Influence of shear effects on the characteristics of axisymmetric wave propagation in a buried fluid-filled pipe

: The acoustic propagation characteristics of axisymmetric waves have been widely used in leak detection of fluid-filled pipes. The related acoustic methods and equipment are gradually coming to the market, but their theoretical research obviously lags behind the field practice, which seriously restricts the breakthrough and innovation of this technology. Based on the fully three-dimensional effect of the surrounding medium, a coupled motion equation of axisymmetric wave of buried liquid-filled pipes is derived in detail, a contact coefficient is used to express the coupling strength between surrounding medium and pipe, then, a general equation of motion was derived which contain the pipe soil lubrication contact, pipe soil compact contact and pipe in water and air. Finally, the corresponding numerical calculation model is established and solved used numerical method. The shear effects of the surrounding medium and the shear effects at the interface between surrounding medium and pipe are discussed in detail. The output indicates that the surrounding medium is to add mass to the pipe wall, but the shear effect is to add stiffness. With the consideration of the contact strength between the pipe and the medium, the additional mass and the pipe wall will resonate at a specific frequency, resulting in a significant increase in the radiation wave to the surrounding medium. The research contents have great guiding effect on the theory of acoustic wave propagation and the engineering application of leak detection technology in the buried pipe.


Introduction
With the development of urbanisation in past decades, the component like pipe has become a necessary means of liquid and gas transportation. The issue of plumbing leakage is a widespread concern in both industry and academy due to its social, environmental and economic utility. It has been estimated that, in mainland China, current direct economic losses caused by underground pipe network leaks exceed 200 billion Yuan (approximately 22.65 billion pound sterling), sometimes the pipe leakage might cause unexpected major source of hidden danger not only for damaging urban biological environment and the security of people's lives and properties. The latest research of the 'track and trace' technology for buried pipelines, which is widely used in the transportation of liquid and gas media such as crude oil and natural gas, is being widely discussed and studied extensively. Acoustic detection methods has attracted more attention because of such non-destructive evaluation does not directly destroy the structure of the original piping system [1][2][3][4].
When the leaking of the pipeline happened, the high-pressure fluid in the pipe will be pressured out of the pipeline and cause unavoidable noise, The acoustic leak detection method is used to detect the events of leaks in different locations of the pipeline, using cross-correlation technique to estimate the delay of leakage noise between two measuring points, then the location of leakage would be calculated [5]. The effectiveness of these methods much depends on the rationality of the selection of the propagation characteristic parameters of the dominant wave in the pipeline.
The studies of the acoustic characteristics and propagation mechanism of pipes is developed with pipe leak detection.
At present, the acoustic leak detection and location of pipelines are usually carried out by time delay estimation method, which depends strongly on wave propagation characteristics. The early research was mainly to solve the wave equation in thin-walled shell pipe, in which the external medium outside pipes is considered vacuum. Fuller, et.al [6] derived the propagation characteristics of n=0 wave within elastic fluid-filled pipes in vacuo defined as "hard" and "soft" shells, then the energy distribution of radial input force and internal pressure fluctuation under various waveforms is studied theoretically [7]. Xu, et.al [8] studied the vibration propagation characteristics of liquid filling pipes in a vacuum environment. Pinnington, et al [9] established an acoustic propagation model of cast iron pipes without considering the dispersion characteristics, after that, studied the n=0 wave propagation characteristic and transfer equation of the pressurised pipes [10,11]. With the development of pipe research, the influence of the medium around the pipes has gradually attracted people's attention. Sinha, et al. [12] studied the numerical results of acoustic wave propagation characteristics of fluid-filled pipes in infinite fluid. Greenspon [13] presented the axisymmetric vibration of thick-wall and thin-walled liquid-filled pipes in water medium. Long et al. [14] put forward a model of acoustic velocity dispersion in the process of acoustic vibration signal propagation and verified by experiments. Zhang, et al. [15] proposed a calculation model for sound velocity under different pipeline embedding conditions. Muggleton et al. [16][17][18] analysed the propagation characteristics of fluid-dominated axisymmetric waves (s=1, n=0) in filled buried pipes. Gao et al. [19,20] developed a general expression for the fluid-dominated wavenumber in a thin-walled fluid-filled pipe surrounded by a layered elastic soil, and the influence of load effect on elastic medium around pipeline is considered.
Current research studies reveal that at low frequencies, the fluid-dominated axisymmetric wave is not only the main carrying waveform of the vibration energy within the buried fluid-filled pipe, but also is an effective signal component which can be used for pipe leakage inspection. This waveform corresponds to the breathing mode of the pipeline, and the current researches on the problem are mainly focused on the metal pipeline. Due to the flexibility of the plastic pipe, the coupling between the pipe and the surrounding medium (mainly soil) is significant, making the influence of the acoustic wave propagation speed and the damping characteristics of the surrounding medium on the energy attenuation more complicated. However, such coupling effect has not been properly addressed in the past; especially the actual contact strength of the pipe-medium interface cannot be considered. With the large-scale use of plastic pipes and the frequent leakage hazards in China's urbanisation construction, it is urgent to carry out related research to avoid unnecessary costs.
In this paper, the coupling vibration equation of "soil-pipe-fluid" is derived in detail, the acoustic wave propagation characteristic model of the buried pipeline is established, and the shear effect of the medium outside the pipe and the shear effect of the interface between the pipe and the medium on the axisymmetric wave of the fluid dominant are discussed.

