[发明专利]一种全分布式考虑流域内水库影响的汇流模拟方法

权利要求书

1.一种全分布式考虑流域内水库影响的汇流模拟方法，其特征在于，用于获取目标流域内各栅格单元在目标时刻t时的泄流量，方法包括如下步骤：

步骤1：基于目标流域的数字高程数据，获取目标流域中各栅格单元的坡度、流向、汇流累计值；基于各栅格单元的流向和汇流累计值，计算各栅格单元的演算次序；然后进入步骤2；

计算演算次序的方法包括如下步骤：

步骤1.1：分别针对汇流累计值等于1的各栅格单元，将其演算次序赋值为1；分别针对汇流累计值不等于1的各栅格单元，初始化其演算次序n＝0；然后进入步骤1.2；

步骤1.2：分别针对各演算次序n＝0的栅格单元，将其作为待处理单元；同时分别针对各待处理单元，执行如下操作，然后进入步骤1.3；

根据与待处理单元相邻各栅格单元的流向，获取水流流向待处理单元的各相邻栅格单元：

若水流流向待处理单元的各栅格单元的演算次序均不等于零，则对水流流向待处理单元的各相邻栅格单元中的演算次序最大值加1获取一累计值，将该累计值作为待处理单元的演算次序；

若至少一个水流流向待处理单元的相邻栅格单元的演算次序等于零，则不改变待处理单元的演算次序；

步骤1.3：判断是否存在演算次序n＝0的栅格单元，是则返回步骤1.2；否则停止循环；

步骤2：根据目标流域内各水库的位置，将各栅格单元区分为水库栅格和非水库栅格；然后进入步骤3；

步骤3：分别针对各水库栅格i，将其作为第一待处理栅格，依次执行步骤3.1.1至步骤3.1.2：

步骤3.1.1：判断第一待处理栅格i在目标时刻t时的库容ResStroage_i，t是否大于或等于其总库容是则根据如下公式获取第一待处理栅格i在目标时刻t时的泄流量ResDrainage_i，t；否则进入步骤3.1.2；

ResDrainage_i，t＝(TotalStorage_i-ComputeStorage_i)*RDC+(ResStroage_i，t-TotalStorage_i)

其中：ComputeStorage_i为第一待处理栅格i的运行库容；RDC为目标流域的水库泄流系数；

步骤3.1.2：判断第一待处理栅格i在目标时刻t时的库容ResStroage_i，t是否大于或等于其运行库容ComputeStorage_i，是则跟据如下公式获取第一待处理栅格i在目标时刻t时的泄流量ResDrainage_i，t，否则取第一待处理栅格i在目标时刻t时的泄流量ResDrainage_i，t＝0；

ResDrainage_i，t＝(ResStroage_i，t-ComputeStorage_i)*RDC

通过执行步骤3.1.1至3.1.2，获取第一待处理栅格对应目标时刻t的泄流量，进而获取各水库栅格分别对应目标时刻t的泄流量；

分别针对演算次序等于1的各非水库栅格，执行如下操作：根据产流模型，获取各非水库栅格在目标时刻t时的产流量，将该产流量作为对应非水库栅格的泄流量；

针对演算次序不等于1的各非水库栅格，按照演算次序由小到大的顺序依次获取各非水库栅格的泄流量，获取对应非水库栅格的泄流量的方法包括如下步骤：

分别针对演算次序不等于1的各非水库栅格，将其作为第二待处理栅格i'，根据如下公式：

O_i′，t＝C₁I_i′，t+C₂I_i′，t-1+C₃O_i′，t-1

获取第二待处理栅格目标时刻t的泄流量O_i′，t，进而获取各非水库栅格分别对应目标时刻t的泄流量；

其中，C₁、C₂、C₃为和第二待处理栅格i'的参考流量以及流速相关的系数，I_i'，t为第二待处理栅格i'目标时刻t时的入流量，I_i'，t-1为第二待处理栅格i'时刻t-1时的入流量，O_i′，t-1为第二待处理栅格i'时刻t-1时的泄流量。