[發明專利]井震標定方法及系統

在審

• 申請號: 201911263986.0
• 申請日: 2019-12-11
• 公開（公告）號: CN112946750A
• 公開（公告）日: 2021-06-11
• 發明（設計）人: 周玉冰;田作基;劉洪昌;張志偉;馬中振;陽孝法;趙永斌;王丹丹;劉亞明
• 申請（專利權）人: 中國石油天然氣股份有限公司
• 主分類號: G01V1/40
• 分類號: G01V1/40;G01V1/28;G01V1/36
• 代理公司: 北京三友知識產權代理有限公司 11127
• 代理人: 趙平;周永君

1.一種井震標定方法，其特征在于，包括：

根據多井井旁地震道數據時頻譜的能量分布特征對所述多井井旁地震道數據進行頻段劃分得到頻段地震道數據；

將預設的地震子波庫中不同地震子波與預設的測井反射系數序列褶積得到合成地震道數據；

根據給定時移范圍平移所述合成地震道數據，計算所述合成地震道數據與所述多井井旁地震道數據和所述頻段地震道數據的相關系數，得到各頻段的時移與子波主頻的相關系數分布圖；

根據所述相關系數分布圖，從所述多井井旁地震道數據中選取與各頻段最大相關系數最接近的子波主頻和對應的時深數據，生成標定結果。

2.根據權利要求1所述的井震標定方法，其特征在于，還包括：

對測井曲線進行校正，得到測井反射系數序列。

3.根據權利要求1所述的井震標定方法，其特征在于，還包括：

根據多井井旁地震道數據的時頻譜構建得到地震子波庫。

4.根據權利要求2所述的井震標定方法，其特征在于，所述對測井曲線進行校正，包括：

對所述測井曲線進行變時窗中值濾波處理；

對中值濾波處理后的所述測井曲線進行約束方波化分層處理，得到方波化的測井反射系數序列。

5.根據權利要求4所述的井震標定方法，其特征在于，所述對所述測井數據進行變時窗中值濾波處理，包括：

選定設定大小的第一時窗，針對所述測井數據中的每個數據點，以每個數據點為起點，對所述第一時窗內的數據點按照大小進行排序生成第一時窗序列，選取每個所述第一時窗序列中的中位數組成第一濾波數據；

執行迭代操作，選定大小大于第一時窗的第二時窗，針對所述第一濾波數據中的每個數據點，以每個數據點為起點，對所述第二時窗內的數據點按照大小進行排序生成第二時窗序列，選取每個所述第二時窗序列中的中位數組成第二濾波數據，用所述第二濾波數據替代所述第一濾波數據，重新選定大小大于第一時窗的第二時窗，直至所述第二濾波數據滿足預設條件。

6.根據權利要求4所述的井震標定方法，其特征在于，所述對中值濾波處理后的所述測井曲線進行約束方波化分層處理，得到方波化的測井反射系數序列，包括：

獲取所述測井曲線中的深度域的n個采樣點以及對應各個點的速度值，將起始采樣點所處平面作為第一個地層界面；其中任意兩個采樣點之間的間隔為設定距離；

對起始采樣點執行迭代操作，計算起始采樣點與相鄰的下一個采樣點之間的速度差的絕對值，若速度差的絕對值小于最小速度差，或者若速度差的絕對值大于等于最小速度差并且后一個采樣點的深度小于預設的巖層最小厚度，則用相鄰兩個采樣點中的前一個采樣點的速度值替代后一個采樣點的速度值，否則對與所述起始采樣點相鄰的下一個采樣點執行所述迭代操作，直至遍歷所有采樣點。

7.根據權利要求1所述的井震標定方法，其特征在于，還包括：

利用小波多尺度變換繪制目的層多井井旁地震道數據時頻譜。

8.根據權利要求7所述的井震標定方法，其特征在于，所述利用小波多尺度變換繪制目的層多井井旁地震道數據時頻譜，包括：

通過平移小波基函數獲取信號的時間信息；

通過縮放小波基函數的寬度獲得信號的頻率特性。

9.根據權利要求1所述的井震標定方法，其特征在于，還包括：

根據多井井旁地震道的時頻譜特征，確定目的層井旁地震道數據的主頻范圍；

根據目的層井旁地震道數據的主頻范圍設置地震子波主頻范圍；

根據所述地震子波主頻范圍，利用Ricker子波構建地震子波庫。