一、需要率定的參數(shù)（對(duì)徑流模擬效果敏感的土壤參數(shù)）

以下參數(shù)，一開(kāi)始沒(méi)有完全準(zhǔn)確的值，需要在一個(gè)范圍內(nèi)率定直到獲得最優(yōu)參數(shù)組合：

infilt：可變下滲曲線方程的冪，通常為0.2~0.6

Ds：非線性基流發(fā)生時(shí)流速占基流最大流速Dsmax的比例

Dsmax：基流最大流速，單位mm/day

Ws：非線性基流發(fā)生時(shí)土壤含水率占飽和含水率的比例

depth：除第一層以外的其他層土壤厚度（即第2、3層），單位m

二、前期準(zhǔn)備工作

1.將打開(kāi)"lvding"文件夾放入產(chǎn)匯流文件夾"F_F"中

2.打開(kāi)start.bat腳本文件（"lvding\code\start.bat"），修改文件路徑（參考之前產(chǎn)流匯流），注意斜杠方向的修改"/"

3.打開(kāi)"NashCoef.m""CreateNewFile.m"腳本文件（lvding文件夾中）及"MinFun.m"腳本文件（lvding\SCEUA文件夾中）

①NashCoef.m：修改第3行的文件路徑（匯流模型輸出的"逐日徑流"數(shù)據(jù)）

②MinFun.m：修改第6、7行的文件路徑（"單位線"數(shù)據(jù)）

③CreateNewFile.m：按照腳本內(nèi)"行數(shù)"修改"4_Sum_soil_param.R"（lvding\code\R\）

"4_Sum_soil_param.R"修改如下：保證這6個(gè)參數(shù)分別在7，8，9，10，37，38行

注意：還需修改"4_Sum_soil_param.R .sample"

4.打開(kāi)"new 3.txt"文本（F:\VIC_modle\result\F_F）

第1行，6個(gè)參數(shù)的取值范圍，需根據(jù)研究區(qū)來(lái)定（查閱文獻(xiàn)），可能需要不斷修改率定

參考文獻(xiàn)：王媛. 基于VIC模型的黃河源區(qū)降雨徑流轉(zhuǎn)化關(guān)系模擬研究[D].青海大學(xué),2022.

三、率定（在Matlab中打開(kāi)以下五個(gè)腳本）

1.在"NashCoef.m"腳本中，將實(shí)際平均徑流更改為研究區(qū)歷年平均徑流

2.在"AutoOptimize.m"腳本中

①在命令行窗口輸入real=0，然后單機(jī)回車，右邊工作區(qū)會(huì)出現(xiàn)一個(gè)real變量

②單擊real打開(kāi)，將實(shí)際徑流數(shù)據(jù)拷貝至表格中

③在命令行窗口輸入"save+文件名"，單機(jī)回車，即可在文件夾中看到該文件，將其拖入resource文件夾中

④更改腳本中的文件名

3.在自動(dòng)率定前，先進(jìn)行手動(dòng)測(cè)試

①運(yùn)行"AutoOptimize.m"腳本

②在命令窗口中輸入MinFu ( [ , , , , ,?] )，括號(hào)中是六個(gè)參數(shù)的值，單機(jī)回車運(yùn)行

4.自動(dòng)率定

①將new 3.txt"內(nèi)容拷貝至Matlab的新腳本中，全選執(zhí)行所選內(nèi)容

②打開(kāi)右邊工作區(qū)出現(xiàn)"V變量"，選擇已經(jīng)組合好的參數(shù)矩陣，將其拷貝至"params.txt"文件中（lvding\resource）

③開(kāi)始自動(dòng)率定，先運(yùn)行"AutoOptimize.m"腳本，再運(yùn)行"ParamsCalculate.m"

Error1：CMD執(zhí)行異常

原因：此路徑下，缺少uh_s文件，之前被誤刪，加上即可，首先要有這個(gè)文件，才能在此命令下刪除

Error2：無(wú)論參數(shù)如何變化，Nash效率系數(shù)一直為1

原因1：我隨意設(shè)置的平均徑流太大，所以nash為1，我將平均徑流修改為真實(shí)合理值，發(fā)現(xiàn)Nash有變化，不再是1，但依然不會(huì)隨著參數(shù)的設(shè)置而變化

在我百思不得其解時(shí)，同學(xué)的一句話點(diǎn)醒了我（感謝他?。。。?strong>重新跑了產(chǎn)流和匯流模型，發(fā)現(xiàn)真的存在問(wèn)題?。。槭裁次揖蜎](méi)有想到去檢查呢，我還很自信的認(rèn)為之前產(chǎn)流匯流都調(diào)試成功了，不用再回頭看/(ㄒoㄒ)/~~

原因2：產(chǎn)流報(bào)錯(cuò)原因是bulk density必須大于0，在之前第一次跑的時(shí)候遇到過(guò)（《學(xué)習(xí)記錄十》），因?yàn)楫?dāng)時(shí)直接在土壤輸出文件中粗暴的把0改為0.1，但代碼沒(méi)改過(guò)程輸出的文件沒(méi)改，所以導(dǎo)致這個(gè)問(wèn)題依然存在

原因3：匯流報(bào)錯(cuò)原因是"單位線文件已存在"，此前這個(gè)問(wèn)題也遇到過(guò)，但率定時(shí)，不用在"F_F_stnloc.txt"中加單位線地址，直接NONE即可，因?yàn)槁识ùa里已經(jīng)有刪除命令

