GB2545044A

GB2545044A - Determining a thickness of a region of wall-or plate-like structure

Info

Publication number: GB2545044A
Application number: GB1613807.5A
Authority: GB
Inventors: Pialucha Tomasz
Original assignee: Guided Ultrasonics Ltd
Current assignee: Guided Ultrasonics Ltd
Priority date: 2016-08-11
Filing date: 2016-08-11
Publication date: 2017-06-07
Anticipated expiration: 2036-08-11
Also published as: US11022436B2; CA3031163A1; GB201613807D0; JP2019525205A; GB2545044B; WO2018029445A1; EP3497402B1; JP7053616B2; GB2552858A; US20200173773A1; CN110088564B; CN110088564A; GB201702166D0; EP3497402A1

Abstract

A method of determining a thickness of a region of wall- or plate-like structure which is thinner than a surrounding region of the structure due to a defect in the structure. The method comprises comparing a measured time-frequency dispersion map (51), obtained by measuring the structure using at least one dispersive guided wave, with a reference time-frequency dispersion map (45, Fig. 19) obtained by modelling the structure. A split or cut-off frequency fc at which the measured and reference time-frequency dispersion maps differ is determined and the thickness of the corroded region is calculated in dependence upon the cut-off frequency. The structure may be modelled by testing the structure at a reference section which is substantially free from defects. The dispersive guided wave may be an antisymmetric or symmetric shear horizontal guided wave. The time-frequency dispersion map may be obtained for a set of different frequencies and a transducer excites the structure using a tone burst to cause fundamental symmetric and first excited mode antisymmetric shear horizontal guided waves to be generated. The tested structure may be a pipe wall. Also claimed is a transducer, comprising an array of spaced transducer elements wherein the transducer spacing is variably adjustable.

Description

Determining a thickness of a region of wall- or plate-like structure Field of the Invention

The present invention relates to determining a thickness of a region of a wall- or platelike structure, such as a wall, plate, pipe or a hollow bar.

Background

Ultrasonic guided waves can be used to inspect pipes, rails, rods, plates and other types of structures for corrosion-type defects, cracks and other types of flaws.

Guided waves can be used to screen structures over long ranges extending, for example, over tens or hundreds of metres. Moreover, guided waves can be used to inspect structures that are coated with insulation, which are difficult to access (for example, due to being buried or being underwater) and/or while they are in use.

