• Geomechanics and Engineering
• /
• v.23 no.4
• /
• pp.327-338
• /
• 2020

The role of precise prediction of subsurface fluids and discrimination among them cannot be ignored in reservoir characterization and petroleum prospecting. A suitable rock physics model should be build for the extraction of valuable information form seismic data. The main intent of current work is to present a rock physics model to analyze the characteristics of seismic wave propagating through a cracked porous rock saturated by a three phase fluid. Furthermore, the influence on wave characteristics due to variation in saturation of water, oil and gas were also analyzed for oil and water as wet cases. With this approach the objective to explore wave attenuation and dispersion due to wave induce fluid flow (WIFF) at seismic and sub-seismic frequencies can be precisely achieved. We accomplished our proposed approach by using BISQ equations and by applying appropriate boundary conditions to incorporate heterogeneity due to saturation of three immiscible fluids forming a layered system. To authenticate the proposed methodology, we compared our results with White's mesoscopic theory and with the results obtained by using Biot's poroelastic relations. The outcomes reveals that, at low frequencies seismic wave characteristics are in good agreement with White's mesoscopic theory, however a slight increase in attenuation at seismic frequencies is because of the squirt flow. Moreover, our work crop up as a practical tool for the development of rock physical theories with the intention to identify and estimate properties of different fluids from seismic data.

• Tunnel and Underground Space
• /
• v.32 no.2
• /
• pp.144-159
• /
• 2022

The present study developed a sequential approach-based numerical simulator for modeling coupled thermal-hydrological-mechanical (THM) processes in the ground and investigated the computational performance of the coupling analysis algorithm. The present sequential approach linked the two different solvers: an open-source numerical code, OpenGeoSys for solving the thermal and hydrological processes in porous media and a commercial code, FLAC3D for solving the geomechanical response of the ground. A benchmark test of the developed simulator was carried out using a THM problem where an analytical solution is given. The benchmark problem involves the coupled behavior (variations in temperature, pore pressure, stress, and deformation with time) of a fully saturated porous medium which is subject to a point heat source. The results of the analytical solution and numerical simulation were compared and the validity of the numerical simulator was investigated.

• Water for future
• /
• v.27 no.3
• /
• pp.107-113
• /
• 1994

The analysis of Richards eq. requires data of the soil water retention function and the unsaturated hydraulic conductivity. The soil water retention function was measured through the use of the developed apparatus and the saturated hydraulic conductivity was measured by the constant head method for each soil sample corresponding to the A, B, C types of SCS. In order to obtain one water retention function and one unsaturated hydraulic conductivity which represent each soil group, van Genuchten's eq. and Mualem's pore-structure model was chosen respectively. Parameters of van Genuchten's eq. are estimated for each soil group using data obtained in the experiments, and estimated values give a basis to analyze the unsaturated flow in the non-measured region efficiently.

• Journal of Soil and Groundwater Environment
• /
• v.21 no.1
• /
• pp.104-110
• /
• 2016

The transport of silver nanoparticles (AgNPs) was investigated through a column packed with sand. A series of column experiments were carried out to evaluate the effect of ionic strength (IS), pH, electrolyte type and clay mineral on mobility of polyvinylpyrrolidone-coated silver nanoparticles (PVP-AgNPs). The deposition of PVP-AgNPs was increased with increasing solution ionic strength and decreasing pH. Furthermore, the depositon of PVP-AgNPs was affected by the electrolyte type (NaCl vs. NaNO₃) and was shown to be greater at NaNO₃ solution. Also, the transport of PVP-AgNPs was greatly increased after the pre-deposition of clay particles on sand. Our results suggest that various environmental factors can influence the mobility of PVP-AgNPs in soil-groundwater systems and should be carefully considered in assessing their environmental risks.

• Geomechanics and Engineering
• /
• v.29 no.1
• /
• pp.53-63
• /
• 2022

