CN110376584B - Water supply pipeline leakage detection method based on ground penetrating radar image characteristic signal identification - Google Patents

Water supply pipeline leakage detection method based on ground penetrating radar image characteristic signal identification Download PDF

Info

Publication number
CN110376584B
CN110376584B CN201910574997.4A CN201910574997A CN110376584B CN 110376584 B CN110376584 B CN 110376584B CN 201910574997 A CN201910574997 A CN 201910574997A CN 110376584 B CN110376584 B CN 110376584B
Authority
CN
China
Prior art keywords
signal
image
radar
water supply
ground penetrating
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.)
Active
Application number
CN201910574997.4A
Other languages
Chinese (zh)
Other versions
CN110376584A (en
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.)
Zhejiang University ZJU
Original Assignee
Zhejiang University ZJU
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 Zhejiang University ZJU filed Critical Zhejiang University ZJU
Priority to CN201910574997.4A priority Critical patent/CN110376584B/en
Publication of CN110376584A publication Critical patent/CN110376584A/en
Application granted granted Critical
Publication of CN110376584B publication Critical patent/CN110376584B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
    • G01S13/88—Radar or analogous systems specially adapted for specific applications
    • G01S13/885—Radar or analogous systems specially adapted for specific applications for ground probing
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
    • G01S13/88—Radar or analogous systems specially adapted for specific applications
    • G01S13/89—Radar or analogous systems specially adapted for specific applications for mapping or imaging
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V3/00—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
    • G01V3/12—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation operating with electromagnetic waves

