GMS:Risk Analysis Wizard
From XMS Wiki - http://wikis.aquaveo.com/xms/
|Building a MODFLOW Model|
|Map to MODFLOW|
|Packages Supported in GMS|
|Saving a MODFLOW Simulation|
|Importing MODFLOW Files|
|Unsupported MODFLOW Features|
|MODFLOW Display Options|
|MODFLOW Post Processing Viewing Options|
|Reading a MODFLOW Simulation|
|Flow:||BCF6, HUF, LPF, UPW|
|Solvers:||PCGN, SAMG, SIP, SOR|
|Automated Parameter Estimation|
|Parameter Estimation Dialog|
|Multiplier Arrays for Parameters|
|Standard MODFLOW Parameters|
|Gaussian Field Generator|
|Risk Analysis Wizard|
The risk analysis wizard is a tool for refining stochastic modeling results. The Risk Analysis Wizard has two options. You can either perform a Capture Zone Analysis or a Threshold Analysis. Capture zone analysis requires a MODFLOW solution set, but threshold analysis can be performed on any solution set type.
Capture Zone Analysis
Capture zone analysis includes running MODPATH for each of the MODFLOW solutions to generate a capture zone for each well or zone code group in the MODFLOW model. These capture zones are combined into one probabilistic capture zone
Threshold analysis is similar to capture zone analysis, but threshold analysis can be performed using any type of 3D data set that has been read into the Project Explorer. You first select a simulation set from the Project Explorer. Next, you set up rules for generating a probabilistic threshold dataset. For example, you might be looking at a specific contaminant and want to know what the probability is that the concentration of this contaminant will be above the EPA level for drinking water. You would generate a rule reflecting this limit. After processing all rules for each simulation, GMS creates a probabilistic threshold dataset.
Risk Analysis Wizard Dialogs
The risk analysis wizard is entered through the right-click menu for a folder in the Project Explorer and choose Risk Analysis... .
Choosing Between Capture Zone and Threshold Analysis
The first step in the risk analysis wizard is to choose what type of simulations you want to process.
The window will list all the solution types as the analysis wizard will only process solutions of the same type. Probabilistic capture zone analysis can only be performed on MODFLOW solutions. Probabilistic threshold analysis can be performed on any solution type.
Capture Zone Analysis Dialog
Capture zone analysis can only be performed on MODFLOW solutions
Analysis Options – Individual wells
Choose this option to create a probabilistic capture zone for each uniquely named well in the MODFLOW model.
Analysis Options – Well groups
Choose this option to create a probabilistic capture zone for each different zone code number.
Particle Starting Locations – Distribute particles on water table surface
Choose this option to distribute particles only on the water table surface. By default, one particle is placed at the xy center of the cell at the water table, but you can use the Particle placement within cells option to change the number of particles placed on the water table in each cell. This option will create a 2D probabilistic capture zone representing the intersection of the entire 3D probabilistic capture zone with the water table.
Particle Starting Locations - Distribute particles within cells
Choose this option to distribute particles within each cell. By default, one particle is placed at the center of each cell, but you can use the Particle placement within cells option to change the number of particles placed on within each cell. This option creates three different probabilistic capture zones. The first represents the 3D capture zone, the other two represent different 2D projections of the 3D capture zone.
Tracking Duration - To end
This option sets MODPATH to move particles through the flow field until they exit the model.
Tracking Duration - Specified duration
This option sets MODPATH to move particles through the flow field until they either exit the model or the time reaches the duration set, which ever comes first.
Particle Placement Within Cells
This option allows you to change the number of particles per cell from the default of one. Increase the number of particles leads to a smoother capture zone, but drastically increases the computation time.
Weight Results Based On Residual Error
In all of the capture zone methods, the algorithms used to synthesize the probability data set can be weighted using observation data. This makes it possible to give more weight to model instances with smaller calibration error when calculating the capture zone probabilities. The weighted head and flow observations can be compared to the computed values to come up with a global error norm, E, for each model run. This error norm can be based on the root mean squared (RMS) error, the sum of the squared weighted residuals, or any other measure selected by the modeler. For capture zone analysis, GMS uses the sum of squared weighted residuals (SSWR) for the error norm. The error norm from each MODFLOW run is used to compute a weight for the given solution using the following equation:
Where Wi is the weight applied to solution i, is a user-defined factor, ME is the mean of the error values from all solutions, Ejk is the error for solution i, and SD is the standard deviation of error values from all solutions. The weights are also normalized as follows
so that the weights sum to unity. Equation 1 was developed to give the greater emphasis to the lower error values and to allow the user to control the relative emphasis given to low vs. high values simply by adjusting the value. The equation also avoids problems when one of the error values is zero, since a zero error value does not result in an infinite weight. We also wanted the equation to scale the weights according to the data being examined. This is done by subtracting the individual SSWR from the mean error and dividing by the standard deviation.
Equation 1 centers the weights on the mean error. The relative weight given to values differing from the mean is biased by the factor. This makes it possible to bias the resulting weight using knowledge of the site and the quality of the observation data.
The figure below shows how the factor in Equation 1 affects the weight applied to a given error. An factor of 1.2 makes the contribution of each SSWR almost linear, whereas an factor of ten gives most of the weight to the lowest 5-10 percent while discounting the other error values. We typically use an value of 2.0.
Threshold Analysis Dialog
Threshold analysis can be performed on any collection of 3D datasets. The threshold analysis dataset is created by using rules. For each rule, you select a dataset (only applies to solutions with multiple data sets.), the greater or less than sign and a value. You can have as many rules as you want. The AND and OR logic options allow you to use: 1st rule AND 2nd rule AND 3rd rule; or 1st rule OR 2nd rule OR 3rd rule. Threshold analysis creates one 3D dataset.
GMS – Groundwater Modeling System
|Modules:||2D Grid · 2D Mesh · 2D Scatter Point · 3D Grid · 3D Mesh · 3D Scatter Point · Boreholes · GIS · Map · Solid · TINs|
|Models:||ADH · ART3D · FEMWATER · MODAEM · MODFLOW · MODPATH · MT3DMS · PEST · RT3D · SEAM3D · PHT3D · SEAWAT · SEEP2D · T-PROGS · UTCHEM · UTEXAS · ZONEBUDGET|