地球物理学报
地毬物理學報
지구물이학보
2003年
3期
421-427
,共7页
罗明秋%刘洪%李幼铭
囉明鞦%劉洪%李幼銘
라명추%류홍%리유명
波场成像%LU分解%谱因式分解
波場成像%LU分解%譜因式分解
파장성상%LU분해%보인식분해
地震波场模拟和偏移成像等有限差分隐格式算法中的重要环节,是实现亥姆霍兹算子表示矩阵H的快速求逆运算. 在螺旋边界条件下,H具有Toepli tz结构的正定厄密矩阵,其快速求逆可由谱法LU分解实现. 本文分析了谱法LU分解对提高计算速度的原理及特点,并着重讨论了在不同类型的介质模型中,采用谱法分解矩阵H时带来的数值误差、误差的分布及其对波场计算的影响. 研究结果表明,对均匀介质而言,矩阵H各列具有相同的非零元素分布,谱法LU分解的误差在吸收边界条件下,不影响波场模拟和成像计算;但对于非均匀介质模型,矩阵H各列具有不同的非零元素分布,谱法LU分解的误差随介质不均匀性程度的增大而增大,势必影响非均匀介质中波场计算. 在波场模拟和成像等有限差分隐格式算法中,采用谱法LU分解完成矩阵求逆时,必须考虑到并尽量减少该方法的误差对波场计算的影响.
地震波場模擬和偏移成像等有限差分隱格式算法中的重要環節,是實現亥姆霍玆算子錶示矩陣H的快速求逆運算. 在螺鏇邊界條件下,H具有Toepli tz結構的正定阨密矩陣,其快速求逆可由譜法LU分解實現. 本文分析瞭譜法LU分解對提高計算速度的原理及特點,併著重討論瞭在不同類型的介質模型中,採用譜法分解矩陣H時帶來的數值誤差、誤差的分佈及其對波場計算的影響. 研究結果錶明,對均勻介質而言,矩陣H各列具有相同的非零元素分佈,譜法LU分解的誤差在吸收邊界條件下,不影響波場模擬和成像計算;但對于非均勻介質模型,矩陣H各列具有不同的非零元素分佈,譜法LU分解的誤差隨介質不均勻性程度的增大而增大,勢必影響非均勻介質中波場計算. 在波場模擬和成像等有限差分隱格式算法中,採用譜法LU分解完成矩陣求逆時,必鬚攷慮到併儘量減少該方法的誤差對波場計算的影響.
지진파장모의화편이성상등유한차분은격식산법중적중요배절,시실현해모곽자산자표시구진H적쾌속구역운산. 재라선변계조건하,H구유Toepli tz결구적정정액밀구진,기쾌속구역가유보법LU분해실현. 본문분석료보법LU분해대제고계산속도적원리급특점,병착중토론료재불동류형적개질모형중,채용보법분해구진H시대래적수치오차、오차적분포급기대파장계산적영향. 연구결과표명,대균균개질이언,구진H각렬구유상동적비령원소분포,보법LU분해적오차재흡수변계조건하,불영향파장모의화성상계산;단대우비균균개질모형,구진H각렬구유불동적비령원소분포,보법LU분해적오차수개질불균균성정도적증대이증대,세필영향비균균개질중파장계산. 재파장모의화성상등유한차분은격식산법중,채용보법LU분해완성구진구역시,필수고필도병진량감소해방법적오차대파장계산적영향.