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

In this paper, we will introduce rock physics modeling technique, which interrelate reservoir properties with seismic properties, and apply the technique to the Donghae-1 gas reservoir. From well-log data analysis, we obtained velocityporosity (Vp-$\phi$) relations for each formation. These relations can used to predict porosity from seismic data. In addition, we analyzed permeability data, which were obtained from core measurements and computational rock physics simulations. We then obtained permeability-porosity ($\kappa-\phi$) relations. Combining $\kappa-\phi$ with Vp-$\phi$ relations, we finally present quantitative Vp-$\kappa$ relations. As to Vp-$\phi$ modeling, we found that the degree of diagenesis and clay contents increase with depth. As to Vp-$\kappa$ relations, though \kappa-\phi relations are almost identical for all formations, we could obtain distinct Vp-$\kappa$ relations due to Vp-$\phi$ variations. In conclusion, the rock physics modeling, which bridges between seismic properties and reservoir properties, can be a very robust tool for quantitative reservoir characterization with less uncertainty.

• 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.13 no.1
• /
• pp.80-87
• /
• 2010

This study was conducted to develop a reservoir modelling workflow to reproduce the heterogeneous distribution of effective permeability that impacts on the performance of SAGD (Steam Assisted Gravity Drainage), the in-situ bitumen recovery technique in the Athabasca Oil Sands. Lithologic facies distribution is the main cause of the heterogeneity in bitumen reservoirs in the study area. The target formation consists of sand with mudstone facies in a fluvial-to-estuary channel system, where the mudstone interrupts fluid flow and reduces effective permeability. In this study, the lithologic facies is classified into three classes having different characteristics of effective permeability, depending on the shapes of mudstones. The reservoir modelling workflow of this study consists of two main modules; facies modelling and permeability modelling. The facies modelling provides an identification of the three lithologic facies, using a stochastic approach, which mainly control the effective permeability. The permeability modelling populates mudstone volume fraction first, then transforms it into effective permeability. A series of flow simulations applied to mini-models of the lithologic facies obtains the transformation functions of the mudstone volume fraction into the effective permeability. Seismic data contribute to the facies modelling via providing prior probability of facies, which is incorporated in the facies models by geostatistical techniques. In particular, this study employs a probabilistic neural network utilising multiple seismic attributes in facies prediction that improves the prior probability of facies. The result of using the improved prior probability in facies modelling is compared to the conventional method using a single seismic attribute to demonstrate the improvement in the facies discrimination. Using P-wave velocity in combination with density in the multiple seismic attributes is the essence of the improved facies discrimination. This paper also discusses sand matrix porosity that makes P-wave velocity differ between the different facies in the study area, where the sand matrix porosity is uniquely evaluated using log-derived porosity, P-wave velocity and photographically-predicted mudstone volume.

• Geophysics and Geophysical Exploration
• /
• v.8 no.1
• /
• pp.1-6
• /
• 2005

Analogue physical modelling using granular materials (i.e., sandbox experiments) has been applied with great success to a number of geological problems at various scales. Such physical experiments can also be simulated numerically with the Discrete Element Method (DEM). In this study, we apply the DEM simulation to the collision between the Indian subcontinent and the Eurasian Plate, one of the most significant current tectonic processes in the Earth. DEM simulation has been applied to various kinds of dynamic modelling, not only in structural geology but also in soil mechanics, rock mechanics, and the like. As the target of the investigation is assumed to be an assembly of many tiny particles, DEM simulation makes it possible to treat an object with large and discontinuous deformations. However, in DEM simulations, we often encounter difficulties when we examine the validity of the input parameters, since little is known about the relationship between the input parameters for each particle and the properties of the whole assembly. Therefore, in our previous studies (Yamada et al.,2002a,2002b,2002c), we were obliged to tune the input parameters by trial and error. To overcome these difficulties, we introduce a numerical biaxial test with the DEM simulation. Using the results of this numerical test, we examine the validity of the input parameters used in the collision model. The resulting collision model is quite similar to the real deformation observed in eastern Asia, and compares well with GPS data and in-situ stress data in eastern Asia.

• Journal of the Korean Society of Groundwater Environment
• /
• v.5 no.4
• /
• pp.192-202
• /
• 1998