Free Motion Equation of fluid-Filled Pipe
In this section, the coupled motion equations of fluid dominant axisymmetric waves in a buried fluid-filled pipes are deduced based on the motion equations of the fluid-filled pipe in vacuum [6]. The soil medium around the pipes is regarded as a homogeneous and isotropic elastic medium which allows both compression wave and shear wave to propagate. According to the current research, the s=1 wave is usually the main carrier of energy in the leakage signal therefore of most interest, and the dynamic damping effect of the pipeline is neglected. Fig.1 shows the cylindrical coordinates of pipes, where u, v, w are the shell displacements in the axial(x), circumferential(θ), and radial(r) directions, respectively. a and h are the pipe radius and the wall thickness respectively and is assumed /1 ha . The internal fluid is assumed to be inviscid, and both the surrounding medium and internal fluid are assumed to be lossless. Fig.1 The coordinate system for a buried fluid-filled pipe For axisymmetric waves (n=0), the rotational motion of the pipe can be neglected, so the circumferential x w u v h a θ displacement and shear stress are both set to zero, free-motion equation of fluid-filled pipe can be described simply according to Donnell-Mushtari shell equation [21] as follow 2  ; () f pa is the internal pressure at the fluid pipe interface.

Motion Equation of Soil Medium
The displacement of soil medium in all directions can be present as ( , , , ) x r t In the column coordinate system, the expansion process of soil medium can be expressed as Here, /0    = . So, the rotating components of soil medium in the axial and radial direction can also be ignored. The rotating component in the direction  is Substituting for r u and According to stresses waves in solids [22], the displacement of a point in the soil medium satisfies the equation of motion Where , mm  are Lame coefficients; m  is the density of the medium;  is the Hamilton differential operator.
Equation of motion shown in Eq. (7) can be expressed as the Bessel equation of cylindrical space outside the pipe Where , GH are the functions of coordinate in the axial (x) and time (t). 01 (), ()

HH
are the Hankel functions of the second kind which describe outgoing waves.
In order to satisfy Eqs. (5) and (9) According to Hook's law, the relationship between stress and strain in the surrounding medium is Where 1 21

