[发明专利]一种基于磁场测量的堤坝渗漏通道检测方法及其装置

权利要求书

1.一种基于磁场测量的堤坝渗漏通道检测方法，其特征在于：包括如下步骤：

S1：获取信号发射机不供电时堤坝上每个磁观测点上的背景磁感应强度矢量以及供电时每个磁观测点上的总磁感应强度矢量；

其中，在堤坝的两侧布设发射电极A和B，所述发射电极A位于水库中并与水接触，所述发射电极A和B之间连接信号发射机，在所述提坝上按照预设线距和点距排布磁观测点；

S2：分别计算每个磁观测点上的总磁感应强度矢量与背景磁感应强度矢量的差值得到每个磁观测点的磁异常矢量；

其中，基于每个磁观测点的磁异常矢量构成提坝的磁异常场；

S3：基于磁异常场中每个磁观测点的磁异常矢量推演出电流密度分布；

其中，电流密度分布不为零处对应为提坝上渗漏通道的位置。

2.根据权利要求1所述的方法，其特征在于：步骤S3是基于磁异常场与电流密度场的关系推演出电流密度分布；

其中，将提坝三维空间离散为M个单元，每个单元具有相同的体积以及每个单元中存在恒定电流密度，磁异常场与电流密度场的关系如下：

B_s＝KJ

式中，B_s为磁异常场矩阵，矩阵大小为3N×1，分别为第1个、第N个磁观测点上的磁异常矢量，N为磁观测点的个数；K为位置矩阵，矩阵大小为3N×3M，分别为第1个磁观测点与第1个单元、第M个单元之间的距离矩阵，分别为第N个磁观测点与第1个单元、第M个单元之间的距离矩阵；J为电流密度场矩阵，矩阵大小为3M×1，J¹、J^M分别为第1个单元、第M个单元上的电流密度；

其中，任意第i个磁观测点的磁异常矢量第i个磁观测点与第k个单元之间的距离矩阵第k个单元上的电流密度J^k如下所示：

式中，分别为磁异常矢量在空间直角坐标系xyz中的x、y、z方向上的分量，μ为自由空间的磁导率，ν为单元体积，r_ik为第i个磁观测点与第k个单元之间的空间距离，分别为第i个磁观测点与第k个单元之间的空间距离在空间直角坐标系xyz中的x、y、z方向上的分量；分别为第k个单元上的电流密度J^k在空间直角坐标系xyz中的x、y、z方向上的分量。

3.根据权利要求2所述的方法，其特征在于：步骤S3中基于磁异常场与电流密度场的关系采用磁异常曲线的特征点法、切线法、二维反演方法或线性三维反演方法推演出电流密度分布。

4.根据权利要求3所述的方法，其特征在于：采用线性三维反演方法推演出电流密度分布的过程如下：

首先，根据最小乘算法原理建立目标函数Φ(J)并得到目标函数的线性逆问题的解决方案；

其中，目标函数如下：

Φ(J)＝[W_d(d_obs-KJ)]^T[W_d(d_obs-KJ)]+λ[W_m(J-J₀)]^T[W_m(J-J₀)]

解决方案的公式如下：

式中，λ为正则化因子，J₀为电流密度的先验模型，d_obs为观测数据，W_d、W_m均为权重矩阵，权重矩阵W_d为观测数据d_obs的误差协方差矩阵的倒数组成的对角矩阵，权重矩阵W_m为光滑度矩阵；

然后，将提坝上的磁异常场矩阵以及位置矩阵代入解决方案的公式中计算出电流密度场矩阵得到提坝上的电流密度分布；

其中，磁异常场矩阵B_s为观测数据d_obs。

5.根据权利要求1所述的方法，其特征在于：在观测点利用磁探头采集磁感应强度矢量。