• 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.8 no.2
• /
• pp.145-155
• /
• 2005

A two-dimensional (2-D) interpretation of MT data has been performed for the purpose of fracture detection for geothermal development. Remote stations have been operated in Kyushu, Japan (480 km apart) as well as in Korea (60 km and 165 km apart in 2002 and 2003 data set, respectively). Apparent resistivity and phase curves calculated by remote processing with the Japan remote data showed enough quality for 2-D inversion for the whole frequency range. Remote reference processing with Korea remote reference data also showed quite good continuity in apparent resistivity and phase curves except some noisy frequency bands; around the power frequency, 60 Hz, and around the dead band $10^{-1}Hz\;Hz\;\~1\;Hz$, where the natural EM signal is known to be very weak. Even though the subsurface showed severe three-dimensional (3-D) characteristics in the survey area so that 2-D inversion by itself could not give enough information for deep geological structures, the 2-D inversion for the 5 survey lines showed several common features. The conductive semi-consolidate mudstone layer is dipping from north to south (about 500 m depth on the south and 200 m on the north most part of the survey area). The boundary between the low (L-2) and high (H-2) resistivity anomalies can be thought as a major fault with strike $N15^{\circ}E$, passing through the sites 206, 112 and 414. The shallow (< 1 km) conductive anomalies (L-4) seem to be fracture zones having strike E-W (at site 105) and $N60^{\circ}W$ (at site 434). And there exists a conductive layer in the western and west-southern part of the survey area in the depth below $2\~3\;km$, for which further investigation is to be needed.

• Geophysics and Geophysical Exploration
• /
• v.7 no.3
• /
• pp.164-173
• /
• 2004

To delineate geothermal water movement at the Pohang geothermal development site, Self-Potential (SP) survey and monitoring were carried out during pumping tests. Before drilling, background SP data have been gathered to figure out overall potential distribution of the site. The pumping test was performed in two separate periods: 24 hours in December 2003 and 72 hours in March 2004. SP monitoring started several days before the pumping tests with a 128-channel automatic recording system. The background SP survey showed a clear positive anomaly at the northern part of the boreholes, which may be interpreted as an up-flow Bone of the deep geothermal water due to electrokinetic potential generated by hydrothermal circulation. The first and second SP monitoring during the pumping tests performed to figure out the fluid flow in the geothermal reservoir but it was not easy to see clear variations of SP due to pumping and pumping stop. Since the area is covered by some 360 m-thick tertiary sediments with very low electrical resistivity (less than 10 ohm-m), the electrokinetic potential due to deep groundwater flow resulted in being seriously attenuated on the surface. However, when we compared the variation of SP with that of groundwater level and temperature of pumping water, we could identify some areas responsible to the pumping. Dominant SP changes are observed in the south-west part of the boreholes during both the preliminary and long-term pumping periods, where 3-D magnetotelluric survey showed low-resistivity anomaly at the depth of $600m\~1,000m$. Overall analysis suggests that there exist hydraulic connection through the southwestern part to the pumping well.

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

Recently, the studies about rock physics model (RPM) in shale reservoir are widely performed. In shale reservoir, the degree of the maturity can be estimated by kerogen and GOR (Gas-Oil Ratio). The researches on the rock physics model of shale reservoir with the amount of kerogen have been actively carried out but not with GOR. Thus, in this study, we analyzed the changes in seismic velocity and density, and AVO (Amplitude Variation with Offset) response depending on changes in GOR and the amount of kerogen. Since the shale consists of plate-like particles, it has vertical transverse isotropy (VTI). Therefore we estimated the seismic velocity and density by using Backus averaging method and analyzed AVO responses based on these estimated properties. The results of analysis showed that the changes in the velocity with the GOR variation are small but the velocity changes with the variation in kerogen amount are relatively larger. In case, GOR 180 (Litre/Litre) which is boundary between heavy oil and light oil, when volume fraction of kerogen increased from 5% to 35%, the P-wave velocity normal to the layering increased 51%. That is, it helps estimating maturity of kerogen through the velocity. Meanwhile, when rates of oil-gas mixture are large, the effect of GOR variation on the velocity change became larger. In case volume fraction of kerogen is 5%, the P-wave velocity normal to the layering was estimated $1.46km/s^2$ in heavy oil (GOR 40) but $1.36km/s^2$ in light oil (GOR 300). The AVO responses analysis showed class 4 regardless of the GOR and amount of kerogen because variation of poisson's ratio is small. Therefore, shale reservoir has possibility to have class 4.

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

