用于波场成像的谱法LU分解

A-Lin

用于波场成像的谱法LU分解

9 C% c5 B* F. a! m0 b4 a. _, I9 q
) Y( v+ |; C3 }; K& E7 D
4 \# f$ N# \* X* M7 Y, x6 H" ~8 H地震波场模拟和偏移成像等有限差分隐格式算法中的重要环节, 是实现亥姆霍兹算子表示矩阵 H 的快速求逆运算. 在螺旋边界条件下, H 具有 Toeplitz 结构的正定厄密矩阵, 其快速求逆可由谱法 LU 分解实现.本文分析了谱法 LU 分解对提高计算速度的原理及特点, 并着重讨论了在不同类型的介质模型中, 采用谱法分解矩阵 H时带来的数值误差、误差的分布及其对波场计算的影响. 研究结果表明, 对均匀介质而言, 矩阵 H 各列具有相同的非零元素分布, 谱法 LU 分解的误差在吸收边界条件下, 不影响波场模拟和成像计算; 但对于非均匀介质模型, 矩阵H 各列具有不同的非零元素分布, 谱法 LU 分解的误差随介质不均匀性程度的增大而增大, 势必影响非均匀介质中波场计算. 在波场模拟和成像等有限差分隐格式算法中, 采用谱法 LU 分解完成矩阵求逆时, 必须考虑到并尽量减少该方法的误差对波场计算的影响。

游客，如果您要查看本帖隐藏内容请回复

fanichicl

研究研究，谢谢楼主分享。