WO2013066682A1 - Inversion-based workflow for processing nuclear density images in high-angle and horizontal wells - Google Patents

Inversion-based workflow for processing nuclear density images in high-angle and horizontal wells Download PDF

Info

Publication number
WO2013066682A1
WO2013066682A1 PCT/US2012/061585 US2012061585W WO2013066682A1 WO 2013066682 A1 WO2013066682 A1 WO 2013066682A1 US 2012061585 W US2012061585 W US 2012061585W WO 2013066682 A1 WO2013066682 A1 WO 2013066682A1
Authority
WO
WIPO (PCT)
Prior art keywords
inversion
formation
density
model
borehole
Prior art date
Application number
PCT/US2012/061585
Other languages
French (fr)
Inventor
Sushil SHETTY
Dzevat Omeragic
Tarek M. Habashy
John C. Rasmus
Jeffrey MILES
Original Assignee
Schlumberger Canada Limited
Services Petroliers Schlumberger
Schlumberger Holdings Limited
Schlumberger Technology B.V.
Prad Research And Development Limited
Schlumberger Technology Corporation
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Schlumberger Canada Limited, Services Petroliers Schlumberger, Schlumberger Holdings Limited, Schlumberger Technology B.V., Prad Research And Development Limited, Schlumberger Technology Corporation filed Critical Schlumberger Canada Limited
Priority to GB1406621.1A priority Critical patent/GB2509450A/en
Priority to US14/352,566 priority patent/US9869792B2/en
Publication of WO2013066682A1 publication Critical patent/WO2013066682A1/en
Priority to NO20140524A priority patent/NO20140524A1/en
Priority to US15/845,515 priority patent/US10809416B2/en

Links

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V5/00—Prospecting or detecting by the use of ionising radiation, e.g. of natural or induced radioactivity
    • G01V5/04—Prospecting or detecting by the use of ionising radiation, e.g. of natural or induced radioactivity specially adapted for well-logging
    • G01V5/08—Prospecting or detecting by the use of ionising radiation, e.g. of natural or induced radioactivity specially adapted for well-logging using primary nuclear radiation sources or X-rays
    • G01V5/12—Prospecting or detecting by the use of ionising radiation, e.g. of natural or induced radioactivity specially adapted for well-logging using primary nuclear radiation sources or X-rays using gamma or X-ray sources
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V5/00—Prospecting or detecting by the use of ionising radiation, e.g. of natural or induced radioactivity
    • G01V5/04—Prospecting or detecting by the use of ionising radiation, e.g. of natural or induced radioactivity specially adapted for well-logging
    • G01V5/08—Prospecting or detecting by the use of ionising radiation, e.g. of natural or induced radioactivity specially adapted for well-logging using primary nuclear radiation sources or X-rays

