(57)【要約】【課題】  液体中に光を照射して液体中の微粒子により発生する動的散乱光を測定する装置において、装置のコストUPを抑制したうえで容易に液体の深さ方向の各位置における動的散乱光を測定する。【解決手段】  レーザ光源１１が出射した広帯域光を光カプラ１５により参照光と測定光に分割し、測定光を液体ＯＢ中に照射して粒子により発生した散乱光を光カプラ１５により参照光と合成させ、回折格子２２で波長により分光して撮像素子２３で受光する。撮像素子２３が出力する信号をセンサ信号取出回路４２がデジタルデータにしてコントローラ５０に出力し、コントローラ５０は撮像素子の位置に対応させた光強度データが記憶する。コントローラ５０は、光強度と対応する撮像素子の位置とを用いてフーリエ変換処理を行うことにより、液体ＯＢの深さ方向における各位置の光強度を算出し、算出した光強度を検出時間に対する関数にして各位置ごとに動的散乱光を算出する。【選択図】  図１
2014092525.tif 000002
【特許請求の範囲】
【請求項１】
光源から出射された光を参照光と測定光に分割し、前記測定光を粒子を含有する液体中に照射して前記粒子にて発生した散乱光を前記参照光と合成させて干渉光とし、前記干渉光の時間に対する光強度変化である動的散乱光を光強度測定手段により測定する動的光散乱測定装置において、
前記光源は広帯域光を出射し、
前記光強度測定手段は、前記干渉光を波長により分光する分光手段と、
前記分光手段により分光された光を撮像素子により受光し、前記撮像素子が出力する受光強度に対応する信号により、前記分光手段の分光方向における前記撮像素子のそれぞれの位置ごとの光強度を設定された時間間隔で検出する分光強度検出手段と、
前記分光強度検出手段により検出された、検出単位ごとの光強度と前記検出された光強度に対応する前記撮像素子の位置とを用いてフーリエ変換処理を行うことにより、前記液体における測定光の照射方向のそれぞれの位置ごとの光強度を算出する光強度データ計算手段と、
前記光強度データ計算手段により算出された前記液体における測定光の照射方向のそれぞれの位置ごとの光強度を、前記分光強度検出手段による検出時間に対する関数にすることにより、前記それぞれの位置ごとの動的散乱光を算出する動的散乱光データ作成手段とを備えたことを特徴とする動的光散乱測定装置。
【請求項２】
光源から出射された光を参照光と測定光に分割し、前記測定光を粒子を含有する液体中に照射して前記粒子にて発生した散乱光を前記参照光と合成させて干渉光とし、前記干渉光の時間に対する光強度変化である動的散乱光を光強度測定手段により測定する動的光散乱測定装置において、
前記光源は波長が一定の周期で変化するレーザ光を出射し、
前記光強度測定手段は、前記干渉光をフォトセンサにより受光し、前記フォトセンサが出力する受光強度に対応する信号により前記干渉光の光強度を設定された時間間隔で検出する受光強度検出手段と、
前記一定の周期を単位として、前記受光強度検出手段により検出された光強度と、前記光強度を検出したタイミングにおける前記光源が出射するレーザ光の波長とを用いてフーリエ変換処理を行うことにより、前記液体における測定光の照射方向のそれぞれの位置ごとの光強度を算出する光強度データ計算手段と、
前記光強度データ計算手段により算出された前記液体における測定光の照射方向のそれぞれの位置ごとの光強度を、前記受光強度検出手段による検出時間に対する関数にすることにより、前記それぞれの位置ごとの動的散乱光を算出する動的散乱光データ作成手段とを備えたことを特徴とする動的光散乱測定装置。
【請求項３】
請求項１または請求項２の動的光散乱測定装置において、
前記動的散乱光データ作成手段により算出された前記液体における測定光の照射方向のそれぞれの位置ごとの動的散乱光を用いて、前記それぞれの位置ごとに前記液体に含有される粒子の径を算出する粒径算出手段を備えたことを特徴とする動的光散乱測定装置。
【請求項４】
請求項２の動的光散乱測定装置において、
前記光源は波長が変化しないレーザ光をも出射でき、
前記光源から波長が変化するレーザ光を出射させ、前記干渉光として前記参照光との合成による正規干渉光を得る場合と、前記光源から波長が変化しないレーザ光を出射させ、前記干渉光として前記参照光との合成によらない仮干渉光を得る場合とのいずれかを選択する選択手段を備え、
前記光強度測定手段は、前記選択手段により仮干渉光を得る場合が選択された場合は、前記受光強度検出手段により設定された時間間隔で検出された光強度を動的散乱光とすることを特徴とする動的光散乱測定装置。
【請求項５】
光源から出射された光を参照光と測定光に分割し、前記測定光を粒子を含有する液体中に照射して前記粒子にて発生した散乱光を前記参照光と合成させて干渉光とし、前記干渉光の時間に対する光強度変化である動的散乱光を光強度測定工程により測定する動的光散乱測定方法において、
前記光源は広帯域光を出射し、
前記光強度測定工程は、前記干渉光を波長により分光する分光ステップと、
前記分光ステップにより分光された光を撮像素子により受光し、前記撮像素子が出力する受光強度に対応する信号により、前記分光ステップの分光方向における前記撮像素子のそれぞれの位置ごとの光強度を設定された時間間隔で検出する分光強度検出ステップと、
前記分光強度検出ステップにより検出された、検出単位ごとの光強度と前記検出された光強度に対応する前記撮像素子の位置とを用いてフーリエ変換処理を行うことにより、前記液体における測定光の照射方向のそれぞれの位置ごとの光強度を算出する光強度データ計算ステップと、
前記光強度データ計算ステップにより算出された前記液体における測定光の照射方向のそれぞれの位置ごとの光強度を、前記分光強度検出ステップによる検出時間に対する関数にすることにより、前記それぞれの位置ごとの動的散乱光を算出する動的散乱光データ作成ステップとを含むことを特徴とする動的光散乱測定方法。
【請求項６】
光源から出射された光を参照光と測定光に分割し、前記測定光を粒子を含有する液体中に照射して前記粒子にて発生した散乱光を前記参照光と合成させて干渉光とし、前記干渉光の時間に対する光強度変化である動的散乱光を光強度測定工程により測定する動的光散乱測定方法において、
前記光源は波長を一定の周期で変化させるレーザ光を出射し、
前記光強度測定工程は、前記干渉光をフォトセンサにより受光し、前記フォトセンサが出力する受光強度に対応する信号により前記干渉光の光強度を設定された時間間隔で検出する受光強度検出ステップと、
前記一定の周期を単位として、前記受光強度検出ステップにより検出された光強度と、前記光強度を検出したタイミングにおける前記光源が出射するレーザ光の波長とを用いてフーリエ変換処理を行うことにより、前記液体における測定光の照射方向のそれぞれの位置ごとの光強度を算出する光強度データ計算ステップと、
前記光強度データ計算ステップにより算出された前記液体における測定光の照射方向のそれぞれの位置ごとの光強度を、前記受光強度検出ステップによる検出時間に対する関数にすることにより、前記それぞれの位置ごとの動的散乱光を算出する動的散乱光データ作成ステップとを含むことを特徴とする動的光散乱測定方法。
【請求項７】
請求項４または請求項５の動的光散乱測定方法において、
前記動的散乱光データ作成ステップにより算出された前記液体における測定光の照射方向のそれぞれの位置ごとの動的散乱光を用いて、前記それぞれの位置ごとに前記液体に含有される粒子の径を算出する粒径算出ステップを含むことを特徴とする動的光散乱測定装置。
【請求項８】
請求項６の動的光散乱測定方法において、
前記光源は波長が変化しないレーザ光をも出射でき、
前記光源から波長が変化するレーザ光を出射させ、前記干渉光として前記参照光との合成による正規干渉光を得る場合と、前記光源から波長が変化しないレーザ光を出射させ、前記干渉光として前記参照光との合成によらない仮干渉光を得る場合とのいずれかを選択する選択ステップを含み、
前記光強度測定工程は、前記選択ステップにより仮干渉光を得る場合が選択された場合は、前記受光強度検出ステップにより設定された時間間隔で検出された光強度を動的散乱光とすることを特徴とする動的光散乱測定方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、微粒子を含有する液体に光を照射して液体中の微粒子により発生する動的散乱光を測定し、測定により得られた動的散乱光のデータを処理して液体中の微粒子の径を検出することが可能な動的光散乱測定装置および動的光散乱測定方法に関する。
【背景技術】
【０００２】
従来から、微粒子を含有する液体にレーザ光を照射して液体中の微粒子により発生する散乱光の強度の時間変化（以下、動的散乱光という）を測定し、この測定結果を基に液体中の微粒子の径を検出することが行われている。液体中の微粒子はブラウン運動をしており、微粒子の動きは粒子径が小さくなるほど速くなるため、微粒子の動く速さと微粒子の径の間には一定の関係がある。この測定方法は、動的散乱光から微粒子の動く速さの度合いである拡散係数を算出し、拡散係数から微粒子の径を求めるものである。以下に測定原理について説明する。
【０００３】
液体中の微粒子の動きを微粒子の拡散と考えると、拡散速度である拡散係数（Ｄ）と微粒子の半径（Ｒ）の間には以下の数１であるアインシュタイン・ストークスの式で示される関係がある。
【数１】
2014092525.tif 000003
ここでＫＢはボルツマン定数、Ｔは絶対温度、ηは液体の粘度である。すなわち、温度を一定にし、微粒子を含有する液体が定まったものであれば、拡散係数（Ｄ）と微粒子の半径（Ｒ）以外の値はすべて定数であり、拡散係数（Ｄ）から微粒子の半径（Ｒ）を求めることができる。そして、拡散係数（Ｄ）は、動的散乱光の測定結果を演算処理することで求めることができる。次にこの演算処理について説明する。
【０００４】
液体中の微粒子にレーザ光を照射し、発生する動的散乱光を測定すると、図９に示すようにランダムに変化する。これは、微粒子はブラウン運動をしているため、光が照射されている箇所に存在する微粒子から発生する散乱光の干渉の仕方が変化するためである。この動的散乱光の測定結果から時間相関関数を計算すると図１０に示すように指数関数的な減衰曲線になる。時間相関関数とは、全データにおける時間ｔのときの値と時間（ｔ＋τ）のときの値との相関係数が時間τの変化とともにどのように変化するかを示したものである。時間τが０のときは同じデータを比較するため相関係数は１であり、時間τが経過するとともに相関係数は小さくなっていくが、これは粒子が拡散していくと考えればいいので微粒子の径が均一であれば時間相関関数は指数関数になる。時間相関関数ｇ（τ）と拡散係数（Ｄ）の間には以下の数２〜数４で表される指数関数の関係がある。
（数２）
ｇ（τ）＝ exp（−Гτ）
（数３）
Г＝ｑ２Ｄ
（数４）
ｑ＝｛４π・ｎ・sin（Θ／２）｝／λ
ここでｎは液体の屈折率、λは照射するレーザ光の波長、Θは散乱光強度を測定する角度である。微粒子を含有する液体が定まったものであり、測定装置の構造が定まっていれば、これらの値は定数であり、数４から計算されるｑは定数である。よって、数３を数２に代入すると、時間相関関数ｇ（τ）から拡散係数（Ｄ）を求めることができることがわかる。そして、拡散係数（Ｄ）が大きくなるほど、すなわち粒子が速く動くほど時間相関関数ｇ（τ）の減衰曲線は急になることがわかる。
【０００５】
上記の説明で、動的散乱光の測定結果から時間相関関数ｇ（τ）を求め、時間相関関数ｇ（τ）から拡散係数（Ｄ）を求め、拡散係数（Ｄ）から微粒子の半径（Ｒ）を求めることができることがわかる。ただし、液体中の微粒子の径は均一ではなく分布をもつので、時間相関関数ｇ（τ）はそれぞれの径の微粒子の時間相関関数ｇ（τ）を加算したものとなる。よって、それぞれの径ごとの微粒子による散乱光強度を全粒子による散乱光強度で除算した値をＧ（Г）とすると、数２の時間相関関数ｇ（τ）の式は次の数５のように表すことができる。
【数５】
2014092525.tif 000004
数５より、ｇ（τ）を逆ラプラス変換すればＧ（Г）を求めることができることがわかる。そして、数１と数３からГは微粒子の半径（Ｒ）の逆数に定数を乗算したものであることがわかるので、Гを微粒子の半径（Ｒ）に変換することができ、図１１に示すように微粒子の半径（Ｒ）に対する存在割合Ｇ（Ｒ）を求めることができる。
【０００６】
上記の説明は、一度微粒子で散乱された光は再度他の微粒子で散乱されないで検出されるという仮定が適用される場合に成り立つ。よって、従来の方法では粒子の濃度が低い液体においては精度良く微粒子の径を測定することができるが、微粒子の濃度が高くなると精度良く微粒子の径を測定することができなくなるという問題がある。この問題を解決する手法として、特許文献１や特許文献２に示されるように、低コヒーレンス光源のマイケルソン型干渉計またはマッハツェンダー型干渉計を用いた装置による動的散乱光の測定が提案されている。この方法は、低コヒーレンス光を測定光と参照光に分割し、測定光のみを液体中に入射させ、液体中で微粒子により散乱された測定光と参照光とを合成させて光強度を検出することにより動的散乱光を測定するものである。これによれば、分割してから合成されるまでの測定光と参照光の光路長が等しくなって測定光と参照光が干渉するポイント（液体中の限られた深さ位置）での動的散乱光のみを測定することができる。よって、この方法を用いれば、微粒子の濃度が高くなっても精度良く微粒子の径を測定することができる。また、特許文献１や特許文献２に示されるように、この方法では参照光の位相を一定の周期で変調させ、動的散乱光をフーリエ変換して得られるデータから変調周期と等しい周期の信号を抽出したうえで時間相関関数を求めれば、S／Nが改善され、粒子の径をさらに精度よく測定することができる。
【先行技術文献】
【特許文献】
【０００７】
【特許文献１】特開２００３−１０６９７９号公報
【特許文献２】特開２０１１−１３１６２号公報
【発明の概要】
【０００８】
しかしながら、低コヒーレンス光源のマイケルソン型干渉計またはマッハツェンダー型干渉計を用いた装置による動的散乱光の測定には次の問題がある。
１）液体中の限られた深さ位置のデータしか得られないため、液体の深さ方向の各点で測定を行いたいときは、測定光の液体に対する出射位置を移動させるための駆動機構と出射位置の移動量を検出する手段が必要になって装置のコストが上がる。
２）測定光の光路長と参照光の光路長が液体中の測定深さ位置でほぼ同一になるように参照光の光路長を調整する必要があり、装置の調整に時間がかかる。
３）参照光の位相を一定の周期で変調させるため、ミラー等を高速で往復運動させる機構等が必要になり、装置のコストが上がる。
【０００９】
本発明はこれらの問題を解消するためなされたもので、その目的は微粒子を含有する液体に光を照射して液体中の微粒子により発生する動的散乱光を測定し、測定により得られた動的散乱光のデータを処理して液体中の微粒子の径を検出することが可能な動的光散乱測定装置および動的光散乱測定方法において、駆動機構や移動量検出手段を設けなくても液体の深さ方向の動的散乱光および微粒子の径が測定可能であり、参照光の光路長を調整する必要がなく、参照光の位相を一定周期で変調させる必要もない動的光散乱測定装置および動的光散乱測定方法を提供することにある。
【００１０】
上記目的を達成するために、本発明の特徴は、光源から出射された光を参照光と測定光に分割し、測定光を粒子を含有する液体中に照射して粒子にて発生した散乱光を参照光と合成させて干渉光とし、干渉光の時間に対する光強度変化である動的散乱光を光強度測定手段により測定する動的光散乱測定装置において、光源は広帯域光を出射し、光強度測定手段は、干渉光を波長により分光する分光手段と、分光手段により分光された光を撮像素子により受光し、撮像素子が出力する受光強度に対応する信号により、分光手段の分光方向における撮像素子のそれぞれの位置ごとの光強度を設定された時間間隔で検出する分光強度検出手段と、分光強度検出手段により検出された、検出単位ごとの光強度と検出された光強度に対応する撮像素子の位置とを用いてフーリエ変換処理を行うことにより、液体における測定光の照射方向のそれぞれの位置ごとの光強度を算出する光強度データ計算手段と、光強度データ計算手段により算出された液体における測定光の照射方向のそれぞれの位置ごとの光強度を、分光強度検出手段による検出時間に対する関数にすることにより、それぞれの位置ごとの動的散乱光を算出する動的散乱光データ作成手段とを備えたことにある。
【００１１】
これによれば、液体中の粒子にて発生した散乱光と参照光とが合成された干渉光が分光手段により波長により分光されて、分光強度検出手段が、それぞれの波長ごとに異なった撮像位置で撮像素子が受光することで、波長に対する光強度を検出する。そして、波長に対する光強度は、液体中の測定光の照射方向のそれぞれの位置で発生する散乱光と参照光とが合成された干渉光を足し合わせたものであり、波長に対する光強度の関係曲線の周期は、液体中における測定光の照射方向のそれぞれの位置に対応する。よって、光強度データ計算手段が撮像素子の位置と光強度の関係またはこの関係から算出した波長に対する光強度の関係をフーリエ変換処理することにより、液体における測定光の照射方向のそれぞれの位置ごとの光強度を算出することができる。そして、動的散乱光データ作成手段が算出されたそれぞれの位置ごとの光強度を、分光された光強度の検出時間に対する関数にすることにより、それぞれの位置ごとの動的散乱光を算出することができる。すなわち、これによれば、駆動機構や移動量検出手段を設けなくても液体の深さ方向のそれぞれの位置ごとの動的散乱光を測定することができる。また、参照光の光路長を調整する必要がないため、装置の調整に時間がかかることがない。さらに、撮像素子の位置と検出した光強度を用いてフーリエ変換処理を行うことでノイズを除去できるので、参照光の位相を一定周期で変調させる必要がない。
【００１２】
また、上記目的を達成するために、本発明の特徴は、光源から出射された光を参照光と測定光に分割し、測定光を粒子を含有する液体中に照射して粒子にて発生した散乱光を参照光と合成させて干渉光とし、干渉光の時間に対する光強度変化である動的散乱光を光強度測定手段により測定する動的光散乱測定装置において、光源は波長が一定の周期で変化するレーザ光を出射し、光強度測定手段は、干渉光をフォトセンサにより受光し、フォトセンサが出力する受光強度に対応する信号により干渉光の光強度を設定された時間間隔で検出する受光強度検出手段と、一定の周期を単位として、受光強度検出手段により検出された光強度と、光強度を検出したタイミングにおける光源が出射するレーザ光の波長とを用いてフーリエ変換処理を行うことにより、液体における測定光の照射方向のそれぞれの位置ごとの光強度を算出する光強度データ計算手段と、光強度データ計算手段により算出された液体における測定光の照射方向のそれぞれの位置ごとの光強度を、受光強度検出手段による検出時間に対する関数にすることにより、それぞれの位置ごとの動的散乱光を算出する動的散乱光データ作成手段とを備えたことにある。
【００１３】
これによれば、広帯域光を出射し干渉光を分光して撮像素子で受光することで波長に対する光強度の関係を検出する替わりに、波長が一定の周期で変化するレーザ光を出射し干渉光をフォトセンサで受光することで波長に対する光強度の関係を検出するようにしたので、同様に目的を達成できる。すなわち、これによっても、駆動機構や移動量検出手段を設けなくても液体の深さ方向のそれぞれの位置ごとの動的散乱光を測定することができる。また、参照光の光路長を調整する必要がないため、装置の調整に時間がかかることがない。さらに、フーリエ変換処理を行うことでノイズを除去できるので、参照光の位相を一定周期で変調させる必要がない。
【００１４】
また、本発明の他の特徴は、動的散乱光データ作成手段により算出された液体における測定光の照射方向のそれぞれの位置ごとの動的散乱光を用いて、それぞれの位置ごとに液体に含有される粒子の径を算出する粒径算出手段を備えたことにある。これによれば、液体の深さ方向のそれぞれの位置ごとに微粒子の径を測定することができる。
【００１５】
また、本発明の他の特徴は、上述の光源から波長が一定の周期で変化するレーザ光を出射する動的光散乱測定装置において、光源は波長が変化しないレーザ光をも出射でき、光源から波長が変化するレーザ光を出射させ、干渉光として参照光との合成による正規干渉光を得る場合と、光源から波長が変化しないレーザ光を出射させ、干渉光として参照光との合成によらない仮干渉光を得る場合とのいずれかを選択する選択手段を備え、光強度測定手段は、選択手段により仮干渉光を得る場合が選択された場合は、受光強度検出手段により設定された時間間隔で検出された光強度を動的散乱光とするようにしたことにある。これによれば、上述のように光源から波長が一定の周期で変化するレーザ光を出射して液体の深さ方向のそれぞれの位置ごとの動的散乱光を測定することができるとともに、液体中の微粒子の濃度が低ければ、光源から波長が変化しないレーザ光を出射して従来技術で説明した方法により動的散乱光を測定することができる。なお、この場合の仮干渉光とは実際は干渉光ではない光であり、液体中の微粒子にて発生した散乱光のみの光である。
【００１６】
さらに、本発明の実施にあたっては、本発明は、動的光散乱測定装置としての発明に限定されることなく、動的光散乱測定方法としての発明にも適用できるものである。
【図面の簡単な説明】
【００１７】
【図１】本発明が適用された第１実施形態による動的光散乱測定装置の全体構成図である。
【図２】第１実施形態による動的光散乱測定装置のコントローラが実行するプログラムのフロー図である。
【図３】第１実施形態による動的光散乱測定装置のコントローラが実行するプログラムのフロー図である。
【図４】フーリエ変換処理により液体中における測定光の照射方向位置に対する光強度を計算することができる原理を示す図である。
【図５】フーリエ変換処理したデータから液体中における測定光の照射方向位置ごとの動的散乱光を算出する処理を視覚的に示した図である。
【図６】本発明が適用された第２実施形態による動的光散乱測定装置の全体構成図である。
【図７】第２実施形態による動的光散乱測定装置のコントローラが実行するプログラムのフロー図である。
【図８】第２実施形態による動的光散乱測定装置のコントローラが実行するプログラムのフロー図である。
【図９】動的散乱光をグラフで示す図である。
【図１０】動的散乱光の時間相関関数をグラフで示す図である。
【図１１】動的散乱光の動時間相関関数から計算された粒子径の分布図である。
【００１８】
（第１実施形態）
以下、本発明が適用された動的光散乱測定装置の第１実施形態について図面を用いて説明する。図１は、本発明が適用された動的光散乱測定装置の全体構成図である。この動的光散乱測定装置は、容器ＶＥに入れられた微粒子を含有する液体ＯＢの深さ方向（液体ＯＢ中の測定光の照射方向）の複数の位置で動的散乱光を測定し、測定結果を処理して液体ＯＢの深さ方向の複数の位置における微粒子の径を算出する装置である。
【００１９】
動的光散乱測定装置は、測定部１０および光ヘッド３０を備えている。測定部１０は、レーザ光源１１、コリメーティングレンズ１２、集光レンズ１３、光カプラ１５、小径のコリメーティングレンズ１９，２１，固定ミラー２０，回折格子２２及び撮像素子２３を有する。そして、光カプラ１５と集光レンズ１３の間および光カプラ１５と小径のコリメーティングレンズ１９，２１との間は、光ファイバー１４，１６，１８で光が伝わるようになっている。また、光ヘッド３０はコリメーティングレンズ３１及びビームエキスパンダーを構成する集光レンズ３２、小径のコリメーティングレンズ３３を備える。そして、測定部１０の光カプラ１５と光ヘッド３０のコリメーティングレンズ３１との間は、光ファイバー１７で光が伝わるようになっている。
【００２０】
レーザ光源１１は広帯域光であるレーザ光を出射する。レーザ光源１１は、例えばスーパールミネセントダイオード（ＳＬＤ）で波長帯域が１００nm程度のものを用いる。コリメーティングレンズ１２は、レーザ光源１１からのレーザ光を平行光に変換し、集光レンズ１３は、コリメーティングレンズ１２からの平行光を集光して光ファイバー１４に入射させる。この場合、集光レンズ１３の焦点距離は、光ファイバー１４内に入射したレーザ光が光ファイバー１４内で全反射するように設定されている。光ファイバー１４に入射したレーザ光は、光カプラ１５に導かれる。光カプラ１５は、光ファイバー１４を介して入射されたレーザ光を２つに分岐させ、一方を固定ミラー２０に通じる光ファイバー１６に入射させ、他方を光ヘッド３０に通じる光ファイバー１７に入射させる。光ファイバー１６を伝わり光ファイバー１６の端から出射したレーザ光は、小径のコリメーティングレンズ１９により平行光になり、固定ミラー２０で反射して同じ光路を戻り、小径のコリメーティングレンズ１９を介して光ファイバー１６を伝わって光カプラ１５に戻る。以下、光カプラ１５で分岐されて固定ミラー２０に入射して反射し光カプラ１５に戻るレーザ光を参照光という。また、光ファイバー１７を伝わり、光ヘッド３０にある光ファイバー１７の端から出射したレーザ光はコリメーティングレンズ３１により平行光になり、集光レンズ３２、小径のコリメーティングレンズ３３により断面径を小さくされて、容器ＶＥに入れられた微粒子を含有する液体ＯＢに入射する。入射したレーザ光は、液体ＯＢ中の微粒子により液体ＯＢの深さ方向のすべての位置で散乱光を発生させ、散乱光の一部は小径のコリメーティングレンズ３３に入射して同じ光路を戻り、集光レンズ３２、コリメーティングレンズ１９を介して光ファイバー１７を伝わって光カプラ１５に戻る。以下、光カプラ１５で分岐されて液体ＯＢに入射し、散乱して光カプラ１５に戻るレーザ光を測定光という。
【００２１】
光カプラ１５まで戻った参照光と測定光は合成され、半分は光ファイバー１４を伝わって光源１１に戻るが、半分は光ファイバー１８を伝わって小径のコリメーティングレンズ２１で断面径の小さい平行光になり、回折格子２２に入射する。以下、光カプラ１５まで戻った参照光と測定光とが合成された光を干渉光という。回折格子２２に入射した干渉光は回折され、ＣＣＤまたはＣＭＯＳでできている撮像素子２３に入射する。回折格子２２で回折した干渉光は、波長によって強め合う回折角度が異なるため波長によって分光される。すなわち、撮像素子２３で受光された光は受光位置と波長が１：１の関係になっており、受光位置から波長を計算することができる。
【００２２】
動的光散乱測定装置は、光源駆動回路４０とセンサ信号取出回路４１とこれらの回路に指令を出力するとともにセンサ信号取出回路４１から入力するデジタルデータを演算処理するコントローラ５０およびコントローラ５０にデータ、条件および指令等を入力する入力装置５２とコントローラ５０で演算処理された結果や測定状況等を表示する表示装置５４を備えている。
【００２３】
光源駆動回路４０は、コントローラ５０から出射指令が入力すると設定された強度のレーザ光がレーザ光源１１から出射するよう、設定された強度の電流および電圧を供給する。そして、コントローラ５０から停止指令が入力すると電流および電圧の供給を停止する。また、センサ信号取出回路４２は、コントローラ５０から作動開始指令が入力すると、設定された時間間隔で撮像素子２３のそれぞれの画素が出力する信号の強度をデジタルデータにしてコントローラ５０に出力する。撮像素子２３のそれぞれの画素が出力する信号の強度は受光した光の強度に対応するため、コントローラ５０に出力されるデータは、それぞれの画素における受光強度である。データを入力したコントローラ５０は入力順にデータをメモリに記憶していくので、メモリの記憶領域からデータに対応する撮像素子２３の画素位置がわかる。
【００２４】
以下、動的光散乱測定装置の電源を入れ、微粒子を含有する液体ＯＢを入れた容器ＶＥを所定の位置に置いた後、コントローラ５０が実行するプログラムに沿って説明する。作業者は容器ＶＥを所定の位置に置いた後、入力装置５２から測定開始の指令を入力する。コントローラ５０は入力装置５２から測定開始の指令が入力すると、図２に示すフローチャートのデータ取込プログラムをステップＳ１０でスタートさせ、図３に示すフローチャートのデータ処理プログラムをステップＳ５０にてスタートさせる。
【００２５】
図２に示すフローチャートのデータ取込プログラムでは、コントローラ５０はステップＳ１２にて変数ｎを１に設定する。変数ｎは、センサ信号取出回路４２からデータを取り込んだ回数を示す。次にコントローラ５０は、ステップＳ１3にて光源駆動回路４０にレーザ光出射の指令を出力し、ステップＳ１４にて、センサ信号取出回路４２に作動開始の指令を出力する。これによりレーザ光源１１から設定された強度のレーザ光が出射され、センサ信号取出回路４２は撮像素子２３のそれぞれの画素が出力する信号の強度をデジタルデータに変換し、設定された時間間隔でコントローラ５０に出力することを繰り返す。次にコントローラ５０はステップＳ１６にて時間計測を開始し、ステップＳ１８にて経過時間が予め設定した微小時間Δｔに変数ｎを乗算した時間に達したか判定し、達していなければ「Ｎｏ」と判定し続けてステップＳ１８の判定処理を繰り返す。そして、経過時間がΔｔ・ｎの時間に達すると「Ｙｅｓ」と判定してステップＳ２０に進み、センサ信号取出回路４２が出力するデジタルデータを記憶する。そして、ステップＳ２２に進み、変数ｎが予め設定したデータ取り込み回数Ｎに達したかを判定し、達していなければ「Ｎｏ」と判定してステップ２４に進み、変数ｎをインクリメントしてステップＳ１８に戻る。このようにして、ステップＳ１８からステップ２４の処理を繰り返し実行することで、コントローラ５０のメモリには微小時間Δｔが経過するごとにセンサ信号取出回路４２が出力するデータが記憶されていく。そして、データを記憶するごとに変数ｎは１づつ増えていくので、変数ｎがデータ取り込み回数Ｎに達したときステップＳ２２で「Ｙｅｓ」と判定してステップＳ２６に進む。そして、ステップＳ２６にて時間計測を停止し、ステップＳ２８にてセンサ信号取出回路４２の作動を停止し、ステップＳ３０にて光源駆動回路４０への指令によりレーザ光出射を停止して、ステップＳ３２にてプログラムの実行を終了する。
【００２６】
図３に示すフローチャートのデータ処理プログラムでは、コントローラ５０はステップＳ５２にて、センサ信号取出回路４２が出力するデジタルデータの１回の取り込みであるデータ群が記憶されているか判定する。図３に示すフローチャートのプログラムは、図２に示すフローチャートのプログラムと同時にスタートするため、スタートしてからデータ群が記憶されるまでは「Ｎｏ」と判定され、ステップＳ５４における図２に示すフローチャートのプログラムが終了したかの判定も「Ｎｏ」と判定されて、ステップＳ５２、Ｓ５４の判定処理が繰り返される。やがて図２に示すフローチャートのプログラムの実行によりデータ群が記憶されると「Ｙｅｓ」と判定されて、ステップＳ５６に進む。コントローラ５０はステップＳ５６にて記憶されたデータ群から予め設定された記憶領域のデータを抽出する。前述のように、記憶領域は撮像素子２３の画素位置に対応するものであり、設定された記憶領域は、対応する画素位置が撮像素子２３における干渉光の分光方向に等間隔にあるものである。または、対応する画素位置から計算される干渉光の波長が等間隔で変化するようにしてもよい。そして、コントローラ５０はステップＳ５８にて、抽出するデータの記憶領域（画素位置）から予め計算されている干渉光の波長と抽出したデータの関係をフーリエ変換処理して、液体ＯＢ中の測定光の照射方向に対する光の強度のデータを生成する。
【００２７】
このフーリエ変換処理について原理を説明する。図４（ａ）に示すように、測定光の照射方向のある位置で測定光が反射した場合を想定する。このとき、図に示すように波長（λ）に対する光の強度（Ｉ）は正弦波状になる。ピーク位置は測定光の波長（λ）と測定光と参照光の光路差（Ｄｅ）がＤｅ＝ｍλの関係にあるとき（位相が合ったとき）であり、ボトム位置はＤｅ＝（ｍ＋１／２）λのとき（位相が１８０°ずれたとき）である（ｍは自然数）。よって、光の強度（Ｉ）を波長（λ）に対する式にすると、次のようになる。
（数６）
Ｉ ∝ ｃｏｓ（２π・Ｄｅ／λ）
この式からわかるように正弦波の周期はＤｅ／λとなり、波長（λ）により変化するが、広帯域光の光の波長は限られた範囲であるので、波長（λ）の変化は微量であり波長（λ）は定数とみなしてよい。よって、測定範囲においては正弦波の周期Ｄｅ／λは一定であり、光路差（Ｄｅ）に対応する値とみなしてよい。そして、測定光が反射する位置は図４（ｂ）に示すように、測定光の照射方向の複数の位置であるが、それぞれの反射位置において、波長（λ）に対する光の強度（Ｉ）は上記のように正弦波となり、それぞれの正弦波の周期は光路差（Ｄｅ）に対応する値である。そして、測定される波長（λ）に対する光の強度（Ｉ）の関係は、図４（ｂ）に示すようにこれらの光路差（Ｄｅ）に対応する周期の正弦波を重ね合わせたものであり、この関係をフーリエ変換すれば、光路差（Ｄｅ）に対応する周期と、光強度（Ｉ）の関係になる。そして、光路差（Ｄｅ）に対応する周期は、上述したようにＤｅ／λであるので、広帯域光の波長（λ）の中央値を乗算すれば、光路差（Ｄｅ）と光強度（Ｉ）の関係を求めることができる。光路差（Ｄｅ）の変化は、液体ＯＢ中の測定光の照射方向の位置の変化（液体ＯＢの深さ方向の位置の変化）であるので、これは液体ＯＢ中の測定光の照射方向の位置と光の強度の関係である。このフーリエ変換処理は、ＦＤ（フーリエドメイン）のＯＣＴの手法を使用すると考えてよい。
【００２８】
コントローラ５０はステップＳ５８の処理が終了すると、ステップＳ５２に戻り、データ群が記憶されているか判定するが、フーリエ変換処理の間に次のデータ群は記憶されているので「Ｙｅｓ」と判定して、再びステップＳ５６とステップＳ５８の処理を行う。このようにして、記憶されたデータ群である干渉光の波長に対する光強度が次々にフーリエ変換処理されて、液体ＯＢ中の測定光の照射方向の位置と光の強度の関係のデータが記憶されていく。そして、記憶されたデータ群がすべてフーリエ変換処理され、処理されるべきデータ群がなくなると、コントローラ５０はステップＳ５２にて「Ｎｏ」と判定し、ステップＳ５４へ進む。そして、処理されるべきデータ群がなくなったときは、既に図２に示すフローチャートのプログラムは終了しているので、ステップＳ５４にて「Ｙｅｓ」と判定され、ステップＳ６０へ進む。
【００２９】
コントローラ５０はステップＳ６０にて、フーリエ変換処理により得られたデータ群である液体ＯＢ中の測定光の照射方向の位置に対する光強度を、設定された位置ごとの時間経過に対する光強度のデータ群に並び替える処理を行う。これは、視覚的に示すと図５に示すように、フーリエ変換処理で得られたそれぞれの変数ｎごとの実線で囲まれたデータ群から点線で囲まれたデータ群を作成する処理である。これにより液体ＯＢ中の測定光の照射方向における各位置ごとの光強度の経過時間に対する変化が得られる。次にコントローラ５０は、ステップＳ６２にて動的散乱光があるデータ群のみを抽出し、さらにその中から後述するデータ処理を行うデータ群を抽出する。これは、液面に測定光が入射した後でしか動的散乱光は発生しないため、光強度が０に近い位置のデータ群は不要として除き、これにより液面の位置がわかるので、液面から予め設定した複数の位置（設定した複数の深さ）に対応したデータ群を抽出する処理である。この結果、液体ＯＢの液面から予め設定したそれぞれの深さ位置ごとに、図９に示す動的散乱光が得られる。
【００３０】
次にコントローラ５０は、ステップＳ６４にて、ステップＳ６２にて抽出した動的散乱光のデータ群の各データ群ごとに時間相関関数を計算する処理を行う。この処理は背景技術において説明したように、図１０に示すように、全データにおける時間ｔのときの値と時間（ｔ＋τ）のときの値との相関係数が時間τの変化とともにどのように変化するかを示す関数を作成するものである。次にコントローラ５０はステップＳ６６にて、ステップＳ６４にて算出した各データ群ごとの時間相関関数から背景技術において示した数５を用いて逆ラプラス変換処理を行い、数１、数３および数４を用いて、粒子径の分布を液体ＯＢの液面から予め設定したそれぞれの深さ位置ごとに算出する。
【００３１】
次にコントローラ５０は、ステップＳ６８にて、ステップＳ６６にて算出された粒子径の分布を、液体ＯＢの液面からの深さ位置に対応させて図１１に示すようなグラフにして表示装置５４に表示させる。作業者はこの表示を見ることにより、液体ＯＢのそれぞれの深さ位置における微粒子の径の分布状態がわかる。なお、液体ＯＢのそれぞれの深さ位置における粒子径の分布のピーク値を算出し、液体ＯＢの深さ方向距離と粒子径との関係を示したグラフを作成して表示装置５４に表示させるようにしてもよい。
【００３２】
上記のように動作する動的光散乱測定装置の実施形態によれば、レーザ光源１１から出射された広帯域光を光カプラ１５により参照光と測定光に分割し、測定光を粒子を含有する液体中ＯＢに照射して粒子にて発生した散乱光を光カプラ１５により参照光と合成させて干渉光とする。そして、干渉光を回折格子２２で回折させることにより波長により分光し、分光された光を撮像素子２３により受光し、撮像素子２３が出力する受光強度に対応する強度の信号をセンサ信号取出回路４２がデジタルデータにしてコントローラ５０出力する。コントローラ５０は入力するデジタルデータを設定された時間間隔で記憶することにより、撮像素子２３の分光方向におけるそれぞれの位置ごとの光強度を設定された時間間隔で検出したデータが記憶される。そして、コントローラ５０は、検出単位ごとの光強度と検出された光強度に対応する撮像素子の位置とを用いてフーリエ変換処理を行うことにより、液体ＯＢにおける測定光の照射方向のそれぞれの位置ごとの光強度を算出する。そして、コントローラ５０は、算出されたそれぞれの位置ごとの光強度を検出時間に対する関数にすることにより、それぞれの位置ごとの動的散乱光を算出する。これにより、駆動機構や移動量検出手段を設けなくても液体ＯＢの深さ方向のそれぞれの位置ごとの動的散乱光を測定することができる。また、参照光の光路長を調整する必要がないため、装置の調整に時間がかかることがない。さらに、撮像素子２３の位置と検出した光強度を用いてフーリエ変換処理を行うことでノイズを除去できるので、参照光の位相を一定周期で変調させる必要がない。
【００３３】
また、上記のように動作する動的光散乱測定装置の実施形態によれば、コントローラ５０は、算出された液体ＯＢにおける測定光の照射方向のそれぞれの位置ごとの動的散乱光を用いて、それぞれの位置ごとに液体に含有される粒子の径を算出する。これにより、液体ＯＢの深さ方向の各位置での微粒子の径を測定することができる。
【００３４】
（第２実施形態）
次に、本発明が適用された動的光散乱測定装置の第２実施形態について説明する。第２実施形態における動的光散乱測定装置は、第１実施形態が光源から広帯域光を出射し分光して波長に対応する光強度を検出したのに対し、光源から波長が一定周期で変化するレーザ光を出射し、干渉光の強度をレーザ光の波長に対応させて検出するものである。すなわち、第１実施形態がＦＤ（フーリエドメイン）におけるＳＤ−ＯＣＴの手法を使用するものであるのに対し、第２実施形態がＦＤ（フーリエドメイン）におけるＳＳ−ＯＣＴの手法を使用するものである。図６は、第２実施形態における本発明が適用された動的光散乱測定装置の全体構成図である。構成が第１実施形態と同じ箇所には同じ番号を付し、違う箇所には新たな番号を付してある。以下、第１実施形態と異なる箇所のみ説明を行う。
【００３５】
動的光散乱測定装置の測定部１０のレーザ光源２５は波長掃引レーザ光源であり、波長が一定の周期で変化するレーザ光を出射する。光源駆動回路４４は、コントローラ５０から掃引出射指令が入力すると、レーザ光源２５から所定強度の波長掃引レーザ光が出射するようにレーザ光源２５に波長掃引を意味する信号と所定強度の電流および電圧を供給する。また、レーザ光源２５は波長掃引範囲のいずれかの波長に波長を固定させたレーザ光を出射することもできる。光源駆動回路４４は、コントローラ５０から固定出射指令が入力すると、レーザ光源２５から所定強度の波長が固定されたレーザ光が出射するようにレーザ光源２５に波長固定を意味する信号と所定強度の電流および電圧を供給する。また、動的光散乱測定装置の測定部１０は光カプラ１５にて参照光と測定光が合成され光ファイバー１８の端から出射した干渉光をフォトディテクタ２４で受光する。フォトディテクタ２４は受光した光の強度に相当する強度の信号を増幅回路４６に出力し、増幅回路４６は入力した信号を設定された増幅率で増幅してＡＤ変換器４８に出力する。ＡＤ変換器４８は、コントローラ５０から作動開始指令が入力すると、設定された時間間隔で入力した信号の強度の瞬時値を表すデジタルデータを作成し、コントローラ５０に出力する。また、動的光散乱測定装置の測定部１０は、コリメーティングレンズ１９と固定ミラー２０との間にシャッタ２６を備える。シャッタ２６は表面が参照光の光軸に対してやや傾いており、シャッタ２６が閉じられると、反射した少量の参照光も光カプラ１５に戻らないようになっている。シャッタ２６は開閉制御回路４９から入力する駆動信号により開閉がなされ、開閉制御回路４９はコントローラ５０から入力する開の指令および閉の指令に基づいて駆動信号をシャッタ２６に出力する。
【００３６】
このように構成された動的光散乱測定装置において、作業者は容器ＶＥを所定の位置に置いた後、入力装置５２から波長掃引レーザ光による測定（すなわち、液体ＯＢの深さ方向位置ごとの粒子径測定）を行うか、波長固定レーザ光による測定（すなわち、液体ＯＢ全体の粒子径測定）を行うか入力し、測定開始の指令を入力する。波長固定レーザ光による測定は、従来技術で説明したように液体ＯＢ中の粒子の濃度が低い場合に有効である。コントローラ５０は入力装置５２から測定開始の指令が入力すると、図７に示すフローチャートのデータ取込プログラムをステップＳ１０でスタートさせ、図８に示すフローチャートのデータ処理プログラムをステップＳ５０にてスタートさせる。
【００３７】
図７に示すフローチャートのデータ取込プログラムは、図２に示すフローチャートのデータ取込プログラムの２箇所を変更したものである。１点はステップＳ１３−１’〜
Ｓ１３−５’として入力された測定の種類に基づいて光源駆動回路４４と開閉制御回路４９に指令を出力するステップを追加した点である。もう１点は、図２に示すフローチャートのデータ取込プログラムにおけるステップＳ１４とステップＳ２８のセンサ信号取出回路をＡ／Ｄ変換器に変えた点である。コントローラ５０はステップＳ１３−１’〜
Ｓ１３−５’にて、波長掃引レーザ光による測定が選択されていれば光源駆動回路４４に掃引出射指令を出力し、開閉制御回路４９に開の指令を出力する。一方、波長固定レーザ光による測定が選択されていれば、光源駆動回路４４に固定出射指令を出力し、開閉制御回路４９に閉の指令を出力する。
【００３８】
また、図８に示すフローチャートのデータ処理プログラムは、図３に示すフローチャートのデータ処理プログラムにステップＳ５１’，５９’として、測定の種類を判定するステップと、波長固定レーザ光による測定が選択された場合の処理のステップを追加し、ステップＳ５７’としてデータ（光強度）に対応する波長を計算するステップを追加したものである。コントローラ５０はステップＳ５１’にて測定の種類を判定し、波長掃引レーザ光による測定が選択されているときは、第１実施形態と同様にステップＳ５２以降の処理を行う。また、波長固定レーザ光による測定が選択されているときはステップＳ５９’に進んで、図７に示すフローチャートのデータ取得プログラムが終了するまで待つ。そして、データ取得プログラムが終了するとステップ６４に進んで時間相関関数の計算を行う。これは、従来技術で説明した方法による粒子径の測定である。一方、ステップＳ５２に進んだ場合は、図３に示すフローチャートのデータ処理プログラムと殆ど同じ処理が行われるが、異なっている点はステップＳ５７’として取得したデータ（光強度）からデータに対応する波長を計算する処理を行う点である。すなわち、第１実施形態ではデータの記憶領域（撮像素子２３の画素位置）から光強度に対応する波長が分かったが、第２実施形態ではこれができないため、Ｓ５７’においてデータ群のそれぞれのデータ（光強度）から波長を求める処理を行う。これは、コントローラ５０にはレーザ光源２５が出射するレーザ光の波長とレーザ強度の特性曲線が予め記憶されており、この特性曲線とデータ群のそれぞれのデータから作成した曲線が最も合うように当てはめることで、それぞれのデータに対応する波長を求めるものである。これ以外のプログラム処理および装置構成は第１実施形態と同じである。
【００３９】
上記のように動作する動的光散乱測定装置の実施形態によれば、広帯域光を出射し干渉光を分光して撮像素子で受光することで波長に対する光強度の関係を検出する替わりに、波長掃引レーザ光源２５から波長が一定の周期で変化するレーザ光を出射し、干渉光をフォトディテクタ２４で受光することで波長に対する光強度の関係を検出するようにしたので、同様に目的を達成できる。すなわち、これによっても、駆動機構や移動量検出手段を設けなくても液体の深さ方向のそれぞれの位置ごとの動的散乱光を測定することができる。また、参照光の光路長を調整する必要がないため、装置の調整に時間がかかることがない。さらに、フーリエ変換処理を行うことでノイズを除去できるので、参照光の位相を一定周期で変調させる必要がない。また、コントローラ５０が液体ＯＢ中の測定光の照射方向におけるそれぞれの位置ごとに動的散乱光を用いて液体に含有される粒子の径を算出するので、液体ＯＢの深さ方向の各位置での微粒子の径を測定することができる。また、波長掃引レーザ光による測定と波長固定レーザ光による測定を選択できるようにしたことから、液体ＯＢの深さ方向のそれぞれの位置ごとの微粒子の径を測定することができるとともに、従来技術で説明した方法により液体ＯＢ全体としての微粒子の径を測定することができる。
【００４０】
以上、本発明の実施形態について説明したが、本発明の実施にあたっては、上記実施形態に限定されるものではなく、本発明の目的を逸脱しない限りにおいて種々の変形も可能である。
【００４１】
上記第１実施形態および第２実施形態では、測定光と参照光への分割および合成を行う構造として、光ファイバーに光を入射させ光カプラにより光の分割と合成を行う構造を用いたが、これに替えて光を断面径の小さい平行光にしてビームスプリッタに入射させ、ビームスプリッタにより光の分割と合成を行う構造にしてもよい。これによっても上記実施形態と同様の効果が得られる。
【００４２】
また、上記第１実施形態および第２実施形態では、測定光と参照光への分割および測定光と参照光の合成を行う干渉計の構造としてマイケルソン型干渉計を用いたが、測定光と参照光への分割および測定光と参照光の合成が可能であれば干渉計はどのような構造でもよい。例えばマッハツェンダー型干渉計を用いてもよい。
【００４３】
また、上記第１実施形態では、回折格子により干渉光を分光させたが、干渉光を波長により分光させることができるならばどのような方法を用いてもよい。例えばプリズムを使用して分光を行うようにしてもよい。
【００４４】
また、上記第１実施形態では、１つのレーザ光源が出射する広帯域光を使用したが、波長範囲を広くしたいときは、波長範囲が異なる複数の光を合成して使用するようにしてもよい。
【００４５】
また、上記第２実施形態では、シャッタの開閉により光カプラにて参照光との合成を行うか否かを選択できるようにしたが、参照光との合成を行うか否かを選択することができればどのような方法を用いてもよい。例えば、レーザ光源と光カプラの間の光ファイバ、光カプラと光ヘッドの間のの光ファイバおよび光カプラとフォトディテクタとの間の光ファイバに光切換スイッチをそれぞれ設け、光切換スイッチ間を別の光ファイバで連結させ、光切換スイッチでの切り換えによりレーザ光源からの出射光をそのまま光ヘッドに導き、光ヘッドからの光をそのままフォトディテクタに導くようにしてもよい。
【００４６】
１０…測定部、１１…レーザ光源、１５…光カプラ、１５…受光センサ、１４，１６，１７，１８…光ファイバー、２０…固定ミラー、２２…回折格子、２３…撮像素子、２４…フォトディテクタ、２５…波長掃引レーザ光源、２６…シャッタ、３０…光ヘッド、４０，４４…光源駆動回路、４２…センサ信号取出回路、４８…Ａ／Ｄ変換器、４９…開閉制御回路、５０…コントローラ、ＯＢ…微粒子含有液体、ＶＥ…容器
【図１】
2014092525.tif 000005
【図２】
2014092525.tif 000006
【図３】
2014092525.tif 000007
【図４】
2014092525.tif 000008
【図５】
2014092525.tif 000009
【図６】
2014092525.tif 000010
【図７】
2014092525.tif 000011
【図８】
2014092525.tif 000012
【図９】
2014092525.tif 000013
【図１０】
2014092525.tif 000014
【図１１】
2014092525.tif 000015