非均匀TI介质P-SV波传播交错网格高阶有限差分数值模拟
|
摘要
给出了在非均匀横向各向同性(TI)介质情况下,四阶时间精度、高阶空间精度的一阶速度-应力P-SV波的波动方程交错网格有限差分解法。首先根据一阶速度(应力)波动方程把速度(应力)对时间的一阶和三阶导数转换为应力(速度)对空间的导数,从而在使用四阶时间精度有限差分格式计算某一时刻的波场时只需要前面两个时间步的波场值;然后在空间上采用高阶有限差分格式以提高数值模拟的精度。数值模拟结果和实测垂直地震剖面(VSP)记录符合得很好,说明该方法是可行的。
A set of stagger-grid finite-difference operations with 4-order temporal accuracy and high order spatial accuracy to one-order velocity-stress P-SV wave equations is presented in heterogeneous transversely isotropic(TI) media.First using the one-order velocity-stress P-SV wave equations,the first and third order temporal derivatives of particle velocity/stress are transformed into spatial derivatives of stress/particle velocity,thus only two former time-step wave fields are needed to compute wave fields for the current time step with the 4-order temporal accuracy finite-differente approximation.Then the high-order spatial finite-difference approximation is used to improve numerical modeling precise.High consistency between the modeling vertical seismic profiles(VSP) records and the field ones demonstrates well feasibility of the present technique.
A set of stagger-grid finite-difference operations with 4-order temporal accuracy and high order spatial accuracy to one-order velocity-stress P-SV wave equations is presented in heterogeneous transversely isotropic(TI) media.First using the one-order velocity-stress P-SV wave equations,the first and third order temporal derivatives of particle velocity/stress are transformed into spatial derivatives of stress/particle velocity,thus only two former time-step wave fields are needed to compute wave fields for the current time step with the 4-order temporal accuracy finite-differente approximation.Then the high-order spatial finite-difference approximation is used to improve numerical modeling precise.High consistency between the modeling vertical seismic profiles(VSP) records and the field ones demonstrates well feasibility of the present technique.
引文
[1]Madariaga R.Dynamic of an expanding circular fault[J].Bull.Seism.Soc.Amer.,1976,66:639-667.
[2]Virieux J.SH-wave propagation in heterogeneous media:ve-locity-stress finite-difference method[J].Geophysics,1984,49(11):1933-1942.
[3]Virieux J.P-SV wave propagation in heterogeneous media:Velocity-stress finite difference method[J].Geophysics,1986,51(4):889-901.
[4]Levander A R.Fourth-order finite-difference P-SV seismo-grams[J].Geophysics,1988,53(11):1 425-1436.
[5]Dai N,Vafidis A,Kanasewich E R.Wave propagation in het-erogeneous,porous media:A velocity-stress,finite-differencemethod[J].Geophysics,1995,60(2):327-340.
[6]董良国,马在田,曹景忠,等.一阶弹性波方程交错网格高阶差分解法[J].地球物理学报,2000,43(3):411-419.
[7]王秀明,张海澜,王东.利用高阶交错网格有限差分法模拟地震波在非均匀孔隙介质中的传播[J].地球物理学报,2003,46(6):842-849.
[8]裴正林,牟永光.非均匀介质地震波传播交错网格高阶有限差分法模拟[J].石油大学学报(自然科学版),2003,27(6):17-21.
[9]裴正林.三维双相各向异性介质中弹性波方程交错网格高阶有限差分法模拟[J].石油大学学报(自然科学版),2004,28(5):23-29.
[10]裴正林,王尚旭.任意倾斜各向异性介质中弹性波波场交错网格高阶有限差分法模拟[J].地震学报,2005,27(4):441-451.
[11]董良国,马在田,曹景忠.一阶弹性波方程交错网格高阶差分解法稳定性研究[J].地球物理学报,2000,43(6):856-864.
[2]Virieux J.SH-wave propagation in heterogeneous media:ve-locity-stress finite-difference method[J].Geophysics,1984,49(11):1933-1942.
[3]Virieux J.P-SV wave propagation in heterogeneous media:Velocity-stress finite difference method[J].Geophysics,1986,51(4):889-901.
[4]Levander A R.Fourth-order finite-difference P-SV seismo-grams[J].Geophysics,1988,53(11):1 425-1436.
[5]Dai N,Vafidis A,Kanasewich E R.Wave propagation in het-erogeneous,porous media:A velocity-stress,finite-differencemethod[J].Geophysics,1995,60(2):327-340.
[6]董良国,马在田,曹景忠,等.一阶弹性波方程交错网格高阶差分解法[J].地球物理学报,2000,43(3):411-419.
[7]王秀明,张海澜,王东.利用高阶交错网格有限差分法模拟地震波在非均匀孔隙介质中的传播[J].地球物理学报,2003,46(6):842-849.
[8]裴正林,牟永光.非均匀介质地震波传播交错网格高阶有限差分法模拟[J].石油大学学报(自然科学版),2003,27(6):17-21.
[9]裴正林.三维双相各向异性介质中弹性波方程交错网格高阶有限差分法模拟[J].石油大学学报(自然科学版),2004,28(5):23-29.
[10]裴正林,王尚旭.任意倾斜各向异性介质中弹性波波场交错网格高阶有限差分法模拟[J].地震学报,2005,27(4):441-451.
[11]董良国,马在田,曹景忠.一阶弹性波方程交错网格高阶差分解法稳定性研究[J].地球物理学报,2000,43(6):856-864.
版权所有：© 2023 中国地质图书馆 中国地质调查局地学文献中心