裂缝和孔洞型储层孔隙模型的理论进展

对有洞的和裂缝型储层的分析已经成为一个热点,因而孔隙模型的研究近年得到了很好的发展.目前已经用双孔隙和三孔隙模型研究这类储层的特性并寻找估计孔隙指数的方法,以便计算含水饱和度.用串联或并联电阻网络模拟储层表明：双孔隙模型适用于基质与非连通孔洞储层以及裂缝和（或）连通孔洞的储层.三孔隙模型更适用于由基质、裂缝和不连通孔洞组成储层的岩石物理评价,对于当前碳酸盐岩和火成岩以及变质岩储层评价具有明显的指导意义.     

Nonlinear simultaneous inversion of pore structure and physical parameters based on elastic impedance

Carbonate reservoirs have complex pore structures, which not only significantly affect the elastic properties and seismic responses of the reservoirs but also affect the accuracy of the prediction of the physical parameters. The existing rockphysics inversion methods are mainly designed for clastic rocks, and the inversion objects are generally porosity and water saturation. The data used are primarily based on the elastic parameters, and the inversion methods are mainly linear approximations. To date, there has been a lack of a simultaneous pore structure and physical parameter inversion method for carbonate reservoirs. To solve these problems, a new Bayesian nonlinear simultaneous inversion method based on elastic impedance is proposed. This method integrates the differential effective medium model of multiple-porosity rocks, Gassmann equation,Amplitude Versus Offset(AVO) theory, Bayesian theory, and a nonlinear inversion algorithm to achieve the simultaneous quantitative prediction of the pore structure and physical parameters of complex porous reservoirs. The forward modeling indicates that the contribution of the pore structure, i.e., the pore aspect ratio, to the AVO response and elastic impedance is second only to that of porosity and is far greater than that of water saturation. The application to real data shows that the new inversion method for determining the pore structure and physical parameters directly from pre-stack data can accurately predict a reservoir's porosity and water saturation and can evaluate the pore structure of the effective reservoir.     

用核磁共振测井资料评价碳酸盐岩等复杂岩性储集层

对于碳酸盐岩等复杂岩性油气藏,由于其储集空间复杂、非均质性强等因素,用常规测井技术难以进行准确描述.核磁共振测井测量的对象是储层空间中的流体,因而可以直接用来划分储集层,而且能提供几乎不受岩性影响的孔隙度和渗透率等参数;同时,由于其T2分布表征了岩石的孔隙结构,所以可以根据T2分布形态判断有效裂缝和溶蚀孔洞.通过多口井的岩心对比和成像测井对比,研究出了一套用MRIL和CMR 的T2分布形态评价储集空间的方法;此外,在特定条件下,根据核磁计算的有效孔隙度和可动流体体积给出了一个计算含油饱和度的公式.利用这些核磁测井技术及其分析方法对车古20古潜山碳酸盐岩储集层裂缝及溶蚀孔洞发育特征进行准确描述,还计算了埕北潜山复杂复杂岩性油气藏的饱和度,对于储量计算具有重要意义.     

基于岩石物理参数转换模型的页岩气含气量评价新方法研究

页岩气在储层中主要以游离气和吸附气的形式存在,其含气量主要通过现场解吸和测井方法计算.由于页岩储层十分致密,并且常伴有微裂缝发育,这就造成了页岩储层的储集空间从纳米级孔隙到微米-毫米级的微裂缝并存的特殊情况.针对微孔和介孔段(50nm以内)的等温吸附方法,以及针对宏孔(50nm以上)的高压压汞方法在表征页岩孔隙结构上都存在着局限性.而核磁(NMR)T2分布本身虽然能够表征孔径分布,但是需要精确刻度横向弛豫时间与孔径的关系.本文在现有方法的基础上,结合上述方法的优势,采用拼接的方式表征页岩整体孔径分布,并将核磁T2分布刻度为孔径分布.基于横向弛豫时间和孔径的关系,结合T2-I转换模型,采用动态参数计算页岩储层含水饱和度和含气量,这为页岩储层含气量评价提供了新的可靠的方法.     

