CN105806949A - Correcting method for radial travel time in ultrasonic CT detection - Google Patents

Correcting method for radial travel time in ultrasonic CT detection Download PDF

Info

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
Application number
CN201610127591.8A
Other languages
Chinese (zh)
Other versions
CN105806949B (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.)
Hangzhou Huaxin Testing Engineering Co Ltd
Original Assignee
Hangzhou Huaxin Testing Engineering Co Ltd
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 Hangzhou Huaxin Testing Engineering Co Ltd filed Critical Hangzhou Huaxin Testing Engineering Co Ltd
Priority to CN201610127591.8A priority Critical patent/CN105806949B/en
Publication of CN105806949A publication Critical patent/CN105806949A/en
Application granted granted Critical
Publication of CN105806949B publication Critical patent/CN105806949B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N29/00Investigating 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/04Analysing solids
    • G01N29/06Visualisation of the interior, e.g. acoustic microscopy
    • G01N29/0654Imaging
    • G01N29/0672Imaging by acoustic tomography
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N29/00Investigating 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/44Processing the detected response signal, e.g. electronic circuits specially adapted therefor
    • G01N29/4463Signal correction, e.g. distance amplitude correction [DAC], distance gain size [DGS], noise filtering
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2291/00Indexing codes associated with group G01N29/00
    • G01N2291/02Indexing codes associated with the analysed material
    • G01N2291/023Solids
    • G01N2291/0232Glass, ceramics, concrete or stone
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2291/00Indexing codes associated with group G01N29/00
    • G01N2291/02Indexing codes associated with the analysed material
    • G01N2291/028Material parameters
    • G01N2291/0289Internal 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

A kind of ray traveltime bearing calibration in ultrasound computed tomography detection
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 point1n
Calculate the modifying factor β as this ray angle linear function1n
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
Δt11·tC1
Δt22·tC2
Δtnn·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-Rn1n.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 function11)-βnn).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
Δt11·tC1
Δt22·tC2
Δtnn·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:
t 1 , t 2 , ... t N , = a 11 a 21 ... a m 1 a 12 a 22 ... a m 2 ... a 1 n a 2 n ... a m n f 1 f 2 ... f m
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 point1n
Calculate the modifying factor β as this ray angle linear function1n
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
Δt11·tC1
Δt22·tC2
Δtnn·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.
CN201610127591.8A 2016-03-07 2016-03-07 A kind of ray traveltime bearing calibration in ultrasound computed tomography detection Active CN105806949B (en)

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)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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

Patent Citations (7)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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