Next Article in Journal
Digitalized Automation Engineering of Industry 4.0 Production Systems and Their Tight Cooperation with Digital Twins
Next Article in Special Issue
Experimental Investigation of the Growth Law of Multi-Fracture during Temporary Plugging Fracturing within a Stage of Multi-Cluster in a Horizontal Well
Previous Article in Journal
Data-Driven Models for Forecasting Failure Modes in Oil and Gas Pipes
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Characteristics Description of Shale Fracture Surface Morphology: A Case Study of Shale Samples from Barnett Shale

1
School of Petroleum Engineering, Xi’an Shiyou University, Xi’an 710065, China
2
Key Laboratory of Unconventional Oil & Gas Development, China University of Petroleum (East China), Ministry of Education, Qingdao 266580, China
3
No.8 Oil Production Plant, Changqing Oilfield, China National Petroleum Corporation (CNPC), Xi’an 710018, China
4
No.2 Oil Production Plant, Changqing Oilfield, China National Petroleum Corporation (CNPC), Qingyang 745100, China
5
Drilling and Production Technology Research Institute, Liaohe Oilfield, China National Petroleum Corporation (CNPC), Panjin 124000, China
6
Oil Production Technology Research Institute, Dagang Oilfield, China National Petroleum Corporation (CNPC), Tianjin 300451, China
7
Directional Well Technology Service Company, Bohai Drilling Engineering Company Limited, China National Petroleum Corporation (CNPC), Tianjin 300457, China
8
School of Aeronautics, Northwestern Polytechnical University, Xi’an 710072, China
*
Author to whom correspondence should be addressed.
Processes 2022, 10(2), 401; https://doi.org/10.3390/pr10020401
Submission received: 14 January 2022 / Revised: 4 February 2022 / Accepted: 16 February 2022 / Published: 18 February 2022

Abstract

:
Shale reservoirs are the hot issue in unconventional resources. The key to the development of shale reservoirs lies in the complex fractures, which are the only path for fluid to migrate from the matrix to the wellbore in shale reservoirs. Therefore, the characteristics of shale fracture surface morphology directly affect fluid migration in shale reservoirs. However, there are few reports about the characteristics of shale fracture surface morphology as the parallel plate model was commonly used to characterize the fracture, neglecting its surface morphology characteristics and leading to great deviation. Thus, description methods were introduced to characterize shale fracture surface morphology with the aim to provide a foundation for the development of shale resources. Three shale samples were fractured by the Brazilian test, and the height distribution of the fracture surface was captured by a three-dimensional profilometer. Then, three-dimensional fracture surface morphology was regarded as a set of two-dimensional profiles, which converted three-dimensional information into two-dimensional data. Roughness, joint roughness coefficient, fractal dimension, tortuosity, and dip angle were employed to characterize shale fracture surface morphology, and their calculation methods were also, respectively, proposed. It was found that roughness, joint roughness coefficient, fractal dimension, tortuosity, and dip angle were all directional, and they varied greatly along with different directions. Roughness, joint roughness coefficient, fractal dimension, tortuosity, absolute dip angle, and overall trend dip angle were among 0.0834–0.2427 mm, 2.5715–10.9368, 2.1000–2.1364, 1.0732–1.1879, 17.7498°–24.5941°, and −3.7223°–13.3042°, respectively. Joint roughness coefficient, fractal dimension, tortuosity, and dip angle were all positively correlated with roughness.

1. Introduction

With the increasing depletion of hydrocarbon resources, unconventional resources such as shale, tight sandstone, and other reservoirs, which were previously considered to be undevelopable, have become hot issues in recent decades [1,2,3]. Compared to conventional low permeability reservoirs, shale reservoirs possess extra-low permeability, usually among nanoscale range, making the development methods for conventional low permeability reservoirs unsuitable for shale reservoirs development [4,5,6]. In recent years, a variety of methods were attempted for the development of shale reservoirs [7,8,9,10]; however, there is still no final conclusion about which method is more suitable. However, there is a consensus that although the shale reservoirs hav extremely low permeability, they posess the common characteristics of developed bedding and natural fracture as well as a high content of brittle minerals. Thus, large-scale volume fracturing technology is usually recommended to transform shale reservoirs [11], forming a volumetric fracture network structure in three-dimensional space, which could communicate the matrix of shale reservoirs and the wellbore, realizing the effective migration of hydrocarbon from the matrix to the wellbore. Therefore, the volumetric fracture network structure is the key to the effective development of shale reservoirs. The volumetric fracture network structure is composed of crisscrossed natural fractures, induced fractures, and hydraulic fractures. According to the fracture morphology, the aforementioned three kinds of fractures can be divided into tensile fractures and shear-slipping fractures [11]. Regardless of the nature and morphology of fractures, the volumetric fracture network structure is essentially composed of fractures, although the fractures are of different shapes.
As the only path for fluid migration from the reservoir matrix to the wellbore, the fracture that directly contacts with reservoir fluid has a binding effect on the migration of fluid [12,13,14]. The characteristics of fracture surface morphology have a direct impact on fluid migration, which requires an accurate and reasonable evaluation of fracture surface morphology. However, in the early studies of fluid migration in fractures, the influence of fracture surface morphology was generally ignored, and the fracture was usually equivalent to a smooth parallel plate model [15,16,17,18]. The equivalent parallel plate model was used for a long time since its’ proposal. However, with the deepening of understanding, it was found that the equivalent model had a great influence on the research results in specific applications. Therefore, in order to understand the characteristics of fluid migration more accurately, the characteristics of the fracture surface morphology cannot be ignored. As for the fluid migration in the fracture [19,20], proppant placement in the fracture [21,22], groundwater seepage [23], and other related situations [24], the actual results that considered the characteristics of the fracture surface morphology were indeed quite different from those of ignoring the morphology characteristics [19,20,21,22,23,24], which also verified that the characteristics of the fracture surface morphology had a significant influence on fluid migration in the fracture. Therefore, various parameters were proposed to describe the characteristics of fracture surface morphology. However, only one single parameter was usually selected for the characteristics description of fracture surface morphology in specific applications [19,20,21,22,23,24].
At present, few reports on the systematic characteristics description of fracture surface morphology are available in the literature, while no reports about the characteristics description of shale fracture surface morphology were found. In order to better serve the development of shale reservoirs, the shale samples from Barnett Shale was taken as the example. The characteristics of shale fracture surface morphology are described by parameters such as roughness, joint roughness coefficient, fractal dimension, tortuosity, and dip angle, so as to comprehensively describe the characteristics of shale fracture surface morphology and provide a solid foundation for accurately describing fluid migration in shale fractures. In addition, this research will also provide a reference for fracture description in other unconventional reservoirs.