The nonlinear oil–water two-phase flow behavior for a horizontal well in triple media carbonate reservoir

Carbonate reservoir is one of the important reservoirs in the world. Because of the characteristics of carbonate reservoir, horizontal well has become a key technology for efficiently developing carbonate reservoir. Establishing corresponding mathematical models and analyzing transient pressure behaviors of this type of well-reservoir configuration can provide a better understanding of fluid flow patterns in formation as well as estimations of important parameters. A mathematical model for a oil–water two-phase flow horizontal well in triple media carbonate reservoir by conceptualizing vugs as spherical shapes are presented in this article. A semi-analytical solution is obtained in the Laplace domain using source function theory, Laplace transformation, and superposition principle. Analysis of transient pressure responses indicates that seven characteristic flow periods of horizontal well in triple media carbonate reservoir can be identified. Parametric analysis shows that water saturation of matrix, vug and fracture system, horizontal section length, and horizontal well position can significantly influence the transient pressure responses of horizontal well in triple media carbonate reservoir. The model presented in this article can be applied to obtain important parameters pertinent to reservoir by type curve matching.     

基于岩石导电效率理论建立碳酸盐岩储层饱和度的计算方法(英文)

由于孔隙型碳酸盐岩储层的复杂孔隙结构和强非均质性,岩石导电效率与含水孔隙度之间理想的线性关系相对复杂或并不存在,故前人基于岩石导电效率理论建立含水饱和度的计算公式并不完全适用于碳酸盐岩储层。基于岩石导电效率理论,推导了岩石导电效率的计算公式,阐明了岩石导电效率与含水孔隙度之间线性关系相对复杂或并不存在的根本原因,发现了岩石导电效率与电阻率之间的幂函数关系,分析了所建立含水饱和度计算公式误差的主控因素,得出了岩石导电效率的计算精度是该方法是否能推广应用的关键。与Archie公式相比,基于岩石导电效率理论建立的含水饱和度计算公式能更准确地计算研究靶区孔隙型碳酸盐岩储层的含水饱和度。在伊拉克某油田和印尼某气田3口井碳酸盐岩储层中的应用,表明当计算岩石导电效率的相对误差不大于0.1时,所计算储层含水饱和度的绝对误差不大于0.1,基本满足孔隙型碳酸盐岩储层精细评价的需求。     

塔河油田碳酸盐岩缝洞型储层的测井识别与评价方法研究

塔河油田奥陶系以碳酸盐岩为主，油气的主要储渗空间为裂缝和溶蚀孔洞，具有很强的非均质性．本文利用常规及成像测井资料，对碳酸盐岩缝洞型储层的识别与评价方法进行研究．为了综合各种测井方法识别裂缝，建立了综合裂缝概率模型，计算综合裂缝概率指示裂缝的发育程度．利用地层微电阻率扫描成像测井资料进行裂缝和溶蚀孔洞的定性、定量解释．定量计算的裂缝参数为：裂缝密度、裂缝长度、裂缝平均宽度、平均水动力宽度、裂缝视孔隙度；定量计算的溶蚀孔洞参数有：面孔率、孔洞密度．根据缝洞型储层孔隙空间类型及其中子孔隙度、补偿密度、声波、双侧向电阻率的测井响应物性特征，建立缝洞型碳酸盐岩储层复杂孔隙介质解释模型，用于确定裂缝、溶蚀孔洞孔隙度和评价储层．     

基于小波变换图像分割技术的电成像测井资料裂缝、孔洞面孔率提取方法

