地球物理学における地下構造モデリング手法
地下構造モデリングは、現代の地球物理学研究の根幹を成すものであり、資源探査、災害対策、インフラ整備計画におけるほぼすべての意思決定は、地表下の状況に関する我々の理解に基づいている。地下構造はほとんどの場合直接観測できないため、地球物理学者は地震波、重力場や磁場の変化、電気応答などの間接的な測定値を用いて地質構造を解釈する。これらのデータから、モデリングプロセスを用いて、2次元断面図または3次元ボリュームのいずれかの形で、妥当な物理的・地質学的表現を構築する。
一般的に、地球物理学における地下構造モデリングは、順モデリングと逆モデリング（またはインバージョン）という2つの主要なアプローチに分けられます。この2つは互いに補完し合う関係にあります。順モデリングは、地下構造モデルが既に仮定されている場合に地球物理学的データを予測するために用いられ、一方、インバージョンは測定データから地下構造モデルを推定することを目的としています。実際には、地球物理学者は、データに適合し、かつ地質学的に整合性のあるモデルが見つかるまで、この2つのアプローチを繰り返し用いることがよくあります。
1. モデリングの基本概念：データからモデルへ
各地球物理学的探査法は、それぞれ特定の物理量を測定します。地震探査では伝播時間または波形を、電気探査では比抵抗または導電率を、電磁探査では誘導応答を、重力探査では重力加速度を、磁気探査では磁場強度を測定します。これらの量は、岩石の特性（密度、磁化率、波速、空隙率、流体含有量、特定の鉱物含有量）と関連付けられます。
しかし、データと岩石特性の関係は一意であることは稀です。多くの異なるモデルが類似したデータ応答を生み出す可能性があります。これは非一意性と呼ばれ、地球物理モデリングにおける古典的な課題です。したがって、より的確な解を得るためには、地質、坑井データ、露頭観測、あるいは複数の手法の統合といった制約条件をモデル化する必要があります。
2. フォワードモデリング
順モデリングとは、想定される地下モデルの地球物理学的応答を計算することです。例えば、特定の密度コントラストを持つ層状モデルを作成すれば、地表における重力異常を予測できます。地震波速度モデルを指定すれば、地震波の到達時間を計算できます。順モデリングは、以下のような場合に重要です。
1. 地質学的仮説の検証（例：火成岩の貫入があるかどうか）。
2. 調査設計（経路間隔、音源周波数などの決定）。
3. 逆解析結果を検証する（物理モデルが現実的であることを確認する）。
計算上、順解析モデルは単純なもの（例えば、層状1次元モデル）から非常に複雑なもの（有限要素法／有限差分法に基づく3次元シミュレーション）まで多岐にわたる。岩石の不均質性、地形、異方性、および3次元効果によって複雑さが増す。
3. インバージョン：地下構造モデリングの中核
インバージョンは、観測データに最もよく一致する地下構造モデルを見つけることを目的としています。インバージョンは通常、最適化問題として解かれます。つまり、測定データとシミュレーションデータの差を最小化すると同時に、モデルが「暴走」したり非現実的になったりしないようにするための追加ルールが適用されます。この差はしばしばミスフィットと呼ばれ、安定化ルールは正則化と呼ばれます。
a. 決定論的逆変換
決定論的逆問題は、その効率性の高さから広く用いられています。例としては、最小二乗法、ガウス・ニュートン法、共役勾配法などが挙げられます。高速で大規模データセットにも適していますが、結果は初期仮定や正則化パラメータに依存します。
b. 確率的（ベイズ的）逆変換
ベイズアプローチでは、モデルを確率分布として扱います。その結果、単一のモデルではなく、複数の可能なモデルとその関連する不確実性の範囲が得られます。マルコフ連鎖モンテカルロ法（MCMC）などの手法は、不確実性の推定を重視する研究でよく用いられますが、計算コストが高いという欠点があります。
c. 正則化と地質学的制約
逆解析は不安定な場合が多いため、平滑化制約（モデルを滑らかにする）、減衰（極端な値を避ける）、疎性制約（明確な境界を促す）などの正則化手法が用いられます。さらに、地震探査データや坑井データから得られる地層境界などの地質学的制約を組み込むことで、モデルの現実性を高めることができます。
4. 地球物理データタイプに基づくモデリング手法
a. 地震モデリング
地震探査は、石油・ガス探査や地殻研究で広く用いられている高解像度手法です。地震モデリングには以下の要素が含まれます。
― 光線追跡法を用いて波の軌跡と伝播時間を計算する。
– 完全な波形をシミュレートするためのフル波形モデリング。
地震波マイグレーションとは、反射されたエネルギーを地表下の実際の位置へと伝達する過程のことである。
逆解析においては、地震トモグラフィー（走時による速度推定）やフルウェーブフォームインバージョン（FWI）が知られており、これらは非常に詳細な速度モデルを作成できるが、初期モデルやデータ品質に影響を受けやすい。
b. 地電気探査および電気抵抗率探査（ERT）によるモデリング
電気抵抗トモグラフィー（ERT）は、電流注入によって生じる電位差測定から地下抵抗率の分布をモデル化する手法です。この方法は、地下水調査、汚染調査、地盤工学において有効です。課題としては、電極接触効果、浅層部の不均質性、高い非線形性などが挙げられます。ERTの逆解析では、通常、正則化を用いてモデルを安定化させ、時間経過に伴う変化（汚染物質の移動や塩水侵入など）を監視するためのタイムラプス逆解析へと拡張することも可能です。
c. 電磁気モデリング（MT、TEM、CSEM）
地磁気地電流法（MT）などの電磁気的手法は、自然の電磁場を利用して深部までの導電率をマッピングするため、地熱や地殻変動の研究において重要です。一方、過渡電磁法（TEM）は、浅層から中層までの調査に適しています。電磁気モデリングは、マクスウェル方程式と強い3次元効果のため計算負荷が高くなる傾向があり、3次元モデリングと逆解析には相当な計算リソースと適切な正則化戦略が必要です。
d. 重力モデリング
重力データは密度変動をマッピングする。モデリング自体は比較的迅速だが、問題は非常に非一意的である。つまり、多くの異なる密度分布が同じ異常を引き起こす可能性がある。そのため、重力解析はほぼ常に地質学的制約、地震データ、または幾何学的境界に依存する。重力モデリングは、堆積盆地、岩塩ドーム構造、または貫入岩体のマッピングによく用いられる。
e. 磁気モデリング
磁気探査法は、磁化率と残留磁化に敏感です。磁気モデリングは、鉱物探査、構造マッピング、火成岩の識別などに有効です。課題としては、地域的な寄与と局所的な寄与の分離、日々の変動の補正、磁化方向の曖昧さなどが挙げられます。磁化率の3次元磁気インバージョンは非常に一般的で、岩石物理学的情報と組み合わせて用いられることもよくあります。
5. マルチメソッド統合と同時逆解析
各手法には感度の違いがあるため、統合は重要な戦略となる。例えば：
地震探査は、速度と音響インピーダンスに敏感である。
重力は密度に影響を受ける。
電磁気学は導電率に敏感である。
結合インバージョンは、複数のデータセットを組み合わせて一意性の欠如を減らし、信頼性を高める手法です。アプローチは2つあります。（1）結合目的関数を用いた同時インバージョン、または（2）ある手法で得られたモデルを別の手法で制約した逐次インバージョンです。地熱などの場合、MT、重力、地震探査を組み合わせた結合インバージョンによって、貯留層、キャップロック、変質帯のより一貫性のある全体像が得られることがあります。
6．主要な課題：解決、調査の深さ、および不確実性
地盤モデリングでは、常に以下の点を考慮する必要があります。
解像度：対象物の詳細度を示す指標。地震探査は通常解像度が高く、重力探査は解像度が低い。
– 調査深度：ERTは浅～中程度、MTは非常に深い調査が可能です。
– 不確実性：データノイズ、モデルの仮定、および一意性の欠如。
最良の方法は、最終モデルだけでなく、適合度、感度、不確実性の範囲も報告することです。これらがなければ、実際には多くの可能性のうちの1つに過ぎないモデルが「正確」であると誤解される危険性があります。
7. 最新の動向
コンピューティング、センサー、アルゴリズムの進歩により、より現実的なモデリングが可能になりつつあります。フル波形インバージョン、大規模3D電磁界インバージョン、およびジョイントインバージョンがますます一般的になっています。一方、機械学習ベースの手法は、解釈、パターン検出、および順方向モデリングを加速するためのサロゲートモデリングアプローチを加速するために使用され始めています。しかし、偏った解釈を避けるためには、機械学習は依然として物理原理と現場での検証によって調整される必要があります。
結論
地球物理学における地下構造モデリング手法は、高度なモデリング、インバージョン、および地質学的知識の統合を組み合わせたものです。地震探査、電気探査、電磁探査、重力探査、磁気探査といった各手法にはそれぞれ長所と短所があるため、手法の選択は対象、深度、および必要な解像度に合わせて行う必要があります。インバージョン、適切な正則化、および複数手法の統合を適切に適用することで、地下構造モデリングは資源探査、リスク軽減、および環境管理のための、より信頼性の高い画像を生成することができます。
ご希望であれば、この記事を特定の状況（例えば、地熱、石油・ガス、鉱物、地質工学探査など）に合わせて修正したり、参考文献を追加したり、引用文献を含むより学術的なバージョンを作成したりすることも可能です。
