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

• The Journal of the Acoustical Society of Korea
• /
• v.24 no.6
• /
• pp.334-342
• /
• 2005

In this Paper. a time-domain geoacoustic inversion was performed using the bulb signals measured during MがU. 04 experiment conducted in the East Sea of Korea in 2004. An obiective function was defined as a direct cross-correlation between the measured and the simulated signals in time domain. The ray theory was used to model the wave propagation in time domain and optimizations were Performed using VFSA (very fast simulated annealing) algorithm. Comparison of inversion results with those from transmission loss matching (an accompanying paper in this issue of the Journal of the Acoustical Society of Korea) shows that Parameters are consistently inverted. Direct time series comparisons between the measured signals and the simulated signals are Presented based on inversion results.

• The Journal of Engineering Geology
• /
• v.19 no.4
• /
• pp.509-515
• /
• 2009

Three earthquakes with local magnitude ($M_L$) greater than 3.0 occurred on April 24, June 2 and September 12 in 1999 nearby the Gyeongju area. Redetermined epicenters were located within the radius of 1 km. We carried out waveform inversion analysis to estimate focal mechanism of June 2 event, and P and S wave polarity and their amplitude ratio analysis to estimate focal mechanisms of April 24 and September 12 events. June 2 and September 12 events had similar fault plane solutions each other. The fault plane solution of April 24 event included those of other 2 events, but its distribution range was relatively broad. Focal mechanisms of those events had a strike slip faulting with a small normal component. P-axes of those events were ENE-WSW which were similar to previous studies on the P-axis of the Korean Peninsula. Considering distances between epicenters, similarities of seismic waves and sameness of polarities of seismic data recorded at common seismic stations, these events might occurred at the same fault. The seismic moment of June 2 event was estimated to be $3.9\;{\times}\;10^{14}\;N{\cdot}m$ and this value corresponded to the moment magnitude ($M_W$) 3.7. The moment magnitude estimated by spectral analysis was 3.8, which was similar to that estimated by waveform inversion analysis. The average stress drop was estimated to be 7.5 MPa. Moment magnitudes of April 24 and September 12 events were estimated to be 3.2 and 3.4 by comparing the spectrum of those events recorded at common single seismic station.

• 한국지구물리탐사학회:학술대회논문집
• /
• 2010.10a
• /
• pp.87-91
• /
• 2010

• Proceedings of the Korean Geotechical Society Conference
• /
• 1996.10a
• /
• pp.203-208
• /
• 1996

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.36 no.3
• /
• pp.203-211
• /
• 2023

A Markov chain Monte Carlo (MCMC) simulation is proposed for probabilistic full waveform inversion (FWI) in a layered half-space. Dynamic responses on the half-space surface are estimated using the thin-layer method when a harmonic vertical force is applied. Subsequently, a posterior probability distribution function and the corresponding objective function are formulated to minimize the difference between estimations and observed data as well as that of model parameters from prior information. Based on the gradient of the objective function, a proposal distribution and an acceptance probability for MCMC samples are proposed. The proposed MCMC simulation is applied to several layered half-space examples. It is demonstrated that the proposed MCMC simulation for probabilistic FWI can estimate probabilistic material properties such as the shear-wave velocities of a layered half-space.

• The Journal of the Acoustical Society of Korea
• /
• v.18 no.8
• /
• pp.32-41
• /
• 1999

In this paper, an inversion method using chirp signals and two near field receivers is proposed. Inversion problems can be formulated into the probabilistic models composed of signals, a forward model and noise. Forward model to simulate chirp signals is chosen to be the source-wavelet-convolution planewave modeling method. The solution of the inversion problem is defined by a posteriori pdf. The wavelet matching technique, using weighted least-squares fitting, estimates the sediment sound-speed and thickness on which determination of the ranges for a priori uniform distribution is based. The genetic algorithm can be applied to a global optimization problem to find a maximum a posteriori solution for determined a priori search space. Here the object function is defined by an L₂norm of the difference between measured and modeled signals. The observed signals can be separated into a set of two signals reflected from the upper and lower boundaries of a sediment. The separation of signals and successive applications of the genetic algorithm optimization process reduce the search space, therefore improving the inversion results. Not only the marginal pdf but also the statistics are calculated by numerical evaluation of integrals using the samples selected during importance sampling process of the genetic algorithm. The examples applied here show that, for synthetic data with noise, it is possible to carry out an inversion for sedimentary layers using the proposed inversion method.

• Geophysics and Geophysical Exploration
• /
• v.13 no.2
• /
• pp.137-143
• /
• 2010

An earthquake with local magnitude $(M_L)$ 3.7 on December 9, 2000 occurred offshore Yeongdeok area, South Korea. In case of applying Chang and Baag (2006) crustal velocity model, the epicenter is $36.4462^{\circ}N\;and\;129.9789^{\circ}E$, which belongs to the inside of the Korean Peninsula Continental Shelf. Although we use the modified model reducing crustal thickness of Chang and Baag (2006) model by 5 km considering the transition from continental crust to oceanic crust in the East Sea, the epicenter was little changed. We carried out the waveform inversion analysis to estimate focal depth and focal mechanism of this event. The focal depth is estimated to be 11 ~ 12 km. The seismic moment is estimated to be $1.0{\times}10^{15}N{\cdot}m$, and this value corresponds to the moment magnitude $(M_W)$ 3.9. The offshore Yeongdeok event including May 29, 2004 offshore Uljin one show typical thrust faulting, and the direction of P-axis is ESE-WNW. The moment magnitude estimated by the spectral analysis is 4.0, which is similar to that by the waveform inversion analysis. Average stress drop is estimated to be 3.4 MPa.

• Journal of the Korean GEO-environmental Society
• /
• v.25 no.2
• /
• pp.5-14
• /
• 2024

Reinforced concrete bridge decks are the first to be damaged by vehicle loads and rain infiltration. Concrete deterioration primarily occurs owing to the corrosion of rebars and other metal components by chlorides used for snow and ice melting. The structural condition and concrete deterioration of the bridge decks within the pavement were evaluated using ground-penetrating radar (GPR) survey data. To evaluate concrete deterioration in bridges, it is necessary to develop GPR data analysis techniques to accurately identify deteriorated locations and rebar positions. GPR exploration involves the acquisition of reflection and diffraction wave signals due to differences in radar wave propagation velocity in geotechnical media. Therefore, a full-waveform inversion (FWI) method was developed to evaluate the deterioration of reinforced concrete bridge decks by estimating the radar wave propagation velocity in geotechnical media using GPR data. Numerical experiments using a GPR velocity model confirmed the deterioration phenomena of bridge decks, such as concrete delamination and rebar corrosion, verifying the applicability of the developed technology. Moreover, using the synthetic GPR data, FWI facilitates the determination of rebar positions and concrete deterioration locations using inverted velocity images.

• 한국지구물리탐사학회:학술대회논문집
• /
• 2011.10a
• /
• pp.37-40
• /
• 2011