SAI-Krylovを用いた高精度かつ高速な3次元過渡電磁場順方向問題モデリングのためのモデル次数削減スペクトル要素法

より高速な仮想センサーで地中を覗く

埋蔵鉱床や地下水、隠れた断層の検出は、多くの場合、地球の電磁場に現れる微妙な揺らぎに依存します。これらの揺らぎを地中の像に変換するには、巨大な3次元モデル上で大量の計算を行う必要があり、非常に時間がかかることがあります。本論文は、実地探査に必要な精度を維持しつつ、計算時間とメモリ消費を大幅に削減する新しい電磁探査シミュレーション手法を紹介します。

電磁探査のシミュレーションが難しい理由

過渡電磁法（TEM）は、地表のループや接地ケーブルに電流パルスを流すことで動作します。電流が遮断されると、地下に渦状の電流が誘導され時間とともに減衰し、岩石の導電率や埋没鉱床に関する手がかりを運びます。受信器が測定する応答を予測するには、マクスウェル方程式を3次元格子上で解く必要があります。従来の手法は、各セル内の変化を線形近似で表す簡単な要素を用いることが多く、実装は容易ですが、微細な変化を捉えようとすると未知数の数が爆発的に増え、何千もの小さな時間刻みを踏むため計算コストが急増します。

少ない要素で精度を上げる格子

著者らは高次のスペクトル要素法を採用し、各格子セルに直線的な変化ではなくより柔軟な関数形状を持たせています。実際には、同じメッシュで電場の変化をより滑らかに表現できるか、あるいは同等の精度をはるかに少ないセルで達成できます。未知数をエッジにどう配置するか、各物理セルを基準立方体へどう写像するか、特殊なガウス–ロバット点と重みをどう用いるかを慎重に設計しています。重要な工夫として「縮約積分（reduced integration）」を導入し、特定の積分の厳密さをわずかに緩めることで導電率行列を厳密に対角化します。これによりスパース性が大きく向上し、メモリ使用量と線形系を解くコストが低下しつつ、目的とする次数での高い精度が保たれます。

物理を損なわずに数式を圧縮する

より高精度な格子を用いても、従来通り時間発展を行うと依然として高コストです。なぜなら時間刻みが変わるたびに巨大行列の再因数分解が必要になるからです。著者らは代わりに問題を連続時間系として書き直し、シフト・インバート型のKrylov部分空間アルゴリズムに基づくモデル次数削減技術を用います。平たく言えば、巨大な全系を電磁場が実際に変化する代表的なパターン群へ射影し、非常に小さな縮約系で表現する手法です。賢いスペクトル変換が重要なモードをクラスタリングし、アルゴリズムの収束を早めます。重要なのは、安価な残差指標を停止規準として用いることで、縮約空間の大きさを実行時に自動決定できる点です。その結果、解法は一度の主要な行列因数分解と少数の逆代入だけで済み、全ての中間時刻を逐次踏むことなく任意の時刻で電磁応答を評価できます。

手法の検証

精度と効率を評価するため、著者らは複数のベンチマークケースを数値実験しました。解析解が知られている単純な半空間では、異なる多項式次と積分戦略を標準的な逆時間ステッピング法と比較しています。2次以上では誤差が概ね1％以下で安定し、新しい縮約モデルでは平均誤差を約0.5％程度まで下げつつ、高次計算でおおむね16〜20倍の速度向上を達成し、メモリ増大を招きませんでした。強い抵抗率コントラストを持つ層状モデルでも、スペクトル要素とKrylovの組合せは全時間窓で約1％の誤差内に収まり、低次要素とより減衰的な時間積分を用いる従来の有限要素法より優れた性能を示しました。最後に、完全な3次元の硫化鉱体シナリオでは、接地線やループ源から誘起された場が複雑な導電体と相互作用して広がり、埋没ターゲットを高空間分解能で描く様子を追跡できました。

地下探査への意義

詳細な3次元電磁モデリングが重要な地球物理探査において、本研究は高忠実度と高速性を両立させる道を示します。高次スペクトル要素と注意深く設計されたモデル次数削減を組み合わせることで、実用的なTEM応答でサブパーセントの誤差を達成しつつ、実績のあるベースラインに比べて計算時間を一桁短縮できることを示しました。現場では、これにより生の観測設計から解釈可能な地下像へのターンアラウンドが速くなり、鉱床探査や地質リスク評価、地下水追跡などを、単純化された近似に頼ることなく物理に基づいた厳密なシミュレーションで実行しやすくなります。

引用: Fan, Y., Lu, K., Huang, Y. et al. Model-order-reduced spectral-element method for high-accuracy and fast 3-D transient electromagnetic forward modeling with SAI-Krylov. Sci Rep 16, 13356 (2026). https://doi.org/10.1038/s41598-026-44053-y

キーワード: 過渡電磁気学, スペクトル要素法, モデル次数削減, 地球物理探査, 数値シミュレーション