• Geophysics and Geophysical Exploration
• /
• v.19 no.3
• /
• pp.111-120
• /
• 2016

In order to enhance the connectivity of fracture network as fluid path in enhanced/engineered geothermal system (EGS), the exact locating of hydraulic fractured zone is very important. Hydraulic fractures can be tracked by locating of microseismic events which are occurred during hydraulic fracture stimulation at each stage. However, since the subsurface velocity is changed due to hydraulic fracturing at each stage, in order to find out the exact location of microseismic events, we have to consider the velocity change due to hydraulic fracturing at previous stage when we perform the mapping of microseimic events at the next stage. In this study, we have modified 3D locating algorithm of microseismic data which was developed by Kim et al. (2015) and have developed 3D velocity update algorithm using occurred microseismic data. Eikonal equation which can efficiently calculate traveltime for complex velocity model at anywhere without shadow zone is used as forward engine in our inversion. Computational cost is dramatically reduced by using Fresnel volume approach to construct Jacobian matrix in velocity inversion. Through the numerical test which simulates the geothermal survey geometry, we demonstrated that the initial velocity model was updated by using microseismic data. In addition, we confirmed that relocation results of microseismic events by using updated velocity model became closer to true locations.

• 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.18 no.2
• /
• pp.64-73
• /
• 2015

Recently, hydraulic fracturing is used in various fields and microseismic monitoring is one of the best methods for judging where hydraulic fractures exist and how they are developing. When locating microseismic events using single vertical well data, distances from the vertical array and depths from the surface are generally decided using time differences between compressional (P) wave and shear (S) wave arrivals and azimuths are calculated using P wave hodogram analysis. However, in field data, it is sometimes hard to acquire P wave data which has smaller amplitude than S wave because microseismic data often have very low signal to noise (S/N) ratio. To overcome this problem, in this study, we developed a grid search algorithm which can find event location using all combinations of arrival times recorded at receivers. In addition, we introduced and analyzed the method which calculates azimuths using S wave. The tests of synthetic data show the inversion method using all combinations of arrival times and receivers can locate events without considering the origin time even using only single phase. In addition, the method can locate events with higher accuracy and has lower sensitivity on first arrival picking errors than conventional method. The method which calculates azimuths using S wave can provide reliable results when the dip between event and receiver is relatively small. However, this method shows the limitation when dip is greater than about $20^{\circ}$ in our model test.

• The Journal of Engineering Geology
• /
• v.19 no.2
• /
• pp.139-144
• /
• 2009

The microseismic monitoring system with wide range of frequency has been operating in real time and it is remotely monitored at indoor and on-site for one year. This system was constructed and established in order to secure the safe and effective operation of the KAERI Underground Research Tunnel(KURT). For one year monitoring work, total 14 events were recorded in the vicinity of the KURT, and the majority of events are regarded as ultramicroseismic earthquake and artificial impacts around the tunnel. The major event is the magnitude 3.4 earthquake which was centered around Gongju city, Chungnam Province. It means that there is no significant evidence of high frequency microseismic event, which is associated with fracture initiation and/or propagation in the rock mass and shotcrete. Three components sensor was applied in order to analyze and define the direction of vibration as well as an epicenter of microseismic origin, and also properly designed and installed in a small borehole. This monitoring system is able to predict the location and timing of fracturing of rock mass and rock fall around an undreground openings as well as analysis on safety of various kinds of engineering structures such as nuclear facilities and other structures.

• Tunnel and Underground Space
• /
• v.18 no.1
• /
• pp.1-9
• /
• 2008

Acoustic emission (AE) and Microseimsic (MS) activities are law-energy seismic events associated with a sudden inelastic deformation such as the sudden movement of existing fractures, the generation of new fractures or the propagation of fractures. These events rapidly increase before major failure and happen within a given rock volume and radiate detectable seismic waves. The main difference between AE and MS signals is that the seismic motion frequencies of AE signals are higher than those of MS signals. As the failure of geotechnical structures usually happens as a high velocity and small displacement, it is nat easy ta determine the precursor and initiation stress level of failure in displacement detection method. To overcame this problem, AE/MS techniques far detection of structure failure and damage have recently adapt in civil engineering. This study deal with the basic theory of AE/MS and state of arts in monitoring technique using AE/MS.