• Geophysics and Geophysical Exploration
• /
• v.11 no.1
• /
• pp.26-33
• /
• 2008

The analysis of wide-angle seismic reflection and refraction data plays an important role in lithospheric-scale crustal structure study. However, it is extremely difficult to develop an appropriate velocity structure model directly from the observed data, and we have to improve the structure model step by step, because the crustal structure analysis is an intrinsically non-linear problem. There are several subjective processes in wide-angle crustal structure modelling, such as phase identification and trial-and-error forward modelling. Because these subjective processes in wide-angle data analysis reduce the uniqueness and credibility of the resultant models, it is important to reduce subjectivity in the analysis procedure. From this point of view, we describe two software tools, PASTEUP and MODELING, to be used for developing crustal structure models. PASTEUP is an interactive application that facilitates the plotting of record sections, analysis of wide-angle seismic data, and picking of phases. PASTEUP is equipped with various filters and analysis functions to enhance signal-to-noise ratio and to help phase identification. MODELING is an interactive application for editing velocity models, and ray-tracing. Synthetic traveltimes computed by the MODELING application can be directly compared with the observed waveforms in the PASTEUP application. This reduces subjectivity in crustal structure modelling because traveltime picking, which is one of the most subjective process in the crustal structure analysis, is not required. MODELING can convert an editable layered structure model into two-way traveltimes which can be compared with time-sections of Multi Channel Seismic (MCS) reflection data. Direct comparison between the structure model of wide-angle data with the reflection data will give the model more credibility. In addition, both PASTEUP and MODELING are efficient tools for handling a large dataset. These software tools help us develop more plausible lithospheric-scale structure models using wide-angle seismic data.

• Geophysics and Geophysical Exploration
• /
• v.11 no.1
• /
• pp.52-59
• /
• 2008

We have developed a new deep-towed marine DC resistivity survey system. It was designed to detect the top boundary of the methane hydrate zone, which is not imaged well by seismic reflection surveys. Our system, with a transmitter and a 160-m-long tail with eight source electrodes and a receiver dipole, is towed from a research vessel near the seafloor. Numerical calculations show that our marine DC resistivity survey system can effectively image the top surface of the methane hydrate layer. A survey was carried out off Joetsu, in the Japan Sea, where outcrops of methane hydrate are observed. We successfully obtained DC resistivity data along a profile ${\sim}3.5\;km$ long, and detected relatively high apparent resistivity values. Particularly in areas with methane hydrate exposure, anomalously high apparent resistivity was observed, and we interpret these high apparent resistivities to be due to the methane hydrate zone below the seafloor. Marine DC resistivity surveys will be a new tool to image sub-seafloor structures within methane hydrate zones.

• Geophysics and Geophysical Exploration
• /
• v.11 no.1
• /
• pp.1-14
• /
• 2008

