OBF-file: Difference between revisions
(12 intermediate revisions by the same user not shown) | |||
Line 2: | Line 2: | ||
<div class="info_versionen"> | <div class="info_versionen"> | ||
'''Versions:''' | '''Versions:''' | ||
* '''current''' v1.9 | * '''current''' v1.9 | ||
* [{{SERVER}}/index.php?title=OBF-file&oldid=7469 v1.8] (up to r1843) | * [{{SERVER}}/index.php?title=OBF-file&oldid=7469 v1.8] (up to r1843) | ||
* [{{SERVER}}/index.php?title=ZIE-file&oldid=7401 v1.7] (as of r1785) | * [{{SERVER}}/index.php?title=ZIE-file&oldid=7401 v1.7] (as of r1785) | ||
Line 19: | Line 19: | ||
The file is divided into several blocks, each of which is optional. Each block is identified by its heading: | The file is divided into several blocks, each of which is optional. Each block is identified by its heading: | ||
* '''Series''': This block defines objective functions that are based on a comparison of simulation time series with a reference time series. | * '''Series''': This block defines objective functions that are based on a comparison of simulation time series with a reference time series. | ||
* '''Values''': This block defines objective functions that are based on a comparison of simulation result values with a reference value ( | * '''Values''': This block defines objective functions that are based on a comparison of simulation result values with a reference value (only for [[SWMM]] applications). | ||
* '''ValueFromSeries''': This block defines objective functions that are based on a comparison of values that are computed from simulated time series with reference values. | * '''ValueFromSeries''': This block defines objective functions that are based on a comparison of values that are computed from simulated time series with reference values. | ||
* '''Aggregate''': Aggregate objective functions | * '''Aggregate''': Aggregate objective functions | ||
Line 32: | Line 31: | ||
*------ | *------ | ||
* | * | ||
*|-----|---------------|---------|---|-------|-------|------------|---------|------------------|------------------|-----------------------------|----------| | *|-----|------------------|---------|---|-------|-------|------------|---------|------------------|------------------|-----------------------------|----------| | ||
*| Opt | Name | *| Opt | Name | Group | D | OpFac | File | SimResult | ObjFunc | Evaluation period | Reference series | CurValue | | ||
*| | | *| | | | | | | | | Start | End | RefName | RefFile | | | ||
*|-----|---------------|---------|---|-------|-------|------------|---------|------------------|------------------|-----------|-----------------|----------| | *|-----|------------------|---------|---|-------|-------|------------|---------|------------------|------------------|-----------|-----------------|----------| | ||
| P | NashSutcliffe | | + | | WEL | B401_1AB | NSE | 30.05.1984 00:00 | 18.06.1984 23:30 | | seriesfile.zre | | | P | Kling-Gupta | | + | | WEL | B401_1AB | KGE | 30.05.1984 00:00 | 18.06.1984 23:30 | | seriesfile.zre | | | ||
| | | S | NashSutcliffe | | + | | WEL | B401_1AB | NSE | 30.05.1984 00:00 | 18.06.1984 23:30 | | seriesfile.zre | | | ||
*|-----|---------------|---------|---|-------|-------|------------|---------|------------------|------------------|-----------|-----------------|----------| | | S | VolumeError | | - | | WEL | B401_1AB | BIAS | 30.05.1984 00:00 | 18.06.1984 23:30 | | seriesfile.zre | | | ||
| S | MeanSquaredError | | - | | WEL | B401_1AB | MSE | 30.05.1984 00:00 | 18.06.1984 23:30 | | seriesfile.zre | | | |||
*|-----|------------------|---------|---|-------|-------|------------|---------|------------------|------------------|-----------|-----------------|----------| | |||
* | * | ||
*Values | *Values | ||
Line 47: | Line 48: | ||
*| | | | | | | | | (Block) | (Column) | RefValue | | | *| | | | | | | | | (Block) | (Column) | RefValue | | | ||
*|-----|---------------|---------|---|-------|-------|-----------|---------|---------------|-------------|----------|----------| | *|-----|---------------|---------|---|-------|-------|-----------|---------|---------------|-------------|----------|----------| | ||
| S | OutfallLoad | | - | | RPT | | | | S | OutfallLoad | | - | | RPT | | AE | OutfallNode | AvgFlow | 0 | | | ||
| S | Storage401 | Storage | - | + | BLZ | B401 | | Speicher | maxVolume | | | | | S | Storage401 | Storage | - | + | BLZ | B401 | | Speicher | maxVolume | | | | ||
| S | Storage701 | Storage | - | + | BLZ | B701 | | Speicher | maxVolume | | | | | S | Storage701 | Storage | - | + | BLZ | B701 | | Speicher | maxVolume | | | | ||
Line 59: | Line 60: | ||
*| | | | | | | | | Start | Ende | ValueType | RefValue | | | *| | | | | | | | | Start | Ende | ValueType | RefValue | | | ||
*|-----|---------------|---------|---|-------|-------|------------|---------|---------|---------|-----------|----------|----------| | *|-----|---------------|---------|---|-------|-------|------------|---------|---------|---------|-----------|----------|----------| | ||
| S | OutfallLoad | | - | | WEL | B401_1AB | | | S | OutfallLoad | | - | | WEL | B401_1AB | AE | | | | | | | ||
| | | P | PeakError | | - | | WEL | S3 _1ZU | SE | | | Max | 29.33 | | | ||
*|-----|---------------|---------|---|-------|-------|------------|---------|---------|---------|-----------|----------|----------| | *|-----|---------------|---------|---|-------|-------|------------|---------|---------|---------|-----------|----------|----------| | ||
* | * | ||
*Aggregate | *Aggregate | ||
Line 104: | Line 95: | ||
Direction of desired optimization: | Direction of desired optimization: | ||
*<code>-</code>: The objective function should be minimized (default) | *<code>-</code>: The objective function should be minimized (default) | ||
*<code>+</code>: The objective function should be maximized | *<code>+</code>: The objective function should be maximized (this causes the objective function values to be multiplied by <code>-1</code> internally). | ||
|- | |- | ||
| '''Op'''erator or '''Fac'''tor<br/>(optional) | | '''Op'''erator or '''Fac'''tor<br/>(optional) | ||
Line 113: | Line 104: | ||
*''number'': The objective function value is multiplied by the given number and then added to the group | *''number'': The objective function value is multiplied by the given number and then added to the group | ||
|- | |- | ||
| '''File''' || The simulation result file containing the simulation result from which the objective function should be calculated. Supported values: | | '''File''' || The file extension of the simulation result file containing the simulation result from which the objective function should be calculated. Supported values: | ||
* <code>WEL</code> (for [[BlueM.Sim]] or [[TALSIM]]) | * <code>WEL</code> (for [[BlueM.Sim]] or [[TALSIM]]) | ||
* <code>KWL</code> (only for [[BlueM.Sim]]) | * <code>KWL</code> (only for [[BlueM.Sim]]) | ||
* <code> | * <code>WBL</code> (only for [[TALSIM]]) | ||
* <code>TEMP.WEL</code> (only for [[TALSIM]]) | * <code>KTR.WEL</code>, <code>TEMP.WEL</code>, <code>CHLO.WEL</code> or any other result file with the ending <code>.WEL</code> (only for [[TALSIM]]) | ||
* <code>ASC</code> (only for [[SMUSI]] and [[SWMM]]) | * <code>ASC</code> (only for [[SMUSI]] and [[SWMM]]) | ||
|- | |- | ||
Line 127: | Line 118: | ||
* <code>AE</code>: Absolute error | * <code>AE</code>: Absolute error | ||
Supported values for comparing series: | Supported values for comparing series: | ||
* <code>MAE</code>: Mean abolute error | |||
* <code>SAE</code>: Sum of absolute errors | |||
* <code>MSE</code>: Mean squared error | |||
* <code>SSE</code>: Sum of squared errors | * <code>SSE</code>: Sum of squared errors | ||
* <code>NSE</code>: Nash-Sutcliffe efficiency{{:Literatur:Nash-Sutcliffe_1970|}} | * <code>NSE</code>: Nash-Sutcliffe efficiency{{:Literatur:Nash-Sutcliffe_1970|}} | ||
:<code>-∞ ≤ NSE ≤ 1</code>. Larger is better. | :<code>-∞ ≤ NSE ≤ 1</code>. Larger is better. | ||
Line 134: | Line 127: | ||
* <code>lnNSE</code>: Logarithmic Nash-Sutcliffe efficiency{{:Literatur:Nash-Sutcliffe_1970|}} | * <code>lnNSE</code>: Logarithmic Nash-Sutcliffe efficiency{{:Literatur:Nash-Sutcliffe_1970|}} | ||
:<code>-∞ ≤ lnNSE ≤ 1</code>. Larger is better. | :<code>-∞ ≤ lnNSE ≤ 1</code>. Larger is better. | ||
:<math>lnNSE=1-\frac{\sum_{t=1}^T\left(ln(Q_{obs}^t)-ln(Q_{sim}^t)\right)^2}{\sum_{t=1}^T\left(ln(Q_{obs}^t)- | :<math>lnNSE=1-\frac{\sum_{t=1}^T\left(ln(Q_{obs}^t+\epsilon)-ln(Q_{sim}^t+\epsilon)\right)^2}{\sum_{t=1}^T\left(ln(Q_{obs}^t+\epsilon)-\overline{ln(Q_{obs}+\epsilon)}\right)^2}</math> | ||
<code>ε</code> is set to 1% of the average value of <code>Q<sub>obs</sub></code> as recommended by {{:Literatur:Pushpalatha_2012}}. | |||
* <code>KGE</code>: Kling-Gupta efficiency{{:Literatur: | * <code>KGE</code>: Kling-Gupta efficiency{{:Literatur:Kling_2012|}} | ||
:<code>-∞ ≤ KGE ≤ 1</code>. Larger is better. | :<code>-∞ ≤ KGE ≤ 1</code>. Larger is better. | ||
:<math>\text{KGE} = 1 - \sqrt{ (r - 1)^2 + (\beta - 1)^2 + (\gamma - 1)^2 }</math> | :<math>\text{KGE} = 1 - \sqrt{ (r - 1)^2 + (\beta - 1)^2 + (\gamma - 1)^2 }</math> | ||
Line 151: | Line 144: | ||
* <code>BIAS</code>: Absolute volume error [%] | * <code>BIAS</code>: Absolute volume error [%] | ||
The following options are '''deprecated''' and only kept for backwards compatibility: | The following options are '''deprecated''' and only kept for backwards compatibility: | ||
* <code>NashSutt</code>: Modified Nash-Sutcliffe efficiency | * <code>NashSutt</code>: Modified Nash-Sutcliffe efficiency | ||
:modified formula with <code>0 ≤ NSE<sub>mod</sub> ≤ ∞</code>. Smaller is better. | :modified formula with <code>0 ≤ NSE<sub>mod</sub> ≤ ∞</code>. Smaller is better. | ||
Line 178: | Line 163: | ||
If '''Start''' and/or '''End''' are left empty, the start and/or end of the entire simulation time period is used. | If '''Start''' and/or '''End''' are left empty, the start and/or end of the entire simulation time period is used. | ||
|- | |- | ||
| '''RefName'''<br/>( | | '''RefName'''<br/>(optional) || Specifies the name of the reference time series. Only necessary if the time series file specified under '''RefFile''' contains multiple series (e.g. [[WEL-format|WEL file]]). | ||
|- | |- | ||
| '''RefFile''' || The filename of the reference time series | | '''RefFile''' || The filename of the file containing the reference time series<br/><span style="color:red">'''Warning:''' The time step size of the reference series must be identical to that of the simulation result! (see [https://github.com/bluemodel/BlueM.Opt/issues/151 #151])</span> | ||
|} | |} | ||
==="Values"=== | ==="Values"=== | ||
Line 245: | Line 229: | ||
|- | |- | ||
| '''ValueType''' || Specifies how to calculate the value from the simulation result series. Supported values: | | '''ValueType''' || Specifies how to calculate the value from the simulation result series. Supported values: | ||
* <code> | * <code>Max</code>: Maximum value of the simulation result | ||
* <code> | * <code>Min</code>: Minimum value of the simulation result | ||
* <code> | * <code>Avg</code>: Average value of the simulation result | ||
* <code> | * <code>Start</code>: First value of the simulation result | ||
* <code> | * <code>End</code>: Last value of the simulation result | ||
* <code> | * <code>Sum</code>: Sum of the simulation result | ||
|- | |- | ||
| '''RefValue''' || The reference value with which the computed simulation result should be compared. | | '''RefValue''' || The reference value with which the computed simulation result should be compared. | ||
|} | |} | ||
==Literature== | ==Literature== |
Latest revision as of 23:43, 4 October 2023
BlueM.Opt | Download | Usage | Development
BlueM.Opt input files: MOD · OPT · OBF · CON
Versions:
This file defines the objective functions for the optimization.
The file must have the same filename as the rest of the dataset's files but have the file extension .OBF
, and be located in the same directory as the rest of the dataset. (In version prior to v1.8, this file had the file extension .ZIE
.)
The file is divided into several blocks, each of which is optional. Each block is identified by its heading:
- Series: This block defines objective functions that are based on a comparison of simulation time series with a reference time series.
- Values: This block defines objective functions that are based on a comparison of simulation result values with a reference value (only for SWMM applications).
- ValueFromSeries: This block defines objective functions that are based on a comparison of values that are computed from simulated time series with reference values.
- Aggregate: Aggregate objective functions
File format
*Objective functions *=================== * *Series *------ * *|-----|------------------|---------|---|-------|-------|------------|---------|------------------|------------------|-----------------------------|----------| *| Opt | Name | Group | D | OpFac | File | SimResult | ObjFunc | Evaluation period | Reference series | CurValue | *| | | | | | | | | Start | End | RefName | RefFile | | *|-----|------------------|---------|---|-------|-------|------------|---------|------------------|------------------|-----------|-----------------|----------| | P | Kling-Gupta | | + | | WEL | B401_1AB | KGE | 30.05.1984 00:00 | 18.06.1984 23:30 | | seriesfile.zre | | | S | NashSutcliffe | | + | | WEL | B401_1AB | NSE | 30.05.1984 00:00 | 18.06.1984 23:30 | | seriesfile.zre | | | S | VolumeError | | - | | WEL | B401_1AB | BIAS | 30.05.1984 00:00 | 18.06.1984 23:30 | | seriesfile.zre | | | S | MeanSquaredError | | - | | WEL | B401_1AB | MSE | 30.05.1984 00:00 | 18.06.1984 23:30 | | seriesfile.zre | | *|-----|------------------|---------|---|-------|-------|------------|---------|------------------|------------------|-----------|-----------------|----------| * *Values *------ * *|-----|---------------|---------|---|-------|-------|-----------|---------|---------------|-------------|----------|----------| *| Opt | Name | Group | D | OpFac | File | SimResult | ObjFunc | TargetElement | TargetValue | | CurValue | *| | | | | | | | | (Block) | (Column) | RefValue | | *|-----|---------------|---------|---|-------|-------|-----------|---------|---------------|-------------|----------|----------| | S | OutfallLoad | | - | | RPT | | AE | OutfallNode | AvgFlow | 0 | | | S | Storage401 | Storage | - | + | BLZ | B401 | | Speicher | maxVolume | | | | S | Storage701 | Storage | - | + | BLZ | B701 | | Speicher | maxVolume | | | *|-----|---------------|---------|---|-------|-------|-----------|---------|-------------- |-------------|----------|----------| * *ValueFromSeries *--------------- * *|-----|---------------|---------|---|-------|-------|------------|---------|---------|---------|----------------------|----------| *| Opt | Name | Group | D | OpFac | File | SimResult | ObjFunc | Evaluation period | Reference value | CurValue | *| | | | | | | | | Start | Ende | ValueType | RefValue | | *|-----|---------------|---------|---|-------|-------|------------|---------|---------|---------|-----------|----------|----------| | S | OutfallLoad | | - | | WEL | B401_1AB | AE | | | | | | | P | PeakError | | - | | WEL | S3 _1ZU | SE | | | Max | 29.33 | | *|-----|---------------|---------|---|-------|-------|------------|---------|---------|---------|-----------|----------|----------| * *Aggregate *--------- * *|-----|--------------|---------|---|----------| *| Opt | Name | Group | D | CurValue | *|-----|--------------|---------|---|----------| | P | TotalStorage | Storage | | | *|-----|--------------|---------|---|----------|
The width of the columns can be adjusted if necessary, but the column separators (|
) must be preserved!
An asterisk (*
) at the beginning of the line causes the line to be ignored.
Description
Common parameters
The blocks "Series", "Values" and "ValuesFromSeries" share a set of parameters.
Parameter Description Opt Specifies, whether the objective function is a primary objective function ( P
) or a secondary objective function (S
).
Primary objective functions are used to calculate the fitness of invidual solutions.
Secondary objective functions are evaluated, but are not used to determine the fitness of individuals. Secondary objective functions are useful for specifiying objective functions that are part of an aggregated objective function or for specifying values of interest that should be evaluated, but not optimized towards.Name Arbitrary name of the objective function Group
(optional)Optional group name, only used for calculating aggregated objective functions. Objective functions with the same group name are aggregated to form one aggregated objective function. Direction
(optional)Direction of desired optimization:
-
: The objective function should be minimized (default)+
: The objective function should be maximized (this causes the objective function values to be multiplied by-1
internally).
Operator or Factor
(optional)Optional, only used for aggregated objective functions
+
: The objective function value is added to the group-
: The objective function value is subtracted from the group- number: The objective function value is multiplied by the given number and then added to the group
File The file extension of the simulation result file containing the simulation result from which the objective function should be calculated. Supported values: SimResult The identifier of the simulation result from which the objective function should be calculated. ObjFunc Specifies which objective function to use. Supported values for comparing values:
SE
: Squared errorAE
: Absolute error
Supported values for comparing series:
MAE
: Mean abolute errorSAE
: Sum of absolute errorsMSE
: Mean squared errorSSE
: Sum of squared errorsNSE
: Nash-Sutcliffe efficiency[1]
-∞ ≤ NSE ≤ 1
. Larger is better.- [math]\displaystyle{ NSE=1-\frac{\sum_{t=1}^T\left(Q_{obs}^t-Q_{sim}^t\right)^2}{\sum_{t=1}^T\left(Q_{obs}^t-\overline{Q_{obs}}\right)^2} }[/math]
lnNSE
: Logarithmic Nash-Sutcliffe efficiency[1]
-∞ ≤ lnNSE ≤ 1
. Larger is better.- [math]\displaystyle{ lnNSE=1-\frac{\sum_{t=1}^T\left(ln(Q_{obs}^t+\epsilon)-ln(Q_{sim}^t+\epsilon)\right)^2}{\sum_{t=1}^T\left(ln(Q_{obs}^t+\epsilon)-\overline{ln(Q_{obs}+\epsilon)}\right)^2} }[/math]
ε
is set to 1% of the average value ofQobs
as recommended by Pushpalatha et al. (2012)[2].KGE
: Kling-Gupta efficiency[3]
-∞ ≤ KGE ≤ 1
. Larger is better.- [math]\displaystyle{ \text{KGE} = 1 - \sqrt{ (r - 1)^2 + (\beta - 1)^2 + (\gamma - 1)^2 } }[/math]
- with
r
: correlation coefficientβ
: bias ratioγ
: variability ratio
NGT
: Relative number of timesteps in which the simulation result is larger than (exceeds) the value of the reference time series [%]SGT
: Sum of the simulation values that are larger than the values of the reference time seriesNLT
: Relative number of timesteps in which the simulation result is smaller than the value of the reference time series [%]SLT
: Sum of the simulation values that are smaller than the values of the reference time seriesDET
: Coefficient of determination (simple linear regression)0 ≤ R² ≤ 1
:
- [math]\displaystyle{ R^2 = \frac{s^2_{xy}}{s^2_x \cdot s^2_y} }[/math]
BIAS
: Absolute volume error [%]
The following options are deprecated and only kept for backwards compatibility:
NashSutt
: Modified Nash-Sutcliffe efficiency
- modified formula with
0 ≤ NSEmod ≤ ∞
. Smaller is better. - [math]\displaystyle{ NSE_{mod}=1-NSE }[/math]
LnNashSutt
: Modified logarithmic Nash-Sutcliffe efficiency
- modified formula with
0 ≤ lnNSEmod ≤ ∞
. Smaller is better. - [math]\displaystyle{ lnNSE_{mod}=1-lnNSE }[/math]
CurValue
(optional)An optional value that is displayed in diagrams to serve as an indication of improvement.
"Series"
Parameter Description Evaluation period
(optional)Specifies the simulation time period for which the objective function should be evaluated.
Format:
dd.mm.yyyy hh:mm
If Start and/or End are left empty, the start and/or end of the entire simulation time period is used.RefName
(optional)Specifies the name of the reference time series. Only necessary if the time series file specified under RefFile contains multiple series (e.g. WEL file). RefFile The filename of the file containing the reference time series
Warning: The time step size of the reference series must be identical to that of the simulation result! (see #151)
"Values"
Parameter Description TargetElement (Block) Specifies in which section (block) the target variable is located. TargetValue (Column) Specifies the target variable within the block. RefValue The reference value with which the simulation result should be compared.
The parameter TargetElement is dependent on the parameter TargetValue and the simulation tool used, only the following combinations are possible:
Application TargetElement (Block) TargetValue (Column) SWMM NodeFlooding HoursFlooded SWMM NodeFlooding FloodVolume SWMM StorageVolume AvgVolume SWMM StorageVolume AvgPctFull SWMM StorageVolume MaxVolume SWMM StorageVolume MaxPctFull SWMM StorageVolume MaxOutflow SWMM OutfallLoad MaxFlow SWMM OutfallLoad FlowVolume SWMM OutfallLoad Pollutant_[no] SWMM Pumping OnlineTime SWMM Pumping TotalEnergy Smusi EntlVolumen SumVol Smusi MaxAbfluss Qmax Smusi EntlFracht CSBspez
The TargetValue Pollutant_[no] allows for the evaluation of up to five substances, [no] can take values of 01, 02, 03, 04 or 05, so e.g. specify Pollutant_02 for the second substance.
"ValueFromSeries"
Parameter Description Evaluation period
(optional)Specifies the simulation time period for which the objective function should be evaluated.
Format:
dd.mm.yyyy hh:mm
If Start and/or End are left empty, the start and/or end of the entire simulation time period is used.ValueType Specifies how to calculate the value from the simulation result series. Supported values: Max
: Maximum value of the simulation resultMin
: Minimum value of the simulation resultAvg
: Average value of the simulation resultStart
: First value of the simulation resultEnd
: Last value of the simulation resultSum
: Sum of the simulation result
RefValue The reference value with which the computed simulation result should be compared.
Literature
- ↑ 1.0 1.1 Nash, J. E. and Sutcliffe, J. V. (1970): River flow forecasting through conceptual models part I — A discussion of principles, Journal of Hydrology, 10 (3), 282–290, DOI:10.1016/0022-1694(70)90255-6.
- ↑ Pushpalatha, Raji, Perrin, Charles, Le Moine, Nicolas, Andréassian, Vazken (2012): A review of efficiency criteria suitable for evaluating low-flow simulations, Journal of Hydrology, Volumes 420–421, 2012, Pages 171-182, ISSN 0022-1694, https://doi.org/10.1016/j.jhydrol.2011.11.055.
- ↑ Kling, H., Fuchs, M., Paulin, M. (2012): Runoff conditions in the upper Danube basin under an ensemble of climate change scenarios, Journal of Hydrology, Volumes 424–425, Pages 264-277, ISSN 0022-1694, https://doi.org/10.1016/j.jhydrol.2012.01.011.