Estimating petrophysical parameters and average mud-filtrate invasion rates using joint inversion of induction logging and pressure transient data

Geophysics ◽  
2011 ◽  
Vol 76 (2) ◽  
pp. E21-E34 ◽  
Author(s):  
Lin Liang ◽  
Aria Abubakar ◽  
Tarek M. Habashy

We introduce an inversion approach for determining the water-based mud-filtrate invasion profile, as well as the formation porosity and horizontal permeability, from the induction logging data. The inversion is constrained by a multiphase fluid flow simulator that simulates the mud-filtrate invasion process to obtain the spatial distributions of the water saturation and the salt concentration, which are in turn transformed into the formation resistivity using a resistivity-saturation formula. By ignoring the diffusion effect, we assume that the mud-filtrate invasion process is mainly convective so that it can be equivalently simulated by providing an average invasion rate and the duration of invasion. The average invasion rate can be directly inverted from the fluid-flow-constrained inversion of induction logging data. We also obtain the mud-filtrate invasion profile, which is consistent with the fluid flow physics. The reconstructed mud-filtrate invasion profile benefits the interpretation of the formation test. When the pressure transient data are available, this approach can be also used to jointly invert both induction logging data and pressure transient data to obtain the mud-filtrate invasion profile, as well as a parametric distribution of the TI-anisotropic formation permeability and porosity. Assuming a vertical well penetrating horizontal formations, the fluid flow problem is solved using an implicit black oil finite-difference simulator with brine tracking option based on a cylindrical, axially symmetric grid, whereas the response of the induction logging tool is simulated using a frequency-domain finite-difference solver based on a Cartesian grid. A Gauss-Newton inversion scheme using the multiplicative regularization technique is used for either the fluid-flow-constrained inversion or the joint inversion. The reliability of the inversion results depends on the accuracy of a priori knowledge of the reservoir, which needs to be confirmed via sensitivity analysis.

Geophysics ◽  
2020 ◽  
Vol 85 (4) ◽  
pp. D105-D119
Author(s):  
Joshua Bautista-Anguiano ◽  
Carlos Torres-Verdín ◽  
Joachim Strobel

The quantitative interpretation of borehole spontaneous potential (SP) measurements via Nernst’s equation often relies on limiting assumptions such as shallow mud-filtrate invasion, negligible streaming potentials, and uniform borehole symmetry. To overcome these limitations while honoring the governing physics of coupled mass transport associated with SP phenomena, we have developed a 3D finite-difference algorithm to simulate borehole SP measurements acquired across water-bearing rocks that incorporates electrochemical, membrane, and electrokinetic SP. The algorithm is based on a mechanistic description of nonequilibrium thermodynamics that enables its coupling with a fluid flow simulator to quantify the effects of continuously varying properties within permeable formations due to mud-filtrate invasion. Numerical modeling of SP measurements acquired under complex petrophysical and geometric conditions enables uncertainty quantification in the estimation of formation-water resitivity, location of bed boundaries, or detection of permeable beds while accounting for shoulder-bed effects, borehole deviation, and borehole eccentricity. Our results indicate that for well trajectories with a relative dip of less than 30°, the assumption of perpendicular beds does not entail significant errors in SP-related calculations, thereby reducing CPU time by a factor of at least 1.76. In vertical wells, SP provides the best resolution possible because deviated wells or dipping beds result in more extended and pronounced shoulder-bed effects. Furthermore, electrokinetic effects can be neglected for commonly used pressure overbalance ranges. In cases in which electrokinetic contributions are not negligible, we conclude that they are more significant when the rock permeability is in the two-figure millidarcy range. Finally, the simulation algorithm enables hypothesis testing to determine the origin and conditions under which SP shale-baseline shifts may occur. The latter shifts can signal vertical variations in salt concentration, which are crucial in the estimation of water saturation and detection of aquifer compartments.


