JP7282580B2 - 地熱発電プラントの坑井制御システム及びその坑井制御方法並びに坑井制御プログラム、地熱発電プラント - Google Patents

Description

本開示は、地熱発電プラントの坑井制御システム及びその坑井制御方法並びに坑井制御プログラム、地熱発電プラントに関するものである。

地熱発電プラントでは、坑井から噴気する地熱流体（主として蒸気及び熱水の混合流体）から分離した蒸気を蒸気タービンに導入して発電を行っている。このため、複数の坑井から噴出した地熱流体を合流して利用するにあたり、各坑井における坑井特性（噴気流体の圧力と流量の特性）に基づいて、噴気される地熱流体の流量等を制御し、発電電力の最大化等を図る坑井管理が行われている（例えば、特許文献１）。

特開平２－２３２４９５号公報

通常では、坑井の運用にあたり、例えば竣工時の噴気試験（初期に行われる試験等）により、図１４のような坑井特性（噴気圧力（坑口圧力）に対して蒸気流量と熱水流量とを関係付けた特性）を初期の坑井特性として取得する。初期の坑井特性を取得した後は、継続運用するために途中で坑井特性の再計測を行わないことが多い。坑井特性を取得するためには、一旦発電設備系統から解列等しなければならず、発電電力の低下等を避けるために、運用に伴う坑井特性の変化状況を把握することはなかった。

しかしながら、坑井は、地熱蒸気の噴出・抽出流量を増加しすぎると坑井の減衰が大きくなる場合があり、地熱蒸気の噴出・抽出流量を絞って圧力を上げすぎると流量が急に低下する場合がある。このため、坑井特性が変化しているにも関わらず、過去に取得した坑井特性に基づいて坑井管理を行った場合には、坑井の制御が不適切となる場合があり、坑井の減衰（短寿命化）を促進させる可能性があった。特定の坑井が減衰してしまうと、発電電力の低下等を避けるために減衰した坑井に代わる新たな坑井の掘削が必要となり、コストの増加を招く可能性もある。このため、坑井管理は、熟練者の豊富な経験のもとで運用されることが多く、坑井の噴気圧力と蒸気流量の適正な調整と、坑井の長寿命化が課題となっている。

また、複数の坑井から噴出された地熱流体の熱水と分離された後の蒸気総流量を用いて発電を実施している運用中において、地熱流体から分離され蒸気タービンにて発電に用いられる蒸気流量（各坑井から噴出された地熱流体から分離した蒸気の総流量）については計測することができる。しかしながら、蒸気タービンで発電に用いられた蒸気流量の総流量から、実際に各坑井から噴出された地熱流体に含まれる各坑井の蒸気流量を精度よく推定することは容易ではない。これは、各坑井から噴出された地熱流体から蒸気を分離する過程等において生ずる圧力変化などに起因して蒸気流量と熱水流量の混合状態が変動するためである。すなわち、各坑井から噴出された地熱流体に含まれる蒸気の総流量と、各坑井から噴出された地熱流体から熱水との混合状態から分離された蒸気の総流量とは必ずしも一致しない。このため、蒸気流量と熱水流量の混合状態の地熱流体から分離した蒸気の総流量に基づいて、各坑井の坑井特性（坑井の坑口圧力に対する熱水流量と蒸気流量の関係）を精度よく推定することは容易ではなく、例えば、各坑井から噴出された地熱流体の蒸気流量と、熱水と分離された後の蒸気流量とが一致するなど仮定条件のもとでの推定を要する。各坑井の運用中には、運用に伴って変化した図１４のような各坑井の坑井特性を正確に把握されていない。このため、正確に把握されていない坑井特性に基づいて坑井管理（坑井から噴気する地熱流体の流量制御）を行っても、精度のよい管理を行えない可能性があった。

また、複数の坑井から噴出された地熱流体を用いて発電を実施している運用中においては、気水分離器で地熱流体から蒸気を生成（分離）する前、すなわち地熱流体の状態では、熱水と蒸気が混合された状態であるため、各坑井から噴出された地熱流体に含まれる蒸気及び熱水の流量を正確に計測することは容易ではない。

本開示は、このような事情に鑑みてなされたものであって、指定期間における坑井特性を反映させ精度よく各坑井を制御することのできる地熱発電プラントの坑井制御システム及びその坑井制御方法並びに坑井制御プログラム、地熱発電プラントを提供することを目的とする。

本開示の第１態様は、予め設定される過去所定期間である指定期間において、複数の坑井から噴出した地熱流体の総流量及び総エンタルピーと、各前記坑井の坑口圧力とに基づいた運転データを蓄積し、蓄積した前記運転データに基づいて前記総流量の変動量及び前記総エンタルピーの変動量を算出し、各前記坑口圧力の基準圧力に対する偏差と前記総流量の前記変動量との関係性及び前記偏差と前記総エンタルピーの前記変動量との関係性に基づいて、各前記坑井における前記坑口圧力の組合せに対応した運定点を複数設定する運定点設定部を備える地熱発電プラントの坑井制御システムである。

複数の坑井から噴出された地熱流体を用いて発電を行う場合、地熱流体から分離され発電に用いられる蒸気流量から、各坑井から噴出された地熱流体に含まれる蒸気流量を精度よく推定することは容易ではない。これは、各坑井から噴出された地熱流体から蒸気を分離する過程等において生ずる圧力変化などに起因して、蒸気流量と熱水流量の混合状態が変動するためである。そこで、地熱流体の総流量及び総エンタルピーと、各坑井の坑口圧力とに基づいて、各坑口圧力の基準圧力に対する偏差と総流量の変動量との関係性及び各坑口圧力の基準圧力に対する偏差と総エンタルピーの変動量との関係性を設定することとした。すなわち、地熱流体の総流量及び総エンタルピーに基づくことで、各坑井の運用に伴う坑井特性の変化を反映することができる。このため、各坑口圧力の基準圧力に対する偏差から、地熱流体の総流量及び総エンタルピーの変動量の関係性を用いて、各坑井の坑口圧力の組合せに対応した複数の運定点を精度よく設定することが可能となる。

上記坑井制御システムにおいて、前記地熱流体の前記総流量と前記総エンタルピーとの関係に対して、前記地熱流体を用いた発電出力に係る情報を対応づけた特性に基づいて、所望の前記地熱流体の総流量状態及び発電出力状態となる前記運定点を所望運定点として特定する所望運定点特定部と、前記所望運定点に対応する各前記坑口圧力を設定する坑口圧力設定部と、を備えることとしてもよい。

上記のような構成によれば、地熱流体の総流量と総エンタルピーの関係に対して発電出力情報に係る情報を対応づけた特性を用いることで、所望の地熱流体の総流量状態及び発電出力が与えられれば、対応する所望運定点（地熱流体の総流量と総エンタルピーの組合せ）を特定することが可能となる。各坑井の坑口圧力を所望運定点に対応した坑口圧力に制御することで、実際の地熱発電プラントの運転状態をより最適化することが可能となる。

上記坑井制御システムにおいて、前記特性は、前記地熱流体の前記総流量と前記総エンタルピーとの関係が前記地熱流体より生成された蒸気流量に基づいて表されており、前記発電出力が前記蒸気流量に対応して表されていることとしてもよい。

上記のような構成によれば、地熱流体より生成された圧力条件と蒸気流量の計測値により、地熱流体の総流量と、総エンタルピーと、発電出力とを関係づけることが可能となる。このため、地熱流体の総流量及び発電出力の所望状態に対応した所望運定点（地熱流体の総流量と総エンタルピーの組合せ）を特定することが可能となる。

上記坑井制御システムにおいて、前記所望運定点特定部は、前記総流量が最小値に近づき、かつ、前記発電出力が最大値に近づくように前記所望運定点を特定することとしてもよい。

上記のような構成によれば、総流量が最小値に近づき、かつ、発電出力が最大値に近づくような所望運定点を特定することができる場合があるため、各坑井から噴出・抽出する地熱流体の総流量が増加して短寿命化することを抑制し、効率的に発電の高出力化を図ることができる。

上記坑井制御システムにおいて、前記所望運定点特定部は、前記総流量が最小値に近づき、かつ、前記発電出力が所定値に近づくように前記所望運定点を特定することとしてもよい。

上記のような構成によれば、総流量が最小値に近づき、かつ、発電出力が所定値に近づくような所望運定点を特定することができるため、各坑井の短寿命化を抑制し、必要とされる発電出力を所定値に維持することで効率的な運用ができる。なお、所定値とは、例えば要求発電出力である。

上記坑井制御システムにおいて、前記所望運定点特定部は、前記総流量が最小値に近づき、かつ、前記発電出力が所定値以上の範囲において前記所定値に近づくように前記所望運定点を特定することとしてもよい。

上記のような構成によれば、総流量が最小値に近づき、かつ、発電出力が所定値以上の範囲において該所定値に近づくような所望運定点を特定することができるため、各坑井の短寿命化を抑制し、必要とされる所定値以上の余剰な発電を抑えた効率的な運転を行うことが可能となる。なお、所定値とは、例えば、地熱発電プラントのブラックアウトを防止するために最低限必要な電力値や、最低限必要な電力値に所定の余裕度（マージン）を加算した値である。

