• Geophysics and Geophysical Exploration
• /
• v.24 no.3
• /
• pp.98-112
• /
• 2021

We followed and extended the algorithm originally made by Das (1995) to calculate the electric potential and field induced by electric current in arbitrary anisotropic layered structure. We confirmed all the theoretical contents and coded the corresponding program to acquire the electric potential and field. Further we extended to forward estimation of apparent resistivity to be attained by electrical resistivity survey on anisotropic layered structure with differing the electrode spacing and azimuth of anisotropy. The effects of anisotropy were reviewed by considering some examples.

• Geophysics and Geophysical Exploration
• /
• v.24 no.2
• /
• pp.45-52
• /
• 2021

With advancements in software technology, more scientists and engineers are employing computer software and programming tools for research. However, several issues can arise in software-based research: environment setting, reproducibility, and loss of source codes. This study investigates the use of Linux containers and version control systems to prevent these problems. Managing research projects using a cloud source-code repository and building a research environment in a Linux container can prevent the abovementioned problems and make research collaboration easier. For researchers with no experience with Linux containers, a repository of project template containing shell scripts for building and running containers has been released.

• Geophysics and Geophysical Exploration
• /
• v.25 no.1
• /
• pp.44-49
• /
• 2022

Closed-form expressions of vector gravity and gravity gradient tensor based on a line segment are derived. If a cylindrical object with axial symmetry is observed from a distance, it is possible to approximate it as a line segment; therefore, it is necessary to compute the gravity and the gravity gradient tensor due to a line source by using closed-form expressions. The gravitational potential for a line segment is defined as a one-dimensional integral, and this integral is differentiated with respect to the Cartesian coordinate system to derive the vector gravity. The expressions of the gravity gradient tensor are derived by differentiating the vector gravity once more in the same coordinate system.

• Geophysics and Geophysical Exploration
• /
• v.25 no.4
• /
• pp.242-251
• /
• 2022

Automatic differentiation automatically calculates the derivatives of a function using the chain rule once the forward operation of a function is defined. Given the recent development of computing libraries that support automatic differentiation, many researchers have adopted automatic differentiation techniques to solve geophysical inverse problems. We analyzed the advantages, disadvantages, and performances of automatic differentiation techniques using the gradient calculations of seismic full waveform inversion objective functions. The gradients of objective functions can be expressed as multiplications of the derivatives of the model parameters, wavefields, and objective functions using the chain rule. Using numerical examples, we demonstrated the speed of analytic differentiation and the convenience of complex gradient calculations for automatic differentiation. We calculated derivatives of model parameters and objective functions using automatic differentiation and derivatives of wavefields using analytic differentiation.

• Geophysics and Geophysical Exploration
• /
• v.27 no.1
• /
• pp.69-83
• /
• 2024

Considering importance and specificity, site investigations for deep geological disposal of Spent Nuclear Fuel require stringent quality control, unlike general geotechnical investigations for tunnels and bridges. In this study, we present a case of selecting geophysical survey method for individual site investigation stage and preparing geophysical survey guideline. The proposed geophysical survey guidelines include procedures, considerations, and quality control for exploration planning, data acquisition, data processing, and interpretation. They comprehensively summarize the contents of airborne electromagnetic survey and seismic reflection survey.

• Geophysics and Geophysical Exploration
• /
• v.15 no.2
• /
• pp.51-56
• /
• 2012

One of the most critical and essential procedures in the interpretation of gravity and magnetic data is to separate the anomaly due to the specific geologic structure from the summation of effects from a broad variety of geologic sources, especially those of different depths. Separation of the residual anomaly from the regional field is the most simple case of the vertical separation. If the anomaly due to a layer of specific depth can be separated or the depth of the separated layer can be quantitatively determined, it may deserve the separation-sounding. We suggest a wavelength filter whose cutoff frequency is determined by log-power spectrum analysis, as a separation-sounding filter. We applied this filter both to synthetic and real gravity data acquired at Heunghae area, and compared the results with those of Jacobsen's upward continuation filter. These showed that the proposed separation-sounding filter could be a useful tool for interpretation of the vertical geologic structure by stripping the gravity effects of geologic sources down to the desired depth.

• Geophysics and Geophysical Exploration
• /
• v.15 no.2
• /
• pp.57-65
• /
• 2012

Various numerical methods in simulation of seismic wave propagation have been developed. Recently an innovative numerical method called as the Spectral Element Method (SEM) has been developed and used in wave propagation in 3-D elastic media. The SEM that easily implements the free surface of topography combines the flexibility of a finite element method with the accuracy of a spectral method. It is generally used a weak formulation of the equation of motion which are solved on a mesh of hexahedral elements based on the Gauss-Lobatto-Legendre integration rule. Variational formulations of velocity-stress motion are newly modified in order to implement ADE-PML (Auxiliary Differential Equation of Perfectly Matched Layer) in wave propagation in 3-D elastic media, because a general weak formulation has a difficulty in adapting CFS (Complex Frequency Shifted) PML (Perfectly Matched Layer). SEM of Velocity-Stress motion having ADE-PML that is very efficient in absorbing waves reflected from finite boundary is verified with simulation of 1-D and 3-D wave propagation.

• Geophysics and Geophysical Exploration
• /
• v.15 no.2
• /
• pp.85-91
• /
• 2012

Multiple Lapse Time Window (MLTW) analysis for obtaining intrinsic attenuation value require numerous data without directional bias to compensate focal mechanism. The first window of MLTW, therefore, shows large deviation in fitting smoothed theoretical curve. The information on the focal mechanism may reduce burdens of number and distribution. This study combined algorithm of computing focal mechanism to DSMC method by Yoshimoto (2000). However, the MLTW method based on the numerous data was not applicable to this study, because of the limited data to the almost same fault plane solution. This study showed that the available data was too insufficient to construct smoothed theoretical curve, although the deviation of the first window was improved. Instead of conventional solution by more data, the study seems to be needed for new constraints to obtain smoothed curve.

• Geophysics and Geophysical Exploration
• /
• v.16 no.3
• /
• pp.139-144
• /
• 2013

A computer program has been developed to estimate time-domain electromagnetic (EM) responses for a onedimensional model with multiple source and receiver dipoles that are finite in length. The time-domain solution can be obtained by applying an inverse fast Fourier transform (FFT) to frequency-domain fields for efficiency. Frequency-domain responses are first obtained for 10 logarithmically equidistant frequencies per decade, and then cubic spline interpolated to get the FFT input. In the case of phases, the phase curve must be made to be continuous prior to the spline interpolation. The spline interpolated data are convolved with a source current waveform prior to FFT. In this paper, only a step-off waveform is considered. This time-domain code is verified with an analytic solution and EM responses for a marine hydrocarbon reservoir model. Through these comparisons, we can confirm that the accuracy of the developed program is fairly high.

• Geophysics and Geophysical Exploration
• /
• v.23 no.4
• /
• pp.223-229
• /
• 2020

Small-loop EM techniques have been used in many geophysical investigations, including shallow engineering and environmental surveys. Even though geometric and frequency soundings have been widely used, there is a debate regarding the effectiveness of frequency sounding, especially when the coil spacing is small. In this study, we analyzed the effectiveness of geometric as well as frequency soundings via the one-dimensional modeling of small-loop EM surveys. The numerical results reveal that geometric sounding can effectively provide underground information. Conversely, the frequency soundings are only effective when the loop spacing is relatively large, that is, when the induction number is large. On the contrary, the frequency soundings fail to provide any information concerning the subsurface layers if the loop spacing is not large.