改进BISQ模型地震波场数值模拟中的边界处理.pdfVIP

第22卷第3期 中国有色金属学报 2012年 3月 、，0lJ22NO．3 TheChineseJournalofNonferrousM etals M ar．20l2 文章编号：1004—0609(2012)03-0976-09 改进BISQ模型地震波场数值模拟中的边界处理 周竹生 一，王 红 ， (1．中南大学 有色金属成矿预测教育部重点实验室，长沙410083 2．中南大学 地球科学与信息物理学院，长沙410083) 摘 要：从改进BISQ模型双相介质所对应的一阶速度一应力运动方程出发，构建2×2Ⅳ阶交错网格有限差分模 拟算法，为了尽可能地减小或消除数值模拟中由人工边界引起的虚假反射，建立完全匹配层(PML)吸收边界的2× 2N阶交错网格有限差分算法。详细地讨论了PML吸收边界条件的构建及其有限差分算法的实现。通过MATLAB 编程进行波场模拟，将加入 PML吸收边界、常规指数衰减吸收边界及未加吸收边界的3种数值模拟结果进行对 比，论证PML吸收边界能十分有效地吸收边界反射。 关键词：改进BISQ模型；双相介质；有限差分；交错网格；完全匹配层(PML) 中图分类号：P631．4 文献标志码：A Boundarytreatmentfornumericalsimulationof seismicwavesbasedonreformulatedBISQmodel ZHOU Zhu．sheng一．WANG Hong， (1．KeyLaboratoryofMetallogenicPredictionofNonferrousMetals，Ministry ofEducation， CentralSouthUniversity,Changsha410083，China； 2．SchoolofGeosciencesandInfo—Physics，CentralSouthUniversity，Changsha410083，China) Abstract：Basedonfirst—ordervelocity—stressequationsofthereformulatedBISQmodelindouble—phasemedia，the simulationalgorithm offinitedifferenceof2-ordertimeand2N—orderspacestaggered—gridwasbuilt．Meanwhile，in ordertominimum effectscausedbytheartificialboundaryinthenumericalsimulation，theconstructionoftheperfectly matchedlayer(PML)absorbingboundaryconditionandtherealizationofthefinite—differencealgorithmwerediscussed indetail．ThearithmeticwasrealizedwithMATLAB．Comparedwiththeconventionaldecayingexponentialabsorbing boun dary andthenon—absorbingboundary，thePMLboun dary canmoreeffectively attenuatereflections，which is supposedbythewavefieldmodeling． Keywords：refomrulatedBISQ model；double—phasemedium；finite—differencemethod