Methodology for Resilience Assessment of Oil Pipeline Network System Exposed to Earthquake

: The oil pipeline network system (OPNS) is an essential part of the critical infrastructure networks (CINs), and is vulnerable to earthquakes. Assessing and enhancing the resilience of the OPNS can improve its capability to cope with earthquakes or to recover the system’s performance quickly after the disturbance. This study deﬁnes the concept of OPNS resilience in the resistive ability, the adaptive ability, and the recovery ability. Then, the quantitative resilience assessment model is established considering the earthquake intensities, the role of safety barriers, the time-variant reliability, and the importance coefﬁcient of each subsystem via a Monte Carlo simulation. Combining the model with GIS technology, a new methodology to evaluate OPNS resilience is proposed, and the resilience partition technology platform is developed, which can visualize the results of the resilience assessment. Finally, a case study is implemented to demonstrate the developed methodology, and a discussion is provided to identify the sensitive variables. The proposed resilience methodology can provide a framework for the probabilistic resilience assessment of OPNS, and could be expanded to other lifeline network systems.


Introduction
The operation and development of cities depend on a continuous energy supply, such as oil and gas, and pipeline transportation provides convenient routes for them [1]. The oil pipeline network system (OPNS), as a part of critical infrastructure networks (CINs), is easily affected by large-scale geological natural disasters, such as earthquakes [2,3]. Once a high-intensity earthquake occurs, it will cause serious damage to the structure and may trigger subsequent technological disasters to increase the severity of accident loss [4]. In the 1987 Ecuador earthquake (Mw 6.9), almost 40 km of the oil pipeline was destroyed, and more than 800 million dollars were lost [5]. The 1994 Northridge earthquake in Colorado ruptured about 50 km of pipeline and caused more than 3500 barrels of oil to be released. The leaked oil triggered more than 100 cases of fires and explosions, resulting in 30 billion dollars of damage and causing severe pollution in the local environment [6].
Currently, the seismic safety research on oil pipelines mainly focuses on seismic design [7,8], vulnerability assessment [9,10], consequence mitigation [11], and mechanical response analysis [12,13]. Makhoul et al. [14] studied the seismic damage results of buried pipelines by using a series of typical fragility models, and proposed a resilient earthquake preparedness strategy for Byblos' infrastructure network. Germoso [15] modeled the seismic responses of buried pipelines and developed a novel dynamic analytic approach by incorporating the nonlinear soil interaction, which enables the analysis of seismic pipeline safety towards an ideal design for its performance and mechanical reliability. Shabarchin [16] quantified the failure probability of oil pipelines owing to wave propagation

•
The resistive ability-refers to the power of the system to resist disaster damage and maintain structural stability [36].

•
The adaptive ability-refers to how a system adapts to the process of disaster evolution and ensures its function is stable [37]. • The recovery ability-the ability of system performance to recover to a certain level after a disaster under the action of recovery measures [38].

Definition of OPNS Resilience
The OPNS is a network system composed of multiple pipeline subsystems and, thus, the resilience of OPNS is affected by each subsystem, and can be quantified by their resilience. However, the importance of each pipe is not the same in the network and these differences must be fully considered to quantify the system resilience accurately.
The importance coefficient of pipe segment (Ic), defined as the influence degree of pipe segment on the overall performance of OPNS, is introduced to optimize the system's resilience. The Ic of the pipeline is a metric that completely takes into account the significant relationship between subsystems and systems, focusing primarily on the state before the disaster, and which will not change with the occurrence of disasters. The larger the Ic is, the higher the influence degree of the pipe is.
As shown in Figure 1, it is assumed that the total flow rate of the pipe network in regular operation is Q and that the flow rate of each pipe section is q i (i = 1,2 . . . n). The Ic of each pipe segment can be calculated according to Equation (1), as follows: a system to prepare and adapt to the changing conditions and to withstand and recover from disruptive events [33][34][35]. The above definition shows three abilities of system resilience, as follows:  The resistive ability-refers to the power of the system to resist disaster damage and maintain structural stability [36].  The adaptive ability-refers to how a system adapts to the process of disaster evolution and ensures its function is stable [37].  The recovery ability-the ability of system performance to recover to a certain level after a disaster under the action of recovery measures [38].

