• Geophysics and Geophysical Exploration
• /
• v.12 no.4
• /
• pp.347-354
• /
• 2009

Physical properties of rocks are strongly dependant on details of pore micro-structures, which can be used for quantifying relations between physical properties of rocks through pore-scale simulation techniques. Recently, high-resolution scan techniques, such as X-ray microtomography and high performance computers make it possible to calculate permeability from pore micro-structures of rocks. We try to extend this simulation methodology to velocity and electrical conductivity. However, the smoothing effect during tomographic inversion creates artifacts in pore micro-structures and causes inaccurate property estimation. To mitigate this artifact, we tried to use sharpening filter and neural network classification techniques. Both methods gave noticeable improvement in pore structure imaging and accurate estimation of permeability and electrical conductivity, which implies that our method effectively removes the smoothing effect in pore structures. However, the calculated velocities showed only incremental improvement. By comparison between thin section images and tomogram, we found that our resolution is not high enough, and it is mainly responsible for the inaccuracy in velocity despite the successful removal of the smoothing effect. In conclusion, our methods can be very useful for pore-scale modeling, since it can create accurate pore structure without the smoothing effect. For accurate velocity estimation, the resolution of pore structure should be at least three times higher than that for permeability simulation.

• Geophysics and Geophysical Exploration
• /
• v.12 no.4
• /
• pp.328-337
• /
• 2009

To investigate the relation between pore fluid conductivity and bulk resistivity of a rock sample it is assumed that electrolyte solution perfectly substitute the pore fluid that occupied the pore space within the sample in general. In this study, it is investigated that how much can the electrolyte solution substitute the pore fluid by repeating the same saturation process. Four kinds of NaCl solutions of 8, 160, 3200, 64000 ${\mu}S$/cm are used. The saturation process has repeated four times for each electrolyte in increasing conductivity order first then four times each in decreasing order. The more the saturation process repeated with the same electrolyte, the more electrolyte solution substitute the pore fluid. Geometric mean of bulk resistivity in increasing and decreasing orders with the same electrolyte solution is assumed to be mostly close to the bulk resistivity with perfect substitution. Bulk resistivity measurements for both increasing and decreasing order differs within 10% to the geometric mean when repeating the saturation process 4 times while maximum 40% difference is observed when single saturation process for each electrolyte solution with increasing order. The modified parallel resistant model can generally represent the relations between pore fluid resistivity and bulk resistivity in the experiment, but more experimental data with various rock samples with different porosity is needed to generalize the model.

• Geophysics and Geophysical Exploration
• /
• v.12 no.4
• /
• pp.316-327
• /
• 2009

A standard for national seismological observatory was proposed on 1999. Since then, Quanterra digitizer has been installed and is operating on almost all of seismic stations which belong to major seismic monitoring organizations. Quanterra digitizer produce and transmit real-time event packet and data packet. Characteristics of event packet and arrival time of each channel's data packet on data center were investigated. Packet selection criteria using signal to noise ratio (hereafter SNR) and signal period from real-time event packet based on 100 samples per second (hereafter sps) velocity data were developed. Estimation of epicenter using time information of the selected event packet were performed and tested. A series of experiment show that event packets were received approximately 3~4 second earlier than data packets and the number of event packet was only 0.3% compare to data packets. Just about 5% against all of event packets were selected as event packet were related P wave of real earthquake. Using the selected event packets we can estimate an epicenter with misfit less than 10 km within 20 sec for local earthquake over magnitude 2.5.

• Geophysics and Geophysical Exploration
• /
• v.12 no.4
• /
• pp.299-308
• /
• 2009

A magentotelluric (MT) survey at the Uiseong area has been performed for the site investigation of pilot scale $CO_2$ sequestration. The purpose of the MT survey is to delineate deeply extended fracture systems that can act as a leakage path of injected $CO_2$ Plume. Since the target area is extremely noisy in electromagentic sense, low frequency data below 1 Hz cannot be used for inversion. Two- and three-dimensional interpretation of the MT data showed a very clear conductive anomaly, which has the direction of $N55\sim65^{\circ}W$ and is extended roughly down to 1.6 km. It have the same direction with the strike-slip faults, the Gaeum and Geumcheon Faults. On the contrary, the eastern part of the survey area shows relatively homogeneous to the depth of 2 km though some small fractures at shallow depths can be found. Test drilling and high-definition borehole surveys should be followed at the eastern part of the survey area and hydraulic fracturing is required for injection of $CO_2$, because mean porosity of the sedimetary rock in the area is only 1.47%.

• Geophysics and Geophysical Exploration
• /
• v.12 no.4
• /
• pp.289-298
• /
• 2009

Optimal data processing parameters were designed to attenuate multiples in seismic data acquired in the south-eastern area of the East Sea, in 2008. Bunch of multiples caused by shallow sea water depth were perceived periodically up to two way travel time of 1,750 ms at every 250 ms over seismic traces. We abbreviated sea bottom multiple as SBM, Peg-leg multiple as PLM, and free-surface multiple as FSM. To attenuate these multiples, seismic data processing flow was constructed including NMO, stack, minimum phase predictive deconvolution filter and wave equation multiple rejections (WEMR). Prevalent multiples were suppressed by predictive deconvolution and remaining multiples were attenuated by WEMR. We concluded that combining deconvolution with WEMR was effective to a seismic data of study area. Derived parameter can be applied to the seismic data processing on adjacent survey area.

