CC BY
33
UDK 622.2
Mechanical properties of sandstone using non-destructive method
Herimitsinjo RAJAOALISON1, Albert ZLOTKOWSKI1», Gérard RAMBOLAMANANA2
lAGH University of Science and Technology, Krakow, Poland
2 Institute and Observatory of Geophysics ofAntananarivo, Antananarivo, Madagascar
The understanding of physical and mechanical properties of rock is considered as critical in drilling, geo-engineering, and construction applications. As an example, the awareness of these rock parameters contributes to avoid or minimizing instability around the wellbore while drilling. The laboratory experiment of understanding of these parameters can be done in two-different ways: static, where the sample subjects to destruction after the test and dynamic, known as non-destruction method. The non-destructive method using ultrasonic waves under a series of different stress conditions, starting from 7 to 56 MPa with incrementation of 7MPa, has been used in this paper in order to characterize the mechanical properties of dry Zbylutôw sandstone at 20 and 80°C. The velocity of primary (P) and secondary (S) waves within these ranges has been recorded in order to understand the behavior of the mechanical properties. The results showed that the Young's modulus, bulk modulus, shear modulus, and Lame's constant of Zbylutôw sandstone have a positive correlation with good coefficient correlation with the increased stress, while the Poisson's ratio showed a negative correlation. Besides, the effect of temperature on the rock parameters is approved by the decrease of primary wave velocity in this two-different temperature range. Such results are necessary when preparing the appropriate mud weight for drilling process, which is related to wellbore instability.
Key words: mechanical properties; non-destructive method; sandstone; stresses; temperature effect
Acknowledgment. The research was supported by the AGH subsidy 16.16.190.779 and AGH-UNESCO.
How to cite this article: Rajaoalison H., Zlotkowski A., Rambolamanana G. Mechanical properties of sandstone using non-destructive method. Journal of Mining Institute. 2020. Vol. 241, p. 113-117. DOI: 10.31897/PMI.2020.1.113
Introduction. All stages of drilling process such as drilling designing, well completion, and production (не является частью процесса бурения) require a better understanding of the mechanical properties of surrounding rock such as elastic modulus and rock compressive strength, in which contributes mainly to the instability of the wellbore as well as inducing fractures [2, 8]. For instance, it is well known that collapsing is due to misunderstanding of these parameters [3, 6]. Measuring of such parameters can be conducted in two different ways: static and dynamic measurement, in which the static measurement can be done in lab where the sample is loaded by stresses until failure occurs, or by in-situ tests, while the dynamic measurement, non-destruction method, is carried out in lab without destroying the samples, in which the sonic waves are used to perform such experiment [17]. Compressional (Vp) and shear waves (Vs) velocities are measured from logs or in lab and used to calculate the dynamic mechanical properties of rock [10, 11].
Numerous studies confirm the difference between the results from these two methods. The static measurement is the more realistic, direct, and very close to the actual engineering needs as it represents the occurred deformation under large amount of stresses, which refers as the underground high stresses condition, while the dynamic measurement is more continuous and easier [5, 9]. Therefore, dynamic moduli are quite different than static moduli but both techniques are required in order to well understand the mechanical parameters of the rock formation. As usual, the static moduli are smaller than the corresponding dynamic. Such difference is depending also on the rocks state. It is larger for weak rocks and decreases with increased confinement. However, Hao Xu et al. has established a relationship between dynamic and static measurement, which can be used in case the static measurement cannot be conducted in the laboratory or missing equipment [16].
While drilling, the in-situ stress distribution around the wellbore is being disturbed and changed, which contributes to the instability of the wellbore. So, understanding of the rock behavior under stress and temperature is a key point to addressing such challenge.
Backgrounds and methods. Mechanical properties of rock are characterized by the elastic as well as inelastic properties. The elastic properties consist of Young's modulus, shear modulus, bulk modulus, and Poisson's ratio; while fracture gradient and formation strength represent the inelastic properties. The elastic properties concern with the resistance to volume and shape deformation and it follows the Hooke's law within small deformation. This deformation is depending on the driven acting force from the geological situation, its nature, and the following elastic moduli such as Young's modulus E, (2) and shear modulus p, (3). The estimation of the dynamic mechanical properties of rock can be done by using the acoustic waves and considering the rock sample as isotropic and linear elastic [1]. Then, the Vp/Vs ratio and density, p, are used to calculate the mechanical properties of the rock sample [8]. The following formula express the relationship between the acoustic waves velocity, density and the mechanical properties of rock:
f S1 1 _| Vl
V = —^, (1)
f K ^ v vp ,
2 1 -
Fp2(1 + v)(1 _ 2V)
E = P „ \-", (1)
(1 _ V)
p = pV2, (3)
X = P( Vp2 _ 1V,2), (4)
AI = pVp. (5)
where Vp is the velocity of P-wave, m/s; Vs is the velocity of ¿'-wave, m/s; p is the bulk density, kg/m3; v is the Poisson's ratio; E is the Young's modulus, Pa; |i is the shear modulus, Pa; X is the Lame's constant, Pa; AI is the acoustic impedance, kgm-1s-1.
Twelve samples from Zbylutow, located in the southern part of Poland, have been used during the experiment. The moisture trapped in the samples was removed using an electrical oven at 70°C for more than 14 hours. Next the density of each sample was measured. The density of dry sample was calculated from the mass (in g) and bulk volume (in cm3) of the sample. Both ends of the samples were polished and lapped parallel up to ± 0.10 mm in order to minimize the error and protect the sample from damage due to higher stress. Also, the cylindrical sides are made straight active with accuracy of ± 0.30 mm over the full length of each sample.
The AVS 1000 apparatus was used to conduct the dynamic measurement (Fig. 1). It is composed of parametric pulsar, tektronix digitizing oscilloscope, and sonic platens containing P and S-wave transducers as well as receivers. The time-of-flight of P-wave and ¿-wave is recorded and their velocities are calculated accordingly as the ratio between the sample length and the travel time [11]. The ultrasonic measurement uses high frequency signals, approximately 1000 kHz for compres-sional and 500 kHz for shear waves. The system was calibrated using aluminum and copper samples with well-known velocities. The arrangement of AVS 1000 is shown by the following diagram.
Results and discussions. The density of Zbylutow sandstone is ranging from 1.34 to 1.48 g/cm3 with a median value of 1.43 g/cm3, whereas the porosity of samples is ranging from 1.80 to 4.95 % with a median value of 4.13 %. The results in this paper represented the median of the investigated properties. The median value has been chosen by the fact that it is not affected by the outlier values. The correlation between the mechanical properties and stress was expressed by the central tendency of the measured values within the ranges of the independent values.
Generator
Pulse/Receiver
Oscilloscope
Switch
(Pore pressure)
(Gas tank)
P
P
Data acquisition
(Radial pressure)
(Axial pressure)
Hand pump
S
S
Fig. 1. Scheme of AVS 1000 apparatus
The measurement of mechanical properties was conducted at two-different temperature and carried out in-situ level of stresses. The result (Fig.2) showed that the velocity of P and ¿'-waves were proportionally increasing with stresses, but only P-waves velocity inversely proportional to temperature. The temperature caused weaknesses to the sample. On the contrary, the increased stress was closing the crack inside the sample [15]. Also, the acoustic impedance profile also showed the increase of compressional velocity of the rock against stresses (Fig.3, f).
As bulk modulus was measuring resistance to compression, then it has been identified that increasing in stress resulted increasing of rock's resistance to compression. As shear modulus represented the ratio of shearing (torsional) stress to shearing strain, then increasing in stress increased the shearing inside the rock. The behavior of Poisson's ratio showed that geometric change of shape under stress was decreasing as long as stress increasing.
Young's modulus, bulk modulus, shear modulus, and Lame modulus of sandstone increase with stresses and decrease with temperatures and the rate of change varies with stress and temperature. The cracks being closed can be used as an explanation of the increase [13, 4], while the pore volume change can be used as an explanation of the decrease. The pore volume is increasing (even small portion) when the temperature raises and the differential expansion of temperature applied to minerals, forming the rock, may generate new cracks especially at the boundary of grains [7, 14].
Conclusions. The dynamic measurement is conducted in this paper in order to understand the impacts of raised stress and temperature to the mechanical properties of the dry Zbylutow sandstone. Based on the results:
1. Young's modulus, bulk modulus, shear modulus, and Lame's constant of Zbylutow sandstone were proportional to increased confining stress and inversely proportional to increased temperature, while the Poisson's ratio is inversely proportional to increased stress and proportional to increased temperature.
u >
6000 -I
5000 -4000 3000 H 2000 1000
0 10 20 30 40 Stress, MPa
=t-„
50
60
P-wave 20 °C P-wave 80 °C
S-wave 20 °C S-wave 80 °C
Fig.2. Velocity of P and S waves versus confining stress
0.42 |0.415
0 0.41 1 'I 0.405 | 0.4 -j
1 0.395 -
^ 0.39 -
«... • '
V..
y =
•••• \\
-0.0007x + 0.4212 R2 = 0.956
0.0005x + 0.4216 R2 = 0.9911
10 20 30 40 Stress (MPa)
50 60
3.5E + 10 3.3E + 10 3.1E + 10 2.9E + 10 -2.7E + 10 -2.5E + 10 -
2.3E + 10 -2.1E + 10 1.9E + 10
0
= 3E + 0.8x + 2E + 10 R2 = 0.9882
y =
2E + 0.8x + 2E + 10 R2 = 0.9989
10
20 30 40 Stress, MPa
50 60
5.1E + 10 4.9E + 10 -4.7E + 10 -4.5E + 10 4.3E + 10 4.1E + 10
0
: 1E + 0.8x + 4E + 10 R2 = 0.9554
9554
.....«-:!...-........
y =
2E + 0.8x + 4E + 10 R2 = 0.9254
10
20 30 40 Stress, MPa
50 60
2
o
K
1.3E + 0.9 ■ 1.2E + 10 ■ 1.1E + 10 ■ 1E + 10 9E + 09 ■ 8E + 09 ■ 7E + 09
: 1E + 0.8x + 7E + 09 R2 = 0.9883
y = 8E + 0.7x + 7E + 09 R2 = 0.9989
10
20 30 40 Stress, MPa
50 60
f
4.3E + 10 -4.2E + 10 -4.1E + 10 -
4E + 10 3.9E + 10 3.8E + 10 H 3.7E + 10 3.6E + 10
0
= 4E + 07x + 4E + 10 R2 = 0.6373
...........•
y
= 2E + 08x + 4E + 10 R2 = 0.8747
10
20 30 40 Stress, MPa
50 60
3
u -ft
M M
=1
o o <i
1.28E + 07 1.26E + 07 1.24E + 07 1 1.22E + 07 1.2E + 07 H 1.18E + 07 1.16E + 07 1.14E + 07 1.12E + 07
= 23002x + 1E + 07 R2 = 0.999
- 28112x + 1E + 07 R2 = 0.9679
0
10
20 30 40 Stress, MPa
50 60
20 °C
i 80 °C
Fig.3. dynamic: a - Poisson's ratio; b - Young modulus; c - bulk modulus; d - shear modulus; e - Lame's constant; f - acoustic impedance versus confining stress at 10 °C and 80 °C
b
a
0
d
c
0
e
2. The acoustic impedance confirmed that raising of stresses resulted increasing of compres-sional velocities, while increasing of temperature decreased the compressional waves velocity through the rock samples.
3. The present result is necessary when formulating the mud weight windows program for drilling and allow the drilling operator to understand the behavior of the sandstone formation against stress and temperature.
REFERENCES
1. Bailin W. Boit's Effective Stress Coefficient Evaluation: Static and Dynamic Approaches. ISRM International Symposium -1nd Asian Rock Mechanics Symposium. International Society for Rock Mechanics and Rock Engineering. 1001.
1. Baoping L., Hongzhi B. Advances in calculation methods for rock mechanics parameters. Petroleum Drilling Techniques. 1005. Vol. 33. N 5, p. 44.
3. Chang C., Zoback M.D., Khaksar A. Empirical relations between rock strength and physical properties in sedimentary rocks. Journal of Petroleum Science and Engineering. 1006. Vol. 51. N 3-4, p. 113-137.
4. Cheng C.H., Toksoz M.N. Inversion of seismic velocities for the pore aspect ratio spectrum of a rock. Journal of Geophysical Research: Solid Earth. 1979. Vol. 84. Iss. B13, p. 7533-7543.
5. Fjar E., Holt R.M., Raaen A.M., Risnes R., Horsrud P. Petroleum related rock mechanics. Elsevier. 1008. Vol. 53, p. 514.
6. Han D., Peng S. Investigation and future research strategy on engineering disasters induced by coal mining at deep depth in China coal industry. China Coal. 2002. Vol. 28. N 2, p. 5-9.
7. Hellwege K-H. Physical properties of rocks: Berlin-Heiddberg. Springer-Verlag. 1982. Vol. 1, p. 476.
8. Guo J., Liu Y. A comprehensive model for simulating fracturing fluid leakoff in natural fractures. Journal of Natural Gas Science and Engineering. 2014. Vol. 21, p. 977-985.
9. Mavko G., Mukerji T., Dvorkin J. Rock physics handbook-Tools for seismic analysis in porous media. Cambridge University Press. 2009, p. 524. DOI: 1017/CB09780511626753
10. Rajaoalison, H., Knez, D., Zlotkowski, A. Changes of dynamic mechanical properties of brine-saturated Istebna sandstone under action of temperature and stress. Przemysl Chemiczny. 2019. Vol. 98. N 5, p. 801-804.
11. Ranjbar-Karami R., Shiri M. A modified fuzzy inference system for estimation of the static rock elastic properties: A case study from the Kangan and Dalan gas reservoirs, South Pars gas field, the Persian Gulf. Journal of Natural Gas Science and Engineering. 2014. Vol. 21, p. 962-976.
12. Rao M.V.M.S., Prasanna Lakshmi K.J. Shear wave propagation in rocks and other lossy media: An experimental study. Curr. Sci. 2003. Vol. 85. N 8, p. 1221-1225.
13. Toksoz M.N., Cheng C.H., Timur A. Velocities of seismic waves in porous rocks. Geophysics. 1976, Vol. 41(4), p.621-645.
14. Zhang J.J., Bentley L.R. Change of bulk and shear moduli of dry sandstone with effective pressure and temperature. CREWES Res Rep. 1999. Vol. 11, p. 1-16.
15. Walsh J.B. The effect of cracks on the compressibility of rock. Journal of Geophysical Research. 1965, Vol. 70(2), p. 381-389.
16. Xu H., Zhou W., Xie R., Da L., Xiao C., Shan Y. and Zhang H. Characterization of rock mechanical properties using lab tests and numerical interpretation model of well logs. Mathematical Problems in Engineering. 2016. Vol. 5. DOI: 10:1155/2016/5967159
17. Yuming S. and Guowei L. Experimental study on dynamic and static mechanics parameters of rocks under formation conditions. Journal of Chengdu University of Technology. 2000, Vol. 27. N 3, p.249-253.
Authors: Herimitsinjo Rajaoalison, PhD candidate, rajaoalison@agh.edu.pl (AGH University of Science and Technology, Krakow, Poland), Albert Zlotkowski, adjunct, azlot@agh.edu.pl (AGH University of Science and Technology, Krakow, Poland), Gérard Rambolamanana, Professor and Leading Researcher, g_rambolamanana@yahoo.fr (Institute and Observatory of Geophysics ofAntananarivo, Antananarivo, Madagascar).
The paper was received on 18 November, 2019.
The paper was accepted for publication on 10 January, 2020.
