页岩气藏压裂水平井试井分析

页岩气藏资源丰富,开发潜力巨大,已成为目前研究的热点.与常规气藏相比,页岩气藏运移机制复杂,流动模式呈非线性,有必要考虑页岩气的吸附解吸,天然微裂缝的应力敏感性,人工裂缝内的非达西流等非线性因素对压裂水平井压力响应的影响. 基于双重介质和离散裂缝混合模型,分别采用Langmuir等温吸附方程描述吸附解吸,渗透率指数模型描述应力敏感,Forchheimer方程描述非达西效应,建立页岩气藏压裂水平井数值试井模型. 运用伽辽金有限元法对模型进行求解.根据试井特征曲线,划分流动阶段,着重分析非线性因素对压力响应的影响.结果表明:页岩气藏压裂水平井存在压裂裂缝线性流、压裂裂缝径向流、地层线性流、系统径向流及封闭边界影响5 种流动阶段.吸附解吸的影响发生窜流之后,Langmuir吸附体积增大,拟压力导数曲线凹槽更加明显,系统径向流出现时间与压力波传播到边界时间均延迟;天然裂缝系统的应力敏感性主要影响试井曲线的晚期段,拟压力和拟压力导数曲线均表现为上翘,应力敏感效应越强,上翘幅度越大;高速非达西效应对早期段影响较大,非达西效应越强,拟压力降幅度越大,试井曲线上翘.与解析解的对比以及矿场实例验证了模型的正确性与适用性.      

考虑诱导缝多段压裂水平井非均质改造压力动态模型研究

为了准确模拟致密油藏水平井大规模压裂形成复杂裂缝网络系统和非均质储层井底压力变化,建立考虑诱导缝矩形非均质储层多段压裂水平井不稳定渗流数学模型,耦合裂缝模型与储层模型得到有限导流裂缝拉普拉斯空间井底压力解,对两种非均质储层模型分别利用数值解、边界元和已有模型验证其准确性.基于压力导数曲线特征进行流动阶段划分和参数敏感性分析,得到以下结果:和常规压裂水平井井底压力导数曲线相比较,理想模式下,考虑诱导缝影响时特有的流动阶段是综合线性流阶段、诱导缝向压裂裂缝“补充”阶段、储层线性流动阶段和拟边界控制流阶段.诱导缝条数的增加加剧了综合线性流阶段的持续时间,降低了流体渗流阻力,早期阶段压力曲线越低;当诱导缝与压裂裂缝导流能力一定时,裂缝导流能力越大,线性流持续时间越长;当所有压裂裂缝不在一个区域时,沿井筒方向两端区域低渗透率弱化了低渗区域诱导缝流体向压裂裂缝“补充”阶段,因此,沿井筒方向两端区域渗透率越低,早期阶段压力曲线越高;当所有压裂裂缝在一个区域时,渗透率变化只影响径向流阶段之后压力曲线形态,外区渗透率越低,早期径向流阶段之后压力曲线越高.通过实例验证,表明该模型和方法的实用性和准确性.     

考虑多重运移机制耦合页岩气藏压裂水平井数值模拟

页岩作为典型的微纳尺度多孔介质,游离气与吸附气共存,传统的达西定律已无法准确描述气体在页岩微纳尺度的运移规律.基于双重介质模型和离散裂缝模型构建页岩气藏分段压裂水平井模型,其中基岩中考虑气体的黏性流、Knudsen 扩散以及气体在基岩孔隙表面的吸附解吸,吸附采用Langmuir等温吸附方程;裂缝中考虑黏性流和Knudsen扩散,在此基础上建立基岩-裂缝双重介质压裂水平井数学模型并采用有限元方法对模型进行求解.结果表明,基岩固有渗透率越小,表面扩散和Knudsen扩散的影响越大,反之则越小;人工裂缝的性质包括条数、开度、半长以及间距,主要影响压裂水平井生产早期,随着人工裂缝参数值的增加,压裂水平井产能增加,累产气量也越大.其次,页岩气藏压裂诱导缝和天然裂缝的发育程度对页岩气藏的产能有很大的影响,水平井周围只有人工裂缝,周围天然裂缝不开启或不发育时,页岩气藏的水平井的产能较低.      

