【地熱】バイナリー地熱発電における発電量の簡易計算方法

ここでは、バイナリー地熱発電における簡易的な出力の計算方法について、解説します。
フラッシュ発電の場合は、今回の計算と異なるため除外します。

バイナリー地熱発電とは?

フロンやアンモニアなど沸点が低い媒体を使用して、タービンを回して発電するシステムです。別名、温泉発電とも言われます。

バイナリー発電システムの構成例

図 バイナリー発電システムの構成例
(経済産業省HPから引用())

特徴としては、比較的低温流体(<120℃)で使用できる、『高温”熱水”で使用できる』ことが挙げられます。
通常のフラッシュ発電では高温高圧の蒸気を必要としますが、バイナリー地熱発電では温泉のような高温の熱水でも活用できます。

水の状態図とエンタルピー

後日更新予定です。他のサイトでも十分理解できるので、リンクを載せておきます。

株式会社大西熱学 顕熱と潜熱

株式会社ミヤワキ 蒸気の基本

蓄熱WEB講座PRO 3-2「潜熱」「顕熱」について

計算方法

この計算は簡易的な計算方法となるため、出力を保証するものではありません。詳細については、別途メーカーに確認してください。

バイナリー発電機の発電効率は、機器や温度条件等によって異なります。地熱流体の温度が高いほどエントロピーが増加し、発電効率は上昇します。また、冷却能力が高いとタービン前後の差圧が上がるため発電機効率が上昇します。諸々踏まえても、一般的に8~15%程度となります。ここでは、E(%)として計算します。

今回は、以下の3パターンに分けて、計算方法を紹介しています。

飽和蒸気単層
二相流
熱水単層

ただし、発電機出口は熱水単層、発電機入口は飽和条件と仮定します。また、計算方法についても2種類用意しました。理解しやすい方で理解して頂ければ幸いです。

過熱蒸気は、計算方法が異なるため注意が必要です。

飽和蒸気のみの場合

バイナリー発電

図 バイナリー発電のシステム図(蒸気ver)

計算方法① 潜熱と顕熱を別に計算

飽和蒸気のみの場合は、「蒸気の潜熱」「熱水になった後の顕熱」の足し算になります。
蒸気流量Q (kg/s)、入口温度T1(℃)、出口温度T2(℃)、温度T1のときの蒸気潜熱h(kJ/kg)とすると、以下の計算式です。熱水量も蒸気量と同じためQ(kg/s)と置けます。

\[ 発電出力(kw) = \bigl\{h + 4.2 \times(T1-T2)\bigr\}\times Q \times E \]

単位計算とすると、以下のようになります。

\[ kw = \{\frac{kJ}{kg} + \frac{kJ}{kg·℃} \times ℃ \}\times \frac{kg}{s} \times 発電効率%\]

また、熱水の比熱は厳密には温度や圧力条件によって変化しますが、ここでは簡易的に4.2(kJ/kg·℃)で一定としています。

計算方法② 比エンタルピーを用いた計算

比エンタルピーを用いた場合、入口の飽和蒸気の比エンタルピーをh(kJ/kg)、出口の熱水の比エンタルピーをh'(kJ/kg)とした、こちらの式です。

\[ 発電出力(kw) = (h-h’)\times Q \times E\]

二相流の場合

バイナリー発電

図 バイナリー発電のシステム図(二相流ver)

二相流は、『蒸気と熱水が混じった流体』です。蒸気と熱水のエネルギーをそれぞれ計算することになります。厳密には、蒸発器と予熱器の出口温度は違いますが、二相流熱交換器などもあることから出口温度は同一であるとします。

計算方法① 蒸気潜熱と熱水顕熱を別に計算

蒸気流量Q1(kg/s)、熱水流量Q2(kg/s)、入口温度T1(℃)、出口温度T2(℃)、温度T1のときの蒸気潜熱h(kJ/kg)とすると、以下の式になります。

\[ 発電出力(kw) = \Bigl[ h\times Q1 + \bigl\{ 4.2 \times(T1-T2) \times (Q1+Q2)\bigr\} \Bigr] \times E \]

単位計算とすると

\[ kw = [\frac{kJ}{kg} \times \frac{kg}{s} + \{\frac{kJ}{kg·℃} \times ℃ \times \frac{kg}{s} \}] \times 発電効率%\]

