CN114624704A - Method for rapidly detecting crack type of asphalt pavement based on ground penetrating radar - Google Patents
Method for rapidly detecting crack type of asphalt pavement based on ground penetrating radar Download PDFInfo
- Publication number
- CN114624704A CN114624704A CN202210146001.1A CN202210146001A CN114624704A CN 114624704 A CN114624704 A CN 114624704A CN 202210146001 A CN202210146001 A CN 202210146001A CN 114624704 A CN114624704 A CN 114624704A
- Authority
- CN
- China
- Prior art keywords
- crack
- ground penetrating
- penetrating radar
- cracks
- amplitude
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 230000000149 penetrating effect Effects 0.000 title claims abstract description 53
- 238000000034 method Methods 0.000 title claims abstract description 44
- 239000010426 asphalt Substances 0.000 title claims abstract description 32
- 238000012545 processing Methods 0.000 claims abstract description 22
- 238000010606 normalization Methods 0.000 claims abstract description 10
- 239000010410 layer Substances 0.000 claims description 54
- 230000008859 change Effects 0.000 claims description 19
- 239000002344 surface layer Substances 0.000 claims description 17
- 239000011384 asphalt concrete Substances 0.000 claims description 9
- 239000002699 waste material Substances 0.000 claims description 8
- 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
- 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
- 239000002689 soil Substances 0.000 claims description 3
- 239000011229 interlayer Substances 0.000 claims description 2
- 238000001514 detection method Methods 0.000 abstract description 15
- 238000004364 calculation method Methods 0.000 abstract description 5
- 238000004836 empirical method Methods 0.000 abstract description 3
- 238000005553 drilling Methods 0.000 description 4
- 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
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 3
- 238000013461 design Methods 0.000 description 2
- 238000011065 in-situ storage Methods 0.000 description 2
- 238000012935 Averaging Methods 0.000 description 1
- 201000004569 Blindness Diseases 0.000 description 1
- 238000009825 accumulation Methods 0.000 description 1
- 230000009286 beneficial 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
- 238000012423 maintenance Methods 0.000 description 1
- 238000011002 quantification Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 230000000087 stabilizing effect Effects 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
- 239000011800 void material Substances 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/88—Radar or analogous systems specially adapted for specific applications
- G01S13/885—Radar or analogous systems specially adapted for specific applications for ground probing
-
- 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 method for rapidly detecting asphalt pavement crack types based on a ground penetrating radar, which can reduce the complexity of theoretical calculation and has high judgment accuracy. The technology carries out crack type identification by calculating the amplitude value of each path of radar data within a certain time range below the isophase axis of interface reflection waves for statistics and carrying out normalization processing, and solves the problem of low precision and low efficiency when judging multiple crack types of long-distance sections of roads by an existing empirical method.
Description
Technical Field
The invention belongs to the technical field of pavement crack detection, and particularly relates to a method for rapidly detecting asphalt pavement crack types based on a ground penetrating radar, which can reduce the complexity of theoretical calculation and has high judgment accuracy.
Background
The asphalt pavement crack types mainly comprise fatigue cracks and reflection cracks, the fatigue cracks generally cause the fatigue cracks of an asphalt surface layer under the action of repeated loading of vehicles and only exist in the asphalt surface layer, and the reflection cracks are formed by reflecting transverse cracks generated by a semi-rigid base layer under the action of temperature and humidity change to the pavement. When external free water continuously enters the structural layer through a pavement crack and is stored at the top of the roadbed and the base layer, the fine aggregate of the base layer and the roadbed soil are flushed and extruded out under the repeated action of the driving load, so that the bottom between layers loses support to form local void, and the structural damage of the road is caused. The traditional detection method judges the type of the crack by core drilling sampling, cannot accurately know the type and the internal development condition of each crack, cannot carry out targeted treatment according to the type of each crack during maintenance design, and enables the design to have certain blindness.
In recent years, ground penetrating radar is widely applied to the field of fine nondestructive detection of deep diseases of roads by virtue of the advantages of continuity, no damage, high speed, high precision and the like. The ground penetrating radar is in a single uniform ideal medium, and the top and bottom points of the crack correspond to a diffracted wave hyperbolic wave group, so that the top and bottom points are easy to identify. However, when the ground penetrating radar propagates in the layered uniform medium, a plurality of diffracted wave hyperbolic wave groups are formed on each medium interface which passes through the ground penetrating radar, so that the diffracted wave hyperbolic 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 judgment of the type of the reflection crack mainly depends 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 complex medium of the road structure interferes with the definition of a hyperbola so as to influence the judgment of the experience; secondly, the detection of the theoretical calculation complexity on the long section of the road is difficult to quickly and accurately detect and identify the type of the crack.
Disclosure of Invention
The invention aims to overcome the defects of the prior art and provides a method for quickly detecting the crack type of the asphalt pavement based on a ground penetrating radar, which is used for solving the problems of low accuracy and complexity of theoretical calculation of the conventional empirical judgment method.
The technical scheme adopted by the invention is as follows:
a method for rapidly detecting the crack type of an asphalt pavement based on a ground penetrating radar comprises the following steps:
step 1), collecting crack underground data by adopting a ground penetrating radar to obtain a oscillogram 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 the reflection amplitude of the electromagnetic waves between the materials of each layer;
step 4), counting the amplitude values of radar data of all crack positions within the time range of 0.7-1.3 ns below the homophase axis of the interface reflection wave between the asphalt concrete surface layer and the semi-rigid base layer in the ground penetrating radar section, and carrying out normalization processing to obtain the maximum value A12Simultaneously obtaining the amplitude intensity variation curve L of the measured line length12；
Step 5), counting the amplitude values of radar data of all crack positions within the time range of 1.7-2.7 ns below the homophase axis of the interface reflection wave between the semi-rigid base and the soil filling subgrade in the ground penetrating radar section, and carrying out normalization processing to obtain the maximum value A23Obtaining the amplitude intensity variation curve L along the length of the measuring line23；
Step 6), along the amplitude intensity variation curve L12Respectively finding the peak value A of the amplitude intensity of the diffraction wave of the crack1,A2,…AiVariation in amplitude intensityCurve L23Respectively finding the peak value point of the amplitude intensity of the diffraction wave of the crack
Step 7), according to A between different cracks12And A23Judging the crack to be a fatigue crack or a reflection crack;
step 8), according to the actual position and the amplitude intensity change curve L of the same crack12And L23And judging the reflection crack to be a vertical crack or an inclined crack according to the corresponding position relation of the amplitude intensity peak point.
Further, in the step 1), a ground penetrating radar is adopted to collect data along the direction perpendicular to the fracture.
Further, in the step 1), the ground penetrating radar uses a 400MHz ground coupling antenna.
Further, in the step 2), the data editing includes a waste track elimination, a direction adjustment, a drift processing, and a track balancing processing.
Further, the step 7) is carried out according to different cracks A12And A23The method for judging the crack as the fatigue crack or the reflection crack comprises the steps of forming a sample library by all cracks detected in the minimum unit of 1km, wherein the number of the cracks in the sample library is not less than 20, randomly selecting 2 cracks in the sample library, respectively naming the cracks as 1 and 2, and when A of the crack 1 is smaller than A, judging the crack as the fatigue crack or the reflection crack12A is not less than crack 223When the crack 1 is a fatigue crack, A of the crack 1 is judged23< A of crack 223Then, the crack 2 is judged to be a reflection crack.
Further, in the step 8), when the actual position and the amplitude intensity variation curve L are in the same12And L23If the peak value of (a) is on a vertical line, the crack is judged to be a vertical crack, otherwise, the crack is judged to be an inclined crack.
The invention has the beneficial effects that: the existing crack type and whether the crack is inclined are judged by experience, quantification is not achieved, and accurate control of the whole process is difficult to achieve. The skill of the artCalculating amplitude value of radar reflected wave data of each crack position and normalizing the amplitude value to find out maximum value A near each crack position12And A23By comparing different cracks A12And A23The numerical relation method is used for identifying the crack type, and has the advantages that: firstly, the method fully utilizes the wave field characteristics of the diffracted waves generated at the top end of the crack, and completely eliminates the influence of artificial subjective factors and personal experience; quantifying the characteristics capable of being compared, and reducing the probability of misjudgment caused by insufficient experience; no matter the judgment of the crack type or the judgment of whether the crack inclines or not, a specific numerical value can be obtained through the statistics of the characteristic quantity in batches, the crack type can be judged and whether the crack inclines or not can be judged by comparing the numerical value with the peak value, and the efficiency is improved by comparing with the traditional experience judgment and the naked eye identification. Compared with the prior art, the method for identifying the crack type by calculating the absolute value of the amplitude of each path of radar data within a certain time range below the isophase axis of the interface reflection wave carries out accumulation summation and averaging, and solves the problem of low precision and low efficiency when the existing empirical method is used for judging the types of a plurality of cracks on the long-distance section of the road.
The frequency of a common radar antenna comprises 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 antennas is shallow, the detection requirement of a deep crack 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 the conventional crack detection, meanwhile, the effective detection depth of the 400MHz antenna is deeper than the maximum depth of the bottom of the conventional crack, the 400MHz antenna is selected to meet both the detection depth requirement and the resolution ratio requirement, and the aim of taking both the detection depth and the resolution ratio can be fulfilled;
in the step 4), because the duration of the homophase axis of the boundary reflected wave between the asphalt concrete surface layer (the dielectric constant is 7.56) and the semi-rigid base layer (the dielectric constant is 9) is within a specific range, selecting data within the range of 0.7-1.3 ns below the homophase axis can ensure that all radar reflected waves of the interface of the asphalt concrete surface layer and the semi-rigid base layer can be included in the time window;
in the step 5), because the semi-rigid base layer (with the dielectric constant of 9) is larger than the dielectric constant of the asphalt concrete surface layer (with the dielectric constant of 7.56), the wavelength of the wavelet in the semi-rigid base layer is larger than that of the asphalt concrete surface layer, so that the frequency of the reflected wave is reduced and the period is increased, and therefore, the time range of 1.7 ns-2.7 ns is selected to include all the radar reflected waves;
in step 7), based on comparing A of different fractures12And A23The method is characterized in that the type of the crack is judged according to the size relationship, and is characterized in that the inventor summarizes rules creatively obtained by mathematical modeling and numerical simulation under the actual scene background, and experimental verification is obtained;
in step 8), the vertex of the diffraction wave caused by the top of the crack is superposed on the reflected wave generated by the interface of the asphalt surface layer and the water-stable layer, so that the amplitude value of the upper layer reflected wave at the diffraction point is maximum, and the transverse position of the diffraction point is correspondingly consistent with the transverse position of the vertex of the crack; the bottom of the crack can also generate diffracted waves, and the vertex of the generated diffracted waves can be superposed on the reflected waves generated by the interface of the water stabilizing layer and the roadbed to cause that the amplitude value of the lower reflected wave at the diffraction point is maximum, and the transverse position of the diffraction point is correspondingly consistent with the transverse position of the bottom end point of the crack, so that when the actual position of the reflection crack and the amplitude intensity change curve L are reflected, the actual position of the crack and the amplitude intensity change curve L are consistent12And L23The peak value of (a) is on a vertical line, and the reflection crack can be judged to be a vertical crack.
In a word, the method for rapidly detecting the type of the asphalt pavement crack based on the ground penetrating radar has the advantages of being scientific, accurate and rapid.
Drawings
FIG. 1 is a cross-sectional view of a road structure layer crack ground penetrating radar actually measured in step 2) of a method for rapidly detecting the type of the asphalt pavement crack based on the ground penetrating radar adopted in embodiment 1;
FIG. 2 is an amplitude intensity variation curve L obtained in step 4) of the method for rapidly detecting the type of the asphalt pavement cracks based on the ground penetrating radar adopted in embodiment 112；
FIG. 3 is a block diagramAmplitude intensity change curve L obtained in step 5) of method for rapidly detecting type of asphalt pavement crack based on ground penetrating radar adopted in example 123；
FIG. 4 is a photograph of a core sample of the location of fracture 1 in example 1;
FIG. 5 is a photograph of a core sample of the location of the fracture 4 in example 1;
FIG. 6 is a cross-sectional view of a road structure layer crack ground penetrating radar actually measured in step 2) of the method for rapidly detecting the type of the asphalt pavement crack based on the ground penetrating radar adopted in the embodiment 2;
FIG. 7 is an amplitude intensity variation curve L obtained in step 4) of the method for rapidly detecting the type of the asphalt pavement cracks based on the ground penetrating radar adopted in example 212；
FIG. 8 is an amplitude intensity variation curve L obtained in step 5) of the method for rapidly detecting the type of the asphalt pavement cracks based on the ground penetrating radar adopted in embodiment 223；
FIG. 9 is a photograph of a core sample of the location of fracture 1 in example 2;
FIG. 10 is a photograph of a core sample of the location of the fracture 2 in example 2;
FIG. 11 is a cross-sectional view of a road structure layer crack ground penetrating radar actually measured in step 2) of the method for rapidly detecting the type of the asphalt pavement crack based on the ground penetrating radar adopted in the embodiment 3;
FIG. 12 is an amplitude intensity variation curve L obtained in steps 4) and 5) of the method for rapidly detecting the type of the asphalt pavement cracks based on the ground penetrating radar adopted in embodiment 312And L23Comparing the images;
FIG. 13 is a photograph of a core sample of the location of fracture 1 in example 3;
FIG. 14 is a photograph of a core sample of the location of the fracture 2 in example 3.
Detailed Description
The core of the invention is to provide a method for rapidly detecting the crack type of the asphalt pavement based on the ground penetrating radar.
Example 1
The embodiment is the crack detection of a certain section of high-speed structural layer of shanxi Jing Wang, and this highway asphalt surface course thickness 15cm, basic unit are 36 cm's the stable basic unit of semi-rigid water. 2.6km of data is collected, 56 cracks are detected, a sample library is established, and 2 cracks are randomly selected from the sample library, namely the crack 1 and the crack 4 in the embodiment. The contents of this embodiment will be further described with reference to the accompanying drawings:
a method for rapidly detecting the crack type of an asphalt pavement based on a ground penetrating radar is implemented according to the following steps:
step 1), adopting a ground penetrating radar ground coupling antenna with a main frequency of 400MHz, and collecting data along the driving direction of a road and perpendicular to the trend of a transverse crack to obtain a oscillogram of an electromagnetic wave time domain signal;
step 2), carrying out data editing on the collected ground penetrating radar data, wherein the data editing comprises waste channel elimination, direction adjustment, drifting treatment, channel equalization, gain adjustment and the like, and obtaining a road structure layer crack ground penetrating radar actual measurement section view as shown in figure 1; wherein the waste channel is removed: the antenna is not placed well or invalid data acquisition is carried out in the antenna moving process, and the scanning is carried out or zero filling is carried out in the radar data before data processing, namely, the scanning is rejected; and (3) drift processing: the data on the radar profile is all positive, all negative or asymmetric in positive and negative half cycles, the data contains a direct current drift amount, and the direct current drift amount needs to be eliminated or suppressed before data processing so as to reduce the influence of the direct current drift amount on the radar data;
step 3), according to the change of the reflection amplitude of the electromagnetic wave among the materials of each layer, 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 surface layer and base layer interface, and the same phase axis of the reflected wave at the depth of 17.32ns is a base layer and roadbed interface;
step 4), according to the content shown in the graph 1, the actually measured section comprises 4 cracks, the amplitude value of the upper layer reflection section in the range of the diffracted wave (including two wings) caused by the cracks is subjected to statistical amplitude value within the time range of (9.58+0.7) ns- (9.58+1.3) ns below the reflection wave coaxial axis of the interface of the asphalt concrete surface layer and the semi-rigid base layer, and the maximum value A is selected after the statistical amplitude value is subjected to normalization processing12To obtainTo amplitude intensity curve L12As shown in fig. 2;
step 5), in the time range of (17.32+1.7) ns- (17.32+2.7) ns below the reflection wave homophase axis of the interface of the semi-rigid base layer and the roadbed, counting the amplitude value of the lower layer reflection section in the range of diffraction waves (including two wings) caused by cracks, normalizing the amplitude, and selecting a selected maximum value A23Obtaining the amplitude intensity variation curve L23As shown in fig. 3;
step 6), along the amplitude intensity variation curve L12Respectively finding the peak value A of the amplitude intensity of the diffraction wave of the crack1,A2,…AiAlong the amplitude intensity curve L23Respectively finding the peak value points of the amplitude intensity of the diffraction wave of the crack
Step 7), selecting the crack 1 and the crack 4 as comparison samples, and comparing a change curve chart 2 and a change curve chart 3 after normalization of the amplitude values of the reflected waves of the upper layer and the lower layer, wherein the change curve chart is shown in figure 2, and A of the crack 1 is shown in figure 2120.7569, A of crack 4120.7389, crack 1A12Not less than 4A of crack12The crack 1 can be judged to be a fatigue crack; shown in FIG. 3, A of crack 1230.7507, A of crack 423Is 1, a crack 1A23< crack 4A23The crack 4 can be determined to be a reflection crack.
The above conclusion is verified by the in-situ core-drilling sampling method, and the crack 1 is confirmed to be a fatigue crack (as shown in fig. 4), and the crack 4 is confirmed to be a reflection crack (as shown in fig. 5).
Example 2
In this embodiment, for crack detection of a certain structural layer at high speed in Shanxi province, 1.8km of data is collected, 37 detected cracks are used to establish a sample library, and 2 cracks are randomly selected from the sample library, i.e., the crack 1 and the crack 2 in this embodiment. The contents of this embodiment will be further described with reference to the accompanying drawings:
a method for rapidly detecting the crack type of an asphalt pavement based on a ground penetrating radar is implemented according to the following steps:
step 1), adopting a ground penetrating radar ground coupling antenna with a main frequency of 400MHz, and collecting data along the driving direction of a road and perpendicular to the trend of a transverse crack to obtain a oscillogram of an electromagnetic wave time domain signal;
step 2), carrying out data editing on the collected ground penetrating radar data, wherein the data editing comprises waste channel elimination, direction adjustment, drifting treatment, channel equalization, gain adjustment and the like, and obtaining a road structure layer crack ground penetrating radar actual measurement section view as shown in fig. 6; wherein the waste channel is removed: the antenna is not placed or invalid data acquisition is carried out in the antenna moving process, and the scanning is carried out or filled with zero in the radar data before data processing, namely the scanning is rejected; and (3) drift processing: the data on the radar profile is all positive, all negative or asymmetric in positive and negative half cycles, the data contains a direct current drift amount, and the direct current drift amount needs to be eliminated or suppressed before data processing so as to reduce the influence of the direct current drift amount on the radar data;
step 3), according to the change of the reflection amplitude of the electromagnetic wave among the materials of each layer, 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 surface layer and base layer interface, and the same phase axis of the reflected wave at the depth of 18.64ns is a base layer and roadbed interface;
step 4), according to the content shown in the graph 6, the actually measured section contains 2 cracks, the statistical amplitude value of the upper layer reflection section amplitude value in the range of the diffracted wave (including two wings) caused by the cracks is carried out in the time range of (9.22+0.7) ns- (9.22+1.3) ns below the reflection wave homophase axis of the interface of the asphalt concrete surface layer and the semi-rigid base layer, and the maximum value A is selected after the statistical amplitude value is subjected to normalization processing12Obtaining the amplitude intensity variation curve L12As shown in fig. 7;
step 5), in the time range of (18.64+1.7) ns- (18.64+2.7) ns below the reflection wave homophase axis of the interface of the semi-rigid base layer and the roadbed, counting the amplitude value of the lower layer reflection section in the range of diffraction waves (including two wings) caused by cracks, normalizing the amplitude, and selecting a selected maximum value A23Obtaining an amplitude intensity variation curve L23As shown in fig. 8;
step 6), along the amplitude intensity variation curve L12Respectively finding the peak value A of the amplitude intensity of the diffraction wave of the crack1,A2,…AiAlong the amplitude intensity curve L23Respectively finding the peak value points of the amplitude intensity of the diffraction wave of the crack
Step 7), selecting the crack 1 and the crack 2 as comparison samples, and comparing a change curve chart 7 after normalization of the amplitude values of the reflected waves of the upper layer and the lower layer with a change curve chart 8, wherein the change curve chart is shown in FIG. 7, and A of the crack 1 is shown in FIG. 7120.3892, A of crack 2120.4084, crack 2A12Not less than 1A of crack12The crack 2 can be judged to be a fatigue crack; shown in FIG. 8, A of crack 1230.9389, A of crack 2230.4766, crack 2A23< crack 1A23The crack 1 can be determined to be a reflection crack.
The above conclusion is verified by the on-site core drilling sampling method, and the fracture 1 is confirmed to be a fatigue fracture (as shown in fig. 9), and the fracture 2 is confirmed to be a reflection fracture (as shown in fig. 10).
Example 3
In this embodiment, fracture detection is performed on a structural layer of a certain high-speed section in wuhan city around, 3.2km data is collected, 63 fractures are detected, a sample library is established, and 2 fractures are randomly selected from the sample library, namely the fracture 1 and the fracture 2 in this embodiment. The contents of this embodiment will be further described with reference to the accompanying drawings:
a method for rapidly detecting the type of a bituminous pavement crack based on a ground penetrating radar is implemented according to the following steps:
step 1), adopting a ground penetrating radar ground coupling antenna with a main frequency of 400MHz, and collecting data along the driving direction of a road and perpendicular to the trend of a transverse crack to obtain a oscillogram of an electromagnetic wave time domain signal;
step 2), carrying out data editing on the collected ground penetrating radar data, wherein the data editing comprises waste channel elimination, direction adjustment, drifting processing, channel equalization, gain adjustment and the like, and obtaining a road structure layer crack ground penetrating radar actual measurement section view as shown in fig. 11; wherein the waste channel is removed: the antenna is not placed or invalid data acquisition is carried out in the antenna moving process, and the scanning is carried out or filled with zero in the radar data before data processing, namely the scanning is rejected; and (3) drift processing: the data on the radar profile is all positive, all negative or asymmetric in positive and negative half cycles, the data contains a direct current drift amount, and the direct current drift amount needs to be eliminated or suppressed before data processing so as to reduce the influence of the direct current drift amount on the radar data;
step 3), according to the change of the reflection amplitude of the electromagnetic wave among the materials of each layer, 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 surface layer and base layer interface, and the same phase axis of the reflected wave at the depth of 17.96ns is a base layer and roadbed interface;
step 4), according to the content shown in the graph 11, the actually measured section comprises 2 cracks, the amplitude value of the upper layer reflection section in the range of the diffraction wave (including two wings) caused by the cracks is subjected to statistical amplitude value and normalized within the time range of (9.86+0.7) ns- (9.86+1.3) ns below the reflection wave event of the interface of the asphalt concrete surface layer and the semi-rigid base layer, and then the maximum value A is selected12Obtaining the amplitude intensity variation curve L12As indicated by line a) in fig. 12;
step 5), in the time range of (17.96+1.7) ns- (17.96+2.7) ns below the same phase axis of the reflection wave of the interface of the semi-rigid base layer and the roadbed, counting the amplitude value of the lower layer reflection section in the range of diffraction waves (including two wings) caused by cracks, and selecting a selected maximum value A after normalization processing of the amplitude value23Obtaining an amplitude intensity variation curve L23As indicated by line b) in fig. 12;
step 6), along the amplitude intensity variation curve L12Respectively finding the peak value A of the amplitude intensity of the diffraction wave of the crack1,A2,…AiAlong the amplitude intensity curve L23Respectively finding the peak value points of the amplitude intensity of the diffraction wave of the crack
Step 7), judging whether the crack is an inclined crack or a vertical crack according to the inclination or plumb position relation of the strongest amplitude values of the top and the bottom of the crack, and judging the amplitude intensity change curve L of the crack 1 and the crack 2 in the graph 1112And L23L of the crack 1, as shown by curves a) and b) in FIG. 12, respectively12And L23Connecting the peak points to obtain a leftward inclined arrow, and judging that the crack 1 is a leftward inclined crack; l of crack 212And L23And connecting the peak points to obtain an arrow in the plumb direction, and judging that the crack 2 is a vertical crack.
The above conclusion is verified by the in-situ core-drilling sampling method, and the fracture 1 is confirmed to be an inclined fracture (as shown in fig. 13) and the fracture 4 is confirmed to be a vertical fracture (as shown in fig. 14).
Claims (6)
1. A method for rapidly detecting the crack type of an asphalt pavement based on a ground penetrating radar is characterized by comprising the following steps:
step 1), collecting crack underground data by adopting a ground penetrating radar to obtain a oscillogram 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 the reflection amplitude of the electromagnetic waves between the materials of each layer;
step 4), counting the amplitude values of radar data of all crack positions within the time range of 0.7-1.3 ns below the homophase axis of the interface reflection wave between the asphalt concrete surface layer and the semi-rigid base layer in the ground penetrating radar section, and carrying out normalization processing to obtain the maximum value A12Obtaining the amplitude intensity variation curve L along the length of the measuring line12；
Step 5), counting the amplitude values of radar data of all crack positions within the time range of 1.7-2.7 ns below the homophase axis of the interface reflection wave between the semi-rigid base and the soil filling subgrade in the ground penetrating radar section, and carrying out normalization processing to obtain the maximum value A23Obtaining the vibration along the length of the measuring lineAmplitude intensity variation curve L23；
Step 6), along the amplitude intensity variation curve L12Respectively finding the peak value A of the amplitude intensity of the diffraction wave of the crack1,A2,…AiAlong the amplitude intensity curve L23Respectively finding the peak value point of the amplitude intensity of the diffraction wave of the crack
Step 7) according to A between different cracks12And A23Judging the crack to be a fatigue crack or a reflection crack;
step 8), according to the actual position and the amplitude intensity change curve L of the same crack12And L23And judging the crack to be a vertical crack or an inclined crack according to the corresponding position relation of the amplitude intensity peak point.
2. The method for rapidly detecting the type of the asphalt pavement crack based on the ground penetrating radar as claimed in claim 1, wherein the method comprises the following steps: in the step 1), data acquisition is carried out by adopting a ground penetrating radar along the direction perpendicular to the direction of the crack.
3. The method for rapidly detecting the type of the asphalt pavement crack based on the ground penetrating radar as claimed in claim 1, wherein the method comprises the following steps: in the step 1), the ground penetrating radar uses a 400MHz ground coupling antenna.
4. The method for rapidly detecting the type of the asphalt pavement crack based on the ground penetrating radar as claimed in claim 1, wherein the method comprises the following steps: in the step 2), the data editing comprises the steps of waste channel elimination, direction adjustment, drifting processing and channel equalization processing.
5. The method for rapidly detecting the type of the asphalt pavement crack based on the ground penetrating radar as claimed in claim 1, wherein the method comprises the following steps: the step 7) is carried out according to different cracks A12And A23The method for judging the crack as a fatigue crack or a reflection crack comprisesForming a sample library by all cracks detected in a minimum unit of 1km, wherein the number of the cracks in the sample library is not less than 20, randomly selecting 2 cracks in the sample library, and respectively naming the cracks as cracks 1 and 2, and when A of the crack 1 is the minimum, randomly selecting 2 cracks in the sample library, respectively naming the cracks 1 and 212A is not less than crack 223When the crack 1 is a fatigue crack, A of the crack 1 is judged23< A of crack 223Then, the crack 2 is judged to be a reflection crack.
6. The method for rapidly detecting the type of the asphalt pavement crack based on the ground penetrating radar as claimed in claim 1, wherein the method comprises the following steps: in the step 8), when the actual position and the amplitude intensity change curve L12And L23If the peak value of (a) is 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 true CN114624704A (en) | 2022-06-14 |
| CN114624704B 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 (4)
| 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 |
| US20190242246A1 (en) * | 2018-02-08 | 2019-08-08 | Saudi Arabian Oil Company | Mapping fracture length using downhole 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 |
-
2022
- 2022-02-17 CN CN202210146001.1A patent/CN114624704B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20190242246A1 (en) * | 2018-02-08 | 2019-08-08 | Saudi Arabian Oil Company | Mapping fracture length using downhole ground penetrating radar |
| 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 |
|---|
| 郭士礼;闫飞;朱培民;李修忠;: "裂缝宽度探地雷达波场响应的数值研究", 地球物理学进展, no. 04, 15 August 2016 (2016-08-15) * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN114624704B (en) | 2024-04-19 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN111476088B (en) | Asphalt pavement water damage identification model construction method, identification method and system | |
| CN106443668B (en) | The evaluation method of highway pavement splitting based on Ground Penetrating Radar | |
| CN112666554A (en) | Method for identifying radar amplitude characteristic crack width of asphalt pavement | |
| Hopman et al. | Repeatability, reproducibility and accuracy of GPR measurements | |
| CN115597530A (en) | Highway medium parameter inversion method | |
| CN112485789A (en) | Asphalt pavement compactness detection method based on three-dimensional ground penetrating radar | |
| CN108333096A (en) | A kind of asphalt concrete pavement porosity detection method based on Ground Penetrating Radar | |
| CN112684440A (en) | Crack evaluation method based on three-dimensional ground penetrating radar detection | |
| Barrett et al. | Understanding the capabilities of GPR for the measurement of ballast fouling conditions | |
| Wang et al. | Automatic asphalt layer interface detection and thickness determination from ground-penetrating radar data | |
| Long et al. | Road distress detection and maintenance evaluation based on ground penetrating radar | |
| CN111398687B (en) | Test method for estimating dielectric constant of asphalt pavement | |
| CN111965711B (en) | GPR image forward-modeling technique pavement reflection crack depth development simulation method | |
| CN114624704A (en) | Method for rapidly detecting crack type of asphalt pavement based on ground penetrating radar | |
| 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. | |
| Walia et al. | Reviewing methods for determination of Dielectric Constant required to Calibrate GPR Study for Asphalt Layers | |
| Dai et al. | Pavement evaluation using ground penetrating radar | |
| Li et al. | Non-destructive survey of pavement layer thicknesses with ground penetrating radar | |
| CN110187191A (en) | A kind of real-time control method for obtaining calibration material dielectric constant | |
| Huang et al. | Research on void signal recognition algorithm of 3D ground-penetrating radar based on the digital image | |
| CN115876661A (en) | Asphalt pavement void ratio evaluation method based on three-dimensional ground penetrating radar image | |
| CN115902884A (en) | Asphalt pavement internal condition evaluation method based on three-dimensional radar typical amplitude value | |
| Batrakov et al. | GPR signals processing for the road pavements monitoring with the weak contrast layers | |
| Lv et al. | Dielectric characteristics of unsaturated loess and the safety detection of the road subgrade based on GPR |
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 |