Next Article in Journal
A Novel Grain-Based DEM Model for Evaluating Surface Integrity in Scratching of RB-SiC Ceramics
Next Article in Special Issue
Using Iron Tailings for Phosphate Removal in Cemented Phosphogypsum (PG) Backfill
Previous Article in Journal
Effect of Reinforcement Size on Mechanical Behavior of SiC-Nanowires-Reinforced 6061Al Composites
Previous Article in Special Issue
Research of Dynamic Tensile Properties of Five Rocks under Three Loading Modes Based on SHPB Device
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

An Evaluation Index of Fracability for Reservoir Rocks Based on Fracture Process Zone

1
Key Laboratory of Shale Gas and Geoengineering, Institute of Geology and Geophysics, Chinese Academy of Sciences, Beijing 100029, China
2
Innovation Academy for Earth Science, CAS, Beijing 100029, China
3
Chinese Society for Rock Mechanics and Engineering, Beijing 100029, China
*
Author to whom correspondence should be addressed.
Materials 2022, 15(23), 8485; https://doi.org/10.3390/ma15238485
Submission received: 25 October 2022 / Revised: 21 November 2022 / Accepted: 24 November 2022 / Published: 28 November 2022

Abstract

:
A reliable evaluation method for the fracability (i.e., ability to generate abundant cracks) of reservoir rocks is a critical issue for maximum hydraulic fracturing efficiency. Most previous fracability indices lacked enough rationality and practicability, and thus could not consistently provide a reliable evaluation. We suggest a new fracability index called crack tolerance, which is represented by the maximum radius of the fracture process zone at the crack tip of a cracked chevron notched Brazilian disk specimen, which corresponds to the critical state for unstable propagation of the notched crack. In experiments and simulations based on the discrete element method, we showed quantitative methods to conveniently determine the value of the crack tolerance and showed that specimens with a greater crack tolerance generated more cracks before rupture and had complex morphologies, which would indicate stronger fracability. The crack tolerance can well characterize the effects of structural and loading conditions, including the grain size heterogeneity, bedding orientation, and environmental temperature, on fracability, and the inherent heterogeneity of rock is the physical basis for it as a fracability evaluation index. Our studies showed the rationality and practicability of this index and provide hints for how to produce abundant complex cracks in reservoirs.

1. Introduction

Human beings’ mining engineering and energy resource exploitation extensively involve the generation and propagation of cracks within rock materials. Hydraulic fracturing is widely used to enhance the fluid conductivity of reservoirs of oil, gas, and geothermal resources. A reliable evaluation of the rock fracability (i.e., ability to generate abundant cracks) is important for hydraulic fracturing [1].
Brittleness, which is generally viewed as a property (or ability) of solid material that ruptures with little appreciable permanent deformation, has long been considered approximately equivalent to fracability, because it shows empirical relevance to the possibility of crack propagation: reservoir comprising brittle rocks usually responds well to stimulation, whereas preexisting and hydraulic fractures tend to heal rather than to propagate in a less brittle reservoir. This is probably attributed to less energy consumed by the ductile deformation of brittle rock materials [2].
In the past decades, a variety of brittleness indices have been developed to evaluate its effect [3,4], which can be classified into several broad categories:
(1)
Based on mineral composition (e.g., [5]), especially the weight or volume proportion of hard minerals such as quartz: a positive correlation seems to exist between the brittleness and mineral contents of rocks. However, such indices do not consider many other factors that also contribute to brittleness, such as grain size and loading conditions.
(2)
Based on elastic parameters (e.g., [6]): for example, rocks with a large Young’s modulus and small Poisson’s ratio are assumed more brittle. However, such indices can be controversial because many laboratories and in situ observations [7,8] contradict this assumption.
(3)
Based on strength: for example, one such index is the ratio of tensile and compressive strengths [9]. Such indices are easily measured, but they lack a physical correlation to brittleness and cracking propagation mechanisms. Thus, these indices may return similar values for various types of rocks with different levels of brittleness.
(4)
Based on characteristics of the stress–strain curve such as the relative stress drop, post-peak modulus, and various combinations [10,11]: these indices characterize rock brittleness well and are widely used in predicting the rockburst proneness. However, high brittleness does not consistently represent strong fracability because brittle rock can also act as a barrier to hydraulic fracturing [12].
In summary, many brittleness indices currently popular in fracability evaluation for reservoir lacks mechanical relevance to the rock cracking process. On the other hand, the evaluation indices used in other areas (e.g., those used to estimate rock cuttability [13]) are usually inapplicable for reservoir fracability evaluation owing to the essential differences of physical meaning between brittleness and fracability. Thus far, few evaluation indices of rock fracability meet the following requirements [3]:
(1)
Has a firm physical basis;
(2)
Consider the heterogeneity of rock material;
(3)
Be convenient to measure;
(4)
Characterize the effects of loading conditions.
To address this issue, we propose a new evaluation index for rock fracability that we call the crack tolerance. See Section 2 for its definition. Section 3 and Section 4 show the experimental measurement of this new index and the corresponding numerical simulation results, respectively, to demonstrate the rationality of the index. Based on these analyses, the effects of several characteristics of the rock materials on the crack tolerance are discussed in Section 5.
This study demonstrated the physical rationality of the crack tolerance as an evaluation index and analyzed the effects of the rock structure and loading conditions on the crack tolerance in an effort to extend our understanding of rock fracability and provide hints for how to produce more cracks in the reservoir.

2. Fracture Process Zone and Crack Tolerance