上記坑井制御システムにおいて、各前記坑井における前記地熱流体の流量及び前記坑井の坑口圧力の少なくともいずれか一方の設定範囲を限定するように設定した坑井評価関数を用いて前記所望運定点を評価し、前記所望運定点を特定することとしてもよい。

上記のような構成によれば、各坑井における地熱流体の流量及び坑井の坑口圧力の少なくともいずれか一方の設定範囲が限定されるため、坑井の特性変化や噴気停止のおそれのある運定点が選定されることを防止でき、各坑井の短寿命化をより一層に抑制することができる。

本開示の第２態様は、坑井から噴出された地熱流体を蒸気と熱水に分離する分離器と、前記分離器により分離された蒸気を用いて発電を行う蒸気タービンと、上記の地熱発電プラントの坑井制御システムと、を備えた地熱発電プラントである。

本開示の第３態様は、予め設定される過去所定期間である指定期間において、複数の坑井から噴出した地熱流体の総流量及び総エンタルピーと、各前記坑井の坑口圧力とに基づいた運転データを蓄積し、蓄積した前記運転データに基づいて前記総流量の変動量及び前記総エンタルピーの変動量を算出し、各前記坑口圧力の基準圧力に対する偏差と前記総流量の前記変動量との関係性及び前記偏差と前記総エンタルピーの前記変動量との関係性に基づいて、各前記坑井における前記坑口圧力の組合せに対応した運定点を複数設定する運定点設定工程を有する地熱発電プラントの坑井制御方法である。

本開示の第４態様は、予め設定される過去所定期間である指定期間において、複数の坑井から噴出した地熱流体の総流量及び総エンタルピーと、各前記坑井の坑口圧力とに基づいた運転データを蓄積し、蓄積した前記運転データに基づいて前記総流量の変動量及び前記総エンタルピーの変動量を算出し、各前記坑口圧力の基準圧力に対する偏差と前記総流量の前記変動量との関係性及び前記偏差と前記総エンタルピーの前記変動量との関係性に基づいて、各前記坑井における前記坑口圧力の組合せに対応した運定点を複数設定する運定点設定処理をコンピュータに実行させるための地熱発電プラントの坑井制御プログラムである。

本開示によれば、指定期間における坑井特性を反映させ精度よく各坑井を制御することができるという効果を奏する。

本開示の第１実施形態に係る坑井制御システムを備えた地熱発電プラントの概略構成を示す図である。 本開示の第１実施形態に係る坑井制御システムを備えた地熱発電プラントが備える機能を示した機能ブロック図である。 本開示の第１実施形態に係る地熱流体の相対坑井特性において総流量の変動と総エンタルピーの変動を示す図である。 本開示の第１実施形態に係る複数の運定点をプロットした図である。 本開示の第１実施形態に係る地熱流体の総流量と総エンタルピーとの関係を示す図である。 本開示の第１実施形態に係る地熱流体の総流量と総エンタルピーとの関係を総蒸気流量に対して示した図である。 本開示の第１実施形態に係る地熱流体の総流量と総エンタルピーとの関係を発電出力に対して示した図である。 本開示の第１実施形態に係る運定点に対して得点を付与した例を示す図である。 本開示の第１実施形態に係る坑井制御システムによる各坑口圧力の設定処理のフローチャートを示した図である。 本開示の第２実施形態に係る運定点に対して得点を付与した例を示す図である。 本開示の第３実施形態に係る運定点に対して得点を付与した例を示す図である。 本開示の第４実施形態に係る坑井評価関数の例を示す図である。 本開示の第４実施形態に係る運定点に対して得点を付与した例を示す図である。 坑井の坑井特性の例を示す図である。

〔第１実施形態〕
以下に、本開示に係る地熱発電プラントの坑井制御システム及びその坑井制御方法並びに坑井制御プログラム、地熱発電プラントの第１実施形態について、図面を参照して説明する。なお、坑井制御システム２０は、坑井２から噴気される地熱流体によって発電を行う地熱発電プラントであれば幅広く適用できるものであって、以下に説明する構成の地熱発電プラント１のみに適用を限定されるものではない。

図１は、本開示の第１実施形態に係る坑井制御システムを備えた地熱発電プラント１の概略構成を示す図である。本実施形態では、坑井２を複数設けるものとして、例えば坑井２を２つ設ける場合について説明する。なお、図１に示す地熱発電プラント１は、例えばフラッシュサイクル型の地熱発電プラントであるが、坑井制御システム２０は、バイナリサイクル型など他の構成の地熱発電プラントであっても同様に適用することが可能である。また、図１では、例えば坑井２を２つ設け、気水分離器（分離器）５を２つ設け、還元井１９を１つ設けることとしているが、上記の構成と構成数量に限らず本開示を適宜適用することが可能である。

本実施形態に係る地熱発電プラント１は、図１に示すように、地熱流体輸送管（二相流体輸送管）３と、噴気流量調整弁（以下、単に「流量調整弁４」という）と、気水分離器（分離器）５と、熱水管６と、蒸気管７と、蒸気タービン８と、復水器１０と、冷却塔１１とを主な構成として備えている。また、地熱発電プラント１には、図２に示すような坑井制御システム（坑井制御装置）２０が適用される。

地熱流体輸送管３は、坑井（地熱坑井）２から噴出された地熱流体を気水分離器５へ導く管である。地下にはマグマ溜りが形成されており、この熱によって地下に浸透した雨水や流入した地下水等が加熱され地熱貯留層が形成される。地熱貯留層には、地熱流体が溜まっており、坑井２によって地上へ取り出される。地熱流体輸送管３は、坑井２を介して地熱貯留層から気水分離器５へ地熱流体を輸送している。なお、地熱流体とは、主として蒸気と熱水からなる二相混合流体である。

流量調整弁（噴気流量調整弁）４は、地熱流体輸送管３上に設けられており、坑井２から気水分離器５へ流入する地熱流体の総流量（噴気量、蒸気と熱水の合計流量）を調整している。なお、流量調整弁４の調整により坑井２の坑口圧力を調整することもできる。

また、地熱流体輸送管３上には、流量調整弁４の地熱流体流れの上流側に、坑口圧力（噴気圧力）を計測する圧力計１６が設けられている（図１における圧力計Ｐ１及びＰ２）。圧力計１６では、坑井２から噴気される地熱流体の圧力を計測している。なお、地熱流体輸送管３上における地熱流体流れの上流側（坑井２出口付近。例えば、圧力計１６と流量調整弁４との間）に、開閉弁１３を設け、地熱流体の導通状態（導通状態または非導通状態）を制御することとしてもよい。

気水分離器５は、地熱流体輸送管３により供給された二相混合流体である地熱流体を、蒸気と熱水に分離する装置（セパレータ）である。気水分離器５によって分離された熱水は熱水管６に導かれ、気水分離器５によって分離された蒸気は蒸気管７へ導かれる。

熱水管６は、気水分離器５によって分離された熱水を還元井１９へ導く管である。還元井１９を介して地下の地熱貯留層に熱水を還すことで、地下の地熱貯留層における地熱流体の枯渇を抑制する。各気水分離器５によって分離された熱水は熱水管６で合流し、ポンプ１４を介して還元井１９へ圧送される。なお、自圧によって還元井１９へ圧送可能な場合には、ポンプ１４を省略したり小容量化してもよい。また、還元井１９に送られる熱水の流量は、熱水管６上に設けた流量計１７によって計測される。なお、本実施形態では、各気水分離器５によって分離された熱水が合流した後の熱水の流量（合計流量）を計測する場合について説明するが、各気水分離器５によって分離された熱水の流量をそれぞれ計測することとしてもよい。また、還元井１９が複数設けられることとしてもよい。

蒸気管７は、気水分離器５によって分離された蒸気を蒸気タービン８へ導く管である。各気水分離器５によって分離された蒸気は、蒸気管７で合流して蒸気タービン８へ供給される。また、蒸気タービン８に送られる蒸気の流量は、蒸気管７上に設けた流量計１８によって計測される。なお、本実施形態では、各気水分離器５によって分離された蒸気が合流した後の蒸気の流量（合計流量）を計測する場合について説明するが、各気水分離器５によって分離された蒸気の流量をそれぞれ計測することとしてもよい。また、蒸気管７には、各気水分離器５によって分離された蒸気が合流した後の蒸気に対して流量調整弁１５を設け、蒸気タービン８への蒸気の供給流量を制御することとしてもよい。

蒸気タービン８は、蒸気管７により供給された蒸気のエネルギーによってタービン翼（不図示）を回転駆動させ、タービン翼の回転軸に接続された発電機９を回転駆動して発電を行う。

復水器１０は、蒸気タービン８においてタービン翼（不図示）を回転駆動させる仕事をし終えた蒸気を冷却水で冷却することで凝縮して復水する装置である。復水された温水は、冷却塔１１へ供給される。

冷却塔１１は、復水器１０において復水された温水を蒸発冷却する装置である。具体的には、冷却塔１１に供給された温水は冷却塔１１上部の散布部から散布される。散布された温水は、冷却塔１１の送風機によって流通する空気と接触することで一部が蒸発し、この蒸発に伴う潜熱によって他の温水が冷やされて水になる。冷やされた水は、冷却塔水として冷却塔１１の水槽に貯水される。水槽に貯水されている冷却塔水は、ポンプ１２を介して冷却水として復水器１０へ供給される。

