Analysis of the influence of P-S wave velocity variation on AVO type in the Dongfang area, Yinggehai Basin

  • Yong Deng ,
  • Shi-you Liu ,
  • Yang-sen Li ,
  • Chen Fu ,
  • Li-xia Shen
Expand
  • Zhanjiang Branch Company, China National Offshore Oil Corporation, Zhanjiang 524057, China

Received date: 2019-02-28

  Revised date: 2019-04-24

  Online published: 2019-10-14

Highlights

Because of the special temperature and pressure environment and formation mechanism in Dongfang area, low-speed mudstone develops especially, while low-speed mudstone and gas-bearing sandstone show relatively similar seismic facies and AVO characteristics, so the effect of identifying sandstone by AVO technique is not ideal. In order to effectively analyze the AVO characteristics of sandstone and mudstone in this area, this paper analyzes the influence of the relative variation of elastic parameters on AVO reflection coefficient based on the simplified formula of Zoeppritz equation. On this basis, the variation rule of gas-bearing sandstone of Huangliu Formation in Dongfang area with depth is clarified. It is clear that the shear wave velocity is the main factor affecting the response characteristics of III and IV types of gas-bearing sandstone in this area. Finally, the same IV type gas-bearing sandstone and low-speed mudstone are analyzed, and the IV type AVO anomaly of low-speed mudstone is more obvious. In order to distinguish the low-speed mudstone and gas-bearing sandstone effectively, the combination of petrophysics and prestack simultaneous inversion can effectively eliminate the impact of low-speed mudstone traps and provide important technical support for reservoir prediction and hydrocarbon detection in the Dongfang area.

Cite this article

Yong Deng , Shi-you Liu , Yang-sen Li , Chen Fu , Li-xia Shen . Analysis of the influence of P-S wave velocity variation on AVO type in the Dongfang area, Yinggehai Basin[J]. Natural Gas Geoscience, 2019 , 30(9) : 1378 -1385 . DOI: 10.11764/j.issn.1672-1926.2019.04.012

0 引言

莺歌海盆地位于南海北部大陆架西北部,为新生代强超压、转换伸展富含气盆地,包括莺东斜坡、莺西斜坡、中央坳陷等一级构造单元[1]。在中央坳陷中浅层天然气勘探利用亮点技术取得比较好的成效,东方1-1气田具有典型的“亮点”型地震响应,存在强振幅反射、频率衰减、低频阴影等地震响应[2]。但随后钻探东方X-6井出现意外情况,低频强振幅亮点反射特征并不是预期的砂岩储层而是低速泥岩[3]。受低速泥岩与含气砂岩反射特征相似的影响,迫使高温高压众多潜在目标钻探搁浅,如何评价这些目标的储层发育情况,避免低速泥岩陷阱成为当前勘探面临的紧迫问题。
相比叠后亮点技术,叠前AVO技术充分考虑地震振幅随偏移距的变化规律,提取的属性参数对流体信息更加敏感,可以有效减少油气检测的多解性,避免特殊岩性体导致的假亮点陷阱[4]。Ostrander[5]最早提出了利用反射系数随入射角变化识别亮点型含气砂岩,克服了传统亮点技术的缺陷。鉴于Zoeppritz方程难以直接用来分析各弹性参数对反射系数的影响,Aki等[6]假设界面两侧弹性性质相差较小条件下,给出了与速度和密度相对变化有关的反射系数近似公式;Shuey[7]推导了泊松比表示的近似式,指出反射系数受界面两侧泊松比之差影响,并提出了梯度和截距的概念;Smith等[8]提出用加权叠加方法估计流体因子和检测气层,实现了AVO对流体性质的检测;为了更好利用AVO属性来识别地层岩性和烃类检测,Rutherford等[9]建立了3类含气砂岩模型,并把AVO异常分为3类,Castagna等[10]后来又将其推广至4类AVO异常;Goodway等[11]利用拉梅常数、弹性模量推导的近似反射系数公式进行岩性识别和油气检测; Kelo等[12]在AVO分析中运用极化方法来识别异常体。在国内,郑晓东[13]提出了幂级数简化反射系数公式;李景叶等[14]从近似式出发,推导出利用多波时移的AVO反演油藏含油饱和度和压力变化的方程;印兴耀等[15]基于Zoeppritz方程近似式,提出直接从弹性波阻抗数据体中提取流体因子进行流体识别的方法。但随着面临的储层类型越来越复杂,AVO在含气砂岩识别上也存在特殊岩性的陷阱,因此亟需探索不同弹性参数对含气砂岩和特殊岩性的AVO特征的影响,以此避免特殊岩性体的勘探风险。
本文主要从Zoeppritz方程出发,根据Aki-Richards近似简化方程分析了纵横波速度及密度变化对AVO反射系数的贡献值影响,进而从本质上分析不同弹性参数变化对亮点型目标的III类、IV类AVO的变化规律。在此基础上结合东方区岩石物理统计规律及其拟合关系,建立东方区AVO的截距和梯度随埋深的变化规律模版,并发现横波速度是影响AVO类型的先决条件。最后对于相同AVO类型的低速泥岩和含气砂岩,分析了其AVO类型的变化特点,这里低速泥岩IV类AVO异常较含气砂岩更为明显,为有效区分低速泥岩和含气砂岩,结合岩石物理分析和叠前同时反演可以有效剔除低速泥岩陷阱的影响,为东方区储层预测和烃类检测提供重要技术支撑。

