Next Article in Journal
Financial Applications on Fractional Lévy Stochastic Processes
Next Article in Special Issue
Investigation on Pore Structure and Permeability of Concrete–Rock Interfacial Transition Zones Based on Fractal Theory
Previous Article in Journal
Ensemble FARIMA Prediction with Stable Infinite Variance Innovations for Supermarket Energy Consumption
Previous Article in Special Issue
Macro- and Micromechanical Assessment of the Influence of Non-Plastic Fines and Stress Anisotropy on the Dynamic Shear Modulus of Binary Mixtures
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Seepage–Fractal Model of Embankment Soil and Its Application

1
School of Civil Engineering, Xuchang University, Xuchang 461000, China
2
College of Civil and Transportation Engineering, Hohai University, Nanjing 210098, China
3
School of Resources and Geosciences, China University of Mining and Technology, Xuzhou 221116, China
4
Yellow River Institute of Hydraulic Research, No.45 Shunhe Road, Zhengzhou 450003, China
*
Author to whom correspondence should be addressed.
Fractal Fract. 2022, 6(5), 277; https://doi.org/10.3390/fractalfract6050277
Submission received: 30 March 2022 / Revised: 17 May 2022 / Accepted: 20 May 2022 / Published: 22 May 2022
(This article belongs to the Special Issue Fractal and Fractional in Geomaterials)

Abstract

:
Over time and across space, the hydraulic conductivity, fractal dimension, and porosity of embankment soil have strong randomness, which makes analyzing seepage fields difficult, affecting embankment risk analysis and early disaster warning. This strong randomness limits the application of fractal theory in embankment engineering and sometimes keeps it in the laboratory stage. Based on the capillary model of porous soil, an analytical formula of the fractal relationship between hydraulic conductivity and fractal dimension is derived herein. It is proposed that the influencing factors of hydraulic conductivity of embankment soil mainly include the capillary aperture, fractal dimension, and fluid viscosity coefficient. Based on random field theory and combined with the embankment parameters of Shijiu Lake, hydraulic conductivity is discretized, and then the soil fractal dimension is approximately solved to reveal the internal relationship between hydraulic gradient, fractal dimension, and hydraulic conductivity. The results show that an increased fractal dimension will reduce the connectivity of soil pores in a single direction, increase the hydraulic gradient, and reduce the hydraulic conductivity. A decreased fractal dimension will lead to consistency of seepage channels in the soil, increased hydraulic conductivity, and decreased hydraulic gradient.

1. Introduction

The development of fractal theory provides a new theoretical basis for the study of soil with a complex internal structure [1,2]. Fractal theory is widely used in research on geotechnical materials as a theoretical method to describe the geometric characteristics of porous media [3,4]. There are extensive experimental data and theoretical research in this field showing that the pore distribution in soil has obvious fractal characteristics [5,6]. At present, the research mainly focuses on the pore–solid model, while there is less research on the seepage–fractal model.
Soil is a loose, porous three-phase medium composed of soil particles, air, and water. Soil particles form the skeleton of soil, while air and water exist in the interconnected pores between soil particles in the soil. In nature, the hydraulic parameters of soil have strong spatial variability at different scales. The pore structure and its connectivity are very complex, and pore scale ranges from millimeters to microns. The process of fluid migration depends on the multiscale pore structure of soil [7,8]. Wang [9] measured the pore size of soil with high-resolution computed tomography (CT) imaging technology, calculated the soil porosity, reconstructed the soil pores in three-dimensional space, and studied the spatial distribution law of soil pores based on their fractal characteristics. Rainfall will introduce changes to the pore structure in the soil surface, and soil particles’ transportation caused by infiltration will clog soil pores and reduce the porosity of surface soil, causing a continuous change in soil pore micro characteristics consequently affecting the mechanical properties of soil [10]. The pore structure of soil is affected by many factors, and soil conditioner can significantly increase the number of macro and small pores to improve the soil structure and porosity [11].
Seepage refers to the phenomenon of pore water in the soil flowing through the interconnected pores in the soil particle skeleton under the action of water [12]. Soil permeability represents the ability to allow water or other fluids to flow through its interior. The study of soil permeability is very important in many engineering fields, including studies on the seepage instability of embankments, waterproof foundation pits, etc. [13,14,15]. There are interconnected pores in the soil particle skeleton, which explains the soil permeability. Therefore, the shape, size, and connection mode of pores have an important impact on the permeability of the soil. Many scholars around the world have done meaningful research on the characteristics of soil permeability; they have explored the internal connection between hydraulic conductivity and other physical parameters and tried to establish a model to determine the hydraulic conductivity [16,17].
Soil particle size distribution refers to the percentage of particles of different sizes in the total particles. It is one of the basic parameters of soil and has a strong impact on the hydraulic properties. Graded entropy can be used to represent the particle size distribution of soil. Through soil particle size, the relationship between hydraulic conductivity and graded entropy can be established. This approach can be applied to predict the hydraulic conductivity of various soils [18].
Darcy’s law describes the linear function between the seepage velocity of water in saturated soil and the hydraulic gradient, also known as the linear seepage law. However, according to test results, this linear function is only applicable under certain conditions. With increased soil viscosity, the linear relationship no longer exists. Therefore, a new method is needed to calculate the hydraulic conductivity of cohesive soil. The plasticity index, average pore diameter, and particle size distribution of cohesive soil affect its permeability, resulting in changes to the hydraulic gradient and seepage velocity. The function describing the relationship between these influencing factors and hydraulic conductivity can be obtained with the single factor test. Scholars have studied the internal relationship between the physical indices and hydraulic conductivity of cohesive soil and have tried to find the optimal method to calculate the hydraulic conductivity [19].
Stress in the soil causes consolidation and deformation, which leads to changes in the pore structure and affects the distribution of hydraulic conductivity. This phenomenon is more obvious in the construction of pile foundations [20]. The pore fractal dimension refers to the size and development degree of pore channels, indicating the distribution characteristics of soil pores [21]. In the process of seepage failure in porous media, movable fine particles are continuously lost from the pore channels, and physical and mechanical parameters such as hydraulic conductivity, porosity, and the nonuniformity coefficient change, resulting in changes to the fractal dimension [22,23,24,25]. Studying the fractal characteristics of porous soil seepage and predicting and controlling the formation and development of seepage failure are significance for the long-term and safe operation of dam engineering. Capillary curvature, pressure, and wall roughness are closely related to soil hydraulic conductivity. When the capillary structure changes, the fractal dimension changes. Therefore, hydraulic conductivity can be expressed as an equation of capillary state and the fractal dimension [26]. Computer tomography and three-dimensional modeling technology can eliminate the influence of geometric size, establish relationships between soil parameters and the fractal dimension, reflect the fractal characteristics of the model in real time and reduce the error [27]. Meanwhile, the shrinkage crack caused by water loss changes the fractal characteristics, and permeability is obviously affected [28]. The expansive agent reduces the cracks of porous media, and leads to the decrease in permeability [29,30]. Soil properties can be explained by fractal theory; however, under the action of many factors, soil physical indices are not ideally fractal [31]. For complex models, the fractional Adams–Bashforth method can quickly obtain the approximate solution, which plays a positive role in the application of fractal theory [32].
At present, the study of seepage by fractal theory is less related to the permeability of embankment soil. Research on hydraulic conductivity, hydraulic gradient, and its influencing factors is still insufficient, and the functional relationship between hydraulic conductivity, hydraulic gradient, and fractal dimension needs to be studied further. Based on the fractal characteristics of embankment soil, an analytical formula of the relationship between hydraulic conductivity, porosity, and fractal dimension is derived. Considering the influence of randomness, the seepage fractal model of embankment soil was established. Combined with engineering practice, the relationship between hydraulic conductivity, hydraulic gradient, and fractal dimension was further studied, and the internal relationship between seepage failure and fractal dimension of porous media soil is revealed.

