2. 中国石油化工股份有限公司科技部油田处, 北京 100728;
3. 中国石油化工股份有限公司石油勘探开发研究院, 北京 100083
2. Science and Technology Department, Sinopec, Beijing 100728, China;
3. Sinopec Petroleum Exploration and Production Research Institute, Beijing 100083, China
在我国, 无论是油气勘探较为成熟的东部地区, 还是勘探程度相对较低的西部各大盆地, 岩性油气藏逐步成为油气勘探的主要对象之一。对我国陆相沉积盆地而言, 岩性油气藏的地层厚度一般低于常规地震勘探的垂向分辨率, 以薄层的形式存在[1]。我国煤系地层厚度多在2~10m, 且常以薄层的形式存在, 反射响应非单一界面产生, 而是顶底反射、层间多次波叠加而成的复合波[2]。因此, 研究薄层地震反透射理论, 对油气、煤炭等资源的开发具有一定的理论及现实意义。
对薄层反射透射的理论研究可追溯到20世纪60年代。布列霍夫斯基赫等[3]从弹性波动力学理论的角度讨论了平面波在层状介质中的反射和透射情况, 建立了传播矩阵方法; 之后WIDESS[4]忽略层间多次波及转换波的影响, 建立薄层反射振幅与薄层厚度的关系式; KENNETT[5-6]给出球面波入射情况下层间反射系数的递归关系; 汪恩华等[7-8]给出纵波入射和全固态介质假设条件下薄层的反透射公式, 并进一步给出单个薄层垂直入射时的纵波反射系数谱计算公式; 21世纪后薄层反射理论研究发展迅速, 郭智奇等[9]结合薄层反射动力学及粘弹性理论研究了薄储层物性与AVO属性的关系。YANG等[10]给出薄层厚度较小时单薄层反透射系数的拟Zoeppritz矩阵方程, 并给出小角度及弱阻抗差前提下的单薄层PP波反射的近似解析表达式, 但该近似式忽略了层间转换波及多次波的影响[11]; 张军华等[12]梳理了薄层叠后地震解释技术, 建立了多个模型并对其效果进行了评述; 彭军等[13]使用基追踪方法识别薄层模型并在实际油气田中检测应用效果; 马跃华等[14]建立了“岩性组合-阻抗组合-地质成因”的多维度薄层分类标准, 并以大港油田的两个应用研究为例, 分析了薄层的地震响应特征; 杨春等[15]结合Gassmann方程与薄层反射透射矩阵方程分析了含油气薄砂储层的地震反射特征。以上研究在很大程度上推动了薄层反透射理论的发展, 然而, 所给出的薄层反透射系数公式有些需要进行复杂的矩阵运算, 有些忽略了层间多次波及转换波的影响, 对薄层AVO分析及反演具有一定的局限性。
由于薄层精准反透射系数方程较为复杂, 无法直接应用于实际的薄层反演, 因此, 需要给出形式较为简洁的近似公式。参照单界面反透射系数近似的研究方法[16-18], 本文讨论了薄层反透射系数与射线参数的关系, 对薄层转换波和非转换波与射线参数的奇偶关系予以证明, 并给出薄层反透射系数的低阶近似形式。在此基础上, 分别对给定的两类典型薄储层模型(强阻抗差薄层、弱阻抗差薄层)的反射系数精确解与低阶近似解进行了误差分析, 讨论了小角度范围内(≤30°)射线参数高次项舍弃的可行性与舍弃后的低阶近似精度。
1 薄层反透射系数与射线参数的关系P波入射时, 薄层反透射系数矩阵可写为:
| $ \left[\begin{array}{llll}{m_{11}} & {m_{12}} & {m_{13}} & {m_{14}} \\ {m_{21}} & {m_{22}} & {m_{23}} & {m_{24}} \\ {m_{31}} & {m_{32}} & {m_{33}} & {m_{34}} \\ {m_{41}} & {m_{42}} & {m_{43}} & {m_{44}}\end{array}\right] \cdot\left[\begin{array}{c}{R_{\mathrm{PP}}} \\ {R_{\mathrm{PS}}} \\ {T_{\mathrm{PP}}} \\ {T_{\mathrm{PS}}}\end{array}\right]=\left[\begin{array}{c}{n_{1}} \\ {n_{2}} \\ {n_{3}} \\ {n_{4}}\end{array}\right] $ | (1) |
其中, RPP, RPS, TPP和TPS分别为薄层PP波和PS波的反射和透射系数; mij为反透射系数矩阵方程组的系数矩阵元素; ni为反透射系数矩阵方程组的常数项列矩阵元素; mij和ni由目标薄层的弹性参数、薄层厚度、入射波频率、入射角决定, 具体表达式见参考文献[19]。
将入射角记为θ, 且3层介质中纵波传播方向与法线的夹角分别记作θ1, θ2, θ3, 横波传播方向与法线的夹角分别记作δ1, δ2, δ3。根据Snell定律, 有:
| $ \frac{\sin \theta}{v_{\mathrm{Pl}}}=\frac{\sin \theta_{1}}{v_{\mathrm{Pl}}}=\frac{\sin \theta_{2}}{v_{\mathrm{P} 2}}=\frac{\sin \theta_{3}}{v_{\mathrm{P} 3}}=\frac{\sin \delta_{1}}{v_{\mathrm{Sl}}}=\frac{\sin \delta_{2}}{v_{\mathrm{S} 2}}=\frac{\sin \delta_{3}}{v_{\mathrm{S} 3}} $ | (2) |
因此, 当sinθ取值-sinθ时, sinθ1, sinθ2, sinθ3, sinδ1, sinδ2, sinδ3相应取值分别为-sinθ1, -sinθ2, -sinθ3, -sinδ1, -sinδ2, -sinδ3; sin2θ1, sin2θ2, sin2θ3, sin2δ1, sin2δ2, sin2δ3相应取值分别为-sin2θ1, -sin2θ2, -sin2θ3, -sin2δ1, -sin2δ2, -sin2δ3; tanθ1, tanθ2, tanθ3, tanδ1, tanδ2, tanδ3相应取值分别为-tanθ1, -tanθ2, -tanθ3, -tanδ1, -tanδ2, -tanδ3; cosθ1, cosθ2, cosθ3, cosδ1, cosδ2, cosδ3, cos2θ1, cos2θ2, cos2θ3, cos2δ1, cos2δ2, cos2δ3取值不变。设此时PP波和PS波的反射和透射系数分别为R′PP, R′PS, T′PP, T′PS, 将它们代入方程(1)并整理可得:
| $ \left[\begin{array}{llll}{m_{11}} & {m_{12}} & {m_{13}} & {m_{14}} \\ {m_{21}} & {m_{22}} & {m_{23}} & {m_{24}} \\ {m_{31}} & {m_{32}} & {m_{33}} & {m_{34}} \\ {m_{41}} & {m_{42}} & {m_{43}} & {m_{44}}\end{array}\right] \cdot\left[\begin{array}{c}{R_{\mathrm{PP}}^{\prime}} \\ {-R_{\mathrm{PS}}^{\prime}} \\ {T_{\mathrm{PP}}^{\prime}} \\ {-T_{\mathrm{PS}}^{\prime}}\end{array}\right]=\left[\begin{array}{l}{n_{1}} \\ {n_{2}} \\ {n_{3}} \\ {n_{4}}\end{array}\right] $ | (3) |
比较(1)式和(3)式, 可以看出:
| $ \begin{array}{l}{R_{\mathrm{PP}}^{\prime}=R_{\mathrm{PP}}} \\ {R_{\mathrm{PS}}^{\prime}=-R_{\mathrm{PS}}} \\ {T_{\mathrm{PP}}^{\prime}=T_{\mathrm{PP}}} \\ {T_{\mathrm{PS}}^{\prime}=-T_{\mathrm{PS}}}\end{array} $ | (4) |
由此可知, 薄层非转换波的反射和透射系数是sinθ(或射线参数)的偶函数, 而转换波的反射和透射系数是sinθ(或射线参数)的奇函数。
利用反射和透射系数的奇偶性可以将其表示成以下sinθ的幂级数形式:
| $ \left\{\begin{array}{l}{R_{\mathrm{PP}}=\sum\limits_{n=0}^{\infty} A R_{2 n} \sin ^{2 n} \theta} \\ {R_{\mathrm{PS}}=\sum\limits_{n=0}^{\infty} A R_{2 n+1} \sin ^{2 n+1} \theta} \\ {T_{\mathrm{PP}}=\sum\limits_{n=0}^{\infty} A T_{2 n} \sin ^{2 n} \theta} \\ {T_{\mathrm{PS}}=\sum\limits_{n=0}^{\infty} A T_{2 n+1} \sin ^{2 n+1} \theta}\end{array}\right. $ | (5) |
式中:AR2n, AR2n+1, AT2n, AT2n+1分别为薄层PP波反射系数、PS波反射系数、PP波透射系数及PS波透射系数的幂级数展开式的系数。入射角较小时, sinnθ随着级数n的增大而迅速衰减, 则近似求解出参数AR2n, AR2n+1, AT2n, AT2n+1, 并根据精度需求对(5)式进行高阶舍弃, 从而获得相应的薄层反射或透射系数的近似公式。
2 薄层反射系数的低阶拟合针对油气及煤资源勘探面临的薄储层模型, 我们分别建立了强、弱阻抗差的薄层模型, 并讨论n≤1情况下两种模型的反射系数低阶近似与精确解的差异。表 1给出了强阻抗差的薄煤层和弱阻抗差的含油气薄砂岩层的模型参数。设定入射纵波主频为50Hz, 计算小角度范围内(入射角θ≤30°)不同厚度下薄层PP波、PS波的反射系数。薄层厚度依次取值为:λ/4, λ/8, λ/10, λ/20, λ/30, λ/40, λ/60, λ/80, λ/100, 0, 其中, λ为纵波波长, 0代表薄层上覆围岩及下伏围岩构成的单界面特例, 顶界面、底界面分别代表薄层顶界面、底界面的反射情况。鉴于薄层反射系数为复反射系数, 而实际地震勘探中接收到的为实地震信号, 所以以下仅对薄层反射系数的实部进行讨论。对PP波反射系数的实部进行sin2θ线性拟合, PS波反射系数的实部进行sinθ线性拟合, 并与薄层顶、底界面反射系数情况进行对比。
|
表 1 薄煤层和薄砂层的模型参数 |
图 1至图 3给出了不同厚度下含油气砂岩模型的RPP关于sin2θ的线性拟合结果。由图 1至图 3可以看出, 随着薄层厚度的减小, 含油气砂岩层的RPP实部随sin2θ的增加其变化趋势发生改变:当薄砂层厚度小于等于λ/10时, RPP的实部随sin2θ的增加而减小, 与薄层顶界面和底界面反射系数的实部随sin2θ的变化规律一致, 说明在此厚度范围内, 顶界面与底界面的阻抗差对反射系数的影响较大; 厚度大于等于λ/8时, RPP的实部随着sin2θ的增加而增大。当入射角θ≤30°时, 薄层反射系数RPP的实部对sin2θ线性程度很高, 特别是厚度范围为(λ/60, λ/20)时; λ/8和λ/10为RPP实部随sin2θ的增加而变化趋势发生改变的界限点, 反射系数的实部随sin2θ变化的线性程度较其它厚度低, 但对sin2θ线性拟合的确定系数R2也大于等于0.910。薄煤层模型呈现规律大体上与薄含油气砂岩一致, 考虑篇幅限制, 我们只给出了含油气砂岩层模型的RPP实部部分厚度的拟合图。图 4给出了两个模型不同厚度及顶、底界面反射系数线性拟合(RPP的实部关于sin2θ的线性拟合, RPS的实部关于sinθ的截距为0的线性拟合)的确定系数R2。由图 4可以看出, 含油气砂岩层模型厚度为λ/10的RPP实部拟合的确定系数最小, 但也高达0.912。说明在小角度范围(入射角≤30°)内, 不论是强阻抗差的煤层模型还是弱阻抗差的含油气砂岩模型, 反射系数在n≤1情况下的拟合程度都很高(R2≥0.912)。
|
图 1 厚度为λ/4(a), λ/8(b)的含油气砂岩模型的RPP关于sin2θ的线性拟合结果 |
|
图 2 厚度为λ/10(a), λ/20(b), λ/30(c), λ/40(d), λ/60(e), λ/80(f), λ/100(g)的含油气砂岩模型的RPP关于sin2θ的线性拟合结果 |
|
图 3 厚度为0(a)及顶(b)、底(c)界面含油气砂岩模型的RPP关于sin2θ的线性拟合结果 |
|
图 4 不同厚度下反射系数的确定系数R2 |
将拟合曲线斜率记为G, 截距记为P, 计算不同厚度下薄层与薄层厚度减至0时上覆围岩与下伏围岩组成的单界面、薄层顶界面、薄层底界面的反射系数关于|ΔG|及|ΔP|的相对变化(图 5)。由图 5可以看出, 对于含油气砂岩层模型, 当薄层厚度较大时(λ/4, λ/8, λ/10), 反射系数的实部更接近于顶界面情况, 说明薄层与上覆围岩介质参数对薄层反射系数的影响较大; 当厚度较小(λ/20, λ/30, λ/40, λ/60, λ/80, λ/100)时, 反射系数的实部更接近于薄层厚度为0时的情况, 即薄层上、下围岩的介质参数共同影响薄层反射系数。由于λ为纵波波长, 而横波波长(λS)满足:λ/8 < λS/4 < λ/4, 因此, 对横波而言仅考虑厚度小于等于λ/8的情况:薄层厚度较大(λ/8, λ/10)时, RPS的实部更接近于顶界面情况; 薄层厚度较小(λ/20, λ/30, λ/40, λ/60, λ/80, λ/100)时, RPS的实部更接近于薄层厚度为0时的情况。对于煤层模型计算结论与含油气砂岩层模型相似。
|
图 5 不同厚度下含油气砂岩模型相对不同界面的RPP, RPS关于|ΔG|及|ΔP|的变化情况 a RPP关于|ΔG|的相对变化; b RPP关于|ΔP|的相对变化; c RPS关于|ΔG|的相对变化 |
利用反射和透射系数的奇偶性可以将其表示成sinθ的幂级数形式, 薄层非转换波的反射和透射系数是sinθ(或射线参数)的偶函数, 而转换波的反射和透射系数是sinθ(或射线参数)的奇函数; 级数n≤1情况下两种薄层模型的反射系数低阶近似与精确解的拟合结果说明, 对薄层PP波反射系数的实部进行sin2θ线性拟合、PS波反射系数的实部进行sinθ线性拟合, 都可以获得较高的拟合精度; 薄层上覆、下伏围岩的介质参数对薄层反射系数的影响较大。
本文仅给出了薄层转换波及非转换波反透射系数与射线参数的奇偶性关系, 公式(5)中参数AR2n, AR2n+1, AT2n, AT2n+1的具体表达式并未给出, 需要进一步研究给出完整的反透射系数表达式, 并根据精度需求对公式(5)进行高阶舍弃获得相应的薄层近似公式。在下一步薄层AVO属性分析中, 可以利用以上结论绘制薄层AVO的PG交会图, 以获得薄层属性。本文仅讨论了煤层和含油气砂岩层模型, 在未来地震勘探中还需要讨论更多薄层模型得到普适性的结论。
| [1] |
张玉芬.影响薄互层地震反射特征因素研究和井震联合地震道多尺度反演[D].武汉: 中国地质大学(武汉), 2000 ZHANG Y F.Study on the characteristics of seismic reflections affecting thin interbeds and multi-scale inversion of well-seismic seismic traces[D].Wuhan: China University of Geosciences (Wuhan), 2000 http://www.wanfangdata.com.cn/details/detail.do?_type=degree&id=Y411764 |
| [2] |
张爱敏, 汪洋, 赵世尊. 不同厚度煤层AVO特征及模型研究[J]. 中国矿业大学学报, 1997, 26(3): 36-41. ZHANG A M, WANG Y, ZHAO S Z. Study on the AVO model and AVO character of coal seams with different thickness[J]. Journal of China University of Mining & Technology, 1997, 26(3): 36-41. DOI:10.3321/j.issn:1000-1964.1997.03.009 |
| [3] |
布列霍夫斯基赫, 杨训仁. 分层介质中的波[M]. 北京: 科学出版社, 1985: 31-36. BREKHOVSKIKH, YANG X R. Waves in layered media[M]. Beijing: Science Press, 1985: 31-36. |
| [4] |
WIDESS M B. How thin is a thin bed?[J]. Geophysics, 1973, 38(6): 1176-1180. DOI:10.1190/1.1440403 |
| [5] |
KENNETT B. Reflections, rays and reverberations[J]. Bulletin of the Seismological Society of America, 1974, 64(6): 1685-1696. |
| [6] |
KENNETT B. Theoretical reflection seismograms for elastic media[J]. Geophysical Prospecting, 1979, 27(2): 301-321. DOI:10.1111/j.1365-2478.1979.tb00972.x |
| [7] |
汪恩华. 适用于AVO分析的层状介质反射系数公式及模型计算[J]. 石油地球物理勘探, 1991, 26(4): 424-436. WANG E H. Reflection coefficient formula and model calculation which are applicable for AVO analysis in the case of layered medium[J]. Oil Geophysical Prospecting, 1991, 26(4): 424-436. |
| [8] |
汪恩华, 贺振华, 李庆忠. 基于薄层的反射系数谱理论与模型正演[J]. 成都理工学院学报, 2001, 28(1): 70-74. WANG E H, HE Z H, LI Q Z. Reflectance spectrum theory and model computation based on thin beds[J]. Journal of Chengdu University of Technology, 2001, 28(1): 70-74. DOI:10.3969/j.issn.1671-9727.2001.01.014 |
| [9] |
郭智奇, 刘财, 冯晅, 等. 薄储层的反射特性及其AVO属性分析[J]. 石油物探, 2009, 48(5): 453-458. GUO Z Q, LIU C, FENG X, et al. Reflection characteristics of thin reservoirs and its AVO attributes analysis[J]. Geophysical Prospecting for Petroleum, 2009, 48(5): 453-458. DOI:10.3969/j.issn.1000-1441.2009.05.004 |
| [10] |
YANG C, WANG Y, WANG Y H. Reflection and transmission coefficients of a thin bed[J]. Geophysics, 2016, 81(5): N31-N39. DOI:10.1190/geo2015-0360.1 |
| [11] |
YANG C, WANG Y, LU J. Weak impedance difference approximations of thin-bed PP-wave reflection response[J]. Journal of Geophysics and Engineering, 2017, 14(4): 1010-1019. DOI:10.1088/1742-2140/aa6dd8 |
| [12] |
张军华, 王庆峰, 张晓辉, 等. 薄层和薄互层叠后地震解释关键技术综述[J]. 石油物探, 2017, 56(4): 459-471. ZHANG J H, WANG Q F, ZHANG X H, et al. Poststack interpretation key techniques for thin layer and thin interbed reservoirs[J]. Geophysical Prospecting for Petroleum, 2017, 56(4): 459-471. DOI:10.3969/j.issn.1000-1441.2017.04.001 |
| [13] |
彭军, 周家雄, 王宇, 等. 基追踪在薄层识别中的研究与应用[J]. 地球物理学进展, 2017, 32(3): 1243-1250. PENG J, ZHOU J X, WANG Y, et al. Base pursuit study and application in the identification thin layer[J]. Progress in Geophysics, 2017, 32(3): 1243-1250. |
| [14] |
马跃华, 周宗良, 李振永, 等. 薄层分类及其地震响应分析——以大港油田两个应用研究为例[J]. 石油物探, 2018, 57(6): 902-913. MA Y H, ZHOU Z L, LI Z Y, et al. Classification and seismic response analysis of thin bed:Two cases study from Dagang oilfield[J]. Geophysical Prospecting for Petroleum, 2018, 57(6): 902-913. DOI:10.3969/j.issn.1000-1441.2018.06.013 |
| [15] |
杨春, 王赟, 陈业全, 等. 双相介质薄砂层的地震反射特征分析[J]. 石油物探, 2018, 57(2): 186-197. YANG C, WANG Y, CHEN Y Q, et al. Characteristic analysis of seismic reflections from a thin sand reservoir in two-phase media[J]. Geophysical Prospecting for Petroleum, 2018, 57(2): 186-197. DOI:10.3969/j.issn.1000-1441.2018.02.003 |
| [16] |
郑晓东. Zoeppritz方程的近似及其应用[J]. 石油地球物理勘探, 1991, 26(2): 129-144. ZHENG X D. Approximation of Zoeppritz equation and its application[J]. Oil Geophysical Prospecting, 1991, 26(2): 129-144. |
| [17] |
郑晓东. 转换波和非转换波与射线参数P的关系[J]. 地球物理学报, 1991, 34(6): 781-787. ZHENF X D. Relationships between converted and non-converted waves and parameter P[J]. Chinese Journal of Geophysics, 1991, 34(6): 781-787. DOI:10.3321/j.issn:0001-5733.1991.06.013 |
| [18] |
杨绍国. Zoeppritz方程的级数表达式及近似[J]. 石油地球物理勘探, 1994, 29(4): 399-412. YANG S G. Series expression and approximation of Zoeppritz equation[J]. Oil Geophysical Prospecting, 1994, 29(4): 399-412. |
| [19] |
张川, 杨春, 王赟. 关于薄层与单界面模型弹性反透射系数的讨论[J]. 煤田地质与勘探, 2015, 43(2): 86-90. ZHANG C, YANG C, WANG Y. Discussion on elastic reflection and transmission coefficients of thin-bed and single interface models[J]. Coal Geology & Exploration, 2015, 43(2): 86-90. DOI:10.3969/j.issn.1001-1986.2015.02.018 |