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

1、天山北麓平原区包气带水分运移机理与数值分析HYDRUS1D软件使用手册目录一、软件介绍二、安装指南三、使用指南四、实例介绍 一、软件介绍HYDRUS1D是国际地下水模型中心公布的，计算包气带水分、盐分运移规律的软件，用它可以解算在不同边界条件制约下的数学模型。若将坐标原点选在地面，取z轴向下为正，则一维饱和非饱和带水分运移基本方程为：理丄|k（6）（也-1）1-Sdtdzdz_式中：e为含水率；h为负压水头；S为植物根系吸水量，对裸露区为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|二）数据输入及运行1、选择计算模块2、模型几何结构属性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、选择输出结果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、岩性参数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、上、下边界条件数据输入10、

10、几何结构确定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、模型运行（三）运行结果1、剖面水头、含水率等计算结果2、边界流量信息LDHORUS-1Df

11、ileYitmriBaimdaryfatcrTinafjninEVacseebG-bondtryInft1LmeLdde泌FriutInforncVaterFl州-nterKL&iv-VatarFldiv-&terFI&iy-YuriihlaBdhj3、参数属性HTDRUS-lIii-IFilmViwTyi口圍風亘j-FlE-pTQEZULLngMildPracAE-EtGesneiryInee沾aridHeais7ZZ/hbeffibefiriorizontaJ.VarticalLl_二=rHa九fTDC45ZGGfonetryInfcliniEInEornPrintInEornsiVfii

12、erFLw-RffllarFLw-WalerTLom-IVqIsrTLom-arLatleBotitgoalfrofilgSoiLFrofiliRun-TimeTineTTin亡Step*vrinaSlapfelmbsifiItdrtiTis間帘帅电羁權冊卿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、实例介绍对奎屯原位试验场包气带裸露区非均质水分运移规律进行模拟。剖面岩性结构从上到下共7层；包括：壤土、粉壤土、粉质黏壤土以及砂质壤土4种岩性类型；水位埋深最大不超过9m；模型运行期为4月1号至11月20号，共234天。选择万方程作为模型计算的方程；上边界为变化的负压边界，下边界为变化的水头边界；初始条件以4月1号观测的负压为准输入模型进行计算；输出每月最后一天的结果信息。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第一步：软件可以进行计算的模块包括水分运移模块、溶质运移模块、热传导模块、植物根系吸收水分模块以及植

15、物根系生长模块等，本次模拟只对包气带水分运移规律进行模拟，故选择水流模块，操作如下图所示：HHmtUS-lD-世检试豔璋*LinlAHiE订匕VliErTrocesiiLEQd-ciilatiorifiesyiltsHindoo心眉圜釦%|S|ffl|?|Jfrfrir&cessinr第二步：选定计算长度单位m,剖面岩性类型4种，从上到下共分为7层,计算深度为9m,操作如下图所示:第二步：选定计算时间单位Days，计算时间从4月1号至11月20号，共234天,边界条件随时间变化记录的最小间隔为1天,故随时间变化的边界条件记录数为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创uDdaryCondjtio厂feiscroLuEi匚辺Data.厂DailyVarjationsBoringTisyS-nerAteibyI匕HYDKUS-ID-爺也姑鮭场*Z回冈JlilaV:sivTrfl-pracqeelngLalcvl.tioilKarultzIViniiA,Hdlpq|Q|知3,sm?Fte-口rue亡三ring第三步：要求输出每月最后一天的计算结果，统计共有9天的计算结果，操作如下图所示：HYDRUS-ID-炬住3DTlcttcopertics-Iterationrallyni.tTiijl.tiedFYLntTie柯.IGon

18、itryTineInformationPrantIVsler-B&uriduyLari-alilbEonriduryCondi3cnLFtdi!e-GrghisSoilForHblp3ptraeeFL第四步：迭代误差限的选择以模型默认的为主，如果软件在运行的过程中出现计算不收敛的情况，可以通过适当调整误差限来达到计算收敛的目的，操作如下图所示：FflHTDRUS-ID-*I二1叵冈剧囁別加smt牺：）fre-pracqeel即菖|Fost-prtc4in.gtjflometryInTimeInformalInfotti-VflrifiL14BduiUpperOptimalIterationna

19、ngeWaterWaterWaterLaverTimeStepRultipLicB.liQnlaterCDiEntcdSoil珀口辽亡Intcorn-alLntdrpalbtinnJbblAE1CH3J6LiverLimitoftheTenicn1DD3近pgrLarrit4ftlwTaiLsinUpp&rTim4StepNultipLicB.tionForMaljo.prbte7L第五步：土壤水力传导模块选择土壤水计算中普遍认可的且通用的wan方程作为本次计算的拟合方程，计算中不考虑水分运移滞后现象，操作如下图所示：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种不同类型土壤岩性参数，可以通过输入实验室土壤颗粒分析结果获得，操作如下图所示：HnKUS-lD-月锻试整场*口冋冈二tileVitvEre-prciCHSsinsQalculitionEtsullEtindwElpi制樹e亘同血iri枸训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第七步：水流运移模块上边界条件据实际情况选择随时间变化的水头边界，下边界选择随时间变化的水头边界，以变化的地下水位确定，初始条件选择压力水头来确定，操作如下图所示：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第十二步：软件运行的准备工作设置就绪，点击0K软件进行运行计算，如下图所示：第十三步：下图为软件运行

28、完毕对话框，从下图中可以看出不同时间的计算结果。-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