CN109541586A - A kind of Ground Penetrating Radar quick detection road engineering concrete cushion thickness analysis method - Google Patents
A kind of Ground Penetrating Radar quick detection road engineering concrete cushion thickness analysis method Download PDFInfo
- Publication number
- CN109541586A CN109541586A CN201811326015.1A CN201811326015A CN109541586A CN 109541586 A CN109541586 A CN 109541586A CN 201811326015 A CN201811326015 A CN 201811326015A CN 109541586 A CN109541586 A CN 109541586A
- Authority
- CN
- China
- Prior art keywords
- module
- ground penetrating
- data
- penetrating radar
- parameter
- 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.)
- Pending
Links
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
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B7/00—Measuring arrangements characterised by the use of electric or magnetic techniques
- G01B7/02—Measuring arrangements characterised by the use of electric or magnetic techniques for measuring length, width or thickness
- G01B7/06—Measuring arrangements characterised by the use of electric or magnetic techniques for measuring length, width or thickness for measuring thickness
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Electromagnetism (AREA)
- General Physics & Mathematics (AREA)
- Computer Networks & Wireless Communication (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
The invention belongs to Ground Penetrating Radar field of measuring technique, a kind of Ground Penetrating Radar quick detection road engineering concrete cushion Thickness Analysis system and analysis method are disclosed, Ground Penetrating Radar quick detection road engineering concrete cushion Thickness Analysis system includes: parameter setting module, data acquisition module, main control module, data de-noising module, thickness measure module, modeling module, data memory module, display module.The present invention is up to the purpose for automatically removing coherent noise by data de-noising module, avoids the trouble for needing artificially to judge singular value based on SVD decomposition method, singular value is directly judged according to root mean square high parameter, not only accurate but also time saving;Simultaneously by modeling module can accurate simulation go out communication process of the radar signal in layered medium, the scope of application covering far field of the model and near field, to provide strong theoretical model using characteristic progress accurate detection of the Ground Penetrating Radar to layered medium.
Description
Technical field
The invention belongs to Ground Penetrating Radar field of measuring technique more particularly to a kind of Ground Penetrating Radar quick detection road engineering are mixed
Solidifying soil padding thickness analysis method.
Background technique
Road engineering refers to the overall process of the planning, design, construction, maintenance and the management work that carry out using road as object
And its engineering entity being engaged in.It is the same with the civil engineering of other any classes, road distance have apparent technology, economy and
The characteristic of management etc..Roadbed is both the main body of route and the basis on road surface, and bears vehicular load jointly with road surface.
Roadbed is divided into embankment formula, road cutting type protecting, cut-fill three classes by the digging situation of filling out of its section.Road shoulder is road surface two sides shoulder
Within area, be used to supporting road surface, for interim parked vehicles or pedestrian's walking.Roadbed earthwork is by the difficulty or ease excavated
It is divided into earth work (loosening the soil, common soil, pan soil three-level) and stonework (soft rock, secondary hard stone, hard stone three-level).However, existing
Ground Penetrating Radar is easy the accuracy of measurement data affected by noise;Existing Ground Penetrating Radar modeling has ignored antenna and believes radar simultaneously
Number influence;The use scope of model is limited, and mainly detected target (layered medium, such as highway) must be positioned at spy land mine
Up to the far field of system, i.e., farther out from radar, the divergence loss of radar signal is increased in this way, reduces signal-to-noise ratio, thus seriously
Affect measurement accuracy.
In conclusion problem of the existing technology is:
Existing Ground Penetrating Radar is easy the accuracy of measurement data affected by noise;Existing Ground Penetrating Radar modeling has ignored simultaneously
Influence of the antenna to radar signal;The use scope of model is limited, and mainly detected target (layered medium, such as highway) must
That is, farther out from radar the far field that ground penetrating radar system must be located at increases the divergence loss of radar signal in this way, thus serious shadow
Measurement accuracy is rung.
The Numerical accuracy that the velocity of electromagnetic wave data of the concrete of the prior art obtain is poor.
It is low that the prior art obtains concrete cushion exact thickness h data accuracy.
Summary of the invention
In view of the problems of the existing technology, the present invention provides a kind of Ground Penetrating Radar quick detection road engineering concrete
Cushion thickness analysis method.
The invention is realized in this way a kind of Ground Penetrating Radar quick detection road engineering concrete cushion Thickness Analysis side
Method, comprising:
Road engineering radar image in Project Areas is detected using profile method by data acquisition module, passes through site road coring
Situation is modified velocity of electromagnetic wave, and the velocity of electromagnetic wave for obtaining the concrete under specific proportion is v2 data;It specifically includes:
1) signal model of the Project Areas road engineering radar image time-frequency overlapping MASK detected are as follows:
Wherein, N is the signal component number of time-frequency overlapped signal, and n (t) is additive white Gaussian noise, siIt (t) is time-frequency weight
The signal component of folded signal, is expressed asA in formulaiIndicate signal point
The amplitude of amount, ai(m) indicate that the symbol of signal component, p (t) indicate molding filter function, TiIndicate the symbol of signal component
Period, fciIndicate the carrier frequency of signal component,Indicate the phase of signal component;
It is rightNonlinear transformation is carried out,
It obtains: y (t) are as follows:
Y (t)=x (t)+n (t);
Wherein, x (t) is digital modulation signals, and n (t) is the impulsive noise of obedience standard S α S distribution, the parsing shape of x (t)
Formula indicates are as follows:
Wherein, N is sampling number, anFor the information symbol of transmission, in MASK signal, an=0,1,2 ..., M-1, M are
Order of modulation, an=ej2πε/M, ε=0,1,2 ..., M-1, g (t) expression rectangle molding pulse, TbIndicate symbol period, fcIt indicates
Carrier frequency, carrier wave initial phaseIt is the equally distributed random number in [0,2 π];
2) the multipath space of n signal is constructed are as follows:
Wherein,Q is sampling number, and K is maximum delay, by
Maximum detectable range Rmax/ c is obtained, wherein xreciIt (t) is reference signal, RmaxFor maximum detectable range, c is the light velocity;
3) then inhibit direct wave and its multipath using principle of least square method, min will be sought | | Ssur-Xref·α||2Conversion
To askIt obtains:
Substitute into αestim, it solves:
Wherein, SsurFor echo channel signal, α is adaptive weight, αestimFor the estimated value of α,For XrefTransposition,
SotherFor remaining echo and noise final in echo channel;
Main control module carries out denoising to Coherent Noise in GPR Record by data de-noising module;Existed by thickness measure module
Surface layer and the bed course for distinguishing Project Areas road in EKKO-Project software by the processing of gain declining profit, measure in velocity of electromagnetic wave
Road evenness and concrete cushion thickness h under v2;
In the DATA REASONING of concrete cushion thickness h, thickness h error parameter need to be handled, be specifically included:
Using inverse Fourier transform, the scattering parameter that frequency domain is measured transforms to time domain, obtains time domain impulse response, reflection
Parameter corresponds to Time Domain Reflectometry response, and configured transmission corresponds to time-domain transmission response;
According to the position of the first two pulse on a timeline in reflex response and transmission response, four time domain choosings are constructed respectively
Pass function;
Using gate function, the reflex response and transmission response of time domain are gated, extract respectively reflex response and
The first two pulse in transmission response;
Time domain impulse after time domain gating is passed through into Fourier transformation respectively, obtains frequency domain gated data;
Information comprising detected concrete bed course in EKKO-Project software in frequency domain gated data, utilizes obtained letter
Breath, according to formula construction compensation factor Fcf(i);Using following formula, compensation factor F is constructedcf(i):
(i=1,2 ...)
G1(i)~G4It (i) is the frequency domain gated data obtained in upper step;
R (i) is ratio factor;
Fcf(i) compensation factor;
Using occlusion compensation formula, the reflection parameters F of error is not coveredS11(i) and configured transmission FS21(i);
The model of ground penetrating radar system is constructed by modeling module;Ground Penetrating Radar measurement is stored by data memory module
Data information;
Further, the Ground Penetrating Radar quick detection road engineering concrete cushion thickness analysis method further include:
By parameter setting module, radar parameter: radar frequency 250MHz is set;When window 36ns;Velocity of electromagnetic wave v1=
0.1m/ns;
The thickness data information of Ground Penetrating Radar measurement is shown by display module.
Further, using following occlusion compensation formula, the reflection parameters F of error is not coveredS11(i) it and transmits
Parameter FS21(i):
(i=1,2 ...).
Data de-noising module denoising method described further includes:
(1) for original Coherent Noise in GPR Record, the data in its preceding k nanosecond is intercepted and determine the standard track of data, k is big
In 0 real number;
(2) it for any scanning road, calculates cross-correlation function within preceding k nanosecond of the scanning road and standard track and takes it most
Big value;
(3) all scanning roads are traversed according to step (2), obtains the first arrival being made of each scanning road cross-correlation function maximum value
Curve L, and then calculate the high D of root mean square of first arrival curve L;
(4) singular value number N corresponding with coherent noise is calculated according to the following formula;N=ROUND { 0.2634D+
1.3086 } wherein: ROUND { } is the function that rounds up;
(5) singular value decomposition is carried out to original Coherent Noise in GPR Record, the top n data component after abnormal value elimination decomposition,
Denoising process is completed after remaining data ingredient is recombined.
Further, modeling module construction method includes:
The course of work of the ground penetrating radar system under Near Field is abstracted, and model parameter is set, depicts thunder
Up to the work block diagram of system;
Using the work block diagram of radar system, the equation of the Ground Penetrating Radar working principle described based on model parameter is obtained;
Model parameter calibration experiments are carried out, based on the equation of Ground Penetrating Radar working principle, using the measured value of different height
Construct one group of equation group about model parameter;
The equation group about model parameter is solved, the model of ground penetrating radar system is obtained;
The model parameter is the characterisitic parameter of radar system, which includes: the characterisitic parameter of transmitting antenna, connects
It receives more between parameter, receiving antenna and the target for the effect that intercouples between the characterisitic parameter of antenna, transmitting antenna and receiving antenna
The parameter of secondary reflection effect;
The equation of the Ground Penetrating Radar working principle based on model parameter description, meets:
Wherein, S (ω) indicates transmitting signal and receives the ratio of signal;A (ω) indicates a certain frequency of radar emission
Signal;B (ω) indicates the signal for the frequency that radar receives;ω indicates angular frequency;T0(ω) indicates transmitting antenna and receives
Intercouple the parameter of effect between antenna;TsIndicate the characterisitic parameter vector of receiving antenna;IN indicates unit matrix;TiIndicate hair
Penetrate the characterisitic parameter vector of antenna;G, G0It is Green's function matrix of the electric dipole in layered medium, wherein G indicates electricity
The launch point and receiving point of magnetic wave, which exist, to be biased, referred to as single biasing Green's function;G0Indicate the launch point and receiving point of electromagnetic wave
It is overlapped, referred to as zero offset Green's function;Rs indicates the parameter matrix of the multiple reflections effect between receiving antenna and target.
The Ground Penetrating Radar quick detection road engineering concrete mats are realized another object of the present invention is to provide a kind of
The computer program of layer thickness analysis method.
Another object of the present invention is to provide a kind of terminal, the terminal, which is at least carried, realizes that the Ground Penetrating Radar is quick
Detect the controller of road engineering concrete cushion thickness analysis method.
Another object of the present invention is to provide a kind of computer readable storage mediums, including instruction, when it is in computer
When upper operation, so that computer executes the Ground Penetrating Radar quick detection road engineering concrete cushion thickness analysis method.
Another object of the present invention is to provide a kind of Ground Penetrating Radar quick detection road engineering concrete cushion thickness point
Analysis system includes:
Parameter setting module: connecting with main control module, for radar parameter to be arranged;
Data acquisition module: connecting with main control module, for detecting Project Areas road engineering radar image, passes through live road
Road coring situation is modified velocity of electromagnetic wave, and the velocity of electromagnetic wave for obtaining the concrete under specific proportion is the data of v2;
Main control module: with parameter setting module, data acquisition module, data de-noising module, thickness measure module, modeling mould
Block, data memory module, display module connection, work normally for controlling modules;
Data de-noising module: connecting with main control module, for carrying out denoising to Coherent Noise in GPR Record;
Thickness measure module: connecting with main control module, for being handled in EKKO-Project software by gain declining profit
Surface layer and the bed course of Project Areas road are distinguished, road evenness and concrete cushion thickness h at velocity of electromagnetic wave v2 are measured;
Modeling module: connecting with main control module, for constructing the model of ground penetrating radar system;
Data memory module: connecting with main control module, for storing the data information of Ground Penetrating Radar measurement;
Display module: connecting with main control module, for showing the thickness data information of Ground Penetrating Radar measurement by display screen.
Another object of the present invention is to provide a kind of equipment for detecting road engineering concrete cushion thickness, the detections
The equipment of road engineering concrete cushion thickness at least carries the Ground Penetrating Radar quick detection road engineering concrete cushion
Thickness Analysis system.
Advantages of the present invention and good effect are as follows:
Coherent Noise in GPR Record is utilized singular value decomposition at a series of unusual with difference by the present invention by data de-noising module
It is worth corresponding data component, in order to individually handle data according to unlike signal or noise type;Then Ground Penetrating Radar number is utilized
According to first arrival lineups approximation detected object surface undulation situation, while the root mean square for calculating first arrival lineups is high;It finally will be equal
Root height is as the foundation of singular value corresponding with coherent noise is judged automatically, to reach the mesh for automatically removing coherent noise
, the trouble for needing artificially to judge singular value based on SVD decomposition method is avoided, surprise is directly judged according to root mean square high parameter
Different value, it is not only accurate but also time saving;Pass through the course of work and feature of the modeling module by ground penetrating radar system under Near Field simultaneously
It is abstracted, the sending and receiving antenna of ground penetrating radar system is indicated that the gain of antenna, loss, phase are prolonged using one group of electric dipole
The characteristics such as the coupling late and between sending and receiving antenna establish physical parameter, and retouch using one group function representation related with frequency
The work block diagram for stating radar system is established according to the physical parameter relation equation measured under different height and is used to solving model
The objective function of parameter, and solved using optimizing, considering exist about in solution procedure between some physical parameters
Beam relationship simplifies the number of the unknown physical parameter of solution accordingly, and initial value is according to the different geometry of electric dipole point
Cloth uses different calculations, and the influence that obtained Ground Penetrating Radar model considers antenna to radar signal being capable of accurate mould
Communication process of the radar signal in layered medium, the scope of application covering far field of the model and near field are drawn up, to utilize spy ground
Radar carries out accurate detection to the characteristic of layered medium and provides strong theoretical model.
Road engineering radar image in Project Areas is detected using profile method by data acquisition module, passes through site road coring
Situation is modified velocity of electromagnetic wave, and the velocity of electromagnetic wave for obtaining the concrete under specific proportion is v2 data;The item of detection
The signal model of mesh area road engineering radar image time-frequency overlapping MASK are as follows:
Wherein, N is the signal component number of time-frequency overlapped signal, and n (t) is additive white Gaussian noise, siIt (t) is time-frequency weight
The signal component of folded signal, is expressed asA in formulaiIndicate signal
The amplitude of component, ai(m) indicate that the symbol of signal component, p (t) indicate molding filter function, TiIndicate the code of signal component
First period, fciIndicate the carrier frequency of signal component,Indicate the phase of signal component;And after carrying out nonlinear transformation, it can obtain
The velocity of electromagnetic wave v2 accurate data of concrete is obtained, is improved compared with the prior art very much.
Main control module of the present invention carries out denoising to Coherent Noise in GPR Record by data de-noising module;Pass through thickness measure
Module distinguishes surface layer and the bed course of Project Areas road in EKKO-Project software by the processing of gain declining profit, measures in electromagnetism
Road evenness and concrete cushion thickness h under wave velocity v2;
In the DATA REASONING of concrete cushion thickness h, thickness h error parameter need to be handled, utilize Fourier's inversion
It changing, the scattering parameter that frequency domain is measured transforms to time domain, obtains time domain impulse response, and reflection parameters correspond to Time Domain Reflectometry response,
Configured transmission corresponds to time-domain transmission response;
According to the position of the first two pulse on a timeline in reflex response and transmission response, four time domain choosings are constructed respectively
Pass function;
Using gate function, the reflex response and transmission response of time domain are gated, extract respectively reflex response and
The first two pulse in transmission response;
Time domain impulse after time domain gating is passed through into Fourier transformation respectively, obtains frequency domain gated data;The above method
Operation can get concrete cushion exact thickness h, and data accuracy improves a lot compared with the prior art.
Detailed description of the invention
Fig. 1 is that the present invention implements the Ground Penetrating Radar quick detection road engineering concrete cushion thickness analysis method stream provided
Cheng Tu.
Fig. 2 is that the present invention implements the Ground Penetrating Radar quick detection road engineering concrete cushion Thickness Analysis system knot provided
Structure block diagram.
In figure: 1, parameter setting module;2, data acquisition module;3, main control module;4, data de-noising module;5, thickness is surveyed
Measure module;6, modeling module;7, data memory module;8, display module.
Specific embodiment
In order to make the objectives, technical solutions, and advantages of the present invention clearer, with reference to embodiments, to the present invention
It is further elaborated.It should be appreciated that the specific embodiments described herein are merely illustrative of the present invention, it is not used to
Limit the present invention.
With reference to the accompanying drawing and specific embodiment is further described application principle of the invention.
As shown in Figure 1, Ground Penetrating Radar quick detection road engineering concrete cushion thickness provided in an embodiment of the present invention point
Analysis system and analysis method the following steps are included:
Radar parameter: radar frequency 250MHz is arranged by parameter setting module in S101;When window 36ns;Velocity of electromagnetic wave v1
=0.1m/ns;
S102 detects Project Areas road engineering radar image using profile method by data acquisition module, passes through site road
Coring situation is modified velocity of electromagnetic wave, and the velocity of electromagnetic wave for obtaining the concrete under specific proportion is v2 data;
S103 main control module carries out denoising to Coherent Noise in GPR Record by data de-noising module;Pass through thickness measure mould
Block distinguishes surface layer and the bed course of Project Areas road in EKKO-Project software by the processing of gain declining profit, measures in electromagnetic wave
Road evenness and concrete cushion thickness h under speed v2;
S104 constructs the model of ground penetrating radar system by modeling module;Ground Penetrating Radar is stored by data memory module to survey
The data information of amount;
S105 shows the thickness data information of Ground Penetrating Radar measurement by display module.
As shown in Fig. 2, Ground Penetrating Radar quick detection road engineering concrete cushion Thickness Analysis system provided by the invention
It include: parameter setting module 1, data acquisition module 2, main control module 3, data de-noising module 4, thickness measure module 5, modeling mould
Block 6, data memory module 7, display module 8.
Parameter setting module 1: it is connect with main control module 3, for radar parameter to be arranged.Wherein, radar frequency 250MHz, when
Window 36ns, velocity of electromagnetic wave v1=0.1m/ns;
Data acquisition module 2: connecting with main control module 3, for detecting Project Areas road engineering radar image, passes through scene
Road coring situation is modified velocity of electromagnetic wave, and the velocity of electromagnetic wave for obtaining the concrete under specific proportion is v2 data;
Main control module 3: with parameter setting module 1, data acquisition module 2, data de-noising module 4, thickness measure module 5,
Modeling module 6, data memory module 7, display module 8 connect, and work normally for controlling modules;
Data de-noising module 4: connecting with main control module 3, for carrying out denoising to Coherent Noise in GPR Record;
Thickness measure module 5: connecting with main control module 3, for passing through at gain declining profit in EKKO-Project software
Reason distinguishes surface layer and the bed course of Project Areas road, measures road evenness and concrete cushion thickness at velocity of electromagnetic wave v2
h;
Modeling module 6: connecting with main control module 3, for constructing the model of ground penetrating radar system;
Data memory module 7: connecting with main control module 3, for storing the data information of Ground Penetrating Radar measurement;
Display module 8: connecting with main control module 3, for showing that the thickness data of Ground Penetrating Radar measurement is believed by display screen
Breath.
Data de-noising module 4 provided by the invention, denoising method is as follows:
(1) for original Coherent Noise in GPR Record, the data in its preceding k nanosecond is intercepted and determine the standard track of data, k is big
In 0 real number;
(2) it for any scanning road, calculates cross-correlation function within preceding k nanosecond of the scanning road and standard track and takes it most
Big value;
(3) all scanning roads are traversed according to step (2), obtains the first arrival being made of each scanning road cross-correlation function maximum value
Curve L, and then calculate the high D of root mean square of first arrival curve L;
(4) singular value number N corresponding with coherent noise is calculated according to the following formula;N=ROUND { 0.2634D+
1.3086 } wherein: ROUND { } is the function that rounds up;
(5) singular value decomposition is carried out to original Coherent Noise in GPR Record, the top n data component after abnormal value elimination decomposition,
Denoising process is completed after remaining data ingredient is recombined.
6 construction method of modeling module provided by the invention is as follows:
Firstly, being abstracted to the course of work of the ground penetrating radar system under Near Field, and model parameter is set, is described
The work block diagram of radar system out;
Secondly, obtaining the Ground Penetrating Radar working principle described based on model parameter using the work block diagram of radar system
Equation;
Then, model parameter calibration experiments are carried out, based on the equation of Ground Penetrating Radar working principle, using the survey of different height
Magnitude constructs one group of equation group about model parameter;
Finally, solving the equation group about model parameter, the model of ground penetrating radar system is obtained.
Model parameter provided by the invention is the characterisitic parameter of radar system, which includes: the spy of transmitting antenna
Intercouple parameter, receiving antenna and the mesh of effect between property parameter, the characterisitic parameter of receiving antenna, transmitting antenna and receiving antenna
The parameter that multiple reflections act between mark.
The equation of Ground Penetrating Radar working principle provided by the invention based on model parameter description, meets:
Wherein, S (ω) indicates transmitting signal and receives the ratio of signal;A (ω) indicates a certain frequency of radar emission
Signal;B (ω) indicates the signal for the frequency that radar receives;ω indicates angular frequency;T0(ω) indicates transmitting antenna and receives
Intercouple the parameter of effect between antenna;TsIndicate the characterisitic parameter vector of receiving antenna;IN indicates unit matrix;TiIndicate hair
Penetrate the characterisitic parameter vector of antenna;G, G0It is Green's function matrix of the electric dipole in layered medium, wherein G indicates electricity
The launch point and receiving point of magnetic wave, which exist, to be biased, referred to as single biasing Green's function;G0Indicate the launch point and receiving point of electromagnetic wave
It is overlapped, referred to as zero offset Green's function;Rs indicates the parameter matrix of the multiple reflections effect between receiving antenna and target.
The invention will be further described combined with specific embodiments below.
The embodiment of the present invention provides Ground Penetrating Radar quick detection road engineering concrete cushion thickness analysis method, comprising:
Road engineering radar image in Project Areas is detected using profile method by data acquisition module, passes through site road coring
Situation is modified velocity of electromagnetic wave, and the velocity of electromagnetic wave for obtaining the concrete under specific proportion is v2 data;It specifically includes:
1) signal model of the Project Areas road engineering radar image time-frequency overlapping MASK detected are as follows:
Wherein, N is the signal component number of time-frequency overlapped signal, and n (t) is additive white Gaussian noise, siIt (t) is time-frequency weight
The signal component of folded signal, is expressed asA in formulaiIndicate signal
The amplitude of component, ai(m) indicate that the symbol of signal component, p (t) indicate molding filter function, TiIndicate the code of signal component
First period, fciIndicate the carrier frequency of signal component,Indicate the phase of signal component;
It is rightNonlinear transformation is carried out,
It obtains: y (t) are as follows:
Y (t)=x (t)+n (t);
Wherein, x (t) is digital modulation signals, and n (t) is the impulsive noise of obedience standard S α S distribution, the parsing shape of x (t)
Formula indicates are as follows:
Wherein, N is sampling number, anFor the information symbol of transmission, in MASK signal, an=0,1,2 ..., M-1, M are
Order of modulation, an=ej2πε/M, ε=0,1,2 ..., M-1, g (t) expression rectangle molding pulse, TbIndicate symbol period, fcIt indicates
Carrier frequency, carrier wave initial phaseIt is the equally distributed random number in [0,2 π];
2) the multipath space of n signal is constructed are as follows:
Wherein,Q is sampling number, and K is maximum delay, by
Maximum detectable range Rmax/ c is obtained, wherein xreciIt (t) is reference signal, RmaxFor maximum detectable range, c is the light velocity;
3) then inhibit direct wave and its multipath using principle of least square method, min will be sought | | Ssur-Xref·α||2Conversion
To askIt obtains:
Substitute into αestim, it solves:
Wherein, SsurFor echo channel signal, α is adaptive weight, αestimFor the estimated value of α,For XrefTransposition,
SotherFor remaining echo and noise final in echo channel;
Main control module carries out denoising to Coherent Noise in GPR Record by data de-noising module;Existed by thickness measure module
Surface layer and the bed course for distinguishing Project Areas road in EKKO-Project software by the processing of gain declining profit, measure in velocity of electromagnetic wave
Road evenness and concrete cushion thickness h under v2;
In the DATA REASONING of concrete cushion thickness h, thickness h error parameter need to be handled, be specifically included:
Using inverse Fourier transform, the scattering parameter that frequency domain is measured transforms to time domain, obtains time domain impulse response, reflection
Parameter corresponds to Time Domain Reflectometry response, and configured transmission corresponds to time-domain transmission response;
According to the position of the first two pulse on a timeline in reflex response and transmission response, four time domain choosings are constructed respectively
Pass function;
Using gate function, the reflex response and transmission response of time domain are gated, extract respectively reflex response and
The first two pulse in transmission response;
Time domain impulse after time domain gating is passed through into Fourier transformation respectively, obtains frequency domain gated data;
Information comprising detected concrete bed course in EKKO-Project software in frequency domain gated data, utilizes obtained letter
Breath, according to formula construction compensation factor Fcf(i);Using following formula, compensation factor F is constructedcf(i):
(i=1,2 ...)
G1(i)~G4It (i) is the frequency domain gated data obtained in upper step;
R (i) is ratio factor;
Fcf(i) compensation factor;
Using occlusion compensation formula, the reflection parameters F of error is not coveredS11(i) and configured transmission FS21(i);
The model of ground penetrating radar system is constructed by modeling module;Ground Penetrating Radar measurement is stored by data memory module
Data information;
The Ground Penetrating Radar quick detection road engineering concrete cushion thickness analysis method further include:
By parameter setting module, radar parameter: radar frequency 250MHz is set;When window 36ns;Velocity of electromagnetic wave v1=
0.1m/ns;
The thickness data information of Ground Penetrating Radar measurement is shown by display module.
Using following occlusion compensation formula, the reflection parameters F of error is not coveredS11(i) and configured transmission FS21
(i):
(i=1,2 ...).
In the above-described embodiments, can come wholly or partly by software, hardware, firmware or any combination thereof real
It is existing.When using entirely or partly realizing in the form of a computer program product, the computer program product include one or
Multiple computer instructions.When loading on computers or executing the computer program instructions, entirely or partly generate according to
Process described in the embodiment of the present invention or function.The computer can be general purpose computer, special purpose computer, computer network
Network or other programmable devices.The computer instruction may be stored in a computer readable storage medium, or from one
Computer readable storage medium is transmitted to another computer readable storage medium, for example, the computer instruction can be from one
A web-site, computer, server or data center pass through wired (such as coaxial cable, optical fiber, Digital Subscriber Line (DSL)
Or wireless (such as infrared, wireless, microwave etc.) mode is carried out to another web-site, computer, server or data center
Transmission).The computer-readable storage medium can be any usable medium or include one that computer can access
The data storage devices such as a or multiple usable mediums integrated server, data center.The usable medium can be magnetic Jie
Matter, (for example, floppy disk, hard disk, tape), optical medium (for example, DVD) or semiconductor medium (such as solid state hard disk Solid
State Disk (SSD)) etc..
The foregoing is merely illustrative of the preferred embodiments of the present invention, is not intended to limit the invention, all in essence of the invention
Made any modifications, equivalent replacements, and improvements etc., should all be included in the protection scope of the present invention within mind and principle.
Claims (10)
1. a kind of Ground Penetrating Radar quick detection road engineering concrete cushion thickness analysis method, which is characterized in that the spy ground
Radar quick detection road engineering concrete cushion thickness analysis method includes:
Road engineering radar image in Project Areas is detected using profile method by data acquisition module, passes through site road coring situation
Velocity of electromagnetic wave is modified, the velocity of electromagnetic wave for obtaining the concrete under specific proportion is v2 data;It specifically includes:
1) signal model of the Project Areas road engineering radar image time-frequency overlapping MASK detected are as follows:
Wherein, N is the signal component number of time-frequency overlapped signal, and n (t) is additive white Gaussian noise, siIt (t) is time-frequency overlapping letter
Number signal component, be expressed asA in formulaiIndicate signal component
Amplitude, ai(m) indicate that the symbol of signal component, p (t) indicate molding filter function, TiIndicate the symbol week of signal component
Phase, fciIndicate the carrier frequency of signal component,Indicate the phase of signal component;
It is rightNonlinear transformation is carried out,
It obtains: y (t) are as follows:
Y (t)=x (t)+n (t);
Wherein, x (t) is digital modulation signals, and n (t) is the impulsive noise of obedience standard S α S distribution, the analytical form table of x (t)
It is shown as:
Wherein, N is sampling number, anFor the information symbol of transmission, in MASK signal, an=0,1,2 ..., M-1, M are modulation
Order, an=ej2πε/M, ε=0,1,2 ..., M-1, g (t) expression rectangle molding pulse, TbIndicate symbol period, fcIndicate carrier wave
Frequency, carrier wave initial phaseIt is the equally distributed random number in [0,2 π];
2) the multipath space of n signal is constructed are as follows:
Xref=[Xref1 Xref2...Xrefn];
Wherein,Q is sampling number, and K is maximum delay, by maximum
Detection range Rmax/ c is obtained, wherein xreciIt (t) is reference signal, RmaxFor maximum detectable range, c is the light velocity;
3) then inhibit direct wave and its multipath using principle of least square method, min will be sought | | Ssur-Xref·α||2It is converted into and asksIt obtains:
Substitute into αestim, it solves:
Wherein, SsurFor echo channel signal, α is adaptive weight, αestimFor the estimated value of α,For XrefTransposition, Sother
For remaining echo and noise final in echo channel;
Main control module carries out denoising to Coherent Noise in GPR Record by data de-noising module;Existed by thickness measure module
Surface layer and the bed course for distinguishing Project Areas road in EKKO-Project software by the processing of gain declining profit, measure in velocity of electromagnetic wave
Road evenness and concrete cushion thickness h under v2;
In the DATA REASONING of concrete cushion thickness h, thickness h error parameter need to be handled, be specifically included:
Using inverse Fourier transform, the scattering parameter that frequency domain is measured transforms to time domain, obtains time domain impulse response, reflection parameters
Corresponding Time Domain Reflectometry response, configured transmission correspond to time-domain transmission response;
According to the position of the first two pulse on a timeline in reflex response and transmission response, four time domain gating letters are constructed respectively
Number;
Using gate function, the reflex response and transmission response of time domain are gated, extract reflex response and transmission respectively
The first two pulse in response;
Time domain impulse after time domain gating is passed through into Fourier transformation respectively, obtains frequency domain gated data;
Information comprising detected concrete bed course in EKKO-Project software in frequency domain gated data, using obtained information,
According to formula construction compensation factor Fcf(i);Using following formula, compensation factor F is constructedcf(i):
(i=1,2 ...)
G1(i)~G4It (i) is the frequency domain gated data obtained in upper step;
R (i) is ratio factor;
Fcf(i) compensation factor;
Using occlusion compensation formula, the reflection parameters F of error is not coveredS11(i) and configured transmission FS21(i);
The model of ground penetrating radar system is constructed by modeling module;The data of Ground Penetrating Radar measurement are stored by data memory module
Information.
2. Ground Penetrating Radar quick detection road engineering concrete cushion thickness analysis method as described in claim 1, feature exist
In the Ground Penetrating Radar quick detection road engineering concrete cushion thickness analysis method further include:
By parameter setting module, radar parameter: radar frequency 250MHz is set;When window 36ns;Velocity of electromagnetic wave v1=0.1m/
ns;
The thickness data information of Ground Penetrating Radar measurement is shown by display module.
3. Ground Penetrating Radar quick detection road engineering concrete cushion thickness analysis method as described in claim 1, feature exist
In not covered the reflection parameters F of error using following occlusion compensation formulaS11(i) and configured transmission FS21(i):
(i=1,2 ...).
4. Ground Penetrating Radar quick detection road engineering concrete cushion thickness analysis method as described in claim 1, feature exist
In the data de-noising module denoising method includes:
(1) for original Coherent Noise in GPR Record, the data in its preceding k nanosecond is intercepted and determine the standard track of data, k is greater than 0
Real number;
(2) it for any scanning road, calculates cross-correlation function within preceding k nanosecond of the scanning road and standard track and takes its maximum
Value;
(3) all scanning roads are traversed according to step (2), obtains the first arrival curve being made of each scanning road cross-correlation function maximum value
L, and then calculate the high D of root mean square of first arrival curve L;
(4) singular value number N corresponding with coherent noise is calculated according to the following formula;N=ROUND { 0.2634D+1.3086 }
Wherein: ROUND { } is the function that rounds up;
(5) singular value decomposition is carried out to original Coherent Noise in GPR Record, the top n data component after abnormal value elimination decomposition will remain
Remaining data component completes denoising process after recombining.
5. Ground Penetrating Radar quick detection road engineering concrete cushion thickness analysis method as described in claim 1, feature exist
In modeling module construction method includes:
The course of work of the ground penetrating radar system under Near Field is abstracted, and model parameter is set, depicts radar system
The work block diagram of system;
Using the work block diagram of radar system, the equation of the Ground Penetrating Radar working principle described based on model parameter is obtained;
Model parameter calibration experiments are carried out, based on the equation of Ground Penetrating Radar working principle, are constructed using the measured value of different height
One group of equation group about model parameter;
The equation group about model parameter is solved, the model of ground penetrating radar system is obtained;
The model parameter is the characterisitic parameter of radar system, which includes: the characterisitic parameter of transmitting antenna, receives day
It is repeatedly anti-between parameter, receiving antenna and the target of effect that intercouple between the characterisitic parameter of line, transmitting antenna and receiving antenna
Penetrate the parameter of effect;
The equation of the Ground Penetrating Radar working principle based on model parameter description, meets:
Wherein, S (ω) indicates transmitting signal and receives the ratio of signal;A (ω) indicates the signal of a certain frequency of radar emission;
B (ω) indicates the signal for the frequency that radar receives;ω indicates angular frequency;T0(ω) is indicated between transmitting antenna and receiving antenna
Intercouple the parameter of effect;TsIndicate the characterisitic parameter vector of receiving antenna;IN indicates unit matrix;TiIndicate transmitting antenna
Characterisitic parameter vector;G, G0It is Green's function matrix of the electric dipole in layered medium, wherein G indicates electromagnetic wave
Launch point and receiving point, which exist, to be biased, referred to as single biasing Green's function;G0It indicates that the launch point of electromagnetic wave is overlapped with receiving point, claims
For zero offset Green's function;Rs indicates the parameter matrix of the multiple reflections effect between receiving antenna and target.
6. Ground Penetrating Radar quick detection road engineering concrete cushion thickness described in a kind of realization Claims 1 to 5 any one
The computer program of analysis method.
7. a kind of terminal, which is characterized in that the terminal, which is at least carried, visits land mine described in realization Claims 1 to 5 any one
Up to the controller of quick detection road engineering concrete cushion thickness analysis method.
8. a kind of computer readable storage medium, including instruction, when run on a computer, so that computer is executed as weighed
Benefit requires Ground Penetrating Radar quick detection road engineering concrete cushion thickness analysis method described in 1-5 any one.
9. Ground Penetrating Radar quick detection road engineering concrete cushion thickness described in a kind of realization Claims 1 to 5 any one
The Ground Penetrating Radar quick detection road engineering concrete cushion Thickness Analysis system of analysis method, which is characterized in that the spy ground
Radar quick detection road engineering concrete cushion Thickness Analysis system includes:
Parameter setting module: connecting with main control module, for radar parameter to be arranged;
Data acquisition module: connecting with main control module, for detecting Project Areas road engineering radar image, is taken by site road
Core situation is modified velocity of electromagnetic wave, and the velocity of electromagnetic wave for obtaining the concrete under specific proportion is the data of v2;
Main control module: with parameter setting module, data acquisition module, data de-noising module, thickness measure module, modeling module,
Data memory module, display module connection, work normally for controlling modules;
Data de-noising module: connecting with main control module, for carrying out denoising to Coherent Noise in GPR Record;
Thickness measure module: connecting with main control module, distinguishes for passing through the processing of gain declining profit in EKKO-Project software
The surface layer of Project Areas road and bed course measure road evenness and concrete cushion thickness h at velocity of electromagnetic wave v2;
Modeling module: connecting with main control module, for constructing the model of ground penetrating radar system;
Data memory module: connecting with main control module, for storing the data information of Ground Penetrating Radar measurement;
Display module: connecting with main control module, for showing the thickness data information of Ground Penetrating Radar measurement by display screen.
10. a kind of equipment for detecting road engineering concrete cushion thickness, which is characterized in that the detection road engineering concrete
The equipment of cushion thickness at least carries Ground Penetrating Radar quick detection road engineering concrete cushion thickness as claimed in claim 9 point
Analysis system.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811326015.1A CN109541586A (en) | 2018-11-08 | 2018-11-08 | A kind of Ground Penetrating Radar quick detection road engineering concrete cushion thickness analysis method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811326015.1A CN109541586A (en) | 2018-11-08 | 2018-11-08 | A kind of Ground Penetrating Radar quick detection road engineering concrete cushion thickness analysis method |
Publications (1)
Publication Number | Publication Date |
---|---|
CN109541586A true CN109541586A (en) | 2019-03-29 |
Family
ID=65845207
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811326015.1A Pending CN109541586A (en) | 2018-11-08 | 2018-11-08 | A kind of Ground Penetrating Radar quick detection road engineering concrete cushion thickness analysis method |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109541586A (en) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110258267A (en) * | 2019-05-27 | 2019-09-20 | 苏州奇摩智能科技有限公司 | A kind of control method, apparatus and system determining asphalt paving thickness based on Ground Penetrating Radar |
CN111595293A (en) * | 2020-05-29 | 2020-08-28 | 山东交通学院 | Surface deformation monitoring method and system based on multi-source monitoring data fusion |
CN111624596A (en) * | 2020-04-21 | 2020-09-04 | 中建路桥集团有限公司 | Road and bridge road acceptance detection system |
CN114061432A (en) * | 2022-01-17 | 2022-02-18 | 中大检测(湖南)股份有限公司 | Geological radar detection system and method |
CN114166111A (en) * | 2021-11-11 | 2022-03-11 | 北京戍宁信息技术有限公司 | Method and system for monitoring synthetic aperture radar of foundation by wide-angle deformation |
Citations (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20090167589A1 (en) * | 2006-06-29 | 2009-07-02 | Igor Bausov | Earth-penetrating radar with inherent near-field rejection |
WO2014153263A1 (en) * | 2013-03-14 | 2014-09-25 | Robert Ernest Troxler | Systems and methods for asphalt density and soil moisture measurements using ground penetrating radar |
CN104133212A (en) * | 2014-07-09 | 2014-11-05 | 中国电子科技集团公司第二十二研究所 | Radar device and method used for real-time detection of road structure layer thickness |
CN104931912A (en) * | 2015-06-12 | 2015-09-23 | 西安电子科技大学 | Time domain compensation method of vector network analyzer |
CN105866750A (en) * | 2016-05-23 | 2016-08-17 | 西安电子科技大学 | Method for detecting multiple GPS (global positioning system) satellite weak echo signals |
CN105939188A (en) * | 2016-06-13 | 2016-09-14 | 西安电子科技大学 | Carrier frequency estimation method for underlay spectrum sharing time-frequency overlapped MASK signal |
CN106501940A (en) * | 2016-12-12 | 2017-03-15 | 湖南工业大学 | A kind of height degree of immersing Head-mounted display control system |
CN106646459A (en) * | 2017-03-07 | 2017-05-10 | 山东农业大学 | Method for quickly detecting thickness of road engineering concrete cushion by using ground penetrating radar |
CN107450054A (en) * | 2017-07-14 | 2017-12-08 | 浙江省交通规划设计研究院 | A kind of adaptive Coherent Noise in GPR Record denoising method |
CN108061920A (en) * | 2017-12-07 | 2018-05-22 | 中国科学院电子学研究所 | The method of Ground Penetrating Radar modeling |
CN108270625A (en) * | 2018-01-30 | 2018-07-10 | 河南质量工程职业学院 | A kind of data calculating control system based on cloud service platform |
CN108470389A (en) * | 2018-02-02 | 2018-08-31 | 武汉城市职业学院 | A kind of intelligent express delivery based on NFC technique is from receive-transmit system |
CN108761446A (en) * | 2018-04-09 | 2018-11-06 | 中国科学院电子学研究所 | The modeling method of frequency stepping Ground Penetrating Radar |
-
2018
- 2018-11-08 CN CN201811326015.1A patent/CN109541586A/en active Pending
Patent Citations (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20090167589A1 (en) * | 2006-06-29 | 2009-07-02 | Igor Bausov | Earth-penetrating radar with inherent near-field rejection |
WO2014153263A1 (en) * | 2013-03-14 | 2014-09-25 | Robert Ernest Troxler | Systems and methods for asphalt density and soil moisture measurements using ground penetrating radar |
CN104133212A (en) * | 2014-07-09 | 2014-11-05 | 中国电子科技集团公司第二十二研究所 | Radar device and method used for real-time detection of road structure layer thickness |
CN104931912A (en) * | 2015-06-12 | 2015-09-23 | 西安电子科技大学 | Time domain compensation method of vector network analyzer |
CN105866750A (en) * | 2016-05-23 | 2016-08-17 | 西安电子科技大学 | Method for detecting multiple GPS (global positioning system) satellite weak echo signals |
CN105939188A (en) * | 2016-06-13 | 2016-09-14 | 西安电子科技大学 | Carrier frequency estimation method for underlay spectrum sharing time-frequency overlapped MASK signal |
CN106501940A (en) * | 2016-12-12 | 2017-03-15 | 湖南工业大学 | A kind of height degree of immersing Head-mounted display control system |
CN106646459A (en) * | 2017-03-07 | 2017-05-10 | 山东农业大学 | Method for quickly detecting thickness of road engineering concrete cushion by using ground penetrating radar |
CN107450054A (en) * | 2017-07-14 | 2017-12-08 | 浙江省交通规划设计研究院 | A kind of adaptive Coherent Noise in GPR Record denoising method |
CN108061920A (en) * | 2017-12-07 | 2018-05-22 | 中国科学院电子学研究所 | The method of Ground Penetrating Radar modeling |
CN108270625A (en) * | 2018-01-30 | 2018-07-10 | 河南质量工程职业学院 | A kind of data calculating control system based on cloud service platform |
CN108470389A (en) * | 2018-02-02 | 2018-08-31 | 武汉城市职业学院 | A kind of intelligent express delivery based on NFC technique is from receive-transmit system |
CN108761446A (en) * | 2018-04-09 | 2018-11-06 | 中国科学院电子学研究所 | The modeling method of frequency stepping Ground Penetrating Radar |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110258267A (en) * | 2019-05-27 | 2019-09-20 | 苏州奇摩智能科技有限公司 | A kind of control method, apparatus and system determining asphalt paving thickness based on Ground Penetrating Radar |
CN111624596A (en) * | 2020-04-21 | 2020-09-04 | 中建路桥集团有限公司 | Road and bridge road acceptance detection system |
CN111595293A (en) * | 2020-05-29 | 2020-08-28 | 山东交通学院 | Surface deformation monitoring method and system based on multi-source monitoring data fusion |
CN114166111A (en) * | 2021-11-11 | 2022-03-11 | 北京戍宁信息技术有限公司 | Method and system for monitoring synthetic aperture radar of foundation by wide-angle deformation |
CN114166111B (en) * | 2021-11-11 | 2023-08-04 | 北京戍宁信息技术有限公司 | Wide-angle deformation monitoring foundation synthetic aperture radar method and system |
CN114061432A (en) * | 2022-01-17 | 2022-02-18 | 中大检测(湖南)股份有限公司 | Geological radar detection system and method |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109541586A (en) | A kind of Ground Penetrating Radar quick detection road engineering concrete cushion thickness analysis method | |
US11546068B2 (en) | Ray tracing technique for wireless channel measurements | |
CN106254010B (en) | A kind of time-varying ocean channel modeling method | |
CN105137424B (en) | Real beam scanning radar angle ultra-resolution method under a kind of clutter background | |
CN103139907A (en) | Indoor wireless positioning method by utilizing fingerprint technique | |
CN104360332A (en) | Atmospheric phase screen extraction method based on ground-based SAR (synthetic aperture radar) interference | |
CN105445701A (en) | Mono-pulse angle estimation method for DDMA-MIMO radar target | |
Klysz et al. | Evaluation of dielectric properties of concrete by a numerical FDTD model of a GPR coupled antenna—parametric study | |
CN104270713A (en) | Passive type moving target track mapping method based on compressed sensing | |
CN102662195A (en) | Underground cavity detection system | |
CN108089147A (en) | Improved shortwave unit localization method | |
CN105372507A (en) | Real-time monitoring method and system for electromagnetic environment | |
CN105372506A (en) | Variable-scale rasterization method and system for regional electromagnetic environment calculation | |
CN104678392A (en) | Foundation SAR (Synthetic Aperture Radar) atmosphere delaying correction method based on stable pixel weighting | |
CN109752715A (en) | A kind of SAR data perfect diffuser detection method and device | |
CN107607945A (en) | A kind of scanning radar forword-looking imaging method based on spatial embedding mapping | |
Ji et al. | A novel Wi-Fi AP localization method using Monte Carlo path-loss model fitting simulation | |
CN107102329A (en) | Ka wave bands TOPS InSAR systems and its method of work | |
CN104808180A (en) | Steady waveform optimizing method for MIMO radar in clutter background | |
CN115113155A (en) | Airborne distributed aperture coherent synthetic radar testing and evaluating method | |
CN104808190B (en) | Improve the sane waveform design method of the worst parameter Estimation performance of MIMO radar | |
CN103176171B (en) | Simulation method of interferometric SAR (synthetic aperture radar) distributed targets | |
Arai et al. | Color radiomap interpolation for efficient fingerprint wifi-based indoor location estimation | |
CN103376443B (en) | Ground penetrating radar terrestrial interference detecting and fast eliminating method | |
CN102711043B (en) | RSS (Root Sum Square)-based passive motion tracking method suitable for sparse node distribution |
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 | ||
WD01 | Invention patent application deemed withdrawn after publication | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20190329 |