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

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

ANALYSIS OF THE WELLBORE STABILITY BASED ON PLASTIC DAMAGE THEORY IN SOFT SHALE FORMATION

软泥页岩在地层高压、流体环境下独有的力学变形特性一直是井壁稳定性控制工程领域的研究热点。针对软泥页岩低强度、高塑性的特点,考虑孔隙渗流、塑性损伤引起的渗透率变化等因素,建立了渗流-- 应力耦合形式的井壁稳定性评价数值模型,研究了流体渗流、泥浆压力、水平应力差的影响。结果表明,对过平衡钻井,当泥浆压力超过地层孔隙压力10% 时,渗流-- 应力耦合效应导致井壁塑性损伤区扩大至少65%,需予以考虑。水平应力差的增大会促进井壁塑性损伤区向最小水平主应力方向集中,且损伤区深度增大。     

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

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

Study on steady surge pressure for Yield-Pseudoplastic fluid in a concentric annulus

Surge pressure caused by pipe moving through the wellbore full of drilling fluid is the main factor affecting wellbore stability, on which the additional mud density is designed. Thus it is desired to predict the surge pressure with high accuracy. While this is closely related with the selection of mud rheological models. Yield-Pseudoplastic model which has three parameters and is famous for its higher accuracy has been used to describe the mud rheological properties in recent years. Based on this model the paper presents a new theoretical model for calculating surge pressure caused by mud viscosity during tripping or casing in a concentric annulus of directional wells under steady laminar situation. For convenience, the paper plots the distribution of the coefficients of surge pressure for different conditions. An example is given. These results provide the basis for controlling pressure surges and tripping velocity.     

STUDY ON STEADY SURGE PRESSURE FOR YIELD-PSEUDOPLASTIC FLUID IN A CONCENTRIC

ＷａｎｇＨａｉｇｅ（汪海阁）；ＬｉｕＸｉｓｈｅｎｇ（刘希圣）（ＲｅｃｅｉｖｅｄＪｕｌｙ８，１９９４；ＣｏｍｍｕｎｉｃａｔｅｄｂｙＣｈｉｅｎＷｅｉｚａｎｇ）Ａｂｓｔｒａｃｔ：Ｓｕｒｇｅｐｒｅｓｓｕｒｅｃａｕｓｅｄｂｙｐｉｐｅｍｏｖｉｎｇｔｈｒｏｕｇｈｔ...     

井筒气液两相流对致密气压裂水平井试井的影响

多段压裂水平井技术是目前开采致密气最常用的方法之一,在致密气压裂水平井试井测试中常常伴随着一定的产水量,井筒气液两相流会增加井筒流体的流动阻力,加大井筒流体流动对试井解释的影响.为了明确井筒气液两相流对致密气藏压裂水平井试井的影响,提高产水致密气压裂水平井的试井解释精度,建立了一种井筒气液两相流与地层渗流耦合的试井模型,采用数值方法对模型进行求解,获得了考虑井筒气液两相流的压裂水平井试井理论曲线、压力场分布及裂缝产量分布.研究结果表明:井筒气液两相流会增加试井理论曲线中压力和压力导数值,造成靠近入窗点的压力扩散要快于远离入窗点的压力扩散,引起靠近入窗点的裂缝产量要高于远离入窗点的裂缝产量.现场实例分析进一步说明,不考虑井筒两相流可能会对产水压裂水平井的试井解释结果产生很大误差,主要表现为水平井筒假设为无限大导流能力会使得拟合得到的表皮系数偏大,将测试点视为入窗点会使得拟合得到的原始地层压力偏小.所建立的考虑井筒两相流的压裂水平井试井模型为产水致密气井试井资料的正确解释提供了重要技术保障.      

Meso-level simulation of gas hydrate dissociation inlow-permeability sediments

This paper presents a meso-level simulation of gas hydrate dissociation in low-permeability marine sediments. Interstitial pores are defined to describe fluid flow and particle movement. The proposed model couples multiphase fluid flow with particle movement to simulate the thermodynamics of gas hydrate dissociation triggered by sharp temperature rises. Hydrates respond quickly to temperature rise in low-permeability sediments. Dissociation causes pore pressure to rise rapidly to equilibrium then steadily increase above equilibrium pressure. Lower permeability sediment builds up greater excess pore pressure as the dissipation of pore pressure is constrained.     

Self-Similarity in Contamination Transport to a Formation Fluid Tester During Cleanup Production

Sampling hydrocarbon fluids from a reservoir during an early phase of formation testing (typically during interruptions in drilling) represent an important step in formation evaluation targeting the characterization of composition and pressure/volume/temperature (PVT) properties of petroleum reserves required for their efficient development. A sampling job is usually accomplished with a tool, the formation tester, which can be conveyed downhole on a wireline cable, drill-pipe, or tubing. This tool has a probe and a pump allowing for the production of reservoir fluid from a small spot of a borehole wall covered by filtercake or mudcake deposited during drilling. The filtercake is built by the invasion of mud filtrate into the formation under pressure overbalance created by the mud column inside the wellbore. For this reason, after drilling, each wellbore is surrounded by a cylindrical zone saturated with mud filtrate and, prior to obtaining a sample of virgin formation fluid, cleanup production must be initiated to overcome the consequences of invasion. A 3D model of flow and contamination transport during cleanup production is considered in this article. The model assumes a piston-like displacement and no filtrate leakage through mudcake. This model leads to a multiscale problem of flow and transport in porous media with geometrical and hydrodynamic singularities. A full sensitivity analysis of its solutions has been carried out versus the initial depth of filtrate invasion and the viscosity contrast between the formation fluid and invaded mud filtrate. The numerical modeling has revealed two phases of cleanup production. During the early phase, the contamination of produced fluid is predominantly affected by a circumferential contamination transport. During the later phase, the evolution of contamination is controlled mainly by mud-filtrate displacement vertically. In the absence of viscosity contrast, both phases of cleanup can be described by exponential laws for the contamination versus the produced volume with exponents matching the empirical correlation (−5/12) and the far-field pseudospherical flow pattern (−2/3) for the early and late phases of cleanup, respectively. In the presence of viscosity contrast, the far-field asymptote holds whereas the circumferential regime of cleanup is strongly affected by the viscosity ratio. The application of solutions found to the cleanup production monitoring is briefly described.     

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

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

