- Home
- A-Z Publications
- Geophysical Prospecting
- Previous Issues
- Volume 66, Issue 2, 2018
Geophysical Prospecting - Volume 66, Issue 2, 2018
Volume 66, Issue 2, 2018
-
-
Post‐stack diffraction imaging in vertical transverse isotropy media using non‐hyperbolic moveout approximations
Authors M. Javad Khoshnavaz, Andrej Bóna and Milovan UrosevicABSTRACTDiffractions carry valuable information about local discontinuities and small‐scale objects in the subsurface. They are still not commonly used in the process of geological interpretation. Many diffraction imaging techniques have been developed and applied for isotropic media, whereas relatively few techniques have been developed for anisotropic media. Ignoring anisotropy can result in low‐resolution images with wrongly positioned or spurious diffractors. In this article, we suggest taking anisotropy into account in two‐dimensional post‐stack domain by considering P‐wave non‐hyperbolic diffraction traveltime approximations for vertical transverse isotropy media, previously developed for reflection seismology. The accuracy of the final images is directly connected to the accuracy of the diffraction traveltime approximations. We quantified the accuracy of six different approximations, including hyperbolic moveout approximation, by the application of a post‐stack diffraction imaging technique on two‐dimensional synthetic data examples.
-
-
-
Modelling and straight‐ray tomographic imaging studies of cross‐hole radio‐frequency electromagnetic data for mineral exploration
Authors Yongxing Li and Richard S. SmithABSTRACTRadio‐frequency electromagnetic tomography (or radio imaging method) employs radio‐frequency (typically 0.1–10 MHz) electromagnetic wave propagation to delineate the distribution of electric properties between two boreholes. Currently, the straight‐ray imaging method is the primary imaging method for the radio imaging method data acquired for mineral exploration. We carried out synthetic studies using three‐dimensional finite‐element modelling implemented in COMSOL Multiphysics to study the electromagnetic field characteristics and to assess the capability of the straight‐ray imaging method using amplitude and phase data separately. We studied four sets of experiments with models of interest in the mining setting. In the first two experiments, we studied models with perfect conductors in homogeneous backgrounds, which show that the characteristics of the electromagnetic fields depend mainly on the wavelength. When the borehole separations are less than one wavelength, induction effects occur; conductors with simple geometries can be recovered acceptably with amplitude data but are incorrectly imaged on the phase tomogram. When the borehole separations are longer than two wavelengths, radiation effects play a major role. In this case, phase tomography provides images with acceptable quality, whereas amplitude tomography does not provide satisfactory results. The third experiment shows that imaging with both original and reciprocal datasets is somewhat helpful in improving the imaging quality by reducing the impact of noise. In the last experiment, we studied models with conductive zones extended into the borehole plane with different lengths, which were not accurately recovered with amplitude tomography. The experiment implies that it is difficult to determine the extent of a mineralised zone that has been intersected by one of the boreholes. Due to the large variation of the wavelength in the radio‐frequency range, we suggest investigating the local electric properties to select an operating frequency prior to a survey. We conclude that straight‐ray tomography with either amplitude or phase data cannot provide high‐quality imaging results. We suggest using more general methods based on full electromagnetic modelling to interpret the data. In circumstances when computational time is critical, we suggest saving time by using either induction methods for borehole separations less than one wavelength or wave‐based methods (only radiation fields are considered) for borehole separation larger than two wavelengths.
-
-
-
The design of visco‐acoustic weighted L1‐error frequency–space wavefield extrapolators
Authors Syed Abdul Salam and Wail A. MousaABSTRACTSeismic imaging is an important step for imaging the subsurface structures of the Earth. One of the attractive domains for seismic imaging is explicit frequency–space (f – x) prestack depth migration. So far, this domain focused on migrating seismic data in acoustic media, but very little work assumed visco‐acoustic media. In reality, seismic exploration data amplitudes suffer from attenuation. To tackle the problem of attenuation, new operators are required, which compensates for it. We propose the weighted L1‐error minimisation technique to design visco‐acoustic f – x wavefield extrapolators. The L1‐error wavenumber responses provide superior extrapolator designs as compared with the previously designed equiripple L4‐norm and L∞‐norm extrapolation wavenumber responses. To verify the new compensating designs, prestack depth migration is performed on the challenging Marmousi model dataset. A reference migrated section is obtained using non‐compensating f – x extrapolators on an acoustic dataset. Then, both compensating and non‐compensating extrapolators are applied to a visco‐acoustic dataset, and both migrated sections are then compared. The final images show that the proposed weighted L1‐error method enhances the resolution and results in practically stable images.
-
-
-
Wavefield interpolation in 3D large‐step Fourier wavefield extrapolation
Authors Gaoxiang Chen, Li‐Yun Fu, Wei Wei and Weijia SunABSTRACTExtrapolating wavefields and imaging at each depth during three‐dimensional recursive wave‐equation migration is a time‐consuming endeavor. For efficiency, most commercial techniques extrapolate wavefields through thick slabs followed by wavefield interpolation within each thick slab. In this article, we develop this strategy by associating more efficient interpolators with a Fourier‐transform‐related wavefield extrapolation method. First, we formulate a three‐dimensional first‐order separation‐of‐variables screen propagator for large‐step wavefield extrapolation, which allows for wide‐angle propagations in highly contrasting media. This propagator significantly improves the performance of the split‐step Fourier method in dealing with significant lateral heterogeneities at the cost of only one more fast Fourier transform in each thick slab. We then extend the two‐dimensional Kirchhoff and Born–Kirchhoff local wavefield interpolators to three‐dimensional cases for each slab. The three‐dimensional Kirchhoff interpolator is based on the traditional Kirchhoff formula and applies to moderate lateral velocity variations, whereas the three‐dimensional Born–Kirchhoff interpolator is derived from the Lippmann–Schwinger integral equation under the Born approximation and is adapted to highly laterally varying media. Numerical examples on the three‐dimensional salt model of the Society of Exploration Geophysicists/European Association of Geoscientists demonstrate that three‐dimensional first‐order separation‐of‐variables screen propagator Born–Kirchhoff depth migration using thick‐slab wavefield extrapolation plus thin‐slab interpolation tolerates a considerable depth‐step size of up to 72 ms, eventually resulting in an efficiency improvement of nearly 80% without obvious loss of imaging accuracy. Although the proposed three‐dimensional interpolators are presented with one‐way Fourier extrapolation methods, they can be extended for applications to general migration methods.
-
-
-
An upwind fast sweeping scheme for calculating seismic wave first‐arrival travel times for models with an irregular free surface
Authors Haiqiang Lan and Ling ChenABSTRACTThe topography‐dependent eikonal equation formulated in a curvilinear coordinate system has recently been established and revealed as being effective in calculating first‐arrival travel times of seismic waves in an Earth model with an irregular free surface. The Lax–Friedrichs sweeping scheme, widely used in previous studies as for approximating the topography‐dependent eikonal equation viscosity solutions, is more dissipative and needs a much higher number of iterations to converge. Furthermore, the required number of iterations grows with the grid refinement and results in heavy computation in dense grids, which hampers the application of the Lax–Friedrichs sweeping scheme to seismic wave travel‐time calculation and high‐resolution imaging. In this paper, we introduce a new upwind fast sweeping solver by discretising the Legendre transform of the numerical Hamiltonian of the topography‐dependent eikonal equation using an explicit formula. The minimisation related to the Legendre transform in the sweeping scheme is solved analytically, which proved to be much more efficient than the Lax–Friedrichs algorithm in solving the topography‐dependent eikonal equation. Several numerical experiments demonstrate that the new upwind fast sweeping method converges and achieves much better accuracy after a finite number of iterations, independently of the mesh size, which makes it an efficient and robust tool for calculating travel times in the presence of a non‐flat free surface.
-
-
-
Detecting underground cavities using microtremor data: physical modelling and field experiment
Authors Yu.I. Kolesnikov and K.V. FedinABSTRACTIn this paper, we study the possibilities of the use of microtremor records in the detection and delineation of near‐surface underground cavities. Three‐dimensional physical modelling data showed that the averaging amplitude spectra of a large number of microtremor records makes it possible to evaluate the frequencies and amplitudes of compressional standing waves generated by microtremor in the space between the ground surface and underground cavities. We illustrate how these parameters can be used to estimate the shape of the underground cavity horizontal projection. If the compressional wave velocity in the enclosing rock is known, it is possible to evaluate the depth to the cavity top using the frequencies of the standing waves. The results of the field experiment confirmed the possibility of underground cavities detection using microtremor data.
-
-
-
Practical and robust experimental determination of c13 and Thomsen parameter δ
Authors Fuyong Yan, De‐Hua Han and Xue‐Lian ChenABSTRACTWe analysed the complications in laboratory velocity anisotropy measurement on shales. There exist significant uncertainties in the laboratory determination of c13 and Thomsen parameter δ. These uncertainties are primarily related to the velocity measurement in the oblique direction. For reliable estimation of c13 and δ, it is important that genuine phase velocity or group velocity be measured with minimum uncertainty. The uncertainties can be greatly reduced if redundant oblique velocities are measured. For industrial applications, it is impractical to make multiple oblique velocity measurements on multiple core plugs. We demonstrated that it is applicable to make multiple genuine oblique group velocity measurements on a single horizontal core plug. The measurement results show that shales can be classified as a typical transversely isotropic medium. There is a coupling relation between c44 and c13 in determining the directional dependence of the seismic velocities. The quasi‐P‐wave or quasi‐S‐wave velocities can be approximated by three elastic parameters.
-
-
-
Basement structural architecture and hydrocarbon conduit potential of polygonal faults in the Cooper‐Eromanga Basin, Australia
Authors David Kulikowski, Khalid Amrouch, Dennis Cooke and Michael Edward GrayABSTRACTA thorough and complete understanding of the structural geology and evolution of the Cooper‐Eromanga Basin has been hampered by low‐resolution seismic data that becomes particularly difficult to interpret below the thick Permian coal measures. As a result, researchers are tentative to interpret the basement fault architecture within the basin, which is largely undefined. To provide a better understanding of the basement fault geometry, all available two‐dimensional seismic lines together with 12 three‐dimensional seismic surveys were structurally interpreted with assistance from seismic attribute analysis. The Upper Cretaceous Cadna‐owie Formation and top Permian reflectors were analysed using a common seismic attribute technique (incoherency) that was used to infer the presence of faults that may have otherwise been overlooked. Detailed basement fault maps for each seismic survey were constructed and used in conjunction with two‐dimensional seismic data interpretation to produce a regional basement fault map. Large north‐northeast–south‐southwest‐striking sinistral strike–slip faults were identified within the Patchawarra Trough appearing to splay from the main northeast–southwest‐striking ridge. These sinistral north‐northeast–south‐southwest‐striking faults, together with field‐scale southeast–northwest‐striking dextral strike–slip faults, are optimally oriented to have potentially developed as a conjugated fault set under a south‐southeast–north‐northwest‐oriented strike–slip stress regime. Geomechanical modelling for a regionally extensive system of Cretaceous polygonal faults was performed to calculate the Leakage Factor and Dilation Tendency of individual faults. Faults that extend into Lower Cretaceous oil‐rich reservoirs with strikes of between 060°N and 140°N and a high to near‐vertical dip angle were identified to most likely be acting as conduits for the tertiary migration of hydrocarbons from known Lower Cretaceous hydrocarbon reservoirs into shallow Cretaceous sediments. This research provides valuable information on the regional basement fault architecture and a more detailed exploration target for the Cooper‐Eromanga Basin, which were previously not available in literature.
-
-
-
3D joint inversion of magnetotelluric and airborne tipper data: a case study from the Morrison porphyry Cu–Au–Mo deposit, British Columbia, Canada
More LessABSTRACTZ‐axis tipper electromagnetic and broadband magnetotelluric data were used to determine three‐dimensional electrical resistivity models of the Morrison porphyry Cu–Au–Mo deposit in British Columbia. Z‐axis tipper electromagnetic data are collected with a helicopter, thus allowing rapid surveys with uniform spatial sampling. Ground‐based magnetotelluric surveys can achieve a greater exploration depth than Z‐axis tipper electromagnetic surveys, but data collection is slower and can be limited by difficult terrain. The airborne Z‐axis tipper electromagnetic tipper data and the ground magnetotelluric tipper data show good agreement at the Morrison deposit despite differences in the data collection method, spatial sampling, and collection date. Resistivity models derived from individual inversions of the Z‐axis tipper electromagnetic tipper data and magnetotelluric impedance data contain some similar features, but the Z‐axis tipper electromagnetic model appears to lack resolution below a depth of 1 km, and the magnetotelluric model suffers from non‐uniform and relatively sparse spatial sampling. The joint Z‐axis tipper electromagnetic inversion solves these issues by combining the dense spatial sampling of the airborne Z‐axis tipper electromagnetic technique and the deeper penetration of the lower frequency magnetotelluric data. The resulting joint resistivity model correlates well with the known geology and distribution of alteration at the Morrison deposit. Higher resistivity is associated with the potassic alteration zone and volcanic country rocks, whereas areas of lower resistivity agree with known faults and sedimentary units. The pyrite halo and ≥0.3% Cu zone have the moderate resistivity that is expected of disseminated sulphides. The joint Z‐axis tipper electromagnetic inversion provides an improved resistivity model by enhancing the lateral and depth resolution of resistivity features compared with the individual Z‐axis tipper electromagnetic and magnetotelluric inversions. This case study shows that a joint Z‐axis tipper electromagnetic–magnetotelluric approach effectively images the interpreted mineralised zone at the Morrison deposit and could be beneficial in exploration for disseminated sulphides at other porphyry deposits.
-
-
-
Current and voltage source induced polarization transients: a comparative consideration
Authors N.O. Kozhevnikov and E.Y. AntonovABSTRACTThe article discusses the excitation of transient induced polarization responses using current and voltage sources. The first method has found a wide application in induced polarization surveys and—directly or indirectly—in the theory of the induced polarization method. Typically, rectangular current pulses are injected into the earth via grounding electrodes, and decaying induced polarization voltage is measured during the pauses between pulses. In this case, only the secondary field is recorded in the absence of the primary field, which is an important advantage of this method. On the other hand, since the current injected into the ground is fully controlled by the source, this method does not allow studying induced polarization by measuring the current in the transmitter line or associated magnetic field. When energising the earth with voltage pulses, the measured quantity is the transient induced polarization current. In principle, this method allows induced polarization studies to be done by recording the transmitter line current, the associated magnetic field, or its rate of change. The decay of current in a grounded transmitter line depends not only on the induced polarization of the earth but also on the polarization of the grounding electrodes. This problem does not occur when induced polarization transients in the earth are excited inductively. A grounded transmitter line is a mixed‐type source; hence, for a purely inductive excitation of induced polarization transients, one should use an ungrounded loop, which is coupled to the earth solely by electromagnetic induction.
-
-
-
Depth from extreme points method for gravity gradient tensor data
Authors Wenna Zhou and Chengying LiuABSTRACTThe ‘depth from extreme points’ method is an important tool to estimate the depth of sources of gravity and magnetic data. In order to interpret gravity gradient tensor data conveniently, formulas for the tensor data form regarding depth from the extreme points method were calculated in this paper. Then, all of the gradient tensor components were directly used to interpret the causative source. Beyond the gzz component, also the gxx and gyy components can be used to obtain depth information. In addition, the total horizontal derivative of the depth from extreme points of the gradient tensor can be used to describe the edge information of geologic sources. In this paper, we investigated the consistency of the homogeneity degree calculated by using the different components, which leads to the calculated depth being confirmed. Therefore, a more integrated interpretation can be obtained by using the gradient tensor components. Different synthetic models were used with and without noise to test the new approach, showing stability, accuracy and speed. The proposed method proved to be a useful tool for gradient tensor data interpretation. Finally, the proposed method was applied to full tensor gradient data acquired over the Vinton Salt Dome, Louisiana, USA, and the results are in agreement with those obtained in previous research studies.
-
-
-
Analysing two‐dimensional effects in central loop transient electromagnetic sounding data using a semi‐synthetic tipper approach
Authors P. Yogeshwar and B. TezkanABSTRACTWe present a simple and feasible approach to analyse and identify two‐dimensional effects in central loop transient electromagnetic sounding data and the correspondingly derived quasi two‐dimensional conductivity models. The proposed strategy is particularly useful in minimising interpretation errors. It is based on the calculation of a semi‐synthetic transient electromagnetic tipper at each sounding and for each observational transient time point. The semi‐synthetic transient electromagnetic tipper is derived from the measured vertical component of the induced voltage and the synthetically calculated horizontal component. The approach is computationally inexpensive and involves one two‐dimensional forward calculation of an obtained quasi two‐dimensional conductivity section. Based on a synthetic example, we demonstrate that the transient electromagnetic tipper approach is applicable in identifying which transient data points and which corresponding zones in a derived quasi two‐dimensional subsurface model are affected by two‐dimensional inhomogeneities. The one‐dimensional inversion of such data leads to false models. An application of the semi‐synthetic transient electromagnetic tipper to field data from the Azraq basin in Jordan reveals that, in total, eight of 80 investigated soundings are affected by two‐dimensional structures although the field data can be fitted optimally using one‐dimensional inversion techniques. The largest semi‐synthetic tipper response occurs in a 300 m‐wide region around a strong lateral resistivity contrast. The approach is useful for analysing structural features in derived quasi two‐dimensional sections and for qualitatively investigating how these features affect the transient response. To avoid misinterpretation, these identified zones corresponding to large tipper values are excluded from the interpretation of a quasi two‐dimensional conductivity model. Based on the semi‐synthetic study, we also demonstrate that a quantitative interpretation of the horizontal voltage response (e.g. by inversion) is usually not feasible as it requires the exact sensor position to be known. Although a tipper derived purely from field data is useful as a qualitative tool for identifying two‐dimensional distortion effects, it is only feasible if the sensor setup is sufficiently accurate. Our proposed semi‐synthetic transient electromagnetic tipper approach is particularly feasible as an a posteriori approach if no horizontal components are recorded or if the sensor setup in the field is not sufficiently accurate.
-
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)