[发明专利]一种城镇化下平原区分布式产流计算方法

权利要求书

1.一种城镇化下平原区分布式产流计算方法，其特征在于，方法步骤如下：

S1：采用ArcGIS软件划分子产水单元；

S2：采用ArcGIS软件统计每个子产水单元中不同土地利用类型的面积；

S3：采用MATLAB软件计算每个子产水单元中不同土地利用类型的产流量。

2.根据权利要求1所述的一种城镇化下平原区分布式产流计算方法，其特征在于，所述S1中，通过ArcGIS将降雨站点和区域边界进行叠加。

3.根据权利要求2所述的一种城镇化下平原区分布式产流计算方法，其特征在于，所述降雨站点通过泰森多边形确认区域的划分，所述区域边界为子流域边界或行政区边界。

4.根据权利要求1所述的一种城镇化下平原区分布式产流计算方法，其特征在于，所述S2中土地利用类型包括水面、水田、旱地与非耕地、及城镇建设用地。

5.根据权利要求1所述的一种城镇化下平原区分布式产流计算方法，其特征在于，所述S3中MATLAB的计算所需数据包括降雨量数据、土地利用类型数据、系数。

6.根据权利要求5所述的一种城镇化下平原区分布式产流计算方法，其特征在于，所述系数包括水面蒸发系数、陆地蒸发系数、水稻需水系数、蒸发量、月降雨日数。

7.根据权利要求1所述的一种城镇化下平原区分布式产流计算方法，其特征在于，所述S3中MATLAB的计算结果包括各子产水单元月产水量，各区域、各类土地利用类型年产水量。

8.根据权利要求1-7任一项所述的一种城镇化下平原区分布式产流计算方法，其特征在于，对于产流量计算，采用月产流模型方法进行计算，首先分别计算各产水单元内每种土地利用类型产流，然后进行汇总，从而计算总的产流量。

a.水面产流：水面产流为降水量与蒸发量之差，即：

R_W＝P-C_E×E

式中，R_W为时段水面产水量；P为时段降水量；C_E为蒸发皿折算系数；E为蒸发量。

b.水田产流：时段水田产水量为：

R_R＝P-α×C_E×E×D_P

式中α为水稻需水系数，D_P为月内降水日数比例。

c.旱地与非耕地产流：采用蓄满产流模型计算旱地与非耕地产流计算模型，其计算式为：

EE=CK×E×WDWM]]>

W_MM＝W_M×(1+B)

A=WMM×(1-(1-WDWM)11+B)]]>

式中，C_K为陆面蒸发折算系数；W_D为土壤初始含水量；W_M为平均蓄水量(土层最大可能缺水量)；E_E为旱地蒸发量；W_MM为蓄水容量曲线的最大值；B为蓄水量曲线指数；R_D为旱地与非耕地产水量。

d.城镇建设用地产流：建设用地的透水性差，其产流可以表示为：

R_C＝C_C×P

式中，R_C为时段(月)建设用地产水量；C_C为城镇径流系数；P为时段(月)降水量。

e.总产水量计算：各分区的总产流深为各下垫面产流之和：

R_T＝R_W+R_R+R_D+R_C

式中，R_T为分区的总产水量。