Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)Wellbore flow and its coupling with porous media has been an active field of research, mainly because wellbore hydraulics (especially in complex geometry wells, like horizontal or multilateral) plays an important role in the pressure and production profiles. An inherent difficulty is associated to the multiscale nature of the flow problem: the cell hosting the well is generally much larger than wellbore radius. As a consequence, average pressure in a cell is not a good approximation of the well pressure. Moreover, at the interface between well and reservoir there is a change in flow pattern, from porous media to pipe flow. In this work a well-reservoir coupling technique is developed and applied to the numerical simulation of horizontal wells in gas reservoirs. This model considers wellbore hydraulics and storage. We consider the 3D single-phase isothermal flow of a compressible fluid in porous media and ID single-phase isothermal flow of a compressible fluid in the wellbore. In both regions, mass conservation, momentum balance and real gas equation of state are the governing equations. Conservative forms are used for continuity and momentum. Flow continuity and pressure equilibrium are considered at the interface between well and reservoir. The nonlinear partial differential equations can be rewritten to allow the coupling of the pressure field between the porous media (reservoir) and the free media (well) by means of an implicit formulation for the finite difference method presented. A hybrid grid technique was applied to improve the description of the flow close to the wellbore. In addition, pressure in reservoir and wellbore are solved simultaneously, while velocities in wellbore are updated at every iterative step. An iterative method to solve numerically the resulting linear system for the primitive variable pressure was developed and it includes the wellbore/reservoir coupling into a finite difference numerical model. This methodology is useful for the analysis of transient pressure and velocity behavior in both reservoir and wellbore including the change in flow pattern from porous media to pipe flow. Effects of permeability, gas specific gravity, formation damage and grid refinement are considered in different sensitivity analysis for several well-reservoir systems. Results for pressure and pressure derivative followed the expected behavior of the reservoir model chosen, both in early and long times. Effects of friction, convection and compressibility in wellbore hydraulics were successfully included in the numerical approach. Copyright 2014, Society of Petroleum Engineers.2 137213972011/11897-6; FAPESP; São Paulo Research Foundation; 2012/19874-8; FAPESP; São Paulo Research FoundationFundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)Aguilar, C., Ozkan, E., Et Al., K.H., Transient behavior of multilateral wells in numerical models: A hybrid analytical-numerical approach (2007) Society of Petroleum Engineers Middle East Oil and Gas Show and Conference, , Paper SPE 104581 presented at the, Kingdom of Bahrain, 11-14 MarchAkgun, F., A finite element model for analyzing horizontal well BHA behavior (2004) Journal of Petroleum Science and Engineering, 42 (1-2), pp. 121-132Al-Hussainy, R., Ramey, H.J., Jr., Crawford, P.B., The flow of real gases through porous media (1966) Journal of Petroleum Technology, 18 (5), pp. 624-636. , SPE-1243-B-PAAl-Kobaisi, M., Ozkan, E., Kazemi, H., A hybrid numerical/analytical model of a finite-conductivity vertical fracture intercepted by a horizontal well (2006) Society of Petroleum Engineers Reservoir Evaluation & Engineering, 9 (4), pp. 345-355. , SPE-92040-PAAl-Mohannadi, N., Ozkan, E., Kazemi, H., Grid-system requirements in numerical modeling of pressure-transient tests in horizontal wells (2007) Society of Petroleum Engineers Journal, 10 (2), pp. 122-131. , SPE-92041-PABeam, R.M., Warming, R.F., Yee, H.C., Stability analysis of numerical boundary conditions and implicit difference approximation for hyperbolic equations (1981) Symposium Held at NASA Ames Research Center, , Moffett Field, CA, 19-20 OctoberBerry, R.A., Notes on the PCICE method: Simplification, generalization and compressibility properties (2006) Journal of Computational Physics, 215 (1), pp. 6-11Bonner, M.P., (2007) Compressible Subsonic Flow on A Staggered Grid, , Master's Thesis, The University of British Columbia, Canada (October 2007)Bourdet, D., Whittle, T.M., Douglas, A., Pirard, Y.M., A new set of type curves simplifies well test analysis (1983) World Oil, pp. 95-106. , May 1983Chaudhry, A.U., (2003) Gas Well Testing Handbook, , Burlington, Massachusetts: ElsevierColebrook, C.F., Turbulent flow in pipes with particular reference to the transition region beteween the smooth and rough pipe laws (1939) Journal Institutational Civil Engineering, 11 (4), pp. 133-156Collins, D.A., Nghiem, L.X., Sharma, R., Agarwal, R.K., Jha, K.N., Fieldscale simulation of horizontal wells with hybrid grids (1991) Society of Petroleum Engineers Symposium on Reservoir Simulation, , Paper SPE 21218-MS presented at the, Anaheim, California, USA, 17-20 FebruaryDing, D.Y., Coupled simulation of near-wellbore and reservoir models (2011) Journal of Petroleum Science and Engineering, 76 (1-2), pp. 21-36Dranchuk, P.M., Abou-Kassem, J.H., Calculation of Z factors for natural gases using equations of state (1975) Journal of Canadian Petroleum Technology, 14 (3), pp. 34-36Dumkwu, F.A., Islam, A.W., Carlson, E.S., Review of well models and assessment of their impacts on numerical reservoir simulation perfomance (2012) Journal of Petroleum Science and Engineering, 82-83, pp. 174-186Ehrl, E., Schueler, S.K., Simulation of a tight gas reservoir with horizontal multifractured wells (2000) Society of Petroleum Engineers European Petroleum Conference, , Paper SPE 21218-MS presented at the, Paris, France, 24-25 OctoberErtekin, T., Abou-Kassem, J.H., King, G.R., (2001) SPE Textbook Series, Basic Applied Reservoir Simulation, 7. , Richardson, Texas, USA: Society of Petroleum EngineersEzulike, O., Igbokoyi, A., Horizontal well pressure transient analysis in anisotropic composite reservoirs - A three dimensional semi-analytical approach (2012) Journal of Petroleum Science and Engineering, 96-97, pp. 120-139Fuentes-Nucamendi, F.A., (1996) A Two-phase Transient Wellbore Model for Simulation of Pressure Data, , PhD's Thesis, Tulsa UniversityGalindo, J., Tiseira, A., Fajardo, P., Navarro, R., Coupling methodology of ID finite difference and 3D finite volume CFD codes based on the method of characteristics (2011) Mathematical and Computer Modelling, 54 (7-8), pp. 1738-1746Gui, P., Cunha, L.B., Cunha, J.C., A numerical two-phase flow model coupling reservoir and multisegment horizontal well (2007) Rocky Mountain Oil & Gas Technology Symposium, , Paper SPE 107989-MS presented at the, Denver, Colorado, USA, 16-18 AprilHawkins, M.F., Jr., A note on the skin effect (1956) Journal of Petroleum Technology, Trans., AIME, 8 (12), pp. 65-66. , SPE-732-GHestenes, M.R., Stiefel, E., Method of conjugate gradients for solving linear systems (1952) Journal of Research of the National Bureau Standards, 49 (6), pp. 409-436Hubbert, M.K., Darcy's law and the field equations of the flow of underground fluids (1956) Journal of Petroleum Technology, 207, pp. 222-239. , SPE-749-GKhoriakov, V., Johansen, A.C., Johansen, T.E., Transient flow modeling of advanced wells (2012) Journal of Petroleum Science and Engineering, 86-87, pp. 99-110Lee, J., Rollins, J.B., Spivey, J.P., (2003) Pressure Transient Testing, , Richardson, USA: Society of Petroleum EngineersLee, A.L., Gonzalez, M.H., Eakin, B.E., The viscosity of natural gases (1966) Journal of Petroleum Technology, 18 (8), pp. 997-1000. , SPE-1340-PAOuyang, L.B., Arbabi, S., Aziz, K., General wellbore flow model for horizontal, vertical, and slanted well completions (1998) Society of Petroleum Engineers Journal, 3 (2), pp. 124-133. , SPE-36608-PAPatankar, S.V., Spalding, D.B., A calculation procedure for heat, mass and momentum transfer in three-dimensional parabolic flows (1972) International Journal of Heat and Mass Transfer, 15 (10), pp. 1787-1806Pedrosa, O.A., Jr., Aziz, K., Use of hybrid grid in reservoir simulation (1986) Society of Petroleum Engineers Reservoir Engineering, 1 (6), pp. 611-621. , SPE-13507-PAPenmatcha, V.R., Aziz, K., Comprehensive reservoir/wellbore model for horizontal wells (1999) Society of Petroleum Engineers Journal, 4 (3), pp. 224-234. , SPE-57194-PAPenmatcha, V.R., Arbabi, S., Aziz, K., Effects of pressure drop in horizontal wells and optimum well length (1999) Society of Petroleum Engineers Journal, 4 (3), pp. 215-223. , SPE-57193-PAPourafshary, P., Varavei, A., Sepehrnoori, K., Podio, A., A compositional wellbore/reservoir simulator to model multiphase flow and temperature distribution (2009) Journal of Petroleum Science and Engineering, 69 (1-2), pp. 40-52Saad, Y., (2003) Iterative Methods for Sparse Linear Systems, 2. , ed. USA: Society of Industrial and Applied MathematicsShterev, K.S., Štefanov, S.K., Pressure based finite volume method for calculation of compressible viscous gas flow (2010) Journal of Computational Physics, 229 (2), pp. 461-480Stone, T.W., Edmunds, N.R., Kristoff, B.J., A comprehensive wellbore/reservoir simulator (1989) Society of Petroleum Engineers Symposium on Reservoir Simulation, , Paper SPE 18419-MS presented at the, Houston, Texas, USA, 06-08 FebruaryStrikwerda, J.C., (2004) Finite Difference Schemes and Partial Differential Equations, 2. , ed. Philadelphia, USA: Society of Industrial and Applied MathematicsVersteeg, H.K., Malalasekera, W., (2007) Introduction to Computational Fluid Dynamics, the Finite Volume Method, 2. , ed. Edinburgh Gate, England: Pearson Education LimitedVestergaard, H., Olsen, H., Sikandar, A.S., Noman, R., The application of unstructured-gridding techniques for full-field simulation of a giant carbonate reservoir developed with long horizontal wells (2008) Society of Petroleum Engineers Reservoir Evaluation & Engineering, 11 (6), pp. 958-967. , SPE 120887-PAVicente, R., Sarica, C., Ertekin, T., A numerical model coupling reservoir and horizontal well flow dynamics: Transient behavior of single-phase liquid and gas flow (2000) Society of Petroleum Engineers/CIM International Conference on Horizontal Well Technology, , Paper SPE 65508-MS presented at the, Calgary, Alberta, Canada, 6-8 NovemberVicente, R., Sarica, C., Ertekin, T., An investigation of horizontal well completions using a two-phase model coupling reservoir and horizontal well flow dynamics (2001) Society of Petroleum Engineers Annual Technical Conference and Exhibition, , Paper SPE 71601-MS presented at the, New Orleans, Louisiana, USA, 30 September-3 OctoberWang, X., The state-of-the-art in natural gas production (2009) Journal of Natural Gas Science and Engineering, 1 (1), pp. 14-24Wang, X., Key issues in high-perfomance natural gas wells (2012) Journal of Natural Gas Science and Engineering, 6, pp. 37-42Zhao, L., (1994) Horizontal Well Simulation with Local Grid Refinement, , PhD's Thesis, Texas A & M University, Texas, USA (August 1994)