Skip to content

Commit f053716

Browse files
kandersolarkandersolar
committed
create sedes2, include in docs
1 parent cc75255 commit f053716

File tree

3 files changed

+206
-0
lines changed

3 files changed

+206
-0
lines changed

‎docs/sphinx/source/reference/effects_on_pv_system_output/spectrum.rst‎

Lines changed: 1 addition & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -7,6 +7,7 @@ Spectrum
77
:toctree: ../generated/
88

99
spectrum.spectrl2
10+
spectrum.sedes2
1011
spectrum.get_example_spectral_response
1112
spectrum.get_reference_spectra
1213
spectrum.calc_spectral_mismatch_field

‎pvlib/spectrum/__init__.py‎

Lines changed: 1 addition & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -1,4 +1,5 @@
11
from pvlib.spectrum.spectrl2 import spectrl2 # noqa: F401
2+
from pvlib.spectrum.sedes2 import sedes2 # noqa: F401
23
from pvlib.spectrum.mismatch import ( # noqa: F401
34
calc_spectral_mismatch_field,
45
spectral_factor_caballero,

‎pvlib/spectrum/sedes2.py‎

Lines changed: 204 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,204 @@
1+
"""
2+
The ``sedes2`` module implements the SEDES2 all-sky spectral irradiance model.
3+
"""
4+
5+
import numpy as np
6+
from scipy.interpolate import make_interp_spline
7+
from pvlib.tools import cosd
8+
9+
SEDES2_COEFFS = np.array([
10+
# wavelength, A1, A2, B1, B2, C1, C2
11+
[320, 1.285724, 0.306791, -0.29613, -0.58516, 0.020632, 0.209150],
12+
[330, 1.235103, 0.262007, -0.28377, -0.53864, 0.010728, 0.206493],
13+
[340, 1.206166, 0.250204, -0.25258, -0.51989, 0.004315, 0.204614],
14+
[350, 1.139737, 0.242676, -0.19222, -0.49821, -0.01184, 0.201329],
15+
[360, 1.091643, 0.244214, -0.13386, -0.48722, -0.02720, 0.200767],
16+
[370, 1.033731, 0.251496, -0.07915, -0.48133, -0.04285, 0.202966],
17+
[380, 0.997179, 0.243862, -0.06550, -0.45039, -0.03607, 0.191915],
18+
[390, 0.997948, 0.227502, -0.08976, -0.40715, -0.01039, 0.173710],
19+
[400, 0.990572, 0.205403, -0.12091, -0.35735, 0.018082, 0.152083],
20+
[410, 0.984024, 0.193105, -0.13671, -0.32748, 0.034395, 0.140702],
21+
[420, 0.971385, 0.177868, -0.15584, -0.29288, 0.051752, 0.127548],
22+
[430, 0.976450, 0.159398, -0.18434, -0.25421, 0.072127, 0.112706],
23+
[440, 0.973204, 0.142079, -0.20773, -0.21836, 0.088689, 0.098569],
24+
[450, 0.979785, 0.129315, -0.22806, -0.19197, 0.103365, 0.087170],
25+
[460, 0.985780, 0.119208, -0.24438, -0.17140, 0.117445, 0.076708],
26+
[470, 0.998610, 0.109176, -0.26163, -0.15113, 0.132599, 0.066066],
27+
[480, 1.005317, 0.099677, -0.27866, -0.13004, 0.147219, 0.055762],
28+
[490, 1.019677, 0.089575, -0.30482, -0.10709, 0.166260, 0.045127],
29+
[500, 1.024404, 0.080517, -0.32229, -0.08750, 0.179513, 0.036474],
30+
[510, 1.031585, 0.069067, -0.34795, -0.06441, 0.196865, 0.025472],
31+
[520, 1.049367, 0.056443, -0.38233, -0.04055, 0.218811, 0.013729],
32+
[530, 1.063939, 0.046316, -0.40907, -0.02121, 0.236117, 0.004201],
33+
[540, 1.071553, 0.038300, -0.42769, -0.00587, 0.248413, -0.00299],
34+
[550, 1.070387, 0.031852, -0.43045, 0.004491, 0.251829, -0.00768],
35+
[560, 1.062834, 0.026342, -0.41879, 0.011999, 0.246651, -0.01046],
36+
[570, 1.045843, 0.024689, -0.37226, 0.009432, 0.223078, -0.00801],
37+
[580, 1.037469, 0.023472, -0.33927, 0.008971, 0.207507, -0.00690],
38+
[590, 1.026083, 0.023298, -0.31410, 0.008151, 0.195734, -0.00518],
39+
[600, 1.040383, 0.015681, -0.34917, 0.024337, 0.218887, -0.01426],
40+
[610, 1.050820, 0.006659, -0.38518, 0.041762, 0.241564, -0.02411],
41+
[620, 1.051636, 0.000294, -0.39171, 0.051032, 0.246386, -0.02902],
42+
[630, 1.040294, -0.00264, -0.36449, 0.050866, 0.230633, -0.02769],
43+
[640, 1.040910, -0.00243, -0.35577, 0.051709, 0.225536, -0.02653],
44+
[650, 1.040678, -0.00316, -0.34746, 0.053757, 0.221069, -0.02611],
45+
[660, 1.065054, -0.00775, -0.38644, 0.068596, 0.246247, -0.03470],
46+
[670, 1.081709, -0.01020, -0.40061, 0.077291, 0.257483, -0.04034],
47+
[680, 1.077241, -0.00697, -0.36968, 0.071587, 0.240555, -0.03716],
48+
[690, 1.040409, -0.00413, -0.28523, 0.052308, 0.187536, -0.02455],
49+
[700, 1.016405, -0.00067, -0.23359, 0.036039, 0.150176, -0.01227],
50+
[710, 1.006519, -0.00416, -0.21335, 0.030742, 0.130581, -0.00725],
51+
[720, 1.015005, -0.00986, -0.20643, 0.033451, 0.120013, -0.00709],
52+
[730, 1.112120, -0.03985, -0.37030, 0.086799, 0.198931, -0.03506],
53+
[740, 1.259638, -0.07938, -0.63633, 0.167885, 0.336038, -0.08023],
54+
[750, 1.359701, -0.10681, -0.82757, 0.227298, 0.435026, -0.11411],
55+
[760, 1.364134, -0.10886, -0.84101, 0.233641, 0.440063, -0.11907],
56+
[770, 1.413504, -0.12491, -0.91952, 0.262684, 0.480399, -0.13497],
57+
[780, 1.472111, -0.14378, -1.00406, 0.291321, 0.524579, -0.14918],
58+
[790, 1.460142, -0.14248, -0.96339, 0.280995, 0.499939, -0.14149],
59+
[800, 1.397082, -0.12613, -0.83251, 0.242547, 0.428312, -0.11892],
60+
[810, 1.303223, -0.09812, -0.64065, 0.184689, 0.325405, -0.08646],
61+
[820, 1.231193, -0.08347, -0.50422, 0.149742, 0.253542, -0.06661],
62+
[830, 1.278968, -0.09801, -0.59564, 0.179143, 0.301940, -0.08288],
63+
[840, 1.394604, -0.12999, -0.82226, 0.248595, 0.424660, -0.12262],
64+
[850, 1.486840, -0.15767, -1.02211, 0.309727, 0.533828, -0.15811],
65+
[860, 1.533058, -0.17332, -1.12535, 0.343352, 0.589583, -0.17738],
66+
[870, 1.548420, -0.17691, -1.14042, 0.350555, 0.597075, -0.18138],
67+
[880, 1.509161, -0.16271, -1.02979, 0.319605, 0.536668, -0.16360],
68+
[890, 1.398192, -0.12470, -0.77108, 0.242984, 0.400871, -0.12150],
69+
[900, 1.176119, -0.06824, -0.34215, 0.120864, 0.176265, -0.05349],
70+
[910, 0.986845, -0.01315, 0.015890, 0.015608, -0.00784, 0.003255],
71+
[920, 0.830406, 0.031590, 0.284686, -0.06127, -0.14019, 0.042905],
72+
[930, 0.611229, 0.097008, 0.607703, -0.15086, -0.28158, 0.082580],
73+
[940, 0.369127, 0.137444, 0.920399, -0.22796, -0.42836, 0.122108],
74+
[950, 0.306382, 0.132255, 1.017929, -0.25108, -0.50619, 0.144856],
75+
[960, 0.427639, 0.084802, 0.857880, -0.20327, -0.46987, 0.132757],
76+
[970, 0.650115, 0.034501, 0.600522, -0.12507, -0.37126, 0.097659],
77+
[980, 0.843689, -0.01411, 0.352464, -0.04375, -0.26576, 0.058199],
78+
[990, 1.018712, -0.05584, 0.115209, 0.032978, -0.16069, 0.019510],
79+
[1000, 1.110714, -0.08242, -0.02662, 0.081820, -0.09732, -0.00507],
80+
[1010, 1.158305, -0.09845, -0.10842, 0.111704, -0.05980, -0.02013],
81+
[1020, 1.187785, -0.10971, -0.17215, 0.134360, -0.02617, -0.03236],
82+
[1030, 1.216623, -0.12039, -0.24681, 0.157773, 0.018206, -0.04635],
83+
[1040, 1.242954, -0.13007, -0.32480, 0.179514, 0.068458, -0.06071],
84+
[1050, 1.242954, -0.13007, -0.32480, 0.179514, 0.068458, -0.06071],
85+
])
86+
87+
88+
def sedes2(spectra_direct_clear, spectra_diffuse_clear, wavelengths,
89+
poa_sky_diffuse, ghi, dni, dhi, solar_zenith, aoi):
90+
"""
91+
Estimate all-sky spectral irradiance using the SEDES2
92+
cloud cover modifiers model.
93+
94+
The SEDES2 model [1]_ estimates all-sky spectral irradiance by
95+
adjusting simulated clear-sky spectra for the effect of cloud cover.
96+
The model was developed for clear-sky spectra simulated with SPECTRL2,
97+
but has since been used with other clear-sky spectral irradiance models
98+
as well.
99+
100+
Parameters
101+
----------
102+
spectra_direct_clear : TYPE
103+
DESCRIPTION.
104+
spectra_diffuse_clear : TYPE
105+
DESCRIPTION.
106+
wavelengths : TYPE
107+
DESCRIPTION.
108+
poa_sky_diffuse : numeric
109+
Plane of array sky diffuse irradiance.
110+
See :term:`poa_sky_diffuse`. [Wm⁻²]
111+
ghi : numeric
112+
Global horizontal irradiance. See :term:`ghi`. [Wm⁻²]
113+
dni : numeric
114+
Direct normal irradiance. See :term:`dni`. [Wm⁻²]
115+
dhi : numeric
116+
Diffuse horizontal irradiance. See :term:`dhi`. [Wm⁻²]
117+
solar_zenith : numeric
118+
Solar zenith angle. See :term:`solar_zenith`. [°]
119+
aoi : numeric
120+
Angle of incidence. See :term:`aoi`. [°]
121+
122+
Returns
123+
-------
124+
out : TYPE
125+
DESCRIPTION.
126+
127+
See also
128+
--------
129+
pvlib.spectrum.spectrl2
130+
131+
Notes
132+
-----
133+
The model lacks a definitive reference, with some implementations
134+
using slightly different coefficients, interpolation schemes,
135+
and normalizations. This implementation follows the description in [2]_.
136+
The coefficients are interpolated to the wavelengths of the input
137+
spectra, extrapolating beyond the range of the coefficients (320-1050 nm)
138+
as needed.
139+
140+
The cloud coverage modifier coefficient values were derived using
141+
equator-facing, fixed-tilt measurements.
142+
143+
References
144+
----------
145+
.. [1] S. Nann and K. Emery, "Spectral effects on PV-device rating,"
146+
Solar Energy Materials and Solar Cells, vol. 27, no. 3,
147+
pp. 189–216, Aug. 1992, :doi:`10.1016/0927-0248(92)90083-2`
148+
.. [2] C. M. Whitaker and J. D. Newmiller, "Photovoltaic module
149+
energy rating procedure. Final subcontract report,". Technical
150+
Report NREL/SR-520-23942, Jan. 1998. :doi:`10.2172/563232`
151+
"""
152+
# Equation numbers refer to Section 2.1.3 of [2]
153+
154+
# arrays are generally 2D, with shape (n_wavelengths, n_timestamps)
155+
156+
poa_sky_diffuse = np.atleast_1d(poa_sky_diffuse)[np.newaxis, :]
157+
ghi = np.atleast_1d(ghi)[np.newaxis, :]
158+
dni = np.atleast_1d(dni)[np.newaxis, :]
159+
dhi = np.atleast_1d(dhi)[np.newaxis, :]
160+
cosZ = np.atleast_1d(cosd(solar_zenith))[np.newaxis, :]
161+
cos_aoi = np.clip(cosd(aoi), a_min=0, a_max=None)
162+
cos_aoi = np.atleast_1d(cos_aoi)[np.newaxis, :]
163+
164+
# interpolate coefficients to match input wavelengths
165+
interp = make_interp_spline(SEDES2_COEFFS[:, 0], SEDES2_COEFFS[:, 1:], k=1)
166+
coef = interp(wavelengths).T
167+
coef = [x[:, np.newaxis] for x in coef] # add dimension for time
168+
A1, A2, B1, B2, C1, C2 = coef
169+
170+
# normalization factors
171+
dhi_clear = np.trapezoid(spectra_diffuse_clear, wavelengths, # 2-44
172+
axis=0)
173+
dni_clear = np.trapezoid(spectra_direct_clear, wavelengths, # 2-45
174+
axis=0)
175+
ndir = (ghi - dhi) / (dni_clear*cosZ) # 2-42
176+
ngh = ghi / (dhi_clear + dni_clear*cosZ) # 2-43
177+
178+
# section 2.1.3.4: rescale spectral irradiance to match broadband values
179+
spectra_direct_clear_rescaled = spectra_direct_clear * ndir # 2-41
180+
# note, [2] is missing a cosZ term here:
181+
spectra_diffuse_clear_rescaled = ( # 2-40
182+
(spectra_direct_clear * cosZ + spectra_diffuse_clear) * ngh
183+
- spectra_direct_clear_rescaled * cosZ
184+
)
185+
186+
# section 2.1.3.3: apply wavelength=dependent cloud cover modifiers
187+
ccm = A1 + A2/cosZ + (B1 + B2/cosZ) * ngh + (C1 + C2/cosZ) * ngh**2 # 2-39
188+
spectra_diffuse = spectra_diffuse_clear_rescaled * ccm # 2-38
189+
spectra_direct = spectra_direct_clear_rescaled * ccm # 3-37
190+
191+
# section 2.1.3.2: convert to POA
192+
spectra_diffuse_poa = spectra_diffuse * poa_sky_diffuse / dhi # 2-36
193+
spectra_direct_poa = spectra_direct * cos_aoi
194+
spectra_global_poa = spectra_direct_poa + spectra_diffuse_poa # 2-35
195+
196+
out = {
197+
'wavelength': wavelengths,
198+
'dhi': spectra_diffuse,
199+
'dni': spectra_direct,
200+
'poa_sky_diffuse': spectra_diffuse_poa,
201+
'poa_direct': spectra_direct_poa,
202+
'poa_global': spectra_global_poa,
203+
}
204+
return out

0 commit comments

Comments
 (0)