Evaluation of Seismic Protection Methods for Buried Fuel Pipelines Subjected to Fault Rupture

Buried steel fuel pipelines are critical lifelines for the society and the economy but are very vulnerable to earthquake-induced ground failure. Traversing seismic areas inevitably results in several pipe-fault crossings. Fault rupture forces a buried pipeline to undergo deformations that could be substantial and heavily endanger its integrity. Due to the grave consequences of a potential pipe failure, mitigating measures are commonly applied at pipe-fault crossings to reduce the effects of a potential fault activation. In this paper, a comparative review of several measures that are used in practice or have been proposed in the technical literature is presented. Numerical analyses are then carried out to compare the effectiveness of commonly used measures and to extract conclusions regarding their applicability. Results indicate that the most efficient among the evaluated measures are pipe placement within culverts and use of flexible joints. Trench backfilling with pumice is a moderately effective measure in terms of pipe protection, while steel grade upgrade, wall thickness increase, and pipe wrapping with geotextile are found to be insufficient protection methods.

section_content: modes in such cases are local buckling of the pipe wall and tensile fracture. The former is caused by compression, while the latter by the concentration of tensile strains that mainly endanger the integrity of girth welds between the adjacent pipe segments. The pipeline structural response due to faulting depends on the fault type. Namely, in case of normal faulting, the soil rupture incline is downwards and the pipe is subjected to bending and tension (Figure 2A ), while in case of reverse faulting, the soil rupture incline is upwards (both soil blocks are heavily compressed) and the pipe undergoes bending and mainly compression (Figure 2B ). Finally, in case of strike-slip faulting, the fault plane is almost vertical and relative block slip take places place along the rupture plane, consequently the pipe is subjected principally to bending and to a smaller extent to either tension or compression (Figure 2C ), depending on the pipeline-fault crossing angle (Abdoun et al., 2009). 

The assessment of pipe behavior due to fault rupture has been of great concern for researchers, engineers, and pipeline operators. Analytical and numerical approaches have been developed for the assessment of the pipe strain-state caused by faulting. Newmark and Hall (1975) were the pioneers in this ongoing effort, presenting a pertinent simplified analytical methodology, considering the pipe as a long cable subjected to relative small displacements caused by a planar fault dividing the ground into two blocks. Then, several researchers extended this analytical procedure by taking into account material non-linearity and non-symmetric conditions regarding the pipe deformation (Kennedy et al., 1977; Kennedy and Kincaid, 1983 ; Wang and Yeh, 1985; Wang and Wang, 1995; Takada et al., 2001; Karamitros et al., 2007 Karamitros et al.,, 2011;; Trifonov and Cherniy, 2010, 2012) and the pipe service loads (Zhang et al., 2016b). 

The analytical solutions remain useful for the rapid assessment of the pipe behavior in the preliminary design stage. However, securing the pipe resilience necessitates a more rigorous numerical approach. The finite element method (FEM) is used in order to properly take into account the geometrical non-linearity of the problem, as well as the non-linearity of the pipeline material and the surrounding soil. Regarding the numerical modeling of the pipe-fault crossing, two approaches are available. 

The first model is that of a Winkler beam, where the pipe is meshed into beam-type finite elements and the surrounding soil with non-linear translational springs in four directions that model the soil resistance due to the pipe movement in the trench (axial, transverse horizontal, and vertical upward and downward). This model is characterized by sufficient reliability, reduced computational cost, and easy implementation for design purposes. However, local buckling and tensile fracture are not directly "captured" by the numerical analysis and must instead be identified by comparing developing strains to code-based limits. This model has been used by researchers [e. g., Ogawa et al. (2004), Joshi et al. (2011), Chaudhary et al. (2013), Uckan et al. (2015), and Liu et al. (2016) ] and is also recommended by pertinent structural codes [e. g., American Lifelines Alliance (2005), CEN (2006), and Indian Institute of Technology Kanpur and Gujarat State Disaster Management Authority (2007) ]. 

