(57)【要約】【課題】複屈折媒体の積層構造解析方法および積層構造解析装置，およびそのプログラムを記録したコンピュータ読み取り可能な記録媒体に関するものである。【解決手段】積層された複屈折媒体を解析する積層構造解析方法において，送信電磁波の電界と複屈折積層媒体の電磁波の伝搬に関するパラメータと複屈折媒体を伝搬する電磁波が反射する層における電磁波の散乱に関するパラメータと電磁波の偏向面の複屈折媒体の主軸に対する角度θをパラメータとして保持し，入力した電磁波の主軸に基づく座標成分について複屈折媒体における伝搬特性および反射面における散乱特性に基づいて複屈折媒体を伝搬し，反射面で反射し，複屈折媒体を伝搬してもどってくる電磁波を求め，前記角度θおよび反射面の深さをパラメータとして変更し，前記角度θおよび反射面の深さおよび受信電磁波の関係を出力する方法。【選択図】図１
2005069785.tif 000002
【特許請求の範囲】
【請求項１】
入力装置と出力装置とＣＰＵとメモリを備え，コンピュータ処理により，積層された複屈折媒体を解析する積層構造解析方法において，
送信電磁波の電界と複屈折積層媒体の電磁波の伝搬に関するパラメータと複屈折媒体を伝搬する電磁波が反射する層における電磁波の散乱に関するパラメータと電磁波の偏向面の複屈折媒体の主軸に対する角度θをパラメータとして保持し，
入力した電磁波の主軸に基づく座標成分について複屈折媒体における伝搬特性および反射面における散乱特性に基づいて複屈折媒体を伝搬し，反射面で反射し，複屈折媒体を伝搬してもどる電磁波を求め，前記角度θおよび反射面の深さをパラメータとして変更し，前記角度θおよび反射面の深さおよび受信電磁波の関係を出力することを特徴とする積層構造解析方法。
【請求項２】
入力電磁波の成分を要素とする送信波マトリックスの生成手段と，該電磁波の複屈折媒体中の伝搬特性を表す伝搬マトリックスを生成する手段と，性質の異なる媒体層の界面における散乱マトリックスを生成する手段と，性質の異なる界面で反射して該複屈折媒体を伝搬して受信される電磁波の成分を要素とするマトリックスを生成する手段とを備え，
送信電磁波の電界を( ＥxT, ＥyT) ，深さｒで反射して受信される受信電磁波の電界を( ＥxR，ＥyR )として，
各層の伝搬マトリックスをＴ，回転マトリックスをＲ（θ），反射面の散乱マトリックスをＳとし，
送信電磁波と受信電磁波の関係を
【数１】
2005069785.tif 000020
で求めることを特徴とする請求項１に記載の積層された積層構造解析方法。
【請求項３】
入力電磁波の成分を要素とする送信波マトリックスの生成手段と，該電磁波の複屈折媒体中の伝搬特性を表す伝搬マトリックスを生成する手段と，性質の異なる媒体層の界面における散乱マトリックスを生成する手段と，性質の異なる界面で反射して該複屈折媒体を伝搬して受信される電磁波の成分を要素とするマトリックスを生成する手段とを備え，
送信電磁波の電界を( ＥxT, ＥyT) ，深さｒの層の電界を( Ｅx ，Ｅy ) ，深さｒで散乱した伝搬する電界を( Ｅx ' ，Ｅ y' ），電磁波の電界を( ＥxR，ＥyR )として，
層ｉの伝搬マトリックスをＴi ，回転マトリックスをＲ（θi ），ｉ層での反射面の散乱マトリックスをＳi とし，
層ｉに入射される電界は層ｉ−１から出力される電界であるとし，反射するまでの深さｒの層ｉの電界を
【数２】
2005069785.tif 000021
で求め，
反射層での電界を
【数３】
2005069785.tif 000022
で求め、
反射した後の電界を
【数４】
2005069785.tif 000023
求めることを特徴とする請求項１に記載の積層構造解析方法。
【請求項４】
入力装置と出力装置とＣＰＵとメモリを備え，コンピュータ処理により，積層された複屈折媒体を解析する積層構造解析装置において，
送信電磁波の電界と複屈折積層媒体の電磁波の伝搬に関するパラメータと複屈折媒体を伝搬する電磁波が反射する層における電磁波の散乱に関するパラメータおよび電磁波の偏向面の複屈折媒体の主軸に対する角度θを保持するパラメータ保持手段と，
電磁波の主軸に基づく座標成分について複屈折媒体における伝搬特性および反射面における散乱特性に基づいて複屈折媒体を伝搬し，反射面で反射し，複屈折媒体を伝搬してもどる電磁波を求める伝搬特性演算手段と，演算結果を保持する演算保持手段とを備え，
前記角度θおよび反射面の深さをパラメータとして変更し，前記角度θおよび反射面の深さおよび受信電磁波の関係を出力することを特徴とする積層構造解析装置。
【請求項５】
入力電磁波の成分を要素とする送信波マトリックスの生成手段と，該電磁波の複屈折媒体中の伝搬特性を表す伝搬マトリックスを生成する手段と，性質の異なる媒体層の界面における散乱マトリックスを生成する手段と，性質の異なる界面で反射して該複屈折媒体を伝搬して受信される電磁波の成分を要素とするマトリックスを生成する手段とを備え，
送信電磁波の電界を( ＥxT, ＥyT) ，深さｒで反射して受信される受信電磁波の電界を( Ｅ xR,Ｅ yR ) として，
層ｉの伝搬マトリックスをＴi ，回転マトリックスをＲ（θi ），層ｉの反射面の散乱マトリックスをＳi とし，
【数５】
2005069785.tif 000024
を求めることを特徴とする請求項４に記載の積層構造解析装置。
【請求項６】
入力電磁波の成分を要素とする送信波マトリックスの生成手段と，該電磁波の複屈折媒体中の伝搬特性を表す伝搬マトリックスを生成する手段と，性質の異なる媒体層の界面における散乱マトリックスを生成する手段と，性質の異なる界面で反射して該複屈折媒体を伝搬して受信される電磁波の成分を要素とするマトリックスを生成する手段とを備え，
送信電磁波の電界を( ＥxT, ＥyT) ，深さｒの層の電界を( Ｅx ，Ｅy ) ，深さｒで散乱した伝搬する電界を( Ｅx ' ，Ｅy ' ），電磁波の電界を( Ｅ xR,Ｅ yR ) として，
各層ｉの伝搬マトリックスをＴi ，回転マトリックスをＲ（θ），反射面の散乱マトリックスをＳとし，
層ｉに入射される電界は層ｉ−１から出力される電界であるとし，反射するまでの深さｒの層ｉの電界を
【数６】
2005069785.tif 000025
で求め，
反射層での電界を
【数７】
2005069785.tif 000026
で求め、
反射した後の電界を
【数８】
2005069785.tif 000027
求めることを特徴とする請求項４に記載の積層構造解析装置。
【請求項７】
積層された複屈折媒体の積層構造解析するプログラムを記録するコンピュータ読み取り可能な記録媒体において，
該プログラムは，送信電磁波の電界と複屈折積層媒体の電磁波の伝搬に関するパラメータと複屈折媒体を伝搬する電磁波が反射する層における電磁波の散乱に関するパラメータと電磁波の偏向面の複屈折媒体の主軸に対する角度θをパラメータとして保持し，
入力した電磁波の主軸に基づく座標成分について複屈折媒体における伝搬特性および反射面における散乱特性に基づいて複屈折媒体を伝搬し，反射面で反射し，複屈折媒体を伝搬してもどる電磁波を求め，前記角度θおよび反射面の深さをパラメータとして変更し，前記角度θおよび反射面の深さおよび受信電磁波の関係を出力することを特徴とするコンピュータ読み取り可能な記録媒体。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は，複屈折媒体の積層構造解析方法および積層構造解析装置，およびそのプログラムを記録したコンピュータ読み取り可能な記録媒体に関するものである。特に，南極の氷床のレーダ観測等の地球物理的な積層構造の解析に有効なものである。
【背景技術】
【０００２】
１．はじめに
アイスレーダは，氷床や氷河の表面から基盤までの厚さの観測や氷体内部の物理構造を探る目的で使用されてきた。アイスレーダ技術は，送信アンテナから氷に向けて鉛直下向きに発射された電波が，氷体内部を伝搬し，反射や散乱を経て受信アンテナに帰還するものである。氷の内部の電波の伝搬にとって，最も重要な点の一つは，電波伝搬媒体が誘電体としての氷結晶であり，周波数に対応した複素誘電率の実数部と虚数部の値によって決定づけられることである。氷結晶は，結晶学的分類では六方晶構造の分子配列をもち，「氷Ｉｈ結晶」と呼ばれている。電波の速度や減衰をはじめ，反射・散乱・複屈折などが，氷Ｉｈ結晶の誘電特性に起因してひき起こされる。氷Ｉｈ結晶の誘電特性として特徴的な点は，結晶構造がＣ軸を対称軸とした一軸対称性をもつことに起因して，誘電率も一軸対称の異方性をもつ。氷Ｉｈ結晶が光の周波数帯で屈折率（誘電率の平方根）に一軸対称性をもつことによって，交差偏光板にかざした氷薄片がさまざまな色を見せることは良く知られている。この現象は複屈折現象によって発生している。いっぽう，氷河や氷床は多結晶の集合構造体である。これを構成する結晶粒は，氷河や氷床で直径約1mm から数cmの大きさであるが，この集合としての方向性（軸性）を示すことにより，結晶に比べて非常に大きな波長である，メガヘルツ帯の電波伝搬の際にも，同様な複屈折現象が発生する。
【０００３】
複屈折現象の観測には，偏波面を相互に平行或いは垂直に保った送・受信アンテナを，鉛直軸を回転軸として，水平面内で360 度回転するような観測を実施する。誘電異方性をもった氷床に入射する直線偏波の電波は，それぞれが直交する通常成分と異常成分に分かれてそれぞれが異なる位相速度をもちつつ伝搬する。反射散乱体では独立に反射散乱し，そして再びそれぞれが異なる位相速度をもちつつ氷床表面に現れる。空気中に現れた２つの成分は楕円偏波として再合成され，それを直線偏波アンテナで受信される。電波は複屈折媒体内を伝搬した結果として，アンテナ回転方位のπ/2を一周期として強弱をもった信号として現れる。
【０００４】
一方，氷床や氷河は，長い年月に降り積もった雪や凝結霜の積層体である。層を構成する結晶は，その堆積時の環境に応じて異なる結晶成長や変形履歴をもつ。したがって，上で述べた多結晶の誘電異方性を持つ氷床を構成する各層位について異なった性質を持っていることが予想される。実際，氷床の内部で発生する電波の反射散乱は，この層毎の多結晶氷としての複素誘電率の差違によって発生する。アイスレーダで観測される氷床内部からの反射エコーの３大原因は，（１）密度変化，（２）酸性度変化または，（３）結晶Ｃ軸方位分布の変化であることが知られている。特に（３）のメカニズムは氷床の中層部や深層部で卓越していることが報告されている。これまでの議論から，必然的に反射係数は方位異方性をももつことが知られている。
【０００５】
氷床内部での電波伝搬についても，上記で述べた反射・散乱と同様に，多結晶氷の誘電異方性を考慮するが重要である。氷床は積層構造をもつ複屈折媒体であり，その各層の誘電異方性はそれぞれに微少量だけ異なる。Hargreaves(1977)は，複屈折媒体としての氷床内部のレーダ波伝搬をはじめてモデル化した（非特許文献１参照）。
【０００６】
図３０は従来の氷床観測装置の構成を示す。図３０において，１２０は送信機である。１２１はパルス変調器である。１２２は電力増幅器１である。１２３は電力増幅器２である。１２４はミキサである。１２５は１７９ＭＨｚの発振器である。
【０００７】
１３０は受信器である。１３１は線型増幅器である。１３２は減衰器である。１３３は高周波増幅器である。１３４はミキサである。１３５はローパスフィルタである。１３６は発振器である。１３７は１ＭＨｚのフィルタ増幅器である。１３８は４ＭＨｚのフィルタ増幅器である。１３９は１４ＭＨｚのフィルタ増幅器である。１４３はＬＯＧ増幅器である。１４１はタイミングパルス発生器である。１４２はデジタルオッシロスコープである。１４０はデータ処理装置であって，受信データをデータ処理するものである。
【非特許文献１】Hargreaves, N D, 1977: The polarization of radio signals in the radio echo sounding of ice sheets, J. Phys. D: Appl. Phys, 10, 1285-1304。
【非特許文献２】Maeno, H., S. Fujita, K. Kamiyama, H. Motoyama, T. Furukawa, and S. Uratsuka, 1995: Relation  Between  Surface Ice Flow and Anisotropic Internal Radio-echoes in The East Queen Maud Land Ice Sheet, Antarctica, Proc. NIPR Symp. Polar Meteorol. Glaciol., 9, 76-86.
【発明の開示】
【発明が解決しようとする課題】
【０００８】
上記のHargreavesのモデルは，氷床の厚さ全体を一つの均一層として取り扱うという特徴をもっている。これは近似的な取り扱いには適していたが，「積層構造をもつ複屈折媒体」というより現実的なモデルにするには不十分である。しかし，積層構造の複屈折のモデル化は試みられたことがなかった。本願発明は複屈折率をもつ媒体の積層構造を解析する複屈折媒体の積層構造解析方法および構造解析装置，およびそのプログラムを記録したコンピュータ読み取り可能な記録媒体を提供することを目的にする。
【課題を解決するための手段】
【０００９】
本発明は，複屈折率をもつ積層媒体の解析には，各単位層の伝搬を記述し，なおかつその積層毎の伝搬の連鎖を示すには，後述する伝搬マトリックスを用いた定式化が有効であることに着目した。一方，近年になって，氷の誘電異方性の正確な値が決定され，積層構造を持ったモデル化が可能となった。これらの条件から，本発明者は，「積層構造をもつ複屈折媒体」の完全モデル化と，氷床内部の電波伝搬シミュレータを開発し，アイスレーダのデータに含有する各層位における氷床内部で発生する結晶方位分布と，氷結晶の統計的性質を知ることを可能にした。
【００１０】
本発明の複屈折率媒体の解析方法は，入力装置と出力装置とＣＰＵとメモリを備え，コンピュータ処理により，積層された複屈折媒体を解析する方法であって，送信電磁波の電界と複屈折積層媒体の電磁波の伝搬に関するパラメータと複屈折媒体を伝搬する電磁波が反射する層における電磁波の散乱に関するパラメータと電磁波の偏向面の複屈折媒体の主軸に対する角度θをパラメータとして保持し，入力した電磁波の主軸に基づく座標成分について複屈折媒体における伝搬特性および反射面における散乱特性に基づいて複屈折媒体を伝搬し，反射面で反射し，複屈折媒体を伝搬してもどってくる電磁波を求め，前記角度θおよび反射面の深さをパラメータとして変更し，前記角度θおよび反射面の深さおよび受信電磁波の関係を出力するようにした。
【００１１】
本発明の複屈折媒体の積層構造解析装置は，入力装置と出力装置とＣＰＵとメモリを備え，コンピュータ処理により，積層された複屈折媒体を解析するものであって，送信電磁波の電界と複屈折積層媒体の電磁波の伝搬に関するパラメータと複屈折媒体を伝搬する電磁波が反射する層における電磁波の散乱に関するパラメータおよび電磁波の偏向面の複屈折媒体の主軸に対する角度θを保持するパラメータ保持手段と，電磁波の主軸に基づく座標成分について複屈折媒体における伝搬特性および反射面における散乱特性に基づいて複屈折媒体を伝搬し，反射面で反射し，複屈折媒体を伝搬してもどる電磁波を求める伝搬特性演算手段と，演算結果を保持する演算保持手段とを備え，前記角度θおよび反射面の深さをパラメータとして変更し，前記角度θおよび反射面の深さおよび受信電磁波の関係を出力する。
【００１２】
本発明の積層された複屈折媒体の積層構造解析するプログラムを記録するコンピュータ読み取り可能な記録媒体は，該プログラムは，送信電磁波の電界と複屈折積層媒体の電磁波の伝搬に関するパラメータと複屈折媒体を伝搬する電磁波が反射する層における電磁波の散乱に関するパラメータと電磁波の偏向面の複屈折媒体の主軸に対する角度θをパラメータとして保持し，入力した電磁波の主軸に基づく座標成分について複屈折媒体における伝搬特性および反射面における散乱特性に基づいて複屈折媒体を伝搬し，反射面で反射し，複屈折媒体を伝搬してもどってくる電磁波を求め，前記角度θおよび反射面の深さをパラメータとして変更し，前記角度θおよび反射面の深さおよび受信電磁波の関係を出力することを特徴とする。
【発明の効果】
【００１３】
本発明によれば，複屈折率をもつ積層構造の解析を容易に行なうことができる。また，実際の観測結果との解析結果との比較も容易に行なうことが可能である。そのため，実際の観測結果に構造解析を容易に行なうことができ，複屈折媒体の積層構造を容易に把握することができる。
【００１４】
例えば，積層構造をレーダにより解析した場合に応用すると，（１）結晶主軸分布の平面内での異方性分布により生じる複屈折現象を確認できる。また，（２）通常波と異常波の位相が逆転するような積層構造を確認できる。さらに，（３）  （２）の現象は，電波の周波数および誘電異方性の割合に依存するので，これを解析することにより積層構造の誘電異方性の強度を解析できる。（４）さらに，一軸対称性の結晶構造をもつ積層構造の場合，散乱行列を考慮して求めた受信電力によりその異方性を確認できる。（５）電気伝導度の異方性が，偏波としての電波伝搬に与える影響は実用上小さいことが確認できる。
【００１５】
さらに，本発明によれば，複屈折率媒体の積層構造の解析結果と実際の観測データに基づく表示を比較することにより，解析に使用するパラメータの入力値を変更することにより解析結果と観測データを近づけ，入力パラメータを基に実際の積層構造を推定することが可能になる。
【発明を実施するための最良の形態】
【００１６】
入力装置により入力された送信電磁波の電界と複屈折積層媒体の電磁波の伝搬に関するパラメータと複屈折媒体を伝搬する電磁波が反射する層における電磁波の散乱に関するパラメータと電磁波の偏向面の複屈折媒体の主軸に対する角度を入力して保持してパラメータがパラメータ保持部に保持される。入力した電磁波の主軸に基づく座標成分について複屈折媒体における伝搬特性および反射面における散乱特性に基づいて複屈折媒体を伝搬し，反射面で反射し，複屈折媒体を伝搬してもどってくる電磁波がコンピュータ処理により求められ，表示出力される。
【実施例１】
【００１７】
図１（ａ），（ｂ），（ｃ），（ｄ）は，本発明の複屈折率をもつ積層体の構造解析方法の説明図である。
【００１８】
図１（ａ）において，１は複屈折媒体であって，層１〜層Ｎをもつものである。１１は送信波である電磁波を出力する送信アンテナであり，１２は複屈折率媒体の距離（深さ）ｒの面で反射し，積層体を伝搬した電磁波を受信する受信アンテナである。１３は基盤であって，複屈折媒体を支持する基盤である。
【００１９】
図１（ｂ）は，送信波の進行方向と偏波面の座標をあらわし，ｘｙ平面は複屈折媒体の表面の面であり，ｒ方向に送信波が進行する。送信波は楕円偏向しているが，電界のｙ成分は小さく，ｘ方向の成分が大きい。
【００２０】
図１（ｃ）は，送信波の座標系と複屈折媒体の主軸を表す座標系との関係を表す。図１（ｃ）において，ａは複屈折媒体１の主軸方向であり，ｂはその垂直方向である。
【００２１】
図１（ｄ）は，送信波の座標系と散乱体の主軸を表す座標系の関係を表す。ａ’は主軸であり，ｂ’はａ’に垂直な軸である。
２．複屈折媒体の積層構造体
2.1.伝搬モデル
以下，アイスレーダにより氷床（南極の氷床等）の層構造を解析する場合を例として説明する。一般にアイスレーダのアンテナのビーム幅は，１０数度以上あるため，例えば，氷床底面への散乱・伝搬を議論するには，斜め入射を考慮した２次元または３次元の取り扱いが必要である。しかし，水平積層氷床内部でのレーダ波の散乱・伝搬を議論する場合には，こうした効果は２次的であるため，ここでは，深さ方向（鉛直）に伝搬する１次元のモデルを構築する。鉛直方向のみの伝搬を考慮することから偏波面は，水平面内に限られるため複素誘電率（誘電率実数部および電気伝導度）の異方性も，水平面内のみに限定し，その原因としての氷結晶のＣ軸対称性を用いＣ軸平均方位が水平面内でガウス分布するという仮定を置いてモデル化する。
【００２２】
図１（ａ），（ｂ）のように，アイスレーダアンテナの向きを中心とした（ｘ−ｙ−ｒ）系と，氷床内の誘電率の方向性を元に（ａ−ｂ−ｒ) 系の２種類の座標系を設定する。r はその座標系にも共通で鉛直下向きとする。いま，氷床に入射されるアンテナからの送信電場を (ＥxT, ＥyT) とベクトルで表すことにする。
【００２３】
このベクトル成分は，ＥXT >ＥyTを想定し一方の偏波成分のみしか存在しない場合には，直線偏波の特性をもつアンテナの，偏波方向と垂直方向の電場成分とを考えることができるが，一般的には，（ＥX ，Ｅy ）をそれぞれ複素数として楕円偏波で扱う。
【００２４】
次に，深さｒの層において，図１（ｂ）に示す様に，氷床の誘電率が水平面内の２軸（ａ, ｂ方向と定義する）に主軸をもつような異方性をもつ伝搬路であるものとする。この時，アンテナを基準とした座標系（ｘ，ｙ）と，氷床の異方性誘電率の主軸による座標系（ａ, ｂ）にθだけ回転しているとする。前にも述べた様に，氷床内の多結晶氷は，複屈折性を持ち，その通常成分と異常成分では，伝搬定数の大きさが異なるため，（ａ, ｂ）成分には位相差が生じ，その結果，楕円偏波のパラメータが変形を受けることになる。伝搬路Ｔを伝搬してきた電場は，ある距離ｒにおいて伝搬路と同じ異方性の主軸による座標系（ａ',ｂ' ）にθだけ回転している散乱体Ｓで散乱される。散乱された電場は，同じ伝搬路を上向へと進むこととなる。以下に示すマトリックス計算では，距離ｒ毎に計算を行い積算して受信強度を得ている。
2.2.  マトリックス
氷床内での電波の伝搬を記述するために用いる電場ベクトルを以下の様にまとめて定義する。
( ＥxT，ＥyT) : 氷床直上にあるアンテナからの送信電波の強度（電場）
( Ｅx ，Ｅv ) :     深さｒにある氷床内で散乱する直前の電波の強度（電場）
( Ｅx ' ，Ｅy ' ） :  深さｒで散乱した直後の電波の強度（電場）
( Ｅ xR ，Ｅ yR ) :   氷床直上にあるアンテナで受信される電波の強度（電場）
次に，電場ベクトル成分の定義と，伝搬マトリックスＴを用いて，レーダから散乱体への電波の往復過程を説明する。模式的に以下のように表す。
【００２５】
往路：  ( ＥXT, ＥyT )→（Ｔ）→( ＥxT, ＥyT）→ (氷床内部散乱体での散乱) 
復路：  ( ＥXR，ＥyR )←（Ｔ）←( Ｅx ' ，Ｅy ' ）← (氷床内部散乱体での散乱) 
ただし，ここで使用する矢印は，レーダから電波が発射された後に電波に発生するイベントの順番を意味している。矢印の反転は，散乱体での散乱後に後方散乱する成分に着目することを意味している。いま，氷床表面から深さｒの位置まで，電波が伝搬するときの伝搬の様子を記述する量として伝搬マトリックス（Propagation Matrix）Ｔ↓を次のように定義する。
【００２６】
【数１】
2005069785.tif 000003
ここでｄaa, ｄbbは複屈折媒体内部での複屈折の主軸方向（電波の進行方向に直交する２軸）の成分であり，
【００２７】
【数２】
2005069785.tif 000004
と表す。ここで，k0は自由空間での伝搬定数であり，k a と kb は，ａ，ｂ偏波の伝搬定数, そしてｒはレーダからの距離である。k a と kb はさらに，伝搬経路での媒質の誘電率および電気伝導度の複素関数として表すことができる。
【００２８】
【数３】
2005069785.tif 000005
となる。ただし，ε0 とμ0 は真空中での誘電率と透磁率，εaa，εbb，σaa，σbbは，それぞれの偏波方向での，比誘電率と電気伝導度成分である。これらによって決定されるk a ,kb の実数部である kareal ,kbreal はそれぞれ位相回転をあらわし kaimaginary,kbimaginaryは減衰を表す。εa とεb とは，複屈折現象の記述にとって最も重要な部分である（誘電率はテンソルである）。氷多結晶の場合には，単位として考える氷の体積中を構成する結晶粒のもつ結晶方位とその誘電異方性をもとに，アンサンブル平均としてのテンソル量を求めることになる。すなわち，氷床氷のもつ結晶ファブリックから決まる。先にも述べた様に，座標系としてのａ軸，ｂ軸は，計算の便利のため，氷床の平面内での異方性の主軸に合わせる。
【００２９】
まず，送信アンテナでの電波発信から，散乱体へ入射するときの，偏波成分の変化は伝搬マトリックスを用いて次の様に定式化される。
【００３０】
【数４】
2005069785.tif 000006
ここでθ1 はアンテナ偏波面座標と氷床の誘電率テンソルの主軸座標のなす角であり，Ｒ( θ1 ）は座標回転角である。
【００３１】
【数５】
2005069785.tif 000007
氷床内部では一般には多層に積層し，かつ，その伝搬マトリックスＴは深さ方向に変化すると考えることができる。この場合には，深さ毎にＴ1,Ｔ2,Ｔ3,... ＴN といった複数のマトリックスに分解して，それぞれに固有の値を取ることで現実的なモデルの生成が可能である。この場合，第ｉ層での電界は次のように表すことができる。
【００３２】
【数６】
2005069785.tif 000008
この計算式は，第ｉ層は，第ｉ−１層から出力される電波を入力信号として前述の伝搬式により電波伝搬を計算することを表わす。回転角θは各層で同じであるとしてθ1 とした。
【００３３】
次に，入射に続く散乱体での散乱プロセスは，以下のように表現することができる。
【００３４】
【数７】
2005069785.tif 000009
Ｓは，以下に示す，散乱マトリックスである。θ2 はアンテナ偏波面と散乱マトリックスの主軸のなす角である。
【００３５】
【数８】
2005069785.tif 000010
散乱される深さからアンテナまでの電波の伝搬は，一般に下向きの伝搬と特性の異なるマトリックスＴ↑を用いて次のように表現できる。
【００３６】
【数９】
2005069785.tif 000011
反射した電磁波が多層構造を伝搬する式は次のように表すことができる。
【００３７】
【数１０】
2005069785.tif 000012
ここで，ｒは散乱面からの距離である。Ｔi 散乱面から表面に向かって数える。
ａ’，ｂ’は，先のアンテナの方位や氷床の誘電テンソルの主軸とは本来独立である。しかしながら，現実の氷床においては，氷床の誘電テンソルは，Ｃ軸分布の対称性をもった選択的異方性によって発生する。伝搬マトリックスも散乱マトリックスも，共にＣ軸分布の異方性によって生じるので，通常は平面内での主軸は共通であるとみなせる。また，散乱マトリックスが回転対称である場合，Ｓb ’a ’，Ｓa ’b ’は等しい。マトリックスの主軸を選択すれば，これらは共にゼロである。Ｓa ’a ’とＳb ’b ’の大きさに関しては，氷床内部での電波散乱に異方性がない場合には，Ｓa ’a ’とＳb ’b ’がほぼ等しいとみなすことができる。氷の内部反射の三大原因（密度揺らぎ，酸性度揺らぎ，結晶主軸方位揺らぎ）のうち，前者２つに起因する反射にはこれがあてはまる。
【００３８】
前述したように, θ2 はアンテナ偏波面と散乱マトリックスの主軸とのなす角であり, 氷床のような場合には, 通常θ1 とθ2 は一致する。例外的に一致しないケースとしては，結晶主軸分布が，何らかの側方からのねじりの力によって層位毎に微少量主軸方位を変えてしまうようなケースを考えることができる。
【００３９】
本来，多層構造で揺らぎながら徐々に変動する伝搬マトリックスを仮定する場合，実際の電波伝搬において「揺らぎ」は，電波の散乱を発生するが，位相の急激な変化等を発生することはないので，本実施例１で伝搬マトリックスＴを一定において近似する。逆に，複数のマトリックスＴi 用いて多層化する数層程度のモデル計算でも，伝搬モデルの各要素がどのようなふるまいをするかを実証することが優先される。こうした多層の場合で，それぞれの層で異なる伝搬マトリックスをとる場合にもそれぞれの層での合成となるので，本実施例１では現実の氷床の記述に最も基本的なモデルとして，Ｔが深さ方向に変化しない単層の場合について説明する。
【００４０】
次に散乱される深さからアンテナまでの電波の伝搬は，一般に下向きの伝搬と特性の異なるマトリックスＴ↑を用いて以下のように表現できる。
【００４１】
一方，これまで述べてきた様な複屈折の原因を考慮すると，上向きと下向きとは双方性(reciprocity) が成立し，Ｔ＝Ｔ↓＝Ｔ↑，かつ，θ＝θ1 ＝θ2 として考えても，大きな矛盾を持たない。このとき，
【００４２】
【数１１】
2005069785.tif 000013
と合成することにより，送信ベクトルから受信ベクトルを求める定式化が可能となる。
2.3 現実に氷床・氷河内部で発生しうる複素誘電率テンソルの範囲の評価
2.3.1 散乱体の異方性
氷Ｉｈ結晶が光の周波数帯で屈折率に一軸対称性をもつことによって，交差偏光板にかざした氷薄片がさまざまな色を見せる事は前述した。これは，氷床や氷河に通常発生する応力形状（圧縮，引っ張り，純粋剪断，単純剪断）に対しても，発生する結晶Ｃ軸の選択的配向もやはり「一軸対称性」をもつと言える。こういった物理的な応力が発生すると散乱体の構造自体も異方性を持つことが推定できる。散乱体である結晶構造集合体が電波の波長に対して十分な体積散乱を起こす構造であるとすると，散乱体の構造が他の異方性を生じさせる要素にくらべ電波伝搬にどの程度反映するかモデル計算を行う必要性がある。本実施例１では，散乱体の異方性の割合を決めるために，基盤または氷床の内部からの反射強度の異方性が約6 dBから10 dB であったとの報告があるので，計算にはこの値を採用した。
2.3.2 偏波成分の割合
一般的にモデル上では任意の楕円偏波を取り扱う事が可能であるが，ここでは最も単純な場合として，ＥxT＞Ｅ yT とする。一般的なアンテナの交差偏波成分の比は，−３０dB以上とされている。アイスレーダに用いられる八木アンテナは，この条件を満たすが，より現実的なパラメータとしてアンテナ自体が完全に１方向の偏波成分で無いことやアンテナ近傍に雪上車の筐体等があることを考慮してＥxTに対してＥyTに−２５ｄＢ偏波成分があるものとした。
2.3.3 誘電率実数部の異方性とそれにより生じる複屈折強度
誘電体中での電波の伝搬速度Ｖe は，誘電率実数部ε' の関数として次式のように近似できる。
【００４３】
【数１２】
2005069785.tif 000014
ここで，ｃは，光の速度である。すなわち，誘電率実数部に異方性をもつ氷床が存在するとすれば伝搬速度の違いによる電界強度の変化が予想される。すなわち，複屈折を生じる結果となることが予想される。実施例１では単層モデルとして南極氷床の鉛直方向（深さ方向）の誘電率実数部の異方性を仮定し，θ方向にはねじれを持たず一定の方位角を保っているとする。現実の氷床の中で発生する複屈折の強度は，水平面内の誘電率テンソル成分の差によって以下のように表現できる。
【００４４】
【数１３】
2005069785.tif 000015
ここで，φは，正常成分と異常成分の位相差である。ｆは電波の周波数，c は真空中の電波の速度である。この式は，任意の２つの深度ｚ１とｚ２間の電波の往復伝搬の間に，位相差がどれだけ出現するかを表した式である。この式により，レーダに使用する電波の周波数が高いほど，φは大きく現れることがわかる。具体例として，上の式に典型的な数字を当てはめてみると，φの具体的範囲がわかる。例えば，氷の誘電率を３．１５とし，誘電異方性を０．０３４という単結晶の値とする。観測周波数は，日本の観測隊が使用をしているＶＨＦ帯の１７９ＭＨｚとすると，１，０００ｍの区間の電磁波の往復に対してのφは，２２．８πとなる。実際の氷床の内部での誘電異方性は，多結晶内部の伝搬であることから，φはこの１／３程度以下になる。また，もし，レーダの周波数を変更する場合には，φと周波数の比例関係を考慮して見積もる必要がある。
2.3.4 電気伝導度の異方性
一般的に誘電体の電気伝導度は誘電率虚数部と相関があり，減衰係数を示す要素として扱われる。また，減衰係数には，温度依存性がある。しかし，これまでの氷床の誘電率虚数部の研究の中で氷床の水平面内に電気伝導度の異方性が存在することは，証明されていない。加えて，氷がそのようなふるまいをするかどうかは，物理的に証明されていない。しかし，本発明の実施例１のモデルは，電気伝導度の異方性を排除することなく計算可能であるので，本実施例１では，電気伝導度の異方性が２０％までの異方性を考慮してどの程度の受信強度に異方性を生じるかについての検討も行った。  
【００４５】
図２は本発明のシステム構成の実施例１を示す図である。図２において，２１はＣＰＵ，２２はメモリ，２３は各種演算手段であって，メモリにロードされたプログラムを表す。２４は入出力インタフェースである。２５は入力装置であって，キーボード等である。２６はディスプレイである。２７はプリンタである。２８は外部記録媒体入力装置であって，磁気ディスク，磁気テープの磁気記憶媒体，ＣＤＲＯＭ，ＤＶＤ等の光記憶媒体等の記録データを入力するための装置である。３０は観測データ記録媒体であって，レーダ装置等により実際に観測した複屈折媒体の観測データを記録するものであり，磁気記録媒体，ＣＤＲＯＭ等の光記録媒体である。
【００４６】
各種演算手段２３において，４１は送信波マトリックス生成手段であり，送信波の電界成分のマトリックスを生成するものである。４２は伝搬マトリックス生成手段であって送信波および反射波が複屈折媒体中を伝搬する状態を表す伝搬マトリックスを生成する手段である。４３は散乱マトリックス生成手段であって，電磁波が反射される面の散乱マトリックスを生成するものである。４４は回転マトリックス生成手段であって，送信波のＸＹ座標軸と複屈折誘電体表面の主軸の座標軸（ａ，ｂ）との座標軸の作る角度に基づく座標回転を表す回転マトリックスを生成するものである。４５は受信波マトリックス生成手段である。４６は伝搬特性演算手段であって，送信マトリックス，回転マトリックス，散乱マトリックス，伝搬マトリックス等に基づいて複屈折媒体中を伝搬して，受信される電磁波の特性を演算する手段である。
【００４７】
４７はその他の演算手段であって，各種図表を生成するための演算手段等である。４８は複屈折率強度演算手段である。
【００４８】
図３は本発明の送信波マトリックス生成，伝搬マトリックス生成，散乱マトリックス生成，回転マトリックス生成，受信波マトリックス生成の構成を示す図である。図３において，２５は入力装置であって，各種パラメータ（送信波（ＥＴ），送信角周波数（ω）等）を入力するものである。５０は入力パラメータ保持部であって，入力装置２５から入力されたパラメータを保持するものである（メモリ上の領域）。４１は送信波マトリックス生成手段であって，送信波のパラメータに基づいて送信波のマトリックス〔ＥＴ〕＝〔ＥxT，ＥyT〕（二行一列の行列であるが，便宜的にこのように表す）の配列を生成するものである。４２は伝搬マトリックス生成手段であって，伝搬マトリックスの要素（ｄaa，０，０，ｄbb）の配列を生成するものである。４３は散乱マトリックス生成手段であって，入力パラメータに基づいて散乱マトリックスＳ＝（Ｓaa，Ｓab，Ｓba，Ｓbb）の配列を生成するものである（散乱マトリックスは二行二列の行列であるが，便宜的にこのように表す）。４４は回転マトリックス生成手段であって，ＸＹ座標軸とａｂ座標軸のなす角θに基づく回転マトリックスＲ（θ）＝（ｃｏｓθ，−ｓｉｎθ，ｓｉｎθ，ｃｏｓθ），回転マトリックスＲ（−θ）＝（ｃｏｓθ，ｓｉｎθ，−ｓｉｎθ，ｃｏｓθ）を生成するものである（Ｒ（θ），Ｒ（−θ）はいずれも二行二列のマトリックスであるが，便宜的にこのように表す）。４５は受信波マトリックス生成手段であって，演算で求められた受信波の成分（ＥxR，ＥyR）を要素とするマトリックスの配列を生成するものである。
【００４９】
５１は送信波マトリックス保持部であって，送信波マトリックス生成手段４１の生成した送信マトリックスを保持するものである。５２は伝搬マトリックス保持部であって，伝搬マトリックス生成手段４２の保持した伝搬マトリックスを保持するものである。５３は散乱マトリックス保持部であって，散乱マトリックス生成手段４３の生成した散乱マトリックスを保持するものである。５４は回転マトリックス保持部であって，回転マトリックス生成手段４４の生成した回転マトリックスを保持するものである。５５は受信波マトリックス保持手段であって，受信波マトリックス生成手段４５の生成した受信波マトリックスを保持するものである。
【００５０】
図４は本発明の伝搬マトリックス要素生成手段の構成である。６０は伝搬マトリックス要素生成手段である。６１はｋ演算部であって，伝搬マトリックスの生成に必要な電磁波のａ方向の偏向成分，ｂ方向の偏向成分の伝搬定数を演算するものである。５０は入力パラメータ保持部である。６２はｋa 演算部であって，入力パラメータ（ε0 ，μ0 ，εa ，σa ，ω）に基づいてｋa を演算するものである。６３はｋb 演算部であって，入力パラメータ（ε0 ，μ0 ，εb ，σb ，ω）に基づいてｋb を演算するものである。
【００５１】
６５は位相回転保持部であって，ｋa 演算部６２とｋb 演算部６３の演算したｋa とｋb の実部である（ｋareal ，ｋbreal ）を保持するものである。ｋa 演算部６２とｋb 演算部６３の演算したｋa とｋb の虚部である（ｋaimaginary，ｋbimaginary）を保持するものである。４２は伝搬マトリックス生成手段であって，伝搬マトリックス（ｄaa，０，０，ｄbb）を演算するものである（伝搬マトリックスは二行二列の行列であるが，便宜的にこのように表す）。５２は伝搬マトリックス保持部であって，伝搬マトリックス生成手段４２の求めた伝搬マトリックス（ｄaa，０，０，ｄbb）を保持するものである。６５は位相回転保持部であって，ｋ演算部６１の生成したｋa とｋb の実数部（ｋareal ，ｋbreal ）を保持するものである。６６は減衰率保持部であって，ｋ演算部６１の生成したｋa とｋb の虚数部（ｋaimaginary，ｋbimaginary）を保持するものである。
【００５２】
図４において，ｋa 演算部６２は，入力パラメータ保持部の保持するパラメータ（ε0 ，μ0 ，εa ，σa ，ω）に基づいて，
【００５３】
【数１４】
2005069785.tif 000016
によりｋa を演算し，ｋa の実部および虚部を求める。また，ｋb 演算部６３は，入力パラメータ保持部の保持する入力パラメータ（ε0 ，μ0 ，εb ，σb ，ω）に基づいて，
【００５４】
【数１５】
2005069785.tif 000017
を演算し，ｋb の実部および虚部を求める。求められたｋa の実部ｋareal とｋb の実部ｋbreal は位相回転保持部６５に保持される。また，求められたｋa の虚部ｋaimaginaryとｋb の虚部ｋbimaginaryは減衰率保持部６６に保持される。
【００５５】
伝搬マトリックス生成手段４２は位相回転保持部６５の保持する（ｋareal ，ｋbreal ）と減衰率保持部６６の保持する（ｋaimaginary，ｋbimaginary）を入力し，さらにｄa ＝ｅｘｐ（−ｊｋ0 ｒ＋ｊｋa ｒ）とｄb ＝ｅｘｐ（−ｊｋ0 ｒ＋ｊｋb ｒ）およびｋa ＝ｋareal ＋ｊｋaimaginary，およびｋb ＝ｋareal ＋ｊｋbimaginaryに基づいて伝搬マトリックス（ｄaa，０，０，ｄbb）を求める。伝搬マトリックス（ｄaa，０，０，ｄbb）は伝搬マトリックス保持部５２に保持される。
【００５６】
図５は本発明の実施例の伝搬特性演算するための構成を示す。図５において，２５は入力装置であって，演算に必要な各パラメータに設定値（シミュレーションするための具体的な数値）を入力するものである。４６は伝搬特性演算手段である。５０は入力パラメータ保持部であって，各パラメータに設定された設定値を保持するものである。４１は送信波マトリックス生成手段である。４２は伝搬マトリックス生成手段である。４３は散乱マトリックス生成手段である。４４は回転マトリックス生成手段である。４５は受信波マトリックス生成手段である。５１は送信波マトリックス保持部である。５２は伝搬マトリックス保持部である。５３は散乱マトリックス保持部である。５４は回転マトリックス保持部である。５５は受信波マトリックス保持部である。
【００５７】
７１は演算結果保持部であって，伝搬特性演算手段４６の演算結果を保持するものである。７３は受信波出力換算手段であって，マトリックス成分（ＥxR，ＥyR）受信電力に換算するものである。８０は各種図表生成手段であって，演算結果をもとに表示する各種図表を生成するものである。９１は複屈折強度演算部であって，前述の式（１１）を演算するものである。
【００５８】
図５の構成の動作を説明する。入力装置２５により構造解析に必要な各パラメータに設定する数値を入力する。各パラメータの数値は入力パラメータ保持部に設定される。送信波マトリックス生成手段４１は，送信波のマトリックスの要素ＥxTとＥyTに入力された値を設定する。ｋ演算部６１において入力パラメータ（ε0 ，μ0 ，εa ，σa ，ω）の入力値に基づいてｋa を演算し，ｋb 演算部において入力パラメータ（ε0 ，μ0 ，εb ，σb ，ω）の設定値に基づいてｋb を演算する。ｋa とｋb の実部は位相回転保持部６５に保持され，虚部は減衰率保持部６６に保持される。伝搬マトリックス生成手段４２は位相回転保持部６５に保持されている位相回転と減衰率保持部６６に保持されている減衰率に基づいて伝搬マトリックスの要素（ｄaa，０，０，ｄbb）の値を演算し，伝搬マトリックスを生成する。散乱マトリックス生成手段は散乱マトリックスの要素（Ｓaa，Ｓab，Ｓba，Ｓbb）に入力された設定値を設定し，散乱マトリックスを生成し，散乱マトリックス保持部４３に保持する。回転マトリックス生成部は回転各θの入力値に基づいて，ｓｉｎθ，ｃｏｓ（θ）を演算し，回転マトリックスＲ（θ）とＲ（−θ）を演算し，回転マトリックス要素に設定する。その回転マトリックスを回転マトリックス保持部５４に保持される。
【００５９】
伝搬特性演算手段４６は，送信波マトリックス保持部５１，伝搬マトリックス保持部５２，散乱マトリックス保持部５３および回転マトリックス保持部５４に保持されているそれぞれのマトリックスを入力し，
【００６０】
【数１６】
2005069785.tif 000018
を演算して受信波の電界成分ＥxR，ＥyRを求め，演算結果保持部７１に保持する。受信波マトリックス生成手段４５は電界成分ＥxR，ＥyRを要素とする受信波マトリックスを生成し受信波マトリックス保持部５５に保持する。
【００６１】
受信波出力換算手段７３は演算結果保持部７１に保持されている演算結果の受信電界をもとに，
ＰxR＝２０ｌｏｇＥxR      （ｄＢ）                (12)
ＰyR＝２０ｌｏｇＥyR      （ｄＢ）                (13)
を演算し，受信波マトリックス生成手段４５はその換算された電力に基づいた受信波マトリックスを生成し，演算結果保持部７１に保持する。また，必要ならば，複屈折強度演算部９１は，
【００６２】
【数１７】
2005069785.tif 000019
に従って複屈折強度φを求め，演算結果保持部７１に保持する。各種図表生成手段８０は出力する形式のグラフを作成し，演算結果出力部８２により出力する。
【００６３】
受信波出力換算手段７３は受信波の各要素の電界を電力に換算する。各種図表生成手段８０は受信電力に関する図表を作成し，演算結果出力部８２により作成した図表を表示出力する。
【００６４】
図６は本発明の入力パラメータ保持部の構成を示す。メモリ上に保持される入力パラメータを示す。図６は入力パラメータ群Ａ，入力パラメータ群Ｂ，入力パラメータ群Ｃが保持されている場合を示している。
【００６５】
図７は本発明の演算結果保持部の構成を示す。図７（ａ）は受信電界保持部の構成を示す。受信電界は入力パラメータ群と角度θ0 ，θ1 ，θ2 ，・・・θn 毎に深さｒ1 ，ｒ2 ，・・・ｒn 毎の式（９）の演算結果を保持する。図７（ａ）ではパラメータ群Ａのみに対する受信電界を示しているが，必要に応じて，他のパラメータ群の演算結果を保持することもできる。
【００６６】
図７（ｂ）は正常成分と異常成分の位相差φについての演算式（１１）の演算結果を保持する位相差保持部を示している。パラメータ群毎に深さの差毎にφを保持する（図７（ｂ）では深さｒ１から深さｒ２での差を示している）。
【００６７】
図８（ａ）は本発明の送信波マトリックス生成手段のフローチャートである。送信波マトリックスの生成は，まず送信波の電界パラメータ（ＥxT，ＥyT）を入力する（Ｓ１）。その電界パラメータにおける電界の成分（ＥxT，ＥyT）を要素とする送信波マトリックスの配列を生成する（Ｓ２）。生成した送信波マトリックスを保持する（Ｓ３）。
【００６８】
図８（ｂ）は散乱マトリックス生成手段のフローチャートである。散乱マトリックスのパラメータ（Ｓaa，Ｓab，Ｓba，Ｓbb）を入力する（Ｓ１）。散乱マトリックス（Ｓaa，Ｓab，Ｓba，Ｓbb）の要素の配列を生成する（Ｓ２）。生成した散乱マトリックスを保持する（Ｓ３）。
【００６９】
図９（ａ）は回転マトリックス生成手段のフローチャートである。回転マトリックスのパラメータθを入力する（Ｓ１）。回転マトリックスの要素（ｃｏｓθ，−ｓｉｎθ，ｓｉｎθ，ｃｏｓθ）を要素とする配列を生成する（Ｓ２）。生成した回転マトリックスＲ（θ，Ｒ（−θ）を保持する（Ｓ３）。
【００７０】
図９（ｂ）はｋ演算部のフローチャートである。パラメータ（ε0 ，μ0 ，εa ，σa ，ω）の入力する（Ｓ１）。式（３ａ）によりｋareal ，ｋaimaginaryを求め，位相回転保持部，減衰率保持部に保持する（Ｓ２）。パラメータ（ε0 ，μ0 ，εb ，σb ，ω）を入力する（Ｓ３）。式（３ｂ）を演算し，てｋbreal ，ｋbimaginaryを求め，位相回転保持部，減衰率保持部に保持する（Ｓ４）。
【００７１】
図１０は伝搬マトリックス保持部のフローチャートである。伝搬マトリックスの要素（ｄaa, ０，０，ｄbb）の配列を生成する（Ｓ１）。ｋa （ｋareal ，ｋaimaginary）およびｋa （ｋbreal ，ｋbimaginary）を位相回転保持部および減衰率保持部から入力して保持する（Ｓ２）。式（２ａ），式（２ｂ）により伝搬マトリックス要素ｄaaとｄbbを求める（Ｓ３）。生成した伝搬マトリックス（ｄaa，０，０，ｄbb）を保持する（Ｓ４）。
【００７２】
図１１は，本発明の伝搬特性演算手段のフローチャートＡである。入力装置により入力された各パラメータの設定値を入力パラメータ保持部に保持し，電界成分の値により，送信波の電界マトリックスの各要素の配列に電界値を設定する（Ｓ１）。入力された散乱マトリックスの各要素の設定値を，散乱マトリックス要素の配列に設定する（Ｓ２）。角度θに従って，ｃｏｓθ，ｓｉｎθを求め，回転マトリックスの各要素の配列に設定する（Ｓ３）。ｋ演算部においてｋa 演算部は，パラメータ（ε0 ，μ0 ，εa ，σa ，ω）の設定値を入力し，ｋa を演算する（Ｓ４）。さらに，計算されたｋareal ，ｋaimaginaryを位相回転保持部，減衰率保持部に保持する（Ｓ５）。また，ｋb 演算部は，パラメータ（ε0 ，μ0 ，εb ，σb ，ω）の設定値を入力し，ｋb を計算する（Ｓ６）。計算されたｋbreal ，ｋbimaginaryを位相回転保持部，減衰率保持部に保持する（Ｓ７）。伝搬特性演算部は位相回転保持部に保持されているｋareal ，ｋbreal ，減衰率保持部ち保持されているｋaimaginary，ｋbimaginaryを入力し，伝搬マトリックスの要素ｄaaとｄbbを求め，伝搬マトリックス（ｄaa，０，０，ｄbb）の配列に設定する（Ｓ８）。各要素に設定された各マトリックス要素（ｄaaとｄbb）に基づいて式（９）により受信波のｘ座標成分とｙ座標成分を求め，受信波マトリックスを生成し，保持する（Ｓ９）。
【００７３】
図１２は本発明の伝搬特性演算手段のフローチャートＢであって，図１１の伝搬特性演算手段のフローチャートＡのステップＳ９の動作の詳細フローチャートである。図１２は深さｒi を変更することによりそれぞれの深さで反射する受信電界を求めるフローチャートである。
【００７４】
図１２のＳ１において，送信電界（ＥxT，ＥyT），角度（θ），深さ（ｒ）の初期設定をする。Ｓ２で座標回転角度θi を定める。Ｓ３，Ｓ４で入力電界（ＥxTi ＥyTi ），測定する反射面の深さｒi を定める。Ｓ５においてその条件で伝搬式（９）を演算する。Ｓ６で演算結果を保持する。Ｓ７で伝搬特性の演算処理を終了するか判断して，終了しないなら，送信電界、深さ、角度を入力する。その際それぞれの値を変更するなら，変更した値を入力し，変更しないなら前回の値より，Ｓ２以降の処理を繰り返す。入力値に基づき深さｒi+1 で反射した受信電界を求める。Ｓ７で処理を終了するなら伝搬特性演算の処理を終了する。
【００７５】
図１３は本発明の実施例２であって，複屈折媒体の伝搬マトリックスが層毎に異なるとして計算する場合の構成を示す。５０は入力パラメータ保持部で，層毎に誘電率，電気伝導度等のパラメータを保持する。４２は伝搬マトリックス生成手段である。５２は伝搬マトリックス保持部であって，層毎に伝搬マトリックスを保持する。５２１は層１の伝搬マトリックス，５２２は層２の伝搬マトリックス，５２３は層Ｎの伝搬マトリックスである。
【００７６】
入力パラメータ保持部５０に保持されている層１パラメータにより伝搬マトリックス生成手段４２は層１の伝搬マトリックスを生成し，伝搬マトリックス保持部５２に保持する。層２パラメータにより伝搬マトリックス生成手段４２は層２の伝搬マトリックスを生成し，伝搬マトリックス保持部５３に保持する。層Ｎパラメータにより伝搬マトリックス生成手段４２は層Ｎの伝搬マトリックスを生成し，伝搬マトリックス保持部５２に保持する。
【００７７】
図１４は本発明の実施例２の伝搬特性演算のためのシステム構成を示す。図１４の構成の動作を説明する。伝搬特性演算手段４６において，送信波マトリックス，散乱マトリックス，回転マトリックス，層毎の伝搬マトリックスにより，伝搬特性の演算式Ａ（５−１），演算式Ｂ（６），演算式Ｃ（８−１）により送信電界に対する受信電界を演算する。伝搬特性演算手段４６において，４６１は演算式Ａ保持部であって，演算式（５−１）を保持するものである。４６２は演算式Ｂ保持部であって，演算式（６）を保持するものである。４６３は演算式Ｃ保持部であって，演算式（８−１）を保持するものである。
【００７８】
７１１は演算式Ａの演算結果を保持する演算式Ａ演算結果保持部である。７１２は演算式Ｂの演算結果を保持する演算式Ｂ演算結果保持部である。７１３は演算式Ｃの演算結果を保持する演算式Ｃ演算結果保持部である。８１は受信電界保持部であって，演算で求められた受信結果を保持するものである。８２は演算結果出力部であって，演算結果を保持するものである。図１４においては，図５にある各種図表生成手段８０，複屈折強度演算部９１，受信波出力換算手段７３等は図示を省略されている。  
【００７９】
図１５は本発明の実施例２の演算式Ａの演算のフローチャートであって，下向きの送信波の電界を求めるフローチャートである（層１に入力されて各層を伝搬して反射面で反射するまで伝搬する電磁波の電界を求めるフローチャート）。Ｓ１において，初期設定をする。Ｓ２において，各種パラメータを設定する。Ｓ３において，θの変更をする（前述の説明ではθは各層においてθ1 で一定であるとしたが，本実施例２では各層毎にθ（Ｘ軸と主軸の間の角度）を変更できるようにしてある。Ｓ４で，層ｉの伝搬マトリックス，回転マトリックスを求める。Ｓ５で層ｉ−１の演算結果を演算式Ａ演算結果保持部から求め，演算式Ａの送信電界とする。Ｓ６で演算式Ａにより層ｉを伝搬する電界を演算する。Ｓ７で層ｉの演算結果を演算式Ａ保持部に保持する。Ｓ８で層Ｎまで演算したか判定する。層Ｎまで演算をしていなければＳ１０でｉ＋１をｉとして演算対象の層を次の層にし，Ｓ３以降の処理を繰り返す。
【００８０】
図１６は本発明の実施例２の演算式Ｂの演算のフローチャートであって，反射面で散乱される電磁波の電界を求めるフローチャートである。Ｓ１で初期設定をする。Ｓ２で各種パラメータを設定する。Ｓ３で散乱マトリックス，回転マトリックスを求める。Ｓ４で層Ｎ−１の演算結果を演算式Ａ演算結果保持部から入力する。Ｓ５で演算結果Ａ保持部から入力した電界を送信電界として演算式Ｂ（式（６））の演算をする。Ｓ６で演算結果を演算式Ｂ保持部に保持する。
【００８１】
図１７は本発明の実施例２の演算式Ｃの演算のフローチャートであって，層Ｎで反射してから各層を伝搬する電磁波の電界を求めるフローチャートである。層Ｎから表面に向かって層を数えるものとする。Ｓ０において，初期設定をする。Ｓ１において，各種パラメータを設定する。Ｓ２において，θの変更をする（前述の説明ではθは各層においてθ1 で一定であるとしたが，本実施例２では各層毎にθ（ｘ軸と主軸の間の角度）を変更できるようにしてある。Ｓ３で，層Ｎから数えてｉ番目の層の伝搬マトリックスを求める。Ｓ４で層Ｎから数えてｉ−１番目の演算結果を演算式Ｃ結果保持部から求め，送信電界とする（反射した後の最初の層の送信電界は演算結果Ｂ保持部の値を送信電界とする）。Ｓ５で演算式Ｃを演算する。Ｓ６で層ｉの演算結果Ｃを保持する。層Ｎで反射してから層１まで伝搬した電磁波の電界を求めたか判定する。層１まで求まっていなかったらｉをｉ＋１としてＳ２以降の処理を繰り返す。Ｓ８で層１まで求まっていたら処理を終了する。
【００８２】
図１８は，本発明のシステム構成の実施例３の構成を示す。図１８は，パラメータ入力により演算した複屈折積層媒体の解析結果と該観測データに基づく観測データを表示出力し，観測データと本発明の装置による解析結果を比較して入力パラメータを変更し，実際の複屈折積層媒体の積層構造を解析することができるようにしたものである。
【００８３】
図１８において，２５は入力装置である。４６は伝搬特性演算手段である。５８は各種マトリックス生成手段であって，入力パラメータをもとに送信波マトリックス，伝搬マトリックス，散乱マトリックス，回転マトリックス，受信波マトリックスを生成するものである。５０は入力パラメータ保持部である。７１は演算結果保持部である。７２は出力データ作成部である。８０は各種図表生成手段である。９６は観測データ入力装置であって，レーダ等により氷床等の実際の複屈折媒体の層に電磁波を照射して観測した観測データを保持する磁気記録媒体である。９７は観測データ保存部であって，実際の観測データを保持するものである。
【００８４】
１００は観測データとシミュレーション結果の比較処理であって，オペレータ判断を表す。１０１は入力パラメータ変更の処理であり，入力装置２５により入力パラメータを変更する操作を表す。
【００８５】
図１８の構成の動作を説明する。最初にシミュレーションの動作について説明する。入力装置２５によりシミュレーションに必要なパラメータを入力する。入力されたパラメータは入力パラメータ保持部５０に保持される。各種マトリックス生成手段５８は送信波マトリックス，伝搬マトリックス，散乱マトリックス，回転マトリックス，受信波マトリックスを生成する（伝搬マトリックスはｋ演算部の演算結果をもとにマトリックス要素の配列に演算値が設定され，受信波マトリックスは演算結果がでた時点でマトリックス要素に具体的な値が設定される）。伝搬マトリックス要素生成手段６０は回転位相，減衰率を求める。伝搬特性演算手段４６は送信波マトリックス，伝搬特性マトリックス，散乱マトリックス，回転マトリックスに基づいて前述の各式（式（５−１），（６），（８−１），式（９）等）の演算をして，受信波の電界ｘ成分，ｙ成分を求め，演算結果保持部７１に保持する。出力データ作成部７２は受信波を電力に換算し，各種図表生成手段８０は，各種図表を生成する。出力装置９８は解析結果に基づく各種図表を表示出力する。
【００８６】
次に，観測データの表示方法について説明する。実際の観測で得られた観測データを保持する記録媒体をもとに観測データ入力装置９６により観測データを入力する。観測データは入力インタフェース２４１観測データ保存部９７に保持される。出力データ作成部７２は観測データを表示する図表を作成する。その図表は出力インタフェース２４２を介して出力装置９８に送られ，出力装置９８により観測データが表示される。
【００８７】
オペレータは観測結果の図表の表示と解析結果の表示を比較し，解析結果が観測結果にできるだけ近づくようにパラメータを変更する等の操作をする。入力装置２５により変更するパラメータを入力し，そのパラメータに基づくシミュレーションが再度なされ，表示され，実際の観測結果と比較される。
【００８８】
図１９は図１８の構成の動作のフローチャートである。Ｓ０で解析処理（シミュレーション）をするか，観測データの処理をするか選択する。解析による処理をする場合には，解析処理が開始され（Ｓ１），パラメータを入力する（Ｓ２）。入力されたパラメータが保持される（Ｓ３）。さらに，伝搬マトリックス要素が演算され（Ｓ４）。各種マトリックス要素が生成され（Ｓ５），伝搬特性が演算される（Ｓ６）。演算結果を保持する（Ｓ７）。出力データを作成し，各種図表を生成する（Ｓ８，Ｓ９）。
【００８９】
Ｓ０で観測データの処理を選択した場合には，観測データに基づく処理が開始される（Ｓ１１）。観測データを入力し，観測データを保持する（Ｓ１２，Ｓ１３）。出力する観測データを選択する（Ｓ１４）。選択した観測データに基づく各種図表を作成し，出力する（Ｓ１５，Ｓ１６）。
【００９０】
オペレータは，パラメータに基づく演算結果と観測データに基づく結果を比較し，パラメータを変更するか判断する（Ｓ１７，Ｓ１８）。Ｓ１９でパラメータを変更する場合には変更パラメータを入力するＳ２以降の処理を繰り返す。パラメータを変更しないのであればＳ２０で表示する観測データを変更するか判断し，変更するのであれば，保持されている観測データの中から，表示するデータを選択し，Ｓ１４以降の処理を繰り返す。入力する観測データを変更するのであればＳ１２以降の処理を繰り返す。
3.伝搬モデルによる構造解析の例
3.1 構造解析の条件
以下，本発明の実施例１の解析装置により解析した結果の例を示す。深さｒの面で散乱層があるとして，二層の場合について，散乱層の深さおよび，送信電界の主軸に対する方位角を変えながら，様々な誘電異方性の場合に受信電界の強度を求めた。
【００９１】
表１（図２８）は，送信波の偏波成分（ＥxT，ＥyT）および氷床の持つ電気的な性質をもつ要素（誘電率実数部，電気伝導度，散乱体の構造）に具体的なパラメータを与える場合分けについて示したものである。表１に従って（１）全く異方性の無い場合，（２）誘電率実数部に異方性を持たせた場合について基本的な複屈折計算を行った。さらに，（３）散乱体に異方性を与えた場合，（４）誘電率実数部と散乱体に異方性がある場合，（５）偏波成分の割合を変えた場合，（６）電気伝導度に異方性を持たせた場合について構造解析を行った。表１の計算番号の先頭または２番目数字の１から６までは，上記のそれぞれの場合を意味する。また，末尾のｘ，ｙは，受信電界の成分ＥxR，ＥyRを意味する。図１７，２０，２１，２２は，表１に従って計算を行った代表的な結果を示す。ここでは，本発明者等が実際に観測したアイスレーダで用いられた周波数１７９ＭＨｚでのレスポンスを求めた。各図において縦軸には氷床表面からの深さ，横軸に電波の偏波面を回転させた方位角θを示す。また，受信電力ＰＲは，次式のとおり受信電界強度ＥＲから求めたものを表示した。
【００９２】
ＰxR＝２０ｌｏｇＥxR      （ｄＢ）                (12)
ＰyR＝２０ｌｏｇＥyR      （ｄＢ）                (13)
3.2 異方性がない場合
図２６は，異方性のない場合に本発明の解析装置う使用して複屈折媒体を解析した結果を示す。１ｘは，送信波のｙ成分が０の直線偏波で誘電率実数部，電気伝導度，散乱体の異方性が無い場合である。この時，氷床とアンテナの向きの関数であるθを変化させてもｘ成分の受信電力は変化しないことが，この計算からも確認される。送信波のｘとｙの偏波成分に１と０．０５５（−２５ｄＢ）の割合を持たせた５１ｘの場合であっても異方性は生じていない。もともとｙ偏波成分のない1Yにおいてｙ偏波成分で受信した電界では，受信信号が０，すなわち交差偏波は生じない事が確認できる。いっぽう，５１ｙは，送信波にｙ成分があるため受信感度があったが，これはもとのｙ成分によるものである。
【００９３】
図２７は，１ｘで示す１方位角における深さ方向での減衰の様子を示したものである。ここでは，周波数，誘電率実数部を変化させても全く同じ変化を示す。これらの結果から，モデルがこれまでの偏波を考慮しない場合のアイスレーダの受信強度の説明と矛盾しない結果を示している。
3.3 誘電率実数部に異方性がある場合
図２０の２Ｘ，２Ｙ，５２Ｘ，５２Ｙは，ａ方向とｂ方向の誘電率実数部の値にわずかな異方性すなわち３．１５，０．１４８とした場合の計算である。２Ｘでは，約７５０ｍと約２２５０ｍの深さで受信電力が大きく低下する現象（ディップ）が現れた。このディップは，１Ｘの場合の受信強度に比べ約５０ｄＢ低い値を示した。また，出現する方位角は，４５°，１３５°，２２５°であって，９０°間隔である。図２１は，周波数毎にこのディップがどう出現するか誘電率実数部の値を，取りうる最大の異方性の差０．０３４を与え計算を行った。ここでは電気伝導度は，一定とした。周波数が高くなるにつれディップの深さ方向の間隔ｄｉは，狭くなる。図２２は，周波数とディップの間隔の関係について誘電率実数部の差の値を０．００８，０．０１７，０．０３４と変えて計算を行った結果である。周波数が高くなるほどｄｉは狭くなり，誘電率の異方性の差が大きいほどｄｉは狭くなる。図２０の２Ｙは，ｙ面の受信電力ＰyRの分布を示した。０°，９０°，１８０°，２７０°に約−３００ｄＢの大きな受信電力ＰyRの低下がおきる。このため他の角度分布がわかりにくいが４極の異方性を持つ構造を持っている。また，約１，５００ｍの深度の全方位角に約−１５０ｄＢの受信電力の低下が見られる。ｘ偏波成分に対しｙ偏波成分を−２５ｄＢ低下させた偏波成分を持たせた場合の異方性計算の結果は，５２Ｘ，５２Ｙで示す。５２Ｘ，５２Ｙは，２Ｘと比較するとディップの電力低下量が変化する。例えば，２ｘの７５０ｍ付近のディップでは，約７５ｄＢの受信電力の低下が５２Ｘでは約２４ｄＢとなる。５２Ｙでは，２Ｙのような極端な特性に対し両偏波成分が干渉した分布状態を示している。
【００９４】
図２２は誘電異方性がある場合に，誘電異方性の強さの異なる場合において，送信電波の周波数とディップｄｉを生じる深さの関係を求めたものである。
【００９５】
3.4 散乱体に異方性がある場合
図２３は散乱マトリックスに異方性がある場合に，本発明の解析装置で解析した結果を示す。散乱マトリックスＳの異方性が等方的ではなく一軸対称性であるとして，その設定を表の５３１Ｘ, ５３１Ｙ, ５３２Ｘ, ５３２ＹのようにそれぞれＳaaに対してＳbbを１０ｄＢ程度（０．３）と６ｄＢ程度（０．５）の低い値の異方性を持つ散乱体を想定して計算した。送信波の偏波成分を表１の様に与えた。図２３の５３１Ｘは，１０ｄＢ程度の異方性を持たせた場合であるが，受信電力ＰxRは，0 °方位と９０°方位の受信強度の絶対値が約１０ｄＢの異方性として出力された。５３２Ｘは，６ｄＢ程度の異方性を持たせた場合で，同様に，受信電力ＰXRは，約６ｄＢの異方性として出力された。５３１Ｙは，４極構造に近い分布となっているが，１８０°，３６０°の受信電力の低下する角がやや９０°，１７０°に接近している。この影響は，偏波成分の割合と散乱体の異方性の影響である。５３２Ｙは，受信電力の低下の影響が和らいだ４極構造である。
【００９６】
3.5   誘電率実数部と散乱体の異方性を組み合わせた場合
ここでは，３．３節と３．４節の計算で示した複屈折の効果と散乱体異方性を組み合わせた。表１の５４１Ｘ, ５４１Ｙ, ５４２Ｘ, ５４２Ｙの様に各要素を設定し計算を行った。送信波も表１の様に偏波成分を与えた。図２４の５４１Ｘは，やや散乱体の異方性の効果が強く７５０ｍ付近に複屈折によるディップがあらわれているが，全体的には２極の方位角依存性を持っている。また，図２０の２Ｘに比較してディップの方位角の４５°と１３５°がそれぞれ６０°と１０５°になり，２２５°と３１５°が２４０°と２８５°となっている。５４２Ｘは，図２０の２Ｘにやや近い分布を示しており，複屈折の影響が大きく４極構造を示している。また，５４１Ｘ，５４２Ｘは，散乱体の一軸対称性の割合が異なるためディップの方位分布が異なる。５４１Ｙ，５４２Ｙは，５３１Ｙ，５２Ｙに近い受信電力の分布を示している。
【００９７】
3.6   電気伝導度に異方性を持たせた場合
図２５は電気伝導度に異方性のある場合に本発明の解析装置を使用した得た演算結果の例である。電気伝導度の異方性の割合を２０％程度，σa に４．８×１０-6，σb に４×１０-6を与えて計算を行った。誘電率実数部の異方性は無く，偏波成分はｘのみ，散乱体の異方性は無く等方的とした。図２５の６Ｘは，0 °の方位の減衰係数が０．００７４ｄＢ／ｍで９０°方向の減衰係数は０．００８８ｄＢ／ｍの２極の異方性を示した。この場合は，減衰量の絶対値は，変化せず減衰係数のみ変化した。６Ｙは，４極構造を示した。
４．観測データとの比較
実際のアイスレーダアンテナから発射される電波は，アンテナの形状，垂直面および水平面指向特性，偏波特性，アンテナ近傍の状況（雪上車筐体や空気と氷床の境界面）によって，放射される電界に様々な影響を及ばせると考えられる。また，受信電界は，鉛直下向きの電界を受信するものとして計算している。よって，指向性のあるアンテナでは，斜めから散乱して返ってくる成分も合わせて受信することとなる。また，空気と氷床表面からの強力な電波反射により氷床内部からの散乱にも影響を及ぼす。ただし，指向性範囲内の氷床の反射面は，起伏がなく一様に平坦であると考えられるので，本モデルの様な平面波仮定が成り立つものと考える。また，氷床の温度特性の影響や地域性もあり，鉛直下向きに散乱帯等の構造が一定ではなく，水平面に対して深さ毎にねじれたり，急激な異方性の構造変化が生じている事も考えられる。本発明の構造解析装置は，異方性の効果を確認するために基本的な応答を求めたものであるが，アンテナの偏波，氷床の誘電率により現実的なパラメータを設定することにより，現実的な現象を説明することが可能である。
【００９８】
図２９は，ドームＦ（ＤＦ８０）（緯度７７°２２’，経度３９°３６’（南極）の地点）で実測した複屈折によるディップを示す。本発明者が実際に観測したデータに基づくものである（非特許文献２参照）。
【００９９】
偏波観測は，２２．５度ずつ１６方位について測定し，受信強度の深さ分布を示した。図２９（ａ）Ｔｘ｜｜Ｒｘは，送受信アンテナを平行に維持した平行偏波での方位角毎の受信強度を，受信機の検出限界−１１５ｄＢｍを考慮して−７０ｄＢｍからおよそ−１１５ｄＢｍまでの観測結果を示した。深さにして６００ｍから２０００ｍまでを示した。図２９（ｂ）Ｔｘ⊥Ｒｘは送受信アンテナを９０°偏波面の方位角を変えた交差偏波での方位角毎の受信強度の観測結果を同様に示す。周波数は１７９ＭＨｚ，パルス幅は１５０ｎｓｅｃ（実測地）を使用し，受信機感度は１ｄＢｍ単位で求めた。Ｄ１，Ｄ２は，受信電力が低下するディップの位置を示す。Ｄ1 は，深さ１２５０ｍから１５００ｍで方位角２２．５°，１１２．５°，２２９．５°に分布している。Ｄ２は，深さ１０００ｍ以下から検出限界まで広く分布し，方位角は６７．５°，１５７．５°，２４７．５°，３３７．５°に分布している。
【０１００】
次に伝搬モデルデータと比較する。伝搬モデル計算を行なった結果と実際の偏波観測で同様の結果が示されれば伝搬モデルの妥当性が証明されたことになり，氷床の誘電特性が特定できる。観測結果は平行偏波，交差偏波のいずれにもディップが発生している。このディップの分布の特徴は，図２９（ａ）のＤ1 の場合，一定の深さにディップが位置し，Ｄ１の方位とは４５°ずれて出現する。この結果は，Ｄ2 は，深さにあまり依存せず，Ｄ１の方位とは４５°ずれて出現する。この結果は，Hargreves （1977）（非特許文献１参照）が反射率に異方性を生じる場合，９０°毎に最小値が出現すると示している事実と同様であり，伝搬モデルで計算した複屈折の効果を示す結果（異方性がある場合の演算結果を示す図１７等）とも類似している。これらは，本発明の解析装置の有効性を示すものである。
【産業上の利用可能性】
【０１０１】
本発明によれば，複屈折率をもつ積層媒体の構造解析を容易に行なうことができる。特に，氷床の構造解析に応用した場合には，鉛直方向に伝搬する電波が，氷媒体の誘電率実数部の設定，電気伝導度の設定，散乱体の構造にどう反応するかを計算して評価することが可能であり，構造解析により，例えば，以下のようなことを評価する等が可能になる。（１）結晶主軸分布の平面内での異方性により氷床を構成する多結晶のマクロ的な誘電異方性が生じ，複屈折現象を起こすこと。（２）通常波と異常波の位相が逆転する際には，電波が互いにうち消し合い，レーダ波強度が大きく低下する現象が発生すること。（３）  （２）の現象は，電波の周波数および誘電異方性の割合に依存し，これを調査すれば氷床の誘電異方性の強度を遠隔探査できる。（４）一軸対称性の結晶構造をもつ氷床の散乱行列を計算することにより一義的に強い受信電力の異方性を示すこと。また，（５）電気伝導度の異方性が偏波としての電波伝搬に与える影響は実用上小さい値であったことが確認された。これらの結果を用いて実際の南極氷床における偏波観測の説明を試み矛盾なく説明できる。
【図面の簡単な説明】
【０１０２】
【図１】本発明の原理説明図である。
【図２】本発明のシステム構成の実施例１を示す図である。
【図３】本発明の各種マトリックス生成のための構成の実施例を示す図である。
【図４】本発明の伝搬マトリックス要素生成手段の構成を示す図である。
【図５】本発明の伝搬特性を求めるための構成を示す図である。
【図６】本発明の入力パラメータ保持部の構成を示す図である。
【図７】本発明の演算結果保持部の構成を示す図である。
【図８】本発明の送信波マトリックス生成手段と散乱マトリックス生成手段のフローチャートを示す図である。
【図９】本発明の回転マトリックス生成手段とｋ演算部のフローチャートを示す図である。
【図１０】本発明の伝搬マトリックス生成手段のフローチャートを示す図である。
【図１１】本発明の伝搬特性演算手段のフローチャートＡを示す図である。
【図１２】本発明の伝搬特性演算手段のフローチャートＢを示す図である。
【図１３】本発明の実施例２を示す図である。
【図１４】本発明の実施例２の伝搬特性演算のためのシステム構成を示す図である。
【図１５】本発明の実施例２の演算式Ａ演算のフローチャートを示す図である。
【図１６】本発明の実施例２の演算式Ｂ演算のフローチャートを示す図である。
【図１７】本発明の演算式Ｃの演算のフローチャートを示す図である。
【図１８】本発明のシステム構成の実施例２を示す図である。
【図１９】本発明の実施例３のフローチャートを示す図である。
【図２０】誘電率実数部に異方性がある場合のシミュレーションの結果を示す図である。
【図２１】周波数毎にこのディップがどう出現するか誘電率実数部の値を，取りうる最大の異方性の差0.034 を与えたシミュレーションの結果を示す図である。
【図２２】周波数とディップの間隔の関係について誘電率実数部の差をパラメータとしてシミュレーションを行なった結果を示す図である。
【図２３】散乱体に異方性がある場合のシミュレーションの結果う示す図である。
【図２４】誘電率実数部と散乱体に異方性がある場合のシミュレーションの結果を示す図である。
【図２５】電気伝導度に異方性がある場合のシミュレーションの結果を示す図である。
【図２６】異方性のない場合のシミュレーションの結果を示す図である。
【図２７】異方性のない場合に，周波数を変えてある方位角における深さ方向での減衰の様子をシミュレーションした結果を示す図である。
【図２８】表１を示す図である。
【図２９】観測データにおける複屈折のディップを示す図である。
【図３０】従来の観測装置の構成を示す図である。
【符号の説明】
【０１０３】
２１：ＣＰＵ
２２：メモリ
２３：各種演算手段（メモリ）
２４: 入出力インタフェース
２５: 入力装置
２６：ディスプレイ
２７：プリンタ
２８：外部記録媒体入力装置
２９：シミュレーションプログラム記録媒体                                        ３０：観測データ記録媒体
４１: 送信波マトリックス生成手段
４２：伝搬マトリックス生成手段
４３：散乱マトリックス生成手段
４４：回転マトリックス生成手段
４５：受信波マトリックス生成手段
４６：伝搬特性演算手段
４７：その他の演算手段
４８：複屈折率強度演算手段
【図１】
2005069785.tif 000028
【図２】
2005069785.tif 000029
【図３】
2005069785.tif 000030
【図４】
2005069785.tif 000031
【図５】
2005069785.tif 000032
【図６】
2005069785.tif 000033
【図７】
2005069785.tif 000034
【図８】
2005069785.tif 000035
【図９】
2005069785.tif 000036
【図１０】
2005069785.tif 000037
【図１１】
2005069785.tif 000038
【図１２】
2005069785.tif 000039
【図１３】
2005069785.tif 000040
【図１４】
2005069785.tif 000041
【図１５】
2005069785.tif 000042
【図１６】
2005069785.tif 000043
【図１７】
2005069785.tif 000044
【図１８】
2005069785.tif 000045
【図１９】
2005069785.tif 000046
【図２０】
2005069785.tif 000047
【図２１】
2005069785.tif 000048
【図２２】
2005069785.tif 000049
【図２３】
2005069785.tif 000050
【図２４】
2005069785.tif 000051
【図２５】
2005069785.tif 000052
【図２６】
2005069785.tif 000053
【図２７】
2005069785.tif 000054
【図２８】
2005069785.tif 000055
【図２９】
2005069785.tif 000056
【図３０】
2005069785.tif 000057
【手続補正書】
【提出日】平成15年10月9日(2003.10.9)
【手続補正１】
【補正対象書類名】図面
【補正対象項目名】図２０
【補正方法】変更
【補正の内容】
【図２０】
2005069785.tif 000058
【手続補正２】
【補正対象書類名】図面
【補正対象項目名】図２３
【補正方法】変更
【補正の内容】
【図２３】
2005069785.tif 000059
【手続補正３】
【補正対象書類名】図面
【補正対象項目名】図２４
【補正方法】変更
【補正の内容】
【図２４】
2005069785.tif 000060
【手続補正４】
【補正対象書類名】図面
【補正対象項目名】図２５
【補正方法】変更
【補正の内容】
【図２５】
2005069785.tif 000061
【手続補正５】
【補正対象書類名】図面
【補正対象項目名】図２６
【補正方法】変更
【補正の内容】
【図２６】
2005069785.tif 000062
【手続補正６】
【補正対象書類名】図面
【補正対象項目名】図２９
【補正方法】変更
【補正の内容】
【図２９】
2005069785.tif 000063