2. Methodology

2.1. Hydraulic Conductivity Based on Pore Radius and Fractal Dimension

Soil is a typical porous medium. Spaces without a solid skeleton are called pores, which are occupied by liquid or gas, and are connected to each other to form channels. Water can flow through these pores in the soil. Therefore, the structure and scale of pores affect the hydraulic characteristics of the soil.
When there is seepage failure of soil, such as with piping, movable fine particles are continuously lost from the pore channels, resulting in changes in pore size and internal structure, and then in the fractal dimension of the soil. Under the action of seepage, some soil particles are brought to other positions by the seepage water, and the pore radius and permeability of soil change. It can be seen that pore size has an effect on the fractal dimension and hydraulic conductivity, resulting in a function that describes the relationship between them.
Hydraulic conductivity reflects the permeability of porous medium. Due to the irregularity and complexity of soil structure, a unified analytical mathematical formula cannot be used to express hydraulic conductivity at present. In engineering, some empirical values or formulas are usually used, or hydraulic conductivity is determined according to indoor geotechnical tests, in situ tests, and numerical test inversion. Considering that soil particles and pores show fractal characteristics from the atomic scale to grain size, the micro pore structure of soil is analyzed, and the analytical relationship between fractal dimension and hydraulic conductivity and porosity is deduced to further explore seepage and fractal characteristics of porous soil media.
In order to describe the characteristic of soil, a cluster of capillaries with different radiuses is introduced to simulate the pores in soil. Take a cylinder in the soil with a radius of R and a length of L. Suppose that in the cylinder, interconnected pores can be replaced by capillaries, and the radius (r) values of capillaries differ from each other. When the capillary radius satisfies r min < r < r max , the total number of capillaries with radius greater than r in the unit section of soil approximately obeys Equation (1) [33]:
M ( r c > r ) = R D r D
where D is the fractal dimension of capillaries in porous media. Since capillaries are composed of interconnected pores, their fractal dimension is equivalent to the fractal dimension of pores in porous media. R and r are the radius of the cylinder and capillary. The fractal dimension of capillaries in soil can be obtained by Equation (2):
D = d ln n ln ( r min / r max )
where d is the fractal dimension of Euclidean space, and n is the porosity of porous media. Obviously, the fractal dimension of a plane in Euclidean space is 2, and the fractal dimension of capillaries in soil meets 1 < D < 2 . The fractal dimension of three-dimensional Euclidean space is 3, and the fractal dimension of capillaries in soil meets 2 < D < 3 .
By differentiating Equation (1), when the capillary radius changes from r to r + dr, the cumulative number of capillaries can be obtained. As shown in Equation (3), it is obvious that the cumulative number of capillaries in the interval decreases with increased pore size:
d M ( r ) = D R D r 1 D d r
According to Poiseuille’s law, when fluid moves in a laminar flow in a horizontal circular pipe, such as capillary, the difference between the volume flow and the pressure at both ends of the pipe, the radius and length of the pipe, and the viscosity coefficient of the fluid meets Equation (4):
q ( r ) = 1 8 η L π r 4 ( p 1 p 2 ) = 1 8 η L π r 4 γ Δ h
where γ and η are the gravity and viscosity coefficient of fluid, respectively; r and L are the radius and length of the pipe, respectively, and Δ h is the head difference.
Assuming that the minimum and maximum capillary apertures in a section are r min and r max , respectively, and the aperture is continuously distributed, the total flow through the section can be obtained by integrating the flow through each capillary in the section in the interval [ r min , r max ] . The total flow is shown in Equation (5):
Q = r min r max q ( r ) d M ( r )
Bringing Equations (3) and (4) into Equation (5), we obtain Equation (6):
Q = π R D D 8 η γ 4 D A J ( r max 4 D r min 4 D )
Based on Darcy’s law, the hydraulic conductivity (k) of soil is shown in Equation (7):
k = v J = Q A J
where ν is the seepage velocity, and J is the hydraulic gradient. Bringing Equation (6) into Equation (7), after sorting, hydraulic conductivity meets Equation (8):
k = π R D D 8 η γ 4 D r max 4 D [ 1 ( r min r max ) 4 D ]
Soil porosity is the ratio of the pores volume to the total volume. In this study, the pores volume in soil can be expressed by the volume of capillary clusters in the interval [ r min , r max ] , and the total soil volume is equal to the volume of a hypothetical cylinder. Therefore, the soil porosity can be expressed by Equation (9):
n = V p V = r min r max π r 2 L d M ( r ) A L
where n is the porosity, V is the total soil volume, and V p is the pore volume in the capillary cluster. By introducing Equation (3) into Equation (9), we obtain Equation (10):
n = π R D D 2 D r max 2 D [ 1 ( r min r max ) 2 D ]
In nature, soil is a typical porous medium. In Equations (8) and (10), the hydraulic conductivity and porosity of soil are functions of the capillary cluster radius and fractal dimension. If the fractal dimension and pore size are different, the hydraulic conductivity and porosity are also different. Strictly speaking, due to the assumption, the relationship between hydraulic conductivity and porosity is difficult to express analytically. Considering the complexity and irregularity of soil pores, the maximum pore radius is much larger than the minimum pore radius. At this time, the relationship between maximum and minimum pore radius can be expressed by Equation (11):
r min r max 0
By introducing Equation (11) into Equations (8) and (10), the hydraulic conductivity and porosity of soil can be expressed by Equations (12) and (13):
k = π R D D 8 η γ 4 D r max 4 D
n = π R D D 2 D r max 2 D
Obviously, the hydraulic conductivity of soil is affected by the fractal dimension, pore size, water viscosity coefficient, and other factors.