The second model is the continuum one, where the pipe is modeled with shell elements, the surrounding soil with 3D solid elements, and the pipe-soil interface with contact elements [e. g., Kokavessis and Anagnostidis (2006), Odina and Tan (2009), Vazouras et al. (2010 Vazouras et al. (, 2012 Vazouras et al. (, 2015)), Zhang et al. (2014, 2016a ), and Trifonov (2015) ]. This model allows the direct investigation of local buckling occurrence and the consideration of trench geometry and detailed fault movement, as well as the use of more detailed soil material laws. Nevertheless, it is associated with significant computational cost (e. g., increased number of degreesof-freedom, convergence difficulties due to contact elements), and therefore its application for design purposes is very limited. 

Pipeline design against permanent ground displacements is carried out in strain terms, considering that the problem is displacement-controlled and taking advantage of pipeline steel ductility. Thus, pertinent structural codes provide upper limits for tensile and compressive strains in order to avoid tensile fracture and local buckling, respectively [e. g., American Lifelines Alliance (2005 ), CEN (2006 ), Canadian Standards Association (2007), and Indian Institute of Technology Kanpur and Gujarat State Disaster Management Authority (2007) ]. Useful information regarding strain-based design and code-based strain limits can indicatively be found in Gresnigt (1986), Mohr (2003), and Liu et al. (2009). The operable strain limits, after (American Lifelines Alliance, 2005) neglecting the internal pressure as a less favorable situation, are adopted within this study and in particular the tensile strain (εt) limit:

and the compressive strain (εc) limit:

where D is the pipe external diameter, t is the pipe wall thickness, and Dmin is the minimal pipe diameter due to possible ovalization. The response of a pipe subjected to fault rupture is highly affected by the pipeline-fault crossing angle (Ha et al., 2010), with respect to the relation of developing pipe axial force and bending moment. In case of β = 90° the pipe behavior is predominantly flexural, while for angle β < 90° the pipe is additionally subjected to tensile and for β > 90° to compressive axial force, leading to pipe elongation and shortening, respectively, as illustrated in Figure 3. 

The potentially severe consequences of faulting on buried steel fuel pipelines in terms of high developing strains have forced engineers to implement preventive measures at pipefault crossings. These measures usually aim either at reducing the pipe-soil friction (e. g., trench backfilling with loose granular soil) or at increasing the pipe strength (e. g., steel grade upgrade). The objective of this paper is to present an overview of such mitigating measures that have been either applied in practice or proposed by researchers. The effectiveness of these measures is then investigated using advanced numerical analysis, and useful remarks are extracted regarding their contribution to pipe protection. For that purpose, a typical buried steel pipeline is considered crossing a strike-slip fault, chosen as a common fault type causing symmetric pipe behavior around the fault trace.

section_content: Pipeline protection against the consequences of faulting necessitates the employment of pertinent measures. It is, as a general rule, advisable to avoid pipeline-fault crossings at the stage of pipeline route selection; however, this is often impossible in seismic regions. In case such crossings are unavoidable, it is recommended to avoid sharp bends that increase constraints to axial displacements and may impose additional forces on the pipeline (CEN, 2006; Gantes and Bouckovalas, 2013). In case there are such restrictions in the route selection procedure that bends are necessary close to the fault crossing, then the anchor length has to be considered regarding the location of bends. Nevertheless, and in addition to the above guidelines, further specific measures are often necessary to minimize the potential of pipeline failure. To provide a comprehensive overview, the measures are categorized depending on the way protection is provided: pipe-soil friction reduction (Friction Reduction Measures), pipe strengthening (Pipe Strengthening Measures), and other measures (Other Measures).
section_name: MiTigaTing MeasUres

section_content: Fault activation forces the pipe to move within the trench, and consequently friction is generated on the pipe-soil interface that results in the development of strains in the pipe wall. Several preventive measures have been proposed, aiming at reducing the soil-induced friction. A detailed description of such measures is provided subsequently and then a summary is presented in Table 1, where each "method" is briefly described, differentiating between methods that have been applied in practice or have been only proposed by researchers. 

• Pipeline wrapping with friction-reducing geotextile (Figure 4A ) aims at reducing pipe-soil friction and increasing the anchor length (Gantes and Bouckovalas, 2013). Monroy-Concha (2013) has carried out experimental tests and found that wrapping the pipe with a double layer geotextile is effective only if the distance between the pipe and the trench wall is smaller than half pipe diameter. • By placing the pipeline within buried prefabricated concrete culverts, which are "sacrificed" during fault rupture, the pipe is retained undeformed (Figure 4B ). The lack of backfill within the culvert practically eliminates the pipe-soil friction. The application of this measure is limited to strike-slip fault crossings. Cost and constructional issues (e. g., fabrication, transportation and installation of culverts, and trench widening) have to be considered. Moreover, pipe placement within culverts, especially for high fault offset magnitude, arises constructional issues that engineers have to bear in mind, such as the required size of culverts, construction and transportation of oversized culverts, extensive excavation for culvert installation, as well as the capability of culverts behaving as they are expected. • Trench backfilling with loose granular soil, mainly pumice, aims at reducing soil resistance and consequently soil-induced friction on the pipe (Gantes and Bouckovalas, 2013). • Construction of a wider trench can contribute toward allowing pipeline deformation to take place over a longer length and therefore avoiding strain concentration (Figure 4C ). 

To evaluate the efficiency of this design approach, it is necessary to consider trench dimensions. The latter cannot be realized in common design procedures, where the beamtype FEM model is implemented. Thus, either research results are necessary for the quantification of the effect of trench widening on soil spring properties (Chaloulos et al., 2015) or using the continuum FEM model is necessary [e. g., Trifonov (2015) ]. • Geocells and geogrids can be used in the trench above the pipeline (Hegde and Sitharam, 2015), or the pipe trench can be backfilled with tire-derived aggregate surrounded by sand (Sim et al., 2012) to reduce pipe strains. • The use of expanded polystyrene (EPS) geofoam as backfill material has also been proposed, taking advantage of EPS' very low weight. EPS has been applied in transportation infrastructure for the protection of culverts and buried pipes against the applied horizontal and vertical forces. Bartlett et al. (2015) have presented experimental and numerical studies, indicating that EPS geofoam backfill can be effective for the protection of buried steel pipes against permanent ground displacements.
section_name: Friction reduction Measures

section_content: There is a second category of measures, aiming at increasing the pipe strength. This task can be achieved by upgrading the steel grade, increasing the pipe wall thickness to reduce the pipe's bending deformation by increasing pipe stiffness (Gantes and Bouckovalas, 2013; Karamanos et al., 2014), or by wrapping the pipeline with composite wraps (Mokhtari and Alavi Nia, 2015; Trifonov and Cherniy, 2016). These measures are briefly listed in Table 2.
section_name: Pipe strengthening Measures

section_content: Seismic countermeasures that cannot be classified either as friction-reducing (Friction Reduction Measures) or pipe strengthening (Pipe Strengthening Measures) are included in this third category. These measures are presented subsequently in detail and are summarized in Table 3. 

• By introducing bellow-type flexible joints between adjacent pipe parts in the fault vicinity, deformations due to fault activation are absorbed as rotations at the joints, and the pipe steel segments remain barely undeformed and consequently unstressed (Melissianos et al., 2016). The behavior of flexible joints in pipes under fault rupture is depicted in Figure 5, indicating that the continuous pipe structural system is modified to "segmented. " Flexible joints of bellow-type have not been used until now in buried high pressure pipes subjected to faulting applications. However, the extensive parametric feasibility study presented in Melissianos et al. (2016) provides encouraging results on the efficiency of this approach. It is noted that the practical application of joints has to be preceded by practical and regulatory provisions regarding several technological issues, such as the pipe-joint welding and the joint protection against corrosion and external damage. • By elevating sufficiently long part of the pipeline on the two sides of fault crossing above the ground, on appropriate supports, the pipe is no longer forced to follow ground movement during fault activation. The efficiency of this approach has been clearly demonstrated in case of the Trans-Alaska pipeline-Denali fault crossing (Honegger and Nyman, 2004 ). 

The solution is applicable in strike-slip fault crossings, and it is the recommended option in case the expected fault offset is significantly high. It has to be noted that a different design procedure (design in operable conditions and against accidental actions) has to be followed for the above-ground segment of the pipeline. • Zhang et al. (2016) have proposed to use a protective device, whose function relies on imposing external hydrostatic pressure on the pipeline in order to avoid local buckling. • By constructing a dog-leg structure, involving a pipe route change with bend of very high radius, the pipeline deformation is allowed to take place within a larger area, especially in case there is great uncertainty on the fault trace location (Besstrashnov and Strom, 2011). The application of this solution may be impeded by environmental considerations, given that a larger area has to be affected due to the route modification. Additionally, increased cost is required due to the need for more extensive excavation. • It has been proposed to create a predefined buckling pattern, consisting of localized deformation of the pipe wall in shape of local buckling at specified sections, in order to control the location and mode of pipe local deformation and increase the fault displacement that can be absorbed (Hasegawa et al., 2014).
section_name: Other Measures

section_content: A typical API5L-X65 high pressure and large diameter natural gas pipeline is considered as the reference pipe (abbreviated as R hereinafter), featuring diameter of D = 914 mm and wall thickness of t = 12. 7 mm. The pipe steel properties are elastic modulus E = 210 GPa, yield stress fy = 448. 5 MPa, yield strain εy = 0. 214%, failure stress fu = 510 MPa, and failure strain εu = 18%. The pipeline is assumed to be corrosion-and defect-free, coated with coal tar, and buried under 1. 30 m of loose cohesionless sand with unit weight γ = 18 kN/m 3 and internal friction angle φ = 36°. The interface angle of pipe-soil friction δ is estimated as:

where f is the coating-dependent factor that relates the internal friction angle φ of the soil to the friction angle at the pipe-soil interface, considered here as f = 0. 90. Keeping in mind that one of the main goals of mitigating measures is the reduction of the soil-induced friction on the pipe, the corresponding maximum axial (Tu) soil friction force is estimated after ALA (American Lifelines Alliance, 2005) as:

where D is the pipe external diameter, α is the adhesion factor [see Appendix B of ALA (American Lifelines Alliance, 2005)], c is the soil cohesion (here taken equal to 0), H is the depth to pipe centerline, and Ko is the coefficient of pressure at rest. Similar relations for the transverse soil springs are also taken from the same code and are not repeated here in the interest of brevity. 

The eventual pipe internal pressure counterbalances the external soil pressure and acts beneficially against local buckling occurrence. Thus, considering a worst-case scenario (e. g., unpressurized pipe due to maintenance), the internal pressure is neglected. The pipeline is intercepted by a strike-slip fault with vertical fault plane, while for the pipe-fault crossing angle β (Figure 1 ) three indicative values are considered, namely β = 75°, β = 90°, and β = 105°. Finally, the code-based strain limits for the pipeline under consideration are obtained from Eqs 1 and 2, equal to 2% for the tensile and 0. 39% for the compressive limit. 

Fuel pipelines are buried below the ground surface for protection mainly against corrosion and accidental or third-party damage. Burial depth in practice ranges from 1. 0 to 3. 0 D at fault crossings. These upper layers of the ground are usually earth fill with inhomogeneous properties, which in case of fault activation may deviate the rupture propagation from the underlying bedrock to the ground surface. In such case, which is very common in nature, uncertainty regarding the exact fault trace location on the ground surface emerges. Seismological, geological, and geotechnical surveys are required to estimate the length over which the fault trace might "appear" on the ground surface. This length may range from a few meters to hundreds of meters on each fault wall. Within the presented study, the length of fault trace uncertainty is assumed to be 100 m on each side of the fault trace. It is also noted that in current engineering practice, fault trace uncertainty is encountered by applying the mitigation measure along this entire length. This approach is adopted also hereinafter by assuming that the fault trace is located in the middle of the pipe modeled length, while alternative trace locations are not examined considering that sufficient pipeline length is modeled on both sides of the fault.
section_name: case sTUDY reference Pipeline

section_content: Numerical modeling of the pipeline-fault crossing is carried out using the commercial FEM software ADINA (ADINA R&D, 2006). The FE model is based on the suggestions of Eurocode 8-Part 4 (CEN, 2006), ALA (American Lifelines Alliance, 2005). A beam-type numerical model is developed, where the pipe is meshed with PIPE elements, which are two-node beamtype elements including supplementary degrees-of-freedom to account for the additional stresses caused by the cross-section ovalization. The pipeline is meshed into 4,000 pipe elements with length 0. 25 m each. Sensitivity analyses have been conducted on the mesh density with respect to accuracy and computational cost efficiency. The modeled pipe length is sufficient in order for attenuation of strains due to fault activation to have fully developed. Hence, boundary conditions at the ends have no effect. Fixity has been applied only on one end for numerical reasons. Furthermore, the numerical model on the use of flexible joints has been validated based on experimental results presented in Melissianos et al. (2017). The pipe strain-state is assessed through stresses and strains that are calculated at integration points of several cross-sections along the pipeline. Strain estimation is essential in order to perform the safety checking, i. e., to check whether the developed strain exceeds the code-based limit for avoiding tensile fracture and/or local buckling. Moreover, sensitivity analyses that were carried out revealed that modeling each meter of pipe with four elements are sufficient to reliably assess the pipe behavior. 

The surrounding soil is modeled with non-linear translational springs, whose properties are estimated according to ALA provisions (American Lifelines Alliance, 2005), which connect pipe nodes to "ground nodes. " It has to be noted that backfill soil properties are assumed to be similar to those of the native soil. This assumption is commonly applied in practice and is based on the construction process of buried pipelines, where the trench is usually backfilled with the soil that was excavated during trenching. An overview of the numerical model is schematically shown in force 1,494 kN/m and yield displacement 114 mm) model the soil resistance to any vertical pipe movement. The fault plane divides the earth crust into two blocks: (i) the moving one is called hanging wall and is modeled by applying the fault movement as imposed displacement on "ground nodes" and (ii) the stable one is called footwall, in which fixities are applied on "ground nodes. " Analyses are performed by considering geometrical and material non-linearity and employing either Newton-Raphson or arc-length (Bathe and Dvorkin, 1983 ) solution algorithms. The number of analysis steps is selected in such a way that numerical convergence is achieved, and at the same time the displacement is applied smoothly in order to be consistent with the quasi-static nature of the problem. Non-seismic and operational loads (e. g., internal pressure, hydraulic actions, etc. ) are not taken into account.
section_name: numerical Modeling

