[發明專利]快速獲取地震勘探觀測系統面元信息的方法及裝置

有效

• 申請號: 201410814786.0
• 申請日: 2015-08-03
• 公開（公告）號: CN104502967A
• 公開（公告）日: 2015-07-29
• 發明（設計）人: 黃興貴;寧克巖
• 申請（專利權）人: 中國石油天然氣集團公司;中國石油集團東方地球物理勘探有限責任公司
• 主分類號: G01V1/28
• 分類號: G01V1/28
• 代理公司: 北京三友知識產權代理有限公司 11127
• 代理人: 黨曉林;李永強

技術領域

本發明涉及地球物理勘探技術領域，尤其是涉及一種快速獲取地震勘探觀測系統面元信息的方法及裝置。

背景技術

地震勘探設計是地震勘探過程中的首要任務之一，在設計過程中需要對地震勘探觀測系統進行充分地論證，用面元的覆蓋次數、炮檢對中點坐標、炮檢對距離以及炮檢連線方位角等數據信息進行全面的分析論證。面元信息計算就是找出地震勘探工區中每個面元的這些屬性值，供技術人員分析使用，因此，面元信息是石油勘探中觀測系統優化設計的重要數據信息。

目前石油勘探已經進入了精細勘探時代，三維和高分辨勘探已成為主要的勘探施工方案，觀測系統的優化設計中就需要計算大量的各種類型的面元信息數據。這些觀測系統一般都有上億個炮檢對，計算出的面元信息達到幾十GB乃至幾百GB的數據量。而采用傳統的面元信息計算方法是：如計算第(i,j)面元的信息，則從炮點開始循環搜索，依據炮檢關系，循環搜索檢點，測試該炮檢對的中點是否在該面元內，若是則計算其他信息，否則繼續下個點。這種計算方法耗費大量的計算時間用于查找符合條件的炮檢對，從而給后續面元分析工作帶來了瓶頸問題，對工作效益造成巨大影響。

發明內容

本發明的目的在于提供一種快速獲取地震勘探觀測系統面元信息的方法及裝置，以提高面元信息的獲取效率和靈活性。

為達到上述目的，一方面，本發明提供了一種快速獲取地震勘探觀測系統面元信息的方法，包括以下步驟：

構建工區內待計算面元的邏輯炮檢點矩形網格；

基于所述邏輯炮檢點矩形網格建立所述工區內待計算面元的面元矩陣；

對所述面元矩陣進并行計算獲取所述旋轉后工區的面元信息。

本發明的獲取地震勘探觀測系統面元信息的方法，所述構建工區內待計算面元的邏輯炮檢矩形網格，具體包括：

將待計算檢點面元的工區按預設角度旋轉，使旋轉后工區內的檢點的主測線方向平行于坐標軸，獲得旋轉后工區內的邏輯檢點；

平移所述邏輯檢點的分布，使所述邏輯檢點的檢線距及檢點距均勻分布；

依據平移后的邏輯檢點的分布，通過添加虛擬檢點的方式使每條檢線上的邏輯檢點相同，獲取邏輯檢點數據結構；

依據上述方法獲取邏輯炮點數據結構；

將邏輯檢點坐標及邏輯炮點坐標對應讀入所述邏輯檢點數據結構及所述邏輯炮點數據結構，并依據所述邏輯檢點的主測線方向重新調整所述邏輯炮點分布，使得調整后的邏輯炮點的炮線方向與所述邏輯檢點的檢線相同，從而獲取工區內待計算面元的邏輯炮檢矩形網格。

本發明的獲取地震勘探觀測系統面元信息的方法，所述面元矩陣為[nx,ny]，其中，nx表示X方向的面元數，ny表示Y方向的面元數，X方向的面元大小表示為dfx，Y方向的面元大小表示為dfy。

本發明的獲取地震勘探觀測系統面元信息的方法，所述對所述面元矩陣進并行計算獲取所述旋轉后工區的面元信息，具體包括：

根據公式isl0＝(i_fex-maxHalfOffset)/dsl計算所述旋轉后工區的第i行第j列面元相關的首炮線isl0；

根據公式isl1＝(i_fex+maxHalfOffset)/dsl計算所述旋轉后工區的第i行第j列面元相關的終炮線isl1；

根據公式isp0＝(j_fey-maxHalfOffset)/dsp計算所述旋轉后工區的第i行第j列面元相關的首炮點isp0；

根據公式isp1＝(j_fey+maxHalfOffset)/dsp計算所述旋轉后工區的第i行第j列面元相關的終炮點isp1；

根據公式irl0＝(i_fex-maxHalfOffset)/drl計算所述旋轉后工區的第i行第j列面元相關的首檢線irl0；

根據公式irl1＝(i_fex+maxHalfOffset)/drl計算所述旋轉后工區的面元的終檢線irl1；

根據公式irp0＝(j_fey-maxHalfOffset)/drp計算所述旋轉后工區的第i行第j列面元相關的首檢點irp0；

根據公式irp1＝(j_fey+maxHalfOffset)/drp計算所述旋轉后工區的第i行第j列面元相關的終檢點irp1；

其中，dsl、dsp分別表示炮線間隔、炮點間隔，drl、drp分別表示檢線間隔、檢點間隔，maxHalfOffset表示最大炮檢距的一半，i_fex、j_fey分別表示第i行第j列面元的橫坐標和縱坐標。