部分饱和砂岩储层地震物理模拟及含气饱和度预测分析

含气饱和度预测是天然气储层地震解释工作的重要目标.本文将岩石物理分析与地震物理模拟技术相结合,构建了部分;饱和砂岩储层物理模型并进行含气饱和度预测分析.物理模型中设置了高孔渗常规砂岩和低孑孔渗致密砂岩两种模拟储层,每种储层都是由具有不同含水饱和度的气-水双相饱和砂体组成.岩石物理分析结果显示在低孔渗致密砂岩中气-水混合流体更加倾向于非均匀的斑块分布,而结合了Brie等效流体公式的Gassmann流体替换理论可以更准确地描述纵波速度随含水饱和度的变化趋势.对物理模型进行地震资料采集处理后,对比了AVO特征和叠前同步反演结果对两种砂岩储层含气饱和度预测能力的差异.AVO特征结果显示,对于混合流体均匀分布的高孔渗砂岩储层,AVO响应曲线和属性变化很难对含气饱和度进行估算;对于混合流体斑块分布的致密砂岩储层,AVO特征可以定性地分辨出储层是否为高、中、低含气情况.反演结果显示,密度及纵横波速度比分别对高孔渗及致密砂岩储层的含气饱和度有着较好的指示能力.     

中深层天然气藏地震预测技术研究

济阳坳陷天然气资源丰富,已在多个层系发现了不同类型的中深层天然气藏,由于气藏地质条件、地震反射复杂,针对不同的气藏需采用不同的地震预测技术。本文在中深层天然气藏地质特点分析的基础上进行了岩石物理、地震反演、地震属性等预测技术研究,结果表明：中层气藏储层物性较好,储层与围岩的岩石物理参数差异明显,应用亮点、AVO、叠前弹性阻抗反演技术可较好地识别预测三角洲砂岩气藏;针对深层低孔低渗砂岩气藏,地震属性及多属性反演是有效的储层预测技术,并预测了孤北地区二叠系有利储集相带分布,结合构造分析,指出了孤北地区深层气勘探方向。     

Numerical simulation of frequency-dependent seismic response and gas reservoir delineation in turbidites: A case study from China

To simulate the frequency-dependent response of turbidite reservoirs in the JZ Area, near the Bohai Sea, China, we apply a diffusive and viscous wave equation (DVWE), which takes into account diffusive and viscous attenuation and velocity dispersion in fluid-bearing poroelastic media. We use a seismic data-driven geological model building approach to produce physical parameter sections, which are then used to numerically synthesize the frequency-dependent seismic response in the DVWE-based simulation. The DVWE-based synthetic section shows the characteristic reflection and geometry of turbidites and delineates the phase delay, instantaneous dominant frequency decrease and magnitude attenuation related to the gas-bearing reservoir. The common frequency sections obtained by instantaneous spectral decomposition of the synthetic section show that a low-frequency shadow (LFS) lies immediately beneath the reservoir. Next, following the implications of the numerical simulation, we then apply LFS and fluid mobility to the data volume. Both hydrocarbon indicators clearly delineate the bright gas reservoir and its spatial distribution. The workflow and methodologies can be expected to be applicable to other frequency-dependent hydrocarbon indicators.     

基于频变AVO反演的页岩储层含气性地震识别技术

结合地震岩石物理技术,研究了叠前频变AVO反演在四川盆地龙马溪组页岩储层含气性识别中的应用.首先,应用Backus平均理论将测井数据粗化为地震尺度储层模型,应用传播矩阵理论进行高精度地震正演及井震标定,分析页岩气储层地震响应特征.其次,基于岩心观测结果,应用Chapman多尺度裂缝理论设计页岩气储层理论模型,研究储层衰减、频散以及对应的地震反射特征.应用该理论模型测试频变AVO反演方法,计算结果表明：对于研究区地层结构和地震数据,区分流体类型的优势频率不是地震子波的主频,还受层间调谐干涉等储层结构因素控制,也进一步说明理论模型测试和标定的重要性.最后,将频变AVO反演技术应用到四川盆地龙马溪组页岩地层,计算得到的频散属性为页岩气储层含气性识别提供依据.     

A new energy density calculation method for reservoir delineation using seismic low-frequency signal

In comparison to high-frequency signals, low-frequency seismic signals suffer less from scattering and intrinsic attenuation during wave propagation, penetrate deeper strata and thus can provide more energy information related to the hydrocarbon reservoirs. Based on the asymptotic representation for the frequency-dependent reflections in the fluid-saturated pore-elastic media, we first derive a novel equation of the reservoir energy density and present an efficient workflow to calculate the reservoir energy density using low-frequency seismic data. Then, within a low-frequency range (from 1 to 30 Hz), we construct an objective function to determine the optimal frequency, using the energy densities calculated from the post-stack seismic traces close to the wells. Next, we can calculate the reservoir energy density using the instantaneous spectra of optimal frequency at the low-frequency end of the seismic spectrum. Tests on examples for synthetic and field data demonstrate that the proposed reservoir energy density can produce high-quality images for the fluid-saturated reservoirs, and it produces less background artefacts caused by elastic layers. This method provides a new way to detect the location of hydrocarbon reservoirs and characterize their spatial distribution.     

