Next Article in Journal
Experimental Study on a Reinforced Concrete Element to Extract the Durability Index with the Automated Visualization
Next Article in Special Issue
Enlarged Head Pressure-Dispersed Anchor Cable for Foundation Pit Engineering Purposes
Previous Article in Journal
Convolutional Neural Network Algorithm Trained with Anteroposterior Radiographs to Diagnose Pre-Collapse Osteonecrosis of the Femoral Head
Previous Article in Special Issue
Probabilistic Analysis of Ground Surface Settlement of Excavation Considering Spatial Variable Modified Cam-Clay Model Parameters
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Stability Conditions in Lignite Open Pits from Romania, Case Study: Oltețu Open Pit

Department of Environmental Engineering and Geology, Faculty of Mining, University of Petrosani, 332006 Petroșani, Romania
*
Author to whom correspondence should be addressed.
Appl. Sci. 2022, 12(19), 9607; https://doi.org/10.3390/app12199607
Submission received: 17 August 2022 / Revised: 20 September 2022 / Accepted: 20 September 2022 / Published: 24 September 2022

Abstract

:
The problem of the slope stability of open pit mines is one of constant interest and great importance, both during the period of operation, but also post-closure. The research focused on the Oltețu open pit (located in Berbesti Mining Basin, Romania) and was directed in such a way as to allow consideration in the stability analyses of natural (predisposing the investigated area to landslides) and anthropogenic (specific to open pit mining) factors and causes as well as their combined effect. The field investigations (observations on the technical condition of the slopes, discussions with the technical personnel from Oltețu open pit, and sampling) were completed with analyses and laboratory tests (physical–mechanical properties of rocks in the composition of the slopes). The stability analyses took into account different hypotheses related to the actual geometry of the working slopes, and a predictive analysis was also carried out for the forecasted evolution of the working fronts and lateral slopes. Following stability analyses, it was found that for most slopes, the stability reserve is insufficient to allow continuing lignite exploitation under safe conditions. The last part of the paper presents the solutions identified by the authors (adoption of new geometries of the working front and lateral slopes) in order to increase the stability reserve to a minimum acceptable level, which would allow the safe continuation of lignite extraction, and, in the end, some practical recommendations are briefly presented.

1. Introduction

Signing the European Green Deal by Romania, doubled by the National Recovery and Resilience Plan, means the abandonment of coal-based energy production by 2025 (at the latest by 2032, according to the first document), a fact that attracts the closure of the current lignite exploitation perimeters. However, as long as the lignite open pits are in operation, it is necessary to ensure the stability conditions of the work fronts and thus ensure the safety and security of the workplace.
The exploitation of lignite in Romania began in a systematic way in the second half of the last century, coal production having over time a significant weight in the national energy mix. In the last 10 years, the lignite extracted in Oltenia Mining Basin has contributed 15–25% to the total electricity production of Romania [1,2,3,4].
The stability of slopes is one of the most important problems that must be solved throughout the life of the open pits. Landslides are phenomena with a high probability of occurrence, even more so in the case of lignite open pits, which operate in predominantly soft rocks (sedimentary rocks and soils), such as clays, sands and marls. The failure and sliding of slopes in an open pit can have more or less serious consequences. Among the consequences of landslides produced over time in lignite open pits are the total or partial damage of machinery (mostly bucket wheel excavators), the degradation or destruction of the nearby habitats involved in the sliding process, the blocking of water courses, the destruction of communication routes, or even the death of people operating the machines or who live in nearby area.
To support the research topic, we mention two large landslides which involved the displacement of important quantities of material by sliding and generated important material losses (destruction of machinery and equipment) and supplementary costs related to the elimination of the effects of these slides: the sliding of the working front of the Alunu open pit in 2019 [5,6] (total destruction of a large-capacity excavator, the loading trolley and the transport lines); and the sliding of the Alunu open pit steps system in 2020 (which led to the cessation of productive activities for more than 2 years). Additionally, several instability phenomena were recorded over time that involved the movement of material from the dumps, from the working slopes or from the natural slopes, damage to households and communication and utility networks, and geostructural changes, mainly as a result of the influence of a complex of factors (geological, hydrogeological, hydrometeorological, climatic and anthropological), among which we mention the following: the sliding occurred on the southern part of the Berbești–Alunu mining perimeter from 2010; about 20 landslides occurred in the Panga mining perimeter (on the northern, southern, or north-eastern parts of the open-pit, and also on the slopes of the external dump) between 2008 and 2016; and the sliding of the West Berbești external dump from 2017 [6,7].
In fact, landslides that affect either the working or final slopes of lignite open pits, or the steps of internal or external dumps, are quite common phenomena (their occurrence and manifestation being signaled all over Europe, for example, in Poland, Czech Republic, Greece, Serbia, Kosovo, Germany, and Spain, but also in other parts of the world), being described in specialized literature (in terms of favoring or triggering factors and causes, unfolding mechanisms, damage caused, and preventive measures for the future) [8,9,10,11,12,13,14,15].
Although the triggers of major landslide can often be technically detected, their prediction and mitigation remain a major challenge, both for practicians and scholars.
The slope angle and height of the slopes are the main geometric elements that can influence their stability/instability [16,17,18,19,20,21,22,23,24,25]. Therefore, the design and compliance with the geometry of the in situ and dump slopes are extremely important elements for maintaining their stability both in the exploitation process and in the closing and reclamation phases of the open pits. In addition, the identification and control of the factors that influence the technical condition of the slopes must represent a permanent concern throughout the life of the open pit. The stability of slopes is influenced by numerous factors, which can be grouped into several categories (geological, mechanical, meteorological, etc.) [16,17,18,19,20,21,22,23,24,25,26,27]. It is difficult to assess the extent to which these factors contribute to triggering a landslide, which is why they must be analyzed and controlled as carefully as possible for each particular situation [25,27,28,29,30,31].
In the conditions where unforeseen sliding phenomena occur, the causes must be analyzed, and proper measures to eliminate them must be taken [29,32,33,34,35,36]. In any situation, monitoring the stability conditions is absolutely necessary, early warning systems being the solution to avoid sliding phenomena that can result in material and human losses [37,38,39,40,41,42,43]. Sometimes, the evolution of major landslides can be subdivided into several sliding models, offering the possibility to analyze the evolution, in terms of space and time, of similar ones [44].
Slope stability in mining perimeters is important for the security of the production processes, mining equipment, and especially for workers. The purpose of the stability analyzes is to evaluate the safety factor or to check the stability reserve of a slope and, based on it, to reduce the geotechnical risks that may occur by identifying and implementing the appropriate stabilization measures [17,19,27,28,45].
The specialized literature has been enriched over time with numerous stability analysis methods developed by researchers in the field, among which we mention the following:
  • The analysis methods based on the limit equilibrium concept (LEM) [17,18,21,22,23,24,25,26,27,28,34] considers the soil both a load and a resistance factor, being commonly used, as it gives accurate results on an acceptable level.
  • The finite element method [21,25,46] is a modern, precise method that can be applied to complex structures. It allows the consideration of hydrological and hydrogeological conditions, the elements of support and improvement of land, and drainage systems. It does not require the imposition of preliminary assumptions, and it provides indications regarding the evolution of the safety factor and deformations during the entire sliding process.
  • Methods based on 3D analyses [47,48,49,50] are applied in the case of slopes with complex geometries, variable structure, and partial loading, with a very well-defined breaking mechanism, and assume the addition of the third dimension in the 2D procedures, which leads to an increase in the number of unknowns. Three-dimensional methods are used quite rarely due to the fact that they have limited applicability.
  • Probabilistic methods [21,25,51,52,53] assume, in principle, the probabilistic analysis/processing of the physical and mechanical properties of the rocks and the probabilistic calculations by numerical or analytical methods.
  • Different combinations of the above-mentioned methods.