2. Acquisition of Shale Fracture Surface Height Distribution

Three shale samples collected from Barnett Shale were fractured using the Brazilian test, and the fractures of the three shale samples are shown in Figure 1. The height distribution of the shale fracture surface was obtained by a three-dimensional optical profilometer, and both the lateral and longitudinal interval for data collection was 0.01 mm. The principle of data acquisition by a three-dimensional optical profilometer was listed as follows: First, a lateral coordinate was fixed, and the longitudinal height distribution under the lateral coordinate was collected successively with the data sampling accuracy taken as the interval until the data acquisition under this lateral coordinate was completed. Then, the lateral coordinate was moved by one longitudinal interval and fixed. The longitudinal height distribution under this lateral coordinate was also collected successively until all the information along this lateral coordinate was collected. The aforementioned steps were repeated until all the height distribution information of the shale fracture surface was collected.
According to its scanning principle, the three-dimensional optical profilometer transformed the morphology of the shale fracture surface into the point set of the profile. However, it could not obtain all the information of the shale fracture surface, as the scanned part of the adjacent profile was the blind area where the information could not be collected. Therefore, the blind area was excluded when calculating the relevant description parameters; in other words, the description parameters of shale fracture surface morphology were calculated according to the information from the profile obtained using a three-dimensional optical profilometer. The specific data acquisition process is shown in Figure 2.

3. Description Parameters of Shale Fracture Surface Morphology and Their Calculation

3.1. Height Distribution Pre-Processing of Shale Fracture Surface Morphology

As the fracture was split by the Brazilian test, in order to prevent the influence of end discontinuity and data acquisition, the central position of the shale fracture surface was firstly located. The rectangular region with 35 mm length and 21 mm width that lay in the central area of shale fracture surface was selected as the object for characteristics description, and the central position was used as the center of the rectangular region.
Prior to calculating the description parameters of shale fracture surface morphology, the height distribution of the shale fracture surface was processed as follows: for the data point set to be studied, the initial point of both lateral direction and longitudinal direction were reset to zero; then, the minimum height of shale fracture surface was regarded as the zero point of height distribution, and the minimum height was subtracted from each data point to obtain a new set of height distribution values for the shale fracture surface. The study focused on the characteristics of shale fracture surface morphology, which was only related to the fracture morphology and was independent of the location and the absolute height of the point. Therefore, the data processing had no effect on the results.
The fracture surface morphologies of three shale samples were reconstructed after the aforementioned pre-processing, and the reconstructed images are shown in Figure 3.

3.2. Characteristics Description of Shale Fracture Surface Morphology

So far, there are no systematic reports available regarding the characteristics of fracture surface morphology, and the characteristics description methods for fracture surface morphology are randomly selected. Generally speaking, roughness and fractal dimension are widely accepted to describe the characteristics of fracture surface morphology in specific studies that involve fluid migration, while joint roughness coefficient is a classical parameter recommended in the field of rock mechanics. However, tortuosity is often adopted when the study is related to both the characteristics description of fracture surface morphology and fluid migration. Besides, dip angle is a parameter used to describe the overall trend of fracture surface morphology. There are few reports that utilize dip angle alone to describe the characteristics of fracture surface morphology.
One single parameter was mostly used to describe the characteristics of fracture surface morphology, but the selection basis was not expounded [19,20,21,22,23,24]. Most of the selection was based on experience; however, it was doubtful whether one single parameter could comprehensively describe the characteristics of fracture surface morphology. Considering the disadvantage of the single description parameter used currently, the characteristics description parameters of fracture surface morphology were summarized, and their calculation methods were reasonably selected according to the scanning principle of the three-dimensional optical profilometer to systematically and comprehensively calculate the characteristics description parameters of shale fracture surface morphology.

3.2.1. Roughness

