裂缝介质弹性波场数值模拟及波场特征NUMERICAL MODELING AND CHARACTERISTIC ANALYSES OF THE ELASTIC WAVEFIELDS IN THE FRACTURE MEDIA
李昂,张丽艳
LI Ang,ZHANG Liyan
摘要(Abstract):
裂缝介质的数值模拟是研究裂缝介质中地震波传播特征及规律的主要手段。以各向异性介质的弹性波动方程为基础,推导出空间交错网格方法一阶导数的任意偶数阶精度差分格式,并利用此差分格式对裂缝介质进行了2.5D三分量正演模拟。模拟结果表明:受裂缝影响弹性波在裂缝介质中传播的波场是非常复杂的,除了存在反射P波、反射P-SV波,还存在横波分裂现象和沿裂缝介质传播的导波等;数值模拟的三分量的方位道集呈现出波场传播的时间和振幅与裂缝的方位存在很大的关联,其方位各向异性特征与实际资料的方位各向异性特征一致,可以指导今后的方位各向异性处理和裂缝预测。
The numerical modeling of Fracture media is the main method to study the propagation characteristics and rules of the seismic waves in underground fracture media. On the basis of the elastic wave equation of the anisotropic media,the scheme of the precision difference is presented for random even order of the first order of the spatial staggered grid formulations,and moreover with the help of the above scheme,the forward modeling for 2. 5D components is conducted for the fracture media. The simulation results show that due to the influences of the fractures,the wave fields of the elastic waves propagating in the fracture media are rather complex: besides there are P- wave and P- SV wave,the following phenomena and so forth occur: the shear wave splitting and guided wave propagated along the fracture media; there are great connection with the fracture azimuths for the time and amplitude of the wave field propagation in the three- component azimuth gathers of the numerical simulation,their azimuth anisotropy characteristics are consistent with those of the real data,thus the future azimuth anisotropy processing and fracture predication can be guided.
关键词(KeyWords):
裂缝介质;弹性波;交错网格;正演模拟;方位各向异性
fracture medium;elastic wave;staggered grid;forward modeling;azimuth anisotropy
基金项目(Foundation): 中国石油天然气股份有限公司重大专项“大庆探区非常规油气地震勘探技术研究”(2012E-2603-04)
作者(Author):
李昂,张丽艳
LI Ang,ZHANG Liyan
参考文献(References):
- [1]郭锐,王尚旭,余刚.基于高频恢复和多尺度分解的裂缝分析方法[J].东北石油大学学报,2014,38(3):49-56.
- [2]袁书坤,陈开远,Hardage B A,等.多分量地震技术在裂缝储层预测中的应用[J].大庆石油地质与开发,2014,33(6):141-145.
- [3]张利萍,潘仁芳,樊太亮.基于P波方位AVO技术的火成岩裂缝监测[J].大庆石油地质与开发,2011,30(4):144-148.
- [4]董良国,马在田,曹景忠,等.一阶弹性波方程交错网格高阶差分解法[J].地球物理学报,2000,43(3):412-418.
- [5]董良国,马在田,曹景忠,等.一阶弹性波方程交错网格高阶差分解法稳定性研究[J].地球物理学报,2000,43(6):856-864.
- [6]常君勇.交错网格法有限差分数值模拟[D].北京:中国地质大学,2007.
- [7]陈可洋,杨微,刘洪林,等.二阶弹性波动方程高精度交错网格波场分离数值模拟[J].物探与化探,2009,33(6):700-703.
- [8]周学明,李庆春,李志华,等.地震波交错网格高阶差分数值模拟研究[J].铁道学报,2011,8(1):1-6.
- [9]李振春,张华,刘庆敏,等.弹性波交错网格高阶有限差分法波场分离数值模拟[J].石油地球物理勘探,2007,42(5):510-515.
- [10]李军,苏云,李录明.各向异性黏弹性介质波场数值模拟[J].大庆石油学院学报,2009,33(6):38-42.
- [11]Liu Y,Sen M K.An implicit staggered-grid finite-difference method for seismic modeling[J].Geophysical Journal International,2009,179(1):459-474.
- [12]裴正林.层状各向异性介质中横波分裂和再分裂数值模拟[J].石油地球物理勘探,2006,41(1):17-25.
- [13]狄帮让,裴正林.岩丘模型弹性波方程正演模拟及波场特征分析[J].石油地球物理勘探,2010,45(6):826-832.
- [14]Pei Zhenglin,Fu Liyun,Sun Weijia,et al.Anisotropic finitedifference algorithm for modeling elastic wave propagation in fractured coalbeds[J].Geophysics,2012,77(1):13-26.
- [15]裴正林.三维各向异性介质中弹性波方程交错网格高阶有限差分法数值模拟[J].石油大学学报(自然科学版),2004,28(5):23-29.
- [16]裴正林.任意起伏地表弹性波方程交错网格高阶有限差分法数值模拟[J].石油地球物理勘探,2004,39(6):629-634.
- [17]裴正林,董玉珊,彭苏萍.裂隙煤层弹性波场方位各向异性特征数值模拟研究[J].石油地球物理勘探,2007,42(6):665-672.
- [18]Liu Y,Sen M K.Scalar wave equation modeling with time-space domain dispersion-relation-based staggered-grid finite-difference schemes[J].Bulletin of the Seismological Society of America,2011,101(1):141-159.
- [19]Liu Y,Sen M K.Finite-difference modeling with adaptive variablelength spatial operators[J].Geophysics,2011,76(4):T79-T89.
- [20]Alterman Z,Karal F C.Propagation of seismic wave in layered media by finite difference methods[J].Bulletin of Seismological Society of America,1968,58(2):367-398.
- [21]田东升,王云专,李义鹏等.单程和双程波动方程叠前深度偏移方法[J].东北石油大学学报,2014,38(4):39-44.