CN105806949A - Correcting method for radial travel time in ultrasonic CT detection - Google Patents
Correcting method for radial travel time in ultrasonic CT detection Download PDFInfo
- Publication number
- CN105806949A CN105806949A CN201610127591.8A CN201610127591A CN105806949A CN 105806949 A CN105806949 A CN 105806949A CN 201610127591 A CN201610127591 A CN 201610127591A CN 105806949 A CN105806949 A CN 105806949A
- Authority
- CN
- China
- Prior art keywords
- ray
- ultrasonic signal
- point
- launch point
- traveltime
- 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
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N29/00—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
- G01N29/04—Analysing solids
- G01N29/06—Visualisation of the interior, e.g. acoustic microscopy
- G01N29/0654—Imaging
- G01N29/0672—Imaging by acoustic tomography
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N29/00—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
- G01N29/44—Processing the detected response signal, e.g. electronic circuits specially adapted therefor
- G01N29/4463—Signal correction, e.g. distance amplitude correction [DAC], distance gain size [DGS], noise filtering
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2291/00—Indexing codes associated with group G01N29/00
- G01N2291/02—Indexing codes associated with the analysed material
- G01N2291/023—Solids
- G01N2291/0232—Glass, ceramics, concrete or stone
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2291/00—Indexing codes associated with group G01N29/00
- G01N2291/02—Indexing codes associated with the analysed material
- G01N2291/028—Material parameters
- G01N2291/0289—Internal structure, e.g. defects, grain size, texture
Landscapes
- Physics & Mathematics (AREA)
- Biochemistry (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Acoustics & Sound (AREA)
- Engineering & Computer Science (AREA)
- Signal Processing (AREA)
- Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
Abstract
The invention provides a correcting method for radial travel time in ultrasonic CT detection. The correcting method is oriented to ultrasonic detection of concrete with a heterogeneity ecderon layer, wherein the reflection echo of an ultrasonic signal on an interface between the ecderon layer and the concrete is collected, the first wave travel time of the reflection echo is utilized to correct the radial travel time and the part transmitted on the ecderon layer in the radial travel time is eliminated, so that the influence of the ecderon layer on the radial travel time after being corrected can be eliminated and the signal transmission time of the ultrasonic signal in the concrete can be accurately reflected. The correcting method provided by the invention has the characteristics that the ultrasonic signal with time-varying characteristics is utilized to overcome the influences of signal attenuation, various echoes and bottom waves in the process of detecting the reflection echo of the signal emitted at each moment.
Description
Technical field
The present invention relates to the concrete nondestructive testing technology realized by ultrasound computed tomography, particularly relate to the ray traveltime bearing calibration in the detection of a kind of ultrasound computed tomography.
Background technology
Inside at concrete structure, likely because there is cavity, crack or being mingled with silt foreign material and produce leakiness district, and because leakage is shaken, isolated or cell that is built on stilts and that formed, or loose shape region is formed because lacking cement, and because being caused loose shape region by accidental injury.These defects destroy the seriality within concrete structure and integrity, and reduce concrete bearing strength and durability to a great extent.Effectively find out the position of drawbacks described above, scope and size in concrete structure, be the precondition carrying out toolability technical finesse for these defects.
The one quantification detection method accurate, efficient formed after the detection of concrete ultrasound computed tomography is concrete Ultrasonic Nondestructive and computerized tomography (CT) be combined with each other.In concrete ultrasound computed tomography detection process, transmitting transducer and reception transducer are separately fixed on two apparent surfaces of concrete structure, by transmitting transducer repeat its transmission ultrasonic pulse, ultrasonic pulse, at concrete structure internal communication, is then received by reception transducer.Ultrasonic pulse signal is in the communication process within concrete structure, when running into rejected regions such as there is honeycomb, cavity or crack, path can only be changed walk around these rejected regions and could continue to propagate, increase when causing walking because transmission path increases and the velocity of sound reduces;Reflection and scattering occurring in the interface of rejected region, cause that ultrasonic wave acoustic energy is decayed, the amplitude of ultrasonic signal is decreased obviously;Because the high fdrequency component in ultrasonic signal decays more significantly, thus ultrasonic signal frequency reduces under the effect of rejected region;Because of the impact of path changing and scattering, reflection, ultrasonic signal is overlapped mutually and makes signal waveform produce distortion.Thus, by gather receive transducer and receive when walking of ultrasonic signal, amplitude, the parameter such as frequency, it is possible to carry out the analysis of concrete quality.On this basis, by adopting ultrasound wave ray that is intensive and that intersect to penetrate each section of concrete structure, then inverting is carried out according to parameters such as when walking of ultrasonic signal, amplitude, frequencies, rebuild the concrete inner structure status image reflecting each section, carry out quantitative measurement and the evaluation of rejected region in combining image process and mode identification technology.
It is the key component that concrete ultrasound computed tomography detects that detection ultrasonic signal carries out inversion imaging at parameter when walking and when passing through to walk concrete structure within.Ultrasonic signal extends transmission path when running into rejected region because detouring, and causes signal to pass through increasing when walking of this position, and the velocity of wave thus corresponding to rejected region is slack-off.Therefore, by measure ultrasonic signal on each section when walking, recycle the wave speed distribution situation on suitable algorithm inversion reckoning section, the position of rejected region, size, distribution and type can be reflected according to the difference of wave speed distribution.
The ultimate principle of detection and inversion imaging when Fig. 1 illustrates away.On the detection section of shown in Fig. 1 a ultrasound computed tomography, S1-Sn is the launch point (launch point such as transmitting transducer) of ultrasonic signal, R1-Rn is the reception point of ultrasonic signal, it is assumed that ultrasonic signal is linearly propagated between each launch point and each reception point, thus as it is shown in figure 1, there is a ray between each launch point and reception point, amounts to N bar ray, and it is 1-N that ray sequence number is compiled.From receiving t when walking that can record every ray1-tN, when walking and ultrasonic signal from launch point along ray propagate receive point used time.Further, this section is decomposed into m image-generating unit, in Fig. 1, has split m=8*8 amount to 64 image-generating units, be numbered 1-m;The quantity of m is determined by imaging precision, and, when the numerical value of m ensures that each image-generating unit is sufficiently small, ultrasound wave can be approximated to be a fixed value in the slowness (slowness is the inverse of velocity of wave) of each image-generating unit, is expressed as f1-fm.Thus, for wherein i-th ray, t when it is walkediCan be expressed as:Wherein ajiRepresent this i-th ray length of process in jth image-generating unit.ajiEngineering construction drawing or the angle of in-site measurement and ray according to concrete structure are readily available, therefore, and t when above formula also just characterizes awayiAnd the operation relation between the slowness (namely velocity of wave) on each image-generating unit, for utilize detectable walk time each image-generating unit of inverting velocity of wave lay the foundation, and the travel-time tomography of ultrasound computed tomography is exactly that velocity of wave distribution on each image-generating unit is showed in the form of images.
For when the walking of whole rays recorded, it is possible to obtain following Matrix Formula:
By length aiiIn the coefficient matrix of composition, owing to every ray only can by the sub-fraction in whole m image-generating units, therefore the most elements in this matrix will value be 0.For the feature that coefficient matrix is sparse, it is possible to use the iterative algorithm such as ART realizes based on the inverting to wave speed distribution of measured value when walking.
T when walking to whole rays1-tNMeasurement, be inverting obtain wave speed distribution the basis of rejected region in concrete structure can be shown real table.In the middle of practice, arriving the time receiving point as when walking being transmitted into Mintrop wave from ultrasonic signal, the degree of accuracy measured when walking needs to reach Microsecond grade, even up to 1/10th microseconds.
Under desirable measuring state, the emitter of ultrasonic signal and receptor should be close to concrete structure surface, and at joint place daubing coupling agent to avoid the occurrence of space.But, in some cases, there is heterogeneous serving in concrete surface, such as the top layer etc. that outer surface swabbing etc. is formed, thus, in this case the measured ray traveltime obtained is strictly speaking, be not ultrasonic signal propagate in concrete structure when walking, but by its being also included within when walking in the middle of serving.Owing to prior art cannot accurately eliminate ultrasonic signal part when walking in serving, and examining and making cuts when can only walk with an estimated value based on serving thickness and prediction velocity of wave, obtain ultrasonic signal in concrete structure when walking, this obviously deviates from mutually with requirement of measurement pinpoint accuracy when walking.
Summary of the invention
For the defect overcoming above-mentioned prior art to exist, it is desirable to provide the ray traveltime bearing calibration in the detection of a kind of ultrasound computed tomography.The present invention is towards the concrete ultrasound detection with heterogeneous serving, wherein gather ultrasonic signal reflection echo at interface between serving and concrete, utilize reflection echo head wave traveltime that ray traveltime is corrected, eliminate the part in serving transmission in ray traveltime, so that the ray traveltime after correction can get rid of the impact of serving, accurately reflect the ultrasonic signal signal transmission time at inside concrete.Utilize during another feature of the present invention and there is the ultrasonic signal of time variation feature, detect each moment launch signal reflection echo process in overcome signal attenuation and the impact of each echo, end ripple time domain superposition.
Ray traveltime bearing calibration in ultrasound computed tomography of the present invention detection, it is characterised in that including:
At the launch point of serving surface configuration ultrasonic signal, launch point launches the ultrasonic signal with time variation characteristic parameter one by one;Wherein, the described time variation characteristic parameter that ultrasonic signal has at each x time of launch point is different;
Point is received in the side surface configuration first group relative with the surface arranging launch point;The described point that receives in first group measures the time of the signal Mintrop wave receiving described ultrasonic signal, as the ray traveltime connecting described launch point and this reception point;
The described surface that launch point is set arranges second group and receives point;Described reception point in second group receives and goes out the described ultrasonic signal reflection echo that interface produces between described serving and xoncrete structure according to described time variation characteristic parameter extraction, and metering receives the time of the signal Mintrop wave of reflection echo as when correcting;
The spacing between thickness and reception point and the launch point of second group according to described serving, obtain the reflection echo ray put that receives of described launch point to second group, and and then draw each reflection echo ray and from launch point to the ray angle received between the ray put of first group;
According to the ray angle calcu-lation modifying factor as this ray angle linear function;Ray traveltime correction it is converted into based on described modifying factor when described correction being walked;
Utilize described ray traveltime correction, apply to revise to the ray traveltime of the ray receiving point of launch point to first group.
Preferably, described time variation characteristic parameter is the mid frequency of ultrasonic signal.It may further be preferable that launch point launches ultrasonic signal with predetermined fixed center frequency, and the mid frequency of the ultrasonic signal of each moment transmitting of launch point is time-varying;Second group receive is put received ultrasonic signal and records the mid frequency of echo-signal;Judge that the skew of mid frequency of the ultrasonic signal of echo-signal mid frequency and the transmitting of particular transmission moment is whether within estimating threshold range, if within threshold range, then determine and receive the reflection echo that the ultrasonic signal that receives of point comprises the ultrasonic signal that this x time is launched, and this reception point is detected first reflection echo when walking as described correction walk time.
Preferably, described time variation characteristic parameter is wavelet-based attribute vector.It may further be preferable that frequency range when being decomposed into N number of by the ultrasonic signal that each for launch point moment launches, after Denoising disposal, calculate ultrasonic signal energy coefficient in frequency range when each;Energy distribution coefficient table is shown as characteristic vector (E1, E2 ... EN), and the characteristic vector of the ultrasonic signal of each moment transmitting of launch point is time-varying;Second group receive is put received ultrasonic signal and decomposes the characteristic vector that acquisition receives;The relatively characteristic vector of the ultrasonic signal that particular moment launches and the degree of association receiving the characteristic vector that point receives;When the two degree of association is more than threshold value, it is determined that receive the reflection echo that the ultrasonic signal that receives of point comprises the ultrasonic signal that this particular moment launches, and this reception point is detected first reflection echo when walking as described correction walk time.
Preferably, the revised ray traveltime of each ray is substituted into following Matrix Formula:
In formula, t1’-tn' for revising the ray traveltime after just;ajiRepresent this i-th ray length of process in jth image-generating unit;Concrete body A is divided m image-generating unit, and ultrasonic signal is approximately fixed value f in the slowness of each image-generating unit1-fm;Utilize the iterative algorithm such as ART realize based on revised walk time the numerical value inverting to wave speed distribution.
Preferably, for arbitrary launch point S1 ultrasonic signal launched, metering connects this launch point and first group n the ray traveltime t receiving point1-tn;
For this launch point S1 ultrasonic signal launched, measure this launch point to second group n and receive t when the correction put is walkedC1-tCn;
The thickness HB of serving can be obtained according to engineering construction drawing or in-site measurement, and each spacing receiving point and this launch point according to second group, obtain reflection echo ray, and then draw each reflection echo ray and from the ray angle α between first group of ray respectively received a little of this launch point1-αn;
Calculate the modifying factor β as this ray angle linear function1-βn;
Based on modifying factor, it is possible to t when correction is walkedC1-tCnIt is converted into ray traveltime correction, it may be assumed that
Δt1=β1·tC1
Δt2=β2·tC2
…
Δtn=βn·tCn
Utilize above ray traveltime correction, this launch point to first group is respectively received whole n ray traveltime t a little1-tnApply to revise, revised ray traveltime t1’-tn' as follows:
t1'=t1-Δt1
t2'=t2-Δt2
…
tn'=tn-Δtn
Revised ray traveltime t1’-tn' as ultrasonic signal within concrete body A actual walk time.
Thus, the ultrasonic signal that the present invention has time variation feature in utilization correctly detects on the basis of reflection echo that each moment launches signal, accurately determine when correction is walked, and then according to above projectional technique, the ray traveltime in xoncrete structure ultrasound computed tomography detection process is carried out rate correction.Ray traveltime after calibrated can get rid of the impact of serving, accurately reflects the ultrasonic signal signal transmission time at inside concrete, thus for setting up good measurement data basis towards the wave speed distribution inverting of inside concrete and imaging.
Accompanying drawing explanation
Below in conjunction with the drawings and specific embodiments, the present invention is further detailed explanation:
The basic principle schematic of detection and inversion imaging when Fig. 1 is to walk;
Fig. 2 is the some position schematic diagram of launch point and reception point in ray traveltime bearing calibration of the present invention;
Fig. 3 is rate correction principle schematic of the present invention.
Detailed description of the invention
By the examples below, technical scheme is done more specifically bright.
The point position schematic diagram of launch point and reception point in ray traveltime bearing calibration of the present invention shown in Figure 2.Detect on section shown in Fig. 2, it is seen that structural member includes concrete body A and serving B.Serving B and concrete body A is mutually heterogeneous, thus forms interface between, and this interface is expressed as A-B.S1-Sn is the launch point (launch point such as transmitting transducer) of ultrasonic signal;R1-Rn is the reception point measured when walking for realizing ultrasonic signal, R1-Rn is arranged at the side surface that the surface at launch point place is relative.Further, in order to correct serving B to R1Impact when-Rn records away, at launch point S1The phase the same side, surface at-Sn place is arranged to the reception point C of correction1-Cn。
In ultrasound computed tomography detection process, launch point S1-Sn launches ultrasonic signal one by one, when each launch point launches ultrasonic signal, and each reception point R1-Rn metering receives the time of signal Mintrop wave, as connecting this launch point and each ray traveltime receiving point.For example, for S1Point R is received to each1Whole n bar rays of-Rn, the ray traveltime that metering obtains is t1-tn。
Ultrasonic signal is incident to interface A-B at a certain angle in communication process, it may occur that reflection, produces reflection echo, and reflection echo can by receiving some C1When-Cn receives and determines that it is walked the time of advent according to Mintrop wave, it is referred to as correction when walking;Such as, for S1The ultrasonic signal launched, receives some C1The correction that-Cn obtains is referred to as t when walkingC1-tCn.And ultrasonic signal is at interface A-B it also occur that refraction, and the ultrasonic signal after refraction continues to propagate in concrete body A, until it reaches with the opposite side surface, surface at launch point place, it is achieved receive some R1The signal of-Rn receives, it is thus achieved that described ray traveltime t1-tn.The described opposite side surface, surface with launch point place also may proceed to reflectance ultrasound ripple signal, and the ultrasonic signal of reflection, as end ripple first, returns to the surface of emitter and received some C1-Cn receives;By that analogy, some C is received1-Cn also can receive ripple etc. at the bottom of second trip echo and secondary, and each time echo can progressively be decayed with end ripple.
Rate correction principle schematic shown in Figure 3, can obtain the thickness HB of serving B, and receive some C according to each according to engineering construction drawing or in-site measurement1-Cn and launch point S1Spacing, it is easy to obtain shown in Fig. 3 from launch point S1Point C is received to each1The reflection echo ray of-Cn, and and then draw each reflection echo ray with from launch point S1Point R is received to each1Ray angle α between the ray of-Rn1-αn.In order to avoid irrelevance is excessive, described ray angle should be maintained in the scope of 3.5-15.5 degree.As seen from Figure 3, launch point S1Point R is received to each1There is fixing progressive error with reflection echo ray in the ray of-Rn, when thickness HB determines, ray angle is poor with this fixed journey system is corresponding, therefore, it can, according to ray angle, calculate the modifying factor β as this ray angle linear function1(α1)-βn(αn).Based on this modifying factor, it is possible to t when correction is walkedC1-tCnIt is converted into ray traveltime correction, it may be assumed that
Δt1=β1·tC1
Δt2=β2·tC2
…
Δtn=βn·tCn
Utilize above ray traveltime correction, to S1Point R is received to each1Whole n ray traveltime t of-Rn1-tnApply to revise, revised ray traveltime t1’-tn' as follows:
t1'=t1-Δt1
t2'=t2-Δt2
…
tn'=tn-Δtn
Revised ray traveltime t1’-tn' as ultrasonic signal within concrete body A actual walk time.And then, it is possible to for launch point S1Each launch point and R in-Sn1Whole N bar rays between each reception point of-Rn, calculate the revised ray traveltime of each ray, then substitute into following Matrix Formula:
In formula, ajiRepresent this i-th ray length of process in jth image-generating unit;Concrete body A is divided m image-generating unit, as in figure 2 it is shown, ultrasound wave is approximately fixed value f in the slowness (slowness is the inverse of velocity of wave) of each image-generating unit1-fm。
Utilize the iterative algorithm such as ART realize based on revised walk time the numerical value inverting to wave speed distribution.Difference according to wave speed distribution can reflect the position of rejected region, size, distribution and type, and wave speed distribution imaging is showed.
Ultrasound computed tomography detects duration, any point of geodesic structure is existed incidence wave (launch point S all simultaneously1-Sn directly provides ultrasonic signal) and various reflection echo, thus each receives superposing of incidence wave that point receives and reflection echo.Wherein, some R is received1-Rn receives and measures launch point S1The incidence wave that-Sn directly provides, because of incidence wave because conducting path is short, thus arrives and receives some R1-Rn's is little when walking, and signal attenuation is also smaller.But, with a reception point R1-Rn receives different with measuring incidence wave, receives some C1-Cn receives and measures the head wave traveltime of the reflection echo produced by interface A-B, the signal attenuation of reflection echo is obvious, and long relative to incidence wave during because walking, it is not easily distinguishable time domain is mixed in together with ripple at the bottom of other each echo and each time, this just gives and receives a some C1T when correction is walked by-CnC1-tCnMeasurement bring obstacle.
For solving this difficult problem, present invention employs feature and there is the ultrasonic signal of time variation as the detection signal of application in ultrasound computed tomography detection and rate correction.Such as, for a certain moment T1 ultrasonic signal launched, the Mintrop wave of its reflection echo arrives and receives some C1Moment be T2;Receive some C1Each echo and the end ripple signal of signal is also launched at T2 reception to other moment, but these echoes or end ripple signal are different from T1 due to x time, therefore the feature of its time variation is necessarily different from reflection echo Mintrop wave produced by the ultrasonic signal that the T1 moment launches.Thus, it is possible to by extracting time variation feature, the reflection echo Mintrop wave accurately obtaining T1 moment ultrasonic signal from the signal that time domain mixes arrives t when the correction receiving some C1 is walkedC1, thus avoiding error detection or inspection does not measure the Mintrop wave of reflection echo.
The mid frequency of ultrasonic signal can be adopted as the parameter of the time variation feature of ultrasonic signal.At particular moment T1, launch point launches ultrasonic signal with predetermined fixed center frequency;For the T2 reception point C1 reflection echo signal received, power spectrum method is adopted to record the mid frequency of echo-signal;Because the decay of reflection echo high frequency components is big, low frequency component decay is little, thus mid frequency can offset.Judge that the skew of mid frequency is whether within estimating threshold range, if within threshold range, determine and receive the reflection echo that the some C1 ultrasonic signal received comprises the ultrasonic signal that the T1 moment launches;T when walking as correction when walking of detection reflection echo firstC1.The mid frequency of the ultrasonic signal launched due to each moment is time-varying, thus each the echo of ultrasonic signal that beyond T1, other moment launches and end ripple be not because possessing the mid frequency within threshold range, because the reflection echo without being taken as the ultrasonic signal that the T1 moment launches is detected.
In order to improve reflection echo detection reliability further, the parameter of the time variation feature of ultrasonic signal can select wavelet-based attribute vector.Frequency range when being decomposed into N number of by ultrasonic signal, after Denoising disposal, calculates ultrasonic signal energy coefficient in frequency range when each;Energy distribution coefficient table is shown as characteristic vector (E1, E2 ... EN), the degree of association of characteristic vector with the characteristic vector of the T2 reception point C1 reflection echo received by comparing the particular moment T1 ultrasonic signal launched, when the two degree of association is more than threshold value, it is determined that receive the reflection echo that the some C1 ultrasonic signal received comprises the ultrasonic signal that the T1 moment launches;T when walking as correction when walking of detection reflection echo firstC1.Time variation because of the ultrasonic signal that each moment launches, thus each the echo of ultrasonic signal that beyond T1, other moment launches and end ripple be not because possessing the characteristic vector within relevance threshold, because the reflection echo without being taken as the ultrasonic signal that the T1 moment launches is detected.
Thus, the ultrasonic signal that the present invention has time variation feature in utilization correctly detects on the basis of reflection echo that each moment launches signal, accurately determine when correction is walked, and then according to above projectional technique, the ray traveltime in xoncrete structure ultrasound computed tomography detection process is carried out rate correction.Ray traveltime after calibrated can get rid of the impact of serving, accurately reflects the ultrasonic signal signal transmission time at inside concrete, thus for setting up good measurement data basis towards the wave speed distribution inverting of inside concrete and imaging.
Above example is merely to illustrate the present invention; and it is not limitation of the present invention; those of ordinary skill about technical field; without departing from the spirit and scope of the present invention; can also make a variety of changes and modification; therefore all equivalent technical schemes fall within scope of the invention, and the scope of patent protection of the present invention should be defined by the claims.
Claims (7)
1. the ray traveltime bearing calibration in a ultrasound computed tomography detection, it is characterised in that including:
At the launch point of serving surface configuration ultrasonic signal, launch point launches the ultrasonic signal with time variation characteristic parameter one by one;Wherein, the described time variation characteristic parameter that ultrasonic signal has at each x time of launch point is different;
Point is received in the side surface configuration first group relative with the surface arranging launch point;The described point that receives in first group measures the time of the signal Mintrop wave receiving described ultrasonic signal, as the ray traveltime connecting described launch point and this reception point;
The described surface that launch point is set arranges second group and receives point;Described reception point in second group receives and goes out the described ultrasonic signal reflection echo that interface produces between described serving and xoncrete structure according to described time variation characteristic parameter extraction, and metering receives the time of the signal Mintrop wave of reflection echo as when correcting;
The spacing between thickness and reception point and the launch point of second group according to described serving, obtain the reflection echo ray put that receives of described launch point to second group, and and then draw each reflection echo ray and from launch point to the ray angle received between the ray put of first group;
According to the ray angle calcu-lation modifying factor as this ray angle linear function;Ray traveltime correction it is converted into based on described modifying factor when described correction being walked;
Utilize described ray traveltime correction, apply to revise to the ray traveltime of the ray receiving point of launch point to first group.
2. the ray traveltime bearing calibration in ultrasound computed tomography according to claim 1 detection, it is characterised in that described time variation characteristic parameter is the mid frequency of ultrasonic signal.
3. the ray traveltime bearing calibration in ultrasound computed tomography according to claim 2 detection, it is characterised in that launch point launches ultrasonic signal with predetermined fixed center frequency, and the mid frequency of the ultrasonic signal of each moment transmitting of launch point is time-varying;Second group receive is put received ultrasonic signal and records the mid frequency of echo-signal;Judge that the skew of mid frequency of the ultrasonic signal of echo-signal mid frequency and the transmitting of particular transmission moment is whether within estimating threshold range, if within threshold range, then determine and receive the reflection echo that the ultrasonic signal that receives of point comprises the ultrasonic signal that this x time is launched, and this reception point is detected first reflection echo when walking as described correction walk time.
4. the ray traveltime bearing calibration in ultrasound computed tomography according to claim 1 detection, it is characterised in that described time variation characteristic parameter is wavelet-based attribute vector.
5. the ray traveltime bearing calibration in ultrasound computed tomography according to claim 4 detection, it is characterized in that, frequency range when being decomposed into N number of by the ultrasonic signal that each for launch point moment launches, calculates ultrasonic signal energy coefficient in frequency range when each after Denoising disposal;Energy distribution coefficient table is shown as characteristic vector (E1, E2...EN), and the characteristic vector of the ultrasonic signal of each moment transmitting of launch point is time-varying;Second group receive is put received ultrasonic signal and decomposes the characteristic vector that acquisition receives;The relatively characteristic vector of the ultrasonic signal that particular moment launches and the degree of association receiving the characteristic vector that point receives;When the two degree of association is more than threshold value, it is determined that receive the reflection echo that the ultrasonic signal that receives of point comprises the ultrasonic signal that this particular moment launches, and this reception point is detected first reflection echo when walking as described correction walk time.
6. the ray traveltime bearing calibration in ultrasound computed tomography according to claim 1 detection, it is characterised in that the revised ray traveltime of each ray is substituted into following Matrix Formula:
In formula, t1’-tn' for revising the ray traveltime after just;ajiRepresent this i-th ray length of process in jth image-generating unit;Concrete body A is divided m image-generating unit, and ultrasonic signal is approximately fixed value f in the slowness of each image-generating unit1-fm;Utilize the iterative algorithm such as ART realize based on revised walk time the numerical value inverting to wave speed distribution.
7. the ray traveltime bearing calibration in ultrasound computed tomography according to claim 1 detection, it is characterised in that for arbitrary launch point S1 ultrasonic signal launched, metering connects this launch point and first group n the ray traveltime t receiving point1-tn;
For this launch point S1 ultrasonic signal launched, measure this launch point to second group n and receive t when the correction put is walkedC1-tCn;
The thickness HB of serving can be obtained according to engineering construction drawing or in-site measurement, and each spacing receiving point and this launch point according to second group, obtain reflection echo ray, and then draw each reflection echo ray and from the ray angle α between first group of ray respectively received a little of this launch point1-αn;
Calculate the modifying factor β as this ray angle linear function1-βn;
Based on modifying factor, it is possible to t when correction is walkedC1-tCnIt is converted into ray traveltime correction, it may be assumed that
Δt1=β1·tC1
Δt2=β2·tC2
…
Δtn=βn·tCn
Utilize above ray traveltime correction, this launch point to first group is respectively received whole n ray traveltime t a little1-tnApply to revise, revised ray traveltime t1’-tn' as follows:
t1'=t1-Δt1
t2'=t2-Δt2
…
tn'=tn-Δtn
Revised ray traveltime t1’-tn' as ultrasonic signal within concrete body A actual walk time.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610127591.8A CN105806949B (en) | 2016-03-07 | 2016-03-07 | A kind of ray traveltime bearing calibration in ultrasound computed tomography detection |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610127591.8A CN105806949B (en) | 2016-03-07 | 2016-03-07 | A kind of ray traveltime bearing calibration in ultrasound computed tomography detection |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105806949A true CN105806949A (en) | 2016-07-27 |
CN105806949B CN105806949B (en) | 2019-04-02 |
Family
ID=56467488
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610127591.8A Active CN105806949B (en) | 2016-03-07 | 2016-03-07 | A kind of ray traveltime bearing calibration in ultrasound computed tomography detection |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105806949B (en) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108732245A (en) * | 2018-05-29 | 2018-11-02 | 辽宁省交通规划设计院有限责任公司 | A method of it is influenced for removal bituminous paving in treatment rates intensity detection |
CN108872388A (en) * | 2018-06-21 | 2018-11-23 | 长江勘测规划设计研究有限责任公司 | A kind of Ultrasonic CT detection system |
CN110927246A (en) * | 2018-09-20 | 2020-03-27 | 广州市市政工程试验检测有限公司 | Sleeve grouting compactness detection method and device based on ultrasonic CT |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2002243850A (en) * | 2001-02-21 | 2002-08-28 | Kajima Corp | Continuous electromagnetic wave type tomography and apparatus thereof |
CN102680579A (en) * | 2012-06-01 | 2012-09-19 | 中国计量学院 | Concrete internal defect detection method based on acoustic emission sensor array |
EP2541243A1 (en) * | 2011-07-01 | 2013-01-02 | Nederlandse Organisatie voor toegepast- natuurwetenschappelijk onderzoek TNO | Non destructive testing apparatus and method using ultrasound imaging |
CN103941289A (en) * | 2014-05-09 | 2014-07-23 | 中国科学院地质与地球物理研究所 | Diffraction point imaging method and equipment |
CN104374830A (en) * | 2014-11-26 | 2015-02-25 | 南京信息工程大学 | Near field phased array structure health monitoring method based on piezoelectric array |
CN104374828A (en) * | 2014-11-18 | 2015-02-25 | 上海岩土工程勘察设计研究院有限公司 | Ultrasonic tomography imaging method of detection on hidden defect |
CN104698086A (en) * | 2015-03-13 | 2015-06-10 | 哈尔滨工程大学 | Nondestructive flaw detection method for one-dimensional stress waves based on time inversion and multipath effects |
-
2016
- 2016-03-07 CN CN201610127591.8A patent/CN105806949B/en active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2002243850A (en) * | 2001-02-21 | 2002-08-28 | Kajima Corp | Continuous electromagnetic wave type tomography and apparatus thereof |
EP2541243A1 (en) * | 2011-07-01 | 2013-01-02 | Nederlandse Organisatie voor toegepast- natuurwetenschappelijk onderzoek TNO | Non destructive testing apparatus and method using ultrasound imaging |
CN102680579A (en) * | 2012-06-01 | 2012-09-19 | 中国计量学院 | Concrete internal defect detection method based on acoustic emission sensor array |
CN103941289A (en) * | 2014-05-09 | 2014-07-23 | 中国科学院地质与地球物理研究所 | Diffraction point imaging method and equipment |
CN104374828A (en) * | 2014-11-18 | 2015-02-25 | 上海岩土工程勘察设计研究院有限公司 | Ultrasonic tomography imaging method of detection on hidden defect |
CN104374830A (en) * | 2014-11-26 | 2015-02-25 | 南京信息工程大学 | Near field phased array structure health monitoring method based on piezoelectric array |
CN104698086A (en) * | 2015-03-13 | 2015-06-10 | 哈尔滨工程大学 | Nondestructive flaw detection method for one-dimensional stress waves based on time inversion and multipath effects |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108732245A (en) * | 2018-05-29 | 2018-11-02 | 辽宁省交通规划设计院有限责任公司 | A method of it is influenced for removal bituminous paving in treatment rates intensity detection |
CN108732245B (en) * | 2018-05-29 | 2020-10-16 | 辽宁省交通规划设计院有限责任公司 | Method for removing influence of asphalt pavement in service bridge strength detection |
CN108872388A (en) * | 2018-06-21 | 2018-11-23 | 长江勘测规划设计研究有限责任公司 | A kind of Ultrasonic CT detection system |
CN108872388B (en) * | 2018-06-21 | 2021-06-22 | 长江勘测规划设计研究有限责任公司 | Ultrasonic CT detection system |
CN110927246A (en) * | 2018-09-20 | 2020-03-27 | 广州市市政工程试验检测有限公司 | Sleeve grouting compactness detection method and device based on ultrasonic CT |
Also Published As
Publication number | Publication date |
---|---|
CN105806949B (en) | 2019-04-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Hua et al. | Time-frequency damage index of Broadband Lamb wave for corrosion inspection | |
Krause et al. | Elastic wave modes for the assessment of structural timber: ultrasonic echo for building elements and guided waves for pole and pile structures | |
Lin et al. | Dijkstra’s algorithm-based ray tracing method for total focusing method imaging of CFRP laminates | |
CN110346453B (en) | Method for rapidly detecting reflection echoes of small defect arrays in concrete structure | |
US8770027B2 (en) | Pulse-echo method by means of an array-type probe and temperature compensation | |
CN104142195A (en) | Device and method for detecting interior initial stress of steel structural member based on ultrasonic method | |
CN105699492A (en) | An ultrasonographic method used for weld seam detection | |
Michaels | Ultrasonic wavefield imaging: Research tool or emerging NDE method? | |
Zielińska et al. | Detection of debonding in reinforced concrete beams using ultrasonic transmission tomography and hybrid ray tracing technique | |
CN104698089A (en) | Ultrasonic relative time propagation technology suitable for inclined crack quantifying and imaging | |
CN106018548A (en) | Ultrasonic online detection method and system for uniformity of solid-liquid two-phase mixture | |
CN105806949A (en) | Correcting method for radial travel time in ultrasonic CT detection | |
CN106596725A (en) | Method for ultrasonic distinguishing of R-region defect of composite material structure | |
CN104142326A (en) | Attenuation coefficient detection method | |
Xu et al. | Quantification of guided mode propagation in fractured long bones | |
Hall | Air-coupled ultrasonic tomographic imaging of concrete elements | |
US10386335B2 (en) | Method for processing signals from an ultrasound probe acquisition, corresponding computer program and ultrasound probe device | |
Hoegh | Ultrasonic linear array evaluation of concrete pavements | |
CN104897778A (en) | Section grid ray average acoustic velocity method for detecting internal defect of concrete | |
Castaings et al. | Air-coupled measurement of plane wave, ultrasonic plate transmission for characterising anisotropic, viscoelastic materials | |
Cawley | Guided waves in long range nondestructive testing and structural health monitoring: Principles, history of applications and prospects | |
Cao et al. | In situ elastic constant determination of unidirectional CFRP composites via backwall reflected multi-mode ultrasonic bulk waves using a linear array probe | |
Belanger et al. | Lamb wave tomography to evaluate the maximum depth of corrosion patches | |
CN113433214B (en) | Composite plate Lamb wave time reversal multi-path imaging damage positioning method | |
Khazanovich et al. | Quantitative ultrasonic evaluation of concrete structures using one-sided access |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |