新安江模型介绍：新安江模型介绍：三水源新安江模型蒸散发计算采用三层模型；产流计算采用蓄满产流模型；用自由水蓄水库结构将总径流划分为地表径流、壤中流和地下径流三种；流域汇流计算采用线性水库。模型结构：模型结构：模型计算：模型计算：在新安江模型中，流域蒸散发计算没有考虑流域内土壤含水量在面上分布的不均匀性，而是按土壤垂向分布的不均匀性将土层分为三层，用三层蒸散发模型计算蒸散发量。参数有流域平均张力水容量WM（mm），上层张力水容量UM(mm)，下层张力水容量LM(mm)，深层张力水容量DM(mm)，蒸散发折算系数KC和深层蒸散发扩散系数C。具体计算为若P+WU>=EP,则EU=EP,EL=0,ED=0;若P+WU<EP,则EU=P+WU;若WL>C*LM,则WL=(EP-EU)WL/LM,ED=0;若WL<C*LM且WL>=C*(EP-EU),则EL=C*(EP-EU),ED=0;若WL<C*LM且WL<C*(EP-EU),则EL=WL,ED=C*(EP-EU)-WL;水源划分中，本小组采用的是三水源划分。三水源用自由水蓄水库结构解决水源划分问题，自由水蓄水库结构考虑了包气带的垂向调蓄作用，按蓄满产流模型计算出总径流量R，先进入自由水蓄水库调蓄，再划分水源。模型参数调整：模型参数调整：1蒸散发能力折算系数KCKC是影响产流量计算最为重要和敏感的参数，产流计算中KC控制着水量平衡，因此，对水量计算是最重要的。KC主要反映流域平均高程与蒸发站高程之间差别的影响和蒸发皿蒸散发于陆面蒸散发间差别的影响。在实际模拟计算中KC值往往变化很大，最后需经模型调试并验证后确定。2流域平均张力水容量WM流域平均张力水容量WM表示流域干旱程度，分为UM,LM,DM。根据经验，南方湿润地区WM约为120~150mm，半湿润地区WM约为150~200mm。3流域蓄水容量―面积分布曲线指数B流域蓄水容量―B值反映划分单元流域张力水蓄水分布的不均匀程度。在一般情况下其取值与单元流域面积有关。在山丘地区，若单元流域面积较小，B=0.1左右；若单元流域面积中等，有几百到一千平方千米，B=0.2~0.3；若单元流域面积有几千平方千米，B=0.4左右。4不透水面积占全流域面积的比例IMIM值可由大比例尺的地形图通过GIS现代际醪饬砍隼础Ｔ谔烊涣饔颍琁M=0.01~0.02。5深层蒸散发扩散系数CC值主要取决于流域内深根植物的覆盖面积，南方多林地区C=0.15~0.20；在北方半湿润地区C=0.09~0.15。6自由水蓄水容量SMSM反映表土蓄水能力，其值受降雨资料时段均化的影响明显。当以日作为时段长时，在图层很薄的山区，其值为10mm或更小；在土深林茂透水性很强的流域，其值可取50mm或更大；一般流域在10~20mm。7自由水蓄水容量―面积分布曲线指数EX自由水蓄水容量―EX反映流域自由水蓄水分布的不均匀程度，一般EX=1.0~1.58自由水蓄水库对地下水和壤中流的日出流系数KG+KIKG的大小反映基岩和深层土壤的渗透性，KI的大小反映表层土的渗透性。9地下水消退系数CGCG可根据枯季地下径流退水规律来推求。若以日作为计算时段长，CG=0.950~0.99810壤中流消退系数CI若无壤中流CI=0，若壤中流丰富，则CI=0.9本小组经过程序的编写，用C++语言编写了新安江模型的计算界面，通过界面填写各种参数的数值调试计算结果，从而达到参数的简单率定。部分程序：部分程序：//三层蒸发模型的时段蒸发计算for(i=0;i<N;i++){W[i]=WU[i]+WL[i]+WD[i];if(WU[i]+P[i]>=EP[i]){EU[i]=EP[i];EL[i]=0;ED[i]=0;}else{EU[i]=WU[i]+P[i];if(WL[i]>=C*WLM){EL[i]=(EP[i]-EU[i])*WL[i]/WLM;ED[i]=0;}elseif(WL[i]>=C*(EP[i]-EU[i])){EL[i]=C*(EP[i]-EU[i]);ED[i]=0;}else{EL[i]=WL[i];ED[i]=C*(EP[i]-WU[i])-EL[i];}}//时段总蒸发量的计算E[i]=EU[i]+EL[i]+ED[i];//时段净雨量的计算PE[i]=P[i]-E[i];//时段初蓄水容量值的计算a[i]=WMM*(1.0-pow((1.0-W[i]/WM),(1.0/(1.0+B))));//时段产流量的计算if(PE[i]>0){if(a[i]+PE[i]<=WMM)R[i]=PE[i]+W[i]-WM+WM*pow