CN113960679A - 一种基于六面体磁探测装置的渗漏检测方法及装置 - Google Patents

一种基于六面体磁探测装置的渗漏检测方法及装置 Download PDF

Info

Publication number
CN113960679A
CN113960679A CN202111254499.5A CN202111254499A CN113960679A CN 113960679 A CN113960679 A CN 113960679A CN 202111254499 A CN202111254499 A CN 202111254499A CN 113960679 A CN113960679 A CN 113960679A
Authority
CN
China
Prior art keywords
magnetic field
magnetic
power supply
field gradient
measuring point
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202111254499.5A
Other languages
English (en)
Other versions
CN113960679B (zh
Inventor
郭玉松
汤井田
王小鹏
胡双贵
郭士明
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Yellow River Engineering Consulting Co Ltd
Original Assignee
Yellow River Engineering Consulting Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Yellow River Engineering Consulting Co Ltd filed Critical Yellow River Engineering Consulting Co Ltd
Priority to CN202111254499.5A priority Critical patent/CN113960679B/zh
Publication of CN113960679A publication Critical patent/CN113960679A/zh
Application granted granted Critical
Publication of CN113960679B publication Critical patent/CN113960679B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V3/00Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
    • G01V3/08Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation operating with magnetic or electric fields produced or modified by objects or geological structures or by detecting devices
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01PMEASURING LINEAR OR ANGULAR SPEED, ACCELERATION, DECELERATION, OR SHOCK; INDICATING PRESENCE, ABSENCE, OR DIRECTION, OF MOVEMENT
    • G01P13/00Indicating or recording presence, absence, or direction, of movement
    • G01P13/02Indicating direction only, e.g. by weather vane
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00Arrangements or instruments for measuring magnetic variables
    • G01R33/02Measuring direction or magnitude of magnetic fields or magnetic flux
    • G01R33/022Measuring gradient
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02ATECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A90/00Technologies having an indirect contribution to adaptation to climate change
    • Y02A90/30Assessment of water resources

Abstract

本发明公开了一种基于六面体磁探测装置的渗漏检测方法及装置,该检测方法包括以下步骤:步骤1:在堤坝的两侧布设供电电极A和B;步骤2:供电装置通过供电电极A、B向堤坝和水中供入电流;步骤3:利用六面体磁探测装置获取通电时每个测点对应的磁场梯度G;步骤4:基于每个测点的磁场梯度G获取渗漏通道的分布,其中,磁场梯度G分布中变化最大的区域为渗漏通道的分布区域。本发明提供了一种全新的手段确定渗漏通道的分布,即根据磁梯度数据来确定渗漏入口的位置和堤坝中渗漏通道的分布,丰富了渗漏检测技术。

Description

一种基于六面体磁探测装置的渗漏检测方法及装置
技术领域
本发明属于勘探地球物理技术领域,具体涉及一种基于六面体磁探测装置的渗漏检测方法及装置。
背景技术
作为水利基础工程建设的组成之一的堤坝,在水利治理体系中占有重要地位。据统计,截止2017年,我国已建成5级以上堤坝30多万公里,已建成各类水库将近10万座。其中,其中大约有30%的坝堤存在不同程度的病险隐患,其表现形式为不同程度下的渗漏。当洪水来临时,由于其高水头渗透作用,平时一些细小隐患很可能引发重大险情。因此对堤坝的结构隐患的探测成为防止堤坝在某一时期突然失效的重要举措。
目前应用于堤坝隐患探测的常用物探方法,包括自然电场法、高密度电阻率法、探地雷达、瞬变电磁法、瑞雷面波法、CT技术以及流场拟合法探测等。中南大学何继善院士提出的流场拟合法,利用水流场与电流场在一定条件下的某些数学物理之间的相似性,在管涌通道出口至堤坝前水域建立人工电流场来拟合渗漏水流场,通过测定电流密度的分布来查明水流场的流向和相对流速取得了良好的效果。但是由于我国堤坝类型众多、结构多样、渗漏情况复杂,所以当前检测技术的探测精度和速度不能满足要求。由于磁场测量具有快速、高效、经济以及磁梯度受地磁场影响较小等特点,故利用磁场测量来确定堤坝渗漏入口和水流方向有着十分重要的意义,可以为堤坝渗漏通道的快速检测提供一种全新的技术手段。
而现有方法中利用磁场测量堤坝渗漏的技术也多是测量坝堤顶部磁场,进而使用复杂算法反演电流密度,且所利用的感应式磁传感器仪器笨重,不方便使用,因此需要提供一种更加简单、快速的方法来确定渗漏通道及水流方向。此外,现有磁场测量堤坝渗漏的技术中常常使用磁异常来反演电流密度分布,从而确定渗漏通道的位置等。为了丰富检测手段,如何充分利用磁技术,开辟其他有效途径检测渗漏通道也是有待研究的。
发明内容
本发明的目的是提供了一种基于六面体磁探测装置的渗漏检测方法及装置,所述方法在水中和堤坝上测量磁场,并根据所获得的磁场计算磁场梯度,根据磁场梯度数据来确定渗漏入口的位置和堤坝中渗漏通道的分布。本发明提供了另一个手段来确定渗漏通道的分布,丰富了渗漏检测技术。
一方面,本发明提供.一种基于六面体磁探测装置的渗漏检测方法,其包括以下步骤:
步骤1:在堤坝的两侧布设供电电极A和B,其中,供电电极A位于堤身内侧渗漏通道的出水口,供电电极B位于堤坝迎水坡一侧,供电电极通过导线与供电装置相连接;
步骤2:供电装置通过供电电极A、B向堤坝和水中供入电流;
步骤3:利用六面体磁探测装置获取通电时每个测点对应的磁场梯度G;
步骤4:基于每个测点的磁场梯度G获取渗漏通道的分布,其中,磁场梯度G分布中变化最大的区域为渗漏通道的分布区域。
关于渗漏通道的位置,使用发电机通过导线向供电电极A和B注入电流信号,电流沿着堤坝渗漏通道产生的优先导电路径引导时,由于通道的存在所产生的电流扰动会产生磁异常,当电流扰动越大,其产生的磁异常越大,同时其梯度变化就越明显。通过处理和分析所获得的磁场梯度数据就能判断通道的分布,而本发明之所以创造性选择用磁场梯度来反应渗漏通道的位置,是考虑到相比于磁场矢量,磁场梯度数据具有九个分量,能更好、更准确的反映渗漏通道的位置,同时能抵抗随机噪声和地磁场的干扰,精度也进一步提高。
可选地,步骤3中所述每个测点对应的磁场梯度G的获取方式为:基于每个测点对应的总磁感应强度并采用中心差分方式得到测点的磁场梯度G。
本发明选择采用中心差分方式得到磁场梯度,并通过验证得知,在本应用领域使用中心差分方式得到的磁场梯度与理论磁场的误差很小,因此,利用中心差分方式得到的磁场梯度是可靠的。
以六面体上的八个磁传感器为例,本发明使用中心差分的方式计算磁场梯度,首先通过六面体装置八个顶点上的磁传感器,测得由渗漏通道、导线、背景及环境所产生的磁场总和。接着用两点差分计算出六面体装置两个相对面之间四条棱边中点的梯度,再通过将四个梯度进行平均计算获得六面体中心点位置处的磁场梯度。该方式不仅可以充分利用八个磁传感器的磁场,同时获得可以计算三个方向电流密度的梯度分量,而且可以很大程度上消除由导线和外部环境引起的外部磁场。
可选地,每个测点上设有八个磁传感器,利用中心差分方式得到测点的磁场梯度G的公式如下:
Figure BDA0003323419750000021
Figure BDA0003323419750000022
Figure BDA0003323419750000031
Figure BDA0003323419750000032
Figure BDA0003323419750000033
Figure BDA0003323419750000034
Figure BDA0003323419750000035
Figure BDA0003323419750000036
Figure BDA0003323419750000037
Figure BDA0003323419750000038
式中,B为磁场矢量,Bx、By、Bz表示在磁观测点上x、y、z方向的磁场分量,
Figure BDA0003323419750000039
表示第i个磁传感器在k方向的磁场分量,d为任意两个磁传感器之间的距离,
Figure BDA00033234197500000310
为偏导数符号,
Figure BDA00033234197500000311
为磁场梯度。
基于八个磁通门磁场传感器计算磁场梯度时每一个分量其计算公式比较对称,计算方便,易于理解。其利用的磁场数据更多,磁场梯度计算流程更加精细,计算出的磁场梯度精度进一步提高,进而电流密度计算更加准确,同时减少了后续数据处理的步骤,提高了探测速度。因此,本发明为了提高磁梯度的精度,鄙弃正四面体结构(正四面体构架结构简单,所测磁场数据较少,计算出的磁场梯度误差较大,最终难以准确反映渗漏点),选择八个磁通门磁场传感器,提高了磁梯度的进度,进而最终提高了渗漏通道检测结果的可靠性。
可选地,所述方法还包括:
获取每个测点上的磁场旋度,再基于所述磁场旋度获得测点上的传导电流密度;
根据所述传导电流密度的汇聚方向得到流向渗漏通道的水流方向,所述汇聚方向与所述水流方向相对应。
本发明通过直接计算磁场旋度来获得传导电流密度,相比于利用磁异常反演电流密度同时也更加高效、快速,同时由于采用了正六面体磁梯度计算使梯度计算容易实施,计算精度更高,在测量中仪器重量可以进一步减小,测量易于实施。而传导电流密度分布与水流方向两者之间的关系在于:当水流由坝坡渗漏入口通过渗漏通道流向下游时,水库中的水流会向渗漏口流动,根据电流场与水流场的相似性,传导电流密度汇聚方向即为水流的流动方向。本发明通过测量磁场的方式得到水中的传导电流密度分布情况,进而利用传导电流密度汇聚方向判断流向渗漏通道的水流流向。
可选地,所述传导电流密度的计算公式如下;
Figure BDA0003323419750000041
式中,
Figure BDA0003323419750000042
表示磁观测点处的磁场旋度,μ0为自由空间的磁导率,Jc为传导电流密度。
可选地,每个测点上的六面体磁探测装置为:八个磁传感器固定在六面体支撑件的八个顶点上,所述六面体支撑件为正六面体。其中,优选六面体支撑架采用环氧板材料,可以在保证支撑架强度时,减少其对磁传感器数据的影响。此外,可以适应长时间野外工作,仪器性能有了大幅度提高。
可选地,所述供电装置所发出的电流激励信号为直流信号或者频率为0.1~380Hz的方波或伪随机信号,供电电极A、B与供电模块之间的供电导线U字排列。
第二方面,本发明提供一种基于所述渗漏通道检测方法的装置,其包括:电流装置、测量模块以及检测模块;
其中,所述电流装置包括供电电极A和B、供电装置,供电电极A位于堤身内侧渗漏通道的出水口,供电电极B位于堤坝迎水坡一侧,供电电极通过导线与供电装置相连接;
测量模块为六面体磁探测装置,其设置于每个测点上,其内设有磁传感器;
检测模块与所述测量模块连接,用于接收所述测量模块中磁传感器采集的数据,并进行渗透检测。
第三方面,本发明提供一种检测装置,其包括:
一个或多个处理器;
存储一个或多个程序的存储器;
所述处理器调用所述程序以执行:
获取通电时每个测点对应的磁场梯度G;
基于每个测点的磁场梯度G获取渗漏通道的分布,其中,磁场梯度G分布中变化最大的区域为渗漏通道的分布区域。
第四方面,本发明提供一种可读存储介质,其存储了计算机程序,所述计算机程序被处理器调用以执行:
获取通电时每个测点对应的磁场梯度G;
基于每个测点的磁场梯度G获取渗漏通道的分布,其中,磁场梯度G分布中变化最大的区域为渗漏通道的分布区域。
有益效果
1.本发明首次将根据磁梯度变化判断渗漏通道的分布,提供了一种全新的思路以及技术来确定渗漏通道的分区区域。使用发电机通过导线向供电电极A和B注入电流信号,电流沿着堤坝渗漏通道产生的优先导电路径引导时,由于通道的存在所产生的电流扰动会产生磁异常,当电流扰动越大,其产生的磁异常越大,同时其梯度变化就越明显。利用磁测梯度进行渗漏通道检测,相比于磁异常数据,其分量更多,数据更多,异常更加明显,同时所需仪器轻便、操作简单、测量效率高,为渗漏检测提供另一种检测手段。
2.本发明进一步的优选方案中,还可以检测出流向渗漏通道的水流流向,本发明具体是根据Maxwell-Ampere定律推导求得磁场梯度中所包含的磁场旋度与传导电流密度之间的数学关系,进一步推导了正六面体磁场梯度计算公式。所以,此发明通过直接计算磁场旋度来获得传导电流密度,,相比于利用磁异常反演电流密度同时也更加高效、快速;同时由于采用了正六面体磁梯度计算使梯度计算容易实施,计算精度更高,在测量中仪器重量可以进一步减小,测量易于实施。而传导电流密度分布与水流方向两者之间的关系在于,当水流由坝坡渗漏入口通过渗漏通道流向下游时,水库中的水流会向渗漏口流动,根据电流场与水流场的相似性,传导电流密度汇聚方向即为水流的流动方向。本发明通过测量磁场的方式得到水中的传导电流密度分布情况,进而利用传导电流密度汇聚方向判断流向渗漏通道的水流流向。相较于利用复杂的方式反演电流密度,此方法能够更加简单快速的获得水中电流密度,从而可以确定流向渗漏通道的水流方向。
附图说明
图1为本发明基于六面体磁梯度装置的堤坝渗漏检测示意图。其中A、B为供电电极,A位于渗漏通道的出水口,B位于堤坝迎水坡一侧;I(ω,t)为供电电流;观测区域位于堤坝上和水库的水中。
图2为正六面体磁梯度测量观测系统,其中1、2、3、4、5、6、7和8分别为三分量磁通门磁传感器,其边长均为d=500mm。
图3为水中的测量区域示意图。其中I为从无穷远到o点电流的大小,z为测点的水深,r为测点到z轴的垂直距离。
图4为水深为50m处,测点沿着与x轴夹角为45°的方向移动时,理论磁场梯度与六面体装置磁场梯度相对误差对比图。
图5为水深为50m处,测点沿着与x轴夹角为45°的方向移动时,理论电流密度与六面体装置磁场电流密度相对误差对比图,其中,Jx和Jy重合;
图6是本发明实施例2提供的一种堤坝渗漏通道检测方法的流程示意图。
具体实施方式
本发明提供了一种堤坝渗漏通道检测方法用于检测堤坝的渗漏情况,其核心是利用磁场梯度判断渗漏通道的分布;此外,还通过直接计算磁场旋度来获得传导电流密度,从而利用传导电流密度汇聚方向判断流向渗漏通道的水流流向。下面将结合实施例对本发明做进一步的说明。
实施例1:
本实施例提供的一种堤坝渗漏通道检测方法,其主要目的在于检测出渗漏通道的分布,具体包括以下几个步骤:
步骤1:在堤坝的两侧布设供电电极A和B,其中,优选供电电极A位于堤身内侧渗漏通道的出水口,供电电极B位于堤坝迎水坡一侧,供电电极通过导线与供电装置相连接。
其中,本实施例还优选供电导线以U型布置,并将供电装置放置较远处,以减小导线磁场对测量数据的影响。
步骤2:供电装置通过供电电极A、B向堤坝和水中供入电流,其电流激励信号为直流信号或者频率为0.1~380Hz的方波或伪随机信号。
步骤3:获取通电时每个测点对应的磁场梯度G。
本实施例中基于每个测点对应的总磁感应强度并采用中心差分方式得到测点的磁场梯度G。具体如下:
将磁场测量装置放置在测点上,使用六面体八个顶点上的磁传感器获得每个传感器上的x、y、z方向的磁场分量,总共得到24个磁场分量;
然后,根据八个磁传感器在x、y、z方向磁场分量并通过中心差分方式得到磁场梯度G,
Figure BDA0003323419750000061
Figure BDA0003323419750000062
Figure BDA0003323419750000071
Figure BDA0003323419750000072
Figure BDA0003323419750000073
Figure BDA0003323419750000074
Figure BDA0003323419750000075
Figure BDA0003323419750000076
Figure BDA0003323419750000077
Figure BDA0003323419750000078
式中,Bx、By、Bz表示在磁观测点上x、y、z方向的磁场分量,
Figure BDA00033234197500000710
表示第i个磁传感器在k方向的磁场分量,d为任意两个磁传感器之间的距离,
Figure BDA0003323419750000079
为偏导数符号,x方向平行堤坝轴向,y方向垂直堤坝并指向堤坝,z方向垂直水面并指向下。
应当理解,本实施例每个测点是利用了8个磁传感器,并使用中心差分方式计算磁场梯度G;其他可行的实施方式中,也可以采用其他手段计算磁场梯度G,本发明对此不进行具体的限定。
步骤4:基于每个测点的磁场梯度G获取渗漏通道的分布,其中,磁场梯度G分布中变化最大的区域为渗漏通道的分布区域。
其中,使用发电机通过导线向供电电极A和B注入电流信号,电流沿着堤坝渗漏通道产生的优先导电路径引导时,由于通道的存在所产生的电流扰动会产生磁异常,当电流扰动越大,其产生的磁异常越大,同时其梯度变化就越明显。通过处理和分析所获得的磁场梯度数据就能判断通道的分布。
实施例2:
本实施例在实施例1的基础上,所述渗漏通道检测方法还可以检测出流向渗漏通道的水流方向,具体是通过如下步骤得到:
获取每个测点上的磁场旋度,再基于所述磁场旋度获得测点上的传导电流密度;
根据所述传导电流密度的汇聚方向得到流向渗漏通道的水流方向,所述汇聚方向与所述水流方向相对应。
其中,所述传导电流密度的计算公式如下;
Figure BDA0003323419750000081
式中,
Figure BDA0003323419750000082
表示磁观测点处的磁场旋度,μ0为自由空间的磁导率,Jc为传导电流密度。即本实施例中利用中心差分方式得到测点的磁场梯度G的参量可以计算出磁场旋度,从而计算出传导电流密度。
为了验证中心差分方式所获得的磁场梯度和电流密度与理论值是否相等,设计了如下模型。
在一个平静的水面上,假设水面为均匀半空间介质,在水面一点上供入电流,并设置该点为坐标原点。则水中的电流密度为:
Figure BDA0003323419750000083
Figure BDA0003323419750000084
Figure BDA0003323419750000085
其中,x、y、z为测点坐标,I为供电电流。则测点处的磁场:
Figure BDA0003323419750000086
Figure BDA0003323419750000087
Bz=0
式中,μ0为自由空间的磁导率,I为供电电流。根据理论推导可到磁梯度公式如下:
Figure BDA0003323419750000091
Figure BDA0003323419750000092
Figure BDA0003323419750000093
Figure BDA0003323419750000094
Figure BDA0003323419750000095
Figure BDA0003323419750000096
Bzx=0
Bzy=0
Bzz=0
则在水深为30m处,测点沿着与x轴夹角为45°的方向移动,共设置测点20个,每个测点间隔5米。通过正演计算可以获得其理论磁场梯度与实际磁场梯度相对误差对比,如图2所示,图中可以看出,通过中心差分方式获得的磁场,与理论磁场误差很小,最大误差约为0.03%,证明中心差分方式获得的数据可靠。
同时,在相同条件下使用中心差分所获得的磁场梯度进而计算电流密度,即:
Figure BDA0003323419750000097
式中,μ0为自由空间的磁导率,Jc为传导电流密度。中心差分计算出的电流密度与理论电流密度的相对误差对比,如图3所示。由图中可知,其最大误差仅为0.03%。由此可以得到通过六面体磁梯度装置可以获得水中的电流密度。
实施例3:
本实施例提供一种基于渗漏通道检测方法的检测装置,其包括:电流装置、测量模块以及检测模块。
本实施例中,电流装置包括供电电极A和B、供电装置,供电电极A位于堤身内侧渗漏通道的出水口,供电电极B位于堤坝迎水坡一侧,供电电极通过导线与供电装置相连接;同时供电导线布置成U型,并将供电装置放置较远处,以减小导线磁场对测量数据的影响。供电装置提供的供电电流信号为直流信号或者频率为0.1~380Hz的方波或伪随机信号。
测量模块设置于每个测点上,即每个测点上设有一个测量模块。本实施例中测量装置为六面体磁探测装置,其包括:六面体支撑件、8个磁传感器、采集箱;
其中,8个磁传感器分别设置在六面体支撑件的8个顶点上,六面体支撑件选用六面体环氧板框架,其各边长度均为500mm,框架底部带有水平装置可以调节仪器达到水平(在水中时通过绳索调平)。8个磁传感器均通过电缆与采集箱连接,采集箱内设有模数转换模块,采集箱采回的数据经过模数转换之后通过USB接口传回电脑终端,并保存。且优选八个磁传感器为磁通门式磁传感器。在测量时,可在堤坝表面测量和船上将仪器放置水中测量,也可将磁梯度测量装置搭载在无人机等航空设备之上,实现磁场矢量及梯度和电流密度的空中快速测量。需要说明的是,采集箱可以是每个测量装置均设有一个,也可以是几个测点的测量装置共享一个,磁传感器与采集箱的通讯方式也据局限于电缆,本发明对此不进行具体的限定。
检测模块用于与所述测量模块连接,用于接收所述测量模块中磁传感器采集的数据,并按照本发明提供的一种堤坝渗漏通道检测方法进行渗透检测。具体是利用每个测点的磁场梯度G获取渗漏通道的分布以及基于所述磁场旋度获得测点上的传导电流密度,再根据所述传导电流密度的汇聚方向得到流向渗漏通道的水流方向。
需要说明的是,检测模块可以是计算机终端,或者是其他具有数据处理分析功能的处理器。
实施例4:
本实施例提供一种检测装置,其包括:
一个或多个处理器;
存储一个或多个程序的存储器;
所述处理器调用所述程序以执行:
获取通电时每个测点对应的磁场梯度G;
基于每个测点的磁场梯度G获取渗漏通道的分布,其中,磁场梯度G分布中变化最大的区域为渗漏通道的分布区域。
在一些实现方式中,程序被处理器调用还执行:
获取每个测点上的磁场旋度,再基于所述磁场旋度获得测点上的传导电流密度;
根据所述传导电流密度的汇聚方向得到流向渗漏通道的水流方向,所述汇聚方向与所述水流方向相对应。
该终端还包括:通信接口,用于与外界设备进行通信,进行数据交互传输。
其中,存储器可能包含高速RAM存储器,也可能还包括非易失性除颤器,例如至少一个磁盘存储器。
如果存储器、处理器和通信接口独立实现,则存储器、处理器和通信接口可以通过总线相互连接并完成相互间的通信。所述总线可以是工业标准体系结构总线,外部设备互联总线或扩展工业标准体系结构总线等。所述总线可以分为地址总线、数据总线、控制总线等。
可选的,在具体实现上,如果存储器、处理器和通信接口集成在一块芯片上,则存储器、处理器即通信接口可以通过内部接口完成相互之间的通信。
各个步骤的具体实现过程请参照前述方法的阐述。
应当理解,在本发明实施例中,所称处理器可以是中央处理单元(CentralProcessing Unit,CPU),该处理器还可以是其他通用处理器、数字信号处理器(DigitalSignal Processor,DSP)、专用集成电路(Application Specific Integrated Circuit,ASIC)、现成可编程门阵列(Field-Programmable GateArray,FPGA)或者其他可编程逻辑器件、分立门或者晶体管逻辑器件、分立硬件组件等。通用处理器可以是微处理器或者该处理器也可以是任何常规的处理器等。存储器可以包括只读存储器和随机存取存储器,并向处理器提供指令和数据。存储器的一部分还可以包括非易失性随机存取存储器。例如,存储器还可以存储设备类型的信息。
实施例5:
本实施例提供一种可读存储介质,其存储了计算机程序,所述计算机程序被处理器调用以执行:
获取通电时每个测点对应的磁场梯度G;
基于每个测点的磁场梯度G获取渗漏通道的分布,其中,磁场梯度G分布中变化最大的区域为渗漏通道的分布区域。
在一些实现方式中,所述计算机程序被处理器调用还执行:
获取每个测点上的磁场旋度,再基于所述磁场旋度获得测点上的传导电流密度;
根据所述传导电流密度的汇聚方向得到流向渗漏通道的水流方向,所述汇聚方向与所述水流方向相对应。
各个步骤的具体实现过程请参照前述方法的阐述。
所述可读存储介质为计算机可读存储介质,其可以是前述任一实施例所述的控制器的内部存储单元,例如控制器的硬盘或内存。所述可读存储介质也可以是所述控制器的外部存储设备,例如所述控制器上配备的插接式硬盘,智能存储卡(Smart Media Card,SMC),安全数字(Secure Digital,SD)卡,闪存卡(Flash Card)等。进一步地,所述可读存储介质还可以既包括所述控制器的内部存储单元也包括外部存储设备。所述可读存储介质用于存储所述计算机程序以及所述控制器所需的其他程序和数据。所述可读存储介质还可以用于暂时地存储已经输出或者将要输出的数据。
基于这样的理解,本发明的技术方案本质上或者说对现有技术做出贡献的部分,或者该技术方案的全部或部分可以以软件产品的形式体现出来,该计算机软件产品存储在一个存储介质中,包括若干指令用以使得一台计算机设备(可以是个人计算机,服务器,或者网络设备等)执行本发明各个实施例所述方法的全部或部分步骤。而前述的可读存储介质包括:U盘、移动硬盘、只读存储器(ROM,Read-OnlyMemory)、随机存取存储器(RAM,RandomAccess Memory)、磁碟或者光盘等各种可以存储程序代码的介质。
需要强调的是,本发明所述的实例是说明性的,而不是限定性的,因此本发明不限于具体实施方式中所述的实例,凡是由本领域技术人员根据本发明的技术方案得出的其他实施方式,不脱离本发明宗旨和范围的,不论是修改还是替换,同样属于本发明的保护范围。

Claims (10)

1.一种基于六面体磁探测装置的渗漏检测方法,其特征在于:包括以下步骤:
步骤1:在堤坝的两侧布设供电电极A和B,供电电极通过导线与供电装置相连接;
步骤2:供电装置通过供电电极A、B向堤坝和水中供入电流;
步骤3:利用六面体磁探测装置获取通电时每个测点对应的磁场梯度G;
步骤4:基于每个测点的磁场梯度G获取渗漏通道的分布,其中,磁场梯度G分布中变化最大的区域为渗漏通道的分布区域。
2.根据权利要求1所述的方法,其特征在于:步骤3中所述每个测点对应的磁场梯度G的获取方式为:基于每个测点对应的总磁感应强度并采用中心差分方式得到测点的磁场梯度G。
3.根据权利要求2所述的方法,其特征在于:每个测点上设有八个磁传感器,利用中心差分方式得到测点的磁场梯度G的公式如下:
Figure FDA0003323419740000011
Figure FDA0003323419740000012
Figure FDA0003323419740000013
Figure FDA0003323419740000014
Figure FDA0003323419740000015
Figure FDA0003323419740000016
Figure FDA0003323419740000017
Figure FDA0003323419740000018
Figure FDA0003323419740000019
Figure FDA00033234197400000110
式中,B为磁场矢量,Bx、By、Bz表示在磁观测点上x、y、z方向的磁场分量,
Figure FDA0003323419740000024
表示第i个磁传感器在k方向的磁场分量,d为任意两个磁传感器之间的距离,
Figure FDA0003323419740000022
为偏导数符号,
Figure FDA0003323419740000023
为磁场梯度。
4.根据权利要求1所述的方法,其特征在于:还包括:
获取每个测点上的磁场旋度,再基于所述磁场旋度获得测点上的传导电流密度;
根据所述传导电流密度的汇聚方向得到流向渗漏通道的水流方向,所述汇聚方向与所述水流方向相对应。
5.根据权利要求3和4所述的方法,其特征在于:所述传导电流密度的计算公式如下;
Figure FDA0003323419740000021
式中,
Figure FDA0003323419740000025
表示磁观测点处的磁场旋度,μ0为自由空间的磁导率,Jc为传导电流密度。
6.根据权利要求3所述的方法,其特征在于:每个测点上的六面体磁探测装置为:八个磁传感器固定在六面体支撑件的八个顶点上,所述六面体支撑件为正六面体,六面体支撑架采用环氧板材料。
7.根据权利要求1所述的方法,其特征在于:所述供电装置所发出的电流激励信号为直流信号或者频率为0.1~380Hz的方波或伪随机信号,供电电极A、B与供电模块之间的供电导线U字排列。
8.一种基于权利要求1-7任一项所述方法的装置,其特征在于:包括:电流装置、测量模块以及检测模块;
其中,所述电流装置包括供电电极A和B、供电装置,供电电极A位于堤身内侧渗漏通道的出水口,供电电极B位于堤坝迎水坡一侧,供电电极通过导线与供电装置相连接;
测量模块为六面体磁探测装置,其设置于每个测点上,其内设有磁传感器;
检测模块与所述测量模块连接,用于接收所述测量模块中磁传感器采集的数据,并进行渗透检测。
9.一种检测装置,其特征在于:包括:
一个或多个处理器;
存储一个或多个程序的存储器;
所述处理器调用所述程序以执行:
获取通电时每个测点对应的磁场梯度G;
基于每个测点的磁场梯度G获取渗漏通道的分布,其中,磁场梯度G分布中变化最大的区域为渗漏通道的分布区域。
10.一种可读存储介质,其特征在于:存储了计算机程序,所述计算机程序被处理器调用以执行:
获取通电时每个测点对应的磁场梯度G;
基于每个测点的磁场梯度G获取渗漏通道的分布,其中,磁场梯度G分布中变化最大的区域为渗漏通道的分布区域。
CN202111254499.5A 2021-10-27 2021-10-27 一种基于六面体磁探测装置的渗漏检测方法及装置 Active CN113960679B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111254499.5A CN113960679B (zh) 2021-10-27 2021-10-27 一种基于六面体磁探测装置的渗漏检测方法及装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111254499.5A CN113960679B (zh) 2021-10-27 2021-10-27 一种基于六面体磁探测装置的渗漏检测方法及装置

Publications (2)

Publication Number Publication Date
CN113960679A true CN113960679A (zh) 2022-01-21
CN113960679B CN113960679B (zh) 2024-01-26

Family

ID=79467529

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111254499.5A Active CN113960679B (zh) 2021-10-27 2021-10-27 一种基于六面体磁探测装置的渗漏检测方法及装置

Country Status (1)

Country Link
CN (1) CN113960679B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114114433A (zh) * 2021-12-08 2022-03-01 黄河勘测规划设计研究院有限公司 全张量梯度的坝体渗漏通道探测装置及方法

Citations (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0981060A2 (de) * 1998-08-20 2000-02-23 Forschungszentrum Jülich Gmbh Verfahren und Vorrichtung zur oberflächennahen Detektion von Stromdichteverteilungen in einem Untergrund
JP2006126108A (ja) * 2004-11-01 2006-05-18 Taisei Corp 漏水検知方法
WO2014196892A1 (en) * 2013-06-04 2014-12-11 Siemens Aktiengesellschaft System for leakage and collapse detection of levees and method using the system
KR101513227B1 (ko) * 2014-10-31 2015-04-17 김용태 비탈면, 댐, 제방 및 교량의 변위 측정을 위한 변위측정 센서모듈 및 이를 이용한 시스템
KR20180013458A (ko) * 2016-07-29 2018-02-07 한국수자원공사 누수 경로 탐지 시스템 및 방법
CN109946370A (zh) * 2019-02-19 2019-06-28 中南大学 一种基于磁场测量的堤坝渗漏通道检测方法及其装置
CN110243559A (zh) * 2019-07-08 2019-09-17 吉林大学 基于流场拟合法的堤坝渗漏探测装置及方法
CN110376275A (zh) * 2019-07-24 2019-10-25 中南大学 一种基于磁场测量的堤坝渗漏连通性检测方法及其装置
CN110632131A (zh) * 2019-10-16 2019-12-31 黄河勘测规划设计研究院有限公司 监测渠道堤防工程渗漏的方法
CN110940725A (zh) * 2019-11-06 2020-03-31 中国地质大学(武汉) 一种堤防内部渗漏通道的探测装置及方法
CN111708088A (zh) * 2020-06-28 2020-09-25 中国矿业大学 一种基于磁梯度张量的瞬变电磁实时动态超前探测方法及系统
CN113514535A (zh) * 2021-06-23 2021-10-19 长江水利委员会长江科学院 可移动式堤坝渗透路径弱磁探测数据采集装置和方法

Patent Citations (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0981060A2 (de) * 1998-08-20 2000-02-23 Forschungszentrum Jülich Gmbh Verfahren und Vorrichtung zur oberflächennahen Detektion von Stromdichteverteilungen in einem Untergrund
JP2006126108A (ja) * 2004-11-01 2006-05-18 Taisei Corp 漏水検知方法
WO2014196892A1 (en) * 2013-06-04 2014-12-11 Siemens Aktiengesellschaft System for leakage and collapse detection of levees and method using the system
KR101513227B1 (ko) * 2014-10-31 2015-04-17 김용태 비탈면, 댐, 제방 및 교량의 변위 측정을 위한 변위측정 센서모듈 및 이를 이용한 시스템
KR20180013458A (ko) * 2016-07-29 2018-02-07 한국수자원공사 누수 경로 탐지 시스템 및 방법
CN109946370A (zh) * 2019-02-19 2019-06-28 中南大学 一种基于磁场测量的堤坝渗漏通道检测方法及其装置
CN110243559A (zh) * 2019-07-08 2019-09-17 吉林大学 基于流场拟合法的堤坝渗漏探测装置及方法
CN110376275A (zh) * 2019-07-24 2019-10-25 中南大学 一种基于磁场测量的堤坝渗漏连通性检测方法及其装置
CN110632131A (zh) * 2019-10-16 2019-12-31 黄河勘测规划设计研究院有限公司 监测渠道堤防工程渗漏的方法
CN110940725A (zh) * 2019-11-06 2020-03-31 中国地质大学(武汉) 一种堤防内部渗漏通道的探测装置及方法
CN111708088A (zh) * 2020-06-28 2020-09-25 中国矿业大学 一种基于磁梯度张量的瞬变电磁实时动态超前探测方法及系统
CN113514535A (zh) * 2021-06-23 2021-10-19 长江水利委员会长江科学院 可移动式堤坝渗透路径弱磁探测数据采集装置和方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
李永丰 等: "《综合探测方法在某堤防渗漏强弱分区中的应用》", no. 1, pages 66 - 71 *
杨天春 等: "《堤坝管涌隐患的电磁勘探法理论分析》", vol. 33, no. 10, pages 14 - 17 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114114433A (zh) * 2021-12-08 2022-03-01 黄河勘测规划设计研究院有限公司 全张量梯度的坝体渗漏通道探测装置及方法
CN114114433B (zh) * 2021-12-08 2023-12-05 黄河勘测规划设计研究院有限公司 全张量梯度的坝体渗漏通道探测装置及方法

Also Published As

Publication number Publication date
CN113960679B (zh) 2024-01-26

Similar Documents

Publication Publication Date Title
US9529110B2 (en) Constructing a reduced order model of an electromagnetic response in a subterranean structure
CN110940725A (zh) 一种堤防内部渗漏通道的探测装置及方法
CN107742015A (zh) 基于任意偶极‑偶极装置的直流激电法三维数值模拟方法
CN110376275B (zh) 一种基于磁场测量的堤坝渗漏连通性检测方法及其装置
JP2010508537A (ja) ストリーマー電気比抵抗探査システム及びこれを用いた河底地盤構造解析方法
US7657391B2 (en) Electromagnetically detecting thin resistive bodies in shallow water and terrestrial environments
CN105279330B (zh) 海面运动舰船湍流尾迹的数值仿真方法
CN109668938B (zh) 堤防渗漏通道三维磁测电阻率探测装置及方法
Lin et al. A para-whole space model for underground magnetic resonance sounding studies
CN113960679A (zh) 一种基于六面体磁探测装置的渗漏检测方法及装置
AU2014201354B2 (en) Systems and methods for measuring water properties in electromagnetic marine surveys
KR20090081146A (ko) 스트리머 전기비저항 탐사 시스템 및 이를 이용한 하저지반구조 해석 방법
CN109557592A (zh) 一种全方位观测的Emn广域电磁法
US20220137249A1 (en) Method of prospecting for three-dimensional bodies using geoelectric tm-polarization techniques
CN111175832A (zh) 一种估算地下水三维流动特性的频域电磁感应勘探方法
US8289025B2 (en) Deriving an electromagnetic field in one direction based on measurement data of one or more sensing elements for measuring an electromagnetic field in another direction
CN115128679A (zh) 一种频率域电磁测深方法、系统及电子设备
CN108089231A (zh) 三维梯度矢量法瞬变电磁仪及其应用
KR20160011357A (ko) 수변 구조물의 파이핑 위치 탐지 시스템 및 그 방법
Ruo et al. The analysis of CSAMT responses of dyke embedded below conductive overburden
Su Magnetic Anomaly Data Detection of Local Marine Geomagnetic Field Model considering Robust Trend Surface Scientific Calculation Algorithm
BR102012030628B1 (pt) Método e aparelho para determinar características de uma região-alvo abaixo de um corpo de água
Vrbancich et al. Testing the limits of AEM bathymetry with a floating TEM system
CN112596108B (zh) Amt剖面探测方法、装置和设备
Huan et al. Test and verification of electric field characteristics in seawater based on ship scaled model

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant