CN109387834B - Dam damage detection method for underground reservoir and electronic equipment - Google Patents
Dam damage detection method for underground reservoir and electronic equipment Download PDFInfo
- Publication number
- CN109387834B CN109387834B CN201810990646.7A CN201810990646A CN109387834B CN 109387834 B CN109387834 B CN 109387834B CN 201810990646 A CN201810990646 A CN 201810990646A CN 109387834 B CN109387834 B CN 109387834B
- Authority
- CN
- China
- Prior art keywords
- ground penetrating
- penetrating radar
- diffracted wave
- diffracted
- angle
- 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
Links
- 238000001514 detection method Methods 0.000 title claims abstract description 20
- 230000000149 penetrating effect Effects 0.000 claims abstract description 145
- 238000003384 imaging method Methods 0.000 claims abstract description 94
- 238000000034 method Methods 0.000 claims abstract description 27
- 238000007781 pre-processing Methods 0.000 claims abstract description 22
- 230000015654 memory Effects 0.000 claims description 16
- 230000005284 excitation Effects 0.000 claims description 10
- 230000006870 function Effects 0.000 claims description 8
- 239000000126 substance Substances 0.000 claims description 3
- 201000010099 disease Diseases 0.000 abstract description 12
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 abstract description 12
- 238000005516 engineering process Methods 0.000 description 7
- 238000012545 processing Methods 0.000 description 4
- 238000011160 research Methods 0.000 description 4
- 230000008859 change Effects 0.000 description 3
- 230000001066 destructive effect Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000001914 filtration Methods 0.000 description 2
- 238000013508 migration Methods 0.000 description 2
- 230000005012 migration Effects 0.000 description 2
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 2
- 238000013473 artificial intelligence Methods 0.000 description 1
- 238000013528 artificial neural network Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 238000010295 mobile communication Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000011435 rock Substances 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
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)
- Geophysics And Detection Of Objects (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
The invention discloses a dam body damage detection method for an underground reservoir and electronic equipment, wherein the method comprises the following steps: acquiring a ground penetrating radar signal of an underground reservoir dam body; preprocessing the acquired ground penetrating radar signal to obtain ground penetrating radar preprocessing data comprising reflected waves and diffracted waves; separating out diffracted waves in the ground penetrating radar preprocessed data; carrying out angle domain compensation on the separated diffracted waves to obtain a primary ground penetrating radar diffracted wave imaging gather; leveling diffracted wave signals in the primary ground penetrating radar diffracted wave imaging channel set to obtain an offset speed model of the ground penetrating radar signals; and performing angle domain compensation on the separated diffracted wave signals again by using an offset velocity model to obtain a final diffracted wave imaging gather of the ground penetrating radar, and superposing along the angle dimension of the final diffracted wave imaging gather to obtain a diffracted wave focusing imaging result. The method improves the identification capability and efficiency of the disease diffraction wave signals of the dam body of the underground reservoir, and timely warns the diseases of the potential safety hazard of the underground reservoir.
Description
Technical Field
The invention relates to the related technical field of an underground reservoir, in particular to a dam body damage detection method of the underground reservoir and electronic equipment.
Background
The ground penetrating radar is an effective device for exploring underground targets developed in recent decades, and has wide application prospects in many fields such as roads, airports, water conservancy, mines, tunnels and the like. In the aspect of artificial dam body disease detection, a large amount of research is also carried out on the ground penetrating radar technology. The research on reservoir dam leakage detection of ground penetrating radar (2005) has been made; in the year 2006, by detecting the characteristic change of the landmine waveform and the continuous change of the same-phase axis of the waveform, the seepage hidden danger of reservoir earth-rock dams in northwest regions is analyzed; the Jiayonmei et al (2012) researches the quality detection of the plastic concrete impervious wall of the reservoir by adopting a ground penetrating radar method; schnewavian et al (2015) have studied the method of detecting a dam cavity by the ground penetrating radar technology. In conclusion, the ground penetrating radar technology is an effective method for detecting disaster bodies such as cavity cracks, and related technology researches relate to mechanism, numerical simulation, data processing algorithm, signal analysis technology, artificial intelligence identification and the like of disease body detection, but only solve the preliminary problems of disease detection such as cavity cracks (Wang Cheng Ming, et al, 2007; Zhang Hai Yan, 2000; Lianping, et al, 2003; Ningyue, 2005; Zhang scholarly, 2000; Zhang Xiaoji, et al, 2008).
At present, the underground space disease detection is mostly based on the diffracted wave of an electromagnetic wave signal as an explanation basis, and the technology can only give out the position of the signal on the radar section where the signal is suddenly changed and the sudden change energy. The method of wavelet analysis and fuzzy neural network can not obtain accurate identification result due to its inherent defects (such as structure is not easy to be determined, and it falls into local minimum point). Therefore, the traditional land radar technology has exploratory and insufficient data processing, and is difficult to effectively solve the problem of detecting the holes and cracks of the dam body of the underground reservoir, and an effective solution is not provided at present for the problem.
Disclosure of Invention
Therefore, it is necessary to provide a method for detecting damage to an underground reservoir dam and an electronic device, aiming at the technical problem that the detection of the cavities and cracks of the underground reservoir dam in the prior art is difficult to effectively solve.
The invention provides a dam body damage detection method for an underground reservoir, which comprises the following steps:
acquiring a ground penetrating radar signal of the dam body of the underground reservoir, wherein the ground penetrating radar signal carries the continuous layered interface inside the dam body and the electrical parameter information of the discontinuous body;
preprocessing the acquired ground penetrating radar signal to obtain ground penetrating radar preprocessing data comprising reflected waves and diffracted waves;
separating out diffracted waves in the ground penetrating radar preprocessed data;
carrying out angle domain compensation on the separated diffracted waves to obtain a primary ground penetrating radar diffracted wave imaging gather;
leveling diffracted wave signals in the primary ground penetrating radar diffracted wave imaging channel set to obtain an offset speed model of the ground penetrating radar signals;
and performing angle domain compensation on the separated diffracted wave signals again by using the offset velocity model to obtain a final diffracted wave imaging gather of the ground penetrating radar, and superposing along the angle dimension of the final diffracted wave imaging gather to obtain a diffracted wave focusing imaging result.
Further, the performing angular domain compensation specifically includes:
the following angle domain compensation formula is adopted for compensation:
wherein M is a diffracted wave multipath focusing result, x is an arbitrary imaging point position in the dam body of the underground reservoir, s is the number of ground penetrating radar excitation sources, alpha is an inclination angle, T is a position of a ground penetrating radar excitation sourceDTo compensate for the coefficients, the form is as follows:
wherein λ represents the ground penetrating radar signal wavelength;
r (x, α) is an angular gather imaging formula of the form:
R(x,α)＝∫∫(θ-α)(t-t0)u(r,t)drdtdθ，
wherein, the parameter is Dirac function, alpha is the angle parameter set by the angle gather, theta, t0The angle parameter and the travel time parameter are respectively calculated according to ray tracing of the speed model, u (r, t) is diffracted wave data in the ground penetrating radar preprocessing data, and r and t respectively correspond to the abscissa and the acquisition time of a receiving channel of the ground penetrating radar.
Further, the leveling up the diffracted wave signals in the primary ground penetrating radar diffracted wave imaging channel set to obtain an offset velocity model of the ground penetrating radar signals specifically includes:
scanning a plurality of different offset speeds, and summing diffracted wave signals of curves corresponding to the different offset speeds to obtain a maximum diffracted wave summation value;
and substituting the offset speed corresponding to the maximum diffracted wave summation value into an angle channel integrated imaging formula R (x, alpha) of the angle domain compensation formula, flattening the diffracted wave signals in the primary ground penetrating radar diffracted wave imaging channel set, and obtaining an offset speed model of the ground penetrating radar signals.
Still further, the curve is calculated as follows:
and substituting the acquired abscissa of the current offset position, the abscissa of the true position of the diffraction point, the initial diffraction point depth, the inclination angle and the preset initial offset speed into an angle domain diffraction wave velocity analysis formula to obtain a relation curve of the diffraction point depth when the diffraction wave offset speed and the inclination angle are alpha.
Still further, the angular domain diffracted wave velocity analysis formula is as follows:
wherein z is0As the initial diffraction point depth, z (α) is the angle of inclinationDepth of diffraction point at α, ρ ═ xm-xd，xmIs the abscissa, x, of the current offset positiondIs the abscissa of the true position of the diffraction point, gamma-vm/v0Deflecting velocity v for diffracted wavesmFrom the initial offset velocity v0The ratio of (a) to (b).
The invention provides an electronic device for detecting damage of an underground reservoir dam body, which comprises:
at least one processor; and the number of the first and second groups,
a memory communicatively coupled to the at least one processor; wherein the content of the first and second substances,
the memory stores instructions executable by the one processor to cause the at least one processor to:
acquiring a ground penetrating radar signal of the dam body of the underground reservoir, wherein the ground penetrating radar signal carries the continuous layered interface inside the dam body and the electrical parameter information of the discontinuous body;
preprocessing the acquired ground penetrating radar signal to obtain ground penetrating radar preprocessing data comprising reflected waves and diffracted waves;
separating out diffracted waves in the ground penetrating radar preprocessed data;
carrying out angle domain compensation on the separated diffracted waves to obtain a primary ground penetrating radar diffracted wave imaging gather;
leveling diffracted wave signals in the primary ground penetrating radar diffracted wave imaging channel set to obtain an offset speed model of the ground penetrating radar signals;
and performing angle domain compensation on the separated diffracted wave signals again by using the offset velocity model to obtain a final diffracted wave imaging gather of the ground penetrating radar, and superposing along the angle dimension of the final diffracted wave imaging gather to obtain a diffracted wave focusing imaging result.
Further, the performing angular domain compensation specifically includes:
the following angle domain compensation formula is adopted for compensation:
wherein M is a diffracted wave multipath focusing result, x is an arbitrary imaging point position in the dam body of the underground reservoir, s is the number of ground penetrating radar excitation sources, alpha is an inclination angle, R (x, alpha) is an angle gather imaging formula, and T isDTo compensate for the coefficients, the form is as follows:
wherein λ represents the ground penetrating radar signal wavelength;
r (x, α) is an angular gather imaging formula of the form:
R(x,α)＝∫∫(θ-α)(t-t0)u(r,t)drdtdθ，
wherein, the parameter is Dirac function, alpha is the angle parameter set by the angle gather, theta, t0The angle parameter and the travel time parameter are respectively calculated according to ray tracing of the speed model, u (r, t) is diffracted wave data in the ground penetrating radar preprocessing data, and r and t respectively correspond to the abscissa and the acquisition time of a receiving channel of the ground penetrating radar.
Further, the leveling up the diffracted wave signals in the primary ground penetrating radar diffracted wave imaging channel set to obtain an offset velocity model of the ground penetrating radar signals specifically includes:
scanning a plurality of different offset speeds, and summing diffracted wave signals of curves corresponding to the different offset speeds to obtain a maximum diffracted wave summation value;
and substituting the offset speed corresponding to the maximum diffracted wave summation value into an angle channel integrated imaging formula R (x, alpha) of the angle domain compensation formula, flattening the diffracted wave signals in the primary ground penetrating radar diffracted wave imaging channel set, and obtaining an offset speed model of the ground penetrating radar signals.
Still further, the curve is calculated as follows:
and substituting the acquired abscissa of the current offset position, the abscissa of the true position of the diffraction point, the initial diffraction point depth, the inclination angle and the preset initial offset speed into an angle domain diffraction wave velocity analysis formula to obtain a relation curve of the diffraction point depth when the diffraction wave offset speed and the inclination angle are alpha.
Still further, the angular domain diffracted wave velocity analysis formula is as follows:
wherein z is0As the initial diffraction point depth, z (α) is the diffraction point depth at an angle of inclination α, ρ ═ xm-xd，xmIs the abscissa, x, of the current offset positiondIs the abscissa of the true position of the diffraction point, gamma-vm/v0Deflecting velocity v for diffracted wavesmFrom the initial offset velocity v0The ratio of (a) to (b).
According to the method, the multipath homing and energy attenuation rules of the diffracted waves are considered, and meanwhile, the focusing performance of the diffracted waves is greatly enhanced by obtaining an accurate offset speed model, so that the identification capability and efficiency of disease diffracted wave signals of the dam body of the underground reservoir are improved, and the diseases of potential safety hazards of the underground reservoir are early warned in time.
Drawings
FIG. 1 is a work flow chart of a dam damage detection method of an underground reservoir according to the invention;
FIG. 2 is a flow chart of the work of the method for detecting damage to an underground reservoir dam body according to the preferred embodiment of the invention;
fig. 3 is a schematic diagram of a hardware structure of an electronic device for detecting damage to an underground reservoir dam according to the present invention.
Detailed Description
The invention is described in further detail below with reference to the figures and specific examples.
Fig. 1 shows a work flow chart of the method for detecting damage to an underground reservoir dam body, which includes:
s101, acquiring a ground penetrating radar signal of the dam body of the underground reservoir, wherein the ground penetrating radar signal carries a continuous layered interface inside the dam body and electrical parameter information of a discontinuous body;
step S102, preprocessing the acquired ground penetrating radar signal to obtain ground penetrating radar preprocessing data comprising reflected waves and diffracted waves;
s103, separating diffracted waves in the ground penetrating radar preprocessed data;
step S104, performing angle domain compensation on the separated diffracted waves to obtain a primary ground penetrating radar diffracted wave imaging gather;
s105, leveling diffracted wave signals in the primary ground penetrating radar diffracted wave imaging channel set to obtain an offset speed model of the ground penetrating radar signals;
and S106, performing angle domain compensation on the separated diffracted wave signals again by using the offset velocity model to obtain a final diffracted wave imaging gather of the ground penetrating radar, and superposing the final diffracted wave imaging gather along the angle dimension to obtain a diffracted wave focusing imaging result.
Specifically, the method for detecting damage to the dam body of the underground reservoir comprises the following steps: acquiring a ground penetrating radar signal of a dam body of the underground reservoir, wherein the ground penetrating radar signal carries electrical parameter information of a continuous layered interface, a cavity, a crack and other discontinuous bodies in the dam body; preprocessing the acquired ground penetrating radar signal, including removing pulse interference, time gain and the like, so that the preprocessed ground penetrating radar signal mainly comprises reflected waves and diffracted waves; separating out diffracted waves in the ground penetrating radar preprocessed data by using a plane wave destructive filtering method, wherein the diffracted waves correspond to discontinuous bodies such as cavities, cracks and the like in the dam body of the underground reservoir; aiming at the separated diffracted waves, obtaining a ground penetrating radar diffracted wave imaging gather by an angle domain compensation type imaging method; by utilizing an angle domain diffracted wave velocity analysis formula, leveling diffracted wave signals concentrated by a diffracted wave imaging channel of the ground penetrating radar to obtain an accurate ground penetrating radar signal offset velocity model; and obtaining a final ground penetrating radar diffracted wave imaging channel set by using the separated diffracted waves and the obtained offset velocity model according to an angle domain compensation type imaging method, and superposing the final ground penetrating radar diffracted wave imaging channel set along the angle dimension of the diffracted wave imaging channel set to obtain a diffracted wave focusing imaging result.
Compared with the prior art, the method and the device have the technical problems that the detection precision is poor due to the fact that multipath propagation and migration velocity modeling of the diffracted waves are not considered when the disease target body is detected, the method and the device can more accurately detect the interior of the dam body of the underground reservoir through diffracted wave separation, angular trace speed analysis and multipath migration imaging, the technical problem that the detection resolution ratio is limited due to the fact that a single data processing mode is used when small-scale holes and cracks are detected in the prior art is solved, and therefore the technical effect of detecting potential small-scale disease bodies with high resolution is achieved.
According to the method, the multipath homing and energy attenuation rules of the diffracted waves are considered, and meanwhile, the focusing performance of the diffracted waves is greatly enhanced by obtaining an accurate offset speed model, so that the identification capability and efficiency of disease diffracted wave signals of the dam body of the underground reservoir are improved, and the diseases of potential safety hazards of the underground reservoir are early warned in time.
In one embodiment, the performing angular domain compensation specifically includes:
the following angle domain compensation formula is adopted for compensation:
wherein M is a diffracted wave multipath focusing result, x is an arbitrary imaging point position in the dam body of the underground reservoir, s is the number of ground penetrating radar excitation sources, alpha is an inclination angle, R (x, alpha) is an angle gather imaging formula, and T isDTo compensate for the coefficients, the form is as follows:
wherein λ represents the ground penetrating radar signal wavelength;
r (x, α) is an angular gather imaging formula of the form:
R(x,α)＝∫∫(θ-α)(t-t0)u(r,t)drdtdθ，
wherein, the parameter is Dirac function, alpha is the angle parameter set by the angle gather, theta, t0Respectively according to the velocity modeAnd (c) calculating an angle parameter and a travel time parameter by ray tracing, wherein u (r, t) is diffracted wave data in the ground penetrating radar preprocessing data, and r and t respectively correspond to the abscissa and the acquisition time of a receiving channel of the ground penetrating radar.
The angle domain compensation formula of the present embodiment is used in step S104 and step S106, and R (x, α) is an existing angle gather imaging formula. In the embodiment, the diffracted wave imaging angle gather is obtained by calculating a traditional angle gather imaging formula and compensating for diffracted wave energy attenuation loss.
In one embodiment, the leveling up the diffracted wave signals in the primary ground penetrating radar diffracted wave imaging channel set to obtain an offset velocity model of the ground penetrating radar signal specifically includes:
scanning a plurality of different offset speeds, and summing diffracted wave signals of curves corresponding to the different offset speeds to obtain a maximum diffracted wave summation value;
and substituting the offset speed corresponding to the maximum diffracted wave summation value into an angle channel integrated imaging formula R (x, alpha) of the angle domain compensation formula, flattening the diffracted wave signals in the primary ground penetrating radar diffracted wave imaging channel set, and obtaining an offset speed model of the ground penetrating radar signals.
In one embodiment, the curve is calculated as follows:
and substituting the acquired abscissa of the current offset position, the abscissa of the true position of the diffraction point, the initial diffraction point depth, the inclination angle and the preset initial offset speed into an angle domain diffraction wave velocity analysis formula to obtain a relation curve of the diffraction point depth when the diffraction wave offset speed and the inclination angle are alpha.
In one embodiment, the angular domain diffracted wave velocity analysis formula is:
wherein z is0As the initial diffraction point depth, z (α) is the diffraction point depth at an angle of inclination α, ρ ═ xm-xd，xmAs the current offset positionAbscissa of (a), xdIs the abscissa of the true position of the diffraction point, gamma-vm/v0Deflecting velocity v for diffracted wavesmFrom the initial offset velocity v0The ratio of (a) to (b).
According to the method, the accurate ground penetrating radar signal speed model is obtained by leveling the diffracted wave signals concentrated in the imaging channel of the ground penetrating radar diffracted wave according to an angle domain diffracted wave speed analysis formula.
Fig. 2 is a flowchart illustrating a method for detecting damage to an underground reservoir dam according to a preferred embodiment of the present invention, including:
step S201, acquiring a ground penetrating radar signal of an underground reservoir dam body;
the ground penetrating radar signal carries electrical parameter information of a continuous layered interface and discontinuous bodies such as cavities and cracks in the dam body;
step S202, preprocessing the acquired ground penetrating radar signal to obtain a preprocessed ground penetrating radar signal which mainly comprises a reflected wave and a diffracted wave;
step S203, separating diffracted waves in the ground penetrating radar preprocessed data by using a plane wave destructive filtering method;
step S204, aiming at the separated diffracted waves, obtaining a ground penetrating radar diffracted wave imaging gather by an angle domain compensation type imaging method;
specifically, the following angle domain compensation formula is adopted for compensation:
wherein M is a diffracted wave multipath focusing result, x is an arbitrary imaging point position in the dam body of the underground reservoir, s is the number of ground penetrating radar excitation sources, alpha is an inclination angle, R (x, alpha) is an angle gather imaging formula, and T isDTo compensate for the coefficients, the form is as follows:
wherein λ represents the ground penetrating radar signal wavelength;
step S205, by utilizing an angle domain diffracted wave velocity analysis formula, leveling diffracted wave signals concentrated by a diffracted wave imaging channel of the ground penetrating radar to obtain an accurate ground penetrating radar signal offset velocity model;
specifically, scanning a plurality of different offset speeds, and summing diffracted wave signals of curves corresponding to the different offset speeds to obtain a maximum diffracted wave summation value;
substituting the offset speed corresponding to the maximum diffracted wave summation value into an angle channel integrated imaging formula R (x, alpha) of the angle domain compensation formula, flattening diffracted wave signals in the primary ground penetrating radar diffracted wave imaging channel set, and obtaining an offset speed model of the ground penetrating radar signals, wherein the curve is calculated in the following mode:
substituting the obtained abscissa of the current offset position, the abscissa of the true position of the diffraction point, the initial diffraction point depth, the inclination angle and the preset initial offset speed into an angle domain diffraction wave velocity analysis formula to obtain a relation curve of the diffraction wave offset speed and the diffraction point depth when the inclination angle is alpha, wherein the angle domain diffraction wave velocity analysis formula is as follows:
wherein z is0As the initial diffraction point depth, z (α) is the diffraction point depth at an angle of inclination α, ρ ═ xm-xd，xmIs the abscissa, x, of the current offset positiondIs the abscissa of the true position of the diffraction point, gamma-vm/v0Deflecting velocity v for diffracted wavesmFrom the initial offset velocity v0The ratio of (a) to (b).
Step S206, obtaining a final ground penetrating radar diffracted wave imaging gather by utilizing the separated diffracted waves and the obtained offset speed model and according to an angle domain compensation type imaging method, namely substituting the offset speed corresponding to the maximum diffracted wave summation value into an angle domain compensation formula of an angle gather imaging formula R (x, alpha), and obtaining the final ground penetrating radar diffracted wave imaging gather through the separated diffracted waves;
and step S207, superposing along the angle dimension of the diffraction wave imaging channel set to obtain a diffraction wave focusing imaging result.
Fig. 3 is a schematic diagram of a hardware structure of an electronic device for detecting damage to an underground reservoir dam according to the present invention, which includes:
at least one processor 301; and the number of the first and second groups,
a memory 302 communicatively coupled to the at least one processor 301; wherein the content of the first and second substances,
the memory 302 stores instructions executable by the one processor to cause the at least one processor to:
acquiring a ground penetrating radar signal of the dam body of the underground reservoir, wherein the ground penetrating radar signal carries the continuous layered interface inside the dam body and the electrical parameter information of the discontinuous body;
preprocessing the acquired ground penetrating radar signal to obtain ground penetrating radar preprocessing data comprising reflected waves and diffracted waves;
separating out diffracted waves in the ground penetrating radar preprocessed data;
carrying out angle domain compensation on the separated diffracted waves to obtain a primary ground penetrating radar diffracted wave imaging gather;
leveling diffracted wave signals in the primary ground penetrating radar diffracted wave imaging channel set to obtain an offset speed model of the ground penetrating radar signals;
and performing angle domain compensation on the separated diffracted wave signals again by using the offset velocity model to obtain a final diffracted wave imaging gather of the ground penetrating radar, and superposing along the angle dimension of the final diffracted wave imaging gather to obtain a diffracted wave focusing imaging result.
In fig. 3, a processor 302 is taken as an example.
The electronic device may further include: an input device 303 and an output device 304.
The processor 301, the memory 302, the input device 303 and the display device 304 may be connected by a bus or other means, and are illustrated as being connected by a bus.
The memory 302 is a non-volatile computer-readable storage medium, and can be used for storing non-volatile software programs, non-volatile computer-executable programs, and modules, such as program instructions/modules corresponding to the method for detecting damage to an underground reservoir dam in the embodiment of the present application, for example, the method flows shown in fig. 1 and fig. 2. The processor 301 executes various functional applications and data processing by running nonvolatile software programs, instructions and modules stored in the memory 302, so as to implement the dam damage detection method of the underground water reservoir in the above embodiment.
The memory 302 may include a storage program area and a storage data area, wherein the storage program area may store an operating system, an application program required for at least one function; the storage data area may store data created according to use of the dam damage detection method of the underground reservoir, and the like. Further, the memory 302 may include high speed random access memory, and may also include non-volatile memory, such as at least one magnetic disk storage device, flash memory device, or other non-volatile solid state storage device. In some embodiments, the memory 302 may optionally include a memory remotely located from the processor 301, and these remote memories may be connected via a network to a device that performs the method for dam damage detection in an underground reservoir. Examples of such networks include, but are not limited to, the internet, intranets, local area networks, mobile communication networks, and combinations thereof.
The input device 303 may receive input of user clicks and generate signal inputs related to user settings and function control of the underground reservoir dam damage detection method. The display device 304 may include a display screen or the like.
The method for detecting dam damage in an underground reservoir according to any of the above method embodiments is performed when the one or more modules are stored in the memory 302 and executed by the one or more processors 301.
In one embodiment, the performing angular domain compensation specifically includes:
the following angle domain compensation formula is adopted for compensation:
wherein M is a diffracted wave multipath focusing result, x is an arbitrary imaging point position in the dam body of the underground reservoir, s is the number of ground penetrating radar excitation sources, alpha is an inclination angle, R (x, alpha) is an angle gather imaging formula, and T isDTo compensate for the coefficients, the form is as follows:
wherein λ represents the ground penetrating radar signal wavelength;
r (x, α) is an angular gather imaging formula of the form:
R(x,α)＝∫∫(θ-α)(t-t0)u(r,t)drdtdθ，
wherein, the parameter is Dirac function, alpha is the angle parameter set by the angle gather, theta, t0The angle parameter and the travel time parameter are respectively calculated according to ray tracing of the speed model, u (r, t) is diffracted wave data in the ground penetrating radar preprocessing data, and r and t respectively correspond to the abscissa and the acquisition time of a receiving channel of the ground penetrating radar.
In one embodiment, the leveling up the diffracted wave signals in the primary ground penetrating radar diffracted wave imaging channel set to obtain an offset velocity model of the ground penetrating radar signal specifically includes:
scanning a plurality of different offset speeds, and summing diffracted wave signals of curves corresponding to the different offset speeds to obtain a maximum diffracted wave summation value;
and substituting the offset speed corresponding to the maximum diffracted wave summation value into an angle channel integrated imaging formula R (x, alpha) of the angle domain compensation formula, flattening the diffracted wave signals in the primary ground penetrating radar diffracted wave imaging channel set, and obtaining an offset speed model of the ground penetrating radar signals.
In one embodiment, the curve is calculated as follows:
and substituting the acquired abscissa of the current offset position, the abscissa of the true position of the diffraction point, the initial diffraction point depth, the inclination angle and the preset initial offset speed into an angle domain diffraction wave velocity analysis formula to obtain a relation curve of the diffraction point depth when the diffraction wave offset speed and the inclination angle are alpha.
In one embodiment, the angular domain diffracted wave velocity analysis formula is:
wherein z is0As the initial diffraction point depth, z (α) is the diffraction point depth at an angle of inclination α, ρ ═ xm-xd，xmIs the abscissa, x, of the current offset positiondIs the abscissa of the true position of the diffraction point, gamma-vm/v0Deflecting velocity v for diffracted wavesmFrom the initial offset velocity v0The ratio of (a) to (b).
The above-mentioned embodiments only express several embodiments of the present invention, and the description thereof is more specific and detailed, but not construed as limiting the scope of the present invention. It should be noted that, for a person skilled in the art, several variations and modifications can be made without departing from the inventive concept, which falls within the scope of the present invention. Therefore, the protection scope of the present patent shall be subject to the appended claims.
Claims (6)
1. A dam body damage detection method for an underground reservoir is characterized by comprising the following steps:
acquiring a ground penetrating radar signal of the dam body of the underground reservoir, wherein the ground penetrating radar signal carries the continuous layered interface inside the dam body and the electrical parameter information of the discontinuous body;
preprocessing the acquired ground penetrating radar signal to obtain ground penetrating radar preprocessing data comprising reflected waves and diffracted waves;
separating out diffracted waves in the ground penetrating radar preprocessed data;
carrying out angle domain compensation on the separated diffracted waves to obtain a primary ground penetrating radar diffracted wave imaging gather;
leveling diffracted wave signals in the primary ground penetrating radar diffracted wave imaging channel set to obtain an offset speed model of the ground penetrating radar signals;
performing angle domain compensation on the separated diffracted wave signals again by using the offset velocity model to obtain a final diffracted wave imaging gather of the ground penetrating radar, and stacking along the angle dimension of the final diffracted wave imaging gather to obtain a diffracted wave focusing imaging result;
the performing angle domain compensation specifically includes:
the following angle domain compensation formula is adopted for compensation:
wherein M is a diffracted wave multipath focusing result, x is an arbitrary imaging point position in the dam body of the underground reservoir, s is the number of ground penetrating radar excitation sources, alpha is an inclination angle, T is a position of a ground penetrating radar excitation sourceDTo compensate for the coefficients, the form is as follows:
wherein λ represents the ground penetrating radar signal wavelength;
r (x, α) is an angular gather imaging formula of the form:
R(x,α)＝∫∫(θ-α)(t-t0)u(r,t)drdtdθ，
wherein is the Dirac function, θ, t0The angle parameter and the travel time parameter are respectively calculated by ray tracing according to a speed model, u (r, t) is diffracted wave data in the ground penetrating radar preprocessing data, and r and t respectively correspond to the abscissa and the acquisition time of a receiving channel of the ground penetrating radar;
flattening diffracted wave signals concentrated in the primary ground penetrating radar diffracted wave imaging channel to obtain an offset speed model of the ground penetrating radar signals, and the method specifically comprises the following steps:
scanning a plurality of different offset speeds, and summing diffracted wave signals of curves corresponding to the different offset speeds to obtain a maximum diffracted wave summation value;
and substituting the offset speed corresponding to the maximum diffracted wave summation value into an angle channel integrated imaging formula R (x, alpha) of the angle domain compensation formula, flattening the diffracted wave signals in the primary ground penetrating radar diffracted wave imaging channel set, and obtaining an offset speed model of the ground penetrating radar signals.
2. The method for detecting damage to an underground reservoir dam according to claim 1, wherein the curve is calculated as follows:
and substituting the acquired abscissa of the current offset position, the abscissa of the true position of the diffraction point, the initial diffraction point depth, the inclination angle and the preset initial offset speed into an angle domain diffraction wave velocity analysis formula to obtain a relation curve of the diffraction point depth when the diffraction wave offset speed and the inclination angle are alpha.
3. The method for detecting damage to an underground reservoir dam according to claim 2, wherein the angular domain diffraction wave velocity analysis formula is:
wherein z is0As the initial diffraction point depth, z (α) is the diffraction point depth at an angle of inclination α, ρ ═ xm-xd，xmIs the abscissa, x, of the current offset positiondIs the abscissa of the true position of the diffraction point, gamma-vm/v0Deflecting velocity v for diffracted wavesmFrom the initial offset velocity v0The ratio of (a) to (b).
4. An electronic device for detecting damage of an underground reservoir dam body is characterized by comprising:
at least one processor; and the number of the first and second groups,
a memory communicatively coupled to the at least one processor; wherein the content of the first and second substances,
the memory stores instructions executable by the one processor to cause the at least one processor to:
acquiring a ground penetrating radar signal of the dam body of the underground reservoir, wherein the ground penetrating radar signal carries the continuous layered interface inside the dam body and the electrical parameter information of the discontinuous body;
preprocessing the acquired ground penetrating radar signal to obtain ground penetrating radar preprocessing data comprising reflected waves and diffracted waves;
separating out diffracted waves in the ground penetrating radar preprocessed data;
carrying out angle domain compensation on the separated diffracted waves to obtain a primary ground penetrating radar diffracted wave imaging gather;
leveling diffracted wave signals in the primary ground penetrating radar diffracted wave imaging channel set to obtain an offset speed model of the ground penetrating radar signals;
performing angle domain compensation on the separated diffracted wave signals again by using the offset velocity model to obtain a final diffracted wave imaging gather of the ground penetrating radar, and stacking along the angle dimension of the final diffracted wave imaging gather to obtain a diffracted wave focusing imaging result;
the performing angle domain compensation specifically includes:
the following angle domain compensation formula is adopted for compensation:
wherein M is a diffracted wave multipath focusing result, x is an arbitrary imaging point position in the dam body of the underground reservoir, s is the number of ground penetrating radar excitation sources, alpha is an inclination angle, T is a position of a ground penetrating radar excitation sourceDTo compensate for the coefficients, the form is as follows:
wherein λ represents the ground penetrating radar signal wavelength;
r (x, α) is an angular gather imaging formula of the form:
R(x,α)＝∫∫(θ-α)(t-t0)u(r,t)drdtdθ，
wherein, the parameter is Dirac function, alpha is the angle parameter set by the angle gather, theta, t0The angle parameter and the travel time parameter are respectively calculated by ray tracing according to a speed model, u (r, t) is diffracted wave data in the ground penetrating radar preprocessing data, and r and t respectively correspond to the abscissa and the acquisition time of a receiving channel of the ground penetrating radar;
flattening diffracted wave signals concentrated in the primary ground penetrating radar diffracted wave imaging channel to obtain an offset speed model of the ground penetrating radar signals, and the method specifically comprises the following steps:
scanning a plurality of different offset speeds, and summing diffracted wave signals of curves corresponding to the different offset speeds to obtain a maximum diffracted wave summation value;
and substituting the offset speed corresponding to the maximum diffracted wave summation value into an angle channel integrated imaging formula R (x, alpha) of the angle domain compensation formula, flattening the diffracted wave signals in the primary ground penetrating radar diffracted wave imaging channel set, and obtaining an offset speed model of the ground penetrating radar signals.
5. The electronic device of claim 4, wherein the curve is calculated as follows:
and substituting the acquired abscissa of the current offset position, the abscissa of the true position of the diffraction point, the initial diffraction point depth, the inclination angle and the preset initial offset speed into an angle domain diffraction wave velocity analysis formula to obtain a relation curve of the diffraction point depth when the diffraction wave offset speed and the inclination angle are alpha.
6. The electronic device of claim 5, wherein the angular domain velocity analysis formula is:
wherein z is0As the initial diffraction point depth, z (α) is the angle of inclinationDepth of diffraction point at α, ρ ═ xm-xd，xmIs the abscissa, x, of the current offset positiondIs the abscissa of the true position of the diffraction point, gamma-vm/v0Deflecting velocity v for diffracted wavesmFrom the initial offset velocity v0The ratio of (a) to (b).
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810990646.7A CN109387834B (en) | 2018-08-28 | 2018-08-28 | Dam damage detection method for underground reservoir and electronic equipment |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810990646.7A CN109387834B (en) | 2018-08-28 | 2018-08-28 | Dam damage detection method for underground reservoir and electronic equipment |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109387834A CN109387834A (en) | 2019-02-26 |
| CN109387834B true CN109387834B (en) | 2020-10-27 |
Family
ID=65417591
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810990646.7A Active CN109387834B (en) | 2018-08-28 | 2018-08-28 | Dam damage detection method for underground reservoir and electronic equipment |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109387834B (en) |
Families Citing this family (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110632667B (en) * | 2019-10-18 | 2020-12-18 | 徐州工程学院 | Hidden collapse column advanced detection method based on shock wave shock condition |
| CN111025383B (en) * | 2019-11-21 | 2021-09-24 | 徐州工程学院 | Method for qualitatively judging water filling condition of tunnel front karst cave based on diffracted transverse waves |
| CN111123264A (en) * | 2019-12-31 | 2020-05-08 | 桂林理工大学 | Ground penetrating radar reverse time migration imaging method considering medium frequency dispersion and attenuation compensation |
| CN112924962B (en) * | 2021-01-29 | 2023-02-21 | 上海匀羿电磁科技有限公司 | Underground pipeline lateral deviation filtering detection and positioning method |
| CN113945987B (en) * | 2021-10-15 | 2022-04-26 | 中国矿业大学(北京) | Method and device for detecting disease geologic body and electronic equipment |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106199704A (en) * | 2016-09-13 | 2016-12-07 | 中国海洋石油总公司 | A kind of Three-dimendimal fusion submarine cable seismic data velocity modeling method |
| CN107817523A (en) * | 2017-10-30 | 2018-03-20 | 中国矿业大学(北京) | The analysis method and device of diffracted wave migration velocity |
| CN108226999A (en) * | 2018-01-19 | 2018-06-29 | 中国石油化工股份有限公司 | The processing method of the small scale fracture hole body information of carbonate rock |
-
2018
- 2018-08-28 CN CN201810990646.7A patent/CN109387834B/en active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106199704A (en) * | 2016-09-13 | 2016-12-07 | 中国海洋石油总公司 | A kind of Three-dimendimal fusion submarine cable seismic data velocity modeling method |
| CN107817523A (en) * | 2017-10-30 | 2018-03-20 | 中国矿业大学(北京) | The analysis method and device of diffracted wave migration velocity |
| CN108226999A (en) * | 2018-01-19 | 2018-06-29 | 中国石油化工股份有限公司 | The processing method of the small scale fracture hole body information of carbonate rock |
Non-Patent Citations (2)
| Title |
|---|
| VTI介质快速偏移速度分析;于彩霞 等;《石油地球物理勘探》;20131031;第48卷(第5期);全文 * |
| 倾角成像道集中反射波和绕射波特征分析及;徐辉 等;《石油物探》;20150331;第54卷(第2期);全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109387834A (en) | 2019-02-26 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109387834B (en) | Dam damage detection method for underground reservoir and electronic equipment | |
| Paik et al. | Image processing-based mine detection techniques: A review | |
| CN104569972B (en) | Plant root system three-dimensional configuration nondestructive testing method | |
| Andersen et al. | Processing and performance of topobathymetric lidar data for geomorphometric and morphological classification in a high-energy tidal environment | |
| CN103941257A (en) | Navigation radar image sea surface wind direction inversion method based on wave number energy spectrum | |
| CN102201125A (en) | Method for visualizing three-dimensional imaging sonar data | |
| Wiehle et al. | Automated waterline detection in the Wadden Sea using high-resolution TerraSAR-X images | |
| KR101426226B1 (en) | Signal processing method of radar | |
| Wang et al. | Qualitative research: The impact of root orientation on coarse roots detection using ground-penetrating radar (GPR) | |
| CN105654477A (en) | Detecting and positioning method for strip-shaped underground object | |
| CN104359422A (en) | Device and method for detecting geometric outline of cavity through borehole camera shooting | |
| Li et al. | Artifact suppression of back-projection algorithm under multiple buried objects | |
| Chen et al. | Feasibility study of using mobile laser scanning point cloud data for GNSS line of sight analysis | |
| Wang et al. | Submarine karst morphology detection method based on multi-frequency ultrasound | |
| CN110221289A (en) | Object detection method for three coordinate Connectors for Active Phased Array Radar | |
| Leng et al. | Satellite derived active-passive fusion bathymetry based on GRU model | |
| DE112016006775B4 (en) | Radar signal processing apparatus and radar signal processing method | |
| CN112666552B (en) | Ground penetrating radar data background clutter self-adaptive removing method | |
| CN113468467A (en) | Data detection and preprocessing method for lightning field synchronous ground-air coupling inversion | |
| CN113945925A (en) | Method and apparatus for mitigating blurring effects in ground-based radar images | |
| Saepuloh et al. | Characterizing surface manifestation of geothermal system under Torrid Zone using Synthetic Aperture Radar (SAR) data | |
| CN114488346B (en) | Underground space abnormal body detection method, device, equipment and medium | |
| Xie et al. | Study on radargram characteristics of the backfill grouting quality evaluation of a shield tunnel using GPR | |
| Hao et al. | Numerical simulation method of ground-penetrating radar identification of piping based on finite-difference-time-domain method | |
| Pambudi et al. | Robust detection for forward-looking GPR in rough-surface clutter environments |
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 |