Research Article |
Corresponding author: Ivan Belozerov ( i.belozerov@narfu.ru ) Academic editor: Aleksandr I. Malov
© 2018 Ivan Belozerov.
This is an open access article distributed under the terms of the Creative Commons Attribution License (CC BY 4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Citation:
Belozerov I (2018) Experimental determination of porosity and permeability properties of terrigenous reservoirs for creation and validation of a digital core model. Arctic Environmental Research 18(4): 141-147. https://doi.org/10.3897/issn2541-8416.2018.18.4.141
|
Digital core modelling is a vital task assessing original-oil-in-place. This technology can be seen as an additional tool for physical experiments capable of providing fast and efficient modelling of porous media. The objective of the paper is to determine experimentally the porosity and permeability properties of rocks and justify the possibility of using them for digital core modelling. The paper also validates feasibility of using the results of lithologic and petrographic surveys of thin sections in digital core modelling. The experimental studies of reservoir conditions allowed us to obtain curves of the dependence between the kerosene permeability of the terrigenous reservoir of the Buff Berea field and the temperature and to determine its main porosity and permeability properties. The paper also validates feasibility of applying the results of lithologic and petrographic surveys of thin sections of the reservoir to form the structure of the pore space of a digital core model by machine learning. The choice of this reservoir stems from the fact that the terrigenous sandstones of Berea Sandstone (USA) are characterised by minimal anisotropy of porosity and permeability properties, relatively high porosity and permeability, as well as uniformly sized grains of the composing rocks and good sorting. Oil industry experts therefore consider samples of these rocks to be most suitable for conducting applied research and testing various technologies. The results obtained were used to select the parameters required for modelling filtration flows in a digital model of the core.
terrigenous reservoirs, porosity and permeability properties, digital core model, permeability change dynamics, supercomputer, modelling of filtration processes
The technology of digital core modelling is currently a popular and developing direction in assessing the world’s original-oil-in-place (
In Russia, this technology is only just starting to develop. The first Russian laboratory for digital core modelling was opened in 2016. The works of Russian researchers on digital core modelling have been published relatively recently (
In the initial stage of creating a digital core model, the geometric model of the pore space can be presented in the form of a close packing of balls, the size of which is selected based on the data obtained during determination of the granulometric size composition of samples, microtomography and other studies. Yet rock particles are rarely shaped like balls (
Various mathematical models and calculation methods may be used for modelling the filtration processes occurring in the reservoir. The most common include: pore network models, model of lattice Boltzmann equations, smoothed particle method, finite volume methods, diffuse boundary models, as well as models based on solving the Navier-Stokes/Stokes equations (
Nevertheless, methods of finite volumes (
When forming a digital core model, methods must be used that allow the shape and size of the grains forming the rock to be determined. Subsequently, based on the presented grain sizes, a porous medium model may be formed including various processes and initial parameters, by which permeability, porosity and other properties are calculated. The initial parameters of a digital core model include the distribution of the particle sizes of rocks, their placement and sorting, chemical interaction parameters, etc.
Various approaches and methods have been used in foreign studies to obtain the initial parameters of the digital core model. For instance, the study (
The study (
Despite the existing advantages of using 3D tomography and X-ray structural analysis of core samples and petrographic thin sections for obtaining initial parameters when building the volumetric structure of rocks, these approaches have some disadvantages. One of their main problems is that they fail to provide a complete description of and substantiation for the micro-processes occurring in the reservoir. This may be partially resolved by increasing the resolution of the images during tomographic studies of core material and accurately determining the granulometric size composition of rocks by laser diffraction.
In order to create and validate the digital core model, the author of this paper proposes using a set of parameters obtained from studying core material and lithologic and petrographic surveys of thin sections. In this case, the digital core model may be presented as a close stochastic packing that will form the microstructure of the core and it is then possible to apply the methods of molecular dynamics. Next, in the process of mathematical modelling of filtering processes, the results obtained on a supercomputer are compared with the results of physical experiments, making it possible to evaluate the effectiveness and accuracy of the model obtained.
An important step in creating the digital core model is to obtain the results of experimental studies of physical samples required for comparison with the results of computational experiments on a supercomputer. In order to obtain the experimental data required for subsequent studies, the main porosity and permeability properties of the terrigenous reservoir, as well as the kerosene permeability coefficients at various temperatures in reservoir conditions, were determined in the paper. The studies were conducted at the Innovative Technological Centre “Arctic Oil and Gas Laboratory Research” of Northern (Arctic) Federal University named after M.V. Lomonosov.
Buff Berea sandstone was selected as the terrigenous reservoir. This reservoir was selected owing to its homogeneity in terms of material composition and its porosity and permeability properties, as well as its wide recognition by oil industry experts in conducting applied research and testing of various technologies.
The study subject is Buff Berea sandstone, representing material for synthesising a digital core model of terrigenous reservoirs that might potentially accumulate hydrocarbons.
Ten standard-sized core samples were selected for study. Regular-shaped samples of rock were taken by drilling, trimming, surfacing and grinding in accordance with the requirements of GOST 26450.0-85 (Rocks. General requirements for selection and preparation of samples for determining reservoir properties. Date of issue: 1986–07–01.). Next, absolute permeability coefficients were determined for stationary filtration and linear gas flow in accordance with the requirements of GOST 26450.2-85 (Method for determination of the absolute gas permeability coefficient for stationary and non-stationary filtration. Date of issue: 1986–07–01.). Nitrogen was used as the gas.
The open porosity coefficients were determined according to GOST 26450.1-85 (Method for determination of the open porosity coefficient by fluid saturation. Date of issue: 1986–07–01.). Non-polar kerosene was used as the saturating fluid.
At the next stage, the samples were loaded into a UIK-5 (7) filtration unit (Fig.
To saturate the samples, the MS-535 manual saturator (Fig.
A UIP-ATM.002 training and research permeameter was used to determine the absolute permeability coefficients (Fig.
The results of determining the coefficients of open porosity and absolute permeability of the samples are presented in Table
Results of determination of coefficients of open porosity and absolute permeability of the samples
Sample No. | Length, cm | Diameter, cm | Absolute permeability coefficient, 10-3 µm2 | Open porosity coefficient, % |
г-1 | 3.068 | 3.003 | 101.01 | 20.11 |
i6 | 2.682 | 2.997 | 125.64 | 20.42 |
i1 | 3.110 | 2.999 | 130.25 | 17.99 |
i2 | 2.722 | 2.996 | 96.43 | 18.89 |
i4 | 2.743 | 3.002 | 85.82 | 17.48 |
i7 | 2.854 | 3.003 | 91.64 | 18.55 |
229 | 2.974 | 2.995 | 121.13 | 19.63 |
1 | 3.580 | 2.995 | 101.032 | 18.07 |
i8 | 2.967 | 2.997 | 97.22 | 18.52 |
2 | 3.630 | 3.003 | 83.64 | 18.79 |
As can be seen in Table
For mathematical modelling of the core, the reservoir temperature is an important factor influencing the values of the permeability coefficient. This is confirmed by studies conducted on core samples of the Buff Berea field. The graph of variance of the permeability of core samples depending on the temperature in reservoir conditions is presented in Figure
In general, the dynamics of changes in kerosene permeability depending on the temperature in reservoir conditions shows a decrease in permeability as the temperature increases for a given collector. By selecting the function giving the most accurate description of the nature of the curves, a quintic polynomial function was chosen and equations were derived for each of the ten dependencies. Next, a generalised expression (1) of the dependence of kerosene permeability on the temperature in reservoir conditions was obtained for these rocks.
y = 0.000001x5 − 0.0002x4 + 0.203x3 − 0.8642x2 + 17.3x − 54.55 (1)
where, у – kerosene permeability of the core, 10-3 µm2; х – temperature, °С.
The resulting dependence may differ significantly for other terrigenous reservoirs. The reason is that the structure of terrigenous rocks has a huge number of pore channels of complex geometric shape and the rocks may differ significantly in their material composition over the entire length of the formation. The presence of capillary forces, contacts of rock-forming minerals with oil, that differ in strength and nature of interaction, and many other factors make it difficult to forecast the properties of terrigenous reservoirs that are characterised by significant differences. This is confirmed in a number of papers (
When creating a close stochastic packing of a digital core model to bring the shapes of the grains of the modelled rock closer to the real ones, along with the results of determining the granulometric size composition of rocks, it is proposed to use the results of lithologic and petrographic surveys of thin sections. These results may later be used to determine the geometric parameters of the modelled particles by machine learning, which will allow extracting context-sensitive values of the lexical units from photographs and apply latent-semantic analysis to form a stochastic packing of the digital core model. This will have a positive impact on the approximation of the parameters of the generated digital core model to natural samples.
In the beginning stages of these studies, the initial image of the petrographic thin section of the terrigenous sandstone Berea Sandstone was transformed in order to establish an algorithm that could make it possible to quickly and accurately determine the characteristics of the microstructure of the rock pore space for digital core modelling. It has been found that, owing to compaction of the rock under study, it was difficult to determine the parameters of the grain boundary section. To solve this problem, the Hough transformation was applied in the OpenCV library of computer analysis algorithms of the thin section under study, making it possible to obtain the image shown in Figure
Next, a threshold transformation of the thin section image was performed. To capture the brightness of all the grain particles, an adaptive threshold transform was used, which considered objects not in one pixel but in its vicinity. The position of the particles was determined by converting the brightness of the pixel in its vicinity through weighted coefficients in accordance with the Gaussian function. The result of the threshold transformation of the thin section is shown in Figure
The threshold transformation made it possible to determine better the position of the grains in the microstructure of the rock under study. The next step was morphological transformation of the image required for determining the area of the voids in the thin section. At this stage, the colour of the rock grains was replaced with black and the voids with white (Fig.
At the last stage, the black and white pixels were counted and the number of white pixels divided by the number of black ones, which corresponds to the petrographic method for determining porosity.
There are many proven software tools for high-quality and fast determination of the structure of the voids in petrographic thin sections using the threshold image transformation algorithm considered by the author. This algorithm may be applied in digital core modelling to form the structure of the pore space of the rock.
This paper presents the results of determining open porosity and absolute permeability (20 definitions), as well as kerosene permeability coefficients, at various temperatures in reservoir conditions (100 definitions). The studies were conducted according to the algorithm of threshold-transformation of petrographic thin sections images proposed by the author.
The results presented in the paper in the form of the established equation of the dependence of kerosene permeability on the temperature in reservoir conditions allows percolation analysis to be applied for forecasting the permeability of the digital core model. The results of lithologic and petrographic surveys of thin sections may be used to determine the geometric parameters of the modelled particles of stochastic packing of the digital core model using machine learning.