Retrieving the Seismic Moment Release Rate of Small Seismic Events Using the Adjoint Method - Theory and Numerical Validation

지구가 흔들릴 때, 아주 작은 진동조차도 지표면이 어떻게 파괴되었는지에 대한 숨겨진 이야기를 담고 있습니다. 지진학자들은 오랫동안 지진이 어디서 시작되었고 얼마나 많은 에너지를 방출했는지는 정확히 짚어낼 수 있었지만, 그 에너지가 시간에 따라 정확히 어떻게 방출되었는지를 이해하는 것은 어려운 난제로 남아 있었습니다. 이 이야기는 소스 시간 함수(source time function)라고 불리는 것, 즉 첫 균열부터 마지막 멈춤까지의 파쇄 리듬을 기록한 것에 의해 전해집니다. 이 리듬을 아는 것은 매우 중요한데, 왜냐하면 그것이 파쇄의 속도, 이동 방향, 그리고 지하 깊은 곳에서 발생한 응력 변화를 드러내기 때문입니다. 그러나 지표면에서 기록된 무질서한 파동으로부터 이 리듬을 추출하는 것은 매우 까다로운 일입니다. 전통적인 방법들은 더 작고 유사한 지진을 찾아 참조용으로 사용하는 방식(주관적이고 부정확할 수 있는 과정)에 의존하거나, 일상적인 용도로 쓰기에는 너무 느린 대규모의 컴퓨터 집약적 시뮬레이션을 실행하려고 시도합니다.

새로운 연구에서, 폴란드 과학 아카데미 지질학 연구소의 보이치에흐 뎀브스키(Wojciech Dębski)와 카밀 바시케비치(Kamil Waśkiewicz) 연구진은 이러한 숨겨진 리듬을 듣는 다른 방법을 제안했습니다. 그들은 소리와 지진파가 시간이 앞으로 흐르는 것과 마찬가지로 시간을 거슬러 뒤로 흐를 수 있다는 단순한 물리적 개념에 기반한 기술을 개발했습니다. 지표면의 기록을 가져와 시간을 역전시킨 후, 컴퓨터 모델을 통해 땅속으로 다시 보냄으로써 파동이 자연스럽게 지진이 시작된 지점으로 다시 집중되도록 하는 것입니다. 이 역전된 파동들이 수렴함에 따라, 에너지 방출의 원래 패턴을 재구성합니다. 연구팀은 이 아이디어를 광업 활동 등으로 유발된 작은 인위적 지진의 컴퓨터 시뮬레이션, 즉 **노이즈가 없는 합성 데이터(noise-free synthetic data)**를 사용하여 테스트했습니다. 그들의 결과는 이 방법이 복잡한 수학적 역산이나 참조 지진을 찾을 필요 없이 상세한 파쇄 타이밍을 성공적으로 복구할 수 있음을 보여줍니다. 이는 복잡한 수학적 문제를 단순한 물리적 과정으로 전환하여, 작은 지진 사건의 역학을 이해하는 더 빠르고 견고한 방법을 제공합니다.

이 새로운 접근법의 핵심은 파동 물리학의 근본적인 특성인 시간 가역성(time reversibility)에 기초합니다. 에너지가 열이나 마찰로 손실되지 않는 완전 탄성 매질에서, 파동의 움직임을 설명하는 방정식은 시간이 앞으로 흐르든 뒤로 흐르든 동일하게 작동합니다. 연못에 돌이 떨어지는 영상을 상상해 보십시오. 만약 영상을 역재생한다면, 잔물결은 다시 안쪽으로 모여들며 돌이 물 밖으로 튀어 오를 것입니다. 실제 세상에서 시간을 실제로 되돌릴 수는 없지만, 수학은 이 과정을 시뮬레이션하는 것을 허용합니다. 연구진은 이 개념을 사용하여 TRSTF라고 부르는 알고리즘을 만들었습니다. 과정은 여러 관측소에서 기록된 지진파를 가져오는 것으로 시작됩니다. 그 후 이 파동들을 시간적으로 역전시켜 지구 지각의 컴퓨터 모델로 다시 주입합니다. 역전된 파동들은 모델을 통과하며 경계면에서 반사되고 복잡한 지질 구조와 상호작용하며, 결국 원래의 진원 위치로 수렴하게 됩니다.

이 역전된 파동들이 진원에서 만날 때, 이들은 서로 간섭합니다. 정확한 지진 발생 위치에서 파동들은 완벽하게 정렬되어 강력하고 집중된 신호를 만들어냅니다. 다른 어떤 위치에서도 파동들은 서로 어긋나 도착하여 서로를 상쇄시키고, 오직 희미한 노이즈만을 남깁니다. 이 집중 효과가 이 방법의 핵심입니다. 초점에 나타나는 신호는 지진 에너지 방출의 직접적인 재구성인 소스 시간 함수입니다. 연구진은 이 기술이 노이즈가 없는 합성 실험에서, 심지어 지구 내부 모델이 완벽하게 알려지지 않은 경우에도 놀라울 정도로 잘 작동한다는 것을 발견했습니다. 컴퓨터 실험에서 그들은 실제 모니터링 설정들을 모사하기 위해 세 가지 서로 다른 시나리오를 시뮬레이션했습니다. 하나는 센서가 광업 활동 수준 근처에 배치된 경우, 다른 하나는 이벤트 위의 지표면에 센서가 있는 경우, 그리고 세 번째는 다른 지표면 배치를 가진 경우였습니다. 모든 경우에서, 이 방법은 원래의 파쇄 패턴 형태를 성공적으로 복구했습니다.

