Large-scale Gauss-Newton inversion of transient controlled-source electromagnetic measurement data using the model reduction framework

Geophysics ◽  
2013 ◽  
Vol 78 (4) ◽  
pp. E161-E171 ◽  
Author(s):  
M. Zaslavsky ◽  
V. Druskin ◽  
A. Abubakar ◽  
T. Habashy ◽  
V. Simoncini
Keyword(s):  
Frequency Domain ◽  
Time Domain ◽  
Large Scale ◽  
Jacobian Matrix ◽  
Computational Cost ◽  
Measurement Data ◽  
Multiple Sources ◽  
Time Data ◽  
Controlled Source ◽  
The Time Domain

Transient data controlled-source electromagnetic measurements are usually interpreted via extracting few frequencies and solving the corresponding inverse frequency-domain problem. Coarse frequency sampling may result in loss of information and affect the quality of interpretation; however, refined sampling increases computational cost. Fitting data directly in the time domain has similar drawbacks, i.e., its large computational cost, in particular, when the Gauss-Newton (GN) algorithm is used for the misfit minimization. That cost is mainly comprised of the multiple solutions of the forward problem and linear algebraic operations using the Jacobian matrix for calculating the GN step. For large-scale 2.5D and 3D problems with multiple sources and receivers, the corresponding cost grows enormously for inversion algorithms using conventional finite-difference time-domain (FDTD) algorithms. A fast 3D forward solver based on the rational Krylov subspace (RKS) reduction algorithm using an optimal subspace selection was proposed earlier to partially mitigate this problem. We applied the same approach to reduce the size of the time-domain Jacobian matrix. The reduced-order model (ROM) is obtained by projecting a discretized large-scale Maxwell system onto an RKS with optimized poles. The RKS expansion replaces the time discretization for forward and inverse problems; however, for the same or better accuracy, its subspace dimension is much smaller than the number of time steps of the conventional FDTD. The crucial new development of this work is the space-time data compression of the ROM forward operator and decomposition of the ROM’s time-domain Jacobian matrix via chain rule, as a product of time- and space-dependent terms, thus effectively decoupling the discretizations in the time and parameter spaces. The developed technique can be equivalently applied to finely sampled frequency-domain data. We tested our approach using synthetic 2.5D examples of hydrocarbon reservoirs in the marine environment.

Solution of large-scale 3D controlled-source electromagnetic modeling problem using efficient iterative solvers

Geophysics ◽  
2021 ◽  
pp. 1-64
Author(s):  
Changkai Qiu ◽  
Changchun Yin ◽  
Yunhe Liu ◽  
Xiuyan Ren ◽  
Hui Chen ◽  
...  
Keyword(s):  
Linear System ◽  
Linear Equation ◽  
Frequency Domain ◽  
Time Domain ◽  
Large Scale ◽  
Equation System ◽  
Iterative Solvers ◽  
Controlled Source ◽  
Linear Equation System ◽  
The Time Domain

With geophysical surveys evolving from traditional 2D to 3D models, the large volume of data adds challenges to inversion, especially when aiming to resolve complex 3D structures. An iterative forward solver for a controlled-source electromagnetic method (CSEM) requires less memory than that for a direct solver; however, it is not easy to iteratively solve an ill-conditioned linear system of equations arising from finite-element discretization of Maxwell’s equations. To solve this problem, we have developed efficient and robust iterative solvers for frequency- and time-domain CSEM modeling problems. For the frequency-domain problem, we first transform the linear system into its equivalent real-number format, and then introduce an optimal block-diagonal preconditioner. Because the condition number of the preconditioned linear equation system has an upper bound of √2, we can achieve fast solution convergence when applying a flexible generalized minimum residual solver. Applying the block preconditioner further results in solving two smaller linear systems with the same coefficient matrix. For the time-domain problem, we first discretize the partial differential equation for the electric field in time using an unconditionally stable backward Euler scheme. We then solve the resulting linear equation system iteratively at each time step. After the spatial discretization in the frequency domain, or space-time discretization in the time domain, we exploit the conjugate-gradient solver with auxiliary-space preconditioners derived from the Hiptmair-Xu decomposition to solve these real linear systems. Finally, we check the efficiency and effectiveness of our iterative methods by simulating complex CSEM models. The most significant advantage of our approach is that the iterative solvers that we adopt have almost the same accuracy and robustness as direct solvers but require much less memory, rendering them more suitable for large-scale 3D CSEM forward modeling and inversion.

Frequency‐domain elastic wave modeling by finite differences: A tool for crosshole seismic imaging

Geophysics ◽  
1990 ◽  
Vol 55 (5) ◽  
pp. 626-632 ◽  
Author(s):  
R. Gerhard Pratt
Keyword(s):  
Frequency Domain ◽  
Time Domain ◽  
Equations Of Motion ◽  
Domain Modeling ◽  
Multiple Sources ◽  
State Equations ◽  
Migration Imaging ◽  
Time Domain Modeling ◽  
Frequency Components ◽  
The Time Domain

The migration, imaging, or inversion of wide‐aperture cross‐hole data depends on the ability to model wave propagation in complex media for multiple source positions. Computational costs can be considerably reduced in frequency‐domain imaging by modeling the frequency‐domain steady‐state equations, rather than the time‐domain equations of motion. I develop a frequency‐domain approach in this note that is competitive with time‐domain modeling when solutions for multiple sources are required or when only a limited number of frequency components of the solution are required.

Frequency-domain wave-equation traveltime inversion with a monofrequency component

Geophysics ◽  
2021 ◽  
Vol 86 (6) ◽  
pp. R913-R926
Author(s):  
Jianhua Wang ◽  
Jizhong Yang ◽  
Liangguo Dong ◽  
Yuzhu Liu
Keyword(s):  
Wave Equation ◽  
Frequency Domain ◽  
Time Domain ◽  
Computational Cost ◽  
Velocity Model ◽  
Receiver Side ◽  
Data Set ◽  
Traveltime Inversion ◽  
Velocity Models ◽  
The Time Domain

Wave-equation traveltime inversion (WTI) is a useful tool for background velocity model building. It is generally formulated and implemented in the time domain, in which the gradient is calculated by temporally crosscorrelating the source- and receiver-side wavefields. The time-domain source-side snapshots are either stored in memory or are reconstructed through back propagation. The memory requirements and computational cost of WTI are thus prohibitively expensive, especially for 3D applications. To partially alleviate this problem, we provide an implementation of WTI in the frequency domain with a monofrequency component. Because only one frequency is used, it is affordable to directly store the source- and receiver-side wavefields in memory. There is no need for wavefield reconstruction during gradient calculation. In such a way, we have dramatically reduced the memory requirements and computational cost compared with the traditional time-domain WTI realization. For practical implementation, the frequency-domain wavefield is calculated by time-domain finite-difference forward modeling and is transformed to the frequency domain by an on-the-fly discrete Fourier transform. Numerical examples on a simple lateral periodic velocity model and the Marmousi model demonstrate that our method can obtain accurate background velocity models comparable with those from time-domain WTI and frequency-domain WTI with multiple frequencies. A field data set test indicates that our method obtains a background velocity model that well predicts the seismic wave traveltime.

JOINT INTERPRETATION OF AIRBORNE ELECTROMAGNETIC DATA IN THE TIME DOMAIN AND FREQUENCY DOMAIN

2018 ◽  
Vol 12 (7-8) ◽  
pp. 76-83
Author(s):  
E. V. KARSHAKOV ◽  
J. MOILANEN
Keyword(s):  
Fourier Transform ◽  
Frequency Domain ◽  
Time Domain ◽  
Frequency Interval ◽  
Single Spectrum ◽  
Simultaneous Inversion ◽  
Homogeneous Half Space ◽  
Time Domain Data ◽  
The Time Domain

Тhe advantage of combine processing of frequency domain and time domain data provided by the EQUATOR system is discussed. The heliborne complex has a towed transmitter, and, raised above it on the same cable a towed receiver. The excitation signal contains both pulsed and harmonic components. In fact, there are two independent transmitters operate in the system: one of them is a normal pulsed domain transmitter, with a half-sinusoidal pulse and a small "cut" on the falling edge, and the other one is a classical frequency domain transmitter at several specially selected frequencies. The received signal is first processed to a direct Fourier transform with high Q-factor detection at all significant frequencies. After that, in the spectral region, operations of converting the spectra of two sounding signals to a single spectrum of an ideal transmitter are performed. Than we do an inverse Fourier transform and return to the time domain. The detection of spectral components is done at a frequency band of several Hz, the receiver has the ability to perfectly suppress all sorts of extra-band noise. The detection bandwidth is several dozen times less the frequency interval between the harmonics, it turns out thatto achieve the same measurement quality of ground response without using out-of-band suppression you need several dozen times higher moment of airborne transmitting system. The data obtained from the model of a homogeneous half-space, a two-layered model, and a model of a horizontally layered medium is considered. A time-domain data makes it easier to detect a conductor in a relative insulator at greater depths. The data in the frequency domain gives more detailed information about subsurface. These conclusions are illustrated by the example of processing the survey data of the Republic of Rwanda in 2017. The simultaneous inversion of data in frequency domain and time domain can significantly improve the quality of interpretation.

scholarly journals Time and Frequency Domain Dynamic Analysis of Offshore Mooring

2021 ◽  
Vol 9 (7) ◽  
pp. 781
Author(s):  
Shi He ◽  
Aijun Wang
Keyword(s):  
Dynamic Analysis ◽  
Frequency Domain ◽  
Time Domain ◽  
Linearization Method ◽  
Euler Method ◽  
Single Component ◽  
Hydrodynamic Drag ◽  
Lumped Mass ◽  
Mooring Lines ◽  
The Time Domain

The numerical procedures for dynamic analysis of mooring lines in the time domain and frequency domain were developed in this work. The lumped mass method was used to model the mooring lines. In the time domain dynamic analysis, the modified Euler method was used to solve the motion equation of mooring lines. The dynamic analyses of mooring lines under horizontal, vertical, and combined harmonic excitations were carried out. The cases of single-component and multicomponent mooring lines under these excitations were studied, respectively. The case considering the seabed contact was also included. The program was validated by comparing with the results from commercial software, Orcaflex. For the frequency domain dynamic analysis, an improved frame invariant stochastic linearization method was applied to the nonlinear hydrodynamic drag term. The cases of single-component and multicomponent mooring lines were studied. The comparison of results shows that frequency domain results agree well with nonlinear time domain results.

A Multiple Harmonic Open-Loop Controller for Hydro/Aerodynamic Force Measurements in Rotating Machinery Using Magnetic Bearings

2002 ◽  
Vol 124 (4) ◽  
pp. 827-834 ◽  
Author(s):  
D. O. Baun ◽  
E. H. Maslen ◽  
C. R. Knospe ◽  
R. D. Flack
Keyword(s):  
Frequency Domain ◽  
Time Domain ◽  
Open Loop ◽  
Rotating Machinery ◽  
Operating Conditions ◽  
Rotor Vibration ◽  
Aerodynamic Forces ◽  
Analog Input ◽  
Input And Output ◽  
The Time Domain

Inherent in the construction of many experimental apparatus designed to measure the hydro/aerodynamic forces of rotating machinery are features that contribute undesirable parasitic forces to the measured or test forces. Typically, these parasitic forces are due to seals, drive couplings, and hydraulic and/or inertial unbalance. To obtain accurate and sensitive measurement of the hydro/aerodynamic forces in these situations, it is necessary to subtract the parasitic forces from the test forces. In general, both the test forces and the parasitic forces will be dependent on the system operating conditions including the specific motion of the rotor. Therefore, to properly remove the parasitic forces the vibration orbits and operating conditions must be the same in tests for determining the hydro/aerodynamic forces and tests for determining the parasitic forces. This, in turn, necessitates a means by which the test rotor’s motion can be accurately controlled to an arbitrarily defined trajectory. Here in, an interrupt-driven multiple harmonic open-loop controller was developed and implemented on a laboratory centrifugal pump rotor supported in magnetic bearings (active load cells) for this purpose. This allowed the simultaneous control of subharmonic, synchronous, and superharmonic rotor vibration frequencies with each frequency independently forced to some user defined orbital path. The open-loop controller was implemented on a standard PC using commercially available analog input and output cards. All analog input and output functions, transformation of the position signals from the time domain to the frequency domain, and transformation of the open-loop control signals from the frequency domain to the time domain were performed in an interrupt service routine. Rotor vibration was attenuated to the noise floor, vibration amplitude ≈0.2 μm, or forced to a user specified orbital trajectory. Between the whirl frequencies of 14 and 2 times running speed, the orbit semi-major and semi-minor axis magnitudes were controlled to within 0.5% of the requested axis magnitudes. The ellipse angles and amplitude phase angles of the imposed orbits were within 0.3 deg and 1.0 deg, respectively, of their requested counterparts.

Time-Domain Nonlinear SSI Analysis of Foundation Sliding Using Frequency-Dependent Foundation Impedance Derived From SASSI

2008 ◽  
Author(s):  
Mansour Tabatabaie ◽  
Thomas Ballard
Keyword(s):  
Frequency Domain ◽  
Linear Systems ◽  
Time Domain ◽  
Soil Structure ◽  
Wave Radiation ◽  
Far Field ◽  
Frequency Dependent ◽  
Nonlinear Springs ◽  
Mat Foundation ◽  
The Time Domain

Dynamic soil-structure interaction (SSI) analysis of nuclear power plants is often performed in frequency domain using programs such as SASSI [1]. This enables the analyst to properly a) address the effects of wave radiation in an unbounded soil media, b) incorporate strain-compatible soil shear modulus and damping properties and c) specify input motion in the free field using the de-convolution method and/or spatially variable ground motions. For structures that exhibit nonlinearities such as potential base sliding and/or uplift, the frequency-domain procedure is not applicable as it is limited to linear systems. For such problems, it is necessary to solve the problem in the time domain using the direct integration method in programs such as ADINA [2]. The authors recently introduced a sub-structuring technique called distributed parameter foundation impedance (DPFI) model that allows the structure to be partitioned from the total SSI system and analyzed in the time domain while the foundation soil is modeled using the frequency-domain procedure [3]. This procedure has been validated for linear systems. In this paper we have expanded the DPFI model to incorporate nonlinearities at the soil/structure interface by introducing nonlinear shear and normal springs arranged in series between the DPFI and structure model. This combination of the linear far-field impedance (DPFI) plus nonlinear near-field soil springs allows the foundation sliding and/or uplift behavior be analyzed in time domain while maintaining the frequency-dependent stiffness and radiation damping nature of the far-field foundation impedance. To check the accuracy of this procedure, a typical NPP foundation mat supported at the surface of a layered soil system and subjected to harmonic forced vibration was first analyzed in the frequency domain using SASSI to calculate the target linear response and derive a linear, far-field DPFI model. The target linear solution was then used to validate two linear time-domain ADINA models: Model 1 consisting of the mat foundation+DPFI derived from the linear SASSI model and Model 2 consisting of the total SSI system (mat foundation plus a soil block). After linear alignment, the nonlinear springs were added to both ADINA models and re-analyzed in time domain. Model 2 provided the target nonlinear solution while Model 1 provided the results using the DPFI+nonlinear springs. By increasing the amplitude of the vibration load, different levels of foundation sliding were simulated. Good agreement between the results of two models in terms of the displacement response of the mat and cyclic force-displacement behavior of the springs validates the accuracy of the procedure presented herein.

IDENTIFICATION OF INTERHARMONICS BASED ON WAVELET PACKET TRANSFORM

2021 ◽  
Vol 3 (1) ◽  
pp. 031-036
Author(s):  
S. A. GOROVOY ◽  
◽  
V. I. SKOROKHODOV ◽  
D. I. PLOTNIKOV ◽  
◽  
...  
Keyword(s):  
Fourier Transform ◽  
Simulation Model ◽  
Frequency Domain ◽  
Time Domain ◽  
Wavelet Packet ◽  
High Accuracy ◽  
Wavelet Packet Transform ◽  
Mathematical Apparatus ◽  
Nonlinear Load ◽  
The Time Domain

This paper deals with the analysis of interharmonics, which are due to the presence of a nonlinear load. The tool for the analysis was a mathematical apparatus - wavelet packet transform. Which has a number of advantages over the traditional Fourier transform. A simulation model was developed in Simulink to simulate a non-stationary non-sinusoidal mode. The use of the wavelet packet transform will allow to determine the mode parameters with high accuracy from the obtained wavelet coefficients. It also makes it possible to obtain information, both in the frequency domain of the signal and in the time domain.

scholarly journals A research on fiber-optic vibration pattern recognition based on time-frequency characteristics

2018 ◽  
Vol 10 (12) ◽  
pp. 168781401881346 ◽  
Author(s):  
Tabi Fouda Bernard Marie ◽  
Dezhi Han ◽  
Bowen An ◽  
Jingyun Li
Keyword(s):  
Pattern Recognition ◽  
Frequency Domain ◽  
Time Domain ◽  
Fiber Optic ◽  
Pattern Recognition Method ◽  
Recognition Method ◽  
Time Frequency ◽  
Mode Decomposition ◽  
Vibration Pattern ◽  
The Time Domain

To detect and recognize any type of events over the perimeter security system, this article proposes a fiber-optic vibration pattern recognition method based on the combination of time-domain features and time-frequency domain features. The performance parameters (event recognition, event location, and event classification) are very important and describe the validity of this article. The pattern recognition method is precisely based on the empirical mode decomposition of time-frequency entropy and center-of-gravity frequency. It implements the function of identifying and classifying the event (intrusions or non-intrusion) over the perimeter to secure. To achieve this method, the first-level prejudgment is performed according to the time-domain features of the vibration signal, and the second-level prediction is carried out through time-frequency analysis. The time-frequency distribution of the signal is obtained by empirical mode decomposition and Hilbert transform and then the time-frequency entropy and center-of-gravity frequency are used to form the time-frequency domain features, that is, combined with the time-domain features to form feature vectors. Multiple types of probabilistic neural networks are identified to determine whether there are intrusions and the intrusion types. The experimental results demonstrate that the proposed method is effective and reliable in identifying and classifying the type of event.

scholarly journals Comparison of machine learning models based on time domain and frequency domain features for faults diagnosis in rotating machines

2018 ◽  
Vol 211 ◽  
pp. 17009
Author(s):  
Natalia Espinoza Sepulveda ◽  
Jyoti Sinha
Keyword(s):  
Machine Learning ◽  
Frequency Domain ◽  
Time Domain ◽  
Intelligent Systems ◽  
Learning Models ◽  
Machine Vibration ◽  
Vibration Data ◽  
Machine Learning Model ◽  
The Time Domain ◽  
Machine Learning Models

The development of technologies for the maintenance industry has taken an important role to meet the demanding challenges. One of the important challenges is to predict the defects, if any, in machines as early as possible to manage the machines downtime. The vibration-based condition monitoring (VCM) is well-known for this purpose but requires the human experience and expertise. The machine learning models using the intelligent systems and pattern recognition seem to be the future avenue for machine fault detection without the human expertise. Several such studies are published in the literature. This paper is also on the machine learning model for the different machine faults classification and detection. Here the time domain and frequency domain features derived from the measured machine vibration data are used separated in the development of the machine learning models using the artificial neutral network method. The effectiveness of both the time and frequency domain features based models are compared when they are applied to an experimental rig. The paper presents the proposed machine learning models and their performance in terms of the observations and results.

Sign in / Sign up

Export Citation Format

Share Document