Multiphase flow model developed for simulating gas hydrate transport in horizontal pipe

The hydrate formation or dissociation in deep subsea flow lines is a challenging problem in oil and gas transport systems. The study of multiphase flows is complex while necessary due to the phase changes (i.e., liquid, solid, and gas) that occur with increasing the temperature and decreasing the pressure. A one-dimensional multiphase flow model coupled with a transient hydrate kinetic model is developed to study the characteristics of the multiphase flows for the hydrates formed by the phase changes in the pipes. The multiphase flow model is derived from a multi-fluid model, while has been widely used in modelling multiphase flows. The heat convection between the fluid and the ambient through the pipe wall is considered in the energy balance equation. The developed multiphase flow model is used to simulate the procedure of the hydrate transport. The results show that the formation of the hydrates can cause hold-up oscillations of water and gas.     

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.     

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.     

波动钻压下阶梯钻柱系统稳定性研究

处于狭长井筒中的钻柱,其动力响应受到钻具组合、内外钻井液流动以及钻井参数等因素的影响,钻柱动力失稳导致的剧烈振动是井壁坍塌和钻具失效的重要原因。考虑到钻杆和钻铤在刚度和线密度上存在很大的差别,论文将钻柱简化为单阶梯输液管柱,钻井液沿着钻柱内部向下泵入并从环空返回地面。耦合考虑钻柱自重、随时间简谐变化的波动钻压、稳定器以及钻井液的水动力和阻尼力,建立了直井中钻柱横向振动的解析模型。利用有限单元法离散为四阶常微分方程后,采用Bolotin法得到临界频率方程确定系统的不稳定区范围,研究了钻压、钻杆长度、稳定器安装位置、钻井液的流速和密度等参数对系统稳定性影响的机理。研究表明：钻压的平均值和波动幅值都是钻柱失稳的驱动因素,而系统的稳定性对处于受拉状态的钻杆的长度变化不敏感。在论文所研究的参数范围内,降低钻井液流速和密度、下移稳定器的安装位置均有助于增强系统的稳定性。     

虚拟质量力对酸性气体-钻井液两相流波速的影响

基于两流体模型、酸性气体和钻井液状态方程,考虑酸性气体与钻井液相间虚拟质量力、粘性剪切力、相间动量交换及狭义相间阻力等条件,建立酸性气体与钻井液两相中压力波传播速度的数学模型,依据小扰动原理,对波速模型求解,得到关于波数K的波速方程。结果表明,在一定范围内,随空隙率、频率的增大,虚拟质量力对波速的影响显著增强;在高空隙率下,压强增大,虚拟质量力对波速的影响减弱;增大流体的密度或不可压缩性,均可使两相压力波速增大;延长气液交换时间或减小波动频率使相间有足够时间进行动量交换,两相压力波波速随之减小。     

One- and two-phase flow during fluid sampling by a wireline tool

Wireline sampling tools withdraw a few litres of fluid from a permeable formation via a small sink probe pressed against the borehole wall. The aim is to recover, quickly and cheaply, a representative native fluid sample. Unfortunately, the formation in the near wellbore region is invaded by mud filtrate, and withdrawal of nonnative fluid initially is inevitable. It is therefore of interest to estimate the proportion of native fluid in the sample stream, as a function of time.Semi-analytical calculations of one- and two-phase sampling flows are presented, for the special case of constant total fluid mobility in the limits of very deep or very shallow invasion. Both the interaction of the initially cylindrically symmetric invasion profile with the spherically symmetric flow and the capillary shock-forming dynamics of two-phase flow are found to control the character of sample composition variation. The wide variety of sample stream composition histories is displayed.     

Front Regime of Heat and Mass Transfer in a Gas Hydrate Reservoir under the Negative Temperature Conditions

The analytical self-similar solution to the nonlinear problem of the front regime of heatand- mass transfer in a gas hydrate reservoir under the negative temperature conditions is obtained. In the initial state the reservoir is assumed to be saturated with a heterogeneous gas hydrate–ice–gas mixture. In particular cases there may be no ice or/and gas. The ice and gas are formed behind the gas hydrate dissociation front. The calculations are presented for a stable hydrate–gas system. The critical curves are constructed in the well-pressure–reservoir-permeability plane. These curves separate the domains of the front regime and the regime of volume gas hydrate dissociation ahead of the front. The velocity of the gas hydrate dissociation front is investigated as a function of various problem parameters. The characteristic temperature and pressure distributions corresponding to various regimes on the diagram are investigated.     

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.     

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

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

Numerical and experimental study on parameters distribution inside a two-phase ejector

The two-phase flow process in an ejector was numerically and experimentally studied using R141b as a working fluid. A modified one-dimensional gas–liquid ejector model was proposed to remedy the defect in the traditional one. Gas–liquid boundary layer regions were discussed and used to close the model. Mac Cormack method is used to discrete controlling equations of gas–liquid two-phase flow in the ejector. The radial distribution of velocity and temperature, the variation of void fraction, the axial velocity variation and the influence of primary steam pressure on the mixing process were predicted with the numerical model. An experimental rig was set up to validate the model by comparing the experimental pressure distribution in the ejector with the calculating one.