裂缝型储层流体识别方法

裂缝型储层的描述包括预测裂缝分布特征和识别裂隙充填物.依据等效介质理论计算的纵波速度随裂缝密度的增大而减小.正演地震记录显示,裂缝介质含气时反射振幅最大,且变化程度比含油或含水时大.叠前方位AVO反演所得的各向异性梯度Bani与裂缝密度成正比,可用于描述有效裂缝发育强度.对于不同的裂缝密度,各向异性梯度Bani与各向同性梯度Biso的比值I(1)fluid近似为常数,且对流体敏感.经裂缝纵横比和背景介质拉梅常数修正后,流体因子Ifluid既不随纵横比变化,又不受背景介质的影响,是裂缝型储层敏感的流体识别因子.在塔里木盆地塔北哈拉哈塘地区热瓦普区块碳酸盐岩储层裂缝发育区域,运用该参数在井点处的流体识别效果与钻井结果一致.     

基于各向异性岩石物理的缝隙流体因子AVAZ反演

裂缝型储层表现出较强的各向异性特征.缝隙中充填不同流体时,裂缝储层的地震响应特征也不相同.本文从各向异性岩石物理模型出发,引入可有效识别缝隙流体的指示因子,并研究缝隙充填流体类型、饱和度以及缝隙纵横比与流体因子的相互关系,进而分析不同流体充填时介质的地震响应特征,并基于AVAZ反演方法估测缝隙流体指示因子.首先对缝隙流体因子的敏感性进行了分析,讨论当缝隙充填不同流体时,缝隙流体因子值的变化特征,同时研究了不同流体类型充填时裂缝储层反射系数随方位角和入射角的变化特征.某工区测井数据和复杂裂缝模型应用表明,基于各向异性岩石物理的缝隙流体因子AVAZ反演方法合理、可靠,且具有良好的抗噪性,即当对合成地震记录添加信噪比不小于1/2的随机噪声时,利用AVAZ反演方法估测所得流体因子值与真实值仍然吻合较好.     

流体饱和度对裂缝—孔隙岩石频变纵波速度及其各向异性的影响

流体饱和度会改变裂缝性储层的纵波速度,从而影响地层速度的频散特性及各向异性程度,导致储层地震响应特征复杂,储层预测多解性强,流体识别难度大.本文根据多相流体饱和裂缝性储层的特点,借助于Norris和KG模型,建立部分饱和裂缝-孔隙等效介质模型,给出频变地震波速度随流体饱和度变化的精确关系式.数值模拟结果表明,当气、水两相共存时,随着含水饱和度的增加,高频段纵波相速度逐渐增大,各向异性程度逐渐减小;低频段纵波相速度逐渐减小,各向异性程度不变;相速度频散及其各向异性程度逐渐增强.组合已有的孔隙弹性理论模型,对实验室人工裂缝-孔隙砂岩岩样的纵波速度进行拟合,计算得到的曲线与实验室测量散点值吻合度较高,表明组合模型在给定参数下的有效性.该研究能够为多相流体饱和裂缝性储层的地震响应特征分析奠定扎实的理论基础,为提高储层预测的确定性和流体识别的准确性提供可靠的理论依据.     

弱各向异性介质中P-SV波弹性阻抗与流体识别因子(英文)

基于Petr Jilek正交各向异性介质中的P-SV波反射系数精确公式,通过近似得到VTI和HTI介质中的P-SV波近似公式,同时结合Kenneth横波弹性阻抗原理,推导建立了VTI和HTI介质中P-SV波的弹性阻抗近似公式;在此基础上构建了弱各向异性介质中的流体识别因子公式,通过理论模型对Castagna和Smith的岩性组合进行了识别;最后详细分析了各向异性参数的取值对P-SV波弹性阻抗及识别因子产生的影响及变化趋势。     

基于挤喷流效应的Russell流体因子推广及应用

