频率域并行化弹性波有限差分方法研究.pdf

·180· 地学与赍源研究新进展 频率域并行化弹性波有限差分方法研究 殷 文1 印兴耀1 吴国忱2 梁 锫1 [1．中国石油大学(华东)地球资源与信息学院，东营257061 2．同挤大学教育部海洋地质重点实验室，上海200092] 摘 要 有限差分方法是溲场数值模捌的一个重要的方法，其在地震波正演模拟及偏移成像等方面得 到了广泛应用。但常规的有限差分法本身存在着数值频散问题．会降低波场模拟的精度与分辨率，为了克服 常规差分算子的数值频散，这里采用了25点有限差分算子，再根据最优化理论求取的优化系数，建立了频率 空同域中弹性波波动方程的差分格式；为了消除边界反射．引入了最佳匹配层．构造带有最佳匹配层的差分方 程。另外，针对有限差分算法存在计算效率低的问题，本文对有限差分进行了并行算法的分析与研究，结出了 相关算法描述。最后由弹性波波动方程和边界条件，蠢过频率域有限差分法．分别利用不同震源对弹性波在 均匀各向同性介质及层状介质中的传播过程进行了数值正演模拟。通过正演模拟，得到了单频波波场、肘同 切片和共炮点道集，为下一步的研究工作如成像、反演打下了研究基础。 关键词正演模拟；有限差分算子；最佳匹配层；并行计算 TheMethodResearchofParallelElasticWave in FiniteDifferences Domain Frequency YinWenlYin WuGuochen2Kail Xingya01 Liang ofGeo—resourcesandInformation，Chinaof 257061 E1．College UniversityPetroleum，Dongying Marine of 200092] 2．Lahorator]ofGeology．MinistryEducation，TongjiUniversity，Shanghai Finitedifferencesisan methodofwave—fieldnumericalsimulation．whichis Abstract important applied insomefields．suchasseismicwa,veforwardsimulation，excursionand widely imaging．ere．However， inherethe ofnumerical whichreducethe andresolution routinefinitedifference problem dispersion precision ofwave—fieldsimulation．Inordertoovercomethenumerical ofroutinefinite—difference dispersion operatort we the coefficientonbasis adopt25一pointoptima[finite·differenceoperator，andget optimized ofoptimization finite—differenceformatofelasticwave in domain．Weintroduce andestablish theory