Dynamic programming methods can be used to determine the critical sliding surface and the stability factor and allow the incorporation of aspects such as spatial variability and time decay of the mechanical parameters. The results obtained by this method may differ substantially from those obtained by static models [54,55].
In this paper, slope stability analysis was performed using the limit equilibrium methods (Fellenius, Bishop, Simplified Janbu and Morgenstern–Price), as they offer satisfactory and credible results.

2. Materials and Methods

2.1. Problem Statement

Berbeşti Mining Basin (Figure 1) is localized within the Amaradia-Tărâia interfluve, the basin also bearing this name in the specialized literature. Geographically, the basin is located in the hilly area between Gorj and Vâlcea counties. The high-capacity mining exploitations (open pits) operating in Berbeşti Mining Basin are Olteţu, Alunu, Panga, and West Berbeşti. These exploitations have produced, since the date of commissioning and up to now, a quantity of approximately 90 million tons of lignite, and approximately 550 million cubic meters of waste rocks resulted to obtain this production. The average lignite/overburden ratio, of approximately 1/6.2 t/m3, is due to the location of these perimeters in a hilly area [56].
For the Oltețu open pit (located in the Alunu exploitation perimeter along with Alunu open pit), part of Govora Electric-Heating Plant (more precisely, the mining division), which is planned to be kept in operation at least until 2025, at the end of 2021, it was considered imperative to conduct a research on the stability of the actual working fronts and lateral slopes, as well as a predictive one, for the 2022 forecasted situation.
The major structural unit of which the lignite deposits in Northern Oltenia are a part is the premontane “Getic Depression”, which took over the function of a sedimentation area, evolving as such in the Paleogene and Neogene [59].
The premontane depression, formed during the Laramic movements, subsequently endured a sedimentation process that began in the Paleogene and during which several discontinuities with non-general character were noted. The sedimentary rocks of the Getic Depression belong to the Paleogene–Quaternary interval, and they are kilometers thick [60,61].
The studies carried out between 1985 and 1992 by I. Andreescu, N. Ţicleanu et al. [62,63,64] led to the division of the Dacian–Romanian deposits into three lithostratigraphic formations (Figure 2). This division represents a local fragmentation of the aforementioned deposits, the Jiu-Motru Formation incorporating the other two, the Berbeşti Formation and the Cîndeşti Formation, with lithostratigraphic differences specific to each area. It can also be noted the comprehensive nature of the formation of Jiu-Motru, which includes both Dacian and Romanian.
The exploitation of lignite in the Oltenia open pits includes a series of activities, such as excavation (mostly carried out in continuous flow with bucket wheel excavators, in working tiers with heights of 25–30 m); the transport of the excavated material (on belt conveyors); the deposition of waste rocks with dumping machines (into internal and external dumps); and the deposition of coal using storage machines.
In the particular case of the Oltețu open pit, at the upper part of the excavation slope, there is a layer of clay that is very unstable, sliding frequently and covering the excavation front (Figure 3).
For this reason, there are losses of time and productivity of the machines, which is reflected in the economic efficiency of the mining activity. Moreover, the continuous sliding of the clay layer can cause the lower layers to move as well, triggering a large-scale slide that can put the excavator in danger.

2.2. Field Observations and Laboratory Tests

The lithological sequence (see Figure 2) found in the structure of the working step in the Olteţu open pit is the one that appears in the Berbești mining basin and consists of yellow-brown clay, sandy marl, and lignite. The layers made up of these rocks alternate and have variable thicknesses within the analyzed perimeter (Figure 4).
The main rock types that participate in the structure of sterile intercalations and the overburden of the lignite deposit are found in various proportions: clay 65–72%; dust 13–21%; sand 5–14% [65].

2.2.1. The Yellow-Brown Clay

The yellow-brown clay layer is the first encountered from the surface of the land, and its thickness varies widely, between 5 and 15 m. In the bedding of this horizon, a layer of sandy marl is found.
After performing the grain size analyses, the specific curves were graphically represented and it was found that they have a similar shape, with reduced inclinations, characteristic of very uneven rocks (the resulting non-uniformity coefficient, U, being variable and much higher than 15), with a large diversity of particle diameters. The compaction capacity is increased, as the small particles tend to fill the gaps between the larger ones and, in theory, under these conditions, the porosity should be reduced. However, cracking, disaggregation, and repeated contraction, swelling, and compaction phenomena that affected the rocks overtime led to a relatively high porosity (49–58%). The high porosity allows the accumulation of large volumes of water, which leads to a significant reduction in the cohesion [65].
The meteorological conditions influence the moisture of the yellow-brown clay, which varies widely, respectively its behavior. The determinations showed that the natural humidity of this clayey rock is about 40% in the conditions of a dry period. When the moisture increases, the amount of adsorbed water also increases. When the clayey rock particles are surrounded by water, and they are separated from each other, the flow phenomenon can occur. This situation defines a fluid state when the rock, due to high moisture, has weak cohesion and practically does not resist shear stress. The swelling potential of this clay was observed while determining the flowing limit when, prior to flowing, a sudden swelling of the material occurred. The increase in volume is caused by a change in the state of stress and in turn, produces displacements or landslides. The flow limit corresponds to a moisture value of approximately 69–73%, which confirms the property of clay to adsorb water in large quantities. Thus, the clay can pass from a solid state to a fluid state, leading to the production of large flow–slide phenomena.
With the decrease in moisture (relatively high moisture, but unsaturated state), this type of rock becomes plastic, and its cohesion increases significantly. The lower the moisture, the harder the rock becomes. In these conditions, the yellow-brown clay contracts and cracks, forming fissures of different lengths and depths and favoring rock disaggregation. These characteristics determine the occurrence of cracks and fissures (transverse cracks and fissures being typically observed), causing the detachment of rocks and the appearance of sliding phenomena named “dry flows” (also described as rotational landslides/dry landslides, rock falls or rolling). These types of landslides occur successively and do not involve large volumes of rocks. In this way, the massif is gradually unloaded, the situation being favorable for its stability.
Locally, on the lateral slopes, sandy intercalations were observed. These intercalations have an inclination that concurs with the inclination of the land or the open pit slope, thus representing detachment or sliding surfaces of clayey rocks. The cracks and fissures that appear in dry periods favors water infiltration during rainy periods. Therefore, there are conditions that favor the appearance of negative geotechnical phenomena, regardless of the situation.
Depending on the degree of moisture (saturation coefficient), the yellow-gray clay is classified as very wet to saturated.
According to current standards [66,67,68] and specialized literature [18,19,21], depending on the value of the plasticity and consistency indexes (Ip higher than 35% and Ic between 0.69 and 0.88), the yellow-brown clay is a clay with very high plasticity and belongs to the class of consistent plastic clayey rocks.
Relying on classifications in the specialized literature [18,19,21], such as the classifications based on grain size composition, the consistency and plasticity index, and the ternary diagram, the classification and description of the rock were made. Thus, it was determined that it is a fat clay (a cohesive soil with relatively low sand, dust, and gravel content, more than 60% clay with a diameter lower than 0.005 mm, with a plasticity index of over 35%). The sand and gravel validate the presence of a sedimentation area.
In conclusion, the analyzed yellow-brown clay presents, due to its characteristics, a risk of sliding in both dry and wet conditions. Slides can occur in various forms/types, from drops, and rolls to flows, depending on the moisture.
According to the studies and observations made over time, the following aspects were found regarding the yellow-brown clay:
In the case of very low moisture, a gradual discharge of the massif takes place;
In the case of high moisture, the risk of landslides in which significant volumes of rocks can be entrained increases;
In the case of difficult conditions, such as a thick yellow-brown clay layer and the unfavorable geometrical elements of the slopes (high slope and large slope angle), the risk of sliding becomes greater, regardless of the moisture of the rocks.