2014 ◽  
Vol 2014 ◽  
pp. 1-5
Author(s):  
Jianhua Zhang ◽  
Zhenhua Liu

The process of drilling mud filtrate invading into a reservoir is time dependant. It causes dynamic invasion profiles of formation parameters such as water saturation, salinity, and formation resistivity. Thus, the responses of a high-definition induction log (HDIL) tool are time dependent. The logging time should be considered as an important parameter during logging interpretation for the purposes of determining true formation resistivity, estimating initial water saturation, and evaluating a reservoir. The time-dependent HDIL responses are helpful for log analysts to understand the invasion process physically. Field examples were illustrated for the application of present method.


2006 ◽  
Vol 9 (01) ◽  
pp. 39-49 ◽  
Author(s):  
Ihsan M. Gok ◽  
M. Onur ◽  
Peter S. Hegeman ◽  
Fikri J. Kuchuk

Summary This paper examines the effect of a mud-filtrate-invaded zone on pressure transients from multiprobe/packer-probe wireline formation testers (WFTs). Invasion zones are modeled as composite zones concentric with the wellbore that have different rock and fluid properties (permeability, porosity, viscosity, and compressibility) from those of the native uninvaded formation. The results show that for multiprobe wireline testers, the sink- (production) and horizontal-probe pressure responses are highly affected by the invaded-region properties, while the vertical-probe pressures are influenced mainly by the properties of the uninvaded zones. For the packer-probe configuration, similar results are obtained (i.e., the vertical-probe pressures are influenced mainly by the properties of the uninvaded zones, while the packer-interval pressures at early times are influenced by the invaded-zone properties). It is shown that if the invaded zones are incorporated into the interpretation process with a 3D r-???-z single-phase, finite-difference model like the one developed in this work, simultaneous matching of spatially available WFT pressure-data sets using nonlinear regression can provide estimates of both invaded- and uninvaded-zone parameters. A synthetic example of a multiprobe test is presented to confirm the theory and procedures developed in this work. Introduction The multiprobe and packer-probe WFTs are used to conduct controlled local production and buildup tests as well as horizontal- and vertical-interference tests (interval pressure transient tests, or IPTT). These tools provide formation-fluid samples and estimates of horizontal and vertical permeabilities and wellbore damage. Further details about WFTs, equipped with packers and multiple probes, can be found in Zimmerman et al. (1990), Goode et al. (1991), Goode and Thambynayagam (1992), Head and Bettis (1993), Pop et al. (1993), Kuchuk et al. (1994), and Onur et al. (2004). It has always been a concern how the pressure transients from these formation testers are affected by the presence of invaded regions around the wellbore. When an oil/gas well is drilled, some of the borehole fluid (mud filtrate) can leak into the formation, displacing the native formation fluid and creating an invaded zone around the wellbore. The invading fluid usually has a viscosity and compressibility that differ from those of the formation fluid. All WFT drawdown and buildup tests and IPTTs are conducted in an openhole environment, in which formation invasion takes place until the mudcake buildups, after which invasion becomes negligible. The invasion depth may vary from a few inches to a few feet, depending on formation and mud properties, as well as drilling parameters. Therefore, one must know how the invaded zone affects parameter estimates from WFT drawdown and buildup tests and IPTTs. Second, can we estimate some of the parameters of invaded zones? For instance, if we can obtain the permeabilities of both invaded and uninvaded zones, this can give us the endpoint effective permeabilities for both water and oil on a much larger scale than that from cores. In general, the process of drilling-fluid invasion is quite complicated because it involves solid and solute (and solvent, in the case of oil-based mud) transport and precipitation as well as multiphase flow, capillary pressure hysteresis, wettability alterations, chemical adsorption, and gravity effects (Ferguson and Klotz 1954; Bailey et al. 2000; Civan 2002). One of the main difficulties is how to model the initial condition set by the invasion process before sampling and pressure transients from WFTs. In addition, modeling the effects of invasion zone(s) on fluid sampling during production and on pressure-transient data from WFTs during pumpout and buildup periods is challenging, where we are faced by two distinct problems. However, if one assumes that capillary and gravity effects are negligible and that permeability and porosity impairments caused by solid and solute invasion and precipitation are minimal, then one may tackle the invasion problem by modeling it as miscible or two-phase immiscible flow. Several researchers (Phelps et al. 1984; Hammond 1991; Akram et al. 1999; Zeybek et al. 2001; Wu et al. 2002) have taken this approach to include the effect of invasion on fluid sampling and pressure-transient data from WFTs. However, all work given in Phelps et al. (1984), Hammond (1991), Akram et al. (1999), Zeybek et al. (2001), and Wu et al. (2002) is limited to a single-layer system. In addition, these references do not address the parameter-estimation problem (mainly permeability) in the presence of invasion. For instance, Akram et al. (1999) numerically simulated the invaded zone by injecting water to create a sharp interface (or shock front) between filtrate and oil in the simulator and then studied the effect of invasion on sampling quality (water cut) vs. time during the pumpout (production) period, which is production from the formation into the wellbore by means of a pump. In the presence of water-based-mud invasion, Zeybek et al. (2001) presented an inversion methodology using a numerical simulator to refine oil/water relative permeabilities by integrating dual-packer pressure and water-cut measurements with openhole array resistivity measurements. During their inversion process, they also updated the absolute horizontal and vertical permeability in the simulator. However, they did not consider vertical-observation probe pressure data in their inversion methodology. Wu et al. (2002) presented miscible and immiscible multiphase sampling modeling for oil- and water-based muds for a dual-probe configuration. Assuming that the capillary pressure and relative permeability curves (and the mud-filtrate characteristics) are known, they presented an inversion methodology of horizontal and vertical permeabilities by matching sampling quality and pressure measurements acquired in the presence of oil- and water-based invaded zones. Their inversion method is based on a neural-network approach coupled with a numerical simulator.