Definitions

  • Nuclear density measurements progressed from designs focused on measuring an average density for vertical wells to a fully azimuthal measurement for inclined wells. While being a shallower measurement than the resistivity, the nuclear density is also affected by bed crossings and adjacent beds, and additionally by standoff and asymmetric ("teardrop") invasion.
  • Embodiments herein relate to apparatus and methods for characterizing a subterranean formation traversed by a wellbore including collecting data from the formation using a tool wherein the tool collects data to form an azimuthal image, characterizing a section of the formation comprising data and images acquired in a high angle wellbore section or horizontal wellbore section using a parametric model, and performing an inversion using apparent densities and volumetric photoelectric factor images to build a formation model wherein the inversion is tailored for high angle wellbore sections and/or horizontal wellbore sections.
  • Figure 1 is a computational grid for fast-forward modeling: 3D grid (left) and r- ⁇ grid (right), tied to the sensor, extending from -45° to 45°.
  • Figure 2 is a plot of LS nuclear density integrated sensitivity functions azimuthal dependence: radial (left) and axial (right).
  • Figure 3 is a plot of SS nuclear density integrated sensitivity functions azimuthal dependence: radial (left) and axial (right).
  • Figure 4 is a schematic drawing of a logging scenario used to derive model for sinusoidal features in COMP image for an "ideal" tool as described in the text.
  • the ideal tool measures the formation density within an infinitesimal volume-of-investigationS, located a radial distance EP L from the borehole wall.
  • Figure 5 is a COMP image corresponding to logging scenario shown in Fig. 4, for an idealized tool with infinitesimal volume-of-investigation and negligible borehole effects. Sensor azimuth ⁇ is referenced to top of hole.
  • Figure 6 are two plots of an EPL calibration for an 8.25 inch ECOSCOPE LWDTM tool.
  • Figure 7 is a schematic drawing of the high-angle well model parameterization, including set of dipping layers with trajectory of varying inclination, strictly up-section or down- section. Arbitrary borehole geometry is represented with dashed lines.
  • Figure 8 is a schematic drawing of the horizontal well model parameterization, including set of layers parallel with the tool at arbitrary azimuth.
  • the borehole is represented with dashed lines; its geometry is defined separately.
  • Figure 9 is a schematic drawing of a borehole model parameterization: (a) cylindrical geometry and (b) arbitrary geometry. Geometry in plane perpendicular to tool axis (upper panel) and plane containing tool axis (lower panel).
  • Figure 10 is a flow chart of an inversion-based workflow for processing data with a high-angle well model.
  • Figure 11 is an inversion-based workflow for processing data using a horizontal well model.
  • Figure 12 is a composite figure for high-angle well pre-processing for noisy synthetic data: bottom COMP squaring (red curve) produced a density profile (blue) (left); Sinusoid extraction from COMP image (dashed black curves) (middle); Adaptive segmentation (solid black lines) (right).
  • Figure 13 is a composite figure for high-angle inversion of noisy synthetic data-set: final layer profile (blue) and COMP bottom quadrant log (measured - red, reconstructed - green), (left); trajectory inclination, initial and final dip, and azimuth (right).
  • Figure 14 is a composite figure for reconstruction of COMP image for high-angle inversion of noisy synthetic data-set with relative dip variation 82° to 87°.
  • Figure 15 is a composite figure for reconstruction of LS image for high-angle inversion of noisy synthetic data-set with relative dip variation 82° to 87°.
  • Figure 16 is a composite figure for reconstruction of SS image for high-angle inversion of noisy synthetic data-set with relative dip variation 82° to 87°.
  • Figure 17 is a composite figure for borehole geometry (stand-off in inches) for high- angle inversion of noisy synthetic data: true geometry (left), initial geometry (middle), and final geometry from the inversion (right).
  • Figure 18 is a composite figure for reconstruction of COMP image for horizontal well inversion of noisy synthetic data-set with constant 89° relative dip.
  • Figure 19 is a composite figure for reconstruction of LS image for horizontal well inversion of noisy synthetic data-set with constant 89° relative dip.
  • Figure 20 is a composite figure for reconstruction of SS image for horizontal well inversion of noisy synthetic data-set with constant 89° relative dip.
  • Figure 21 is a composite figure for horizontal well inversion of noisy synthetic data (rel. dip 89°): Measured COMP image (top), final borehole-corrected image (middle), and final formation model in curtain-section plane (bottom). The green curves represent trajectory and the borehole surfaces.
  • Figure 22 is a reconstruction of COMP image for horizontal well inversion of noisy synthetic data-set with sinusoidal relative dip in the range 88°-92°.
  • Figure 23 is a horizontal well inversion of noisy synthetic data with variable rel. dip 88°-92°: Measured COMP image (top), final borehole-corrected image (middle), and final formation model in curtain-section plane (bottom). The green curves represent trajectory and the borehole surfaces.
  • Figure 24 is a composite figure illustrating high-angle well pre-processing of the field data: bottom quadrant COMP squaring (red curve) produced a density profile (blue) (left);
  • Figure 25 is a high-angle inversion of field data: Inverted density profile (blue), bottom quadrant COMP measurements (red) reconstructed (green); Trajectory inclination, initial and final apparent dip, and relative azimuth (right).
  • Figure 26 is a composite figure of reconstruction of COMP image for high-angle well inversion of interval from North Sea field data-set.
  • Figure 27 is a composite figure of reconstruction of LS images for high-angle well inversion of interval from North Sea field data-set.
  • Figure 28 is a composite figure of reconstruction of SS images for high-angle well inversion of interval from North Sea field data-set.
  • Figure 29 is a comparison plot of high-angle inversion of an interval from the North Sea field data-set: Average error per sector for each channel (left) Standard deviation of error per sector for each channel (right).
  • Figure 30 is a comparison plot of borehole geometry for high-angle well inversion of an interval from North Sea field data-set: density-caliper (left); initial geometry (middle), and inversion-based caliper (right) in inches.
  • Figure 31 is a composite figure of a high-angle inversion workflow on North Sea field data set: Measured COMP image (top), final borehole-corrected image (middle), and final formation model in curtain-section plane (bottom).
  • the black curves represent trajectory and the borehole surfaces.
  • Figure 32 is composite figure for high-angle well pre-processing for gas shale field data-set: bottom quadrant COMP squaring (red curve) produced a density profile (blue) (left); Sinusoid extraction from COMP image (dashed black curves) (middle); Adaptive segmentation (solid black lines) (right).
  • Figure 33 is a composite of a high-angle inversion of gas shale field data-set. Inverted density profile (blue), bottom quadrant COMP: measurements (red) reconstructed (green);
  • Figure 34 is a composite of a reconstruction of COMP image for high-angle well inversion of interval from a gas shale field data-set.
  • Figure 35 is composite of a reconstruction of LS images for high-angle well inversion of interval from a gas shale field data-set.
  • Figure 36 is a composite of a reconstruction of SS images for high-angle well inversion of interval from a gas field data-set.
  • Figure 37 is plot of a HA inversion of an interval from a gas shale field data-set:
  • Figure 38 is a comparison of borehole geometry for high-angle well inversion of an interval from a gas shale field data-set: density-caliper (left); initial geometry (middle), and inversion-based caliper (right) in inches.
  • Figure 39 is a composite figure of a HA inversion workflow on gas shale field data set: measured COMP image (top), final borehole-corrected image (middle) and final formation model in curtain-section plane (bottom).
  • the black solid curve is the tool trajectory and the dashed black curves are the borehole surfaces in the curtain-section plane.
  • Figure 40 is a composite figure for reconstruction of COMP image for horizontal well inversion of interval from North Sea field data-set.
  • Figure 41 is a plot of horizontal well inversion of an interval from North Sea field data- set: Average error per sector for each channel (left) Standard deviation of error per sector for each channel (right).
  • Figure 42 is a comparison of borehole geometry for horizontal well inversion of an interval from North Sea field data-set: density-caliper (left); initial geometry (middle), and inversion-based caliper (right) in inches.
  • Figure 43 is a composite plot of a final formation model in curtain-section plane reconstructed by HZ well workflow (bottom), corresponding measured COMP image and borehole-corrected image (middle).
  • composition used/disclosed herein can also comprise some components other than those cited.
  • model-based inversion incorporating the new second- order 3D sensitivity functions.
  • This inversion is capable of processing density images in HA as well as HZ wells, solving for a ID layered formation model, formation dip and azimuth, and 3D borehole geometry with locally homogeneous mud properties. Lateral variations in the model parameters are captured using an adaptive trajectory segmentation based on relative dip.
  • the inversion-based workflows are driven by measurement sensitivities, ensuring that the interpretation is fully consistent with the measurements from all sectors and channels.
  • the algorithm is primarily intended for 3D petrophysical interpretation in high-angle and horizontal wells, and for real-time well placement.
  • the methodology and model are primarily intended for 3D petrophysical interpretation in high-angle and horizontal wells, and for real-time well placement.
  • parameterization are fairly general and are compatible with models used for interpretation of other measurements (such as resistivity and directional EM) in HA, HZ wells.
  • a fast-forward model for the gamma-gamma density tool response was developed based on second-order 2D (axisymmetric) sensitivity functions.
  • the original 2D sensitivity functions have been expanded into a full 3D model of spatial sensitivity.
  • the simulation time is reduced from hours using MCNP to milliseconds for the sensitivity approach, making it applicable for commercial petrophysics.
  • Due to use of second-order 3D sensitivity functions the accuracy of fast- forward model is comparable to measurement accuracy in wide range of high- angle and horizontal well scenarios, enabling use in rigorous physics-based inversion workflows to reduce ambiguity and minimize subjectivity in the interpretation process.
  • the inversion-based approach promises to represent a step change in quantitative interpretation of nuclear measurements in HA/HZ wells, producing more accurate density profile, formation dip and azimuth, with accurate boundary positions and layer densities free of boundary and borehole effects, resulting in improved porosity and reserves estimates.
  • the parametric model used for interpretation includes a multi-layer dipping formation, mud properties, arbitrary borehole geometry, and 3D well trajectory. Lateral variations in formation and borehole parameters are captured by defining the models in discrete trajectory segments using an adaptive segmentation based on the local relative dip. Measurement sensitivities are used to design the flexible and robust inversion-based workflow for determining optimum parameter values from all the available measurements. A sliding window is used to enforce consistency of models between adjacent segments. The result of the inversion is the accurate layer thicknesses, shoulder-bed corrected layer densities, formation dip and azimuth in each segment. The inversion also produces a borehole-corrected formation density image and a robust caliper that takes into account the layered formation as background. The algorithm is primarily intended for 3D petrophysical interpretation in high-angle and horizontal wells. The methodology and model parameterization is general and compatible with models used for interpretation of other measurements (such as resistivity and directional EM), making feasible future integration in joint multi-physics interpretation.
  • Other measurements such as resistivity and directional
  • the horizontal well workflow (for relative dip 88°-92°) takes into account non-crossed layers and enables building the 2D curtain section model as superposition of ID formation models defined in very short segments along the trajectory, with layer boundaries in each segment oriented parallel to the trajectory.
  • the inversion uses a sliding window scheme to enforce model consistency between adjacent segments.
  • Inversion-based approach for estimating the borehole geometry based on two model parameterizations removes sinusoidal artifacts in conventional density caliper results, which are due to coupling of boundary and borehole effects.
  • the inversion-based caliper is more robust because it uses all the information available, and is more accurate because it models the true geometry as well as the layered background, effectively decoupling the borehole effect from boundary effect.
  • the ECOSCOPE LWDTM tool integrates a full suite of formation evaluation measurements surpassing the classic "triple-combo" service, with many of the measurements co- located and closer to the bit compared with the previous generation.
  • the tool has two density measurements: a gamma-gamma density from a
  • PNG pulsed neutron generator
  • Nuclear gamma-gamma density measurements are one of the key formation evaluation measurements and are used to infer formation porosity from the knowledge of lithology and pore fluids.
  • the gamma-gamma measurement uses a gamma-ray source and two gamma-ray detectors generating three independent measurements: a long-spacing (LS) and short- spacing (SS) apparent electron density and the apparent volumetric photoelectric factor (UAPP), used to infer the lithology.
  • the density measurement is resolved into 16 azimuthal sectors each spanning 22.5°.
  • the radial sensitivity extends radially approximately 6 inches into the formation, while the axial sensitivity covers 14-18 inches. Therefore, in HAHZ wells, the effective true vertical resolution will vary between these two factors depending on the relative angle between the trajectory and the formation layers.
  • the LS and SS based electron densities are combined in a "spine-and-ribs" algorithm to compute the compensated (COMP) density for each sector, quadrant and 16-sector average.
  • the compensation is also used along with the apparent mud density to compute a standoff distance which is added to the stabilizer size to create a 3D caliper, henceforth referred to as the density caliper.
  • the differential (DRHO) density defined as the difference between COMP and LS densities, is used as a quality indicator for the COMP measurement.
  • the fast-forward model uses precomputed first and second order sensitivity functions defined for the 3D grid shown in Fig. 1 to compute measurement responses for density and photoelectric factor properties defined on the same grid.
  • the radial grid extends out to 9 inches from the tool, with the cell size gradually increasing from 1/8 inch to 1 inch, total 22 cells.
  • the uniform grid in axial direction has 19 cells.
  • the spatial variation of integrated radial and axial sensitivity for LS and SS is shown in Figs. 2 and 3, respectively. Including the full 3D sensitivity is especially important for the LS channel.
  • the forward modeling is about million times faster than MCNP simulation used to generate them. It takes milliseconds to compute the responses for a single log-point and sector.
  • the model has representative accuracies of 0.02 g/cc, 0.05 g/cc, and 0.03 g/cc for LS, SS, and COMP, respectively.
  • the representative accuracy of the UAPP model is 1 b/cc for non-barite muds, while for heavy barite muds, the model is only qualitative. In configurations with standoff up to 0.25in, the accuracy of the model for all channels is comparable to the measurement accuracy, better than 0.02g/cc.
  • the error in approximation of COMP remains the same for standoffs up to linch, while it nearly doubles for LS over the same range.
  • the SS accuracy decrease by a factor of two for an intermediate standoff of 0.5 in.
  • the parametric model includes separate components for describing the formation and borehole along a given wellbore. Complete and general model parameterization enables easy modeling, sensitivity analysis and interfacing to the inversion. In general, model
  • parameterization should not be related to measurements, and is shared for interpretation of different measurements to build a common model.
  • Isocontours of the COMP image are computed using the "marching squares” algorithm, which is the 2D adaptation of the “marching cubes” algorithm for extracting isosurfaces from a 3D scalar field.
  • isocontours are computed for isovalues spanning the density range in the image, in increments of the measurement accuracy, 0.015 gfcc.
  • 3 ⁇ 4 is the location of the sinusoid in MD
  • a k is the sinusoid amplitude
  • k is the sinusoid phase.
  • Isocontours will not be exactly sinusoidal if the image contains noise, or if assumptions underlying the parametric model break down, for example, tool crossing a non- planar layer boundary, variable trajectory inclination and azimuth, or significant
  • the high-angle well model assumes that the trajectory can be discretized into segments such that the formation in each segment is ID layering with constant dip and azimuth.
  • the trajectory azimuth is assumed to be constant in each segment, and the trajectory inclination is assumed to be "down-section", with layers crossed at relative dip less than 90°, or "up-section", when layers are crossed at relative dip greater than 90°.
  • Fig. 7 shows the formation for a typical down-section segment when viewed in the "curtain-section" plane, a vertical surface defined by the trajectory.
  • the free parameters for the segment are the apparent formation dip ⁇ ⁇ , TVD location payers Q ⁇ i a y er b oun d ar y ? layer densities pj ayers , and layer photoelectric factors PEFj ayers .
  • the layer boundary locations may also be parameterized by their locations in MD, lj ayers , which is a more stable parameterization due to the fact that the initial values for l ayers are constrained from log-squaring, whereas the initial values of zj ayers are coupled to lj ayers via the initial value of the apparent formation dip ⁇ ⁇ .
  • the segment sizes, number of layers and relative azimuth are determined and fixed from the data in a pre-processing step.
  • the horizontal well model also assumes discrete segments such that the formation in each segment is ID layering with dip and azimuth, and such that the trajectory azimuth is constant in each segment.
  • the horizontal well model constrains the layer boundaries in each segment to be parallel to the trajectory, i.e., constant relative dip 90°.
  • Fig. 8 shows the model for a typical horizontal well segment when viewed in the curtain section.
  • the model parameters for each segment are the distance zj ayers in TVD from the first log-point to each layer boundary, layer densities p j layers , and layer photoelectric factors PEFj ayers .
  • the number of segments, segment sizes and number of layers J ⁇ ayers in each segment are determined and fixed from the data in a pre-processing step.
  • the horizontal well model is applicable to intervals with variable relative dip close to 90°, i.e., with borehole parallel or near parallel to the surrounding formation layers.
  • segments must be chosen sufficiently small (not below the LS axial resolution) so that the model error in each segment - due to variable trajectory inclination or layer boundaries not parallel to the trajectory - corresponds to change in responses smaller than the data noise.
  • the mud properties, p muc i and PEF mud are defined in the same trajectory segments used for the formation model, and are assumed to be uniform in each segment.
  • the geometry of the borehole is defined in terms of the radial distance from the tool axis to the borehole wall, denoted by r(9,MD), where ⁇ is the tool azimuth.
  • r(9,MD) the corresponding stand-off s(9,MD) is defined as where r too i is the tool radius.
  • r(9,MD) may be parameterized using one of the following models.
  • the radial distance to the borehole wall is assumed to be piecewise-constant function of both ⁇ and MD, defined by a discrete radius for each sector and each log-point.
  • the harmonic geometry can describe an arbitrary borehole surface.
  • the radial distance to the borehole wall is a harmonic function of ⁇ :
  • the high-angle pre-processing workflow consists of the following steps:
  • Log-Squaring - A log-squaring algorithm is applied to the bottom quadrant density (ROBB) log to construct a layered density profile along MD.
  • a photoelectric factor for each layer is obtained by sampling the UAPP bottom quadrant log at the layer mid-points.
  • Sinusoid Extraction The amplitude and phase of sinusoidal features corresponding to layer boundaries are extracted from the COMP image either manually, or by an automatic method. For each layer boundary, the extracted phase and amplitude are converted to relative dip and relative azimuth.
  • Trajectory Segmentation An adaptive algorithm is used to discretize the trajectory such that the COMP image for each segment includes at least half a sinusoidal feature, which improves sensitivity to dip.
  • the algorithm also constrains segment end-points to coincide with layer mid-points, so that each segment includes at least half of a layer, improving sensitivity to properties of truncated layers.
  • Initial Model The output from log-squaring and sinusoid extraction is used to define the number of layers, their properties and boundaries, apparent formation dip, and relative azimuth for each segment.
  • the layer boundaries from log-squaring must be shifted in MD in order to obtain their locations referenced to the tool axis.
  • the initial borehole geometry is defined as a cylinder with diameter equal to the bit-size. Alternatively, the initial borehole geometry may also be defined using the density caliper.
  • the initial mud density and mud photoelectric factor are specified by the user.
  • a subset of model parameters is optimized from a subset of channels which are selected based on their sensitivity to the parameters being optimized, resulting in a more robust workflow.
  • Gauss-Newton optimization with line search, adaptive regularization and parameter constraints is used to minimize the cost function, defined to be the weighted Z 2 -norm error between the measurements and the modeled responses over all log-points in a processing window (see below for definition of processing window).
  • the weights in the cost function are the reciprocal of the measurement noise level for the channel used (0.015 g/cc for SS, LS, COMP, and 0.5 b/e- for UAPP). Measurements from all sectors are included with equal weights; no special weighting is applied to sectors from the bottom of the hole.
  • the model parameters are inverted from measurements in a processing window that includes all log- points in the segment, as well as log-points from buffer regions on each side of the segment. Formation layers from the previous and subsequent segments that influence measurements in the buffer regions are added to the formation model in each segment. The properties and thicknesses of the layers included from the previous segment are held fixed during the inversion.
  • the purpose of the post-processing workflow is to visualize the inversion results, and to compute various measures of the quality of the results.
  • the formation model is visualized as a sequence of ID multi-layer dipping formations in the curtain-section plane or as the inversion- based borehole-corrected image, which is obtained by projecting the layer properties in each segment onto a cylinder with radius r too/ +EPL, where EPL is the effective penetration length of the COMP channel.
  • EPL is the effective penetration length of the COMP channel.
  • the uncertainties in the final model parameters for each segment may be estimated from the model covariance matrix.
  • Fig. 11 The steps for processing a horizontal interval are shown in Fig. 11 and are grouped in a similar manner to the high-angle workflow, with a pre-processing workflow, an inversion workflow, and a post-processing workflow.
  • segment sizes must be chosen so as to balance the opposing requirements of a small window size for minimizing model error, and a large window size for improving inversion robustness. Based on tests with perfect synthetic data described earlier, we use segment size of 2 ft in MD.
  • the initial parameters for the borehole model are defined using the same procedure described for the high-angle pre-processing workflow.
  • the density images are processed segment by segment using a sliding window scheme analogous to the one used for the high-angle well workflow described earlier.
  • the 8.25 inch tool size was modeled over the 200 ft trajectory section with inclination varying linearly from 77° to 82°.
  • Gaussian noise with zero mean and standard deviation of 0.02 g/cc which is comparable to the amplitude of real measurement noise, was added to the synthetic measurements for all channels.
  • the initial mud properties were set to their true values.
  • Fig. 12 The output from the pre-processing workflow is shown in Fig. 12.
  • Log-squaring and sinusoid extraction are used to build an initial segmented multi-layer formation model.
  • the initial densities in the thin layers have errors due to shoulder-bed effect.
  • the inverted, fully shoulder-bed corrected profile is shown in Fig 13.
  • Fig. 14 shows the measured, initial, and final COMP density images, as well as reconstruction errors for the initial guess and the solution.
  • the initial error shows large systematic discrepancies due to errors in the initial dip, initial layer densities and boundary locations obtained from log-squaring and sinusoid extraction. These discrepancies are corrected by the inversion.
  • Figs. 15 and 16 show the measured, initial, and final images for the LS and SS densities, and corresponding error images. The initial error shows similar systematic
  • Fig. 17 shows the standoff from the inversion.
  • the final errors appear random for all channels and the mean and standard deviation of the error for each sector.
  • the errors are comparable to the measurement noise, which confirms the accuracy of the inversion for this data-set.
  • Fig. 19 shows the borehole-corrected image and final formation model in the curtain-section view.
  • Fig. 24 shows the output from pre-processing, as a result of which the trajectory has been discretized into segments with a ID multi-layer formation in each segment.
  • the initial borehole was assumed to be a cylinder with diameter equal to the bit- size of 8.5 inch, and the initial mud density and mud PEF were set to 1.4 g/cc and 10 b/e-, respectively.
  • Fig. 25 shows the reconstructed COMP bottom quadrant log, and final layer profile, dip, and azimuth. Note the shoulder-bed corrections in thin layers and corrections to the initial dip of approximately 1°.
  • Figs. 26-28 show the reconstructed images for the COMP, LS and SS density channels, as well as the corresponding errors.
  • COMP image reconstruction shows that errors are reduced in all quadrants due to optimization of the layer properties, thicknesses, and dip. Final errors are less than 5 %, and are larger in the upper quadrant relative to the bottom quadrant, possibly due to inaccuracies in the mud properties or lower accuracy of fast- forward model for the large standoff.
  • LS and SS reconstruction is similar, but in addition, significant reduction in the initial error in the upper quadrant is seen due to optimization of the borehole geometry.
  • Fig. 29 shows the average and standard deviation of the reconstruction error for each sector.
  • the final average error is less than 0.03 g/cc for all channels over all sectors, with close to zero average error in the bottom quadrant.
  • the final standard deviation of the error is less than 0.03 g/cc for all channels in all sectors.
  • Fig. 30 shows the final borehole geometry, which was parameterized using a harmonic model with two harmonics in the expansion. The same figure also shows the geometry from the density caliper. The density caliper has sinusoidal artifacts whose locations correlate with layers in the COMP image. These artifacts are not present in the geometry from the inversion.
  • Fig. 31 shows the final borehole-corrected image and final formation model in the curtain-section view. Field data-set from a gas shale well
  • Fig. 32 shows the output from preprocessing, where the trajectory has been discretized into segments with associated initial ID multi-layer formation in each segment.
  • the initial borehole geometry was assumed to be a cylinder with diameter equal to the bit-size (8.5 inch), and the mud density and mud PEF were set to 1.4 g/cc and 50 b/e-, respectively.
  • Fig. 33 shows the reconstructed layered density profile, dip, and azimuth.
  • the reconstructed dip is corrected up to 2° compared to initial guess.
  • Fig. 34 shows the measured, initial, and final images for the COMP channel, as well as the initial and final errors.
  • Fig. 35-36 show the measured, initial, and final images for the LS and SS channels, as well as the initial and final errors.
  • Fig. 37 shows the average and standard deviation of the initial and final errors for each sector. The inversion significantly reduces the mismatch between modeling and the measurements.
  • Fig. 38 shows the final borehole geometry.
  • Fig. 39 shows the final borehole-corrected image and final formation model in the curtain- section plane.
  • Fig. 40 The measured COMP image, the initial guess, inversion reconstructed image and corresponding misfit are shown in Fig. 40.
  • Fig. 41 shows the average and standard deviation of the initial and final errors for each sector.
  • the inversion significantly reduces the mismatch between modeling and the measurements.
  • Borehole geometry in terms of standoff is shown in Fig 42.
  • Initial geometry was cylinder with 8.5 inch diameter. Harmonic parameterization with two harmonics in the expansion was used.
  • Figure 43 bottom track shows the final formation model in the curtain-section plane.
  • the green curve is the tool trajectory and the dashed green curves are the borehole surfaces.

