基于核磁共振孔隙划分的致密油藏自发渗吸原油可动性研究

自发渗吸驱油是致密油藏提高采收率的有效手段,但不同的孔隙划分方法会导致原油可动性精细定量表征存在差异性.基于此,以鄂尔多斯盆地延长组致密油藏为研究对象,开展了四种典型致密岩心的自发渗吸驱油实验,利用基于核磁共振分形理论的流体分布孔隙精细划分方法,区分了致密砂岩岩心孔隙类型,明确了不同类型岩心孔隙结构对原油可动性和自发渗吸驱油速率的控制特征.研究结果表明不同类型岩心自发渗吸模拟油动用程度介于22.07%～33.26%,核磁共振T2谱双峰型岩心自发渗吸模拟油动用程度高于单峰型岩心;不同类型致密砂岩岩心中流体分布孔隙可初步划分出P1和P2两种类型, P1类型孔隙则可进一步划分出P1-1, P1-2和P1-3三种亚类型孔隙;致密砂岩岩心中P1和P2类孔隙中模拟油均有不同程度的动用, P1类孔隙作为致密岩心中主要孔隙,尤其是P1类孔隙中P1-2和P1-3类孔隙的数量决定了自发渗吸模拟油动用程度;P1-1, P1-2和P1-3类孔隙结构差异性对自发渗吸模拟油动用程度起决定性作用,较小尺寸孔径孔隙较大的孔隙结构差异性不仅提升了自发渗吸模拟油动用程度,而且提升了自发渗吸驱油速率;流体可动性指数较高的P...     

Relative Permeability of Gas for Unconventional Reservoirs

Relative permeability of gas gains great significance in exploring unconventional gas. This paper developed a universal relative permeability model of gas, which is applicable for unconventional gas reservoirs such as coal, tight sandstone and shale. The model consists of the absolute relative permeability of gas and the gas slippage permeability. In the proposed model, the effects of water saturation and mean pore pressure on gas slippage permeability are taken into account. Subsequently, the evaluation of the model with existing model is done and then the validation of the model is made with data of tight sandstones, coals and shales from published literatures. The modeling results illustrate that a strong power-law relationship between relative permeability of gas and water saturation and the contribution of gas slippage permeability to relative permeability is determined by water saturation and mean pore pressure simultaneously. Furthermore, a sensitivity analysis of the impact of the parameters in the model is conducted and their effects are discussed.     

Drawdown Induced Changes in Permeability of Coalbeds: A New Interpretation of the Reservoir Response to Primary Recovery

A model for pore pressure-dependent cleat permeability is presented for gas-desorbing, linear elastic coalbeds under uniaxial strain conditions experienced in producing reservoirs. In the model, changes in the cleat permeability of coalbeds, which are idealised to have a bundled matchstick geometry, is controlled by the prevailing effective horizontal stresses normal to the cleats. Variations in the effective horizontal stresses under uniaxial strain conditions are expressed as a function of pore pressure reduction during drawdown, which includes a cleat compression term and a matrix shrinkage term that have competing effects on cleat permeability. A comprehensive analysis has revealed that the shape of the stress – pore pressure curve is predominantly determined by the magnitude of recovery pressure and rebound pressure relative to the initial reservoir pressure. A total of five possible scenarios have been identified with regard to response of the horizontal stress function to reservoir drawdown. When applied to four coalbed wells at two separate sites in the fairway of the San Juan basin, the model predictions at one site, where the three wells have shown increased absolute permeability during gas production, are in excellent agreement with the published pore pressure dependent permeability changes that were obtained independently from history matching the field production data. At a separate site the model correctly predicts, at least qualitatively, a strong permeability rebound at lower drawdown pressures that has been inferred through history matching the production data. An analysis of the effects of initial reservoir pressure on the response of effective horizontal stress to drawdown was carried out, with reference to the range of pressure likely to be encountered in the San Juan basin. The implications of this in terms of pore pressure dependent permeability are discussed.     

致密砂岩逆向渗吸作用距离实验研究

