• Geophysics and Geophysical Exploration
• /
• v.7 no.3
• /
• pp.191-196
• /
• 2004

Seismic waveform inversion can be solved by using the classical Gauss-Newton method, which needs to construct the huge Hessian by the directly computed Jacobian. The property of Hessian mainly depends upon a source and receiver aperture, a velocity model, an illumination Bone and a frequency content of source wavelet. In this paper, we try to invert the Marmousi seismic data by controlling the huge Hessian appearing in the Gauss-Newton method. Wemake the two kinds of he approximate Hessian. One is the banded Hessian and the other is the approximate Hessian with automatic gain function. One is that the 1st updated velocity model from the banded Hessian is nearly the same of the result from the full approximate Hessian. The other is that the stability using the automatic gain function is more improved than that without automatic gain control.

• Geophysics and Geophysical Exploration
• /
• v.2 no.4
• /
• pp.167-173
• /
• 1999

Since the resolution of the 2-D hole-to-hole seismic traveltime tomography is affected by the limited ray transmission angle, various methods were used to improve the resolution. Linear traveltime interpolation(LTI) ray tracing method was chosen for forward-modeling method. Inversion results using the LTI method were compared with those using the other ray tracing methods. As an inversion algorithm, SIRT method was used. In the iterative non-linear inversion method, the cost of ray tracing is quite expensive. To reduce the cost, each raypath was stored and the inversion was performed from this information. Using the proposed method, fast convergence was achieved. Inversion results are likely to be affected by the initial velocity guess, especially when the ray transmission angle was limited. To provide a good initial guess for the inversion, generalized regression neural network(GRNN) method was used. When the transmitted raypath angle is not limited or the geological model is very complex, the inversion results are not affected by initial velocity model very much. Since the raypath angles, however, are limited in most geophysical tomographic problems, the enhancement of resolution in tomography can be achieved by providing a proper initial velocity model by another inversion algorithm such as GRNN.

• Geophysics and Geophysical Exploration
• /
• v.7 no.1
• /
• pp.14-18
• /
• 2004

High levels of ambient noise and safety factors often limit the use of 'active-source' seismic methods for geotechnical investigations in urban environments. As an alternative, shear-wave velocity-depth profiles can be obtained by treating the background microtremor wave field as a stochastic process, rather than adopting the traditional approach of calculating velocity based on ray path geometry from a known source. A recent field test in Melbourne demonstrates the ability of the microtremor method, using only Rayleigh waves, to resolve a velocity inversion resulting from the presence of a hard, 12 m thick basalt flow overlying 25 m of softer alluvial sediments and weathered mudstone. Normally the presence of the weaker underlying sediments would lead to an ambiguous or incorrect interpretation with conventional seismic refraction methods. However, this layer of sediments is resolved by the microtremor method, and its inclusion is required in one-dimensional layered-earth modelling in order to reproduce the Rayleigh-wave coherency spectra computed from observed seismic noise records. Nearby borehole data provided both a guide for interpretation and a confirmation of the usefulness of the passive Rayleigh-wave microtremor method. Sensitivity analyses of resolvable modelling parameters demonstrate that estimates of shear velocities and layer thicknesses are accurate to within approximately $10\%\;to\;20\%$ using the spatial autocorrelation (SPAC) technique. Improved accuracy can be obtained by constraining shear velocities and/or layer thicknesses using independent site knowledge. Although there exists potential for ambiguity due to velocity-thickness equivalence, the microtremor method has significant potential as a site investigation tool in situations where the use of traditional seismic methods is limited.

• Economic and Environmental Geology
• /
• v.49 no.3
• /
• pp.225-242
• /
• 2016

Rock physics serves as a useful tool for seismic reservoir characterization and monitoring by providing quantitative relationships between rock properties and seismic data. Rock physics models can predict effective moduli for reservoirs with different mineral components and pore fluids from well-log data. The distribution of reservoirs and fluids for the entire seismic volume can also be estimated from rock physics models. The first part of this report discusses the Voigt, Reuss, and Hashin-Shtrikman bounds for effective elastic moduli and the Gassmann fluid substitution. The second part reviews various contact models for moderate- to high-porosity sands. In the third part, constant-cement model, known to work well for the sand that gradually loses porosity with deteriorating sorting, was applied to the well-log data from an oil field in the North Sea. Lastly, the rock physics template constructed from the constant-cement model and the results from the prestack inversion of 2D seismic data were combined to predict the lithology and fluid types for the sand reservoir of this oil field.

• Geophysics and Geophysical Exploration
• /
• v.5 no.1
• /
• pp.33-45
• /
• 2002

