상수 Q와 안정 확률분포를 갖는 지진 파동 전파

초록

본 논문은 품질인자 Q가 주파수에 독립적인 경우, 1차원 지진 파동이 시간에 대한 분수 차수 미분 방정식으로 기술됨을 보인다. 이 방정식은 열방정식과 파동방정식 사이를 연결하며, Cauchy와 Signalling 문제의 기본 해는 Wright 유형 전체함수로 표현된다. 해는 완전 점성 유체와 완전 탄성 고체의 극한 형태 사이의 중간 특성을 보이며, 거의 탄성적인 경우를 중점적으로 분석한다. 또한, 이러한 기본 해는 분수 차수에 의해 결정되는 안정 확률분포와 직접적인 연관성을 가진다.

상세 분석

논문은 먼저 선형 점탄성 이론을 Riemann‑Liouville·Caputo 분수 미분 체계에 연결한다. 점탄성 매질의 크리프 컴플라이언스 J(t)가 초기값 J₀와 파워‑법칙 항 t^ν (0<ν≤1)으로 전개될 때, 라플라스 변환을 적용하면 파동‑확산 방정식
∂^{2β}w/∂t^{2β}=D ∂²w/∂x², 2β=2−ν (1/2≤β<1)
이 도출된다. β=1/2이면 전통적인 열방정식, β=1이면 고전적인 파동방정식이 회복된다. 이 식은 “분수 확산‑파동 방정식”이라 명명되며, 시간 차수가 1보다 작아 비국소적 메모리 효과를 내포한다.

Cauchy 문제(전역 초기값)와 Signalling 문제(반무한 경계) 각각에 대해 라플라스 변환 후 2차 공간 미분 방정식을 풀어 Green 함수의 변환형을 구한다. 역변환을 수행하면 두 기본 해는 동일한 similarity 변수 r=|x|/(√D t^β) 를 갖는 두 개의 보조 함수 F(r;β)와 M(r;β) 로 표현된다. 이 보조 함수들은 복소 적분 형태의 Wright 전체함수이며, r→0,∞ 에서의 asymptotic은 각각 급격한 피크와 장거리 꼬리를 제공한다. 특히, β가 1에 가까울수록 해는 매우 좁고 뾰족해져 거의 탄성 매질에서 관측되는 지진 펄스와 일치한다.

또한 저자들은 Green 함수 사이에 2β x G_c = t G_s 라는 reciprocity 관계를 증명한다. 이를 이용해 G_c와 G_s 를 서로 변환할 수 있으며, 이는 해석적 연속성 및 수치적 구현에 유용하다.

마지막으로, Signalling 문제의 Green 함수 G_s(x,t;β) 를 확률밀도 함수와 동일시한다. β에 따라 안정 확률분포의 지수 α=2β 가 결정되며, β=1/2 일 때는 Lévy‑α‑stable 분포(α=1)와 일치한다. 따라서 분수 차수는 물리적 감쇠(Q)와 확률적 안정성 사이의 매개변수 역할을 한다. Q와 ν 사이의 관계 Q⁻¹=tan(πν/2) 로부터, 실제 지진 매질에서 Q≈10³ 정도이면 ν≈0.001 수준으로 거의 탄성에 가깝고, 이에 대응하는 β≈0.9995 로서 해는 거의 파동형태를 유지한다.

이와 같이 논문은 분수 미분 방정식이 지진 파동 전파의 물리적, 수학적, 확률적 측면을 통합적으로 설명함을 보여준다.