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.     

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.     

Injection of A Hydrate-Forming Gas into A Snow Layer Saturated with the Same Gas

The problem of injection of a hydrate-forming gas (methane) into a snow layer whose pores are initially saturated with the same gas is solved. Self-similar solutions describing the temperature and pressure fields and the snow, hydrate, and gas distributions in the layer are constructed. It is shown that, depending on the initial thermobaric state of the snow–methane system and the rate of gas injection, three characteristic zones can be distinguished in the filtration region: a near zone, in which snow is completely converted into hydrate and, consequently, the hydrate layer is saturated with gas; an intermediate zone, in which gas, snow, and hydrate are in phase equilibrium; far zone filled with gas and snow. It is shown that the length of the heated zone decreases with increasing initial snow content in the layer and with decreasing injected gas pressure. It is also shown that the length of the region of hydrate formation increases with increasing permeability. It is noted that the heating of the intermediate zone occurs more rapidly.     

沉积物中水合物形成机理及分解动力学研究进展

对沉积物中水合物形成机理的研究进展作了述评. 简述了纯水合物生成的机理及研究进展, 包括成核机理、铠甲效应、记忆效应等现象和机理;阐述了低可溶气体如甲烷生成水合物的难点, 以及主要因素如孔隙尺寸对沉积物中水合物生成的影响研究进展; 简要介绍了目前关于水合物形成和分解动力学的代表性数学模型;并对今后的研究方向进行了阐述.      

Breakdown of under-water gas hydrate deposits

At the floor and upper sediment layers of the submarine volcano Håkon Mosby gas (methane) hydrate deposits have been discovered. Measurements of the water-dissolved methane concentration at different depths near the volcano have shown that its level is higher than the average level in the Ocean. In order to understand the effects observed in concentration distribution, a model is proposed for describing the propagation of methane formed as a result of decomposition of depth and floor hydrate over the hydrosphere. Within the framework of physical modeling the kinetics of methane dissolution and its distribution in the hydrosphere are traced numerically with account for transfer anisotropy in different directions and undercurrents. The dissolution rate and the mass transfer directions are determined and the concentration fields in the floor layers are calculated. It is shown that methane is rapidly dispersed in the hydrosphere and, thus, significant volumes with high methane concentration cannot be formed. Practical recommendations are given for hunting hydrate deposits in water areas. The model proposed is applicable to describing the propagation of any admixture in different anisotropic media.     

Methane Hydrate and Two-Dimensional Fluid Transport Model: Comparison with Blake Ridge Chlorinity Measurements

The potential benefits and risks of natural methane hydrate occurrences have raised the need to understand the processes governing hydrate formation and dissociation over the last decades, and models of increasing complexity have been developed for this purpose. We propose a formulation of a multi-dimensional methane hydrate model that couples the established chemistry of hydrate formation to the fluid flow in marine sediments undergoing compaction. The numerical model applies the Finite Volume Method to construct a segregated solver for the coupled system. The solution of the sequence of individual processes is stabilised with an adaptive Picard iteration in each time step. We implement the model based on the OpenFOAM library and extend the functionality by generalising the diffusion term to take into account variable porosity of the sediment matrix. The advantage of this robust and efficient scheme is that the formulation is conservative by construction, allowing an accurate solution of mass transport when strong concentration gradients develop in the vicinity of phase boundaries. We validate the model using data from the Blake Ridge hydrate province. Good agreement is found for the pore water chlorinity, a proxy for hydrate formation, except for regions where heterogeneous hydrate formation results in highly variable measurements.     

Numerical investigation of the decomposition of gas hydrates coexisting with gas in natural reservoirs

The problem of methane hydrate decomposition in a reservoir saturated with a gas and hydrate mixture is investigated numerically. The results of the numerical simulation and an analytic solution obtained in the linear approximation are compared. It is shown that for high-permeability rocks the convective heat transfer in the near-well space of the reservoir predominates over the conductive transfer. This makes the use of intra-well heaters ineffective. It is found that an increase in the reservoir and well pressures and a decrease in the permeability suppress the formation of an extended hydrate dissociation region. Critical diagrams of existence of the frontal decomposition regime are constructed.     

