Next Article in Journal
Transient Dynamics of a Fractional Fisher Equation
Previous Article in Journal
Changes in Pore Structure and Fractal Characteristics of Solvents Pretreated High-Rank Coal under Supercritical CO2
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Characterization of the Fine-Scale Evolution of Damage in Shale under the Influence of Two-Way Stress Differences Based on CT Images and Fractal Theory—The Example of the Anba Dyke in the Wufeng–Longmaxi Formation

1
College of Resource and Environmental Engineering, Guizhou University, Guiyang 550025, China
2
Guizhou Energy Industry Research Institute Co., Ltd., Guiyang 550025, China
3
College of Civil Engineering, Guizhou University, Guiyang 550025, China
*
Author to whom correspondence should be addressed.
Fractal Fract. 2024, 8(3), 142; https://doi.org/10.3390/fractalfract8030142
Submission received: 29 December 2023 / Revised: 2 February 2024 / Accepted: 27 February 2024 / Published: 29 February 2024

Abstract

:
To better understand the influence of different levels of two-way stress differences on the development of damage in Anchang diametral laminar shale in the northern Qianbei area, a numerical model of laminar shale with a representative fine-scale structure was established by using RFPA3D-CT. A triaxial compression test was conducted; a three-dimensional mesoscale fracture box dimension algorithm based on digital images was generated by using MATLAB R2020b; and the fractal characteristics were quantitatively analyzed. The results showed that under the influence of the horizontal stress ratio and two-way stress, the greater the two-way stress is, the more notable the plastic characteristics of specimen damage are, and the higher the residual strength is. The specimens with lower two-way stress exhibited obvious brittle damage characteristics. The difficulty degree of complex fracture network formation increased with the increase in the horizontal tension ratio, and the degree of increase in the fracture network complexity gradually decreased. At a horizontal stress ratio of 1.25, the fractal dimension was the highest, which indicates that the cracks were the most pronounced. Fracture formation after specimen damage was the most common phenomenon. Under the condition of a lower horizontal stress ratio, a large number of fracture structures could be generated in shale specimens after damage, promoting the expansion of natural fractures.

1. Introduction

As the Chinese economy continues to develop at a rapid pace, the supply and demand for oil and natural gas resources are notably increasing [1,2,3,4,5]; thus, there is a pressing need to strengthen national endeavors in oil and gas exploration and development. Shale gas, classified as an unconventional energy source [6,7,8], exhibits tremendous potential as a feasible substitute for traditional oil and gas resources. As a result, it has attracted considerable interest from scholars both domestically and internationally [9,10,11,12]. China possesses abundant supplies of shale gas, and many other nations have intensified their research on and use of shale gas in recent years [13,14,15]. China’s recoverable shale gas volume ranges from about 11.5 × 1012 to 36.1 × 1012 m3 and is among the largest worldwide [16].
Zou [17] et al. noted that hydraulic fractures can intersect natural fractures at high levels of differential stress. Rickman et al. [18] postulated that the presence of a substantial clay mineral content impedes the development of intricate fracture networks, which are more prone to forming in brittle shale characterized by a high mineral content. Olsen et al. [19] demonstrated that shale formations exhibiting low Poisson’s ratios and extremely low permeability are more prone to undergoing complex fracturing. Hou et al. [20] and collaborators further revealed, via comprehensive laboratory experiments, that rocks enriched with a substantial quantity of brittle minerals demonstrate an increased predisposition to the formation of intricate fracture networks in comparison to other rock types. This rock type not only contains discontinuous surfaces that can be activated and connected into a complex fracture network but also provides the large fracture lengths required for notable reservoir transformations [21,22].
However, most of the above studies have focused on shale samples, as well as shale outcrops containing a small number of discontinuities, which do not suitably represent actual layered shale reservoirs. The extension geometry of shallow shale fractures, rather than that of deep shale horizons, has been investigated in several shale outcrop experiments; the latter structures are widely found in the northern part of the Qianbei region. Deep shale horizons are located above 3500 m, and they exhibit different geomechanical properties from shallow shale horizons. These rocks have higher horizontal stress and temperatures above 110 °C [23], and in such deep shale horizons, the confluence of elevated fracture pressure, high density, and diminished brittleness culminates in restricted fracture width, constrained fracture propagation, sand accumulation, and diminished flow conductivity, thus challenging reservoir modification, making it difficult to implement fracturing treatments, and potentially leading to a decrease in shale gas production. To date, there has been a dearth of research investigating approaches to estimating and regulating fracture extension and reservoir remodeling volume in deep formations with elevated stress concentration levels. Therefore, delving into the mechanisms governing fracture extension in deep shale formations is imperative for the precise estimation of and enhancement in reservoir transformation volume. Numerical simulation helps to identify crack extension paths. However, current two-dimensional numerical simulation methods are flawed and inaccurate in assessing reservoir volume modification. Through the simplification of the complex network of fractures within a three-dimensional domain to multiple fractures projected onto a two-dimensional horizontal plane, these methodologies concentrate exclusively on fracture extension at a designated elevation, neglecting the interplay between fractures and structural discontinuities [24].
Presently, a deficiency persists in the utilization of three-dimensional (3D) numerical models incorporating CT scanning within extant simulation investigations aimed at scrutinizing the intricate internal composition of laminated shale formations. Consequently, this investigation employs cutting-edge CT scanning, advanced imaging technology, and the state-of-the-art RFPA3D-CT finite element computation method, thereby facilitating the construction of a sophisticated three-dimensional numerical model to aptly capture the intricacies of weak lamination surfaces. As a result, a comprehensive and accurate 3D numerical model specifically tailored to laminated shale formations in the Zheng’an region is successfully formulated. In order to examine the damage evolution, damage morphology, and mechanical characteristics of laminated shale specimens sourced from the Zheng’an region, a rigorous true triaxial compression test was executed. The primary aim of this experiment was to evaluate the influence of biaxial stress differentials on the mechanical attributes, damage morphology, and fracture intricacies of shale formations.
The utilization of the fractal dimension, serving as an index of fracture network effects, permits the identification of optimal fracturing conditions for the adjustment of fracture networks within hydraulic fracturing operations in shale oil reservoirs. This investigation has the potential to provide a theoretical framework for hydraulic fracturing endeavors within deep, high-stress shale formations manifesting diverse disparities.

