Sealing the base and slopes of the landfill creates an impermeable sealing barrier, protecting the ground against the penetration of leachates and landfill gases into the lower layers of the ground and groundwater, as well as draining the resulting leachate to the treatment system. High-density polyethylene (HDPE) geomembranes are one of the synthetic materials used as artificial sealing barriers in municipal waste landfills. One of the disadvantages of HDPE geomembrane is its smooth surface, which results in a low value of interface shear strength obtained for multilayered liner systems. This fact was especially noticed after the slope-stability failure of a Class I hazardous waste landfill at Kettleman Hills in California [12, 16]. The failure was caused by insufficient shear strength between layers of mixed storage seals. The interaction of smooth geomembrane and compacted clay layer was characterized by a very low residual value of the interface friction angle equalled 8°. Nowadays, geomembranes with textured surfaces are produced to prevent slippage along phases with mixed seal systems. Factors that affect landfill stability can be divided into internal (geological) and external (geo-environmental). More detailed they can be assessed as [7]: engineering properties of waste, structural features of the waste body (leachate and landfill gas) and dynamic engineering geological processes (earthquake, rainfall, leachate, excavation, overloaded).
Important elements affecting the stability of the landfill are the geometric dimensions (height, width and inclination of the slope of the waste massif), as well as their physical and mechanical parameters. Municipal waste collected in landfills is a material that is very diverse in terms of morphology and density. Fresh municipal waste has a density in the range of 0.4−1.0 Mg/m3, while landfill waste has a density of 0.8−1.2 Mg/m3 [25]. The most common ranges of strength parameters are about 20−35° for the angle of internal friction and 15−40 kPa for the cohesion resistance [26]. The slightly different values of strength parameters are given by Dixon et al. [3]: 15−42° for the angle of internal friction and 0–28 kPa for cohesion intercept. The shear strength of municipal waste varies over time, which is mainly related to its compression and decomposition of organic substances. The period of about 1.5–3 years after the ending of deposition corresponds to a change in the intensity of the processes of bio-decomposition [6]. Attention should also be paid to the gradual decrease in the strength parameters of waste due to the progressing decomposition of municipal waste.
Many publications have been assigned to the issue of the stability of landfills, e.g. [2, 8, 9, 10, 13]. Both the limit equilibrium methods based on a cylindrical (circular) slip surface and the finite element method (FEM) can be used to analyze the stability of landfills. In the most used limit equilibrium methods, the factor of safety (F) determined from the ratio of the stabilizing resistances and the destabilizing effect of actions is to be greater than the permissible value of the stability factor, which in the case of landfills should be taken in the range of 1.2 to 1.3, depending on the importance of the facility and threats to the adjacent areas. The slopes of municipal landfills with the factor
The construction of municipal waste landfills in Poland is currently regulated by the Act of 14 December 2012 on waste and the Regulation of the Minister of the Environment of 30 April 2013 on landfills [15, 19]. The regulation [15] stated that the operation of the landfill should ensure “geotechnical stability of the stored waste”, however, the method of ensuring this has not been defined. The slopes of the landfill in the post-operational phase should also be subject to the assessment of stability “determined by geotechnical methods”.
The aim of the work is to verify the stability of municipal waste stored in a landfill of a specific structure, assuming variables related to the waste massif, such as height, width of the crest and slope inclination of the massif. Two different materials were considered as a mineral sealing layer: compacted highly plastic clay and compacted fly ash that meets the conditions of the material for the construction of the sealing layer [20, 23]. The analysis may be helpful in determining a procedure for waste placement storage.
In the ground, under the landfill and its side walls, there should be a natural geological barrier in the form of a continuous layer of soil with a permeability coefficient k ≤ 10−9 m/s [2, 11, 14, 25]. In the case of the absence of a suitable natural geological barrier, an artificial barrier is made. The natural or artificial barrier is accompanied by a synthetic geomembrane. The base of the municipal waste landfill and its slopes are also equipped with a drainage system for leachate. According to [15], the minimum thickness of the natural geological barrier should not be less than 1.0 m, and the artificial mineral barrier should be at least 0.5 m thick. Drainage layers are designed from soil materials with a permeability coefficient
The landfill was assumed as a sub-level in the excavation, where the maximum height of the waste is equal to the height of the excavation slope. The slope of the excavation is made of fine sand. Variable geometrical parameters of the municipal waste massif were assumed, such as the height of the waste massif
Scheme of municipal waste storage at the landfill
The sealing layers of the base and slopes of the municipal waste landfill were adopted in accordance with applicable legal regulations and literature recommendations. The drainage layer is planned of medium-dense sand with 0.5 m, while the mineral sealing layer with 1 m is made of compacted stiff highly plastic clay. It should be noted that the optimum moisture content (
Cross-section through a single-sealing layer of slope and base of the landfill made of compacted clay or compacted fly ash
The values of geotechnical parameters of municipal waste and geosynthetic materials forming the landfill base were taken after [13]. The strength parameters of the synthetic layers were given as interface contact parameters, which were presented as peak strength values at maximum shearing resistance. Earlier authors’ research [18, 24] showed that the main failure mechanism took place in the base of the filled landfill, so peak values were used for calculation. The peak strength parameters are generally used in landfill base stability analyses when the residual values are used for the calculation of the stability of multilayer surface sealings [17]. Sliding resistance is not taken into consideration. Interface unit weights were taken as an arithmetic mean of the weights of two adjacent materials. The fly ash parameters were given after [21, 22] for the material characterized by the lowest hydraulic conductivity that was established at
Parameters of materials used for calculations of waste landfill stability
Layer number | Material | γ (kN/m3) | ||
---|---|---|---|---|
I | Municipal waste | 10.20 | 30.0 | 3.0 |
II | Medium sand |
16.68 | 33.6 | – |
III | Medium sand + Non-woven geotextile | 9.02 | 27.0 | 14.0 |
IV | Non-woven geotextile + Textured geomembrane HDPE | 5.29 | 24.0 | 0.0 |
V | Textured geomembrane HDPE + Compacted clay | 14.91 | 19.0 | 9.3 |
VI | Compacted clay |
20.60 | 17.5 | 30.1 |
Explanation:
Parameters of materials used for alternative calculations of waste landfill stability
Layer number | Material | |||
---|---|---|---|---|
V | Textured geomembrane HDPE + Compacted fly ash | 11.77 | 12.0 | 10.0 |
VI | Compacted fly ash R=0.97 | 14.32 | 40.0 | 42.0 |
Explanation:
The soil and water conditions in the subsoil were assumed to be simple. Near-surface formations are non-cohesive sandy soils in the form of medium sands at medium dense, also lying deeper in the subgrade. Below the non-cohesive formations, there are glacial clayey soils in the form of stiff sandy clays and clays. The presence of groundwater and leachate levels is not assumed in the analyzed subsoil.
Landfill slope stability is typically assessed using limit equilibrium methods. The most used methods are the classic limit equilibrium methods: Fellenius, Bishop, Janbu, or Morgenstern-Price. When performing calculations according to the recommendations of Eurocode 7 [4], it should be considered that the Eurocode imposes the assumption of horizontal forces between vertical stripes, which excludes the use of the Fellenius method. In the Fellenius method, zero shear and normal forces are assumed between the calculation blocks, which results in lower values of the obtained stability factors. Additionally, the heterogeneity of municipal waste deposited in the landfill increases the range of generated errors, so the Fellenius method can only be used for an approximate forecast of the stability of landfill slopes [8].
The most unfavourable circular slip surfaces of 45 construction variants of municipal waste massifs at two different constructions of landfill slope and base sealing layers, were analyzed. The considerations were carried out according to approach 3 (DA3) of Eurocode 7 [1, 4], approved according to the National Annex for checking the state of equilibrium (stability) and determining the degree of utilization. Stability calculations were also made considering the values of safety factors, i.e., using the characteristic values of parameters and actions.
The value of the degree of utilization (utilization factor) for the ultimate limit state GEO according to Eurocode 7 [1, 4] is given by the formula (1):
Comparatively, the results are presented as values of the factor of safety (F):
The structure stability analysis was performed using the GEO5 numerical program (Slope Stability module), considering the limit equilibrium methods: Fellenius/Petterson, Bishop, Spencer, Janbu and Morgenstern-Price, assuming a circular slip surface. The calculations were carried out several times, looking for the slip surface with the lowest factor of safety, called critical slip surface [5].
The calculation results are presented, depending on the geometry of the slope and the calculation method in Table 3 for compacted clay as a material in sealing, and Table 4 – for fly ash as a part of sealing. The results are shown as the degree of utilization for the limit state (
Percentage utilization for the limit state
Geometrical parameters | Percentage utilization |
||||||
---|---|---|---|---|---|---|---|
Bishop | Fellenius/Petterson | Spencer | Janbu | Morgenstern-Price | |||
α=20° | B=10 m | H=5 m | 94.1/1.33 | ||||
H=10 m | |||||||
H=30 m | |||||||
H=50 m | 72.3/1.73 | 73.0/1.71 | 71.6/1.75 | 71.4/1.75 | 71.3/1.75 | ||
B=50 m | H=5 m | 94.1/1.33 | |||||
H=10 m | |||||||
H=30 m | |||||||
H=50 m | 70.8/1.77 | 72.3/1.73 | 70.8/1.76 | 70.8/1.76 | 70.8/1.76 | ||
α=25° | B=10 m | H=5 m | 94.1/1.33 | ||||
H=10 m | |||||||
H=30 m | |||||||
H=50 m | 87.7/1.42 | 89.1/1.40 | 87.3/1.43 | 87.2/1.43 | 87.2/1.43 | ||
B=50 m | H=5 m | 94.1/1.33 | |||||
H=10 m | |||||||
H=30 m | |||||||
H=50 m | 88.8/1.41 | 90.9/1.38 | 88.9/1.41 | 88.9/1.41 | 88.9/1.41 | ||
α=30° | B=10 m | H=5 m | 94.1/1.33 | ||||
H=10 m | |||||||
H=30 m | 97.0/ |
98.4/1.27 | |||||
H=50 m | – | – | – | – | – | ||
B=50 m | H=5 m | 94.1/1.33 | |||||
H=10 m | |||||||
H=30 m | 102.0/1.22 | 105.4/1.19 | 102.1/1.22 | 102.1/1.22 | 102.1/1.22 | ||
H=50 m | 107.5/1.16 | 110.3/1.13 | 107.6/1.16 | 107.6/1.16 | 107.6/1.16 | ||
α=45° | B=10 m | H=5 m | 99.5/1.26 | 106.2/1.18 | 99.8/1.25 | 99.8/1.25 | 99.2/1.26 |
H=10 m | 133.6/0.94 | 137.7/0.91 | 133.1/0.94 | 132.8/0.94 | |||
H=30 m | – | – | – | – | – | ||
H=50 m | – | ||||||
B=50 m | H=5 m | 105.3/1.19 | 110.4/1.13 | 105.7/1.18 | 105.2/1.19 | ||
H=10 m | 126.5/0.99 | 132.7/0.94 | 127.0/0.98 | 126.9/0.98 | 127.0/0.98 | ||
H=30 m | 158.8/0.79 | 165.4/0.76 | 159.2/0.79 | 159.3/0.78 | 159.2/0.79 | ||
H=50 m | 171.5/0.73 | 178.4/0.70 | 171.8/0.73 | 171.8/0.73 |
Percentage utilization for the limit state (
Geometrical parameters | Percentage utilization |
||||||
---|---|---|---|---|---|---|---|
Bishop | Fellenius/Petterson | Spencer | Janbu | Morgenstern-Price | |||
α=30° | B=10 m | H=5 m | 94.1/1.33 | ||||
H=10 m | |||||||
H=30 m | 97.3/ |
98.4/1.27 | |||||
H=50 m | – | – | – | – | – | ||
B=50 m | H=5 m | 94.1/1.33 | |||||
H=10 m | |||||||
H=30 m | 102.0/1.22 | 105.4/1.19 | 102.1/1.22 | 102.1/1.22 | 102.1/1.22 | ||
H=50 m | 107.5/1.16 | 110.3/1.13 | 107.6/1.16 | 107.6/1.16 | 107.6/1.16 | ||
α=45° | B=10 m | H=5 m | 99.5/1.26 | 106.2/1.18 | 99.8/1.25 | 99.8/1.25 | 99.2/1.26 |
H=10 m | 133.6/0.94 | 137.7/0.91 | 133.1/0.94 | 132.8/0.94 | |||
H=30 m | – | – | – | – | – | ||
H=50 m | – | – | – | – | – | ||
B=50 m | H=5 m | 105.3/1.19 | 110.4/1.13 | 105.7/1.18 | 105.2/1.19 | ||
H=10 m | 126.5/0.99 | 132.7/0.94 | 127.0/0.98 | 126.9/0.98 | 127.0/0.98 | ||
H=30 m | 158.8/0.79 | 165.4/0.76 | 159.2/0.79 | 159.3/0.78 | 159.2/0.79 | ||
H=50 m | 171.5/0.73 | 178.4/0.70 | 171.8/0.73 | 171.8/0.73 |
In the vast majority of calculation cases, slightly lower values of the factors of safety
Analyzing the calculations made in accordance with DA3 of Eurocode 7 [4], it was found that with the assumed structure and geometric dimensions of the landfill, the waste mass can be considered stable at the storage height
Stability calculations considering the factor of safety may be more or less rigorous compared to the calculations according to Eurocode 7 (DA3), depending on the permissible value of the stability factor adopted. Assuming that the factor of safety should be
In the case of the analyzed slope inclination
It should be noted that the location of the critical slip lines varies depending on the geometrical dimensions of the waste body, and is generally independent of the adopted calculation method. However, the location of the slip line is affected by the type of layer sealing the slope and the base of the landfill. In the case of classic clay sealing, the slope inclination
Examples of slip surfaces generated with the GEO5 program for classic clay mineral sealing layer: a), b), c) the stability of the structure is preserved, d) the stability of the structure was not preserved
In the case of sealing with fly ash as the mineral layer of the sealing, similar courses of circular slip lines (Fig. 4) and mostly identical values of safety factors and utilization degrees were obtained, within the scope of calculations performed (Tables 3 and 4). Variation of the slip line was observed in the case of the slope inclination
Examples of slip surfaces generated with the GEO5 program for fly ash as a mineral sealing layer: a) the stability of the structure is preserved, b) the stability of the structure was not preserved
Calculations of the slope stability of municipal waste stored in landfills were made in accordance with the recommendations of approach 3 (DA3) for the ultimate limit state GEO of Eurocode 7 and by analyzing safety factors. Evaluation of slope stability using both methods is comparable if the permissible values of factors of stability are greater than 1.2. If
Municipal waste stored in a sub-level landfill is generally stable if the inclination of the waste slope is
The given values of geometrical parameters of the stored waste should be treated as indicative only, due to the large diversity of physical and mechanical parameters of municipal waste and their heterogeneity.
Using the Fellenius/Petterson method can lead to an underestimation of the factor of safety and an overstatement of the degree of utilization, and consequently to incorrect assessment of the safety of the structure. Other methods of assessing structure stability – the Bishop, Janbu, or Morgenstern-Price methods give comparable or the same results.
The use of various stability assessment methods (the Fellenius/Petterson, Bishop, Janbu, or Morgenstern-Price methods) leads to very similar circular slip lines. Changing the material of the mineral sealing layer can lead to a change in the course of the circular slip line.