• Title/Summary/Keyword: Numerical inversion method

Search Result 148, Processing Time 0.029 seconds

2차원 전기비저항 모니터링 자료의 시간경과 역산 (Time-lapse Inversion of 2D Resistivity Monitoring Data)

  • 김기주;조인기;정재형
    • 지구물리와물리탐사
    • /
    • 제11권4호
    • /
    • pp.326-334
    • /
    • 2008
  • 전기비저항 탐사법은 지하 천부의 전기비저항 영상화에 널리 사용되어 왔으며, 최근 신속한 자동측정이 가능해짐에 따라 모니터링이 용이하게 되었다. 이 연구에서는 전기비저항 모니터링 자료의 해석을 위한 시간경과 역산법을 개발하였다. 개발된 역산법은 기준모델에 비하여 변화가 큰 모델변수에는 약한 제한을, 변화가 미미한 모델변수에는 강한 제한을 가하는 방법이다. 수치시험을 통하여 시간경과 역산 영상은 시간적으로 변화한 영역을 보다 정확하고 뚜렷하게 영상화함을 확인하였다. 또한 저수지 누수문제에 개발된 시간경과 역산법을 적용하여 누수구간을 탐지할 수 있었으며, 이 누수구간은 측정 기간내에 크게 변화하지 않았음을 확인할 수 있었다.

Depth Scaling Strategy Using a Flexible Damping Factor forFrequency-Domain Elastic Full Waveform Inversion

  • Oh, Ju-Won;Kim, Shin-Woong;Min, Dong-Joo;Moon, Seok-Joon;Hwang, Jong-Ha
    • 한국지구과학회지
    • /
    • 제37권5호
    • /
    • pp.277-285
    • /
    • 2016
  • We introduce a depth scaling strategy to improve the accuracy of frequency-domain elastic full waveform inversion (FWI) using the new pseudo-Hessian matrix for seismic data without low-frequency components. The depth scaling strategy is based on the fact that the damping factor in the Levenberg-Marquardt method controls the energy concentration in the gradient. In other words, a large damping factor makes the Levenberg-Marquardt method similar to the steepest-descent method, by which shallow structures are mainly recovered. With a small damping factor, the Levenberg-Marquardt method becomes similar to the Gauss-Newton methods by which we can resolve deep structures as well as shallow structures. In our depth scaling strategy, a large damping factor is used in the early stage and then decreases automatically with the trend of error as the iteration goes on. With the depth scaling strategy, we can gradually move the parameter-searching region from shallow to deep parts. This flexible damping factor plays a role in retarding the model parameter update for shallow parts and mainly inverting deeper parts in the later stage of inversion. By doing so, we can improve deep parts in inversion results. The depth scaling strategy is applied to synthetic data without lowfrequency components for a modified version of the SEG/EAGE overthrust model. Numerical examples show that the flexible damping factor yields better results than the constant damping factor when reliable low-frequency components are missing.

Development of an Inversion Analysis Technique for Downhole Testing and Continuous Seismic CPT

  • Joh, Sung-Ho;Mok, Young-Jin
    • 한국지반공학회지:지반
    • /
    • 제14권3호
    • /
    • pp.95-108
    • /
    • 1998
  • 지반의 S파 및 P파의 깊이에 따른 변화를 원위치에서 측정하기 위하여 다운흘 시험 (downhole testing)과 SCPT (seismic CPT) 등이 널리 사용되어 왔다. 다운홀 시험과 SCPT는 경제성, 운용의 용이성, 발진원의 단순성 등의 측면에서 효율적이기 때문에, 현재 지반조사에서 그 사용빈도가 더욱 증가하고 있는 추세이다. 특히 최근에는 다운흘과 SCPT의 자료 분석을 자동화하기 위한 노력의 일환으로 interval measurements의 기법이 활용되고 있는데, 현재 이에 대한 적절한 역산해석 (inversion analysis) 기법이 없는 형편이다. 따라서, 본 논문에서는 다운홀이나 SCPT의 interval measurements를 분석하기 위한 새로운 역산해석 기법을 제안하였다. 제안된 역산해석 기법의 정모델링(forward modeling)에서는 탄성파의 전파를 Snell의 법칙에 의거하여 굴절.반사되는 현상을 고려하였곡, 역산해석을 위해서는 최대공산법 (maximum likelihood method)을 적용하였다. 그리고, 본 논문에서 제안한 역산해석 기법의 검증을 위하여, 하나의 S파 주상도를 가정하고 이에 대하여 다운흘 시험을 모사하였다. 이론적으로 수행한 다운홀 시험 결과에 대하여 기존의 비 역산해석 방법과 본 논문에서 제안한 역산해 석 기법에 의해서 S파 주상도를 추정하였는데, 그 결과 본 논문에서 제시한 역산기법이 가장 정확한 결과를 도출하였으며, 다운홀 시험과 SCPT을 자동화하는데 효율적으로 적용이 될 수 있음을 입증하였다.

  • PDF