坑井制御システム２０は、所望の運定点（地熱流体の総流量と総エンタルピーの組合せ）に対して各坑井２を制御する。複数の各坑井２から噴出・抽出された地熱流体は、流量調整弁４及び気水分離器５を介して熱水と蒸気とに分離されるが、この分離過程で発生する圧力変化（例えば、流量調整弁４での減圧フラッシュ）によって、分離の前後での蒸気流量の比率は同じにならない場合がある。具体的には、気水分離器５での分離過程で発生する減圧作用により、熱水流量が減少するとともに蒸気流量は増加する。すなわち、複数の各坑井２から噴出・抽出された地熱流体の蒸気流量と、気水分離器５で分離された後の蒸気流量とは一致しない場合がある。熱水も同様に気水分離器５で分離される前後で流量は一致しない場合がある。このため、気水分離器５における分離後の蒸気流量及び熱水流量に基づいても、各坑井２から噴出された蒸気流量及び熱水流量（地熱流体に含まれる蒸気流量及び熱水流量）の絶対値を推定することは容易ではない。すなわち、各坑井２の坑井特性（坑口圧力に対する蒸気流量と熱水流量との関係）を実際の状況（現時点における状態）に対応して把握することは容易ではない。すなわち、例えば竣工時の噴気試験で取得した初期の坑井特性を用いて、実際の状況を把握されていない（運用に伴う変化が正確に反映されていない）各坑井特性に基づいて各坑井２を制御したとしても、適切に制御できない場合がある。そこで、本実施形態に係る坑井制御システム２０では、各坑井２の坑井特性を図１４に示したような噴気圧力に対する蒸気流量と熱水流量の関係を坑井特性で表したものとして把握するのではなく、地熱流体の総流量及び総エンタルピーと各坑井２の坑口圧力の状態とに基づいて推定した相対的な坑井２の特性に基づいて各坑井２を精度よく制御する。

坑井制御システム２０は、例えば、図示しないＣＰＵ（中央演算装置）、ＲＡＭ（Ｒａｎｄｏｍ Ａｃｃｅｓｓ Ｍｅｍｏｒｙ）等のメモリ、及びコンピュータ読み取り可能な記録媒体等から構成されている。後述の各種機能を実現するための一連の処理の過程は、プログラムの形式で記録媒体等に記録されており、このプログラムをＣＰＵがＲＡＭ等に読み出して、情報の加工・演算処理を実行することにより、後述の各種機能が実現される。なお、プログラムは、ＲＯＭやその他の記憶媒体に予めインストールしておく形態や、コンピュータ読み取り可能な記憶媒体に記憶された状態で提供される形態、有線または無線による通信手段を介して配信される形態等が適用されてもよい。コンピュータ読み取り可能な記憶媒体とは、磁気ディスク、光磁気ディスク、ＣＤ－ＲＯＭ、ＤＶＤ－ＲＯＭ、半導体メモリ等である。

図２は、坑井制御システム２０が備える機能を示した機能ブロック図である。図２に示されるように、坑井制御システム２０は、運定点設定部２１と、所望運定点特定部２２と、坑口圧力設定部２３とを備えている。

坑井特性によれば、図１４に示すように噴気圧力（坑口圧力）が決定されると、蒸気流量と熱水流量が一義的に特定される。すなわち、各坑井２では、坑口圧力に対して、噴出された地熱流体に含まれる蒸気流量と熱水流量が対応している。ここで、本実施形態では、ｉ番目の坑井２の蒸気流量をＱｓｉとし、熱水流量をＱｂｉとすると、ｉ番目の坑井２から噴出された地熱流体の流量Ｑｉは、蒸気流量と熱水流量の和で表される（Ｑｉ＝Ｑｓｉ＋Ｑｂｉ）。また、ｉ番目の坑井２から噴出された地熱流体のエンタルピーＨｉは、蒸気の飽和比エンタルピーをｈｓｉとし、熱水の飽和比エンタルピーをｈｂｉとすると、蒸気流量と蒸気の飽和比エンタルピーとの積である蒸気エンタルピーと、熱水流量と熱水の飽和比エンタルピーとの積である熱水エンタルピーとの和で表される（Ｈｉ＝ｈｓｉ・Ｑｓｉ＋ｈｂｉ・Ｑｂｉ）。このため、各坑井２から噴出された地熱流体の総流量Ｑは、各坑井２から噴出された地熱流体の流量Ｑｉの総和として定義され、各坑井２から噴出された地熱流体の総エンタルピーＨは、各坑井２から噴出された地熱流体のエンタルピーＨｉの総和として定義される。すなわち、各坑井２の坑井特性と、地熱流体の総エンタルピーＨと地熱流体の総流量Ｑとは互いに関係している。このように、各坑井２の坑井特性を、地熱流体の総エンタルピーＨと地熱流体の総流量Ｑ（坑井２における熱水と蒸気の合計流量）へ変換することで、相対的な坑井特性を表すことができる。

実際の地熱発電プラントにおいて地熱流体の総流量Ｑ及び地熱流体の総エンタルピーＨを求める際は、後述の（６）式及び（７）式を用いて計測算定する。すなわち、総蒸気流量Ｑｓｓ及び総熱水流量Ｑｂｓの計測値を用いて、地熱流体の総流量Ｑ及び地熱流体の総エンタルピーＨを算出する。このため、各坑井２の個々の坑井特性を正確に把握していなくても、各坑井の個々の坑井特性と関係性を変換して保有する地熱流体の総流量Ｑ及び地熱流体の総エンタルピーＨを取得することができ、これらに基づいて相対的に坑井特性を反映して最適な運定点（各坑井２の坑口圧力）を特定することが可能となる。ここで気水分離器５により分離後の総蒸気流量Ｑｓｓ及び総熱水流量Ｑｂｓは流量計で計測される蒸気及び熱水の流量であり、ｈｓｓ、ｈｂｓは各流量計における図示しない圧力計、温度計による計測結果をもとに得られる物性値である蒸気、熱水の比エンタルピーである。なお実際のプラントにおいて流量計が複数ある場合は、それぞれの流量計におけるエンタルピーを個別に算出したうえで積算するものとする。

運定点設定部２１は、複数の坑井２から噴出した地熱流体の総流量及び総エンタルピーの計測算定値と各坑井２の坑口圧力の状態とに基づいて設定した相対坑井特性に対して、各坑口圧力の基準圧力に対する偏差と総流量の変動量との関係性及び偏差と総エンタルピーの変動量との関係性を用い、各坑井２の制御可能範囲内における各坑口圧力の組合せに対応した複数の運定点を設定する。

具体的には、運定点設定部２１には、各坑口圧力の基準圧力に対する偏差と総流量の変動量との関係性及び各坑口圧力の基準圧力に対する偏差と総エンタルピーの変動量との関係性が格納される。なお、該関係性については、運定点設定部２１において自動的に設定することとしてもよいし、地熱発電プラント１の運転員等によって設定されてもよいし、他の装置等における演算結果として取得されることとしてもよい。

まず、各坑井２の噴気圧力（坑口圧力）Ｐｉ（本実施形態の例ではｉ＝１、２）について、基準圧力をＰｉ０とし、基準圧力Ｐｉ０からの変動量をΔＰｉとすると、運定点での噴気圧力（坑口圧力）Ｐｉは、以下（１）式で表される。

（１）式において、基準圧力Ｐｉ０は、各坑井２に対してそれぞれ設定され、通常運転される圧力範囲内であれば任意に設定可能である。そして、各坑井２の坑口圧力が基準圧力Ｐｉ０である場合における地熱流体の総流量を基準総流量Ｑ０（例えば単位は、ｋｇ／ｓで示される）、総エンタルピーを基準総エンタルピーＨ０（例えば単位は、ｋＷで示される）とすると、各坑口圧力が任意の状態における地熱流体の総流量Ｑ及び総エンタルピーＨは、以下（２）式及び（３）式で表される。

ここで、（２）式におけるΔＱ（ΔＰｉ）は、総流量Ｑの基準総流量Ｑ０からの変動量であり、各坑井２の坑口圧力の変動量ΔＰｉの関数となる。また、（３）式におけるΔＨ（ΔＰｉ）は、総エンタルピーＨの基準総エンタルピーＨ０からの変動量であり、各坑井２の坑口圧力の変動量ΔＰｉの関数となる。すなわち、（２）式及び（３）式で示される関係を図示すると、図３のような相対坑井特性となる。すなわち、各坑井２の坑口圧力が基準圧力Ｐｉ０である場合の地熱流体の総流量（基準総流量Ｑ０）と総エンタルピー（基準総エンタルピーＨ０）を基準運定点（白抜き丸印）として、各坑口圧力が基準圧力Ｐｉ０から変動（ΔＰ１、ΔＰ２）すると、総流量Ｑ及び総エンタルピーＨがそれぞれΔＱ及びΔＨだけ変動する。すなわち、各坑口圧力が任意の状態での運定点（丸囲み×印）は、基準運定点から移動した位置に表示できる。