2.2. Randomness of Hydraulic Conductivity of Embankment Soil

Due to the differences in internal structure and pore distribution, soil hydraulic conductivity has a certain spatial variability, resulting in different values throughout the soil. In embankment engineering, due to the influence of time and space, the hydraulic conductivity of soil shows strong randomness, which cannot be ignored. It has a significant impact on the safe operation and early risk warning of embankments. Deterministic research cannot meet the requirements of the development of modern science and technology. More engineering applications need to consider the randomness of hydraulic conductivity in order to meet the needs of design and risk protection.
The randomness of hydraulic conductivity is affected by many factors, such as pore structure, plasticity index, and compactness. The relationship between these factors is complex and difficult to study at the same time. Some scholars try to use statistical methods to study the random characteristics of hydraulic conductivity, testing multiple groups of soil samples, counting the mean and variance, and obtaining an approximate expression of the probability density function through regression analysis. The concept of randomness is introduced in these studies. Soil hydraulic conductivity is regarded as a random variable that follows a certain distribution law. It is generally believed that the log normal distribution has good adaptability and can reflect the spatial structure of soil hydraulic conductivity better, and it has been applied in some studies [34,35,36,37].
In embankment engineering, when measuring the hydraulic conductivity of multiple groups of soil samples and counting the mean and variance, the mean and variance of logarithmic hydraulic conductivity can be expressed by Equation (14):
σ ln k 2 = ln ( 1 + σ k 2 μ ¯ k 2 ) μ ¯ ln k = 1 2 σ ln k 2 + ln ( μ ¯ k )
In the discrete random field, the hydraulic conductivity of the jth soil element can be expressed by Equation (15):
ln k j = μ ¯ ln k + σ ln k G j
where G j is the local average of the standard normal random field in the jth soil element. There is a certain correlation between the points in the random field, which can be described by the correlation function. The exponential correlation function can represent the random characteristics of embankment soil better as shown in Equation (16):
ρ ( τ ) = exp ( 2 θ | τ | )
The fluctuation scale ( θ ) can be regarded as the maximum distance beyond which points in the random field are almost irrelevant. Linear distance ( τ ) between any points in the random field is closely related to correlation. The randomness of hydraulic conductivity is affected by the coefficient of variation (COV) and θ . With increased COV, randomness increases. With increased θ , the correlation between points is enhanced, randomness is reduced, and the soil tends to be uniform.
As shown in Figure 1, a hydraulic conductivity random field is generated based on the random field model in this paper. The shape of the random field is a 40 m square. There are two soil layers in the random field with the same thickness. Assume that the mean of soil hydraulic conductivity is 1 × 10−5 cm/s and the fluctuation scale of the random field in any direction is fixed at 3 m. Different COV values result in changes in the dispersion degree and spatial structure of hydraulic conductivity, while the mean of hydraulic conductivity of soil is the same. In order to show this spatial distribution characteristic, different COV values were selected to generate the hydraulic conductivity random field: COV of 0.1 and 0.3 (Figure 1a) and 0.2 and 0.3 (Figure 1b). There are 1600 elements in the random field, and the element size is a 1 m square. Hydraulic conductivity is expressed in gray scale; from pure black to pure white, the scale represents hydraulic conductivity from minimum to maximum. According to the principle of minimum potential energy, groundwater mainly flows through light colored units.
In Figure 1a, the gray levels of the upper soil units are close to each other, the overall color difference is small, the distribution is relatively uniform, and the gray levels of the units are far away from pure white or pure black, which shows that the gray value is concentrated in the middle part. Correspondingly, in this case, the value of hydraulic conductivity changes little, and the deviation from its mean value is small. The COV of upper soil is 0.1, which means that the deviation of unit hydraulic conductivity from its mean value is small. There are obvious differences in the gray level of subsoil units; it is close to pure black or pure white, and the gray difference increases. Based on the definition of COV, the greater the value, the higher the deviation of hydraulic conductivity from its mean value. Relatively, the unit can be assigned a larger or smaller value of hydraulic conductivity. In Figure 1a, the random field reflects the spatial distribution of hydraulic conductivity, and similar results can be seen in Figure 1b. In this model, the dispersion degree and spatial structure of hydraulic conductivity of embankment soil will change, although the mean values are equal.
Obviously, it is difficult to characterize the spatial variation of hydraulic conductivity in the fractal dimension, which limits the application of fractal theory in engineering, particularly embankment engineering. Therefore, embankment soil is divided into several basic units, the soil hydraulic conductivity random field is generated, the fractal dimension of each soil unit is calculated, and then the influence of fractal dimension on hydraulic gradient is studied, which has positive significance for engineering practice.

