PEST中應(yīng)用SVD-Assist和Parallel PEST

如果估算的參數(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


最后編輯于
?著作權(quán)歸作者所有,轉(zhuǎn)載或內(nèi)容合作請聯(lián)系作者
【社區(qū)內(nèi)容提示】社區(qū)部分內(nèi)容疑似由AI輔助生成,瀏覽時請結(jié)合常識與多方信息審慎甄別。
平臺聲明:文章內(nèi)容(如有圖片或視頻亦包括在內(nèi))由作者上傳并發(fā)布,文章內(nèi)容僅代表作者本人觀點(diǎn),簡書系信息發(fā)布平臺,僅提供信息存儲服務(wù)。

相關(guān)閱讀更多精彩內(nèi)容

友情鏈接更多精彩內(nèi)容