FilesExpand file tree

 eta

/

plot_fig3A_hsu_soiling_example.py

Copy path

Blame

More file actions

Blame

More file actions

Latest commit

 

History

History

History

79 lines (70 loc) · 2.6 KB

 eta

/

plot_fig3A_hsu_soiling_example.py

Copy path

Top

File metadata and controls

• Code
• Blame

79 lines (70 loc) · 2.6 KB

Raw

Copy raw file

Download raw file

Open symbols panel

Edit and raw actions

1

2

3

4

5

6

7

8

9

10

11

12

13

14

15

16

17

18

19

20

21

22

23

24

25

26

27

28

29

30

31

32

33

34

35

36

37

38

39

40

41

42

43

44

45

46

47

48

49

50

51

52

53

54

55

56

57

58

59

60

61

62

63

64

65

66

67

68

69

70

71

72

73

74

75

76

77

78

79

"""

HSU Soiling Model Example

=========================

Example of soiling using the HSU model.

"""

# %%

# This example shows basic usage of pvlib's HSU Soiling model [1]_ with

# :py:func:`pvlib.soiling.hsu`.

#

# References

# -----------

# .. [1] M. Coello and L. Boyle, "Simple Model For Predicting Time Series

# Soiling of Photovoltaic Panels," in IEEE Journal of Photovoltaics.

# doi: 10.1109/JPHOTOV.2019.2919628

#

# This example recreates figure 3A in [1]_ for the Fixed Settling

# Velocity case.

# Rainfall data comes from Imperial County, CA TMY3 file

# PM2.5 and PM10 data come from the EPA. First, let's read in the

# weather data and run the HSU soiling model:

import pathlib

from matplotlib import pyplot as plt

from pvlib import soiling

import pvlib

import pandas as pd

# get full path to the data directory

DATA_DIR = pathlib.Path(pvlib.__file__).parent / 'data'

# read rainfall, PM2.5, and PM10 data from file

imperial_county = pd.read_csv(DATA_DIR / 'soiling_hsu_example_inputs.csv',

index_col=0, parse_dates=True)

rainfall = imperial_county['rain']

depo_veloc = {'2_5': 0.0009, '10': 0.004} # default values from [1] (m/s)

rain_accum_period = pd.Timedelta('1h') # default

cleaning_threshold = 0.5

tilt = 30

pm2_5 = imperial_county['PM2_5'].values

pm10 = imperial_county['PM10'].values

# run the hsu soiling model

soiling_ratio = soiling.hsu(rainfall, cleaning_threshold, tilt, pm2_5, pm10,

depo_veloc=depo_veloc,

rain_accum_period=rain_accum_period)

# %%

# And now we'll plot the modeled daily soiling ratios and compare

# with Coello and Boyle Fig 3A:

daily_soiling_ratio = soiling_ratio.resample('d').mean()

fig, ax1 = plt.subplots(figsize=(8, 2))

ax1.plot(daily_soiling_ratio.index, daily_soiling_ratio, marker='.',

c='r', label='hsu function output')

ax1.set_ylabel('Daily Soiling Ratio')

ax1.set_ylim(0.79, 1.01)

ax1.set_title('Imperial County TMY')

ax1.legend(loc='center left')

daily_rain = rainfall.resample('d').sum()

ax2 = ax1.twinx()

ax2.plot(daily_rain.index, daily_rain, marker='.',

c='c', label='daily rainfall')

ax2.set_ylabel('Daily Rain (mm)')

ax2.set_ylim(-10, 210)

ax2.legend(loc='center right')

fig.tight_layout()

fig.show()

# %%

# Here is the original figure from [1]_ for comparison:

#

# .. image:: ../_images/Coello_Boyle_2019_Fig3.png

# :alt: Figure 3A from the paper showing a simulated soiling signal.

#

# Note that this figure shows additional timeseries not calculated here:

# modeled soiling ratio using the 2015 PRISM rainfall dataset (orange)

# and measured soiling ratio (dashed green).