对于裂缝、溶蚀孔洞发育的碳酸盐岩缝洞储层,如何从测井资料中提取裂缝、溶蚀孔洞信息是评价储层有效性的关键问题.为了从电成像测井静态图像上准确地分割出清晰的裂缝、溶蚀孔洞子图像进而提取其参数信息,本文提出了一种基于小波变换模极大值图像分割技术的电成像测井资料缝洞面孔率提取方法.以钮扣电极电导率曲线为对象,先消除井壁凹凸不平导致的地层背景噪声的影响,利用小波变换模极大值图像分割方法得到包含裂缝和溶蚀孔洞目标的子图像,根据子图像提取裂缝-孔洞总面孔率、裂缝面孔率、孔洞面孔率等信息.本文利用塔里木盆地奥陶系碳酸盐岩地层的电成像测井数据提取了缝洞面孔率参数,还利用同井岩心CT扫描图像计算的平均缝洞面孔率、双侧向电阻率、常规测井资料三孔隙度模型计算的相对连通缝洞孔隙度进行了对比,并进行了试油验证.对比表明,电成像测井裂缝-孔洞总面孔率、裂缝面孔率、孔洞面孔率与岩心CT扫描图像平均缝洞面孔率、双侧向电阻率、相对连通缝洞孔隙度、试油结果均有较好的一致性.这一方面验证了采用本文方法提取的裂缝-孔洞总面孔率、裂缝面孔率、孔洞面孔率的合理性,另一方面表明所提取参数可用于指示缝洞型碳酸盐岩储层的渗透性和有效性.     

Modeling and interpretation of Q logs in carbonate rock using a double porosity model and well logs

Attenuation data extracted from full waveform sonic logs is sensitive to vuggy and matrix porosities in a carbonate aquifer. This is consistent with the synthetic attenuation (1 / Q) as a function of depth at the borehole-sonic source-peak frequency of 10 kHz. We use velocity and densities versus porosity relationships based on core and well log data to determine the matrix, secondary, and effective bulk moduli. The attenuation model requires the bulk modulus of the primary and secondary porosities. We use a double porosity model that allows us to investigate attenuation at the mesoscopic scale. Thus, the secondary and primary porosities in the aquifer should respond with different changes in fluid pressure. The results show a high permeability region with a Q that varies from 25 to 50 and correlates with the stiffer part of the carbonate formation. This pore structure permits water to flow between the interconnected vugs and the matrix. In this region the double porosity model predicts a decrease in the attenuation at lower frequencies that is associated with fluid flowing from the more compliant high-pressure regions (interconnected vug space) to the relatively stiff, low-pressure regions (matrix). The chalky limestone with a low Q of 17 is formed by a muddy porous matrix with soft pores. This low permeability region correlates with the low matrix bulk modulus. A low Q of 18 characterizes the soft sandy carbonate rock above the vuggy carbonate.This paper demonstrates the use of attenuation logs for discriminating between lithology and provides information on the pore structure when integrated with cores and other well logs. In addition, the paper demonstrates the practical application of a new double porosity model to interpret the attenuation at sonic frequencies by achieving a good match between measured and modeled attenuation.     

Unsteady inter-porosity flow modeling for a multiple media reservoir

The paper deals with unsteady inter-porosity flow modeling of underground fluid in a multiple media reservoir. Assuming spherical vugs, symmetrically distributed pressure, negligible inter-porosity flow between matrix and vug systems and centrifugal flow of the fluid from matrix blocks or vugs to fractures, and treating media directly connected with wellbore as the fracture system, we establish and solve a model of unsteady inter-porosity flow for dual and triple porosity media reservoirs. We provide simulated graphs of pressure and pressure derivative log-log type curves, and analyze the transient flow process and characteristics of type curves affected by different parameters. The new type curves of unsteady inter-porosity flow modeling are evidently different in shape and characteristics from those of pseudo-steady inter-porosity flow modeling. The location of dimensionless pressure of unsteady inter-porosity is lower than that of pseudo-steady inter-porosity, which indicates that unsteady inter-porosity flow accelerates an energy supplement during production. Qualitatively, the unsteady inter-porosity flow modeling reduces the classical V-shaped response. We also estimated parameters from well test data in real applications using this model.     

CEC比值法在凝灰质砂岩储层测井评价中的应用——以海-塔盆地X凹陷凝灰质砂岩储层为例

