應用自適應變網格緊致差分的地震波數值模擬及逆時偏移

陳亮 黃建平 王自穎 慕鑫茹

摘要：緊致有限差分（ＣｏｍｐａｃｔＤｉｆｆｅｒｅｎｃｅ，ＣＤ）格式相對于中心有限差分（ＦｉｎｉｔｅＤｉｆｆｅｒｅｎｃｅ，ＦＤ）格式具有計算精度高、邊界處理簡單等優(yōu)點，其中，規(guī)則網格四階ＣＤ 格式所需計算節(jié)點數少，計算精度高，但使用該格式求解聲波方程時需求解大型稀疏對角矩陣方程組，當模型較大時，該格式計算效率較低，同時用于存儲和計算對角矩陣方程組的內存需求也較大。自適應變網格策略能夠通過優(yōu)化網格剖分，有效降低模型網格數，從而達到提高計算效率、降低內存占用的目的。文中推導了求解二維聲波方程的時間二階中心差分、空間四階ＣＤ 的離散格式，并應用于正演模擬；通過引入自適應變網格策略，推導了修正聲波方程，發(fā)展了高效、高精度的自適應變網格緊致差分逆時偏移成像方法。數值實驗表明：該方法可實現復雜構造的高效、高分辨率成像；與規(guī)則細網格剖分方法相比，在保證精度的前提下，該方法可將Ｍａｒｍｏｕｓｉ模型的成像效率提高約３３．７％。

關鍵詞：緊致差分，自適應變網格，數值模擬，逆時偏移，計算效率，計算精度

中圖分類號：Ｐ６３１文獻標識碼：Ａ doi：１０．１３８１０／ｊ．ｃｎｋｉ．ｉｓｓｎ．１０００-７２１０．２０２３．０３．０１７