Landscapes

  • Engineering & Computer Science (AREA)
  • Remote Sensing (AREA)
  • Physics & Mathematics (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Electromagnetism (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Environmental & Geological Engineering (AREA)
  • Geology (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Geophysics (AREA)
  • Radar Systems Or Details Thereof (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

本发明属于市政工程城市供水管道漏损的无损检测领域，公开了一种基于探地雷达图像特征信号识别的供水管道漏损检测方法。步骤一，获取供水管道的探地雷达原始图像数据；步骤二，对获取的原始图像数据进行去直流、零时校正、振幅增强和带通滤波处理；步骤三，结合双曲线信号特征和多次波形成原理以及它们在图像上的位置信息，观察图像中各信号的特征，识别出图像中的管道漏损信号和干扰信号，以对图像中管道漏损信息的准确识别。本发明改变了传统的探地雷达图像分析方式，不但可以识别出雷达图像上的管道漏损信号，也可以排除干扰信号，还可以精确定位漏损的具体位置和深度，从而提高探地雷达在供水管道漏损检测中的效率和准确率。The invention belongs to the field of non-destructive detection of leakage of urban water supply pipelines in municipal engineering, and discloses a leakage detection method of water supply pipelines based on ground penetrating radar image feature signal recognition. Step 1, obtain the original image data of the ground penetrating radar of the water supply pipeline; Step 2, perform DC removal, zero-time correction, amplitude enhancement and band-pass filtering processing on the obtained original image data; Step 3, combine hyperbolic signal characteristics and multiple times. Wave formation principle and their position information on the image, observe the characteristics of each signal in the image, identify the pipeline leakage signal and interference signal in the image, and accurately identify the pipeline leakage information in the image. The invention changes the traditional ground penetrating radar image analysis method, not only can identify the pipeline leakage signal on the radar image, but also can eliminate the interference signal, and can also accurately locate the specific position and depth of the leakage, so as to improve the accuracy of the ground penetrating radar. Efficiency and accuracy in leak detection of water supply pipelines.

Description

Water supply pipeline leakage detection method based on ground penetrating radar image characteristic signal identification
Technical Field
The invention belongs to the field of nondestructive detection of urban water supply pipeline leakage in municipal engineering, and particularly relates to a water supply pipeline leakage detection method based on ground penetrating radar image characteristic signal identification.
Background
As a nondestructive detection method, the ground penetrating radar has the advantages of high efficiency, visual image, difficulty in external interference and the like, and has good application in the fields of archaeology, underground cavity detection, geological forecast and the like, but the application and research in the aspect of leakage detection of municipal water supply pipelines are few. Fig. 1 shows an imaging principle of a ground penetrating radar, which is based on an electromagnetic wave propagation theory and detects a target object by using reflection of high-frequency pulse electromagnetic waves on the premise of differences in electrical properties (electrical conductivity and dielectric constant) of media. The larger the dielectric constant difference between the target object and the soil body is, the larger the reflected power of the target object is, and the clearer the radar imaging of the target object is. Table 1 shows the relative dielectric constant of a common material at 100MHz, and it can be seen from the table that the dielectric constant of water is 81, the dielectric constant of a general soil body is 16, the difference is large, theoretically, when the electromagnetic wave reaches the interface of dry soil and wet soil, the reflection coefficient is large enough, and the radar can receive the electromagnetic wave reflected by the radar, so that a signal generated by leakage of a water supply pipeline appears on a radar image. Therefore, the ground penetrating radar has an advantage in the aspect of detecting the leakage of the water supply pipeline theoretically. However, in actual engineering, due to the diversity and complexity of underground media, interference signals similar to leakage signals always appear on the images, so that the leakage signals on the radar images are difficult to accurately identify, and the application effect of the ground penetrating radar in the aspect of pipeline leakage detection is poor. For example, the casing often contains rocks of varying sizes, each of which produces interference signals on the radar image similar to the leakage signals. The main problems of the identification of the leakage signal in the existing image are as follows:
(1) the images contain a large amount of information, i.e. there are many complex, similar signals in the images;
(2) the image signal has strong explanation subjectivity, namely the theoretical basis of the image signal is insufficient, and the image signal is greatly dependent on the detection personnel with rich experience;
(3) the image signal is underutilized, i.e. the detector usually only focuses on where the signal appears in the image, but does not analyze the signal itself.
Therefore, how to accurately identify the leakage signals on the image from the ground penetrating radar image and eliminate the non-leakage signals on the image becomes a key problem for detecting the pipeline leakage by applying the ground penetrating radar. The problem is solved, whether the pipeline leaks or not and the leakage condition can be detected more efficiently and accurately, and the popularization of the ground penetrating radar in the aspect of pipeline leakage detection is facilitated.
TABLE 1 relative dielectric constant of common substances at 100MHz
Figure BDA0002111862160000011
Figure BDA0002111862160000021
Disclosure of Invention
Aiming at the problems in the prior art, the invention provides a water supply pipeline leakage detection method based on ground penetrating radar image characteristic signal identification.
The invention is realized in such a way, and the method for detecting the leakage of the water supply pipeline based on the ground penetrating radar image characteristic signal identification is characterized by comprising the following steps:
acquiring original image data of a ground penetrating radar of a water supply pipeline;
step two, performing DC removal, zero-time correction, amplitude enhancement and band-pass filtering processing on the acquired original image data;
and thirdly, observing the characteristics of each signal in the image by combining the hyperbolic signal characteristics and the multiple wave forming principle and the position information of the hyperbolic signal characteristics and the multiple waves on the image, and identifying the pipeline leakage signal and the interference signal in the image so as to accurately identify the pipeline leakage information in the image.
Further, the acquisition of the raw image data of the ground penetrating radar in the first step includes the following steps:
(1) according to the known pipeline buried depth and pipe diameter, determining the center frequency of the used radar antenna through a ground penetrating radar horizontal resolution and vertical resolution formula; the relation between the horizontal resolution and the vertical resolution of the ground penetrating radar is shown as the following formula:
Figure BDA0002111862160000022
Figure BDA0002111862160000023
Figure BDA0002111862160000024
wherein: rf、RνRespectively the horizontal resolution and the vertical resolution of the ground penetrating radar; h is the embedding depth of the object; lambda is the wavelength of the ground penetrating radar electromagnetic wave when the ground penetrating radar electromagnetic wave is transmitted in an underground medium; c is the propagation rate of the electromagnetic wave in vacuum; f is probeCenter frequency of the radar antenna; ε is the relative permittivity of the subsurface medium;
(2) and according to the position and the trend of the pipeline, the radar antenna is arranged right above the pipeline and moves along the axial direction of the pipeline to finish the acquisition of the original image data.
Further, the second step is as follows:
firstly, zeroing a direct current component of an acquired radar original image signal by using a gather-DC-Shift module in reflexw one-dimensional filtering of radar general processing software to achieve the purposes of removing direct current and zero point drift;
then, a Move start time module in reflexw static correction is used for selecting a first negative peak value or a positive peak value of a direct wave as a zero-time correction point, and the time of the position is set to be 0;
finally, an Energy decay module in the reflexw gain is used for amplifying the amplitude of the weak signals in depth, and then a bandpassbatterfworth module in the reflexw one-dimensional filtering is used for selecting frequency signals in a specific range, so that high-frequency signals can pass through to the maximum extent, and low-frequency signals are attenuated and suppressed.
Further, the third step is as follows:
firstly, according to the formation process of a hyperbolic signal, the size and the depth position of an object can be calculated by a hyperbolic equation extracted from a radar image; the calculated size or depth position of the object can exclude non-leakage signals which do not meet the actual situation on the size or depth position;
secondly, according to the formation principle of multiple waves, when the electromagnetic waves reach an object for the first time, the first reflection occurs, when the reflected waves return to the ground, a part of energy can penetrate through the ground to be received by the radar for the first time, the other part of energy can reach the object again through the ground reflection to form the second reflection and reception, and the second reflection and reception are sequentially performed until the reflected energy is attenuated to the extent that the radar cannot receive, so that a plurality of vertically distributed hyperbolas generated by the same object can appear on an image, and then interference signals are eliminated;
finally, according to the position relation between the water supply pipeline and the leakage area, the hyperbolic signal of the pipeline leakage on the image is overlapped with the banded signal of the water supply pipeline, so that the banded signal of the pipeline is discontinuous, and the isolated hyperbolic signal in the image can be eliminated.
Further, the step three interference signals include: signals generated by air, stone, plastic and metal existing underground on radar images
Further, when the object is a cylindrical object, the relationship between the hyperbolic equation and the size and depth of the object is shown in the formula:
Figure BDA0002111862160000031
Figure BDA0002111862160000032
Figure BDA0002111862160000041
Figure BDA0002111862160000042
wherein: t is tn、t0At positions x and x for the antenna, respectively0Travel time of the electromagnetic wave; r is the radius of the cylindrical object; v is the propagation velocity of the radar electromagnetic wave in the underground medium; x is the number of1、y1The abscissa and the ordinate are respectively the central point of the hyperbola; p and q are respectively the half-axis length of a hyperbolic real axis and the half-axis length of a virtual axis; h is0Burying the object deeply; and t is the electromagnetic wave travel time when the antenna is right above the object.
In conclusion, the beneficial effects of the invention are as follows: the leakage signal identification method is applied to the ground penetrating radar water supply pipeline image, can identify the pipeline leakage signal on the image, eliminates non-leakage signals such as stones and the like on the image, and further more accurately and efficiently positions the pipeline leakage position. The method changes the analysis mode of the traditional ground penetrating radar image signals in the aspect of leakage signal identification, and can accurately judge non-leakage signals such as pipeline leakage signals and stones and the like according to hyperbolic signal characteristics, a multiple wave forming principle and position information of the hyperbolic signal characteristics and the multiple waves on the image. For example, in fig. 7, the dashed box is a pipe leakage signal, and the solid box is a stone non-leakage signal. Comparing the signal characteristics of the two signals can clearly see the multiple wave phenomenon of the leakage signal and the phenomenon that the water supply pipeline signal is overlapped to make the pipeline strip signal discontinuous, but the non-leakage signal does not have the two characteristics.
Drawings
Fig. 1 is a schematic diagram of a ground penetrating radar imaging principle provided by an embodiment of the present invention.
FIG. 2 is a flow chart of a water supply pipeline leakage detection method based on ground penetrating radar image characteristic signal identification according to an embodiment of the present invention.
FIG. 3 is a schematic diagram of a hyperbolic curve forming process provided by an embodiment of the present invention.
Fig. 4 is a schematic diagram illustrating a multiple forming principle according to an embodiment of the present invention.
FIG. 5 is a schematic view of a water supply pipeline leakage model according to an embodiment of the present invention.
Fig. 6 is a schematic diagram of radar acquisition provided by an embodiment of the present invention.
Fig. 7 is a cross-sectional view of a radar including leakage and non-leakage signals according to an embodiment of the present invention.
Detailed Description
In order to make the objects, technical solutions and advantages of the present invention more apparent, the present invention is further described in detail with reference to the following embodiments. It should be understood that the specific embodiments described herein are merely illustrative of the invention and are not intended to limit the invention.
As shown in fig. 2, the invention provides a method for detecting leakage of a water supply pipeline based on ground penetrating radar image characteristic signal identification, which is characterized by comprising the following steps:
acquiring original image data of a ground penetrating radar of a water supply pipeline;
step two, performing DC removal, zero-time correction, amplitude enhancement and band-pass filtering processing on the acquired original image data;
and thirdly, observing the characteristics of each signal in the image by combining the hyperbolic signal characteristics and the multiple wave forming principle and the position information of the hyperbolic signal characteristics and the multiple waves on the image, and identifying the pipeline leakage signal and the interference signal in the image so as to accurately identify the pipeline leakage information in the image.
According to the method, the characteristic performance of the pipeline leakage signal is summarized by analyzing the formation principle of the pipeline leakage signal in the ground penetrating radar image, and the characteristic performance is distinguished from the common stone non-leakage signal, so that the judgment basis of the pipeline leakage signal is enhanced, and the leakage position of the water supply pipeline is judged more efficiently and accurately.
In a preferred embodiment of the present invention, the obtaining of the ground penetrating radar raw image data comprises:
(1) according to the known pipeline buried depth and pipe diameter, determining the center frequency of the used radar antenna through a ground penetrating radar horizontal resolution and vertical resolution formula; the relation between the horizontal resolution and the vertical resolution of the ground penetrating radar is shown in a formula:
Figure BDA0002111862160000051
Figure BDA0002111862160000052
Figure BDA0002111862160000053
wherein: rf、RνRespectively the horizontal resolution and the vertical resolution of the ground penetrating radar; h is the embedding depth of the object; lambda is the wavelength of the ground penetrating radar electromagnetic wave when the ground penetrating radar electromagnetic wave is transmitted in an underground medium; c is the propagation rate of the electromagnetic wave in vacuum; f is the center frequency of the ground penetrating radar antenna; ε is the relative permittivity of the subsurface medium.
(2) And according to the position and the trend of the pipeline, the radar antenna is arranged right above the pipeline and moves along the axial direction of the pipeline to finish the acquisition of the original image data.
In the preferred embodiment of the invention, the acquired raw image data is subjected to de-DC, zero-time correction, amplitude enhancement and band-pass filtering. The method comprises the following specific steps:
firstly, for an acquired radar original image signal, a direct current component of the signal is set to be zero by using a solenoid-DC-Shift module in reflexw one-dimensional filtering of radar general processing software, so that the purposes of removing direct current and zero point drift are achieved. And then, selecting a first negative peak value or a positive peak value of the direct wave as a zero-time correction point by using a Move start time module in reflexw static correction, and setting the time of the position to be 0. Finally, an Energy decay module in the reflexw gain is used for amplifying the amplitude of the weak signals in depth, and then a bandpassbatterfworth module in the reflexw one-dimensional filtering is used for selecting frequency signals in a specific range, so that high-frequency signals can pass through to the maximum extent, and low-frequency signals are attenuated and suppressed.
In a preferred embodiment of the present invention, the third step is specifically as follows:
first, as shown in fig. 3, according to the formation process of the hyperbolic signal, the hyperbolic equation extracted from the radar image can calculate the size and depth position of the object. Taking a cylindrical object as an example, the relationship between the hyperbolic equation and the size and depth position of the object is shown in the formula:
Figure BDA0002111862160000061
Figure BDA0002111862160000062
Figure BDA0002111862160000063
Figure BDA0002111862160000064
wherein: t is tn、t0At positions x and x for the antenna, respectively0Travel time of the electromagnetic wave; r is the radius of the cylindrical object; v is the propagation velocity of the radar electromagnetic wave in the underground medium; x is the number of1、y1The abscissa and the ordinate are respectively the central point of the hyperbola; p and q are respectively the half-axis length of a hyperbolic real axis and the half-axis length of a virtual axis; h is0Burying the object deeply; and t is the electromagnetic wave travel time when the antenna is right above the object.
The calculated size or depth position of the object can exclude non-leakage signals which do not meet the actual situation on the size or depth position;
next, as shown in fig. 4, according to the formation principle of multiple waves, when an electromagnetic wave reaches an object for the first time, a first reflection occurs, when the reflected wave returns to the ground, a part of energy penetrates through the ground and is received by the radar for the first time, and another part of energy is reflected by the ground and reaches the object again, so as to form a second reflection and reception, and the second reflection and reception are sequentially performed until the reflected energy is attenuated until the radar cannot receive, so that a plurality of vertically distributed hyperbolas generated by the same object appear on an image. Because the dielectric constant of the leakage area is greatly different from the surrounding medium, the energy of the first reflected wave is strong, a remarkable multiple wave phenomenon can be generated on the image, and the dielectric constant of the interfering object such as stone is slightly different from the surrounding medium, the multiple wave phenomenon on the image is not obvious, so that the non-leakage signal of the interfering object such as stone can be eliminated;
finally, according to the position relation between the water supply pipeline and the leakage area, the hyperbolic signal of the pipeline leakage on the image is overlapped with the banded signal of the water supply pipeline, so that the banded signal of the pipeline is discontinuous, and the isolated hyperbolic signal in the image can be eliminated.
The application of the principles of the present invention will now be described in further detail with reference to the accompanying drawings.
As shown in fig. 5 to 7, the method for detecting water supply pipeline leakage based on ground penetrating radar image characteristic signal identification provided by the embodiment of the invention comprises the following steps:
(1) and establishing a water supply pipeline leakage model. The method comprises the following steps: the buried depth of the pipeline is 1.0m, the pipe diameter is 100mm, the pipe is filled with water, a round hole with the diameter of 5mm is formed in the front end of the pipeline, dry sand with uniform particles is filled in the model, and a square concrete block with the side length of about 35cm and the thickness of about 5cm is arranged above the rear end of the pipeline. Fig. 5 is a schematic view of a water supply pipeline leakage model after 3 hours of leakage, in fig. 5, 1 is a square concrete block, 2 is a water supply pipeline, and 3 is a spherical leakage area.
(2) And according to the position and the trend of the pipeline, the radar antenna is arranged right above the pipeline and moves along the axial direction of the pipeline to finish the acquisition of the original image data. Fig. 6 shows data acquisition by radar, wherein 4 in fig. 6 is a water supply pipeline, 5 is a ground penetrating radar antenna, and 6 is a moving path of the antenna.
(3) After the original data are obtained, the radar profile is subjected to basic processing, including direct current removal, zero time correction, amplitude enhancement and band-pass filtering.
(4) Fig. 7 is a cross-sectional view of the ground penetrating radar with pipeline leakage obtained after processing, wherein pipeline leakage signals are shown in a dotted line frame in fig. 7, and interference signals of a disturbing object, namely a stone block, are shown in a solid line frame. Comparing the signal characteristics of the two signals can clearly see the multiple wave phenomenon of the leakage signal and the phenomenon that the water supply pipeline signal is overlapped to make the pipeline strip signal discontinuous, and the interference signal does not have the two characteristics.
The above description is only for the purpose of illustrating the preferred embodiments of the present invention and is not to be construed as limiting the invention, and any modifications, equivalents and improvements made within the spirit and principle of the present invention are intended to be included within the scope of the present invention.

Claims (5)

1.一种基于探地雷达图像特征信号识别的供水管道漏损检测方法，其特征在于，该方法包括：1. a water supply pipeline leakage detection method based on ground penetrating radar image feature signal identification, is characterized in that, the method comprises: 步骤一，获取供水管道的探地雷达原始图像数据；Step 1, obtain the original image data of the ground penetrating radar of the water supply pipeline; 步骤二，对获取的原始图像数据进行去直流、零时校正、振幅增强和带通滤波处理；所述步骤二具体如下：Step 2: Perform DC removal, zero-time correction, amplitude enhancement and band-pass filtering processing on the acquired original image data; the second step is as follows: 首先对采集到的雷达原始图像信号，利用雷达通用处理软件reflexw一维滤波中的Subtract-DC-Shift模块将信号的直流分量置零，达到去直流和去零点漂移的目的；Firstly, for the collected radar original image signal, use the Subtract-DC-Shift module in the one-dimensional filtering of the radar general processing software reflexw to zero the DC component of the signal, so as to achieve the purpose of removing DC and removing zero drift; 而后利用reflexw静校正中的Move start time模块选择直达波的第一负峰值或正峰值作为零时校正点，将该位置的时间设置为0；Then use the Move start time module in the reflexw static correction to select the first negative peak or positive peak of the direct wave as the zero-time correction point, and set the time of this position to 0; 最后利用reflexw增益中的Energy decay模块将深度的微弱信号进行振幅放大，再通过reflexw一维滤波中的bandpassbutterworth模块对特定范围内的频率信号进行选择，最大限度的让高频信号通过，并对低频信号进行衰减和抑制；Finally, the energy decay module in the reflexw gain is used to amplify the amplitude of the weak signal in depth, and then the bandpassbutterworth module in the reflexw one-dimensional filter is used to select the frequency signal within a specific range, so as to maximize the passage of high-frequency signals, and to reduce low-frequency signals. Signal attenuation and suppression; 步骤三，结合双曲线信号特征和多次波形成原理以及它们在图像上的位置信息，观察图像中各信号的特征，识别出图像中的管道漏损信号和干扰信号，以对图像中管道漏损信息的准确识别。Step 3: Combine the characteristics of hyperbolic signals and the formation principle of multiples and their position information on the image, observe the characteristics of each signal in the image, and identify the pipeline leakage signal and interference signal in the image, so as to detect pipeline leakage in the image. Accurate identification of damaged information. 2.如权利要求1所述的基于探地雷达图像特征信号识别的供水管道漏损检测方法，其特征在于，所述步骤一中探地雷达原始图像数据的获取包括如下步骤：2. The method for detecting leakage of water supply pipelines based on GPR image feature signal identification as claimed in claim 1, wherein in the step 1, the acquisition of the original image data of GPR comprises the following steps: (1)根据已知的管道埋深和管径大小，通过探地雷达水平分辨率和垂直分辨率公式确定所使用的雷达天线的中心频率；其中，探地雷达水平分辨率和垂直分辨率的关系见下式：(1) According to the known buried depth and diameter of the pipeline, the center frequency of the radar antenna used is determined by the horizontal resolution and vertical resolution formula of the ground penetrating radar; among them, the horizontal resolution and vertical resolution of the ground penetrating radar The relationship is as follows:
Figure FDA0002955719040000011
Figure FDA0002955719040000011
Figure FDA0002955719040000012
Figure FDA0002955719040000012
Figure FDA0002955719040000013
Figure FDA0002955719040000013
 其中：Rf、Rv分别为探地雷达的水平分辨率和垂直分辨率；h为物体的埋置深度；λ为探地雷达电磁波在地下介质传播时的波长；c为电磁波在真空中的传播速率；f为探地雷达天线的中心频率；ε为地下介质的相对介电常数；Among them: R f and R v are the horizontal resolution and vertical resolution of the ground penetrating radar respectively; h is the buried depth of the object; λ is the wavelength of the ground penetrating radar electromagnetic wave propagating in the underground medium; c is the electromagnetic wave in the vacuum. Propagation rate; f is the center frequency of the ground penetrating radar antenna; ε is the relative permittivity of the underground medium; (2)根据管道位置和走向，将雷达天线置于管道正上方并沿管道轴向方向移动完成原始图像数据的采集。(2) According to the position and direction of the pipeline, the radar antenna is placed directly above the pipeline and moved along the axial direction of the pipeline to complete the acquisition of the original image data. 3.如权利要求1所述的基于探地雷达图像特征信号识别的供水管道漏损检测方法，其特征在于，所述步骤三具体如下：3. The water supply pipeline leakage detection method based on the recognition of ground penetrating radar image feature signal as claimed in claim 1, is characterized in that, described step 3 is as follows: 首先，根据双曲线信号的形成过程，从雷达图像中提取的双曲线方程可以计算出物体的尺寸大小和深度位置；由计算得到的物体尺寸或深度位置可以排除在尺寸大小或深度位置上不符合实际情况的非漏损信号；First, according to the formation process of the hyperbolic signal, the hyperbolic equation extracted from the radar image can calculate the size and depth position of the object; the size or depth position of the object obtained from the calculation can be excluded from the size or depth position that does not conform to the non-missing signal of the actual situation; 其次，根据多次波的形成原理，当电磁波第一次到达物体会发生第一次反射，该反射波在回到地面时一部分能量会穿透地面被雷达第一次接收，另一部分能量会经地面反射再次到达该物体上，形成第二次反射、接收，依次进行直到反射能量衰减至雷达无法接收，从而在图像上会出现由同一物体产生的垂直分布的多个双曲线，进而排除干扰信号；Secondly, according to the formation principle of multiple waves, when the electromagnetic wave reaches the object for the first time, the first reflection will occur. When the reflected wave returns to the ground, a part of the energy will penetrate the ground and be received by the radar for the first time, and the other part of the energy will be received by the radar for the first time. The ground reflection reaches the object again, forming a second reflection and reception, which are carried out in sequence until the reflected energy is attenuated to the point where the radar cannot receive it, so that multiple hyperbolas with vertical distribution generated by the same object will appear on the image, thereby eliminating the interference signal. ; 最后，根据供水管道与漏损区域之间的位置关系，可知图像上管道漏损的双曲线信号应与供水管道的带状信号发生重叠而使得管道带状信号不连续，进而可以排除图像中孤立的双曲线信号。Finally, according to the positional relationship between the water supply pipeline and the leakage area, it can be seen that the hyperbolic signal of pipeline leakage on the image should overlap with the band-shaped signal of the water supply pipeline, so that the pipeline band-shaped signal is discontinuous, which can eliminate the isolation in the image. the hyperbolic signal. 4.如权利要求1所述的基于探地雷达图像特征信号识别的供水管道漏损检测方法，其特征在于，所述步骤三干扰信号包括：地下存在的空气、石块、塑料和金属在雷达图像上产生的信号。4. The method for detecting leakage of water supply pipelines based on GPR image feature signal identification as claimed in claim 1, wherein the interference signal in step 3 comprises: air, stones, plastics and metals existing in the ground in the radar signal generated on the image. 5.如权利要求3所述的基于探地雷达图像特征信号识别的供水管道漏损检测方法，其特征在于，当物体为圆柱形物体时，其双曲线方程与物体尺寸大小、深度位置关系见公式所示：5. The water supply pipeline leakage detection method based on the feature signal recognition of ground penetrating radar image as claimed in claim 3, it is characterized in that, when the object is a cylindrical object, its hyperbolic equation and object size, depth position relationship see The formula shows:
Figure FDA0002955719040000021
Figure FDA0002955719040000021
Figure FDA0002955719040000022
Figure FDA0002955719040000022
Figure FDA0002955719040000023
Figure FDA0002955719040000023
Figure FDA0002955719040000024
Figure FDA0002955719040000024
 其中：tn、t0分别为天线在位置x和x0时电磁波的旅行时间；R为圆柱形物体的半径；v为雷达电磁波在地下介质中的传播速率；x1、y1为分别为双曲线中心点的横坐标和纵坐标；p、q分别为双曲线实轴的半轴长和虚轴的半轴长；h0为物体埋深；t为天线在物体正上方时电磁波旅行时间。Among them: t n , t 0 are the travel time of the electromagnetic wave when the antenna is at the position x and x 0 respectively; R is the radius of the cylindrical object; v is the propagation speed of the radar electromagnetic wave in the underground medium; x 1 , y 1 are respectively The abscissa and ordinate of the center point of the hyperbola; p and q are the semi-axis length of the real axis of the hyperbola and the semi-axis length of the imaginary axis, respectively; h 0 is the buried depth of the object; t is the electromagnetic wave travel time when the antenna is directly above the object .
CN201910574997.4A 2019-06-28 2019-06-28 Water supply pipeline leakage detection method based on ground penetrating radar image characteristic signal identification Active CN110376584B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910574997.4A CN110376584B (en) 2019-06-28 2019-06-28 Water supply pipeline leakage detection method based on ground penetrating radar image characteristic signal identification

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910574997.4A CN110376584B (en) 2019-06-28 2019-06-28 Water supply pipeline leakage detection method based on ground penetrating radar image characteristic signal identification

Publications (2)

Publication Number Publication Date
CN110376584A CN110376584A (en) 2019-10-25
CN110376584B true CN110376584B (en) 2021-09-14

Family

ID=68251255

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910574997.4A Active CN110376584B (en) 2019-06-28 2019-06-28 Water supply pipeline leakage detection method based on ground penetrating radar image characteristic signal identification

Country Status (1)

Country Link
CN (1) CN110376584B (en)

Families Citing this family (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112285790B (en) 2020-03-30 2021-06-18 中国科学院地质与地球物理研究所 Method and device for determining attenuation coefficient of electromagnetic wave field and medium
CN111551927B (en) * 2020-05-19 2023-01-31 上海圭目机器人有限公司 Underground pipeline diameter measuring method based on three-dimensional ground penetrating radar
CN111812609B (en) * 2020-06-28 2023-12-01 中国地质调查局西安地质调查中心(西北地质科技创新中心) Geological radar signal recovery method with multistage filtering nested amplitude gain
CN111853557B (en) * 2020-07-28 2021-06-08 中国水利水电科学研究院 An automatic identification method and system for detecting pipeline leakage based on ground penetrating radar
CN112180452B (en) * 2020-09-23 2023-09-29 中国建筑第八工程局有限公司 Underground pipeline buried depth estimation method based on ground penetrating radar and three-dimensional velocity spectrum
CN112346049A (en) * 2020-09-28 2021-02-09 浙江大学 Pipeline leakage detection method and device based on three-dimensional image energy density attribute of ground penetrating radar
CN113688692A (en) * 2021-07-29 2021-11-23 浙江大学 Water supply pipeline leakage detection method based on time-frequency scale characteristics of ground penetrating radar
CN116360001A (en) * 2023-04-04 2023-06-30 中国地质科学院地球物理地球化学勘查研究所 Method for eliminating low-frequency unshielded ground penetrating radar multiple signal interference in frozen soil region
CN119164336A (en) * 2024-10-15 2024-12-20 山东聚合管业有限公司 A deformation detection method and device for PE pure water diversion pipe
CN119987384A (en) * 2025-04-16 2025-05-13 江苏可一网络有限公司 A radar vehicle control system for surveying underground urban roads

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109031290A (en) * 2018-07-27 2018-12-18 山东农业大学 A kind of quality determining method of land engineering
CN109188383A (en) * 2018-08-16 2019-01-11 广西有色勘察设计研究院 Rich water gravel layer shield duct piece grouting behind shaft or drift lining geologic radar detection scheme
JP2019015699A (en) * 2017-07-11 2019-01-31 富士通株式会社 Ranging device, water level measuring system, and ranging method
CN109343022A (en) * 2016-12-30 2019-02-15 北京师范大学 Methods of Estimating Soil Moisture Content Between Layers

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1808087A (en) * 2006-01-27 2006-07-26 重庆大学 Digital and intelligent method for positioning and detecting water leakage of underground pipeline and instrument therefor
CN103675922B (en) * 2013-12-13 2016-03-30 南京工业大学 Operation period underground pipeline pipe diameter measuring method based on ground penetrating radar
US20160281498A1 (en) * 2014-02-18 2016-09-29 Halliburton Energy Services, Inc. Methods for obtaining data from a subterranean formation
US20160370491A1 (en) * 2015-06-21 2016-12-22 Industrea Mining Technology Pty Limited Drill sensor system and method
JP6880850B2 (en) * 2017-03-13 2021-06-02 富士通株式会社 Distance measuring device, water level measurement system and distance measuring method
US10234552B1 (en) * 2018-06-27 2019-03-19 University Of South Florida Precise infrastructure mapping using full-waveform inversion of ground penetrating radar signals

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109343022A (en) * 2016-12-30 2019-02-15 北京师范大学 Methods of Estimating Soil Moisture Content Between Layers
JP2019015699A (en) * 2017-07-11 2019-01-31 富士通株式会社 Ranging device, water level measuring system, and ranging method
CN109031290A (en) * 2018-07-27 2018-12-18 山东农业大学 A kind of quality determining method of land engineering
CN109188383A (en) * 2018-08-16 2019-01-11 广西有色勘察设计研究院 Rich water gravel layer shield duct piece grouting behind shaft or drift lining geologic radar detection scheme

Also Published As

Publication number Publication date
CN110376584A (en) 2019-10-25

Similar Documents

Publication Publication Date Title
CN110376584B (en) Water supply pipeline leakage detection method based on ground penetrating radar image characteristic signal identification
CN104020495B (en) Automatic underground pipeline parameter recognizing method on basis of ground penetrating radar
US11900630B2 (en) Method for detecting leakage of water supply pipe based on ground-penetrating radar three-dimensional image attribute analysis
CN106446919B (en) A fast detection method for ground-penetrating radar hyperbolic targets
CN109343022B (en) Method for estimating interlayer soil water content
CN109387834B (en) Dam damage detection method for underground reservoir and electronic equipment
CN107402386B (en) A method for estimating the radius and buried depth of underground metal circular pipes based on BP neural network
Bi et al. Multi-frequency GPR data fusion and its application in NDT
CN108333634A (en) A kind of Ground Penetrating Radar buried target localization method based on frequency spectrum remnants conspicuousness detection methods
CN108845317B (en) Frequency domain reverse time migration method based on layered medium Green function
De Coster et al. Fusion of multifrequency GPR data freed from antenna effects
Yang et al. Small‐Scale Void‐Size Determination in Reinforced Concrete Using GPR
CN111929678A (en) Method for measuring water content of shallow surface soil
CN108614307B (en) Inverted arch quality detection method and system by using comprehensive geophysical prospecting technology
CN112180452B (en) Underground pipeline buried depth estimation method based on ground penetrating radar and three-dimensional velocity spectrum
Lu et al. Ground penetrating radar applications in mapping underground utilities
CN108646229A (en) Underground column reflector inclination angle detection method
Gharamohammadi et al. A robust whitening algorithm to identify buried objects with similar attributes in correlation-based detection
Alshamy et al. Evaluation of GPR Detection for buried objects material with different depths and scanning angles
CN108303745A (en) A kind of inversion method of the buried cable detection based on electromagnetic wave saturating ground technology
CN107356985A (en) The method that earthen ruins archaeological investigation is carried out using GPR wave impedance inversion
JP2528148B2 (en) Method and device for detecting underground buried objects
Li et al. Distorted hyperbola extraction of GPR data using edge detection for buried non-metallic pipes
CN112346049A (en) Pipeline leakage detection method and device based on three-dimensional image energy density attribute of ground penetrating radar
CN115013739B (en) A method for locating underground water pipe leakage based on ground penetrating radar and wavelet entropy

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