(57)【要約】【課題】太陽光発電システムの全体システムとしての効率はシステム出力係数により表されることが多いが、その損失分のうちどの程度が太陽電池温度上昇による損失であるかは明らかでなかった。そのためシステムの諸損失の解明ができずシステムの効率向上に結びつかなかった。【解決手段】「システム出力係数」に「温度上昇による損失」を加えた効率を算出する２つの方法（「月別温度係数」による方法、「温度補正後システム出力係数」による方法）を開発し、両方法で算出した結果がよく似ていることが明らかとなった。またいずれの方法も、任意の日射強度・太陽電池温度のＩ−Ｖカーブの作成方法として「理論式によるＩ−Ｖカーブ作成法」すれば効果的であることがわかった、【選択図】        図１
2005051048.tif 000002
【特許請求の範囲】
【請求項１】
対象太陽光発電システムの月ごとのシステム出力係数（％）（＝（月間実測システム発電量（ＫＷｈ））/ (月間受光面日射量（ＫＷｈ/ m2）)/太陽電池定格容量（ＫＷ/ＫＷ/ m2）＊１００）を算出し、
そして、太陽電池の温度上昇による損失（％）（＝（１００−月別温度係数））を算出するために、月別温度係数（％）（＝月間計算システム発電量（ＫＷｈ）/月間受光面日射量（ＫＷｈ/ m2）/太陽電池定格容量（ＫＷ/ＫＷ/ m2）＊１００）計算のために必要な月間計算システム発電量の算出にさらに必要な、任意の日射強度・太陽電池温度の太陽電池Ｉ−Ｖカーブを作成する方法として、（1）標準太陽電池のＩ−Ｖカーブを変換する方法（「実用的Ｉ−Ｖカーブ作成法」）又は、（2）太陽電池基本式を解く方法（「理論式によるＩ−Ｖカーブ作成法」）を使い、温度上昇による損失（％）を算出し、システム出力係数に該温度上昇による損失を加えた値により、
太陽光発電システムの発電量を評価することを特徴とする評価方法
【請求項２】
対象太陽光発電システムの月ごとのシステム出力係数（％）（＝（月間実測システム発電量（ＫＷｈ））/ (月間受光面日射量（ＫＷｈ/ m2）)/太陽電池定格容量（ＫＷ/ＫＷ/ m2）＊１００）を算出し、
そして、太陽電池の温度上昇による損失（％）（＝（１００−月別温度係数））を該システム出力係数に加えた値により太陽光発電システムの発電量を評価することを特徴とする方法において、
月別温度係数（％）（＝（月間計算システム発電量（ＫＷｈ）/月間受光面日射量（ＫＷｈ/ m2）/太陽電池定格容量（ＫＷ/ＫＷ/ m2）＊１００）の算出に必要な月間計算システム発電量の算出過程で必要な、任意の日射強度・太陽電池温度の発電量算出に必要なＩ―Ｖカーブ作成方法において、
｛01｝電圧Ｖ、電流Ｉ、日射強度１kW/m2 での光起電流IL、飽和電流温度係数Co、接合定数n、並列抵抗Rsh、直列抵抗Rs、太陽電池モジュール温度T （絶対温度）を含んだ
関数：Func(V,I,IL,Co,n,Rsh,Rs,T) = IL - CoT3 exp(-qEg/nk0T)*(exp( q*(V+Rs*I)/(n*k0*T) )-1) - (V+Rs*I)/ Rsh - I  を作成し、つぎに、
｛02｝該関数Func(V,I,IL,Co,n,Rsh,Rs,T)を変数Vで微分した
関数：Div(V,I,IL,Co,n,Rsh,Rs,T)を作成し、
｛03｝太陽電池の基準状態（モジュール温度Ta(298K(ta=25℃))、日射強度Ea（１kW/m2 ））での仕様値、短絡電流Isca、最適電流Iopa 、最適電圧Vopa 、開放電圧Vocaの点P1(0,Isca),P2(Vopa,Iopa),P3(Voca,0)を選択し、
｛04｝前記関数Func(V,I,IL,Co,n,Rsh,Rs,T)= 0のTに基準状態の温度Ta(298K),直列抵抗Rsに基準温度での仕様値Rsa,および前記P1,P2,P3の点の値を代入し、IL,Co,n,Rshを未知数とする
関係式:Func(0,Isca,IL,Co,n,Rsh,Rsa,Ta)= 0,
関係式:Func(Voca,0,IL,Co,n,Rsh,Rsa,Ta)= 0,
関係式:Func(Vopa,Iopa,IL,Co,n,Rsh,Rsa,Ta)= 0を作成し、
｛05｝前記関数Div(V,I,IL,Co,n,Rsh,Rs,T)= 0 に、基準状態の温度Ta(298K)の直列抵抗Rsに値Rsa および前記点P2の値(Vopa,Iopa) を代入して、IL,Co,n,Rshを未知数とする、
関係式:Div(Vopa,Iopa,IL,Co,n,Rsh,Rsa,Ta)= 0を作成し、つぎに、
｛06｝前記４つの関係式:Func(0,Isca,IL,Co,n,Rsh,Rsa,Ta)= 0,Func(Voca,0,IL,Co,n,Rsh,Rsa,Ta)= 0,Func(Vopa,Iopa,IL,Co,n,Rsh,Rsa,Ta)= 0,Div(Vopa,Iopa,IL,Co,n,Rsh,Rsa,Ta)= 0 を満たす解A(ILa,Coa,na,Rsha)を、非線形解法のプログラムによって算出し、つぎに、
｛07｝太陽電池の基準状態と同様に、第２の太陽電池モジュール温度Tb(313K(tb=40℃))・日射強度Eb（１kW/m2 ）)での仕様値である、短絡電流Iscb、最適電流Iopb 、最適電圧Vopb 、開放電圧Vocbの点P1(0,Iscb),P2(Vopb,Iopb),P3(Vocb,0)を選択し、
｛08｝前記関数Func(V,I,IL,Co,n,Rsh,Rs,T)= 0のTに第２の温度Tb(313K),直列抵抗Rsに該温度での仕様値Rsb,および前記{07}のP1,P2,P3の点の値を代入し、IL,Co,n,Rshを未知数とする
関係式:Func(0,Iscb,IL,Co,n,Rsh,Rsb,Tb)= 0,
関係式:Func(Vocb,0,IL,Co,n,Rsh,Rsb,Tb)= 0,
関係式:Func(Vopb,Iopb,IL,Co,n,Rsh,Rsb,Tb)= 0を作成し、
｛09｝前記関数Div(V,I,IL,Co,n,Rsh,Rs,T)= 0 に、第２の温度Tb(313K)の直列抵抗Rsに第２の温度での値Rsb および前記{07}の点P2の値(Vopb,Iopb) を代入して、IL,Co,n,Rshを未知数とする、
関係式:Div(Vopb,Iopb,IL,Co,n,Rsh,Rsb,Tb)= 0を作成し、つぎに、
｛10｝前記４つの関係式:Func(0,Iscb,IL,Co,n,Rsh,Rsb,Tb)= 0,Func(Vocb,0,IL,Co,n,Rsh,Rsb,Tb)= 0,Func(Vopb,Iopb,IL,Co,n,Rsh,Rsb,Tb)= 0,Div(Vopb,Iopb,IL,Co,n,Rsh,Rsb,Tb)= 0 を満たす解A(Ilb,Cob,nb,Rshb)を、非線形解法のプログラムによって、算出し、つぎに、
｛11｝太陽電池の基準状態と同様に第３の太陽電池モジュール温度Tc(328K(tc=55℃))・日射強度Ec（１kW/m2 ）)での仕様値である、短絡電流Iscc、最適電流Iopc 、最適電圧Vopc 、開放電圧Voccの点P1(0,Iscc),P2(Vopc,Iopc),P3(Vocc,0)を選択し、
｛12｝前記関数Func(V,I,IL,Co,n,Rsh,Rs,T)= 0のTに第３の温度Tc(328K),直列抵抗Rsに該温度での仕様値Rsc,および前記{11}のP1,P2,P3の点の値を代入し、IL,Co,n,Rshを未知数とする
関係式:Func(0,Iscc,IL,Co,n,Rsh,Rsc,Tc)= 0,
関係式:Func(Vocc,0,IL,Co,n,Rsh,Rsc,Tc)= 0,
関係式:Func(Vopc,Iopc,IL,Co,n,Rsh,Rsc,Tc)= 0を作成し、
｛13｝前記関数Div(V,I,IL,Co,n,Rsh,Rs,T)= 0 に、第３の温度Tc(328K)の直列抵抗Rsに第３の温度での値Rsc および前記{11}の点P2の値(Vopc,Iopc) を代入して、IL,Co,n,Rshを未知数とする、
関係式:Div(Vopc,Iopc,IL,Co,n,Rsh,Rsc,Tc)= 0を作成し、つぎに、
｛14｝前記４つの関係式:Func(0,Iscc,IL,Co,n,Rsh,Rsc,Tc)= 0,Func(Vocc,0,IL,Co,n,Rsh,Rsc,Tc)= 0,Func(Vopc,Iopc,IL,Co,n,Rsh,Rsc,Tc)= 0,Div(Vopc,Iopc,IL,Co,n,Rsh,Rsc,Tc)= 0 を満たす解A(Ilc,Coc,nc,Rshc)を、非線形解法のプログラムによって算出し、つぎに、
｛15｝評価すべき太陽電池の実測値日射強度Ej,太陽電池モジュール温度tj(摂氏:絶対温度Tj=tj＋273)およびこの条件での発生電圧Vj−発生電流Ijの各値を取り込み、
｛16｝基準状態での前記温度ta(摂氏25℃:絶対温度Ta(K)=ta(℃)＋273) における前記｛06｝の解A(ILa,Coa,na,Rsha)、前記温度tb(摂氏:Tb=tb＋273) における前記｛10｝の解B(ILb,Cob,nb,Rshb)、前記温度tc(摂氏:Tc=tc＋273) における前記｛14｝の解C(ILc,Coc,nc,Rshc)および、入力値Rsa、Rsb、Rsa、Rsc のそれぞれ(IL,Co,n,Rsh,Rs)に関して３点について曲線補間をして、実測温度tj(摂氏:Tj=tj＋273)での特性値M(ILm,Com,nm,Rshm,Rsm) を算出し、つぎに、
｛17｝ILmを実測された日射強度EjによりIL'm=Ilm*Ej/Ea により補正した後、前記関係式：Func(V,I,IL,Co,n,Rsh,Rs,T)=0にIL'm,Com,nm,Rshm,Rsmを代入して、Func(V,I,IL'm,Com,nm,Rshm,Rsm,Tj)=0 を作成し、電圧(V)-電流(I)の関係(約40〜50点)を非線形解法のプログラムによって求め、電圧(V)-電流(I)の関係またはそれを結んだI-V カーブ,P-V カーブを作成し、該カーブを使用すること
を特徴とする太陽光発電量評価方法。
【請求項３】
対象太陽光発電システムの月ごとのシステム出力係数（％）（＝（月間実測システム発電量（ＫＷｈ））/ (月間受光面日射量（ＫＷｈ/ m2）)/太陽電池定格容量（ＫＷ/ＫＷ/ m2）＊１００）を算出し、
そして、太陽電池の温度上昇による損失（％）（＝（１００−月別温度係数））を算出するために、月別温度係数（％）（＝（月間計算システム発電量（ＫＷｈ）/月間受光面日射量（ＫＷｈ/ m2）/太陽電池定格容量（ＫＷ/ＫＷ/ m2）＊１００）を計算し、システム出力係数に該損失を加えた値により評価することに特徴があり、本計算に必要な月間計算システム発電量の算出過程に必要な、任意の日射強度・太陽電池温度の発電量算出のためのＩ−Ｖカーブ作成方法において、
｛18｝電圧Ｖ、電流Ｉ、日射強度１kW/m2 での光起電流IL、飽和電流温度係数Co、接合定数n、並列抵抗Rsh、直列抵抗Rs、太陽電池モジュール温度T （絶対温度）を含んだ
関数：Func(V,I,IL,Co,n,Rsh,Rs,T) = IL - CoT3 exp(-qEg/nk0T)*(exp( q*(V+Rs*I)/(n*k0*T) )-1) - (V+Rs*I)/ Rsh - I  を作成し、つぎに、
｛19｝該関数Func(V,I,IL,Co,n,Rsh,Rs,T)を変数Vで微分した
関数：Div(V,I,IL,Co,n,Rsh,Rs,T)を作成し、
｛20｝太陽電池の基準状態（モジュール温度Ta(298K(ta=25℃))、日射強度Ea（１kW/m2 ））での仕様値、短絡電流Isca、最適電流Iopa 、最適電圧Vopa 、開放電圧Vocaの点P1(0,Isca),P2(Vopa,Iopa),P3(Voca,0)およびこれら３点とは近接しない任意のP4(V4a,I4a)を選択し、
｛21｝前記関数Func(V,I,IL,Co,n,Rsh,Rs,T)= 0のTに基準状態の温度Ta(298K),前記P1,P2,P3,P4の点の値を代入し、IL,Co,n,Rsh,Rsを未知数とする
関係式:Func(0,Isca,IL,Co,n,Rsh,Rs,Ta)= 0,
関係式:Func(Voca,0,IL,Co,n,Rsh,Rs,Ta)= 0,
関係式:Func(Vopa,Iopa,IL,Co,n,Rsh,Rs,Ta)= 0、
関係式:Func(V4a,I4a,IL,Co,n,Rsh,Rs,Ta)= 0を作成し、
｛22｝前記関数Div(V,I,IL,Co,n,Rsh,Rs,T)= 0 に、基準状態の温度Ta(298K)の前記点P2の値(Vopa,Iopa) を代入して、IL,Co,n,Rsh,Rsを未知数とする、
関係式:Div(Vopa,Iopa,IL,Co,n,Rsh,Rs,Ta)= 0を作成し、つぎに、
｛23｝前記５つの関係式:Func(0,Isca,IL,Co,n,Rsh,Rs,Ta)= 0,Func(Voca,0,IL,Co,n,Rsh,Rs,Ta)= 0,Func(Vopa,Iopa,IL,Co,n,Rsh,Rs,Ta)= 0,Div(Vopa,Iopa,IL,Co,n,Rsh,Rs,Ta)= 0、Func(V4a,I4a,IL,Co,n,Rsh,Rs,Ta)= 0 を満たす解A(ILa,Coa,na,Rsha,Rsa)を、非線形解法のプログラムによって算出し、つぎに、
｛24｝太陽電池の基準状態と同様に、第２の太陽電池モジュール温度Tb(313K(tb=40℃))・日射強度Eb（１kW/m2 ）)での仕様値である、短絡電流Iscb、最適電流Iopb 、最適電圧Vopb 、開放電圧Vocbの点P1(0,Iscb),P2(Vopb,Iopb),P3(Vocb,0), P4(V4b,I4b)を選択し、
｛25｝前記関数Func(V,I,IL,Co,n,Rsh,Rs,T)= 0のTに第２の温度Tb(313K),および前記{07}のP1,P2,P3,P4の点の値を代入し、IL,Co,n,Rshを未知数とする
関係式:Func(0,Iscb,IL,Co,n,Rsh,Rs,Tb)= 0,
関係式:Func(Vocb,0,IL,Co,n,Rsh,Rs,Tb)= 0,
関係式:Func(Vopb,Iopb,IL,Co,n,Rsh,Rs,Tb)= 0、
関係式:Func(V4b,I4b ,IL,Co,n,Rsh,Rs,Tb)= 0を作成し、
｛26｝前記関数Div(V,I,IL,Co,n,Rsh,Rs,T)= 0 に、第２の温度Tb(313K)の前記{24}の点P2の値(Vopb,Iopb) を代入して、IL,Co,n,Rsh,Rsを未知数とする、
関係式:Div(Vopb,Iopb,IL,Co,n,Rsh,Rs,Tb)= 0を作成し、つぎに、
｛27｝前記５つの関係式:Func(0,Iscb,IL,Co,n,Rsh,Rs,Tb)= 0,Func(Vocb,0,IL,Co,n,Rsh,Rs,Tb)= 0,Func(Vopb,Iopb,IL,Co,n,Rsh,Rs,Tb)= 0,Div(Vopb,Iopb,IL,Co,n,Rsh,Rs,Tb)= 0, Func(V4b,I4b ,IL,Co,n,Rsh,Rs,Tb)= 0 を満たす解A(Ilb,Cob,nb,Rshb,Rsb)を、非線形解法のプログラムによって、算出し、つぎに、
｛28｝太陽電池の基準状態と同様に第３の太陽電池モジュール温度Tc(328K(tc=55℃))・日射強度Ec（１kW/m2 ）)での仕様値である、短絡電流Iscc、最適電流Iopc 、最適電圧Vopc 、開放電圧Voccの点P1(0,Iscc),P2(Vopc,Iopc),P3(Vocc,0), P4(V4c,I4c)を選択し、
｛29｝前記関数Func(V,I,IL,Co,n,Rsh,Rs,T)= 0のTに第３の温度Tc(328K),前記{28}のP1,P2,P3,P4の点の値を代入し、IL,Co,n,Rsh,Rsを未知数とする
関係式:Func(0,Iscc,IL,Co,n,Rsh,Rs,Tc)= 0,
関係式:Func(Vocc,0,IL,Co,n,Rsh,Rs,Tc)= 0,
関係式:Func(Vopc,Iopc,IL,Co,n,Rsh,Rs,Tc)= 0,
関係式:Func(V4b,I4b,IL,Co,n,Rsh,Rs,Tc)= 0を作成し、
｛30｝前記関数Div(V,I,IL,Co,n,Rsh,Rs,T)= 0 に、第３の温度Tc(328K)の前記{28}の点P2の値(Vopc,Iopc) を代入して、IL,Co,n,Rsh,Rsを未知数とする、
関係式:Div(Vopc,Iopc,IL,Co,n,Rsh,Rs,Tc)= 0を作成し、つぎに、
｛31｝前記５つの関係式:Func(0,Iscc,IL,Co,n,Rsh,Rs,Tc)= 0,Func(Vocc,0,IL,Co,n,Rsh,Rs,Tc)= 0,Func(Vopc,Iopc,IL,Co,n,Rsh,Rs,Tc)= 0,Div(Vopc,Iopc,IL,Co,n,Rsh,Rs,Tc)= 0, Func(V4b,I4b,IL,Co,n,Rsh,Rs,Tc)= 0 を満たす解A(Ilc,Coc,nc,Rshc,Rsc)を、非線形解法のプログラムによって算出し、つぎに、
｛32｝評価すべき太陽電池の実測値日射強度Ej,太陽電池モジュール温度tj(摂氏:絶対温度Tj=tj＋273)およびこの条件での発生電圧Vj−発生電流Ijの各値を取り込み、
｛33｝基準状態での前記温度ta(摂氏25℃:絶対温度Ta(K)=ta(℃)＋273) における前記｛23｝の解A(ILa,Coa,na,Rsha,Rsa)、前記温度tb(摂氏:Tb=tb＋273) における前記｛27｝の解B(ILb,Cob,nb,Rshb,Rsb)、前記温度tc(摂氏:Tc=tc＋273) における前記｛31｝の解C(ILc,Coc,nc,Rshc,Rsc)のそれぞれ(IL,Co,n,Rsh,Rs)に関して３点について曲線補間をして、実測温度tj(摂氏:Tj=tj＋273)での特性値M(ILm,Com,nm,Rshm,Rsm) を算出し、つぎに、
｛34｝ILmを実測された日射強度EjによりIL'm=Ilm*Ej/Ea により補正した後、前記関係式：Func(V,I,IL,Co,n,Rsh,Rs,T)=0にIL'm,Com,nm,Rshm,Rsmを代入して、Func(V,I,IL'm,Com,nm,Rshm,Rsm,Tj)=0 を作成し、電圧(V)-電流(I)の関係(約40〜50点)を非線形解法のプログラムによって求め、電圧(V)-電流(I)の関係またはそれを結んだI-V カーブ,P-V カーブを作成して、該カーブを使用すること
を特徴とする太陽光発電量評価方法。
【請求項４】
対象太陽光発電システムの月ごとのシステム出力係数（％）（＝（月間実測システム発電量（ＫＷｈ））/ (月間受光面日射量（ＫＷｈ/ m2）)/太陽電池定格容量（ＫＷ/ＫＷ/ m2）＊１００）を算出し、
そして、太陽電池の温度上昇による損失（％）（＝（１００−月別温度係数））を算出するために、月別温度係数（％）（＝（月間計算システム発電量（ＫＷｈ））/月間受光面日射量（ＫＷｈ/ m2）/太陽電池定格容量（ＫＷ/ＫＷ/ m2）＊１００）を計算し、システム出力係数に本損失を加えた値により評価することを特徴としており、月間計算システム発電量算出過程にお任意の日射強度・太陽電池温度の発電量の算出に必要な、Ｉ−Ｖカーブ作成方法において、
｛35｝電圧Ｖ、電流Ｉ、日射強度１kW/m2 での光起電流IL、飽和電流温度係数Co、接合定数n、並列抵抗Rsh、直列抵抗Rs、太陽電池モジュール温度T （絶対温度）を含んだ
関数：Func(V,I,IL,Co,n,Rsh,Rs,T) = IL - CoT3 exp(-qEg/nk0T)*(exp( q*(V+Rs*I)/(n*k0*T) )-1) - (V+Rs*I)/ Rsh - I  を作成し、つぎに、
｛36｝太陽電池の基準状態（モジュール温度Ta(298K(ta=25℃))、日射強度Ea（１kW/m2 ））でのメーカーから与えられるｎ個の点P1(V1a,I1a),P2(V2a,I2a),P3(V3a,I3a),・・・PN(VNa,INa)を、前記の関係式:Func(V,I,IL,Co,n,Rsh,Rs,T)= 0に入れてｎ個の関係式Func(V,I,IL,Co,n,Rsh,Rs,Ta)= 0を作成し、
｛37｝該ｎ個の関係式に最も適したIla,Coa,na,Rsha,Rsaを最小自乗問題解法などの数学的な手法で解いて、つぎに、
｛38｝太陽電池の基準状態と同様に、第２の太陽電池モジュール温度Tb(313K(tb=40℃))・日射強度Eb（１kW/m2 ）)でのメーカーから与えられるｎ個の点P1(V1b,I1b),P2(V2b,I2b),P3(V3b,I3b),・・・PN(VNb,Inb)を、前記の関係式:Func(V,I,IL,Co,n,Rsh,Rs,T)= 0に入れてｎ個の関係式Func(V,I,IL,Co,n,Rsh,Rs,Tb)= 0を作成し、
｛39｝該ｎ個の関係式に最も適したIlb,Cob,nb,Rshb,Rsbを最小自乗問題解法などの数学的な手法で解いて、つぎに、
｛40｝太陽電池の基準状態と同様に、第３の太陽電池モジュール温度Tb(328K(tb=55℃))・日射強度Eb（１kW/m2 ）)でのメーカーから与えられるｎ個の点P1(V1c,I1c),P2(V2c,I2c),P3(V3c,I3c),・・・PN(VNc,Inc)を、前記の関係式:Func(V,I,IL,Co,n,Rsh,Rs,T)= 0に入れてｎ個の関係式Func(V,I,IL,Co,n,Rsh,Rs,Tc)= 0を作成し、
｛41｝該ｎ個の関係式に最も適したIlc,Coc,nc,Rshc,Rscを最小自乗問題解法などの数学的な手法で解いて、つぎに、
｛42｝評価すべき太陽電池の実測値日射強度Ej,太陽電池モジュール温度tj(摂氏:絶対温度Tj=tj＋273)およびこの条件での発生電圧Vj−発生電流Ijの各値を取り込み、
｛43｝基準状態での前記温度ta(摂氏25℃:絶対温度Ta(K)=ta(℃)＋273) における前記｛37｝の解A(ILa,Coa,na,Rsha,Rsa)、前記温度tb(摂氏:Tb=tb＋273) における前記｛39｝の解B(ILb,Cob,nb,Rshb,Rsb)、前記温度tc(摂氏:Tc=tc＋273) における前記｛41｝の解C(ILc,Coc,nc,Rshc,Rsc) のそれぞれ(IL,Co,n,Rsh,Rs)に関して３点について曲線補間をして、実測温度tj(摂氏:Tj=tj＋273)での特性値M(ILm,Com,nm,Rshm,Rsm) を算出し、つぎに、
｛44｝ILmを実測された日射強度EjによりIL'm=Ilm*Ej/Ea により補正した後、前記関係式：Func(V,I,IL,Co,n,Rsh,Rs,T)=0にIL'm,Com,nm,Rshm,Rsmを代入して、Func(V,I,IL'm,Com,nm,Rshm,Rsm,Tj)=0 を作成し、電圧(V)-電流(I)の関係(約40〜50点)を非線形解法のプログラムによって求め、電圧(V)-電流(I)の関係またはそれを結んだI-V カーブ,P-V カーブを作成し、該カーブを使用すること
を特徴とする太陽光発電量評価方法。
【請求項５】
対象太陽光発電システムの一定時間（１０分〜１時間程度）の温度補正後システム出力係数（＝（実測出力（ＫＷｈ））/日射量（ＫＷｈ/ m2）/太陽電池面積（m2）＊（２５℃における最大電力（Ｐmax）（ＫＷ））/（測定温度における最大電力（Ｐmax）（ＫＷ））＊（太陽電池定格容量（ＫＷ/ＫＷ/m2））/基準日射量（ＫＷｈ/m2）（＝１）/太陽電池面積（m2）＊１００）の計算において、２５℃におけるＰmaxおよび測定温度おけるＰmaxの計算に必要な、任意の日射強度・太陽電池温度の太陽電池Ｉ−Ｖカーブを作成する方法として、（1）標準太陽電池のＩ−Ｖカーブを変換する方法（「実用的Ｉ−Ｖカーブ作成法」）又は、（2）太陽電池基本式を解く方法（「理論式によるＩ−Ｖカーブ作成法」）を使い算出し、太陽光発電システムの発電量を評価することを特徴とする評価方法
【請求項６】
対象太陽光発電システムの一定時間（１０分〜１時間程度）の温度補正後システム出力係数（＝（実測出力（ＫＷｈ））/日射量（ＫＷｈ/ m2）/太陽電池面積（m2）＊（２５℃における最大電力（Ｐmax）（ＫＷ））/（測定温度における最大電力（Ｐmax）（ＫＷ））＊（太陽電池定格容量（ＫＷ/ＫＷ/m2））/基準日射量（ＫＷｈ/m2）（＝１）/太陽電池面積（m2）＊１００）の、１〜３時間程度の平均値により評価することを特徴とする評価方法において、２５℃と測定温度におけるＰmaxを算出する過程において、
｛45｝電圧Ｖ、電流Ｉ、日射強度１kW/m2 での光起電流IL、飽和電流温度係数Co、接合定数n、並列抵抗Rsh、直列抵抗Rs、太陽電池モジュール温度T （絶対温度）を含んだ
関数：Func(V,I,IL,Co,n,Rsh,Rs,T) = IL - CoT3 exp(-qEg/nk0T)*(exp( q*(V+Rs*I)/(n*k0*T) )-1) - (V+Rs*I)/ Rsh - I  を作成し、つぎに、
｛46｝該関数Func(V,I,IL,Co,n,Rsh,Rs,T)を変数Vで微分した
関数：Div(V,I,IL,Co,n,Rsh,Rs,T)を作成し、
｛47｝太陽電池の基準状態（モジュール温度Ta(298K(ta=25℃))、日射強度Ea（１kW/m2 ））での仕様値、短絡電流Isca、最適電流Iopa 、最適電圧Vopa 、開放電圧Vocaの点P1(0,Isca),P2(Vopa,Iopa),P3(Voca,0)を選択し、
｛48｝前記関数Func(V,I,IL,Co,n,Rsh,Rs,T)= 0のTに基準状態の温度Ta(298K),直列抵抗Rsに基準温度での仕様値Rsa,および前記P1,P2,P3の点の値を代入し、IL,Co,n,Rshを未知数とする
関係式:Func(0,Isca,IL,Co,n,Rsh,Rsa,Ta)= 0,
関係式:Func(Voca,0,IL,Co,n,Rsh,Rsa,Ta)= 0,
関係式:Func(Vopa,Iopa,IL,Co,n,Rsh,Rsa,Ta)= 0を作成し、
｛49｝前記関数Div(V,I,IL,Co,n,Rsh,Rs,T)= 0 に、基準状態の温度Ta(298K)の直列抵抗Rsに値Rsa および前記点P2の値(Vopa,Iopa) を代入して、IL,Co,n,Rshを未知数とする、
関係式:Div(Vopa,Iopa,IL,Co,n,Rsh,Rsa,Ta)= 0を作成し、つぎに、
｛50｝前記４つの関係式:Func(0,Isca,IL,Co,n,Rsh,Rsa,Ta)= 0,Func(Voca,0,IL,Co,n,Rsh,Rsa,Ta)= 0,Func(Vopa,Iopa,IL,Co,n,Rsh,Rsa,Ta)= 0,Div(Vopa,Iopa,IL,Co,n,Rsh,Rsa,Ta)= 0 を満たす解A(ILa,Coa,na,Rsha)を、非線形解法のプログラムによって算出し、つぎに、
｛51｝太陽電池の基準状態と同様に、第２の太陽電池モジュール温度Tb(313K(tb=40℃))・日射強度Eb（１kW/m2 ）)での仕様値である、短絡電流Iscb、最適電流Iopb 、最適電圧Vopb 、開放電圧Vocbの点P1(0,Iscb),P2(Vopb,Iopb),P3(Vocb,0)を選択し、
｛52｝前記関数Func(V,I,IL,Co,n,Rsh,Rs,T)= 0のTに第２の温度Tb(313K),直列抵抗Rsに該温度での仕様値Rsb,および前記{07}のP1,P2,P3の点の値を代入し、IL,Co,n,Rshを未知数とする
関係式:Func(0,Iscb,IL,Co,n,Rsh,Rsb,Tb)= 0,
関係式:Func(Vocb,0,IL,Co,n,Rsh,Rsb,Tb)= 0,
関係式:Func(Vopb,Iopb,IL,Co,n,Rsh,Rsb,Tb)= 0を作成し、
｛53｝前記関数Div(V,I,IL,Co,n,Rsh,Rs,T)= 0 に、第２の温度Tb(313K)の直列抵抗Rsに第２の温度での値Rsb および前記{07}の点P2の値(Vopb,Iopb) を代入して、IL,Co,n,Rshを未知数とする、
関係式:Div(Vopb,Iopb,IL,Co,n,Rsh,Rsb,Tb)= 0を作成し、つぎに、
｛54｝前記４つの関係式:Func(0,Iscb,IL,Co,n,Rsh,Rsb,Tb)= 0,Func(Vocb,0,IL,Co,n,Rsh,Rsb,Tb)= 0,Func(Vopb,Iopb,IL,Co,n,Rsh,Rsb,Tb)= 0,Div(Vopb,Iopb,IL,Co,n,Rsh,Rsb,Tb)= 0 を満たす解A(Ilb,Cob,nb,Rshb)を、非線形解法のプログラムによって、算出し、つぎに、
｛55｝太陽電池の基準状態と同様に第３の太陽電池モジュール温度Tc(328K(tc=55℃))・日射強度Ec（１kW/m2 ）)での仕様値である、短絡電流Iscc、最適電流Iopc 、最適電圧Vopc 、開放電圧Voccの点P1(0,Iscc),P2(Vopc,Iopc),P3(Vocc,0)を選択し、
｛56｝前記関数Func(V,I,IL,Co,n,Rsh,Rs,T)= 0のTに第３の温度Tc(328K),直列抵抗Rsに該温度での仕様値Rsc,および前記{55}のP1,P2,P3の点の値を代入し、IL,Co,n,Rshを未知数とする
関係式:Func(0,Iscc,IL,Co,n,Rsh,Rsc,Tc)= 0,
関係式:Func(Vocc,0,IL,Co,n,Rsh,Rsc,Tc)= 0,
関係式:Func(Vopc,Iopc,IL,Co,n,Rsh,Rsc,Tc)= 0を作成し、
｛57｝前記関数Div(V,I,IL,Co,n,Rsh,Rs,T)= 0 に、第３の温度Tc(328K)の直列抵抗Rsに第３の温度での値Rsc および前記{55}の点P2の値(Vopc,Iopc) を代入して、IL,Co,n,Rshを未知数とする、
関係式:Div(Vopc,Iopc,IL,Co,n,Rsh,Rsc,Tc)= 0を作成し、つぎに、
｛58｝前記４つの関係式:Func(0,Iscc,IL,Co,n,Rsh,Rsc,Tc)= 0,Func(Vocc,0,IL,Co,n,Rsh,Rsc,Tc)= 0,Func(Vopc,Iopc,IL,Co,n,Rsh,Rsc,Tc)= 0,Div(Vopc,Iopc,IL,Co,n,Rsh,Rsc,Tc)= 0 を満たす解A(Ilc,Coc,nc,Rshc)を、非線形解法のプログラムによって算出し、つぎに、
｛59｝評価すべき太陽電池の実測値日射強度Ej,太陽電池モジュール温度tj(摂氏:絶対温度Tj=tj＋273)およびこの条件での発生電圧Vj−発生電流Ijの各値を取り込み、
｛60｝基準状態での前記温度ta(摂氏25℃:絶対温度Ta(K)=ta(℃)＋273) における前記｛50｝の解A(ILa,Coa,na,Rsha)、前記温度tb(摂氏:Tb=tb＋273) における前記｛54｝の解B(ILb,Cob,nb,Rshb)、前記温度tc(摂氏:Tc=tc＋273) における前記｛58｝の解C(ILc,Coc,nc,Rshc)および、入力値Rsa、Rsb、Rsa、Rsc のそれぞれ(IL,Co,n,Rsh,Rs)に関して３点について曲線補間をして、実測温度tj(摂氏:Tj=tj＋273)での特性値M(ILm,Com,nm,Rshm,Rsm) を算出し、つぎに、
｛61｝ILmを実測された日射強度EjによりIL'm=Ilm*Ej/Ea により補正した後、前記関係式：Func(V,I,IL,Co,n,Rsh,Rs,T)=0にIL'm,Com,nm,Rshm,Rsmを代入して、Func(V,I,IL'm,Com,nm,Rshm,Rsm,Tj)=0 を作成し、電圧(V)-電流(I)の関係(約40〜50点)を非線形解法のプログラムによって求め、電圧(V)-電流(I)の関係またはそれを結んだI-V カーブ,P-V カーブを作成し、該カーブを使用すること
を特徴とする太陽光発電量評価方法。
【請求項７】
対象太陽光発電システムの一定時間（１０分〜１時間程度）の温度補正後システム出力係数（＝（実測出力（ＫＷｈ））/日射量（ＫＷｈ/ m2）/太陽電池面積（m2）＊（２５℃における最大電力（Ｐmax）（ＫＷ））/（測定温度における最大電力（Ｐmax）（ＫＷ））＊（太陽電池定格容量（ＫＷ/ＫＷ/m2））/基準日射量（ＫＷｈ/m2）（＝１）/太陽電池面積（m2）＊１００）の、１〜３時間程度の平均値により評価することを特徴とする評価方法において、２５℃と測定温度におけるＰmaxを算出する過程で必要なＩ-Ｖカーブ作成において、
｛62｝電圧Ｖ、電流Ｉ、日射強度１kW/m2 での光起電流IL、飽和電流温度係数Co、接合定数n、並列抵抗Rsh、直列抵抗Rs、太陽電池モジュール温度T （絶対温度）を含んだ
関数：Func(V,I,IL,Co,n,Rsh,Rs,T) = IL - CoT3 exp(-qEg/nk0T)*(exp( q*(V+Rs*I)/(n*k0*T) )-1) - (V+Rs*I)/ Rsh - I  を作成し、つぎに、
｛63｝該関数Func(V,I,IL,Co,n,Rsh,Rs,T)を変数Vで微分した
関数：Div(V,I,IL,Co,n,Rsh,Rs,T)を作成し、
｛64｝太陽電池の基準状態（モジュール温度Ta(298K(ta=25℃))、日射強度Ea（１kW/m2 ））での仕様値、短絡電流Isca、最適電流Iopa 、最適電圧Vopa 、開放電圧Vocaの点P1(0,Isca),P2(Vopa,Iopa),P3(Voca,0)およびこれら３点とは近接しない任意のP4(V4a,I4a)を選択し、
｛65｝前記関数Func(V,I,IL,Co,n,Rsh,Rs,T)= 0のTに基準状態の温度Ta(298K),前記P1,P2,P3,P4の点の値を代入し、IL,Co,n,Rsh,Rsを未知数とする
関係式:Func(0,Isca,IL,Co,n,Rsh,Rs,Ta)= 0,
関係式:Func(Voca,0,IL,Co,n,Rsh,Rs,Ta)= 0,
関係式:Func(Vopa,Iopa,IL,Co,n,Rsh,Rs,Ta)= 0、
関係式:Func(V4a,I4a,IL,Co,n,Rsh,Rs,Ta)= 0を作成し、
｛66｝前記関数Div(V,I,IL,Co,n,Rsh,Rs,T)= 0 に、基準状態の温度Ta(298K)の前記点P2の値(Vopa,Iopa) を代入して、IL,Co,n,Rsh,Rsを未知数とする、
関係式:Div(Vopa,Iopa,IL,Co,n,Rsh,Rs,Ta)= 0を作成し、つぎに、
｛67｝前記５つの関係式:Func(0,Isca,IL,Co,n,Rsh,Rs,Ta)= 0,Func(Voca,0,IL,Co,n,Rsh,Rs,Ta)= 0,Func(Vopa,Iopa,IL,Co,n,Rsh,Rs,Ta)= 0,Div(Vopa,Iopa,IL,Co,n,Rsh,Rs,Ta)= 0、Func(V4a,I4a,IL,Co,n,Rsh,Rs,Ta)= 0 を満たす解A(ILa,Coa,na,Rsha,Rsa)を、非線形解法のプログラムによって算出し、つぎに、
｛68｝太陽電池の基準状態と同様に、第２の太陽電池モジュール温度Tb(313K(tb=40℃))・日射強度Eb（１kW/m2 ）)での仕様値である、短絡電流Iscb、最適電流Iopb 、最適電圧Vopb 、開放電圧Vocbの点P1(0,Iscb),P2(Vopb,Iopb),P3(Vocb,0), P4(V4b,I4b)を選択し、
｛69｝前記関数Func(V,I,IL,Co,n,Rsh,Rs,T)= 0のTに第２の温度Tb(313K),および前記{51}のP1,P2,P3,P4の点の値を代入し、IL,Co,n,Rshを未知数とする
関係式:Func(0,Iscb,IL,Co,n,Rsh,Rs,Tb)= 0,
関係式:Func(Vocb,0,IL,Co,n,Rsh,Rs,Tb)= 0,
関係式:Func(Vopb,Iopb,IL,Co,n,Rsh,Rs,Tb)= 0、
関係式:Func(V4b,I4b ,IL,Co,n,Rsh,Rs,Tb)= 0を作成し、
｛70｝前記関数Div(V,I,IL,Co,n,Rsh,Rs,T)= 0 に、第２の温度Tb(313K)の前記{68}の点P2の値(Vopb,Iopb) を代入して、IL,Co,n,Rsh,Rsを未知数とする、
関係式:Div(Vopb,Iopb,IL,Co,n,Rsh,Rs,Tb)= 0を作成し、つぎに、
｛71｝前記５つの関係式:Func(0,Iscb,IL,Co,n,Rsh,Rs,Tb)= 0,Func(Vocb,0,IL,Co,n,Rsh,Rs,Tb)= 0,Func(Vopb,Iopb,IL,Co,n,Rsh,Rs,Tb)= 0,Div(Vopb,Iopb,IL,Co,n,Rsh,Rs,Tb)= 0, Func(V4b,I4b ,IL,Co,n,Rsh,Rs,Tb)= 0 を満たす解A(Ilb,Cob,nb,Rshb,Rsb)を、非線形解法のプログラムによって、算出し、つぎに、
｛72｝太陽電池の基準状態と同様に第３の太陽電池モジュール温度Tc(328K(tc=55℃))・日射強度Ec（１kW/m2 ）)での仕様値である、短絡電流Iscc、最適電流Iopc 、最適電圧Vopc 、開放電圧Voccの点P1(0,Iscc),P2(Vopc,Iopc),P3(Vocc,0), P4(V4c,I4c)を選択し、
｛73｝前記関数Func(V,I,IL,Co,n,Rsh,Rs,T)= 0のTに第３の温度Tc(328K),前記{11}のP1,P2,P3,P4の点の値を代入し、IL,Co,n,Rsh,Rsを未知数とする
関係式:Func(0,Iscc,IL,Co,n,Rsh,Rs,Tc)= 0,
関係式:Func(Vocc,0,IL,Co,n,Rsh,Rs,Tc)= 0,
関係式:Func(Vopc,Iopc,IL,Co,n,Rsh,Rs,Tc)= 0,
関係式:Func(V4b,I4b,IL,Co,n,Rsh,Rs,Tc)= 0を作成し、
｛74｝前記関数Div(V,I,IL,Co,n,Rsh,Rs,T)= 0 に、第３の温度Tc(328K)の前記{72}の点P2の値(Vopc,Iopc) を代入して、IL,Co,n,Rsh,Rsを未知数とする、
関係式:Div(Vopc,Iopc,IL,Co,n,Rsh,Rs,Tc)= 0を作成し、つぎに、
｛75｝前記５つの関係式:Func(0,Iscc,IL,Co,n,Rsh,Rs,Tc)= 0,Func(Vocc,0,IL,Co,n,Rsh,Rs,Tc)= 0,Func(Vopc,Iopc,IL,Co,n,Rsh,Rs,Tc)= 0,Div(Vopc,Iopc,IL,Co,n,Rsh,Rs,Tc)= 0, Func(V4b,I4b,IL,Co,n,Rsh,Rs,Tc)= 0 を満たす解A(Ilc,Coc,nc,Rshc,Rsc)を、非線形解法のプログラムによって算出し、つぎに、
｛76｝評価すべき太陽電池の実測値日射強度Ej,太陽電池モジュール温度tj(摂氏:絶対温度Tj=tj＋273)およびこの条件での発生電圧Vj−発生電流Ijの各値を取り込み、
｛77｝基準状態での前記温度ta(摂氏25℃:絶対温度Ta(K)=ta(℃)＋273) における前記｛67｝の解A(ILa,Coa,na,Rsha,Rsa)、前記温度tb(摂氏:Tb=tb＋273) における前記｛71｝の解B(ILb,Cob,nb,Rshb,Rsb)、前記温度tc(摂氏:Tc=tc＋273) における前記｛75｝の解C(ILc,Coc,nc,Rshc,Rsc)のそれぞれ(IL,Co,n,Rsh,Rs)に関して３点について曲線補間をして、実測温度tj(摂氏:Tj=tj＋273)での特性値M(ILm,Com,nm,Rshm,Rsm) を算出し、つぎに、
｛78｝ILmを実測された日射強度EjによりIL'm=Ilm*Ej/Ea により補正した後、前記関係式：Func(V,I,IL,Co,n,Rsh,Rs,T)=0にIL'm,Com,nm,Rshm,Rsmを代入して、Func(V,I,IL'm,Com,nm,Rshm,Rsm,Tj)=0 を作成し、電圧(V)-電流(I)の関係(約40〜50点)を非線形解法のプログラムによって求め、電圧(V)-電流(I)の関係またはそれを結んだI-V カーブ,P-V カーブを作成し、該カーブを使用すること
を特徴とする太陽光発電量評価方法。
【請求項８】
対象太陽光発電システムの一定時間（１０分〜１時間程度）の温度補正後システム出力係数（＝（実測出力（ＫＷｈ））/日射量（ＫＷｈ/ m2）/太陽電池面積（m2）＊（２５℃における最大電力（Ｐmax）（ＫＷ））/（測定温度における最大電力（Ｐmax）（ＫＷ））＊（太陽電池定格容量（ＫＷ/ＫＷ/m2））/基準日射量（ＫＷｈ/m2）（＝１）/太陽電池面積（m2）＊１００）の、１〜３時間程度の平均値により評価することを特徴とする評価方法において、２５℃と測定温度におけるＰmaxを算出する過程に必要なＩ-Ｖカーブ作成において、
｛76｝電圧Ｖ、電流Ｉ、日射強度１kW/m2 での光起電流IL、飽和電流温度係数Co、接合定数n、並列抵抗Rsh、直列抵抗Rs、太陽電池モジュール温度T （絶対温度）を含んだ
関数：Func(V,I,IL,Co,n,Rsh,Rs,T) = IL - CoT3 exp(-qEg/nk0T)*(exp( q*(V+Rs*I)/(n*k0*T) )-1) - (V+Rs*I)/ Rsh - I  を作成し、つぎに、
｛77｝太陽電池の基準状態（モジュール温度Ta(298K(ta=25℃))、日射強度Ea（１kW/m2 ））でのメーカーから与えられるｎ個の点P1(V1a,I1a),P2(V2a,I2a),P3(V3a,I3a),・・・PN(VNa,INa)を、前記の関係式:Func(V,I,IL,Co,n,Rsh,Rs,T)= 0に入れてｎ個の関係式Func(V,I,IL,Co,n,Rsh,Rs,Ta)= 0を作成し、
｛78｝該ｎ個の関係式に最も適したIla,Coa,na,Rsha,Rsaを最小自乗問題解法などの数学的な手法で解いて、つぎに、
｛79｝太陽電池の基準状態と同様に、第２の太陽電池モジュール温度Tb(313K(tb=40℃))・日射強度Eb（１kW/m2 ）)でのメーカーから与えられるｎ個の点P1(V1b,I1b),P2(V2b,I2b),P3(V3b,I3b),・・・PN(VNb,Inb)を、前記の関係式:Func(V,I,IL,Co,n,Rsh,Rs,T)= 0に入れてｎ個の関係式Func(V,I,IL,Co,n,Rsh,Rs,Tb)= 0を作成し、
｛80｝該ｎ個の関係式に最も適したIlb,Cob,nb,Rshb,Rsbを最小自乗問題解法などの数学的な手法で解いて、つぎに、
｛81｝太陽電池の基準状態と同様に、第３の太陽電池モジュール温度Tb(328K(tb=55℃))・日射強度Eb（１kW/m2 ）)でのメーカーから与えられるｎ個の点P1(V1c,I1c),P2(V2c,I2c),P3(V3c,I3c),・・・PN(VNc,Inc)を、前記の関係式:Func(V,I,IL,Co,n,Rsh,Rs,T)= 0に入れてｎ個の関係式Func(V,I,IL,Co,n,Rsh,Rs,Tc)= 0を作成し、
｛82｝該ｎ個の関係式に最も適したIlc,Coc,nc,Rshc,Rscを最小自乗問題解法などの数学的な手法で解いて、つぎに、
｛83｝評価すべき太陽電池の実測値日射強度Ej,太陽電池モジュール温度tj(摂氏:絶対温度Tj=tj＋273)およびこの条件での発生電圧Vj−発生電流Ijの各値を取り込み、
｛84｝基準状態での前記温度ta(摂氏25℃:絶対温度Ta(K)=ta(℃)＋273) における前記｛78｝の解A(ILa,Coa,na,Rsha,Rsa)、前記温度tb(摂氏:Tb=tb＋273) における前記｛80｝の解B(ILb,Cob,nb,Rshb,Rsb)、前記温度tc(摂氏:Tc=tc＋273) における前記｛82｝の解C(ILc,Coc,nc,Rshc,Rsc) のそれぞれ(IL,Co,n,Rsh,Rs)に関して３点について曲線補間をして、実測温度tj(摂氏:Tj=tj＋273)での特性値M(ILm,Com,nm,Rshm,Rsm) を算出し、つぎに、
｛85｝ILmを実測された日射強度EjによりIL'm=Ilm*Ej/Ea により補正した後、前記関係式：Func(V,I,IL,Co,n,Rsh,Rs,T)=0にIL'm,Com,nm,Rshm,Rsmを代入して、Func(V,I,IL'm,Com,nm,Rshm,Rsm,Tj)=0 を作成し、電圧(V)-電流(I)の関係(約40〜50点)を非線形解法のプログラムによって求め、電圧(V)-電流(I)の関係またはそれを結んだI-V カーブ,P-V カーブを作成し、該カーブを使用すること
を特徴とする太陽光発電量評価方法。
【請求項９】
請求項５、６、７、８において評価に使うデータの条件は、（1）太陽電池温度２５℃近辺で、日射強度が８５０Ｗ/m2程度以上または（2）日射強度が１０００Ｗ/m2
程度で、１〜３時間程度の安定した値が得られたデータなどを使い発電量を評価することを特徴とする太陽光発電量評価方法。
【請求項１０】
請求項１または請求項２または請求項３または請求項４または請求項５または請求項６または請求項７または請求項８または請求項９の太陽光発電評価方法で処理する処理プログラムを記録することを特徴とするコンピュータ読み取り可能な記録媒体。
【請求項１１】
請求項１０記載の記録媒体を動作しうるコンピュータからなることを特徴とする太陽電池出力評価装置・システム。
【発明の詳細な説明】
【技術分野】
【０００１】
太陽電池は、図3のごとく、太陽の光エネルギーを電気エネルギーに直接変換するものである。すなわち光電効果の一種である光起電力効果を応用し、起電力を発生させるものであり、太陽電池中に適当なエネルギーを持った光（光子）が入射すると、自由な電子と正孔が発生する。太陽電池は半導体中のpn接合近傍に達した電子と正孔は、それぞれｎ型半導体側、ｐ型半導体側に拡散し、両電極部に集まるので、電力が取り出せ、電圧および電流が発生するというわけである。
太陽電池は、結晶シリコン系、非結晶シリコン系、化合物系に大別できる。  結晶シリコン系太陽電池、特に単結晶太陽電池は製造工程が複雑で、製造に大量の電力を必要とするため、コスト低減の研究が進められている。最近、素子構造の工夫などにより、20％を超える変換効率を達成した報告もある。製造工程の少し簡単な多結晶シリコン太陽電池では、太陽電池変換効率は実用上では10％〜15％程度である。非結晶シリコン系太陽電池（アモルファス太陽電池）では、製造工程が簡単で製造エネルギーが少なく、またシリコン材料が少なくてすむため、低コスト太陽電池として有望である。さらに、薄膜で各種の基盤上に形成できるため、広い応用範囲が期待される。効率は実用上では10％程度である。
本発明は、このような太陽電池を使った発電システム（太陽光発電システム）において、フィールドに設置した太陽光発電システムの発電出力、発電量が工場での太陽電池の測定値又は仕様値（それぞれの太陽電池の工場試験値又はその型式の太陽電池の定格値）に比較してどの程度であるかなどを短期間に正確かつ的確に評価する方法、およびそのプログラムを記録したコンピュータ読み取り可能な記録媒体とその評価装置・システムに関する。
【背景技術】
【０００２】
太陽光発電システムは、地球環境問題解決の高まりのなか、太陽電池とそのシステム価格の低減、逆潮流に対する電力会社の受入体制の整備、国の助成措置などにより最近の普及が著しい。また、新しい高機能の太陽電池も開発されている。しかし最近では、実際に設置した太陽電池の発電出力や発電量が本来の仕様値ほど得られていないのではないかという疑問を、太陽光発電システムの設置者、設置工事業者などからよく聞く。本発明者らが実際に太陽電池のこれら出力・発電量を調査した結果でも、出力・発電量がかなり小さい事例をよく見かける。
国などの実施するフィールド試験では、「システム出力係数」がよく使われている。「システム出力係数」とは、［ある期間の太陽光発電システムの発電量（KWH）］を［同じ期間のそのシステムに照射された日射量（KWH/ m2）＊基準日射強度での太陽電池定格出力（KW/KW/ m2）］で除した値である。また従来から実施されている方法として、一定期間の発電量を、同一期間に測定した日射計による太陽電池受光面の日射量に太陽電池面積を乗じた値で除した「太陽電池変換効率」を求めて評価する方法もある。また、一定期間の積算値で評価するのでなく、それぞれの時点における測定した太陽電池出力を日射強度・太陽電池温度から概略求めた出力値で除して評価することも行なわれている。しかしこれらの方法では、文字通り概略値の評価であり、評価に長期の期間がかかったり、的確な評価値が得られなかったりする。また、これらの技術では評価精度が低いため、出力・発電量不足がみられても、その原因解明とその対策の実施をすることが難しかった。
【０００３】
そこで最近では、「（太陽電池）I-Vカーブ測定器」が開発され市販されている。この装置では、評価すべき太陽電池を、その接続されたインバータ・負荷などの機器から切離し、模擬負荷（コンデンサ負荷、電子負荷など）を接続し、その負荷の値を高速に自動切替・変化することにより、0.5秒程度以内に40〜50組程度の発生電圧（Ｖ）−発生電流（Ｉ）値を測定して電圧−電流カーブ（I-Vカーブ）を得、同時に当該太陽電池の温度（以下太陽電池温度又はモジュール温度という）と当該太陽電池受光面の日射強度を得て、これらの値を使い基準状態のI-Vカーブにもどし、メーカー仕様値と比較して、太陽電池の出力を評価するものである。すなわち、この装置の評価の方法は、評価すべき時刻の日射強度・太陽電池温度において測定したI-Vカーブを基準状態（日射強度１kW／ｍ2、太陽電池温度25℃）のI-Vカーブに変換式を使って変換し、太陽電池メーカーが示している仕様値（基準状態における短絡電流（Isc）、開放電圧（Voc）、最適電圧（Vop）、最適電流（Iop）など）、実際は最大出力電力Pmax（Iop＊Vop）と比較・評価するものである。そして、評価対象の太陽電池の出力（電力）とその評価結果などをパソコン・マイコンで出力するものである。
太陽電池の出力・発電量の評価では、この「I-Vカーブ測定器」が重要な地位を占めているにもかかわらず、この種装置で市販されている機種は国内外に意外に少なく、またこれら機種のハード面の原理および評価ソフトは類似している。一般的に「I-Vカーブ測定器」による評価方法・装置の課題は次のとおりである。
(1) この装置の評価ソフトの中核部分は、測定時の日射強度、太陽電池温度における発生電圧Ｖ−発生電流Ｉの値（I-Vカーブの値）を、基準状態（日射強度１kW／ｍ2、モジュール温度25℃）の電圧−電流値（I-Vカーブの値）に、変換（換算）する方法（変換式）である。従来から使われている変換式では（図５の表の左欄又は中欄の式）であり、変換式を適用できる日射強度の条件（範囲）が0.8kW／ｍ2程度以上で、日射も安定している場合などに限定されており、実際には評価ができない日射強度条件の場合も多かった。すなわち、汎用的な精度の高い変換式が一般には確立していなかった。
(2) 太陽電池や太陽光発電システムのそれぞれの時刻における出力の評価の他に、一定期間の発電量すなわちKWhで正確に評価する方法が確立していない。
(3) 一般に日射の変動は激しく、0.5秒以内に最高から最低まで日射強度が変動することもみられる。したがって、より高速に電圧−電流の値（I-Vカーブ）を測定する技術が必要であり、この点でも従来の技術では十分とはいえない。
(4) 太陽電池出力を評価する場合、運転中の太陽電池を連系されている太陽電池から切離すことなく太陽電池の出力・発電量を計測し、評価することが望まれる。
【０００４】
本発明では以上の事実を考慮しながら、太陽光発電システムの発電量評価が、太陽電池を太陽光発電システムから切離すことなく計測でき、発電量評価することができる方法を中心に考えるため、ここでは「Ｉ−Ｖカーブ測定器」は使用しない方法とする。すなわち、連系状態で評価する方法とする。すなわち、実際の太陽光発電システムは常に運転状態にあり、停止して系統から切離して出力・発電量を評価することは難しい面が多い。そのため、「I-Vカーブ測定器」を使わないで実施できる評価方法が望まれている。一般に、太陽光発電システムを連系状態のままで評価する方法としては前述のように、ある時刻のシステムの出力（発電）と、その時の日射強度と太陽電池温度から計算した出力（発電）の比を求めて評価することが基本となる。そのためには、次の主要技術が必要となるが、いずれも未だに確立していない。
【０００５】
太陽電池特性値から任意の日射強度・太陽電池温度におけるＩ−Ｖカーブを作成する技術
このようなＩ−Ｖカーブが描ければ計算出力（発電）が正確に算出でき基本的には太陽光発電システムの発電量の正確な評価が可能となる。
【０００６】
（２）激しい日射変動のもとでも、また日射強度や太陽電池温度の大きさに影響されにくく安定した評価ができる技術
激しい日射変動のもとなどでも安定した評価ができることが必要である。
【０００７】
（３）その他
できるだけ短期間の測定で評価ができること、および太陽光発電システムがＭＰＰＴ運転できているかを評価できる技術も必要である。
【０００８】
ここで評価に関する従来の技術と課題を十分理解するために、太陽電池の基本的な技術を述べるとともに、太陽電池の特性と特性曲線、および評価の方法に関して説明する。
図５は太陽電池の出力特性曲線である、電圧−電流曲線（I-Vカーブ）を示しており、横軸は電圧Ｖ、縦軸は電流Ｉを示している。また破線は電圧−電力曲線（P-Vカーブ）を示しており、横軸は電圧Ｖ、縦軸は電力Ｐ（Ｖ＊Ｉ）である。
太陽電池の出力電力は、電圧Ｖと電流Ｉの積であり、この太陽電池の出力が最大となる点が、最大出力点（Pmax）と呼ばれている。最大出力点Pmaxにおける電流Ｉおよび電圧Ｖのそれぞれの値は、最適電流Iopおよび最適電圧Vopと呼ばれている。そして、電流が零のときの電圧Ｖを解放電圧Vocとまた、電圧が零のときの電流Ｉを短絡電流Iscという。
太陽電池は、その個々の特性の違いにより、また受光する光の日射強度や太陽電池温度により、その発生電圧と発生電流の曲線であるI-Vカーブが異なってくる。そこで太陽電池の特性値は、基準日射強度（１kW／ｍ2）で基準太陽電池温度（25℃）の状態（「基準状態」という）における短絡電流Isc、最適電流Iop、最適電圧Vop、開放電圧Vocおよび、太陽電池直列抵抗Rs、温度が１℃変化したときの短絡電流Iscの変動値α（Ａ／℃）、温度が１℃変化したときの開放電圧Vocの変動値β（Ｖ／℃）、曲線補正因子Ｋで示される。なお、太陽電池の特性値は、日射強度１kW／ｍ2で太陽電池55℃の状態（「動作状態」ともいう）での短絡電流Isc、最適電流Iop、最適電圧Vop、開放電圧Vocなどで示されることもある。
また、太陽電池の出力特性値は、前述のごとく一般に、日射強度が１kW／ｍ2で太陽電池温度が25℃（場合によっては55℃を含むことがある）の状態がメーカーから与えられる。このため、太陽電池温度が25℃や55℃以外の場合や、日射強度が１kW／ｍ2以外の場合には、評価対象の太陽電池の出力や発電量が太陽電池本来の出力や発電量に比較してどの程度であるかわからない。つまり、太陽電池温度が基準温度と異なる値、例えば38℃で、日射強度が基準日射強度と異なる値、例えば700Ｗ／ｍ2の場合には、I-Vカーブのデータは基準状態のI-Vカーブを日射強度700Ｗ／ｍ2、太陽電池温度38℃に変換できなければ評価対象の太陽電池の出力がどの程度の評価値であるかわからない。このように、太陽電池の出力・発電量の評価では、基準状態のI-Vカーブを任意の日射強度、任意の太陽電池温度のI-Vカーブに変換する式、あるいは任意の日射強度、任意の太陽電池温度で測定したI-Vカーブを基準状態（日射強度１kW／ｍ2、太陽電池温度25℃）のI-Vカーブに正確に変換する式が不可欠である。前記のとおり日射強度がごく限られた条件（800〜850W／ｍ2程度以上で安定しているとき）においては、基準状態のI-Vカーブのデータに戻す式として、従来からJIS C 8913などで示されている変換式が一般に適用されている（図４の表の左欄又は中欄の変換式。これらの式はIEAでも示されている式で、米国製の「I-V測定器」にも適用されている。）しかし、このように、ごく限られた日射強度条件などで適用でき、精度に問題がある変換式では太陽電池出力評価には実質的に十分役立たない。また、この変換式の適用にあたって基準状態での特性値であるα，β，Ｒｓ，Ｋを種々の日射強度・太陽電池温度条件のもとでそのまま式の中に使えるかどうかという課題もある。
【０００９】
また、発電電力量の評価に必要な、基準状態の特性値Isc，Iop，Vop，Vocなどから基準状態のI-Vカーブを描く方法も1部にはあったが（図７）、理論式によるI-Vカーブ作成方法はその適用と応用に関しては、十分には未だ確立されていない。
さらに、太陽電池の特性値である直列抵抗Rsに関しては、太陽電池温度による影響が大きいため、太陽電池温度の関数として表わす式が知られている。例えば、(財)日本品質保証機構では、
Rs＝｛１＋3.3717＊10-3（Ｔ−298）＋9.7058＊10-5（Ｔ−298）2｝＊Rs'
の式（図６）により、太陽電池温度25℃（絶対温度298K）の直列抵抗値Rs'から太陽電池温度Ｔ（絶対温度）での直列抵抗Rsの値を推定することを提案している。しかし、この式はある限定された太陽電池などを使って測定した実験式であり汎用的な式とは言い難いとともに、部分的な解決の試みでしかない。
【発明の開示】
【発明が解決しようとする課題】
【００１０】
上記「従来の技術」の中で技術の現状と共に一般的な課題についても述べたが、ここでは、本発明が解決しようとする課題を中心に整理して述べる。
その第１は、前述のように測定した日射強度・太陽電池温度における太陽電池の電圧−電流カーブ（I-Vカーブ）を基準状態（日射強度１kW／ｍ2、太陽電池温度25℃）に変換して比較評価する正確かつ汎用的な変換式がまだ十分確立してないという課題である。またその逆に、基準状態での特性値（Isc，Iop，Vop，Voc，α，β，Rs，Ｋ）から測定時の日射強度・太陽電池温度における電圧−電流カーブ（I-Vカーブ）を求める正確かつ汎用的な式の確立もできていない。
具体的事例で示すと、現在一般に使用されている任意の日射強度・太陽電池温度条件のI-Vカーブを基準状態のI-Vカーブに変換する変換式（図４の表の左又は中央の欄の変換式）は、前述のように日射強度が800〜850Ｗ／ｍ2程度以下の場合には、計算精度が悪いという欠点などがある。このため、測定時の日射強度が800〜850Ｗ／ｍ2程度より大きい場合でなければ、実測したI-Vカーブを標準状態のI-Vカーブに変換することができないという課題が指摘されている。すなわち、「I-Vカーブ測定器」を使った評価では、測定した測定時日射強度・太陽電池温度条件におけるI-V，P-Vカーブと太陽電池の基本特性値から求めたI-V，P-Vカーブを測定時条件で十分な精度で比較・評価することができなかった。本発明は、この課題を解決するため、精度・汎用性の高いI-Vカーブの作成方法と、これを使った評価の方法をを提案するものである。
【００１１】
第２は、基準状態の特性値（Isc，Iop，Vop，Vocなど）から基準状態のI-Vカーブを描く方法が十分確立していないことである。そのため、測定時日射強度・太陽電池温度条件のI-Vカーブを正確・汎用的に描くことができなかったため、実際に設置された太陽電池の出力やその積算値である発電量を正確・汎用的に評価できなかった。また、その方法が十分確立していないため、太陽光発電システムの年間発電量のシミュレーション計算が正確にできず、太陽光発電システムの設計・運用が十分に実施できにくかったということもある。また、図４の下欄に記載の式にかわる、より理論的な方法を探る動きもあった。例えば、変換式として太陽電池基本式を使い、その基本特性値（Io，Rs，Rshなど）に日射強度・太陽電池温度を使った関数であらわし、この値を基本特性式に適用し解いてI-Vカーブを求める方法である。この方法では、太陽電池の特性値である直列抵抗Rs、並列抵抗Rshなどは太陽電池温度による影響が大きいため、これらの値を太陽電池温度の関数としてあらわしているが、Rs，Rshなどは、前記のとおり特定種類の太陽電池などを使った限られた条件のもとで測定して得られた値である。そのため、汎用的で正確な方法とは言い難く、これらの値を使って求めた計算結果にも課題があった。
他の理論的な方法として、太陽電池基本式を非線形の方程式として解く方法（電気学会論文１（伊賀：「太陽電池の光照射状態での電圧−電流特性を用いたI-Vカーブ作成法とその活用」、電学論116巻、10号、1996））が発表されているが、前述の太陽電池特性値（Isc，Iop，Vop，Voc）から基準状態のI-Vカーブを描く方法、基準状態以外の条件（太陽電池温度40℃、55℃など）のI-Vカーブの作成方法、Rsの求め方、基本特性値（IL，Co，ｎ，Rsh，Rsなど）の曲線補間の方法などにまだ課題があった。
【００１２】
系統と連系したままで太陽光発電システムの評価が安定して的確にできるには、どのような課題があり（Ｐmax運転の評価方法、運転電圧での評価方法など）、どのように解決すべきかという課題があった。
【００１３】
日射変動が激しい等の条件の場合などにおいても、どのような方法なら安定した的確な評価ができるかという課題があった。
【００１４】
また、どのような日射条件と太陽電池温度条件の日・時なら安定した評価ができるかも明らかでなかった。
【００１５】
太陽光発電システムの効率を向上させるため、上記に評価方法をどのように適用して諸損失を分別・解明するかという課題もあった。
【課題を解決するための手段】
【００１６】
請求項１記載の太陽光発電システムの発電量評価方法は、
対象太陽光発電システムの月ごとのシステム出力係数（％）（＝（月間実測システム発電量（ＫＷｈ））/ (月間受光面日射量（ＫＷｈ/ m2）)/太陽電池定格容量（ＫＷ/ＫＷ/ m2）＊１００）を算出し、
そして、太陽電池の温度上昇による損失（％）（＝（１００−月別温度係数））を算出するために、月別温度係数（％）（＝（月間計算システム発電量（ＫＷｈ））/月間受光面日射量（ＫＷｈ/ m2）/太陽電池定格容量（ＫＷ/ＫＷ/ m2）＊１００）計算に必要な月間計算システム発電量算出に必要な、任意の日射強度・太陽電池温度の太陽電池Ｉ−Ｖカーブを作成する方法として、（１）標準太陽電池のＩ−Ｖカーブを変換する方法（「実用的Ｉ−Ｖカーブ作成法」）又は、（２）太陽電池基本式を解く方法（「理論式によるＩ−Ｖカーブ作成法」）を使い、温度上昇による損失（％）を算出し、
該システム出力係数に該温度上昇による損失を加えた値により太陽光発電システムの発電量を評価することを特徴とする。
【００１７】
請求項２記載の太陽光発電システムの発電量評価方法は、
対象太陽光発電システムの月ごとのシステム出力係数（％）（＝（月間実測システム発電量（ＫＷｈ））/ (月間受光面日射量（ＫＷｈ/ m2）)/太陽電池定格容量（ＫＷ/ＫＷ/ m2）＊１００）を算出し、
そして、太陽電池の温度上昇による損失（％）（＝（１００−月別温度係数））を該システム出力係数に加えた値により、太陽光発電システムの発電量を評価することを特徴とする方法であり、
月別温度係数（％）（＝（月間計算システム発電量（ＫＷｈ））/月間受光面日射量（ＫＷｈ/ m2）/太陽電池定格容量（ＫＷ/ＫＷ/ m2）＊１００）計算に必要な月間計算システム発電量の算出過程において、任意の日射強度・太陽電池温度の発電量算出のためのＩ−Ｖカーブ作成において、
｛01｝電圧Ｖ、電流Ｉ、日射強度１kW/m2 での光起電流IL、飽和電流温度係数Co、接合定数n、並列抵抗Rsh、直列抵抗Rs、太陽電池モジュール温度T （絶対温度）を含んだ
関数：Func(V,I,IL,Co,n,Rsh,Rs,T) = IL - CoT3 exp(-qEg/nk0T)*(exp( q*(V+Rs*I)/(n*k0*T) )-1) - (V+Rs*I)/ Rsh - I  を作成し、つぎに、
｛02｝該関数Func(V,I,IL,Co,n,Rsh,Rs,T)を変数Vで微分した
関数：Div(V,I,IL,Co,n,Rsh,Rs,T)を作成し、
｛03｝太陽電池の基準状態（モジュール温度Ta(298K(ta=25℃))、日射強度Ea（１kW/m2 ））での仕様値、短絡電流Isca、最適電流Iopa 、最適電圧Vopa 、開放電圧Vocaの点P1(0,Isca),P2(Vopa,Iopa),P3(Voca,0)を選択し、
｛04｝前記関数Func(V,I,IL,Co,n,Rsh,Rs,T)= 0のTに基準状態の温度Ta(298K),直列抵抗Rsに基準温度での仕様値Rsa,および前記P1,P2,P3の点の値を代入し、IL,Co,n,Rshを未知数とする
関係式:Func(0,Isca,IL,Co,n,Rsh,Rsa,Ta)= 0,
関係式:Func(Voca,0,IL,Co,n,Rsh,Rsa,Ta)= 0,
関係式:Func(Vopa,Iopa,IL,Co,n,Rsh,Rsa,Ta)= 0を作成し、
｛05｝前記関数Div(V,I,IL,Co,n,Rsh,Rs,T)= 0 に、基準状態の温度Ta(298K)の直列抵抗Rsに値Rsa および前記点P2の値(Vopa,Iopa) を代入して、IL,Co,n,Rshを未知数とする、
関係式:Div(Vopa,Iopa,IL,Co,n,Rsh,Rsa,Ta)= 0を作成し、つぎに、
｛06｝前記４つの関係式:Func(0,Isca,IL,Co,n,Rsh,Rsa,Ta)= 0,Func(Voca,0,IL,Co,n,Rsh,Rsa,Ta)= 0,Func(Vopa,Iopa,IL,Co,n,Rsh,Rsa,Ta)= 0,Div(Vopa,Iopa,IL,Co,n,Rsh,Rsa,Ta)= 0 を満たす解A(ILa,Coa,na,Rsha)を、非線形解法のプログラムによって算出し、つぎに、
｛07｝太陽電池の基準状態と同様に、第２の太陽電池モジュール温度Tb(313K(tb=40℃))・日射強度Eb（１kW/m2 ）)での仕様値である、短絡電流Iscb、最適電流Iopb 、最適電圧Vopb 、開放電圧Vocbの点P1(0,Iscb),P2(Vopb,Iopb),P3(Vocb,0)を選択し、
｛08｝前記関数Func(V,I,IL,Co,n,Rsh,Rs,T)= 0のTに第２の温度Tb(313K),直列抵抗Rsに該温度での仕様値Rsb,および前記{07}のP1,P2,P3の点の値を代入し、IL,Co,n,Rshを未知数とする
関係式:Func(0,Iscb,IL,Co,n,Rsh,Rsb,Tb)= 0,
関係式:Func(Vocb,0,IL,Co,n,Rsh,Rsb,Tb)= 0,
関係式:Func(Vopb,Iopb,IL,Co,n,Rsh,Rsb,Tb)= 0を作成し、
｛09｝前記関数Div(V,I,IL,Co,n,Rsh,Rs,T)= 0 に、第２の温度Tb(313K)の直列抵抗Rsに第２の温度での値Rsb および前記{07}の点P2の値(Vopb,Iopb) を代入して、IL,Co,n,Rshを未知数とする、
関係式:Div(Vopb,Iopb,IL,Co,n,Rsh,Rsb,Tb)= 0を作成し、つぎに、
｛10｝前記４つの関係式:Func(0,Iscb,IL,Co,n,Rsh,Rsb,Tb)= 0,Func(Vocb,0,IL,Co,n,Rsh,Rsb,Tb)= 0,Func(Vopb,Iopb,IL,Co,n,Rsh,Rsb,Tb)= 0,Div(Vopb,Iopb,IL,Co,n,Rsh,Rsb,Tb)= 0 を満たす解A(Ilb,Cob,nb,Rshb)を、非線形解法のプログラムによって、算出し、つぎに、
｛11｝太陽電池の基準状態と同様に第３の太陽電池モジュール温度Tc(328K(tc=55℃))・日射強度Ec（１kW/m2 ）)での仕様値である、短絡電流Iscc、最適電流Iopc 、最適電圧Vopc 、開放電圧Voccの点P1(0,Iscc),P2(Vopc,Iopc),P3(Vocc,0)を選択し、
｛12｝前記関数Func(V,I,IL,Co,n,Rsh,Rs,T)= 0のTに第３の温度Tc(328K),直列抵抗Rsに該温度での仕様値Rsc,および前記{11}のP1,P2,P3の点の値を代入し、IL,Co,n,Rshを未知数とする
関係式:Func(0,Iscc,IL,Co,n,Rsh,Rsc,Tc)= 0,
関係式:Func(Vocc,0,IL,Co,n,Rsh,Rsc,Tc)= 0,
関係式:Func(Vopc,Iopc,IL,Co,n,Rsh,Rsc,Tc)= 0を作成し、
｛13｝前記関数Div(V,I,IL,Co,n,Rsh,Rs,T)= 0 に、第３の温度Tc(328K)の直列抵抗Rsに第３の温度での値Rsc および前記{11}の点P2の値(Vopc,Iopc) を代入して、IL,Co,n,Rshを未知数とする、
関係式:Div(Vopc,Iopc,IL,Co,n,Rsh,Rsc,Tc)= 0を作成し、つぎに、
｛14｝前記４つの関係式:Func(0,Iscc,IL,Co,n,Rsh,Rsc,Tc)= 0,Func(Vocc,0,IL,Co,n,Rsh,Rsc,Tc)= 0,Func(Vopc,Iopc,IL,Co,n,Rsh,Rsc,Tc)= 0,Div(Vopc,Iopc,IL,Co,n,Rsh,Rsc,Tc)= 0 を満たす解A(Ilc,Coc,nc,Rshc)を、非線形解法のプログラムによって算出し、つぎに、
｛15｝評価すべき太陽電池の実測値日射強度Ej,太陽電池モジュール温度tj(摂氏:絶対温度Tj=tj＋273)およびこの条件での発生電圧Vj−発生電流Ijの各値を取り込み、
｛16｝基準状態での前記温度ta(摂氏25℃:絶対温度Ta(K)=ta(℃)＋273) における前記｛06｝の解A(ILa,Coa,na,Rsha)、前記温度tb(摂氏:Tb=tb＋273) における前記｛10｝の解B(ILb,Cob,nb,Rshb)、前記温度tc(摂氏:Tc=tc＋273) における前記｛14｝の解C(ILc,Coc,nc,Rshc)および、入力値Rsa、Rsb、Rsa、Rsc のそれぞれ(IL,Co,n,Rsh,Rs)に関して３点について曲線補間をして、実測温度tj(摂氏:Tj=tj＋273)での特性値M(ILm,Com,nm,Rshm,Rsm) を算出し、つぎに、
｛17｝ILmを実測された日射強度EjによりIL‘m=Ilm*Ej/Ea により補正した後、前記関係式：Func(V,I,IL,Co,n,Rsh,Rs,T)=0にIL‘m,Com,nm,Rshm,Rsmを代入して、Func(V,I,IL‘m,Com,nm,Rshm,Rsm,Tj)=0 を作成し、電圧(V)-電流(I)の関係(約40〜50点)を非線形解法のプログラムによって求め、電圧(V)-電流(I)の関係またはそれを結んだI-V カーブ,P-V カーブを作成し、該カーブを使用すること
を特徴とする。
【００１８】
請求項３記載の太陽光発電システムの発電量評価方法は、
対象太陽光発電システムの月ごとのシステム出力係数（％）（＝（月間実測システム発電量（ＫＷｈ））/ (月間受光面日射量（ＫＷｈ/ m2）)/太陽電池定格容量（ＫＷ/ＫＷ/ m2）＊１００）を算出し、
そして、太陽電池の温度上昇による損失（％）（＝（１００−月別温度係数））を算出するために、月別温度係数（％）（＝（月間計算システム発電量（ＫＷｈ））/月間受光面日射量（ＫＷｈ/ m2）/太陽電池定格容量（ＫＷ/ＫＷ/ m2）＊１００）を計算して該システム出力係数に加算した値により評価することに特徴があり、月間計算システム発電量算出に必要な、任意の日射強度・太陽電池発電量算出におけるＩ−Ｖカーブ作成において、
｛18｝電圧Ｖ、電流Ｉ、日射強度１kW/m2 での光起電流IL、飽和電流温度係数Co、接合定数n、並列抵抗Rsh、直列抵抗Rs、太陽電池モジュール温度T （絶対温度）を含んだ
関数：Func(V,I,IL,Co,n,Rsh,Rs,T) = IL - CoT3 exp(-qEg/nk0T)*(exp( q*(V+Rs*I)/(n*k0*T) )-1) - (V+Rs*I)/ Rsh - I  を作成し、つぎに、
｛19｝該関数Func(V,I,IL,Co,n,Rsh,Rs,T)を変数Vで微分した
関数：Div(V,I,IL,Co,n,Rsh,Rs,T)を作成し、
｛20｝太陽電池の基準状態（モジュール温度Ta(298K(ta=25℃))、日射強度Ea（１kW/m2 ））での仕様値、短絡電流Isca、最適電流Iopa 、最適電圧Vopa 、開放電圧Vocaの点P1(0,Isca),P2(Vopa,Iopa),P3(Voca,0)およびこれら３点とは近接しない任意のP4(V4a,I4a)を選択し、
｛21｝前記関数Func(V,I,IL,Co,n,Rsh,Rs,T)= 0のTに基準状態の温度Ta(298K),前記P1,P2,P3,P4の点の値を代入し、IL,Co,n,Rsh,Rsを未知数とする
関係式:Func(0,Isca,IL,Co,n,Rsh,Rs,Ta)= 0,
関係式:Func(Voca,0,IL,Co,n,Rsh,Rs,Ta)= 0,
関係式:Func(Vopa,Iopa,IL,Co,n,Rsh,Rs,Ta)= 0、
関係式:Func(V4a,I4a,IL,Co,n,Rsh,Rs,Ta)= 0を作成し、
｛22｝前記関数Div(V,I,IL,Co,n,Rsh,Rs,T)= 0 に、基準状態の温度Ta(298K)の前記点P2の値(Vopa,Iopa) を代入して、IL,Co,n,Rsh,Rsを未知数とする、
関係式:Div(Vopa,Iopa,IL,Co,n,Rsh,Rs,Ta)= 0を作成し、つぎに、
｛23｝前記５つの関係式:Func(0,Isca,IL,Co,n,Rsh,Rs,Ta)= 0,Func(Voca,0,IL,Co,n,Rsh,Rs,Ta)= 0,Func(Vopa,Iopa,IL,Co,n,Rsh,Rs,Ta)= 0,Div(Vopa,Iopa,IL,Co,n,Rsh,Rs,Ta)= 0、Func(V4a,I4a,IL,Co,n,Rsh,Rs,Ta)= 0 を満たす解A(ILa,Coa,na,Rsha,Rsa)を、非線形解法のプログラムによって算出し、つぎに、
｛24｝太陽電池の基準状態と同様に、第２の太陽電池モジュール温度Tb(313K(tb=40℃))・日射強度Eb（１kW/m2 ）)での仕様値である、短絡電流Iscb、最適電流Iopb 、最適電圧Vopb 、開放電圧Vocbの点P1(0,Iscb),P2(Vopb,Iopb),P3(Vocb,0), P4(V4b,I4b)を選択し、
｛25｝前記関数Func(V,I,IL,Co,n,Rsh,Rs,T)= 0のTに第２の温度Tb(313K),および前記{07}のP1,P2,P3,P4の点の値を代入し、IL,Co,n,Rshを未知数とする
関係式:Func(0,Iscb,IL,Co,n,Rsh,Rs,Tb)= 0,
関係式:Func(Vocb,0,IL,Co,n,Rsh,Rs,Tb)= 0,
関係式:Func(Vopb,Iopb,IL,Co,n,Rsh,Rs,Tb)= 0、
関係式:Func(V4b,I4b ,IL,Co,n,Rsh,Rs,Tb)= 0を作成し、
｛26｝前記関数Div(V,I,IL,Co,n,Rsh,Rs,T)= 0 に、第２の温度Tb(313K)の前記{24}の点P2の値(Vopb,Iopb) を代入して、IL,Co,n,Rsh,Rsを未知数とする、
関係式:Div(Vopb,Iopb,IL,Co,n,Rsh,Rs,Tb)= 0を作成し、つぎに、
｛27｝前記５つの関係式:Func(0,Iscb,IL,Co,n,Rsh,Rs,Tb)= 0,Func(Vocb,0,IL,Co,n,Rsh,Rs,Tb)= 0,Func(Vopb,Iopb,IL,Co,n,Rsh,Rs,Tb)= 0,Div(Vopb,Iopb,IL,Co,n,Rsh,Rs,Tb)= 0, Func(V4b,I4b ,IL,Co,n,Rsh,Rs,Tb)= 0 を満たす解A(Ilb,Cob,nb,Rshb,Rsb)を、非線形解法のプログラムによって、算出し、つぎに、
｛28｝太陽電池の基準状態と同様に第３の太陽電池モジュール温度Tc(328K(tc=55℃))・日射強度Ec（１kW/m2 ）)での仕様値である、短絡電流Iscc、最適電流Iopc 、最適電圧Vopc 、開放電圧Voccの点P1(0,Iscc),P2(Vopc,Iopc),P3(Vocc,0), P4(V4c,I4c)を選択し、
｛29｝前記関数Func(V,I,IL,Co,n,Rsh,Rs,T)= 0のTに第３の温度Tc(328K),前記{28}のP1,P2,P3,P4の点の値を代入し、IL,Co,n,Rsh,Rsを未知数とする
関係式:Func(0,Iscc,IL,Co,n,Rsh,Rs,Tc)= 0,
関係式:Func(Vocc,0,IL,Co,n,Rsh,Rs,Tc)= 0,
関係式:Func(Vopc,Iopc,IL,Co,n,Rsh,Rs,Tc)= 0,
関係式:Func(V4b,I4b,IL,Co,n,Rsh,Rs,Tc)= 0を作成し、
｛30｝前記関数Div(V,I,IL,Co,n,Rsh,Rs,T)= 0 に、第３の温度Tc(328K)の前記{28}の点P2の値(Vopc,Iopc) を代入して、IL,Co,n,Rsh,Rsを未知数とする、
関係式:Div(Vopc,Iopc,IL,Co,n,Rsh,Rs,Tc)= 0を作成し、つぎに、
｛31｝前記５つの関係式:Func(0,Iscc,IL,Co,n,Rsh,Rs,Tc)= 0,Func(Vocc,0,IL,Co,n,Rsh,Rs,Tc)= 0,Func(Vopc,Iopc,IL,Co,n,Rsh,Rs,Tc)= 0,Div(Vopc,Iopc,IL,Co,n,Rsh,Rs,Tc)= 0, Func(V4b,I4b,IL,Co,n,Rsh,Rs,Tc)= 0 を満たす解A(Ilc,Coc,nc,Rshc,Rsc)を、非線形解法のプログラムによって算出し、つぎに、
｛32｝評価すべき太陽電池の実測値日射強度Ej,太陽電池モジュール温度tj(摂氏:絶対温度Tj=tj＋273)およびこの条件での発生電圧Vj−発生電流Ijの各値を取り込み、
｛33｝基準状態での前記温度ta(摂氏25℃:絶対温度Ta(K)=ta(℃)＋273) における前記｛23｝の解A(ILa,Coa,na,Rsha,Rsa)、前記温度tb(摂氏:Tb=tb＋273) における前記｛27｝の解B(ILb,Cob,nb,Rshb,Rsb)、前記温度tc(摂氏:Tc=tc＋273) における前記｛31｝の解C(ILc,Coc,nc,Rshc,Rsc)のそれぞれ(IL,Co,n,Rsh,Rs)に関して３点について曲線補間をして、実測温度tj(摂氏:Tj=tj＋273)での特性値M(ILm,Com,nm,Rshm,Rsm) を算出し、つぎに、
｛34｝ILmを実測された日射強度EjによりIL‘m=Ilm*Ej/Ea により補正した後、前記関係式：Func(V,I,IL,Co,n,Rsh,Rs,T)=0にIL‘m,Com,nm,Rshm,Rsmを代入して、Func(V,I,IL‘,Com,nm,Rshm,Rsm,Tj)=0 を作成し、電圧(V)-電流(I)の関係(約40〜50点)を非線形解法のプログラムによって求め、電圧(V)-電流(I)の関係またはそれを結んだI-V カーブ,P-V カーブを作成し、該カーブを使用すること
を特徴とする。
【００１９】
請求項４記載の太陽光発電システムの発電量評価方法は、
対象太陽光発電システムの月ごとのシステム出力係数（％）（＝（月間実測システム発電量（ＫＷｈ））/ (月間受光面日射量（ＫＷｈ/ m2）)/太陽電池定格容量（ＫＷ/ＫＷ/ m2）＊１００）を算出し、
そして、太陽電池の温度上昇による損失（％）（＝（１００−月別温度係数））を算出するために、月別温度係数（％）（＝（月間計算システム発電量（ＫＷｈ））/月間受光面日射量（ＫＷｈ/ m2）/太陽電池定格容量（ＫＷ/ＫＷ/ m2）＊１００）を計算し、該システム出力係数に加算した値により評価することを特徴とする方法において、月間計算システム発電量算出に必要な、任意の日射強度・太陽電池温度の発電量算出のためのＩ−Ｖカーブ作成において、
｛35｝電圧Ｖ、電流Ｉ、日射強度１kW/m2 での光起電流IL、飽和電流温度係数Co、接合定数n、並列抵抗Rsh、直列抵抗Rs、太陽電池モジュール温度T （絶対温度）を含んだ
関数：Func(V,I,IL,Co,n,Rsh,Rs,T) = IL - CoT3 exp(-qEg/nk0T)*(exp( q*(V+Rs*I)/(n*k0*T) )-1) - (V+Rs*I)/ Rsh - I  を作成し、つぎに、
｛36｝太陽電池の基準状態（モジュール温度Ta(298K(ta=25℃))、日射強度Ea（１kW/m2 ））でのメーカーから与えられるｎ個の点P1(V1a,I1a),P2(V2a,I2a),P3(V3a,I3a),・・・PN(VNa,INa)を、前記の関係式:Func(V,I,IL,Co,n,Rsh,Rs,T)= 0に入れてｎ個の関係式Func(V,I,IL,Co,n,Rsh,Rs,Ta)= 0を作成し、
｛37｝該ｎ個の関係式に最も適したIla,Coa,na,Rsha,Rsaを最小自乗問題解法などの数学的な手法で解いて、つぎに、
｛38｝太陽電池の基準状態と同様に、第２の太陽電池モジュール温度Tb(313K(tb=40℃))・日射強度Eb（１kW/m2 ）)でのメーカーから与えられるｎ個の点P1(V1b,I1b),P2(V2b,I2b),P3(V3b,I3b),・・・PN(VNb,Inb)を、前記の関係式:Func(V,I,IL,Co,n,Rsh,Rs,T)= 0に入れてｎ個の関係式Func(V,I,IL,Co,n,Rsh,Rs,Tb)= 0を作成し、
｛39｝該ｎ個の関係式に最も適したIlb,Cob,nb,Rshb,Rsbを最小自乗問題解法などの数学的な手法で解いて、つぎに、
｛40｝太陽電池の基準状態と同様に、第３の太陽電池モジュール温度Tb(328K(tb=55℃))・日射強度Eb（１kW/m2 ）)でのメーカーから与えられるｎ個の点P1(V1c,I1c),P2(V2c,I2c),P3(V3c,I3c),・・・PN(VNc,Inc)を、前記の関係式:Func(V,I,IL,Co,n,Rsh,Rs,T)= 0に入れてｎ個の関係式Func(V,I,IL,Co,n,Rsh,Rs,Tc)= 0を作成し、
｛41｝該ｎ個の関係式に最も適したIlc,Coc,nc,Rshc,Rscを最小自乗問題解法などの数学的な手法で解いて、つぎに、
｛42｝評価すべき太陽電池の実測値日射強度Ej,太陽電池モジュール温度tj(摂氏:絶対温度Tj=tj＋273)およびこの条件での発生電圧Vj−発生電流Ijの各値を取り込み、
｛43｝基準状態での前記温度ta(摂氏25℃:絶対温度Ta(K)=ta(℃)＋273) における前記｛37｝の解A(ILa,Coa,na,Rsha,Rsa)、前記温度tb(摂氏:Tb=tb＋273) における前記｛39｝の解B(ILb,Cob,nb,Rshb,Rsb)、前記温度tc(摂氏:Tc=tc＋273) における前記｛41｝の解C(ILc,Coc,nc,Rshc,Rsc) のそれぞれ(IL,Co,n,Rsh,Rs)に関して３点について曲線補間をして、実測温度tj(摂氏:Tj=tj＋273)での特性値M(ILm,Com,nm,Rshm,Rsm) を算出し、つぎに、
｛44｝ILmを実測された日射強度EjによりIL‘m=Ilm*Ej/Ea により補正した後、前記関係式：Func(V,I,IL,Co,n,Rsh,Rs,T)=0にIL‘m,Com,nm,Rshm,Rsmを代入して、Func(V,I,IL‘m,Com,nm,Rshm,Rsm,Tj)=0 を作成し、電圧(V)-電流(I)の関係(約40〜50点)を非線形解法のプログラムによって求め、電圧(V)-電流(I)の関係またはそれを結んだI-V カーブ,P-V カーブを作成し、該カーブを使用すること
を特徴とする。
【００２０】
請求項５記載の太陽光発電システムの発電量評価方法は、
対象太陽光発電システムの一定時間（１０分〜１時間程度）の温度補正後システム出力係数（＝（実測出力（ＫＷｈ））/日射量（ＫＷｈ/ m2）/太陽電池面積（m2）＊（２５℃における最大電力（Ｐmax）（ＫＷ））/（測定温度における最大電力（Ｐmax）（ＫＷ））＊（太陽電池定格容量（ＫＷ/ＫＷ/m2））/基準日射量（ＫＷｈ/m2）（＝１）/太陽電池面積（m2）＊１００）の計算において、２５℃におけるＰmaxおよび測定温度おけるＰmaxの計算に必要な、任意の日射強度・太陽電池温度の太陽電池Ｉ−Ｖカーブを作成する方法において、標準太陽電池のＩ−Ｖカーブを変換する方法（「実用的Ｉ−Ｖカーブ作成法」）又は、太陽電池基本式を解く方法（「理論式によるＩ−Ｖカーブ作成法」）を使い算出し、太陽光発電システムの発電量を評価することを特徴とする。
【００２１】
請求項６記載の太陽光発電システムの発電量評価方法は、
対象太陽光発電システムの一定時間（１０分〜１時間程度）の温度補正後システム出力係数（＝（実測出力（ＫＷｈ））/日射量（ＫＷｈ/ m2）/太陽電池面積（m2）＊（２５℃における最大電力（Ｐmax）（ＫＷ））/（測定温度における最大電力（Ｐmax）（ＫＷ））＊（太陽電池定格容量（ＫＷ/ＫＷ/m2））/基準日射量（ＫＷｈ/m2）（＝１）/太陽電池面積（m2）＊１００）の、１〜３時間程度の平均値により評価することを特徴とする評価方法であり、２５℃と測定温度におけるＰmaxを算出する過程で、任意の日射強度・太陽電池温度のＩ−Ｖカーブ作成において、
｛45｝電圧Ｖ、電流Ｉ、日射強度１kW/m2 での光起電流IL、飽和電流温度係数Co、接合定数n、並列抵抗Rsh、直列抵抗Rs、太陽電池モジュール温度T （絶対温度）を含んだ
関数：Func(V,I,IL,Co,n,Rsh,Rs,T) = IL - CoT3 exp(-qEg/nk0T)*(exp( q*(V+Rs*I)/(n*k0*T) )-1) - (V+Rs*I)/ Rsh - I  を作成し、つぎに、
｛46｝該関数Func(V,I,IL,Co,n,Rsh,Rs,T)を変数Vで微分した
関数：Div(V,I,IL,Co,n,Rsh,Rs,T)を作成し、
｛47｝太陽電池の基準状態（モジュール温度Ta(298K(ta=25℃))、日射強度Ea（１kW/m2 ））での仕様値、短絡電流Isca、最適電流Iopa 、最適電圧Vopa 、開放電圧Vocaの点P1(0,Isca),P2(Vopa,Iopa),P3(Voca,0)を選択し、
｛48｝前記関数Func(V,I,IL,Co,n,Rsh,Rs,T)= 0のTに基準状態の温度Ta(298K),直列抵抗Rsに基準温度での仕様値Rsa,および前記P1,P2,P3の点の値を代入し、IL,Co,n,Rshを未知数とする
関係式:Func(0,Isca,IL,Co,n,Rsh,Rsa,Ta)= 0,
関係式:Func(Voca,0,IL,Co,n,Rsh,Rsa,Ta)= 0,
関係式:Func(Vopa,Iopa,IL,Co,n,Rsh,Rsa,Ta)= 0を作成し、
｛49｝前記関数Div(V,I,IL,Co,n,Rsh,Rs,T)= 0 に、基準状態の温度Ta(298K)の直列抵抗Rsに値Rsa および前記点P2の値(Vopa,Iopa) を代入して、IL,Co,n,Rshを未知数とする、
関係式:Div(Vopa,Iopa,IL,Co,n,Rsh,Rsa,Ta)= 0を作成し、つぎに、
｛50｝前記４つの関係式:Func(0,Isca,IL,Co,n,Rsh,Rsa,Ta)= 0,Func(Voca,0,IL,Co,n,Rsh,Rsa,Ta)= 0,Func(Vopa,Iopa,IL,Co,n,Rsh,Rsa,Ta)= 0,Div(Vopa,Iopa,IL,Co,n,Rsh,Rsa,Ta)= 0 を満たす解A(ILa,Coa,na,Rsha)を、非線形解法のプログラムによって算出し、つぎに、
｛51｝太陽電池の基準状態と同様に、第２の太陽電池モジュール温度Tb(313K(tb=40℃))・日射強度Eb（１kW/m2 ）)での仕様値である、短絡電流Iscb、最適電流Iopb 、最適電圧Vopb 、開放電圧Vocbの点P1(0,Iscb),P2(Vopb,Iopb),P3(Vocb,0)を選択し、
｛52｝前記関数Func(V,I,IL,Co,n,Rsh,Rs,T)= 0のTに第２の温度Tb(313K),直列抵抗Rsに該温度での仕様値Rsb,および前記のP1,P2,P3の点の値を代入し、IL,Co,n,Rshを未知数とする
関係式:Func(0,Iscb,IL,Co,n,Rsh,Rsb,Tb)= 0,
関係式:Func(Vocb,0,IL,Co,n,Rsh,Rsb,Tb)= 0,
関係式:Func(Vopb,Iopb,IL,Co,n,Rsh,Rsb,Tb)= 0を作成し、
｛53｝前記関数Div(V,I,IL,Co,n,Rsh,Rs,T)= 0 に、第２の温度Tb(313K)の直列抵抗Rsに第２の温度での値Rsb および前記の点P2の値(Vopb,Iopb) を代入して、IL,Co,n,Rshを未知数とする、
関係式:Div(Vopb,Iopb,IL,Co,n,Rsh,Rsb,Tb)= 0を作成し、つぎに、
｛54｝前記４つの関係式:Func(0,Iscb,IL,Co,n,Rsh,Rsb,Tb)= 0,Func(Vocb,0,IL,Co,n,Rsh,Rsb,Tb)= 0,Func(Vopb,Iopb,IL,Co,n,Rsh,Rsb,Tb)= 0,Div(Vopb,Iopb,IL,Co,n,Rsh,Rsb,Tb)= 0 を満たす解A(Ilb,Cob,nb,Rshb)を、非線形解法のプログラムによって、算出し、つぎに、
｛55｝太陽電池の基準状態と同様に第３の太陽電池モジュール温度Tc(328K(tc=55℃))・日射強度Ec（１kW/m2 ）)での仕様値である、短絡電流Iscc、最適電流Iopc 、最適電圧Vopc 、開放電圧Voccの点P1(0,Iscc),P2(Vopc,Iopc),P3(Vocc,0)を選択し、
｛56｝前記関数Func(V,I,IL,Co,n,Rsh,Rs,T)= 0のTに第３の温度Tc(328K),直列抵抗Rsに該温度での仕様値Rsc,および前記{55}のP1,P2,P3の点の値を代入し、IL,Co,n,Rshを未知数とする
関係式:Func(0,Iscc,IL,Co,n,Rsh,Rsc,Tc)= 0,
関係式:Func(Vocc,0,IL,Co,n,Rsh,Rsc,Tc)= 0,
関係式:Func(Vopc,Iopc,IL,Co,n,Rsh,Rsc,Tc)= 0を作成し、
｛57｝前記関数Div(V,I,IL,Co,n,Rsh,Rs,T)= 0 に、第３の温度Tc(328K)の直列抵抗Rsに第３の温度での値Rsc および前記{55}の点P2の値(Vopc,Iopc) を代入して、IL,Co,n,Rshを未知数とする、
関係式:Div(Vopc,Iopc,IL,Co,n,Rsh,Rsc,Tc)= 0を作成し、つぎに、
｛58｝前記４つの関係式:Func(0,Iscc,IL,Co,n,Rsh,Rsc,Tc)= 0,Func(Vocc,0,IL,Co,n,Rsh,Rsc,Tc)= 0,Func(Vopc,Iopc,IL,Co,n,Rsh,Rsc,Tc)= 0,Div(Vopc,Iopc,IL,Co,n,Rsh,Rsc,Tc)= 0 を満たす解A(Ilc,Coc,nc,Rshc)を、非線形解法のプログラムによって算出し、つぎに、
｛59｝評価すべき太陽電池の実測値日射強度Ej,太陽電池モジュール温度tj(摂氏:絶対温度Tj=tj＋273)およびこの条件での発生電圧Vj−発生電流Ijの各値を取り込み、
｛60｝基準状態での前記温度ta(摂氏25℃:絶対温度Ta(K)=ta(℃)＋273) における前記｛50｝の解A(ILa,Coa,na,Rsha)、前記温度tb(摂氏:Tb=tb＋273) における前記｛54｝の解B(ILb,Cob,nb,Rshb)、前記温度tc(摂氏:Tc=tc＋273) における前記｛58｝の解C(ILc,Coc,nc,Rshc)および、入力値Rsa、Rsb、Rsa、Rsc のそれぞれ(IL,Co,n,Rsh,Rs)に関して３点について曲線補間をして、実測温度tj(摂氏:Tj=tj＋273)での特性値M(ILm,Com,nm,Rshm,Rsm) を算出し、つぎに、
｛61｝ILmを実測された日射強度EjによりIL‘m=Ilm*Ej/Ea により補正した後、前記関係式：Func(V,I,IL,Co,n,Rsh,Rs,T)=0にIL‘m,Com,nm,Rshm,Rsmを代入して、Func(V,I,IL‘m,Com,nm,Rshm,Rsm,Tj)=0 を作成し、電圧(V)-電流(I)の関係(約40〜50点)を非線形解法のプログラムによって求め、電圧(V)-電流(I)の関係またはそれを結んだI-V カーブ,P-V カーブを作成し、該カーブを使用すること
を特徴とする。
【００２２】
請求項７記載の太陽光発電システムの発電量評価方法は、
対象太陽光発電システムの一定時間（１０分〜１時間程度）の温度補正後システム出力係数（＝（実測出力（ＫＷｈ））/日射量（ＫＷｈ/ m2）/太陽電池面積（m2）＊（２５℃における最大電力（Ｐmax）（ＫＷ））/（測定温度における最大電力（Ｐmax）（ＫＷ））＊（太陽電池定格容量（ＫＷ/ＫＷ/m2））/基準日射量（ＫＷｈ/m2）（＝１）/太陽電池面積（m2）＊１００）の、１〜３時間程度の平均値により評価することを特徴とする評価方法であり、２５℃と測定温度におけるＰmaxを算出する過程で、任意の日射強度・太陽電池温度のＩ−Ｖカーブ作成おいて
｛62｝電圧Ｖ、電流Ｉ、日射強度１kW/m2 での光起電流IL、飽和電流温度係数Co、接合定数n、並列抵抗Rsh、直列抵抗Rs、太陽電池モジュール温度T （絶対温度）を含んだ
関数：Func(V,I,IL,Co,n,Rsh,Rs,T) = IL - CoT3 exp(-qEg/nk0T)*(exp( q*(V+Rs*I)/(n*k0*T) )-1) - (V+Rs*I)/ Rsh - I  を作成し、つぎに、
｛63｝該関数Func(V,I,IL,Co,n,Rsh,Rs,T)を変数Vで微分した
関数：Div(V,I,IL,Co,n,Rsh,Rs,T)を作成し、
｛64｝太陽電池の基準状態（モジュール温度Ta(298K(ta=25℃))、日射強度Ea（１kW/m2 ））での仕様値、短絡電流Isca、最適電流Iopa 、最適電圧Vopa 、開放電圧Vocaの点P1(0,Isca),P2(Vopa,Iopa),P3(Voca,0)およびこれら３点とは近接しない任意のP4(V4a,I4a)を選択し、
｛65｝前記関数Func(V,I,IL,Co,n,Rsh,Rs,T)= 0のTに基準状態の温度Ta(298K),前記P1,P2,P3,P4の点の値を代入し、IL,Co,n,Rsh,Rsを未知数とする
関係式:Func(0,Isca,IL,Co,n,Rsh,Rs,Ta)= 0,
関係式:Func(Voca,0,IL,Co,n,Rsh,Rs,Ta)= 0,
関係式:Func(Vopa,Iopa,IL,Co,n,Rsh,Rs,Ta)= 0、
関係式:Func(V4a,I4a,IL,Co,n,Rsh,Rs,Ta)= 0を作成し、
｛66｝前記関数Div(V,I,IL,Co,n,Rsh,Rs,T)= 0 に、基準状態の温度Ta(298K)の前記点P2の値(Vopa,Iopa) を代入して、IL,Co,n,Rsh,Rsを未知数とする、
関係式:Div(Vopa,Iopa,IL,Co,n,Rsh,Rs,Ta)= 0を作成し、つぎに、
｛67｝前記５つの関係式:Func(0,Isca,IL,Co,n,Rsh,Rs,Ta)= 0,Func(Voca,0,IL,Co,n,Rsh,Rs,Ta)= 0,Func(Vopa,Iopa,IL,Co,n,Rsh,Rs,Ta)= 0,Div(Vopa,Iopa,IL,Co,n,Rsh,Rs,Ta)= 0、Func(V4a,I4a,IL,Co,n,Rsh,Rs,Ta)= 0 を満たす解A(ILa,Coa,na,Rsha,Rsa)を、非線形解法のプログラムによって算出し、つぎに、
｛68｝太陽電池の基準状態と同様に、第２の太陽電池モジュール温度Tb(313K(tb=40℃))・日射強度Eb（１kW/m2 ）)での仕様値である、短絡電流Iscb、最適電流Iopb 、最適電圧Vopb 、開放電圧Vocbの点P1(0,Iscb),P2(Vopb,Iopb),P3(Vocb,0), P4(V4b,I4b)を選択し、
｛69｝前記関数Func(V,I,IL,Co,n,Rsh,Rs,T)= 0のTに第２の温度Tb(313K),および前記{51}のP1,P2,P3,P4の点の値を代入し、IL,Co,n,Rshを未知数とする
関係式:Func(0,Iscb,IL,Co,n,Rsh,Rs,Tb)= 0,
関係式:Func(Vocb,0,IL,Co,n,Rsh,Rs,Tb)= 0,
関係式:Func(Vopb,Iopb,IL,Co,n,Rsh,Rs,Tb)= 0、
関係式:Func(V4b,I4b ,IL,Co,n,Rsh,Rs,Tb)= 0を作成し、
｛70｝前記関数Div(V,I,IL,Co,n,Rsh,Rs,T)= 0 に、第２の温度Tb(313K)の前記{68}の点P2の値(Vopb,Iopb) を代入して、IL,Co,n,Rsh,Rsを未知数とする、
関係式:Div(Vopb,Iopb,IL,Co,n,Rsh,Rs,Tb)= 0を作成し、つぎに、
｛71｝前記５つの関係式:Func(0,Iscb,IL,Co,n,Rsh,Rs,Tb)= 0,Func(Vocb,0,IL,Co,n,Rsh,Rs,Tb)= 0,Func(Vopb,Iopb,IL,Co,n,Rsh,Rs,Tb)= 0,Div(Vopb,Iopb,IL,Co,n,Rsh,Rs,Tb)= 0, Func(V4b,I4b ,IL,Co,n,Rsh,Rs,Tb)= 0 を満たす解A(Ilb,Cob,nb,Rshb,Rsb)を、非線形解法のプログラムによって、算出し、つぎに、
｛72｝太陽電池の基準状態と同様に第３の太陽電池モジュール温度Tc(328K(tc=55℃))・日射強度Ec（１kW/m2 ）)での仕様値である、短絡電流Iscc、最適電流Iopc 、最適電圧Vopc 、開放電圧Voccの点P1(0,Iscc),P2(Vopc,Iopc),P3(Vocc,0), P4(V4c,I4c)を選択し、
｛73｝前記関数Func(V,I,IL,Co,n,Rsh,Rs,T)= 0のTに第３の温度Tc(328K),前記のP1,P2,P3,P4の点の値を代入し、IL,Co,n,Rsh,Rsを未知数とする
関係式:Func(0,Iscc,IL,Co,n,Rsh,Rs,Tc)= 0,
関係式:Func(Vocc,0,IL,Co,n,Rsh,Rs,Tc)= 0,
関係式:Func(Vopc,Iopc,IL,Co,n,Rsh,Rs,Tc)= 0,
関係式:Func(V4b,I4b,IL,Co,n,Rsh,Rs,Tc)= 0を作成し、
｛74｝前記関数Div(V,I,IL,Co,n,Rsh,Rs,T)= 0 に、第３の温度Tc(328K)の前記{72}の点P2の値(Vopc,Iopc) を代入して、IL,Co,n,Rsh,Rsを未知数とする、
関係式:Div(Vopc,Iopc,IL,Co,n,Rsh,Rs,Tc)= 0を作成し、つぎに、
｛75｝前記５つの関係式:Func(0,Iscc,IL,Co,n,Rsh,Rs,Tc)= 0,Func(Vocc,0,IL,Co,n,Rsh,Rs,Tc)= 0,Func(Vopc,Iopc,IL,Co,n,Rsh,Rs,Tc)= 0,Div(Vopc,Iopc,IL,Co,n,Rsh,Rs,Tc)= 0, Func(V4b,I4b,IL,Co,n,Rsh,Rs,Tc)= 0 を満たす解A(Ilc,Coc,nc,Rshc,Rsc)を、非線形解法のプログラムによって算出し、つぎに、
｛76｝評価すべき太陽電池の実測値日射強度Ej,太陽電池モジュール温度tj(摂氏:絶対温度Tj=tj＋273)およびこの条件での発生電圧Vj−発生電流Ijの各値を取り込み、
｛77｝基準状態での前記温度ta(摂氏25℃:絶対温度Ta(K)=ta(℃)＋273) における前記｛67｝の解A(ILa,Coa,na,Rsha,Rsa)、前記温度tb(摂氏:Tb=tb＋273) における前記｛71｝の解B(ILb,Cob,nb,Rshb,Rsb)、前記温度tc(摂氏:Tc=tc＋273) における前記｛75｝の解C(ILc,Coc,nc,Rshc,Rsc)のそれぞれ(IL,Co,n,Rsh,Rs)に関して３点について曲線補間をして、実測温度tj(摂氏:Tj=tj＋273)での特性値M(ILm,Com,nm,Rshm,Rsm) を算出し、つぎに、
｛78｝ILmを実測された日射強度EjによりIL‘m=Ilm*Ej/Ea により補正した後、前記関係式：Func(V,I,IL,Co,n,Rsh,Rs,T)=0にIL‘m,Com,nm,Rshm,Rsmを代入して、Func(V,I,IL‘m,Com,nm,Rshm,Rsm,Tj)=0 を作成し、電圧(V)-電流(I)の関係(約40〜50点)を非線形解法のプログラムによって求め、電圧(V)-電流(I)の関係またはそれを結んだI-V カーブ,P-V カーブを作成し、該カーブを使用すること
を特徴とする。
【００２３】
請求項８記載の太陽光発電システムの発電量評価方法は、
対象太陽光発電システムの一定時間（１０分〜１時間程度）の温度補正後システム出力係数（＝（実測出力（ＫＷｈ））/日射量（ＫＷｈ/ m2）/太陽電池面積（m2）＊（２５℃における最大電力（Ｐmax）（ＫＷ））/（測定温度における最大電力（Ｐmax）（ＫＷ））＊（太陽電池定格容量（ＫＷ/ＫＷ/m2））/基準日射量（ＫＷｈ/m2）（＝１）/太陽電池面積（m2）＊１００）の、１〜３時間程度の平均値により評価することを特徴とする評価方法であり、２５℃と測定温度におけるＰmaxを算出する過程で、任意の日射強度・太陽電池温度のＩ−Ｖカーブ作成において、
｛76｝電圧Ｖ、電流Ｉ、日射強度１kW/m2 での光起電流IL、飽和電流温度係数Co、接合定数n、並列抵抗Rsh、直列抵抗Rs、太陽電池モジュール温度T （絶対温度）を含んだ
関数：Func(V,I,IL,Co,n,Rsh,Rs,T) = IL - CoT3 exp(-qEg/nk0T)*(exp( q*(V+Rs*I)/(n*k0*T) )-1) - (V+Rs*I)/ Rsh - I  を作成し、つぎに、
｛77｝太陽電池の基準状態（モジュール温度Ta(298K(ta=25℃))、日射強度Ea（１kW/m2 ））でのメーカーから与えられるｎ個の点P1(V1a,I1a),P2(V2a,I2a),P3(V3a,I3a),・・・PN(VNa,INa)を、前記の関係式:Func(V,I,IL,Co,n,Rsh,Rs,T)= 0に入れてｎ個の関係式Func(V,I,IL,Co,n,Rsh,Rs,Ta)= 0を作成し、
｛78｝該ｎ個の関係式に最も適したIla,Coa,na,Rsha,Rsaを最小自乗問題解法などの数学的な手法で解いて、つぎに、
｛79｝太陽電池の基準状態と同様に、第２の太陽電池モジュール温度Tb(313K(tb=40℃))・日射強度Eb（１kW/m2 ）)でのメーカーから与えられるｎ個の点P1(V1b,I1b),P2(V2b,I2b),P3(V3b,I3b),・・・PN(VNb,Inb)を、前記の関係式:Func(V,I,IL,Co,n,Rsh,Rs,T)= 0に入れてｎ個の関係式Func(V,I,IL,Co,n,Rsh,Rs,Tb)= 0を作成し、
｛80｝該ｎ個の関係式に最も適したIlb,Cob,nb,Rshb,Rsbを最小自乗問題解法などの数学的な手法で解いて、つぎに、
｛81｝太陽電池の基準状態と同様に、第３の太陽電池モジュール温度Tb(328K(tb=55℃))・日射強度Eb（１kW/m2 ）)でのメーカーから与えられるｎ個の点P1(V1c,I1c),P2(V2c,I2c),P3(V3c,I3c),・・・PN(VNc,Inc)を、前記の関係式:Func(V,I,IL,Co,n,Rsh,Rs,T)= 0に入れてｎ個の関係式Func(V,I,IL,Co,n,Rsh,Rs,Tc)= 0を作成し、
｛82｝該ｎ個の関係式に最も適したIlc,Coc,nc,Rshc,Rscを最小自乗問題解法などの数学的な手法で解いて、つぎに、
｛83｝評価すべき太陽電池の実測値日射強度Ej,太陽電池モジュール温度tj(摂氏:絶対温度Tj=tj＋273)およびこの条件での発生電圧Vj−発生電流Ijの各値を取り込み、
｛84｝基準状態での前記温度ta(摂氏25℃:絶対温度Ta(K)=ta(℃)＋273) における前記｛78｝の解A(ILa,Coa,na,Rsha,Rsa)、前記温度tb(摂氏:Tb=tb＋273) における前記｛80｝の解B(ILb,Cob,nb,Rshb,Rsb)、前記温度tc(摂氏:Tc=tc＋273) における前記｛82｝の解C(ILc,Coc,nc,Rshc,Rsc) のそれぞれ(IL,Co,n,Rsh,Rs)に関して３点について曲線補間をして、実測温度tj(摂氏:Tj=tj＋273)での特性値M(ILm,Com,nm,Rshm,Rsm) を算出し、つぎに、
｛85｝ILmを実測された日射強度EjによりIL‘m=Ilm*Ej/Ea により補正した後、前記関係式：Func(V,I,IL,Co,n,Rsh,Rs,T)=0にIL‘m,Com,nm,Rshm,Rsmを代入して、Func(V,I,IL‘m,Com,nm,Rshm,Rsm,Tj)=0 を作成し、電圧(V)-電流(I)の関係(約40〜50点)を非線形解法のプログラムによって求め、電圧(V)-電流(I)の関係またはそれを結んだI-V カーブ,P-V カーブを作成し、該カーブを使用すること
を特徴とする。
【００２４】
請求項９記載の太陽光発電システムの発電量評価方法は、
請求項５、６、７、８において評価に使う測定データの条件は、太陽電池温度２５℃近辺で日射強度が８５０Ｗ/m2程度以上または、日射強度が１０００Ｗ/m2
程度で、１〜３時間程度の安定した値が得られたデータで発電量を評価することを特徴とする。
【００２５】
請求項１０記載のコンピュータ読み取り可能な記憶媒体は、
請求項１または請求項２または請求項３または請求項４または請求項５または請求項６または請求項７または請求項８または請求項９の太陽光発電評価方法で処理する処理プログラムを記録することを特徴とする。
【００２６】
請求項１１記載の太陽光発電評価装置・システムは、
請求項１０記載の記録媒体を動作しうるコンピュータからなることを特徴とする。
【００２７】
ここで本発明に関する用語・機能・技術などについて説明する。
【００２８】
（１）  「太陽電池温度」は「太陽電池モジュール温度」、「モジュール温度」ともいわれ、通常は太陽電池モジュールのセルに埋込んだ熱電対により測定する。なお、太陽電池はその構成の段階により（太陽電池）セル→（太陽電池）モジュール→（太陽電池）アレイという名称で呼ばれている。
【００２９】
（２）   太陽電池温度に使っている記号については、小文字で示したｔ(℃)は摂氏を、大文字で示したＴ(Ｋ(ケルビン))は絶対温度を示す。すなわちＴ(Ｋ)＝ｔ(℃)＋273である。大文字Ｔは主に太陽電池基本式の中で使われており、その他の所では小文字ｔがよく使われている。
【００３０】
（３）  太陽電池出力・発電量に関しては、太陽電池の瞬間の発生電圧Ｖと瞬間の発生電流Ｉの積を太陽電池出力（単位Ｗ又はKW）と言い、その時間積分値が発電量（電力量）（単位Wh又はKWh）である。
【００３１】
（４）  太陽電池基本特性式は次の式である。
Ｉ＝IL−Co＊Ｔ3＊exp（−mqEg/nKoT）＊（exp（q（Ｖ＋RsI）/nKoT）−１）−（Ｖ＋RsI）/Rsh
ここに各記号は次のとおりである。
Ｉ  ：出力電流［Ａ］        Co  ：飽和電流温度係数
Ｖ  ：出力電圧［Ｖ］        Eg  ：エネルギーギャップ［eV］
IL  ：光起電流［Ａ］        Ｔ  ：太陽電池素子絶対温度［Ｋ］
IO  ：飽和電流［Ａ］        Ko  ：ボルツマン定数［J/K］
Rs  ：直流抵抗［Ω］        ｑ  ：電子の電荷量［℃］
Rsh ：並列抵抗［Ω］        ｎ  ：接合定数    
m：１モジュールを構成するセル数（本発明の中では説明省略個所あり）
上記式は半導体の基本特性に基づく式で理論的式である。この式からI-Vカーブ上の各点を求めるには、この式の左辺のＩを右辺に移項し
Func（Ｖ，Ｉ，IL，Co，ｎ，Rsh，Rs，Ｔ）
＝IL−Co＊Ｔ3＊exp（−qEg/nKoT）＊（exp q（Ｖ＋RsI）/nKoT−１）−（Ｖ＋RsI）/Rsh−Ｉ
の関数をつくり、基本特性値（IL，Co，ｎ，Rsh，Rs）を代入のうえ、Ｖ，Ｉの関係を非線形解法プログラムで解いてI-Vカーブを求めることができる。
【００３２】
（５）  本発明では特性値については次のとおり使い分けている。
(ａ) （太陽電池）基本特性値…IL，Co，ｎ，Rsh，Rs
(ｂ) （太陽電池）特性値…Isc，Iop，Vop，Voc，α，β，Rs，Ｋ
なお、上記のとおりRsは(1)，(2)の両方で使われている。
【００３３】
（６）   変換式については次のとおりである。
(ａ) 「実用的I-Vカーブ変換式」
Ｉ1＝Ｉ2＋Isc＊（Ｅ1／Ｅ2−１）＋α＊（ｔ1−ｔ2）
Ｖ1＝Ｖ2＋β＊（ｔ1−ｔ2）−Rs＊（Ｉ1−Ｉ2）−Ｋ＊Ｉ1＊（ｔ1−ｔ2）
(ｂ) 「実用的I-Vカーブ変換式」の逆変換式（逆の適用）
（(1)の式のＶ2，Ｉ2について式を解いて変形したものである）
Ｉ2＝Ｉ1＋Isc＊（Ｅ2−Ｅ1）／Ｅ2＋α（ｔ2−ｔ1）
Ｖ2＝Ｖ1＋β＊（ｔ2−ｔ1）−Rs＊（Ｉ2−Ｉ1）−Ｋ＊Ｉ1＊（ｔ2−ｔ1）
ここで、(1)，(2)は一般に知られているJIS8913，8914，8919の式と形の上では似ているが、以下で述べるように異なる新しい優れた式である。特に(1)の式については特願平6-2626で出願の上、電気学会論文２等で発表している。
また、これらの式で使っている記号は基準状態での電圧値、電流値、日射強度、太陽電池温度をそれぞれＶ2，Ｉ2，Ｅ2，Ｔ2または任意の条件（測定時条件）での各値をＶ1，Ｉ1，Ｅ1，Ｔ1とする。
また、Isc：短絡電流（Ａ）
Iop：最適電流（Ａ）
Vop：最適電圧（Ｖ）
Voc：開放電圧（Ｖ）
α：温度が１℃変化した時のIscの変動値（A/℃）
β：温度が１℃変化した時のVocの変動値（V/℃）
Rs：モジュールの直列抵抗（Ω）
Ｋ：曲線補正因子（Ω/℃）である。
そして、図４の下欄(1)，(2)式が「実用的I-Vカーブ変換式」に、また図４の右欄の(3)，(4)式が「実用的I-Vカーブ変換式」の逆の適用に相当する。
基準状態の電圧−電流値と測定時日射強度・太陽電池温度条件における電圧−電流値の変換式の一覧は前述のとおり図４に示してある。この図は電気学会論文３（伊賀：「実用的I-Vカーブ作成法を使った太陽電池日射計」、電学論Ｄ、117巻10号、1997）に掲載されたもので右端の欄および下欄が発明者の式で下欄は前記のように特願平6-2626で論文発表の前に特許申請がされている。一般には、まだ左端又は中央の欄の式が使われ、前記「I-Vカーブ測定器」でも同様に使われている。
【００３４】
ここで変換式(1)，(2)の適用方法と特徴について、さらに説明を加える。
本発明で使用している(1)，(2)式とJIS補正式（(3)'，(4)'）、CEI IECの式、IEEEの式は一見形の上ではよく似ているように見えるが、図４のように電圧Ｖ、電流Ｉの日射強度・太陽電池温度条件が全く異なっているだけでなくIsc，α，β，Rs，Ｋの日射強度・太陽電池温度条件が全く異なっている。発明者の適用している式（(1)，(2)式）では、電圧、電流をはじめIsc，α，β，Rs，Ｋはすべて基準状態（日射強度１KW／ｍ2、太陽電池温度25℃）における値を使っており、このような適用方式は他に見られない。
そして、このような適用をすることにより、本明細書に記載しているように、精度と汎用性よく測定時条件の電圧−電流特性曲線（I-Vカーブ）、電圧−電力特性曲線（P-Vカーブ）を描け、太陽電池の評価や発電量の予測計算に役立つ効果がある。
そこで、この式（(1)，(2)式）を適用することにより精度と適用性よく、測定時の日射強度・太陽電池温度条件のI-Vカーブ，P-Vカーブが求まるかの理由を次に示す。
（ａ）  太陽電池メーカーから公表されるのは基準状態のパラメータや特性値であり、その他の日射強度・太陽電池温度におけるこれらの値はほとんど示されない。本発明は基準状態の値を使うことを基本としているため、本発明は適用性が広い。
また、他の方法のように、種々の日射強度・太陽電池温度条件のパラメータや特性値を混合して適用することは演算精度・取扱上好ましくない。
（ｂ）  太陽電池の特性を比較・評価する場合、測定値を加工して基準値と比較するより、基準値を測定時条件に変換して測定値と比較する方が直接的で好ましい評価ができる。
【００３５】
太陽光発電システムの出力・発電量の評価をする場合、その測定方法からは、大別して,次の２つに分類できる。
（１）太陽光発電システムを連系運転状態から切離して測定
太陽電池をシステムから切離し「I-Vカーブ測定器」を接続し、模擬負荷抵抗値を高速に変えることにより連続的に発生電圧と発生電流値のｎ組の値（例えば数１０組の値）の関係を測定し、太陽電池の電圧−電流カーブ（I-Vカーブ）を作成し、前記基準状態のＩ−Ｖカーブに変換するなどして、評価する。
（２）太陽光発電システムを連系運転状態のままで電圧・電流を測定
連係状態のままで日射強度、太陽電池温度、出力電圧・電流などを測定して評価する。
本発明では（2）方法に関する発明である。
【発明の効果】
【００３６】
請求項１〜４により算出した、「システム出力係数」＋「太陽電池温度上昇による損失」の値により、対象太陽光発電システムの諸損失が安定かつ正確に把握できるため、諸損失の解明からシステムの効率向上に結びつく。
請求項１では、評価のために必要な太陽電池Ｉ−Ｖカーブ作成方法に種々の方法が使用できるため、適用性の広い汎用性な評価方法といえる。
請求項２により、太陽電池特性値としてメーカーから、短絡電流Isc、最適電流Iop 、最適電圧Vop 、開放電圧Vocのほかに、太陽電池直列抵抗Ｒsが与えられることにより、Ｉ−Vカーブが作成でき、評価ができる。
請求項３、４により、太陽電池特性値としてメーカーから、短絡電流Isc、最適電流Iop 、最適電圧Vop 、開放電圧Vocと、Ｉ−Ｖカーブ上のこれらの点以外のの１点以上の点がが与えられれば、直列抵抗Ｒsが与えられなくてもＩ―Ｖカーブが作成でき,評価ができる。
請求項５〜８は上記の、「システム出力係数」＋「太陽電池温度上昇による損失」の値が直接得られる方法であるとともに、１〜３時間程度の測定値により得られ、短時間の測定値による評価ができる。
請求項５では、評価のために必要な太陽電池Ｉ−Ｖカーブ作成方法に種々の方法が使用できるため、適用性の広い汎用性な評価方法といえる。
請求項６により、太陽電池特性値としてメーカーから、短絡電流Isc、最適電流Iop 、最適電圧Vop 、開放電圧Vocのほかに、太陽電池直列抵抗Ｒsが与えられることにより、Ｉ−Vカーブが作成でき、評価ができる。
請求項７、８により、太陽電池特性値としてメーカーから、短絡電流Isc、最適電流Iop 、最適電圧Vop 、開放電圧Vocと、Ｉ−Ｖカーブ上のこれらの点以外のの１点以上の点がが与えられれば、直列抵抗Ｒsが与えられなくてもＩ―Ｖカーブが作成でき,評価ができる。
請求項９では基準状態に近い日射強度・太陽電池温度条件のデーターを使うことによりメーカーから得られた特性値に誤差があってもその誤差にあまり影響されずに正確な評価が実施できる。
請求項１０、１１により、本評価方法を実施できる装置・システムが製作でき、実用化に結びつけることができる。
【発明を実施するための最良の形態】
【００３７】
つぎに、本発明の実施形態を主な図面に基づき説明する。
図１は太陽光発電システムの発電量評価の例であり、「システム出力係数」に「温度上昇による損失（温度特性部分）」を加えた値（請求項１〜４関係）は年間を通してあまり大きい変化がないことがわかる。すなわち、この値により対象とする太陽光発電システムの諸損失の割合がわかる。その発電量の測定位置を変えることにより、諸損失要素の解明につながる。また、×○印（請求項５〜９関係）は上記値と大差がないことがわかり、日射強度・太陽電池温度条件さえ整えば１〜３時間程度の測定でも評価が可能なことがわかる。
図２は理論式によるＩ−Ｖカーブ作成法のフロー図である。太陽電池温度２５℃の太陽電池特性値（Isc,Vop,Iop,Voc）と２５℃のＲs、α、β、Ｋが与えられた場合、および２５℃、４０℃、５５℃における電池特性値（Isc,Vop,Iop,Voc）と２５℃のＲsが与えられた場合の任意日射強度・温度のＩ−Ｖカーブ，Ｐ−Ｖカーブを作成する方法を示している。なお図では示していないが、太陽電池特性値（Isc,Vop,Iop,Voc）と、Ｉ−Ｖカーブ上のこれらの点以外の１以上の点が与えられれば、同様にＩ−Ｖカーブを作成することができる。
本発明の請求項１〜９は、各種の入力条件から、ニュートン・ラプソン法にもとづく解法によりＩ−Ｖカーブを作成し、評価する方法を示している。また、請求項１および５は「実用的I-Vカーブ作成方法」による方法を含んいる。また、任意の日射強度・太陽電池温度におけるＩ−Vカーブの作成方法であれば、本請求項で記述した方法に限定されずこの評価方法が適用可能である。
図８は請求項５〜９の内容に関して説明したもので、温度補正後システム出力係数（基準（温度）出力係数）の算出方法を示したものである。図ではＩ−Vカーブ作成方法として代表的な「理論式によるＩ−Ｖカーブ作成法」で示しているが、上記のようにいろいろな方法が適用可能である。図の（１）式で基準温度（２５℃）の太陽電池変換効率（「温度補正後変換効率」）と基準温度（２５℃）の出力係数の算出方法を示すものである。
【図面の簡単な説明】
【００３８】
【図１】太陽光発電システムの発電量評価結果の例である。
【図２】理論式によるＩ−Ｖカーブ作成法のフロー図である。
【図３】太陽電池の発電原理である。
【図４】測定時日射強度・太陽電池温度条件におけるI-Vカーブの値を基準状態のI-Vカーブの値への変換式（(3)，(4)式および(3)'、(4)'）およびその逆変換式（図の下欄(1)，(2)式）の説明の図である（電気学会論文３：「実用的I-Vカーブ作成法を使った太陽電池日射計、電学論Ｄ、117巻10号、1997より」。なお、前記のとおりこの表では添字１は測定時日射強度・太陽電池温度条件を添字２は基準状態（日射強度１kw/m2、太陽電池温度２５℃）をあらわしており、本発明本文の中の添字１．２の使い方は逆になっているので注意を要する。
【図５】太陽電池出力特性曲線（I-V，P-Vカーブ）の説明図である。
【図６】(財)日本品質保証機構などが提案する太陽電池基本特性値（Iph，Io，Rs，Rsh）と太陽電池絶対温度（Tp）の関係式である。この式により異なる温度のこれら基本特性値が計算により求まる。
【図７】太陽電池特性値から測定時日射強度・太陽電池温度条件のＩ−Ｖカーブの作成法（「実用的Ｉ−Ｖカーブ作成法」）である。この方法を使うことにより、太陽電池基本式を解かなくても任意の日射強度・太陽電池温度のＩ−Ｖカーブが得られる。
【図８】温度補正後システム出力係数（基準温度の出力係数）の算出フロー図である。
【図１】
2005051048.tif 000003
【図２】
2005051048.tif 000004
【図３】
2005051048.tif 000005
【図４】
2005051048.tif 000006
【図５】
2005051048.tif 000007
【図６】
2005051048.tif 000008
【図７】
2005051048.tif 000009
【図８】
2005051048.tif 000010