2-D Numerical Simulation of Natural Gas Hydrate Decomposition Through Depressurization by Fully Implicit Method

Natural gas hydrate, as a potential energy resource, deposits in permafrost and marine sediment with large quantities. The current exploitation methods include depressurization, thermal stimulation, and inhibitor injection. However, many issues have to be resolved before the commercial production. In the present study, a 2-D axisymmetric simulator for gas production from hydrate reservoirs is developed. The simulator includes equations of conductive and convective heat transfer, kinetic of hydrate decomposition, and multiphase flow. These equations are discretized based on the finite difference method and are solved with the fully implicit simultaneous solution method. The process of laboratory-scale hydrate decomposition by depressurization is simulated. For different surrounding temperatures and outlet pressures, time evolutions of gas and water generations during hydrate dissociation are evaluated, and variations of temperature, pressure, and multiphase fluid flow conditions are analyzed. The results suggest that the rate of heat transfer plays an important role in the process. Furthermore, high surrounding temperature and low outlet valve pressure may increase the rate of hydrate dissociation with insignificant impact on final cumulative gas volume.     

热水驱替开采Ⅱ类水合物藏规律研究

天然气水合物广泛分布于陆地冻土带和深海地层,资源潜力巨大,其中Ⅱ类水合物藏占有重要地位。为加强对Ⅱ类水合物藏开采规律的认识,结合实际水合物藏参数,使用数值模拟方法研究了热水驱替开采Ⅱ类水合物藏的动态规律,并与降压法的开采效果进行了对比分析。结果表明:①热水驱替开采Ⅱ类水合物藏时,产气速率和分解气速率首先快速上升,然后以较快速度下降至趋于相对稳定;累产气和累分解气上升较快;气体采出程度和水合物分解程度均处于较高水平(>60%)。②热水驱替对Ⅱ类水合藏的开采具有一定的适应性,与降压法开采相比,热水驱替方式下储层水合物的分解更彻底,气体采出程度、水合物分解程度也更优,但具有较低的累积气水比,产水量较大。     

天然气水合物稳定带内流体压裂计算的程序耦合方法

海洋天然气水合物稳定带气烟囱结构中存在被水合物充填的裂隙, 表明在自然条件下沉积物中曾发生过流体压裂以及相关的流体流动和水合物形成。在水合物稳定带内实施人为的流体压裂工程, 并联合其他方法(如降压或注热)进行水合物开采, 有望提高开采效率。水合物稳定带内, 无论是自然条件下发生的流体压裂过程, 还是人为实施的流体压裂工程, 都存在水合物反应和沉积物裂隙变形之间的耦合响应。当前, 已有不少数值程序对水合物反应与沉积物弹塑性变形的耦合过程进行了定量研究, 但尚没有数值程序能够计算水合物反应和离散裂隙变形之间的耦合过程。文章将TOUGH+Hydrate程序、IC-FERST和Solidity两者的耦合程序进行了进一步耦合, 为水合物稳定带内的流体压裂计算提供了一种耦合计算方法, 同时通过一个算例初步验证了该耦合计算方法的可行性。验证结果表明, 该耦合计算方法经进一步改进后有望应用于定量研究水合物稳定带内的裂隙变形和水合物反应过程。     

海底天然气水合物开发的前景研究——基于经济可行性的视角

海底天然气水合物的开发利用对于国家能源安全与环境保护具有重要意义。然而,作为一种新型深海能源,海底天然气水合物至今仍未实现商业化生产。除开发技术等因素外,天然气水合物开发投入、产出的成本问题也是制约海底天然气水合物大规模开发的另一关键因素。文章从经济可行性的视角分析海底天然气水合物目前开采成本、未来开采成本的变化以及未来的价格优势,从而对天然气水合物开发的前景进行相关研究。研究结果显示,虽然目前天然气水合物的开采成本要高于常规天然气,但随着今后开采技术的改进,天然气水合物开采的总成本、边际成本和平均成本都会相应降低,当天然气水合物的价格低于常规天然气价格时,天然气水合物的商业化生产便可能实现。     

