本发明公开了一种区域降雨人工模拟系统的降雨过程设计方法，该方法包括获取研究流域的空间地理信息，并将其按设定比例缩小得到的模型流域导入室内降雨区中；基于区域降雨人工模拟系统的空间剖分结构将室内降雨区划分成若干区块，并将模型流域位于的区块标记为流域子块；采用研究流域中雨量站测量的日降雨量及雨量站与流域子块的权重，计算流域子块的日降雨量；采用流域子块的日降雨量、雨量雨力关系模型、降雨历时及每个流域子块中模型流域的占比，计算模型流域的室外小时降雨量；计算模型流域在室外的净雨量；计算模型流域在室外的瞬时单位线；根据瞬时单位线构建室外净雨的时段单位线；之后根据时段单位线计算模型流域的室内净雨量。
1.一种区域降雨人工模拟系统的降雨过程设计方法，其特征在于，包括：获取研究流域的矢量信息，并将其按设定比例缩小得到的模型流域导入室内区域降雨人工模拟系统的覆盖区域中；
按照区域降雨人工模拟系统的空间剖分结构将模型流域剖分成若干区块，将模型流域对应的区块标记为流域子块；
采用研究流域中雨量站测量的日降雨量及雨量站与流域子块的权重，计算流域子块的日降雨量；
采用流域子块的日降雨量、雨量雨力关系模型、降雨历时及每个流域子块中模型流域的占比，计算模型流域的室外小时降雨量；
其中，I0为室外初损值；为室外平均下渗率；tc为产流历时；P后为室外降雨后期不产流的雨量；D为模型流域的室外小时降雨量；
根据模型流域的室外净雨量和室外径流量，计算模型流域在室外的瞬时单位线u(0，t)：其中，Q为室外径流量；n′为反应模型流域调蓄能力的参数；Γ(n′)为n′的伽马函数；k′为线性水库的调蓄系数；e为自然对数；t为时间变量； 分别为室外径流量Q的一阶和二阶原点矩； 分别为室外净雨量R的一阶和二阶原点矩；
根据模型流域在室外的瞬时单位线，构建模型流域在室外净雨为z mm的时段单位线：其中，q(Δt，t)为室外时段单位线；Δt为净雨时段；F为模型流域面积；
当单位时段内的小时净雨量为K个单位时，根据室外净雨的时段单位线计算模型流域第j时段末的室内净雨量R′j：q′＝K5*q(Δt，t)
其中， K5为q(Δt，t)与模型流域的室内时段单位线q′的换算比例尺，K5＜1；Q′i为i时段末的室内理论流量，i＝1，2，···，l为流量过程线时段数；q′i-j+1为i-j+1时段的室内时段单位线，i-j+1＝1，2，···n为单位线时段数；R′j为模型流域第j时段末的室内净雨量，j＝1，2，···，m为净雨时段数；
根据模型流域的模型流域第j时段末的室内净雨量R′j，计算模型流域的室内时段降雨量D′j：其中，I′0为室内初损值； 为室内平均下渗率；tc为产流历时；P′后为室内降雨后期不产流的雨量；
根据计算的所有室内时段降雨量，生成室内模型流域正上方人工降雨装置的降雨过程控制文件。
2.根据权利要求1所述的区域降雨人工模拟系统的降雨过程设计方法，其特征在于，所述流域子块的日降雨量的计算公式为：其中，Hj为第j个流域子块的日降雨量；wij为第j个流域子块到第i个雨量站的权重；pi为第i个雨量站的日降雨量；c为雨量站的总数量。
3.根据权利要求2所述的区域降雨人工模拟系统的降雨过程设计方法，其特征在于，所述wij的计算公式为：其中，rij为第j个流域子块到第i个雨量站的距离；b为权重指数；b＝0为平方法，b＝1为线性反比法，b＝2为RDS法。
4.根据权利要求1-3任一所述的区域降雨人工模拟系统的降雨过程设计方法，其特征在于，所述采用流域子块的日降雨量、雨量雨力关系模型、降雨历时及每个流域子块中模型流域的占比，计算模型流域的室外小时降雨量进一步包括：采用流域子块的日降雨量、雨量雨力关系模型、降雨历时，计算流域子块在降雨历时内的小时降雨量：其中，Aj为第j个流域子块在降雨历时内的小时降雨量；Hj为第j个流域子块的日降雨量；Sj为第j个流域子块的雨量雨力关系模型；a和b均为雨量雨力关系参数；ε为残差；T为降雨历时；n为暴雨衰减系数；
根据流域子块的小时降雨量和每个流域子块中模型流域的占比，计算模型流域的室外小时降雨量：其中，D为模型流域的室外小时降雨量；q为每个流域子块中模型流域的占比；B为流域子块的面积；Bj为第j个流域子块中模型流域的面积。
5.根据权利要求4所述的区域降雨人工模拟系统的降雨过程设计方法，其特征在于，所述计算流域子块在降雨历时内的小时降雨量进一步包括：根据研究流域的日降雨量，构建雨量雨力关系模型：
根据计算的雨力Sj、流域子块的日降雨量Hj和率定的衰减系数n，计算一天内的降雨历时T：根据雨力Sj和降雨历时T，流域子块在降雨历时内的小时降雨量为：
6.根据权利要求1所述的区域降雨人工模拟系统的降雨过程设计方法，其特征在于，所述 和 的计算公式分别为：其中，Rl为第l时段的模型流域的室外净雨量；Ql为第l时段的模型流域的室外径流量；
7.根据权利要求1-3、5、6任一所述的区域降雨人工模拟系统的降雨过程设计方法，其特征在于，所述研究流域的矢量信息的生成方法包括：获取研究流域的地形数据和雨量站的位置，并将地形数据和雨量站的位置导入ArcGis中生成研究流域的矢量信息。
一种区域降雨人工模拟系统的降雨过程设计方法
技术领域
[0001] 本发明涉及水文水资源领域，具体涉及一种区域降雨人工模拟系统的降雨过程设计方法。
背景技术
[0002] 区域水文过程试验是当前水文试验研究的主要发展方向之一，其中如何根据试验研究目标科学设定室内区域水文过程试验的人工降雨过程是其关键问题之一。当前的人工降雨试验主要利用点尺度的降雨观测资料，选定典型降雨雨强和降雨历时，开展人工降雨试验，实验过程强调降雨区域内的均匀度指标。然而由于实际降雨过程存在时间和空间变化，区域水文过程试验的研究对象强调指定区域内降雨、地形、径流等过程的相似性、一致性。由于两者之间的关注点不一致，从而致使现有的室内水文试验中对人工降雨过程的设定方法难以用于区域水文过程试验研究中的人工降雨过程设计。
发明内容
[0003] 针对现有技术难以实现区域水文过程实验中难以模拟降雨时空变异性等不足，本发明提供了一种能够模拟降雨时空变异性的区域降雨人工模拟系统的降雨过程设计方法。
[0004] 为了达到上述发明目的，本发明采用的技术方案为：
[0005] 提供一种区域降雨人工模拟系统的降雨过程设计方法，其包括：
[0006] 获取研究流域的矢量信息，并将其按设定比例缩小得到的模型流域导入室内降雨区中；
[0007] 按照区域降雨人工模拟系统的空间剖分结构将模型流域剖分成若干区块，将模型流域对应的区块标记为流域子块；
[0008] 采用研究流域中雨量站测量的日降雨量及雨量站与流域子块的权重，计算流域子块的日降雨量；
[0009] 采用流域子块的日降雨量、雨量雨力关系模型、降雨历时及每个流域子块中模型流域的占比，计算模型流域的室外小时降雨量；
[0010] 计算模型流域在室外的净雨量R：
[0011]
[0012] 其中，I0为室外初损值； 为室外平均下渗率；tc为产流历时；P后为室外降雨后期不产流的雨量；D为模型流域的室外小时降雨量；
[0013] 根据模型流域的室外净雨量和室外径流量，计算模型流域在室外的瞬时单位线u(0，t)：
[0014]
[0015]
[0016] 其中，Q为室外径流量；n′为反应模型流域调蓄能力的参数；Γ(n′)为n′的伽马函数；k′为线性水库的调蓄系数；e为自然对数；t为时间变量； 分别为室外径流量Q的一阶和二阶原点矩；MR(1)、MR(2)分别为室外净雨量R的一阶和二阶原点矩；
[0017] 根据模型流域在室外的瞬时单位线，构建室外净雨为z mm的时段单位线：
[0018]
[0019]
[0020] 其中，q(Δt，t)为室外时段单位线；Δt为净雨时段；F为模型流域面积；
[0021] 当单位时段内的小时净雨量为K个单位时，根据室外净雨的时段单位线计算模型流域第j时段末的室内净雨量Rj′：
[0022]
[0023] 其中， K5为q(Δt，t)与模型流域的室内时段单位线q′的换算比例尺，K5＜1；Q′i为i时段末的室内理论流量，i＝1，2，···，l为流量过程线时段数；q′i-j+1为i-j+1时段的室内时段单位线，i-j+1＝1，2，···n为单位线时段数；R′j为模型流域第j时段末的室内净雨量，j＝1，2，···，m为净雨时段数；
[0024] 根据模型流域的模型流域第j时段末的室内净雨量R′j，计算模型流域的室内时段降雨量D′j：
[0025]
[0026] 其中，I′0为室内初损值； 为室内平均下渗率；tc为产流历时；P′后为室内降雨后期不产流的雨量；
[0027] 根据计算的所有室内时段降雨量，生成室内模型流域正上方人工降雨装置的降雨过程控制文件。
[0028] 进一步地，所述流域子块的日降雨量的计算公式为：
[0029]
[0030] 其中，Hj为第j个流域子块的日降雨量；wij为第j个流域子块到第i个雨量站的权重；pi为第i个雨量站的日降雨量；c为雨量站的总数量。
[0031] 进一步地，所述wij的计算公式为：
[0032]
[0033] 其中，rij为第j个流域子块到第i个雨量站的距离；b为权重指数；b＝0为平方法，b＝1为线性反比法，b＝2为RDS法。
[0034] 进一步地，所述采用流域子块的日降雨量、雨量雨力关系模型、降雨历时及每个流域子块中模型流域的占比，计算模型流域的室外小时降雨量进一步包括：
[0035] 采用流域子块的日降雨量、雨量雨力关系模型、降雨历时，计算流域子块在降雨历时内的小时降雨量：
[0036]
[0037] 其中，Aj为第j个流域子块在降雨历时内的小时降雨量；Hj为第j个流域子块的日降雨量；Sj为第j个流域子块的雨量雨力关系模型；a和b均为雨量雨力关系参数；ε为残差；T为降雨历时；n为暴雨衰减系数；
[0038] 根据流域子块的小时降雨量和每个流域子块中模型流域的占比，计算模型流域的室外小时降雨量：
[0039]
[0040] 其中，D为模型流域的室外小时降雨量；q为每个流域子块中模型流域的占比；B为流域子块的面积；Bj为第j个流域子块中模型流域的面积。
[0041] 进一步地，所述计算流域子块在降雨历时内的小时降雨量进一步包括：
[0042] 根据研究流域的日降雨量，构建雨量雨力关系模型：
[0043] Sj＝aHj+b+ε；
[0044] 根据计算的雨力Sj、流域子块的日降雨量Hj和率定的衰减系数n，计算一天内的降雨历时T：
[0045]
[0046] 根据雨力Sj和降雨历时T，流域子块在降雨历时内的小时降雨量为：
[0047]
[0048] 进一步地，所述 MR(1)和MR(2)的计算公式分别为：
[0049]
[0050]
[0051] 其中，Rl为第l时段的模型流域的室外净雨量；Ql为第l时段的模型流域的室外径流量； l＝1，2···，n-1；Δt为净雨时段。
[0052] 进一步地，所述研究流域的矢量信息的生成方法包括：
[0053] 获取研究流域的地形数据和雨量站的位置，并将地形数据和雨量站的位置导入ArcGis中生成研究流域的矢量信息。
[0054] 本发明的有益效果为：本方案利用地理信息系统进行研究流域和模型流域的空间对应分析，基于区域降雨人工模拟系统的空间剖分结构，通过构建研究流域和模型流域的时段单位线，实现区域降雨人工模拟系统的降雨过程设计；采用该方法可以有效满足室内区域水文过程试验研究对模拟降雨过程的设计要求，实现室内水文过程较好的拟合野外研究区域，具有较高的可靠度、准确性和实用性。
附图说明
[0055] 图1为研究流域在室内区域降雨量的预测方法的流程图。
[0056] 图2为模型流域导入室内降雨区的示意图。
[0057] 图3为S曲线的示意图。
具体实施方式
[0058] 下面对本发明的具体实施方式进行描述，以便于本技术领域的技术人员理解本发明，但应该清楚，本发明不限于具体实施方式的范围，对本技术领域的普通技术人员来讲，只要各种变化在所附的权利要求限定和确定的本发明的精神和范围内，这些变化是显而易见的，一切利用本发明构思的发明创造均在保护之列。
[0059] 参考图1，图1示出了研究流域在室内区域降雨量的预测方法的流程图；如图1所示，该方法100包括步骤101至步骤110。
[0060] 在步骤101中，获取研究流域的矢量信息，并将其按设定比例缩小得到的模型流域导入室内区域降雨人工模拟系统的覆盖区域中，其中，区域降雨人工模拟系统为：有效降雨面积26m×40m，共分为110个独立控制单元(其中共分为10个大区，每个大区又分为11个小区)；每个降雨单元雨强变化范围10-200mm/h；单元内均匀度大于0.8；雨强变化频率可以做到5分钟调整一次；通过输入降雨控制文件，可以实现人工降雨。
[0061] 如图2所示，40×26的区块为室内降雨区，其中位于室内降雨区中的弧形区域为模型流域。在进行缩小时，必须保证研究流域的矢量信息在x、y轴上的长宽比例不变的情况下对研究流域的矢量信息进行缩放得到模型流域，模型流域的水平比例尺为K1，垂直比例尺为K2。
[0062] 其中设定比例的选择主要跟室内降雨区的面积相关，在缩放过程中，为了更好地在室内模拟野外流域降雨过程，优选缩小后的模型流域尽量填充整个室内降雨区。
[0063] 实施时，本方案优选研究流域的矢量信息的生成方法包括：获取研究流域的地形数据和雨量站的位置，并将地形数据和雨量站的位置导入ArcGis中生成研究流域的矢量信息。
[0064] 在步骤102中，按照区域降雨人工模拟系统的空间剖分结构将模型流域剖分成若干区块，将模型流域对应的区块标记为流域子块；实施时，本方案优选将室内降雨区共划分成110个区块，找到模型流域所在区域，并将有模型流域位于的区块标记为流域子块，并将所有的流域子块依次编上序号。
[0065] 在步骤103中，采用研究流域中雨量站测量的日降雨量及雨量站与流域子块的权重，计算流域子块的日降雨量；具体地，流域子块的日降雨量的计算公式为：
[0066]
[0067] 其中，Hj为第j个流域子块的日降雨量；wij为第j个流域子块到第i个雨量站的权重；pi为第i个雨量站的日降雨量；c为雨量站的总数量。
[0068] 所述wij的计算公式为：
[0069]
[0070] 其中，rij为第j个流域子块到第i个雨量站的距离；b为权重指数；b＝0为平方法，b＝1为线性反比法，b＝2为RDS法。
[0071] 在步骤104中，采用流域子块的日降雨量、雨量雨力关系模型、降雨历时及每个流域子块中模型流域的占比，计算模型流域的室外小时降雨量；
[0072] 在本发明的一个实施例中，所述采用流域子块的日降雨量、雨量雨力关系模型、降雨历时及每个流域子块中模型流域的占比，计算模型流域的室外小时降雨量进一步包括步骤201和步骤202：
[0073] 在步骤201中，采用流域子块的日降雨量、雨量雨力关系模型、降雨历时，计算流域子块在降雨历时内的小时降雨量：
[0074]
[0075] 其中，Aj为第j个流域子块在降雨历时内的小时降雨量；Hj为第j个流域子块的日降雨量；Sj为第j个流域子块的雨量雨力关系模型；a和b均为雨量雨力关系参数；ε为残差；T为降雨历时；n为暴雨衰减系数。
[0076] 在步骤201中计算流域子块在降雨历时内的小时降雨量的方法包括：
[0077] 首先通过对研究流域的日降雨量进行分析，构建雨量雨力关系模型：
[0078] Sj＝aHj+b+ε
[0079] 由于雨量雨力关系模型中的流域子块的日降雨量Hj可以在103步骤中求解出来，则本方案通过构建的雨量雨力关系模型就可以直接计算出第j个流域子块的雨力Sj。
[0080] 之后，根据计算的雨力Sj、流域子块的日降雨量Hj和率定的衰减系数n，计算一天内的降雨历时T：
[0081]
[0082] 之后，根据雨力Sj和降雨历时T，流域子块在降雨历时内的小时降雨量为：
[0083]
[0084] 在步骤202中，根据流域子块的小时降雨量和每个流域子块中模型流域的占比，计算模型流域的室外小时降雨量：
[0085]
[0086] 其中，D为模型流域的室外小时降雨量；q为每个流域子块中模型流域的占比；B为流域子块的面积；Bj为第j个流域子块中模型流域的面积。
[0087] 在步骤105中，计算模型流域在室外的净雨量R：
[0088]
[0089] 其中，I0为室外初损值； 为室外平均下渗率；tc为产流历时；P后为室外降雨后期不产流的雨量；D为模型流域的室外小时降雨量。
[0090] 在步骤106中，根据模型流域的室外净雨量和室外径流量，计算模型流域在室外的瞬时单位线u(0，t)：
[0091]
[0092]
[0093] 其中，Q为室外径流量；n′为反应模型流域调蓄能力的参数；Γ(n′)为n′的伽马函数；k′为线性水库的调蓄系数；e为自然对数；t为时间变量； 分别为室外径(1) (2)流量Q的一阶和二阶原点矩；MR 、MR 分别为室外净雨量R的一阶和二阶原点矩。
[0094] 实施时，本方案优选所述 MR(1)和MR(2)的计算公式分别为：
[0095]
[0096]
[0097] 其中，Rl为第l时段的模型流域的室外净雨量；Ql为第l时段的模型流域的室外径流量； l＝1，2···，n-1；Δt为净雨时段。
[0098] 在步骤107中，根据模型流域在室外的瞬时单位线，构建模型流域在室外净雨为z mm的时段单位线：
[0099]
[0100]
[0101] 其中，q(Δt，t)为室外时段单位线；Δt为净雨时段；F为模型流域面积。
[0102] 在本发明的一个实施例中，根据模型流域在室外的瞬时单位线，构建模型流域在室外净雨为z mm的时段单位线进一步包括：
[0103] 首先对得到的模型流域在室外的瞬时单位线进行求导，得到S曲线：
[0104]
[0105] 由于n′、k′已知，以不同的t代入上式积分，就可以得到如图3所示的S曲线。以t＝0为起点的S(t)曲线向后平移一个Δt时段，即可得到S(t-Δt)曲线，两条曲线的纵坐标差：
[0106] U(Δt，t)＝S(t)-S(t-Δt)
[0107] 上式U(Δt，t)即为时段为Δt的无因次时段单位线，将无因次时段单位线换算成时段为Δt，净雨为zmm的时段单位线为：
[0108]
[0109] 其中，q(Δt，t)为室外时段单位线；Δt为净雨时段；F为模型流域面积。
[0110] 在步骤108中，当单位时段内的小时净雨量为K个单位时，根据室外净雨的时段单位线计算模型流域第j时段末的室内净雨量R′j：
[0111]
[0112] 其中， K5为q(Δt，t)与模型流域的室内时段单位线q′的换算比例尺，K5＜1；Q′i为i时段末的室内理论流量，i＝1，2，···，l为流量过程线时段数；q′i-j+1为i-j+1时段的室内时段单位线，i-j+1＝1，2，···n为单位线时段数；R′j为模型流域第j时段末的室内净雨量，j＝1，2，···，m为净雨时段数；
[0113] 当i＝1的时候，j＝1；i＝2的时候，是j＝1和j＝2时两个结果的和；i＝3时，就是j＝1、j＝2、j＝3这三个的和。
[0114] 在步骤109中，根据模型流域的模型流域第j时段末的室内净雨量R′j，计算模型流域的室内时段降雨量Dj′：
[0115]
[0116] 其中，I′0为室内初损值； 为室内平均下渗率；tc为产流历时；P′后为室内降雨后期不产流的雨量。
[0117] 在步骤110中，根据计算的所有室内时段降雨量，生成室内模型流域正上方人工降雨装置的降雨过程控制文件。
[0118] 其中的控制文件为每个流域子块的人工降雨量，其具体获取方法为：按照分区情况将降雨(所有室内时段降雨量)分配到每一个区块上，利用面积权重来分配得到每个区域的降雨，并采用每一个降雨插值分区的室内降雨量实施人工降雨。
[0119] 下面采用Nash效率、相关系数和相对误差对本方案构建的时段单位线的效果进行说明：
[0120] 将生成的室内模型流域正上方人工降雨装置的降雨过程控制文件，输入到人工降雨控制平台中实施降雨，并对降雨过程中的雨强和指定观测点的流量过程进行监测，记录完整的数据。
[0121] 采用以下计算公式计算模型流域的室内理论流量：
[0122]
[0123] 其中，Q′理i为i时段末的室内理论流量，i＝1，2，···，l为流量过程线时段数；q′i-j+1为i-j+1时段的室内时段单位线，i-j+1＝1，2，···n为单位线时段数；R′j为模型流域第j时段末的室内净雨量，j＝1，2，···，m为净雨时段数；
[0124] 当i＝1的时候，j＝1；i＝2的时候，是j＝1和j＝2时两个结果的和；i＝3时，就是j＝1、j＝2、j＝3这三个的和。
[0125] 下面采用Nash效率、相关系数和相对误差对Q′理i与实验得到实测流量过程线Q′n进行校验：
[0126] Nash效率的计算公式为：
[0127]
[0128] 其中，Qi为采用本方案的方法计算得到的降雨量，Qi为上面计算的Q′理i，qi为室内实测降雨量，为室内实测降雨量的平均值。
[0129] 相关系数的计算公式为：
[0130]
[0131] rxy为相关系数；n为系列的样本数；X、Y分别代表实测系列和模拟系列的数值。
[0132] 相对误差的计算公式为：
[0133]
[0134] 其中，Dv为模相对误差(％)；F0为室内实测降雨量的均值；R为采用本方案的方法计算得到的降雨量均值。
[0135] 将在模拟人工降雨过程中记录的雨强和指定观测点的流量带入Nash效率、相关系数和相对误差的计算公式中，可以得到Nash效率接近于1，相关系数rxy位于0.8-1.0之间；相对误差绝对值接近于零。
[0136] 通过上述对对Q′理与实验得到实测流量过程线Q′n的校验可以看出，本方案得到的区域降雨人工模拟系统的降雨过程设计方法能够到到试验效果。