연구는 또한 지구의 음속 모델의 오류가 결과에 어떤 영향을 미치는지 탐구했습니다. 현실 세계에서 지질학자들은 모든 지점에서의 정확한 지진파 속도를 알지 못하며, 추정치에 의존합니다. 연구진은 모델에 무작위 오류와 블록 형태의 변화를 도입하여, 속도 모델이 약 90% 정도만 정확한 상황을 시뮬레이션함으로써 자신들의 방법을 테스트했습니다. 결과는 부정확한 속도 모델의 사용이 재구성된 소스 시간 함수에 계통적 오차를 유발한다는 것을 나타냈습니다. 이 방법은 무작위로 흩어진 오류보다 조직된 블록 형태의 오류에 더 민illary(민감)했습니다. 그러나 이러한 오류에도 불구하고, 복구된 신호는 식별 가능한 상태를 유지했으며 유용한 정보를 제공했습니다. 이는 이 방법이 불완전한 모델을 사용하더라도 유용성을 유지하지만, 정확도가 모델의 부정확함으로부터 자유롭지는 않음을 보여줍니다.

가장 놀라운 발견 중 하나는 센서의 배치가 결과의 품질에 어떻게 영향을 미치는지였습니다. 연구진은 연구 대상 구역의 경계에 가까운 곳에 센서를 배치하는 것이 오히려 재구성을 돕는다는 것을 발견했습니다. 이는 경계로부터의 반사가 추가적인 가상 센서 역할을 하여 시스템에 더 많은 정보를 제공하기 때문입니다. 센서가 모델의 상단 가장자리 근처에 배치되었던 테스트 케이스 중 하나에서는, 원시 데이터가 간섭으로 인해 더 왜곡되어 보였음에도 불구하고 재구성된 신호가 다른 설정보다 더 날카롭고 정확했습니다. 이 반직관적인 결과는 이 방법이 파동장의 복잡성을 활용하여, 반사를 노이즈로 취급하는 대신 이를 유리하게 이용한다는 점을 시사합니다.

연구진은 또한 재구성된 신호의 구체적인 특성, 예를 들어 파쇄가 얼마나 빨리 시작되었고 얼마나 오래 지속되었는지를 측정했습니다. 그들은 재구성된 신호가 원래 신호보다 약간 더 넓어졌는데, 이는 아마도 시뮬레이션에 사용된 센서의 수가 제한적이었기 때문일 것이라고 판단했습니다. 하지만 전체적인 형태와 타이밍은 보존되었습니다. 파쇄 지속 시간을 추정하는 능력은 매우 중요한데, 이는 과학자들이 사건 동안 방출된 스트로크 드롭(stress drop, 응력 강하)을 계산하는 데 도움을 주기 때문입니다. 이 정보는 광업이나 유체 주입 등에 의해 유발된 지진의 역학을 이해하는 데 필수적입니다. 연구는 이 새로운 방법이 기존 기술과 대등하거나 더 나은 정확도로 이러한 매개변수들을 추정할 수 있음을 보여주었으나, 훨씬 적은 계산량으로 이를 수행했습니다.

수천 번의 시뮬레이션과 최적의 적합점을 찾기 위한 복잡한 최적화 루프를 실행해야 하는 전체 파형 역산(full-waveform inversion)과 달리, 이 새로운 방법은 모델을 통해 역전된 파동을 단 한 번 통과시키는 것만을 요구합니다. 이는 계산 효율성을 높여, 작은 모델의 경우 표준 데스크톱 컴퓨터에서, 지역적 규모의 경우 대형 슈퍼컴퓨터에서 실행할 수 있게 합니다. 또한 참조 지진을 선택할 필요가 없으므로, 전통적인 방법에서 인간의 편향과 불확실성을 초래할 수 있는 단계를 피할 수 있습니다. 연구진은 이 접근법이 특히 소규모 지진 사건에 매우 적합하다고 강조하는데, 이는 진원이 단일 점으로 취급될 만큼 작기 때문입니다. 그들은 이 방법이 지면이 완벽한 탄성체처럼 행동한다는 것(즉, 에너지를 흡수하지 않음)과 진원이 파동의 파장에 비해 작다는 것을 가정하고 있음을 인정합니다. 이는 많은 광업 관련 사건들에 대해 합리적인 가정이지만, 연구진은 에너지 손실이 상당한 더 복적인 실제 환경에서 이 방법이 어떻게 작동하는지 이해하기 위해 추가적인 연구가 필요하다고 언급했습니다.

논문은 이 기술이 지진학의 표준 도구가 될 가능성을 강조하며 마무리됩니다. 소스 시간 함수를 회수하는 빠르고 신뢰할 수 있으며 물리적으로 직관적인 방법을 제공함으로써, 이 기술은 이전에는 상세한 특성 파악이 너무 어려웠던 작은 지진 사건들의 정기적인 분석을 위한 길을 열어줍니다. 연구진은 이 방법이 기존의 모니터링 시스템에 통합되어 유발 지진의 본질에 대한 즉각적인 통찰력을 제공할 수 있다고 제안합니다. 현재의 연구는 컴퓨터 시뮬레이션에 국한되어 있지만, 이론적 토대는 견고하며 결과는 유망합니다. 저자들에 따르면 다음 단계는 실제 지진 데이터를 사용하여 이 방법을 테스트하고, 그 한계와 불확실성에 대한 이해를 정교화하는 것입니다. 현재로서 이 연구는 시간 역전이 지구의 과거를 들여다보는 강력한 렌즈가 되어, 지면이 어떻게 파괴되는지에 대한 숨겨진 이야기를 밝혀낼 수 있다는 것을 보여주는 개념 증명(proof of concept) 단계에 있습니다.