[논문 리뷰] Analytical Solution for Wave Propagation in Stratified Acoustic/Porous Media. Part II: the 3D Case

Julien Diaz, Abdelaâziz Ezziani|ArXiv.org|2008. 07. 25.

Seismic Imaging and Inversion Techniques참고 문헌 9인용 수 3

한 줄 요약

이 논문은 Cagniard-de Hoop 방법을 사용하여 이중층 음향/다공성 탄성체에서 3차원 파동 전파를 분석적으로 해결한다. 이는 이전의 2차원 연구를 3차원으로 확장한 것으로, 반사파 및 투과파(각각 P, S, Ps 유형 포함)에 대한 명시적 표현을 유도하여 경계에서의 파동 거동을 정밀하게 모델링할 수 있게 하며, 점원천과 가우시안 시간 함수를 사용한 수치 시뮬레이션을 통해 검증된다.

ABSTRACT

We are interested in the modeling of wave propagation in an infinite bilayered acoustic/poroelastic media. We consider the biphasic Biot's model in the poroelastic layer. The first part is devoted to the calculation of analytical solution in two dimensions, thanks to Cagniard de Hoop method. In this second part we consider the 3D case.

연구 동기 및 목표

• 3차원 다층 구조로 이루어진 음향층 위에 다공성 탄성층이 있는 매질에서 일시적 파동 전파에 대한 완전한 분석적 해를 개발하기 위해.
• 다공성 탄성체 파동 문제에 대해 Cagniard-de Hoop 방법을 2차원에서 3차원으로 확장하여 모든 파동 유형을 독립적으로 계산할 수 있도록 하기 위해.
• 모든 파동 모드(P, S, Ps, Pf 포함)에 대해 그린 함수와 변위 성분에 대한 명시적 표현을 유도하기 위해.
• 점원천에 제5도 도함수를 갖는 가우시안 시간 함수를 사용한 수치 시뮬레이션을 통해 분석적 해를 검증하기 위해.
• 특히 복잡한 다층 구조에서 파동 전파를 시뮬레이션하는 수치 코드의 기준이 되는 해를 제공하기 위해.

제안 방법

• 다공성 탄성층에서는 Biot의 방정식을, 음향층에서는 표준 파동 방정식을 사용하여 3차원 파동 문제를 수립한다.
• 유도된 편미분 방정식계를 주파수-파수 도메인에서 해결하기 위해 Cagniard-de Hoop 기법을 적용한다.
• 경로 적분과 안장점 방법을 사용하여 라플라스-푸리에 변환을 역행렬화함으로써 변위와 압력에 대한 그린 함수를 유도한다.
• 변분적 접근을 사용하여 경계를 통과하는 파동의 최적 전파 시간을 최소화함으로써 투과파(예: Ps)의 허구적 및 실제 도착 시간을 계산한다.
• 최적의 경계 반사점에 대한 4차 다항방정식을 풀어 파동 경로 기하학을 결정한다.
• 중점 구적법을 사용하여 수치적 적분을 수행하고, 결과를 소스 시간 함수와 콘볼루션하여 시간 도메인 신호를 산출한다.

실험 결과

연구 질문

• RQ1 다공성 탄성체에서의 파동 전파에 대해 Cagniard-de Hoop 방법을 2차원에서 3차원으로 어떻게 확장할 수 있는가?
• RQ2 3차원 음향/다공성 탄성체 이중층 매질에서 반사파 및 투과파장(P, S, Ps, Pf)에 대한 분석적 표현은 무엇인가?
• RQ3 다양한 파동 유형(Ps 및 Pf 포함)의 파동 속도와 도착 시간은 파수 및 경계 기하학에 따라 어떻게 달라지는가?
• RQ4 허구적 속도는 Cagniard-de Hoop 방법에서 어떤 역할을 하는가? 그리고 이는 시간 도메인 해의 정확성에 어떤 영향을 미치는가?
• RQ5 어떻게 하면 알려진 소스 함수와 수신기 구성 조건에 대해 분석적 해를 수치적으로 검증할 수 있는가?

주요 결과

• 분석적 해는 3차원에서 P, S, Ps, Pf 모든 파동 유형을 별도로 정확하게 계산하며, 시간 도메인 신호에서 명확한 도착 시간이 관측된다.
• 500 m 상단에 위치한 점원천과 중심 주파수 15 Hz인 가우시안 소스 함수의 제5도 도함수를 사용하여 수치적 검증이 수행되었다.
• 음향층과 다공성 탄성층 양쪽 모두에서 경계로부터 533 m 떨어진 곳에 위치한 수신기들은 Fig. 9와 10에서 확인되듯이 Ps 파동 및 Pf 파동을 포함한 파동 도착 시간이 명확하게 분리되어 있음을 보였다.
• 허구적 도착 시간 함수 $\tilde{t}_0(q_y)$ 는 $\mathbb{R}^+$ 에서 짝함수이면서 증가함수이며, $\mathbb{R}^+$ 에서 $\mathbb{R}^+$ 로의 전단사 함수이므로 Cagniard-de Hoop 방법에서 안정적인 역행렬화를 보장한다.
• 이 방법은 모든 파동 모드의 군파면과 진폭을 정확히 포착하며, Ps 파동은 그 군속도가 낮아 P 및 S 파동보다 후행하여 도착한다.
• 수신기에서 계산된 변위 성분은 입사파, 반사파, 투과파에 해당하는 명확한 파aket을 보이며, 이는 해의 물리적 일관성을 확인시킨다.