[論文レビュー] L0+L1+L2 mixed optimization: a geometric approach to seismic imaging and inversion using concepts in topology and semigroup

August Lau, Chuan Yin|arXiv (Cornell University)|Jul 12, 2010

Topological and Geometric Data Analysis参考文献 7被引用数 4

ひとこと要約

本稿では、データの条件付けとミグレーション精度の向上を図るために、L1、L0、L2ノルムを逐次的に適用する幾何最適化フレームワークを提案する。L1はスパイク除去に、L0（ベッチ数を介して）はトポロジカルな複雑性の最小化に、L2の半群は幾何フリーなデータ適合を定義するために用いられ、最小限のアーティファクトで最適な速度における画像化を実現する。

ABSTRACT

The mathematical interpretation of L0, L1 and L2 is needed to understand how we should use these norms for optimization problems. The L0 norm is combinatorics which is counting certain properties of an object or an operator. This is the least amplitude dependent norm since it is counted regardless of the magnitude. The L1 norm could be interpreted as minimal geometric description. It is somewhat sensitive to amplitude information. In geophysics, it has been used to edit outliers like spikes in seismic data. This is a good application of L1 norm. The L2 norm could be interpreted as the numerically simplest solution to fitting data with a differential equation. It is very sensitive to amplitude information. Previous application includes least square migration. In this paper, we will show how to combine the usage of L0 and L1 and L2. We will not be optimizing the 3 norms simultaneously but will go from one norm to the next norm to optimize the data before the final migration.

研究の動機と目的

• トポロジカル、幾何的、最小二乗の原則を統合した、地震画像化のためのマルチノルム最適化フレームワークの構築。
• 従来のL2ベースのミグレーションが非地質的で正弦波的なアーティファクトを生じるという限界の解消。
• 最適なミグレーション速度選定の基準としてのトポロジカルな単純さ（L0）を用いた、速度解析と画像化の改善。
• L1（スパイク除去）、L0（トポロジカルな最小化）、L2（半群に基づく適合）を組み合わせることで、より頑健で地質学的に妥当な結果が得られることの実証。
• 合成および実データ2次元地震データセットを用いた手法の妥当性の検証により、トポロジカルな複雑性を最小化する速度値で画像化が向上することの確認。

提案手法

• 最初にL1変分ノルムを適用し、全 variation を最小化することで、スパイクのような特異イベントを除去し、平坦なイベントを保持するとともに、滑らかでない特徴を抑制する。
• ベッチ数を介してL0ノルムを用い、ミグレーテッド画像内のトポロジカルな複雑性（例：穴の数）を定量化し、その数を最小化することで最適なミグレーション条件を同定する。
• 拡散マップを通じてL2半群ノルムを用い、フーリエ変換やウェーブレット変換に依存せずに、内在するデータ幾何を定義することで、正弦波的アーティファクトを回避する。
• 最適化の順序を定める：L1でデータの条件付け、L0でトポロジカルな単純化、L2で最終的なミグレーションにより幾何的忠実性を確保する。
• 異なる速度で半群ミグレーションを実行し、画像の系列を生成した後、速度に応じてベッチ数（例：B1）を計算し、複雑性が最小になる点を特定する。
• Betti数のトレンドとミグレーション速度を比較することで手法の妥当性を検証し、最小値が正確な速度推定に対応することを確認する。

実験結果

リサーチクエスチョン

• RQ1 L1変分ノルムは、滑らかなイベントを歪めることなく、未ミグレーション地震データ内の誤ったスパイクや特異イベントを効果的に抑制できるか？
• RQ2 ベッチ数で測定されるL0によるトポロジカルな複雑性の最小化は、合成および実データセットにおいて、正しいミグレーション速度を信頼性高く示唆するか？
• RQ3 L2半群に基づくミグレーションは、標準的な最小二乗ミグレーションで一般的に見られる非地質的正弦波的アーティファクトを回避できるか？
• RQ4 L1、L0、L2ノルムを逐次的に適用することで、従来の手法と比較して、画像品質および速度解析がどの程度向上するか？
• RQ5 ミグレーテッド画像で観察される最小のベッチ数は、真の地下速度モデルを信頼性高く示す指標となるか？

主な発見

• L1変分ノルムは、図1(b)に示すように、実データからスパイク様の特異性を効果的に除去した一方で、平坦で滑らかなイベントを保持した。
• 3つのdiffractorsを含む合成データセットにおいて、ベッチ数B1は1500 m/sのミグレーション速度で最小値を示し、真の速度と一致した。これは、トポロジカルな単純さが有効な最適化基準であることを確認した。
• 実2次元スタック断面において、ベッチ数B1は約1800 m/sで最小値を示し、これが画像化に最適な速度であると示唆された（図7参照）。
• 半群ミグレーションの系列は、真の値から逸脱するにつれて、次第に複雑化する干渉パターンを生じさせた。これは、トポロジカルな複雑性が増加する傾向を裏付けた。
• L1（ノイズ抑制）、L0（トポロジカル最適化）、L2（幾何学的注意を払った適合）の組み合わせにより、標準的なL2手法よりも解釈可能で地質学的に妥当な画像が得られた。
• L2ミグレーションにおける半群の使用により、フーリエ変換などの事前定義された変換に依存せず、幾何学的注意を払った適合が可能となり、人工的な正弦波的アーティファクトが低減された。