Analytical Modeling of Flow Behavior for Wormholes in Naturally Fractured–Vuggy Porous Media

Acidizing technology has been widely applied when developing naturally fractured–vuggy reservoirs. So testing and evaluating acidizing wells’ pressure behavior become necessary for further improving the wells’ performance. Analyzing transient pressure data can estimate some key reservoir parameters. Generally speaking, carbonate minerals are usually composed of dolomite and calcite which are easy to be dissolved by hydrochloric acid which is often used to react with the rock to create a high conductivity channel, namely wormhole. Pressure transient behavior in fractured–vuggy reservoirs has been studied for many years; however, the models of acidizing wells with wormholes were not reported in previous studies. This article presented an analytical model for wormholes in naturally fractured–vuggy carbonate reservoirs, and wormholes solutions were obtained through point sink integral method. The results were validated accurately by comparing with previous results and numerical simulation. Then in this paper, type curves were established to recognize the flow characteristics, and flow was divided into six flow regimes comprehensively. The calculative results showed that the characteristics of type curves were influenced by inter-porosity flow factor, wormhole number, fluids capacitance coefficient. We also showed that the pressure behavior was affected by the angles between wormholes, and the pressure depletion increased as the angle decreased, because the wormholes were closer, their interaction became stronger. At the end, a reservoir example was showed to demonstrate the methodology of new type curve analysis.     

THE TRANSIENT ELLIPTIC FLOW OF POWER-LAW FLUID IN FRACTAL POROUS MEDIA

1　ＴｈｅＦｌｏｗＭｏｄｅｌｏｆＰｏｗｅｒ＿ＬａｗＦｌｕｉｄｉｎＲａｄｉｃａｌＦｒａｃｔａｌＲｅｓｅｒｖｏｉｒＴｈｅｔｒａｎｓｉｅｎｔｆｌｏｗｏｆｐｏｗｅｒ＿ｌａｗｆｌｕｉｄｉｎｒａｄｉｃａｌｆｒａｃｔａｌｒｅｓｅｒｖｏｉｒｉｓｓｔｕｄｉｅｄｉｎＲｅｆ.[1 ] ,ａｎｄａｎａｌｙｔｉｃａｌｓｏｌｕｔｉｏｎｏｆＬａｐｌａｃｅｓｐａｃｅｉｓｄｅｒｉｖｅｄ .ＩｎＲｅｆ.[2 ] ,ｔｈｅｔｒａｎｓｉｅｎｔｅｌｌｉｐｔｉｃａｌｆｌｏｗｉｓｒｅｓｅａｒｃｈｅｄｏｎｍｏｄｅｌｏｆｅｘｐａｎｄｉｎｇｒｅｃｔａｎｇｌｅ .Ｔ…     

One-Dimensional Matrix-Fracture Transfer in Dual Porosity Systems with Variable Block Size Distribution

Most of the developed models for fractured reservoirs assume ideal matrix block size distribution. This assumption may not be valid in reality for naturally fractured reservoirs and possibly lead to errors in prediction of production from the naturally fractured reservoirs especially during a transient period or early time production from the matrix blocks. In this study, we investigate the effect of variable block size distribution on one- dimensional flow of compressible fluids in fractured reservoirs. The effect of different matrix block size distributions on the single phase matrix-fracture transfer is studied using a recently developed semi-analytical approach. The proposed model is able to simulate fluid exchange between matrix and fracture for continuous or discrete block size distributions using probability density functions or structural information of a fractured formation. The presented semi-analytical model demonstrates a good accuracy compared to the numerical results. There have been recent attempts to consider the effect of variable block size distribution in naturally fractured reservoir modeling for slightly compressible fluids with a constant viscosity and compressibility. The main objective of this study is to consider the effect of variable block size distribution on a one-dimensional matrix-fracture transfer function for single-phase flow of a compressible fluid in fractured porous media. In the proposed semi-analytical model, the pressure variability of viscosity and isothermal compressibility is considered by solving the nonlinear partial differential equation of compressible fluid flow in the fractured media. The closed form solution provided can be applied to flow of compressible fluids with variable matrix block size distribution in naturally fractured gas reservoirs.     

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

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

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

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

致密砂岩垂直裂缝井CO2混相驱试井模型研究

向地层注入CO2可以有效地提高致密砂岩原油采收率,常规的试井解释数学模型不能满足致密砂岩无限导流垂直裂缝井CO2混相驱试井解释的需求．因此,基于渗流力学基本原理建立考虑应力敏感影响的无限导流垂直裂缝井CO2混相驱试井解释数学模型,利用Laplace变换、摄动变换和Stehfest数值反演的方法进行求解,编程绘制典型特征曲线并进行敏感性分析．研究表明：该模型典型特征曲线共划分为八个流动阶段．由于应力敏感效应的影响,径向流阶段内、外区压力导数曲线不再是0.5和0.5M13水平线,而是呈“上翘”的曲线,并且应力敏感系数越大,曲线“上翘”越明显;混相区压力导数曲线符合幂律型变化指数规律且高于(1-n)/(3-n)斜率直线;内区、混相区半径和M12的变化都会使得外区压力曲线升高;通过该模型可以有效地对致密砂岩压裂井CO2混相驱试井资料进行解释．     

