Optimization and robustness of cost-efficient seismic arrays for Newtonian noise cancellation at the Einstein Telescope

기술 요약: 아인슈타인 망원경의 뉴턴 노이즈 제거를 위한 비용 효율적인 지진계 어레이의 최적화 및 강건성

문제 정의
제3세대 지하 중력파 검출기를 목표로 하는 아인슈타인 망원경(Einstein Telescope, ET)은 3 Hz까지의 저주파수 대역에서 중력파를 검출하고자 한다. 이 저주파수 대역에서는 뉴턴 노이즈(Newtonian Noise, NN)가 지배적인 노이즈원이 될 것으로 예상된다. NN은 주변 암석의 밀도 변화에 의해 발생하는 중력적 힘에서 기인하며, 주로 지진체파(P파 및 S파)에 의해 발생한다. 이를 완화하기 위해서는 보조 지진계 어레이를 사용하여 지진파장을 측정하고, 이를 통해 인터로미터 출력에서 NN을 예측하여 제거해야 한다. 그러나 센서를 배치하기 위해 시추공(borehole)을 뚫는 것은 상당한 비용 상승 요인이 된다. 따라서 여러 개의 센서를 하나의 시추공에 배치하거나 기존 인프라를 활용하는 구성을 탐색함으로써, 비용을 최소화하면서도 완화 성능을 극대화할 수 있는 센서 위치 최적화가 과제이다.

방법론
저자들은 ET의 한 코너(네 개의 상관된 거울)에 대한 최적의 지진계 위치를 결정하기 위해 Ref. [12]에서 확립된 최적화 프레임워크를 사용한다. 이 최적화는 정적 위너 필터(stationary Wiener filters)를 사용하여 NN 신호의 평균 잔차를 최소화한다.

• 최적화 알고리즘: 초기화를 위해 입자 군집 최적화(Particle Swarm Optimization, PSO)를 사용하고, 이후 미세 조정을 위해 디퍼렌셜 옵티마이저인 Adam을 사용한다.
• 성능 지표: 완화 계수 $M(f)$는 네 개의 거울 중 최대 잔차의 제곱근 역수로 정의되며, 이는 진폭 감소를 나타내는 보수적인 지표 역할을 한다.
• 테스트된 시나리오:
  1. 벤치마크: 시추공당 하나의 지진계를 배치한 어레이로, 총 시추공 수($N = 20$ ~ $70$)를 변화시킨다.
  2. 다중 센서 시추공: 각 시추공에 $n = 1, 2, 3, 5, 10$개의 지진계가 포함된 $N=20$개의 시추공 어레이. 최적화 제약 조건은 시추공 내에서 $x$와 $y$ 위치는 가변적으로 허용하되, $z$ 위치는 독립적으로 가변하도록 허용한다.
  3. 인프라 통합: 코너로부터 최대 5km까지의 ET 터널(암) 내부에 100m 간격으로 고정된 위치의 지진계 100개를 추가하며, 이들의 특정 위치는 추가 최적화하지 않는다.
  4. 강건성 분석: 모든 세 공간 방향에 대해 무작위 가우시안 변위($\sigma = 50$ m)를 도입하여 이러한 구성들의 안정성을 테스트한다.
• 매개변수: 시뮬레이션은 신호 대 잡음비(SNR) 15, P파 분율 $p=0.2$, 그리고 특정 파동 속도($c_P = 6000$ m/s, $c_S = 4000$ m/s)를 가정한다. 어레이는 특히 $f_{opt} = 10$ Hz에 맞춰 최적화되며, 이들의 "광대역" 성능은 1~10 Hz 구간에서 평가된다.

주요 기여

1. 비용 효율적인 다중 센서 구성: 본 연구는 단일 센서가 포함된 다수의 시추공을 사용하는 어레이보다, 하나의 시추공에 여러 개의 센서를 배치하는 것이 훨씬 더 많은 시추공을 사용하는 것과 유사한 완화 능력을 달달할 수 있음을 보여준다. 이는 시추 비용 문제를 해결하기 위해 각 시추공의 효용성을 극대화한다.
2. 인프라 활용: 기존 ET 터널 인프라 내에 지진계를 배치함으로써 얻는 성능 이득을 정량화한다. 이 방식은 추가적인 시추가 필요 없으며 광대역 성능을 크게 향상시킨다.
3. 강건성 정량화: 무작위 배치 오류(최대 50 m)가 완화 성능에 미치는 영향을 상세히 분석하여, 더 큰 규모의 어레이와 다중 센서 구성이 단일 센서 어레이보다 안정성을 더 잘 유지함을 보여준다.

결과

• 광대역 성능: 10 Hz에 최적화된 어레이는 낮은 주파수(1–10 Hz)에서도 잘 일반화된다.
  • $N=20$개의 시추공에 대해, 시추공당 센서 수($n$)를 늘리면 완화 계수가 크게 향상된다. 구체적으로, $n=5$ (총 100개 센서)는 $n=1$인 $N=60$개의 시추공 어레이와 유사한 성능을 달성한다.
  • $N=20$ 구성에 터널 지진계를 추가하면 성능이 더욱 향착된다. $n=10$인 시추공과 터널 센서를 결합하면, $f \gtrsim 3$ Hz에서 완화 계수가 10을 초과하며, 피크 값은 약 34에 달한다.
• 강건성:
  • 무작위 변위는 높은 주파수(변위가 파장의 더 큰 비율을 차지하는 구간)에서 성능을 가장 크게 저하시킨다.
  • 더 큰 어레이가 더 높은 안정성을 보인다. 연구된 가장 큰 구성($N=50$ 시추공, $n=10$ 센서/시추공, 터널 센서 포함)의 경우, 무작위 변위가 있어도 $f \gtrsim 3$ Hz에서 완화 계수가 15 이상을 유지한다.
  • 터널 센서가 포함된 $N=20, n=5$ 구성의 경우, 변위 변화에도 불구하고 $f \gtrsim 3$ Hz에서 완화 계수가 7 이상을 유지한다.
  • $N=20, n=3$의 비교적 단순한 구성(터널 센서 포함)조차도 무작위 변위에도 불구하고 $f \gtrsim 3$ Hz에서 완화 계수 6 이상을 유지한다.

의의 및 주장
본 논문은 높은 주파수(예: 10 Hz)를 위해 지진 어레이를 최적화하는 것이 3 Hz까지의 효과적인 광대역 NN 완화를 제공한다고 주장한다. 주요 의의는 ET를 위한 비용 효율적인 전략을 제시했다는 점에 있다. 즉, 단일 센서가 포함된 방대한 수의 시축공을 뚫는 대신, 여러 개의 센서를 갖춘 더 적은 수의 시추공을 사용하고 기존 ET 인프라 내의 센서를 보완적으로 활용하는 방안을 제안한다.

저자들은 이러한 구성이 "유망한 광대역 성능"을 제공하며 "지진계 위치 변화에 특히 안정적"이라고 결론짓는다. 이들은 이러한 결과가 다중 센서 시추공의 기술적 타당성과 ET 인프라로의 센서 통합에 대한 추가 연구를 촉발하는 동기가 된다고 밝힌다. 또한 광대역 최적화가 더 나은 절충안을 제공할 수 있으나, 여기서 테스트된 단일 주파수 최적화 전략도 효과적임을 언급한다. 향후 연구로는 실제 지질 환경에서의 유한 요소 시뮬레이션(finite-element simulations)과 융합 어레이(지진계, 변형률계, 기울기계의 결합)에 대한 연구가 제안되었다.