考虑单元劈裂的流-固耦合连续-非连续方法及定向水力压裂模拟

为了更真实地模拟水力压裂过程中的岩石变形、裂缝扩展及流体流动,在自主开发的拉格朗日元与离散元耦合的连续-非连续方法的基础上,发展了一种流-固耦合方法。在该方法中,裂缝可沿四边形单元对角线和单元边界扩展,流体流动满足立方定律。通过与单一裂缝非稳态渗流模型及KGD模型的理论解进行对比,验证了该方法的正确性。由定向射孔水力压裂的模拟结果可以发现,(1)距离射孔越远,流体压力越小;随着时间的增加,裂缝中流体压力降低。(2)随着射孔角度的增加,裂缝起裂和扩展过程中的流体压力及转向距离增加;随着x方向水平应力的增加,裂缝起裂和扩展过程中的流体压力增加;两个方向水平应力之差越大,裂缝转向距离越小。(3)随着时间的增加,裂缝区段数目的增速变慢,这与裂缝体积增加变快有关。     

考虑非达西渗流的单层非均质地基固结分析

各类地质作用形成的非均质粘性土地基在工程界中广泛存在,同时非均质粘性土中渗流在低水力坡降下会出现非达西渗流的现象。但考虑非达西渗流的非均质粘性土地基一维固结理论尚不完善,应用达西定律计算低水力坡降情况下非均质粘性土中的渗流,会存在较大的误差。因此通过假定粘性土压缩模量随深度线性增加且渗透系数随深度指数减小来考虑地基的非均质性,并考虑实际中变荷载作用建立了非达西渗流下非均质地基固结模型。在此基础上,运用有限差分法获其数值解,并与相应的达西渗流解析解进行对比,验证了本文解的正确性;分别研究了非达西渗流参数m、i_1和非均质地基参数α、μ、p对土体孔压、地基沉降、地基平均固结度的影响。结果表明:m和i_1值越大,非均质地基中孔压消散越慢,土体固结速率也随之变慢;非均质粘性土地基参数的影响各异,α越大或μ越大、p越小的情况下,固结速率越快,且随着α值增大,非均质土体的最终沉降值越小;而m、i_1、μ、p对地基最终沉降量没有影响。非达西渗流和地基的非均质性都会影响土体的固结性状,故此两者在固结计算过程中不可忽略。     

复合射孔爆燃气体压裂裂缝起裂扩展研究

结合线弹性断裂力学的裂缝尖端应力强度因子判据,建立了复合射孔爆燃气体压裂裂缝的起裂扩展模型,通过建立与多个变量相关的缝内气体压力分布函数,利用迭代法实现了模型的数值求解,获得了缝内气体压力分布随时间的动态变化规律,并分析了不同特征参数对裂缝起裂扩展与止裂过程的影响。实例计算结果表明:(1)随着裂缝扩展的进行,爆燃气体流动尖端与裂缝尖端经历了由重合到不重合再到重合的过程;(2)地应力越大,裂缝起裂扩展越困难,爆燃气体有效致裂作用时间越短,最终得到的裂缝扩展长度也越小;(3)初始裂缝越长,裂缝更容易起裂扩展,爆燃气体能量利用率越高,裂缝扩展更长;(4)岩石断裂韧性的改变对裂缝起裂、止裂和裂缝扩展长度没有明显的影响;(5)升压速率越小,爆燃气体有效致裂作用时间越长,最终裂缝扩展也更长,但对裂缝起裂压力与止裂压力几乎没有影响。     

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

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

动边界双重介质油藏低速非达西渗流试井模型

