基于GHI、DHI、DNI的POA计算问题求助（pvlib库）

我用pvlib通过GHI、DNI、DHI、太阳方位角、天顶角计算POA辐照度，晴天结果正常，但阴天时波动极大，求解决思路。

我的代码如下：

i=24
for j in range(i):
    Y = X[j]
    M = ['ghi','dni','dhi','zenith','azimuth','period_end','air_temp']
    values = list( map(Y.get, M) )
    df_poa = pvlib.irradiance.get_total_irradiance(
        surface_tilt=10,  # 与水平面夹角10度
        surface_azimuth=0,  # 朝南                                            
        dni=values[1],
        ghi=values[0],
        dhi=values[2],
        solar_zenith=values[3],
        solar_azimuth=values[4],
        model='isotropic')
    OUT = dict(df_poa)
    Final_out = list(OUT.values())
    POA_Global = Final_out[0]

可能的原因及解决方法 #

• 低太阳高度角的数值敏感性：阴天时太阳天顶角通常接近90°，此时cos(zenith)的数值对天顶角的微小变化极其敏感，一点点测量误差就会导致POA结果大幅波动。
  解决：对天顶角大于85°的数据直接标记为无有效辐照度（设POA为0），或者用平滑算法处理输入的天顶角数据。

• DNI数据可靠性问题：阴天时DNI几乎为0，但实际测量数据可能带有噪声（比如微小的正负波动），isotropic模型对DNI的微小变化会放大输出波动。
  解决：过滤DNI小于阈值（比如5 W/㎡）的数据，将其设为0；或者更换散射模型，比如haydavies或reindl，这类模型对低DNI场景的鲁棒性更强。

• 太阳方位角定义不匹配：pvlib中太阳方位角的定义是从南向顺时针度量，如果你的输入数据是从北向逆时针度量的标准方位角，会导致方位角差值的余弦计算异常，尤其在低太阳高度时影响更明显。
  解决：统一方位角定义，将输入的方位角转换为pvlib兼容的格式：pvlib_azimuth = (azimuth + 180) % 360（如果输入是北向为0的方位角）。

• 输入数据噪声：阴天时辐照强度本身很低，传感器的测量噪声占比更高，会直接导致POA计算结果波动。
  解决：对GHI、DHI做滑动平均处理（比如取前后3个时间点的均值），减少瞬时噪声的影响。

代码优化建议 #

可以简化取值逻辑，增加异常处理，示例如下：

i=24
# 定义阈值
ZENITH_THRESHOLD = 85
DNI_THRESHOLD = 5

for j in range(i):
    Y = X[j]
    # 直接取值更清晰
    ghi = Y['ghi']
    dni = max(Y['dni'], 0)  # 确保DNI非负
    dhi = Y['dhi']
    zenith = Y['zenith']
    azimuth = Y['azimuth']
    
    # 处理低太阳高度情况
    if zenith >= ZENITH_THRESHOLD:
        POA_Global = 0
        continue
    
    # 过滤无效DNI
    if dni < DNI_THRESHOLD:
        dni = 0
    
    # 更换更鲁棒的模型
    df_poa = pvlib.irradiance.get_total_irradiance(
        surface_tilt=10,
        surface_azimuth=0,
        dni=dni,
        ghi=ghi,
        dhi=dhi,
        solar_zenith=zenith,
        solar_azimuth=azimuth,
        model='haydavies')  # 替代isotropic
    
    POA_Global = df_poa['poa_global']  # 直接取列名更可靠

内容的提问来源于stack exchange，提问作者chandrakant singh