• Geophysics and Geophysical Exploration
• /
• v.15 no.4
• /
• pp.227-234
• /
• 2012

In this study, we make simulation of density log using a Monte Carlo N-Particle (MCNP) algorithm to make an analysis on density logging under different borehole environments, since density logging is affected by various borehole conditions like borehole size, density of borehole fluid, thickness and type of casing, and so on. MCNP algorithm has been widely used for simulation of problems of nuclear particle transportation. In the simulation, we consider the specific configuration of a tool (Robertson Geologging Co. Ltd) that Korea institute of geoscience and mineral resources (KIGAM) has used. In order to measure accurate bulk density of a formation, it is essential to make a calibration and correction chart for the tool under considerations. Through numerical simulation, this study makes calibration plot of the density tool in material with several known bulk densities and with boreholes of several different diameters. In order to make correction charts for the density logging, we simulate and analyze measurements of density logging under different borehole conditions by considering borehole size, density of borehole fluid, and presence of casing.

• Geophysics and Geophysical Exploration
• /
• v.12 no.3
• /
• pp.255-262
• /
• 2009

Conventionally, vertical-seismic-profiling (VSP) survey that provides high-resolution information has mainly performed to obtain the exact depth of the gas hydrate-bearing sediment, which is one of the key factors in the development of the gas hydrate. In this study, we extracted interval velocities and created corridor stacks from the first domestic zero-offset VSP data, which were acquired with three component receivers at UBGH09 borehole in Ulleung Basin where gas hydrate exists. Then we compared the corridor stacks with a CMP stacked section from surface seismic data. First of all, we converted the signals recorded with three component receivers to true vertical and horizontal components by phase rotation, and divided the data into direct waves and reflected waves by wavefield separation processing. The trend of the interval velocity extracted from the zero-offset VSP was similar to that of the sonic log obtained at the same borehole. Because the interval velocity of the gas hydrate-bearing sediment above the BSR was high, and it decreased suddenly through the BSR, we could infer that free gas is accumulated below the BSR. The results of comparing the corridor stacks to the CMP stacked section of the surface seismic data showed that most reflection events agreed well with those in the surface CMP stacked section and that the phase-rotated VSP data corresponded better with the surface seismic data than the VSP data without phase rotation. In addition, by comparing a corridor stack produced from the transverse component with the CMP stacked section of the surface seismic data, we could identify PS mode-converted reflections in the CMP stacked section.

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

• Geophysics and Geophysical Exploration
• /
• v.11 no.2
• /
• pp.137-147
• /
• 2008

Fault rupture directivity of the Odaesan earthquake, which was inferred to be the main cause of the high PGAvalue (> 0.1 g) unusually observed at the near-source region, was analyzed by using the data from the nearby (R < 100 km) dense seismic stations. The Boatwright's method (2007) was adopted for this purpose in which the azimuth and takeoff angle of the unilateral rupture directivity function could be estimated based on the relative peak ground-motions of seismic stations resulting from the nature of the rupture directivity. In this study, the approximate values of the relative peak ground-motions was derived from the difference between the log residuals of the point-source spectral model (Boore, 2003) for the main and secondary events based on the Random Vibration Theory. In this derivation, the spectral difference for a frequency range between the source corner frequencies of main and secondary events was considered to reflect only the effect of the fault directivity. The inversion result of the model parameters for the fault directivity function showed that the fault-plane of NWW-SEE direction dipping steeply to the North with high rupture velocity near upward in SE direction is responsible for the observed high level of ground-motion at the near-source region.

• Geophysics and Geophysical Exploration
• /
• v.19 no.4
• /
• pp.187-199
• /
• 2016

In order to delineate the soft ground distribution, an integrated geostatistical analysis was performed using the MASW (Multichannel Analysis of Surface Wave) which has the information of overall region and CPTu (Piezo Cone Penetration Test) which provides the direct information of the measuring point of the ground. MASW results were known to have close relationship with the ground stiffness. This correlation was confirmed through the comparison of MASW data obtained from two survey lines to the laboratory test with extracted soil samples. 3D physical property distribution in the study area was acquired by geostatistical integrated analysis with the data of tip resistance ($q_c$) and pore pressure (u) from the CPTu obtained at 6 points within the study area. The integrated analysis was conducted by applying the COSGSIM (Sequential Gaussian Co-Simulation) technology which can carry out the simulation in accordance with the spatial correlation between the MASW results and both tip resistance and pore pressure. Besides the locations of CPTu, borehole investigations were also conducted at two different positions. As a result, the N value of SPT and borehole log could be secured, so these data were used for the analysis of the geotechnical engineering accuracy of the integrated analysis result. For the verification of reliability of the 3D distribution of tip resistance and pore pressure secured through integrated analysis, the geotechnical information gained from the two drilling areas was compared, and the result showed extremely high correlation.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.42 no.1
• /
• pp.35-44
• /
• 2022

The V_S30 map is used as a key variable for site amplification in the ShakeMap, which predicts ground motion at any site. However, no V_S30 map considering Korean geology and geomorphology has been developed yet. To develop a proxy-based V_S30 map, we used 1,101 V_S profiles obtained from a geophysical survey and collected proxy layers of geological and topographical information for the Korean Peninsula. Then, V_S30 prediction models were developed using linear regression analysis for each geological age considering the distribution of V_S30. As a result, models depending on geomorphology were suggested per each geologic group, including Quaternary, Fill, Ocean, Mesozoic group and Precambrian. Resolution of map is doubled from that of V_S30 map by U.S. Geological Survey (USGS). Standard deviation of residual in natural log of proxy-based V_S30 map is 0.233, whereas standard deviation of slope-based USGS V_S30 map is 0.387. Therefore, the proxy-based V_S30 map developed in this study is expected to have less uncertainty and to contribute to predicting more accurately the ground motion amplitude.

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