地球物理学におけるMT法の基礎
地磁気地電流法（MT法）は、岩石の電気的特性、特に比抵抗（またはその逆数である導電率）に基づいて地球の地下構造を調査するために広く用いられている受動的な地球物理学的探査手法です。人工的なエネルギー源を必要とする能動的な地球物理学的探査手法（例えば、地震探査や直流比抵抗法）とは異なり、MT法は太陽風と磁気圏および電離圏との相互作用、ならびに大気中の雷活動から生じる自然の電磁場を利用します。MT法の主な利点は、数百メートルから数十キロメートル、あるいは数百キロメートルにも及ぶ深度まで探査できることであり、地熱、地殻変動、鉱物資源探査の研究において非常に重要な手法となっています。
地磁気地電流法の基本原理
概念的に言えば、MT法は、地球表面における電場（E）と磁場（H）の自然な変動を時間関数として測定する手法です。これらの電磁場の変動は地球内部に浸透し、地下の物質と相互作用します。岩石はそれぞれ鉱物組成、多孔性、流体含有量、温度、変質度などによって異なる抵抗率を持つため、地表で記録される応答には、地下の抵抗率分布に関する情報が含まれています。
周波数領域における電場と磁場の関係は、インピーダンステンソル（Z）によって表される。
E(ω) = Z(ω) · H(ω)
ここでωは角周波数です。インピーダンスは、地球が磁気信号を電気信号に「変換」する仕組みを表します。このインピーダンスから、見かけの抵抗率や位相などの重要なパラメータが導き出され、これらはMTにおける解釈の基礎となります。
電磁場の自然発生源
MT信号は一般的に、主に2つの周波数帯域から発生します。
1. 高周波（約1～10.000Hz）：地球規模の雷活動において支配的な周波数であり、電磁波（スフェリックス）を発生させる。この周波数帯は、浅層から中層までの調査に有用である。
2. 低周波（約0,0001～1Hz）：電離層および磁気圏における電流の変動（地磁気脈動）に由来する。低周波はより深く浸透できるため、地殻から上部マントルまでの構造をマッピングするのに適している。
信号の浸透深度は、周波数と媒体の抵抗率によって決まります。周波数が低いほど、信号はより深く浸透します。簡単に言うと、この概念は表皮効果として知られています。
表皮深さと調査深度の概念
表皮深さ（δ）は、電磁波の振幅が著しく減少する実効的な深さを表す指標です。おおよそ次のようになります。
δ ≈ 500 √(ρ / f)
ここで、δはメートル単位、ρは抵抗率（オームメートル）、fは周波数（Hz）です。この式は、抵抗率の高い岩石では信号がより深く浸透し、導電性の高い岩石では浸透が浅くなることを示しています。
例えば、抵抗率が100オームメートル、周波数が1Hzの場合、表皮深さは約5000メートルです。しかし、周波数が0,01Hzになると、深さは約50kmにまで増加します。これが、MT法が地域調査や大規模地熱システムにおいてしばしば選択される理由です。
現場でのMTデータ取得
MT測定は、観測地点（ステーション）に電場・磁場センサーを設置することによって行われます。データは、対象深度やその場所のノイズ品質に応じて、数時間から数日間にわたる時系列データとして記録されます。
測定対象となる構成要素は通常以下のとおりです。
– 電場（Ex、Ey）：地面に埋め込まれた2組の非分極性電極を使用して測定され、一定の長さ（例えば50～200m）でx方向とy方向に双極子を形成します。
– 磁場（Hx、Hy、Hz）：磁力計（通常、中高周波数では誘導コイル、低周波数ではフラックスゲート）を使用して測定します。
インピーダンステンソル解析においては、センサーの向き（南北方向と東西方向）を決定することが重要です。さらに、電極と土壌との接触状態、湿度条件、設置の安定性も、電界信号の品質に大きく影響します。
データ処理：時系列からインピーダンスまで
MTデータ処理の段階は一般的に以下の通りです。
1. 周波数領域への変換：フーリエ変換法を用いて時系列データを周波数スペクトルに変換します。
2. インピーダンス推定：各周波数帯域で安定したZテンソルを得るために統計的手法を用いて実施する。
3. ノイズフィルタリング：ノイズは、人間の活動（送電線、列車、工場）、センサーを振動させる風、または電極の接触不良から発生する可能性があります。
品質向上に重要な手法の一つに、リモートリファレンシングがあります。これは、対象エリアと電磁波の影響を受けにくい場所の2か所で同時にデータを記録する手法です。2つの場所間の相関関係を利用することで、局所的なノイズの影響を軽減し、より信頼性の高いインピーダンス推定値を得ることができます。
主要パラメータ：見かけ比抵抗と位相
インピーダンスから計算すると、次のようになります。
– 見かけ比抵抗（ρa）：特定の周波数の波が観測する「平均」比抵抗を表します。
– 位相（φ）：電界と磁界の間の位相差を示し、媒体の誘導特性に関連しています。
初期解釈では、周波数に対するρaとφの曲線が分析される。一般的な傾向としては、高周波数は浅い深度を、低周波数は深い深度を表す。曲線の急激な変化は、抵抗層／導電層の境界、変質帯、または流体の存在を示している可能性がある。
モデル寸法：1D、2D、3D
MT解析の結果は、地質構造の複雑さによって左右される。
– 1次元モデルは、比抵抗が深さ（水平層）のみで変化すると仮定しています。初期調査や単純な領域に適しています。
2Dモデルでは、比抵抗は横方向の一方向で深度とともに変化すると仮定し、他の方向は均質であると仮定します。地熱研究や細長い堆積盆地で広く用いられています。
– 3Dモデルは、あらゆる方向における比抵抗の変化を考慮します。これは複雑な地質構造を最も現実的に捉えることができますが、高密度のデータ、膨大な計算量、そして適切な逆解析手法が必要となります。
データから地下抵抗率モデルを導出するプロセスは、インバージョンと呼ばれます。MTインバージョンは一意ではなく、同じデータを複数のモデルで説明できる可能性があります。したがって、MT解釈は、地質情報、その他の地球物理学的データ（重力、磁気、地震など）、およびパラメータ制約（現実的な抵抗率の限界など）によって裏付けられる必要があります。
MT法の適用
MT法は、以下のような幅広い用途に利用されています。
1. 地熱探査：MT法は、粘土層（粘土の変質による導電性ゾーン）、上昇流経路、および比較的抵抗率の高いキャップ/貯留岩のマッピングに非常に効果的です。
2. 地殻構造と地殻の研究：下部地殻における主要な断層帯、縫合線、プレート境界、および導電性構造をマッピングする。
3. 鉱物探査：塊状硫化物や流体関連鉱化帯などの導電性物質の検出。
4. 堆積盆地および炭化水素盆地：比抵抗に影響を与える堆積層の厚さ、基盤岩、および岩相の変化を特定するのに役立ちます。
限界と課題
MTは強力ではあるものの、いくつかの限界がある。
– 文化的なノイズの影響を受けやすい：送電線、パイプ、フェンス、工業施設などが信号に干渉する可能性がある。
– 静的シフト：浅い不均質性（砂利、薄い導電層など）による電場の局所的な歪みで、見かけ比抵抗曲線の位相を大きく変化させることなくシフトする現象。これには、TDEM/CSAMTデータとの統合や特定の逆解析手法など、特別な処理が必要となる。
– 非一意性と解像度：MT法は、微細な幾何学的構造よりも比抵抗のコントラストに敏感です。解像度は深度が深くなるほど低下します。
閉鎖
地磁気地電流法（MT法）は、自然の電磁場を利用して、浅層から深層までの地下抵抗率分布をマッピングする受動的な地球物理学的探査手法です。地表における電場と磁場を測定し、それらをインピーダンス、見かけ抵抗率、位相に変換することで、MT法は地質構造、流体帯、変質帯、地殻境界に関する重要な知見を提供します。文化的ノイズやインバージョン結果の非一意性といった課題に直面しながらも、MT法は地熱探査、地殻構造研究、および様々な地殻調査において依然として主要な手法となっています。MT法の成功の鍵は、優れた調査設計、質の高いデータ取得、堅牢なデータ処理、そして地質学的理解やその他の補助データと統合された解釈にあります。
