
==== Front
Sci Rep
Sci Rep
Scientific Reports
2045-2322
Nature Publishing Group UK London

39237621
71351
10.1038/s41598-024-71351-0
Article
Quantitative classification evaluation model for tight sandstone reservoirs based on machine learning
Song Xinglei 12
Feng Congjun fengcj@nwu.edu.cn

12
Li Teng 345
Zhang Qin 6
Pan Xinhui 12
Sun Mengsi 7
Ge Yanlong 12
1 grid.412262.1 0000 0004 1761 5538 State Key Laboratory of Continental Dynamics, Northwest University, Xi’an, 710069 China
2 https://ror.org/00z3td547 grid.412262.1 0000 0004 1761 5538 Department of Geology, Northwest University, No. 229, Taibai North Road, Xi’an, 710069 Shaanxi China
3 https://ror.org/040c7js64 grid.440727.2 0000 0001 0608 387X School of Petroleum Engineering, Xi’an Shiyou University, Xi’an, 710065 China
4 https://ror.org/03m01yf64 grid.454828.7 0000 0004 0638 8050 Engineering Research Center of Development and Management for Low to Ultra-Low Permeability Oil & Gas Reservoirs in West China, Ministry of Education, Xi’an, 710065 China
5 Xi’an Key Laboratory of Tight Oil (Shale Oil) Development, Xi’an, 710065 China
6 https://ror.org/02awe6g05 grid.464414.7 0000 0004 1765 2021 PetroChina Research Institute of Petroleum Exploration & Development, Beijing, 100083 People’s Republic of China
7 https://ror.org/01dyr7034 grid.440747.4 0000 0001 0473 0092 School of Petroleum Engineering and Environmental Engineering, Yan’an University, Yan’an, 716000 China
5 9 2024
5 9 2024
2024
14 2071222 4 2024
27 8 2024
© The Author(s) 2024
2024
https://creativecommons.org/licenses/by-nc-nd/4.0/ Open Access This article is licensed under a Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International License, which permits any non-commercial use, sharing, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if you modified the licensed material. You do not have permission under this licence to share adapted material derived from this article or parts of it. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by-nc-nd/4.0/.
Tight sandstone reservoirs are a primary focus of research on the geological exploration of petroleum. However, many reservoir classification criteria are of limited applicability due to the inherent strong heterogeneity and complex micropore structure of tight sandstone reservoirs. This investigation focused on the Chang 8 tight reservoir situated in the Jiyuan region of the Ordos Basin. High-pressure mercury intrusion experiments, casting thin sections, and scanning electron microscopy experiments were conducted. Image recognition technology was used to extract the pore shape parameters of each sample. Based on the above, through grey relational analysis (GRA), analytic hierarchy process (AHP), entropy weight method (EWM) and comprehensive weight method, the relationship index Q1 between initial productivity and high pressure mercury injection parameters and the relationship index Q2 between initial productivity and pore shape parameters are obtained by fitting. Then a dual-coupled comprehensive quantitative classification prediction model for tight sandstone reservoirs was developed based on pore structure and shape parameters. A quantitative classification study was conducted on the target reservoir, analyzing the correlation between reservoir quality and pore structure and shape parameters, leading to the proposal of favourable exploration areas. The research results showed that when Q1 ≥ 0.5 and Q2 ≥ 0.5, the reservoir was classified as type I. When Q1 > 0.7 and Q2 > 0.57, it was classified as type I1, indicating a high-yield reservoir. When 0.32 < Q1 < 0.47 and 0.44 < Q2 < 0.56, was classified as type II. When 0.1 < Q1 < 0.32 and 0.3 < Q2 < 0.44, it was classified as type III. Type I reservoirs exhibit a zigzag pattern in the northwest part of the study area. Thus, the northwest should be prioritized in actual exploration and development. Additionally, the initial productivity of tight sandstone reservoirs showed a positive correlation with the porosity, permeability, sorting coefficient, coefficient of variation, and median radius. Conversely, it demonstrated a negative correlation with the median pressure and displacement pressure. The perimeters of pores, their circularity, and the length of the major axis showed a positive correlation with the porosity, permeability, sorting coefficient, coefficient of variation, and median radius. On the other hand, they exhibited a negative correlation with the median pressure and displacement pressure. This study quantitatively constructed a new classification and evaluation system for tight sandstone reservoirs from the perspective of microscopic pore structure, achieving an overall model accuracy of 93.3%. This model effectively predicts and evaluates tight sandstone reservoirs. It provides new guidance for identifying favorable areas in the study region and other tight sandstone reservoirs.

Keywords

Tight sandstone
Pore structure
Quantitative evaluation
High-pressure mercury injection
Image recognition
Machine learning
Subject terms

Geology
Geophysics
Petrology
http://dx.doi.org/10.13039/501100017596 Natural Science Basic Research Program of Shaanxi Province Grant No. 2020JQ-798 Grant No. 2017JM4013 Feng Congjun Sun Mengsi issue-copyright-statement© Springer Nature Limited 2024
==== Body
pmcIntroduction

With the depletion of conventional oil and gas reservoirs, tight oil reservoirs have gradually become a hot topic and a focal point for exploration and development, both domestically and internationally1–9. However, tight sandstone oil reservoirs exhibit complex reservoir characteristics, primarily manifested in their deep burial depths, wide distribution, and complex depositional processes. The reservoirs exhibit characteristics of low porosity, poor permeability, and high heterogeneity. The dominant pores are micro- and nano-scale, with narrow and dispersed throats, and are unfavorable for the migration and accumulation of oil and gas10–15. These factors necessitate considering the interdependent influences of multiple factors when classifying and evaluating tight sandstone reservoirs, which affects the accuracy of reservoir evaluation and hinders the selection of high-quality reservoirs. Therefore, the rapid and effective classification and evaluation of tight sandstone reservoirs has long been a focal point of scholarly research.

The quality of the reservoir is a key factor that determines the oil and gas production capacity. The classification and evaluation of reservoirs are central to reservoir studies and play a significant role in oilfield development. With the continuous advancement of oilfield development technologies, reservoir classification and evaluation methods have become increasingly diverse, gradually evolving from qualitative to quantitative research and from macro-parameter to micro-parameter evaluation. At present, both domestic and international scholars classify reservoirs using two main methods. The first is the traditional classification and evaluation method, which directly uses indicators such as the lithology, physical properties, pore structure, sedimentary facies, and oil and production experiments for classification. For example, Wei et al. classified the tight sandstone reservoirs of the Sha Creek Formation in the central Sichuan Basin based on the transverse relaxation (T2) distribution of nuclear magnetic resonance16. Xu et al. studied the characteristics and controlling factors of tight sandstone using thin-section casting, scanning electron microscopy, X-ray diffraction (XRD), and spontaneous imbibition experiments17. Wu et al. analyzed the logging response characteristics using core data and electric imaging logging data and identified the reservoir type with the highest industrial production in the study area18. Zhang et al. established classification criteria for the third member of the Quan Formation based on mercury injection curves, core physical properties, and sedimentary facies characteristics19. Talib et al. quantitatively characterized tight oil and gas reservoirs through rock physics experiments and seismic inversion profiles20.

The second approach to reservoir classification involves initially choosing evaluation parameters that align with the geological conditions of the target area. Subsequently, machine learning techniques such as GRA the AHP, the EWM, and fuzzy analysis are employed to assign weight coefficients to each evaluation parameter. Finally, the reservoir is comprehensively scored. For example, Fang et al. proposed an automatic classification and verification method for reservoir types based on k-means clustering and Bayesian discriminant theory, using core logging and logging data from coring wells, combined with physical characteristics such as reservoir deposition and diagenesis21. Li et al. classified the Fuyu reservoir using GRA, Q clustering analysis, and discriminant analysis22. Wang et al.combined AHP and EWM, used the multi-factor superposition method, and established a new reservoir classification and evaluation method23. Fan et al. quantified the weight of evaluation parameters’ contribution to production by combining the relationships between variables and directional good production using the GRA24. Niu et al. proposed a new machine learning framework (GCA-CE-MGPK) for shale reservoirs, achieving efficient and accurate multi-scale evaluation of shale reservoirs25. In summary, traditional classification and evaluation methods are costly, inefficient and require extensive experimental data. They are mainly suitable for specific regions, making them inadequate for large-scale reservoir evaluation and prediction. Although machine learning techniques can improve efficiency and reduce costs, their accuracy often depends on the optimization of various mathematical methods, leading to high subjectivity in some models and lower overall precision, failing to meet the practical needs of production. Moreover, previous studies have primarily focused on evaluating single factors, lacking the integration of macro and micro perspectives. Based on these, this study combined multiple machine learning methods to directly link actual oilfield production data with micro-scale pore shape and structure parameters, effectively integrating macro and micro parameters.

Given the significant influence of subjective factors on the classification criteria for the quantitative evaluation of conventional reservoirs, adopting a new method for reservoir evaluation is essential. This study focuses on the Chang 8 tight sandstone reservoir in the Jiyuan area of the Ordos Basin, extracting pore shape parameters from 52 rock samples. Combined with the experimental data of high pressure mercury injection and the actual initial production capacity of the oil field. Through GRA, AHP, EWM and comprehensive weight method, the relationship index Q1 between initial productivity and high pressure mercury injection parameters and the relationship index Q2 between initial productivity and pore shape parameters are obtained by fitting. Then a dual-coupled comprehensive quantitative classification prediction model for tight sandstone reservoirs was developed based on pore structure and shape parameters. A quantitative classification study was conducted on the target reservoir, analyzing the correlation between reservoir quality and pore structure and shape parameters, leading to the proposal of favourable exploration areas. This method effectively combined the subjectivity-influenced AHP with the objectivity-influenced EWM to calculate the comprehensive weight coefficient, mitigating the impact of subjective factors and enhancing the model's accuracy. Validation results indicate that the model has an overall accuracy of 93.3%. Therefore, it was an effective tool for predicting and classifying tight sandstone reservoirs. It is significant for further exploration in the study area and other similar reservoirs.

Geological setting