Well Test Solutions for Vertically Fractured Injection Wells

The transient behavior of a vertically fractured pressure response due to the presence of an infinite-conductivity vertical fracture is determined by solving the diffusivity equation in elliptical coordinates. The solution is then extended to a composite elliptical system to provide for the different fluid banks present during water injection. The validity of the analytical solutions presented is demonstrated by comparing limiting forms with those available elsewhere in the literature. Computational issues which became evident during the verification stage of our work are also discussed. The solutions have been developed in the Laplace domain to facilitate the addition of fissures and variable rate production (i.e. wellbore storage).A pressure transient test for tracking the advancement of a water front during the early stages of waterflooding is described. We utilize the composite elliptical model developed herein to provide for two distinct regions in which the flow behavior resulting from an induced fracture is elliptical rather than radial. A relationship between the increasing elliptical distance to the waterfront and the resulting change in the apparent (total) skin factor is obtained. Through the analysis of successive falloff tests, this relationship may be used to monitor the advancement of the front provided the cumulative volume of injected water is known. The fluid saturations and the mobilities of the swept and unswept regions are assumed unknown and are obtained from the test analysis.Finally, we present methods for computing the Mathieu functions necessary in solving the diffusivity equation in elliptical coordinates. Mathieu functions are utilized in many applications involving elliptical geometry and we feel the efficient evaluation of these functions is an important contribution of this work.     

Time-Dependent Shape Factors for Uniform and Non-Uniform Pressure Boundary Conditions

Matrix–fracture transfer functions are the backbone of any dual-porosity or dual-permeability formulation. The chief feature within them is the accurate definition of shape factors. To date, there is no completely accepted formulation of a matrix–fracture transfer function. Many formulations of shape factors for instantly-filled fractures with uniform pressure distribution have been presented and used; however, they differ by up to five times in magnitude. Based on a recently presented transfer function, time-dependent shape factors for water imbibing from fracture to matrix under pressure driven flow are proposed. Also new matrix–fracture transfer pressure-based shape factors for instantly-filled fractures with non-uniform pressure distribution are presented in this article. These are the boundary conditions for a case for porous media with clusters of parallel and disconnected fractures, for instance. These new pressure-based shape factors were obtained by solving the pressure diffusivity equation for a single phase using non-uniform boundary conditions. This leads to time-dependent shape factors because of the transient part of the solution for pressure. However, approximating the solution with an exponential function, one obtains constant shape factors that can be easily implemented in current dual-porosity reservoir simulators. The approximate shape factors provide good results for systems where the transient behavior of pressure is short (a case commonly encountered in fractured reservoirs).     

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

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

Transient Pressure Behavior of Reservoirs with Discrete Conductive Faults and Fractures

Fractures and faults are common features of many well-known reservoirs. They create traps, serve as conduits to oil and gas migration, and can behave as barriers or baffles to fluid flow. Naturally fractured reservoirs consist of fractures in igneous, metamorphic, sedimentary rocks (matrix), and formations. In most sedimentary formations both fractures and matrix contribute to flow and storage, but in igneous and metamorphic rocks only fractures contribute to flow and storage, and the matrix has almost zero permeability and porosity. In this study, we present a mesh-free semianalytical solution for pressure transient behavior in a 2D infinite reservoir containing a network of discrete and/or connected finite- and infinite-conductivity fractures. The proposed solution methodology is based on an analytical-element method and thus can be easily extended to incorporate other reservoir features such as sealing or leaky faults, domains with altered petrophysical properties (for example, fluid permeability or reservoir porosity), and complicated reservoir boundaries. It is shown that the pressure behavior of discretely fractured reservoirs is considerably different from the well-known Warren and Root dual-porosity reservoir model behavior. The pressure behavior of discretely fractured reservoirs shows many different flow regimes depending on fracture distribution, its intensity and conductivity. In some cases, they also exhibit a dual-porosity reservoir model behavior.     

Non-linear modeling and dynamic analysis of hydraulic control valve; effect of a decision factor between experiment and numerical simulation

