(57)【要約】【課題】耐用年数を確保しながら発電量の減少を極力抑えること。【解決手段】風車の荷重データを用いて、所定期間における風車の疲労損傷度を評価するための疲労等価荷重を算出する疲労等価荷重算出部１１と、疲労等価荷重と風車の耐用年数に応じて決定される基準荷重とを比較し、その差分が所定の閾値を超えていた場合に、現在採用されている運転制限の条件設定値を該差分に応じて更新する設定値更新部１２とを具備する風車の運転制限調整装置１０を提供する。【選択図】図２
2010048239.tif 000002
【特許請求の範囲】
【請求項１】
風速が所定の条件設定値を超えた場合に、風車の運転を停止する運転制限を行う風車に適用され、該運転制限の条件設定値を更新する運転制限調整装置であって、
風車の荷重データを用いて、所定期間における風車の疲労損傷度を評価するための疲労等価荷重を算出する疲労等価荷重算出手段と、
前記疲労等価荷重と風車の耐用年数から決定される基準荷重とを比較し、その差分が所定の閾値を超えていた場合に、現在採用されている運転制限の条件設定値を該差分に応じて更新する設定値更新手段と
を具備する風車の運転制限調整装置。
【請求項２】
前記運転制限における条件設定値が予め設定された風向毎に設定されている場合に、前記疲労等価荷重算出手段は、該風向毎に前記疲労等価荷重を算出し、
前記設定値更新手段は、該疲労等価荷重と該基準荷重とを風向毎に比較し、その差分に応じて現在採用されている運転制限の該条件設定値を風向毎に更新する請求項１に記載の風車の運転制限調整装置。
【請求項３】
前記荷重データは時系列荷重データであり、
前記疲労等価荷重は、該所定期間よりも短い所定の分割時間毎に、該分割時間内における前記時系列荷重データを統計的に処理することで、該分割時間における疲労等価荷重を求め、この疲労等価荷重を前記所定期間において積分することで、該所定期間における疲労等価荷重を算出する請求項１または請求項２に記載の風車の運転制限調整装置。
【請求項４】
前記疲労等価荷重は、前記分割時間毎にレインフローカウント法を採用した以下の式を用いて算出される請求項３に記載の風車の運転制限調整装置。
【数１】
2010048239.tif 000006
ここで、leq,kは該所定期間における処理開始からk番目の前記分割時間における疲労等価荷重、liは疲労荷重スペクトルのi番目の荷重レンジ、niは疲労荷重スペクトルのi番目の荷重レンジのサイクル数、neqは、前記第２の期間における基準等価サイクル数、mは関連材料のS-N曲線の勾配である。
【請求項５】
複数の風況条件と推定荷重とを関連付けたデータベースを有し、
前記運転制限の初期条件設定値は、
風車の建設地域における風況を解析し、
該風況に対応する推定荷重を該データベースから取得し、
取得した推定荷重と前記基準荷重との比較結果に基づいて決定される請求項１から請求項４のいずれかに記載の風車の運転制限調整装置。
【請求項６】
風速が所定の条件設定値を超えた場合に、風車の運転を停止する運転制限を行う風車に適用され、該運転制限の条件設定値を更新する運転制限調整方法であって、
風車の荷重データを用いて、所定期間における風車の疲労損傷度を評価するための疲労等価荷重を算出する過程と、
前記疲労等価荷重と風車の耐用年数から決定される基準荷重とを比較し、その差分が所定の閾値を超えていた場合に、現在採用されている運転制限の条件設定値を該差分に応じて更新する過程と
を具備する風車の運転制限調整方法。
【請求項７】
風速が所定の条件設定値を超えた場合に、風車の運転を停止する運転制限を行う風車に適用され、該運転制限の条件設定値を更新する運転制限調整をコンピュータに実現させるための風車の運転制限調整プログラムであって、
風車の荷重データを用いて、所定期間における風車の疲労損傷度を評価するための疲労等価荷重を算出する処理と、
前記疲労等価荷重と風車の耐用年数から決定される基準荷重とを比較し、その差分が所定の閾値を超えていた場合に、現在採用されている運転制限の条件設定値を該差分に応じて更新する処理と
を有する風車の運転制限調整プログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、風車の運転制限調整装置及び方法並びにプログラムに関するものである。
【背景技術】
【０００２】
ウィンドファーム内に設置された風車は、風車の配置と風況条件次第で、風車間の相互干渉による疲労荷重が増加し、設計荷重を超過する場合がある。対策を講じずにそのまま運転した場合、設計時に考慮された耐用年数（典型的な風車では、例えば、２０年）を待たずして、風車の故障若しくは事故が生じる可能性がある。そのため、風車干渉の影響が大きい風向において、通常のカットアウト風速（典型的な風車では２５ｍ／ｓ）以下の風速で運転停止（シャットダウン）することにより運転時間を制限し、疲労荷重を低減させる対策が一般的にとられる。これを運転制限と呼ぶ。運転制限が行われると運転時間が減少するため、当然の如く、発電量も減少してしまう。
【０００３】
そこで、発電減少量を最小限にとどめるために、事前検討を行い、風車を停車させる運転制限の条件設定値を極力適切な値に決定する方法が提案されている（例えば、特許文献１参照）。このように、事前検討を行うことで、過剰な運転制限を避け、所望の発電量を確保することが期待できる。
【非特許文献１】”Recent Developments in Turbine Site-Suitability Assessment”，T.R. Camp and D.V. Witcher, 2003 EWEC AbstractBook, European Wind Energy Conference and Exhibition, 16-19 June 2003 in Madrid, Spain
【発明の開示】
【発明が解決しようとする課題】
【０００４】
しかしながら、事前検討で予測した風況等と実際の風況等とが異なることも多く、この場合には、事前検討において決定した運転制限の条件設定値が、現実には最適な値に設定されていないおそれがあり、例えば、発電量確保の点からみて過剰、あるいは、疲労荷重低減の点からみて過小となっており、所望の発電量を得られないという問題があった。
【０００５】
本発明は、上記問題を解決するためになされたもので、耐用年数を確保しながら発電量の減少を極力抑えることのできる風車の運転制限調整装置及び方法並びにプログラムを提供することを目的とする。
【課題を解決するための手段】
【０００６】
上記課題を解決するために、本発明は以下の手段を採用する。
本発明は、風速が所定の条件設定値を超えた場合に、風車の運転を停止する運転制限を行う風車に適用され、該運転制限の条件設定値を更新する運転制限調整装置であって、風車の荷重データを用いて、所定期間における風車の疲労損傷度を評価するための疲労等価荷重を算出する疲労等価荷重算出手段と、前記疲労等価荷重と風車の耐用年数から決定される基準荷重とを比較し、その差分が所定の閾値を超えていた場合に、現在採用されている運転制限の条件設定値を該差分に応じて更新する設定値更新手段とを具備する風車の運転制限調整装置を提供する。
【０００７】
このような構成によれば、実際に風車を建設した後、各風車の荷重を定期的に計測し、この荷重から所定期間における疲労等価荷重を算出し、この疲労等価荷重を風車の耐用年数から決定される基準荷重と比較することで、現在採用されている運転制限の条件設定値の見直しを行うので、風車の建設後において、その疲労損傷度に応じて運転制限の条件設定値を適切な値に調整することが可能となる。これにより、運転制限によって風車が停止される頻度を低減させることができ、耐久年数を確保しながら、運転制限による発電量の低下を抑制することが可能となる。
【０００８】
上記風車の運転制限調整装置において、前記運転制限における条件設定値が予め設定された風向毎に設定されている場合に、前記疲労等価荷重算出手段は、該風向毎に前記疲労等価荷重を算出し、前記設定値更新手段は、該疲労等価荷重と該基準荷重とを風向毎に比較し、その差分に応じて現在採用されている運転制限の該条件設定値を風向毎に更新することとしてもよい。
【０００９】
このような構成によれば、風向毎に運転制限の条件設定値を見直すことができ、運転制限の調整精度を向上させることが可能となる。
【００１０】
上記風車の運転制限調整装置において、前記荷重データは時系列荷重データであり、前記疲労等価荷重は、該所定期間よりも短い所定の分割時間毎に、該分割時間内における前記時系列荷重データを統計的に処理することで、該分割時間における疲労等価荷重を求め、この疲労等価荷重を前記所定期間において積分することで、該所定期間における疲労等価荷重を算出することとしてもよい。
【００１１】
このように、疲労等価荷重は、例えば、分割時間毎の荷重データを統計的に処理して求めた疲労等価荷重を所定期間において積分することで求められる。
【００１２】
上記風車の運転制限調整装置において、例えば、前記疲労等価荷重は、前記分割時間毎にレインフローカウント法を採用した以下の式を用いて算出される。
【数１】
2010048239.tif 000003
ここで、leq,kは該所定期間における処理開始からk番目の前記分割時間における疲労等価荷重、liは疲労荷重スペクトルのi番目の荷重レンジ、niは疲労荷重スペクトルのi番目の荷重レンジのサイクル数、neqは、前記第２の期間における基準等価サイクル数、mは関連材料のS-N曲線の勾配である。
【００１３】
上記風車の運転制限調整装置は、複数の風況条件と推定荷重とを関連付けたデータベースを有し、前記運転制限の初期条件設定値は、風車の建設地域における風況を解析し、該風況に対応する推定荷重を該データベースから取得し、取得した推定荷重と前記基準荷重との比較結果に基づいて決定されることとしてもよい。
【００１４】
このように、複数の風況条件と推定荷重とを関連付けたデータベースを予め用意しておくことで、事前検討において運転制限の初期条件設定値を決定する場合に、推定した風況に最も近い風況条件に対応する推定荷重をデータベースから取得し、この推定荷重に基づいて風車の疲労損傷度を予測することが可能となる。これにより、容易に、かつ、迅速に各風車の疲労損傷度を推定することができる。
【００１５】
本発明は、風速が所定の条件設定値を超えた場合に、風車の運転を停止する運転制限を行う風車に適用され、該運転制限の条件設定値を更新する運転制限調整方法であって、風車の荷重データを用いて、所定期間における風車の疲労損傷度を評価するための疲労等価荷重を算出する過程と、前記疲労等価荷重と風車の耐用年数から決定される基準荷重とを比較し、その差分が所定の閾値を超えていた場合に、現在採用されている運転制限の条件設定値を該差分に応じて更新する過程とを具備する風車の運転制限調整方法を提供する。
【００１６】
本発明は、風速が所定の条件設定値を超えた場合に、風車の運転を停止する運転制限を行う風車に適用され、該運転制限の条件設定値を更新する運転制限調整をコンピュータに実現させるための風車の運転制限調整プログラムであって、風車の荷重データを用いて、所定期間における風車の疲労損傷度を評価するための疲労等価荷重を算出する処理と、前記疲労等価荷重と風車の耐用年数から決定される基準荷重とを比較し、その差分が所定の閾値を超えていた場合に、現在採用されている運転制限の条件設定値を該差分に応じて更新する処理とを有する風車の運転制限調整プログラムを提供する。
【発明の効果】
【００１７】
本発明によれば、耐用年数を確保しながら発電量の減少を極力抑えることができるという効果を奏する。
【発明を実施するための最良の形態】
【００１８】
以下に、本発明に係る風車の運転制限調整装置及び方法並びにプログラムの一実施形態について、図面を参照して説明する。
図１は、本実施形態に係る風車の運転制限調整装置の概略構成を示したブロック図である。図１に示すように、本実施形態に係る風車の運転制限調整装置１０は、コンピュータシステム（計算機システム）であり、ＣＰＵ（中央演算処理装置）１、ＲＡＭ（Random Access Memory）等の主記憶装置２、ＨＤＤ（Hard Disk Drive）等の補助記憶装置３、キーボードやマウスなどの入力装置４、及びモニタやプリンタなどの出力装置５、外部の機器と通信を行うことにより情報の授受を行う通信装置６などで構成されている。
補助記憶装置３には、各種プログラム（例えば、風車の運転制限調整プログラム）が格納されており、ＣＰＵ１が補助記憶装置３から主記憶装置２にプログラムを読み出し、実行することにより種々の処理を実現させる。
【００１９】
この風車の運転制限調整装置１０は、各風車に設けられていてもよいし、ウィンドファームに建設されている全て或いは一部の風車の制御を一元的に管理・監視する制御室に設けられていても良い。
【００２０】
図２は、風車の運転制限調整装置１０が備える機能を展開して示した機能ブロック図である。図２に示されるように、風車の運転制限調整装置１０は、風速が所定の条件設定値を超えた場合に、風車の運転を停止する運転制限を行う風車に適用され、該運転制限の条件設定値を更新する装置であって、疲労等価荷重算出部（疲労等価荷重算出手段）１１と、設定値更新部（設定値更新手段）１２と、設定値記憶部１３とを主な構成として備えている。
【００２１】
疲労等価荷重算出部１１は、風車に取り付けられた荷重センサから出力される荷重データを取得し、この荷重データを用いて所定期間における風車の疲労損傷度を評価するための疲労等価荷重を算出する。
【００２２】
設定値更新部１２は、疲労等価荷重算出部１１によって算出された疲労等価荷重と運転制限の該初期条件設定値を決定する際に用いられた推定疲労荷重とを比較し、その差分が所定の閾値を超えていた場合に、該差分に応じて、設定値記憶部１３に格納されている現在の運転制限の条件設定値を更新する。
【００２３】
本実施形態に係る風車の運転制限調整装置１０は、事前検討によって求められた風車の運転制限の初期条件設定値を風車の疲労損傷度合に応じて適切な値に更新していくものである。
以下、風車の運転制限の初期条件設定値を決定する事前検討における処理手順を説明した後、本実施形態に係る風車の運転制限調整装置１０により実現される処理手順について図３を参照して説明する。ここでは、複数の風車が建設されるウィンドファームを前提に説明する。
【００２４】
〔事前検討〕
まず、事前検討では、ウィンドファームに建設が予定されている風車毎に、風況条件を計算する風況解析が行われる（図３のステップＳＡ１）。風況解析は、例えば、実際の計測値または数値計算、若しくはこれらを組み合わせることで行われる。風況解析では、例えば、風速風向出現頻度分布、乱れ度分布、吹上角、ウィンドシア、風車相互干渉による乱れ度分布（等価乱れ度）等が求められる。ここで、ウィンドシアとは、風車の鉛直方向における風速の分布である。
【００２５】
次に、上記風況解析で求められた値を入力情報として用いて、各風車において所定の風向毎に所定期間（例えば、１年間）にわたる疲労荷重が推定される（図３のステップＳＡ２）。ここで、推定された疲労荷重を推定荷重という。また、推定荷重は、全方位を所定の分割数で分割した方位（例えば１６方位）においてそれぞれ算出される。
推定荷重の算出方法は、例えば、以下に示す３つの方法を任意に選択して採用することが可能である。
〔方法１〕
上記風況解析の解析結果を入力として、建設が予定されている全ての風車に対して一台毎に時系列応答計算を行うことで推定荷重を算出する。この方法では、風車台数が１００台を超えるような大規模なウィンドファームになると、膨大な処理時間を要するという欠点がある。
【００２６】
〔方法２〕
ウィンドファーム内に建設予定の風車の中で、荷重が最も大きいと推定できる風車、換言すると、風車干渉の影響が最も大きい風車を経験的に選出し、この風車の推定荷重を上記方法１を用いて算出するとともに、算出した推定荷重を他の全ての風車にも適用する。
なお、この方法を用いた場合には、経験によって最も荷重が大きいと考えられる風車を選定するので、風車の選定が誤っていた場合には、その結果、求められる運転制限の条件設定値も適切な値に設定されない可能性があり、信頼性が低いという欠点がある。
【００２７】
〔方法３〕
様々な風況条件と推定荷重とを関連付けたデータベースを予め用意しておき、各風車の荷重計算を行う際には、上記風況解析で得られた風況結果に最も近い推定荷重を上記データベースの中から取得することで、各風車の推定荷重を得ることとしてもよい。
【００２８】
具体的には、図４と図５に示すような、空気密度、ウィンドシアのべき指数、吹上角および乱れ度等の風況解析で求められる各種パラメータの組み合わせからなる荷重を風速ごとに予め用意しておき、各風車における所定期間の推定荷重は、各風車の風況解析結果に最も近い風況パラメータの風速毎の推定荷重を該データベースから取得し、風速風向出現頻度分布と組み合わせて積分することにより算出する。なお、風況解析結果に一致するケースが存在しない場合には、風況解析結果に近いケースの推定荷重を内挿することで求めることとしてもよい。
なお、上記べき指数は、ウィンドシアを規定する以下の（１）式で用いられる数値であり、以下の（１）式におけるαに相当する。
Ｖ（ｚ）＝Ｖｏ（Ｚ／Ｚｏ）α      （１）
【００２９】
各風車について推定荷重が求められると、続いて、運転制限の条件設定値が風車毎に決定される（図３のステップＳＡ３）。
運転制限の初期条件設定値は、以下の処理を風車毎に行うことで、決定される。
まず、ステップＳＡ２において算出された全風向の推定荷重の合計と風車の耐用年数から決定される既定の基準荷重とを比較し、全風向の推定荷重の合計が該基準荷重を超えていなければ運転制限は不要であると判断する。
【００３０】
一方、全風向の推定荷重の合計が該基準荷重を超えていた場合には、以下の手順により運転制限の初期条件設定値が決定される。
まず、風車干渉によって荷重が大きい風向、すなわち、運転制限が必要な風向を算定するために、上記ステップＳＡ２において算出された風向毎の推定荷重と予め設定されている風向毎の基準荷重（例えば、設計風況による荷重の１．２倍等）とを比較し、推定荷重が基準荷重を超えている風向を運転制限の対象風向として特定する。
【００３１】
次に、運転制限の対象風向において、運転制限の条件設定値としてカットアウト風速よりも低い風速、例えば、通常のカットアウト風速が２５ｍ／ｓであった場合には、２４ｍ／ｓ等を仮設定し、それ以上の風速出現頻度をゼロとして、上記ステップＳＡ２において実施した推定荷重の算出処理を再度実施し、全風向における推定荷重の合計を再度算出する。
続いて、再度算出した全風向における推定荷重の合計と上記基準荷重とを比較し、全風向における推定荷重の合計が基準荷重を超えていた場合には、運転制限の条件設定値を更に小さい値に仮設定し、推定荷重の再計算を行う。そして、全風向における推定荷重の合計が基準荷重以下となるまで、上記運転制限の条件設定値の調整を繰り返し行い、全風向における推定荷重の合計が基準荷重以下となったときの仮設定値を運転制限の初期条件設定値として決定する。
【００３２】
このようにして、各風車について運転制限の初期条件設定値が決定されると、これらの値が各風車と対応付けられて設定値記憶部２３に格納される。なお、このとき、上述したように、運転制限の初期条件設定値は、風車毎、かつ、風向毎に設定され、設定値記憶部２３に格納される。
上記処理までが事前検討で行われる処理内容である。
【００３３】
〔風車の運転制限調整〕
続いて、風車が実際に建設された後は、本実施形態に係る風車の運転制限調整装置によって、上記運転制限の初期条件設定値が評価され、定期的に、この条件設定値が最適な値になるように見直されることとなる。
【００３４】
以下、風車の運転制限調整方法について説明する。なお、以下の処理は、風車毎にそれぞれ行われる。
まず、図２に示した疲労等価荷重算出部２１は、風車に取り付けられている荷重センサから出力される荷重データを取得する（図３のステップＳＡ４）。この荷重データは、例えば、風車のブレードや主軸、主軸受、ナセルおよびタワー等に設置された歪ゲージによって取得された時系列データを所定の較正式を用いて荷重データに変換したものである。風車の運転制限調整装置１０が風車の外部に設けられていた場合には、該荷重データは無線または有線の通信媒体を介して本装置に送られてくるものとする。
【００３５】
疲労等価荷重算出部２１は、ステップＳＡ４において取得された荷重データ（時系列データ）を用いて、所定の期間（以下「分割期間」という）、例えば、１０分間毎に、かつ、所定の風向毎に、疲労等価荷重を算出する（図３のステップＳＡ５）。疲労等価荷重は、例えば、レインフローカウント法を採用した以下の（２）式を用いて算出される。
【００３６】
【数２】
2010048239.tif 000004
【００３７】
ここで、leq,kは該所定期間における処理開始からk番目の前記分割時間における疲労等価荷重、liは疲労荷重スペクトルのi番目の荷重レンジ、niは疲労荷重スペクトルのi番目の荷重レンジのサイクル数、neqは、前記第２の期間における基準等価サイクル数、mは関連材料のS-N曲線の勾配である。
【００３８】
続いて、上記疲労等価荷重を所定の期間（例えば、１年間）にわたって積分することで、年間における疲労等価荷重を風向毎に算出する。年間の疲労等価荷重は以下の（３）式を用いて算出される。
【００３９】
【数３】
2010048239.tif 000005
【００４０】
上記（３）式において、Leqは年間の疲労等価荷重、Ｎeqは、年間の基準等価サイクル数である。
【００４１】
次に、風向毎に算出された年間の疲労等価荷重を風車の耐用年数で決定される基準荷重とそれぞれ比較し、これらの比較結果から事前検討において決定された初期条件設定値が適切であるか否かを判断し、適切でなかった場合には、該初期条件設定値を更新する（図３のステップＳＡ６，ＳＡ７）。
【００４２】
具体的には、風向毎に算出された疲労等価荷重と事前検討において用いられた風向毎の基準荷重との差分を算出する。
続いて、上記差分が予め設定されている所定の閾値（たとえば、疲労等価荷重のプラスマイナス１０％）を超えているか否かを判断し、超えていた場合には、その差分に応じて当該風向における初期条件設定値を変更する。
【００４３】
例えば、疲労等価荷重が基準荷重よりも大きい風向に対しては、初期条件設定値を予め設定されている所定量小さく修正し、疲労等価荷重が基準荷重よりも小さい風向に対しては、初期条件設定値を予め設定されている所定量大きく修正する。所定量は、任意に設定できる。一例としては、１ｍ／ｓが挙げられる。
【００４４】
図６に、各風向における疲労等価荷重と基準荷重とを比較した一例を示す。図６において、実線は疲労等価荷重、破線は基準荷重を示している。図６においては、北を０°、東を９０°、南を１８０°、西を２７０°と設定し、全１６方位についてそれぞれ疲労等価荷重、基準荷重を算出する場合を想定している。東及び西において、疲労等価荷重と基準荷重との差分が閾値を超えている。
具体的には、東では疲労等価荷重が基準荷重を上回っているので、初期条件設定値を所定量小さく修正し、西では疲労等価荷重が基準荷重を下回っているので、初期条件設定値を所定量大きく修正する。
【００４５】
このようにして、各風車において初期条件設定値の見直しが行われると、設定値記憶部１３に格納されている運転制限の初期条件設定値は修正後の値に更新され、以降の運転制限には、変更後の条件設定値に基づいて行われることとなる。
そして、上述したステップＳＡ４からＳＡ６の処理を所定期間毎（例えば、１年毎）に繰り返し行うことにより、実際の風車の疲労損傷度に応じて運転制限の条件設定値を適切な値に更新することが可能となる。
【００４６】
以上、説明してきたように、本実施形態に係る風車の運転制限調整装置及び方法並びにプログラムによれば、風車を建設した後において、風車の荷重データを定期的に取得し、この荷重データに基づいて風車の疲労等価荷重を求め、この疲労等価荷重を用いて現在採用されている運転制限の条件設定値の見直しを行うので、風車の運用に応じて運転制限の条件設定値を適切な値に変更することが可能となる。これにより、運転制限によって風車が停止される頻度を低減させることができ、耐久年数を確保しながら、運転制限による発電量の低下を抑制することが可能となる。
【００４７】
なお、本実施形態では、事前検討において決定された運転制限の初期条件設定値を設定記憶部１３に格納することとし、事前検討自体は風車の運転制限調整装置において行われないものとしたが、例えば、事前検討における一連の処理についても風車の運転制限調整装置によって行われることとしてもよい。
【図面の簡単な説明】
【００４８】
【図１】本発明の一実施形態に係る風車の運転制限調整装置の概略構成を示した図である。
【図２】本発明の一実施形態に係る風車の運転制限調整装置の機能を展開して示した機能ブロック図である。
【図３】本発明の一実施形態に係る風車の運転制限調整装置によって実現される運転制限調整方法について説明するためのフローチャートである。
【図４】事前検討における推定荷重の算出方法について説明するための図である。
【図５】事前検討における推定荷重の算出方法について説明するための図である。
【図６】疲労等価荷重と基準荷重とを比較して示した分布図である。
【符号の説明】
【００４９】
１  ＣＰＵ
２  主記憶装置
３  補助記憶装置
４  入力装置
５  出力装置
６  通信装置
１１  疲労等価荷重算出部
１２  設定値更新部
１３  設定値記憶部
【図１】
2010048239.tif 000007
【図２】
2010048239.tif 000008
【図３】
2010048239.tif 000009
【図４】
2010048239.tif 000010
【図５】
2010048239.tif 000011
【図６】
2010048239.tif 000012