section_content: The case study aims at evaluating the effectiveness of alternative mitigating measures against the consequences of faulting on buried steel pipelines. To that effect, the commonly adopted measures in design are examined, namely pipe wrapping with geotextiles, pipe placement within culverts, pumice backfilling, steel grade upgrade, and wall thickness increase. These measures have been implemented in several pipeline-fault crossings (Besstrashnov and Strom, 2011; Gantes and Bouckovalas, 2013; Chenna et al., 2014). Additionally, numerical results presented in Melissianos et al. (2015 Melissianos et al. (, 2016) ) indicate that the integration of flexible joints is a promising solution for the protection of buried pipes, considering on the one hand that joints are commercial products and on the other hand the absence of significant constructional requirements. Thus, the introduction of joints will also be examined and compared to the above commonly used measures. The examined measures are the following, considering that each is applied along the entire fault trace uncertainty length:

• Pipe Wrapping with Friction-Reducing Geotextile (FR1)

The pipeline is wrapped with friction-reducing geotextile in order to reduce the developing friction on the pipe-soil interface (Table 4 ). The interface angle δ′ of pipe-soil friction is taken equal to half of the one in case of the reference pipeline (R), i. e., δ′ = 0. 5 δ, which is an assumption often adopted in design practice. Thus, the axial spring force is reduced accordingly (Figure 7 ). Numerical implementation of pipe wrapping with friction-reducing geotextile is carried out by modifying the properties of the axial springs only.
section_name: Preventive Measures