多波时移地震AVO反演研究

数值模拟了油藏含油饱和度与有效压力变化时移地震AVO的响应,确定利用时移地震AVO区分油藏参数的变化、实现油藏定量解释的可行性.从Aki等 AVO近似方程出发,详细推导了P_P波和P_S转换波时移地震AVO计算公式.结合岩石物理近似关系和本文推导的时移地震AVO计算公式,推导了利用多波时移地震AVO反演油藏含油饱和度和压力变化的方程.数据试验表明,文中推导的多波时移地震AVO方程能较好地反演油藏含油饱和度变化和有效压力变化,实现油藏定量解释.     

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.     

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.     

Study on attribute characterization for reservoir dynamic monitoring by seismic

Study on characterizing reservoir parameters dynamic variations by time-lapse seismic attributes is the theoretical basis for effectively distinguishing reservoir parameters variations and conducting time-lapse seismic interpretation,and it is also a key step for time-lapse seismic application in real oil fields. Based on the rock physical model of unconsolidated sandstone,the different effects of oil saturation and effective pressure variations on seismic P-wave and S-wave velocities are calculated and analyzed. Using numerical simulation on decoupled wave equations,the responses of seismic amplitude with different offsets to reservoir oil saturation variations are analyzed,pre-stack time-lapse seismic attributes differences for oil saturation and effective pressure variations of P-P wave and P-S converted wave are calculated,and time-lapse seismic AVO (Amplitude Versus Offset) response rules of P-P wave and P-S converted wave to effective pressure and oil saturation variations are compared. The theoretical modeling study shows that it is feasible to distinguish different reservoir parameters dynamic variations by pre-stack time-lapse seismic information,including pre-stack time-lapse seismic attributes and AVO information,which has great potential in improving time-lapse seismic interpreta-tion precision. It also shows that the time-lapse seismic response mechanism study on objective oil fields is especially important in establishing effective time-lapse seismic data process and interpreta-tion scheme.     

基于频变AVO技术对多尺度裂缝内流体属性反演与识别（英文）

通过地震数据获取裂缝储藏中流体的性质并对流体类型进行识别,是地震勘探岩性反演的重要问题之一。由于地震波的速度、储层的密度等弹性参数对某些流体不具有很强的敏感性,使只依赖振幅信息进行流体识别的传统AVO方法面临困境。作为传统叠前振幅反演的一个拓展,频变AVO（FDAVO）技术进一步考虑了振幅对频率的依赖关系,将这种依赖关系与地下裂缝结构、流体填充对应起来,能带来更丰富的流体信息。利用该技术,本文提出了一种基于地震数据参数化Chapman模型的贝叶斯反演新方法（BIDCMP）,它包含两步算法,即,FDAVO反演储层的非弹性属性和贝叶斯框架下的流体识别。首先,通过匹配观测数据和模型数据,构造差函数反演裂缝储层非弹性参数。随后,在贝叶斯框架下,使用马尔科夫随机场（MRF）作为先验模型,联合多参数场识别流体。本方法在计算过程中,除综合考虑了弹性参数场、测井资料等常规信息外,还特别地加人了第一步中反演得的非弹性参数的约束,从而充分利用了流体粘性差异,最后在最大后验概率（MAP）准则下输出最佳岩性一流体识别结果。分别对合成地震记录和模拟岩性—流体剖面验证本文方法的有效性,结果证明本文方法获得的流体识别结果准确可信。     

Frequency‐ and angle‐dependent poroelastic seismic analysis for highly attenuating reservoirs

We extend the frequency‐ and angle‐dependent poroelastic reflectivity to systematically analyse the characteristic of seismic waveforms for highly attenuating reservoir rocks. It is found that the mesoscopic fluid pressure diffusion can significantly affect the root‐mean‐square amplitude, frequency content, and phase signatures of seismic waveforms. We loosely group the seismic amplitude‐versus‐angle and ‐frequency characteristics into three classes under different geological circumstances: (i) for Class‐I amplitude‐versus‐angle and ‐frequency, which corresponds to well‐compacted reservoirs having Class‐I amplitude‐versus‐offset characteristic, the root‐mean‐square amplitude at near offset is boosted at high frequency, whereas seismic energy at far offset is concentrated at low frequency; (ii) for Class‐II amplitude‐versus‐angle and ‐frequency, which corresponds to moderately compacted reservoirs having Class‐II amplitude‐versus‐offset characteristic, the weak seismic amplitude might exhibit a phase‐reversal trend, hence distorting both the seismic waveform and energy distribution; (iii) for Class‐III amplitude‐versus‐angle and ‐frequency, which corresponds to unconsolidated reservoir having Class‐III amplitude‐versus‐offset characteristic, the mesoscopic fluid flow does not exercise an appreciable effect on the seismic waveforms, but there exists a non‐negligible amplitude decay compared with the elastic seismic responses based on the Zoeppritz equation.     

