[發明專利]一種基于復電阻率參數的礦體位置確定方法及系統

有效

• 申請號: 202011301082.5
• 申請日: 2020-11-19
• 公開（公告）號: CN112394420B
• 公開（公告）日: 2023-10-20
• 發明（設計）人: 張志勇;劉祜;汪碩;黃笑;王寒冰;王春陽
• 申請（專利權）人: 核工業北京地質研究院
• 主分類號: G01V3/38
• 分類號: G01V3/38
• 代理公司: 北京高沃律師事務所 11569
• 代理人: 杜陽陽

1.一種基于復電阻率參數的礦體位置確定方法，其特征在于，包括：

獲取Cole-Cole模型；

將所述Cole-Cole模型重參數化，得到重參數化的Cole-Cole模型；

根據所述重參數化的Cole-Cole模型，利用所述重參數化的Cole-Cole模型中的復電阻率值確定所述重參數化的Cole-Cole模型的靈敏度；

利用所述靈敏度，采用阻尼最小二乘法反演理論進行反演，得到所述重參數化的Cole-Cole模型中的參數反演值；

根據所述參數反演值確定礦體位置。

2.根據權利要求1所述的基于復電阻率參數的礦體位置確定方法，其特征在于，所述將所述Cole-Cole模型重參數化，得到重參數化的Cole-Cole模型，具體包括：

將所述Cole-Cole模型中的時間常數替換為設定時間常數，將所述Cole-Cole模型中的極化率替換為相位角，得到重參數化的Cole-Cole模型。

3.根據權利要求2所述的基于復電阻率參數的礦體位置確定方法，其特征在于，所述將所述Cole-Cole模型中的時間常數替換為設定時間常數，將所述Cole-Cole模型中的極化率替換為相位角，得到重參數化的Cole-Cole模型，具體包括：

根據如下公式替換所述Cole-Cole模型中的時間常數和所述Cole-Cole模型中的極化率：

其中，m₀表示極化率，表示相位角，a(n)表示第一參數，b(n)表示第二參數，n表示模型參數的序號，τ_ρ表示時間常數，ρ表示電阻率，表示設定時間常數，表示相位，Re()表示取復數的實部，Im()表示取復數的虛部，i表示復數的虛部，c表示頻率相關系數。

4.根據權利要求1所述的基于復電阻率參數的礦體位置確定方法，其特征在于，所述根據所述重參數化的Cole-Cole模型，利用所述重參數化的Cole-Cole模型中的復電阻率值確定所述重參數化的Cole-Cole模型的靈敏度，具體包括：

改變所述重參數化的Cole-Cole模型中的任一參數，計算參數改變后的復電阻率；

根據所述參數改變后對應的復電阻率確定所述重參數化的Cole-Cole模型的靈敏度。

5.根據權利要求1所述的基于復電阻率參數的礦體位置確定方法，其特征在于，所述利用所述靈敏度，采用阻尼最小二乘法反演理論進行反演，得到所述重參數化的Cole-Cole模型中的參數反演值，具體包括：

根據所述靈敏度，確定所述重參數化的Cole-Cole模型參數反演的初值；

以單個Cole-Cole模型表示復電阻率正演模型，利用相對偏差確定所述重參數化的Cole-Cole模型參數反演的初值處的泰勒展開式；所述相對偏差為所述復電阻率正演模型理論正演場值和實測場值之間的擬合程度；

根據所述泰勒展開式確定擬合誤差；

判斷所述擬合誤差是否小于設定擬合誤差，得到第一判斷結果；

若所述第一判斷結果表示為是，則確定所述重參數化的Cole-Cole模型參數反演的初值為所述重參數化的Cole-Cole模型中的參數反演值；

若所述第一判斷結果表示為否，則根據所述擬合誤差確定異常體電阻率模型參數；

根據所述異常體電阻率模型參數和所述重參數化的Cole-Cole模型參數反演的初值之和更新所述重參數化的Cole-Cole模型參數反演的初值，并返回步驟“以單個Cole-Cole模型表示復電阻率正演模型，利用相對偏差確定所述重參數化的Cole-Cole模型參數反演的初值處的泰勒展開式”。