ここで、基準総流量Ｑ０からの変動量ΔＱ（ΔＰｉ）及び基準総エンタルピーＨ０からの変動量ΔＨ（ΔＰｉ）は、それぞれ各坑口圧力の変動量ΔＰｉの関数となっている。これらの関数は、坑井２の運転とともに（例えば経時的に）変化する可能性がある。このため、変動量ΔＱ及び変動量ΔＨを後述する（４）式及び（５）式のように係数を用いて表し、各係数を地熱発電プラント１の運転状態に合わせて更新することで、各坑井２の坑井特性の変化を反映させる。具体的には、指定期間（過去所定期間）において、地熱発電プラント１が運転中に特定される地熱流体の総流量Ｑ及び総エンタルピーＨと、各坑井２の坑口圧力とに基づいた運転データを蓄積する。そして、蓄積した運転データに基づいて（４）式及び（５）式の係数を都度適正化し、適正化した各係数を用いて変動量ΔＱ及び変動量ΔＨを算出する。すなわち、総流量の変動量ΔＱ及び総エンタルピーの変動量ΔＨの関数を、運用中の地熱発電プラント１から取得が可能で、複数に取得した情報（総流量の計測算定値、総エンタルピーの計測算定値、各坑口圧力の状態）に基づいて運定点を設定することで、指定期間における坑井２の相対坑井特性を把握することができる。そして、総流量の変動量ΔＱ及び総エンタルピーの変動量ΔＨを指定期間における坑井２の特性に適合するように算出することができる。すなわち、運用や時間経過に伴って各坑井の坑井特性が変化したとしても、各坑井の坑井特性自体を再現することは困難であっても、指定期間における相対坑井特性を反映して総流量の変動量ΔＱ及び総エンタルピーの変動量ΔＨを算出できる。

なお、総流量Ｑの計測算定値は、後述する（６）式を用いて、気水分離器５での分離後における総蒸気流量（Ｑｓｓ）の計測値と総熱水流量（Ｑｂｓ）の計測値を合計した値として算出される。また、総エンタルピーＨの計測算定値は、後述する（７）式を用いて、総蒸気流量（Ｑｓｓ）の計測値及び総熱水流量（Ｑｂｓ）の計測値と、気水分離器５では圧力により飽和温度が決定されるため、気水分離器５内の圧力により決まる物性値である蒸気の飽和比エンタルピー（ｈｓｓ）及び熱水の飽和比エンタルピー（ｈｂｓ）とを用いて算出される。

総流量の変動量ΔＱ及び総エンタルピーの変動量ΔＨについては、本実施形態では２つの坑井２で構成されることから、例えば、以下（４）式及び（５）式で表すことができる。

ここで、（４）式におけるａ１、ａ２は、それぞれΔＰ１及びΔＰ２の係数である。また、（５）式におけるｂ１、ｂ２は、それぞれΔＰ１及びΔＰ２の係数である。そして、総流量の変動量ΔＱ及び総エンタルピーの変動量ΔＨを（４）式及び（５）式で表した場合には、後述するように、計測した総蒸気流量及び総熱水流量に基づいて算出した地熱流体の総流量Ｑ及び総エンタルピーＨから、基準総流量及び基準総エンタルピーとの差が変動量ΔＱ及びΔＨとして算出される。また、総蒸気流量及び総熱水流量を計測した時点における各坑口圧力（Ｐ１、Ｐ２）を取得して、基準圧力との差をΔＰ１及びΔＰ２として算出する。すなわち、指定期間内において計測を行った時点に対して、地熱流体の総流量Ｑに係る運転データ情報と、総エンタルピーＨに係る運転データ情報と、各坑口圧力（Ｐ１、Ｐ２）に係る運転データ情報を蓄積することによって、指定期間において、ΔＱ、ΔＨ、ΔＰ１、ΔＰ２の組合せを複数パターン取得する。

そして、指定期間において取得した情報（ΔＱ、ΔＨ、ΔＰ１、ΔＰ２の組合せ）に対して例えば最小二乗法等によるフィッティング手法を適用することで、（４）式及び（５）式における各係数（ａ１、ａ２、ｂ１、ｂ２）を特定する。このようにして、実際の地熱発電プラント１の運転状態を示す総蒸気流量及び総熱水流量の計測値を、総流量の変動量ΔＱ及び総エンタルピーの変動量ΔＨを算出する式（（４）式及び（５）式における各係数）に反映させる。指定期間とは、例えば最近の１ヶ月から６ヶ月や最近の１年間などであり、予め設定される期間（過去所定期間）である。すなわち、地熱流体の総流量Ｑの計測算定値及び総エンタルピーＨの計測算定値と、計測した時点における各坑口圧力（Ｐ１、Ｐ２）とに基づいて、各係数が指定期間における坑井２の特性に適合したものとすることができるため、相対坑井特性を反映して、総流量の変動量ΔＱ及び総エンタルピーの変動量ΔＨを算出することが可能となる。なお、坑井２の数量が２より多い場合（例えば、坑井２の数量＝ｎ）は、ΔＱの係数をａｎ・ΔＰｎまで増加して加算し、ΔＨの係数をｂｎ・ΔＰｎまで増加して加算する。このように、坑井２の数に限定されず、運用に伴って変化した各坑井２の坑井特性を相対的に加味して、各坑井２の圧力変化に対する総流量の変動量ΔＱ及び総エンタルピーの変動量ΔＨをでき、運定点を精度よく推定することが可能となる。

なお、総流量の変動量ΔＱを表す関数及び総エンタルピーの変動量ΔＨを表す関数については、上記の（４）式及び（５）式に限らず適用することができる（例えば任意の次数の関数等）。また、総流量の変動量ΔＱを表す関数及び総エンタルピーの変動量ΔＨを任意の関数として、機械学習等によって求めることとしてもよい。すなわち、指定期間における坑井２の特性を反映して、各坑口圧力（Ｐ１、Ｐ２）から地熱流体の総流量の変動量ΔＱ及び総エンタルピーの変動量ΔＨを算出可能であれば、上記に限定されず適用することができる。

このように、各坑口圧力の変動量ΔＰｉに基づいて、指定期間の状態における坑井２の特性に適合するように総流量の変動量ΔＱ及び総エンタルピーの変動量ΔＨを算出することができるため、基準運定点（基準総流量Ｑ０及び基準総エンタルピーＨ０）に対して、算出した総流量の変動量ΔＱ及び総エンタルピーの変動量ΔＨを加算することで、任意の坑口圧力（Ｐ１、Ｐ２）における地熱流体の総流量Ｑ（運定点）及び総エンタルピーＨ（運定点）を算出することができる（図３）。

運定点設定部２１は、各坑井２の制御可能範囲内における各坑口圧力（Ｐ１、Ｐ２）の組合せに対応した複数の運定点を設定する。制御可能範囲とは、特性に応じて坑井２毎に設定されているものであり、流量調整弁４を用いて制御可能な坑口圧力の範囲である。なお、制御可能範囲は、実際に弁操作で実現可能な圧力範囲のうち、坑井として運用可能な範囲を指している。運定点設定部２１は、後述する所望運定点特定部２２における所望運定点の特定のために、各坑井２の制御可能範囲内において取り得る様々な運定点を設定する。すなわち、指定期間（過去所定期間）において、地熱発電プラント１が運転中に特定される地熱流体の総流量Ｑ及び総エンタルピーＨと、各坑井２の坑口圧力とに基づいた運転データを蓄積した中から、運定点設定部２１は、実際の運転状態として取り得る運定点のパターンを、係数が更新された（４）式及び（５）式を用いて設定する。

具体的には、運定点設定部２１は、各坑井２の制御可能範囲において、坑井２毎に任意の坑口圧力（Ｐ１、Ｐ２）を選定する。このようにして、運定点設定部２１は、各坑口圧力の組合せを複数組設定する。そして、設定した各坑口圧力の組合せでの総流量の変動量ΔＱ及び総エンタルピーの変動量ΔＨを（４）式及び（５）式を用いて算出する。そして、総流量の変動量ΔＱ及び総エンタルピーの変動量ΔＨと、基準運定点（基準総流量Ｑ０及び基準総エンタルピーＨ０）とに基づいて、各坑口圧力の組合せにおける運定点（総流量Ｑ及び総エンタルピーＨ）を算出する。運定点設定部２１は、運定点算出の処理を、設定した複数の各坑口圧力の組合せ毎に実行する。すなわち、運定点設定部２１では、各坑井の坑口圧力（Ｐ１、Ｐ２）と、該坑口圧力に対応した運定点（総流量Ｑ及び総エンタルピーＨ）との組み合わせが、複数パターン設定されることとなる。このように設定した複数の運定点を図示したものを図４に示す。図４における各×印が、設定された複数の運定点（運定点群）を示す。すなわち、図４に示す運定点群のそれぞれには、各坑井の坑口圧力の組合せが対応づけられている。後述する所望運定点特定部２２では、このように設定された複数の運定点のうち、所望の条件に適合した運定点を所望運定点として特定する。

所望運定点特定部２２は、複数の坑井２から噴出した地熱流体の総流量Ｑと総エンタルピーＨとの関係が地熱流体を用いた発電出力毎に表された特性に基づいて、所望の地熱流体の総流量状態及び発電出力状態となる運定点を所望運定点として特定する。すなわち、所望運定点特定部２２は、運定点設定部２１において設定された複数の運定点の中から、所望の運定点（所望の運転状態における総流量と総エンタルピーの組合せ）を特定する。

所望運定点特定部２２において用いる特性は、地熱流体の総流量Ｑと総エンタルピーＨとの関係が地熱流体より生成された蒸気流量に基づいて表されており、発電出力Ｌが蒸気流量に基づいて表されている関係性である。このような特性について説明する。なお、以降の説明では、本実施形態での一例として複数の各気水分離器（セパレータ）５での地熱流体を蒸気と熱水に分離するセパレータ圧力を等しいとする仮定条件のもとでの簡易な算出方法を示すものであり、例えば複数の各セパレータ圧力が異なる場合や、地熱エネルギーをバイナリサイクル方式で利用する場合などでは、別途各々に合せたモデルを設定する。地熱流体の総流量Ｑ及び総エンタルピーＨは、気水分離器５で分離後の総蒸気流量Ｑｓｓ及び総熱水流量Ｑｂｓを用いて以下（６）式及び（７）式で表される。

（７）式において、ｈｓｓは気水分離器５内の圧力Ｐｓにおける蒸気の飽和比エンタルピー（例えば単位は、ｋＪ／ｋｇで示される）であり、飽和温度は圧力Ｐｓで決定されるため、ｈｂｓは気水分離器５内の圧力Ｐｓにおける熱水の飽和比エンタルピー（例えば単位は、ｋＪ／ｋｇで示される）であり、圧力Ｐｓにより決まる物性値である。ここで、本実施形態では、気水分離器５は２つ設けることとしているが、各気水分離器５において出口から蒸気タービン８の入口までの圧力損失は比較的小さいため、圧力損失を省略（もしくは一定）とすると、各気水分離器５内の圧力は蒸気タービン８の入口の圧力と等しいと考えることができる。このため、蒸気タービン８の入口の圧力に基づいて特定される蒸気の飽和比エンタルピーｈｓｓ及び熱水の飽和比エンタルピーｈｂｓは各気水分離器５内の圧力Ｐｓで決定する値に等しいと考えることができ、地熱流体の総流量Ｑ及び総エンタルピーＨは、分離後の総蒸気流量Ｑｓｓ及び総熱水流量Ｑｂｓを用いて上記（６）式及び（７）式のように表すことができる。

（６）式及び（７）式を用いて、Ｑｂｓを消去すると、以下の関係式を得る。

（８）式は、地熱流体の総流量Ｑ及び総エンタルピーＨがわかれば、総蒸気流量Ｑｓｓを算出可能なことを示している。（８）式をさらに変形すると、以下の関係を得る。

ここで、（９）式において総蒸気流量Ｑｓｓを一定とすると、総エンタルピーＨは、傾きをｈｂｓとする総流量Ｑの一次関数で表される。すなわち、（９）式は、図５のような直線関係（線形）のグラフを描くことができる。ここで、図５では、総蒸気流量Ｑｓｓをパラメータとして、ある固定値（一定値）としている。このため、総蒸気流量Ｑｓｓをより大きな固定値とした場合には一次関数の切片（（ｈｓｓ－ｈｂｓ）・Ｑｓｓ）が大きくなり、図６に示すように、縦軸（総エンタルピー）が増加する方向へ一次関数が平行移動する。一方で、総蒸気流量Ｑｓｓをより小さな固定値とした場合には一次関数の切片（（ｈｓｓ－ｈｂｓ）・Ｑｓｓ）が小さくなり、図６に示すように、縦軸（総エンタルピー）が減少する方向へ一次関数が平行移動する。すなわち、（９）式の関係から、地熱流体の総流量Ｑと総エンタルピーＨとの関係が、パラメータである地熱流体より生成された総蒸気流量Ｑｓｓに基づいて表されたこととなる。

一方で、蒸気タービン８では、地熱流体より生成された蒸気（総蒸気流量Ｑｓｓ）によって発電を行っており、発電出力（発電量）Ｌ（例えば単位は、ｋＷで示される）は、以下（１０）式で表される。

（１０）式において、Δｈｔは蒸気タービン８の有効熱落差（例えば単位は、ｋＪ／ｋｇで示される）で、蒸気タービン８の出口蒸気の運用により決定される。すなわち、（１０）式では、発電出力Ｌが蒸気タービン８に供給された総蒸気流量Ｑｓｓによって示されている。すなわち、蒸気タービン８の有効熱落差Δｈｔが一定とみなせる条件範囲では、発電出力Ｌと総蒸気流量Ｑｓｓとは比例関係にあるため、図６に示すような総蒸気流量Ｑｓｓと一次関数の関係は、図７に示すように、発電出力Ｌと一次関数の関係と考えることができる。すなわち、総蒸気流量Ｑｓｓを増加させると、総エンタルピーＨが増加するとともに、発電出力Ｌも増加する。一方で、総蒸気流量Ｑｓｓを減少させると、総エンタルピーＨが減少するとともに、発電出力Ｌも減少する。

すなわち、（９）式及び（１０）式を整理すると、図７に示すような、地熱流体の総流量Ｑと総エンタルピーＨとの関係が地熱流体を用いた発電出力Ｌ毎に表された特性を得ることができる。すなわち、坑井２の短寿命化の要因の一つとなる地熱流体の総流量Ｑの状態（Ｑ＝Ｑｓｓ＋Ｑｂｓ）と、地熱発電プラント１における発電出力状態（Ｌ）において、所望の運定点を特定することが可能となる。

本実施形態では、所望運定点特定部２２は、各坑井２の制御可能範囲において、総流量Ｑが最小値に近づき、かつ、発電出力Ｌが最大値に近づくように所望運定点を特定する。地熱発電プラント１の各坑井２の運用として、地熱流体の総流量Ｑが多すぎると、地熱貯留層に貯留されている地下流体の水源が乏しい場合などでは、地熱流体の減少を促進し、坑井２の寿命を余分に短くしてしまう可能性がある。しかしながら、地熱発電プラント１としては、発電出力Ｌはできるだけ大きい方が好ましい。そこで、本実施形態では、総流量Ｑが最小値に近づき、かつ、発電出力Ｌが最大値に近づくような運定点を特定する。このため、所望運定点特定部２２は、図７に示すような地熱流体の総流量Ｑと総エンタルピーＨとの関係が地熱流体を用いた発電出力Ｌ毎に表された特性に対して、運定点設定部２１にて設定した複数の運定点をプロットし、所望の運定点を特定する（図８）。

所望運定点特定部２２は、所望の運定点を特定可能なように、得点付けによる評価を行う。本実施形態では、総流量Ｑが最小値に近づき、かつ、発電出力Ｌが最大値に近づくように運定点を特定する。このため、例えば本実施形態に係る所望運定点特定部２２は、例えば、図４で示した運定点群を図７に対してプロットすることで、図８を得る。図８に示すように、プロットした運定点群のうち総流量Ｑが最小値となる運定点を最も高得点となるように設定し、総流量Ｑが多くなるにつれて所定の間隔で得点を低くして設定する。図８の例では、プロットした運定点群のうち総流量が最小値となる運定点に最高得点として１０点を付け、総流量Ｑが多くなるにつれて所定の間隔で９点、８点等と点数が低下するように付ける。また、所望運定点特定部２２は、プロットした運定点群のうち発電出力Ｌが最大値となる運定点を最も高得点となるように設定し、発電出力Ｌが少なくなるにつれて所定の間隔で得点を低くして設定する。図８の例では、プロットした運定点群のうち発電出力Ｌが最大値となる運定点に最高得点として１０点を付け、発電出力Ｌが少なくなるにつれて所定の間隔で９点、８点等と点数が低下するように付ける。

図８のように得点を付与した場合、総流量の最小値へ近づく得点と発電出力Ｌが最大値へ近づく得点を合算して評価する。例えば丸囲み×印で示したＷ１が最も高得点となる。このため、所望運定点特定部２２は、総流量Ｑが最小値に近づき、かつ、発電出力Ｌが最大値に近づく最適な運定点としてＷ１を特定する。このようにして、所望運定点特定部２２は、最適となる所望の運定点を特定する。なお、本実施形態では、得点を付与することにより所望運定点を特定しているが、所望の運転状態（総流量が最小値に近づき、かつ、発電出力が最大値に近づく状態）の運定点を特定することができれば、得点付けの方法に限定されず適用することができる。

坑口圧力設定部２３は、所望運定点に対応する各坑口圧力（Ｐ１ｓｅｔ、Ｐ２ｓｅｔ）を設定する。運定点設定部２１は、各坑口圧力の組合せと運定点とを対応させている。このため、所望運定点特定部２２において所望運定点として特定された運定点Ｗ１には、すでに対応する各坑口圧力（Ｐ１、Ｐ２）が対応づけられている。このため、坑口圧力設定部２３は、所望運定点に対応づけられた各坑口圧力の組合せを、所望の各坑口圧力（Ｐ１ｓｅｔ、Ｐ２ｓｅｔ）として設定する。なお、坑口圧力設定部２３は、各坑口圧力を設定後、実際に流量調整弁４を制御して各坑井２の坑口圧力を設定した坑口圧力（Ｐ１ｓｅｔ、Ｐ２ｓｅｔ）となるように調整し、地熱発電プラント１の実際の運定点を所望の運定点へ近づけることとしてもよい。なお、坑口圧力設定部２３によって設定した各坑口圧力に基づいて、地熱発電プラント１の運転員等が流量調整弁４を制御することとしてもよい。

