Owned by: China Association for Science and Technology
Sponsored by: China Coal Society
Published by: Springer Nature
About issue

The International Journal of Coal Science & Technology is a peer-reviewed open access journal. It focuses on key topics of coal scientific research and mining development, serving as a forum for scientists to present research findings and discuss challenging issues.


Coverage includes original research articles, new developments, case studies and critical reviews in all aspects of scientific and engineering research on coal, coal utilizations and coal mining. Among the broad topics receiving attention are coal geology, geochemistry, geophysics, mineralogy, and petrology; coal mining theory, technology and engineering; coal processing, utilization and conversion; coal mining environment and reclamation and related aspects.


The International Journal of Coal Science & Technology is published with China Coal Society, who also cover the publication costs so authors do not need to pay an article-processing charge.


The journal operates a single-blind peer-review system, where the reviewers are aware of the names and affiliations of the authors, but the reviewer reports provided to authors are anonymous.


  • A forum for new research findings, case studies and discussion of important challenges in coal science and mining development

  • Offers an international perspective on coal geology, coal mining, technology and engineering, coal processing, utilization and conversion, coal mining environment and reclamation and more

  • Published with the China Coal Society

Show More
Editors-in-Chief
Suping Peng, Shimin Liu
Managing Editor
Wanjie Wang
Associate Editors
Bo Hyun Kim, Dongjie Xue, Pedram Roghanchi, Wu Xiao, Zhiqiang Wu
Publishing model
Open Access. Learn about publishing OA with us
Home > Volumes and issues > Volume 10, issue 10

Stability analysis of a slope containing water-sensitive mudstone considering different rainfall conditions at an open-pit mine

Research Article

Open Access

Published: 06 October 2023

0 Accesses

International Journal of Coal Science & Technology Volume 10, article number 64, (2023)

Abstract

Mudstone, as a typical soft rock with wide distribution, has been endangering the slopes containing mudstone by its water-sensitivity of swelling and weakening strength when encountering water. To comprehensively understand the water-sensitivity of mudstone and reveal its influence on slope stability, we took the working slope containing water-sensitive mudstone of Shengli No.1 open-pit coal mine in Xilinhot, Inner Mongolia, China, as an example. Mudstone samples taken from the working slope were remodeled and saturated, and then triaxial tested to obtain the effective cohesion and effective internal friction angle. The filter paper method was used to obtain the soil–water characteristic curve of unsaturated mudstone. The pore structure of mudstone samples with different water contents were analyzed using the mercury intrusion porosimetry tests combined with the fractal dimension. The total pore content of the mudstone sample with lower water content is greater than that of the mudstone sample with higher water content. The mesopores are more in the mudstone sample with lower water content, while the small pores are more in the mudstone sample with higher water content. The variation of water content will change the complexity of mudstone pore structure. The higher the water content, the simpler the mudstone pore structure and the smoother the pore surface. Numerical calculations were conducted on the stability of the working slope under different rainfall conditions. The effective saturation on the mudstone layer surface changed and the plastic strain all occurred on the mudstone steps under different rainfall conditions. The key to preventing landslide of the slope containing water-sensitive mudstone in Shengli No.1 open-pit coal mine is to control the deformation and sliding of the mudstone layer.

1.Introduction

Slope stability has always been a concern for many scholars. Compared with other natural or artificial slopes, open-pit coal mine slopes characterized by high height, long strike, and generally unmaintained (Shu 2009; Wu et al. 2022). Therefore, once a landslide occurs, it is often characterized by a large sliding area, a large volume of the slide and a large surface drop (Li et al. 2022). Rainfall infiltration is one of the main factors causing slope instability. When surface water infiltration encounters water-sensitive soils or soft rocks, the strength of water-sensitive soils or soft rocks will be significantly reduced and accompanied by a significant increase in deformation, which makes the slope equilibrium state disrupted and then causes slope disaster accidents such as landslides (Wang et al. 2018; Ye et al. 2015; Long et al. 2019). Mudstone, as a widely distributed water-sensitive soft rock in China (He et al. 2021; Zhu et al. 2021), endangers the safe production and operation of open-pit mines in various regions(Zhao et al. 2020), and is one of the important challenges for mining engineering and geotechnical engineering research.

Current research on water-sensitive soils and soft rocks focuses on three aspects: deformation, strength and pore structure. The deformation characteristics of water-sensitive soils and soft rocks are mainly divided into collapsibility and expansibility. Among them, collapsible water-sensitive soils and soft rocks are mainly represented by wet-collapsible loess, while expansive soils and soft rocks mainly include expansive mudstone and clay. As the large deformation characteristics of water-sensitive soils and soft rocks brought great harm to engineering construction, scholars have carried out a lot of research on the deformation characteristics of water-sensitive soils and soft rocks. Shao et al. (2018) used a WG tripe triple high pressure consolidator to perform graded wetting tests on remodeled loess with different dry densities, and found that the critical collapse saturation degree of remodeled loess increased with the increasing initial dry density and decreased as the vertical pressure increased. Xie et al. (2018) conducted single oedometer-collapse tests on Malan, upper Lishi and lower Lishi loess soils, and observed that the collapse potential-pressure curves may present bimodal characteristics. Kong et al. (2018) conducted swelling deformation tests on yellow–brown and magenta mudstones from a representative landslide site in the Yanji section of the Jilin-Hunchun high-speed railway using a consolidation instrument, and found that the swelling strain of both mudstones increased with decreasing initial water content and vertical pressure. Dai et al. (2021) proposed a decoupling analysis method for the seepage and swelling of red-bed mudstone, predicting the process of water infiltration through the fissures into the mudstone and the temporal and spatial variation of swelling caused by the slow water absorption of the mudstone.

The weakening of the strength of water-sensitive soils and soft rocks when encountering water is very significant, which seriously affects the stability of tunnels, slopes, and other engineering projects. Many scholars have studied the changes in the strength of water-sensitive soils and soft rocks when encountering water. Xu et al. (2009) observed that the strength of mudstone weakened by water at different stress stages by experiments. Liu and Li (2012) carried out direct shear tests on natural and saturated schist and concluded that the cohesion, internal friction angle, residual cohesion and residual internal friction angle of saturated schist were reduced by 36.7%, 4.6%, 66.7% and 11.3%, respectively, compared to the natural state. Li et al. (2019) pointed out that the uniaxial compressive strength and elastic modulus of mudstone gradually decreased and stabilized with the increase of clay mineral content. Yao et al. (2021) found that the shear strength, cohesion and internal friction angle of mudstone samples decreased exponentially, negative exponentially and linearly, respectively, with increasing water content.