Landscapes

  • Physics & Mathematics (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Geophysics (AREA)
  • Geophysics And Detection Of Objects (AREA)
  • Image Processing (AREA)

Abstract

Methods and apparatus for characterizing a subterranean formation traversed by a wellbore including collecting data from the formation using a tool wherein the tool collects data to form an azimuthal image, characterizing a section of the formation comprising data and images acquired in a high angle wellbore section or horizontal wellbore section using a parametric model, and performing an inversion using apparent densities and volumetric photoelectric factor images to build a formation model wherein the inversion is tailored for high angle wellbore sections and/or horizontal wellbore sections.

Description

INVERSION-BASED WORKFLOW FOR PROCESSING NUCLEAR DENSITY IMAGES IN HIGH-ANGLE AND HORIZONTAL WELLS
PRIORITY
[0001] This application claims priority to United States Provisional Patent Application Serial Number 61/553810, filed October 31, 2011, which is incorporated by reference herein.
BACKGROUND
[0002] Formation evaluation in high-angle and horizontal (HA, HZ) wells is still a challenge due to measurement complexity, environmental effects and lack of adequate interpretation answer products. The effects of bed boundary crossings at high angles, influence of adjacent bed above and below the wellbore, as well as anisotropy were recognized on all the measurements. In case of resistivity interpretation, these effects made it impossible for the petrophysicists to use any of resistivity measurements as a direct representation of Rt over the intervals where the effects are present. It is necessary to model the tool response in order to determine the individual layer Rt, or to compare directly the offset well values for correlation and infer an Rt for petrophysics.
Development of directional measurements significantly increased the potential to accurately characterize formation structure near the wellbore, correct the effects of adjacent boundaries and made possible use of inversion to build more accurate 2D or 3D formation models.
[0003] Nuclear density measurements progressed from designs focused on measuring an average density for vertical wells to a fully azimuthal measurement for inclined wells. While being a shallower measurement than the resistivity, the nuclear density is also affected by bed crossings and adjacent beds, and additionally by standoff and asymmetric ("teardrop") invasion.
[0004] Traditionally, modeling of the nuclear tool responses relies on Monte Carlo Nuclear Particle (MCNP) modeling codes that are computationally intensive and not suitable for log analysis. They are primarily used in the tool design process to optimize the sensitivity and accuracy of the measurement and to develop cross-plots, correction charts and case studies to better understand responses in complex scenarios. [0005] The first quantitative inversion-based interpretation of nuclear density images was developed using a fast- forward model based on linear approximation and 2D flux sensitivity functions. The methodology ignores the borehole and determines the initial formation model using only the bottom quadrant data where the standoff is generally negligible. Based on the dip extracted from image sinusoids, the technique is able to build a ID layered earth model of formation density in high-angle wells. However, its applicability is limited due to the limitations of 2D (axisymmetric) assumption, and the loss of information from not using data from all sectors. The same constraint severely limits the applicability to horizontal wells, where the image usually does not contain sinusoids, and use of bottom sector data is not sufficient in automated interpretation workflows due to variation of dip and effects of adjacent (non-crossed) layers.
SUMMARY
[0006] Embodiments herein relate to apparatus and methods for characterizing a subterranean formation traversed by a wellbore including collecting data from the formation using a tool wherein the tool collects data to form an azimuthal image, characterizing a section of the formation comprising data and images acquired in a high angle wellbore section or horizontal wellbore section using a parametric model, and performing an inversion using apparent densities and volumetric photoelectric factor images to build a formation model wherein the inversion is tailored for high angle wellbore sections and/or horizontal wellbore sections.
FIGURES
[0007] Figure 1 is a computational grid for fast-forward modeling: 3D grid (left) and r-θ grid (right), tied to the sensor, extending from -45° to 45°.
[0008] Figure 2 is a plot of LS nuclear density integrated sensitivity functions azimuthal dependence: radial (left) and axial (right).
[0009] Figure 3 is a plot of SS nuclear density integrated sensitivity functions azimuthal dependence: radial (left) and axial (right).
[00010] Figure 4 is a schematic drawing of a logging scenario used to derive model for sinusoidal features in COMP image for an "ideal" tool as described in the text. At a given log- point 0, and sensor azimuth(), the ideal tool measures the formation density within an infinitesimal volume-of-investigationS, located a radial distance EP L from the borehole wall.
[00011] Figure 5 is a COMP image corresponding to logging scenario shown in Fig. 4, for an idealized tool with infinitesimal volume-of-investigation and negligible borehole effects. Sensor azimuth Θ is referenced to top of hole.
[00012] Figure 6 are two plots of an EPL calibration for an 8.25 inch ECOSCOPE LWD™ tool.
[00013] Figure 7 is a schematic drawing of the high-angle well model parameterization, including set of dipping layers with trajectory of varying inclination, strictly up-section or down- section. Arbitrary borehole geometry is represented with dashed lines.
[00014] Figure 8 is a schematic drawing of the horizontal well model parameterization, including set of layers parallel with the tool at arbitrary azimuth. The borehole is represented with dashed lines; its geometry is defined separately.
[00015] Figure 9 is a schematic drawing of a borehole model parameterization: (a) cylindrical geometry and (b) arbitrary geometry. Geometry in plane perpendicular to tool axis (upper panel) and plane containing tool axis (lower panel).
[00016] Figure 10 is a flow chart of an inversion-based workflow for processing data with a high-angle well model.
[00017] Figure 11 is an inversion-based workflow for processing data using a horizontal well model.
[00018] Figure 12 is a composite figure for high-angle well pre-processing for noisy synthetic data: bottom COMP squaring (red curve) produced a density profile (blue) (left); Sinusoid extraction from COMP image (dashed black curves) (middle); Adaptive segmentation (solid black lines) (right).
[00019] Figure 13 is a composite figure for high-angle inversion of noisy synthetic data-set: final layer profile (blue) and COMP bottom quadrant log (measured - red, reconstructed - green), (left); trajectory inclination, initial and final dip, and azimuth (right).
[00020] Figure 14 is a composite figure for reconstruction of COMP image for high-angle inversion of noisy synthetic data-set with relative dip variation 82° to 87°. [00021] Figure 15 is a composite figure for reconstruction of LS image for high-angle inversion of noisy synthetic data-set with relative dip variation 82° to 87°.
[00022] Figure 16 is a composite figure for reconstruction of SS image for high-angle inversion of noisy synthetic data-set with relative dip variation 82° to 87°.
[00023] Figure 17 is a composite figure for borehole geometry (stand-off in inches) for high- angle inversion of noisy synthetic data: true geometry (left), initial geometry (middle), and final geometry from the inversion (right).
[00024] Figure 18 is a composite figure for reconstruction of COMP image for horizontal well inversion of noisy synthetic data-set with constant 89° relative dip.
[00025] Figure 19 is a composite figure for reconstruction of LS image for horizontal well inversion of noisy synthetic data-set with constant 89° relative dip.
[00026] Figure 20 is a composite figure for reconstruction of SS image for horizontal well inversion of noisy synthetic data-set with constant 89° relative dip.
[00027] Figure 21 is a composite figure for horizontal well inversion of noisy synthetic data (rel. dip 89°): Measured COMP image (top), final borehole-corrected image (middle), and final formation model in curtain-section plane (bottom). The green curves represent trajectory and the borehole surfaces.
[00028] Figure 22 is a reconstruction of COMP image for horizontal well inversion of noisy synthetic data-set with sinusoidal relative dip in the range 88°-92°.
[00029] Figure 23 is a horizontal well inversion of noisy synthetic data with variable rel. dip 88°-92°: Measured COMP image (top), final borehole-corrected image (middle), and final formation model in curtain-section plane (bottom). The green curves represent trajectory and the borehole surfaces.
[00030] Figure 24 is a composite figure illustrating high-angle well pre-processing of the field data: bottom quadrant COMP squaring (red curve) produced a density profile (blue) (left);
Sinusoid extraction from COMP image (dashed black curves) (middle); Adaptive segmentation (solid black lines) (right). [00031] Figure 25 is a high-angle inversion of field data: Inverted density profile (blue), bottom quadrant COMP measurements (red) reconstructed (green); Trajectory inclination, initial and final apparent dip, and relative azimuth (right).
[00032] Figure 26 is a composite figure of reconstruction of COMP image for high-angle well inversion of interval from North Sea field data-set.
[00033] Figure 27 is a composite figure of reconstruction of LS images for high-angle well inversion of interval from North Sea field data-set.
[00034] Figure 28 is a composite figure of reconstruction of SS images for high-angle well inversion of interval from North Sea field data-set.
[00035] Figure 29 is a comparison plot of high-angle inversion of an interval from the North Sea field data-set: Average error per sector for each channel (left) Standard deviation of error per sector for each channel (right).
[00036] Figure 30 is a comparison plot of borehole geometry for high-angle well inversion of an interval from North Sea field data-set: density-caliper (left); initial geometry (middle), and inversion-based caliper (right) in inches.
[00037] Figure 31 is a composite figure of a high-angle inversion workflow on North Sea field data set: Measured COMP image (top), final borehole-corrected image (middle), and final formation model in curtain-section plane (bottom). The black curves represent trajectory and the borehole surfaces.
[00038] Figure 32 is composite figure for high-angle well pre-processing for gas shale field data-set: bottom quadrant COMP squaring (red curve) produced a density profile (blue) (left); Sinusoid extraction from COMP image (dashed black curves) (middle); Adaptive segmentation (solid black lines) (right).
[00039] Figure 33 is a composite of a high-angle inversion of gas shale field data-set. Inverted density profile (blue), bottom quadrant COMP: measurements (red) reconstructed (green);
Trajectory inclination, initial and final apparent dip, and relative azimuth (right).
[00040] Figure 34 is a composite of a reconstruction of COMP image for high-angle well inversion of interval from a gas shale field data-set. [00041] Figure 35 is composite of a reconstruction of LS images for high-angle well inversion of interval from a gas shale field data-set.
[00042] Figure 36 is a composite of a reconstruction of SS images for high-angle well inversion of interval from a gas field data-set.
[00043] Figure 37 is plot of a HA inversion of an interval from a gas shale field data-set:
Average error per sector for each channel (left) Standard deviation of error per sector for each channel (right).
[00044] Figure 38 is a comparison of borehole geometry for high-angle well inversion of an interval from a gas shale field data-set: density-caliper (left); initial geometry (middle), and inversion-based caliper (right) in inches.
[00045] Figure 39 is a composite figure of a HA inversion workflow on gas shale field data set: measured COMP image (top), final borehole-corrected image (middle) and final formation model in curtain-section plane (bottom). The black solid curve is the tool trajectory and the dashed black curves are the borehole surfaces in the curtain-section plane.
[00046] Figure 40 is a composite figure for reconstruction of COMP image for horizontal well inversion of interval from North Sea field data-set.
[00047] Figure 41 is a plot of horizontal well inversion of an interval from North Sea field data- set: Average error per sector for each channel (left) Standard deviation of error per sector for each channel (right).
[00048] Figure 42 is a comparison of borehole geometry for horizontal well inversion of an interval from North Sea field data-set: density-caliper (left); initial geometry (middle), and inversion-based caliper (right) in inches.
[00049] Figure 43 is a composite plot of a final formation model in curtain-section plane reconstructed by HZ well workflow (bottom), corresponding measured COMP image and borehole-corrected image (middle). DETAILED DESCRIPTION
[00050] At the outset, it should be noted that in the development of any such actual
embodiment, numerous implementation— specific decisions must be made to achieve the developer's specific goals, such as compliance with system related and business related constraints, which will vary from one implementation to another. Moreover, it will be appreciated that such a development effort might be complex and time consuming but would nevertheless be a routine undertaking for those of ordinary skill in the art having the benefit of this disclosure. In addition, the composition used/disclosed herein can also comprise some components other than those cited. In the summary of the invention and this detailed
description, each numerical value should be read once as modified by the term "about" (unless already expressly so modified), and then read again as not so modified unless otherwise indicated in context. Also, in the summary of the invention and this detailed description, it should be understood that a concentration range listed or described as being useful, suitable, or the like, is intended that any and every concentration within the range, including the end points, is to be considered as having been stated. For example, "a range of from 1 to 10" is to be read as indicating each and every possible number along the continuum between about 1 and about 10. Thus, even if specific data points within the range, or even no data points within the range, are explicitly identified or refer to only a few specific, it is to be understood that inventors appreciate and understand that any and all data points within the range are to be considered to have been specified, and that inventors possessed knowledge of the entire range and all points within the range.
[00051] The statements made herein merely provide information related to the present disclosure and may not constitute prior art, and may describe some embodiments illustrating the invention.
[00052] In this application, we present a model-based inversion incorporating the new second- order 3D sensitivity functions. This inversion is capable of processing density images in HA as well as HZ wells, solving for a ID layered formation model, formation dip and azimuth, and 3D borehole geometry with locally homogeneous mud properties. Lateral variations in the model parameters are captured using an adaptive trajectory segmentation based on relative dip. The inversion-based workflows are driven by measurement sensitivities, ensuring that the interpretation is fully consistent with the measurements from all sectors and channels.
[00053] The algorithm is primarily intended for 3D petrophysical interpretation in high-angle and horizontal wells, and for real-time well placement. The methodology and model
parameterization are fairly general and are compatible with models used for interpretation of other measurements (such as resistivity and directional EM) in HA, HZ wells.
[00054] A fast-forward model for the gamma-gamma density tool response was developed based on second-order 2D (axisymmetric) sensitivity functions. Here, the original 2D sensitivity functions have been expanded into a full 3D model of spatial sensitivity. The simulation time is reduced from hours using MCNP to milliseconds for the sensitivity approach, making it applicable for commercial petrophysics. Due to use of second-order 3D sensitivity functions, the accuracy of fast- forward model is comparable to measurement accuracy in wide range of high- angle and horizontal well scenarios, enabling use in rigorous physics-based inversion workflows to reduce ambiguity and minimize subjectivity in the interpretation process. The inversion-based approach promises to represent a step change in quantitative interpretation of nuclear measurements in HA/HZ wells, producing more accurate density profile, formation dip and azimuth, with accurate boundary positions and layer densities free of boundary and borehole effects, resulting in improved porosity and reserves estimates.
[00055] Workflows are developed for the interpretation of logging-while-drilling density images in high-angle and horizontal wells. The key component of the workflow is the 3D parametric inversion using a robust Gauss-Newton optimization engine and a new fast-forward model based on second-order 3D sensitivity functions.
[00056] The parametric model used for interpretation includes a multi-layer dipping formation, mud properties, arbitrary borehole geometry, and 3D well trajectory. Lateral variations in formation and borehole parameters are captured by defining the models in discrete trajectory segments using an adaptive segmentation based on the local relative dip. Measurement sensitivities are used to design the flexible and robust inversion-based workflow for determining optimum parameter values from all the available measurements. A sliding window is used to enforce consistency of models between adjacent segments. The result of the inversion is the accurate layer thicknesses, shoulder-bed corrected layer densities, formation dip and azimuth in each segment. The inversion also produces a borehole-corrected formation density image and a robust caliper that takes into account the layered formation as background. The algorithm is primarily intended for 3D petrophysical interpretation in high-angle and horizontal wells. The methodology and model parameterization is general and compatible with models used for interpretation of other measurements (such as resistivity and directional EM), making feasible future integration in joint multi-physics interpretation.
[00057] The workflow is tuned for scenarios with wellbore trajectory nearly parallel to layer boundaries, where ambiguity in interpretation is increased due to difficulty in determining dip, lateral changes in layer properties and influence of stand-off and nearby non-crossed boundaries.
[00058] The methodology is validated on synthetic thin layer models with variable dip and azimuth. Several examples are presented showing successful application of workflow to field data. The inversion-derived caliper is compared to the conventional density caliper to illustrate the improvement.
[00059] We include in the analysis layered formation models simultaneously with the 3D borehole geometry to better capture variable standoff and properties of thin layers, essential for improved quantitative interpretation especially in horizontal wells. Variations in formation parameters along the trajectory in high angle wells (with relative dip up to 88°) were captured by defining the ID formation models in discrete trajectory segments using an adaptive segmentation algorithm based on image sinusoids so as to improve inversion robustness in each segment.
[00060] The horizontal well workflow (for relative dip 88°-92°) takes into account non-crossed layers and enables building the 2D curtain section model as superposition of ID formation models defined in very short segments along the trajectory, with layer boundaries in each segment oriented parallel to the trajectory. For both high-angle and horizontal well workflows, the inversion uses a sliding window scheme to enforce model consistency between adjacent segments.
[00061] Inversion-based approach for estimating the borehole geometry based on two model parameterizations removes sinusoidal artifacts in conventional density caliper results, which are due to coupling of boundary and borehole effects. The inversion-based caliper is more robust because it uses all the information available, and is more accurate because it models the true geometry as well as the layered background, effectively decoupling the borehole effect from boundary effect.
[00062] Inversion results have been validated using noisy synthetic data and have been successfully applied to several field data sets. Results for field data-sets show good consistency of the final formation model from one segment to the next. Measurement reconstruction for all channels is very good and the final reconstruction errors are typically comparable to the measurement accuracy.
NUCLEAR DENSITY MEASUREMENTS
[00063] The ECOSCOPE LWD™ tool integrates a full suite of formation evaluation measurements surpassing the classic "triple-combo" service, with many of the measurements co- located and closer to the bit compared with the previous generation. In addition to propagation resistivities, the tool has two density measurements: a gamma-gamma density from a
conventional chemical source and a neutron-gamma density derived from the pulsed neutron generator (PNG). The PNG also enables the spectroscopy, neutron, and sigma measurements.
[00064] Nuclear gamma-gamma density measurements are one of the key formation evaluation measurements and are used to infer formation porosity from the knowledge of lithology and pore fluids. The gamma-gamma measurement uses a gamma-ray source and two gamma-ray detectors generating three independent measurements: a long-spacing (LS) and short- spacing (SS) apparent electron density and the apparent volumetric photoelectric factor (UAPP), used to infer the lithology. The density measurement is resolved into 16 azimuthal sectors each spanning 22.5°. The radial sensitivity extends radially approximately 6 inches into the formation, while the axial sensitivity covers 14-18 inches. Therefore, in HAHZ wells, the effective true vertical resolution will vary between these two factors depending on the relative angle between the trajectory and the formation layers.
[00065] The LS and SS based electron densities are combined in a "spine-and-ribs" algorithm to compute the compensated (COMP) density for each sector, quadrant and 16-sector average. The compensation is also used along with the apparent mud density to compute a standoff distance which is added to the stabilizer size to create a 3D caliper, henceforth referred to as the density caliper. The differential (DRHO) density, defined as the difference between COMP and LS densities, is used as a quality indicator for the COMP measurement.
[00066] The images are used to identify bed boundaries and estimate dip, enabling optimal well placement, steering the well to more productive zones. Normally the bottom quadrant density is used for petrophysical interpretations. However, in HAHZ wells their use is not straightforward since they often need to be combined with non-azimuthal gamma-ray and resistivity
measurements, which, unlike the bottom quadrant density response, are also affected by layers to the side and above the borehole. Although shallow, the density response can be affected when approaching bed boundaries.
Fast forward modeling of nuclear density
[00067] The fast-forward model uses precomputed first and second order sensitivity functions defined for the 3D grid shown in Fig. 1 to compute measurement responses for density and photoelectric factor properties defined on the same grid. There are 5 cells in the azimuthal direction, with grid nodes located at ± 5°, ±15° and ± 45° with respect to the detector face normal. The radial grid extends out to 9 inches from the tool, with the cell size gradually increasing from 1/8 inch to 1 inch, total 22 cells. The uniform grid in axial direction has 19 cells. The spatial variation of integrated radial and axial sensitivity for LS and SS is shown in Figs. 2 and 3, respectively. Including the full 3D sensitivity is especially important for the LS channel.
[00068] The forward modeling is about million times faster than MCNP simulation used to generate them. It takes milliseconds to compute the responses for a single log-point and sector. The model has representative accuracies of 0.02 g/cc, 0.05 g/cc, and 0.03 g/cc for LS, SS, and COMP, respectively. The representative accuracy of the UAPP model is 1 b/cc for non-barite muds, while for heavy barite muds, the model is only qualitative. In configurations with standoff up to 0.25in, the accuracy of the model for all channels is comparable to the measurement accuracy, better than 0.02g/cc. The error in approximation of COMP remains the same for standoffs up to linch, while it nearly doubles for LS over the same range. The SS accuracy decrease by a factor of two for an intermediate standoff of 0.5 in. MODEL PARAMETERIZATION
[00069] The parametric model includes separate components for describing the formation and borehole along a given wellbore. Complete and general model parameterization enables easy modeling, sensitivity analysis and interfacing to the inversion. In general, model
parameterization should not be related to measurements, and is shared for interpretation of different measurements to build a common model.
[00070] We consider separately the high angle well scenario from the horizontal well scenario. In the first, the tool crosses multiple boundaries, while in the latter the tool is nearly parallel to formation boundaries and may include non-crossed boundaries. The models are compatible and could be easily merged. In both scenarios the borehole models are the same, defined in discrete trajectory segments. To decide which scenario best applies to a set of data, the COMP image is used to form a parametric model for sinusoidal features.
Methodology for Automatic Dip Estimation
[00071] Although the parametric model was derived for the case of the tool crossing a single planar boundary, it should be clear that if the tool crosses a sequence of planar boundaries, then the corresponding COMP image will be characterized by a sequence of sinusoidal contours, each with a different amplitude and phase depending on the local orientation of the layer boundary relative to the trajectory. Also, for a real tool with finite volume of investigation, we expect the sharp density jump in Fig. 5 to be smeared out over a thin "transition zone" near each layer boundary, with zone thickness dependent on the dimensions of the volume-of-investigation and on the physics of the measurement. Assuming that isocontours of the COMP image in each transition zone are described by the sinusoidal model we can estimate dip and azimuth of layer boundaries.
Sinusoid Extraction
[00072] Isocontours of the COMP image are computed using the "marching squares" algorithm, which is the 2D adaptation of the "marching cubes" algorithm for extracting isosurfaces from a 3D scalar field. To ensure that isocontours are extracted uniformly from the image, isocontours are computed for isovalues spanning the density range in the image, in increments of the measurement accuracy, 0.015 gfcc. Phase and amplitude for the extracted contours are determined using the least-squares method, by fitting each contour, denoted y (Θ) for £ = 1, N1 , to a sinusoidal model of the form γ άε1(θ)≡
AiCos(0 - + Ck. Here ¾ is the location of the sinusoid in MD, Ak is the sinusoid amplitude, and k is the sinusoid phase.
[00073] Isocontours will not be exactly sinusoidal if the image contains noise, or if assumptions underlying the parametric model break down, for example, tool crossing a non- planar layer boundary, variable trajectory inclination and azimuth, or significant
uncompensated borehole effects.
High-angle well model
[00074] The high-angle well model assumes that the trajectory can be discretized into segments such that the formation in each segment is ID layering with constant dip and azimuth. The trajectory azimuth is assumed to be constant in each segment, and the trajectory inclination is assumed to be "down-section", with layers crossed at relative dip less than 90°, or "up-section", when layers are crossed at relative dip greater than 90°. Fig. 7 shows the formation for a typical down-section segment when viewed in the "curtain-section" plane, a vertical surface defined by the trajectory.
[00075] The free parameters for the segment are the apparent formation dip βα, TVD location payers Q^ iayer boundary? layer densities pjayers, and layer photoelectric factors PEFjayers . The layer boundary locations may also be parameterized by their locations in MD, ljayers, which is a more stable parameterization due to the fact that the initial values for l ayers are constrained from log-squaring, whereas the initial values of zjayers are coupled to ljayers via the initial value of the apparent formation dip βα. The segment sizes, number of layers
Figure imgf000015_0001
and relative azimuth are determined and fixed from the data in a pre-processing step.
[00076] The segmented nature of the high-angle model allows it to capture variations in dip and layer properties along the trajectory. The efficiency of the high-angle parameterization degrades for relative dip close to 90°, where a separate approach is needed.
Horizontal well model
[00077] The horizontal well model also assumes discrete segments such that the formation in each segment is ID layering with dip and azimuth, and such that the trajectory azimuth is constant in each segment. However, the horizontal well model constrains the layer boundaries in each segment to be parallel to the trajectory, i.e., constant relative dip 90°. Fig. 8 shows the model for a typical horizontal well segment when viewed in the curtain section. The model parameters for each segment are the distance zjayers in TVD from the first log-point to each layer boundary, layer densities pj layers, and layer photoelectric factors PEFjayers . The number of segments, segment sizes and number of layers J^ayers in each segment are determined and fixed from the data in a pre-processing step.
[00078] The horizontal well model is applicable to intervals with variable relative dip close to 90°, i.e., with borehole parallel or near parallel to the surrounding formation layers. For all these intervals, segments must be chosen sufficiently small (not below the LS axial resolution) so that the model error in each segment - due to variable trajectory inclination or layer boundaries not parallel to the trajectory - corresponds to change in responses smaller than the data noise.
[00079] To estimate the contribution from model error, perfect synthetic data was generated using the fast-forward model and inverted using the horizontal well model with different segment sizes. For segment sizes 2 ft, the average reconstruction error in each segment was less than 0.015 g/cc for relative dip 88° to 92°.
Borehole model
[00080] The mud properties, pmuci and PEFmud, are defined in the same trajectory segments used for the formation model, and are assumed to be uniform in each segment. The geometry of the borehole is defined in terms of the radial distance from the tool axis to the borehole wall, denoted by r(9,MD), where Θ is the tool azimuth. Given r(9,MD), the corresponding stand-off s(9,MD) is defined as
Figure imgf000016_0001
where rtooi is the tool radius. For each segment, r(9,MD) may be parameterized using one of the following models.
Cylindrical Geometry
[00081] The borehole in each segment is assumed to be a cylinder parameterized by the borehole radius rBn and the stand-off Shot at 9=180°, as shown in Fig. 9. The radial distance of the borehole wall from the tool axis is τ(β, MD) = d2sin20 + TgH + dcosO, and d≡ rBH— ( t00l + sbot). Staircase Geometry
[00082] The radial distance to the borehole wall is assumed to be piecewise-constant function of both Θ and MD, defined by a discrete radius for each sector and each log-point.
Harmonic Geometry
[00083] The harmonic geometry can describe an arbitrary borehole surface. The radial distance to the borehole wall is a harmonic function of Θ:
MD) = a0 (MD)
Figure imgf000017_0001
[00084] Here, the harmonic coefficients a0, a„ and b„ for n=\ ,2,...N, are piecewise-constant functions defined at each log point. For the data-sets considered in this paper, N=2 was found to be sufficient.
INVERSION-BASED WORKFLOWS FOR PROCESSING NUCLEAR DENSITY IMAGES IN HIGH-ANGLE
WELLS
[00085] The steps for processing a high-angle interval are shown in Fig. 10. Pre-processing workflow for high-angle wells
[00086] The high-angle pre-processing workflow consists of the following steps:
[00087] Log-Squaring - A log-squaring algorithm is applied to the bottom quadrant density (ROBB) log to construct a layered density profile along MD. A photoelectric factor for each layer is obtained by sampling the UAPP bottom quadrant log at the layer mid-points.
[00088] Sinusoid Extraction - The amplitude and phase of sinusoidal features corresponding to layer boundaries are extracted from the COMP image either manually, or by an automatic method. For each layer boundary, the extracted phase and amplitude are converted to relative dip and relative azimuth.
[00089] Trajectory Segmentation - An adaptive algorithm is used to discretize the trajectory such that the COMP image for each segment includes at least half a sinusoidal feature, which improves sensitivity to dip. The algorithm also constrains segment end-points to coincide with layer mid-points, so that each segment includes at least half of a layer, improving sensitivity to properties of truncated layers.
[00090] Definition of Initial Model - The output from log-squaring and sinusoid extraction is used to define the number of layers, their properties and boundaries, apparent formation dip, and relative azimuth for each segment. The layer boundaries from log-squaring must be shifted in MD in order to obtain their locations referenced to the tool axis. The initial borehole geometry is defined as a cylinder with diameter equal to the bit-size. Alternatively, the initial borehole geometry may also be defined using the density caliper. The initial mud density and mud photoelectric factor are specified by the user.
Inversion workflow for high-angle wells
[00091] As shown in Fig. 10, in each step of the inversion workflow, a subset of model parameters is optimized from a subset of channels which are selected based on their sensitivity to the parameters being optimized, resulting in a more robust workflow.
[00092] Gauss-Newton optimization with line search, adaptive regularization and parameter constraints is used to minimize the cost function, defined to be the weighted Z2-norm error between the measurements and the modeled responses over all log-points in a processing window (see below for definition of processing window). The weights in the cost function are the reciprocal of the measurement noise level for the channel used (0.015 g/cc for SS, LS, COMP, and 0.5 b/e- for UAPP). Measurements from all sectors are included with equal weights; no special weighting is applied to sectors from the bottom of the hole.
[00093] Inversion workflow steps are as follows:
1. Get initial model parameters for a segment from the output of pre-processing.
2. Assume cylindrical borehole geometry, and invert for the borehole radius and mud properties using the SS, DRHO, and UAPP channels for all sectors and log-points in the processing window.
3. Invert for the geometry of the borehole using the staircase or harmonic parameterization using the SS and DRHO channels for that log-point or sector. 4. Keep mud properties and borehole geometry fixed from previous steps. Invert for formation parameters using LS, COMP, and UAPP channels for all sectors and log-points in the processing window.
5. Repeat steps 2-4 until all parameters are converged to a specified tolerance, or until the specified maximum number of repetitions is reached.
Sliding window scheme
[00094] To enforce consistency of models between adjacent segments, for each segment, the model parameters are inverted from measurements in a processing window that includes all log- points in the segment, as well as log-points from buffer regions on each side of the segment. Formation layers from the previous and subsequent segments that influence measurements in the buffer regions are added to the formation model in each segment. The properties and thicknesses of the layers included from the previous segment are held fixed during the inversion.
Post-processing workflow for high-angle wells
[00095] The purpose of the post-processing workflow is to visualize the inversion results, and to compute various measures of the quality of the results. The formation model is visualized as a sequence of ID multi-layer dipping formations in the curtain-section plane or as the inversion- based borehole-corrected image, which is obtained by projecting the layer properties in each segment onto a cylinder with radius rtoo/+EPL, where EPL is the effective penetration length of the COMP channel. For quality control we can also visualize the inversion outputs, individual and average misfit for all measurements and sectors. If required, the uncertainties in the final model parameters for each segment may be estimated from the model covariance matrix.
INVERSION-BASED WORKFLOWS FOR PROCESSING NUCLEAR DENSITY IMAGES IN HORIZONTAL WELLS
[00096] The steps for processing a horizontal interval are shown in Fig. 11 and are grouped in a similar manner to the high-angle workflow, with a pre-processing workflow, an inversion workflow, and a post-processing workflow.
Pre-processing workflow for horizontal wells
Trajectory Segmentation [00097] The segment sizes must be chosen so as to balance the opposing requirements of a small window size for minimizing model error, and a large window size for improving inversion robustness. Based on tests with perfect synthetic data described earlier, we use segment size of 2 ft in MD.
Definition of initial model
[00098] The assumption that the relative dip in each segment is close to 90° allows the approximate layer locations and densities to be derived from the width and density of each stripe in the COMP image by squaring the azimuthal response using the following algorithm:
1. Compute ρ(θ) by averaging the COMP image over the MD interval spanning the segment.
2. Extract phase by computing centroid of ρ(θ). Shift ρ(θ) to remove phase and symmetrize shifted response about 9=0°.
3. Square symmetrized response in the interval from θ=[0°,180°], to obtain number of layers in segment, boundary locations zjayer and layer densities pj layer.
4. Use the UAPP image to get layer photoelectric factors PEFjayer.
5. Map layer boundaries from Θ to TVD: zjayer =-(rtooi+EP ) cos$Jayerl cos a,.
The initial parameters for the borehole model are defined using the same procedure described for the high-angle pre-processing workflow.
Inversion workflow for horizontal wells
[00099] The steps in the horizontal well inversion workflow are similar to steps of the high- angle inversion workflow:
1. Obtain initial model parameters for a segment from the pre-processing output.
2. Invert for mud parameters and radius of cylindrical borehole using SS, UAPP, and DRHO channels.
3. Invert for borehole geometry parameters using SS and DRHO channels.
4. Invert for formation parameters using LS, UAPP, and COMP channels.
5. Repeat steps 2-4 until parameters converge, or until specified number of repetitions is reached. [000100] Here also, the density images are processed segment by segment using a sliding window scheme analogous to the one used for the high-angle well workflow described earlier.
Post-processing workflow for horizontal wells
[000101] The details of the horizontal well post-processing are similar to the high-angle well post-processing, with visualization of the formation model in the curtain-section plane, borehole- corrected image, and reconstructed images for all channels.
VALID A TION ON SYNTHETIC DA TA EXAMPLES
High-angle well example: Noisy data-set with variable rel. dip from 82° to 87°
[000102] We tested the workflow on a formation model with dip βΩ=5° and zero azimuth (ar=0°), with multiple layers of thicknesses 0.5 ft, 1 ft and 2 ft, and layer densities alternating between 1.95 g/cc and 2.6 g/cc, and photoelectric factor alternating between 1.49 b/e- and 2.9 b/e-. The borehole is a cylindrical,
Figure imgf000021_0001
inch), and the drilling mud is fresh water.
[000103] The 8.25 inch tool size was modeled over the 200 ft trajectory section with inclination varying linearly from 77° to 82°. Gaussian noise with zero mean and standard deviation of 0.02 g/cc, which is comparable to the amplitude of real measurement noise, was added to the synthetic measurements for all channels. The initial mud properties were set to their true values.
[000104] The output from the pre-processing workflow is shown in Fig. 12. Log-squaring and sinusoid extraction are used to build an initial segmented multi-layer formation model. The initial densities in the thin layers have errors due to shoulder-bed effect. The inverted, fully shoulder-bed corrected profile is shown in Fig 13.
[000105] Fig. 14 shows the measured, initial, and final COMP density images, as well as reconstruction errors for the initial guess and the solution. The initial error shows large systematic discrepancies due to errors in the initial dip, initial layer densities and boundary locations obtained from log-squaring and sinusoid extraction. These discrepancies are corrected by the inversion. Figs. 15 and 16 show the measured, initial, and final images for the LS and SS densities, and corresponding error images. The initial error shows similar systematic
discrepancies as the initial error for the COMP channel. In addition, due to the sensitivity of the LS and SS channels to borehole effects, there are large systematic discrepancies near the top of the hole due to the error in the initial borehole geometry. These discrepancies are absent from the final error due to the reconstruction of the borehole geometry.
[000106] The standoff from the inversion is shown in Fig. 17. The final errors appear random for all channels and the mean and standard deviation of the error for each sector is shown in Fig. 18. The errors are comparable to the measurement noise, which confirms the accuracy of the inversion for this data-set. Fig. 19 shows the borehole-corrected image and final formation model in the curtain-section view.
Horizontal well example 1: Noisy data-set with constant rel. dip 89°
[000107] We tested the horizontal well workflow on a formation model with multiple thin layers of thickness 0. 75 ft and 0.5 ft, and layer densities alternating between 1.95 g/cc and 2.6 g/cc, photoelectric factor alternating between 1.49 b/e- and 2.9 b/e-, and zero dip (βΩ=αΓ=0°). The borehole is cylindrical, rBH=4.5 inch, the bottom stand-off is zero, and the drilling mud is fresh water. The 8.25 inch tool size was modeled over the 200 ft of trajectory at constant inclination of 89°. Gaussian noise with standard deviation equal to the accuracy of the fast- forward model was added to the measurements.
[000108] The reconstruction of COMP, LS and SS density images by the inversion is shown in Figs. 18-20. The results in curtain section view are shown in Fig. 21, along with the final borehole and shoulder-bed corrected density image. We observe consistent reconstruction of boundaries. Boundaries are detected before crossing, and as expected they are more accurate for lower beds than for upper beds, where there is a possibility of ambiguity between the boundary and large standoff in the 9.0 inch diameter borehole.
Horizontal well example 2: Noisy data-set with sinusoidal rel. dip 88°-92°
[000109] We tested the horizontal well workflow on a 320 ft long section of trajectory with variable inclination from 88° to 92°, with the build/drop rate 1° per 20ft, in a test formation with zero dip (βΩ=αΓ=0°) traversing a sequence of layers 0.5, 0.75, 1 and 2 ft thick, with density alternating between 1.95 g/cc and 2.6 g/cc, and photoelectric factor alternating between 1.49 b/e- and 2.9 b/e-. The borehole is a cylinder with
Figure imgf000022_0001
inch, the tool is 8.25 inch size, with zero bottom stand-off
Figure imgf000022_0002
inch), and the drilling mud is fresh water. Gaussian noise with zero mean and standard deviation equal to the accuracy of the fast-forward model was added to the measurements for each channel. The results of the inversion are shown in Figs. 22-23.
FIELD EXAMPLES
High-angle well example: North Sea data set
[000110] We apply the high-angle workflow to a 180 ft interval from a data-set acquired in a North Sea field, for which the COMP image is characterized by a sequence of down-section sinusoidal features. Fig. 24 shows the output from pre-processing, as a result of which the trajectory has been discretized into segments with a ID multi-layer formation in each segment. The initial borehole was assumed to be a cylinder with diameter equal to the bit- size of 8.5 inch, and the initial mud density and mud PEF were set to 1.4 g/cc and 10 b/e-, respectively.
[000111] Fig. 25 shows the reconstructed COMP bottom quadrant log, and final layer profile, dip, and azimuth. Note the shoulder-bed corrections in thin layers and corrections to the initial dip of approximately 1°.
[000112] Figs. 26-28 show the reconstructed images for the COMP, LS and SS density channels, as well as the corresponding errors. COMP image reconstruction shows that errors are reduced in all quadrants due to optimization of the layer properties, thicknesses, and dip. Final errors are less than 5 %, and are larger in the upper quadrant relative to the bottom quadrant, possibly due to inaccuracies in the mud properties or lower accuracy of fast- forward model for the large standoff. LS and SS reconstruction is similar, but in addition, significant reduction in the initial error in the upper quadrant is seen due to optimization of the borehole geometry.
[000113] Fig. 29 shows the average and standard deviation of the reconstruction error for each sector. The final average error is less than 0.03 g/cc for all channels over all sectors, with close to zero average error in the bottom quadrant. The final standard deviation of the error is less than 0.03 g/cc for all channels in all sectors.
[000114] Fig. 30 shows the final borehole geometry, which was parameterized using a harmonic model with two harmonics in the expansion. The same figure also shows the geometry from the density caliper. The density caliper has sinusoidal artifacts whose locations correlate with layers in the COMP image. These artifacts are not present in the geometry from the inversion. Fig. 31 shows the final borehole-corrected image and final formation model in the curtain-section view. Field data-set from a gas shale well
[000115] We apply the high-angle workflow to a 100 ft down-section interval of the data-set acquired in the Marcellus gas shale well in NE Pennsylvania. Fig. 32 shows the output from preprocessing, where the trajectory has been discretized into segments with associated initial ID multi-layer formation in each segment. The initial borehole geometry was assumed to be a cylinder with diameter equal to the bit-size (8.5 inch), and the mud density and mud PEF were set to 1.4 g/cc and 50 b/e-, respectively.
[000116] Fig. 33 shows the reconstructed layered density profile, dip, and azimuth. The reconstructed dip is corrected up to 2° compared to initial guess.
[000117] Fig. 34 shows the measured, initial, and final images for the COMP channel, as well as the initial and final errors. Fig. 35-36 show the measured, initial, and final images for the LS and SS channels, as well as the initial and final errors. Fig. 37 shows the average and standard deviation of the initial and final errors for each sector. The inversion significantly reduces the mismatch between modeling and the measurements. Fig. 38 shows the final borehole geometry. Fig. 39 shows the final borehole-corrected image and final formation model in the curtain- section plane.
Horizontal Well Example: North Sea data set
[000118] We apply the horizontal well workflow to 350 ft long data interval from a North Sea field. The initial mud density and mud photoelectric factor were set to 1.4 g/cc and 10 b/e-, respectively.
[000119] The measured COMP image, the initial guess, inversion reconstructed image and corresponding misfit are shown in Fig. 40. Fig. 41 shows the average and standard deviation of the initial and final errors for each sector. The inversion significantly reduces the mismatch between modeling and the measurements. Borehole geometry in terms of standoff is shown in Fig 42. Initial geometry was cylinder with 8.5 inch diameter. Harmonic parameterization with two harmonics in the expansion was used. The inversion-based caliper does not have sinusoidal artifacts at MD= 100- 180ft.
[000120] Figure 43 bottom track shows the final formation model in the curtain-section plane. The green curve is the tool trajectory and the dashed green curves are the borehole surfaces.