3. Case Study

In recent years, within the general trend of global warming, extreme weather has occurred frequently. In order to ensure the stability of lake embankments, the embankments of Shijiu and Gucheng Lakes in Nanjing City, Jiangsu Province, are reinforced. During the reinforcement process, the hydraulic conductivity of many embankment soil samples was measured. Based on this, the embankment section of Shijiu Lake was selected for analysis. It is located at K19 + 400 stake in Lishui District, near Lianhewei drainage station (shown in Figure 2). In Lishui District, the width of the embankment top is generally 5–7 m, and some embankment tops are earth roads. The lower part of the upstream slope has dry masonry and riprap toe protection, and the upper part has mortar masonry or concrete protection. The downstream slope is basically natural grassland or shrub, and there are berms on the downstream slope of some embankment sections; there is a mortar masonry wave wall approximately 0.3–1.5 m in size near the water.
Due to the length of time the embankment has been there, the upstream–downstream slopes have been eroded to varying degrees. In studying the randomness of hydraulic conductivity, the analysis model was established according to the engineering practice without considering the erosion effect on both sides. In this section, the embankment soil is not layered and is composed of heavy silty loam.
Figure 3 shows the embankment model of Shijiu Lake. According to the engineering practice, the generalized model of the lake embankment was established. The width of the embankment top and bottom was 6 and 48 m, the vertical distance from the top to the foundation was 8 m, and the model was 35 m long in the z-direction. The model was divided into 26,880 elements and 59,976 nodes. The dimensions of the elements in the y-and z-directions were 0.5 and 1 m, and the cell size was not fixed in the x-direction. Six sections were chosen to study the variation of soil hydraulic conductivity and fractal dimension: z = 5, 10, 15, 20, 25, and 30.
In this study, through on-site borehole sampling, indoor permeability tests were carried out on several samples of embankments in Lishui District. A total of 54 undisturbed soil samples were obtained in the soil layers inside and outside the impervious core of the embankment section. The undisturbed soil sample was 20 cm long. Each sample was measured three times, and then the hydraulic conductivity of embankment soil samples in Lishui District was obtained. The mean value measured in the test was compared with the hydraulic conductivity provided in the survey report and analyzed to determine the hydraulic conductivity of each soil layer as shown in Table 1.
Based on the data in Table 1, a random field of hydraulic conductivity was generated, with a mean of 3.32 × 10−6 cm/s and COV of 0.3, and then the fractal dimension of soil element was calculated according to the relationship between soil fractal dimension and hydraulic conductivity. There were three important boundary conditions: (1) the upstream slope was permeable, and the total head was 6 m; (2) the downstream slope was permeable, and the total head was 0 m; (3) without considering the influence of unsaturated zone, the bottom of the embankment was impervious.

4. Results and Discussion

