[發明專利]一種地震數據處理方法及裝置有效
| 申請號： | 201711259624.5 | 申請日： | 2017-12-04 |
| 公開（公告）號： | CN108107470B | 公開（公告）日： | 2019-06-11 |
| 發明（設計）人： | 詹仕凡;李磊;陶春峰;萬忠宏;崔京彬;冉賢華 | 申請（專利權）人： | 中國石油天然氣集團公司;中國石油集團東方地球物理勘探有限責任公司 |
| 主分類號： | G01V1/30 | 分類號： | G01V1/30 |
| 代理公司： | 北京三友知識產權代理有限公司 11127 | 代理人： | 李輝;陳剛 |
| 地址： | 100007 *** | 國省代碼： | 北京;11 |
| 權利要求書： | 查看更多 | 說明書： | 查看更多 |
| 摘要： | |||
| 搜索關鍵詞： | 地震道數據 子波 殘差 地震數據處理 預設 疊加 疊前地震數據 希爾伯特變換 數據減去 原始地震 閾值時 分辨率 申請 | ||
1.一種地震數據處理方法，其特征在于，提供有目的層位的原始疊前地震道數據；所述方法包括：
對所述原始疊前地震道數據進行希爾伯特變換，得到希爾伯特數據；
確定所述原始疊前地震道數據對應的目標子波；其中，所述目標子波根據預設子波模型、中心時間和目標主頻確定；所述目標主頻根據所述預設子波模型、所述原始疊前地震道數據和所述中心時間確定；所述中心時間按照下述方法確定：確定所述希爾伯特數據中多個極大值數據點以及所述極大值數據點對應的采樣時間，其中，所述極大值數據點與采樣時間一一對應，并計算所述多個極大值數據點分別與所述目的層位之間的距離，所述極大值數據點與所述目的層位之間的距離為所述極大值數據點對應的采樣時間與所述目的層位對應的預設采樣時間之間的差值的絕對值，并將所述多個極大值數據點中與所述目的層位之間的距離最小的極大值數據點作為目標極大值數據點，并將所述目標極大值數據點對應的采樣時間作為所述中心時間；
將所述原始地震道數據減去所述目標子波，得到殘差地震道數據；
當所述殘差地震道數據中的振幅的絕對值的最大值大于或等于預設振幅閾值時，將所述殘差地震道數據作為新的原始疊前地震道數據，并生成所述新的原始疊前地震道數據對應的新的目標子波和新的殘差地震道數據，直至所述新的殘差地震道數據中振幅的絕對值的最大值小于所述預設振幅閾值，并將所述新的目標子波以及在所述新的目標子波之前生成的目標子波進行疊加，并將疊加的結果作為目標疊前地震道數據。
2.根據權利要求1所述的一種地震數據處理方法，其特征在于，所述目標主頻按照下述方法確定：
根據所述預設子波模型和所述中心時間，確定所述原始疊前地震道數據對應的初始子波模型；其中，所述初始子波模型中子波的相位為零；
確定所述目標極大值數據點對應的瞬時頻率，并根據所述瞬時頻率確定目標主頻范圍；
遍歷所述目標主頻范圍，以獲取初始主頻，并將所述初始主頻作為所述初始子波模型的主頻，得到所述原始疊前地震道數據對應的初始子波，并計算所述初始子波與所述原始疊前地震道數據的相似度；其中，所述目標主頻范圍中的初始主頻與相似度一一對應；所述目標主頻范圍中的多個初始主頻對應多個相似度；
將所述多個相似度中數值最大的相似度對應的主頻作為所述目標主頻。
3.根據權利要求2所述的一種地震數據處理方法，其特征在于，所述計算所述初始子波與所述原始疊前地震道數據的相似度，包括：
計算所述初始子波與所述原始疊前地震道數據之間的內積，并將內積結果作為所述相似度。
4.根據權利要求2所述的一種地震數據處理方法，其特征在于，采用下述公式確定目標主頻范圍：
(fk-fT)≤f≤(fk+fT)
其中，f表示所述目標主頻范圍內的任一主頻，fk表示所述瞬時頻率。
5.根據權利要求1所述的一種地震數據處理方法，其特征在于，采用下述公式表征所述預設子波模型：
其中，gγ(x)表示所述預設子波模型，u表示子波的中心時間，ζ表示子波的主頻，表示子波的相位，x表示采樣時間。
6.根據權利要求1所述的一種地震數據處理方法，其特征在于，所述預設振幅閾值為所述原始疊前地震道數據中的振幅的絕對值的最大值乘以5％。
該專利技術資料僅供研究查看技術是否侵權等信息，商用須獲得專利權人授權。該專利全部權利屬于中國石油天然氣集團公司;中國石油集團東方地球物理勘探有限責任公司，未經中國石油天然氣集團公司;中國石油集團東方地球物理勘探有限責任公司許可，擅自商用是侵權行為。如果您想購買此專利、獲得商業授權和技術合作，請聯系【客服】
本文鏈接：http://www.szxzyx.cn/pat/books/201711259624.5/1.html，轉載請聲明來源鉆瓜專利網。
