A direct imaging method for inverse scattering problem of biharmonic wave with phased and phaseless data

기술적 요약: 역 바이하모닉 산란 문제를 위한 직접 영상화 방법

문제 제기
본 논문은 측정 데이터를 이용하여 $\mathbb{R}^2$ 내의 확장된 장애물 $D$의 형상과 위치를 결정하는 역 바이하모닉 산란 문제를 다룬다. 음향, 탄성, 전자기파와 같은 2 차 미분 연산자와 관련된 더 일반적인 역 산란 문제와 달리, 본 연구는 4 차 편미분 방정식인 바이하모닉 파동 방정식에 초점을 맞춘다. 본 연구는 위상 정보가 포함된 위상 데이터 (phased data) 와 필드의 크기 (modulus) 만이 측정된 위상 없는 데이터 (phaseless data) 모두를 고려한다. 장애물은 판 이론에 기반한 네 가지 경계 조건 중 하나를 따르며, 이는 고정 (clamped), 단순 지지 (simply supported), 롤러 지지 (roller supported), 또는 자유 (free) 조건이다.

방법론
저자들은 역시간 이동 (Reverse Time Migration, RTM) 기법에 기반한 직접 영상화 방법을 제안한다. 반복적 방법과 달리, 이 접근법은 산란체에 대한 사전 기하학적 가정을 요구하지 않으며, 각 반복 단계에서 직접 산란 문제를 수치적으로 풀 필요도 없다.

방법론의 핵심은 영역 $\Omega$ 내의 샘플링 점 $z$에 대해 특정 영상화 함수 ($I_j(z)$) 를 구성하는 것이다. 이러한 함수들은 샘플링 점이 장애물 경계 근처에 있는지 아니면 멀리 있는지에 따라 뚜렷한 특성을 보이도록 설계되었다.

1. 위상 데이터 (IP-1 에서 IP-10): 저자들은 다양한 측정 데이터의 조합을 활용하여 10 개의 영상화 함수를 정의한다. 여기에는 다음이 포함된다:

   • 산란장 ($u^{sc}$) 과 그 법선 도함수 ($\partial_n u^{sc}$).
   • 굽힘 모멘트 ($M u^{sc}$) 와 횡방향 힘 ($N u^{sc}$).
   • 원거리 패턴 ($u^\infty$).
   • 점 소스와 평면파에 의해 생성된 데이터.

   영상화 함수는 바이하모닉 파동 방정식의 기본 해가 아닌 헬름홀츠 방정식 ($\Phi_\kappa$) 의 기본 해를 사용하여 구성된다. 이러한 공식화는 계산 복잡도를 줄이기 위한 것이다. 이러한 함수들은 측정 경계 ($\Gamma_r, \Gamma_s$) 또는 단위 원 ($S$) 에 대한 적분을 포함하며, 측정 데이터와 헬름홀츠 기본 해 및 그 도함수를 결합한다.

2. 위상 없는 데이터 (IP-11): 총장의 크기 $|u|$ 만이 이용 가능한 경우, 특정 영상화 함수 $I_{11}(z)$ 가 구성된다. 이 함수는 입사장에 의해 정규화된 총장의 제곱 크기에서 입사장의 제곱 크기를 뺀 차이에 의존한다.

이론적 분석 및 주요 기여
본 논문은 제안된 영상화 함수에 대한 엄밀한 수학적 분석을 제공한다:

• 분해능 분석: 저자들은 임의의 샘플링 점 $z$에 대해 영상화 함수 $I_j(z)$ 가 주요 항과 잡음/오차 항으로 분해될 수 있음을 증명한다. 주요 항은 헬름홀츠 기본 해의 허수 부분에서 유도된 경계 데이터를 갖는 바이하모닉 파동 방정식을 만족하는 보조 해 $\psi(\cdot, z)$ 의 원거리 패턴의 제곱 크기의 적분에 비례한다.
• 대조 메커니즘: 분석에 따르면, 영상화 함수는 샘플링 점 $z$가 장애물 경계 근처에 있을 때 높은 대조도를 보이며, $z$가 멀어짐에 따라 감소한다. 이러한 거동은 헬름홀츠 방정식 기본 해의 허수 부분 (베셀 함수) 의 특성에 기인하며, 이는 소스와 관측 점이 일치할 때 최대값을 갖는다.
• 공식화의 독창성: 중요한 기여는 바이하모닉 기본 해가 아닌 헬름홀츠 방정식 (또는 음향 평면파) 의 기본 해에만 의존하는 영상화 함수를 유도한 것이다. 저자들은 2 차 미분 연산자의 높은 차수 때문에 음향, 탄성, 전자기파에 대한 이전 RTM 연구들이 바이하모닉 경우에 직접 확장될 수 없다고 지적한다.
• 경계 조건에 대한 강건성: 이론적 프레임워크는 네 가지 유형의 경계 조건 (고정, 단순 지지, 롤러 지지, 자유) 모두에 적용 가능하도록 수립되었으며, 각각에 대한 특정 점근적 거동이 유도되었다.

수치적 결과
본 논문은 이론적 발견을 검증하기 위한 수치 실험을 제시한다:

• 단일 장애물 재구성: 알고리즘은 11 개의 영상화 함수 (10 개의 위상 데이터, 1 개의 위상 없는 데이터) 를 사용하여 노이즈가 없는 상태에서 단위 원을 성공적으로 재구성한다.
• 노이즈 강건성: 연날리기 모양의 장애물을 이용한 실험은 5% 및 10% 의 노이즈 수준에 대한 알고리즘의 강건성을 입증한다. 재구성된 형상은 정확하고 뚜렷하게 유지된다.
• 다중 장애물: 이 방법은 노이즈가 존재하는 상황에서도 두 개의 서로 다른 장애물 (원형과 연날리기 모양) 을 동시에 효과적으로 위치시키고 형상을 파악하는 것으로 나타났다.

의의 및 주장
본 논문은 제안된 직접 영상화 방법이 반복적 방법에 비해 두 가지 주요 이점을 제공한다고 주장한다:

1. 산란체에 대한 사전 기하학적 가정의 필요성을 제거한다.
2. 직접 산란 문제를 반복적으로 푸는 계산적 부담을 피한다.

저자들은 그들의 접근 방식이 산란장부터 위상 없는 총장에 이르는 단일 유형의 측정 데이터를 사용하면서 더 간단한 헬름홀츠 기본 해에 의존하는 바이하모닉 역 산란 문제에 RTM 기반 직접 영상화를 적용한 최초의 사례임을 강조한다. 분해능 분석은 이 방법이 장애물 경계를 식별하기 위한 명확한 대조도를 제공함을 확인시켜 주며, 수치 실험은 복잡한 기하학적 구조를 재구성하는 효율성과 노이즈 환경에서의 안정성을 검증한다. 본 연구는 직접 영상화 방법의 적용 범위를 2 차 파동 방정식에서 4 차 바이하모닉 파동 방정식으로 확장한다.