(完整版)Hydrus-1D中文說明

1、天山北麓平原區(qū)包氣帶水分運移機理與數(shù)值分析HYDRUS1D軟件使用手冊目錄一、軟件介紹二、安裝指南三、使用指南四、實例介紹 一、軟件介紹HYDRUS1D是國際地下水模型中心公布的，計算包氣帶水分、鹽分運移規(guī)律的軟件，用它可以解算在不同邊界條件制約下的數(shù)學(xué)模型。若將坐標原點選在地面，取z軸向下為正，則一維飽和非飽和帶水分運移基本方程為：理丄|k（6）（也-1）1-Sdtdzdz_式中：e為含水率；h為負壓水頭；S為植物根系吸水量，對裸露區(qū)為0。二、安裝指南第一步：安裝界面；2RfDIITi-ID51dtHYDRVS-1DVMecwWtwH別如冷聲IWi帥，|如刊IH臨The.崎仍打ht-tod1

2、Jraw*ra*nrf!mbdrtbAmd1btLKf.r阿0眄|呻郵怕c崗RewjM*嶺廚陽電第二步：選擇安裝目錄；bircniiilb泌”士tEwidfdiridrcdvtriiwgunijhx朗aB4HiHrTRJ5-li臥出宵iMdioIjhvhe.I第三步：修改存放名稱；個工JLBVjhLTZ*加削FoHh.HYDRUS-SDwW朗蝕mr如博diwFft那Fc*i-4TBft.iTsrtcrgjhj!3冋Ml：kJCkl-H*JCfiWu第四步：安裝進度顯示條；Tie3*3iTcVyciusld*f/r/dnaid.d*l：!CTDU1Pl臥CfilIIYDRUS-IDIg|第五步：

3、完成安裝；三、使用指南一）初始界面軟件開始界面如圖所示：Ietdrus-id-展試UTlgfx1114VaawFrc-procsssiiiAC-alcilGtr：KASTiltsFViniiDMHalp間喝I副國電舊血I丫丨別呵IPre-prc3l口二P3t-rrnce35LILC田#Ma.inFmeeu呂./GepmeiryInu-n-3.t.ion耆#TinaHitfoftlDn昌FrlntInfornatLan1ji|WatorJLcm-Itoi4.ticoiCritori4.WaterFl酬-SoilM/dr泗lieFropertyNodel1IWitcrFLoh-SoilHyxhr4