Numerical problems that are usually ignored in the dynamic analysis of hydraulic control valves are described, and an analysis of the effects of such problems on the numerical modeling is provided. Previous studies have ignored the effects of changes in the flow coefficient in the orifice, the solenoid force along the spool movement in the valve and an ascending tendency of pressure during reach to the steady state. To eliminate these problems, it was studied a method to substantiate the non-linearity of the pressure loss caused by passing between the orifice and port as well as that caused by interaction with the solenoid. Moreover, the movement of the spool and spring is expressed using the time-delay-element (TDE). The proposed numerical model has been used in the Bond graphs method of a hydraulic control valve and the simulation results have been shown to be accurate. It is known that differences between simulated and experimental results can have a considerable impact on the function of actual systems. The contribution of each parameter is measured separately for the transient state and steady state. Analysis standard observed the first peak value, pressure increase to the steady state and the settling time in the response results.     

分形在油气田开发中的应用

分形用于测井资料分析，为油藏描述提供了新的工具，使得利用油藏数值模拟来模拟复杂结构的油气藏变得更加便利；分形油藏上的试井分析解决了一些以前难以解释的问题；用分形来描述裂缝和孔洞结构，提高了酸化压裂的设计水平；用分形理论来描述两相流，揭示了部分粘性指进的触发机理．总之，分形理论在油气田开发中的应用，促进了渗流力学基础理论的发展．     

A new computation model of fractured horizontal well coupling with reservoir

Based on Green's functions and Newman's product principle, pressure drop formula was derived for considering simultaneous production of fractures and horizontal wellbore in unsteady state. A reservoir/fractured horizontal well coupling model is developed for finite conductivity condition that can be solved by the combination of quasi‐Newton method and PSO (Particle Swarm Optimization) algorithm. The solution of a practical example shows that lots of factors can affect the productivity of the fractured horizontal well. The number of the fractures has an optimizing range. Different fractures have different flow rates and the production in the central fracture is the lowest. The distribution of production rate in the wellbore has a wave‐like shape due to the influence of fractures. The closer fractures are the lower flow rate in well segments. Considering the horizontal wellbore production, the production rates in the fractures are asymmetric and the corresponding pressure drop curve is smoother. Copyright © 2010 John Wiley & Sons, Ltd.     

Unknown Rate History Calculation from Down-hole Transient Pressure Data Using Wavelet Transform

In the oilfield, rate metering is less satisfying compared with pressure measurement. Detailed rate history for individual well is not always available. The common practice is that few of the wells are measured together through manifolds and daily, monthly rate or total cumulative production is measured with surface separation equipments. Rate allocation algorithm is utilized to assign production to each well, but the results are inaccurate due to low measurement frequency. For well testing, another problem is that the rate measured from the surface is very often inconsistent with the pressure measured from down-hole. These problems together make it a very challenging task for engineer to interpret transient pressure from the well down-hole gauges. This article introduced wavelet transform method to calculate unknown rate history from the down-hole transient pressure data. The Haar wavelet was chosen to process transient pressure data, where flow events can be identified with the amplitude of pressure transform. For the reservoir with constant reservoir-well parameters, the amplitude of pressure transforms is proportional to the change of rate. Based on this discovery, an algorithm for calculating unknown rate history from cumulative production has been developed. As this method is independent of the developed reservoir model, reservoir heterogeneity, skin, and wellbore storage have no effect on this method. For other complex reservoirs with variable reservoir-well parameters, such as gas reservoir and multiphase reservoir, this algorithm can also be applied resulting in small errors. In practice, the developed method can identify the offset of the flowing rate and the pressure, so synchronization of the rate and pressure can be achieved, which greatly increased the confidence of the pressure data analysis. Both theory and case studies were tested to demonstrate the robustness of this method in practice.     

Displacement of Cr(III)–Partially Hydrolyzed Polyacrylamide Gelling Solution in a Fracture in Porous Media

During waterflooding of a fractured formation, water may channel through the fracture or interconnected network of fractures, leaving a large portion of oil bearing rock unswept. One remedial practice is injection of a gelling solution into the fracture. Such placement of a gelling mixture is associated with leak-off from the fracture face into the adjoining matrix. Design of a gel treatment needs understanding of the flow of gelling mixture in and around the fracture. This flow is addressed here for Cr(III)–partially hydrolyzed polyacrylamide formulation through experiments and conceptual model. A fractured slab was used to develop a lab-model, where the flow along the fracture and simultaneous leak-off into the matrix can be controlled. Also, the fracture and matrix properties had to be evaluated individually for a meaningful analysis of the displacement of gelling solution. During this displacement, the gelling fluid leaked off from the fracture into the matrix as a front, resulting in a decreasing velocity (and pressure gradient) along the fracture. With pressure in the fracture held constant with time, the leak-off rate decreased as the viscous front progressed into the matrix. The drop in leak-off rate was rapid during the initial phase of displacement. A simple model, based on the injection of a viscous solution into the dual continua, could explain the displacement of Cr(III)–polyacrylamide gelling mixture through the fractured slab. This study rules out any major complication from the immature gelling fluid, e.g., build-up of cake layer on the fracture face. The model, due to its simplicity may become useful for quick sizing of gel treatment, and any regression-based evaluation of fluid properties in a fracture for other applications.     

