[發明專利]一種基于互相關道的散射波成像方法及系統

有效

• 申請號: 201710236627.0
• 申請日: 2017-04-12
• 公開（公告）號: CN108693560B
• 公開（公告）日: 2020-04-07
• 發明（設計）人: 程磊磊;姜宇東;王躍
• 申請（專利權）人: 中國石油化工股份有限公司;中國石油化工股份有限公司石油物探技術研究院
• 主分類號: G01V1/34
• 分類號: G01V1/34
• 代理公司: 北京思創畢升專利事務所 11218
• 代理人: 孫向民;廉莉莉

1.一種基于互相關道的散射波成像方法，其特征在于，該方法包括：

輸入地震數據道集、速度模型以及初始化成像參數；

正演計算檢波點到成像點的旅行時；

將地震數據道集轉化為互相關道集；

利用互相關道集進行散射波成像；

輸出散射波成像結果；

其中，將地震數據道集轉化為互相關道集包括：

讀入地震數據道集，道集中數據道數為N；

對于第i道數據道，求取第j道數據與第i道數據道間互相關函數，其中i＝1,2,...,N-1，j＝i+1,i+2,...,N；

改變i值，重復求取第j道數據與第i道數據道間互相關函數，直到任意兩道間的互相關函數計算完畢，得到互相關道集；

其中，利用公式c_ij(t)＝∫s_i(τ)·s_j(τ+t)dτ求取第j道數據與第i道數據道間互相關函數，其中，c_ij(t)表示互相關道集，s_i(t)和s_j(t)分別表示第i和j道地震數據；

其中，利用互相關道集進行散射波成像包括：

讀入任意一道互相關數據c_ij(t)，確定生成該互相關道數據的地震數據的兩個檢波器坐標r_i和r_j；

計算成像范圍內任意成像點x到兩個檢波器間的旅行時差Δt_ij(x)；

對于任意成像點x，利用計算出的旅行時差Δt_ij(x)，在互相關數據道中提取對應時刻的采樣點值；

將采樣點值乘以振幅校正因子1/(r_i-x)(r_j-x)，輸出和疊加到成像點x處；

重復上述步驟，直到相關道集中所有數據道計算完畢。

2.根據權利要求1所述的基于互相關道的散射波成像方法，其特征在于，正演計算檢波點到成像點的旅行時包括：

根據地震數據道集的檢波點坐標和成像范圍，采用旅行時正演算法，計算檢波點到所有成像點的旅行時。

3.一種基于互相關道的散射波成像系統，其特征在于，所述系統包括：

存儲器，存儲器存儲地震數據道集、速度模型、初始化成像參數以及可執行指令；

處理器，所述處理器調用存儲器存儲的數據，執行所述可執行指令完成以下步驟：

正演計算檢波點到成像點的旅行時；

將地震數據道集轉化為互相關道集；

利用互相關道集進行散射波成像；

輸出散射波成像結果；

其中，將地震數據道集轉化為互相關道集包括：

讀入地震數據道集，道集中數據道數為N；

對于第i道數據道，求取第j道數據與第i道數據道間互相關函數，其中i＝1,2,...,N-1，j＝i+1,i+2,...,N；

改變i值，重復求取第j道數據與第i道數據道間互相關函數，直到任意兩道間的互相關函數計算完畢，得到互相關道集；

其中，利用公式c_ij(t)＝∫s_i(τ)·s_j(τ+t)dτ求取第j道數據與第i道數據道間互相關函數，其中，c_ij(t)表示互相關道集，s_i(t)和s_j(t)分別表示第i和j道地震數據；

其中，利用互相關道集進行散射波成像包括：

讀入任意一道互相關數據c_ij(t)，確定生成該互相關道數據的地震數據的兩個檢波器坐標r_i和r_j；

計算成像范圍內任意成像點x到兩個檢波器間的旅行時差Δt_ij(x)；

對于任意成像點x，利用計算出的旅行時差Δt_ij(x)，在互相關數據道中提取對應時刻的采樣點值；

將采樣點值乘以振幅校正因子1/(r_i-x)(r_j-x)，輸出和疊加到成像點x處；

重復上述步驟，直到相關道集中所有數據道計算完畢。

4.根據權利要求3所述的基于互相關道的散射波成像系統，其特征在于，正演計算檢波點到成像點的旅行時包括：

根據地震數據道集的檢波點坐標和成像范圍，采用旅行時正演算法，計算檢波點到所有成像點的旅行時。