凝灰质砂岩储层一直是测井评价中的难点.凝灰质的存在使储层的孔隙结构和物性变化都很大,其含量直接影响储层参数的求取,对储层饱和度的影响不容忽视.本文以海-塔盆地X凹陷的凝灰质砂岩储层为例进行测井评价研究.基于泥质和凝灰质的测井响应差异利用粒子群和细菌觅食的混合优化算法计算储层中两者的含量,然后通过阳离子交换量(CEC)的实验数据验证了凝灰质具有导电性,进而利用CEC与电阻率之间的关系得到凝灰质电阻率的计算方法,并且应用到饱和度的计算,最终得到了一种新的计算凝灰质砂岩储层饱和度的方法——CEC比值法,并取得了良好的应用效果.     

地层模量分解及在流体识别中的应用

储层流体识别是确定油气水分布,合理布设井位,提高钻井成功率的关键之一.本文基于流体饱和孔隙介质岩石物理模型,对地震反演的地层体积模量进行分解,获得孔隙流体体积模量,并依据油、气、水(尤其是气-油、气-水)模量的显著差异进行识别.文中简要分析了Gassmann模型和Kuster-Toksz模型的特征,详细讨论了孔隙形态和饱和度对弹性模量的影响,提出了联合Kuster-Toksz方程和Gassmann方程的体积模量分解方法.该方法通过Kuster-Toksz方程从测井数据中反演地层骨架固体和干骨架的弹性模量,再利用Gassmann方程对地层体积模量进行分解,既考虑了孔隙形态,又充分利用了Gassmann方程的易用性.理论模型结果表明方法是可行的.方法应用于西部地区某气田,流体识别与地层含气性预测结果与钻井基本一致,进一步证实了方法的有效性.     

复杂孔隙储层三维岩石物理模版

复杂孔隙储层往往同时发育孔缝洞等多种孔隙类型,这种孔隙结构的复杂性使得岩石的速度与孔隙度之间的相关性很差.经典的二维岩石物理模版只研究弹性参数与孔隙度和饱和度之间的定量关系,而不考虑孔隙结构的影响,用这样的模版来预测复杂孔隙储层的物性参数时带来很大偏差.本文首先证明多重孔隙岩石的干骨架弹性参数可以用一个等效孔隙纵横比的单重孔隙岩石物理模型来模拟;进而基于等效介质岩石物理理论和Gassmann方程,建立一个全新的三维岩石物理模版,用它来建立复杂孔隙岩石的弹性性质与孔隙扁度及孔隙度和饱和度之间的定量关系;在此基础上,预测复杂储层的孔隙扁度、孔隙度以及孔隙中所包含的流体饱和度.实际测井和地震反演数据试验表明,三维岩石物理模版可有效提高复杂孔隙储层参数的预测精度.     

Characterizing moldic and vuggy pore space in karst aquifers using borehole-wall,slabbed-core and thin-section images

Carbonate aquifers are prolific and important sources of potable water in many parts of the world owing to enlarged dissolution features that enhance porosity and interconnectivity. To better understand the variations of pore space in different karst aquifers, image and geospatial analyses are used to analyze pore attributes (i.e., pore area and perimeter) in images of vuggy aquifers. Pore geometry and 2D porosity derived from images of the moldic Castle Hayne and vuggy Biscayne aquifers are analyzed at three scales of observation: borehole televiewer, core and thin-section. The Castle Hayne and Biscayne aquifers are the foci of this study because the pore spaces that control the hydrologic properties in each of these aquifers are markedly different even though both of these carbonate reservoirs are prolific aquifers. Assessments of pore area, perimeter and shape index (a measure of shape complexity) indicate that pore geometries and pore complexities vary as a function pore type and scale of observation. For each aquifer type, the areas, perimeters and complexities of pores are higher at the larger scale of observation (e.g., borehole) than the smaller scale of observation (e.g., thin section). When the complexity of the moldic pores is compared to the complexity of vuggy pores, the results indicate that moldic pores are generally more complex than vuggy pores at the same scale of observation. Whereas estimates of 2D porosity from the borehole televiewer image of the vuggy aquifer are higher than those derived from the moldic aquifer, the range of 2D porosities is larger in core and thin section images for the vuggy aquifer than the moldic aquifer. A model for the development of pores is presented that suggests that the coalescence of small pores with simple shapes leads to the growth of larger pores with more complex shapes. The model suggests that the younger Biscayne aquifer is a more mature karst than the Castle Hayne aquifer.     