在油气勘探与开发过程中,寻求能够从地震资料中直接识别储层油气的流体指示因子至关重要.由于现存多数流体指示因子都是在Biot理论假设前提下建立起来的,因此在双相孔隙介质条件下不能有效地识别流体.为此,本文基于前人提出的双孔介质统一波动理论,考虑岩石裂缝间挤喷流效应,在经典流体指示因子基础上构建了一种新的流体指示因子.该流体因子能够较好地反映岩石内孔隙流体的变化对波传播产生的影响.在实际资料的应用中,采用多种方法分析、对比该流体识别因子与以Biot理论为基础的传统流体因子的优劣,理论分析与实际资料的验证表明该流体因子对于储层中油气的检测有较高的精度和灵敏度.     

流体因子关系分析以及新流体因子的构建

目前流体因子(ρf),泊松阻抗和,λρ,μρ等属性常用来进行流体识别和岩性预测.这些方法都能在流体识别中发挥作用,但是由于信息量比较单一,所以在地质情况比较复杂的时候不是很有效.本文以流体因子(ρf)和泊松阻抗(PI)为基础,构建了新的流体因子公式,提高了流体识别和岩性预测的敏感性.将构建的新流体因子属性应用到数值模拟...     

背景为各向异性的含裂缝岩石频散和衰减计算方法研究

裂缝广泛分布于各类储层岩石中,并且会显著提高储层的渗流能力.因此,裂缝的评价和表征对于提高油气产能具有重要意义.由于裂缝与背景介质之间的波致流会显著影响地震波的频散和衰减特性,所以地震勘探是评价裂缝性储层的有效手段.裂缝地震定量表征的前提是要基于含裂缝岩石中波致流对频散和衰减的影响建立含裂缝岩石物理特性与地震性质的关系.然而,目前相关的理论研究大部分基于各向同性背景这一假设,难以有效应用于常见的各向异性储层.本文针对背景为各向异性的含裂缝岩石提出了频散和衰减的计算方法.该方法首先将含裂缝岩石中的各向异性背景介质等效为层状背景介质;然后,通过分析不同频率下层状含裂缝岩石中的流体压力分布,理论计算了两个特定的中间频率并求解得到两个中间频率下的弹性参数;进一步,以计算得到的两个特定中间频率以及高低频极限下的弹性参数为基础,应用数值方法求解得到弛豫函数中的未知参数,最终实现了背景为各向异性含裂缝岩石中频散和衰减的理论模拟.通过将理论预测结果与实验测量和数值模拟结果进行对比,验证了该方法在背景为各向异性含不同分布裂缝岩石中的有效性.本文提出的方法考虑了常见的各向异性背景对含裂缝岩石频散和衰减的影...     

地震资料中的多尺度断裂与流体响应特征：数值模拟分析

流体在断裂和岩石骨架间的交换被认为是影响岩石弹性参数各向异性的主要原因,理论研究表明断裂尺度同样对弹性参数的各向异性也有影响.为了说明两者对各向异性影响以实现多尺度断裂裂隙的识别,本文在等效介质模型的基础上,应用数值分析的方法研究速度和衰减(1/Q)随多尺度断裂、频率和流体因子变化规律.结果表明介质弹性参数是频率依赖的,并且参数中存在衰减项,而这种频率依赖性与介质物性参数中的断裂尺度及流体性质存在一定的联系;当断裂定向分布时,参数结果显示为各向异性;不同断裂尺度具有不同的波速频散特性,剪切波分裂程度依赖于频率,断裂尺度起着控制作用,高频时对小尺度的敏感,低频段对大尺度敏感.在地震频段Thomsen参数随着频率的增大而减小,随着断裂尺寸的增大而减小.因此地震数据可能区分断裂和微裂隙引起各向异性,从而可测量断裂尺度.     

利用叠前地震数据预测裂缝储层的应用研究

