[發明專利]基于時空變系數加權的優化疊加方法及系統

在審

• 申請號: 201710632145.7
• 申請日: 2017-07-28
• 公開（公告）號: CN109307889A
• 公開（公告）日: 2019-02-05
• 發明（設計）人: 許璐;李博;王小青
• 申請（專利權）人: 中國石油化工股份有限公司;中國石油化工股份有限公司石油物探技術研究院
• 主分類號: G01V1/36
• 分類號: G01V1/36
• 代理公司: 北京思創畢升專利事務所 11218
• 代理人: 孫向民;廉莉莉

本發明提出了一種基于時空變系數加權的優化疊加方法及系統，該方法包括：由原始數據疊加求取參考序列；計算參考序列和疊前道集間的局部時空變相關系數；確定時空變的加權系數；將時空變的加權系數應用于對應的疊前道集；將加權濾波后的疊前道集進行疊加，得到最終疊加剖面。本發明介紹了一種利用疊前道集與參考序列集之間的時空變相關系數作為加權系數的優化疊加方法，加權系數濾波處理后的疊前道集中的噪聲干擾得到有效壓制，有效信號得到增強，提高了地震成像道集和疊加剖面的信噪比和整體質量。

技術領域

本發明屬于地震勘探資料處理技術領域，涉及高精度地震成像及地震解釋處理技術，具體地涉及基于時空變系數加權的優化疊加方法及系統。

背景技術

面對復雜地區的勘探開發，獲得正確的構造認識和油氣有利圈閉對地震解釋工作的精度要求越來越高。由于得到共反射點道集的過程是在疊前實現的，因此其信噪比一般比較低，低信噪比成像道集的直接疊加會降低成像疊加剖面的質量，困擾著傳統解釋工作的進行。不能被現有方法有效利用的噪聲干擾是造成信噪比低的原因之一，并且以隨機噪聲居多，可以采用疊前去除隨機噪聲方法來消除或壓制，如f-k域傾向濾波、f-x域預測濾波、中值濾波等，但是這些濾波方法沒有方向特性，易造成對傾斜同相軸振幅的壓制，同時也會模糊小的斷層和裂縫等不連續結構。

結構濾波是一種具有地層傾向的濾波方法，能夠用于疊前資料的噪聲壓制，并且具有保護邊緣結構的良好特性，但如何有效地對地層傾向進行估計是該方法的關鍵和難點。

發明內容

本發明介紹了一種利用疊前道集與參考序列集之間的時空變相關系數作為加權系數的優化疊加方法，加權系數濾波處理后的疊前道集中的噪聲干擾得到有效壓制，有效信號得到增強，提高了地震成像道集和疊加剖面的信噪比和整體質量，對于推動更高精度的地震成像具有重要的現實意義。

針對地震成像道集中存在的噪聲干擾問題，本發明介紹了一種利用疊前道集與參考序列集之間的時空變相關系數作為加權系數的優化疊加方法，加權系數濾波處理后的疊前道集中的噪聲干擾得到有效壓制，有效信號得到增強，提高了地震成像道集和疊加剖面的信噪比和整體質量。

根據本發明的一個方面，提供一種基于時空變系數加權的優化疊加方法，該方法包括：

由原始數據疊加求取參考序列；

計算參考序列和疊前道集間的局部時空變相關系數；

確定時空變的加權系數；

將時空變的加權系數應用于對應的疊前道集；

將加權濾波后的疊前道集進行疊加，得到最終疊加剖面。

進一步地，所述參考序列的求取公式為：

式中，N為疊加道數，j為采樣點數，α_i,j(t)為疊前道集數據。

進一步地，參考序列和疊前道集間的局部時空變相關系數表示為：

c₁＝[λ²I+S(A^TA-λ²I)]^-1SA^Tb,

c₂＝[λ²I+S(B^TB-λ²I)]^-1SB^Ta.

式中，A和B分別表示由a和b元素組成的對角矩陣算子，S為整形算子，a和b為兩個離散信號的矢量表示，λ為正則化控制算子，c₁和c₂的積的各個分量即為兩個信號的局部時空變相關系數。

進一步地，時空變的加權系數的優化疊加公式表示為：