Quantitative Formation Permeability Evaluation from Stoneley Waves

2000 ◽  
Vol 3 (02) ◽  
pp. 109-117 ◽  
Author(s):  
A. Brie ◽  
T. Endo ◽  
D.L. Johnson ◽  
F. Pampuri

Summary The borehole Stoneley wave has been shown to be sensitive to fluid mobility, the ratio of permeability, and viscosity. The phenomenon is well described using the Biot theory and the effect of the mudcake was modeled as an elastic membrane. An inversion technique, which uses both slowness and attenuation of the Stoneley wave over a range of frequency to evaluate mobility, was proposed. This paper describes the implementation of an interpretation methodology based on this technique. The error analysis shows that an accurate determination of the fluid mobility requires that some critical parameters, such as the mud slowness, mud attenuation, and the pore-fluid modulus, be precisely determined. An interpretation procedure is proposed to determine these parameters with good accuracy. The fluid mobility can then be determined without the need for external calibration with another measurement. The intrinsic permeability of the rock can then be derived, knowing the various fluid components, their relative permeabilities, and their respective viscosity. When using the proposed methodology within the applicable limits, the Stoneley wave can provide a continuous estimation of the formation permeability along the well. Core measurements are not required for calibration although they can be used for verification. This technique will find applications in reservoir engineering optimization of well production through better placement of the perforated intervals. Introduction Permeability information is essential for oil and gas production, once reserves have been identified and evaluated, to optimize well completion and field development. Permeability is needed to determine the optimal perforated interval with respect to the reservoir boundaries and the water table. More generally, permeability is needed for:completion and production optimization to maximize production while minimizing water cut,production prediction and planning to maximize hydrocarbon recovery, anddefinition of drainage pattern. Although the absolute value of permeability in the reservoir is usually considered to be the most important, the variations of permeability along the well are equally important. Permeability, however, is one of the most difficult measurements to get in an oil well. Direct measurements either provide only a few points along the well, as is the case with well testing or wireline testers, or provide measurements under different conditions in the case of core measurements. With indirect measurements, permeability is inferred from a different property (porosity, nuclear magnetic resonance, or geochemical logs) using models and assumptions. As the models are not exact, the uncertainty attached to the results is high. Another technique, the study of invasion profiles, give only qualitative information about permeability. The Stoneley wave is the only technique to provide a continuous, direct measurement of permeability along the well. However, although the principle of the measurement has been known for quite some time, obtaining a reliable and accurate measurement of permeability from Stoneley waves has proved difficult. At low frequency, the Stoneley mode becomes the tube wave and propagates as a piston-like compression of the borehole fluid in the borehole. When the borehole crosses permeable zones or permeable fractures, some fluid movement occurs between the borehole and the formation. This results in some energy loss, hence attenuation, and a slowing down of the wave, hence increased Stoneley wave slowness (Fig. 1). Fractures and permeable zones have different characteristics and affect the Stoneley wave in different ways. In particular, in the case of permeable fractures, the strong, localized impedance contrasts also cause reflections of the Stoneley wave that appear as chevron patterns on a variable density log (VDL) display. Specific techniques are used to evaluate permeable fractures with the Stoneley wave.1–3 The objective of the present work is to evaluate the permeability of nonfractured reservoirs, i.e., essentially, distributed permeability from the pore space. In effect, the formation parameter the Stoneley wave measured is not exactly the formation permeability, but rather the fluid mobility (i.e., the ratio of permeability to fluid viscosity, ko/µp). The permeability is evaluated in millidarcy units and the viscosity in centipoise, so that the mobility is usually given in md/cp. As the viscosity of water is about 1 cp, mobility and permeability take the same value in water-bearing reservoirs. The pore fluid viscosity is normally known with sufficient accuracy, hence a measurement of permeability can be obtained from the Stoneley wave. There have been many attempts at evaluating permeability from the Stoneley wave. Rosenbaum proposed measuring permeability with the Stoneley wave as early as 1974.4 In 1984, Williams et al. showed conclusive correlations between permeability and Stoneley attenuation in field logs.5 Following this idea, many log interpreters tried to develop empirical correlations between the Stoneley wave energy and permeability. But calibration with other information was needed in all cases, and the reported successes were, unfortunately, followed by disappointing results. Theoretical models based on the Biot poroelastic theory were developed by Schmitt et al.6 and Chang et al.7 Comparisons with laboratory experiments made by Winkler et al.8 provided a validation of these models. Based on these models, a simplified algorithm was developed that estimates permeability from the difference between the Stoneley slowness measured at a given frequency and the slowness calculated using a purely elastic, nonpermeable formation model. This method, called the S-Se technique, is currently proposed commercially.3 It produces useful indications in hard rocks but is often unsuccessful in soft rocks. Tang et al. in 1984, proposed a simplification of the Biot-Rosenbaum model and developed an inversion technique for the Stoneley wave amplitude.9 Later, they studied the effect of the presence of the tool in the borehole.10 Cassell et al. in 1994 proposed a simple technique to extract a permeability indicator from the variation of the Stoneley attenuation with frequency.11