Coupling motion equations of pipe-medium
When single axisymmetric waves are considered, the load distribution at the pipe-medium interface is shown in Fig.   2. It assumes that the pipe and medium are in constant contact during the course of movement, the radial displacement of the surrounding medium at the pipe interface is assumed the same as that of the pipe-wall, ur=w. The contact stress of pipe and medium in the radial direction are considered as identical, which expressed as rr  . The frictional stress along the pipe surface equals the shear stress in the axial direction, which expressed as x  . The coupled motion equations given by Eq. (1) can be written as 2 11 13 According to Donnell-Mushtari shell equation [23], the displacement of pipe can be expressed  According to reference [6], for the liquid cannot withstand shear force, the liquid pressure inside the pipe can be expressed directly as the normal displacement of the pipe wall.
Where, is the internal fluid radial wavenumber which can be expressed , here, is the () Substituting Eqs. (18) and (19) Solving Eq. (20) can be obtained 22 The fluid loading term, FL [7] and the surrounding medium loading matrix, SL, are given by

Wave characteristic
For the s=1 wavenumber in buried fluid-filled pipes, 22 1 L kk , so where  stands for the surrounding medium loading and pipe parameters which can be used to evaluate the influence of soil load on the pipe wall displacement.  refers to fluid and pipe parameters which can be used to evaluate the influence of fluid load on the pipe wall displacement. By means of a complex modulus of elasticity Ep (  and  always complex), it is found from Eq. (26) that k1 is always complex indicating the s=1 wave decays as it propagates.
Then  and  are described as the measures of the loading effects of surrounding medium and fluid on the pipe wall respectively.
By Eq. (27),  can be obtained directly, but  which is related to the unknown wavenumber k1 cannot be solved directly. When the pipe is placed in a different medium, the equations can be simplified by boundary conditions.

Lubrication Contact of Pipe-Medium
On the condition of lubrication contact, the contact coefficient 0 It can be seen from Eq. (28) that the propagation of k1 wave will be delayed as it propagates caused by the effect of the pipe wall (i.e., a complex ) and additional damping of the surrounding medium (i.e., a complex ), although there is no frictional damping between pipe and surrounding medium. (1 ) 1 SL is a function of the complex k1, is a complex value. So, s1 wave attenuation is attributed to both material losses along the pipe wall (i.e., a complex ) and radiation losses due to the added damping of the surrounding medium (i.e., a complex ).

Pie in Air
For an air medium, the loading effects of air on the pipe wall can be neglected, the contact coefficient is considered zero. Then, T=0, SL=0, and The Neldes-Mead method [24] is used to solve Eq. (31). In the optimization progress, the termination condition is set as Where, n is the iterations, , re im kk are the centre of the simplex in the current step,  is the tolerance.
In the calculation process, the derivative of Bessel function can be deal based on its property. 01 01 0 2 0 For the arguments of the Bessel or Hankel functions are derived from the square root, it is important to choose the sign of the root. The method to choose the sign of the root in reference [18] is be used. If the real part is larger, the partial wave can be considered homogeneous and must propagate away from the shell, so the positive square root is chosen. On the other hand, if the imaginary part is larger, the partial wave can be considered inhomogeneous and must decay away from the shell, so the negative square root is selected.
The material properties of the fluid, pipe and surrounding medium are shown in Table 1. Considering the efficiency and convergence, the wavenumbers are calculated up to 1 Hz. The thickness/radius ratio of pipe is 0.125, and the plate compressional wave speed is 1725 m/s.  In order to eliminate the effects of surrounding medium and pipe interface friction, the "lubrication contact" assume is used here. Fig.3 gives the wavenumber for the soil with different shear modulus. It can be seen that from Fig.3(a), as the previous theoretical analysis, the effect of the pipe and surrounding medium is used to substantially increase the real part of the s=1 wavenumber from the free-field value kf . This solution is similar to the reference [18]. As the shear effect of the surrounding medium increased, the real part of the wavenumber gradually decreased. The overall loading influences of the surrounding medium are to add mass to the pipe wall, but the shear effect is to add stiffness, and additional stiffness increases with the increase of the shear effect of surrounding medium.

Shear effects at the pipe/medium interface
The shear effects at the pipe/medium interface exist in buried pipe systems, but it is normally ignored because the coupling relationship between pipes and surrounding medium is not clear. Fig.7 and Fig.8 show the wavenumber for s=1 wave with different contact coefficients. In calculation, the shear wave speed of surrounding medium is set to 300 m/s. As shown in Fig.7, the real part of wavenumber for s=1 wave much larger than the free-field value kf with different contact coefficients. Whether the shear effect of the interface between surrounding medium and pipe wall has a great influence on the real part at high frequency(kfa>0.09). Once these shear effects are included, the influence of the contact coefficient on the real part is mainly reflected in the middle frequency band ( (0.02,0.09) f ka  ). In Fig.8, it can be seen that with the consideration of the shear effect of the interface, more waves can be radiated into the surrounding medium.
The attenuation increased as the frequency increased, and has a local peak in the middle frequency band. The frequency corresponding to these peaks increased with the contact coefficient increased. The local peak of attenuation is caused by the resonance between the additional mass of the surrounding medium and the pipe wall, resulting in a dramatic increase in the wave intensity of radiation to the surrounding medium. At the high frequency, the friction force between surrounding medium and pipe wall seems not enough for surrounding medium to vibrate with the pipe wall, so the results of attenuation gradually approaching the value in the case of lubrication contact.
In engineering applications, the actual wave number and attenuation should be obtained according to the actual pipe-soil coupling situation, and the delay estimation should be carried out according to the leakage signal dispersion in order to obtain the accurate leakage location results. Medium and low-frequency signal is generally be used, and the contact coefficient which has a certain influence on the propagation characteristics of the wave in this frequency band should not be ignored.

Field test
This section presents some numerical results of wavenumber from actual plastics water pipes. Two field tests were carried out for pipes in different buried environments. The details of the experimental setup and analysis for test 1 can be found in Ref. [17,20]. The test 2 selected the water supply network in Southwest Jiaotong University, the sensor is set in the pipe well, and fire hydrant discharge signal as leakage signal. The distance from pipe well to hydrant discharge is 12 m. The soil and pipe coupling parameters cannot be obtained in both tests, the theoretical calculation and experimental comparison are not presented in this paper.   Fig.9 Real part of measured wavenumber Fig.10 Loss of measured wavenumber

Conclusion
Axisymmetric waves in thin-walled fluid-filled plastic pipe surrounded by an infinite elastic medium which can sustain both longitudinal and shear waves have been studied. The contact coefficient has been introduced to describe the contact strength of pipe and surrounding medium. Then a general expression for the fluid-dominated wavenumber has been presented in buried fluid-filled plastic pipe.
For axisymmetric waves, the fluid loading dominates the vibration of the plastic pipe wall. The overall loading effects of the surrounding medium are to add mass to the pipe wall, but the shear effect is to add stiffness, which increases with the shear effect of surrounding medium. The shear effect of the surrounding medium also influences the attenuation of the wave.
The added mass of the surrounding medium will resonate with the pipe wall at a specific frequency under the shear effects at the pipe/medium interface, resulting in the change of the propagation characteristics of the wave near the frequency. At higher frequency, the influence of shear effects at the interface on the propagation characteristics is not obvious. The wavenumber can be solved by compact contact theory and the attenuation will approximate the lubricating contact state with the frequency increases.

Availability of data and materials
This paper includes all the data and materials which can be shared.

Competing interests
Not applicable.

Funding
National Nature Science Foundation of China (No.11774378).

Authors' contributions
Ping Lu completed the equation derivation, calculation and writing of the paper.
Xiaozhen Sheng guided the research ideas and results of the paper.
Yan Gao proposed the research idea of this article.
Ruichen Wang revised the language of the paper.