The Ordos Basin is a large, multi-cycle, cratonic basin that formed on the crystalline basement during the Paleoproterozoic–Mesoproterozoic. The Ordos Basin, the second-largest sedimentary basin in China, has experienced five significant stages of sedimentary evolution. These stages include the middle to late Proterozoic rift valley, the early Paleozoic shallow marine platform, the late Paleozoic nearshore plain, the Mesozoic inland lake basin, and Cenozoic peripheral subsidence. This basin is known for its substantial reserves of oil and gas. The Ordos Basin extends across five provinces and regions, namely, Shaanxi, Gansu, Shanxi, Ningxia, and Inner Mongolia. Geographically, it stretches from the Yin Mountains in the north to the Qinling Mountains in the south, and from the Liupan Mountains in the west to the Lvliang Mountains in the east. The basin’s total area is 25 × 104 km2, with favorable areas covering 9.9 × 104 km2. The estimated resource volume is 6.2 × 1012 m3, indicating significant exploration and development potential. Based on the basin’s geological nature, tectonic evolution, and structural pattern, the Ordos Basin can be divided into six primary tectonic units: the northern Shaanxi slope, the Tianhuan Depression, the western thrust fault zone, the Yimeng Uplift, the Weihebei Uplift, and the western Shanxi fold belt. The Jiyuan area, located in the central-western part of the Ordos Basin, covers a total area of 1302 km2 (Fig. 1a, c). This area spans the two primary tectonic units of the northern Shaanxi slope and the Tianhuan Depression, exhibiting a gently inclined monocline structure towards the west. Since the Mesozoic, the basin has developed thick fluvio-lacustrine deposits. In the Cenozoic, rift valleys were formed around the basin due to fault subsidence. The overall geological conditions are relatively complex, posing challenges for exploration. However, the area is rich in oil and gas resources, indicating favourable exploration prospects26–29. The proven petroleum geological reserves in this area amount to 800 × 106 t, with annual crude oil production of 700 × 104 t, making it the oilfield with the largest reserves and production levels in the Ordos Basin from the Mesozoic. Existing exploration results indicate that the Chang 8 oil-bearing formation is one of the most favourable hydrocarbon accumulation zones in the Jiyuan area, with a proven favourable oil-bearing area of 1500 km2.Fig. 1 (a) Location of the study area(modified from Tong29), (b) columnar diagram of the Chang 8 formation, (c) well location distribution map of the study area.

The Chang 8 reservoir is located in the lower part of the Upper Triassic Yan’an Formation. It is primarily composed of grey sandstone and dark black mudstone interbeds. These sedimentary microfacies are predominantly characterized by subaqueous distributary channels and underwater distributary bays, indicating a deposition pattern typical of a shallow-water deltaic environment (Fig. 1b). Based on the thin-section identification of the study area (Fig. 2), the lithology of the Chang 8 reservoir is predominantly composed of fine-grained feldspathic sandstone, feldspathic lithic sandstone, and a small amount of feldspar sandstone. The detrital components in the study area mainly consist of quartz, feldspar, and detritus. The ranges of contents are as follows: the quartz content is 20.1% to 58.6%, with an average of 31.21%; the feldspar content is 23.56% to 57.62%, with an average of 34.43%; and the detritus content is 6.25% to 29.45%, with an average of 21.38%.Fig. 2 Triangular diagram and detrital composition diagram of the study area. (a) Triangular classification diagram of the sandstone in the Chang 8 reservoir, (b) histogram of the relative content of detrital components in the Chang 8 reservoir.

Materials and statistical methods

Materials and experiments

Materials

In this study, 52 drilling core samples were obtained from the Chang 8 reservoir in Jiyuan, Ordos Basin, with all samples exhibiting a fine sandstone lithology. The samples underwent oil washing, gas permeability measurements, and the weight method for porosity calculation, allowing the determination of the reservoir’s petrophysical parameters (Table 1). The samples' average porosity was 8.23%, between 2.41 and 13.6%. The average permeability was 0.18 × 10–3 µm2, ranging between 0.01 × 10–3 µm2 and 1.10 × 10–3 µm2. Subsequently, thin-section casting and scanning electron microscopy experiments were conducted, resulting in 300 photographs. Additionally, high-pressure mercury intrusion was performed on the 52 samples to obtain the micropore throat characteristic parameters. Table 1 Petrophysical parameters of the experimental samples.

Well	Core number	Diameter (cm)	Length (cm)	Porosity (%)	Permeability (10–3 µm2)	Lithology	
JY-1	1	2.50	5.20	9.29	0.39	FLFS	
JY-2	2	2.48	5.20	9.32	0.14	FLFS	
JY-3	3	2.50	5.10	10.36	0.94	FLFS	
JY-4	4	2.49	5.00	10.43	0.19	FLFS	
JY-5	5	2.50	5.20	11.54	0.23	FLFS	
JY-6	6	2.50	5.20	9.08	0.17	FLFS	
JY-7	7	2.50	5.20	9.97	0.13	FLFS	
JY-8	8	2.48	5.10	8.53	0.12	FLFS	
JY-9	9	2.50	5.00	9.09	0.10	FLFS	
JY-10	10	2.50	5.20	9.38	0.10	FLFS	
JY-11	11	2.50	5.20	12.25	0.26	FLFS	
JY-12	12	2.50	5.20	9.40	0.30	FLFS	
JY-13	13	2.48	5.00	10.10	0.21	FLFS	
JY-14	14	2.50	5.20	10.80	0.21	FGFS	
JY-15	15	2.50	5.20	11.85	0.78	FLFS	
JY-16	16	2.50	5.10	12.80	0.64	FLFS	
JY-17	17	2.50	5.20	11.55	0.16	FLFS	
JY-18	18	2.50	5.20	12.40	0.32	FLFS	
JY-19	19	2.49	5.20	10.10	0.20	FLFS	
JY-20	20	2.50	5.20	13.60	1.10	FLFS	
JY-21	21	2.50	5.20	7.07	0.09	FLFS	
JY-22	22	2.50	5.10	5.65	0.08	FLFS	
JY-23	23	2.50	5.20	6.68	0.05	FLFS	
JY-24	24	2.49	5.20	7.50	0.09	FGFS	
JY-25	25	2.49	5.20	9.55	0.22	FLFS	
JY-26	26	2.50	5.10	8.87	0.09	FLFS	
JY-27	27	2.48	5.20	7.16	0.10	FLFS	
JY-28	28	2.49	5.20	6.49	0.08	FLFS	
JY-29	29	2.50	5.10	6.98	0.06	FLFS	
JY-30	30	2.50	5.20	6.15	0.05	FLFS	
JY-31	31	2.50	5.20	7.35	0.05	FLFS	
JY-32	32	2.50	5.20	8.20	0.09	FLFS	
JY-33	33	2.50	5.20	7.00	0.18	FFLS	
JY-34	34	2.49	5.20	8.17	0.09	FFLS	
JY-35	35	2.50	5.20	6.32	0.08	FFLS	
JY-36	36	2.50	5.20	8.17	0.12	FFLS	
JY-37	37	2.50	5.10	7.22	0.08	FFLS	
JY-38	38	2.49	5.00	9.99	0.05	FFLS	
JY-39	39	2.50	5.00	6.69	0.02	FFLS	
JY-40	40	2.50	5.00	6.52	0.06	FFLS	
JY-41	41	2.50	5.00	10.40	0.11	FFLS	
JY-42	42	2.50	5.00	7.04	0.02	FFLS	
JY-43	43	2.49	5.00	9.81	0.10	FFLS	
JY-44	44	2.50	5.20	2.41	0.01	FFLS	
JY-45	45	2.50	5.20	3.64	0.03	FFLS	
JY-46	46	2.50	5.20	9.01	0.43	FFLS	
JY-47	47	2.50	5.00	3.09	0.02	FFLS	
JY-48	48	2.50	5.20	4.78	0.02	FFLS	
JY-49	49	2.50	5.20	5.48	0.05	FFLS	
JY-50	50	2.49	5.20	4.91	0.04	FFLS	
JY-51	51	2.50	5.20	5.78	0.01	FFLS	
JY-52	52	2.50	5.10	3.47	0.02	FGFS	
FLFS fine-grained lithic feldspar sandstone, FFLS fine-grained feldspar lithic sandstone, FGFS fine-grained feldspar sandstone.

High pressure mercury intrusion and scanning electron microscopy

High pressure mercury intrusion experiment was used to evaluate the micropore throat characteristics of reservoirs quantitatively. This is achieved by observing the pressure changes during mercury injection into the pores, analyzing the characteristics of the capillary pressure curves, and studying the relationship between the intrusion volume of mercury and these characteristics30,31. In this experiment, the Auto Pore IV 9530 fully automated mercury porosimeter was utilized, with a pore diameter measurement range of 3 nm to 1100 μm. Continuous mercury injection was employed, with volume accuracy of less than 0.1 μl for both injection and withdrawal. The experimental procedure followed the national standard GB/T29171-2012, and the maximum mercury injection pressure reached 95.39 MPa.

Scanning electron microscopy (SEM) allows for high-resolution morphological observation and analysis of samples, as well as structural and compositional characterization. It also enables direct observation of the development characteristics of the micro-pore throats in the reservoir32–34. The experiment employed the Japanese Electron JSM-7500F field emission scanning electron microscope, which achieves a secondary-electron image resolution of 1 nm and magnification ranging from 20 to 300,000 times.

Pore parameter extraction technology

The ImageJ software, initially developed by Wayne Rasband at the National Institutes of Health in the United States, is a powerful open-source image processing system written in Java. It was initially applied in the fields of biomedical and agricultural sciences35. Recently, an increasing number of scholars have used it to identify and extract reservoir pores and fracture features36–39. In this study, the ImageJ software was used to process 210 scanning electron microscope images, extracting various pore parameters, including the perimeter, circularity, major axis length, aspect ratio, and solidity.

Statistical methodology

GRA

GRA is to address infinite space problems using finite sequences. It aims to evaluate the correlations between various factors within a system and determine the significance of each factor to the target function. This approach helps to avoid the subjective process of manually assigning weights to factor indicators40. In recent years, GRA has been applied in production forecasting and development plan optimization for tight sandstone reservoirs41–44. The specific steps are as follows.Determine the initial sequence:1 X0={X0(1),X0(2),…,X0(n)},

2 Xi={Xi(1),Xi(2),…,Xi(n)},(i=1,2,…,m),

where X0 is the reference sequence, Xi is the comparative sequence, i is the number of comparative sequences, m is the number of independent variables, and n is the number of samples.

Normalize the data using the extreme value method:3 Xi¯(k)=Xi(k)-minXi(k)maxXi(k)-minXi(k)(i=1,2,3…,m;k=1,2,3,…,n)

Calculate the gray correlation coefficient:4 γ(x0(k),xi(k))=Δmin+ρΔmaxΔik+ρΔmax

5 Δmin=miniminkx0(k)-xi(k)

6 Δmax=maximaxkx0(k)-xi(k)

7 Δik=|x0(k)-xi(k)|

Obtain the gray correlation coefficient matrix:8 γ01(1)γ02(1)⋯γ0m(1)γ01(2)γ02(2)⋯γ0m(2)⋮⋮⋮⋮γ01(n)γ02(n)⋯γ0m(n)

where ρ is the resolution coefficient, which takes values between 0 and 1. A smaller resolution coefficient indicates greater differences between the correlation coefficients and stronger discriminatory power. Usually, ρ is set to 0.5.

Determine the correlation degree. Represent the correlation strength between the series using the average of the n correlation coefficients:9 εoi=1n∑k=1nγoi(K)

where εoi represents the correlation degree between the i-th comparative sequence and the reference sequence.

Determine the weights and rank the correlation degrees. Normalize the correlation degrees to obtain the weight Wi of each comparative sequence:10 Wi=ε0i∑i=1mε0i

AHP

AHP is a methodology that categorizes the factors within a complex problem into interconnected and prioritized levels. This approach facilitates the process of making decisions based on multiple criteria. It is primarily used to determine the weighting coefficients for comprehensive evaluations45–47. The process is as follows. Construction of a judgment matrix: a judgment matrix is constructed to compare the importance of different factors:11 AW=λmaxW

where A is the matrix of pairwise comparisons, W is the weight vector, and λmax is the maximum eigenvalue.

Calculation of weights: the weight vector W is determined using the sum-product method.

Consistency check:12 Ic=λmax-nn-1

13 IR=λ′max-nn-1

14 IcR=IcIR

where n is the number of elements, Ic is the consistency index, IR is the random consistency index, IcR is the consistency ratio, and λ′max is the average of the maximum eigenvalues.

If IcR < 0.10, the consistency of the judgment matrix is considered acceptable.

EWM

EWM is an objective weighting approach that comprehensively examines the underlying patterns and informational value of unprocessed data. It can determine the uncertainty in variables through entropy values, where larger information content corresponds to smaller uncertainty and smaller entropy, and vice versa. The entropy weighting method is characterized by high accuracy and strong objectivity, and many scholars have applied it to oilfield production with good results 48,49. The basic steps are as follows.Normalize the data and calculate the information entropy:15 Ei=-∑j=1NXijlog2Xijlog2N

where Ei is the information entropy of the ith indicator, Xij is the value of the ith indicator on the jth sample, and N is the number of samples.

Calculate the weights:16 Wi=1-Ei∑j=1M(1-Ej)

where Wi is the weight of the ith indicator, Ei is the information entropy of the ith indicator, and M is the number of indicators.

Comprehensive weight coefficient

Weight coefficients can be used to classify and evaluate the reservoir quality effectively, and several methods are currently available to determine the weight coefficients. These include GRA, the expert evaluation method, Q clustering analysis and discriminant analysis, and factor analysis50–52. In this research, a comprehensive weight analysis methodology that integrated AHP and EWM was employed. The key advantage of this approach lies in its amalgamation of the subjective AHP analysis and the objective numerical analysis of EWM. This combination helps to mitigate the influence of subjective factors to a certain extent, thereby enhancing the reliability of the data.17 wcomprehensive=wiAHP×wiEWM∑i=1mwiAHP×wiEWM

where WiAHP is the weight coefficient obtained from the AHP method, and WiEWM is the weight coefficient obtained from the EWM method.

Results and discussion

Evaluation parameter selection

Tight sandstone reservoirs are influenced by deposition, tectonics, and diagenesis.. These reservoirs demonstrate significant heterogeneity and an intricate distribution of micropore throats. The pore structure plays a crucial role in governing the storage and flow behaviour of the reservoir, where the different shape parameters of the pores govern the micropore structure of the rock formation53–57. Considering the characteristics above, this study aimed to provide a quantitative characterization of the reservoir by assessing three key aspects: the pore structure, the physical properties, and the pore shape parameters. Twelve parameters were selected to establish the relationship between the initial production capacity index and the pore structure and shape parameters. The actual initial production capacity of the oilfield was used as the indicator.

Sensitivity parameter selection for pore structure characteristics

The selected 52 samples were subjected to high-pressure mercury intrusion experiments using an Auto Pore IV 9530 automatic mercury porosimeter. The sorting coefficient varied between 1.5 and 2.74, with an average of 2.10. The coefficient of variation ranged between 13.94 and 17.32, with a mean value of 15.54. With an average value of 13.86 MPa, the median pressure varied between 10.5 and 18.79 MPa. The average displacement pressure was 1.23 MPa, ranging between 0.09 and 2.57 MPa. The median radius had a mean value of 0.09 μm and varied from 0.05 to 0.15 μm. With a mean value of 84.52%, the maximum mercury saturation varied from 62.77 to 93.76%. With an average of 34.90%, the mercury withdrawal efficiency varied between 16.7 and 46.6%. Overall, the pore structure of the reservoir in the study area was poor, with uneven sorting and poor connectivity among the pore throats, indicating strong heterogeneity. Correlation analysis was conducted on the initial production and mercury intrusion parameters (Fig. 3), and it was found that the correlation between the production capacity and permeability and porosity was the strongest, with correlation coefficients (R2) of 0.91 and 0.75, respectively. This is mainly because porosity plays a crucial role in determining the size of the pore space within a reservoir, while permeability governs its flow capacity. In the context of tight sandstone reservoirs, the reservoir quality often depends on favourable pore permeability. The sorting coefficient and coefficient of variation provide insights into the uniformity of the distribution of the pore throat sizes. Higher values of these parameters indicate an improved pore structure and increased reservoir productivity. The median radius and median pressure indicate the pore permeability of the reservoir. A larger median radius and smaller median pressure indicate a larger pore space and stronger flow capacity, resulting in a larger oil production capacity. Therefore, the median radius positively correlates with production, while the median pressure is inversely correlated. The displacement pressure is inversely correlated with production (R2 = 0.65). This is because displacement pressure refers to the capillary pressure corresponding to the largest connected pore, and a higher displacement pressure means a higher capillary pressure, making it more difficult for fluid to flow through. This indicates that tight oil has poor flow capacity in the reservoir and is more difficult to accumulate and extract. In conclusion, the initial production capacity is sensitive to the porosity, permeability, sorting coefficient, coefficient of variation, median pressure, median radius, and displacement pressure.Fig. 3 Relationship between initial production and porosity, permeability, selectivity coefficient, coefficient of variation, median pressure, median radius, and displacement pressure.

Selection of pore-shape-sensitive parameters

A total of 210 high-resolution SEM images were captured for the 52 samples. The rock core pores were identified and extracted using ImageJ, obtaining pore shape parameters such as the perimeter, circularity, major axis length, aspect ratio, and solidity (Fig. 4, Table 2). The average values of the identified pore shape parameters for each sample were then calculated. It was found that the pore perimeters of the 52 samples varied between 40.3 and 486.2 μm, with a mean value of 250.5 μm. The circularity ranged between 0.11 and 0.96, with a mean value of 0.31. The major axis lengths of the circumscribed ellipses spanned from 42.52 to 221.19 μm, with an average of 111.67 μm. The aspect ratios ranged from 1.14 to 2.92, and the average value was 2.32. The solidity values ranged between 0.09 and 0.89, with an average of 0.67. In general, the pore shape parameters of the tight sandstone reservoirs exhibited a wide range of variation, with relatively large average perimeters, average major axis lengths of the circumscribed ellipses, aspect ratios, and solidity, and with small average circularity (Fig. 5). This indicates that the pore shapes in tight sandstone are diverse, predominantly irregular and elongated, with few circular pores. Pearson correlation analysis was conducted between the most sensitive parameters for the prioritized pore structure characteristics and the extracted pore shape parameters (Fig. 6). The absolute value of the correlation coefficient always lies between −1 and 1. In this context, a value closer to 1 indicates a stronger positive relationship between the two independent variables, a value closer to -1 indicates a stronger negative relationship between the independent variables, and a value closer to 0 indicates a weak relationship between the variables. A significant and strong correlation (R2 > 0.5) observed between the different shape parameters of the pores and the mercury injection parameters. This suggests that the shape parameters of the pores play a crucial role in determining the pore structures of tight sandstone reservoirs. In general, the perimeter, circularity, and major axis length of the pores displayed a positive correlation with the porosity (Φ), permeability (K), sorting coefficient (Sp), coefficient of variation (Dr), and median radius (R50). Conversely, they exhibited a negative correlation with the median pressure (P50) and displacement pressure (Pd). On the other hand, the aspect ratio and solidity of the pores were inversely proportional to the porosity, permeability, sorting coefficient, coefficient of variation, and median radius. However, they were positively correlated with the median pressure and displacement pressure. Among them, there was a strong positive correlation (R2 = 0.914) between the perimeter and porosity and a relatively strong negative correlation (R2 = –0.766) with the displacement pressure. A larger pore perimeter results in a greater contact area between the reservoir fluid and the solid, facilitating fluid infiltration and storage. Circularity was strongly positively correlated with permeability (R2 = 0.927) and negatively correlated with the displacement pressure (R2 = –0.604). This is because larger circularity indicates a closer approximation to circular pores, which typically exhibit a uniform distribution, resulting in improved connectivity and fluid flow. The major axis length was strongly positively correlated with the permeability and porosity because the major axis length of the circumscribed ellipses of pores affects the connectivity and fluid flow path within the pores. A larger major axis length indicates better connectivity between pores, resulting in a more direct fluid flow path and higher permeability. Moreover, a longer major axis length corresponds to a larger pore size and higher porosity. The aspect ratio exhibited a strong negative correlation with the permeability and selectivity coefficient (R2 = –0.866, R2 = –0.754, respectively) and a strong positive correlation with the displacement pressure (R2 = 0.652). As the aspect ratio increases, the pores become narrower and more uneven, resulting in longer and narrower flow channels, making fluid flow more difficult. As a result, the displacement pressure increases, the selectivity coefficient decreases, and the permeability decreases. Solidity exhibited a strong negative correlation with permeability (R2 = –0.862) and a positive correlation with the displacement pressure (R2 = 0.574). As the solidity increases, the pore shape becomes more concave, and the roundness deteriorates, making fluid flow between the pores more difficult. In conclusion, it can be observed that the perimeter, circularity, major axis of the circumscribed ellipse, aspect ratio, and solidity of the pores are sensitive to various parameters of mercury intrusion.Fig. 4 Visualization of pore extraction results for rock samples. (A) Pore identification (sample no. 1), (B) pore extraction (sample no. 1), (C) pore identification (sample no. 10), (D) pore extraction (sample no. 10), (E) pore identification (sample no. 25), (F) pore extraction (sample no. 25).

