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

• 한국지구물리탐사학회:학술대회논문집
• /
• 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.

• The Journal of Engineering Geology
• /
• v.13 no.1
• /
• pp.129-137
• /
• 2003

We interpreted the seismic signal characteristics from crosswell seismic tomography in the underground cavity like abandoned mines. The first arrival time delay and amplitude attenuation showed clearly at the low velocity zone of cavity and fracture. Also ray density decreased by detour of raypath. As a result of the amplitude spectrum analysis of fresh rock and low velocity zone, there were no noticeable differences of the amplitude up to about 1000Hz frequency, but indicated that the one passed around cavity decreased about 7dB at 2000Hz, and 20dB at 3000Hz. It was possible to compare the signal characteristics between two media by extracting the signal data from the fresh rock zone and the underground cavity through the seismic crosswell tomography.

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

• Geophysics and Geophysical Exploration
• /
• v.9 no.1
• /
• pp.30-36
• /
• 2006

Japan's first pilot-scale $CO_2$ sequestration experiment has been conducted in Nagaoka, where 10400 t of $CO_2$ have been injected in an onshore aquifer at a depth of about 1100 m. Among various measurements conducted at the site for monitoring the injected $CO_2$, we conducted time-lapse crosswell seismic tomography between two observation wells to determine the distribution of $CO_2$ in the aquifer by the change of P-wave velocities. This paper reports the results of the crosswell seismic tomography conducted at the site. The crosswell seismic tomography measurements were carried out three times; once before the injection as a baseline survey, and twice during the injection as monitoring surveys. The velocity tomograms resulting from the monitoring surveys were compared to the baseline survey tomogram, and velocity difference tomograms were generated. The velocity difference tomograms showed that velocity had decreased in a part of the aquifer around the injection well, where the injected $CO_2$ was supposed to be distributed. We also found that the area in which velocity had decreased was expanding in the formation up-dip direction, as increasing amounts of $CO_2$ were injected. The maximum velocity reductions observed were 3.0% after 3200 t of $CO_2$ had been injected, and 3.5% after injection of 6200 t of $CO_2$. Although seismic tomography could map the area of velocity decrease due to $CO_2$ injection, we observed some contradictions with the results of time-lapse sonic logging, and with the geological condition of the cap rock. To investigate these contradictions, we conducted numerical experiments simulating the test site. As a result, we found that part of the velocity distribution displayed in the tomograms was affected by artefacts or ghosts caused by the source-receiver geometry for the crosswell tomography in this particular site. The maximum velocity decrease obtained by tomography (3.5%) was much smaller than that observed by sonic logging (more than 20%). The numerical experiment results showed that only 5.5% velocity reduction might be observed, although the model was given a 20% velocity reduction zone. Judging from this result, the actual velocity reduction can be more than 3.5%, the value we obtained from the field data reconstruction. Further studies are needed to obtain more accurate velocity values that are comparable to those obtained by sonic logging.

• Geophysics and Geophysical Exploration
• /
• v.6 no.2
• /
• pp.57-63
• /
• 2003

For the dectection of small cavity in the hard rock, we investigated the feasibility of crosswell travel-time tomography and Kirchhoff migration technique. In travel-time tomography, first arrival anomaly caused by small cavity was investigated by numerical modeling based on the knowledge of actual field information. First arrival delay was very small (<0.125 msec) and detectable receiver offset range was limited to 4m with respect to $1\%$ normalized first arrival anomaly. As a consequence, it was turned out that carefully designed survey array with both sufficient narrow spatial spacing and temporal (<0.03125 msec) sampling were required for small cavity detection. Also, crosswell Kirchhoff migration technique was investigated with both numerical and real data. Stack section obtained by numerical data shows the good cavity image. In crosswell seismic data, various unwanted seismic events such as direct wave and various mode converted waves were alto recorded. To remove these noises und to enhance the diffraction signal, combination of median and bandpass filtering was applied and prestack and stacked migration images were created. From this, we viewed the crosswell migration technique as one of the adoptable method for small cavity detection.

• 한국지구물리탐사학회:학술대회논문집
• /
• 2003.11a
• /
• pp.321-327
• /
• 2003

[ $CO_2$ ] sequestration in oil reservoirs can be one of the most effective strategies for long-term removal of greenhouse gas from atmosphere. This paper presents an advantage of the localized nonlinear approximation of integral equation solutions for inverting crosswell electromagnetic data, which are observed as a part of pilot project of $CO_2$ flooding at the Lost Hills oil field in central California, U.S.A. To monitor the migration of $CO_2$, we have used 2-D cylindrically symmetric and 2.5-D tomographic inversion methods. These two schemes produce nearly the same images if the borehole separation is large compared with the skin depth. However, since the borehole separation is much less than five skin depths in this $CO_2$ injection experiment, the 2.5-D model seems to be more reliable than the 2-D model. In fact, the pre-injection 2.5-D image is more successfully compared with induction logs observed in the two wells than the 2-D model. From the time-lapse crosswell imaging, we can confirm the replacement of brine with $CO_2$ makes a decrease of conductivity.

• Geophysics and Geophysical Exploration
• /
• v.7 no.1
• /
• pp.51-55
• /
• 2004

This paper presents a fast and stable imaging scheme using the localized nonlinear (LN) approximation of integral equation (IE) solutions for inverting electromagnetic data obtained in a crosswell survey. The medium is assumed to be cylindrically symmetric about a source borehole, and to maintain the symmetry a vertical magnetic dipole is used as a source. To find an optimum balance between data fitting and smoothness constraint, we introduce an automatic selection scheme for a Lagrange multiplier, which is sought at each iteration with a least misfit criterion. In this selection scheme, the IE algorithm is quite attractive for saving computing time because Green's functions, whose calculation is a most time-consuming part in IE methods, are repeatedly re-usable throughout the inversion process. The inversion scheme using the LN approximation has been tested to show its stability and efficiency, using both synthetic and field data. The inverted image derived from the field data, collected in a pilot experiment of water-flood monitoring in an oil field, is successfully compared with that derived by a 2.5-dimensional inversion scheme.

• 한국지구과학회:학술대회논문집
• /
• 2003.09a
• /
• pp.140-141
• /
• 2003

• Geophysics and Geophysical Exploration
• /
• v.6 no.3
• /
• pp.138-142
• /
• 2003

Crosswell seismic data were acquired in three sections crossing a tunnel of 3 different types; one was empty, another was ailed by sand, and the other was filled by rock debris. Both the P- and S-wave first arrivals were picked and the traveltime tomography was conducted to generate the P- and S- wave velocity tomograms on the all three sections. Among six tomograms, only one tomogram shows a low velocity zone that can be interpreted as a tunnel image. The tomogram is the P wave velocity image of a section that crosses an empty tunnel. The result of numerical analysis for the spatial resolution of the traveltime tomography was consistent to this finding.