インバージョン 基礎の基礎

            138.873238   38.988033    -73.4
            138.873184   38.989834    -73.1
            138.873130   38.991635    -73.1
            138.873077   38.993436    -73.5
            

             255：olon=14700.0
             254：olon=14100.0
             253：olon=13500.0
             252：olon=12900.0
             251：olon=12300.0
           

            254←zone54だから基準経度は東経141度である。
            315.8  4317.6   -73.4
            315.8  4317.8   -73.1
            315.8  4318.0   -73.1
            315.8  4318.2   -73.5
            315.8  4318.4   -74.2
            315.8  4318.6   -75.1
           

• サンプルデータ１　サンプルデータ２
  
• プログラム
  
• Y軸の中心にあたる点のIGRFの値を計算しておく。サンプルデータの場合は143度15分、42度0分（1点であれば京都大学のサイトが便利）
  計算結果の例（とりあえず2000年で計算した）
  

               北緯 42.000 度   東経 143.250 度    標高: 0.0m での 2000 年 の地磁気要素 (IGRF-11)
               全磁力 (F)	 48821.3 nT
               偏角 (D)	 -8.599度
               伏角 (I)	 55.904度
               北向き (X)	 27060.4nT
               東向き (Y)	 -4091.9nT
               下向き (Z)	 40429.1nT
               水平 (H)	 27368.0nT
               

• 一般的に磁化したプリズムの組み合わせで磁気異常を説明する場合
  
  • 磁化強度を一定にして上面深度を変える
    
  • 上面深度を一定にして磁化強度を変える
    
  の2種類を考える。どちらがよりふさわしいかは、たとえばそのエリアの地質分布などを参考に考える。
  参考までに両者のモデルが一致することはほとんどない。
  
• モデルはいずれの場合も打ち切りによる影響を小さくするために両端には広めにプリズムを置く。
  
• 特に顕著な磁気異常の正負のピークがプリズムの両端に来るようにプリズムを置く。
  
• 推定するパラメータが多くなりすぎると計算が大変なので、モデルはアバウトに作ると良い。
  
  • ただし、あくまでも元データの分解能に対して判断するもの。
    一般に未知数のほうが観測データを大きく上回る場合が多い。このときに、モデルに様々な制約条件を加えていく。
    制約を加える代わりに解析モデルを粗いメッシュに切って未知数個数を減らそうとする方法は、元のデータに対しておかしな制約を加えたことに相当し、物理的な意味がわからなくなる恐れがある。
    （中塚正,2013：磁気探査入門講座より）
    
  • ここでは、こんな複雑な解析は扱わないのですが、基本姿勢として。
    求めるもの、無視すること、をきちんと把握して計算を始めることが大切。