Table 2 Extracted pore shape parameters.

Core number	Φ (%)	K (10–3 µm2)	Perimeter(µm)	Circularity	Major axis (µm)	Aspect ratio	Solidity	
1	9.29	0.39	312.34	0.63	112.89	1.61	0.35	
2	9.32	0.14	322.51	0.31	129.11	2.28	0.71	
3	10.36	0.94	352.42	0.89	142.23	1.48	0.14	
4	10.43	0.19	342.81	0.22	148.68	2.17	0.71	
5	11.54	0.23	378.41	0.34	161.37	2.33	0.61	
6	9.08	0.17	272.46	0.29	105.38	2.37	0.76	
7	9.97	0.13	368.14	0.23	139.40	2.12	0.77	
8	8.53	0.12	325.47	0.24	149.89	2.48	0.77	
9	9.09	0.10	330.97	0.23	130.97	2.55	0.77	
10	9.38	0.10	314.48	0.29	123.78	2.29	0.78	
11	12.25	0.26	420.48	0.38	183.94	1.98	0.65	
12	9.40	0.30	310.64	0.32	159.89	2.05	0.63	
13	10.10	0.21	310.04	0.29	139.07	2.32	0.76	
14	10.80	0.21	337.37	0.28	125.77	2.36	0.73	
15	11.85	0.78	385.59	0.75	175.33	1.31	0.22	
16	12.80	0.64	416.32	0.79	195.34	1.62	0.26	
17	11.55	0.16	399.94	0.31	184.82	2.36	0.62	
18	12.40	0.32	288.56	0.61	107.45	1.72	0.26	
19	10.10	0.20	295.68	0.36	132.39	2.22	0.58	
20	13.60	1.10	486.24	0.96	221.19	1.14	0.09	
21	7.07	0.09	122.71	0.19	93.94	2.68	0.88	
22	5.65	0.08	100.77	0.11	82.38	2.59	0.85	
23	6.68	0.05	100.16	0.18	72.42	2.69	0.72	
24	7.50	0.09	120.20	0.17	82.97	2.67	0.83	
25	9.55	0.22	238.05	0.31	120.53	2.21	0.61	
26	8.87	0.09	327.73	0.11	130.77	2.73	0.83	
27	7.16	0.10	219.06	0.24	92.28	2.32	0.77	
28	6.49	0.08	179.13	0.12	85.46	2.71	0.84	
29	6.98	0.06	138.49	0.19	42.52	2.62	0.81	
30	6.15	0.05	164.49	0.12	73.29	2.73	0.89	
31	7.35	0.05	180.40	0.28	87.58	2.41	0.74	
32	8.20	0.09	245.46	0.35	102.76	2.19	0.62	
33	7.00	0.18	218.19	0.48	84.79	1.97	0.38	
34	8.17	0.09	267.47	0.36	114.27	2.16	0.56	
35	6.32	0.08	179.18	0.21	73.20	2.24	0.68	
36	8.17	0.12	264.94	0.43	118.24	2.18	0.62	
37	7.22	0.08	217.58	0.25	83.13	2.11	0.45	
38	9.99	0.05	299.56	0.13	110.26	2.42	0.75	
39	6.69	0.02	232.39	0.17	77.71	2.57	0.82	
40	6.52	0.06	261.47	0.23	115.18	2.78	0.83	
41	10.40	0.11	366.41	0.24	130.67	2.42	0.72	
42	7.04	0.02	228.36	0.17	86.91	2.61	0.79	
43	9.81	0.10	271.22	0.21	134.88	2.66	0.79	
44	2.41	0.01	40.31	0.18	67.51	2.69	0.83	
45	3.64	0.03	51.32	0.12	50.68	2.87	0.82	
46	9.01	0.43	219.38	0.61	80.18	1.59	0.39	
47	3.09	0.02	68.02	0.12	56.38	2.62	0.84	
48	4.78	0.02	119.68	0.16	61.94	2.75	0.88	
49	5.48	0.05	142.86	0.21	76.23	2.46	0.79	
50	4.91	0.04	179.53	0.18	85.38	2.53	0.81	
51	5.78	0.01	226.21	0.14	100.06	2.92	0.84	
52	3.47	0.02	66.69	0.22	63.77	2.69	0.76	

Fig. 5 Distribution of pore shape parameters. (a) Distribution range of pore perimeter and major axis, (b) distribution range of pore circularity, solidity, and aspect ratio.

Fig. 6 Correlations between pore structural parameters and pore shape parameters.

Reservoir classification evaluation

Quantitative classification prediction formula

Based on the results of the GRA, AHP, and EWM, a comprehensive quantitative classification prediction formula was constructed using the superposition principle. This formula was then used to classify and evaluate tight sandstone reservoirs.18 Q=∑i=1nai·bi,N

where Q is the productivity index, ai is the dimensionless weight coefficients of various parameters, bi,N is the dimensionless normalized parameters, and n is the number of parameters.

Determination of weight coefficients

In this study, the initial production rate directly reflecting the reservoir quality was taken as the fundamental sequence. Seven sensitive parameters, namely, the porosity, permeability, sorting coefficient, coefficient of variation, median pressure, median radius, and displacement pressure, were considered as sub-sequences. The principles and steps of GRA were employed to determine the weights of various parameters, thereby assessing the sensitivity of each factor to the initial production rate (Table 3). Combining the correlation degree between the sensitive parameters determined by the gray correlation method and the initial productivity. Then, the parameters were compared in pairs, and values were assigned based on the 9-point scale method. The judgment matrix was obtained by pairwise comparisons of the seven sensitive parameters (Table 4). Subsequently, the weight coefficients were determined using the weighted product method within the AHP (Table 5). Formula (14) shows that the judgment matrix IcR = 0.093 is less than 0.1, meeting the consistency requirements. Subsequently, the EWM analysis method was employed to conduct an objective analysis of each sensitive parameter, resulting in objective weight indices. The comprehensive weight coefficients were calculated using Eq. (17) (Table 5). The formula for the initial productivity and the mercury intrusion sensitivity parameter can be obtained as follows:19 Q1=0.228R50,N+0.217KN+0.191ΦN+0.143Sp,N+0.116Dr,N+0.077P50,N+0.028Pd,N

where ΦN is the normalized porosity, KN is the normalized permeability, SP,N is the normalized sorting coefficient, Dr,N is the normalized coefficient of variation, P50,N is the normalized median pressure, R50,N is the normalized median radius, and Pd,N is the normalized displacement pressure. Table 3 Correlation and ranking of factors with initial yield.

Factor	Φ (%)	R50 (µm)	Sp	Dr	P50	K (10–3 µm2)	Pd (MPa)	
Correlation coefficient (R2)	0.927	0.896	0.894	0.87	0.83	0.824	0.755	
Rank	1	2	3	4	5	6	7	

Table 4 Judgment matrix.

Element	Pd (MPa)	K (10–3 µm2)	P50 (MPa)	Dr	SP	R50 (µm)	Φ (%)	
Pd (MPa)	1	1/3	1/3	1/3	1/3	1/5	1/5	
K (10–3 µm2)	3	1	1/3	1/3	1/3	1/3	1/3	
P50 (MPa)	3	3	1	1/3	1/3	1/3	1/3	
Dr	3	3	3	1	1/3	1/3	1/3	
SP	3	3	3	3	1	1/3	1/3	
R50 (µm)	5	3	3	3	3	1	1/3	
Φ (%)	5	3	3	3	3	3	1	

Table 5 Comprehensive weight coefficients for production capacity.

Method	Pd (MPa)	K (10–3 µm2)	P50 (MPa)	Dr	SP	R50 (µm)	Φ (%)	
AHP	0.038	0.066	0.091	0.120	0.158	0.222	0.306	
EWM	0.087	0.392	0.101	0.115	0.108	0.123	0.074	
Comprehensive weight coefficient	0.028	0.217	0.077	0.116	0.143	0.228	0.191	

Then, using the mercury intrusion parameter as the fundamental sequence, five sensitive parameters related to the pore shape, namely, the perimeter, circularity, major axis length, aspect ratio, and solidity, were considered sub-sequences. The correlation between the mercury intrusion parameters and the pore-shape-sensitive parameters was determined using GRA. The comprehensive weight coefficients for each mercury intrusion parameter were calculated using a combination of the AHP and the EWM (Table 6). Based on these weight coefficients, the correlation formulas between each mercury intrusion parameter and the pore shape parameters were obtained as follows:20 ΦN=0.323CN+0.306MN+0.164PN+0.107AN+0.098SN

21 KN=0.649CN+0.145PN+0.111MN+0.101SN+0.094AN

22 SP,N=0.332MN+0.247AN+0.176CN+0.126SN+0.120PN

23 Dr,N=0.378AN+0.301SN+0.139MN+0.113CN+0.071PN

24 P50,N=0.441AN+0.250SN+0.136MN+0.127PN+0.043CN

25 R50,N=0.311MN+0.182PN+0.179CN+0.177AN+0.150SN

26 Pd,N=0.419AN+0.272SN+0.156MN+0.097PN+0.054CN

Table 6 Comprehensive weight coefficients for mercury pressure parameters.

Method	Perimeter (µm)	Circularity	Major axis (µm)	Aspect ratio	Solidity	
Comprehensive weight coefficient (Φ)	0.164	0.323	0.306	0.107	0.098	
Comprehensive weight coefficient (K)	0.145	0.649	0.111	0.094	0.101	
Comprehensive weight coefficient (SP)	0.120	0.176	0.332	0.247	0.126	
Comprehensive weight coefficient (Dr)	0.071	0.113	0.139	0.378	0.301	
Comprehensive weight coefficient (P50)	0.127	0.043	0.136	0.441	0.250	
Comprehensive weight coefficient (R50)	0.182	0.179	0.311	0.177	0.150	
Comprehensive weight coefficient (Pd)	0.097	0.054	0.156	0.419	0.272	

Combined with Formula (19), the relationship between the initial productivity and pore shape parameters can be obtained:27 Q2=0.286CN+0.231MN+0.206AN+0.147SN+0.13PN

where PN is the normalized perimeter, CN is the normalized circularity, MN is the normalized major axis, AN is the normalized aspect ratio, and SN is the normalized solidity.

Classification scheme and feature evaluation

