CN114624704B - Asphalt pavement crack type rapid detection method based on ground penetrating radar - Google Patents
Asphalt pavement crack type rapid detection method based on ground penetrating radar Download PDFInfo
- Publication number
- CN114624704B CN114624704B CN202210146001.1A CN202210146001A CN114624704B CN 114624704 B CN114624704 B CN 114624704B CN 202210146001 A CN202210146001 A CN 202210146001A CN 114624704 B CN114624704 B CN 114624704B
- Authority
- CN
- China
- Prior art keywords
- crack
- ground penetrating
- penetrating radar
- amplitude
- amplitude intensity
- 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
- 230000000149 penetrating effect Effects 0.000 title claims abstract description 54
- 239000010426 asphalt Substances 0.000 title claims abstract description 27
- 238000001514 detection method Methods 0.000 title claims abstract description 22
- 238000000034 method Methods 0.000 claims abstract description 30
- 238000010606 normalization Methods 0.000 claims abstract description 8
- 239000010410 layer Substances 0.000 claims description 55
- 230000008859 change Effects 0.000 claims description 36
- 239000002344 surface layer Substances 0.000 claims description 16
- 239000002699 waste material Substances 0.000 claims description 11
- 239000011384 asphalt concrete Substances 0.000 claims description 9
- 230000008878 coupling Effects 0.000 claims description 5
- 238000010168 coupling process Methods 0.000 claims description 5
- 238000005859 coupling reaction Methods 0.000 claims description 5
- 238000010586 diagram Methods 0.000 claims description 5
- 230000008030 elimination Effects 0.000 claims description 5
- 238000003379 elimination reaction Methods 0.000 claims description 5
- 239000000463 material Substances 0.000 claims description 5
- 239000011229 interlayer Substances 0.000 claims description 2
- 238000012545 processing Methods 0.000 abstract description 7
- 238000004364 calculation method Methods 0.000 abstract description 5
- 238000004836 empirical method Methods 0.000 abstract description 3
- 230000008569 process Effects 0.000 description 4
- 238000005070 sampling Methods 0.000 description 4
- 230000009471 action Effects 0.000 description 3
- 238000005259 measurement Methods 0.000 description 3
- 238000013461 design Methods 0.000 description 2
- 239000000758 substrate Substances 0.000 description 2
- 201000004569 Blindness Diseases 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000005336 cracking Methods 0.000 description 1
- 230000001186 cumulative effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 201000010099 disease Diseases 0.000 description 1
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 1
- 238000005553 drilling Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000012423 maintenance Methods 0.000 description 1
- 238000011002 quantification Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 239000002689 soil Substances 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
- 239000011800 void material Substances 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
- 238000004804 winding Methods 0.000 description 1
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
-
- E—FIXED CONSTRUCTIONS
- E01—CONSTRUCTION OF ROADS, RAILWAYS, OR BRIDGES
- E01C—CONSTRUCTION OF, OR SURFACES FOR, ROADS, SPORTS GROUNDS, OR THE LIKE; MACHINES OR AUXILIARY TOOLS FOR CONSTRUCTION OR REPAIR
- E01C23/00—Auxiliary devices or arrangements for constructing, repairing, reconditioning, or taking-up road or like surfaces
- E01C23/01—Devices or auxiliary means for setting-out or checking the configuration of new surfacing, e.g. templates, screed or reference line supports; Applications of apparatus for measuring, indicating, or recording the surface configuration of existing surfacing, e.g. profilographs
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N27/00—Investigating or analysing materials by the use of electric, electrochemical, or magnetic means
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Radar, Positioning & Navigation (AREA)
- Chemical & Material Sciences (AREA)
- General Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Civil Engineering (AREA)
- Architecture (AREA)
- Chemical Kinetics & Catalysis (AREA)
- Electrochemistry (AREA)
- Structural Engineering (AREA)
- Computer Networks & Wireless Communication (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Radar Systems Or Details Thereof (AREA)
- Road Repair (AREA)
Abstract
The invention belongs to the technical field of pavement crack detection, and particularly relates to a quick detection method for asphalt pavement crack types based on a ground penetrating radar, which can reduce theoretical calculation complexity and has high judgment accuracy. The method carries out crack type identification by calculating the amplitude value of each path of radar data within a certain time range below the interface reflection wave phase axis and carrying out normalization processing, and solves the problems of low precision and low efficiency when the conventional empirical method is used for judging a plurality of crack types of a long-distance section of a road.
Description
Technical Field
The invention belongs to the technical field of pavement crack detection, and particularly relates to a quick detection method for asphalt pavement crack types based on a ground penetrating radar, which can reduce theoretical calculation complexity and has high judgment accuracy.
Background
The asphalt pavement crack types mainly comprise fatigue cracks and reflection cracks, wherein the fatigue cracks generally cause fatigue cracking of the asphalt pavement layer under the repeated loading action of a vehicle, the fatigue cracks exist only in the asphalt pavement layer, and the reflection cracks are formed by reflecting transverse cracks generated by the semi-rigid base layer under the action of temperature and humidity change to the pavement. When external free water continuously enters the inside of the structural layer through pavement cracks and is accumulated at the tops of the roadbeds and the base layers, the base layer fine aggregate and roadbed soil are scoured and extruded under the repeated action of running load, so that the bottom of each layer is out of support to form local void, and the structural damage of the road is caused. The traditional detection method judges the type of the crack through core drilling sampling, so that the type and internal development condition of each crack cannot be accurately known, and the specific treatment cannot be carried out according to the type of each crack during maintenance design, so that the design has certain blindness.
In recent years, the ground penetrating radar has the advantages of continuous nondestructive, high speed, high precision and the like, and is widely applied to the field of fine nondestructive detection of deep diseases of highways. In the ground penetrating radar, in a single uniform ideal medium, the top and bottom points of the crack correspond to a diffraction wave hyperbolic wave group, so that the top and bottom points of the ground penetrating radar are easy to identify. However, when the ground penetrating radar propagates in the layered uniform medium, multiple diffraction wave hyperbola wave groups are formed at each medium interface passing through the ground penetrating radar, so that the diffraction wave hyperbola wave groups corresponding to the top end point and the bottom end point of the crack are similar in form, and the bottom end point of the crack cannot be judged. At present, the type judgment of the reflection crack mainly relies on an empirical method and a theoretical calculation method, but because the asphalt pavement belongs to a multi-phase discrete random medium layered structure, firstly, the definition of a hyperbola is interfered by a complex medium of the road structure, so that the judgment of experience is affected; secondly, the complexity of theoretical calculation is difficult to quickly and accurately detect and identify the crack type in road long section detection.
Disclosure of Invention
The invention aims to overcome the defects of the prior art, and provides a quick detection method for the type of asphalt pavement cracks based on a ground penetrating radar, which is used for solving the problems of low accuracy and complexity of theoretical calculation of the conventional experience judgment method.
The technical scheme adopted by the invention is as follows:
a method for rapidly detecting crack type of asphalt pavement based on a ground penetrating radar comprises the following steps:
Step 1), acquiring underground data of a crack by adopting a ground penetrating radar to obtain a waveform diagram of an electromagnetic wave time domain signal;
step 2), editing the collected ground penetrating radar data;
Step 3), determining interlayer interfaces between the asphalt surface layer and the base layer and between the base layer and the roadbed according to the change of reflection amplitude of electromagnetic waves among materials of each layer;
Step 4), carrying out statistics and normalization treatment on amplitude values of radar data of each crack position in a time range of 0.7-1.3 ns below an interface reflection wave phase axis between an asphalt concrete surface layer and a semi-rigid base layer in a ground penetrating radar section to obtain a maximum value A 12, and simultaneously obtaining an amplitude intensity change curve L 12 along the length of a measuring line;
Step 5), carrying out statistics and normalization treatment on amplitude values of radar data of each crack position in a time range of 1.7-2.7 ns below an interface reflection wave phase axis between a semi-rigid base layer and a filling road base in a ground penetrating radar profile, and obtaining a maximum value A 23 to obtain an amplitude intensity change curve L 23 along the length of a measuring line;
Step 6), find the peak points of the amplitude intensity A 1,A2,…Ai of the crack diffraction wave along the amplitude intensity change curve L 12, and find the peak points of the amplitude intensity of the crack diffraction wave along the amplitude intensity change curve L 23
Step 7), judging the crack as a fatigue crack or a reflection crack according to the size relation of A 12 and A 23 among different cracks;
And 8) judging the reflection crack as a vertical crack or an inclined crack according to the corresponding position relation between the actual position of the same crack and the amplitude intensity peak points in the amplitude intensity change curves L 12 and L 23.
Further, in the step 1), a ground penetrating radar is adopted to collect data along the direction perpendicular to the crack.
Further, in the step 1), the ground penetrating radar uses a ground coupling antenna of 400 MHz.
Further, in the step 2), the data editing includes the waste channel elimination, the direction adjustment, the drift processing and the channel equalization processing.
Further, the method for judging that the crack is a fatigue crack or a reflection crack according to the size relation of different cracks A 12 and A 23 in the step 7) is that all the detected cracks in the minimum unit of 1km form a sample library, the number of the cracks in the sample library is not less than 20, 2 cracks are randomly selected in the sample library and are respectively named as a crack 1 and a crack 2, when A 12 of the crack 1 is more than or equal to A 23 of the crack 2, the crack 1 is judged to be the fatigue crack, and when A 23 of the crack 1 is less than A 23 of the crack 2, the crack 2 is judged to be the reflection crack.
Further, in the step 8), when the peak values of the actual position and the amplitude intensity change curves L 12 and L 23 are on a vertical line, the crack is determined to be a vertical crack, otherwise, the crack is determined to be an oblique crack.
The invention has the beneficial effects that: the existing crack type and whether the inclination is judged through experience, quantification is not achieved, and accurate control of the whole process is difficult to achieve. The technique counts the amplitude value of radar reflected wave data at each crack position and performs normalization processing on the amplitude value, finds out the maximum values A 12 and A 23 near each crack position, and performs crack type identification by comparing the numerical relations of different cracks A 12 and A 23, and has the advantages that: ① The method fully utilizes the wave field characteristics of diffracted waves generated at the top of the crack, and completely eliminates the influence of artificial subjective factors and personal experience; ② The characteristics which can be compared are quantified, so that the probability of erroneous judgment caused by insufficient experience is reduced; ③ Whether the crack type is judged or whether the crack is inclined or not is judged, specific numerical values can be obtained through batch feature quantity statistics, the crack type can be judged and whether the crack is inclined or not can be judged by comparing the numerical values with the peak positions, and the efficiency is improved compared with the traditional experience judgment and naked eye identification. Compared with the prior art, the method for carrying out the crack type identification by calculating the absolute value of the amplitude of each path of radar data in a certain time range below the phase axis of the interface reflected wave, and carrying out the cumulative summation and average value obtaining, and the method solves the problems of low precision and low efficiency when the conventional empirical method is used for judging the types of a plurality of cracks of a long-distance section of a road.
The frequency of the common radar antenna is 200MHz, 400MHz, 900MHz and 1.5GHz, the frequency of the 200MHz antenna is too low, the resolution ratio cannot meet the precision requirement of crack detection, the detection depth of the 900MHz and 1.5GHz antenna is shallow, the detection requirement of deeper cracks at the bottom cannot be completely met, the resolution ratio of electromagnetic waves emitted by the 400MHz radar antenna can completely meet the resolution ratio requirement of conventional crack detection, meanwhile, the effective detection depth of the 400MHz antenna is deeper than the maximum depth of the bottoms of conventional cracks, the detection depth requirement can be met by selecting 400MHz, and the purpose of combining the two requirements can be achieved;
In the step 4), as the duration of the phase axis of the boundary reflection wave between the asphalt concrete surface layer (dielectric constant is 7.56) and the semi-rigid base layer (dielectric constant is 9) is within a specific range, the data within the range of 0.7-1.3 ns below the phase axis is selected to ensure that radar reflection waves at the interface between the asphalt concrete surface layer and the semi-rigid base layer can be completely contained within the time window;
In step 5), since the semi-rigid substrate (dielectric constant 9) is larger than the dielectric constant of the asphalt concrete surface layer (dielectric constant 7.56), the wavelength of the wavelet in the semi-rigid substrate is larger than that in the asphalt concrete surface layer, so that the frequency of the reflected wave is reduced and the period is increased, and the whole radar reflected wave can be included by selecting the time range of 1.7 ns-2.7 ns;
In the step 7), judging the type of the crack according to the size relation of A 12 and A 23 of different cracks, wherein the inventor carries out mathematical modeling and numerical simulation under the actual scene background to creatively obtain a rule summary, and obtains the verification of an experiment;
In step 8), the peak of the diffracted wave caused by the top of the crack is superimposed on the reflected wave generated by the interface between the asphalt surface layer and the water-stable layer, so that the amplitude value of the reflected wave at the upper layer is the largest at the diffraction point, and the transverse position of the diffraction point is correspondingly consistent with the transverse position of the crack peak; the bottom of the crack also generates diffraction waves, and the vertex of the generated diffraction waves is superposed on the reflected waves generated by the interface between the water-stable layer and the roadbed, so that the amplitude value of the lower reflected waves at the diffraction point is maximum, and the transverse position of the diffraction point corresponds to the transverse position of the bottom end point of the crack, so that when the actual position of the reflected crack and the peak value of the amplitude intensity change curves L 12 and L 23 are on the same plumb line, the reflected crack can be judged to be a vertical crack.
In a word, the asphalt pavement crack type rapid detection method based on the ground penetrating radar provided by the invention has the advantages of being scientific, accurate and rapid.
Drawings
FIG. 1 is a measured cross-sectional view of a road structural layer crack ground penetrating radar obtained in step 2) of a method for rapidly detecting a type of a crack of a bituminous pavement based on the ground penetrating radar, which is adopted in example 1;
FIG. 2 is an amplitude intensity variation curve L 12 obtained in step 4) of a method for rapidly detecting the type of a crack in a bituminous pavement based on a ground penetrating radar, which is adopted in example 1;
FIG. 3 is an amplitude intensity variation curve L 23 obtained in step 5) of a method for rapidly detecting a crack type of a bituminous pavement based on a ground penetrating radar, which is adopted in example 1;
FIG. 4 is a photograph of a sample of the drill core at the location of fracture 1 in example 1;
FIG. 5 is a photograph of a sample of the drill core at the location of fracture 4 in example 1;
FIG. 6 is a cross-sectional view of a pavement structure layer crack ground penetrating radar obtained in step 2) of a pavement crack type rapid detection method based on a ground penetrating radar, which is adopted in example 2;
FIG. 7 is a graph L 12 of the amplitude intensity variation curve obtained in the step 4) of the method for rapidly detecting the type of the crack of the asphalt pavement based on the ground penetrating radar, which is adopted in the embodiment 2;
FIG. 8 is a graph L 23 of the amplitude intensity variation curve obtained in step 5) of the method for rapidly detecting the type of the crack of the asphalt pavement based on the ground penetrating radar, which is adopted in the embodiment 2;
FIG. 9 is a photograph of a sample of the drill core at the location of fracture 1 in example 2;
FIG. 10 is a photograph of a sample of the drill core at the location of fracture 2 in example 2;
FIG. 11 is a cross-sectional view of a pavement structure layer crack ground penetrating radar obtained in step 2) of a pavement crack type rapid detection method based on a ground penetrating radar used in example 3;
FIG. 12 is a graph comparing the amplitude intensity curves L 12 and L 23 obtained in steps 4) and 5) of a method for rapidly detecting the type of a crack in a bituminous pavement based on a ground penetrating radar, which is adopted in example 3;
FIG. 13 is a photograph of a sample of the drill core at the location of fracture 1 in example 3;
FIG. 14 is a photograph of a sample of the drill core at the location of fracture 2 in example 3.
Detailed Description
The invention aims at providing a quick detection method for the type of asphalt pavement cracks based on a ground penetrating radar.
Example 1
The embodiment is characterized in that cracks of a structural layer of the Shaanxi Jing Wang Gaosu section are detected, the thickness of an asphalt surface layer of the highway is 15cm, and a base layer is a 36cm semi-rigid water-stable base layer. 2.6km of data are collected, the number of detected cracks is 56, a sample library is built, and 2 cracks are randomly selected from the sample library, namely, the cracks 1 and 4 in the embodiment. The contents of this embodiment are further described below with reference to the accompanying drawings:
A method for rapidly detecting crack type of asphalt pavement based on a ground penetrating radar is implemented according to the following steps:
Step 1), acquiring data along the driving direction of a road and perpendicular to the trend of a transverse crack by adopting a ground coupling antenna of a ground penetrating radar with a main frequency of 400MHz to obtain a waveform diagram of an electromagnetic wave time domain signal;
Step 2), data editing is carried out on the collected ground penetrating radar data, including waste channel elimination, direction adjustment, drift treatment, channel equalization, gain adjustment and the like, as shown in fig. 1, so as to obtain a road structural layer crack ground penetrating radar actual measurement sectional view; wherein, reject the waste channel: the antenna is not placed or invalid data acquisition is carried out in the antenna moving process, and the scanning is required to be proposed or zero-filled from radar data before the data processing, namely, the scanning is used for removing waste channels; drift treatment: the phenomenon that the data on the radar profile are all positive, all negative or asymmetric in positive and negative half cycles appears, the data contain direct current drift amount, and the direct current drift amount needs to be eliminated or suppressed before the data are processed so as to reduce the influence on the radar data;
Step 3), determining that the same phase axis at the depth of 4.62ns is a ground direct wave, the same phase axis of the reflected wave at the depth of 9.58ns is a boundary surface between the surface layer and the base layer, and the same phase axis of the reflected wave at the depth of 17.32ns is a boundary surface between the base layer and the roadbed according to the change of reflection amplitude of the electromagnetic wave among materials of each layer;
Step 4), according to the content shown in fig. 1, the actually measured section contains 4 cracks, in the ns- (9.58+1.3) ns time range below the interface reflection wave phase axis of the asphalt concrete surface layer and the semi-rigid base layer (9.58+0.7), the amplitude value of the upper reflection section in the range of diffraction waves (including two wings) caused by the cracks is counted and normalized, and then the maximum value A 12 is selected, so as to obtain an amplitude intensity change curve L 12, as shown in fig. 2;
Step 5), in the ns- (17.32+1.7) ns time range below the phase axis of the reflected wave at the interface between the semi-rigid base layer and the roadbed, counting the amplitude value of the lower layer reflected section in the range of the diffracted wave (including two wings) caused by the crack, normalizing the amplitude, and then selecting a selected maximum value A 23 to obtain an amplitude intensity change curve L 23, as shown in figure 3;
Step 6), find the peak points of the amplitude intensity A 1,A2,…Ai of the crack diffraction wave along the amplitude intensity change curve L 12, and find the peak points of the amplitude intensity of the crack diffraction wave along the amplitude intensity change curve L 23
Step 7), selecting a crack 1 and a crack 4 as comparison samples, comparing a change curve graph2 and a change curve graph 3 after normalization of amplitude values of reflected waves of an upper layer and a lower layer, wherein in the curve graph2, A 12 of the crack 1 is 0.7569, A 12 of the crack 4 is 0.7389, and crack 1A 12 is more than or equal to crack 4A 12, so that the crack 1 can be judged to be a fatigue crack; as shown in fig. 3, crack 1, a 23, is 0.7507, crack 4, a 23, is 1, crack 1a 23 < crack 4a 23, and crack 4 can be determined to be a reflective crack.
After the above conclusion is verified by the on-site core sampling method, the crack 1 is proved to be a fatigue crack (as shown in fig. 4), and the crack 4 is a reflection crack (as shown in fig. 5).
Example 2
In the embodiment, the Shanghai Shaanxi high-speed crack detection of a certain section of structural layer is performed, 1.8km of data is acquired, the number of detected cracks is 37, a sample library is built, and 2 cracks are randomly selected from the sample library, namely, the cracks 1 and the cracks 2 in the embodiment. The contents of this embodiment are further described below with reference to the accompanying drawings:
A method for rapidly detecting crack type of asphalt pavement based on a ground penetrating radar is implemented according to the following steps:
Step 1), acquiring data along the driving direction of a road and perpendicular to the trend of a transverse crack by adopting a ground coupling antenna of a ground penetrating radar with a main frequency of 400MHz to obtain a waveform diagram of an electromagnetic wave time domain signal;
Step 2), data editing is carried out on the collected ground penetrating radar data, including waste channel elimination, direction adjustment, drift treatment, channel equalization, gain adjustment and the like, and as shown in fig. 6, a road structural layer crack ground penetrating radar actual measurement sectional view is obtained; wherein, reject the waste channel: the antenna is not placed or invalid data acquisition is carried out in the antenna moving process, and the scanning is required to be proposed or zero-filled from radar data before the data processing, namely, the scanning is used for removing waste channels; drift treatment: the phenomenon that the data on the radar profile are all positive, all negative or asymmetric in positive and negative half cycles appears, the data contain direct current drift amount, and the direct current drift amount needs to be eliminated or suppressed before the data are processed so as to reduce the influence on the radar data;
Step 3), determining that the same phase axis at the depth of 5.01ns is a ground direct wave, the same phase axis of the reflected wave at the depth of 9.22ns is a boundary surface between the surface layer and the base layer, and the same phase axis of the reflected wave at the depth of 18.64ns is a boundary surface between the base layer and the roadbed according to the change of reflection amplitude of the electromagnetic wave among materials of each layer;
Step 4), according to the content shown in fig. 6, the actually measured section contains 2 cracks, in the time range of ns- (9.22+1.3) ns below the interface reflection wave phase of the asphalt concrete surface layer and the semi-rigid base layer (9.22+0.7), the amplitude value of the upper reflection section in the range of diffraction waves (including two wings) caused by the cracks is counted and normalized, and then the maximum value A 12 is selected, so as to obtain an amplitude intensity change curve L 12, as shown in fig. 7;
Step 5), in the ns- (18.64+1.7) ns time range below the phase axis of the reflected wave of the semi-rigid base layer and the roadbed interface, counting the amplitude value of the lower layer reflected section in the range of the diffracted wave (including two wings) caused by the crack, and normalizing the amplitude, and then selecting a selected maximum value A 23 to obtain an amplitude intensity change curve L 23, as shown in fig. 8;
Step 6), find the peak points of the amplitude intensity A 1,A2,…Ai of the crack diffraction wave along the amplitude intensity change curve L 12, and find the peak points of the amplitude intensity of the crack diffraction wave along the amplitude intensity change curve L 23
Step 7), selecting a crack 1 and a crack 2 as comparison samples, comparing a change curve graph 7 and a change curve graph 8 after normalization of amplitude values of reflected waves of an upper layer and a lower layer, wherein in the curve graph 7, A 12 of the crack 1 is 0.3892, A 12 of the crack 2 is 0.4084, and the crack 2A 12 is more than or equal to the crack 1A 12, so that the crack 2 can be judged to be a fatigue crack; as shown in fig. 8, crack 1, a 23, 0.9389, crack 2, a 23, 0.4766, and crack 2a 23 < crack 1a 23, can be determined to be reflective crack 1.
After the above conclusion is verified by the on-site core sampling method, the crack 1 is proved to be a fatigue crack (as shown in fig. 9), and the crack 2 is proved to be a reflection crack (as shown in fig. 10).
Example 3
In the embodiment, the cracks of a structural layer of the Wuhan winding city at a high speed are detected, 3.2km of data are collected, the number of detected cracks is 63, a sample library is built, and 2 cracks are randomly selected from the sample library, namely, the cracks 1 and the cracks 2 in the embodiment. The contents of this embodiment are further described below with reference to the accompanying drawings:
A method for rapidly detecting crack type of asphalt pavement based on a ground penetrating radar is implemented according to the following steps:
Step 1), acquiring data along the driving direction of a road and perpendicular to the trend of a transverse crack by adopting a ground coupling antenna of a ground penetrating radar with a main frequency of 400MHz to obtain a waveform diagram of an electromagnetic wave time domain signal;
Step 2), data editing is carried out on the collected ground penetrating radar data, including waste channel elimination, direction adjustment, drift treatment, channel equalization, gain adjustment and the like, and as shown in fig. 11, a road structural layer crack ground penetrating radar actual measurement sectional view is obtained; wherein, reject the waste channel: the antenna is not placed or invalid data acquisition is carried out in the antenna moving process, and the scanning is required to be proposed or zero-filled from radar data before the data processing, namely, the scanning is used for removing waste channels; drift treatment: the phenomenon that the data on the radar profile are all positive, all negative or asymmetric in positive and negative half cycles appears, the data contain direct current drift amount, and the direct current drift amount needs to be eliminated or suppressed before the data are processed so as to reduce the influence on the radar data;
Step 3), determining that the same phase axis at the depth of 4.85ns is a ground direct wave, the same phase axis of the reflected wave at the depth of 9.86ns is a boundary surface between the surface layer and the base layer, and the same phase axis of the reflected wave at the depth of 17.96ns is a boundary surface between the base layer and the roadbed according to the change of reflection amplitude of the electromagnetic wave among materials of each layer;
Step 4), according to the content shown in fig. 11, the actually measured section contains 2 cracks, in the ns- (9.86+1.3) ns time range below the interface reflection wave phase axis of the asphalt concrete surface layer and the semi-rigid base layer (9.86+0.7), the amplitude value of the upper reflection section in the diffraction wave (including two wings) range caused by the cracks is counted and normalized, and then the maximum value A 12 is selected, so as to obtain an amplitude intensity change curve L 12, as shown by a line a) in fig. 12;
Step 5), in the ns- (17.96+1.7) ns time range below the phase axis of the reflected wave of the semi-rigid base layer and the roadbed interface, counting the amplitude value of the lower layer reflected section in the range of the diffracted wave (including two wings) caused by the crack, and normalizing the amplitude, and then selecting a selected maximum value A 23 to obtain an amplitude intensity change curve L 23, as shown by a line b) in fig. 12;
Step 6), find the peak points of the amplitude intensity A 1,A2,…Ai of the crack diffraction wave along the amplitude intensity change curve L 12, and find the peak points of the amplitude intensity of the crack diffraction wave along the amplitude intensity change curve L 23
Step 7), judging whether the crack is an inclined crack or a vertical crack according to the position relation of the top and bottom strongest amplitude values of the crack as inclination or plumb, wherein the amplitude intensity change curves L 12 and L 23 of the crack 1 and the crack 2 in fig. 11 are respectively shown as the curves a) and b) in fig. 12, connecting the peak points of L 12 and L 23 of the crack 1 to obtain a left inclined arrow, and judging that the crack 1 is the crack inclined to the left; and connecting the peak points of the L 12 and the L 23 of the crack 2 to obtain an arrow in the plumb direction, and judging that the crack 2 is a vertical crack.
After the above conclusion is verified by the on-site core sampling method, the crack 1 is proved to be an inclined crack (as shown in fig. 13), and the crack 4 is proved to be a vertical crack (as shown in fig. 14).
Claims (6)
1. The quick detection method for the asphalt pavement crack type based on the ground penetrating radar is characterized by comprising the following steps of:
Step 1), acquiring underground data of a crack by adopting a ground penetrating radar to obtain a waveform diagram of an electromagnetic wave time domain signal;
step 2), editing the collected ground penetrating radar data;
Step 3), determining interlayer interfaces between the asphalt surface layer and the base layer and between the base layer and the roadbed according to the change of reflection amplitude of electromagnetic waves among materials of each layer;
step 4), carrying out statistics and normalization treatment on amplitude values of radar data of each crack position in a time range of 0.7-1.3 ns below an interface reflection wave phase axis between an asphalt concrete surface layer and a semi-rigid base layer in a ground penetrating radar section, and obtaining a maximum value A 12 to obtain an amplitude intensity change curve L 12 along the length of a measuring line;
Step 5), carrying out statistics and normalization treatment on amplitude values of radar data of each crack position in a time range of 1.7-2.7 ns below an interface reflection wave phase axis between a semi-rigid base layer and a filling road base in a ground penetrating radar profile, and obtaining a maximum value A 23 to obtain an amplitude intensity change curve L 23 along the length of a measuring line;
Step 6), find the peak points of the amplitude intensity A 1,A2,…Ai of the crack diffraction wave along the amplitude intensity change curve L 12, and find the peak points of the amplitude intensity of the crack diffraction wave along the amplitude intensity change curve L 23
Step 7), judging the crack as a fatigue crack or a reflection crack according to the size relation of A 12 and A 23 among different cracks;
And 8) judging the crack to be a vertical crack or an inclined crack according to the corresponding position relation between the actual position of the same crack and the amplitude intensity peak points in the amplitude intensity change curves L 12 and L 23.
2. The method for rapidly detecting the type of the asphalt pavement crack based on the ground penetrating radar according to claim 1, which is characterized by comprising the following steps: in the step 1), a ground penetrating radar is adopted to collect data along the direction perpendicular to the crack.
3. The method for rapidly detecting the type of the asphalt pavement crack based on the ground penetrating radar according to claim 1, which is characterized by comprising the following steps: in the step 1), the ground penetrating radar uses a ground coupling antenna of 400 MHz.
4. The method for rapidly detecting the type of the asphalt pavement crack based on the ground penetrating radar according to claim 1, which is characterized by comprising the following steps: in the step 2), the data editing comprises waste channel elimination, direction adjustment, drift treatment and channel equalization treatment.
5. The method for rapidly detecting the type of the asphalt pavement crack based on the ground penetrating radar according to claim 1, which is characterized by comprising the following steps: the step 7) is to judge that the crack is a fatigue crack or a reflection crack according to the size relation of different cracks A 12 and A 23, wherein the method comprises the steps that a sample library is formed by all the detected cracks in a minimum unit of 1km, the number of the cracks in the sample library is not less than 20, 2 cracks are randomly selected in the sample library and are respectively named as a crack 1 and a crack 2, when A 12 of the crack 1 is more than or equal to A 23 of the crack 2, the crack 1 is judged to be the fatigue crack, and when A 23 of the crack 1 is less than A 23 of the crack 2, the crack 2 is judged to be the reflection crack.
6. The method for rapidly detecting the type of the asphalt pavement crack based on the ground penetrating radar according to claim 1, which is characterized by comprising the following steps: in the step 8), when the peak values of the actual position and the amplitude intensity change curves L 12 and L 23 are on a vertical line, the crack is judged to be a vertical crack, otherwise, the crack is judged to be an inclined crack.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210146001.1A CN114624704B (en) | 2022-02-17 | 2022-02-17 | Asphalt pavement crack type rapid detection method based on ground penetrating radar |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210146001.1A CN114624704B (en) | 2022-02-17 | 2022-02-17 | Asphalt pavement crack type rapid detection method based on ground penetrating radar |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN114624704A CN114624704A (en) | 2022-06-14 |
| CN114624704B true CN114624704B (en) | 2024-04-19 |
Family
ID=81900725
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202210146001.1A Active CN114624704B (en) | 2022-02-17 | 2022-02-17 | Asphalt pavement crack type rapid detection method based on ground penetrating radar |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN114624704B (en) |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109440610A (en) * | 2018-12-18 | 2019-03-08 | 江苏中路工程技术研究院有限公司 | A kind of pavement crack width, depth determination method and quick maintenance method based on Ground Penetrating Radar |
| CN112666554A (en) * | 2020-12-17 | 2021-04-16 | 江苏中路工程技术研究院有限公司 | Method for identifying radar amplitude characteristic crack width of asphalt pavement |
| CN112684440A (en) * | 2020-12-17 | 2021-04-20 | 江苏中路工程技术研究院有限公司 | Crack evaluation method based on three-dimensional ground penetrating radar detection |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US10392929B1 (en) * | 2018-02-08 | 2019-08-27 | Saudi Arabian Oil Company | Mapping fracture length using downhole ground penetrating radar |
-
2022
- 2022-02-17 CN CN202210146001.1A patent/CN114624704B/en active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109440610A (en) * | 2018-12-18 | 2019-03-08 | 江苏中路工程技术研究院有限公司 | A kind of pavement crack width, depth determination method and quick maintenance method based on Ground Penetrating Radar |
| CN112666554A (en) * | 2020-12-17 | 2021-04-16 | 江苏中路工程技术研究院有限公司 | Method for identifying radar amplitude characteristic crack width of asphalt pavement |
| CN112684440A (en) * | 2020-12-17 | 2021-04-20 | 江苏中路工程技术研究院有限公司 | Crack evaluation method based on three-dimensional ground penetrating radar detection |
Non-Patent Citations (1)
| Title |
|---|
| 裂缝宽度探地雷达波场响应的数值研究;郭士礼;闫飞;朱培民;李修忠;;地球物理学进展;20160815(04);全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN114624704A (en) | 2022-06-14 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108981557B (en) | Detection method for simultaneously measuring diameter of steel bar in concrete and thickness of protective layer of steel bar | |
| CN110346454B (en) | Concrete shallow surface layer ultrasonic surface wave detection method based on arrayed ultrasonic seismic sources | |
| CN105716969A (en) | Asphalt pavement roughening effectiveness evaluation method based on three-dimensional laser detection technology | |
| CN112666554A (en) | Method for identifying radar amplitude characteristic crack width of asphalt pavement | |
| CN106443668B (en) | The evaluation method of highway pavement splitting based on Ground Penetrating Radar | |
| CN114624704B (en) | Asphalt pavement crack type rapid detection method based on ground penetrating radar | |
| Hopman et al. | Repeatability, reproducibility and accuracy of GPR measurements | |
| CN115597530A (en) | Highway medium parameter inversion method | |
| CN112684440A (en) | Crack evaluation method based on three-dimensional ground penetrating radar detection | |
| CN104100262A (en) | Method for identifying properties of reservoir fluid by three-dimensional quantitative fluorescence parameter charts | |
| CN108333096A (en) | A kind of asphalt concrete pavement porosity detection method based on Ground Penetrating Radar | |
| CN112835124B (en) | Crack effectiveness evaluation method based on imaging logging and array acoustic logging data | |
| CN111398687B (en) | Test method for estimating dielectric constant of asphalt pavement | |
| CN106337354B (en) | The evaluation method of highway pavement isolation based on Ground Penetrating Radar | |
| Jin-sung et al. | Technique for Detecting Subsurface Cavities of Urban Road Using Multichannel Ground-penetrating Radar Equipment. | |
| CN110715963A (en) | Method for detecting ice content of concrete | |
| Huang et al. | Research on void signal recognition algorithm of 3D ground-penetrating radar based on the digital image | |
| CN115825411A (en) | Crack size evaluation method | |
| CN115876661A (en) | Asphalt pavement void ratio evaluation method based on three-dimensional ground penetrating radar image | |
| CN112014881B (en) | Water flooding speed prediction method based on time lapse earthquake | |
| Chen et al. | Automatic detection of asphalt layer thickness based on Ground Penetrating Radar | |
| Wang et al. | Designing a skid-resistant and durable asphalt mixture based on the stress concentration distribution rate | |
| CN110082758B (en) | Pavement interlayer void detection and maintenance method | |
| CN114397656A (en) | Dynamic asphalt pavement thickness adjusting method based on three-dimensional ground penetrating radar | |
| Batrakov et al. | GPR signals processing for the road pavements monitoring with the weak contrast layers |
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 |