• 제목/요약/키워드: domain-inversion

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

시간영역에서 가우스뉴튼법을 이용한 탄성파 파형역산 (Time Domain Seismic Waveform Inversion based on Gauss Newton method)

  • 신동훈;박창업
    • 한국지구물리탐사학회:학술대회논문집
    • /
    • 한국지구물리탐사학회 2006년도 공동학술대회 논문집
    • /
    • pp.131-135
    • /
    • 2006
  • 본 논문에서는 가우스 뉴튼법을 이용한 중합전 탄성파 자료의 파형역산에 관한 연구를 수행하였다. 탄성파 파형역산에 가우스 뉴튼법을 적용하는 방법은 80년대에 제시되었으나 최근 들어서야 활발히 연구가 진행되고 있는데 이는 연산 능력과 기억용량의 한계에 기인한 것이다. 이를 극복하기 위해 본 연구에서는, 파동 전파 수치모의와 역산과정에서 각각 다른 크기의 격자간격을 사용하고, 필요한 시간영역의 파동전파 모사와 가상 진원의 근사를 통해 편미분 파형을 계산하였으며, 효과적으로 슈퍼컴퓨터를 활용하기 위해 병렬처리 기법을 사용하였다. 수치모의를 통해, 가우스 뉴튼법을 이용한 파형 역산의 수렴속도가 빠르고 정확한 것을 알 수 있었으며, 이를 통해 본 연구에서 제시한 방법의 실제 탄성파 자료를 이용한 역산에의 적용가능성을 확인하였다.

  • PDF

PROBLEMS IN INVERSE SCATTERING-ILLPOSEDNESS, RESOLUTION, LOCAL MINIMA, AND UNIQUENESSE

  • Ra, Jung-Woong
    • 대한수학회논문집
    • /
    • 제16권3호
    • /
    • pp.445-458
    • /
    • 2001
  • The shape and the distribution of material construction of the scatterer may be obtained from its scattered fields by the iterative inversion in the spectral domain. The illposedness, the resolution, and the uniqueness of the inversion are the key problems in the inversion and inter-related. The illposedness is shown to be caused by the evanescent modes which carries and amplifies exponentially the measurement errors in the back-propagation of the measured scattered fields. By filtering out all the evanescent modes in the cost functional defined as the squared difference between the measured and the calculated spatial spectrum of the scattered fields from the iteratively chosen medium parameters of the scatterer, one may regularize the illposedness of the inversion in the expense of the resolution. There exist many local minima of the cost functional for the inversion of the large and the high-contrast scatterer and the hybrid algorithm combining the genetic algorithm and the Levenberg-Marquardt algorithm is shown to find efficiently its global minimum. The resolution of reconstruction obtained by keeping all the propating modes and filtering out the evanescent modes for the regularization becomes 0.5 wavelength. The super resolution may be obtained by keeping the evanescent modes when the measurement error and instance, respectively, are small and near.

  • PDF

ERROR ESTIMATIES FOR A FREQUENCY-DOMAIN FINITE ELEMENT METHOD FOR PARABOLIC PROBLEMS WITH A NEUMANN BOUNDARY CONDITION

  • Lee, Jong-Woo
    • 대한수학회보
    • /
    • 제35권2호
    • /
    • pp.345-362
    • /
    • 1998
  • We introduce and anlyze a naturally parallelizable frequency-domain method for parabolic problems with a Neumann boundary condition. After taking the Fourier transformation of given equations in the space-time domain into the space-frequency domain, we solve an indefinite, complex elliptic problem for each frequency. Fourier inversion will then recover the solution of the original problem in the space-time domain. Existence and uniqueness of a solution of the transformed problem corresponding to each frequency is established. Fourier invertibility of the solution in the frequency-domain is also examined. Error estimates for a finite element approximation to solutions fo transformed problems and full error estimates for solving the given problem using a discrete Fourier inverse transform are given.

  • PDF

노후된 콘크리트 구조물의 안전도 평가를 위한 초음파기법의 주파수 및 시간영역 해석에 관한 연구 (A Study on Frequency and Time Domain Interpretation for Safety Evaluation of old Concrete Structure)

  • 서백수;손권익
    • 터널과지하공간
    • /
    • 제15권5호
    • /
    • pp.352-358
    • /
    • 2005
  • 콘크리트 비파괴 검사를 위하여 배면공동모형과 교량공동모형에 대하여 시간영역 탐사와 주파수영역 탐사를 실시하였다. 시간영역 탐사는 초동주시 역산법을 이용하여 토모그래피를 작성하여 공동의 여부를 해석하였다. 주파수영역 탐사는 시간영역 기록을 푸리에 변환에 의한 주파수영역에서의 최대 주파수를 분석하여 해석하였다.

$Ti:LiNbO_3$ 도파로 제작을 위한 열처리 과정 동안 강유전 도메인 특성에 미치는 영향 (Ferroelectric domain inversion in $LiNbO_3$ crystal plate during heat treatment for Ti in-diffusion)

  • 양우석;이형만;권순우;김우경;이한영;윤대호
    • 한국결정성장학회지
    • /
    • 제15권3호
    • /
    • pp.124-127
    • /
    • 2005
  • [ $Ti:LiNbO_3$ ], 광도파로 제작을 위해 큐리온도$(T_c)$ 아래에서 백금박스 내에서 알곤 과 산소 분위기 내에서 열처리 과정 동안 기판 표면의 강유전 도메인 특성 변화를 관찰하였다. 열처리 된 $LiNbO_3$ 기판의 +Z면의 경우 전체적으로 약 $1.6{\mu}m$ 두께로 도메인 반전이 이루어 졌으며, 표면에서 etch hillock이 관찰되었다. $LiNbO_3$ 결정 표면의 Li 이온이 외부로 확산 되는 영향을 감소시킬 수 있는 환경에 있는 기판 면에서 하나의 도메인이 관찰되었으며, 이때 결정 표면에서의 식각특성, 결정성 및 양이온 분포변화에 관하여 X-선 회절, AFM 및 SIMS를 이용하여 분석하였다.

공간영역에서 모멘트방법과 개선된 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

Phase inversion of seismic data

  • Kim, Won-Sik;Shin, Chang-Soo;Park, Kun-Pil
    • 한국지구물리탐사학회:학술대회논문집
    • /
    • 한국지구물리탐사학회 2003년도 Proceedings of the international symposium on the fusion technology
    • /
    • pp.459-463
    • /
    • 2003
  • Waveform inversion requires extracting a reliable low frequency content of seismic data for estimating of the low wave number velocity model. The low frequency content of the seismic data is usually discarded or neglected because of the band-limited response of the source and the receivers. In this study, however small the spectral of the low frequency seismic data is, we assume that it is possible to extract a reliable phase information of the low frequency from the seismic data and use it in waveform inversion. To this end, we exploit the frequency domain finite element modeling and source-receiver reciprocity to calculate the $Frech\`{e}t$ derivative of the phase of the seismic data with respect to the earth model parameter such as velocity, and then apply a damped least squares method to invert the phase of the seismic data. Through numerical example, we will attempt to demonstrate the feasibility of our method in estimating the correct velocity model for prestack depth migration.

  • PDF

Loop-loop EM inversion and its applicability to subsurface exploration

  • Sasaki, Yutaka
    • 한국지구물리탐사학회:학술대회논문집
    • /
    • 한국지구물리탐사학회 2006년도 공동학술대회 논문집
    • /
    • pp.3-6
    • /
    • 2006
  • There are three types of frequency-domain loop-loop EM induction method, depending on the loop separation and their location relative to the ground surface: horizontal-loop EM (HLEM), fixed small-loop EM, and helicopter-borne EM (HEM) methods. Multidimensional inversion provides tomographic images of the subsurface resistivity structure and thus enhances the interpretational accuracy of loop-loop EM data. HLEM method is shown to be effective for exploring groundwater resources in weathered and fractured crystalline basement terrains in semi-arid regions. Also, HEM method is useful for locating weak zones in landslide areas. The applicability of inversion to small-loop EM data depends solely on the S/N ratio. The quadrature response of small-loop EM data can only give the equivalent conductivity of a homogenous half-space model, and thus the in-phase component is essential in inverting EM data. However, the in-phase response is much lower and decreases more rapidly with decreasing frequency than the quadrature response. Further work is needed to obtain conductivity-depth images from small-loop EM data.

  • 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