(57)【特許請求の範囲】
【請求項１】
加熱炉を含む複数の熱処理炉からなる、金属ストリップの連続熱処理設備での通板の中央速度を制御するための速度制御方法であって、
前記連続熱処理設備の各熱処理炉での設備能力、および、前記金属ストリップの板幅と板厚を含む製造仕様を含む制約条件に基づいて、各熱処理炉の数式モデルを用いて各熱処理炉の理論制約速度を算出し、前記各熱処理炉の理論制約速度の最低速度を求めることにより、前記金属ストリップが現在から所定距離だけ進行する間の制約速度スケジュールであるネック速度スケジュールを作成するネック速度スケジュール作成手順と、
前記ネック速度スケジュール作成手順により作成されたネック速度スケジュールをベースにして、予め用意した複数パターンの速度変更タイミング及び速度変更レートに基づいて速度スケジュールを作成し、それぞれの速度スケジュールについて前記連続熱処理設備に含まれる各熱処理炉を模擬したシミュレータにより板温を予測するシミュレーション手順と、
前記シミュレーション手順により得られた板温および目標板温並びに速度変更レートを要素として含む評価関数を使用して、速度変更タイミング及び速度変更レートを決定する中央速度探索手順とを実行することを特徴とする連続熱処理設備の速度制御方法。
【請求項２】
通板の中央速度の上限値及びネック速度スケジュールに対するバイアス値をオペレータから入力可能としておき、
前記シミュレーション手順では、入力された中央速度の上限値及びネック速度スケジュールに対するバイアス値を考慮した上で、予め用意された複数パターンの速度変更タイミング及び速度変更レートに基づいて速度スケジュールを作成することを特徴とする請求項１に記載の連続熱処理設備の速度制御方法。
【請求項３】
加熱炉を含む複数の熱処理炉からなる、金属ストリップの連続熱処理設備での通板の中央速度を制御するための速度制御装置であって、
前記連続熱処理設備の各熱処理炉での設備能力、および、前記金属ストリップの板幅と板厚を含む製造仕様を含む制約条件に基づいて、各熱処理炉の数式モデルを用いて各熱処理炉の理論制約速度を算出し、前記各熱処理炉の理論制約速度の最低速度を求めることにより、前記金属ストリップが現在から所定距離だけ進行する間のネック速度スケジュールを作成するネック速度スケジュール作成手段と、
前記ネック速度スケジュール作成手段により作成されたネック速度スケジュールをベースにして、予め用意した複数パターンの速度変更タイミング及び速度変更レートに基づいて速度スケジュールを作成し、それぞれの速度スケジュールについて前記連続熱処理設備に含まれる各熱処理炉を模擬したシミュレータにより板温を予測するシミュレーション手段と、
前記シミュレーション手段により得られた板温および目標板温並びに速度変更レートを要素として含む評価関数を使用して、速度変更タイミング及び速度変更レートを決定する中央速度探索手段とを備えたことを特徴とする連続熱処理設備の速度制御装置。
【請求項４】
通板の中央速度の上限値及びネック速度スケジュールに対するバイアス値をオペレータが入力できる入力手段をさらに備え、
前記シミュレーション手段では、入力された中央速度の上限値及びネック速度スケジュールに対するバイアス値を考慮した上で、予め用意された複数パターンの速度変更タイミング及び速度変更レートに基づいて速度スケジュールを作成することを特徴とする請求項３に記載の連続熱処理設備の速度制御装置。
【請求項５】
加熱炉を含む複数の熱処理炉からなる、金属ストリップ連続熱処理設備での通板の中央速度を制御するための処理をコンピュータに実行させるコンピュータプログラムあって、
前記連続熱処理設備の各熱処理炉での設備能力、および、前記金属ストリップの板幅と板厚を含む製造仕様を含む制約条件に基づいて、各熱処理炉の数式モデルを用いて各熱処理炉の理論制約速度を算出し、前記各熱処理炉の理論制約速度の最低速度を求めることにより、前記金属ストリップが現在から所定距離だけ進行する間の制約速度スケジュールであるネック速度スケジュールを作成するネック速度スケジュール作成処理と、
前記ネック速度スケジュール作成処理により作成されたネック速度スケジュールをベースにして、予め用意した複数パターンの速度変更タイミング及び速度変更レートに基づいて速度スケジュールを作成し、それぞれの速度スケジュールについて前記連続熱処理設備に含まれる各熱処理炉を模擬したシミュレータにより板温を予測するシミュレーション処理と、
前記シミュレーション処理により得られた板温および目標板温並びに速度変更レートを要素として含む評価関数を使用して、速度変更タイミング及び速度変更レートを決定する中央速度探索処理とをコンピュータに実行させることを特徴とするコンピュータプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、連続溶融めっき設備や連続焼鈍設備等の連続熱処理設備での通板の中央速度を制御するための速度制御方法、装置、及びコンピュータプログラムに関する。
【背景技術】
【０００２】
例えば鉄鋼業等の金属製造プラントで、金属ストリップの連続溶融めっき設備や連続焼鈍設備等の連続熱処理設備において、加熱炉入側から合金化炉出側までの間、または加熱炉入側から一連の熱処理炉出側までの間の通板の中央速度を制御する場合、設備能力や製造仕様といった制約条件に基づいて制約速度スケジュール（本明細書では「ネック速度スケジュール」と称する）を決定する必要がある。制約条件が複数存在する場合は、それら制約条件をすべて満たすように、各制約条件から求められる理論制約速度の最低速度をネック速度スケジュールとして抽出する。
【０００３】
ネック速度スケジュールが決定されたならば、それをベースにして適宜な加減速レートで速度変更させるのであるが、固定的に決められている加減速レートで速度変更した場合、板温が上下限範囲を外れる等して操業を不安定にする結果をまねくおそれもある。
【０００４】
現在のところ、ネック速度スケジュールを画面表示する等のガイダンスレベルのシステム構築はなされているものの、最終的な速度決定はオペレータの判断によるものとされているのが実情である。
【０００５】
連続熱処理設備での通板速度を制御する技術として、例えば特許文献１には、板温測定値に基づいてライン速度の加減速レートを自動設定することにより、連続焼鈍炉に通板される鋼帯に板温外れを生じさせることなく、適正な加減速レートで鋼帯のライン速度を変更するようにしたライン速度の加減速レート自動設定システムが開示されている。
【０００６】
また、特許文献２には、炉負荷燃焼、熱処理方法による制約に基づいて中央速度の制約速度を算出する工程と、加減速率の決定を行う工程とを行うようにした連続焼鈍ラインの中央速度制御方法が開示されている。
【０００７】
また、特許文献３には、推定された板温変化量とストリップの目標板温との偏差が許容範囲外である場合には、加熱炉の炉温制御及びストリップの通板速度の修正制御を行うようにした連続焼鈍炉における板温制御方法が開示されている。
【０００８】
【特許文献１】特開２００３−２２６９１１号公報
【特許文献２】特開平６−３３０１８２号公報
【特許文献３】特開平２−２５８９３３号公報
【発明の開示】
【発明が解決しようとする課題】
【０００９】
しかしながら、特許文献１に開示されたライン速度の加減速レート自動設定システムでは、制約条件に基づくネック速度スケジュールが考慮されていないという問題がある。
【００１０】
また、特許文献２に開示された連続焼鈍ラインの中央速度制御方法では、加減速率を決定するに際して、実験から知見された各熱処理炉の板温変動しない限界速度（律速線）を用いるが、あくまでも実験に基づくものであるため、厳密性に欠けるという問題がある。
【００１１】
また、特許文献３に開示された板温制御方法では、推定された板温変化量とストリップの目標板温との偏差が許容範囲外である場合に、加熱炉の炉温制御及びストリップの通板速度の修正制御を行うものであり、板温変化量と速度変更レートを同時に計算して最適操業条件を決定する技術ではない。
【００１２】
本発明は上記のような点に鑑みてなされたものであり、制約条件に基づくネック速度スケジュールを考慮しつつ板温が上下限範囲を外れることなく、厳密性を保ちながら、板温変化量と速度変更レートを同時に計算して可及的に良好な操業条件を決定することを目的とする。
【課題を解決するための手段】
【００１３】
本発明による連続熱処理設備の速度制御方法は、加熱炉を含む複数の熱処理炉からなる、金属ストリップの連続熱処理設備での通板の中央速度を制御するための速度制御方法であって、前記連続熱処理設備の各熱処理炉での設備能力、および、前記金属ストリップの板幅と板厚を含む製造仕様を含む制約条件に基づいて、各熱処理炉の数式モデルを用いて各熱処理炉の理論制約速度を算出し、前記各熱処理炉の理論制約速度の最低速度を求めることにより、前記金属ストリップが現在から所定距離だけ進行する間の制約速度スケジュールであるネック速度スケジュールを作成するネック速度スケジュール作成手順と、前記ネック速度スケジュール作成手順により作成されたネック速度スケジュールをベースにして、予め用意した複数パターンの速度変更タイミング及び速度変更レートに基づいて速度スケジュールを作成し、それぞれの速度スケジュールについて前記連続熱処理設備に含まれる各熱処理炉を模擬したシミュレータにより板温を予測するシミュレーション手順と、前記シミュレーション手順により得られた板温および目標板温並びに速度変更レートを要素として含む評価関数を使用して、速度変更タイミング及び速度変更レートを決定する中央速度探索手順とを実行する点に特徴を有する。
また、本発明による連続熱処理設備の速度制御方法の他の特徴とするところは、通板の中央速度の上限値及びネック速度スケジュールに対するバイアス値をオペレータから入力可能としておき、前記シミュレーション手順では、入力された中央速度の上限値及びネック速度スケジュールに対するバイアス値を考慮した上で、予め用意された複数パターンの速度変更タイミング及び速度変更レートに基づいて速度スケジュールを作成する点にある。
本発明による連続熱処理設備の速度制御装置は、加熱炉を含む複数の熱処理炉からなる、金属ストリップの連続熱処理設備での通板の中央速度を制御するための速度制御装置であって、前記連続熱処理設備の各熱処理炉での設備能力、および、前記金属ストリップの板幅と板厚を含む製造仕様を含む制約条件に基づいて、各熱処理炉の数式モデルを用いて各熱処理炉の理論制約速度を算出し、前記各熱処理炉の理論制約速度の最低速度を求めることにより、前記金属ストリップが現在から所定距離だけ進行する間のネック速度スケジュールを作成するネック速度スケジュール作成手段と、前記ネック速度スケジュール作成手段により作成されたネック速度スケジュールをベースにして、予め用意した複数パターンの速度変更タイミング及び速度変更レートに基づいて速度スケジュールを作成し、それぞれの速度スケジュールについて前記連続熱処理設備に含まれる各熱処理炉を模擬したシミュレータにより板温を予測するシミュレーション手段と、前記シミュレーション手段により得られた板温および目標板温並びに速度変更レートを要素として含む評価関数を使用して、速度変更タイミング及び速度変更レートを決定する中央速度探索手段とを備えた点に特徴を有する。
また、本発明による連続熱処理設備の速度制御装置の他の特徴とするところは、通板の中央速度の上限値及びネック速度スケジュールに対するバイアス値をオペレータが入力できる入力手段をさらに備え、前記シミュレーション手段では、入力された中央速度の上限値及びネック速度スケジュールに対するバイアス値を考慮した上で、予め用意された複数パターンの速度変更タイミング及び速度変更レートに基づいて速度スケジュールを作成する点にある。
本発明によるコンピュータプログラムは、加熱炉を含む複数の熱処理炉からなる、金属ストリップ連続熱処理設備での通板の中央速度を制御するための処理をコンピュータに実行させるコンピュータプログラムあって、前記連続熱処理設備の各熱処理炉での設備能力、および、前記金属ストリップの板幅と板厚を含む製造仕様を含む制約条件に基づいて、各熱処理炉の数式モデルを用いて各熱処理炉の理論制約速度を算出し、前記各熱処理炉の理論制約速度の最低速度を求めることにより、前記金属ストリップが現在から所定距離だけ進行する間の制約速度スケジュールであるネック速度スケジュールを作成するネック速度スケジュール作成処理と、前記ネック速度スケジュール作成処理により作成されたネック速度スケジュールをベースにして、予め用意した複数パターンの速度変更タイミング及び速度変更レートに基づいて速度スケジュールを作成し、それぞれの速度スケジュールについて前記連続熱処理設備に含まれる各熱処理炉を模擬したシミュレータにより板温を予測するシミュレーション処理と、前記シミュレーション処理により得られた板温および目標板温並びに速度変更レートを要素として含む評価関数を使用して、速度変更タイミング及び速度変更レートを決定する中央速度探索処理とをコンピュータに実行させる点に特徴を有する。
【発明の効果】
【００１４】
本発明によれば、制約条件に基づくネック速度スケジュールを考慮しつつ、熱処理炉を模擬したシミュレータによりシミュレーションを行うことで厳密性を保ちながら、板温変化量と速度変更レートを同時に計算して可及的に良好な操業条件を決定することができる。
【発明を実施するための最良の形態】
【００１５】
以下、添付図面を参照して、本発明の好適な実施形態について説明する。図２に概要のみを示すが、鋼材等の金属ストリップ（ストリップ）１００が加熱炉２００、均熱炉３００、及び合金化炉４００に連続的に通板される連続熱処理設備の例である。
【００１６】
図１に、本発明を適用した連続熱処理設備の速度制御装置の構成を示す。本実施形態の速度制御装置は、図２に示した連続熱処理設備での通板の中央速度を制御するものである。
【００１７】
１は入力部であり、通板の中央速度の上限値及びネック速度スケジュールに対するバイアス値をオペレータが入力するためのものである。図３に、入力部１の一例としてオペレータ介入操作盤を示す。ここでは、任意の上限値［ｍｐｍ］、及び、０〜−１００［ｍｐｍ］の範囲で任意のバイアス値を入力することができる。
【００１８】
２はネック速度スケジュール作成部であり、連続熱処理設備の各熱処理炉２００、３００、４００での制約条件に基づいて、現在からストリップ１００がｘ［ｍ］進行する間のネック速度スケジュールを作成する。
【００１９】
３は加熱炉シミュレータであり、ネック速度スケジュール作成部２により作成されたネック速度スケジュールをベースにして、入力部１から入力された中央速度の上限値及びネック速度スケジュールに対するバイアス値を考慮した上で、予め用意された複数パターンの速度変更タイミング及び速度変更レートに基づいて速度スケジュールを作成し、それぞれの速度スケジュールについて加熱炉２００の出側板温の応答性をシミュレーションにより予測する。
【００２０】
４は中央速度探索部であり、加熱炉２００の出側板温の応答性及び速度変更レートを要素として含む評価関数を使用して、最適と考えられる速度変更タイミング及び速度変更レートを決定する。そして、決定された速度変更タイミング及び速度変更レートに基づいて連続熱処理設備での通板の速度制御を実行する。
【００２１】
５は表示装置であり、例えばネック速度スケジュール作成部２により作成されるネック速度スケジュールや、加熱炉シミュレータ３により作成される複数パターンの速度スケジュールを可視的に画像表示する。これにより、オペレータに対する操業ガイダンス機能の向上が実現される。
【００２２】
図４は、本実施形態における連続熱処理設備の速度制御の処理動作を説明するためのフローチャートである。
【００２３】
まず、ステップＳ１において、ネック速度スケジュール作成部２は、加熱炉２００、均熱炉３００、及び合金化炉４００での理論制約速度を算出する。
【００２４】
例えば加熱炉２００での設備仕様上の加熱能力に基づく理論制約速度ＶＣHFは、下式（１）のように適宜な数式モデルｆ1を用いて算出される。
ＶＣHF＝ｆ1（ＴＰＨ，ρ，ＷＤ，ＴＨ，ＶＲ）・・・（１）
ただし、ＴＰＨ：Ｔ／Ｈ［ｔｏｎ／ｈｏｕｒ］、ρ：比重［ｔｏｎ／ｍ3］、ＷＤ：板幅［ｍ］、ＴＨ：板厚［ｍ］、ＶＲ：速度達成率［−］である。
【００２５】
また、均熱炉３００での設備仕様上の均熱時間に基づく理論制約速度ＶＣRSFは、下式（２）のように適宜な数式モデルｆ2を用いて算出される。
ＶＣRSF＝ｆ2（ＰＬRSF，ＴｒｅｆRSF）・・・（２）
ただし、ＰＬRSF：均熱炉パス長［ｍ］、ＴｒｅｆRSF：均熱時間［ｓｅｃ］である。
【００２６】
また、合金化炉４００での設備仕様上の合金化時間に基づく理論制約速度ＶＣHOは、下式（３）のように適宜な数式モデルｆ3を用いて算出される。
ＶＣHO＝ｆ3（ＰＬHO，ＴｒｅｆHO）・・・（３）
ただし、ＰＬHO：合金化炉パス長［ｍ］、ＴｒｅｆHO：合金化時間［ｓｅｃ］である。
【００２７】
次に、ステップＳ２において、ネック速度スケジュール作成部２は、上記ステップＳ１で算出された各熱処理炉２００、３００、４００での理論制約速度ＶＣHF、ＶＣRSF、ＶＣHOに基づいて、現在からストリップ１００がｘ［ｍ］進行する間のネック速度スケジュールを作成する。
【００２８】
具体的には、図５（ａ）に示すように、各熱処理炉２００、３００、４００での理論制約速度ＶＣHF、ＶＣRSF、ＶＣHOをすべて満たすように、理論制約速度ＶＣHF、ＶＣRSF、ＶＣHOの最低速度をネック速度スケジュールＶＣminとして抽出する。
【００２９】
次に、ステップＳ３において、加熱炉シミュレータ３は、上記ステップＳ２で作成されたネック速度スケジュールＶＣminをベースにして、入力部１から入力された中央速度の上限値及びネック速度スケジュールに対するバイアス値を考慮した上で、予め用意された複数パターンの速度変更タイミング及び速度変更レートに基づいて速度スケジュールＶを作成し、それぞれの速度スケジュールＶについて加熱炉２００の出側板温の応答性をシミュレーションにより予測する。
【００３０】
具体的には、図５（ｂ）に示すように、上記ステップＳ２で作成されたネック速度スケジュールＶＣminを、入力部１から入力されたバイアス値分だけ下げるとともに、入力部１から入力された上限値を超える部分をカットする（図中点線を参照）。そして、予め用意された複数パターンの速度変更タイミング及び速度変更レートに基づいて、複数パターンの速度スケジュールＶを作成する。図５（ｂ）には、速度変更レートを変更している様子の一例を示す。
【００３１】
本実施形態では、加熱炉シミュレータ３により予測する加熱炉２００の出側板温の応答性として、図６に示すように、複数パターンの速度スケジュールＶそれぞれについて、加熱炉２００の出側板温の目標板温に対する誤差予測値ΔＴＳ（ｔ）［℃］を求める。加熱炉２００の出側板温の目標板温は、不図示のプロセスコンピュータにより生産スケジュールに従って決定されている。
【００３２】
図７には加熱炉シミュレータ３の概要を示す。本実施形態の加熱炉シミュレータ３では、下式の１〜５で表わされるラジアントチューブ温度モデル７０１、炉温モデル７０２、ハースロール温度モデル７０３、板温モデル７０４、炉壁温度モデル７０５が構築されており、これら数式モデルを用いて加熱炉２００での熱処理を模擬演算し、複数パターンの速度スケジュールＶそれぞれについて加熱炉２００の出側板温の誤差予測値ΔＴＳ（ｔ）［℃］を求める。
【００３３】
【数１】
0004833614.tif 000002
【００３４】
なお、本実施形態では、熱処理炉２００、３００、４００のうち加熱炉２００の板温変化が操業上で重要であることから、加熱炉２００を模擬したシミュレータを構築するようにしたが、例えば各熱処理炉２００、３００、４００についてシミュレータを構築するようにしてもよい。
【００３５】
図８は、ステップＳ３において複数パターンの速度変更タイミング及び速度変更レートに基づいて速度スケジュールＶを作成し、それぞれの速度スケジュールＶについて加熱炉２００の出側板温の応答性をシミュレーションにより予測するための処理動作を詳細に説明するためのフローチャートである。
【００３６】
例えば、速度変更タイミングのパターンα（ｉ）がｉ＝０〜Ｎまで、速度変更レートのパターンβ（ｊ）がｊ＝０〜Ｍまで用意されているものとする。
【００３７】
最初は、ｉ＝０から始め（ステップＳ３０１）、速度変更タイミングパターンα（ｉ）を設定し（ステップ３０２）、ｊ＝０の速度変更レートパターンβ（ｉ）を設定してΣΔ
ＴＳ（ｔ）2を求めてデータベースへ格納する（ステップＳ３０３〜Ｓ３０５）。
【００３８】
次に、ｊに１を加算して（ステップＳ３０６）、ｊが最大となるまで、すなわちｊ＝ＭとなるまでステップＳ３０２〜Ｓ３０５を繰り返す（ステップＳ３０７）。
【００３９】
その後、ｉに１を加算して（ステップＳ３０８）、ｉが最大となるまで、すなわちｉ＝ＮとなるまでステップＳ３０２〜Ｓ３０８を繰り返す（ステップＳ３０９）。
【００４０】
以上述べた処理により、複数パターンの速度変更タイミング（ｉ＝０〜Ｎ）及び速度変更レート（ｊ＝０〜Ｍ）に基づいて速度スケジュールＶを作成し、それぞれの速度スケジュールＶについてΣΔＴＳ（ｔ）2を取得することができる。
【００４１】
次に、ステップＳ４において、中央速度探索部４は、加熱炉２００の出側板温の応答性（誤差予測値ΔＴＳ（ｔ））及び速度変更レートを要素として含む評価関数Ｊsを使用して、最適と考えられる速度変更タイミング及び速度変更レートを決定する。
【００４２】
具体的には、上記ステップＳ３で算出されたΣΔＴＳ（ｔ）2と、速度変更レートＶra
te［ｍｐｍ／ｓｅｃ］と、適宜な重みｗ［−］とにより、下式（４）により構成される評価関数Ｊsが最小の値をとる速度変更タイミング及び速度変更レートを決定する。
Ｊs＝ΣΔＴＳ（ｔ）2＋ｗ・（１／Ｖrate）・・・（４）
【００４３】
その後、ステップＳ４により決定された速度変更タイミング及び速度変更レートに基づいて連続熱処理設備での通板の速度制御を実行する。
【００４４】
以上述べたように、制約条件に基づくネック速度を考慮しつつ、加熱炉２００を模擬したシミュレータによりシミュレーションを行うことで厳密性を保ちながら、板温変化量と通板の中央速度との相関の最適化を図ることができる。
【００４５】
なお、以上述べた実施形態の速度制御装置は具体的にはコンピュータシステム或いは装置により構成されるものである。したがって、上述した機能を実現するソフトウェアのプログラムコードを記録した記憶媒体をシステム或いは装置に供給し、そのシステム或いは装置のコンピュータ（又はＣＰＵやＭＰＵ）が記憶媒体に格納されたプログラムコードを読み出し実行することによっても、達成されることは言うまでもない。
【００４６】
この場合、記憶媒体から読み出されたプログラムコード自体が上述した実施形態の機能を実現することになり、プログラムコード自体及びそのプログラムコードを記憶した記憶媒体は本発明を構成することになる。プログラムコードを供給するための記憶媒体としては、例えば、フレキシブルディスク、ハードディスク、光ディスク、光磁気ディスク、ＣＤ−ＲＯＭ、ＣＤ−Ｒ、磁気テープ、不揮発性のメモリカード、ＲＯＭ等を用いることができる。
【図面の簡単な説明】
【００４７】
【図１】本実施形態の連続熱処理設備の中央速度制御装置の構成を示す図である。
【図２】本実施形態の連続熱処理設備の概要を説明するための図である。
【図３】オペレータ介入操作盤を示す図である。
【図４】連続熱処理設備の速度制御の処理動作を説明するためのフローチャートである。
【図５】進行距離と速度との関係を示す図であり、（ａ）がネック速度スケジュールを作成する様子を示す図であり、（ｂ）が速度スケジュールを作成する様子を示す図である。
【図６】時間と板温との関係を示す図である。
【図７】加熱炉シミュレータの概要を示す図である。
【図８】複数パターンの速度変更タイミング及び速度変更レートに基づいて速度スケジュールを作成し、それぞれの速度スケジュールについて加熱炉出側板温の応答性をシミュレーションにより予測するための処理動作を説明するためのフローチャートである
【符号の説明】
【００４８】
１      入力部
２      ネック速度スケジュール作成部
３      加熱炉シミュレータ
４      中央速度探索部
５      表示装置
１００  金属ストリップ
２００  加熱炉
３００  均熱炉
４００  合金化炉
【図１】
0004833614.tif 000003
【図２】
0004833614.tif 000004
【図３】
0004833614.tif 000005
【図４】
0004833614.tif 000006
【図５】
0004833614.tif 000007
【図６】
0004833614.tif 000008
【図７】
0004833614.tif 000009
【図８】
0004833614.tif 000010