Pressure recovery processes in fractured reservoirs

The analysis of the pressure recovery in an oil field following the shutting in of a well is one of the basic methods of determining or refining the reservoir parameters. For fractured reservoirs the process is complicated by the strong dependence of their flow-capacity characteristics on the state of stress. It is shown that, qualitatively, the pressure recovery processes in fractured and fractured porous media can be correctly described using the integral conservation laws. Numerical solutions are derived in order to estimate the validity of the results obtained.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 5, pp. 77–83, September–October, 1990.     

Pressure Buildup During Supercritical Carbon Dioxide Injection From a Partially Penetrating Borehole into Gas Reservoirs

Injecting CO₂ into a subsurface formation causes a buildup of pressure in the vicinity of the injection well. While a large injection rate can reduce the cost associated with injection, an indefinitely large injection rate can result in excessive formation damage. To obtain an optimal injection rate without exceeding the safe pressure limits, one will like to have some knowledge of the transient pressure buildup characteristics resulting from a particular injection rate. While elaborate numerical simulations can provide reliable pressure buildup predictions, they require extensive knowledge about the formation, which is normally not available at the start of an injection process. To alleviate this problem, using some simplifying assumptions, we have developed a solution to predict the transient buildup of pressure resulting from injection of supercritical carbon dioxide from a partially penetrating well into a gas reservoir. The solution in space and time is first obtained in the Fourier–Laplace transform space, and then inverted back into real space (in cylindrical coordinates) and time. We use the solution to study pressure transient characteristics for different formation permeabilities and anisotropy ratios. Results obtained using the solution compared well with those from numerical simulations.     

Simulation of Pressure Transient Behavior for Asymmetrically Finite-Conductivity Fractured Wells in Coal Reservoirs

Based on Fick’s law in matrix and Darcy flow in cleats and hydraulic fractures, a new semi-analytical model considering the effects of boundary conditions was presented to investigate pressure transient behavior for asymmetrically fractured wells in coal reservoirs. The new model is more accurate than previous model proposed by Anbarci and Ertekin, SPE annual technical conference and exhibition, New Orleans, 27–30 Sept 1998 because new model is expressed in the form of integral expressions and is validated well through numerical simulation. (1) In this paper, the effects of parameters including fracture conductivity, coal reservoir porosity and permeability, fracture asymmetry factor, sorption time constant, fracture half-length, and coalbed methane (CBM) viscosity on bottomhole pressure behavior were discussed in detail. (2) Type curves were established to analyze both transient pressure behavior and flow characteristics in CBM reservoir. According to the characteristics of dimensionless pseudo pressure derivative curves, the process of the flow for fractured CBM wells was divided into six sub-stages. (3) This paper showed the comparison of transient steady state and pseudo steady state models. (4) The effects of parameters including transfer coefficient, wellbore storage coefficient, storage coefficient of cleat, fracture conductivity, fracture asymmetry factor, and rate coefficient on the shape of type curves were also discussed in detail, indicating that it is necessary to keep a bigger fracture conductivity and fracture symmetry for enhancing well production and reducing pressure depletion during the hydraulic fracturing design.     

Transient pressure of percolation through one dimension porous media with threshold pressure gradient

IntroductionManyexperimentsshow[1,2,3]thatthepercolationinlowPermeabilityisnotfittoDarcy'sLaw.Theparticularcharacteristheexistenceofthresholdpressuregradient[1],inotherwords:thefluidcanflow,onlywhentherealpressuregradientislargerthanthresholdpressuregradient.Thedevelopmentinthisresearchingfieldhasn'tmuchprogressbecauseofquestionsineconomicsandteelmologysincethisconceptwasformedin1951.ThisfieldisbecominganewwarmresearchingPOintasChinadevelopedthelowPermeabilityreservoirsinlastfewyears.Thego…     

Prediction of transient displacement, velocity, and force on projectiles penetrating cohesive soils

A closed-form solution is developed for predicting transient displacement, velocity, and force as well as residual displacement and duration for both flat and pointed projectiles penetrating cohesive soils. Because transient response can be represented as a 3 parameter space of nondimensional numbers, displacement, velocity, and force are each presented in a single graph. Experimental test data for different soils and penetrator contours are used to demonstrate the validity of both transient and residual responses. Of great importance in this solution and soil dynamics in general is the advancement of a hypothesis for the dynamic interaction of intergranular stress, pore water pressure, and pore air pressure in the resisting force that retards penetration. Because fundamental dynamic pore pressure theory is used to develop this analysis, basic soil properties (cohesive yield strength, density, and degree of saturation) are used to characterize the cohesive soil.