Based on the indices Q1, which relate initial productivity to high-pressure mercury intrusion sensitivity parameters, and Q2, which relate initial productivity to pore shape parameters, a classification and evaluation scheme for the Chang 8 tight sandstone reservoir have been determined. As depicted in Fig. 7, Q1 for type III reservoirs ranges from 0.1 to 0.31, and Q2 ranges between 0.3 and 0.44. For type II reservoirs, Q1 ranges from 0.32 to 0.47, and Q2 ranges from 0.44 to 0.56. For type I reservoirs, Q1 ≥ 0.5 and Q2 ≥ 0.5. Moreover, type I reservoirs can be further divided into type I1, comprising high-yield reservoirs, and type I2, comprising high-quality reservoirs, with Q1 > 0.7 and Q2 > 0.57 indicating type I1 high-yield reservoirs. Type I reservoirs are considered optimal for the Chang 8 formation, with 15 out of 52 samples belonging to this type, accounting for 28.8%. The characteristics associated with this type of reservoir include favourable pore permeability, featuring an average porosity of 11.1% and permeability of 0.4 × 10–3 µm2. Additionally, these reservoirs possess a low displacement pressure of 0.62 MPa, a low median pressure of 11.79 MPa, and a relatively high median radius of 0.12 µm. The reservoir exhibits good pore throat selectivity, characterized by a large sorting coefficient (2.5) and variation coefficient (16.43). The average pore perimeter of the reservoir is relatively long (360.30 µm), with good circularity (0.50) and a small aspect ratio (1.92). This indicates that the pore shape is more regular and almost circular. Generally, type II displays moderate petrophysical characteristics, characterized by an average porosity of 8.43% and permeability of 0.1 × 10–3 µm2. Within this classification, 19 samples contribute to 36.54% of the dataset. Compared to type I, this reservoir type has a somewhat higher average displacement pressure and median pressure (1.11 MPa and 13.48 MPa, respectively). The median radius is lower (0.10 µm), and the average sorting coefficient and coefficient of variation are 2.41 and 16.18, respectively, indicating moderate sorting. The average pore perimeter of this reservoir type is smaller than that of type I (261.61 µm), with smaller circularity (0.26) and a larger aspect ratio (2.41). Compared to type I, the pores of type II reservoirs exhibit irregular and more elongated shapes. Type III exhibits poorer petrophysical properties, with an average permeability of 0.06 × 10–3 μm2 and porosity of 5.7%, significantly lower than those of type I and type II. There were 18 samples belonging to this type, accounting for 34.62%. This reservoir type has an average displacement pressure of 1.89 MPa and a median pressure of 16.1 MPa, greater than type II. The median radius is the smallest (0.07 µm). The average sorting coefficient and coefficient of variation are 1.81 and 14.7, respectively, indicating poor pore throat sorting. The average pore perimeter is the smallest (147.37 µm), with the poorest circularity (0.19) and the largest aspect ratio (2.56). This indicates that the pores of type III reservoirs are more elongated and slender, making them unfavorable for fluid flow and leading to poor reservoir permeability. In summary, it can be observed that as the reservoir quality deteriorates, the pore structure becomes increasingly worse, and the pore shapes become more complex and variable.Fig. 7 Comprehensive quantitative classification prediction model for the research area of the Chang 8 reservoir.

According to the distribution maps of the well locations and sedimentary microfacies (Figs. 1c, 8), it is observed that type I reservoir wells are mostly found in the northwest of the research region, within the subaqueous distributary channels, exhibiting a zigzag pattern. Most type II reservoir wells are located in the study area's centre, mainly within the middle portions of the subaqueous distributary channel's lateral sand bodies. On the other hand, the relatively poor type III reservoir wells are scattered around the type II reservoirs, with most of them located in the marginal areas adjacent to the interdistributary bay and the edge of the channel’s lateral sand bodies. Therefore, in practical exploration and development, the high-quality reservoirs (type I) in the study area's northwest part should be prioritised.Fig. 8 Planar distribution map of comprehensive quantitative classification for the research area of the Chang 8 reservoir.

Additionally, the main reason for the high productivity of type I1 reservoirs is the higher content of dissolved pores in type I reservoirs. According to Table 7 and Fig. 9, samples 3, 15, 16, and 20 from type I reservoirs exhibit significant development of feldspar dissolution pores, intergranular pores, and a small number of rock particles that dissolve pores. The average absolute contents of feldspar dissolution and intergranular pores are 1.2% and 5.15%, respectively. The average face rate is 0.8%, higher than the other samples. The greater the development of feldspar dissolution and intergranular pores, the larger the flow channels and storage space they provide, thus improving the reservoir’s porosity and permeability, resulting in high-productivity reservoirs. The pore shape parameters of samples 3, 15, 16, and 20 were compared with those of the other samples (Table 2). It was found that these four samples have longer pore perimeters and major axes, larger shape factor (roundness) coefficients, and relatively smaller aspect ratios and concavity. This indicates that high-productivity reservoirs (type I1) have larger pore perimeters, an increased contact area between the pores and reservoir fluids, higher pore circularity, and more circular shapes favourable for fluid flow and storage. Furthermore, as shown in Fig. 8, the four high-productivity wells (JY-3, JY-15, JY-16, JY-20) are all located on the main channel of the subaqueous distributary channel. Therefore, from a macro perspective, thicker sand bodies may be another reason for their high productivity. Table 7 Absolute porosity content of type I reservoir samples.

Sample number	Absolute porosity content (%)	Face rate (%)	
Intergranular pores	Feldspar dissolution pores	Rock debris dissolution pores	Intercrystalline pores	
1	4.5	0.7	0.2	0.3	0.6	
3	4.8	1	0.1	0.4	0.9	
4	3.4	0.5	0.2	0.2	0.4	
5	3.6	0.5	0.1	0.2	0.4	
7	3.5	0.8	0.2	0.1	0.6	
11	3.2	0.4	0.1	0.2	0.3	
12	3.7	0.3	0.2	0.1	0.5	
13	2.7	0.4	0.1	0.2	0.6	
14	3.3	0.8	0.1	0.1	0.1	
15	4.9	1.2	0.1	0.4	0.8	
16	5.2	1.1	0.1	0.3	0.7	
17	3.2	0.6	0.2	0.2	0.3	
18	3.8	0.7	0.1	0.3	0.5	
19	3.0	0.6	0.1	0.1	0.2	
20	5.7	1.5	0.1	0.3	0.8	

Fig. 9 Porosity structure of type I1 reservoir. (A) Intergranular pores, developed dissolution pores (sample no. 3), (B) feldspar dissolution pores (sample no. 20), (C) rock fragment dissolution pores (sample no. 15), (D) intergranular pores, locally developed dissolution pores (sample no. 16).

Model validation

In order to verify the model, 15 coring wells in Jiyuan Chang 8 reservoir were selected. High-pressure mercury intrusion tests, scanning electron microscopy, and thin-section casting experiments were conducted on corresponding samples to extract the pore shape parameters. Next, the comprehensive indices Q1 and Q2, for reservoir categorization, were determined using the GRA, the AHP, and the EWM. Finally, the accuracy of the classification results was compared with that of the existing oil test parameters. As shown in Fig. 10, three wells were classified as type I reservoirs, with an average initial yield of 5.73 t/d. Six wells were classified as type II reservoirs, with an average initial yield lower than type I at 3.52 t/d. One well was misclassified, deviating from the expected value. Five wells were classified as type III reservoirs, with the lowest average initial yield of 1.32 t/d. The quantitative evaluation of the comprehensive parameters matched the actual production capacity results, demonstrating a high matching rate of 93.3%. Compared to conventional models by other scholars for tight sandstone reservoirs, this model establishes a direct connection between actual oilfield production data, microscale pore shape parameters, and pore structure parameters, leading to quantitative reservoir classification evaluation 58–60. It demonstrates higher and more stable classification accuracy.Fig. 10 Comparative analysis of the integrated quantitative classification prediction for the Chang 8 reservoir.

Conclusions

Tight sandstone reservoirs display significant heterogeneity and intricate microscopic pore structures, which impact the accuracy of reservoir assessment. This study employed scanning electron microscopy, thin section analysis, and high-pressure mercury intrusion data as samples. It utilized image recognition technology and machine learning methods to develop a novel classification and evaluation system for tight sandstone reservoirs based on microscopic pore structures. This method utilizes minimal experimental data, is cost-effective, demonstrates relatively high model accuracy, and is particularly suitable for tight sandstone reservoirs. The research conclusions are as follows:By analyzing high pressure mercury parameters, scanning electron microscopy images, and thin sections of the study area in the Chang 8 reservoir, a comprehensive quantitative classification prediction model for tight sandstone reservoirs was established. The model was constructed using twelve sensitive parameters: porosity, permeability, sorting coefficient, coefficient of variation, median pressure, median radius, displacement pressure, pore perimeter, circularity, major axis length, aspect ratio, and solidity, all extracted using image recognition technology.

The case study based on the comprehensive quantitative classification prediction model showed that Q1 ≥ 0.5 and Q2 ≥ 0.5 corresponded to type I reservoirs, while Q1 > 0.7 and Q2 > 0.57 corresponded to type I1 high-yield reservoirs. When 0.32 < Q1 < 0.47 and 0.44 < Q2 < 0.56, a type II reservoir was identified. When 0.1 < Q1 < 0.32 and 0.3 < Q2 < 0.44, a type III reservoir was identified. Additionally, the presence of high-content dissolution pores, intergranular pores, and larger pore perimeters, as well as higher pore circularity, were the main factors contributing to high-yield reservoirs (type I1). The model was validated, achieving an overall accuracy of 93.3%, which indicates its effectiveness in predicting the classification and evaluation of tight reservoirs.

Reservoir quality is influenced by the pore structure characteristics and shape parameters. In tight sandstone reservoirs, the productivity is positively correlated with the porosity, permeability, sorting coefficient, coefficient of variation, and median radius, but negatively correlated with the median pressure and displacement pressure. The perimeter, circularity, and major axis length of the pores are positively correlated with the porosity, permeability, sorting coefficient, coefficient of variation, and median radius, but negatively correlated with the median pressure and displacement pressure.

Type I reservoir wells were primarily found in the northwest of the research region, within the subaqueous distributary channels, exhibiting a zigzag pattern. The majority of type II reservoir wells were located in the study area's center, mostly within the middle portions of the subaqueous distributary channel’s lateral sand bodies. In contrast, the relatively inferior type III reservoir wells were dispersed among the type II reservoirs, primarily situated in the marginal zones bordering the interdistributary bay and the periphery of the channel’s lateral sand bodies. Therefore, in terms of practical exploration and development, priority should be given to the superior reservoirs (type I) in the northwestern sector of the research region.