次に、上述の坑井制御システム２０による各坑口圧力の設定処理について図９を参照して説明する。図９に示すフローは、例えば、地熱発電プラント１が稼働中において、所定の制御周期で繰り返し実行される。なお、図９に示すフローは、坑井制御を実施する場合において、運転員等によって開始指令が入力された場合に実行することとしてもよい。

まず、指定期間（過去所定期間）において、地熱発電プラント１が運転中に特定される地熱流体の総流量Ｑ及び総エンタルピーＨと、各坑井２の坑口圧力とに基づいた運転データを蓄積した中から、各坑口圧力（Ｐ１、Ｐ２）の組合せに対応した運定点を複数設定する（Ｓ１０１）。例えば、図４や図８で各×印で示した運定点群を設定する。

次に、地熱流体の総流量Ｑと総エンタルピーＨとの関係が地熱流体を用いた発電出力Ｌ毎に表された特性に基づいて、設定した各運定点に得点を付与する（Ｓ１０２）。

次に、得点に基づいて運定点から最適となる所望運定点を特定する（Ｓ１０３）。例えば、図８で丸囲み×印で示したＷ１を特定する。

次に、特定した所望運定点Ｗ１に基づいて、対応する各坑井２の坑口圧力（Ｐ１ｓｅｔ、Ｐ２ｓｅｔ）を設定する（Ｓ１０４）。

なお、本実施形態では、発電出力Ｌと総流量Ｑの得点付けにおける最高得点を同じ値としているが、得点付けの方法は一例であり、上記方法に限定されない。総流量Ｑの減少への得点と、発電出力Ｌの増大への得点につき、地熱発電プラント１としてより重視したい得点へ重みを付与してもよい。すなわち、総流量Ｑの増加及び発電出力Ｌの減少の所定の間隔でさらに重みの勾配を付加した得点として、大きく低くなるように設定してもよい。総流量Ｑの増加及び発電出力Ｌの減少のいずれかに異なる重みを持たせて、最高得点を異なる値としてもよい。すなわち、発電出力Ｌを重視する場合には、発電出力に関する得点付けの最高得点を、総流量Ｑに関する得点付けの最高得点に対して高くすることとしてもよい。

以上説明したように、本実施形態に係る地熱発電プラントの坑井制御システム及びその坑井制御方法並びに坑井制御プログラム、地熱発電プラントによれば、指定期間における坑井特性を反映させ精度よく各坑井を制御することが可能となる。複数の坑井２から噴出・抽出された地熱流体を用いて発電を行う場合、気水分離器５などで地熱流体から分離され発電に用いられる蒸気流量により、実際に各坑井２から噴出・抽出された地熱流体に含まれる蒸気流量を精度よく推定することは容易ではない。これは、複数の各坑井２から噴出された地熱流体から蒸気を気水分離器５などで分離する過程等において生じる圧力変化などに起因して蒸気流量と熱水流量の混在状態が変動するためである。すなわち、精度よく複数の各坑井２の坑井特性を把握することは容易ではなかった。そこで、地熱流体の総流量Ｑの計測算定値及び総エンタルピーＨの計測算定値と、複数の各坑井２の噴気圧力（坑口圧力）Ｐｉの状態とに基づいて、各坑口圧力の基準圧力Ｐｉ０に対する変動量ΔＰｉと総流量の変動量ΔＱとの関係性及び各坑口圧力の基準圧力に対する変動量ΔＰｉと総エンタルピーの変動量ΔＨとの関係性を指定期間内における計測値に基づいて設定することとした。すなわち、地熱流体の総流量Ｑの計測算定値及び総エンタルピーＨの計測算定値に基づくことで、運用に伴い変化した坑井２の坑井特性を反映することができる。すなわち、各坑口圧力の基準圧力に対する偏差から、地熱流体の総流量の変動量ΔＱ及び総エンタルピーの変動量ΔＨを精度よく求めることができる。そして、これらの関係性を用いて、各坑井２の制御可能範囲内における各坑口圧力の組合せに対応した複数の運定点を求めることで、各坑井２の制御可能範囲内で取り得る運定点を推定することが可能となる。

また、地熱流体の総流量Ｑと総エンタルピーＨがパラメータとして発電出力Ｌ毎に表された特性を用いることで、所望の地熱流体の総流量状態及び発電出力状態が与えられれば、対応する所望運定点（地熱流体の総流量Ｑと総エンタルピーＨの組合せ）を特定することが可能となる。各坑井２の坑口圧力を所望運定点に対応した坑口圧力（Ｐ１ｓｅｔ、Ｐ２ｓｅｔ）に制御することで、実際の地熱発電プラント１の運転状態をより最適化することが可能となる。

また、得点付けによる評価を行い、総流量Ｑが最小値に近づき、かつ発電出力Ｌが最大値に近づくような所望運定点Ｗ１を特定することができるため、各坑井２の短寿命化を抑制し、効率的に発電の高出力化を図ることができる。

〔第２実施形態〕
次に、本開示の第２実施形態に係る地熱発電プラントの坑井制御システム及びその坑井制御方法並びに坑井制御プログラム、地熱発電プラントについて説明する。
上述した第１実施形態では、所望運定点特定部２２は、総流量Ｑが最小値に近づき、かつ、発電出力Ｌが最大値に近づくように所望運定点を特定する場合について説明していたが、本実施形態では、所望運定点特定部２２は、総流量Ｑが最小値に近づき、かつ、発電出力Ｌが所定値に近づくように所望運定点を特定する。以下、本実施形態に係る地熱発電プラントの坑井制御システム及びその坑井制御方法並びに坑井制御プログラム、地熱発電プラントについて、第１実施形態と異なる点について主に説明する。

本実施形態における所望運定点特定部２２は、各坑井２の制御可能範囲において、総流量Ｑが最小値に近づき、かつ、発電出力Ｌが所定値に近づくように所望運定点を特定する。地熱発電プラント１では、要求発電出力等によって、発電出力が指定される場合がある。このため、要求発電出力等によって指定される発電出力Ｌを所定値として、所定値に近づくように所望運定点を特定する。

所望運定点特定部２２は、所望の運定点を特定可能なように、得点付けを行う。本実施形態では、総流量Ｑが最小値に近づき、かつ、発電出力Ｌが所定値に近づくように運定点を特定する。このため、所望運定点特定部２２は、図１０に示すように、プロットした運定点群のうち総流量Ｑが最小値となる運定点を最も高得点となるように設定し、総流量Ｑが多くなるにつれて所定の間隔で得点を低くして設定する。図１０の例では、プロットした運定点群のうち総流量Ｑが最小値となる運定点に最高得点として１０点を付け、総流量Ｑが多くなるにつれて所定の間隔で９点、８点等と点数が低下するように付ける。また、所望運定点特定部２２は、プロットした運定点群のうち発電出力Ｌが所定値となる運定点を最も高得点となるように設定し、発電出力Ｌが所定値から離れるにつれて所定の間隔で得点を低くして設定する。図１０の例では、プロットした運定点群のうち発電出力Ｌが所定値となる運定点に最高得点として１０点を付け、発電出力Ｌが所定値から離れるにつれて所定の間隔で９点等と点数が低下するように付ける。

図１０のように得点を付与した場合、総流量が最小値へ近づくほど高得点となるように付与された得点と、発電出力Ｌが所定値へ近づくほど高得点となるように付与された得点とを合算して評価する。例えば丸囲み×印で示したＷ２が最も高得点となる。このため、所望運定点特定部２２は、総流量Ｑが最小値に近づき、かつ、発電出力Ｌが所定値に近づく最適な運定点としてＷ２を特定する。このようにして、所望運定点特定部２２は、最適となる所望の運定点を特定する。そして、坑口圧力設定部２３では、所望運定点に対応する各坑口圧力（Ｐ１ｓｅｔ‘、Ｐ２ｓｅｔ‘）を設定する。

以上説明したように、本実施形態に係る地熱発電プラントの坑井制御システム及びその坑井制御方法並びに坑井制御プログラム、地熱発電プラントによれば、総流量Ｑが最小値に近づき、かつ発電出力Ｌが所定値に近づくような所望運定点を特定することができるため、各坑井２の短寿命化を抑制し、効率的に発電出力Ｌを要求発電出力等によって指定される所定値に維持することができる。

〔第３実施形態〕
次に、本開示の第３実施形態に係る地熱発電プラントの坑井制御システム及びその坑井制御方法並びに坑井制御プログラム、地熱発電プラントについて説明する。
上述した第２実施形態では、所望運定点特定部２２は、総流量Ｑが最小値に近づき、かつ、発電出力Ｌが所定値に近づくように所望運定点を特定する場合について説明していたが、本実施形態では、所望運定点特定部２２は、各坑井２の制御可能範囲において、総流量Ｑが最小値に近づき、かつ、発電出力Ｌが所定値以上の範囲において所定値に近づくように所望運定点を特定する。以下、本実施形態に係る地熱発電プラント１の坑井制御システム及びその坑井制御方法並びに坑井制御プログラム、地熱発電プラント１について、第１実施形態及び第２実施形態と異なる点について主に説明する。