Definition of OPNS Resilience
The OPNS is a network system composed of multiple pipeline subsystems and, thus, the resilience of OPNS is affected by each subsystem, and can be quantified by their resilience. However, the importance of each pipe is not the same in the network and these differences must be fully considered to quantify the system resilience accurately.
The importance coefficient of pipe segment (Ic), defined as the influence degree of pipe segment on the overall performance of OPNS, is introduced to optimize the system's resilience. The Ic of the pipeline is a metric that completely takes into account the significant relationship between subsystems and systems, focusing primarily on the state before the disaster, and which will not change with the occurrence of disasters. The larger the Ic is, the higher the influence degree of the pipe is.
As shown in Figure 1, it is assumed that the total flow rate of the pipe network in regular operation is Q and that the flow rate of each pipe section is qi (i = 1,2…n). The Ic of each pipe segment can be calculated according to Equation (1), as follows: 1 1 where R is the resilience of the OPNS, and Res is the pipeline subsystem resilience.
Ic i ×Res (2) where R is the resilience of the OPNS, and Res is the pipeline subsystem resilience. The measurement of system performance is assumed to be system reliability for convenience. As shown in Figure 3, the pipeline reliability varies over time. Before the occurrence of an earthquake, the initial state reliability of pipes is the maximum value P 0 , and this decreases immediately to P 1 due to the disaster (earthquake) occurring at time t 0 . With the evolution of primary accident scenarios and the adaption of the pipeline, reliability will further reduce to the minimum value P 2 at time t 1 . In addition, the recovery strategy is prepared during the time from t 0 to t 1 and applied at time t 1 , and the pipeline reliability begins to recover at time t 2 and fully recovers to P 3 at time t 3 .

Resilience Metrics
The general approach to quantifying resilience is shown in Equation (3), as follows: where P(t) is the performance of a system, t 0 is the time when the disaster event occurs, and t f is the time to complete restoration of system performance. The numerator of Equation (3) indicates the accumulation of system performance P(t) between t 0 (disruption) and t 3 (fully recovered). The denominator represents the accumulation of initial system performance P(t 0 ) between t 0 and t 3 . Although the resilience metric is illustrated by the case of system  The measurement of system performance is assumed to be system reliability for convenience. As shown in Figure 3, the pipeline reliability varies over time. Before the occurrence of an earthquake, the initial state reliability of pipes is the maximum value P0, and this decreases immediately to P1 due to the disaster (earthquake) occurring at time t0. With the evolution of primary accident scenarios and the adaption of the pipeline, reliability will further reduce to the minimum value P2 at time t1. In addition, the recovery strategy is prepared during the time from t0 to t1 and applied at time t1, and the pipeline reliability begins to recover at time t2 and fully recovers to P3 at time t3.

Resilience Metrics
The general approach to quantifying resilience is shown in Equation (3), as follows:   The measurement of system performance is assumed to be system reliability for convenience. As shown in Figure 3, the pipeline reliability varies over time. Before the occurrence of an earthquake, the initial state reliability of pipes is the maximum value P0, and this decreases immediately to P1 due to the disaster (earthquake) occurring at time t0. With the evolution of primary accident scenarios and the adaption of the pipeline, reliability will further reduce to the minimum value P2 at time t1. In addition, the recovery strategy is prepared during the time from t0 to t1 and applied at time t1, and the pipeline reliability begins to recover at time t2 and fully recovers to P3 at time t3.

Resilience Metrics
The general approach to quantifying resilience is shown in Equation (3), as follows: However, to assess the resilience of the OPNS exposed to an earthquake, a more robust approach is needed to define system performance on a resilience metric. Combined with the logical relation of resilience in Section 2.2, the OPNS resilience (Res) is quantified as the function among resistive, adaptive, and recovery abilities of the system, as follows in Equation (4).
where Re is the quantification of resistive ability, which is defined as the ability of the pipe section to resist external disaster (earthquake) disturbance to a certain extent to ensure normal function. It can be estimated as the ratio of residual reliability to the initial reliability at the time of the disruptive event [39] as follows in Equation (5): Ad is the quantification of adaptive ability, which is defined as the power of a pipe segment to effectively adapt to the evolution process of primary accident scenarios caused by disasters (earthquakes) to avoid the functional reliability reduction. It can be quantified as the ratio of the lowest reliability to residual reliability during the period before recovery action starts, as follows in Equation (6): Rec is the quantification of recovery ability, which is defined as the ability of the pipe section to reach the minimum performance level after emergency repair quickly, and can be quantified as the slope of the recovery and estimated using the following Equation (7): Therefore, based on Figure 3 and Equations (3)-(7), the resilience metrics can be transformed as follows in Equation (8):

Quantification of Resistive Ability
According to the definition of resistive ability in Section 2.3, the R 1 can be quantified through Equation (9), as follows: where P damage is the probability of pipeline damage under disasters. For a pipeline, it can be considered that the damage to the pipe is random and independent along its length, and that it obeys the Poisson distribution. Therefore, the P damage can be calculated by Equation (10), as follows: where L is the length of a pipe segment, and Rr is the damage rate of the pipe exposed to disaster, which refers to the number of damaged areas per kilometer of a pipeline under a certain kind of disaster. For an earthquake, it contains pipe diameter, materials, site conditions, soil liquefaction degree, and other factors, and it can be obtained as follows in Equations (11) and (12) [40]: where C p is the influence coefficient of the pipeline, C d is the influence coefficient of diameter, C g is the influence coefficient of geological topography, C l is the influence coefficient of soil liquefaction, R 0 is the standard damage rate of the pipeline, PGA is the peak ground motion acceleration, PGV is the peak ground motion velocity, and the quantification of these coefficients can be found in [40].

Quantification of Adaptive Ability
In the adaptation stage, the pipeline structure damage caused by a disaster (earthquake) can lead to oil leakage and may occur in the subsequent accident, such as a fire and explosion, which will further reduce the reliability of a system to the minimum level (P 2 in Figure 3) before it finally reaches a stable state before restoration.
The involvement of safety barriers can prevent and mitigate these subsequent scenarios, which will increase R 2 and reduce the pipeline reliability loss, thus, improving the adaptive ability and restraining the reliability decline. A technical safety barrier can be generically defined as a physical means to prevent, mitigate, or control undesired events or accidents [41,42]. This paper classified technical safety barriers into the following two categories [43]: • Active barriers, which require automatic/external activation [44], such as PSVs and emergency shutdowns; • Passive barriers which are already in place and do not require external activation [45] (e.g., PFP and flexible connections).
The quantification of safety barrier performance mainly considers two indicators, namely reliability and effectiveness [46], and, thus, the probability of failure on demand PFD and effectiveness η [47] were used to characterize the performance of different safety barriers in various scenarios in this article, as is shown in the following Equation (13): where PFD 0 is the baseline value for barrier probability of failure on demand, η 0 is the baseline value for barrier effectiveness, φ is the modification factor that is used to determine a tailored value of PFD, and η, which can be interpreted as the likelihood that barrier systems are impaired or damaged by natural hazards, hence, a higher value (i.e., closer to 1) indicates a higher probability that the barrier will fail in providing a successful protection action. Table 1 shows the common types of safety barriers in OPNS [47][48][49]. Then, P 2 can be calculated as follows in Equation (14): where P e is the probability of subsequent accidents triggered by disaster, and the event tree representing the evolution process of disaster must be built to calculate P e . For the OPNS, the parameters of each part of the segment and the barrier resettlement are different and it is necessary to create an event tree. Then, the probability can be quantified using the following Equation (15): where P safe is the probability that the final state is safe in the development of the event tree, and x is the number of the safe state in the evolution of the event tree. Figure 4 shows the event tree as an example. The logical operators described as gates can be obtained in Figure 5, adapted from a previous study [50].
where Psafe is the probability that the final state is safe in the development of the event tree, and x is the number of the safe state in the evolution of the event tree. Figure 4 shows the event tree as an example.  The logical operators described as gates can be obtained in Figure 5, adapted from a previous study [50].

