含高速屏蔽层的地震数据采集及分析
|
摘要
高速层与围岩之间波阻抗差异大，反射系数也大，形成强反射界面，若采用常规的勘探方法，纵波向下入射遇到高速层时，临界角比较小，很小的入射角就产生全反射，地震波难以透射下去达到勘探所要求达到的目的层，所以高速层会对下传的地震波的能量产生_定的屏蔽作用。而且透射不下去的那部分能量会产生了各种干扰波，使得折射波及多次反射波的能量特别强。本文根据Snell定理和Zoeppritz方程，通过射线追踪进行正演模拟，对含高速屏蔽层地区的地震数据采集工作做了一定的研究。
首先，在不同地层参数条件下，较为系统的分析了P波和SV波反射透射系数随入射角变化的规律。从而得出可针对高速层反射透射的特点，相应对其采集方法和采集技术进行改进。
其次，对广角纵波反射勘探作了一定的研究，当地下存在高速层时，采用长排列接收到广角反射波时可以获得较强的能量。但长排列时若用常规的动校正方法就不能满足勘探的精度要求，所以要采用高阶的非双曲时距方程。本文通过理论公式的推导还得出不能准确进行动校正的原因是因为地震波沿射线路径的平均速度和均方根速度都随炮检距的增加而变大。
再次，根据地震波的反射和透射特点知道在一定的地质条件下可以获得相对较强的反射SV波，即在含高速屏蔽层的地质条件下可进行PS转换波勘探。本文通过对一元四次方程和一元三次方程的求解推导出对均匀水平界面可以准确求出共转换点的位置，并对单参数、双参数和三参数时距方程进行了分析和对比，还对PS转换波的AVO特点进行了一定的分析。
The discrepancy of high velocity layer and wall rock in wave impedance or in reflectivity is palpable, so it will come into being a strong reflector. If we adopt a normal method of seismic exploration, the high velocity will shield the energy of down wave, because the critical angle is very small and the full reverberation is easy, also it is difficulty for seismic wave strike the interface which we will explore. Moreover, the energy of the interference wave, such as multiple wave and refractive wave, are very strong. Based on Snell law and Zoeppritz equation and used the forward modeling in ray theory, this master's thesis have do some research on gathering the seismic data in the high velocity shielded layer area.
Firstly we analysis the reflection coefficient and transmission coefficient of P-wave or SV-wave vary with the incidence angle under different condition. So we can ameliorate ways and means in seismic exploration according to the peculiarity of the high velocity shielded layer.
Secondly this thesis study the wide-angel reflection, we can receive stronger energy with a long layout in the high velocity layer area. In this condition, it is inapplicable to NMO correction with normal method, so non-hyperbolic time-distance equation is required. In addition, this thesis also deduce the V_(av) or V_r will be accelerate along with the off-set.
Finally based on the peculiarity of reflection and transmission we can receive comparatively stronger reflected SV wave, so it is feasible to explore with PS converted wave in the high velocity shielded layer area. This thesis has sought exactly the converted point in homogeneous medium by solving the cubic equation and quartic equation, it also analysis and compare one-parameter, two-parameter and three-parameter time-distance equation of PS converted wave. Moreover, it analysis the peculiarity of AVO.
首先，在不同地层参数条件下，较为系统的分析了P波和SV波反射透射系数随入射角变化的规律。从而得出可针对高速层反射透射的特点，相应对其采集方法和采集技术进行改进。
其次，对广角纵波反射勘探作了一定的研究，当地下存在高速层时，采用长排列接收到广角反射波时可以获得较强的能量。但长排列时若用常规的动校正方法就不能满足勘探的精度要求，所以要采用高阶的非双曲时距方程。本文通过理论公式的推导还得出不能准确进行动校正的原因是因为地震波沿射线路径的平均速度和均方根速度都随炮检距的增加而变大。
再次，根据地震波的反射和透射特点知道在一定的地质条件下可以获得相对较强的反射SV波，即在含高速屏蔽层的地质条件下可进行PS转换波勘探。本文通过对一元四次方程和一元三次方程的求解推导出对均匀水平界面可以准确求出共转换点的位置，并对单参数、双参数和三参数时距方程进行了分析和对比，还对PS转换波的AVO特点进行了一定的分析。
The discrepancy of high velocity layer and wall rock in wave impedance or in reflectivity is palpable, so it will come into being a strong reflector. If we adopt a normal method of seismic exploration, the high velocity will shield the energy of down wave, because the critical angle is very small and the full reverberation is easy, also it is difficulty for seismic wave strike the interface which we will explore. Moreover, the energy of the interference wave, such as multiple wave and refractive wave, are very strong. Based on Snell law and Zoeppritz equation and used the forward modeling in ray theory, this master's thesis have do some research on gathering the seismic data in the high velocity shielded layer area.
Firstly we analysis the reflection coefficient and transmission coefficient of P-wave or SV-wave vary with the incidence angle under different condition. So we can ameliorate ways and means in seismic exploration according to the peculiarity of the high velocity shielded layer.
Secondly this thesis study the wide-angel reflection, we can receive stronger energy with a long layout in the high velocity layer area. In this condition, it is inapplicable to NMO correction with normal method, so non-hyperbolic time-distance equation is required. In addition, this thesis also deduce the V_(av) or V_r will be accelerate along with the off-set.
Finally based on the peculiarity of reflection and transmission we can receive comparatively stronger reflected SV wave, so it is feasible to explore with PS converted wave in the high velocity shielded layer area. This thesis has sought exactly the converted point in homogeneous medium by solving the cubic equation and quartic equation, it also analysis and compare one-parameter, two-parameter and three-parameter time-distance equation of PS converted wave. Moreover, it analysis the peculiarity of AVO.
引文
[1] 李正文，贺振华等．勘查技术工程学[M]．北京：地质出版社，2003．
[2] 黄德济，贺振华等．地震勘探资料数字处理[M]．北京：地质出版社，1990．
[3] 贺振华，黄德济等．复杂油气藏地震波场特征方法理论及应用[M]．四川科学技术出版社，1999．
[4] 朱介寿等，地震学中的计算方法[M]．北京：地质出版社，1988．
[5] 陆基孟．地震勘探原理[M]．山东东营：石油大学出版社，2004．
[6] Cerveny V，Molotko I．A著：刘福田，胡戈等译：地震学中的射线方法[M]．北京：地质出版社，1986．
[7] 何樵登，地震勘探原理和方法[M]．北京：地质出版社，1986．
[8] 傅淑芳等：地震学教程[M]．北京：地震出版社，1980．
[9] 王志，贺振华，黄德济等．高速屏蔽层下弱反射层地震勘探[J]．勘探地球物理进展，2002，25(5)，23-27．
[10] 王志，贺振华等．广角反反射特征研究及高斯射线正演模拟分析[J]．成都理工人学学报，2002(1)：56-81．
[11] 王熙明，王志，贺振华等．广角反射方式下观测系统设计与优化[J]．物探化探计算技术，2002，24(4)，298-303．
[12] 李智宏，胡中平，王辉明等．含高速层地震物理模型的广角地震波场研究[J]．勘探地球物理进展，2005，28(5)，345-348．
[13] 姬小兵，尚应军，张帆．山地地震勘探采集技术研究[J]．天然气地球科学，2005，16(1)：114-116．
[14] 郭晓龙，欧阳永林，耿晶等．鄂尔多斯盆地苏里格气田地震转换波解释[J]，天然气地球科学，2005，16(5)：650-653．
[15] 刘洋，魏修成．转换波地震勘探中的最大炮检距设计[J]．勘探地球物理进展，2005，28(3)：247-251．
[16] 刘洋，魏修成．转换波三参数速度分析及动校正方法[J]．石油地球物理探，2005，40(5)，504-509．
[17] 薛冈，王良书等．大炮检距地震资料动校正方法比较[J]．石油地球物理探，2003，38(2)，151-155．
[18] Vlastislav Cerveny. Seismic Ray Theory[M]. The Pitt Building, Trumpington Street, Cambridge, United Kingdom, 2001.
[19] 王建花，李庆忠等．浅层强反射界面的能量屏蔽作用[J]．石油地球物理探，2003，38(6)， 589-502.
[20] Vloritz M, Fliendner, Robert S White. Using Wide-angle seismic data for basalt and sub-basalt imaging[J]. In: Annual Meeting Abstrats[C]. society of Exploration Geophysicists, 1999.
[21] 刘宏伟，董幼瑞．陵北玄武岩发育区广角—转换波地震勘探[J]．物探化探计算技术，2004，26(3)，199-201
[22] 裴正林，牟永光．火成岩区地震波传播规律研究[J]．石油物探，2004，43(5)，433-437．
[23] Purnell G W. Imaging beneath a high velocity layer using converted waves [J]. Geophysics, 1992, 57(4): 1444-1452.
[24] Virieux J. PSV wave propagation in heterogeneous media: Velocity stress finite difference method[J]. Geophysics, 1986, 51(4): 889-9014.
[25] 郑晓东．Zoeppritz方程的近似及其应用[J]．石油地球物理勘探，1991，26(2)：129．
[26] 王小六，李振春等．广角地震采集综述[J]．勘探地球物理进展，2004，27(5)，321-326．
[27] 胡中平，管路平等．高速屏蔽层下广角地震波场分析及成像方法[J]．地球物理学报，2004，47(1)，88-94．
[28] Richard J Castle. A theory of normal moveout. Geophysics, 1994, 59(6): 983-999.
[29] Alkhalifah T, TsvankinI. Velocity analysis for transversely isotropic media. Geophysics, 1995, 60(5): 1550-1566.
[30] Tsvankin Ⅰ and Thomsen I. Non-hyperbolic reflection moveout in anisotropic media. Geophysics, 1994, 59: 1290-1304.
[31] Sun C Wet al. Optimized 6th order NMO correction for long-offset seismic data. Expanded Abstracts of 72nd SEG Mtg, 2002.
[32] 刘洋．双曲时距方程对地震资料处理的影响[J]．石油地球物理勘探，2002 39(2)，134-138
[33] 夏洪瑞，葛川庆等．广角反射高阶NMO校正中存在的问题及对策[J]．江汉石油科技，2005，15(1)，14-19．
[34] Castle R J. A theory of normal moveout. Ceophysics, 1994, 59(6): 983-999.
[35] 刘洋，魏修成．转换波地震勘探的若干问题与对策[J]．勘探地球物理进展，2003，26(4)，247-251．
[36] Stewart R R, Gaiser J E, BrownR J, etal. Converted wave seismic Exploration: Methods[J]Geophysics, 2002, 67(5): 1348-1363.
[37] Cary P W, EatonDWS. A simple method for resolving large converted wave(PSV) statics[J]. Geophysics, 1993, 58(3): 429-433
[38] 朱海波．射线正演研究与采集参数的设计[J]．成都理工大学硕士论文，2004．
[39] 蒋先艺．基于二维与三维复杂结构模型正演的地震数据采集设计方法研究．成都理工大 学博士论文，2003．
[40] 王立明．提高多波速度分析精度的研究．长安大学硕士论文，2005．
[41] 李录明，罗省贤．PSV波转换波速度分析及解释方法[J]．石油地球物理勘探，1995，30(1)：66-74．
[42] 许士勇，马在田．快速有效的转换波共转换点叠加技术[J]．地球物理学报，2002，45(4)：557-568．
[43] 郭向宇，凌云，魏修成．PS转换波共转换点的几种计算方法及实际应用[J]．石油物探，2002，41(2)：141-148．
[44] 毕丽飞，李振春等．转换波处理中共转换点的定位[J]．勘探地球物理进展，2005，28(5)，353-356．
[45] Tessmer G, Behle A. Common reflection point data stacking technique for converted waves[J]. Geophysical Prospecting, 1988, 36: 671-688.
[46] Leon Thomsen. Converted wave reflection seismology over in homogeneous, anisotropic media[J]. Geophysics, 1999, 64(3), 678-690.
[47] Tessmer G. Processing of PS reflection data applying a common conversion point stacking technique[J]. Geophysical Prospecting, 1990, 38: 267-286.
[48] 苑春方，彭苏萍等．水平界面上PSV转换波转换点的精确解[J]．地球物理学报．2005，48(5)，1179-1184．
[49] Ilya Tsvankin and Leon Thornsen Nonhyperbolic reflection moveout in anisotropie media[J]. GEOPHYSICS, VOL.59, NO.8 (AUGUST 1994); P 1290-1304.
[50] 郭刚明，郭树祥等．四阶速度分析与动校正方法及应用效果分析[J]．油气地质与采收率，2003，10(2)，38-39．
[51] Zhang Wenpo, Guo Ping. Study and Practice of Wide-Angle Seismic Data Processing APPLIED GEOPHYSICS, Vol.1, No.1(July 2004), 31-37.
[52] 徐常练，许云等．速度随炮检距变化(VVO)分析．石油地球物理勘探，1998，33(6)，733-74
[53] 张璐，印成耀等．双相介质中AVO正演模拟[J]．地球物理学进展，2005，20(2)，319-322
[54] 阮爱国，李家彪等．海底天然气水合物层界面反射[J]．AVO数值模拟，2006，49(6)，1826-1834
[55] 孙鹏远，孙建国．P-SV波AVO分析[J]．石油地球物理勘探，2003，38(2)，131-135
[2] 黄德济，贺振华等．地震勘探资料数字处理[M]．北京：地质出版社，1990．
[3] 贺振华，黄德济等．复杂油气藏地震波场特征方法理论及应用[M]．四川科学技术出版社，1999．
[4] 朱介寿等，地震学中的计算方法[M]．北京：地质出版社，1988．
[5] 陆基孟．地震勘探原理[M]．山东东营：石油大学出版社，2004．
[6] Cerveny V，Molotko I．A著：刘福田，胡戈等译：地震学中的射线方法[M]．北京：地质出版社，1986．
[7] 何樵登，地震勘探原理和方法[M]．北京：地质出版社，1986．
[8] 傅淑芳等：地震学教程[M]．北京：地震出版社，1980．
[9] 王志，贺振华，黄德济等．高速屏蔽层下弱反射层地震勘探[J]．勘探地球物理进展，2002，25(5)，23-27．
[10] 王志，贺振华等．广角反反射特征研究及高斯射线正演模拟分析[J]．成都理工人学学报，2002(1)：56-81．
[11] 王熙明，王志，贺振华等．广角反射方式下观测系统设计与优化[J]．物探化探计算技术，2002，24(4)，298-303．
[12] 李智宏，胡中平，王辉明等．含高速层地震物理模型的广角地震波场研究[J]．勘探地球物理进展，2005，28(5)，345-348．
[13] 姬小兵，尚应军，张帆．山地地震勘探采集技术研究[J]．天然气地球科学，2005，16(1)：114-116．
[14] 郭晓龙，欧阳永林，耿晶等．鄂尔多斯盆地苏里格气田地震转换波解释[J]，天然气地球科学，2005，16(5)：650-653．
[15] 刘洋，魏修成．转换波地震勘探中的最大炮检距设计[J]．勘探地球物理进展，2005，28(3)：247-251．
[16] 刘洋，魏修成．转换波三参数速度分析及动校正方法[J]．石油地球物理探，2005，40(5)，504-509．
[17] 薛冈，王良书等．大炮检距地震资料动校正方法比较[J]．石油地球物理探，2003，38(2)，151-155．
[18] Vlastislav Cerveny. Seismic Ray Theory[M]. The Pitt Building, Trumpington Street, Cambridge, United Kingdom, 2001.
[19] 王建花，李庆忠等．浅层强反射界面的能量屏蔽作用[J]．石油地球物理探，2003，38(6)， 589-502.
[20] Vloritz M, Fliendner, Robert S White. Using Wide-angle seismic data for basalt and sub-basalt imaging[J]. In: Annual Meeting Abstrats[C]. society of Exploration Geophysicists, 1999.
[21] 刘宏伟，董幼瑞．陵北玄武岩发育区广角—转换波地震勘探[J]．物探化探计算技术，2004，26(3)，199-201
[22] 裴正林，牟永光．火成岩区地震波传播规律研究[J]．石油物探，2004，43(5)，433-437．
[23] Purnell G W. Imaging beneath a high velocity layer using converted waves [J]. Geophysics, 1992, 57(4): 1444-1452.
[24] Virieux J. PSV wave propagation in heterogeneous media: Velocity stress finite difference method[J]. Geophysics, 1986, 51(4): 889-9014.
[25] 郑晓东．Zoeppritz方程的近似及其应用[J]．石油地球物理勘探，1991，26(2)：129．
[26] 王小六，李振春等．广角地震采集综述[J]．勘探地球物理进展，2004，27(5)，321-326．
[27] 胡中平，管路平等．高速屏蔽层下广角地震波场分析及成像方法[J]．地球物理学报，2004，47(1)，88-94．
[28] Richard J Castle. A theory of normal moveout. Geophysics, 1994, 59(6): 983-999.
[29] Alkhalifah T, TsvankinI. Velocity analysis for transversely isotropic media. Geophysics, 1995, 60(5): 1550-1566.
[30] Tsvankin Ⅰ and Thomsen I. Non-hyperbolic reflection moveout in anisotropic media. Geophysics, 1994, 59: 1290-1304.
[31] Sun C Wet al. Optimized 6th order NMO correction for long-offset seismic data. Expanded Abstracts of 72nd SEG Mtg, 2002.
[32] 刘洋．双曲时距方程对地震资料处理的影响[J]．石油地球物理勘探，2002 39(2)，134-138
[33] 夏洪瑞，葛川庆等．广角反射高阶NMO校正中存在的问题及对策[J]．江汉石油科技，2005，15(1)，14-19．
[34] Castle R J. A theory of normal moveout. Ceophysics, 1994, 59(6): 983-999.
[35] 刘洋，魏修成．转换波地震勘探的若干问题与对策[J]．勘探地球物理进展，2003，26(4)，247-251．
[36] Stewart R R, Gaiser J E, BrownR J, etal. Converted wave seismic Exploration: Methods[J]Geophysics, 2002, 67(5): 1348-1363.
[37] Cary P W, EatonDWS. A simple method for resolving large converted wave(PSV) statics[J]. Geophysics, 1993, 58(3): 429-433
[38] 朱海波．射线正演研究与采集参数的设计[J]．成都理工大学硕士论文，2004．
[39] 蒋先艺．基于二维与三维复杂结构模型正演的地震数据采集设计方法研究．成都理工大 学博士论文，2003．
[40] 王立明．提高多波速度分析精度的研究．长安大学硕士论文，2005．
[41] 李录明，罗省贤．PSV波转换波速度分析及解释方法[J]．石油地球物理勘探，1995，30(1)：66-74．
[42] 许士勇，马在田．快速有效的转换波共转换点叠加技术[J]．地球物理学报，2002，45(4)：557-568．
[43] 郭向宇，凌云，魏修成．PS转换波共转换点的几种计算方法及实际应用[J]．石油物探，2002，41(2)：141-148．
[44] 毕丽飞，李振春等．转换波处理中共转换点的定位[J]．勘探地球物理进展，2005，28(5)，353-356．
[45] Tessmer G, Behle A. Common reflection point data stacking technique for converted waves[J]. Geophysical Prospecting, 1988, 36: 671-688.
[46] Leon Thomsen. Converted wave reflection seismology over in homogeneous, anisotropic media[J]. Geophysics, 1999, 64(3), 678-690.
[47] Tessmer G. Processing of PS reflection data applying a common conversion point stacking technique[J]. Geophysical Prospecting, 1990, 38: 267-286.
[48] 苑春方，彭苏萍等．水平界面上PSV转换波转换点的精确解[J]．地球物理学报．2005，48(5)，1179-1184．
[49] Ilya Tsvankin and Leon Thornsen Nonhyperbolic reflection moveout in anisotropie media[J]. GEOPHYSICS, VOL.59, NO.8 (AUGUST 1994); P 1290-1304.
[50] 郭刚明，郭树祥等．四阶速度分析与动校正方法及应用效果分析[J]．油气地质与采收率，2003，10(2)，38-39．
[51] Zhang Wenpo, Guo Ping. Study and Practice of Wide-Angle Seismic Data Processing APPLIED GEOPHYSICS, Vol.1, No.1(July 2004), 31-37.
[52] 徐常练，许云等．速度随炮检距变化(VVO)分析．石油地球物理勘探，1998，33(6)，733-74
[53] 张璐，印成耀等．双相介质中AVO正演模拟[J]．地球物理学进展，2005，20(2)，319-322
[54] 阮爱国，李家彪等．海底天然气水合物层界面反射[J]．AVO数值模拟，2006，49(6)，1826-1834
[55] 孙鹏远，孙建国．P-SV波AVO分析[J]．石油地球物理勘探，2003，38(2)，131-135