(57)【要約】【課題】単一モード光ファイバで構成される光ファイバ伝送路において単一モード光ファイバの実効ラマン利得係数を片端から、かつ簡単に評価できる単一モード光ファイバの実効ラマン利得係数評価装置を提供することにある。【解決手段】単一モード光ファイバ伝送路１６に波長λpの励起光パワーＰP1と波長λのＯＴＤＲ光をＷＤＭカップラーで合波して入射したときに、単一モード光ファイバ１６の任意の位置ｚでの後方散乱光強度Ｓ１（λ，ｚ，Ｐp1）と、励起パワーをＰP２に変化させた場合の後方散乱光強度Ｓ２（λ，ｚ，Ｐp2）とを測定する後方散乱光強度測定器１１とそれらの後方散乱光強度差と励起光波長λpでの損失係数、励起光パワーの差（Ｐp1−Ｐp2）を用い、単一モード光ファイバ１６の実効ラマン利得係数を演算する解析装置及び／又は調整用コンピュータ装置とを有するようにした。【選択図】図１
2008009036.tif 000002
【特許請求の範囲】
【請求項１】
励起光パワーＰp1が供給された単一モード光ファイバ伝送路の任意の位置ｚからの後方散乱光Ｓ１(λ,z,Pp1) (単位はdB)と、励起光パワーＰp2の場合の前記位置ｚからの後方散乱光Ｓ２(λ,z,Pp２)とを測定し、その差分値を計測し、前記励起光の波長λpにおける損失係数αpの値との相関において、実効ラマン利得係数ｇＲ/Ａeff  (ｇR：ラマン利得、Ａeff：実効断面積)として評価することを特徴とする単一モード光ファイバ伝送路の実効ラマン利得係数の評価方法。
【請求項２】
励起光パワーＰp1が供給された単一モード光ファイバ伝送路の任意の位置ｚからの後方散乱光Ｓ１(λ,z,Pp1) (単位はdB)、励起光パワーＰp2の場合の前記位置ｚからの後方散乱光Ｓ２(λ,z,Pp２)とを計測し、その差分値を計測し、および前記励起光の波長λpにおける損失係数αpの値を入力し、
【数１】
2008009036.tif 000022
を用いて演算し、実効ラマン利得係数ｇＲ/Ａeff  (ｇR：ラマン利得、Ａeff：実効断面積)を評価することを特徴とする請求項１記載の単一モード光ファイバ伝送路の実効ラマン利得係数の評価方法。
【請求項３】
請求項２において、実効ラマン利得係数ｇＲ/Ａeffと励起光パワーＰp1が供給された単一モード光ファイバ伝送路の任意の位置Ｚ0からの後方散乱光Ｓ１(λ,z,Pp1)と励起光パワーＰp2の場合の前記位置Ｚ0からの後方散乱光Ｓ２(λ,z,Pp２)の差分値Ｓ１(λ,z,Pp1)−Ｓ２(λ,z,Pp２)が相互作用長Ｌeffに比例する数（２）の関係から比例定数を求めて演算し、その値から実効ラマン利得係数ｇＲ/Ａeffを評価することを特徴とする単一モード光ファイバ伝送路の実効ラマン利得係数の評価方法。
【数２】
2008009036.tif 000023
【請求項４】
波長λpの励起光パワーＰpを供給した単一モード光ファイバ伝送路の位置ｚにおける波長λでの後方散乱光強度Ｓ(λ,z,Pp)を測定する機能と、該後方散乱光強度波形の演算処理を行なう機能とを有し、励起光パワーＰp1が供給された単一モード光ファイバ伝送路の任意の位置ｚからの後方散乱光Ｓ１(λ,z,Pp1) (単位はdB)と、励起光パワーＰp2の場合の前記位置ｚからの後方散乱光Ｓ２(λ,z,Pp２)とを計測し、その差分値を測定し、前記励起光の波長λpにおける損失係数αpの値との相関において演算し、実効ラマン利得係数ｇＲ/Ａeff  (ｇR：ラマン利得、Ａeff：実効断面積)として評価することを特徴とする単一モード光ファイバ伝送路の実効ラマン利得係数の評価装置。
【請求項５】
波長λpの励起光パワーＰpを供給した単一モード光ファイバ伝送路の位置ｚにおける波長λでの後方散乱光強度Ｓ(λ,z,Pp)を測定する機能と、該後方散乱光強度波形の演算処理を行なう機能とを有し、
励起光パワーＰp1が供給された単一モード光ファイバ伝送路の任意の位置ｚからの後方散乱光Ｓ１(λ,z,Pp1) (単位はdB)、励起光パワーＰp2の場合の前記位置ｚからの後方散乱光Ｓ２(λ,z,Pp２)を測定する工程および前記励起光の波長λpにおける損失係数αpの値を入力して
【数３】
2008009036.tif 000024
を用いて演算する工程を備え、実効ラマン利得係数ｇＲ/Ａeff  (ｇR：ラマン利得、Ａeff：実効断面積)を評価することを特徴とする単一モード光ファイバ伝送路の実効ラマン利得係数の評価装置。
【請求項６】
波長λpの励起光パワーＰpを供給した単一モード光ファイバ伝送路の位置ｚにおける波長λでの後方散乱光強度Ｓ(λ,z,Pp)を測定する機能と、該後方散乱光強度波形の演算処理を行なう機能とを有し、
励起光パワーＰp1が供給された単一モード光ファイバ伝送路の任意の位置ｚからの後方散乱光Ｓ１(λ,z,Pp1) (単位はdB)、励起光パワーＰp2の場合の前記位置ｚからの後方散乱光Ｓ２(λ,z,Pp２)を測定する工程および前記励起光の波長λpにおける損失係数αpの値を入力し、相互作用長Ｌeffに比例する数（４）の関係から比例定数を求めて演算し、その値から実効ラマン利得係数ｇＲ/Ａeffを評価することを特徴とする単一モード光ファイバ伝送路の実効ラマン利得係数の評価装置。
【数４】
2008009036.tif 000025
【請求項７】
請求項３を実施するために、前記励起波長λpにおける損失係数αpを測定するため手段として、波長λpの後方散乱光強度を測定する機能と、該後方散乱光強度波形の演算処理を行なう機能とを有することを特徴とする単一モード光ファイバ伝送路の実効ラマン利得係数の評価装置。
【請求項８】
請求項３を実施するために、ＯＴＤＲ装置、励起光源、信号と励起光とを合波するＷＤＭ合分波器、波形解析装置及びこれらを調整するコンピュータとからなる単一モード光ファイバ伝送路の実効ラマン利得係数の評価装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は単一モード光ファイバの実効ラマン利得係数の評価法及び評価装置に関し、特に、接続された複数の光ファイバで構成される光ファイバ伝送路の実効ラマン利得係数の評価法、及び装置に関する。
【背景技術】
【０００２】
昨今、インターネットサービスの急速な普及により、トラフィック量が急増しており、そのトラフィックを処理するために波長多重技術（ＷＤＭ）が用いられている。また、低損失光ファイバ、低損失波長帯を利用した波長域の開拓及び増幅技術の開発に伴い、光ファイバによる長距離伝送化が、さらに低コストで効率的な伝送を行うための低損失無中継伝送システムが検討され、光ファイバ伝送路を増幅媒体とする広帯域な光増幅技術の適用も考検討されている。
【０００３】
このような単一モード光ファイバ伝送路を用いた通信システムにおいて、分布ラマン増幅（Distributed Raman Amplifier：DRA）技術の開発が進んでいる。ラマン増幅とは石英ガラスで構成される光ファイバに信号光と、信号光より周波数が高い励起光とを同時に入力すると、石英ガラス中の誘導ラマン散乱現象により、励起光のエネルギーの一部が信号光に移る、即ち信号光が増幅される現象である。
【０００４】
分布ラマン増幅とは、光ファイバ伝送路自体を増幅媒体として用い、その伝送路に励起光を入射することによりラマン増幅効果を得る方法である。分布ラマン増幅を適用した光ファイバ伝送システムでは、伝送路の伝搬損失がラマン増幅で補償されるため、伝送可能距離を伸すことができる。
【０００５】
以下、従来の光伝送システムにおける測定方法について説明する。ラマン利得係数とは、各ファイバについて送信光のパワー１Ｗに対し、受信側で得られた利得を示すパラメータとして用いられる。
【０００６】
従来の実効ラマン利得係数の評価技術について説明する。伝送路自体を増幅媒体とする分布ラマン増幅のラマン利得係数の測定は、従来は、図６のような構成で行われてきた。図６は、従来のラマン利得係数測定方法例の構成図である。ここでは、ラマン利得係数の測定に用いられる信号光をテスト光と呼ぶ。伝送路の一方の端にテスト光源３１を備え、他方の端に励起光波長とテスト光波長を合分波するＷＤＭ(wavelength division multiplexer)カプラ３２を備える。ＷＤＭカプラ３２の励起光波長ポートに励起光源３３を備え、テスト光波長ポートにはテスト光パワーを測定するための受光器３４を備える。
【０００７】
受光器３４としては、たとえば光スペクトラムアナライザ、光パワーメータ等を用いる。テスト光を伝送路３６に入力し、励起光源出力を停止した状態で受光器にて検出されるテスト光パワーＰ１を測定する。次に励起光源を出力させた状態で受光器にて検出されるテスト光パワーＰ２を測定する。Ａは測定端３５を示す。デシベル表示でＰ２からＰ１を減算することで、テスト光が受けたラマン利得を求める。この利得を励起光出力パワーおよび実効断面積で除算することで、実効ラマン利得係数を算出する。
【０００８】
図７は、特許文献１に示されている実効ラマン利得係数を測定する別の方法を示す。
特許文献１では、光ファイバ伝送路４７に分布ラマン増幅を適用し信号光を増幅伝送する工程と、前記光ファイバ伝送路への励起光の出力を調整する励起光調整工程と、前記励起光を出力および停止した状態について前記光ファイバの伝播損失を時間領域光反射法（ＯＴＤＲ４５）によるテスト光によって測定する測定工程と、前記二状態における伝播損失の差分をもとに前記光ファイバのラマン利得を算出する算出工程を備える。
【０００９】
この方法は、光ファイバの一端に励起光が供給された前記光ファイバのラマン利得を測定する測定装置において、前記励起光を出力した状態における第１の戻り光パワーと、前記励起光を停止した状態における第２の戻り光パワーとの比が一定となる、前記光ファイバの他端を除く点における前記第１及び第２の戻り光パワーに基づいて、前記光ファイバの全長で生じるラマン利得を測定することを特徴とした方法である。尚、４１は励起光源、４２はＷＤＭ、４３は波長フィルタ、４４は光アッテネータ、Ａは測定端４６である。
【００１０】
しかしながら、この方法では、励起光のパワーを調整することが必要となり、測定におけるあいまいさがあり、また、任意の点での実効ラマン利得係数の分布が評価できない。
【特許文献１】特開２００４−２４０４６１号公報
【発明の開示】
【発明が解決しようとする課題】
【００１１】
実効ラマン利得係数測定の際、従来手法では、伝送路の両端に測定器、光源、及び作業者を配置して作業をしなければならないので、作業性が悪かった。このため、伝送路の一方の端での作業のみで実効ラマン利得係数を簡便に測定する手段が求められていた。
【００１２】
また、ＯＴＤＲを用いた片端で測定する方法も提案されているが、その評価法においては、励起光を入れた場合と入れない場合において、戻り光パワーが一定になるように励起光入力を調整する必要があり、測定精度の向上が求められていた。
【課題を解決するための手段】
【００１３】
主な第１の解決手段は、励起光パワーＰp1が供給された単一モード光ファイバ伝送路の任意の位置ｚからの後方散乱光Ｓ１(λ,z,Pp1) (単位はdB)と、励起光パワーＰp2の場合の前記位置ｚからの後方散乱光Ｓ２(λ,z,Pp２)との差値と、前記励起光の波長λpにおける損失係数αpの値との相関において、実効ラマン利得係数ｇＲ/Ａeff  (ｇR：ラマン利得、Ａeff：実効断面積)として評価することを特徴とする単一モード光ファイバ伝送路の実効的ラマン利得係数の評価方法とする。
【００１４】
更には、波長λpの励起光パワーＰp1が供給された単一モード光ファイバ伝送路の任意の位置ｚからの後方散乱光Ｓ１(λ,z,Pp1)、励起光パワーＰp2の場合の前記位置ｚからの後方散乱光Ｓ２(λ,z,Pp２)および前記励起光の波長λpにおける損失係数αpの値を用い、当該単一モード光ファイバ伝送路の任意の位置ｚにおける実効ラマン利得係数ｇＲ/Ａeff  式（１）に示す関係式より評価することにより、前記課題を解決する手段としている。
【数１】
2008009036.tif 000003
【００１５】
主な第２の解決手段は、実効ラマン利得係数ｇＲ/Ａeff  とを評価する際、波長λpの励起光パワーＰp1が供給された単一モード光ファイバ伝送路の任意の位置ｚからの後方散乱光Ｓ１(λ,z,Pp1)と励起光パワーＰp2の場合の前記位置ｚからの後方散乱光Ｓ２(λ,z,Pp２)が相互作用長Ｌeffが式（２）で示す関係式により記述できることより、この比例定数を求めることによって実効ラマン利得係数ｇＲ(z)/Ａeff  を求めることを、課題を解決する手段としている。
【数２】
2008009036.tif 000004
【００１６】
また、測定装置としても、ＯＴＤＲ装置、励起光源、信号光と励起光とを合波するＷＤＭ合分波器、波形解析装置及び／又はこれらを調整するコンピュータから構成される、きわめて簡単な装置構成で、前記課題を解決することができる。
【発明の効果】
【００１７】
本発明によれば、従来評価不可能であった敷設後の複数の光ファイバで構成される単一モード光ファイバ伝送路の実効ラマン利得係数を片端からかつ短時間に評価できるといった効果を奏する。
【００１８】
また、本発明を用いることにより、現場（中継点又は工場など）での位置で、実効ラマン利得係数が測定可能であり、機能としても簡易化されているので持ち運びなど利便性があり、簡単な装置構成で、簡単な評価法で実効ラマン利得係数を測定できる。
【発明を実施するための最良の形態】
【００１９】
以下に本発明の第１実施例について、図を用いて説明する。
【００２０】
（実施例１）
本発明の第１の実施例では、単一モード光ファイバもしくは光伝送路の任意の位置ｚにおける実効ラマン利得係数の評価手順について説明する。
図１は本発明による単一モード光ファイバ伝送路１６の実効ラマン利得係数の評価装置の構成を示す概略図である。本発明による単一モード光ファイバ伝送路１６の実効ラマン利得係数の評価装置は、後方散乱光強度測定装置１１、後方散乱光強度波形解析装置１２、励起用光源１３、光合分波器１４により構成される。
【００２１】
また、図２は、コンピュータを装備した場合の本発明１実施例装置の概略図である。普通には、図１の後方散乱光強度波形解析装置１２に演算機能を装備しているが、さらに、コンピュータ１７を設置し、測定した数値の記録、演算式による演算機能、各装置の調整機能を備える装置とすることもできる。
【００２２】
後方散乱光強度測定装置１１は、通常のＯＴＤＲ（Optical Time Domain Reflectometer）と同様の機能を有し、少なくとも１波長以上の測定光源を有する。後方散乱光強度測定装置１１では、長さＬｋｍの被測定単一モード光ファイバ、もしくは光ファイバ伝送路に波長λの測定パルス光を、測定端Ａに入射し、また、同時に、波長λpの連続光の励起光Ｐp(W)を入射し、該単一モード光ファイバもしくは単一モード光ファイバ伝送路の位置ｚからの後方散乱光強度、Ｓ(λ,z,Pp)（単位：dB）を測定する。
【００２３】
一方、後方散乱光強度波形解析装置では、以下に示す手順により、当該単一モード光ファイバ、もしくは単一モード光ファイバ伝送路の位置ｚにおける実効ラマン利得係数を評価する。散乱光までの距離ｚの位置は、ＯＴＤＲからの信号光である光パルスを被測定光ファイバに入射させ、各散乱点から戻ってくる反射光をＯＴＤＲで検出し、入射光と反射光の時間差から評価する。
【００２４】
まず、入射端から距離ｚにおける信号光パワーＰs(z)は次式の結合方程式を解くことによって求めることができる。
【数３】
2008009036.tif 000005
【数４】
2008009036.tif 000006
【００２５】
ここで、Ａeffは励起光と信号光とが相互作用する実効的な断面積であり、ｇRはラマン利得係数である。また、αpおよびαsはそれぞれ励起波長および信号波長における損失係数を表している。境界条件としてＰP(0)＝ＰPを与えて、式（４）を解くと
【数５】
2008009036.tif 000007
式（５）を式（４）に代入して、境界条件としてPs(0)=Psを与えて、実効ラマン係数ｇR/Ａeffが定数であると仮定して、信号光パワーＰs(z)を求めると、
【数６】
2008009036.tif 000008
【００２６】
ここで、
【数７】
2008009036.tif 000009
は相互作用長を表しており、式（７）のように距離ｚの関数として定義されている。
【数８】
2008009036.tif 000010
したがって、式（６）は、
【数９】
2008009036.tif 000011
式（８）は、任意の距離ｚにおける励起光に増幅された信号光のパワーを表す。
【００２７】
さて、ここで、ＯＴＤＲの測定原理に戻って距離ｚから後方に散乱される信号光について考える。
ｚ＝ｚでの信号光パワーＰs(z)は後方にαsＢ(z)Ｐs(z)分だけ散乱されるが、その散乱パワーと励起光とが相互作用をすることにより、式（８）と同様な増幅作用を経験しながら、入射端のほうに散乱パワーは増幅されながら伝搬する。したがって、任意の位置ｚから後方に散乱される後方散乱光パワーＰs(z)は次式で記述できる。
【数１０】
2008009036.tif 000012
ここで、αは散乱係数、Ｂ(z)は次式で定義される捕獲率を示す。
【数１１】
2008009036.tif 000013
【００２８】
また、励起光がない場合には、式（９）でＰp(0)＝0とおくと、よく知られた方程式が得られる。
【数１２】
2008009036.tif 000014
【００２９】
そこで、これらの関係式を利用した実効ラマン利得係数を求める方法を述べる。励起波長λpの励起パワーＰp1を供給した単一モード光ファイバあるいは単一モード光ファイバ伝送路にＯＴＤＲで、励起光と同じ片端から測定された後方散乱光強度をＳ１(λ,z,Pp1) (=10log(P(z))，単位：dB)とすると次式で表せる。
【数１３】
2008009036.tif 000015
【００３０】
一方、励起光がパワーをＰp2と変化させた場合、測定される後方散乱光強度をＳ２(λ,z,Pp２)は、
【数１４】
2008009036.tif 000016
したがって、式（１２）と式（１３）より、励起パワーを変化させた場合の後方散乱光強度を測定することにより、次式より実効ラマン利得係数をもとめることができる。
【数１５】
2008009036.tif 000017
【００３１】
したがって、実効ラマン利得係数は次式で求めることができる。
【数１６】
2008009036.tif 000018
式（１５）より、励起パワーを変化させて測定した後方散乱光と励起波長λpでの損失係数αpを測定することにより実効ラマン利得係数を評価することができる。
【００３２】
また、一方、式（１４）から、励起光パワーを変化させた場合の後方散乱光強度の差（Ｓ１(λ,z,Pp1)−Ｓ２(λ,z,Pp２)）は相互作用長Ｌeffに正比例することがわかる。したがって、後方散乱光強度の差と相互作用長Ｌeffの関係から、比例定数Ｋは最小二次近似法を用いて次式で求めることができる。
【数１７】
2008009036.tif 000019
これより、実効ラマン利得係数は、次式から簡単に求めることができる。
【数１８】
2008009036.tif 000020
【００３３】
（実施例２）
以下では、本発明による単一モード光ファイバ伝送路の実効ラマン利得係数の評価例について図表を用いて説明する。
後方散乱光強度測定装置は波長１５５５ｎｍの測定光源を有する。本実施では長さ２２ｋｍの１．３μｍ帯零分散ファイバ（ファイバＡ）と長さ２２ｋｍのノンゼロ分散シフトファイバ（ファイバＢ）を用いた。
【００３４】
図３に、ファイバＢに対して、励起光を入射した場合としない場合についての後方散乱光波形を示す。測定した光ファイバのパラメータは表１に示す。また、測定波長は１５５５ｎｍで、ＯＴＤＲは計の平均化時間は５分である。図３の破線は、励起光を入射していない場合の波形である。また、実線は、励起光を入射した場合の波形で、励起パワーがＰP2＝２４ｍＷとＰP1＝３３ｍＷの場合について示す。図からわかるように、励起光を入射していない場合には、後方散乱光強度はファイバ長とともに対数的に直線的に減少しているのがわかる。一方、励起光を入射すると、ＯＴＤＲ波形が変化しているのが明確に観測できる。この差が、ラマン散乱の影響である。
【００３５】
この時の、波長λ＝１５５５ｎｍにおける、測定端Ａからの後方散乱光強度波形、それぞれＳ１(λ,ｚ,Pp1)及びＳ２(λ,ｚ,PP2)を測定し、後方散乱光強度波形解析装置を用い、実施例１に記載の手順により被測定光ファイバ伝送路の実効ラマン利得係数を評価した。
【００３６】
表１に、図１に示した評価系で用いた各単一モード光ファイバの特性を示す。尚、表中の特性は、各単一モード光ファイバを接続前に、従来技術により当該光ファイバの一端、もしくは全長の平均値として評価した値を示す。
【表１】
2008009036.tif 000021
【００３７】
図４にファイバＡおよびファイバＢに対して波長λ＝１５５５ｎｍにおける実効ラマン利得係数を、前記、関係式（１６）を用いて評価した結果を示す。図４より、ファイバＡ及びＢの実効ラマン利得係数は、従来技術による評価結果と良く一致していることが分かる。
【００３８】
図３よりわかるように、ファイバＡとファイバＢとの実効ラマン利得係数にちがいがあることが明確に分かる。これは、ファイバＢの方が、比屈折率差ΔがファイバＡに比べて大きくかつモードフィールド径が小さいことによる。また、ここで評価された値は、Ｓｔｕａｒｔ  Ｇｒａｙによる研究発表(Stuart Gray; Raman gain measurements in optical fibers, SOFM, pp. 151-154, September 2000）で示されている値とよく一致しているのがわかる。
【００３９】
図５は、横軸（ｘ）に相互作用長（Leff）と縦軸（ｙ）に後方散乱パワー差との相関関係を示し、ファイバＢに関して、関係式（１４）に基づいて、Ｓ１(λ,z,Pp1)−Ｓ２(λ,z,Pp２)とＬeffの関係を示す。この図の傾きＫは式（１６）で与えられるので、式（１７）から実効ラマン係数を求めることができる。この図では、直線状を形成し、次の一次方程式が成り立つ。
ｙ＝−０．４４０３２＋０．１４１０７ｘ                      (１８)
【産業上の利用可能性】
【００４０】
上記のように、本発明によれば、従来評価不可能であった敷設後の複数の光ファイバで構成される単一モード光ファイバ伝送路の実効ラマン利得係数を片端からかつ短時間に評価でき、また、本発明を用いることにより、現場（中継点又は工場など）での位置で、実効ラマン利得係数が測定可能であり、機能としても簡易化されているので持ち運びなど利便性があり、簡単な装置構成で、簡単な評価法で実効ラマン利得係数を測定できる。
【図面の簡単な説明】
【００４１】
【図１】本発明による単一モード光ファイバ伝送路の実効ラマン利得係数評価装置の構成を示す概略図である。
【図２】単一モード光ファイバ伝送路の実効ラマン利得係数評価装置にコンピュータを備えた本発明の第１実施例の装置の概略図である。
【図３】本発明の第２実施例におけるファイバＢの、励起光パワーＰpを３３ｍＷ、２６ｍＷおよび０ｍＷの場合の後方散乱光波形を示す図である。
【図４】本発明の第２実施例におけるファイバＡおよびファイバＢの実効ラマン利得係数の評価結果を示す図である。
【図５】本発明の第２実施例におけるファイバＢのラマン利得係数の評価法を示す図である。
【図６】従来のラマン利得及びラマン利得効率測定方法例の構成図である。
【図７】従来の別のラマン利得及びラマン利得効率測定方法例の構成図である。
【符号の説明】
【００４２】
１１  後方散乱光強度測定装置
１２  後方散乱光強度波形解析装置
１３  励起用光源
１４  ＷＤＭカップラー
１５  測定端Ａ
１６  被測定単一モード光ファイバもしくは光ファイバ伝送路
１７  コンピュータ
３１  テスト光源
３２  励起用光源
３３  ＷＤＭカップラー
３４  受光器
３５  測定端Ａ
３６  被測定単一モード光ファイバもしくは光ファイバ伝送路
【図１】
2008009036.tif 000026
【図２】
2008009036.tif 000027
【図３】
2008009036.tif 000028
【図４】
2008009036.tif 000029
【図５】
2008009036.tif 000030
【図６】
2008009036.tif 000031
【図７】
2008009036.tif 000032