Measurement and Modelling of Gas Condensate Flow Around Rock Perforation

In the petroleum production industry wells are mostly cased and perforated in the producing formation. Perforation characteristics such as size, length, number of perforation tunnels and their arrangements as well as fluid and rock properties determine fluid flow behaviour in the wellbore region, hence, well productivity. Flow of gas and condensate around a perforation tunnel (including the damaged zone) has been studied by performing steady-state core experiments and simulating the results numerically, using a finite element modelling approach. The model allows for the changes in fluid properties and accounts for the coupling of the two phases and the inertial effect using a fractional flow based correlation. The results indicated that different sets of thickness-permeability (h−k) values obtained from matching single-phase flow performance could be assigned to the damaged zone around perforation to represent the two-phase flow performance. The status of the tip of the perforation for two extreme cases of totally closed and fully open was investigated and found to have a minimal effect on the performance of the system.     

Modeling Water Leak-off Behavior in Hydraulically Fractured Gas Shale under Multi-mechanism Dominated Conditions

Fracturing-fluid leak-off in fractured gas shale is a complex process involving multiple pore/fluid transports and interactions. However, water leak-off behavior has not been modeled comprehensively by considering the multi-pores and multi-mechanisms in shale with existing simulators. In this paper, we present the development of a comprehensive multi-mechanistic, multi-porosity, and multi-permeability water/gas flow model that uses experimentally determined formation properties to simulate the fracturing-fluid leak-off of hydraulically fractured shale gas wells. The multi-mechanistic model takes into account water transport driven by hydraulic convection, capillary and osmosis, gas transport caused by hydraulic convection, and salt ion transport caused by advection and diffusion. The multi-porosity includes hydraulic fracture millipores, organic nanopores, clay nanopores, and other inorganic micropores. The multi-permeability model accounts for all the important processes in shale system, including gas adsorption on the organics’ surface, multi-mechanistic clay/other inorganic mineral mass transfer, inorganic mineral/hydraulic fracture mass transfer, and injection from a hydraulically fractured wellbore. The dynamic water saturation and pressure profiles within clay and other inorganic matrices are compared, revealing the leak-off behavior of water in rock media with different physicochemical properties. In sensitivity analyses, cases with different clay membrane efficiency, volume proportion of source rock, connate water salinity, and saturation are considered. The impacts of shale properties on water fluxes through wellbore, hydraulic fracture and matrix, and the total injection and leak-off volumes of the well during the treatment of hydraulic fracturing are investigated. Results show that physicochemical properties in both organic and inorganic matrices affect the water leak-off behavior.     

The Dual-Reciprocity Boundary Element Analysis for Hydraulically Fractured Shale Gas Reservoirs Considering Diffusion and Sorption Kinetics

This work presents a new application of boundary element method (BEM) to model fluid transport in unconventional shale gas reservoirs with discrete hydraulic fractures considering diffusion, sorption kinetics and sorbed-phase surface diffusion. The fluid transport model consists of two governing partial differential equations (PDEs) written in terms of effective diffusivities for free and sorbed gases, respectively. Boundary integral formulations are analytically derived using the fundamental solution of the Laplace equation for the governing PDEs and Green’s second identity. The domain integrals arising due to the time-dependent function and nonlinear terms are transformed into boundary integrals employing the dual-reciprocity method. This transformation retains the domain-integral-free, boundary-integral-only character of standard BEM approaches. In the proposed solution, the free- and sorbed-gas flow in the shale matrix is solved simultaneously after coupling the fracture flow equation of free gas. Well production performance under the effect of relaxation phenomenon due to delayed responses of sorbed gas under nonequilibrium sorption condition is rigorously captured by imposing the zero-flux condition at fracture–matrix interface for the sorbed-gas transport equation. The validity of proposed solution is verified using several case studies through comparison against a commercial finite-element numerical simulator.

     

Pressure Transient Analysis for Multi-stage Fractured Horizontal Wells in Shale Gas Reservoirs

This article presents the PTA on the multi-stage fractured horizontal well in shale gas reservoirs incorporating desorption and diffusive flow in the matrix. Currently, most PTA models are simply based on Darcy flow both in natural fractures and matrix without considering the mechanisms of desorption and diffusion in shale matrix. Source function and Laplace transform with the numerical discrete method are employed to solve the mathematical model. The solution is presented in the Laplace domain so that the wellbore storage effect and skin factor can be easily incorporated by convolution. Type curves are plotted with Stehfest algorithm and different flow regimes are identified. The presented model could be used to interpret pressure signals more accurately for shale gas reservoirs.     

Landfill Gas Flow: Collection by Horizontal Wells

