一种基于DRAINMOD数值模拟暗管排水系统布置最优参数确定方法

一种基于drainmod数值模拟暗管排水系统布置最优参数确定方法
技术领域
1.本发明涉及农田暗管布设技术领域，具体是涉及一种基于drainmod数值模拟暗管排水系统布置最优参数确定方法。


背景技术：

2.近年来，农村规模化种植的发展以及机械化水平的提高，与现有农田水利田间工程配套不相适应；加之气象条件多变，农田排水不畅导致作物生长季节渍害严重以及收获季节农业机械无法下田作业的情况，继而造成粮食欠收的不良现象。农田末级排水明沟承担农田排涝、降渍任务，田块长度均取决于末级排水明沟的间距，现有的明沟排水系统不利于机械化作业，阻碍了农业现代化发展的步伐。针对目前农业规模化、机械化和现代化发展要求，采用农田暗管降渍，将排涝、降渍系统分开，可为满足农田灌溉、排涝、降渍、农机作业等综合要求的高标准农田建设创造条件，以便形成连片土地，提高农机作业效率，进而促进农业生产。开展面向农业机械化的田间暗管排水工程模式研究，确定满足现代农业发展需要的暗管管径、埋深、间距、土地平整度等参数很有必要，为了确定合适的暗管埋深、间距、半径等参数，保证排水效果，需要开发更为高效和准确的计算机模拟方法。


技术实现要素：