本実施形態における所望運定点特定部２２は、各坑井２の制御可能範囲において、総流量Ｑが最小値に近づき、かつ、発電出力Ｌが所定値以上の範囲において所定値に近づくように所望運定点を特定する。例えば、地熱発電プラント１では、接続されている系統側の事故等によって、所内単独運転を行う場合がある。所内単独運転とは、地熱発電プラント１内における各機器への電力供給を自ら発電電力でまかなう状態である。すなわち、地熱発電プラント１内において最低限必要な電力が設定されており、発電出力Ｌがこの最低限必要な電力を下回ると、地熱発電プラント１内における各機器が止まりブラックアウトとなる。このため、地熱発電プラント１では、最低限必要な電力（必要蒸気量）に対して所定の余裕度（マージン）を加算した電力値を発電出力Ｌの所定値として、地熱流体の生産を行う。一方で、実際に発電に必要な地熱流体に対する過剰分は系外排出される等の処理がなされるため地熱流体の生産量はなるべく低い方が好ましい。そこで、所望運定点特定部２２は、総流量Ｑが最小値に近づき、かつ、発電出力Ｌが所定値以上の範囲において所定値に近づくように所望の運定点を特定する。なお、所定値を下回る範囲では、地熱発電プラント１がブラックアウトとなってしまうため、選定不可領域と設定する。また、蒸気タービン８に供給する総蒸気流量Ｑｓｓのうち一部を、蒸気タービン８の入口側において設けた放出管（蒸気管７から蒸気管７を流通する蒸気の一部を大気に放出する管）によって放出し、発電出力Ｌを調整することとしてもよい。この場合には、放出管に設けた放出弁によって、放出する蒸気の量を調節することとすればよい。

所望運定点特定部２２は、所望の運定点を特定可能なように、得点付けを行う。本実施形態では、総流量Ｑが最小値に近づき、かつ、発電出力Ｌが所定値以上の範囲において所定値に近づくように運定点を特定する。このため、所望運定点特定部２２は、図１１に示すように、プロットした運定点群のうち総流量Ｑが最小値となる運定点を最も高得点となるように設定し、総流量Ｑが多くなるにつれて所定の間隔で得点を低くして設定する。図１１の例では、プロットした運定点群のうち総流量Ｑが最小値となる運定点に最高得点として１０点を付け、総流量が多くなるにつれて所定の間隔で９点、８点等と点数が低下するように付ける。また、所望運定点特定部２２は、プロットした運定点群のうち発電出力Ｌが所定値となる運定点を最も高得点となるように設定し、発電出力Ｌが所定値以上の範囲において増加するにつれて所定の間隔で得点を低くして設定する。なお、発電出力Ｌが所定値を下回る範囲において選定されないように例えば得点を極めて低く設定する。図１１の例では、プロットした運定点群のうち発電出力Ｌが所定値となる運定点に最高得点として１０点を付け、発電出力Ｌが所定値から増加するにつれて所定の間隔で９点等と点数が低下するように付ける。なお、発電出力Ｌが所定値を下回る運定点には選定不可領域として、例えば０点やマイナス等の点数を付けるか、もしくはそのような運定点は棄却されるようにする。

図１１のように得点を付与した場合、総流量が最小値へ近づくほど高得点となるように付与された得点と、発電出力Ｌが所定値以上の範囲において所定値に近づくほど高得点となるように付与された得点とを合算して評価する。例えば丸囲み×印で示したＷ３が最も高得点となる。このため、所望運定点特定部２２は、総流量Ｑが最小値に近づき、かつ、発電出力Ｌが所定値以上の範囲において所定値に近づく最適な運定点としてＷ３を特定する。このようにして、所望運定点特定部２２は、最適となる所望の運定点を特定する。そして、坑口圧力設定部２３では、所望運定点に対応する各坑口圧力（Ｐ１ｓｅｔ‘‘、Ｐ２ｓｅｔ‘‘）を設定する。

以上説明したように、本実施形態に係る地熱発電プラントの坑井制御システム及びその坑井制御方法並びに坑井制御プログラム、地熱発電プラントによれば、総流量Ｑが最小値に近づき、かつ発電出力Ｌが所定値以上の範囲において該所定値に近づくような所望運定点を特定することができるため、各坑井２の短寿命化を抑制し、発電出力Ｌが最低限必要な電力である所定値以上としながら、余剰な発電を抑えた運転を行うことができ、地熱発電プラント１のブラックアウトを防止することが可能となる。

〔第４実施形態〕
次に、本開示の第４実施形態に係る地熱発電プラントの坑井制御システム及びその坑井制御方法並びに坑井制御プログラム、地熱発電プラントについて説明する。
本実施形態では、所望運定点特定部２２は、各坑井２における地熱流体の流量Ｑｉ及び噴気圧力Ｐｉの少なくともいずれか一方の設定範囲を限定するように設定した坑井評価関数を用いて、各坑井２での制御可能範囲（流量調整弁４を用いて制御可能な範囲）を更に制限することで所望運定点を限定したものにおいて、第１実施形態から第３実施形態で説明した所望運定点を特定する。
本実施形態では、一例として坑井評価関数を用いて、各坑井２での制御可能範囲を更に制限することで所望運定点を限定したものにおいて、第１実施形態で示した、総流量Ｑが最小値に近づき、発電出力Ｌが最大値に近づくように所望運定点を特定するものを説明する。以下、本実施形態に係る地熱発電プラントの坑井制御システム及びその坑井制御方法並びに坑井制御プログラム、地熱発電プラントについて、第１実施形態と異なる点について主に説明する。

本実施形態における所望運定点特定部２２は、各坑井２での制御可能範囲を更に制限するものにおいて、総流量Ｑが最小値に近づき、かつ、発電出力Ｌが最大値に近づくような所望運定点を特定する。各坑井２での制御可能範囲を更に制限する理由として、各坑井２の噴気圧力Ｐｉを低くして流量Ｑｉが大きくなった場合に坑井２の過熱化や枯渇減衰が大きくなる場合がある。このため、各坑井２での坑井特性に対応した流量Ｑｉ及び噴気圧力Ｐｉの少なくともいずれか一方の範囲を設定する。流量Ｑｉ及び噴気圧力Ｐｉの少なくともいずれか一方の範囲を適切に設定しないと（制御可能範囲外として考慮すべき状況としては）、例えば、制御可能範囲の流量を超える多量の地熱流体を噴気し続けた場合には、地熱貯留層に貯留されている地下流体の水源が乏しい場合などでは地熱流体の減少を促進し、坑井２の寿命を余分に短くしてしまうおそれがある。また、例えば、制御可能範囲の噴気圧力を下回る圧力の地熱流体を噴気し続けた場合には、複数の坑井２のフィードポイント間のバランスが変化して各坑井２の噴気特性が変化するおそれがある。また、例えば、制御可能範囲の噴気圧力を超える圧力の地熱流体を噴気し続けた場合には、地熱流体の噴気が停止するおそれがある。これら状況を防止するために、各坑井２の運用を制御可能範囲にて制限する。すなわち、所望運定点を特定する際には各坑井２での制御可能範囲外の運転条件が含まれる運定点が選択されないように、各坑井２での流量Ｑｉ及び噴気圧力Ｐｉの少なくともいずれか一方に対応した坑井評価関数を与える。

図１２は一例として、各坑井２のうちの１つの坑井２の噴気圧力Ｐｉに対応して設定された坑井評価関数を示す図である。なお、坑井評価関数は、地熱流体の流量Ｑｉに対応して設定されてもよい。図１２では、初期の噴気特性試験等で得られた坑井特性やトレーサ試験（貯留層の流動調査試験）の結果、または坑井２の実運用状況などに基づいて、噴気圧力Ｐｉに対する制御可能範囲を設定可能範囲として設け、設定可能範囲の下限値Ｐｉａと上限値Ｐｉｄが設定される。また、坑井の寿命に関する経験的な情報や実運用状況などに基づいて、噴気圧力Ｐｉに対する制御可能範囲をさらに制限した（設定可能範囲の下限値Ｐｉａと上限値Ｐｉｄとの間の範囲に上下限値が設定された）設定推奨範囲を設け、設定推奨範囲の下限値Ｐｉｂと上限値Ｐｉｃが設定される。本実施形態では坑井評価関数は、例えば設定推奨範囲の下限値Ｐｉｂと上限値Ｐｉｃの間で１とし、制御可能範囲内の下限値ＰｉａとＰｉｂの間、及びＰｉｃと上限値Ｐｉｄの間では０から１の間の正の値とし、制御可能範囲外である下限値Ｐｉａより低い圧力、及び上限値Ｐｉｄより大きな圧力では０とする。
坑井評価関数で、噴気圧力Ｐｉを設定推奨範囲の下限値Ｐｉｂと上限値Ｐｉｃの間に選定することで、坑井評価関数は１となり、各坑井２の運用を続けた際に坑井特性に大きな影響が発生しないことを示している。また、噴気圧力Ｐｉを設定推奨範囲のＰｉｂからＰｉａへと噴気圧力が低下する間、及び設定推奨範囲のＰｉｃからＰｉｄへと噴気圧力が上昇する間では、坑井評価関数は１から０へと漸次減少し、各坑井２の運用を続けた際に坑井特性に大きな影響が発生し易くなることを示している。