We present a method for interpreting seismic records with arrivals and waveforms having characteristics which could be generated by extremely inhomogeneous velocity structures, such as non-typical oceanic crust, decollement at subduction zones, and seamounts in oceanic regions, by comparing them with synthetic waveforms. Recent extensive refraction and wide-angle reflection surveys in oceanic regions have provided us with a huge number of high-resolution and high-quality seismic records containing characteristic arrivals and waveforms, besides first arrivals and major reflected phases such as PmP. Some characteristic waveforms, with significant later reflected phases or anomalous amplitude decay with offset distance, are difficult to interpret using only a conventional interpretation method such as the traveltime tomographic inversion method. We find the best process for investigating such characteristic phases is to use an interactive interpretation method to compare observed data with synthetic waveforms, and calculate raypaths and traveltimes. This approach enables us to construct a reasonable structural model that includes all of the major characteristics of the observed waveforms. We present results here with some actual observed examples that might be of great help in the interpretation of such problematic phases. Our approach to the analysis of waveform characteristics is endorsed as an innovative method for constructing high-resolution and high-quality crustal structure models, not only in oceanic regions, but also in the continental regions.

• Geophysics and Geophysical Exploration
• /
• v.20 no.3
• /
• pp.146-162
• /
• 2017

Multibeam and seismic data in the northwestern part of the Ulleung Basin were analyzed to study stratigraphy and evolutionary history of submarine canyon. A detailed analysis reveals that the sedimentary sequences in this area consist of four stratigraphic units separated by erosional unconformities. On the continental slope, these units are dominated by well-stratified facies with some slope failures, whereas these units show well-stratified and chaotic facies toward the basin floor. Generally, the sediment thickness is relatively thin on the slope, whereas thick sediment accumulation occurs on the base of slope and basin floor. Based on seismic characteristics and distribution, the deposition of each units are well correlated with the evolutionary history of the submarine canyon. Unit 1 directly overlying the acoustic basement has thin sediment layer on the slope, whereas its thickness gradually increase toward the basin floor. Compared to other units, Unit 2 is relatively thick accumulations on the slope and contains some slope failures related to faults systems. The mass transport sediments due to slope failures, mainly deposited on the base of slope as a submarine fan. The width and depth of submarine canyon increase due to dominant of the erosional process rather than the sediment deposition. Unit 3 is thin accumulation on the slope around the submarine canyon. Toward the basin floor, its thickness gradually increases. Unit 4 is characterized by thin layers including slides and slumps on the slope, whereas it formed thick accumulations at the base of slope as a submarine fan. The increase in the width and depth of submarine canyon results from the dominant of the erosional process and slope failures around the submarine canyon. Consequently, the formation of sedimentary units combined with the development of submarine canyon in this area is largely controlled by the amounts of sediment supply originated from slope failures, regional tectonic effects and sea-level fluctuations.

• Economic and Environmental Geology
• /
• v.51 no.3
• /
• pp.249-277
• /
• 2018

This study aims to review the research trend of articles published in the journal of Economic and Environmental Geology(EEG) over a period of 50 years from 1968 to 2017. Based on the literature review, Research theme as analytical criteria was adopted. 1,886 academic articles were analyzed and classified by decade, by change period of journal name, and by volume(in sets of 10 volumes). Results indicate that from the foundation of journal to the present, the article in the field of petrology and geochemistry covers the highest proportion as 30.4%, the next best things are in the fields of metallic ore deposit(18.3%) and geophysics and exploration(14.9%). Based on the change period of journal name, the field of petrology and geochemistry have the most obvious difference in the number of articles. Until 1993 as before the change of journal name from Mining Geology to Economic and Environmental Geology, the coverage percentage in the field of petrology and geochemistry is 21.5%(N=121). However, after 1994 changed to Economic and Environmental Geology, the coverage percentage in the field of petrology and geochemistry is 34.1%(N=452). According to the analytical result in case of volume and research theme, the number of articles in the field of petrology and geochemistry have the abrupt increase from the volume number in 21 to 30. However, the number of articles in the field of metallic and nonmetallic ore deposit have the decreasing tendency toward the recent volume. In order to explore again the 50 years history of Economic and Environmental Geology, we mainly analyzed the articles quantitatively using the criteria as research theme. Based on these analyses, this review aims to offer the beginning for the exploration of the development subject in the future.