지진 반사 데이터 마이그레이션의 수학적 원리와 구현

초록

본 논문은 반사 지진학에서 마이그레이션 과정을 수학적으로 정밀히 분석한다. 평면 반사체의 정상 이동(NMO) 식을 기울어진 반사체에 적용하면 이동곡선의 최고점이 타원 궤적을 그임을 보이고, 실제 지하 구조는 비균일하고 3차원적 변화를 갖기에 파동 방정식을 기반으로 한 전파 모델이 필요함을 제시한다. WKB 근사와 정적 위상법을 이용해 깊이‑속도 변화를 포함한 파동장의 진폭·위상을 도출하고, 일정 속도 가정하에 Stolt 마이그레이션과 정적 위상법이 동등함을 증명한다.

상세 분석

이 논문은 반사 지진 데이터의 마이그레이션을 두 단계로 나누어 체계적으로 전개한다. 첫 번째 단계에서는 가장 단순한 2차원 모델—수평면에 평탄한 반사층을 가진 경우—을 고찰하고, 소스‑리시버 간 거리 x와 도착 시간 t 사이의 관계를 피타고라스 정리를 이용해 정상 이동(NMO) 방정식 (t^{2}=t_{0}^{2}+x^{2}/v^{2}) 로 도출한다. 여기서 (t_{0}=2d/v) 는 두 번 왕복한 수직 반사 시간이며, 이 식은 하이퍼볼라 형태의 이동곡선을 만든다.

다음으로, 반사면이 일정 각도 θ 로 기울어졌을 때를 고려한다. 기울어진 면에 대해 소스‑리시버 쌍마다 반사점이 달라지므로, 좌표 변환 (\bar{x}=x-2d\sin\theta,\ \bar{t}{0}=2d\cos\theta/v) 을 도입해 동일한 하이퍼볼라 형태를 유지한다. 이때 하이퍼볼라의 최고점 좌표는 ((-2d\sin\theta,\ \bar{t}{0})) 가 되며, θ 가 변함에 따라 최고점이 그리는 궤적은 ((x_{m}^{2}/(2d)^{2})+(t_{m}^{2}v^{2}/(2d)^{2})=1) 이라는 타원 방정식으로 나타난다. 즉, 기울어진 반사면은 이동곡선의 피크를 타원 궤적으로 이동시킨다.

그러나 실제 지하 구조는 다수의 층, 비선형 기울기, 단층 및 깊이에 따른 속도 v(z) 변화를 포함한다. 따라서 단일 x‑t 관계만으로는 충분하지 않으며, 3차원 파동 전파를 기술하는 파동 방정식
\