scholarly journals Effects of full transmitting-current waveforms on transient electromagnetics: Insights from modeling the Albany graphite deposit

Geophysics ◽  
2019 ◽  
Vol 84 (4) ◽  
pp. E255-E268 ◽  
Author(s):  
Sihong Zeng ◽  
Xiangyun Hu ◽  
Jianhui Li ◽  
Colin G. Farquharson ◽  
Peter C. Wood ◽  
...  
Keyword(s):  
Real Data ◽  
Forward Modeling ◽  
Northern Ontario ◽  
Data Set ◽  
Rock Types ◽  
Full Waveform ◽  
Transient Electromagnetic ◽  
Steady Stage ◽  
Geologic Model ◽  
Off Time

In transient electromagnetic (TEM) methods, the full transmitting-current waveform, not just the abrupt turn-off, can have effects on the measured responses. A 3D finite-element time-domain forward-modeling solver was used to investigate these effects. This was motivated by an attempt to match, via forward-modeling, real data from the Albany graphite deposit in northern Ontario, Canada. Initial modeling results for homogeneous half-spaces illustrate the effects that a full waveform can have on TEM responses, especially the durations of the steady stage and turn-off time. For the Albany data set, a geophysical conductivity model was developed from a geologic model that itself had been constructed predominantly from drillhole information. The conductivities of the various geologic units in the model were first estimated based on typical conductivity values for the respective rock types, then adjusted to fit the measured TEM data as closely as possible. We found that the TEM responses differed significantly from the pure step-off response and that incorporating the effects of the full waveform (particularly the linear ramp turn-off) greatly improved the match between observed and computed responses, especially for the early measurement times. In addition, this Albany example illustrates the presence of sign changes in TEM data caused primarily by localized conductivity targets.

A new parameter set for anisotropic multiparameter full-waveform inversion and application to a North Sea data set

Geophysics ◽  
2016 ◽  
Vol 81 (4) ◽  
pp. U25-U38 ◽  
Author(s):  
Nuno V. da Silva ◽  
Andrew Ratcliffe ◽  
Vetle Vinje ◽  
Graham Conroy
Keyword(s):  
North Sea ◽  
Waveform Inversion ◽  
Real Data ◽  
Model Parameters ◽  
Full Waveform Inversion ◽  
Radiation Patterns ◽  
Data Set ◽  
Full Waveform ◽  
Seismic Image ◽  
Central North Sea

Parameterization lies at the center of anisotropic full-waveform inversion (FWI) with multiparameter updates. This is because FWI aims to update the long and short wavelengths of the perturbations. Thus, it is important that the parameterization accommodates this. Recently, there has been an intensive effort to determine the optimal parameterization, centering the fundamental discussion mainly on the analysis of radiation patterns for each one of these parameterizations, and aiming to determine which is best suited for multiparameter inversion. We have developed a new parameterization in the scope of FWI, based on the concept of kinematically equivalent media, as originally proposed in other areas of seismic data analysis. Our analysis is also based on radiation patterns, as well as the relation between the perturbation of this set of parameters and perturbation in traveltime. The radiation pattern reveals that this parameterization combines some of the characteristics of parameterizations with one velocity and two Thomsen’s parameters and parameterizations using two velocities and one Thomsen’s parameter. The study of perturbation of traveltime with perturbation of model parameters shows that the new parameterization is less ambiguous when relating these quantities in comparison with other more commonly used parameterizations. We have concluded that our new parameterization is well-suited for inverting diving waves, which are of paramount importance to carry out practical FWI successfully. We have demonstrated that the new parameterization produces good inversion results with synthetic and real data examples. In the latter case of the real data example from the Central North Sea, the inverted models show good agreement with the geologic structures, leading to an improvement of the seismic image and flatness of the common image gathers.

Frequency‐domain acoustic‐wave modeling and inversion of crosshole data: Part II—Inversion method, synthetic experiments and real‐data results

