[논문 리뷰] Seismic inverse scattering in the `wave-equation' approach

Christiaan C. Stolk, Maarten V. de Hoop|arXiv (Cornell University)|2001. 12. 17.

Seismic Imaging and Inversion Techniques참고 문헌 17인용 수 25

한 줄 요약

이 논문은 파동방정식 접근법을 사용하여 지구물리학적 역산산산란 문제를 위한 미세국소 분석 프레임워크를 개발한다. 주로 매끄러운 배경에서 파동 전파를 모델링하기 위해 이중 제곱근(DSR) 방정식을 사용한다. 각도 변환을 통해 생성되는 공통 이미지 점 군집이 잡음 없이 작동함을 증명하고, DSR 방정식을 기반으로 한 가짜미분 연산자(annihilators)를 구성하여 비수평 광선 조건 하에서 신뢰성 있는 데이터 중복성 활용이 가능함을 보여준다. 이는 엄밀한 수학적 근거를 바탕으로 한 지구물리학적 이미징 기법이다.

ABSTRACT

Seismic data are commonly modeled by a high-frequency single scattering approximation. This amounts to a linearization in the medium coefficient about a smooth background. The discontinuities are contained in the medium perturbation. The wave solutions in the background medium admit a geometrical optics representation. Here we describe the wave propagation in the background medium by a one-way wave equation. Based on this we derive the double-square-root equation, which is a first order pseudodifferential equation, that describes the continuation of seismic data in depth. We consider the modeling operator, its adjoint and reconstruction based on this equation. If the rays in the background that are associated with the reflections due to the perturbation are nowhere horizontal, the singular part of the data is described by the solution to an inhomogeneous double-square-root equation. We derive a microlocal reconstruction equation. The main result is a characterization of the angle transform that generates the common image point gathers, and a proof that this transform contains no artifacts. Finally, pseudodifferential annihilators based on the double-square-root equation are constructed. The double-square-root equation approach is used in seismic data processing.

연구 동기 및 목표

• 파동방정식 접근법을 사용하여 지구물리학적 역산산산란 문제를 위한 엄밀한 수학적 프레임워크를 개발한다.
• 공통 이미지 점 군집에서 사용되는 각도 변환의 특성을 규명하고, 이가 잡음 없음을 증명한다.
• 이중 제곱근(DSR) 방정식을 기반으로 하향 연장(seismic data downward continuation)을 위한 미세국소 재구성 방정식을 유도한다.
• DSR 프레임워크 내에서 데이터 중복성을 활용하는 가짜미분 연산자(annihilators)를 구성한다.
• 정규화 연산자가 가짜미분 연산자임을 보장하는 조건을 설정하여 매질의 변화를 미세국소적으로 재구성할 수 있도록 한다.

제안 방법

• 일방향 파동방정식을 사용하여 배경 매질 내의 파동 전파를 모델링하고, 이를 통해 일阶 가짜미분 방정식으로서의 이중 제곱근(DSR) 방정식을 도출한다.
• 광선이 어디서도 수평이 아닐 경우, 지구물리학적 데이터의 특이 부분이 비균일 DSR 방정식의 해임을 유도한다.
• 미세국소 분석을 통해 모델링 연산자의 캐논리컬 관계를 특성화하고, 관련된 푸리에 적분 연산자의 가역성을 증명한다.
• 주어진 매질의 속도 c₀에 따라 조정되는 수정된 DSR 연산자(eAWE)를 정의하며, 이는 주성분 기호가 1인 가역적인 푸리에 적분 연산자로 작용한다. 이를 통해 미세국소 재구성을 가능하게 한다.
• 특히 차분 연산자 M = r − s를 활용하여 데이터 중복성을 활용하는 연산자의 복합화를 통해 가짜미분 연산자(annihilators)를 구성한다.
• 수정된 연산자에 정규화된 역연산자를 적용하여, 모든 차수에서 데이터를 제거하는 완전한 연산자 W = K fM eK∗ 를 도출한다.

실험 결과

연구 질문

• RQ1 공통 이미지 점 군집에서 사용되는 각도 변환이 파동방정식 접근법 하에서 엄밀히 잡음 없이 특성화될 수 있는가?
• RQ2 매질의 변화가 존재하는 상황에서 이중 제곱근 방정식이 지구물리학적 데이터를 모델링하는 데 유효한 조건은 무엇인가?
• RQ3 DSR 방정식에서 데이터 중복성을 활용하는 가짜미분 연산자(annihilators)를 어떻게 구성할 수 있는가?
• RQ4 광선이 어디서도 수평이 아닐 경우, 재구성 연산자의 미세국소적 구조는 어떻게 되는가?
• RQ5 파동방정식 접근법에서 정규화 연산자가 가짜미분 연산자임을 입증할 수 있는가? 이를 통해 미세국소 역재구성을 가능하게 하는가?

주요 결과

• 공통 이미지 점 군집을 생성하는 각도 변환이 파동방정식 접근법 하에서 잡음 없이 작동함을 증명하여, 기존 키르히호프 기반 방법에서 제기된 우려를 해결한다.
• 이중 제곱근 방정식이 매끄러운 배경 매질 내에서 지구물리학적 데이터의 하향 연장을 모델링하는 일阶 가짜미분 방정식임을 입증한다.
• 미세국소 재구성 연산자 eAWE가 주성분 기호가 1인 가역적인 푸리에 적분 연산자임을 증명하여, 미세국소 동치성 이하에서 편차를 정확히 복원할 수 있음을 보장한다.
• 모든 차수에서 데이터를 제거하는 가짜미분 연산자 W = K fM eK∗ 가 구성되었으며, 이는 배경 매질 c₀에 의존한다.
• 연산자 W[c0]는 데이터 일관성에 대한 정량적 측도를 제공하는 파동방정식 버전의 차분 일관성 기능(differential semblance functional)의 역할을 한다.
• 증명은 모델링 연산자의 캐논리컬 관계의 미세국소 가역성에 기반하며, 반사 광선이 어디서도 수평이 아닐 조건 하에서 성립한다.