Experimental study on the effect of rock pressure on sandstone permeability circular. Proceeding from JOURNAL

10.31897/PMI.2022.24 Abstract. The results of laboratory studies to determine the effect of effective stress on the permeability of sandstone are presented. During the test, the samples were subjected to a stepwise increase or decrease of the effective stress (at a constant pore pressure) in a specified step. The values of rock permeability at different values of effective stress were determined, and the influence of the grain size of the reservoir rock matrix on the character of the change in the sandstone permeability coefficient was also established. During the test, a decrease in permeability was observed with an increase in effective stress. It was found that as a result of gradual loading/unloading of the sandstone sample, the original permeability values were not restored, which indicates the beginning of the formation of residual strains in the rock. This effect should be taken into account when modeling field development because in the process of reserves extraction the effective stress acting on the reservoir rock skeleton changes, which results in a significant chang in rock permeability. The results of laboratory studies showed that the deviation of permeability in medium-grained sandstones relative to the initial value was greater than in medium- and fine-grained sandstones. The pressure sensitivity coefficient and constant of material, which are used in empirical relationships between permeability and effective stress, were numerically estimated. At the same time, the constant of material showed no such convergence, which indicates that the values of this parameter are individual for each


Introduction.
Designing the development of an oil and gas field is a complex process, consisting of several stages. Creating hydrodynamic and geological models of a field is one of the important stages. Model adequacy influences accuracy and reliability of predictions for key figures that determine efficiency of oil and gas fields' development, completeness of reserves depletion.
In order to create hydrodynamic and geological model of a formation, it is necessary to carry out a large amount of research work, including laboratory investigation. One of the main types of research is filtration studies, since fluid flow through the pore space is one of the key processes occurring in the formation. Fluid movement in porous media depends on many parameters: porosity and permeability of the reservoir, pressure gradients in the formation, formation pressure, effective stress, etc. [1]. A large number of works are dedicated to the simulation of fluid flow through porous media. Many existing models are based on Darcy equation [2][3][4].
The main task in selecting a model is to describe the investigated physical process accurately and precisely, so modern fluid flow models are multi-level systems of equations with many variables, which have to take into account a number of factors [5][6][7]:  Parameters that take into account well imperfection by nature and degree of drilling-in.  Calculations by Dupuis formula assuming a radial feed contour geometry. As the practice of developing the oil fields shows, the shape of the feed contour is not always circular. Proceeding from this fact and the fact that the shape of the feed contour is one of the factors influencing fluid inflow, its geometry has to be taken into account [8].
 When simulating the process of fluid flow it is necessary to take into account changes in permeability, as it is one of the main parameters characterizing permeability and porosity properties of reservoir rocks and well productivity. The value of this parameter is influenced by many processes, including change of the stress-strain state (SSS) of rocks [9][10][11][12], which depends on physical and mechanical properties of rocks [13][14][15], rock and pore pressures.
