A New Method to Determine the Lower Limit of Reservoir Physical Properties—Corrected Minimum Flow Pore-throat Radius Method

. The lower limit of reservoir physical properties is an important parameter for identifying reservoirs and determining effective thickness in reserves evaluation, and is also an important basis for selecting perforated test intervals in oilfield exploration and development. There are many methods to determine the lower limit of reservoir physical properties, and the minimum flow pore throat radius method is one of the commonly used methods. The method uses 0.1μm as the minimum flow pore-throat radius, and uses this to calibrate the lower limit of reservoir physical properties. However, according to the water film theory, the minimum radius of the reservoir's flowing pore throat is not a definite value, but varies with the displacement dynamics. Therefore, there is no exact basis for using 0.1μm as the minimum flow pore-throat radius, so it needs to be corrected. To this end, a new method for determining the lower limit of reservoir physical properties—the corrected minimum flow pore-throat radius method is proposed. The correction method comprehensively considers the factors of oil and gas accumulation dynamics, and determines the lower limit of reservoir physical properties by obtaining the minimum flow pore-throat radius value suitable for oil and gas accumulation dynamics. A case study of Chang 6 3 reservoir in A Oilfield shows that the minimum flow pore radius of oil and gas determined by the correction method is 0.08 μm, and the lower limit of reservoir physical properties (porosity 9.1%, permeability 0.117 × 10 -3 μm 2 ). The traditional method has a minimum flow pore-throat radius of 0.1 μm and a lower limit of reservoir physical properties (porosity of 9.8% and permeability of 0.133 × 10 -3 μm 2 ). Due to full consideration of the impact of oil and gas accumulation dynamics, the minimum flow pore-throat radius determined by the correction method is more reliable than the traditional method, and the lower limit of the reservoir physical property calibrated by it has practical significance.


1Introduction
The lower limit of reservoir physical properties is the minimum effective porosity and minimum permeability that can store and percolate fluid. It is usually expressed by a certain value of porosity or permeability [1] . The lower limit of reservoir physical properties is an important parameter for identifying reservoirs and determining effective thickness in reserves evaluation, and it is also an important basis for selecting perforated test intervals in oilfield exploration and development [2][3] . In view of the basicity and necessity of the lower limit of reservoir physical properties, research on the lower limit of reservoir physical properties is the focus of oil and gas reservoir geology and engineering. For a long time, many scientists and technicians have conducted a lot of research on the lower limit of reservoir physical properties, and put forward more methods to determine the lower limit of reservoir physical properties  . Among them, the minimum flow pore-throat radius method is one of the commonly used methods to determine the lower limit of reservoir physical properties [24][25][26][27][28][29][30][31] . The method believes that the pores and throats of rocks are spaces and channels for oil and gas storage and flow. Whether oil and gas can flow out of rocks under a certain pressure difference depends on the thickness of the throat. That is, the radius of the rock throat is the key factor that determines whether oil and gas can flow out of the rock under a certain pressure difference. This minimum throat radius that can store oil and gas and allow oil and gas to seep is the minimum flow porethroat radius of oil and gas [24,25,32] . After determining the minimum flow pore-throat radius, the correlation curve of pore-throat radius and porosity and permeability can be drawn according to the principle of statistical analysis, and the corresponding porosity and permeability lower limit values can be calibrated according to the minimum flow pore-throat radius [26][27] . At present, the minimum flow pore-throat radius method generally adopts 0.1μm as the minimum flow pore-throat radius of oil and gas. It is also believed that 0.1μm is equivalent to the thickness of the water film attached to the surface of the water-wet clastic rock. The pore-throat radius is less than this value. Oil and gas are not flowing, and the corresponding porosity and permeability values are the lower limits of reservoir physical properties [24][25][26][27][28][29][30][31] .
The advantage of the minimum flow pore-throat radius method is that the lower limit of the reservoir physical properties and the minimum flow pore-throat radius are linked together, and the lower limit of the reservoir physical properties corresponding to the specific pore-throat radius (0.1μm) is calibrated. The process is simple and clear, easy to master. However, this method also has obvious shortcomings. The theoretical basis for determining the minimum flow porethroat radius as a fixed value (0.1 μm) is not sufficient. According to previous studies on water film theory [33][34] , the thickness of the water film attached to the rock surface is not fixed, and the thickness of the water film is a function of the displacement pressure. As the displacement power increases, the thickness of the water film gradually decreases, and the corresponding minimum flow pore-throat radius of oil and gas must also decrease, which shows that the minimum flow porethroat radius should not be a certain value. From this point of view, it is not appropriate to use 0.1μm as the minimum flow pore-throat radius, and there is a certain deviation in the lower limit of the physical properties calibrated by this. Therefore, how to calculate the minimum flow pore-throat radius value more accurately? It is the key to objectively calibrate the lower limit of reservoir physical properties using the minimum flow pore-throat radius method. To this end, this paper proposes a new method to determine the lower limit of reservoir physical properties-the corrected minimum flow pore-throat radius method (hereinafter referred to as the correction method).

