CN111007572A - Automatic identification method, device and system for road underground cavity - Google Patents

Automatic identification method, device and system for road underground cavity Download PDF

Info

Publication number
CN111007572A
CN111007572A CN201911163987.8A CN201911163987A CN111007572A CN 111007572 A CN111007572 A CN 111007572A CN 201911163987 A CN201911163987 A CN 201911163987A CN 111007572 A CN111007572 A CN 111007572A
Authority
CN
China
Prior art keywords
reflected wave
area
reflection coefficient
amplitude
value
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
CN201911163987.8A
Other languages
Chinese (zh)
Other versions
CN111007572B (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.)
Beijing Zhongke blueprints Technology Co.,Ltd.
Zhongke yuntu Technology Co., Ltd
Original Assignee
Beijing Zhongke Blueprints Technology 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 Beijing Zhongke Blueprints Technology Co Ltd filed Critical Beijing Zhongke Blueprints Technology Co Ltd
Priority to CN201911163987.8A priority Critical patent/CN111007572B/en
Publication of CN111007572A publication Critical patent/CN111007572A/en
Application granted granted Critical
Publication of CN111007572B publication Critical patent/CN111007572B/en
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/12Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation operating with electromagnetic waves
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Systems 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/88Radar or analogous systems specially adapted for specific applications
    • G01S13/885Radar or analogous systems specially adapted for specific applications for ground probing

Landscapes

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

Abstract

The invention provides an automatic identification method, a device and a system for road underground cavities, wherein the method comprises the following steps: acquiring sensing data of the ground penetrating radar of each area on the road, wherein the sensing data is used for representing the amplitude of a reflected wave based on the ground penetrating radar; calculating the reflected wave energy value of each area according to the sensing data of the ground penetrating radar, and determining a first area of which the reflected wave energy value is greater than a first preset threshold value; calculating a reflection coefficient of a reflected wave in the first region and a Fourier spectrum amplitude of the reflected wave according to the sensing data of the first region, and further calculating a cavity response value of the reflected wave in the first region according to the reflection coefficient of the reflected wave in the first region and the Fourier spectrum amplitude of the reflected wave; and judging the position corresponding to the first area with the cavity response value larger than the second preset threshold value as an underground cavity. The method can improve the accuracy and reliability of the detection of the underground cavities of the road.

Description