1 纵横波速度及密度变化对AVO的影响

当平面简谐波入射两层半无限弹性介质的分界面上时,产生反射纵波、透射纵波、转换反射纵波和转换透射横波,因此反射系数并不是简单的由界面上下介质的波阻抗关系求取,Zoeppritz方程很好地描述了弹性纵波在介质分界面上的反射和透射问题,则其方程具体表示如下:
s i n θ 1 c o s φ 1 - s i n θ 2 c o s φ 2 - c o s θ 1 s i n φ 1 - c o s θ 2 - s i n φ 2 s i n 2 θ 1 V P 1 V S 1 c o s 2 φ 1 ρ 2 V S 2 2 V P 1 ρ 1 V S 1 2 V P 2 s i n 2 θ 2 - ρ 2 V S 2 V P 1 ρ 1 V S 1 c o s 2 φ 2 c o s 2 φ 1 - V S 1 V P 1 s i n 2 φ 1 - ρ 2 V P 2 ρ 1 V P 1 c o s 2 φ 2 - ρ 2 V S 2 ρ 1 V P 1 s i n 2 φ 2 R p p R p s T p p T p s = - s i n θ 1 - c o s θ 1 s i n 2 θ 1 - c o s 2 φ 1
式中:Vp 1Vs 1ρ 1Vp 2Vs 2ρ 2分别为界面上下纵波速度,m/s、横波速度,m/s、密度,kg/m3θ 1θ 2φ 1φ 2分别为纵波入射、透射角、横波反射、透射角;Rpp、Rps为纵横波反射系数;Tpp、Tps为纵横波透射系数。
Aki-Richards假设相邻地层弹性参数变化较小条件,即ΔV pV p、ΔV s /Vs、Δρ/ρ和其他值相比为小值,所以可略去它们的高次项,将Zoeppritz方程进一步简化,得到了纵波反射系数近似方程:
R θ R 0 + G s i n 2 θ + F t a n 2 θ - s i n 2 θ
这里:
R 0 = 1 2 Δ V p V p + Δ ρ ρ
G = 1 2 Δ V p V p - 2 V s V p 2 2 Δ V s V s + Δ ρ ρ
F = 1 2 Δ V p V p
式中: Δ V p = V p 2 - V p 1为反射界面纵波速度变化量; V p = V p 2 + V p 1 / 2为反射截面纵波速度平均值; Δ ρ = ρ 2 - ρ 1表示反射界面密度的变化量; ρ = ρ 2 + ρ 1 / 2为反射界面密度平均值; Δ V s = V s 2 - V s 1为反射界面横波速度变化量; Δ V s = V s 2 + V s 1 / 2为反射界面横波速度平均值。
对密度、速度变化量分别进行合并同类项,可以得到纵波反射系数随密度、速度的变化率相对关系,即:
R θ 1 2 1 - 4 V s V p 2 s i n 2 θ Δ ρ ρ + 1 2 1 c o s 2 θ Δ V p V p - 4 V s V p 2 s i n 2 θ Δ V s V s
其中:密度项系数表示为 1 2 1 - 4 V s V p 2 s i n 2 θ,纵波速度项系数表示为 1 2 1 c o s 2 θ,横波速度项系数表示为 - 4 V s V p 2 s i n 2 θ
由上式可以看出,AVO反射系数主要可分成密度、纵波速度和横波速度变化率贡献,其弹性参数相对变化关系对比如图1所示。通常砂岩密度小于泥岩密度,因此对于密度项系数始终为负值(角度0o~45o范围),并随着入射角的增大,密度项系数绝对值越来越小,在大角度影响作用大。其次对于纵波速度项系数取决于上下地层纵波速度相对关系,当上覆地层纵波速度大于下伏地层纵波速度时,其纵波速度项系数为负值,随角度的增加纵波速度项系数绝对值越来越大;反之,纵波速度项系数为正值,随角度增加其贡献数值越来越大,纵波速度在法向入射角影响较大。而横波速度对AVO反射系数贡献度也取决于上下地层横波速度相对关系,其横波速度项系数同纵波速度规律正好相反,主要对AVO中远角度影响较大。因此AVO反射系数是密度、纵波速度和横波速度相对变化的综合响应结果,只是对于不同角度范围,各个弹性参数的贡献程度各不相同。
图1 不同弹性参数对AVO反射系数的贡献度对比分析