2.2.2. The sandy Marls

In the area of the working fronts, there are three layers of sandy marls located, as follows: in the bed of the first lignite layer (open pit base), in the roof of the first lignite layer (respectively in the bed of the second lignite layer) with variable thickness between 0.3 and 1.5 m, and in the roof of the second lignite layer with a thickness between 3 and 4 m.
After carrying out the grain size analysis, it was found that the sandy marl is composed mostly of sand, dust, some elements of gravel, and clay (5–10%), having an average uniformity.
The natural moisture of the sandy marl varies from 16.93% to 56%, and it influences the behavior and the resistance characteristics in the sense of increasing the shearing resistance with the decreasing of the moisture. In case of high moisture, due to the content of sand and dust, the rock is easily breakable, and in cases of low and very low moisture, the rock is hard, compact, and it has a high cohesion (between 1.38 and 1.72 daN/cm2) [65].

2.2.3. The Lignite

The first and second lignite layers (I, II a and II b, see Figure 2) represent the main (with economic importance) lignite layers in the Oltețu perimeter by quantity. They are thick and extend over almost the entire mining field.
Within the Oltețu mining perimeter, the first layer varies between 2.3 and 3.2 m thick. Qualitatively and qualitatively, it represents the main layer of lignite [56].
The second lignite layer is complex from a lithological point of view. It consists of two coal banks (layers II a and b) separated by a thin intercalation of sandy marl rocks. The second layer varies from 2.6 to 5 m in thickness [56].
The two lignite layers are separated in the roof and in the bed by sandy marl layers with varying thicknesses.
Generally, the physical–mechanical properties of lignite are favorable in terms of slope stability, but their values also vary as a result of the partially woody structure and clay inclusions.
Field research also led to the following findings:
  • The geometric non-uniformity of the excavation fronts is caused by the morphology of the land, as the excavation process is completed in one step. This way of excavation leads to variations of the geometric elements, both on the width and on the depth (expansion) of the working front.
  • Phenomena of fragmentation and detachment of the rocks from the slopes in conditions of low moisture and plastic flow tendencies of the clay from the upper part of the step under conditions of its saturation.
  • Local manifestation of instability, affecting the entire height of the step, with variable extent depending on the structure of the massif and the hydro-meteorological conditions.
  • Under the conditions of the initial geometry (December 2021) of the excavation fronts, the risk of massive landslides, which would involve large volumes of rocks and affect the safety of machinery and personnel, is relatively low.
  • Taking into account the terrain morphology, the difficulties of ensuring a relatively uniform geometry of the excavation fronts will be accentuated with their advance.
  • In the areas with increased amplitudes (the western part of the perimeter), taking into account the functional technological parameters of the E-02 excavator (ERC 1400-30/7 type) and possible instability phenomena, the design of appropriate exploitation technologies must be imposed so that the geometry of the slopes ensures their stability.
In order to carry out the stability analyses, three sampling campaigns were carried out (22 October; 12 November, and 2 December 2021) from all the rock layers encountered, and the samples were analyzed within the Soil Mechanics Laboratory (LMP) of the University of Petroşani (Figure 5). To verify the results, some tests were demanded from the Călan GeoLogic laboratory, especially for the first layer from the surface, consisting of yellow-brown clay (Table 1).
Following the analyses, it was clear that the yellow-brown clay was the least resistant rock. It has the highest affinity for water, and therefore, the deterioration of these characteristics is likely under unfavorable hydro-meteorological conditions. In this situation, the state of tension in the massif changes, and structural deformations occur. This conclusion is also confirmed by field observations. Stability analyses were performed, taking into account these aspects.

3. Results and Discussion

As shown in the introduction, there are a number of factors that influence the stability of slopes. For the case study analyzed in the paper, the main factors that can be considered are the hydrostatic pressure generated by the water from rock pores, seismic shocks (the analyzed area being characterized by a seismic intensity of 8/MSK scale, the average return period being about half a century [70]) and the vibrations generated by machines.
Following the in situ documentation, it was found that the vibrations transmitted mainly by the E-02 excavator to the slope, under the conditions where the yellow-brown clay layer was at natural moisture (approximately 40%, after a period without of significant precipitations) favored its fragmentation and the appearance of “dry flow” type slides (rolling, falling of “lumps” of clay) characterized as rotational and successive slides. These types of slides involve small volumes of rocks, gradually unloading the massif, so the mining activity can continue under relatively safe and satisfactory conditions.
It is difficult to appreciate the influence of the vibrations on slope stability in conditions of high moisture of the rocks (mainly clayey rocks) when such simulations cannot be carried out in the laboratory [65].
The authors considered that the vibrations produced by the heavy tonnage machines (excavators) are much lower in intensity and amplitude compared to those originating from natural earthquakes such that they cannot trigger liquefaction-type slides or slides caused by the thixotropic behavior of the material.
The maximum value for the seismic acceleration in the studied area is a = 0.2 g, the return period of the earthquake is over 2 centuries (about 225 years according to [70]), and the exceeding probability is low. Knowing these and the fact that the vibrations contribute to a small degree to triggering a major landslide, in the performed stability analyses, the seismic acceleration was considered to be a = 0.1 g.

3.1. Analysis of the Stability of In Situ Slopes

The stability analyses were carried out in seven sections, including two longitudinal (L1 and L2) and five cross (P0, P1, P2, P4 and P5), and in each section, the position of the working fronts was taken into account at the level of December 2021 and quarters I-IV of 2022, which led to the analysis of 13 geotechnical profiles, drawn on the situation plan with the projection of the open-pit evolution (Figure 6). The stability analyses for the lateral slopes were performed only in the case of the slopes located on the southwest side, as slopes on the north-eastern side have heights that do not exceed 15 m, so they do not raise stability problems.
The cases of lateral slopes in sections P0, P1, and P5 were analyzed from the perspective of ascending slopes, considering that their height increases and substantial changes in the geometry of the steps are required. We mention that, at the request of the research team, the mining operator made available eight sections L1, L2, P0–P5, essential for the development of the stability study. Of the provided sections, the P3 cross section was not used, being outside the interest area.
The stability calculations were carried out in three variants, respectively:
  • Drained slopes, unaffected by external factors (Fs1);
  • Undrained slopes, where pore water pressure can induce landslides (Fs2);
  • Undrained slopes and under the influence of seismic shocks, caused by the vibrations generated by the excavator and conveyor belts and/or the occurrence of an earthquake (Fs3).
