• Structural Engineering and Mechanics
• /
• 제15권6호
• /
• pp.735-752
• /
• 2003

The quasi-static and dynamic responses of a linear viscoelastic circular beam on Winkler foundation are studied numerically by using the mixed finite element method in transformed Laplace-Carson space. This element VCR12 has 12 independent variables. The solution is obtained in transformed space and Schapery, Dubner, Durbin and Maximum Degree of Precision (MDOP) transform techniques are employed for numerical inversion. The performance of the method is presented by several quasi-static and dynamic example problems.

• 지구물리와물리탐사
• /
• 제21권3호
• /
• pp.139-149
• /
• 2018

전기비저항 모니터링은 지하 매질의 시간적 변화를 파악하기 위한 효과적인 물리탐사법이다. 전기비저항 모니터링은 최근 정교한 자동측정 장비와 통신기술의 발달에 힘입어 지하의 시공간적 물성변화의 탐지를 위하여 널리 사용되고 있다. 이 연구에서는 지하 매질의 의미있는 시공간적 변화대를 효과적으로 해석할 수 있는 4차원 역산 알고리듬을 개발하였다. 4차원 역산의 문제점인 시간축을 따른 지나친 평활화 제한 문제를 극복하기 위하여 자료오차에 대한 시공간 모델제한을 일정하게 유지하도록 라그랑지 곱수를 자동으로 조정하였다. 또한 강조 모델제한자를 도입하여 보다 선명하게 지하의 시공간 변화대를 영상화하고자 하였다. 시간경과 모델에 대한 수치자료에 대하여 개발된 4차원 역산을 수행하여 개발된 역산 알고리듬의 타당성을 검토하였다.

• Structural Engineering and Mechanics
• /
• 제69권5호
• /
• pp.557-566
• /
• 2019

The reservoir basin bedrock produced significant impact on the long-term service safety of super-high arch dams. It was important for accurately identifying geomechanical parameters and its evolution process of reservoir basin bedrock. The deformation modulus mechanism research methods of reservoir basin bedrock deformation modulus for super-high arch dams was carried out by finite element numerical calculation of the reservoir basin bedrock deformation and in-situ monitoring data analysis. The deformation modulus inversion principle of reservoir basin bedrock in a wide range was studied. The convergence criteria for determining the calculation range of reservoir basin of super-high arch dams was put forward. The implementation method was proposed for different layers and zones of reservoir basin bedrock. A practical engineering of a super-high arch dam was taken as the example.

• 지구물리와물리탐사
• /
• 제12권1호
• /
• pp.99-104
• /
• 2009

본 연구에서는 천부 토양층의 2차원 불균질 S파 단면을 결정하기 위해 천부 굴절법 탐사로부터 얻은 SH파 자료에 대한 파형역산 기법을 제안한다. 2차원 매질에서 SH파의 전파를 모사하기 위해 2.5차원 파동 방정식을 사용한다. 합성탄성파를 계산하기 위해 공간축으로 4차, 시간축으로 2차 근사한 staggered grid 유한차분법을 사용하여 파동 방정식을 푼다. 계산된 파형과 측정 파형의 잔여오차로 정의되는 목적함수를 하이브리드 발견적 탐색기법에 의해 최소화한다. 2차원 지하구조 모형은 각기 다른 심도 경계면을 갖는 블록과 블록 내부의 S파 속도에 의해 매개화한다. 수치실험은 암영층파 불규칙한 경계를 갖는 모델에 대해 백색잡음을 추가한 합성 SH파 자료를 이용하여 수행하였다. 지표 굴절법 자료로부터 암영층을 포함한 구조를 적절한 계산시간 내에 영상화할 수 있었다

• 지구물리와물리탐사
• /
• 제10권4호
• /
• pp.229-240
• /
• 2007

이 연구에서는 지하 매질의 속도 이방성을 규명하고 이를 고려하여 지하 단면을 해석하는 탄성파 이방성 토모그래피(seismic anisotropic tomography) 알고리듬을 개발하였다. 이는 지하 매질을 대칭축이 기울어진 횡등방성(Tilted Transversely Isotropic) 매질로 가정하여 지하의 속도 구조와 더불어 이방성 구조를 함께 파악할 수 있도록 고안되었다. 개발된 알고리듬은 이방성 토모그래피에서 흔히 나타날 수 있는 역산의 안정성 문제를 해결하기위해여 역산 매개 변수를 역속도비, 역속도 그리고 이방성 대칭축 방향으로 설정하였고, 의사 베타 변환(pseudo-beta transform)을 통해 이방성 계수의 범위를 제한하였으며, ACB (Active Constraint Balancing) 법을 이용하여 효율적으로 역산 제한조건을 적용하였다. 특히, 기존의 등방성 토모그래피에서 사용하던 프레넬대(Fresnel volume)를 고려한 역산 기법을 이방성 토모그래피로 확장하여 파선 토모그래피에 비해 넓은 범위의 전파각 영역을 확보함으로써 이방성 토모그래피의 안정성을 높일 수 있었다. 개발된 이방성 토모그래피 알고리듬을 복잡한 수치 모델에 적용하여 알고리듬의 타당성을 검증하였으며 현장 자료에 적용하여 시추 결과 및 지질 정보들과 일치하는 향상된 결과를 얻었다. 개발된 이방성 토모그래피 알고리듬을 통해 지하매질의 속도 정보와 이방성 정보를 함께 고려함으로써 보다 정확한 지하 구조의 해석이 가능할 것으로 기대된다.

