ペロブスカイト太陽光発電システムを数式で理解する：基礎・発電量計算・係数・シミュレーション・API・データセット完全ガイド

想定読者:
- 太陽光発電・蓄電池・PPA・GXサービスの開発者
- pvlib / SAM / PVGIS / Solar Pro / PVsyst などで発電量シミュレーションを扱うエンジニア
- ペロブスカイト太陽電池を「ニュース」ではなく「計算モデル」として理解したい人
- 再エネSaaS / API / BPO / 提案自動化に関わるPdM・事業開発・技術責任者

タグ候補:
perovskite, solar, pvlib, renewable-energy, simulation, python, photovoltaic, energy

TL;DR：ペロブスカイト太陽光の本質は「高効率」よりも「設置制約を破る計算対象」である

ペロブスカイト太陽電池は、単に「次世代の高効率太陽電池」ではない。
工学的には、軽量・薄膜・曲面・壁面・低耐荷重屋根・建材一体型という、従来の結晶シリコン太陽光が入りにくかった場所を発電面に変える技術である。

ただし、シミュレーションの難易度は上がる。

従来の屋根置き太陽光では、概算なら「日射量 × 容量 × パフォーマンスレシオ」でかなり見通せた。ペロブスカイトでは、そこに次の論点が強く入る。

軽量・薄膜ゆえの温度挙動
フィルム・壁面・曲面設置による入射角損失
都市部・壁面・低高度日射で増える部分影・反射・低照度特性
湿気・熱・紫外線・酸素・イオン移動による劣化モデル
単接合か、シリコンとのタンデムかによる電気モデルの違い
実証段階から量産初期へ移る過程での保証・認証・バンカビリティ

