(57)【要約】【課題】時間的に離間した、ガスの存在又は濃度の計測値を処理する方法を提供し、及びこの方法により、設計条件を変更することなく応答の振舞いを改善すること。【解決手段】ガスセンサ装置（１００）に使用するための本発明の方法は、指数関数により時間特性を表現できる時間的に離間した計測値を処理する方法であって、計測値フィルタ（１０９）を使用する。時間的に離間した計測値を処理する方法を特定するために、時間間隔（Δt）で記録された検出器ユニット（１０８）の計測値が、同様に指数関数に従うが変更された時定数を有する、時間的に離間した出力関数（ｚ(t)）に変換される。出力値を計算するために、第１計測値と、第２計測値及び第１計測値間の加重差との合計が計算され、この合計値は出力値として使用される。【選択図】図２
2007298507.tif 000002
【特許請求の範囲】
【請求項１】
第１時定数を有する第１指数関数により時間特性を表現できる時間的に離間した計測値を処理する方法であって、
規定時間間隔で計測値を検出し、互いに対して前記規定時間間隔を有する、少なくとも第１時間における第１計測値及び第２時間における第２計測値を記憶する工程と、
前記第１時定数とは異なる第２時定数を有する第２指数関数により時間特性を表現できる時間的に離間した出力値を生成するために前記計測値をフィルタリングする工程と、
前記出力値を出力する工程とを具備し、
前記フィルタリング工程は、前記第１計測値と、前記第２計測値及び前記第１計測値間の加重差との合計であって、前記出力値を形成する前記合計を計算する工程を具備することを特徴とする時間的に離間した計測値の処理方法。
【請求項２】
前記加重差は、該加重差に一定係数を積算することにより形成されることを特徴とする請求項１記載の処理方法。
【請求項３】
前記加重差は、該加重差に調整可能な係数を積算することにより形成されることを特徴とする請求項１又は２記載の処理方法。
【請求項４】
前記調整可能な係数は、前記第２計測値及び前記第１計測値間の差の値に依存して決定されることを特徴とする請求項３記載の処理方法。
【請求項５】
時間的に離間した計測値を提供し、第１指数関数により時間特性を表現できると共に、第１時定数を有する、センサ用のフィルタであって、
該フィルタは、
互いに対して規定された時間間隔を有する、少なくとも第１時間における第１計測値及び第２時間における第２計測値を記憶する記憶ユニットと、
前記第１時定数とは異なる第２時定数を有する第２指数関数により時間特性を表現できる時間的に離間した出力値を生成する計算ユニットとを具備し、
該計算ユニットは、前記第１計測値と、前記第２計測値及び前記第１計測値間の加重差との合計であって、前記出力値を形成する前記合計を計算するよう構成されていることを特徴とするガスセンサ用フィルタ。
【請求項６】
放射を放出する少なくとも１個の放射源、計測される少なくとも１種類の検体を含む試料ガスで充填され得るガス計測室、前記放射を検出し、検体の存在及び濃度の一方又は両方に依存した検出器信号を発生する少なくとも１個の検出器、及び前記検出器信号を記録すると共に出力信号を発する制御ユニットを具備し、
該制御ユニットは、請求項１ないし４のうちいずれか１項記載の方法を実行するよう構成されたフィルタユニットを具備することを特徴とするガスセンサ装置。
【請求項７】
二酸化炭素を検出すること、二酸化炭素の濃度を決定することの一方又は両方を実行するよう設定されていることを特徴とする請求項６記載のガスセンサ装置。
【請求項８】
放射を放出する少なくとも１個の放射源、計測される少なくとも１種類の検体を含む試料ガスで充填され得るガス計測室、前記放射を検出し、検体の存在及び濃度の一方又は両方に依存した検出器信号を発生する少なくとも１個の検出器、及び前記検出器信号を記録すると共に出力信号を発する制御ユニットを具備するガスセンサ装置を作動させる方法において、
規定時間間隔で計測値を検出し、互いに対して前記規定時間間隔を有する、少なくとも第１時間における第１計測値及び第２時間における第２計測値を記憶する工程と、
前記第１時定数とは異なる第２時定数を有する第２指数関数により時間特性を表現できる時間的に離間した出力値を生成するために前記計測値をフィルタリングする工程と、
前記出力値を出力する工程とを具備し、
前記フィルタリング工程は、前記第１計測値と、前記第２計測値及び前記第１計測値間の加重差との合計であって、前記出力値を形成する前記合計を計算する工程を具備することを特徴とするガスセンサ装置の作動方法。
【請求項９】
二酸化炭素を検出すること、二酸化炭素の濃度を決定することの一方又は両方を特徴とする請求項８記載の作動方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、時間特性が指数関数で記載できる時間的に離間した計測値の処理方法、計測値フィルタの使用方法に関する。特に本発明は、ガスセンサ装置におけるこの種の使用方法、及び関連するガスセンサ装置に関する。
【背景技術】
【０００２】
公知のガスセンサは、放射を放出する放射源、計測される少なくとも１種類の検体を含む試料ガスで充填され得るガス計測室、及び検体の存在、濃度に依存した出力信号を発生する検出器を具備する。この種のガスセンサ装置は、例えば二酸化炭素又はメタン等の幅広い検体を検出するものとして知られている。従来のガスセンサは、赤外線（ＩＲ）波長範囲の放射を吸収する多くの極性ガスの特性に基づいている。赤外線光は、極性分子の双極子モーメントと相互作用することにより、回転及び振動の励起により分子の励起状態を生じさせることができる。赤外線光の熱エネルギーはこのようにしてガスに移転され、ガスを通過する赤外線光の強度が同様に減少する。励起状態に従い、例えば二酸化炭素の場合は4.25μmで問題のガスの特性である波長で吸収が生ずる。
【０００３】
特に二酸化炭素の検出は、現在では多数の用途領域において次第に重要になりつつある。例えば自動車の分野において、二酸化炭素の検出は、必要に応じ、すなわち二酸化炭素濃度の増加がある場合にのみ、適当な換気扇の起動により新鮮な空気の供給を開始するために、加熱及び空調のエネルギー効率を増加させるよう室内の空気の二酸化炭素含有量を監視するのに使用可能である。さらに、最新の自動車空調システムはまた、冷媒として二酸化炭素を使用しているので、自動車分野の二酸化炭素ガスセンサは、故障の際の漏れ二酸化炭素に関連する監視機能を満たすことができる。特に自動車分野において、ガスセンサは、頑丈さ、信頼性及び小型可能性の観点で最も高い要求事項を満たさねばならない。さらに、センサの応答時間は、安全用途の所定の制限値を超えてはならない。
【０００４】
特許文献１に開示されたような公知の検出器の出力信号は、次式に従った指数関数により濃度の急激な変化に対する応答として一般的に記述することができる。
【数１】
2007298507.tif 000003
ここで、ｙ(t)は時間tで検出器が放出する検出器信号を表わし、y0は検出器信号が近似する最終値を表わし、τは指数関数の時定数を表わす。
【０００５】
ガスセンサの応答の振舞いを特徴付けるために、実際には最も多くの部分で使用される、零における指数関数への正接を表わす時定数τではなく、検出器信号が最終値y0の90％に達する時間を表わすいわゆるt90時間である。通常、10秒未満のt90時間は、ガスセンサの安全用途で要求される。
【特許文献１】独国特許出願第１０２００５０５５８６９．７号明細書
【特許文献２】独国特許出願第１０２００５０５５８６０．７号明細書
【発明の開示】
【発明が解決しようとする課題】
【０００６】
しかし、以下で詳細に検討される種々の設計の境界条件は、公知のガスセンサの場合、t90時間が長すぎる結果をもたらすことが多い。
【０００７】
例えば、汚染からの保護の必要があって入口に膜を使用することは、原理上、内方への拡散遅延のため、ガスセンサの応答をさらに遅延してしまう。ガスセンサの小型化への厳しい要求により試料ガス用入口の寸法がさらに小さくなるので、同様に応答時間が増加してしまう。
【０００８】
センサの感度を向上させるために、特許文献２では回転対称である小部屋（cell）を使用することが提案されている。しかし、これは、ガス混合時間の延長、延いては応答時間の増加を意味する。このため、公知のガスセンサは、t90時間が安全用途の要求する制限値を超えるという不都合を有する。
【０００９】
他方、不十分な感度を受け入れなければならない等の不都合があるので、設計条件は可能な限り変更しないようにすべきである。
【００１０】
従って、本発明は、時間的に離間した計測値を処理する方法を特定し、この方法により、設計条件を変更することなく応答の振舞いを改善することを目的とする。
【課題を解決するための手段】
【００１１】
この目的は、独立請求項に記載された事項により達成される。本発明の有利な展開は従属請求項に記載される。
【００１２】
本発明は、ガスセンサ装置の改善された時間特性がフィルタの使用により達成できるという技術思想に基づいている。ここで、時間間隔Δtで記録された検出器ユニットの計測値が、同様に指数関数に従うが変更された時定数を有する、時間的に離間した出力関数に変換される。ガスセンサの指数階段関数応答及び作動時間離間モードのため、付加的な差部が計測信号に容易に導入できるので、このようなセンサシステムはより速い時定数を有する。
【００１３】
本発明によれば、出力値を計算するために、第１計測値と、第２計測値及び第１計測値間の加重差との合計が計算され、この合計値は出力値として使用される。これは、出力値が、以前の計測値と、現在の計測値及び加重係数Ｄを積算した以前の計測値の間の差との合計として計算される。この関係は以下の式（２）で記述される。
【数２】
2007298507.tif 000004
ここで、ｚ(ti)は出力値の時間離間順で現在の出力である値を表わし、ｙi-1は第１時間の検出器値を表わし、ｙiは、時間間隔Δtだけ第１時間と異なる第２時間の検出器値を表わし、Ｄは加重係数を表わす。
【００１４】
加重係数Ｄは、τ2を０に設定できる固定値にしてもよい。或いは、加重係数Ｄもまた、差Δyすなわち現在の計測値及び以前の計測値間の差に依存した関数にしてもよい。このように調整された加重係数を使用することにより、例えば新たな時定数τ2を設定することが可能であり、さらに例えば小さな変化がある場合に、ノイズが過度に増加しないようにすることが可能である。
【００１５】
例えば一定の加重係数としてＤ＝0.5が使用される場合、平均化と等しくなり、ノイズを低減する。
【００１６】
本発明は、ガス濃度の急激な増加に対するより速いセンサ応答が設計計測に頼る必要なく保証されるという利点がある。
【００１７】
本発明による計測値処理の有利な特性は特に、例えば自動車分野において、二酸化炭素の漏れを監視し客室内部の空気品質を確認するための二酸化炭素検出用のガスセンサ装置に活用できる。もちろん、本発明に従った原理はまた、他の任意のガス検出に関連して使用することができ、計測信号が第１時定数で指数関数に従う全てのセンサで重要であり、変更された時定数で指数関数に従う出力値が得られる。
【発明を実施するための最良の形態】
【００１８】
以下、添付図面に示された利点のある構成を参照して本発明を詳細に説明する。図面において、本発明に従った事項と同様又は対応する詳細には同じ参照符号が付される。
【００１９】
本発明のガスセンサ装置の構造及び時間離間して計測値を処理する本発明の作動方法を、図面を参照して詳細に説明する。
【００２０】
図１に示されるように、検出器ユニット１０８は、ここではほぼ階段関数として仮定できるガス濃度Ｃ(t)を検出する。しかし、言うまでもないが、検出器ユニット１０８用の入力信号はガス濃度とする必要はない。むしろ、本発明の原理によれば、任意のセンサの時間特性が影響を受ける。
【００２１】
検出器ユニット１０８は、第１指数関数ｙ(t)に従い、階段関数と比べて時定数τ1だけ遅延した応答を有する時間離間した検出器信号を供給する。一般式（１）から類推して、検出器信号の式は以下の式（３）で与えられる。
【数３】
2007298507.tif 000005
ここで、本発明に従った検出器信号は、時間間隔Δtで得られた時間離間した計測値ｙ(ti)からなると定義される。
【００２２】
より速いセンサ特性を達成できるようにするために、検出器信号は、本発明に従った出力信号として直接使用されないが、検出器信号ｙ(t)を出力信号ｚ(t)に変換するフィルタユニット１０９に供給される。図１に示されるように、出力信号ｚ(t)は、より短い過渡回復時間τ2を有するのが好適である。
【００２３】
検出器信号の出力信号への変換のための計算ルールは、本発明によれば式（２）により与えられる。
【００２４】
図２は、本発明に従ったフィルタユニットを使用するガスセンサ装置の概観を示す。図２に示されるように、本発明のガスセンサ装置１００は、本明細書では広帯域赤外線放射源である放射源１０２を具備する。原理上、ガスセンサ装置１００はいわゆるＮＤＩＲ（非分散型赤外線）センサである。赤外線放射源１０２を除き、基本部品は、ガス計測室、波長フィルタ１０６、及び検出器ユニット１０８としての赤外線検出器である。検出されるガス成分用にチェックされる試料ガス１１０は、入口及び出口１１２，１１４により符号化されているように、ガス計測室１０４内に送り込まれ又は拡散する。求められているガスの存在、濃度は、上述したように、赤外線範囲の特定波長の吸収により、光電的に決定することができる。図２に示されているように、求められているガスの濃度は、階段関数により時間範囲で近似できる。
【００２５】
放出される赤外線放射１１６は、ガス計測室１０４を通って検出器ユニット１０８に伝達される。検出器ユニット１０８に配置されているのは、検出されるガス分子が吸収する波長範囲を通ってのみ可能にする光フィルタである。他のガス分子は特定波長で通常、いかなる光も吸収しないので、検出器ユニット１０８に到達する放射量に影響を与えない。適当な全ての赤外線検出器は検出器ユニットとして使用可能であり、本発明に従った方法は任意のタイプの検出器に適応できる。
【００２６】
例えば、検出器は焦電要素、赤外線熱電対又はフォトダイオードであってもよい。いずれの場合においても、各要求事項に従って適当な検出器を選択することができる。このため、熱電対検出器とも称されるサーモコラム（thermocolumn）が選択されたスペクトル範囲で放射を特に高く一様な吸収するという利点を提供するのに対し、フォトダイオードは比較的低コストの部品を代表する利点を提供する。最後に、焦電センサは、非常に高い感度レベル及び小型化の製造可能性という利点を有する。
【００２７】
本発明によれば、赤外線信号は、所望の信号から熱的背景信号を濾波することができるように、放射源１０２によりパルス化される。このため、検出器ユニットにより供給される計測値は、図２に示されるように第１時定数τ1を有する指数関数を満たす時間離間した値ｙ(t)の形態で現れる。
【００２８】
コントローラ（制御ユニット）１２０は、一方で放射源１０２を始動し、他方で検出器ユニット１０８の検出器信号を受信すると共に本発明の原理に従ってさらにそれら信号を処理する。特に、コントローラは、検出器信号ｙ(t)を短い時定数τ2で出力信号ｚ(t)に変換するフィルタユニット１０９を具備する。図２に概略的に示されるように、フィルタユニット１０９は、フィルタ関数に従って変換が行われる計測値及び計算ユニット１２４を記憶するための記憶ユニット１２２を具備する。しかし、言うまでもないが、これら２ユニットもまた、例えば、コントローラ１２０で既に使用されたマイクロコントローラ等の記憶ユニット及び計算ユニットで実施できる。
【００２９】
図３は、図１及び図２に静的に表わされる時間関数Ｃ(t)、ｙ(t)及びｚ(t)を詳細に示す。このグラフを参照すると、加重係数Ｄは式（２）で以下のように推定される。
【００３０】
最初に、本発明は電流値及び最後の計測値のみが利用可能であるという基本条件から始まることに留意されたい。図３に示されるように、個々の計測値はそれぞれ、例えば、2.5秒等の規定された時間間隔Δtにより互いに分離される。ガス濃度Ｃ(t)は0〜100％の時間間隔Δt内で急激に変化するという以下の計算式が仮定される。ここで、図３に示される信号曲線は、明確化のためにｙ0=100に規格化されることを考慮されたい。曲線ｙ(t)は上述の式（３）に従い、既に一般的に公知であるように、時定数τ1は以下の式（４）によりt90時間で相関される。
【数４】
2007298507.tif 000006
【００３１】
これは、本発明に従った短いt90時間を得るために、過渡回復時間τが同様に短縮化されねばならない。出力値関数ｚ(t)として短縮されたτ2を有するこのような関数は、図３に概略的に表わされる。関数ｚ(t)は以下の式（５）に従う。
【数５】
2007298507.tif 000007
【００３２】
本発明によれば、得られる最終値ｙ0は、変化しないため、式（３）から式（５）の最終値ｙ0に対応する。基本指数関数特性は、式（１）にあるように同じ関数項で同様に提供される。本発明によれば、検出器ユニットにより達成されるτ1のみが式（５）の短い時間τ2で置換される。以下に推定されるように、前に計測された値ｙ(ti-1)を使用して現在の計測値ｙ(ti)から現在の出力値ｚ(ti)への変換は、式（２）により実行される。
【００３３】
以下の定義はこれのために使用される。すなわち、t＝ti-1は計測が実行される第１時間であり、t＋Δt＝tiは第２の計測が実行される、ti-1によりΔt後の第２時間である。２つの時間での計測値の差はΔyで表わされる。
【００３４】
この命名法及び式（３）及び式（５）を使用すると、加重係数Ｄの式は以下のように推定される。第１及び第２のそれぞれの時間で計測された第１及び第２の計測値については、以下のように書くことができる：
【数６】
2007298507.tif 000008
【数７】
2007298507.tif 000009
これにより、２計測値の差Δyは以下の通りになる。
【数８】
2007298507.tif 000010
次式と同様に式（５）から得られた式（６）及び式（８）が式（２）に代入されると、Ｄについての一般式が対応する解決の後、以下の結果になる。
【数９】
2007298507.tif 000011
【００３５】
実質的に境界線の場合は、τ2が0と仮定される場合、すなわち階段関数Ｃ(t)が正確に反復する理想的センサが達成される場合である。
【００３６】
0に接近するτ2についてＤ式の極限値が形成され、以下の簡略化した式が最早時間ｔに依存しないＤについての結果になる。
【数１０】
2007298507.tif 000012
例えばτ1に9.8秒が代入され、Δtに2.5秒が代入されると、加重係数Ｄは4.17の値になる。
【００３７】
任意のτ2に対しては、さらに変換が妥当であるので、Ｄの計算が可能である。第１式（６）は変形される結果、時間ｔはｙ0及びｙ1の関数として以下のように表わされる。
【数１１】
2007298507.tif 000013
両辺の対数をとると、
【数１２】
2007298507.tif 000014
式（１１）及び式（１２）の項ｅ−ｔ／τ1及びｔを式（９）に代入すると、加重係数Ｄについて以下の一般式が得られる。
【数１３】
2007298507.tif 000015
このため、出力値は、前の計測値に依存する加重係数に積算される。
【００３８】
図４及び図５は検出器の実際の計測値曲線４０１と比較して加重係数Ｄの異なる２つの値を使用した結果のシミュレーションを示し、図５は図４の拡大した詳細を表わす。ここで、曲線４０２はτ2→０の極限値にほぼ対応する加重係数Ｄ＝4.17に対応するので、曲線４０２は検出器信号４０１より事実上短い時定数を有する。しかし、シミュレーション結果はこの場合に若干行き過ぎを示すので、より小さいＤ値を使用する調節は理にかなっている。曲線４０３は、このような加重係数Ｄ＝2.23の出力値を表わす。
【００３９】
自動車用途でガスセンサを作動させる際の更なる基本的条件は、時間当たり1000ppmの濃度の変化で立ち上がる縁を少なくとも依然として検出可能でなければならないという安全要求事項である。他方、第１及び第２の計測値の差Δyが非常に小さい場合には、強すぎるノイズを防止するために、可能ならば加速された出力信号への変換が行われるべきでない。
【００４０】
Δyの絶対値に依存して種々の加重係数Ｄが意味をなすことができることが示される。これらの関係は図６に示される。或る規定値Δyminを下回る値の場合、より遅く応答して過度のノイズ及び行き過ぎを防止するセンサを構成するよう低いＤ値が便宜的に選択される。
【００４１】
規定値Δyminより上では、反応時間が早くなることが目的であるので、この最小値より上で、より高いＤ値を使用してフィルタリングが行われる。センサが遅く応答する傾向がある程度まで、最小値Δyは、以下の式（１４）を使用して近似的に計算することができる。
【数１４】
2007298507.tif 000016
【００４２】
時間当たり1000ppmの縁に到達し、２計測点間の間隔が2.5秒の場合、ｙs＝2500ppm及びτ1＝22.5で564ppmのΔymin値という結果となる。実際の時定数値τ1が依然として低い、例えば４０秒の場合、Δyminは335ppmとなる。
【００４３】
図６に示されるように、Δy＝Δyminについて本発明に従うと最も遅くとも加速されたシステム応答に導くＤ値が得られなければならない。これは、単純階段関数により例えば曲線６０１として示されるようにここで達成できるのみならず、曲線６０２から６０５を参照して以下に説明される、より複雑な依存性Ｄ（Δymin）を選択可能である。
【００４４】
階段関数６０１は、Δyをノイズ６１０（ここでは例示として値２τが選択される）から識別されるまで、一定値Ｄ＝0.5を使用する。Ｄ＝0.5を使用するフィルタリングは単純平均化に対応する。Δy＞Δyminの場合、Ｄ＝１とＤ＝１／（１−ｅ−Δt／τ1）との間の値はＤについて入力されるので、出力値は減少した時定数τ2を有する。
【００４５】
曲線６０２は、τ1の後、新たなＤ値がΔyに依存した直線的立ち上がりの形態でゆっくりと近似される場合の概略を示す。
【００４６】
このような直線依存性はまた、曲線６０３の概略としてΔy＝０で開始し、遅くともΔy＝Δyminで得られなければならないＤ＝１／（１−ｅ−Δt／τ1）の値である限り立ち上がるかもしれない。原理上、曲線６０５の曲線進行により示されるように、本発明に従ったフィルタ１０９の計算ユニット１２４において任意の依存性が実現可能である。
【００４７】
（図６の領域６０６により特徴付けられる）Ｄ＞１／（１−ｅ−Δt／τ1）のＤ値に対して、システムは強い行き過ぎを示すのに対し、領域６０７において、すなわち１＜Ｄ＜１／（１−ｅ−Δt／τ1）に対して、過渡回復時間における本発明による収縮が得られる。領域６０８において、すなわち１〜0.5のＤ値に対して、システムは、検出器が実際に規定する過渡回復時間よりも遅くなり、Ｄ＝0.5未満の値（領域６０９）に対して、センサは、曲線６０４により示されるように前に計測された値を保つ（hangs on）。
【００４８】
特にガスセンサに関連して本発明に従った計測値フィルタリングを使用すると、急激に増加するガス濃度値に対するセンサシステム全体のより速い応答を達成することができる。上記は通しで説明されているＮＤＩＲ二酸化炭素センサの場合であるが、最終値への指数近時が行われ、計測値が時間間隔の形態で現れる全センサシステムに適用可能であることは明らかである。
【００４９】
さらに、上述の通りに実行された計算は、初期値０から始まる計測信号が正の立ち上がりで最終値ｙ0に接近するという仮定に専ら基づいている。しかし、本発明に従った原理は、以下の一般式に従った指数崩壊の場合にも類推により移行することができる。
【数１５】
2007298507.tif 000017
この場合、ｂは崩壊信号の開始値を表わし、ａはｔ→∞の場合に信号が接近する最終値を表わす。
【００５０】
最後に、唯一の放射源及び唯一の検出器を有するガスセンサ装置は常に上述のように仮定される。当然、例えば２個の放射源又は２個以上の検出器等の全ての公知の基準原理をも使用することができる。例えば、図７は赤外線放射源の経年劣化の兆候を無くすために使用できる、２個の放射源を有するこのような装置を概略的に示す。
【図面の簡単な説明】
【００５１】
【図１】本発明のフィルタを有する検出器ユニットを示す概略図である。
【図２】本発明のフィルタ方法を使用するガスセンサユニットを示す概略図である。
【図３】時間の関数として実際に現れる検出器信号、出力値、及び濃度の変化の時間特性を示すグラフである。
【図４】種々の加重係数Ｄの時間特性を示すグラフである。
【図５】図４とは異なるメモリでの詳細を示すグラフである。
【図６】２つの検出器信号間の差に依存した種々の加重係数を示すグラフである。
【図７】２個の放射源を有するガスセンサ装置を示す概略図である。
【符号の説明】
【００５２】
１０２  放射源
１０４  ガス計測室
１０８  検出器
１０９  フィルタ
１１０  試料ガス
１２２  記憶ユニット
１２４  計算ユニット
ｙ(t)   検出器信号
ｚ(t)   出力信号
【図１】
2007298507.tif 000018
【図２】
2007298507.tif 000019
【図３】
2007298507.tif 000020
【図４】
2007298507.tif 000021
【図５】
2007298507.tif 000022
【図６】
2007298507.tif 000023
【図７】
2007298507.tif 000024