所望運定点特定部２２は、所望の運定点を特定可能なように、得点付けを行う。本実施形態では、各坑井２における地熱流体の流量Ｑｉ及び噴気圧力Ｐｉの少なくともいずれか一方に基づいて、各坑井２での制御可能範囲を更に制限する範囲内において、総流量Ｑが最小値に近づき、かつ、発電出力Ｌが最大値に近づくような所望運定点を特定する。
例えば、図１２で示すように、噴気圧力Ｐｉが設定推奨範囲の下限値Ｐｉｂと上限値Ｐｉｃの間にあり、総流量Ｑが最小値に近づき、かつ、発電出力Ｌが最大値に近づくように所望の運定点を特定する。このため、所望運定点特定部２２は、まず図１３に示すように、プロットした運定点群のうち、噴気圧力Ｐｉが設定推奨範囲にあるかどうかを設定する坑井評価関数により得点付けをする。坑井評価関数による得点付けは、噴気圧力Ｐｉが設定推奨範囲の下限値Ｐｉｂと上限値Ｐｉｃの間となる運定点を最も高得点の１となるように設定し、噴気圧力Ｐｉが小さくなりＰｉｂとＰｉａの間では１から０へと連続的に得点が低くなるように設定し、更に噴気圧力Ｐｉが小さくなりＰｉａよりも低い圧力では最低点の０となるよう設定する。また、噴気圧力Ｐｉが大きくなりＰｉｃとＰｉｄの間では１から０へと連続的に得点が低くなるように設定し、更に噴気圧力Ｐｉが大きくなりＰｉｄよりも高い圧力では最低点の０となるよう設定する。また、所望運定点特定部２２は、プロットした運定点群のうち、総流量Ｑが最小値となる運定点を最も高得点として１０点となるように設定し、総流量Ｑが多くなるにつれて所定の間隔で得点を低くして設定する。また、所望運定点特定部２２は、プロットした運定点群のうち発電出力Ｌが最大値となる運定点を最も高得点として１０点となるように設定し、発電出力Ｌが少なくなるにつれて所定の間隔で得点を低くして設定する。

図１３のように得点を付与した場合、総流量が最小値へ近づくほど高得点となるように付与した得点と、発電出力Ｌが最大値へ近づくほど高得点となるように付与した得点を合算して、合算した値に対して噴気圧力Ｐｉの選定による坑井評価関数による得点を積算することによって最終得点を算出し評価する。例えば図１３の例では、丸囲み×印で示したＷ４が最も高得点となる。このため、所望運定点特定部２２は、各坑井２での制御可能範囲を更に制限する範囲内において、総流量Ｑが最小値に近づき、かつ、発電出力Ｌが最大値に近づく最適となる運定点としてＷ４を特定する。なお、本実施形態では、丸囲み×印で示したＷ４の最適点の選定にあたっては、Ｗ４の近傍の運定点に対して、坑井評価関数の得点を追加評価することによって、Ｗ４の最適点を選定している。このようにして、所望運定点特定部２２は、最適となる所望の運定点を特定する。また、設定推奨範囲の下限値Ｐｉｂと上限値Ｐｉｃは、各坑井の特性から各坑井毎での経験的に基づいて設定されてもよい。噴気圧力Ｐｉを設定可能範囲のＰｉｂから下限値Ｐｉａへと噴気圧力が低下する間隔、及び設定可能範囲のＰｉｃから上限値Ｐｉｄへと噴気圧力が上昇する間隔は、設定推奨範囲とは無関係に狭くしてもよく、または実質的に間隔を無しとして設定推奨範囲をより明確に限定してもよい。

このように、各坑井２における地熱流体の流量Ｑｉ及び噴気圧力Ｐｉの少なくともいずれか一方の制御可能範囲に対応した設定可能範囲と、坑井の寿命等を加味して設定可能範囲を限定して設定された設定推奨範囲とが設けられた坑井評価関数を用いることによって、設定推奨範囲に含まれる運定点を高く評価することができる。なお、坑井評価関数において、設定可能範囲に対して設定推奨範囲の方が評価が高くなるように設定されている。設定推奨範囲に含まれる運定点が高評価されやすくなることによって、坑井の寿命等を加味してより最適な運定点を特定することが可能となる。

以上説明したように、本実施形態に係る地熱発電プラントの坑井制御システム及びその坑井制御方法並びに坑井制御プログラム、地熱発電プラントによれば、坑井の特性変化や噴気停止のおそれのある運定点が選定されることを防止できるため、各坑井２の短寿命化をより一層に抑制することができる。

また、制御可能範囲内において、安全性を確保することが可能な噴気圧力Ｐｉが下限値Ｐｉｂと上限値Ｐｉｃの間からなる制御推奨範囲を設定してあるので、制御推奨範囲内における坑井評価関数は、各坑井２の短寿命化をより一層に抑制できる適切な運定点が選定されやすくなり、地熱発電プラント１の信頼性を向上することができる。

本開示は、上述の実施形態のみに限定されるものではなく、発明の要旨を逸脱しない範囲において、種々変形実施が可能である。なお、各実施形態を組み合わせることも可能である。

１ ：地熱発電プラント
２ ：坑井
３ ：地熱流体輸送管
４ ：流量調整弁
５ ：気水分離器（分離器）
６ ：熱水管
７ ：蒸気管
８ ：蒸気タービン
９ ：発電機
１０ ：復水器
１１ ：冷却塔
１２ ：ポンプ
１３ ：開閉弁
１４ ：ポンプ
１５ ：流量調整弁
１６ ：圧力計
１７ ：流量計
１８ ：流量計
１９ ：還元井
２０ ：坑井制御システム
２１ ：運定点設定部
２２ ：所望運定点特定部
２３ ：坑口圧力設定部

Claims (10)

1. 予め設定される過去所定期間である指定期間において、複数の坑井から噴出した地熱流体の総流量及び総エンタルピーと、各前記坑井の坑口圧力とに基づいた運転データを蓄積し、
   蓄積した前記運転データに基づいて前記総流量の変動量及び前記総エンタルピーの変動量を算出し、
   各前記坑口圧力の基準圧力に対する偏差と前記総流量の前記変動量との関係性及び前記偏差と前記総エンタルピーの前記変動量との関係性に基づいて、各前記坑井における前記坑口圧力の組合せに対応した運定点を複数設定する運定点設定部を備える地熱発電プラントの坑井制御システム。
2. 前記地熱流体の前記総流量と前記総エンタルピーとの関係に対して、前記地熱流体を用いた発電出力に係る情報を対応づけた特性に基づいて、所望の前記地熱流体の総流量状態及び発電出力状態となる前記運定点を所望運定点として特定する所望運定点特定部と、
   前記所望運定点に対応する各前記坑口圧力を設定する坑口圧力設定部と、
   を備える請求項１に記載の地熱発電プラントの坑井制御システム。
3. 前記特性は、前記地熱流体の前記総流量と前記総エンタルピーとの関係が前記地熱流体より生成された蒸気流量に基づいて表されており、前記発電出力が前記蒸気流量に対応して表されている請求項２に記載の地熱発電プラントの坑井制御システム。
4. 前記所望運定点特定部は、前記総流量が最小値に近づき、かつ、前記発電出力が最大値に近づくように前記所望運定点を特定する請求項２または３に記載の地熱発電プラントの坑井制御システム。
5. 前記所望運定点特定部は、前記総流量が最小値に近づき、かつ、前記発電出力が所定値に近づくように前記所望運定点を特定する請求項２または３に記載の地熱発電プラントの坑井制御システム。
6. 前記所望運定点特定部は、前記総流量が最小値に近づき、かつ、前記発電出力が所定値以上の範囲において前記所定値に近づくように前記所望運定点を特定する請求項２または３に記載の地熱発電プラントの坑井制御システム。
7. 各前記坑井における前記地熱流体の流量及び前記坑井の坑口圧力の少なくともいずれか一方の設定範囲を限定するように設定した坑井評価関数を用いて前記所望運定点を評価し、前記所望運定点を特定する請求項４から６のいずれか１項に記載の地熱発電プラントの坑井制御システム。
8. 坑井から噴出された地熱流体を蒸気と熱水に分離する分離器と、
   前記分離器により分離された蒸気を用いて発電を行う蒸気タービンと、
   請求項１から７のいずれか１項に記載の地熱発電プラントの坑井制御システムと、
   を備えた地熱発電プラント。
9. 予め設定される過去所定期間である指定期間において、複数の坑井から噴出した地熱流体の総流量及び総エンタルピーと、各前記坑井の坑口圧力とに基づいた運転データを蓄積し、
   蓄積した前記運転データに基づいて前記総流量の変動量及び前記総エンタルピーの変動量を算出し、
   各前記坑口圧力の基準圧力に対する偏差と前記総流量の前記変動量との関係性及び前記偏差と前記総エンタルピーの前記変動量との関係性に基づいて、各前記坑井における前記坑口圧力の組合せに対応した運定点を複数設定する運定点設定工程を有する地熱発電プラントの坑井制御方法。
10. 予め設定される過去所定期間である指定期間において、複数の坑井から噴出した地熱流体の総流量及び総エンタルピーと、各前記坑井の坑口圧力とに基づいた運転データを蓄積し、
    蓄積した前記運転データに基づいて前記総流量の変動量及び前記総エンタルピーの変動量を算出し、
    各前記坑口圧力の基準圧力に対する偏差と前記総流量の前記変動量との関係性及び前記偏差と前記総エンタルピーの前記変動量との関係性に基づいて、各前記坑井における前記坑口圧力の組合せに対応した運定点を複数設定する運定点設定処理をコンピュータに実行させるための地熱発電プラントの坑井制御プログラム。

Images