14.07.2013 Views

User's guide of Proceessing Modflow 5.0

User's guide of Proceessing Modflow 5.0

User's guide of Proceessing Modflow 5.0

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

148 Processing <strong>Modflow</strong><br />

If a parameter b is relative-limited, The relative change between optimization<br />

iterations i-1 and i is defined as<br />

(b i&1 & b i )<br />

b i&1<br />

b i&1 /b i if |b i&1 | > |b i |, or<br />

b i /b i&1<br />

3.6.5 PEST (Inverse Modeling)<br />

if |b i | > |b i&1 |<br />

(3.55)<br />

The absolute value <strong>of</strong> this relative change must be less than RELPARMAX. If a<br />

parameter upgrade vector is calculated such that the relative adjustment for one or more<br />

relative-limited parameters is greater than RELPARMAX, the magnitude <strong>of</strong> the upgrade<br />

vector is reduced such that this no longer occurs.<br />

If parameter b is factor-limited, the factor change between optimization iterations i-1<br />

and i is defined as<br />

(3.56)<br />

This factor change must be less than FACPARMAX. If a parameter upgrade vector is<br />

calculated such that the factor adjustment for one or more factor-limited parameters is<br />

greater than FACPARMAX, the magnitude <strong>of</strong> the upgrade vector is reduced such that<br />

this no longer occurs. It is important to note that a factor-limit will not allow a parameter<br />

to change sign. If a parameter must be free to change sign during an optimization<br />

process, it must be relative-limited; furthermore RELPARMAX must be set at greater<br />

than unity or the change <strong>of</strong> sign will be impossible. Similarly, if a parameter's upper or<br />

lower limit is zero, it cannot be factor-limited and RELPARMAX must be at least unity.<br />

Suitable values for RELPARMAX and FACPARMAX can vary enormously from<br />

case to case. If you are unsure <strong>of</strong> how to set these parameters, a value <strong>of</strong> 5 for each <strong>of</strong><br />

them is <strong>of</strong>ten suitable. For highly non-linear problems, these values are best set low. If<br />

they are set too low, however, the estimation process can be very slow. An inspection <strong>of</strong><br />

the PEST run record (by pressing the ESC key) will <strong>of</strong>ten show whether you have set<br />

these values too low, for PEST records the maximum parameter factor and relative<br />

changes are recorded on this file at the end <strong>of</strong> each optimization iteration. If these<br />

changes are always at their upper limits and the estimation process is showing no signs<br />

<strong>of</strong> instability, it is quite possible that RELPARMAX and/or FACPARMAX are to low<br />

and could be increased.<br />

Note that FACPARMAX can never be less than 1; RELPARMAX can be less than<br />

1 as long as no parameter's upper and lower bounds are <strong>of</strong> opposite sign. (If necessary,<br />

use OFFSET to shift the parameter domain so that it does not include zero.)<br />

- FACORIG is a criterion for modifying RELPARMAX and FACPARMAX. If, in the

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!