CN103149274B - A kind of concrete defect detection method - Google Patents

A kind of concrete defect detection method Download PDF

Info

Publication number
CN103149274B
CN103149274B CN201310031762.3A CN201310031762A CN103149274B CN 103149274 B CN103149274 B CN 103149274B CN 201310031762 A CN201310031762 A CN 201310031762A CN 103149274 B CN103149274 B CN 103149274B
Authority
CN
China
Prior art keywords
signal
linear
concrete
echoed
large aperture
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.)
Expired - Fee Related
Application number
CN201310031762.3A
Other languages
Chinese (zh)
Other versions
CN103149274A (en
Inventor
李长征
张碧星
师芳芳
吴先梅
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Institute of Acoustics CAS
Original Assignee
Institute of Acoustics CAS
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Institute of Acoustics CAS filed Critical Institute of Acoustics CAS
Priority to CN201310031762.3A priority Critical patent/CN103149274B/en
Publication of CN103149274A publication Critical patent/CN103149274A/en
Application granted granted Critical
Publication of CN103149274B publication Critical patent/CN103149274B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Abstract

The present invention relates to a kind of concrete defect detection method.Method comprises: the linear FM signal in broadband when the diverse location sequential firing in space is wide; Receive echoed signal; For the echoed signal received, broad pulse linear FM signal in time domain is compressed in the shorter time period, according to locus and phase relation, relevant process of making a reduction is carried out to the echoed signal of diverse location, form the large aperture image of equivalence, based on described large aperture image recognition concrete defect.Pulse compression technique is combined with synthetic aperture imaging detection method by the present invention, improves signal to noise ratio (S/N ratio) and the resolution of detection signal, is conducive to the detection of the especially less defect of concrete defect.

Description

A kind of concrete defect detection method
Technical field
The present invention relates to Acoustic detection field, particularly relate to a kind of concrete defect detection method.
Background technology
Current Ultrasonic Detection is widely used in engineering geophysics field, is especially applied in the defects detection of the mass concretes such as dam, bridge, tunnel.The accuracy of ultrasonic testing results is directly connected to safety and the quality evaluation result of engineering.
Pulse echo method is the detection that the mode utilizing single-shot list to receive carries out defect, and ultrasound wave propagation distance in tested component is longer, energy significantly decay after echoed signal be often submerged in noise, cause flaw indication not easy to identify.Sound beam focusing to specific region, thus is strengthened the signal to noise ratio (S/N ratio) of echoed signal by ultrasonic phase array.Phased-array technique is useful in raising penetration range and signal to noise ratio (S/N ratio), but this technology needs complicated hardware and software control system, in complex dielectrics array signal wave beam control in also there is certain difficulty.Meanwhile, Site Detection condition is comparatively complicated, and the transducer of larger aperture is laid exists certain difficulty.
There is not the signal to noise ratio (S/N ratio) and resolution that can improve Ultrasonic Detection at present, again the convenient scheme implemented.
Summary of the invention
For the problems referred to above of the prior art, the invention provides a kind of concrete defect detection method, described method comprises: the linear FM signal in transmitting transducer broadband when the diverse location sequential firing of concrete surface is wide; Receiving transducer receives echoed signal, and is compressed in the shorter time period by the broad pulse linear FM signal in the time domain in the echoed signal received; Carry out coherence stack process according to locus and phase relation to corresponding to the echoed signal after the compression of described diverse location, form the large aperture image of equivalence; Based on described large aperture image recognition concrete defect.
Preferably, broad pulse linear FM signal in time domain in the echoed signal received is compressed in the shorter time period by described receiving transducer, be specially: described broad pulse linear FM signal and reference signal carry out related calculation by receiving transducer, thus described broad pulse linear FM signal was compressed in the shorter time period.
Preferably, described reference signal obtains in the following manner: transmitting transducer directly docks with receiving transducer, and launch described wide time broadband linear FM signal give described receiving transducer, the signal that receiving transducer receives is described reference signal.
Pulse compression technique is combined with synthetic aperture imaging detection method by the embodiment of the present invention, implements comparatively simple, and improves signal to noise ratio (S/N ratio) and the resolution of detection signal, be conducive to the detection of the especially less defect of concrete defect.
Accompanying drawing explanation
Be described in detail specific embodiment of the invention scheme below with reference to accompanying drawings, advantage of the present invention will highlight more accordingly.In the accompanying drawings:
Fig. 1 is the detection scene schematic diagram of the embodiment of the present invention;
Fig. 2 be the embodiment of the present invention relate to containing decaying concrete test block/test specimen schematic diagram;
Fig. 3 is the Received signal strength schematic diagram received through 29 measuring points corresponding to transmitting LFM signal of the embodiment of the present invention;
Fig. 4 is the Received signal strength schematic diagram received through 29 measuring points corresponding to transmitting short pulse signal of the embodiment of the present invention;
Fig. 5 be the embodiment of the present invention about short pulse receiving signal and LFM compressed signal, for the 18th measuring point signal to noise ratio (S/N ratio) contrast schematic diagram;
Fig. 6 is the SAFT imaging results schematic diagram about short pulse excitation of the embodiment of the present invention;
Fig. 7 is the SAFT imaging results schematic diagram about LFM excitation of the embodiment of the present invention.
Embodiment
Below by embodiment, technical scheme of the present invention is described in further detail.
In view of the defect of existing scheme, consider pulse compression technique to be combined with synthetic aperture imaging detection method.
The diverse location sequential firing of described Synthetic Aperture Focussing Imaging (SAFT) in space also receives echoed signal, carries out coherence stack process according to locus and phase relation to the echoed signal of diverse location, thus forms the large aperture of equivalence.The method effectively can improve the signal to noise ratio (S/N ratio) of Ultrasonic Detection.
Described pulse compression technique is the linear FM signal in the broadband when transmitting terminal is launched wide, thus increase the energy transmitted, by the principle of correlation reception, the broad pulse linear FM signal in time domain is compressed within the very short time period at receiving end, both ensure that the resolution of compressed signal, improve again energy and the signal to noise ratio (S/N ratio) of signal largely, investigation depth is increased greatly.At present, pulse compression technique concentrates on underwater sound field and medical research field, and in engineering geophysics field, application is less.
Generally, technical scheme of the present invention realizes in the following manner:
1) choose applicable Site Detection condition and the transmitting transducer that resonance frequency band is wider and receiving transducer.
2) formula (1) is utilized to edit linear frequency modulation coded signal S code, f 0for FM signal centre frequency, sampling time interval is Δ t, and amplitude is A, and time width is T, and modulating bandwidth is B, and signal chirp rate is coded signal length (sampling number) is M 0, the centre frequency f of signal 0should be consistent with the centre frequency of transducer, the bandwidth of linear FM signal should be selected in transducer-6dB bandwidth range.
For ensureing that signal possesses good pulse compression characteristics, the pulse width T of linear FM signal and the product of bandwidth B should be greater than 10, and are less than 100.
S code ( iΔt ) = A cos [ 2 π ( f 0 - B 2 ) iΔt + μ ( iΔt ) 2 2 ] , - - - ( 1 )
0≤i Δ t≤T, i are positive integer, and Δ t is sampling interval.
3) by coded signal S code(i Δ t) lead-in signal generator is used for excitation, and transmitting transducer directly docks with receiving transducer, record Received signal strength S ref, S refas with reference to function, consider the delay of transducer, S refsampling length M 1be greater than M 0.
4) according to ultrasonic field technical norms for measurement, adopt reflectometry to detect, adopt the transducer apparatus of bistatic, T-R transducer distance is less than 2cm, in the position of the N number of equally spaced measuring point of concrete surface mark, and S codet () lead-in signal generator excitation transducer detects, record the echoed signal e of N number of measuring point i(t) (1≤i≤N), and record measuring point coordinate (x i, y i).For ensure launch wide time band signal be all received, Received signal strength length (sampling number) M 2should meet: M 2be greater than M 1with Δ Lsample_rate/V estimatesum.Wherein, Δ L is detected concrete thickness, and sample_rate is sample frequency, V estimatefor wave velocity estimation value, wherein coded signal should be consistent with the sample frequency sample_rate of collection signal.
5) formula (2) is utilized, reference function S refwith Received signal strength e icarry out related calculation, obtain N number of compressed signal if M 3=M 2+ M 1,
S c i ( m ) = Σ k = 1 M 1 S ref ( k ) e i ( k + m ) , 1 ≤ m ≤ M 3 , 1 ≤ i ≤ N , - - - ( 2 )
If (k+m) is greater than M in computation process 2, get e i(k+m)=0.
6) do not consider the effect of distance of sending and receiving transducer, the Gridding length of employing two-dimensional imaging is
Δ l=v p/ sample_rate, if concrete long l 1, high l 2, then transverse grid adds up to h number=l 1/ Δ l, longitudinal grid number is v number=l 2/ Δ l, obtains the image array I that initial value is 0 i(m, n), 1≤m≤h number, 1≤n≤V number.
Adopt cycle calculations, n1 is increased to M from 1 3, step-length is 1, if space networks lattice point (m, n) meets [(m Δ l-Δ l) 2+ (n Δ l-Δ l) 2] 1/2=k Δ t/2v p, (3)
Then, image array I i ( m , n ) = I i ( m , n ) + S c i ( nl ) , 1 ≤ i ≤ N - - - ( 4 )
M, nl, n are integer, 1≤nl≤M 3, v pfor concrete velocity of longitudinal wave.
7) to the compressed signal of all measuring points carry out the computing of (4) formula, and superpose, final image matrix can be obtained I sum ( m , n ) = Σ i = 1 N I i ( m , n ) - - - ( 5 ) .
According to the cross-sectional sizes of defect, set up image array I sumthe brightness of each grid (m, n) of (m, n) is by compressed signal ? the amplitude of moment ripple determines. that theoretic ripple is from focus to the travel-time of defect impact point.The compressed signal peak value of all measuring points is all by grid (m corresponding to defect impact point aim, n aim) carry out amplitude superposition, result is (m aim, n aim) place's amplitude maximum, on the contrary, zero defect grid place amplitude is less.The amplitude of all image lattices is by with different look chi imagings, and defect will show with bright band region.
Below in conjunction with accompanying drawing, the embodiment of the present invention is described in more details, so that those skilled in the art understand its principle and implementation detail better.
Fig. 1 is for detecting schematic diagram.For implementing described detection, make the concrete test block as Fig. 2, the long 35cm of this test block, high 29cm, thick 10cm.There is the rectangle defect of a long 10cm, high 3cm centre.Linear FM signal encourage, time wide 20 μ s, bandwidth range 250kHz-750kHz.
Particularly, following steps are carried out:
(1) according to test block size, test block is divided into transverse direction 176, vertical direction 145 grids.Tested by reality, test block longitudinal wave velocity is 3970m/s.Select transducer centre frequency 500kHz ,-6dB bandwidth 500kHz.
(2) 3 in technical scheme is utilized) obtain reference signal.
(3) sending and receiving Distance Between Transducers is 1cm, the left hand edge 3cm of first measuring point distance test block, successively the transducer of synchronizing moving transmitting-receiving to the right, and adjacent measuring point distance 1cm, detect 29 measuring points altogether, detection signal is as Fig. 3.
(4) after carrying out same-amplitude compacting to direct-path signal, utilize formula (2) to obtain the pulse compression signal of 29 measuring points, as Fig. 4, the signal to noise ratio (S/N ratio) of compressed signal improves degree in table 1.
(5) utilize detection scheme 6) and 7) obtain the imaging results shown in Fig. 5, by Fig. 5, imaging results and actual defects position comparatively identical.
For encouraging Detection results to contrast with short pulse, be the excitation of 500kHz short pulse by the frequency in 2 cycles, test mode is identical with impulse compression method test mode.Received signal strength is as Fig. 4.After the signal compression of the 18th measuring point that contrast LFM detects, comparatively short pulse detection signal-to-noise ratio improves about 8dB, as Fig. 5.
For contrast SAFT image effect, the signal amplitude within 12dB is only calculated during SAFT imaging, Fig. 6 and Fig. 7 is respectively the Received signal strength of short pulse excitation and the SAFT imaging results of LFM excitation compressed signal, two kinds of methods all can carry out imaging to the defective locations of inside concrete, the defect imaging position of Fig. 6 is a little less than the defective locations of Fig. 7, and the bright band re-imaging length of two figure all can show the lateral length of defect.Noise spot in Fig. 6 apparently higher than noise spot in Fig. 7, will show that impulse compression method can significantly improve the signal to noise ratio (S/N ratio) of defects detection.
To sum up, visible key point of the present invention and advantage as follows:
1) launch linear FM signal when having wide, broadband character, have good compression property, after compression, the signal to noise ratio (S/N ratio) of signal is higher.
2) transmit wide time, broadband character determines that the energy of signal is comparatively large, can penetrate the tested component of long distance.
3) basic waveform of linear FM signal is sinusoidal wave, and the length of signal does not have anything to limit,
Its pulse compression ratio can manual control, is is easily transmitted and received by existing sonac.
4) the detected image noise point calculated is few, easily to the identification of defect.
Those skilled in the art should recognize further, the software unit that the method step in conjunction with embodiment disclosed herein description can use hardware, processor to perform, or the combination of the two is implemented.Software unit can be placed in the storage medium of other form any known in random access memory (RAM), internal memory, ROM (read-only memory) (ROM), electrically programmable ROM, electrically erasable ROM, register, hard disk, moveable magnetic disc, CD-ROM or technical field.
It is to be noted, these are only present pre-ferred embodiments, not be used for limiting practical range of the present invention, the technician with professional knowledge base can realize the present invention by above embodiment, therefore every any change according to making within the spirit and principles in the present invention, amendment and improvement, all cover by the scope of the claims of the present invention.Namely, above embodiment is only in order to illustrate technical scheme of the present invention and unrestricted, although with reference to preferred embodiment to invention has been detailed description, those of ordinary skill in the art is to be understood that, can modify to technical scheme of the present invention or equivalent replacement, and not depart from the spirit and scope of technical solution of the present invention.

