- Home
- A-Z Publications
- Geophysical Prospecting
- Previous Issues
- Volume 53, Issue 1, 2005
Geophysical Prospecting - Volume 53, Issue 1, 2005
Volume 53, Issue 1, 2005
-
-
Residual dip moveout in VTI media
Authors Tariq Alkhalifah and Maarten V. De HoopABSTRACTDip‐moveout (DMO) correction is often applied to common‐offset sections of seismic data using a homogeneous isotropic medium assumption, which results in a fast execution. Velocity‐residual DMO is developed to correct for the medium‐treatment limitation of the fast DMO. For reasonable‐sized velocity perturbations, the residual DMO operator is small, and thus is an efficient means of applying a conventional Kirchhoff approach. However, the shape of the residual DMO operator is complicated and may form caustics. We use the Fourier domain for the operator development part of the residual DMO, while performing the convolution with common‐offset data in the space–time domain. Since the application is based on an integral (Kirchhoff) method, this residual DMO preserves all the flexibility features of an integral DMO. An application to synthetic and real data demonstrates effectiveness of the velocity‐residual DMO in data processing and velocity analysis.
-
-
-
Electromagnetic fields in a steel‐cased borehole
Authors Ki Ha Lee, Hee Joon Kim and Toshihiro UchidaABSTRACTThe development of an electromagnetic numerical modelling scheme for a magnetic dipole in an arbitrary casing segment in an inhomogeneous conductivity background has been difficult, due to the very high electrical conductivity and magnetic permeability contrasts between the steel casing and the background medium. To investigate the effect of steel casing efficiently, we have developed an accurate but simple finite‐element modelling scheme to simulate electromagnetic fields in a medium of cylindrically symmetric conductivity structures. In order to preserve the cylindrical symmetry in the resulting electromagnetic fields, a horizontal loop current source is used throughout. One of the main advantages of the approach is that the problem is scalar when formulated using the azimuthal electric field, even if the casing is both electrically conductive and magnetically permeable. Field calculations have been made inside the cased borehole as well as in another borehole which is not cased. Careful analyses of the numerical modelling results indicate that the anomaly observed in a cross‐borehole configuration is sensitive enough to be used for tomographic imaging.
-
-
-
Traveltime computation by wavefront‐orientated ray tracing
Authors Radu Coman and Dirk GajewskiABSTRACTFor multivalued traveltime computation on dense grids, we propose a wavefront‐orientated ray‐tracing (WRT) technique. At the source, we start with a few rays which are propagated stepwise through a smooth two‐dimensional (2D) velocity model. The ray field is examined at wavefronts and a new ray might be inserted between two adjacent rays if one of the following criteria is satisfied: (1) the distance between the two rays is larger than a predefined threshold; (2) the difference in wavefront curvature between the rays is larger than a predefined threshold; (3) the adjacent rays intersect. The last two criteria may lead to oversampling by rays in caustic regions. To avoid this oversampling, we do not insert a ray if the distance between adjacent rays is smaller than a predefined threshold. We insert the new ray by tracing it from the source. This approach leads to an improved accuracy compared with the insertion of a new ray by interpolation, which is the method usually applied in wavefront construction. The traveltimes computed along the rays are used for the estimation of traveltimes on a rectangular grid. This estimation is carried out within a region bounded by adjacent wavefronts and rays. As for the insertion criterion, we consider the wavefront curvature and extrapolate the traveltimes, up to the second order, from the intersection points between rays and wavefronts to a gridpoint. The extrapolated values are weighted with respect to the distances to wavefronts and rays. Because dynamic ray tracing is not applied, we approximate the wavefront curvature at a given point using the slowness vector at this point and an adjacent point on the same wavefront. The efficiency of the WRT technique is strongly dependent on the input parameters which control the wavefront and ray densities. On the basis of traveltimes computed in a smoothed Marmousi model, we analyse these dependences and suggest some rules for a correct choice of input parameters. With suitable input parameters, the WRT technique allows an accurate traveltime computation using a small number of rays and wavefronts.
-
-
-
Sub‐basalt decoupled walkaway VSP imaging
More LessABSTRACTWalkaway vertical seismic profile (VSP) data acquired in basalt‐covered areas can be used to improve knowledge of the sub‐basalt structure. A synthetic example and a case study from the North Atlantic (UK) show that elastic two‐way downward‐continuation migration combined with the stationary‐phase principle is well suited to the processing of VSP data. Vector data are processed using decoupled elastic migration algorithms in both isotropic and anisotropic media. To illustrate the value of decoupled imaging equations, conventional PP imaging is carried out on the enhanced VSP data and compared with the decoupled scheme. Decoupled vector migration operates directly on the displacement vector, and uses various wave modes. Downgoing waves are migrated to image basalt lava flows and measure their anisotropy. Upgoing waves are used for high‐resolution sub‐basalt imaging.
-
-
-
Simulation of full‐waveform log in saturated cracked formations using Hudson's approach
Authors E. Kazatchenko, M. Markov and A. MousatovABSTRACTWe study the propagation of elastic waves that are generated in a fluid‐filled borehole surrounded by a cracked transversely isotropic medium. In the model studied the anisotropy and borehole axes coincide. To obtain the effective elastic moduli of a cracked medium we have applied Hudson's theory that enables the determination of the overall properties as a function of the crack orientation in relation to the symmetry axis of the anisotropic medium. This theory takes into account the hydrodynamic mechanism of the elastic‐wave attenuation caused by fluid filtration from the cracks into a porous matrix.
We have simulated the full waveforms generated by an impulse source of finite length placed on the borehole axis. The kinematic and dynamic parameters of the compressional, shear and Stoneley waves as functions of the matrix permeability, crack orientation and porosity were studied. The modelling results demonstrated the influence of the crack‐system parameters (orientation and porosity) on the velocities and amplitudes of all wave types. The horizontally orientated cracks result in maximal decrease of the elastic‐wave parameters (velocities and amplitudes).
Based on the fact that the shear‐ and Stoneley‐wave velocities in a transversely isotropic medium are determined by different shear moduli, we demonstrate the feasibility of the acoustic log to identify formations with close to horizontal crack orientations.
-
-
-
Interpretation of magnetic data using analytic signal derivatives
By Ahmed SalemABSTRACTThis paper develops an automatic method for interpretation of magnetic data using derivatives of the analytic signal. A linear equation is derived to provide source location parameters of a 2D magnetic body without a priori information about the nature of the source. Then using the source location parameters, the nature of the source can be ascertained. The method has been tested using theoretical simulations with random noise for two 2D magnetic models placed at different depths with respect to the observation height. In both cases, the method gave a good estimate for the location and shape of the sources. Good results were obtained on two field data sets.
-
-
-
Depth and structural index from normalized local wavenumber of 2D magnetic anomalies
Authors Ahmed Salem and Richard SmithABSTRACTRecent improvements in the local wavenumber approach have made it possible to estimate both the depth and model type of buried bodies from magnetic data. However, these improvements require calculation of third‐order derivatives of the magnetic field, which greatly enhances noise. As a result, the improvements are restricted to data of high quality. We present an alternative method to estimate both the depth and model type using the first‐order local wavenumber approach without the need for third‐order derivatives of the field. Our method is based on normalization of the first‐order local wavenumber anomalies and provides a generalized equation to estimate the depth of some 2D magnetic sources regardless of the source structure. Information about the nature of the sources is obtained after the source location has been estimated. The method was tested using synthetic magnetic anomaly data with random noise and using three field examples.
-
-
-
Airborne gravimetry using a strapped‐down LaCoste and Romberg air/sea gravity meter system: a feasibility study
Authors J. Verdun and E.E. KlingeléABSTRACTMarine gravimeters mounted on stabilized platforms are commonly used in aircraft to perform airborne gravity measurements. The role of the stabilized platform is to level the sensor mechanically, whatever the aircraft attitude. However, this compensation is generally insufficient due to the sensitivity of modern gravity sensors. Correcting the offlevel error requires that an offlevel correction calculated from positioning data be added to gravimeter measurements, which complicates not only the processing, but also the assessment of precision and resolution.
This paper is a feasibility study describing the levelling of a completely strapped‐down LaCoste and Romberg gravimeter for airborne gravimetry operation, by means of GPS positioning data. It focuses on the calculation of the sensor offlevel correction needed for the complete gravity data processing. The precision of the offlevel correction that can be achieved, in terms of GPS data precision and gravity wavelengths, is theoretically studied and estimated using the gravity and GPS data acquired during the Alpine Swiss French airborne gravity survey carried out in 1998 over the French Western Alps. While a 1 cm precision of GPS‐determined baseline coordinates is sufficient to achieve a 5 mGal precision of the offlevel correction, we maintain that this precision has to reach 1 mm to ensure a 1 mGal precision of the offlevel correction at any wavelength.
Without a stabilized platform, the onboard instrumentation becomes significantly lighter. Furthermore, the correction for the offlevel error is straightforward and calculated only from GPS data. Thus, the precision and the resolution of airborne gravity surveys should be estimated with a better accuracy.
-
-
-
Inversion of nuclear well‐logging data using neural networks
ABSTRACTThis work looks at the application of neural networks in geophysical well‐logging problems and specifically their utilization for inversion of nuclear downhole data. Simulated neutron and γ‐ray fluxes at a given detector location within a neutron logging tool were inverted to obtain formation properties such as porosity, salinity and oil/water saturation. To achieve this, the forward particle‐radiation transport problem was first solved for different energy groups (47 neutron groups and 20 γ‐ray groups) using the multigroup code EVENT. A neural network for each of the neutron and γ‐ray energy groups was trained to re‐produce the detector fluxes using the forward modelling results from 504 scenarios. The networks were subsequently tested on unseen data sets and the unseen input parameters (formation properties) were then predicted using a global search procedure. The results obtained are very encouraging with formation properties being predicted to within 10% average relative error. The examples presented show that neural networks can be applied successfully to nuclear well‐logging problems. This enables the implementation of a fast inversion procedure, yielding quick and reliable values for unknown subsurface properties such as porosity, salinity and oil saturation.
-
-
-
The effect of orthorhombic anisotropy and its implication for oil recovery and reservoir exploitation
More LessABSTRACTResults of an experimental study of shear‐ and compressional‐wave propagation in an orthorhombically anisotropic medium are presented. The experiments were performed on a physical model consisting of two sets of fractures. The first set consisted of orientated rubber inclusions simulating weak material‐filled cracks. The second set consisted of a system of closely spaced parallel fractures simulated by thin plates of epoxy resin, superimposed orthogonally on the first set. Three cases of fracture orientations within the model were identified and studied. Case 1 is analogous to a jointed fracture reservoir with one vertical set of fluid‐filled cracks or fractures and one non‐filled horizontal set. This case is referred to as JFV. Case 2 is analogous to a double fracture reservoir with one horizontal set of fluid‐filled fractures or cracks and one non‐filled vertical set. This is referred to as DFH. Case 3 is analogous to a double fracture reservoir with two vertical sets of fractures or cracks, with only one fluid‐filled. Case 3 is referred to as DFV. A pulse transmission method was performed on all three modelled cases along the three principal axes. A directional variation in the compressional‐ and shear‐wave velocities, as well as distinct shear‐wave splitting, was observed. The elastic constants for each case were determined and differences between them were noted and compared with the controlled results of both layered (transverse isotropy, TI) and vertically fractured (azimuthally anisotropic models, VF) media. The differences in elastic moduli and velocities indicate the potential of recognizing the different fracture orientations and suggest an approach to designing a method of drilling to further enhance oil recovery and reservoir exploitation.
-
-
-
Dense 3D residual moveout analysis as a tool for HTI parameter estimation
Authors Evgenii A. Kozlov and Dmitry Y. VarivodaABSTRACTThree‐dimensional residual moveout analysis is the basic step in velocity model refinement. The analysis is generally carried out using horizontal and/or vertical semblances defined on a sparse set of in‐lines or cross‐lines with densely sampled source–receiver offsets. An alternative approach, which we call dense residual moveout analysis (DRMA), is to use all the bins of a three‐dimensional survey but sparsely sampled offsets. The proposed technique is very fast and provides unbiased and statistically efficient estimates of the residual moveout. Indeed, for the sparsest possible offset distribution, when only near‐ and far‐angle stacks are used, the variance of the residual moveout estimate is only 1.4 times larger than the variance of the least‐squares estimate obtained using all offsets.
The high performance of DRMA makes it a useful tool for many applications, of which azimuthal velocity analysis is considered here. For a horizontal transverse isotropy (HTI) model, a deterministic procedure is proposed to define, at every point of residual moveout estimation, the azimuthal angle of the HTI axis of symmetry, the Thomsen anisotropy coefficients, and the interval (or root‐mean‐square) velocities in both the HTI isotropy and symmetry planes. The procedure is not restricted by DRMA assumptions; for example, it is also applicable to semblance‐based residual moveout estimates.
The high resolution of the technique is illustrated by azimuthal velocity analysis over an oilfield in West Siberia.
-
-
-
Non‐minimum‐phase wavelet estimation using second‐ and third‐order moments
By Wenkai LuABSTRACTThis paper presents a new algorithm for estimating non‐minimum‐phase seismic wavelets by using the second‐ and higher‐order statistics (HOS) of the wavelets. In contrast to many, if not most, of the HOS‐based methods, the proposed method does not need to assume that subsurface seismic reflectivity is a non‐Gaussian, statistically independent and identically distributed random process. The amplitude and phase spectra of the wavelets are estimated, respectively, using the second‐order statistics (SOS) and third‐order moment (TOM) of the wavelets, which will, in turn, be derived from the HOS of the seismic traces. In our approach, the wavelets can be ‘calculated’ from seismic traces efficiently; no optimization or inversion is necessarily required. Very good results have been obtained by applying this method to both synthetic and real‐field data sets.
-
Volumes & issues
-
Volume 72 (2023 - 2024)
-
Volume 71 (2022 - 2023)
-
Volume 70 (2021 - 2022)
-
Volume 69 (2021)
-
Volume 68 (2020)
-
Volume 67 (2019)
-
Volume 66 (2018)
-
Volume 65 (2017)
-
Volume 64 (2015 - 2016)
-
Volume 63 (2015)
-
Volume 62 (2014)
-
Volume 61 (2013)
-
Volume 60 (2012)
-
Volume 59 (2011)
-
Volume 58 (2010)
-
Volume 57 (2009)
-
Volume 56 (2008)
-
Volume 55 (2007)
-
Volume 54 (2006)
-
Volume 53 (2005)
-
Volume 52 (2004)
-
Volume 51 (2003)
-
Volume 50 (2002)
-
Volume 49 (2001)
-
Volume 48 (2000)
-
Volume 47 (1999)
-
Volume 46 (1998)
-
Volume 45 (1997)
-
Volume 44 (1996)
-
Volume 43 (1995)
-
Volume 42 (1994)
-
Volume 41 (1993)
-
Volume 40 (1992)
-
Volume 39 (1991)
-
Volume 38 (1990)
-
Volume 37 (1989)
-
Volume 36 (1988)
-
Volume 35 (1987)
-
Volume 34 (1986)
-
Volume 33 (1985)
-
Volume 32 (1984)
-
Volume 31 (1983)
-
Volume 30 (1982)
-
Volume 29 (1981)
-
Volume 28 (1980)
-
Volume 27 (1979)
-
Volume 26 (1978)
-
Volume 25 (1977)
-
Volume 24 (1976)
-
Volume 23 (1975)
-
Volume 22 (1974)
-
Volume 21 (1973)
-
Volume 20 (1972)
-
Volume 19 (1971)
-
Volume 18 (1970)
-
Volume 17 (1969)
-
Volume 16 (1968)
-
Volume 15 (1967)
-
Volume 14 (1966)
-
Volume 13 (1965)
-
Volume 12 (1964)
-
Volume 11 (1963)
-
Volume 10 (1962)
-
Volume 9 (1961)
-
Volume 8 (1960)
-
Volume 7 (1959)
-
Volume 6 (1958)
-
Volume 5 (1957)
-
Volume 4 (1956)
-
Volume 3 (1955)
-
Volume 2 (1954)
-
Volume 1 (1953)