Spatial Analysis of Leakage Accidents of Oil Pipeline River Crossing Sections

All published articles of this journal are available on ScienceDirect.

RESEARCH ARTICLE

Spatial Analysis of Leakage Accidents of Oil Pipeline River Crossing Sections

The Open Civil Engineering Journal 30 Nov 2016 RESEARCH ARTICLE DOI: 10.2174/1874149501610010836

Abstract

Oil pipeline leakages often bring catastrophic impacts on rivers along the pipeline, drinking water sources, and farm land irrigation etc. Based on the distinct spatial-temporal attributes in the analysis of oil spill, diffusion and consequences thereof, we present the theoretical concept and reasoning for leakage accidents spatial analysis through combining GIS and accident aftermath model. Specifically, we aim to achieve quantitative analysis and visualization of oil spill accidents through several GIS functions, which include buffer analysis, overlay analysis, spatial statistical analysis and visualized analysis. Taking a specific leakage accident of oil pipeline river-crossing section as an example, we have analyzed the polluted areas under different diffusion time and the hydrological conditions. In addition, we have evaluated populations and farm land areas affected in terms of drinking water and irrigation as well as their spatial distributions. The results demonstrate that our method can achieve quantitative analysis and visualized representations of oil spills to support policy-makers in making contingency plans.

Keywords: GIS, Oil pipelines, Oil spills, Quantitative analysis, River, Spatial analysis.

1. INTRODUCTION

Environment is the foundation of survival and the development of mankind. Environmental pollution has become one of the most severe threats against the survival of humans. Among different types of pollution, water environment contamination directly endangers life and health of people [1]. The oil pipeline leakage is one of the main causes of water pollution. In recent years, numerous researches have been conducted across the world in water pollution accidents and Oil pipeline leakages. Wu and Alireda et al. [2, 3] carried out risk evaluations of offshore crude oil pipeline based on Bayesian Network model. Joanna evaluated the reliability and risks associated with the port oil pipeline transportation system [4]. Zhang and Zhao et al. [5-7] proposed early-warning models for river pollution contingencies based on 1D, 2D dynamics water quality model, and EFDC model, respectively. By using six-point and four-point implicit difference scheme for numerical solutions, Xu and Wu et al. [8, 9] established engineering-oriented models suitable for emergency management of sudden water pollution incidents. Through using Cartesian cut cell method, Cheng and Wu et al. [10, 11] achieved the computation of pollutant diffusion in complex river courses as well as its GIS visualization. Li and Zhang et al. [12-15] constructed water quality early-warning system after integrating GIS with water pollutant diffusion models. Anjuman and Wang et al. [16-19] undertook quantitative analyses over the impacts of oil pipeline leakages on surface and underground water. Song and Yang et al. [20, 21] statistically analyzed the causes for oil & gas pipeline failure as well as their damage on environment. Despite emerging studies on pollutant diffusion in river oil spills, few existing literature touch upon water source contamination caused by such accidents or statistical analysis and visualization of their further impacts on supply of drinking and irrigation water. To fill this void, here we present a spatial analysis approach for leakage accidents of the river crossing sections of oil pipelines.

2. MODEL FOR ACCIDENT AFTERMATH ANALYSIS

2.1. Computation Model for Leakage Amounts

During oil pipeline leakages, the maximum oil spill is comprised of spill during the responding time and maximum spill from isolated sections after valve closure, as expressed in Equation (1).

(1)

Where, W denotes the maximum oil spill amounts (kg); Q represents the leaking rate within the responding time (kg/s); t denotes the respond time(s); Wrem represents the maximum spill amounts from isolated sections after valve closure (kg).

2.1.1. Leaking Rate within the Responding Time

For small hole leaking, the leaking rate within the responding time can be calculated through Bernoulli [22].

(2)

Where Q denotes the liquid leaking rate (kg/s); p1 denotes the liquid density (kg/m3); δ1 is the liquid leakage coefficient, frequently from 0.6 to 0.64; p represents the pressure in the container (Pa); p 0 represents the environmental pressure (Pa), A 0 denotes the surface area of the crack(m2) by generally taking the diameter of 20% to 100%; g is the gravity acceleration rate (9.8 m/s2); h1 denotes the height of liquid level at the leakage site (m). Liquid leaking rate depends on the pressure difference between the container inside and outside.

For pipeline rupture, leakage rate within responding time equals to the pipeline oil transportation rate, as expressed in Equation (3).

(3)

Where Q denotes the liquid leaking rate (kg/s); Qtra represents the oil transportation rate within the pipeline (kg/s).

2.1.2. Maximum Oil Spill from Isolated Sections