叠前地震资料储层预测技术是在Zoeppritz方程基础上发展起来的,通过处理地震数据随着不同入射角地震反射属性,得到地震属性随着入射角变化而改变,研究分析得到反映岩性变化的纵波速度、横波速度、泊松比和截距与梯度剖面,预测裂缝储层的发育及分布.同时根据不同方位角地震资料属性,计算得到不同方位角目的层的属性差异,使用各向异性椭圆公式作拟合,求出背景趋势A和各向异性因子B,利用最大振幅包络方位和对应θ,求出裂缝发育优势方向,及B/A各向异性因子,实现对裂缝储层预测.     

基于地震资料的三种岩性流体预测方法对比分析

利用地震资料识别储层流体是储层预测的一项重要内容,对储层流体的识别是通过流体识别因子来实现的.首先引入三种流体识别因子并介绍其预测流体的方法原理,在此基础上以湿砂岩和含气砂岩为例,对比分析了不同流体识别因子对流体的敏感度,然后通过一个具体实例进一步验证了不同方法的特点.分析结果表明,第三种方法对流体的敏感度最高,能更准确地表征岩石孔隙中流体的性质.     

Fluid identification based on P-wave anisotropy dispersion gradient inversion for fractured reservoirs

Fluid identification in fractured reservoirs is a challenging issue and has drawn increasing attentions. As aligned fractures in subsurface formations can induce anisotropy, we must choose parameters independent with azimuths to characterize fractures and fluid effects such as anisotropy parameters for fractured reservoirs. Anisotropy is often frequency dependent due to wave-induced fluid flow between pores and fractures. This property is conducive for identifying fluid type using azimuthal seismic data in fractured reservoirs. Through the numerical simulation based on Chapman model, we choose the P-wave anisotropy parameter dispersion gradient (PADG) as the new fluid factor. PADG is dependent both on average fracture radius and fluid type but independent on azimuths. When the aligned fractures in the reservoir are meter-scaled, gas-bearing layer could be accurately identified using PADG attribute. The reflection coefficient formula for horizontal transverse isotropy media by Rüger is reformulated and simplified according to frequency and the target function for inverting PADG based on frequency-dependent amplitude versus azimuth is derived. A spectral decomposition method combining Orthogonal Matching Pursuit and Wigner–Ville distribution is used to prepare the frequency-division data. Through application to synthetic data and real seismic data, the results suggest that the method is useful for gas identification in reservoirs with meter-scaled fractures using high-qualified seismic data.     

Estimation of Dry Fracture Weakness,Porosity, and Fluid Modulus Using Observable Seismic Reflection Data in a Gas-Bearing Reservoir

Fracture detection and fluid identification are important tasks for a fractured reservoir characterization. Our goal is to demonstrate a direct approach to utilize azimuthal seismic data to estimate fluid bulk modulus, porosity, and dry fracture weaknesses, which decreases the uncertainty of fluid identification. Combining Gassmann’s (Vier. der Natur. Gesellschaft Zürich 96:1–23, 1951) equations and linear-slip model, we first establish new simplified expressions of stiffness parameters for a gas-bearing saturated fractured rock with low porosity and small fracture density, and then we derive a novel PP-wave reflection coefficient in terms of dry background rock properties (P-wave and S-wave moduli, and density), fracture (dry fracture weaknesses), porosity, and fluid (fluid bulk modulus). A Bayesian Markov chain Monte Carlo nonlinear inversion method is proposed to estimate fluid bulk modulus, porosity, and fracture weaknesses directly from azimuthal seismic data. The inversion method yields reasonable estimates in the case of synthetic data containing a moderate noise and stable results on real data.     

针对深层流体识别的两项弹性阻抗反演与Russell流体因子直接估算方法研究