The evaluation results of the quantitative classification of tight sandstone reservoirs using machine learning are generally consistent with previous multiparameter conventional evaluation studies. However, this approach effectively integrates macroscopic and microscopic parameters, resulting in higher model accuracy, easier operation, and lower costs. It is particularly suitable for large-scale quality assessments of tight sandstone reservoirs, offering essential guidance for further exploration in the study area and other similar reservoirs.

Abbreviations

AHP Analytic hierarchy process

GRA Grey relational analysis

EWM Entropy weight method

XRD X-ray diffraction

SEM Scanning electron microscopy

FLFS Fine-grained lithic feldspar sandstone

FFLS Fine-grained feldspar lithic sandstone

FGFS Fine-grained feldspar sandstone

GCA-CE-MGPK Grey correlation analysis, clustering ensemble, and the Kriging model combined with macro geological parameters

Acknowledgements

This research was sponsored by Natural Science Basic Research Plan in Shaanxi Province of China (Grant No. 2017JM4013; Grant No. 2020JQ-798).

Author contributions

Xinglei Song: Investigation, Formal analysis, Conceptualization, Data Curation, Writing-Original Draft; Congjun Feng: Writing-Review & Editing, Supervision, Funding acquisition,Methodology; Teng Li: Investigation, Resources, Data Curation, Writing-Review & Editing; Qin Zhang: Investigation, Resources, Data Curation; Xinhui Pan: Supervision, Project administration; Mengsi Sun: Supervision, Writing-Review & Editing, Project administration; Yanlong Ge: Investigation, Resources, Data Curation. All authors have read and agreed to the published version of the manuscript.

Data availability

The data that support the findings of this study are available from the corresponding author upon reasonable request.

Competing interests

The authors declare no competing interests.

Publisher's note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
==== Refs
References

