• 제목/요약/키워드: model inversion

검색결과 506건 처리시간 0.025초

주기적 송신원 추출과 참조 송신원 부분집합을 이용한 완전 파형 역산 (Full Waveform Inversion using a Cyclic-shot Subsampling and a Reference-shot Subset)

  • 조상훈;하완수
    • 지구물리와물리탐사
    • /
    • 제22권2호
    • /
    • pp.49-55
    • /
    • 2019
  • 본 연구에서는 주기적 송신원 추출 기법을 사용한 완전 파형 역산 시 목적함수의 안정적인 수렴을 위해 참조 송신원 부분집합을 사용하는 방법을 제안하였다. 완전 파형 역산은 반복적인 파동 전파 모델링을 통해 수행되며, 송신원 개수가 증가할수록 계산 시간이 증가하게 된다. 완전 파형 역산의 계산량을 줄이기 위한 기법들 중 하나로, 주기적 송신원 추출 기법을 사용할 수 있지만 이 경우 역산 초기부터 목적함수가 진동하며 수렴하기 때문에 수렴 판별에 문제가 생기게 된다. 이러한 문제를 해결하고자 본 연구에서는 주기적 송신원 추출 기법을 이용해 모델을 갱신하되, 고정된 참조 송신원 부분집합을 이용해 목적함수를 계산하는 방법을 제안하였다. Marmousi 속도 모델을 이용한 완전 파형 역산 예제를 통해 참조 송신원 부분집합을 이용하면 주기적 송신원 추출 기법을 사용하더라도 목적함수가 안정적으로 수렴할 수 있음을 확인하였다.

시간영역 유도분극 자료의 Cole-Cole 역산 (Spectral Inversion of Time-domain Induced Polarization Data)

  • 김연정;조인기
    • 지구물리와물리탐사
    • /
    • 제24권4호
    • /
    • pp.171-179
    • /
    • 2021
  • 시간영역 유도분극 탐사 자료로부터 Cole-Cole 변수를 추정하는 2차원 역산법을 개발하였다. 모든 유도분극 과도 전위 자료를 역산하여 전기비저항, 충전성, 완화 시간 및 주파수 승수 등의 2D Cole-Cole 변수를 추정하였다. 개발된 역산법은 2단계로 구성된다. 우선 음의 겉보기 충전성 문제를 피하기 위하여 측정된 유도분극 반응을 전류 주입 중 겉보기 전기비저항으로 변환하였다. 1단계 역산에서는 시간에 따라 항상 증가하는 전기비저항을 추정하는 4차원 역산을 통하여 각 역산 블록에서의 전기비저항 시계열 모델을 구축하였다. 2단계 역산에서는 4차원 역산에서 얻어진 전기비저항 시계열 자료를 역산하여 Cole-Cole 변수를 추정하였다. 이때 격자 탐색법을 통하여 참값에 근접한 초기 모델을 설정하는 방법을 통하여 신속한 역산이 가능하였다. 마지막으로 수치 자료에 대한 역산 실험을 통해 개발된 알고리즘이 Cole-Cole 지하 모델을 효과적으로 영상화할 수 있음을 확인하였다.

탄성파 지오토모그래피를 이용한 지하터널 탐사 (Tunnel Detetion Using Seismic Geotomography)

  • 서백수;현병구
    • 터널과지하공간
    • /
    • 제3권1호
    • /
    • pp.50-53
    • /
    • 1993
  • Geotomography, which reconstructs underground structures, is very important task in recent geophysical data processing. In this study, the field data acquired by U.S.Army was used for tomographic inversion and the result was compared with the tomographic inversion and the result was compared with the tomogram from theoretical model data.

  • PDF

Geostatistics for Bayesian interpretation of geophysical data

  • Oh Seokhoon;Lee Duk Kee;Yang Junmo;Youn Yong-Hoon
    • 한국지구물리탐사학회:학술대회논문집
    • /
    • 한국지구물리탐사학회 2003년도 Proceedings of the international symposium on the fusion technology
    • /
    • pp.340-343
    • /
    • 2003
  • This study presents a practical procedure for the Bayesian inversion of geophysical data by Markov chain Monte Carlo (MCMC) sampling and geostatistics. We have applied geostatistical techniques for the acquisition of prior model information, and then the MCMC method was adopted to infer the characteristics of the marginal distributions of model parameters. For the Bayesian inversion of dipole-dipole array resistivity data, we have used the indicator kriging and simulation techniques to generate cumulative density functions from Schlumberger array resistivity data and well logging data, and obtained prior information by cokriging and simulations from covariogram models. The indicator approach makes it possible to incorporate non-parametric information into the probabilistic density function. We have also adopted the MCMC approach, based on Gibbs sampling, to examine the characteristics of a posteriori probability density function and the marginal distribution of each parameter. This approach provides an effective way to treat Bayesian inversion of geophysical data and reduce the non-uniqueness by incorporating various prior information.

  • PDF

수신함수를 이용한 관측소 하부의 지진파 속도구조 (Crustal structure beneath broadband seismic station using receiver function)

  • 박윤경;전정수;김성균
    • 한국지진공학회:학술대회논문집
    • /
    • 한국지진공학회 2003년도 춘계 학술발표회논문집
    • /
    • pp.45-49
    • /
    • 2003
  • The velocity structure beneath the CHNB broadband station is determined by receiver function analysis using by from teleseismic P waveforms. The detailed broadband receiver functions are obtained by stacking method for source-equalized vertical, radial and tangential components of teleseismic P waveforms. A time domain inversion uses the stacked radial receiver function to determine vertical P wave velocity structure beneath the station. The crustal velocity structures beneath the stations are estimated using the receiver function inversion method in the case at the crustal model parameterized by many thin, flat-tying, homogeneous layers. The result of crust at model inversion shows the crustal velocity structure beneath the CHNB station varies smoothly with increasing depth, and there are six discontinuity around 2.5km, 6.25km, 12.5km, 22.5km and 27.5km depth, with Moho discontinuity at about 32.5km depth.

  • PDF

수신함수를 이용한 관측소 하부의 지진파 속도구조 (2) (Crustal structure beneath broadband seismic station using receiver function (2))

  • 박윤경;전정수;김성균
    • 한국지진공학회:학술대회논문집
    • /
    • 한국지진공학회 2003년도 추계 학술발표회논문집
    • /
    • pp.3-7
    • /
    • 2003
  • The velocity structure beneath the CHNB broadband station is determined by receiver function analysis using by from teleseismic P waveforms. The detailed broadband receiver functions are obtained by stacking method for source-equalized vertical, radial and tangential components of teleseismic P waveforms. A time domain inversion uses the stacked radial receiver function to determine vertical P wave velocity structure beneath the station. The crustal velocity structures beneath the stations are estimated using the receiver function inversion method in the case at the crustal model parameterized by many thin, flat-lying, homogeneous layers. Events divide into 4 groups. four azimuths corresponding to events in group a(southwest), b(south), c(southeast), d(northeast). The result of crust at model inversion shows the crustal velocity structure beneath the CHNB station varies smoothly with increasing depth. The conard discontinuity lies around 18 km and moho discontinuity lies range from 30 to 34 km.

  • PDF

Trajectory Guidance and Control for a Small UAV

  • Sato, Yoichi;Yamasaki, Takeshi;Takano, Hiroyuki;Baba, Yoriaki
    • International Journal of Aeronautical and Space Sciences
    • /
    • 제7권2호
    • /
    • pp.137-144
    • /
    • 2006
  • The objective of this paper is to present trajectory guidance and control system with a dynamic inversion for a small unmanned aerial vehicle (UAV). The UAV model is expressed by fixed-mass rigid-body six-degree-of-freedom equations of motion, which include the detailed aerodynamic coefficients, the engine model and the actuator models that have lags and limits. A trajectory is generated from the given waypoints using cubic spline functions of a flight distance. The commanded values of an angle of attack, a sideslip angle, a bank angle and a thrust, are calculated from guidance forces to trace the flight trajectory. To adapt various waypoint locations, a proportional navigation is combined with the guidance system. By the decision logic, appropriate guidance law is selected. The flight control system to achieve the commands is designed using a dynamic inversion approach. For a dynamic inversion controller we use the two-timescale assumption that separates the fast dynamics, involving the angular rates of the aircraft, from the slow dynamics, which include angle of attack, sideslip angle, and bank angle. Some numerical simulations are conducted to see the performance of the proposed guidance and control system.