裂缝性油藏中基质岩块的渗透率一般很低,大量岩心测试实验证实在基质岩块内的液体渗流和在一定含水饱和度下的气体渗流将偏离达西渗流,往往出现低速非达西渗流,表现出启动压力梯度以及流体流动边界不断向外扩展等特殊现象。本文充分考虑启动压力梯度与动边界的影响,建立了微可压缩双重介质油藏低速非达西渗流的试井数学模型,对时间和空间变量...     

页岩气藏压裂缝网扩展数值模拟

为探究页岩气藏水力压裂复杂裂缝网络的形成机理,开展了缝网扩展的数值模拟研究.考虑应力阴影和天然裂缝作用,建立了井筒和裂缝中流体流动模型,利用位移不连续方法求解应力与位移不连续量,然后构建了压力与裂缝宽度的迭代方程,并采用牛顿迭代法求解.通过比较数值解经典模型解析解,验证了模型和迭代解法的正确性.多簇裂缝同步扩展时裂缝间距越小,压裂液分配到各条裂缝越不均匀,靠近井筒跟部的裂缝的分流量越大,从而裂缝宽度越大;考虑天然裂缝作用时,逼近角越小或者应力各向异性越弱,水力裂缝越容易发生转向扩展,裂缝网络越复杂.      

页岩气藏的双重介质-离散裂缝模型

考虑页岩气藏开发中渗流的多尺度效应,提出了一个基于裂缝-孔隙双重介质的离散裂缝模型．在该模型中,基质、天然裂缝和人工压裂裂缝采用各自控制方程独立计算,不同介质之间通过流量交换相互关联．为分析模型可靠性,分别和基于渗透率粗化及压裂裂缝导流能力无穷大的模型对比．数值算例显示,伴随着网格细分,该模型与精确渗透率粗化模型具有相同计算精度,两者收敛速度均较快,但该模型易推广到多相流动问题,而等压模型对产量将有所高估．研究了地质参数和工艺参数对气井产量的影响规律．计算结果表明天然裂缝渗透率及基质孔隙扩散系数对产气速率有着重要影响,产气速率伴随着人工压裂裂缝导流能力、长度以及数目的增加而增加,但是增加幅度会逐步趋缓．     

Numerical Simulation of Single-Phase and Multiphase Non-Darcy Flow in Porous and Fractured Reservoirs

A numerical method as well as a theoretical study of non-Darcy fluid flow through porous and fractured reservoirs is described. The non-Darcy behavior is handled in a three-dimensional, multiphase flow reservoir simulator, while the model formulation incorporates the Forchheimer equation for describing single-phase or multiphase non-Darcy flow and displacement. The non-Darcy flow through a fractured reservoir is handled using a general dual-continuum approach. The numerical scheme has been verified by comparing its results against those of analytical methods. Numerical solutions are used to obtain some insight into the physics of non-Darcy flow and displacement in reservoirs. In addition, several type curves are provided for well-test analyses of non-Darcy flow to demonstrate a methodology for modeling this type of flow in porous and fractured rocks, including flow in petroleum and geothermal reservoirs.     

A New Non-Darcy Flow Model for Low-Velocity Multiphase Flow in Tight Reservoirs

The pore and pore-throat sizes of shale and tight rock formations are on the order of tens of nanometers. The fluid flow in such small pores is significantly affected by walls of pores and pore-throats. This boundary layer effect on fluid flow in tight rocks has been investigated through laboratory work on capillary tubes. It is observed that low permeability is associated with large boundary layer effect on fluid flow. The experimental results from a single capillary tube are extended to a bundle of tubes and finally to porous media of tight formations. A physics-based, non-Darcy low-velocity flow equation is derived to account for the boundary layer effect of tight reservoirs by adding a non-Darcy coefficient term. This non-Darcy equation describes the fluid flow more accurately for tight oil reservoir with low production rate and low pressure gradient. Both analytical and numerical solutions are obtained for the new non-Darcy flow model. First, a Buckley–Leverett-type analytical solution is derived with this non-Darcy flow equation. Then, a numerical model has been developed for implementing this non-Darcy flow model for accurate simulation of multidimensional porous and fractured tight oil reservoirs. Finally, the numerical studies on an actual field example in China demonstrate the non-negligible effect of boundary layer on fluid flow in tight formations.     

