An Experimental Investigation of Shale Mechanical Properties Through Drained and Undrained Test Mechanisms
Md. Aminul Islam and Paal Skalle
Shale mechanical properties are evaluated from laboratory tests after a complex workflow that covers tasks from sampling to testing. Due to the heterogeneous nature of shale, it is common to obtain inconsistent test results when evaluating the mechanical properties. In practice, this variation creates errors in numerical modeling when test results differ significantly, even when samples are from a similar core specimen. This is because the fundamental models are based on the supplied test data and a gap is, therefore, always observed during calibration. Thus, the overall goal of this study was to provide additional insight regarding the organization of the non-linear model input parameters in borehole simulations and to assist other researchers involved in the rock physics-related research fields. To achieve this goal, the following parallel activities were carried out: (1) perform triaxial testing with different sample orientations, i.e., 0°, 45°, 60°, and 90°, including the Brazilian test and CT scans, to obtain a reasonably accurate description of the anisotropic properties of shale; (2) apply an accurate interpretative method to evaluate the elastic moduli of shale; (3) evaluate and quantify the mechanical properties of shale by accounting for the bed- dings plane, variable confinement pressures, drained and undrained test mechanisms, and cyclic versus monotonic test effects. The experimental results indicate that shale has a significant level of heterogeneity. Post failure analysis confirmed that the failure plane coincides nicely with the weak bedding plane. The drained Poisson’s ratios were, on average, 40% or lower than the undrained rates. The drained Young’s modulus was approximately 48% that of the undrained value.
These mechanical properties were significantly impacted by the bedding plane orientation. Based on the Brazilian test, the predicted tensile strength perpendicular to the bedding plane was 12% lower than the value obtained using the standard isotropic correlation test. The cyclic tests provided approximately 6 % higher rock strength than those predicted by the monotonic tests.
Most borehole stability problems occur during drilling in overburden shale or shale-like materials, such as mudstones. Borehole instabilities are typically associated with high pore pressures (PPs) in shale located immediately above the hydrocarbon reservoir. Worldwide, there is an increasing tendency to go toward deeper high-pressure/high-temperature (HPHT) reservoirs. Consequently, the drilling window margin (collapse–fracture) is reduced, and more emphasis is placed on borehole stability predictions to avoid borehole collapse or fracture. Moreover, field observations imply that borehole collapse is strongly dependent on the wellbore orientation, especially in areas where the formation is strongly bedded and where the stress field is non-isotropic. Under such circumstances, accurate information on the rock strength and rock failure behavior in overburden shale is crucial to improve drilling safety. Knowledge of the mechanical properties of shale is vital to implement any 3D shale anisotropy borehole instability model (Crook et al. 2002; Søreide et al. 2009).
As an example, to implement the modified Cam Clay model for the shale-related borehole instability assessment, mechanical properties and matrix anisotropic parameters, which consist of 21 attributes in total, are required (Søreide et al. 2009). Most of the parameters are related to the mechanical characterization of shale. However, a wide range in their magnitude was observed specifically for the Young’s modulus, Poisson’s ratio, friction angle, and the shear and bulk moduli. An accurate interpretation, together with sound engineering judgment, appear necessary to evaluate and quantify the mechanical properties of shale.
This paper has addressed these concerns by interpreting the results obtained from extensive laboratory testing of shale. It has previously been assumed that shale is not a reservoir rock and is, thus, not interesting in terms of hydrocarbon production (Hofmann and Johnson 2006). In general, the reservoir sections are the primary target, so shale samples from the overburden sections are limited. Some previous publications on the mechanical properties of tested shale were reported by Closmann and Bradley (1979), Steiger and Leung (1992), and Hoang and Abousleiman (2010). However, the scarcity of deep shale testing samples has lead to a lack of published shale elasticity data.
Meanwhile, the interest in shale has increased over the past few years, which was understandably triggered by recognition by the oil industry that the primary migration of oil is incomplete in some shale. Potential challenges to the understanding of shale are shale geochemistry, anisotropy, petrophysical properties, and the response of the seal to pressure as a barrier. The growing interest in shale-gas and caprock integrity for CO2 storage has resulted in a rising demand for fundamental rock property data. Performing work on stimulation wells to create a complex fracture network in a shale-gas field is a difficult task, particularly the drilling of horizontal wells. However, this is the only viable method at this stage to develop a drainage strategy in the shale-gas reservoirs. Fluid flow modeling to develop an analytical solution in this particular event is one current issue where accurate shale property knowledge is essential. Moreover, infill drilling is a practice used by the industry to penetrate interbedded shale-sand formations, requiring the most accurate rock mechanical models to correctly assess the mud weight operational window. Accurate study of all issues requires a database of shale characteristics and its mechanical properties.
The prime consideration when evaluating and modeling borehole stability problems in overburden shale or in an interbedded reservoir is the lack of relevant test data to describe shale properties accurately. In reality, cutting a real shale specimen from an interbedded sand-shale layer is a challenging, costly, and time-consuming task. Experimental data for real shale specimens at the maintained downhole conditions are rare, while data on weak shales (i.e., Pierre-1) are fairly abundant. It is, therefore, a well-accepted approach to execute experimental investigations on outcrop shale and use the data to provide the necessary material data sets for the fundamental borehole stability model. In practice, fundamental models are calibrated against field cases and readjusted. Many investigators (Johnston 1987; Schmitt et al. 1994; Horsrud et al. 1994; Horsrud 1998; Spaar et al. 1995; Holt et al. 1996; Claesson and Bohloli 2002; Crook et al. 2002; Søreide et al. 2009; Fjær et al. 2008) have, instead, used outcrop shales to calibrate borehole stability models. The present work deals with data sets for Pierre-1 shale at room temperature and at the high confining pressure present at downhole conditions. It was our intention to compile published shale outcrops test data. This was impossible due to the inconsistency in the reported information and evaluation techniques. It was also unclear as to which types of outcrop batch were used.
Shale is a rock with laminated structure, and the bedding plane orientation determines the mechanical properties of each shale sample. It is a well-accepted statement and a well-defined concept in soil and rock mechanics textbooks that the shale bedding plane leads to borehole instability. Jaeger and Cook (1979) reported the laboratory test results of mechanical properties of layered rocks during failure. A comprehensive study, together with field evidence on the plane of weakness and anisotropy, were introduced and reported in the oil industry by Bradley (1979), Maury and Santarelli (1988), and Aadnøy and Chenevert (1987). Aadnøy and Chenevert (1987) investigated the effects of wellbore inclination, anisotropic stresses, and anisotropic rock strength by modeling highly inclined boreholes. Under certain conditions, the rock would fail along the planes of weakness. Because of the geomechanical properties of the weak plane (alignment of phyllosilicates due to overburden digenesis), slip surfaces may exhibit a significantly higher potential to fail compared with stronger rock units, such as limestone and sandstone. However, this statement is valid only for particular geological formations and may not be applicable in other geological arenas. For example, the same author performed a separate study (Aadnøy et al. 2009) on a tightly folded structure in British Columbia, and it was observed there that the planes of weakness in bedded rocks may lead to borehole collapse. However, in a three- dimensional space, there were combinations along wellbore inclinations and azimuths, instead of along the weak planes, where the sample failed. The bedding plane is, thus, a key issue in this study. In general, the critical parameters that cause instability are the planes of weakness, anisotropy, the relative normal stress values on the borehole, the well trajectory, and the relative angle between the borehole and the bedding plane (Bradley 1979; Aadnøy and Chenevert 1987; Aadnøy et al. 2009; Crook et al. 2002; Fjær et al. 2008; Søreide et al. 2009; Islam et al. 2009).
As discussed earlier, many researchers have studied the properties of clay and shale. Selected test results of mechanical shale properties under drained conditions are published and reported by Closmann and Bradley (1979), Horsrud et al. (1998), Crook et al. (2002), Søreide et al. (2009), and Hoang and Abousleiman (2010). These previous investigations can also be useful for this study to some degree. However, determination of the mechanical properties of shale under drained and undrained test conditions with variable confinement pressure is currently insufficient. To overcome this limitation, this study performed an extensive test program to explore shale heterogeneity and its directional properties. The following parallel activities were performed:
- Characterization and evaluation of shale heterogeneity through a drained and undrained stress path. The mechanical properties and post failure material behavior were determined under different sample orientations, confinement pressures, and loading steps.
- Determination of the estimated and quantified mechanical properties of shale.
- Proposition of proper data sets for a borehole stability model.
Both compressibility and triaxial tests were utilized, with emphasis on evaluating the elasticity, yielding, and failure response as a function of the confining pressure and bedding plane orientation. The petrophysical testing included porosity, permeability, and bulk-and solid-density measurements. A total of three triaxial tests with different sample orientations (i.e., 0°, 45°, 60°, and 90°), including Brazilian tests and CT scans, were performed to obtain a reasonably accurate description of shale. Moreover, the cyclic stress-dependent rock stiffness and strength were evaluated and compared with the monotonic test results. The method used to achieve relatively consistent test results is briefly discussed at a later stage in this paper.
2 Materials and Methods
2.1 Material and Samples
The outcrop Pierre-1 shale was used in this study for the following reasons:
- The scarcity of real shale specimens.
- The difficulty in obtaining intact samples of real shale specimens due to splitting along the bedding planes during the coring process.
- Mechanically, Pierre-1 shale exhibits plastic and anisotropic behavior that could be used in borehole simulation modeling of real shale.
The outcrop shale block was preserved by wax when received. A core barrel was used to drill the core samples from a large block stored in laboratory-grade oil (MarcolTM) to prevent fluid losses. Both the drilling and the subsequent end-grinding of the plugs to obtain sufficiently plane and parallel end surfaces were performed whilst surrounded in the same oil. The prepared samples were subsequently stored in MarcolTM until they were assembled in the test cell to avoid desiccation effects.
The samples were equilibrated to the brine (3.5 w% NaCl) prior to ramping up stresses to consolidation level, so as the Pierre is a outcrop material with initial low stress levels, this should assume saturation close to 1 during the tests. Porosity was determined from fluid loss by heating typically to about 105°C.
The test samples were 38mm in diameter and had a length of approximately two times the diameter. A summary of the test specimens and other necessary information are presented in Table 1. The mineralogical and petrophysical properties of the Pierre-1 shale used are presented in Tables 2, 3, and 4.
We used samples from two neighboring blocks. All triaxial samples had densities between 2.32 and 2.34 g/cm3, with an average density of 2.33 ± 0.005 g/cm3. We also completed separate density and porosity measurements on the two blocks from which we drilled the samples. From the loss of pore fluid measurement, the densities and porosities were determined to be 2.33g/cm3 and 23.2 % (block 2) and 2.34 g/cm3 and 22.4 % (block 3), respectively, i.e., they were fairly similar.
2.2 Sample Selection
CT scans were performed to check for possible heterogeneities in the test samples. To illustrate the heterogeneity effect of Pierre-1 shale, two examples of CT scans (radial slices) are presented in Fig. 1a, b. Lighter areas generally correspond to denser (higher absorption) material. The very light color on one of the scans is likely calcite (Fig. 1a). Such samples were discarded from the consolidated isotropic undrained triaxial compression tests (CIU). The other image contains a less light, gray region/layer (Fig. 1b), which is not calcite but likely represents a more silty region. Such samples were used in the consolidated isotropic drained triaxial compression tests (CID).
Fig. 1 CIU and CID tests for Pierre-1 shale; a,b heterogeneous nature of shale. The very light color on one of the scans is probably calcite (a). A less light, gray region/layer (Fig. 1b) is not calcite but probably a more silty region.
2.3 Bedding Plane and Loading
The bedding plane orientation in the triaxial test is denoted by the angle h, which is measured clockwise from the loading direction relative to the bedding plane (Fig. 2). The monotonic CIU tests were performed with angles of 0°,45°, 60°, and 90° between the load axis and the bedding plane. The same is true for the CID tests, with the angles of 0° and 90°. The postfailure samples are presented in Fig. 3. Through a postmortem analysis of the tests samples, it is observed that, in most cases, the failure plane transverses nicely through the weak bedding plane.
Fig. 2 Principle drawing including bedding angle and loading; a drilled normal to bedding (h = 0°), b drilled parallel to bedding (h = 90°), and c drilled inclined to normal to bedding (h = inclined). The angle h is measured clockwise from the loading direction relative to the bedding plane.
Fig. 3 Post failure views of the Pierre-1 shale samples performed through monotonic loading; a loading normal to bedding (θ=0°), b loading 45° to the normal to bedding (θ=45°), c loading 60° to the normal to bedding (θ= 60°), and d loading parallel to bedding (θ= 90°).
2.4 Brazilian Test
The tensile strength of rocks is among the most important parameters influencing rock deformability, rock crushing, and blasting results. To calculate the tensile strength from the indirect tensile (Brazilian) test, one must know the principal tensile stress, in particular, at the disc center, where a crack initiates. This stress can be calculated by an analytical solution given by Claesson and Bohloli (2002).
A Brazilian test was conducted to estimate the indirect tensile strength at both θ= 90° and 0° conditions. One sample (length=19.55 mm and diameter=37.70 mm) with a bulk density of 2.32g/cm3 showed a rock tensile strength of 0.62 and 0.76 MPa at θ= 90° and 0°, respectively. This measured tensile strength through θ= 0° is 12% lower than the result of the standard isotropic correlation given by Claesson and Bohloli (2002). The test results also indicated that the tensile strength is higher for vertical than for horizontal bedding.
2.5 Triaxial Testing
A schematic drawing is presented in Fig. 4 to illustrate the methods of collecting data using the triaxial test.
Fig. 4 Principle sketch of the interior of a triaxial cell, showing the loading piston with fluid ports, optional sintered plates for fluid distribution, a radial strain jig (suspended from the specimen) measuring two orthogonal diameters, and an axial strain jig measuring the change in sample length between the pistons. It is also possible to measure axial strain directly on the sample, to measure radial strain in more than two directions, or measure the change in circumference by a chain around the sample (after Fjær et al. 2008).
2.6 Drained and Undrained Test Mechanisms
During drilling in low permeable shale, a typical undrained stress relaxation response is observed at the borehole wall. With extended time, excess PP will dissipate and equilibrate with the surrounding pressure. A redistribution of stresses might provide a time-delayed failure in the shale formation. In this study, fluid permeability was measured with the pulse-decay technique at a constant confining pressure (15 MPa). This is a standard pulse-decay method on samples subjected to hydrostatic stress. We used brine as the contacting pore fluid. The measurements were performed carefully to avoid effects from temperature variations and other factors influencing the results.
Due to the low shale permeability, steady-state type measurement techniques are not generally applied. Instead, techniques based on transient measurements are used (Horsrud et al. 1998). We used a technique where a thin (6- mm) circular disk (38 mm in diameter) was placed into a pressure cell, where PP and hydrostatic confining pressure could be applied separately. After loading to predetermined levels of PP and confining pressure, ample time was given for the sample to consolidate. The PP was increased by 0.5 MPa on one side of the sample and reduced by the same amount on the other side (see Fig. 5). The decay in this pressure difference was then recorded by a differential pressure transducer, and the permeability was calculated from this time response.
Fig. 5 Illustration of the pressure transient method used for permeability estimation (Horsrud et al. 1994).
From the permeability (k) measurements within a confining pressure of 15 MPa and a PP of 5 MPa, k was 14 and 49 ηD at θ=0° and 90°, respectively. Because shale is expected to be anisotropic, it is reasonable to expect the permeability to be anisotropic and to expect the largest values parallel to the bedding. Islam et al. (2009) reported that the magnitude of the consolidation time for the finite diffusive length adjacent to the borehole wall was only in the range of hours and days for samples with a permeability of 100 ηD. Therefore, drained tests were necessary for time-delayed material failure analysis. The drained stiffness was somewhat lower than the undrained stiffness, as drained stiffness accounts for effective stress.
Generally, in drained tests, the outlets through the pistons will be open such that the pore fluid pressure can be maintained at any prescribed value. During the hold periods, the PP equilibration is expected to take place inside the sample. During the testing of shale, the drained condition means that the test must run sufficiently slowly to avoid unacceptable pressure build-up when the sample deforms. In this case, a strain rate of 2x10-8s-1 was maintained until failure. How long the drained test will take is fairly difficult to estimate and basically depends on the strain rate used during the triaxial loading of the specific sample. The standard procedure is to estimate this strain rate from the consolidation behavior in the initial hydrostatic phase, which depends on the specific sample, confining stress, and PP, as well as how effective the PP drainage is working for a given sample. The desired strain rate for the drained condition is often 10 times slower than that for the undrained condition. For example, based on an assumed strain rate of 2x10-7s-1 (and a 10-mm axial strain at failure) and 24 h initial consolidation, 68 h was required to complete an undrained test (sample drilled parallel to the bedding; θ= 90°), whereas for the drained test, it would be approximately 10 times longer. A brief discussion is included on this topic in the next section. From a theoretical point of view, drained conditions mean that the PP in the sample is an independent variable, while for undrained tests, it is a dependent variable. In real-life borehole stability simulations, both drained and undrained rock mechanical properties are necessary because they depend on the time scale. In the case of a short borehole simulation, undrained rock properties can be used, whereas for long simulations, the drained rock properties are applied.
2.7 Triaxial Testing of Shales
The triaxial tests were run as consolidated drained and undrained tests in a servo-controlled load frame. Consolidated undrained (CIU) tests were selected for the testing of shales (Steiger and Leung 1992; Nakken et al. 1989; Horsrud et al. 1998). The procedure was adopted from soil testing (e.g., Head 1984). This test consisted of three distinct phases:
- Loading to the predetermined level of confining pressure and PP.
- Consolidating, maintaining a constant confining pressure, and allowing drainage of the pore fluid against a constant PP.
- Undrained axial loading under a constant axial displacement rate beyond failure of the sample.
The internal instrumentation of the test sample is shown in Fig. 4. In addition to the measurements of the external load, pressure, and deformations, the PP on both ends of the sample were recorded, and acoustic wave trends in both the axial and radial directions were also recorded.
CIU tests were performed at angles of 0°, 45°, 60°, and 90° between the sample axis and the bedding plane. In addition, two CID tests were performed to obtain the effective stress parameters. For triaxial testing, all of the samples were consolidated up to 10 MPa of PP. The confining pressure was ramped up to 20, 25, and 30 MPa to simulate a closer range of downhole effective stress conditions.
The samples were left for 24h to establish consolidation. We monitored the strain (deformation) so as to evaluate when consolidation would be sufficient. In general, 24 h of consolidation for our sample geometries is considered to be a sufficient amount of time. The time dependency of the consolidation phase was also used to determine the loading rate during the main triaxial phases (which were run in a constant strain mode). In some undrained tests, the strain rate was 1–2x10-7s-1. For the drained tests, the rate was approximately 2x10-8s-1. In addition, the size of the test sample and the time required to perform a test depended highly on the permeability of the shale (Nes et al. 1998; Islam et al. 2009; Loret et al. 2001) and the fluid viscosity.
A typical response from the triaxial loading of the consolidation type is shown in Fig. 6. The drained test was extremely time consuming up to sample failure (see Fig. 6e). For a drained test, approximately 20 days was required to reach yielding with a strain rate of 2x10-8 s-1. This is the main reason for running the triaxial loading part in the undrained mode. Concepts from soil mechanical testing can be applied to assist in determining when consolidation is complete and also when the appropriate displacement rate is achieved in the undrained part of the test (Katsube et al. 1991). Determination of the strain rate in the undrained part was also based on the consolidation response of the sample to ensure that PP equilibrium was achieved throughout the sample. To reduce the test time as much as possible, both axial and radial drainage of the sample were included.
Fig. 6 Triaxial testing of Pierre-1 shale subjected to consolidation, triaxial loading, and development of PP; a the timeline shows the development of PP and failure load with variable confinement pressure, at sample θ=0°, b effect of bedding plane on PP and failure load, c variable confining pressure to estimate failure load and to develop PP against total test time at θ=90°, d material failure and PP response under constant confinement pressure with different bedding plane orientations, e failure response against total test time for CID for two samples, θ=0° and 90°, f PP and failure response against axial deformation subjected to variable confinement pressure. The two dotted circles indicate abnormal strain rate and unusual PP development.
In the case of CIU tests conditions, a total test time of 2–4 days, on average, was required, including consolidation of approximately 20–25 h to failure. The bedding plane sample at 60° took a maximum of 5 days to fail (Fig. 6d). From the material stiffness analysis, the strongest sample was determined to be the sample drilled parallel to the bedding and would normally take a longer time to fail.
The PP increased up to a certain level with increasing axial stress, after which the PP tended to decrease until the sample failed. The response of PP during triaxial testing is presented in Fig. 6. Referring to Fig. 6a, with 30 MPa of confining pressure, the net PP increased 9 MPa, but at the failure position, the PP dropped 7 MPa and decreased further. This decrease indicates that the sample dilates before failure. In general, induced cracks or microfractures were created before the material reached failure. Basically, induced cracks or microfractures increase the pore volume (Fjær et al. 2008). The increased pore volume led to material dilatancy, consequently reducing the PP. In general, the failure position is controlled by increased pore volume, which reduces the PP. The development of PP under different bedding planes and confining pressures is presented in Fig. 6a–f. From this analysis, PP changed due to the change of the bedding plane. The confining pressure was quantified. The PP development versus confining pressure response is shown in Fig. 6f. This analysis indicates that the PP development reached a maximum at 0° bedding and was the lowest for 90° bedding. For the CID condition, PP was constant (see Fig. 6e). Thus, the actual effective stress-related parameters can be obtained through such tests.
PP development is linked with the effective stress to define material failures. According to the Mohr–Coulomb failure envelope, a material with 0° or 45° bedding to normal will reach the failure line prior to other orientations. One anomaly is indicated by a dotted circle at 90° bedding in Fig. 6e. The anomaly may be due to a contribution in the pore-pressure fluid volume or uncertainities related to existing cracks. The stiffness analysis trends are presented in Fig. 6b, which shows that samples drilled parallel to the bedding are the strongest. A summary of the rock failure related data applicable at different sample orientations is presented in Table 1. The summary implies that net PP developments during CIU tests are not similar. The bedding plane and material heterogeneity are two key concerns to control the material failure and the development of PP.
3 Elastic Moduli of Shale—Theory
The elastic moduli of a material are primarily described by its Young’s modulus (E) and Poisson’s ratio (ν). These two parameters are interpreted in this study using the triaxial test results through a stress–strain trend analysis. Several simplifications have been introduced to calculate the elastic properties of shale from the triaxial tests. One key simplification is that the elastic stress–strain response is assumed to be linear. This is generally not the case for both the hydrostatic and triaxial phases. Thus, when presenting the elastic moduli, it is essential to define how the interpretation has been made. The following alternatives are applicable to estimate the Young’s modulus (Wood 1990; Fjær et al. 2008):
- Initial modulus, given as the initial slope of the stress–strain curve.
- Secant modulus, measured at a fixed percentage of the peak stress.
- Tangent modulus, given at a specific percentage of the peak stress.
- Average modulus, given within a specific maximum and minimum stress level.
Of these alternatives, the tangent modulus at 50% of the peak stress value was chosen in this study simply because it is the most widely used method. The Young’s modulus was determined as the tangent modulus measured at 50% of the peak stress (see Fig. 7a) and is, therefore, denoted as E50. The Poisson’s ratio was determined from the inclination of a straight line that passes from the origin to the point of the curve that corresponds to 50 % of the peak stress in an εvol versus ε1 diagram (see Fig. 7b), and is, therefore, denoted as ν50.
According to the American Society for Testing and Materials (ASTM) and International Society for Rock Mechanics (ISRM) standards, it is recommended to use the deformation rates at 50% of the peak stress level for the determination of ν. However, if the curve is strongly non-linear, complete information can only be given if the entire curve is presented. Shale is an anisotropic material, thus, the Poisson’s ratio and the Young’s modulus are not the proper parameters to describe the mechanical behavior. The elastic stiffness tensor, defined by Fjær et al. (2008), is required for an accurate description of the shale.
Fig. 7 Elastic moduli estimation techniques of shale; a Young’s modulus, b Poisson’s ratio. The symbol Ψ is the dilatancy angle and ν is the Poisson’s ratio. In practice, the ductile region may be very small.
4 Elastic Moduli—Results and Discussion
4.1 Strain Measurements
We have carefully leveled the strike and slip directions of the test samples. One arm of the cantilever was placed in the strike direction (along the bedding), and another arm was placed in the slip direction (perpendicular to the bedding). Radial deformation was measured by two pairs of calipers (pairs #1 and #2) mounted orthogonal to one another. One of the pairs always measured deformation along the strike direction.
Because we performed several tests along different bedding inclinations, a summary of sample orientations and caliper positions to measure the radial deformations are presented below:
- Perpendicular to bedding (θ= 0°): pair #1 measures deformation along the strike direction.
- Angle to bed (θ=45°): pair #1 measures deformation along the dip direction.
- Angle to bed (θ=60°): pair #1 measures deformation along the dip direction.
- Parallel to bedding (θ=90°): pair #1 measures deformation along the slip direction.
4.2 Elastic Moduli of Shale—the CIU Case
During the CIU testing of shale, the confining stress was generally kept constant on a total stress basis. However, the effective confining stress is affected by the PP response. Consequently, the anisotropic effective stress parameters were difficult to obtain because the radial effective stress changed. In the CIU test conditions, the stiffness of the water was also taken into account and, because there was no volumetric change (because water is incompressible), the Young’s modulus was expected to be greater than for the CID case. In general, for partially undrained conditions, the stiffness is higher and the volumetric strain is reduced (Holt et al. 2011; Rozhko 2011).
As mentioned earlier, the Young’s modulus normally appeared to be greater for the undrained conditions. In addition, the higher the confining stresses, the higher the Young’s modulus for the same type of test and the same angle to the bedding. For similar confining pressures, the Young’s modulus was greater (the material is stiffer) for loading parallel to the bedding and smaller for loading perpendicular to the bedding. The stress–strain curves are presented in Fig. 8 to show how the material stiffness varied due to the effects of confining pressure and of the bedding plane.
The Young’s moduli and the Poisson’s ratios were calculated from Fig. 8, in both drained and undrained test conditions. The determined elastic parameters, including both E50 and ν50, together with the bulk modulus, shear modulus, and Lame modulus, are summarized in Table 5. The effect of the temperature on these elastic moduli is not included in this paper. Similar studies (Horsrud et al. 1998; Søreide et al. 2009) indicate that, at a downhole temperature of 120° C, a significant reduction in both stiffness (approx. -25%) and strength (approx. -35%) of the shale is observed. According to Søreide et al. (2009), the downhole undrained stiffness for North Sea shale at 140°C is 57 % of the stiffness at room temperature. One possible countermeasure of the temperature increase is to improve the saturation of the sample, weaken the rock frame, and/or increase the fluid modulus, which would reduce capillary effects. This can contribute significantly to the strength of rocks with small pores (Schmitt et al. 1994; Papamichos et al. 1997). In general, it was found that, compared with sandstone, the response of shales were more dependent on temperature and less dependent on pressure (Holt et al. 1996; Horsrud et al. 1998).
Fig.8 Stress–strain response subjected to different confining pressures and bedding angles; a θ=0°, b θ =90°, c θ =45°, and d θ=60°. At low confinement pressure, natural microfractures cause non-linear stress–strain curve response (marked with the dotted ellipse in a).
Figure 8a represents the stress–strain response at θ=0° bedding subjected to the effects of confinement. The material stiffness under CID and CIU tests differed significantly. Under CIU tests, the stiffness curves increased non-linearly up to a certain level and then decreased. A possible reason for this behavior is that the PP develops asymptotically up to a certain level close to the peak stress. The CIU test subjected to drilling normal to bedding would, thus, provide a higher development response of the PP. The developed PP under this test did not account for the effective stress because the volumetric strain was constant. However, for the CID test, the effective stresses play a role in determining the rock strength and stiffness. The initial non-linear trend in Fig. 8a (marked by a dotted elipse) may be due to the presence of induced cracks or natural microfractures. Through the postfailure analysis (Fig. 3) of these samples, it was inferred that the additional induced cracks or microfractures can further contribute to increased pore volume prior to failure. The pore fluid may penetrate through the microfracture network and further reduce the PP, resulting in reduced rock strength. Moreover, the natural properties of shale in terms of crack or microfracture orientation may also produce various nonlinear behaviors. Depending on the loading direction, the fracture network/crack can close or open. This variance leads to a decrease or increase in pore volume in response to the PP. A similar observation was found at 90 bedding (Fig. 8b). The concave-type CIU test results were mostly defined by the PP increase through the triaxial phase.
The volumetric-axial strain curve (Fig. 9) was used to calculate the Poisson’s ratio. The results are presented in Table 5. For the undrained tests, the value of the Poisson’s ratios ranged between 0.34 and 0.75, while for the drained case, it decreased to 0.2. This behavior is illustrated and discussed later in this paper.
4.3 Elastic Moduli of Shale—the CID Case
The drained test for low permeability rocks is essential to obtain the effective stress parameters. These effective stress parameters are important for borehole stability modeling, especially in shale. Directional stiffness parameters (Et and Ep) are essential for the shale anisotropy borehole stability model. The normal or transversal stiffness parameter (Et) was determined at sample θ=0° (load normal or transversal to the bedding), and Ep was estimated when θ=90° (loading parallel to the bedding). Søreide et al. (2009) showed analytically that the stiffness for the samples with loading parallel to the bedding (Ep) is higher than for loading perpendicular to the bedding (Et). In our case, the estimated values of Et and Ep were 0.65 and 1.55 GPa, respectively (see Fig. 10a). The estimated mean relative difference of the E between these two samples (in absolute value) was approximately 58 %. A significant difference between these two moduli indicated a strong heterogeneous nature.
Fig. 9 Poisson’s ratio from axial strain versus volumetric strain through CID tests; a h=0 , b h=90 . The dotted circle indicates the strain rate effect due to a temporary stop in the experiment.
Figure 10e presents the behavior of the Poisson’s ratios for samples at θ=0° and 90°. At θ=0°; the observed Poisson’s ratios were found to be similar. The two measured radial strains will always be equal because the sample plane acts as a ring and, thus, carries similar material behavior, which produces the same deformations in all directions in that particular plane. In the case of θ=90°, the material behavior will be anisotropic, which means that one radial deformation is parallel to the bedding and the other is perpendicular to the bedding. These two measured radial strains will never be equal, and the Poisson’s ratios will, therefore, be different. Figure 10e, f supports these conclusions.
The shear modulus is calculated based on the estimated Young’s moduli and Poisson’s ratios, and the resulting calculations are presented in Fig. 10c, d. The lower limit of the shear modulus for θ=0° was 0.24 GPa, but for θ=90°, it was 0.62 GPa (on average). In Fig. 10a, b, there is an irregularity, marked by dotted ellipses, caused by stopping the test at a high strain rate. This irregularity may affect the quality of the results. One point that should be carefully addressed is how to select the measured data points. The trend of the Young’s modulus at the horizontal bedding is almost linear, but it is non-linear at θ=90°.
Fig. 10 Shale elastic moduli at drained test conditions; a variation of E with axial strain, b variation of E with deviatoric stress, c, d variation of G with axial strain and deviatoric stress, and, finally, e, f variation of Poisson’s ratios with axial strain and deviatoric stress. The dotted ellipses indicate the abnormal strain rates due to a temporary stop in the experiment.
4.4 Elastic Moduli of Shale—the CIU Case with Constant Confining Pressure
At constant confining pressure (i.e., 25 MPa), the E and the ν were calculated and are presented in Fig. 11. The response of the E and m decreased exponentially with increasing axial strain or deviatoric stress. For example, in Fig. 11a, for 10 % axial strain, E for the 90°, 60°, 45°, and 0° samples are 2.1, 1.7, 1.5, and 1.4 GPa, respectively. The samples showed the highest stiffness at θ= 90° and the lowest stiffness at θ=0°. It is recommended to choose the E data point based on the expected axial strain or deviatoric stress condition.
Fig. 11 Estimation of Young’s modulus in CIU tests at different constant confinement pressures.
Two Poisson’s ratios were defined, ν1 and ν2. The values of ν1 vary significantly with respect to all bedding plane samples. These values are increasing exponentially but, in some cases, also increasing gradually versus changing deviatoric stress (Fig. 12a, b). For θ=0°, the magnitudes of ν1 and ν2 followed a non-linear trend and varied from 0.4 to 0.7. At other sample orientations (except θ=0°), ν1 showed extremely high values between 0.6 and 0.8. However, a significant reduction occurred in the interpreted ν2, with values varying between 0.25 and 0.4. At θ=0° and θ=90°, the Poisson’s ratios ν1=0.52 and 0.68, respectively, while ν2 was calculated to be 0.53 and 0.35, respectively.
Fig. 12 Estimation of Poisson’s ratios in CIU tests at different constant confinement pressures.
4.5 Effects of Variable Confining Pressure on E (for Both CIU and CID Tests)
The confining pressure influences material stiffness and shear strength. The higher the confining pressure, the stiffer the material and the higher the shear strength. The material stiffness was evaluated by changing the confining pressure, and the interpreted results are presented in Fig. 13. The variation of rock stiffness was clearly observed from this analysis.
The material stiffness for the CIU cases gradually decreased with increasing axial strain (Fig. 13a) and with increasing deviatoric stresses (Fig. 13b). Estimates of E for the CIU cases are confusing because it is difficult to find a linear trend, as in the CID cases. The anomaly of the CID curves (dotted circles) is due to a pause in the experiment at a high strain rate and reloading the sample to run the test again. This analysis helped us to choose E according to the projection of confining pressure and material deformation state.
Fig. 13 a–h Variation of the Young’s modulus subjected to variable confinement pressures.
4.6 Effects of Variable Confining Pressure on ν (for Both CIU and CID Tests)
For homogeneous materials, it is generally accepted that the higher the confining pressure, the higher the Poisson’s ratio. Figure 14a shows that, for samples at θ=0°, both ν1 and ν2 increased non-linearly with increasing confining pressure. The changing values are minor compared with other bedding orientations.
On the other hand, a large gap between ν1 and ν2 was found for the sample drilled at θ=90° (Fig. 14c, d), and similar trends were observed for θ=45° and 60°. At θ=90°, values of 0.3 and 0.7 were obtained for ν1 and ν2, respectively, and this difference is significant. In the undrained test, the Poisson’s ratios observed in the tested shale rocks tended to be larger than 0.5. This observation is supported by Aadnøy and Chenevert (1987), who showed that, for laminated materials, the Poisson’s ratio can be higher than 0.5 and is, in fact, dependent on the elasticity modulus ratio (a ratio between the vertical and horizontal moduli).
Fig. 14 a–h Quantifying the effects of confining pressure on the Poisson’s ratio.
In the part B of Experimental Investigation of Shale Mechanical Properties:
5 Cyclic Versus Monotonic Test Effects on Rock Strength
6 Rock Deformation Versus Bedding Plane
7 Shale Anisotropy,
8 Pore Pressure Response,
9 Interpretation of the Mohr–Coulomb Failure Parameters from the Triaxial Tests
10 Potential Applications
Md. A. Islam (&) P. Skalle
Department of Petroleum Engineering and Applied Geophysics, Norwegian University of Science and Technology, Trondheim, Norway
- Aadnøy BS, Chenevert ME (1987) Stability of highly inclined boreholes. In: Proceedings of the IADC/SPE drilling conference, New Orleans, Louisiana, 15–18 March 1987, SPE 16052
- Aadnøy BS, Hareland G, Kustamsi A (2009) Borehole failure related to bedding plane. In: Proceedings of the ARMA conference, Asheville, North Carolina, 28 June–1 July 2009
- Bradley WB (1979) Failure of inclined boreholes. J Energy Resour Technol 101:232–239
- Claesson J, Bohloli B (2002) Brazilian test: stress field and tensile strength of anisotropic rocks using an analytical solution. Int J Rock Mech Min Sci 39:991–1004
- Closmann PJ, Bradley WB (1979) The effect of temperature on tensile and compressive strengths and Young’s modulus of oil shale. In: Proceedings of the SPE-AIME 52nd annual fall technical conference and exhibition, Denver, Colorado, 9–12 October 1977, SPE 6734
- Crook AJL, Yu J-G, Willson SM (2002) Development of an orthotropic 3D elastoplastic material model for shale. In: Presented at the SPE/ISRM rock mechanics conference, San Antonio, Texas, 20–23 October 2002
- Fjær E, Holt RM, Horsrud P, Raaen AM, Risnes R (2008) Petroleum related rock mechanics, 2nd edn. Elsevier, Amsterdam, p 267
- Fjær E, Holt RM, Nes O-M (2011) The transition from elastic to non- elastic behavior. In: Presented at the 45th US rock mechanics/ geomechanics symposium, San Francisco, CA, 26–29 June 2011
- Head KH (1984) Soil laboratory testing, vol I–III. ELE Int., London
- Hoang SK, Abousleiman YN (2010) Openhole stability and solids production simulation of emerging gas shales using anisotropic thick wall cylinders. In: Proceedings of the IADC/SPE Asia Pacific drilling technology conference and exhibition, Ho Chi Minh City, Vietnam, 1–3 November 2010, IADC/SPE 135865
- Horsrud P (1998) Estimating mechanical properties of shale from empirical correlations. J SPE 16(2):68–73
- Horsrud P, Holt RM, Sonstebo EF, Svano G, Bostrom B (1994) Time dependent borehole stability: laboratory studies and numerical simulation of different mechanisms in shale. In: Presented at the Eurock SPE/ISRM rock mechanics in petroleum engineering conference, Delft, The Netherlands, 29–31 August 1994
- Horsrud P, Sønstebø EF, Bøe R (1998) Mechanical and petrophysical properties of North Sea shales. Int J Rock Mech Min Sci 35(8):1009–1020
- Holt RM, Sønstebø EF, Horsrud P (1996) Acoustic velocities of North Sea shales. In: Proceedings of the 58th EAGE conference and technical exhibition, Amsterdam, The Netherlands, 3–7 June 1996
- Holt RM, Fjær E, Nes O-M, Alassi HT (2011) A shaly look at brittleness. In: Presented at the 45th US rock mechanics/ geomechanics symposium, San Francisco, California, 26–29 June 2011
- Hofmann CE, Johnson RK (2006) Estimation of formation pressures from log-derived shale properties. J Petrol Technol 17:717–722
- Islam MdA (2011) Modeling and prediction of borehole collapse pressure during underbalanced drilling in shale. Ph.D. thesis at NTNU, Norway. ISBN 978-82-471-2599-1
- Islam MdA, Skalle P, Faruk ABM, Pierre B (2009) Analytical and numerical study of consolidation effect on time delayed borehole stability during underbalanced drilling in shale. In: Proceedings of the KIPCE 09/SPE drilling conference, Kuwait City, Kuwait, 14–16 December 2009, SPE 127554
- Jaeger JC, Cook NGW (1979) Fundamentals of rock mechanics, 3rd edn. Chapman and Hall, New York
- Johnston DH (1987) Physical properties of shale at temperature and pressure. Geophysics 52:1391–1401
- Katsube TJ, Mudford BS, Best ME (1991) Petrophysical characteristics of shales from the Scotian shelf. Geophysics 56:1681–1689
- Lockner DA, Stanchits SA (2002) Undrained poroelastic response of sandstones to deviatoric stress change. J Geophys Res 107(B12): 2353. doi:10.1029/2001JB001460
- Loret B, Rizzi E, Zerfa Z (2001) Relations between drained and undrained moduli in anisotropic poroelasticity. J Mech Phys Solids 40:2593–2619
- Maury V, Santarelli FJ (1988) Borehole stability: a new challenge for an old problem. In: Cundall PA, Sterling RL, Starfield AM (eds) Key questions in rock mechanics. Balkema, Rotterdam, pp 453–460. ISBN 90 6191 8359
- Nakken SJ, Christensen TL, Marsden R, Holt RM (1989) Mechanical behaviour of clays at high stress levels for well bore stability applications. In: Maury V, Fourmaintraux D (eds) Rock at great depth. Balkema, Rotterdam, pp 141–148
- Nes O-M, Sønstebø EF, Holt RM (1998) Dynamic and static measurements on mm-size shale samples. In: Proceedings of Eurock 98, Trondheim, Norway, Trondheim, 8–10 July 1998, pp 23–32, vol 2, SPE/ISRM 47200
- Niandou H, Shao JF, Henry JP, Fourmaintraux D (1997) Laboratory investigation of the mechanical behavior of tournemire shale. Int J Rock Mech Min Sci 34(1):3–16
- Papamichos E, Brignoli M, Santarelli FJ (1997) An experimental and theoretical study of a partially saturated collapsible rock. Mech Cohes Frict Mater 2:251–278
- Rozhko AY (2011) Capillary phenomena in partially-saturated rocks: theory of effective stress. In: Presented at the 45th US rock mechanics/geomechanics symposium, San Francisco, California, 26–29 June 2011
- Sarout J, Molez L, Gue´guen Y, Hoteit N (2007) Shale dynamic properties and anisotropy under triaxial loading: experimental and theoretical investigations. Phys Chem Earth 32:896–906
- Schmitt L, Forsans T, Santarelli FJ (1994) Shale testing and capillary phenomena. Int J Rock Mech Min Sci Geomech Abstr 31: 411–427
- Skempton AW (1954) The pore-pressure coefficients A and B. Geotechnique 4:143–147
- Sønstebø EF, Horsrud P (1996) Effects of brines on mechanical properties of shales under different test conditions. In: Proceedings of Eurock 96, Turin, Italy, 2–5 September 1996. Balkema, Rotterdam, pp 91–98
- Søreide OK, Bostrøm B, Horsrud P (2009) Borehole stability simulations of an HPHT field using anisotropic shale modeling. In: Proceedings of the ARMA conference, Asheville, North Carolina, June 28–July 1 2008
- Spaar JR, Ledgerwood LW, Hughes C, Goodman H, Graff RL, Moo TJ (1995) Formation compressive strength estimates for predicting drillability and PDC bit selection. In: Presented at the SPE/ IADC drilling conference, Amsterdam, The Netherlands, 28 February–2 March 1995, SPE/IADC 29397
- Steiger RP, Leung PK (1992) Quantitative determination of the mechanical properties of shales. SPE Drill Eng 7:181–185 Wood DM (1990) Soil behaviour and critical state soil mechanics. Cambridge University Press, Cambridge.
Open Access This article is distributed under the terms of the Creative Commons Attribution License which permits any use, distribution, and reproduction in any medium, provided the original author(s) and the source are credited.