Quantification of Recovery Ability
Many factors, such as daily usage of hazardous substances, equipment value, the degree of change in system performance, and the cost of recovery need to be considered to quantify the recovery ability. Due to the duration of an earthquake disaster being concise and aftershocks making the post-earthquake recovery mechanism very complicated, it is  where Psafe is the probability that the final state is safe in the development of the event tree, and x is the number of the safe state in the evolution of the event tree. Figure 4 shows the event tree as an example.  The logical operators described as gates can be obtained in Figure 5, adapted from a previous study [50].

Quantification of Recovery Ability
Many factors, such as daily usage of hazardous substances, equipment value, the degree of change in system performance, and the cost of recovery need to be considered to quantify the recovery ability. Due to the duration of an earthquake disaster being concise and aftershocks making the post-earthquake recovery mechanism very complicated, it is

Quantification of Recovery Ability
Many factors, such as daily usage of hazardous substances, equipment value, the degree of change in system performance, and the cost of recovery need to be considered to quantify the recovery ability. Due to the duration of an earthquake disaster being concise and aftershocks making the post-earthquake recovery mechanism very complicated, it is reasonable to simplify the recovery ability in the cross-regional level energy transmission system [51]. By assuming the damage of each pipe section in the OPNS subject to the Poisson distribution with the parameter Rr along the pipeline [52], the distance D i between two consecutive rupture points in each pipe section system in the OPNS exposed to earthquake can be calculated according to Equation (16), as follows: where UD is a random variable subject to the uniform distribution of 0-1, and D i are randomly generated according to Equation (15). For a single pipeline, when m-1 is the maximum number of rupture points. Typically, the oil station at both ends of the pipeline is equipped with corresponding emergency teams to repair the damaged parts of the pipe in time after the disaster occurs. The repair performance of the pipe is defined as the ratio of the number of rupture points repaired per unit of time to the total number of them. Assuming that all emergency teams will immediately fix the damaged parts when the pipe is damaged and the rescue time needed to repair each rupture point is subject to N~(6h,3h) according to the HAZUS-MHMR4 [53], the performance recovery curve of each pipe section at different times can be obtained by a Monte Carlo simulation, which is shown in Figure 6.
Based on Figure 6, the time begins to recover (t 2 in Figure 3) and the time to fully recover (t 3 in Figure 3) can be obtained as the average time through the simulation, and a shorter t 3 indicates that more rupture points can be repaired per unit time and the pipe section can reach the minimum operation level earlier. emergency teams to repair the damaged parts of the pipe in time after the disaster occurs. The repair performance of the pipe is defined as the ratio of the number of rupture points repaired per unit of time to the total number of them. Assuming that all emergency teams will immediately fix the damaged parts when the pipe is damaged and the rescue time needed to repair each rupture point is subject to N~(6h,3h) according to the HAZUS-MHMR4 [53], the performance recovery curve of each pipe section at different times can be obtained by a Monte Carlo simulation, which is shown in Figure 6. Based on Figure 6, the time begins to recover (t2 in Figure 3) and the time to fully recover (t3 in Figure 3) can be obtained as the average time through the simulation, and a shorter t3 indicates that more rupture points can be repaired per unit time and the pipe section can reach the minimum operation level earlier.

Methodology of OPNS Resilience Exposed to Earthquake and the Resilience Partitioning Platform
In the present study, the quantification result of the pipe section was visualized using ArcGIS 10.9 by utilizing natural breaks to separate each index into five assessment grades. The higher the score for each index, the stronger each system's abilities are, and the more resilient the OPNS is. The combination of resilience model and GIS technology adopts the source code programming method. This method takes ArcGIS as an integration platform and uses the Java programming language and SQL statements, which connect resilience as an analysis module with ArcGIS to provide a unified user interface. Combined with the proposed model framework, an OPNS resilience partitioning platform is developed.

Methodology of OPNS Resilience Exposed to Earthquake and the Resilience Partitioning Platform
In the present study, the quantification result of the pipe section was visualized using ArcGIS 10.9 by utilizing natural breaks to separate each index into five assessment grades. The higher the score for each index, the stronger each system's abilities are, and the more resilient the OPNS is. The combination of resilience model and GIS technology adopts the source code programming method. This method takes ArcGIS as an integration platform and uses the Java programming language and SQL statements, which connect resilience as an analysis module with ArcGIS to provide a unified user interface. Combined with the proposed model framework, an OPNS resilience partitioning platform is developed. The spatial analysis and visualization function of the rasterized ArcGIS can partition the subsystem's resistive ability, adaptive ability, recovery ability, and resilience of OPNS accurately, efficiently, and conveniently, and realize the visualization of results partitioning. The framework for the OPNS resilience partition platform is shown in Figure 7. The spatial analysis and visualization function of the rasterized ArcGIS can partition the subsystem's resistive ability, adaptive ability, recovery ability, and resilience of OPNS accurately, efficiently, and conveniently, and realize the visualization of results partitioning. The framework for the OPNS resilience partition platform is shown in Figure 7. Combining the developed partitioning platform and the quantitative resilience assessment model of OPNS, a six-step methodology is further proposed (see Figure 8) to analyze the resilience of OPNS under an earthquake. First, data is collected related to the OPNS, including geographic details, the kind of safety barrier, the pipe diameter, length, and material of each pipe segment. Second, we confirm the earthquake frequency and intensity parameters in the region where each pipe segment is located according to the catastrophes that need to be evaluated. The next step ensures the maintenance plan ac- Combining the developed partitioning platform and the quantitative resilience assessment model of OPNS, a six-step methodology is further proposed (see Figure 8) to analyze Sustainability 2023, 15, 972 9 of 20 the resilience of OPNS under an earthquake. First, data is collected related to the OPNS, including geographic details, the kind of safety barrier, the pipe diameter, length, and material of each pipe segment. Second, we confirm the earthquake frequency and intensity parameters in the region where each pipe segment is located according to the catastrophes that need to be evaluated. The next step ensures the maintenance plan according to the severity of the disaster, and Step 4 calculates three assessment indexes of pipeline resilience combined with the determined parameters, respectively.
Step 5 calculates the resilience of the pipeline segment, and the OPNS resilience is then obtained after weighting. Finally, the last step achieves the partition and visualization of the results based on the OPNS resilience partitioning platform.

Data Collection
This paper uses the buried OPNS in Guangdong Province, China, as a case study. Analysis reveals that the total length of this network is 1000.776 km, the flow direction of the material is from west to east, and the entire flow of the pipe network is 700 m 3 /h. Table  2 and Figure 9 detail each pipe section's parameters and network topology.

Data Collection
This paper uses the buried OPNS in Guangdong Province, China, as a case study. Analysis reveals that the total length of this network is 1000.776 km, the flow direction of the material is from west to east, and the entire flow of the pipe network is 700 m 3 /h. Table 2 and Figure 9 detail each pipe section's parameters and network topology.
According to history, the seismic intensity of this province is basically between VI and VII, and the recurrence period of an earthquake is between 80 years and 100 years, respectively. Based on the Chinese seismic intensity scale (GB/T 17742-2020) [54], the PGA = 1.35 m/s 2 and PGV = 0.12 m/s earthquake was selected as the natural disaster case study based on the most adverse principle. Assuming partial soil liquefaction for all pipes, Table 3 summarizes the several types of safety barriers that can be set up for each pipe section.

Data Collection
This paper uses the buried OPNS in Guangdong Province, China, as a case study. Analysis reveals that the total length of this network is 1000.776 km, the flow direction of the material is from west to east, and the entire flow of the pipe network is 700 m 3 /h. Table  2 and Figure 9 detail each pipe section's parameters and network topology.

Calculation of Resilience Assessment Indexes
Set the number of simulations in the quantification of Rec to 10,000, and prepare for restoration by assuming a suspended time (t 0 -t 1 ) of 0.5 weeks [28]. The event trees required by the calculation of P accident are listed in Figures 10-13, and Figure 14a depicts the resilience curve of pipe subsystems for the reliability level at various times. The results of the resilience indexes of each pipeline are summarized in Table 4, and Table 5 Table 4. Quantification results of resilience indexes in the pipeline network.

