Research Article  Open Access
Numerical Modeling of the Conductivity of the Particle Monolayer with Reduced Size
Abstract
Fractures filled with a proppant monolayer play an important role in the hydraulic fracture network. Predicting the conductivity of these fractures is the basis of fracture network optimization. However, little attention has been paid to the conductivity of the proppant monolayer. The change of conductivity under various conditions is currently not fully understood. Therefore, in this paper, the conductivity variation under different conditions are simulated. The reduction of particle size was calculated by existing analytical models. The permeability variation was calculated through computational fluid dynamics (CFD) combined with COMSOL Multiphysics. The controlling factors of conductivity under a proppant monolayer were identified. Simulation results indicate that elastic parameters, closure pressure, and proppant distribution have significant influence on conductivity, while creep parameters, such as rock viscosity and time, have limited influence on conductivity. Moreover, the changes in permeability, porosity, and tortuosity with variation of embedment were analyzed. Results indicated that with an increase in embedment, the permeability and porosity decrease as expected. The main reduction (nearly half) emerges in the first 20% of proppant embedment. Furthermore, the permeability of a single particle deviates largely from the prediction of CarmanKozeny (CK) equation. The tortuosity of proppant particle increases with a decrease in particle size due to embedment. A modification of the CarmanKozeny equation is proposed to address this influence.
1. Introduction
The reserve of oil and gas in conventional and unconventional resources is abundant [1, 2]. Unconventional oil and gas reservoirs have become the focus of exploration and exploitation in petroleum industry. Horizontal drilling and multistage fracturing are widely used in the exploitation of these reservoirs [3]. Accurately measuring and understanding these technologies are critical for ensuring successful reservoir exploitation [4]. For hydraulic fracturing, it consists of initiating, propagating, and maintaining a fracture from the wellbore to pay zones. Proppants are pumped into the created fracture along with the fracturing fluid. When the hydraulic pressure is removed from a well, the surfaces of fractures compress onto the proppants, creating a highpermeability pathway in deep rock formations through which natural gas and oil flow more freely.
The performance of fractured wells is controlled by two types of parameters. One type of parameters is formation parameters, such as porosity, permeability, and mechanical properties [5]. The second type of parameters is fracture parameters, such as fracture length and fracture conductivity [6, 7]. Fracture conductivity reflects the transport capacity of the high permeable channel, which limits the amount of stimulation achieved from the fracturing treatments of many wells [8]. It is defined as the product of fracture width and the apparent permeability of the proppant pack. Results summarized from over 80 field studies show the benefit of increasing fracture conductivity [9]. As lowquality reservoirs become the increasingly popular focus, increasing the conductivity of hydraulic fracture is one of the possible strategies in the future to improve the fracturing response [10]. Any measures that intended to increase the fracture conductivity will increase the cost. The International Energy Agency (IEA) indicates that the new worldwide oil and gas infrastructure is projected to cost US $20 trillion between 2011 and 2035 [11]. However, no matter whether the project is onshore or offshore, many oil and gas projects experienced significant cost overruns [12]. Under the condition of low oil and gas prices, improving the evaluation efficiency of oil and gas reservoir quality and optimizing the technical parameters are two effective ways to cut the cost. Therefore, it is worth optimizing the fracture conductivity in order to reduce the costs and improve postfrac performance.
Slick water fracturing is widely used in the development of tight reservoirs [13]. A complex network of narrow secondary fractures is generated to extend the contact area between the network of flow channels and formation. Some of these narrow fractures are only propped by a proppant monolayer [14]. Initially, the industry believes that a monolayer is difficult to achieve. However, the situation changes in unconventional reservoirs; many hydraulically activated shear fractures propagate deeply in the reservoir. These fractures could only open at a limited width and would only be poorly propped by a monolayer of proppants. Therefore, the conductivity of the proppant monolayer is an important factor in the design and optimization of hydraulic fracturing.
But little attention has been paid in previous studies. The dominant factors that affect the conductivity of the proppant monolayer are not clear. Besides, the prediction of fracture conductivity is the foundation of fracture optimization. In order to predict the fracture conductivity with the proppant monolayer, a good understanding of the changes of fracture width and permeability of the proppant pack is required. Proppant embedment is the main cause of fracture width reduction, which can be solved by the models available in the literatures. On the calculation of permeability of proppant pack, CarmanKozeny (CK) equation was widely used [15–18]. The semiempirical CarmanKozeny (CK) equation is the most famous model that describes the permeabilityporosity relation. However, this relation has many limitations. It assumes that the porous medium is a packing of identical spheres and the tortuosity is an empirical constant. Then, when there is only a proppant monolayer, whether the equation is still applicable is unclear. Furthermore, after proppant embedding into rock, the geometry of the proppant is no longer a circular shape, as shown in Figure 1. How will this affect the permeability and its deviation from the prediction of CK equation?
Therefore, in this paper, the conductivity of the proppant monolayer was studied. First, computational fluid dynamic theory is used to study the influencing factors controlling conductivity under a proppant monolayer. Then, the permeability of the proppant monolayer with different levels of embedment was studied. A modification of CK equation is proposed based on numerical results. The models and results can provide an insight into the features and factors controlling the conductivity of the proppant monolayer.
2. Numerical Methodology
A computational fluid dynamics (CFD) approach was utilized to model the conductivity reduction due to proppant embedment. A series of numerical simulation were performed to study the variation of fracture conductivity with the proppant monolayer. CFD has been widely used in studying fluid flow [14, 17, 19, 20]. In this work, COMSOL Multiphysics was used. If the structure of the medium is determined, CFD can solve the steadystate Navier–Stokes equations to obtain the velocity field and then compute the permeability. Incompressible fluid and laminar flow model were used, as shown in Table 1.