Automatic identification method, device and system for road underground cavity
Technical Field
The invention relates to the field of road detection, in particular to an automatic identification method, device and system for underground cavities of roads.
Background
Due to vehicle vibration, road surface water seepage, underground pipeline water seepage and other reasons, the urban road frequently has the problems of cracking, deformation, settlement, collapse and the like. By means of a geophysical prospecting method, urban road detection is carried out regularly, and hidden danger of road collapse can be early warned in advance. The ground penetrating radar method is a road detection geophysical prospecting method which can simultaneously meet the requirements of rapidness, no damage and high resolution.
At present, radar road detection data are interpreted manually, and in the face of mass data generated by rapid road detection, the manual method is low in efficiency and cannot submit an interpretation result in time, interpretation accuracy rates of different personnel are greatly different, and reliability is low.
Disclosure of Invention
The invention aims to provide a method, a device and a system for automatically identifying underground road cavities, which aim to solve the problems of low accuracy and low reliability of the manual interpretation of radar road detection data in the prior art.
According to a first aspect of the invention, a method for automatically identifying a road underground cavity comprises the following steps: acquiring sensing data of the ground penetrating radar of each area on the road; the sensing data is used for representing the amplitude of a ground penetrating radar-based reflected wave; calculating the reflected wave energy value of each area according to the sensing data of the ground penetrating radar, and determining a first area of which the reflected wave energy value is greater than a first preset threshold value; calculating a reflection coefficient of a reflected wave in the first region and a Fourier spectrum amplitude of the reflected wave according to the sensing data of the first region, and further calculating a cavity response value of the reflected wave in the first region according to the reflection coefficient of the reflected wave in the first region and the Fourier spectrum amplitude of the reflected wave; and judging the position corresponding to the first area with the cavity response value larger than the second preset threshold value as an underground cavity.
Further, the step of calculating the cavity response value of the reflected wave in the first region according to the reflection coefficient of the reflected wave in the first region and the fourier spectrum amplitude of the reflected wave includes:
and calculating the void response value of the reflected wave in the first area according to the weighted summation of the reflection coefficient of the reflected wave in the first area and the Fourier spectrum amplitude of the reflected wave.
Further, before the step of calculating the hole response value of the reflected wave in the first region according to the weighted summation of the reflection coefficient of the reflected wave in the first region and the fourier spectrum amplitude of the reflected wave, the method comprises:
calculating the reflection coefficient of the reflected wave in each first area, and calculating the average value of the reflection coefficients of the reflected wave in each first area;
determining areas with the reflection coefficients larger than the average value in each first area as second areas, wherein the first areas comprise the second areas;
calculating a cavity response value of the reflected wave in the first area according to the weighted summation of the reflection coefficient of the reflected wave in the first area and the Fourier spectrum amplitude of the reflected wave, and determining the position corresponding to the first area with the cavity response value larger than a second preset threshold value as an underground cavity, wherein the step comprises the following steps:
and calculating the cavity response value of the reflected wave in the second area according to the weighted summation of the reflection coefficient of the reflected wave in the second area and the Fourier spectrum amplitude of the reflected wave, and judging the position corresponding to the second area with the cavity response value being greater than a second preset threshold value as an underground cavity.
Further, the step of acquiring sensing data of the ground penetrating radar is followed by: and filtering the sensing data of the acquired ground penetrating radar.
Further, the reflected wave energy value E is calculated as follows:
Figure BDA0002285420820000021
wherein N is the number of sampling points, PiThe amplitude of the reflected wave of the ith sampling point is obtained;
the reflection coefficient R is calculated as follows:
Figure BDA0002285420820000031
wherein, P0The peak value of the direct wave; the fourier spectrum amplitude F (ω) of the reflected wave is calculated as follows:
Figure BDA0002285420820000032
f (t) is the time domain signal of the amplitude of the reflected wave, ω is the amplitude signal frequency of the reflected wave, and t is the amplitude signal of the reflected waveSampling time;
the hole response value K of the reflected wave is calculated in the following way that K is α R + β F (omega), α is a preset reflection coefficient weight, and β is a preset Fourier spectrum amplitude weight.
According to a second aspect of the present invention, an apparatus for automatically identifying a road underground cavity, comprises: the system comprises an acquisition module, a detection module and a processing module, wherein the acquisition module is used for acquiring sensing data of the ground penetrating radar in each area on a road, and the sensing data is used for representing the amplitude of a reflected wave based on the ground penetrating radar;
the processing module is used for calculating the reflected wave energy value of each area according to the sensing data of the ground penetrating radar and determining a first area with the reflected wave energy value being larger than a first preset threshold value; calculating a reflection coefficient of a reflected wave in the first region and a Fourier spectrum amplitude of the reflected wave according to the sensing data of the first region, and further calculating a cavity response value of the reflected wave in the first region according to the reflection coefficient of the reflected wave in the first region and the Fourier spectrum amplitude of the reflected wave;
and the identification module is used for judging the position corresponding to the first area with the cavity response value larger than a second preset threshold value as an underground cavity.
Further, the processing module is further configured to calculate a void response value of the reflected wave in the first region according to a weighted summation of a reflection coefficient of the reflected wave in the first region and a fourier spectrum amplitude of the reflected wave.
Further, the processing module is further configured to calculate a reflection coefficient of the reflected wave in each first area, and calculate an average value of the reflection coefficients of the reflected waves in each first area; determining an area with a reflection coefficient larger than the average value in each first area as a second area; calculating a hole response value of the reflected wave in the second area according to the weighted summation of the reflection coefficient of the reflected wave in the second area and the Fourier spectrum amplitude of the reflected wave, and judging the position corresponding to the second area with the hole response value larger than a second preset threshold value as an underground hole; wherein the first region comprises a second region.
Further, the automatic identification device for the underground cavity of the road further comprises: the filtering module is used for filtering the sensing data of the acquired ground penetrating radar;
the reflected wave energy value E is calculated as follows:
Figure BDA0002285420820000041
wherein N is the number of sampling points, PiThe amplitude of the reflected wave of the ith sampling point is obtained;
the reflection coefficient R is calculated as follows:
Figure BDA0002285420820000042
wherein, P0The peak value of the direct wave; the fourier spectrum amplitude F (ω) of the reflected wave is calculated as follows:
Figure BDA0002285420820000043
f (t) is a time domain signal of the amplitude of the reflected wave, omega is the amplitude signal frequency of the reflected wave, and t is the amplitude signal sampling time of the reflected wave;
the hole response value K of the reflected wave is calculated in the following way that K is α R + β F (omega), α is a preset reflection coefficient weight, and β is a preset Fourier spectrum amplitude weight.
According to a third aspect of the invention, the system for automatically identifying the road underground cavity comprises the device for automatically identifying the road underground cavity and a ground penetrating radar, wherein the device for automatically identifying the road underground cavity is in communication connection with the ground penetrating radar.
According to the method, the device and the system for automatically identifying the road underground cavity, the first area with the energy accumulated value of the ground penetrating radar sensing data larger than the first preset threshold value is determined, the cavity response value of the first area is obtained through calculation according to the reflection coefficient of the first area and the Fourier spectrum amplitude, the position corresponding to the first area with the cavity response value larger than the second preset threshold value is determined as the underground cavity, and the accuracy and the reliability of radar road detection data interpretation are improved based on the analysis of the difference of the physical properties of the cavity and the surrounding medium at the interface of the underground cavity.
Other characteristic features and advantages of the invention will become apparent from the following description of exemplary embodiments, which is to be read in connection with the accompanying drawings.
Drawings
The accompanying drawings, which are incorporated in and constitute a part of the specification, illustrate embodiments of the invention and together with the description, serve to explain the principles of the invention. In the drawings, like reference numerals are used to indicate like elements. The drawings in the following description are directed to some, but not all embodiments of the invention. For a person skilled in the art, other figures can be derived from these figures without inventive effort.
FIG. 1 is a flow chart of an embodiment of a method for automatically identifying road underground cavities according to the present invention;
FIG. 2 is a flow chart of another embodiment of the method for automatically identifying underground cavities in roads of the present invention;
FIG. 3 is a schematic diagram of features of the holes identified by the automatic identification method shown in FIG. 2 in a radar image;
FIG. 4 is a schematic diagram of a Fourier magnitude spectrum of a region shown by a dashed box at a hole position in FIG. 3, wherein a peak between an abscissa 200 and an abscissa 400 is a characteristic of the hole in the Fourier spectrum;
the dotted line in FIG. 5 marks the features of the homogeneous medium in the radar image;
FIG. 6 is a Fourier magnitude spectrum at the dashed line of the uniform medium in FIG. 5, wherein the peaks between the abscissa 200 and the abscissa 400 are the characteristic of the non-voids in the Fourier spectrum;
FIG. 7 is a block diagram of an embodiment of an automatic road underground cavity recognition device according to the present invention;
FIG. 8 is a block diagram of an embodiment of an automatic identification system for underground cavities in a road according to the present invention.
Detailed Description
In order to make the objects, technical solutions and advantages of the embodiments of the present invention clearer, the technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are some, but not all, embodiments of the present invention. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention. It should be noted that the embodiments and features of the embodiments in the present application may be arbitrarily combined with each other without conflict.
The research of the invention finds that: the difference of the physical properties of the cavity and the surrounding medium at the interface of the underground cavity is huge, strong reflection is generated when electromagnetic waves are transmitted, the strong reflection is mainly based on high-frequency components, and a high-frequency peak value can be formed in a frequency spectrum. When the electromagnetic wave is transmitted in a uniform medium, strong reflection cannot be generated, and corresponding high-frequency response cannot exist in a frequency spectrum. The difference between the electromagnetic physical properties of the cavity and the surrounding medium is huge, the reflection coefficient obtained through calculation is larger than that of a normal area, and meanwhile the underground cavity of the road can be accurately identified by combining Fourier frequency spectrum. The following detailed description is made in conjunction with the embodiments illustrated in the various figures.
As shown in FIG. 1, the invention relates to an automatic identification method of road underground cavities, which comprises the following steps:
step 101: acquiring sensing data of the ground penetrating radar of each area on the road, wherein the sensing data is used for representing the amplitude of a reflected wave based on the ground penetrating radar;
step 102: calculating the reflected wave energy value of each area according to the sensing data of the ground penetrating radar, and determining a first area of which the reflected wave energy value is greater than a first preset threshold value; calculating a reflection coefficient of a reflected wave in the first region and a Fourier spectrum amplitude of the reflected wave according to the sensing data of the first region, and further calculating a cavity response value of the reflected wave in the first region according to the reflection coefficient of the reflected wave in the first region and the Fourier spectrum amplitude of the reflected wave;
step 103: and judging the position corresponding to the first area with the cavity response value larger than the second preset threshold value as an underground cavity.
According to the embodiment of the automatic identification method of the road underground cavity, the first area with the energy accumulated value of the ground penetrating radar sensing data larger than the first preset threshold value is determined, the cavity response value of the first area is obtained through calculation according to the reflection coefficient of the first area and the Fourier spectrum amplitude, the position corresponding to the first area with the cavity response value larger than the second preset threshold value is determined as the underground cavity, and the accuracy and the reliability of radar road detection data interpretation are improved based on the analysis of the difference of the physical properties of the cavity and the surrounding medium at the interface of the underground cavity.
As shown in fig. 2, the present invention provides an automatic identification method for underground road cavities, which is a preferred implementation of the embodiment of the method shown in fig. 1, and the explanation of the embodiment shown in fig. 1 can be applied to this embodiment, and the automatic identification method of this embodiment includes:
step 201: carrying out data loading to acquire sensing data of the ground penetrating radar;
step 202: filtering the ground penetrating radar data to remove interference clutter; the specific filtering manner can be various, such as mean filtering, FK filtering (Frequency Wave-number domain filtering), etc.;
step 203: calculating an energy accumulated value for the filtered data; the energy accumulation value can be calculated in various ways, such as summing the squares of the amplitude values of all points of each data track; the reflected wave energy value E is calculated as follows:
Figure BDA0002285420820000071
wherein N is the number of sampling points, PiThe amplitude of the reflected wave of the ith sampling point is obtained;
step 204: determining a region with an energy value smaller than or equal to a first preset threshold as a non-cavity region; the size of the first preset threshold is set according to actual needs, and is not limited herein;
step 205: determining an area with an energy value larger than a first preset threshold value as a first area, and calculating a reflection coefficient; calculating the average value of the reflection coefficients of all the first areas; the reflection coefficient refers to the ratio of the amplitude value domain direct wave amplitude value of each point of each channel of data; such asThe reflection coefficient R is calculated as follows:
Figure BDA0002285420820000072
wherein, P0The peak value of the direct wave;
step 206: determining a first area with a reflection coefficient less than or equal to the average value of the reflection coefficients as a non-void area;
step 207: calculating a fourier spectrum of a region (as a second region, the first region includes the second region) in which the reflection coefficient is larger than the average value; for the area (namely a second area) with the reflection coefficient larger than the average value, carrying out weighted summation on the reflection coefficient and the Fourier spectrum amplitude to obtain a void response value; such as: the fourier spectrum magnitude F (ω) of the reflected wave is calculated as follows:
Figure BDA0002285420820000073
f (t) is a time domain signal of the amplitude of the reflected wave, omega is the amplitude signal frequency of the reflected wave, and t is the amplitude signal sampling time of the reflected wave;
the calculation method of the hole response value K of the reflected wave is as follows, wherein K is α R + β F (omega), α is a preset reflection coefficient weight, β is a preset Fourier spectrum amplitude weight, the specific values of α and β can be set according to actual needs, including but not limited to the fact that the weights of the two are both 1, and unequal weights can be selected according to actual experimental conditions.
Step 208: judging the position of the void response value less than or equal to a second preset threshold value as a non-underground void; the size of the second preset threshold is set according to actual needs, and is not limited herein;
step 209: and judging the position with the cavity response value larger than the second preset threshold value as an underground cavity.
The experimental results of the automatic identification method for the road underground cavity in the embodiment can be seen in fig. 3-6, and it can be seen that the accuracy and reliability of interpretation of radar road detection data can be effectively improved by the automatic identification method for the road underground cavity in the embodiment.
According to the embodiment of the automatic identification method of the road underground cavity, after the first area is determined, the reflection coefficient of each data is calculated for the whole radar data, the size of the reflection coefficient of each data is compared with the average value of the reflection coefficients, the second area with the reflection coefficient larger than the average value of the reflection coefficients is defined, FFT operation is conducted on the data of the second area, the Fourier frequency spectrum of each data is calculated, weighted summation calculation is conducted on the reflection coefficient and the Fourier frequency spectrum, the cavity response value is obtained, the position with the cavity response value larger than the threshold value can be determined as the underground cavity area, and the accuracy and reliability of radar road detection data interpretation can be effectively improved.
As shown in fig. 7, the present invention further provides an automatic identification device for underground road cavities, which is a corresponding device embodiment to the method embodiments shown in fig. 1 and 2, and the explanation of the embodiments shown in fig. 1 to 6 can be applied to this embodiment, and the automatic identification device includes:
an obtaining module 701, configured to obtain sensing data of a ground penetrating radar in each area on a road, where the sensing data is used to characterize an amplitude of a reflected wave based on the ground penetrating radar;
the processing module 702 is configured to calculate a reflected wave energy value of each region according to the sensing data of the ground penetrating radar, and determine a first region in which the reflected wave energy value is greater than a first preset threshold; calculating a reflection coefficient of a reflected wave in the first region and a Fourier spectrum amplitude of the reflected wave according to the sensing data of the first region, and further calculating a cavity response value of the reflected wave in the first region according to the reflection coefficient of the reflected wave in the first region and the Fourier spectrum amplitude of the reflected wave;
the identifying module 703 is configured to determine a position corresponding to the first area where the cavity response value is greater than the second preset threshold as an underground cavity.
Further, the processing module 702 is further configured to calculate a void response value of the first region according to a weighted summation of the reflection coefficient of the first region and the fourier spectrum amplitude.
Further, the processing module 702 is further configured to calculate a reflection coefficient of the reflected wave in each first area, and calculate an average value of the reflection coefficients of the reflected waves in each first area; determining an area with a reflection coefficient larger than the average value in each first area as a second area; calculating a hole response value of the reflected wave in the second area according to the weighted summation of the reflection coefficient of the reflected wave in the second area and the Fourier spectrum amplitude of the reflected wave, and judging the position corresponding to the second area with the hole response value larger than a second preset threshold value as an underground hole; wherein the first region comprises a second region.
Further, the automatic identification device for the underground road cavity may further include: and a filtering module 704, configured to perform filtering processing on the sensing data of the ground penetrating radar.
As shown in fig. 8, the present invention further provides an automatic identification system for underground road cavities, which comprises the automatic identification device shown in fig. 7. The explanation of the embodiment shown in fig. 1-7 can be applied to this embodiment, and the automatic identification system of this embodiment includes the automatic identification device for the underground cavity of the road and the ground penetrating radar, and the automatic identification device is connected to the ground penetrating radar in communication.
The automatic identification system for the road underground cavity has the corresponding technical effects of the automatic identification device for the road underground cavity, and the details are not repeated herein.
The above-described aspects may be implemented individually or in various combinations, and such variations are within the scope of the present invention.
Finally, it should be noted that: the above examples are only for illustrating the technical solutions of the present invention, and are not limited thereto. Although the present invention has been described in detail with reference to the foregoing embodiments, it will be understood by those of ordinary skill in the art that: the technical solutions described in the foregoing embodiments may still be modified, or some technical features may be equivalently replaced; and such modifications or substitutions do not depart from the spirit and scope of the corresponding technical solutions of the embodiments of the present invention.

