Petrophysical inversion of borehole array-induction logs: Part II — Field data examples

Geophysics ◽  
2006 ◽  
Vol 71 (5) ◽  
pp. G261-G268 ◽  
Author(s):  
Carlos Torres-Verdín ◽  
Faruk O. Alpak ◽  
Tarek M. Habashy

We describe the application of Alpak et al.’s (2006) petrophysical inversion algorithm to the interpretation of borehole array induction logs acquired in an active North American gas field. Layer-by-layer values of porosity and permeability were estimated in two closely spaced vertical wells that penetrated the same horizontal rock formation. The wells were drilled with different muds and overbalance pressures, and the corresponding electromagnetic induction logs were acquired with different tools. Rock-core laboratory measurements available in one of the two wells were used to constrain the efficiency of gas displacement by water-based mud during the process of invasion. Estimated values of porosity and permeability agree well with measurements performed on rock-core samples. In addition to estimating porosity and permeability, the petrophysical inversion algorithm provided accurate spatial distributions of gas saturation in the invaded rock formations that were not possible to obtain with conventional procedures based solely on the use of density and resistivity logs.

Geophysics ◽  
2006 ◽  
Vol 71 (4) ◽  
pp. F101-F119 ◽  
Author(s):  
Faruk O. Alpak ◽  
Carlos Torres-Verdín ◽  
Tarek M. Habashy

We have developed a new methodology for the quantitative petrophysical evaluation of borehole array-induction measurements. The methodology is based on the time evolution of the spatial distributions of fluid saturation and salt concentration attributed to mud-filtrate invasion. We use a rigorous formulation to account for the physics of fluid displacement in porous media resulting from water-base mud filtrate invading hydrocarbon-bearing rock formations. Borehole array-induction measurements are simulated in a coupled mode with the physics of fluid flow. We use inversion to estimate parametric 1D distributions of permeability and porosity that honor the measured array-induction logs. As a byproduct, the inversion yields 2D (axial-symmetric) spatial distributions of aqueous phase saturation, salt concentration, and electrical resistivity. We conduct numerical inversion experiments using noisy synthetic wireline logs. The inversion requires a priori knowledge of several mud, petrophys-ical, and fluid parameters. We perform a systematic study of the accuracy and reliability of the estimated values of porosity and permeability when knowledge of such parameters is uncertain. For the numerical cases considered in this paper, inversion results indicate that borehole electromagnetic-induction logs with multiple radial lengths of investigation (array-induction logs) enable the accurate and reliable estimation of layer-by-layer absolute permeability and porosity. The accuracy of the estimated values of porosity and permeability is higher than 95% in the presence of 5% measurement noise and 10% uncertainty in rock-fluid and mud parameters. However, for cases of deep invasion beyond the radial length of investigation of array-induction logging tools, the estimation of permeability becomes unreliable. We emphasize the importance of a sensitivity study prior to inversion to rule out potential biases in estimating permeability resulting from uncertain knowledge about rock-fluid and mud properties.


2020 ◽  
Vol 21 (3) ◽  
pp. 9-18
Author(s):  
Ahmed Abdulwahhab Suhail ◽  
Mohammed H. Hafiz ◽  
Fadhil S. Kadhim

   Petrophysical characterization is the most important stage in reservoir management. The main purpose of this study is to evaluate reservoir properties and lithological identification of Nahr Umar Formation in Nasiriya oil field. The available well logs are (sonic, density, neutron, gamma-ray, SP, and resistivity logs). The petrophysical parameters such as the volume of clay, porosity, permeability, water saturation, were computed and interpreted using IP4.4 software. The lithology prediction of Nahr Umar formation was carried out by sonic -density cross plot technique. Nahr Umar Formation was divided into five units based on well logs interpretation and petrophysical Analysis: Nu-1 to Nu-5. The formation lithology is mainly composed of sandstone interlaminated with shale according to the interpretation of density, sonic, and gamma-ray logs. Interpretation of formation lithology and petrophysical parameters shows that Nu-1 is characterized by low shale content with high porosity and low water saturation whereas Nu-2 and Nu-4 consist mainly of high laminated shale with low porosity and permeability. Nu-3 is high porosity and water saturation and Nu-5 consists mainly of limestone layer that represents the water zone.


Geophysics ◽  
2018 ◽  
Vol 83 (6) ◽  
pp. D187-D202 ◽  
Author(s):  
Elsa Maalouf ◽  
Carlos Torres-Verdín

Detecting vertical transversely isotropic (VTI) formations and quantifying the magnitude of anisotropy are fundamental for describing organic mudrocks. Methods used to estimate stiffness coefficients of VTI formations often provide discontinuous or spatially averaged results over depth intervals where formation layers are thinner than the receiver aperture of acoustic tools. We have developed an inversion-based method to estimate stiffness coefficients of VTI formations that are continuous over the examined depth interval and that are mitigated for spatial averaging effects. To estimate the coefficients, we use logs of frequency-dependent compressional, Stoneley, and quadrupole/flexural modes measured with wireline or logging-while-drilling (LWD) instruments in vertical wells penetrating horizontal layers. First, we calculate the axial sensitivity functions of borehole sonic modes to stiffness coefficients; next, we use the sensitivity functions to estimate the stiffness coefficients of VTI layers sequentially from frequency-dependent borehole sonic logs. Because sonic logs exhibit spatial averaging effects, we deaverage the logs by calculating layer-by-layer slownesses of formations prior to estimating stiffness coefficients. The method is verified with synthetic models of homogeneous and thinly bedded formations constructed from field examples of organic mudrocks. Results consist of layer-by-layer estimates of [Formula: see text], [Formula: see text], [Formula: see text], [Formula: see text], and [Formula: see text]. We observe three sources of error in the estimated coefficients: (1) bias error originating from deaveraging the sonic logs prior to the sequential inversion, (2) error propagated during the sequential inversion, and (3) error associated with noisy slowness logs. We found that the relative bias and uncertainty of the estimated coefficients are largest for [Formula: see text] and [Formula: see text] because borehole modes exhibit low sensitivity to these two coefficients. The main advantage of our method is that it mitigates spatial averaging effects of sonic logs, while at the same time it detects the presence of anisotropic layers and yields continuous estimations of stiffness coefficients along the depth interval of interest.


Geophysics ◽  
2002 ◽  
Vol 67 (6) ◽  
pp. 1753-1768 ◽  
Author(s):  
Yuji Mitsuhata ◽  
Toshihiro Uchida ◽  
Hiroshi Amano

Interpretation of controlled‐source electromagnetic (CSEM) data is usually based on 1‐D inversions, whereas data of direct current (dc) resistivity and magnetotelluric (MT) measurements are commonly interpreted by 2‐D inversions. We have developed an algorithm to invert frequency‐Domain vertical magnetic data generated by a grounded‐wire source for a 2‐D model of the earth—a so‐called 2.5‐D inversion. To stabilize the inversion, we adopt a smoothness constraint for the model parameters and adjust the regularization parameter objectively using a statistical criterion. A test using synthetic data from a realistic model reveals the insufficiency of only one source to recover an acceptable result. In contrast, the joint use of data generated by a left‐side source and a right‐side source dramatically improves the inversion result. We applied our inversion algorithm to a field data set, which was transformed from long‐offset transient electromagnetic (LOTEM) data acquired in a Japanese oil and gas field. As demonstrated by the synthetic data set, the inversion of the joint data set automatically converged and provided a better resultant model than that of the data generated by each source. In addition, our 2.5‐D inversion accounted for the reversals in the LOTEM measurements, which is impossible using 1‐D inversions. The shallow parts (above about 1 km depth) of the final model obtained by our 2.5‐D inversion agree well with those of a 2‐D inversion of MT data.


2011 ◽  
Vol 58-60 ◽  
pp. 1926-1931
Author(s):  
Fei Yu Lian ◽  
Qing Li

In this paper, we proposed a punctuated target identification method for underground homogeneous medium based on dielectric constant inversion algorithm. Its main idea is that, initially, Hough transform is used to locate the target characterized by hyperbola in a radar profile map, then a layer-by-layer waveform inversion algorithm is used to invert the dielectric constant on the location where the target lies. To guarantee the correctness of inversion, the transmission beam method is adopted to obtain initial parameters and calibration factors required by inversion. Through numerical analysis of the dielectric constant of the target, the classification of multiple targets in an underground medium can be confirmed. This method not only overcomes the failure of the traditional phase-comparison method to distinguish different kinds of targets in the same region, but also overcomes the limitation of the image processing method, in which it only classifies the target in a coarse-grained manner. Experimental results show that this method has many advantages, such as fine-grained classification, high precision, and multiple target identification, in the identification of underground targets.


2018 ◽  
Vol 58 (1) ◽  
pp. 339
Author(s):  
Brenton Richards ◽  
Alexander Côté

Over the past decade, there has been a paradigm shift in the exploitation strategy in North American tight gas plays from vertical to horizontal wells. This shift has yet to occur in Australia. The Cooper Basin has vast amounts of contingent and prospective tight gas resources that have yet to be unlocked commercially. These resources continue to be developed primarily with hydraulic fracture stimulated vertical wells. Operators have yet to challenge the status quo and test the Cooper Basin tight gas potential with a drilled, completed and tested horizontal well. There are many advantages to horizontal well developments, from the ability to target a specific high graded reservoir unit to increased capital efficiency. Operators need to break away from convention and take a new approach to Cooper Basin tight gas exploration and development in the quest to demonstrate commerciality. A review of the inherent challenges in Cooper Basin gas field developments and the current exploitation strategies employed in analogous tight gas plays have been integrated to produce a pragmatic workflow to identify potential reservoir units that would benefit from a change in development strategy.


2000 ◽  
Vol 40 (1) ◽  
pp. 355
Author(s):  
C.J. Shield

Water saturation (Sw) values calculated from resistivity or induction logs are often higher than those measured from core-derived capillary pressure (Pc) measurements. The core-derived Sw measurements are commonly applied for reservoir simulation in preference to the log-derived Sw calculations. As it is economically and logistically impractical to core every hydrocarbon reservoir, a method of correlating the core-derived Sw to resistivity/induction logs is required. Two-dimensional resistivity modelling is applied to dual laterolog data to ascertain the applicability of this technique.The Griffin and Scindian/Chinook Fields, offshore Western Australia, have been producing hydrocarbons since 1994 from two early-to-middle Cretaceous reservoirs, the clean quartzose sandstones of the Zeepaard Formation and the overlying glauconitic, quartzose sandstones of the Birdrong Formation. Routine and special core analysis of cores recovered from wells intersecting these two reservoirs creates an excellent data set with which to correlate the good quality wireline log data.A strong relationship is noted between the modelled water saturation from resistivity logs, and the irreducible water saturation measured from core capillary pressure data. Correlation between the core-derived permeability and the invasion diameter calculated from the modelled laterolog data is shown to produce a locally applicable means of estimating permeability from the resistivity modelling results.The evaluation of these data from the Griffin and Scindian/Chinook Fields provides a method for reducing appraisal and development well analysis costs, through the closer integration of core and wireline log data at an earlier stage of the field appraisal phase.


2021 ◽  
Vol 61 (2) ◽  
pp. 720
Author(s):  
Kasia Sobczak ◽  
Heinz-Gerd Holl ◽  
Andrew Garnett

The Upper Jurassic Walloon Coal Measures of the Surat Basin (Queensland) host some of the most prominent coal seam gas (CSG) resources in Australia. The Walloon Coal Measures are directly overlain by the Springbok Sandstone formation, historically referred to as a regional aquifer. An increasing number of studies and industry models suggest relatively limited hydraulic connectivity within the formation and between it and the underlying coal measures, due to extreme lithological heterogeneity. Accurate evaluation of the permeability, as well as lateral and vertical continuity of the lithological units within the Springbok Sandstone, is critical in reservoir models that form the basis of reasonable aquifer protection practices and impact prediction. This study presents a wireline log-based workflow applied to identify permeable zones within the Springbok Sandstone in 31 CSG wells across the Surat Basin that allows robust estimations of porosities and Klinkenberg permeabilities. The workflow primarily utilises spontaneous potential, density, neutron and resistivity logs, and was developed by integrating current industry practices implemented by operators on a local scale to identify risk (permeable) zones in the vicinity of targeted coal seams. The results of this case study indicate that permeable zones within the interval are volumetrically minor (on average 25% N/G) and likely isolated, with Klinkenberg permeabilities rarely exceeding 10–20mD. This evidence for low hydraulic connectivity, as well as significant local variations in the character of the Springbok Sandstone, suggests that the definition of the formation as a regional, continuous aquifer and the way it is modelled needs to be revised.


2013 ◽  
Vol 385-386 ◽  
pp. 474-477
Author(s):  
Yu Feng Lu ◽  
Min Yang ◽  
San Chuan Li ◽  
Bo Lv

Daqing Wuzhan gas field is a deformable low permeability gas filed. Based on the lab study of changing law of porosity and permeability with confining pressure, a numerical simulation model of deformable low permeability gas reservoir is established and is solved by using IMPES method. Comparing with the result of eclipse shows the method is valid. Field example shows the dynamic geological reserves of Wuzhan gas field is low and the new wells should be drilled at the area where there is production potential and larger effective thickness. Comparison of results obtained from different models with field data shows the result considering medium deformation is more reasonable.


Sign in / Sign up

Export Citation Format

Share Document