ULTRASONIC MEASUREMENT OF SURFACE PROFILE AND AVERAGE DIAMETER OF A TUBE

Abstract

Disclosed is a system and method for ultrasonic measurement of the average diameter and surface profile of a tube. A calibration block is used to calibrate the average tube diameter, and a correction is applied to account for any temperature difference of the couplant between calibration and test measurements. By using a linear probe, or a single probe with a finely pitched helicoidal scan, errors in diameter measurement due to presence of surface pits may be compensated.

Claims

1. A method of measuring a tube in a test mode, the tube having a tube outer surface, the method comprising the steps of: providing a measurement assembly configured to be mounted on the tube outer surface and to hold at least one ultrasonic probe, the probe is configured to emit an ultrasonic beam in a beam plane and to receive response signals from the tube outer surface, the probe being coupled to the tube outer surface with a test couplant made of a couplant material; causing relative rotation between the tube and the probe, the rotation being about a tube axis oriented along a tube axial direction, wherein an angle is an angle of rotation, and wherein the angle is between 0 and 360; measuring ultrasonic time of flight measurements, TOF(k), between the probe and the tube outer surface, wherein k is a sampling number and =k , wherein is an angular rotation increment between each ultrasonic time of flight measurement; converting the ultrasonic time of flight measurements to distance measurements, q(k), using a known test couplant sound velocity; calculating a Fourier transform of q(k); deriving harmonics, S(n), where n=0, 2, 3, . . . N, where S(n) is a n.sup.th harmonic of a tube outer surface profile, N is a highest harmonic number and n=0 corresponds to an average circle having an average diameter of the tube outer surface; calculating an inverse Fourier transform of S(n) to derive the tube outer surface profile; measuring a test couplant temperature; measuring a calibration time of flight in a calibration mode wherein the probe is coupled to a calibration block with a calibration couplant made of the couplant material and having a calibration couplant temperature; determining the average diameter based on the calibration couplant temperature and the test couplant temperature; and, calculating a roundness profile, wherein the roundness profile is a deviation of the tube outer surface profile from the average circle, the deviation being measured as a function of the angle .

2. The method of claim 1 wherein the probe is supported by a first support frame connected to a first contact bar and a second support frame connected to a second contact bar, wherein the first contact bar and the second contact bar are in contact with the tube outer surface.

3. The method of claim 2 wherein the calibration mode comprises the steps of: acoustically coupling the probe to the calibration block with the calibration couplant, wherein the probe is supported by the first support frame and the second support frame, and wherein the first contact bar and the second contact bar are in contact with a calibration block surface such that a calibration beam path between the probe and the calibration block surface corresponds to a nominal tube diameter; measuring the calibration time of flight in the calibration couplant from the probe to the calibration block surface and back to the probe; measuring the calibration couplant temperature; and, determining a calibration couplant sound velocity from a known couplant material sound velocity temperature dependence.

4. The method of claim 3 wherein the beam plane intersects the tube axis, wherein the first contact bar and the second contact bar are in contact with the tube outer surface at a first contact line and a second contact line respectively, wherein a first plane joining the first contact line to the tube axis makes an angle with the beam plane, and wherein a second plane joining the second contact line to the tube axis makes an angle with the beam plane.

5. The method of claim 4 wherein the average diameter is calculated in accordance with the equation $D_1=\left({\mathrm{TOF}}_1{V}_1-{\mathrm{TOF}}_0{V}_0\right).Math.\frac{\cos .Math..Math.}{1-\cos .Math..Math.}+D_0$ wherein D.sub.1 is the average diameter, D.sub.0 is the nominal tube diameter, TOF.sub.0 is the calibration time of flight, TOF.sub.1 is an average time of flight, wherein the average time of flight is the average of the ultrasonic time of flight measurements, V.sub.0 is the calibration couplant sound velocity and V.sub.1 is the test couplant sound velocity at the test couplant temperature.

6. The method of claim 2 wherein the tube outer surface has at least one surface pit thereon, the surface pit having a pit axial length, wherein the first contact bar has a bar axial length and the second contact bar has the bar axial length, and wherein the bar axial length is greater than the pit axial length and the linear probe length is greater than the pit axial length.

7. The method of claim 6 wherein the step of causing relative rotation includes simultaneously translating the tube in the axial direction thereby forming a helicoidal tube scan.

8. The method of claim 7 wherein the ultrasonic probe is a single probe and the helicoidal tube scan has a helicoidal pitch which is less than half the pit axial length.

9. The method of claim 6 wherein the ultrasonic probe is a linear probe having a plurality of elements configured as a plurality of transmission apertures and having a linear probe direction with a linear probe length oriented substantially parallel to the axial direction, wherein the plurality of transmission apertures measures a plurality of aperture times of flight between each aperture and the tube outer surface, and wherein there is a minimum aperture time of flight corresponding to each rotation angle .