中国致密油储量丰富, 但多数致密储层波及效率低, 衰竭开发效果较差. 逆向渗吸是致密油藏注水开发过程中的一种重要的提高采收率途径, 目前许多学者主要针对致密油藏渗吸采收率及其影响因素开展研究, 而对于渗吸作用距离(表征致密油藏渗吸作用范围)研究较少. 本文采用CT在线扫描装置建立了致密岩心逆向渗吸作用距离量化方法, 明确了逆向渗吸的作用范围, 进一步研究了流体压力、含水饱和度、岩心渗透率和表面活性剂对逆向渗吸作用距离的影响, 阐明了逆向渗吸作用距离与渗吸采收率的关系, 为提高致密油藏采收率提供指导. 研究结果表明, 渗透率为0.3 mD的致密岩心逆向渗吸作用距离尺度仅为1.25 ~ 1.625 cm; 5 MPa条件下渗透率为0.302 mD的岩心逆向渗吸作用距离为1.375 cm. 在本实验条件下, 流体压力和初始含水饱和度对致密岩心逆向渗吸作用距离的影响较小, 而渗透率和表面活性剂对致密岩心逆向渗吸作用距离的影响显著, 渗透率为0.784 mD的岩心逆向渗吸作用距离相较于渗透率为0.302 mD的岩心提高2.63倍. 逆向渗吸作用距离是渗吸采收率表征的重要参数, 决定了逆向渗吸作用的波及范围.      

Unstable Immiscible Displacements in Oil-Wet Rocks

Displacement of a viscous fluid by a lower viscosity immiscible fluid (such as waterflood of a viscous oil) in a porous medium is unstable. The displacement front generates viscous fingers which lead to low oil recovery efficiency. These fingers are much smaller in width than typical reservoir simulation grid blocks, and capturing their effect in reservoir simulation is important. A dimensionless scaling group (viscous finger number) had been suggested in the past, which has a power-law relationship with the breakthrough recovery and cumulative recovery in unstable core floods. The relative permeability used in large grid block simulations had been modified to so-called pseudo-relative permeability on the basis of the dimensionless group, thus incorporating the effect of fingers in waterflood predictions. However, the previous proposed models were constructed from experiments in only water-wet rocks. This paper extends the recent viscous fingering models to oil-wet systems. Sandstone cores were treated to alter the wettability to oil-wet. Adverse viscosity water floods were performed in oil-wet cores. Viscosity ratio, velocity and diameter were varied. It is shown that the previously developed viscous finger number does not work for the oil-wet experiments. The correlating dimensionless number is modified for oil-wet systems; it is also different from the dimensionless group identified by Peters and Flock (Soc Petroleum Eng, 1981. doi: 10.2118/8371-PA) for oil-wet cores. A pseudo-relative permeability model has been developed for oil-wet cores. Corefloods have been matched by the new pseudo-relative permeability model to determine the model parameters. This pseudo-relative permeability model can be used in reservoir simulations of water and polymer floods in viscous oil-wet 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.     

A Numerical Study of Microscale Flow Behavior in Tight Gas and Shale Gas Reservoir Systems

Various attempts have been made to model flow in shale gas systems. However, there is currently little consensus regarding the impact of molecular and Knudsen diffusion on flow behavior over time in such systems. Direct measurement or model-based estimation of matrix permeability for these “ultra-tight” reservoirs has proven unreliable. The composition of gas produced from tight gas and shale gas reservoirs varies with time for a variety of reasons. The cause of flowing gas compositional change typically cited is selective desorption of gases from the surface of the kerogen in the case of shale. However, other drivers for gas fractionation are important when pore throat dimensions are small enough. Pore throat diameters on the order of molecular mean free path lengths will create non-Darcy flow conditions, where permeability becomes a strong function of pressure. At the low permeabilities found in shale gas systems, the dusty-gas model for flow should be used, which couples diffusion to advective flow. In this study we implement the dusty-gas model into a fluid flow modeling tool based on the TOUGH+ family of codes. We examine the effects of Knudsen diffusion on gas composition in ultra-tight rock. We show that for very small average pore throat diameters, lighter gases are preferentially produced at concentrations significantly higher than in situ conditions. Furthermore, we illustrate a methodology which uses measurements of gas composition to more uniquely determine the permeability of tight reservoirs. We also describe how gas composition measurement could be used to identify flow boundaries in these reservoir systems. We discuss how new measurement techniques and data collection practices should be implemented in order to take advantage of this method. Our contributions include a new, fit-for-purpose numerical model based on the TOUGH+ code capable of characterizing transport effects including permeability adjustment and diffusion in micro- and nano-scale porous media.     

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

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

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

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

Construction of Dual Pore 3-D Digital Cores with a Hybrid Method Combined with Physical Experiment Method and Numerical Reconstruction Method

In order to further improve the accuracy of digital core modeling, a new hybrid method was proposed to construct the dual pore 3-D digital core with high precision. First of all, 3-D macro-pore digital cores were constructed by micro-CT. Secondly, based on the high-resolution 2-D scanning electron microscopy images of rock cores, micro-pore digital cores were constructed by simulated annealing method. And then, a superposition method was used to construct the digital core which could describe different pore characteristics. Finally, pore structures of digital cores were compared, and lattice Boltzmann method was used to analyze the percolation properties. The results show that the carbonate dual pore digital core constructed by the new hybrid method has a high accuracy, which can capture the pore properties of both macro-pore and micro-pore, and whose permeability simulation results are in good agreement with the experimental measurements. In addition, the new hybrid method is not only accurate and reliable, but also high efficient and economic, and can be applied to all kinds of reservoir modeling.     