先ほどと同様に熱水の比熱は、簡易的に4.2(kJ/kg·℃)で一定としています。

計算方法② 比エンタルピーを用いた計算

別の計算方法として、乾き度比エンタルピーを用いて簡単に計算する方法もあります。まずは、乾き度と比エンタルピーを求めます。

乾き度は、流量全体のうちの蒸気量です。式は以下の通りです。

\[ 乾き度X = \frac{Q1}{(Q1+Q2)} \]

飽和蒸気の比エンタルピーをhG、飽和熱水の比エンタルピーをhLとしたとき、二相流の比エンタルピーhは以下のように計算します。

\[ h(kJ/kg) = X \times h_G +(1-X)\times h_L \]

以上の結果より、出口の熱水の比エンタルピーをh'(kJ/kg)としたとき、以下のように発電量を計算します。

\[ 発電出力(kw) = (h-h’)\times Q \times E\]

二相流の比エンタルピーについては、私が作成した以下のサイトで飽和圧力と乾き度から自動計算することも可能ですので、ご利用ください。

蒸気・熱水二相流の物性値

熱水のみの場合

バイナリー発電

図 バイナリー発電のシステム図(熱水ver)

熱水の場合は、熱水流量Q (kg/s)、入口温度T1(℃)、出口温度T2(℃)としたとき、以下の計算式になります。

\[ 発電出力(kw) = 4.2 \times(T1-T2) \times Q \times E \]

先ほどと同様に熱水の比熱は、簡易的に4.2(kJ/kg·℃)で一定としています。

計算例

例)総流量 10t/h乾き度 0.3坑口圧力 0.5MPaGの地熱流体が噴出しているとします。このとき、想定される発電出力はいくつになるでしょうか?また、このとき発電効率Eは12%、出口の熱水温度T2は90℃とします。

計算方法①の場合

まず、蒸気流量Q1と熱水流量Q2をそれぞれ求めます。

\[ Q1(t/h) = 10 \times 0.3 = 3 \] \[ Q1(kg/s) = 3 \times 1000 \div 3600 = 0.83\]
\[ Q2(t/h) = 10 \times 0.7 = 7 \] \[ Q2(kg/s) = 7 \times 1000 \div 3600 = 1.94\]

次に蒸気表より潜熱h(kJ/kg)、飽和温度T(℃)について調べます。

\[ 潜熱 = 飽和蒸気比エンタルピー - 飽和熱水比エンタルピー\] \[ h(kJ/kg) = 2756.23 - 670.876 = 2085.36\] \[ 温度T1(℃) = 158.919 \]

以上より、発電量は以下のようになります。

\[ 発電出力(kw) = \Bigl[ h\times Q1 + \bigl\{ 4.2 \times(T1-T2) \times (Q1+Q2)\bigr\} \Bigr] \times E \] \[ = \Bigl[ 2085.36 \times 0.83 + \bigl\{ 4.2 \times(158.9-90) \times (2.77)\bigr\} \Bigr] \times 12% \] \[ \doteqdot (1730.8 + 801.6) \times 12% \doteqdot 303.9\]

計算方法②の場合

まず、二相流の比エンタルピーを導出します。

\[ h(kJ/kg) = X \times h1 +(1-X)\times h2 \] \[ = 0.3 \times 2756.23 +(1-0.3)\times 670.876 \] \[ =  1296.483 \]

飽和熱水の比エンタルピーは蒸気表より以下します。

\[ h'(kJ/kg) = 376.968 \]

以上より、発電量を算出すると

\[ 発電出力(kw) = (h-h’)\times Q \times E\] \[   = (1296.483 -376.968 )\times 2.77 \times 12%\] \[   = 305.6 \]

熱水の比エンタルピーを飽和蒸気表を使ったこと、所々で数字を丸めたことでで、①の計算と若干数字がずれていますが、おおよそ同じ値となりました。

自動計算サイト

自動計算サイトを作成しました。ご活用ください。
ただし、出力を保証するものではございませんので、ご利用の際には十分お気を付けください。

https://www.isenthalpiccalc.com/geo/binary

関連記事

最新情報をチェックしよう!