Simplified analytical solutions are developed for the dynamic analyses of an axially loaded pile foundation embedded in a transverse-isotropic, fluid-filled, poro-visco-elastic soil with rigid substratum. The pile is modeled as a viscoelastic Rayleigh-Love rod, while the surrounding soil is regarded as a transversely isotropic, liquid-saturated, viscoelastic, porous medium of which the mechanical behavior is represented by the Boer's poroelastic media model and the fractional derivative model. Upon the separation of variables, the frequency-domain responses for the impedance function of the pile top, and the vertical displacement and the axial force along the pile shaft are gained. Then by virtue of the convolution theorem and the inverse Fourier transform, the time-domain velocity response of the pile head is derived. The presented solutions are validated, compared to the existing solution, the finite element model (FEM) results, and the field test data. Parametric analyses are made to show the effect of the soil anisotropy and the excitation frequency on the pile-soil dynamic responses.

• The Journal of Engineering Geology
• /
• v.14 no.3 s.40
• /
• pp.301-311
• /
• 2004

This study is suggested a new dielectric tracer test method to understand geological structure of porous media and groundwater flow to use the dielectric constant which is one of electrical special quality of various geological materials. To measure their parameters, tracer material is made an ethanol mixing liquid(EML) having a same specific gravity of water. Also, soil materials are prepared a dielectric tracer test using the FDR system that could measure dielectric constant for saturated standard sand and river sand layers which have different initial porosity. To compare with their results, we discussed with the concentration variation of saline water having a saline concentration $3\%$ which is general tracer material by using the electro multi-meter system in the laboratory or field test. In two tracer experiment results, EML tracer test could confirm definitely EML concentration variation from each saturated soil layer as standard and river sands. However, tracer test of saline water $3\%$ concentration could not confirm permeating movement of water by degree of salinity change because these are settled at lower part column in a whole column area continuously. These causes are that specific gravity of saline water is heavier than water. That is, it could know that deposition of saline water is composed of lower part of soil column continuously independently of the direction of water into saturated soil material.

• Journal of the Korean Society of Safety
• /
• v.11 no.4
• /
• pp.143-150
• /
• 1996

This is an explicit-Implicit, finite element analysis for linear as well as nonlinear hygrothermal stress problems. Additional features, such as moisture diffusion equation, crack element and virtual crack extension(VCE ) method for evaluating J-integral are implemented in this program. The Linear Elastic Fracture Mechanics(LEFM) Theory is employed to estimate the crack driving force under the transient condition for and existing crack. Pores in materials are assumed to be saturated with moisture in the liquid form at the room temperature, which may vaporize as the temperature increases. The vaporization effects on the crack driving force are also studied. The Ideal gas equation is employed to estimate the thermodynamic pressure due to vaporization at each time step after solving basic nodal values. A set of field equations governing the time dependent response of porous media are derived from balance laws based on the mixture theory Darcy's law Is assumed for the fluid flow through the porous media. Perzyna's viscoplastic model incorporating the Von-Mises yield criterion are implemented. The Green-Naghdi stress rate is used for the invariant of stress tensor under superposed rigid body motion. Isotropic elements are used for the spatial discretization and an iterative scheme based on the full newton-Raphson method is used for solving the nonlinear governing equations.

• Journal of the Korean Society for Nondestructive Testing
• /
• v.28 no.2
• /
• pp.112-118
• /
• 2008

In the present study, the dependence of ultrasonic phase velocity on porosity and frequency in cancellous bone was predicted using the Biot model and the modified Biot-Attenborough (MBA) model for propagation in fluid-saturated porous media. It was also compared with previously published measurements in human and bovine cancellous bones in vitro. It was shown that the phase velocity in cancellous bone decreased with increasing porosity and frequency The dependence of phase velocity on propagation angle in cancellous bone as predicted using the Schoenberg model together with the Biot model and tile MBA model which were modified to include the effect of angle. The theoretical models used in the present study advance our understanding of the interaction between ultrasound and cancellous bone and can be expected to be usefully employed for the diagnosis of osteoporosis.

• Journal of Soil and Groundwater Environment
• /
• v.26 no.6
• /
• pp.1-17
• /
• 2021

Self-potential (SP) surveys measure naturally occurring differences in electrical potential in the absence of artificial sources and have been applied to various fields since the first application in mineral explorations. Among various causes of SP occurrences, streaming potential is generated by the flow of groundwater, and makes SP surveys suitable for the exploration of groundwater table fluctuation, fractures, sinkholes and landslide occurrences. Recently, there has been many studies that applied SP surveys to monitor water leakage through dikes and embankments. In this review paper, we first review the characteristics and theoretical backgrounds of streaming potential in saturated or unsaturated porous media to introduce it in the embankment among various application field. After the review of the background theory, we review the past cases of field SP surveys on dikes and embankments and also the characteristics of field streaming potential data in the surveys. Further, by analyzing past studies of qualitative as well as quantitative interpretation of SP survey data, we show the possibility of quantitative interpretation of streaming potential data obtained on dikes and embankments. Consequently, it is hope that this review paper helps researches on SP surveys on dikes and embankments, and provides basis for interpretation methods of the SP data to identify leaked area and further leakage rate (or permeability).

• Proceedings of the Korea Water Resources Association Conference
• /
• 2023.05a
• /
• pp.62-62
• /
• 2023

강우가 지표면 아래로 침투할 때 초기에는 투수층이 불포화 상태이어서 부압이 작용하면서 침투할 것이다. Richards 식(Richards, 1931)을 써서 불포화 투수층의 침투를 모의할 수 있다. 강우가 지속되는 동안 하상 아래 어느 구간은 포화 상태가 되어 Richards 식을 더 이상 사용할 수 없다. 하지만 현재까지의 연구는 Richards 식을 사용하여 침투를 모의하는 오류를 범하고 있다. 강우에 의한 침투를 예측할 때 지표면에서의 침투율 q_b 가 필요한 데 현존하는 연구에서는 Horton 식(Horton, 1941)을 사용하여 초기 침투율 f_o 와 장시간 후 침투율 f_c 와 시간에 따라 지수함수로 감소하는 계수 k 의 3가지 계수값을 실험이나 현장 관측값에서 찾아서 쓰고 있다. 그런데, 이 계수값은 강우강도 r_i 가 클수록 침투율 q 가 커지는 물리 현상을 반영하지 못하는 한계가 있다. 본 연구에서 먼저 포화 투수층에서의 침투를 모의하는 식을 개발하였다. 지표면 아래에서 불포화 투수층에는 Richards식을 사용하고 포화 투수층에는 개발한 식을 사용하여 침투를 모의하였다. 또한 지표면에서의 침투율 q_b 를 구하는 공식을 개발하였다. 하상에서의 침투율의 최대값은 $q_{bmax}=-{\lambda}{\sqrt{2g(s-b)}}$ 일 것이다. 여기서 λ 는 투수층의 공극율, s 는 유출수면의 위치, b 는 지표면의 위치이다. 지표면에서의 침투율의 최소값 q_bmin 은 지표면 바로 아래 지점에서의 침투율일 것이다. 지표면에서의 침투율 q_b 로 q_bmax 와 q_bmin 사이의 적절한 값을 선택한다. 강우강도를 r_i 라고 하면 지표면 위 유출수의 연속방정식은 다음과 같다: $s-b={\int}(r_i-{\mid}q_b{\mid})dt$. 즉, 유출수면의 위치 s 는 강우강도 r_i 가 클수록 또는 지표면에서의 유출율의 크기 |q_b| 가 작을수록 크다. 또한 지표면에서의 침투율 q_b 와 지표면 아래에서의 침투율 q 는 s - b 가 클수록 크다. 따라서, 강우강도 r_i 가 클수록 침투율 q_b, q 가 큰 현상이 잘 반영되었다. 강우-침투-유출 모형실험을 수행하여 강우강도에 따라 침투율과 유출량이 다른 현상을 관측하여 수치실험 결과와 비교·검증하였다.