No.
P 0 P 1 P 2 P 3 t 0 t 1 t 2 t 3 Re Ad Rec   Figure 15 shows the visualization results of grade division for each index. From Figure 15a, the distribution of the Re in each pipe section in the OPNS tends to decrease from west to east when facing an earthquake of VII intensity, and the quantification of pipes 12, 13, and 14 is relatively low. This is because the length of these pipes is longer than others and spans a more extensive geological range, resulting in a broader range of disturbances affected by the earthquake, and it is even more challenging to maintain the stability of its functions in the event of disasters.
As for the adaptive ability, Figure 15b shows that all the pipes in this network can adapt to the accident evolution under earthquake disasters to a certain extent and reduce accidents after damages caused by subsequent disasters due to the involvement of safety barriers. However, pipes 3, 7, 11, and 14 are on a low level compared with other pipes, which means that these pipes are at greater risk of fire or explosion accidents influenced by the system performance and further damage to the structure caused by the earthquake. Therefore, setting up an emergency team near these sections in the deployment of an emergency rescue force can reduce the system's deterioration risk and prevent further expansion of subsequent accidents.
Regarding recovery ability, it can be seen in Figure 15c that the spatial distribution trend is the same as that of Re (a gradually decreasing trend from west to east). Furthermore, the Rec level for pipes 12, 13, and 14 are also lower than the others according to Figures 14b and 15c, mainly because a longer pipeline will generate more rupture points and requires more recovery time under the same repair efficiency. Therefore, by increasing the number of emergency teams and the daily maintenance frequency of pipes 12, 13, and 14, the repair rate of the pipeline section can be further improved so that the OPNS can quickly reach the lowest performance level to meet the daily needs after the disaster. As for the adaptive ability, Figure 15b shows that all the pipes in this network can adapt to the accident evolution under earthquake disasters to a certain extent and reduce accidents after damages caused by subsequent disasters due to the involvement of safety barriers. However, pipes 3, 7, 11, and 14 are on a low level compared with other pipes, which means that these pipes are at greater risk of fire or explosion accidents influenced by the system performance and further damage to the structure caused by the earthquake. Therefore, setting up an emergency team near these sections in the deployment of an emergency rescue force can reduce the system's deterioration risk and prevent further expansion of subsequent accidents.
Regarding recovery ability, it can be seen in Figure 15c that the spatial distribution trend is the same as that of Re (a gradually decreasing trend from west to east). Furthermore, the Rec level for pipes 12, 13, and 14 are also lower than the others according to Figures 14b and 15c, mainly because a longer pipeline will generate more rupture points and requires more recovery time under the same repair efficiency. Therefore, by increasing the number of emergency teams and the daily maintenance frequency of pipes 12, 13, and 14, the repair rate of the pipeline section can be further improved so that the OPNS can quickly reach the lowest performance level to meet the daily needs after the disaster. Table 5 shows the quantification of resilience under an earthquake of each pipe calculated by Equation (8). Figure 15d deploys the differentiated and expressed results spatially.  Table 5 shows the quantification of resilience under an earthquake of each pipe calculated by Equation (8). Figure 15d deploys the differentiated and expressed results spatially.

Calculation of Network Resilience
From Figure 15d, it can be obtained that resilience still decreases gradually from west to east. For single pipe sections, the resilience of pipes 12, 13, and 14 is lower than others due to the low Re and Rec, which means that they are more likely to be damaged and cause fires, explosions, or other dangerous accidents in the face of the VII earthquake. Furthermore, the length of these pipes is too long, and it takes longer to recover to the minimum performance level, which turned these pipes into objects requiring special attention in OPNS.
When it comes to the OPNS, the quantification of resilience is 0.7753, as in Equation (2). The result means that this network can resist an earthquake of seismic intensity VII, adapt to the evolution of the post-earthquake scene after the quake, and simultaneously avoid all kinds of derivative accidents. Meanwhile, it would take nearly five and a half weeks to recover to the optimal overall performance level to meet the city's daily needs.

