Кейс: фотодиод в солнечной панели на крыше частного дома даёт меньше мощности вечером, чем прогнозировалось — какие факторы (угол падения, спектр света, температура, засорение поверхности, деградация материалов) нужно проанализировать и как их количественно учесть?
Кратко — какие факторы анализировать и как их учесть количественно. 1) Геометрия и угол падения - Эффективная падающая радиация на модуль: учитываем косинусный закон и коэффициент уголовой чувствительности (IAM): Geff=Gbeamcosθ⋅IAM(θ)+Gdiffuse⋅Fdiff
G_{\text{eff}} = G_{\text{beam}}\cos\theta\cdot IAM(\theta) + G_{\text{diffuse}}\cdot F_{diff} Geff=Gbeamcosθ⋅IAM(θ)+Gdiffuse⋅Fdiff
где θ\thetaθ — угол падения прямого луча на поверхность, IAM(θ)IAM(\theta)IAM(θ) — измеренная/модельная Incidence Angle Modifier (напр. приближённо IAM(θ)=1−a1(secθ−1)−a2(secθ−1)2IAM(\theta)=1-a_1(\sec\theta-1)-a_2(\sec\theta-1)^2IAM(θ)=1−a1(secθ−1)−a2(secθ−1)2), FdiffF_{diff}Fdiff — фактор для рассеянной составляющей. - Измерения: позиция солнца (время, широта, долгота), измерить или смоделировать θ\thetaθ, измерить или вычислить Gbeam,GdiffuseG_{\text{beam}},G_{\text{diffuse}}Gbeam,Gdiffuse (пиранометр или метео‑данные). 2) Спектр света (спектральный диз‑мэтч) - Квантитативно через спектральную реактивность (EQE/SR) модуля: Isc=q∫0∞Φ(λ) SR(λ) dλ
I_{sc} = q\int_0^\infty \Phi(\lambda)\,SR(\lambda)\,d\lambda Isc=q∫0∞Φ(λ)SR(λ)dλ
и фактор спектрального несоответствия относительно эталона: M=∫E(λ) SR(λ) dλ∫Eref(λ) SR(λ) dλ
M = \frac{\int E(\lambda)\,SR(\lambda)\,d\lambda}{\int E_{ref}(\lambda)\,SR(\lambda)\,d\lambda} M=∫Eref(λ)SR(λ)dλ∫E(λ)SR(λ)dλ
где E(λ)E(\lambda)E(λ) — измеренный спектр (или модель AM). Тогда IscI_{sc}Isc и мощность масштабируются на MMM. - Измерения: портативный спектророметр или использовать эталонный калиброванный приёмник (reference cell) и/или модель AM (вечером спектр изменяется: больше длинноволновой доли — влияние зависит от SR модуля). 3) Температура - Температурный эффект обычно моделируется через температурный коэффициент мощности γP\gamma_PγP (обычно −0.3-0.3−0.3 … −0.5-0.5−0.5%/°C для кремниевых модулей). В практической форме: P(T)=PSTC[1+γP (Tcell−25∘C)]
P(T) = P_{STC}\big[1+\gamma_P\,(T_{cell}-25^\circ C)\big] P(T)=PSTC[1+γP(Tcell−25∘C)]
- Оценка температуры ячейки через NOCT: Tcell≈Tamb+NOCT−20800 G
T_{cell}\approx T_{amb} + \frac{NOCT-20}{800}\,G Tcell≈Tamb+800NOCT−20G
- Можно также через параметры IV: Isc≈Isc,refGGref[1+αI(Tcell−25)],Voc≈Voc,ref+βV(Tcell−25)+nVTlnGGref
I_{sc}\approx I_{sc,ref}\frac{G}{G_{ref}}\big[1+\alpha_I (T_{cell}-25)\big],\quad V_{oc}\approx V_{oc,ref} + \beta_V (T_{cell}-25) + nV_T\ln\frac{G}{G_{ref}} Isc≈Isc,refGrefG[1+αI(Tcell−25)],Voc≈Voc,ref+βV(Tcell−25)+nVTlnGrefG
где αI\alpha_IαI — коэфф. тока (прибл. +0.03+0.03+0.03–0.060.060.06%/°C), βV\beta_VβV — коэфф. напряжения (прибл. −0.25-0.25−0.25–−0.40-0.40−0.40%/°C); точные значения — из даташита. 4) Засорение поверхности (soiling) - Влияет как на прозрачность покрытия: вводят множитель TsoilingT_{soiling}Tsoiling (доля переданной энергии, 0..10..10..1). Garriving=Geff⋅Tsoiling
G_{\text{arriving}} = G_{\text{eff}}\cdot T_{soiling} Garriving=Geff⋅Tsoiling
- Измерить практически: зафиксировать IscI_{sc}Isc до и после очистки; Tsoiling≈Isc,soiled/Isc,cleanT_{soiling}\approx I_{sc,\text{soiled}}/I_{sc,\text{clean}}Tsoiling≈Isc,soiled/Isc,clean при прочих равных условиях. - Можно оценить зависимость по времени (накопление пыли) и скомпенсировать в моделях. 5) Деградация материалов и повреждения - Долгосрочная деградация: годовой процент деградации ddd: P(t)=P0(1−d)t или приближённо P(t)≈P0(1−dt)
P(t)=P_0(1-d)^t\ \text{или приближённо}\ P(t)\approx P_0\big(1-dt\big) P(t)=P0(1−d)tилиприближённоP(t)≈P0(1−dt)
- Острая деградация (PID, SPF, микротрещины, выход из строя диодов) определяется через измерение IV‑характеристики, термографии, эл‑люминесценцию. - Количественно: сравнить номинальные Isc,Voc,FF,PmaxI_{sc},V_{oc},FF,P_{max}Isc,Voc,FF,Pmax с текущими; снижение IscI_{sc}Isc обычно указывает на засор/покрытие/обрыв, снижение VocV_{oc}Voc и FF — на деградацию/повреждение. Как собрать данные и вычислить ожидаемую мощность - Комбинированная модель: P≈A⋅Gglobal⋅M⋅IAM(θ)⋅Tsoiling⋅ηSTC[1+γP (Tcell−25)]
P \approx A\cdot G_{\text{global}}\cdot M\cdot IAM(\theta)\cdot T_{soiling}\cdot \eta_{STC}\big[1+\gamma_P\,(T_{cell}-25)\big] P≈A⋅Gglobal⋅M⋅IAM(θ)⋅Tsoiling⋅ηSTC[1+γP(Tcell−25)]
где AAA — площадь модуля, ηSTC\eta_{STC}ηSTC — КПД при STC. - Или через IV: 1) вычислить IscI_{sc}Isc с учётом G,GG,GG,G-масштабирования, MMM, TsoilingT_{soiling}Tsoiling и температурного коэфф.: формула выше; 2) вычислить VocV_{oc}Voc с температурой и логарифмическим сдвигом от изменения освещённости; 3) оценить FFFFFF (с учётом температуры и деградации) и получить Pmax=ImpVmp≈IscVocFFP_{max}=I_{mp}V_{mp}\approx I_{sc}V_{oc}FFPmax=ImpVmp≈IscVocFF. Практическая последовательность диагностики 1. Измерить профиль освещённости: Gbeam,Gdiffuse,G_{beam},G_{diffuse},Gbeam,Gdiffuse, спектр (или использовать reference cell). 2. Определить геометрию (угол, ориентация) и вычислить θ\thetaθ, IAMIAMIAM. 3. Замерить модульную и амб. температуру, вычислить TcellT_{cell}Tcell. 4. Измерить IV‑кривую (вычленить Isc,Voc,FF,PmaxI_{sc},V_{oc},FF,P_{max}Isc,Voc,FF,Pmax). 5. Очистить панель и повторить замеры для оценки TsoilingT_{soiling}Tsoiling. 6. Сравнить измеренное PmaxP_{max}Pmax с рассчитанным по модели; если разница остаётся — проводить визуальный осмотр, термографию, EL‑съёмку для поиска деградации/шунтов/диодов. Коротко о вечернем эффекте - Вечером важны: меньше G (основная причина), большой угол падения (снижение по косинусу и IAM), изменённый спектр (влияет в зависимости от типа ячейки), иногда понижение температуры (обычно благоприятно для мощности через повышение VocV_{oc}Voc — но проигрывает снижение освещённости), и накопление пыли/росы на поверхности. Количественно оцените все перечисленные факторы по формулам выше и сравните с фактической IV‑кривой. Используйте даташит модуля (температурные коэффициенты, SR/EQE при наличии), портативный пиранометр/спектрометр, термометр/термопару для модуля и инструмент для IV‑замера — тогда расчёты будут точными.
1) Геометрия и угол падения
- Эффективная падающая радиация на модуль: учитываем косинусный закон и коэффициент уголовой чувствительности (IAM):
Geff=Gbeamcosθ⋅IAM(θ)+Gdiffuse⋅Fdiff G_{\text{eff}} = G_{\text{beam}}\cos\theta\cdot IAM(\theta) + G_{\text{diffuse}}\cdot F_{diff}
Geff =Gbeam cosθ⋅IAM(θ)+Gdiffuse ⋅Fdiff где θ\thetaθ — угол падения прямого луча на поверхность, IAM(θ)IAM(\theta)IAM(θ) — измеренная/модельная Incidence Angle Modifier (напр. приближённо IAM(θ)=1−a1(secθ−1)−a2(secθ−1)2IAM(\theta)=1-a_1(\sec\theta-1)-a_2(\sec\theta-1)^2IAM(θ)=1−a1 (secθ−1)−a2 (secθ−1)2), FdiffF_{diff}Fdiff — фактор для рассеянной составляющей.
- Измерения: позиция солнца (время, широта, долгота), измерить или смоделировать θ\thetaθ, измерить или вычислить Gbeam,GdiffuseG_{\text{beam}},G_{\text{diffuse}}Gbeam ,Gdiffuse (пиранометр или метео‑данные).
2) Спектр света (спектральный диз‑мэтч)
- Квантитативно через спектральную реактивность (EQE/SR) модуля:
Isc=q∫0∞Φ(λ) SR(λ) dλ I_{sc} = q\int_0^\infty \Phi(\lambda)\,SR(\lambda)\,d\lambda
Isc =q∫0∞ Φ(λ)SR(λ)dλ и фактор спектрального несоответствия относительно эталона:
M=∫E(λ) SR(λ) dλ∫Eref(λ) SR(λ) dλ M = \frac{\int E(\lambda)\,SR(\lambda)\,d\lambda}{\int E_{ref}(\lambda)\,SR(\lambda)\,d\lambda}
M=∫Eref (λ)SR(λ)dλ∫E(λ)SR(λ)dλ где E(λ)E(\lambda)E(λ) — измеренный спектр (или модель AM). Тогда IscI_{sc}Isc и мощность масштабируются на MMM.
- Измерения: портативный спектророметр или использовать эталонный калиброванный приёмник (reference cell) и/или модель AM (вечером спектр изменяется: больше длинноволновой доли — влияние зависит от SR модуля).
3) Температура
- Температурный эффект обычно моделируется через температурный коэффициент мощности γP\gamma_PγP (обычно −0.3-0.3−0.3 … −0.5-0.5−0.5%/°C для кремниевых модулей). В практической форме:
P(T)=PSTC[1+γP (Tcell−25∘C)] P(T) = P_{STC}\big[1+\gamma_P\,(T_{cell}-25^\circ C)\big]
P(T)=PSTC [1+γP (Tcell −25∘C)] - Оценка температуры ячейки через NOCT:
Tcell≈Tamb+NOCT−20800 G T_{cell}\approx T_{amb} + \frac{NOCT-20}{800}\,G
Tcell ≈Tamb +800NOCT−20 G - Можно также через параметры IV:
Isc≈Isc,refGGref[1+αI(Tcell−25)],Voc≈Voc,ref+βV(Tcell−25)+nVTlnGGref I_{sc}\approx I_{sc,ref}\frac{G}{G_{ref}}\big[1+\alpha_I (T_{cell}-25)\big],\quad
V_{oc}\approx V_{oc,ref} + \beta_V (T_{cell}-25) + nV_T\ln\frac{G}{G_{ref}}
Isc ≈Isc,ref Gref G [1+αI (Tcell −25)],Voc ≈Voc,ref +βV (Tcell −25)+nVT lnGref G где αI\alpha_IαI — коэфф. тока (прибл. +0.03+0.03+0.03–0.060.060.06%/°C), βV\beta_VβV — коэфф. напряжения (прибл. −0.25-0.25−0.25–−0.40-0.40−0.40%/°C); точные значения — из даташита.
4) Засорение поверхности (soiling)
- Влияет как на прозрачность покрытия: вводят множитель TsoilingT_{soiling}Tsoiling (доля переданной энергии, 0..10..10..1).
Garriving=Geff⋅Tsoiling G_{\text{arriving}} = G_{\text{eff}}\cdot T_{soiling}
Garriving =Geff ⋅Tsoiling - Измерить практически: зафиксировать IscI_{sc}Isc до и после очистки; Tsoiling≈Isc,soiled/Isc,cleanT_{soiling}\approx I_{sc,\text{soiled}}/I_{sc,\text{clean}}Tsoiling ≈Isc,soiled /Isc,clean при прочих равных условиях.
- Можно оценить зависимость по времени (накопление пыли) и скомпенсировать в моделях.
5) Деградация материалов и повреждения
- Долгосрочная деградация: годовой процент деградации ddd:
P(t)=P0(1−d)t или приближённо P(t)≈P0(1−dt) P(t)=P_0(1-d)^t\ \text{или приближённо}\ P(t)\approx P_0\big(1-dt\big)
P(t)=P0 (1−d)t или приближённо P(t)≈P0 (1−dt) - Острая деградация (PID, SPF, микротрещины, выход из строя диодов) определяется через измерение IV‑характеристики, термографии, эл‑люминесценцию.
- Количественно: сравнить номинальные Isc,Voc,FF,PmaxI_{sc},V_{oc},FF,P_{max}Isc ,Voc ,FF,Pmax с текущими; снижение IscI_{sc}Isc обычно указывает на засор/покрытие/обрыв, снижение VocV_{oc}Voc и FF — на деградацию/повреждение.
Как собрать данные и вычислить ожидаемую мощность
- Комбинированная модель:
P≈A⋅Gglobal⋅M⋅IAM(θ)⋅Tsoiling⋅ηSTC[1+γP (Tcell−25)] P \approx A\cdot G_{\text{global}}\cdot M\cdot IAM(\theta)\cdot T_{soiling}\cdot \eta_{STC}\big[1+\gamma_P\,(T_{cell}-25)\big]
P≈A⋅Gglobal ⋅M⋅IAM(θ)⋅Tsoiling ⋅ηSTC [1+γP (Tcell −25)] где AAA — площадь модуля, ηSTC\eta_{STC}ηSTC — КПД при STC.
- Или через IV:
1) вычислить IscI_{sc}Isc с учётом G,GG,GG,G-масштабирования, MMM, TsoilingT_{soiling}Tsoiling и температурного коэфф.: формула выше;
2) вычислить VocV_{oc}Voc с температурой и логарифмическим сдвигом от изменения освещённости;
3) оценить FFFFFF (с учётом температуры и деградации) и получить Pmax=ImpVmp≈IscVocFFP_{max}=I_{mp}V_{mp}\approx I_{sc}V_{oc}FFPmax =Imp Vmp ≈Isc Voc FF.
Практическая последовательность диагностики
1. Измерить профиль освещённости: Gbeam,Gdiffuse,G_{beam},G_{diffuse},Gbeam ,Gdiffuse , спектр (или использовать reference cell).
2. Определить геометрию (угол, ориентация) и вычислить θ\thetaθ, IAMIAMIAM.
3. Замерить модульную и амб. температуру, вычислить TcellT_{cell}Tcell .
4. Измерить IV‑кривую (вычленить Isc,Voc,FF,PmaxI_{sc},V_{oc},FF,P_{max}Isc ,Voc ,FF,Pmax ).
5. Очистить панель и повторить замеры для оценки TsoilingT_{soiling}Tsoiling .
6. Сравнить измеренное PmaxP_{max}Pmax с рассчитанным по модели; если разница остаётся — проводить визуальный осмотр, термографию, EL‑съёмку для поиска деградации/шунтов/диодов.
Коротко о вечернем эффекте
- Вечером важны: меньше G (основная причина), большой угол падения (снижение по косинусу и IAM), изменённый спектр (влияет в зависимости от типа ячейки), иногда понижение температуры (обычно благоприятно для мощности через повышение VocV_{oc}Voc — но проигрывает снижение освещённости), и накопление пыли/росы на поверхности. Количественно оцените все перечисленные факторы по формулам выше и сравните с фактической IV‑кривой.
Используйте даташит модуля (температурные коэффициенты, SR/EQE при наличии), портативный пиранометр/спектрометр, термометр/термопару для модуля и инструмент для IV‑замера — тогда расчёты будут точными.