Claims

Claims
1. A method for characterizing a subterranean formation traversed by a wellbore, comprising: collecting data from the formation using a density tool wherein the tool collects data to form an azimuthal image; characterizing a section of the formation comprising data and images acquired in a high angle wellbore section or horizontal wellbore section using a parametric model; and performing an inversion using apparent densities and volumetric photoelectric factor images to build a formation model wherein the inversion is tailored for high angle wellbore sections and/or horizontal wellbore sections.
2. The method of claim 1, further comprising post-processing comprising visualizing the inversion built model.
3. The method of claim 1, wherein the formation model further comprises borehole geometry and layer photoelectric models.
4. The method of claim 1, further comprising pre-processing the data.
5. The method of claim 4, wherein pre-processing is selected from the group consisting of log- squaring, sinusoid extraction from an image, trajectory segmentation, and a combination thereof.
6. The method of claim 5, wherein the log is the azimuthal compensated density image acquired at a single measurement depth.
7. The method of claim 5, wherein the log is the bottom quadrant density from the compensated image.
8. The method of claim 5, wherein the sinusoid extraction comprised fitting the contour using sine function and finding amplitude and phase.
9. The method of claim 5, wherein the sinusoid extraction comprises processing to calculate dip and azimuth of a layer boundary from sinusoid phase and amplitude.
10. The method of claim 5, wherein the sinusoid extraction comprises a marching-squares contour extraction and log-squaring.
11. The method of claim 5, trajectory segmentation comprises computation of size of inversion processing window using computed relative dip.
12. The method of claim 11, wherein a window size is inversely proportional to relative dip.
13. The method of claim 5, where trajectory segmentation comprises selection between high angle and horizontal inversion workflows based on relative dip.
14. The method of claim 13, wherein the relative dip is about 87° to about 93°.
15. The method of claim 1, wherein the performing an inversion comprises inverting for borehole geometry and/or mud density and/or photoelectric factor using a short spacing density and/or differential density and/or photeoelectric factor measurements.
16. The method of claim 15, wherein inverting the borehole geometry is based on parametric model of the borehole geometry and functional form of borehole radius or standoff between the sensor and the formation.
17. The method of claim 15, wherein inverted measurements are an azimuthal sector of an image or full image.
18. The method of claim 15, wherein parameterization is using harmonic representation of borehole radius.
19. The method of claim 15, wherein standoff is inverted for each sector independently.
20. The method of claim 1, wherein the performing an inversion comprises inverting formation parameters using long spacing density, compensated density, and/or apparent volumetric photoelectric factor azimuthal images.
21. The method of claim 1, wherein the performing an inversion is selected from the group consisting of inverting formation layer thickness, densities, photoelectric factors, layer dip, or a combination thereof.
22. The method of claim 1, wherein the inversion is a Gauss-Newton inversion method.
23. The method of claim 1, wherein the tool is acquiring data while azimuthally rotating.
24. The method of claim 1, wherein the tool is logging the formation while drilling.
PCT/US2012/061585 2011-10-31 2012-10-24 Inversion-based workflow for processing nuclear density images in high-angle and horizontal wells WO2013066682A1 (en)

