scholarly journals Three-Dimensional Anisotropic Inversions for Time-Domain Airborne Electromagnetic Data

Minerals ◽  
2021 ◽  
Vol 11 (2) ◽  
pp. 218
Author(s):  
Yang Su ◽  
Changchun Yin ◽  
Yunhe Liu ◽  
Xiuyan Ren ◽  
Bo Zhang ◽  
...  

Rocks and ores in nature usually appear macro-anisotropic, especially in sedimentary areas with strong layering. This anisotropy will lead to false interpretation of electromagnetic (EM) data when inverted under the assumption of an isotropic earth. However, the time-domain (TD) airborne EM (AEM) inversion for an anisotropic model has not attracted much attention. To get reasonable inversion results from TD AEM data, we present in this paper the forward modeling and inversion methods based on a triaxial anisotropic model. We apply three-dimensional (3D) finite-difference on the secondary scattered electric field equation to calculate the frequency-domain (FD) EM responses, then we use the inverse Fourier transform and waveform convolution to obtain TD responses. For the regularized inversion, we calculate directly the sensitivities with respect to three diagonal conductivities and then use the Gauss–Newton (GN) optimization scheme to recover model parameters. To speed up the computation and to reduce the memory requirement, we adopt the moving footprint concept and separate the whole model into a series of small sub-models for the inversion. Finally, we compare our anisotropic inversion scheme with the isotropic one using both synthetic and field data. Numerical experiments show that the anisotropic inversion has inherent advantages over the isotropic ones, we can get more reasonable results for the anisotropic earth structures.

Geophysics ◽  
2012 ◽  
Vol 77 (4) ◽  
pp. WB59-WB69 ◽  
Author(s):  
Leif H. Cox ◽  
Glenn A. Wilson ◽  
Michael S. Zhdanov

Time-domain airborne surveys gather hundreds of thousands of multichannel, multicomponent samples. The volume of data and other complications have made 1D inversions and transforms the only viable method to interpret these data, in spite of their limitations. We have developed a practical methodology to perform full 3D inversions of entire time- or frequency-domain airborne electromagnetic (AEM) surveys. Our methodology is based on the concept of a moving footprint that reduces the computation requirements by several orders of magnitude. The 3D AEM responses and sensitivities are computed using a frequency-domain total field integral equation technique. For time-domain AEM responses and sensitivities, the frequency-domain responses and sensitivities are transformed to the time domain via a cosine transform and convolution with the system waveform. We demonstrate the efficiency of our methodology with a model study relevant to the Abitibi greenstone belt and a case study from the Reid-Mahaffy test site in Ontario, Canada, which provided an excellent practical opportunity to compare 3D inversions for different AEM systems. In particular, we compared 3D inversions of VTEM-35 (time-domain helicopter), MEGATEM II (time-domain fixed-wing), and DIGHEM (frequency-domain helicopter) data. Our comparison showed that each system is able to image the conductive overburden and to varying degrees, detect and delineate the bedrock conductors, and, as expected, that the DIGHEM system best resolved the conductive overburden, whereas the time-domain systems most clearly delineated the bedrock conductors. Our comparisons of the helicopter and fixed-wing time-domain systems revealed that the often-cited disadvantages of a fixed-wing system (i.e., response asymmetry) are not inherent in the system, but rather reflect a limitation of the 1D interpretation methods used to date.


2021 ◽  
Vol 11 (1) ◽  
Author(s):  
Souma Jinno ◽  
Shuji Kitora ◽  
Hiroshi Toki ◽  
Masayuki Abe

AbstractWe formulate a numerical method on the transmission and radiation theory of three-dimensional conductors starting from the Maxwell equations in the time domain. We include the delay effect in the integral equations for the scalar and vector potentials rigorously, which is vital to obtain numerically stable solutions for transmission and radiation phenomena in conductors. We provide a formalism to connect the conductors to any passive lumped-parameter circuits. We show one example of numerical calculations, demonstrating that the new formalism provides stable solutions to the transmission and radiation phenomena.


2018 ◽  
Vol 2018 ◽  
pp. 1-11
Author(s):  
Chen Xia ◽  
Chengzhi Qi ◽  
Xiaozhao Li

Transmitting boundaries are important for modeling the wave propagation in the finite element analysis of dynamic foundation problems. In this study, viscoelastic boundaries for multiple seismic waves or excitations sources were derived for two-dimensional and three-dimensional conditions in the time domain, which were proved to be solid by finite element models. Then, the method for equivalent forces’ input of seismic waves was also described when the proposed artificial boundaries were applied. Comparisons between numerical calculations and analytical results validate this seismic excitation input method. The seismic response of subway station under different seismic loads input methods indicates that asymmetric input seismic loads would cause different deformations from the symmetric input seismic loads, and whether it would increase or decrease the seismic response depends on the parameters of the specific structure and surrounding soil.


Geophysics ◽  
2015 ◽  
Vol 80 (6) ◽  
pp. K25-K36 ◽  
Author(s):  
Michael S. McMillan ◽  
Christoph Schwarzbach ◽  
Eldad Haber ◽  
Douglas W. Oldenburg

Author(s):  
Changkun Wei ◽  
Jiaqing Yang ◽  
Bo Zhang

In this paper, we propose and study the uniaxial perfectly matched layer (PML) method for three-dimensional time-domain electromagnetic scattering problems, which has a great advantage over the spherical one in dealing with problems involving anisotropic scatterers. The truncated uniaxial PML problem is proved to be well-posed and stable, based on the Laplace transform technique and the energy method. Moreover, the $L^2$-norm and $L^{\infty}$-norm error estimates in time are given between the solutions of the original scattering problem and the truncated PML problem, leading to the exponential convergence of the time-domain uniaxial PML method in terms of the thickness and absorbing parameters of the PML layer. The proof depends on the error analysis between the EtM operators for the original scattering problem and the truncated PML problem, which is different from our previous work (SIAM J. Numer. Anal. 58(3) (2020), 1918-1940).


Author(s):  
Yoshiyuki Inoue ◽  
Md. Kamruzzaman

The LNG-FPSO concept is receiving much attention in recent years, due to its active usage to exploit oil and gas resources. The FPSO offloads LNG to an LNG carrier that is located close to the FPSO, and during this transfer process two large vessels are in close proximity to each other for daylong periods of time. Due to the presence of neighboring vessel, the motion response of both the vessels will be affected significantly. Hydrodynamic interactions related to wave effects may result in unfavorable responses or the risk of collisions in a multi-body floating system. Not only the motion behavior but also the second order drift forces are influenced by the neighboring structures due to interactions of the waves among the structures. A study is made on the time domain analysis to assess the behavior and the operational capability of the FPSO system moored in the sea having an LNG carrier alongside under environmental conditions such as waves, wind and currents. This paper presents an analysis tool to predict the dynamic motion response and non-linear connecting and mooring forces on a parallel-connected LNG-FPSO system due to non-linear exciting forces of wave, wind and current. Simulation for the mooring performance is also investigated. The three-dimensional source-sink technique has been applied to obtain the radiation forces and the transfer function of wave exciting forces on floating multi-bodies. The hydrodynamic interaction effect between the FPSO and the LNG carrier is included to calculate the hydrodynamic forces. For the simulation of a random sea and also for the generation of time depended wind velocity, a fully probabilistic simulation technique has been applied. Wind and current loads are estimated according to OCIMF. The effects of variations in wave, wind and current loads and direction on the slowly varying oscillations of the LNG and FPSO are also investigated in this paper. Finally, some conclusions are drawn based on the numerical results obtained from the present time domain simulations.


2019 ◽  
Vol 871 ◽  
pp. 350-376 ◽  
Author(s):  
Agnès Maurel ◽  
Kim Pham ◽  
Jean-Jacques Marigo

We study the propagation of water waves over a ridge structured at the subwavelength scale using homogenization techniques able to account for its finite extent. The calculations are conducted in the time domain considering the full three-dimensional problem to capture the effects of the evanescent field in the water channel over the structured ridge and at its boundaries. This provides an effective two-dimensional wave equation which is a classical result but also non-intuitive transmission conditions between the region of the ridge and the surrounding regions of constant immersion depth. Numerical results provide evidence that the scattering properties of a structured ridge can be strongly influenced by the evanescent fields, a fact which is accurately captured by the homogenized model.


Geophysics ◽  
1999 ◽  
Vol 64 (1) ◽  
pp. 278-288
Author(s):  
Chengshu Wang

I consider a new dip‐moveout (DMO) processing technique in the Radon domain called Radon DMO. The Radon DMO operator directly maps data from the NMO-corrected time domain to the DMO wavefield in the Radon domain. The method is built upon a process that transforms a single NMO-corrected trace into multiple traces spread along hyperbolas in the Radon domain. These hyperbolas are a linear Radon map of the DMO ellipses in the time domain. In this paper, I introduce the amplitude‐preserving Radon DMO and compare some examples of Radon DMO and Fourier DMO for both synthetic and real data. I also show the better frequency preservation properties of the Radon DMO method. Three‐dimensional data are often irregularly sampled with respect to fold, azimuth, and offset. Population deficiencies are exaggerated in the common‐offset domain. Radon DMO does not require that input traces belong to one common‐offset bin as does the Fourier method. Input traces can be organized from multiple offset bins grouping to perform Radon DMO, which is well used in 3-D surveys. Some synthetic and real data examples show these properties.


Sign in / Sign up

Export Citation Format

Share Document