波场模拟中的数值频散分析与校正策略.pdf

第 2 0 卷 第 1 期 地 球 物 理 学 进 展 Vol .20 No .1 2005 年3 月( 页码:58 ～65) PROGRESS I N GEOPHYSI CS Marc h 2005 波 场 模 拟 中 的数 值 频 散 分 析 与校 正 策 略 1 , 2 1 吴国忱 , 王华忠 ( 1 . 同济大 学教 育部海 洋地质 重点 实验室 , 上海 200092 ; 2 . 石 油大 学( 华东) 地 球资源 与信息 学院, 东 营 25706 1) 摘 要 波动方程有限差分法正演模拟, 对认识地震 波传播 规律 、进行地震属 性研 究、地震 资料地 质解释 、储层评价 等, 均具有重要的理论和实际意义. 但有限差分法本身固有存在着数值频散 问题, 数值频散在 正演模拟 中是 一种严重 的干扰, 会 降低波场模拟的精度与分辨率. 针对TI 介质波场模拟的交错 网格 有限差分方 法, 本文从 空间网格 离散 、时 间网格 离散和算子近似等三个方面对其产生的数值频散进行 了分析, 并结合其他 学者的研 究成 果给 出 了TI 介质波 场模拟中压制数值频散的方法与策略: 在 已知介质频散关 系时, 对差分算 子可实施算 子校正; 通过提 高差分 方程 的阶 数来提 高波场模拟精度; 采 用流体力学中守恒式方程 的通量校 正传输 方法来压制波场模拟 中的数值频散; 在 实际正 演模拟时, 采 用交错 网格 高阶有限差分方程, 不仅在 空间上采 用 高阶差分, 而且在时 间上也要采 用 高阶差分, 否则只 在单一方向上( 空间或时间) 提 高方程 的阶数对压制数值频散也不会取得理想的效果. 关键词 有限差分法, 交错 网格, 波场模拟, 数值频散, TI 介质 中图分类号 P315 ,P631 文献标识码 A 文章编号 1004- 2903( 2005) 01-0058-08 Analysis of numerical dispersion in wave-field si mulation WU Guo-chen1 , 2 , WA NG Hua-zhong 1 ( 1 . Laborat ory of Mari ne Geol ogy under Educ ation Mi nistr y , Tongj i Uni versit y , Shanghai 200092 , Chi na ; 2 . Colle ge of Earth Resource and I nf or mation , Uni versit y of Petroleu m , Dong Yi ng 257061 , Chi na) Abstract It is well kno wn t hat fi nite- difference si mul ati on of el astic wave equation is very i mportant for us to recog- nize t he regul ar pattern i n propagation of seis mi c wave 、perfor mresearch of seis mi c attri b