강비선형 전파 반전용 직접 엔벨로프 프레셰 도함수 기반 다중 스케일 지진 엔벨로프 반전

초록

전통적인 엔벨로프 반전은 파동 프레셰 도함수에 의존해 약산란 가정에 제한된다. 저자는 에너지 산란 이론을 이용해 체인 룰 없이 직접적인 엔벨로프 프레셰 도함수를 도출하고, 이를 다중 스케일 윈도우‑평균 엔벨로프(WAE)와 결합한 새로운 반전 프레임워크(MS‑DEI)를 제시한다. 강한 대비를 가진 염전·암석 구조에 대해 전통적 FWI와 공동 최적화함으로써 저주파 정보를 효과적으로 복원한다. 1‑D 및 2‑D SEG/EAGE 염전 모델 실험에서 성공을 입증한다.

상세 분석

본 논문은 기존 엔벨로프 반전이 파동 프레셰 도함수(웨이브폼 프레셰 도함수)를 이용해 비선형 미스핏 함수를 선형화하고, 그 과정에서 Born 근사(약산란) 전제에 의존한다는 근본적인 한계를 지적한다. 약산란 가정은 파동이 배경 매질에 대해 작은 변동만을 가질 때만 유효하므로, 염전·암석과 같이 강한 임피던스 대비를 가진 구조에 대해서는 파동 전파가 비선형적으로 왜곡되어 선형화가 무의미해진다. 저자는 이러한 문제를 회피하기 위해 엔벨로프 자체를 에너지 산란 관점에서 재정의한다. 에너지(진폭 제곱)는 위상 정보를 무시하고 진폭만을 다루므로, 파동의 복합 간섭을 평균화한 형태로 강경계 반사를 선형적으로 기술할 수 있다. 이를 기반으로 ‘직접 엔벨로프 프레셰 도함수(Direct Envelope Fréchet Derivative, DEFD)’를 유도했으며, 핵심은 체인 룰을 사용하지 않고 에너지 가상 소스 연산자(Virtual Source Operator, VSO)를 직접 구성한다는 점이다. 약산란(볼륨 산란)에서는 VSO가 배경 그린 함수와 속도 교란의 제곱에 비례하는 형태이지만, 강경계 산란(반사)에서는 경계 조건을 만족하는 Kirchhoff 적분을 통해 얻은 VSO가 주파수 독립적인 반사 계수를 포함한다. 이렇게 얻은 DEFD는 기존 파동 프레셰 도함수에 비해 6차 역속도 의존성을 가지며, 강한 대비 구조에서도 민감도 커널이 정확히 정의된다.

또한 저자는 엔벨로프 데이터를 다중 스케일 윈도우‑평균 엔벨로프(WAE) 형태로 변환한다. WAE는 시간 창을 조절함으로써 저주파 성분을 강조하고, 장파장 속도 구조를 효과적으로 추출한다. WAE와 전통적인 파동 데이터의 잔차를 동시에 최소화하는 공동 목적함수를 설계해, 저주파 정보는 WAE가, 고주파 세부 구조는 FWI가 담당하도록 하였다. 최적화는 L‑BFGS와 같은 2차 미분 기반 알고리즘으로 수행되며, 각 반복에서 DEFD와 파동 프레셰 도함수를 각각 역전파해 그라디언트를 계산한다.

수치 실험에서는 1‑D 두꺼운 염전 층 모델과 2‑D SEG/EAGE 염전 모델을 대상으로 검증하였다. 초기 모델을 단순 1‑D 선형 속도 프로파일로 설정했음에도 불구하고, 전통적 FWI는 강한 사이클‑스키핑과 지역 최소점에 빠져 복원에 실패한다. 반면 MS‑DEI와 FWI를 공동으로 적용한 경우, WAE가 제공하는 저주파 정보가 초기 모델을 빠르게 교정하고, 이후 파동 데이터가 고주파 세부 구조를 정밀하게 복원한다. 결과적으로 염전 경계와 내부 속도 대비가 정확히 재구성되었으며, 전통적 방법 대비 잔차 감소와 모델 오류가 현저히 낮았다.

본 연구는 엔벨로프 반전의 비선형성을 근본적으로 해소하고, 강비선형(Strong‑Nonlinear) 전파 문제에 적용 가능한 새로운 민감도 연산자를 제공함으로써, 저주파 정보가 부족한 현장 데이터에서도 안정적인 대규모 구조 반전이 가능함을 입증한다.