A gray scale diagram of hydraulic conductivity and fractal dimension of six sections in the embankment are shown in Figure 4. Figure 4a shows the distribution law of the hydraulic conductivity of each unit in Section 1. In the section, the gray values of the units in the middle and upper left areas are similar, showing darker colors, and the hydraulic conductivity in this area was low. On both sides of the section, the gray values of the units were close to each other and the color was light, which means the hydraulic conductivity of the soil was high. It can be seen that dark and light soil units appeared as irregular blocks. In the junction area of blocks, there was no sudden change in gray level, but an obvious transition state, gradually changing from pure white to pure black.
Figure 4b shows the distribution law of the fractal dimension of each unit in Section 1. In this section, the gray values of the units in the middle and upper left regions were similar, showing lighter colors. The unit fractal dimension in this region was large; according to the statistics, the soil fractal dimension was mainly distributed in the interval (1.7, 1.8), which was close to two. In a plane, the soil fractal dimension was distributed in the interval (1, 2); therefore, in this part of the region, it was large. On both sides of the section, the gray values of the units were also close to each other, but the color in this part was dark. According to the statistics, in this part of the area, the soil fractal dimension was mainly distributed in the interval (1.55, 1.6), which means that it was small. At the same time, it can be seen that the dark and light units appeared as irregular blocks. In the junction area of the blocks, there was no sudden change in the gray, but an obvious shift from pure white to pure black.
Comparing Figure 4a,b, the unit distribution area with dark color in Figure 4a was close to the unit distribution area with light color in Figure 4b, and the unit distribution area with light color in Figure 4a was close to the unit distribution area with medium dark color in Figure 4b, indicating that the unit gray levels of hydraulic conductivity and fractal dimension present opposite characteristics. In Figure 4a, the soil units with low hydraulic conductivity were represented by dark color, and the units with high hydraulic conductivity were represented by light color. In Figure 4b, the soil units with small fractal dimension were represented by dark color, and the units with large fractal dimension were represented by light color. In a unit at the same point, when the hydraulic conductivity was low, the corresponding fractal dimension was large, and when the hydraulic conductivity was high, the corresponding fractal dimension was small.
The reasons for this phenomenon are as follows: The fractal dimension of two-dimensional Euclidean space is 2. For the soil distributed in this space, considering the influence of pore structure, the fractal dimension is 1 < D < 2 . Assuming that the structure of pores in the soil tends to be consistent, a single-direction seepage channel is formed, and groundwater can only flow through this channel. Based on the meaning of fractal, the fractal dimension of soil is infinitely close to 1. In this case, the seepage path of water is a straight line. In the flow process, it is rarely hindered by soil particles or soil particle agglomerates, the flow rate is high, and the hydraulic conductivity of soil can reach the maximum value.
Correspondingly, it is assumed that there is no rule governing the structure of pores. At any point, the seepage channel can point in any direction and water can flow in any direction. In this case, the fractal dimension of the soil is infinitely close to 2, and the seepage path of water is no longer a straight line but a curve related to the head difference. In the process of flow, soil particles or soil particle agglomerates form the greatest obstruction to seepage, the hydraulic conductivity and flow rate are reduced and can reach the minimum value. In nature, the two extreme working conditions are almost impossible. In the section, when the fractal dimension of soil increases, the hydraulic conductivity decreases, and when the fractal dimension decreases, the hydraulic conductivity increases.
Figure 4c,d show similar laws. Comparing Figure 4a,c, due to the randomness of hydraulic conductivity, although the laws are consistent, its distribution in the two sections differs. In Figure 4e,g, the dark areas are widely distributed, which indicates that most of the unit hydraulic conductivity is low in the two sections. In Figure 4f,h, there are many units with light color, which means that the fractal dimension of most units is large. Consistent results can also be seen in Figure 4i–l. Obviously, hydraulic conductivity can be explained by fractal theory; however, it is not an ideal fractal, which leads to the fractal dimension of soil properties always being greater than 1.5 [31,38]. Similar explanations can be found in this study.
Figure 5 shows the variations in hydraulic conductivity with fractal dimension in different sections. In the sections, hydraulic conductivity decreases with increased fractal dimension. For the same fractal dimension, some of the unit hydraulic conductivity is not unique, and its value changes within a range. The reason is that the unit permeability is affected by the maximum pore diameter and porosity, and these are random variables, which will affect the fractal dimension and hydraulic conductivity values.
The fractal dimension values in each section are in the interval [1.5, 1.8], which is far from 1.0. In this condition, it is difficult to form a single-direction seepage channel in the soil, and water can flow in multiple directions. The fractal dimension is close to 2.0, which means that the seepage channel can point in almost any direction in the plane, and the water can flow in any direction in the soil. The seepage path of water is not a straight line but is in a complex and changeable state.
In each section, the unit hydraulic conductivity is mostly distributed in the interval [2 × 10−6, 4 × 10−6], and only a few hydraulic conductivity units are distributed outside this interval. This phenomenon conforms to the objective law. In this research, the mean of hydraulic conductivity is 3.32 × 10−6 cm/s, COV is 0.3, and the scale of fluctuations is 6 and 3 m, respectively. When the hydraulic conductivity random field is discretized, COV is within a certain range, and the value after discretization deviates little from the mean value. With increased COV, the value after discretization will be distributed in a larger range. In Figure 5b,c,e, the discrete points show a slight downward bending phenomenon, which means that the unit permeability may decrease slowly with increased fractal dimension.
In order to study the relationship between the fractal dimension and hydraulic gradient of soil, four soil units were selected near the upstream slope angle of each section, which were in the same position. The random seepage field of the embankment was calculated by ABAQUS software, and the hydraulic gradient at each point was obtained. The variation in hydraulic gradient with fractal dimension in different sections is shown in Figure 6. With increased fractal dimension, the hydraulic gradient increased accordingly, because the soil fractal dimension was inversely proportional to the hydraulic conductivity; when the fractal dimension increased the unit permeability decreased, which means that water flows the same distance in the soil, encounters more resistance, and needs more potential energy, resulting in increased hydraulic gradient.

