滑动最小二乘法求解地震波波动方程-中国科技论文在线.PDFVIP

中国科技论文在线 第 20卷 第 4期 地 球 物 理 学 进 展 VoI_2O No．4 2005年 12月 (页码 ：920～924) PR()GRESS 1N GE0PHYS1CS Dec． 2005 滑动最小二乘法求解地震波波动方程 贾晓峰， 胡天跃 (北京大学地球与空间科学学院，北京 100871) 摘 要 本文首次将 固体力学领域 中的滑动最小二乘拟合法用于求解地震波波动方程．与有限差分法相 比，该方法 所采用的拟合思想使得解在全空间域更加连续，同时仍然保持 高精度的特性 ；与有限元法相 比，该方法绕过变分原 理，经过算法的有限差分化之后 ，其所需的计算量及数据存储量均可大大降低．薄膜震动 的数值算例表 明了该方法的 上述特 点．最后 ，通过对几个模 型的合成地震模拟试验 ，进一步证 明了滑动最小二乘法用于波动方程数值模拟 的可行 性 ． 关键词 滑动最小二乘，有 限差分 ，有限元 ，波动方程 ，模拟 ，成本 中图分类号 P315，P641 文献标识码 A 文章编号 1004—2903(2005)04-0920—05 Solvingseismicwaveequationbymovingleastsquares(M LS)method JIA Xiao—feng，HU Tian—yue (SchoolofEarthandSpaceSciences，PekingUniversity，Beijing100871，China) Abstract Inthispaper，themovingleastsquares (MLS)fittingcriterioncommonlyusedinsolidmechanicsisdis— cussedforsolving seismicwaveequations．Thefittingin MLS insteadofinterpolation infinitedifferencemethod (FDM)makesthesolutionmorecontinuousin thewholespacedomain withoutany lossofaccuracy．Ontheother hand，althoughtheabsenceofvariationalprinciplewhichisnecessaryinfiniteelementmethod(FEM)savelittlecom— putationalvolume，thecostcouldfurtherbecutdowngreatlybyadifferentiatedmethod．Anexampleofvibrantfilm showsthemeritsoftheMISmethod．Furthermore，somesyntheticseismogramsobtainedbyM ISmethodproveits effectivityforwaveequationseismicmodeling． Keywords movingleastsquares，finitedifference，finiteelement，waveequation，modeling，cost 算性能要求较高．因此在这两种方法不断被改进 的 0 引 言 同时，也需要对各种新方法进行有效的尝试． 地震波模拟与成像是地震学特别是勘探地震学 滑动最小二乘拟合方法 (MLS)最初是在