Next Article in Journal
Pore-Fracture System Distribution Heterogeneity by Using the T2 Spectral Curve under a Centrifugal State
Next Article in Special Issue
Effects of Compositions and Fractal Pores on CO2 Adsorption in Lacustrine Shale
Previous Article in Journal
Antibacterial Properties of Dandelion Extract-Based PVA/CTS/DAN/CuNP Composite Gel
Previous Article in Special Issue
The Strike-Slip Fault System and Its Influence on Hydrocarbon Accumulation in the Gudong Area of the Zhanhua Depression, Bohai Bay Basin
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Method for Predicting Bound Water Saturation in Tight Sandstone Reservoirs Using Morphology and Fractal Models

1
The Second Exploration Team of Shandong Coalfield Geological Bureau, Jining 272000, China
2
College of Earth Sciences & Engineering, Shandong University of Science and Technology, Qingdao 266590, China
3
Department of Mechanical, Materials and Manufacturing Engineering, Faculty of Engineering, University of Nottingham, Nottingham NG7 2RD, UK
*
Authors to whom correspondence should be addressed.
Processes 2024, 12(9), 1811; https://doi.org/10.3390/pr12091811
Submission received: 21 July 2024 / Revised: 13 August 2024 / Accepted: 19 August 2024 / Published: 26 August 2024

Abstract

:
The nuclear magnetic resonance T2 spectrum was used to identify the T2 cut-off value, which is the key to determining the irreducible water saturation of a reservoir. In this paper, the saturation and centrifugal T2 spectra of sandstone and coal samples were used to explore the correlation between each parameter and the T2 cut-off value, using a single fractal dimension, a multifractal dimension and a spectral morphology discrimination method. The conclusions are as follows: (1) The T2 spectra of nine sandstone samples in this paper can be divided into four types. Type A is represented by sample 2, wherein the T2 spectrum shows a bimodal state and the area of the right T2 spectrum (2.5~100 ms) is larger than that of the left T2 spectrum (T2 < 2.5 ms), indicating that the sample has good pore connectivity and belongs to the macroporous development sample. The B-type T2 spectrum is unimodal, and the pore connectivity is poor, indicating that it is a large-pore development sample. The T2 spectrum of the C-type sample is unimodal, and the pore connectivity is very poor, indicating that it is a mesoporous development sample. The T2 spectrum of the D-type sample shows a single peak state, and the main T2 is distributed within 0.1~2.5 ms. The pore connectivity is very poor, which indicates that it belongs to the small pore development type sample. (2) The single fractal model shows that, compared with other single fractal parameters, D2 increases with the increase in the T2 cut-off value, but the correlation is weak. Therefore, it is not feasible to predict the T2 cut-off value using the single fractal dimension parameter. (3) The multifractal model shows that D−10D10 increases linearly with the increase in D−10D0, but there is no obvious linear correlation between D0D10 and D−10D10, indicating that the low pore volume area in this kind of sample controls the overall heterogeneity of pore distribution. (4) The related parameters affecting the T2 cut-off value include D−10D10, D−10/D10, D−10D0, TM and D2. Therefore, based on the above five parameters, a T2 cut-off value prediction model is constructed. The T2 cut-off value calculated by the model is highly consistent with the experimental value, which proves the reliability of the model.

1. Introduction

With the continuous intensification of energy exploration and development and energy consumption, conventional oil and gas resources in China are no longer sufficient to meet the needs of social development, and the development of unconventional oil and gas resources is imperative. Tight sandstone gas has become one of the most important unconventional energy sources in China, and the development of pore and fracture structures in tight sandstone reservoirs is a key factor restricting the development of tight sandstone gas. Low-field nuclear magnetic resonance (LF-NMR) technology has the advantages of being non-destructive, fast and highly precise, and is widely used in fluid identification, pore fracture structure characterization and other fields [1]. Nuclear magnetic resonance logging curves can be used to distinguish movable fluid and immovable fluid, and the key parameter used to distinguish these is the T2 cut-off value. The T2 cut-off value is a key parameter in nuclear magnetic resonance logging, which reflects the boundary between the movable fluid and the bound fluid in the reservoir. In nuclear magnetic resonance measurements, T2 relaxation time (also known as transverse relaxation time) is closely related to the fluidity of fluids in pores. Longer T2 values correspond to larger pores and more fluid flowing, while shorter T2 values indicate smaller pores. Based on saturated and centrifugal T2 spectra, a T2 cut-off value can be determined, which distinguishes the movable fluid from the confined fluid. The T2 cut-off value is selected as a point where the T2 spectrum undergoes significant changes during the centrifugation process. The T2 values before this point correspond to the movable fluid, and the T2 values after this point correspond to the confined fluid.
This value is related to the measurement accuracy of irreducible water saturation and permeability [2,3]. At present, there are two methods to determine the T2 cut-off value. On the one hand, the T2 cut-off value under laboratory conditions can be obtained by calculating the T2 spectral integral curve obtained after saturation and the centrifugation of water [4]. However, this method has the disadvantages of a high cost, a long cycle and only limited amounts of data. On the other hand, according to empirical inference, the T2 cut-off values of sandstone and limestone are 33 and 93 ms, respectively; however, the T2 cut-off values of different research areas, target layers and reservoir types vary greatly, resulting in large calculation errors [5].
At present, the more commonly used method to predict the T2 cut-off value is to use the saturated T2 spectrum obtained from the NMR logging curve for conversion. After studying the relationship between the T2 spectrum and the T2 cut-off value of a large number of saturated samples, scholars proposed a method to calculate the variable T2 cut-off value by using the morphological differences in the T2 spectrum [6,7]. On this basis, the normal distribution function is used to fit the T2 spectrum of the centrifugal bound water, and then the variable T2 cut-off value is obtained. Subsequently, to quantitatively characterize the complexity of pore fracture structure, scholars have used fractal theory to characterize the heterogeneity of pore distribution. The results show that there is a linear relationship between the multifractal parameters (D0D10) and the T2 cut-off value. Therefore, a T2 cut-off prediction model based on multifractal parameters is proposed [8,9].
In summary, the relevant literature has established a T2 cut-off prediction model to achieve the quantitative prediction of irreducible water saturation in specific reservoirs [10,11,12]. However, there are still the following problems in this field. Firstly, the correlation results show that there is a relationship between T2 spectral morphology and the T2 cut-off value, but there is a lack of research on representative morphological parameters. At the same time, the relationship between multifractal parameters and T2 cut-off values has been studied, but there are few studies on the correlation between single fractal parameters and T2 cut-off values. Therefore, for specific reservoir types, the applicability of different prediction methods for the T2 cut-off value of the same reservoir needs to be further studied.
Therefore, this paper takes the sandstone samples collected in the Qinshui Basin as an example, and uses low-field nuclear magnetic resonance technology to obtain the saturation and centrifugal T2 spectra of different samples, as well as the T2 cut-off value of each sample. On this basis, from a morphological point of view, the parameters including peak value, peak position and peak area corresponding to the hole are counted, and the relationship between the morphological parameters and the T2 cut-off value is discussed. At the same time, single and multifractal models are used to realize the quantitative characterization of the heterogeneity of the pore fracture distribution of the target sample, and the correlation between the multifractal parameters and the pore fracture parameters is discussed. On this basis, the key parameters are selected from the morphological and fractal parameters, and the T2 cut-off value is calculated using the single fractal dimension, multifractal dimension and spectral morphology discrimination methods, to construct the T2 cut-off value prediction model.

2. Experimental Test and Theoretical Analysis

2.1. Experimental Test

The target samples were collected from the Ordos Basin in the western part of the North China Plate. The region is a typical superimposed basin on the edge of the craton, with a total area of about 37 × 104 km2. The basin is composed of the western margin thrust belt, Tianhuan sag, Yishan slope, western Shanxi bending fold belt, Yimeng uplift and Weibei uplift tectonic belt (Hu et al., 2020) [8]. The Benxi Formation in the Qinshui Basin is mainly composed of alternating marine and continental mudstone, sandy mudstone, claystone, limestone and limonite layers interbedded with coal lines, with a thickness ranging from 10 to 50 m. In this formation, tight sandstone is an important component, but its scale and thickness are relatively small, and the sand bodies are relatively independent of each other. The tight sandstone of the Benxi Formation is mainly composed of fine-grained sandstone, with small and tightly arranged particles and relatively low porosity and permeability; hence, it is called tight sandstone. This sandstone typically has high mechanical strength and stability. Due to the lower porosity and permeability of tight sandstone, its reservoir properties are relatively poor. However, through reservoir modification techniques such as fracturing to increase production, the physical properties of tight sandstone reservoirs can be significantly improved, thereby increasing the efficiency of oil and gas resource extraction. To study the pore distribution of tight sandstone reservoirs at different depths, 11 tight sandstone samples of Benxi Formation were collected from different exploration wells, with each having a sampling depth of 3000~3900 m.
Low-field nuclear magnetic resonance (LF-NMR) was used [13]. Cylinders with a diameter of 25 mm and a length of 30 mm were prepared from sandstone samples. Before the test, the sample was placed in a drying oven and the temperature was adjusted to 105 °C for about 12 h. Then, the dried samples were weighed with a balance. The samples were saturated with distilled water in vacuum (i.e., no confining pressure) for 48 h. After saturation, the wet sample was weighed and the porosity (%) of the sample was calculated. After measuring the T2 spectrum, the sample was centrifuged at a certain speed to obtain the centrifugal T2 spectrum [14,15].

2.2. Model Calculation and Parameter Acquisition

