• Tunnel and Underground Space
• /
• v.23 no.6
• /
• pp.493-505
• /
• 2013

This numerical study investigates seismicity and fault slip induced by fluid injection in deep geothermal reservoir with pre-existing fractures and fault. Particle Flow Code 2D is used with additionally implemented hydro-mechanical coupled fluid flow algorithm and acoustic emission moment tensor inversion algorithm. The output of the model includes spatio-temporal evolution of induced seismicity (hypocenter locations and magnitudes) and fault deformation (failure and slip) in relation to fluid pressure distribution. The model is applied to a case of fluid injection with constant rates changing in three steps using different fluid characters, i.e. the viscosity, and different injection locations. In fractured reservoir, spatio-temporal distribution of the induced seismicity differs significantly depending on the viscosity of the fracturing fluid. In a fractured reservoir, injection of low viscosity fluid results in larger volume of induced seismicity cloud as the fluid can migrate easily to the reservoir and cause large number and magnitude of induced seismicity in the post-shut-in period. In a faulted reservoir, fault deformation (co-seismic failure and aseismic slip) can occur by a small perturbation of fracturing fluid (<0.1 MPa) can be induced when the injection location is set close to the fault. The presented numerical model technique can practically be used in geothermal industry to predict the induced seismicity pattern and magnitude distribution resulting from hydraulic stimulation of geothermal reservoirs prior to actual injection operation.

• Geophysics and Geophysical Exploration
• /
• v.20 no.4
• /
• pp.216-225
• /
• 2017

Recently, steel casing became an interesting issue when applying controlled-source electromagnetic (EM) method to various fields because effects of steel casing on EM responses are not negligible. This study employed an approach that approximates the steel casing as a series of electric dipole sources in order to develop the numerical algorithm for the efficient simulation of EM responses in the presence of steel casing. After verifying the validity of the developed algorithm, we analyze effects of steel casing on EM responses with the synthetic model simulating geothermal reservoir environment. The analysis showed that the effects of steel casing on EM responses are localized near the casing and increase as the transmitter becomes close to the casing. In addition, through the analysis on the EM responses by the injection of clean water, we confirm that the effects of casing are negligible when interpreting the after-injection data acquired using the transmitter located far enough from the casing. Considering the difference in EM responses between before and after injection in inversion, the effects of the casing can be neglected although after-injection data shows considerable difference due to the close distance between the transmitter and casing. To investigate this kind steel casing effect, the precise analysis on EM responses should be preceded. The algorithm introduced in this study will contribute to the reliable calculations of EM responses distorted by the conductive steel casing.

• Geophysics and Geophysical Exploration
• /
• v.4 no.3
• /
• pp.70-79
• /
• 2001

The integral equation method is a powerful tool for numerical electromagnetic modeling. But the difficulty of this technique is the size of the linear equations, which demands excessive memory and calculation time to invert. This limitation of the integral equation method becomes critical in inverse problem. The conventional Born approximation, where the electric field in the anomalous body is approximated by the background field, is very rapid and easy to compute. However, the technique is inaccurate when the conductivity contrast between the body and the background medium is large. Quasi-linear, quasi-analytical and extended Born approximations are novel approaches to 3-D EM modeling based on the linearization of the integral equations for scattered EM field. These approximation methods are much less time consuming than full integral equation method and more accurate than conventional Born approximation. They we, however, still approximate methods for 3-D EM modeling. Iterative series methods such as modified Born, quasi-linear and quasi-analytical can be used to increase the accuracy of various approximation methods. Comparisons of numerical performance against a full integral equation and various approximation codes show that the iterative series methods are very accurate and almost always converge. Furthermore, they are very fast and easy to implement on a computer. In this study, extended Born series method is developed and it shows more accurate result than that of other series methods. Therefore, Iterative series methods, including extended Born series, open principally new possibilities for fast and accurate 3-D EM modeling and inversion.

• Journal of the Korean Geotechnical Society
• /
• v.21 no.4
• /
• pp.71-81
• /
• 2005

Surface wave techniques were initially based on 2-D plane waves and were later improved to the techniques based the 3-D based cylindrical waves. However, body-wave interference, near-field effect and limited technology in surface wave measurements restricted the use of 3-D cylindrical waves to the 1-D evaluation of subgrade stiffness. In this study, by the numerical simulation of SASW measurements, the dispersion properties of surface waves including vertical, horizontal Rayleigh waves and Love waves were thoroughly investigated in the 3-D domain, and a new filter criteria to minimize the near-field effect was established, which led to CAP (common-array-profiling)-SASW technique. The CAP-SASW technique enabled the evaluation of subgrade stiffness fur a specific subgrade segment, not for a whole section of measurement array. Therefore, a contour plot of subgrade stiffness with a ground-truth quality can be obtained by the CAP-SASW technique. The procedure proposed in this study was verified by comparing the shear-wave velocity profiles with the shear-wave velocity profiles of downhole testing at two geotechnical sites.

• Journal of the Korean earth science society
• /
• v.28 no.7
• /
• pp.936-946
• /
• 2007

Analytical geopotential field in balance with the sectoral mode (the first symmetric mode with respect to the equator) of the Rossby-Haurwitz wave on the inclined rotation axis was derived in presence of superrotation background flow. The balanced field was obtained by inverting the divergence equation with the time derivative being zero. The inversion consists of two steps, i.e., the evaluation of nonlinear forcing terms and the finding of analytical solutions based on the Poisson's equation. In the second step, the forcing terms in the from of Legendre function were readily inverted due to the fact that Legendre function is the eigenfunction of the spherical Laplacian operator, while other terms were solved either by introducing a trial function or by integrating the Legendre equation. The balanced field was found to be expressed with six zonal wavenumber components, and shown to be of asymmetric structure about the equator. In association with asymmetricity, the advantageous point of the balanced field as a validation method for the numerical model was addressed. In special cases where the strength of the background flow is a half of or exactly the same as the rotation rate of the Earth it was revealed that one of the zonal wavenumber components vanishes. The analytical balanced field was compared with the geopotential field which was obtained using a spherical harmonics spectral model. It was found that the normalized difference lied in the order of machine rounding, indicating the reliability of the analytical results. The stability of the sectoral mode of Rossby-Haurwitz wave and the associated balanced field was discussed, comparing with the flrst antisymmetric mode.

• Geophysics and Geophysical Exploration
• /
• v.17 no.3
• /
• pp.129-136
• /
• 2014

Electromagnetic (EM) methods are generally divided into frequency-domain EM (FDEM) and time-domain EM (TDEM) methods, depending on the source waveform. The FDEM and TDEM fields are mathematically related by the Fourier transformation, and the TDEM field can thus be obtained as the Fourier transformation of FDEM data. For modeling in time-domain, we can use fast frequency-domain modeling codes and then convert the results to the time domain with a suitable numerical method. Thus, frequency-to-time transformations are of interest to EM methods, which is generally attained through fast Fourier transform. However, faster frequency-to-time transformation is required for the 3D inversion of TDEM data or for the processing of vast air-borne TDEM data. The diffusion expansion method (DEM) is one of smart frequency-to-time transformation methods. In DEM, the EM field is expanded into a sequence of diffusion functions with a known frequency dependence, but with unknown diffusion-times that must be chosen based on the data to be transformed. Especially, accuracy of DEM is sensitive to the diffusion-time. In this study, we developed a method to determine the optimum range of diffusion-time values, minimizing the RMS error of the frequency-domain data approximated by the diffusion expansion. We confirmed that this method produces accurate results over a wider time range for a homogeneous half-space and two-layered model.

• Korean Journal of Optics and Photonics
• /
• v.8 no.6
• /
• pp.445-449
• /
• 1997

The complex refractive indices of $Ge_2Se_2Te_5$ which show reversible phase change between the crystalline phase and an amorphous one depending upon the annealing process have been determined in the spectral range of 0.7-4.5 eV. The $Ge_2Se_2Te_5$ films were DC sputter deposited on the crystalline silicon substrate. The spectro-ellipsometry data of a thick film were analyzed following the modelling procedure where the quantum mechanical dispersion relation were used for the complex refractive indices of both the cryastalline phase $Ge_2Se_2Te_5$ and and amorphous phase $Ge_2Se_2Te_5$, respectively. On the other hand, with the surface micro-roughness layer whose effective thickness was determined from AFM analysis, the spectro-ellipsometry data were numerically inverted to yield the complex refractive index of $Ge_2Se_2Te_5$ at each wavelength. With these set of complex refractive indices, the reflectance spectra were calculated and those spectra obtained from the numerical inversion showed better agreement with the experimental reflection spectra for both the cryastalline phase and an amorphous phase. Finally, the thin $Ge_2Se_2Te_5$ film which has the optimum thickness of 26 nm as the medium for optical recording was also analyzed and the quantitative result of the film thickness and the surface microroughness has been reported.

• Geophysics and Geophysical Exploration
• /
• v.11 no.2
• /
• pp.85-92
• /
• 2008

Magnetotelluric (MT) methods are widely applied as an effective exploration technique to geothermal surveys. Two-dimensional (2-D) analysis is frequently used to investigate a complicated subsurface structure in a geothermal region. A 2-D finite-element method (FEM) is usually applied to the MT analysis, but we must pay attention to the accuracy of so-called auxiliary fields. Rodi (1976) proposed an algorithm of improving the accuracy of auxiliary fields, and named it as the MOM method. Because it introduces zeros into the diagonal elements of coefficient matrix of the FEM total equation, a pivoting procedure applied to the symmetrical band matrix makes the numerical solution far less efficient. The MOM method was devised mainly for the inversion analysis, in which partial derivatives of both electric and magnetic fields with respect to model parameters are required. In the case of forward modeling, however, we do not have to resort to the MOM method; there is no need of modifying the coefficient matrix, and the auxiliary fields can be elicited from the regular FEM solution. The computational efficiency of the MOM method, however, can be greatly improved through a sophisticated rearrangement of the total equation.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.34 no.3
• /
• pp.167-174
• /
• 2021

Recent investigation into the integrity of nuclear containment buildings has highlighted the importance of developing an elaborate diagnostic method to evaluate the distribution and size of cavities inside concrete walls. As part of developing such a method, this paper presents a finite element approach to modeling elastic waves propagating in the containment building walls of a nuclear power plant. We introduce a perfectly matched layer (PML) wave-absorbing boundary to limit the large-scale nuclear containment wall to the region of interest. The formulation results in a semi-discrete form with symmetric damping and stiffness matrices. The transient elastic wave equations for a mixed unsplit-field PML were solved for displacement and stresses in the time domain. Numerical results show that the sensitivity of displacement, velocity, acceleration, and stresses is large depending on the size and location of the cavity. The dynamic response of the wall slightly differs depending on the existence of the containment liner plate. The results of this study can be applied to a full-waveform inversion approach for characterizing cavities inside a containment wall.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.26 no.3C
• /
• pp.143-155
• /
• 2006

SPT-Uphole tomography method was introduced for the evaluation of near subsurface shear wave velocity (Vs) distribution map. In SPT-Uphole method, SPT (Standard Penetration Test) which is common in geotechnical site investigation was used as a source and several surface geophones in line were used as receivers. Vs distribution map which is the triangular shape around the boring point can be developed by tomography inversion. To obtain the exact travel time information of shear wave component, a procedure using the magnitude summation of vertical and horizontal components was used based on the evaluation of particle motion at the surface. It was verified that proposed method could give reliable Vs distribution map through the numerical study using the FEM (Finite Element Method) model. Finally, SPT-Uphole tomography method was performed at the weathered soil site where several boring data with SPT-N values are available, and the feasibility of proposed method was verified in the field.