• Geophysics and Geophysical Exploration
• /
• v.9 no.1
• /
• pp.10-18
• /
• 2006

We measured the changes in P-wave velocity that occur when injecting $CO_2$ in gaseous, liquid, and supercritical phases into water-saturated anisotropic sandstones. P-wave velocities were measured in two cylindrical samples of Tako Sandstone, drilled along directions normal and parallel to the bedding plane, using a piezo-electric transducer array system. The velocity changes caused by $CO_2$ injection are typically -6% on average, with maximum values about -16% for the case of supercritical $CO_2$ injection. P-wave velocity tomograms obtained by the differential arrival-time method clearly show that $CO_2$ migration behaviour is more complex when $CO_2$ flows normal to the bedding plane than when it flows parallel to bedding. We also found that the differences in P-wave velocity images were associated both with the $CO_2$ phases and with heterogeneity of pore distribution in the rocks. Seismic images showed that the highest velocity reduction occurred for supercritical $CO_2$ injection, compared with gaseous or liquid $CO_$ injection. This result may justify the use of the seismic method for $CO_2$ monitoring in geological sequestration.

• Geophysics and Geophysical Exploration
• /
• v.9 no.1
• /
• pp.37-43
• /
• 2006

In 2006, the Cooperative Research Centre for Greenhouse Gas Technologies (CO2CRC) plans to undertake (subject to receiving the necessary approvals) a Pilot program for $CO_2$ storage within a depleted gas reservoir. The Otway Basin Pilot Program (OBPP) aims to demonstrate that subsurface $CO_2$ storage is both economically and environmentally sustainable in Australia. This will be the first $CO_2$ storage program in the world to utilise a depleted gas reservoir and, hence, the experience gained will be a valuable addition to the range of international $CO_2$ storage programs that are underway or being planned. A key component of the OBPP is the design of an appropriate geophysical monitoring strategy that will allow the subsurface migration of the $CO_2$ plume to be tracked and to verify that containment has been successful. This paper presents the results from modelling the predicted gravity response to $CO_2$ injection into the Otway Basin reservoir, where the goal was to determine minimum volumes of $CO_2$ that may be detectable using non-seismic geophysical techniques. Modelling results indicate that gravity measurements at 10 m spacing within the existing observation well and the planned $CO_2$ injection well would provide excellent vertical resolution, even for the smallest $CO_2$ volume modelled (10000 tonnes), but resolving the lateral extent of the plume would not be possible without additional wells at closer spacing.

• Geophysics and Geophysical Exploration
• /
• v.18 no.3
• /
• pp.97-104
• /
• 2015

KIGAM has surveyed most of the Korean territory since 1982 using airborne gamma-ray spectrometry, and complete the nationwide scale map in the near future. However, since the duration of survey is too long and the conditions of survey is not consistent, the data does not have physical consistency. In addition, the window counts (count/sec) were recorded instead of potassium, uranium and thorium radioelement concentrations. Thus, the data could not be registered to the International Atomic Energy Agency (IAEA) radioelement datum. This limits the usefulness of the data and it is not possible to easily combine surveys into regional compilations or make quantitative interpretations between different survey areas. To solve these problems, we undertook a test baseline survey over Jincheon-Eumseong area, to level the different two sets of data and to map radioelement concentrations. This survey confirms to IAEA radioelement baseline. The method and procedures of data leveling prepared by this study improve the usefulness and usability of the radiometric data, and make it enable to compile the nationwide scale radioelement concentration maps.

• Geophysics and Geophysical Exploration
• /
• v.18 no.3
• /
• pp.154-161
• /
• 2015

For efficiently designing neutron induced gamma spectroscopy sonde, Monte Carlo simulation is employed to understand a dominant location of thermal neutron and classify the formation elements from the energy peak of capture gamma spectrum. A pulsed neutron generator emitting 14 MeV neutron particles was used as a source, and flux of thermal neutron was calculated from the twelve detectors arranged at each 10 cm intervals from the source. Design for reducing borehole effects using shielding materials was also applied to numerical sonde model. Moreover, principal elements and quantities of numerical earth models were verified through the energy spectrum analysis of capture gamma detected from a gamma detector. These results can help to enhance the signal-to-noise ratio, and determine an optimal placement of capture gamma detectors of neutron induced gamma spectroscopy sonde.

• Geophysics and Geophysical Exploration
• /
• v.18 no.3
• /
• pp.143-153
• /
• 2015

Polymer flooding for enhancing hydrocarbon production injects into a reservoir polymer solution that is viscous. It is very important to monitor the behavior pattern of the polymer solution in order to evaluate the effectiveness of polymer flooding. To monitor the distribution of polymer solution and thus fluid substitution within the reservoir, we first construct seismic and resistivity rock physics models (RPMs), which are functions of reservoir parameters such as rocks and type of fluid, fluid saturation. For the seismic and resistivity RPMs, responses of seismic and electromagnetic (EM) tomography are numerically simulated as polymer injection, using two dimensional (2D) staggered-grid finite difference elastic modeling and 2.5D finite element EM modeling algorithms, respectively. In constructing RPM for EM tomography, three different reservoir rocks are considered: clean-sand, dispersed shale-sand, and sand-shale lamination rocks. The polymer solution is assumed to have 2 wt% of polymer as normally generated, while water is freshwater or saltwater. Further, neutron logging is also considered to check its sensitivity to polymer flooding. The techniques discussed in the paper are important in monitoring not only hydrocarbon but also geothermal reservoirs.

• Geophysics and Geophysical Exploration
• /
• v.18 no.3
• /
• pp.125-132
• /
• 2015

The water saturation ($S_w$), which is very important to estimate hydrocarbon reserves in the reservoir, has been determined from resistivity index (RI) by using the Archie's formula. However, in many cases, it has been reported that n is not constant for a given formation and it could be varied with water saturation. In addition, it frequently happens that the line obtained by linear regression analysis on log-log scale does not pass through the origin. In order to overcome these drawbacks, we suggested a modified Archie's formula, which can handle almost all the RI vs. $S_w$ cross-plots whether the trend is straight or curved and whether it passes through the origin or not. We also demonstrated that how to determine conductivity exponent ${\mu}$, critical water saturation $S_c$, and saturation distribution factor b in the laboratory to use the modified Archie's formula. Since the modified Archie's formula takes into account pore structure, pore water distribution, and wettability of reservoir such as clean sandstone and carbonate rocks, it might improve field applicability.

• Geophysics and Geophysical Exploration
• /
• v.10 no.1
• /
• pp.1-18
• /
• 2007

Non-linear elastic wavefield inversion is a powerful method for estimating elastic parameters for physical constraints that determine subsurface rock and properties. Here, I introduce six elastic-wave velocity models by reconstructing elastic-wave velocity variations from real data and a 2D elastic-wave velocity model. Reflection seismic data information is often decoupled into short and long wavelength components. The local search method has difficulty in estimating the longer wavelength velocity if the starting model is far from the true model, and source frequencies are then changed from lower to higher bands (as in the 'frequency-cascade scheme') to estimate model elastic parameters. Elastic parameters are inverted at each inversion step ('simultaneous mode') with a starting model of linear P- and S-wave velocity trends with depth. Elastic parameters are also derived by inversion in three other modes - using a P- and S-wave velocity basis $('V_P\;V_S\;mode')$; P-impedance and Poisson's ratio basis $('I_P\;Poisson\;mode')$; and P- and S-impedance $('I_P\;I_S\;mode')$. Density values are updated at each elastic inversion step under three assumptions in each mode. By evaluating the accuracy of the inversion for each parameter set for elastic models, it can be concluded that there is no specific difference between the inversion results for the $V_P\;V_S$ mode and the $I_P$ Poisson mode. The same conclusion is expected for the $I_P\;I_S$ mode, too. This gives us a sound basis for full wavelength elastic wavefield inversion.

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

In this study, we propose a joint inversion method, using genetic algorithms, to estimate an S-wave velocity structure for deep sedimentary layers from receiver functions and surface-wave phase velocity observed at several sites. The method takes layer continuity over a target area into consideration by assuming that each layer has uniform physical properties, especially an S-wave velocity, at all the sites in a target area in order to invert datasets acquired at different sites simultaneously. Numerical experiments with synthetic data indicate that the proposed method is effective in reducing uncertainty in deep structure parameters when modelling only surface-wave dispersion data over a limited period range. We then apply the method to receiver functions derived from earthquake records at one site and two datasets of Rayleigh-wave phase velocity obtained from microtremor array surveys performed in central Tokyo, Japan. The estimated subsurface structure is in good agreement with the results of previous seismic refraction surveys and deep borehole data. We also conclude that the proposed method can provide a more accurate and reliable model than individual inversions of either receiver function data only or surface-wave dispersion data only.

• Geophysics and Geophysical Exploration
• /
• v.10 no.1
• /
• pp.29-36
• /
• 2007

Recent seismic surveys have shown that the presence of methane hydrate (MH) in sediments has significant influence on seismic attenuation. I have used vertical seismic profile (VSP) data from a Nankai Trough exploratory well, offshore Tokai in central Japan, to estimate compressional attenuation in MH-bearing sediments at seismic frequencies of 30-110 Hz. The use of two different measurement methods (spectral ratio and centroid frequency shift methods) provides an opportunity to validate the attenuation measurements. The sensitivity of attenuation analyses to different depth intervals, borehole irregularities, and different frequency ranges was also examined to validate the stability of attenuation estimation. I found no significant compressional attenuation in MH-bearing sediments at seismic frequencies. Macroscopically, the peaks of highest attenuation in the seismic frequency range correspond to low-saturation gas zones. In contrast, high compressional attenuation zones in the sonic frequency range (10-20 kHz) are associated with the presence of methane hydrates at the same well locations. Thus, this study demonstrated the frequency-dependence of attenuation in MH-bearing sediments; MH-bearing sediments cause attenuation in the sonic frequency range rather than the seismic frequency range As a possible reason why seismic frequencies in the 30-110 Hz range were not affected in MH-bearing sediments, I point out the effect of thin layering of MH-bearing zones.

• Geophysics and Geophysical Exploration
• /
• v.10 no.1
• /
• pp.50-59
• /
• 2007

The advent of towed geo-electrical array surveying on water and land has resulted in datasets of magnitude approaching that of airborne electromagnetic surveying and most suited to 1D inversion. Robustness and complete automation is essential if processing and reliable interpretation of such data is to be viable. Sharp boundaries such as river beds and the top of saline aquifers must be resolved so use of smoothness constraints must be minimised. Suitable inversion algorithms must intelligently handle low signal-to-noise ratio data if conductive basement, that attenuates signal, is not to be misrepresented. A noise-level aware inversion algorithm that operates with one elastic thickness layer per electrode configuration has been coded. The noise-level aware inversion identifies if conductive basement has attenuated signal levels so that they are below noise level, and models conductive basement where appropriate. Layers in the initial models are distributed to span the effective depths of each of the geo-electric array quadrupoles. The algorithm works optimally on data collected using geo-electric arrays with an approximately exponential distribution of quadrupole effective depths. Inversion of data from arrays with linear electrodes, used to reduce contact resistance, and capacitive-line antennae is plausible. This paper demonstrates the effectiveness of the algorithm using theoretical examples and an example from a salt interception scheme on the Murray River, Australia.