[发明专利]阵列感应视电导率权系数计算方法及装置

权利要求书

1.一种阵列感应视电导率权系数计算方法，其特征在于，包括：

获取各子线圈系视电导率的采样点深度值；

将所述采样点深度值代入最小二乘法的寻优函数，并根据所述寻优函数构建各个所述子线圈系视电导率的权系数线性方程组；

从所述权系数线性方程组提取权系数的实部方程组和虚部方程组的系数矩阵；

对所述实部方程组的系数矩阵和常数项向量进行分块；

根据分块的所述系数矩阵和所述常数项向量计算对应子线圈系的权系数的实部；

根据所述子线圈系的权系数的实部计算当前的总几何因子；

计算所述总几何因子与实部目标几何因子的差值；

若所述差值小于预设阈值，则对所述权系数的实部进行归一化处理得到归一化的权系数的实部；

根据所述权系数的实部计算实部的总几何因子；

以所述实部的总几何因子减去实部目标几何因子得到差值，并将所述差值与函数的乘积赋给权系数虚部目标几何因子；其中，a为常数；r为预设柱坐标系的极径坐标；

根据所述虚部目标几何因子得到虚部方程组的常数项向量；

根据所述常数项向量和所述系数矩阵，利用最小二乘法构建权系数的虚部方程组；

根据所述虚部方程组计算权系数的虚部；

根据所述实部和所述虚部生成各个所述子线圈系视电导率的权系数。

2.根据权利要求1所述的阵列感应视电导率权系数计算方法，其特征在于，所述对所述实部方程组的系数矩阵和常数项向量进行分块的步骤，包括：

按照子线圈系视电导率的采样点个数，对所述实部方程组的系数矩阵和常数项向量进行分块。

3.根据权利要求1所述的阵列感应视电导率权系数计算方法，其特征在于，所述根据分块的所述系数矩阵和所述常数项向量计算对应子线圈系的权系数的实部的步骤，包括：

设定第一临时变量p＝0，并给定M，M-1，...M-x_p值，并将M-x₀赋值给第二临时变量m；所述第一临时变量用于记录计算子线圈系权系数的次数；其中，M表示当前需要的最长源矩的子阵列编号；

选取分块的系数矩阵中M到M-x₀行与M到M-x₀列的子矩阵，分块的常数项向量中M到M-x₀行的子向量，以根据所述子矩阵和所述子向量构建关于M、M-1、...、M-x₀号子线圈系权系数的实部线性方程组；

求解所述实部线性方程组得到M、M-1、...、M-x₀号子线圈系权系数的实部。

4.根据权利要求3所述的阵列感应视电导率权系数计算方法，其特征在于，在所述求解所述实部的线性方程组得到M、M-1、...、M-x₀号子线圈系权系数的实部的步骤之前，还包括：

利用松弛不完全LU分解法和近似逆法对所述实部线性方程组进行预处理，以降低所述实部线性方程组的方程条件数。

5.根据权利要求3所述的阵列感应视电导率权系数计算方法，其特征在于，所述方法还包括：

若所述差值大于预设阈值，则所述第一临时变量p自加1，所述第二临时变量m自减1；

根据当前总的几何因子和实部目标几何因子之差建立第m、m-1、...、m-x_p号子线圈系权系数的实部线性方程组；

根据所述实部线性方程组计算第m、m-1、...、m-x_p号子线圈系权系数的实部。

6.一种阵列感应视电导率权系数计算装置，其特征在于，包括：

采样点深度值获取模块，用于获取各子线圈系视电导率的采样点深度值；

权系数线性方程组构建模块，用于将所述采样点深度值代入最小二乘法的寻优函数，并根据所述寻优函数构建各个所述子线圈系视电导率的权系数线性方程组

提取模块，用于从所述权系数线性方程组提取权系数的实部方程组和虚部方程组的系数矩阵；

实部计算模块，用于对所述实部方程组的系数矩阵和常数项向量进行分块；根据分块的所述系数矩阵和所述常数项向量计算对应子线圈系的权系数的实部；根据所述子线圈系的权系数的实部计算当前的总几何因子；计算所述总几何因子与实部目标几何因子的差值；若所述差值小于预设阈值，则对所述权系数的实部进行归一化处理得到归一化的权系数的实部；

虚部方程组构建模块，用于根据所述权系数的实部计算实部的总几何因子；以所述实部的总几何因子减去实部目标几何因子得到差值，并将所述差值与函数的乘积赋给权系数虚部目标几何因子；其中，a为常数；r为预设柱坐标系的极径坐标；

根据所述虚部目标几何因子得到虚部方程组的常数项向量；根据所述常数项向量和所述系数矩阵，利用最小二乘法构建权系数的虚部方程组；

虚部计算模块，用于根据所述虚部方程组计算权系数的虚部；

权系数生成模块，用于根据所述实部和所述虚部生成各个所述子线圈系视电导率的权系数。