Numerous researchers have revealed that the propagation of macroscopic cracks within rock under tension is attributed to progressive generation, interaction, and nucleation of micro-cracks from the macroscopic crack tips as follows [14]. When the imposed tensile load is small, only a few independent micro-cracks can arise around each crack tip (Figure 1a). As the tensile load increases, the distribution range of the micro-cracks expands, and their density increases. They interact with each other and coalesce (Figure 1b) to cause a gradual macroscopic propagation of the preexisting crack (Figure 1c,d). These micro-cracks indicate nonlinear deformation in the region around a crack tip preceding crack unstable propagation, which is referred to as a fracture process zone (FPZ) [14].
Crack propagation in tensile mode is most common in hydraulic fracturing because the effect of hydraulic pressure imposed on the crack surface approximates remote tensile stress in nature; additionally, rocks have a much lower tensile strength than compressive and shear strengths. Thus, cracks easily propagate driven by an injected fluid. The principal stresses at a tensile crack tip can be described as [15]
σ 1 σ 2 = K I 2 π r cos θ 2 1 + sin θ 2 cos θ 2 1 sin θ 2
where σ1 and σ2 are maximum and intermediate principal stresses, KI is tensile stress intensity factor, r and θ are polar radius and polar angle for polar coordinate system from the tip. Note that the minimum principal stress not listed here equals to zero. The range of FPZ (i.e., its size) is calculated based on the hypothesis that nonlinear deformation occurs within a region around crack tip when the local stress state satisfies a certain criterion (e.g., tensile strength criterion for rock materials, von Mises criterion for metal materials). The tensile cracks are assumed to propagate parallel to their own plane (i.e., θ = 0) when the σ1 reaches the tensile strength of the rock (σt), because the critical state of crack propagation is attained, which corresponds to the maximum size of the FPZ:
σ t = K IC 2 π r c cos 0 1 + sin 0
which leads to
rc = (KIC/σt)2/(2π),
where KIC is the tensile fracture toughness, and rc is the maximum FPZ size. In this context, the FPZ is represented by a circle centered on a fracture tip [14] (Figure 1e), and rc is the radius of a circular FPZ.
A large rc would indicate that micro-cracks are distributed within a large FPZ in front of a preexisting crack tip. It would also suggest a considerable number of micro-cracks within the FPZ because a preexisting crack will not propagate until the micro-crack density is high enough to reach a critical level [16]. Therefore, rc may characterize the maximum number of micro-cracks generated in the preparation stage for macroscopic crack propagation. In other words, rc can be used to indicate the ability of a rock to tolerate micro-cracks before crack unstable propagation. For this reason, we refer to rc as the crack tolerance. The crack morphology may also depend on the crack tolerance because a large rc would indicate an extensive distribution of micro-cracks, which would likely result in irregular and branch cracks.
The concepts of the FPZ and rc derive from the propagation process of a single crack with specific boundary conditions. Nevertheless, this process represents the inherent mechanical rule of crack generation within rocks because each crack started as an FPZ. Based on this understanding, the crack tolerance may reflect the potential of a given rock stratum to generate abundant cracks. Recent studies [17,18] have shown that rock specimens with a larger FPZ produce more fragments, which suggests greater fracability and provides evidence supporting our hypothesis. The maximum FPZ radius has exhibited dependence on the structure [19,20] and loading conditions [21] of rock. Thus, we conducted experiments and numerical simulations to analyze their effects on crack tolerance and demonstrate its rationality as an evaluation index.

3. Experiments

3.1. Specimens

We used marble, shale, and sandstone collected from Xishan, Beijing for experiments because marble was observed in some geothermal reservoirs, and shale and sandstone are representative lithologies comprising oil and gas reservoirs. The marble was divided into types A and J (Table 1): marble A totally constituted by calcite had a greater average grain size and was more heterogeneous as defined by Han et al. [22], and marble J mainly consisting of dolomite had an equigranular texture. The microscopy observation and X-ray diffraction (XRD) analysis showed that the shale with fine grains consisted of quartz (55.4%), plagioclase (6.2%), and clay minerals (38.4%, brown grains in Table 1). The quartz and clay minerals were alternatively layered. The sandstone consisted of quartz (69.5%), plagioclase (22.1%), and potassium feldspar (8.4%), and these xenomorphic grains have similar sizes (~2 mm). Most plagioclase grains experienced sericitization.

3.2. Experimental Methodology

The cracked chevron notched Brazilian disk (CCNBD) test involves the formation of FPZs at the two tips of a prefabricated notched crack, which is analogous to a natural crack. Therefore, the CCNBD test is applicable to evaluating crack tolerance. According to Equation (3), quantifying the crack tolerance requires determining the tensile fracture toughness KIC and tensile strength σt, which are measured by the CCNBD and Brazilian disk (BD) tests, respectively, as recommended by the International Society for Rock Mechanics (ISRM) and American Society of Testing Materials (ASTM).
The notched crack of each CCNBD specimen was created by a 1 mm thick circular diamond saw. To ensure cutting accuracy, the expected locations of the circular center and the initial and final chevron notched cracks were marked on each disk. We measured the actual values of the parameters shown in Figure 2a,b and confirmed that the dimensionless parameters α1 and αB of all CCNBD specimens were within the valid range (Figure 2c). The method reported by Fowell et al. [23] was used to calculate the KIC:
KIC = PmaxY*min/BR1/2,
where Pmax is the peak applied axial load in the CCNBD test and Y*min is the critical dimensionless stress intensity value. This is determined by
Y* = u · exp (v · α1),
where u and v are geometric constants that are determined by α0 and αB as reported by Fowell et al. [23].
The thickness (B’) and diameter (D’) of the BD specimens were set identical to those of the CCNBD specimens to eliminate the size effect on the calculated crack tolerance. The B’-to-D’ ratio was within the range recommended by the ASTM of 0.2–0.75 [24]. The σt was calculated as follows:
σt = 2PmaxBD’,
where Pmax is the peak applied axial load in the BD test.
Each CCNBD or BD test (Figure 3a,b) was performed at a constant displacement rate of 0.06 mm/min by an MTS servo-control testing machine (series CMT) with a maximum loading force of 100 kN. This machine is equipped with an SNAS GDS-300 environmental chamber controlled by a WK650 controller (Figure 3c,d). These apparatuses permit environmental temperatures within the chamber up to 200 °C by electrical heaters (Figure 3b). To investigate the effect of temperature, several sandstone specimens were placed in the chamber at 75 or 125 °C for 1 h before the tests began, so that the notched crack propagated within rocks under higher temperatures. Other tests were performed at room temperature (~25 °C). The bedding planes of the shale specimens were set perpendicular (horizontal) or parallel (vertical) to the notched cracks to analyze the effect of the bedding orientation.
BD tests were conducted on at least three specimens in parallel with the same lithology, bedding orientation, and temperature, and the average strength was taken as the tensile strength for the corresponding set of conditions. The KIC of each CCNBD specimen and the above average σt were used in Equation (3) to calculate the crack tolerance.

3.3. Experimental Results