Transport in Porous Media - Collection of landfill gas by horizontal perforated wells is studied. The problem combines flow through porous media in the landfill and unobstructed pipe flow in the...     

Microvisual Study of Multiphase Gas Condensate Flow in Porous Media

Gas condensate reservoirs constitute a significant portion of hydrocarbon reserves worldwide. The liquid drop-out in these reservoirs may lead to recovery problems such as near wellbore permeability impairment and uncertainty in the actual location of the target condensate. Such technical issues can be addressed through improved understanding of the formation of condensate and the multiphase flow of gas/condensate/water in the reservoir as characterized by relative permeability curves. The appropriate relative permeability curves in turn can be used in reservoir simulators to assist in optimization of field development. This paper reports results of experiments conducted in micromodels, in support of possible core flow tests, using reservoir fluids under reservoir conditions. In particular, visualizations of condensate formation with and without connate water are presented and the differences between the two cases as well as the possible implications for the relative permeability measurements are discussed. Furthermore, the flow of gas and condensate at different force ratios (capillary and Bond numbers) are presented. It is postulated that a single dimensionless number may not be sufficient to characterize the multiphase flow in gas condensate reservoirs. The physical mechanisms occurring under various field conditions are examined in the light of these observations.     

Well Test Solutions for Vertically Fractured Injection Wells

The transient behavior of a vertically fractured pressure response due to the presence of an infinite-conductivity vertical fracture is determined by solving the diffusivity equation in elliptical coordinates. The solution is then extended to a composite elliptical system to provide for the different fluid banks present during water injection. The validity of the analytical solutions presented is demonstrated by comparing limiting forms with those available elsewhere in the literature. Computational issues which became evident during the verification stage of our work are also discussed. The solutions have been developed in the Laplace domain to facilitate the addition of fissures and variable rate production (i.e. wellbore storage).A pressure transient test for tracking the advancement of a water front during the early stages of waterflooding is described. We utilize the composite elliptical model developed herein to provide for two distinct regions in which the flow behavior resulting from an induced fracture is elliptical rather than radial. A relationship between the increasing elliptical distance to the waterfront and the resulting change in the apparent (total) skin factor is obtained. Through the analysis of successive falloff tests, this relationship may be used to monitor the advancement of the front provided the cumulative volume of injected water is known. The fluid saturations and the mobilities of the swept and unswept regions are assumed unknown and are obtained from the test analysis.Finally, we present methods for computing the Mathieu functions necessary in solving the diffusivity equation in elliptical coordinates. Mathieu functions are utilized in many applications involving elliptical geometry and we feel the efficient evaluation of these functions is an important contribution of this work.     

Pore-Level Modeling of Gas and Condensate Flow in Two- and Three-Dimensional Pore Networks: Pore Size Distribution Effects on the Relative Permeability of Gas and Condensate

We present a mechanistic model of retrograde condensation processes in two- and three-dimensional capillary tube networks under gravitational forces. Condensate filling-emptying cycles in pore segments and gas connection–isolation cycles are included. With the pore-level distribution of gas and condensate in hand, we determine their corresponding relative permeabilities. Details of pore space and displacement are subsumed in pore conductances. Solving for the pressure field in each phase, we find a single effective conductance for each phase as a function of condensate saturation. Along with the effective conductance for the saturated network, the relative permeability for each phase is calculated. Our model porous media are two- and three-dimensional regular networks of pore segments with distributed size and square cross-section. With a Monte Carlo sampling we find the optimum network size to avoid size effects and then we investigate the effect of network dimensionality and pore size distribution on the relative permeabilities of gas and condensate.     

Effect of Hydrogen Sulfide on Condensation and Flow of Gas Condensate Fluids in Porous Media

Condensation and flow experiments were conducted at subsurface conditions in a glass micromodel using reservoir fluids with and without the hydrogen sulfide component. It has been noted that the formation of the condensing phase as well as modes of condensate flow are similar for both fluids. Furthermore, an additional condensate transport mechanism, termed lamella flow, was observed with the sour fluid. It has been concluded that core flow experiments conducted with sweet reservoir fluid should reproduce the flow of sour fluid to a large extent.     

Gas Condensate Relative Permeability for Well Calculations

