[发明专利]一种基于复电阻率参数的矿体位置确定方法及系统

权利要求书

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模型参数反演的初值处的泰勒展开式”。