2. Geological and Tectonic Background of the Study Area

The Anchang syncline is located in the Wuling depression between the southeastern edge of the Middle and Upper Yangzi Massifs and the Xuefeng Mountain basement detachment orogenic belt (Figure 1), which is the western side of the Sichuan eastern compartmentalized trough fold belt. Tectonically, it is located north of Guizhou Province and is part of the northeastern tectonic deformation area of the Yangzi quasi-platform, Qianbei Plateau, convex Zunyi fault, and Fenggang north-northeast tectonic deformation zone. The northern region of Guizhou Province exhibits notable structural features characterized by prominently developed synclines and fault structures. These folds predominantly exhibit a consistent northeastern to north-northeastern orientation, and they are mainly divided trough folds, with narrow and tightly closed, steeply dipping troughs and banded distributions on the surface. The backslopes are wide and gradually sloping in the form of a box and often exhibit the characteristics of compound structures, while internal extrusion tectonics, such as lenticular bodies, rubbing folds, compressive faults, and slip tectonics, are very well developed. The fold axis is generally an NNE–NE arc-like curve. Fractures in the fold belt are extremely common, with the main strike occurring in the NE–NNE direction, mostly in the core and wings of the dorsal slope, and the fractures are mostly retrograde in nature, with a small number of positive fractures also present. The Anchang dorsal slope is located in the weak–medium deformation zone in northern Yunnan–Guizhou, which is a negative tectonic unit formed against the geotectonic background of the Early Paleozoic era, whereas the positive tectonic attributes of the southern neighboring East Yunnan–Guizhou–China uplift are exactly the opposite, but the formation and development processes of these two structures are complementary and closely related. Deformational tectonics in the deformation area were mainly formed during the Yanshan main fold orogeny in the Early Cretaceous period, and the typical tectonic type is the syncline fold-and-rupture combination with folds, which is a straight torsion tectonic type.

3. Introduction to RFPA3D