section_content: The pipeline segment along the 200-m uncertainty length is placed within concrete culverts without backfill soil (Table 4 ). The lack of backfilling is numerically modeled by removing the soil springs along this particular length. The friction between the pipe bottom and the trench base is neglected as the contact surface is very small.
section_name: • Use of Culverts (FR2)

section_content: The trench is backfilled with pumice (Table 4 ) having unit weight γ = 8 kN/m 3, cohesion c = 0 kPa, and internal friction angle φ = 33°. The soil properties after ALA (American Lifelines Alliance, 2005) are presented in Figure 7. Numerical modeling of the pumice backfill as preventive measure is carried out by recalculating the axial, transverse horizontal (lateral), and upward soil spring properties along the fault trace uncertainty length. Downward soil springs are estimated using the properties of the native soil and are therefore not modified. Additionally, it is assumed that the trench is wide enough so that its boundaries are not to significantly modify the lateral spring properties. Even though this assumption might be considered as a drawback in lateral spring calculation, the effect of trench boundaries is identified in Appendix B of ALA guidelines (American Lifelines Alliance, 2005) but is not quantified. To address this task, there are very recent research efforts that try to quantify the effect of trench dimensions (Chaloulos et al., 2015 (Chaloulos et al.,, 2017)).
section_name: • Pumice Backfill (FR3)

section_content: The steel grade is upgraded to API5L-X80 with properties: elastic modulus E = 210 GPa, yield stress fy = 530 MPa, yield strain εy = 0. 25%, failure stress fu = 621 MPa, and failure strain εu = 18% (Table 4 ). Numerical modeling of steel grade upgrade (ST1) is carried out by modifying the steel properties of the reference pipelines (R), while soil spring properties are not modified (Figure 7 ).
section_name: • Steel Grade Upgrade (ST1)

section_content: The pipe wall thickness is increased from 12. 70 to 19. 05 mm, which is a commercially available thickness for pipes with diameter equal to 914 mm (Table 4 ). Numerical modeling of wall thickness increase (ST2) is performed through the increase of the thickness of the pipe cross-section of the reference pipeline (R), while soil spring properties are not modified (Figure 7 ).
section_name: • Wall Thickness Increase (ST2)

section_content: Hinged bellow-type flexible joints are placed along the pipeline around the fault vicinity, adopting the configuration proposed in Melissianos et al. (2016) and taking also into account the fault trace uncertainty. Hinged bellows are numerically modeled as generic flexible joints (EJMA, 2008; Peng and Peng, 2009) represented by a rotational spring at the center point without modeling the joint length. The relative axial and lateral displacement capabilities are restricted through appropriate numerical constraints. The adopted modeling approach for pipelines with hinged flexible joints under fault has been verified in Melissianos et al. (2017). 