3.发明目的：本发明提供一种基于drainmod数值模拟暗管排水系统布置最优参数确定方法，能够高效和准确确定暗管排水系统布置的最优参数。
4.技术方案：为解决上述问题，本发明提供一种基于drainmod数值模拟暗管排水系统布置最优参数确定方法，包括以下步骤：
5.(1)模型构建：设计暴雨，以暴雨雨后n至n+1天地下水降至田面以下最小埋深为目标，建立如下目标函数：
6.h＝min hw(d,h,p,s,t,f,γ,h1,h2,v,m,σ,r)
ꢀꢀꢀ
(1)
7.式中，h为地下水最小埋深；hw为地下水埋深；d为暗管管径；h为暗管埋深；p为暗管间距；s为地表平整度；t为排水控制时段；f为土壤吸力；γ为土壤体积含水率；h1为最大地表蓄水深度；h2为kirkham积水深度；v为土壤可排空体积；m为潜水上升通量；σ为排水模数；r为暗管有效半径；
8.(2)设置目标函数求解的约束条件；
9.(3)采用大系统试验选优方法，结合约束条件对目标函数求解，具体包括以下步骤：
10.(3.1)数据准备：测定田块面积a、田面高程z1、农沟水位z2；按照暴雨统计资料，推算设计暴雨过程；测定当地土壤水分特征参数值；设定排水系统特性参数；测定当地旱作物生长的盛期要求的地下水埋深数值；测定当地农沟水面离地面距离；
11.(3.2)根据暗管管径d、暗管埋设深度h、暗管间距p、地表平整度s四种因素，每种因
素分别设置为j种水平，参照大系统试验选优理论构建4因素j水平正交表；
12.(3.3)对正交表中每一行试验方案，结合约束条件，采用drainmod模型进行数值模拟；
13.(3.4)根据步骤(3.3)的模拟结果构造极差表，通过正交分析，获得理论最优暗管管径d
*
、暗管埋设深度h
*
、暗管间距p
*
、地表平整度s
*
；
14.(3.5)将步骤(3.4)获得的理论最优参数组合，代入目标函数并结合约束条件，通过drainmod模型重新进行数值模拟，获得地下水最小埋深h。
15.进一步的，步骤(2)约束条件包括：
16.(2.1)地下水埋深约束：
17.h≥0.8
ꢀꢀꢀ
(2)
18.(2.2)排水控制时段约束：
19.n≤t≤n+1
ꢀꢀꢀ
(3)
20.(2.3)暗管埋深约束：
21.h≥hwꢀꢀꢀ
(4)
22.h≤h
农-0.1
ꢀꢀꢀ
(5)
23.式中，h
农
为农沟水面离地面距离，hw为旱作物生长的盛期要求的地下水埋深；
24.进一步的，步骤(3.3)包括具体以下步骤：
25.(3.3.1)设置drainmod模型的输入参数，输入参数包括土壤水分特征参数、土壤排水特性参数、排水系统特征参数；
26.(3.3.2)对于正交表中每一个方案，结合地下水埋深约束、排水控制时段约束、暗管埋深约束条件，采用drainmod模型，分别模拟暴雨雨后n至n+1天土壤地下水埋深变化过程，获得控制排水时段t内地下水埋深变化过程，以及控制排水时段t末的地下水埋深h
w,i
，i为正交表中试验方案编号(i＝1,2,
…
,j2)。
27.进一步的，步骤(1)所述的n取值为2。
28.进一步的，步骤(3.1)所述的当地土壤水分特征参数值包括：土壤吸力f、土壤体积含水率γ；
29.进一步的，步骤(3.1)所述的土壤排水特性参数包括：土壤可排空体积v、潜水上升通量m。
30.进一步的，步骤(3.1)所述的排水系统特性参数包括：不透水层深度c、排水模数σ、最大地表蓄水深度h1、kirkham积水深度h2、控制排水堰深h3、排水控制时段t。
31.一种计算机可读存储介质，所述计算机可读存储介质存储有计算机程序，所述计算机程序被处理器执行时实现上述的方法步骤。
32.一种计算机设备，包括存储器、处理器及在所述存储器上存储并可运行的程序，所述程序被处理器执行时实现上述方法的步骤。
33.有益效果：本发明相对于现有技术，其显著优点：1、将大系统试验选优理论应用于解决暗管排水系统布置最优参数确定问题，针对多因素多水平非线性数学模型求解，可高效获取暗管布置优化参数，节省模型求解工作量，提高模型求解精度；2、结合平原河网地区实际，设置了满足平原河网地区要求的地下水埋深约束、排水控制时段约束以及满足旱作物正常生育需要和自流排水要求的暗管埋深约束，提高了模型求解结果的实用性。
附图说明
34.图1所示为本发明的流程图。
具体实施方式
35.下面结合附图对本发明的技术方案作进一步说明。
36.如图1所述，本发明所述的一种基于drainmod数值模拟暗管排水系统布置最优参数确定方法，包括以下步骤：
37.第一步：进行模型的构件；
38.以20年一遇设计暴雨，暴雨过后2-3天地下水降至田面以下最小埋深为目标，建立如下目标函数：
39.h＝min hw(d,h,p,s,t,f,γ,h1,h2,v,m,σ,r)
ꢀꢀꢀꢀ
(1)
40.式中，h为地下水最小埋深；hw为地下水埋深；d为暗管管径；h为暗管埋深；p为暗管间距；s为地表平整度；t为排水控制时段；f为土壤吸力；γ为土壤体积含水率；h1为最大地表蓄水深度；h2为kirkham积水深度；v为土壤可排空体积；m为潜水上升通量；σ为排水模数；r为暗管有效半径；
41.第二步，结合平原河网地区的实际情况，设计目标函数的约束条件；
42.1、地下水埋深h约束：
43.h≥0.8
ꢀꢀꢀ
(2)
44.2、排水控制时段t约束：
45.48≤t≤72
ꢀꢀꢀ
(3)
46.3、满足旱作物正常生育需要的暗管埋深h约束：
47.h≥hwꢀꢀꢀ
(4)
48.4、满足自流排水要求的暗管埋深h约束：
49.h≤h
农-0.1
ꢀꢀꢀ
(5)
50.式中，hw为旱作物生长的盛期要求的地下水埋深(m)；h
农
为农沟水面离地面距离(m)。
51.第三步，进行模型求解；
52.1、进行模型数据的准备，具体包括：
53.(1)测定田块面积a、田面高程z1、农沟水位z2；
54.(2)收集典型降雨资料，设计20年一遇设计暴雨过程；所述的设计暴雨为专有名词，按照过去若干年的长系列暴雨统计资料，按照水文学统计方法，排频计算出来；
55.(3)测定当地土壤水分特征参数值，包括：土壤吸力f、土壤体积含水率γ；测定土壤排水特性参数，包含：土壤可排空体积v、潜水上升通量m；
56.(4)设定排水系统特性参数，包含：不透水层深度d、排水模数σ、最大地表蓄水深度h1、kirkham积水深度h2、控制排水堰深h3、排水控制时段t；
57.(5)测定当地旱作物生长的盛期要求的地下水埋深hw值，以及测定当地农沟水面离地面距离h
农
值。
58.2、采用大系统试验选优方法，根据约束条件对目标函数求解，具体包括以下步骤：
59.(1)考虑设置暗管管径d、暗管埋设深度h、暗管间距p、地表平整度s四种因素；并将
暗管管径d设置为d1、d2、d3三种不同水平，将暗管埋设深度h设置为h1、h2、h3三种不同水平，将暗管间距p设置为p1、p2、p3三种不同水平，将地表平整度s设置为s1、s2、s3三种不同水平；
60.(2)参照大系统试验选优理论，构建4因素3水平正交表，如下表1：
61.表1 4因素3水平l9(34)正交表
[0062][0063]
(3)对正交表中的每一行试验方案，分别开展drainmod数值模拟，结合约束条件，获取控制排水时段t末的地下水埋深，具体包括以下步骤：
[0064]
(a)设置drainmod模型输入参数，如表2所示：
[0065]
表2 drainmod模型主要输入参数
[0066][0067]
(b)对正交表中每一个方案，在20年一遇设计暴雨下，考虑地下水埋深约束、排水控制时段约束、暗管埋深约束条件，采用drainmod模型，分别模拟20年一遇设计暴雨雨后2-3天土壤地下水埋深变化过程，获得该控制排水时段t内地下水埋深变化过程，以及控制排水时段t末的地下水埋深h
w,i
，其中i为正交表中试验方案编号，i＝1,2,
…
,9。
[0068]
(4)根据每一行试验方案数值模拟结果，构造极差表。通过正交分析，获得全部组合(34)的理论最优暗管管径d
*
、暗管埋设深度h
*
、暗管间距p
*
、地表平整度s
*
。
[0069]
(5)将此理论最优参数组合，代入目标函数(1)，结合约束条件(2)-(5)，通过drainmod模型重新进行数值模拟，获得地下水最小埋深h。
[0070]
此外，本发明还提供一种计算机可读存储介质，所述计算机可读存储介质存储有计算机程序，所述计算机程序被处理器执行时实现上述的方法步骤。一种计算机设备，包括存储器、处理器及在所述存储器上存储并可运行的程序，所述程序被处理器执行时实现上述方法的步骤。