各向异性双相介质弹性波场褶积算法数值模拟 - 地震工程学报.pdf

第 33卷 第 4期 西 北 地 震 学 报 Vo1．33 No．4 2011年 12月 N0RTHWESTERN SEISMOLOGICALJOURNAL Dec．，2O11 各向异性双相介质弹性波场褶积算法数值模拟 宫 猛，李信 富 (中国地质大学(北京)地球物理与信息技术学院，北京 100083) 摘 要：从各向同性介质中波场数值模拟的褶积微分算子法 出发 ，推导出了各向异性双相介质 中波 场传播数值计算的褶积新算法。将常见的二阶微分 Biot波动方程用等效的一阶速度一应力双 曲 方程表示，其中未知的波场向量包括 固相和流体的速度分量和应力分量，由此对方程的时间项使用 交错 网格差分方法计算，而对空间项则采用褶积微分算法进行求解。对各 向异性双相介质在单层 介质模型和双层介质模型中的波场特征进行 了研究。研究的结果显示，在两层介质分界面上 当地 震波产生反射时能观i贝l到两类纵波和横波，并且在衰减系数大的介质里慢纵波很难见到 。 关键词 ：各向异性介质 ：弹性波 ；数值模拟 ；褶积微分算子 中图分类号 ：P315．31 文献标识码 ：A 文章编号 ：1000—0844(2011)04—0313—06 NumericalSimulationofElasticW avePropagation inHeter0gene0us Two—phaseM ediawithConvolutionalDifferentiatorM ethod GONG Meng，L1Xin—fu (Schoolo GeophysicsandInformationTechnology，ChinaUniversityofGeosciences。Beijing 100083，China) Abstract：A new numericalsimulationmethodwithconvolutionalalgorithm forelasticwavesprop— agationinheterogeneoustwo—phasemediaisderivedfrom theconvolutionaldifferentialsimulation forthepropagationinhomogeneousmedia，andthesecond—orderBiotwaveequationisexpressed asfirst—ordervelocity—stresshyperbolicequations，inwhich theunknownwavefieldvectorsare velocityand stresscomponentsin both fluidand solidphases．Thetemporalterm in thewave equationiscomputedbystager—gridfinitedifferencemethodandthespatialterm iscomputedby thenew convolutionaldifferentiator．Thewavefieldcharactersofheterogeneoustwo—phasemedia arestudiedinbothsingleandtwo-layeredmodels．Theresultsshow thatonthesubsurfaceinthe two—layeredmodel，twotypesoffastP—waveandonetypeofS—wavecanbeobservedwhenseis— micwaveisreflected。meanwhileinthemediawith1argeattenuationcoefficienttheslow P—wave ishardtobeobserved． Keywords：Heterogeneousmedia；Elasticwave；Numericalsimulation；Convolutionaldifferen