Claims (10)

1. An automatic identification method for underground road cavities is characterized by comprising the following steps:
acquiring sensing data of the ground penetrating radar of each area on the road, wherein the sensing data is used for representing the amplitude of a reflected wave based on the ground penetrating radar;
calculating the reflected wave energy value of each area according to the sensing data of the ground penetrating radar, and determining a first area of which the reflected wave energy value is greater than a first preset threshold value;
calculating a reflection coefficient of a reflected wave in the first region and a Fourier spectrum amplitude of the reflected wave according to the sensing data of the first region, and further calculating a cavity response value of the reflected wave in the first region according to the reflection coefficient of the reflected wave in the first region and the Fourier spectrum amplitude of the reflected wave;
and judging the position corresponding to the first area with the cavity response value larger than the second preset threshold value as an underground cavity.
2. The method for automatically identifying an underground cavity of a road according to claim 1, wherein: calculating a hole response value of the reflected wave in the first area according to the reflection coefficient of the reflected wave in the first area and the Fourier spectrum amplitude of the reflected wave, wherein the step comprises the following steps:
and calculating the void response value of the reflected wave in the first area according to the weighted summation of the reflection coefficient of the reflected wave in the first area and the Fourier spectrum amplitude of the reflected wave.
3. The method for automatically identifying an underground cavity of a road according to claim 2, wherein: the method comprises the following steps of calculating the void response value of the reflected wave in the first area according to the weighted summation of the reflection coefficient of the reflected wave in the first area and the Fourier spectrum amplitude of the reflected wave, wherein the method comprises the following steps:
calculating the reflection coefficient of the reflected wave in each first area, and calculating the average value of the reflection coefficients of the reflected wave in each first area;
determining areas with the reflection coefficients larger than the average value in each first area as second areas, wherein the first areas comprise the second areas;
calculating a cavity response value of the reflected wave in the first area according to the weighted summation of the reflection coefficient of the reflected wave in the first area and the Fourier spectrum amplitude of the reflected wave, and determining the position corresponding to the first area with the cavity response value larger than a second preset threshold value as an underground cavity, wherein the step comprises the following steps:
and calculating the cavity response value of the reflected wave in the second area according to the weighted summation of the reflection coefficient of the reflected wave in the second area and the Fourier spectrum amplitude of the reflected wave, and judging the position corresponding to the second area with the cavity response value being greater than a second preset threshold value as an underground cavity.
4. The method of claim 3, wherein the step of obtaining sensed data of the ground penetrating radar of each area on the road is followed by:
and filtering the sensing data of the acquired ground penetrating radar.
5. The method for automatically identifying an underground cavity of a road according to any one of claims 1 to 4, wherein:
the reflected wave energy value E is calculated as follows:
Figure FDA0002285420810000023
wherein N is the number of sampling points, PiThe amplitude of the reflected wave of the ith sampling point is obtained;
the reflection coefficient R is calculated as follows:
Figure FDA0002285420810000021
wherein, P0The peak value of the direct wave; the fourier spectrum amplitude F (ω) of the reflected wave is calculated as follows:
Figure FDA0002285420810000022
f (t) is the time domain signal of the amplitude of the reflected wave, ω is the amplitude of the reflected waveThe signal frequency t is the amplitude signal sampling time of the reflected wave;
the hole response value K of the reflected wave is calculated in the following way that K is α R + β F (omega), α is a preset reflection coefficient weight, and β is a preset Fourier spectrum amplitude weight.
6. An automatic identification device of road underground cavity, characterized by, includes:
the system comprises an acquisition module, a detection module and a processing module, wherein the acquisition module is used for acquiring sensing data of the ground penetrating radar in each area on a road, and the sensing data is used for representing the amplitude of a reflected wave based on the ground penetrating radar;
the processing module is used for calculating the reflected wave energy value of each area according to the sensing data of the ground penetrating radar and determining a first area with the reflected wave energy value being larger than a first preset threshold value; calculating a reflection coefficient of a reflected wave in the first region and a Fourier spectrum amplitude of the reflected wave according to the sensing data of the first region, and further calculating a cavity response value of the reflected wave in the first region according to the reflection coefficient of the reflected wave in the first region and the Fourier spectrum amplitude of the reflected wave;
and the identification module is used for judging the position corresponding to the first area with the cavity response value larger than a second preset threshold value as an underground cavity.
7. The apparatus for automatically identifying an underground cavity of a road according to claim 6, wherein: the processing module is further configured to calculate a void response value of the reflected wave in the first region according to a weighted summation of a reflection coefficient of the reflected wave in the first region and a fourier spectrum amplitude of the reflected wave.
8. The apparatus according to claim 7, wherein the processing module is further configured to calculate a reflection coefficient of the reflected wave in each of the first areas, and calculate an average value of the reflection coefficients of the reflected waves in each of the first areas; determining an area with a reflection coefficient larger than the average value in each first area as a second area; calculating a hole response value of the reflected wave in the second area according to the weighted summation of the reflection coefficient of the reflected wave in the second area and the Fourier spectrum amplitude of the reflected wave, and judging the position corresponding to the second area with the hole response value larger than a second preset threshold value as an underground hole; wherein the first region comprises a second region.
9. The apparatus for automatically identifying an underground cavity of a road according to any one of claims 6 to 8, further comprising:
the filtering module is used for filtering the sensing data of the acquired ground penetrating radar;
the reflected wave energy value E is calculated as follows:
Figure FDA0002285420810000031
wherein N is the number of sampling points, PiThe amplitude of the reflected wave of the ith sampling point is obtained;
the reflection coefficient R is calculated as follows:
Figure FDA0002285420810000032
wherein, P0The peak value of the direct wave; the fourier spectrum amplitude F (ω) of the reflected wave is calculated as follows:
Figure FDA0002285420810000033
f (t) is a time domain signal of the amplitude of the reflected wave, omega is the amplitude signal frequency of the reflected wave, and t is the amplitude signal sampling time of the reflected wave;
the hole response value K of the reflected wave is calculated in the following way that K is α R + β F (omega), α is a preset reflection coefficient weight, and β is a preset Fourier spectrum amplitude weight.
10. An automatic identification system for underground road cavities, characterized in that it comprises an automatic identification device for underground road cavities according to any one of claims 6 to 9 and a ground penetrating radar, said automatic identification device being connected in communication with the ground penetrating radar.
CN201911163987.8A 2019-11-22 2019-11-22 Automatic identification method, device and system for road underground cavity Active CN111007572B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201911163987.8A CN111007572B (en) 2019-11-22 2019-11-22 Automatic identification method, device and system for road underground cavity

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911163987.8A CN111007572B (en) 2019-11-22 2019-11-22 Automatic identification method, device and system for road underground cavity

Publications (2)

Publication Number Publication Date
CN111007572A true CN111007572A (en) 2020-04-14
CN111007572B CN111007572B (en) 2020-11-10

Family

ID=70113248

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911163987.8A Active CN111007572B (en) 2019-11-22 2019-11-22 Automatic identification method, device and system for road underground cavity

Country Status (1)

Country Link
CN (1) CN111007572B (en)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114076943A (en) * 2021-11-17 2022-02-22 广州市市政工程试验检测有限公司 Rapid automatic comprehensive detection method for urban road underground cavity
CN114137517A (en) * 2022-02-07 2022-03-04 北京中科蓝图科技有限公司 Ground penetrating detection method and device for road and ground penetrating radar device
CN115267143A (en) * 2022-09-28 2022-11-01 江苏筑升土木工程科技有限公司 Road cavity defect detection system and detection method
CN115793086A (en) * 2023-02-07 2023-03-14 武汉新楚光电科技发展有限公司 Optical cable laying environment underground cavity judgment method and system based on optical fiber sensing

Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS6022681A (en) * 1983-07-19 1985-02-05 Japan Radio Co Ltd Method and device for wave height radar observation
JPH0587945A (en) * 1991-09-27 1993-04-09 Doro Hozen Gijutsu Center Cavity inspection method for paved road
JP5719075B1 (en) * 2014-10-06 2015-05-13 ジオ・サーチ株式会社 Cavity thickness exploration method
US20160202346A1 (en) * 2013-06-15 2016-07-14 Howard University Using An MM-Principle to Enforce a Sparsity Constraint on Fast Image Data Estimation From Large Image Data Sets
CN105891896A (en) * 2016-04-25 2016-08-24 湖南科技大学 Feature information recognition and analysis method for underground mined area
CN109001728A (en) * 2018-06-26 2018-12-14 石家庄铁道大学 Method and device based on disease inside Ground Penetrating Radar detection armored concrete
CN109521421A (en) * 2018-01-27 2019-03-26 河南工业大学 A kind of Ground Penetrating Radar thin layer object recognition and detection method
CN109615692A (en) * 2018-12-12 2019-04-12 中国矿业大学(北京) A kind of underground piping surrounding objects interpolation modeling algorithm based on Ground Penetrating Radar
CN110319916A (en) * 2019-06-06 2019-10-11 浙江大学 A kind of limited pond low frequency expanding method based on the measurement of waters interference of reflected wave

Patent Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS6022681A (en) * 1983-07-19 1985-02-05 Japan Radio Co Ltd Method and device for wave height radar observation
JPH0587945A (en) * 1991-09-27 1993-04-09 Doro Hozen Gijutsu Center Cavity inspection method for paved road
US20160202346A1 (en) * 2013-06-15 2016-07-14 Howard University Using An MM-Principle to Enforce a Sparsity Constraint on Fast Image Data Estimation From Large Image Data Sets
JP5719075B1 (en) * 2014-10-06 2015-05-13 ジオ・サーチ株式会社 Cavity thickness exploration method
CN105891896A (en) * 2016-04-25 2016-08-24 湖南科技大学 Feature information recognition and analysis method for underground mined area
CN109521421A (en) * 2018-01-27 2019-03-26 河南工业大学 A kind of Ground Penetrating Radar thin layer object recognition and detection method
CN109001728A (en) * 2018-06-26 2018-12-14 石家庄铁道大学 Method and device based on disease inside Ground Penetrating Radar detection armored concrete
CN109615692A (en) * 2018-12-12 2019-04-12 中国矿业大学(北京) A kind of underground piping surrounding objects interpolation modeling algorithm based on Ground Penetrating Radar
CN110319916A (en) * 2019-06-06 2019-10-11 浙江大学 A kind of limited pond low frequency expanding method based on the measurement of waters interference of reflected wave

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114076943A (en) * 2021-11-17 2022-02-22 广州市市政工程试验检测有限公司 Rapid automatic comprehensive detection method for urban road underground cavity
CN114137517A (en) * 2022-02-07 2022-03-04 北京中科蓝图科技有限公司 Ground penetrating detection method and device for road and ground penetrating radar device
CN115267143A (en) * 2022-09-28 2022-11-01 江苏筑升土木工程科技有限公司 Road cavity defect detection system and detection method
CN115793086A (en) * 2023-02-07 2023-03-14 武汉新楚光电科技发展有限公司 Optical cable laying environment underground cavity judgment method and system based on optical fiber sensing
CN115793086B (en) * 2023-02-07 2023-06-06 武汉新楚光电科技发展有限公司 Optical cable laying environment underground cavity judging method and system based on optical fiber sensing

Also Published As

Publication number Publication date
CN111007572B (en) 2020-11-10

Similar Documents

Publication Publication Date Title
CN111007572B (en) Automatic identification method, device and system for road underground cavity
CN111476088B (en) Asphalt pavement water damage identification model construction method, identification method and system
CN104020495B (en) Automatic underground pipeline parameter recognizing method on basis of ground penetrating radar
CN110376584B (en) Water supply pipeline leakage detection method based on ground penetrating radar image characteristic signal identification
KR101917374B1 (en) Apparatus and method for processing 3d ground penetrating radar signal
EP0943101A1 (en) Method of detecting and classifying objects by means of radar
KR101865016B1 (en) Cavity detection from ground penetrating radar data and cavity detection apparatus
CN112666554A (en) Method for identifying radar amplitude characteristic crack width of asphalt pavement
CN102621531A (en) Rainfall interference suppression method based on X-band radar images
CN108710888A (en) A kind of Coherent Noise in GPR Record method for registering
Wang et al. Automatic asphalt layer interface detection and thickness determination from ground-penetrating radar data
CN111679275A (en) Underground pipeline identification method based on ground penetrating radar
Economou et al. Attenuation analysis of real GPR wavelets: The equivalent amplitude spectrum (EAS)
CN113514833B (en) Sea surface arbitrary point wave direction inversion method based on sea wave image
CN110991507A (en) Road underground cavity identification method, device and system based on classifier
CN111007464B (en) Road underground cavity identification method, device and system based on optimal weighting
CN107656270A (en) The measurement apparatus and measuring method of a kind of contactless buried pipe track forces cun
CN106443674A (en) Ground penetrating radar wave velocity estimation method based on diffraction, imaging and minimum entropy technology
CN102518422A (en) Method for detecting and identifying current stress of oil field downhole casing
CN108303745A (en) A kind of inversion method of the buried cable detection based on electromagnetic wave saturating ground technology
CN111691876B (en) Method, device and storage medium for imaging adjacent well by using acoustic logging
Shen et al. An alternative method for surface current extraction from X-band marine radar images
RU2803396C1 (en) Method for detecting objects and determining their location in real time using distributed optic fibre interferometric vibration sensors
CN112464777B (en) Intelligent estimation method for vertical distance of optical fiber vibration source
JP2551952B2 (en) Invisible object detection method

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
TA01 Transfer of patent application right

Effective date of registration: 20200827

Address after: No.86-a2112, Wanxing Road, Changyang, Fangshan District, Beijing 102400

Applicant after: Beijing Zhongke blueprints Technology Co.,Ltd.

Applicant after: Zhongke yuntu Technology Co., Ltd

Address before: No.86-a2112, Wanxing Road, Changyang, Fangshan District, Beijing 102400

Applicant before: Beijing Zhongke blueprints Technology Co.,Ltd.

TA01 Transfer of patent application right
GR01 Patent grant
GR01 Patent grant