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

• Journal of the korean society of oceanography
• /
• v.32 no.4
• /
• pp.163-170
• /
• 1997

Some CMP gathers acquired from shallow marine seismic reflection survey in offshore Korea do not show the hyperbolic trend of moveout. It originated from so-called swell effect of source and streamer, which are towed under rough sea surface during the data acquisition. The observed time deviations of NMO-corrected traces can be entirely ascribed to the swell effect. To correct these time deviations, a residual statics is introduced using Genetic Algorithms (GA) into the swell correction. A new class of global optimization methods known as GA has recently been developed in the field of Artificial Intelligence and has a resemblance with the genetic evolution of biological systems. The basic idea in using GA as an optimization method is to represent a population of possible solutions or models in a chromosome-type encoding and manipulate these encoded models through simulated reproduction, crossover and mutation. GA parameters used in this paper are as follows: population size Q=40, probability of multiple-point crossover P$_c$=0.6, linear relationship of mutation probability P$_m$ from 0.002 to 0.004, and gray code representation are adopted. The number of the model participating in tournament selection (nt) is 3, and the number of expected copies desired for the best population member in the scaling of fitness is 1.5. With above parameters, an optimization run was iterated for 101 generations. The combination of above parameters are found to be optimal for the convergence of the algorithm. The resulting reflection events in every NMO-corrected CMP gather show good alignment and enhanced quality stack section.

• 한국신재생에너지학회:학술대회논문집
• /
• 2006.06a
• /
• pp.377-381
• /
• 2006

A seismic survey for gas hydrate have performed over the East sea by the KIGAM since 1997. General indicator of gas hydrate in seismic data is commonly inferred from the BSR(Bottom Simulating Reflector) that occurred parallel to the sea floor, amplitude decrease at the top of the BSR, amplitude blanking at the bottom of the BSR, decrease of the interval velocity and the reflection phase reversal at the BSR. In this paper we had analyzed optimum parameters of the field data to detect the 9as hydrate. Shot delay correction is applied 95ms, spherical divergence correction is applied velocity library 3, bandpass filter is applied 25-30-115-120Hz deconvolution operator length is applied 60ms, lag is 6ms and accurate velocity analysis NMO correction, stack is performed. Geobit 2.11.0 developed by the KIGAM was used for all data processing. Processing results say that the BSR occurred parallel to the sea floor were shown at 3,150m/s of two way travel time from the sea floor through shot point 5,000-5,610, and identified the interval velocity decrease around BSR and the reflection phase reversal corresponding to the reflection at the sea floor.

• The Journal of Engineering Geology
• /
• v.16 no.2 s.48
• /
• pp.161-169
• /
• 2006

A short seismic reflection survey was performed to map the basement and the upper units in the Gomso Bay. This research was mainly aimed at clarifying the basement by improving the signal-to-noise ratio in data processing steps. The strategies employed in this research included enhancement of the signal interfered with large-amplitude noise, through pre- and post-stack processing such as time-variant filtering, bad trace edit, careful muting after f-k filter and NMO correction. The subsurface structure mapped from this survey mainly consists of the top of basement and the upper three units, which were well correlated to the result from the previously conducted MT survey. Furthermore seismic section clarifies approximately 30m deep subhorizontal event of the top of the basement, which was not shown in the central portion of the MT section due to data qualify.

• The Journal of Engineering Geology
• /
• v.28 no.3
• /
• pp.465-473
• /
• 2018

The application of singular value decomposition (SVD) filtering is examined for attenuation of the ground-roll in land seismic data. Prior to the SVD computation to seek singular values containing the highly correlatable reflection energy, processing steps such as automatic gain control, elevation and refraction statics, NMO correction, and residual statics are performed to enhance the horizontal correlationships and continuities of reflections. Optimal parameters of SVD filtering are effectively chosen with diagnostic display of inverse NMO (INMO) corrected CSP (common shot point) gather. On the field data with dispersion of ground-roll overwhelmed, continuities of reflection events are much improved by SVD filtering than f-k filtering by eliminating the ground-roll with preserving the low-frequency reflections. This is well explained in the average amplitude spectra of the f-k and SVD filtered data. The reflectors including horizontal layer of the reservoir are much clearer on the stack section, with laminated events by SVD filtering and subsequent processing steps of spiking deconvolution and time-variant spectral whitening.

