贴体坐标系二维声波方程SBP有限差分法的稳定性分析.pdfVIP

第 55卷第 1期 石 油 物 探 Vo1．55。No．1 2016年 1月 GE0PHYSICALPR0SPECTING F0R PETROLEUM Jan．，2016 王颖，周辉，盛善波．贴体坐标系二维声波方程SBP有限差分法的稳定性分析[J]．石油物探，2016，55(1)：33—40 W ANG Ying，ZHOU Hui，SHENG Shanbo．StabilityanalysisofSBP finitedifferenceschemefortwo-dimensionalacoustic waveequationinboundary-conforminggrids[J]．GeophysicalProspectingforPetroleum，2016，55(1)：33—40 贴体坐标系二维声波方程 SBP有限差分法的 稳定性分析 王 颖 ，周 辉 ，盛善波 (1．中国石油大学(北京)油气资源与探测国家重点实验室，中国石油大学(北京)CNPC物探重点实验室， 中国石油大学(北京)，北京 102249；2．中国石油天然气勘探开发公司，北京 100034) 摘要：有限差分方法(FiniteDifferenceMethod，FDM)是波动方程正演数值模拟领域应用最为广泛的方法之一，然而，当模 拟区域不规则或者地表起伏不平时，规则网格有限差分法求解波动方程会产生阶梯状近似，影响模拟的精度。借助贴体网 格技术，将不规则的物理区域转换为规则的计算域，给出了贴体坐标系下的二维声波方程及其二阶精度的分部求和(Sum- marionbyParts，SBP)有限差分离散格式 ，采用Fourier谱分析方法分析了该离散格式的稳定性，得到了贴体网格二维声波 方程SBP有限差分方法的稳定性条件。数值实验结果表明：① 当时间采样间隔的选取满足稳定性条件时，贴体网格 SBP 有限差分的数值计算过程是稳定的；②与贴体网格中心差分方法相比，贴体网格SBP有限差分方法的稳定性更好。 关键词：贴体坐标；声波方程；分部求和有限差分；中心差分 ；稳定性条件 中图分类号 ：P631 文献标识码 ：A ， 文章编号i1000—1441(2016)01—0033—08 DOI：10．3969／j．issn．1000—1441．2016．01．005 StabilityanalysisofSBP finitedifferenceschemefortwo—dimensional acousticwaveequation inboundary--conform inggrids WANG Ying ，ZHOU Hui，SHENG Shanbo。 (1．StateKeyLaboratoryofPetroleumResourceandProspecting，CNPCKeyL拍ofGeophysicalProspecting，ChinaUni— versityofPetroleum，Beijing102249，China；2．China NationalOil GasExplorationandDevelopmentCorporation，Bei— jing100034，China) Abstract：Traditionally，finitedifferencemethodiSwidelyusedasafastandaccuratemethodfornumericalsimulationandmi— grationofwaveequation．However，finitedifferencemethodfacesobstacleswhentherearesurfacetopographyandirregular interfaces．Thestaircaseapproximationscausedbyregulargridsinfluencetheaccuracyoffinite