A summary of the reference pipe and the measures under investigation is presented in Table 4 in order to provide a comprehensive overview of the case study. 

Moreover, the soil spring properties for each examined case are illustrated in Figure 7.
section_name: • Flexible Joints (OT1)

section_content: The maximum tensile and compressive strains obtained from the numerical analyses with respect to the normalized fault offset (Δ/D) for crossing angle β = 75° are shown in Figure 8, where additionally the horizontal lines indicate the ALA (American Lifelines Alliance, 2005) strain limits. It is observed that the maximum developing tensile strains εt,max increase with fault offset increase, while the corresponding compressive strains εc,max increase only up to fault offset equal roughly to Δ/D = 1, which is attributed to the fact that for low fault offset the pipe behavior is mainly bending. Then, as the fault offset increases further, the pipe is subjected mainly to tension and consequently compressive strains due to bending are offset. 

Regarding the efficiency of the examined alternative mitigating measures, it is primarily noted that the efficiency of pipe strengthening measures, namely steel grade upgrade (ST1) and wall thickness increase (ST2), is very limited in terms of strain reduction. Pipe wrapping with geotextile (FR1) does not protect the pipe neither against tensile fracture nor against local buckling. Even though the motivation for wrapping the pipeline with geotextile is to reduce the developing friction on the pipe surface, in practice this measure proves insufficient, taking into account that the pipe surface (including the coating) is already relatively smooth and thus the reduction of friction is limited, while transverse soil-pipe interaction remains unaffected. Trench backfilling with pumice (FR3) contributes to strain reduction and may be considered as a moderately effective approach. The most effective measures are the use of culverts (FR2) and the integration of flexible joints (OT1). Comparing these two approaches, it is shown that for this crossing angle: (i) flexible joints are effective only for low fault offset regarding tensile strains because as the fault offset increases the hinged layout of the joint cannot relief the axial strains due to overall tension, (ii) flexible joints are significantly effective in protection against local buckling as for low fault offset already compressive strains are practically eliminated, and (iii) culverts provide the maximum tensile strain reduction and at the same time do not "allow" the development of compression. 

In case the pipe crosses the fault plane perpendicularly (β = 90°), it is predominately subjected to bending. The maximum tensile and compressive strains with respect to the normalized fault offset (Δ/D) are depicted in Figure 9 for all mitigating measures. It is observed again that steel grade upgrade (ST1), wall thickness increase (ST2), and pipe wrapping with geotextile (FR1) are not effective. In particular, the compressive strains exceed the code-based strain limit for very low offset. Moreover, use of culverts (FR2), trench backfilling with pumice (FR3), and the integration of flexible joints (OT1) can affectively protect the pipe by minimizing the potential of tensile failure and/or local buckling through the limitation of the pipe developing strains. 

The design of a pipe-fault crossing with angle β > 90° is not recommended by pertinent structural codes [e. g., Eurocode 8-Part 4 (CEN, 2006), ALA (American Lifelines Alliance, 2005)] and is avoided as much as possible in design in order to prevent the pipeline from being subjected to significant compression. The resulting tensile and compressive strains with respect to the fault offset (Δ/D) are illustrated in Figure 10. It is again observed that steel grade upgrade (ST1), wall thickness increase (ST2), and pipe wrapping with geotextile (FR1) are insufficient solutions. Pumice backfill (FR3) can be seen as a viable solution only for low fault offset magnitudes. On the other hand, pipe placement within culverts (FR2) and integration of flexible joints (OT1) remain the most efficient solutions in terms of reducing the developing pipe strains. 

It is noted that the case of pipe with flexible joints (OT1) differs from the other cases in terms of developing strains, and therefore compressive strains for this case (Figure 10B ) vary with respect to fault offset, and the pertinent strain curve is not strictly increasing. Flexible joints act as "internal hinges" and thus the adjacent pipe parts are relatively "rotated. " For low fault offset the rotation of joints is small, and compression is developed. Then, compression is progressively minimized as joints are rotated until strains are practically eliminated. Compression is again developed for higher fault offset, due to the significant pipe deformation in the trench. However, it is of great importance to point out that in FR2 case the pipe will undergo significant lateral deformation for large fault offset magnitude that can be judged as a state of global instability. This fact is indicatively illustrated in Figure 11 for fault offset Δ/D = 4, where the pipeline transverse horizontal displacement is presented on the vertical axis and the distance from the fault trace on the horizontal axis. It is shown that the pipe lateral displacement is significant, triggering the question whether this displacement can take place within the culverts. Similarly, in case of integrating joints (OT1), large relative rotations are encountered at the joints, which may exceed the joint rotational capacity and/or exhaust the available trench width. The pipe deformation of OT1 case is also presented in Figure 11. This is at the expense of large relative rotations at the flexible joints, which must not exceed the joints' rotational capacity. 

As a final comment on the use of geotextile, it is derived from Figures 8-10 that within a specific range of fault displacement for all crossing angles, a relatively minor strain amplification is detected if friction-reducing geotextile is used. This is attributed to the fact that even though friction springs are softer due to the geotextile and consequently developing axial strains are lower, transverse springs remain unaffected, thus bending action dominates the response, and hence bending strains are higher. Therefore, longitudinal strains, which are the summation of axial and bending strains, are slightly higher in some cases.
section_name: PresenTaTiOn anD DiscUssiOn OF resUlTs

section_content: The protection of buried steel fuel pipelines subjected to seismic fault rupture has been reviewed and investigated. The mitigating measures commonly used in practice and/or investigated in the literature have been presented, briefly identifying their advantages and disadvantages and discussing constructional considerations of each approach. Most such measures aim either at minimizing the soil-induced friction along the pipeline or at strengthening the pipeline. Even though significant experimental, analytical, and numerical research has been conducted over the past decades and until now, pertinent research on mitigating measures can been seen as limited. Thus, the application of these seismic countermeasures is mostly based on engineering judgment and experience. 