Spill of remaining oil in the pipeline after valve closure depends on the leakage site pressure, which in turn depends on the elevation of leakage site. In Fig. (1), red segments of the curve represent the pipeline sections with potential leakage. The leak detection system based on negative pressure wave could locate the leakage site, as shown in Fig. (1). Thus, the maximum oil spill from isolated sections can be calculated through Equation (4).

Fig. (1). Diagram for calculating maximum leakage from isolable pipeline sections.
(4)

Where R is the pipeline radius (m); p denotes the oil density within the pipeline (kg/m3); Li (i = 1, 2, 3...) represents the pipeline length (m).

2.2. Water Pollution Prediction Model

If only considering leaked oil’s process of plug-flow transfer, dispersion and attenuation at the longitudinal direction, we can use the one-dimensional analysis model to establish a prediction model of transient emission from point-source [23]. The basic format of one-dimensional water quality model can be expressed as Equation (5).

(5)

Considering some factors that emission from oil spill is not stably continuous and concentration at any point is a function of time, we can derive a concentration prediction model as expressed in Equation (6).

(6)

where C is the pollutant concentration at the predicted site (mg/L); x is longitudinal distance between predicted site and leakage site (m); t is the predicted time (s); M denotes the total amount of oil spill (kg); h represents the average water depth of the river (m); σx is the longitudinal dispersion coefficient of the river (m2/s); u is the average longitudinal flow rate of the river (m/s); k denotes the pollutant degradation rate inside the river. As crude oil is hard to degrade, here k = 0.

Among them, river longitudinal dispersion coefficient is obtained through empirical formula based on the flat and straight river stretches. The calculation formula is as Equation (7).

(7)

Where I is the river bed gradient; g is the acceleration rate of gravity, set as 9.8m/s2.

3. PRINCIPLES OF SPATIAL ANALYSIS

Methods for leakage accident spatial analysis is mainly comprised of buffer analysis, overlay analysis, spatial statistical analysis and visualized analysis.

3.1. Buffer Analysis

Buffer analysis is the distance analysis based on the topological relationships among spatial targets (point, line, and plane). Its core idea is to specify spatial targets and define their surrounding areas, the sizes of which depend on the areas’ radii [24]. By using water quality prediction model to calculate the pollutant concentration in the water, we can derive the river pollution buffer analysis model as in Equation (8).

(8)

Where A is the polluted river area; x 0, y 0 represent the point coordinate of leakage accident (m); r denotes the maximum diffusion distance of oil film (m); R+represents downstream river area of the leakage accident site.

3.2. Overlay Analysis

Overlay analysis seeks to generate new features or establish spatial correlations between geographical objects by overlaying two or more geographical feature layers. Once we have the polluted river areas, we can overlay it with the layer of water withdrawal sites to identify the contaminated ones. Its mathematical model is as Equation (9).

(9)

Where C represents the assemblage of environmentally sensitive sites that are contaminated; A denotes the region that oil film flows through; B represents the assemblage of environmentally sensitive sites.

3.3. Spatial Statistical Analysis

Spatial Statistical Analysis is also known as Spatial Data Statistical Analysis. At its essence, the method seeks to establish the statistical relationship among data based on spatial locations as well as identify spatial dependence. Correlation, self-correlation within data related to geographical positions. Through overlay analysis, we obtain the polluted water withdrawal sites and the areas affected by accidents, then combine water withdrawal sites with topological relationship between population and farm land through spatial statistical analysis. Its mathematical model can be expressed as Equation (10).

(10)

Where S denotes areas affected by the accidents; Myi represents the population whose drinking water sources are polluted (person); Mgi represents the area of farm land whose irrigation water sources are polluted (Km2), i = 1, 2, 3...

3.4. Visualization Analysis

Visualization refers to the conversion of symbol or data into intuitive geometric figures that are convenient for researchers to monitor the simulation and computation processes [25]. This method primarily consists of visualized representations of map data, geographic information, spatial analysis results as well as interactive demonstration of graphs and data of analysis results.

4. EXEMPLARY CALCULATION

As shown in Fig. (2), the river crossing oil pipeline in our study has a diameter of D610mm and an operation pressure of 10MPa. The nearest valve chambers at upstream and downstream of river crossing section are RTU(Remote Terminal Unit)-based. Moreover, a leak detection system based on negative pressure wave is installed in the pipeline. For the purpose of spatial analysis of leakage accidents, we consider a hypothetical situation that the pipeline is ruptured due to illegal operation of a sand dredger, and regard it as small hole leaking. The river hydrological parameters are shown in Table 1.