Modulus defect,velocity dispersion and attenuation in partially-saturated reservoirs of Jurassic sandstone,Indus Basin,Pakistan

Partially saturated reservoirs are one of the major sources of seismic wave attenuation, modulus defect and velocity dispersion in real seismic data. The main attenuation and dispersion phenomenon is wave induced fluid flow due to the heterogeneity in pore fluids or porous rock. The identification of pore fluid type, saturation and distribution pattern within the pore space is of great significance as several seismic and petrophysical properties of porous rocks are largely affected by fluid type, saturation and fluid distribution pattern. Based on Gassmann-Wood and Gassmann- Hill rock physics models modulus defect, velocity dispersion and attenuation in Jurassic siliclastic partially-saturated rocks are studied. For this purpose two saturation patterns - uniform and patchy - are considered within the pore spaces in two frequency regimes i.e., lower frequency and higher frequency. The results reveal that at low enough frequency where saturation of liquid and gas is uniform, the seismic velocity and bulk modulus are lower than at higher frequency where saturation of fluid mixture is in the form of patches. The velocity dispersion and attenuation is also modeled at different levels of gas saturation. It is found that the maximum attenuation and velocity dispersion is at low gas saturation. Therefore, the dispersion and attenuation can provide a potential way to predict gas saturation and can be used as a property to differentiate low from high gas saturation.     

An insightful parametrization for the flatlander's interpretation of time‐lapsed seismic data

An approximation is developed that allows mapped 4D seismic amplitudes and time‐shifts to be related directly to the weighted linear sum of pore pressure and saturation changes. The weights in this relation are identified as key groups of parameters from a petroelastic model and include the reservoir porosity. This dependence on groups of parameters explains the inherent non‐uniqueness of this problem experienced by previous researchers. The proposed relation is of use in 4D seismic data feasibility studies and inversion and interpretation of the 4D seismic response in terms of pore pressure and water saturation changes. A further result is drawn from analysis of data from the North Sea and West Africa, which reveals that the relative interplay between the effects of pore pressure and saturation changes on the seismic data can be simplified to the control of a single, spatially variant parameter C_S/C_P. Combining these results with those from published literature, we find that C_S/C_P = 8 appears to be a generality across a range of clastic reservoirs with a similar mean porosity. Using this C_S/C_P value, an in situ seismic‐scale constraint for the rock stress sensitivity component of the petroelastic model is constructed considering this component carries the largest uncertainty.     

Facies-constrained prestack seismic probabilistic inversion driven by rock physics

Seismic Rock physics plays a bridge role between the rock moduli and physical properties of the hydrocarbon reservoirs. Prestack seismic inversion is an important method for the quantitative characterization of elasticity, physical properties, lithology and fluid properties of subsurface reservoirs. In this paper, a high order approximation of rock physics model for clastic rocks is established and one seismic AVO reflection equation characterized by the high order approximation(Jacobian and Hessian matrix) of rock moduli is derived. Besides, the contribution of porosity, shale content and fluid saturation to AVO reflectivity is analyzed. The feasibility of the proposed AVO equation is discussed in the direct estimation of rock physical properties. On the basis of this, one probabilistic AVO inversion based on differential evolution-Markov chain Monte Carlo stochastic model is proposed on the premise that the model parameters obey Gaussian mixture probability prior model. The stochastic model has both the global optimization characteristics of the differential evolution algorithm and the uncertainty analysis ability of Markov chain Monte Carlo model. Through the cross parallel of multiple Markov chains, multiple stochastic solutions of the model parameters can be obtained simultaneously, and the posterior probability density distribution of the model parameters can be simulated effectively. The posterior mean is treated as the optimal solution of the model to be inverted.Besides, the variance and confidence interval are utilized to evaluate the uncertainties of the estimated results, so as to realize the simultaneous estimation of reservoir elasticity, physical properties, discrete lithofacies and dry rock skeleton. The validity of the proposed approach is verified by theoretical tests and one real application case in eastern China.     

