• Geophysics and Geophysical Exploration
• /
• v.10 no.2
• /
• pp.131-137
• /
• 2007

With consideration of the first arrival picking methodology and inherent errors in picking process, I propose, from the computerization point of view, a practical algorithm for picking and error computation. The proposed picking procedure consists of 2-step; 1) picking the first coherent peak or trough events, 2) derive a line which approximates the record in the interval prior to the pick, and set the intercept time of the line as the first break. The length of fitting interval used in experiment, is few samples less than 1/4 width of the arriving wavelet. A quantitative measure of the error involved in first arrival picking is defined as the time length that needed to determine if an event is the first arrival or not. The time length is expressed as a function of frequency bandwidth of the signal and the S/N ratio. For 3 sets of cross-well seismic data, first breaks are picked twice, by manually, and by the proposed method. And at the same time, the error bound for each trace is computed. Experiment results show that good performance of the proposed picking method, and the usefulness of the quantitative error measure in pick-quality evaluation.

• 한국지구물리탐사학회:학술대회논문집
• /
• 2008.10a
• /
• pp.69-72
• /
• 2008

EZTOMO is a crosswell seismic tomography software system. The system has capability of event picking, raytracing, inversion, error analysis, and visualization of the processing results. Waveform of the first arrival signal has been utilized to select the event of the first motion, and uncertainty measured in estimation of the first breaks has been utilized to improve the inversion process.

• Tunnel and Underground Space
• /
• v.29 no.4
• /
• pp.243-261
• /
• 2019

Microseismic monitoring technologies have been recognized for its superiority over traditional methods and are used in domestic and overseas underground mines. However, the complex gangway layout of underground mines in Korea and the mixed structure of excavated space and rock masses make it difficult to estimate the microseismic propagation and to determine the arrival time of microseismic wave. In this paper, experimental studies were carried out to determine the source location according to various arrival picking methods and dimensional conditions. The arrival picking methods used were FTC (First Threshold Cross), Picking window, AIC (Akaike Information Criterion), and 2-D and 3-D source generation experiments were performed, respectively, under the 2-D sensor array. In each experiment, source location algorithm used iterative method and genetic algorithm. The iterative method was effective when the sensor array and source generation were the same dimension, but it was not suitable to apply when the source generation was higher dimension. On the other hand, in case of source location using RCGA, the higher dimensional source location could be determined, but it took longer time to calculate. The accuracy of the arrival picking methods differed according to the source location algorithms, but picking window method showed high accuracy in overall.

• Geophysics and Geophysical Exploration
• /
• v.2 no.1
• /
• pp.17-25
• /
• 1999

A traveltime tomography has been carried out by transforming electromagnetic data in frequency domain to wave-like domain. The transform uniquely relates a field satisfying a diffusion equation to an integral of the corresponding wavefield. But direct transform of frequency domain magnetic fields to wave-field domain is ill-posed problem because the kernel of the integral transform is highly damped. In this study, instead of solving such an unstable problem, it is assumed that wave-fields in transformed domain can be approximated by sum of ray series. And for further simplicity, reflection and refraction energy compared to that of direct wave is weak enough to be neglected. Then first arrival can be approximated by calculating the traveltime of direct wave only. But these assumptions are valid when the conductivity contrast between background medium and the target anomalous body is low enough. So this approach can only be applied to the models with low conductivity contrast. To verify the algorithm, traveltime calculated by this approach was compared to that of direct transform method and exact traveltime, calculated analytically, for homogeneous whole space. The error in first arrival picked by this study was less than that of direct transformation method, especially when the number of frequency samples is less than 10, or when the data are noisy. Layered earth model with varying conductivity contrasts and inclined dyke model have been successfully imaged by applying nonlinear traveltime tomography in 30 iterations within three CPU minutes on a IBM Pentium Pro 200 MHz.

• Geophysics and Geophysical Exploration
• /
• v.6 no.2
• /
• pp.71-76
• /
• 2003

Accurate estimation of the first arrival travel time is an essential task to obtain a high resolution velocity tomogram. Accuracy of the travel time estimation may be influenced by two factors; geological and mechanical. A serious mechanical factor is the source firing control problems. We found the control problems in the records generated by tome impulsive borehole sources. The problems are; irregular firing control and uncertainty in estimation of the absolute firing-times shown in records. Definitely, the time difference will introduce an error to the first arrival times, and accordingly; it will cause some distortion in the resulting velocity tomogram. A method to determine the firing time is suggested here. The method determines the optimum onset time by comparing the horizontal and the NMO velocity with various amount of delay time adjustment.

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