자력검층을 이용한 파일 심도 예측 (Pile Depth Prediction by Magnetic Logging)

  • 김진후
    • 한국지반공학회:학술대회논문집
    • /
    • 한국지반공학회 2000년도 가을 학술발표회 논문집
    • /
    • pp.231-236
    • /
    • 2000
  • In order to predict depth of the pile forward modeling and inversion of magnetic logging data was conducted by using a finite line of dipoles model. The horizontal component as well as the vertical component of magnetic fields can be measured in the borehole, and the magnetic anomalies can be obtained by subtracting the Earth's magnetic field from the measurement. The magnetic anomalies of the pile are considered as vector sum of induced magnetization due to the Earth's magnetic field and remnant magnetization possessed by steel strings in the pile. The magnetic anomalies are used as input data for inversion from which the length, the magnetic moment per unit length, and the dip angle of the pile can be obtained. From the inversion of synthetic noisy data, and the data obtained from the field model test it is found that the driving depth of the pile can be determined as close to the order of measuring interval (5∼10㎝). It is also found that the resultant magnetic anomalies due to an individual steel string in the pile are almost same as those due to a group of steel strings located at the center of the pile. The magnetic logging method also can be used for locating reinforced bars, pipes, and steel casings.

  • PDF

파장별 복사강도를 사용한 온도 역계산의 적용 (Application of Temperature Inversion by Using Spectral Radiation Intensities)

  • 양수석;송태호
    • 대한기계학회논문집B
    • /
    • 제24권4호
    • /
    • pp.533-542
    • /
    • 2000
  • Analytical experiments to determine the line-of-sight temperature distribution is conducted by using spectral radiation intensities. For this study, fourteen narrow bands of $25cm^{-1}$ interval in $CO_2\;4.3{\mu}m$ band ($2,050cm^{-1}$ to $2375cm^{-1}$) are selected. The applied system is a one-dimensional gas slab filled with 100% $CO_2$ gas at 1 atm. Two types of temperature profile are tested; parabolic and boundary layer types. Three kinds of radiation calculation are used in the iteration procedure for the temperature inversion; LBL(Line by Line), SNB(Statistical Narrow Band) and WNB(WSGGM. based Narrow Band) models. The LBL solution shows perfect agreement while some error of temperature prediction is caused by radiation modeling error when using SNB and WNB models. The inversion result shows that the WNB model may be used more accurately in spectral remote sensing techniques than the traditional SNB model.

3-D gravity terrain inversion for high resolution gravity data analysis

  • Lee Heuisoon;Park Gye-Soon;Kwon Byung-Doo;Oh Seok Hoon;Yang Junmo
    • 한국지구물리탐사학회:학술대회논문집
    • /
    • 한국지구물리탐사학회 2003년도 Proceedings of the international symposium on the fusion technology
    • /
    • pp.379-382
    • /
    • 2003
  • In gravity data correction process, mass effect of the upper part of base level is removed with Bouguer density. Usually, Bouguer density is estimated as a mean density in the field area. But, this may causes a serious problem when ore body is in the area. To overcome this problem, we tried to apply a new method mixing up mass corrections and inversion (3DGTI). 3-D Gravity Terrain Inversion (3DGTI) includes information of topography and distribution of Bouguer density. For this method does not remove the mass effect above base level, it is no longer useless to use Bouguer density. Numerical model tests have shown that the 3DGIT successfully retrieves the anomalous subsurface density distribution of both surface and deeper layers. Model tests shows that this method shows better results than those of conventional one, especially when main target is ore body. The inversion result well delineates the three-dimensional shape of the intruded granite body and basement.

  • PDF