• Geophysics and Geophysical Exploration
• /
• v.13 no.1
• /
• pp.109-117
• /
• 2010

To enable high-performance computing (HPC) for applications with large datasets using a $Sony^{(R)}$ $PLAYSTATION^{(R)}$ 3 ($PS3^{TM}$) video game console, we configured a hybrid system consisting of a $Windows^{(R)}$ PC and a $PS3^{TM}$. To validate this system, we implemented the real-time multiplet identifier (RTMI) application, which identifies multiplets of microearthquakes in terms of the similarity of their waveforms. The cross-correlation computation, which is a core algorithm of the RTMI application, was optimised for the $PS3^{TM}$ platform, while the rest of the computation, including data input and output remained on the PC. With this configuration, the core part of the algorithm ran 69 times faster than the original program, accelerating total computation speed more than five times. As a result, the system processed up to 2100 total microseismic events, whereas the original implementation had a limit of 400 events. These results indicate that this system enables high-performance computing for large datasets using the $PS3^{TM}$, as long as data transfer time is negligible compared with computation time.

• Proceedings of the Korean Geotechical Society Conference
• /
• 2008.03a
• /
• pp.199-210
• /
• 2008

Acoustic emission(AE)/Microseimsic(MS) activities are low-energy seismic events associated with a sudden inelastic deformation such as the sudden movement of existing fractures, the generation of new fractures or the propagation of fractures. These events rapidly increase before major failure and happen within a given rock volume and radiate detectable seismic waves. The main difference between AE and MS signals is that the seismic motion frequencies of AE signals are higher than those of MS signals. As the failure of geotechnical structures usually happens as a high velocity and small displacement, it is not easy to determine the precursor and initiation stress level of failure in displacement detection method. To overcome this problem, AE/MS techniques for detection of structure failure and damage have recently adopt in civil engineering. In this study, AE/MS monitoring system, which consist of sensor, data acquisition and operation program, is constructed with domestic technology. To verify and optimize the developed system, we are now carrying out the field application at an underground research laboratory and the developed AE/MS monitoring will be used in detecting of seismic events with various scales.

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

The hypocenter distribution of microseismic events generated by hydraulic fracturing for shale gas development provides essential information for understanding characteristics of fracture network. In this study, we evaluate how inaccurate velocity models influence the inversion results of two widely used location programs, hypoellipse and hypoDD, which are developed based on an iterative linear inversion. We assume that 98 stations are densely located inside the circle with a radius of 4 km and 5 artificial hypocenter sets (S0 ~ S4) are located from the center of the network to the south with 1 km interval. Each hypocenter set contains 25 events placed on the plane. To quantify accuracies of the inversion results, we defined 6 parameters: difference between average hypocenters of assumed and inverted locations, $d_1$; ratio of assumed and inverted areas estimated by hypocenters, r; difference between dip of the reference plane and the best fitting plane for determined hypocenters, ${\theta}$; difference between strike of the reference plane and the best fitting plane for determined hypocenters, ${\phi}$; root-mean-square distance between hypocenters and the best fitting plane, $d_2$; root-mean-square error in horizontal direction on the best fitting plane, $d_3$. Synthetic travel times are calculated for the reference model having 1D layered structure and the inaccurate velocity model for the inversion is constructed by using normal distribution with standard deviations of 0.1, 0.2, and 0.3 km/s, respectively, with respect to the reference model. The parameters $d_1$, r, ${\theta}$, and $d_2$ show positive correlation with the level of velocity perturbations, but the others are not sensitive to the perturbations except S4, which is located at the outer boundary of the network. In cases of S0, S1, S2, and S3, hypoellipse and hypoDD provide similar results for $d_1$. However, for other parameters, hypoDD shows much better results and errors of locations can be reduced by about several meters regardless of the level of perturbations. In light of the purpose to understand the characteristics of hydraulic fracturing, $1{\sigma}$ error of velocity structure should be under 0.2 km/s in hypoellipse and 0.3 km/s in hypoDD.