• Geophysics and Geophysical Exploration
• /
• v.14 no.3
• /
• pp.191-202
• /
• 2011

Reserve-time migration (RTM) using a two-way wave equation is one of the most accurate migration techniques. RTM has been conducted by assuming that subsurface media are isotropic. However, anisotropic media are commonly encountered in reality. Conventional isotropic RTM may yield inaccurate results for anisotropic media. In this paper, we develop RTM algorithms for vertical transversely isotropic media (VTI) and tilted transversely isotropic media (TTI). For this, the pseudo-acoustic wave equations are used. The modeling algorithms are based on the high-order finite-difference method (FDM). The RTM algorithms are composed using the cross-correlation imaging condition or the imaging condition using virtual sources. By applying the developed RTM algorithms to the Hess VTI and BP TTI models, we could obtain better images than those obtained by the conventional isotropic RTM.

• Journal of the Korean Geophysical Society
• /
• v.2 no.2
• /
• pp.123-134
• /
• 1999

Magnetotelluric responses may be affected by strong anisotropy of the high-conductivity layers (HCL) in the upper mantle or lower crust. We have studied two-dimensional anisotropy MT modelling to examine the effect of high anisotropic media. Electrical properties of a homogeneous anisotropic body are defined by a symmetric conductivity tensor and the problem is described by coupled diffusion equation in the frequency domain. In two-dimensional anisotropic environments, diagonal elements of the impedance tensor have higher values than those in isotropic environments. In some cases, TM mode phases reach more than 90°and apparent resistivities decrease for some frequency range because of telluric distortion. GB decomposition may be used to recover regional responses, but can be affected by the regional anisotropic effect. Considering these results, BC87 dataset was interpreted with a modified anisotropic model.

• Journal of Broadcast Engineering
• /
• v.10 no.3
• /
• pp.383-391
• /
• 2005

In this paper, a new anisotropic diffusion based on Alvarez, Lions, and Morel (ALM) diffusion model is proposed for the suppression of blocking artifact caused by discrete cosine transform (DCT) based image compression. The proposed diffusion model, which incorporates a 'rate control parameter' (RCP), makes it possible to reduce blocking artifacts while to preserve the edge. The RCP controls the rate between isotropic and anisotropic diffusion. Isotropic diffusion is encouraged to eliminate the blocking artifacts in a block boundary of a smooth region, while anisotropic diffusion is encouraged to keep the edge or texture sharp in edge and a block boundary within an edge region. Additionally, to avoid oversmoothness of the texture region, a 'speed control parameter' (SCP), which makes diffusion process slow in the texture region, is employed.

• Proceedings of the Optical Society of Korea Conference
• /
• 2007.02a
• /
• pp.96-97
• /
• 2007

• Tunnel and Underground Space
• /
• v.13 no.6
• /
• pp.465-475
• /
• 2003

Groundwater inflow into the caverns constructed in fractured rock mass was simulated by numerical modeling, NAPSAC (DFN, discrete fracture network model) and NAMMU (CPM, continuous porous media model), a finite-element software package for groundwater flow in 3D fractured media developed by AEA Technology, UK. The input parameters for modeling were determined on surface fracture survey, core logging and single hole hydraulic test data. In order to predict the groundwater inflow more accurately, the anisotropic hydraulic conductivity was considered. The anisotropic hydraulic conductivities were calculated from the fracture network properties. With a minor adjustment during model calibration, the numerical modeling is able to reproduce reasonably groundwater inflows into cavern and the travel length and times to the ground surface along the flow paths in the normal, dry and rainy seasons.

• Korean Journal of Optics and Photonics
• /
• v.5 no.1
• /
• pp.25-30
• /
• 1994

We derive a set of formuale which show one-to-one correspondence between the the unitary Jones matrices of transparent anisotropic media and the rotational transformations on the Poincare sphere. By using the formuale one can determine the vectorial representation of the rotational transformation on the Poincare sphere which specifies the direction of the axis and the angle of the rotation in terms of the three parameters specific to the corresponding unitary Jones matrix, and conversely the the three parameters of the uniatry Jones matrix in terms of the vectorial representation of the corresponding rotational transformation on the Poincare sphere. To understand the polarization transmission characteristics of an optical system consisting of transparent linear anisotropic media, start with the Jones calculus to get the unitary Jones matrix for the whole system and then convert it to a rotational transformation on the Poincare sphere, from which we can intuitively understand the effect of the optical system on the polarization state of the light passing through the system.

• 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.

• 한국지구물리탐사학회:학술대회논문집
• /
• 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.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.5 no.1
• /
• pp.55-78
• /
• 2001

Given the anisotropic Poisson equation $-{\nabla}{\cdot}{\mathcal{K}}{\nabla}p=f$, one can convert it into a system of two first order PDEs: the Darcy law for the flux $u=-{\mathcal{K}{\nabla}p$ and conservation of mass ${\nabla}{\cdot}u=f$. A very natural mixed finite volume method for this system is to seek the pressure in the nonconforming P1 space and the Darcy velocity in the lowest order Raviart-Thomas space. The equations for these variables are obtained by integrating the two first order systems over the triangular volumes. In this paper we show that such a method is really a standard finite element method with local recovery of the flux in disguise. As a consequence, we compare two approaches in analyzing finite volume methods (FVM) and shed light on the proper way of analyzing non co-volume type of FVM. Numerical results for Dirichlet and Neumann problems are included.

• Structural Engineering and Mechanics
• /
• v.59 no.6
• /
• pp.951-972
• /
• 2016

The purpose of this study is to illustrate the propagation of the shear waves (SH-waves) in a prestressed hetrogeneous orthotropic media overlying a pre-stressed anisotropic porous half-space with self weight. It is considered that the compressive initial stress, mass density and moduli of rigidity of the upper layer are space dependent. The proposed model is solved to obtain the different dispersion relations for the SH-wave in the elastic-porous medium of different properties. The effects of compressive and tensile stresses along with the heterogeneity, porosity, Biot's gravity parameter on the dispersion of SH-wave are shown numerically. The wave analysis further indicates that the technical parameters of upper and lower half-space affect the wave velocity significantly. The results may be useful to understand the nature of seismic wave propagation in geophysical applications and in the field of earthquake and material science engineering.