4、iLicPdr-rnetiers1JHaWtrFIciMi-BauridAryCondiLion-zF*(V-iriatliBDunduryCand.3-tiansP|Soilfrafila-Gnphicd.Eiitar岡|SciLProfile-SmnaryFoorJfelp,prcssfl|IfLM|二）數(shù)據(jù)輸入及運行1、選擇計算模塊2、模型幾何結(jié)構(gòu)屬性nHTDRUS-1D-黏土JQ岡Fllc卑刑Frc-proc&ee3Calculfl.tidnEaeuI1cXiridowIfelp聘口|綱%|s|m|?|(o|HGllFrbCAECEF口3l-pro匚亡rsixifGeonttryInt

5、oirn&tsnnrinbInEcrnbtianlYlsr1.PrnnlGe口Et:ryDifDTBarti口nVaterIVaierTL-dilLengthUnitenmanij*n3、時間信息WBmo期ProfilaTnEorn-tionfrater71c*-口皿&叱fFlwiesaridHeadsSailKydrinLicPlropMrliEhJfmbercfEoiLlVl.tc2*3&lsKtimbarofLdybirEfarIHaeeJO4、選擇輸出結(jié)果HTDRU5-ID-祐Jnlxi|1FilaViavFrflprfleqeelngCalinilAtionKbulLeWindowH

6、alpn圍樹口|顧皂口e竜z=Pit-prQUddLng訝訝自呂TO圈BoundaryPrintTinitsrrmiDpiLooisI_LeTtlInfarmatInforri3.ilgflaw-ItcarMa.ldPracAEtiKGeomeiryInf4rnatTinaInfarn-atianPrintWaterWaterWaterWaterEduxl血ryCmiliSoilPrafil*-Gr&phic-iSolLPrfalt-E-umniiyQutpuiVjjnliErqEPrinibpdFtldEF&rMtlp,prtzs715、迭代步長L&werLptirridlItcratloh

7、HengeDptimiLTlarilLon.KangaLoerTimetepFliHtipiicationVpperFlow-F1d1-lateruterTiraeStephlullLplj:bIjodEnformitIiEornEFl州-FIdiy-IatarTriwrViuri4bL4BothSoilProfileS1Profilefluiniui-DJfijrriDsrot11aratlcue?aterContent.TolermceInterrialIuterpliLCJiFaliLsLinitdIksI&nsi&r：L4-006LCCOOUpperLiinitoftheTension

8、irrvi6、土壤水運移方程7、巖性參數(shù)8、上、下邊界條件確定nHTDJtUS-lD-黏土1FileJfiTrc-proaeeeingCdl2ijlati4iLIsYirid&wPferpiroc戸：5*：咬Va.terFlc-wBDimiliiTyf口nd.ixinxMannrrDcaEcEGeonetryInternrinal3iEcorn&.tii:PrintInfnrina-I.LIVatsrILcm-ItWaterFlow-SoWaitfnfl-S&nffll&rTLc-BdarL-sblEBQUTidaSoiLFrofila-SoilFrnfile-UpparBouiid.-ury

9、CanditiqlrCoustantPtrxw口tpTia汩AtnospticflcECviittiSurfL-rAlnospticsrLCBCthSurfa.ac童皿01VirKbLcrrcssuorJia4.dVriibLcFrszjareMad/FLoH-LoverBoodaryConditionCjuisimtConsint常VaricLbLc-VAriibLc*PressureHea.d1jZ3HTadHedsJ-Luk7LuhFreeDrsinaeeHeepDrainatpd.gPb.CMoriaontslDrbinsOForHelp,pressFL9、上、下邊界條件數(shù)據(jù)輸入10、

10、幾何結(jié)構(gòu)確定11、初始條件zhKootAezDxz1st10E8870111120.05-0.3H4350111131111IIi11140:00501111EII”D-II111160.250150111170.3D.201111BUJ:UJbu111190.：.30-11110045035II1111110.6D.4011111PIIT114.IIi111130.：.5011111也0G5055II11111!E0.7D.601111匚二二工二二二CancelEfcvlous|fieztL=lk12、模型運行（三）運行結(jié)果1、剖面水頭、含水率等計算結(jié)果2、邊界流量信息LDHORUS-1Df

11、ileYitmriBaimdaryfatcrTinafjninEVacseebG-bondtryInft1LmeLdde泌FriutInforncVaterFl州-nterKL&iv-VatarFldiv-&terFI&iy-YuriihlaBdhj3、參數(shù)屬性HTDRUS-lIii-IFilmViwTyi口圍風亙j-FlE-pTQEZULLngMildPracAE-EtGesneiryInee沾aridHeais7ZZ/hbeffibefiriorizontaJ.VarticalLl_二=rHa九fTDC45ZGGfonetryInfcliniEInEornPrintInEornsiVfii

12、erFLw-RffllarFLw-WalerTLom-IVqIsrTLom-arLatleBotitgoalfrofilgSoiLFrofiliRun-TimeTineTTin亡Step*vrinaSlapfelmbsifiItdrtiTis間簾帥電羈權(quán)冊卿9KW擁RRWWHMHTime|daysg+-土:-.ForHpLrptfpwFL5、水均衡信息HU-US1Tdnerrj.ODOOSixb-fefdontivuxi.1At&b.LL900001401.3LI0U0E401W-valun.aLL3B421+01.35427E401In-fLauLL/TJooooaE+oa.0U0U0E-F

13、00hHsanLL.433611+CIJ.433E7E4-口J.TopffiunLL/T.EZS口占3-口：3Ear.ffluK1T-.1359dI-03四、實例操作流程介紹1、實例介紹對奎屯原位試驗場包氣帶裸露區(qū)非均質(zhì)水分運移規(guī)律進行模擬。剖面巖性結(jié)構(gòu)從上到下共7層；包括：壤土、粉壤土、粉質(zhì)黏壤土以及砂質(zhì)壤土4種巖性類型；水位埋深最大不超過9m；模型運行期為4月1號至11月20號，共234天。選擇萬方程作為模型計算的方程；上邊界為變化的負壓邊界，下邊界為變化的水頭邊界；初始條件以4月1號觀測的負壓為準輸入模型進行計算；輸出每月最后一天的結(jié)果信息。2、軟件操作流程：-SinmlateIk-lr

14、nlDLysc-n4tryIrifoim4.114：TimaInfQkTriiticin.FriuiiAtomiixtionIatarFldw-rtirataterFlSoilK：int4rFl-SaiLKVaterFl州-BoundiV&riibl亡BoruidaryCSoilPrflEili-Grap；SilProFile-SnrriruSoLuteTrarisp&rtTf-2nspartIHootVaterUptake廠HootGrcvth.Fost-proctssnncEz訂Fz呂口gE第一步：軟件可以進行計算的模塊包括水分運移模塊、溶質(zhì)運移模塊、熱傳導(dǎo)模塊、植物根系吸收水分模塊以及植

15、物根系生長模塊等，本次模擬只對包氣帶水分運移規(guī)律進行模擬，故選擇水流模塊，操作如下圖所示：HHmtUS-lD-世檢試豔璋*LinlAHiE訂匕VliErTrocesiiLEQd-ciilatiorifiesyiltsHindoo心眉圜釦%|S|ffl|?|Jfrfrir&cessinr第二步：選定計算長度單位m,剖面巖性類型4種，從上到下共分為7層,計算深度為9m,操作如下圖所示:第二步：選定計算時間單位Days，計算時間從4月1號至11月20號，共234天,邊界條件隨時間變化記錄的最小間隔為1天,故隨時間變化的邊界條件記錄數(shù)為234天,操作如下圖所示：IT-!TieEnfarMaticiii

16、!and.HaAs:MainFrocadTigDiEcratiz-atiddIniiL-alTineFrlntInforJLawCcaacJJinlineInitL-slHlninijTi001FlowVeriatliE:HwinuniTie電SispSailFrflfilSailPra:lFlovFlavWaterWaterLepineiryInTinaInfarnit4rofriras-ViriabliEurLilNy(.口:S.econdsNinutas:Kfl.nr3荷BiyzriunbercisoroLciE-tl広.遷.Rtdi.M11jl.；BonridiryCond.L-tio

17、nzTlrid-V-urjakda創(chuàng)uDdaryCondjtio廠feiscroLuEi匚辺Data.廠DailyVarjationsBoringTisyS-nerAteibyI匕HYDKUS-ID-爺也姑鮭場*Z回岡JlilaV:sivTrfl-pracqeelngLalcvl.tioilKarultzIViniiA,Hdlpq|Q|知3,sm?Fte-口rue亡三ring第三步：要求輸出每月最后一天的計算結(jié)果，統(tǒng)計共有9天的計算結(jié)果，操作如下圖所示：HYDRUS-ID-炬住3DTlcttcopertics-Iterationrallyni.tTiijl.tiedFYLntTie柯.IGon

18、itryTineInformationPrantIVsler-B&uriduyLari-alilbEonriduryCondi3cnLFtdi!e-GrghisSoilForHblp3ptraeeFL第四步：迭代誤差限的選擇以模型默認的為主，如果軟件在運行的過程中出現(xiàn)計算不收斂的情況，可以通過適當調(diào)整誤差限來達到計算收斂的目的，操作如下圖所示：FflHTDRUS-ID-*I二1叵岡劇囁別加smt犧：）fre-pracqeel即菖|Fost-prtc4in.gtjflometryInTimeInformalInfotti-VflrifiL14BduiUpperOptimalIterationna

19、ngeWaterWaterWaterLaverTimeStepRultipLicB.liQnlaterCDiEntcdSoil珀口遼亡Intcorn-alLntdrpalbtinnJbblAE1CH3J6LiverLimitoftheTenicn1DD3近pgrLarrit4ftlwTaiLsinUpp&rTim4StepNultipLicB.tionForMaljo.prbte7L第五步：土壤水力傳導(dǎo)模塊選擇土壤水計算中普遍認可的且通用的wan方程作為本次計算的擬合方程，計算中不考慮水分運移滯后現(xiàn)象，操作如下圖所示：G二匸1!；!HXThIII-.-.I-r：-lI-il:-IIIITn.I

20、1-.7iiilIV”1iHYDKnS-lD-舷歯倔SoxlEtydraiiLicModelz/z/EIEE1E1B囪E10GeonctryIrinaInEcaPrintInfiWoiet-FLmIValsrTLWaterTLmHydraulicNeddSnuflePorositya-*cutGdnuchtari-Nualitn|7YilhAir-EnlrrValueqt廠HadifiAdvuiGanuchtanBroftks-Ccarftj-廣4sillagrDOEj.)1：iVaitfTlwF血/QTLdllsW1SoilTrof:IWISo:LProf:-Btl&I-T&rcsity/D

21、uai.-PenrseobilltyHodoLi-廣虹遼bQroaLlyd)uizuj-jdualvhilGsjnuchtm廣al-iorc-ityCyiabiJf_irmabilfe,vstergTtfi廣Du.41poirc-ELtyLc-irmbilc：!hitzd.gee二t廣Eu.4t；cifTTCabilltyiJdiL2T-s.ticvarcccno-ti&c廠DuiJp腫t陀曲訂ilyiBerkeandvan-?htec_K廣Lcok-apTabLcs-bb-stereEisr*Iohystr52EIfcsteresLEin.retejitiaicurvJtsteresLEin

22、.retejitiancufzeandjloei廣Lr：iticLLdryingcnCInitiall?11.ngcmForKitLppressTL第六步：包氣帶4種不同類型土壤巖性參數(shù)，可以通過輸入實驗室土壤顆粒分析結(jié)果獲得，操作如下圖所示：HnKUS-lD-月鍛試整場*口冋岡二tileVitvEre-prciCHSsinsQalculitionEtsullEtindwElpi制樹e亙同血iri枸訓(xùn)FrBrprwxTimg|FxLjnptjiaEjigTaterFlowParaacdzersSciL刁險iralHetvorkrtMiLctLon廠lenperatureDependenceIf

23、eLpatQrAlphan.I10.0780.433.61.150.249G0.520.0570.4551.410.1DH0.530.(0.4341.230.01GB0.540.065J/lbJ.b：lLiFarltslp.jtaeePL第七步：水流運移模塊上邊界條件據(jù)實際情況選擇隨時間變化的水頭邊界，下邊界選擇隨時間變化的水頭邊界，以變化的地下水位確定，初始條件選擇壓力水頭來確定，操作如下圖所示：rrt-procesEanaMainFro匚InEcornTineInfurnatioFrantZnErnatLWaterTlow-ItWaitfnfl-Sonffll&rTLom-SoiVoitf

24、riwfl-Bdfari*blpSonrida.SnLFrofile-rFrc*ieutaMaFLcccAtnoapheric:BCthSjrfactLferAinosptieTLICmilbSurfia.eeRunOJLavorSoimiisiryCanLtL3nant廠Conslant肓MiriitLeVwklLLcPressureHeadILwUpperBughryLaniLiiQonsiantPres:Cojiz-tantELozVriatLePressureTieadVariJJ-erressijreHe汕T*LaFracBiapDrain-agA廠eepaeeFaceKoriion

25、tiLDrliius第八步：輸入隨時間變化的上、下邊界壓力水頭值，操作如下圖所示：nHTDK0S-1D-1原隹試螫壕*dEIil4V:Fr4-pracacsln：ralcultionKbeuJIeIVind?*1HdlpTi.rhTpGVL丄11-0.0837381.01?II:4-11ll330.5921481.014Q-II4-4711ll&:0.6098531.01甘-0J730d6蘋770.3340631B；日-D.J07796Q.9Bu.J5:2.；1.Ul10I-J馬:.?11110217320.991212-O.J4O121a.99in13-1w:a.!i.14M0.026：.9

26、s15150.J93G39.95K1匚uijcqI1PravLOUEHglji.AddLiri4DdlgtaLin3DfoilltEFre-procesEijig一日IE1IE1ProfileInfornstLuaVaterFlow-Bomt-ldTyFIjhcsenilKe&dsSoil.draulLeFrprti45Rim7LmeInfnrTatinnVl4teLFlow-ITIVlmtfifFlov-HTOEUS-IDguideIVlilerFl-Be珂】ViriabLsEogdoiDoyou-3111torunPROFILE-3IipLlC&.tl4Tt?SoilProfiLeSoil

27、FrofLLe-ZZlat1:-0GE873:J-I12II)1113010.45824601111蚪.15-0.343Q1111E020.407333J111I6a.25-0.47670111210.30637501112aa.35-0.01233301112gIM-Hi;.ll;-1111”ID.450.732011121105-0Z3454G0111212.55-0.73709101113130c-073963G3-1EH0.65-07210201113150.7-07147270-13T第十二步：軟件運行的準備工作設(shè)置就緒，點擊0K軟件進行運行計算，如下圖所示：第十三步：下圖為軟件運行

28、完畢對話框，從下圖中可以看出不同時間的計算結(jié)果。-Hydras-IDCalculation.:原桂式醴場2制219,333342497,68E-ffi-b51E-92b66E-01-1.62JH0EPS219.6667325002E-fl3-.5BE-02.StKEKHJG5E-01-1.62測1.75220.02250215E-83-.49E-02.e0E*8064E-01-1.6200l.?5320.566042506-.28E-fl3-.51E-02h63E-01-1h54.991.75221.888B22S9S-,16E-fiG-.52E-02,8BE*fl0.61E-01-1,54.

29、BBK?5221.5SS042512-50E-B3-.54E-0260E-01-1.48測1.7532515-.33E8356E-03.60E*e053E-01-1H8.09l.?5333.888042519-57E-fi3-.61E-02.60E*flOE7E-01-i獨B01.75224_eeee42523-.8E-fl3-.tRE-02-eaE*eo5GE-B1-i.3fci.何225.008042527-76E-3-.76E-02.StKEKH54E-01-1.30測1.75226.042531-21E-83-.78E-02.ee(E*8053E-01-1.3600l.?5337.66

30、6042535.28E-04-.7BE-B2.60E*W52E-01-1.43.091.75228_888B42539J.8E-B3-.76E-02-saE+fla51E-01-1.51.R1何229.000O4254328E-fl3-.73E-02胴E*fl0.5BE-01-1.5脳1.75E3B.666032546.37E-8369E-03.60E*e043E-01-1.6?.09l.?5TimeItUJtCuniulaipEoTnpSullDOtEuBnthTnphflDothBot231.008032549.35E-B3-.66E-02.00E048E-01-1專001.75333.8

31、88022S51.3E-e9-WE偵.80E*8e47E-01-1.80.89l.?522553_38E-3-59E-B2-saE+fla.4?E-01-1.B7.R234.000022555-.55E-02胴E*fl046E-01-1.94脳1.75IfiBaltineIseclI.SQQQODDDD00PressEntertaGontiniiiE第十四步：不同埋深點包氣帶剖面壓力水頭隨時間的變化曲線，如下圖所示：PrintInforrifraterF譽-WaterJLawWaterFlor-WitcrFla廿VirifiLlaBdijtSol!FrafnleQefaulttrintLtgm

32、eiryIjlTinQInfarniKorizantalObservationNodes:PressureHeadsFsrKelj,.pressFt第十五步：不同埋深點包氣帶剖面含水率隨時間的變化曲線，如下圖所示：AEATtlKAbdsPrintLnfom込Flov-WaterFlov-NilerFl-fl&4WTiQtrylnfTintrInfdrn.1_nn站d加Flow-fViriatleBui|g|SoilFrofLLaWSoilPrcfLLe囲就HDrixonlsl7LITi4FPti-p-roaGsirigYeri匚alV*aterContentT-FileV：=H?LHefaoi

33、ltErint土-輕-:口I匸m巨二|CloseKim第十六步：包氣帶剖面壓力水頭變化曲線，如下圖所示PrintVla.lex壯tearWalerVl4.tcarFlov-ViriableB?ujSoilFroflLflSoilProfLLeesen.ilKe&ds/回回回EElalglIElKgw哥U-uriFrocesstGeOTetryInftTineInfortna-InfoniincarnwntilYtrtacalProfileInformation:PressureHeadlipfaultPrintfreri,otierideaForJfelp,pressFl第十七步：包氣帶剖面含

34、水率變化曲線，如下圖所示：atDkus-to-匚歸匡Prft-prflaesEingN-iiriFrocesst(reometryInf(lameInforri-iiPrintInfDmiVlilerFIoy-YlittiOfn.w-lrtilerFJ.q-W也七牡Flv-VariableBlubSoilFrofiUSoilProfLLe/IGJ回回EE0IE1IE1ForJfelp,pressFlH&riiantialProfile-InfdrTialioiVIktarCantmlYei-taaalDep-IIlIc二沙/r叵兇Prepfi：cessingQMo.LD.KotrizDn.tLVerticalTineZJActualSurfictPlmc討czzGeonetrzInfflTineInfarm釈FtlntInforricWaterFlor-ActualSurfa