利用注聚合物泥浆伤害岩电实验参数变化评价特低渗透储层

针对鄂尔多斯盆地陕北斜坡长6特低渗透储层受聚合物泥浆伤害,微电极电阻率曲线在渗透层上的正幅度差异不明显,直观指示油气层和水层的深、中,浅探测电阻率在常规储层的有序排列基本消失.通过注聚合物泥浆驱替前后岩石电阻率与孔隙度及含水饱和度实验关系分析,驱替后岩样电阻率显著增大,岩样电阻率随含水饱和度增大而减小的基本规律不复存在.归纳其岩性系数a、电阻率系数b比驱替前大得多,饱和度指数n由正变负,反映出注聚合物泥浆驱替破坏了储层孔隙结构及其基本特征,测井中造成了范围较小的高侵和特高侵地层电阻率带.从而提出在较为致密的低渗砂、致密砂聚合物泥浆伤害储层中,利用微电极负差异及其电阻率曲线不规则增高变化划分特低渗透油层有效厚度,并以实例阐明了利用岩电实验参数变化评价特低渗透储层的方法.     

基于核磁共振实验的低渗透砂岩岩电参数分类及应用——以东营凹陷南坡沙四段为例

低渗透砂岩储层孔隙结构复杂,储层有效性识别及饱和度准确计算难度较大.笔者以东营凹陷南坡沙四段（Es₄）低渗透砂岩为研究对象,根据压汞、物性、薄片及核磁等资料,将研究区孔隙结构分为三大类、五小类.在岩样孔隙结构分类基础上,明确了孔隙结构类型与岩电参数之间存在确定的关系,而核磁共振T₂谱定量特征参数在一定程度上能够表征孔隙结构类型及其细节信息,通过提取T₂谱中T₂几何平均值（T_2g）、T₂均值（T₂）、峰度（K_G）、可动流体分量（S_mf）及区间孔隙分量等孔隙结构参数,建立了基于核磁T₂谱特征参数的孔隙结构识别图版,显示核磁T₂谱孔隙结构参数对不同类型的储层有较好的识别效果,进而探讨核磁孔隙结构参数和岩电参数之间的关系,结果表明,T₂谱峰度值与孔隙胶结指数（m）值相关性较高,进一步确定了岩电参数m的核磁计算公式.最终,将该套方法应用于研究区井筒剖面中,有效地提高了饱和度计算精度,也为东营凹陷南坡低渗透砂岩油藏储量估算与高效开发提供了依据.     

碳酸盐岩储层孔隙结构对电阻率的影响研究

碳酸盐岩储层孔隙类型多样,各种孔隙的尺寸变化范围可以跨越几个数量级,孔隙结构非常复杂,这种复杂孔隙结构和不均匀分布的多元孔隙空间使得储层电性呈现明显非阿尔奇特性.为了了解影响电阻率变化的控制因素,本次研究选取中三叠世雷口坡组的8块全直径碳酸盐岩岩样,开展了核磁共振、岩电实验、孔渗实验、压汞实验及薄片等实验,并利用数字图像分析法定量分析了孔隙结构特征.研究结果表明：①孔隙度是影响电阻率高低的重要因素,但并非唯一因素,除孔隙度以外,孔隙尺寸和数量、孔隙网络复杂程度远比吼道大小对电阻率的影响大;②在孔隙度一定的条件下,胶结指数m随储层中孤立大孔隙占比的增多而增大,当孔隙度增大到一定程度后,胶结指数m又随大孔隙占比的增多而减小,微裂缝起重要沟通作用;③在给定孔隙度时,以简单大孔隙为主的岩样表现为胶结指数m值较大,而以复杂孔隙网络、细小孔隙为主的岩样表现为胶结指数m值较小,具分散、孤立大孔隙的岩样,胶结指数m值最高;④依据孔隙几何参数与电阻率和胶结指数之间的关系,可以利用测井资料间接判别储层类型,从而提高储层有效性和含水饱和度评价精度.