We investigated the geochemistry and environmental isotopes of granite-bedrock groundwater in the Yeongcheon diversion tunnel which is located about 300 m below the land surface. The hydrochemistry of groundwaters belongs to the Ca-HCO$_3$type, and is controlled by flow systems and water-rock interaction in the flow conduits (fractures). The deuterium and oxygen-18 data are clustered along the meteoric water line, indicating that the groundwater are commonly of meteoric water origin and are not affected by secondary isotope effects such as evaporation and isotope exchange. Tritium data show that the groundwaters were mostly recharged before pre-thermonuclear period and have been mixed with younger surface water flowing down rapidly into the tunnel along fractured zones. Based on the mass balance and reaction simulation approaches, using both the hydrochemistry of groundwater and the secondary mineralogy of fracture-filling materials, we have modeled the low-temperature hydrogeochemical evolution of groundwater in the area. The results of geochemical simulation show that the concentrations of Ca$\^$2＋/, Na$\^$＋/ and HCO$_3$and pH of waters increase progressively owing to the dissolution of reactive minerals in flow paths. The concentrations of Mg$\^$2＋/ and K$\^$＋/ frist increase with the dissolution, but later decrease when montmorillonite and illitic material are precipitated respectively. The continuous adding of reactive minerals, namely the progressively larger degrees of water/rock interaction, causes the formation of secondary minerals with the following sequence: first hematite, then gibbsite, then kaolinite, then montmorillonite, then illtic material, and finally microcline. During the simulation all the gibbsite is consumed, kaolinite precipitates and then the continuous reaction converts the kaolinite to montmorillonite and illitic material. The reaction simulation results agree well with the observed, water chemistry and secondary mineralogy, indicating the successful applicability of this simulation technique to delineate the complex hydrogeochemistry of bedrock groundwaters.

• The Journal of Engineering Geology
• /
• v.31 no.4
• /
• pp.507-522
• /
• 2021

Specific survey objectives often cannot be met using only one geophysical method, as each method's results are influenced by the specific physical properties of subsurface materials. In particular, areas susceptible to geological hazards require investigation using more than one method in order to reduce risks to life and property. Instead of analyzing the results from each method separately, this work develops a four-quadrant criterion for classifying areas of levees as safe or weak. The assessment is based on statistically determined thresholds of seismic velocity (P-wave velocity from seismic refraction and S-wave velocity from multichannel analysis of surface waves) and electrical resistivity. Thresholds are determined by subtracting the standard deviation from the mean during performance testing of this correlation technique applied to model data of four horizontal and inclined fracture zones. Compared with results from the crossplot of resistivity and P-wave velocity, crossplot analysis using resistivity and S-wave velocity data provides more reliable information on the soil type, ground stiffness, and lithological characteristics of the levee system. A loose and sandy zone (represented by low S-wave velocity and high resistivity) falling within the second quadrant is interpreted to be a weak zone. This interpretation is well supported by the N values from standard penetrating test for the central core.

• Economic and Environmental Geology
• /
• v.50 no.6
• /
• pp.487-495
• /
• 2017

A granite rock body, called 'Bukbawi', located on a mountaineering trail at Mt. Palgong Provincial Park is popular among the public because it resembles a percussion instrument. If someone hits the specific surface area of this rock body, people can hear drum-like sound. Such phenomenon may be geologically associated with exfoliation process of the granite body or miarolitic cavity developed after gasses escaped during formation of granite. To understand better the inner structure causing drum-like sound, we carried out a non-destructive ground-penetrating radar survey. In this study, as our primary target is very close to the surface, we utilized 1 GHz antennas to produce high-resolution near-surface images. In order to construct 3-D internal images, the measurements were conducted along a pre-defined grid. The processed radargrams revealed that the locations associated with 'drum' sound coincide with strong reflections. In addition, both reflection patterns of fracture and cavity were observed. To further quantify the observed reflections, we simulated GPR scans from a synthetic fracture in a granite body, filled with different materials. The simulated results suggest that both exfoliation process and miarolitic cavity may have contributed to the 'drum' phenomena. Furthermore, the radargrams showed a well-developed cavity signature where two major reflection planes were crossed. Thus, our study is an example of non-destructive geophysical studies that can promote Earth Science in the broader community by examining geological structures attracting the public.