• Economic and Environmental Geology
• /
• v.46 no.6
• /
• pp.485-494
• /
• 2013

$CO_2$ Monitoring, Mitigation and Verification (MMV) is the essential part in the Carbon Capture and Storage (CCS) project in order to assure the storage permanence economically and environmentally. In large-scale CCS projects in the world, the seismic time-lapse survey is a key technology for monitoring the behavior of injected $CO_2$. In this study, we developed a basic process procedure for 3-D seismic baseline data from the Aquistore project, Estevan, Canada. Major target formations of Aquistore CCS project are the Winnipeg and the Deadwood sandstone formations located between 1,800 and 1,900 ms in traveltime. The analysis of trace energy and similarity attributes of seismic data followed by spectral decomposition are carried out for the characterization of $CO_2$ injection zone. High trace energies are concentrated in the northern part of the survey area at 1,800 ms and in the southern part at 1,850 ms in traveltime. The sandstone dominant regions are well recognized with high reflectivity by the trace energy analysis. Similarity attributes show two structural discontinuities trending the NW-SE direction at the target depth. Spectral decomposition of 5, 20 and 40 Hz frequency contents discriminated the successive E-W depositional events at the center of the research area. Additional noise rejection and stratigraphic interpretation on the baseline data followed by applying appropriate imaging technique will be helpful to investigate the differences between baseline data and multi-vintage monitor data.

• Geophysics and Geophysical Exploration
• /
• v.16 no.3
• /
• pp.163-170
• /
• 2013

The amplitude variation with offset of seismic data can detect fluids in the sediment and resolve the petrophysical properties of hydrocarbons in the subsurface. We analyzed and described the amplitude variation in deep sea seismic data obtained from the Ulleung Basin, East Sea. By inspecting seismic CDP-offset and CDP-angle gathers which show a bright reflection event, we decided a target zone for amplitude variation analysis. From the seismic angle gather at the middle of Ulleung Basin, we recognized amplitude increase or decrease versus offset on the intercept-gradient curve. Using the product attribute and Poisson's ratio change attribute computed in terms of intercept with gradient, the top and the base of gas saturated sediments were described. The area of amplitude variation suggestive of the presence of gas saturated sediments is shown at the depth of 3 s traveltime. Anomalous features of seismic amplitude in the Ulleung Basin were classified by the crossplot of intercept and gradient. The background trend of crossplot between intercept and gradient shows an inverse proportional relation that is common for wet sediments. Anomalous amplitudes of Class III fall into the first and the third quadrants on crossplots. We inferred regional gas/water saturated area with the horizontal dimension of 150 m in the Ulleung Basin by cross-section with respect to cross-plot anomaly.

• 한국지구물리탐사학회:학술대회논문집
• /
• 2006.06a
• /
• pp.125-130
• /
• 2006

Both surface seismic and far-offset VSP data were recorded alongtwo mutually perpendicular profiles in the Pungam basin. The first-arrival times were simultaneously inverted using the tomography method. For the surface data, seismic energy was generated by a 5-kg sledgehammer at 48 stations and detected by 21 surface geophones at 3 m intervals and one 3-component geophone in test borehole for the purpose of static corrections. For the VSP data, seismic waves generated by the sledgehammer on the ground were detected by a 3-component borehole geophone in a depth range of $9{\sim}99\;m$. Delay times of the hammer data were corrected using the seisgun data before the inversion to yield velocity tomograms. The tomograms indicates that the soil layer with velocities less than 750 m/s averages 1.8 m thick. The velocity varies from 5353 m/s at the depth range of $31{\sim}40\;m$ to 4262 m/s at the depth range of $65{\sim}73\;m$. Compared with core samples, the relatively large variation in velocity may due to lithology changes and fracture effects with depth.

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

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

• Economic and Environmental Geology
• /
• v.47 no.6
• /
• pp.625-633
• /
• 2014

The Aquistore project is the world's first commercial capture, transportation, utilization and storage project of post-combustion $CO_2$ from a coal-fired thermo electric power plant, and the proposed storage is a saline aquifer at a depth of about 3,500 m. Deep saline aquifer, compared to hydrocarbon reservoir, provides the great volumetric potential for storage of $CO_2$ anywhere in the world, therefore the research results from the project may be exported globally to other sites. Geological $CO_2$ storage characterization for saline aquifer instead of hydrocarbon reservoir needs to estimate the geophysical properties of subsurface geology. This study calculated the geophysical property of water-saturated formation by applying amplitude variation analysis developed from oil and gas exploration. We correlated horizon tops at the well logs to seismic traveltime of 1,815 and 1,857 ms as Winnipeg and Deadwood formations. Gradient analysis from seismic traces showed correlation coefficient of 45 - 81 % on amplitude variation with respect to incident angle. Crossplot of intercept and gradient shows the inverse proportional trend which represents typical water saturated sediments. Product attribute of intercept and gradient described the base of wet sediment. Poisson's ratio change attribute increased at the top of target area satisfying with wet sediment and decreased at the top of basement in a dry rock bed.

• Geophysics and Geophysical Exploration
• /
• v.13 no.2
• /
• pp.153-158
• /
• 2010