The results of the stability analyses are presented in Table 2 for the working front slopes in sections L1 and L2, and for lateral sections P2 and P4 (it is planned that in the IV quarter of 2022, the E02 excavator will be withdrawn in the approximately same position as December 2021, and the excavation will be carried out in the lateral slope, corresponding to cross sections P2 and P4, thus becoming the working front) taking into account the advance of the open pit, in Table 3, for the lateral slopes, P0, P1, and P5 (in accordance with the semester in which they materialize) and in Figure 7, Figure 8, Figure 9 and Figure 10.
The stability analyses were performed using the SLOPE software, a software specialized in geotechnical problems, considering that a possible landslide can occur along a curved sliding surface.
This software performs stability analyses based on limit equilibrium theory (LEM), the safety factor, or the stability coefficient being given by the ratio between the sum of the moments of the resistance forces and the sum of the moments of the forces generating the slide. For this purpose, based on the recommendations of the literature [17,18,21,34], the methods of Fellenius, Bishop, Simplified Janbu and Morgenstern–Price were used. The values presented in Table 2 and Table 3 are the lower ones as determined for each section and quarter.
As a result of carrying out the stability analyses, the following conclusions can be drawn:
  • The stability analyses were performed for the initial conditions of the slopes only in the L1 (December 2021) and L2 (December 2021) sections, using the stratigraphic columns provided by the mining operator, and confirm the situation observed during the field investigations. Thus, in the area of L1 section, for the configuration of the slope from December 2021, instability phenomena of the clay and marl layers were manifested. The stability analysis carried out led to a value of 1.01 for the stability factor under conditions of no external stress on the slope (limit equilibrium state) and to sub-unit values for the case where pore water pressure and/or the effect of vibrations were manifested. The slope in L2 section was stable in December 2021, and the analyses performed indicate a stability reserve of 22% under no external stress, equilibrium limit under the conditions of the manifestation of pore water pressure, and instability under the conditions of the simultaneous effect of water and vibrations (shocks).
  • The other sections, namely cross sections P0–P2, P4 and P5, present the hypothetical position of the lateral slopes at the end of the four quarters of 2022 (as they materialize).
  • The stability factor depends on the height, inclination, and structure of the slope, a significant influence on the physical state (static and dynamic) coming from the weight of the clay layer in the structure of the slope. As for the sliding surface, in conditions where the thickness of the clay layer on the upper part exceeds 9–10 m, it is transmitted only through this layer and, as a rule, the stability coefficient is sub-unitary, indicating clay instability phenomena. Under the conditions of a clay layer of less thickness, the sliding surface is transmitted either through the layers of clay and sandy marl or through the entire step.
  • Stability is ensured for slopes between 18 and 26 m high, and slope angles of up to 60°, in conditions of no external stress. For all other cases, namely the presence of water and seismic shocks (except for the slope in section L2, III quarter), the stability factor is sub-unitary or approaches the equilibrium limit.
Considering the advancement of work fronts in areas with an ascending slope, in most cases, it is necessary to build a system of steps, which can ensure, on the one hand, the technical possibility of excavation (considering the fact that the theoretical excavation height of the ERC 1400-30/7 bucket wheel excavator is of maximum 30 m), and on the other hand, maintaining stability conditions. Such situations are recorded in the case of cross-sections P5, P1 and P0, corresponding to the stage of the excavation works in the I, II, and III quarters of 2022, where the distance between the bed of the first lignite layer (open-pit base) and the surface of the land is 40, 41, and, respectively, 38 m. For these cases, the stability of the designed (forecasted) slopes (with the geometry resulting from the provided sections) was analyzed, and the general conclusion is that such slopes would be unstable in all the three calculation hypotheses of the stability factor.

3.2. Determination of the Stable Geometry of the Working Front and Lateral Slopes

To ensure the stability conditions of the slopes, several geometries of the working steps and the lateral slopes were analyzed, determining their height and inclination for an imposed value of the stability factor.
In order not to affect in an unjustified manner the excavator’s efficiency indicators, the dimensioning took into account the fact that the working slopes have a short lifetime, so a stability reserve of 5–10% can be considered sufficient, and in the case of side slopes, which remain in place from a few months up to a year, a stability reserve of 20–25% was considered sufficient.

3.2.1. Stable Geometry of Working Slopes

For this case, section L1 was taken into account (Figure 11), and the following values of geometrical elements were established in order to ensure a stability reserve of 5%:
H = 25 m;
α = 52°;
Yellow-brown clay layer thickness: 9 m;
Fs1 = 1.05.
The thickness of the yellow-brown clay layer is lower than 10 m (the maximum allowed height), so it is not necessary to create sub-steps in this case.

3.2.2. Stable Geometry of Lateral Slopes

For this case, section P1 was taken into account, where the thickness of the rock package between the lignite layer no. I bed (elevations +345 to +334 m above sea level) and the land surface (elevations +350 to 380 m above sea level) records maximum values, respectively, 41–46 m at a distance of 70–90 m from the outcrop area of the lignite layers.
In these conditions, it is necessary to adopt a technology of excavation in steps and sub-steps to ensure stability conditions.
The first step (from the base of the open pit) up to +370 m surface level should be of maximum 26 m (crossing the lignite layers, the sandy marl between and above the lignite layers and the yellow-brown clay at the top on maximum thickness of 7.6 m).
The second step, up to +380 m surface level, will be excavated exclusively in the yellow-brown clay layer, and will have a maximum height of 10 m.
When the yellow-brown clay layer exceeds 10 m thickness (above the elevation of +380 m surface level), it is necessary to divide the second (upper) step into sub-steps so as not to exceed the maximum allowed height of 10 m. The upper sub-step will have a variable thickness, between 3 and 10 m, as the excavation advances toward the planned configuration (in the final stage, the slope will be divided in two steps, one of 26 m (the lower one) in height and one of 20 m (the upper one) divided in two sub steps of 10 m in height each).
a. Lower step with lignite layers (Figure 12):
H = 26 m;
α = 55°;
Yellow-brown clay layer thickness: 7.6 m;
Fs1 = 1.25.
b. Upper step and sub-steps in yellow-brown clay (Figure 13):
Height of sub-steps H = 10 m;
The height of the sub-step system (upper step) Hgen = 20 m;
The width of the intermediate berm B = 60 m;
The width of the safety berm b = 30 m;
Inclination of sub-steps slopes αind = 45°;
General inclination of slope αgen = 25°;
The stability factor for sub-steps Fs1ind = 1.15;
The stability factor for general slope (sub-step system) Fs1gen = 1.26.

3.3. Recommendations for Ensuring the Stability Reserve

  • Framing the excavation fronts in the geometric parameters that ensure the necessary stability reserve, depending on the type and lifetime of the slopes, according to the elements resulting from the determination of the stable geometry.
  • Preliminary uncovering (stripping) of areas with pronounced non-uniformity of the land morphology, thus ensuring the geometry of the slopes in conditions of stability.
  • Taking into account the tendency of the yellow-brown clay layer to disaggregate and the possibility of plastic failure, manifested by phenomena of “flow” on the slope both in a dry state as well as in the wet state, the maximum thickness in the structure of the excavation steps must be of maximum 10 m.
  • A permanent monitoring of the deformation phenomena of the front and lateral slopes (cracks and fissures, subsidence or swelling areas, and other signs that reflect changes in the state of tension in the massif).
  • Detailing the geological structure from the advance of the open pit, either by recovering the data of the initial exploration drillings or by executing new drillings that provide more precise information regarding the stratigraphy, tectonics, and hydrogeology of the area.

4. Conclusions

As a result of this research, the following observations were made regarding the characteristics and behavior of rocks:
  • The lignite layers in the perimeter of the Oltețu quarry are buried in sedimentary formations, being intercalated by clayey and marly rock layers;
  • According to the grain size analyses, these rocks contain variable amounts of sand (generally fine sand), dust, and gravel;
  • The geometry of the working fronts is uneven, being conditioned by the geomorphology of the region;
  • The variation of the physical and mechanical properties of rocks is accentuated and highly influenced by the water content (moisture);
  • The mechanical characteristics (cohesion and internal friction angle) decrease with increasing moisture, being dependent on hydro-meteorological conditions;
  • The great variety of physical–mechanical characteristics required statistical processing and the selection of values considered representative for the stability analyzes;
  • Considering the reduced resistance characteristics and the plastic behavior of the yellow-brown clay, the first layer from the land surface, it can be stated that it determines the occurrence of instability phenomena, such as sliding, disaggregation, dry flow, and plastic failure, regardless of the rock moisture;
  • The stratigraphy is variable and not known precisely until the time of excavation (relatively small number of prospecting drillings on the advancement direction).