전파형 역산을 이용한 시추공 영상의 분해능 (Resolution Limits of Cross-Well Seismic Imaging Using Full Waveform Inversion)

  • 조창수;이희일;서정희
    • 지구물리와물리탐사
    • /
    • 제5권1호
    • /
    • pp.33-45
    • /
    • 2002
  • 주시 토모그래피의 한계인 분해능을 극복하기 위하여 새로운 영상화 기법이 요구되며 그 중의 하나가 파형 역산이다. 파형 역산은 위상뿐만 아니라 파동의 진폭을 동시에 이용하므로 지하구조를 고해상으로 영상화할 수 있는 기법이다. 그러나 파형역산은 전파와 역전파의 모형반응 계산이 요구되므로 많은 계산 시간이 요구된다. 본 연구에서는 파형역산 기법에서 효율적인 합성 파동장 계산을 위하여 속도-응력법을 이용하였다. 시추공 영상화 기법들의 분해능을 알아 보기 위하여 수치모형에 적용, 비교하여 파형역산과 주시 토모그래피의 분해능을 살펴보았다. 파형역산의 분해능 한계는 Schuster가 유도한 구조보정의 분해능 한계와 유사함을 알 수 있었다. 시추공의 기하학적인 문제로 인한 커버리지의 부족의 문제는 VSP자료를 적용함으로서 해결할 수 있어 수평적인 분해능이 향상되었다. 또한 구현된 알고리듬의 현장적용성을 평가하기 위하여 실제와 유사한 이론모형에 적용해보았으며 이 때 발생하는 비선형성을 줄이기 위해 초기치로 주시역산 토모그램을 적용하여 좋은 결과를 얻었다.