• Geophysics and Geophysical Exploration
• /
• v.17 no.4
• /
• pp.231-241
• /
• 2014

Recently, single channel seismic survey for engineering purpose have been used widely taking advantage of simple processing. However it is very difficult to obtain high fidelity subsurface image by single channel seismic due to insufficient fold coverage. Recently, seismic waveform inversion in multi channel seismic survey is utilized for accurate subsurface imaging even in complex terrains. In this paper, we propose the seismic waveform inversion algorithm for velocity model building using a single channel seismic data. We utilize the Gauss-Newton method and assume that subsurface model is 1-Dimensional. Seismic source estimation technique is used and offset effect is also corrected by removing delay time by offset. Proposed algorithm is verified by applying modified Marmousi2 model, and applied to field data set obtained in port of Busan.

• Geophysics and Geophysical Exploration
• /
• v.17 no.4
• /
• pp.216-230
• /
• 2014

Based on the interpretation of 3D seismic profiles acquired in the northwestern area of the Ulleung Basin, East Sea, the shallow sediments consist of five seismic units separated by regional reflectors. An anticline is present in the study area that documents activity of many faults. Bottom simulating reflectors are characterized by high RMS amplitude. Acoustic blanking with low RMS amplitude is distinctively recognized in the gas hydrate stability zone. Seismic attribute analysis shows that if gas hydrates are underlain by free gas, the high reflection strength and the low instantaneous frequency are displayed below the boundary between them. Whereas, if not, the reflection strength is low and instantaneous frequency is high continuously below the gas hydrate zone. Based on the spectral decomposition of the bottom simulating reflector, the high envelope at the specific high frequency range indicates the generation of the tuning effect due to the lower free gas content. Four models for the occurrence of the gas hydrate are suggested considering the slope of sedimentary layers as well as the presence of gas hydrate or free gas.

• Geophysics and Geophysical Exploration
• /
• v.17 no.4
• /
• pp.187-201
• /
• 2014

Compared with the separated inversion of electromagnetic (EM) and seismic data, a joint inversion using both EM and seismic data reduces the uncertainty and gives the opportunity to use the advantage of each data. Seismic fullwaveform inversion allows velocity information with high resolution in complicated subsurface. However, it is an indirect survey which finds the structure containing oil and gas. On the other hand, marine controlled-source EM (mCSEM) inversion can directly indicate the oil and gas using different EM properties of hydrocarbon with marine sediments and cap rocks whereas it has poor resolution than seismic method. In this paper, we have developed a joint EM inversion algorithm using a cross-gradient technique. P-wave velocity structure obtained by full-waveform inversion using plane wave encoding is used as structure constraints to calculate the cross-gradient term in the joint inversion. When the jointinversion algorithm is applied to the synthetic data which are simulated for subsea reservoir exploration, images have been significantly improved over those obtained from separate EM inversion. The results indicate that the developed joint inversion scheme can be applied for detecting reservoir and calculating the accurate oil and gas reserves.

• Geophysics and Geophysical Exploration
• /
• v.16 no.2
• /
• pp.79-90
• /
• 2013

Geophysical investigation in non-destructive testing is economically less expensive than boring testing and providing geotechnical information over wide-area. But, it provides only limited geotechnical information, which is hardly used to the design. Accordingly, we performed electrical resistivity experiments on large scale of soil model to analyze the correlation between electrical resistivity response and soil water contents. The soils used in the experiments were the Jumunjin standard sand and weathered granite soil. Each soil particle size distribution and coefficient of uniformity of experimental material obtained in the experiments were maintained in a state of the homogeneous. The specifications of the model used in this study is $160{\times}100{\times}50$(cm) of acrylic, and each soil was maintained at the height 30 cm. The water content were measured using the 5TE sensors (water contents sensors) which is installed 7 ~ 8 cm apart vertically by plugging to floor. The results of the resistivity behavior pattern for Jumunjin standard sand was found to be sensitive to the water content, while the weathered granite soil was showing lower resistivity over the time, and there was no significant change in behavior pattern observed. So, it results that the Jumunjin standard sand's particle current conduction was better than the weathered granite soil's particle through contact with the distilled water. This lab test was also compared with the result of a test bed site composed of similar weathered soil. It was confirmed that these experiments were underlying research of non-destructive investigation techniques to improve the accuracy to estimate the geotechnical parameter.

• Geophysics and Geophysical Exploration
• /
• v.16 no.2
• /
• pp.71-78
• /
• 2013

The one-way wave equation migration is much more computationally efficient comparing with reverse time migration and it can provide better image than the migration algorithm based on the ray theory. We have developed the prestack depth migration module adopting (GS) propagator designed for vertical transverse isotropic media. Since GS propagator considers the higher-order term by expanding the Taylor series of the vertical slowness in the thin slab of the phase-screen propagator, the GS migration can offer more correct image for the complex subsurface with large lateral velocity variation or steep dip. To verify the validity of the developed GS migration module, we analyzed the accuracy with the order of the GS propagator for VTI media (GSVTI propagator) and confirmed that the accuracy of the wavefield propagation with the wide angles increases as the order of the GS propagator increases. Using the synthetic seismic data, we compared the migration results obtained from the isotropic GS migration module with the anisotropic GS migration module. The results show that the anisotropic GS migration provides better images and the improvement is more evident on steeply dipping structures and in a strongly anisotropic medium.