[發(fā)明專利]基于貝葉斯和級(jí)數(shù)反演理論的AVO反演方法及系統(tǒng)

有效

• 申請(qǐng)?zhí)枺?/td>: 201910571024.5
• 申請(qǐng)日: 2019-06-28
• 公開(kāi)（公告）號(hào): CN112147682B
• 公開(kāi)（公告）日: 2022-08-02
• 發(fā)明（設(shè)計(jì)）人: 時(shí)磊;劉俊州;韓磊;王震宇;溫立峰;薛明喜
• 申請(qǐng)（專利權(quán)）人: 中國(guó)石油化工股份有限公司;中國(guó)石油化工股份有限公司石油勘探開(kāi)發(fā)研究院
• 主分類號(hào): G01V1/30
• 分類號(hào): G01V1/30
• 代理公司: 北京思創(chuàng)畢升專利事務(wù)所 11218
• 代理人: 孫向民;廉莉莉

1.一種基于貝葉斯和級(jí)數(shù)反演理論的AVO反演方法，其特征在于，包括：

根據(jù)Zoeppritz方程和Gassmann理論，獲得基于入射角近似的流體項(xiàng)高階AVO表達(dá)式；

根據(jù)級(jí)數(shù)理論，計(jì)算所述流體項(xiàng)高階AVO表達(dá)式的參數(shù)，獲得流體項(xiàng)高階AVO公式；

根據(jù)貝葉斯理論獲取反演目標(biāo)函數(shù)；

根據(jù)所述流體項(xiàng)高階AVO公式與所述反演目標(biāo)函數(shù)，構(gòu)建疊前縱波高階AVO反演公式；

其中，所述流體項(xiàng)高階AVO表達(dá)式為：

其中，

其中，

當(dāng)時(shí)，為一階線性AVO表達(dá)式；當(dāng)時(shí)，為二階非線性AVO表達(dá)式；當(dāng)時(shí)，為三階非線性AVO表達(dá)式；

其中，所述根據(jù)級(jí)數(shù)理論，計(jì)算所述流體項(xiàng)高階AVO表達(dá)式的參數(shù)，獲得流體項(xiàng)高階AVO公式包括：

獲得縱波反射系數(shù)的一階、二階、三階展開(kāi)式；

根據(jù)實(shí)際地震記錄與所述一階展開(kāi)式，計(jì)算一階參數(shù)；

根據(jù)所述一階參數(shù)與所述二階展開(kāi)式，計(jì)算二階參數(shù)；

根據(jù)所述二階參數(shù)與所述三階展開(kāi)式，計(jì)算三階參數(shù)；

根據(jù)所述一階參數(shù)、所述二階參數(shù)、所述三階參數(shù)，計(jì)算所述流體項(xiàng)高階AVO表達(dá)式的參數(shù)；

將所述流體項(xiàng)高階AVO表達(dá)式的參數(shù)代入所述流體項(xiàng)高階AVO表達(dá)式中，獲得所述流體項(xiàng)高階AVO公式。

2.根據(jù)權(quán)利要求1所述的基于貝葉斯和級(jí)數(shù)反演理論的AVO反演方法，其中，所述一階展開(kāi)式為：

3.根據(jù)權(quán)利要求1所述的基于貝葉斯和級(jí)數(shù)反演理論的AVO反演方法，其中，所述二階展開(kāi)式為：

4.根據(jù)權(quán)利要求1所述的基于貝葉斯和級(jí)數(shù)反演理論的AVO反演方法，其中，所述三階展開(kāi)式為：

5.根據(jù)權(quán)利要求1所述的基于貝葉斯和級(jí)數(shù)反演理論的AVO反演方法，其中，所述流體項(xiàng)高階AVO表達(dá)式的參數(shù)為：

6.根據(jù)權(quán)利要求1所述的基于貝葉斯和級(jí)數(shù)反演理論的AVO反演方法，其中，所述反演目標(biāo)函數(shù)為：

其中，m為反演目標(biāo)函數(shù)，d為觀測(cè)數(shù)據(jù)，G為觀測(cè)數(shù)據(jù)和模型數(shù)據(jù)之間的映射算子，G^T為G的轉(zhuǎn)置矩陣，C_d為噪音協(xié)方差矩陣，α為權(quán)重系數(shù)，Q為取決于所選擇的先驗(yàn)分布類型的正則化項(xiàng)。

7.一種基于貝葉斯和級(jí)數(shù)反演理論的AVO反演系統(tǒng)，其特征在于，該系統(tǒng)包括：

存儲(chǔ)器，存儲(chǔ)有計(jì)算機(jī)可執(zhí)行指令；

處理器，所述處理器運(yùn)行所述存儲(chǔ)器中的計(jì)算機(jī)可執(zhí)行指令，執(zhí)行以下步驟：

根據(jù)Zoeppritz方程和Gassmann理論，獲得基于入射角近似的流體項(xiàng)高階AVO表達(dá)式；

根據(jù)級(jí)數(shù)理論，計(jì)算所述流體項(xiàng)高階AVO表達(dá)式的參數(shù)，獲得流體項(xiàng)高階AVO公式；

根據(jù)貝葉斯理論獲取反演目標(biāo)函數(shù)；

根據(jù)所述流體項(xiàng)高階AVO公式與所述反演目標(biāo)函數(shù)，構(gòu)建疊前縱波高階AVO反演公式；

其中，所述流體項(xiàng)高階AVO表達(dá)式為：

其中，

其中，

當(dāng)時(shí)，為一階線性AVO表達(dá)式；當(dāng)時(shí)，為二階非線性AVO表達(dá)式；當(dāng)時(shí)，為三階非線性AVO表達(dá)式；

其中，所述根據(jù)級(jí)數(shù)理論，計(jì)算所述流體項(xiàng)高階AVO表達(dá)式的參數(shù)，獲得流體項(xiàng)高階AVO公式包括：

獲得縱波反射系數(shù)的一階、二階、三階展開(kāi)式；

根據(jù)實(shí)際地震記錄與所述一階展開(kāi)式，計(jì)算一階參數(shù)；

根據(jù)所述一階參數(shù)與所述二階展開(kāi)式，計(jì)算二階參數(shù)；

根據(jù)所述二階參數(shù)與所述三階展開(kāi)式，計(jì)算三階參數(shù)；

根據(jù)所述一階參數(shù)、所述二階參數(shù)、所述三階參數(shù)，計(jì)算所述流體項(xiàng)高階AVO表達(dá)式的參數(shù)；

將所述流體項(xiàng)高階AVO表達(dá)式的參數(shù)代入所述流體項(xiàng)高階AVO表達(dá)式中，獲得所述流體項(xiàng)高階AVO公式。