• Geophysics and Geophysical Exploration
• /
• v.18 no.2
• /
• pp.74-84
• /
• 2015

To yield physically meaningful images through elastic reverse-time migration, the wavefield separation which extracts P- and S-waves from reconstructed vector wavefields by using elastic wave equation is prerequisite. For expanding the application of the elastic reverse-time migration to anisotropic media, not only the anisotropic modelling algorithm but also the anisotropic wavefield separation is essential. The anisotropic wavefield separation which uses pseudo-derivative filters determined according to vertical velocities and anisotropic parameters of elastic media differs from the Helmholtz decomposition which is conventionally used for the isotropic wavefield separation. Since applying these pseudo-derivative filter consumes high computational costs, we have developed the efficient anisotropic wavefield separation algorithm which has capability of parallel computing by using GPUs (Graphic Processing Units). In addition, the highly efficient anisotropic elastic reverse-time migration algorithm using MPI (Message-Passing Interface) and incorporating the developed anisotropic wavefield separation algorithm with GPUs has been developed. To verify the efficiency and the validity of the developed anisotropic elastic reverse-time migration algorithm, a VTI elastic model based on Marmousi-II was built. A synthetic multicomponent seismic data set was created using this VTI elastic model. The computational speed of migration was dramatically enhanced by using GPUs and MPI and the accuracy of image was also improved because of the adoption of the anisotropic wavefield separation.

• Tunnel and Underground Space
• /
• v.10 no.2
• /
• pp.165-172
• /
• 2000

When tunnels or underground structures are constructed in anisotropic rock mass, designers and constructors have to consider the anisotropic characteristics in rock mass because their physical and mechanical properties are depended on the anisotropic angles(${\beta}$). In this study, therefore, we have first investigated the mechanical behavior of the gneiss specimen from lab. tests, and then have analysed the behavior of specimens for to the transversely isotropic model in elastic medium using the FLAC program. The results of this study were summarized as follows; 1) In the result of the variation tests, in general, the properties of strength were depended on the angle of inclination in spite of the hard rock. And except for the shear strength test, the lowest and peak stress were appeared at 60$^{\circ}$ and 90$^{\circ}$respectively. 2) The results of specimen modeling analysis using FDM well indicated the mechanical behaviors of the specimen of transversely isotropic model.

• Journal of the Korean earth science society
• /
• v.28 no.2
• /
• pp.227-239
• /
• 2007

The MT data at the Okchon Belt show peculiar phase responses exceeding $90^{\circ}$. A reasonable explanation is that those responses are due to an electrical anisotropy structure which is composed of a narrow anisotropic block and an anisotropic layer. Considering the dominant anisotropic strikes of the block (NE-trend) and the layer (NW-trend) inferred from the MT data, if existing, the electrical anisotropy in the Okchon Belt was probably produced by the deformations in the pre-Jurassic period, since the NE-trending shearing or thrusting should create alternating bands of metamorphic rocks and fractures with NE-trending. Correlation of the structural strike of 2-D block with the latest EW-trending deformation events demonstrates that the geometrical structure of the anisotropic block was formed by the latest Daebo and Bulgugsa orogeny.

• Economic and Environmental Geology
• /
• v.36 no.1
• /
• pp.49-58
• /
• 2003

To analyze the anisotropic characteristics of fractured aquifer, variations of streaming potential were measured during and after pumping over several wells at the two test sites. Surface electrical resistivity survey, normal resistivity logging, and slug test were performed at the wells to identify the hydrogeological structure. Applying the results to the recently suggested model, the aquifer of the two test sites showed confined characteristics. Anisotropic direction appeared in using equi-potential maps from self-potential monitoring results matched well with the results of the hydrogeological test. The self-potential monitoring method adopted in this study would be useful for providing a more reliable information on the anisotropy of aquifer in the pumping test at single well.

• Geophysics and Geophysical Exploration
• /
• v.18 no.4
• /
• pp.232-240
• /
• 2015

Magnetotelluric (MT) survey investigates electrical structure of subsurface by measuring natural electromagnetic fields on the earth surface. For the accurate interpretation of MT data, the precise three-dimensional (3-D) modeling algorithm is prerequisite. Since MT responses are affected by electrical anisotropy of medium, the modeling algorithm has to incorporate the electrical anisotropy especially when analyzing time-lapse MT data sets, for monitoring engineered geothermal system (EGS) reservoir, because changes in different-vintage MT-data sets are small. This study developed a MT modeling algorithm for the simulation MT responses in the presence of electrical anisotropy by improving a pre-existing staggered-grid finite-difference MT modeling algorithm. After verifying the developed algorithm, we analyzed the effect of vertical transversely isotropic (VTI) anisotropy on MT responses. In addition, we are planning to extend the applicability of the developed algorithm which can simulate not only the horizontal transversely isotropic (HTI) anisotropy, but also the tiled transversely isotropic (TTI) anisotropy.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.22 no.2
• /
• pp.189-196
• /
• 2009

