• Geophysics and Geophysical Exploration
• /
• v.5 no.1
• /
• pp.18-22
• /
• 2002

As a single arrival traveltime, maximum energy arrival traveltime has been known as the most proper operator for Kirchhoff migration. In case of the model having the simple structure, both the first arrival traveltime and the maximum energy arrival traveltime can be used as the correct operators for Kirchhoff migration. However for some model having the complex and high velocity contrast structure, the migration using the first arrival traveltime can't give the correct depth section. That is, traveltime to be required in Kirchhoff migration is the maximum energy traveltime, but, needs considerably more calculation time than that of first arrival. In this paper, we propose the method for calculating the traveltime approximated to the maximum energy arrival using one-way wave equation. After defining the WAS(Wrap Around Suppression) factor to be used for calculating the first arrival traveltime using one-way wave equation as the function of lateral grid interval and depth and considering the delay time of source wavelet. we calculate the traveltime approximated to the maximum energy arrival. to verify the validity of this traveltime, we applied this to the migraion for simple structure and complex structure and compared the depth section with that obtained by using the first arrival traveltime.

• Tunnel and Underground Space
• /
• v.11 no.2
• /
• pp.120-121
• /
• 2001

The purpose of this study is to investigate concrete flaw using seismic first arrival and various inversion method. Seismic wave propagation was calculated using finite element method in theoretical modelling and tomogram was made using various inversion methods in theoretical and experimental modelling. Five steps of seismic first arrival were selected from FEM results and these data were used to calculate seismic velocity section. According to the results, exact seismic first arrival picking method was proposed and experimental modelling was conducted.

• Geophysics and Geophysical Exploration
• /
• v.22 no.1
• /
• pp.12-20
• /
• 2019

The first-arrival traveltime calculation method based on the damped wave equation overcomes the shortcomings of ray-tracing methods. Since this algorithm needs to solve the damped wave equation, numerical modeling is essential. However, it is not desirable to use the finite-difference method (FDM), which has good computational efficiency, for simulating the land seismic data because of irregular topography. Thus, the finite-element method (FEM) which requires higher computational cost than FDM has been used to correctly describe the irregular topography. In this study, we computed first-arrival traveltimes in an irregular topographic model using FDM incorporating embedded boundary method (EBM) to overcome this problem. To verify the accuracy and efficiency of the proposed algorithm, we compared our results with those of FEM. As a result, the proposed method using EBM not only provided the same accuracy as the FEM but also showed the improved computational efficiency.

• Geophysics and Geophysical Exploration
• /
• v.14 no.2
• /
• pp.146-151
• /
• 2011

A method to evaluate and to reduce the source- and receiver- consistent noise in a cross-well travel time data was proposed. These systematic noises, which can cause some serious effects on the result of a travel time tomography, can be considered as the source and receiver statics. The method evaluates the statics through a curve-fitting of the first arrival travel times in the common source and common receiver gathers. Feasibility study was conducted on a synthetic data which simulates the cross-well travel time tomography to detect a small scale tunnel in a uniform background medium. First arrival travel times at a given source and receiver points are computed by a raytracing method, and the source consistent- and receiver consistent noises are added to the record. In case of the added noise with rms amounting to 25% of the maximum expected anomalous travel time delays, it is confirmed that the method successfully extracted the noise at the 7th step of iteration.

• Geophysics and Geophysical Exploration
• /
• v.12 no.3
• /
• pp.233-238
• /
• 2009

An inversion method regulated by the error in the measurement of the first arrival time was developed, and we conducted a feasibility study by applying the method to a real cross-well seismic data. The inversion is a two-step regulation process; 1) derive the measurement error bound based on the resolution of the velocity image want to derive, and exclude the records whose picking error is larger than the error bound, 2) set the travel time residual to zero if the residual is less than the measurement error. This process prevents the trivial residuals are accumulated and contribute to the velocity-model update. Comparison of two velocity images, one by using all records and another by using the regulate inversion method, shows that the later velocity image exhibits less numerical artefacts, and it also indicates that, according to the Fermat's principle, the latter image is a more feasible velocity model.

• Geophysics and Geophysical Exploration
• /
• v.20 no.1
• /
• pp.1-11
• /
• 2017