Method principle
Pore throats in reservoirs vary in size and geometry, and pore throats of different orders interweave to form an intricately interconnected system. The complexity of the reservoir pore-throat network, especially the order of pore-throat distribution, will inevitably have a direct and significant effect on the seepage response of the reservoir fluid. When the displacement power is constant, only the part of the pore-throat network with pore-throat size greater than a certain limit will seep, while other parts of the pore-throat network below this limit cannot seep. That is, the fluid seepage in the pore-throat network is selective. This phenomenon can be understood as, as the pore-throat size of the reservoir decreases, the liquid particles are continuously strengthened by the capillary force and the molecular force of the surrounding solid interface, and the seepage resistance increases accordingly. When the pore-throat radius reaches a certain limit, the displacement power and the seepage resistance are in equilibrium. When the pore-throat radius is greater than this limit, the displacement power is greater than the seepage resistance, and the liquid flows. When the pore-throat radius is smaller than this limit, the displacement power is less than the seepage resistance, and the liquid does not flow. Therefore, the pore-throat radius in the displacement equilibrium state is the minimum flow pore-throat radius.
As the displacement power increases, the minimum flow pore-throat radius of the reservoir gradually decreases. When the minimum flow pore-throat radius of the reservoir is reduced to just the full use of oil and gas in the reservoir, the corresponding minimum flow porethroat radius is the minimum flow pore-throat radius of the oil and gas, and the lower limit of physical properties calibrated based on this is the effective lower limit. Because the production process is limited by many factors, it is difficult to accurately determine whether the oil and gas are fully utilized and the minimum flow porethroat radius when the oil and gas is fully utilized. However, considering that the process of oil and gas accumulation and production is the reciprocal process of oil drainage and water flooding, the minimum flow porethroat radius of oil and gas can be obtained according to the oil and gas accumulation dynamics.
The mercury intrusion curve is a curve reflecting the change of the saturation of the non-wet phase with the displacement pressure during the displacement of the wet-phase fluid by the non-wet phase fluid. The mercury intrusion process can be regarded as the accumulation and filling process of oil drainage, so the mercury intrusion curve actually reflects the correspondence between oil saturation and accumulation dynamics. In other words, the mercury inlet pressure corresponding to the original oil saturation read on the mercury pressure curve is the accumulation power, and the pore-throat radius corresponding to this accumulation power is the minimum pore-throat radius of oil and gas. Based on the determination of the minimum flow pore-throat radius, and then based on the principle of statistical analysis, the correlation curve of pore-throat radius and porosity and permeability is drawn, and the lower limit of reservoir physical properties can be calibrated according to the minimum flow pore-throat radius. Taking the Chang 6 3 reservoir in the A oilfield of the Ordos Basin as an example, the lower limit of the physical properties of the reservoir is determined by the correction method.

Geological overview
The structure of the A oil field belongs to the southwest of the Yishan slope in the Ordos Basin. It is located in Huachi and Qingyang in Gansu Province, with an area of 2600km 2 . More than 300 exploration and evaluation wells have been drilled. The main production layer of the oil field is the Chang 63 oil group of the Upper Triassic Yanchang Formation in the Triassic system, which belongs to the gravity flow sedimentation of deep lakesemi-deep lake facies, with an average thickness of 47m and a sand-to-land ratio of 0.52. The surface of the area belongs to the loess plateau landform, the terrain is undulating, the ground elevation is about 1150-1650m, and the relative height difference is about 500m. The structure of Chang 6 period is relatively simple. The overall structure is a gentle west-dipping monoclinic with an inclination angle of less than 1 degree. Analysis of core physical property data of 8288 samples from 77 cored wells shows that the porosity of Chang 63 reservoir in A oilfield is distributed between 4 and 15%, with an average porosity of 9.1%. The permeability distribution is between 0.01 and 0.8 × 10 -3 μm 2 , and the average permeability is 0.152 × 10 -3 μm 2 . The reservoir belongs to ultra-low porosity-ultra-low permeability reservoir.

Characteristics of mercury intrusion curve
Based on core observation and sample collection, five rock samples were selected for mercury intrusion experiment. The characteristics of the rock samples are shown in Table 1. It can be seen that the porosity of the rock samples is distributed between 6.86 and 13.51%, with an average of 9.97%; the permeability is distributed between 0.035 and 0.203 × 10 -3 μm 2 , with an average of 0.137 × 10 -3 μm 2 ; It can reflect the physical properties of Chang 63 low porosity and low permeability reservoir. According to the experimental data of mercury intrusion, a mercury intrusion curve (Figure 1) was drawn and the characteristics of mercury intrusion parameters of rock samples were calculated ( Table 2). It can be seen that the displacement pressure of the five rock samples is distributed between 0.78 and 2.85MPa, with an average of 1.82MPa. The median pressure is distributed between 2.61 and 14.31MPa, with an average of 8.51MPa. The maximum pore-throat radius is distributed between 0.258 and 0.943μm, with an average of 0.524μm. The median throat radius is between 0.051 and 0.282μm, with an average of 0.135μm. The maximum mercury saturation is between 83.94 and 92.13%, with an average of 87.58%. The mercury withdrawal efficiency is distributed between 26.83 and 32.6%, with an average of 30.36%. Overall, the displacement pressure and median pressure of the Chang 63 reservoir are higher, the median throat radius is lower (average 0.135μm), the reservoir throat is small, the maximum mercury saturation is high (average 87.58%), and mercury is withdrawn Low efficiency (average 30.36%).

