• Journal of the Korean Geotechnical Society
• /
• v.20 no.8
• /
• pp.147-158
• /
• 2004

Laboratory dynamic tests are carried out to assess the liquefaction potential of saturated sands in most countries. However, simple results such as the maximum cyclic shear stress and the number of cycles at initial liquefaction are used in the experimental assessment of liquefaction potential, even though various results can be obtained from the dynamic test. In addition, it seemed to be inefficient because more than three dynamic tests with different stress ratio have to be carried out to draw a liquefaction resistance experimental curve. To improve the present assessment method fur liquefaction potential, a new critical resistible characteristic far soil liquefaction is proposed and verified through conventional cyclic triaxial tests with Jumunjin sand. In the proposed method, various experimental data such as effective stress path, stress-strain relationship, and the change of excess pore water pressure can be used in the determination of cumulative plastic shear strains at every 1/4 cycle. Especially, the critical cumulative plastic shear strain to initiate liquefaction can be defined in a specific point called a phase change point in the effective stress path and it can be calculated from a hysteric curve of stress-strain relationship up to this point. Through this research, it is found that the proposed cumulative plastic shear strain can express the dissipated energy to resist dynamic loads and consider the realistic soil dynamic behavior of saturated sands reasonably. It is also found that the critical plastic shear strain can be used as a registible index of soils to represent the critical soil dynamic state, because it seems to include no effect of large deformation.

• Journal of Soil and Groundwater Environment
• /
• v.9 no.1
• /
• pp.79-86
• /
• 2004

Hydraulic conductivity along rock fracture is mainly dependent on fracture geometries such as orientation, aperture, roughness and connectivity. Therefore, it needs to consider fracture geometries sufficiently on a fracture model for a numerical analysis to calculate permeability coefficient in a fracture. This study performed new type of numerical analysis using a homogenization analysis method to calculate permeability coefficient accurately along single fractures with several fracture models that were considered fracture geometries as much as possible. First of all, fracture roughness and aperture variation due to normal stress applied on a fracture were directly measured under a confocal laser scaning microscope (CLSM). The acquired geometric data were used as input data to construct fracture models for the homogenization analysis (HA). Using the constructed fracture models, the homogenization analysis method can compute permeability coefficient with consideration of material properties both in microscale and in macroscale. The HA is a new type of perturbation theory developed to characterize the behavior of a micro inhomogeneous material with a periodic microstructure. It calculates micro scale permeability coefficient at homogeneous microscale, and then, computes a homogenized permeability coefficient (C-permeability coefficient) at macro scale. Therefore, it is possible to analyze accurate characteristics of permeability reflected with local effect of facture geometry. Several computations of the HA were conducted to prove validity of the HA results compared with the empirical equations of permeability in the previous studies using the constructed 2-D fracture models. The model can be classified into a parallel plate model that has fracture roughness and identical aperture along a fracture. According to the computation results, the conventional C-permeability coefficients have values in the range of the same order or difference of one order from the permeability coefficients calculated by an empirical equation. It means that the HA result is valid to calculate permeability coefficient along a fracture. However, it should be noted that C-permeability coefficient is more accurate result than the preexisting equations of permeability calculation, because the HA considers permeability characteristics of locally inhomogeneous fracture geometries and material properties both in microscale and macroscale.

• Journal of the Korean Society for Marine Environment & Energy
• /
• v.12 no.3
• /
• pp.217-226
• /
• 2009

Marine geological storage of $CO_2$ is regarded as one of the most promising options to response climate change. Marine geological storage of $CO_2$ is to capture $CO_2$ from major point sources, to transport to the storage sites and to store $CO_2$ into the marine geological structure such as deep sea saline aquifer. Up to now, process design for this $CO_2$ marine geological storage has been carried out mainly on pure $CO_2$. Unfortunately the captured $CO_2$ mixture contains many impurities such as $N_2$, $O_2$, Ar, $H_2O$, $SO_x$, $H_2S$. A small amount of impurities can change the thermodynamic properties and then significantly affect the compression, purification and transport processes. In order to design a reliable $CO_2$ marine geological storage system, it is necessary to perform numerical process simulation using thermodynamic equation of state. The purpose of the present paper is to compare and analyse the relevant equations of state including PR, PRBM, RKS and SRK equation of state for $CO_2-N_2$ mixture. To evaluate the predictive accuracy of the equation of the state, we compared numerical calculation results with reference experimental data. In addition, optimum binary parameter to consider the interaction of $CO_2$ and $N_2$ molecules was suggested based on the mean absolute percent error. In conclusion, we suggest the most reliable equation of state and relevant binary parameter in designing the $CO_2-N_2$ mixture marine geological storage process.

• Journal of the Korean Institute of Gas
• /
• v.17 no.4
• /
• pp.58-69
• /
• 2013

This paper presents production data analysis for two production wells located in the shale gas field, Canada, with the proper analysis method according to each production performance characteristics. In the case A production well, the analysis was performed by applying both time and superposition time because the production history has high variation. Firstly, the flow regimes were classified with a log-log plot, and as a result, only the transient flow was appeared. Then the area of simulated reservoir volume (SRV) analyzed based on flowing material balance plot was calculated to 180 acres of time, and 240 acres of superposition time. And the original gas in place (OGIP) also was estimated to 15, 20 Bscf, respectively. However, as the area of SRV was not analyzed with the boundary dominated flow data, it was regarded as the minimum one. Therefore, the production forecasting was conducted according to variation of b exponent and the area of SRV. As a result, estimated ultimate recovery (EUR) increased 1.2 and 1.4 times respectively depending on b exponent, which was 0.5 and 1. In addition, as the area of SRV increased from 240 to 360 acres, EUR increased 1.3 times. In the case B production well, the formation compressibility and permeability depending on the overburden were applied to the analysis of the overpressured reservoir. In comparison of the case that applied geomechanical factors and the case that did not, the area of SRV was increased 1.4 times, OGIP was increased 1.5 times respectively. As a result of analysis, the prediction of future productivity including OGIP and EUR may be quite different depending on the analysis method. Thus, it was found that proper analysis methods, such as pseudo-time, superposition time, geomechanical factors, need to be applied depending on the production data to gain accurate results.