Russell流体因子是判识深层储层含流体特征的敏感指示参数.针对深层地震资料缺失大角度入射信息以及信噪比普遍较低的特点,开展了基于贝叶斯理论框架的两项弹性阻抗反演与Russell流体因子直接估算方法研究.首先,在饱和孔隙弹性介质理论的指导下,推导了包含Russell流体因子和纵波阻抗参数的两项反射系数近似公式,并且进一步得到了新的两项弹性阻抗方程.其次通过模型的近似精度分析,验证了新弹性阻抗方程满足小角度近似精度要求.最后提出了贝叶斯弹性阻抗反演和Russell流体因子的直接估算方法,通过模型试算表明直接估算方法不仅充分利用了弹性阻抗反演的高抗噪性与实用性,并且有效的避免了间接累计误差对流体因子计算精度的影响,实际应用则进一步验证了直接估算的Russell流体因子可以提高深层储层流体识别的可靠性.     

Numerical analysis of seismic wave propagation in fluid-saturated porous multifractured media

Elastic wave propagation and attenuation in porous rock layers with oriented sets of fractures, especially in carbonate reservoirs, are anisotropic owing to fracture sealing, fracture size, fracture density, filling fluid, and fracture strike orientation. To address this problem, we adopt the Chapman effective medium model and carry out numerical experiments to assess the variation in P-wave velocity and attenuation, and the shear-wave splitting anisotropy with the frequency and azimuth of the incident wave. The results suggest that velocity, attenuation, and anisotropy vary as function of azimuth and frequency. The azimuths of the minimum attenuation and maximum P-wave velocity are nearly coincident with the average strike of the two sets of open fractures. P-wave velocity is greater in sealed fractures than open fractures, whereas the attenuation of energy and anisotropy is stronger in open fractures than sealed fractures. For fractures of different sizes, the maximum velocity together with the minimum attenuation correspond to the average orientation of the fracture sets. Small fractures affect the wave propagation less. Azimuth-dependent anisotropy is low and varies more than the other attributes. Fracture density strongly affects the P-wave velocity, attenuation, and shear-wave anisotropy. The attenuation is more sensitive to the variation of fracture size than that of velocity and anisotropy. In the seismic frequency band, the effect of oil and gas saturation on attenuation is very different from that for brine saturation and varies weakly over azimuth. It is demonstrated that for two sets of fractures with the same density, the fast shear-wave polarization angle is almost linearly related with the orientation of one of the fracture sets.     

Azimuthal anisotropy analysis of walkaround vertical seismic profiling vertical seismic profiling: a case study from Saudi Arabia

Understanding fracture orientations is important for optimal field development of fractured reservoirs because fractures can act as conduits for fluid flow. This is especially true for unconventional reservoirs (e.g., tight gas sands and shale gas). Using walkaround Vertical Seismic Profiling (VSP) technology presents a unique opportunity to identify seismic azimuthal anisotropy for use in mapping potential fracture zones and their orientation around a borehole. Saudi Aramco recently completed the acquisition, processing and analysis of a walkaround VSP survey through an unconventional tight gas sand reservoir to help characterize fractures. In this paper, we present the results of the seismic azimuthal anisotropy analysis using seismic traveltime, shear‐wave splitting and amplitude attenuation. The azimuthal anisotropy results are compared to the fracture orientations derived from dipole sonic and image logs. The image log interpretation suggests that an orthorhombic fracture system is present. VSP data show that the P‐wave traveltime anisotropy direction is NE to SW. This is consistent with the cemented fractures from the image log interpretation. The seismic amplitude attenuation anisotropy direction is NW to SE. This is consistent with one of the two orientations obtained using transverse to radial amplitude ratio analysis, with the dipole sonic and with open fracture directions interpreted from image log data.     

瞬时频域衰减属性及其在碳酸盐岩裂缝检测中的应用

