油藏模型含油饱和度变化的地震振幅响应特征
吴宝年1,2,3
1.中国石化 石油勘探开发研究院,北京 100083
2.中国石化集团国际石油勘探开发有限公司,北京 100083
3.中国石油大学(北京) 地球科学学院,北京 102249

作者简介: 吴宝年(1976-),男,高级工程师,博士后在读,研究方向为含油气盆地分析及物探解释技术。

摘要

针对储层流体地震检测技术不容易识别油层(不含气)的现状,深入研究油、水弱弹性参数差异条件下的地震响应同油水饱和度变化之间的关系,对提升油层地震检测成功率有重要意义。因此,利用矢量波场数值模拟技术,对不同含油饱和度条件下的油藏模型进行波场正演,利用正演结果探索性地研究仅含油—水混合流体油藏的含油饱和度变化引起的地震响应特征,为油层地震检测和定量解释提供参考。研究结果显示:仅含油—水混合流体的油藏在不同含油饱和度时的AVO响应特征基本一致,利用常规手段较难有效区分含油与含水,但通过压制储层反射中的非油气响应,突出油气响应的振幅校正处理后,可以明显看到储层含油饱和度变化引起的地震响应,且不同含油饱和度范围的地震响应特征不同,呈明显的非线性特征,因此,地震资料携带了储层流体信息,利用地震资料检测储层含油性及含油量是可行的,但对储层含油性和饱和度的合理解释需要以正确认识储层含油饱和度变化同地震响应之间的非线性规律为前提基础。

关键词: 含油饱和度; 地震检测; 振幅校正; 振幅响应
中图分类号:P631.4 文献标志码:A 文章编号:1000-8918(2015)06-1271-07
Seismic amplitude response characteristics of oil-bearing saturability variation of the oil pool model
WU Bao-Nian1,2,3
1.Exploration and Production Research Institute,SINOPEC,Beijing 100083,China
2.SINOPEC International Petroleum Exploration and Production Corporation,Beijing 100083,China
3.College of Geosciences,China University of Petroleum,Beijing 102249,China
Abstract

To tackle the status that the reservoir fluid seismic detection technique cannot easily recognize oil bed (containing no gas),it is of great significance to study the relationship between the seismic response and the oil-water saturability under the condition of the difference in oil and water weak elastic parameters so as to improve the success rate of the oil bed seismic detection.Utilizing the vector wave field digital simulation technique,the author performed wave field forward processing of the oil pool under the condition of different oil-bearing saturabilities and,using the inversion results,tentatively studied the seismic response characteristics caused by the oil-bearing saturability variation in the oil pool that contains only oil-water mixed fluid,so as to provide reference for the seismic detection and quantitative interpretation of the oil bed.According to the results obtained,the AVO response characteristics of the oil reservoir that contains only oil-water mixed fluid are basically consistent under the condition of different oil-bearing saturabilities,and hence it is difficult to use routine means to distinguish oil from water.Nevertheless,through suppressing non-oil and gas response in the reflection of the oil reservoir and making prominent the oil and gas response by amplitude correction,the seismic response resulting from the oil-bearing saturability of the oil reservoir can be clearly observed,and the seismic response characteristics in different oil-bearing saturability ranges are different and exhibit evident nonlinear features.In this way,the seismic data carry fluid information of the oil reservoir,and the utilization of the seismic data to detect the oil-bearing potential and oil-bearing quantity is feasible.However,the reasonable interpretation of the oil-bearing potential and saturability should be based on the prerequisite of the correct understanding of the nonlinear relationship between the oil-bearing saturability variation of the oil reservoir and the seismic response.

Keyword: oil-bearing saturability; seismic detection; amplitude correction; amplitude response

利用地震勘探进行油气预测的方法、技术众多[1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13], 但在油气勘探实践中, 这些技术方法成功应用的实例绝大多数是对气层的检测[13, 14], 而对油层的识别效果则较差。其主要原因是油和水的弹性参数(如体积模量和密度)差别很小, 不能引起足够大且易识别的地震反射差异, 尤其是当储层中油水混合流体的饱和度变化时, 更不易区分其中的差异。因此, 针对储层流体地震检测技术不容易识别油层(不含气)的现状, 深入研究弱弹性参数差异条件下的地震响应同油水饱和度变化之间的关系, 对提升油层地震检测成功率有重要意义。

笔者利用地震波场数值模拟技术, 使用矢量波动方程对利用Boit-Gassmann方程建立的油藏地质模型在含油饱和度变化时进行地震波场正演, 根据正演结果分析油藏储层在不同含油饱和度时的反射P-P波和P-S波振幅响应, 研究和探讨储层含油饱和度变化引起的P-P波和P-S波振幅振变化特征, 为油层地震检测和定量解释提供参考。

1 油藏岩石物理模型建立

为便于分析, 建立简化的双层油藏地质模型, 上层为泥岩盖层, 下层为砂岩储层。泥岩盖层纵、横波速度和密度保持不变, 分别为2 320 m/s, 1 031 m/s和2.15 g/cm3。砂岩储层为纯石英砂岩, 孔隙度30%, 顶面埋深600 m, 含油饱和度范围为0~1.0。储层砂岩基质模量为石英的体积模量和剪切模量, 干岩石骨架模量由Murphy等研究干岩石骨架模量与孔隙度的关系给出的最佳经验拟合公式计算[15], 孔隙流体体积模量用Wood方程计算[16]。储层在不同含油饱和度时的体积模量和剪切模量用Boit-Gassmann方程计算, 纵、横波速度和密度分别由

vp=[K+(4μ)/3]/ρ, vs=μ/ρρ=ρd+ϕρf, ρf=ρwSw+ρhc(1-Sw)1

计算。式中:vpvs分别为纵、横波速度; ρ Kμ 分别为流体饱和多孔岩石的密度、体积模量和剪切模量; ρ dρ f分别是干燥岩石、孔隙流体密度; ρ wρ hc分别是水、烃类密度; Sw为含水饱和度; ϕ 为孔隙度。

图1 模型砂岩速度随含油饱和度So的变化关系

计算出的砂岩储层在不同含油饱和度时的纵、横波速度和密度见图1图2。从图中可看到模型储层随含油饱和度增加, 纵波速度呈下降趋势, 横波速度则略微增加, 密度呈缓慢下降趋势。

图2 模型砂岩密度随含油饱和度So的变化关系

2 不同含油饱和度油藏模型地震波场正演

鉴于矢量波动方程正演能很好地保持地震波场传播过程中的动力学信息, 因此使用矢量波动方程对建立的油藏模型进行波场正演, 含油饱和度从0~1, 步长变化为0.1, 计算方法使用伪谱法[17]。伪谱法没有数值频散, 计算精度高, 能达到有限差分法的极限值[18]等优点。波场计算使用的波动方程为弹性介质运动平衡方程的位移方程, 其矢量形式

(λ+μ)(·U)+μ2U+ρF=ρ2Ut2,  (2)

式中:λ μ 为拉梅系数, ∇2= 2x2+ 2y2+ 2z2, U为位移矢量, ρ 为密度, F为体力, t为时间。

由式(2)可得到以速度和位移形式表示的运动平衡方程

2ut2=vp22ux2+2vxy+2wxz+vs22uy2+2uz2-2vxy-2wxz+Fx2vt2=vp22uxy+2vy2+2wyz+vs22vx2+2vz2-2uxy-2wyz+Fy2wt2=vp22uxz+2vyz+2wz2+vs22wx2+2wy2-2uxz-2vyz+Fz3

式中:uvw分别为xyz方向的位移分量; FxFyFz分别为xyz方向的体力分量。

求解方程组(3)得到的是纵、横波混合的矢量波场, 为便于分析, 需要从混合波场中分离出纵波场和横波场。对混合波场分别求散度和旋度, 可得到分解的纵波场和横波场。由于弹性波场可以用一个标量位的梯度和一个矢量位的旋度之和表示。设位移矢量u的标量和矢量位函数分别为φ ψ ; 力位函数F的标量和矢量位函数分别为Φ Ψ 。于是位移矢量u和力位函数F分别可表示为

u=up+us=φ+×ψ, (4)F=Fp+Fs=Φ+×Ψ(5)

式中:upFpusFs分别满足∇× up=∇× (∇φ )=0, ∇× Fp=∇× (∇Φ )=0, ∇· us=∇· (∇× ψ )=0, ∇· Fs=∇· (∇× Ψ )=0。

将式(4)、式(5)代入式(3)求解, 可同时得到混合波场、分解的纵波场和横波场。

波场计算使用主频50 Hz的雷克子波作震源函数[17], 边界条件采用吸收边界, 道间距20 m, 时间采样间隔1 ms。图3为含油饱和度0.80(油藏模型So=0.80), 炮点坐标(640 m, 40 m), 记录长度1.8 s的单炮正演记录。

图3 油藏模型So=0.8时的单炮地震记录(a— 混合波水平分量; b— 混合波垂直分量; c— P波水平分量; d— P波垂直分量; e— S波水平分量; f— S波垂直分量)

3 储层含油饱和度变化振幅响应特征
3.1 P-P波AVO特征分析及振幅校正

为便于说明, 将油藏模型含油饱和度从0~1的过程简化为含油饱和度分别为0、0.3、0.6、0.9代表的纯水层、含油水层、含水油层、油层四种状态进行分析说明。

从振幅随偏移距的变化关系(如图4)可看到以上四种状态的储层P-P波垂直分量振幅随偏移距增大均呈减小趋势, 且变化特征一致, 不容易区分。

图4 含油饱和度不同时模型储层反射P-P波原始振幅极大值随偏移距的变化关系

由于储层地震反射振幅是岩性、物性、含流体性等多因素的综合响应, 储层中油水混合流体的地震响应相对岩性等其他因素引起的地震响应弱很多, 直接从储层地震反射中分析油水饱和度变化引起的振幅变化, 会很困难, 而且多解性强。因此, 需要压制储层非油气因素引起的地震响应, 突出油气因素引起的地震响应, 并在此基础上进行储层含油饱和度变化引起的地震振幅响应特征分析。

能引起储层反射振幅变化的因素很多, 除含油气外, 地震波入射角、岩性、物性、含水、波前扩散等都可以引起振幅的变化。为消除或压制这些非油气因素引起的振幅变化, 以储层含油饱和度为0代表的纯水层为基准, 计算振幅校正系数, 计算公式

ri=A0/Ai, (6)

式中:ri为第i道地震波振幅校正系数, A0为零偏移距道的地震波振幅, Ai为第i道地震波振幅。

然后, 用振幅校正系数与储层反射记录相乘, 得到压制储层非油气因素振幅响应之后的地震记录。此时, 地震记录中储层反射振幅随偏移距的变化基本反映储层中的烃类引起的振幅变化。图5为含油饱和度分别为0、0.3、0.6、0.9时储层反射P-P波垂直分量振幅校正前后对比, 图中含油饱和度为0的储层反射记录经振幅校正后无AVO现象。

图5 含油饱和度不同时的振幅校正前、后P-P波垂直分量振幅记录

3.2 含油饱和度变化P-P波AVO特征分析

图6为含油饱和度为0, 0.3, 0.6, 0.9时的P-P波作压制储层非油气因素振幅响应处理后的垂直分量振幅随偏移距的变化关系。从图中可看到振幅校正后的P-P波垂直分量振幅随偏移距的变化在不同含油饱和度具有明显不同的特征。当含油饱和度不为零时, P-P波垂直分量振幅随偏移距增大而减小, 并随含油饱和度增加, 振幅衰减幅度加大, 且不同含油饱和度对应的AVO曲线形态不同, 说明含油因素引起的振幅变化同含油饱和度的变化是非线性关系, 如果利用线性刻度进行基于振幅的储层含油饱和度解释, 则会有一定的风险。但根据储层反射压制非油气响应之后保留的烃类响应的AVO特征可以判断油层的含油饱和度范围。另外, 不同偏移距范围内的烃类AVO响应特征不同, 如含水油层(So=0.6), 在近— 中远偏移距范围内无AVO现象, 远偏移距之后振幅快速衰减。

图6 含油饱和度不同时模型储层反射振幅校正后P-P波垂直分量振幅随偏移距变化关系

分析振幅校正后的储层反射记录偏移距分别为0 m(地震波入射角为0)、500 m(地震波入射角为23° )、1 000 m(地震波入射角为40° )代表的近、中、远偏移距振幅同含油饱和度之间的变化关系, 如图7所示。从图中可看到, 振幅校正后的储层反射记录近、中偏移距P-P波垂直分量振幅随含油饱和度的增加呈非线性变化, 二者变化规律一致, 而远偏移距P-P波垂直分量振幅随含油饱和度的增加呈近线性减小。在近— 中偏移距道, 含油引起的P-P波垂直分量振幅随含油饱和度增加, 当0≤ So< 0.4时呈线性减小; 0.4≤ So< 0.7时, 为常数, 不随含油饱和度的变化而变化; 当So> 0.7时, 也呈线性减小, 但斜率小于So< 0.4时的线性下降斜率, 且二者相交于So=0.5的位置。对于远偏移距道, 含油引起的P-P波垂直分量振幅随含油饱和度增加, 呈近线性减小。利用这一现象, 可以根据时移地震近、中、远偏移距道集中的油藏反射振幅变化判断油藏含油饱和度。

图7 不同偏移距道的振幅校正后P-P波振幅垂直分量随So的变化关系

3.3 含油饱和度变化P-S波AVO特征分析

对正演记录中的P-S波作振幅校正, 图8为含油饱和度分别为0, 0.3, 0.6, 0.9时的P-S波振幅校正后垂直分量振幅随偏移距的变化关系。从图中可看到振幅校正后的储层反射P-S波振幅不随含油饱和度的变化而变化, 且无AVO现象。其原因是:振幅校正后的储层地震反射记录已消除非油气因素引起的储层反射振幅变化, 反映的是储层中烃类流体信息, 而横波不能在流体中传播, 因此不会出现流体因素引起的横波振幅变化。

图8 含油饱和度不同时的P-S波振幅校正后垂直分量随偏移距的变化关系

4 结论与认识

通过分析和对比不同含油饱和度的油藏地质模型正演合成地震记录中的储层反射振幅特征, 可以看到, 油层中油水混合流体引起的地震响应相对岩性、地震波入射角改变等其他因素引起的地震响应弱很多, 直接根据储层反射振幅的变化区分储层含油与含水及定量解释饱和度难度较大, 多解性强。研究通过压制储层非油气因素引起的地震响应, 突出储层含油饱和度变化引起的地震响应后, 可明显看到储层含油饱和度变化引起的地震响应, 且不同含油饱和度范围具有不同的特征。因此, 地震资料携带了储层流体信息, 利用地震资料检测储层含油性及含油量是可行的, 但对储层含油性和饱和度的合理解释需要以正确认识储层含油饱和度变化同地震响应之间的非线性规律为前提基础。

致谢:本论文在完成过程中得到了中国石化石油勘探开发研究院金之钧院士的指导和帮助, 在此表示感谢!

The authors have declared that no competing interests exist.

参考文献
[1] 郑晓东. AVO理论和方法的一些新进展[J]. 石油地球物理勘探, 1992, 27(3): 305-317. [本文引用:1]
[2] 胡中平. 用P-SV波AVO特性判别油气异常[J]. 石油地球物理勘, 1995, 30(5): 653-658. [本文引用:1]
[3] Goodway B, Chen T, Downton J. Improved AVO fluid detection and lithology discrimination using Lame petrophysical parameters from P and S inversions[C]//Expand ed Abstracts of 67th SEG Annual Internet Meeting, 1997: 183-186. [本文引用:1]
[4] 陈小宏, 牟永光. 四维地震油藏监测技术及其应用[J]. 石油地球物理勘探, 1998, 33(6): 707-705. [本文引用:1]
[5] Connolly P. Elastic impedance[J]. The Leading Edge, 1999, 18(4): 438-452. [本文引用:1]
[6] 马劲风. 地震勘探中广义弹性阻抗的正反演[J]. 地球物理学报, 2003, 46(1): 118-124. [本文引用:1]
[7] 黄绪德. 地震勘探直接找油气在国外的最新发展[J]. 勘探地球物理进展, 2004, 27(3): 218-227. [本文引用:1]
[8] Hampson D P, Russell B H. Simultaneous inversion of pre-stack seismic data[C]//Expand ed Abstracts of 75th SEG Annual Internet Meeting, 2005: 1633-1637. [本文引用:1]
[9] 李景叶, 陈小宏, 金龙. 时移地震AVO反演在油藏定量解释中的应用[J]. 石油学报, 2005, 26(3): 68-73. [本文引用:1]
[10] Contreras A, Torres-Verdin C, Fasnacht T. AVA simultaneous inversion of partially stacked seismic amplitude data for the spatial delineation of lithology and fluid units of deepwater hydrocarbon reservoirs in the central Gulf of Mexico[J]. Geophysics, 2006, 71(4): E41-E48. [本文引用:1]
[11] 宁忠华, 贺振华, 黄德济. 基于地震资料的高灵敏度流体识别因子[J]. 石油物探, 2006, 45(3): 239-241. [本文引用:1]
[12] 管路平. 地震叠前反演与直接烃类指示的探讨[J]. 石油物探, 2008, 47(3): 228-234. [本文引用:1]
[13] 孙鹏远. AVO技术研究的新进展[J]. 勘探地球物理进展, 2005, 28(6): 432-438. [本文引用:2]
[14] 王振国, 陈小宏, 王学军, . AVO方法检测油气应用实例分析[J]. 石油地球物理勘探, 2007, 42(2): 194-197. [本文引用:1]
[15] 云美厚, 易维启, 庄红艳. 砂岩的弹性模量与孔隙率、泥质含量、有效压力和温度的经验关系[J]. 石油地球物理勘探, 2001, 36(3): 308-314. [本文引用:1]
[16] Wood A W. A Textbook of Sound[M]. New York: The Maclillan Corporation, 1955: 360. [本文引用:1]
[17] 吴宝年, 吴肃琴, 黄义, . 交错网格在伪谱法弹性波场数值模拟中的应用[J]. 石油物探, 2012, 51(5): 440-445. [本文引用:2]
[18] Fornberg B. The pseudospectral method: Comparisons with finite differences for the elastic wave equatio[J]. Geophysics, 1987, 52(4): 438-501. [本文引用:1]