Reservoir-forming power and minimum flow pore-throat radius
The original oil saturation is an important basis for determining the accumulation power. Since the rock samples taken are not in the original state, the original oil saturation cannot be measured. To this end, based on the logging interpretation results, the original oil saturation of the five rock samples is counted, and then combined with the mercury intrusion test curve to obtain the accumulation power and the minimum flow pore-throat radius (Table 3). It can be seen from Table 3 that the oil saturation of the five rock samples is distributed between 44.3 and 74.4%, with an average of 56.1%. The accumulation power is distributed between 8.97 and 9.31MPa, with an average of 9.18MPa. This shows that although the mercury intrusion curves of different rock samples are different in shape ( Figure 1) and the oil saturations are also very different, the accumulation power is relatively stable (all around 9MPa). If factors such as thin Chang 63 reservoir, gentle underground structure, and interpretation error of oil saturation logging are considered, the reservoir-forming power of Chang 63 reservoir in A oilfield is a certain value. For this reason, based on various reservoir-forming dynamics values, the reservoir-forming power of Chang 63 reservoir calculated by arithmetic average method is 9.18 MPa, and the corresponding minimum flow pore radius is 0.08 μm.

lower limit of reservoir physical properties
On the basis of the determination of the minimum flow pore-throat radius of oil and gas in Chang 63 oil layer, according to the principle of statistical analysis, the correlation curves of pore-throat radius and porosity and permeability are drawn (Figure 2-Figure 3). According to the minimum flow pore-throat radius, the corresponding lower limit of porosity and permeability is calculated. Figure 2 is a graph of the intersection of the median radius and porosity of the reservoir pore throats of five mercury-injected samples. It can be seen that the two are logarithmically related and the fitting relationship is as following, POR = 3.219×ln(R 50 ) + 17. 19 (1) The correlation coefficient is: R² = 0.810.  Figure 3 is the intersection graph of the median porethroat radius and permeability of five mercury-injected samples. It can be seen that the two are logarithmically related, and the fitting relationship is as following, PERM = 0.075×ln(R 50 ) + 0.306 (2) The correlation coefficient is: R² = 0.781. According to the minimum flow pore-throat radius of 0.08μm for oil and gas, the lower limit of physical properties of Chang 6 3 reservoir in A oilfield is calculated according to formula (1) and (2): porosity is 9.1%, permeability is 0.117 × 10 -3 μm 2 . If calculated according to the minimum flow pore-throat radius of 0.1μm for oil and gas, the lower limit of the physical properties of Chang 6 3 reservoir in A oilfield is: porosity 9.8%, permeability 0.133 × 10 -3 μm 2 . It can be seen from the comparison that the absolute error of the porosity is 0.7% and the relative error is 7.1% between the lower limit of the physical property of the reservoir calibrated according to 0.1 μm and the lower limit of the physical property of the reservoir calibrated by the correction method. The absolute error of permeability is 0.016 × 10 -of the reservoir physical property calibrated by it is of guiding significance.

Conclusion and understanding
The minimum flow pore-throat radius of reservoir oil and gas is not a definite value, but a function of reservoir-forming power. As the reservoir-forming power increases, the minimum flow pore-throat radius of oil and gas decreases. The traditional method of using 0.1μm as the minimum flow pore-throat radius of oil and gas lacks a precise theoretical basis and needs to be corrected. The mercury intrusion curve reflects the accumulation process of oil drainage and the corresponding relationship between oil saturation and accumulation dynamics. The mercury intrusion curve can be used to determine the reservoir accumulation dynamics and the minimum flow pore-throat radius.
The mercury intrusion test and logging interpretation oil saturation analysis of 5 samples of Chang 6 3 reservoir in A oilfield show that the reservoir forming power of Chang 6 3 reservoir is 9.18 MPa, the corresponding minimum flow pore radius of oil and gas is 0.08 μm, and the lower limit of the reservoir physical properties is: porosity 9.1%, permeability 0.117 × 10 -3 μm 2 . While the traditional method takes the minimum flow pore-throat radius of 0.1μm, the lower limit of the calibration reservoir physical properties is: porosity 9.8%, permeability 0.133 × 10 -3 μm 2 . Because the influence of reservoir forming dynamics is fully considered, the minimum flow pore-throat radius determined by the correction method is more reliable than the traditional method, and the lower limit of reservoir physical properties calibrated by this method has practical significance.