10. The method of claim 9 wherein each of the plurality of aperture times of flight is set equal to the minimum aperture time of flight corresponding to each rotation angle , and wherein a reconstructed surface profile is thereby derived, wherein the reconstructed surface profile is a profile of the tube outer surface in the absence of the at least one pit and a reconstructed tube diameter is the average tube diameter in the absence of the at least one pit.

11. The method of claim 10 wherein a true surface profile at an axial location is calculated from the reconstructed surface profile, the plurality of aperture times of flight and the minimum aperture time of flight, wherein the true surface profile is an angular profile of the tube outer surface at the axial location, including the at least one pit.

12. The method of claim 11 wherein the true surface profile at the axial location is calculated in accordance with the equation
tsp()=rsp()V.sub.1.Math.[t()t.sub.min()] wherein tsp() is a radius of the true surface profile, rsp() is a radius of the reconstructed surface profile, V.sub.1 is the test couplant sound velocity, t() is a one of the plurality of aperture times of flight corresponding to an aperture at the axial location, and t.sub.min() is the minimum aperture time of flight.

13. The method of claim 1 wherein the couplant material is water.

14. A measurement apparatus for measuring a tube having a tube outer surface, the apparatus comprising, at least one ultrasonic probe, the probe is configured to emit an ultrasonic beam in a beam plane and to receive response signals from the tube outer surface, the probe being coupled to the tube outer surface with a test couplant made of a couplant material; a measurement assembly configured to be mounted on the tube outer surface, to hold the probe and to cause a relative rotation between the probe and the tube outer surface, the rotation being about a tube axis along a tube axial direction, wherein an angle is an angle of rotation, and wherein the angle is between 0 and 360; a test acquisition unit configured to measure ultrasonic time of flight measurements, TOF(k), between the probe and the tube outer surface, wherein k is a sampling number and =k , wherein is an angular rotation increment between each ultrasonic time of flight measurement, the test acquisition unit further configured to convert the ultrasonic time of flight measurements to distance measurements, q(k), using a known test couplant sound velocity; a Fourier transform calculation module configured to calculate a Fourier transform of q(k) and to derive harmonics, S(n), where n=0, 2, 3, . . . N, where S(n) is a n.sup.th harmonic of a tube outer surface profile, N is a highest harmonic number, and n=0 corresponds to an average circle having an average diameter of the tube outer surface; an inverse Fourier transform calculation module configured to calculate an inverse Fourier transform of S(n) to derive the tube outer surface profile; a temperature compensator configured to measure a test couplant temperature; a calibration system configured to measure a calibration time of flight wherein the probe is coupled to a calibration block with a calibration couplant made of the couplant material and having a calibration couplant temperature; an average diameter calibrator configured to calculate the average diameter based on the calibration couplant temperature and the test couplant temperature; and, a surface profile calculator configured to calculate a roundness profile, wherein the roundness profile is a deviation of the tube outer surface profile from the average circle, the deviation being measured as a function of .

15. The measurement apparatus of claim 14 wherein the Fourier transform is a fast Fourier transform and the inverse Fourier transform is an inverse fast Fourier transform.

16. The measurement apparatus of claim 14, wherein the tube is configured to rotate about the tube axis, and the measurement assembly is configured to hold the probe still.

17. The measurement apparatus of claim 14 wherein the measurement assembly further comprises a first support frame connected to a first contact bar and a second support frame connected to a second contact bar, wherein the first contact bar and the second contact bar are in contact with the tube outer surface at a first contact line and a second contact line respectively, wherein a first plane joining the first contact line to the tube axis makes an angle with the beam plane, and wherein a second plane joining the second contact line to the tube axis makes an angle with the beam plane.

18. The measurement apparatus of claim 17 further comprising a calibration module, the calibration module comprising: a calibration block, wherein during a calibration mode, the probe is coupled to the calibration block with the calibration couplant, wherein the probe is supported by the first support frame and the second support frame, and wherein the first contact bar and the second contact bar are in contact with a calibration block surface such that a calibration beam path in the beam plane between the probe and the calibration block surface corresponds to a nominal tube diameter; a calibration acquisition unit configured to measure the calibration time of flight from the probe to the calibration block surface and back to the probe; and, the temperature compensator configured to measure the calibration couplant temperature and to determine a calibration couplant sound velocity from a known couplant material sound velocity temperature dependence.

19. The measurement apparatus of claim 18 further comprising an average diameter calibrator configured to calculate a calibrated average diameter in accordance with the equation $D_1=\left({\mathrm{TOF}}_1{V}_1-{\mathrm{TOF}}_0{V}_0\right).Math.\frac{\cos .Math..Math.}{1-\cos .Math..Math.}+D_0$ wherein D.sub.1 is the average diameter, D.sub.0 is the nominal tube diameter, TOF.sub.0 is the calibration time of flight, TOF.sub.1 is an average time of flight, wherein the average time of flight is the average of the ultrasonic time of flight measurements, V.sub.0 is the calibration couplant sound velocity and V.sub.1 is the test couplant sound velocity at the test couplant temperature.