基于AVO梯度、走时及速度的纵波方位角裂缝检测已有较多研究和应用.由于裂缝及所含流体的综合作用,碳酸岩裂缝性储层的地震响应具有高频衰减、频带低移及方位各向异性的特点,因此应用衰减属性进行裂缝预测具有重要价值.基于小波变换计算的瞬时频域振幅谱及以此为基础提取的衰减属性,由于在时域和频域都有较好的分辨率能力,能够准确地刻画出碳酸岩裂缝性储层的非均质性.结合三维地震纵波数据的分方位角处理和分析技术,可以由方位衰减属性估算裂缝相对密度.对Z区下古生界潜山碳酸岩储层进行的裂缝预测表明,由衰减梯度和85%能量比频率两种属性导出的裂缝发育密度具有很好的相关性,且与实际情况符合.     

富含煤层气储层的叠前地震反演预测方法研究——以山西晋东南测区为例

叠前地震反演技术能够将地震中振幅、偏移距和入射角等多种有效信息与油气的敏感参数结合起来,在有效识别油气藏"甜点"方面虽发挥了重要的作用,但对于煤层气的有效预测仍需研究.针对煤层气与石油天然气的赋存地质条件不同,本文以山西沁水盆地煤层气地震资料为例,借鉴叠前地震反演技术的多参数预测思想,通过模型正演和实测资料处理,得到纵波、横波速度和密度等不同的数据体,利用数据体和弹性参数关系式,采用AkiRichard近似公式与叠前宽角度反演联合的方法,寻找出富含煤层气储层的敏感参数.从而实现富含煤层气储层的有效预测.经与测井曲线结果对比,吻合率较高,说明了该方法具有可行性与实用性.     

基于贝叶斯线性AVAZ的TTI介质裂缝参数反演

裂缝储层岩石物理参数的准确获得对地下裂缝预测具有重要意义,而叠前方位AVA地震反演是获得裂缝岩石物理参数的有效手段.假设地下岩石为倾斜横向各向同性（TTI）介质,本文从裂缝岩石物理等效模型的构建出发,从测井数据中估计出纵横波相对反射系数和裂缝柔度参数.通过推导含裂缝柔度的方位各向异性反射系数公式,基于贝叶斯反演框架建立了P波线性AVAZ反演方法.合成地震数据应用表明基于贝叶斯理论的TTI介质裂缝柔度反演方法具有一定抗噪性,可以降低裂缝柔度估测的不确定性,为地下裂缝预测提供有力的依据.     

Seismic amplitude inversion for the transversely isotropic media with vertical axis of symmetry

Transverse isotropy with a vertical axis of symmetry is a common form of anisotropy in sedimentary basins, and it has a significant influence on the seismic amplitude variation with offset. Although exact solutions and approximations of the PP-wave reflection coefficient for the transversely isotropic media with vertical axis of symmetry have been explicitly studied, it is difficult to apply these equations to amplitude inversion, because more than three parameters need to be estimated, and such an inverse problem is highly ill-posed. In this paper, we propose a seismic amplitude inversion method for the transversely isotropic media with a vertical axis of symmetry based on a modified approximation of the reflection coefficient. This new approximation consists of only three model parameters: attribute A, the impedance (vertical phase velocity multiplied by bulk density); attribute B, shear modulus proportional to an anellipticity parameter (Thomsen's parameter ε−δ); and attribute C, the approximate horizontal P-wave phase velocity, which can be well estimated by using a Bayesian-framework-based inversion method. Using numerical tests we show that the derived approximation has similar accuracy to the existing linear approximation and much higher accuracy than isotropic approximations, especially at large angles of incidence and for strong anisotropy. The new inversion method is validated by using both synthetic data and field seismic data. We show that the inverted attributes are robust for shale-gas reservoir characterization: the shale formation can be discriminated from surrounding formations by using the crossplot of the attributes A and C, and then the gas-bearing shale can be identified through the combination of the attributes A and B. We then propose a rock-physics-based method and a stepwise-inversion-based method to estimate the P-wave anisotropy parameter (Thomsen's parameter ε). The latter is more suitable when subsurface media are strongly heterogeneous. The stepwise inversion produces a stable and accurate Thomsen's parameter ε, which is proved by using both synthetic and field data.