Geophysics ◽  
1995 ◽  
Vol 60 (3) ◽  
pp. 796-809 ◽  
Author(s):  
Zhong‐Min Song ◽  
Paul R. Williamson ◽  
R. Gerhard Pratt
Keyword(s):  
Frequency Domain ◽  
Real Data ◽  
Forward Modeling ◽  
Two Dimensions ◽  
Inversion Method ◽  
Inversion Algorithm ◽  
Data Set ◽  
Wave Inversion ◽  
Full Wave ◽  
And Inversion

In full‐wave inversion of seismic data in complex media it is desirable to use finite differences or finite elements for the forward modeling, but such methods are still prohibitively expensive when implemented in 3-D. Full‐wave 2-D inversion schemes are of limited utility even in 2-D media because they do not model 3-D dynamics correctly. Many seismic experiments effectively assume that the geology varies in two dimensions only but generate 3-D (point source) wavefields; that is, they are “two‐and‐one‐half‐dimensional” (2.5-D), and this configuration can be exploited to model 3-D propagation efficiently in such media. We propose a frequency domain full‐wave inversion algorithm which uses a 2.5-D finite difference forward modeling method. The calculated seismogram can be compared directly with real data, which allows the inversion to be iterated. We use a descents‐related method to minimize a least‐squares measure of the wavefield mismatch at the receivers. The acute nonlinearity caused by phase‐wrapping, which corresponds to time‐domain cycle‐skipping, is avoided by the strategy of either starting the inversion using a low frequency component of the data or constructing a starting model using traveltime tomography. The inversion proceeds by stages at successively higher frequencies across the observed bandwidth. The frequency domain is particularly efficient for crosshole configurations and also allows easy incorporation of attenuation, via complex velocities, in both forward modeling and inversion. This also requires the introduction of complex source amplitudes into the inversion as additional unknowns. Synthetic studies show that the iterative scheme enables us to achieve the theoretical maximum resolution for the velocity reconstruction and that strongly attenuative zones can be recovered with reasonable accuracy. Preliminary results from the application of the method to a real data set are also encouraging.

scholarly journals Accelerated full-waveform inversion using dynamic mini-batches

2020 ◽  
Vol 221 (2) ◽  
pp. 1427-1438 ◽  
Author(s):  
Dirk Philip van Herwaarden ◽  
Christian Boehm ◽  
Michael Afanasiev ◽  
Solvi Thrastarson ◽  
Lion Krischer ◽  
...  
Keyword(s):  
Waveform Inversion ◽  
Real Data ◽  
Complete Data ◽  
Control Group ◽  
Full Waveform Inversion ◽  
African Plate ◽  
Data Set ◽  
Full Waveform ◽  
Spatial Coverage ◽  
Quasi Newton

SUMMARY We present an accelerated full-waveform inversion based on dynamic mini-batch optimization, which naturally exploits redundancies in observed data from different sources. The method rests on the selection of quasi-random subsets (mini-batches) of sources, used to approximate the misfit and the gradient of the complete data set. The size of the mini-batch is dynamically controlled by the desired quality of the gradient approximation. Within each mini-batch, redundancy is minimized by selecting sources with the largest angular differences between their respective gradients, and spatial coverage is maximized by selecting candidate events with Mitchell’s best-candidate algorithm. Information from sources not included in a specific mini-batch is incorporated into each gradient calculation through a quasi-Newton approximation of the Hessian, and a consistent misfit measure is achieved through the inclusion of a control group of sources. By design, the dynamic mini-batch approach has several main advantages: (1) The use of mini-batches with adaptive size ensures that an optimally small number of sources is used in each iteration, thus potentially leading to significant computational savings; (2) curvature information is accumulated and exploited during the inversion, using a randomized quasi-Newton method; (3) new data can be incorporated without the need to re-invert the complete data set, thereby enabling an evolutionary mode of full-waveform inversion. We illustrate our method using synthetic and real-data inversions for upper-mantle structure beneath the African Plate. In these specific examples, the dynamic mini-batch approach requires around 20 per cent of the computational resources in order to achieve data and model misfits that are comparable to those achieved by a standard full-waveform inversion where all sources are used in each iteration.

