河北博野某地热系统现场阻垢试验及阻垢效果评价

2022-10-24谭现锋卜宪标郑慧铭

河北工程大学学报(自然科学版) 2022年3期

刘肖，谭现锋*，张丰，张茜，卜宪标，郑慧铭

(1.山东省鲁南地质工程勘察院(山东省地质矿产勘查开发局第二地质大队)，山东济宁 272100；2.山东省地热清洁能源勘查开发工程研究中心，山东济宁 272100；3.中国科学院广州能源研究所，广东广州 510640)

随着国家碳排放目标的提出，地热资源作为一种绿色高效能源有助于“双碳”目标实现，它主要源于地球内部物质中放射性同位素衰变和地球形成演化过程中能量的积累[1-2]，可用于发电、供暖、制冷和干燥等[3-8]。在地热资源开发利用过程中，当地热流体从热储层向井口运移过程中，由于温度和压力发生变化，导致矿物溶解度发生变化，使得部分矿物达到过饱和状态而沉淀，堵塞井孔或管道，严重制约了地热能的可持续高效开发利用[9-10]。沉淀物的主要类型为钙垢、硅垢等，钙垢一般发生在闪蒸点附近，对地热井而言结垢现象在井筒、井上管道设备中均有可能出现，使地热能的利用效率大大降低[10-11]。钙垢在中低温系统中比较常见，也是影响地热系统正常高效运行的主要问题之一。羊八井地热电站[12]、西藏那曲地热电站[8,10]、甘孜DZK02地热井[8,10,13]、台湾3 MW清水地热电站[14]、Los Humeros 地热井[15]、Mahanagdong地热井[16]均因碳酸钙结垢问题影响了地热田产能以及被迫关闭。

为了提高地热资源的开发利用效益，避免或者抑制钙垢，国内外学者在钙垢形成过程及其碳酸钙溶解沉淀规律、钙垢发生的深度和范围、闪蒸趋势及结垢位置预测、地热水阻垢室内试验及机理分析等方面均进行了大量研究。Rahmani[17]采用WATCH软件对伊朗两口地热井进行绝热沸腾计算，结果表明当温度为180 ℃时，CO32-突增，流体中方解石过饱和程度最大，基于结垢速率，建议采取周期性机械除垢或者注入抑制剂的方法来除垢；Wangyal[18]采用井筒模拟程序HOLA评价沸腾是否发生及发生的深度，结果表明该深度与井口压力关系密切；Patzay[19]对匈牙利南部高温地热井结垢的影响因素进行了分析；Lee[20]对台湾清水地热井的闪蒸位置进行了模拟研究，并基于此进行了阻垢试验；Ramos-Candelaria[21]模拟了在不同井口压力、井口质量流量下的闪蒸位置，并设计了阻垢系统；郑尚基[22]以丰顺地热电站的地热水为研究对象，采用地热水动态结垢趋势分析系统和多功能岩心驱替装置，开展了地热电厂管线和地下储层的结垢、阻垢模拟实验，研究了其结垢趋势和阻垢防垢效果。Li[23]以川西高原康定地热田为例，对地热田方解石结垢和腐蚀进行了定量评估。由此可以看出针对地热系统碳酸钙结垢情况，如何对结垢情况进行抑制及避免十分必要。

尽管针对地热系统结垢问题已有理论研究和阻垢试验研究，但是多集中在阻垢的某一个环节进行研究分析，缺少对阻垢过程的全流程研究。因此本文以河北博野某结垢地热系统为例，在现场考察、结垢样品取样测试、阻垢工艺设计及现场加注装置研发加工、阻垢现场试验、阻垢效率评价等全过程基础上，开展不同工况下阻垢剂的阻垢性能试验，对阻垢效果进行定量评价，分析阻垢机理，为实现地热资源可持续高效开发利用，减少钙垢的产生提供了理论依据。

1 试验场地基本情况

1.1 地热井基本概况

河北省博野县地热区位于渤海湾盆地冀中坳陷高阳低凸起边缘，构造位置处于华北板块(1级)、渤海湾盆地(2级)、冀中坳陷(3级)、高阳低凸起(4级构造单元)的中东部边缘，高阳低凸起与饶阳凹陷交界处(图1)。本次现场阻垢试验选定地热井于2019年11月成井，成井深度3 860 m。井管结构如图2所示，0.00～550.00 m套管直径406 mm，513.66～1 950.00 m套管直径273 mm，1 909.86～3 275.00 m套管直径177.80 mm，3 275.00～3 860.00 m裸眼直径152.00 mm。本井自上而下揭露地层为第四系平原组，新近系明华镇组，上第三系馆陶组，古近系东营组，下第三系沙河街组、孔店组，蓟县系雾迷山组。地热井的取水层位为雾迷山组白云岩含水层，利用井段3 275.00～3 860.00 m。含水层总厚度38.10 m，共11层，其中二类裂隙缝层2层，总厚度9.10 m，三类裂隙缝层9层，总厚度29.00 m，二类、三类裂缝含水层为本井主要产水层。

1.2 供暖系统运行情况

2019—2020年采暖季该井作为回灌井使用，2020—2021年采暖季作为开采井使用。该井泵室管直径406 mm，深度513.66 m，电潜泵下深230 m，泵管直径159 mm，供暖系统如图3所示。井口的地热水直接进入敞口水箱，然后流进热水罐。热水罐内的地热水为常压，其内部的地热水经地热水泵加压后进入板换，在板换内地热水和建筑供暖的二次水换热。板换内换完热后的地热水分为两路，一路流入回灌井完成回灌；另一路通过旁侧管进入生产井重复利用。旁侧管的直径140 mm，深度400 m，见图3。地热井出水温度95 ℃左右，流量约100 m3/h，，供暖面积约18×104m2。该地热井所用电潜泵的耐温不超过100 ℃，而井下温度近116 ℃，因此通过设置旁侧管的方式，将部分38.5 ℃的回水掺入开采流体中，降低泵前流体温度，以防止电潜泵被高温烧坏。井口地热水的化学成分见表1。溶解性总固体5 122 mg/L，水化学类型为HCO3·Cl-Na型。

1.3 结垢情况分析

2020—2021年采暖季供暖期间，地热供暖系统结垢现象严重，在井口管道、敞口水箱、热水罐、阀门、地热水泵、地热水输运管线和板换位置均有垢生成，结垢最为严重的是热水罐，地热水泵以及两者之间的管道和阀门，导致地热水泵半个月就需要清垢一次。图4是管道和阀门结垢情况，图4(a)中结垢位置位于图3中的A点，为地热井出井口后的水平运输管道上，图4(b)是热水罐和地热水泵之间管道上的阀门。在井口处安装了压力和温度监测点以及取样口。

图1 冀中坳陷及饶阳凹陷构造位置图[24]Fig.1 Structural location map of Jizhong depression and Raoyang depression[24]

图2 井管结构示意图Fig.2 Schematic diagram of well casing structure

注：1—电潜泵，2—井管，3—泵管，4—敞口水箱， 5—热水罐，6—地热水泵，7—板式换热器，8—供暖循环泵，9— 建筑，10—回灌管线，11—回灌井，12—旁侧管，A—取样口1， B—取样口2，C—阻垢剂加注装置。图3 地热供暖系统图Fig.3 Geothermal heating system diagram

对图中的垢样进行晶型结构分析，分析仪器为D8 Advance衍射仪。XRD谱图显示两处结垢样品为同一晶型结构，未检测到其他晶型结构(图5)，垢质主要成分为方解石和文石，质量比为45∶55，为碳酸钙结垢，敞口水箱和热水罐的生锈导致了图4结垢颜色的不同。

对图4(a)管道结垢样品中选定的区域进行二次电子样品形貌及能谱分析，分析仪器为扫描电子显微镜/Merlin Compact，加速电压EHT=10.00 kV，工作距离WD=9.4 mm，Mag=500X(即放大倍数500倍)，检测器为 SE2探测器(Signal A=SE2)。图6(c)—(f)为Ca、O、C、Mg在区域上分布图，图中亮点为元素含量高的地方，阴影部位为元素含量较低区域。结合能谱分析结果(表2)，管道内结垢样品主要元素为Ca、O、C、Mg，Ca含量41.95%、O含量45.18%、C含量8.75%、Mg含量1.55%。

图4 结垢情况[25]Fig.4 Scaling situation[25]

对图4(b)阀门结垢样品中选定的区域进行二次电子样品形貌分析，分析仪器为扫描电子显微镜/Merlin Compact，加速电压EHT=15.00 kV，工作距离WD=7.9 mm，Mag=100X(即放大倍数100倍)，检测器为 SE2探测器(Signal A=SE2)。图7(c)—(f)为Ca、O、Mg、Fe在区域上分布图，图中亮点为元素含量高的地方，阴影部分为元素含量较低区域。结合能谱分析结果(表2)，管道内结垢样品主要元素为Ca、O、C、Mg，Ca含量37.90%、O含量48.09%、C含量10.82%、Mg含量1.75%。

图5 结垢样品XRD图谱Fig.5 XRD pattern of scaling samples