In the CCNBD test, the marble A specimens with stronger heterogeneity had greater tensile strength, fracture toughness, and crack tolerance than the marble J specimens, with relatively homogeneous small grains (Figure 4a). White patches indicating FPZs [25] appeared in front of the notched crack tips (Figure 5a) as the peak loads of the marble A specimens were approached. The patches corresponded to the sparkling areas on the rupture surface (Figure 5b), which may imply breaking cleavages of grains. However, such patches were not observed for the marble J specimens (Figure 5c,d), and neither were the discernible sparkling areas. Furthermore, the main cracks in the marble A specimens had branches causing more fragments (Figure 5b) while the crack in the marble J specimens propagated along a straight path (Figure 5c). These phenomena suggest that a more heterogeneous grain size corresponds to a larger crack tolerance and thus a stronger ability for crack generation.
The mean crack tolerance of the shale specimens was less with a vertical bedding orientation than with a horizontal orientation (Figure 4b). The tensile strength and fracture toughness displayed similar variation trends with bedding orientation. Similar results can be acquired based on the data from Wang [26]. With a vertical orientation, the main crack of the specimen propagated along the bedding planes (Figure 6a), which generated a smooth rupture surface (Figure 6b). In contrast, with a horizontal orientation, the main crack spanned across bedding planes, and the path with steps was more irregular (Figure 6c,d). This is because the main crack was offset or even bifurcated when it encountered a bedding plane. The branch cracks were captured by bedding planes and then propagated along them, thereby their morphologies were smooth.
The crack tolerance value of the sandstone specimens consistently declined as the environmental temperature rises from 25 °C to 125 °C, while the tensile strength and fracture toughness exhibited V-shaped trends within this temperature range (Figure 4c). It is difficult to identify changes in crack morphology with the rising temperature with the naked eye (Figure 7a,c,e). According to the edge of their rupture surface, we speculated that the main cracks in the specimens at 125 °C may propagate along less curved paths than the specimens at lower temperatures did (Figure 7b,d,f). The variations of the crack tolerance value and crack morphologies imply that high temperatures possibly reduce rock fracability.
The crack tolerance preliminarily showed an ability to address the aforementioned problems of previous fracability indices. Firstly, this index has a firm physical basis derived from the FPZ size, representing the nonlinearity of deformation due to micro-crack generation before the macroscopic propagation of the crack. Additionally, the formation of the FPZ is the inherent mechanical behavior of heterogeneous rock materials, and the FPZ size highly depends on the degree of heterogeneity as the previous [27] and present experiment results revealed. From the aspect of practicability, the crack tolerance value can be determined conveniently in the laboratory because BD and CCNBD tests are very common rock mechanical tests, and the small-size specimens they use can be easily obtained from cores. Finally, this index may characterize the effects of structure and loading conditions on fracability to an extent, as the tests on the shale and sandstone showed.

4. Numerical Analysis with the Discrete Element Method

4.1. Particle Flow Code

To test the rationality of the crack tolerance as an evaluation index for rock fracability, we adopt particle flow code in two dimensions (PFC2D), which is widely used for discrete element method (DEM). Rock was modeled as a dense packing of non-uniform-sized and inter-bonded circular particles using this method, and its mechanical behavior relied on the microscale properties and constitutive relations of the bonded contacts between the particles. Following Newton’s laws of motion, the force acting at each contact were updated with the particle movements during the simulation process, and the breakage of bonds representing crack generation [28] occurred when a component of the contact force satisfied a certain criterion.
We used the experimental results for marble A and J as examples for the DEM simulation because the marble contained polygonal minerals >1 mm in size, which allowed us to implement a grain-based model (GBM) with a polygon-tessellation grain boundaries [29]. Such a model takes the mineral grain texture into account (Figure 8a), making a simulation more vivid. The modeling method for GBM of marble refers to [27]. Soft-bonded [30] and smooth-joint models [31] were employed to express bonded and unbonded behaviors characteristic of intra-grain and inter-grain contacts (Figure 8b–d), respectively.

4.2. Model Setup and Parameter Calibration

Based on the grain size distributions in Table 1, four circular 75 mm GBMs were created representing the marble A and J specimens in the BD and CCNBD tests. Each model comprised ~20,000 circular basic particles with a 0.2–0.3 mm radius. Since even a single mineral crystal is anisotropic along different atomic lattices [32], we set the strength and deformation parameters of the soft-bonded contacts to follow the Weibull distribution, and the shape parameter representing heterogeneity was set to 3 and 5 for marble A and J, respectively. A small value for the shape parameter indicates strong heterogeneity.
Before conducting the simulation of the CCNBD tests, the microscale parameters of the particles and contacts required iteratively calibrating through trial and error referring to the BD test results and the previous work [27,33,34]. The GBMs of the BD specimens were positioned between two stiff walls representing the loading end and platform of a compression machine, and the walls moved toward each other at the same constant velocity to result in a quasi-static loading rate. The calibration completes until the simulated load–displacement curves and crack morphology fit well with the observations in the BD tests (Figure 9). The calibrated microscale parameters (Table 2) were used to simulate the marble specimens in the CCNBD test.

4.3. Simulation Results

In the numerical simulations of the CCNBD test, when the applied load reached a certain level, micro-cracking was initiated near the notched crack tips of the specimens (Figure 10a and Figure 11a). At the peak loads (Pmax) of the marble A and J specimens, the micro-cracks around the crack tips tended to coalesce to form new macroscopic cracks (Figure 10b and Figure 11b). After that, the notched crack propagated dramatically, which caused a rapid post-peak drop in the applied load and the specimen to rupture (Figure 10c and Figure 11c). Therefore, the preparation stage for dramatic propagation of a notched crack can be defined as from the initiation of micro-cracking to the reaching of the peak load, during which micro-cracks generate to develop the FPZ. The FPZ is the area near the crack tip with a dense micro-crack distribution when the peak load is reached that stays in the critical state of macroscopic rupture. As mentioned in Section 2, the crack tolerance is characterized by the size of the FPZ.
As the simulation results (Figure 10b and Figure 11b) showed, the micro-crack density at the vicinity of notched crack tip was especially high in the whole specimen, owing to the nonlinear deformation brought by stress concentration. The FPZ was assumed as a tip-centered circle with radius of rc [14] that covered the area with high micro-crack density. Therefore, with increasing distance (radius) Rf from the tip and the diminishing intensity of the stress concentration, the deformation transitions from nonlinear inside the FPZ to quasilinear outside the FPZ, and thus the micro-crack density outside the FPZ declined to the background density of the rock [35].
To simplify the analysis, we assumed that the distributions of micro-crack inside and outside the FPZ are uniform but have different density. On the basis of this, the total micro-crack number N within a certain circular statistical range with Rf radius can be formulated as
N R = π ,         R f < r c π r c 2 d d 0 + π R f 2 d 0 ,         R f r c
where d and d0 are the average micro-crack density inside and outside the FPZ, respectively. N displayed a positively correlation with Rf; however, the curves of N deflected when Rf increased to rc that defined the boundary of the FPZ (Figure 12). This is because the micro-crack density d inside the FPZ can be up to ~15 times as great as the background density d0 [36]; the increasing rate of N will decelerate once the statistical range extend outside the FPZ. Such a deflection became more identifiable with the increasing ratio of d/d0. Thus, the crack tolerance of the specimens can be determined by the radius corresponding to the deflection point.
As the above calculation predicts, the micro-crack number N in Figure 10b and Figure 11b for the marble A and J specimens increased with Rf, and the N–Rf curves deflected at radii of 7 and 4 mm, respectively (Figure 13), which were closed to the mean values of the crack tolerance of A (~8 mm) and J (~3 mm) measured by experiments. These results demonstrate that calculating the crack tolerance using the tensile fracture toughness and average tensile strength in Equation (3) leads to reliable results. The simulations also showed that marble A had a greater crack tolerance than marble J, and the FPZ of the former contained more micro-cracks than that of the latter preceding specimen rupture (Figure 13c,d). Correspondingly, the GBM of marble A generates 1473 micro-cracks in the whole loading process, more than that of marble J (1404). These results showed that a greater crack tolerance can represent a stronger ability to generate micro-cracks.
The strike angle of micro-cracks (Figure 14) for the marble J mainly distributed in the range of 70–100°, which was narrower than that, 60–110°, for the marble A. This result suggests that more micro-cracks deviated from the loading direction (90°) in the marble A. The coalescence of such micro-cracks with various strike angles resulted in macroscopic cracks that propagated along irregular even branched paths and radiated from the notched crack tips (Figure 5 and Figure 10). Otherwise, the macro-cracks will develop primarily parallel to the loading direction, and thus their morphologies were less complex, as the marble J specimen showed. Therefore, these observations confirmed the assumption that fracability can be characterized by the crack tolerance.

5. Discussion

5.1. Effect of Grain Size on Crack Tolerance

The grain size greatly influences the cracking behavior of rocks [37]. Regarding rock consisting of grains with various size, the grain size greatly differs across parts of the rock specimen, so does the microscopic strength, which enhances the rock heterogeneity. This is why the heterogeneity index (shape parameter for the Weibull distribution) of marble A was set as smaller than that for marble J in the GBMs.
The microscopic strength in different parts of a strongly heterogeneous rock specimen can distribute in a wide range, so a small increment in the stress near a crack tip can easily cause micro-cracking within such rocks. Therefore, the initiation of micro-cracking was earlier, i.e., corresponding to a smaller ratio between the applied load and the peak load, in marble A than in marble J (Figure 10 and Figure 11). However, cracking can also be arrested easily because it probably encounters stronger local parts soon. Therefore, the rupture of strongly heterogeneous specimens will not occur until there are adequate micro-cracks to create FPZs and macroscopic cracks. In summary, strong heterogeneity strengthens the crack tolerance of rocks.
For rock specimens with a homogeneous grain size distribution, its microscopic strength in different parts can be generally closed to a certain level. Thus, only a few micro-cracks arise before the stress near the notched crack tips reaches that strength level. Once the strength is reached, the cracks propagate dramatically, which causes a rapid rupture. Macroscopic cracks spanning across the specimen form nearly instantaneously following the generation of a small FPZ. For these reasons, relatively homogeneous marble J had a smaller crack tolerance than marble A.

5.2. Effects of the Bedding Orientation and Environmental Temperature

With a vertical bedding orientation, micro-cracking naturally initiates within bedding planes in front of crack tips and propagates along them because the tensile strength of shale bedding planes is usually much smaller than that of layers between planes. Therefore, the tensile strength and fracture toughness were relatively low. Since the micro-cracks are limited to thin bedding planes, the corresponding crack tolerance is also small. With a horizontal bedding orientation, layers comprising various mineral grains make the distribution of the microscopic tensile strength near each tip more heterogeneous, so micro-cracking may be scattered among the layers. The behavior of micro-cracks can be complex when they cross bedding planes: cracks may branch along bedding planes and result in curved macroscopic cracks. Although hydraulic fracturing involves with many factors affecting the interaction between cracks and bedding planes [38], it is generally recognized that a crack propagating along a bedding plane is the most unfavorable situation for generating complex fracture networks [39,40].
Thermal treatment of sandstone leads to dehydration and the thermal expansion of minerals, which promotes the brittle–ductile transition of minerals [41]. Dehydration occurs at ~100 °C, which is when absorbed water escapes from the mineral surface [42], and enhances the friction and bonding strength between minerals. Thermal expansion takes effect when the temperature exceeds 100 °C and closes preexisting micro-cracks in rocks [43], which enhances the tensile strength of mineral grains and boundaries within a certain temperature range. Owing to these effects arising exceeding 100 °C, the sandstone at 125 °C had higher tensile strength and fracture toughness than that at 75 °C (Figure 4c). However, the dehydration and thermal expansion of minerals reduce the structural heterogeneity and rock fracability. In addition, fracability also weakens when sandstone transitions from brittle to ductile [44]. These observations suggest that increasing the temperature from 25 °C to 125 °C should be unfavorable for crack generation in sandstone. Such a negative effect of temperature was also observed in Longmaxi shale [45], a commercial shale gas reservoir in Chongqing, China. Considering that the downhole temperature, especially in deep and geothermal wells, can be much higher than the surface temperature, more attention should be given to the effect of temperature over a broader range on rock fracability.

5.3. Implications in Hydraulic Fracturing