To investigate the velocity structure in the central and southern parts of the Korean peninsula, a 299-km NW-SE seismic refraction profile KCRT-2008was obtained across major tectonic boundaries. Seismic waves were generated by detonating 250 ~ 1500 kg explosives at depths of 50 ~ 100 m in eight drill holes located at intervals of 21 ~ 113 km. The seismic signals were detected by 4.5 Hz geophones at a nominal interval of 500 m. The first-arrival times were inverted to derive a velocity tomogram. The raypaths indicate several mid-crust interfaces including those at approximate depths of 2 ~ 3, 11 ~ 13, and 20 km. The Moho discontinuity with refraction velocity of 7.7 to 8.1 km/s has a maximum depth of 34.5 km under the central portion of the peninsula. The Moho becomes shallower as the Yellow Sea and the East Sea are approached on the west and east coasts of the peninsula, respectively. The depth of the 7.6 km/s velocity contour varies from 31.3 km to 34.4 km. The velocity tomogram shows the existence of a 129 km wide low-velocity zone centered at 7.2 km depth under the Okchon fold belt and Gyeonggi massif and low-velocity(< 5.4 km/s) rocks in the Gyeongsang sedimentary basin with a maximum thickness of 2.6 km

• Geophysics and Geophysical Exploration
• /
• v.24 no.2
• /
• pp.53-66
• /
• 2021

Velocity model building is an essential procedure in seismic data processing. Conventional techniques, such as traveltime tomography or velocity analysis take longer computational time to predict a single velocity model and the quality of the inversion results is highly dependent on human expertise. Full-waveform inversions also depend on an accurate initial model. Recently, deep neural network techniques are gaining widespread acceptance due to an increase in their integration to solving complex and nonlinear problems. This study investigated cases of seismic velocity model building using deep neural network techniques by classifying items according to the neural networks used in each study. We also included cases of generating training synthetic velocity models. Deep neural networks automatically optimize model parameters by training neural networks from large amounts of data. Thus, less human interaction is involved in the quality of the inversion results compared to that of conventional techniques and the computational cost of predicting a single velocity model after training is negligible. Additionally, unlike full-waveform inversions, the initial velocity model is not required. Several studies have demonstrated that deep neural network techniques achieve outstanding performance not only in computational cost but also in inversion results. Based on the research results, we analyzed and discussed the characteristics of deep neural network techniques for building velocity models.

• Geophysics and Geophysical Exploration
• /
• v.11 no.1
• /
• pp.1-14
• /
• 2008

We present a method for interpreting seismic records with arrivals and waveforms having characteristics which could be generated by extremely inhomogeneous velocity structures, such as non-typical oceanic crust, decollement at subduction zones, and seamounts in oceanic regions, by comparing them with synthetic waveforms. Recent extensive refraction and wide-angle reflection surveys in oceanic regions have provided us with a huge number of high-resolution and high-quality seismic records containing characteristic arrivals and waveforms, besides first arrivals and major reflected phases such as PmP. Some characteristic waveforms, with significant later reflected phases or anomalous amplitude decay with offset distance, are difficult to interpret using only a conventional interpretation method such as the traveltime tomographic inversion method. We find the best process for investigating such characteristic phases is to use an interactive interpretation method to compare observed data with synthetic waveforms, and calculate raypaths and traveltimes. This approach enables us to construct a reasonable structural model that includes all of the major characteristics of the observed waveforms. We present results here with some actual observed examples that might be of great help in the interpretation of such problematic phases. Our approach to the analysis of waveform characteristics is endorsed as an innovative method for constructing high-resolution and high-quality crustal structure models, not only in oceanic regions, but also in the continental regions.

• Economic and Environmental Geology
• /
• v.44 no.3
• /
• pp.229-238
• /
• 2011

S-wave, which provides lithology and pore fluid information, plays a key role in estimating gas-hydrate saturation. In general, P- and S-wave velocities increase in the presence of gas-hydrate and the P-wave velocity decreases in the presence of free gas under the gas-hydrate layer. Whereas there are very small changes, even slightly increases, in the S-wave velocity in the free gas layer because S-wave is not affected by the pore fluid when propagating in the free gas layer. To verify those velocity properties of the BSR (bottom-simulating reflector) depth in the gas-hydrate prospect area in the Ulleung Basin, P- and S-wave velocity profiles were derived from multi-component ocean-bottom seismic data which were acquired by Korea Institute of Geoscience and Mineral Resources (KIGAM) in May 2009. OBS (ocean-bottom seismometer) hydrophone component data were modeled and inverted first through the traveltime inversion method to derive P-wave velocity and depth model of survey area. 2-D multichannel stacked data were incorporated as an initial model. Two horizontal geophone component data, then, were polarization filtered and rotated to make radial component section. Traveltimes of main S-wave events were picked and used for forward modeling incorporating Poisson's ratio. This modeling provides S-wave profiles and Poisson's ratio profiles at every OBS site. The results shows that P-wave velocities in most OBS sites decrease beneath the BSR, whereas S-wave velocities slightly increase. Consequently, Poisson's ratio decreased strongly beneath the BSR indicating the presence of a free gas layer under the BSR.