Feasibility study on joint exploitation of methane hydrate with deep geothermal energy
-
摘要: 随着全球能源消耗不断增加，天然气水合物和地热资源具有储量丰富、清洁高效等优势成为世界研究的热点，中国南海海域同时具有丰富的水合物资源和地热资源。由此，提出了联合深层地热资源开采浅部水合物的方法，通过向深层地热储层注入海水，海水在深层地热中吸收热量后循环至浅部水合物储层，结合降压法和注热法促使水合物分解。利用数值模拟对联合法的可行性进行评估，并对地层热物性、开采参数等储层敏感性进行分析。模拟结果表明：联合法能够有效地利用深层地热将海水加热，海水进入水合物层时的温度保持约为50℃，与注热法和降压法相比具有更高的产气量，具有良好的可行性；注入速度、井底压力、地层导热系数和地温梯度对联合法开采效果具有显著影响；注入速度和井底压力对前期的产气效果影响较大，而较大的地层导热系数有利于海水与地层的换热；地温梯度小于0.025 m/℃时，联合方法的换热性能极大减弱，甲烷累计产气量大幅度降低，联合方法的商业价值降低，可行性减弱。Abstract: With the drastic increase in energy consumption, both the natural gas hydrates and geothermal resources have become research focuses in the world due to their enormous reserves. Then the method to exploit shallow gas hydrates together with deep geothermal resources becomes attractive. In this method, seawater will be injected into the deep geothermal reservoirs and then bring into the shallow hydrate reservoirs after circulation and absorbing enough heat from the deep geothermal reservoir. Depressurization and thermal recovery technique are used to encourage the decomposition of gas hydrates. In this paper, the feasibility of the joint exploitation method is evaluated through numerical simulation, and the thermal properties of hydrate bearing sediment, exploiting parameters and reservoir sensitivities studied. Results show that effective utilization of geothermal resources to heat seawater may enable the temperature of seawater entering the hydrate layer to maintain on a level of about 50℃, and higher gas production will achieved comparing to the method of heat injection and depressurization techniques. It is indeed a method with good feasibility. Some factors, such as injecting rate, bottom hole pressure, thermal conductive factor of formation and geothermal gradient, have a significant impact on the exploiting results. In addition, the injection rate and bottom hole pressure may bring greatly influence to gas production in the early stage, while the thermal conductivity of larger formation has a favorable contribution to the heat exchange between the seawater and formation. Results also suggest that the performance of heat transferring of the method be largely attenuated and the cumulative gas production of methane be substantially reduced within the area with a geothermal gradient lower than 0.025m/℃. In such a circumstance, the commercial value of the deposits and their feasibility of exploitation will decrease.
-
Keywords:
- geothermal /
- natural gas hydrate /
- depressurization /
- thermal recovery /
- feasibility
-
天然气水合物是一种储量丰富的清洁能源，广泛分布于深海沉积物和陆地的永久冻土层，其全球储量相当于目前已探明常规能源的两倍以上[1]。目前开采水合物的方法主要包括热激发法(或注热法、热激法)、降压法、抑制剂注入法和置换法。注热法、降压法和抑制剂注入法通过破坏水合物平衡条件或者改变水合物储层的环境，使水合物处于热力学不稳定状态发生分解。降压法在开采过程中不需向水合物储层注入热量，是一种经济环保的开采方法。但随水合物不断分解吸收热量，引起储层温度降低导致水合物的二次生成，并且固体水合物的分解破坏了水合物储层的地层结构，容易引起沉降坍塌等地质问题[2, 3]。相比于降压法，注热法通过向水合物储层注入热流体加热储层，热流体为水合物分解提供了足够的热量，保证水合物能够持续分解，是一种较为理想的水合物开采方式。但是注热法存在能源消耗大、经济效益不佳和热效率低等问题[4-6]。Wang建立了多孔介质水合物解离过程中传质和传热的解析模型，分析了降压法、注热法和降压法与注热法相结合的3种开采方式传质与传热机制，证实了降压法与注热法相结合的方式能够提高水合物的分解速度[7-9]。CO2置换法是向水合物储层中注入CO2气体，由于二氧化碳水合物比甲烷水合物稳定，二氧化碳分子会置换出甲烷分子。Ohgaki等第一次通过实验证实了CO2置换开采甲烷水合物的可行性[10]。甲烷水合物和二氧化碳水合物的分解热分别为54.49和57.98kJ/mol，根据热力学理论CO2置换甲烷水合物的吉布斯自由能为负，化学反应能够自发进行, 但是CO2置换反应过程却十分缓慢[11-13]。地热能是一种环境友好、资源丰富的可再生清洁能源。根据2010年世界地热大会统计，截至到2009年，大约有70多个国家对地热直接利用，年利用总量达72622GW。据2015年世界地热大会报道，中国地热发电的装机容量达到27MW[14]。我国南海珠江口盆地具有丰富的水合物资源和地热资源，2013年，中国在珠江口盆地东部海域进行钻探，成功获得天然气水合物样品[15]。2017年, 中国在南海神狐海域成功实施天然气水合物试采，试采稳步持续了60天, 天然气总产量达到30.9×104m3[16]。珠江口盆地是典型的“热”盆地，蕴藏着丰富的地热资源，其地热梯度范围为24.7~60.8℃/km, 平均值为37.9±7.4℃/km[17]。在南海珠江口盆地利用深层地热资源开采水合物似乎具有良好的可行性，目前，诸多学者提出了利用地热资源开采水合物的方法。宁伏龙等首先提出利用地热开采海底天然气水合物的设想，通过泵将海水从注入井压入干热岩层，被加热的海水沿注水管和套管的孔隙上返至水合物储层[18]。窦斌等对利用海底地热开采水合物的方法进行了可行性分析，通过块裂介质模型模拟显示出随着地热系统的运行越来越有利于水合物的开采[19]。Liu提出了CO2置换法联合地热开采水合物的方法，将CO2从热交换井注入地热储层，由地热储层加热的CO2向上流入水合物储层，以注热法与CO2置换相结合的方式完成了天然气水合物的开采，同时实现了CO2地质埋存。但是由于换热井为垂直井，导致整个系统的换热性能受到限制，同时CO2气源也会限制该方法在实际过程中的实施[20]。
目前，针对于利用地热开采水合物的研究尚不完善，诸多研究只是提出利用地热开采水合物的设想，开采过程中存在热流体漏失和如何保证被加热的流体回流至采出端等问题没有给出详细的解决方法，而对于利用地热开采水合物的可行性评价以及影响因素分析尚不成熟。因此，本文提出一种联合深层地热利用注热法和降压法开采海域水合物的方法，系统地阐述了联合深层地热海域水合物开采的工艺设计，通过数值模拟方法对其可行性进行评估，并对联合开采方法的敏感性进行分析。
1. 联合深层地热海域甲烷水合物开采工艺设计
中国南海区域拥有丰富地热和天然气水合物资源，诸多学者提出了利用深层地热资源开采水合物的设想，基于现有地热水合物联合开采的基础，设计了一套通过向换热井注入海水经深层地热储层吸收热量上返至水合物储层的联合开采方法。本文利用海水作为循环介质，海水来源广、成本低，并且含有大量的盐，渗入到水合物储层中有利传热，加快水合物的分解。联合开采方法系统示意如图 1所示，海水从换热井注入，流经地热储层被加热，上返至浅层采出井水平段进入水合物储层，生产井的水平段与采出井的水平段位于同一水合物层位。图 1中红色管柱代表换热井，换热井主要包括换热井水平段和换热井垂直段两部分，换热井水平段位于目标地热储层，是海水换热的主要场所。蓝色管柱代表采出井，采出井主要用于将加热的海水收集并运送到水合物储层，主要由水平段和垂直段组成，水平段进行射孔完井，垂直段底端与换热井水平段相连形成类似“U”型结构。经过地热储层加热的海水直接进入采出井，由于整个循环过程都是在油管内进行，避免了海水漏失到地热储层的问题。海水在采出井井筒内上返时会与上覆岩层发生热交换，上返井段内的油管使用绝热材料以减缓海水上返时的热损失，图 1中黑色管柱为使用绝热材料井段。生产井保持较低的井底压力从而实现降压开采的效果，冯景春等通过室内物理实验对双水平井降压与注温水联合开采法进行研究，研究结果表明水平井相比于垂直井开采水合物，双水平井能效比更高，因此，生产井主要设置成水平井[21]。
![]()
图 1 联合开采方法系统a:换热井和采出井结构, b:生产井井底与采出井水平段Figure 1. Schematic illustration of joint exploiting method systema:Heat exchange well and production well, b:The bottom hole of production well and the horizontal section of the production well相比于常规水合物开采方法，本文提出的联合开采方法结合了注热法和降压法2种方法的优势。李淑霞等利用实验和仿真模拟探究了注热法结合降压法开采水合物，结果表明注热法与降压法结合具有更高的能量利用效率[22]。联合开采方法克服了降压法因长期开采导致温度降低生成二次水合物的问题，又避免了注热法能源消耗严重等问题；由于利用深层地热加热海水，不需要地面的加热设备，节省了地面设备的投入，节约了大量的化石能源，同时海水来源广，产出水可以回注到地下重复利用；整个循环在油管内进行，避免了海水在储层漏失的问题。综上所述，本文提出的联合开采方式是一种经济环保、高效节能开采措施，具有良好的可行性。
2. 水合物开采数值模拟分析
2.1 数值模型描述
本次研究利用CMG-STARS模拟器完成了联合深层地热开采海域水合物的数值模拟分析，数值模型考虑了3个组分，分别为液相水、气相甲烷和固体甲烷水合物。由于利用了注热法开采水合物，模型忽略了开采水合物过程中冰的形成。化学反应(1)表示了水合物的形成和分解过程，本文采用Kim第一次定量研究的水合物分解动力学方程(2)[23]，水合物的形成动力学方程采用了Vysniauskas-Bishinoi模型，表达式为方程式(3)[24]，甲烷水合物分解热为47.5 ~ 62.9kJ/mol[25]。
基于中国南海海域的地质特征建立数值模型，为了便于模拟研究将换热井和采出井用一口井等效井(INJ-1)替代，即用一口井(INJ-1)完成海水经地热储层加热上返至水合物储层的整个过程。图 2为数值模拟中井位设置，整个联合开采方法模拟系统由等效井(INJ-1)、和生产井(PRO，PRO-1)组成，生产井井底与等效井井底距离为54m。为了提高计算精度和计算速度，整个模型被剖分为41×41×40即67240个网格，平面单个网格尺寸为6m×6m，水合物储层纵向网格尺寸为3m，而下覆岩层纵向网格尺寸为120m，地热储层纵向网格尺寸为50m。根据GMGS3-W19站位资料显示海底水深为1273.9m，海底温度为3.75℃，地温梯度为0.045℃/m，水合物储层主要赋存于海底以下135~170m[26]，不含有自由气，初始含气饱和度为0，初始含水饱和度为0.6，非水合物储层初始含水饱和度为1。模型没有考虑水合物储层上部的覆盖层，模型顶部深度为1408m，地热储层底部深度为4388 m，具体模型参数设置如表 1和表 2。模型渗透率为40mD，孔隙度为0.3[27]，模型外边界为封闭边界。
表 1 水合物动力学参数设置Table 1. Parameters of hydrate kinetic properties参数 CH4·nH2O Kd0, gmole/(s·Pa·m2) 1.24×105 Kf0, gmole/(s·Pa·m2) 2.9×103 Adec, m2/m3 3.75×105 ΔE, J/mol 81084.2 H, J/mol 54490.0 表 2 数值模型参数设置Table 2. Parameter for numerical models参数 数值 单位 顶部深度 1408 m 水合物储层厚度 30 m 岩石体积热容量 2.12×106 J/(m3·K) 水合物热熔 1600 J/(kg·K) 岩石导热系数 5 W/(m·K) 水合物导热系数 3.92 W/(m·K) 水的导热系数 0.60 W/(m·K) 气体导热系数 0.03 W/(m·K) 绝缘层导热系数 0.2 W/(m·K) 水合物储层初始含水饱和度 0.6 － 水合物储层初始含气饱和度 0 － 水合物储层初始水合物饱和度 0.4 － 深度，d - m 初始地层压力 9800*d Pa 岩石密度 2650 kg/m3 CH4·nH2O密度 919.7 kg/m3 $$ {\text{C}}{{\text{H}}_4}\left( g \right) + n{{\text{H}}_2}{\text{O}}\;\;f\;\;{\text{C}}{{\text{H}}_4} \cdot n{{\text{H}}_2}{\text{O}}\left( s \right) $$ (1) $$ {v_d} = K_d^0{A_{dec}}\left( {{P_e}-{P_g}} \right)\exp \left( {-\frac{{\Delta E}}{{RT}}} \right) $$ (2) $$ {v_f} = K_f^0{A_{dec}}\left( {{P_e}-{P_g}} \right)\exp \left( {-\frac{{\Delta E}}{{RT}}} \right) $$ (3) 水合物的形成和分解过程可以用化学反应(1)(2)和(3)来表示，式中vd，水合物分解速率；Kd0, 水合物分解速率常数，gmole/(s·Pa·m2); Adec，单位体积内反应面积，m2/m3; Pe, 气相平衡压力，Pa; Pg, 气相分压力，Pa; △E, 水合物反应活化能，J/mol; R, 气体常数，J/(molK); T, 温度, K; vf, 水合物形成速率；Kf0, 水合物分解速率常数，gmole/(s·Pa·m2)。
2.2 可行性分析
井网结构布置如图 2所示，两口生产井水平段长度为246m，生产井井底压力保持在2.5MPa，注入井保持注入速度为200m3/d，注入温度为5℃。为了充分证明联合开采方法有良好的可行性，与注热法和降压法进行对比。由于海水的初始温度对注热法开采水合物具有较大的影响，初始温度越高注热开采效果越佳，但过高的初始温度需要消耗大量的能源。综合考虑，本次模拟注热法海水注入温度为30℃。
图 3展示了换热井水平段所在网格地层温度的分布，对比原始地层温度分布和开采1年、5年和10年时的地层温度分布，海水刚进入水平段时温差较大，根据傅里叶定律可知，在相同条件下温差越大热流密度越大，即单位时间内通过单位面积的热量越多，这导致了水平段入口端地层温度要比出口端地层温度下降快。水平段入口端地层温度在10年内下降了46℃，而随着海水在水平段内不断地吸热升温，与地热储层岩体的温差逐渐变小，10年内地层温度下降约为21℃。图 3可以看出在整个生产过程中地热储层的温度能够保持在100℃以上，海水进入采出井时的温度达到100℃以上，深层地热能够为海水加热提供稳定的热源。图 4曲线展示了采出井在水合物储层入口端所在网格地层温度变化，可以看出尽管海水在采出井内上返时发生热损失，进入水合物储层时的温度仍保持在50℃左右，储层最高温度达到55℃，能够满足注热法开采水合物的要求。图 5表示了开发10年时换热井水平段温度的空间分布，从图中可以看出随着海水不断与换热井周围岩体进行热交换，换热井两侧20m内岩体的温度逐渐降低，形成了一个40m宽的温度降低区域，井筒内的温度稳定在111℃，20m外的岩体温度仍能保持在原始状态，这说明地热储层能够稳定地为水合物开采提供热量，不会因换热而引起整个地热储层温度的严重降低，地热储层蕴含的热量远远高于海水换热所需要的热量，因此, 整个开采过程中拥有充足的能量补充。图 6表示了开发10年时水合物储层的压力空间分布，随着气体水合物和水的产出，整个水合物储层的平均压力由14 MPa逐渐降低至约3.2MPa，由于气体的产生和海水的补充，水合物储层压力没有降低到生产井井底压力，生产井井底周围的压力保持在2.7MPa，采出井井底周围的压力保持在3.8MPa，存在的压力差有利于海水流向生产井。图 7、图 8和图 9分别表示开采10年时水合物储层的温度分布，温度沿着采出井向生产井井底逐渐降低，生产井井底温度达到33℃，注热法开发10年时生产井井底的温度只有10℃，降压法井底周围的温度下降到1.5℃，海水流向生产井井底的过程中对水合物储层进行加热，为水合物开采提供热量。图 10为甲烷累计产气量和产气速度曲线，生产初期甲烷累计产气量曲线基本重合，水合物的分解主要通过降压法，而随着水合物的分解吸收了大量热量，联合开采方式能够持续不断地为水合物分解提供热量，使水合物保持稳定的产气速度，拥有较高的甲烷累计产量，而对比降压法，随着开采的进行地层温度逐渐降低，产气速度下降，导致最终的累计产气量较低。
![]()
图 5 换热井水平段温度分布Figure 5. Horizontal temperature distribution of the heat exchange well![]()
图 7 联合法水合物储层温度分布Figure 7. Temperature distribution of a hydrate reservoir by the joint method![]()
图 8 注热法水合物储层温度分布Figure 8. Temperature distribution of hydrate reservoir by the thermal method![]()
图 9 降压法水合物储层温度分布Figure 9. Temperature distribution of hydrate reservoir by depressurization method综合而言，整个联合开采过程中，地热储层为水合物开采提供一个稳定的热源，保证了水合物分解所需要的热量，具有较高的甲烷累计产量，同时联合法不需要地面加热设备，极大地降低了能源消耗，节省了开采成本，具有良好的可行性。
3. 敏感性分析
3.1 注入速度
保持其余参数不变，注入速度分别设置为200、300、400和500m3/d。图 11展示了生产10年时不同注入速度条件下沿着水平段的地层温度分布。注入速度为200m3/d时，水平段入口端温度下降至85℃，这说明地热储层与海水进行了充分换热，由曲线可以看出海水经入口端流至采出井井底过程中不断被加热，水平段50m后的地层温度基本维持不变，这说明海水的温度经加热基本与地层温度持平，可以看出海水进入采出井井底的温度大约113℃。注入速度为500m3/d时，水平段入口端温度为76℃，海水进入采出井井底温度大约110℃。这是因为海水在入口端温度较低，海水与地热储层温差较大，地热储层与海水换热效果明显，导致了入口端地层温度较低。注入速度越大，入口端的地层温度下降越大，但是不同速度条件下海水到达采出井井底时的温度差异较小，这是因为换热井水平段较长，海水能够在地热储层进行充分换热升温。图 12所示为不同注入速度条件下采出井水平段所在水合物储层入口端的地层温度变化，图中可以明显地看出注入速度越小，入口端的地层温度越高，注入速度为200m3/d时，入口处地层温度保持在48℃，而注入速度为500m3/d时，入口处地层温度为39℃。图 13和图 14分别为甲烷的累计产气量曲线和产气速度曲线，由图 13可以看出，不同注入速度条件下的累计产气量差距较小，注入速度为200和500 m3/d所对应的甲烷累计产气量分别为3.32×107和3.25×107m3/d。生产初期产气速度急剧增大，达到峰值后迅速下降直到生产速度保持平稳。注入速度为200m3/d时的最大产气速度高于500m3/d，而达到平衡时的产气速度低于500m3/d，这是因为在生产初期主要靠降压法生产，压力差越大产气速度越大，500m3/d补充地层能量充足导致初期生产效果不佳，而生产后期水合物的分解主要通过注热法，因此产气量效果优于200m3/d。由于本次的模拟区域较小，高的注入速度具有较大的波及区域的优势没有得到体现，而较低的注入速度就能提供水合物分解所需的热量，并且低的注入速度加热地层的效果优于高的注入速度，这导致了本次模拟中低的注入速度具有更大的累计产气量。
3.2 井底压力
本文探究了生产井井底压力为2.5、3.5、4.5和5.5 MPa对水合物开采的影响，图 15和图 16为甲烷累计产气量和产气速度曲线。在生产初期，生产井井底周围的水合物受到生产井降压的影响，水合物的相平衡状态被破坏，分解成水和甲烷气体，在压差的作用下迅速的流入生产井，由图 15和图 16可以看出井底压力越小，生产井的前期的产气速度越快，井底压力为2.5MPa时的最大产气速度达到1.23×105m3/d，井底压力为5.5MPa时的最大产气速度仅仅为1.8×104m3/d，最大产气速度随井底压力的减小而增大。井底压力主要影响前期的产气速度，产气速度到达峰值之后急剧下降而趋于平稳，后期产气速度趋于平缓时，不同井底压力的产气速度十分接近。这是因为前期海水的驱替前缘还未波及至生产井，生产井井底周围水合物的分解主要通过降压法，井底压力越低水合物分解速度越快，进而产气速度越大。但是水合物分解吸收了大量热量，导致采出井井底地层温度逐渐降低，产气速度迅速下降。随着海水驱替前缘向生产井井底扩散，及时补充了降压法所消耗的地层热量，与此同时地层能量也得到了补充，降压法的作用减弱，甲烷气体主要通过海水加热促使水合物分解产生，因此, 后期水合物产气速度保持平稳且相接近。
3.3 地层导热系数
维持其余参数不变，地层的导热系数设置为3、5、7和9 W/(m·K)。图 17是开发10年时地层导热系数对地热储层换热的影响。由图 17可以看出，地层的导热系数越大，开发10年时的地热储层温度下降越快，这意味着海水进入采出井井底的温度越低，这是因为较大的导热系数有利于地热储层与海水之间的传热，在相同的开发时间内地热储层温度下降较快。图 18采出井水平段所在水合物储层入口端的地层温度变化，开发10年时3、5、7和9W/(m·K)时入口端的地层温度分别为41、47、50和53℃，较大的导热系数具有较好的换热性能，海水流经地热储层时提取热能能力强，同时上返至水合物储层时也能更快的加热水合物储层。图 19显示出3、5、7和9W/(m·K)时的累计产气量分别为3.24×107、3.32×107、3.38×107和3.42×107m3。从图 19和图 20可以看出不同的导热系数甲烷的累计产气量和产气速度较为接近, 这说明地层的导热系数为3 W/(m·K)时，海水经加热能够为水合物的分解提供足够的热量，维持稳定的开采速度。
3.4 地温梯度
设置地温梯度为0.025、0.035、0.045和0.050 m/℃，其余参数保持不变，探究地温梯度对联合开采性能的影响。图 21是开发10年时地温梯度对地热储层换热的影响。地温梯度为0.025、0.035、0.045、0.050m/℃时，海水进入采出井井底的温度分别为65、90、114和127℃，地温梯度越大，海水升温越高。这是因为随着地温梯度的增大，井筒中的海水与地热储层间的温差增大，从而使海水从地热储层提取了更多的热量。采出井水平段所在水合物储层入口端的地层温度变化图 22所示，地层温度为0.05m/℃时，开发10年时，水平段入口端地层温度被加热至52℃，同样情况地温梯度为0.025m/℃时，水平段入口端地层温度仅仅只有28℃。这是因为海水到达采出井井底的温度过低，同时海水上返时与上覆地层发生热损失，因此, 海水进入水合物层的温度就会较低，对水合物储层的热量补充效果变差。这种影响可以从图 23和图 24看出，地温梯度为0.05m/℃时，甲烷气体的累计产量为3.46×107m3，而地温梯度为0.025m/℃时，甲烷累计产气量仅仅为2.54×107m3，地温梯度对甲烷产气量具有较大的影响，呈现地温梯度越大，甲烷的累计产气量越大和产气速度越大的规律。在本次模拟中，地温梯度对联合开采方法的效果具有较大的影响，地温梯度小于0.025m/℃时，联合方法的换热性能极大减弱，甲烷累计产气量大幅度降低，进而导致联合方法的商业价值降低，可行性减弱。因此，实施联合开采方法需要充分考虑地温梯度的影响。
4. 结论
本文提出了一种联合深层地热海域甲烷水合物开采方法，以海水作为携热介质将地热能输送到水合物储层，促进水合物的分解。相比于注热法和降压法，联合法即不需要地面加热，极大地降低了能源消耗，节省了开采成本，又能弥补降压法开采导致生成二次水合物的问题。同时，通过建立三维数值模型对联合法的可行性进行评价以及进行敏感性分析，模拟结果如下：
(1) 联合开采方式能够有效地利用地热资源将海水加热，海水进入采出井时的温度达到100℃以上，进入水合物层时的温度仍保持50℃左右，使水合物保持稳定的产气速度，与注热法和降压法相比具有更高的产气量，联合开采方式具有良好的可行性；
(2) 注入速度主要影响水合物储层压力和温度的变化，生产初期水合物的开采主要通过降压法，生产后期主要通过注热法。注入速度越大，初期的产气速度越小，平稳时的产气速度越大。注入速度为200和500m3/d所对应的甲烷累计产气量分别为3.32×107和3.25×107 m3/d，不同注入速度条件下的累计产气量差距较小；
(3) 井底压力主要影响前期的产气速度，产气速度到达峰值之后急剧下降而趋于平稳。井底压力越小，生产井的前期的产气速度越快，井底压力为2.5 MPa时的最大产气速度达到1.23×105m3/d，井底压力为5.5MPa时的最大产气速度仅仅为1.8×104 m3/d，最大产气速度随井底压力的减小而增大。后期产气速度趋于平缓时，不同井底压力的产气速度十分接近；
(4) 地层的导热系数影响海水在地热储层和水合物储层的换热性能，导热系数越大，海水与地层的换热性能越好，地热储层温度下降越快，海水对水合物层的加热效果越佳。联合开采方法对地层的导热性没有较高的要求，导热系数为3W/(m·K)时，就能有较好的换热效果；
(5) 地温梯度影响对联合法开采具有较大的影响，地温梯度小于0.025m/℃时，联合方法的换热性能极大减弱，甲烷累计产气量大幅度降低，联合方法的商业价值降低，可行性减弱。地温梯度越大越有利于水合物的开采，呈现地温梯度越大，甲烷的累计产气量越大和产气速度越大的规律。地温梯度为0.05m/℃时，甲烷气体的累计产量为3.46×107 m3，而地温梯度为0.025m/℃时，甲烷累计产气量仅仅为2.54×107 m3，相比于前者减少26.6%。
-
![]()
图 1 联合开采方法系统
a:换热井和采出井结构, b:生产井井底与采出井水平段
Figure 1. Schematic illustration of joint exploiting method system
a:Heat exchange well and production well, b:The bottom hole of production well and the horizontal section of the production well
![]()
图 5 换热井水平段温度分布
Figure 5. Horizontal temperature distribution of the heat exchange well
![]()
图 7 联合法水合物储层温度分布
Figure 7. Temperature distribution of a hydrate reservoir by the joint method
![]()
图 8 注热法水合物储层温度分布
Figure 8. Temperature distribution of hydrate reservoir by the thermal method
![]()
图 9 降压法水合物储层温度分布
Figure 9. Temperature distribution of hydrate reservoir by depressurization method
表 1 水合物动力学参数设置
Table 1 Parameters of hydrate kinetic properties
参数 CH4·nH2O Kd0, gmole/(s·Pa·m2) 1.24×105 Kf0, gmole/(s·Pa·m2) 2.9×103 Adec, m2/m3 3.75×105 ΔE, J/mol 81084.2 H, J/mol 54490.0 
下载: 导出CSV
表 2 数值模型参数设置
Table 2 Parameter for numerical models
参数 数值 单位 顶部深度 1408 m 水合物储层厚度 30 m 岩石体积热容量 2.12×106 J/(m3·K) 水合物热熔 1600 J/(kg·K) 岩石导热系数 5 W/(m·K) 水合物导热系数 3.92 W/(m·K) 水的导热系数 0.60 W/(m·K) 气体导热系数 0.03 W/(m·K) 绝缘层导热系数 0.2 W/(m·K) 水合物储层初始含水饱和度 0.6 － 水合物储层初始含气饱和度 0 － 水合物储层初始水合物饱和度 0.4 － 深度，d - m 初始地层压力 9800*d Pa 岩石密度 2650 kg/m3 CH4·nH2O密度 919.7 kg/m3
下载: 导出CSV
-
[1] Klauda J B, Sandler S I. Global distribution of methane hydrate in ocean sediment[J]. Energy & Fuels, 2005, 19(2): 459-470.
[2] Konno Y, Masuda Y, Akamine K, et al. Sustainable gas production from methane hydrate reservoirs by the cyclic depressurization method[J]. Energy Conversion & Management, 2016, 108:439-445. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=c1294ffa686fe21ff5d8962c0a33e6d5
[3] Falser S, Uchida S, Palmer A C, et al. Increased gas production from hydrates by combining depressurization with heating of the wellbore[J]. Energy Fuels, 2012, 26(10):6259-6267. doi: 10.1021/ef3010652
[4] Islam M R. A new recovery technique for gas production from Alaskan gas hydrates[J]. Journal of Petroleum Science & Engineering, 1991, 11(4):267-281.
[5] 唐良广, 李刚, 冯自平, 等.热力法开采天然气水合物的数学模拟[J].天然气工业, 2006, 26(10):105-107. doi: 10.3321/j.issn:1000-0976.2006.10.033 TANG Liangguang, LI Gang, FENG Ziping, et al. Mathematic modeling on thermal recovery of natural gas hydrate[J]. Natureal Gas lndustry, 2006, 26(10):105-107. doi: 10.3321/j.issn:1000-0976.2006.10.033
[6] Fitzgerald G C, Castaldi M J, Schicks J M. Methane hydrate formation and thermal based dissociation behavior in silica glass bead porous media[J]. Industrial & Engineering Chemistry Research, 2014, 53(16):6840-6854. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=57fea6a51ecb6751922735da7f6627b4
[7] Wang Y, Feng J C, Li X S, et al. Experimental and modeling analyses of scaling criteria for methane hydrate dissociation in sediment by depressurization[J]. Applied Energy, 2016, 181:299-309. doi: 10.1016/j.apenergy.2016.08.023
[8] Wang Y, Feng J C, Li X S, et al. Fluid flow mechanisms and heat transfer characteristics of gas recovery from gas-saturated and water-saturated hydrate reservoirs[J]. International Journal of Heat & Mass Transfer, 2018, 118:1115-1127. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=d2a68854081ac571a8bb81c7ffc68998
[9] Wang Y, Feng J C, Li X S, et al. Analytic modeling and large-scale experimental study of mass and heat transfer during hydrate dissociation in sediment with different dissociation methods[J]. Energy, 2015, 90:1931-1948. doi: 10.1016/j.energy.2015.07.029
[10] Ohgaki K, Takano K, Sangawa H, et al. Methane exploitation by carbon dioxide from gas hydrates -phase equilibria for CO2-CH4 mixed hydrate system[J]. Journal of Chemical Engineering of Japan, 1996, 29(3):478-483. doi: 10.1252/jcej.29.478
[11] Yezdimer E M, Cummings P T, Chialvo A A. Determination of the Gibbs free energy of gas replacement in SI clathrate hydrates by molecular simulation[J]. Journal of Physical Chemistry, 2002, 106(34):7982-7987. doi: 10.1021/jp020795r
[12] Uddin M, Coombe D A, Law H S, et al. Numerical studies of gas hydrate formation and decomposition in a geological reservoir[J]. Journal of Energy Resources Technology, 2008, 130(3):032501-1-14. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=3bc22ce759981cba290b6b6f63eeae88
[13] Uddin M, Coombe D, Wright F. Modeling of CO2-hydrate formation in geological reservoirs by injection of CO2 gas[J]. Journal of Energy Resources Technology Transactions of the Asme, 2008, 130(3):102-108. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=56c4d61bb0450e5712a9170354987ab5
[14] Bertani R. Geothermal power generation in the world 2010-2014 update report[J]. Geothermics, 2016, 60(1):31-43. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=008286d0fc7e28ebead100b697d2be77
[15] Liu C, Meng Q, He X, et al. Characterization of natural gas hydrate recovered from Pearl River Mouth basin in South China Sea[J]. Marine & Petroleum Geology, 2015, 61(61):14-21.
[16] 李彦龙.我国海域天然气水合物试开采圆满完成并取得历史性突破[J].海洋地质与第四纪地质, 2017, 37(5):34. http://hydz.chinajournal.net.cn/WKD/WebPublication/paperDigest.aspx?paperID=4f0c459a-7f39-4e5f-9697-dcf201702507 LI Yanlong. Successful completion and historic breakthrough of natural gas hydrate trial production in China′s offshore areas[J]. Marine Geology & Quaternary Geology, 2017, 37(5):34. http://hydz.chinajournal.net.cn/WKD/WebPublication/paperDigest.aspx?paperID=4f0c459a-7f39-4e5f-9697-dcf201702507
[17] Li Y, Jiang Z, Jiang S, et al. Heat flow and thermal evolution of a passive continental margin from shelf to slope-A case study on the Pearl River Mouth Basin, northern South China Sea[J]. Journal of Asian Earth Sciences, 2017.
[18] 宁伏龙, 蒋国盛, 汤凤林, 等.利用地热开采海底天然气水合物[J].天然气工业, 2006, 26(12):136-138. doi: 10.3321/j.issn:1000-0976.2006.12.038 NING Fulong, JIANG Guosheng, TANG Fenglin, et al.Utilizing geothermal energy to exploit marine gas hydrate[J]. Natureal Gas lndustry, 2006, 26(12):136-138. doi: 10.3321/j.issn:1000-0976.2006.12.038
[19] 窦斌, 秦明举, 蒋国盛, 等.利用地热开采南海天然气水合物的技术研究[J].海洋地质前沿, 2011(10):49-52. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=hydzdt201110007 DOU Bin, QIN Mingju, JIANG Guosheng, et al. A discussion on technology for gas hydrates production in the South China Sea using geothermal as an energy source[J].Marine Geology Frontiers, 2011(10):49-52. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=hydzdt201110007
[20] Liu Y, Hou J, Zhao H, et al. A method to recover natural gas hydrates with geothermal energy conveyed by CO2[J]. Energy, 2018, 144:265-278. doi: 10.1016/j.energy.2017.12.030
[21] 冯景春, 李小森, 王屹, 等.三维实验模拟双水平井联合法开采天然气水合物[J].现代地质, 2016, 30(4):929-936. doi: 10.3969/j.issn.1000-8527.2016.04.023 FENG Jingchun, LI Xiaosen, WANG Yi, et al.Three-dimension experimental investigation of hydrate dissociation by the combined method with dual horizontal wells[J]. Geoscience, 2016, 30(4):929-936. doi: 10.3969/j.issn.1000-8527.2016.04.023
[22] 李淑霞, 王炜, 陈月明, 等.多孔介质中天然气水合物注热开采影响因素实验研究[J].海洋地质前沿, 2011, 27(6):49-54. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=hydzdt201106008 LI Shuxia, WANG Wei, CHEN Yueming, et al. Experimental study on influence factors of hot-brine stimulation for dissociation of ngh in porous medium[J].Marine Geology Frontiers, 2011, 27(6):49-54. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=hydzdt201106008
[23] Dornan P, Alavi S, Woo T K. Free energies of carbon dioxide sequestration and methane recovery in clathrate hydrates[J]. Journal of Chemical Physics, 2007, 127(12):52. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=d4e9f3915f857fc6a5a8cf372591a699
[24] Vysniauskas A, Bishnoi P R. A kinetic study of methane hydrate formation[J]. Chemical Engineering Science, 1983, 38(7):1061-1072. doi: 10.1016/0009-2509(83)80027-X
[25] Anderson G K. Enthalpy of dissociation and hydration number of methane hydrate from the Clapeyron equation[J]. Journal of Chemical Thermodynamics, 2004, 36(12):1119-1127. doi: 10.1016/j.jct.2004.07.005
[26] 张伟, 梁金强, 何家雄, 等.南海北部神狐海域GMGS1和GMGS3钻探区天然气水合物运聚成藏的差异性[J].天然气工业, 2018, 38(3):138-149. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=trqgy201803018 ZHANG Wei, LIANG Jinqiang, HE Jiaxiong, et al. Differences in natural gas hydrate migration and accumulation between GMGS1 and GMGS3 drilling areas in the Shenhu area, northern South China Sea[J]. Natural Gas Industry, 2018, 38(3):138-149. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=trqgy201803018
[27] 万义钊, 吴能友, 胡高伟, 等.南海神狐海域天然气水合物降压开采过程中储层的稳定性[J].天然气工业, 2018, 38(4):117-128. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=trqgy201804015 WAN Yizhao, WU Nengyou, HU Gaowei, et al. Reservoir stability in the process of natural gas hydrate production by depressurization in the Shenhu area of the South China Sea[J]. Natural Gas Industry, 2018, 38(4):117-128. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=trqgy201804015
-
期刊类型引用(11)
1. 田振环，王厚杰，王威，史经昊. 海上地热能开发现状及其对中国的启示. 海洋地质前沿. 2024(06): 1-12 . 
百度学术
2. 王胜坡. 海洋可燃冰资源开发新技术现状及展望. 化学工程师. 2024(09): 71-75 . 百度学术
3. 陈野，赵笑寒，姜志晨，王佳惠，江鑫，卞琦，刘松. 基于SWOT的中国南海可燃冰资源商业开采前景浅析. 天然气与石油. 2024(06): 35-41 . 百度学术
4. 欧芬兰，于彦江，寇贝贝，陈靓. 水合物藏的类型、特点及开发方法探讨. 海洋地质与第四纪地质. 2022(01): 194-213 . 本站查看
5. 刘乐乐，万义钊，李承峰，张永超，刘昌岭，吴能友. 天然气水合物储层有效绝对渗透率现场测试进展. 海洋地质前沿. 2022(11): 40-55 . 百度学术
6. 王维希，张春生，吴颜雄，张审琴，夏晓敏. 联合深海地热开采天然气水合物技术展望. 现代化工. 2021(09): 17-21 . 百度学术
7. 冯轩，翟亚若，王久星，韩金虎，陈映赫，马麟. 置换法联合压裂开采天然气水合物技术. 现代化工. 2021(12): 22-26 . 百度学术
8. 吴能友，李彦龙，万义钊，孙建业，黄丽，毛佩筱. 海域天然气水合物开采增产理论与技术体系展望. 天然气工业. 2020(08): 100-115 . 百度学术
9. 卜庆涛，刘圣彪，胡高伟，刘昌岭，万义钊. 含水合物沉积物声学特性——实验模拟与数值模拟的对比分析. 海洋地质前沿. 2020(09): 56-67 . 百度学术
10. 陈强，胡高伟，李彦龙，万义钊，刘昌岭，吴能友，刘洋. 海域天然气水合物资源开采新技术展望. 海洋地质前沿. 2020(09): 44-55 . 百度学术
11. 韩笑，刘姝，万青翠，陈玲玲，张郑挥，李波. 热激法开采天然气水合物研究进展. 油气储运. 2019(08): 849-855 . 百度学术
其他类型引用(10)
计量
- 文章访问数: 2842
- HTML全文浏览量: 445
- PDF下载量: 30
- 被引次数: 21