Two-phase flow for a horizontal well penetrating a naturally fractured reservoir with edge water injection

This paper examines the two-phase flow for a horizontal well penetrating a naturally fractured reservoir with edge water injection by means of a fixed streamline model. The mathematical model of the vertical two-dimensional flow or oil-water for a horizontal well in a medium with double-porosity is established, and whose accurate solutions are obtained by using the characteristic method. The saturation distributions in the fractured system and the matrix system as well as the formula of the time of water free production are presented. All these results provide a theoretical basis and a computing method for oil displacement by edge water from naturally fractured reservoirs.     

TWO－PHASE FLOW FOR A HORIZONTAL WELL PENETRATING A NATURALLY FRACTURED RESERVOIR WITH EDGE WATER INJECTION

ＴＷＯ－ＰＨＡＳＥＦＬＯＷＦＯＲＡＨＯＲＩＺＯＮＴＡＬＷＥＬＬＰＥＮＥＴＲＡＴＩＮＧＡＮＡＴＵＲＡＬＬＹＦＲＡＣＴＵＲＥＤＲＥＳＥＲＶＯＩＲＷＩＴＨＥＤＧＥＷＡＴＥＲＩＮＪＥＣＴＩＯＮＧｕｏＤａｌｉ（郭大立）ＬｉｕＣｉｑｕｎ（刘慈群）（Ｒｅｃｅｉｖｃ...     

TWO－PHASE FLOW FOR A HORIZONTAL WELL PENETRATING A NATURALLY FRACTURED RESERVOIR WITH EDGE WATER INJECTION

ＴＷＯ－ＰＨＡＳＥＦＬＯＷＦＯＲＡＨＯＲＩＺＯＮＴＡＬＷＥＬＬＰＥＮＥＴＲＡＴＩＮＧＡＮＡＴＵＲＡＬＬＹＦＲＡＣＴＵＲＥＤＲＥＳＥＲＶＯＩＲＷＩＴＨＥＤＧＥＷＡＴＥＲＩＮＪＥＣＴＩＯＮＧｕｏＤａｌｉ（郭大立）ＬｉｕＣｉｑｕｎ（刘慈群）（Ｒｅｃｅｉｖｅ...     

Radiation effect on mixed convection from a horizontal surface in a porous medium

The effect of thermal radiation on the non-Darcy mixed convection flow over a non-isothermal horizontal surface immersed in a saturated porous medium has been studied. The wall temperature is assumed to have a power-law variation with the distance measured from the leading edge of the plate. The non-linear coupled parabolic partial differential equations governing the flow have been solved numerically using a finite-difference scheme. For some particular cases, the self-similar solution has also been obtained. The heat transfer is found to be strongly influenced by the radiative flux number, buoyancy parameter, variation of wall temperature, non-Darcy parameter and the nature of the free stream velocity.     

低渗透多孔介质渗流动边界模型的解析与数值解

考虑启动压力梯度的低渗透多孔介质非达西渗流模型属于强非线性动边界问题, 分别利用相似变量变换方法和基于空间坐标变换的有限差分方法, 对内边界变压力情况下、考虑启动压力梯度的一维低渗透多孔介质非达西渗流动边界模型进行了精确解析与数值求解研究. 研究结果表明:该动边界模型存在唯一的精确解析解, 且所求得的精确解析解可严格验证数值解的正确性;且当启动压力梯度值趋于零时, 非达西渗流动边界模型的精确解析解将退化为达西渗流情况下的精确解析解. 由求解结果作出的非零无因次启动压力梯度下的地层压力分布曲线表现出紧支性特点, 其与达西渗流模型的有显著不同. 因此, 研究低渗透多孔介质中非稳态渗流问题时, 应该考虑动边界的影响. 研究内容完善了低渗透多孔介质的非达西渗流力学理论, 为低渗透油气藏开发的试井解释与油藏数值模拟技术提供了理论基础.