• Journal of Soil and Groundwater Environment
• /
• v.13 no.1
• /
• pp.1-12
• /
• 2008

The purpose of this study is to characterize the hydrogeochemical characteristics of hot spring waters and to interpret the source of noble gases and the geochemical environment of the hot spring waters distributed along the eastern area of the Korean peninsula. For this purpose, We carried out the chemical, stable isotopic and noble gas isotopic analyses for eleven hot spring water and fourteen hot spring gas samples collected from six hot spring sites. The hot spring waters except the Osaek hot spring water show the pH range of 7.0 to 9.1. However, the Osaek $CO_2$-rich hot spring water shows a weak acid of pH 5.7. The temperature of hot spring waters in the study area ranges from $25.7^{\circ}C$ to $68.3^{\circ}C$. Electrical conductivity of hot spring waters varies widely from 202 to $7,130{\mu}S/cm$. High electrical conductivity (av., $3,890{\mu}S/sm$) by high Na and Cl contents of the Haeundae and the Dongrae hot spring waters indicates that the hot spring waters were mixed with seawater in the subsurface thermal system. The type of hot springs in the viewpoint of dissolved components can be grouped into three types: (1) alkaline Na-$HCO_3$ type including sulfur gas of the Osaek, Baekam, Dukgu and Chuksan hot springs, and (2) saline Na-Cl type of the Haeundae and Dongrae hot springs, and (3) weak acid $CO_2$-rich Na-$HCO_3$ type of Osaek hot spring. Tritium ratios of the Haeundae and the Dongrae hot springs indicate different residence time in their aquifers of older water of $0.0{\sim}0.3$ TU and younger water of $5.9{\sim}8.8$ TU. The ${\delta}^{18}O$ and ${\delta}D$ values of hot spring waters indicate that they originate from the meteoric water, and that the values also reflect a latitude effect according to their locations. $^3He/^4He$ ratios of the hot spring waters except Osaek $CO_2$-rich hot spring water range from $0.1{\times}10^{-6}$ to $1.1{\times}10^{-6}$ which are plotted above the mixing line between air and crustal components. It means that the He gas in hot spring waters was originated mainly from atmosphere and crust sources, and partly from mantle sources. The Osaek $CO_2$-rich hot spring water shows $3.3{\times}10^{-6}$ in $^3He/^4He$ ratio that is 2.4 times higher than those of atmosphere. It provides clearly a helium source from the deep mantle. $^{40}Ar/^{36}Ar$ ratios of hot spring water are in the range of an atmosphere source.

• Journal of Soil and Groundwater Environment
• /
• v.10 no.6
• /
• pp.20-31
• /
• 2005

Three dimensional Monte-Carlo simulations of non-ergodic transport of a lion-reactive solute plume by steady-state groundwater flow under a uniform mean velocity in isotropic heterogeneous aquifers were conducted. The log-normally distributed hydraulic conductivity, K(x), is modeled as a random field. Significant efforts are made to reduce tile simulation uncertainties. Ensemble averages of the second spatial moments of the plume and plume centroid variances were simulated with 1600 Monte Carlo runs for three variances of log K, ${\sigma}_Y^2=0.09,\;0.23$, and 0.46, and three dimensionless lengths of line plume sources normal to the mean velocity. The simulated second spatial moment and the plume centroid variance in longitudinal direction fit well to the first order theoretical results while the simulated transverse moments are generally larger than the first order results. The first order theoretical results significantly underestimated the simulated dimensionless transverse moments for the aquifers of large ${\sigma}_Y^2$ and large dimensionless time. The ergodic condition for the second spatial moments is far from reaching in all cases simulated, and transport In transverse directions may reach ergodic condition much slower than that in longitudinal direction. The evolution of the contaminant transported in a heterogeneous aquifer is not affected by the shape of the initial plume but affected mainly by the degree of the heterogeneity and the size of the initial plume.

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

• Geophysics and Geophysical Exploration
• /
• v.3 no.1
• /
• pp.1-6
• /
• 2000