Priority Applications (4)

Application Number Priority Date Filing Date Title
GB1406621.1A GB2509450A (en) 2011-10-31 2012-10-24 Inversion-based workflow for processing nuclear density images in high-angle and horizontal wells
US14/352,566 US9869792B2 (en) 2011-10-31 2012-10-24 Inversion-based workflow for processing nuclear density images in high-angle and horizontal wells
NO20140524A NO20140524A1 (en) 2011-10-31 2014-04-23 INVERSION-BASED WORKFLOW TO TREAT CORE DENSITY IMAGES IN STEEL-ANGLED AND HORIZONTAL WELLS
US15/845,515 US10809416B2 (en) 2011-10-31 2017-12-18 Inversion-based workflow for processing nuclear density images in high-angle and horizontal wells

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
US201161553810P 2011-10-31 2011-10-31
US61/553,810 2011-10-31

Related Child Applications (2)

Application Number Title Priority Date Filing Date
US14/352,566 A-371-Of-International US9869792B2 (en) 2011-10-31 2012-10-24 Inversion-based workflow for processing nuclear density images in high-angle and horizontal wells
US15/845,515 Continuation US10809416B2 (en) 2011-10-31 2017-12-18 Inversion-based workflow for processing nuclear density images in high-angle and horizontal wells