This paper addresses several issues related to the modeling and experimental design of relative permeabilities used for simulating gas condensate well deliverability. Based on the properties of compositional flow equations, we make use of the fact that relative permeability ratio k _rg/k _ro is a purely thermodynamic variable, replacing saturation, when flow is steady-state. The key relation defining steady-state flow in gas condensate wells is relative permeability k _rg as a function of k _rg/k _ro. Consequently, determination of saturation and k _r as a function of saturation is not important for this specific calculation. Once the k _rg=f(k _rg/k _ro) relationship is experimentally established and correlated with capillary number (N _c), accurate modeling of condensate blockage is possible. A generalized model is developed for relative permeability as the function of k _rg/k _ro and capillary number. This model enables us to link the immiscible or rock curves with miscible or 'straight-line curves by a transition function dependent on the capillary number. This model is also extended to the case of high-rate, inertial gas flow within the steady-state condensate blockage regionand locally at the wellbore. We have paid particular attention to the effect of hysteresis on the relation k _rg=f(k _rg/k _ro), based on our observation that many repeated cycles of partial/complete imbibition and drainage occur in the near-well region during the life of a gas condensate well. Finally, the composite effect of condensate blockage is handled using a Muskat pseudopressure model, where relative permeabilities are corrected for the positive effect of capillary number dependence and the negative effect of inertial high velocity flow. Special steady-state experimental procedures have been developed to measure k _rg as a function of k _rg/k _ro and N _c. Saturations, though they can be measured, are not necessary. An approach for fitting steady-state gas condensate relative permeability data has been developed and used for modeling relative permeability curves.     

Finite-Volume Discretisations for Flow in Fractured Porous Media

Over the last decades, finite-volume discretisations for flow in porous media have been extended to handle situations where fractures dominate the flow. Successful discretisations have been based on the discrete fracture-matrix models to yield mass conservative methods capable of explicitly incorporating the impact of fractures and their geometry. When combined with a hybrid-dimensional formulation, two central concerns are the restrictions arising from small cell sizes at fracture intersections and the coupling between fractures and matrix. Focusing on these aspects, we demonstrate how finite-volume methods can be efficiently extended to handle fractures, providing generalisations of previous work. We address the finite-volume methods applying a general hierarchical formulation, facilitating implementation with extensive code reuse and providing a natural framework for coupling of different subdomains. Furthermore, we demonstrate how a Schur complement technique may be used to obtain a robust and versatile method for fracture intersection cell elimination. We investigate the accuracy of the proposed elimination method through a series of numerical simulations in 3D and 2D. The simulations, performed on fractured domains containing permeability heterogeneity and anisotropy, also demonstrate the flexibility of the hierarchical framework.     

Flow Field Around Single and Tandem Piers

The present study provides a comparison between the flow pattern around two circular piers in tandem and a single pier set up on a moderately rough flat bed in a laboratory flume. Velocities are measured by an Acoustic Doppler Velocimeter (ADV). The contours of the time-averaged velocity components, Reynolds shear stresses, turbulence intensities and turbulence kinetic energy at different planes are presented. Streamlines and vectors are used to study the flow features. The analysis of power spectra around the piers is also presented. The results show that the presence of downstream pier changes the flow structure to a great extent, particularly in the near-wake region. Within the gap between the two piers, a stronger and substantial upflow is shaped. However, a weaker transverse-deflection is formed in comparison with that in the single pier. Near the bed, the velocity of flow approaching the downstream pier decreases to 0.2–0.3 times of the approach mean velocity (U ₀) due to the sheltering effect of the upstream pier. In the wake of downstream pier, the flow structure is completely different from the one in the wake of single pier. In comparison with the single pier, the values of turbulence kinetic energy and turbulence intensities show a considerable decrease around the tandem piers. In tandem piers, the high values of turbulence characteristics are found near the downstream pier. There is a recirculation zones just upstream of the sheltered pier close to the bed and another behind that pier near the free surface. The results show a decrease in the strength of vortical structure in the wake of tandem piers in comparison with single pier. It is shown that the formation of flow with different Reynolds number along the flow depth due to the effect of bed roughness, as well as pier spacing, can influence the type of flow regime of tandem case. In addition to enhancing the flow structure indulgence, the present detailed measurements can also be used for verification of numerical models.     

Simulation of Pressure Transient Behavior for Asymmetrically Finite-Conductivity Fractured Wells in Coal Reservoirs

Based on Fick’s law in matrix and Darcy flow in cleats and hydraulic fractures, a new semi-analytical model considering the effects of boundary conditions was presented to investigate pressure transient behavior for asymmetrically fractured wells in coal reservoirs. The new model is more accurate than previous model proposed by Anbarci and Ertekin, SPE annual technical conference and exhibition, New Orleans, 27–30 Sept 1998 because new model is expressed in the form of integral expressions and is validated well through numerical simulation. (1) In this paper, the effects of parameters including fracture conductivity, coal reservoir porosity and permeability, fracture asymmetry factor, sorption time constant, fracture half-length, and coalbed methane (CBM) viscosity on bottomhole pressure behavior were discussed in detail. (2) Type curves were established to analyze both transient pressure behavior and flow characteristics in CBM reservoir. According to the characteristics of dimensionless pseudo pressure derivative curves, the process of the flow for fractured CBM wells was divided into six sub-stages. (3) This paper showed the comparison of transient steady state and pseudo steady state models. (4) The effects of parameters including transfer coefficient, wellbore storage coefficient, storage coefficient of cleat, fracture conductivity, fracture asymmetry factor, and rate coefficient on the shape of type curves were also discussed in detail, indicating that it is necessary to keep a bigger fracture conductivity and fracture symmetry for enhancing well production and reducing pressure depletion during the hydraulic fracturing design.     

