[發明專利]基于歸一化函數的正演電磁數據擬地震剖面顯示方法

在審

• 申請號: 202011475977.0
• 申請日: 2020-12-14
• 公開（公告）號: CN112817052A
• 公開（公告）日: 2021-05-18
• 發明（設計）人: 田紅軍;胡曉穎;古志文;周印明;尹雪波;譚章坤;杜蛟;張劍;葉恒;王雨辰
• 申請（專利權）人: 四川中成煤田物探工程院有限公司
• 主分類號: G01V3/08
• 分類號: G01V3/08;G01V3/38
• 代理公司: 成都華復知識產權代理有限公司 51298
• 代理人: 龐啟成

1.基于歸一化函數的正演電磁數據擬地震剖面顯示方法，其特征在于，包括：

S1：建立反映地層厚度與電阻率對應關系的層狀地電模型；

S2：設置電磁數據采集參數；

S3：根據水平層狀介質水平電偶極子電磁場公式計算在笛卡爾坐標系下每個發射電流頻率所對應的非均勻半空間場值數據；

S4：設置均勻半空間電阻率參數，根據均勻半空間水平電偶極子電磁場公式計算在笛卡爾坐標系下每個發射電流頻率所對應的均勻半空間場值數據；

S5：將所述非均勻半空間場值數據實部值與所述均勻半空間場值數據實部值相比，得到歸一化場值數據；

S6：以偏移距為橫軸、對數發射電流頻率為縱軸，繪制所述歸一化場值數據的同相軸剖面顯示圖。

2.根據權利要求1所述的顯示方法，其特征在于：步驟S2包括設置發射電流、電偶極子長度；步驟S3包括改變發射電流頻率，根據水平層狀介質水平電偶極子電磁場公式分別計算柱坐標系下的場值，再根據柱坐標系下的場值轉換為笛卡爾坐標系下的場值。

3.根據權利要求2所述的顯示方法，其特征在于：對于長導線源：步驟S3將所述長導線源分解成多個偶極子源分別計算，得到多個中間非均勻場值，再將多個所述中間非均勻場值累加得到最終非均勻半空間場值；步驟S4將所述長導線源分解成多個偶極子源分別計算，得到多個中間均勻場值，再將多個所述中間均勻場值累加得到最終均勻半空間場值；步驟S5將所述最終非均勻半空間場值實部值與所述最終均勻半空間場值實部值相比，得到歸一化場值數據。

4.根據權利要求3所述的顯示方法，其特征在于：所述場值為電場強度沿偶極方向水平分量EX、電場強度垂直偶極方向水平分量EY、電場強度垂直分量EZ、磁場強度沿偶極方向水平分量HX、磁場強度垂直偶極方向水平分量HY、磁場強度垂直分量HZ中的任一種。

5.根據權利要求4所述的顯示方法，其特征在于：采用801個濾波系數進行漢克爾變換實現含有0階和1階貝塞爾函數積分的電磁場分量計算。

6.根據權利要求5所述的顯示方法，其特征在于：采用逆插值樣條方法求解非均勻半空間場值，采用平均值法進行多次迭代求取均勻半空間場值。

7.根據權利要求1-6任一項所述的顯示方法，其特征在于：步驟S6對發射頻率10^-1～10⁴Hz的對數域進行N等分，從-1起取N+1個采樣點作為縱軸刻度。

8.根據權利要求7所述的顯示方法，其特征在于：步驟S6繪制同相軸剖面顯示圖之前還包括對歸一化場值放大指定倍數。

9.一種計算機設備，包括存儲器、處理器及存儲在存儲器上并可在處理器上運行的計算機程序，其特征在于，所述處理器執行所述計算機程序時實現權利要求1至8任一所述方法。

10.一種計算機可讀存儲介質，其特征在于，所述計算機可讀存儲介質存儲有執行權利要求1至8任一所述方法的計算機程序。