Based on the saturated T2 spectrum, the following three morphological parameters are proposed:
The main peak position (TM) represents the T2 value corresponding to the highest spectral area, in ms.
The percentage of pore volume of smaller pores (SPVP) and larger pores (LPVP) represents the percentage of spectral area in a certain T2 spectrum range, in %; relevant literature has studied the relationship between the T2 value and pore size. However, the conversion coefficients between the two parameters are not consistent. Based on the research results of Yao et al. (2010), the pore fracture system can be divided into small pores (0–102 nm, T2 < 2.5 ms) and large pores (102–104 nm, 2.5 ms < T2) using the T2 spectrum and the response characteristics of pores with different pore sizes [4,16,17,18].
Based on the saturated T2 spectrum, two kinds of fractal parameters are proposed, as follows:
The first type is the single fractal parameter, which is calculated using the single fractal model. The calculation formula is as follows:
lg ( V p ) = ( 3 D w ) lg ( T 2 ) + ( D w 3 ) lg T 2 max
In the formula, Vp is the pore volume percentage of the T2 spectrum under a saturated water condition, in %; T2 max is the maximum transverse relaxation time, in ms; and DS is the fractal dimension in a saturated water state, which is dimensionless.
The relevant literature shows that the fractal curve can be divided into two sections using the T2 cut-off value of saturated water [19,20]. The slope of the curve with a T2 value of less than 2.5 ms corresponds to the heterogeneity of pore distribution with a smaller pore size (D1), and the slope of the curve with a T2 value greater than 2.5 ms corresponds to the heterogeneity of pore distribution with a larger pore size (D2).
The second type is the multifractal parameter, which is calculated using the multifractal model. The specific calculation formula is shown in Zhang et al. (2020) [1]. The multifractal characteristics of NMR data under different test methods (reflecting the inhomogeneity of pore size distribution) are studied using a box counting method [21,22,23]. When analyzing the volume probability of the T2 spectrum in the interval [A, B], it is necessary to determine the scale and measure. The expression of proportion and measure is as follows:
ε = 2 k L
Δ n i ¯ = Δ n i Δ n min i = 1 N ( Δ n i Δ n min )
The T2 value is divided into several equal-length intervals for NMR testing, and the interval scale is represented by this. The analysis range of NMR data is 0~1000 ms, as detailed in Zhang et al. (2020) [1]. The multifractal dimension can be divided into two languages, namely the singular fractal spectrum and the generalized fractal spectrum. The characteristic parameters of f(a) mainly include amin, amax, a0, a0amax, amina0 and a. This is a set of main languages describing multifractal local features, called a multifractal spectrum [24,25]. From the perspective of information theory, another set, q~D(q), is introduced, which is called the generalized fractal dimension [26].

3. Results and Discussion

3.1. NMR T2 Spectrum Distribution and Type Division

Unlike unconventional reservoirs such as coal and shale, tight sandstone reservoirs do not develop micropores, and the T2 spectrum is dominated by single peaks [9]. The T2 spectra of nine sandstone samples can be divided into the following four types. Type A is represented by sample 2, wherein the T2 spectrum shows a bimodal state, and the area of the right T2 spectrum (2.5~100 ms) is larger than that of the left T2 spectrum (T2 < 2.5 ms), indicating that the sample has good pore connectivity and belongs to the macroporous development sample (Figure 1a). Type B is represented by sample 8, wherein the T2 spectrum is unimodal. The area of the right T2 spectrum is significantly larger than that of the left T2 spectrum, indicating that the pore connectivity of this type of sample is poor and that it belongs to the macroporous development sample (Figure 1b). The C-type sample is represented by sample 10. The T2 spectrum shows a single peak state, and the main T2 is distributed in 1~50 ms, indicating that the pore connectivity of this type of sample is very poor, and that it belongs to the mesoporous development type sample (Figure 1c). The D-type sample is represented by sample 52, and the T2 spectrum is unimodal, and the main T2 is distributed in 0.1~2.5 ms, indicating that the pore connectivity of this type of sample is very poor and that it belongs to the small pore development sample (Figure 1d). The results in Table 1 show that the porosity and permeability of type A and B samples are larger than those of type C and D samples. This is because the large pores in the pores and fractures of the former are developed and the pore connectivity is better.

3.2. Correlation between Single Fractal Parameters and T2 cut-off Value

Based on Equation (1), the single fractal dimension values of D1 (small pores), D2 (large pores) and DT (total pores) of the same sample can be calculated using the saturated T2 spectrum of the sample. Figure 2a shows that D1 is not in the range of 2~3, which indicates that the distribution characteristics of small pores have no fractal significance. The D2 and DT of the sandstone samples are 2.77~2.97 and 2.06~2.45, respectively, indicating that the macropores and total pores in the sandstone samples have fractal characteristics. Figure 2b shows that D2 increases with the increase in pore volume percentage, but the correlation is weak, while the correlation between DT and pore volume parameters is weak. Compared with small pores, there is no significant correlation between the volume percentage of large pores and the single fractal dimension value (Figure 2c), indicating that small pores are the key factors restricting the heterogeneity of reservoir pore distribution. Figure 2d shows that D2 increases with the increase in T2 cut-off value, but the correlation is weak. Therefore, it is not feasible to predict the T2 cut-off value using the single fractal dimension parameter.