Fig. (2). Environment overview of the analyzed area.
Table 1.
River Hydrological Parameters.
Season River Flow (m3/s) Flow Rate(m/s) Depth(m) Surface width(m) River bed gradient(‰)
High-flow 6861 3.77 10.1 182 6.3
Medium-flow 1989 2.12 6.9 136 6.3
Low-flow 43 0.57 1.2 63 6.3

4.1. Oil Film Diffusion Distance

4.1.1. Calculation of Leakage

The responding time between occurrence of the accident and closure of valve chambers at two ends is not more than 3 minutes. According to the leakage calculation Equation (1)-(4), p1 = 0.9kg/m2, δ1 = 0.6, p = 1 × 107pa, p 0 = 1.01 × 105pa, A 0 = 0.61 m2, g = 9.8m/s2, h1 = 3m, we can derive liquid leaking rate Q as 190kg/s, the total leakage within the 3 minutes responding time as 57080 kg. Meanwhile, based on the pipeline elevation parameters, we can determine the maximum leakage of isolated sections as 38860 kg. Therefore, the maximum oil spill of this accident is 95940 kg.

Fig. (3). River stretches with excessive pollutants after 30/60/90 min in medium-flow season.

4.1.2. River Stretches with Excessive Pollutants at Different Times

Here we assume that the pipeline’s owner launches a rapid rescue response and dispatches repair workers to the accident site within 90 min when leakage occurs. Combining water quality prediction model with buffer analysis based on GIS platform, we can determine the river stretches with excessive pollutants in medium-flow season with oil film diffusion times being 30 min, 60 min and 90 min respectively. The visualized results are shown in Fig. (3).

4.1.3. River Stretches with Excessive Pollutants Under Different Hydrological Conditions

Combining water quality prediction model with buffer analysis based on GIS, we can determine the river stretches with excessive pollutants in low-flow, medium-flow or high-flow seasons with the oil film diffusion time being 90 min. The visualized results are shown in Fig. (4).

Fig. (4). River stretches with excessive pollutants after 90 min in low-flow, medium-flow or high-flow seasons.

4.2. Areas with Contaminated Water Sources along River Stretches that Oil Film Flows through

Based on excessively polluted areas, we can use buffer analysis to determine riverine regions that oil film flows through (namely contaminated regions) at different time points. In the next step, we undertake overlay analysis by superposition of affected riverine region layer and water withdrawal sites layer. On the basis of the supply-demand topological relationships between water withdrawal sites and residential areas plus farm land, we can compute the areas mostly affected by the accident. The results are visualized, as shown in Figs. (5-7). Through spatial statistical analysis, we also obtain the data of areas with affected water sources along river stretches that oil film flows through, as demonstrated in Table 2, Figs. (8-10).

Fig. (5). Affected areas after 30/60/90 min in low-flow season.
Fig. (6). Affected areas after 30/60/90 min in medium-flow season.
Fig. (7). Affected areas after 30/60/90 min in high-flow season.
Fig. (8). Diffusion distance at different time under the same hydrological condition.
Fig. (9). Diffusion distance at 90 min under different hydrological condition.
Table 2.
Statistics of areas affected by water source pollution along river stretches that oil film flows through.
Hydrology Condition Affected objects Influence 30min 60min 90min
Low-flow season Resident Drinking water 0 0 0
Farmland Irrigation 0 0 8.4 Km2
Medium-flow season Resident Drinking water 0 0 5370 Person
Farmland Irrigation 8.4 Km2 8.4 Km2 8.4 Km2
High-flow season Resident Drinking water 0 5370 Person 9390 Person
Farmland Irrigation 8.4 Km2 8.4 Km2 25.1 Km2
Fig. (10). Area of regions with excessive pollutants under different conditions.

4.3. Analysis and Discussion

Collectively, Fig. (8) demonstrates areas with excessive water pollutant and related statistical data at different times under similar hydrological conditions. Between 30min and 90min, the oil film central diffusion distance increases from 3.82Km to 11.43Km; the oil film central concentration declines from 4095g/L to 1362g/L; Fig. (10) shows the area of regions with excessive river pollutants increases from 0.208 Km2 to 0.320 Km2. Thus, in a pipeline leakage accident oil film diffusion distance increases and maximum water pollutant concentration decreases with the passage of time. Fig. (9) shows areas with excessive water pollutant and related statistical data at the same time point but under different hydrological conditions. With the hydrological condition changing from low-flow season, to medium-flow and high-flow seasons, the oil film central diffusion distance increases from 3.07Km to 20.35Km; the oil film central concentration declines from 5487g/L to 655g/L. The size of regions with excessive river pollutants increases from 0.327 Km2 to 0.393 Km2 in Fig. (10). The results clearly show that low-flow season is characterized by the lowest oil film diffusion distance and a high maximum water pollutant concentration, indicating the difficulty of pollutant diffusion at this season. In contrast, pollutant diffusion is the easiest in high-flow season due to the highest pollutant diffusion distance and low maximum water pollutant concentration.

