• 한국지구물리탐사학회:학술대회논문집
• /
• 2007.06a
• /
• pp.203-208
• /
• 2007

Due to the long tectonic history and the very complex geologic formations in Korea, the anisotropic characteristics of subsurface material may often change very greatly and locally. The algorithms for the travel time computation commonly used, however, may not give sufficiently precise results particularly for the complex and strong anisotropic model, since they are based on the two-dimensional (2D) earth and/or weak anisotropy assumptions. This study is intended to develope a three-dimensional (3D) modeling algorithm to precisely calculate the first arrival time in the complex anisotropic media. We assume 3D TTI (tilted transversely isotropy) medium having the arbitrary symmetry axis. The algorithm includes the 2D non-linear interpolation scheme to calculate the traveltimes inside the grid and the 3D traveltime mapping to fill the 3D model with first arrival times. The weak anisotropy assumption, moreover, can be overcome through devising a numerical approach of the steepest descent method in the calculation of minimum traveltime, instead of using approximate solution.

• Geophysics and Geophysical Exploration
• /
• v.9 no.3
• /
• pp.241-249
• /
• 2006

Due to the long tectonic history and the very complex geologic formations in Korea, the anisotropic characteristics of subsurface material may often change very greatly and locally. The algorithms commonly used, however, may not give sufficiently precise computational results of traveltime data particularly for the complex and strong anisotropic model, since they are based on the two-dimensional (2D) earth and/or weak anisotropy assumptions. This study is intended to develope a three-dimensional (3D) modeling algorithm to precisely calculate the first arrival time in the complex anisotropic media. Considering the complex geology of Korea, we assume 3D TTI (tilted transversely isotropy) medium having the arbitrary symmetry axis. The algorithm includes the 2D non-linear interpolation scheme to calculate the traveltimes inside the grid and the 3D traveltime mapping to fill the 3D model with first arrival times. The weak anisotropy assumption, moreover, can be overcome through devising a numerical approach of the steepest descent method in the calculation of minimum traveltime, instead of using approximate solution. The performance of the algorithm developed in this study is demonstrated by the comparison of the analytic and numerical solutions for the homogeneous anisotropic earth as well as through the numerical experiment for the two layer model whose anisotropic properties are greatly different each other. We expect that the developed modeling algorithm can be used in the development of processing and inversion schemes of seismic data acquired in strongly anisotropic environment, such as migration, velocity analysis, cross-well tomography and so on.

• Geophysics and Geophysical Exploration
• /
• v.19 no.3
• /
• pp.111-120
• /
• 2016

In order to enhance the connectivity of fracture network as fluid path in enhanced/engineered geothermal system (EGS), the exact locating of hydraulic fractured zone is very important. Hydraulic fractures can be tracked by locating of microseismic events which are occurred during hydraulic fracture stimulation at each stage. However, since the subsurface velocity is changed due to hydraulic fracturing at each stage, in order to find out the exact location of microseismic events, we have to consider the velocity change due to hydraulic fracturing at previous stage when we perform the mapping of microseimic events at the next stage. In this study, we have modified 3D locating algorithm of microseismic data which was developed by Kim et al. (2015) and have developed 3D velocity update algorithm using occurred microseismic data. Eikonal equation which can efficiently calculate traveltime for complex velocity model at anywhere without shadow zone is used as forward engine in our inversion. Computational cost is dramatically reduced by using Fresnel volume approach to construct Jacobian matrix in velocity inversion. Through the numerical test which simulates the geothermal survey geometry, we demonstrated that the initial velocity model was updated by using microseismic data. In addition, we confirmed that relocation results of microseismic events by using updated velocity model became closer to true locations.