Volume 48 Number 4 2017
EG15077Relationship between bulk mineralogy and induced polarisation responses in iron oxide-copper-gold and porphyry copper mineralisation, northern Chile
We have studied the correlation between bulk mineralogy and induced polarisation (IP) responses in iron oxide-copper-gold (IOCG) and porphyry copper concealed mineralised systems in northern Chile. Data analyses show that bulk sulphide concentration is the first order factor that explains variation in the IP effect, even for the case of magnetite-related IOCG mineralisation.
EG15074Divergence correction schemes in finite difference method for 3D tensor CSAMT in axial anisotropic media
This paper presents a study on divergence correction in the finite difference method for 3D tensor controlled-source audio-frequency magnetotellurics (CSAMT) in an axial anisotropic media. We show how divergence correction can be performed more efficiently and the study results can be used to help accelerate a whole 3D tensor CSAMT forward modelling and inversion.
EG15095Short-offset grounded-wire TEM method for efficient detection of mined-out areas in vegetation-covered mountainous coalfields
The objective of this study was to overcome the limitations in using the loop-source TEM method and realise the efficient detection of mined-out areas using the short-offset grounded-wire TEM (SOTEM) method in vegetation-covered mountainous coalfields. The feasibility of SOTEM is verified by the drilling result in the field test.
EG15059Ground resistivity method and DCIP2D forward and inversion modelling to identify alteration at the Midwest uranium deposit, northern Saskatchewan, Canada
Forward and inverse modelling are used to demonstrate that the resistivity method can detect alteration above the Midwest uranium deposit in Canada. The inverse section of the field data is similar to the inverse section of a synthetic model when alteration is present, but dissimilar when there is no alteration.
EG15078Near-surface high resolution imaging of a metallogenic zone in the northern fringe of Dalma volcanics in eastern India using electrical resistivity tomography
An ERT survey was conducted over a defined conducting zone through a magnetotelluric survey in the Dalma volcanics, India. The resistivity models resolved four major zones and the results corroborate well with the geological succession. The resistivity values of less than 1 Ωm at a depth of ~25 m correspond to the deposit.
We obtained the shallow S-wave velocity profile distributions in the Bandung Basin, Indonesia, from the inversion of the observed Rayleigh wave phase velocity spectra using a microtremor survey. The profiles will be used for calculation of the average S-wave velocity of the upper 30 m and site amplification factors for seismic hazard evaluation in the area.
Subsurface-property estimation from remote geophysical measurements is subject to non-Gaussian errors. We propose to minimise the misfit with a robust error measure, which is based on a generalised Gaussian distribution. A suboptimal solution is proposed through a mixed-norm functional combination of the l1 and l2 norms.
EG15097Simultaneous seismic interpolation and denoising based on sparse inversion with a 3D low redundancy curvelet transform
This paper introduces a low redundancy curvelet transform which can reduce the redundancy to 10 for three dimensional data and simultaneously interpolate and denoise. Numerical experiments on synthetic and field data demonstrate that the low redundancy transform can provide reliable results while at the same time improving the computational efficiency.
EG15112Improved reservoir characterisation using fuzzy logic platform: an integrated petrophysical, seismic structural and poststack inversion study
In this study, petrophysical properties, seismic structural and poststack seismic inversion results are integrated using the fuzzy logic AND operator to characterise the Tensleep Sandstone Formation at Powder River Basin, Wyoming, USA. The 3D seismic cube and logging datasets of 17 wells are used for enhancing the seismic and petrophysical evaluation.
Benefitting from the forward modelling studies, the phenomenon of spatial mispositioning in apparent reflections caused by the symmetry axis orientation of the TI overburden was demonstrated. Numerical manipulations show the imaging problems can be well described by the TERF method and properly handled by the RAKTM method.
A new parameterisation using P-wave velocity, Poisson’s ratio and density is proposed for multi-parameter full waveform inversion (FWI) for isotropic elastic media. By introducing Poisson’s ratio, the virtual source for the P-wave velocity generates additional S-waves in the partial derivative wavefields, which increases the spatial resolution of FWI.
EG16031Separation of prestack seismic diffractions using an improved sparse apex-shifted hyperbolic Radon transform
We have developed an improved sparse apex-shifted hyperbolic Radon transform with Fourier-kernel Stolt-based modelling and imaging operators to separate seismic diffractions from reflections. The time axis is stretched to overcome the limitation of the multi-velocity model and a sparsity-promoting inversion is implemented to enhance the resolution of the Radon panel.
EG16020Identification of small-scale discontinuities based on dip-oriented gradient energy entropy coherence estimation
We have proposed a dip-oriented gradient energy entropy (DOGEE) coherence method to detect subtle faults and structural features. Theoretical synthetic datasets and a real dataset show that our proposed method is not only robust to noise, but also improves the clarity of small-scale discontinuities.
EG15107Analysis of crooked-line 2D seismic reflection data recorded in areas with complex surface and subsurface conditions
A crooked-line 2D seismic reflection dataset recorded in the northern part of the Eastern Pannonian Basin was processed to obtain a high-resolution seismic image for the area of interest. Refraction interferometry was used to enhance the first arrivals and to compute static corrections with good accuracy.
Depth variation of scattering attenuation by applying single events to multiple lapse time window analysis (MLTWA) was recently shown; however, this application inevitably combined several events due to insufficient data. This study demonstrated that a flexible range window can be applied successfully to solely single events in MLTWA.
EG15021Radiogenic heat production in Rudeis Formation, Lower Miocene, Belayim marine oil field, Gulf of Suez, Egypt
This study describes radiogenic heat production (RHP) estimated from two sets of well logging data recorded in Rudeis Formation, Belayim marine oil field. Cross-correlation of the estimated RHP (AR and ABR) yielded a moderate fit of 0.75. The three-dimensional (3D)-slicing of the reservoir shows that the northern and north-western parts of the study area have higher RHP compared to other parts.
EG16010The research on the buried public monumental complexes of Lupiae (Lecce) by geophysical prospecting
This paper presents a case study that involves a geophysical survey employing the surface 3D ground penetrating radar (GPR) techniques, in order to archaeologically characterise the areas located in the south-western sector of the historical centre of Lecce (Italy). GPR measurements significantly contributed to reconstruction of the complex subsurface properties in these modern urban areas.
The peer-reviewed and edited version of record published online before inclusion in an issue
Multiple-point geostatistical simulation (MPS) was applied to develop 3D ore models matched to surrounding geological information accompanying aeromagnetic data using training image. The present study proposes a method for reducing the uncertainty of the 3D ore model, applying MPS to create probabilistic ore models and analysing the correlation between the models and geophysical data.
Building on the concepts of cohesion degree and local relaxation, we propose an integrated hierarchical equilibrium parallel finite-element reverse time migration (HEP-FE-RTM) algorithm, which has the distinct advantage in that parallel efficiency does not decrease as the number of processors increases.
EG17014Amplification characteristics at Iedang Reservoir dam sites determined using H/V spectral ratio with background noise, S-wave and coda wave energy
Amplification characteristics are critical in evaluating the reliability in seismic design for engineering and seismic source and crustal attenuation characteristics for seismology. Horizontal-to-vertical spectral ratios were analysed using six seismograms observed at four sites near Iedang Reservoir in Korea. Amplification characteristics were compared using the S-wave, coda wave and background noise of each seismogram.
EG16132Least-squares Kirchhoff migration with non-smooth regularisation strategy for subsurface imaging
A fast non-smooth regularisation method is proposed to solve the difficulties involved in the classical least-squares Kirchhoff migration. It not only accounts for the irregular and incomplete data sampling, but also compensates for the anomalous ray coverage and multipathing problem. Numerical experiments show that the method works well.
EG17068Geophysical strata rating (GSR) as an aid in carbonate reservoir characterisation: an example from the South Pars gas field, Persian Gulf Basin
In this study, the geophysical strata rating (GSR), which is an empirical measure of rock competency, is calculated from petrophysical data. The GSR is then extended to the whole South Pars gas field in the framework of 3D seismic data through an acoustic impedance poststack seismic inversion.
The Marine Vibrator Joint Industry Project (MVJIP), sponsored by Shell, Total and ExxonMobil, has been ongoing for four years. This paper discusses the history, motivation, technology and latest developments from the project.
EG17077Estimation of reservoir properties from seismic data through a Markov Chain Monte Carlo-AVA inversion algorithm
We formulate the amplitude versus angle (AVA) inversion in terms of a Markov Chain Monte Carlo algorithm and apply it for reservoir characterisation and litho-fluid facies prediction in offshore Nile Delta. A blind test, based on available well log information, demonstrates the applicability of the proposed method and the reliability of the results.
In laboratories, core cannot always be cut exactly along the axis of symmetry (normal to the bedding plane), which leads to a minor core cutting rotation error. This paper investigates the sensitivity of Thomson’s anisotropic parameters, epsilon (ε), gamma (γ) and delta (δ), to the cutting error.
EG17139Robust data processing of noisy marine controlled-source electromagnetic data using independent component analysis
We have developed a new controlled-source electromagnetic (CSEM) data processing scheme that uses a robust method based on independent component analysis to extract interpretable datasets from noisy marine CSEM data. The results show that the coherent and environment noise is reduced sufficiently. It makes interpretation of the resulting data possible.
EG17055A constrained spectral inversion method based on compressive sensing in order to distinguish high-quality shale
Distinguishing high-quality shale effectively is a difficult exploration problem. We propose a constrained spectral inversion method based on compressive sensing to handle this challenge. This method can help us improve the resolution and continuity of the profiles and enhance our ability to discover high-quality shale.
EG17048Seismic facies-controlled prestack simultaneous inversion of elastic and petrophysical parameters for favourable reservoir prediction
This paper describes the development of a deterministic prestack inversion for estimating elastic and petrophysical parameters. The Gassmann equation is used to construct the relationship between the seismic data and petrophysical parameters. Seismic facies constraints were introduced to improve the accuracy. The very fast-simulated annealing method is used to quickly find the optimal solutions.
EG17057Channel edge detection using 2D complex shearlet transform: a case study from the South Caspian Sea
Channel detection plays a significant role in seismic interpretation. Shearlet transform as a multi-scale and multi-directional transformation is capable of detecting anisotropic singularities. We applied complex-valued shearlet edge measure to synthetic and real seismic time-slices from the South Caspian Sea. The proposed algorithm outperformed both Sobel and Canny edge detectors.
EG17049Synthetic modelling and analysis of CSEM full-field apparent resistivity response combining EM induction and IP effect for 1D medium
Full-field apparent resistivity response combining EM induction and IP effect for layered earth were modelled using the Cole–Cole model. The influence of chargeability is greater than that of the other two parameters, time constant and frequency dependent coefficient. The layered model was distorted seriously by using the inversion algorithm ignoring IP effect.
This paper demonstrates the magnitude of change that results when different information (relating to calibration factors, time zones and time changes, height, geodetic datums and gravity datums) is introduced at various stages of the gravity reduction process. The differences range from below 0.01 mGal to over 16.0 mGal.
EG16156Estimation of elastic anisotropy from three-component ultrasonic measurements using laser Doppler interferometry
Laser Doppler interferometers were previously employed to detect ultrasonic waves propagating in different directions and to estimate elastic anisotropy from these measurements. Our numerical simulations and laboratory measurements show that the recorded wavefield contains converted PS-waves, which need to be taken into consideration to obtain robust estimates of anisotropy.
Measurements of natural time-fluctuating magnetic and electric fields at the surface of Earth produce magnetotelluric data, most of which can be expressed as quantities invariant to the directions of field measurement. Such invariants give information on geological structure.
EG17129Study of the Nankai seismogenic fault using dynamic wave propagation modelling of digital rock from the Nobeoka Fault
We quantified effective VP, VS and the ratio of VP to VS of the 3D digital rock models with different crack-filling minerals. By comparing the elastic properties derived from the digital rock models with the seismic velocities around the Nankai seismogenic fault, we characterised the evolution process of the seismogenic fault.
EG16154The archaeological site of Sagalassos (Turkey): exploring the mysteries of the invisible layers using geophysical methods
The ancient Roman city of Sagalassos (Turkey) is covered by layers of eroded soil that have preserved many secrets waiting to be revealed. A geophysical campaign was planned to highlight the buried structure. Geophysics revealed evidence of a clay quarry, a number of tombs related to the Byzantine period and defensive walls.
EG17110The relationship between normalised horizontal-to-vertical spectral ratios (HVSRs) of microtremors and the F distribution
This paper shows that a method used to identify a peak of horizontal-to-vertical spectral ratios (HVSRs) of microtremors, which was proposed in 2001, is unreliable, and clarifies the reason for the unreliability. Because of the resemblance between the distribution of normalised HVSRs and the F distribution, we frequently fail to detect the deviation between them.
EG16123A new method to enhance the characterisation of seismically thin beds based on the generalised S transform maximum modulus
Thin beds are common in underground hydrocarbon reservoirs, but rarely detected and characterised. In this work, we propose a new method to improve characterisation of thin beds and we are able to find and characterise more thin beds accurately. This new method will enhance our ability to discover more oil and gas.
In this paper, the correction coefficients of a fluxgate magnetic tensor gradiometer and the magnetic properties of an aircraft are determined. A recursive method is used to combine the magnetic interference model and the error model into a unified calibration model. The method is significant as there will be a greater use of airborne magnetic tensor gradiometers using a wide range of aircraft.
EG16134Velocity analysis of CMP gathers acquired by an array GPR system ‘Yakumo’: results from field application to tsunami deposits
Velocity analysis is applied to explore the capability of a multi-static ground penetrating radar (GPR) system to estimate a vertical profile of dielectric constants of stratified tsunami layers. Both laboratory and field experiments demonstrate that multi-static GPR can accurately delineate the layered geological structures.
In this study, the deficiency of the Kozeny–Carman assumption is investigated and alternative derived equations for the Kozeny–Carman equation are suggested, including equations where the grain size is replaced with the pore size and with varying tortuosity. Relationships for the permeability of shaly sand reservoir are introduced.
EG17104A simple inversion algorithm to estimate a linearly increasing velocity model for microseismic monitoring
In this study, we suggest a simple algorithm to estimate 1D velocity gradient for microseismic monitoring. The proposed algorithm is based on the analytically-derived ray formulas for a linearly increasing velocity model and the Gauss–Newton method. Numerical examples show that the proposed algorithm is robust to picking errors and initial guess.
We propose an elastic full-waveform inversion method based on P- and S-wave mode separation to mitigate the crosstalk artefacts between P- and S-wave modes by deriving the gradient formulas with respect to various wave modes using a P- and S-wave mode separated first-order velocity-stress wave equation.
Considering that the focal transform can flexibly perform order reduction or elevation of multiples and surface-related multiple elimination (SRME) can effectively separate primaries and multiples, a practical and novel method based on both the focal transform and SRME is proposed to separate different order multiples in this paper.
EG17085Three-dimensional gravity inversion using automatic regional-residual components separation with application to the Iwaki area in Fukushima Prefecture, Japan
We propose a simple method for analysing gravity data to define 3D subsurface structures based on automatic separation of gravity anomalies into residual gravity and trend surfaces. This is achieved by assuming that the gravity anomaly trend surface is extremely smooth prior to analysis. The new method was applied to an area around Iwaki City, Fukushima Prefecture, Japan, and the results evaluated using a boot-strap simulation.
A hydrogeophysically and geologically complex bitumen deposit in western Iran was investigated by means of resistivity and induced polarisation (IP) methods. The results show the high potential of IP data to improve bitumen detection. It was also proven that saline water can sometimes be useful in target identification in mineral explorations.
The Gilmore Fault Zone marks a distinct geophysical and complex geological boundary. This work improves on existing models of this fault zone and presents a revised potential field model. These models of the subsurface morphology of this major thrust and/or strike-slip fault system are being used to refine tectonic models for the evolution of the Eastern Lachlan Orogen.
We derive two formulas of forward modelling for least-squares reverse-time migration based on Born and Kirchhoff approximations. Analysis unveils an iω difference exists between the two formulas. Consequently, their seismic images have different shape/phase: the Born approximation produces anti-symmetric images while the Kirchhoff approximation gives symmetric images. Numerical examples demonstrate these features between the two formulas.
We present a compensation model of aircraft effects for full magnetic gradient tensor data acquired in a towed bird. Training flight tests were designed so that aircraft compensation parameters could be estimated. The feasibility of the compensation method was verified by modelling and by a flight simulation test.
EG17016Topographic correction of magnetic data on rugged topography with application to Río Blanco-Los Bronces and El Teniente porphyry copper districts, Southern Andes, Chile
This paper describes a semiquantitative, 3D-modelling based methodology that has been developed to remove topographic effects on magnetic data. Through synthetic and real cases in the Andes, we show that airborne magnetic surveys over rugged topography generate magnetic signatures that could be of similar amplitude with the signal of geological interest.
EG17012Hybrid absorbing boundary condition for piecewise smooth curved boundary in 2D acoustic finite difference modelling
Flexible computational domains with irregular absorbing boundaries can improve the efficiency of seismic modelling. We propose a hybrid absorbing boundary condition for piecewise smooth curved boundary in mesh-free discretisation based on radial-basis-function-generated finite difference modelling. Modelling examples in both homogeneous and heterogeneous media with flexibly shaped computational domain demonstrate the effectiveness of our proposed method.
We propose a modification of the reverse time migration imaging condition by reducing the modelling time to about half of the recording time. Results showed that the proposed modification has no significant effect on the illumination of reflectors, while increasing the CPU and memory allocation performances by about 25% and 50%, respectively.
A fast marching eikonal solver is employed in the isotropic and VTI concepts. Instead of using a linear eikonal equation, we consider a nonlinear approximation because it is more realistic and accurate than the former. Overall, anisotropic approach is stable, fast and generates high-quality images with accurate details in deep structures.
In this paper, we propose an adaptive denoising method based on data-driven signal mode decomposition, where the noise is represented by the residual/last mode. The proposed approach adaptively extracts the noise component depending on the data statistics rather than defining a fixed priori threshold.
EG17081Imaging the crustal structure along the KCRT-2008 profile in the Southern Korean Peninsula using regularisation
In this study, high-resolution subsurface structures of the southern Korean Peninsula were delineated. To acquire images, we applied seismic data processing tools that are commonly used in oil exploration to irregularly sampled data after regularisation by matching pursuit interpolation.
EG16072Enhancement of fault interpretation using multi-attribute analysis and artificial neural network (ANN) approach: a case study from Taranaki Basin, New Zealand
This paper describes an improved and efficient workflow that refines the art of structural interpretation from 3D seismic data. The research aims to streamline the interpretation workflow using better data conditioning scheme and defining a meta-attribute (fault probability cube) using an artificial neural network. A robust pathway is provided for delineating structural details of complex geological terrain.
Global lightning networks catalogue the time and location of up to four million lightning strikes per day. We use lightning network data to predict sferic arrival times, vertical electric field amplitudes and arrival azimuths at an audio-frequency magnetotelluric sounding near Heathcote, Victoria.
EG17007Application of full waveform inversion algorithms to seismic data lacking low-frequency information from a simple starting model
The low-frequency information and correct starting velocity are important for full waveform inversion (FWI). However, obtaining low-frequency information and accurate starting velocity from the field seismic exploration is difficult. This paper suggests a 4-phase FWI to invert the correct velocity model when a dataset lacks low-frequency information and accurate starting velocity.
EG17054Estimating high hydraulic conductivity locations through a 3D simulation of water flow in soil and a resistivity survey
We propose a simple method for estimating high hydraulic conductivity locations. The proposed method uses the 3D simulations of soil water flow and resistivity survey during a groundwater recharge experiment. Results of numerical and field experiments indicate that the proposed method estimates the high hydraulic conductivity locations more precisely compared with 3D inversion of in-line data.
A reduced order modelling method is introduced for wave propagation modelling by using the mode shapes of the model. The numerical accuracy, computational performance and boundary conditions of the proposed method were investigated. According to the results, the proposed method substantially improves the computational efficiency of the reverse time migration.
EG16115Regional and residual gravity anomaly separation using the singular spectrum analysis-based low pass filtering: a case study from Nagpur, Maharashtra, India
We present a singular spectrum analysis (SSA)-based algorithm for regional and residual gravity anomaly separation. Initially, we tested the method on synthetic data and applied gravity profiles from the Umred coalfield, Nagpur, Maharashtra, India. Final 2D crustal models obtained from the filtered regional anomaly agree well with regional geology and borehole information.
Two consistent methods of comparing borehole induction logs with models from inversion of AEM data have been developed: one in model space and one in data space. The two methods are applied to the Broken Hill Managed Aquifer Recharge project conducted by Geoscience Australia on AEM data from the SkyTEM system.
EG16048Study of Scholte wave dispersion curves and modal energy distribution using a wavefield numerical simulation method
In this paper, we examine the characteristics of Scholte wave dispersion curves and their modal energy distribution by combining analytical solutions of the dispersion equation with numerical modelling. The modelling results showed that our approach can facilitate Scholte wave exploration in water-covered areas.
EG160583D seismic attributes for structural mapping and enhancement of deep gold mining: a case study from the West Wits Line goldfields, South Africa
In this study, we use sophisticated seismic attributes to interpret the 3D reflection seismic data to: (1) assess and mitigate the risks posed by deep mining activities and (2) improve the resource evaluation of the gold-bearing quartz pebble conglomerate horizons (reefs) in the world’s deepest gold mines (South Africa).
In this paper, we present a scheme to incorporate 3D controlled-source audio frequency magnetotelluric (CSAMT) topographic distortions into the 3D inversion instead of correcting them. This approach has been verified by comparison with 2D FEM CSAMT solutions and synthetic inversion examples. The field example also illustrates the effectiveness of our approach.
EG16079Three-dimensional tensor controlled-source audio-frequency magnetotelluric inversion using LBFGS
We have shown in this paper that directly using a magnetotelluric method to invert the data of the tensor controlled-source audio-frequency magnetotelluric (CSAMT) will obtain an incorrect result, the inversion result of tensor CSAMT is more reliable than that of the traditional CSAMT, and the limited-memory Broyden-Fletcher-Goldfarb-Shanno (LBFGS) method is more efficient than the nonlinear conjugate gradient (NLCG) method for tensor CSAMT.
EG16097Azimuthal seismic responses from shale formation based on anisotropic rock physics and reflectivity method: a case study from south-west China
We have constructed an anisotropic rock physics model at the seismic scale for shales, and applied it to the Longmaxi Shale in the Sichuan Basin, south-west China. The simplified reflectivity method is proposed to calculate seismic responses for amplitude variation with offset and azimuth (AVAz) analysis.
EG16038Mapping electrical structures in the southern Great Khingan Range, north-east China, through two-dimensional magnetotelluric sounding
Magnetotelluric measurements were carried out in the southern Great Khingan Range, Inner Mongolia, China, from 2012 to 2014. After 2D inversion, the Paleozoic Hunnitu and Gadasu depressions were discovered and inferred to possibly contain Linxi strata. The Mesozoic Xiretu depression was also discovered, which possibly contains Cretaceous hydrocarbon-bearing strata.
EG16100Regularisation parameter adaptive selection and its application in the prestack AVO inversion
In this paper, based on L-curve criterion, we propose an improved method for the adaptive acquisition of regularisation parameters for arbitrary norm condition. A detailed derivation of the proposed method is described. Numerical experiments confirm that the proposed method is more accurate and robust than its main competitor, generalised cross-validation.
We present a novel approach for surface-consistent residual phase corrections based on migrated gathers to improve the SNR and resolution of migrated images. This method includes the following new aspects: the residual statics and phase are estimated by dual parameter cross-correlation, and determined after migration.
We have conducted numerical and field experiments to investigate the applicability of electrode configurations and line layouts commonly used for 2D resistivity surveys to 3D inversion. We propose that parallel lines are useful to highlight areas of particular interest where further detailed work with an intersecting line could be carried out.
EG16046Noise reduction of grounded electrical source airborne transient electromagnetic data using an exponential fitting-adaptive Kalman filter
GREATEM field data usually includes a mixed variety of noises, which makes the exponential decaying signal too difficult to identify. This paper presents an exponential fitting-adaptive Kalman filter (EF-AKF) to remove mixed electromagnetic noises, while preserving the signal characteristics. As a new method, the EF-AKF can be used for denoising exponential decaying signals.
EG16022Subsurface attenuation estimation using a novel hybrid method based on FWE function and power spectrum
A new Q estimation approach, called improved frequency weighted exponential (IFWE), is presented by combining the advantage of the FWE method and the power spectrum. Tests of synthetic and field data show that the IFWE is more robust and the bandwidth selection for the IFWE is more tolerant than the FWE.
EG16139Trans-dimensional Bayesian inversion of airborne electromagnetic data for 2D conductivity profiles
We apply a novel trans-dimensional Bayesian approach using a wavelet parameterisation to airborne electromagnetic (AEM) inversions using data from the Broken Hill region. This approach allows exploration of a range of plausible subsurface conductivity models and provides more robust uncertainty estimates while accounting for potential non-uniqueness.
EG16094An implicit spatial and high-order temporal finite difference scheme for 2D acoustic modelling
We propose a new finite difference (FD) scheme that is implicit in space (adopting the Taylor-expansion-based or least-square based FD coefficients) and has high-order accuracy in time (adopting a combination of rhombus and cross stencils). Dispersion analysis, stability analysis and modelling examples validate the superiority of our method.
To evaluate the crustal structure of a site proposed for ocean drilling, high-resolution dense seismic surveys were carried out in the Izu rear arc, Japan. Our results suggest that use of a dense grid configuration is important in volcanic rear-arc situations, in order to recognise heterogeneous crustal structure, such as sediments from different origins.
The artefacts caused by random fluctuations in the background velocity model will contaminate the final image, which makes the work of interpretation difficult to implement in reverse-time migration. The imaging condition proposed in this paper, which exploits the pseudo-Wigner distribution function, can reduce this noise, and also help improve the continuity of complex layers.
To alleviate the effect of phase in waveform classification, the adaptive phase k-means is introduced for unsupervised seismic facies analysis. This method improves the traditional k-means algorithm by using an adaptive phase distance for waveform similarity measure, and is thus robust to phase variations caused by horizon interpretation.
EG16057Investigating subsurface structures of Gachsaran oil field in Iran using 2D inversion of magnetotelluric data
In 2012, the National Iranian Oil Company conducted an electromagnetic survey in the Gachsaran oil field. The Bahr’s skew and Mohr diagrams were used to perform the dimensionality analysis. Magnetotelluric data were modelled using the smoothness-constrained least-squares method. The resulting model revealed the main anticline and overthrust zone in the region.
EG16011Effect of ground undulation and mounted vehicle velocity variation on stepped frequency continuous wave GPR data
This paper describes a method to generate dataset on stepped frequency continuous wave ground penetrating radar for land mine detection. A novel ground undulation and ground bounce removal filter has been proposed which mitigates the effects of vehicle velocity variations. Kirchhoff and F-K migration algorithms are applied to the filtered output.
Residual noise remains in airborne time domain electromagnetic profiles after preprocessing the data, and this noise affects the exploration of targets. An approach to reduce this noise based on the minimum noise fraction has been proposed in this paper.
EG16136Comparative analyses of seismic site conditions and microzonation of the major cities in Gangwon Province, Korea
To determine the seismic site conditions and microzonation of Chuncheon, Wonju, and Gangneung in Korea, shear-wave velocities were derived at 313 sites by the ESPAC method. The proxy-based Vs30 indicated that these cities were mainly categorised into NEHRP classes B, C, and D, with a minor proportion of A.
EG16044MOPSO: a new computing algorithm for joint inversion of Rayleigh wave dispersion curve and refraction traveltimes
This study introduces the multi-objective particle swarm optimisation (MOPSO) strategy as a new tool for joint inversion of Rayleigh wave dispersion curve and refraction traveltimes. The proposed algorithm was tested on two synthetic datasets and an experimental dataset. The results showed that the applied joint inversion technique can considerably reduce uncertainties of the inverted models.
EG16037Multi-scale full waveform inversion for areas with irregular surface topography in an auxiliary coordinate system
We implement full waveform inversion in an auxiliary coordinate system to improve inversion quality of near-surface regions with strong elevation and velocity variation. Furthermore, a time-domain multi-scale decomposition method and an optimised encoding strategy are introduced to the inversion frame to promote the practical application of our method.
EG16056Applying multi-wavelength, multi-polarimetric and interferometric airborne radar data for geological mapping in densely forested north-western Tasmania
Airborne multi-wavelength, multi-polarimetric and interferometric radar data (NASA’s AIRSAR/TOPSAR) was processed and applied for geological mapping within densely forested north-western Tasmania. The results successfully compared with geophysical, topographic and vegetation classification information and showed the potential of such data to assist in the upgrading of existing published geological mapping.
EG16042Regional magnetic lithologies and structures as controls on porphyry copper deposits: evidence from Iran
The Tertiary magmatic belts host the porphyry copper deposits in Iran. The combination of upward continuation and edge enhancement filters were used to find and predict the pattern of distribution of these deposits. The results prove that the deposits locate on regional magnetic lithologies and structures.
EG15114Seismic data decomposition and reconstruction with sparse Gaussian beams and sparse optimisation method
Sparse Gaussian beam decomposition and reconstruction is an ill-posed inverse problem. We construct a zero-norm constrained minimisation model, develop a dip-angle scanning strategy and adopt a quasi-Newton method to find the solution. This method can recover full wave field data and can be used for sparse Gaussian beam migration.
We have developed a new interpolation process for insufficient 3D field data that applies 2D curvelet transform-based projection onto convex sets (POCS) to the kx-ky transformed data of each time slice of 3D data. Additionally, to acquire accurate interpolation results, we have designed the preparation process to render the input data irregularly distributed with small-sized gaps.
EG16076New constraints on the neotectonic stress pattern of the Flinders and Mount Lofty Ranges, South Australia
The results of this study suggest a regional E–W orientation of the maximum horizontal stress in the Flinders and Mount Lofty Ranges. Focal mechanism solutions of earthquakes suggest that the majority of events occur in a thrust faulting stress regime. However, our data also indicate the presence of strike-slip and normal faulting stress regimes in the region, which has not been suggested extensively before this study.
EG15123Traveltimes and conversion-point positions for P-SV converted wave propagation in a transversely isotropic medium: numerical calculations and physical model studies
This study uses ultrasonic physical modelling to verify that Fermat’s minimum-time principle is better than the anisotropic non-hyperbolic moveout and conversion-point (CP) equations for calculating the traveltime and CP position of a P-SV wave reflected from a strong vertical transversely isotropic medium.
The outliers-out (OlO) stack is a novel seismic data stacking technique that is based on excluding a specific number of outliers from each time sample before stacking. The OlO stacking parameters are automatically obtained from the data according to the statistical distribution of amplitudes. The superiority of the OlO stack has been verified by experiments.
These articles have been peer reviewed and accepted for publication. They are still in production and have not been edited, so may differ from the final published form.
Comparison of the Projection Onto Convex Sets and Iterative Hard Thresholding methods for seismic data interpolation and denoising
The Most Read ranking is based on the number of downloads in the last 60 days from papers published on the CSIRO PUBLISHING website within the last 12 months. Usage statistics are updated daily.
Exploration Geophysics (Online Early)Philip Heath
Exploration Geophysics (Online Early)Frederick E. M. Lilley
Seismic facies-controlled prestack simultaneous inversion of elastic and petrophysical parameters for favourable reservoir predictionExploration Geophysics (Online Early)Sheng Zhang, Handong Huang, Baoheng Zhu, Huijie Li, Lihua Zhang
Velocity analysis of CMP gathers acquired by an array GPR system ‘Yakumo’: results from field application to tsunami depositsExploration Geophysics (Online Early)Honghua Wang, Hai Liu, Jie Cui, Xiangyun Hu, Motoyuki Sato
Robust data processing of noisy marine controlled-source electromagnetic data using independent component analysisExploration Geophysics (Online Early)Naoto Imamura, Tada-nori Goto, Takafumi Kasaya, Hideaki Machiyama
Synthetic modelling and analysis of CSEM full-field apparent resistivity response combining EM induction and IP effect for 1D mediumExploration Geophysics (Online Early)Weiqiang Liu, Pinrong Lin, Qingtian Lü, Yong Li, Jianhua Li
A new method to enhance the characterisation of seismically thin beds based on the generalised S transform maximum modulusExploration Geophysics (Online Early)Mingfei Fan, Shenghe Wu, Guangyi Hu, Jingjing Qu, Jiajia Zhang
Exploration Geophysics (Online Early)Fangyu Li, Bo Zhang, Sumit Verma, Kurt J. Marfurt
Amplification characteristics at Iedang Reservoir dam sites determined using H/V spectral ratio with background noise, S-wave and coda wave energyExploration Geophysics (Online Early)Jun Kyoung Kim, Soung Hoon Wee, Seong Hwa Yoo, Kwang Hee Kim
Application of full waveform inversion algorithms to seismic data lacking low-frequency information from a simple starting modelExploration Geophysics (Online Early)Hyunggu Jun, Jungkyun Shin, Changsoo Shin
Exploration Geophysics (Online Early)Rhys Hawkins, Ross C. Brodie, Malcolm Sambridge
Exploration Geophysics (Online Early)Deepika Venkataramani, Robert J. Musgrave, David A. Boutelier, Alistair C. Hack, William J. Collins
Exploration Geophysics (Online Early)Changhong Lin, Handong Tan, Wangwang Wang, Tuo Tong, Miao Peng, Mao Wang, Weihua Zeng
3D seismic attributes for structural mapping and enhancement of deep gold mining: a case study from the West Wits Line goldfields, South AfricaExploration Geophysics (Online Early)Nomqhele Z. Nkosi, Musa S. D. Manzi, Oleg Brovko, Raymond J. Durrheim
Mapping electrical structures in the southern Great Khingan Range, north-east China, through two-dimensional magnetotelluric soundingExploration Geophysics (Online Early)Weijun Zhao, Shuwang Chen, Qiuhong Ding
Exploration Geophysics (Online Early)Chunxiao Xiu, Xiaohong Meng, Lianghui Guo, Sheng Zhang, Xingdong Zhang
Exploration Geophysics (Online Early)Kunpeng Wang, Handong Tan, Changhong Lin, Jianlong Yuan, Cong Wang, Jing Tang
Enhancement of fault interpretation using multi-attribute analysis and artificial neural network (ANN) approach: a case study from Taranaki Basin, New ZealandExploration Geophysics (Online Early)Priyadarshi Chinmoy Kumar, Animesh Mandal
Subsurface attenuation estimation using a novel hybrid method based on FWE function and power spectrumExploration Geophysics (Online Early)Jingnan Li, Shangxu Wang, Dengfeng Yang, Genyang Tang, Yangkang Chen
Study of Scholte wave dispersion curves and modal energy distribution using a wavefield numerical simulation methodExploration Geophysics (Online Early)Hua Wu, Guang-Zhou Shao, Qing-Chun Li