3.3. Correlation between Multifractal Parameters and T2 cut-off Value

As described in Section 2.2, the multifractal dimension values of all sandstone samples can be calculated using the saturated T2 spectrum curve (Figure 3). The results show that the q~D(q) spectra of all sandstone samples are obvious inverse S-typed, which is also another typical feature of the sample pore size distribution conforming to multifractal characteristics. The spectral line can effectively characterize the complexity of the pore size distribution characteristics at different stages, to reveal the local difference characteristics of the mosaic in the whole. The results show that the D0 of all samples is one, indicating that each partition box contains pores, corresponding to the Euclidean dimension of one-dimensional distribution. D1 is the information entropy dimension, which is related to the entropy value of the system and can be used to characterize the heterogeneity of the overall pore distribution. The higher the value, the higher the pore distribution range and the stronger the heterogeneity. Figure 4 shows that there is a linear positive correlation between the singular fractal spectrum parameters and the generalized fractal spectrum parameters. Except for the correlation coefficient between D0D10 and a0a10, which is 0.61, the other correlation coefficients are greater than 0.9. Therefore, only the generalized fractal parameters are used to characterize the multifractal characteristics.
According to the above, D0D10 characterizes the heterogeneity of pore volume development area distribution, and D−10D0 characterizes the heterogeneity of pore volume underdevelopment area distribution. Taking Figure 1d as an example, when T2 is less than 2.5 ms, the pore volume ratio of small pores is greater than 80%; therefore, this part of the region can be characterized by D0D10. On the contrary, when T2 is greater than 2.5 ms, the corresponding macropore volume ratio is less than 20%; therefore, this part of the region can be characterized by D−10D0. It can be seen from Figure 5a that the correlation between D0D10 and D−10D0 is weak, indicating that the distribution heterogeneity of the low value area and the high value area of pore volume in this kind of sample is strong. At the same time, Figure 5b,c show that, with the increase in D−10D0, D−10D10 increases linearly, but there is no obvious linear correlation between D0D10 and D−10D10, which indicates that the low value area of pore volume in this kind of sample controls the overall heterogeneity of pore distribution.
Figure 6a shows that there is no correlation between the cut-off values of D10 and T2. Figure 6b shows that there is a linear positive correlation between the cut-off values of D−10 and T2, wherein the R2 value is 0.75. Figure 6c,d show that there is a linear positive correlation between T2 and D−10D10 and D−10/D10 cut-off values, with R2 being 0.68 and 0.76, respectively. Figure 6e shows that there is a linear positive correlation between D−10D0 and T2 cut-off value, wherein the R2 value is 0.75, which indicates that the heterogeneity parameter of the pore volume distribution of regional pore size can be used to predict the T2 cut-off value.

3.4. Correlation between T2 Spectral Morphological Parameters and T2 cut-off Values

Peak position, peak value and T2 spectral area ratio are three important parameters used to distinguish T2 spectral types. Figure 7a shows that the correlation between TM and the volume percentage of small pores and large pores is weak. Figure 7b shows that the correlation between TM and the single fractal dimension D1 is weak, but the linear fitting degree with D2 is high, and the R2 value can reach 0.45. Figure 7c shows that the correlation between TM and the multifractal dimension is very poor. The above results show that the physical meaning of the morphological parameters is completely different from those of the fractal parameters.
In addition, Figure 8 shows that, except for the single fractal parameter D2, the pore volume percentage has a weak correlation with other single and multifractal parameters. This also shows that the pore volume percentage is different from the fractal parameters and can be used as one of the parameters for T2 cut-off value prediction.
When the T2 cut-off value of the sample is 0~10 ms, the T2 cut-off value decreases with the increase in the TM value (Figure 8a), but its linear correlation is much weaker than with the multifractal parameters (Figure 6c,d). In addition, there is no significant correlation between the volume percentage of small and large pores and the T2 cut-off value (Figure 9a–c). It is worth noting that this result is not consistent with the results of Zhang et al. (2023), which were mainly affected by the number of samples and the type of reservoir. As a result, the correlation between the morphological parameters and the T2 cut-off value is weak [9].

3.5. T2 cut-off Value Prediction Model Based on Multifractal Parameters

In the above content, the parameters affecting the T2 cut-off value are found by using the single fractal and multifractal methods and the morphological discrimination method. The results show that the relevant parameters affecting the T2 cut-off value include D−10D10, D−10/D10, D−10D0, TM and D2. Therefore, taking D−10D10, D−10/D10, D−10D0, TM and D2 as the prediction parameters for the T2 cut-off value of sandstone samples, a multiple linear regression analysis is carried out using SPSS software v19.0 (SPSS Inc., Chicago, IL, USA), and the prediction model for the T2 cut-off value is obtained. The calculation formula is as follows:
T2 cut-off value =7.386 − 3.30 × (D−10D10) − 0.29 × D−10/D10 + 5.40 × (D−10D0) − 0.24 × TM − 1.14 × D2
The results show that the positive correlation coefficient of the model can reach 0.88. The model is used to predict the T2 cut-off value of the samples in the study area. The basic information of the samples to be tested is shown in Table 2 and Figure 10. The comparison shows that the error range between the T2 cut-off value calculated using Equation (4) and the T2 cut-off value measured using centrifugal experiment NMR is less than 5% (Figure 10), indicating that the T2 cut-off values calculated by the above two methods are highly consistent. Therefore, it can be considered that the model established in this paper has good applicability.

4. Conclusions

In this paper, the parameters including peak value, peak position and peak area corresponding to small holes are counted, and the relationship between the morphological parameters and the T2 cut-off value is discussed. At the same time, the single and multifractal models are used to realize the quantitative characterization of the heterogeneity of the pore fracture distribution of the target sample, and the correlation between the multifractal parameters and the pore fracture parameters is discussed. On this basis, the key parameters are selected from the morphological and fractal parameters, and the T2 cut-off value is calculated using the single fractal dimension, multifractal dimension and spectral morphology discrimination methods, in order to construct the T2 cut-off value prediction model. The following conclusions are drawn.
(1)
The T2 spectra of sandstone samples in the study area can be divided into four types. The A-type T2 spectrum shows a bimodal state, and the area of the right T2 spectrum is larger than that of the left T2 spectrum, indicating that the samples have good pore connectivity and belong to the macroporous development type samples. The B-type T2 spectrum is unimodal, and the pore connectivity of the samples is poor, indicating that they belong to the macroporous development sample. The T2 spectrum of the C-type samples is unimodal, and the pore connectivity is very poor, indicating that they belong to the mesoporous development sample. The T2 spectrum of the D-type samples is unimodal, and the main T2 is distributed within 0.1~2.5 ms, and the pore connectivity is very poor, indicating that they belong to the small pore development sample.
(2)
The single model shows that small pore volume is the key factor restricting the non-uniformity of reservoir pore distribution. Compared with other single fractal parameters, D2 increases with the increase in the T2 cut-off value, but the correlation is weak. Therefore, it is not feasible to predict the T2 cut-off value using the single fractal dimension parameter.
(3)
The multifractal model shows that D−10D10 increases linearly with the increase in D−10D0, but there is no obvious linear correlation between D0D10 and D−10D10, indicating that the low pore volume area in this kind of sample controls the overall heterogeneity of pore distribution.
(4)
The relevant parameters affecting the cut-off value of T2 include D−10D10, D−10/D10, D−10D0, TM and D2. Therefore, the T2 cut-off value prediction model is constructed based on the above five parameters. The T2 cut-off value calculated by the model is highly consistent with the experimental value, which proves the reliability of the model.

Author Contributions

Methodology, D.Z.; software, T.T.; validation, Y.S. and J.Z.; formal analysis, Y.H.; writing—review and editing, Z.Q. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

The original contributions presented in the study are included in the article, further inquiries can be directed to the corresponding author.

Conflicts of Interest

Authors Di Zhang, Tian Tian and Yong Shi were employed by the Second Exploration Team of Shandong Coalfield Geological Bureau. The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest. The Second Exploration Team of Shandong Coalfield Geological Bureau had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. Zhang, J.; Wei, C.; Ju, W.; Qin, Z.; Ji, Y.; Quan, F.; Hu, Y. Microscopic distribution and dynamic variation of water under stress in middle and high rank coal samples. J. Nat. Gas Sci. Eng. 2020, 79, 103369. [Google Scholar] [CrossRef]
  2. Bai, Y.; Wang, S.; Xu, B.; Chen, L.; Li, Y.; Dong, Z. Experimental study on the relationship between irreducible water saturation and micro pore throat in tight sandstone. China Offshore Oil Gas 2022, 34, 65–71. [Google Scholar]
  3. Tan, L.; Cai, M.; Cai, D.; Zhang, C. Study on Bound Water Saturation of Tight Sandstone Based on Nuclear Magnetic Experiment. Contemp. Chem. Ind. 2022, 51, 1066–1070. [Google Scholar]
  4. Yao, Y.; Liu, D.; Che, Y.; Tang, D.; Tang, S.; Huang, W. Petrophysical characterization of coals by low-field nuclear magnetic resonance (NMR). Fuel 2010, 89, 1371–1380. [Google Scholar] [CrossRef]
  5. Zhu, L.; Zhang, C.; Shi, W.; Wei, Y.; Chen, Y.; Liu, J.; Gou, C.; Cheng, Y. Study on the Method of Prediction of Irreducible Water Saturation by Combining Mercury Intrusion and NMR Logging Data. Sci. Technol. Eng. 2016, 16, 22–29. [Google Scholar]
  6. Wang, Z.; Zhang, C.; Xiao, C.; Chen, X.; Song, F. Experimental study on T2 cut-off value of low permeability reservoir and its logging application. Geophys. Prospect. Pet. 2004, 19, 652–655. [Google Scholar]
  7. Shao, W.; Ding, Y.; Xiao, F. On the Method of Determing T2 Cutoff Value with the T2 Spectrum Characteristics. Well Logging Technol. 2009, 33, 430–435. [Google Scholar]
  8. Hu, Y.; Guo, Y.; Zhang, J.; Shangguan, J.; Li, M.; Quan, F.; Li, G. A method to determine nuclear magnetic resonance T2 cut-off value of tight sandstone reservoir based on multifractal analysis. Energy Sci. Eng. 2020, 8, 1135–1148. [Google Scholar] [CrossRef]
  9. Junjian, Z.; Guo, Y.; Qin, Z.; Chongtao, W.; Qinhong, H.; Vandeginste, V.; Heyao, M.; Peng, Y.; Pengfei, Z. Predicting irreducible water saturation of unconventional reservoirs by using NMR T2 spectra: Methods of morphological division and fractal models. Acta Geol. Sinaca 2023. [Google Scholar] [CrossRef]
  10. Lai, F.; Li, Z.; Zhang, W.; Dong, H.; Kong, F.; Jiang, Z. Investigation of Pore Characteristics and Irreducible Water Saturation of Tight Reservoir Using Experimental and Theoretical Methods. Energy Fuels 2018, 32, 3368–3379. [Google Scholar] [CrossRef]
  11. Meng, M.; Ge, H.; Shen, Y.; Ji, W.; Ji, W.; Wang, Q. Rock Fabric of Tight Sandstone and Its Influence on Irreducible Water Saturation in Eastern Ordos Basin. Energy Fuels 2023, 37, 3685–3696. [Google Scholar] [CrossRef]
  12. Lai, F.; Li, Z.; Dong, H.; Jiang, Z.; Mao, G.T. Micropore structure characteristics and water distribution in a coalbed methane reservoir. Aust. J. Earth Sci. 2019, 66, 741–750. [Google Scholar] [CrossRef]
  13. Sun, Z.; Li, M.; Yuan, S.; Hou, X.; Bai, H.; Zhou, F.; Liu, X.; Yang, M. The flooding mechanism and oil recovery of nanoemulsion on the fractured/non-fractured tight sandstone based on online LF-NMR experiments. Energy 2024, 291, 130226. [Google Scholar] [CrossRef]
  14. Zhang, T.; Tang, M.; Yuan, L.; Zhu, G.; Zhang, C.; Liu, Y.; Li, Y.; Wang, W.; Yang, X. Experimental study on stress-dependent multiphase flow in ultra-low permeability sandstone during CO2 flooding based on LF-NMR. Energy 2023, 278, 127874. [Google Scholar] [CrossRef]
  15. Hou, X.; Zhu, Y.; Chen, S.; Wang, Y.; Liu, Y. Investigation on pore structure and multifractal of tight sandstone reservoirs in coal bearing strata using LF-NMR measurements. J. Pet. Sci. Eng. 2020, 187, 106757. [Google Scholar] [CrossRef]
  16. Dai, Q.; He, Y.; Da, Z.; Zhu, Y.; Wang, R.; Yu, H.; Zhou, Y.; Wang, G.; Li, Y.; Wu, H. Evaluating Pore Structure of Shale Reservoirs through Pore Space Multifractal Characterization. Geofluids 2023, 2023, 3042994. [Google Scholar] [CrossRef]
  17. Sun, Y.; Zhai, C.; Zhao, Y.; Xu, J.; Cong, Y.; Zheng, Y.; Tang, W. Multifractal Analysis and Neural Network Prediction of Pore Structures in Coal Reservoirs Based on NMR T2 Spectra. Energy Fuels 2021, 35, 11306–11318. [Google Scholar] [CrossRef]
  18. Zheng, S.; Yao, Y.; Liu, D.; Cai, Y.; Liu, Y.; Li, X. Nuclear magnetic resonance T2 cut-offs of coals: A novel method by multifractal analysis theory. Fuel 2019, 241, 715–724. [Google Scholar] [CrossRef]
  19. Zheng, S.; Sang, S.; Yao, Y.; Liu, D.; Liu, S.; Wang, M.; Feng, G. A multifractal-based method for determination NMR dual T2 cut-offs in coals. J. Pet. Sci. Eng. 2022, 214, 110488. [Google Scholar] [CrossRef]
  20. Ge, X.; Fan, Y.; Zhu, X.; Chen, Y.; Li, R. Determination of nuclear magnetic resonance T2 cut-off value based on multifractal theory: An application in sandstone with complex pore structure. Geophysics 2015, 80, D11–D21. [Google Scholar] [CrossRef]
  21. Liu, H.; Xie, H.; Wu, F.; Li, C.; Gao, R. A novel box-counting method for quantitative fractal analysis of three-dimensional pore characteristics in sandstone. Int. J. Min. Sci. Technol. 2024, 34, 479–489. [Google Scholar] [CrossRef]
  22. Wu, M.; Wang, W.; Shi, D.; Song, Z.; Li, M.; Luo, Y. Improved box-counting methods to directly estimate the fractal dimension of a rough surface. Meas. J. Int. Meas. Confed. 2021, 177, 109303. [Google Scholar] [CrossRef]
  23. Ficker, T. Normalized multifractal spectra within the box-counting method. Czechoslov. J. Phys. 2000, 50, 389–403. [Google Scholar] [CrossRef]
  24. Xie, Y.; Michael, Z.; Li, C. Anisotropic characteristics of acoustic emission and the corresponding multifractal spectrum during progressive failure of shale under cyclic loading. Int. J. Rock Mech. Min. Sci. 2023, 165, 105364. [Google Scholar] [CrossRef]
  25. Liu, M.; Xie, R.; Guo, J.; Jin, G. Characterization of Pore Structures of Tight Sandstone Reservoirs by Multifractal Analysis of the NMR T2 Distribution. Energy Fuels 2018, 32, 12218–12230. [Google Scholar] [CrossRef]
  26. Fang, S.; Liang, Y.; Liu, Z.; Yang, W. Heterogeneity of tight sandstone reservoirs based on fractal theory: The Xu-6 member of Xujiahe Formation in Guang’an area, central Sichuan Basin. Arab. J. Geosci. 2021, 14, 1515. [Google Scholar]
Figure 1. Typical T2 spectral distribution of sandstone samples: (a) Typical T2 spectral distribution of type A; (b) Typical T2 spectral distribution of type B; (c) Typical T2 spectral distribution of type C; (d) Typical T2 spectral distribution of type D.
Figure 1. Typical T2 spectral distribution of sandstone samples: (a) Typical T2 spectral distribution of type A; (b) Typical T2 spectral distribution of type B; (c) Typical T2 spectral distribution of type C; (d) Typical T2 spectral distribution of type D.
Processes 12 01811 g001
Figure 2. Correlation of univariate fractal dimension values with pore volume parameters and T2 cut-off: (a) fractal dimension curves; (b,c) correlation analyses of fractal dimension values with pore volume percentages; (d) correlation analysis of fractal dimension values with T2 cut-off.
Figure 2. Correlation of univariate fractal dimension values with pore volume parameters and T2 cut-off: (a) fractal dimension curves; (b,c) correlation analyses of fractal dimension values with pore volume percentages; (d) correlation analysis of fractal dimension values with T2 cut-off.
Processes 12 01811 g002
Figure 3. Pore structure characterization of samples based on multi-fractal models: (a) multiple-fractal model of DS; (b) multiple-fractal model of D.
Figure 3. Pore structure characterization of samples based on multi-fractal models: (a) multiple-fractal model of DS; (b) multiple-fractal model of D.
Processes 12 01811 g003
Figure 4. Correlation analysis of singular fractal parameters with generalized fractal parameters: (a) D−10D0 and a−10a0; (b) D0D10 and a0a10; (c) D−10D10 and a−10a10; (d) D10 and a10.
Figure 4. Correlation analysis of singular fractal parameters with generalized fractal parameters: (a) D−10D0 and a−10a0; (b) D0D10 and a0a10; (c) D−10D10 and a−10a10; (d) D10 and a10.
Processes 12 01811 g004
Figure 5. Generalized fractal parameter correlation analysis: (a) D−10D0 and D0D10; (b) D−10D0 and D−10D10; (c) D0D10 and D−10D10.
Figure 5. Generalized fractal parameter correlation analysis: (a) D−10D0 and D0D10; (b) D−10D0 and D−10D10; (c) D0D10 and D−10D10.
Processes 12 01811 g005
Figure 6. Correlation analysis of multiple fractal dimensions with T2 cut-offs: (a) Correlation analysis of D10 with T2 cut-offs; (b) Correlation analysis of D−10 with T2 cut-offs; (c) Correlation analysis of D−10D10 with T2 cut-offs; (d) Correlation analysis of D−10/D10 with T2 cut-offs; (e) Correlation analysis of D−10D0 with T2 cut-offs; (f) Correlation analysis of D0D10 with T2 cut-offs.
Figure 6. Correlation analysis of multiple fractal dimensions with T2 cut-offs: (a) Correlation analysis of D10 with T2 cut-offs; (b) Correlation analysis of D−10 with T2 cut-offs; (c) Correlation analysis of D−10D10 with T2 cut-offs; (d) Correlation analysis of D−10/D10 with T2 cut-offs; (e) Correlation analysis of D−10D0 with T2 cut-offs; (f) Correlation analysis of D0D10 with T2 cut-offs.
Processes 12 01811 g006aProcesses 12 01811 g006b
Figure 7. TM correlation analysis with pore structure parameters and fractional dimension values: (a) Correlation analysis of Pore volume percentage with TM; (b) Correlation analysis of Single fractal dimension D1 with TM; (c) Correlation analysis of D−10/D10 with TM.
Figure 7. TM correlation analysis with pore structure parameters and fractional dimension values: (a) Correlation analysis of Pore volume percentage with TM; (b) Correlation analysis of Single fractal dimension D1 with TM; (c) Correlation analysis of D−10/D10 with TM.
Processes 12 01811 g007
Figure 8. Correlation of pore volume percentage of small pores with multiple and single weight fractal dimensions, respectively: (a) Correlation of pore volume percentage with D−10/D10 and D−10D10; (b) Correlation of pore volume percentage with D2 and DT.
Figure 8. Correlation of pore volume percentage of small pores with multiple and single weight fractal dimensions, respectively: (a) Correlation of pore volume percentage with D−10/D10 and D−10D10; (b) Correlation of pore volume percentage with D2 and DT.
Processes 12 01811 g008
Figure 9. Correlation analysis of T2 cut-off values with morphological parameters: (a) Correlation analysis of T2 cut-off values with peak position; (b) Correlation analysis of T2 cut-off values with pore volume percentage of small pore volume; (c) Correlation analysis of T2 cut-off values with pore volume percentage of large pore volume.
Figure 9. Correlation analysis of T2 cut-off values with morphological parameters: (a) Correlation analysis of T2 cut-off values with peak position; (b) Correlation analysis of T2 cut-off values with pore volume percentage of small pore volume; (c) Correlation analysis of T2 cut-off values with pore volume percentage of large pore volume.
Processes 12 01811 g009
Figure 10. Calculation process of T2 cut-off values of samples to be tested and similarity test.
Figure 10. Calculation process of T2 cut-off values of samples to be tested and similarity test.
Processes 12 01811 g010
Table 1. Sample base information sheet.
Table 1. Sample base information sheet.
Sample NumberΦwater (%)Φnitrogen (%)Permeability
(mD)
T2 cut-off (ms)Saturated Water Swi (%)<2.5 (ms)2.5~100 (ms)Main Peak Position (ms)
23.756.551.01732.9842.910.410.5114.97
83.253.500.11694.8227.960.180.7421.09
103.243.200.10077.7971.190.330.624.78
115.745.410.38718.3964.360.250.686.00
122.622.800.11164.6463.280.460.503.02
163.553.200.03526.6763.150.310.666.00
4510.249.273.03581.3011.490.130.8618.82
522.222.010.01460.4966.380.920.0714.97
589.779.030.57305.6551.930.270.727.54
Table 2. Test model sample information sheet.
Table 2. Test model sample information sheet.
SampleD−10D10D−10/D10D−10D0TMD2T2C-NMR
12.545.072.171.362.951.60
21.743.301.500.392.991.60
35.0312.614.464.262.982.68
44.049.453.523.802.992.03
54.438.784.002.412.993.10
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Zhang, D.; Tian, T.; Shi, Y.; He, Y.; Zhang, J.; Qin, Z. Method for Predicting Bound Water Saturation in Tight Sandstone Reservoirs Using Morphology and Fractal Models. Processes 2024, 12, 1811. https://doi.org/10.3390/pr12091811

AMA Style

Zhang D, Tian T, Shi Y, He Y, Zhang J, Qin Z. Method for Predicting Bound Water Saturation in Tight Sandstone Reservoirs Using Morphology and Fractal Models. Processes. 2024; 12(9):1811. https://doi.org/10.3390/pr12091811

Chicago/Turabian Style

Zhang, Di, Tian Tian, Yong Shi, Yaomiao He, Junjian Zhang, and Zhenyuan Qin. 2024. "Method for Predicting Bound Water Saturation in Tight Sandstone Reservoirs Using Morphology and Fractal Models" Processes 12, no. 9: 1811. https://doi.org/10.3390/pr12091811

APA Style

Zhang, D., Tian, T., Shi, Y., He, Y., Zhang, J., & Qin, Z. (2024). Method for Predicting Bound Water Saturation in Tight Sandstone Reservoirs Using Morphology and Fractal Models. Processes, 12(9), 1811. https://doi.org/10.3390/pr12091811

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop