高阶高密度三维多波多分量弹性波波场分离正演数值模拟
详细信息 本馆镜像全文    |  推荐本文 | | 获取馆网全文
摘要
针对二维全波场弹性波模拟无法准确表征真实的三维地球介质波场响应,为此首次开展了三维弹性波全波场波场分离正演数值模拟研究,该方法是全面认识和研究地震波在三维地球介质中传播规律的重要手段。从三维一阶速度—应力弹性波波动方程出发,采用高阶(16阶)交错网格有限差分法和内外侧镶边法PML吸收边界条件,并结合了基于散度和旋度原理的三维空间高阶波场分离算子。以三维地质模型为例,设计模拟了野外地表多分量地震资料的采集、零偏移距和非零偏移距VSP记录的采集,开展了三维弹性波全波场高精度、高密度正演数值模拟,并对计算结果进行了三维立体显示。从理论分析和数值试验均表明,文中方法的模拟结果中包括3个多分量混合波场、1个纯纵波波场和3个纯横波波场,在三维弹性波波场快照和模拟记录中均可以容易识别和追踪出各个地层反射界面的反射波和转换波等复杂波场,且计算结果的精度较高,边界吸收较好,采用的三维弹性波波场分离模拟方法还可为实际多波多分量地震资料波场分离方法的准确有效性和三维多波多分量弹性波逆时偏移成像提供方法指导,该方法研究在多波多分量地震勘探中具有一定的推广价值。
For the reason that the 2D full elastic wave simulation method cannot correctly depict the real 3D earth medium’s wave re sponses,3D full elastic wave field separating forward numerical simulation is carried out,which is the important mean of comprehen sively recognize and study the seismic wave’s propagating law in 3D earth medium.This paper started from 3D first order velocitystress elastic wave equation,used high-order(16th) staggered grid finite difference method and PML absorbing boundary condition of inside and outside bordering method and combined with 3D high-order spatial wave separating operator based on divergence and curl law.Taking 3D geologic model as an example,this paper designed and simulated acquisition of the field surface multi-compo nent seismic data,zero offset VSP records and nonzero offset VSP records,and then carried out the high-accuracy high-density for ward numerical simulation of 3D full elastic wave fields,and used 3D visualization on the computational results.The theoretical anal ysis and the numerical experiments showed that the simulation results contained three multi-component hybrid wave field,one pure compress wave field and three pure shear wave field.Every formational interface’s complex wave field such as reflective waves and transformation waves could be obviously recognized and traced within 3D elastic wave snapshots and numerical records.The comput ed results had perfect numerical accuracy and boundary absorbing effect.The adopted 3D elastic wave field separating method could also guide for the accurate and effective of the practical multi-wave multi-component seismic data wave field separating method and the reverse-time migration imagery.This method could be extended in multi-wave multi-component seismic prospecting.
引文
[1]陈可洋.高阶弹性波波动方程正演模拟及逆时偏移成像研究[D].黑龙江大庆:大庆石油学院,2009.
    [2]张文生,宋海斌.三维正交各向异性介质三分量高精度有限差分正演模拟[J].石油地球物理勘探,2001,36(4):422-432.
    [3]裴正林.三维各向异性介质中弹性波方程交错网格高阶有限差分法数值模拟[J].石油大学学报(自然科学版),2004,28(5):23-29.
    [4]裴正林.三维各向同性介质弹性波方程交错网格高阶有限差分法模拟[J].石油物探,2005,44(4):308-315.
    [5]陈可洋.标量声波波动方程高阶交错网格有限差分法[J].中国海上油气,2009,21(4):232-236.
    [6]陈可洋.一阶速度—应力Biot双相各向同性介质弹性波波场分离数值模拟[J].计算物理,2011,28(3):404-412.
    [7]陈可洋,杨微,刘洪林,等.二阶弹性波动方程高精度交错网格波场分离数值模拟[J].物探与化探,2009,33(6):700-704.
    [8]陈可洋.井间弹性波波场散射特征数值模拟分析[J].岩性油气藏,2011,23(3):91-96.
    [9]刘洪林,陈可洋,杨微,等.高阶交错网格有限差分法纵横波波场分离数值模拟[J].地球物理学进展,2010,25(3):877-884.
    [10]陈可洋.高精度数值模拟多分量数据的采集[J].油气地球物理,2010,8(2):44-47.
    [11]陈可洋.基于交错网格的弹性波波场分离数值模拟方法[J].中国石油勘探,2010,14(3):33-37.
    [12]裴正林,牟永光.地震波传播的数值模拟[J].地球物理学进展,2004,19(4):933-941.
    [13]Alford R M,Kelly K R,Boore D M.Accuracy of finite-differ ence modeling of the acoustic wave equation[J].Geophysics,1974,39(6):834-842.
    [14]Marfurt K J.Accuracy of finite difference and finite elementmodeling of the scalar and elastic equations[J].Geophysics,1984,49(5):533-549.
    [15]Dablain M A.The application of high-differencing to the scalarwave equation[J].Geophysics,1986,51(1):54-66.
    [16]Igel H,Morap,Riolletb.Anisotropic wave propagation throughfinite-difference grids[J].Geophysics,1995,60(4):1203-1216.
    [17]刘金平,张向君,董建平.波动方程显式短算子三维正演模拟[J].吉林大学学报(地球科学版),2004,34(1):138-141.
    [18]王润秋,张明,费建博,等.精细积分法三维地震波正演模拟[J].勘探地球物理进展,2006,29(6):394-397.
    [19]何兵寿,魏修成,刘洋.三维波动方程的数值频散关系及其叠前和叠后数值模拟[J].石油大学学报(自然科学版),2001,25(1):67-71.
    [20]熊晓军,贺振华,黄德济.三维波动方程正演及模型应用研究[J].石油物探,2005,44(6):554-556.
    [21]熊高君,贺振华,张毅祥.三维混合延拓一步法波动方程正演模拟[J].成都理工学院学报,1999,26(1):48-51.
    [22]陈可洋.三维快速高精度地震波正演数值模拟方法及其应用[J].天然气勘探与开发,2011,34(3):12-15.
    [23]康玮,程玖兵.横向各向同性介质拟声波方程及其在逆时偏移中的应用[J].地球物理学报,2012,55(3):1033-1045.
    [24]王泽利,何樵登.用邦德矩阵法推导各向异性介质的弹性常数[J].长春科技大学学报,2000,30(4):404-409.
    [25]陈可洋.地震波数值模拟中差分近似的各向异性分析[J].石油物探,2010,49(1):19-22.
    [26]陈可洋.网格剖分及其精度和计算量分析[J].内陆地震,2011,25(1):12-20.
    [27]陈可洋.地震波数值模拟中差分阶数的边界效应[J].高原地震,2011,23(1):20-23.
    [28]陈可洋.边界吸收中镶边法的评价[J].中国科学院研究生院学报,2010,27(2):170-175.
    [29]杨微,陈可洋.加权吸收边界条件的优化设计[J].石油物探,2009,48(3):244-246,251.
    [30]陈可洋.完全匹配层吸收边界条件研究[J].石油物探,2010,49(5):472-477.
    [31]陈可洋.声波完全匹配层吸收边界条件的改进算法[J].石油物探,2009,48(1):76-79.
    [32]陈可洋.完全匹配层吸收边界条件在弹性波波场分离数值模拟中的应用[J].石油工业计算机应用,2010,65(1):10-12.
    [33]Fei T,Larner K.Elimination of numerical dispersion in finite-difference modeling and migration by flux-corrected transport[J].Geophysics,1995,60(6):1830-1842.
    [34]董良国,李培明.地震波传播数值模拟中的频散问题[J].天然气工业,2004,24(6):53-56.
    [35]吴国忱,王华忠.波场模拟中的数值频散分析与校正策略[J].地球物理学进展,2005,20(1):58-65.
    [36]宁刚,熊章强,陈持.波动方程有限差分正演模拟误差来源分析[J].物探与化探,2008,32(2):203-206.
    [37]董良国,马在田,曹景忠,等.一阶弹性波方程交错网格高阶差分解法[J].地球物理学报,2000,43(3):411-419.
    [38]董良国,马在田,曹景忠.一阶弹性波方程交错网格高阶差分解法稳定性研究[J].地球物理学报,2000,43(6):856-864.
    [39]陈可洋.基于拉普拉斯算子的叠前逆时噪声压制方法[J].岩性油气藏,2011,23(5):87-95.
    [40]陈可洋.高精度地震纯波震源数值模拟[J].岩性油气藏,2012,24(1):84-91.
    [41]陈可洋.地震波数值模拟中优化的通量校正传输方法[J].内陆地震,2012,26(2):169-179.

版权所有:© 2023 中国地质图书馆 中国地质调查局地学文献中心