공간영역에서 모멘트방법과 개선된 Newton 알고리즘을 이용한 반복 역산란 방법 (Iterative Inversion Using Moment Method and Improved Newton`s AIgorithmin the Configuration Domain)

  • 박천석;김정혜;라정웅
    • 전자공학회논문지A
    • /
    • 제31A권5호
    • /
    • pp.39-49
    • /
    • 1994
  • In this paper, An inversion technique to reconstruct permittivity profiles of 2-D inhomogeneous dielectric objects by iterativeprocess using the moment-methodand improved newton's algoritham is presented. In order to reduce the noise effect in the scattered fieldon the reconstructed permittivity profiles, the cell size of inversescattering is made be larger than that of forward scattering. Performing numerical calculations of dielectric scatterer it is demonstrated that this inversion is able to reconstruct dielectric objectshaving large size and inhomogeneous characteristics, which is insentive tothe noise effect in the scattered field on the reconstructed result.

  • PDF

탄성파 수직반사자료의 역산 연구 (A Study on Inversion of Seismic Normal Reflection Data)

  • 양동우;양승진;장성형
    • 자원환경지질
    • /
    • 제29권5호
    • /
    • pp.629-637
    • /
    • 1996
  • In this paper a numerical experiment is conducted to determine the low acoustic impedance of a thin oil or gas reservoir from a seismogram by using the generalized linear inversion method. The seismograms used are normal incident synthetic seismograms containing p-wave primary reflections, multiples, and peg-leg multiples on the layers consisting of oil-, gas-, water-filled sandstone incased in shales. In this experiment the acoustic impedance, the location of reservoir boundary, thickness, and source wavelet are assumed initially and revised iteratively by the least-squares-error technique until the difference between the seismogram and calculated one is very small. This experiment shows that the acoustic impedance and thickness, about 10 m thick, can be determined by the inversion.

  • PDF

전기비저항 탐사자료의 3차원 역산 (Three-dimensional Inversion of Resistivity Data)

  • 이명종;김정호;조성준;정승환;송윤호
    • 지구물리와물리탐사
    • /
    • 제2권4호
    • /
    • pp.191-201
    • /
    • 1999
  • 현재까지 전기비저항 탐사자료의 해석은 대부분 지하구조가 1차원 또는 2차원으로 가정하는데 근간을 두고 있다. 그러나, 국내와 같이 지질구조가 복잡하며 지형의 기복이 심한 경우, 이와 같은 가정을 사용은 전기비저항 탐사자료의 해석에 오류를 야기할 가능성이 높다. 본 연구에서는 3차원적으로 복잡하게 변화하는 지하구조 및 지형의 기복에 의한 지형효과를 정확히 계산할 수 있는 유한요소법에 의한 3차원 전기비저항 모델링 알고리듬을 근간으로 하는 역산 알고리듬을 개발하였다. 이 알고리듬은 ACB(Active Constraint Balancing)법을 채용함으로써 역산의 분해능 향상을 기하였다. 간단한 3차원 지하구조 및 지형기복에 대한 수치모형 실험결과 지표의 불규칙한 기복을 포함하여 역산을 수행함으로써 정확한 지하구조 영상을 획득할 수 있음을 보였으며, 국내의 교량건설 예정부지에서 수행된 격자망을 이룬 전기비저항 탐사자료에 대한 3차원 역산 결과 매우 신뢰도 높은 지하구조에 대한 영상을 도출할 수 있었다.

  • PDF

전자파의 과도적 산란 특성에 관한 새로운 수치해석 연구 (On the Study of New Numerical Analysis in the Transient Electromagnetic Wave Scattering)

  • 이강호;이상회;김정기
    • 한국통신학회논문지
    • /
    • 제15권1호
    • /
    • pp.1-8
    • /
    • 1990
  • 본 논문에서는 유전체 원통주의 전자파 과도 응답 산란 특성을 새로운 수치 해석 방법으로 연구하낟. 수치 해석에 사용하는 경계적분 방정식(BIE)의 기본 수식은 가중 잔차법으로 부터 시작된다. 경계 적분 방정식을 확장 경계 조건과 표면 경계 조건을 물체의 표면점 내부와 외부에 적용시켜 두개의 연립방정식을 만든다. 수치해석으로 경계요소법을 사용하고, 이 방법에는 직접법과 간접법 두 가지가 있다. 그래서 연산자를 메트릭스 역변환하는 간접법을 수치해석방법으로 사용한다. 결과값은 다른 연구 결과값과 좋은 일치를 보인다.

  • PDF

출력재정의를 통한 비선형제어 기법의 미사일 오토파일롯 응용 (Application of nonlinear control via output redefinition to missile autopilot)

  • 류진훈;탁민제
    • 제어로봇시스템학회:학술대회논문집
    • /
    • 제어로봇시스템학회 1996년도 한국자동제어학술회의논문집(국내학술편); 포항공과대학교, 포항; 24-26 Oct. 1996
    • /
    • pp.1496-1499
    • /
    • 1996
  • A nonlinear tracking control technique developed for the control of nonlinear systems has been applied to the autopilot design of missile system. The difficulties in the application of inversion based control methods such as input-output feedback linearization and sliding mode control due to nonminimum phase characteristics are discussed. To avoid the stability problem associated with unstable zero dynamics, the input-output feedback linearization is applied with output-redefinition method to normal acceleration control. The output-redefinition method gives an indirect way to apply the nonlinear controls to nonminimum phase plants by redefining the plant output such that the tracking control of the modified output ensures the asymptotic tracking of the original output. The numerical simulation shows satisfactory results both for nominal and for slightly perturbed missile systems adopting the sliding mode control technique. However, the robustness problem in this method is briefly investigated and verified with the simulation.

  • PDF

과냉각을 동반한 동결과정의 수치해석 (Numerical Analysis for Cooling and Freezing Processes with Subcooling)

  • 윤정인;김재돌;김성규
    • 설비공학논문집
    • /
    • 제8권4호
    • /
    • pp.451-462
    • /
    • 1996
  • In this study, which focuses on ice storage, a fundamental study in cooling and solidification was performed, including the interesting phenomena of density inversion, supercooling and dendritic ice. A numerical study was performed for natural convection and ice formation considering existence of subcooling and dendritic ice were analyzed numerically by using finite difference method and boundary fixing method. In the mesh, the solid fraction was introduced with adding as a term to the energy conservation equation. A flow in the dendrite was modelled as a flow in a porous medium, and the momentum conservation equation was modified to incorporate resistance forces involved in flows through porous media. A numerical solution of the time dependencies of dendrite area and dense ice front was successfully obtained, and the numerical results were good agreement with experimental results. Based on this methodology, a discussion was made of phenomena and characteristics of cooling and freezing processes under various conditions.

  • PDF