The stability analyses were carried out for the areas of interest of the mining operator, respectively, for the forecasted slopes in 2022 according to the preliminary production plan. The results of the analyses show that the most unfavorable stability conditions occur in the case of front structures where the clay layer at the top of the slopes exceeds 10 m in thickness. In the presence of other disturbing factors (precipitation, vibrations, and seismic shocks), in most cases, large-scale instability phenomena can be triggered.
It is recommended to provide slope angles of maximum 52° and step heights of 25–26 m for front slopes. In the case of the lateral slopes, in order to ensure the reserve of stability of the step configured in the clay layer, it is necessary to divide it into sub-steps of 10 m each, and to reduce the slope angle to the value of 45°, resulting in a general slope angle of 25°.
By studying the behavior of the rocks that appear in the geometry of the working and final slopes depending on the thickness of the layers, their location in the lithological structure, and the existing meteorological conditions allowed the authors to draw some general conclusions regarding the necessary measures that must be taken in order to safely continue the lignite extraction activity in Oltețu open pit.
However, as the general geology (in particular the lithological structure), as well as other natural and anthropogenic factors, is similar throughout the entire Berbești Mining Basin, these results can provide valuable information to mining developers and engineers operating the other three open pits on how to deal with the top layer of clay so as to avoid unwanted events (large landslides capable of producing significant economic or even human losses).
The present paper has the merit of addressing a geotechnical phenomenon in the mining perimeters of Oltenia that was not particularly taken into account in the stage of designing the geometry of the open-pit slopes, namely the dry or wet flow of the yellow clay from the upper part of the open pit. Such flows generate changes in the state of stresses and strains and can lead to the triggering of slope failure along the entire height. Thus, detailed laboratory studies were carried out and, following stability analyses, the main failure mechanisms were defined. Based on the results obtained, it was proposed to modify the geometry of the working steps and implicitly the excavation technology, so that the lignite exploitation can be carried out in safe conditions, even if the lithological structure of the slope is less known. The results obtained are valuable both for mining specialists (the risk of slope failure will be greatly reduced), but also for students who are preparing in the field of mining and geotechnics (by including these results in the academic training process).
An interesting, and at the same time challenging, direction of continuing this research is represented by the introduction within the existing models of structural weakening coefficients of the massifs, when we are dealing with rocks with a behavior similar to that of the yellow-brown clay described in this paper, or even introducing this type of behavior into dynamic simulations.
Therefore, the paper makes a definite contribution to mining practice (especially the extraction by surface mining techniques of coal deposits, confined in sedimentary formations, weakly cohesive and prone to alteration and modification of resistance characteristics under the influence of external factors, mainly hydro-meteorological), but it also serves as a useful study material and starting point for improving existing software and stability models.

Author Contributions