Publications (1)

Publication Number Publication Date
WO2013066682A1 true WO2013066682A1 (en) 2013-05-10

Family

ID=48192624

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/US2012/061585 WO2013066682A1 (en) 2011-10-31 2012-10-24 Inversion-based workflow for processing nuclear density images in high-angle and horizontal wells

Country Status (4)

Country Link
US (2) US9869792B2 (en)
GB (1) GB2509450A (en)
NO (1) NO20140524A1 (en)
WO (1) WO2013066682A1 (en)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2015051350A1 (en) * 2013-10-04 2015-04-09 Schlumberger Canada Limited Inversion-based workflow for consistent interpretation of nuclear density images in horizontal wells
WO2015131016A1 (en) * 2014-02-28 2015-09-03 Schlumberger Canada Limited Automatic method for three-dimensional structural interpretation of borehole images acquired in high-angle and horizontal wells
US9366135B2 (en) 2013-10-08 2016-06-14 Exxonmobil Upstream Research Company Automatic dip picking from wellbore azimuthal image logs
US10451769B2 (en) 2015-01-26 2019-10-22 Schlumberger Technology Corporation Method for determining petrophysical properties from logging measurements
US20230102131A1 (en) * 2018-05-11 2023-03-30 Schlumberger Technology Corporation Geologic formation neutron porosity system

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB2509450A (en) 2011-10-31 2014-07-02 Schlumberger Holdings Inversion-based workflow for processing nuclear density images in high-angle and horizontal wells
GB2511744B (en) * 2013-03-11 2020-05-20 Reeves Wireline Tech Ltd Methods of and apparatuses for identifying geological characteristics in boreholes
US9811882B2 (en) * 2014-09-30 2017-11-07 Electronics And Telecommunications Research Institute Method and apparatus for processing super resolution image using adaptive preprocessing filtering and/or postprocessing filtering
US10222498B2 (en) 2015-05-15 2019-03-05 Weatherford Technology Holdings, Llc System and method for joint inversion of bed boundaries and petrophysical properties from borehole logs
CN110443283A (en) * 2019-07-11 2019-11-12 长江大学 A kind of discriminant classification method and system of tight gas reservoir horizontal well

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20090093962A1 (en) * 2003-10-03 2009-04-09 Ridvan Akkurt System and methods for t1-based logging
US20100004867A1 (en) * 2008-07-01 2010-01-07 Schlumberger Technology Corporation Forward models for gamma ray measurement analysis of subterranean formations
US20100312478A1 (en) * 2009-06-04 2010-12-09 Tabanou Jacques R Log processing in highly deviated wellbores
US20110077867A1 (en) * 2008-03-11 2011-03-31 Michael Evans Method and apparatus for determining formation and fluid properties
US20110168879A1 (en) * 2007-02-05 2011-07-14 Schlumberger Technology Corporation Nuclear Tool

