• 한국지구물리탐사학회:학술대회논문집
• /
• 2006.06a
• /
• pp.9-14
• /
• 2006

Estimations of depth, magnetic orientation, and strength of dipole moments aid discrimination between unexploded ordnance (UXO) and non-UXO using magnetic surveys. Such estimations may be hindered by geologic noise, magnetic clutter, and overlapping tails of nearby dipole fields. An improved method of inversion for anomalies of single or multiple dipoles with arbitrary polarization was developed to include intra-inversion filtering and estimation of background field gradients. Data interpolated to grids are flagged so that only nodes nearest to measurement stations are used. To apply intra-inversion filtering to such data requires a gapped filter. Moreover, for data with significant gaps in coverage, or along the edges or corners of survey areas, intra-inversion filters must be appropriately modified. To that end, edge-adaptive and gapped gradient-nulling filters have been designed and tested. Applications are shown for magnetic field data from Chongcho Lake, Sokcho, Korea and the U. S. Army's Aberdeen Proving Ground in Maryland.

• Geophysics and Geophysical Exploration
• /
• v.11 no.4
• /
• pp.343-349
• /
• 2008

We have constructed a car-borne magnetic exploration system, in which a car drags a non-magnetic cart on which a magnetometer is installed. In the total magnetic field measured as a vectorial sum in this system, are included the magnetic field generated by the car itself. This magnetic field, doing the role of a magnetic noise, should be eliminated. For this purpose, we have set up a measurement condition to get the same effect as if we have put the car in one point and thereafter measured the magnetic field around it. In this case, if there is any magnetically anomalous body in the area, we can consider all the remaining magnetic field to have been generated by the car itself, once the geomagnetic field eliminated. We tried to invert the magnetic field considered to have been generated by the car and succeeded to derive the magnetic moment and the direction of the induced and remanent magnetic field of the car respectively. Once the magnetic moment and the direction of the induced and remanent magnetic field have been calculated, the magnetic field generated by them in specific points can be directly and analytically calculated. This result can be used in the future to eliminate the magnetic field generated by the car itself doing the role of a magnetic noise during the procedure of reduction of the measured magnetic exploration data by the car-borne magnetic exploration system.

• 한국지구물리탐사학회:학술대회논문집
• /
• 2006.06a
• /
• pp.141-151
• /
• 2006

A blank test was done to calculatee the car itself's magnetic effect as noise and to eliminate it from the data set of total magnetic intensity(=magnetic flux density) exploration in a car-borne magnetic exploration system. To calculate the induced magnetic intensity(= magnetization) and the remanent magnetic intensity(= magnetization) of the car itself, we have installed the magnetometer on a fixed point and measured the magnetic intensity letting the car move around the magnetometer, and we have changed the data set into an analogous data set as if acquired in the condition that we have parked the car on the same fixed point and measured the magnetic intensity moving the magnetometer around the magnetometer. Through an inversion with the later data set as input, we have calculated the magnetic center and the magnetic moments of the induced magnetic intensity(= magnetization) and the remanent magnetic intensity(= magnetization) of the car itself with the two centers coincided because of some barriers of the inversion algorithm that we have used in this study. On the other hand, we have extracted the magnetic anomaly by reducing i. e. vectorially eliminating the induced magnetic intensity(= magnetization) and the remanent magnetic intensity(= magnetization) of the car itself calculated forwardly, from the magnetic exploration data set acquired by the car-borne magnetic exploration system.

• Geophysics and Geophysical Exploration
• /
• v.2 no.1
• /
• pp.43-53
• /
• 1999

Numerical modeling and inversion for electromagnetic exploration methods are essential to understand behaviour of electromagnetic fields in complex subsurface. In this study, a finite element method was adopted as a numerical scheme for the 2.5-dimensional forward problem. And a finite element equation considering linear conductivity variation was proposed, when 2.5-dimensional differential equation to couple eletric and magnetic field was implemented. Model parameters were investigated for near-field with large source effects and far-field with responses dominantly by homogeneous half-space. Numerical responses by this study were compared with analytic solutions in homogeneous half-space. Blocky inversion model was modified to be applied to the forward calculation in this study and it was also adopted in the inversion algorithm. Resolution for isolated bodies were investigated to confirm possibility and limitation of inversion for electromagnetic exploration data.

• Geophysics and Geophysical Exploration
• /
• v.17 no.2
• /
• pp.95-103
• /
• 2014

This paper presents a practical inversion method for recovering a three-dimensional (3D) resistivity model and static shifts simultaneously. Although this method is based on a Gauss-Newton approach that requires a sensitivity matrix, the computer time can be greatly reduced by implementing a simple and effective procedure for updating the sensitivity matrix using the Broyden's algorithm. In this research, we examine the approximate inversion procedure and the weighting factor ${\beta}$ for static shifts through inversion experiments using synthetic MT data. In methods using the full sensitivity matrix constructed only once in the iteration process, a procedure using the full sensitivity in the earlier stage is useful to produce the smallest rms data misfit. The choice of ${\beta}$ is not critical below some threshold value. Synthetic examples demonstrate that the method proposed in this paper is effective in reconstructing a 3D resistivity structure from static-shifted MT data.

• Geophysics and Geophysical Exploration
• /
• v.1 no.3
• /
• pp.155-160
• /
• 1998

EM tomography technique has been developed. The algorithm used the extended Born approximations for forward modeling and reconstructed a conductivity image by a smoothness constraint least squares inversion method. Observed data, the vertical components of secondary magnetic fields, were simulated with the 3-D integral equation code. The results showed that the location of anomalous body could be imaged very well, but conductivity of the body was lower than real one and the vertical resolution was much higher than the horizontal resolution.

• Geophysics and Geophysical Exploration
• /
• v.27 no.2
• /
• pp.129-143
• /
• 2024

In magnetotelluric (MT) surveys, small inhomogeneities near the surface cause a static shift in which apparent resistivities shift regardless of frequency. As the static shift in MT data leads to errors in subsurface structure interpretation, many studies have been conducted over the past few decades to mitigate or remove the distortions it caused. The most representative method involves removing static shifts from the data before inversion. Conversely, static shifts can be corrected during inversion or included in the inversion process. In addition, other geophysical data can be used to remove static shifts. However, the correction methods are limited to one-dimensional (1D) static responses, and limitations remain in two- or three-dimensional (2D or 3D) interpretation of distorted MT data owing to static shifts. This study provides a foundation for future studies on static shift by analyzing several previously published methods.

• Geophysics and Geophysical Exploration
• /
• v.15 no.2
• /
• pp.66-74
• /
• 2012

This paper presents development of a three-dimensional marine controlled-source electromagnetic (mCSEM) modeling algorithm and its application to a salt and reservoir model to examine detectability of mCSEM for a reservoir under complex subsurface structures. The algorithm is based on the finite difference method, and employs the secondary field formulation for an accurate and fast calculation of modeling responses. The algorithm is verified for a two-layer model by comparing solutions not only with analytic solutions but also with those from other 3D modeling algorithm. We calculate and analyze electric and magnetic fields and their normalized responses for a salt and reservoir model due to three sources located at boundaries between a salt, a reservoir, and background. Numbers and positions of resistive anomalies are informed by normalized responses for three sources, and types of resistive anomalies can be informed when there is a priori information about a salt by seismic exploration.

• Geophysics and Geophysical Exploration
• /
• v.11 no.3
• /
• pp.184-196
• /
• 2008

We perform MT soundings in Bendigo, the northern part of Victoria, Australia, to investigate the deep subsurface geologic structure. The primary purpose of this survey is to figure out whether the discontinuity such as faults extends northward. The time series of MT signal were measured over 11 days at 71 measurement stations together with at remote reference, which help enhance the quality of impedance estimation and its interpretation. The impedances are estimated by robust processing using remote reference technique and then inverted with 2D MT 2D inversion. We can see that known faults are clearly imaged in MT 2D inversion. Comparing resistivity images from MT 2D inversion with interpreted boundary from reflection seismic exploration, two interpretations match well each other.

• Geophysics and Geophysical Exploration
• /
• v.2 no.1
• /
• pp.17-25
• /
• 1999

A traveltime tomography has been carried out by transforming electromagnetic data in frequency domain to wave-like domain. The transform uniquely relates a field satisfying a diffusion equation to an integral of the corresponding wavefield. But direct transform of frequency domain magnetic fields to wave-field domain is ill-posed problem because the kernel of the integral transform is highly damped. In this study, instead of solving such an unstable problem, it is assumed that wave-fields in transformed domain can be approximated by sum of ray series. And for further simplicity, reflection and refraction energy compared to that of direct wave is weak enough to be neglected. Then first arrival can be approximated by calculating the traveltime of direct wave only. But these assumptions are valid when the conductivity contrast between background medium and the target anomalous body is low enough. So this approach can only be applied to the models with low conductivity contrast. To verify the algorithm, traveltime calculated by this approach was compared to that of direct transform method and exact traveltime, calculated analytically, for homogeneous whole space. The error in first arrival picked by this study was less than that of direct transformation method, especially when the number of frequency samples is less than 10, or when the data are noisy. Layered earth model with varying conductivity contrasts and inclined dyke model have been successfully imaged by applying nonlinear traveltime tomography in 30 iterations within three CPU minutes on a IBM Pentium Pro 200 MHz.