Numerical studies of methane production from Class 1 gas hydrate accumulations enhanced with carbon dioxide injection

Class 1 gas hydrate accumulations are characterized by a permeable hydrate-bearing interval overlying a permeable interval with mobile gas, sandwiched between two impermeable intervals. Depressurization-induced dissociation is currently the favored technology for producing gas from Class 1 gas hydrate accumulations. The depressurization production technology requires heat transfer from the surrounding environment to sustain dissociation as the temperature drops toward the hydrate equilibrium point and leaves the reservoir void of gas hydrate. Production of gas hydrate accumulations by exchanging carbon dioxide with methane in the clathrate structure has been demonstrated in laboratory experiments and proposed as a field-scale technology. The carbon dioxide exchange technology has the potential for yielding higher production rates and mechanically stabilizing the reservoir by maintaining hydrate saturations. We used numerical simulation to investigate the advantages and disadvantages of using carbon dioxide injection to enhance the production of methane from Class 1 gas hydrate accumulations. Numerical simulations in this study were primarily concerned with the mechanisms and approaches of carbon dioxide injection to investigate whether methane production could be enhanced through this approach. To avoid excessive simulation execution times, a five-spot well pattern with a 500-m well spacing was approximated using a two-dimensional domain having well boundaries on the vertical sides and impermeable boundaries on the horizontal sides. Impermeable over- and under burden were included to account for heat transfer into the production interval. Simulation results indicate that low injection pressures can be used to reduce secondary hydrate formation and that direct contact of injected carbon dioxide with the methane hydrate present in the formation is limited due to bypass through the higher permeability gas zone.     

腰英台气田天然气水合物防治

松辽盆地南部长岭断陷腰英台深层构造腰深1井井筒及输气管线发生水合物堵塞现象,分析认为天然气水合物是导致堵塞的原因。天然气水合物是在一定温度和压力条件下由天然气中烃分子与游离水结合而形成。计算出腰深1井水合物添加甲醇的合理量,通过从油管和套管添加甲醇,使该井恢复了正常生产,日产气量（6～12）×10^4m^3。防堵的关键是注入甲醇和连续平稳生产。     

海洋天然气水合物开采方法及产量分析

海洋天然气水合物的巨大储量刺激了世界各国能源部门努力研究如何从天然气水合物储层生产天然气。根据水合物形成的条件,只有当水合物处在其相平衡条件以外,水合物才能分解。因此,水合物的开采方法只能为热熔法、抑制剂刺激法、减压法和地面分解法。为了对天然气水合物储层中气体的生产有个定量的评估,本文以水合物开采井为例,运用数学方法推导了水合物井中气体的产生量。结果表明,在天然气水合物储层中,天然气释放量是井内水合物分解温度、压力及水合物层气体渗透性的敏感函数。该函数可以用于天然气水合物井气体开采量的计算及对水合物储层可开采性评价。     

Permafrost-associated natural gas hydrate occurrences on the Alaska North Slope

In the 1960s Russian scientists made what was then a bold assertion that gas hydrates should occur in abundance in nature. Since this early start, the scientific foundation has been built for the realization that gas hydrates are a global phenomenon, occurring in permafrost regions of the arctic and in deep water portions of most continental margins worldwide. In 1995, the U.S. Geological Survey made the first systematic assessment of the in-place natural gas hydrate resources of the United States. That study suggested that the amount of gas in the gas hydrate accumulations of northern Alaska probably exceeds the volume of known conventional gas resources on the North Slope. Researchers have long speculated that gas hydrates could eventually become a producible energy resource, yet technical and economic hurdles have historically made gas hydrate development a distant goal. This view began to change in recent years with the realization that this unconventional resource could be developed with existing conventional oil and gas production technology. One of the most significant developments was the completion of the BPXA-DOE-USGS Mount Elbert Gas Hydrate Stratigraphic Test Well on the Alaska North Slope, which along with the Mallik project in Canada, have for the first time allowed the rational assessment of gas hydrate production technology and concepts. Almost 40 years of gas hydrate research in northern Alaska has confirmed the occurrence of at least two large gas hydrate accumulations on the North Slope. We have also seen in Alaska the first ever assessment of how much gas could be technically recovered from gas hydrates. However, significant technical concerns need to be further resolved in order to assess the ultimate impact of gas hydrate energy resource development in northern Alaska.     