Conceptualization, M.L., F.F. and I.-M.A.; methodology, M.L., F.F. and I.-M.A.; software, M.L., F.F. and I.-M.A.; validation, M.L., F.F. and I.-M.A.; formal analysis, M.L., F.F. and I.-M.A.; investigation, M.L., F.F. and I.-M.A.; resources, M.L., F.F. and I.-M.A.; data curation, M.L., F.F. and I.-M.A.; writing—original draft preparation, M.L., F.F. and I.-M.A.; writing—review and editing, M.L., F.F. and I.-M.A.; visualization, M.L., F.F. and I.-M.A.; supervision, M.L., F.F. and I.-M.A.; project administration, M.L., F.F. and I.-M.A.; funding acquisition, M.L., F.F. and I.-M.A. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Informed Consent Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Fodor, D.; Diaconescu, P. Evolution of lignite exploitation in the coal pits of Rovinari mining basin. Min. Rev. 2012, 18, 22–27. Available online: https://www.upet.ro/revistaminelor/archive/rmpdf2012/nr3en.html (accessed on 20 July 2022).
  2. Huidu, E. Exploitation of Coal Deposits through Open Pit Mining Works; Măiastra: Târgu Jiu, Romania, 2012; pp. 34–35. (In Romanian)
  3. Apostu, I.M.; Lazăr, M.; Faur, F. Implications of coal disappearance from the energy mix in Romania. In Proceedings of the 10th edition of the International Symposium on Occupational Health and Safety—SESAM 2021, Petrosani, Romania, 18 October 2021; Găman, G.A., Ed.; 2022; Volume 354. [Google Scholar] [CrossRef]
  4. Tomescu, C.; Cioclea, D.; Boanta, C.; Morar, M. Restructuring of coal mining in Romania between the climate crisis and the energy transition. In Proceedings of the 10th edition of the International Symposium on Occupational Health and—SESAM 2021, Petrosani, Romania, 18 October 2021; Găman, G.A., Ed.; 2022; Volume 354. [Google Scholar] [CrossRef]
  5. Lazăr, M.; Faur, F.; Rotunjanu, I.; Rada, C.; Apostu, I.M. The influence of hydrometeorological conditions on the stability of slopes from Alunu mining perimeter. An. Univ. Petroşani Min. Eng. 2020, 21, 86–99. Available online: https://www.upet.ro/annals/mining/pdf/Annals%20Mining%202020%20v1.pdf (accessed on 2 July 2021).
  6. Rada, C.; Lazăr, M.; Faur, F.; Apostu, I.M. Comparative study regarding the causes of land instability phenomena in the area of the Getic Subcarpathians. In Proceedings of the 9th edition of the International Multidisciplinary Symposium “UNIVERSITARIA SIMPRO 2021”: Quality and Innovation in Education, Research and Industry—The Success Triangle for a Sustainable Economic, Social and Environmental Development, Petrosani, Romania, 27–28 May 2021; Lazăr, M., Faur, F., Popescu-Stelea, M., Eds.; 2021; Volume 342. [Google Scholar] [CrossRef]
  7. Faur, F.; Lazăr, M.; Apostu, I.M.; Rada, C.; Moisuc-Hojda, D. Assessment of the Stability State and the Risk of Landslides within Berbești Mining Basin (Romania) Post Closure. Inż. Min.-J. Polish Min. Eng. Soc. 2022. submitted. [Google Scholar]
  8. Steiakakis, E.; Kavouridis, K.; Monopolis, D. Large scale failure of the external waste dump at the “South Field” lignite mine, Northern Greece. Eng. Geol. 2009, 104, 269–279. [Google Scholar] [CrossRef]
  9. Ylber, M. Evaluation of landslides and engineering measures on lignite open pit slope in south east Sibovc-Kosovo Coal Basin. In Proceedings of the Landslide and Flood Hazard Assessment, 1st Regional Symposium on Landslides in the Adriatic-Balkan Region and 3rd Workshop of the Japanese-Croatian Project “Risk Identification and Land-Use Planning for Disaster Mitigation of Landslides and Floods in Croatia”, Zagreb, Croatia, 6–9 March 2013; Arbanas, S.M., Arbanas, Z., Eds.; City of Zagreb Emergency Management Office: Zagreb, Croatia, 2013; Volume 1. [Google Scholar]
  10. Bednarczyk, Z. Slope instabilities in Polish open-pit mines. In Landslides and Engineered Slopes. Experience, Theory and Practice, 1st ed.; CRC Press Taylor and Francis Group: New York, NY, USA, 2016; 9p. [Google Scholar] [CrossRef]
  11. Steiakakis, C.; Apostolou, E.; Papavgeri, G.; Agioutantis, Z.G.; Schilizzi, P. Large moving landslide inside a lignite mine in northern Greece. In Landslides and Engineered Slopes. Experience, Theory and Practice; CRC Press Taylor and Francis Group: New York, NY, USA, 2016. [Google Scholar] [CrossRef]
  12. Steiakakis, E.; Thalassinakis, J.; Agioutantis, Z.G.; Galetakis, M. Failure process of the external waste dump of a lignite mine using FEM. In Proceedings of the 13th International Symposium Continuous Surface Mining (ISCSM 2016), Belgrade, Serbia, 12–14 September 2016. [Google Scholar]
  13. Baumann, A.B.G. Landslides from Space—Amyntaio Lignite Coal Mine Landslide (10th June 2017); Zenodo, European Organization for Nuclear Research: Meyrin, Switzerland, 2017. [Google Scholar] [CrossRef]
  14. Vanneschi, C.; Eyre, M.; Burda, J.; Žižka, L.; Francioniac, M.; Coggan, J.S. Investigation of landslide failure mechanisms adjacent to lignite mining operations in North Bohemia (Czech Republic) through a limit equilibrium/finite element modelling approach. Geomorphology 2018, 320, 142–153. [Google Scholar] [CrossRef]
  15. Chen, T.; Shu, J.; Han, L.; Li, B. Landslide mechanism and stability of an open-pit slope: The Manglai open-pit coal mine. Nat. Haz. 2022, submitted. [CrossRef]
  16. Hoek, E.; Bray, J.W. Rock Slope Engineering, 3rd ed.; Institute of Mining and Metallurgy: London, UK, 1981; pp. 67–69. [Google Scholar]
  17. Wyllie, D.C.; Mah, C.W. Rock Slope Engineering. Civil and Mining, 4th ed.; Spon Press—Taylor & Francis Group: New York, NY, USA, 2004; pp. 51–54. [Google Scholar]
  18. Duncan, J.M.; Wright, S.G.; Brandon, T.L. Soil Strenght and Soil Stability, 2nd ed.; John Wiley & Sons, Inc.: Hoboken, NJ, USA, 2014; pp. 81–132. [Google Scholar]
  19. Derek, H.C. Landslides in Practice: Investigation, Analysis, and Remedial/Preventative Options in Soils; John Wiley & Sons, Inc.: Hoboken, NJ, USA, 2005; pp. 170–219. [Google Scholar]
  20. Soltanmohammadi, M.; Osanloo, A.S.; Malekzadeh, S.B. Selection of practical bench height in open pit mining using a multi-criteria decision-makingsolution. J. Geol. Min. Res. 2010, 2, 48–59. [Google Scholar]
  21. Stănciucu, M. Stability of Natural and Artificial Slopes; Technical: Bucharest, Romania, 2018; pp. 11–150. (In Romanian) [Google Scholar]
  22. Abramson, L.W.; Lee, T.S.; Sharma, S.; Boyce, G.M. Slope Stability and Stabilization Methods; John Wiley & Sons, Inc.: New York, NY, USA, 2002; pp. 378–382. [Google Scholar]
  23. Huang, Y.H. Slope Stability Analysis by the Limit Equilibrum Methods: Fundamentals and Methods; ASCE Press: New York, NY, USA, 2014; pp. 229–264. [Google Scholar]
  24. Salgado, R. The Engineering of Foundation, Slopes and Retaining Structures, 2nd ed.; CRC Press: Boca Raton, FL, USA, 2022; pp. 894–917. [Google Scholar]
  25. Dana, H.; Kakaie, R.; Rafiee, R.; Bafghi, A. Effects of geometrical and geomechanical properties on slope stability of open-pit mines using 2D and 3D finite difference methods. J. Min. Environ. 2018, 9, 941–957. [Google Scholar]
  26. Luo, S.; Huang, D.; Peng, J.; Tomás, R. Influence of permeability on the stability of dual-structure landslide with different deposit-bedding interface morphology: The case of the three Gorges Reservoir area. China. Eng. Geol. 2022, 296, 106480. [Google Scholar] [CrossRef]
  27. Kumar, M.; Raj Kumar, M. Factors Affecting Slope Stability of an Opencast Mine: A Brief Study. Int. J. Innov. Sci. Res. Technol. 2022, 7, 878–882. [Google Scholar] [CrossRef]
  28. Wang, C.Y.; Zhang, F.; Han, W.D. Sensitivity Analysis of Slope Stability Influence Factors Based on BP Neural Network. Adv. Mat. Res. 2014, 1010, 1544–1547. [Google Scholar] [CrossRef]
  29. Zevgolis, I.; Deliveris, A.; Koukouzas, N. Slope failure incidents and other stability concerns in surface lignite mines in Greece. J. Sustain. Min. 2019, 18, 182–197. [Google Scholar] [CrossRef]
  30. Ye, S.; Huang, A. Sensitivity Analysis of Factors Affecting Stability of Cut and Fill Multistage Slope Based on Improved Grey Incidence Model. Soil Mec. Found. Eng. 2020, 57, 8–17. [Google Scholar] [CrossRef]
  31. Han, D.; Danqing, S. Investigation of failure prediction of open-pit coal mine landslides containing complex geological structures using the inverse velocity method. Nat. Hazards 2022, 111, 2819–2854. [Google Scholar] [CrossRef]
  32. Islam, S.; Hosain, M.; Islam, M.S.; Ahmed, D.; Hoque, F.; Karim, S.U.; Islam, A. Disaster due to slope failure in the Pahartoli area of the Chittagong city, Bangladesh. Int. J. Sci. Eng. Res. 2015, 6, 246–251. [Google Scholar]
  33. Ortiz, R.; Silva, R.; Michalak, N. Application of the Response Surface Methodology to 3DEC Analysis of Open Pit Slopes. In Proceedings of the 8th South American Congress on Rock Mechanics, Buenos Aires, Argentina, 15–18 November 2015. [Google Scholar]
  34. Dongping, D.; Li, L.; Zhao, L. Limit equilibrium method (LEM) of slope stability and calculation of comprehensive factor of safety with double strength-reduction technique. J. Mount. Sci. 2017, 14, 2311–2324. [Google Scholar] [CrossRef]
  35. Kolapo, P.; Oniyide, G.O.; Said, K.O.; Lawal, A.I.; Onifade, M.; Munemo, P. An Overview of Slope Failure in Mining Operations. Mining 2022, 2, 350–384. Available online: https://www.mdpi.com/2673-6489/2/2/19 (accessed on 11 June 2022). [CrossRef]
  36. Kavvadas, M.; Roumpos, C.; Servou, A.; Paraskevis, N. Geotechnical Issues in Decommissioning Surface Lignite Mines—The Case of Amyntaion Mine in Greece. Mining 2022, 2, 278–296. [Google Scholar] [CrossRef]
  37. Crosta, G.; Agliardi, F.; Show, C.R.; Alberti, S.; Dei Cas, L. Long-term evolution and early warning strategies for complex rockslides by real-time monitoring. Landslides 2017, 14, 1615–1632. [Google Scholar] [CrossRef]
  38. Aditiana, A.; Kubotab, T.; Shinoharab, Y. Comparison of GIS-based landslide susceptibility models using frequency ratio, logistic regression, and artificial neural network in a tertiary region of Ambon, Indonesia. Geomorphology 2018, 318, 101–111. [Google Scholar] [CrossRef]
  39. Pecoraro, G.; Calvello, M.; Piciullo, L. Monitoring strategies for local landslide early warning systems. Landslides 2019, 16, 213–231. [Google Scholar] [CrossRef]
  40. Intrieri, E.; Carlà, T.; Gigli, G. Forecasting the time of failure of landslides at slope-scale: A literature review. Earth-Sci. Rev. 2019, 193, 333–349. [Google Scholar] [CrossRef]
  41. Chen, M.; Jiang, Q. An early warning system integrating time-of-failure analysis and alert procedure for slope failures. Eng. Geol. 2020, 272, 105629. [Google Scholar] [CrossRef]
  42. Mohmmed, M. A review on slope monitoring and application methods in open pit mining activities. Int. J. Sci. Technol. Res. 2021, 10, 181–186. Available online: http://www.ijstr.org/final-print/feb2021/A-Review-On-Slope-Monitoring-And-Application-Methods-In-Open-Pit-Mining-Activities.pdf (accessed on 6 June 2022).
  43. Shang, L.; Nguyen, H.; Bui, X.N.; Vu, T.H.; Costache, R.; Hanh, L.T.M. Toward state-of-the-art techniques in predicting and controlling slope stability in open-pit mines based on limit equilibrium analysis, radial basis function neural network, and brainstorm optimization. Acta Geotech. 2022, 17, 1295–1314. [Google Scholar] [CrossRef]
  44. Wang, J.; Schweizer, D.; Liu, Q.B.; Su, A.; Hu, X.; Blum, P. Three-dimensional landslide evolution model at the Yangtze River. Eng. Geol. 2021, 292, 106275. [Google Scholar] [CrossRef]
  45. Khan, M.I.; Wang, S. Slope Stability Analysis to Correlate Shear Strength with Slope Angle and Shear Stress by Considering Saturated and Unsaturated Seismic Conditions. Appl. Sci. 2021, 11, 4568. [Google Scholar] [CrossRef]
  46. Arinze, E.; Okafor, C. Finite Element Method of Stability Analysis and Stabilization of Gully Erosion Slopes—A Study of the Otampa Gu. J. Civil Environ. Eng. 2017, 7, 8. [Google Scholar]
  47. Boujaj, A.; Lahcen, B.; Ouadif, L.; Khadija, B. A methodology based on GIS for 3D slope stability analysis. Int. J. Eng. Technol. 2016, 8, 2259–2264. [Google Scholar] [CrossRef]
  48. Albataineh, N. Slope Stability Analysis Using 2D and 3D Methods. Master’s Thesis, University of Akron, Akron, OH, USA, 2006. [Google Scholar]
  49. Michalowski, R.L. Limit analysis and stabilioty charts for 3D slope failure. J. Geotech. Geoenviron. Eng. 2010, 136, 583–593. [Google Scholar] [CrossRef]
  50. Yu, G.; Xie, M.; Liang, J.; Farooq, A.; Williams, E.J. A GIS-based 3D slope stability analysis method based on the assumed normal stress on the slip surface. Sci. Rep. 2020, 10, 4442. [Google Scholar] [CrossRef]
  51. Rosenblueth, E. Two Point Estimates in Probabilities. Appl. Math. Model. 1981, 5, 329–335. [Google Scholar] [CrossRef]
  52. Luo, W. Application of Rosenblueth moment estimationmethod into probability analysis of slope stability. Ch. J. Rock Mech. Eng. 2003, 22, 232–235. [Google Scholar]
  53. Low, B.K. Practical probabilistic approach using spreadsheet. In Uncertainty in the Geologic Environment: From theory to Practise; Shackelford, C.D., Nelson, P.P., Roth, M.J.S., Eds.; ASCE: New York, NY, USA, 1996; pp. 1284–1302. [Google Scholar]
  54. Liu, X.; Tang, H.; Xiong, C.-R.; Liu, Q.B. A new method for reliability analysis of dynamic slope stability with considering energy-time distribution. Rock Soil Mech. 2015, 36, 1428–1443, 1454. [Google Scholar] [CrossRef]
  55. Liu, F.; Yang, T.; Zhou, J.; Deng, W.; Yu, Q.L.; Zhang, P.; Cheng, G. Spatial variability and time decay of rock mass mechanical parameters: A landslide study in the Dagushan open-pit mine. Rock Mech. Rock Eng. 2020, 53, 3031–3053. [Google Scholar] [CrossRef]
  56. Govora Electric-Heating Plant; Technical Documentation; Mining Division: Berbești, Vâlcea, Romania, 2021. (In Romanian)
  57. Ministry of Energy. Sector of Mineral Energy Resources—Lignite Mining Activity. Available online: http://energie.gov.ro/wp-content/uploads/2016/08/ACTIVITATEA-MINIERA-2016-2030-22072016.pdf (accessed on 1 August 2022). (In Romanian)
  58. Google Earth Pro. Available online: https://www.google.com/intl/ro/earth/versions/ (accessed on 9 August 2022).
  59. Fodor, D.; Dican, N. Open pit exploitation of coal deposits from Berbeşti mining basin. Min. Rev. 2013, 19, 2–11. Available online: https://www.upet.ro/revistaminelor/archive/rmpdf2013/nr2en.html (accessed on 8 July 2022).
  60. Ielenicz, M.; Pătru, I.G.; Clius, M. The Romanian Subcarpathians; University of Bucharest: Bucharest, Romania, 2005; pp. 71–74. (In Romanian) [Google Scholar]
  61. Chiriță, R.V. Research on Geomorphological Changes Generated by Mining Activities in the Berbeşti Mining Basin and Their Impact on the Environment. Doctoral Thesis, University of Petrosani, Petrosani, Romania, September 2019. (In Romanian). [Google Scholar]
  62. Andreescu, I.; Ţicleanu, N.; Pană, I.; Pauliuc, S.; Pelin, M.; Barus, T. Stratigraphie des dépôts pliocènes à charbons, zone est d’Oltenie (secteur Olt-Jiu). An. Univ. Bucharest Geol. 1985, 62, 87–96. [Google Scholar]
  63. Ţicleanu, N.; Roman, S.; Givulescu, R.; Biţoianu, C. Considérations paléofloristiques et pétrographiques sur la génèse des lignites de l’ouest d’Olténie. Inst. Geol. Geof. DS 1985, 4, 71–87. [Google Scholar]
  64. Ţicleanu, N. The Genetic Study of the Main Neogene Coal Deposits in Romania Based on the Characteristic Paleophytocoenoses, with Special Regard to Oltenia. Doctoral Thesis, University of Bucharest, Bucharest, Romania, 1992. (In Romanian). [Google Scholar]
  65. Apostu, I.M.; Rada, C.; Lazăr, M.; Faur, F.; Sîli, N. Investigation of the Causes and Factors Generating Land Instability in the Berbeşti Mining Basin. Min. Rev. 2022, 28, 24–36. [Google Scholar] [CrossRef]
  66. SR EN ISO 14688-1; Geotechnical Investigations and Tests. Identification and Classification of Lands. Part 1: Identification and Description. International Organization for Standardization: Geneva, Switzerland, 2017.
  67. SR EN ISO 14688-2; Geotechnical Investigations and Tests. Identification and Classification of Lands. Part 2: Principles for a Classification. International Organization for Standardization: Geneva, Switzerland, 2017.
  68. SR EN ISO 17892-12; Geotechnical Investigations and Tests. Soil Laboratory Tests. Part 12: Determination of Liquidity and Plasticity Limits. International Organization for Standardization: Geneva, Switzerland, 2018.
  69. Geologic Site (Călan, Hunedoara, Romania). Test Reports No. 014533/from 23.11.2021. Available online: https://www.geologicsite.ro/ (accessed on 12 November 2021).
  70. Technical University of Civil Engineering Bucharest. (Bucharest, Romania). Seismic Design Code P 100-1—2013. Available online: https://www.mdlpa.ro/userfiles/reglementari/Domeniul_I/I_22_P100_1_2013.pdf (accessed on 10 August 2022). (In Romanian).
