Download User Guide - WRc Websites

Transcript
STORMPAC 4.1 User Guide
C4
DISAGGREGATION PROCEDURE
C4.1
Summary of the Algorithm
The disaggregation procedure can be summarised as follows. A 1-hour rainfall depth is read in
from a file, which contains the hourly series to be disaggregated. A 5-minute series is simulated
using the fitted model. The 5-minute simulated series is summed and the 1-hour total of the
simulated series compared to the 1-hour total that was read in. The simulated series is discarded
if the absolute difference between the simulated 1-hour total and the total read in exceeds
0.05 mm. The process is repeated until the totals are in agreement to within 0.05 mm. The
simulated series then represents a possible realisation of 5-minute data representative of the 1hour value that was read in. Following this procedure, a record of 1-hour rainfall depths can be
disaggregated.
In the above procedure, some additional rules are applied. These include using overlapping 5minute values from the previous disaggregated hour provided the total of these overlapping
values do not exceed the 1-hour total to be disaggregated. This allows for some influence, due to
overlapping cells, from the previous hour. In addition, the disaggregated series are scaled to
achieve an exact match to the hourly series (although clearly this is only a very minor adjustment,
given they fall to within 0.05 mm of the hourly total).
It should be noted that the disaggregation procedure is not accounted for when fitting the model –
it is a separate procedure that could be applied to any model capable of simulating 5-minute
series. Therefore it is likely that some (probably small) bias would be introduced into the
distributional properties of the resultant disaggregated series.
C4.2
Tests
The Farnborough data were aggregated to hourly values and then disaggregated using the
above procedure. The distribution of the resultant disaggregated series was compared to the
historical 5-minute distribution using quantile plots; the results are shown in Figures C3-5.
Figure C3 gives the overall result, where the data for all months have been pooled. From this
figure it can be seen that the distribution of the disaggregated series compares favourably to the
distribution of the historical 5-minute series. The exception is a slight curvature in the mid to lower
part of the plot, which would translate into a slight over-estimation in part of the distribution tail.
The plots for January and July show a similar pattern, i.e. a slight over-estimation for some of the
distribution tail (Figures C4 and C5).
89
© WRc plc 2009