Many researchers examine changes in physical and mechanical properties in different geological conditions [16][17][18]. The authors of [10] performed a laboratory study of changes in permeability of IBS (interbedded shale) samples in pressure and temperature conditions similar to those that exist during SAGD (Steam Assisted Gravity Drainage) operations and observed that Young modulus is temperature dependent, whereas Poisson coefficient and compressibility of IBS samples are not.
A number of studies have been carried out to determine the durability characteristics of sedimentary rocks [19][20][21]. Custom designs and equipment have been used to investigate deformation dynamics, sandstone creep, pore collapse and rock destruction [22][23][24].
In [9] the samples were tested for uniaxial compression and tension with the determination of Poisson coefficient and Young modulus. The samples were divided into two groups: from waterflooded and non-waterflooded formations. The research results showed that in samples taken from waterflooded formations, the values of Poisson coefficient increase with a simultaneous significant decrease in the modulus of elasticity and durability limit in uniaxial compression and tension.
Work [14] notes the effect of saturating fluid on the basic physical and mechanical properties of the rock. During the experiments, the samples were divided into three groups: dry, saturated by oil and formation water (3 % KCl).
The analysis of conducted studies has shown that fluid saturation of sandstone leads to a decrease in its durability properties. The greatest decrease in uniaxial compressive and tensile durability limit is observed with water saturation. The presence of fluid (oil or formation water) in the pore space of a rock leads to a decrease in the dynamic modulus of elasticity. Presence of oil in the reservoir rock has the greatest influence on this characteristic.
A large number of works have been devoted to identifying the effect of pore and rock pressure on the SSS [25][26][27]. When describing the dependence of pore pressure and effective stress, the contribution of K.Terzaghi [28] should be taken into account. He proposed a model to calculate the clay permeability coefficient, which describes the influence of fluid flow and changes in SSS: where ij -effective stress; Sij -full stresses; Ppor -pore pressure.
The model proposed by Terzaghi was further developed by M.Biot. He derived the relationship between the rock skeleton stress and the pressure of the fluid that fills the pore space. Biot poroelasticity model has been reviewed in many works [29][30][31]. Biot's law is presented in [32-34]: Most of modern fluid flow models do not fully take into account rock mechanics and the relationship between the main parameters characterizing the SSS with the main permeability and porosity properties. This aspect is relevant in modern oil and gas industry, as more precise correction of fluid flow law in pore space when designing hydrocarbon field development will allow to predict technical and economic indicators of development more accurately.
Since permeability is one of the main parameters characterizing the permeability and porosity properties of reservoir rocks, there is a need to establish a relationship between effective stress and rock permeability.
The objective of this research is to determine the relationship between effective stress and rock permeability for accurate prediction of the parameters of hydrocarbon field development depending on the stage. In order to achieve the goal, the following tasks were set: determination of the dependence of reservoir rock permeability on effective stress, the relationship between reservoir rock granularity and the nature of permeability changes.
Methodology. Most oil and gas fields in Russia are confined to terrigenous rocks with complex geological structures, so sandstone from the West Siberian field was used as samples: 1-4 -light gray sandstone, medium-and fine-grained, weakly mica, with clay cement (dense, firm, massive); 5-6 -light gray sandstone, medium-grained, weakly mica, with clay cement (dense, firm, massive).
Before the test, the samples were prepared on stonecutting equipment to ensure parallelism of the sample surfaces. The deviation from parallelism was measured with an indicator mounted on a tripod in two mutually perpendicular directions; the diameter dimension should be no more than 0.2 mm (GOST 21153.8-88).
After checking the parallelism of the sample surfaces, the next preparation step was the cleaning of the rock samples, which included:  cleaning of samples from residual fluid (water, oil, etc.), which was achieved by washing, leaking or contacting with various solvents;  drying of all samples to reduce the natural moisture content (each rock sample was dried until its mass became constant).
The last preparation step involves measuring the height, diameter and mass of each sample. The diameter was measured at three locations along the height of the sample (in the middle and at the end surfaces) in two mutually perpendicular directions. The difference in diameter for these measurements did not exceed 0.5 mm. The arithmetical mean of all measurements (GOST 21153.8-88) was taken as the calculated value. Six samples were examined in total.
The following methodology was chosen for the experiment: 1. The sample was tested in the gas permeability analyzer PIK-PP, which allows to determine the permeability values in the range from 0.01-10 −3 to 10 000 μm 2 as well as to set rock pressure up to 70 MPa. The operating principle of the analyzer is based on Boyle -Marriott law, where the basic parameter of rock permeability is determined by helium flow rate through the tested sample at constant pressure (Ppor = 0.1 MPa). To obtain more accurate results, the permeability was measured three times at each point; the test time per sample was 10 min.
2. Measuring the initial permeability value k1 at 10 MPa, t = 25 °C. 3. Stepwise increase in rock pressure (with permeability measurements at each loading stage) in steps of 10 to 40 MPa (maximum rock pressure value was determined based on bedding conditions).
4. When the critical rock pressure was reached, the lateral pressure was reduced to 10 MPa with a step of 10 MPa. 5. After completing stage 4, the sample was removed from the PIK-PP gas analyzer, replaced by another sample and stages 1-4 were repeated.
6. To determine the effect of rock pressure magnitude on rock permeability, dependence plots were constructed (Fig.1).
Discussion of the results. The analysis of the research results allowed establishing a correlation between one of the parameters characterizing the effective stress and the change in permeability. The results of the studies on the effect of the effective stress (at the assumed 1   ) on the permeability coefficient are shown in Fig.1 and Table 1.
In most models describing fluid flow in pore space, permeability is assumed to be a constant, although the results of these and other studies indicate that permeability values may vary depending on the stress acting on the rock skeleton. Consequently, permeability is a function of stress. It should be noted that depending on permeability and composition of reservoir rock, the effect of rock pressure on permeability values will vary. In some works, equations are given to establish the relationship between effective stress and rock permeability, but there is a need for additional definition of the parameter, as it depends on the type of rock [35]: where k -permeability at effective stress Peff; k0 -permeability at pressure P0; γ -pressure sensitivity coefficient.
The dependence [36] can be used to determine the relationship between effective stress and permeability: where p -constant of the material. According to studies, the pressure sensitivity coefficient γ is in the range of 1.37·10 −3 -18.1·10 −3 MPa −1 for different sandstones [35,36]. The measured values of the pressure sensitivity coefficient γ given in Table 2 are very close to those given in [36], while the material constant p differs from that presented in this paper.  (5) and (6) are individual for each field. Consequently, determination of the relationship between effective stress and permeability must be carried out at the stage of pilot field operation, that is, at the initial stage of oil or gas field development, in order to provide a more accurate description of processes occurring in the reservoir and to accurately predict the main development parameters.
Conclusion. Conducted studies allowed to establish relationships between effective stress and change in permeability, which allow more accurate prediction of key development indicators (e.g., well productivity). In addition, analysis of the obtained dependencies showed that the determination of the relationship between sandstone permeability and rock pressure can be derived from the following dependencies:  (3), if pore pressure changes significantly (e.g., by increasing or decreasing water injection pressure), the reservoir rock skeleton may be destroyed, which will lead to irreversible permeability change. The permeability change at each point was determined relative to the initial value obtained at the beginning of the test (k1 at Peff = 10 MPa). In the process of increasing or decreasing the effective stress acting on the sample, hysteresis (permeability value did not return to its initial value) was observed. In samples 1, 2 and 4 the permeability decrease was uneven. For example, a sharp change of permeability in samples 2 and 4 is observed in the pressure range of 20-30 MPa that can be explained by possible repacking of pore space in the samples.
In hydrodynamic modelling, permeability should be a function of stress. The proposed equations (5) and (6) can be used as a simplified relationship for express estimation in development monitoring. However, there is a need for additional definition of pressure sensitivity coefficient γ and constant of material p, as these parameters are individual for each rock type. Despite the fact that the pressure sensitivity coefficient γ from the presented laboratory studies is very close to the values given in [36], it is recommended to perform laboratory rock samples tests in order to identify the values of this parameter for a specific type of material.
The nature of the change in the permeability coefficient of sandstone with changes in SSS depends on the skeleton grain size, type and composition of the cementing material. The initial permeability values of the sample do not recover when the effective stress decreases that indicates the beginning of residual deformation of the rock. As shown in fig. 1, medium-grained sandstones undergo greater irreversible deformation than medium-and fine-grained sandstones. For medium-grained sandstones, the permeability range is 0.97-2.1 %, for medium-and fine-grained sandstones it is 0.63-0.85 %.