米国DOEは、ペロブスカイトの主な課題として、安定性・耐久性、モジュールスケール効率、製造性、実証・バンカビリティを挙げている。これは、発電量計算でもまったく同じ順番で効いてくる。研究セルの変換効率だけを見ても、事業性は決まらない。(The Department of Energy's Energy.gov)

1. ペロブスカイト太陽電池とは何か

ペロブスカイト太陽電池、特に金属ハライドペロブスカイト太陽電池は、ABX3 型の結晶構造を持つ材料を光吸収層に使う太陽電池である。代表例として、MAPbI3 のような鉛ハライド系材料が知られている。DOEの解説では、ペロブスカイトは結晶シリコンと比べて非常に薄い活性層で太陽光を吸収でき、溶液塗布や蒸着など多様な製造プロセスに対応し得る技術として整理されている。(The Department of Energy's Energy.gov)

1.1 基本構造

典型的な単接合ペロブスカイト太陽電池は、概念的には次の層で構成される。

光入射側
  ↓
透明電極
電子輸送層 ETL
ペロブスカイト光吸収層
正孔輸送層 HTL
金属電極
封止層

実際の構成は、n-i-p構造、p-i-n構造、カーボン電極型、メソポーラス型、フレキシブル基板型、ガラス基板型などに分かれる。

エンジニア視点で重要なのは、材料名よりも次の4点である。

観点	シミュレーション上の意味
光吸収層が薄い	入射角・反射・低照度特性の影響を受けやすい場合がある
バンドギャップを調整しやすい	タンデム構成に向く
塗布・印刷・蒸着に対応し得る	大面積化・歩留まり・面内ばらつきが重要
湿気・熱・酸素・光で劣化しやすい材料系がある	劣化率をシリコンと同じに置くと危険

ペロブスカイトの強みは、セル単体の効率だけではない。
「今まで発電面として扱えなかった場所」を計算対象に入れることにある。

2. 現在地：研究効率、国内政策、量産初期の実装状況

2026年時点で、ペロブスカイトは研究開発だけの技術ではなく、実証・商用初期・量産準備へ移行している。

NRELのBest Research-Cell Efficienciesチャートでは、ペロブスカイト単接合およびペロブスカイト/シリコンタンデムの研究セル効率が高い水準に到達しており、ハイブリッドタンデム領域では30%台半ばの記録が示されている。研究セル効率という意味では、すでに「高効率候補」ではなく「高効率領域に入った技術」と見てよい。(NREL)

日本では、次世代型太陽電池戦略の中で、2040年に約20GWという導入目標、発電コスト目標として2025年20円/kWh、2030年14円/kWh、2040年10〜14円/kWh以下といった水準が示されている。

NEDOのグリーンイノベーション基金では、2030年度までにペロブスカイト単接合で14円/kWh以下、タンデム型で住宅用12円/kWh、変換効率30%という目標が掲げられている。(NEDO グリーンイノベーション基金)

国内メーカーでは、積水化学がフィルム型ペロブスカイト太陽電池「SOLAFIL」関連の事業開始を公表し、2026年度は限定生産、2027年度に100MW級製造ラインを優先する方針を示している。公開資料では、100MW級ライン、2027年稼働計画、2030年のGW級製造も言及されている。(積水化学株式会社)

また、初期モジュール仕様例として、最大出力110W、寸法約1,450mm × 985mm × 30mm、重量6.4kg/枚、面密度4.5kg/m²のフィルム型モジュールが示されている。単純計算すると、この仕様例の面積当たり定格効率は約7.7%である。ただしこれは特定初期製品仕様からの導出値であり、ペロブスカイト技術全体の効率を代表する値ではない。

area_m2 = 1.45 * 0.985
p_stc_w = 110
eta = p_stc_w / (1000 * area_m2)
print(round(eta * 100, 2))  # 約7.7%

ここで重要なのは、研究セル効率と量産初期モジュール効率を混同しないことだ。

研究セル効率:
  材料・デバイスの到達可能性を見る指標

量産モジュール効率:
  面積、封止、配線、歩留まり、耐久性、安全性を含む製品指標

システム発電量:
  設置面、日射、温度、影、入射角、PCS、劣化、運用を含む事業指標

この3つを混ぜると、提案も投資判断も壊れる。

3. ペロブスカイトの発電原理を数式で見る

太陽電池の基本は、光で発生したキャリアを電流として取り出すことである。
セルの変換効率は、標準試験条件では次のように表せる。

η = P_{max} / (G_{STC} × A)

また、電気的には次式でも表せる。

η = (V_{oc} × J_{sc} × FF) / G_{STC}

記号	意味
`η`	変換効率
`Pmax`	最大出力
`G_STC`	標準日射強度。通常1000 W/m²
`A`	セルまたはモジュール面積
`Voc`	開放電圧
`Jsc`	短絡電流密度
`FF`	フィルファクター

フィルファクターは次で定義される。

FF = (V_{mp} × I_{mp}) / (V_{oc} × I_{sc})

3.1 単接合セルの理想と現実

単接合セルでは、1つの吸収層が広い波長範囲の光を受ける。
しかし、バンドギャップより低いエネルギーの光は吸収されず、バンドギャップより大きいエネルギーの余剰分は主に熱として失われる。

この限界を緩める発想がタンデムである。

3.2 タンデムセルの基本式

ペロブスカイト/シリコンタンデムでは、ペロブスカイトを上側セル、シリコンを下側セルとして使うことが多い。

2端子タンデム

2端子タンデムでは、上側セルと下側セルが直列接続される。したがって電流は小さい方に制約される。

J_{tandem} = min(J_{top}, J_{bottom})

電圧は足し合わされる。

V_{tandem} = V_{top} + V_{bottom}

したがって概念的には、

P_{tandem} ≈ J_{match} × (V_{top} + V_{bottom}) × FF

となる。

2端子タンデムのシミュレーションでは、電流マッチングが支配的になる。
上側ペロブスカイト層のバンドギャップ、厚み、透過率、反射、スペクトル変動が、下側シリコンセルの発電量まで左右する。

4端子タンデム

4端子タンデムでは、上側セルと下側セルを電気的に別々に取り出せる。

P_{total} = P_{top} + P_{bottom}

4端子は電流マッチング制約が緩い一方、光学設計、配線、インバータ、制御、コストが複雑になる。

ソフトウェア実装では、2端子は「結合された1つの特殊モジュール」、4端子は「2つのPVサブシステム」として扱うと整理しやすい。

4. I-Vカーブと単一ダイオードモデル

太陽電池の電気モデルでは、単一ダイオードモデルが基本になる。

I = I_L - I_0 \left[\exp\left(\frac{V + I R_s}{n N_s V_T}\right)-1\right] - \frac{V + I R_s}{R_{sh}}

記号	意味
`I_L`	光生成電流
`I_0`	ダイオード逆飽和電流
`R_s`	直列抵抗
`R_sh`	並列抵抗
`n`	理想係数
`N_s`	直列セル数
`V_T`	熱電圧

pvlibのsinglediode()は、この単一ダイオード方程式を解くために、Lambert W法やBishop法などのアルゴリズムを提供している。(Pvlib Python)

ペロブスカイトで難しいのは、I-V特性が単純な静的モデルに収まりにくい場合があることだ。
イオン移動、ヒステリシス、界面状態、時間依存劣化が効く場合、同じ日射・温度でも、スキャン方向や履歴でI-Vカーブが変わることがある。

高精度モデルでは、状態変数を導入する。

I(t) = f(V(t), G(t), T(t), ξ(t))

\frac{dξ}{dt} = \frac{ξ_{eq}(V,G,T) - ξ}{τ}

ここでξは、イオン分布や界面状態を表す抽象的な内部状態である。
商用シミュレーションではここまで入れないことが多いが、研究・検証用途では重要になる。

5. 発電量シミュレーションの基本パイプライン

ペロブスカイト太陽光発電システムの発電量計算は、次の9段階に分解できる。

1. 位置・時刻
2. 気象データ
3. 水平面日射の分解
4. 傾斜面・壁面・曲面への変換
5. 影・反射・入射角補正
6. モジュール温度計算
7. DC発電量計算
8. インバータ・AC変換・制御
9. 劣化・停止・経済性評価

従来PVと同じに見えるが、ペロブスカイトでは特に4〜6と9が重い。

6. 日射計算：GHI、DNI、DHI、POA

気象データとしてよく使う日射量は次の3つである。

記号	意味
`GHI`	全球水平面日射量
`DNI`	直達法線面日射量
`DHI`	散乱水平面日射量

傾斜面、壁面、曲面で使うのは、最終的にはPOA、つまりPlane of Array irradianceである。

G_{POA} = G_{beam} + G_{sky} + G_{ground}

より具体的には、

G_{POA} = DNI × cos(θ) + DHI × F_{sky} + ρ × GHI × F_{ground}

記号	意味
`θ`	太陽光線とモジュール法線の入射角
`ρ`	地表面反射率、アルベド
`F_sky`	空を見る形状係数
`F_ground`	地面を見る形状係数

JPPAの壁面太陽光発電量推定手法ガイドラインでは、水平面日射を直達・散乱に分離し、ErbsモデルやPerezモデルなどを用いて任意方位・傾斜面の日射量へ変換する流れが整理されている。

6.1 ペロブスカイトで壁面・垂直面が重要になる理由

ペロブスカイトの社会実装では、軽量性を活かした低耐荷重屋根、壁面、曲面、既存建物への後付けが重要になる。NEDOも、日本では平地面積が限られるため、耐荷重の小さい屋根や建物壁面など、従来型太陽電池では導入しにくい場所への適用を重視している。(NEDO グリーンイノベーション基金)

しかし、壁面PVのシミュレーションは難しい。

JPPAのガイドラインでは、壁面PVでは都市部の近傍建物による遮蔽、天空率の低下、高入射角、低照度、反射、ミスマッチ、MPPTの難しさなどに注意が必要とされている。

つまり、ペロブスカイトは「設置できる場所」を増やす。
同時に、「雑な発電量計算では外しやすい場所」を増やす。

ここが反直感ポイントである。

7. 入射角損失：IAMの計算

太陽光がモジュール面に斜めから入ると、反射や光路変化により有効吸収量が落ちる。これを入射角補正、IAM、Incidence Angle Modifierと呼ぶ。

代表的なASHRAEモデルは次で表される。

IAM = 1 - b_0 \left(\frac{1}{cosθ} - 1\right)

JPPAガイドラインでは、このASHRAEモデルに加えて、Martin-Ruizモデルも整理されている。

IAM(θ) =
\frac{1 - exp[-cos(θ)/a_r]}{1 - exp[-1/a_r]}

ペロブスカイトでは、基板、封止、表面フィルム、電極、光取り込み構造が従来シリコンモジュールと異なるため、IAM係数をそのまま流用するのは危険である。

実務では次の順で扱う。

初期検討:
  ASHRAEまたはMartin-Ruizで仮置き

設計検討:
  メーカー実測または角度別透過・反射データで補正

実証後:
  実測発電量からIAM係数を逆推定して更新

8. モジュール温度モデル

太陽電池の出力はセル温度に依存する。
一般にセル温度が上がると出力は下がる。

簡易モデルでは、

P_{dc}(t) =
P_{STC}
× \frac{G_{POA}(t)}{1000}
× [1 + γ × (T_{cell}(t) - 25)]

と書ける。

記号	意味
`P_STC`	標準条件での定格出力
`G_POA`	アレイ面日射
`γ`	出力温度係数
`T_cell`	セル温度

NREL SAMの技術リファレンスでは、Sandia Array Performance Modelの温度式として、次の形が示されている。(NREL)

T_m = G × e^{a + b v_w} + T_a

T_c = T_m + \frac{G}{1000} × ΔT

記号	意味
`T_m`	モジュール温度
`T_c`	セル温度
`G`	POA日射量
`v_w`	風速
`T_a`	外気温
`a, b, ΔT`	設置方式・モジュール構造ごとの係数

SAMの例では、Polymer/Thin Film/SteelのOpen Rack条件で a=-3.58, b=-0.113, ΔT=3 が示されている。(NREL)

ここで注意したい。

フィルム型ペロブスカイトを金属屋根に貼る、壁面に貼る、建材一体で使う、空気層なしで使う、という条件では、温度モデルは大きく変わる。
「薄膜だから温度が低い」と決めつけてもいけないし、「シリコンと同じ係数でよい」とも言えない。

温度モデルは、ペロブスカイト発電量計算の隠れた支配変数である。

9. DC発電量モデル

発電量シミュレーションで使いやすい形にまとめると、DC出力は次のように書ける。

P_{dc}(t) =
P_{STC}
× \frac{G_{POA}(t)}{1000}
× f_{temp}(t)
× f_{IAM}(t)
× f_{spectral}(t)
× f_{lowlight}(t)
× f_{shade}(t)
× f_{mismatch}(t)
× f_{degradation}(t)

それぞれの意味は次の通り。

係数	意味	ペロブスカイトでの注意点
`f_temp`	温度補正	フィルム、壁面、金属屋根、封止で変わる
`f_IAM`	入射角補正	壁面・曲面で効きやすい
`f_spectral`	スペクトル補正	バンドギャップ・タンデムで重要
`f_lowlight`	低照度補正	壁面・曇天・都市影で重要
`f_shade`	影補正	建物密集地で重要
`f_mismatch`	ミスマッチ損失	部分影・面内ばらつきで重要
`f_degradation`	劣化補正	最重要。固定値を安易に置かない

これを年間発電量にすると、

E_{dc,year} = \sum_{t=1}^{N} P_{dc}(t) × Δt

AC発電量は、

E_{ac,year} =
\sum_{t=1}^{N}
P_{dc}(t)
× η_{inv}(t)
× f_{clip}(t)
× f_{availability}(t)
× Δt

となる。

10. 劣化モデル：ペロブスカイトで最も危険な仮定

ペロブスカイトの実務シミュレーションで、最も危険なのは劣化率を雑に置くことだ。

従来の結晶シリコンPVでは、年率劣化を一定値で置く単純モデルがよく使われる。

P(t) = P_0 × (1 - d)^t

しかしペロブスカイトでは、熱、湿気、酸素、紫外線、電界、イオン移動、封止劣化、界面劣化などの支配因子が複雑に絡む。MRS Bulletinのレビューでも、ペロブスカイト太陽電池の安定性、標準化、ISOSプロトコル、IEC試験、商用段階への移行に関する課題が整理されている。(Springer)

加速劣化試験では、しばしばArrhenius型の考え方を使う。

k(T) = k_0 × exp\left(-\frac{E_a}{k_B T}\right)

光強度依存も加えるなら、

k(T,I) = k_0 × I^m × exp\left(-\frac{E_a}{k_B T}\right)

ただし、MRS Bulletinのレビューは、こうした加速モデルが有効なのは、支配的な劣化メカニズムやストレス因子が分かっている場合に限られると指摘している。(Springer)

実務では、次のように分けるのが安全である。

段階	劣化率の扱い
初期概算	複数シナリオで感度分析
実証前	メーカー保証値、認証試験、屋外暴露試験値を確認
実証中	実測発電量から月次・季節別に逆推定
商用判断	P50/P90、交換費、保証、停止リスクを財務モデルに入れる

1つの劣化率で「20年の発電量」を断定するのは、まだ危ない。
ペロブスカイトでは、劣化率は単なる係数ではなく、事業性そのものを左右する設計変数である。

11. LCOE、NPV、回収年数の計算

発電コスト、つまりLCOEは次のように表せる。

LCOE =
\frac{
CAPEX + \sum_{t=1}^{N}\frac{OPEX_t + Replacement_t - Residual_t}{(1+r)^t}
}{
\sum_{t=1}^{N}\frac{E_t}{(1+r)^t}
}

記号	意味
`CAPEX`	初期投資
`OPEX`	運用保守費
`Replacement`	交換費
`Residual`	残存価値
`E_t`	年間発電量
`r`	割引率
`N`	評価年数

ペロブスカイトでLCOEを計算するときは、CAPEX単価だけでなく、次を分ける必要がある。

モジュール費
施工費
下地処理費
接着・固定・防水費
封止・端部処理費
PCS・配線費
足場・高所作業費
保守点検費
交換費
撤去・リサイクル費
保証・保険コスト

ここでまた反直感な論点がある。

ペロブスカイトは、モジュール単価が安いから勝つとは限らない。
むしろ初期市場では、従来PVを設置できない場所で発電価値を作れるから勝つ可能性が高い。

例えば、耐荷重が足りず結晶シリコンPVを置けない屋根に、軽量フィルム型PVを置けるなら、比較対象は「安いシリコンPV」ではない。
比較対象は「発電ゼロ」である。

この場合、LCOEの評価軸は変わる。

従来PVが置ける場所:
  シリコンPV vs ペロブスカイトPV のコスト・発電量比較

従来PVが置けない場所:
  発電ゼロ vs ペロブスカイトPV の価値比較

この比較軸を間違えると、ペロブスカイトの価値を過小評価する。

12. 自家消費・蓄電池・PPAとの接続

需要家向けの導入判断では、単に発電量を出すだけでは足りない。
発電が「いつ」出て、需要が「いつ」あるかを見る必要がある。

自家消費率は、

SCR =
\frac{\sum_t min(PV_t, Load_t)}
{\sum_t PV_t}

自給率は、

SSR =
\frac{\sum_t min(PV_t, Load_t)}
{\sum_t Load_t}

蓄電池を入れると、

PV_t + Discharge_t - Charge_t - Export_t = SelfConsumption_t

という時系列制約になる。

PPAでは、さらに次のような収益・費用を分ける。

需要家側:
  電気料金削減額
  PPA単価支払い
  契約期間
  再エネ価値
  BCP価値

PPA事業者側:
  CAPEX
  OPEX
  発電量
  PPA売上
  余剰売電
  保守リスク
  劣化・交換リスク

ペロブスカイトでは、壁面・軽量屋根・曲面など、発電プロファイルが通常の南面屋根PVと異なる可能性がある。
そのため、年間発電量だけでなく、時間帯別発電量 × 需要カーブで見る必要がある。

13. 係数・パラメーター早見表

以下は、ペロブスカイト太陽光シミュレーションで最低限管理すべきパラメーターである。

区分	パラメーター	単位	代表的な扱い	注意点
日射	`GHI`	W/m²	気象データから取得	水平面値のまま使わない
日射	`DNI`, `DHI`	W/m²	分離モデルまたはデータ取得	Erbs/Perez等の選択で差が出る
面	方位角	deg	設計値	壁面・曲面で重要
面	傾斜角	deg	設計値	垂直面は90度
光学	IAM係数	-	ASHRAE/Martin-Ruiz等	ペロブスカイト実測が望ましい
温度	`a,b,ΔT`	-	Sandiaモデル等	フィルム・屋根一体で要検証
電気	`P_STC`	kW	モジュール仕様	面積効率と混同しない
電気	温度係数 `γ`	1/℃	データシート	シリコン値の流用禁止
電気	I-Vパラメーター	各種	実測またはフィッティング	ヒステリシスに注意
劣化	年率劣化 `d`	%/年	感度分析	固定値断定は危険
経済性	CAPEX	円	個別見積	初期市場は相場未形成
経済性	OPEX	円/年	保守条件	交換・点検費を忘れない
経済性	LCOE	円/kWh	割引現在価値で算出	発電量劣化を反映

14. 相場数値・政策目標・実装初期値

ペロブスカイトの「相場」を語るときは注意が必要である。
2026年5月時点では、結晶シリコンPVのように全国の施工相場が厚く形成されている段階ではない。したがって、公開情報から使える数値は、主に政策目標、量産計画、初期モジュール仕様、補助制度である。

論点	数値	意味
日本の導入目標	2040年約20GW	次世代型太陽電池戦略上の目標
発電コスト目標	2025年 20円/kWh、2030年 14円/kWh、2040年 10〜14円/kWh以下	政策上のコスト目標
NEDO単接合目標	2030年度 14円/kWh以下	GI基金の技術目標
NEDOタンデム目標	2030年度 12円/kWh、変換効率30%	タンデム型の技術目標
積水化学100MWライン	2027年度予定	量産初期の重要マイルストーン
初期モジュール仕様例	110W、約4.5kg/m²	特定公開仕様からの参考値

政策目標はMETIおよびNEDOの公開資料に基づく。
積水化学の量産計画および初期モジュール仕様例は、同社および関連会議資料に基づく。(積水化学株式会社)

重要なのは、これらを「施工単価」として読まないことだ。
例えば、100MWラインへの投資額を、設置単価に単純変換してはいけない。それは工場投資であり、顧客が支払うシステム価格ではない。

実務モデルでは、当面は次のように感度分析するのがよい。

CAPEX:
  Low / Base / High の3ケース

発電量:
  P50 / P90 の2ケース

劣化:
  保守的 / 中位 / 楽観 の3ケース

工期:
  通常 / 建物制約あり / 高所・特殊施工あり

価値:
  自家消費価値 / PPA価値 / BCP価値 / PR価値 を分ける

15. シミュレーションソフト・API・データセット

15.1 ソフトウェア・API比較

ツール	主用途	ペロブスカイトでの使い方
pvlib	PythonによるPV計算	カスタムモジュール・IAM・温度・単一ダイオード実装
NREL SAM / PySAM	詳細な性能・財務モデル	PV＋蓄電池＋経済性評価
PVWatts API	簡易PV発電量API	初期概算、比較用
PVGIS API	欧州中心だがグローバルな日射・PV性能	API経由の概算、地点比較
PVsyst	商用PV設計・影・損失分析	高精度設計、銀行向け資料
Solar Pro	3D CAD影・I-V・国内気象データ活用	日本案件、詳細影解析
HOMER	マイクログリッド・蓄電池・経済性	離島、BCP、ハイブリッド電源
OpenDSS	配電系統解析	大量導入、電圧、逆潮流、系統影響

PVGISは、極域を除く地点で太陽放射とPVシステム性能を評価でき、APIも提供している。(Joint Research Centre)
NSRDBはGOESやHimawari系の衛星データを含み、アジア向けにはHimawariデータも提供されている。(デベロッパーネットワーク)
PVWattsは、世界各地の系統連系PVの発電量を簡易推定するためのNRELツールで、API版も提供されている。(PVWatts Calculator)
PySAMは、SAMの性能モデルと財務モデルをPythonから扱うためのライブラリである。(PySAM)

PVsyst、Solar Pro、HOMER、OpenDSSは、それぞれ詳細PV設計、3D影解析、マイクログリッド経済性、配電系統解析に使われる代表的なツールである。(PVSyst)

15.2 日本向けデータセット

データセット	用途	注意点
NEDO METPV-20	時別日射量	国内PV発電量計算の基礎
NEDO MONSOLA-20	月別日射量	概算・比較向き
気象庁データ	気温・風速・日照等	温度モデルに必要
JWA等の高解像度気象データ	商用高精度計算	ライセンス確認が必要
NSRDB / Himawari	衛星ベース日射	海外・広域比較向き
PVGIS	APIによる地点別概算	欧州中心だが広域対応

JPPAガイドラインでは、NEDOのMETPV-20が2010〜2018年の時別データを含む1kmメッシュベースの日射量データとして整理されている。

16. Pythonで最小限のペロブスカイト発電量モデルを書く

以下は、ペロブスカイト専用モデルではなく、ペロブスカイトをカスタムPVとして扱うための最小実装である。

import numpy as np
import pandas as pd

def sandia_cell_temperature(poa_wm2, temp_air_c, wind_ms,
                            a=-3.58, b=-0.113, delta_t=3.0):
    """
    Sandia系の簡易セル温度モデル。
    デフォルト値は thin film / steel / open rack 系の例。
    実案件では設置方式・メーカー実測で更新する。
    """
    poa_wm2 = np.asarray(poa_wm2)
    temp_module = poa_wm2 * np.exp(a + b * wind_ms) + temp_air_c
    temp_cell = temp_module + (poa_wm2 / 1000.0) * delta_t
    return temp_cell


def ashrae_iam(theta_deg, b0=0.05):
    """
    ASHRAE IAMモデル。
    theta_deg: 入射角[deg]
    b0: 実測またはメーカー値で更新すべき係数
    """
    theta_rad = np.deg2rad(theta_deg)
    cos_theta = np.cos(theta_rad)
    cos_theta = np.maximum(cos_theta, 1e-6)

    iam = 1.0 - b0 * ((1.0 / cos_theta) - 1.0)
    return np.clip(iam, 0.0, 1.0)


def perovskite_dc_power(
    p_stc_kw,
    poa_wm2,
    temp_cell_c,
    theta_deg,
    gamma_pdc=-0.0025,
    b0=0.05,
    spectral_factor=1.0,
    lowlight_factor=1.0,
    shade_factor=1.0,
    mismatch_factor=0.98,
    degradation_factor=1.0
):
    """
    ペロブスカイトPVのDC出力を、簡易損失チェーンで計算する。

    注意:
    gamma_pdc, b0, spectral_factor, lowlight_factor, degradation_factor は
    ペロブスカイト一般の真値ではない。
    実測・データシート・実証結果で更新する。
    """
    poa_ratio = np.maximum(poa_wm2, 0.0) / 1000.0
    temp_factor = 1.0 + gamma_pdc * (temp_cell_c - 25.0)
    iam_factor = ashrae_iam(theta_deg, b0=b0)

    p_dc = (
        p_stc_kw
        * poa_ratio
        * temp_factor
        * iam_factor
        * spectral_factor
        * lowlight_factor
        * shade_factor
        * mismatch_factor
        * degradation_factor
    )

    return np.maximum(p_dc, 0.0)


def ac_energy_kwh(p_dc_kw, inverter_eff=0.96, dt_hours=1.0):
    """
    DC→AC変換後の電力量。
    クリッピングやPCS制御を入れる場合はここを拡張する。
    """
    p_ac_kw = p_dc_kw * inverter_eff
    return np.sum(p_ac_kw * dt_hours)

このコードのポイントは、ペロブスカイトを「特別な魔法の発電源」として扱っていないことだ。
むしろ、通常PVモデルの各係数を、ペロブスカイト固有の実測値で更新できる形にしている。

実務で重要なのは、モデルを複雑にすることではない。

どの係数が仮置きか
どの係数がメーカー値か
どの係数が実測で更新済みか
どの係数が事業性を左右するか

を明示することである。

17. ペロブスカイト向けシミュレーションAPIを設計するなら

再エネSaaSや業務システムに組み込むなら、ペロブスカイトPVは単なる「PV種別」ではなく、設置面と劣化モデルを持つ資産として扱うべきである。

17.1 入力JSONの例

{
  "site": {
    "latitude": 35.6812,
    "longitude": 139.7671,
    "timezone": "Asia/Tokyo"
  },
  "surface": {
    "type": "wall",
    "azimuth_deg": 180,
    "tilt_deg": 90,
    "area_m2": 500,
    "sky_view_factor": 0.65,
    "albedo": 0.2
  },
  "module": {
    "technology": "perovskite_film",
    "p_stc_kw": 35.0,
    "module_efficiency": 0.08,
    "temperature_coefficient": -0.0025,
    "iam_model": "ashrae",
    "iam_b0": 0.05
  },
  "losses": {
    "soiling": 0.98,
    "mismatch": 0.98,
    "availability": 0.99
  },
  "degradation": {
    "model": "scenario",
    "annual_rate_cases": {
      "conservative": 0.03,
      "base": 0.015,
      "optimistic": 0.005
    }
  },
  "inverter": {
    "efficiency": 0.96,
    "dc_ac_ratio": 1.1
  },
  "load_profile": {
    "type": "hourly_csv",
    "unit": "kWh"
  },
  "economic": {
    "capex_yen": 50000000,
    "opex_yen_per_year": 500000,
    "discount_rate": 0.04,
    "analysis_years": 20,
    "electricity_price_yen_per_kwh": 25
  }
}

ここで重要なのは、surfaceを明示していることだ。

ペロブスカイトは、屋根、壁、曲面、低耐荷重屋根、建材一体、仮設、車庫、倉庫外壁など、設置面の多様性が価値になる。
したがってAPI上も、単にcapacity_kwだけでなく、surfaceを第一級オブジェクトにするべきである。

17.2 出力JSONの例

{
  "summary": {
    "annual_ac_kwh_p50": 42000,
    "annual_ac_kwh_p90": 36000,
    "self_consumption_rate": 0.82,
    "self_sufficiency_rate": 0.18,
    "lcoe_yen_per_kwh": 18.5,
    "simple_payback_years": 12.4
  },
  "diagnostics": {
    "dominant_loss": "high_incidence_angle",
    "critical_assumption": "degradation_rate",
    "needs_field_validation": true
  },
  "loss_breakdown": {
    "iam_loss_pct": 12.0,
    "temperature_loss_pct": 4.0,
    "shade_loss_pct": 8.0,
    "mismatch_loss_pct": 2.0,
    "inverter_loss_pct": 4.0
  },
  "scenario_comparison": {
    "conservative": {
      "lcoe_yen_per_kwh": 25.0,
      "payback_years": 18.0
    },
    "base": {
      "lcoe_yen_per_kwh": 18.5,
      "payback_years": 12.4
    },
    "optimistic": {
      "lcoe_yen_per_kwh": 14.2,
      "payback_years": 9.8
    }
  }
}

出力で最も大事なのは、単一の発電量ではなく、診断を返すことだ。

この案件は日射不足なのか
入射角損失なのか
影なのか
温度なのか
劣化仮定なのか
施工費なのか

これが分からない発電量APIは、意思決定に使いにくい。

18. システムエンジニアリング上の論点

ペロブスカイト太陽光発電システムは、セル材料の話だけでは終わらない。
商用システムとしては、次の設計が必要になる。

領域	論点
構造	風荷重、接着、剥離、下地、曲面追従
防水	端部処理、貫通部、雨水侵入、屋根保証
電気	配線、絶縁、接地、短絡、直流アーク
熱	放熱、熱だまり、金属屋根温度
劣化	湿気、酸素、紫外線、熱、電界
安全	鉛含有、破損時回収、リサイクル
保守	点検方法、交換単位、遠隔監視
金融	保証、保険、性能保証、バンカビリティ
制度	JET評価、補助対象、建築基準、防火

環境・健康リスクについては、多くのペロブスカイト組成が鉛を含むため、破損・水濡れ・廃棄・リサイクル時の管理が論点になる。レビュー論文でも、鉛含有ペロブスカイト太陽電池の環境・健康リスクが重要課題として扱われている。(PMC)

国内の社会実装支援では、2025年度の環境省事業において、ペロブスカイト太陽電池の社会実装モデル創出支援が公募され、補助対象はJETによる性能評価等を満たすフィルム型ペロブスカイト太陽電池とされている。(環境省)

つまり、実務導入では「軽くて貼れる」だけでは足りない。

貼れる
発電する
剥がれない
漏水しない
安全である
劣化を説明できる
保証できる
金融・保険・補助制度に乗る

ここまで揃って、初めてシステムになる。

19. よくある誤解

誤解1：ペロブスカイトはシリコンPVをすぐ置き換える

短期的には、すべてを置き換えるというより、シリコンPVが苦手な場所を補完すると見る方が現実的である。

特に、低耐荷重屋根、壁面、曲面、建材一体、都市部の未利用面が重要になる。

誤解2：変換効率が高ければ発電量も必ず多い

発電量は、効率だけでは決まらない。

発電量 = 効率 × 日射 × 設置面 × 温度 × 影 × 入射角 × 劣化 × 稼働率

壁面や曲面では、入射角・影・反射・低照度が大きく効く。
研究セル効率だけで案件の発電量を語るのは危険である。

誤解3：薄くて軽いから施工費も必ず安い

軽いことは大きな利点だが、防水、接着、端部処理、屋根保証、足場、安全対策、交換性が追加論点になる。

軽量化は施工費を下げる可能性がある。
しかし、初期市場では特殊施工費や検証コストが上乗せされる可能性もある。

誤解4：年間発電量だけ出せば導入判断できる

需要家案件では、年間発電量だけでは不十分である。
必要なのは、時間帯別発電量、需要カーブ、自家消費率、余剰、蓄電池、電気料金単価、PPA単価、補助金、劣化、保証をつないだ評価である。

20. 実務での最小入力セット

ペロブスカイトPV案件を概算する場合、最低限必要な入力は次の通りである。

入力	必須度	理由
緯度・経度	必須	太陽位置・日射計算
設置面の方位・傾斜	必須	POA日射に直結
設置面積	必須	容量上限の推定
モジュール定格・効率	必須	P_STC計算
設置方式	必須	温度モデルに影響
周辺遮蔽	高	壁面・都市部で重要
IAM係数	高	高入射角で重要
温度係数	高	発電量に影響
劣化シナリオ	高	LCOEに直結
需要カーブ	需要家案件では必須	自家消費率計算
電気料金単価	需要家案件では必須	経済効果計算
CAPEX/OPEX	必須	LCOE/NPV計算

最小努力で始めるなら、次の3つから集める。

1. 設置面:
   面積、方位、傾斜、影

2. モジュール:
   定格、重量、効率、温度係数、保証

3. 需要:
   30分値または月別使用量、料金単価

この3つがない状態で回収年数を語ると、ほぼ占いになる。占いはQiitaではなく別媒体でやるべきだ。

21. ペロブスカイト発電量シミュレーションの設計チェックリスト

□ 研究セル効率と製品モジュール効率を分けたか
□ 設置面の方位・傾斜・天空率を入れたか
□ GHIをそのまま使わずPOAに変換したか
□ 直達・散乱・反射を分けたか
□ 壁面・曲面の入射角損失を見たか
□ 近傍建物の影を評価したか
□ 温度モデルを設置方式に合わせたか
□ 温度係数をデータシート値で確認したか
□ 劣化率を単一固定値で断定していないか
□ P50/P90や保守的ケースを出したか
□ 自家消費率と自給率を分けたか
□ 蓄電池を入れる場合、充放電制約を入れたか
□ LCOEに交換費・保守費・撤去費を入れたか
□ 鉛・廃棄・リサイクル・安全対策を確認したか
□ JET評価、補助対象、保証条件を確認したか

22. エネがえる文脈での実装アイデア

ペロブスカイト太陽光は、再エネシミュレーションSaaS/APIにとって非常に面白い対象である。理由は、単に新技術だからではない。
試算の難易度が上がるほど、標準化された計算基盤の価値が上がるからである。

例えば、エネがえるのような経済効果シミュレーション基盤であれば、次のような拡張が考えられる。

ペロブスカイトPV入力:
  - 設置面種別: 屋根 / 壁面 / 曲面 / 金属屋根 / 建材一体
  - モジュール種別: フィルム型 / ガラス型 / タンデム型
  - 面密度: kg/m²
  - IAMモデル
  - 温度モデル
  - 劣化シナリオ

発電量計算:
  - NEDO METPV-20
  - 高解像度日射データ
  - POA変換
  - 壁面補正
  - 影補正

経済性:
  - 自家消費
  - PPA
  - 補助金
  - 蓄電池併設
  - 劣化・交換費
  - LCOE / NPV / 回収年数

出力:
  - 発電量
  - 損失内訳
  - 支配要因
  - 追加確認事項
  - P50/P90
  - 稟議用レポート

特に価値が出るのは、次の3つである。

22.1 設置できない屋根を「発電候補」に変える

低耐荷重屋根では、従来PVが構造上難しい場合がある。
ペロブスカイトはここに入れる可能性がある。

この場合、シミュレーションは単なる発電量計算ではなく、

構造制約 → 設置可否 → 発電量 → 自家消費 → 経済性

をつなぐ意思決定支援になる。

22.2 壁面PVの発電量を説明可能にする

壁面PVでは、なぜ発電量が少ないのか、どの面が効くのか、午前・午後・季節でどう変わるのかを説明する必要がある。

単一の年間kWhでは足りない。
損失内訳と時間帯別の可視化が必要になる。

22.3 補助金・PPA・蓄電池を含めた比較にする

ペロブスカイトは初期市場で補助制度や実証事業と絡みやすい。
そのため、設備単価だけでなく、補助金、PPA契約、自家消費単価、余剰売電、蓄電池併設を含めた比較が必要になる。

ここはExcel手作業では壊れやすい。
API化・BPO化・標準レポート化の余地が大きい。

23. 実装ロードマップ

Phase 1：概算モデル

目的:
  技術検討・営業初期提案

入力:
  住所、面積、方位、傾斜、概算容量

モデル:
  METPV-20 / PVGIS / PVWatts相当
  簡易POA
  簡易温度
  固定損失
  劣化3シナリオ

出力:
  年間発電量
  概算LCOE
  自家消費率
  注意点

Phase 2：設計モデル

目的:
  提案・稟議・PPA検討

入力:
  3D影、需要カーブ、料金単価、モジュール仕様、施工費

モデル:
  IAM
  近傍影
  温度モデル
  インバータ
  劣化シナリオ
  蓄電池制御

出力:
  月別・時間帯別発電量
  損失内訳
  P50/P90
  NPV/LCOE
  比較表

Phase 3：実証・運用モデル

目的:
  実測検証・性能保証・金融・保険

入力:
  実測発電量、日射、温度、風速、異常ログ

モデル:
  デジタルツイン
  パラメーター逆推定
  劣化率推定
  異常検知

出力:
  実測対予測
  劣化推定
  保守アラート
  保証判定
  レポート

24. FAQ

Q1. ペロブスカイト太陽電池はもう実用化されていますか？

研究段階だけではなく、実証・商用初期・量産準備へ進んでいる。国内では積水化学がフィルム型ペロブスカイト太陽電池の事業開始と2027年度100MW級ラインの方針を公表している。ただし、結晶シリコンPVのような広範な施工相場が確立した段階とは分けて見るべきである。(積水化学株式会社)

Q2. ペロブスカイトはシリコンより高効率ですか？

単接合同士で単純比較するより、ペロブスカイト/シリコンタンデムで30%台の高効率を狙える点が重要である。NRELの研究セル効率チャートでは、タンデム領域で高い効率記録が示されている。(NREL)

Q3. 発電量計算は従来PVと同じでよいですか？

基本パイプラインは同じだが、壁面・曲面・低耐荷重屋根などで使う場合、入射角損失、影、低照度、温度、劣化の扱いが重要になる。特に壁面PVでは、JPPAガイドラインでも高入射角、近傍遮蔽、天空率、MPPTなどへの注意が整理されている。

Q4. ペロブスカイトの劣化率は何%で置けばよいですか？

現時点では、案件や製品仕様に依存するため、固定値で断定しない方がよい。メーカー保証、認証試験、屋外実証、実測データを確認し、保守的・中位・楽観の複数シナリオで評価するのが安全である。

Q5. ペロブスカイトはどこに向いていますか？

従来PVが置きにくい場所に向く。具体的には、低耐荷重屋根、壁面、曲面、建材一体、都市部の未利用面などである。すでにシリコンPVを安く大量に置ける広い屋根では、短期的にはシリコンPVの方が有利なケースも多い。

Q6. 発電コストの目標はどの程度ですか？

日本の政策資料では、次世代型太陽電池について2025年20円/kWh、2030年14円/kWh、2040年10〜14円/kWh以下といった目標が示されている。NEDOのGI基金でも、2030年度に単接合14円/kWh以下、タンデム型12円/kWhなどの目標が掲げられている。

Q7. ペロブスカイトPVのシミュレーションに最適なソフトはどれですか？

目的による。Pythonで柔軟に組むならpvlib、性能・財務モデルならSAM/PySAM、API概算ならPVGISやPVWatts、詳細設計ならPVsystやSolar Pro、系統影響ならOpenDSSが候補になる。ペロブスカイト固有のパラメーターは、既存ツールの標準DBにない場合が多いため、カスタム入力が前提になる。

25. まとめ：ペロブスカイトPVは「材料技術」ではなく「設置面を再定義するシステム技術」である

ペロブスカイト太陽光発電を理解するうえで、最も重要な認識はこれである。

ペロブスカイトの価値は、
シリコンPVを安く置き換えることだけではない。

これまで発電面として扱えなかった屋根・壁・曲面・建物外皮を、
発電資産として計算可能にすることにある。

そのため、発電量シミュレーションも変わる。

従来のように、

容量 × 年間日射 × PR

で終わらせるには無理がある。

これから必要になるのは、

設置面
日射
影
入射角
温度
劣化
需要
料金
蓄電池
補助金
保証

をつないだ、説明可能なシステムモデルである。

ペロブスカイトは、太陽電池の材料を変えるだけではない。
発電量シミュレーション、PPA設計、建築外皮、都市エネルギー、BPO、API、保証、金融の接続点を増やす。

つまり、エンジニアにとってはかなり面白い。
ただし、雑に扱うと一瞬で壊れる。
そこがいい。

参考文献・一次情報

DOE Solar Energy Technologies Office「Perovskite Solar Cells」(The Department of Energy's Energy.gov)
NREL Best Research-Cell Efficiencies Chart(NREL)
経済産業省次世代型太陽電池戦略関連資料
NEDO グリーンイノベーション基金次世代型太陽電池関連資料(NEDO グリーンイノベーション基金)
積水化学フィルム型ペロブスカイト太陽電池関連公表資料(積水化学株式会社)
JPPA「壁面設置太陽光発電システム発電量推定手法ガイドライン」
NREL SAM Photovoltaic Model Technical Reference(NREL)
pvlib documentation(Pvlib Python)
PVGIS overview and API documentation(Joint Research Centre)
NSRDB Developer API(デベロッパーネットワーク)
NREL PVWatts / PySAM documentation(PVWatts Calculator)
MRS Bulletin review on perovskite commercialization and stability(Springer)