• Geophysics and Geophysical Exploration
• /
• v.4 no.4
• /
• pp.133-144
• /
• 2001

An effective seismic reflection technique for mapping the cavities and bedrock surface in carbonate rocks is described. The high resolution seismic reflection images were successfully registered by using the hydrophones employed in the stream-water driven trench, and were effectively focused by applying optimal data processing sequences. The strategy included enhancement of the signal interfered with the large-amplitude scattering noise, through pre- and post stack processing such as time-variant filtering, bad-trace editing, residual statics, velocity analysis, and careful muting after NMO (normal moveout) correction. The major reflections including the bedrock surface were mapped with the desired resolution and were correlated to the seismic crosshole tomographic data. Shallow major reflectors could be identified and analyzed on the AGC (auto gain control)-applied field records. Three subhorizontal layers were identified with their distinct velocities; overburden (<3000 m/s), sediments (3000-4000 m/s), limestone bedrock (>4000 m/s). Taking into account of no diffraction effects in the field records, gravel-rich overburdens and sediments are considered to be well sorted. Based on the images mapped consistently on the whole survey line and seismic velocity increasing with depth, this area probably lacks in sizable cavities (if any, no air-filled cavities).

• Geophysics and Geophysical Exploration
• /
• v.13 no.1
• /
• pp.31-39
• /
• 2010

The goal of automatic velocity analysis is to extract accurate velocity from voluminous seismic data with efficiency. In this study, we developed an efficient automatic velocity analysis algorithm by using bootstrapped differential semblance (BDS) and Monte Carlo inversion. To estimate more accurate results from automatic velocity analysis, the algorithm we have developed uses BDS, which provides a higher velocity resolution than conventional semblance, as a coherency estimator. In addition, our proposed automatic velocity analysis module is performed with a conditional initial velocity determination step that leads to enhanced efficiency in running time of the module. A new optional root mean square (RMS) velocity constraint, which prevents picking false peaks, is used. The developed automatic velocity analysis module was tested on a synthetic dataset and a marine field dataset from the East Sea, Korea. The stacked sections made using velocity results from our algorithm showed coherent events and improved the quality of the normal moveout-correction result. Moreover, since our algorithm finds interval velocity ($\nu_{int}$) first with interval velocity constraints and then calculates a RMS velocity function from the interval velocity, we can estimate geologically reasonable interval velocities. Boundaries of interval velocities also match well with reflection events in the common midpoint stacked sections.

• Geophysics and Geophysical Exploration
• /
• v.19 no.1
• /
• pp.11-19
• /
• 2016

Two representative residual static methods of traveltime decomposition and stack-power maximization are discussed in terms of application to land seismic data. For the model data with synthetic shot/receiver statics (time shift) applied and random noises added, continuities of reflection event are much improved by stack-power maximization method, resulting the derived time-shifts approximately equal to the synthetic statics. Optimal parameters (maximum allowable shift, correlation window, iteration number) for residual statics are effectively chosen with diagnostic displays of CSP (common shot point) stack and CRP (common receiver point) stack as well as CMP gather. In addition to removal of long-wavelength time shift by refraction statics, prior to residual statics, processing steps of f-k filter, predictive deconvolution and time variant spectral whitening are employed to attenuate noises and thereby to minimize the error during the correlation process. The reflectors including horizontal layer of reservoir are more clearly shown in the variable-density section through repicking the velocities after residual statics and inverse NMO correction.

• Journal of the Korean Geophysical Society
• /
• v.2 no.1
• /
• pp.45-56
• /
• 1999