2 现场阻垢试验及分析

2.1 阻垢剂加注设备研发

根据现场实地考察，结合试验场地供暖管网及地热井情况，为解决阻垢剂流量计量、加注深度计量及调节控制绞盘、测量仪表设计安装等问题，自主研发了一套阻垢剂地面加注装置(图8)。其中隔膜式计量泵两台，一备一用，功率90 W，流量150 L；加注管线直径6 mm，316 L材质；阻垢剂加注箱为直径20 cm、高1.0 m的亚克力材质管；阀门和管线均耐腐蚀性。图3中的C位置即为阻垢剂加注装置。

2.2 阻垢剂选择

针对地热系统碳酸钙结垢及相应阻垢剂注入情况，不同地热系统采用不同阻垢剂进行处理，减少碳酸钙结垢形成主要是减少碳酸钙形成的量或抑制碳酸钙晶体的生长，添加的化学抑制剂主要包括有机磷酸盐[26]、腐植酸[27]、氨基三亚甲基膦酸(ATMP)、多胺聚醚亚甲基膦酸酯(PAPEMP)[28]、乙二胺四乙酸(EDTA)[29]、金属离子[30]、聚丙烯酸(PAA)[28,31]、聚马来酸(PMA)[32]和聚甲基丙烯酸(PMAA)、聚丙烯酸-2-丙烯酰胺-2-甲基丙烷磺酸)(PAA-PAMPS)[33]、。这些抑制剂吸附在生长的碳酸钙晶体上并与生长部位结合，抑制晶体生长并破坏方解石晶体的规则形状，导致晶体稳定性减弱。目前纳米气泡技术应用于酸性地热流体中也取得了显著效果[34-35]。根据对国内外阻垢剂市场的调研，本次选择了四种阻垢剂进行现场试验，分别为阻垢剂1126、1131、4200、ZC WT-1100，前三种为进口阻垢剂，主要成分分别为丙烯酸共聚物，磺化离子聚合物和聚丙烯酸钠，第四种为国产阻垢剂羧酸均聚物。阻垢剂1126在含有二价阳离子的水中，比如钙和铁的水中，具有很好的防垢性能，在广泛的pH条件下具有优异的热稳定性和配方能力。阻垢剂1131是一种水性压克力助剂，能有效地抑制方解石(CaCO3)，重晶石(BaSO4)以及硫酸钙(CaSO4)等矿物质结垢问题，表现出极好的二价离子兼容性、极好的高温稳定性、极好的与防冻剂兼容稳定性等。阻垢剂4200对硫酸钡、硫酸钙、碳酸钙和其他低溶解度盐类有高效的阻垢效果。阻垢剂ZC WT-1100对碳酸钙、硫酸钡、硫酸钙结垢具有较好阻垢效果。

表1 井口地热水的离子浓度Tab.1 Ion concentration of geothermal water

表2 区域元素含量百分比一览表Tab.2 List of regional element content percentages

图6 管道结垢样品二次电子样品形貌图Fig.6 Morphology of secondary electronic sample of pipeline scaling sample

图7 阀门结垢样品二次电子样品形貌图Fig.7 Morphology of secondary electronic sample of valve scaling sample

图8 阻垢剂加注装置示意图Fig.8 Schematic diagram of scale inhibitor filling device

表3 不同阻垢剂在不同工况下阻垢效果一览表Tab.3 List of scale inhibition effect of different scale inhibitors under different working conditions

2.3 现场试验方案

本次试验通过泵房内旁侧管阀门(图3中C点)加入阻垢剂，在抽水井井口(图3中A点)和回水回灌后管道(图3中B点)进行取样。四种阻垢剂分别开展了不同工况下的试验来研究阻垢剂最佳浓度、最佳阻垢效果。结合已有研究成果，土耳其Germencik地热田采用膦酸酯和聚丙烯酸酯阻垢剂，阻垢剂用量6 mg/L可以达到阻垢效果；菲律宾BacMan地热田采用膦基羧酸共聚物，用量8 mg/L[11,16,20]。通过对比研究博野场地水质情况，结垢情况等，选定本次试验最低浓度10 mg/L。试验开始前，按照1∶3稀释阻垢剂置于阻垢剂箱，通过调节计量泵的阀门达到调节阻垢剂浓度的目的。试验开始前分别在取样口1、2取样，作为试验开始时的初始值。试验开始后，每30 min同时在取样口1、2进行取样，现场采用EDTA滴定法测试样品中的钙离子浓度。试验结束后，选取关键节点数据送往实验室进行检测。