莺歌海盆地“暗点”型气藏地震识别技术研究

莺歌海盆地乐东区深层与东方区具有类似的成藏条件,是高温高压领域获取下一个重大勘探突破的最现实区域.但是该区已钻遇的“亮点”型储层物性差,普遍表现为低孔-特低渗特征,而区域基础规律研究表明,区域潜在甜点储层在地震上应该表现为“暗点”型特征.但是,“暗点”型储层在叠后地震资料上无法有效识别,如何准确预测“暗点”型储层是乐东区勘探新突破的重点攻关方向.本文针对“暗点”识别难的问题,提出了基于叠前地震资料的储层预测技术.从岩石物理出发系统性地梳理该区“暗点”型储层的地球物理特征,总结出其具有“叠加无反射、近道弱波峰、远道弱波谷、道集极性反转”的规律,根据该规律建立了Fstack属性重构算法,利用该算法开展甜点储层筛选并进行“暗点”储层预测,结合道集资料、反演结果等地球物理信息,最终确定潜在“暗点”型气层的分布范围.钻井结果证实该技术对乐东区深层潜在“暗点”型储层具有很好的预测效果.     

Seismic low-frequency-based calculation of reservoir fluid mobility and its applications

Low frequency content of seismic signals contains information related to the reservoir fluid mobility. Based on the asymptotic analysis theory of frequency-dependent reflectivity from a fluid-saturated poroelastic medium, we derive the computational implementation of reservoir fluid mobility and present the determination of optimal frequency in the implementation. We then calculate the reservoir fluid mobility using the optimal frequency instantaneous spectra at the low-frequency end of the seismic spectrum. The methodology is applied to synthetic seismic data from a permeable gas-bearing reservoir model and real land and marine seismic data. The results demonstrate that the fluid mobility shows excellent quality in imaging the gas reservoirs. It is feasible to detect the location and spatial distribution of gas reservoirs and reduce the non-uniqueness and uncertainty in fluid identification.     

有效压力与饱和度变化的气藏多波时移地震AVO响应

天然气在开发过程中,储层有效压力和含气饱和度均会发生变化,研究有效压力和含气饱和度的变化对地震响应特征的影响,在基于时移地震的剩余气分布预测研究中具有重要意义。天然气和石油的声学性质有着明显的差异,油藏时移地震的研究成果不能直接应用于气藏,因此需要开展气藏的时移地震研究。利用Shapiro模型表征干岩石弹性模量随有效压力的变化,借助Batzle-Wang方程描述流体速度随压力的变化关系,联合Gassmann理论进行流体替代,表征饱和流体岩石速度随含气饱和度的变化,建立了饱和流体岩石速度随有效压力和饱和度变化的岩石物理模型。基于该模型,对不同含气饱和度和不同有效压力下的气藏储层模型进行了多波时移地震叠前振幅变化（AVO）模拟。结果表明多波时移地震AVO技术可以有效地区分有效压力变化和含气饱和度变化,为进一步开展气藏多波时移地震流体监测提供了理论参考依据。      

球面波PP反射系数的频变特征研究

与平面波反射系数相比,球面波反射系数可以更精确地描述实际地震波的反射特征.近些年关于球面波的研究主要聚焦于球面波反射系数随入射角的变化规律,很少对球面波反射系数随频率的变化(频变)做详细研究.为了更全面地了解球面波的反射机制,本文研究了两层弹性介质中球面波PP反射系数(幅值和相位)的频变规律.文中基于经典的Sommerfeld积分构造球面波PP反射系数,通过自适应的Gauss-Kronrod求积算法对其进行稳定的数值计算.数值试验发现,对于不同的介质参数,球面波反射系数表现出了复杂的频变规律.尤其是当平面波反射系数为零时,对应的球面波反射系数是非零的,且球面波反射系数的相位随频率增加在高频趋近于90°或-90°,即此时球面反射波相对于入射波会有90°的相位旋转.对四类AVO模型的测试表明,球面波反射系数与平面波反射系数在临界角附近和低频时差异很大.     

Analysis of reservoir heterogeneity-induced amplification effect on time-lapse seismic responses of fluid substitution: A physical modelling study

