• Geophysics and Geophysical Exploration
• /
• v.9 no.3
• /
• pp.207-224
• /
• 2006

It has been widely known that the seismic piezo-cone penetration test (SCPTu) is one of the most useful techniques for investigating the geotechnical characteristics such as static and dynamic soil properties. As practical applications in Korea, SCPTu was carried out at two sites in Busan and four sites in Incheon, which are mainly composed of alluvial or marine soil deposits. From the SCPTu waveform data obtained from the testing sites, the first arrival times of shear waves and the corresponding time differences with depth were determined using the cross-over method, and the shear wave velocity $(V_S)$ profiles with depth were derived based on the refracted ray path method based on Snell's law. Comparing the determined $V_S$ profile with the cone tip resistance $(q_t)$ profile, both trends of profiles with depth were similar. For the application of the conventional CPTu to earthquake engineering practices, the correlations between $V_S$ and CPTu data were deduced based on the SCPTu results. For the empirical evaluation of $V_S$ for all soils together with clays and sands which are classified unambiguously in this study by the soil behavior type classification index $(I_C)$, the authors suggested the $V_S-CPTu$ data correlations expressed as a function of four parameters, $q_t,\;f_s,\;\sigma'_{v0}$ and $B_q$, determined by multiple statistical regression modeling. Despite the incompatible strain levels of the downhole seismic test during SCPTu and the conventional CPTu, it is shown that the $V_S-CPTu$ data correlations for all soils, clays and sands suggested in this study is applicable to the preliminary estimation of $V_S$ for the soil deposits at a part in Korea and is more reliable than the previous correlations proposed by other researchers.

• Journal of the Korean earth science society
• /
• v.37 no.5
• /
• pp.277-285
• /
• 2016

We introduce a depth scaling strategy to improve the accuracy of frequency-domain elastic full waveform inversion (FWI) using the new pseudo-Hessian matrix for seismic data without low-frequency components. The depth scaling strategy is based on the fact that the damping factor in the Levenberg-Marquardt method controls the energy concentration in the gradient. In other words, a large damping factor makes the Levenberg-Marquardt method similar to the steepest-descent method, by which shallow structures are mainly recovered. With a small damping factor, the Levenberg-Marquardt method becomes similar to the Gauss-Newton methods by which we can resolve deep structures as well as shallow structures. In our depth scaling strategy, a large damping factor is used in the early stage and then decreases automatically with the trend of error as the iteration goes on. With the depth scaling strategy, we can gradually move the parameter-searching region from shallow to deep parts. This flexible damping factor plays a role in retarding the model parameter update for shallow parts and mainly inverting deeper parts in the later stage of inversion. By doing so, we can improve deep parts in inversion results. The depth scaling strategy is applied to synthetic data without lowfrequency components for a modified version of the SEG/EAGE overthrust model. Numerical examples show that the flexible damping factor yields better results than the constant damping factor when reliable low-frequency components are missing.

• Geophysics and Geophysical Exploration
• /
• v.16 no.4
• /
• pp.225-233
• /
• 2013

The accuracy of the automatic NMO velocity analysis, which is used for an effective and objective NMO velocity analysis, is highly affected by the velocity resolution of the velocity spectrum. In this study, we have developed an automatic NMO velocity algorithm, where the velocity spectra are created using high-resolution bootstrapped differential semblance (BDS), and the velocity analysis on CMP gathers is performed in parallel with MPI. We also compared the velocity models from the developed automatic NMO velocity algorithm with high-resolution BDS to those from BDS. To verify the developed automatic velocity analysis module we created synthetic seismic data from a velocity model including horizon layers. We confirmed that the developed automatic velocity analysis module estimated velocity more accurately. In addition, NMO velocity which yielded a CMP stacked section, where the coherency of the events were improved, was estimated when the developed module was applied to a marine field data set.

• Geophysics and Geophysical Exploration
• /
• v.9 no.1
• /
• pp.93-101
• /
• 2006

In this study, we propose a joint inversion method, using genetic algorithms, to determine the shear-wave velocity structure of deep sedimentary layers from receiver functions and surface-wave phase velocity. Numerical experiments with synthetic data indicate that the proposed method can avoid the trade-off between shear-wave velocity and thickness that arises when inverting the receiver function only, and the uncertainty in deep structure from surface-wave phase velocity inversion alone. We apply the method to receiver functions obtained from earthquake records with epicentral distances of about 100 km, and Rayleigh-wave phase velocities obtained from a microtremor array survey in the Kanto Plain, Japan. The estimated subsurface structure is in good agreement with the previous results of seismic refraction surveys and deep borehole data.

• Geophysics and Geophysical Exploration
• /
• v.6 no.2
• /
• pp.95-100
• /
• 2003

In this study, we developed reflection tomography inversion algorithm using Straight Ray Technique (SRT) which can calculate travel time easily and fast for complex geological structure. The inversion process begins by setting the initial velocity model as a constant velocity model that hat only impedance boundaries. The inversion process searches a layer-interface structure model that is able to explain the given data satisfactorily by inverting to minimize data misfit. For getting optimal solution, we used Gauss-Newton method that needed constructing the approximate Hessian matrix. We also applied the Marquart-Levenberg regularization method to this inversion process to prevent solution diverging. The ability of the method to resolve typical target structures was tested in a synthetic salt dome inversion. Using the inverted velocity model, we obtained the migration image close to that of the true velocity model.

• Geophysics and Geophysical Exploration
• /
• v.20 no.4
• /
• pp.232-240
• /
• 2017