1. Ledingham Glen W Santigo, pool. California: Geological note AAPG Bull. 1947 31 11 2063 2067
Ledingham, & Glen, W. Santigo, pool. California: Geological note. AAPG Bull. 31(11), 2063–2067 (1947).
2. Wang QR Tao SZ Guan P Progress in research and exploration & development of shale oil in continental basins in China Nat. Gas Geosci. 2020 31 3 417 427
Wang, Q. R., Tao, S. Z. & Guan, P. Progress in research and exploration & development of shale oil in continental basins in China. Nat. Gas Geosci. 31(3), 417–427 (2020).
3. Zou CN Zhang GS Yang Z Tao SZ Hou LH Zhu RK Yuan XJ Ran QQ Li DH Wang ZP Geological concepts, characteristics, resource potential and key techniques of unconventional hydro-carbon: On unconventional petroleum geology Pet. Explor. Dev. 2013 40 4 385 399 10.1016/S1876-3804(13)60053-1
Zou, C. N. et al. Geological concepts, characteristics, resource potential and key techniques of unconventional hydro-carbon: On unconventional petroleum geology. Pet. Explor. Dev. 40(4), 385–399 (2013).10.1016/S1876-3804(13)60053-1
4. Zou CN Yang Z Zhu RK Zhang GS Hou LH Wu ST Tao SZ Yuan XJ Dong DZ Wang YM Progress in China’s unconventional oil & gas exploration and development and theoretical technologies Acta Geol. Sin. 2015 89 6 979 1007
Zou, C. N. et al. Progress in China’s unconventional oil & gas exploration and development and theoretical technologies. Acta Geol. Sin. 89(6), 979–1007 (2015).
5. Zhou N Lu S Wang M Huang W Xiao D Jiao C Wang J Tian W Zhou L Chen F Liu W Wang Z Limits and grading evaluation criteria of tight oil reservoirs in typical continental basins of China Petrol. Explor. Dev. 2021 48 05 1089 1100 10.1016/S1876-3804(21)60093-9
Zhou, N. et al. Limits and grading evaluation criteria of tight oil reservoirs in typical continental basins of China. Petrol. Explor. Dev. 48(05), 1089–1100 (2021).10.1016/S1876-3804(21)60093-9
6. Zhao W Hu S Hou L Yang T Li X Guo B Yang Z Types and resource potential of continental shale oil in China and its boundary with tight oil Petrol. Explor. Dev. 2020 47 01 1 11 10.1016/S1876-3804(20)60001-5
Zhao, W. et al. Types and resource potential of continental shale oil in China and its boundary with tight oil. Petrol. Explor. Dev. 47(01), 1–11 (2020).10.1016/S1876-3804(20)60001-5
7. Sun L Zou C Jia A Wei Y Zhu R Wu S Guo Z Development characteristics and orientation of tight oil and gas in China Petrol. Explor. Dev. 2019 46 06 1073 1087 10.1016/S1876-3804(19)60264-8
Sun, L. et al. Development characteristics and orientation of tight oil and gas in China. Petrol. Explor. Dev. 46(06), 1073–1087 (2019).10.1016/S1876-3804(19)60264-8
8. Xiang F Zhi Y Weipeng Y Xuguang G Yanxiong W Juntian L Classification evaluation criteria and exploration potential of tight oil resources in key basins of China J. Nat. Gas Geosci. 2019 4 6 309 319 10.1016/j.jnggs.2019.11.002
Xiang, F. et al. Classification evaluation criteria and exploration potential of tight oil resources in key basins of China. J. Nat. Gas Geosci. 4(6), 309–319 (2019).10.1016/j.jnggs.2019.11.002
9. Gao X Chen J XU R Zhen Z Zeng X Chen X Cui L Research progress and prospect of the materials of bipolar plates for proton exchange membrane fuel cells (PEMFCs)[J] International Journal of Hydrogen Energy 2024 50 711 743 10.1016/j.ijhydene.2023.09.005
Gao, X., Chen, J., Xu, R., Zhen, Z., Zeng, X., Chen, X. & Cui, L. Research progress and prospect of the materials of bipolar plates for proton exchange membrane fuel cells (PEMFCs)[J]. International Journal of Hydrogen Energy. 50, 711–743 (2024).10.1016/j.ijhydene.2023.09.005
10. Wang J Wu S Li Q Guo Q An investigation into pore structure fractal characteristics in tight oil reservoirs: A case study of the Triassic tight sandstone with ultra-low permeability in the Ordos Basin, China Arab. J. Geosci. 2020 13 18 961 10.1007/s12517-020-05928-0
Wang, J., Wu, S., Li, Q. & Guo, Q. An investigation into pore structure fractal characteristics in tight oil reservoirs: A case study of the Triassic tight sandstone with ultra-low permeability in the Ordos Basin, China. Arab. J. Geosci. 13(18), 961 (2020).10.1007/s12517-020-05928-0
11. Gao, H., Cao, J., Wang, C., He, M., Dou, L., Huang, X. & Li, T. Comprehensive characterization of pore and throat system for tight sandstone reservoirs and associated permeability determination method using SEM, rate controlled mercury and high pressure mercury. J. Petrol. Sci. Eng. 174 (2018).
12. Gao H Wang C Cheng Z Li T Dou L Zhao K Xue J Luo K Effect of pressure pulse stimulation on imbibition displacement within a tight sandstone reservoir with local variations in porosity Geoenergy Sci. Eng. 2023 226 211811 10.1016/j.geoen.2023.211811
Gao, H. et al. Effect of pressure pulse stimulation on imbibition displacement within a tight sandstone reservoir with local variations in porosity. Geoenergy Sci. Eng. 226, 211811 (2023).10.1016/j.geoen.2023.211811
13. Wang C Gao H Gao Y Fan H Influence of pressure on spontaneous imbibition in tight sandstone reservoirs Energy Fuels 2020 34 8 9275 9282 10.1021/acs.energyfuels.0c00639
Wang, C., Gao, H., Gao, Y. & Fan, H. Influence of pressure on spontaneous imbibition in tight sandstone reservoirs. Energy Fuels 34(8), 9275–9282 (2020).10.1021/acs.energyfuels.0c00639
14. Wang C Li T Gao H Zhao J Gao Y Quantitative study on the blockage degree of pores due to asphaltene precipitation in low-permeability reservoirs with NMR technique J. Petrol. Sci. Eng. 2018 163 703 711 10.1016/j.petrol.2017.11.021
Wang, C., Li, T., Gao, H., Zhao, J. & Gao, Y. Quantitative study on the blockage degree of pores due to asphaltene precipitation in low-permeability reservoirs with NMR technique. J. Petrol. Sci. Eng. 163, 703–711 (2018).10.1016/j.petrol.2017.11.021
15. Gao H Zhugeng B Wang C Cao J He M Dou L Effects of pore structure and salinity on the imbibition of shale samples using physical simulation and NMR technique: A case from Chang 7 shale, Ordos basin Simulation. 2021 97 2 167 173 10.1177/0037549719830194
Gao, H. et al. Effects of pore structure and salinity on the imbibition of shale samples using physical simulation and NMR technique: A case from Chang 7 shale, Ordos basin. Simulation. 97(2), 167–173 (2021).10.1177/0037549719830194
16. Wei H Xie R Guo J Liu J Xu C Wang Y Classification of tight sandstone reservoirs based on the nuclear magnetic resonance T2 distribution: A case study on the Shaximiao Formation in Central Sichuan, China Energy Fuels 2022 36 10803 10812 10.1021/acs.energyfuels.2c01612
Wei, H. et al. Classification of tight sandstone reservoirs based on the nuclear magnetic resonance T2 distribution: A case study on the Shaximiao Formation in Central Sichuan, China. Energy Fuels 36, 10803–10812 (2022).10.1021/acs.energyfuels.2c01612
17. Xu J Yu C Qi D Yang R Yu H Zhang W Wang T Characteristics and controlling factors of tight gas sandstones from the Upper Shanxi and Lower Shihezi Formations in the Northern Sulige Area, Ordos Basin, China Energy Fuels 2023 37 20 15712 15729 10.1021/acs.energyfuels.3c02518
Xu, J. et al. Characteristics and controlling factors of tight gas sandstones from the Upper Shanxi and Lower Shihezi Formations in the Northern Sulige Area, Ordos Basin, China. Energy Fuels 37(20), 15712–15729 (2023).10.1021/acs.energyfuels.3c02518
18. Wu X Su Y Zhang C Xin Y Chen X Li N Huang R Tang B Zhao X A novel evaluation method of dolomite reservoir using electrical image logs: The Cambrian dolomites in Tarim Basin, China Geoenergy Sci. Eng. 2024 233 212509 10.1016/j.geoen.2023.212509
Wu, X. et al. A novel evaluation method of dolomite reservoir using electrical image logs: The Cambrian dolomites in Tarim Basin, China. Geoenergy Sci. Eng. 233, 212509 (2024).10.1016/j.geoen.2023.212509
19. Zhang Q Zhang Y Wang B Yin S Wu X Yuan C Comprehensive evaluation and reservoir classification in the Quan 3 member of the Cretaceous Quantou Formation in the Fuxin Uplift, Songliao Basin Front. Earth Sci. 2022 10 1016924 10.3389/feart.2022.1016924
Zhang, Q. et al. Comprehensive evaluation and reservoir classification in the Quan 3 member of the Cretaceous Quantou Formation in the Fuxin Uplift, Songliao Basin. Front. Earth Sci. 10, 1016924 (2022).10.3389/feart.2022.1016924
20. Talib M Durrani MZA Palekar AH Sarosh B Rahman SA Quantitative characterization of unconventional (tight) hydrocarbon reservoir by integrating rock physics analysis and seismic inversion: A case study from the Lower Indus Basin of Pakistan Acta Geophys. 2022 70 6 2715 2731 10.1007/s11600-022-00885-6
Talib, M., Durrani, M. Z. A., Palekar, A. H., Sarosh, B. & Rahman, S. A. Quantitative characterization of unconventional (tight) hydrocarbon reservoir by integrating rock physics analysis and seismic inversion: A case study from the Lower Indus Basin of Pakistan. Acta Geophys. 70(6), 2715–2731 (2022).10.1007/s11600-022-00885-6
21. Fang X Zhu G Yang Y Li F Feng H Quantitative method of classification and discrimination of a porous carbonate reservoir integrating k-means clustering and Bayesian theory Acta Geol. Sin. (Beijing) 2023 97 1 176 189 10.1111/1755-6724.14941
Fang, X., Zhu, G., Yang, Y., Li, F. & Feng, H. Quantitative method of classification and discrimination of a porous carbonate reservoir integrating k-means clustering and Bayesian theory. Acta Geol. Sin. (Beijing) 97(1), 176–189 (2023).10.1111/1755-6724.14941
22. Li Y Niu D Zhang Y Wu H Fu H Wang Z Microscopic pore-throat grading evaluation in a tight oil reservoir using machine learning: A case study of the Fuyu oil layer in Bayanchagan area, Songliao Basin central depression Earth Sci. Inform. 2021 14 2 601 617 10.1007/s12145-021-00594-6
Li, Y. et al. Microscopic pore-throat grading evaluation in a tight oil reservoir using machine learning: A case study of the Fuyu oil layer in Bayanchagan area, Songliao Basin central depression. Earth Sci. Inform. 14(2), 601–617 (2021).10.1007/s12145-021-00594-6
23. Wang Z Tang H Hou Y Shi H Li J Yang T Feng Y Meng W Quantitative evaluation of unconsolidated sandstone heavy oil reservoirs based on machine learning Geol. J. (Chichester, England). 2023 58 6 2321 2341
Wang, Z. et al. Quantitative evaluation of unconsolidated sandstone heavy oil reservoirs based on machine learning. Geol. J. (Chichester, England). 58(6), 2321–2341 (2023).
24. Fan J Shi J Wan X Xie Q Wang C Classification evaluation method for Chang 7 oil group of Yanchang formation in Ordos Basin J. Pet. Explor. Prod. Te. 2021 12 825 834
Fan, J., Shi, J., Wan, X., Xie, Q. & Wang, C. Classification evaluation method for Chang 7 oil group of Yanchang formation in Ordos Basin. J. Pet. Explor. Prod. Te. 12, 825–834 (2021).
25. Niu D Li Y Zhang Y Sun P Wu H Fu H Wang Z Multi-scale classification and evaluation of shale reservoirs and 'sweet spot' prediction of the second and third members of the Qingshankou Formation in the Songliao Basin based on machine learning J. Petrol Sci. Eng. 2022 216 110678 10.1016/j.petrol.2022.110678
Niu, D. et al. Multi-scale classification and evaluation of shale reservoirs and “sweet spot” prediction of the second and third members of the Qingshankou Formation in the Songliao Basin based on machine learning. J. Petrol Sci. Eng. 216, 110678 (2022).10.1016/j.petrol.2022.110678
26. Li C Liu G Cao Z Sun M You Y Liu N Oil charging pore throat threshold and accumulation effectiveness of tight sandstone reservoir using the physical simulation experiments combined with NMR J. Petrol. Sci. Eng. 2022 208 109 338
Li, C. et al. Oil charging pore throat threshold and accumulation effectiveness of tight sandstone reservoir using the physical simulation experiments combined with NMR. J. Petrol. Sci. Eng. 208, 109–338 (2022).
27. Li S Yao J Mou W Luo A Wang Q Deng X Chu M Li Y Yan C The dissolution characteristics of the Chang 8 tight reservoir and its quantitative influence on porosity in the Jiyuan area, Ordos Basin, China J. Nat. Gas Geosci. 2018 3 2 95 108 10.1016/j.jnggs.2018.04.002
Li, S. et al. The dissolution characteristics of the Chang 8 tight reservoir and its quantitative influence on porosity in the Jiyuan area, Ordos Basin, China. J. Nat. Gas Geosci. 3(2), 95–108 (2018).10.1016/j.jnggs.2018.04.002
28. Song X Gao H Feng C Yi P Wang C Li T Analysis of the influence of micro-pore structure on oil occurrence using nano-CT scanning and nuclear magnetic resonance technology: An example from Chang 8 tight sandstone reservoir, Jiyuan, Ordos Basin Processes 2023 11 11274 10.3390/pr11041127
Song, X. et al. Analysis of the influence of micro-pore structure on oil occurrence using nano-CT scanning and nuclear magnetic resonance technology: An example from Chang 8 tight sandstone reservoir, Jiyuan, Ordos Basin. Processes 11, 11274 (2023).10.3390/pr11041127
29. Tong Q He D Xia Z Huang J Li Y Xu F Guo S Research on sand body architecture at the intersection of a bidirectional sedimentary system in the Jiyuan area of Ordos Basin Sci. Rep. 2023 13 12261 10.1038/s41598-023-28302-y 37507439
Tong, Q. et al. Research on sand body architecture at the intersection of a bidirectional sedimentary system in the Jiyuan area of Ordos Basin. Sci. Rep. 13, 12261 (2023).37507439 10.1038/s41598-023-28302-y
30. Fu S Fang Q Li A Li Z Han J Dang X Han W Accurate characterization of full pore size distribution of tight sandstones by low-temperature nitrogen gas adsorption and high-pressure mercury intrusion combination method Energy Sci. Eng. 2021 9 1 80 100 10.1002/ese3.817
Fu, S. et al. Accurate characterization of full pore size distribution of tight sandstones by low-temperature nitrogen gas adsorption and high-pressure mercury intrusion combination method. Energy Sci. Eng. 9(1), 80–100 (2021).10.1002/ese3.817
31. Li P Sun W Wu B Huang R Gao Y Yan J Huang H Occurrence characteristics and main controlling factors of movable fluids in Chang 81 reservoir, Maling Oilfield, Ordos Basin, China J. Petrol. Explor. Prod. Technol. 2018 9 1 17 29 10.1007/s13202-018-0471-2
Li, P. et al. Occurrence characteristics and main controlling factors of movable fluids in Chang 81 reservoir, Maling Oilfield, Ordos Basin, China. J. Petrol. Explor. Prod. Technol. 9(1), 17–29 (2018).10.1007/s13202-018-0471-2
32. Li C Chen G Li X Zhou Q Sun Z The occurrence of tight oil in the Chang 8 lacustrine sandstone of the Huaqing area, Ordos Basin, China: Insights into the content of adsorbed oil and its controlling factors J. Nat. Gas Geosci. 2022 7 1 27 37 10.1016/j.jnggs.2021.11.001
Li, C., Chen, G., Li, X., Zhou, Q. & Sun, Z. The occurrence of tight oil in the Chang 8 lacustrine sandstone of the Huaqing area, Ordos Basin, China: Insights into the content of adsorbed oil and its controlling factors. J. Nat. Gas Geosci. 7(1), 27–37 (2022).10.1016/j.jnggs.2021.11.001
33. Gong Y Liu K Pore throat size distribution and oiliness of tight sands-A case study of the Southern Songliao Basin, China J. Petrol. Sci. Eng. 2020 184 106508 10.1016/j.petrol.2019.106508
Gong, Y. & Liu, K. Pore throat size distribution and oiliness of tight sands-A case study of the Southern Songliao Basin, China. J. Petrol. Sci. Eng. 184, 106508 (2020).10.1016/j.petrol.2019.106508
34. Liu Y Dong X Chen Z Hou Y Luo Q Chen Y A novel experimental investigation on the occurrence state of fluids in microscale pores of tight reservoirs J. Petrol. Sci. Eng. 2021 196 107656 10.1016/j.petrol.2020.107656
Liu, Y. et al. A novel experimental investigation on the occurrence state of fluids in microscale pores of tight reservoirs. J. Petrol. Sci. Eng. 196, 107656 (2021).10.1016/j.petrol.2020.107656
35. Sandhya N Baviskar A A quick & automated method for measuring Am. Biol. Teach. 2011 73 9 554 556 10.1525/abt.2011.73.9.9
Sandhya, N. & Baviskar, A. A quick & automated method for measuring. Am. Biol. Teach. 73(9), 554–556 (2011).10.1525/abt.2011.73.9.9
36. Curtis ME Cardott BJ Sondergeld CH Development to for organic porosity in the Woodford shale with increasing thermal maturity Int. J. Coal Geol. 2012 26 31 26 30 10.1016/j.coal.2012.08.004
Curtis, M. E., Cardott, B. J. & Sondergeld, C. H. Development to for organic porosity in the Woodford shale with increasing thermal maturity. Int. J. Coal Geol. 26(31), 26–30 (2012).10.1016/j.coal.2012.08.004
37. Keller LM Schuetz P Erni R Characterization of multi-scale micro-structural features in opalinus clay Microporous Mesoporous Mater. 2013 83 84 90
Keller, L. M., Schuetz, P. & Erni, R. Characterization of multi-scale micro-structural features in opalinus clay. Microporous Mesoporous Mater. 83, 84–90 (2013).
38. Jin L Mathur R Rother G Evolution of porosity and geochemistry in Mar cell us formation black shale during weathering Chem. Geol. 2013 50 51 56
Jin, L. et al. Evolution of porosity and geochemistry in Mar cell us formation black shale during weathering. Chem. Geol. 50, 51–56 (2013).
39. Rine JM Smart E Dorsey W Comparison of porosity distribution with in selected north American shale units by SEM examination of argon-ion-milled samples Electron Microsc. Shale Hydrocarbon Reserv. AAPG Memoir. 2013 102 137 152
Rine, J. M. et al. Comparison of porosity distribution with in selected north American shale units by SEM examination of argon-ion-milled samples. Electron Microsc. Shale Hydrocarbon Reserv. AAPG Memoir. 102, 137–152 (2013).
40. Zhao JY An XP Wang J Fan JM Kang XM Tan XQ Li WQ A quantitative evaluation for well pattern adaptability in ultra-low permeability oil reservoirs: A case study of Triassic Chang 6 and Chang 8 reservoirs in Ordos Basin Pet. Explor. Dev. 2018 45 3 482 488 10.1016/S1876-3804(18)30055-7
Zhao, J. Y. et al. A quantitative evaluation for well pattern adaptability in ultra-low permeability oil reservoirs: A case study of Triassic Chang 6 and Chang 8 reservoirs in Ordos Basin. Pet. Explor. Dev. 45(3), 482–488 (2018).10.1016/S1876-3804(18)30055-7
41. Dong Q Dai Yin Y Ya Zhou Z Fine classification of ultra-low permeability reservoirs around the Placanticline of Da Qing oilfield (PR of China) J. Petrol. Sci. Eng. 2019 174 1042 1052 10.1016/j.petrol.2018.12.008
Dong, Q., Dai Yin, Y. & Ya Zhou, Z. Fine classification of ultra-low permeability reservoirs around the Placanticline of Da Qing oilfield (PR of China). J. Petrol. Sci. Eng. 174, 1042–1052 (2019).10.1016/j.petrol.2018.12.008
42. Gao Y Liu L Liu H Zheng F Wu L Zhou J Qu D Application of an analytic hierarchy process to hydro-carbon accumulation coefficient estimation Petrol. Sci. 2010 7 3 337 346 10.1007/s12182-010-0075-z
Gao, Y. et al. Application of an analytic hierarchy process to hydro-carbon accumulation coefficient estimation. Petrol. Sci. 7(3), 337–346 (2010).10.1007/s12182-010-0075-z
43. Liu Y Zhang X Shi J A reservoir quality evaluation approach for tight sandstone reservoirs based on the gray correlation algorithm: A case study of the Chang 6 layer in the W area of the as oilfield, Ordos Basin Energy Explor. Exploit. 2021 39 4 1027 1056 10.1177/0144598721998510
Liu, Y. et al. A reservoir quality evaluation approach for tight sandstone reservoirs based on the gray correlation algorithm: A case study of the Chang 6 layer in the W area of the as oilfield, Ordos Basin. Energy Explor. Exploit. 39(4), 1027–1056 (2021).10.1177/0144598721998510
44. Shi B Chang X Yin W Li Y Mao L Quantitative evaluation model for tight sandstone reservoirs based on statistical methods—A case study of the Triassic Chang 8 tight sandstones, Zhenjing area, Ordos Basin, China J. Petrol. Sci. Eng. 2019 173 601 616 10.1016/j.petrol.2018.10.035
Shi, B., Chang, X., Yin, W., Li, Y. & Mao, L. Quantitative evaluation model for tight sandstone reservoirs based on statistical methods—A case study of the Triassic Chang 8 tight sandstones, Zhenjing area, Ordos Basin, China. J. Petrol. Sci. Eng. 173, 601–616 (2019).10.1016/j.petrol.2018.10.035
45. Liu B The analytic hierarchy process for the reservoir evaluation in Chaoyanggou oilfield Adv. Petrol. Explor. Dev. 2014 6 46 50
Liu, B. The analytic hierarchy process for the reservoir evaluation in Chaoyanggou oilfield. Adv. Petrol. Explor. Dev. 6, 46–50 (2014).
46. Shang YZ Application of analytical hierarchy process in the low-grade oil reservoirs evaluation Daqing Petrol. Geol. Oilfield Dev. 2014 33 55 59
Shang, Y. Z. Application of analytical hierarchy process in the low-grade oil reservoirs evaluation. Daqing Petrol. Geol. Oilfield Dev. 33, 55–59 (2014).
47. Xi Y Li Y Liu J Wu S Lu N Liao G Wang Q Application of analytic hierarchy process in mineral prospecting prediction based on an integrated geology—aerogeophysics—geochemistry model Minerals 2023 13 7 978 10.3390/min13070978
Xi, Y. et al. Application of analytic hierarchy process in mineral prospecting prediction based on an integrated geology—aerogeophysics—geochemistry model. Minerals 13(7), 978 (2023).10.3390/min13070978
48. Lai F Luo H Qin D Crushability evaluation of shale gas reservoir based on analytic hierarchy process Spec. Oil Gas Reserv. 2018 25 3 154 159
Lai, F. et al. Crushability evaluation of shale gas reservoir based on analytic hierarchy process. Spec. Oil Gas Reserv. 25(3), 154–159 (2018).
49. Elhaj, M.A., Imtiaz, S. A., Naterer, G. F. & Zendehboudi, S. Production optimization of hydrocarbon reservoirs by entropy generation minimization. J. Nat. Gas Sci. Eng. 83, 103538 (2020).
50. Szabo NP Cluster analysis of core measurements using heterogeneous data sources: An application to complex Miocene reservoirs J. Petrol. Sci. Eng. 2019 178 575 585 10.1016/j.petrol.2019.03.067
Szabo, N. P. et al. Cluster analysis of core measurements using heterogeneous data sources: An application to complex Miocene reservoirs. J. Petrol. Sci. Eng. 178, 575–585 (2019).10.1016/j.petrol.2019.03.067
51. Oliveira GP Santos MD Roque WL Constrained clustering approaches to identify hydraulic flow units in petroleum reservoirs J. Petrol. Sci. Eng. 2020 186 106732 10.1016/j.petrol.2019.106732
Oliveira, G. P., Santos, M. D. & Roque, W. L. Constrained clustering approaches to identify hydraulic flow units in petroleum reservoirs. J. Petrol. Sci. Eng. 186, 106732 (2020).10.1016/j.petrol.2019.106732
52. Jia A Wei Y Jin Y Progress in key technologies for evaluating marine shale gas development in China Petrol. Explor. Dev. 2016 43 6 1035 1042 10.1016/S1876-3804(16)30120-3
Jia, A., Wei, Y. & Jin, Y. Progress in key technologies for evaluating marine shale gas development in China. Petrol. Explor. Dev. 43(6), 1035–1042 (2016).10.1016/S1876-3804(16)30120-3
53. Xiao L Bi L Yi T Lei Y Wei Q Pore structure characteristics and influencing factors of tight reservoirs controlled by different provenance systems: A case study of the Chang 7 members in Heshui and Xin'anbian of the Ordos Basin Energies 2023 16 34108 10.3390/en16083410
Xiao, L., Bi, L., Yi, T., Lei, Y. & Wei, Q. Pore structure characteristics and influencing factors of tight reservoirs controlled by different provenance systems: A case study of the Chang 7 members in Heshui and Xin’anbian of the Ordos Basin. Energies 16, 34108 (2023).10.3390/en16083410
54. Dong J Huang Z Chen J Li T Zhao J Pan Y Qu T Pore structure and fractal characteristics of tight sandstone: A case study for Huagang Formation in the Xihu Sag, East China Sea Basin, China Energies 2023 16 20134 10.3390/en16042013
Dong, J. et al. Pore structure and fractal characteristics of tight sandstone: A case study for Huagang Formation in the Xihu Sag, East China Sea Basin, China. Energies 16, 20134 (2023).10.3390/en16042013
55. Gao J Wei H Zhou R Ren Q Tian T Ning B Ding X Zhang Z Study on the coupling law between pore-scale fluid flow capacity and pore-throat configuration in tight sandstone reservoirs Geofluids 2023 2023 1 1693773
Gao, J. et al. Study on the coupling law between pore-scale fluid flow capacity and pore-throat configuration in tight sandstone reservoirs. Geofluids 2023(1), 1693773 (2023).
56. Zhang R Jiang S Zhang L Wang H Zhang T Yu R Zhang L Chima FU Microscopic pore structures and their controlling factors of the lower carboniferous Luzhai Shale in Guizhong depression, China Geofluids 2023 2023 8890709 10.1155/2023/8890709
Zhang, R. et al. Microscopic pore structures and their controlling factors of the lower carboniferous Luzhai Shale in Guizhong depression, China. Geofluids 2023, 8890709 (2023).10.1155/2023/8890709
57. Du M Yang Z Yang S Feng C Wang G Jia N Li H Shi X Study on the quantitative characterization and heterogeneity of pore structure in deep ultra-high pressure tight glutenite reservoirs Minerals 2023 13 6015 10.3390/min13050601
Du, M. et al. Study on the quantitative characterization and heterogeneity of pore structure in deep ultra-high pressure tight glutenite reservoirs. Minerals 13, 6015 (2023).10.3390/min13050601
58. Wu BH Xie RH Xiao LZ Guo JF Jin GW Fu JW Integrated classification method of tight sandstone reservoir based on principal component analysis-simulated annealing genetic algorithm-fuzzy cluster means Petrol. Sci. 2023 20 5 2747 2758 10.1016/j.petsci.2023.04.014
Wu, B. H. et al. Integrated classification method of tight sandstone reservoir based on principal component analysis-simulated annealing genetic algorithm-fuzzy cluster means. Petrol. Sci. 20(5), 2747–2758 (2023).10.1016/j.petsci.2023.04.014
59. Lu X Xing X Hu K Zhou B Classification and evaluation of tight sandstone reservoirs based on MK-SVM Processes. 2023 11 9 2678 10.3390/pr11092678
Lu, X., Xing, X., Hu, K. & Zhou, B. Classification and evaluation of tight sandstone reservoirs based on MK-SVM. Processes. 11(9), 2678 (2023).10.3390/pr11092678
60. Qiu X Ding L Liu J Yan Z Bao Y Tan C Quantitative evaluation of reservoir quality of tight oil sandstones in Chang 7 member of Ordos Basin Front. Earth Sci. 2023 10 1046489 10.3389/feart.2022.1046489
Qiu, X. et al. Quantitative evaluation of reservoir quality of tight oil sandstones in Chang 7 member of Ordos Basin. Front. Earth Sci. 10, 1046489 (2023).10.3389/feart.2022.1046489