南海神狐海域天然气水合物饱和度的数值模拟分析

珠江口盆地神狐海域是天然气水合物钻探和试验开采的重点区域,大量钻探取心、测井与地震等综合分析表明不同站位水合物的饱和度、厚度与气源条件存在差异。本文利用天然气水合物调查及深水油气勘探所采集的测井和地震资料建立地质模型,利用PetroMod软件模拟地层的温度场、有机质成熟度、烃源岩生烃量、流体运移路径以及不同烃源岩影响下的水合物饱和度,结果表明:生物成因气分布在海底以下1500 m范围内的有机质未成熟地层,而热成因气分布在深度超过2300 m的成熟、过成熟地层。水合物稳定带内生烃量难以形成水合物,形成水合物气源主要来自于稳定带下方向上运移的生物与热成因气。模拟结果与测井结果对比分析表明,稳定带下部生物成因气能形成的水合物饱和度约为10%,在峡谷脊部的局部区域饱和度较高;相对高饱和度(>40%)水合物形成与文昌组、恩平组的热成因气沿断裂、气烟囱等流体运移通道幕式释放密切相关,W19井形成较高饱和度水合物的甲烷气体中热成因气占比达80%,W17井热成因气占比为73%,而SH2井主要以生物成因为主,因此,不同站位甲烷气体来源占比不同。     

Investigation of the hydrate reservoir production under different depressurization modes

Abstract

Large reserves of natural gas hydrates exist, and the depressurization method has the greatest potential for gas hydrate reservoir recovery. Currently, the most commonly adopted depressurization simulation method is a constant bottom-hole pressure production scheme. This study proposes a new depressurization mode with decreasing bottom-hole pressure. The production characteristic was numerically investigated using this method. The results show the following: (1) As the depressurization exponent (n) decreases, the development effect improves, and production indexes including cumulative gas production/dissociation and gas-water ratio increase. However, the reservoir energy consumption is higher and the hydrate reformation is more severe. (2) Compared to the proposed depressurization mode, the hydrate production index of the constant bottom-hole pressure production (n?=?0) is better. However, the hydrate reservoir energy consumption is higher and the hydrate reformation is more severe using constant bottom-hole pressure production. (3) To achieve a balance between production and reservoir energy consumption during depressurization production, the bottom-hole pressure should be controlled by selecting a suitable depressurization exponent between n_min and n_max, which can be determined through numerical simulations.     

Gas hydrate within the Winona Basin, offshore western Canada

In western Canada gas hydrates have been thought to exist primarily in the Cascadia accretionary prism off southern Vancouver Island, British Columbia (BC). We present evidence for the existence of gas hydrate in folds and ridges of the Winona Basin up to 40 km seaward from the foot of the continental slope off northern Vancouver Island. The occurrence of a bottom-simulating reflector (BSR) observed in a number of vintage seismic reflection profiles is strongly correlated to faulted, and folded sedimentary ridges and buried folds. The observed tectonic structures of the Winona Basin are within the rapidly evolving Juan de Fuca - Cascadia - Queen Charlotte triple junction off BC. Re-processing of multi-channel data imaged mildly to strongly deformed sediments; the BSR is confined to sediments with stronger deformation. Changes in the amplitude character of sediment-reflections above and below the depth of the base of gas hydrate stability zone were also used as an indicator for the presence of gas hydrate. Additionally, regional amplitude and frequency reduction below some strong BSR occurrences may indicate free gas accumulations. Gas hydrate formation in the Winona Basin appears strongly constrained to folds and ridges and thus correlated to deeper-routed fluid-advection regimes. Methane production from in situ microbial activities as a source of gas to form gas hydrates, as proposed to be a major contributor for gas hydrates within the accretionary prism to the south, appears to be insufficient to produce the widespread gas hydrate occurrences in the Winona Basin. Potential reasons for the lack of sufficient in situ gas production may be that sedimentation rates are 5-100 times higher than those in the accretionary prism so that available organic carbon moves too quickly through the gas hydrate stability field. The confinement of BSRs to ridges and folds within the Winona Basin results in an areal extent of gas hydrate occurrences that is a factor of five less than what is expected from regional gas hydrate stability field mapping using water-depth (pressure) as the only controlling factor only.     

Pre- and post-drill comparison of the Mount Elbert gas hydrate prospect, Alaska North Slope

In 2006, the United States Geological Survey (USGS) completed a detailed analysis and interpretation of available 2-D and 3-D seismic data, along with seismic modeling and correlation with specially processed downhole well log data for identifying potential gas hydrate accumulations on the North Slope of Alaska. A methodology was developed for identifying sub-permafrost gas hydrate prospects within the gas hydrate stability zone in the Milne Point area. The study revealed a total of 14 gas hydrate prospects in this area.In order to validate the gas hydrate prospecting protocol of the USGS and to acquire critical reservoir data needed to develop a longer-term production testing program, a stratigraphic test well was drilled at the Mount Elbert prospect in the Milne Point area in early 2007. The drilling confirmed the presence of two prominent gas-hydrate-bearing units in the Mount Elbert prospect, and high quality well logs and core data were acquired. The post-drill results indicate pre-drill predictions of the reservoir thickness and the gas-hydrate saturations based on seismic and existing well data were 90% accurate for the upper unit (hydrate unit D) and 70% accurate for the lower unit (hydrate unit C), confirming the validity of the USGS approach to gas hydrate prospecting. The Mount Elbert prospect is the first gas hydrate accumulation on the North Slope of Alaska identified primarily on the basis of seismic attribute analysis and specially processed downhole log data. Post-drill well log data enabled a better constraint of the elastic model and the development of an improved approach to the gas hydrate prospecting using seismic attributes.     

Numerical simulations of depressurization-induced gas production from gas hydrate reservoirs at the Walker Ridge 313 site,northern Gulf of Mexico

In 2009, the Gulf of Mexico (GOM) Gas Hydrates Joint-Industry-Project (JIP) Leg II drilling program confirmed that gas hydrate occurs at high saturations within reservoir-quality sands in the GOM. A comprehensive logging-while-drilling dataset was collected from seven wells at three sites, including two wells at the Walker Ridge 313 site. By constraining the saturations and thicknesses of hydrate-bearing sands using logging-while-drilling data, two-dimensional (2D), cylindrical, r-z and three-dimensional (3D) reservoir models were simulated. The gas hydrate occurrences inferred from seismic analysis are used to delineate the areal extent of the 3D reservoir models. Numerical simulations of gas production from the Walker Ridge reservoirs were conducted using the depressurization method at a constant bottomhole pressure. Results of these simulations indicate that these hydrate deposits are readily produced, owing to high intrinsic reservoir-quality and their proximity to the base of hydrate stability. The elevated in situ reservoir temperatures contribute to high (5–40 MMscf/day) predicted production rates. The production rates obtained from the 2D and 3D models are in close agreement. To evaluate the effect of spatial dimensions, the 2D reservoir domains were simulated at two outer radii. The results showed increased potential for formation of secondary hydrate and appearance of lag time for production rates as reservoir size increases. Similar phenomena were observed in the 3D reservoir models. The results also suggest that interbedded gas hydrate accumulations might be preferable targets for gas production in comparison with massive deposits. Hydrate in such accumulations can be readily dissociated due to heat supply from surrounding hydrate-free zones. Special cases were considered to evaluate the effect of overburden and underburden permeability on production. The obtained data show that production can be significantly degraded in comparison with a case using impermeable boundaries. The main reason for the reduced productivity is water influx from the surrounding strata; a secondary cause is gas escape into the overburden. The results dictate that in order to reliably estimate production potential, permeability of the surroundings has to be included in a model.     

Production potential and stability of hydrate-bearing sediments at the site GMGS3-W19 in the South China Sea: A preliminary feasibility study

