Sand Production Modelled with Darcy Fluid Flow Using Discrete Element Method

Authors: M. N. Nwodo, Y. P. Cheng, N. H. Minh

Abstract:

In the process of recovering oil in weak sandstone formations, the strength of sandstones around the wellbore is weakened due to the increase of effective stress/load from the completion activities around the cavity. The weakened and de-bonded sandstone may be eroded away by the produced fluid, which is termed sand production. It is one of the major trending subjects in the petroleum industry because of its significant negative impacts, as well as some observed positive impacts. For efficient sand management therefore, there has been need for a reliable study tool to understand the mechanism of sanding. One method of studying sand production is the use of the widely recognized Discrete Element Method (DEM), Particle Flow Code (PFC^3D) which represents sands as granular individual elements bonded together at contact points. However, there is limited knowledge of the particle-scale behavior of the weak sandstone, and the parameters that affect sanding. This paper aims to investigate the reliability of using PFC^3D and a simple Darcy flow in understanding the sand production behavior of a weak sandstone. An isotropic tri-axial test on a weak oil sandstone sample was first simulated at a confining stress of 1MPa to calibrate and validate the parallel bond models of PFC^3D using a 10m height and 10m diameter solid cylindrical model. The effect of the confining stress on the number of bonds failure was studied using this cylindrical model. With the calibrated data and sample material properties obtained from the tri-axial test, simulations without and with fluid flow were carried out to check on the effect of Darcy flow on bonds failure using the same model geometry. The fluid flow network comprised of every four particles connected with tetrahedral flow pipes with a central pore or flow domain. Parametric studies included the effects of confining stress, and fluid pressure; as well as validating flow rate – permeability relationship to verify Darcy’s fluid flow law. The effect of model size scaling on sanding was also investigated using 4m height, 2m diameter model. The parallel bond model successfully calibrated the sample’s strength of 4.4MPa, showing a sharp peak strength before strain-softening, similar to the behavior of real cemented sandstones. There seems to be an exponential increasing relationship for the bigger model, but a curvilinear shape for the smaller model. The presence of the Darcy flow induced tensile forces and increased the number of broken bonds. For the parametric studies, flow rate has a linear relationship with permeability at constant pressure head. The higher the fluid flow pressure, the higher the number of broken bonds/sanding. The DEM PFC^3D is a promising tool to studying the micromechanical behavior of cemented sandstones.

Keywords: Discrete Element Method, fluid flow, parametric study, sand production/bonds failure.

Digital Object Identifier (DOI): doi.org/10.5281/zenodo.1111616

References:

[1] H. Rahmati, et al., “Review of sand production prediction models,” Journal of Petroleum Engineering, vol. 2013, article ID 864981, pp. 1-17, Nov. 2012.
[2] V. Fattahpour, M. Moosavi, and M. Mehranpour, “An experimental investigation on the effect of rock strength and perforation size on sand production,” Journal of Petroleum Science and Engineering, vol. 2012, no. 86-87, pp. 172-189, Mar. 2012.
[3] K. Bjorlykke, Petroleum Geoscience: from Sedimentary Environments to Rock Physics. Berlin: Springer, 2010.
[4] H. H. Vaziri, E. Lemoine, I. D. Palmer, J. McLennan, and R. Islam, “How can sand production yield a several-fold increase in productivity: experimental and field data,” at SPE Annu. Technical Conf. and Exhibition, Texas, 2000, pp. 1-4.
[5] P. G. Ranjith, M. S. A. Perera, W. K. G. Perera, S. K. Choi, and E. Yasar, “Sand production during the extrusion of hydrocarbons from geological formations: a review,” Journal of Petroleum Science and Engineering, vol. 124, pp. 72-82, Oct. 2014.
[6] H. H. Vaziri, E. M. Lemoine, and Y. Xiao, “Quantification of sand production induced improvement in productivity index,” Canadian Geotechnical J., vol. 39, pp. 1088-1102, Sept. 2002.
[7] A. Nouri, H. Vaziri, E. Kuru, and R. Islam, “A comparison of two sanding criteria in physical and numerical modelling of sand production,” Journal of Petroleum Science and Engineering, vol. 2006, no. 50, pp. 55-70, Oct. 2005.
[8] Y. Wang, and E. Papamichos, “Sand prediction by different criteria and validation through a hollow cylinder test,” at the 46th U.S. Rock Mechanics/Geomechanics Symposium, Chicago, 2012, pp. 24-27.
[9] Y. H. Wang, and S. C. Leung, “A particulate-scale investigation of cemented sand behavior,” Canadian Geotechnical J., vol. 45, pp. 29-44, Feb. 2008.
[10] N. Climent, M. Arroyo, C. O’Sullivan, and A. Gens, “Sand production simulation coupling DEM with CFD,” European Journal of Environmental and Civil Engineering, vol. 18, no. 9, pp. 983-1008, May, 2014.
[11] Z. Y. Zhou, A. B. Yu, and S. K. Choi, “Numerical simulation of the liquid-induced erosion in a weakly bonded sand assembly,” Powder Technology, vol. 211, pp. 237-249, Apr. 2011.
[12] B. Ikporo, and O. Sylvester, “Effect of sand invasion on oil well production: a case study of Garon field in the Niger Delta,” The International Journal of Engineering and Science, vol. 4, no. 5, pp. 64-72, 2015.
[13] P. A. Cundall, “A computer model for simulating progressive large scale movement in blocky rock systems,” in Proceedings of the Symposium of International Society of Rock Mechanics, France, vol. 1, paper no. II-8, 1971.
[14] P. A. Cundall, “Formulation of a three-dimensional distinct element model – Part I. A scheme to detect and represent contacts in a system composed of many polyhedral blocks,” International Journal of Rock Mechanics, Mineral Science and Geomechanics, vol. 25, no. 3, pp. 107-116, 1988.
[15] R. Hart, P. A. Cundall, and J. Lemos, “Formulation of a three-dimensional distinct element model – Part II. Mechanical calculations for motion and interaction of a system composed of many polyhedral blocks,” International Journal of Rock Mechanics, Mineral Science and Geomechanics, vol. 25, no. 3, pp. 117-125, 1988.
[16] R. M. O’Connor, J. R. Torczynski, D. S. Preece, J. T. Klosek, and J. R. Williams, “Discrete element modelling of sand production,” International Journal of Rock Mechanics and Mining Sciences, vol. 34, no. 231, pp. 3-4, 1997.
[17] L. Li, and R. M. Holt, “Particle scale reservoir mechanics,” Oil and Gas Science and Technology, vol. 57, no. 5, pp. 525-538, 2002.
[18] L. Li, E. Papamichos, and P. Cerasi, “Investigation of sand production mechanisms using DEM with fluid flow,” Multiphysics Coupling and Long Term Behavior in Rock Mechanics, London, 2006.
[19] D. O. Potyondy, “Simulating stress corrosion with a bonded-particle model for rock,” International Journal of Rock Mechanics and Mining Sciences, vol. 44, pp. 677-691, 2006.
[20] E. Papamichos, “Sand production physical and experimental evidence,” at the French 18th ALERT School-Geomechanical and Structural Issues in Energy Production, Aussois, 2006.
[21] S. M. Haeri, A. Seiphoori, and A. Rahmati, “The behavior of a limy cemented gravely sand under static loading-case study of Tehran alluvium,” Electronic Journal of Geotechnical Engineering, vol. 13, pp. 1-10, 2008.
[22] L. Y. G. Cheung, C. O’Sullivan, and M. R. Coop, “Discrete element method simulations of analogue reservoir sandstones,” International Journal of Rock Mechanics and Mining Sciences, vol. 63, pp. 93-103, July, 2013.