Geophysics ◽  
1986 ◽  
Vol 51 (1) ◽  
pp. 137-147 ◽  
Author(s):  
Jeffry L. Stevens ◽  
Steven M. Day

We apply an iterative, linearized inversion method to Stoneley waves recorded on acoustic logs in a borehole. Our objective is to assess inversion of Stoneley wave phase and group velocity as a practical technique for shear velocity logging in slow formations. Indirect techniques for shear logging are of particular importance in this case because there is no shear head wave arrival. Acoustic logs from a long‐spaced sonic tool provided high‐quality, low‐noise data in the 1 to 10 kHz band for this experiment. A shear velocity profile estimated by inversion of a 60 ft (18 ⋅ 3 m) section of full‐wave acoustic data correlates well with the P‐wave log for the section. The inferred shear velocity ranges from 60 to 90 percent of the sound velocity of the fluid. Formal error estimates on the shear velocity are everywhere less than 5 percent. Moreover, application of the same inversion method to synthetic waveforms corroborates these error estimates. Finally, a synthetic acoustic waveform computed from inversion results is an excellent match to the observed waveform. On the basis of these results, we conclude that Stoneley‐wave inversion constitutes a practical, indirect, shear‐logging technique for slow formations. Success of the shear‐logging method depends upon availability of high‐quality, low‐noise waveform data in the 1 to 4 kHz band. Given good prior estimates of compressional velocity and density of the borehole fluid, only rough estimates of borehole radius and formation density and compressional velocity are required. The existing inversion procedure also yields estimates of formation Q inferred from spectral amplitudes of Stoneley waves. This extension of the method is promising, since amplitudes of Stoneley waves in a slow formation are highly sensitive to formation Q. Attenuation caused by formation Q dominates over attenuation caused by fluid viscosity if the viscosity is less than about [Formula: see text]. However, Stoneley‐wave amplitudes are also sensitive to gradients in shear velocity in the direction of propagation. In some cases, correction for the effects of shear‐velocity gradients is required to obtain the formation Q from Stoneley‐wave attenuation.


2019 ◽  
Vol 14 (1) ◽  
pp. 52-58 ◽  
Author(s):  
A.D. Nizamova ◽  
V.N. Kireev ◽  
S.F. Urmancheev

The flow of a viscous model fluid in a flat channel with a non-uniform temperature field is considered. The problem of the stability of a thermoviscous fluid is solved on the basis of the derived generalized Orr-Sommerfeld equation by the spectral decomposition method in Chebyshev polynomials. The effect of taking into account the linear and exponential dependences of the fluid viscosity on temperature on the spectral characteristics of the hydrodynamic stability equation for an incompressible fluid in a flat channel with given different wall temperatures is investigated. Analytically obtained profiles of the flow rate of a thermovisible fluid. The spectral pictures of the eigenvalues of the generalized Orr-Sommerfeld equation are constructed. It is shown that the structure of the spectra largely depends on the properties of the liquid, which are determined by the viscosity functional dependence index. It has been established that for small values of the thermoviscosity parameter the spectrum compares the spectrum for isothermal fluid flow, however, as it increases, the number of eigenvalues and their density increase, that is, there are more points at which the problem has a nontrivial solution. The stability of the flow of a thermoviscous fluid depends on the presence of an eigenvalue with a positive imaginary part among the entire set of eigenvalues found with fixed Reynolds number and wavenumber parameters. It is shown that with a fixed Reynolds number and a wave number with an increase in the thermoviscosity parameter, the flow becomes unstable. The spectral characteristics determine the structure of the eigenfunctions and the critical parameters of the flow of a thermally viscous fluid. The eigenfunctions constructed in the subsequent works show the behavior of transverse-velocity perturbations, their possible growth or decay over time.


2019 ◽  
Vol 11 (19) ◽  
pp. 5283 ◽  
Author(s):  
Gowida ◽  
Moussa ◽  
Elkatatny ◽  
Ali

Rock mechanical properties play a key role in the optimization process of engineering practices in the oil and gas industry so that better field development decisions can be made. Estimation of these properties is central in well placement, drilling programs, and well completion design. The elastic behavior of rocks can be studied by determining two main parameters: Young’s modulus and Poisson’s ratio. Accurate determination of the Poisson’s ratio helps to estimate the in-situ horizontal stresses and in turn, avoid many critical problems which interrupt drilling operations, such as pipe sticking and wellbore instability issues. Accurate Poisson’s ratio values can be experimentally determined using retrieved core samples under simulated in-situ downhole conditions. However, this technique is time-consuming and economically ineffective, requiring the development of a more effective technique. This study has developed a new generalized model to estimate static Poisson’s ratio values of sandstone rocks using a supervised artificial neural network (ANN). The developed ANN model uses well log data such as bulk density and sonic log as the input parameters to target static Poisson’s ratio values as outputs. Subsequently, the developed ANN model was transformed into a more practical and easier to use white-box mode using an ANN-based empirical equation. Core data (692 data points) and their corresponding petrophysical data were used to train and test the ANN model. The self-adaptive differential evolution (SADE) algorithm was used to fine-tune the parameters of the ANN model to obtain the most accurate results in terms of the highest correlation coefficient (R) and the lowest mean absolute percentage error (MAPE). The results obtained from the optimized ANN model show an excellent agreement with the laboratory measured static Poisson’s ratio, confirming the high accuracy of the developed model. A comparison of the developed ANN-based empirical correlation with the previously developed approaches demonstrates the superiority of the developed correlation in predicting static Poisson’s ratio values with the highest R and the lowest MAPE. The developed correlation performs in a manner far superior to other approaches when validated against unseen field data. The developed ANN-based mathematical model can be used as a robust tool to estimate static Poisson’s ratio without the need to run the ANN model.


Geophysics ◽  
2020 ◽  
Vol 85 (4) ◽  
pp. MR201-MR212
Author(s):  
Zhi-Qiang Yang ◽  
Tao He ◽  
Chang-Chun Zou

Velocity dispersion is a common phenomenon for fluid-charged porous rocks and carries important information on the pore structure and fluid in reservoir rocks. Previous ultrasonic experiments had measured more significant non-Biot velocity dispersion on saturated reservoir sandstones with increasing pore-fluid viscosity. Although wave-induced local squirt-flow effect could in theory cause most of the non-Biot velocity dispersion, its quantitative prediction remains a challenge. Several popular models were tested to predict the measured velocities under undrained conditions, but they either underestimated the squirt-flow effect or failed to simultaneously satisfy P- and S-wave velocity dispersions (especially for higher viscosity fluids). Based on the classic double-porosity theory that pore space is comprised of mainly stiff/Biot’s porosity and minor compliant porosity, an effective “wet frame” was hypothesized to account for the squirt-flow effect, whose compliant pores are filled with a hypothesized fluid with dynamic modulus. A new dynamic elastic model was then introduced by extending Biot theory to include the squirt-flow effect, after replacing the dry-frame bulk/shear moduli with their wet-frame counterparts. In addition to yielding better velocity predictions for P- and S-wave measurements of different fluid viscosities, the new model is also more applicable because its two key tuning parameters (i.e., the effective aspect ratio and porosity of compliant pores) at in situ reservoir pressure could be constrained with laboratory velocity measurements associated with pore-fluid viscosities.


2019 ◽  
Vol 12 (10) ◽  
pp. 5503-5517 ◽  
Author(s):  
Pascal Hedelt ◽  
Dmitry S. Efremenko ◽  
Diego G. Loyola ◽  
Robert Spurr ◽  
Lieven Clarisse

Abstract. The accurate determination of the location, height, and loading of sulfur dioxide (SO2) plumes emitted by volcanic eruptions is essential for aviation safety. The SO2 layer height is also one of the most critical parameters with respect to determining the impact on the climate. Retrievals of SO2 plume height have been carried out using satellite UV backscatter measurements, but, until now, such algorithms are very time-consuming. We have developed an extremely fast yet accurate SO2 layer height retrieval using the Full-Physics Inverse Learning Machine (FP_ILM) algorithm. This is the first time the algorithm has been applied to measurements from the TROPOMI instrument onboard the Sentinel-5 Precursor platform. In this paper, we demonstrate the ability of the FP_ILM algorithm to retrieve SO2 plume layer heights in near-real-time applications with an accuracy of better than 2 km for SO2 total columns larger than 20 DU. We present SO2 layer height results for the volcanic eruptions of Sinabung in February 2018, Sierra Negra in June 2018, and Raikoke in June 2019, observed by TROPOMI.


2010 ◽  
Vol 02 (03) ◽  
pp. 585-615 ◽  
Author(s):  
BORIS D. PLYSHCHENKOV ◽  
ANATOLY A. NIKITIN

Numerical experiments based on Pride's model of electrokinetic phenomena have shown that electromagnetic Stoneley waves as well as pressure Stoneley waves are most sensitive to permeability variations. A new way for quantitative evaluation of any value of formation permeability is presented. It is based on simultaneous measurement of pressure field and axial component of electric field excited by an acoustic source in fluid-filled borehole with help from a set of receivers in borehole. Frequency dependence of ratio of the complex-valued amplitudes of the electric Stoneley wave to the pressure Stoneley wave obtained as a result of plane waves decomposition of pressure field and mentioned component of electric field carries important information about permeability. The ratio of the real part of this ratio to its imaginary part is very sensitive to permeability variations. The approximate analytical expressions for this ratio derived for open and sealed pores on borehole wall are base for construction of a new way of quantitative evaluation of formation permeability.


Author(s):  
Sergey V. Kuznetsov

Stoneley wave velocity variation is analyzed by solving the modified Scholte secular equation for velocity of Stoneley waves, allowing to find dependency of the Stoneley wave velocity on the Wiechert parameter and construct a set of inequalities that confines region of existence for the appropriate root of the secular equation. Numerical analysis for Stoneley wave velocity dependence on the Wiechert parameter for both auxetics (materials with negative Poisson’s ratio) and nonauxetics revealed the presence of (i) asymptotes indicating degeneracy of Stoneley waves into the corresponding Rayleigh waves; and (ii) common extremums relating to degeneracy of Stoneley waves into the corresponding bulk shear waves.


1988 ◽  
Vol 6 (2) ◽  
pp. 118-125
Author(s):  
D.J. Norris

TCPL has recently carried out ‘reconnaissance or exploration’ 3-D surveys, in three different blocks, each designed to solve a different type of problem. In each case a considerable improvement in data quality and the resulting structural/stratigraphic interpretation was achieved. The Kupe South structure is a wrench-induced feature cross-cut by numerous small-medium faults. Stratigraphic changes across the prospect produce a variable quality seismic event at the top reservoir level. The Pataka Prospect comprises a narrow horst block trend within the Oakura fault zone, offshore New Plymouth. Accurate determination of potential reserves required a reliable interpretation of the fault configuration, and the amount of displacement of the reservoir horizon by the critical faults. The Waitara Prospect is affected by a ‘no-data’ zone possibly associated with volcanics in the near surface. It was necessary to define the extent and nature of the no-data zone and to calculate the effect of the interpreted volcanics on the time structure map. Prior to the surveys we modelled the effects of such critical parameters as sail-line separation, final interpolation spacing and the dimensions of the 3-D grids using existing 2-D data. Good results were obtained with a wider line spacing than is strictly required for true 3-D. Whilst the Reconnaissance 3-D method has not removed all of the difficulties with interpretation, a considerable improvement was obtained in data quality and ease of interpretation.


Sign in / Sign up

Export Citation Format

Share Document