• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.35 no.8
• /
• pp.957-964
• /
• 2011

This investigation presents a finite element method to obtain the transmission properties of bulk elastic waves in piezoelectric band gap structures(phonon crystals) for varying frequencies and modes. To this end, periodic boundary conditions are imposed on a three-dimensional model while both in-plane and out-of-plane modes are included. In particular, the mode decoupling characteristics between in-plane and out-of-plane modes are identified for each electric poling direction and the results are incorporated in the finite element modeling. Through numerical simulations, the proposed modeling method was found to be a useful, effective one for analyzing the wave characteristics of various types of piezoelectric phononic band gap structures.

• 한국신재생에너지학회:학술대회논문집
• /
• 2007.06a
• /
• pp.489-492
• /
• 2007

The P-wave velocity at the formation which contains gas hydrate varies very wide upon gas hydrate existence. These features on seismic shot gather can not be simulated normally by numerical modeling of homogeneous medium so that we need that of random inhomogeneous medium instead. We, in this study generated random inhomogeneous medium using gaussian ACF, exponential ACF and von Karman ACF and that we supposed the random inhomogeneous medium be gas hydrate formation to execute numeric modeling. The modeling result shows the typical effect by scattering caused by random hydrate formation as is observed from seismic shot gather where hydrate exist.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.20 no.6
• /
• pp.523-531
• /
• 2008

In hole type CPTu equipment which combines the concepts of inhole test method and piezocone test method was newly developed in order to evaluate the dynamic properties of marine soils. It is possible to perform inhole type CPTu without any additional source device because the source and receiver are contained inside the cone rod, which is different from the conventional seismic cone system. In this study, laboratory tests using kaolinite as soft soil and numerical simulations using finite element method were carried out to investigate the effects of several parameters including test methods and soil conditions on the test results from inhole type CPTu and to find out the optimum test method. It was found that it is necessary to maintain the length of swing arm as well as the distance between source and receiver consistently to obtain the rigorous test results. The laboratory test and numerical results also reveal that contrary to the input wave frequency, the water content of soil layer and the disturbance due to the installation of swing arm apparently affect the shear wave velocity.

• 한국지구과학회:학술대회논문집
• /
• 2002.09a
• /
• pp.50-51
• /
• 2002

• Geophysics and Geophysical Exploration
• /
• v.13 no.1
• /
• pp.80-87
• /
• 2010

This study was conducted to develop a reservoir modelling workflow to reproduce the heterogeneous distribution of effective permeability that impacts on the performance of SAGD (Steam Assisted Gravity Drainage), the in-situ bitumen recovery technique in the Athabasca Oil Sands. Lithologic facies distribution is the main cause of the heterogeneity in bitumen reservoirs in the study area. The target formation consists of sand with mudstone facies in a fluvial-to-estuary channel system, where the mudstone interrupts fluid flow and reduces effective permeability. In this study, the lithologic facies is classified into three classes having different characteristics of effective permeability, depending on the shapes of mudstones. The reservoir modelling workflow of this study consists of two main modules; facies modelling and permeability modelling. The facies modelling provides an identification of the three lithologic facies, using a stochastic approach, which mainly control the effective permeability. The permeability modelling populates mudstone volume fraction first, then transforms it into effective permeability. A series of flow simulations applied to mini-models of the lithologic facies obtains the transformation functions of the mudstone volume fraction into the effective permeability. Seismic data contribute to the facies modelling via providing prior probability of facies, which is incorporated in the facies models by geostatistical techniques. In particular, this study employs a probabilistic neural network utilising multiple seismic attributes in facies prediction that improves the prior probability of facies. The result of using the improved prior probability in facies modelling is compared to the conventional method using a single seismic attribute to demonstrate the improvement in the facies discrimination. Using P-wave velocity in combination with density in the multiple seismic attributes is the essence of the improved facies discrimination. This paper also discusses sand matrix porosity that makes P-wave velocity differ between the different facies in the study area, where the sand matrix porosity is uniquely evaluated using log-derived porosity, P-wave velocity and photographically-predicted mudstone volume.

• Geophysics and Geophysical Exploration
• /
• v.13 no.1
• /
• pp.24-30
• /
• 2010

This study utilises repeated numerical tests to understand the effects of variable near-surface conditions on time-lapse seismic surveys. The numerical tests were aimed at reproducing the significant scattering observed in field experiments conducted at the Naylor site in the Otway Basin for the purpose of $CO_2$ sequestration. In particular, the variation of elastic properties of both the top soil and the deeper rugose clay/limestone interface as a function of varying water saturation were investigated. Such tests simulate the measurements conducted in dry and wet seasons and to evaluate the contribution of these seasonal variations to seismic measurements in terms of non-repeatability. Full elastic pre-stack modelling experiments were carried out to quantify these effects and evaluate their individual contributions. The results show that the relatively simple scattering effects of the corrugated near-surface clay/limestone interface can have a profound effect on time-lapse surveys. The experiments also show that the changes in top soil saturation could potentially affect seismic signature even more than the corrugated deeper surface. Overall agreement between numerically predicted and in situ measured normalised root-mean-square (NRMS) differences between repeated (time-lapse) 2D seismic surveys warrant further investigation. Future field studies will include in situ measurements of the elastic properties of the weathered zone through the use of 'micro Vertical Seismic Profiling (VSP)' arrays and very dense refraction surveys. The results of this work may impact on other areas not associated with $CO_2$ sequestration, such as imaging oil production over areas where producing fields suffer from a karstic topography, such as in the Middle East and Australia.

• Geophysics and Geophysical Exploration
• /
• v.23 no.1
• /
• pp.38-49
• /
• 2020

Full-waveform inversion (FWI) is an optimization process of fitting observed and modeled data to reconstruct high-resolution subsurface physical models. In acoustic FWI (AFWI), pressure data acquired using a marine streamer has mainly been used to reconstruct the subsurface P-wave velocity models. With recent advances in marine seismic-acquisition techniques, acquiring multi-component data in marine environments have become increasingly common. Thus, AFWI strategies must be developed to effectively use marine multi-component data. Herein, we proposed an AFWI strategy using horizontal and vertical particle-acceleration data. By analyzing the modeled acoustic data and conducting sensitivity kernel analysis, we first investigated the characteristics of each data component using AFWI. Common-shot gathers show that direct, diving, and reflection waves appearing in the pressure data are separated in each component of the particle-acceleration data. Sensitivity kernel analyses show that the horizontal particle-acceleration wavefields typically contribute to the recovery of the long-wavelength structures in the shallow part of the model, and the vertical particle-acceleration wavefields are generally required to reconstruct long- and short-wavelength structures in the deep parts and over the whole area of a given model. Finally, we present a sequential-inversion strategy for using the particle-acceleration wavefields. We believe that this approach can be used to reconstruct a reasonable P-wave velocity model, even when the pressure data is not available.

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

• Geophysics and Geophysical Exploration
• /
• v.12 no.2
• /
• pp.173-182
• /
• 2009

In April 2008, KIGAM carried out an ocean-bottom seismometer (OBS) survey in the central Ulleung Basin where strong bottom simulating reflectors (BSRs) were revealed from previous surveys and some gas-hydrate samples were retrieved by direct sampling. The purpose of this survey is to estimate the velocity structure near the BSR in the gas hydrate prospect area using wide-angle seismic data recorded on the ocean-bottom seismometers. Along with the OBS survey, a 2-D seismic survey was performed whereby stratigraphic and preliminary velocity information was obtained. Two methods were applied to wide-angle data for estimating P wave velocity; one is velocity analysis in the $\tau$-p domain and the other is seismic traveltime inversion. A 1-D interval velocity profile was obtained by the first method, which was refined to layered velocity structure by the latter method. A layer stripping method was adopted for modeling and inversion. All velocity profiles at each OBS site clearly show velocity reversal at BSR depths due to the presence of gas hydrates. In addition, we could confirm high velocity in the column/chimney structure.

• Geophysics and Geophysical Exploration
• /
• v.14 no.1
• /
• pp.98-104
• /
• 2011

Numerical simulation in exploration geophysics provides important insights into subsurface wave propagation phenomena. Although elastic wave simulations take longer to compute than acoustic simulations, an elastic simulator can construct more realistic wavefields including shear components. Therefore, it is suitable for exploration of the responses of elastic bodies. To overcome the long duration of the calculations, we use a Graphic Processing Unit (GPU) to accelerate the elastic wave simulation. Because a GPU has many processors and a wide memory bandwidth, we can use it in a parallelised computing architecture. The GPU board used in this study is an NVIDIA Tesla C1060, which has 240 processors and a 102 GB/s memory bandwidth. Despite the availability of a parallel computing architecture (CUDA), developed by NVIDIA, we must optimise the usage of the different types of memory on the GPU device, and the sequence of calculations, to obtain a significant speedup of the computation. In this study, we simulate two- (2D) and threedimensional (3D) elastic wave propagation using the Finite-Difference Time-Domain (FDTD) method on GPUs. In the wave propagation simulation, we adopt the staggered-grid method, which is one of the conventional FD schemes, since this method can achieve sufficient accuracy for use in numerical modelling in geophysics. Our simulator optimises the usage of memory on the GPU device to reduce data access times, and uses faster memory as much as possible. This is a key factor in GPU computing. By using one GPU device and optimising its memory usage, we improved the computation time by more than 14 times in the 2D simulation, and over six times in the 3D simulation, compared with one CPU. Furthermore, by using three GPUs, we succeeded in accelerating the 3D simulation 10 times.