Description Usage Arguments Value Examples

Simulates a reservoir for a given inflow time series and assuming Standard Operating Policy (meet target at all times, unless constrained by available water in reservoir plus incoming flows) or an optimised policy deived using `sdp_supply`

.

1 2 |

`Q` |
vector or time series object. Net inflow totals to the reservoir. Mm^3 (Million cubic meters). |

`target` |
numerical constant, or a time series or vector of the target releases. Must be the same length as Q is given as a vector or time series. Mm^3 (Million cubic meters). |

`capacity` |
numerical. The reservoir capacity. Should be same volumetric unit as Q. Mm^3 (Million cubic meters). |

`surface_area` |
numerical. The reservoir surface area at full capacity. Must be in square kilometers (km^2), or Mm^2. |

`max_depth` |
numerical. The maximum water depth of the reservoir at maximum capacity. Must be in meters. If omitted, the depth-storage-area relationship will be estimated from surface area and capacity only. |

`evap` |
vector or time series object of length Q, or a numerical constant. Evaporation from losses from reservoir surface. Varies with level if depth and surface_area parameters are specified. Recommended units: meters, or kg/m2 * 10 ^ -3. |

`double_cycle` |
logical. If TRUE the Q and R time series will be replicated and placed end-to-end to double the simulation. Recommended if the critical period occurs at the end of the sequence. |

`plot` |
logical. If TRUE (the default) the storage and release time series are plotted. |

`S_initial` |
numerical. The initial storage as a ratio of capacity (0 <= S_initial <= 1). The default value is 1. |

`policy` |
list. The output of the SDP function. If omitted, Standard Operating Policy is assumed. |

Returns the no-fail storage capacity and corresponding storage behaviour time series.

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 | ```
# simulate a reservoir assuming standard operating policy, then compare with SDP-derived policy
#trained on historical flows.
# DEFINE RESERVOIR SPECS AND MODEL INPUTS
res_cap <- 1500 #Mm3
targ <- 150 #Mm3
area <- 40 #km2
max_d <- 40 #m
ev = 0.2 #m
Q_pre1980 <- window(resX$Q_Mm3, end = c(1979, 12), frequency = 12)
Q_post1980 <- window(resX$Q_Mm3, start = c(1980, 1), frequency = 12)
# SIMULATE WITH SOP
layout(1:3)
simSOP <- simRes(Q_post1980, capacity = res_cap, target = targ,
surface_area = area, max_depth = max_d, evap = ev)
# TRAIN SDP POLICY ON HISTORICAL FLOWS
policy_x <- sdp_supply(Q_pre1980, capacity = res_cap, target = targ,
surface_area = area, max_depth = max_d, evap = ev, Markov = TRUE, plot = FALSE, S_disc = 100)
# SIMULATE WITH SDP-DERIVED POLICY
simSDP <- simRes(Q_post1980, capacity = res_cap, target = targ,
surface_area = area, max_depth = max_d, evap = ev, policy = policy_x)
``` |

swd-turner/reservoir documentation built on Sept. 7, 2017, 8:11 p.m.

Embedding an R snippet on your website

Add the following code to your website.

For more information on customizing the embed code, read Embedding Snippets.