## Services on Demand

## Article

## Indicators

## Related links

- Cited by Google
- Similars in Google

## Share

## Journal of the South African Institution of Civil Engineering

*Print version* ISSN 1021-2019

### J. S. Afr. Inst. Civ. Eng. vol.53 no.1 Midrand Apr. 2011

**TECHNICAL PAPER**

**The failure probability of welded steel pipelines in dolomitic areas **

**D T Nel; J Haarhoff **

**ABSTRACT**

This paper considers aspects related to the nature of dolomite, sinkholes, the risk classification of dolomitic land, as well as factors affecting the failure of pipelines in dolomitic areas. The information and data presented are used to derive equations that may be used to predict the probability of failure of steel pipelines in dolomitic areas subject to sinkhole formation. Consideration is also given to other factors that may influence the failure of pipelines in dolomitic areas.

**Key words:** dolomitic, sinkholes, pipeline, failure, probability

**INTRODUCTION**

Sections of water distribution pipelines are often installed in dolomitic areas subject to sinkhole formation. The possibility of losing a bulk water pipeline through sinkhole formation, however remote, has an obvious and large impact on the reliability of water supply. This paper develops a methodology to estimate the failure probability of steel pipelines associated with sinkhole formation.

The assessment of the probability of sinkhole formation, and the related sinkhole diameter, draws on the substantial body of experience accumulated in South Africa in the past years. The data and guidelines had to be adapted for the purposes of this analysis. This had to be complemented with a structural and spatial analysis of a steel pipe intersecting a sinkhole. The workability of the approach is demonstrated with a hypothetical example. The paper closes with a critical assessment of the methodology and how it could be improved and refined.

**THE NATURE OF DOLOMITE **

Dolomitic land refers to areas underlain directly or at shallow depth by dolomite, which is a sedimentary rock, has a calcareous composition, is chemically formed and consists of the double carbonate of calcium and magnesium (Lurie 1977).

Chemical composition of dolomite:

Dolomitic rock has a number of unique features (Brink 1979):

It is a compact, impervious rock with a porosity of approximately 0,3%.

It has a highly developed network of joints, tension-fractures and faults that allows water to percolate easily through the rock mass.

It is easily dissolved by carbon dioxide in the groundwater. The carbon dioxide in conjunction with groundwater forms a weak carbonic acid, which dissolves the dolomite to bicarbonates.

The main insoluble residues left from the weathering of dolomite are chert (SiCy, iron from the layers of ferruginous dolomite, and wad.

In the presence of a small concentration of manganese, the dolomite is found in the form of wad. Wad (or manganiferous earth) is an insoluble and highly compressible material that consists of manganese and iron oxides with minor impurities. A lowering of the water table may produce significant ground movement at the surface due to compression of the wad. This ground movement manifests itself as a doline, being a shallow enclosed depression in the ground surface. A doline can be from a few metres across to more than 1 km in length.

**SINKHOLES **

**Types of sinkholes **

Sinkholes (also referred to as sinks) are a type of land subsidence since they involve a vertical downward movement of the land surface. Three different types of sinkholes are distinguished (Tharp 1999; Waltham* et al* 2005):

Solution sinkholes (Figure 1) form where dolomitic bedrock is exposed at the land surface and subjected to weathering by dissolution. Surface water collects in natural depressions and the bedrock slowly dissolves to form a sinkhole.

Cover subsidence sinkholes (Figure 2) are similar to solution sinkholes except that the soluble bedrock is covered by a thin layer of soil or sediment. Surface water infiltration dissolves the dolomite to form cavities where the bedrock is most intensely fractured, and the overlying sediment gradually moves downward into the expanding cavity.

Cover collapse sinkholes (Figure 3) form when surface materials suddenly drop into a subsurface cavity. A cavity will form slowly over time as groundwater moves along fractures in soluble bedrock and enlarges the cavity through dissolution. The actual collapse can occur in two different ways. Firstly, when a cavity gets sufficiently large, the roof of the cavity becomes too thin to support the weight of the overlying rock, sediment or imposed loads, so it collapses into the cavity. Secondly, caves filled with groundwater are sometimes able to support the weight of overlying sediment, but if groundwater levels are lowered, the overlying sediment will first erode and then collapse into the dewatered cavity. The final breakthrough of a cover collapse sinkhole can occur suddenly and it may have catastrophic consequences.

**Sinkhole formation **

Sinkholes result from the hollowing out or formation of a void below the earth's surface. Sinkholes can form naturally as a consequence of normal geological processes, or they may have anthropogenic causes. According to Schoning (1990) anthropogenic sinkhole formation requires three conditions - the right geotechnical conditions, inappropriate development relative to the geotechnical conditions and adequate rainfall. Anthropogenic causes such as the construction of roads, forming of drainage ditches, township development and the associated services, groundwater extraction, groundwater recharge, etc, may cause sinkholes to form.

Buttrick (1992) discussed various factors that may influence sinkhole formation and the size of the sinkhole that will form, namely:

The receptacle within the bedrock or overburden that can receive the mobilised sinkhole material.

The throat size of the conduit that feeds the mobilised soil material into the receptacle.

The blanketing layer of material that overlies the dolomite bedrock.

A mobilising agency that will induce mobilisation of the material in the blanketing layer through the throat and into the receptacle.

The soil angle of internal friction.

The presence and influence of intrusive layers within the blanketing layer that will reduce the mobilisation potential of materials in the blanketing layers.

The position of the water table.

To gain a better understanding of the mechanisms of sinkhole formation, Goodings* et al* (2002) constructed and tested forty-nine physical models of sinkhole development in a laboratory. The models simulated weakly cemented sand overlying cavities in karst limestone, but did not consider the influence of groundwater. The models were divided into two groups - those with nothing above the weakly cemented sand layer, and those with an additional layer of loose, uncemented soil above the cemented sand layer. For the experiments with no uncemented soil overburden, the main findings were:

The parameters critical to predicting failure were the unit weight of the cemented sand (g

_{c}), the thickness of the cemented sand overlying the karst cavity (H_{c}), the true cohesion of the cemented sand (c') and the diameter of the underlying karst cavity (D).Failure was manifested as an intact plug of soil falling into the cavity below. Three modes of failure were observed:

In models with thin layers (H

_{c}/D < 0,25) failure occurred as a breakthrough plug that left a hole in the top of the weakly cemented layer but with little overhang on the inclines.In models with H

_{c}/D = 0,25, failure occurred as a breakthrough plug, but there was some overhang observed on the inclines.In models with thick layers (H

_{c}/D>0,31) failure occurred when a plug of soil dropped into the underlying cavity, but the hole did not break through to the top of the weakly cemented soil layer.

A sinkhole stability chart was developed that can be used to predict full-scale failures under certain conditions, as shown in Figure 4. The x-axis value relates to the ratio of the depth of overburden and the diameter of the underlying cavity, whereas the y-axis value deals with the characteristics of the soil above the cavity where a sinkhole may possibly form. N is a scale factor with N = 1 for full scale models. The stability chart is for failure conditions and includes no factors of safety. If the plotted value of a specific scenario being investigated falls on or above the failure envelope, failure is predicted.

For the experiments with an uncemented soil overburden above the weakly cemented sand, the H_{c}/D ratio was kept at or below 0,25. The main finding was that the additional pressure of the uncemented soil bearing down on the breakthrough plug never exceeded the weight of a cone of sand of diameter D' and height 1,25D', where D' is the diameter of the top of the inclined breakthrough plug.

**Sinkhole geometry**

Three distinct forms related to sinkholes have been recognised, namely shaft, undercut and bowl forms (Hyatt* et al* 1996):

Shaft forms are steep-sided, with a flat bottom.

Undercut forms have a distinct overhang on at least three of the four sides.

Bowl forms are simple depressions with side slopes less than or equal to 90 degrees.

Hyatt* et al* (1996) developed a threedimensional volume ratio index* (VRI)* which compares the actual volume of the sinkhole

(*V _{S}*) to the volume of an elliptical cone

*(V*having the same major radii, minor radii and depth as the sinkhole that has formed, as defined by Equation 3.

_{ec})where

VRI | is the volume ratio index Vs is the volume of the sinkhole that formed |

V _{ec} | is the volume of an elliptical cone with the same major radii, minor radii and depth as the sinkhole that formed |

a | is the major radii of the sinkhole that has formed (length / 2) |

b | is the minor radii of the sinkhole that has formed (width / 2) |

d | is the depth of the sinkhole that has formed |

The* VRI* provides an indication of the geometric form to which the sinkhole is volumetrically most similar, as shown in Figure 5.

Figure 6 reflects the results of a study related to sinkholes that formed in the Dougherty Plain at Albany, Georgia in the USA (Hyatt* et al* 1996). For every sinkhole, the sinkhole* VRI* value was plotted against the diameter of the sinkhole. The results indicated that the generic form of small sinkholes (less than 2 m in diameter) varied greatly. As the sinkhole diameter increased, the* VRI* values converged to an ideal bowl form.

The parameter used to describe the surface shape of sinkholes is the ratio of the minor to the major surface dimension of the sinkhole. Sinkholes tend to show a circular or sub-elliptical shape, and the following ratios have been reported:

Bruno

et al(2008); predominantly in the range 0,62 to 1,0.Hyatt

et al(1996); predominantly from 0,9 to 1,0, but sinkholes with ratios from 0,48 to 0,9 were also encountered.Brinkmann

et al(2008); within a range of 0,153 to 0,954, with an average ratio of 0,715.

Pilecki* et al* (2006) presented the relationship given by Equation 4, which suggests that the greater the thickness of the layer of loose overburden, the greater the area of the sinkhole at the ground surface. Above a critical overburden thickness, however, the relationship breaks down and the sinkhole diameter remains constant, due to the socalled chimney action process.

where

D_{sh} | is the diameter of the sinkhole on the ground surface in m |

Z | is the thickness of the loose overburden in m |

θ | is the angle of internal friction of the loose overburden |

r_{s} | is the radius of the shaft that can receive the sinkhole material |

Aspects pertaining to sinkholes are widely reported in the literature, but data related to the sizes of sinkholes that occurred is limited. Table 1 provides a summary of the dimensions of sinkholes as reflected in a number of sources in the literature. Hyatt* et al* (1996), Bruno* et al* (2008) and Brinkmann *et al* (2008) indicated that the sizes of sinkholes follow a positively skewed lognormal distribution.

In order to establish a statistical distribution of the average sinkhole sizes that could be expected, the data pertaining to the average sinkhole diameter or length and the maximum sinkhole sizes given in Table 1 were used. It effectively translates to a sample size of 1 393 data points, comprising average and maximum sinkhole sizes only. The shaded data was excluded due to insufficient data, or since the large sinkholes sizes reported may be dolines. The software programme, EasyFit Professional Version 5.1, was used to select the numerical function that provides a good fit in respect of the average sinkhole diameter data. The Kolmogorov-Smirnov test was applied to test the goodness of fit of the data compared to a range of hypothesised distribution functions. The Kolmogorov-Smirnov statistic is based on the largest vertical difference between the theoretical and empirical cumulative distribution function. The parameters of the sinkhole size probability density function and the value of the Kolmogorov-Smirnov statistic are given in Table 2. A plot of the Table 1 data used and the fitted cumulative distribution function related to the anticipated sinkhole sizes are shown in Figure 7. The variables in the lognormal distributions are to the base e.

Although the data used comprised a limited data set, it does provide an indication of the probability that a certain sinkhole size may be exceeded. However, this distribution cannot be generally applied to all areas, due to the fact that each area is unique in relation to the factors and conditions that may determine the actual size of the sinkhole that may form. What will add great value is that data related to sinkhole sizes should be gathered in order to provide more reliable sinkhole frequency distributions for specific areas.

**RISK CLASSIFICATION OF DOLOMITIC LAND **

The probability of sinkhole formation exihibits great spatial variations, influenced by the stratigraphy and geological history of an area, as well as by anthropogenic processes. The risk classification of dolomitic areas entails mapping the areas where sinkholes are likely to occur, and should ideally also estimate their size and probability of occurrence.

Galve* et al* (2009) described different methods that can be applied to assess the susceptibility of dolomitic land subject to sinkhole formation, namely:

Deterministic models (based on numerical models or stability analysis) are used to assess the degree of stability at different points. These models are based on numerical supposition, rely on data that is often difficult and expensive to obtain, does not take account of the complexity of subsidence processes and cannot be applied to large areas.

Direct mapping of susceptibility zonations comprise maps produced by experts with a good knowledge of the study area and the phenomena, based on specific criteria used.

Susceptibility maps based on the spatial distribution of sinkholes. A number of approaches exist:

Heuristic methods that base susceptibility assessments on the establishment of a scoring system to a group of conditioning factors.Sinkhole density maps.

Sinkhole susceptibility zonations determined by the distance of each point in an area to the nearest sinkhole and expressing it in terms of a nearest neighbour index.

Probabilistic methodologies derive the probability of sinkhole formation from the analysis of statistical relationships between the known sinkholes and a group of factors used to predict sinkhole formation.

Thomas* et al* (1999) described a number of site characterisation methods which may be used to locate sinkholes or subsurface voids. It is stated that the reliability of the methods to locate existing sinkholes is generally good, but that the reliability associated with finding subsurface voids is dependent on the number of probes utilised in the area under investigation. These methods can comprise deterministic models, direct mapping of zones subject to sinkhole formation or the preparation of susceptibility maps.

A number of heuristic methods have been developed. The method of scenario supposition for the stability evaluation of dolomitic areas was developed to characterise dolomitic areas in respect of the risk of sinkhole formation (Buttrick 1992; Buttrick* et al *2001; Department of Public Works 2004). In terms of this method, dolomitic areas are divided into eight inherent risk classes. Each inherent risk class is a function of the size of the sinkhole that may develop, as well as the anticipated number of ground movement events that may occur per unit area per year, based on the assumption that infrastructure services in dolomitic areas are poorly maintained. The framework of this method is outlined in Table 3.

Buttrick (2001) related the number of anticipated ground movement events that may occur to an area's inherent risk categorisation, as reflected in Table 4.

The above review on sinkhole probability and sizes provides a number of useful preliminary pointers. The size distribution of sinkhole sizes shown in Table 1 with its mathematical description in Figure 7 should be further reinforced with an analysis of actual sinkhole data. The risk classification framework of Buttrick (2001) provides an estimate of the probability of sinkhole formation, as well as a rough indication of which sizes can be expected. Before the analysis of sinkhole sizes causing pipeline failures can be done, the factors affecting the failure of pipelines in dolomitic areas have to be considered. Following this, one can narrow down the problem by establishing the critical sinkhole size which would pose a real danger to a welded steel pipe collapsing. This is the topic in the remainder of the paper.

**FACTORS AFFECTING THE FAILURE OF PIPELINES IN DOLOMITIC AREAS **

**The soil angle of internal friction **

The soil angle of internal friction plays an important role in determining the size and shape of the sinkhole that may form, influenced by the thickness of the blanketing layer, intrusive layers within the blanketing layer and the position of the water table.

Table 5 summarises the angle of internal friction of various soil types (Spangler* et al* 1984; Buttrick 1992; Pilecki* et al* 2006). The angle of internal friction of a specific soil may be further influenced by the degree of soil compaction, the soil moisture content, etc.

**The pipeline depth of cover **

Assuming that the geometric form of sinkholes is predominantly conical in shape (VRI = 1), it follows that the unsupported length of pipe spanning across a sinkhole will increase as the pipeline depth of cover decreases. The depth of cover over a pipeline is governed by a range of factors, such as local soil conditions, pipeline loading, water table, third party activities in the area, etc. Table 6 reflects data used by Rand Water related to the minimum depth of cover when pipelines are installed (Turnbull 1996).

An envelope line intersecting or being below the data in Table 6 is given by Equation (5).

where

h_{cover} | is the pipeline depth of cover in m |

D | is the diameter of the pipeline in m |

**The pipeline wall thickness **

Bulk water distribution pipelines are designed to safely accommodate internal pressure, external loads, handling, buckling, circumferential and longitudinal deflections, and to span between supports where required. Depending on the pipeline coating or lining that will be applied, special consideration has to be given to limit the strain in the pipeline material.

Analysis of Rand Water's steel pipeline data (Reyneke 2007) indicated that pipelines are generally conservatively designed. The pipeline design stress is limited to 55% of the yield stress, while steel pipelines that have a cement mortar lining are designed so that the strain in the pipeline does not exceed a value of 800x10^{-6}. Figure 8 provides data pertaining to the D/t ratio versus pipe diameter relationship of Rand Water steel pipelines installed from 1907 to date.

The conservative envelope line around the data of Figure 8 is described by Equation (6).

where

D | is the diameter of the pipeline in m |

t | is the pipe wall thickness in m |

**The pipeline safe span distance **

AWWA (2003) provided a table that may be used to determine the safe span distances for a simply supported steel pipe, completely filled with water and supported at intervals on 120 degree contact saddles. Table 7 illustrates that a relationship exists between the pipeline safe span distance, the pipe diameter* (D)* and the pipe wall thickness* (t). *The* D/t* ratio has been determined for each row and column intersection, and this value is reflected together with their corresponding pipeline safe span distance* (L _{s-spar} )* in Table 7.

Exponential interpolation between the values of Table 7, using the D/t ratio in Equation (6), allows a direct plot of the maximum safe span distance against the pipe diameter, as shown in Figure 9. A trendline was fitted to these data points and Equation (7) defines this relationship mathematically.

where

L_{s-span} | is the pipeline safe span distance |

D | is the diameter of the pipeline in m |

**SINKHOLE DIAMETERS, POSITIONS AND THE NUMBER OF EVENTS INFLUENCING THE FAILURE OF STEEL PIPELINES **

**Critical sinkhole diameter for a pipe running through its centre **

A sinkhole that may form under specific conditions and that is large enough in diameter will leave a certain length of unsupported pipe across the sinkhole (*L _{u-span}*), as shown in Figure 10. The following assumptions are made related to the derivations that follow:

A sinkhole is conical in shape, and symmetrical about its vertical axis.

Only steel pipelines with welded joints are considered.

The throat size of a sinkhole, which can vary greatly, is not taken into account in terms of its potential to affect the failure of the pipeline passing through the sinkhole.

Equation (8) is derived by considering the geometrical relationships within Figure 10 and expresses the diameter of the sinkhole as a function of the pipe diameter, pipeline unsupported span distance across the sinkhole, soil angle of internal friction, as well as the pipeline depth of cover.

where

D_{sh} | is the diameter of the sinkhole at the ground surface in m |

L _{u_span} | is the pipeline unsupported span distance across the cavity caused by the sinkhole in m |

h_{cover} | is the pipeline depth of cover in m |

D | is the diameter of the pipeline in m |

θ | is the soil angle of internal friction in degrees |

Equation (9) is obtained by substituting *h _{cover}* from Equation (5) into Equation (8), as well as setting

*L*and substituting it into Equation (8). Equation (9) provides the maximum sinkhole diameter that can be spanned by a pipeline across its centrewithout structural collapse.

_{u_span}= L_{s_span}Equation (9) has been applied to a range of pipeline diameters and soils with different angles of internal friction. The results thereof are shown in Figure 11, illustrating that sinkholes with a surface diameter less than 15 m pose no risk to pipelines larger than 600 mm in diameter where the soil angle of internal friction is within the range 15º __<__ θ __<__ 45º.

**The zone of influence of sinkholes larger than the critical sinkhole diameter **

In the previous section, it was assumed that the pipeline runs through the centre of the sinkhole. In this section, the analysis will be broadened to cases where the pipeline may intersect a sinkhole at any position. The objective is to estimate the width of the pipeline failure strip within the sinkhole, associated with sinkholes larger than the critical sinkhole diameter that will cause pipeline failure. Figure 12 depicts a sinkhole in different positions relative to the pipeline.

For the purpose of the derivation that follows, it is assumed that the sinkhole approaches and crosses the pipeline at right angles. The first critical point occurs as the sinkhole is moving towards the pipeline (from left to right) and is just in the position where the pipe reaches the limit of its unsupported span distance,* L _{s_span}*. As the

sinkhole moves further towards the right, through a distance of* X _{fs}/2,* it will reach a point where the pipeline unsupported span distance

*(L*will be a maximum. From here, if the sinkhole moves through a further horizontal distance

_{s_span})*X*to the right of the pipe it will again reach a point where the pipe will reach the limit of its safe unsupported span distance. Hereafter the pipe will again be able to support itself without collapse. Consider triangle TUV that lies on the horizontal plane that coincides with the centre line of the pipeline in Figure 12. It follows that:

_{fs}/2 The unsupported span distance, L_{u_span} can be determined from Equation (8):

Equation (12) is obtained through substitution of Equations (11) into (10), and simplifying by substituting for *h _{cover}* and

*L*using Equations (5) and (7) respectively, and solving for

_{s_span}*X*. Equation (12) is used to calculate the width of the strip within the sinkhole where a pipe would fail.

_{fs}where

X_{fs} | is the width of the pipeline failure strip in the sinkhole in m |

**The number of sinkhole events likely to occur along a length of pipeline within a dolomitic area **

Consider a pipeline of length *Ls*, traversing land where the sinkhole failure rate is *λ**s*. The number of sinkholes that are likely to occur along this length of pipeline can be calculated as follows:

where

Ns | is the number of sinkholes that may form along the section of pipeline considered per year |

X_{fs} | is the width of the pipeline failure strip in the sinkhole in m |

Ls | is the length of pipeline passing through the dolomitic area subject to sinkhole formation in m |

λs | is the sinkhole failure rate in the dolomitic area within which the pipeline is installed, expressed as the number of events / km^{2} / year |

**The simultaneous failure of parallel pipelines **

If a sinkhole occurs, it may cause the failure of more than one pipeline, since large water distribution pipelines often run side-by-side within pipeline servitudes. The spacing of adjacent bulk water distribution pipelines laid in parallel in a servitude are governed by factors such as the availability of land, diameter of the pipeline, width of the trench, local soil conditions, and required working strips. Table 8 provides a guideline in respect of requirements to be met in terms of the spacing of pipelines to facilitate installation, repair and maintenance work to take place (Turnbull 1996).

A trendline providing a good fit to the pipeline spacing data pertaining to the minimum spacing of parallel pipelines is given by Equation (14).

where

Y_{ps} | is the pipeline spacing of parallel pipelines (between their centre lines) in m |

D | is the diameter of the pipeline in m |

The number of parallel pipeline failures that may be associated with a particular sinkhole that has occurred, can be determined as follows:

where

n _{pp} | is the number of parallel pipeline failures that may be anticipated under specific soil conditions, size of sinkhole, pipe diameter, and the spacing of parallel pipelines |

**PIPELINE FAILURE RATE DUE TO SINKHOLES **

The work undertaken in the previous sections provides a number of important contributions to the proposed methodology, namely:

Equation (7) provides a relationship between the pipe diameter and the pipeline safe span distance.

Equation (8) expresses the diameter of the sinkhole as a function of the pipe diameter, pipeline unsupported span distance across the sinkhole, soil angle of internal friction, as well as the pipeline depth of cover.

Equation (9) provides a relationship between the diameter of the pipeline, the soil angle of internal friction and the sinkhole diameter that will cause pipeline failure.

It was shown that for steel pipelines larger than 600 mm in diameter, there is no need to consider sinkholes with diameters less than 15 m if the soil angle of internal friction is within a range 15º

<θ<45º. In such instances the middle three columns of Table 3 can be disregarded.Equation (12) allows the width of the strip within the sinkhole where a pipe would fail to be calculated.

Equation (13) makes it possible to estimate the number of sinkholes, for any given sinkhole size, that could cause pipeline failure along a certain length of pipeline.

Equation (15) makes it possible to determine the number of parallel pipeline failures that may be anticipated under specific soil conditions, size of sinkhole, pipe diameter, and the spacing of parallel pipelines

**The proposed methodology **

What remains is to develop a method to calculate the pipeline probability of failure that takes account of the sinkhole size and its probability of occurrence. The proposed methodology will utilise the dolomitic risk classification framework developed by Buttrick (2001) related to the relationships that exist between the eight inherent risk classes, the sinkhole sizes and the sinkhole failure rates, as reflected in Table 9. The methodology to be developed in this paper will be developed for a pipeline that falls within one inherent risk class only, while it will also only consider sinkholes greater than 15 m in diameter.

The number of sinkholes that may occur along the length of a pipeline falling totally within an inherent risk class is determined as follows:

The base sinkhole diameter is denoted by

D_{sh(15) }The sinkhole event failure rate, assuming that all sinkholes that will occur are equal to or greater than

Dis equal to_{sh(i5)},λs. The value ofλsis given in Table 9.For sinkholes equal to or greater than

D, where_{sh,}Dis larger in size than_{sh}D, the sinkhole event failure rate,_{sh(15)}λ, is calculated as follows:_{Dsh}

where

λ_{Dsh} | is the sinkhole event failure rate associated with sinkhole sizes equal to or greater than D_{sh} (_{Dsh} > D_{sh(15)}), expressed as the number of sinkholes occurring/km^{2}/year | |

λs | is the sinkhole event failure rate associated with sinkhole sizes equal to D_{sh(i5)}, expressed as the number of sinkholes occurring/km^{2}/year |

F(D_{sh}) | is the probability that sinkholes equal to or greater than D s h will occur | |

F(D_{sh(i5)}) | is the probability that sinkholes equal to or greater than 15 m will occur |

*F(D _{sh})* and

*F(D*should be based on sinkhole failure rate data in the area under investigation in order to obtain the relationship between the sinkhole size and its probability of exceedance. In this paper, the relationships will be based on the lognormal distribution function of which the parameters are provided in Table 2.

_{sh(15)}) The number of sinkholes larger than* Dsh *that may occur along the length of pipeline is determined using Equation (17).

where

NS_{(Dsh)} | is the number of sinkholes equal to or greater than Dsh associated with inherent risk class i, measured as the number of events along the length of pipeline Ls/ year | |

X_{fs(DSh)} | is the sinkhole failure strip width in m | |

Ls | is the length of pipeline Ls in m |

The pipeline failure rate associated with sinkhole formation is calculated using Equation (18).

where

λP_{Dsh} | is the pipeline failure rate, measured as the number of failures/km/ year |

**A HYPOTHETICAL EXAMPLE **

A hypothetical example will be used to illustrate the application of the methodology based on the following input parameters:

1. Pipelines installed in a dolomitic area where the sinkhole failure rate associated with sinkholes greater than 15 m in diameter will be equal to 5 events/km

^{2}/year.2. The soil angle of internal friction may be 20º, 40º, 60º or 80º respectively.

3. The pipeline depth of cover and the spacing of parallel pipelines are defined by Equations 5 and 14 respectively.

4. The installed pipelines will consist of the following diameters: 0,9 m; 1,2 m; 1,5 m; 1,8 m; 2,1 m and 2,4 m.

5. The probability that a specific sinkhole diameter will be exceeded is governed by the data presented in Table 2.

The methodology developed will be applied to the above input parameters, exploring different scenarios, in order to determine the following:

The critical sinkhole diameter that may cause different pipeline diameters to fail for a range of soil types.

The width of the pipeline failure strip within the sinkhole, associated with sinkholes larger than the critical sinkhole diameter.

The pipeline failure rate.

The maximum number of parallel pipeline failures that may occur.

**The critical sinkhole diameter causing pipeline failure **

The application of Equation 9, utilising the relevant input parameters, is used to plot the graphs shown in Figure 13. Figure 13 illustrates the following:

The critical sinkhole diameter causing pipeline failure increases as the pipe diameter increases.

For a specific pipeline installed, the larger the soil angle of internal friction, the smaller the sinkhole causing pipeline failure.

**The width of the pipeline failure strip within sinkholes larger than the critical sinkhole diameter **

The application of Equation 12, utilising the relevant input parameters, is used to plot the graphs of the sinkhole diameter versus the pipeline failure strip width for* D* = 0,9 m and *D* = 2,4 m diameter respectively, as shown in Figure 14. Figure 14 illustrates the following:

The width of the pipeline failure strip increases as the sinkhole diameter increases.

For a specific pipe diameter and a specific sinkhole that may cause pipeline failure, the width of the pipeline failure strip increases as the soil angle of internal friction increases.

If a specific sinkhole has formed, and for pipelines installed in an area with a similar soil angle of internal friction, the width of the failure strip increases as the pipe diameter decreases.

**The pipeline failure rate **

The relevant input parameters, the data obtained in respect of the pipeline failure strip width, as well as Equations 16 and 18, are utilised to determine the pipeline failure rate. The results are presented in Figure 15 for* D *= 0,9 m and* D* = 2,4 m diameter respectively. Figure 15 illustrates the following:

For a specific pipeline installed, the pipeline failure rate increases as the sinkhole diameter increases until it reaches a local maximum value, after which it again decreases.

For a specific pipe diameter and a specific sinkhole that may cause pipeline failure, the pipeline failure rate increases as the soil angle of internal friction increases.

If a specific sinkhole has formed, and for a pipeline installed in an area with a similar soil angle of internal friction, the pipeline failure rate increases as the pipe diameter decreases.

For each of the curves plotted in Figure 15, it is noted that the maximum failure rate that occurs is characterised by unique values in respect of the pipe diameter, the diameter of the sinkhole that causes pipeline failure, the soil angle of internal friction and the pipeline failure rate. As a result thereof, the characteristics of all such points were determined in respect of all pipe diameter and soil angle of internal friction value combinations. The results are presented in Figures 16 and 17 respectively. The significance of these two figures relate to the fact that for a given pipe diameter installed in an area with a particular soil angle of internal friction value, the maximum pipeline failure rate that may occur, as well as the associated sinkhole size that will cause it, can be predicted.

A question that arises is, how critical are these pipeline failure rates when compared to normal age-related pipeline failure rates. Nel (2010) provided a relationship to calculate the pipeline failure rate for welded steel pipes, based on an analysis of published steel pipeline failure rate data, and although it was derived based on pipeline failure rate data for pipelines __<__ 800 mm in diameter. This relationship is used to determine the failure rates, and it is also plotted in Figure 16. It is noted from Figure 16 that pipeline failure rates associated with sinkholes in dolomitic areas may exceed the pipeline failure rates associated with normal pipeline failure rates, even with the sinkhole failure rate taken as low as 5 events/km^{2}/year. The probability of welded steel pipeline failure rates associated with sinkhole formation is therefore a significant factor that should be considered when assessing bulk water distribution reliability.

**The maximum number of parallel pipeline failures that may occur within a sinkhole **

The relevant input parameters, the data obtained in respect of the pipeline failure strip width, as well as Equations 12, 14 and 15 are utilised to determine the number of parallel pipeline failures that may occur within the sinkhole failure strip area. The results are presented in Table 10. Table 10 highlights the following:

For a specific sinkhole causing pipeline failure, the number of simultaneous parallel pipeline failures decreases as the pipeline diameter increases.

For a specific sinkhole causing pipeline failure, the number of simultaneous parallel pipeline failures increases as the soil angle of internal friction increases.

**CONSIDERATION OF OTHER FACTORS THAT MAY INFLUENCE THE FAILURE OF PIPELINES IN DOLOMITIC AREAS **

There are other factors that may adversely affect the failure of pipelines within dolomitic areas, such as for instance:

The provision of services such as water pipes, sewer pipes, storm water pipes and the construction of buildings. Buttrick

et al(2008) provided data related to the failure of pipelines attributed to dolomitic ground movement events, with reference to a study area located to the immediate south and southwest of central Pretoria (Gauteng Province, South Africa). The data presented highlighted the fact that 98,9% of the events that occurred in the study area could be attributed to leaking water services.Pipelines installed in areas where the pipeline may be subject to variable differential settlement along its length will experience bending stresses at the edge and centre of the depression (American Lifelines Alliance 2001). The movement of the soil around the pipeline will subject the pipeline to additional tensile and compressive stresses at the edge and the centre of the depression. The induced stresses may cause failure of the pipeline, and high compressive stresses may cause buckling of pipelines with high D/t ratios, while tensile stresses may cause failure of the pipeline material and failure of welded joints. This may imply that a pipeline spanning an area formed by a sinkhole may fail sooner than predicted.

Schoning (1990) and Buttrick

et al(2008) reported that rainfall and inadequate storm water disposal tend to increase the frequency of sinkhole formation.The loss of water from a leaking pipeline may increase the sinkhole size and thus enhance pipeline failure.

The support conditions at the edge of the sinkhole may be unstable due to the potential development of shear zones around the collapsed sinkhole, thus increasing the unsupported pipeline distance.

Soil that may be present on top of the pipeline after the sinkhole has formed will increase the load carried by the unsupported pipeline and enhance pipeline failure.

Structures that may be built around the pipeline will add an additional load to an unsupported pipeline.

The control and channelling of scour water from pipelines will require detailed planning in dolomitic areas subject to sinkhole formation, since scour water may trigger sinkhole formation that could stimulate pipeline failure.

The dynamic loads that may be induced onto a pipeline during sinkhole formation may also be a significant factor that should be considered.

**CONCLUSIONS **

The assessment of the probability of sinkhole formation and the related sinkhole diameter draws on the substantial body of experience accumulated in South Africa in the past years. Specific data and guidelines had to be adapted for the purposes of this analysis, and were complemented with a structural and spatial analysis of a steel pipe intersecting a sinkhole.

The methodology developed provides a means to predict the probability of failure of welded steel pipelines installed in areas subject to sinkhole formation. The methodology relies on the soil angle of internal friction to be known, as well as the cumulative distribution function with respect to the size of sinkholes that may be expected in an area.

The workability of the approach was demonstrated with a hypothetical example that highlighted the following aspects worth noting:

Pipelines larger than 600 mm in diameter will most likely not fail as a result of sinkholes less than 15 m in diameter. This finding is significant since it leads to a considerable simplification of the application of the risk classification framework of Buttrick (2001).

The maximum pipeline failure rates for pipelines between 0,9 m to 2,4 m in diameter and 20º

<θ<80º, are associated with sinkholes in the range 18 m to 37 m in diameter.Pipeline failure rates associated with sinkholes in dolomitic areas may exceed the pipeline failure rates associated with normal pipeline failure rates, even with the sinkhole failure rate taken as low as 5 events/km

^{2}/year. The probability of welded steel pipeline failure rates associated with sinkhole formation is therefore a significant factor that should be considered when assessing bulk water distribution reliability.

The methodology developed is conservatively related to the fact that sinkholes may be more concentrated in areas where leaking wet services or ponding surface water are encountered. Furthermore, the pipeline safe span distance is based upon static forces that the pipeline is subjected to, while the dynamic forces associated with sinkhole formation are not considered. Further refinements to this methodology, as well as better and more reliable data, will undoubtedly ensure that the level of confidence related to estimating welded steel pipeline failure rates in dolomitic areas associated with sinkhole formation can be improved.

**REFERENCES **

American Lifelines Alliance 2001. Guidelines for the design of buried steel pipe, ASCE, 43. [ Links ]

AWWA (American Water Works Association) 2003. Steel pipe - A guide for design and installation. Manual of water supply practices, M11, 74. [ Links ]

Brink, A B A 1979.* Engineering geology of southern Africa,* Vol 1. Pretoria: Building Publications. [ Links ]

Brinkmann, R, Parise, M & Dye, D 2008. Sinkhole distribution in a rapidly developing urban environment: Hillsborough County, Tampa Bay area, Florida. *Engineering Geology,* 99(3-4): 169-184. [ Links ]

Bruno, E, Calcaterra, D & Parise, M 2008. Development and morphometry of sinkholes in coastal plains of Apulia, southern Italy. Preliminary sinkhole susceptibility assessment.* Engineering Geology,* 99: 198-209, doi:10.1016/j.enggeo.2007.11.017. [ Links ]

Buttrick, D B 1992. Characterisation and appropriate development of sites on dolomite, Unpublished PhD thesis, Pretoria: University of Pretoria. [ Links ]

Buttrick, D B, Van Schalkwyk, A, Kleywegt, R J & Watermeyer, R B 2001. Proposed method for dolomite land hazard and risk assessment in South Africa.* Journal of the South African Institution of Civil Engineering,* 43(2): 27-36. [ Links ]

Buttrick, D B, Watermeyer, R B, Gerber, A P, Pieterse, N & Trollip, N Y G 2008. The influence of urbanisation on sinkhole and doline formation, unpublished paper. [ Links ]

Department of Public Works, South Africa 2004. Appropriate development of infrastructure on dolomite. Guidelines for consultants, Report PW 344, www.info.gov.za [ Links ]

Galve, J P, Gutierrez, F, Remondo, J, Lucha, P & Cendrero, A 2009. Evaluating and comparing methods of sinkhole susceptibility mapping in the Ebro Valley evaporite karst (NE Spain). *Geomorphology,* 111(3-4): 160-172, doi:10.1016/j. geomorph.2009.04.017. [ Links ]

Goodings, G J, Waleed, A & Abdulla, B 2002. Stability charts for predicting sinkholes in weakly cemented sand over karst limestone.* Engineering Geology,* 65: 179-184. [ Links ]

Hyatt, J A & Jacobs, P M 1996. Distribution and morphology of sinkholes triggered by flooding following Tropical storm Alberto at Albany, Georgia, USA. *Geomorphology,* 17: 305-316. [ Links ]

Kaufmann, O & Quinif, Y 1999. Cover-collapse sinkholes in the "Tournaisis" area, southern Belgium. *Engineering Geology,* 52: 15-22. [ Links ]

Lurie, J 1977.* South African geology for mining, metallurgical, hydrological and civil engineering. *Johannesburg: McGraw Hill. [ Links ]

Nel, D T 2010.* Factors that may compromise bulk water distribution reliability.* Doctoral thesis, University of Johannesburg, South Africa. [ Links ]

Pilecki, Z & Baranowski, A 2006. Estimation of dimension of a rectangular type sinkhole activated by abandoned shafts. Warsaw: Institute of Geophysics of the Polish Academy of Sciences, Report M-29 (395), http://agp.igf.edu.pl, 18 July 2008, 16:10. [ Links ]

Reyneke, P L 2007. Rand Water pipeline data. Excel spreadsheet prepared by Planning Department, Rand Water Head Office, Johannesburg. [ Links ]

Schoning, W L 1990.Verspreiding van sinkgate en versakkings in die dolomietgebiede suid van Pretoria. (Distribution of sinkholes and subsidence in the dolomite areas south of Pretoria). Unpublished MSc dissertation, Pretoria: University of Pretoria, pp 111-112. [ Links ]

Sinkhole.org 2008. Types of sinkholes. www.sinkhole.org. [ Links ]

Spangler, M G & Handy, R L 1984.* Soil engineering,* 4th ed. New York: Harper & Row. [ Links ]

Tharp, T M 1999. Mechanics of upward propagation of cover-collapse sinkholes.* Engineering Geology,* 52: 23-33. [ Links ]

Thomas, B & Roth, M J S 1999. Evaluation of site characterization methods for sinkholes in Pennsylvania and New Jersey.* Engineering Geology,* 52: 147-142. [ Links ]

Turnbull, E D 1996. A study into Rand Water's pipeline avenues and servitudes. Report of the Rand Water (South Africa) Pipelines Department, Rand water, Johannesburg. [ Links ]

Waltham, T, Bell, F G & Culshaw, M G 2005.* Sinkholes and subsidence. - Karst and cavernous rocks in engineering and construction.* Berlin: Springer verlag. [ Links ]

**Contact details:**

PO Box 112 Meyerton 196"

South Africa

T: +27 (0)16 362 0772 F: +27 (0)86 550 3544

E: daclc@absamail.co.za

**Contact details:**

PO Box 524 Auckland Park 2006

South Africa

T: +27 (0)11 559 2148 F: +27 (0)11 559 2395

E: jhaarhoff@uj.ac.za

DANIE NEL obtained his bachelor and master's degrees in civil engineering at the then Ranc Afrikaans University and his D Ing degree at the University of Johannesburg. He spent his career in the municipal and bulk water supply environment and is currently an independent consultant. He has vast experience in many of the disciplines related to the provision, operation and maintenance of civil engineering infrastructure, bulk water supply infrastructure projects, risk assessment related to bulk water distribution, management, pipeline risk assessment in dolomitic areas, as wel as bulk sanitation services. |

JOHANNES HAARHOFF graduated as civil anc environmental engineer at the University of Stellenbosch and Iowa State University. Prof Haarhoff spent the first part of his career in municipal employ, civil construction, design and project management before joining the University of Johannesburg in 1990, where he eads the Water Research Group. This group focuses on both applied and fundamental aspects of drinking water treatment and supply |