WMS:Stochastic Modeling

There is always a great deal of uncertainty in hydrologic and hydraulic modeling and the parameters that are used to develop solutions. Despite this, a typical flood plain boundary is black and white in that you are either in or out of the flood plain. A good engineer might be able to dispute a flood plain boundary by performing a hydrologic/hydraulic analysis with a set of equally probable parameters that results in a difference in the flood plain delineation. Until recently, computer programs lacked the ability to consider multiple probable answers and report a probabilistic floodplain boundary, but with the Stochastic Modeling tools in WMS this is possible using a combination of HEC-1 for hydrologic analysis, HEC-RAS for 1D hydraulic river modeling, and the WMS flood plain delineation tools.

You can "connect" the results of HEC-1 to a developed HEC-RAS model and then run them as many times consecutively, with the results of the HEC-1 analysis feeding the boundary conditions for an HEC-RAS model. Certain parameters (at this point only basin CN and precipitation within HEC-1, and Manning's roughness within HEC-RAS) can be varied within a range of reasonable values using Monte Carlo or Latin Hypercube simulations in order to create a number of simulations. The results of each HEC-RAS model can then be used to delineate a series of flood plains. The combination of all floodplains can then be examined in order to derive a "probabilistic" flood plain where a region flooded by 100% of the model simulation combinations can be distinguished from an area that is flooded by only 50% of the models as shown in the figure below:

File:WMSImage277.jpg

The following steps outline the process for developing a stochastic model of floodplain boundaries using WMS, HEC-1, and HEC-RAS.

Step 1. Develop a hydrologic model with HEC-1

File:Image85.jpg

The HEC-1 Interface tools in WMS can be used to develop a working HEC-1 model. It is important that this model be running effectively (i.e. calibrated and/or adjusted to give credible results) prior to using it for the Stochastic Modeling simulation.

The only parameters at this point that can be varied within a range of probable answers are rainfall and CN (curve number for the different basins. This is done by setting the precipitation or CN to be a negative number in their respective dialogs. This negative number is a key number and should be unique for each stochastic variable you create. When running the stochastic simulation WMS will substitute the simulation specific parameter for the defined key. You can then setup a stochastic variable for HEC-1 in the Stochastic Run Parameters dialog. A key value (matching the key you defined in the materials property) starting value, min value, max value, standard deviation and distribution type.

File:WMSImage310.jpg

Step 2. Develop a working model in HEC-RAS

File:Image84.jpg

The HEC-RAS interface tools in WMS can be used to developing a working HEC-RAS model. It is important that this model be running effectively (i.e. calibrated and/or adjusted to give credible results) prior to using it for the Stochastic Modeling simulation.

File:WMSImage311.jpg

The only parameter at this point that can be varied within a range of probable answers are Manning's coefficients for the different material types. This is done by setting the roughness to be a negative number in the HEC-RAS Materials dialog. This negative number is a key number and should be unique for each stochastic variable you create. When running the stochastic simulation WMS will substitute the simulation specific parameter for the defined key. You can then setup a stochastic variable for HEC-RAS in the Stochastic Run Parameters dialog. A key value (matching the key you defined in the materials property) starting value, minimum value, maximum value, standard deviation and distribution type.

Step 3. Establish appropriate flood plain delineation parameters for your area

The floodplain delineation portion of the stochastic modeling uses the results from each HEC-RAS model to develop a floodplain for each run. The floodplain delineation is the same model used by WMS to perform individual floodplain delineations.

The only parameters at this point that can be varied within a range of probable answers the search radius of the flood plain delineation. This is done by setting the search radius in the Floodplain Delineation Options be a negative number. This negative number is a key number and should be unique for each stochastic variable you create. When running the stochastic simulation WMS will substitute the simulation specific parameter for the defined key. You can then setup a stochastic variable for Floodplain in the Stochastic Run Parameters dialog. A key value (matching the key you defined in the materials property) starting value, min value, max value, standard deviation and distribution type.

Step 4. Set up stochastic simulation parameters

There are several simulation parameters that control a stochastic simulation. They are defined in the Stochastic Run Parameters dialog shown below. Each section of this dialog is discussed below.

File:Image86.jpg

Simulation Type

The simulation type and number of simulations can be set. In a Monte Carlo simulation, each specified input variable is randomly varied within a specified minimum and maximum value a given number of times. If only a few simulations are run it may not be guaranteed to fully explore the parameter space. A Latin Hypercube simulation, on the other hand, divides the range into intervals and insures that parameters are chosen from each interval. With this kind of simulation you are more likely to explore the parameter space with fewer simulations.

Stochastic Models

HEC-1 or TR-20 for hydrologic modeling, HEC-RAS for hydraulic modeling, and Floodplain Delineation are the only currently available models for stochastic modeling. For each model you include a basin input file and solution files directory needs to be defined. For HEC-1, TR-20, or HEC-RAS you should select the input file of the already created model. These models will have key values (negative numbers) for the input parameters that will be defined as stochastic variables. The current floodplain delineation options will be saved in the flood run file.

You can add stochastic variables for any of the models. Each stochastic variable requires a key value (a negative number that has been entered in place of a parameter such as precipitation), a type, a starting value, a minimum value, a maximum value, a standard deviation, and a distribution. The distribution can be either normal or uniform and optionally defined as log.

Step 5. Assign boundary conditions between models and run the stochastic simulation

The final step before running a stochastic simulation is to link any required models. In particular you have to link how the hydrologic modeling results are used as boundary conditions for the HEC-RAS hydraulic model. When you choose the Run Stochastic Model option the Model Steering Dialog appears. For the hydrologic model you are using (HEC-1) you must assign the appropriate hydrograph (basin or outlet) to the river reach and section in the hydraulic model.

The option to Assign Hydrographs Automatically can be used providing a drainage coverage and centerline coverage are being used to set up the hydrologic and hydraulic models.

After identifying the hydrograph and the river, reach, and cross section station select the Assign Hydrograph BC button to link the models for this point. Continue until all of the appropriate model locations are linked.

When selecting OK your model simulation will run.

File:Image144.jpg

Step 6. Run the simulation

File:Image83.jpg

Once the simulation begins running you will see the parameters selected for each model run, as well as the status of each run. The Read solution on exit toggle is on by default and will cause that all model solutions (hydrographs, water surface elevations, and floodplain delineations) are read when the simulations are completed.

Step 7. Post Process the results

After finishing a stochastic simulation there are two primary results read back into WMS for each simulation: hydrographs from the HEC-1 model, and the floodplain depths and water surface elevations for each run.

Hydrographs

A series of hydrographs are loaded for each hydrograph station and can be viewed in the normal way hydrographs are viewed.

File:Image91.jpg

Floodplains

Each floodplain delineation results in a water surface elevation and a flood depth dataset. Each pair of datasets are organized in a folder underneath the TIN in the Project Explorer. You can set the contour options for a TIN and select the dataset you wish to be active and displayed from the Project Explorer.

File:Image330.gif

File:Image137.jpg

Probabilistic Floodplain Boundary

A final dataset showing the probability of flooding is also created from all of the individual floodplain datasets. The probabilistic flood plain indicates for each vertex on the TIN the percentage of model runs that resulting in inundation at the point.

File:Image135.jpg

File:Image136.jpg Close up view

Annual Exceedance Probability (AEP) Map

When using the stochastic hydrologic/hydraulic/floodplain delineation tools it is possible to generate an annual exceedance probability map. This is done by generating inundation maps that consider the range of all possible floods for all return periods. The result is a map that identifies the annual exceedance probability of flooding for every TIN vertex. For example if a point is flooded 10 times in 1000 simulations then it would represent the .01 probability. The Return Period → Feature Objects... command allows you to generate contours from the AEP map for specified return periods. You should remember that such a map does not represent a solution from a single set of input parameters, but is rather the composite of several hundred or thousand simulations.

File:WMSImage312.jpg

Related Topics

HEC-1 Modeling

HEC-RAS Modeling

Flood Plain Delineation