20. The measurement apparatus of claim 14 wherein the tube is translated in the axial direction and the roundness profile and the average diameter are calculated for a series of axial positions.

21. The measurement apparatus of claim 20 wherein the probe is a linear probe having a plurality of elements configured as a plurality of transmission apertures and having a linear probe direction with a linear probe length oriented substantially parallel to the axial direction, wherein the plurality of transmission apertures measures a plurality of aperture times of flight between each aperture and the tube outer surface, and wherein there is a minimum aperture time of flight corresponding to each rotation angle .

22. The measurement apparatus of claim 21 wherein the tube outer surface has at least one surface pit thereon, the surface pit having a pit axial length, wherein the first contact bar has a bar axial length and the second contact bar has the bar axial length, and wherein the bar axial length is greater than the pit axial length and the linear probe length is greater than the pit axial length.

23. The measurement apparatus of claim 22 wherein the test acquisition unit sets each of the plurality of aperture times of flight to be equal to the minimum aperture time of flight corresponding to each rotation angle , and wherein a reconstructed surface profile is thereby derived, wherein the reconstructed surface profile is a profile of the tube outer surface in the absence of the at least one pit and a reconstructed tube diameter is the average tube diameter in the absence of the at least one pit.

24. The measurement apparatus of claim 23 wherein a true surface profile calculator calculates a true surface profile at an axial location from the reconstructed surface profile, the plurality of aperture times of flight and the minimum aperture time of flight, wherein the true surface profile is an angular profile of the tube outer surface at the axial location, including the at least one pit.

25. The measurement apparatus of claim 24 wherein the true surface profile calculator calculates the true surface profile at the axial location in accordance with the equation
tsp()=rsp()V.sub.i.Math.[t()t.sub.min()] wherein tsp() is a radius of the true surface profile, rsp() is a radius of the reconstructed surface profile, V.sub.1 is the test couplant sound velocity, t() is a one of the plurality of aperture times of flight corresponding to an aperture at the axial location, and t.sub.min() is the minimum aperture time of flight.

26. The measurement apparatus of claim 14 wherein the couplant material is water.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

[0020] FIG. 1A is a schematic representation of an ovality measurement.

[0021] FIG. 1B is a schematic representation of a roundness measurement.

[0022] FIG. 2A is a schematic representation of a test measurement system for roundness profiling of a tube according to the present disclosure.

[0023] FIG. 2B is a schematic representation of parameters used for determination of the average diameter of a tube.

[0024] FIG. 3 shows a phase representation of the amplitude of undulations for the first 4 harmonics, n=0, 1, 2, and 3.

[0025] FIG. 4 shows a phase representation of the amplitude of undulations for the first 3 harmonics, n=0, 1, and 2, together with circular geometric representations of the tube profiles corresponding to each harmonic.

[0026] FIG. 5A is a diagram of a calibration block according to the present disclosure.

[0027] FIG. 5B is a schematic representation of a calibration system according to the present disclosure.

[0028] FIG. 6A is a graph showing the temperature dependence of sound velocity in water.

[0029] FIG. 6B is a graph showing the dependence of measured tube radius on the water temperature difference between calibration and measurement.

[0030] FIG. 7 is a diagram illustrating the principle of compensation of the measurement of the average diameter.

[0031] FIG. 8 is a photograph of a representative pit on a tube surface.

[0032] FIG. 9 is a schematic representation of a measurement system showing the effect of a pit in the tube surface.

[0033] FIG. 10 is a side view of a measurement system in the presence of a pit.

[0034] FIG. 11 is an illustration of geometry for measuring the depth of a pit.

[0035] FIG. 12 is a flowchart of a method of measuring a roundness profile and an average diameter of a tube according to the present disclosure.

[0036] FIG. 13 is a schematic representation of a measurement system for measuring a roundness profile and an average diameter of a tube according to the present disclosure.

DETAILED DESCRIPTION OF PREFERRED EMBODIMENT

