2

我正在为我的硕士论文使用 pvlib 库。当我在下午 4 点以后运行该示例时,它通常会报告交流电源 -0.02。有人知道为什么吗?我正在使用下面的代码:

import pandas as pd
import numpy as np

# pvlib imports
import pvlib

from pvlib.pvsystem import PVSystem
from pvlib.location import Location
from pvlib.modelchain import ModelChain

# load some module and inverter specifications
sandia_modules = pvlib.pvsystem.retrieve_sam('SandiaMod')
cec_inverters = pvlib.pvsystem.retrieve_sam('cecinverter')

sandia_module = sandia_modules['Canadian_Solar_CS5P_220M___2009_']
cec_inverter = cec_inverters['ABB__MICRO_0_25_I_OUTD_US_208_208V__CEC_2014_']
location = Location(latitude=49.0205559, longitude=12.057453900000041)
system = PVSystem(surface_tilt=20, surface_azimuth=200,
              module_parameters=sandia_module,
              inverter_parameters=cec_inverter)
mc = ModelChain(system, location)
python_native_dt = datetime.datetime.now()
weather = pd.DataFrame([[1050, 1000, 100, 30, 5]],
                   columns=['ghi', 'dni', 'dhi', 'temp_air', 'wind_speed'],
                   index=[pd.Timestamp(pytz.timezone('Etc/GMT+2').localize(python_native_dt))])

mc.run_model(times=weather.index, weather=weather)
print(mc.ac)

执行mc.ac将导致:2018-06-05 16:20:19.117017-02:00 -0.02 dtype: float64

4

2 回答 2

1

-0.02是您选择的逆变器在输入直流功率低于其激活阈值时消耗的能量。

为了提高可重复性并帮助我们找到答案,我建议您指定一个确切的时间,而不是依赖datetime.datetime.now(). 使用index=[pd.Timestamp('2018-06-05 16:20:19.117017-02:00')],我得到2018-06-05 16:20:19-02:00 13.660678

我建议您确认mc.aoimc.solar_position与您的天气输入一致。它们源自时间索引,用于计算阵列辐照度平面。

如果那没有帮助... pvlib 和 pandas 的版本是什么?请注意,该示例还需要import pytzandimport datetime才能运行。

于 2018-06-07T16:59:44.810 回答
0

G'day 马克西米利安,

我怀疑这是天气数据框。我注意到来自 NWP 中心的大多数天气文件(我正在使用 ECMWF ERA5)代表 DNI GHI DHI 的累积。对于 ERA5 数据,我注意到原始数据的负功率值。

通过除以 60(secs)*60(minutes) 从累积转换为 W/s 后,它起作用了。

于 2018-07-05T23:19:40.957 回答