• Economic and Environmental Geology
• /
• v.44 no.2
• /
• pp.161-170
• /
• 2011

The multiphase flow simulator, CHEMPS, was developed based on the fractional flow approach reported in the petroleum engineering literature considering fully three phase flow in physically and chemically heterogeneous media. It is a extension of MPS developed by Suk and Yeh (2008) to include the effect of wettability on the migration of NAPL. The fractional flow approach employs water, total liquid saturation and total pressure as the primary variables. Most existing models are limited to two-phase flow and specific boundary conditions when considering physically heterogeneous media. In addition, these models focused mainly on the water-wet media. However, in a real system, variations in wettability between water-wet and oil-wet media often occur. Furthermore, the wetting of porous media by oil can be heterogeneous, or fractional, rather than uniform due to the heterogeneous nature of the subsurface media and the factors that affect the wettability. Therefore, in this study, the chemically heterogeneous media considering fractional wettability as well as physically heterogeneous media were simulated using CHEMPS. In addition, the general boundary conditions were considered to be a combination of two types of boundaries of individual phases, flux-type and Dirichlet type boundaries.

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

We have developed a random heterogeneous velocity model with bimodal distribution in methane hydrate-bearing Bones. The P-wave well-log data have a von Karman type autocorrelation function and non-Gaussian distribution. The velocity histogram has two peaks separated by several hundred metres per second. A random heterogeneous medium with bimodal distribution is generated by mapping of a medium with a Gaussian probability distribution, yielded by the normal spectral-based generation method. By using an ellipsoidal autocorrelation function, the random medium also incorporates anisotropy of autocorrelation lengths. A simulated P-wave velocity log reproduces well the features of the field data. This model is applied to two simulations of elastic wane propagation. Synthetic reflection sections with source signals in two different frequency bands imply that the velocity fluctuation of the random model with bimodal distribution causes the frequency dependence of the Bottom Simulating Reflector (BSR) by affecting wave field scattering. A synthetic cross-well section suggests that the strong attenuation observed in field data might be caused by the extrinsic attenuation in scattering. We conclude that random heterogeneity with bimodal distribution is a key issue in modelling hydrate-bearing Bones, and that it can explain the frequency dependence and scattering observed in seismic sections in such areas.

• Journal of the Korean Geosynthetics Society
• /
• v.18 no.3
• /
• pp.69-77
• /
• 2019

This paper presents a new analytical model for estimating the free vibration of tapered friction piles. The governing differential equation for the free vibration of statically axially-loaded piles embedded in non-homogeneous soil is derived. The equation is numerically integrated by the Runge-Kutta method, and then the eigenvalue of natural frequency is determined by the Regula-Falsi scheme. For a cylindrical non-tapered pile, the computed natural frequencies compare well with the available data from literature. Numerical examples are presented to investigate the effects of the tapering, the skin friction resistance, the end condition of the pile, the vertical compressive loading, and the soil non-homogeneity on the natural frequency and mode shape of tapered friction piles.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.40 no.4
• /
• pp.373-379
• /
• 2016

It is known that it is expedient to represent the distribution of the properties of a bone with complex heterogeneity as B-spline volume functions. For B-spline-based representation, the pixel values of CT images are interpolated by B-spline volume functions. However, the CT images of a bone are three-dimensional and very large, and hence a large amount of memory and long computation time for the interpolation are required. In this study, a method for resolving these problems is proposed. In the proposed method, the B-spline volume interpolation problem is simplified by using the uniformity of pixel spacing of the image and the properties of B-spline basis functions. This results in a reduction in computation time and the amount of memory used. The proposed method was implemented and it was verified that the computation time and the amount of memory used were reduced.

• The Korean Journal of Petroleum Geology
• /
• v.1 no.1 s.1
• /
• pp.47-52
• /
• 1993

A method and results of computations are presented for the 2-D seismic migration process in the frequency-wavenumber domain for the laterally and vertically inhomogeneous medium. In order to take the intrinsic attenuation effect into account in the migration process the complex-valued wave velocity is used in the wavefield extrapolation operator, improving the generalized frequency-wavenumber migration technique. The imaginary part of the complex-valued wave velocity includes the seismic quality factor Q value. In derivation of the solution of the wave equation for the medium of inhomogeneous wave velocity and anelasticity, the inhomogeneous medium is mathematically converted to an equivalent system which consists of a homogeneous medium of averaged slowness and an inhomogeneous distribution of hypothetical wave source. The strength of the hypothetical wave source depends on the deviation of squared slowness from the averaged value of the medium. Results of numerical computation using the technique show more distinct geologic images than those using the convensional generalized frequency-wavenumber migration. Especially, the obscured images due to the wave attenuation by anelasticity are restored to show sharp boundaries of structures. The method will be useful in the imaging of the reflection data obtained in the regions of possible petroleum or natural gas reservoir and of fractured zone.

• Journal of the Korean GEO-environmental Society
• /
• v.23 no.8
• /
• pp.29-36
• /
• 2022

The importance of subsurface information is becoming crucial in urban area due to increase of underground construction. The position of underground facilities should be identified precisely before excavation work. Geophyiscal exporation method such as ground penetration radar (GPR) can be useful to investigate the subsurface facilities. GPR transmits electromagnetic waves to the ground and analyzes the reflected signals to determine the location and depth of subsurface facilities. Unfortunately, the readability of GPR signal is not favorable. To overcome this deficiency and automate the GPR signal processing, deep learning technique has been introduced recently. The accuracy of deep learning model can be improved with abundant training data. The ground is inherently heteorogeneous and the spacially variable ground properties can affact on the GPR signal. However, the effect of ground heterogeneity on the GPR signal has yet to be fully investigated. In this study, ground heterogeneity is simulated based on the fractal theory and GPR simulation is carried out by using gprMax. It is found that as the fractal dimension increases exceed 2.0, the error of fitting parameter reduces significantly. And the range of water content should be less than 0.14 to secure the validity of analysis.

• 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.8 no.1
• /
• pp.7-17
• /
• 2005

The determination of seismic velocities in refractors for near-surface seismic refraction investigations is an ill-posed problem. Small variations in the computed time parameters can result in quite large lateral variations in the derived velocities, which are often artefacts of the inversion algorithms. Such artefacts are usually not recognized or corrected with forward modelling. Therefore, if detailed refractor models are sought with model based inversion, then detailed starting models are required. The usual source of artefacts in seismic velocities is irregular refractors. Under most circumstances, the variable migration of the generalized reciprocal method (GRM) is able to accommodate irregular interfaces and generate detailed starting models of the refractor. However, where the very-near-surface environment of the Earth is also irregular, the efficacy of the GRM is reduced, and weathering corrections can be necessary. Standard methods for correcting for surface irregularities are usually not practical where the very-near-surface irregularities are of limited lateral extent. In such circumstances, the GRM smoothing statics method (SSM) is a simple and robust approach, which can facilitate more-accurate estimates of refractor velocities. The GRM SSM generates a smoothing 'statics' correction by subtracting an average of the time-depths computed with a range of XY values from the time-depths computed with a zero XY value (where the XY value is the separation between the receivers used to compute the time-depth). The time-depths to the deeper target refractors do not vary greatly with varying XY values, and therefore an average is much the same as the optimum value. However, the time-depths for the very-near-surface irregularities migrate laterally with increasing XY values and they are substantially reduced with the averaging process. As a result, the time-depth profile averaged over a range of XY values is effectively corrected for the near-surface irregularities. In addition, the time-depths computed with a Bero XY value are the sum of both the near-surface effects and the time-depths to the target refractor. Therefore, their subtraction generates an approximate 'statics' correction, which in turn, is subtracted from the traveltimes The GRM SSM is essentially a smoothing procedure, rather than a deterministic weathering correction approach, and it is most effective with near-surface irregularities of quite limited lateral extent. Model and case studies demonstrate that the GRM SSM substantially improves the reliability in determining detailed seismic velocities in irregular refractors.

• Proceedings of the Korean Society of Soil and Groundwater Environment Conference
• /
• 2005.04a
• /
• pp.72-75
• /
• 2005

본 연구에서는 야외규모 TCE (trichloroethylene) 질량전이역학을 평가하기 위해 TCE의 대수층 내 유입 및 용해상 거동을 고해상도 수치 모사를 실시하였다. 공간적으로 불균질한 대수층 $10{\times}10{\times}10m$ 도메인 내부의 유입된 TCE 분포를 모사하기 위해 본 연구를 통해 개발된 정성적 침투모델(percolation model)이 이용되었다. 이를 기초로 하여 연계되어진 (coupled) 지하수 유동 및 용해상 거동 장기 (long-term) 모사를 실시하였다. 지엽적으로 일어나는 질량전이는 기존 연구를 통하여 발표되어진 실험실규모 연구에 기초하였다. 지하수가 도메인을 지나 흘러나가는 경계면에서 측정되어진 용해상 (aqueous phase) TCE의 질량선속 (mass flux)을 통해 실질 야외규모 질량전이 상수가 계산되었다. 관찰된 바 야외규모 질량전이 상수는 실험실 연구를 통해 측정된 값에 비하여 휠씬 작은 값을 보였으며 이는 지하수 유속 및 TCE의 불균질한 분포에 기인한다. 야외규모 질량전이 상수는 평균 지하수 유속에 직접 비례하는 것으로 관찰되었고 이는 기존 실험실 연구를 통해 알려진 평균 지하수 유속의 0.7승이라는 결과와 대조되는 것이다. 또한 모사를 통해 관찰된 야외규모 질량전이 상수는 상대 TCE 질량의 고갈상수 승에 비례함을 보였다. 이러한 고갈상수는 TCE가 측방으로 퍼지는 현상이 강한 대수층, 즉 저투수성 층의 발달이 양호한 대수층, 에서는 1보다 작은 값을 갖고 그렇지 않은 대수층, 즉 저투수성 층의 발달이 미약한 대수층, 에서는 대체적으로 1보다 높은 값을 갖는 것으로 관측되었다. 이는 DNAPL의 측방 퍼짐이 강한 대수층에서는 용해로 인한 시간에 따른 오염원 부근에서의 농도 감소가 미약하기 때문인 것이며, 그와 반대로 DNAPL의 측방 퍼짐이 약한 대수층에서는 시간이 지남에 따라 용해에 의해 지속적으로 오염원 부근에서의 농도가 감소 또는 소멸함으로 인하여 측정되는 용해상 DNAPL의 질량 선속 역시 계속적으로 감소되는 것으로 밝혀졌다.

• Journal of Biomedical Engineering Research
• /
• v.43 no.6
• /
• pp.369-374
• /
• 2022

Purpose: Focused ultrasound is an emerging technology for treating the brain locally in a noninvasive manner. In this study, we have investigated the influence of skull properties on simulating transcranial pressure field. Methods: A 3D computational model of transcranial focused ultrasound was constructed using female and male CT data to solve for intracranial pressure. For heterogeneous model, the acoustic properties were calculated from CT Hounsfield units based on a porosity. The homogeneous model assigned constant acoustic properties for the single-layered skull. Results: A computational model was validated against empirical data. The homogeneous models were then compared with the heterogeneous model, resulted in 10.87% and 7.19% differences in peak pressure for female and male models respectively. For the focal volume, homogeneous model demonstrated more than 94% overlap compared with the heterogeneous model. Conclusion: Homogeneous model can be constructed using MR images that are commonly used for the segmentation of the skull. We propose the possibility of the homogeneous model for the simulating transcranial pressure field owing to comparable focal volume between homogeneous model and heterogeneous model.