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

技術要約：アインシュタイン・テレスコープにおけるニュートンノイズ除去のための、コスト効率の高い地震計アレイの最適化とロバスト性

問題提起
提案されている第三世代の地下重力波干渉計であるアインシュタイン・テレスコープ（ET）は、3 Hzまでの周波数で重力波を検出することを目指している。この低周波数帯域では、ニュートンノイズ（NN）が支配的なノイズ源になると予想される。NNは、周囲の岩石の密度変動によって生じる重力的な力に起因し、主に地震波（P波およびS波）によって引き起こされる。これを軽減するためには、地震波場を測定するための補助的な地震計のアレイが必要であり、これによりNNを予測し、干渉計の出力から差し引くことが可能になる。しかし、センサーを設置するためのボーリング（掘削）を行うことは、多大なコスト要因となる。課題は、複数のセンサーを搭載したより少ない数のボーリング孔を利用する構成や、既存のインフラを活用する構成を検討することで、軽減性能を最大化しつつコストを最小限に抑えるよう、センサーの位置を最適化することにある。

手法
著者らは、ETのフルコーナー（4つの相関するミラー）に対して最適な地震計配置を決定するために、既報[12]で確立された最適化フレームワークを用いている。この最適化は、定常ウィーナーフィルタを用いた後のニュートンノイズ信号の平均残留値を最小化するように設計されている。

• 最適化アルゴリズム: プロセスでは、初期化のために粒子群最適化（PSO）を用い、その後の微調整のためにディファレンシャル最適化アルゴリズムであるAdamを使用する。
• 性能指標: 軽減係数 $M(f)$ は、4つのミラーのうちの最大残留値の平方根の逆数として定義され、振幅減少の保守的な指標として機能する。
• テストされたシナリオ:
  1. ベンチマーク: 1つのボーリング孔につき1つの地震計を配置するアレイで、総ボーリング孔数（$N = 20$ から $70$）を変化させる。
  2. マルチセンサー・ボーリング孔: 各ボーリング孔に $n = 1, 2, 3, 5, 10$ 個の地震計を含む $N = 20$ 個のボーリング孔のアレイ。最適化の制約として、ボーリング孔内の $x$ および $y$ 位置は変化させることができるが、$z$ 位置はボーリング孔内で独立して変化させる。
  3. インフラの統合: コーナーから5 kmまでの干渉計トンネル（アーム）内に、位置を固定した100個の地震計（100 m間隔）を追加する。これについては、特定の場所のさらなる最適化は行わない。

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. コスト効率の高いマルチセンサー構成: 本研究は、単一のボーリング孔内に複数の地震計を配置することで、はるかに多くのボーリング孔（各孔に1つのセンサーのみ）を持つアレイと同等の軽減能力が得られることを示している。これは、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までのニュートンノイズに対する効果的なブロードバンド軽減を提供する、と主張している。主な意義は、ETにとってコスト効率の高い戦略を特定したことにある。すなわち、多数のシングルセンサー・ボーリング孔を掘削するのではなく、複数のセンサーを備えたより少ない数のボーリング孔を使用し、既存のETインフラ内のセンサーで補完することを提案している。

著者らは、このような構成が「有望なブロードバンド性能」を提供し、「地震計の位置の変動に対して特に安定している」と結論付けている。これらの知見は、マルチセンサー・ボーリング孔の技術的実現可能性の検討や、ニュートンノイズ軽減の実行可能な戦略としてのETインフラへのセンサー統合に向けた研究を促進するものである。また、ブロードバンド最適化の方がより良いトレードオフをもたらす可能性があるものの、ここでテストされた単一周波数最適化戦略も有効であるとも述べている。今後の課題としては、現実的な地質における有限要素法シミュレーションや、融合アレイ（地震計、歪計、および傾斜計の組み合わせ）の研究が挙げられている。