• 한국지구물리탐사학회:학술대회논문집
• /
• 한국지구물리탐사학회 2007년도 공동학술대회 논문집
• /
• pp.124-129
• /
• 2007

For seismic imaging and inversion, the inverted image depends on how we define the objective function. ${\ell}^1$-norm is more robust than ${\ell}^2$-norm. However, it is difficult to apply the Newton-type algorithm directly because the partial derivative for ${\ell^1$-norm has a singularity. In our paper, to overcome the difficulties of singularities, Huber function given by hybrid ${\ell}^1/{\ell}^2$-norm is used. We tested the robustness of our new object function with several noisy data set. Numerical results show that the new objective function is more robust to band limited spiky noise than the conventional object function.

• 한국지구물리탐사학회:학술대회논문집
• /
• 한국지구물리탐사학회 2003년도 Proceedings of the international symposium on the fusion technology
• /
• pp.402-409
• /
• 2003

Inversion of traveltime requires an efficient algorithm for computing the traveltime as well as its $Frech\hat{e}t$ derivative. We compute the traveltime of the head waves using the damped wave solution in the Laplace domain and then present a new algorithm for calculating the $Frech\hat{e}t$ derivative of the head wave traveltimes by exploiting the numerical structure of the finite element method, the modem sparse matrix technology, and SWEET algorithm developed recently. Then, we use a properly regularized steepest descent method to invert the traveltime of the Marmousi-2 model. Through our numerical tests, we will demonstrate that the refraction tomography with large aperture data can be used to construct the initial velocity model for the prestack depth migration.

• 지구물리와물리탐사
• /
• 제12권3호
• /
• pp.233-238
• /
• 2009

시추공 주시 토모그래피의 주요 요소 중 하나인 초동주시의 측정오차 정보를 사용한 주시 토모그래피 법을 제안하고, 이 역산 방법을 실제 시추공 탐사 자료에 적용하는 실험을 수행 하였다. 제안된 역산방법은 2단계 제어 처리 공정이다. 1) 성취하고자하는 속도 영상의 분해능에 기초하여 측정 주시의 허용 오차한계를 설정하고 측정오차가 오차한계 이상인 기록은 역산에서 제외한다. 2) 반복역산 과정에서 계산된 주시 차 중에서 해당기록의 측정오차 보다 작은 주시 차를 0(영) 처리함으로써 무의미한 수치가 누적되어 속도 개선에 기여함을 방지한다. 모든 기록을 사용한 속도영상과 제어역산으로 도출한 속도영상을 비교한 결과에 의하면, 제어 역산으로 도출한 속도영상에서 수치잡음이 상대적으로 적게 표출되었다. 그리고 Fermat 원리에 의하면, 이 영상이 보다 현실가능성이 높은 속도모델임을 확인하였다.

• 자원환경지질
• /
• 제33권1호
• /
• pp.41-50
• /
• 2000

Geo-radar survey which has the advantage of high-resolution and relatively fast survey has been widely used for engineering and environmental problems. Three-dimensional effects have to be considered in the interpretation of geo-radar for high-resolution. However, there exists a trouble on the analysis of the three dimensional effects. To solve this problem an efficient three dimension numerical modeling algorithm is needed. Numerical radar modeling in three dimensional case requires large memory and long calculating time. In this paper, a finite difference method time domain solution to Maxwell's equations for simulating electromagnetic wave propagation in three dimensional media was developed to make economic algorithm which requires smaller memory and shorter calculating time. And in using boundary condition Liao absorption boundary. The numerical result of cross-hole radar survey for tunnel is compared with real data. The two results are well matched. To prove application to three dimensional analysis, the results with variation of tunnel's incident angle to survey cross-section and the result when the tunnel is parallel to the cross-section were examined. This algorithm is useful in various geo-radar survey and can give basic data to develop dat processing and inversion program.

• 대한전기학회:학술대회논문집
• /
• 대한전기학회 1996년도 추계학술대회 논문집 학회본부
• /
• pp.130-132
• /
• 1996

The power flow calculations(PFc) are the most important and powerful tools in power systems engineering. The conventional power now problem is solved generally with numerical methods such as Newton-Raphson. The conventional numerical method generally have some convergency problem, which is sensitive to initial value, and numerical stability problem concerned with matrix inversion. This paper presents a new power flow calculation algorithm based on the genetic algorithm(GA) which can overcome the disadvantages mentioned above. Some case studies with IEEE 6 bus system also presented to show the performance of proposed algorithm.