Study on Voids and Seepage Characteristics within Rock Fracture after Shear Dislocation Viewing from CT Test and Numerical Modeling

: In rock mass engineering, stress balance changes often cause the relative slip of fractures along a wall surface, impacting the seepage behavior of ﬂ uid in the fractures. Using computer tomography (CT) scanning, spatial models of fractures with dislocations ranging from 0 to 10 mm were created to explore the relationship between changes in fracture dislocation and changes in ﬂ uid ﬂ ow behavior, respectively. The spatial fractal dimension of cavity distribution within the fractures was calculated using a thin-plate ﬁ lling approach to characterize the complexity of the fracture cavity distribution. The ﬂ uid ﬂ ow within the dislocation fractures was then simulated using COM-SOL


Introduction
Numerous fractures of different scales exist in subsurface rock mass, which controls engineering rock mass' physical, hydrological, and mechanical properties, thus affecting the deformation of surrounding rock, the internal stress distribution, and the seepage characteristics of the surrounding rock mass.Stress balance changes during rock engineering often result in relative sliding along fracture surfaces, altering the fluid permeation behavior within the fractures.For instance, in deep rock engineering applications such as enhanced geothermal system development [1], unconventional oil and gas reservoir exploitation [2], and carbon dioxide sequestration [3], the injection of high-pressure fluids into deep rock mass is employed to expand interconnected fracture pathways, which ultimately constructs a fracture network to enhance energy recovery or increase the storage capacity of deep rock strata.The injection of high-pressure fluids induces changes in the stress field in the rock mass, leading to relative sliding along the fracture surfaces and consequently impacting the flow characteristics.The complexity of fluid flow within the fractures arises from alterations in the stress environment caused by high-pressure fluid injection.As a result, fluid flow within the fractures transitions from linear to nonlinear states, affecting the solute and heat transfer properties within the fractures [4].Therefore, investigating the changes in permeability characteristics of fluid flow in rock fractures resulting from slip displacement is important.
Studying the deformation and permeability characteristics of a single fracture is a prerequisite for understanding the permeability properties of fractured rock masses [5,6].Early research on rock fractures assumed that the upper and lower surfaces of fractures were two smooth parallel plates.A cubic law was derived by simplifying the solution to the Navier-Stokes equation, indicating that the fracture's volume flow rate is proportional to the cube of the pressure gradient and the fracture aperture [7].However, natural fractures are not simple parallel plate structures; they possess a complex rough structure, heterogeneous void distribution, and tortuous geometric shapes.These factors lead to deviations from the cubic law in permeability [8,9], especially when fractures undergo deformation and contact influenced by various factors.This prompts changes in fluid flow velocity or direction within the channels, causing nonlinear fluid flow behaviors such as local eddy currents and turbulence [10].In contrast to the cubic law, the Forchheimer law describes the nonlinear relationship between pressure gradient and volume flow rate and has been verified in various contexts [11][12][13].Rong et al. [14] carried out non-matching flow tests on fractures and found that the Forchheimer equation is highly sensitive to changes in fracture aperture.Zhou et al. [15] proposed a quantitative analysis of the global inertial permeability for the nonlinear term in the Forchheimer equation, characterizing non-Darcy flow in porous media.Therefore, the Forchheimer equation provides an accurate model for describing the nonlinear seepage of mismatched fractures and better understanding the seepage characteristics of fractures under complex conditions.
For a single fracture, its void distribution exhibits complex characteristics closely related to the roughness structure of the fracture surface [16][17][18].Scholars usually parameterize the roughness structure of fracture surfaces and analyze its influence on the permeability characteristics of fractures.Parameters include joint roughness coefficients and plane fractal dimensions.By wavelet analysis, Zou et al. [19] divided the fracture surface roughness into primary and secondary roughness.They found that the secondary roughness significantly influences the nonlinear behavior of seepage.Zhang et al. [20] utilized 3D printing technology to create precise samples and conducted seepage tests and numerical simulations.They found that fracture roughness and aperture distribution affect the seepage and pressure fields.On the other hand, characterizing the cavity structure of fractures is also a research focus.Chen et al. [21], through improvements to the cubic covering method, filled the fracture space with ultra-thin square plates, describing the anisotropy of fracture space based on the difference in fractal geometric dimensions between the fracture space and parallel plates with equivalent apertures.However, when fractures undergo shear or slip dislocation under complex stress, the void distribution of fractures undergoes intricate changes with dislocation distance, and different stages have varying impacts on seepage [22,23].The irregular and unpredictable void distribution of fracture spaces poses a challenge to predicting permeability [10,24].Researchers have conducted numerous numerical simulations and experimental analyses to explore the seepage characteristics of non-matching fractures [25].Dang et al. [26] concluded that asperities' influence on aperture changes during the shear process can not be ignored, significantly affecting the seepage characteristics of the fluid in sheared fractures.Xiao et al. [27] summarized coupled shear seepage experiments on fractures, finding that a relatively low-seepage water pressure in the existing test system has a non-negligible impact on the understanding of seepage behavior in a high-stress environment.However, the challenge of maintaining specimen sealing under fracture dislocation conditions in pressure chambers with large shear displacements remains to be addressed and remains a challenge for current research development.
Most current experimental techniques face challenges in directly observing the details of fluid flow in actual rock fractures, especially in profoundly understanding the relationship between local geometric characteristics and seepage behavior in rock fractures [28].This includes non-linear seepage behaviors such as groove flow [19], eddy current [29], and reflux [30].Research on observing seepage behavior in fractures typically falls into two methods: destructive observation and non-destructive observation.Among nondestructive observation techniques, computer tomography technology (CT) can accurately obtain data on the distribution of voids in fractures without causing artificial damage to the fracture surfaces [16,[31][32][33].Using CT technology, Watanabe et al. [34,35] analyzed the flow characteristics of two granite single fracture specimens under confining pressure conditions.They found differences in fracture porosity and permeability under stress influence, but image noise led to insufficient data points.Additionally, based on pore distribution data and numerical simulation methods, the study further investigated the permeability characteristics of actual fractures.This research conclusively confirms the dominant flow behavior within the fractures.
Most of the above studies on fracture seepage are based on rock fractures in the original fracture state.However, there are still insufficient studies on the seepage law of mismatched fractures formed by rock dislocations on both sides of fractures.By combining CT experiments and numerical simulation methods, a more in-depth understanding of the evolution of fluid seepage behavior within non-matching fractures can be achieved.In this study, CT technology was utilized to obtain data on the distribution of voids in single granite fractures at different dislocation distances (0-10 mm).The fracture cavity model corresponding to the dislocation state is established, and the numerical simulation of dislocation fracture seepage is carried out.The parameters are quantified based on the fracture surface morphology data obtained by laser scanning to explore the internal causes of fracture space changes.At the same time, the changes in cavity interconnection channels are investigated through the fracture space structure parameters.Combined with the numerical simulation of seepage velocity field changes, the purpose is to understand further the influence of fracture dislocation state on fluid seepage characteristics.

Sample Preparation
The Alashan Block is located in the western part of Inner Mongolia, China, and is part of the North China Craton.This region features extensive granite coverage, including Gallidonian granite, Carboniferous granite, and Permian granite.Geological surveys suggest that the Alashan Block may be a favorable candidate for the geological disposal of high-level radioactive waste.For this study, core samples of banded granite from below the surface to a depth of 500 m were collected in the region.The granite is composed of approximately 45% plagioclase, 18% orthoclase, 25% quartz, and about 12% biotite.The grain size distribution of biotite, plagioclase, quartz, and orthoclase is approximately 0.8-1.5 mm, 1.3-3.0mm, 1.5-4.0mm, and 2.0-8.0 mm, respectively.The average density of the rock is around 2.65 g/cm 3 , with average shear and compressional wave velocities of 3151.8 m/s and 5860.4 m/s, respectively.The banded granite core samples were processed into standard specimens with a diameter of 50 mm and a height of 100 mm.Subsequently, the Brazilian splitting method was employed to split the specimens along the axial direction, obtaining granite samples with rough fractures, as shown in Figure 1a.
This study focuses on the variation in void distribution caused by rough fractures along a single direction, generated by 0-10 mm shear dislocation, and its impact on seepage characteristics.To ensure the stability of the samples after dislocation during the CT experiments and reduce the artificial dislocation error, semi-cylindrical metal shims with thicknesses ranging from 1 to 10 mm were manufactured based on the geometric characteristics of the samples, as shown in Figure 1b.Two shims of each thickness were prepared and attached to both ends of the samples using glue, resulting in granite fracture samples under 0-10 mm dislocation states.

Fracture Surface Scanning
To further understand the relationship between the rough structure of the fracture surface and the change in the void when the fracture dislocation occurs, three-dimensional surface morphology data were initially obtained using a 3D laser scanner [36].The 3D laser scanner employed in this study is the Revopoint MINI, featuring a scanning accuracy of 50 µm and a scanning rate of 10 frames per second, as illustrated in Figure 2. The process of acquiring surface morphology data is outlined below: Before scanning, the 3D laser scanner underwent calibration and adjustment to ensure proper height and angle, enabling the complete recording of the specimen's morphology.Once the scanner was fixed, the specimen was positioned on the tray at an appropriate location.A black accessory that the instrument could not recognize was used to stabilize the specimen, ensuring balance.The scanning program was initiated, where the tray drove the specimen's rotation, collecting point cloud data of the surface morphology of the fracture through scans at various angles.The comprehensive point cloud data of the fracture surface were imported into Geomagic Wrap v. 2021 software for post-processing tasks such as filtering, registration, and merging.This process aimed to create a precise three-dimensional model of the rough fracture surface.

CT Tests
Based on the prepared fracture samples, CT experiments were conducted to obtain cavity distribution data of 0-10 mm dislocated fracture samples.The CT scanning experiments relied on a real-time high-energy CT scanning system provided by the Institute of Geology and Geophysics, Chinese Academy of Sciences.This system comprises an X-ray source, detector, rotating platform, pressurization system, and computer, as illustrated in Figure 3.The samples under different shear dislocation conditions were subjected to slice scanning to obtain precise information about the three-dimensional fracture cavities without damaging the specimens, which facilitated the analysis of fracture spatial structure distribution.The test did not consider the effect of stress to minimize the influence of damage accumulation on the dislocation fracture surface under stress conditions on the exploration of nonlinear seepage laws.
The experimental steps of the CT scanning experiment are as follows: firstly, determining equipment parameters such as imaging range, scanning time, and radiation dose; secondly, fixing the samples in different shear dislocation states using elastic gasket and tape; placing the samples on the scanning bed, and adjusting position as needed to cover the observation area; subsequently, initiating the scanning program, with the X-ray source and detector starting to rotate, capturing a projection image every 0.4°, and recording a series of two-dimensional slice images with a slice thickness of 50 µm.Parameters such as scanning speed, overlap rate, interval, and sampling frequency could be adjusted during the scanning process.After completing the scan, preprocessing steps were applied to eliminate noise and artifacts, including image denoising, enhancement, and filtering using three-dimensional visualization and data analysis software (AVZIO v.

Fracture Cavity Characterization
The characteristics of rock fractures typically include the morphology of fracture surfaces and the diversity of void changes.Quantifying these characteristics can aid in analyzing the impact of fracture features on seepage.Parameters such as mechanical aperture, fractal dimension, roughness angle, and relative spatial heterogeneity fractal dimension can effectively describe the complexity of fracture surface roughness and the distribution of space void changes.This enables a better understanding of how the structural changes in fracture space, resulting from the uplift of contacted fracture surfaces due to dislocation, influence fluid flow behavior.
Based on the point cloud data obtained from three-dimensional laser scanning of the upper and lower surfaces of the fracture, high-precision grids were fitted with grid intervals of 0.1 mm in the X and Y directions.Therefore, the mechanical aperture of the fracture can be calculated using Equation (1) [37]: where ei represents the mechanical aperture at the i-th grid node along the fracture dislocation, and zi + and zi − denote the spatial heights of the upper and lower fracture surfaces at the corresponding grid node, respectively.Under stress conditions, fracture surfaces typically experience contact closure.In reality, the mechanical aperture of natural fractures cannot be less than zero, and under stress, fracture surfaces tend to come into contact and close.In this study, an aperture size smaller than 50 µm is considered closed.In this closed state, the contribution to fluid flow can be negligible and is treated as a stagnant fluid region [30].To provide a more intuitive description of the changes in fracture aperture with increasing fault distance, the average mechanical aperture (em) is calculated by taking the mean of all mechanical aperture data.
The aperture distribution of fractures reflects the uneven distribution of voids caused by the contact uplift of asperities on the fracture surface.To accurately assess the influence of asperities on fracture aperture, scholars have quantified and investigated the impact of fracture surface roughness on fracture permeability [38][39][40].The contact uplift of the fracture surface is primarily related to the friction of the asperities on the fracture surface.Therefore, the roughness angle α of the fracture surface is used to characterize the absolute tilt parameter of the asperities along the dislocation direction, exploring the reasons for the variation in the void distribution with dislocation.The fracture surface roughness angle along the single direction of dislocation can be determined by Equation ( 2) [41].
where zi and yi represent the coordinate values on the fracture surface at the corresponding grid node, and i is the characterization of the grid nodes of the single profile along the shear direction.The degree of tortuosity of the profile along the single direction is determined by the above equation, based on which the effect of the tortuosity angle of the fracture surface on the change of the fracture dislocation aperture is analyzed.The roughness angle αi of the fracture surface describes a single scalar value and is insufficient to characterize the unevenness of the contact and cavity formed by two rough fracture surfaces.Chen et al. [21] made improvements using the cube covering method, filling the fracture space with ultra-thin square plates.The anisotropy of the fracture space is described based on the difference in fractal dimensions between parallel plates with the same opening, as shown in Figure 5.The surface roughness of natural fractures typically follows a self-similar fractal distribution, and the fracture space exhibits self-similarity.The ultra-thin square plate filling method describes interconnected pore space's three-dimensional irregular cavity distribution.The spatial heterogeneity fractal dimension can be obtained by calculating the number of ultra-thin square plates at different side lengths.This parameter reflects the changes in spatial dimension, connectivity, and symmetry.This quantitative method can better understand the complexity and heterogeneity of the fracture space cavity distribution.Based on the above research, the application of spatial fractional dimension can be further expanded to describe the complexity of fracture space cavity distribution under non-matching conditions of fracture surfaces.The specific method is as follows: Based on the high-precision mesh fitting of the actual surface morphology mentioned before, ultra-thin square plates with a width of δ and a thickness of ξ are filled at the mesh nodes.Therefore, the number of thin square plates covering the fracture space at the mesh nodes can be determined by Equation (3).
where the different elevations of the four mesh nodes are represented by z(i,j) and + and − denote the upper and lower surfaces, respectively.Here, i and j represent the position indices along the width and length directions of the ultra-thin square plate, and INT(X) is the function that returns the largest integer not greater than X.Therefore, the total number of plates filling the entire interconnected fracture space at a given width can be calculated using Equation (4).
By continuously changing the size of the mesh nodes through Kriging interpolation, the corresponding adjustment of the ultra-thin square plate width δ is made to match the mesh nodes.The total number of ultra-thin square plates also changes accordingly.The relationship between the two can be plotted in a double-logarithmic coordinate system, and the slope of the curve is considered as −D, representing the calculated spatial fractal dimension, as shown in Equation (5).
where k is a proportionality constant, and D represents the corresponding spatial heterogeneous fractional dimension, which describes the heterogeneity of the distribution of interconnected cavity channels in the fracture space.Based on this, by comparing parallel plates with the same mechanical aperture, DΔ = D − D ref can be obtained.This represents the relative heterogeneous fractional dimension and is used to characterize the complexity of the fracture spatial void distribution.

Seepage Theoretical Model
The Navier-Stokes equation is generally used to describe fluid motion [38,39].The steady flow of incompressible fluids can be expressed as Equation (6).The N-S and continuity Equations ( 6) and ( 7) are widely used to analyze single fracture flow characteristics [42].
u 0   (7) where u represents the velocity vector in m/s, p is the pressure in Pa, ρ is the fluid density in kg/m 3 , and µ is the fluid's dynamic viscosity in Pa•s.
With the increase in shear dislocation, the distribution of voids in the fracture space continuously changes, and the fluid in the permeability process is more likely to enter a nonlinear flow state locally.The calculation of nonlinear fluid flow is commonly based on the Forchheimer empirical equation, as shown in Equation ( 8) [43,44].
where ▽P is the pressure gradient, A is the linear term coefficient or viscous term coefficient, and B is the nonlinear term coefficient or inertia coefficient.When the influence of fluid inertia can be neglected at low flow rates, the nonlinear term in the second part (BQ 2 ) disappears, and the Forchheimer law reduces to Darcy's law.The coefficients A and B are mainly related to the geometric characteristics of the fracture surface tomography [45].Combining Darcy's law, the permeability of the fracture is calculated through Equation ( 9) [46]: where K represents the permeability of the fracture, m 2 ; S is the cross-sectional area of the fracture, m 2 .The dimensionless Reynolds number (Re) is commonly used as an indicator for determining fluid flow regimes, defined as the ratio of inertial forces to viscous forces.For fluid within fractures, Re can be calculated using Equation (10) [22]: In the equation, bh represents the hydraulic aperture, and w is the width of the fracture along the direction of the pressure gradient.

Numerical Models
To investigate the changes in the seepage field of fractured rock masses with displacements ranging from 0 to 10 mm, the single-fracture seepage simulations were conducted based on the fracture cavity model morphology obtained from CT scans under different dislocation conditions.The simulations involved solving the N-S and continuity equations to get the seepage velocity field for dislocation ranging from 0 to 10 mm.The corresponding mesh models are presented in Figure 6.The analysis of cavity variations was combined with the simulation results to understand the causes of complex fluid seepage behavior.In the simulation process, the fluid was assumed to be an incompressible Newtonian fluid with constant density (1 × 10 3 kg/m 3 ) and dynamic viscosity (1 × 10⁻ 3 Pa•s).The simulation focused on single-phase, saturated, and steady-state conditions for the single-fracture seepage.Considering granite's low porosity and permeability, fluid exchange between the granite matrix and fractures was neglected [46].Previous experiments and observations have indicated that roughness distribution on the fracture surfaces leads to complex fluid seepage behavior [19].Solving the threedimensional N-S equation remains challenging.Under certain conditions, where the fluid is in laminar or creeping flow, i.e., Re < 2000, the fracture permeability can be estimated using Darcy's law and the Forchheimer equation [45,47].
COMSOL Multiphysics (https://www.comsol.com/) is a finite element-based multiphysics simulation software widely used for simulating fluid flow in rock fractures.Initially, the models of fractured cavities with various dislocation states are imported into the software.The Y-axis is set along the length of the fracture, the X-axis along the width, and the Z-axis represents the direction of the dislocation-induced uplift, as illustrated in Figure 7.The color filling indicates the local mechanical aperture size of the dislocationinduced fracture.Subsequently, a network is established for the dislocation-induced fracture space model.COMSOL's built-in mesh generator is employed to generate a discretized mesh automatically suitable for the solver.The mesh is optimized based on the pore size, and the optimal mesh size is selected based on the computation capacity of the workstation.Boundary conditions for the spatial model are then determined according to the physical area.Uniform inlet and outlet boundary conditions are applied to all dislocation fracture models, with the fluid injection flow rate set as the inlet condition (10 mL/s) and the outlet boundary condition set to 0 Pa.The fluid flows primarily along the direction of fracture length under the fluid pressure difference.Other boundaries are treated as no-slip boundary conditions, consistent with many laboratory experimental conditions [12,18].
This study uses COMSOL's laminar flow module to solve the N-S equations.When fluid seepage is simulated, the mesh conditions at the inlet and outlet boundaries are further optimized.A fine mesh with a size of 0.001 mm is used at the inlet and outlet boundaries, while a coarse mesh with a size of 0.1 mm is applied internally.The total number of tetrahedral elements is typically on the order of 10 6 .The seepage field for different dislocation states is exported, and the dimensionless Reynolds number Re for fracture flow under different dislocation states is calculated as lower than 300.Thus, the flow regime in the dislocation-induced fractures can be assumed to be linear flow.The permeability of dislocated fracture models can then be calculated by Equation ( 9).

Statistics of Fractures' Local Mechanical Aperture Based on CT Data
Figure 8 shows the distribution pattern of local mechanical aperture data with the fracture dislocation distance varying from 0 to 10 mm as a visualization of the variation in the distribution pattern with dislocation distance.The horizontal axis represents the mechanical aperture size, while the vertical axis of each dislocation reflects the amount of corresponding aperture data.It can be observed that with the increase in fracture dislocation distance, the distribution of local mechanical aperture data undergoes a significant evolution from a relatively concentrated distribution in a small aperture to a more extensive local aperture distribution.For instance, when the fracture dislocation is at a distance of 1 mm, the distribution of local mechanical aperture data is relatively narrow, displaying a predominant normal distribution.The local aperture values are primarily concentrated around 1~2 mm.As the fracture dislocation distance increases, the distribution of local aperture data significantly widens.When the fracture dislocation distance reaches 10 mm, most local aperture values are between 1.5 and 3.5 mm. Figure 9 presents the variation in the average mechanical aperture with fracture dislocation distance under normal stress conditions and stress-free conditions.It is evident that, with the increase in fracture dislocation distance, the average mechanical aperture exhibits nonlinear growth, with the growth rate gradually slowing down, following a logarithmic relationship.This behavior aligns with the experimental and numerical simulation results of Li et al. [47] under normal stress conditions, represented by the red and blue curves in Figure 9. Furthermore, it can be observed that the values of calculated average mechanical aperture em do not entirely align with the fitted curve, reflecting the influence of the actual fracture roughness on the average fracture aperture.Compared with the fitting results obtained by Zou et al. [48] through shear experiments, the average mechanical aperture under stress-free conditions is notably lower than the values under certain normal stress conditions.Considering the curve in the graph variations of the average mechanical aperture em with fracture dislocation distance, it can be observed that when dislocation occurs on the fracture surface, the change in fracture aperture generally maintains a certain level.As the fracture dislocation distance increases, the growth of the average mechanical aperture em gradually slows down.However, for certain fracture dislocation distances, such as 7 mm and 10 mm, the average mechanical aperture em deviates from the curve's growth trend.To better understand the variation in fracture apertures, the roughness angle α is calculated along the length of the rock core in Figure 10. Figure 10a presents a two-dimensional vertical height map of fracture sections along the radial direction of the core, obtained by stacking fracture surface segments taken every 2 mm in the horizontal direction.It shows that the undulation range of the fracture surface is about 4 mm, consistent with the trend of aperture changes when the fracture is uplifted.Figure 10b  Fracture displacement involves altering the void distribution within a fracture and affects the connectivity of void spaces.In this study, MATLAB was employed to code the ultra-thin square plate filling method from Section 2.3 to calculate the relative fractal dimension DΔ of the fracture void space for fracture displacements ranging from 1 to 10 mm.This was carried out to assess spatial complexity changes resulting from fracture displacement variations.The ultra-thin square plate filling method is closely related to the selected thickness of the plates.This research chose four thicknesses, namely 0.1 mm, 0.2 mm, 0.3 mm, and 0.4 mm [5].Subsequently, a log-log plot was established between the number of ultra-thin square plates and square plate width, and the curve equation was fitted.The coefficient of the linear term in the curve fitting equation represents the calculated spatial fractal dimension.The spatial fractal dimension D was then subtracted from the fractal dimension of void space with the same mechanical aperture, and the relative spatial fractal dimension was obtained.The relationship between the relative spatial fractal dimension DΔ and fracture dislocation distance is shown in Figure 11, where some of the fracture cavity complex variations are labeled in the figure.At the beginning of the curve, the distribution of fracture voids approaches that of parallel plates with the same mechanical aperture, indicating a reduction in void connectivity and a decrease in fracture complexity.As fracture displacement changes, it can be observed that the relative spatial fractal dimension of the fracture exhibits nonlinear growth.Specifically, at 3 mm and 7 mm fracture displacements, the relative spatial fractal dimension DΔ deviates from the curve, smaller than the predicted relative spatial fractal dimensions.Figure 6 shows that the 3 mm dislocation fracture model has more contacts and widespread aperture distribution, resulting in a relatively complex fracture void distribution.Regarding the spatial fractal dimension variation in the model of 7 mm dislocation, combined with the 7 mm fracture void model in Figure 6 and the average mechanical aperture changes in Figure 7, it is inferred that the decrease in DΔ is due to a local reduction in fracture aperture, leading to an increase in void connectivity and an overall increase in complexity.

Seepage Velocity Field Analysis of Dislocation Fracture
Figure 12 presents the curve of the fracture permeability K as a function of fracture displacement d.The graph shows that with the increase in fracture displacement d, the fracture permeability K shows a nonlinear growth trend, and the growth trend slows down.This trend is consistent with the results obtained by Li et al. [47] in experiments and simulations of fracture shear seepage under normal stress conditions.Under stress conditions, fracture shear mainly occurs in the compaction stage, and the fracture aperture changes relatively little with fracture displacement, resulting in nonlinear growth in permeability K, which eventually tends to flatten.When fractures experience fracture displacement under stress-free conditions, the average mechanical aperture of the fractures is more significant than that under stress conditions.Similarly, the change in permeability K is more important than that under stress conditions, especially after the fracture displacement exceeds 5 mm, where the permeability is twice that of the shear permeability under stress conditions.Figure 12 shows that when the fracture displacement is 7 mm, the fracture permeability K deviates from the curve, and the permeability is lower than that under previous fracture displacement states.This is closely related to the increase in fracture complexity mentioned in Section 4.1.Further analysis of the fracture flow velocity field is conducted to understand better the relationship between permeability K and the void space variation caused by fracture displacement.Table 1 summarizes the parameters analyzed for the flow velocity field.The occurrence of nonlinear phenomena is closely related to the contact state and shape of the fracture surfaces.For example, when the fracture displacement is 3 mm, the contact area of the fracture surface is more significant than that at a 2 mm fracture displacement and exhibits an irregular distribution.The increase in the relative spatial fractal dimension DΔ leads to fluid accumulation and pressure at smaller apertures, forming vortices after contact points.This is consistent with the understanding of nonlinear behavior in flow presented by Lee et al. [10].The right-upward causes the vortex at 40 mm along the Y-axis inclined geometric structure with a 10° inclination, and it is influenced by the boundary shape, resulting in a flow direction change and the formation of a vortex.Regarding the deviation of the permeability K from the growth trend in Figure 11 when the fracture displacement d is 7 mm, in addition to the increased complexity indicated by the relative spatial fractal dimension DΔ, the development area of the vortex is much larger than that in the flow velocity field at other fracture displacements.This is because at this fracture displacement, the larger aperture at 40 mm along the Y-axis occupies a larger area.Considering the variation in the fracture flow velocity field with fracture displacement d in Figure 13, it can be seen that the influence range of the vortex gradually expands with increasing fracture displacement, and it is affected by the boundary shape and fracture surface structure, impacting the development of permeability K.Although the Reynolds number Re is lower than 300 in Table 1, the fluid velocity fields in Figure 13 show a century non-linear behavior.To better understand the nonlinear changes in the process of fracture seepage, the nonlinear factor E is introduced as in [49]: This parameter quantifies the proportion of nonlinear behavior in the overall seepage process, providing a comprehensive perspective for understanding the complexity of fracture seepage.When the nonlinear factor E is set to 0.1, the fluid transitions from linear to nonlinear states [50], and the corresponding critical hydraulic gradient is defined as the critical hydraulic gradient Jc [51].The Forchheimer equation is employed in describing the nonlinear flow characteristics of fracture seepage.The parameters of this equation have distinct meanings, where A represents the intrinsic permeability of the fracture, and B is considered a crucial parameter for describing nonlinear flow [15,18,45].Therefore, an analysis of the variations in the critical hydraulic gradient Jc and the inertial coefficient B with the fracture displacement is presented in Figure 14.
Figure 14 shows that with the increase in fracture displacement d, the critical hydraulic gradient Jc and the inertial coefficient B exhibit a gradually decreasing power-law relationship.The rate of decrease gradually slows down, and the nonlinear factor E under different fracture displacements is consistently greater than 0.1.This indicates that fracture seepage predominantly manifests as nonlinear seepage in this study.Figure 14b illustrates the variation in the inertial parameter B of the Forchheimer equation with fracture displacement d.The curve's form aligns with the results of other researchers and is smaller than the seepage under shear conditions when the fracture experiences stress.This is because the average aperture change when the fracture undergoes displacement is more significant than the aperture change under shear conditions with stress, making fracture seepage less prone to becoming nonlinear, manifested in the curve as a reduction in the inertial coefficient B. As the fracture displacement d increases, the fluid becomes less inclined to exhibit nonlinearity.It is noteworthy that, as indicated by the flow field distribution, even in cases of relatively small nonlinear coefficients, complex fluid behaviors such as vortices may still occur during the seepage process [10], as shown in Figure 13d.This suggests that the nonlinear coefficient does not necessarily preclude the occurrence of other complex flow phenomena during the seepage process; instead, it is used to quantify the proportion of nonlinear phenomena in the overall seepage.Therefore, considering both the nonlinear coefficient and the velocity field distribution comprehensively enhances our understanding of the complexity of fracture seepage.

Conclusions
In this study, CT scans were employed to investigate the void distribution of actual fractures under different fracture dislocation distance states.Spatial variations of void during the fracture dislocation distance process were analyzed based on parameters such as average mechanical aperture, relative spatial fractal dimension, and fracture surface roughness angle distribution.Subsequently, numerical fluid flow simulations were conducted based on the fracture void models obtained from CT experiments at different fracture dislocation distances.The impact of fracture dislocation on fracture permeability and the nonlinear behavior of fluid flow was analyzed.The main conclusions are as follows: In summary, this study reveals the influence of fracture dislocation on fluid flow characteristics through a combination of CT experiments and numerical simulations.This research holds crucial significance for understanding the nonlinear behavior of fracture flow and the impact of fracture dislocation on permeability, providing a more profound perspective for fracture flow studies.Further studies are still needed on the effects of different fracture morphologies and the seepage-flow coupling of fracture shear processes.

Figure 1 .
Figure 1.Granite fracture and sample end gasket: (a) granite samples with a single fracture surface; (b) diagram of gasket for controlling dislocation.

Figure 2 .
Figure 2. Fracture surface morphology reconstruction by three-dimensional laser scanning.
2019).Different shear dislocation states fracture cavity spatial models were established, as shown in Figure 4. Subsequently, post-processing was performed based on the established three-dimensional cavity spatial models to enhance model quality.

Figure 3 .
Figure 3.In situ loading high-energy CT real-time scanning system (in Institute of Geology and Geophysics, Chinese Academy of Sciences).

Figure 4 .
Figure 4. CT reconstruction of rock fracture space with different shear dislocations.

Figure 5 .
Figure 5. Ultra-thin square plate filling method for the spatial fractional dimension of the void within a fracture.

Figure 6 .
Figure 6.Grid models for the fractures at different dislocations from 1 mm to 10 mm.

Figure 7 .
Figure 7.The simulation model of seepage in dislocation fracture, where the color fills in the figure characterize the magnitude of local aperture.

Figure 8 .
Figure 8. Statistical analysis of fractures' local mechanical aperture.

Figure 9 .
Figure 9.The relationship between the average mechanical aperture of the fracture and the dislocation distance [47].

Figure 10 .
Figure 10.Roughness characteristics of a fracture surface: (a) stacking diagram of fracture surface topography in the direction of fracture length; (b) absolute roughness angle distribution of fracture surface along the dislocation direction.

Figure 11 .
Figure 11.The relationship between fracture spatial fractal dimension and dislocation distance.(The red arrow leads to a comparison of fracture openings).

Figure 13
Figure 13 illustrates the distribution of the flow velocity field in fractures at various fracture displacements, with 400 streamlines set at the inlet.A grayscale map of the local aperture distribution is filled on the lower surface to provide a more intuitive understanding of the influence of fracture displacement on fluid flow.Figure 13 shows that the flow velocity within the fracture gradually decreases with increasing fracture displacement and exhibits a distinct preferential flow phenomenon, with streamlines avoiding areas of smaller local apertures.Additionally, as the fault displacement d gradually increases from 1 mm to 4 mm, the contact between fracture surfaces along the Y-axis at approximately 70 mm gradually occurs, increasing the contact area.When the fracture displacement d increases from 5 mm to 10 mm, the contact region along the Y-axis at 70 mm becomes a lower aperture, and the aperture distribution becomes complex and variable.Considering the flow velocity field at various fracture displacements, it is observed that the mechanical apertures are more extensive along the Y-axis at 40 mm compared to other regions.When the fracture displacement d increases from 3 mm to 10 mm, vortices and other nonlinear phenomena occur in the fracture at approximately 40 mm and 70 mm along the Y-axis.The occurrence of nonlinear phenomena is closely related to the contact state and shape of the fracture surfaces.For example, when the fracture displacement is 3 mm, the contact area of the fracture surface is more significant than that at a 2 mm fracture displacement and exhibits an irregular distribution.The increase in the relative spatial fractal dimension DΔ leads to fluid accumulation and pressure at smaller apertures, forming vortices after contact points.This is consistent with the understanding of nonlinear behavior in flow presented by Lee et al.[10].The right-upward causes the vortex at 40 mm along the Y-axis inclined geometric structure with a 10° inclination, and it is influenced by the boundary shape, resulting in a flow direction change and the formation of a vortex.Regarding the deviation of the permeability K from the growth trend in Figure11when the fracture displacement d is 7 mm, in addition to the increased complexity indicated by the relative spatial fractal dimension DΔ, the development area of the vortex is much larger than that in the flow velocity field at other fracture displacements.This is because at this fracture

Figure 13
Figure 13 illustrates the distribution of the flow velocity field in fractures at various fracture displacements, with 400 streamlines set at the inlet.A grayscale map of the local aperture distribution is filled on the lower surface to provide a more intuitive understanding of the influence of fracture displacement on fluid flow.Figure 13 shows that the flow velocity within the fracture gradually decreases with increasing fracture displacement and exhibits a distinct preferential flow phenomenon, with streamlines avoiding areas of smaller local apertures.Additionally, as the fault displacement d gradually increases from 1 mm to 4 mm, the contact between fracture surfaces along the Y-axis at approximately 70 mm gradually occurs, increasing the contact area.When the fracture displacement d increases from 5 mm to 10 mm, the contact region along the Y-axis at 70 mm becomes a lower aperture, and the aperture distribution becomes complex and variable.Considering the flow velocity field at various fracture displacements, it is observed that the mechanical apertures are more extensive along the Y-axis at 40 mm compared to other regions.When the fracture displacement d increases from 3 mm to 10 mm, vortices and other nonlinear phenomena occur in the fracture at approximately 40 mm and 70 mm along the Y-axis.The occurrence of nonlinear phenomena is closely related to the contact state and shape of the fracture surfaces.For example, when the fracture displacement is 3 mm, the contact area of the fracture surface is more significant than that at a 2 mm fracture displacement and exhibits an irregular distribution.The increase in the relative spatial fractal dimension DΔ leads to fluid accumulation and pressure at smaller apertures, forming vortices after contact points.This is consistent with the understanding of nonlinear behavior in flow presented by Lee et al.[10].The right-upward causes the vortex at 40 mm along the Y-axis inclined geometric structure with a 10° inclination, and it is influenced by the boundary shape, resulting in a flow direction change and the formation of a vortex.Regarding the deviation of the permeability K from the growth trend in Figure11when the fracture displacement d is 7 mm, in addition to the increased complexity indicated by the relative spatial fractal dimension DΔ, the development area of the vortex is much larger than that in the flow velocity field at other fracture displacements.This is because at this fracture

Figure 13 .
Figure 13.Distribution of flow velocity field under different dislocation distances of fractures.

Figure 14 .
Figure 14.The relationship between critical hydraulic gradient and inertia coefficient with fracture dislocation distance: (a) variation with the critical hydraulic gradient of dislocation distance; (b) variation with the inertial coefficient B of dislocation distance [47].

( 1 )
Fracture dislocation significantly influences aperture distribution.The aperture frequency distribution obtained from CT experiments exhibits a normal distribution.As the fracture dislocation distance increases, the numerical range of the local aperture gradually transitions from relatively narrow to broader, and the maximum local aperture of fractures increases.The average mechanical aperture (em) growth follows a logarithmic function and is closely related to the roughness angle of the fracture surface.(2) The average mechanical aperture and the distribution of local apertures during fracture dislocation affect the complexity of void connectivity.With increasing fracture dislocation distance, the relative spatial fractal dimension (DΔ) of fractures shows a logarithmic relationship, indicating a decrease in void connectivity and a reduction in the complexity of the fluid flow pathways.(3) The fracture permeability increases nonlinearly with the increase in dislocation distance.Under the same injection flow rate, the seepage flow in the fracture changes from linear flow to nonlinear flow, and the dominant flow linear flow in dislocation fracture is very obvious.When the dislocation distance reaches 5 mm, nonlinear behaviors such as eddy currents occur, and the influence range of eddy currents gradually expands with the increase in dislocation distance.The fracture surface morphology structure as the boundary of seepage has a significant influence on the seepage behavior.It is also found that the change in dislocation fracture permeability with dislocation distance is closely related to the change in spatial fractal dimension.(4) With the increase in dislocation distance, the critical hydraulic gradient Jc and the inertia coefficient B, which distinguish linear and nonlinear flows, show a power function decrease and gradually tend to a stable value, indicating that the fluid in the fracture is more likely to produce nonlinear flow with the increase in dislocation distance.The degree of distribution of the cavity is gradually reduced, and the fluid presents a dominant flow, which enhances the nonlinear ability of fluid in the fracture.The inertia coefficient B and the critical hydraulic gradient Jc in the Forchheimer equation can quantitatively describe the nonlinear seepage characteristics of the fracture.

Table 1 .
Seepage velocity field parameters for different dislocation fractures.