An introduction to and an overview of the capabilities of long range ultrasonic guided wave testing systems can be found in M. Lowe and P. Cawley: “Long Range Guided Wave Inspection Usage - Current Commercial Capabilities and Research Directions” (2006) (http://www3.imperial.ac.Uk/pls/portallive/docs/1/55745699.PDF).

Guided waves can propagate in a structure in a number of different modes including torsional modes, and flexural and compressional modes. Torsional modes, in the form of shear horizontal waves, can be particularly useful for inspecting pipes while they are in use because little of the guided wave energy is coupled into the external media inside and outside the pipe.

Although shear horizontal waves can be used to measure average pipe thickness, shear horizontal waves have not been considered in industry to be particularly suited to detecting and determining the size of defects, such as reduced wall thickness arising from corrosion of an inner or outer diameter of the wall.

Lamb and shear horizontal waves near their cut-off frequencies to determine the presence pipe corrosion using transmission or reflection coefficients, or mode conversion phenomenon have been investigated in the laboratory.

For example, the use of mode cut off to identify and to obtain a thickness of a corroded region has been investigated experimentally as described in J. L. Rose & J. N. Barshinger: “Using Ultrasonic Guided Wave Mode Cutoff for Corrosion Detection and Classification”, IEEE ULTRASONICS SYMPOSIUM, page 851 (1998), W. Zhu, J. L. Rose, J. N. Barshinger, V. S. Agarwal: “Ultrasonic Guided Wave NDT for Hidden Corrosion Detection”, Research in Nondestructive Evaluation, pages 205 to 225 (1998) and D. Tuzzeo & F. Lanza di Scalea: “Noncontact Air-Coupled Guided Wave Ultrasonics for Detection of Thinning Defects in Aluminum Plates” volume 10, pages 61 to 77 (2001).

Mode conversion has also been observed experimentally as described in Nurmalia, N. Nakamura, H. Ogi, and M. Hirao: “Detection of Shear Horizontal Guided Waves Propagating in Aluminum Plate with Thinning Region”, Japanese Hournal of Applied Physics, volume 50, pages 07HC17-1 to 07HC17- 5 (2011), Nurmalia, N. Nakamura, H. Ogi, M. Hirao, K. Nakahata: “Mode conversion behaviour of SH guided wave in a tapered plate”, NDT&E International, volume 45, pages 156 to 161 (2012) and P. Belanger: “High order shear horizontal modes for minimum remnant thickness”, Ultrasonics, volume 54, pages 1079 to 1087 (2014).

Notwithstanding these laboratory investigations, these approaches have not been adopted in real-world guided waves inspection systems, not least because of a number of unknown variables which makes identification and characterisation of defects in a real-world structure (which can have any size of defects) more difficult.

Summary

According to a first aspect of the present invention there is provided a method of determining a thickness (“a remaining thickness”) of a region of a wall- or plate-like or -based structure (such as pipe, rectangular conduit or cylindrical vessel) which is thinner than a thickness of a surrounding region of the structure due to a cavity in the structure. The method comprises comparing a measured time-frequency dispersion map for at least one dispersive guided wave obtained by measuring the structure using guided waves with a reference time-frequency dispersion map obtained by modelling the cavity-free structure, determining a cut-off frequency at which the measured time-frequency dispersion map and the reference time-frequency dispersion map differ and calculating the remaining thickness in dependence upon the cut-off frequency.

The at least one dispersive guided wave is at least one shear horizontal guided wave of the higher order than one fundamental symmetric shear horizontal guided wave (SH0), for example, a first-excited mode antisymmetric shear horizontal guided wave (AHi).

The method may comprise obtaining the measured time-frequency dispersion map by, for a set of different frequencies, causing a transmitter transducer to excite the structure using a tone burst at a given frequency, receiving a time-varying signal from a receiver transducer and storing the time-varying signal.

The time-varying signal may include signals corresponding to one or more guided waves. The one or more guided waves may include at least one fundamental symmetric shear horizontal guided wave (SH0) and/or at least one a first-excited mode antisymmetric shear horizontal guided wave (AHi). The one or more guided waves may include at least one transmitted wave. The at least one transmitted wave may include periodically-detected wave (i.e. a wave which travels around a loop in a structure (such as pipe circumference) and which is detected each time the wave returns to a receiver transducer thereby resulting in two or more waves). The one or more guided waves may include at least one reflected wave.

Causing the transmitter transducer to excite the structure using the tone burst at the given frequency may cause a fundamental symmetric shear horizontal, SH0, guided wave and a first-excited mode antisymmetric shear horizontal, AHi, guided wave to be generated in the structure.

Modelling the structure may comprise generating the reference time-frequency dispersion map in dependence upon a thickness of the structure. Generating the reference time-frequency dispersion map may comprise calculating a time of flight as a function of frequency or vice versa dependent on guided wave mode number, a distance between a transmitter transducer and a receiver transducer, a thickness of a wall or plate of the structure and bulk shear velocity.

Modelling the structure may comprise, for a reference section of the structure (for example, a section of a pipe which is substantially free from defects or at least free from defects which are of interest), causing the transmitter transducer to excite the structure using a reference tone burst, receiving a reference time-varying signal from a receiver transducer and determining a distance between a transmitter transducer and a receiver transducer and a thickness of a wall or plate of the structure in dependence upon the reference time-varying signal.

The cavity is preferably a surface cavity, i.e. the cavity extends from a surface (which maybe an outer surface or an inner surface) of the structure. For example, the cavity may be at the surface of a pipe.

The structure is preferably elongate, i.e. has a length which is greater than its transverse dimension (e.g. diameter). The structure preferably comprises a wall or plate which has a substantially uniform thickness. The wall or plate may have a curvature in a transverse direction (i.e. around a pipe circumference) which is preferably smooth, i.e. there are no abrupt changes in thickness or in orientation.

The wall- or plate-like structure may be a pipe comprising a pipe wall.

The method may be implemented in software, e.g. in a (programmable) computer system comprising memory and at least one processor. The method may be implemented, at least in part, in hardware. For example, an application-specific integrated circuit (ASIC) or field-programmable gate array (FPGA) maybe used to carry out at least some of the steps of the method.

According to a second aspect of the present invention there is provided a computer program which, when executed by a computer system, causes the computer to perform the method.

According to a third aspect of the present invention there is provided a computer-readable medium, which may be a non-transitory computer-readable medium, storing the computer program.

According to a third aspect of the present invention is provided apparatus configured to compare a measured time-frequency dispersion map for at least one dispersive guided wave obtained by measuring the structure using guided waves with a reference time-frequency dispersion map obtained by modelling the cavity-free structure, to determine a cut-off frequency at which the measured time-frequency dispersion map and the reference time-frequency dispersion map differ and to calculating the thickness of the thinner region in dependence upon the cut-off frequency.

The apparatus may be a computer system comprising memory and at least one processor, and the at least one processor configured to compare a measured time-frequency dispersion map for at least one dispersive guided wave obtained by measuring the structure using guided waves with a reference time-frequency dispersion map obtained by modelling the structure, to determine a cut-off frequency at which the measured time-frequency dispersion map and the reference time-frequency dispersion map differ and to calculating the thickness of the thinner region in dependence upon the cut-off frequency.

For a set of different frequencies, the apparatus is configured to cause a transmitter transducer to excite the structure using a tone burst at a given frequency, to receive a time-vaiying signal from a receiver transducer and to store the time-varying signal so as to generate the measured time-frequency dispersion map.

The apparatus may be configured to cause a transducer to be mechanically adjusted so as to access a different range of wavelengths. The apparatus maybe configured to cause spacing of elements in the transducer to change. For example, the apparatus may be configured to change spacing of elements in the transducer by an actuator, such as electric motor.

According to a fourth aspect of the present invention is provided a guided wave inspection system comprising the apparatus (e.g. the computer system), a transmitter comprising a transmitter transducer for exciting guided wave modes in the structure and a receiver comprising a receiver transducer for obtaining a time-varying signal.

The transmitter and receiver transducers may comprise respective transducers.

The transmitter transducer and receive transducer may comprise an electromagnetic acoustic transducer (EMAT) or first and second electromagnetic acoustic transducers (EMATs).

The transducer(s) may (each) comprise an array of spaced-apart transducer elements which are arranged such that spacing between the transducer elements is variably adjustable. For example, the transducer may comprise at least one linear bearing and at least one motor, for example in the form of a stepper motor. One or more transducer elements may be mounted on or to the linear bearing.

According to a fifth aspect of the present invention is provided a transducer for guided wave inspection comprising an array of spaced-apart transducer elements which are arranged such that spacing between the transducer elements is variably adjustable.

The transducer may be an electromagnetic acoustic transducer (EMAT) and the transducer elements comprise permanent magnets.

The transducer is configured such that the spacing is variable so as to take a value between a lower limit and an upper limit. The lower limit may lie in the range l mm to to mm and the upper limit may lie in the range 10 mm to 20 mm.

Brief Description of the Drawings

Certain embodiments of the present invention will now be described, by way of example, with reference to the accompanying drawings, in which:

Figure 1 is a cross-sectional view of a pipe wall and a defect in the pipe wall;

Figure 2 is a schematic block diagram of a guided wave shear horizontal sensing system which can be used to screen a plate or pipe;

Figure 3 is a partial transverse cross-section of a steel pipe and illustrates an arrangement in which a transmitter and a receiver of a sensing system are placed at different positions around the circumference of a pipe next to the outside of pipe;

Figure 4 is schematic time-domain plot of a fundamental shear horizontal mode (SH0) tone burst which arrives at a receiver 100.0 ps after being generated by a transmitter in the arrangement shown in Figure 3;

Figure 5 is a partial transverse cross-section of a steel pipe and illustrates an arrangement in which a transmitter and a receiver of a sensing system are placed at different positions around the circumference of a pipe next to the outside of pipe similar to that shown in Figure 3;

Figure 6 is schematic time-domain plot of a dispersive, anti-symmetric shear horizontal mode (AHi) tone burst and a corresponding Fourier transformation of a windowed signal centred around 103.6 ps and illustrating a dominant frequency of the windowed signal as being 254.5 kHz;

Figure 7 illustrates a first transmitter-receiver arrangement, similar to that shown in Figure 5, in which a receiver is interposed between a transmitter and a pipe defect; Figure 8 is a schematic time-domain plot of a transmitted dispersive, anti-symmetric shear horizontal mode (AHi) and its reflection from a defect in the first arrangement shown in Figure 7;

Figure 9 is a schematic plot of a fast Fourier Transform of a reference AHi signal shown in Figure 8;

Figure 10 is a schematic plot of a fast Fourier Transform of a signal shown in Figure 8 reflected from a defect, where fc is a cut-off frequency of the reflected signal;

Figure 11 illustrates a second transmitter-receiver arrangement, similar to that shown in Figure 5, in which a defect is interposed between a transmitter and a receiver;

Figure 12 is a schematic time-domain plot of a transmitted dispersive AHi signal in the second arrangement shown in Figure 11;

Figure 13 is a schematic plot of a fast Fourier Transform of a signal shown in Figure 12 transmitted through the defect, where fc is a cut-off frequency of the reflected signal

Figure 14 illustrates a transmitter-receiver arrangement around the circumference of a pipe;

Figure 15 is a schematic plot of non-dispersive SH0 wavelets travelling around the circumference of a pipe in clockwise and anticlockwise directions;

Figure 16 is a schematic time-frequency plot of dispersive AHi mode travelling around the circumference of a pipe;

Figure 17a is a process flow diagram of a first part of method of determining depth of a defect;

Figure 17b is a process flow diagram of a second part of method of determining depth of a defect

Figure 18 is schematic time-domain plot of SH0 and AHi modes travelling around the circumference of a pipe without a defect;

Figure 19 is a schematic plot of frequency against time of the SHo and AHi modes travelling around a circumference of a pipe without a defect;

Figure 20 illustrates a transmitter-receiver arrangement around the circumference of a pipe and a support-type defect;

Figure 21 is a schematic plot of frequency against time of the SHo and AHi modes travelling around a circumference of a pipe with a defect, such that the AHi passes through the defect and reflection and transmission occurs, wherein the system measures the cut-off frequency fc and, thus, calculates minimum remaining wall thickness Lmin;

Figure 22 is a schematic plot of wavelength against frequency for SH0, AHi and SHi modes travelling circumferentially around a pipe having a wall thickness of 10 mm and illustrating a probe frequency range of 180 to 325 kHz and a wavelength range of 11 to 41 mm to detect defects extending 1 to 5 mm into the pipe wall; and Figure 23 is a schematic diagram of a three-element sensor showing a magnet arrangement in which a centre magnet is fixed and adjacent magnets can be adjusted so as to vary spacing.

Detailed Description of Certain Embodiments

Defect in pipe

The purpose of a guided wave sensing system herein described is to detect and estimate the depth of support-type defects and other similar types of defects and flaws found in pipes used for carrying, for example, oil, gas, petrochemicals, other types of processed chemicals, water or coolants. These types of pipes can be found in a wide-range of industries and applications such as oil and gas extraction platforms, refineries, chemical plants, pipelines, power stations and other types of industrial plant and infrastructure.

Referring to Figure 1, a pipe defect 1 can take the form of a cavity or pit formed in the wall 2 of a pipe 3 resulting in a substantial reduction of wall thickness, L. A defect 1 can extend lengthwise and/or circumferentially by a distance, d, ranging from a few millimetres to a few hundreds of millimetres on the inside 4 or outside 3 of the pipe 3. The defect 1 can have a depth Lmax which can be X % of the original wall thickness, where ο < X < too. For example, a defect 1 can be relatively shallow, e.g. ο % < X % < 50 %, or relatively deep, 50 % < X % <100 %. The cavity can be filled with material which is different to that of the pipe wall 2. For example, the pipe wall 2 may be formed from steel and the cavity may be filled with rust.

Ultrasonic guided wave shear horizontal sensing system

Referring to Figure 2, a sensing system 6 employing guided waves is shown which can be used for detecting and sizing defects.

The system 6 includes a transmitter 7 and a receiver 8 comprising respective transducers 9,10 capable of generating and detecting shear horizontal (SH) guided waves in a pipe wall 2. The transducers 9,10 take the form of electromagnetic acoustic transducers (EMATs). However, other forms of transduction such as piezoelectric transduction can be implemented as well.

The system 6 includes a signal generator 11 capable of generating rf signals 12 having a suitable frequency, which is usually of the order of tens or hundreds of kilohertz (kHz), and a suitable shape, such as, for example, a n-cycle suitably-windowed tone burst, where n is a positive number equal to or greater than 1, preferably an integer or half integer, preferably taking a value in the range 3 > n > 10, and where a suitable windowing function can be a Gaussian function. The signal generator 11 feeds the rf signal 12 to the transmitter transducer 9 which converts the signal 12 into a guided wave 13 in the pipe wall 2, for example, by Lorentz force and/or magnetostriction transduction. The guided wave 13 travels along the pipe wall 2 and may be transmitted and reflected at a discontinuity.

The receiver transducer 10 converts a (received) guided wave 13 into an electrical signal 15. The receiver transducer 10 feeds the electrical signal 15 to a signal receiver 16. The signal receiver 16 may include an amplifier (not shown) and an analogue-to-digital converter (not shown) which generates a digitized signal 17 of the electrical signal 15.

The system 6 includes a computer system 18 for controlling the signal generator 11 and processing signals 17 received from the signal receiver 16. The computer system 18 includes at least one processor 19 and memory 20 which stores software 21 for implementing a defect-detection and defect-sizing program.

The system 6 may also include a display 22, for example, to allow measurements to be displayed and one or more user input devices 23. The display 22 may take the form of a touch display which also provides user input.

The sensing system 6 is capable of automatically measuring pipe diameter and the pipe wall thickness, L. These two basic parameters are then used by the computer program in detecting defects and measuring their size.

Calculating the effective path length (the distance) between transmitter and receiver Referring also to Figure 3 and 14, the transducers 9,10 are placed against or suitably close to the steel pipe wall 2 in an annular section 30 of the pipe 3 (in other words, at the same longitudinal position, z, along the pipe) at different angular positions (in other words, at different angular positions, Θ, around the circumference of the pipe, as shown in Figure 14).

Referring also to Figure 4, the transmitter 7 generates a fundamental shear horizontal mode (SH0) guided wave tone burst propagating in the direction of receiver 8. The system 6 measures the time taken by the wave 13 to travel from the transmitter transducer 9 to the receiver transducer 10.

The effective path length (the distance), Dtr, between transmitter 7 and the receiver 8 can be calculated as:

(1) where c is the speed of the SH0 mode in steel and T0f is the time of flight between transmitter 7 and the receiver 8.

For example, the transmitter 7 generates a short tone burst SH0 signal at a time = 0.00 ps. The receiver 8 picks up the tone burst at a time = 100.0 ps. The speed, c, of the SH0 mode in steel is 3250.0 m/s. If the system 6 measures the time of flight T0/to be 100.0 ps, then, using equation (1) above, the system 6 can calculate the effective path length (the distance), Dtr, between the transmitter 7 and receiver 8 to be 325.0 mm.

Calculating pipe wall thickness

Referring to Figure 5, the transducers 9,10 are placed some distance apart and facing each other along the circumference of the steel pipe.

Referring also to Figure 6, the transmitter 7 generates a dispersive shear horizontal mode SHn (where n is a positive integer) tone burst wave propagating in the direction of the receiver 8. The time taken by the wave to travel from transmitter transducer 9 to receiver transducer 10 and the dominant frequency of the received signal are measured by the system 6 (Figure 2).

The group velocity for the measured frequency of the SHn mode can be obtained using:

(2) where Dtr is the effective path length (the distance), between the transmitter transducer 9 to the receiver transducer 10 and 71,/-is the time of flight taken for the wave to travel between the transmitter transducer 9 to the receiver transducer 10.

Using the dispersion equation, the thickness, L, of the pipe wall 2 can be calculated using

(3) where υ is the group velocity of the SHn tone burst,/is the dominant frequency of the SHn tone burst, c is the bulk shear velocity in steel, and n is mode number, where n = 1, 2, ...

For example, the AHi tone burst signal travels from the transmitter 7 to the receiver 8. The transducers 9,10 are placed at the distance Dtr = 185.5 mm apart.

The system 6 (Figure 2) selects the part of the AHi tone burst within a time window 31 centred at Tof = 103.6 ps. The system 6 takes the fast Fourier transform (FFT) of the windowed signal and finds the dominant frequency f = 254.5 kHz.

Using equation (2) above, the group velocity, υ, of the AHi mode (that is, for n=i) associated with this frequency can be obtained as υ = 185.5 / 103.6 = 1790.5 m/s. Bulk shear velocity, c, in steel is 3250 m/s. Using equation (3) above the system 6 can calculate the pipe wall thickness, L, as being 7.65 mm.

Calculating depth of the defect using reflected waves

Referring to Figure 7, the transducers 9,10 are positioned some distance apart and facing each other along the circumference of the pipe 3. The receiver transducer 10 is situated between the transmitter transducer 9 and the defect 1.

Referring also to Figure 8, the transmitter 7 generates an AHi mode tone burst. The travelling AHi wave is picked by the receiver 8. The wave 13 travels further and is reflected by the defect 1. The reflected wave 14 is subsequently picked up by the receiver 8. The time domain view of the signal is shown in Figure 8.

Referring to Figures 9 and 10, the Fourier transforms 35, 36 of the original AHi signal and the reflected signal are shown.

The reflected signal 14 has a different frequency composition from that of the original signal 13. Specifically, higher frequency content will be significantly attenuated above a specific cut-off frequency, fc.

The system 6 calculates a minimum wall thickness, Lmm, under the defect 1 using:

(4) where c is the bulk shear velocity in steel and fc is the cut-off frequency of the reflected signal.

For example, the shear wave velocity, c, in steel is 3250.0 m/s. If the cut-off frequency, fc, is measured as 280.0 kHz (see Figure 10), then, using equation (4) above, the system 6 determines the minimum wall thickness, Lmin, to be 5.8 mm. If the pipe wall thickness is 7.7 mm then the depth of the defect is 7.7 - 5.8 = 1.9 mm.

Calculating depth of the defect using transmitted waves

Referring to Figure 11, the transducers 9,10 are positioned some distance apart and facing each other along the circumference of the steel pipe 2. The defect 1 is situated between the transmitter transducer 9 and the receiver transducer 10.

Referring to Figure 12, the transmitter 7 generates a dispersive AHi mode tone burst 13. The travelling wave crosses the defective zone and the transmitted wave 13’ is picked by the receiver 8.

Referring also to Figure 13, the Fourier transform 38 of the transmitted signal 13 of Figure 12 is shown. The wave travelling through the defect area loses some of the lower frequency components 39. Below a certain cut-off frequency 40, the signal is attenuated.

The system 6 (Figure 2) calculates the minimum wall thickness, Lm!n, under the defect 1 as:

(5) where c is the bulk shear velocity in steel and fc is the cut-off frequency of the reflected signal.

For example, the bulk shear wave velocity in steel is assumed to be 3250.0 m/s. If the system 6 (Figure 2) finds that fc is 325.0 kHz (Figure 13), then, using equation (5) above, the system 6 calculates the minimum wall thickness Lm!n to be 5.0 mm. If the pipe wall thickness is 7.7 mm then the depth of the defect is 7.7 - 5.0 = 2.7 mm.

Defect sizing process

Referring to Figure 14, an arrangement of a transmitter 7 and a receiver 8 around a pipe 3 is shown.

Pipe circumference

The transmitter 7 is arranged to generate shear waves in the pipe wall 2 propagating around the pipe circumference. The transmitter 7 generates a SH0 mode tone burst wave which is picked up by the receiver 8.

Referring also to Figure 15, multiple, round-trip wavelets 131,132,133,134 travelling circumferentially in the clockwise and anticlockwise directions are detected by the system 6 and, based on time of flight, are used to calculate accurately the effective path length (the distance), Dtr, between transmitter transducer 9 and the receiver transducer 10 and also pipe circumference, Circ·

Both variables Dtr and Circ are measured by the system 6 (Figure 2) to a precision equal to or less than 0.1 mm.

Pipe wall thickness

Using the dispersion properties of higher-order modes, the system 6 (Figure 2) accurately measures the pipe wall thickness using a time-frequency decomposition of measured signals. For practical purposes, the system 6 (Figure 2) uses the non-dispersive SH0 mode and the first dispersive mode AHi, that is, modes number n = 0 and n=i of equation (3).

Referring also to Figure 16, the system 6 (Figure 2) can generate a predicted (or “theoretical”) time-frequency map 40 (herein also referred to as “setting-up a time-frequency map”) for the pipe given the measured values of transmitter-receiver distance, Dtr, pipe circumference Circ, and the pipe wall thickness, L, as will now be described in more detail.

Setting up a time-frequency map

Referring also to Figure 14, the transmitter 7 and receiver 8 are set up on the pipe 3 around the pipe circumference. Circumferential distances A, B between the transmitter 7 to receiver 8 can be measured along top and bottom paths respectively.

The pipe circumference C can be calculated as: C = A + B (6)

For example, referring also to Figure 15, using the time-of-flights for first and second wavelets 131,132 and the bulk shear wave velocity of 3250 ms-1, the circumferential distances for the top and bottom paths can be determined to be A = 57.1 x 3250 = 185.5 mm and B = 150.2 x 3250 = 488.0 mm respectively and so the circumference C = A + B = 673.5 mm.

Equation (3) above can be re-arranged to give:

(7)

Therefore, the time-of-arrival of a wavelet can be calculated as:

(8) where D is the path length between transmitter 7 and receiver 8. D can be the top path A, the bottom path B or the multiple circumferential paths, i.e.: D = A, B, A+C, B+C, A+2C, B+2C,.... (9)

For example, to generate time-frequency curve of the AHi mode (i.e. for mode n = 1) for the top path, the system 6 (Figure 2) uses equation (8) above and the values of n, D, A and c, as shown in Figure 16. For n = 1, D = A = 185.5 mm, and c = 3250 ms-1:

(10)

The system 6 generates AHi mode curves for other circumferential paths A, B, A+C, B+C, A+2C, B+2C etc. using equation (8) and obtains a family of time-frequency curves as shown in Figure 16.

The position of the AHi curves 4I1, 4i2 4¾ on this time-frequency map 40 is determined by the pipe wall thickness, L. For larger values of pipe wall thickness, L, the AHi curves are shifted down, while for smaller values of pipe wall thickness, L, the AHi curves are shifted up.

The system 6 (Figure 2) calculates pipe wall thickness, L, using a number of points using equation (3) along the AHi lines on the time-frequency map 40. Typically, the wall thickness is measured to within a small fraction of a millimetre.

In general, time-frequency curves can be calculated for any SHn mode in the pipe 3, using equation (8), for n= o, 1, 2,... etc.

The system 6 (Figure 2) generates families of time-frequency curves associated with SHo, AHi and higher modes and creates the reference time-frequency map (Figure 19) which is then used by the system 6 (Figure 2) to detect and evaluate defects 1.

Detecting support-type defects

Referring to Figures 2,14 and 17 to 19, to detect and determine the depth of a defect 1, a multiple step process is used.

First, the system 6 is used to inspect a defect-free section of the pipe 3.

The transmitter 7 and receiver 8 are positioned in a reference section of a pipe 3 along its circumference (step Si) and the computer system 18 causes the transmitter 7 to transmit SH0, AHi or higher modes’ tone bursts at one or more frequencies (step S2).

If necessary, higher-order modes can be chosen.

The computer system 18 receives measurements of SH0, AHi or higher modes from the receiver 8 (step S3). An illustrative time-domain response is shown in Figure 18 which includes received SH0 wavelets 431, 432,433 and AHi wavelets 441,442, 443.

The computer system 18 uses the measurements to calculate the transmitter-receiver distance, Dtr, the pipe circumference, Circ, and wall thickness, L (steps S4, S5 & S6).

In many practical applications, the system 6 firstly generates only SH0 mode using frequency tone burst below the AHi cut-off to accurately calculate the transmitter-receiver distance, Dtr and the pipe circumference, GrC· Subsequently the system 6 uses higher frequencies to generate higher modes in order to measure wall thickness, L.

Referring in particular to Figure 19, the computer system 18 generates a reference time -frequency map 45 (step S7). The time -frequency map 45 includes lines 461,462, 463 corresponding to SH0 modes and curves 471, 472, 473 corresponding to AHi modes. It will be appreciated that the time-frequency map 45 need not be actually displayed, but can simply be stored in memory.

The system 6 is able to probe (or “access”) a range 48 of frequencies between first and second frequencies fA, ίβ· Thus, the system is about to probe sections 49,50 of the lines 461, 462, 463 and curves 471, 472, 473 (i.e. the sections 49, 50 of the lines and curves in the hatched areas). The system 6 maybe able to sweep the frequency range 48 electronically. However, as will be explained in more detail later, the system 6 may use a different set of transducers with a different spacing of elements or transducers in which it is possible to controllably vary spacing of elements so as to access wider wavelength range and, thus, access a broader aggregate frequency range.

Having established a predicted time -frequency map 45, the system 6 is used then to inspect a section of the pipe 3 for defects.

Referring also to Figure 20, the transmitter 7 and receiver 8 are positioned in a section of the pipe 3 under test along its circumference (step S8). Although defects 1 can be located in lower parts of the pipe 3, e.g. due to pipe supports, they can be found anywhere around the pipe 3.

Starting with a first frequency, e.g. at the lower end of the range, the system 6 probes the pipe 3 across a range of frequencies (step S9 & S10)

The computer system 18 causes the transmitter 7 to generate the required modes, in this case, a SH0 wave and an AHi wave using a tone burst at a frequency^ (step S10). The computer system 18 receives and stores measurements of SH0 and AHi from the receiver 8 (steps S11 & S12).

The computer system 18 may continue to probe the pipe at different frequencies until measurements for the complete range of frequencies have been acquired (steps S13 & S14), in case when, it might be necessary to analyse higher order modes. Thus, the computer system 18 effectively constructs a measured time-frequency map 51.

Referring also to Figure 21, some AHi curves appear in different positions on the time-frequency map. As the dispersive mode wave passes over the defect 1, a frequency split may take place between the reflected and transmitted waves so that specific patterns can be observed on the time-frequency map.

These different patterns, in particular the cut-off frequency split between the reflected and transmitted signals, are used by the system 6 to calculate the maximum depth,

Dmax, of the defect 1.

The computer system 18 identifies the frequency split and calculates the cut-off frequency fc (see Figure 21). The minimum remaining wall thickness is calculated using:

(11) where c is the bulk shear velocity in steel and fc is the frequency at which the split between reflected and transmitted signal takes place. The computer system 18 can then calculate the maximum depth, Dmax, of the defect 1 using:

Dmax = L - Lmin (l2)

For example, the pipe wall thickness L is 7.7 mm and the split frequency, fc, is 325.0 kHz. Using equation (6) above, the minimum wall thickness Lm!n is 5.0 mm and, using equation (7) above, the system 6 calculates the maximum depth of the defect Dmax to be 7.7 - 5.0 = 2.7 mm.

The computer system 18 compares the measured time-frequency map 51 with the predicted time-frequency map 45 to identify missing and/or shifted portions 52 of AHi curves 53 (step S15) and, if found, identifies the split frequency fc (step S16).

Using the split frequency fc, the computer system 18 is able to calculate the minimum remaining wall thickness (step S17) and, this, the maximum depth (step S18).

Transducer

To detect and successfully size defects, the transducers 9,10 cover a sufficiently wide range of frequencies and corresponding wavelengths. Such a requirement can pose a challenge for multi-element type transducers, such as periodic permanent magnets EMAT transducers, which might only support a narrow range of the wavelengths.

Referring to Figure 22, a dispersion map is shown for shear horizontal Lamb waves propagating circumferentially in a 203.2-mm (8-inch) diameter pipe having a wall thickness of 10 mm. The wave dispersion is computed as a relationship of wave frequency versus its wavelength for the SH0, ALL and SLL modes.

Using equation (7) above, for a pipe having a wall thickness, L, of 10 mm, if the system 6 is to be able to measure a maximum defect depth, Dmax, in a range between 1 and 5 mm, then the transducers 9,10 should be able to generate and detect ALL waves in a frequency range of 180 to 325 kHz.

This frequency range translates into a range of AHi mode wavelengths of between 11 and 41 mm. Similar analysis can be done for other mode types such as SHi, AH2 etc.

The required frequency range can be readily obtained using modern electronic circuitry. However, wide wavelength range can be difficult to achieve when using multielement probes.

To be able to sweep across a wide range of the wavelengths, the transducers 9,10 can incorporate a mechanism allowing smooth adjustment of the magnet spacing.

Figure 23 shows the magnet arrangement for the three-element sensor comprising first, second and third magnets 61, 62, 63. The second magnet 62 is interposed between the first and third magnets 61, 63 and is fixed. The first and third magnets 61, 63 are moveable so as to vary spacing, s. The spacing, s, be variable, for example, between 1 mm and 20 mm.

This or similar arrangements can be readily incorporated into an EMAT or other form of multi-element ultrasonic transducer to achieve variable spacing.

It will be appreciated that various modifications may be made to the embodiments hereinbefore described. Such modifications may involve equivalent and other features which are already known in the design, manufacture and use of ultrasonic guided wave systems and component parts thereof and which may be used instead of or in addition to features already described herein. Features of one embodiment may be replaced or supplemented by features of another embodiment. A single transducer, particularly a piezoelectric transducer, can be used which serves as both transmitter and receiver. EMAT transducers need not be used. Instead, piezoelectric or magnetostrictive or other transducers can be used.

The wall- or plate-like structure can be a wall or plate, a pipe or a vessel, such as a cylindrical or rectangular tank with closed ends, or hollow bar or hollow box, such as a channel or conduit. The structure can be formed from a metal, metal alloy, plastic or any solid material. The type of material may affect the type of transducer used.

The defect need not be a corrosion-based defect.

Although claims have been formulated in this application to particular combinations of features, it should be understood that the scope of the disclosure of the present invention also includes any novel features or any novel combination of features disclosed herein either explicitly or implicitly or any generalization thereof, whether or not it relates to the same invention as presently claimed in any claim and whether or not it mitigates any or all of the same technical problems as does the present invention. The applicants hereby give notice that new claims may be formulated to such features and/or combinations of such features during the prosecution of the present application or of any further application derived therefrom.

Claims

Claims

1. A method of determining a thickness of a region of a wall- or plate-like structure which is thinner than a thickness of a surrounding region of the structure due to a cavity in the structure, the method comprising: comparing a measured time-frequency dispersion map for at least one dispersive guided wave obtained by measuring the structure using guided waves with a reference time-frequency dispersion map obtained by modelling the structure; determining a split frequency at which the measured time-frequency dispersion map and the reference time-frequency dispersion map differ; and calculating the thickness of the thinner wall thickness region in dependence upon the split frequency.
2. A method according to claim l, wherein the at least one dispersive guided wave is at least one antisymmetric shear horizontal guided wave or at least one symmetric shear horizontal guided wave
3. A method according to claim 1 or 2, further comprising: obtaining the measured time-frequency dispersion map by, for a set of different frequencies: causing a transmitter transducer to excite the structure using a tone burst at a given frequency; receiving a time-varying signal from a receiver transducer; and storing the time-varying signal.
4. A method according to claim 3, wherein causing the transmitter transducer to excite the structure using the tone burst at the given frequency causes a fundamental symmetric shear horizontal, SH0, guided wave and a first-excited mode antisymmetric shear horizontal, AH,, guided wave to be generated in the structure.
5. A method according to any preceding claim, wherein modelling the structure comprises generating the reference time-frequency dispersion map in dependence upon a thickness of the structure.
6. A method according to claim 5, wherein generating the reference time-frequency dispersion map comprises calculating a time of flight as a function of frequency or vice versa dependent on guided wave mode number, a distance between a transmitter transducer and a receiver transducer, a thickness of a wall or plate of the structure and bulk shear velocity.
7. A method according to any preceding claim, wherein modelling the structure comprises: for a reference section of the structure; causing the transmitter transducer to excite the structure using reference tone bursts; receiving reference time-varying signals from a receiver transducer; and determining a distance between a transmitter transducer and a receiver transducer and a thickness of a wall or plate of the structure in dependence upon the reference time-varying signal.
8. A method according to any preceding claim, wherein the wall- or plate-like structure is a pipe comprising a pipe wall.
9. A computer program which, when executed by a computer, causes the computer system to perform a method according to any one of claims 1 to 8.
10. A computer-readable medium, optionally a non-transitory computer-readable medium, storing a computer program according to claim 9.
11. Apparatus configured to compare a measured time-frequency dispersion map for at least one dispersive guided wave obtained by measuring the structure using guided waves with a reference time-frequency dispersion map obtained by modelling the structure, to determine a cut-off frequency at which the measured time-frequency dispersion map and the reference time-frequency dispersion map differ and to calculating the thickness of the thinner region in dependence upon the cut-off frequency.
12. Apparatus according to claim 11, wherein, for a set of different frequencies, to obtain the measured time-frequency dispersion map, the apparatus is configured to cause a transmitter transducer to excite the structure using a tone burst at a given frequency, to receive a time-varying signal from a receiver transducer and to store the time-varying signal. 13· Apparatus according to claim 12, further comprising: causing a transducer to be mechanically adjusted so as to access a different range of wavelengths.
14. Apparatus according to claim 13, wherein causing the transducer to be mechanically adjusted so as to access the different range of wavelengths comprises causing spacing of elements in the transducer to change.
15. Apparatus according to any one of claims 11 to 14, which is a computer system comprising memory and at least one processor.
16. A guided wave inspection system comprising: apparatus according to any one of claims 11 to 15; a transmitter comprising a transmitter transducer for exciting guided wave modes in the structure; and a receiver comprising a receiver transducer for generating a time-varying signal.
17. A system according to claim 15, wherein the transmitter and receiver transducers comprise respective transducers.
18. A system according to claim 16 or 17, wherein the transmitter transducer and receive transducer comprise an electromagnetic acoustic transducer or first and second electromagnetic acoustic transducers.
19. A transducer for guided wave inspection, comprising: an array of spaced-apart transducer elements which are arranged such that spacing between the transducer elements is variably adjustable.
20. A transducer according to claim 19, wherein the transducer is an electromagnetic acoustic transducer and the transducer elements comprise permanent magnets.
21. A transducer according to claim 19 or 20, which is configured such that the spacing is variable so as to take a value between a lower limit and an upper limit, wherein the lower limit lies in a range between 1 mm and 10 mm and the upper limit lies in a range between 10 mm and 20 mm.