Roughness is not only the most basic parameter to describe the characteristics of fracture surface morphology but also the most widely verified parameter [25]. During the study of nonlinear flow through rock fracture networks, it was shown that the surface roughness that contributes to complex void geometries and streamlines structures could significantly reduce the critical Re for the flow regime transition [25]. Roughness belongs to microscopic geometrical error, and it refers to the asperity of the surface morphology with both small spacing and small peaks/valleys. The smaller the roughness, the smoother the surface morphology will be. The commonly recommended parameters for calculating roughness include standard deviation, skewness, kurtosis, arithmetic average height, the maximum height of the peak, etc. [26].
As for the description of fracture surface morphology, standard deviation is usually employed to calculate fracture roughness, which can be divided into overall standard deviation and sample standard deviation [27]. Taking the principle of data acquisition into account, the roughness of profiles along both lateral and longitudinal directions were calculated, and the arithmetic square root of the two was regarded as the roughness of shale fracture surface morphology. Sample standard deviation was introduced to characterize the roughness of shale fracture surface morphology, and the specific methods are shown in Equations (1)–(3):
ε x = i = 1 m Z i x Z x ¯ 2 m 1 ,
ε y = i = 1 n Z i y Z y ¯ 2 n 1 ,
ε = ε x ε y ,
The roughness of shale fracture surface morphology is shown in Table 1.
As shown in Table 1, the roughness of the shale fracture morphology was different, among 0.0834 mm–0.2427 mm. For the same shale fracture surface morphology, roughness varied greatly along different directions, with a maximum difference of seven times. Therefore, there was heterogeneity of roughness along different strikes.

3.2.2. Joint Roughness Coefficient

The joint roughness coefficient is a parameter recommended for characterizing the joint morphology in the field of rock mechanics [28]. Joint refers to those with fractures and no obvious relative displacement on either side of the fracture surface, e.g., with an essential class of fracture structure. However, joints and fractures are not easily distinguished in the specific application of joint roughness coefficient [29]. In the experimental study of fluid flow through fractured granite, a joint roughness coefficient was employed to qualitatively describe the joint, and specimens with fairly smooth joints subjected to low confining pressures showed a linear laminar flow region against fluid pressures [30]. Barton firstly proposed the concept of joint roughness coefficient (JRC) through the summarization of rough fracture profiles [31], and ten standard profiles (JRC values among 0–20) were proposed. The actual value of the joint roughness coefficient could be acquired by comparing the actual profile of fracture surface to that of the standard joint roughness coefficient. This method has been widely applied in related descriptions of fracture surface morphology since it was proposed and recommended, and many empirical equations used for the calculation of fracture surface parameters were developed according to their relationship to fit the joint roughness coefficient [32,33,34,35].
The most widely used statistical parameter Z2 was adopted to calculate the joint roughness coefficient of fracture profiles [33], while the Z2 along the lateral direction was listed as follows,
Z 2 x = 1 n i = 1 n ( z i + 1 z i ) 2 x i + 1 x i 1 2 ,
The equations of the relationship between Z2 and joint roughness coefficient developed by Tse and Cruden were chosen to calculate the joint roughness coefficient of the profile along the lateral direction [32], and it was formulated as,
J R C 3 D x = 32.2 + 32.47 1 m i = 1 m lg Z 2 x i ,
Both the calculations of the statistical parameter Z2y and the joint roughness coefficient JRC3Dy along the longitudinal direction were similar to that of the lateral direction. The joint roughness coefficient of shale fracture surface morphology is shown in Equation (6).
J R C 3 D = J R C 3 D x J R C 3 D y ,
The JRC3D of shale fracture surface morphology is shown in Table 2.
As shown in Table 2, the joint roughness coefficient of different samples varied greatly, among 2.5715–10.9368. The difference between the joint roughness coefficient along the lateral direction and that along the longitudinal direction was significant, up to two to three times.

3.2.3. Fractal Dimension

