지구 자유진동을 위한 등방성 확장 버거스 모델 이완 텐서 명시식 및 스펙트럼 역전

초록

본 논문은 등방성 확장 버거스 모델(q‑EBM)의 이완 텐서를 명시적으로 유도하고, 이를 이용해 구형 지구(단위 구)에서 자유진동 고유값을 계산한다. 즉시항(instantaneous)과 기억항(memory)으로 분리된 연산자를 분석하여, 즉시항 고유값과 동일한 고유함수를 갖지만 이동된 위치에 있는 고유값 군(C‑ev)을 도출한다. 또한 C‑ev의 구조를 밝히고, 관측된 C‑ev로부터 q‑EBM 파라미터를 복원하는 역전 공식을 제시한다.

상세 분석

이 논문은 지구 내부의 비탄성(anelastic) 거동을 묘사하기 위해 기존의 Burgers 모델을 확장한 quasi‑static extended Burgers model(q‑EBM)을 채택한다. 저자는 먼저 일반적인 비균질 등방성 q‑EBM에 대해 이완 텐서 G(t)를 명시적으로 구성하는 방법을 제시한다. 핵심은 변형률 텐서를 부피·전단 성분으로 분해하고, 각 성분에 대해 4차 텐서 I_m, J_m을 이용해 직교 투영을 수행함으로써 복잡한 연립 미분‑적분 방정식을 블록 행렬 형태로 변환하는 것이다. 이 과정에서 얻어진 L_0, L_1 행렬은 각각 부피와 전단 모드에 대응하며, 이들의 고유값 −τ_j, −κ_j는 음의 실수이며, 해당 고유벡터는 서로 직교한다. 이를 통해 이완 텐서의 (1,1) 성분 g_{00}(t)와 g^{0}_{00}(t)를 각각 exp(−τ_j t), exp(−κ_j t)의 선형 결합으로 표현한다(정리 2.3, 2.5).

다음 단계에서는 q‑EBM 연산자를 “즉시항 + 기억항” 형태로 재구성한다. 즉시항은 고전적인 탄성 텐서 C_0에 의해 정의되며, 기억항은 위에서 도출한 G(t)의 시간미분을 이용한 컨볼루션 형태이다. 초기 변형률이 영이라는 가정 하에, 시간 라플라스 변환을 적용하면 z∈ℂ(복소수 파라미터)와 연관된 다항식 P(z)=det(zI−A(z))의 근이 자유진동 고유값이 된다. 여기서 A(z)는 즉시항 고유값과 기억항의 가중합으로 구성된 (n+1)×(n+1) 행렬이며, 고유함수는 즉시항 고유함수와 동일함을 증명한다.

특히 저자는 “클러스터 고유값(C‑ev)”이라는 개념을 도입한다. C‑ev는 즉시항 고유값 군을 중심으로, 기억항에 의해 일정량 이동된 복합 고유값 집합이다. 이들은 다항식 P(z)의 실근 혹은 복소쌍으로 나타나며, 각 클러스터는 동일한 고유함수를 공유한다. 클러스터 구조는 τ_j, κ_j와 비례 상수 (v_j^0)^2, (q_j^0)^2에 의해 완전히 규정된다.

마지막으로 논문은 관측된 C‑ev 집합으로부터 q‑EBM 파라미터(λ_i, μ_i, η_i)를 복원하는 역전 공식을 제시한다. 역전 과정은 고유값 클러스터의 위치와 폭을 이용해 τ_j, κ_j를 추정하고, 이를 다시 L_0, L_1의 고유분해와 연결시켜 η_i와 탄성계수를 구한다. 이 절차는 기존에 존재하던 수치적 추정법보다 해석적으로 명확하며, 파라미터 식별에 필요한 최소 데이터 양을 이론적으로 제시한다.

전반적으로 이 연구는 (1) 등방성 q‑EBM의 이완 텐서를 명시적으로 구성, (2) 자유진동 고유값을 즉시항 고유함수와 연결, (3) 클러스터 고유값 구조를 분석, (4) 역전 공식을 도출함으로써 지구 내부 비탄성 모델링에 새로운 해석적 도구를 제공한다는 점에서 큰 의미를 가진다.