(57)【特許請求の範囲】
【請求項１】
複数枚のランナ羽根を有するランナと、このランナの外周側に円形翼列状に配置された複数枚のガイドベーンとを具備する水力機械において、上記ランナの外縁部に近接した静止部に設置されランナ振動の信号を発生する非接触振動センサと、上記非接触振動センサの上記ランナ振動信号から、Ｚｒ・Ｎ（Ｚｒは上記ランナ羽根の枚数であり、Ｎは回転同期周波数である）の周波数成分及びその高次調波成分を抽出してそれぞれの周波数成分信号を発生する周波数分析装置と、上記ランナ円周上の所定の位置が上記静止部の所定の位置に達した時にトリガー信号を発生するトリガー装置と、上記トリガー信号に対する上記抽出された各周波数成分信号の初期位相角を検出する初期位相角検出装置と、上記周波数分析装置の各周波数成分信号の大きさと上記初期位相角検出装置によって検出された各周波数成分信号の初期位相角とに基づきランナ振動モードの変位パターンを算出しそれらを合算する演算装置とを具備することを特徴とする水力機械。
【請求項２】
上記ガイドベーンの枚数をＺｇ、回転軸回りの角度座標をランナ回転方向にφ、上記周波数成分信号のｍ次成分における片振幅の大きさをＡｍ、時間をｔ、上記初期位相角をΨｍとした時に、上記演算装置は、上記ランナ振動モードの変位パターンＸｍを、
ｋ＝ｍ・Ｚｒ−Ｚｇ
Ｘｍ＝Ａｍ・ｓｉｎ［２π・Ｚｇ・Ｎ・ｔ−ｋ｛φ−Ｚｇ・Ψｍ／ｋ
（Ｚｇ＋ｋ）｝］
から算出し、ランナの実振動モードパターンＸｒを、
【数１】
0003756531.tif 000011
として合算することを特徴とする請求項１に記載の水力機械。
【請求項３】
上記非接触振動センサは上記ランナの上カバー及び下カバーに夫々設置されていることを特徴とする請求項１に記載の水力機械。
【請求項４】
上記演算装置の出力を判別するパターン判別装置と、このパターン判別装置の出力に基づき水力機械を停止させる制御装置とを具備することを特徴とする請求項１に記載の水力機械。
【発明の詳細な説明】
【０００１】
【産業上の利用分野】
本発明は、水力機械に係り、特にランナの振動応力を低減した水車やポンプやポンプ水車等の水力機械に関する。
【０００２】
【従来の技術】
近年水車やポンプ水車は、高落差化及び高速化が進み、これに伴って種々の振動問題が発生するようになっている。特に、揚水発電所に使用されるポンプ水車は、経済性の観点から高速化が著しく、ランナに非常に大きな振動が発生する傾向にある。
【０００３】
図９は、一般的なフランシス形のポンプ水車の構造を示したもので、図示を省略した発電電動機に直結された主軸１の下端には、ランナ２が連結され、このランナ２は、クラウン３とバンド４とこれらの間に挟持された複数枚のランナ羽根５とから構成される。ランナ２は上面及び下面が夫々、上カバー６及び下カバー７によって覆われている。ランナ２の外周部には複数枚のガイドベーン８が円形翼列状に配置され、これらのガイドベーン８の外周には複数枚のステーベーン９が円形翼列状に配置されている。これらのステーベーン９の外周にはケーシング１０が配置されている。また、ランナ２の下方には吸い出し管１１が連接されている。
【０００４】
このような構成のフランシス形ポンプ水車は、水車運転の場合にはケーシング１０内の高圧水が、ステーベーン９を介してガイドベーン８に流入し、このガイドベーン８によって流量制御されてランナ２に流入し、このランナ２を回転駆動した後に、吸い出し管１１を通って排出される。
ランナ羽根５とガイドベーン８との翼列干渉に起因する振動がランナ２及び静止部に発生し、特に高落差高速機では振動周波数が高くなるため、ランナ２が共振してランナ羽根５に高い繰返し変動応力が作用し、材料疲労限界を越えると、ランナ２の疲労破壊に至るといった問題がある。
【０００５】
そこで、最近では新設の高落差高速ポンプ水車の営業運転開始に際しては、現地試運転時に、運転中のランナに作用する変動応力を実測して、作用応力が疲労強度上、問題のないレベルであることを確認する試験が行われる。
このランナ応力測定は、一般にランナに歪みゲージを取付けて、作用応力を直接に測定する方式が採られている。しかしながら、歪みゲージ及びゲージリードを滑らかなランナ流水面に高速流水に耐えられるように強固に取付けることは、困難であり、特に実働応力測定の必要性が特に大きい高落差高速機ほどゲージが剥がれ易いといった問題があった。更に、歪みゲージによる測定は、ランナの限られた箇所の応力測定であり、ランナ全体の振動状態の情報を得ることはできないといった問題もある。
【０００６】
一方、ランナ振動は、ポンプ水車の運転状態を監視する上で極めて重要な測定ファクターであり、従来は、ランナに近接した静止部にギャップセンサを設置して、このギャップセンサによってこのランナ振動の測定を行い、この測定信号から主要周波数成分の大きさを監視することが行われていた。しかしながら、このような従来の方法では、静止部側からのランナ振動の測定結果が一体、ランナのどのような振動モード現象を表しているのかが曖昧であり、何か異常が発生した場合にも、振動状態が通常と異なるといった情報が得られるに過ぎず、その原因や異常の程度を判断することができないといった問題がある。
【０００７】
他方、ポンプ水車の高落差高速化に伴ってランナの振動現象に関する研究も勧められ、日本機械学会論文集４９巻４３９号の論文「静止側の分布励振源による羽根付回転円板の振動」には、羽根の付いた円板が静止翼列と干渉を起こした時に、円板には特定のスピンする節直径振動モードが励振されることが開示されている。また、この論文には、羽根付き回転円板と見做せるランナにも節直径振動モードが励振される可能性が示され、しかもモードの節直径数は、ランナ羽根枚数とガイドベーン枚数との組合わせによって限定されることが明らかにされた。更に、ランナ振動を静止部側から観測すると、ランナから観測した周波数とは異なる周波数が観測されることも明らかにされ、それまで不明であったポンプ水車の静止部側に現れる特定の振動周波数の意味が解釈されるようになった。
【０００８】
上記の論文によると、ランナとガイドベーンとの翼列干渉により節直径振動が励振される条件は、次式で与えられる。
ｈ・Ｚｇ±ｎ＝ｍ・Ｚｒ            （１）
ここで、Ｚｇはガイドベーン枚数
Ｚｒはランナ羽根枚数
ｈ、ｍ、ｎは、夫々任意の正の整数
なお、ｎの前に付された複号±は、スピン振動モードのスピン方向を示し、記号＋はランナの回転方向と同方向を表し、記号−はランナと逆の回転方向を表している。
【０００９】
式（１）を満足するｈ、ｍ、ｎの組合わせが存在すれば、ランナにはｎ節直径のスピン振動モードが励振され、回転同期周波数をＮとすると、ランナで観測したその振動周波数ｆｒは、次式で表される。
ｆｒ＝ｈ・Ｚｇ・Ｎ                （２）
また、振動を静止部で観測した時の振動周波数ｆｓは、次式で表される。
ｆｓ＝ｍ・Ｚｒ・Ｎ                （３）
また、上記の論文によるとランナから観測した振動Ｘｍは次式（４）で表される。
【数２】
0003756531.tif 000002
この時には、静止部で観測した振動Ｘｍは次式（５）で表される。
【数３】
0003756531.tif 000003
ここで、ω＝２π・ｆｒ
Ω＝２π・Ｎ
ｔは時間
φはランナに固定された座標における円周角
θは静止部に固定された座標における円周角
従って、ランナが単一の節直径振動モードで振動している場合には、静止部で測定した振動から式（３）及び（５）を参照して、Ａｍとｍとを求め、更にｈを仮定することによって式（１）からｎが求まるので、式（４）からランナの振動を算出することが可能になる。
また、高落差ポンプ水車では、ランナが共振することを避けるために、上述のように実機にセンサを設置して振動を検出する代りに、最近では、高落差ポンプ水車の設計の際に、相似模型による実落差模型試験を実施し、実機ランナが定格速度点で共振しないことを模型で確認した上で実機製作を行う方法がとられている。
【００１０】
図１０は横軸に回転速度をとり、縦軸にランナ羽根先端付根に作用する変動応力の片側振幅をとって、ランナの共振特性と使用落差との関係を示したグラフである。縦線は、このランナを２００ｍ級から７００ｍ級までの各落差に適用する場合に水力設計上で必要となる回転速度を示している。このグラフから分かるように、ランナは一般に１次共振ピークＰ１及び２次共振ピークＰ２を持った振動特性であるため、模型試験の結果、定格回転速度点付近に１次または２次共振ピークＰ１、Ｐ２が認められた場合には、定格回転速度がこれらの１次及び２次共振ピークＰ１、Ｐ２から十分に離れるように、ランナのクラウン形状やバンド形状を修正して共振点を移動させていた。
【００１１】
【発明が解決しようとする課題】
ところが、上述した静止部側に振動センサを設置して振動の式からランナの振動モードを算出する方式式は、ランナが単一の節直径振動モードで振動している場合には確かに有効であるが、しかしながら、実際のランナは上記の式（４）で表される複数の節直径振動モードが混在しているため、実際のランナの振動モードを算出することができない。即ち、複数の節直径振動モードが混在している場合には、各振動モードを加算して重ね合わせる必要がある。しかしながら、上記の式（４）には、モードを加算するために必要なモードの空間位相が示されておらず、加算計算を行うことができない。
【００１２】
また、上述した模型試験の結果に基づき定格回転速度を共振点から十分に離す方式は、定速度運転のポンプ水車にとっては有効であるが、しかしながら、可変速運転のポンプ水車の場合には必ずしも万能とはいえない。
この点を詳述すると、近年は、パワーエレクトロニクスの進歩により、大容量の周波数変換装置が開発され、揚水入力可変幅や発電運転効率上昇幅を大きくとれるポンプ水車の可変速運転が実用段階に入っている。このような可変速機では、図１１に示したように中心回転数Ｎ０から±ΔＮの可変速幅の中で運転され、この可変速幅±ΔＮは通常中心回転数Ｎ０の約±１０％に設定される。可変速幅
がこのように大きい可変速機にあっては、模型試験結果に基づくランナのクラウン形状やバンド形状の修正によっては可変速範囲から共振ピークを追い出すことが極めて困難である。
【００１３】
特に、４００ｍ級以上の高落差のポンプ水車では、図１０に示したように運転点が共振点に近付き、５００ｍ級以上になると、一次共振点Ｐ１を越えるため、４００ｍ級以上の高落差可変速ポンプ水車では、図１１に示したようにランナのクラウン形状やバンド形状の修正によってランナの固有振動数を調整して、振動特性を特性Ａから特性ＢまたはＣに変化させても、運転可変速範囲に共振点が入ってしまう。
そこで、本発明の一般的な目的は、ランナの振動応力を低減することができる水車やポンプやポンプ水車等の水力機械を提供することである。
特に、請求項１に記載の発明の目的は、静止部からの振動測定によって、運転中のランナの振動モードを正確に測定することができる水力機械を提供することである。
請求項５及び６に記載の発明の目的は、共振ピークにおいてもランナの変動応力が十分に小さい可変速水力機械を提供することである。
【００１４】
【課題を解決するための手段】
この目的を達成するために請求項１に記載の発明は、複数枚のランナ羽根を有するランナと、このランナの外周側に円形翼列状に配置された複数枚のガイドベーンとを具備する水力機械において、上記ランナの外縁部に近接した静止部に設置されランナ振動の信号を発生する非接触振動センサと、上記非接触振動センサの上記ランナ振動信号から、Ｚｒ・Ｎ（Ｚｒは上記ランナ羽根の枚数であり、Ｎは回転同期周波数である）の周波数成分及びその高次調波成分を抽出してそれぞれの周波数成分信号を発生する周波数分析装置と、上記ランナ円周上の所定の位置が上記静止部の所定の位置に達した時にトリガー信号を発生するトリガー装置と、上記トリガー信号に対する上記抽出された各周波数成分信号の初期位相角を検出する初期位相角検出装置と、上記周波数分析装置の各周波数成分信号の大きさと上記初期位相角検出装置によって検出された各周波数成分信号の初期位相角とに基づきランナ振動モードの変位パターンを算出しそれらを合算する演算装置とを具備することを特徴とするものである。
【００１５】
この構成にあっては、上記ガイドベーンの枚数をＺｇ、回転軸回りの角度座標をランナ回転方向にφ（図２参照、ただし、φはランナ回転方向が正である）、上記周波数成分信号のｍ次成分における片振幅の大きさをＡｍ、時間をｔ、上記初期位相角をΨｍとした時に、上記演算装置は、上記ランナ振動モードの変位パターンＸｍを、
ｋ＝ｍ・Ｚｒ−Ｚｇ
Ｘｍ＝Ａｍ・ｓｉｎ［２π・Ｚｇ・Ｎ・ｔ−ｋ｛φ−Ｚｇ・Ψｍ／ｋ
（Ｚｇ＋ｋ）｝］
から算出し、ランナの実振動モードパターンＸｒを、
【数４】
0003756531.tif 000004
として合算することが好ましい。
更に、上記非接触振動センサは上記ランナの上カバー及び下カバーに夫々設置されていることが望ましい。
また、上記演算装置の出力を判別するパターン判別装置と、このパターン判別装置の出力に基づき水力機械を停止させる制御装置とを更に具備することが望ましい。
【００１６】
【作用】
請求項１に記載の発明の作用は以下の通りである。
非接触振動センサは運転中のランナの振動を検出し、振動の信号を発生する。この振動信号はＺｒ・Ｎの周波数成分及びその高次調波成分を含み、周波数分析装置は、これらのＺｒ・Ｎの周波数成分及びその高次調波成分を抽出して、それぞれの周波数成分信号を発生する。トリガー装置はランナ円周上の所定の位置が静止部の所定の位置に達した時にトリガー信号を発生する。
初期位相角検出装置はトリガー信号に対する上記抽出された各周波数成分信号の初期位相角を検出し、演算装置は周波数分析装置の各周波数成分信号の大きさと初期位相角検出装置によって検出された各周波数成分信号の初期位相角とに基づきランナ振動モードの変位パターンを算出しそれらを合算することによって、ランナの実振動モードパターンを合成する。
【００１７】
次に、請求項１に記載の発明の測定原理を数式を使用して詳述する。
上述の式（４）及び式（５）に±ｎ＝ｋを代入すると、夫々次式（６）及び（７）が得られる。
Ｘｍ＝ＡｍＳＩＮ（ω・ｔ−ｋ・φ）                  （６）
Ｘｍ＝ＡｍＳＩＮ｛（ω＋ｋ・Ω）ｔ−ｋ・θ｝        （７）
この式（７）は、静止部からランナ振動を測定した場合の一般式である。非接触振動センサはトリガー装置からの静止部のθ＝０の位置でトリガー信号によって、トリガーをかけて、実際のランナ振動を観測すると、初期位相角Ψｍを含んだ次式が得られる。
Ｘｍ＝ＡｍＳＩＮ｛（ω＋ｋ・Ω）ｔ＋Ψｍ−ｋ・θ｝    （８）
これは次式（９）に変形される。
0003756531.tif 000005
一方、式（６）が表す振動と式（７）が表す振動は、完全に同一であるので、式（９）に対応する式（６）相当の式は、次式（１０）で表すことができる。
0003756531.tif 000006
この式（１０）は更に次式（１１）に変形することができる。
0003756531.tif 000007
この式（１１）にω＝２π・ｆｒ及びΩ＝２π・Ｎを代入すると、次式（１２）が得られる。
0003756531.tif 000008
この式（１２）におけるｈは、ガイドベーン側の励振源としての水圧分布や流速分布の高調波次数を示すが、これまでの多くのランナ実働応力測定の結果から考察すると、ランナ変動応力の周波数成分にｈ≧２のｈ・Ｚｇ・Ｎ成分が見られたことはない。従って少なくともポンプ水車に関しては、ｈ＝１と見做して差支えない。そこで、ｈ＝１を上式（１２）に代入すると、次式（１３）が得られる。
0003756531.tif 000009
この式（１３）のＺｇ・Ψｍ／ｋ（Ｚｇ＋ｋ）は、ランナに励振されているｍ次節直径振動成分の空間座標における初期位相角を示している。
このように、式（１３）は、ランナに励振される複数の節直径振動モードの各々がランナ上のどこに存在するかを示す空間座標における位相角情報を有しているので、周波数分析装置の各周波数成分に対応する夫々の節直径成分を同一のφ座標の上で加算することができ、実態振動モードＸｒは、次式（１４）から求めることができる。
【数５】
0003756531.tif 000010
なお、式（１３）に示された振動式は、周波数分析装置の各周波数成分信号の大きさと初期位相角検出装置によって検出された各周波数成分信号の初期位相角とに基づき演算装置によって算出される。
【００１９】
【実施例】
以下に本発明による水力機械の実施例を図９及び図１１と同部分には同一符号を付して示した図１乃至図８を参照して説明する。
図１及び図２は本発明の第１の実施例を示したもので、この実施例のポンプ水車は、図２に明示したようにランナ羽根枚数Ｚｒ＝６、ガイドベーン枚数Ｚｇ＝２０である。
ランナ２の外縁部に近接した上カバー６などの静止部に非接触振動センサ１２が、設置されランナ振動の信号を発生する。また、主軸１の外周囲に対向するようにトリガー発生装置１３が配置されている。このトリガー発生装置１３は、例えば主軸１の外周に貼設した信号片とこれを検出するセンサとから構成され、ランナの回転に伴いランナ上の所定位置が静止部の所定位置に対向したことを検出して、トリガー信号Ｔｒを発生する。
【００２０】
非接触振動センサ１２には振動計１４が接続され、この振動計１４は非接触振動センサ１２からの振動信号を増幅する。振動計１４は周波数分析装置１５に接続され、この周波数分析装置１５は振動計１４で増幅された非接触振動センサ１２からの振動信号を受け、この振動信号から、ランナ羽根枚数Ｚｒに回転同期周波数Ｎを乗じたＺｒ・Ｎの周波数成分及びその高次調波成分を抽出する。なお、Ｚｒ・Ｎの基調成分とその高次調波成分とを以下ではｍ・Ｚｒ・Ｎ成分と総称する。
周波数分析装置１５は、抽出したｍ・Ｚｒ・Ｎ成分信号を位相角検出装置１６に送出する。この位相角検出装置１６はトリガー発生装置１３からのトリガー信号Ｔｒと周波数分析装置１５からのｍ・Ｚｒ・Ｎ成分信号とを受け、トリガー信号に対するｍ・Ｚｒ・Ｎ成分の各成分ごとに初期位相角を求める。演算装置１７は、ｍ・Ｚｒ・Ｎ成分の各成分ごとに振動モードを計算して割り付け、割り付けた各振動モードの変位パターンを合成して、ランナの実振動モードパターンを合成計算する。
【００２１】
次に、第１の実施例の作用を説明する。
非接触振動センサ１２は運転中のランナ２の振動を検出し、振動の信号を発生する。周波数分析装置１５は、振動計１４によって増幅された振動信号の周波数を分析し、この振動信号からｍ・Ｚｒ・Ｎ成分信号を抽出する。なお、図３はこの周波数分析装置１５による周波数分析結果を示したグラフであり、ｍ・Ｚｒ・Ｎ成分信号の各々について、横軸に周波数をとり縦軸に振幅をとったグラフ、及びｍ＝１、２、３の成分信号について横軸に時間ｔをとり縦軸に振幅をとったグラフが夫々示されている。
トリガー装置１３はランナ円周上の所定の位置が静止部の所定の位置に達した時にトリガー信号Ｔｒを発生し、初期位相角検出装置１６はトリガー信号Ｔｒに対する各ｍ・Ｚｒ・Ｎ成分信号の初期位相角Ψｍを検出し、演算装置１７は周波数分析装置１５の各ｍ・Ｚｒ・Ｎ成分信号の大きさと初期位相角検出装置１６によって検出された各周波数成分信号の初期位相角Ψｍとに基づき、式（１３）に従ってランナ振動モードの変位パターンを算出しそれらを合算することによって、ランナの実振動モードパターンを合成する。この時、式（１３）の計算に必要なガイドベーン枚数Ｚｇやランナ羽根枚数Ｚｒや回転同期周波数Ｎなどの情報は、予め演算装置１７に入力されている。
【００２２】
図４は、ランナの変動応力分布を求める方法を示したもので、有限要素法などの数値解析手法を用いて、ランナ外縁部振動モード１８Ａと構造解コード１８Ｂとからランナの変動応力分布１８Ｃを計算することができる。
図５は、第１の実施例の変形例を示したもので、上カバー６の非接触振動センサ１２の他に、下カバー７にも非接触振動センサ１９が設置されている。この非接触振動センサ１９の出力信号は、振動計２０を介して周波数分析装置１５に入力される。その他の構成は図１と同一である。
Ｚｒ＝６、Ｚｇ＝２０のポンプ水車におけるランナのクラウン側及びバンド側の両方の振動計測結果の一例を次に説明する。
計測の結果、ｍ・Ｚｒ・Ｎ成分としては、ｍ＝２、３、４、５の各成分が検出された。これらの成分は、ｋ＝ｍ・Ｚｒ−Ｚｇより、それぞれｋ＝−８、−２、４、１０となる。これらの４成分について式（１３）を適用して夫々求めた８節直径モードパターンと２節直径モードパターンと４節直径モードパターンと１０節直径モードパターンが、図６の（ａ）、（ｂ）、（ｃ）、（ｄ）に示されている。また、これらを合算した合成実振動モードパターンが図６の（ｅ）に示されている。この図６（ｅ）の合成実振動モードパターンから分かるように、クラウン３とバンド４は、羽根５の間の部分が振動し、羽根５の部分がほとんど振動しておらず、自然な振動モードが得られている。
【００２３】
図７は、空中においてランナ単体を打撃することによって求めた空中固有振動モードの一例であり、図６の水中でのランナ振動モードとは多少異なるが、基本的なパターンは良く類似しており、運転中のランナ振動モードがランナ単体の空中固有振動モードと基本的に同一であることが検証された。なお、打撃試験を実施すると、数十もの空中固有振動モードが観測され、そのどのモードで振動しているかを特定することは一般に非常に困難である。しかし、本発明による振動モード測定を実施することによって、適確かつ簡単に特定するできることが実証された。
【００２４】
図８は、本発明の第２の実施例を示したもので、演算装置１７の出力にはパターン判別装置２１が接続され、このパターン判別装置２１の出力にはポンプ水車制御装置２２が接続されている。その他の構成は図５と同一である。
パターン判別装置２１は、演算装置１７の出力であるランナの振動モードパターンを識別し、ランナ振動モードが異常かどうかを判別すると共に、モード異常が重大事故か否かを判別する。例えば、振動モードパターンがランナ羽根が伸び縮みするモードである場合には、ランナ羽根とクラウンやバンドとの接合部に割れが入っているので、パターン判別装置２１はポンプ水車の急停止信号をポンプ水車制御装置２２に送出し、これに応じてポンプ水車制御装置２２は、ポンプ水車を急停止すると共に、警報を発生する。
【００３６】
【発明の効果】
以上の説明から明らかなように本発明によれば、運転中のランナの振動を静止部から測定した測定値に基づき、複数の節直径振動モードを重ね合わせてランナの実振動モードパターンを合成計算することができる。
【図面の簡単な説明】
【図１】  本発明による水力機械の第１の実施例を概略的に示した断面図。
【図２】  第１の実施例のポンプ水車の水平断面図。
【図３】  横軸に周波数及び時間をとり、縦軸に振幅をとったグラフ。
【図４】  ランナの応力分析手法を示した説明図。
【図５】  第１の実施例の変形例を概略的に示した断面図。
【図６】  運転中のランナの振動状態を示した説明図。
【図７】  ランナの空中固有振動モードを示した説明図。
【図８】  本発明による水力機械の第２の実施例を概略的に示した断面図。
【図９】  従来のフランシス形ポンプ水車を概略的に示した断面図。
【図１０】  回転速度と変動応力との関係を示したグラフ。
【図１１】  回転速度と変動応力との関係を示したグラフ。
【符号の説明】
２  ランナ
３  クラウン
４  バンド
５  ランナ羽根
８  ガイドベーン
１２  非接触振動センサ
１３  トリガー装置
１５  周波数分析装置
１６  位相角検出装置
１７  演算装置
１９  非接触振動センサ
【図１】
0003756531.tif 000012
【図２】
0003756531.tif 000013
【図３】
0003756531.tif 000014
【図４】
0003756531.tif 000015
【図５】
0003756531.tif 000016
【図６】
0003756531.tif 000017
【図７】
0003756531.tif 000018
【図８】
0003756531.tif 000019
【図９】
0003756531.tif 000020
【図１０】
0003756531.tif 000021
【図１１】
0003756531.tif 000022