Rockfall Process Analysis in Three Dimensions (RFPA3D), an extension of Rockfall Process Analysis in Two Dimensions (RFPA2D) [25,26], is a simulation tool tailored for modeling fracture and damage phenomena in quasi-brittle materials such as rock. RFPA3D adopts the finite element method (FEM) incorporating eight-node isoparametric elements for stress analysis, enabling the depiction of the nonhomogeneous nature of rock materials and intricate characterization of rock mechanical parameter distributions. Furthermore, RFPA3D employs an elastic–brittle intrinsic model to accurately portray the nonlinear damage-phase behavior exhibited by rock.
Equation (1) in RFPA3D is frequently employed to quantitatively assess the maximum tensile strain threshold, which serves as a criterion for identifying the initiation of fracture generation and the occurrence of tensile failure. Furthermore, Equation (2) is employed as the secondary damage threshold, supplying valuable information regarding the occurrence of shear mode damage in a given element [27,28].
ε 3 σ t 0 E
F = σ 1 ε 3 1 + sin φ 1 sin φ σ c 0
An elastic–brittle damage constitutive model is used to describe the stress–strain relationship and consider its residual strength. Each unit is considered elastic, and before reaching the damage threshold, the stress–strain curve of each unit is considered linear elastic. The determination of damage thresholds in RFPA3D relies on the utilization of both the maximum compressive stress (or strain) criterion and the Mohr–Coulomb criterion. Upon reaching the prescribed strength criterion, the onset of damage accumulation ensues within the element, ultimately resulting in the stepwise degradation of the element’s elastic modulus as the damage process advances. The constitutive relationship of damaged components can be expressed as follows:
σ = E ε = E 0 ( 1 D ) ε
where E is the modulus of elasticity of the unit after damage has occurred, E0 is the initial modulus of elasticity of the unit, and D is the damage parameter. With D = 0, the unit exhibits no damage; with D = 1, the unit exhibits a complete damage state; and with 0 < D < 1, the unit exhibits damage. Notably, according to Equation (3), the modulus of elasticity is E = 0 if damage parameter D = 1, which does not correspond to the actual damage to the rock specimen. For this reason, a relatively small E value of 1.0 × 10−5 is specified in the tests, where E, E0, and D are scalar quantities, because RFPA3D assumes that the unit and damage are isotropic.
Figure 2 elucidates the ontological interconnection between the maximum compressive stress (or strain) criterion and microscopic cells during the triaxial compression process, and the corresponding damage variables are given as follows [29,30,31]:
D = { 0 ε > ε t 0 1 σ t r ε E 0 ε t 0 ε ε t u 1 ε < ε t u
In Equation (4), σtr is the residual compressive strength of the unit; σt0 is the ultimate compressive strength of the unit; σrt = λσt0; λ is the residual compressive strength coefficient, and 0 < λ < 1. Moreover, εt0 is the compressive strain of the unit at the elastic limit, which is the strain threshold for compression damage, and εtu is the ultimate compressive strain of the unit when it is completely destroyed. Finally, εtu = ηεt0, where η is the ultimate compressive strain coefficient.
In multiaxial stress scenarios, the influence of the supplementary principal stresses on the maximum or minimum principal strain is duly considered. As the compressive strain surpasses the specified damage strain threshold (εt0), the substitution of the equivalent strain for the compressive strain (ε) in Equation (4) becomes appropriate. The equivalent strain can be expressed as follows:
ε ¯ = ε 1 2 + ε 2 2 + ε 3 2
where ε1, ε2, and ε3 are the principal strains along the three directions. The angle brackets are functions defined as follows:
ε i = { ε i ( ε i 0 , i = 1 , 2 , 3 ) 0 ( ε i > 0 , i = 1 , 2 , 3 )
The stress–strain relationship derived in terms of elastic damage theory can be expressed as follows:
σ i j = { 2 G ε i j + λ δ i j ε k k ε ¯ > ε t 0 σ r t ε ¯ E 0 ( 2 G ε i j + λ δ i j ε k k ) = σ r t ε ¯ ( ε i j 1 + υ + δ i j ε k k ( 1 + υ ) ( 1 2 υ ) ) ε t 0 ε ¯ ε t u 0 ε ¯ < ε t u
where G is the shear modulus of the material, with G = E 0 2 ( 1 + v ) ; η is the Poisson’s ratio of the material; λ is the residual coefficient, with λ = E 0 v ( 1 + v ) ( 1 2 v ) ; and when i = j, σij = 1. Otherwise, σij = 0.

4. CT Scanning Experiments and Numerical Image Characterization

4.1. CT Scanning Process

The shale specimens are transformed into digital images utilizing CT scanning and digital image processing techniques. Grayscale and color data are used to differentiate the geometric and spatial distributions of mesoscopic components. Through the application of a predefined threshold to partition each medium, an image is generated, portraying the heterogeneous nature of the material distribution [32]. As shown in Figure 3a, the scanning test apparatus used was an X-ray 3D microscope-equipped nanoVoxel-4000 system. The primary mechanism of CT scanning is the interaction of radiation with matter. An item scatters a high percentage of input photons as a ray travels through it. This attenuates the intensity of the ray, which is then translated by imaging tools into a CT image (Figure 3b). To obtain the scanned slice shown in Figure 3d, the specimen was scanned every 0.1 mm from top to bottom (Figure 3c). The image has a resolution of 1000 × 1000 pixels and measures 100 mm × 100 mm. Different rock types are represented by various shades of gray. The high-density rock areas, mostly comprising quartz, appear lighter. In contrast, the low-density areas, mainly consisting of laminar structures, appear black. The shale matrix, which has a density between the two, appears gray. The image in Figure 3d was imported into ImageJ2 software for the purpose of segmentation, facilitating the measurement of grayscale values associated with individual pixels along a designated cut line. This segmentation allowed for the creation of a two-dimensional bar distribution of grayscale values, as shown in Figure 3e. The thresholds for binarization were determined as 85 and 100.
The digital representation in Figure 3f illustrates the nonuniformity of the fine-scale structure of the shale specimen after image processing. The experiment employs high-resolution CT scan images to improve the detection of color features occurring at the intersections of multiple material types. This enables more accurate material zoning and parameter assignment in the digital image characterization process.

4.2. Numerical Modeling

The images were numbered sequentially and converted into model data. The data were then stacked along the z-axis and imported into RFPA3D to reconstruct a 3D structural model. A cubic volume unit was used to measure finite quantities. To characterize the different fine-grained materials in the shale specimen, the image was converted into multiple finite element cubes. Given the nonuniformity of the material, it is proposed that the mechanical properties of the fine-grained cells follow the Weibull distribution, expressed as Equation (8) [27].
f ( x ) = m β ( x β ) m 1 e ( x β ) m , x 0
where the variable x denotes the mechanical property parameter associated with the fine-scale rock unit, encompassing characteristics such as modulus of elasticity, strength, Poisson’s ratio, and others. The parameter β represents the average value of each mechanical property parameter for the fine-scale unit, while the coefficient m quantifies the material’s homogeneity, reflecting its degree of uniformity.
The Monte Carlo method was utilized to assign values to the material components based on color and to input the corresponding nonuniformity coefficients according to the mechanical parameters listed in Table 1. Figure 3b shows the numerical model that incorporates the actual fine-scale structure of the material. To ensure consistency, each sample was processed through digital image processing to maintain the same scale. The specimens used in the simulation were 100 mm × 100 mm × 100 mm in size, and a total of five specimens were numerically simulated. Considering the disparity between the tensile strength and compressive strength of brittle materials, RFPA3D-Basic software incorporates a tailored Mohr–Coulomb criterion as the unit damage strength criterion. The loading process involves gradually applying the load up to the predetermined perimeter pressure value, followed by the application of axial displacement once the system perimeter pressure has stabilized. The initial displacement was set to 0.001 mm, with subsequent increments of 0.0005 mm per step. Loading continued until the specimen reached a damaged state, as shown in Figure 4.
The analysis of stress distribution in the Zheng’an shale gas reservoirs of the Anchang slant region, aimed at real shale gas extraction, entailed the investigation of five distinct regional well groups. In order to guarantee precise simulation, five sets of two-dimensional stress loading values were defined, and the corresponding experimental parameters are reported in Table 2.

5. Experimental Results and Discussion

5.1. Mechanical Characterization

In the course of triaxial compaction, the stress–strain curves of shale samples manifested three distinct phases, as illustrated in Figure 5: During the elastic stage, the stress–strain relationship demonstrates linearity with a consistent slope, alongside the gradual advancement of microscopic ruptures. In the nonlinear deformation stage, the slope of the curve decreases as the loading level exceeds approximately 70–80% of the peak strength, indicating that the stress continues to cause specimen damage and decrease the force. The elastic potential energy stored is released, resulting in irreversible deformation, which in turn causes the compaction of some cracks. The stress in the shale specimens of Groups 1, 2, and 5 sharply decreases after the peak phase; the internal microscopic cracks expand through the specimens and form macroscopic main cracks; and the specimens are destabilized and damaged. Similarly, as the post-peak curves decrease almost perpendicularly, irreversible plastic distortion occurs, and the load carrying capacity gradually decreases to 0. The specimens of Groups 3 and 4 still provide a certain load-bearing capacity after failure. Cracks do not penetrate all the way through the specimens at peak stress. When the loading process is continued, new cracks are formed and continue to branch out; the time when the main fracture emerges and when the specimen emits the loudest sounds, therefore, lags behind the ultimate load.

5.2. Analysis of the Rupture Process

The damage to the cells during specimen loading is shown in Figure 6a, where the red color indicates the damaged cells. The isotropic stress slightly affects the specimen, and little major damage occurs during the elastic phase. In nonlinear deformation, fine-scale damage begins to occur at two sites, one corresponding to particle cementation of the laminated surface and one corresponding to the shale matrix, and the number of units associated with internal damage continues to increase. Surface cracks are developed in some of the damaged elements.
In Figure 6d, individual circles symbolize acoustic emission events, with the circle’s diameter correlating with the magnitude of the emitted acoustic emission energy. Each circle’s center demarcates the position of the damaged unit, distinguishing tensile damage in blue and shear damage in red. During the elastic phase, localized tensile stress concentrations arise at the ends of the joints, accompanied by an evenly distributed pressure load across the laminae, which is manifested as scattered micropoint fractures. This indicates that the laminae are fragile and can be easily damaged. As loading proceeds, the energy gradually accumulates, and when the accumulated energy reaches a certain extent, a plethora of sizable red acoustic emission circles are identified within the specimen, indicating the presence of compressive damage within the macroscopic crack. The fracture zone within the high-stress region notably plays a substantial role in the energy release process during specimen fracturing, as illustrated in Figure 6b, and the dominant damage type near the interface between the sample and platen is shear damage. Compressive damage of the fine-grained shale unit leads to the initiation and expansion of fine-grained holes and cracks because the laminar part of the shale specimen cannot withstand high compressive stresses, and the compressive stresses first reach the compressive strength.
In the post-peak phase, the lamellae exhibit the formation of extensive and relatively shallow damage zones, accompanied by interconnected damage among particles. As the number of fine-scale damage points through the joints surpasses a predetermined threshold, the macroscopic damage becomes irreversibly established. The laminae represent the most vulnerable locations, with the most severe damage to the members at the junctions and cracks at the joints. The sequence of damage development in the specimen progresses from the inside to the outside and finally to the surface of the sample.
Under high peripheral pressure, a number of crack zones is generated due to shear damage, and the crack zones extend beyond the region of localized deformation in the direction of the applied load, as depicted in Figure 6c. Under controlled tensile stress, tilted monoplane tensile cracks begin to appear inside the specimen along the nodal direction and extend along the weak surface, and elongation occurs approximately parallel to the maximum principal stress.

6. Fractal Characterization of Shale Fracture Structures

6.1. Calculation of the Fractal Dimension of a Three-Dimensional Fracture Structure

The fractal dimension provides a quantitative measure to assess the intricacy and structural characteristics of shale fractures [35,36]. Different research objects exhibit fractal characteristics within different scale ranges, so different methods can be used to calculate the fractal dimension in a targeted manner [37]. Lopes and Betrouni classified these methods into three categories, one of which includes the box dimension method [38,39]. Because the algorithm for determining the box dimension value is simple, is more suitable for analyzing digital images, and has a clearer physical interpretation, it is commonly utilized in the analysis of crack morphology [40,41,42].
In the model constructed in RPFA3D, since the shale specimen comprises multiple cubic units of the same size, each pixel has a corresponding grayscale value [43,44]. In space, cubes of size m×n×h pixel points are divided (as shown in Figure 7). Therefore, a cubic box can be selected to cover the crack model in three dimensions. The bisection method is chosen to determine the side lengths of the small cubes.
δ k = 1 C k 1 | 0 < k ln min ( m , n , h ) ln c | + 1 k = 0 , 1 , 2 , 3
In RPFA3D, the fracture distribution map of each slice can be generated and exported. The spatial arrangement of cracks in the image was analyzed using the principles of 3D reconstruction to assess their shapes. In Figure 8, the cracks in the postprocessed slices obtained with RFPA3D are generated via a custom procedure to determine the coordinates of the pixels with cracks.
The size of the cube directly impacts the number of small cubes (Nk) required to encompass the crack space. Following image preprocessing procedures, a binary representation of the crack image is produced. The total number of cubes necessary to fully cover the entire crack space is calculated by multiplying it by the constant C. For any value of c, a series of data (1/ck−1, Nk) can be obtained. In the double logarithmic coordinate system, the straight line represented by (1/ck−1, Nk) reflects the fractal characteristics. In the model constructed by RPFA3D, since the shale specimen comprises multiple cubic units of the same size, a cubic box is selected to cover the three-dimensional fracture model.
lg c k 1 = lg ( k )
D = lim r 0 lg N ( k ) lg ( k )

6.2. Effect of the Horizontal Stress Ratio

According to the definition of the cube coverage method, an increase in the complexity of cracks in shale specimens leads to a higher fractal dimension [45,46]. From the fracture network morphology and the fractal dimensions of the shale specimens under each horizontal stress difference condition (Figure 9), the fractal dimensions of the fractal networks based on the cube coverage method range from about 2.339 to 2.517. By analyzing the fractal dimensions of the damaged fracture networks under varying horizontal stresses, it can be determined that within the same specimen, there is significant variation in the fractal size of the cracks after failure under the different horizontal stress levels. For horizontal stress ratios below 1.282, an incremental growth in the horizontal stress ratio is associated with a gradual reduction in the rate of change in the fractal dimension. Within the range of 1.278 to 1.331 for the horizontal stress ratio, a notable decline in the rate of change in the fractal dimension is observed. This observation implies that elevating the horizontal stress ratio within this interval hinders the evolution of a sophisticated fracture network, leading to a gradual decline in the intricacy of the fracture network. At a horizontal stress ratio of 1.25, the fractal dimension is the highest, and most cracks are formed after the specimens are damaged. At lower horizontal stress ratios, a large number of fracture structures are formed after shale sample damage, and further expansion of natural fractures is promoted.
Compared with the results for a horizontal tension of 1.25, the fractal dimension for a horizontal tension of 1.282 is lower, and relatively few cracks are formed after specimen damage. However, there are cracks along the direction of the laminar surface, which is favorable for achieving connectivity with natural fractures and promotes the formation of a single crack (Figure 10). At a high horizontal stress ratio, no single crack produces fine crack branches, except for the main crack, and the fractures do not extend along the laminar surface. Thus, with the increase in the horizontal stress ratio, the fractal dimension of the crack decreases. At tension levels of 1.331 and 1.355 in the horizontal direction, the manipulation of crack orientation poses a notable challenge. At this time, it is more difficult for the cracks to become connected and for them to communicate with more natural cracks and to induce the expansion of the natural cracks. Moreover, the final reconstructed crack network showed that the main fractures perpendicular to the minimal horizontal main tension are more likely to be formed at this time. In other words, when the horizontal tension is relatively low, shale damage is more prone to forming a complex network of fractures, while an excessive horizontal stress ratio promotes the creation of a single giant crack after damage, which is unfavorable for the development of a complex fracture network.

7. Conclusions

Through shale true triaxial compression simulation experiments, fractal theory and 3D CT reconstruction technology were utilized to study the fracture network complexity, and the conclusions below were drawn.
(1) The macrofracture pattern is consistent with the cumulative development of microfractures. The damage to the specimen spreads from the interior to the exterior and eventually extends to the surface. Influenced by the horizontal stress ratio and two-way stress, the higher the two-way stress is, the more notable the plastic characteristics are, and the higher the residual stress of the sample is after damage. Moreover, all the specimens with lower two-way stress exhibit obvious brittle damage characteristics.
(2) Compared with the structure obtained with the 2D numerical simulation, a 3D structure provides a more intuitive and precise depiction of the internal fracture morphology and cell behavior throughout deformation and damage mechanisms.
(3) Fine-scale nonhomogeneity and laminar defects exert a substantial influence on molding the stress distribution, crack propagation trajectory, damage mechanism, and overarching fracturing progression. The strength of shale material exhibits obvious nonuniformity and statistical characteristics, i.e., the initial defects, such as gaps and other areas of low strength, are destroyed first, followed by high-strength regions.
As the horizontal tensile stress ratio increases, the complexity of developing intricate crack networks also escalates, and the extent of the increase in the fracture network complexity gradually decreases. At a horizontal stress ratio of 1.25, the highest fractal dimension is observed, which suggests that the cracks produced after specimen damage are the most well developed. Under a reduced horizontal stress ratio, shale specimens display abundant fracture patterns after damage, thereby promoting the propagation of pre-existing natural fractures.

Author Contributions

Writing—original draft preparation, B.L., Conceptualization, B.L. and R.Y.; methodology, B.L. and Z.W.; software, Z.W.; validation, H.J. and X.L.; formal analysis, B.L. and X.L.; writing—review and editing, B.L., R.Y. and H.J.; project administration, Z.W. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China, grant numbers 51964007, 51774101, and 52104080; Guizhou Science and Technology Fund, grant numbers [2020]4Y046, [2019]1075, and [2018]1107; and Key technologies and engineering tests for shale gas benefit development in Guizhou Province Qiankehe Strategic Prospecting [2022]ZD005.

Data Availability Statement

This study entailed the use of publicly available datasets. The details and sources of these datasets are provided in the relevant literature cited within the article.

Conflicts of Interest

Authors Baofeng Lan and Haishen Jiang were employed by the company Guizhou Energy Industry Research Institute Co., Ltd. The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

References

  1. Fa, G.F.; Wang, Z.Q.; Gao, F.; Wang, Z.S.; Li, Z.Y. Global Natural Gas Production and Development Status, Supply and Demand Situation Analysis and Prospects. In Proceedings of the International Field Exploration and Development Conference 2021; Springer: Singapore, 2022; Volume 8, p. 73. [Google Scholar]
  2. Zheng, T.; Liu, X.; Yang, Z.; Luo, Y.; Zhao, X. Identification of Seepage Mechanisms for Natural Gas Huff-n-Puff and Flooding Processes in Hydrophilic Reservoirs with Low and Ultra-Low Permeabilities. J. Energy Resour. Technol. 2021, 143, 063004. [Google Scholar] [CrossRef]
  3. Chen, G.W. On China’s Oil and Gas Security in Light of “the One Belt and One Road Initiative”. J. Jiangnan Soc. Univ. 2017. [Google Scholar]
  4. Hu, T.; Wu, G.; Xu, Z.; Pang, X.; Liu, Y.; Yu, S. Potential resources of conventional, tight, and shale oil and gas from Paleogene Wenchang Formation source rocks in the Huizhou Depression. Adv. Geo-Energy Res. 2022, 6, 402–414. [Google Scholar] [CrossRef]
  5. Guoxin, L.I.; Zhu, R.; Zhang, Y.; Chen, Y.; Cui, J.; Jiang, Y.; Kunyu, W.U.; Sheng, J.; Xian, C.; Liu, H. Geological characteristics, evaluation criteria and discovery significance of Paleogene Yingxiongling shale oil in Qaidam Basin, NW China. Pet. Explor. Dev. 2022, 49, 21–36. [Google Scholar]
  6. Liu, H.; Meng, S.; Su, J.; Zhang, G.; Chen, L. Reflections and suggestions on the development and engineering management of shale gas fracturing technology in China. Nat. Gas Ind. B 2019, 6, 539–545. [Google Scholar] [CrossRef]
  7. Dong, H.; Guang, Y.; Zhi, Y.; Tian-Quan, Y.; Rong, B.; Yu-Cong, L.; Hong-Ming, D. New understanding and development potential of tight oil exploration and development in Sichuan Basin. Nat. Gas Geosci. 2019, 30, 1212–1221. [Google Scholar]
  8. Cui, H.; Wu, Z.; Li, L.; Wang, J.; Song, S.; Sun, S.; Zuo, Y.; Liu, H.; Lou, Y. Numerical Test Study of the Microscale Failure Modes and Fractal Analysis of Lower Cambrian Shale Based on Digital Images. Adv. Civ. Eng. 2020, 2020, 8819711. [Google Scholar] [CrossRef]
  9. Jia, J.T. China achieved fruitful results in oil-shale gas-coalbed methane exploration and development in 2021. China Geol. 2022, 5, 355–356. [Google Scholar]
  10. Wenke, C.; Wei, L.; Hailong, L. Effect of formation strength anisotropy on wellbore shear failure in bedding shale. J. Pet. Sci. Eng. 2022, 208, 109183. [Google Scholar] [CrossRef]
  11. Chong, Z.; Yao, Q.; Li, X.; Zhu, L.; Tang, C. Effect of rock brittleness on propagation of hydraulic fractures in shale reservoirs with bedding-lanes. Energy Sci. Eng. 2020, 8, 2352–2370. [Google Scholar] [CrossRef]
  12. Wei, T.; Cheng, Z.; Xu, Y.; Jizhao, X.; Yong, S.; Yuzhou, C.; Yangfeng, Z.; Yu, W. Dynamic Brazilian splitting experiment of bedding shale based on continuum-discrete coupled method. Int. J. Impact Eng. 2022, 168, 104289. [Google Scholar]
  13. Wu, Z.; Yang, Y.; Zuo, Y.; Meng, X.; Wang, W.; Lei, W. Damage evolution characteristics of 3D reconstructed bedding-containing shale based on CT technology and digital image processing. Acta Geophys. 2023, 1–17. [Google Scholar] [CrossRef]
  14. Pu, B.; Wang, F.Q.; Dong, D.Z.; Sun, J.B. Challenges of terrestrial shale gas exploration and development from Chang 7 shale in the Ordos Basin. Arab. J. Geosci. 2021, 14, 644. [Google Scholar] [CrossRef]
  15. Sun, C.; Nie, H.; Dang, W.; Chen, Q.; Zhang, G.; Li, W.; Lu, Z. Shale Gas Exploration and Development in China: Current Status, Geological Challenges, and Future Directions. Energy Fuel 2021, 35, 6359–6379. [Google Scholar] [CrossRef]
  16. Yong, C.; Yufeng, G.; Guodong, C.; Wanlin, Y.; Falong, W. On Resource Potential, Distribution and Characteristics of Shale Gas in China. J. Kunming Metall. Coll. 2017, 33, 17–24. [Google Scholar]
  17. Zou, Y.; Huang, M.Z.; Liu, F.Y.; Yang, B.C.; Wang, L.Q.; Wang, F.; Yu, X.H.; Wan, L.; Wan, X.D.; Xu, X.Y. Absence of DICER1, CTCF, RPL22, DNMT3A, TRRAP, IDH1 and IDH2 hotspot mutations in patients with various subtypes of ovarian carcinomas. Biomed Rep. 2015, 3, 33–37. [Google Scholar] [CrossRef] [PubMed]
  18. Rickman, R.; Mullen, M.; Petre, E.; Grieser, B.; Kundert, D. A Practical Use of Shale Petrophysics for Stimulation Design Optimization: All Shale Plays Are Not Clones of the Barnett Shale. In Proceedings of the SPE Annual Technical Conference & Exhibition, Denver, CO, USA, 21–24 September 2008; p. SPE-115258-MS. [Google Scholar]
  19. Olsen, T.N.; Bratton, T.R.; Thiercelin, M.J. Quantifying Proppant Transport for Complex Fractures in Unconventional Formations. Soc. Pet. Eng. 2009, SPE-119300-MS. [Google Scholar] [CrossRef]
  20. Bing, H.; Ruxin, Z.; Yijin, Z.; Weineng, F.; Yeerfulati, M.; Mian, C. Analysis of hydraulic fracture initiation and propagation in deep shale formation with high horizontal stress difference. J. Pet. Ence Eng. 2018, 170, S973227303. [Google Scholar]
  21. Li, X.; Guo, P.; Hu, Y.; Li, S.; Yang, W. Experiment and numerical simulation of hydraulic fracturing in lacustrine shale: Taking the Ordos Basin as an example. Oil Gas Geol. 2023, 44, 1009–1019. [Google Scholar]
  22. Li, Z. Description, Identification and Distribution Evaluation of Fractures in Deep Tight Sandstone Gas Reservoirs: A Case Study of Xuer Gas Reservoir in Xinchang Gas Field; Chengdu University of Technology: Chengdu, China, 2014. [Google Scholar]
  23. Hou, B.; Diao, C.; Li, D. An experimental investigation of geomechanical properties of deep tight gas reservoirs. J. Nat. Gas Sci. Eng. 2017, 47, 22–33. [Google Scholar] [CrossRef]
  24. Tan, P.; Jin, Y.; Han, K.; Hou, B.; Chen, M.; Guo, X.; Gao, J. Analysis of hydraulic fracture initiation and vertical propagation behavior in laminated shale formation. Fuel 2017, 206, 482–493. [Google Scholar] [CrossRef]
  25. Li, G.; Tang, C.A.; Li, L.C. Three-Dimensional Micro Flow-Stress-Damage (FSD) Model And Application In Hydraulic Fracturing In Brittle And Heterogeneous Rocks. Key Eng. Mater. 2011, 452–453, 581–584. [Google Scholar] [CrossRef]
  26. Wang, Y.; Gong, B.; Tang, C. Numerical simulation of fracture mechanism of columnar jointed basalts with different shapes under compression. Mech. Eng. 2022, 44, 88–99. [Google Scholar] [CrossRef]
  27. Xin, L.; Meng, S.; Liang, Z.; Tang, C.; Tao, J.; Tang, J. Microscale Crack Propagation in Shale Samples Using Focused Ion Beam Scanning Electron Microscopy and Three-Dimensional Numerical Modeling. Pet. Sci. 2023, 20, 1488–1512. [Google Scholar]
  28. Lou, Y.; Wu, Z.; Sun, W.; Yin, S.; Wang, A.; Liu, H.; Zuo, Y. Study on failure models and fractal characteristics of shale under seepage-stress coupling. Energy Sci. Eng. 2020, 8, 1634–1649. [Google Scholar] [CrossRef]
  29. Tang, M.; Wu, Z.; Wang, A.; Zuo, Y.; Lou, Y.; Liu, H.; Sun, W. Study on the Microscopic Fracture Process and Acoustic Emission of Shale Based on Digital Image. Geofluids 2021, 2021, 8874918. [Google Scholar] [CrossRef]
  30. Wu, Z.; Tang, M.; Zuo, Y.; Lou, Y.; Wang, W.; Liu, H.; Sun, W. Acoustic emission-based numerical simulation of tectonic stress field for tectoclase prediction in shale reservoirs of the northern Guizhou area, China. Energy Geosci. 2021, 3, 436–443. [Google Scholar] [CrossRef]
  31. Liu, H.; Mukherjee, B.; Zuo, Y.; Roy, P.N.S.; Wu, Z.; Sun, W.; Zheng, L.; Pan, C.; Lin, J. Fractal dimension used as a proxy to understand the spatial distribution for carlin-type gold deposits. Ore Geol Rev. 2023, 158, 105534. [Google Scholar] [CrossRef]
  32. Shi, K.; Zuo, Y.; Sun, W.; Liu, H.; Lin, J.; Li, T.; Chen, B.; Zheng, L. Fracture Mechanical Properties of Three—Dimensional Irregular Sandstone Particles Based on Digital-Image. Geotech. Geol. Eng. 2022, 40, 4991–5005. [Google Scholar] [CrossRef]
  33. Wu, Z.; Song, H.; Li, L.; Zhou, Z.; Zuo, Y.; Sun, W.; Liu, H.; Lou, Y. Study on the Damage Evolution Process and Fractal of Quartz-Filled Shale under Thermal-Mechanical Coupling. Geofluids 2021, 2021, 8843120. [Google Scholar] [CrossRef]
  34. Wu, Z.; Lou, Y.; Yin, S.; Wang, A.; Liu, H.; Sun, W. Acoustic and fractal analyses of the mechanical properties and fracture modes of bedding-containing shale under different seepage pressures. Energy Sci. Eng. 2020, 8, 3638–3656. [Google Scholar] [CrossRef]
  35. Zhang, P.; Lu, S.; Li, J.; Xue, H.; Wang, S. Quantitative characterization of microscopic pore structure for shales using scanning electron microscopy. J. China Univ. Pet. (Ed. Nat. Sci.) 2018, 42, 19–28. [Google Scholar]
  36. Cai, J.; Wei, W.; Hu, X.; Liu, R.; Wang, J. Fractal characterization of dynamic fracture network extension in porous media. Fractals 2017, 25, 1750023. [Google Scholar] [CrossRef]
  37. Li, J.; Zhao, W.; Liu, S. Description of Shale Pore Fracture Structure Based on Multi-fractal Theory. Adv. Pet. Explor. Dev. 2019, 18, 1–9. [Google Scholar]
  38. Lehamel, M.; Hammouche, K. Comparative Evaluation of Various Fractal Dimension Estimation Methods. In Proceedings of the 2012 Eighth International Conference on Signal Image Technology and Internet Based Systems, Sorrento, Italy, 25–29 November 2012; IEEE: Piscataway, NJ, USA, 2012; pp. 259–266. [Google Scholar]
  39. Panigrahy, A.M.N.K. Quantitative texture measurement of gray-scale images: Fractal dimension using an improved differential box counting method. Measurement 2019, 147, 106859. [Google Scholar] [CrossRef]
  40. Li, X.; Wei, W.; Wang, L. A new method for evaluating the pore structure complexity of digital rocks based on the relative value of fractal dimension. Mar. Pet. Geol. 2022, 141, 105694. [Google Scholar] [CrossRef]
  41. Vishal, V.; Chandra, D. Mechanical response and strain localization in coal under uniaxial loading, using digital volume correlation on X-ray tomography images. Int. J. Rock Mech. Min. Sci. 2022, 154, 105103. [Google Scholar] [CrossRef]
  42. Li, B.; Yu, S.; Yang, L.; Zhu, W.; Xue, Y.; Feng, D.; Wang, C.; Chen, Y. Multiscale fracture characteristics and failure mechanism quantification method of cracked rock under true triaxial compression. Eng. Fract. Mech. 2022, 262, 108257. [Google Scholar] [CrossRef]
  43. Ai, T.; Zhang, R.; Zhou, H.W.; Pei, J.L. Box-counting methods to directly estimate the fractal dimension of a rock surface. Appl. Surf. Sci. 2014, 314, 610–621. [Google Scholar] [CrossRef]
  44. Zhou, H.W.; Xue, D.J.; Jiang, D.Y. On fractal dimension of a fracture surface by volume covering method. Surf. Rev. Lett. 2014, 21, 212–481. [Google Scholar] [CrossRef]
  45. Sui, L.; Ju, Y.; Yang, Y.; Yang, Y.; Li, A. A quantification method for shale fracability based on analytic hierarchy process. Energy 2016, 115, 637–645. [Google Scholar] [CrossRef]
  46. Liang, X.; Hou, P.; Xue, Y.; Yang, X.; Gao, F.; Liu, J.; Taylor, T. A fractal perspective on fracture initiation and propagation of reservoir rocks under water and nitrogen fracturing. Fractals 2021, 29, 2150189. [Google Scholar] [CrossRef]
Figure 1. Surface geological map of the Zheng’an block.
Figure 1. Surface geological map of the Zheng’an block.
Fractalfract 08 00142 g001
Figure 2. Elastic–brittle damage ontological relationships of model cells under uniaxial compression.
Figure 2. Elastic–brittle damage ontological relationships of model cells under uniaxial compression.
Fractalfract 08 00142 g002
Figure 3. CT scanning and digital image characterization process. (a) X-ray 3D microscope-equipped nanoVoxel-4000 system; (b) CT scanning schematic; (c) CT scanning frequency; (d) 2D CT slice image; (e) grayscale value scanning map; (f) constructed laminated shale model.
Figure 3. CT scanning and digital image characterization process. (a) X-ray 3D microscope-equipped nanoVoxel-4000 system; (b) CT scanning schematic; (c) CT scanning frequency; (d) 2D CT slice image; (e) grayscale value scanning map; (f) constructed laminated shale model.
Fractalfract 08 00142 g003
Figure 4. Schematic of the stress loading method.
Figure 4. Schematic of the stress loading method.
Fractalfract 08 00142 g004
Figure 5. Stress–strain curve analysis.
Figure 5. Stress–strain curve analysis.
Fractalfract 08 00142 g005
Figure 6. Mechanical behavior of shale during damage—an example of Group 4. (a) Elastic modulus damage map. (b) Maximum principal stress cloud map. (c) Y-axis displacement vector map. (d) Acoustic emission evolution map.
Figure 6. Mechanical behavior of shale during damage—an example of Group 4. (a) Elastic modulus damage map. (b) Maximum principal stress cloud map. (c) Y-axis displacement vector map. (d) Acoustic emission evolution map.
Fractalfract 08 00142 g006
Figure 7. Three-dimensional fracturing image and meshing for fractal dimension calculation. (a) Three-dimensional crack images. (b) Schematic of the cubic coverage method for network segmentation.
Figure 7. Three-dimensional fracturing image and meshing for fractal dimension calculation. (a) Three-dimensional crack images. (b) Schematic of the cubic coverage method for network segmentation.
Fractalfract 08 00142 g007
Figure 8. Three-dimensional crack characterization process (a), and crack structure in slices (bd) at different locations.
Figure 8. Three-dimensional crack characterization process (a), and crack structure in slices (bd) at different locations.
Fractalfract 08 00142 g008
Figure 9. Horizontal stress ratio versus fractal dimension.
Figure 9. Horizontal stress ratio versus fractal dimension.
Fractalfract 08 00142 g009
Figure 10. Cloud map of the final damage displacement field and elastic modulus damage map of the shale specimen.
Figure 10. Cloud map of the final damage displacement field and elastic modulus damage map of the shale specimen.
Fractalfract 08 00142 g010
Table 1. Material parameters of rock samples [33,34].
Table 1. Material parameters of rock samples [33,34].
MaterialElastic Modulus/GPaCompressive Strength/MPaPoisson’s RatioCompression-Tension RatioInternal Friction Angle/(°)
Black shale51.61450.221435
Bedding30.91160.311330
Quartz96.03750.081560
Table 2. Surrounding pressure settings.
Table 2. Surrounding pressure settings.
MaterialGroup 1Group 2Group 3Group 4Group 5
Maximum horizontal stress (MPa)45.8464.9671.3980.8454.66
Minimum horizontal stress (MPa)36.8247.9555.8263.0841.07
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Lan, B.; Yang, R.; Wu, Z.; Jiang, H.; Li, X. Characterization of the Fine-Scale Evolution of Damage in Shale under the Influence of Two-Way Stress Differences Based on CT Images and Fractal Theory—The Example of the Anba Dyke in the Wufeng–Longmaxi Formation. Fractal Fract. 2024, 8, 142. https://doi.org/10.3390/fractalfract8030142

AMA Style

Lan B, Yang R, Wu Z, Jiang H, Li X. Characterization of the Fine-Scale Evolution of Damage in Shale under the Influence of Two-Way Stress Differences Based on CT Images and Fractal Theory—The Example of the Anba Dyke in the Wufeng–Longmaxi Formation. Fractal and Fractional. 2024; 8(3):142. https://doi.org/10.3390/fractalfract8030142

Chicago/Turabian Style

Lan, Baofeng, Ruidong Yang, Zhonghu Wu, Haishen Jiang, and Xinzheng Li. 2024. "Characterization of the Fine-Scale Evolution of Damage in Shale under the Influence of Two-Way Stress Differences Based on CT Images and Fractal Theory—The Example of the Anba Dyke in the Wufeng–Longmaxi Formation" Fractal and Fractional 8, no. 3: 142. https://doi.org/10.3390/fractalfract8030142

APA Style

Lan, B., Yang, R., Wu, Z., Jiang, H., & Li, X. (2024). Characterization of the Fine-Scale Evolution of Damage in Shale under the Influence of Two-Way Stress Differences Based on CT Images and Fractal Theory—The Example of the Anba Dyke in the Wufeng–Longmaxi Formation. Fractal and Fractional, 8(3), 142. https://doi.org/10.3390/fractalfract8030142

Article Metrics

Back to TopTop