[發明專利]利用陸上地震資料進行波形層析的速度建模方法

有效

• 申請號: 201710507204.8
• 申請日: 2017-06-28
• 公開（公告）號: CN109143350B
• 公開（公告）日: 2021-02-19
• 發明（設計）人: 王杰;邵文潮;劉定進;胡光輝
• 申請（專利權）人: 中國石油化工股份有限公司;中國石油化工股份有限公司石油物探技術研究院
• 主分類號: G01V1/30
• 分類號: G01V1/30
• 代理公司: 北京思創畢升專利事務所 11218
• 代理人: 孫向民;廉莉莉

1.一種利用陸上地震資料進行波形層析的速度建模方法，其特征在于，該方法包括：

步驟1，基于已知的地質資料，把工區劃分成不同區塊；

步驟2，在互相關波形層析框架下構建陸上資料速度建模梯度場；

步驟3，基于不同區塊和陸上資料速度建模梯度場，進行分區塊波數域濾波以及整個區域高斯濾波；

步驟4，更新層析模型，迭代次數達到預設值則終止迭代得到最終層析結果，若不滿足轉到步驟2，此時陸上資料速度建模梯度場由更新后的層析模型求解；

其中，把工區劃分成不同區塊包括：

基于深度域偏移剖面構建地質構造傾角；

基于深度域偏移剖面構建主要層位信息；

根據構造傾角及層位信息，把工區分為不同的區塊，每個區塊對應不同的主要反射層構造傾角；

其中，采用自動追蹤層位技術提取深度域偏移剖面主要反射層層位信息，或者采用手動解釋方式來獲取層位信息；

其中，步驟3包括：

將不同區塊映射到梯度場，把梯度場分為不同的區塊，針對不同區塊采用相應的波數域濾波算子進行濾波；

對不同區塊的梯度場按順序整合成完整梯度場并進行高斯濾波；

其中，對不同區塊的梯度場按順序整合成完整梯度場進行高斯平滑濾波：

d_n＝Gs_n

式中，s_n為整合梯度場，d_n表示高斯濾波梯度場，G為高斯濾波算子，分量形式為：

其中L_x和L_z為x和z方向的相關長度。

2.根據權利要求1所述的速度建模方法，其特征在于，采用局部平面波分解算法對工區已知的深度域地震偏移剖面估計構造傾角場，其中地震偏移剖面滿足的局部平面波差分方程為：

式中，P(x,z)為深度域偏移剖面場，σ(x,z)為局部構造傾角場。

3.根據權利要求1所述的速度建模方法，其特征在于，構建觀測地震資料和人工合成地震資料的歸一化互相關場；利用歸一化互相關場構建陸上資料速度建模梯度場。

4.根據權利要求3所述的速度建模方法，其特征在于，采用歸一化互相關目標泛函完成波形層析速度建模，歸一化互相關目標泛函為：

式中，τ為延遲時間，x_r為檢波器，x_s為震源，d_cal(t,x_r；x_s)為人工合成地震記錄，d_obs(t+τ,x_r；x_s)是時延觀測地震記錄，W(τ)為衰減時窗。

5.根據權利要求1所述的速度建模方法，其特征在于，在每個區塊找到最大構造傾角θ_i，利用不同區塊的最大構造傾角和各自相應的波數域濾波算子分別對不同區塊梯度噪音壓制，不同區塊采用的噪音壓制算子公式為：

式中，表示第i區塊梯度場，Pⁱ為傅立葉正變換，Tⁱ表示波數選擇算子其是θ_i的函數，(P^-1)ⁱ為傅里葉反變換，為第i區塊噪音壓制后的梯度場。

6.根據權利要求1所述的速度建模方法，其特征在于，步驟4中采用的模型更新計算式為：

m_n+1＝m_n+αd_n，

其中m_n+1為第n+1迭代次數的層析模型；m_n為第n次迭代的層析模型；α為更新步長；d_n為梯度。