전이방향성 이방성 매질의 여행시간 계산을 위한 효율적인 티아이 이코날 솔버

초록

본 논문은 전이방향성(transversely isotropic, TTI) 매질의 이코날 방정식을 η 파라미터에 대한 1차 섭동 전개로 근사하여, 복잡한 4차 방정식 해석 대신 선형·이차 PDE만으로 정확한 여행시간을 빠르게 계산하는 방법을 제시한다.

상세 분석

이 논문은 전이방향성 이방성(TTI) 매질에서 발생하는 이코날 방정식의 4차 비선형성을 해소하기 위해, η(비타원성) 파라미터를 기준점 η=0 주변에서 섭동 전개한다는 핵심 아이디어를 도입한다. 기존 방법은 각 격점마다 quartic 방정식을 풀어야 하며, 올바른 P‑wave 해를 선택하는 과정에서 수치적 불안정성이 발생한다. 저자들은 먼저 2차원 TTI 이코날 방정식을 전통적인 1차 차분 스킴으로 이산화하고, τ(i,j)≈τ⁽⁰⁾+τ⁽¹⁾η+τ⁽²⁾η² 형태의 다항식 근사식을 가정한다. 이때 τ⁽⁰⁾는 η=0인 경우, 즉 TEA(equivalent elliptical anisotropy) 방정식의 해이며, 이는 기존 fast‑sweeping 방법으로 효율적으로 구할 수 있다. 이후 η에 대한 1차·2차 항을 차례로 대입해 방정식을 η의 차수별로 정리하면, 각각 f₀(τ⁽⁰⁾)=1, f₁(τ⁽⁰⁾,τ⁽¹⁾)=0, f₂(τ⁽⁰⁾,τ⁽¹⁾,τ⁽²⁾)=0 형태의 연립식이 얻어진다. 저자들은 이 연립식을 풀어 τ⁽¹⁾와 τ⁽²⁾에 대한 명시적 식을 도출했으며, 이 식들은 모두 기존 격점값과 파라미터(v₀, vₙₘₒ, θ)만을 이용한 간단한 연산으로 구현된다. 특히, η가 공간적으로 급격히 변해도 섭동 전개가 이산화 단계에서 적용되므로 매끄러운 η 필드에 대한 제한이 사라진다. 수렴 속도를 높이기 위해 Shanks 변환을 적용해 τ⁽⁰⁾,τ⁽¹⁾,τ⁽²⁾를 조합한 가속된 근사값을 제시한다. 실험에서는 균일 TTI 모델, VTI Marmousi 모델, BP TTI 모델을 대상으로 정확도와 계산 시간을 비교했으며, 기존 4차 방정식 기반 솔버 대비 오차는 12% 이하로 감소하면서 실행 시간은 3050% 정도 절감되었다. 따라서 제안된 방법은 높은 정확도와 낮은 연산 비용을 동시에 만족시키는 실용적인 TTI 이코날 솔버로 평가된다.