scholarly journals Inverting magnetotelluric data with distortion correction—stability, uniqueness and trade-off with model structure

2020 ◽  
Vol 222 (3) ◽  
pp. 1620-1638 ◽  
Author(s):  
M Moorkamp ◽  
A Avdeeva ◽  
Ahmet T Basokur ◽  
Erhan Erdogan
Keyword(s):  
Real Data ◽  
Distortion Correction ◽  
Impedance Tensor ◽  
Strong Impact ◽  
Magnetotelluric Data ◽  
Data Set ◽  
Near Surface ◽  
Transient Electromagnetic ◽  
Electromagnetic Measurements ◽  
Galvanic Distortion

SUMMARY Galvanic distortion of magnetotelluric (MT) data is a common effect that can impede the reliable imaging of subsurface structures. Recently, we presented an inversion approach that includes a mathematical description of the effect of galvanic distortion as inversion parameters and demonstrated its efficiency with real data. We now systematically investigate the stability of this inversion approach with respect to different inversion strategies, starting models and model parametrizations. We utilize a data set of 310 MT sites that were acquired for geothermal exploration. In addition to impedance tensor estimates over a broad frequency range, the data set also comprises transient electromagnetic measurements to determine near surface conductivity and estimates of distortion at each site. We therefore can compare our inversion approach to these distortion estimates and the resulting inversion models. Our experiments show that inversion with distortion correction produces stable results for various inversion strategies and for different starting models. Compared to inversions without distortion correction, we can reproduce the observed data better and reduce subsurface artefacts. In contrast, shifting the impedance curves at high frequencies to match the transient electromagnetic measurements reduces the misfit of the starting model, but does not have a strong impact on the final results. Thus our results suggest that including a description of distortion in the inversion is more efficient and should become a standard approach for MT inversion.

scholarly journals Source-independent efficient wavefield inversion

2020 ◽  
Vol 222 (1) ◽  
pp. 697-714
Author(s):  
Chao Song ◽  
Tariq Alkhalifah
Keyword(s):  
Waveform Inversion ◽  
Computational Cost ◽  
Real Data ◽  
Velocity Model ◽  
Data Set ◽  
Full Waveform ◽  
Source Estimation ◽  
Reference Trace ◽  
New Formulation ◽  
Wavefield Inversion

SUMMARY Full-waveform inversion (FWI) is an effective tool to retrieve a high-resolution subsurface velocity model. The source wavelet accuracy plays an important role in reaching that goal. So we often need to estimate the source function before or within the inversion process. Source estimation requires additional computational cost, and an inaccurate source estimation can hamper the convergence of FWI. We develop a source-independent waveform inversion utilizing a recently introduced wavefield reconstruction based method, which we refer to as efficient wavefield inversion (EWI). In EWI, we essentially reconstruct the wavefield by fitting it to the observed data as well as a wave equation based on iterative Born scattering. However, a wrong source wavelet will induce errors in the reconstructed wavefield, which may lead to a divergence of this optimization problem. We use a convolution-based source-independent misfit function to replace the conventional data fitting term in EWI to formulate a source-independent EWI (SIEWI) objective function. By convolving the observed data with a reference trace from the predicted data and convolving the predicted data with a reference trace from the observed data, the influence of the source wavelet on the optimization is mitigated. In SIEWI, this new formulation is able to mitigate the cycle-skipping issue and the source wavelet uncertainty simultaneously. We demonstrate those features on the Overthrust model and a modified Marmousi model. Application on a 2-D real data set also shows the effectiveness of the proposed method.

EXPONENTIATED HALF-LOGISTIC LOMAX DISTRIBUTION WITH PROPERTIES AND APPLICATION

