3차원 층상 음향·다공성 매질 파동 전파의 정확한 해법

초록

본 논문은 무한히 넓은 두 층(음향층과 비오트 다공성층)으로 구성된 매질에서 3차원 파동 전파를 해석적으로 풀기 위해 Cagniard‑de Hoop 기법을 확장한다. 비오트 이론에 기반한 다공성층의 복합 파동(빠른·느린 압축파와 전단파)과 음향층의 단일 압축파 사이의 경계조건을 엄격히 적용하여, 라플라스·푸리에 변환 후 복소 적분을 통해 Green 함수 형태의 해를 도출한다. 결과는 파동의 반사·투과 계수와 전파 시간‑지연을 명시적으로 제공하며, 2D 해와 일관성을 확인한다.

상세 분석

이 연구는 기존 2차원 해법을 3차원으로 일반화함으로써, 층상 구조에서 실제 현장 적용이 가능한 해석적 도구를 제공한다는 점에서 큰 의의를 가진다. 먼저, 음향층은 비점음원에 대한 표준 파동 방정식(헬름홀츠 방정식)으로 기술하고, 다공성층은 비오트의 전이 방정식 두 개(고체 변위와 유체 흐름)으로 모델링한다. 여기서 중요한 점은 비오트 매개변수(탄성계수, 비투과성, 동점성 등)를 복소수 형태로 취급해, 저주파·고주파 영역 모두에서 정확한 파동 속도와 감쇠를 반영한다는 것이다.

경계조건은 두 층 사이의 연속성(압력·법선 변위 연속)과 응력·유량 연속을 동시에 만족하도록 설정한다. 이를 라플라스 변환(t→s)과 2차원 푸리에 변환(x,y→k₁,k₂)으로 변환하면, 각 파동 모드에 대한 특성 방정식이 복소 평면에서 다중 근을 갖는 형태가 된다. 저자들은 Cagniard‑de Hoop 기법을 이용해 복소 적분 경로를 실시간 파동 전파 시간 변수 τ로 매핑함으로써, 복소 근들을 실제 물리적 파동 전파 경로와 연결한다. 이 과정에서 빠른 압축파(P₁), 느린 압축파(P₂), 전단파(S) 각각에 대한 시간‑지연 함수와 진폭 함수를 명시적으로 도출한다.

특히 3차원에서는 파동 전파가 구면파 형태를 띠므로, 방사형 거리 r=√(x²+y²+z²)와 각도 θ, φ에 대한 의존성을 포함해야 한다. 저자들은 구면 좌표계에서의 푸리에 역변환을 수행할 때, 베셀 함수와 구면 조화 함수를 활용해 각 파동 모드의 방사 패턴을 정확히 표현한다. 결과적으로, 반사·투과 계수는 입사 각도와 주파수에 따라 복소수 형태로 제공되며, 이는 전통적인 수치 해법에서 발생하는 불안정성을 회피한다.

또한, 해의 수렴성을 보장하기 위해 복소 평면에서의 브랜치 컷 선택과 적절한 경로 변형을 상세히 논의한다. 이는 특히 느린 압축파가 강한 감쇠를 보이는 고주파 영역에서 중요한데, 적절한 브랜치 선택 없이는 물리적으로 의미 없는 해가 도출될 수 있다. 저자들은 이러한 수학적 세부사항을 체계적으로 정리함으로써, 해가 모든 파라미터 범위에서 일관되게 적용될 수 있음을 증명한다.

마지막으로, 2D 해와의 비교 실험을 통해 3D 해가 정확히 2D 해의 축대칭 특수 경우를 포함한다는 점을 확인한다. 이는 모델의 일반성과 확장성을 뒷받침한다. 전반적으로, 본 논문은 복잡한 층상 매질에서 3차원 파동 전파를 해석적으로 다루는 최초의 시도 중 하나이며, 향후 지진 탐사·해양 음향·비파괴 검사 등 다양한 분야에 직접 활용될 수 있는 기반을 제공한다.