The FPZ indicates nonlinear deformation (i.e., micro-crack generation) within a rock, which originates from rock heterogeneity. Nonlinear deformation diminishes and transitions into linear deformation with decreasing rock heterogeneity, which would reduce the FPZ size (i.e., crack tolerance). Complex cracks barely form in the absence of micro-cracks and their interactions. Because the formation of the FPZ is intrinsic to heterogeneous rock, the association between the FPZ size and heterogeneity is the physical basis for the crack tolerance as an evaluation index of fracability. The crack tolerance can reflect the effects of structural and environmental factors because they influence rock heterogeneity [46,47].
A large crack tolerance indicates many micro-cracks within the FPZ, which would tend to cause an irregular morphology of macroscopic cracks and wider zones of micro-cracks along both sides of cracks. These characteristics allow more cracks and pores in the rocks to connect with the main cracks during crack propagation driven by fluid. This increases the volume of cracks, which enhances the fluid conductivity of rock to exploit oil, gas, and geothermal resources.
The present study mainly provided laboratory observations on the effects of three factors to support the rationality of crack tolerance in reservoir fracability evaluation. However, the complexity of crack networks is dependent on various factors. Further investigations that consider other environmental and structural effects (e.g., magnitude and direction of crustal stress) are required to test this fracability index. On the other hand, the reliability and practical applicability of this index need to be further tested using rock specimens collected from reservoirs being exploited.

6. Conclusions

Aiming to the challenge of lacking a reliable index of fracability evaluation for hydraulic fracturing, we suggest that the crack tolerance, i.e., the maximum radius of the FPZ, may be used to characterize the fracability.
Crack tolerance originates from the unique heterogeneity of rock and inherent rules of crack generation, and thus it has a clear physical meaning and firm mechanical basis.
This index can be conveniently quantified in the laboratory using BD and CCNBD tests. We showed that the crack tolerance is positively correlated with the grain size heterogeneity and is negatively correlated with the environmental temperature (25–125 °C). The crack tolerance of shale specimens was greater with a horizontal bedding orientation than with a vertical orientation. Thus, crack tolerance can well characterize the effects of certain rock properties and environmental factors.
In summary, crack tolerance is promising in serving as a reliable evaluation index of rock fracability in terms of rationality and practicability, which has significant engineering implications for efficient hydraulic fracturing.

Author Contributions

Conceptualization, J.N. and H.C.; methodology, J.N., M.Z. and H.C.; validation, J.N.; formal analysis, H.C.; investigation, H.C., M.Z. and J.N.; resources, J.N.; writing—original draft preparation, H.C.; writing—review and editing, J.N., M.Z. and H.C.; visualization, H.C. and M.Z.; funding acquisition, J.N. and H.C. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China, grant numbers 42107184 and 41602330.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

All the data required to evaluate the conclusions of this study are present in the paper. The authors will provide additional data related to this paper upon request.

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclature

List of symbols
a0Initial chevron notched crack length
a1Final chevron notched crack length
ASTMAmerican Society of Testing Materials
BThickness of CCNBD specimen
BThickness of BD specimen
BDBrazilian disk
CCNBDCracked chevron notched Brazilian disk
dAverage micro-crack density inside FPZ
DDiameter of CCNBD specimen
DDiameter of BD specimen
d0Average micro-crack density outside FPZ
DEMDiscrete element method
FPZFracture process zone
GBMGrain-based model
ISRMInternational Society for Rock Mechanics
KITensile stress intensity factor
KICTensile fracture toughness
NTotal micro-crack number
PFCParticle flow code
PmaxPeak applied axial load in CCNBD test
PmaxPeak applied axial load in BD test
RRadius of CCNBD specimen
rPolar radius
rcRadius of circular FPZ
RfRadius of circular statistical range
RsSaw radius
uGeometric constants
vGeometric constants
XRDX-ray diffraction
Y*minCritical dimensionless stress intensity
Greek symbols
α0Dimensionless parameter
α1Dimensionless parameter
αBDimensionless parameter
θPolar angle
σ1Maximum principal stress
σ2Intermediate principal stress
σtTensile strength