According to the calculation results and statistical data of areas with contaminated water sources shown in Table 2, the river stretch flowed through by oil film is the shortest in the low-flow season, in which affected area is kept at the minimum. Within 90min, none of the drinking water sources is polluted and only irrigation water for 8.4 Km2 vegetable plantation field is contaminated. The impacted area is the largest in the high-flow season, when oil film flows through the longest river stretch. Within 90min, drinking water sources for 9390 people are contaminated; irrigation water sources for 25.1 Km2 vegetable/fruit plantation fields are polluted. In the meantime, from variation in affected area during low-flow season shown in Table 2, we can infer that the size of affected area mainly depends on the number of water withdrawal sites near the oil film pollution. If no drinking/irrigation water withdrawal site exists in the stretch flowed through by oil film, the impact of pollution will mainly fall onto the river instead of resident areas and farm land.

Taken together, in a pipeline leakage accident, oil-film-polluted area largely depends on hydrological conditions and responding time of emergency repair by the pipeline owner. Moreover, in regard to areas with affected water sources along the river stretch flowed through by oil-film, their sizes mainly depend on the number of environment sensitive sites such as water withdrawal sites as well as population and farm land associated with affected water withdrawal sites.

CONCLUSION

Our spatial analysis method takes advantage of GIS to reveal the evolving process and spatial features of leakage accidents of pipeline river crossing sections. So our method can be used to enhance the scientific objectiveness and intuitiveness of pipeline/environment safety management, it could have invaluable applications in areas such as environment/ecosystems protection and water supply security.

  1. We have established a GIS spatial analysis method for leakage accidents of pipeline river crossing sections. The method is based on combination of the analysis model of accident aftermath with several core GIS functions, including buffer analysis, overlay analysis, spatial statistical analysis and visualized analysis. With this method, we have achieved orderly quantitative analysis and visualized representations of the complex environment, process of the beginning and development of river oil spills.
  2. We have elaborated the principles and procedures of the spatial analysis for leakage accidents of pipeline river crossing sections. We have also demonstrated the feasibility of the method by conducting calculation and analysis on an example. The results show that we can use this method to determine the maximum distance of oil film diffusion under different diffusion time and hydrological conditions as well as affected water withdrawal sites. By analyzing the topological relationships between water withdrawal sites and residential areas/farm land, we can ascertain the population number and farm land area impacted by the leakage accident as well as their spatial distributions.
  3. Compared to ocean oil spills, river oil spills are characterized by complex spatial relationship due to more complex hydrological conditions. As a result, the size of polluted region is under the influences of a diverse range of factors. Further studies can be undertaken from these perspectives to enhance the precision of accident aftermath analysis.

CONFLICT OF INTREST

The authors confirm that this article content has no conflict of interest.

ACKNOWLEDGEMENTS

This paper is sponsored by the funding project of Geographic National Condition Monitoring Engineering Research Center of Sichuan Province (GC201517), Planned Project of Science and Technology Support of Sichuan Province (2015SZ0046) and Open Fund Project of Key Laboratory of Structural Engineering of Institution of Higher Education in Sichuan Province-Methods Study of Earthquake Safety Evaluation of Oil and Gas Pipeline on Basis of Environmental Risks. Sichuan Bureau of Surveying, Mapping and Geoinformation Support Issues (J2014ZC12).

REFERENCES