Claims (3)

1. a concrete defect detection method, is characterized in that, described method comprises:
The linear FM signal in transmitting transducer broadband when the diverse location sequential firing of concrete surface is wide;
Receiving transducer receives echoed signal, and is compressed in the shorter time period by the broad pulse linear FM signal in the time domain in the echoed signal received;
Carry out coherence stack process according to locus and phase relation to corresponding to the echoed signal after the compression of described diverse location, form the large aperture image of equivalence;
Based on described large aperture image recognition concrete defect;
Described linear FM signal is determined by following formula:
S code ( iΔt ) = A cos [ 2 π ( f 0 - B 2 ) iΔt + μ ( iΔt ) 2 2 ]
Wherein, described f 0for described linear FM signal centre frequency, described Δ t is sampling time interval, and described A is amplitude, and described B is modulating bandwidth, and 0≤i Δ t≤T, i are positive integer, and described T is time width, and the product of T and B is greater than 10, and is less than 100, and signal chirp rate is coded signal length is M 0;
Echoed signal after described compression is determined by following formula:
S c i ( m ) = Σ k = 1 M 1 S ref ( k ) e i ( k + m ) , 1 ≤ m ≤ M 3 , 1 ≤ i ≤ N
Wherein, described S reffor Received signal strength, M 3=M 1+ M 2, M 1for signal S refsampling length, M 1be greater than M 0, M 2for Received signal strength length, N is the number of the equidistant point position at concrete surface mark, described e ifor the broad pulse linear FM signal in the time domain in the described echoed signal received, if k+m>M in computation process 2, then e i(k+m)=0;
Describedly carry out coherence stack process according to locus and phase relation to the echoed signal after the compression corresponding to described diverse location, the large aperture image forming equivalence specifically comprises:
Utilize the echoed signal after described compression, obtain image array, overlap-add procedure is carried out to described image array, obtain the large aperture image of equivalence; Described image array is determined by following formula;
I nl i ( m , n ) = I i nl - 1 ( m , n ) + S c i ( nl ) ,
Wherein, the Gridding length adopting two-dimensional imaging is Δ l=v p/ sample_rate, if concrete long l 1, high l 2, then transverse grid adds up to h number=l 1/ Δ l, longitudinal grid number is v number=l 2/ Δ l, image array I i(m, n) initial value is 0,1≤i≤N, 1≤m≤h number, 1≤n≤v number, m, nl, n are integer, 1≤nl≤M 3, v pfor concrete velocity of longitudinal wave; Adopt cycle calculations, n1 is increased to M from 1 3, step-length is 1, if space networks lattice point (m, n) meets [(m Δ l-Δ l) 2+ (n Δ l-Δ l) 2] 1/2=k Δ t/2v p; Then described large aperture image is determined by following formula:
I sum ( m , n ) = Σ i = 1 N I i ( m , n ) .
2. method according to claim 1, it is characterized in that, broad pulse linear FM signal in time domain in the echoed signal received is compressed in the shorter time period by described receiving transducer, be specially: described broad pulse linear FM signal and reference signal carry out related calculation by receiving transducer, thus described broad pulse linear FM signal was compressed in the shorter time period.
3. method according to claim 2, it is characterized in that, described reference signal obtains in the following manner: transmitting transducer directly docks with receiving transducer, and launch described wide time broadband linear FM signal give described receiving transducer, the signal that receiving transducer receives is described reference signal.
CN201310031762.3A 2013-01-28 2013-01-28 A kind of concrete defect detection method Expired - Fee Related CN103149274B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201310031762.3A CN103149274B (en) 2013-01-28 2013-01-28 A kind of concrete defect detection method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201310031762.3A CN103149274B (en) 2013-01-28 2013-01-28 A kind of concrete defect detection method

Publications (2)

Publication Number Publication Date
CN103149274A CN103149274A (en) 2013-06-12
CN103149274B true CN103149274B (en) 2015-10-07

Family

ID=48547478

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201310031762.3A Expired - Fee Related CN103149274B (en) 2013-01-28 2013-01-28 A kind of concrete defect detection method

Country Status (1)

Country Link
CN (1) CN103149274B (en)

Families Citing this family (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR102193234B1 (en) * 2013-11-25 2020-12-18 현대모비스 주식회사 Device for increasing the detection distance of ultrasonic sensors of vehicle and method thereof
CN104535615B (en) * 2015-01-18 2018-03-13 无锡双马钻探工具有限公司 A kind of cross-correlation vortex thermal imaging imperfection detection and chromatography imaging method and system
CN106546661B (en) * 2016-10-24 2018-07-27 河海大学常州校区 A kind of interference ultrasonic synthetic aperture is to inside concrete imaging method
CN109521099B (en) * 2018-01-25 2021-12-28 四川升拓检测技术股份有限公司 Mobile nondestructive testing method based on phase difference array of sound pick-up
CN108415020A (en) * 2018-02-14 2018-08-17 中北大学 A kind of improved time-domain ultrasonic signal synthetic aperture algorithm
CN110243945B (en) * 2019-07-23 2021-04-30 大连理工大学 Ultrasonic TOFD blind area suppression method based on synthetic aperture focusing and mode converted wave
CN111044569B (en) * 2019-12-27 2023-01-24 中铁建大桥工程局集团第二工程有限公司 Tunnel concrete structure defect detection method

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6009755A (en) * 1996-11-08 2000-01-04 Mitsubishi Denki Kabushiki Kaisha Ultrasonic transceiver displaying modified B scope
CN1731214A (en) * 2005-08-25 2006-02-08 中国船舶重工集团公司第七一五研究所 Synthetic aperture three-dimensional acoustic imaging method
CN1996047A (en) * 2006-12-07 2007-07-11 中国人民解放军海军工程大学 Single receiving array synthetic-aperture sonar interferometry
CN101644763A (en) * 2008-12-12 2010-02-10 中国科学院声学研究所 Real-time signal-processing method of synthetic aperture sonar based on cluster processor and system
CN101825709A (en) * 2009-12-08 2010-09-08 中国科学院声学研究所 Underwater high-resolution side-looking acoustic imaging system and method thereof
CN102426384A (en) * 2011-09-06 2012-04-25 赵永贵 Method for detecting underground goaf and karst distribution
CN102788845A (en) * 2012-09-04 2012-11-21 黄河水利委员会黄河水利科学研究院 Barker coding excitation ultrasonic detection method of concrete structure defect

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6009755A (en) * 1996-11-08 2000-01-04 Mitsubishi Denki Kabushiki Kaisha Ultrasonic transceiver displaying modified B scope
CN1731214A (en) * 2005-08-25 2006-02-08 中国船舶重工集团公司第七一五研究所 Synthetic aperture three-dimensional acoustic imaging method
CN1996047A (en) * 2006-12-07 2007-07-11 中国人民解放军海军工程大学 Single receiving array synthetic-aperture sonar interferometry
CN101644763A (en) * 2008-12-12 2010-02-10 中国科学院声学研究所 Real-time signal-processing method of synthetic aperture sonar based on cluster processor and system
CN101825709A (en) * 2009-12-08 2010-09-08 中国科学院声学研究所 Underwater high-resolution side-looking acoustic imaging system and method thereof
CN102426384A (en) * 2011-09-06 2012-04-25 赵永贵 Method for detecting underground goaf and karst distribution
CN102788845A (en) * 2012-09-04 2012-11-21 黄河水利委员会黄河水利科学研究院 Barker coding excitation ultrasonic detection method of concrete structure defect

Also Published As

Publication number Publication date
CN103149274A (en) 2013-06-12

Similar Documents

Publication Publication Date Title
CN103149274B (en) A kind of concrete defect detection method
CN102809610B (en) Phased array ultrasonic testing method based on improved dynamic depth focusing
Prado et al. Lamb mode diversity imaging for non-destructive testing of plate-like structures
CN102539530B (en) Method for measuring echo reduction/reflection coefficient of underwater sound passive material based on time reversal focusing
CN103462643B (en) Shear wave speed measurement method, device and system
CN104007176B (en) Full-wave field detection system and method of complex geotechnical engineering medium
JP2012225651A (en) Underwater detection device, fish finder, underwater detection method and program
CN103235039A (en) Ultrasonic TOFD (Time of Flight Diffraction) detection method for weld defects based on linear frequency modulation technology
Lin et al. Excitation waveform design for Lamb wave pulse compression
Ganguli et al. Synthetic aperture imaging for flaw detection in a concrete medium
KR101922522B1 (en) Sound speed imaging using shear waves
CN104634876A (en) Method for detecting inclusions in metal material by virtue of ultrasonic scanning microscope
CN110361453B (en) Array type ectopic receiving ultrasonic synthetic aperture focusing method based on nondestructive testing
CN105699486B (en) A kind of detection method of inclination surface crack inclination angle degree and depth
CN110346453A (en) Defect minispread reflection echo rapid detection method in a kind of concrete structure
CN109991590B (en) System and method for testing low-frequency emission characteristic of transducer in pressure tank in limited space
CN104698089A (en) Ultrasonic relative time propagation technology suitable for inclined crack quantifying and imaging
CN109374748A (en) Antidirection finding imaging method when acoustic emission source
CN104605890A (en) Shear wave crest value waveform correction method, device and system and application thereof
CN102788845A (en) Barker coding excitation ultrasonic detection method of concrete structure defect
CN110927252B (en) Targeted shear wave elastography detection method
RU2422769C1 (en) Method of ultrasound echo-pulse thickness measurement
CA3110818A1 (en) Continuous wave ultrasound or acoustic non-destructive testing
CN104777227A (en) Internal concrete defect detection method based on principles of three views
Frazier et al. Acoustic imaging of objects buried in soil

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20151007

Termination date: 20220128

CF01 Termination of patent right due to non-payment of annual fee