References

  1. Mullen, M.J.; Enderlin, M.B. Fracability index—More than rock properties. In Proceedings of the SPE Annual Technical Conference and Exhibition, San Antonio, TX, USA, 8–10 October 2012; pp. 1–10. [Google Scholar]
  2. Wang, S.; Li, X.; Yao, J.; Gong, F.; Li, X.; Du, K.; Tao, M.; Huang, L.; Du, S. Experimental investigation of rock breakage by a conical pick and its application to non-explosive mechanized mining in deep hard rock. Int. J. Rock Mech. Min. Sci. 2019, 122, 104063. [Google Scholar] [CrossRef]
  3. Meng, F.; Wong, L.N.Y.; Zhou, H. Rock brittleness indices and their applications to different fields of rock engineering: A review. J. Rock Mech. Geotech. Eng. 2021, 13, 221–247. [Google Scholar] [CrossRef]
  4. Zhang, D.; Ranjith, P.G.; Perera, M.S.A. The brittleness indices used in rock mechanics and their application in shale hydraulic fracturing: A review. J. Pet. Sci. Eng. 2016, 143, 158–170. [Google Scholar] [CrossRef]
  5. Jarvie, D.M.; Hill, R.J.; Ruble, T.E.; Pollastro, R.M. Unconventional shale-gas systems: The Mississippian Barnett Shale of north-central Texas as one model for thermogenic shale-gas assessment. AAPG Bull. 2007, 91, 475–499. [Google Scholar] [CrossRef]
  6. Rickman, R.; Mullen, M.J.; Petre, J.E.; Grieser, W.V.; Kundert, D. A practical use of shale petrophysics for stimulation design optimization: All shale plays are not clones of the Barnett shale. In Proceedings of the SPE Annual Technical Conference and Exhibition, Denver, CO, USA, 21–24 September 2008; pp. 1–9. [Google Scholar]
  7. Bai, M. Why are brittleness and fracability not equivalent in designing hydraulic fracturing in tight shale gas reservoirs. Petroleum 2016, 2, 1–19. [Google Scholar] [CrossRef] [Green Version]
  8. Nasehi, M.J.; Mortazavi, A. Effects of in-situ stress regime and intact rock strength parameters on the hydraulic fracturing. J. Pet. Sci. Eng. 2013, 108, 211–221. [Google Scholar] [CrossRef]
  9. Hucka, V.; Das, B. Brittleness determination of rocks by different methods. Int. J. Rock Mech. Min. Sci. Geomech. Abstr. 1974, 11, 389–392. [Google Scholar] [CrossRef]
  10. Meng, F.; Zhou, H.; Zhang, C.; Xu, R.; Lu, J. Evaluation methodology of brittleness of rock based on post-peak stress–strain curves. Rock Mech. Rock Eng. 2015, 48, 1787–1805. [Google Scholar] [CrossRef]
  11. Xia, Y.; Zhou, H.; Zhang, C.; He, S.; Gao, Y.; Wang, P. The evaluation of rock brittleness and its application: A review study. Eur. J. Environ. Civ. Eng. 2022, 26, 239–279. [Google Scholar] [CrossRef]
  12. Jin, X.; Shah, S.N.; Roegiers, J.-C.; Zhang, B. Fracability evaluation in shale reservoirs—An integrated petrophysics and geomechanics approach. In Proceedings of the SPE Hydraulic Fracturing Technology Conference, The Woodlands, TX, USA, 4–6 February 2014; pp. 1–14. [Google Scholar]
  13. Wang, S.; Sun, L.; Li, X.; Zhou, J.; Du, K.; Wang, S.; Khandelwal, M. Experimental investigation and theoretical analysis of indentations on cuboid hard rock using a conical pick under uniaxial lateral stress. Geomech. Geophys. Geo Energ. Geo Resour. 2022, 8, 34. [Google Scholar] [CrossRef]
  14. Atkinson, B.K. Fracture Mechanics of Rock; Academic Press: London, UK, 1987. [Google Scholar]
  15. Lawn, B.R.; Wilshaw, T.R. Fracture of Brittle Solids; Cambridge University Press: New York, NY, USA, 1975. [Google Scholar]
  16. Saleur, H.; Sammis, C.G.; Sornette, D. Discrete scale invariance, complex fractal dimensions, and log-periodic fluctuations in seismicity. J. Geophys. Res. Solid Earth 1996, 101, 17661–17677. [Google Scholar] [CrossRef]
  17. Ghamgosar, M.; Bahaaddini, M.; Erarslan, N.; Williams, D.J. A new experimental approach to quantify microfractures in the Fracture Process Zone (FPZ) under various loading conditions. Eng. Geol. 2021, 283, 106024. [Google Scholar] [CrossRef]
  18. Ghamgosar, M.; Erarslan, N.; Williams, D.J. Experimental investigation of fracture process zone in rocks damaged under cyclic loadings. Exp. Mech 2017, 57, 97–113. [Google Scholar] [CrossRef]
  19. Guha Roy, D.; Singh, T.N.; Kodikara, J. Influence of joint anisotropy on the fracturing behavior of a sedimentary rock. Eng. Geol. 2017, 228, 224–237. [Google Scholar] [CrossRef]
  20. Zietlow, W.K.; Labuz, J.F. Measurement of the intrinsic process zone in rock using acoustic emission. Int. J. Rock Mech. Min. Sci. 1998, 35, 291–299. [Google Scholar] [CrossRef]
  21. Erarslan, N. Microstructural investigation of subcritical crack propagation and Fracture Process Zone (FPZ) by the reduction of rock fracture toughness under cyclic loading. Eng. Geol. 2016, 208, 181–190. [Google Scholar] [CrossRef]
  22. Han, Z.; Zhang, L.; Azzam, R.; Zhou, J.; Wang, S. A statistical index indicating the degree and mechanical effects of grain size heterogeneity in rocks. Eng. Geol. 2021, 293, 106292. [Google Scholar] [CrossRef]
  23. Fowell, R.J.; Xu, C.; Dowd, P.A. An update on the fracture toughness testing methods related to the cracked chevron-notched Brazilian disk (CCNBD) specimen. Pure Appl. Geophys. 2006, 163, 1047–1057. [Google Scholar] [CrossRef]
  24. American Society of Testing Materials. Standard Test Method for Splitting Tensile Strength of Intact Rock Core Specimens; American Society of Testing Materials: West Conshohocken, PA, USA, 2016; pp. 1–5. [Google Scholar] [CrossRef]
  25. Wong, L.N.Y.; Einstein, H.H. Crack coalescence in molded gypsum and Carrara marble: Part 2—Microscopic observations and interpretation. Rock Mech. Rock Eng. 2009, 42, 513–545. [Google Scholar] [CrossRef]
  26. Wang, H. Experimental Study of Shale Rock Mechanical Properties and Volume Fracturing. Ph.D. Thesis, University of Chinese Academy of Sciences, Beijing, China, 2015. [Google Scholar]
  27. Chen, H.; Niu, J.; Zhai, M. Characteristics of the fracture process zone for reservoir rock with various heterogeneity. Energies 2022, 15, 8332. [Google Scholar] [CrossRef]
  28. Potyondy, D.O.; Cundall, P.A. A bonded-particle model for rock. Int. J. Rock Mech. Min. Sci. 2004, 41, 1329–1364. [Google Scholar] [CrossRef]
  29. Potyondy, D.O. A grain-based model for rock: Approaching the true microstructure. In Proceedings of the Rock Mechanics in the Nordic Countries 2010, Kongsberg, Norway, 9–12 June 2010; pp. 225–234. [Google Scholar]
  30. Itasca Consulting Group Inc. PFC2D—Particle Flow Code in 2 Dimensions; Itasca Consulting Group Inc: Minneapolis, MN, USA, 2019. [Google Scholar]
  31. Pierce, M.; Cundall, P.; Potyondy, D.; Ivars, D.M. A synthetic rock mass model for jointed rock. In Proceedings of the 1st Canada—U.S. Rock Mechanics Symposium, Vancouver, BC, Canada, 27–31 May 2007. [Google Scholar]
  32. Whitney, D.L.; Broz, M.; Cook, R.F. Hardness, toughness, and modulus of some common metamorphic minerals. Am. Miner. 2007, 92, 281–288. [Google Scholar] [CrossRef]
  33. Wang, S.; Zhou, J.; Zhang, L.; Han, Z.; Zhang, F. Parameter studies on the mineral boundary strength influencing the fracturing of the crystalline rock based on a novel Grain-Based Model. Eng. Fract. Mech. 2021, 241, 107388. [Google Scholar] [CrossRef]
  34. Zhou, J.; Lan, H.; Zhang, L.; Yang, D.; Song, J.; Wang, S. Novel grain-based model for simulation of brittle failure of Alxa porphyritic granite. Eng. Geol. 2019, 251, 100–114. [Google Scholar] [CrossRef]
  35. Moore, D.E.; Lockner, D.A. The role of microcracking in shear-fracture propagation in granite. J. Struct. Geol. 1995, 17, 95–114. [Google Scholar] [CrossRef]
  36. Janssen, C.; Wagner, F.; Zang, A.; Dresen, G. Fracture process zone in granite: A microstructural analysis. Int. J. Earth. Sci. 2001, 90, 46–59. [Google Scholar] [CrossRef]
  37. Eberhardt, E.; Stimpson, B.; Stead, D. Effects of grain size on the initiation and propagation thresholds of stress-induced brittle fractures. Rock Mech. Rock Eng. 1999, 32, 81–99. [Google Scholar] [CrossRef]
  38. Zhang, Q.; Zhang, X.; Sun, W. A review of laboratory studies and theoretical analysis for the interaction mode between induced hydraulic fractures and pre-existing fractures. J. Nat. Gas. Sci. Eng. 2021, 86, 103719. [Google Scholar] [CrossRef]
  39. Lin, C.; He, J.; Li, X.; Wan, X.; Zheng, B. An experimental investigation into the effects of the anisotropy of shale on hydraulic fracture propagation. Rock Mech. Rock Eng. 2017, 50, 543–554. [Google Scholar] [CrossRef]
  40. Zhao, Y.; Zhang, Y.; Wang, C.; Liu, Q. Hydraulic fracturing characteristics and evaluation of fracturing effectiveness under different anisotropic angles and injection rates: An experimental investigation in absence of confining pressure. J. Nat. Gas. Sci. Eng. 2022, 97, 104343. [Google Scholar] [CrossRef]
  41. Sha, S.; Rong, G.; Tan, J.; He, R.; Li, B. Tensile strength and brittleness of sandstone and granite after high-temperature treatment: A review. Arab. J. Geosci. 2020, 13, 598. [Google Scholar] [CrossRef]
  42. Zhang, W.; Sun, Q.; Hao, S.; Geng, J.; Lv, C. Experimental study on the variation of physical and mechanical properties of rock after high temperature treatment. Appl. Eng. 2016, 98, 1297–1304. [Google Scholar] [CrossRef]
  43. Sirdesai, N.N.; Mahanta, B.; Ranjith, P.G.; Singh, T.N. Effects of thermal treatment on physico-morphological properties of Indian fine-grained sandstone. Bull. Eng. Geol. Env. 2019, 78, 883–897. [Google Scholar] [CrossRef]
  44. Zuo, J.; Zhou, H.; Xie, H. Fracture characteristics of sandstone under thermal effects. Eng. Mech. 2008, 25, 124–130. [Google Scholar]
  45. Niu, J. Research on the Fracture Mechanics Index of Rock Fracability Evaluation. Ph.D. Thesis, University of Chinese Academy of Sciences, Beijing, China, 2015. [Google Scholar]
  46. Chen, B.; Barboza, B.R.; Sun, Y.; Bai, J.; Thomas, H.R.; Dutko, M.; Cottrell, M.; Li, C. A review of hydraulic fracturing simulation. Arch. Comput. Method E. 2022, 29, 1–58. [Google Scholar] [CrossRef]
  47. Liu, P.; Ju, Y.; Gao, F.; Ranjith, P.G.; Zhang, Q. CT identification and fractal characterization of 3-D propagation and distribution of hydrofracturing cracks in low-permeability heterogeneous rocks. J. Geophys. Res. Sol. Ea. 2018, 123, 2156–2173. [Google Scholar] [CrossRef]
