System and method for rapid acoustic assessment of pipe obstructions
11668819 · 2023-06-06
Assignee
Inventors
Cpc classification
G01S15/104
PHYSICS
International classification
Abstract
A method for locating and quantifying obstructions in a pipe is described. The method comprises the steps of emitting at one end of the pipe, by means of a loudspeaker, an emitted signal comprising a wave train at a first frequency, the waves being of an acoustic type; receiving at the same end of the pipe, by means of a microphone, a reflected signal, resulting from the reflection of the wave train from obstructions in the pipe; determining a position of each of the obstructions according to a delay measured between the wave train of the emitted signal and the wave train of the reflected signal received by the microphone; and determining, for each of the obstructions, its degree of obstruction by extrapolating the energies of the emitted, reflected and transmitted signals. A system comprising the loudspeaker, the microphone, a processor, and a memory for performing the above method is also described.
Claims
1. A method for locating and quantifying obstructions in a pipe, the method comprising: emitting, at an end of the pipe, through a loudspeaker, an emitted signal comprising a wave train, the waves being of an acoustic type; receiving at the same end of the pipe, through a microphone, a reflected signal resulting from a reflection of the wave train on the obstructions in the pipe; determining a position of each of the obstructions based on a time delay measured between the wave train of the emitted signal and the wave train of the reflected signal received by the microphone; determining, for a most upstream obstruction, its degree of obstruction by performing the steps of: a) extrapolating an energy of the emitted signal and an energy of the reflected signal at the position of the obstruction; b) calculating a reflection coefficient C.sub.D=E.sub.r/E.sub.i, where E.sub.i is the extrapolated emitted signal energy and E.sub.r is the extrapolated reflected signal energy; c) calculating a ratio a/A of a partly obstructed section and a complete section of the pipe, the ratio a/A being determined by a function of the reflection coefficient C.sub.D; and determining, for subsequent obstructions, their respective degrees of obstruction by repeating steps a) through c) for each one, and by replacing the extrapolated emitted signal energy E.sub.i with a transmitted energy E.sub.t once the upstream obstruction has been passed.
2. The method according to claim 1, wherein the step of emitting the emitted signal comprises a step of forming the wave train with a number n of cycles, where n≥3.
3. The method according to claim 1, comprising: determining a number of close proximity obstructions m by estimating, from the reflected signal, a number of overlapped cycles n.sub.t greater than the number of emitted cycles n in the wave train, such that m is a rounded-up integer of a ratio (n.sub.t/n).
4. The method according to claim 1, wherein the wave train is a first wave train emitted at a first frequency, the first frequency corresponding approximately to a first wavelength of 3 times a diameter of the pipe.
5. The method according to claim 4, wherein the emitted signal comprises a second wave train at a second frequency corresponding approximately to a second wavelength of 10 times the diameter of the pipe, the first and second wave trains permitting to detect different types of obstructions; wherein the method further comprises a step of corroborating the respective positions of the obstructions and their degrees of obstruction using the results obtained with the first and second wave trains.
6. The method according claim 1, wherein the reflected signal comprises a first signal reflected from the obstructions of the pipe, downstream of the microphone, and a second signal reflected from the loudspeaker, upstream of the microphone, the method comprising: eliminating the second signal reflected by emitting through the loudspeaker, a wave train equal to but in opposite phase with the second signal reflected from the loudspeaker.
7. The method according to claim 1, wherein the function that determines the ratio a/A with respect to the reflection coefficient C.sub.D is given by a/A=P.sub.1×C.sub.D+P.sub.2, P.sub.1 being a first constant between 0.1 and 2 and P.sub.2 being a second constant between −1 and 1, according to the frequency used.
8. The method according to claim 7, wherein P.sub.1 is equal to 1 and P.sub.2 is equal to 0.
9. The method according to claim 1, wherein the transmitted energy E.sub.t corresponds to the extrapolated emitted signal energy E.sub.i multiplied by [1−C.sub.D.sup.C.sup.
10. The method according to claim 9, wherein the value of C.sub.trans is equal to 2.55 when the detected element is a distinct obstruction.
11. The method according to claim 9, wherein the value of C.sub.trans is equal to 1 when the detected element is a close proximity obstruction.
12. The method according to claim 9, wherein the value of C.sub.trans is equal to 6 when the detected element is a cavity.
13. The method according to claim 1, wherein the step of extrapolating comprises a step of applying an attenuation factor on the emitted and reflected signals, based on a distance of a given obstruction from the end of the pipe.
14. The method according to claim 13, comprising a step of simulating a complete blockage of the pipe at a predetermined distance from the end of the pipe, to determine the attenuation factor to be applied in the step of extrapolating the emitted signal and the reflected signal to the position of the obstructions.
15. The method according to claim 1, comprising: detecting a cavity in the pipe by identifying, in the reflected signal, a wave train in anti-phase relative to the wave train of the emitted signal.
16. The method according to claim 5, comprising: identifying a bend in the pipe when the microphone receives the reflected signal of the first wave train but does not receive the reflected signal of the second wave train.
17. A portable system for locating and quantifying obstructions in a pipe, the system comprising: a loudspeaker, configured to emit at an end of the pipe an emitted signal comprising a wave train, the waves being of an acoustic type; a microphone, configured to receive, at the same end of the pipe, a reflected signal, resulting from a reflection of the wave train on the obstructions in the pipe; a processor and a memory in communication with the processor, the processor and the memory being operatively connected to the microphone and loudspeaker, the memory comprising instructions executable by the processor to: determine a position of each of the obstructions based on a measured time delay between the emitted signal wave train and the reflected signal wave train received by the microphone; and determine, for a most upstream obstruction, its degree of obstruction by performing the steps of: a) extrapolating an emitted signal energy and a reflected signal energy at the position of the obstruction; b) calculating a reflection coefficient C.sub.D=E.sub.r/E.sub.i, where E.sub.i is the extrapolated emitted signal energy and E.sub.r is the extrapolated reflected signal energy; and c) calculating a ratio a/A between a partly obstructed section and a complete section of the pipe, the ratio a/A being determined by a function of the reflection coefficient Co; wherein the processor determines, for subsequent obstructions, their respective degree of obstruction by repeating steps a) to c) for each one and by replacing the extrapolated emitted signal energy E.sub.i with a transmitted energy E.sub.t once the upstream obstruction has been passed.
18. The system according to claim 17, wherein the loudspeaker and the microphone are disposed in a housing, and wherein the processor and the memory are offset from the housing.
19. The system, according to claim 18, wherein the housing is positioned at the end of the pipe, with the loudspeaker and microphone aligned on a central axis of the pipe.
20. The system, according to claim 17, wherein the loudspeaker is located approximately 3 cm upstream of the microphone relative to the end of the pipe.
Description
BRIEF DESCRIPTION OF THE DRAWINGS
(1) A detailed description of the invention will be given below with reference to the following drawings:
(2)
(3)
(4)
(5)
(6)
(7)
(8)
(9)
(10)
(11)
DETAILED DESCRIPTION
(12) In the description and figures that follow, the same reference numbers refer to similar elements of the invention. Further, in order not to unduly clutter the figures, it is possible that a figure may not contain all of the reference numbers of the elements therein. Thus, it is also possible that certain elements or components are referenced in only one figure. The elements thus referenced can be easily inferred in the other figures presented. The embodiments, geometric configurations, materials and/or dimensions shown in the figures or described in the present disclosure are only indicative, and show possible embodiments, presented as examples, and should not be construed as limitations of the invention.
(13) Referring to
(14) The system 10 includes means for emitting signals, such as a loudspeaker 100 and means for receiving signals, such as a microphone 110. The loudspeaker 100 is configured to emit an emitted signal comprising a wave train, or a train of waves, the waves being of an acoustic type. The loudspeaker used may be a miniature loudspeaker, resistant to outdoor use thanks to moisture protection properties. It can also be used in a range of outdoor temperatures, for example between −40° C. and +80° C. According to one possible embodiment, the loudspeaker is dimensioned so that its cone or the circumference of its membrane is equal to or smaller than the diameter of the pipe to be inspected.
(15) The microphone 110 is configured to receive a reflected signal, resulting from the reflection of the wave train emitted by the loudspeaker on the obstructions in the pipe. The microphone used may be of the low noise, pre-polarized measurement microphone type. In some cases, the microphone may contain a built-in preamplifier. The microphone is also preferably sized so that it can be inserted, in whole or in part, at the end of the pipe where the loudspeaker is located.
(16)
(17) The emitted signal as well as the reflected signal are digitized by the acquisition software. The signals are recorded at a sampling rate, for example 51.2 kHz. The data can be recorded and stored in the computer memory 130 for analysis. According to one possible embodiment, the analyzed signals may correspond to the result of a given number of time averages, for example 5 time averages, synchronized from the signal transmitted via the audio card. The purpose of the time averages is to denoise the signal measured by the microphone.
(18) In one embodiment, the loudspeaker 100 and microphone 110 may be inserted into a housing 112. This housing 112 may be made of PVC, rubber, plastic, or any relatively rigid material. The housing allows the loudspeaker 100 and microphone 110 to be held in an aligned position at a given distance from each other. In one possible embodiment, the distance between the loudspeaker and the microphone is between 1 and 10 cm, and preferably about 3 cm. In one embodiment, the housing may comprise a PVC tube 122 about 5 cm in diameter and 30 cm in length, and a rubber adapter 124 attached to the tube 122, the adapter 124 having an enlarged diameter relative to the tube 122, so as to exactly match the diameter of the pipe 20. The tube 122 allows for better handling and grip of the portable system 10. In use, the portable system 10 is positioned at a proximal end 200 of the pipe 20, such that the loudspeaker 100 and microphone 110 are aligned with the centerline 210 of the pipe. The configuration of the housing is for illustrative purposes only, and other configurations are possible. The housing may allow for increased portability of the system, so that the system can be easily and quickly installed and uninstalled at the proximal end of the pipe to be inspected.
(19)
(20)
(21) The method involves emitting an acoustic signal comprising at least a number of cycles n at a frequency f. The frequency f is determined by the ratio of a propagation speed V of the wave in the homogeneous medium determined by a wavelength λ of the signal.
(22)
In the illustrated embodiment, the homogeneous medium in which the signal propagates is air, for which the propagation speed V of the wave is equal to 340 m/s. The wavelength λ is determined with respect to the diameter d of the pipe. In order to ensure that the wave propagating inside the pipe is relatively flat, i.e., a wave whose wave fronts are relatively perpendicular to a direction of wave propagation, its wavelength λ must be at least 2 times the diameter d of the pipe.
(23) In order to apply the developed models for predicting the degree of obstruction, two frequencies f1 and f2 were selected. Considering for example a pipe of 10 cm diameter (d=10 cm), a first wavelength λ1 of about 2 to 4 times the diameter, and preferably 3 times the diameter of the pipe (λ1=3d=30 cm) and a second wavelength λ2 of about 8 to 12 times the diameter, and preferably 10 times the diameter of the pipe (λ2=10d=100 cm=1 m) can be selected. In this case, the corresponding frequencies f1 and f2 are calculated such that:
(24)
(25) The value of the first wavelength λ1, corresponding to about 3 times the pipe diameter, allows more spatial resolution of obstructions, especially in the case of smooth pipes such as PVC pipes. By spatial resolution, we mean that the position and distinction of individual obstructions is better identified than with other wavelengths. In addition, the value of λ1 makes it possible to identify a succession of close proximity obstructions up to a critical distance d.sub.cr of
(26)
(27) The value of the second wavelength λ2, corresponding to about 10 times the diameter of the pipe, performs better in the case of pipes with rough, cracked or irregular walls and in the case where the energy dissipated through the wall is greater, such as for pipes with concrete walls. The value of the second wavelength λ2 can also be used to detect certain deformations in the pipe, such as a bend or a cavity.
(28) Thus, according to an embodiment, the signal emitted by the loudspeaker 100 comprises at least a first wave train at a first frequency f1. As explained above, this frequency f1 may be between 850 and 1700 Hz, and preferably correspond to about 1100 Hz, when the pipe to be inspected has a diameter of about 10 cm and comprises air. In one possible embodiment, the signal emitted by the loudspeaker 100 comprises a second wave train at a second frequency f2. This frequency f2 may be between 283 and 425 Hz, and preferably corresponds to about 340 Hz, for the pipe described above. The emitted signal comprising the two wave trains at frequencies f1 and f2 will be reflected by the various obstructions clogging the pipe, resulting in first and second reflected signals, recorded by the microphone 110. These first and second reflected signals can then be corroborated, for example by comparing them, to confirm the position of the obstructions and their degree of obstruction. By degree of obstruction, we mean the ratio of obstructed sections vs. complete section of the drain. For example, it can be, according to a cross section of the drain at a distance x from the proximal end, the ratio of the area occupied by a limestone clump to the complete area of the drain. It has been shown that the use of two wave trains at two different frequencies allows better detection of all the irregularities in the pipe, compared to when a single frequency wave train is used. This is because either of the two frequencies may be more or less sensitive to the composition of the pipe walls, the spatial distribution of obstructions, or the presence of bends in the pipe. In other embodiments, it is possible to use a signal comprising more than two wave trains, at frequencies different from the two frequencies selected above. For the detection of calcareous obstructions in dam drains, transmitting two wave trains at wavelengths of about 3 and 10 times the diameter of the drain, has been found to be advantageous in locating and determining the degree of calcareous obstructions more accurately than when a single wavelength is used, without significantly increasing system complexity or analysis time.
(29) The number of cycles n of each emitted wave train may be determined as a compromise between i) a sufficient number of peaks to account for attenuation of the first and last peaks by the response time of the loudspeaker, and ii) a limited number of cycles to optimize detection of close proximity obstructions whose minimum critical distance d.sub.cr to identify two consecutive obstructions is given by
(30)
In the described embodiment, the optimized value of the number of cycles n is equal to or greater than 3
(31) Note that the values selected for d, λ1, λ2 and n are examples only, and other values may be selected, depending on the applications for which the system is to be used.
(32) In the embodiment shown in
(33)
(34) In order to determine a degree of obstruction of the first obstruction 220 encountered by the emitted signal (the obstruction furthest upstream or closest to the proximal end of the pipe 200), the energy of the emitted and reflected signals must be extrapolated. The energy E.sub.y of the acoustic wave train of a signal is determined by the product of an effective amplitude of the acoustic wave train over n cycles in A.sub.RMS Pascals and the pulse duration T in seconds:
E.sub.y=A.sub.RMS*T
(35) The subscript y denotes the signal that is being measured. With reference to
(36) E.sub.0 is the measured energy of the signal emitted at the proximal end of the pipe 200,
(37) E.sub.r0 is the measured energy of the signal reflected from the obstruction at the proximal end of the pipe 200.
(38) The value of this energy can also be extrapolated to the threshold of the obstruction 220:
(39) E.sub.i is the extrapolated energy of the emitted signal at the threshold of the obstruction 220,
(40) E.sub.r is the extrapolated energy of the signal reflected by the obstruction at the threshold of the obstruction 220,
(41) E.sub.t is the extrapolated energy of the signal transmitted through the obstruction 220, in the case of a partial obstruction.
(42) The extrapolated energy at the threshold of the obstruction can be evaluated by simulating the attenuation α of the wave as the acoustic pulse propagates along the pipe. The attenuation α is generated by friction, damping and leakage through the walls of the pipe 20. The calculation is based on the principle that the dissipated energy e.sub.F is proportional to the magnitude of its energy E by the attenuation coefficient α.
(43) Thus, the signal energy at a distance x from the proximal end of the pipe E.sub.(x) is given by:
E(x)=E.sub.0*e.sup.−αx
(44) In the embodiment shown in
(45) In one embodiment, calibration of the attenuation coefficient α yielded the following values:
(46) TABLE-US-00001 Diameter Frequency Attenuation Material (cm) (Hz) “α” PVC 7.6 1100 0.026 PVC 7.6 340 0.016 Fibro-ciment 10 1100 0.033 Fibro-ciment 10 340 0.037
(47) Returning to
(48) The reflection coefficient C.sub.D can be calculated by taking the ratio of the extrapolated energy of the signal reflected by the obstruction at the obstruction threshold E.sub.r to the extrapolated energy of the transmitted signal at the obstruction threshold Ei:
(49)
(50) The transmission coefficient C.sub.T can be calculated as the ratio of the extrapolated energy of the signal transmitted through the obstruction E.sub.t to the extrapolated energy of the signal emitted at the obstruction threshold Ei:
(51)
(52) This transmission coefficient C.sub.T is used to extrapolate the incident and reflected waves to the next obstruction. The following formula can be used to determine C.sub.T as a function of C.sub.D:
C.sub.T=(1−C.sub.D.sup.C.sup.
Where C.sub.trans is a coefficient of transmissibility obtained by a convergence algorithm based on data from laboratory tests (calibration). The value of C.sub.trans is between 1 and 6, depending on the type of detected element. By detected element, we mean an element among close proximity obstructions, distinct obstructions, a bend in the pipe, or a cavity or a widening and/or narrowing of the pipe section.
(53) This transmissibility coefficient C.sub.trans is determined by a calibration carried out during various laboratory tests, with obstructions positioned at known distances, on a calibration pipe having similar properties (wall roughness, diameter, material) to those of the pipe to be inspected. The different values of C.sub.trans will be detailed in the following description.
(54) Thus, when an obstruction has been identified, i.e. its position has been determined by the measured delay between the emitted signal wave train and the reflected signal wave train, its degree of obstruction can be determined. This degree of obstruction can be identified by a ratio a/A, where a represents the obstructed area and A represents the complete cross-sectional area of the pipe. This degree of obstruction is therefore a number between 0 and 1, where 0 represents no obstruction in the pipe and 1 represents total obstruction of the pipe.
(55) The degree of obstruction a/A can be determined from a function of the reflection coefficient C.sub.D. This function can be a polynomial function or any other function that relates the value of a/A to that of C.sub.D. For the most upstream obstruction in the pipe, i.e., the obstruction closest to the proximal end, the ratio a/A can be determined by a function of the reflection coefficient C.sub.D. For subsequent obstructions, their respective degree of obstruction is determined by repeating the steps described below, i.e., a) extrapolating emitted signal energy and a reflected signal energy at the position of the obstruction; b) calculating the degree of reflection C.sub.D, but replacing the extrapolated emitted signal energy E.sub.i with the transmitted energy E.sub.t once the upstream obstruction is passed, as described above; and finally c) calculating the ratio a/A is determined by a function of the reflection coefficient C.sub.D.
(56) According to one possible embodiment, where the pipe is an air-filled barrier drain, this function is a 1st degree polynomial function, or linear function, expressed by:
a/A=P.sub.1×C.sub.D+P.sub.2, where
P.sub.1 is a first constant between 0.1 and 2; and
P.sub.2 is a second constant between −1 and 1.
(57) The values of P.sub.1 and P.sub.2 depend mainly on the frequency of the emitted signal. The various tests carried out show that these constants are not very sensitive to the diameter of the pipe or the surface conditions, since these parameters are already taken into account in the determination of the reflection coefficient C.sub.D.
(58) Thus, according to this embodiment, P.sub.1=1 and P.sub.2=0, so that the degree of obstruction a/A corresponds to the reflection coefficient C.sub.D.
(59) In some cases, as shown in
(60) One way of suppressing this spurious reflected signal 114 is to emit a wave train equal to but in phase opposition to the spurious reflected signal through the loudspeaker, the so-called cancellation wave train 115. Among other things, a delay may be applied before emitting this wave train, in order to reflect the distance between the microphone and the loudspeaker. Thus, when the reflection bounces off the loudspeaker, the loudspeaker reacts with a mirror wave of this reflection in opposite phase so that the spurious reflected signal is cancelled, as shown in
(61) As shown in
(62) The phenomenon is caused by the approach of obstructions whose critical distance d.sub.cr between two consecutive obstructions is too short compared to the number of cycles n emitted, times the wavelength λ of the emitted signal.
(63) Knowing that
(64)
the criterion of chained reflections occurs when a distance D.sub.x that separates two consecutive obstructions is smaller than d.sub.cr. We talk about close proximity obstructions.
(65) For example, for an emitted signal with wavelength λ=340 Hz and a number of cycles n=3, chained reflections will occur if consecutive obstructions are less than 1.5 m apart.
(66) The calculation of the number of close proximity obstructions m in an area of chained reflections can be established from a total number of cycles of the reflected signal n.sub.t and the number of cycles of the emitted signal n, such that m is the rounded-up integer of the ratio (n.sub.t/n), with n.sub.t>n.
(67) In one embodiment and as illustrated in
(68) On the contrary, if the distance D.sub.x between two consecutive obstructions is greater than d.sub.cr, the obstructions are said to be distinct.
(69)
(70) Thus, in order to be able to determine the degree of obstruction a/A of several distinct consecutive obstructions, it is preferable to determine the transmissibility coefficient C.sub.trans for each obstruction. To do this, it is possible to start sequentially from the first obstruction. Then, the reflection coefficient C.sub.D and transmission coefficient C.sub.T can be calculated for each previous obstruction so that the incident and reflected energies can be extrapolated to the next obstruction. This procedure can be applied to all m.sup.th obstructions identified, as follows:
(71)
(72) For example, in the case illustrated in
C.sub.D1=E.sub.i1/E.sub.r1 and C.sub.D2=E.sub.i2/E.sub.r2.
(73) In some cases, as shown in
(74) When the emitted signal encounters a cavity 240, the reflected signal received by the loudspeaker 100 exhibits a 180-degree phase change. This antiphase 119 is illustrated in
(75) In one embodiment and as illustrated in
(76) In other embodiments, as shown in
(77) The three graphs shown in
(78) It should be noted that graphs 10-B and 10-C show a disparity 231. Indeed, the acoustic wave of the signal emitted at the frequency f.sub.1=1100 Hz (graph 10-B) detects the bend 230 while the acoustic wave of the signal emitted at the frequency f.sub.2=340 Hz (graph 10-C) does not detect the bend.
(79) Thus, a bend 230 in a pipe may be identified by corroborating the reflected signal from the first wave train emitted at a frequency f.sub.1 with the reflected signal from the second wave train emitted at a frequency f.sub.2. For example, this corroboration could involve superimposing the reflected signals from the first and second wave trains emitted at respective frequencies of f.sub.1 and f.sub.2 in order to identify the mismatch 231, since the reflected signal from the first wave train shows a series of pulses, while the reflected signal from the second wave train shows no significant pulses.
(80) In one embodiment, the computer processor 130 may automate all of the steps of the above method. An algorithm developed as part of signal processing and obstruction prediction software can identify and characterize the reflection, identifying in the signal measured by the microphone: A beginning and an end of a pulse of the emitted and reflected signal, The position of the first obstruction x, An effective amplitude A.sub.RMS, and A duration of the pulse T.
(81) Once these characteristics are identified, the software can identify each reflected signal according to 4 categories for a signal emitted at a frequency f.sub.2: 1. Distinct obstructions if the number of cycles of the reflected signal is less than or equal to the number of cycles of the emitted signal, 2. Close proximity and chained obstructions if the number of cycles of the reflected signal is higher or equal to the number of cycles of the emitted signal+1 cycle, 3. Cavities, or 4. Pipe outlet.
(82) In a subsequent step, the software may also allow the degree of obstruction and the coefficient of transmissibility of each obstruction to be determined.
(83) Although certain advantages have been described, the person skilled in the art may discover other advantages and/or features inherent in the invention that have not been explicitly described. Furthermore, although certain configurations and embodiments have been described herein, it is appreciated that they are by way of example only and should not be taken to limit the scope of the invention.