The concept of fractal dimension was proposed based on the self-similarity of geometric features, which is an important metric parameter to determine the disorder of geometrical morphology and the irregularity of complex geometries [36]. Compared to the integer dimension reflecting the static characteristics of the object, the fractal dimension represents the dynamic change process of the object. If it is extended to the dynamic behaviors and phenomena of nature, the fractal dimension is a correlated characterization of the whole system that is composed of small local characteristics in natural phenomena. That is, for an object, it can accurately reflect its irregularity and complexity using the dimension scale of the measured non-integer value; as a result, the dimension of the non-integer value is regarded as the fractal dimension [37]. During the experimental analysis of single-phase flow through rough fracture replicas, the compensation of turbulence effect developed earlier for rough fractures was compared to the smooth parallel plate model and was longer for the fractures with higher fractal dimensions. Transmissivity values showed a decrease with increasing fractal dimensions and normal loading and also exhibited directional-dependent behavior. The percentages of water-invaded wet planar areas showed a tendency to decrease with increasing fractal dimension [38].
For irregular profiles of complex fracture surfaces, the yard stick method and covering method are commonly used [39]. However, for fracture surfaces, it is impossible to directly cover the rough surface with two-dimensional Euclidean geometry with a certain scale.
Currently, the reported methods for calculating the fractal dimension of fracture surface include the triangular prism surface area method, projection covering method, and cube covering method [39]. The first two methods possess the defect that the calculation of rough surface area is approximate, which leads to the deviation of calculation results. In contrast, the cube covering method can effectively avoid the calculation deviation caused by approximation. Besides, the fractal dimension estimated by the cube covering method is the pure geometric fractal dimension. Each calculation step has an accurate method, and there is no approximate process during the calculation; therefore, the calculated fractal dimension is close to the true value.
The specific process of the cube covering method was provided as follows [40]: assuming there is a square grid on the plane XOY, and the length of its basic element unit is δ, the four corners of the square grid correspond to the four heights h(i, j), h(i, j + 1), h(i + 1, j), and h(i + 1, j + 1) (i, j is among (1, n − 1), and n is the number of measuring points each side), respectively. Cubes with side length δ were employed to cover the fracture surface, and the number of cubes in the covering area δ × δ was calculated, that is, the number of cubes Ni,j covering the fracture surface in the grid (i, j) was provided as follows,
N i , j = I N T δ 1 [ max ( h ( i , j ) , h ( i , j + 1 ) , h ( i + 1 , j ) , h ( i + 1 , j + 1 ) min ( h ( i , j ) , h ( i , j + 1 ) , h ( i + 1 , j ) , h ( i + 1 , j + 1 ) ) ] + 1
Therefore, the total number of cubes required to cover the entire fracture surface was
N δ = i , j = 1 n 1 N i , j ,
Then, the whole fracture surface was covered again by changing the measuring dimension, and the total number of cubes for covering the whole surface was calculated using Equation (8). According to the fractal theory, the following relationship could be obtained,
N δ ~ δ D ,
Moreover, the logarithms of N and δ at both ends of Equation (9) were, respectively, calculated. The relationship between them was fitted and plotted, and the slope was fractal dimension D.
The fractal dimension of shale fracture surface morphology is shown in Table 3.

3.2.4. Tortuosity

The concept of tortuosity was primarily proposed in the study of fluid seepage in porous media, and it is defined as the ratio of the fluid actual seepage path to its apparent length. Tortuosity represents the complexity of the fluid seepage path, so it essentially describes the complexity of the fluid seepage channel, which is also known as the channel tortuosity. Tortuosity has a negative effect on the permeability of the porous media, which means that the fluid seepage resistance increases with the increase of tortuosity [41]. Additionally, tortuosity is closely related to the physical properties of porous media, such as porosity and permeability [41].
The methods to study tortuosity mainly include casting thin sections, a capillary pressure curve, and empirical formula [42,43]. Currently, there are two main methods used to calculate the tortuosity of fracture surface morphology, which are, respectively, based on the profile and area of fracture surface morphology. Based on the acquisition principle of fracture surface height distribution, the calculation based on the profile of fracture surface morphology was selected to calculate the tortuosity of shale fracture surface morphology. For the characteristics of shale fracture surface morphology, three-dimensional distribution, both the lateral and longitudinal tortuosity along fluid seepage direction, were calculated, respectively, and the root mean square of the two tortuosities was taken as the tortuosity of shale fracture surface morphology. The specific equations are provided as follows:
τ x = 1 m i = 1 m j = 1 N x 1 Δ x 2 + z i + 1 , j z i , j 2 l x ,
τ y = 1 n i = 1 n j = 1 N y 1 Δ y 2 + z i , j + 1 z i , j 2 l y ,
τ = τ x τ y ,
The tortuosity of shale fracture surface morphology is shown in Table 4.

3.2.5. Dip Angle

The concept of dip angle was proposed during the description of fracture structure, and then dip angle was gradually adopted for the description of fracture structure. Dip angle is a parameter used to describe the trend of the fracture surface; however, the application of dip angle in fracture description is less than that of other parameters. There is almost no report of taking dip angle alone to describe the characteristics of fracture surface morphology, and it is generally used in combination with other parameters. During the study of gas flow in shale micro-fracture, dip angle combined with other parameters were employed to fit shale gas flow, and the results showed that dip angle had an inhibitory effect on gas flow in shale micro-fracture [44].
The angle of the profile can be divided into the angle with an upward trend and the angle with a downward trend along the profile; the influence of the two on fluid seepage is significantly different. The angle of the upward trend has an inhibitory effect on fluid seepage, while the angle of the downward trend has a positive effect. Therefore, the dip angle is divided into absolute dip angle and overall trend dip angle. The absolute dip angle refers to the absolute value accumulation of the angle along the profile, and upward and downward trends are not distinguished. The absolute angle represents the complexity of the profile. The overall trend dip angle refers to the superposition of the angles along the profile, distinguishing the angle of the upward trend from the angle of the downward trend. The overall trend dip angle represents the overall trend of the profile.
Based on the scanning principle of a three-dimensional optical profilometer, the dip angle of each profile was calculated successively, and the dip angle of the whole shale fracture surface morphology could be obtained by weighted averaging the dip angles of all profiles.
The absolute dip angle was calculated by Equation (13), and the arithmetic square root of the absolute dip angle product in both directions was taken as the absolute dip angle of shale fracture surface morphology. The details are shown in Table 5.
α = tan 1 ( 1 m 1 N x 1 i = 1 N x 1 z i + 1 z i Δ x ) ,
As can be seen in Table 5, the absolute dip angle of shale fracture surface morphology showed different characteristics along both the lateral and longitudinal directions. The absolute dip angle of shale sample 1 was basically the same along both directions, which were all around 24°. The absolute dip angle of shale sample 2 along the lateral direction was about 3° higher than that along the longitudinal direction, which was 19.406° and 16.235°, respectively. The absolute dip angle of shale sample 3 along the lateral direction was about 10° smaller than that along the longitudinal direction, which was 19.539° and 29.861°, respectively. In general, the absolute dip angle along the lateral direction was 19.406°–24.256°, and the absolute angle along the longitudinal direction was 16.235°–29.861°. The difference between shale samples along the lateral direction was about 5°, while it was about 13° along the longitudinal direction.
Meanwhile, the overall trend dip angle was calculated using Equation (14). The arithmetic square root of the overall trend dip angle product in both directions was taken as the overall trend dip angle of shale fracture surface morphology. The details are shown in Table 6.
β = tan 1 ( 1 m 1 N x 1 i = 1 N x 1 z i + 1 z i Δ x ) ,
The overall trend dip angle of shale fracture surface morphology is shown in Table 6. The overall trend dip angle of different shale samples varied greatly in both the lateral and longitudinal direction, and that of shale sample 3 was minimal, about 2.4°, while that of shale sample 1 was the maximum, about 12°. The difference between the overall trend dip angle of shale samples along the lateral direction was smaller, about 2.5°, while that along the longitudinal direction was larger, about 22.5°, roughly nine times larger than that along the lateral direction.

4. Relationship between Roughness and Other Description Parameters of Shale Fracture Surface Morphology

Summary of Shale Fracture Surface Morphology Characteristics

In the above study, the description parameters including roughness, joint roughness coefficient, fractal dimension, tortuosity, and dip angle of shale fracture surface morphology were systematically calculated. The characteristics of shale fracture surface morphology are summarized in Table 7.
Roughness is the most basic, simple, and widely used parameter to describe the characteristics of fracture surface morphology. The joint roughness coefficient, fractal dimension, tortuosity, absolute dip angle, and overall trend dip angle were positively correlated with roughness, but the physical meanings described by each parameter were different, and the sensitivity of each parameter to roughness varied greatly.

5. Conclusions

(1)
The shale samples collected from Barnett Shale were fractured using the Brazilian test, and the information regarding shale fracture surface morphology was obtained. According to the collecting principle of shale fracture surface morphology, the calculation methods of roughness, joint roughness coefficient, fractal dimension, tortuosity, and dip angle were proposed.
(2)
Roughness, joint roughness coefficient, fractal dimension, tortuosity, and dip angle were all directional, and they varied greatly along different directions.
(3)
The parameters of roughness, joint roughness coefficient, fractal dimension, tortuosity, absolute dip angle, and overall trend dip angle were among 0.0834 mm–0.2427 mm, 2.5715–10.9368, 2.1000–2.1364, 1.0732–1.1879, 17.7498°–24.5941°, and −3.7223°–13.3042°, respectively.
(4)
Joint roughness coefficient, fractal dimension, tortuosity, and dip angle were all positively correlated with roughness. However, the sensitivity of each parameter to roughness varied greatly.

Author Contributions

T.S. and Z.Z. conceived and designed the experiments; J.P. performed the experiments; G.Q., H.W., and R.G. calculated the parameters; G.Q. and J.S. wrote the paper; K.Z. drew the figures. All authors have read and agreed to the published version of the manuscript.

Funding

The study was funded by Youth Program of National Natural Science Foundation of China (51804255), the Natural Science Basic Research Program of Shaanxi (2019JQ-488), the Opening Fund of Key Laboratory of Unconventional Oil & Gas Development (China University of Petroleum (East China) (19CX05005A-6), Ministry of Education and the Fundamental Research Funds for the Central Universities (19CX05005A-6).

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Acknowledgments

The authors would like to thank Randy Doly Hazllet, for his help in obtaining the shale samples. We appreciate the valuable comments from the editors and anonymous reviewers.

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclature

ε, εx, εyroughness of the whole shale fracture surface morphology, along the lateral and longitudinal direction, respectively, mm
Zix, Ziymean height of ith profile along the lateral and longitudinal direction, respectively, mm
Z x ¯ , Z y ¯ mean height of the profile along the lateral and longitudinal direction, respectively, mm
ziheight of an ith point on the shale fracture surface, mm
nnumber of points on shale fracture surface
Z2root mean square first derivative values
JRC2Dtwo-dimensional joint roughness coefficient
JRC3Dthree-dimensional joint roughness coefficient
Lthe length of shale fracture
xithe position of an ith point on the profile
INTingetral function
lx,lythe apparent length of a profile in the lateral and longitudinal direction, respectively, mm
τ, τx, τytortuosity of the whole shale fracture surface morphology, along the lateral and longitudinal direction, respectively
Δx, Δyscanning interval along the lateral and longitudinal direction, respectively
m, nnumber of profiles along the lateral and longitudinal direction, respectively
θdip angle of shale fracture surface morphology
mnumber of profiles along lateral direction
Nxnumber of points on a profile along the longitudinal direction

References

  1. Mendhe, V.A.; Kamble, A.D.; Bannerjee, M.; Mishra, S.; Mukherjee, S.; Mishra, P. Evaluation of shale gas reservoir in Barakar and barren measures formations of north and south Karanpura Coalfields, Jharkhand. J. Geol. Soc. India 2016, 88, 305–316. [Google Scholar] [CrossRef]
  2. Wang, Y.; Liang, C.; Sun, X. Shale oil reservoir characteristics and enrichment in the Jiyang depression, Bohai Bay Basin. East China J. Earth Sci. 2017, 28, 977–986. [Google Scholar] [CrossRef]
  3. Jiang, Z.-X.; Li, Z.; Li, F.; Pang, X.-Q.; Yang, W.; Liu, L.-F.; Jiang, F.-J. Tight sandstone gas accumulation mechanism and development models. Pet. Sci. 2015, 12, 587–605. [Google Scholar] [CrossRef] [Green Version]
  4. Chen, X.; Chen, L.; Jiang, S.; Liu, A.; Luo, S.; Li, H.; Li, P.; Chen, P. Evaluation of Shale Reservoir Quality by Geophysical Logging for Shuijingtuo Formation of Lower Cambrian in Yichang Area, Central Yangtze. J. Earth Sci. 2021, 32, 766–777. [Google Scholar] [CrossRef]
  5. Ji, W.; Song, Y.; Rui, Z.; Meng, M.; Huang, H. Pore characterization of isolated organic matter from high matured gas shale reservoir. Int. J. Coal Geol. 2017, 174, 31–40. [Google Scholar] [CrossRef]
  6. Shan, C.; Zhang, T.; Wei, Y.; Zhang, Z. Shale gas reservoir characteristics of Ordovician-Silurian formations in the central Yangtze area, China. Front. Earth Sci. 2017, 11, 184–201. [Google Scholar] [CrossRef]
  7. Jia, B.; Tsau, J.S.; Barati, R. A review of the current progress of CO2 injection EOR and carbon storage in shale oil reservoirs. Fuel 2019, 236, 404–427. [Google Scholar] [CrossRef]
  8. Junira, A.; Biancardi, S.; Ganjdanesh, R.; Sepehrnoori, K.; Yu, W.; Ambrose, R. In situ fluid and injection gas compatibility of huff-n-puff gas EOR in a low-permeability shale-sand-carbonate reservoir. J. Pet. Sci. Eng. 2022, 208, 109587. [Google Scholar] [CrossRef]
  9. Tu, J. Study of Surfactant-Based Shale Oil EOR Under High Confining Pressure Conditions. In SPE Annual Technical Conference and Exhibition; Society of Petroleum Engineers: Calgary, AB, Canada, 2019. [Google Scholar] [CrossRef]
  10. Li, L.; Sheng, J.J. Experimental study of core size effect on CH4 huff-n-puff enhanced oil recovery in liquid-rich shale reservoirs. J. Nat. Gas Sci. Eng. 2016, 34, 1392–1402. [Google Scholar] [CrossRef] [Green Version]
  11. Zhong, J.; Liu, S.; Ma, Y.; Yin, C.; Liu, C.; Li, Z.; Liu, X.; Li, Y. Macro-fracture mode and micro-fracture mechanism of shale. Pet. Explor. Dev. 2015, 42, 269–276. [Google Scholar] [CrossRef]
  12. Odling, N.E.; Gillespie, P.; Bourgine, B.; Castaing, C.; Chiles, J.P.; Christensen, N.P.; Fillion, E.; Genter, A.; Olsen, C.; Thrane, L.; et al. Variations in fracture system geometry and their implications for fluid flow in fractures hydrocarbon reservoirs. Pet. Geosci. 1999, 5, 373–384. [Google Scholar] [CrossRef]
  13. Yin, P.; Zhao, G.-F. Numerical simulation of fluid flow through deformable natural fracture network. Geomech. Geophys. Geo-Energy Geo-Resour. 2016, 2, 343–363. [Google Scholar] [CrossRef] [Green Version]
  14. Alvarez, L.L.; Guimarães, L.J.D.N.; Gomes, I.F.; Beserra, L.; Pereira, L.C.; de Miranda, T.S.; Maciel, B.; Barbosa, J.A. Impact of Fracture Topology on the Fluid Flow Behavior of Naturally Fractured Reservoirs. Energies 2021, 14, 5488. [Google Scholar] [CrossRef]
  15. Kishida, K.; Kobayashi, K.; Hosoda, T.; Fujita, T.; Kishi, H.; Kuzuha, Y.; Nobuto, J. Development of grout injection model to single fracture in considering inertia term and its application on parallel plate experiments. J. Soc. Mater. Sci. Jpn. 2012, 61, 245–252. [Google Scholar] [CrossRef] [Green Version]
  16. Klimczak, C.; Schultz, R.A.; Parashar, R.; Reeves, D.M. Cubic law with aperture-length correlation: Implications for network scale fluid flow. Hydrogeol. J. 2010, 18, 851–862. [Google Scholar] [CrossRef]
  17. Brown Stephen, R. Transport of fluid and electric current through a single fracture. J. Geophys. Res. Solid Earth 1989, 94, 9429. [Google Scholar] [CrossRef]
  18. Zhou, B.; Luo, X.; Loggia, D.; Hou, P.; Wang, Z. Experiments and Characterization of Oil Migration in a Single Fracture. Acta Geol. Sin. 2006, 80, 454. [Google Scholar] [CrossRef]
  19. Ouyang, Z.; Liu, D.; Cai, Y.; Yao, Y. Investigating the Fractal Characteristics of Pore-Fractures in Bituminous Coals and Anthracites through Fluid Flow Behavior. Energy Fuels 2016, 30, 10348–10357. [Google Scholar] [CrossRef]
  20. Stoll, M.; Huber, F.; Trumm, M.; Enzmann, F.; Meinel, D.; Wenka, A.; Schill, E.; Schäfer, T. Experimental and numerical investigations on the effect of fracture geometry and fracture aperture distribution on flow and solute transport in natural fractures. J. Contam. Hydrol. 2019, 221, 82–97. [Google Scholar] [CrossRef]
  21. Raimbay, A.; Babadagli, T.; Kuru, E.; Develi, K. Quantitative and visual analysis of proppant transport in rough fractures. J. Nat. Gas Sci. Eng. 2016, 33, 1291–1307. [Google Scholar] [CrossRef]
  22. Li, J.; Liu, P.; Kuang, S.; Yu, A. Visual lab tests: Proppant transportation in a 3D printed vertical hydraulic fracture with two-sided rough surfaces. J. Pet. Sci. Eng. 2021, 196, 107738. [Google Scholar] [CrossRef]
  23. Ni, X.-D.; Niu, Y.-L.; Wang, Y.; Yu, K. Non-Darcy Flow Experiments of Water Seepage through Rough-Walled Rock Fractures. Geofluids 2018, 2018, 8541421. [Google Scholar] [CrossRef] [Green Version]
  24. Verberg, R.; Ladd, A. Simulation of chemical erosion in rough fractures. Phys. Rev. E 2002, 65, 056311. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Liu, R.; Li, B.; Jiang, Y. Critical hydraulic gradient for nonlinear flow through rock fracture networks: The roles of aperture, surface roughness, and number of intersections. Adv. Water Resour. 2016, 88, 53–65. [Google Scholar] [CrossRef]
  26. Horváth, R.; Czifra, Á.; Drégelyi-Kiss, Á. Effect of conventional and non-conventional tool geometries to skewness and kurtosis of surface roughness in case of fine turning of aluminium alloys with diamond tools. Int. J. Adv. Manuf. Technol. 2015, 78, 297–304. [Google Scholar] [CrossRef]
  27. Bland, M. Estimating Mean and Standard Deviation from the Sample Size, Three Quartiles, Minimum, and Maximum. Int. J. Stats Med. Res. 2015, 4, 57–64. [Google Scholar] [CrossRef]
  28. Li, Y.; Zhang, Y. Quantitative estimation of joint roughness coefficient using statistical parameters. Int. J. Rock Mech. Min. Sci. 2015, 77, 27–35. [Google Scholar] [CrossRef] [Green Version]
  29. Mo, P.; Li, Y. Estimating the three-dimensional joint roughness coefficient value of rock fractures. Bull. Eng. Geol. Environ. 2019, 78, 857–866. [Google Scholar] [CrossRef] [Green Version]
  30. Ranjit, P.G. An experimental study of single and two-phase fluid flow through fractured granite specimens. Environ. Earth Sci. 2010, 59, 1389–1395. [Google Scholar] [CrossRef]
  31. Barton, N.; Choubey, V. The shear strength of rock joints in theory and practice. Rock Mech. Rock Eng. 1977, 10, 1–54. [Google Scholar] [CrossRef]
  32. Tse, R.; Cruden, D.M. Estimating joint roughness coefficients. Int. J. Rock Mech. Min. Sci. Geomech. Abstr. 1979, 16, 303–307. [Google Scholar] [CrossRef]
  33. Tatone, B.S.A.; Grasselli, G. A new 2D discontinuity roughness parameter and its correlation with JRC. Int. J. Rock Mech. Min. Sci. 2010, 47, 1391–1400. [Google Scholar] [CrossRef]
  34. Mirzaghorbanali, A.; Nemcik, J.; Aziz, N. Effects of shear rate on cyclic loading shear behaviour of rock joints under constant normal stiffness conditions. Rock Mech. Rock Eng. 2014, 47, 1931–1938. [Google Scholar] [CrossRef]
  35. Sanei, M.; Faramarzi, L.; Goli, S.; Fahimifar, A.; Rahmati, A.; Mehinrad, A. Development of a new equation for joint roughness coefficient (JRC) with fractal dimension: A case study of Bakhtiary Dam site in Iran. Arab. J. Geosci. 2015, 8, 465–475. [Google Scholar] [CrossRef]
  36. Schweinhart, B. Fractal dimension and the persistent homology of random geometric complexes. Adv. Math. 2020, 372, 107291. [Google Scholar] [CrossRef]
  37. Liu, Y.; Wang, Y.; Chen, X.; Zhang, C.; Tan, Y. Two-stage method for fractal dimension calculation of the mechanical equipment rough surface profile based on fractal theory. Chaos Solitons Fractals 2017, 104, 495–502. [Google Scholar] [CrossRef]
  38. Develi, K.; Babadagli, T. Experimental and visual analysis of single-phase flow through rough fracture replicas. Int. J. Rock Mech. Min. Sci. 2015, 73, 139–155. [Google Scholar] [CrossRef]
  39. Zhou, H.; He, X.; Kwasniewski, M. Fractal Dimension of Rough Surface Estimated by the Cubic Covering Method. Tribology 2000, 20, 455–459. [Google Scholar] [CrossRef]
  40. Zuo, J.; Wei, X.; Pei, J.; Zhao, X. Investigation of meso-failure behaviors of Jinping marble using SEM with bending loading system. J. Rock Mech. Geotech. Eng. 2015, 7, 593–599. [Google Scholar] [CrossRef] [Green Version]
  41. Matyka, M.; Khalili, A.; Koza, Z. Tortuosity-porosity relation in porous media flow. Phys. Rev. E 2008, 78, 026306. [Google Scholar] [CrossRef] [Green Version]
  42. Ezzatabadipour, M.; Zahedi, H. A Novel Method for Streamline-Based Tortuosity Calculation and Investigation of Obstacles Shape Effect on Tortuosity in Porous Media with Random Elliptical Obstacles Using Lattice Boltzmann Method. Transp. Porous Media 2021, 136, 103–124. [Google Scholar] [CrossRef]
  43. Ziarani, A.S.; Aguilera, R. Pore-throat radius and tortuosity estimation from formation resistivity data for tight-gas sandstone reservoirs. J. Appl. Geophys. 2012, 83, 65–73. [Google Scholar] [CrossRef]
  44. Qu, G.; Qu, Z.; Hazlett, R.; Freed, D.; Mustafayev, R. Geometrical description and permeability calculation about shale tensile micro-fractures. Pet. Explor. Dev. 2016, 43, 124–130. [Google Scholar] [CrossRef]
Figure 1. Fractures of shale samples. (a) Shale fracture 1; (b) Shale fracture 2; (c) Shale fracture 3.
Figure 1. Fractures of shale samples. (a) Shale fracture 1; (b) Shale fracture 2; (c) Shale fracture 3.
Processes 10 00401 g001
Figure 2. Data acquisition principle of the three-dimensional optical profilometer.
Figure 2. Data acquisition principle of the three-dimensional optical profilometer.
Processes 10 00401 g002
Figure 3. Shale fracture surface morphologies (the subfigures (ac) are the reconstructed images of the rectangular region with 35 mm length and 21 mm width on the three shale surface morphology, and they are the zone for the calculation of characteristics description).
Figure 3. Shale fracture surface morphologies (the subfigures (ac) are the reconstructed images of the rectangular region with 35 mm length and 21 mm width on the three shale surface morphology, and they are the zone for the calculation of characteristics description).
Processes 10 00401 g003aProcesses 10 00401 g003b
Table 1. Roughness of shale fracture surface morphology.
Table 1. Roughness of shale fracture surface morphology.
Shale Sample No.εx, mmεy, mmε, mm
10.30890.19070.2427
20.03060.22710.0834
30.09330.10800.1004
Table 2. JRC3D of shale fracture surface morphology.
Table 2. JRC3D of shale fracture surface morphology.
Shale Sample No.JRC3DxJRC3DyJRC3D
19.216312.977910.9368
23.97981.66162.5715
39.57344.27716.3989
Table 3. Fractal dimension of shale fracture surface morphology.
Table 3. Fractal dimension of shale fracture surface morphology.
Shale Sample No.Fractal Dimension D
12.1364
22.1000
32.1187
Table 4. Tortuosity of shale fracture surface morphology.
Table 4. Tortuosity of shale fracture surface morphology.
Shale Sample No.τxτyτ
11.17431.20151.1879
21.08431.06211.0732
31.09161.21491.1516
Table 5. Absolute dip angle of shale fracture surface morphology.
Table 5. Absolute dip angle of shale fracture surface morphology.
Shale Sample No.αx, °αy, °α, °
124.25624.93724.5941
219.40616.23517.7498
319.53929.86124.1548
Table 6. Overall trend dip angle about shale fracture surface morphology.
Table 6. Overall trend dip angle about shale fracture surface morphology.
Shale Sample No.βx, °βy, °β, °
18.53120.74813.3042
27.037−1.969−3.7223
35.0197.4326.1075
Table 7. Characteristics of shale fracture surface morphology.
Table 7. Characteristics of shale fracture surface morphology.
Shale Sample No.123
Roughness, mm0.24270.08340.1004
Joint roughness coefficient10.93682.57156.3989
Fractal dimension2.13642.10002.1187
Tortuosity1.18791.07321.1516
Absolute dip angle, °24.594117.749824.1548
Overall trend dip angle, °13.3042−3.72236.1075
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Qu, G.; Shi, T.; Zhang, Z.; Su, J.; Wei, H.; Guo, R.; Peng, J.; Zhao, K. Characteristics Description of Shale Fracture Surface Morphology: A Case Study of Shale Samples from Barnett Shale. Processes 2022, 10, 401. https://doi.org/10.3390/pr10020401

AMA Style

Qu G, Shi T, Zhang Z, Su J, Wei H, Guo R, Peng J, Zhao K. Characteristics Description of Shale Fracture Surface Morphology: A Case Study of Shale Samples from Barnett Shale. Processes. 2022; 10(2):401. https://doi.org/10.3390/pr10020401

Chicago/Turabian Style

Qu, Guanzheng, Tao Shi, Zheng Zhang, Jian Su, Haitao Wei, Rui Guo, Jiao Peng, and Kairui Zhao. 2022. "Characteristics Description of Shale Fracture Surface Morphology: A Case Study of Shale Samples from Barnett Shale" Processes 10, no. 2: 401. https://doi.org/10.3390/pr10020401

APA Style

Qu, G., Shi, T., Zhang, Z., Su, J., Wei, H., Guo, R., Peng, J., & Zhao, K. (2022). Characteristics Description of Shale Fracture Surface Morphology: A Case Study of Shale Samples from Barnett Shale. Processes, 10(2), 401. https://doi.org/10.3390/pr10020401

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

Article Metrics

Back to TopTop