[發(fā)明專利]基于歸一化互相關(guān)的Marchenko成像方法及系統(tǒng)

在審

• 申請?zhí)枺?/td>: 202011081247.2
• 申請日: 2020-10-10
• 公開（公告）號: CN114428343A
• 公開（公告）日: 2022-05-03
• 發(fā)明（設(shè)計(jì)）人: 韓冬;李博;崔月
• 申請（專利權(quán)）人: 中國石油化工股份有限公司;中國石油化工股份有限公司石油物探技術(shù)研究院
• 主分類號: G01V1/36
• 分類號: G01V1/36
• 代理公司: 北京知舟專利事務(wù)所(普通合伙) 11550
• 代理人: 郭韞

1.一種基于歸一化互相關(guān)的Marchenko成像方法，其特征在于：所述方法以初始格林函數(shù)和校正后的地震數(shù)據(jù)作為輸入，以迭代反演的方法對Marchenko方程進(jìn)行數(shù)值求解得到上行格林函數(shù)和下行格林函數(shù)，并對所得到的上行格林函數(shù)和初始格林函數(shù)進(jìn)行歸一化互相關(guān)處理得到Marchenko成像結(jié)果。

2.根據(jù)權(quán)利要求1所述的基于歸一化互相關(guān)的Marchenko成像方法，其特征在于：所述方法包括：

步驟一：獲取初始聚焦函數(shù)；

步驟二：多次迭代獲得Marchenko聚焦函數(shù)；

步驟三：根據(jù)所述Marchenko聚焦函數(shù)獲得上行格林函數(shù)和下行格林函數(shù)；

步驟四：進(jìn)行歸一化互相關(guān)得到Marchenko成像結(jié)果。

3.根據(jù)權(quán)利要求2所述的基于歸一化互相關(guān)的Marchenko成像方法，其特征在于：所述步驟一的操作包括：

利用下式計(jì)算獲得初始聚焦函數(shù)：

G_d(x′₀，x_i，-t)為初始格林函數(shù)的逆時(shí)變換。

4.根據(jù)權(quán)利要求3所述的基于歸一化互相關(guān)的Marchenko成像方法，其特征在于：所述步驟二的操作包括：

利用下式計(jì)算獲得經(jīng)過多次迭代后的Marchenko上行聚焦函數(shù)及Marchenko下行聚焦函數(shù)

其中，R(x₀,x'₀,t)為校正后的地震數(shù)據(jù)，是在x′₀點(diǎn)處激發(fā)x₀處接收得到的地表反射響應(yīng)；

為初始上行聚焦函數(shù)；

第k-1次迭代所得的上行聚焦函數(shù)；

為第k次迭代所得的逆時(shí)Marchenko尾波函數(shù)；

為第k次迭代所得的Marchenko尾波函數(shù)；

為第k次迭代所得的下行聚焦函數(shù)；

為第k次迭代所得的上行聚焦函數(shù)；

θ為窗函數(shù)，表達(dá)式如下：

x_i為地下聚焦點(diǎn)空間位置；

x₀為地表震源點(diǎn)的空間位置。

5.根據(jù)權(quán)利要求4所述的基于歸一化互相關(guān)的Marchenko成像方法，其特征在于：所述步驟三的操作包括：

利用下式計(jì)算獲得上行格林函數(shù)和下行格林函數(shù)：

上式中，G^-(x_i,x₀,t)和G⁺(x_i,x₀,t)分別表示地表x₀點(diǎn)處激發(fā)并在地下x_i接收得到的上行格林函數(shù)和下行格林函數(shù)，和分別表示在地下x_i處聚焦的上行聚焦函數(shù)和下行聚焦函數(shù)，其值分別采用和的值。

6.根據(jù)權(quán)利要求2所述的基于歸一化互相關(guān)的Marchenko成像方法，其特征在于：所述步驟四的操作包括：

利用下式對初始格林函數(shù)和上行格林函數(shù)進(jìn)行歸一化互相關(guān)得到Marchenko成像結(jié)果：

7.根據(jù)權(quán)利要求2所述的基于歸一化互相關(guān)的Marchenko成像方法，其特征在于：

在所述步驟一之前還包括：數(shù)據(jù)準(zhǔn)備步驟；

所述數(shù)據(jù)準(zhǔn)備步驟的操作包括：

對地震數(shù)據(jù)進(jìn)行校正獲得校正后的地震數(shù)據(jù)；

以背景速度v₀作為輸入，求解程函方程得到走時(shí)場t，提取出初始時(shí)間t_d，計(jì)算獲得初始格林函數(shù)。