5. Conclusions

This study focused on constructing and applying a seepage–fractal model of embankment soil. The analytical formula of fractal dimension was deduced, and a random field of hydraulic conductivity was generated to calculate the fractal dimension of soil units. With the seepage–fractal model, the internal relationship between embankment hydraulic conductivity, hydraulic gradient, and fractal dimension was revealed. There are four main conclusions:
(1)
The proposed seepage–fractal model of embankment soil is suitable for porous soil media under laminar flow. The influencing factors of hydraulic conductivity mainly include pore size, fractal dimension, and fluid viscosity coefficient, and fractal dimension is the main factor;
(2)
Hydraulic conductivity is inversely proportional to fractal dimension. Increased fractal dimension will reduce the connectivity of soil pores in a single direction, increase the seepage resistance of water, and reduce the hydraulic conductivity. Decreased fractal dimension will lead to consistency of seepage channels in the soil, limited seepage direction of water, decreased resistance in the seepage direction, and increased hydraulic conductivity;
(3)
Increased fractal dimension leads to decreased hydraulic conductivity, increased potential energy consumption through the same seepage path, and increased hydraulic gradient. When the seepage resistance increases further, the seepage path changes, and the water will bypass units with high fractal dimension and flow through units with low fractal dimension;
(4)
In the seepage–fractal model of embankment soil, the fractal dimension and hydraulic conductivity had significant field characteristics. Units with similar attributes formed agglomerates within which the soil interacts, and the attributes tended to be consistent. On the outside, the interaction between clusters shows obvious transition.

Author Contributions

Investigation, methodology, writing—original draft preparation, X.Z.; data curation, B.Y.; conceptualization, S.Y.; funding acquisition, Z.S.; writing—review and editing, project administration, D.F. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by [Key Scientific Research Project of Colleges and Universities in Henan Province] grant number [22B570002 and 21A410003] And [National Natural Science Foundation of China] grant number [U2243210] And [Natural Science Foundation of Henan] grant number [222300420281].

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The datasets generated during the current study are available from the corresponding author on reasonable request.

Acknowledgments

The authors would like to thank all the anonymous referees for their constructive comments and suggestions.

Conflicts of Interest

No conflict of interest exits in the submission of this manuscript, and the manuscript was approved by all authors for publication.

