층상 다공성 매질 지진파를 통한 탄성파라미터 추정

초록

본 연구는 Biot 이론을 기반으로 한 전파 모델링과 전파 반사율 시뮬레이션을 이용해 층상 다공성 매질의 지진파 형태를 전방 문제로 해결하고, 파라미터 민감도와 일반화 최소제곱 역문을 결합한 전파 전파역전(FWI) 절차를 개발한다. 합성 데이터 실험을 통해 공극률과 응집 파라미터가 가장 민감하게 반응함을 확인했으며, 투과성은 지진 응답에 거의 영향을 미치지 않는다는 결론을 얻었다. 다중 파라미터 역전에서는 파라미터 간 강한 결합으로 인한 비식별성이 발생하므로, 암석 물리학 법칙에 기반한 복합 파라미터 혹은 매질 교란(예: 가스 주입) 형태로의 역전을 제안한다.

상세 요약

본 논문은 Biot 이론을 층상 구조에 적용한 전파 모델링을 구현함으로써, 유체가 포화된 다공성 매질에서 발생하는 1차 및 2차 파동(압축파와 전단파)의 전파 특성을 정밀하게 재현한다. 구현된 프로그램은 전통적인 reflectivity 방법을 확장하여, 각 층의 탄성계수, 밀도, 공극률, 투과성, 응집 파라미터(비탄성 계수) 등을 입력 변수로 받아 복합 전파 응답을 계산한다. 전방 모델링 단계에서 파라미터별 민감도 연산자를 도출했으며, 이는 파라미터 변화가 지진파 진폭·위상에 미치는 영향을 정량화한다. 민감도 분석 결과, 공극률(porosity)과 응집 파라미터(consolidation parameter, 흔히 비탄성 계수로 표현)는 전파 속도와 감쇠에 크게 기여해 정역학적·동역학적 응답을 모두 크게 변화시킨다. 반면 투과성(permeability)은 저주파 대역에서 거의 영향을 미치지 않으며, 고주파에서도 파형 변형에 미치는 기여가 제한적이다. 이러한 특성은 역전 과정에서 투과성을 독립적으로 추정하기 어려움을 의미한다.

역전 알고리즘은 일반화 최소제곱(Generalized Least Squares, GLS) 프레임워크를 채택하고, 앞서 정의한 민감도 연산자를 Jacobian 행렬에 삽입한다. 반복적인 선형화와 업데이트 과정을 통해 모델 파라미터를 점진적으로 보정하며, 각 반복마다 합성 파형과 관측 파형 사이의 L2-norm 오차를 최소화한다. 이때 파라미터 간 상관관계가 높은 경우(예: 공극률과 고체 탄성계수, 또는 유체 압축성계수와 고체 압축성계수) Jacobian의 조건수가 크게 악화되어 역전이 수렴하지 않거나 비물리적 해가 도출될 위험이 있다. 저자들은 이러한 문제를 완화하기 위해 두 가지 전략을 제시한다. 첫째, 암석 물리학 관계식(예: Gassmann 방정식)을 이용해 파라미터를 결합한 복합 파라미터(예: 효과 탄성계수)를 직접 역전한다. 둘째, 매질의 절대값이 아니라 교란(perturbation) 형태, 즉 가스 주입에 따른 공극률 변화와 같은 차분값을 역전함으로써 민감도가 높은 파라미터만을 효과적으로 추출한다.

합성 실험에서는 알려진 모델을 기반으로 잡음이 추가된 파형을 생성하고, 역전 알고리즘을 적용하였다. 결과는 공극률과 고체·유체 탄성계수를 정확히 복원했으며, 투과성은 초기값에 크게 의존해 회복이 어려웠다. 또한 파라미터 결합이 강한 경우(예: 고체 탄성계수와 유체 압축성계수)에는 복합 파라미터 접근법이 개별 파라미터 역전보다 안정적인 결과를 제공함을 확인했다. 전반적으로, 본 연구는 Biot 기반 전파 모델링과 민감도 기반 GLS 역전이 다공성 매질의 주요 탄성 파라미터를 추정하는 데 유효함을 입증했으며, 파라미터 결합 문제를 해결하기 위한 실용적인 방안을 제시한다.