Positive Effect of Flow Velocity on Gas–Condensate Relative Permeability: Network Modelling and Comparison with Experimental Results

Positive velocity dependency of relative permeability of gas–condensate systems, which has been observed in many different core experiments, is now well acknowledged. The above behaviour, which is due to two-phase flow coupling in condensing systems at low interfacial tension (IFT) conditions, was simulated using a 3D pore network model. The steady-dynamic bond network model developed for this purpose was also equipped with a novel anchoring technique, which was based on the equivalent hydraulic length concept adopted from fluid flow through pipes. The available rock data on the co-ordination number, capillary pressure, absolute permeability, porosity and one set of measured relative permeability curves were utilised to anchor the capillary, volumetric and flow characteristics of the constructed network model to those properties of the real core sample. Then the model was used to predict the effective permeability values at other IFT and velocity levels. There is a reasonable quantitative agreement between the predicted and measured relative permeability values affected by the coupling rate effect.     

A Triple Pore Network Model (T-PNM) for Gas Flow Simulation in Fractured,Micro-porous and Meso-porous Media

Transport in Porous Media - In this study, a novel triple pore network model (T-PNM) is introduced which is composed of a single pore network model (PNM) coupled to fractures and micro-porosities....     

The Effect of Asphalt Precipitation on Flow Behavior and Production of a Fractured Carbonate Oil Reservoir During Gas Injection

The first field data, collected over an 11 year period, are presented which indicate the possible effect of asphalt precipitation on the permeability and injectivity index of a fractured carbonate oil reservoir. The asphalt aggregates were formed during enhanced oil recovery by injection of a rich gas into the reservoir. The data indicate that, while at the initial stages of the operations the permeability and injectivity index decrease, at later times they appear to oscillate with the process time, with apparent oscillations' periods that depend on the heterogeneity of the reservoir. Two classes of plausible mechanisms that give rise to such oscillatory behavior are discussed. One relies on the changes in the structure of the reservoir's fractures, while the other one is based on asphalt precipitation in the reservoir. Computer simulations of flow and precipitation of asphalt aggregates in a pore network model of the reservoir are carried out. The results appear to support our proposition that asphalt formation and precipitation in the reservoir are the main mechanism for the observed behavior of the injectivity index. We also develop a stochastic continuum model that accurately predicts the time-dependence of the reservoir's permeability and injectivity index during the gas injection process.     

Equivalent Permeability of Fractured Media Incorporating Tortuosity and Nonlinear Flow

Transport in Porous Media - Permeability of fractured media plays an important role in many hydrogeologic applications. In this study, we propose an approach to quantify the equivalent permeability...     

管道内差压驱动机器人相关流场数值模拟研究

管道内流场对机器人的驱动力是设计管内机器人外形尺寸的基本依据,本文用数值方法计算了管道内检测机器人所受的差压驱动力。在合理提出一些基本假设后,用一阶迎风和中心差分格式离散管道内检测机器人附近流场的控制方程,用SIMPLE算法求得了不同入口流速下机器人附近的流场分布,以及流场对机器人的驱动力。结果表明雷诺数为1875时,机器人下游流场变为湍流;当雷诺数为60000时,机器人下游流场变为非定常流,出现周期性流动;计算数据还表示驱动力只与管内平均流速有关而与运行压力无关。     

高速列车紊态外流场的数值模拟研究

高速列车是近地运行的细长、庞大物体,它的空气绕流问题有其特殊性,本文以不可压缩粘性流体的Navier-Stokes方程和k-ε两方程紊流模型为基础,采用有限元方法求解了高速列车三维紊态外流场,针对有限元法应用于流场计算时常出现的问题,采用分离式解法,非对称矩阵一维变带宽压缩存储及带宽极小化等方法,最大限度地降低计算存储量;并采用罚函数法,集中质量矩阵,缩减积分法,带参数迭代法以及 引入松弛因子等技术,提出了一套用有限元法计算非线性问题的求解方法,提高了收敛速度的计算严谨,计算方法和计算结果对列车空气动力学的深入研究有一定的帮助。