References

  1. Peli, T. Multiscale fractal theory and object characterization. JOSA A 1990, 7, 1101–1112. [Google Scholar] [CrossRef]
  2. Perfect, E.; Kay, B.D. Fractal theory applied to soil aggregation. Soil Sci. Soc. Am. J. 1991, 55, 1552–1558. [Google Scholar] [CrossRef]
  3. Wang, L.; Zeng, X.; Yang, H.; Lv, X.; Guo, F.; Shi, Y.; Hanif, A. Investigation and application of fractal theory in cement-based materials: A review. Fractal Fract. 2021, 5, 247. [Google Scholar] [CrossRef]
  4. Dıaz-Zorita, M.; Perfect, E.; Grove, J.H. Disruptive methods for assessing soil structure. Soil Tillage Res. 2002, 64, 3–22. [Google Scholar] [CrossRef]
  5. Tyler, S.W.; Wheatcraft, S.W. Fractal scaling of soil particle-size distributions: Analysis and limitations. Soil Sci. Soc. Am. J. 1992, 56, 362–369. [Google Scholar] [CrossRef]
  6. Romero, E.; Simms, P.H. Microstructure investigation in unsaturated soils: A review with special attention to contribution of mercury intrusion porosimetry and environmental scanning electron microscopy. Geotech. Geol. Eng. 2008, 26, 705–727. [Google Scholar] [CrossRef]
  7. Gimenez, D.; Perfect, E.; Rawls, W.J.; Pachepsky, Y. Fractal models for predicting soil hydraulic properties: A review. Eng. Geol. 1997, 48, 161–183. [Google Scholar] [CrossRef]
  8. Ju, X.; Jia, Y.; Li, T.; Gao, L.; Gan, M. Morphology and multifractal characteristics of soil pores and their functional implication. Catena 2021, 196, 104822. [Google Scholar] [CrossRef]
  9. Wang, J.; Qin, Q.; Guo, L.; Feng, Y. Multi-fractal characteristics of three-dimensional distribution of reconstructed soil pores at opencast coal-mine dump based on high-precision CT scanning. Soil Tillage Res. 2018, 182, 144–152. [Google Scholar] [CrossRef]
  10. Yang, M.; Fu, Y.; Li, G.; Ren, Y.; Li, Z.; Ma, G. Microcharacteristics of soil pores after raindrop action. Soil Sci. Soc. Am. J. 2020, 84, 1693–1704. [Google Scholar] [CrossRef]
  11. Yang, Y.; Wu, J.; Zhao, S.; Han, Q.; Pan, X.; He, F.; Chen, C. Assessment of the responses of soil pore properties to combined soil structure amendments using X-ray computed tomography. Sci. Rep. 2018, 8, 1–10. [Google Scholar] [CrossRef] [Green Version]
  12. Chang, D.S.; Zhang, L.M. Extended internal stability criteria for soils under seepage. Soils Found. 2013, 53, 569–583. [Google Scholar] [CrossRef] [Green Version]
  13. Liu, Q.Q.; Li, J.C. Effects of water seepage on the stability of soil-slopes. Procedia IUTAM 2015, 17, 29–39. [Google Scholar] [CrossRef] [Green Version]
  14. Yang, Y.; Sun, G.; Zheng, H. Modeling unconfined seepage flow in soil-rock mixtures using the numerical manifold method. Eng. Anal. Bound. Elem. 2019, 108, 60–70. [Google Scholar] [CrossRef]
  15. Wilson, G.V.; Periketi, R.K.; Fox, G.A.; Dabney, S.M.; Shields, F.D.; Cullum, R.F. Soil properties controlling seepage erosion contributions to streambank failure. Earth Surf. Processes Landf. J. Br. Geomorphol. Res. Group 2007, 32, 447–459. [Google Scholar] [CrossRef]
  16. Hu, S.; Tian, C.; Gan, Y. Determination and calculation of soil permeability coefficient. Trans. Chin. Soc. Agric. Eng. 2011, 27, 68–72. [Google Scholar]
  17. Elhakim, A.F. Estimation of soil permeability. Alex. Eng. J. 2016, 55, 2631–2638. [Google Scholar] [CrossRef] [Green Version]
  18. O’Kelly, B.C.; Nogal, M. Determination of soil permeability coefficient following an updated grading entropy method. Geotech. Res. 2020, 7, 58–70. [Google Scholar] [CrossRef]
  19. Kozlowski, T.; Ludynia, A. Permeability coefficient of low permeable soils as a single-variable function of soil parameter. Water 2019, 11, 2500. [Google Scholar] [CrossRef] [Green Version]
  20. Lu, M.M.; Xie, K.H.; Guo, B. Consolidation theory for a composite foundation considering radial and vertical flows within the column and the variation of soil permeability within the disturbed soil zone. Can. Geotech. J. 2010, 47, 207–217. [Google Scholar] [CrossRef]
  21. Dathe, A.; Eins, S.; Niemeyer, J.; Gerold, G. The surface fractal dimension of the soil–pore interface as measured by image analysis. Geoderma 2001, 103, 203–229. [Google Scholar] [CrossRef]
  22. Ghanbarian-Alavijeh, B.; Millán, H.; Huang, G. A review of fractal, prefractal and pore-solid-fractal models for parameterizing the soil water retention curve. Can. J. Soil Sci. 2011, 91, 1–14. [Google Scholar] [CrossRef]
  23. Sun, X.; She, D.; Wang, H.; Fei, Y.; Gao, L. Modelling soil hydraulic properties with an improved pore-solid fractal (PSF) model through image analysis. Eur. J. Soil Sci. 2022, 73, e13156. [Google Scholar] [CrossRef]
  24. Wise, W.R. A new insight on pore structure and permeability. Water Resour. Res. 1992, 28, 189–198. [Google Scholar] [CrossRef]
  25. Wang, L.; Song, X.; Yang, H.; Wu, B.; Mao, W. Pore structural and fractal analysis of the effects of MgO reactivity and dosage on permeability and F–T resistance of concrete. Fractal Fract. 2022, 6, 113. [Google Scholar] [CrossRef]
  26. Su, H.; Zhang, Y.; Xiao, B.; Huang, X.; Yu, B. A fractal-monte carlo approach to model oil and water two-phase seepage in low-permeability reservoirs with rough surfaces. Fractals 2021, 29, 2150003. [Google Scholar] [CrossRef]
  27. Chang, B.; Du, C.; Sun, M.; Lin, Y.; Wang, Y.; Chu, X.; Zhang, L.; He, J. Mesoscopic Seepage Simulation and Analysis of Unclassified Tailings Pores Based on 3D Reconstruction Technology. ACS Omega 2021, 6, 14309–14316. [Google Scholar] [CrossRef]
  28. Ye, W.; Hu, J.; Ma, F. Centrifuge model study on the influence of desiccation cracks on the seepage behavior of upstream clay anti-seepage system subjected to abrupt flood. Bull. Eng. Geol. Environ. 2021, 80, 5075–5090. [Google Scholar] [CrossRef]
  29. Wang, L.; Li, G.; Li, X.; Guo, F.; Tang, S.; Xiao, L.; Hanif, A. Influence of reactivity and dosage of MgO expansive agent on shrinkage and crack resistance of face slab concrete. Cem. Concr. Compos. 2022, 126, 104333. [Google Scholar] [CrossRef]
  30. Yang, H.M.; Zhang, S.M.; Lei, W.; Chen, P.; Shao, D.K.; Tang, S.W.; Li, J.Z. High-ferrite Portland cement with slag: Hydration, microstructure, and resistance to sulfate attack at elevated temperature. Cem. Concr. Compos. 2022, 130, 104560. [Google Scholar] [CrossRef]
  31. Burrough, P.A. Multiscale sources of spatial variation in soil. I. The application of fractal concepts to nested levels of soil variation. J. Soil Sci. 1983, 34, 577–597. [Google Scholar] [CrossRef]
  32. Shah, K.; Arfan, M.; Mahariq, I.; Ahmadian, A.; Salahshour, S.; Ferrara, M. Fractal-Fractional Mathematical Model Addressing the Situation of Corona Virus in Pakistan. Results Phys. 2020, 19, 103560. [Google Scholar] [CrossRef] [PubMed]
  33. Yu, B. Analysis of Flow in Fractal Porous Media. Appl. Mech. Rev. 2008, 61, 050801. [Google Scholar] [CrossRef]
  34. Freeze, R.A. A stochastic-conceptual analysis of one-dimensional groundwater flow in nonuniform homogeneous media. Water Resour. Res. 1975, 11, 725–741. [Google Scholar] [CrossRef]
  35. Hoeksema, R.J.; Kitanidis, P.K. Analysis of the spatial structure of properties of selected aquifers. Water Resour. Res. 1985, 21, 563–572. [Google Scholar] [CrossRef]
  36. Sudicky, E.A. A natural gradient experiment on solute transport in a sand aquifer: Spatial variability of hydraulic conductivity and its role in the dispersion process. Water Resour. Res. 1986, 22, 2069–2082. [Google Scholar] [CrossRef]
  37. Huang, H.; Hu, B.X.; Wen, X.H.; Shirley, C. Stochastic inverse mapping of hydraulic conductivity and sorption partitioning coefficient fields conditioning on nonreactive and reactive tracer test data. Water Resour. Res. 2004, 40, 1–16. [Google Scholar] [CrossRef]
  38. Shou, D.; Fan, J.; Ding, F. A difference-fractal model for the permeability of fibrous porous media. Phys. Lett. A 2010, 374, 1201–1204. [Google Scholar] [CrossRef]
