AI要約
背景
表面波トモグラフィーは、地球の地殻や上部マントルの構造を研究する上で重要な役割を果たしてきました。従来の一般的な手法は、まず観測された走時データから2次元の位相速度（または群速度）マップを作成し、次に各地点の分散データのインバージョンにより1次元の横波（S波）速度プロファイルを求めるという、2段階のプロセスを踏みます。しかし、この手法や従来の直接反転法では、多くの場合、大圏コース（Great Circle）伝播を仮定しており、複雑な構造を持つ地殻浅部などでは、波の伝播経路の歪みが無視できず、推定結果にバイアスが生じる可能性がありました。一方で、波形全体を用いるフルウェーブフォームインバージョン（FWI）やアジョイントトモグラフィーは高精度ですが、特に高周波数帯域において計算コストが非常に高いという課題があります。手法
本論文で提案された手法は、中間的な速度マップを作成せず、分散データから直接3次元のS波速度構造を求める単一ステップのインバージョン法です。主な特徴は以下の通りです。Fast Marching Method( FMM) によるレイトレーシング: 各周期において、FMMを用いて震源・受信点間の走時と波の伝播経路を計算します。これにより、不均質メディアにおける大圏コースからの経路の逸脱を考慮することが可能になります。
反復的な更新: 反転プロセスにおいて、新しく得られた3次元速度モデルに基づいて、波の経路とデータ感度行列（カーネル）を逐次更新します。
ウェーブレット基底を用いた疎性制約付き反転: モデルパラメータをウェーブレット領域で表現します。ウェーブレット変換の多重解像度特性を利用し、データ制約が強い地域では高解像度を、制約が弱い地域では長波長の特徴を保持するようにします。
L1ノルム正則化: ウェーブレット係数がまばら（スパース）であると仮定し、L1 ノルム正則化を適用します。これにより、従来の滑らかさを強制する正則化よりも、データの分布に適合した適応的な解像度が得られます。
計算アルゴリズム: L1ノルムベースの最小化問題を解くために、反復再重合最小二乗法（Iteratively Reweighted Least Squares: IRLS）とLSQRアルゴリズムを組み合わせて使用しています。
物性間の関係: P波速度（α）と密度（ρ）は、地殻に関する経験的な関係式（Brocher, 2005）を用いてS波速度（β）に関連付けられており、これらによる分散への影響も感度行列に含まれています。
結果
台湾の台北盆地における、環境ノイズ相関法から抽出された短周期（0.5〜3.0秒）のレイリー波位相速度データにこの手法を適用しました。計算の収束: 反復的な反転により、走時残差の標準偏差は1.44秒から0.91秒に減少し、残差の平均はほぼゼロになりました。
速度構造の解明: 盆地内の低速度域と山地（大屯火山群など）の高速度域が明瞭に描き出されました。これらの結果は地表の地質構造と良好に一致しています。
従来手法との比較: 従来の2段階反転法の結果と概ね一致しましたが、データ網羅性が低い地域では、本手法の方がアーティファクトが少なく、より安定した結果が得られました。
解像度テスト: チェッカーボードテストにより、データ制約が良い地域では水平・垂直方向に高い解像度を持つことが確認されました。
考察
提案手法は、計算効率と精度のバランスに優れており、特に構造が複雑な地殻浅部の調査に有効です。利点: 2次元速度マップの構築が不要であり、3次元参照モデルの組み込みが容易です。また、ウェーブレット基底の採用により、データの密度に応じてモデルの細かさが自動的に調整されます。
限界と仮定: 現状では高周波近似（レイ理論）に基づいており、滑らかな不均質構造を想定しています。
今後の展望: 有限頻度効果（Finite-frequency effects）を取り入れた感度カーネルの導入や、実体波（P波・S波）の走時データとの共同反転（Joint Inversion）への拡張が期待されます。これにより、地殻浅部におけるP波・S波両方の速度構造をより正確に制約できる可能性があります。
こちらも、数式を追うのに時間がかかりました。パスを計算する点や、1次元のインバージョンを介さないという点が新鮮です。
0 件のコメント:
コメントを投稿