It was necessary to devise new techniques to overcome and enhance the resolution limits of traveltime tomography. Waveform inversion has been one of the methods for giving very high resolution result. High resolution image could be acquired because waveform inversion used not only phase but amplitude. But waveform inversion was much time consuming Job because forward and backward modeling was needed at each iteration step. Velocity-stress method was used for effective modeling. Resolution limits of imaging methods such as travel time inversion, acoustic and elastic waveform inversion were investigated with numerical models. it was investigated that Resolution limit of waveform inversion was similar tn resolution limit of migration derived by Schuster. Horizontal resolution limit could be improved with increased coverage by adding VSP data in cross hole that had insufficient coverage. Also, waveform inversion was applied to realistic models to evaluate applicability and using initial guess of travel time tomograms to reduce non-linearity of waveform inversion showed that the better reconstructed image could be acquired.

• Journal of the Korean earth science society
• /
• v.21 no.6
• /
• pp.675-682
• /
• 2000

25 seismic shot gathers were obtained to study the two dimensional subsurface shear wave velocities in a landfilled area near the Keum river estuary. Borehole(BH#1 and BH#2) tests at two sites were made in the same area. Standard Penetration Tests were also performed at the same time. The 2-D shear wave velocity structure resulted from the inversion of the seismic data shows that the subsurface of the studied area consists of the upper 1${\sim}$3 meter thick layer(200 m/sec${\sim}$700 m/sec), the middle 5${\sim}$8 m thick low velocity layer(100 m/sec${\sim}$400 m/sec), and the lower layer of 1000m/sec or higher shear wave velocities. The thickness of the low velocity layer decreases from the BH #1 site to the BH #2 site. The depth to the basement also decreases toward the BH #2 site. The examination of the S wave velocity structure, the description of the geologic contents, and the Standard Penetration Test values indicate that the middle layer of low shear wave velocity may be related to the clay content of the layer. On the other hand, the Standard Penetration test values increase with depth, showing no significant relationship with the geologic contents of the subsurface. This study shows that the inversion of surface waves can be effective in the study of the shear wave velocity, especially in the area where low velocity layers can be found. The method of inversion of surface waves also can be used as a viable technique to overcome the limit of the seismic refraction method.

• Geophysics and Geophysical Exploration
• /
• v.8 no.2
• /
• pp.163-169
• /
• 2005

Recently, the velocity stack domain is having an attention as a very useful domain for various processing in seismic data processing. In order to be used in many applications, the velocity stack should be obtained through an inversion method and the used inversion should have properties like the robustness to noise and the parsimony of velocity stack result. Iteratively Reweighted Least-Squares (IRLS) method is the one of the inversion methods that have such properties. This paper describes the theoretical background, implementation of the method, and examines the characteristics and limits of the IRLS method.

• 한국지구물리탐사학회:학술대회논문집
• /
• 2007.12a
• /
• pp.85-85
• /
• 2007

• Economic and Environmental Geology
• /
• v.29 no.5
• /
• pp.629-637
• /
• 1996

In this paper a numerical experiment is conducted to determine the low acoustic impedance of a thin oil or gas reservoir from a seismogram by using the generalized linear inversion method. The seismograms used are normal incident synthetic seismograms containing p-wave primary reflections, multiples, and peg-leg multiples on the layers consisting of oil-, gas-, water-filled sandstone incased in shales. In this experiment the acoustic impedance, the location of reservoir boundary, thickness, and source wavelet are assumed initially and revised iteratively by the least-squares-error technique until the difference between the seismogram and calculated one is very small. This experiment shows that the acoustic impedance and thickness, about 10 m thick, can be determined by the inversion.

• Proceedings of the Earthquake Engineering Society of Korea Conference
• /
• 2003.03a
• /
• pp.45-49
• /
• 2003

The velocity structure beneath the CHNB broadband station is determined by receiver function analysis using by from teleseismic P waveforms. The detailed broadband receiver functions are obtained by stacking method for source-equalized vertical, radial and tangential components of teleseismic P waveforms. A time domain inversion uses the stacked radial receiver function to determine vertical P wave velocity structure beneath the station. The crustal velocity structures beneath the stations are estimated using the receiver function inversion method in the case at the crustal model parameterized by many thin, flat-tying, homogeneous layers. The result of crust at model inversion shows the crustal velocity structure beneath the CHNB station varies smoothly with increasing depth, and there are six discontinuity around 2.5km, 6.25km, 12.5km, 22.5km and 27.5km depth, with Moho discontinuity at about 32.5km depth.