Fig.1 Comparative of the contribution of different parameters to the reflection coefficient of AVO

通常在中浅层含气砂岩AVO多以III类和IV类亮点型储层为主,含气砂岩往往表现为低纵波速度、低密度的岩石物理特征,因此AVO反射系数的贡献度均属于纵波速度和密度小于上覆地层的情况,密度项贡献度为负值,随着角度增加绝对值逐渐减小,纵波速度项贡献度为负值,随着角度增加其绝对值逐渐增加,这时AVO类型就取决于上下地层横波速度的相对关系(图2)。当上覆地层横波速度小于下伏地层横波速度时,横波速度变化率对反射系数贡献为负值,随着角度的增加,其绝对值越来越大,AVO整体表现为III类AVO现象;当上覆地层横波速度大于下伏地层横波速度时,横波速度变化率对反射系数贡献为正值,随着角度的增加,其值越来越大,AVO整体则表现为IV类AVO现象。
图2 III类、IV类AVO反射系数的贡献度对比分析

Fig.2 Comparative of the contribution of AVO reflection coefficients of class III and IV

2 南海某气田实例研究

在分析不同弹性参数相对变化对AVO反射系数影响的基础上,本文以东方区为例对该区上中新统黄流组含气砂岩AVO变化规律进行分析。东方区上中新统黄流组含气砂岩相对泥岩表现为低速、低密、低阻抗的岩石物理特征(图3),其中图中深度范围均属黄流组。在黄流组浅层,砂岩纵波速度和密度相对泥岩差异较大,气水差异相对明显;随着埋深增加,砂岩纵波速度和密度与泥岩差异减小,气水差异相对略小;而横波速度随深度变化较大,在黄流组浅层,砂岩横波速度大于泥岩横波速度,随着埋深的增加,砂岩横波速度小于泥岩横波速度,由于横波速度不受流体影响,气水差异整体较小。
图3 东方区黄流组密度、纵波速度及横波速度随深度变化关系

Fig. 3 The relationship of density,p-wave and s-wave velocity with depth in Huangliu Formation of Dongfang area

通过岩石物理统计分析,该区纵波速度和密度对AVO反射系数的贡献只有一种情况:在纵波速度贡献度方面,上覆盖层纵波速度大于下伏地层纵波速度,纵波速度对AVO反射系数贡献值为负值,并且随着角度增加,其绝对值逐渐增加;在密度贡献度方面,上覆盖层密度大于下伏地层密度,密度对AVO反射系数贡献值也为负值,并且随着角度增加,其绝对值逐渐减小。而该区横波速度对AVO反射系数的贡献有2种情况:当上覆盖层横波速度小于下伏地层速度,横波速度对AVO反射系数贡献值为负值,并且随着角度增加,其绝对值逐渐增加,这时含气砂岩AVO往往表现为III类AVO异常;当上覆地层横波速度大于下伏地层速度时,横波速度对AVO反射系数贡献值为正值,并且随着角度的增加,其贡献值逐渐增加,这时含气砂岩AVO表现为IV类AVO异常(图4)。因此,横波速度是影响该区含气砂岩III类、IV类AVO响应特征的主要因素。
图4 东方区AVO截距、梯度随深度变化关系

Fig. 4 The relationship between AVO intercept, gradient and depth in Dongfang area

此外,东方区上中新统黄流组是由于快速沉降和快速沉积形成了欠压实低速泥岩(X-6井),低速泥岩同含气砂岩(X-4井)相似,也表现为低纵波速度、低阻抗特征(图5)。从含气砂岩和低速泥岩AVO类型上看,由于含气砂岩和低速泥岩的横波速度均低于上覆盖层的横波速度,因此含气砂岩和低速泥岩都表现为IV类AVO特征(图6),并且低速泥岩横波速度相对盖层横波速度更低,IV类AVO的梯度更大,异常特征更明显。
图5 东方区X-4井,X-6井测井曲线

Fig.5 Well logging curve of Wells X-4 and X-6 in Dongfang area

图6 含气砂岩和低速泥岩AVO特征分析

Fig. 6 Analysis of AVO characteristics of gas-bearing sandstone and low-speed mudstone

为了有效区分低速泥岩和含气砂岩的IV类AVO异常,首先通过该区岩石物理分析发现,低速泥岩和含气砂岩都表现为低阻抗特征,在叠后地震剖面上都表现为强振幅波谷特征,这就给该区利用亮点技术发现油气藏带来一定困难[图7(a)]。但相对于含气砂岩,低速泥岩表现为高纵横波速度比特征,含气砂岩表现为低纵横波速度比特征[图7(b)],这就为利用叠前同时反演区分低速泥岩和含气砂岩提供了有效的解决途径。通过X-4井—X-6井叠前同时反演的纵横波速度比发现(图8),X-4井目的层表现为低纵横波速度比异常特征,而X-6井目的层段表现为高纵横波速度比异常特征,为低速泥岩概率的可能性较大,因此对于东方区低速泥岩陷阱,借助叠前反演的纵横波速度比可以有效降低该区勘探风险,为该区的储层预测和烃类检测提供重要技术支撑。
图7 不同流体砂岩和低速泥岩岩石物理交会分析

Fig. 7 Rocks physical analysis of sandstone and low-velocity mudstone

图8 东方区X-4—X-6井连井地震剖面及叠前反演结果

Fig.8 Seismic profile and prestack inversion results of Wells X-4-X-6 in DongFang area

3 结论

(1) 首先根据Aki-Richards近似公式,分析了不同纵横波速度及密度相对变化对AVO反射系数的影响,密度在大角度影响作用大,纵波速度在法向入射角影响较大,而横波速度对中远角度影响较大。
(2) 通过对东方区岩石物理趋势拟合,明确了东方区黄流组含气砂岩随深度的变化规律,其中横波速度是影响该区含气砂岩III类、IV类AVO响应特征的主要因素,若盖层为硬质泥岩,则含气砂岩表现为IV类AVO特征;相反则表现为III类AVO特征。
(3) 对于相同类型的低速泥岩和含气砂岩,可以借助叠前同时反演的纵横波速度比进行区分,通过反演证实该技术方法的有效性,为东方区储层预测和烃类检测提供重要技术支撑。
1
Deng Yong, Li Yangsen, Liu Shiyou,et al. Formation and identification of dark-spot gas reservoirs in Yinggehai Basin[J]. Progress in Geophysics, 201833(6): 2535-2540.

邓勇, 李洋森, 刘仕友,等. 莺歌海盆地乐东区暗点型气藏形成原因及识别方法[J]. 地球物理学进展, 201833(6): 2535-2540.

2
Tong Chuanxin, Wang Zhenfeng, Li Xushen.Pooling conditions of gas reservoirs in the Dongfang 1-1 Gas Field, Yinggehai Basin[J]. Natural Gas Industry, 2012, 32(8): 11-15.

童传新, 王振峰, 李绪深. 莺歌海盆地东方1-1气田成藏条件及其启示[J]. 天然气工业, 2012, 32(8):11-15.

3
Pei Jianxiang, Pan Guangchao, Zhu Peiyuan, et al. Identification of low velocity mudstone in the mid-deep overpressure zone,Yinggehai Basin[J].Oil Geophysical Prospecting,2016, 51(2): 361-370.

裴健翔, 潘光超, 朱沛苑,等. 莺歌海盆地超压带低速泥岩陷阱的甄别及对勘探的启示[J]. 石油地球物理勘探, 2016, 51(2): 361-370.

4
Yao Shufan, Ding Wenlong, Zhao Gang, et al. The relationship between diagenesis and AVO anomaly of gas-bearing sandstone in the ZJ Formation of the Pearl River Mouth Basin, China[J]. Geophysical Prospecting for Petroleum, 2018, 57(3): 443-451.

姚淑凡, 丁文龙, 赵刚,等. 珠江口盆地深水区珠江组成岩作用对含气砂岩AVO特征的影响研究[J]. 石油物探,2018, 57(3): 443-451.

5
Ostrander W J. Plane-wave reflection coefficients for gas sands at nonnormal angles of incidence[J]. Geophysics,1984,49(10): 1637-1648.

6
Aik K, Richards P.Quantitative seismology:Theory and me-thods[J]. San Francisco : Freeman19801:557 .

7
Shuey R T.A simplification of the Zoeppritz equations[J].Geo-physics198550(9): 609-614.

8
Smith G C, Gidlow P M. Weighted stacking for rock property estimation and detection of gas[J]. Geophysical Prospecting, 1987, 35(9): 993-1014.

9
Rutherford S R, Williams R H. Amplitude-versus-offset variations in gas sands[J]. Geophysics, 1989, 54(6): 680-688.

10
Castagna J P, Swan H W, Foster D J. Framework for AVO gradient and intercept interpretation[J]. Geophysics,1998, 63(3): 948-956.

11
Goodway W, Chen T, Downton J. Improved AVO fluid detection and lithology discrimination using Lame parameters; λρ, μρ, & λ/μ fluid stack, from P and S inversions[J]. SEG Technical Program Expanded Abstracts, 1997,16(1):2067.

12
Keho T, Lemanski S, Ripple R, et al. The AVO hodogram: Using polarization to identify anomalies[J]. The Leading Edge, 2001,20(11):1214-1224.

13
Zheng Xiaodong.Some new development of AVO techni-que[J]. Oil Geophysical Prospecting, 1992, 27(3):307-317.

郑晓东.AVO理论和方法的一些新进展[J].石油地球物理勘探1992, 27(3):307-317.

14
Li Jingye, Chen Xiaohong. The method for impedance inversion of time-lapse seismic data[J].Progress in Exploration Geophysics, 2003, 26(1): 41-44.

李景叶,陈小宏. 时移地震波阻抗反演方法研究[J]. 勘探地球物理进展, 2003, 26(1): 41-44.

15
Yin Xingyao, Zhang Shixin, Zhang Fanchang, et al. Utilizing russell approximation based elastic wave impedance inversion to conduct reservoir description and fluid identification[J]. Oil Geophysical Prospecting, 2010, 45(3):373-380.

印兴耀, 张世鑫, 张繁昌,等. 利用基于Russell近似的弹性波阻抗反演进行储层描述和流体识别[J]. 石油地球物理勘探, 2010, 45(3):373-380.

Outlines

/