$CO_2$-EOR (Carbon Dioxide-Enhanced Oil Recovery), one of the enhanced oil recovery methods, helps to not only enhance the production of oil, but also store carbon dioxide in underground. However, if micro fractures occur when during the injection of $CO_2$, it is difficult to make permanent storage of $CO_2$ in reservoir and can cause contamination of groundwater and soil. Therefore, in this study, we performed microseismic monitoring to investigate the occurrence of fractures during the $CO_2$ injection at the Meruap oil reservoir, Indonesia. To pick the first arrivals of microseismic events, Improved MER (Modified Energy Ratio) method was used. After picking the first arrivals, hodogram analysis was carried out by using the data recorded at three component geophones to calculate the back azimuth of events. Finally, locations of microseismic events were decided by using the results of first arrival picking and hodogram analysis. Estimated locations showed that all microseismic events were occurred at surface and any fracture did not occur around the reservoir. Moreover, by analyzing noise characteristic, we confirmed that almost of picked first arrivals were due to the repetitive mechanical noise.

• Geophysics and Geophysical Exploration
• /
• v.5 no.3
• /
• pp.193-198
• /
• 2002

In crosswell ray tomography, the resultant velocity structure could be affected by source static, first-arrival-time picking errors, convergence to a local minimum due to an inappropriate initial velocity model and etc. In the paper, I propose an algorithm that automatically correct the souce static among these error-prone factors. The algorithm automatically corrects source static using the picking times' differences along the source direction. The application of the algorithm to real data produces a quite satisfactory result. Tile algorithm seems to be helpful for users to apply the souce static correction consistently and to acquire accurate velocity structure.

• The Journal of Engineering Geology
• /
• v.29 no.2
• /
• pp.153-162
• /
• 2019

Rapid variations in the geometry (i.e., thickness) of the refractor and low velocities affect greatly the imaging of the reflectors of land seismic data. Conventional solutions to obtain the weathering models utilizes first break picking method, which requires time consuming steps and causes the human error in picking the first arrivals. A new interferometric approach (interferometric refraction statics, IRS) which utilizes the first arrival signal (S/N enhanced by refraction convolution stack) instead of first break picking, is tested in this study to the synthetic data from the velocity structure provided by surface geophysics (refraction, MASW) and borehole geophysics (tomography, SPS logging) for the Cheongju granitic bodies. The results of IRS approach are found to be better than the ones from conventional first break picking in terms of continuities and horizontal resolution of the reflectors. The unresolved long-wavelength statics in brute stack are much removed by IRS weathering correction and the overlying refractors (the base of shallow weathering zone) are incidentally delineated in the refraction convolution stack.

• Economic and Environmental Geology
• /
• v.28 no.5
• /
• pp.481-485
• /
• 1995

This paper presents a new approaching method to determine the velocity and geometry of shallow subsurface from seismic refraction events. After picking the first breaks from seismic refraction data, we assume that field refraction seismogram can be replaced by the unit delta function having time shift of first break. Time curves are generated by shooting ray tracing. The partial derivatives seismogram for a damped least squares method is computed analytically at each step of the forward ray tracing. The technique is successfully tested on synthetic and real data. It has the advantage of real full waveform inversion, which is robust at low frequency band even if the initial guess is far from the true model.

• Geophysics and Geophysical Exploration
• /
• v.25 no.4
• /
• pp.177-188
• /
• 2022

Recently, the application of distributed acoustic sensors (DAS), which can replace geophones and seismometers, has significantly increased along with interest in micro-seismic monitoring technique, which is one of the CO₂ storage monitoring techniques. A significant amount of temporally and spatially continuous data is recorded in a DAS monitoring system, thereby necessitating fast and accurate data processing techniques. Because event detection and seismic phase picking are the most basic data processing techniques, they should be performed on all data. In this study, a machine learning-based P, S wave phase picking algorithm was developed to compensate for the limitations of conventional phase picking algorithms, and it was modified using a transfer learning technique for the application of DAS data consisting of a single component with a low signal-to-noise ratio. Our model was constructed by modifying the convolution-based EQTransformer, which performs well in phase picking, to the ResUNet structure. Not only the global earthquake dataset, STEAD but also the augmented dataset was used as training datasets to enhance the prediction performance on the unseen characteristics of the target dataset. The performance of the developed algorithm was verified using K-net and KiK-net data with characteristics different from the training data. Additionally, after modifying the trained model to suit DAS data using the transfer learning technique, the performance was verified by applying it to the DAS field data measured in the Pohang Janggi basin.