Statistical estimation of frictional coefficients of faults based on a structural dataset in the Tuy Hoa–Vung Tau Region, Viet Nam
South-Vietnam Geological Mapping Division, 200 Ly Chinh Thang Street, Ward 9, District 3, Ho Chi Minh City 700000, Vietnam
2 Faculty of Geology and Petroleum Engineering, Ho Chi Minh University of Technology (HCMUT), 268 Ly Thuong Kiet Street, District 10, Ho Chi Minh City 700000, Vietnam
3 Vietnam National University Ho Chi Minh City, Linh Trung Ward, Thu Duc District, Ho Chi Minh City 700000, Vietnam
4 PetroVietnam University, 762 Cach Mang Thang Tam Street, Long Toan Ward, Ba Ria City, Ba Ria – Vung Tau Province 790000, Vietnam
5 Department of Petroleum Engineering, Universitas Islam Riau, Jl. Kaharuddin Nasution 113, Pekanbaru Riau 28284, Indonesia
6 Faculty of Geology, University of Science, 227 Nguyen Van Cu Street, District 5, Ho Chi Minh City, 700000, Vietnam
* Corresponding author: email@example.com
Accepted: 18 March 2021
This study estimates the shear friction coefficient from shear friction angles for the prediction of slip tendencies in the Tuy Hoa–Vung Tau region of Southern Vietnam. A dataset consisting of measured data of 355 fractured planes, striations, and unconformities in coastal areas as well as 239 offshore faults was analyzed based on the principles of statistical probability. As a result, 138 friction angles for the onshore and offshore faults were calculated based on shear fracture conjugate pairs. The goodness-of-fit test was used to define the probability distribution of the friction angles, which had a normal distribution. The acceptable average of friction angles for the onshore region with a reliability of more than 95% were in the range of 25.8–31.5°, which corresponds to frictional coefficients of 0.48–0.61. The acceptable average friction angles for the offshore region were relatively low at 23–31°, which corresponds to a frictional coefficient of 0.42–0.60. Owing to the heterogeneity of the fault system, the median value (19.12°) should be used as the lowest threshold value for slipping faults at all conditions. The recommended applicable average friction angles are 28.65° and 27° for the onshore and offshore regions, respectively. The estimation of the frictional coefficients is highly reliable, and it can be applied to other subsurface resource exploitation projects within the study area.
© D. Linh et al., published by IFP Energies nouvelles, 2021
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
In the process of geological mapping and mineral exploration, the identification of striations in the dip direction and dip angle of the fault, which can lead to slip or overturn, is essential for determining the slip friction angle and the cohesion of the fault. Through field investigations, information can be gathered on fracture surface properties, including the continuity of length, surface roughness, width and depth, secondary mineral characteristics in the fracture, and the influence of water on mineral deposits .
The coefficient of friction, μ, is a characteristic of the first-order functional friction. The coefficient of static friction is also known as the coefficient of internal friction, basic friction, or residual friction. It is an essential parameter that controls fault motion but it is not related to stick-slip behavior. The coefficients of static friction (μs) and kinetic friction (μk) depend on the nature of the component surface of the contact rock (Fig. 1). The value of the coefficient of shear friction is less than that of static friction (μk < μs). While performing stability analyses, crustal rocks are often assumed to act as a Mohr–Coulomb material in its response to shear stress, whereby its shear strength is expressed in terms of the cohesion and angle of the slickenlines . Tectonic stress can cause large earthquakes and crustal faults to weaken or slide steadily. Byerlee’s frictional coefficient (0.6 < μ < 0.85) is widely applied to the naturally sliding plane of a fault , which can extend several kilometers in depth through the upper crust, from Earth’s surface toward the point of origin of an earthquake.
Static and kinetic friction.
The shear strength of a flat, clean fracture with no scratches and an adhesive force of zero can be determined using the static friction angle. The friction angle of the rock material is related to the size, shape, and composition of the particles exposed along the fault surface. Therefore, a fine-grained rock with a high mica content will have a small friction angle while a coarse-grained rock will have a large friction angle. The measured friction angles in sedimentary rocks range from 25° to 30°, while in intrusive and metamorphic rocks, they range from 30° to 35°. Table 1 lists the friction angle values of each rock type [4, 5].
Classification of the friction angle by rock composition.
Yukutake et al.  identified that accurately estimating the slip friction coefficient of a fault to evaluate the slip tendency using the frictional resistance of a fault is difficult. They determined the slip friction coefficient by analyzing the conjugate angles of faults based on the focal mechanisms of dip-slip earthquakes. Instead, they applied the coefficient of friction, according to  friction law. However, the results were not convincing when applied to specific areas to estimate slip tendencies and design water injection pressure thresholds. Thus, this method cannot be used in regions with different tectonic and petrological conditions.
At conditions with soft soil and pre-existing fault zones, reactivated faults exhibit a low frictional coefficient of ~0.3 or even 0.1, as observed from dip-slip earthquake focal mechanisms . The frictional reactivation angle for faults in coal mines ranges from 8° to 18°, . Tamura and Kambayashi  presented a fault analysis to estimate the frictional coefficient.
The tilt and direct shear tests are standard methods used to determine the friction angle [10–12]. In terms of geological form, e.g., in a heterogeneous and discontinuous granite in Bukit Pengorak (Malaysia), the basic friction angle of granite was determined to be 30° and the maximum friction angle changed from 47° to 80°, depending on the roughness of the surface .
Furthermore, the coefficient of shear friction depends on many factors, such as the filling material, degree of mineral variation, hydrothermal activity, and orientation of fracture relative to the modern stress field. A study on the Alpine Fault in New Zealand estimated the shear frictional coefficient to be 0.12–0.16 at temperatures from 25 °C to 210 °C, corresponding to an effective normal stress (σn′) of 31.2–93.6 MPa. Boulton et al.  showed that spaces between the fault planes were filled with clay particles, which significantly affected the shear friction.
Jeanne et al.  examined the distribution law of the friction coefficient. They applied the law of friction to evaluate the effect of fracture stability based on seismic data collected during carbon capture and storage. Simulation results showed that the distribution of the static and dynamic frictional coefficients (before failure) varied along the fault, with intensities of 1.88–5.88. The critical values of the frictional coefficients (minimum and maximum) along the fault zone could be predicted to control the degree of stress intensity accumulation before failure.
Assuming that the tectonic stress field and fault direction are stable, the shear stress on a fault plane (i) will not change at a particular moment for at least thousands of years. Therefore, faults slip when the friction angle exceeds the internal friction angle. In principle, the coefficient of shear friction is relatively stable because it reflects the inherent nature of the fault. Thus, we must consider the effect of stress to investigate fault reactivation.
A pioneering method proposed by Anderson  to determine the paleostress direction assumes that conjugate pairs of shear fractures evolved according to the Mohr–Coulomb standard. Pairs of conjugate fractures are two fault planes or microscopic fractures that are generated simultaneously in the newly formed structure. The dihedral angles between the two fault planes vary from 30° to 89° for brittle deformation and are > 90° for ductile deformation. Depending on the nature of the mechanical and physical properties of the brittle deformation, i.e., brittle-ductile or ductile, this angle varies from 30° to 150°. Usually, this dihedral angle in brittle granite ranges from 55° to 62° and can sometimes be close to 90°. When the friction angle is more extensive, this corresponding angle approaches 90°. The conjugate fracture pairs determined using statistical analyses of the outcrop in this study exhibit shear and normal stresses at a point on the largest Mohr circle.
Assumptions of shear friction angles with different values can lead to errors in the design of drilling wells or when determining the water injection volume threshold for oil recovery. For instance, in the White Tiger oil field, offshore Vietnam, the shear friction angles of two wells (A and B) in a fractured granite reservoir changed from 48° to 57° at depths from 4000–4300 m . While examining the stress state of the oil and gas fields in the Cuu Long and Nam Con Son basins, Binh et al.  estimated that the frictional coefficient of a strike-slip normal fault in the basin was ~0.5, which corresponds to a shear friction angle of 26.6°. The lower threshold of the shear friction angle was also set at 16.7° for analysis using the TENSOR program. Thus, the accuracy of the shear friction angle is essential for forecasting the slip tendency and rock stability.
Based on a study of geothermal energy production, Gan and Elsworth  suggested that injecting fluid into a reservoir at high temperatures increased the pore-pressure and widened fracture apertures to facilitate fault motion, causing subterranean fault reactivation. Therefore, frictional properties influence the instability of faults and associated seismic events.
The internal friction angle depends on the rock type, porosity, and bonding. The weaker the rock, the lower the shear friction angle. Lama and Vutukuri  found that higher Young’s elastic modulus values of the rock result in higher shear friction angles. Son et al.  calculated the shear friction angle in the Oligocene sandstone reservoir (4012–4077 m) of the White Tiger oil field to be 28°–45° based on Poisson’s coefficient and the geophysical log Vp /Vs curve obtained from the well.
Sato  used a numerical method to estimate the slip friction coefficient based on the statistical probability distribution theory of points on the Mohr circle. The standard Mohr circle, where σ3 = 0, σ1 = 1, and σ2 = the stress ratio (from the inversion of striations or earthquake focal mechanisms), was used to obtain the shear and normal stresses along the sliding fault plane for any direction and dip angle in space within a specific tectonic stress field [23, 24]. Sato  argued that the points on a Mohr circle represent a density probability distribution close to the Mohr circle, which declines in the direction of the stress vector. Furthermore, by statistically obtaining the decrease in the vector direction, the slip friction coefficient can be calculated.
Energy projects are currently being developed in Tuy Hoa–Vung Tau region, wherein geothermal exploitation, oil and gas field exploitation, underground storage, and dam construction for hydroelectric power will be conducted in onshore and offshore regions. Thus, determining the shear friction coefficient thresholds from the shear friction angles is vital for predicting the slip tendencies, which can identify the impact of human activity or natural disasters. At present, no study has calculated the coefficients of friction in this region; instead, most studies use the assumed values of  or similar rock compositions from previous studies. The shear friction coefficient threshold is a primary parameter that can be used to assist managers in decision-making regarding the pre-feasibility and operational stages, as well as reducing costs by forecasting and minimizing risk. Human activity can change the stress balance and cause faults to slide or reactivate, thereby increasing the reservoir interconnections and enabling efficient resource exploitation.
The Tuy Hoa–Vung Tau region includes a territory in south-central Vietnam (Fig. 2) and its continental shelf, and is located in the western South China Sea (center of the Southeast Asian plate). This area is likely to have formed and developed on the Precambrian crystalline basement, followed by the development of an active continental margin setting during the late Mesozoic. The region was influenced by extension and lateral strike-slip movement in the early Cenozoic, followed by domain uplift with basalt volcanism during the late Cenozoic. From a modern tectonic perspective, the study region includes continental crust structural units (i.e., the Kon Tum, Da Lat, and Can Tho zones) as well as extensional, thinned, and differentiated transitional crust forming numerous oil-bearing Cenozoic basins.
Map of the geoblock, stress field, and human impacts on the environment.
From the Paleoproterozoic to the Quaternary, the study area experienced at least 14 other geodynamic settings, which can be reconstructed by their characteristic petrotectonic assemblages and are outlined in Table 2. The study area is located between three features (Fig. 2): the Indochina continental crust, Cuu Long-Con Son structural domain, and Paracel Islands-Borneo continental-transitional crust terrane.
In the early Cenozoic, the study area formed part of the southeastern interior of the Eurasian plate, which was influenced by two significant dynamic processes:
collision between India and Asia in the west, where the northwest-striking Mae Ping and Three Pagodas faults cut through the northwestern part of Indochina; and
the Red River fault, oriented in a northeast–southeast direction and the South China Sea divergent processes, which began in the Late Cretaceous [25, 26] and continued to spread at a maximum of 32–16 Ma, creating the new South China Sea oceanic plate .
The study area is a part of Indochina where basaltic eruption was accompanied by intensive massif uplift since the late Cenozoic. Here, extensional activity, sliding, and thermal subsidence meet a transitional continental crust massif. The formation mechanism resulted in two primary structural parts:
pre-Cenozoic crystallized basement, consisting of plutonic, volcanic, and sedimentary rocks and constituting the Indochina continental block; and
an overburden layer, consisting of mainly late Cenozoic terrestrial clastic sediment, interbedded with basalt rock, which filled in the Phu Khanh and Cuu Long basins, forming part of the transitional continental crust.
The coastal area from Tuy Hoa developed 21 major fault systems, including four faults at level I, seven faults at level II, and ten faults at level III. The faults are developed in four main directions: NW–SE, NE–SW, E–W, or sub-E–W, and N–S or sub-N–S (Tab. 3).
Morphological and dynamic characteristics of major faults in the coastal area from Tuy Hoa to Vung Tau.
Unconformable contacts are generally referred to as unconformities (Tab. 4).
The Cenozoic tectonic evolution in the study area and adjacent areas includes three main periods with six sub-periods (Tab. 5).
History of the Cenozoic tectonic evolution.
The compressive axis of the tectonic stress field construction has not changed substantially from the late Miocene to Quaternary, although it is different from the previous period in the geodynamic regime (i.e., the early to middle Miocene, when compression occurred in a NEN–SWS direction, and from the middle Eocene to early sub-Miocene, when compression was toward the sub-latitudes).
Measurements were taken of fractured planes, striations, and unconformities in outcrops of the Late Cretaceous granite in the Deo Ca complex, the late Miocene Basalt in the Dai Nga complex, and the Cretaceous rhyolites in the Nha Trang complex along the coastal areas from Tuy Hoa, Ke Ga, and Vung Tau to Phu Quy Island. A geo-compass and measuring tape were used to determine the morphology and scale of the features according to the following procedure:
measure the strike of fractures containing striations;
measure the dip direction and dip of the fractures or microfaults;
determine and measure the plunge and trend of the slip vector on the microfault surface;
determine the movement sense of the microfault or use indicators  (i.e., record the kinematic type as strike-slip, normal, reverse fault, etc.);
measure the rake angle or pitch of the microfault using the geo-compass or a degree ruler;
observe and find the angle of the conjugate fault system, which can be identified by the direction of the sliding pinch on the intersection point between two faults (this angle is twice the friction angle of the fault slip); and
record all measurements in a geological notebook.
The data collected from the steps described in the previous section were processed as follows:
the recorded data was processed to define the stress tensor at an exposed outcrop (from several hundred to thousands of meters in length), which has at least four microfaults measured that are distinct in orientation in the same rock type (e.g., granite, gabbro, sandstone, basalt, etc.);
the Win-Tensor  software was used to separate the structures into different stress tensors or different stress states based on the population of striations and slickenlines. The reduced stress tensor, corresponding to every sub-set microfault, was determined and categorized using four parameters: the stress ratio (ϕ) and the orientations (trend/plunge) of the three principal stress axes (largest, smallest, and intermediate);
MohrPlotter 3D  was used to analyze striations and slickenlines under the stress tensor that governed them, with a normalized tensor of σ1 = 1, σ2 = ϕ, and σ3 = 0, and to determine the shear stress and normal stress on each microfault at any orientation; and
the angle of the set of faults (represented by a cluster of small circles) was then measured to calculate the friction angle by measuring from the origin coordinates of the Mohr circle (see ). An error of ~6% was considered acceptable .
In total, we measured 355 fractured planes, striations, and unconformities in 196 outcrops along the coastal areas from Tuy Hoa, Ke Ga, and Vung Tau to Phu Quy Island. Residual friction angles of 34 granite samples in Cai Mep Port were measured using triaxial tests at depths of 90–173 m .
The offshore dataset included in situ stresses from 13 wells, which were measured in eight EXtended Leak-Off Tests (XLOT), and 23 Formation Micro-Imager (FMI). Twelve earthquake focal mechanisms were analyzed from 315 earthquakes detected between 2002 and May 2020. During this period, earthquakes mainly occurred along the NNE–SSW fault zone in and adjacent to the White Tiger–Rong fields (Fig. 2). We also interpreted 239 faults based on 14 000 km of regional offshore seismic data from the Cuu Long Basin .
Owing to the geological heterogeneity of fault systems, the coefficient of shear friction changes randomly. As the shear friction coefficient profiles were determined by the probability distribution function, a Monte Carlo simulation allows the use of these distributions and independent or dependent random variables . In this study, the normal distribution was used as the variable (x), with the mean and variance, σ 2, based on the following function:
The Mohr circle diagram demonstrates the relationship between the stress state and rheology of rocks in the crust and lithosphere. This diagram shows the stress ratio, which predicts the growth of new ruptures and the reactivation of existing faults.
The internal friction angle can be described as the shear strength due to friction according to the effective normal stress prescribed by the Mohr–Coulomb failure condition. The internal friction angle was determined using the shear and normal stresses at various levels (Fig. 3). According to the Mohr–Coulomb circle, the slope of the curve related to the shear strength and compressional pressure is the internal friction angle.
Mohr–Coulomb failure criterion in the effective stress domain. (a) Shear plane under triaxial compression. (b) Mohr circles and an attenuation contour cut linearly (tangential line of the Mohr circles), where is the angle of internal friction.
The following assumptions were made to calculate the coefficients of shear friction:
brittle deformation must comply with the Mohr–Coulomb failure standard;
the shear friction coefficient of an existing fault before the occurrence of critical sliding was constant for a specific tectonic regime. The faults could be strike-slip, normal, or reverse types;
the normalized Mohr circle was based on the stress ratio in the region corresponding to the tectonic regime. Simplified principal stress intensities were σ3 = 0, σ1 = 1, and σ2 = the stress ratio. For example, if the stress ratio was 0.7, then σ2 = 0.7 on the Mohr circle; and
the tensors of the principal stress axes (σ1, σ2, and σ3) were similar to the regional stress tensors.
This study analyzed the conjugate fault system and the ratio of normal to shear stress based on measured fault profiles to achieve a robust approach to understand the shear friction coefficient. The process of estimating the friction coefficient proceeded as follows:
Step 1: Conduct field measurements on sliding surfaces: dip directions, dip angles, and the directions of the motion vectors (Fig. 4).
Fracture parameters and stress state analysis.
Step 2: Classify the striation data according to petrographic compositions:
+ Late Cretaceous granite from the Deo Ca complex.
+ Intrusive basalt from the late Miocene Dai Nga complex, Pliocene to early Pleistocene Tuc Trung complex, and middle Pleistocene Xuan Loc complex in Ba Ria–Vung Tau province.
+ Cretaceous rhyolites from the Nha Trang complex.
Step 3: Build a contour diagram of fractures/striations using the Stereonet 10.0 software  to find the maximum values. Each maximum value corresponding to a specified dip direction and dip angle that represented a group of striations was identified.
Step 4: Estimate the conjugate angles between striation pairs (conjugate deformation). Two fractures/striations were considered to be conjugates when they formed an angle of 30°–89° . We used the Dips 8.0 software package to determine the striation pairs.
Step 5: Calculate the static friction angles based on the measured conjugate angles using the equation in :
Step 6: Based on the dataset’s friction angles, we performed statistical analyses to estimate the range of values.
The friction angles were determined using an appropriate probability distribution function in two steps:
Step 1: The type of random distribution was hypothesized by analyzing the available data. The proposed distribution assumed that the friction angle types conform to normal, log-normal, and beta distributions.
Step 2: The goodness-of-fit test was used to evaluate the distribution type by fitting the theoretical distribution to the hypothetical distribution at a specified significance level.
Twelve earthquake focal mechanisms were measured by the Vietnam Institute of Geophysics using the technique reported in  technique to analyze 24 conjugate angles based on conjugate fault pairs (Fig. 5). Thus, 24 friction angles were calculated (Tab. 6).
Analysis of the friction angles from 12 earthquake focal mechanisms.
Estimation of friction angles based on earthquake focal mechanisms.
Based on the Mohr circle failure standard, the calculated conjugate angles of fault pairs are represented by points on the Mohr circle. From the 11 fault systems in the White Tiger structure and its vicinity, 41 pairs of conjugate angles were determined and plotted on the Mohr circle. The ratios of shear to normal stresses in the fault profiles in the White Tiger field and adjacent areas were calculated based on data from 239 faults (Tab. 7).
Estimation of friction angles based on pairs of conjugated faults in the Oligocene sandstone and fractured granite reservoirs in the White Tiger field.
The profiles show that the friction angle of the fault system varied from 5° to 58°. The standard deviation and mean values were 12.23° and 27.65°, respectively. The dataset consisted of 73 samples and had a p-value of < 0.005, indicating statistical significance. Therefore, the probability analysis had a very high confidence interval, confirming that the average value of the friction angle was between 25.8° and 31.5°. A histogram of the data exhibited a normal distribution, confirming that the results were reliable (Fig. 6).
Summary of the statistical parameters for onshore friction angles.
The points on the graph were almost linear, indicating that the dataset had a normal distribution. In addition, the p-value was significant (α = 0.05), confirming that the data had a normal distribution (Fig. 7). For the onshore region, the low threshold was 24° (Fig. 6). Figure 7 shows a friction angle of 24° in the 20th percentile based on the diagram. Therefore, ~20% of the friction angles were equal to or less than this value while ~80% of the friction angles were greater than 24°.
Threshold of the onshore friction angle value.
The histogram of a dataset consisting of 65 calculated friction angles for a wide range of values (1–59°) had a normal distribution (Fig. 8). The p-value was 0.016, or less than α = 5%, indicating that the dataset was statistically significant. The mean and standard deviation were 15.95° and 26.93°, respectively. The average friction angle was in the range of 23–31° with high reliability.
Summary of the statistical parameters for offshore friction angles.
The normal probability distribution plot shows that the points are almost linear, indicating that the data have a normal distribution. In addition, the p-value for the test was 0.016, which was smaller than the significance level (α = 0.05). This confirms that the data have a normal distribution (Fig. 9).
Threshold of the offshore friction angle value.
According to  theory of friction, the initial friction angle was measured at 16.7°, which is the lower limit for sliding faults in all conditions. However, in the White Tiger field and adjacent areas, the lowest value of the median was 19.12°. Hence, this value should be used as the lower threshold (Fig. 8). Figure 9 shows that a slip friction angle of 19.12° was in the 38th percentile. Thus, ~38% of the friction angles were equal to or less than 19.12° and ~62% of the friction angles were greater than 19.12°.
One of the main challenges that occurs when performing reliability analysis is the determination of the range of acceptable friction angles. From this perspective, a high confidence level (95%) was selected as the threshold for the mean values in each surveyed area (Fig. 10):
For the onshore region: The threshold value ranged from 25.8° to 31.5° with an average of 28.65°.
For the offshore region: The threshold value ranged from 23° to 31° with an average of 27°. The initial friction angle chosen was the lowest median value (19.12°), corresponding to a frictional coefficient of 0.35, which is considered a low threshold for slipping faults in all conditions.
Threshold of the friction angle at a 95% confidence level.
The chart in Figure 10 shows that the friction angles in the offshore region were, on average, 1.65° lower than those of the onshore region. Thus, the equivalent offshore coefficient of friction was 0.51 while that of the onshore region was 0.55.
As estimating the friction angle of a sliding fault is difficult , previous studies have assumed a wide range for the coefficient of friction (0.6–0.85), equivalent to a friction angle of 31–40° . When compared with the results of this study, the lowest coefficient of friction obtained using the method proposed in  was equal to the upper threshold found for the offshore area of Vietnam. Therefore, applying the friction angle of  caused uncertainty in forecasting the sliding tendency in the Tuy Hoa–Vung Tau region.
Sibson  suggested that reactivated inverse faults have high slope angles, indicating that the optimal angle, θ* (equal to half of the conjugate angle), is 60° or greater. When the reactivated fault has a low slope angle, θ* is in the range of 20–30°, which reflects the frictional coefficient anomaly. This evidence is related to water injection and formation damage caused by human activity. The four earthquake focal mechanisms with conjugate fault plane angles in the range of 20–30° may have been caused by the impact of water injection, hydraulic fracturing, or abnormal pore pressures due to production activities.
The analysis in  suggests that any fault with a friction angle of ϕ < 20° will slip if the optimal θ* angle is 30–40°. If the friction angle of a fault is < 15°, the fault will slip within the direction range of 15° < θ < 60°. Faults with friction angles > 22° will not slip in any direction.
For in situ stress and borehole pressures in the Cuu Long Basin, friction angles of 48–57° obtained for wells A and B in the brittle fractured granite reservoir indicate deformation at depths of 4000–4300 m . This result is inaccurate because the value of 48°–57° was outside the allowable confidence interval, which led to incorrect well drilling design and determinations of the injected water volume threshold. This result was locally calculated and is not an expected value for most of the Cuu Long Basin oilfield.
For the stress field state at the wells in the Chu Long and Nam Con Son basins , a frictional coefficient of 0.5 was used to calculate the maximum horizontal stress for the normal strike-slip fault, which is equivalent to a friction angle of 26.6°. The assumed frictional coefficient was appropriate because the results were similar to those obtained in this study. A comparison of the Leak-Off Test (LOT) and hydraulic fracturing in several boreholes indicates that the maximum horizontal stress was 60–64 MPa. Thus, using a friction angle of 27° and a stress ratio of 0.51 was appropriate for the LOT test conditions.
By investigating the role of the friction coefficient in terms of the optimal slip angle of a fault, King et al.  found that a lower friction angle results in an optimal angle, θ, between the two slip plane faults closer to 90°. In contrast, a higher friction angle results in a smaller optimal angle, θ. The friction angle is important for assessing sensitivity to changes in the Coulomb stress and the effect that present-day tectonic stress fields have on strike-slip faults with changes in the coefficient of friction. The lowest friction angle is usually 8° to 18° for reactivation faults in coal mines .
Computations using the Coulomb v3.3 software package on changes in stress and deformation have suggested that the San Andres fault could cause earthquakes, with a frictional coefficient of 0.4–0.45, corresponding to a friction angle of 21.8–24.2° . This is a lower threshold for the sliding friction angle of faults for normal earthquakes. If the value is abnormally higher or lower than the range of frictional slip fault angle, there are two possibilities:
Fluid injection can be used to increase pore pressure or create hydraulic stimulation ;
The fault anomalies are filled with clay minerals or coal, or have been mineralized.
Therefore, the prediction of the lower limit of the friction angle of major fault systems in the Bach Ho oilfield indicated that slippage can occur. Hence, the consideration of appropriate volumes of water injection and drilling design are important parameters for preventing slip. For estimations, a lower shear friction angle threshold of 23°, or a frictional coefficient of 0.42 can be used for this area.
Abnormal faults with small friction angles can cause slip in high pore pressure conditions. Clay that fills in a fault plane can also reduce the frictional resistance, leading to slip with very small friction angles of 1–2° [6, 40, 41]. The data in Tables 2 and 3 indicate that some conjugate faults with very low shear friction angles can cause slip because of clay mineral infill along the faults.
Barton and Choubey  conducted experiments to determine the static friction coefficients of different rocks, showing that sedimentary rocks slide on magmatic rocks (or vice versa) due to the residual friction angle of the fissure or natural fracture at an angle of 28–32°. These values lowered the upper threshold for the range of the friction angle to 23–31°.
The modern stress field in the coastal area of Tuy Hoa–Vung Tau was assumed to have an average coefficient of friction, i.e., μ = 0.5 (friction angle of 26.6°) at a depth of 5 km, which almost coincides with the offshore value obtained in this study. Calculation results show high reliability corresponding to the vertical, maximum horizontal, and minimum stresses of 127, 158, and 90.5 MPa, respectively. The estimated pore pressure was 49.1 MPa .
Multi-stage triaxial compression tests for the confinement conditions at 20 °C in natural granite fragments show that the average friction angle at normal temperatures is 29.3° , which has a deviation of 0.65° from the onshore friction angle estimated in this study. At 80 °C, the friction angle was reduced to 27.87°, which is similar to the friction angle of 27° estimated for the offshore region at a depth of 4000 m.
Dung et al.  used a critical frictional coefficient of μ = 0.6 to determine the direction of open fractures in the granite basement of the Cuu Long Basin. The results showed two main directions of the critical stress fractures: (1) north–south, with a dip direction toward the east or west and a dip angle of ~75°, and (2) northwest–southeast, with a dip direction toward the northeast or southwest and a dip angle of ~75°. The frictional coefficient (μ = 0.6, corresponding to a friction angle of 31°) of the basement rock coincides with the upper threshold of the average value estimated in this study.
The estimation of statistical probability is an essential tool used in modern data analysis. The main goal of this study was to statistically process a friction angle dataset. The ranges of acceptable threshold values (minimum, mean, and maximum) at a 95% confidence level confirm the normal distributions of data from each of the studied regions. The estimates can be summarized as follows:
Onshore region: The threshold value ranged from 25.8° to 31.5°, with an average of 28.65°, according to a frictional coefficient of 0.55.
Offshore region: The threshold value ranged from 23° to 31°, with an average of 27°, according to a frictional coefficient of 0.51.
The recommended initial friction angle is the lowest median value (19.12°) obtained for the offshore area, which is higher than the 16.7° angle obtained when using the method of . The estimations of the friction coefficients and friction angles using statistical analyses are highly reliable and can contribute significantly to other underground resource exploitation projects in the study area.
The study was limited by any measurement error associated with the fracture morphology and from the sliding angles of kinematic faults. Errors and confusion regarding the identification of the conjugate angles in outcrops may have also occurred.
Future research should focus on building and classifying kinematic fractures in rock types with different ages. An improved identification method for matching pairs of conjugate fractures in outcrops is also required. Furthermore, we must develop an isometric map of the distribution of frictional coefficients, friction angles, and fracture strengths to further aid in the design of underground wells and resource exploitation.
This research is funded by the Office of National Science and Technology Research Programs, the Ministry of Science and Technology under grant number KC.09.38/16-20. We appreciate Vietnam–Russia Joint Venture, Hyosung Vina Petrochemical Co., Ltd, and Hoan Long – Hoan Vu JOC for providing the valuable data. In addition, this research is funded by the Ministry of Science and Technology of Vietnam, under grant number NĐT.48.KR/18. We would like to thank General Department of Geology and Minerals of Vietnam, Ministry of Natural Resource and Environment, Ho Chi Minh City University of Technology (HCMUT), and University of Science VNU-HCM for allowing the use of facilities during this study.
- Wyllie D.C., Norrish N.I. (1996) Rock strength properties and their measurement, TRB Special Report 247, pp. 372–390. [Google Scholar]
- Wyllie D.C. (2014) Rock fall engineering, CRC Press. [Google Scholar]
- Byerlee J. (1978) Friction of rocks, in: Byerlee J.D., Wyss M. (eds), Rock friction and earthquake prediction, Birkhäuser, Basel, pp. 615–626. [Google Scholar]
- Barton N. (1973) Review of a new shear-strength criterion for rock joints, Eng. Geol. 7, 4, 287–332. [Google Scholar]
- Jaeger J.C., Cook N.G., Zimmerman R. (2009) Fundamentals of rock mechanics, John Wiley & Sons. [Google Scholar]
- Yukutake Y., Takeda T., Yoshida A. (2015) The applicability of frictional reactivation theory to active faults in Japan based on slip tendency analysis, Earth Planet Sci. Lett. 411, 188–198. [Google Scholar]
- Middleton T.A., Copley A. (2014) Constraining fault friction by re-examining earthquake nodal plane dips, Geophys. J. Int. 196, 2, 671–680. [Google Scholar]
- Alber M. (2013) Strength of faults – A concern for mining engineers? in: ISRM International Symposium-EUROCK 2013, International Society for Rock Mech Rock Eng. [Google Scholar]
- Tamura H., Kambayashi Y. (2016) Estimation of the coefficient of static friction of surface by analyzing photo images, in: Czarnowski I., Caballero A., Howlett R., Jain L. (eds), Intelligent Decision Technologies 2016, Springer, Cham, pp. 15–26. [Google Scholar]
- Hencher S.R. (1989) Laboratory direct shear testing of rock discontinuities, Ground Eng. 22, 2, 24–31. [Google Scholar]
- Hencher S.R. (1995) Interpretation of direct shear tests on rock joints, in: Proceedings 35th US Symposium on Rock Mechanics, Lake Tahoe, pp. 99–106. [Google Scholar]
- Hencher S.R. (2012) Discussion of Alejano, Gonzalez and Muralha (2012), Rock. Mech. Rock. Eng. 45, 6, 1137–1139. [Google Scholar]
- Ghani R.A., Goh T.L., Hariri A.M., Baizura Y.N. (2011) Field and laboratory-based approach for the determination of friction angle of geological discontinuities of Malaysian granites, ASEAN J. Sci. Technol. Dev. 28, 2, 151–155. [Google Scholar]
- Boulton C., Barth N.C., Moore D.E., Lockner D.A., Townend J., Faulkner D.R. (2018) Frictional properties and 3-D stress analysis of the southern Alpine Fault, New Zealand, J. Struct. Geol. 114, 43–54. [Google Scholar]
- Jeanne P., Rutqvist J., Foxall W., Rinaldi A.P., Wainwright H.M., Zhou Q., Birkholzer J., Layland-Bachmann C. (2017) Effects of the distribution and evolution of the coefficient of friction along a fault on the assessment of the seismic activity associated with a hypothetical industrial-scale geologic CO2 sequestration operation, Int. J. Greenh. Gas Control 66, 254–263. [Google Scholar]
- Anderson E.M. (1951) The dynamics of faulting and dyke formation with applications to Britain, Oliver and Boyd. [Google Scholar]
- Binh N.T.T., Tokunaga T., Okui A. (2004) In-situ stress and pore pressure fields in the North Cuu Long Basin, offshore Vietnam, in: SPE Asia Pacific Conference on Integrated Modelling for Asset Management, Society of Petroleum Engineers. [Google Scholar]
- Binh N.T.T., Tokunaga T., Goulty N.R., Son H.P., Van Binh M. (2011) Stress state in the Cuu Long and Nam Con Son basins, offshore Vietnam, Mar. Pet. Geol. 28, 5, 973–979. [Google Scholar]
- Gan Q., Elsworth D. (2014) Analysis of fluid injection-induced fault reactivation and seismic slip in geothermal reservoirs, J. Geophys. Res. Solid Earth 119, 4, 3340–3353. [Google Scholar]
- Lama R.D., Vutukuri V.S. (1978) Handbook on mechanical properties of rocks-testing techniques and results-volume iii, Vol. 3, No. 2, Trans Tech Publications. [Google Scholar]
- Son H.V., Khanh P.D., Dung P.T. (2015) Analyzing the effect of the stress field on the wellbore stability, Bachelor thesis, Ho Chi Minh University of Technology, Vietnam. [Google Scholar]
- Sato K. (2016) A computerized method to estimate the friction coefficient from the orientation distribution of Meso-scale faults, J. Struct. Geol. 89, 44–53. [Google Scholar]
- Gephart J.W., Forsyth D.W. (1984) An improved method for determining the regional stress tensor using earthquake focal mechanism data: application to the San Fernando earthquake sequence, J. Geophys. Res. Solid Earth 89, B11, 9305–9320. [Google Scholar]
- Michael A.J., Geller R.J. (1984) Linear moment tensor inversion for shallow thrust earthquakes combining first-motion and surface wave data, J. Geophys. Res. Solid Earth 89, B3, 1889–1897. [Google Scholar]
- Holloway N.H. (1982) North Palawan block, Philippines – Its relation to Asian mainland and role in evolution of South China Sea, Am. Assoc. Pet. Geol. Bull. 66, 9, 1355–1383. [Google Scholar]
- Ru K., Pigott J.D. (1986) Episodic rifting and subsidence in the South China Sea, Am. Assoc. Pet. Geol. Bull. 70, 9, 1136–1155. [Google Scholar]
- Taylor B., Hayes D.E. (1983) Origin and history of the South China Sea basin, GMS 27, 23–56. [Google Scholar]
- Angelier J. (1994) Fault slip analysis and Paleostress reconstruction, in: Continental deformation, P.L. Hancock (ed), Pergamon Press, Oxford, pp. 53–100. [Google Scholar]
- Delvaux D., Sperner B. (2003) New aspects of tectonic stress inversion with reference to the TENSOR program, Geol. Soc. Spec. Publ., London 212, 1, 75–100. [Google Scholar]
- Allmendinger R.W., Cardozo N., Fisher D.M. (2012) Structural geology algorithms: Vectors and tensors, Cambridge University Press. [Google Scholar]
- Angelier J. (1989) From orientation to magnitudes in paleostress determinations using fault slip data, J. Struct. Geol. 11, 1–2, 37–50. [Google Scholar]
- Hyosung Co. Ltd (2018) LPG Underground Storage Project in Mined Rock Cavern in Cai Mep – Geological Investigation Report, Hyosung Co. Ltd. [Google Scholar]
- Thuc P.D. (2018) Researching and proposing engineering and technical solutions to improve oil recovery coefficient at the last stage of fractured basement reservoir in White Tiger field. Part II: Evaluation of the effectiveness of maintenance reservoir pressure, exploitation status of each region, existence and causes, Petroleum Magazine, Vietsovpetro. [Google Scholar]
- Wyllie D.C., Mah C.W. (2004) Rock slope engineering: Civil and Mining, CRC Press. [Google Scholar]
- Allmendinger R.W. (2018) Stereonet 10, accessed 1, 11, 2018. [Google Scholar]
- Gzovsky M.V. (1975) Basics of tectonophysics, Science Moscow, p. 535. [Google Scholar]
- Sibson R.H. (1985) A note on fault reactivation, J. Struct. Geol. 7, 6, 751–754. [Google Scholar]
- King G.C., Stein R.S., Lin J. (1994) Static stress changes and the triggering of earthquakes, Bull. Seismol. Soc. Am. 84, 3, 935–953. [Google Scholar]
- Toda S., Stein R.S., Sevilgen V., Lin J. (2011) Coulomb 3.3 Graphic-rich deformation and stress-change software for earthquake, tectonic, and volcano research and teaching – user guide, US Geological Survey Open-File Report 1060, p. 63. [Google Scholar]
- Yin Z.M., Ranalli G. (1995) Estimation of the frictional strength of faults from inversion of fault-slip data: a new method, J. Struct. Geol. 17, 9, 1327–1335. [Google Scholar]
- Sibson R.H. (2000) Fluid involvement in normal faulting, J. Geodyn. 29, 3–5, 469–499. [Google Scholar]
- Barton N.C., Choubey V.V. (1976) The shear strength of rock joints in theory and practice, Rock. Mech. Rock. Eng. 10, 1–54. [Google Scholar]
- Loc N.L., Linh D.V., Van D.Q., Sang N.T., Tan V.T., Ngoc P.Q., Dien T.N. (2020) Study on recent stress change in the Tuy Hoa – Vung Tau marine region arising from oil and gas exploitation activities, J. Geol. 117–133. [Google Scholar]
- Kim T., Zhang N., Jeon S. (2013) A study on shear characteristics of a smooth rock surface under different thermal, hydro and mechanical conditions, in: ISRM International Symposium-EUROCK 2013, International Society for Rock Mech Rock Eng. [Google Scholar]
- Dung T.Q., Dung P.T., Ha L.T., Thuan N.V. (2020) Applying geomechanical model and critical stress to identify open fractures in basement reservoir, Petrovietnam J. 3, 30–37. [Google Scholar]
Mohr–Coulomb failure criterion in the effective stress domain. (a) Shear plane under triaxial compression. (b) Mohr circles and an attenuation contour cut linearly (tangential line of the Mohr circles), where is the angle of internal friction.
|In the text|