[發明專利]基于廊帶數據約束的大地電磁剖面三維結構化反演方法

有效

• 申請號: 202110929060.1
• 申請日: 2021-08-13
• 公開（公告）號: CN113391362B
• 公開（公告）日: 2021-10-29
• 發明（設計）人: 王緒本;楊鈺菡;羅威;李德偉;高永才;陳進超
• 申請（專利權）人: 成都理工大學
• 主分類號: G01V3/38
• 分類號: G01V3/38
• 代理公司: 成都東恒知盛知識產權代理事務所(特殊普通合伙) 51304
• 代理人: 羅江

本發明公開基于廊帶數據約束的大地電磁剖面三維結構化反演方法，包括以下步驟：S1：獲取大地電磁的二三維電性剖面，確定目標異常體，分析目標異常體深度與橫向分布情況；S2：根據目標異常體的分布情況，計算目標異常體的影響范圍；S3：基于目標異常體的影響范圍，確定廊帶區域的范圍及加設廊帶輔助測點的位置與數量；形成廊帶三維觀測模式，進行三維結構化反演，得到大地電磁剖面三維反演結果；本發明通過分析二三維實測電性剖面的異常體分布情況，確定剖面其異常體影響范圍，通過加設廊帶輔助測點構成局部三維廊帶觀測方式，進行結構化反演，可以有限約束旁側異常體的影響，提高目標剖面的可靠性，從而實現剖面三維精細化處理與解釋。

技術領域

本發明屬于數據處理技術領域，尤其涉及一種基于廊帶數據約束的大地電磁剖面三維結構化反演方法。

背景技術

由于體積效應存在，大地電磁的二、三維電性剖面易受旁側異常體影響。不存在于實測剖面上的旁側異常體，其響應范圍可能會擴散到二、三維實測電性剖面。而基于二維的地下電性結構假設計算的二維反演處理方法和只有單剖面數據參與反演的三維反演處理方法，會使該剖面的二、三維反演結果上產生假異常，與真實模型存在較大差異。影響解釋人員對該剖面的解釋分析，錯誤刻畫目標體形態。

發明內容

本發明提供基于廊帶數據約束的大地電磁剖面三維結構化反演方法，旨在解決上述存在的問題。

本發明是這樣實現的，基于廊帶數據約束的大地電磁剖面三維結構化反演方法，包括以下步驟：

S1：獲取大地電磁的二、三維電性剖面，確定目標異常體，分析二、三維電性剖面中沿剖面測線方向的異常體深度與橫向分布情況；

S2：根據目標異常體的分布情況，計算目標異常體的影響范圍；

S3：基于目標異常體的影響范圍，確定廊帶區域的范圍及加設廊帶輔助測點的位置與數量；形成廊帶三維觀測模式，進行三維結構化反演，得到大地電磁剖面三維反演結果；

作為優選的，在步驟S1中，具體包括：

S11：選取測線方向橫向距離最大的異常體作為刻畫廊帶范圍的目標體，目標異常體的長即x剖面方向、寬即y垂直于剖面方向、深度為z方向；

S12：基于有限差分三維正演算法對目標異常體進行三維正演計算，得到第一正演結果；

S13：使用相同正演空間、網格、頻率參數進行均勻背景電阻率下的三維正演計算，得到第二正演結果；

S14：第一正演結果減去第二正演結果，得到處理后數據，即目標異常的分布情況。

進一步的，在步驟S12的三維正演計算中，頻率范圍與所述二、三維電性剖面數據的采集頻段一致；

網格為規則矩形網格剖分，網格間隔需滿足正演異常體其占間隔3個及以上，橫向擴展網格以1.5的幾何增長因子擴展15個，垂向網格按首層厚度以幾何增長因子1.5設置35個；

首層厚度計算公式為：

其中，h為首層厚度，ρ_top為首層電阻率值，f_high為最高頻率。

進一步的，在步驟S12中，進行三維正演計算使用的正演方程組包括：