본 연구에서는 복잡한 지질구조에 대해서도 신속하고 효율적으로 주시를 계산할 수 있는 Straight Ray Technique(SRT)을 이용한 반사주시 토모그래피 역산 알고리듬을 개발하였다. 역산을 위한 초기 속도모델은 지층경계면에 임피던스 변화를 갖는 상속도 모델을 사용하였다. 실제 속도모델의 반사주시와 초기 속도모델의 반사주시 차이를 계산하여 각각의 요소마다 주시의 오차를 줄이는 방법인 가우스-뉴튼 알고리듬을 이용하여 역산온 수행하였다. 자코비안의 요소는 파선이 지나가는 거리함수로 구성되며, 이를 최소자승형태의 근사 헤시안 행렬로 구성하여 역산을 수행하였다. 역산시 해가 수렴할 수 있도록 근사 헤시안 행렬의 대각성분에 일정한 감쇠인자를 더하였다. 역산된 속도모델을 이용하여 Kirchhoff구조보정을 실시한 결과 실제 속도모델구조에 근사한 단면영상을 얹을 수 있었다.
In this study, we developed reflection tomography inversion algorithm using Straight Ray Technique (SRT) which can calculate travel time easily and fast for complex geological structure. The inversion process begins by setting the initial velocity model as a constant velocity model that hat only impedance boundaries. The inversion process searches a layer-interface structure model that is able to explain the given data satisfactorily by inverting to minimize data misfit. For getting optimal solution, we used Gauss-Newton method that needed constructing the approximate Hessian matrix. We also applied the Marquart-Levenberg regularization method to this inversion process to prevent solution diverging. The ability of the method to resolve typical target structures was tested in a synthetic salt dome inversion. Using the inverted velocity model, we obtained the migration image close to that of the true velocity model.