Seismic refracrion and reflection surveys were conducted along an E-W trending track of 482 m long in Ilwall-dong, Pohang. End-on spread was employed as source-receiver configuration with 2 m for both geophone interval and offset. Seismic data were acquired using 24 channels at every shot fired every 2 m along the track. Refraction data were interpreted using equations for multi-horizontal layers. Reflection data were processed in the sequence of trace edit, gain control, CMP sorting, NMO correction, mute, common offset gathering, and filtering to produce a single fold seismic section. There are two layers in shallow subsurface of the study area. Upper layer has the P-wave velocities ranging from 267 to 566 m/s and is interpreted as a layer of unconsolidated sediments. Lower layer has P-wave velocities of 1096-3108 m/s and is interpreted as weathered rock to hard rock. Most of the lower layer classified as soft rock. Upper layer has lateral variations in both P-wave velocity and thickness. The upper layer in the eastern part of the seismic line is 3-5 m thick and has P-wave velocity of 400 m/s in average. The upper layer in the western part is 8-10 m thick and has P-wave velocity of 340 m/s in average. The eastern part is interpreted as unconsolidated beach sand, while the western part is interpreted as infilled soil to develop a construction site. Three fault systems of high angle are imaged in seismic reflection section. It is interpreted that the area between these fault systems are relatively safe. Large buildings should be located in the safe ground condition of no fault and footings should be designed to be in the basement rock of 3-10 m deep below the surface.

• Geophysics and Geophysical Exploration
• /
• v.2 no.1
• /
• pp.26-32
• /
• 1999

Among the various seismic data processing sequences, the velocity analysis is the most time consuming and man-hour intensive processing steps. For the production seismic data processing, a good velocity analysis tool as well as the high performance computer is required. The tool must give fast and accurate velocity analysis. There are two different approches in the velocity analysis, batch and interactive. In the batch processing, a velocity plot is made at every analysis point. Generally, the plot consisted of a semblance contour, super gather, and a stack pannel. The interpreter chooses the velocity function by analyzing the velocity plot. The technique is highly dependent on the interpreters skill and requires human efforts. As the high speed graphic workstations are becoming more popular, various interactive velocity analysis programs are developed. Although, the programs enabled faster picking of the velocity nodes using mouse, the main improvement of these programs is simply the replacement of the paper plot by the graphic screen. The velocity spectrum is highly sensitive to the presence of the noise, especially the coherent noise often found in the shallow region of the marine seismic data. For the accurate velocity analysis, these noise must be removed before the spectrum is computed. Also, the velocity analysis must be carried out by carefully choosing the location of the analysis point and accuarate computation of the spectrum. The analyzed velocity function must be verified by the mute and stack, and the sequence must be repeated most time. Therefore an iterative, interactive, and unified velocity analysis tool is highly required. An interactive velocity analysis program, xva(X-Window based Velocity Analysis) was invented. The program handles all processes required in the velocity analysis such as composing the super gather, computing the velocity spectrum, NMO correction, mute, and stack. Most of the parameter changes give the final stack via a few mouse clicks thereby enabling the iterative and interactive processing. A simple trace indexing scheme is introduced and a program to nike the index of the Geobit seismic disk file was invented. The index is used to reference the original input, i.e., CDP sort, directly A transformation techinique of the mute function between the T-X domain and NMOC domain is introduced and adopted to the program. The result of the transform is simliar to the remove-NMO technique in suppressing the shallow noise such as direct wave and refracted wave. However, it has two improvements, i.e., no interpolation error and very high speed computing time. By the introduction of the technique, the mute times can be easily designed from the NMOC domain and applied to the super gather in the T-X domain, thereby producing more accurate velocity spectrum interactively. The xva program consists of 28 files, 12,029 lines, 34,990 words and 304,073 characters. The program references Geobit utility libraries and can be installed under Geobit preinstalled environment. The program runs on X-Window/Motif environment. The program menu is designed according to the Motif style guide. A brief usage of the program has been discussed. The program allows fast and accurate seismic velocity analysis, which is necessary computing the AVO (Amplitude Versus Offset) based DHI (Direct Hydrocarn Indicator), and making the high quality seismic sections.