運(yùn)行成功?。?！

四、一些注意事項(xiàng)

1.首先要保證匯流和產(chǎn)流是能跑通，沒(méi)有報(bào)錯(cuò)的，才能進(jìn)行率定，否則看似率定出了nash系數(shù)，但其實(shí)是錯(cuò)的，一定要多檢查?。?！再放一次這幾句點(diǎn)醒我的話?。。?/strong>

2.檢查輸出單位是否是"立方米/秒"和"mm"，如果你的"make_convolution.f"中和下圖一樣，那就是已經(jīng)轉(zhuǎn)換過(guò)了，此文件不用再次修改

3.自動(dòng)率定原本代碼只能用逐日實(shí)測(cè)徑流率定，輸入逐月無(wú)法率定，需要改動(dòng)代碼

4.要在一開(kāi)始就設(shè)定好VIC模型的預(yù)熱期、率定期、驗(yàn)證期，進(jìn)行VIC模擬時(shí)得分時(shí)段進(jìn)行模擬。由于剛開(kāi)始模擬計(jì)算時(shí)，流域的初始狀態(tài)未知，所以要設(shè)定一個(gè)預(yù)熱期，以便提高模擬的準(zhǔn)確性；然后在率定期對(duì)相關(guān)參數(shù)進(jìn)行率定，確定最優(yōu)參數(shù)組合；最后用驗(yàn)證期對(duì)模型模擬效果進(jìn)行驗(yàn)證。不要?jiǎng)傞_(kāi)始就傻乎乎一股腦把所有年份一起模擬了，比如我......

5.不是每個(gè)流域都可以用VIC，有些流域人類活動(dòng)干擾很大，用VIC模擬效果不好，比如我正要做的流域......聽(tīng)天由命咯，不行就得換模型了/(ㄒoㄒ)/~~

6.關(guān)于VIC的文獻(xiàn)不是很多，以下幾篇碩博論文中VIC講的比較詳細(xì)可參考

[1]李秀清. 基于VIC模型的丹江流域水文模擬及水資源管理對(duì)策[D].西北大學(xué),2021.

[2]孟鵬飛. 基于水文模擬定量評(píng)估漳河上游流域徑流變化歸因[D].河北工程大學(xué),2020.

[3]王媛. 基于VIC模型的黃河源區(qū)降雨徑流轉(zhuǎn)化關(guān)系模擬研究[D].青海大學(xué),2022.

[4]蔡潔連. 基于VIC模型的衛(wèi)星降水產(chǎn)品在贛江流域的應(yīng)用研究[D].南寧師范大學(xué),2022.

自我反省及感謝Part? ?多看文獻(xiàn)，多思考，多檢查，不要想當(dāng)然?。。。?/span>

這周在做率定的時(shí)候遇到了一些問(wèn)題，在咸魚(yú)找人解疑，要么太貴，要么回復(fù)慢。回想起學(xué)習(xí)VIC的整個(gè)歷程，遇到不少問(wèn)題，每次遇到的問(wèn)題可能很簡(jiǎn)單，但自己總是不開(kāi)竅，需要?jiǎng)e人點(diǎn)一下才懂，并且剛解決完一個(gè)問(wèn)題就接著遇到下一個(gè)問(wèn)題，在自閉中輪回折磨。所以在跑到率定這一步，我非常想立刻跑出來(lái)，完成它?。?！因?yàn)槁识ㄒ恢背霈F(xiàn)問(wèn)題，我連著幾個(gè)晚上失眠到凌晨4點(diǎn)......最后不得已還是去打擾了一位同學(xué)，正巧他告訴我學(xué)習(xí)這個(gè)模型已有一年已經(jīng)跑通，并且愿意幫助我，帶病為我語(yǔ)音講解，沒(méi)有任何推脫，在他的點(diǎn)撥下，我最后才得以跑通！?。》浅８兄x?。。∽Ｋ厴I(yè)順利工作順利??！! 也感謝某位有償解疑答惑的大佬，在別人收200-450不等的時(shí)候，他每次只收20塊，一句回復(fù)就能解決我的問(wèn)題?。?！

因?yàn)閂IC模型學(xué)習(xí)資料網(wǎng)絡(luò)上真的太少了，相關(guān)文獻(xiàn)里也有步驟，但比較簡(jiǎn)略，沒(méi)有具體操作步驟，我的信息渠道也很閉塞，在某魚(yú)找人問(wèn)一個(gè)小問(wèn)題動(dòng)輒就收費(fèi)幾百，所以在此分享記錄我的VIC學(xué)習(xí)歷程，希望對(duì)也在苦苦學(xué)習(xí)這個(gè)模型的同學(xué)有用，大家一起交流共同進(jìn)步?。?！如果你在看這些文章的時(shí)候，發(fā)現(xiàn)我的操作步驟有問(wèn)題，或者對(duì)我遇到的問(wèn)題有好的解決方案，也希望可以告知我！?。≡u(píng)論私聊均可?。?！

OK，算是把VIC模型通學(xué)了一遍，我現(xiàn)在的狀態(tài)還是一知半解，下面我就開(kāi)始好好讀文獻(xiàn)，分好預(yù)熱/率定/驗(yàn)證期，不要盲干了，理解這個(gè)模型后，再?gòu)念^開(kāi)始跑一遍！?。?/p>

標(biāo)簽：