Thermodynamic Conditions of Formation of CO<Subscript>2</Subscript> Hydrate in Carbon Dioxide Injection into a Methane Hydrate Reservoir

The mechanism of replacement of methane by carbon dioxide in the hydrate in the process of CO₂ injection into a reservoir with formation of fronts of methane hydrate dissociation and carbon dioxide hydrate generation is investigated. It is found that such a replacement regime can be implemented in both low- and high-permeability reservoirs. It is shown that in the highintensity injection regime the heat flux from the well does not affect propagation of the fronts of methane hydrate dissociation and carbon dioxide hydrate generation. In this case the replacement regime is maintained by only the heat released at formation of carbon dioxide hydrate. An increase in the injection pressure may lead to suppression of methane hydrate dissociation and termination of the replacement reaction. The critical diagrams of existence of the regime of conversion of methane hydrate to carbon dioxide hydrate are constructed.     

Formation of hydrate in injection of liquid carbon dioxide into a reservoir saturated with methane andwater

Injection of liquid carbon dioxide into a depleted natural gas field is investigated. A mathematical model of the process which takes into account forming CO₂ hydrate and methane displacement is suggested. An asymptotic solution of the problem is found in the one-dimensional approximation. It is shown that three injection regimes can exist depending on the parameters. In the case of weak injection, liquid carbon dioxide boils up with formation of carbon-dioxide gas. The intense regime is characterized by formation of CO₂ hydrate or a mixture of CO₂ and CH₄ hydrates. Critical diagrams of the process which determine the parameter ranges of the corresponding regimes are plotted.     

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

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

天然气水合物降压开采沉积物压缩效应数值模拟研究

深海天然气水合物降压开采过程中,沉积物的压缩会改变储层的物理力学特性,进而对天然气的开采效果产生显著影响.为揭示沉积物压缩效应下井周围储层物理力学特性演化规律,本文建立了考虑沉积物压缩效应的理论模型,通过COMSOL模拟研究了不同初始固有渗透率、初始水合物饱和度和井底压力条件下的降压开采中生产井周围储层的物理力学特性演化规律以及开采效果.结果表明:受沉积物压缩的影响,水合物分解区的渗透率随着与井筒距离的增加先增加后减少;产气与产水速率由零立即上升至峰值,然后迅速下降,并且考虑沉积物压缩时的产气与产水速率比不考虑时低;在水合物完全分解区,渗透率的大小与有效应力成负相关关系,未分解区渗透率的大小与水合物饱和度成负相关关系;井底压力越小,有效应力越大,生产井周围储层的渗透率下降越明显;初始水合物饱和度对产气与产水的影响存在拐点,饱和度拐点位于0.25与0.35之间,高水合物饱和度并不代表储层开采效果好,产气速率的高低还与储层的渗透率有关,高水合物饱和度储层的渗透率较低,产气速率较低;储层初始固有渗透率较高时显著促进了开采效果,但储层变形量较大增加了储层的不稳定性.     

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

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

含水合物粉质黏土压裂成缝特征实验研究

水力压裂技术是一种重要的油气井增产、增注措施,已经广泛应用于页岩油气等非常规资源的商业开采中.目前对于粉质黏土水合物沉积物的水力压裂成缝能力尚不清楚.本文采用南海水合物沉积层的粉质黏土制备沉积物试样,并与实验室配制的粉细砂土沉积物对比,分析粉质黏土沉积物的水力成缝能力及主控因素.实验结果表明含水合物和冰的沉积物破裂压力较高,这与粉质黏土沉积物特殊的应力-应变特征和渗透性有关.当沉积物应变高于6%时, 试样强度迅速上升, 呈现应变强化的特征,对水力拉伸裂缝的扩展具有一定的阻碍作用. 粉质黏土沉积物粒径细小, 渗透性差,难以通过渗透作用传递压力, 提高了沉积层的破裂压力. 此外,粉质黏土水合物沉积层裂缝扩展存在明显延迟效应,说明裂缝扩展受到流体压力和热应力的共同影响. 适当延长注入时间,保持流体与沉积层充分接触, 会起到分解水合物、降低破裂压力的作用.该研究成果有利于深入理解水力裂缝在水合物沉积层中的扩展规律,对探索压裂技术在水合物沉积层开发中的应用具有重要意义.      

