Analytical Modeling of Multistage Hydraulically Fractured Horizontal Wells Producing in Multilayered Reservoirs with Inter-Layer Pure-Planar Crossflow Using Source/Sink Function Method

This study presents a comprehensive analytical modeling technology to model transient behaviors of multilayered reservoirs with inter-layer pure-planar crossflow induced by multi-stage hydraulically fractured horizontal well (MHFHW). The objective of this study is to develop an analytical model for multilayered reservoirs in conjunction with complex MHFHW and to achieve not only accurate and efficient computation, but also well-organized solutions expressed in a systematically integrated manner. The consideration of inter-layer crossflow across adjacent layers sets up the foundation for successful modeling of multilayered reservoirs. Source/sink function method (SSFM) is applied to describe fluid flow. Unsteady-state pressure or production rate solutions of MHFHW with the advantages of fast computation, accurate, and stable solutions are achieved. Comparative and consistent outcomes generated by this work and widely applied industry software have largely enhanced our technical confidence. More importantly, innovatively defined modified dimensionless terms that integrate systematic well-reservoir geometry information, as well as rock/fluid properties of each layer, have been newly applied to regulate the new modified dimensionless rate decline curve. This new technique sheds light on the reservoir characterization practice for complicated reservoir systems. Theoretical results in terms of transient pressure and rate were generated by the proposed multilayered model (SSFM-ML) for five scenarios of general concern, under various reservoir and well parameters, which were examined and discussed to demonstrate technical robustness. Not only does this study give solutions to the targeted multiple layered reservoirs, but it also provides insights into modeling three-dimensional fluid flow in heterogeneous reservoir with complex well configurations. It is recommended that future research should be conducted for more complicated two- and three-dimensional reservoirs, using the similar strategy of developing new type curves through adopting other new forms of modified dimensionless rate and time terms.