Discussion
The method proposed in this paper is used to quantitatively evaluate the resilience of OPNS under different intensities. Considering that the OPNS seismic fortification intensity is VI and above, the content of Section 3.1 shows that the Re of the pipeline has tended to 0 when the earthquake scene is above VIII degrees, indicating that the pipeline has been completely destroyed and that the research value is limited, so an earthquake with an intensity of VI~VIII is selected for the initial disaster scenarios. Assume that the ground motion parameters in the case area are the same for simplification and that all the subsystems can be fully recovered. The resilience evaluation results under different intensities are shown in Table 6.  Figure 16 depicts the resilience curve of pipe subsystems for the reliability level at various times under different disaster scenarios. When PGA exceeds 0.2 g, almost all subsystems cannot achieve a complete performance-recovery within 3 months, indicating that the earthquake of intensity VIII is fatal to the OPNS, which will not only cause severe damage to the structure but also exponentially prolong the repair time of the subsystem. As a result, the OPNS cannot effectively obtain its original state recovery to fulfill the daily demand in this scenario. Furthermore, when the PGA is less than 0.2 g, although the performance of each subsystem will be degraded to a certain extent when earthquake damage comes, the system can be repaired quickly and restored to the initial state due to the high efficiency of the repair scheme, so that OPNS can be maintained at a high resilience level.
As for the three abilities of resilience, considering that the Ad is reflected in the system and prevents the scenario from evolving to an unsafe state through the configuration of the safety barrier, the change in the Ad is mainly related to the configuration and quantity of the barrier. Furthermore, according to Section 3.2 and Figure 3, the involvement of barriers can improve P 2 significantly by mitigating and preventing accident scenarios, making the system more adaptable when facing an earthquake.
Moreover, the Re and Rec decrease with the increase in PGA, as shown in Figure 17a,b, indicating that Re and Rec depend on the severity of earthquake. Interestingly, the Rec decreases rapidly before Re, mainly due to the restoration time vastly increasing when the PGA increases from 0.1 g to 0.15 g. As a result, shortening the restoration time by applying a more efficient strategy can be a great measure of resilience enhancement.
Additionally, with increasing the severity of earthquake, the resilience rapidly decreases and then approaches the minimal resilience value, as shown in Figure 17c,d. The minimal resilience is close to 0.5 for the assumption that all the systems can achieve a complete restoration without considering the repair time, indicating that the OPNS can resist the disturbance of disasters to a certain extent and keep its performance stable in the seismic scenarios with VII and below. When the intensity is above VII, the resilience decreases rapidly due to the possibility, and the degree of structural damage to the system increase sharply, which will prolong the system's recovery time and make it difficult to restore to a specific performance level in time. As for the three abilities of resilience, considering that the Ad is reflected in the system and prevents the scenario from evolving to an unsafe state through the configuration of the safety barrier, the change in the Ad is mainly related to the configuration and quantity of the barrier. Furthermore, according to Section 3.2 and Figure 3, the involvement of barriers can improve P2 significantly by mitigating and preventing accident scenarios, making the system more adaptable when facing an earthquake.
Moreover, the Re and Rec decrease with the increase in PGA, as shown in Figure  17a,b, indicating that Re and Rec depend on the severity of earthquake. Interestingly, the Rec decreases rapidly before Re, mainly due to the restoration time vastly increasing when the PGA increases from 0.1g to 0.15g. As a result, shortening the restoration time by applying a more efficient strategy can be a great measure of resilience enhancement. Furthermore, the inflection point for the rest of the pipes appears at PGV = 0.2 g, whereas the resilience of pipes 12, 13, and 14 starts to decline quickly at PGV = 0.15 g and slowly at PGV = 0.2 g (see Figure 17c), indicating that these pipes are more sensitive to changes in PGV and are more susceptible to losing resilience when confronted with earthquakes. This result is similar to the above assessment, demonstrating that the most vulnerable part in this OPNS is pipes 12, 13, and 14. As such, the OPNS can be made more resilient by enhancing the resilience of these parts.
Finally, comparing the change curves of Re and Rec with PGA (see Figure 17c,d), it can be seen that the variation trend of resilience is compatible with the variation trend of Re, indicating that Re can significantly influence the variation trend of resilience. This result demonstrates that the resistive ability is considered the primary index that contributes to system resilience [39], and reasonably improving the system's Re will be the most efficient strategy to increase resilience. Additionally, with increasing the severity of earthquake, the resilience rapidly decreases and then approaches the minimal resilience value, as shown in Figure 17c,d. The minimal resilience is close to 0.5 for the assumption that all the systems can achieve a complete restoration without considering the repair time, indicating that the OPNS can resist the disturbance of disasters to a certain extent and keep its performance stable in the seismic scenarios with VII and below. When the intensity is above VII, the resilience decreases rapidly due to the possibility, and the degree of structural damage to the system increase sharply, which will prolong the system's recovery time and make it difficult to restore to a specific performance level in time.
Furthermore, the inflection point for the rest of the pipes appears at PGV = 0.2g, whereas the resilience of pipes 12, 13, and 14 starts to decline quickly at PGV = 0.15g and slowly at PGV = 0.2g (see Figure 17c), indicating that these pipes are more sensitive to changes in PGV and are more susceptible to losing resilience when confronted with earthquakes. This result is similar to the above assessment, demonstrating that the most vulnerable part in this OPNS is pipes 12, 13, and 14. As such, the OPNS can be made more resilient by enhancing the resilience of these parts.
Finally, comparing the change curves of Re and Rec with PGA (see Figure 17c,d), it can be seen that the variation trend of resilience is compatible with the variation trend of Re, indicating that Re can significantly influence the variation trend of resilience. This result demonstrates that the resistive ability is considered the primary index that contributes to system resilience [39], and reasonably improving the system's Re will be the most efficient strategy to increase resilience.