Figure 1. Gray scale diagram of hydraulic conductivity: (a) COV of hydraulic conductivity (k) = 0.1 and 0.3; (b) COV of k = 0.2 and 0.3.
Figure 1. Gray scale diagram of hydraulic conductivity: (a) COV of hydraulic conductivity (k) = 0.1 and 0.3; (b) COV of k = 0.2 and 0.3.
Fractalfract 06 00277 g001
Figure 2. Location of Shijiu Lake.
Figure 2. Location of Shijiu Lake.
Fractalfract 06 00277 g002
Figure 3. Embankment model of Shijiu Lake: (a) dimensions of section; (b) mesh generation of section; (c) location of section in embankment model; (d) mesh generation of embankment.
Figure 3. Embankment model of Shijiu Lake: (a) dimensions of section; (b) mesh generation of section; (c) location of section in embankment model; (d) mesh generation of embankment.
Fractalfract 06 00277 g003
Figure 4. Distribution of hydraulic conductivity (k) and fractal dimension (D) of sections 1−6 in the embankment: (a) k in section 1; (b) D in section 1; (c) k in section 2; (d) D in section 2; (e) k in section 3; (f) D in section 3; (g) k in section 4; (h) D in section 4; (i) k in section 5; (j) D in section 5; (k) k in section 6; (l) D in section 6.
Figure 4. Distribution of hydraulic conductivity (k) and fractal dimension (D) of sections 1−6 in the embankment: (a) k in section 1; (b) D in section 1; (c) k in section 2; (d) D in section 2; (e) k in section 3; (f) D in section 3; (g) k in section 4; (h) D in section 4; (i) k in section 5; (j) D in section 5; (k) k in section 6; (l) D in section 6.
Fractalfract 06 00277 g004
Figure 5. Variation in the hydraulic conductivity with fractal dimension in different sections: (a) section 1; (b) section 2; (c) section 3; (d) section 4; (e) section 5; (f) section 6.
Figure 5. Variation in the hydraulic conductivity with fractal dimension in different sections: (a) section 1; (b) section 2; (c) section 3; (d) section 4; (e) section 5; (f) section 6.
Fractalfract 06 00277 g005
Figure 6. Variation in hydraulic gradient with fractal dimension in different sections: (a) point 1; (b) point 2; (c) point 3; (d) point 4.
Figure 6. Variation in hydraulic gradient with fractal dimension in different sections: (a) point 1; (b) point 2; (c) point 3; (d) point 4.
Fractalfract 06 00277 g006
Table 1. Parameters of embankment model.
Table 1. Parameters of embankment model.
SoilkCOV θ x θ y
Heavy silty loam3.32 × 10−6 cm/s0.36 m3 m
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Zhao, X.; Yang, B.; Yuan, S.; Shen, Z.; Feng, D. Seepage–Fractal Model of Embankment Soil and Its Application. Fractal Fract. 2022, 6, 277. https://doi.org/10.3390/fractalfract6050277

AMA Style

Zhao X, Yang B, Yuan S, Shen Z, Feng D. Seepage–Fractal Model of Embankment Soil and Its Application. Fractal and Fractional. 2022; 6(5):277. https://doi.org/10.3390/fractalfract6050277

Chicago/Turabian Style

Zhao, Xiaoming, Binbin Yang, Shichong Yuan, Zhenzhou Shen, and Di Feng. 2022. "Seepage–Fractal Model of Embankment Soil and Its Application" Fractal and Fractional 6, no. 5: 277. https://doi.org/10.3390/fractalfract6050277

APA Style

Zhao, X., Yang, B., Yuan, S., Shen, Z., & Feng, D. (2022). Seepage–Fractal Model of Embankment Soil and Its Application. Fractal and Fractional, 6(5), 277. https://doi.org/10.3390/fractalfract6050277

Article Metrics

Back to TopTop