[發明專利]被動源地震頻率諧振勘探方法

有效

• 申請號: 201811587566.3
• 申請日: 2018-12-25
• 公開（公告）號: CN110954943B
• 公開（公告）日: 2023-05-23
• 發明（設計）人: 薛愛民
• 申請（專利權）人: 北京派特森科技股份有限公司
• 主分類號: G01V1/28
• 分類號: G01V1/28
• 代理公司: 北京市蘭臺律師事務所 11354
• 代理人: 張峰

本發明公開了一種被動源地震頻率諧振勘探方法，包括如下步驟：步驟1，利用檢波器采集地下介質的響應信號形成地震時間序列數據；步驟2，將所述步驟1采集的數據經過傅里葉變換成頻率域數據；步驟3，對經過所述步驟2的在同一檢波點位置的數據進行頻率域疊加，形成頻率域振幅疊加數據；步驟4，將經過所述步驟3處理后的頻率域數據依照標準井參數變換成深度域數據；步驟5，對經過步驟4處理的數據進行處理得到成像數據Imagesubgt;(d)/subgt;，所述成像數據Imagesubgt;(d)/subgt;是隨深度變化的視波阻抗比率或者視波阻抗。該方法可以利用地震波諧振原理對地下介質進行空間和屬性成像。

技術領域

本發明涉及地震勘探領域，特別涉及一種被動源地震頻率諧振勘探方法。

背景技術

常規地震勘探技術是在時間域進行勘探的技術。除面波勘探利用近地表面波的速度頻散特性對地下介質進行空間和屬性成像以外，其它所有的地震勘探技術均由波場走時確定勘探目的物的空間位置和屬性特征。常規地震勘探和面波地震勘探技術已經發展的比較成熟，廣泛應用于資源環境、工程地質等地學領域。

發明內容

本發明的目的在于提供一種被動源地震頻率諧振勘探方法，可以利用地震波諧振原理對地下介質進行空間和屬性成像。

為實現上述目的，本發明提供了一種被動源地震頻率諧振勘探方法，包括如下步驟：步驟1，利用檢波器采集地下介質的響應信號形成地震時間序列數據；步驟2，將步驟1采集的數據經過傅里葉變換成頻率域數據；步驟3，對經過步驟2的在同一檢波點位置的數據進行頻率域疊加，形成頻率域振幅疊加數據Amp_(f)；步驟4，將經過步驟3處理后的頻率域數據經過速度和標準井參數校正轉換到深度域數據：

如果所述步驟1的檢波器采集的是單一水平分量數據，深度域頻譜振幅Γ_(d)＝Amp_(d)為單一分量譜振幅；如果所述步驟1的檢波器采集的是多分量數據，將兩個水平分量進行合成得到深度域頻譜振幅垂直分量單獨形成振幅譜Γ_(d)，經過標準井參數校正的振幅譜隨后轉換到深度域的數據屬性為視波阻抗比率數據，經過積分變換成為視波阻抗數據；步驟5，對經過步驟4得到的數據進行圖像可視化處理后得到視波阻抗比率或視波阻抗成像數據Image(d)，或將數據直接輸出文本型數據以便于應用其它可視化軟件進行圖像處理。

優選地，對地震時間序列數據進行野值剔除、自動增益及幾何位置加載處理。

優選地，在步驟3之前，對所述頻率域數據進行噪音壓制和對不同頻率的吸收衰減按需要進行或不進行補償處理。

優選地，對每個檢波點處采集的數據重復步驟2至步驟5獲得整個剖面的視波阻抗比率或視波阻抗成像數據。

優選地，步驟1中，采集地下介質響應信號的最低頻率應該小于F_min＝0.5·V/D_max，最大頻率應該大于F_max＝0.5·V/D_min，其中D_min為測量的最小深度，D_max為測量的最大深度，V為待測地下介質平均速度。

優選地，步驟1中，采集地下介質響應信號的時間長度應大于10/F_min，數據采集過程中的數據采樣時間間隔為ΔT，且

優選地，對數據進行自動增益數據處理時，其窗口時長大于1/F_min。

優選地，步驟2中所需傅里葉變換的時間域數據的時間長度大于1/F_min。

優選地，步驟4中的標準井參數為實際的近測點處的地質參數或根據勘探精度要求人為制作的標準井地質參數，根據公式

U_I·Г＝U_M