1 X. Han, and T. Huang, "Statistical analysis of sudden water pollution accidents", Water Resour. Prot., vol. 1, pp. 84-86, 2010. 2 W.S. Wu, C.F. Yang, J.C. Chang, and P. Chateau, "Risk assessment by integrating interpretive structural modeling and Bayesian network, case of offshore pipeline project", Reliab. Eng. Syst. Saf., vol. 142, pp. 515-524, 2015.
[CrossRef Link] 3 A. Alireda, K. Faisal, A. Ayhan, and M. Haddara, "Risk assessment of offshore crude oil pipeline failure", J. Loss Prev. Process Ind., vol. 37, pp. 101-109, 2015.
[CrossRef Link] 4 S. Joanna, "Reliability and risk evaluation of a port oil pipeline transportation system in variable operation conditions", Int. J. Press. Vessels Piping, vol. 87, pp. 81-87, 2010.
[CrossRef Link] 5 B. Zhang, Y. Qin, Q. Hao, and X. Wang, "The research of temporal-spatial simulation and dynamic regulation system of water quality on sudden water pollution accidents", Geomatics Info. Sci. Wuhan Univ., vol. 9, pp. 1108-1112, 2013. 6 Y. Zhao, Y. Wang, W. Zhang, and M. Ju, "Two-dimensional numerical simulation model for oil spill in river", Yangtze River, vol. 15, pp. 81-84, 2012. 7 Z. Yu, "J, Xia, and H. Ren, “Water pollution accident emergency response and early-warning model in gansu-ningxia-inner Mongolia section of the yellow river", Yellow River, vol. 4, pp. 37-40, 2014. 8 J. Xu, X. Xie, Y. Fu, and Y. Li, "Spatio-temporal simulation and visualization of sudden water pollution accident", Zhongguo Nongcun Shuili Shuidian, vol. 8, pp. 52-55, 2012. 9 D. Wu, Q. Huang, H. Sun, and Y. Song, "A water pollution diffusion model on emergency response and its visualization in a GIS platform", Geomatics Inf. Sci. Wuhan Univ., vol. 2, pp. 131-134, 2009. 10 J. Cheng, K. Wang, X. Feng, and Y. Wang, "GIS-based water quality simulation and visualization applied in a watershed", J. Hydraul. Eng., vol. 11, pp. 1352-1360, 2014. 11 P. Wu, "Geographic information system-based environmental impact assessment visualization of complex river course based on Cartesian cut cell method", J. Comput. Appl., vol. 3, pp. 780-784, 2014. 12 G. Li, T. Zhou, B. Zeng, J. Shen, and Z. Li, "Design and development of risk assessment system for accidental water environment pollution in three gorges reservoir region", J. Northwest A&F Univ., vol. 8, pp. 221-228, 2014.(Nat. Sci. Ed.) 13 Y. Zhang, W. Luo, A. Lei, and X. Jiang, "Research and application of river oil spill model in advanced warning system of water quality of three gorges", Environ. Sci. Technol., vol. 1, pp. 128-131, 2010. 14 Z. Xiao, C. Mo, and X. Lei, "Research and application of water pollution early warning system based on GIS", Water Resour. Power, vol. 5, pp. 139-141, 2011. 15 B. Chen, D. Tan, and L. Song, "Application research of GIS technology in suddenly-happened water pollution accident emergency response system", J. Yangtze River Sci. Res. Inst., vol. 1, pp. 29-32, 2010. 16 S. Anjuman, S. Rehan, and T. Solomon, "Risk analysis for oil & gas pipelines: A sustainability assessment approach using fuzzy based bow-tie analysis", J. Loss Prev. Process Ind., vol. 25, pp. 505-523, 2012.
[CrossRef Link] 17 L. Wang, X. Yu, and L. Hu, "Environmental risk assessment on crude oil pipeline leakage to surface water", Environ. Sci. and Manag., vol. 12, pp. 193-194, 2012. 18 H. Ma, Y. Li, and Z. Wu, "Environmental risk evaluation of leaked products pipeline", Oil Gas Storage Transp., vol. 11, pp. 801-804, 2011. 19 Q. Li, G. Liu, and J. Zhang, "The groundwater environmental risks assessment based on fuzzy fault tree analysis method to the oil pipeline", Environ. Sci. Technol., vol. 6L, pp. 328-333, 2013. 20 X. Song, G. Liu, and W. Wang, "The influence of oil pipeline leakage on environmental pollution", Oil Gas Storage Transp., vol. 4, pp. 297-298, 2011. 21 S. Yang, H. Jin, S. Yu, Y. Chen, J. Hao, and Z. Zhai, "Environmental hazards and contingency plans along the proposed China–Russia Oil Pipeline route, Northeastern China", Cold Reg. Sci. Technol., vol. 64, pp. 271-278, 2010.
[CrossRef Link] 22 E. Hu, "Practical technology, methods, and case of environmental risk assessment", China Environ. Sci., vol. 1, pp. 272-287, 2009. 23 Z. Peng, T. Yang, X. Liang, and Z. Gu, "Mathematical model of water environment and its application", Chemical Industry., vol. 1, pp. 19-52, 2007. 24 K. Qin, Theories and Methods of Spatial Analysis in GIS., Wuhan University Press: Wuhan, 2010. 25 A. Tang, and X. Yang, "Experimental course of spatial analysis in GIS", Science, vol. 2, pp. 151-192, 2012.