Abstract
Long-distance pipelines are threatened by a variety of natural geological hazards. A stress monitoring system driven by the strain-stress solution algorithm was proposed; it can achieve real-time maximum axial stress measurement by installing vibrating wire gauges (VWGs) on the surface of the pipe. To verify the effectiveness of the algorithm, a large-scale pipe mechanical loading experiment combined with a finite element model (FEM) was conducted. The results show that VWGs were reliable with a relative error of 1.19%∼7.98% compared with resistance strain gauges (SGs). The FEM was also reliable with a maximum relative error of 4.04% compared with theoretical analysis. When the reasonable combination modes of VWGs were chosen utilizing the least square method, the error of the pipe stress detection algorithm could be controlled within the range of −13.33∼16.66%. This pipeline stress monitoring technology can meet the requirement of 24-hour dynamic monitoring of the underground pipeline’s mechanical state, realizing the early warning of geohazards.
1. Introduction
1.1. Background
Pipelines are key elements of the energy supply chain. The geological conditions along the long-distance oil and gas pipeline are complicated, and the pipelines are subjected to geological hazards with high risks such as landslide, fault, earthquake, and permafrost [1, 2]. The geological hazards may result in pipeline deformation, stress concentration, upheaval buckling, creep, and even rupture which seriously threaten the structural integrity and long-term safe operation of the pipeline [3]. Although the types of geological hazards are different, the impact of natural geohazards on the pipeline is ultimately reflected in the mechanical state and bearing capacity of the pipeline. If the excessive load applied to a pipeline exceeds the pipeline limit, the pipeline may fail. In addition, even though stress concentration is within the pipeline limit, accidents may also occur if extra control or elimination measures are not taken.
The pipe’s stress and strain responses under various kinds of loads can be investigated through analysis of pipe-soil interactions. From the perspective of soil, as a kind of granular material, the soil has complex mechanical properties due to involved physical and chemical processes [4–6]. Many types of constitutive models including the elastoplastic model, superplastic model, and thermo-hydro-mechanical (THM) model were developed to investigate the behaviors of various types of soil such as loess, freezing soil, and saline soil. In recent years, the thermodynamic theory was applied in the development of constitutive models [7]. Differing from the critical state soil mechanics theory, some concepts such as the flow rule, yield function, dissipation potential function, and hardening criterion included in the classical elastoplastic model can be avoided. New concepts of particle entropy, particle temperature, migration coefficient, and potential energy density function were introduced [8]. Bai et al. [9] proposed a new coupled THM model established under the framework of granular thermodynamics. The influencing factors such as the temperature, saturation, microstructure of the soil skeleton, and phase transition between the liquid and solid/gas phases were reflected in the model. Differing from the classical effective principle based on linear elastic porous media, a generalized effective principle considering the impact of the stress path, temperature path, and soil structure was derived.
From the perspective of the pipeline itself, there are two pipe strength-checking methods based on stress or strain presented in pipeline design standards. The stress-based design criteria are based on limit stress criteria, which ensure that the pipe stress or equivalent stress generated by an external load is less than the minimum yield stress (SMYS) of the pipe itself [10]. The strain-based design criteria are based on limit state design and displacement-control load, which allow that pipeline strain to exceed the specified yield strain resulting in plastic deformation occurring in the pipeline; however, the pipeline still meets the operation requirements [11].
The pipeline strain induced by deformation is a direct reflection of the mechanical state of the pipeline [12]. Commonly pipeline strain is monitored to calculate the pipe’s stress component which is the most intuitive parameter to examine the safety and stability state of the pipeline [13, 14]. The high-precision measurement of strain value during pipeline operation will be crucial to its strength checking and safety assessment [15]. Because of a characteristic of concealment and temporal hysteresis in the response of pipeline foundation soils to external disturbance, no matter how many considerations have been made in the engineering design stage, it is still impossible to accurately predict all the complex and changeable engineering geological conditions. There will still be great uncertainties and safety risks in the pipeline operation. However, most pipeline accidents attributed to geological hazards have a gradient process in the accumulation of pipe strain and stress, which gives us the possibility of identifying the risk in the early stage and adopting prevention and mitigation measures promptly. Thus, a comprehensive stress monitoring system is required and mandatory, which could make deformation measurements sufficiently sensitive and prompt the detection of the approach to operational tolerance limits.
Monitoring of structural response under various loading conditions has been applied in other geotechnical engineering fields, such as the monitoring of slope, pile structure, tunnel, and dam. Inadequate parameters may lead to uneconomical or sometimes faulty design. Therefore, the application of performance monitoring systems to engineering promotes the transformation from code-based design to performance-based design. Suhail et al. [16] summarized the application of fiber-optic sensor (FOS) technology in monitoring pile performance including pile-soil friction distributions, strain variation with depth, and load transfer behavior. The advantages of FOSs, such as high accuracy, high-temperature durability, and long life cycle, and shortcomings of strain and temperature discrimination were presented. Yuan et al. [17] utilized transparent soil technology combined with particle image velocimetry (PIV) to study the soil deformation processes and displacement trends around a laterally loaded pile. In the experiment, six electrical resistance strain gauges were attached to the model pile body to measure the deformation of the pile. With the rapid development of urbanization, road building and mining have increased the problem of slope stability. Yang et al. [18] summarized the slope deformation monitoring technologies including distributed fiber-optic strain sensor (DFOSS). It is noted that due to the monitoring results being easy to be interfered with by external factors, more studies should be carried out in monitoring data analysis and theoretical models. Some scholars carried out a sensitivity analysis to investigate the effect of influencing factors such as internal friction angle, ratio of the slope, and cohesion on slope stability. The results can provide a reference for the design of the slope and customized slope stability monitoring program. Moreover, postconstruction monitoring and reconnaissance will provide an important feedback loop to the design process [19–21].
Direct monitoring of pipeline strain could be achieved via several technologies. The most notable ones are in-line inspection (ILI), vibrating wire gauges (VWGs), fiber Bragg grating (FBG), and distributed strain sensing (DSS) with fiber optics [22–25]. The ILI tools provide an important means of monitoring pipelines in the permafrost zone, which may experience thaw settlement, frost heave, upheaval buckling buoyancy, and slope instability. For example, the annual ILI tools run over the first 336 km of the Norman Well Pipeline and Trans-Alaska Pipeline System (TAPS) segments in permafrost terrain [26]. A settlement of about 4.57 m at MDX 200 of TAPS was detected using ILI tool [20]. Pipeline bending strain values can be obtained based on the inertial measurement unit (IMU) assembled in the pipeline ILI tool. This technology was employed by the PetroChina Pipeline Company to identify high thawing settlement risk areas along with China-Russia Crude Oil Pipelines (CRCOPs) [27].
Compared with the ILI tool, which provides only a snapshot in time of the strains in the pipeline [28], fiber-optic cable sensor installed along the pipeline and vibrating wire point gauges directly attached to the pipeline could provide strain measurement with a higher frequency. A fiber-optic geotechnical monitoring system has been used on Sakhalin–Khabarovsk–Vladivostok (SKV) gas pipeline, with a total length of approximately 26 km, which crosses 32 Active Tectonic Faults (ATF) in seismically active areas to track pipeline performance. In addition, a total length of 92 km is monitored with strain and temperature fiber-optic sensor cables for ground movement detection [29]. A ground movement monitoring and leak detection system based on fiber-optic monitoring has been applied on the Peru LNG pipeline [30]. The fiber-optic technology can successfully detect pipeline strain with high accuracy and reliability; however, it is complex and costly from the perspective of installation and data processing. It requires a complex installation geometry that is not paired well with other variables’ detection on other abnormal pipeline behavior such as leak detection, ground movement monitoring, and ROW intrusion.
The VWGs could reflect axial strains according to the functional relationship between the vibration frequency of the metal string and the tension force. It has an advantage in reliability compared with SGs, convenient installation compared with FBG or DSS, and high frequency compared with ILI. At the same time, according to the theory of material mechanics and Hooke’s Law, the maximum axial stress could be acquired by only knowing three axial strains of the pipe section. Thus, the actual stress and strain variation could be obtained by installing the VWGs on the pipeline monitoring sections. Based on the monitoring results and stress or strain criteria, the degree of pipeline damage caused by geohazards could be evaluated quantitatively, and preventative mitigation measures could be adopted before the pipeline accidents occur. The purpose of this paper is to propose a pipeline stress monitoring system driven by the strain-stress solution algorithm. To verify the effectiveness of the algorithm, a large-scale mechanical loading experiment combined with FEM establishment was carried out for the first time.
1.2. Outline of a Pipe Stress Monitoring System
The pipe stress monitoring system is composed of VWGs, data acquisition controllers, and a monitoring center (Figure 1). The distributed acquisition mode is adopted by arranging main controllers and sensors near the monitoring pipe segments. 3∼5 VWGs attached to the pipe section are used to acquire pipe axial strain, and the maximum axial stress could be calculated according to the stress detection algorithm. The data acquisition controllers carry out data acquisition and signal processing. The data is communicated to a PC to reduce the transmission distance of analog quantity and improve the antijamming ability and stability of the whole system. Finally, the stress software installed in the monitoring center performs analysis and early warning realizing the 24-hour monitoring of the pipeline.

2. Materials and Methods
This study facilitates a multisource data coupling analysis for the pipeline’s mechanical state (Figure 2): the pipe mechanical loading experiment, the FEM, the algorithm derived, and the basic theory of material mechanics. In particular, it is the first time that the algorithm is published and the pipe mechanical loading experiment for this purpose is carried out.

(a)

(b)
2.1. Theoretical Basis and Algorithm Derivation
It is assumed that the pipe is in a state of linear elasticity and small deformation. Compared with the inner pressure and bending moment , the lateral force and torque T are smaller and could be neglected. Therefore, it is considered that all measuring points in the pipe are in a state of bidirectional stress and are only affected by the axial stress and the hoop stress . The solution formulas of axial stress , hoop stress , and equivalent stress are listed as follows. Formula (3) is applicable under the condition that the hoop stress is the first principal stress and the axial stress is the third principal stress (compression).
On the monitoring section, 3 VWGs are arranged, and the coordinates of the measuring points are , , and , as depicted in Figure 3. The axis and axis define the plane perpendicular to the pipeline axis. The measured axial strains are , , and . The axial force of the pipe is F, and the bending moment is . The angle between the bending neutral axis and the axis is .

The projection of bending moment on axis and axis is shown in (4). According to the generalized Hooke Law, the axial strains of the three measuring points are shown in (5).where is the hoop stress (Pa), is the internal pressure (Pa), is the outer diameter of the pipe (mm), is the wall thickness of the pipe (mm), is the elastic modulus of the pipe (mm), 2.03 × 1011 Pa and is Poisson’s ratio of the pipe, 0.3.
According to the material mechanics principle, the axial stresses of the three measuring points are shown in the following equation:where I is the axial moment of inertia of the pipe section (mm4), is the cross-sectional area of the pipe (mm2), .
There are three unknown parameters, F, M1, and M2, so three pipe strain measuring points are required. The strains are expressed in the following equations:
According to the equations above, the maximum and minimum stress and strain can be derived by the following equations:
2.2. Large-Scale Pipeline Compression Experiment
A test pipe with a length of 8 m, a diameter of 1219 mm, a wall thickness of 26.4 mm, and a steel grade of X80 was installed on the mechanical loading platform. A driving arm and a confining arm were hinged to the pipe, and a hydraulic oil cylinder provided power to the driving arm with a movement of a constant velocity of 600 mm/min (Figure 4). Each movement was 20 mm, and each load was kept for more than 1 minute. The pipe hoop stress was generated by the constant internal water pressure. The axial stress was generated by the force resulting in tension strain in the outer side and compression strain in the inner side. Before the experiment, the welding of pipe flanges, installation of VWGs, and deployment of data acquisition system were completed. The trial operation included a pipe water pressure test, operation of hydraulic cylinder, data acquisition, and communication system test. The whole experiment took 1 hour and 9 minutes from loading to pipe rupture. Five monitoring sections in the test pipe were determined. For each monitoring section, five VWGs were installed at the 9, 10.5, 12, 13.5, and 3 o’clock positions, and one SG was installed at the 3 o’clock position for comparison (Figure 5).

(a)

(b)

(a)

(b)
2.3. Finite Element Simulation
The corresponding finite element model was established based on the experimental conditions with 332360 units and 381308 points (Figure 6). Five monitoring sections were also set in the same positions consistent with the experiment in the field. The x axis was along the driving arm direction, the y axis was along the vertical direction, and the z axis was along the pipe direction. In the model, zero displacement constraints in x, y, and z directions were applied to the restrained end, and zero displacement constraints in x and y directions were applied to the loading end. The hydraulic oil cylinder thrust applied an acting force on the nodes of the loading end in the form of a “concentrated force,” while the internal water pressure stressed the nodes of the pipe inner wall in the form of a “uniform distribution force.” The internal water pressure remained unchanged, and the force of oil cylinder thrust was the main load. The bolts in the field were ignored, and it is assumed that the driving arm and confining arm were directly connected with the pipe.

(a)

(b)
Under most geological conditions, pipelines are subjected to displacement load. Take frost heave displacement as an example; as illustrated in Figure 7, in the pipeline sections with frost-susceptible foundation soils and favorable hydrothermal conditions for frost heaving, the frost bulb generated by a cold pipeline may push the pipe upward while the adjacent pipes anchored in the frozen ground tend to hold the pipe down, causing bending of the pipe. If the bending strain becomes larger than the ultimate strain of the pipe, bucking or rupture failures of the pipe may happen. The pipe’s loading condition was mainly simplified to a “concentrated force” applied by hydraulic oil cylinder thrust. The additional bending stress can be induced by the applied load in this way, and soil around the pipe was ignored. Although it is inconsistent with the field, the pipe’s mechanical states can be reflected without considering the soil around the pipe at a lower cost.

3. Results
3.1. Distributions of Pipe Displacement and Stress
The variation in hydraulic oil cylinder thrust with time is shown in Figure 8. The oil cylinder provided power to the driving arm with a movement of a constant velocity of 600 mm/min. Each movement was 20 mm, and each load was kept for more than 1 minute. From 350 s to 1005 s, the load kept at 1577 kN, so it was observed that there was an obvious load stabilization stage during the loading process. The distributions of overall displacement, pipe displacement, and pipe stress at 651 s with cylinder thrust of 1577 kN and inner pressure of 10.9 MPa are shown in Figure 9. Under the movement of the driving arm, pipe axial stress values on the outer side were positive and the ones on the inner side were negative, illustrating that the outer side was tensile and the inner side was compressed. The displacements of each part of the pipe were distributed symmetrically along the pipe centerline in the z direction. The distributions of displacement and stress as well as symmetry were all under common sense. At the same time, the distributions of six stress components in the pipe are also given in Figure 10.

(a)

(b)

(a)

(b)

(c)

(d)

(a)

(b)

(c)

(d)

(e)

(f)
3.2. Verification of Results of FEM
According to the nonlinear FEM established, the calculated stress results of five pipe sections are listed in Table 1. All the maximum shear stresses in each section were less than 90% of pipe yield strength, illustrating that the pipe was in a linear elastic state. According to the basic theory of material mechanics, the results of FEM were compared with theoretical analysis (Table 2). The comparison results show the calculated values of the two methods were similar, and the maximum relative error was only 4.04%. Thus, the FEM was reliable and could be used to verify the effectiveness of the algorithm combined with the experiment results.
3.3. Comparison of Results between VWGs and SGs
The monitoring results of the VWGs installed at the 3 o’clock position in four pipe sections were compared with SGs which were located at the same positions (Figure 11). Since the monitoring results of section 1 were unreasonable, they were not used for analysis here. The comparison results show that there was a similar variation trend with a maximum relative error of 7.98% and a minimum of 1.19%. This illustrated that the VWGs were reliable in pipe axial strain monitoring with high relevance.

(a)

(b)

(c)

(d)
3.4. Optimization in Combination Modes of VWGs
According to the pipeline stress detection algorithm derived above, the maximum axial stress of the pipe section can be obtained by only knowing three axial strains. In this experiment, there were 5 VWGs arranged on each pipe section, so there were 10 different combination modes in total if 3 were selected for one combination. To evaluate the reliability of the algorithm and optimize the combination modes of VWGs, the results of 10 different combination modes in 5 monitoring sections were substituted into the algorithm to compute the stress for comparison with the FEM results. Take pipe section 1 as an example; the comparison of results between stress algorithm and FEM are listed in Table 3. Before the screening, the relative error of the shear stress was −13.33∼48.11%.
As can be seen from Figure 12, when the combination modes from the first to the sixth were adopted, the relative error of the maximum axial stress ranged from −7.79 to 26.74%. Among all sections, the error of section 4 was relatively small, ranging from −0.48 to 2.82%, and the minimum error was −0.09%. When the third, fifth, sixth, eighth, ninth, and tenth combination modes were adopted, the relative error of the minimum axial stress ranged from −6.00 to 31.48%, and the error of section 4 was relatively small, with a minimum value of −6.00%. The relative error of the maximum shear stress ranged from −13.26 to 14.59% when the second to the sixth combination modes were used. Among all sections, the error of section 5 was relatively small, ranging from −1.46 to 0.72%, and the minimum error was 0.34%. It is important to note that in the stress analysis of a pipeline, the strength check is mainly judged by the maximum shear stress, and the maximum and minimum axial stress are only the intermediate results of the pipeline stress detection algorithm. Thus, the optimal combination modes for the maximum shear stress should be adopted to evaluate pipe safety. Based on the comparison of VWGs monitoring results and FEM simulation results, the error of gauging point S3 was large. The reason is that gauging point S3 is located near the neutral axis. Under the same load condition, the monitoring strain value of this point is much smaller than those of other points (S1, S2, S4, S5), which easily causes large errors. Therefore, the curve fluctuates up and down when combination modes including gauging point S3 are plotted.

(a)

(b)

(c)
4. Discussion
According to the basic theory of material mechanics, under the condition of linear elasticity and small deformation, the strain obeys the assumption of a plane section, so the results of measuring points S1∼S5 in the same section should present a straight line in the ideal state. The monitoring results of five sections could be drawn in one figure (Figure 13(a)). If a measuring point deviated too far from the straight line, this indicated that it was not reliable and should be excluded. In addition, when the force applied on multiple sections was similar, the measuring points of each section should be approximately collinear in the ideal state (Figure 13(b)). If the straight line of one section deviated too far, this indicated that the monitoring result of this section was not reliable and should be excluded.

(a)

(b)
Take the strain results of pipe section 1 as an example. Firstly, the fitting line was drawn utilizing the least square method (Figure 14). Then, the distance from each point to the line was squared, and the point with the greatest deviation was gotten rid of. It was S4 in this case. After the S4 was excluded, the initial ten combination modes could be reduced to four combination modes. Finally, the results of the remaining points, S1, S2, S3, and S5, were substituted into the algorithm to compute the stress values of the four combination modes. The relative error of maximum shear stress could be reduced from −13.33∼48.11% to −13.33∼2.26%. With the same method, the results of S2 in section 2 and S3 in sections 3∼5 deviated from the fitting line and were not substituted into the algorithm in calculating the stress. The relative errors of the maximum shear stress in sections 2∼5 were further reduced to 11.82∼16.66%, 6.94∼11.66%, −9.32∼6.9%, and −1.46∼8.99%, respectively. The overall relative error of the five sections was determined as −13.33∼16.66% which was similar to −13.26∼14.59% presented initially above.

In addition, when the error was negative, it indicated that the calculation result of the algorithm was smaller compared with the real condition. On the contrary, a positive value means that it was larger. The former will lead to pipe strength warning with a delay, which may cause engineering accidents. The latter will cause the warning to be too sensitive, resulting in wasting the bearing capacity of the pipe. However, the former consequence was more serious.
5. Conclusions
(1)According to the theory of material mechanics and Hook’s Law, a pipe stress detection algorithm, which could compute the maximum axial stress by only knowing three axial strains of the pipe section, was derived.(2)To verify the effectiveness of the stress detection algorithm, a large-scale pipeline compression experiment was conducted on a mechanical loading platform, and a corresponding FEM was established. The VWGs were reliable in pipe axial strain monitoring compared with SGs attached at the same location with a relative error of 1.19%∼7.98%. The FEM was also credible compared with theoretical analysis with a maximum relative error of 4.04%. Thus, the algorithm could be verified by the experiment combined with the FEM.(3)The least square method could be utilized in excluding unreliable monitoring points of VWGs, and combination mode numbers could be reduced from 10 to 4. When the reasonable combination modes of VWGs were chosen, the error of the pipeline stress detection algorithm could be controlled within a range of −13.33∼16.66% with enough accuracy for engineering application.(4)Limited by measuring accuracy of VWGs, the measuring point at pipe crown is not recommended to be arranged because this measuring point is easy to cause a large measuring error.Data Availability
The data used to support the findings of this study are available from the corresponding author upon request.
Conflicts of Interest
The authors declare no conflicts of interest.
Acknowledgments
This research was funded by the Chinese Academy of Sciences (CAS) Strategic Priority Research Program (Grant no. XDA20100103).