[논문 리뷰] Sensitivity Analysis of an Inverse Problem for the Wave Equation with Caustics

Gang Bao, Hai Zhang|arXiv (Cornell University)|2012. 11. 27.

Numerical methods in inverse problems참고 문헌 6인용 수 4

한 줄 요약

이 논문은 경계 측정값에서 파동방정식의 속도장 복원이라는 역문제에 대해 접선 군집이 존재하는 경우조차도 국소 안정성 추정을 수립한다. 지오데식 X선 변환에 행렬값 가중치를 도입하고, 푸리에 적분 연산자 이론을 활용하여, 딜리클레-노이만 맵에 대한 작은 변화가 접선 정규성 조건 하에서 속도장에 대한 작은 변화를 유도한다는 것을 증명한다.

ABSTRACT

The paper investigates the sensitivity of the inverse problem of recovering the velocity field in a bounded domain from the boundary dynamic Dirichlet-to-Neumann map (DDtN) for the wave equation. Three main results are obtained: (1) assuming that two velocity fields are non-trapping and are equal to a constant near the boundary, it is shown that the two induced scattering relations must be identical if their corresponding DDtN maps are sufficiently close; (2) a geodesic X-ray transform operator with matrix-valued weight is introduced by linearizing the operator which associates each velocity field with its induced Hamiltonian flow. A selected set of geodesics whose conormal bundle can cover the cotangent space at an interior point is used to recover the singularity of the X-ray transformed function at the point; a local stability estimate is established for this case. Although fold caustics are allowed along these geodesics, it is required that these caustics contribute to a smoother term in the transform than the point itself. The existence of such a set of geodesics is guaranteed under some natural assumptions in dimension greater than or equal to three by the classification result on caustics and regularity theory of Fourier Integral Operators. The interior point with the above required set of geodesics is called "fold-regular"; (3) assuming that a background velocity field with every interior point fold-regular is fixed and another velocity field is sufficiently close to it and satisfies a certain orthogonality condition, it is shown that if the two corresponding DDtN maps are sufficiently close then they must be equal.

연구 동기 및 목표

• 기하학적 구조가 단순하지 않은 경우, 특히 고리성이 존재할 때 역파동 전파 문제의 민감도를 다루는 것.
• 고리성과 접선 특이성을 배제하는 표준 '단순한 계량' 가정을 초월하여 안정성 결과를 확장하는 것.
• 새로운 행렬값 가중치를 가진 X선 변환을 도입하여 접선 고리성이 존재하는 경우의 국소 안정성 추정을 위한 프레임워크를 개발하는 것.
• 해밀턴 흐름에 특이성이 유도될 때 속도장 복원 과정의 정규성 조건을 규명하는 것.
• 경계 딜리클레-노이만 맵에 대한 작은 변화가 복원된 속도장에 대한 작은 변화를 유도하는 조건을 설정하는 것.

제안 방법

• 속도장에서 유도된 해밀턴 흐름으로의 사상 연산자를 선형화하여, 행렬값 가중치를 가진 지오데식 X선 변환을 도출한다.
• 접선 정규성 내부 점의 개념을 도입하여, 제어된 접선 고리성이 있는 지오데식이 코탄젠트 공간을 덮을 수 있음을 보장한다.
• 가우시안 비임의 모드를 사용하여 미세지역 파라미트릭스를 구성하고, 진동적 적분 추정을 통해 그 상호작용을 분석한다.
• 푸리에 적분 연산자 이론을 적용하여 X선 변환의 정규성 분석을 수행하고, 접선 정규성 조건 하에서 안정성 추정을 유도한다.
• 점근 전개와 오차 한계를 사용하여 가우시안 비임 성분과 그 도함수의 $ L^2 $-노름 차이를 추정한다.
• 주파수의 역수인 $ \rho = \frac{1}{\text{frequency}} $ 에서의 $ \frac{1}{\rho} $-유형 감쇠를 통해 편향된 비와 원래 비의 차이를 유계로 제한함으로써 X선 변환의 $ L^2 $-안정성을 확립한다.

실험 결과

연구 질문

• RQ1 속도장이 접선 고리성을 유도할 경우, 역파동방정식 문제에 대한 안정성 추정을 수립할 수 있는가?
• RQ2 속도장에 대한 기하학적 조건은 무엇이어야 하며, 딜리클레-노이만 맵이 속도장을 소수의 변화까지 유일하게 결정할 수 있는가?
• RQ3 접선 고리성이 존재할 경우, 파동 전파의 역문제에서 정규성과 안정성에 어떤 영향을 미치는가?
• RQ4 행렬값 가중치를 가진 지오데식 X선 변환은 고리성이 존재하는 경우에도 속도장의 특이성을 복원할 수 있는가?
• RQ5 두 속도장이 그에 해당하는 딜리클레-노이만 맵이 충분히 가까운 경우, 어떤 조건이 두 속도장이 동일함을 보장하는가?

주요 결과

• 두 비포착 속도장이 경계 근처에서 일정하고, 그에 해당하는 딜리클레-노이만 맵이 충분히 가까운 경우, 유도된 산산각 관계는 반드시 동일해야 한다.
• 행렬값 가중치를 가진 지오데식 X선 변환을 도입하였으며, 접선 정규점에서 속도장에 대한 국소 안정성 추정을 수립하였다.
• 지오데식 상에서 접선 고리성이 허용되며, 이는 변환에서 점 자체보다 더 매끄러운 항을 기여하는 조건을 만족할 경우에 한한다.
• 자연적인 가정 하에 고리성의 분류와 푸리에 적분 연산자의 정규성 이론을 통해, 차원 $ d \geq 3 $ 에서 접선 정규성 조건을 만족하는 지오데식의 집합이 존재함이 보장된다.
• 모든 내부 점이 접선 정규인 고정 배경 속도장이 존재할 때, 다른 속도장이 충분히 가까우며 수직 조건을 만족한다면, 딜리클레-노이만 맵이 동일한 경우 속도장도 동일하다.
• 편향된 비와 원래 비의 성분 간 $ L^2 $-노름 차이는 $ \rho = \frac{1}{\text{frequency}} $ 에서 $ \frac{1}{\rho} $ 의 감쇠를 보이며, 이는 파라미트릭스 구성의 안정성을 확인한다.