Normal moveout correction is one of the main procedures of seismic reflection data processing and a crucial pre-processing step for AVO analysis. Unfortunately, stretch phenomenon, which is the intrinsic problem of NMO correction, degrades the quality of stack section and reliability of AVO analysis. Although muting is applied to resolve this problem, it makes far-offset traces more useful to develop an advanced NMO correction technique without stretch. In this paper, easy and detailed explanations are provided on the definition and methodology of NMO correction, and then the cause of stretch is explained with its characteristics. A graphical explanation for NMO correction is given for the intuitive understanding of stretch phenomenon. Additionally, the theoretical formulation is derived to quantitatively understand the NMO correction. Through explaining the muting process to remove NMO stretch, the limitations of conventional methods are investigated and the need for a new resolution comes to discussion. We describe a stretch-free NMO correction based on inverse theory among many different stretch-free NMO corrections. Finally, the stretch-free NMO correction is verified through synthetic example and real data.

• Geophysics and Geophysical Exploration
• /
• v.11 no.3
• /
• pp.242-249
• /
• 2008

Recently, AVO analysis has been widely used in oil exploration with seismic subsurface section as a direct indicator of the existence of the gas. In the case of the deep reservoirs like the gas reservoirs in the East-sea, it is often difficult to observe AVO responses in CMP gathers even though the bright spots are shown in the stacked section. Because the reservoir becomes more consolidated as its depth deepens, P-wave velocity does not decrease significantly when the pore fluid is replaced by the gas. Thus the difference in Poisson's ratio, which is a key factor for AVO response, between the reservoir and the layer above it does not increase significantly. In this study, we analyzed the effects of Poisson's ratio difference on AVO response with a variety of Poisson's ratios for the upper and lower layers. The results show that, as the difference in Poisson's ratio between the upper and lower layers decreases, the change in the reflection amplitude with incidence angle decreases and AVO responses become insignificant. To consider the limitation of AVO responses shown in the gas reservoir in East-sea, the velocity model was made by simulation Gorae V structure with seismic data and well logs. The results of comparing AVO responses observed from the synthetic data with theoretical AVO responses calculated by using material properties show that the amount of the change in reflection amplitude with increasing incident angle is very small when the difference in Poisson's ratio between the upper and lower layers is small. In addition, the characteristics of AVO responses were concealed by noise or amplitude distortion arisen during preprocessing. To overcome such limitations of AVO analysis of the data from deep reservoirs, we need to acquire precisely reflection amplltudes In data acquisition stage and use processing tools which preserve reflection amplitude in data processing stage.

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

• Geophysics and Geophysical Exploration
• /
• v.10 no.1
• /
• pp.19-28
• /
• 2007

In this study, we propose a joint inversion method, using genetic algorithms, to estimate an S-wave velocity structure for deep sedimentary layers from receiver functions and surface-wave phase velocity observed at several sites. The method takes layer continuity over a target area into consideration by assuming that each layer has uniform physical properties, especially an S-wave velocity, at all the sites in a target area in order to invert datasets acquired at different sites simultaneously. Numerical experiments with synthetic data indicate that the proposed method is effective in reducing uncertainty in deep structure parameters when modelling only surface-wave dispersion data over a limited period range. We then apply the method to receiver functions derived from earthquake records at one site and two datasets of Rayleigh-wave phase velocity obtained from microtremor array surveys performed in central Tokyo, Japan. The estimated subsurface structure is in good agreement with the results of previous seismic refraction surveys and deep borehole data. We also conclude that the proposed method can provide a more accurate and reliable model than individual inversions of either receiver function data only or surface-wave dispersion data only.

• Geophysics and Geophysical Exploration
• /
• v.9 no.1
• /
• pp.60-66
• /
• 2006

Three differing sandstones, two synthetic and one field sample, have been tested ultrasonically under a range of confining pressures and pore pressures representative of in-situ reservoir pressures. These sandstones include: a synthetic sandstone with calcite intergranular cement produced using the CSIRO Calcite In-situ Precipitation Process (CIPS); a synthetic sandstone with silica intergranular cement; and a core sample from the Otway Basin Waarre Formation, Boggy Creek 1 well, from the target lithology for a trial $CO_2$ pilot project. Initial testing was carried on the cores at "room-dried" conditions, with confining pressures up to 65 MPa in steps of 5 MPa. All cores were then flooded with $CO_2$, initially in the gas phase at 6 MPa, $22^{\circ}C$, then with liquid-phase $CO_2$ at a temperature of $22^{\circ}C$ and pressures from 7 MPa to 17 MPa in steps of 5 MPa. Confining pressures varied from 10 MPa to 65 MPa. Ultrasonic waveforms for both P- and S-waves were recorded at each effective pressure increment. Velocity versus effective pressure responses were calculated from the experimental data for both P- and S-waves. Attenuations $(1/Q_p)$ were calculated from the waveform data using spectral ratio methods. Theoretical calculations of velocity as a function of effective pressure for each sandstone were made using the $CO_2$ pressure-density and $CO_2$ bulk modulus-pressure phase diagrams and Gassmann effective medium theory. Flooding the cores with gaseous phase $CO_2$ produced negligible change in velocity-effective stress relationships compared to the dry state (air saturated). Flooding with liquid-phase $CO_2$ at various pore pressures lowered velocities by approximately 8% on average compared to the air-saturated state. Attenuations increased with liquid-phase $CO_2$ flooding compared to the air-saturated case. Experimental data agreed with the Gassmann calculations at high effective pressures. The "critical" effective pressure, at which agreement with theory occurred, varied with sandstone type. Discrepancies are thought to be due to differing micro-crack populations in the microstructure of each sandstone type. The agreement with theory at high effective pressures is significant and gives some confidence in predicting seismic behaviour under field conditions when $CO_2$ is injected.