In random network models, sizes for pores and throats are distributed according to a truncated Weibull distribution. As a result, parameters defining the shape of the distribution are critical for the characteristic of the network. In this paper, an algorithm to distribute pores and throats in random network was established to more representatively describe the topology of porous media. First, relations between Weibull parameters and the distribution of dimensionless throat siz es were studied and a series of standard curves were obtained. Then, by analyzing the capillary pressure curve of the core sample, frequency distribution histogram of throat sizes was obtained. All the sizes were transformed to dimensionless numbers ranged from 0 to 1. Curves of the core were compared to the standard curves, and truncated Weibull parameters could be determined according an inverse algorithm. Finally, aspect ratio and average length of throats were adjust ed to simultaneously fit the porosity and the capillary pressure curves and the whole network was established. The predicted relative permeability cu rves were in good agreement with the experimental data of cores, indicating the validity of the algorithm.
Network models have been widely used to study fluid flow in porous media. There are two general kinds of networks: networks that representatively describe the topology of real core samples and random networks which only capture a few characteristic of porous media. The former ones can be used to accurately predict experimental data such as relative permeabilities and oil recovery (Øren et al., 1998), but constructing this kind of networks is costly. By contrast, random networks are founded with very low cost and mainly used to investigate non-Newtonian flow (Lopez et al. 2003; Taha, 2007) and microscopic mechanism of enhanced recovery of oil such as polymer flooding (Sun et al., 2011) or surfactant flooding (Bo et al., 2003, SPE84906). What’s more, wetting condition plays a significant role in determining transport properties in porous media (Romero-Zeron, 2009; Rostami et al., 2010; Kasin, 2011) and random networks can qualitatively describe wettability effects (Zhao et al., 2010) if wetting conditions can be properly or quantitatively expressed (Øren et al., 2003; Feng et al., 2012) and random networks are representative. So, parameter optimization to construct valid networks is of great importance.
In network models, pores and throats are approximated as cylindrical ducts of arbitrary shapes and different sizes of cross-sections (Patzek, 2000). This is suitable for studying seepage characteristics under the influence of different parameters through stochastic network model. However, no method has been studied and proposed for assigning pore throat parameters of the network model so that the constructed model could represent the pore throat parameters and seepage characteristics of the actual core for specific oilfields.
In this paper, a new algorithm is established to calculate the truncated Weibull parameters which determine the sizes of throats and play critical roles in determining the shape of capillary pressure curves and relative permeability curves.
Sizes of throats are distributed according to truncated Weibull function, which can be written as:
r = ( r max − r min ) ( − δ ln [ x ( 1 − e − 1 / δ ) + e − 1 / δ ] ) 1 / γ + r min (1)
In the formula, r max and r min are the inscribed radius of the biggest and smallest throats respectively. x is a random number ranged from 0 to 1. δ and γ are Weibull parameters.
Assume that there are (n + 1) throats in random network. If network is sufficiently big, x in Equation (1) are statistically evenly valued as 0 n ， 1 n ， 2 n , and so on with the biggest value n n . Sizes of the throats are transformed to dimensionless numbers between 0 and 1:
r D = r r max (2)
As a result, r D max = r max r max = 1 and r D min = r min r max ≈ 0 . Number the throats in
order of size. Using Equation (1) and Equation (2), the dimensionless radius of all throats are
{ r D 0 = ( 1 − r D min ) ( − δ ln [ 0 n ( 1 − e − 1 / δ ) + e − 1 / δ ] ) 1 / γ + r D min r D 1 = ( 1 − r D min ) ( − δ ln [ 1 n ( 1 − e − 1 / δ ) + e − 1 / δ ] ) 1 / γ + r D min ⋮ r D n = ( 1 − r D min ) ( − δ ln [ n n ( 1 − e − 1 / δ ) + e − 1 / δ ] ) 1 / γ + r D min (3)
The dimensionless average throat size can be calculated by
R ¯ D = ∑ i = 0 n r D i n + 1 = 1 − r D min n + 1 ∑ i = 0 n ( − δ ln [ i n ( 1 − e − 1 / δ ) + e − 1 / δ ] ) 1 / γ + r D min (4)
Relationship of R ¯ D versus Weibull parameters are depicted in semi logarithmic coordinates, as shown in
In
Classify the range from 0 to 1 evenly into ten small intervals: 0 to 0.1, 0.1 to 0.2 and so on. Values between 0 and 0.1 are all considered as the value of 0.05, et al., for reasons of convenience. Calculate all dimensionless throat sizes with Equation (3) and distribution frequency are obtained by adding up those dimensionless throat sizes which are valued in the same small intervals.
The frequency distribution curves under different γ are shown in
with δ fixed as 0.2 and in
Peak dimensionless throat sizes, r D , p e a k , are the dimensionless throat sizes that correspond to the biggest frequency. From
That is to say, more throats are assigned small values in these conditions.
Capillary pressure curves of the rock core sample are acquired by mercury intrusion method and porosity are calculated after measurement of void space and bone volume of the core in laboratory. According to the capillary pressure curves, frequency distribution histograms of throat sizes are obtained. P S min , the pressure corresponds to the minimum wetting phase saturation and P T , the threshold capillary pressure can be read from the mercury intrusion curves (Qin & Li, 2006).
The maximum and minimum throat size, r max and r min , can be calculated by replacing P c in Equation (5) with P T and P S min respectively:
r = 2 σ cos θ P c (5)
In which σ is the surface tension of the mercury-air system used in the core experiment and θ is mercury-rock contact angle.
Transform the frequency distribution histogram of throat sizes into frequency distribution curve of dimensionless throat sizes with Equation (2) through which peak dimensionless throat size can be read. Dimensionless average throat size R ¯ D can be obtained by analysis of frequency distribution histogram (Qin & Li, 2006) and calculation with Equation (2). The following steps are carried out to determine the parameters γ and δ ,and the flowchart is shown in
1) γ is assigned a random initial value between 1 and 10.
2) Calculate δ iteratively using Equation (4) or by
3) With current δ and γ , the standard frequency distribution curve are determined (like those in
4) Repeat (2) and (3) until the relative error between r D , c p and r D , p e a k are allowable.
5) δ and γ are determined and throat size distribution of the random network is obtained together with the upper and lower limits of the throat size.
In networks, aspect ratio, the ratio between the pore radius and connecting throat radius, is used to distribute pore sizes, which can be expressed as
α = r p o r e r t h r o a t (6)
According to spatial correlation principle, size of pore is larger than that of all its connecting throats. Known the throat distribution and aspect ratio (Lu et al., 1997; Sun, 2002), pore size can be obtained by
r p = max ( α ∑ i = 1 n C r i n C , max ( r i ) ) (7)
in which n c is the numbers of connecting throats of the pore.
In random networks, pores are distributed in regular node. Throat length is the distance between its connecting pores minus the radius of the two pores. So, the grid of networks determines the length of throats and further the porosity. If the grid lines are denser, pores are closer to each other and the porosity of the network is larger while the grid lines are looser, pores are farther away from each other, resulting in long throats, more bone volume and low porosity.
Core sample is selected. Experimental and treated data are listed in
| Porosity | 0.387 |
|---|---|
| Throat size | 1.93 μm - 8.96 μm |
| Oil-water tension in waterflooding | 30 mN/m |
| Average coordinate number | 4.27 |
According to the algorithm, δ = 0.18 and γ = 2.7 . Adjust the density of grid lines to fit the porosity with fine tuning coordinate numbers to simultaneously fit the capillary curves (dotted lines in
Simulate water flooding with Valvatne’s modeling code. Relative permeability curves with the random network are dotted in
An algorithm is established to determine the truncated Weibull parameters for distribution of throats and pores in random network models using this method, founded networks can more accurately and representatively describe the topology of rock cores.
In random network models, the two Weibull parameters determine the sizes of throats and play critical roles in determining the shape of capillary pressure curves and relative permeability curves. Knowing the dimensionless peak throat size and average throat size, Weibull parameters can be determined with iterations. Density of grid lines in random network has great impact on porosity. The denser the grid lines are, the larger porosity of the network is.
In future, pore throat parameters assigning methods based on artificial intelligence algorithms might be an interesting direction for further research.
The authors declare no conflicts of interest regarding the publication of this paper.
Shi, F., Wu, J. D., Chang, T., Sun, C. J., & Wu, X. Z. (2022). An Algorithm to Determine the Truncated Weibull Parameters for Distribution of Throats and Pores in Random Network Models. Journal of Geoscience and Environment Protection, 10, 46-53. https://doi.org/10.4236/gep.2022.1010004