[發(fā)明專利]基于二階偏微分波動方程的PML吸收邊界條件的方法

有效

• 申請?zhí)枺?/td>: 201310606742.4
• 申請日: 2013-11-25
• 公開（公告）號: CN103616721A
• 公開（公告）日: 2014-03-05
• 發(fā)明（設(shè)計）人: 秦臻;宋建勇;馬曉宇;宋雪娟
• 申請（專利權(quán)）人: 中國石油天然氣股份有限公司
• 主分類號: G01V1/28
• 分類號: G01V1/28
• 代理公司: 北京三友知識產(chǎn)權(quán)代理有限公司 11127
• 代理人: 賈磊

1.一種基于二階偏微分波動方程的PML吸收邊界條件的方法，其特征在于，包括以下步驟：

加載某一采樣時刻的震源的縱波波場；

計算所述采樣時刻的縱波波場在三維空間某一方向的二階偏導(dǎo)數(shù)，并根據(jù)所述方向的二階偏導(dǎo)數(shù)對所述方向的完全匹配層區(qū)域內(nèi)的縱波波場的時間一階導(dǎo)數(shù)波場進(jìn)行吸收衰減；

重復(fù)上一步驟，分別對所述采樣時刻的縱波波場在三維空間另外兩個方向的完全匹配層區(qū)域內(nèi)的縱波波場的時間一階導(dǎo)數(shù)波場進(jìn)行吸收衰減；

重復(fù)以上步驟，直至完成所述縱波波場所有采樣時刻的處理。

2.根據(jù)權(quán)利要求1所述的方法，其特征在于，在所述在整個三維數(shù)值模擬計算區(qū)域?qū)λ隹v波波場進(jìn)行更新之后，重復(fù)以上步驟，直至完成所述縱波波場所有采樣時刻的處理之前，還包括：

在整個三維數(shù)值模擬計算區(qū)域?qū)λ隹v波波場進(jìn)行更新。

3.根據(jù)權(quán)利要求2所述的方法，其特征在于，在所述加載某一采樣時刻的震源的縱波波場之前，還包括：

預(yù)先將縱波波場的三維速度模型在三維空間上進(jìn)行網(wǎng)格離散化，并設(shè)置波動方程數(shù)值模擬所需參數(shù)；

預(yù)先將經(jīng)過網(wǎng)格離散化處理后的三維速度模型的整個三維空間數(shù)值模擬區(qū)域劃分為主區(qū)域和完全匹配層區(qū)域，并在所述完全匹配層區(qū)域的每個方向各引入兩個完全匹配層輔助波場，所述主區(qū)域?yàn)榭v波波場傳播時能量不被吸收衰減的區(qū)域；

預(yù)先對所述完全匹配層區(qū)域進(jìn)行參數(shù)設(shè)置。

4.根據(jù)權(quán)利要求3所述的方法，其特征在于，所述波動方程數(shù)值模擬所需參數(shù)包括采樣網(wǎng)格數(shù)、三維空間方向的采樣步長、三維空間方向的采樣點(diǎn)數(shù)、時間采樣點(diǎn)數(shù)、時間采樣步長、地震子波類型和主頻；所述主區(qū)域的波動方程為：

∂Q∂t=υ2(∂2P∂x2+∂2P∂y2+∂2P∂z2)+f(t),]]>

其中，P為縱波波場，Q為引入的所述縱波波場P對時間變量t的一階導(dǎo)數(shù)波場，υ為聲波波速，x、y、z為三維空間方向，f(t)為地震子波。