Injection of liquid carbon dioxide into a reservoir partially saturated with methane hydrate

A mathematical model is proposed to describe methane–carbon dioxide replacement in gas hydrate by injecting liquid carbon dioxide into a porous medium initially saturated with methane and its hydrate. Self-similar solutions of the axisymmetric problems are constructed that describe the distribution of the main parameters of the reservoir. It is shown that there exist solutions according to which the process can occur both with and without boiling of carbon dioxide. Diagrams of the existence of each type of solution are constructed.     

Analytical solution of the nonlinear problem of gas hydrate dissociation in a formation

The problem of methane hydrate dissociation in a formation at high pressure gradients, when the flow in the near-well region is described by a nonlinear equation, is considered. In the quasistationary approximation, an analytical solution of the problem representable in the form of an implicit function is obtained. It is shown that during dissociation at high pressure gradients ice may be formed at fairly high initial temperatures. A critical diagram of hydrate dissociation regimes is plotted.     

水合物沉积物力学性质的实验装置和研究进展

天然气水合物是一种高效、洁净和储量巨大的新型能源,一般蕴含于砂岩、粘土以及其它土质的沉积物中.对水合物沉积物力学性质的实验研究,是水合物地层中基础稳定性分析和水合物开发评价重要的热点课题之一.本文首先介绍了水合物沉积物合成与分解实验、物性参数测量技术以及力学性质实验装置的主要组成部分和特点,然后介绍了目前国内外在水合物沉积物合成和分解及力学性质实验一体化装置和实验研究两个方面所取得的成果,最后指出在实验装置、测试技术和实验研究方面存在的问题以及今后研究的重点和方向.     

Diffusive dispersion of natural gas hydrates in ocean sediments

A mathematical model of the process of dissolving and dispersion of natural methane hydrates in ocean sediments is proposed. A condition which must be satisfied on the solution boundary of the hydrates in the region of their stable thermodynamic state is derived. The calculated solution rates are several centimeters per year. It is shown that the temperature perturbations initiated by the transition of a gas from the hydrate into the dissolved state are small and the process can be considered isothermal. Moscow. Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No. 1, pp. 170–173, January–February, 1999. The work was carried out with financial support from the Russian Foundation for Fundamental Research (project No. 96-01-00521).     

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.     

Effect of decomposition of a gas hydrate on the gas recovery from a reservoir containing hydrate and gas in the free state

The problem of gas recovery in the process of methane hydrate dissociation in a reservoir saturated with gas-hydrate mixture is considered. The mathematical model of hydrate decomposition into gas and water is generalized to include the negative temperature interval and takes ice formation into account. The solution of the problem is represented in the self-similar approximation. It is shown that there exists a transition hydrate decomposition regime in which water and ice are formed simultaneously. A comparative analysis of the recovery is carried out on the basis of relations derived for the masses of recovered gas in different hydrate dissociation regimes. It is shown that an anomalous increase in the recovered gas volumes is observed in the transition hydrate dissociation regime.__________Translated from Izvestiya Rossiiskoi Academii Nauk, Mekhanika Zhidkosti i Gaza, No. 1, 2005, pp. 132–142. Original Russian Text Copyright © 2005 by Tsypkin.     

石英砂粒径对水合物沉积物力学性质的影响

利用自主研发的水合物沉积物原位合成与力学性质测试的高压低温三轴仪,通过多级加荷的试验方法,以不同粒径的砂粒作为沉积物骨架进行三轴压缩试验,得到了剪切过程的应力-应变关系曲线,以及不同粒径尺寸沉积物的强度,还有剪切过程中的体积变化关系。结果表明:含水合物沉积物强度随着沉积物粒径尺寸的增大而增强;在降压剪切过程中,所有粒径的水合物沉积物式样均有明显的剪缩现象。