如果估算的參數(shù)有100個,那么每進(jìn)行一輪迭代,PEST要運(yùn)行100次,如果參數(shù)有成千上萬個的話,實(shí)際計算時間過程非常長甚至不可行。PEST中SVD-Assist和Parallel PEST就是為了減少運(yùn)行計算時間而設(shè)置的。SVD-Assist選項(xiàng)識別方向模型中的不敏感參數(shù)、然后移除這些參數(shù),Parallel PEST可以對模型進(jìn)行并行計算從而加速模型計算過程。
SVD-Assist得以實(shí)現(xiàn),其內(nèi)部原理主要包含3個基本步驟:
- 1)為了計算矩陣,PEST首先為每個參數(shù)運(yùn)行一次MODFLOW。計算的結(jié)果用來創(chuàng)建超參數(shù),這些超參數(shù)來源于初始指定的參數(shù)集。
- 2)運(yùn)行SVDAPREP來創(chuàng)建新的PEST控制文件。MODFLOW|Parameter Estimation——SVD-Assist Options下的選項(xiàng)是有關(guān)SVDAPREP的選項(xiàng)。其中最重要的選項(xiàng)是Specify # super param(默認(rèn)是No),當(dāng)這個選項(xiàng)為No時,寫入到SVD文件中的信息將用來指定超參數(shù)的數(shù)量。
- 3)使用SVDAPREP寫入的控制文件來運(yùn)行PEST程序。這樣,PEST在每次迭代時,模型運(yùn)行次數(shù)會顯著減小。PEST運(yùn)行結(jié)束時,PAR文件包含的時超參數(shù)的值,而基礎(chǔ)參數(shù)的值包含在BPA文件中,因此,導(dǎo)入最優(yōu)參數(shù)值時,應(yīng)選擇 .bpa文件
當(dāng)再次運(yùn)行模型時,可以通過指定Jacobian和SVD文件,而讓上述第1)步驟省去,具體方式為MODFLOW | Parameter Estimation——SVD-Assist Options——Specify Jacobian file——找到上次模型運(yùn)行產(chǎn)生的文件路徑+文件名(.jco文件),填入到Jacobian file path,找到文件路徑+文件名(.svd文件),填入到SVD file path