Geophysics ◽  
2012 ◽  
Vol 77 (5) ◽  
pp. ID9-ID22 ◽  
Author(s):  
Lin Liang ◽  
Aria Abubakar ◽  
Tarek M. Habashy

We present a fluid-flow constrained inversion approach for integrating controlled-source electromagnetic data and production data. In this approach, we assumed that the reservoir model has been well defined from a priori knowledge obtained from other independent measurements such as seismic and/or well-logs. Our objective was to reconstruct the permeability distribution and the shape and location of the flooding front. A finite-difference reservoir simulator was used to model the water-flooding process to simulate the time-dependent production data as well as the temporal and spatial distributions of water saturation and salt concentration, which were then transformed into the reservoir resistivity distribution using a petrophysical relationship. A finite-difference frequency-domain electromagnetic forward-modeling code was then employed to simulate the controlled-source electromagnetic response. The permeability distribution was reconstructed using a multiplicative-regularized Gauss-Newton algorithm for jointly inverting controlled-source electromagnetic and production data. From the water saturation distribution, we can identify the shape and location of the fluid-front resulting from the recovery process.


2010 ◽  
Vol 7 ◽  
pp. 182-190
Author(s):  
I.Sh. Nasibullayev ◽  
E.Sh. Nasibullaeva

In this paper the investigation of the axisymmetric flow of a liquid with a boundary perpendicular to the flow is considered. Analytical equations are derived for the radial and axial velocity and pressure components of fluid flow in a pipe of finite length with a movable right boundary, and boundary conditions on the moving boundary are also defined. A numerical solution of the problem on a finite-difference grid by the iterative Newton-Raphson method for various velocities of the boundary motion is obtained.


2014 ◽  
Author(s):  
Shafaruniza Mahadi ◽  
Farah Suraya Md Nasrudin ◽  
Faisal Salah ◽  
Zainal Abdul Aziz

Sign in / Sign up

Export Citation Format

Share Document