This paper presents the subparametric finite element model formulated by partial-linear layerwise theory for the analysis of laminate composites. The proposed model is based on refined approximations of two dimensional plane for orthotropic thick laminate plate as well as thin case. Three dimensional problem can be reduced to two dimensional case by assuming piecewise linear variation of in-plane displacement and a constant value of out-of-plane displacement across the thickness. The integrals of Legendre polynomials are chosen to define displacement fields and Gauss-Lobatto numerical integration is implemented in order to directly obtain maximum values occurred at the nodal points of each layer without other extrapolation techniques. The validity and characteristics of the proposed model have been tested by using orthotropic multilayered plate problem as compared to the values available in the published references. In this study, the convergence test has been carried out to determine the optimal layer model in terms of central deflection and stresses. Also, the distribution of displacements and stresses across the thickness has been investigated as the number of layer is increased.

• Proceedings of the Korean Society for Rock Mechanics Conference
• /
• 2000.09a
• /
• pp.129-137
• /
• 2000

This paper presents groundwater flow characteristics in discontinuous rock mass using fracture network program(NAPSAC) by statistical approach. Equivalent permeability coefficients are estimated from borehole data around Mabuk test tunnel site and fracture map on the arch of the tunnel. The reliability of fracture network model is obtained from determination of input data for statistical fracture network analysis from the real data(data of fracture network, data of hydraulic tests). The variation of permeability and mean anisotropic permeability coefficients are calculated from the realized model by increasing the size. As a result of analysis, a strong anisotropy of permeability is observed according to the direction of the fracture sets around the test tunnel.

• Tunnel and Underground Space
• /
• v.10 no.3
• /
• pp.378-386
• /
• 2000

This paper presents groundwater flow characteristics in discontinuous rock mass using fracture network program(NAPSAC) by statistical approach. Equivalent permeability coefficients are estimated from borehole data around Mabuk test tunnel site and fracture map on the arch of the tunnel. The reliability of fracture network model is obtained from determination of input data for statistical fracture network analysis from the real data(data of fracture network, data of hydraulic tests). The variation of permeability and mean anisotropic permeability coefficients are calculated from the realized model by increasing the size. As a result of analysis, a strong anisotropy of permeability is observed according to the direction of the fracture sets around the test tunnel.

• Computational Structural Engineering
• /
• v.8 no.3
• /
• pp.135-141
• /
• 1995

Any arbitrarily shaped laminated composite shells of revolution can be sum of the conical shell elements. Therefore, finite element model of conical shell element will be developed in this study. To verify consistency and validity of this model, natural vibrations of this model is compared with the analytical solution of cylindrical shell. Herein, an extensive parametric study is presented to assess the modeling capability of this model in class of laminated composite cylinders. It is seen that the proposed model provides highly accurate results with analytical solution. Once development of this conical shell element is done, any arbitrarily shaped composite shells of revolution can be easily analyzed.

• Geophysics and Geophysical Exploration
• /
• v.10 no.4
• /
• pp.229-240
• /
• 2007

In this study, seismic anisotropic tomography algorithm was developed for imaging the seismic velocity anisotropy of the subsurface. This algorithm includes several inversion schemes in order to make the inversion process stable and robust. First of all, the set of the inversion parameters is limited to one slowness, two ratios of slowness and one direction of the anisotropy symmetric axis. The ranges of the inversion parameters are localized by the pseudobeta transform to obtain the reasonable inversion results and the inversion constraints are controlled efficiently by ACB(Active Constraint Balancing) method. Especially, the inversion using the Fresnel volume is applied to the anisotropic tomography and it can make the anisotropic tomography more stable than ray tomography as it widens the propagation angle coverage. The algorithm of anisotropic tomography is verified through the numerical experiments. And, it is applied to the real field data measured at limestone region and the results are discussed with the drill log and geological survey data. The anisotropic tomography algorithm will be able to provide the useful tool to evaluate and understand the geological structure of the subsurface more reasonably with the anisotropic characteristics.