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

초록

본 논문은 무한히 넓은 두 층으로 구성된 다공성 매질에서 Biot 이론에 기반한 파동 전파를 3차원으로 해석한다. Cagniard‑De Hoop 기법을 적용해 라플라스‑푸리에 변환 후 복소 평면에서 적분 경로를 변형함으로써, 평면 인터페이스를 통과하는 전파·반사·투과 파동의 정확한 해를 도출한다. 결과는 빠른·느린 압축파와 전단파 각각에 대한 Green 함수 형태로 제시된다.

상세 분석

이 연구는 Biot의 2상 연속체 모델을 기반으로, 층상 구조를 갖는 무한 매질에서 발생하는 파동 전파 현상을 3차원적으로 해석한다는 점에서 학문적·실용적 의의를 가진다. 기존 2차원 해법은 복잡한 경계조건과 파동 모드 결합을 충분히 포착하지 못했으나, 본 논문은 Cagniard‑De Hoop(CDH) 기법을 3차원으로 확장함으로써 이러한 한계를 극복한다. CDH 기법은 라플라스 변환 후 복소 시간 평면에서 적절한 경로(Cagniard 경로)를 선택해 역변환을 수행함으로써, 시간‑공간 영역에서의 해를 직접적으로 얻을 수 있게 한다. 논문은 먼저 Biot 방정식을 라플라스‑푸리에 변환하여 6차원 행렬 형태의 고유값 문제를 만든다. 여기서 고유값은 빠른 압축파(P_f), 느린 압축파(P_s), 전단파(S) 세 종류로 분리되며, 각 파동은 복소 파수와 전파 속도로 표현된다. 층간 경계에서는 연속성 조건(변위·응력·압력 연속)을 적용해 반사·투과 계수를 행렬 형태로 도출한다. 특히, 3차원에서는 파동의 편향각과 방위각이 추가 자유도를 제공하므로, 반사·투과 계수는 각도 의존성을 갖는 복합 행렬로 나타난다. 이후 Cagniard 경로를 이용해 복소 파수 적분을 실시간 변수 t와 거리 r에 대한 실함수 적분으로 변환한다. 이 과정에서 각 파동 모드별로 특이점(극점·분기점)을 정확히 처리해, 급격한 전이 현상(예: 앞선 파와 후속 파)도 수식적으로 기술한다. 최종적으로 얻어진 Green 함수는 시간‑공간 영역에서의 파동장(변위·압력·응력)을 명시적으로 제공하며, 이는 수치 시뮬레이션(예: FEM, FDTD)과 비교해 높은 정확도와 계산 효율성을 보인다. 논문은 또한 파라미터(다공성, 투과성, 비탄성 계수)의 변화가 파동 전파에 미치는 영향을 정량적으로 분석하고, 실험적 데이터와의 일치성을 검증한다. 이러한 분석은 지구물리학적 탐사, 지진공학, 해양 구조물 설계 등 다공성 매질을 다루는 분야에서 실시간 파동 예측 및 역문제 해결에 직접적인 활용 가능성을 제시한다.