- Home
- A-Z Publications
- Geophysical Prospecting
- Previous Issues
- Volume 64, Issue 5, 2016
Geophysical Prospecting - Volume 64, Issue 5, 2016
Volume 64, Issue 5, 2016
-
-
Sensitivity analysis and application of time‐lapse full‐waveform inversion: synthetic testing and field data example from the North Sea, Norway
Authors Hadi Balhareth and Martin LandrøABSTRACTTime‐lapse refraction can provide complementary seismic solutions for monitoring subtle subsurface changes that are challenging for conventional P‐wave reflection methods. The utilization of refraction time lapse has lagged behind in the past partly due to the lack of robust techniques that allow extracting easy‐to‐interpret reservoir information. However, with the recent emergence of the full‐waveform inversion technique as a more standard tool, we find it to be a promising platform for incorporating head waves and diving waves into the time‐lapse framework. Here we investigate the sensitivity of 2D acoustic, time‐domain, full‐waveform inversion for monitoring a shallow, weak velocity change (−30 m/s, or −1.6%). The sensitivity tests are designed to address questions related to the feasibility and accuracy of full‐waveform inversion results for monitoring the field case of an underground gas blowout that occurred in the North Sea. The blowout caused the gas to migrate both vertically and horizontally into several shallow sand layers. Some of the shallow gas anomalies were not clearly detected by conventional 4D reflection methods (i.e., time shifts and amplitude difference) due to low 4D signal‐to‐noise ratio and weak velocity change. On the other hand, full‐waveform inversion sensitivity analysis showed that it is possible to detect the weak velocity change with the non‐optimal seismic input. Detectability was qualitative with variable degrees of accuracy depending on different inversion parameters. We inverted, the real 2D seismic data from the North Sea with a greater emphasis on refracted and diving waves’ energy (i.e., most of the reflected energy was removed for the shallow zone of interest after removing traces with offset less than 300 m). The full‐waveform inversion results provided more superior detectability compared with the conventional 4D stacked reflection difference method for a weak shallow gas anomaly (320 m deep).
-
-
-
Offset‐variable density improves acoustic full‐waveform inversion: a shallow marine case study
Authors Akela Silverton, Michael Warner, Joanna Morgan and Adrian UmplebyABSTRACTWe have previously applied three‐dimensional acoustic, anisotropic, full‐waveform inversion to a shallow‐water, wide‐angle, ocean‐bottom‐cable dataset to obtain a high‐resolution velocity model. This velocity model produced an improved match between synthetic and field data, better flattening of common‐image gathers, a closer fit to well logs, and an improvement in the pre‐stack depth‐migrated image. Nevertheless, close examination reveals that there is a systematic mismatch between the observed and predicted data from this full‐waveform inversion model, with the predicted data being consistently delayed in time. We demonstrate that this mismatch cannot be produced by systematic errors in the starting model, by errors in the assumed source wavelet, by incomplete convergence, or by the use of an insufficiently fine finite‐difference mesh. Throughout these tests, the mismatch is remarkably robust with the significant exception that we do not see an analogous mismatch when inverting synthetic acoustic data. We suspect therefore that the mismatch arises because of inadequacies in the physics that are used during inversion. For ocean‐bottom‐cable data in shallow water at low frequency, apparent observed arrival times, in wide‐angle turning‐ray data, result from the characteristics of the detailed interference pattern between primary refractions, surface ghosts, and a large suite of wide‐angle multiple reflected and/or multiple refracted arrivals. In these circumstances, the dynamics of individual arrivals can strongly influence the apparent arrival times of the resultant compound waveforms. In acoustic full‐waveform inversion, we do not normally know the density of the seabed, and we do not properly account for finite shear velocity, finite attenuation, and fine‐scale anisotropy variation, all of which can influence the relative amplitudes of different interfering arrivals, which in their turn influence the apparent kinematics. Here, we demonstrate that the introduction of a non‐physical offset‐variable water density during acoustic full‐waveform inversion of this ocean‐bottom‐cable field dataset can compensate efficiently and heuristically for these inaccuracies. This approach improves the travel‐time match and consequently increases both the accuracy and resolution of the final velocity model that is obtained using purely acoustic full‐waveform inversion at minimal additional cost.
-
-
-
3‐D shallow‐water seismic survey evaluation and design using the focal‐beam method: a case study offshore Abu Dhabi
Authors Tomohide Ishiyama and Gerrit BlacquièreABSTRACTFor 3‐D shallow‐water seismic surveys offshore Abu Dhabi, imaging the target reflectors requires high resolution. Characterization and monitoring of hydrocarbon reservoirs by seismic amplitude‐versus‐offset techniques demands high pre‐stack amplitude fidelity. In this region, however, it still was not clear how the survey parameters should be chosen to satisfy the required data quality. To answer this question, we applied the focal‐beam method to survey evaluation and design. This subsurface‐ and target‐oriented approach enables quantitative analysis of attributes such as the best achievable resolution and pre‐stack amplitude fidelity at a fixed grid point in the subsurface for a given acquisition geometry at the surface. This method offers an efficient way to optimize the acquisition geometry for maximum resolution and minimum amplitude‐versus‐offset imprint. We applied it to several acquisition geometries in order to understand the effects of survey parameters such as the four spatial sampling intervals and apertures of the template geometry. The results led to a good understanding of the relationship between the survey parameters and the resulting data quality and identification of the survey parameters for reflection imaging and amplitude‐versus‐offset applications.
-
-
-
Normal moveout velocity for pure‐mode and converted waves in layered orthorhombic medium
Authors Igor Ravve and Zvi KorenABSTRACTWe study the azimuthally dependent hyperbolic moveout approximation for small angles (or offsets) for quasi‐compressional, quasi‐shear, and converted waves in one‐dimensional multi‐layer orthorhombic media. The vertical orthorhombic axis is the same for all layers, but the azimuthal orientation of the horizontal orthorhombic axes at each layer may be different. By starting with the known equation for normal moveout velocity with respect to the surface‐offset azimuth and applying our derived relationship between the surface‐offset azimuth and phase‐velocity azimuth, we obtain the normal moveout velocity versus the phase‐velocity azimuth. As the surface offset/azimuth moveout dependence is required for analysing azimuthally dependent moveout parameters directly from time‐domain rich azimuth gathers, our phase angle/azimuth formulas are required for analysing azimuthally dependent residual moveout along the migrated local‐angle‐domain common image gathers. The angle and azimuth parameters of the local‐angle‐domain gathers represent the opening angle between the incidence and reflection slowness vectors and the azimuth of the phase velocity ψphs at the image points in the specular direction. Our derivation of the effective velocity parameters for a multi‐layer structure is based on the fact that, for a one‐dimensional model assumption, the horizontal slowness and the azimuth of the phase velocity ψphs remain constant along the entire ray (wave) path. We introduce a special set of auxiliary parameters that allow us to establish equivalent effective model parameters in a simple summation manner. We then transform this set of parameters into three widely used effective parameters: fast and slow normal moveout velocities and azimuth of the slow one. For completeness, we show that these three effective normal moveout velocity parameters can be equivalently obtained in both surface‐offset azimuth and phase‐velocity azimuth domains.
-
-
-
Modified image algorithm to simulate seismic channel waves in 3D tunnel model with rugged free surfaces
Authors Hui Li, Peimin Zhu, Guangzhong Ji and Qiang ZhangABSTRACTThe existence of rugged free‐surface three‐dimensional tunnel conditions in the coal seams, caused either by geological or mining processes, will inevitably influence wave propagation characteristics when the seismic waves go through the coal mines. Thus, a modified image algorithm has been developed to account for seismic channel waves propagating through this complicated topography with irregular free surfaces. Moreover, the seismic channel waves commonly exhibit damped and dispersive signatures, which is not only because of their own unique sandwich geometry of rock–coal–rock but also because of the viscoelastic behavior of coal. Considering the complexity of programming in three‐dimensional tunnel models with rugged free surfaces, an optimized vacuum grid search algorithm, enabling to model highly irregular topography and to compute efficiently, is also proposed when using high‐order staggered finite‐difference scheme to simulate seismic channel wave propagations in viscoelastic media. The numerical simulations are implemented to investigate the accuracy and stability of the method and the impact of coal's viscoelastic behavior on seismic channel wave propagation characteristics. The results indicate that the automatic vacuum grid search algorithm can be easily merged into high‐order staggered finite‐difference scheme, which can efficiently be applied to calculate three‐dimensional tunnel models with rugged free surfaces in the viscoelastic media. The simulation also suggests that the occurrence of a three‐dimensional tunnel with free surfaces has a remarkable influence on the seismic channel wave propagation characteristics and elastic energy distribution.
-
-
-
Review paper: Exploration geophysics for intrusion‐hosted rare metals
Authors M.D. Thomas, K.L. Ford and P. KeatingABSTRACTIgneous intrusions, notably carbonatitic–alkalic intrusions, peralkaline intrusions, and pegmatites, represent significant sources of rare‐earth metals. Geophysical exploration for and of such intrusions has met with considerable success. Examples of the application of the gravity, magnetic, and radiometric methods in the search for rare metals are presented and described. Ground gravity surveys defining small positive gravity anomalies helped outline the shape and depth of the Nechalacho (formerly Lake) deposit within the Blatchford Lake alkaline complex, Northwest Territories, and of spodumene‐rich mineralization associated with the Tanco deposit, Manitoba, within the hosting Tanco pegmatite. Based on density considerations, the bastnaesite‐bearing main ore body within the Mountain Pass carbonatite, California, should produce a gravity high similar in amplitude to those associated with the Nechalacho and Tanco deposits. Gravity also has utility in modelling hosting carbonatite intrusions, such as the Mount Weld intrusion, Western Australia, and Elk Creek intrusion, Nebraska.
The magnetic method is probably the most successful geophysical technique for locating carbonatitic–alkalic host intrusions, which are typically characterized by intense positive, circular to sub‐circular, crescentic, or annular anomalies. Intrusions found by this technique include the Mount Weld carbonatite and the Misery Lake alkali complex, Quebec. Two potential carbonatitic–alkalic intrusions are proposed in the Grenville Province of Eastern Quebec, where application of an automatic technique to locate circular magnetic anomalies identified several examples. Two in particular displayed strong similarities in magnetic pattern to anomalies accompanying known carbonatitic or alkalic intrusions hosting rare‐metal mineralization and are proposed to have a similar origin.
Discovery of carbonatitic–alkalic hosts of rare metals has also been achieved by the radiometric method. The Thor Lake group of rare‐earth metal deposits, which includes the Nechalacho deposit, were found by follow‐up investigations of strong equivalent thorium and uranium peaks defined by an airborne survey. Prominent linear radiometric anomalies associated with glacial till in the Canadian Shield have provided vectors based on ice flow directions to source intrusions. The Allan Lake carbonatite in the Grenville Province of Ontario is one such intrusion found by this method. Although not discovered by its radiometric characteristics, the Strange Lake alkali intrusion on the Quebec–Labrador border is associated with prominent linear thorium and uranium anomalies extending at least 50 km down ice from the intrusion. Radiometric exploration of rare metals hosted by pegmatites is evaluated through examination of radiometric signatures of peraluminous pegmatitic granites in the area of the Tanco pegmatite.
-
-
-
Estimating azimuthal stress‐induced P‐wave anisotropy from S‐wave anisotropy using sonic log or vertical seismic profile data
Authors Olivia Collet, Boris Gurevich and Guy DuncanABSTRACTMost sedimentary rocks are anisotropic, yet it is often difficult to accurately incorporate anisotropy into seismic workflows because analysis of anisotropy requires knowledge of a number of parameters that are difficult to estimate from standard seismic data. In this study, we provide a methodology to infer azimuthal P‐wave anisotropy from S‐wave anisotropy calculated from log or vertical seismic profile data. This methodology involves a number of steps. First, we compute the azimuthal P‐wave anisotropy in the dry medium as a function of the azimuthal S‐wave anisotropy using a rock physics model, which accounts for the stress dependency of seismic wave velocities in dry isotropic elastic media subjected to triaxial compression. Once the P‐wave anisotropy in the dry medium is known, we use the anisotropic Gassmann equations to estimate the anisotropy of the saturated medium. We test this workflow on the log data acquired in the North West Shelf of Australia, where azimuthal anisotropy is likely caused by large differences between minimum and maximum horizontal stresses. The obtained results are compared to azimuthal P‐wave anisotropy obtained via orthorhombic tomography in the same area. In the clean sandstone layers, anisotropy parameters obtained by both methods are fairly consistent. In the shale and shaly sandstone layers, however, there is a significant discrepancy between results since the stress‐induced anisotropy model we use is not applicable to rocks exhibiting intrinsic anisotropy. This methodology could be useful for building the initial anisotropic velocity model for imaging, which is to be refined through migration velocity analysis.
-
-
-
Field testing of modular borehole monitoring with simultaneous distributed acoustic sensing and geophone vertical seismic profiles at Citronelle, Alabama
Authors T.M. Daley, D.E. Miller, K. Dodds, P. Cook and B.M. FreifeldABSTRACTA modular borehole monitoring concept has been implemented to provide a suite of well‐based monitoring tools that can be deployed cost effectively in a flexible and robust package. The initial modular borehole monitoring system was deployed as part of a CO2 injection test operated by the Southeast Regional Carbon Sequestration Partnership near Citronelle, Alabama. The Citronelle modular monitoring system transmits electrical power and signals, fibre‐optic light pulses, and fluids between the surface and a reservoir. Additionally, a separate multi‐conductor tubing‐encapsulated line was used for borehole geophones, including a specialized clamp for casing clamping with tubing deployment. The deployment of geophones and fibre‐optic cables allowed comparison testing of distributed acoustic sensing. We designed a large source effort (>64 sweeps per source point) to test fibre‐optic vertical seismic profile and acquired data in 2013. The native measurement in the specific distributed acoustic sensing unit used (an iDAS from Silixa Ltd) is described as a localized strain rate. Following a processing flow of adaptive noise reduction and rebalancing the signal to dimensionless strain, improvement from repeated stacking of the source was observed. Conversion of the rebalanced strain signal to equivalent velocity units, via a scaling by local apparent velocity, allows quantitative comparison of distributed acoustic sensing and geophone data in units of velocity. We see a very good match of uncorrelated time series in both amplitude and phase, demonstrating that velocity‐converted distributed acoustic sensing data can be analyzed equivalent to vertical geophones. We show that distributed acoustic sensing data, when averaged over an interval comparable to typical geophone spacing, can obtain signal‐to‐noise ratios of 18 dB to 24 dB below clamped geophones, a result that is variable with noise spectral amplitude because the noise characteristics are not identical. With vertical seismic profile processing, we demonstrate the effectiveness of downgoing deconvolution from the large spatial sampling of distributed acoustic sensing data, along with improved upgoing reflection quality. We conclude that the extra source effort currently needed for tubing‐deployed distributed acoustic sensing vertical seismic profile, as part of a modular monitoring system, is well compensated by the extra spatial sampling and lower deployment cost as compared with conventional borehole geophones.
-
-
-
Properties of low‐frequency trapped mode in viscous‐fluid waveguides
Authors Anatoly A. Nikitin, Boris D. Plyushchenkov and Arkady Yu. SegalABSTRACTWe derived the velocity and attenuation of a generalized Stoneley wave being a symmetric trapped mode of a layer filled with a Newtonian fluid and embedded into either a poroelastic or a purely elastic rock. The dispersion relation corresponding to a linearized Navier–Stokes equation in a fracture coupling to either Biot or elasticity equations in the rock via proper boundary conditions was rigorously derived. A cubic equation for wavenumber was found that provides a rather precise analytical approximation of the full dispersion relation, in the frequency range of 10−3 Hz to 103 Hz and for layer width of less than 10 cm and fluid viscosity below 0.1 Pa· s [100 cP]. We compared our results to earlier results addressing viscous fluid in either porous rocks with a rigid matrix or in a purely elastic rock, and our formulae are found to better match the numerical solution, especially regarding attenuation. The computed attenuation was used to demonstrate detectability of fracture tip reflections at wellbore, for a range of fracture lengths and apertures, pulse frequencies, and fluid viscosity.
-
-
-
On the determination of marine magnetotelluric receiver orientation using land magnetic observatory data
More LessABSTRACTMarine magnetotelluric measurements using “free‐fall’’ instruments without effective compasses suffer from the problem of unknown orientation of the receivers at the seafloor. While past works indicate that marine magnetotelluric orientation of the instruments can be estimated by reference to land deployments of known orientation using the transfer tensor method, there is limited published information on how this is implemented in practice. We document this method and propose a set of new time‐ and frequency‐domain approaches to solve this orientation problem of the seafloor receivers. We test these methodologies in onshore and offshore magnetotelluric data whose orientations are well known and apply these techniques to marine magnetotelluric data with unknown orientation. For the controlled tests, both time‐ and frequency‐domain approaches produce overall comparable results. To investigate the effects of the subsurface structure distribution on the orientation process, a dimensionality analysis of a controlled dataset is carried out. In subsequent analysis using the available disoriented marine magnetotelluric data from offshore Brazil and from the Vassouras magnetic observatory on the mainland for remote referencing, frequency‐domain methods yield approximate orientation angles among themselves with low standard deviation each. Time‐domain results are consistent for most cases but differ from frequency‐domain results for some situations.
-
-
-
From noise correlation to resonant ultrasound spectroscopy in rock acoustics
ABSTRACTThis article introduces an alternative experimental procedure for measuring the elastic properties of a solid material at laboratory scale, using both the principles of passive seismic interferometry and resonance ultrasound spectroscopy. We generate noise into the studied sample with a pneumatic air blow gun, and we cross‐correlate the signals recorded with two passive piezoelectric sensors put in soft contact with the sample surface. Resonance phenomena are induced in the sample, but in contrast with conventional resonance ultrasound spectroscopy experiments, we have no control over the injected frequencies that are sent all together within the noise spectrum. The spectrum of the correlogram is a good approximation of the resonance spectrum of the vibrating sample and can be inverted in terms of the elastic moduli of the constituent material of the sample.
The experimental procedure is validated on samples made of standard materials (here, aluminium and Plexiglas) by consistently comparing the inverted elastic velocities with the velocities independently measured with the conventional technique of ultrasonic pulse transmission. Moreover, we got similar positive results on dry rock samples, such as Vilhonneur limestone. These encouraging preliminary results open up promising prospects for monitoring fluid substitution in rock samples using the technique described in this paper.
-
-
-
Research Note: Compact Depth from Extreme Points: a tool for fast potential field imaging
Authors Jamaledin Baniamerian, Maurizio Fedi and Behrooz OskooiABSTRACTWe propose a fast method for imaging potential field sources. The new method is a variant of the “Depth from Extreme Points,” which yields an image of a quantity proportional to the source distribution (magnetization or density). Such transformed field is here transformed into source‐density units by determining a constant with adequate physical dimension by a linear regression of the observed field versus the field computed from the “Depth from Extreme Points” image. Such source images are often smooth and too extended, reflecting the loss of spatial resolution for increasing altitudes. Consequently, they also present too low values of the source density. We here show that this initial image can be improved and made more compact to achieve a more realistic model, which reproduces a field consistent with the observed one. The new algorithm, which is called “Compact Depth from Extreme Points” iteratively produces different source distributions models, with an increasing degree of compactness and, correspondingly, increasing source‐density values. This is done through weighting the model with a compacting function. The compacting function may be conveniently expressed as a matrix that is modified at any iteration, based on the model obtained in the previous step. At any iteration step the process may be stopped when the density reaches values higher than prefixed bounds based on known or assumed geological information. As no matrix inversion is needed, the method is fast and allows analysing massive datasets. Due to the high stability of the “Depth from Extreme Points” transformation, the algorithm may be also applied to any derivatives of the measured field, thus yielding an improved resolution. The method is investigated by application to 2D and 3D synthetic gravity source distributions, and the imaged sources are a good reconstruction of the geometry and density distributions of the causative bodies. Finally, the method is applied to microgravity data to model underground crypts in St. Venceslas Church, Tovacov, Czech Republic.
-
-
-
Research Note: The mean sensitivity depth of the electrical resistivity method
By S.L. ButlerABSTRACTUntil now, a simple formula to estimate the depth of investigation of the electrical resistivity method that takes into account the positions of all of the electrodes for a general four‐electrode array has not been available. While the depth sensitivity function of the method for a homogeneous infinite half‐space is well known, previous attempts to use it to characterize the depth of investigation have involved calculating its peak and median, both of which must be determined numerically for a general four‐electrode array. I will show that the mean of the sensitivity function, which has not been considered previously, does admit a very simple mathematical formula. I compare the mean depth with the median and peak sensitivity depths for some common arrays. The mean is always greater than or equal to the median that is always greater than the peak. All three measures give reasonable estimates to the depths of actual structures for most circumstances. I will further show that, for 1D soundings, the use of the mean sensitivity depth as the pseudo‐depth assigns an apparent resistivity to a given pseudo‐depth that is consistent between different arrays. One consequence of this is that smoother depth soundings are obtained as “clutches,” caused by a change in the depth sensitivity due to moving the potential electrodes, are effectively removed. I expect that the mean depth formula will be a useful “rule of thumb” for estimating the depth of investigation before the resistivity structure of the ground is known.
-
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)