[0037] FIG. 2A is a schematic representation of a test measurement system 1 according to the present disclosure, wherein test measurement system 1 is used for measurement of roundness and average diameter of a tube 4 having a tube axis 18. Note that the deviations from roundness of tube 4 illustrated in FIG. 2 are not to scale, but have been exaggerated for the purpose of presentation. An ultrasonic probe 10 emits an ultrasonic beam 12 which measures time of flight (TOF) in a couplant fluid 16 between ultrasonic probe 10 and the outer surface of tube 4. Ultrasonic probe 10 may be a single probe or may be a linear probe array having the direction of the linear array substantially parallel to the axis of tube 4. In an embodiment, probe 10 may be a linear phased array probe having 168 elements and an active length of 139 mm, but other probe types are possible and all are within the scope of the present disclosure. Ultrasonic beam 12 may therefore be a single beam or a linearly scanned beam with the scanning direction substantially parallel to axis 18. Support frames 6a and 6b support ultrasonic probe 10 by contacting tube 4 at lines of contact 8a and 8b respectively, wherein lines of contact 8a and 8b are substantially parallel to axis 18 and are defined by the contact of contact bars 9a and 9b respectively with the surface of tube 4. Lines of contact 8a and 8b each subtend an angle with respect to the direction of ultrasonic beam 12. Tube 4 is rotated about axis 18 as illustrated by an arrow 14, the rotation being effected by rollers 7a and 7b, and the rotation of tube 4 being represented by a rotation angle which varies from 0 to 2 during a single revolution of tube 4. Couplant fluid 16, which in an embodiment may be water, is retained in an enclosed space between probe 10, support frames 6a and 6b and tube 4. With knowledge of the sound velocity in couplant fluid 16, a time of flight measurement between ultrasonic probe 10 and the outer surface of tube 4 may be converted to a displacement measurement m(). The roundness profile of tube 4 is represented by s(), which is the dependence on angle of the deviation between the outer surface of tube 4 and an average circle 2.

[0038] FIG. 2B is a schematic representation of parameters used for measurement of average circle 2, wherein average circle 2 is a circle having a diameter D.sub.1, which is the average diameter of tube 4. A line 106a is tangential at a line of contact 8a on average circle 2 and a line 106b is tangential at a line of contact 8b on average circle 2, wherein lines of contact 8a and 8b both subtend angle with respect to the direction of ultrasonic beam 12. Note that, under the assumption that roundness deviations between tube 4 and average circle 2 are small, lines of contact 8a and 8b are close to, but not exactly equivalent to, lines of contact 8a and 8b shown in FIG. 2A. Lines 106a and 106b intersect at an intersection 19, and a distance h is the distance between intersection 19 and the active surface of probe 10. A displacement q() is the sum of displacement measurement m() and distance h, and corresponds to the distance from intersection 19 to the surface of tube 4.

[0039] In order to determine diameter Di, displacement q() must be determined for all values of angle from 0 to 2. It should be noted that distance h has no effect on the roundness measurement, but it does influence measurement of average circle diameter D.sub.1. The value of distance h is not precisely known in practice. However by using a calibration block having the nominal tube diameter, an accurate measurement of average circle diameter D.sub.1 may be made. The principle of the calibration is disclosed below with reference to equations (7), (8) and (9).

[0040] Displacement measurement q() can be related with roundness profile s() and the weighted roundness profile at the position of the two points of contact, s() and s(+), using the following equation:

[00001]$\begin{array}{cc}q\left(\right)=s\left(\right)-\frac{s\left(-\right)+s\left(+\right)}{2.Math..Math.\cos .Math..Math.}& \left(1\right)\end{array}$

[0041] The functions q() and s() can be represented by a sum of Fourier components Q(n) and S(n), where n takes integer values known as harmonics, each harmonic representing the number of undulations per revolution of tube 4. FIG. 3 is a phase representation of the normalized amplitude of undulations for the first four harmonics, n=0, 1, 2 and 3. For simplicity in each case the initial phase is set to 0.

[0042] FIG. 4 shows the phase representation of the first three harmonics, n=0, 1, and 2, together with circular geometric representations of the tube profiles corresponding to each harmonic. In the circular representations, the circles with solid lines are representations of average circle 2, while figures with broken lines represent profiles of actual tube 4.

[0043] The first harmonic (n=0) corresponds to a perfect circle for which the amplitude is a constant value and the deviation s() from roundness is zero all around the tube. The second harmonic (n=1) corresponds to a perfect circle which is not centered with respect to the axis of rotation. The dotted circle in the circular representation of harmonic n=1 shows an offset of tube 4 in the Y axis. The third harmonic (n=2) corresponds to tube 4 having an elliptical (oval) cross-section. The dotted ellipse in the circular representation of harmonic n=2 is centered with respect to the average circle 2. Higher harmonics (n>2) (not shown in FIG. 4) represent higher order variations in s(), representing roundness deviations which are non-oval.

[0044] Using the Fourier components Q(n) and S(n), for n=0, 2, . . . N, where N is the highest harmonic number, a harmonic sensitivity G(n) can be defined as:

[00002]$\begin{array}{cc}G\left(n\right)=\frac{Q\left(n\right)}{S\left(n\right)}& \left(2\right)\end{array}$

where G(n) is the sensitivity of probe measurement Q(n) to the tube profile variation S(n) for the n.sup.th harmonic. G(n) is given by:

[00003]$\begin{array}{cc}G\left(n\right)=1-\frac{\cos .Math..Math.n.Math..Math.}{\cos .Math..Math.}& \left(3\right)\end{array}$

Therefore:

[0045] [00004]$\begin{array}{cc}S\left(n\right)=\frac{Q\left(n\right)}{G\left(n\right)}=\frac{Q\left(n\right)}{1-\frac{\cos \left(n.Math..Math.\right)}{\cos .Math..Math.}}& \left(4\right)\end{array}$

Note that since Q(n) is generally a complex number, S(n) is a complex number in phase with Q(n).

[0046] In order to calculate Q(n), a series of ultrasonic measurements q(k) is made as tube 4 is rotated ( varies from 0 to 2), where k is the sampling number and =k , where is the angular increment between samples. Q(n) is then derived as

Q(n)=FFT{q(k)}(5)

for k=1, 2, . . . K, where K is the number of samples in one revolution of the tube, and where FFT{q(k)} is the Fast Fourier Transform (FFT) of the measurements q(k). The harmonics S(n) may then be calculated from equations (4) and (5), and finally the roundness profile of tube 4 is obtained using the Inverse Fast Fourier Transform of all harmonics S(n):

s(k)=FFT.sup.1{S(n)}(6)

where FFT .sup.1{S(n)} is the Inverse Fast Fourier Transform of S(n).

[0047] Note that not all harmonics can be used because those harmonics for which G(n) is very small may make the calculation unstable in the presence of signal noise. In particular, the second harmonic (n=1) cannot be used because G(1) is zero (see equation (3)). However, the second harmonic represents the displacement of tube 4 from the center of the nominal circle, which is not important for determination of roundness and average diameter. To avoid the problem of zero sensitivity for the second harmonic, the sensitivity

[00005]$\frac{1}{G\left(1\right)}=0$

for n=1 is set to for the Fourier Transform calculation.

[0048] In general, in order to obtain a precise profile of the outside surface of tube 4, it is necessary to have a large number of samples, meaning a small angular increment between samples. It is also necessary to avoid harmonics n for which sensitivity factor G(n) is small or zero. This may be achieved by careful selection of the angle subtended by lines of contact 8a and 8b with respect to the direction of ultrasonic beam 12. In an embodiment, may be equal to 48. In a second embodiment, may be equal to 42. However any other value of a may be advantageous, and all such values are within the scope of the present disclosure.

[0049] In order to avoid instabilities due to harmonics having small or zero sensitivity factor G(n), 1/G(n) may be set to zero for the Fourier Transform calculation for all values of G(n) less than a lower sensitivity limit. In an embodiment, the lower sensitivity limit may be 0.3, but other values of the lower sensitivity limit are possible, and all such values are within the scope of the present disclosure.

[0050] The Fourier Transform calculation may be performed for harmonics n=0, 1, 2, . . . N, where N is the highest order harmonic allowed by the system design. For example, if G(n)>0.3, then according to Eq. 3, 1-cos(n)/cos()>0.3. This relationship may be satisfied for some combinations of and N. In general, the value of N should be as large as possible, and in any event should be greater than or equal to 6. Another design factor to consider is the size of the measuring mechanism. The larger the value of angle , the larger the size of the measuring mechanism. Therefore there is a compromise between the size of the mechanism and the highest achievable harmonic number N. In an embodiment, with =48, the highest order harmonic may be N=13, but other values of the highest order harmonic are possible, and all such values are within the scope of the present disclosure.

[0051] It is important to note that the distances, S(n), Q(n), s(k) and q(k), and the angle , are all referenced to the diameter D.sub.1 of average circle 2. With the assumption that the deviations from average circle 2 are small, only the first harmonic n=0 corresponds to the diameter of average circle 2. However, the diameter of average circle 2 cannot be accurately measured with test measurement system 1 because the exact length of the couplant column between probe 10 and the surface of tube 4 is difficult to measure, and the distance h (see FIG. 2B) is unknown. The inventors of the present disclosure have discovered that the accuracy of the average diameter measurement can be significantly improved by performing a calibration with a calibration system 3 using a calibration block 40 as shown in FIGS. 5A and 5B. Calibration block 40 is configured to allow stable and repeatable location on calibration block 40 of the measurement system comprising probe 10, support frames 6a and 6b and contact bars 9a and 9b. Calibration block 40 is further configured so that when the measurement system is located on calibration block 40, the acoustic path length from probe 10 to the surface of calibration block 40 accurately corresponds to a length L, wherein length L is the acoustic path length from probe 10 to a nominal circle 2a having a nominal diameter D.sub.0 (see FIG. 5B). Note that nominal circle 2a, as defined by the geometry of calibration system 3, represents the nominal diameter of tube 4. However, nominal diameter D.sub.0 may be different from average diameter D.sub.1, which is the average tube diameter as measured with test measurement system 1. As illustrated in FIG. 5B, calibration block 40 has a cross-section which is part of nominal circle 2a. However, those skilled in the art may devise various forms and constructions of calibration block 40 in which the acoustic path length from probe 10 to the surface of calibration block 40 accurately corresponds to the length L, and all such forms and constructions are within the scope of the present disclosure.

