• Geophysics and Geophysical Exploration
• /
• v.20 no.3
• /
• pp.163-175
• /
• 2017

In land seismic exploration, irregular surface topography and weathering layer in near surface distorts the reflected signals of data. Therefore, typical land seismic data should be compensated for this distortion by static correction. To perform the static correction, near-surface velocity is required, which can be obtained by seismic refraction survey. However, land seismic data is often acquired in a limited form of geometry depending on the equipment availability, accessibility condition, and permission for the survey site. In this situation, refraction analysis should be performed using reflection data because it is impossible to acquire refraction-oriented data due to limited source and receiver geometry. In this study, we aimed to analyze the reliability of the results obtained by refraction traveltime tomography when using reflection data with a limited number of sources and receivers from irregular surface topography. By comparing the inversion result from irregular topography with that from flat surface, we found that the surface topography affects the reliability of the inversion results to some degree. We also found that the number of sources has little effect on the inversion results unless the number of sources are very small. On the other hand, we observed that velocity distortion occurred in the overlapped part of receiver arrays when using a limited number of receivers, and therefore suggested the size of the least overlapping ratio to avoid the velocity distortion. Finally, we performed numerical tests for the model which simulates the surface topography and acquisition geometry of the survey region and verified the reliability analysis of inversion results. We identified reliable areas and suspicious area of the inverted velocity model by applying the analysis results to field data.

• 한국지구물리탐사학회:학술대회논문집
• /
• 2000.09a
• /
• pp.132-146
• /
• 2000

Seismic refraction survey was performed for 10 lines along NE-SW and NW-SE directions above Nampoong gallery at Makyo-ri, Dogye, Samcheok, Kangwon-do. 48 geophones were laid in line with the interval of 1m, and a 5Kg hammer was used as a source at 5 points for each line. Data processing was done using reciprocal time method, GRM, and traveltime tomography which utilizes wavefront expansion method for forward process and STRT for inversion. The result shows that the first layer has its lower boundary between 3.49m and 8.88m. The P-wave velocity of the first and the second layer were estimated as 270${\~}$360m/s and 1550${\~}$1940m/s respectively. When the boundary of the first and second layer is smooth enough and the velocity difference is large enough, GRM has little advantage over reciprocal time method. The result of reciprocal method and traveltime tomography shows consistency. The northeast part of the boundary has syncline structure, which is similar to the topography above. This implies that the collapse of the cavities of Nampoong gallery result in the subsidence of the ground surface. The subsidence is in progress across the Youngdong railroad, therefore a proper reinforcement work is required.

• Geophysics and Geophysical Exploration
• /
• v.12 no.1
• /
• pp.105-113
• /
• 2009

The interpretation of observed waveform characteristics identified in refraction and wide-angle reflection data increases confidence in the crustal structure model obtained. When calculating traveltimes and raypaths, wavefront methods on a regular grid based on graph theory are robust even with complicated structures, but basically compute only first arrivals. In this paper, we develop new algorithms to compute traveltimes and raypaths not only for first arrivals, but also for fast and later reflection arrivals, later refraction arrivals, and converted waves between P and S, using the modified wavefront method based on slowness network nodes mapped on a multi-layer model. Using the new algorithm, we can interpret reflected arrivals, Pg-later arrivals, strong arrivals appearing behind Pn, triplicated Moho reflected arrivals (PmP) to obtain the shape of the Moho, and phases involving conversion between P and S. Using two models of an ocean-continent transition zone and an oceanic ridge or seamount, we show the usefulness of this algorithm, which is confirmed by synthetic seismograms using the 2D Finite Difference Method (2D-FDM). Characteristics of arrivals and raypaths of the two models differ from each other in that using only first-arrival traveltime data for crustal structure analysis involves risk of erroneous interpretation in the ocean-continent transition zone, or the region around a ridge or seamount.

• 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.1 no.1
• /
• pp.49-56
• /
• 1998

Seismic tomography has been widely used as high resolution subsurface imaging techniques in engineering applications. Although most of the techniques have been using travel time inversion, waveform method is being driven forward owing to the progress of computational environments. Although full-waveform inversion method has been known as the best method in terms of model resolving power without high-frequency restriction and weak scattering approximation, it has practical disadvantage that it is apt to get stuck in local minimum if the initial guess is far from the actual model and it consumes so much time to calculate. In this study, 2-D full-waveform inversion algorithm in acoustic medium is developed, which uses result of traveltime tomography as initial model. From the application on synthetic data, it is proved that this approach can efficiently reduce the problem of conventional approaches: our algorithm shows much faster convergence rate and improvement of model resolution. Result of application on physical modeling data also shows much improvement. It is expected that this algorithm can be applicable to real data.

• Journal of Industrial Technology
• /
• v.27 no.A
• /
• pp.3-7
• /
• 2007

There are many techniques to calculate the exact position of deep seated tunnel. Especially, tomography method has been used generally in present days. This method has been performed mainly by wave traveltime. Because of short interval of two measuring boreholes, it was very hard to interpret the exact tunnel position. To solve this problem, seismic amplitude method was tried to detect exact pososition of tunnel in this study.

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

• Journal of the Korean Geophysical Society
• /
• v.8 no.2
• /
• pp.81-87
• /
• 2005

Acquisition and processing of vector seismic waves were conducted through simultaneous generation of P, SH, and SV waves and receiving those waves using three-component geophones. Test data were received by 24 8-Hz geophones at an interval of 2 m along a 94-m profile. The data were recorded for 512 ms with sampling intervals of 0.2 ms. Raw data indicate that both reflected and refracted P waves are strongly recorded on the vertical component while SH waves are significant on the transverse horizontal component. On the inline horizontal component, both direct P and converted PS waves are recorded. First arrivals of P and SH waves were detected simultaneously on the vertical and transverse horizontal axes, respectively. The recorded vector data were separately inverted using traveltime tomography to yield P- and SH-wave sections. Using those two velocity sections, Poisson's ratios were able to be obtained effectively.

• Geophysics and Geophysical Exploration
• /
• v.4 no.3
• /
• pp.84-88
• /
• 2001

For the elastic migation, the velocity errors between the initial velocity model and true velocity model seriously affect the migrated images. The assumption of an initial velocity model, thus, is one of the critical factor for the successful migration. In case of applying the layered earth model as an initial velocity model, the layer boundary having large velocity contrast can not be defined well with conventional traveltime calculation algolithms and we have the difficulties for expressing the characteristics of the real subsurface. Smooth Background Model (SBM) we have applied as an initial velocity model in our study is characterized to be linearly varying the velocity with the depth, which can express the velocity variation in the subsurface properly. Thus it can properly be applied to traveltime calculation algolithms such as Vidale's method. In this study, Kirchhoff operator for prestack migration was used and the absolute amplitude obtained by modeling was applied as a weighted value to consider the true amplitude for initial model. Initial velocity model for migration was determined by using stacking velocity and we applied this model to real data.

• 한국지구물리탐사학회:학술대회논문집
• /
• 1999.08a
• /
• pp.137-155
• /
• 1999

Recently, Televiewer(Borehole Acoustic Scanner(Televiewer)) has come to be widely used specially for the general engineering construction design. The Televiewer tool using a focussed acoustic beam is to detect the amplitude and traveltime of each reflected acoustic signal at the wall, resulting in the amplitude- and traveltime image respectively. Fractures can be well detected, because they easily scatter the acoustic energy due to the highly narrow beam. In addition, the drilling work will rough the borehole wall so that the acoustic energy can be scattered simply due to the roughness of the wall. Thus, the amplitude level can be directed associated with the elastic properties(impedance) and the hardness of the rock as well. Meanwhile, the traveltime image provides an information about the borehole shape and can be converted to a high precision 3D caliper log(max. 288 arms). In this paper, based on the high resolution of Televiewer images, general evaluation methods are illustrated to derive very reliable rock parameters.