The pore structure is one of the most important factors determining the mechanical behavior and physical properties of soils and soft rocks (Enrique and Paul 2008; Zhang et al. 2023). Sun et al. (2020) conducted environmental scanning electron microscope (ESEM) and MIP tests on bentonite samples with two different initial dry densities at various suctions along wetting and drying paths, and pointed out that slight increase in fractal dimension of bentonite with increasing suction and with increasing dry density. Wang et al. (2020a, b) observed that loess collapse is mainly caused by the destruction of large ang medium pores, where the destruction of medium pores contributed more significantly to the loess collapse. Wei et al. (2020) pointed out that the swelling and dispersion of clay mineral aggregations especially the degraded illite and illite/smectite mixed layer, and the loss of high suction initiate structural failure under loading and wetting conditions. Ge et al. (2021) confirmed the significant differences in pore size distribution and structural arrangement between samples with the same vertical effective stress following different loading and saturation paths by one-dimensional (1D) compression tests, MIP tests and scanning electron microscopy (SEM) tests on compacted loess from Xi’an.

In summary, scholars have conducted a large number of studies on deformation, strength and pore structure of water-sensitive soils and soft rocks, but there are few studies on SWCC of water-sensitive mudstones. SWCC is used to describe relationship between matrix suction and volumetric water content of unsaturated soils and to predict the permeability of unsaturated soils (Kristo et al. 2019). In addition, studies on the pore structure of water-sensitive soils and soft rocks have mostly focused on research objects such as loess and swelling soils, and there are fewer studies on the changes in the pore structure of mudstone in contact with water.

On October 1, 2020, a landslide occurred in the upper 4 steps of the slope in the Shengli No. 1 open-pit coal mine due to continuous rainfall, with a sliding volume of about 220,000 m3 and a surface subsidence of about 3 m, damaging the transportation channel and affecting normal production. In this study, the working slope containing water-sensitive mudstone of Shengli No.1 open-pit coal mine was used as the engineering background. The shear strength parameters of saturated–unsaturated mudstone were obtained by triaxial tests and filter paper method. The pore structure of mudstone with different water contents was studied through MIP testing. The stability of the working slope under different rainfall conditions was studied by numerical simulation method. This study is of great significance for comprehensively understanding the water-sensitive characteristics of mudstone, revealing the impact of water-sensitive mudstone on slope stability, and for disaster prevention and control of slopes containing mudstone weak layers.

2.Materials and methods

2.1 Experimental materials

The mudstone in this study was taken from the mudstone step in the working slope of Shengli No.1 open-pit coal mine, as shown in Fig. 1. The Shengli No.1 open-pit coal mine is located in the southwest of Shengli Coalfield in Inner Mongolia Autonomous Region, China. It is 6 km south of the northern suburbs of Xilinhot City. The geographical coordinates of the mining area are 115° 56′ 51″–116° 01′ 50″ east longitude and 43° 58′ 16″–44° 03′ 10″ north latitude.

Fig. 1
figure 1

Sampling location of mudstone

The mudstone sampled is dark gray in color and is prone to disintegration and expansion when exposed to water. The mudstone was tested for water content, density, free expansion, liquid limit, plastic limit and saturated permeability coefficient. The basic physical properties of the mudstone are shown in Table 1.

Table 1 Basic physical properties of mudstone

Test items

Test results

Natural water content (%)

19.2

Natural Density (g/cm3)

1.94

Free expansion rate (%)

48.5

Liquid limit (%)

60.26

Plastic limit (%)

29.91

Plasticity index

30.35

Saturated permeability coefficient (cm/s)

7.67 × 10−7

As can be seen from Table 1, the mudstone is weakly expansive, with a large plasticity index and poor permeability. X-ray whole rock diffraction and clay separation diffraction tests were performed on the mudstone. The scanning speed was 2°/min and the sampling step width was 0.02°. The mineral composition content of the mudstone is shown in Tables 2 and 3.

Table 2 Mineral composition content of mudstone

Constituent minerals

Quartz

Potassium feldspar

Plagioclase

Clay minerals

Content (%)

43.0

3.9

5.6

47.5

Table 3 Clay mineral composition content of mudstone

Relative content of clay minerals (%)

Ratio of mixed-layer (%S)

Smectite

Illite/smectite mixed layer

Illite

Kaolinite

Chlorite

Chlorite/smectite mixed layer

Illite/smectite mixed layer

Chlorite/smectite mixed layer

73

11

16

60

According to Tables 2 and 3, the mudstone is mainly composed of clay minerals and quartz. The proportion of clay minerals reached 47.5%. The main clay minerals are illite/smectite mixed layer, kaolinite and illite. The illite/smectite mixed layer is a water-sensitive mineral that accounts for a relatively high proportion in the mudstone. Therefore, the mudstone is highly water-sensitive and is prone to disintegration and expansion when exposed to water.

2.2 Shear strength test of saturated–unsaturated mudstone

The saturation of rock and soil mass at different positions in the slope often varies due to the influence of groundwater and rainfall infiltration. For the same rock and soil mass, there are significant differences in shear strength under different saturation. The shear strength of rock and soil mass directly determines the stability of the slope, so it is necessary to study the shear strength of rock and soil mass with different saturation.

2.2.1 Triaxial tests

The gas in the pores inside the saturated mudstone is completely extruded, and the total stress consists of the effective stress borne by the mudstone skeleton and the pore water pressure borne by the pore water. According to the Mohr Coulomb strength theory and the Terzaghi principle of effective stress, the shear strength of saturated mudstone is:

$$\tau_{\text{f}} = c^{{\prime }} + \left( {\sigma - u_{\text{w}} } \right)\tan \varphi^{{\prime }}$$
(1)

where τf is shear strength of the mudstone, c′ is effective cohesion of the mudstone, σ is normal total stress of the mudstone, φ′ is effective internal friction angle of the mudstone, and uw pore water pressure in the mudstone.

The consolidated undrained (CU) triaxial tests for saturated remodeled mudstone at different confining pressures (σ3) was performed by the GDS triaxial test system (Fig. 2) to measure the effective cohesion and effective internal friction angle of the mudstone. The confining pressures (σ3) are 50 kPa, 100 kPa, 150 kPa, 200 kPa, 250 kPa and 300 kPa respectively.

Fig. 2
figure 2

GDS triaxial test system

2.2.2 Measurement of matrix suction by filter paper method

Unsaturated mudstone contains pores that are not filled with water and is usually considered to be a four-phase system of solid phase (mudstone particles), liquid phase (pore water), gas phase (pore gas), and air–water interface (shrinkage film) (Fredlund and Morgenstern 1977). Shrinkage film affects the mechanical properties of mudstone by pulling the mudstone particles together through surface tension (Kristo et al. 2019). The stress state of unsaturated mudstone is represented by two independent stress state variables, net normal stress (σ − ua) and matric suction (ua − uw) (Fredlund and Morgenstern 1977; Matyas and Radhakrishna 1968). The shear strength of unsaturated mudstone based on the Mohr–Coulomb strength theory and the Terzaghi principle of effective stress can be expressed as (Bishop and Blight 1963):

$$\tau_{\text{f}} = c^{\prime } + \left[ {\left( {\sigma - u_{\text{a}} } \right) + \chi \left( {u_{\text{a}} - u_{\text{w}} } \right)} \right]\tan \varphi^{\prime }$$
(2)

where ua is pore-air pressure in the mudstone, and χ is the parameter related to the saturation of the mudstone. When the mudstone is saturated, χ = 1, and when the mudstone is dry, χ = 0.

The shear strength of unsaturated mudstone based on the dual stress variable theory can also be expressed as (Fredlund et al. 1978):

$$\tau_{\text{f}} = c^{{\prime }} + \left( {\sigma - u_{\text{a}} } \right)\tan \varphi^{{\prime }} + \left( {u_{\text{a}} - u_{\text{w}} } \right)\tan \varphi^{b}$$
(3)

where φ′ is the internal friction angle corresponding to the net normal stress (σ − ua), and φb is the internal friction angle corresponding to the matrix suction (ua − uw).

The matrix suction of the mudstone was measured and the SWCC of the mudstone was plotted by the filter paper method (as shown in Fig. 3) to obtain the relationship between matrix suction and volumetric water content of the mudstone. The filter paper used in the test was "Shuangquan" No. 203, and the equation of the matrix suction curve (Wang et al. 2003) was as follows:

$$\lg h_{\text{m}} = \left\{ {\begin{array}{*{20}c} { - 0.0767\omega_{\text{fp}} + 5.493,} & {\omega_{\text{fp}} \le 47\% } \\ { - 0.0120\omega_{\text{fp}} + 2.470,} & {\omega_{\text{fp}} > 47\% } \\ \end{array} } \right.$$
(4)

where hm is matrix suction of mudstone and ωfp is mass water content of filter paper.

Fig. 3
figure 3

Matrix suction measurement test

2.3 MIP tests

MIP tests are based on capillary law of non-wetting liquid (e. g. mercury) to measure the pore diameter of rock and soil mass, as shown in Eq. (5) (Washburn 1921).

$$D = - \frac{{4\sigma_{{{\text{Hg}}}} \cos \theta_{{{\text{Hg}}}} }}{P}$$
(5)

where D is pore diameter of rock and soil mass, σHg is surface tension of mercury, θHg is contact angle of mercury with the rock and soil mass, and P is intrusion pressure. In this study, the σHg = 0.484 N/m and θHg = 130°.

The effect of different water contents on the pore structure of mudstone was tested by MIP. The dry density of the prepared mudstone samples is 1.4 g/cm3, and the mass water content is 30% and 48% (saturated water content), respectively. The prepared mudstone samples were dried using the freeze-drying method, and the freeze-drying time was 24 h.

To quantitatively describe the irregularity and surface roughness of the micro pore structure of rock and soil mass, the fractal theory has been introduced into the study of the microstructure of rock and soil mass. The fractal dimension value of the pore structure of the rock and soil mass should be between 2 and 3. The larger the fractal dimension value is, the rougher the pore surface is, and the more complex the pore structure is Zhang and Bing (2015). In this study, the thermodynamic model was used to quantitatively describe the complexity of the pore structure and the roughness of the pore surface in mudstones with different water contents. The thermodynamic model is based on the thermodynamic relationship of porous media during mercury intrusion testing. It is considering the balance between the surface energy increased by mercury intrusion and the work done by external forces on mercury (Zhang and Li 1995):

$$\ln \left( {W_{\text{n}} } \right) = \ln \left( {Q_{\text{n}} } \right) + {\text{C}}$$
(6)
$$W_{\text{n}} = \sum\limits_{i = 1}^{n} {\overline{p}_{i} } \Delta V$$
(7)
$$Q_{\text{n}} = r_{\text{n}}^{{2 - D_{\text{z}} }} V_{\text{n}}^{{D_{\text{z}} /3}}$$
(8)

Substituting Eqs. (7) and (8) into Eq. (6) leads to:

$$\ln \left( {W_{\text{n}} /r_{\text{n}}^{2} } \right) = D_{\text{z}} \ln \left( {V_{\text{n}}^{1/3} /r_{\text{n}} } \right) + {\text{C}}^{{\prime }}$$
(9)

where Wn is accumulated surface energy, Qn is a function of pore radius and pore volume, C and C’ are constants, i is the ith intrusion step, n is the total intrusion number, \(\overline{p}_{i}\) is the average intruded mercury pressure for ith intrusion, ΔV is the intruded mercury volume for ith intrusion, r is pore radius, and Dz is fractal dimension calculated from the thermodynamic model.

2.4 Numerical simulation

In this study, the working slope of Shengli No.1 open-pit coal mine was taken as the engineering background. The working slope has a height of 220 m and contains four layers of coal seams. The typical profile of the working slope is shown in Fig. 4. The physical and mechanical parameters of the rock strata in the slope are shown in Table 4, where the effective cohesion and effective internal friction angle of the mudstone are taken according to the results of the triaxial tests in Sect. 3.1.

Fig. 4
figure 4

Typical profile of the working slope in Shengli No.1 open-pit coal mine

Table 4 Physical and mechanical parameters of rock strata

Rock strata

Thickness (m)

Bulk modulus (GPa)

Shear modulus (GPa)

Cohesion (kPa)

Internal friction angle (°)

Saturated permeability coefficient (cm/s)

Density (kg/m3)

Sandy mudstone

25

2.22

1.02

28

26

2 × 10–6

1950

Mudstone

41

1.72

1.45

20.24

17

7.67 × 10–7

1940

Coal (5-1#)

3

0.85

0.38

33

28

1 × 10–5

1350

Mudstone

7.5

1.72

1.45

20.24

17

7.67 × 10–7

1940

Coal (5-2#)

5.5

0.85

0.38

33

28

1 × 10–5

1350

Sandy mudstone

31.23

2.22

1.02

28

26

2 × 10–6

1950

Coal (6#)

20

0.85

0.38

33

28

1 × 10–5

1350

Fine sandstone

79.77

2.7

1.6

118

39

1.2 × 10–5

2200

Coal (7#)

7

1

0.46

43

31

1 × 10–5

1400

Fine sandstone

50

2.7

1.6

118

39

1.2 × 10–5

2200

2.4.1 Numerical model and simulation schemes

In this study, COMSOL Multiphysics software was used to numerically analyze the stability of the working slope under different rainfall conditions. COMSOL Multiphysics is a software that can simulate and calculate multi-physics field. The two-dimensional numerical model of the working slope was established according to Fig. 4 and imported into COMSOL Multiphysics software. The coupling of the stress field and seepage field of the slope model was realized and numerically calculated with the help of the solid mechanics interface and Darcy's law interface in the software.

The numerical model of the working slope is 1840 m long and 270 m high. 26,738 domain units and 10,455 boundary units were divided in COMSOL Multiphysics software for meshing the numerical model of the working slope, as shown in Fig. 5.

Fig. 5
figure 5

Meshing of working slope numerical model

Xilinhot is a semi-arid grassland climate zone with cold winter and hot summer. The annual average precipitation is 294.74 mm. The rainy season is in June, July, and August of each year. The precipitation in the rainy season accounts for about 70% of the annual precipitation. The annual maximum precipitation is 511.0 mm, the monthly maximum precipitation is 281.8 mm, and the daily maximum precipitation is 86.8 mm. To study the stability changes of the working slope under different rainfall conditions, and considering the limited influence of a single rainfall on slope stability (Pan et al. 2020; Wang et al. 2020a, b; Tang et al. 2018; Cho 2016; Qian et al. 2021), we designed different rainfall simulation schemes for the combination of antecedent rainfall and major rainfall according to the classification standard of rainfall intensity and the historical rainfall conditions in the Xilinhot area. Among them, the antecedent rainfall is all with a rainfall intensity of 10 mm/d and a rainfall duration of 5 days. The time interval between the antecedent rainfall and the major rainfall is 5 days. The major rainfall simulation schemes are shown in Table 5.

Table 5 Rainfall simulation scheme

Rainfall scheme

Rainfall intensity (mm/d)

Rainfall grade

Rainfall duration (d)

Total rainfall (mm)

I

20

Medium rain

1

20

II

45

Heavy rain

1

45

III

90

Rainstorm

1

90

IV

10

Medium rain

4.5

45

V

10

Medium rain

9

90

VI

30

Heavy rain

3

90

2.4.2 Numerical model initial conditions

In this study, the safety factors of the working slope under different rainfall conditions were calculated by the strength reduction method. The shear strength parameters of the rock and soil were reduced by a custom program, and the results of the reduced calculation were substituted into the working slope numerical model for calculation through parameter scanning. The infiltration process of rainwater in the working slope was described using the Richards equation (Richards 1931):

$$\frac{\partial }{\partial x}\left[ {k_{ux} \left( h \right)\frac{\partial H}{{\partial x}}} \right] + \frac{\partial }{\partial y}\left[ {k_{uy} \left( h \right)\frac{\partial H}{{\partial y}}} \right] + \frac{\partial }{\partial z}\left[ {k_{uz} \left( h \right)\frac{\partial H}{{\partial z}}} \right] + S = \frac{\partial \theta }{{\partial t}}$$
(10)

where kux(h), kuy(h) and kuz(h) are the unsaturated permeability coefficients of the rock and soil in x, y and z, h is the pressure head, H is the total head, S is the rainfall flow, and t is the time. In this study, we assumed that the rock and soil are isotropic, that is, kux(h) = kuy(h) = kuz(h) = ku(h).

According to the SWCC of unsaturated rock and soil, if C = ∂θ/∂h, then Eq. (10) becomes:

$$\frac{\partial }{\partial x}\left[ {k_{u} \left( h \right)\frac{\partial H}{{\partial x}}} \right] + \frac{\partial }{\partial y}\left[ {k_{u} \left( h \right)\frac{\partial H}{{\partial y}}} \right] + \frac{\partial }{\partial z}\left[ {k_{u} \left( h \right)\frac{\partial H}{{\partial z}}} \right] + S = C\frac{\partial h}{{\partial t}}$$
(11)

where C is the specific water capacity and is the slope of the SWCC. The equation for the specific water capacity (C) is defined in the Van Genuchten model:

$$C = \frac{\partial \theta }{{\partial h}} = \left\{ {\begin{array}{*{20}c} {\frac{am}{{1 - m}}\left( {\theta_{s} - \theta_{r} } \right)Se^{\frac{1}{m}} \left( {1 - Se^{\frac{1}{m}} } \right)^{m} ,} & {h < 0} \\ {0,} & {h \ge 0} \\ \end{array} } \right.$$
(12)

where Se is the effective saturation, which is given by the equation:

$$Se = \frac{{\theta - \theta_{r} }}{{\theta_{s} - \theta_{r} }} = \left\{ {\begin{array}{*{20}c} {\frac{1}{{\left[ {1 + \left| {ah} \right|^{n} } \right]^{m} }},} & {h < 0} \\ {1,} & {h \ge 0} \\ \end{array} } \right.$$
(13)

In addition, the unsaturated permeability coefficients of the rock and soil (ku(h)) in Eq. (10) can be given by the equation:

$$k_{u} \left( h \right) = k_{r} k_{s}$$
(14)

where kr is relative permeability coefficient.

The Van Genuchten model defines the relative permeability coefficient (kr) as:

$$k_{r} = \left\{ {\begin{array}{*{20}c} {Se^{0.5} \left[ {1 - \left( {1 - Se^{\frac{1}{m}} } \right)^{m} } \right]^{2} } & {h < 0} \\ {1,} & {h \ge 0} \\ \end{array} } \right.$$
(15)

The a, m, and n in Eqs. (12), (13), and (15) are the fitting parameters of the Van Genuchten model, which were fitted based on the test results in Sect. 3.2.

The displacement boundary and initial water pressure boundary of the slope numerical model were set as follows:

(1) Set a fixed constraint at the undersurface of the slope numerical model (y = 0), specifying that the displacement at the undersurface of the model is 0 in all directions.

(2) Set displacement constraints on both sides of the slope numerical model (x = 0, x = 1840), specifying that the x-direction displacement on both sides of the model is 0.

(3) Set the water table in the slope numerical model as follows:

$$\left\{ {\begin{array}{*{20}l} {y = - 0.08x + 180,} \hfill & {x \le {1690}} \hfill \\ {y = 45,} \hfill & {1690 < x \le 1840} \hfill \\ \end{array} } \right.$$
(16)

(4) Set water pressure boundaries on both sides of the slope numerical model. The water pressure boundary on the left side of the model (x = 0, 0 ≤ y ≤ 180) is [9800 × (180—y)] Pa, and the water pressure boundary on the right side of the model (x = 1840, 0 ≤ y ≤ 45) is [9800 × (45- y)] Pa.

3.Results and discussion

3.1 Effective shear strength of the mudstone

The failure of the saturated remolded mudstone sample after the consolidated undrained (CU) triaxial tests is shown in Fig. 6. The saturated remolded mudstone samples under different confining pressures (σ3) all exhibited significant shear expansion failure after CU triaxial tests.

Fig. 6
figure 6

Failure of mudstone samples

The stress–strain curves of saturated remolded mudstone samples under different confining pressures are shown in Fig. 7a. The stress–strain curves of saturated remolded mudstone samples can be roughly divided into three stages:

  1. (1)

    When the axial strain ε < 1%, the stress–strain curves are approximately linear, and the saturated remodeled mudstone samples are in an elastic deformation state.

  2. (2)

    When the axial strain is 1% < ε < 3%, the stress–strain curves become nonlinear curves, and the plastic proportion of the saturated remolded mudstone samples increases. When confining pressure (σ3) is 100 kPa and 150 kPa, the rate of increase of deviatoric stress (q) with the increase axial strain ε gradually slows down, and when the deviator stress (q) increases to a peak, it begins to decrease. When the confining pressure (σ3) is 50 kPa, 200 kPa, 250 kPa and 300 kPa, the rate of increase of deviatoric stress (q) gradually decreases but no peak occurs.

  3. (3)

    When axial strain ε > 3%, the deviator stress (q) continues to slowly increase with the increase of axial strain (ε).

Fig. 7
figure 7

Triaxial test results under different confining pressures

As shown in Fig. 7b, the peak deviatoric stress (q) of the saturated remolded mudstone samples exhibits a linear relationship with the variation of confining pressure.

The variation of pore water pressure (uw) with axial strain ε for saturated remodeled mudstone samples under different surrounding pressures is shown in Fig. 7c. The pore water pressure-strain curves of saturated remolded mudstone samples can also be roughly divided into three stages:

  1. (1)

    When the axial strain ε < 1%, the pore water pressure (uw) increases linearly with the increase of axial strain (ε).

  2. (2)

    When the axial strain is 1% < ε < 3%, the increase rate of pore water pressure (uw) gradually slows down with the increase of axial strain (ε).

  3. (3)

    When axial strain ε > 3%, the pore water pressure (uw) tends to be stable with the increase of axial strain (ε).

As shown in Fig. 7d, the peak pore water pressure of the saturated remolded mudstone samples shows a linear relationship with the variation of confining pressure.

The Mohr–Coulomb failure envelope of the mudstone is plotted according to the Mohr–Coulomb strength theory and Terzaghi effective stress principle, as shown in Fig. 8. The intercept of the Mohr–Coulomb failure envelope is the effective cohesion c′ of the mudstone, and the slope is the tangent value of the effective internal friction angle φ′ of the mudstone. From Fig. 8, it can be seen that the effective cohesion (c′) is 20.24 kPa, and the effective internal friction angle (φ′) is 17°.

Fig. 8
figure 8

Mohr–Coulomb failure envelope of the mudstone

3.2 SWCC of the mudstone

In this study, Logistic curve was selected to fit the filter paper test data after comparing numerous curves in Origin software. The Logistic curve equation is as follows:

$$y = \frac{{A_{1} - A_{2} }}{{1 + \left( {x/x_{0} } \right)^{p} }} + A_{2}$$
(17)

where A1, A2, x0 and p are the fitting parameters.

The matrix suction is 0 when the mudstone is saturated. However, there is no 0 point in the logarithmic axis, therefore the matrix suction when the mudstone is saturated is taken as 1 kPa. The curve fitted to the test data based on the logistic curve equation is shown in Fig. 9a. The fitting parameters A1 = 68.86114, A2 = 16.97179, x0 = 54.70808 and p = 1.40178. The fitting degree R2 = 0.9805 indicates a good fitting effect.

Fig. 9
figure 9

SWCC of the mudstone

The SWCC can be divided into three zones as the volume water content of the mudstone decreases: boundary effect zone, transition zone and residual zone. By solving the coordinates of feature points A and B in Fig. 9a, it was found that the air entry value (AEV) of the mudstone is 13.57 kPa, the saturated volume water content (θs) is 67.2%, the residual volume water content (θr) is 20.79%, and the residual matrix suction (ψr) is 193.97 kPa.

When the matrix suction is less than 13.57 kPa, the mudstone is in a saturated or near-saturated state, and its internal pores are filled with water and the pore water is interconnected. When the matrix suction gradually increases to more than 13.57 kPa, air gradually enters the internal pores of the mudstone. As the further increase of the matrix suction, the water and air in the internal pores of the mudstone are in an unconnected state, and the matrix suction generated by the shrinkage film significantly increases with the decrease of the volume water content of the mudstone. When the matrix suction increases to 193.97 kPa, the air in the internal pores of mudstone is interconnected, and only a small portion of water exists in discontinuous small pores. At this time, the mudstone is close to a dry state, and the change in matrix suction has little effect on the volume water content of mudstone.

It is worth noting that although the Logistic curve equation fits the test data well, the left curve of the boundary effect zone of the SWCC is higher than the saturated volume water content of the mudstone, and the fitting parameters have no physical significance. Further improvement is still needed for this fitting curve.

Many empirical models for predicting SWCC by fitting fewer test data have been proposed by many scholars on the basis of a large number of experiments, due to the high cost, complex operation, or long cycle, and high workload of matrix suction measurement methods. Among them, the more classic models include the Gardner model (Gardner 1958), Van Genuchten model (Van Genuchten 1980), and Fredlund & Xing model (Fredlund and Xing 1994), as shown in Table 6. Here, a, b, m and n are the fitting parameters of each SWCC model.

Table 6 Classical SWCC models

SWCC models

Model equations

Fitting parameters

Gardner model

\(\theta = \theta_{r} + \frac{{\theta_{s} - \theta_{r} }}{{1 + a\psi^{b} }}\)

a, b

Van Genuchten model

\(\theta = \theta_{r} + \frac{{\theta_{s} - \theta_{r} }}{{\left[ {1 + \left( {a\psi } \right)^{n} } \right]^{m} }}\)

a, m, n

Fredlund & Xing model

\(\theta = C\left( \psi \right)\frac{{\theta_{s} }}{{\left\{ {\ln \left[ {e + \left( {\frac{\psi }{a}} \right)^{n} } \right]} \right\}^{m} }}\)

\(C\left( \psi \right) = 1 - \frac{{\ln \left( {1 + \frac{\psi }{{\psi_{r} }}} \right)}}{{\ln \left( {1 + \frac{{10^{6} }}{{\psi_{r} }}} \right)}}\), \(0 \le \psi \le 10^{6} \,{\text{kPa}}\)

a, m, n

The SWCCs of mudstone based on three classical SWCC models fitted to the test data are shown in Fig. 9b, and the model fitting parameters and degree are shown in Table 7.

Table 7 Fitting parameters and degree of classical SWCC models for mudstone

SWCC models

Fitting parameters

Fitting degree R2

Gardner model

a = 8.12442 × 10–4, b = 1.82221

0.97357

Van Genuchten model

a = 0.03546, m = 0.40914, n = 2.82345

0.97699

Fredlund & Xing model

a = 24.07855, m = 0.37678, n = 4.08817

0.9928

As shown in Fig. 9b and Table 7, all three typical SWCC models are well fitted, and the fitting degree reach above 0.97. Compared with the Logistic curve, all three typical SWCCs pass the saturated point of the mudstone. The Gardner and Van Genuchten models fit the boundary effect zone and transition zone of the SWCC better, but the residual zone deviates from the actual data. The Fredlund & Xing model fits the boundary effect zone, transition zone, and residual zone of the SWCC well. Therefore, the Fredlund & Xing model is more suitable for fitting and predicting the SWCC of the mudstone.

3.3 Pore size distributions of the mudstones with different water content

Figure 10 shows the pore size distribution (PSD) curves of the mudstone samples with different water contents. It can be seen from Fig. 10a that the total pore volume of the mudstone samples with 30% mass water content is larger than that of the mudstone samples with 48% mass water content. The shape of cumulative pore volume curves was similar, and the slope of the curves changes at pore diameters of approximately 0.1 μm, 1.5 μm and 25 μm. To visually analyze the pore content of mudstone samples with different water contents, the pores of mudstone are divided into micropores (< 0.1 μm), small pores (0.1–1.5 μm), mesopores (1.5–25 μm), and macropores (> 25 μm) according to their diameters.

Fig. 10
figure 10

PSDs of mudstone samples with different water contents

Figure 10b shows the pore size density functions of mudstone samples with different water contents are bimodal, in which the first peak is in the small pores and the second peak is in the mesopores. There is a significant difference in the two peak values of pore size density functions for mudstone samples with different water contents. The pore size density functions of the mudstone samples with different water contents in the micropores and macropores do not differ significantly, indicating that changes in water content have a relatively small impact on micropores and macropores. The most probable pore diameter indicates that the content of pores with that diameter is the highest in mudstone samples. The most probable pore diameter of the mudstone sample with 30% mass water content is 19.10 μm, while the most probable pore diameter of the mudstone sample with 48% mass water content is 6.75 μm. This may be because when the water content increases, the clay minerals in the mudstone expand with water and occupy the space of the mesopores, leading to the shrinkage of the mesopores into small pores.

The fractal dimensions of mudstone samples with different water contents are shown in Fig. 11. According to Eq. (9), the slope of the fitting line in Fig. 11 is the fractal dimension (Dz) of the thermodynamic model. The fractal dimension of mudstone sample with 30% mass water content is 2.87955, while the fractal dimension of mudstone sample with 48% mass water content is 2.84311. This indicates that the change of water content changes the complexity of the pore structure and the roughness of the pore surface inside the mudstone samples. The higher the water content, the simpler the pore structure and the smoother the pore surface inside the mudstone, which causes changes in physical and mechanical properties such as increased permeability and decreased strength of the mudstone.

Fig. 11
figure 11

Fractal dimension of mudstone samples with different water contents

3.4 Stability of working slope under different rainfall conditions

3.4.1 Changes in slope saturation

According to the initial conditions of the slope numerical model in Sect. 2.4.2, the rock strata below the water table are saturated, and the initial effective saturation of the rock strata above the water table is set to 30%, as shown in Fig. 12a.

Fig. 12
figure 12

Effective saturation of rock strata in the working slope under different rainfall conditions

The effective saturation of the rock strata in the working slope under different rainfall conditions is shown in Fig. 12b–i. Due to the overall low permeability of the rock strata in the working slope, there is no hydraulic connection between rainwater and groundwater after rainfall, and only the effective saturation of the surface rock strata in the slope changed. The focus of this study is on the water-sensitive mudstone in the working slope, so the mudstone layer was locally amplified. It can be seen from the local enlargement of the mudstone layer that the effective saturation of the mudstone layer varies under different rainfall conditions. The maximum effective saturation of the mudstone layer surface is 54.80% after the antecedent rainfall. The saturation of the mudstone layer surface decreased slightly after 5 days of the antecedent rainfall, i.e., before the major rainfall, with a maximum effective saturation of 54.60%. When the rainfall intensity is 20 mm/d, the maximum effective saturation of the mudstone layer surface increases to 70.44% after 24 h of rainfall. When the rainfall intensity is 45 mm/d, some areas of the mudstone layer surface are close to saturation after 24 h of rainfall, and the maximum effective saturation increases to 98.24%. When the rainfall intensity is 90 mm/d, after 24 h of rainfall, the mudstone layer surface is fully saturated and the mudstone layer shows obvious saturated and unsaturated transition zones. These results indicate that the greater the rainfall intensity, the greater the effective saturation of the mudstone layer surface when the rainfall duration is the same.

Comparing the simulation results of Scheme II and Scheme IV both with a total rainfall of 45 mm, it can be seen that when the rainfall intensity is 10 mm/d, after 4.5 days of rainfall, the maximum effective saturation of the mudstone layer surface is 93.31%, which is smaller than that in Scheme II. This is due to the fact that the longer the duration of the rainfall, the larger the rainfall infiltration amount when the rainfall total is the same, resulting in a slight decrease in the mudstone layer surface saturation. The mudstone layer surface in Scheme III, V and VI, in which the total amount of rainfall is 90 mm, is fully saturated, and all of them show obvious saturated and unsaturated transition zones. In addition, the simulation results of Scenarios IV and V, in which the rainfall intensity is 10 mm/d, show that when the rainfall intensity is the same, the longer the rainfall duration is, the greater the effective saturation of the mudstone layer surface is.

3.4.2 Changes in slope stability

The safety factors and the critical slip surfaces position of the working slope under different rainfall conditions were obtained based on the calculated non convergence or plastic zone connection as the slope instability criterion. The equivalent plastic strains of the working slope under different rainfall conditions at the time of reduction failure are shown in Fig. 13. The safety factors of the working slope under different rainfall conditions are shown in Fig. 14.

Fig. 13
figure 13

Equivalent plastic strain of the working slope under different rainfall conditions

Fig. 14
figure 14

Safety factors of the working slope under different rainfall conditions

Figure 13 shows that the plastic zones of the working slope at the time of reduction failure under different rainfall conditions appears on the shallow layer of mudstone steps, indicating that the critical slip surfaces of the working slope under different rainfall conditions are all located on the shallow layer of mudstone steps, and the mudstone steps are prone to shallow sliding. Due to the large size of the working slope, only the mudstone weak layer in the upper part of the slope generates plastic strain under rainfall, without forming a plastic zone that runs through the entire slope. Therefore, the influence of the rainfall scheme simulated in this study on the safety factors of the slope is small. However, it can still be used to explore the changes in the stability of a slope containing mudstone weak layers under different rainfall conditions.

As shown in Fig. 14, the safety factors of the working slope decrease under different rainfall conditions. The safety factor of the working slope decreased from 1.180 to 1.170 at the end of the antecedent rainfall. The plastic strain value at the mudstone step 5 days after the end of the antecedent rainfall, i.e., before the major rainfall, is less than the plastic strain value at the mudstone step at the end of the antecedent rainfall, and the safety factor of the working slope increases to 1.171. Comparing the results of each simulation scheme of major rainfall, it can be seen that when the rainfall duration is the same, the greater the rainfall intensity is, the greater the plastic strain in the shallow layer of the mudstone step is, and the smaller the safety factor of the working slope is. When the rainfall intensity is the same, the longer the rainfall duration, the greater the plastic strain in the shallow layer of the mudstone step, and the smaller the safety factor of the working slope. These indicate that the total amount of rainfall is the key factor affecting the stability of the slope. When the total amount of rainfall is the same, the longer the rainfall duration, the smaller the plastic strain in the shallow layer of the mudstone step, and the greater the safety factor of the working slope. It can be seen that the short duration strong rainfall is more threatening to the stability of the slope than the long duration weak rainfall.

The permeability of the rock strata on the working slope of Shengli No.1 open-pit coal mine is poor, and only the water content on the surface of the slope increases under rainfall infiltration. The water-sensitivity of mudstone leads to a significant weakening of its strength and a significant increase in deformation when encountering water, while other rock strata in the slope are non-water-sensitive (in fact, sandy mudstone in the working slope is also weakly water-sensitive. and the weak water-sensitivity of the sandy mudstone is not considered in this study in order to highlight the influence of the weak layer of water-sensitive mudstone on the stability of the slope). After encountering water, the strength will not change significantly. Therefore, the plastic strain occurs only in the mudstone steps after rainfall.

In addition, the fact that rainwater is not contacted with groundwater also results in the absence of a plastic zone that runs through the entire slope. Of course, if the rainfall is large enough and the rainwater keeps infiltrating downward, causing the increase of water content and decrease of shear strength of non-water-sensitive rock strata in a large area of the working slope, as well as the rise of groundwater level, the working slope may form a plastic zone penetrating the slope as a whole and large-scale sliding instability will occur. However, the Xilinhot area is a semi-arid grassland climate with low rainfall in the rainy season and throughout the year, making it unsuitable for large-scale rainfall. Therefore, the key to preventing landslide of the slope containing mudstone weak layer in Shengli No.1 open-pit coal mine is to control the deformation and sliding of the mudstone weak layer.

4.Conclusions

With the working slope containing water-sensitive mudstone of Shengli No.1 open-pit coal mine, as an example, this study investigated the water sensitive characteristics of mudstone and analyzed the stability of the working slope. Firstly, the shear strength of saturated–unsaturated mudstone was investigated using the GDS triaxial test system and the filter paper method. Secondly, the pore structure of mudstone samples with different water contents was analyzed by the MIP tests combined with the fractal dimension. Finally, the stability of the working slope under different rainfall conditions was numerically calculated. The following conclusions can be drawn from the analysis and discussion of the experiments and numerical calculations:

  1. (1)

    The saturated remolded mudstone samples under different confining pressures all exhibited significant shear expansion failure after CU triaxial tests. The curves of deviatoric stress and pore water pressure with the increase of axial strain can be divided into three stages: rapid linear growth, slow growth and stabilization. The peak values of deviatoric stress and pore water pressure are linearly positively correlated with confining pressure. The effective cohesion of the mudstone is 20.24 kPa, and the effective internal friction angle is 17°.

  2. (2)

    The SWCC of the mudstone fitted by the Logistic equation shows that the AEV of the mudstone is 13.57 kPa, the saturated volume water content is 67.2%, the residual volume water content is 20.79%, and the residual matrix suction is 193.97 kPa. The Fredlund & Xing model is more suitable for fitting and predicting the SWCC of the mudstone among the three classic SWCC models.

  3. (3)

    The total pore content of the mudstone sample with lower water content is greater than that of the mudstone sample with higher water content. The mesopores are more in the mudstone sample with lower water content, while the small pores are more in the mudstone sample with higher water content. This may be because when the water content increases, the clay minerals in the mudstone expand with water and occupy the space of the mesopores, leading to the shrinkage of the mesopores into small pores.

  4. (4)

    The fractal dimensions of mudstone samples with different water contents show that the higher the water content, the simpler the pore structure and the smoother the pore surface inside the mudstone, which causes changes in physical and mechanical properties such as increased permeability and decreased strength of the mudstone.

  5. (5)

    Under different rainfall conditions, the effective saturation on surface of the working slope increases, the safety factor decreases, and the plastic zone appears on the mudstone steps. The key to preventing landslide of the slope containing mudstone weak layer in Shengli No.1 open-pit coal mine is to control the deformation and sliding of the mudstone weak layer.

References

[1] Bishop AW, Blight GE (1963) Some aspects of effective stress insaturated and partly saturated soils. Geotechnique 13(3):177–197
[2] Cho SE (2016) Stability analysis of unsaturated soil slopes considering water-air flow caused by rainfall infiltration. Eng Geol 211:184–197
[3] Dai ZJ, Guo JH, Yu F, Zhou Z, Li J, Chen SX (2021) Long-term uplift of high-speed railway subgrade caused by swelling effect of red-bed mudstone: case study in southwest china. Bull Eng Geol Environ 80(1)
[4] Enrique R, Paul HS (2008) Microstructure investigation in unsaturated soils: a review with special attention to contribution of mercury intrusion porosimetry and environmental scanning electron microscopy. Geotech Geol Eng 26:705–727
[5] Fredlund DG, Morgenstern NR (1977) Stress state variables for unsaturated soils. J Geotech Eng Div 103(5):447–466
[6] Fredlund DG, Morgenstern NR, Widger RA (1978) The shear strength of unsaturated soils. Can Geotech J 15(3):313–321
[7] Fredlund DG, Xing A (1994) Erratum : equations for the soil-water characteristic curve. Can Geotech J 31(6):1023–1025
[8] Gardner WR (1958) Some steady-state solutions of the unsaturated moisture flow equation with application to evaporation from a water table. Soil Sci 85(4):228–232
[9] Ge MM, Pineda JA, Sheng DC, Burton GJ, Ning L (2021) Microstructural effects on the wetting-induced collapse in compacted loess. Comput Geotech 138 (Article ID): 104359.
[10] He ZM, Zhong W, Liu ZF, Li Q, Shu QH (2021) Stability analysis of carbonaceous mudstone coarse-grained soil embankment slope based on improved green-ampt infiltration model. J Cent South Univ (sci Technol) 52(7):2179–2187
[11] Kong LW, Zeng ZX, Bai W, Wang M (2018) Engineering geological properties of weathered swelling mudstones and their effects on the landslides occurrence in the yanji section of the jilin-hunchun high-speed railway. Bull Eng Geol Environ
[12] Kristo C, Rahardjo H, Satyanaga A (2019) Effect of hysteresis on the stability of residual soil slope. Int Soil Water Conserv Res 7(3):226–238
[13] Li G, Yang X, Wang D, Wang Y, Yu X (2022) Stabi1ity of inner dump slope under coal pillar support case study in an open-pit coal mine. Int J Coal Sci Technol 9:25
[14] Li GC, Sun CL, He JT, Sun YT, Dong YX, Zhao HS (2019) Macro and meso scalesimulation study of the strength-weakening property of soft mudstone affected by water. J China Univ Min Technol 48(05):935–942
[15] Liu J, Li JP (2012) Study of shear strength water sensitivity of intact schist from guzhu expressway. Rock Soil Mech 33(06):1719–1723
[16] Long AF, Chen KS, Ji YX (2019) Experimental study on wetting-drying cycles of red clay slopes under different rainfall intensities. Chin J Geotech Eng 41(z2):193–196
[17] Matyas EL, Radhakrishna HS (1968) Volume change characteristics of partially saturated soils. Géotechnique 18(4):432–448
[18] Pan Y, Wu G, Zhao Z, He L (2020) Analysis of rock slope stability under rainfall conditions considering the water-induced weakening of rock. Comput Geotech 128 (Article ID): 103806
[19] Qian J, Sun D, Li G, Wu Y (2021) Study on the influence of various rainfall types on the stability of high and steep slopes. Geofluids 2021 (Article ID): 7906573
[20] Richards LA (1931) Capillary conduction of liquids through porous mediums. Physics 1(5):318–333
[21] Shao XX, Zhang HY, Tan Y (2018) Collapse behavior and microstructural alteration of remolded loess under graded wetting tests. Eng Geol 233:11–22
[22] Shu JS (2009) Study on key factors influencing slope stability and integration of mining and slope control in surface mines. Dissertation, China University of Mining & Technology
[23] Sun HQ, Mašín D, Najser J, Neděla V, Navrátilová E (2020) Fractal characteristics of pore structure of compacted bentonite studied by esem and mip methods. Acta Geotech 15(6):1655–1671
[24] Tang G, Huang J, Sheng D, Sloan SW (2018) Stability analysis of unsaturated soil slope under random rainfall patterns. Eng Geol 245:322–332
[25] Van Genuchten TM (1980) A closed-form equation for predicting the hydraulic conductivity of unsaturated soils. Soil Sci Soc Am J 44(5):892–898
[26] Wang JS, Wang SM, Hong ML, Li YZ, Wu G, Li L, Zhang JL (2018) Correlation analysis between clay mineral composition and shear strength. J Southwest Jiaotong Univ 53(5):1033–1038
[27] Wang L, Li XA, Li LC, Hong B, Yao W, Lei HN, Zhang C (2020a) Characterization of the collapsible mechanisms of malan loess on the chinese loess plateau and their effects on eroded loess landforms. Hum Ecol Risk Assess 26(9):2541–2566
[28] Wang Y, Chai J, Cao J, Qin Y, Zhang X (2020b) Effects of seepage on a three-layered slope and its stability analysis under rainfall conditions. Nat Hazards 102(3):1269–1278
[29] Wang Z, Yang JX, Juan KJ, An JY, Dao LY (2003) Application of filter paper method in field measurement of matric suction. Chin J Geotech Eng 25(4):405–408
[30] Washburn EW (1921) Note on a method of determining the distribution of pore sizes in a porous material. Proc Natl Acad Sci USA 7(4):115–116
[31] Wei YN, Fan W, Yu B, Deng LS, Wei TT (2020) Characterization and evolution of three-dimensional microstructure of malan loess. Catena 192 (Article ID): 104585
[32] Wu H, Chen Y, Lv H, Xie Q, Chen Y, Gu J (2022) Stability analysis of rib pillars in highwall mining under dynamic and static loads in open-pit coal mine. Int J Coal Sci Technol 9:38
[33] Xie WL, Li P, Zhang MS, Cheng TE, Wang Y (2018) Collapse behavior and microstructural evolution of loess soils from the loess plateau of china. J Mountain Sci
[34] Xu XL, Zhang N, Li YS (2009) Experimental study of permeability and strength weakening of mudstone considering effect of water under typical stress stages. Chin J Rock Mech Eng 28(S1):3089–3094
[35] Yao QL, Wang WN, Yang SY, Fang J, Zhou BJ (2021) Direct shear and acoustic emission characteristics of sandy mudstone under the effect of moisture content. J China Coal Soc 46(09):2910–2922
[36] Ye WJ, Zhao ZP, Yang GS, Xi JM, Zhang YY (2015) Influence of soil moisture on loess slope spalling hazards. China J Highway Transp 28(07):18–24
[37] Zhang BQ, Li SF (1995) Determination of the surface fractal dimension for porous media by mercury porosimetry. Ind Eng Chem Res 34(4):1383–1386
[38] Zhang C, Bai Q, Han P, Wang L, Wang X, Wang F (2023) Strength weakening and its micromechanism in water–rock interaction, a short review in laboratory tests. Int J Coal Sci Technol 10:10
[39] Zhang Y, Bing H (2015) Experimental study of the effect of freezing-thawing cycles on porosity characters of silty clay by using mercury intrusion porosimetry. J Glaciol Geocryol 37(1):169–174
[40] Zhao H, Wang D, Ma M, Zheng K (2020) Parameter inversion and location determination of evolutionary weak layer for open-pit mine slope. Int J Coal Sci Technol 7(4):714–724
[41] Zhu JH, Yu RG, Han SX, Tong YM, Zhang HY, Zhen ZL (2021) Strength deterioration of mudstone with different initial dry densities under dry-wet cycles. J China Railway Soc 43(10):109–117

About this article

Cite this article

Yang, G., Chen, Y., Liu, X. et al. Stability analysis of a slope containing water-sensitive mudstone considering different rainfall conditions at an open-pit mine.Int J Coal Sci Technol 10, 64 (2023).
  • Received

    31 May 2023

  • Revised

    28 July 2023

  • Accepted

    03 August 2023

  • DOI

    https://doi.org/10.1007/s40789-023-00619-z

  • Share this article

    Copy to clipboard

For Authors

Explore