[0052] Calibration according to the present disclosure is performed by measuring a calibration time of flight in couplant fluid 16 when probe 10 is coupled to calibration block 40. Note that calibration block 40 is not rotated during the calibration measurement. Probe 10 is then coupled to tube 4 which is rotated, and the measured average diameter D.sub.1 (first harmonic) is compensated according to the calibration time of flight. The calibration allows an accurate measurement of the first harmonic n=0 corresponding to the actual average diameter D.sub.1 of tube 4. The method of compensating the first harmonic according to the calibration time of flight is described below in connection with FIG. 7 and equations (7)-(11).

[0053] The inventors of the present disclosure have observed that compensation of the first harmonic may be inaccurate if the temperature of couplant fluid 16 during calibration with calibration system 3 differs from the temperature of couplant fluid 16 during the measurement with test measurement system 1. FIG. 6A shows a curve 28 representing the dependence of sound velocity in water on the temperature of the water. Curve 28 is well known in the art. FIG. 6B shows calculations of the radius of average circle 2 obtained from the calibrated first harmonic when measuring a tube having a nominal radius of 8.075 inches. A line 30 shows the calibrated radius of average circle 2 when the water temperature is 20 C. for both the calibration and the measurement. A curve 32 shows the calculated radius of average circle 2 as a function of the temperature difference between the calibration temperature at 20 C. and the measurement temperature. It may be seen from FIG. 6B that a 5 C. difference in water temperature causes an error of approximately 1% in the calculated radius. In the method of the present disclosure, this error is substantially eliminated by measuring the water temperature during calibration and during measurement, and using the velocity dependence of curve 28 (FIG. 6A) to correctly convert the time of flight measurements to displacement measurements (see equations (7), (8) and (9) below). Note that only the first harmonic is corrected for couplant temperature. Correction of higher harmonics is not required.

[0054] FIG. 7 illustrates the principle of compensation of the measurement of the first harmonic n=0 with calibration block 40. FIG. 7 shows a superposition of measurement parameters for calibration with calibration block 40 (having diameter D.sub.0 and radius R.sub.0) and subsequent measurement of tube 4 (having average diameter D.sub.1 and average radius R.sub.1). In FIG. 7, line 106a is tangential at line of contact 8a to both average circle 2 and calibration block 40, and a line 106b is tangential at line of contact 8b to both average circle 2 and calibration block 40. Note that axis 18 of average circle 2 may not coincide with the corresponding axis 18a of calibration block 40. Distances q.sub.0 and q.sub.1 are measured from intersection 19 to average circle 2 and calibration block 40 respectively. Therefore, the distance q.sub.0-h represents a sound path 12a from the active surface of probe 10 to the surface of calibration block 40, and the distance q.sub.1-h represents the sound path 12 from the active surface of probe 10 to the surface of tube 4. Sound paths 12 and 12a, and the difference q between them, can be calculated with equations (7), (8) and (9) below:

q.sub.0h=TOF.sub.0V.sub.0/2 (7)

q.sub.1h=TOF.sub.1V.sub.1/2 (8)

q=q.sub.1q.sub.0=(TOF.sub.1V.sub.1TOF.sub.0V.sub.0)/2 (9)

where TOF.sub.0 is the complete time of flight in the sound path q.sub.0-h, TOF.sub.1 is the complete time of flight in the sound path q.sub.1-h, wherein TOF.sub.1 is derived as the average time of flight measured during a complete rotation of tube 4, V.sub.0 is the sound velocity in couplant fluid 16 during calibration with calibration block 40, and V.sub.1 is the sound velocity in couplant fluid 16 during measurement of tube 4. Note that V.sub.0 and V.sub.1 may be different because the temperature of couplant fluid 16 may be different for the two measurements. The following equation may be derived from the geometrical relationships shown in FIG. 7:

[00006]$\begin{array}{cc}{D}_1=2.Math..Math..Math..Math.q.Math.\frac{\cos .Math..Math.}{1-\cos .Math..Math.}+D_0,& \left(10\right)\end{array}$

where D.sub.0 is the diameter of calibration block 40 (corresponding to the nominal diameter of tube 4) and D.sub.1 is the average diameter of tube 4 (corresponding to the first harmonic n=0 that needs to be compensated).