Numerical evaluation of the efficiency of alternative mitigating measures, namely pipe wrapping with geotextile, use of culverts, trench backfilling with pumice, steel grade upgrade, pipe wall thickness increase, and integration of flexible joints, has been presented. A beam-type FEM model has been developed, and aspects on the numerical modeling of each measure are highlighted. A typical buried steel pipeline subjected to strike-slip fault rupture is selected as a case study to examine the efficiency of the measures in terms of strain reduction. Three characteristic pipe-fault crossing angles are examined to demonstrate the effect of this critical design parameter on the effectiveness of protection measures. The uncertainty regarding the exact fault trace location has been also considered in the analysis by considering that each measure is applied over a sufficiently long part of the pipeline to cover this uncertainty. The main conclusions of this investigation can be summarized as follows:

• Pipeline wrapping with friction-reducing geotextile is also not an effective mean of protection, as developing strains are similar to those of the unprotected pipe. • The use of prefabricated concrete culverts can decisively reduce the potential of pipe failure, either tensile fracture or local buckling by significantly reducing pipe strains. This is attributed to the lack of backfill within the culverts, where the pipe is placed. The use of culverts is associated to increased cost, as well as several constructional issues. In case the pipe is subjected to bending and tension (pipe-fault crossing angle β < 90°), pipe developing compression is practically zero. Thus, the effectiveness of culverts in such case is decisive in terms that compression is not developed. In case the pipe is subjected to compression (pipe-fault crossing angle β > 90°), the capability of culverts to follow the fault movement may be questioned, as well as the capability of pipe deformation within the culverts. • Trench backfilling with light-weighted pumice is a moderately effective way to protect the pipe, provided that the soil material used in the construction field is well-graded and the backfilling process is carried out following good engineering practice. This approach could be selected in case the use of culverts or the integration of flexible joints is not a viable solution. • Pipe strengthening measures, namely steel grade upgrade and pipe wall thickness increase, cannot provide sufficient strain reduction and their efficiency is limited, regardless the pipefault crossing angle. Results show that these design approaches are only meaningful if the expected magnitude of fault rupture is relatively low. • The integration of hinged bellow-type flexible joints is a very promising alternative for the protection of buried steel pipelines against seismic fault activation for all considered pipefault crossing angles, provided that the technological issues mentioned in Section "Other Measures" are addressed in a satisfactory manner. However, decline of crossing angle leads to tensional rather than bending pipe behavior, which reduces the joints' efficiency. Moreover, the joint rotation capability for high fault offset has to be examined in detail. 

Finally, it is noted that the selection of a mitigating measure has to be carried out within a techno-economic framework. Therefore, considering the most effective measures among the examined ones, designers have to take into account, among other issues: (i) the manufacturing and transportation cost of culverts at the construction site and the need for excavation of greater amount of soil, compared to the case of a typical trench, (ii) the availability and cost of well-graded pumice at remote construction sites and the requirement for special backfilling process, and (iii) the production cost of customized flexible joints and the joints' compliance with the pipe operable conditions. Engineers therefore have to achieve a balance between safety and economy within the specific requirements defined by the pipeline operator and the regulatory authorities. 

aUThOr cOnTriBUTiOns CG supervised the theoretical background and organized the research process. VM carried out the numerical analyses and prepared the results.
section_name: sUMMarY anD cOnclUsiOn

section_content: Part of the present work has been cofinanced by the European Union (European Social Fund-ESF) and Hellenic National Funds through the Operational Program "Education and Lifelong Learning" (NSRF 2007-2013)-Research Funding Program "Aristeia II, " project "ENSSTRAM-Novel Design Concepts for Energy Related Steel Structures using Advanced Materials, " grant number 4916.
section_name: FUnDing

section_name: Conflict of Interest Statement:

section_content: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
section_name: Conflict of Interest Statement: