• Proceedings of the Korea Water Resources Association Conference
• /
• 2017.05a
• /
• pp.538-538
• /
• 2017

우리나라의 현행 수질오염총량제도에서 비점오염원에 대한 배출량은 배출원별로 유량조건에 따라 배출계수를 적용하여 산정하고 있으나, 강우 및 유역특성에 따라 크게 좌우되는 비점오염 배출특성을 반영하는 것에는 어려움이 있다. 한편, 역산모델링은 결과를 바탕으로 원인 값을 추정하는 방식으로 유량 및 수질측정결과, 점오염원 배출량 등 실측결과를 바탕으로 비점오염 배출량을 추정할 수 있다. 본 연구에서는 역산 모델링 기법을 이용하여 경안천 유역에서 배출되는 총인의 비점오염 배출량을 추정하고, 부하지속곡선을 통해 비점배출기준을 산정하고자 한다. 경안천 유역의 수질 및 유량측정자료, 하수처리장 현황, 오염원 자료 등을 수집하여, 점오염원 배출량을 추정하고 역산 모델링을 수행하였으며, 비점배출총량을 산정하였다. 본 연구의 결과는 현행 수질오염 총량제도의 배출계수에 대한 재검토와 향후 유량조건에 따른 차별적인 비점오염원 관리 필요성에 대한 기초자료로 활용될 수 있을 것으로 사료된다.

• Proceedings of the KSEEG Conference
• /
• 2002.04a
• /
• pp.198-200
• /
• 2002

경제적이며 정확한 3차원 전자탐사 모델링을 위해 위해 Habashy et al. (1993)에 의해 제안된 국소 비선형 근사(localized nonlinear approximation)를 이용하여 전자탐사 모델링 알고리듬을 개발하였다. 전자탐사 수치모델링시 많은 계산시간 및 기억용량을 필요로 하는 Green 텐서 적분을 정확하고 빠르게 계산하기 위해, 단일 미소체를 이용한 공간파수 영역에서의 Green 텐서 적분 알고리듬을 제안하였다. 더욱이 Green 텐서의 송수신 방향 및 상반성을 고려하여 각각의 미소체에 의한 전체 미소체에의 Green 텐서 적분을 한 개의 미소체에 의한 전체 미소체에의 Green 텐서 적분 값으로 구하게 하므로 매우 적은 기억용량 만으로 Green 텐서 적분 행렬을 구성할 수 있어, 역산법에 효과적으로 적용할 수 있다. 이 수치 모델링 알고리듬을 기본으로 하여 평활화 제한을 가한 최소자승 역산 알고리듬을 개발하였다. 이 역산 알고리듬을 지표 전자탐사 및 시추공-지표 전자탐사 등에 적용하여 PC에서도 빠르게 3차원 전자탐사 역산이 수행됨을 보였다.

• Geophysics and Geophysical Exploration
• /
• v.10 no.2
• /
• pp.138-146
• /
• 2007

This paper describes 2.5D induced polarization (IP) modeling and inversion algorithms using complex resistivity. The complex resistivity method has merits for acquiring more valuable information about hydraulic parameters and pore fluid than the conventional IP methods. The IP modeling and inversion algorithms are developed by allowing complex arithmetic in existing DC modeling and inversion algorithms. The IP modeling and inversion algorithms use a 2.5D DC finite-element algorithm and a damped least-squares method with smoothness constraints, respectively. The accuracy of the IP modeling algorithm is verified by comparing its responses of two synthetic models with two different approaches: linear filtering for a three-layer model and an integral equation method for a 3D model. Results from these methods are well matched to each other. The inversion algorithm is validated by a synthetic example which has two anomalous bodies, one is more conductive but non-polarizable than the background, and the other is polarizable but has the same resistivity as the background. From the inverted section, we can cleary identify each anomalous body with different locations. Furthermore, in order to verify its efficiency to the real filed example, we apply the inversion algorithm to another three-layer model which includes phase anomaly in the second layer.

• Proceedings of the KSEEG Conference
• /
• 2002.04a
• /
• pp.140-142
• /
• 2002

• Proceedings of the KSEEG Conference
• /
• 2000.04b
• /
• pp.155-157
• /
• 2000

• Geophysics and Geophysical Exploration
• /
• v.5 no.4
• /
• pp.291-298
• /
• 2002

We have developed a two-dimensional (2-D) magnetotelluric (MT) inversion algorithm, which can include topographic effects in inversion. We use the finite element method (FEM) to incorporate topography into forward calculation. Topography is implemented simply by moving nodes of rectangular elements in z-direction according to the elevation of air-earth interface. In the inversion process, we adopt a spatially variable Lagrangian multiplier algorithm in the smoothness-constrained least-squares inversion. The inversion algorithm developed in this study reconstructs subsurface resistivity structure quite well when topography variation exists. Also, it turns out to be effective in both resolution and stability from a model study and field data application.

• Geophysics and Geophysical Exploration
• /
• v.16 no.1
• /
• pp.27-35
• /
• 2013

In this study, we developed the trial product on 3D magnetic inversion tentatively named 'KMag3D'. Also, we briefly introduced its own function and graphic user interface on which especially focused through the development in the form of user manual. KMag3D is consisted of two fundamental frame for the 3D magnetic inversion. First, algebraic reconstruction technique was selected as a 3D inversion algorithm instead of least square method conventionally used in various magnetic inversion. By comparison, it was turned out that algebraic reconstruction algorithm was more effective and economic than that of least squares in aspect of both computation time and memory. Second, for the effective determination of the 3D initial and a-priori information model required in the execution of our algorithm, we proposed the practical technique based on the assemblage of 2D forward modeling and inversion results for individual user-selected 2D profiles. And in succession, initial and a-priori information model were constructed by appropriate interpolation along the strke direction. From this, we concluded that our technique is both suitable and very practical for the application of 3D magentic inversion problem.

• Geophysics and Geophysical Exploration
• /
• v.26 no.4
• /
• pp.268-274
• /
• 2023

Finite source inversion is performed with a Green's function matrix and geodetic coseismic displacement. Conventionally, the Green's function matrix is constructed using the Okada model (Okada, 1985). However, for more realistic earthquake simulations, recent research has widely adopted the physics-based model, which can consider various material properties such as elasticity, viscoelasticity, and elastoplasticity. We used the physics-based software PyLith, which is suitable for earthquake modeling. However, the PyLith does not provide a mesh generator, which makes it difficult to perform finite source inversions that require numerous subfaults and observation points within the model. Therefore, in this study, we developed CPInterface (COMSOL-PyLith Interface) to improve the convenience of finite source inversion by combining the processes of creating a numerical model including sub-faults and observation points, simulating earthquake modeling, and constructing a Green's function matrix. CPInterface combines the grid generator of COMSOL with PyLith to generate the Green's function matrix automatically. CPInterface controls model and fault information with simple parameters. In addition, elastic subsurface anomalies and GPS observations can be placed flexibly in the model. CPInterface is expected to enhance the accessibility of physics-based finite source inversions by automatically generating the Green's function matrix.

• The Journal of Engineering Geology
• /
• v.14 no.1
• /
• pp.21-28
• /
• 2004

Three-dimensional (3-D) resistivity inversion including a topographic effect can be considered theoretically to be the technique of acquiring the most accurate image in the interpretation of resistivity data, because it includes characteristic image that the actual subsurface structure is 3-D. In this study, a finite-element method was used as the numerical method in modeling, and the efficiency of Jacobian calculation has been maximized with sensitivity analysis for the destination block in inversion process. Also, during the iterative inversion, the resolution of inversion can be improved with the method of selecting the optimal value of Lagrange multiplier yielding minimum RMS(root mean square) error in the parabolic equation. In this paper, we present synthetic examples to compare the difference between the case which has the toprographic effect and the other case which has not the effect in the inversion process.

• Geophysics and Geophysical Exploration
• /
• v.23 no.1
• /
• pp.38-49
• /
• 2020

Full-waveform inversion (FWI) is an optimization process of fitting observed and modeled data to reconstruct high-resolution subsurface physical models. In acoustic FWI (AFWI), pressure data acquired using a marine streamer has mainly been used to reconstruct the subsurface P-wave velocity models. With recent advances in marine seismic-acquisition techniques, acquiring multi-component data in marine environments have become increasingly common. Thus, AFWI strategies must be developed to effectively use marine multi-component data. Herein, we proposed an AFWI strategy using horizontal and vertical particle-acceleration data. By analyzing the modeled acoustic data and conducting sensitivity kernel analysis, we first investigated the characteristics of each data component using AFWI. Common-shot gathers show that direct, diving, and reflection waves appearing in the pressure data are separated in each component of the particle-acceleration data. Sensitivity kernel analyses show that the horizontal particle-acceleration wavefields typically contribute to the recovery of the long-wavelength structures in the shallow part of the model, and the vertical particle-acceleration wavefields are generally required to reconstruct long- and short-wavelength structures in the deep parts and over the whole area of a given model. Finally, we present a sequential-inversion strategy for using the particle-acceleration wavefields. We believe that this approach can be used to reconstruct a reasonable P-wave velocity model, even when the pressure data is not available.