Figure 1. Berbeşti Mining Basin [57,58].
Figure 1. Berbeşti Mining Basin [57,58].
Applsci 12 09607 g001
Figure 2. Lithostratigraphy of Berbeşti Mining Basin [61].
Figure 2. Lithostratigraphy of Berbeşti Mining Basin [61].
Applsci 12 09607 g002
Figure 3. Intense fragmentation and (dry) flows of the clay upper layer.
Figure 3. Intense fragmentation and (dry) flows of the clay upper layer.
Applsci 12 09607 g003
Figure 4. (a) The alternation of rock layers; (b) the discontinuity of rocks with variable thicknesses (the pictures were taken at different moments in time, in different weather conditions and from different distances and angles relative to the working front).
Figure 4. (a) The alternation of rock layers; (b) the discontinuity of rocks with variable thicknesses (the pictures were taken at different moments in time, in different weather conditions and from different distances and angles relative to the working front).
Applsci 12 09607 g004
Figure 5. Laboratory tests: (a) direct shearing tests on yellow-brown clay and sandy marl; (b) shearing tests on lignite (mandatory planes); (c) uniaxial compression and traction (Brazilian method) on lignite [65].
Figure 5. Laboratory tests: (a) direct shearing tests on yellow-brown clay and sandy marl; (b) shearing tests on lignite (mandatory planes); (c) uniaxial compression and traction (Brazilian method) on lignite [65].
Applsci 12 09607 g005
Figure 6. Evolution of Oltețu open pit and the longitudinal and cross sections [56].
Figure 6. Evolution of Oltețu open pit and the longitudinal and cross sections [56].
Applsci 12 09607 g006
Figure 7. Section L1, Fs1 = 1.01 (December 2021).
Figure 7. Section L1, Fs1 = 1.01 (December 2021).
Applsci 12 09607 g007
Figure 8. Section L2, Fs1 = 1.22 (December 2021).
Figure 8. Section L2, Fs1 = 1.22 (December 2021).
Applsci 12 09607 g008
Figure 9. Section P2, IV quarter, Fs1 = 0.80.
Figure 9. Section P2, IV quarter, Fs1 = 0.80.
Applsci 12 09607 g009
Figure 10. Section P4, IV quarter, Fs1 = 1.01.
Figure 10. Section P4, IV quarter, Fs1 = 1.01.
Applsci 12 09607 g010
Figure 11. Section L1–Stability analysis for the proposed design.
Figure 11. Section L1–Stability analysis for the proposed design.
Applsci 12 09607 g011
Figure 12. Section P1—Stability analysis of the lower excavation step.
Figure 12. Section P1—Stability analysis of the lower excavation step.
Applsci 12 09607 g012
Figure 13. Section P1—Stability analysis of the sub-steps system (excavated in yellow-brown clay).
Figure 13. Section P1—Stability analysis of the sub-steps system (excavated in yellow-brown clay).
Applsci 12 09607 g013
Table 1. Physical–mechanical properties of the rocks.
Table 1. Physical–mechanical properties of the rocks.
Type of RockLayer Thickness
[m]
Values Determined in the Laboratory
LMP
GeoLogic Lab. [69]
Values Determined from Documentations [56]Statistically
Determined Values
w
[%]
γnatsat
[daN/m3]
c
[daN/cm2]
φ
[°]
c
[daN/cm2]
φ
[°]
c
[daN/cm2]
φ
[°]
Yellow-brown clay5–1546.15
39.66
1845/1920
1782/1825
0.23–0.63
0.33
0–10
11.18
0.22150.3110
Sandy marl3–422.35
16.13
1850/1970
1911/2000
0.43
1.72
32
22.82
0.3220.50.4521
Lignite (IInd layer)2.6–5.036.251180/13000.7 */1.6 **18 */34 **0.3/1/1.230/19/381.126.5
Sandy marl0.3–1.522.35
16.13
1850/1970
1911/2000
0.43
1.72
32
22.82
0.3220.50.4521
Lignit (Ist layer)2.3–3.236.2511870.7 */1.6 **18 */34 **0.3301.126.5
Sandy marl-22.35
16.13
1850/1970
1911/2000
0.43
1.72
32
22.82
0.3220.50.4521
w—natural moisture; γnatsat—natural/saturated volumetric weight; c—cohesion; φ—internal friction angle. * Determined by shearing strength. ** Determined based on tensile and compressive strengths.
Table 2. Results of stability analyses for working fronts (sections L1, L2, P2, and P4).
Table 2. Results of stability analyses for working fronts (sections L1, L2, P2, and P4).
SectionStageHeight
h
[m]
Slope Angle
α [°]
Thickness of the Yellow-Brown Clay Layer [m]Stability FactorObservations on the Transmission Mode of the Potential Slide Surface
Fs1Fs2Fs3
L1December 20212365101.010.810.71Sliding through layers of clay and very little sandy marl
I quarter28.556130.980.790.67Sliding through layer of clay
II quarter31.370140.740.590.53Sliding through layers of clay and sandy marl
III quarter28619.51.050.840.73Sliding through layers of clay and sandy marl
L2December 2021264771.220.980.81The potential sliding surface is transmitted over the entire height of the step exits in front of the foot of the slope
I quarter2164110.910.730.71Sliding through layers of clay and sandy marl
II quarter2064101.120.900.79Sliding through layers of clay and sandy marl
III quarter1846101.541.231.08The potential sliding surface is transmitted over the entire height of the step and exits in front of the foot of the slope
P2IV quarter276913.50.800.640.56Sliding through the clay layer and very little through the sandy marl
P4IV quarter2257101.010.810.70Sliding through the clay layer, very little through the sandy marl, tangent to the layer no. II of lignite
Table 3. Results of the stability analysis for the lateral slopes.
Table 3. Results of the stability analysis for the lateral slopes.
SectionStageHeight
h
[m]
Slope Angle
α [°]
Thickness of the Yellow-Brown Clay Layer [m]Stability FactorObservations on the Transmission Mode of the Potential Slide Surface
Fs1Fs2Fs3
P5I quarter4061180.650.520.45Sliding through layers of clay and sandy marl
P1II quarter4163280.520.420.36Sliding through the clay layer and very little through the sandy marl
P0III quarter384825.50.670.530.44Sliding through layer of clay
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Lazar, M.; Faur, F.; Apostu, I.-M. Stability Conditions in Lignite Open Pits from Romania, Case Study: Oltețu Open Pit. Appl. Sci. 2022, 12, 9607. https://doi.org/10.3390/app12199607

AMA Style

Lazar M, Faur F, Apostu I-M. Stability Conditions in Lignite Open Pits from Romania, Case Study: Oltețu Open Pit. Applied Sciences. 2022; 12(19):9607. https://doi.org/10.3390/app12199607

Chicago/Turabian Style

Lazar, Maria, Florin Faur, and Izabela-Maria Apostu. 2022. "Stability Conditions in Lignite Open Pits from Romania, Case Study: Oltețu Open Pit" Applied Sciences 12, no. 19: 9607. https://doi.org/10.3390/app12199607

APA Style

Lazar, M., Faur, F., & Apostu, I. -M. (2022). Stability Conditions in Lignite Open Pits from Romania, Case Study: Oltețu Open Pit. Applied Sciences, 12(19), 9607. https://doi.org/10.3390/app12199607

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop