基于螺旋线上谱因式分解的地震波场隐式辛算法
|
摘要
均匀介质、复杂各向同性介质和各向异性介质中的地震波传播过程 ,可用统一形式的标量声波方程描述 .考虑到在无损耗条件下 ,地震波方程描述了地震波场这一个无穷维的哈密顿体系随时间的演化过程 ,该过程为一个单参数连续辛变换 ,因而可以在其哈密顿形式表述下导出其辛格式 .与显式辛算法相比 ,隐式辛格式对应的隐式辛几何算法具有无条件稳定的特点 ,可以允许较大的计算步长 .但是由于隐式算法不可避免地面临高阶矩阵的求逆 ,其每一步的计算速度较慢 .为实现矩阵快速求逆 ,文中采用了螺旋边界条件下谱因式分解的方法 .在螺旋边界条件下 ,需要求逆的矩阵化为带状矩阵 ,而且其各列非零元素的位置和大小具有非常好的相似性 ,因而可以采用谱因式分解的方法实现快速LU分解 .文中采用二阶精度的隐式蛙跳辛格式和谱因式分解方法 ,计算了常速度、层状介质和Marmousi模型中的波场 .计算表明 ,隐式辛算法不失为波场计算的一种好方法
Based on the common equation of seismic wave field and its Hamiltonian formalism, an implicit symplectic method (ISM) is proposed for seismic wave modeling in this paper. Combined with helix boundary conditions, the spectral factorization method is applied to give a fast way for the ISM modeling. The propagation of seismic wave can be described by the acoustic equation of a single scalar quantity in the constant, complex isotropy and anisotropy media. With the assumption of no dissipation, the scalar equation gives the symplectic evolution of seismic wave field, which is an infinite dimensional Hamiltonian system. Therefore, its Hamiltionian formalism and symplectic schemes can be given. Compared with explicit symplectic method (ESM), ISM is unconditionally stable method. Consequently, a long time step and a long time calculation are possible for ISM, which may be impossible for ESM in the same conditions. Unfortunately, it is unavoidably to calculate the inverse of large matrix in each step of ISM, which is time consuming process. Fortunately, however, the matrix has striped shape with a helical boundary conditions. And more, the nonzero units in each row of the matrix not only are located in the same relative positions but have similar values. Therefore, spectral factorization method can be used to get the inversion of the matrix in a very quick way. Based on spectral factorization method, an ISM scheme with second order in time is applied to get the wave fields in constant velocity, layer velocity and Marmousi velocity models. As the results indicates, with the help of spectral factorization, ISM is a quite good method in seismic wave modeling.
Based on the common equation of seismic wave field and its Hamiltonian formalism, an implicit symplectic method (ISM) is proposed for seismic wave modeling in this paper. Combined with helix boundary conditions, the spectral factorization method is applied to give a fast way for the ISM modeling. The propagation of seismic wave can be described by the acoustic equation of a single scalar quantity in the constant, complex isotropy and anisotropy media. With the assumption of no dissipation, the scalar equation gives the symplectic evolution of seismic wave field, which is an infinite dimensional Hamiltonian system. Therefore, its Hamiltionian formalism and symplectic schemes can be given. Compared with explicit symplectic method (ESM), ISM is unconditionally stable method. Consequently, a long time step and a long time calculation are possible for ISM, which may be impossible for ESM in the same conditions. Unfortunately, it is unavoidably to calculate the inverse of large matrix in each step of ISM, which is time consuming process. Fortunately, however, the matrix has striped shape with a helical boundary conditions. And more, the nonzero units in each row of the matrix not only are located in the same relative positions but have similar values. Therefore, spectral factorization method can be used to get the inversion of the matrix in a very quick way. Based on spectral factorization method, an ISM scheme with second order in time is applied to get the wave fields in constant velocity, layer velocity and Marmousi velocity models. As the results indicates, with the help of spectral factorization, ISM is a quite good method in seismic wave modeling.
引文
[ 1 ] 杨孔庆 ,罗明秋 ,李幼铭等 .地震波传播的微分几何描述 .地球物理学进展 ,1 998,1 3( 3) :32— 37.WANGKong Qing,LUOMing Qiu ,LIYou Ming,Differentialgeometricdescribingonseismicwavepropagation .ProgressinGeophysics (inChinese) ,1 998,1 3( 3) :32— 37.
[2 ] 罗明秋 ,刘 洪 ,李幼铭等 .地震波场的几何描述 .现代地质 ,1 999,1 3(增刊 ) :72— 77.LUOMing Qiu,LIUHong ,LIYou Ming ,etal.Geometricdescriptionofseismicwavefield ,Geoscience———Jour nalofgraduateschool (inChinese) ,1 999,1 3(Suppl.) :72— 77.
[3] BunemanO .AdvantagesofHamiltonFormulationinComputerSimulations.In :Tabor,ed.,MathematicalMethodsinHydrodynamicsandIntegrabilityofDynamicalSystems.Amer.Inst.Phys.USA ,1 981 .
[4 ] FengK .Ondifferenceschemesandsymplecticgeometry .In :FengK .ed,Proc 1 984 ,Beijing ,SympDiffGeometryandDiffEquations.Beijing :SciencePress,1 985,4 2 — 58.
[5] FengK ,QinM .Z .Hamiltonalgorithmsandacomparativenumericalstudy.ComputPhysComm ,1 991 ,65:1 73—1 87.
[6] FengK ,QinMZ .ThesymplecticmethodsfortheComputationofHamiltonequation .In :ZhuYL ,GuoBen yu ,ed .ProcConfonNumericalMethodsforPDE’s.Berlin :Springer,1 987.
[7] QinM .Z ,ZhangM .Q .Muti stagesymplecticschemesoftwokindofHamiltoniansystermsofwaveequations.ComputerMathApplic,1 990 ,1 9( 1 0 ) :51— 62 .
[8] 罗明秋 ,刘 洪 ,李幼铭 .地震波传播的哈密顿表述及辛几何算法 .地球物理学报 ,2 0 0 1 ,4 4( 1 ) :1 2 0— 1 2 8LUOMing Qiu ,LIUHong,LIYou Ming .Hamiltoniandescriptionofseismicwavepropagationandsymplecticmethods.ChineseJ.Geophys.(inChinese) ,2 0 0 1 ,4 4( 1 ) :1 2 0— 1 2 8.
[9] 刘 洪 ,罗明秋 ,李幼铭 .共反射点轨迹的哈密顿方法 .地球物理学报 ,1 999,4 2 ( 5) :685— 694 .LIUHong,LUOMing Qiu,LIYou Ming.Hamiltonmethodindescribingcommonreflectionpointtrajectory .Chi neseJ.Geophys.(inChinese) ,1 999,4 2 ( 5) :685— 694 .
[1 0 ] LuoMingqiu ,LiYouming ,LiuHong .TheSymplecticGeometricDescriptionandAlgorithmofSeismicWavePropa gation ,SEG 69thInternationalMeeting ,ExpandedAbstracts,1 999,1 852 — 1 854.
[1 1 ] 罗明秋 .反射地震运动学的哈密顿方法 [硕士论文 ].兰州 :兰州大学 ,1 998.LUOMing Qiu .Hamiltonmethodinkinematicsreflectiveseismology [M .Sc.thesis](inChinese) ,Lanzhou :LanzhouUniversity ,1 998.
[1 2 ] KummerB ,BehleA ,DorauF .Hybridmodelingofelastic wavepropagationintwo dimensionallaterallyinhomoge neousmedia.Geophysics,1 987,52 ( 2 ) :765— 771 .
[1 3] KuoJT ,YengYC ,PecholcsP ,etal.Anoteontheinfluenceoftheelasticityofwavepropagationinanacoustic/elasticmedia.Comp .andMathAppl.,1 985,1 1 :887— 896.
[1 4 ] MuftiIR .seismicmodelingintheimplicitmode.Geophysicalprospecting,33:61 9— 656.
[1 5] 刘 洪 ,罗明秋 ,李幼铭 .横向线性变速介质中的共炮检距波场延拓 .地球物理学报 ,1 999,4 2 ( 6) :80 9— 81 7.LIUHong ,LUOMing Qiu,LIYou Ming.Commonoffsetcontinuationinlaterallinearlyvariedvelocitymedia.ChineseJ.Geophys.(inChinese) ,1 999,4 2 ( 6) :80 9— 81 7.
[1 6] LuoMingqiu ,LiYouming ,LiuHong,AmplitudePreservingOffsetContinuationinLateralVariedVelocityMedia,SEG 69thInternationalMeeting .ExpandedAbstracts,1 999,1 933— 1 936.
[1 7] ClaerboutJ.Multidimensionalrecursivefiltersviaahelixwithapplicationtovelocityestimationand 3 Dmigration ,SEG 68thInternationalMeeting.ExpandedAbstracts,1 998,1 995— 1 998.
[1 8] RickettJ.Implicit3 Ddepthmigrationbywavefieldextrapolationwithhelicalboundarycondition ,SEG 68thInterna tionalMeeting.ExpandedAbstracts,1 998,1 1 2 4 — 1 1 2 7.
[1 9] ClaerboutJ.Multidimensionalrecursivefiltersviaahelix .Geophysics,1 998,63:1 532— 1 541 .
[2 0 ] LuoMingqiu ,GaoHongwei,LiuHong,etal.AcomparisonontheLUdecompositionofLaplacianmatrixwithspec tralandgraphicmethods,SEG 70 thInternationalMeeting .ExpandedAbstracts,2 0 0 0 .
[2 ] 罗明秋 ,刘 洪 ,李幼铭等 .地震波场的几何描述 .现代地质 ,1 999,1 3(增刊 ) :72— 77.LUOMing Qiu,LIUHong ,LIYou Ming ,etal.Geometricdescriptionofseismicwavefield ,Geoscience———Jour nalofgraduateschool (inChinese) ,1 999,1 3(Suppl.) :72— 77.
[3] BunemanO .AdvantagesofHamiltonFormulationinComputerSimulations.In :Tabor,ed.,MathematicalMethodsinHydrodynamicsandIntegrabilityofDynamicalSystems.Amer.Inst.Phys.USA ,1 981 .
[4 ] FengK .Ondifferenceschemesandsymplecticgeometry .In :FengK .ed,Proc 1 984 ,Beijing ,SympDiffGeometryandDiffEquations.Beijing :SciencePress,1 985,4 2 — 58.
[5] FengK ,QinM .Z .Hamiltonalgorithmsandacomparativenumericalstudy.ComputPhysComm ,1 991 ,65:1 73—1 87.
[6] FengK ,QinMZ .ThesymplecticmethodsfortheComputationofHamiltonequation .In :ZhuYL ,GuoBen yu ,ed .ProcConfonNumericalMethodsforPDE’s.Berlin :Springer,1 987.
[7] QinM .Z ,ZhangM .Q .Muti stagesymplecticschemesoftwokindofHamiltoniansystermsofwaveequations.ComputerMathApplic,1 990 ,1 9( 1 0 ) :51— 62 .
[8] 罗明秋 ,刘 洪 ,李幼铭 .地震波传播的哈密顿表述及辛几何算法 .地球物理学报 ,2 0 0 1 ,4 4( 1 ) :1 2 0— 1 2 8LUOMing Qiu ,LIUHong,LIYou Ming .Hamiltoniandescriptionofseismicwavepropagationandsymplecticmethods.ChineseJ.Geophys.(inChinese) ,2 0 0 1 ,4 4( 1 ) :1 2 0— 1 2 8.
[9] 刘 洪 ,罗明秋 ,李幼铭 .共反射点轨迹的哈密顿方法 .地球物理学报 ,1 999,4 2 ( 5) :685— 694 .LIUHong,LUOMing Qiu,LIYou Ming.Hamiltonmethodindescribingcommonreflectionpointtrajectory .Chi neseJ.Geophys.(inChinese) ,1 999,4 2 ( 5) :685— 694 .
[1 0 ] LuoMingqiu ,LiYouming ,LiuHong .TheSymplecticGeometricDescriptionandAlgorithmofSeismicWavePropa gation ,SEG 69thInternationalMeeting ,ExpandedAbstracts,1 999,1 852 — 1 854.
[1 1 ] 罗明秋 .反射地震运动学的哈密顿方法 [硕士论文 ].兰州 :兰州大学 ,1 998.LUOMing Qiu .Hamiltonmethodinkinematicsreflectiveseismology [M .Sc.thesis](inChinese) ,Lanzhou :LanzhouUniversity ,1 998.
[1 2 ] KummerB ,BehleA ,DorauF .Hybridmodelingofelastic wavepropagationintwo dimensionallaterallyinhomoge neousmedia.Geophysics,1 987,52 ( 2 ) :765— 771 .
[1 3] KuoJT ,YengYC ,PecholcsP ,etal.Anoteontheinfluenceoftheelasticityofwavepropagationinanacoustic/elasticmedia.Comp .andMathAppl.,1 985,1 1 :887— 896.
[1 4 ] MuftiIR .seismicmodelingintheimplicitmode.Geophysicalprospecting,33:61 9— 656.
[1 5] 刘 洪 ,罗明秋 ,李幼铭 .横向线性变速介质中的共炮检距波场延拓 .地球物理学报 ,1 999,4 2 ( 6) :80 9— 81 7.LIUHong ,LUOMing Qiu,LIYou Ming.Commonoffsetcontinuationinlaterallinearlyvariedvelocitymedia.ChineseJ.Geophys.(inChinese) ,1 999,4 2 ( 6) :80 9— 81 7.
[1 6] LuoMingqiu ,LiYouming ,LiuHong,AmplitudePreservingOffsetContinuationinLateralVariedVelocityMedia,SEG 69thInternationalMeeting .ExpandedAbstracts,1 999,1 933— 1 936.
[1 7] ClaerboutJ.Multidimensionalrecursivefiltersviaahelixwithapplicationtovelocityestimationand 3 Dmigration ,SEG 68thInternationalMeeting.ExpandedAbstracts,1 998,1 995— 1 998.
[1 8] RickettJ.Implicit3 Ddepthmigrationbywavefieldextrapolationwithhelicalboundarycondition ,SEG 68thInterna tionalMeeting.ExpandedAbstracts,1 998,1 1 2 4 — 1 1 2 7.
[1 9] ClaerboutJ.Multidimensionalrecursivefiltersviaahelix .Geophysics,1 998,63:1 532— 1 541 .
[2 0 ] LuoMingqiu ,GaoHongwei,LiuHong,etal.AcomparisonontheLUdecompositionofLaplacianmatrixwithspec tralandgraphicmethods,SEG 70 thInternationalMeeting .ExpandedAbstracts,2 0 0 0 .
版权所有：© 2023 中国地质图书馆 中国地质调查局地学文献中心