맥스웰형 점탄성 유체가 다공성 매질 파동 전파에 미치는 영향

초록

본 연구는 분수 파생 마크스웰 점탄성 유체가 포어엘라스틱 매질에 포화될 때, 파동 전파 속도와 감쇠에 미치는 영향을 분석한다. 분수 미분 기반 응력‑변형률 관계로부터 점성 소산 함수를 도출하고, 이를 이용해 P‑파와 S‑파의 분산·감쇠 식을 얻었다. 수치 예시에서 점탄성 유체는 전통적 점성 유체와 달리 전파 속도 진동과 소산‑탄성 전이 현상을 보이며, 실험적 관측과 일치한다.

상세 분석

이 논문은 전통적인 비점탄성(뉴턴) 유체 대신, 분수 차수(α)와 시간 상수(τ)로 정의되는 마크스웰형 점탄성 유체를 도입함으로써 포어엘라스틱 매질 내 파동 전파 메커니즘을 재정의한다. 핵심은 응력‑변형률 관계 σ = E ε + η D^α ε (여기서 D^α는 분수 미분 연산자)에서 유도된 점성 소산 함수 Ψ = η |ω|^α sin(πα/2)이다. 이 함수는 전통적인 점성 손실(η |ω|)과 달리 주파수에 대한 비선형 의존성을 갖으며, α가 0에 가까울수록 탄성에, 1에 가까울수록 점성에 가까운 거동을 보인다.

소산 함수를 포어엘라스틱 방정식에 삽입하면, 복소 파수 k = k_r + i k_i가 주파수 ω에 대해 다음과 같이 변한다. P‑파와 S‑파 각각에 대해 실재 파수 k_r는 전파 속도 v = ω/k_r를, 허수 파수 k_i는 감쇠 계수 α_att = k_i를 제공한다. 분석 결과, α가 0.3~0.7 사이일 때, 특히 중간 주파수 대역(10–100 Hz)에서 v가 비주기적으로 진동하고, 감쇠는 급격히 감소한다. 이는 소산‑탄성 전이 현상으로, 유체가 저주파에서는 점성 손실을, 고주파에서는 탄성 저장을 주도한다는 것을 의미한다.

또한, 유체 속도 u_f = (φ/η_eff)∇p와 같은 Darcy‑type 관계를 분수 차수와 결합해 유도함으로써, 실험에서 관측된 ‘속도 진동 현상’을 정량적으로 재현한다. 여기서 η_eff = η |ω|^{α-1} sin(πα/2)는 유동 저항을 주파수 의존적으로 변환한다.

수치 시뮬레이션에서는 전형적인 석유·가스 저류층 파라미터(다공성 φ = 0.15, 전단 탄성률 μ = 3 GPa 등)를 사용했으며, α와 τ를 변화시켜 파동 전파 특성을 매핑했다. 결과는 다음과 같다. (1) α가 작을수록 P‑파와 S‑파의 전파 속도가 고주파에서 급격히 상승한다. (2) τ가 클수록 전이 주파수 대역이 낮아져, 저주파에서도 탄성 효과가 두드러진다. (3) 감쇠 계수는 α가 0.5 근처에서 최소값을 보이며, 이는 점탄성 유체가 최적의 에너지 전달 효율을 갖는 조건을 시사한다.

이러한 발견은 기존의 비점탄성 유체 모델이 설명하지 못한 ‘속도 진동’과 ‘감쇠 감소’를 물리적으로 설명한다는 점에서 학술적 의미가 크다. 또한, 비정상적인 유동 거동을 보이는 비전통적 유전(예: 셰일 가스, 탄산염 저장소)에서 파동 기반 탐사 해석에 직접 적용 가능하다.