[發明專利]一種面向分層大地的視在電阻率測算方法、系統及存儲介質

在審

• 申請號: 202211326767.4
• 申請日: 2022-10-25
• 公開（公告）號: CN115718325A
• 公開（公告）日: 2023-02-28
• 發明（設計）人: 吉旺威;林藝;李國城;李光茂;鄭服利;鐘少泉;何昊
• 申請（專利權）人: 廣東電網有限責任公司廣州供電局
• 主分類號: G01V3/08
• 分類號: G01V3/08;G01V3/38
• 代理公司: 廣州粵高專利商標代理有限公司 44102
• 代理人: 劉俊

1.一種面向分層大地的視在電阻率測算方法，其特征在于，包括以下步驟：

S1、根據溫納四極法，初始化設置極距a，對電流極輸入測試電流I，獲取對應極距a的測試電壓ΔV，計算視在電阻率測量值ρ_am；

S2、引入歐拉變換法獲得兩個電壓極間的地表電位差ΔV′(a)，根據地表電位差ΔV′(a)計算分層大地對應極距a的視在電阻率計算值ρ_as；其中，采用歐拉變換法對地表電位差ΔV′(a)進行計算時，對極距a對應的參數以預設的間隔及預設的數量進行采樣；

S3、計算測量值ρ_am和計算值ρ_as間的相對誤差ε，根據相對誤差ε是否超過預設閾值確認是否更新極距a：若超過預設閾值，則確認更新極距a并轉入步驟S1；若不超過預設閾值，則不更新極距a，結束計算并輸出視在電阻率計算值ρ_as。

2.根據權利要求1所述的一種面向分層大地的視在電阻率測算方法，其特征在于，所述步驟S1中，測試電流I為0.1mA-10mA。

3.根據權利要求1所述的一種面向分層大地的視在電阻率測算方法，其特征在于，所述步驟S1中，所述極距a為1m以上。

4.根據權利要求1所述的一種面向分層大地的視在電阻率測算方法，其特征在于，所述步驟S2中，所述歐拉變換法包括以下步驟：

S2.1、對應極距a構建關于地表電位差ΔV′(a)的模型：

其中，λ為積分變量，表示電磁場傳播系數；ρ₁表示水平分層大地第1層大地的電阻率；J₀為第一類零階貝塞爾函數；R₁為表達地表特征阻抗的函數，其形式如下：

其中，R_i表示水平分層大地第i層的特征阻抗；ρ_i表示第i層大地的電阻率；h_i表示第i層大地的分層厚度；m表示水平分層大地最深層所在層數；

S2.2、對式(1)進行漢克爾變換，所述漢克爾變換結果表達式為：

其中，J₀為第一類零階貝塞爾函數；

S2.3、引入歐拉變換法，所述式(3.1)和式(3.2)變換結果為：

其中,λ＝e^-y；a＝e^x；x,y∈(-∞,∞)；

S2.4、按間隔Δ對式(4)中e^x-y進行采樣，所述式(4)輸出為離散信號如下：

其中，H_n為濾波系數；Δ_n為采樣點位置；此時積分核函數f＝R₁；

S2.5、采取90個采樣點形式，對式(5)進行截斷求和，更新地表電位差ΔV′(a)模型，其表達式為：

其中，Δ_k為預設的位置系數，H_k為預設的權重系數。

5.根據權利要求4所述的一種面向分層大地的視在電阻率測算方法，其特征在于，所述步驟S2.4中，所述采樣間隔Δ小于0.5倍積分核函數f的截止頻率。

6.根據權利要求4所述的一種面向分層大地的視在電阻率測算方法，其特征在于，所述步驟S2.5中預設的Δ_k和預設的H_k如下表所示：

7.根據權利要求1所述的一種面向分層大地的視在電阻率測算方法，其特征在于，所述步驟S3中，所述相對誤差ε表達式為：