The 1 $km^2$ area studied is located in Sukchun-ri, Hwasung-koon, the southern part of Kyeonggi-do. Even though this site has been known as a contaminated area caused by seawater intrusions, geophysical and geochemical surveys have never been carried out at the site to determine the extent of the seawater contamination and to investigate whether the seawater intrusion is in progress. The purpose of this study is to determine the extent of seawater contamination and a preferred channel of the seawater intrusion using geophysical methods such as DC resistivity surveys with Schlumberger array and a dipole-dipole array. In order to determine whether the seawater intrusion is in progress in the area, DC resistivity monitoring with Schlumberger array was performed. According to the resistivity map obtained from the inversion of the resistivity data measured with Schlumberger array, the study area is divided into two districts as relatively lowly resistive (less than 30 ohm-m) and highly resistive (more than 30 ohm-m) areas. The distribution of the lowly resistive area is consistent with the distribution of the layer composed of clay minerals, and the resistivity of this layer decreases slowly as approaching to the old seashore. Hydrogeological analysis shows that the clay layer within a distance of about 200 m from the seashore has been already contaminated by sea-water and its electric conductivity is 8 times higher than that of the sand layer covered by the clay layer. According to the results of the 2-dimensional DC resistivity surveys with a dipole-dipole array, there are two preferred channels of the seawater intrusion in the site, and both the channels are in the NW-SE direction from the old seashore. The lowly resistive zone in the southern channel extends to a depth of 80 m. The DC resistivity monitoring with Schlumberger array was carried out along the preferred channel which has the low resistivity Bone (fracture zone) that extended to a depth of 80 m. The time series of apparent resistivity, measured at a distance of 260 m from the old coast line, fluctuates with a period of 12 hours. From these observations, it can be concluded that the seawater intrusion caused by tidal action is still in progress along the fractured zone interpreted by the DC resistivity surveys with a dipole-dipole array.

• Economic and Environmental Geology
• /
• v.56 no.4
• /
• pp.421-433
• /
• 2023

The final disposal of spent nuclear fuel(SNF) from nuclear power plants takes place in a deep geological repository. The metal canister encasing the SNF is made of cast iron and copper, and is engineered to effectively isolate radioactive isotopes for a long period of time. The SNF is further shielded by a multi-barrier disposal system comprising both engineering and natural barriers. The deep disposal environment gradually changes to an anaerobic reducing environment. In this environment, sulfide is one of the most probable substances to induce corrosion of copper canister. Stress-corrosion cracking(SCC) triggered by sulfide can carry substantial implications for the integrity of the copper canister, potentially posing a significant threat to the long-term safety of the deep disposal repository. Sulfate can exist in various forms within the deep disposal environment or be introduced from the geosphere. Sulfate has the potential to be transformed into sulfide by sulfate-reducing bacteria(SRB), and this converted sulfide can contribute to the corrosion of the copper canister. Bentonite, which is considered as a potential material for buffering and backfilling, contains oxidized sulfate minerals such as gypsum(CaSO4). If there is sufficient space for microorganisms to thrive in the deep disposal environment and if electron donors such as organic carbon are adequately supplied, sulfate can be converted to sulfide through microbial activity. However, the majority of the sulfides generated in the deep disposal system or introduced from the geosphere will be intercepted by the buffer, with only a small amount reaching the metal canister. Pyrite, one of the potential sulfide minerals present in the deep disposal environment, can generate sulfates during the dissolution process, thereby contributing to the corrosion of the copper canister. However, the quantity of oxidation byproducts from pyrite is anticipated to be minimal due to its extremely low solubility. Moreover, the migration of these oxidized byproducts to the metal canister will be restricted by the low hydraulic conductivity of saturated bentonite. We have comprehensively analyzed and summarized key research cases related to the presence of sulfates, reduction processes, and the formation and behavior characteristics of sulfides and pyrite in the deep disposal environment. Our objective was to gain an understanding of the impact of sulfates and sulfides on the long-term safety of high-level radioactive waste disposal repository.