本次试验通过分析加入阻垢剂前后Ca2+浓度变化差异来计算阻垢率，见式(1)：

(1)

其中，IE 是阻垢剂阻垢率(%)；Ct0和Ct1为注入阻垢剂前后取样点B的Ca2+浓度(mg/L)；Cwb为井底Ca2+浓度(mg/L)。

在井底、井口和旁侧管中的总水量和钙离子浓度是保持守恒的，因为井底的钙离子浓度可以通过公式(2)、公式(3)计算获得：

MW,wh=MW,st+MW,wb

(2)

Cwh1MW,wh=Ct1MW,t+CwbMW,wb

(3)

其中，MW,wh、MW,t和MW,wb分别为井口、旁侧管和井底流量，(kg/s)；Cwh1为取样点A Ca2+浓度(mg/L)。MW,wh、MW,t和MW,wb是通过超声波流量计测量得到。

2.4 试验结果分析

在2020—2021年供暖季开展了现场阻垢试验，开采井抽水量在110 m3/h左右、旁侧管流量26～30 m3/h，出水温度95 ℃左右，图9为不同阻垢剂的阻垢率变化曲线图，从图中可以看出阻垢率随着阻垢剂浓度增大而增大，阻垢剂浓度越高，阻垢率越大，阻垢效果越好。通过现场试验发现，4种阻垢剂的阻垢率均能达到80%；达到85%的有2种阻垢剂，为阻垢剂1131在阻垢剂加入量25.61 mg/L和阻垢剂4200在阻垢剂加入量48.82 mg/L时；达到92%的有1种阻垢剂(表3)。在同等阻垢剂浓度情况下，阻垢剂1131的阻垢效果最好，阻垢剂1131在浓度35.95 mg/L时，阻垢率为94.49%，四种阻垢剂中ZC WT-1100阻垢剂相比于其他三种阻垢剂阻垢效果较差。

图9 不同阻垢剂阻垢率曲线图Fig.9 Graph of scale inhibition rate of different scale inhibitors

阻垢剂通常是分散剂，可防止水垢紧紧粘附在管道和设备表面[35]。在碳酸钙垢生成过程中，添加到地热水中的阻垢剂主要通过与水中离子的络合增溶、晶格畸变、凝聚与分散等机理(图10)来达到抑制垢和消除垢的目的[22]。针对现场阻垢试验，筛选出的高效阻垢剂通过吸附水中的钙离子于结垢晶体表面后抑制地热水结垢的正常生长，从而使其发生畸变。并且阻垢剂还可与碳酸钙垢发生反应，破坏碳酸钙的晶格结构，产生不规则的非晶格颗粒，使得垢质变松软，从而达到除垢的效果。

图10 阻垢机理示意图[34]Fig.10 Schematic diagram of scale inhibition mechanism[34]

实验开始前均进行了本底值测试，待整个工况稳定后开始加入阻垢剂。对4种阻垢剂在不同取样口处Ca2+浓度进行分析(图11)，实验过程中各种阻垢剂浓度由小到大进行调整，随着阻垢剂的加入取样口处Ca2+浓度整体均表现出上升趋势。对于每个浓度的阻垢剂，A、B两个取样口的Ca2+浓度在3 h后可以达到稳定，因此每种阻垢剂在每个工况下需要4～5 h的测试时间才能达到稳定的Ca2+浓度。阻垢剂1126在阻垢剂浓度11.30 mg/L情况下，井口Ca2+浓度达到稳定需3 h，稳定在93.27 mg/L，此时取样口B Ca2+浓度稳定在81.10 mg/L，加大阻垢剂浓度到30.43 mg/L时，取样口A Ca2+浓度为101.38 mg/L，此时取样口B Ca2+浓度达到93.27 mg/L，此后稳定在该浓度直到试验结束。阻垢剂4200在阻垢剂浓度48.82 mg/L情况下，取样口A Ca2+浓度为99.35 mg/L，取样口B Ca2+浓度稳定在93.27 mg/L，阻垢剂ZC WT-1100在阻垢剂浓度38.56 mg/L情况下，取样口A Ca2+浓度为101.38 mg/L，此时取样口B Ca2+浓度稳定在93.27 mg/L，三种阻垢剂最后稳定浓度一致，这与图9中的阻垢率曲线图相一致，三种阻垢剂最大阻垢率均在80%左右。而阻垢剂1131，在相同阻垢剂浓度情况下，取样口Ca2+浓度稳定值较高，在阻垢剂浓度35.95 mg/L情况下，取样口B Ca2+浓度达到95.30 mg/L，此时阻垢率为94.49%，阻垢效果最好。