According to the preliminary geological data of gas hydrate bearing-sediments (GHBS) at site GMGS3-W19 in the third Chinese expedition to drill gas hydrates in 2015, a production model using three different recovery pressures was established to assess the production feasibility from both production potential and geomechanical response. The simulation results show that for this special Class 1 deposit, it is a little hard for gas production rate to reach the commercial extraction rate because the degree of hydrate dissociation is limited due to the low reservoir permeability and the permeable burdens. However, the free gas accumulating in the lower part of the GHBS can significantly increase gas-to-water ratio. It also generates many secondary hydrates in the GHBS at the same time. Decreasing the well pressure can be beneficial to gas recovery, but the recovery increase is not obvious. In term of geomechanical response of the reservoir during the gas recovery, the permeable burdens are conducive to reduction of the sediment deformation, though they don't facilitate the gas recovery rate. In addition, significant stress concentration is observed in the upper and lower edges of GHBS around the borehole during depressurization because of high pressure gradient, and the greater the well pressure drop, the more obvious the phenomenon. Yield failures and sand production easily take place in the edges. Therefore, in order to achieve the purpose of safe, efficient and long-term gas production, a balance between the production pressure and reservoir stability should be reached at the hydrate site. The production pressure difference and sand production must be carefully controlled and the high stress concentration zones need strengthening or sand control treatment during gas production. Besides, the sensitivity analyses show that the hydrate saturation heterogeneity can affect the production potential and geomechanical response to some extent, especially the water extraction rate and the effective stress distribution and evolution. Increasing GHBS and its underlying free gas formation permeabilities can enhance the gas production potential, but it probably introduces geomechanical risks to gas recovery operations.     

In-situ gas hydrate hydrate saturation estimated from various well logs at the Mount Elbert Gas Hydrate Stratigraphic Test Well, Alaska North Slope

In 2006, the U.S. Geological Survey (USGS) completed detailed analysis and interpretation of available 2-D and 3-D seismic data and proposed a viable method for identifying sub-permafrost gas hydrate prospects within the gas hydrate stability zone in the Milne Point area of northern Alaska. To validate the predictions of the USGS and to acquire critical reservoir data needed to develop a long-term production testing program, a well was drilled at the Mount Elbert prospect in February, 2007. Numerous well log data and cores were acquired to estimate in-situ gas hydrate saturations and reservoir properties.Gas hydrate saturations were estimated from various well logs such as nuclear magnetic resonance (NMR), P- and S-wave velocity, and electrical resistivity logs along with pore-water salinity. Gas hydrate saturations from the NMR log agree well with those estimated from P- and S-wave velocity data. Because of the low salinity of the connate water and the low formation temperature, the resistivity of connate water is comparable to that of shale. Therefore, the effect of clay should be accounted for to accurately estimate gas hydrate saturations from the resistivity data. Two highly gas hydrate-saturated intervals are identified - an upper ∼43 ft zone with an average gas hydrate saturation of 54% and a lower ∼53 ft zone with an average gas hydrate saturation of 50%; both zones reach a maximum of about 75% saturation.     

Mount Elbert Gas Hydrate Stratigraphic Test Well, Alaska North Slope: Overview of scientific and technical program

The Mount Elbert Gas Hydrate Stratigraphic Test Well was drilled within the Alaska North Slope (ANS) Milne Point Unit (MPU) from February 3 to 19, 2007. The well was conducted as part of a Cooperative Research Agreement (CRA) project co-sponsored since 2001 by BP Exploration (Alaska), Inc. (BPXA) and the U.S. Department of Energy (DOE) in collaboration with the U.S. Geological Survey (USGS) to help determine whether ANS gas hydrate can become a technically and commercially viable gas resource. Early in the effort, regional reservoir characterization and reservoir simulation modeling studies indicated that up to 0.34 trillion cubic meters (tcm; 12 trillion cubic feet, tcf) gas may be technically recoverable from 0.92 tcm (33 tcf) gas-in-place within the Eileen gas hydrate accumulation near industry infrastructure within ANS MPU, Prudhoe Bay Unit (PBU), and Kuparuk River Unit (KRU) areas. To further constrain these estimates and to enable the selection of a test site for further data acquisition, the USGS reprocessed and interpreted MPU 3D seismic data provided by BPXA to delineate 14 prospects containing significant highly-saturated gas hydrate-bearing sand reservoirs. The “Mount Elbert” site was selected to drill a stratigraphic test well to acquire a full suite of wireline log, core, and formation pressure test data. Drilling results and data interpretation confirmed pre-drill predictions and thus increased confidence in both the prospect interpretation methods and in the wider ANS gas hydrate resource estimates. The interpreted data from the Mount Elbert well provide insight into and reduce uncertainty of key gas hydrate-bearing reservoir properties, enable further refinement and validation of the numerical simulation of the production potential of both MPU and broader ANS gas hydrate resources, and help determine viability of potential field sites for future extended term production testing. Drilling and data acquisition operations demonstrated that gas hydrate scientific research programs can be safely, effectively, and efficiently conducted within ANS infrastructure. The program success resulted in a technical team recommendation to project management to drill and complete a long-term production test within the area of existing ANS infrastructure. If approved by stakeholders, this long-term test would build on prior arctic research efforts to better constrain the potential gas rates and volumes that could be produced from gas hydrate-bearing sand reservoirs.     

Hydrate bearing clayey sediments: Formation and gas production concepts

Hydro-thermo-chemo and mechanically coupled processes determine hydrate morphology and control gas production from hydrate-bearing sediments. Force balance, together with mass and energy conservation analyses anchored in published data provide robust asymptotic solutions that reflect governing processes in hydrate systems. Results demonstrate that hydrate segregation in clayey sediments results in a two-material system whereby hydrate lenses are surrounded by hydrate-free water-saturated clay. Hydrate saturation can reach ≈2% by concentrating the excess dissolved gas in the pore water and ≈20% from metabolizable carbon. Higher hydrate saturations are often found in natural sediments and imply methane transport by advection or diffusion processes. Hydrate dissociation is a strongly endothermic event; the available latent heat in a reservoir can sustain significant hydrate dissociation without triggering ice formation during depressurization. The volume of hydrate expands 2-to-4 times upon dissociation or CO₂CH₄ replacement. Volume expansion can be controlled to maintain lenses open and to create new open mode discontinuities that favor gas recovery. Pore size is the most critical sediment parameter for hydrate formation and gas recovery and is controlled by the smallest grains in a sediment. Therefore any characterization must carefully consider the amount of fines and their associated mineralogy.     

The mechanism and verification analysis of permafrost-associated gas hydrate formation in the Qilian Mountain,Northwest China

Muri Basin in the Qilian Mountain is the only permafrost area in China where gas hydrate samples have been obtained through scientific drilling. Fracture-filling hydrate is the main type of gas hydrate found in the Qilian Mountain permafrost. Most of gas hydrate samples had been found in a thin-layer-like, flake and block group in a fracture of Jurassic mudstone and oil shale, although some pore-filling hydrate was found in porous sandstone. The mechanism for gas hydrate formation in the Qilian Mountain permafrost is as follows: gas generation from source rock was controlled by tectonic subsidence and uplift--gas migration and accumulation was controlled by fault and tight formation--gas hydrate formation and accumulation was controlled by permafrost. Some control factors for gas hydrate formation in the Qilian Mountain permafrost were analyzed and validated through numerical analysis and laboratory experiments. CSMGem was used to estimate the gas hydrate stability zone in the Qilian permafrost at a depth of 100–400 m. This method was used to analyze the gas composition of gas hydrate to determine the gas composition before gas hydrate formation. When the overlying formation of gas accumulation zone had a permeability of 0.05 × 10⁻¹⁵ m² and water saturation of more than 0.8, gas from deep source rocks was sealed up to form the gas accumulation zone. Fracture-filling hydrate was formed in the overlap area of gas hydrate stability zone and gas accumulation zone. The experimental results showed that the lithology of reservoir played a key role in controlling the occurrence and distribution of gas hydrate in the Qilian Mountain permafrost.