[논문 리뷰] The reconstruction of sound speed in the Marmousi model by the boundary control method

Ivan Ivanov, Μ. I. Belishev외 1명|arXiv (Cornell University)|2016. 09. 24.

Numerical methods in inverse problems참고 문헌 14인용 수 8

한 줄 요약

이 논문은 합성 음향 경계 데이터를 사용하여 복잡한 Marmousi 모델에서 음속을 복원하기 위해 경계 제어 방법(Boundary Control Method, BCM)을 적용하는 것을 보여준다. 이 방법은 매우 비정규적인 음향 광선 경로가 존재하는 영역에서도 상대 오차 10–15% 이내로 평균화된 음속 프로파일을 복원하며, 고해상도 반복 복원 방법을 위한 견고한 초깃값을 제공한다.

ABSTRACT

We present the results on numerical testing of the Boundary Control Method in the sound speed determination for the acoustic equation on semiplane. This method for solving multidimensional inverse problems requires no a priory information about the parameters under reconstruction. The application to the realistic Marmousi model demonstrates that the boundary control method is workable in the case of complicated and irregular field of acoustic rays. By the use of the chosen boundary controls, an `averaged' profile of the sound speed is recovered (the relative error is about $10-15\%$). Such a profile can be further utilized as a starting approximation for high resolution iterative reconstruction methods.

연구 동기 및 목표

• 비정규적인 음향 광선 경로를 가진 현실적인 복잡한 지질 모델에서 경계 제어 방법(Boundary Control Method, BCM)을 이용한 음속 복원의 가능성 테스트.
• 사전에 음속 프로파일에 대한 가정 없이 경계 측정값만을 사용하여 BCM의 성능 평가.
• 고해상도 복원 기법을 위한 후속 반복 복원 기법의 초깃값으로 사용할 수 있는 안정적이고 저해상도의 음속 복원 능력 평가.
• 스칼라 곱 평가의 수치 오차가 더 깊은 깊이에서 복원 정확도에 미치는 영향 조사.
• 지반 탐사에서 흔히 볼 수 있는 매우 비정규적인 파동 전파 환경을 BCM이 어떻게 다룰 수 있는지 탐구.

제안 방법

• BCM은 경계의 일부에 시간 도메인 경계 제어를 적용하여 파동을 자극하고, 시스템의 입력-출력 행동을 코딩한 응답 연산자 $ R^{2T} $ 를 측정한다.
• 이 방법은 파동 해의 완비성과 기하광학 원리를 기반으로 하며, 파동 불연속성과 그들이 광선 沿해 전파되는 방식에 대한 정보를 추출한다.
• 핵심 수학적 도구로는 경계 데이터로부터 스칼라 곱 $ igra c^{-2} u^f u^g igracket $ 를 계산하여 광선 좌표계에서의 파동 영상 복원에 사용한다.
• 복원 과정은 복원된 파동 및 좌표 함수 영상으로부터 카르테시안 좌표계를 반지름 기반(광선 기반) 좌표계로 변환함으로써 진행된다.
• 카르테시안 좌표계와 광선 좌표계 간의 좌표 변환의 자코비안을 계산하여 음속을 복원하며, 광선 관내에서 $ c = rac{1}{ ext{Jacobian}} $ 의 관계를 사용한다.
• 경계 제어에 31개의 공간적 및 32개의 시간적 텐트형 함수 기저를 사용하며, 행렬 요소는 수치 적분을 통해 계산하고 선형 시스템은 LAPACK 루틴을 사용하여 해를 구한다.

실험 결과

연구 질문

• RQ1 경계 제어 방법(Boundary Control Method, BCM)은 음속 구조에 대한 사전 가정 없이 Marmousi와 같은 현실적이고 복잡한 모델에서 음속을 복원할 수 있는가?
• RQ2 경계 제어 수가 제한되고 해상도가 낮을 경우 BCM이 음속의 매끄럽고 평균화된 버전을 얼마나 정확하게 복원하는가?
• RQ3 특히 선형 시스템의 조건수 증가로 인해 깊이가 깊어질수록 수치 적분 정확도가 복원 과정에 미치는 영향은 무엇인가?
• RQ4 매우 비정규적이고 복잡한 음향 광선 패턴이 존재하는 환경에서 BCM이 얼마나 안정적이고 효과적으로 유지되는가?
• RQ5 낮은 해상도의 복원 결과가 후속 고해상도 반복 복원 기법을 위한 신뢰할 수 있는 초깃값으로 사용될 수 있는가?

주요 결과

• BCM은 대부분의 탐사 영역에서 상대 오차 약 10–15% 이내로 Marmousi 모델의 음속을 매끄럽고 평균화된 프로파일로 성공적으로 복원하였다.
• 매우 비정규적이고 복잡한 음향 광선 장이 존재하는 영역에서도 복원이 안정적이고 정확하게 유지되어, 파동면 왜곡에 대한 방법의 견고성을 입증하였다.
• 더 높은 정확도의 내부 적분을 사용한 가짜 복원은 -2km를 초과하는 깊이에서 정확도를 크게 향상시켰으며, 이는 스칼라 곱 평가의 수치 오차가 깊이 해상도를 제한한다는 것을 시사한다.
• 시스템 행렬 $ A^{ heta} $ 의 조건수는 $ heta^4 $ 의 비율로 증가하여 최대 탐사 시간 근처에서 $ 10^5 $ 에서 $ 10^6 $ 수준에 도달하며, 이는 불안정성 문제를 야기하고 깊이 해상도를 제한한다.
• 이 방법은 고해상도 반복 복원 기법을 위한 후속 단계에서 사용할 수 있는 안정적이고 저해상도의 음속 복원을 제공한다.
• 경계에서 슬라이딩 서포트를 사용함으로써 활성 제어 수를 줄였지만 공간 해상도를 유지하였으며, 이는 수치적 안정성과 효율성을 향상시켰다.