Figure 1. (ad) Formation of a fracture process zone (FPZ), shaded in red, and the propagation of macro fractures under the tensile stress σy (denoted by red arrows). (e) Schematic of an FPZ with the size rc.
Figure 1. (ad) Formation of a fracture process zone (FPZ), shaded in red, and the propagation of macro fractures under the tensile stress σy (denoted by red arrows). (e) Schematic of an FPZ with the size rc.
Materials 15 08485 g001
Figure 2. (a) Orthographic and (b) side profiles of a marble CCNBD specimen; (c) valid range for dimensionless parameters α1 and αB (outlined in gray) [23] and the distribution of parameter values for all of the prepared CCNBD specimens. Geometric parameters: diameter D = 75 mm, radius R = 37.5 mm, thickness B = 30 mm, saw radius Rs = 25 mm, initial chevron notched crack length a0 = 8.45 mm, and final chevron notched crack length a1 = 23.5 mm.
Figure 2. (a) Orthographic and (b) side profiles of a marble CCNBD specimen; (c) valid range for dimensionless parameters α1 and αB (outlined in gray) [23] and the distribution of parameter values for all of the prepared CCNBD specimens. Geometric parameters: diameter D = 75 mm, radius R = 37.5 mm, thickness B = 30 mm, saw radius Rs = 25 mm, initial chevron notched crack length a0 = 8.45 mm, and final chevron notched crack length a1 = 23.5 mm.
Materials 15 08485 g002
Figure 3. (a) CCNBD test on a marble A specimen at room temperature; (b) BD test on a sandstone specimen in the chamber at 125 °C; (c) MTS servo-control testing machine, SNAS GDS-300 environmental chamber, and (d) WK650 temperature controller.
Figure 3. (a) CCNBD test on a marble A specimen at room temperature; (b) BD test on a sandstone specimen in the chamber at 125 °C; (c) MTS servo-control testing machine, SNAS GDS-300 environmental chamber, and (d) WK650 temperature controller.
Materials 15 08485 g003
Figure 4. Mean tensile strength (brown), tensile fracture toughness (green), and crack tolerance (sienna) of the (a) marble A (diamond) and J (triangle), (b) shale (circle) with horizontal and vertical orientations and (c) sandstone (square).
Figure 4. Mean tensile strength (brown), tensile fracture toughness (green), and crack tolerance (sienna) of the (a) marble A (diamond) and J (triangle), (b) shale (circle) with horizontal and vertical orientations and (c) sandstone (square).
Materials 15 08485 g004
Figure 5. Crack morphologies of (a,b) marble A and (c,d) J. The red dashed ellipses in (a) denote white patches around the notched crack tips.
Figure 5. Crack morphologies of (a,b) marble A and (c,d) J. The red dashed ellipses in (a) denote white patches around the notched crack tips.
Materials 15 08485 g005
Figure 6. Crack morphologies of the shale CCNBD specimens with (a,b) horizontal and (c,d) vertical bedding orientations.
Figure 6. Crack morphologies of the shale CCNBD specimens with (a,b) horizontal and (c,d) vertical bedding orientations.
Materials 15 08485 g006
Figure 7. Crack morphologies of the sandstone specimens in the CCNBD tests at (a,b) 25 °C, (c,d) 75 °C, and (e,f) 125 °C.
Figure 7. Crack morphologies of the sandstone specimens in the CCNBD tests at (a,b) 25 °C, (c,d) 75 °C, and (e,f) 125 °C.
Materials 15 08485 g007
Figure 8. (a) Grain-based model for marble and (b) the implement of inter-grain soft-bonded and intra-grain smooth-joint contacts. Schematic of bonded and unbonded state for (c) soft-bonded and (d) smooth-joint contacts. (b) is the zoomed-in view of the dashed box in (a).
Figure 8. (a) Grain-based model for marble and (b) the implement of inter-grain soft-bonded and intra-grain smooth-joint contacts. Schematic of bonded and unbonded state for (c) soft-bonded and (d) smooth-joint contacts. (b) is the zoomed-in view of the dashed box in (a).
Materials 15 08485 g008
Figure 9. (a) Load–displacement relationship of experimental and numerical results for marble A and J in the BD test. (b,c) Crack morphologies of the marble A and J specimens in the numerical simulation of BD test and the corresponding experiment results.
Figure 9. (a) Load–displacement relationship of experimental and numerical results for marble A and J in the BD test. (b,c) Crack morphologies of the marble A and J specimens in the numerical simulation of BD test and the corresponding experiment results.
Materials 15 08485 g009
Figure 10. Simulated cracking process of marble A CCNBD specimen. (a) Cracking initiation at notched crack tips (~66% Pmax); (b) reaching Pmax; (c) specimen rupture.
Figure 10. Simulated cracking process of marble A CCNBD specimen. (a) Cracking initiation at notched crack tips (~66% Pmax); (b) reaching Pmax; (c) specimen rupture.
Materials 15 08485 g010
Figure 11. Simulated cracking process of marble J CCNBD specimen. (a) Cracking initiation at notched crack tips (~70% Pmax); (b) reaching Pmax; (c) specimen rupture.
Figure 11. Simulated cracking process of marble J CCNBD specimen. (a) Cracking initiation at notched crack tips (~70% Pmax); (b) reaching Pmax; (c) specimen rupture.
Materials 15 08485 g011
Figure 12. Variation trend of total micro-crack number inside tip-centered circles with various radii Rf. The three curves correspond to cases under different ratios of d/d0.
Figure 12. Variation trend of total micro-crack number inside tip-centered circles with various radii Rf. The three curves correspond to cases under different ratios of d/d0.
Materials 15 08485 g012
Figure 13. Micro-crack distributions near the notched crack tips in the (a) marble A and (b) J specimens at the peak applied load in the CCNBD simulation. (c,d) Total number of micro-cracks in tip-centered circles with different radii. The red and blue curves in (a,b) indicate intra- and inter-grain micro-cracks. The dashed lines in (c,d) are fitting lines for the distribution of the solid circles.
Figure 13. Micro-crack distributions near the notched crack tips in the (a) marble A and (b) J specimens at the peak applied load in the CCNBD simulation. (c,d) Total number of micro-cracks in tip-centered circles with different radii. The red and blue curves in (a,b) indicate intra- and inter-grain micro-cracks. The dashed lines in (c,d) are fitting lines for the distribution of the solid circles.
Materials 15 08485 g013
Figure 14. Strike angle distribution of micro-cracks in the CCNBD simulations for (a) marble A and (b) J, corresponding to Figure 10b and Figure 11b. The dashed boxes define the main distribution ranges.
Figure 14. Strike angle distribution of micro-cracks in the CCNBD simulations for (a) marble A and (b) J, corresponding to Figure 10b and Figure 11b. The dashed boxes define the main distribution ranges.
Materials 15 08485 g014
Table 1. Microstructures of the rock specimens.
Table 1. Microstructures of the rock specimens.
Microscopy ObservationQualitative DescriptionAverage Size (mm)Grain Size Distribution
Materials 15 08485 i001Medium-coarse grains3.43–5 mm: 70%
1–3 mm: 30%
Materials 15 08485 i002Medium grains1.51–2 mm
Materials 15 08485 i003Fine grains<0.1/
Materials 15 08485 i004Coarse grains21–3 mm
Table 2. Calibrated microscale parameters of marble A and J specimens.
Table 2. Calibrated microscale parameters of marble A and J specimens.
ElementsMicroscale ParametersMarble AMarble J
ParticlesDensity (kg/m3)26902690
Effective modulus (GPa)7070
Normal to shear stiffness ratio1.51.5
Friction coefficient0.60.6
Soft-bonded intra-grain contactEffective modulus (GPa)15.010.0
Normal to shear stiffness ratio2.02.0
Friction coefficient0.60.6
Tensile strength (MPa)28.018.0
Cohesion (MPa)112.072.0
Friction angle (°)4545
Softening factor0.10.1
Softening tensile strength factor0.70.7
Smooth-joint inter-grain contactTensile strength coefficient0.30.3
Cohesion coefficient0.80.8
Friction angle coefficient0.60.6
Friction adjustment coefficient0.50.5
Normal stiffness coefficient0.80.8
Shear stiffness coefficient0.150.15
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Chen, H.; Niu, J.; Zhai, M. An Evaluation Index of Fracability for Reservoir Rocks Based on Fracture Process Zone. Materials 2022, 15, 8485. https://doi.org/10.3390/ma15238485

AMA Style

Chen H, Niu J, Zhai M. An Evaluation Index of Fracability for Reservoir Rocks Based on Fracture Process Zone. Materials. 2022; 15(23):8485. https://doi.org/10.3390/ma15238485

Chicago/Turabian Style

Chen, Hongran, Jingrui Niu, and Mengyang Zhai. 2022. "An Evaluation Index of Fracability for Reservoir Rocks Based on Fracture Process Zone" Materials 15, no. 23: 8485. https://doi.org/10.3390/ma15238485

APA Style

Chen, H., Niu, J., & Zhai, M. (2022). An Evaluation Index of Fracability for Reservoir Rocks Based on Fracture Process Zone. Materials, 15(23), 8485. https://doi.org/10.3390/ma15238485

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