Family Cites Families (17)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE4104358A1 (en) 1991-02-13 1992-08-20 Implex Gmbh IMPLANTABLE HOER DEVICE FOR EXCITING THE INNER EAR
US5831935A (en) * 1996-03-05 1998-11-03 Chevron U.S.A. Inc. Method for geophysical processing and interpretation using seismic trace difference for analysis and display
US5960371A (en) * 1997-09-04 1999-09-28 Schlumberger Technology Corporation Method of determining dips and azimuths of fractures from borehole images
US6308136B1 (en) * 2000-03-03 2001-10-23 Baker Hughes Incorporated Method of interpreting induction logs in horizontal wells
WO2004076816A1 (en) * 2003-02-27 2004-09-10 Schlumberger Surenco Sa Estimation of formation characteristics in wells
US7814036B2 (en) * 2003-06-19 2010-10-12 Haliburton Energy Services, Inc. Processing well logging data with neural network
US7279677B2 (en) 2005-08-22 2007-10-09 Schlumberger Technology Corporation Measuring wellbore diameter with an LWD instrument using compton and photoelectric effects
GB2453478B (en) * 2006-07-25 2011-08-17 Exxonmobil Upstream Res Co Method for determining physical properties of structures
EP2238477A4 (en) * 2007-12-19 2016-08-24 Exxonmobil Upstream Res Co Gamma ray tool response modeling
US9557440B2 (en) * 2009-07-17 2017-01-31 Baker Hughes Incorporated Radial waves in a borehole and stoneley waves for measuring formation permeability and electroacoustic constant
US8433518B2 (en) 2009-10-05 2013-04-30 Schlumberger Technology Corporation Multilevel workflow method to extract resistivity anisotropy data from 3D induction measurements
CA2787149A1 (en) * 2010-01-22 2011-07-28 Schlumberger Canada Limited Real-time formation anisotropy and dip evaluation using tri-axial induction measurements
AU2010357606B2 (en) * 2010-07-16 2014-03-13 Halliburton Energy Services, Inc. Efficient inversion systems and methods for directionally-sensitive resistivity logging tools
US8754651B2 (en) 2010-11-15 2014-06-17 Schlumberger Technology Corporation System and method for imaging properties of subterranean formations
US9557433B2 (en) * 2011-03-23 2017-01-31 Seismic Global Ambient, Llc Fracture imaging methods employing skeletonization of seismic emission tomography data
US9239403B2 (en) * 2011-08-29 2016-01-19 Hallibburton Energy Services, Inc. Apparatus and methods of controlling recordation of resistivity-related readings in determining formation resistivity
GB2509450A (en) 2011-10-31 2014-07-02 Schlumberger Holdings Inversion-based workflow for processing nuclear density images in high-angle and horizontal wells

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20090093962A1 (en) * 2003-10-03 2009-04-09 Ridvan Akkurt System and methods for t1-based logging
US20110168879A1 (en) * 2007-02-05 2011-07-14 Schlumberger Technology Corporation Nuclear Tool
US20110077867A1 (en) * 2008-03-11 2011-03-31 Michael Evans Method and apparatus for determining formation and fluid properties
US20100004867A1 (en) * 2008-07-01 2010-01-07 Schlumberger Technology Corporation Forward models for gamma ray measurement analysis of subterranean formations
US20100312478A1 (en) * 2009-06-04 2010-12-09 Tabanou Jacques R Log processing in highly deviated wellbores

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2015051350A1 (en) * 2013-10-04 2015-04-09 Schlumberger Canada Limited Inversion-based workflow for consistent interpretation of nuclear density images in horizontal wells
US10914861B2 (en) 2013-10-04 2021-02-09 Schlumberger Technology Corporation Inversion-based workflow for consistent interpretation of nuclear density images in horizontal wells
US9366135B2 (en) 2013-10-08 2016-06-14 Exxonmobil Upstream Research Company Automatic dip picking from wellbore azimuthal image logs
WO2015131016A1 (en) * 2014-02-28 2015-09-03 Schlumberger Canada Limited Automatic method for three-dimensional structural interpretation of borehole images acquired in high-angle and horizontal wells
EP3111041A4 (en) * 2014-02-28 2017-12-06 Services Pétroliers Schlumberger Automatic method for three-dimensional structural interpretation of borehole images acquired in high-angle and horizontal wells
US10466375B2 (en) 2014-02-28 2019-11-05 Schlumberger Technology Corporation Automatic method for three-dimensional structural interpretation of borehole images acquired in high-angle and horizontal wells
US10451769B2 (en) 2015-01-26 2019-10-22 Schlumberger Technology Corporation Method for determining petrophysical properties from logging measurements
US20230102131A1 (en) * 2018-05-11 2023-03-30 Schlumberger Technology Corporation Geologic formation neutron porosity system
US11774631B2 (en) * 2018-05-11 2023-10-03 Schlumberger Technology Corporaton Geologic formation neutron porosity system

Also Published As

Publication number Publication date
US20180136362A1 (en) 2018-05-17
US20140286539A1 (en) 2014-09-25
GB2509450A (en) 2014-07-02
US10809416B2 (en) 2020-10-20
NO20140524A1 (en) 2014-05-06
GB201406621D0 (en) 2014-05-28
US9869792B2 (en) 2018-01-16

Similar Documents

Publication Publication Date Title
US10809416B2 (en) Inversion-based workflow for processing nuclear density images in high-angle and horizontal wells
EP3111041B1 (en) Automatic method for three-dimensional structural interpretation of borehole images acquired in high-angle and horizontal wells
US9482775B2 (en) Real-time formation anisotropy and dip evaluation using tri-axial induction measurements
Grana et al. Quantitative log interpretation and uncertainty propagation of petrophysical properties and facies classification from rock-physics modeling and formation evaluation analysis
US10451765B2 (en) Post-well reservoir characterization using image-constrained inversion
US10914861B2 (en) Inversion-based workflow for consistent interpretation of nuclear density images in horizontal wells
US8527204B2 (en) Volume of investigation based density image processing
US20140257703A1 (en) Real-Time Formation Anisotropy And Dip Evaluation Using Multiaxial Induction Measurements
CN106154322A (en) Log curve correction method and apparatus
US20200033501A1 (en) Improved methods relating to quality control
Wang et al. Seismic attributes for characterization and prediction of carbonate faulted karst reservoirs in the Tarim Basin, China
Shetty et al. Inversion-based workflows for interpretation of nuclear density images in high-angle and horizontal wells
Shetty et al. 3D Parametric Inversion for Interpretation of Logging-While-Drilling Density Images in High-Angle and Horizontal Wells1
CN114542055A (en) Method and device for focusing and processing array induction logging data of inclined well or inclined stratum
George et al. Challenges and key learning for developing tight carbonate reservoirs
BR102020012943A2 (en) COMPUTATIONAL SCRIPT FOR IMAGE TREATMENT AND ITS APPLICATION IN EASY IMAGE DETERMINATION METHOD
CN110850504B (en) Shale density parameter prestack inversion method based on uranium curve quasi-impedance constraint
Maharaja Global net-to-gross uncertainty Assessment at reservoir appraisal stage
Boesing et al. High-Angle Formation Response Modeling for More Accurate Petrophysics–Challenges and Experiences from Multiple Cases
Bolt et al. A Methodology for Portraying Three-Dimensional Positional Uncertainty Using Along-Hole Depth, Inclination and Azimuth Measurement Accuracies
Zhou et al. Quantitative formation evaluation in high angle and horizontal wells-A step change in reservoir characterization
EP2113795A1 (en) Volume Of Investigation Based Image Processing
Stromberg et al. Reservoir Quality, Net-to-Gross, and Fluid Identification in Laminated Reservoirs from a new generation of NMR logging tools. Examples from the Gharif Formation, Southern Oman.
Zhou Newly developed technology for array resistivity logging
Fishburn et al. SPWLA 39th Annual Logging Symposium, May 26-29, 1998 PRACTICAL INVERSION OF HIGH-DEFINITION INDUCTION LOGS USING A PRIORI INFORMATION

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 12845442

Country of ref document: EP

Kind code of ref document: A1

ENP Entry into the national phase

Ref document number: 1406621

Country of ref document: GB

Kind code of ref document: A

Free format text: PCT FILING DATE = 20121024

WWE Wipo information: entry into national phase

Ref document number: 1406621.1

Country of ref document: GB

WWE Wipo information: entry into national phase

Ref document number: 14352566

Country of ref document: US

NENP Non-entry into the national phase

Ref country code: DE

122 Ep: pct application non-entry in european phase

Ref document number: 12845442

Country of ref document: EP

Kind code of ref document: A1