[0055] Note that in equations (9) and (10) the dependence on unknown distance h has been eliminated, allowing an accurate determination of average diameter D.sub.1.

[0056] From Equations (9) and (10), the following final expression is derived:

[00007]$\begin{array}{cc}{D}_1=\left({\mathrm{TOF}}_1{V}_1-{\mathrm{TOF}}_0{V}_0\right).Math.\frac{\cos .Math..Math.}{1-\cos .Math..Math.}+D_0& \left(11\right)\end{array}$

[0057] Equation (11) allows calibration of the measurement of average diameter D.sub.1 by using calibration block 40 and taking into account the effect of temperature change.

[0058] It should be noted that the equipment and methods of the present disclosure are generally applied to tubes having substantial axial length, and that the average diameter and roundness profile should therefore be measured as a function of axial position. Such measurements may be made by axially translating tube 4 during the measurement and/or by axially scanning ultrasonic beam 12 from ultrasonic probe 10, wherein probe 10 is configured as a linear array probe. In making such measurements, it is assumed that the average diameter and roundness profile of the tube change slowly in the axial direction relative to the axial lengths of ultrasonic probe 10 and contact bars 9a and 9b. While this assumption is generally true, errors can occur if there are small pits in the surface of tube 4. FIG. 8 is a photograph of a surface pit which is typically a result of the method of production of the tube. The dimension p representing the axial length of the pit is typically 1 to 10 mm. The depth of a pit is typically 1 to 3 mm.

[0059] FIG. 9 illustrates a situation where the outside surface of tube 4 has a pit 20 whose dimension in the axial direction is less than the axial length of probe 10 and contact bars 9a and 9b. As illustrated in FIG. 9, when tube 4 is rotated so that pit 20 is under ultrasonic beam 12, ultrasonic beam 12 may be sufficiently focused to probe the depth of pit 20. However, when tube 4 rotates so that pit 20 is under either contact bar 9a or contact bar 9b, the contact bar, being axially longer than pit 20, will ride over pit 20 and the contact bar radial position will be unaffected by pit 20. This situation violates the three point contact assumption inherent in equation (1), and therefore diameter D.sub.1 of average circle 2 will be overestimated at the axial location of pit 20.

[0060] FIG. 10 shows a side view of tube 4, ultrasonic probe 10, support frame 6a and contact bar 9a. Ultrasonic probe 10 is longer in the axial direction than pit 20, and emits multiple ultrasonic beams 12, some of which probe the depth of pit 20 and others do not. Among ultrasonic beams 12 there is an ultrasonic beam 12a which probes a maximum depth d of pit 20 and which therefore has a maximum time of flight t.sub.max between probe 10 and tube 4. Also among ultrasonic beams 12 there is an ultrasonic beam 12b which does not probe pit 20 and which has a minimum time of flight t.sub.min between probe 10 and tube 4. Note that under the assumption that the roundness profile varies slowly in the axial direction, there will be multiple beams 12b having substantially the same minimum time of flight, corresponding to all those beams 12 which do not probe pit 20. Furthermore, under the assumption that the roundness profile varies slowly in the axial direction, the time of flight of those beams 12, such as beam 12a, which do probe pit 20, may be replaced by minimum time of flight t.sub.min in order to obtain a reconstructed surface profile rsp() which represents the profile of tube 4 in the absence of pit 20.

[0061] In summary, roundness calculation errors due to the presence of pit 20 may be corrected by using a linear array probe 10 emitting multiple ultrasonic beams 12, determining a minimum time of flight of ultrasonic beams 12, and calculating displacement measurement q() based on the minimum time of flight for all axial positions along the length of probe 10.

[0062] As an alternative to using a linear probe array, a single probe or a single aperture of a probe array may be used in conjunction with a helicoidal scan, in which case the times of flight of many adjacent scans may be axially compared and the minimum time of flight used to calculate displacement measurement q(). The pitch of the helicoidal scan should be less than half the axial dimension of a typical pit in tube 4.

[0063] FIG. 11 illustrates how the depth profile of the pit may be determined as tube 4 is rotated. A true tube surface profile tsp() at each axial location, which takes the pit profile into account, may then be calculated from the equation:

tsp()=rsp()V.sub.1[t()t.sub.min()](12)

wherein tsp() is the radius of the true surface profile, rsp() is the radius of the reconstructed surface profile, V.sub.1 is the acoustic velocity in couplant 16, t() is the time of flight corresponding to an aperture at the axial location, and t.sub.min() is the minimum time of flight corresponding to beam 12b measured by linear probe 10 at angle .

[0064] As shown in FIG. 11, the maximum depth d of pit 20 will be measured at an angle .sub.p, wherein t(.sub.p)=t.sub.max. In equation (12), variation of V.sub.1 with temperature is not important since very accurate measurement of profile tsp() is not required.

[0065] Note that FIG. 10 illustrates the thickness D of tube 4, which may be measured as a function of angle and axial position by means of the time of flight difference between ultrasonic beams 12 reflected from the outer surface of tube 4 and those reflected from inner surface of tube 4 (not shown), and with knowledge of the acoustic velocity in the material of tube 4. Thus, use of the ultrasonic method has the advantage that both the roundness and the tube thickness may be measured.

[0066] FIG. 12 is a flowchart of a method of measuring the roundness profile and average diameter of a tube according to the present disclosure. Steps 50, 52 and 54 are the steps of a calibration mode 90. In step 50, probe 10 is coupled to calibration block 40 so that the acoustic path length in couplant 16 between probe 10 and calibration block 40 is equal to the nominal acoustic path length L between probe 10 and nominal circle 2a. In step 52, calibration time of flight TOF.sub.0 is measured between probe 10 and the surface of calibration block 40, and the temperature of couplant 16 is measured at the time of calibration. In step 54, the calibration time of flight and couplant temperature are stored for use in subsequent test measurements.

[0067] A test mode measurement 92 on a tube workpiece begins at step 56, in which probe 10 is coupled to tube 4. In step 57 the couplant temperature for the test mode measurement is measured and stored. In step 58, a series of ultrasonic measurements q(k) is made as tube 4 is rotated 360. The measurements are made as a function of axial position with linear probe 10 having its axis parallel to the axis of tube 4, and for each axial position q(k) is derived from the minimum time of flight train of all apertures of probe 10. In step 60, the n.sup.th harmonics Q(n) are derived from the Fast Fourier Transform of q(k), and the corresponding harmonics of the tube profile, S(n), are then calculated from equation (4). In step 62, the inverse FFT is applied to S(n) to derive the reconstructed surface profile s(k) of tube 4, which is the profile in the absence of surface pits.

[0068] Steps 64 and 66 are steps of calibration and of corrections applied to s(k) to account for temperature differences and the presence of pits. In step 64, the average tube diameter D.sub.1 is calculated and temperature compensated in accordance with equation (11), using the calibration time of flight TOF.sub.0, the average time of flight TOF.sub.1, and the different sound velocities at the test couplant temperature and the calibration couplant temperature. In step 66, the true surface profile, taking into account the presence of pits, is calculated from equation (12) based on the difference between the measured TOF at each angle and axial position, and the minimum TOF at that angle for all apertures of linear probe 10.

[0069] FIG. 8 is a schematic representation of a measurement system 100 for measuring the roundness profile and average diameter of a tube according to the present disclosure. Measurement system 100 comprises calibration system 3 with probe 10 coupled to calibration block 40. A calibration acquisition unit 77 acquires ultrasound signals from probe 10 and transmits calibration TOF.sub.0 to an average diameter calibrator 84. A temperature compensation unit 72 obtains a calibration couplant temperature measurement of couplant 16 during the calibration measurement, and outputs a calibration sound velocity V.sub.0 with input from a sound velocity vs temperature table 74.

[0070] Measurement system 100 further comprises test measurement system 1, which is deployed subsequent to the calibration measurement, and from which linear probe 10 transmits data to a test acquisition unit 76 and to a linear probe TOF unit 78. Linear probe TOF unit 78 determines times of flight for all apertures of linear probe 10 and transmits minimum TOF, t.sub.min(), to test acquisition unit 76. Temperature compensation unit 72 obtains a test couplant temperature measurement of couplant 16 during the test measurement, and outputs a test sound velocity V.sub.1 with input from sound velocity vs temperature table 74. Test acquisition unit 76 acquires a series of ultrasonic measurements q(k) as tube 4 is rotated 360, the values of q(k) being based on the minimum TOF, t.sub.min(), for all axial positions along the length of probe 10. A FFT calculation module 80 then performs a Fast Fourier Transform to obtain the harmonics S(n) of the profile of tube 4, and an Inverse FFT calculation module 82 performs an Inverse Fast Fourier Transform to obtain the reconstructed surface profile s(k) of tube 4, which is the profile in the absence of surface pits. With input from linear probe TOF unit 78, a true surface profile calculator 86 uses equation (12) to produce the true surface profile, which is the surface profile of tube 4 taking into account the presence of pits.

[0071] An average diameter calibrator 84 calibrates the average diameter D.sub.1 (the first harmonic of s(k)) by means of equation (11), using calibration TOF.sub.0 from calibration system 3 and calibration and test sound velocities V.sub.0 and V.sub.1 from temperature compensation unit 72.

[0072] Therefore the output from measurement system 100 is the true surface profile and calibrated average diameter D.sub.1.

[0073] Although the present invention has been described in relation to particular embodiments thereof, it can be appreciated that various designs can be conceived based on the teachings of the present disclosure, and all are within the scope of the present disclosure.