In order to improve the accuracy of microseismic epicenter location through the inversion techniques using P and S wave first arrivals, field experiments of microseismic monitoring were performed using borehole 3-component geophones. The direction of epicenter was estimated from the hodograms of P-wave first arrivals through the weight drop experiments in which the $\times$ component of 3-component geophone was aligned to the magnetic north. The picking of S wave first arrival was possible in the polarization filtered data even if S waves are difficult to be identified in raw data. The inversion technique using only P wave first arrival times can often converge to the local minimum when the initial values for epicenter are largely apart from the true epicenter, so that the correct solution can not be found. To solve this problem, the epicenter determination method using differences between P and S wave arrival times was used to estimate proper initial values of epicenter. The inversion result using only P-wave first arrival times which started from the estimated initial values showed the improved accuracy of the epicenter location.

• Geophysics and Geophysical Exploration
• /
• v.13 no.4
• /
• pp.330-335
• /
• 2010

In order to estimate the confidence level of the velocity distribution shown in a velocity image reconstructed from a travel-time tomography, the ray coverage and the inversion characteristics of the system matrix were investigated. The targets of the analysis is the first arrival travel-time, the raypath information, and the resulting velocity model. The ray coverage, degree of ray and model coupling, was estimated by the number of rays and total ray length in a velocity grid, and information regarding the resolution and uncertainties involved in the reconstructed velocity model was derived from the results of the SVD analysis of the system matrix that relates the data space (first arrival travel times) to the model space (velocity distribution in tomogram).

• Geophysics and Geophysical Exploration
• /
• v.14 no.2
• /
• pp.140-145
• /
• 2011

In order to analyze the influence of the noise on a cross-well traveltime tomography to detect a small scale low velocity body in a homogeneous medium, the first arrival travel times were computed one a tunnel model by a finite-difference ray tracing scheme. Three different types and four different intensity levels of white noises were added to the computed first arrival travel times, and velocity tomograms were constructed using an iterative inversion method (SIRT). Tomograms with the noise intensity up to 10% of the maximum traveltime delay in the tunnel model, showed the exact location of the tunnel. However, the velocity shown at the tunnel location was not close to air velocity but only slightly less than the velocity of the background medium. The additive random noise showed significantly less degree of influence on the resulting tomogram than the source- and receiver consistent noise.

• Geophysics and Geophysical Exploration
• /
• v.5 no.3
• /
• pp.150-156
• /
• 2002

In this study, a software for crosswell seismic tomography is developed. The software consists of first arrival picking and adjusting module, crosswell traveltime tomography module, and imaging module. This software allows saying the picked first arrival times into the header of seismic data, and using this data directly to the input of crosswell seismic tomography. With an imaging module, velocity structures and ray path can be imaged directly from the output of the tomography module. Because it is developed on the basis of the SU under the Linux and the GUI environment for user, this software can be carried out directly the first arrival picking, inversion and tomogram for crosswell tomography data in the field. Therefore, this software can be improved the applicability of site investigation by tomography method.

• Journal of the Korean Geophysical Society
• /
• v.9 no.1
• /
• pp.37-45
• /
• 2006

Velocity structures were defined in the vicinty of the 140-m deep test borehole in the pungam basin through simultaneous inversion of surface seismic refraction and far-ofset VSP traveltime data. Seismicenergy generated at the surface by a seisgun was recorded both at 42 surface locations at 3-m intervalsalong the profiles in the N20E and its orthogonal directions and at 71 m depth in the borehole. Forthe ofset VSP study, seismic energy was generated by a 5 kg sledgehamer at the surface in the horizontal ofset range of -19.5∼+19.5 m from the borehole. The seismic signals were detected at 9∼99 m depths with 1∼2 m intervals and recorded for 204 ms per shot. After shot static corrections,first-arrival times picked from both the surface refraction and borehole records were simultaneouslyinverted to yield velocity tomograms. The tomograms indicate that a 1.5 m thick soil layer with velocities les than 500 m/s overlies basements having a velocity range of 3,067 ∼5,717 m/s. Within the basements,∼4 m and deeper than 71 m. The high-velocit yzones may be due to conglomerates intercalated with sandstones and siltstones. No evidence for large-scale fracture zones or faults is detected near the borehole