Seismic wave propagation through a fluid-saturated poroelastic layer might be strongly affected by media heterogeneities. Via incorporating controlled laboratory simulation experiments, we extend previous studies of time-lapse seismic effects to evaluate the wave scattering influence of the heterogeneous nature of porous permeable media and the associated amplification effects on 4D seismic response characteristics of reservoir fluid substitution. A physical model consisted of stratified thin layers of shale and porous sandstone reservoir with rock heterogeneities was built based on the geological data of a real hydrocarbon-saturated reservoir in Northeast China. Multi-surveys data of good quality were acquired by filling poroelastic reservoir layers with gas, water and oil in sequence. Experimental observations show that reservoir heterogeneity effect causes significantly magnified abnormal responses to the fluid-saturated media. Specifically, reflection signatures of the gas-filled reservoir are dramatically deviated from those of the liquid fluid-filled reservoir, compared with ones of the homogeneous media. By removing the influences unrelated to reservoir property alterations, 4D seismic estimates of travel-time and frequency-dependent characteristic are reasonably consistent with fluid variations. Nevertheless, strong 4D amplitude difference anomalies might not correspond to the regions where fluid variations occur. We also find that 4D seismic difference attributes are evident between oil- and water-filled models, whereas significant between oil- and gas-filled models. Meanwhile, rock physics modelling results reveal the predicted 4D seismic differences are obviously smaller than those calculated from seismic observations. The results in this paper, therefore, implicate that the effect of a reservoir's heterogeneous nature might be beneficial for hydrocarbons detection as well as monitoring small variations in pore fluids.     

基于流体替换技术的地震AVO属性气藏识别(英文)

传统上,油藏地球物理工程师是基于测井数据进行流体替换,计算油藏饱和不同流体时的弹性参数,并通过地震正演模拟分析油藏饱和不同流体时的地震响应,从而进行油气藏识别研究。该研究方案为油藏研究提供了重要的弹性参数和地震响应信息,但这些信息仅限于井眼位置。对于实际油藏条件,地下储层参数都是随位置变化而变化的,如孔隙度、泥质含量和油藏厚度等,因此基于传统流体替换方案得到的流体变化地震响应信息对于油气藏识别具有很大的局限性。研究通过设定联系油藏弹性参数与孔隙度、矿物组分等参数的岩石物理模型,并基于三层地质模型,进行地震正演模拟与AVO属性计算。得到油藏孔隙度、泥质含量和储层厚度变化时地震AVO属性,并建立了饱和水储层和含气储层对应AVO属性(包括梯度与截距)之间的定量关系。建立的AVO属性之间的线性关系可以实现基于地震AVO属性直接进行流体替换。最后,应用建立的流体替换前后AVO属性之间线性方程,对模拟地震数据直接进行流体替换,并通过流体替换前后AVO属性交汇图分析实现了气藏识别。     

Petrophysical evaluation and its application to AVO based on conventional and CMR-MDT logs

常规测井为 AVO 分析提供了基础的资料, 成为联系岩石物理与地震资料的桥梁。然而如果储层存在有复杂的流体系统, 如地层被严重地层流体侵入、电阻率响应低及盐水矿化度复杂等的现象, 则常规测井无法提供高质量的测井资料, 导致得出错误的弹性计算结果, 使 AVO 结果与地震资料不吻合。中国渤海湾地区第三系裂缝性储层复杂, 我们利用常规测井和核磁共振测井与模块地层动态测试相结合的组合仪完成了地层评价和储层描述。研究结果表明岩石物理学家利用上述方法技术可以获得诸如空隙度、渗透率、含水饱和度、束缚流体以及空隙压力等重要的储层参数并进一步综合应用这些结果和以实验室测量数据为基础的岩性分析结果进行在地震域岩石物理研究和 AVO 分析。     

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.     

基于AVO反演的频变流体识别方法

研究表明流体引起衰减与频散往往表现为频变AVO现象.一些频散地震属性,例如纵波频散,已经证实为可靠的碳氢指示因子.为了更有效地识别流体,基于f-μ-ρ近似构建了新的流体因子D_f,即频变流体项.该属性的反演首先需要连续小波变换(CWT)谱分解得到不同频带地震数据,通过去相关与先验约束来保证反演结果可靠性.模型试算证实了频变反射系数近似公式的精度可靠性,D_f可以识别出强衰减介质所引起的频散现象.实际数据试算中,D_f可以较好地识别储层孔隙流体,尤其对于气层,具有较好的指示效果.该流体因子将Gassmann流体项的高孔隙流体敏感性与叠前数据丰富的振幅频率信息相结合,反演效果与岩石物理认识相符.此研究有助于利用衰减频散现象借助AVO反演实现流体识别.