裂缝性低渗透油藏流-固耦合理论与数值模拟

根据裂缝性低渗油藏的储层特征,建立适合裂缝性砂岩油藏渗流的等效连续介质模型。将渗流力学与弹塑性力学相结合,建立裂缝性低渗透油藏的流-固耦合渗流数学模型,并给出其数值解.通过数值模拟对一实际井网开发过程中孔隙度、渗透率的变化以及开发指标进行计算,并和刚性模型以及双重介质模型的计算结果进行了分析比较.     

致密砂岩气藏纳微观结构及渗流特征

随着常规油气资源的减少,致密砂岩气藏逐渐成为勘探开发的热点.然而致密砂岩气高效产出机理研究还相对滞后,已经成为目前制约致密砂岩气大规模经济有效开发的瓶颈.致密砂岩由于其特殊的地质特征, 基质微观结构复杂,天然裂缝一定程度发育, 投产一般需要水力压裂等增产措施,气体在其中渗流存在跨尺度效应.研究揭示出致密砂岩气在多孔介质中的流动分为解吸、扩散、渗流等几个方式,包括浓度场下的扩散、压力场中的渗流等.综述了致密砂岩存在的非常规地质特征及其研究描述方法,并依据克努森数大小, 对致密砂岩气的渗流状态进行流态区域划分,为更好的理解和模拟致密砂岩气的流动状态提供了理论依据,并指出建立一种高效、简洁的微观结构精细描述方法,以及解决致密砂岩气在多尺度条件下的传质优化等是今后研究的重点方向.      

A Computational Approach to Determine Average Reservoir Pressure in a Coalbed Methane (CBM) Well Flowing Under Dominant Matrix Shrinkage Effect

Desorption of gas from coal matrix alters the pore volume of fracture network. Consequently, cleat porosity and permeability of reservoir changes as pressure depletes. The method of standard pressure analysis calculations produces incorrect results in the case of coalbed methane reservoirs producing under dominant matrix shrinkage effect. The change in cleat porosity and permeability due to shrinkage of coal matrix following gas desorption with pressure depletion invalidates the underlying assumptions made in the derivation of diffusivity equation. Consequently, equations of pseudo-steady state commonly used in conventional reservoirs no longer remain valid as the porosity and permeability values change with pressure depletion. In this paper, effort has been made to describe pseudo-steady-state flow in coalbed methane reservoirs in the form of a new equation that accounts for pressure dependency of cleat porosity and permeability due to shrinkage of coal matrix. The concept of Al-Hussainy et al. (1966) has been extended to define a new pseudo-pressure function which assimilates within itself the pressure dependence of porosity and permeability Palmer and Mansoori (1998). Equation has been used to relate the cleat porosity with pressure. The equation-based computational method suggested in this paper finds its usefulness in estimating average reservoir pressure for any known flowing bottom hole pressure and thus reducing the frequency of future pressure buildup tests. The new equation is also useful in predicting reservoir pressure under the situation when coal matrix shrinks below desorption pressure. The equation used in the computational method has been validated with the help of numerical simulator CMG-GEM.     

A Semi-Analytical Model for Two Phase Immiscible Flow in Porous Media Honouring Capillary Pressure

This article describes a semi-analytical model for two-phase immiscible flow in porous media. The model incorporates the effect of capillary pressure gradient on fluid displacement. It also includes a correction to the capillarity-free Buckley–Leverett saturation profile for the stabilized-zone around the displacement front and the end-effects near the core outlet. The model is valid for both drainage and imbibition oil–water displacements in porous media with different wettability conditions. A stepwise procedure is presented to derive relative permeabilities from coreflood displacements using the proposed semi-analytical model. The procedure can be utilized for both before and after breakthrough data and hence is capable to generate a continuous relative permeability curve unlike other analytical/semi-analytical approaches. The model predictions are compared with numerical simulations and laboratory experiments. The comparison shows that the model predictions for drainage process agree well with the numerical simulations for different capillary numbers, whereas there is mismatch between the relative permeability derived using the Johnson–Bossler–Naumann (JBN) method and the simulations. The coreflood experiments carried out on a Berea sandstone core suggest that the proposed model works better than the JBN method for a drainage process in strongly wet rocks. Both methods give similar results for imbibition processes.     