2019 ◽  
Vol XVI (2) ◽  
pp. 1-11
Author(s):  
Farrukh Jamal ◽  
Hesham Mohammed Reyad ◽  
Soha Othman Ahmed ◽  
Muhammad Akbar Ali Shah ◽  
Emrah Altun
Keyword(s):  
Real Data ◽  
Continuous Model ◽  
Model Parameters ◽  
Data Set ◽  
Lomax Distribution ◽  
Mathematical Properties ◽  
Proposed Model ◽  
Probability Weighted Moment ◽  
Record Statistics ◽  
Maximum Likelihood Criterion

A new three-parameter continuous model called the exponentiated half-logistic Lomax distribution is introduced in this paper. Basic mathematical properties for the proposed model were investigated which include raw and incomplete moments, skewness, kurtosis, generating functions, Rényi entropy, Lorenz, Bonferroni and Zenga curves, probability weighted moment, stress strength model, order statistics, and record statistics. The model parameters were estimated by using the maximum likelihood criterion and the behaviours of these estimates were examined by conducting a simulation study. The applicability of the new model is illustrated by applying it on a real data set.

Evaluation for estimating of the PDF and the CDF of Generalized Inverted Exponential Distribution with Application in Industry

2020 ◽  
pp. 507-522
Author(s):  
Parisa Torkaman
Keyword(s):  
Least Squares ◽  
Exponential Distribution ◽  
Mean Squared Error ◽  
Weighted Least Squares ◽  
Real Data ◽  
Minimum Variance ◽  
Cumulative Distribution ◽  
Estimation Methods ◽  
Data Set ◽  
Better Than

The generalized inverted exponential distribution is introduced as a lifetime model with good statistical properties. This paper, the estimation of the probability density function and the cumulative distribution function of with five different estimation methods: uniformly minimum variance unbiased(UMVU), maximum likelihood(ML), least squares(LS), weighted least squares (WLS) and percentile(PC) estimators are considered. The performance of these estimation procedures, based on the mean squared error (MSE) by numerical simulations are compared. Simulation studies express that the UMVU estimator performs better than others and when the sample size is large enough the ML and UMVU estimators are almost equivalent and efficient than LS, WLS and PC. Finally, the result using a real data set are analyzed.

HCVS: Pinpointing Chromatin States Through Hierarchical Clustering and Visualization Scheme

2019 ◽  
Vol 14 (2) ◽  
pp. 148-156
Author(s):  
Nighat Noureen ◽  
Sahar Fazal ◽  
Muhammad Abdul Qadir ◽  
Muhammad Tanvir Afzal
Keyword(s):  
Hierarchical Clustering ◽  
Real Data ◽  
Cell Types ◽  
Computational Scheme ◽  
Data Set ◽  
Chromatin States ◽  
Functional Regions ◽  
Visualization Strategy ◽  
Hidden States ◽  
Next Generation Sequencing Ngs

Background: Specific combinations of Histone Modifications (HMs) contributing towards histone code hypothesis lead to various biological functions. HMs combinations have been utilized by various studies to divide the genome into different regions. These study regions have been classified as chromatin states. Mostly Hidden Markov Model (HMM) based techniques have been utilized for this purpose. In case of chromatin studies, data from Next Generation Sequencing (NGS) platforms is being used. Chromatin states based on histone modification combinatorics are annotated by mapping them to functional regions of the genome. The number of states being predicted so far by the HMM tools have been justified biologically till now. Objective: The present study aimed at providing a computational scheme to identify the underlying hidden states in the data under consideration. </P><P> Methods: We proposed a computational scheme HCVS based on hierarchical clustering and visualization strategy in order to achieve the objective of study. Results: We tested our proposed scheme on a real data set of nine cell types comprising of nine chromatin marks. The approach successfully identified the state numbers for various possibilities. The results have been compared with one of the existing models as well which showed quite good correlation. Conclusion: The HCVS model not only helps in deciding the optimal state numbers for a particular data but it also justifies the results biologically thereby correlating the computational and biological aspects.

Sign in / Sign up

Export Citation Format

Share Document