Actually, CFD simulation of fluid flow inside a field scale geometry will require large computational efforts, which exceeds the capacity of most computers. In order to improve the computational efficiency, three computational domains were used as the physical model, as shown in Figure 2. One is a simple regular lattice called facecentered cubic (FCC) packing. Fluid flow in a FCC packing was demonstrated to be consistent with the CarmanKozeny equation [21]. Therefore, it was used to verify the model and numerical solution of COMSOL Multiphysics. The second one is a single particle, which is the focus of this paper. The third one includes multiple particles, which are used to study the influence of proppant distribution on conductivity.
(a) FCC packing
(b) Single particle
(c) Multiple particles
The CarmanKozeny equation is a wellknown formula and a relatively simple way to predict permeability. The permeability can then be expressed as where is the permeability, m^{2}; is the porosity, dimensionless; is the tortuosity, dimensionless; and is the specific surface area per pore volume, 1/m. If the porous medium is a packing of identical spheres, then is given by where is the diameter of particles, m. in equation (2) does not consider the influence of embedment. In later calculation, the definition of is used to calculate its value. Equation (2) can be substituted into equation (1) to yield
For random packing of spheres with the same radius, is commonly used and results in the following wellknown equation:
For the numerical simulations, it is assumed that the flow was threedimensional, singlephase, and laminar. Each domain also contains an entrance and an exit on each side. A constant pressure boundary was used at the inlet. A constant zero static pressure was used at the outlet. Wall boundaries were set on the other surfaces. Grid size affects the accuracy of the numerical simulation. We conducted a series of calculations with different grid sizes to check the computational stability and accuracy. Each set of comparison processes (data in the same figure) used the same mesh to ensure the convergence and computational speed. Simulation results indicate that the numerical permeability obtained by this mesh is 919 μm^{2}, the analytical permeability calculated by CK equation is 928 μm^{2}, and the difference is 0.9%.
3. Numerical Results and Discussions
In the reference [15], analytical models were derived to calculate proppant embedment, proppant deformation, change in fracture aperture, and fracture conductivity. And then, this model has been improved by Zhang and Hou [16] with consideration of the rock creep deformation. Here, we would like to give a brief introduction of their models.
Proppant embedment, proppant deformation, and the change of fracture width were derived as follows:
Then, considering creep deformation, these parameters change to where is the deformation of the proppant, mm; is the diameter of the proppant, mm; is the thickness of rock, mm; is the distance coefficient; is the closure pressure, MPa; is Poisson’s ratio of the proppant, mm; is Poisson’s ratio of rock, mm; is the elastic modulus of the proppant, MPa; is the elastic modulus of rock, MPa; is the value of embedment, mm; is the change in fracture aperture, mm; is the proppant viscosity, mPa·s; and is the rock viscosity, mPa·s.
This paper focuses on the effect of different factors on the conductivity of the proppant monolayer. Therefore, analytical models [15, 16] were used to calculate the embedment depth and fracture residual width under different conditions. Then, fracture conductivity was evaluated by a CFD method.
3.1. Effect of Elastic Parameters
To investigate the effect of elastic parameters on fracture conductivity, a normalized parameter can be defined as follows: where is Young’s modulus (1 is proppant; 2 is rock), Pa. Simulation tests were carried out with six different closure pressure values under different .
Figure 3 shows the relationship between and normalized conductivity at different stress conditions. The conductivity after proppant embedment and deformation is normalized by the conductivity of fracture with a width that equals the proppant diameter (0.85 mm). As shown in the figure, the conductivity decreases when increases for all the cases. With an increase in the proppant elastic modulus or a decrease in the rock elastic modulus, the fracture conductivity shows an approximately linear decrease. When the closure pressure is 10 MPa, the conductivity decrease 6.66% with an increase in . However, when the closure is 60 MPa, the conductivity decrease 26.91% with an increase in . Under high closure pressure, the conductivity is more sensitive to . For some reservoirs, such as shale, the rock elastic modulus will decline when encountered with fracturing fluid [22]. If the proppant with high elastic modulus is used, the will become higher too, which would undoubtedly result in sharp decrease in fracture conductivity.
3.2. Effect of Creep Parameters
Figures 4 and 5 show the relationship among rock viscosity, closure pressure, and normalized conductivity. As the rock viscosity increases, the conductivity increases. It is found that the conductivity increases relatively obvious when the rock viscosity increases from 1000 mPa·s to 4000 mPa·s and gradually levels off when the rock viscosity is above 4000 mPa·s.
As shown in the figure, the conductivity increases by 4.22% as the rock viscosity increases from 1000 mPa·s to 4000 mPa·s under 10 MPa closure pressure and by 16.62% when the closure pressure is 60 MPa. In the case of low closure pressure, the increase of rock viscosity does not make a big difference to the increment of conductivity.
However, the influence of closure pressure on conductivity is significant. As shown in Figures 4 and 5, the conductivity decreases by 44.1% as the closure pressure increases from 10 MPa to 60 MPa with 1000 mPa·s rock viscosity and by 35.78% when the rock viscosity is 10000 mPa·s.
Besides, as shown in Figures 6 and 7, the conductivity decreases by 11.93% as the time increases from 100 d to 600 d at rock viscosity of 1000 mPa·s and by 8.07% when the rock viscosity is 10000 mPa·s. The time does not have a significant influence on conductivity.
3.3. Effect of Diameter and Distribution of the Proppant
Figures 8 and 9 show the relationship among proppant diameter, closure pressure, and normalized conductivity. As the proppant diameter increases, the conductivity increases sharply. The relationship between the particle size and conductivity satisfies the power function.
As shown in Figures 8 and 9, the conductivity increases by 49.00% as the closure pressure increases from 10 MPa to 60 MPa with a proppant diameter of 0.2 mm and by 28.87% when the proppant diameter is 0.85 mm.
Several proppant distribution cases were simulated as shown in Figure 10. When considering proppant embedment, the changes of physical models are illustrated in Figure 11. Simulation results are shown in Figure 12. Embedment ratio is defined as the fraction of the proppant grain embedment to the initial fracture width (proppant diameter). As the proppant distribution changes, the conductivity varies significantly. The patterns of A and E are better than the other ones. The results demonstrate that the discontinuous proppant placement is better than continuous placement under different embedment ratios.
(a)
(b)
(c)
(d)
(e)
3.4. Modification of CarmanKozeny Equation
CarmanKozeny equation can be used to calculate the conductivity as evidenced in many analytical models, such as reference [15]. However, as described before, the CarmanKozeny equation does not consider the particle size change resulted from embedment. A modification of the CK equation is proposed based on numerical results.
Three particle sizes were simulated, ranging from 212 μm to 840 μm (70 mesh to 20 mesh). The results are shown in Figure 13. The permeability is normalized against the permeability of a channel with the maximum initial fracture width, which is the initial permeability of 20mesh proppantsupported fracture. The permeability can be affected by the proppant size and the embedment ratio. With an increase in proppant size and a decrease in the embedment ratio, the permeability increases. For the proppant monolayer, the permeability reduces to its half when the embedment ratio increases to nearly 20%. Proppant has two embedment mechanisms, elastic embedment and creep embedment [16]. The fracture will keep closing under creep deformation.
Furthermore, results show that with an increase in the embedment ratio, the porosity decreases and the specific surface area remains constant, as shown in Figure 14. When the porosity decreases to a certain value, the decreased amplitude reduces. After the proppant embedment increases to 90%, the porosity reduces to 20%. From equation (1), the decreased porosity will also lead to a decrease in permeability. However, the prediction of the CK equation is still larger than the numerical results.
Therefore, tortuosity induces the difference. Previous researches showed that tortuosity will change with porosity [23]. Under each different situation, they meet a specific relationship. Four equations are used to fit the relationship between porosity and tortuosity. Their results are shown in Figure 15.
Results are summarized in Table 2. It can be seen that all the equations’ calculations fit the numerical results well with a . Equation (11) achieves the best fitting with a . However, when these results were used in the permeability calculation, only equation (11) matches well with the results, as shown in Figure 16.
4. Conclusions
We have conducted a numerical study on the fracture conductivity change caused by proppant embedment. Simulation results enhance the understanding of the conductivity variation under proppant monolayer condition. The following conclusions can be summarized according to the present study: (1)Simulation results indicate that elastic parameters, closure pressure, and proppant distribution have significant influence on conductivity, while creep parameters, such as rock viscosity and time, have limited influence on conductivity(2)The numerical results showed that the first 20% of proppant embedment can cause a significant permeability reduction(3)For a single proppant particle, with an increase in the embedment ratio, the porosity decreases and the specific surface area keeps constant. After the proppant embedment increases to 90%, there is still 20% porosity left. Meanwhile, tortuosity of the proppant monolayer increases with the increase in embedment(4)The numerical simulation was able to capture the flow characteristics around the proppant. The permeability of a single particle deviates largely from the prediction of the CK equation. A new correlation between porosity and tortuosity was proposed to solve this deviation(5)The placement of the proppant is very complex, which brings great difficulty to the prediction of the fracture conductivity. One should note that in a real fracture there is no perfect placement with uniform size throughout a complex fracture network. The models in this paper focus on the base case of the proppant monolayer. Any displacement of the proppant monolayer could be a combination of single particle cases. However, the models consider single phase, laminar flow, single particle, and their combination. Any deviation from the assumptions can make the prediction of the models different from real cases. It is recommended to develop a more comprehensive model by the method given in this paper, such as multilayers and twophase flow
Nomenclature
:  Constant 
:  Constant 
:  Constant 
:  Diameter of particles, m 
:  Diameter of the proppant, mm 
:  Thickness of rock, mm 
:  Constant 
:  Elastic modulus of the proppant, MPa 
:  Elastic modulus of rock, MPa 
:  Dimensionless Young’s modulus 
:  Value of embedment, mm 
:  Permeability, m^{2} 
:  Distance coefficient 
:  Closure pressure, MPa 
:  Specific surface area, 1/m 
:  Tortuosity, dimensionless 
:  Porosity, dimensionless 
:  Change in fracture aperture, mm 
:  Deformation of the proppant, mm 
:  Poisson’s ratio of the proppant 
:  Poisson’s ratio of rock 
:  Proppant viscosity, mPa·s 
:  Rock viscosity, mPa·s. 
Conflicts of Interest
The authors declare that there is no conflict of interest regarding the publication of this paper.
Acknowledgments
The authors would like to thank the financial support of National Natural Science Foundation of China (nos. 51525404, 51374178, and 51604232), the Young Scholars Development Fund of SWPU (201599010084), and the National Science and Technology Major Project (no. 2016ZX05002002).
References
 C. Zou, G. Zhai, G. Zhang et al., “Formation, distribution, potential and prediction of global conventional and unconventional hydrocarbon resources,” Petroleum Exploration and Development, vol. 42, no. 1, pp. 14–28, 2015. View at: Publisher Site  Google Scholar
 C. McGlade, J. Speirs, and S. Sorrell, “Unconventional gas – a review of regional and global resource estimates,” Energy, vol. 55, pp. 571–584, 2013. View at: Publisher Site  Google Scholar
 N. R. Warpinski, M. J. Mayerhofer, M. C. Vincent, C. L. Cipolla, and E. P. Lolon, “Stimulating unconventional reservoirs: maximizing network growth while optimizing fracture conductivity,” Journal of Canadian Petroleum Technology, vol. 48, no. 10, pp. 39–51, 2009. View at: Publisher Site  Google Scholar
 G. E. King, “Hydraulic fracturing 101: what every representative, environmentalist, regulator, reporter, investor, university researcher, neighbor and engineer should know about estimating frac risk and improving frac performance in unconventional gas and oil wells,” in SPE Hydraulic Fracturing Technology Conference, p. 80, The Woodlands, TX, USA, 2012. View at: Publisher Site  Google Scholar
 Y. Shen, M. Meng, T. Liu, H. Ge, and Y. Zhang, “Impact of petrophysical properties on hydraulic fracturing and development in tight volcanic gas reservoirs,” Geofluids, vol. 2017, Article ID 5235140, 13 pages, 2017. View at: Publisher Site  Google Scholar
 C. Dai, L. Xue, W. Wang, and X. Li, “Analysis of the influencing factors on the well performance in shale gas reservoir,” Geofluids, vol. 2017, Article ID 7818346, 12 pages, 2017. View at: Publisher Site  Google Scholar
 J. Wang, X. Wang, and W. Dong, “A semianalytical model for multiplefractured horizontal wells with SRV in tight oil reservoirs,” Geofluids, vol. 2017, Article ID 2632896, 15 pages, 2017. View at: Publisher Site  Google Scholar
 C. E. Cooke Jr., “Conductivity of fracture proppants in multiple layers,” Journal of Petroleum Technology, vol. 25, no. 9, pp. 1101–1107, 1973. View at: Publisher Site  Google Scholar
 M. C. Vincent, “Proving it  a review of 80 published field studies demonstrating the importance of increased fracture conductivity,” in SPE Annual Technical Conference and Exhibition, p. 21, San Antonio, TX, USA, 2002. View at: Google Scholar
 M. C. Vincent, “The next opportunity to improve hydraulicfracture stimulation,” Journal of Petroleum Technology, vol. 64, no. 3, pp. 118–127, 2012. View at: Publisher Site  Google Scholar
 IEA, World Energy Outlook, IEA, Paris, France, 2011.
 O. J. Olaniran, P. E. D. Love, D. Edwards, O. A. Olatunji, and J. Matthews, “Cost overruns in hydrocarbon megaprojects: a critical review and implications for research,” Project Management Journal, vol. 46, no. 6, pp. 126–138, 2015. View at: Publisher Site  Google Scholar
 T. T. Palisch, M. Vincent, and P. J. Handren, “Slickwater fracturing: food for thought,” SPE Production & Operations, vol. 25, no. 3, pp. 327–344, 2010. View at: Publisher Site  Google Scholar
 A. Khanna, A. Kotousov, J. Sobey, and P. Weller, “Conductivity of narrow fractures filled with a proppant monolayer,” Journal of Petroleum Science and Engineering, vol. 100, pp. 9–13, 2012. View at: Publisher Site  Google Scholar
 K. Li, Y. Gao, Y. Lyu, and M. Wang, “New mathematical models for calculating proppant embedment and fracture conductivity,” SPE Journal, vol. 20, no. 3, pp. 496–507, 2015. View at: Publisher Site  Google Scholar
 J. Zhang and J. Hou, “Theoretical conductivity analysis of surface modification agent treated proppant II – channel fracturing application,” Fuel, vol. 165, pp. 28–32, 2016. View at: Publisher Site  Google Scholar
 P. Sanematsu, Y. Shen, K. Thompson et al., “Imagebased Stokes flow modeling in bulk proppant packs and propped fractures under high loading stresses,” Journal of Petroleum Science and Engineering, vol. 135, pp. 391–402, 2015. View at: Publisher Site  Google Scholar
 M. M. M. Shamsi, S. F. Nia, and K. Jessen, “Dynamic conductivity of proppantfilled fractures,” Journal of Petroleum Science and Engineering, vol. 151, pp. 183–193, 2017. View at: Publisher Site  Google Scholar
 A. Trykozko, M. Peszynska, and M. Dohnalik, “Modeling nonDarcy flows in realistic porescale proppant geometries,” Computers and Geotechnics, vol. 71, pp. 352–360, 2016. View at: Publisher Site  Google Scholar
 J. Zhang, L. Ouyang, D. Zhu, and A. D. Hill, “Experimental and numerical studies of reduced fracture conductivity due to proppant embedment in the shale reservoir,” Journal of Petroleum Science and Engineering, vol. 130, pp. 37–45, 2015. View at: Publisher Site  Google Scholar
 A. A. Zick and G. M. Homsy, “Stokes flow through periodic arrays of spheres,” Journal of Fluid Mechanics, vol. 115, no. 1, pp. 13–26, 1982. View at: Publisher Site  Google Scholar
 O. M. Akrad, L. M. Jennifer, and P. Manika, “The effects of fracturing fluids on shale rock mechanical properties and proppant embedment,” in SPE Annual Technical Conference and Exhibition, p. 12, Denver, CO, USA, 2011. View at: Publisher Site  Google Scholar
 B. Ghanbarian, A. G. Hunt, R. P. Ewing, and M. Sahimi, “Tortuosity in porous media: a critical review,” Soil Science Society of America Journal, vol. 77, no. 5, pp. 1461–1477, 2013. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2018 Yuxuan Liu et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.