不同水合物饱和度下注热水开采实验研究

为研究不同水合物藏饱和度对注热开采动态的影响,采用自制的一维天然气水合物(NGH)开采模拟实验装置,模拟地层多孔介质的条件,首先在填砂模型中生成不同饱和度的NGH,然后进行注热盐水分解实验。结果表明:不同饱和度的NGH注热分解产气都可划分为三个阶段,不同的是NGH饱和度越高,水合物分解阶段的产气速率越大,且该阶段持续的时间越长;NGH饱和度越高,注热分解阶段电阻率增大的幅度越大,系统各测点温度升高的幅度越小。注热分解过程中产水速率围绕注水速率而上下波动,且NGH饱和度越高,产水速率波动幅度越大;在实验研究范围内,随初始水合物饱和度的升高,注热开采的能量效率逐渐升高。因此,从能量效率来说,高饱和度的水合物藏更适宜于注热开采。     

海洋天然气水合物降压开采地层井壁力学稳定性分析

降压开采天然气水合物使其分解会导致储层孔隙度、渗透率、孔隙压力和岩层骨架有效应力发生改变, 同时降低沉积物的胶结程度, 使地层的抗剪强度和承载能力降低, 从而引起井壁失稳、海底滑坡、海底面沉降等工程问题. 为此, 在地下多相非等温数值模拟软件TOUGH+Hydrate框架内, 基于扩展的三维Biot固结理论, 考虑水合物分解相变、传热(T)、流动(H)、岩土体变形(M)等过程及其相互耦合作用, 建立了新的水合物开采传热-流动-力学(THM)耦合数学模型, 并开发有限元程序对其进行数值求解. 以中国南海神狐海域GMGS1航次SH2站位水合物储层条件为研究对象, 构建了垂直井降压开采THM耦合地层井壁稳定性分析模型, 预测了水合物开采过程中储层温-压-力场和水合物分解区的演化规律, 揭示了地层优势出砂区域和海底面沉降趋势. 结果表明: 储层降压导致地层有效应力增大, 进而引起井周地层发生沉降, 且地层的沉降主要发生在降压开采前期, 最大沉降位置位于井壁周围, 向储层内部延伸地层沉降量快速减小; 水合物分解导致井周地层力学强度降低, 加剧了储层的沉降; 井筒降压造成射孔段井壁应力集中最为明显, 从而造成井壁破坏的潜在风险, 这些区域正是水合物开采出砂防治的关键区域.      

海洋天然气水合物降压开采地层井壁力学稳定性分析

降压开采天然气水合物使其分解会导致储层孔隙度、渗透率、孔隙压力和岩层骨架有效应力发生改变,同时降低沉积物的胶结程度,使地层的抗剪强度和承载能力降低,从而引起井壁失稳、海底滑坡、海底面沉降等工程问题.为此,在地下多相非等温数值模拟软件TOUGH+Hydrate框架内,基于扩展的三维Biot固结理论,考虑水合物分解相变、传热(T)、流动(H)、岩土体变形(M)等过程及其相互耦合作用,建立了新的水合物开采传热-流动-力学(THM)耦合数学模型,并开发有限元程序对其进行数值求解.以中国南海神狐海域GMGS1航次SH2站位水合物储层条件为研究对象,构建了垂直井降压开采THM耦合地层井壁稳定性分析模型,预测了水合物开采过程中储层温-压-力场和水合物分解区的演化规律,揭示了地层优势出砂区域和海底面沉降趋势.结果表明:储层降压导致地层有效应力增大,进而引起井周地层发生沉降,且地层的沉降主要发生在降压开采前期,最大沉降位置位于井壁周围,向储层内部延伸地层沉降量快速减小;水合物分解导致井周地层力学强度降低,加剧了储层的沉降;井筒降压造成射孔段井壁应力集中最为明显,从而造成井壁破坏的潜在风险,这些区域正是水合物开采出砂防治的关键区域.     

降压开采模拟试验的水合物分解阵面演化过程

水合物分解阵面演化过程与开采安全性和产气效率密切相关,是开采原位监测的重要组成部分。在玻璃砂样品中进行了甲烷水合物降压开采模拟试验,探讨了水合物饱和度对渗流阵面和水合物分解阵面演化过程的影响,结合已有理论模型,分析了水合物分解阵面传播速度的关键影响因素。结果表明:渗流阵面和水合物分解阵面的传播距离均与时间平方根呈近似线性关系,传播速度均随水合物饱和度的增加而减小;水合物分解阵面的传播速度随多孔介质的有效渗透率和降压幅度的增加而变快,随孔隙率的增加而变慢,粗砂质地层更有利于水合物降压分解阵面的传播。     

南海水合物黏土沉积物力学特性试验模拟研究

利用自行研制的含水合物沉积物合成、分解与力学性质测量一体化试验设备,以南海水合物区域的海底粉质黏土作为骨架,制备含水合物沉积物样品,并对其进行了三轴压缩试验研究,获得了水合物分解前后的应力应变曲线和抗剪强度特性. 结果表明:在水合物饱和度0%~45% 的范围内,水合物沉积物的应力应变曲线均表现为弹塑性变形,存在明显的应变硬化现象;抗剪强度、内摩擦角和黏聚力随水合物饱和度的增加而增加. 相对而言,内摩擦角随饱和度增加幅度较小,其他参数在水合物饱和度超过25% 时,呈陡然增高趋势;水合物分解后导致抗剪强度最大可降低为初始的1/4,不同初始饱和度条件下水合物完全分解后沉积物的抗剪强度基本相等,并大于同等围压条件下初始不含水合物的沉积物抗剪强度.      

含水合物沉积物三轴剪切试验与损伤统计分析

天然气水合物开采诱发水合物分解,削弱水合物地层强度,可能导致地层滑动和生产平台倒塌等工程地质灾害,对水合物开采安全性构成严重威胁.深入理解含水合物沉积物力学性质并建立合理的本构关系模型是水合物开采安全性评价的前提条件.在自主研发的含水合物沉积物力学性质测试实验装置上,采用饱和水海砂沉积物气体扩散法制备了含水合物沉积物样品,并开展了系列的排水三轴剪切试验,通过时域反射技术实现了样品中水合物饱和度的实时在线测量;基于复合材料的罗伊斯(Reuss)应力串联模型和沃伊特(Voigt)应变并联模型提出了含水合物沉积物等效弹性模量的细观力学混合律模型,结合损伤统计理论和摩尔-库伦破坏准则改进了含水合物沉积物的本构关系模型.结果表明:随着水合物饱和度的增加和有效围压的减小,应力-应变曲线由应变硬化型变为应变软化型,割线模量和峰值强度均随水合物饱和度与有效围压的增加而提高,黏聚力受水合物饱和度影响明显,而内摩擦角基本不变;提出的等效弹性模量细观力学混合律模型与改进的本构关系模型均具有良好的适用性,模型参数少且物理意义明确.      

含水合物沉积物多相渗流特性研究进展

天然气水合物作为一种储量大、无污染的清洁能源近些年受到了广泛关注.近20年来,中国进行了较大范围的陆海域天然气水合物储层勘探与储量预测. 2017年,中国地质调查局牵头对南海神狐海域的天然气水合物进行了基于降压渗流原理的试验性开采.国内外已进行的水合物试采工程面临着气体产量低、出砂较多等问题,其最主要的原因之一是开发过程中沉积物内复杂多相渗流机理尚不明晰.本文综述了平行毛细管模型、Kozeny模型等广泛应用于天然气水合物开发渗流分析的理论模型,对比分析了水合物开发多尺度渗流过程模拟方法,简述了国内外含水合物沉积物渗透率测试、渗流过程中沉积物物性演变以及水合物开采室内模拟等方面的渗流实验进展,总结了矿场尺度的天然气水合物储层开采过程中产气数值模拟手段,展望了多相渗流模型、储层原位含水合物样品室内测试及结构与物性演化、矿场尺度数值模拟与水平井压裂技术等应用研究的未来方向与挑战.     

含水合物沉积物多相渗流特性研究进展

天然气水合物作为一种储量大、无污染的清洁能源近些年受到了广泛关注. 近20年来,中国进行了较大范围的陆海域天然气水合物储层勘探与储量预测.2017年,中国地质调查局牵头对南海神狐海域的天然气水合物进行了基于降压渗流原理的试验性开采.国内外已进行的水合物试采工程面临着气体产量低、出砂较多等问题,其最主要的原因之一是开发过程中沉积物内复杂多相渗流机理尚不明晰.本文综述了平行毛细管模型、Kozeny模型等广泛应用于天然气水合物开发渗流分析的理论模型,对比分析了水合物开发多尺度渗流过程模拟方法,简述了国内外含水合物沉积物渗透率测试、渗流过程中沉积物物性演变以及水合物开采室内模拟等方面的渗流实验进展,总结了矿场尺度的天然气水合物储层开采过程中产气数值模拟手段,展望了多相渗流模型、储层原位含水合物样品室内测试及结构与物性演化、矿场尺度数值模拟与水平井压裂技术等应用研究的未来方向与挑战.      

有效应力升降过程中不同饱和度水合物沉积物渗透率实验研究

为了研究天然气水合物开采过程中水合物饱和度和有效应力变化对含水合物沉积层渗透率变化规律的影响,利用自主研制的水合物沉积物合成与三轴渗流实验一体化装置,进行了有效应力升降过程中不同饱和度水合物沉积物渗透性实验。结果表明:有效应力升降过程中,不同饱和度水合物沉积物的渗透率与有效应力均呈负指数规律变化,并表现出在低有效应力阶段渗透率变化的幅度大于高有效应力阶段;有效应力升降过程中,试样产生的不可恢复变形使渗透率不能完全恢复,造成渗透率永久损失;有效应力不变时,水合物沉积物的渗透率与水合物饱和度呈负指数规律变化,且曲线的斜率随饱和度增加由大变小;水合物的饱和度越大,最大渗透率损害率越大,渗透率恢复的程度越差。     

天然气水合物开发多物理场特征及耦合渗流研究进展与建议

天然气水合物作为一种非常规的清洁能源,在全球分布广、资源量大.自20世纪90年代以来,加拿大、美国、日本、中国已经先后进行了陆域及海域的水合物试采,但发现出砂、单井日产气量低、稳产时间短等问题,试采产量远不能满足商业化开发的需求,其中核心问题是对水合物开发过程中的相变、多相多组分多场耦合渗流特征的认识不够明晰.本文根据天然气水合物开发过程中涉及的渗流场、温度场、化学场、力学场等多场耦合特征,重点综述水合物生成/分解对各物理场主要特征参数的影响,包括水合物储层的孔隙度、水合物饱和度、渗透率、相对渗透率等基础物性参数及其动态演变,天然气水合物的导热系数、比热容、热扩散系数以及水合物生成/分解热等热力学参数,天然气水合物生成、分解动力学特征,纯水合物以及含水合物沉积物的力学性质等,最后阐述了天然气水合物开发渗流中的多场耦合关系及相互作用,提出了今后水合物开发多物理场特征及耦合渗流的科学研究、技术开发的有关建议.     

The Experimental and Numerical Studies on Gas Production from Hydrate Reservoir by Depressurization

A set of experimental system to study hydrate dissociation in porous media is built and some experiments on hydrate dissociation by depressurization are carried out. A mathematical model is developed to simulate the hydrate dissociation by depressurization in hydrate-bearing porous media. The model can be used to analyze the effects of the flow of multiphase fluids, the kinetic process and endothermic process of hydrate dissociation, ice-water phase equilibrium, the variation of permeability, convection and conduction on the hydrate dissociation, and gas and water productions. The numerical results agree well with the experimental results, which validate our mathematical model. For a 3-D hydrate reservoir of Class 3, the evolutions of pressure, temperature, and saturations are elucidated and the effects of some main parameters on gas and water rates are analyzed. Numerical results show that gas can be produced effectively from hydrate reservoir in the first stage of depressurization. Then, methods such as thermal stimulation or inhibitor injection should be considered due to the energy deficiency of formation energy. The numerical results for 3-D hydrate reservoir of Class 1 show that the overlying gas hydrate zone can apparently enhance gas rate and prolong life span of gas reservoir.     

天然气水合物开发多物理场特征及耦合渗流研究进展与建议

天然气水合物作为一种非常规的清洁能源, 在全球分布广、资源量大. 自20世纪90年代以来, 加拿大、美国、日本、中国已经先后进行了陆域及海域的水合物试采, 但发现出砂、单井日产气量低、稳产时间短等问题, 试采产量远不能满足商业化开发的需求, 其中核心问题是对水合物开发过程中的相变、多相多组分多场耦合渗流特征的认识不够明晰. 本文根据天然气水合物开发过程中涉及的渗流场、温度场、化学场、力学场等多场耦合特征, 重点综述水合物生成/分解对各物理场主要特征参数的影响, 包括水合物储层的孔隙度、水合物饱和度、渗透率、相对渗透率等基础物性参数及其动态演变, 天然气水合物的导热系数、比热容、热扩散系数以及水合物生成/分解热等热力学参数, 天然气水合物生成、分解动力学特征, 纯水合物以及含水合物沉积物的力学性质等, 最后阐述了天然气水合物开发渗流中的多场耦合关系及相互作用, 提出了今后水合物开发多物理场特征及耦合渗流的科学研究、技术开发的有关建议.      

The simulation of gas production from oceanic gas hydrate reservoir by the combination of ocean surface warm water flooding with depressurization

A new method is proposed to produce gas from oceanic gas hydrate reservoir by combining the ocean surface warm water flooding with depressurization which can efficiently utilize the synthetic effects of thermal, salt and depressurization on gas hydrate dissociation. The method has the advantage of high efficiency, low cost and enhanced safety. Based on the proposed conceptual method, the physical and mathematical models are established, in which the effects of the flow of multiphase fluid, the kinetic process of hydrate dissociation, the endothermic process of hydrate dissociation, ice-water phase equilibrium, salt inhibition, dispersion, convection and conduction on the hydrate dissociation and gas and water production are considered. The gas and water rates, formation pressure for the combination method are compared with that of the single depressurization, which is referred to the method in which only depressurization is used. The results show that the combination method can remedy the deficiency of individual producing methods. It has the advantage of longer stable period of high gas rate than the single depressurization. It can also reduce the geologic hazard caused by the formation deformation due to the maintaining of the formation pressure by injected ocean warm water.     

Modeling of Gas Hydrate Formation Upon Injection of Carbon Dioxide Gas into A Natural Reservoir

This paper presents the results of numerical modeling of gas hydrate formation upon injection of carbon dioxide into a finite-length reservoir saturated with methane and water. It is shown that at different stages, hydrate formation can occur on both the frontal surface and in a reservoir region of finite length. The effects of pressure at the reservoir boundaries and the effects of the permeability and initial water saturation of the reservoir on the hydrate formation process were studied. The dependences of the time of the complete conversion of water into gas hydrate in the entire reservoir on the injection pressure and reservoir permeability were obtained.     

水平井降压法和热激法水合物开采对海底边坡稳定性的影响

天然气水合物广泛赋存在深海沉积物孔隙中, 被认为是具有巨大开发潜力的未来绿色能源之一, 引起全球的关注. 深海水合物开采将造成含水合物储层的强度劣化, 可能产生孔压积聚, 诱发海床失稳. 本文基于边坡稳定极限平衡分析框架, 引入考虑水合物开采热-流-化学耦合过程的数值分析模型, 研究水合物开采对海底边坡稳定性的影响. 采用TOUGH+HYDRATE热-流-化学耦合分析程序, 模拟了采用水平井降压法和热激法开采深海水合物的过程, 分析了水合物分解锋面扩展和瞬态孔压演变的规律, 并通过SLOPE/W程序采用极限平衡分析方法计算水合物开采过程及停采后的海底边坡安全系数, 分析开采井位置和开采方法对海底边坡稳定性的影响. 研究表明, 对于存在致密盖层的细砂储层陡坡, 单水平井降压开采过程中, 由于孔压降低, 土体有效应力增加, 边坡稳定性显著提高, 当开采井布设在坡体中部时, 边坡稳定性提高最为明显; 停采后, 由于水合物分解导致土体黏聚强度降低, 且孔压逐渐回升到静水压状态, 导致边坡稳定性下降, 最危险滑弧通过水合物分解区. 若采用双水平井热激法开采, 开采过程与停采后的最危险滑弧始终通过水合物分解区, 由于开采过程中温度升高, 井周孔压显著上升, 导致边坡安全系数明显下降, 存在诱发滑坡的风险.      

Injection of carbon dioxide into a reservoir saturated with methane and water

This paper presents a mathematical model for the injection of carbon dioxide into a natural gas reservoir saturated with methane and water accompanied by the formation of carbon dioxide hydrate in an extended region. The dependence of the coordinates of the boundaries of the region of phase transitions on the pressure of the injected gas and the initial parameters of the reservoir are investigated. It is established that the velocity of the near boundary of the region of hydrate formation decreases with increasing water saturation and initial temperature of the reservoir and the velocity of the far boundary of the region of phase transitions increases with increasing pressure of the injected gas and reservoir permeability. It is shown that at high initial temperatures of the reservoir, a regime is possible in which replacement of methane by carbon dioxide without hydrate formation occurs at the far interface, and at the near interface, water is completely incorporated into gas hydrate.