[發(fā)明專利]基于黏聲方程的正演模擬方法及系統(tǒng)

在審

• 申請?zhí)枺?/td>: 201910913281.2
• 申請日: 2019-09-25
• 公開（公告）號: CN112558146A
• 公開（公告）日: 2021-03-26
• 發(fā)明（設(shè)計）人: 許璐;李博;韓冬
• 申請（專利權(quán)）人: 中國石油化工股份有限公司;中國石油化工股份有限公司石油物探技術(shù)研究院
• 主分類號: G01V1/28
• 分類號: G01V1/28
• 代理公司: 北京思創(chuàng)畢升專利事務(wù)所 11218
• 代理人: 孫向民;廉莉莉

1.一種基于黏聲方程的正演模擬方法，其特征在于，包括：

步驟1：在地震數(shù)據(jù)的有效頻帶范圍內(nèi)，給定L個應(yīng)變擬合參數(shù)τ_σl的初始值，其中L表示擬合標準線性體的個數(shù)，l＝1…L；

步驟2：在所述有效頻帶范圍內(nèi)，基于所述L個應(yīng)變擬合參數(shù)τ_σl的初始值計算當前計算頻率下的積分函數(shù)值；

步驟3：基于所述積分函數(shù)值，根據(jù)擬合參數(shù)表達式獲得當前Q值對應(yīng)的應(yīng)力擬合參數(shù)τ；

步驟4：針對Q值場數(shù)據(jù)體內(nèi)的所有Q值點，重復(fù)步驟1至步驟3，得到擬合參數(shù)數(shù)據(jù)體；

步驟5：根據(jù)所述擬合參數(shù)數(shù)據(jù)體和時間空間域控制方程，獲得黏聲介質(zhì)波場傳播算子，進行波場的正演模擬。

2.根據(jù)權(quán)利要求1所述的基于黏聲方程的正演模擬方法，其特征在于，根據(jù)以下公式(4)計算所述積分函數(shù)值：

其中，F(xiàn)表示積分函數(shù)，w為當前計算頻率。

3.根據(jù)權(quán)利要求2所述的基于黏聲方程的正演模擬方法，其特征在于，所述擬合參數(shù)表達式為：

其中，w_a、w_b為有效頻帶范圍的邊界值，Q₀表示當前Q值。

4.根據(jù)權(quán)利要求3所述的基于黏聲方程的正演模擬方法，其特征在于，所述時間空間域控制方程為：

其中，P為波場，v為介質(zhì)速度，t為旅行時，Ψ為輔助變量，s為震源。

5.根據(jù)權(quán)利要求4所述的基于黏聲方程的正演模擬方法，其特征在于，利用高階有限差分算法在時間空間域求解所述時間空間域控制方程，獲得所述黏聲介質(zhì)波場傳播算子。

6.一種基于黏聲方程的正演模擬系統(tǒng)，其特征在于，所述系統(tǒng)包括：

存儲器，存儲有計算機可執(zhí)行指令；

處理器，所述處理器運行所述存儲器中的計算機可執(zhí)行指令，執(zhí)行以下步驟：

步驟1：在地震數(shù)據(jù)的有效頻帶范圍內(nèi)，給定L個應(yīng)變擬合參數(shù)τ_σl的初始值，其中L表示擬合標準線性體的個數(shù)，l＝1…L；

步驟2：在所述有效頻帶范圍內(nèi)，基于所述L個應(yīng)變擬合參數(shù)τ_σl的初始值計算當前計算頻率下的積分函數(shù)值；

步驟3：基于所述積分函數(shù)值，根據(jù)擬合參數(shù)表達式獲得當前Q值對應(yīng)的應(yīng)力擬合參數(shù)τ；

步驟4：針對Q值場數(shù)據(jù)體內(nèi)的所有Q值點，重復(fù)步驟1至步驟3，得到擬合參數(shù)數(shù)據(jù)體；

步驟5：根據(jù)所述擬合參數(shù)數(shù)據(jù)體和時間空間域控制方程，獲得黏聲介質(zhì)波場傳播算子，進行波場的正演模擬。

7.根據(jù)權(quán)利要求6所述的基于黏聲方程的正演模擬系統(tǒng)，其特征在于，根據(jù)以下公式(4)計算所述積分函數(shù)值：

其中，F(xiàn)表示積分函數(shù)，w為當前計算頻率。

8.根據(jù)權(quán)利要求7所述的基于黏聲方程的正演模擬系統(tǒng)，其特征在于，所述擬合參數(shù)表達式為：

其中，w_a、w_b為有效頻帶范圍的邊界值，Q₀表示當前Q值。

9.根據(jù)權(quán)利要求8所述的基于黏聲方程的正演模擬系統(tǒng)，其特征在于，所述時間空間域控制方程為：

其中，P為波場，v為介質(zhì)速度，t為旅行時，Ψ為輔助變量，s為震源。

10.根據(jù)權(quán)利要求9所述的基于黏聲方程的正演模擬系統(tǒng)，其特征在于，利用高階有限差分算法在時間空間域求解所述時間空間域控制方程，獲得所述黏聲介質(zhì)波場傳播算子。