储层含水条件下致密砂岩/页岩无机质纳米孔隙气相渗透率模型

页岩及致密砂岩储层富含纳米级孔隙,且储层条件下页岩孔隙(尤其无机质孔隙)及致密砂岩孔隙普遍含水,因此含水条件下纳米孔隙气体的流动能力的评价对这两类气藏的产能分析及生产预测具有重要意义.本文首先基于纳米孔隙内液态水及汽态水热力学平衡理论,量化了储层孔隙含水饱和度分布特征;进一步在纳米孔隙单相气体传质理论的基础上,考虑了孔隙含水饱和度对气体流动的影响;最终建立了含水饱和度与气相渗透率的关系曲线. 基于本文岩心孔隙分布特征,计算结果表明:储层含水饱和度对气体流动能力的影响不容忽视,在储层含水饱和度20%的情况下,气相流动能力与干燥情况相比将降低约10%;在含水饱和度40% 的情况下,气相流动能力将降低约20%.      

Pore Network Modelling of Electrical Resistivity and Capillary Pressure Characteristics

In order to model petrophysical properties of hydrocarbon reservoir rocks, the underlying physics occurring in realistic rock pore structures must be captured. Experimental evidence showing variations of wetting occurring within a pore, and existence of the so-called 'non-Archie' behaviour, has led to numerical models using pore shapes with crevices (for example, square, elliptic, star-like shapes, etc.). This paper presents theoretical derivations and simulation results of a new pore space network model for the prediction of petrophysical properties of reservoir rocks. The effects of key pore geometrical factors such as pore shape, pore size distribution and pore co-ordination number (pore connectivity) have been incorporated into the theoretical model. In particular, the model is used to investigate the effects of wettability and saturation history on electrical resistivity and capillary pressure characteristics. The petrophysical characteristics were simulated for reservoir rock samples. The use of the more realistic grain boundary pore (GBP) shape allows simulation of the generic behaviour of sandstone rocks, with various wetting scenarios. The predictions are in close agreement with electrical resistivity and capillary pressure characteristics observed in experiments.     

An Experimental and Numerical Study of Relative Permeability Estimates Using Spatially Resolved $$T_1$$- z NMR

Relative permeability is a key characteristic describing flow properties of petroleum reservoirs, aquifers and water retention of soils. Various laboratory methods, typically categorised as steady-state, unsteady-state and centrifuge are used to measure relative permeability and may lead to different results. In recent years, 1D MRI, NMR \(T_2\) and \(T_1\) profiling have been applied for the characterisation of rock cores. It has been shown that spatially resolved NMR in conjunction with centrifuge technique may provide high-quality capillary pressure curves. Combining Burdine and Brooks–Corey models enables estimation of relative permeability from capillary pressure curves. This approach assumes a strong relationship between capillary pressure and relative permeability known to be complex. Here we compare a generalised approach of Green, which relies on saturation profiles set by various capillary drainage techniques, to a NMR relaxation approach. Comparisons are performed experimentally and numerically using three sandstone rocks to test the influence of rock morphology. The numerical part includes simulation of a centrifuge capillary drainage by applying morphological drainage transforms on high-resolution 3D tomograms. \(T_1\) responses along the sample are simulated using a random walk technique. The NMR relaxation-based approach is then compared to LBM simulated relative permeability and to experiment. The study confirms the applicability of NMR relaxation methods for relative permeability estimation of water-wet rocks and validates a numerical approach against experiment.     

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.     

具有井简储集的变形介质双渗模型的压力分析

为更好地研究碳酸盐油藏和低渗油的渗流问题，引入渗透率模数，考虑应力敏感地层中介质的变形，介质的双孔隙度、双渗秀率特征，同时考虑井筒储集的影响，建立新的数学模型。渗透率依赖于孔隙压力变化的流动方程是强非线性的，模型采用Douglas—Jones预估-校正法获得了无限大地层及有界封闭地层的数值解，形成了新的理论图版，并利用这些图版对模型中的有关参数进行了敏感性分析。     

具有井筒储集的变形介质双渗模型的压力分析

为更好地研究碳酸盐油藏和低渗油的渗流问题,引入渗透率模数,考虑应力敏感地层中介质的变形,介质的双孔隙度、双渗秀率特征,同时考虑井筒储集的影响,建立新的数学模型。渗透率依赖于孔隙压力变化的流动方程是强非线性的,模型采用Douglas-Jones预估-校正法获得了无限大地层及有界封闭地层的数值解,形成了新的理论图版,并利用这些图版对模型中的有关参数进行了敏感性分析。