Conclusions
The present study aimed to develop a systematic quantitative resilience assessment approach for OPNS exposed to earthquakes. Through the case study application, the following conclusions are obtained: (1) The resistive ability, the adaptive ability, the recovery ability, and the importance coefficient of each subsystem play a significant role in the resilience of OPNS. The damage severity and recovery time may increase with seismic intensity, resulting in the reduced resistive ability and recovery ability of the investigated system. Furthermore, removing safety barriers can significantly exacerbate the decline in reliability, reducing the adaptive ability of the system. Considering that the importance of each subsystem is different, underestimating the importance coefficient and the indexes mentioned above will lead to an inaccurate assessment of OPNS resilience; (2) The case study demonstrates that the recovery ability is the most sensitive ability index among the three indexes in resilience. In contrast, the resistive ability is recognized as the primary index contributing to system resilience, since it can influence the fluctuation trend of resilience; (3) Under different seismic intensities, the OPNS resilience can vary considerably. In an earthquake scenario with an intensity below VII, the pipe network belongs to the resilient OPNS. When the intensity exceeds VII, the resilience decreases rapidly with the increase in PGA, and the OPNS will not belong to the resilience system when the PGA is more significant than 0.2 g; (4) By combining the resilience assessment model with GIS technology, it is possible to visualize resilience partitioning and increase the application efficiency of this methodology. The outcome of this research can be used to support the decision-making of regional planning, emergency management, disaster prevention, and mitigation in OPNS.
In order to improve and evolve the proposed methodology, future research should adopt the assessment parameters of system resilience most coinciding with the facts. For example, the system's functionality should be considered more comprehensively, including the users' satisfaction, the network connectivity, and the efficiency of oil transportation. In addition, calculating the ground motion parameters using the ground motion prediction equation is meaningful for assessing OPNS resilience. Further efforts are also required to quantify the effect of resilience-strengthening strategies, which can optimize the system's resilience and reduce the impact of the disaster.
Author Contributions: J.M. developed the methodology and its application, analyzed the results, drafted the manuscript, and finalized it; G.C. and Y.Z. administrated and supervised this study, acquired funding and data, and provided suggestions on the manuscript; T.Z. provided support for the assessment model and contributed to manuscript writing; L.Z. and J.Z. provided suggestions on the manuscript and polished the manuscript. All authors have read and agreed to the published version of the manuscript.

Data Availability Statement:
The data presented in this study are available on request from the corresponding author and are not publicly available due to confidentiality.