地面灌溉中地表水流运动的模拟方法

地面灌溉中地表水流运动的模拟方法
【专利摘要】本发明公开了一种地面灌溉中地表水流运动的模拟方法，属于地面灌溉领域。该模拟方法包括：首先，选择向量形式的全水动力学方程作为模拟目标畦田地面灌溉中地表水流运动的地面灌溉模型。然后，对向量形式的全水动力学方程进行时空离散化处理，得到全水动力学方程的时空离散表达式；并将时空离散表达式中的和分别用q和h来表达。最后，将全水动力学方程的时空离散表达式整理成以q和h为未知数的地表水流运动模拟方程，通过求解该方程，获取目标畦田中每一测点处的q和h，从而对目标畦田地面灌溉中地表水流运动进行模拟。该模拟方法计算简便，有效缩短了模拟计算的时间，且计算精度高，适用于大规模地面灌溉中地表水流运动的模拟。
【专利说明】
地面灌溉中地表水流运动的模拟方法
技术领域
[0001] 本发明涉及地面灌溉领域，特别涉及一种地面灌溉中地表水流运动的模拟方法。
【背景技术】
[0002] 地面灌溉是一种常见的农田灌溉方法，其利用灌水沟、畦或格田等，将水引入农 田，并使水在重力和毛细管作用下渗入土壤，为田间作物补充水分。进行地面灌溉时，为保 证田间作物得到充分灌溉，并尽可能节约水资源，需要对水流在目标畦田中的灌溉性能进 行预先评价，以对灌入水量及流速等进行调整，避免水资源的浪费。因此，提供一种地面灌 溉的评价方法是十分必要的。
[0003] 现有技术提供了一种地面灌溉中地表水流运动的模拟方法，该方法采用具有向量 形式的全水动力学方程作为地面灌溉模型，根据此来描述地面灌溉中地表水流运动过程。 在求解具有向量形式的全水动力学方程时，首先对其进行时空离散化处理，得到时空离散 表达式：
[0004]
[0005] 式中，F为数值通量，U为因变量向量，S为源项向量；At和Ax分别为数值计算的时 间步长和空间步长;nt为时间离散迭代步；i为空间离散节点；i+1/2为空间离散节点i与i+1 之间的中点；i-1/2为空间离散节点i与i_l之间的中点。
[0006] 其中，在求解全水动力学方程的上述时空离散表达式时，关键是计算
印 ，现有技术通常采用经典Roe模拟方法进行计算，该经典Roe模拟方法的表达式如下所
示：
[0007]
[0008]
[0009]
[0010]
[0011]
[0021 ]其中，u为水流沿所述目标畦田畦长方向的垂向均布流速，单位为m/s; h为地表水 深，单位为m; q为水流沿所述目标畦田畦长方向的单宽流量，单位为m3/ (s · m)。
[0022]根据上述各计算公式能够求解得到关于q和h的方程，通过该方程绘制h相对于q的 变化曲线，即获得地表水深相对于水流沿目标畦田畦长方向的单宽流量变化情况，以完成 地面灌溉中地表水流运动的模拟，从而实现对地面灌溉的评价。
[0023]发明人发现现有技术至少存在以下技术问题：
[0024]采用现有技术对地面灌溉中地表水流运动进行模拟计算时，虽然经典Roe模拟方 法的模拟精度较高，但其表达式极为复杂，导致求解效率较低，无法有效完成对大规模地面 灌溉水流运动的模拟。

【发明内容】

[0025] 本发明实施例所要解决的技术问题在于，提供了一种计算简单快捷、且精度高的 地面灌溉中地表水流运动的模拟方法，具体技术方案如下：
[0026] 本发明实施例提供了一种地面灌溉中地表水流运动的模拟方法，该模拟方法包括 以下步骤：
[0027]步骤a、选择向量形式的全水动力学方程作为模拟目标畦田地面灌溉中地表水流 运动的地面灌溉模型。
[0028]所述向量形式的全水动力学方程为：
[0029]
[0030]
[0031]
[0032]
[0033] 其中，F为物理通量，U为因变量向量，S为源项向量，t为所述目标畦田以灌溉起始 时刻作为原点的时间坐标，单位为s，x为所述目标畦田中任意测点与所述目标畦田的畦首 之间的距离，单位为m; u为水流沿所述目标畦田畦长方向的垂向均布流速，单位为m/s; h为 地表水深，单位为m;q为水流沿所述目标畦田畦长方向的单宽流量，单位为m3/(s · m)，且q =hu; g为重力加速度，单位为m/s2; i。为实测得到的所述目标畦田地表水入渗率;b为实测得 到的所述目标畦田地表相对高程，单位为m;n为已知的所述目标畦田的畦面糙率系数，单位 为 s/m1/3。
[0034]步骤b、对所述向量形式的全水动力学方程进行时空离散化处理，得到所述全水动 力学方程的时空离散表达式。
[0035]所述全水动力学方程的时空离散表达式为：
[0036]
[0037] 其中，A t为针对所述t得到的时间步长；△ X为针对所述X得到的空间步长;nt为针 对所述t得到的时间离散节点；i为针对所述X得到的空间离散节点；i+1/2为空间离散节点i 与i+Ι之间的中点；i-1/2为空间离散节点i与i-Ι之间的中点。
[0038]步骤c、提供利用所述q和所述h来表达所述\~的第一表达式。
[0039]所述第一表达式的为：
[0040]

[0041] 步骤d、分别提供利用所述q和所述h来表达所过 的第二表达式和 第三表达式。
[0042] 所述第二表达式为：
[0048] 步骤e、将所述第一表达式、所述第二表达式和所述第三表达式代入所述全水动力 学方程的时空离散表达式，得到以所述q和所述h作为未知数的地表水流运动的模拟方程， 通过对所述地表水流运动的模拟方程进行求解，获取所述目标畦田中每一测点处的所述q 和所述h，从而对地面灌溉中地表水流运动进行模拟。
[0049] 具体地，作为优选，所述目标畦田地表水入渗率i。通过Kostiako V入渗公式计算得 至1J，所述Kostiakov入渗公式为：
[0050]
[0051 ] 其中，kin为土壤入渗参数，单位为m/sa;a为无量纲入渗指数参数;Tin为入渗受水时 间，单位为s。
[0052] 具体地，作为优选，所述畦面糙率系数η通过在所述目标畦田的畦首实测畦首地表 水深值，并根据所述畦首地表水深值估算得到。
[0053] 具体地，作为优选，所述时间步长为I-IOs。
[0054] 具体地，作为优选，所述空间步长为l-10m。
[0055] 具体地，作为优选，在所述步骤e中，在求解所述地表水流运动的模拟方程时，将水 流首次进入畦田时的单宽流量作为求解所述q的初始条件。
[0056] 本发明实施例提供的技术方案带来的有益效果是：
[0057] 本发明实施例提供的地面灌溉中地表水流运动的模拟方法，通过对向量形式的全 水动力学方程进行时空离散，获取该全水动力学方程的时空离散表达式，并将该表达式转 换成以单宽流量q和地表水深h作为未知数的地表水流运动的模拟方程。通过求解该模拟方 程，即可获取目标畦田中每一测点处的q和h，进而可获得目标畦田中所有测点处的q和h，进 而通过所求解得到的数据来对地面灌溉中地表水流运动进行模拟。本发明实施例所使用的 上述时空离散表达式简便可靠，能在保证计算精确度的前提下，有效缩短q和h的计算时间， 提高了地面灌溉中地表水流运动模拟数据的计算效率，而且，该方法能够适用于农田等流 速较缓的地面灌溉中地表水流运动的模拟，并可用于大规模地面灌溉中地表水流运动的模 拟。
【附图说明】
[0058] 为了更清楚地说明本发明实施例中的技术方案，下面将对实施例描述中所需要使 用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于 本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他 的附图。
[0059] 图1是本发明实施例提供的地面灌溉中某一时刻下地表水流运动的状态示意图；
[0060] 图2是本发明又一实施例提供的河北冶河灌区、灌溉时间为40min时的地面灌溉中 地表水流运动的模拟结果示意图；
[0061] 图3是本发明又一实施例提供的河北冶河灌区、灌溉时间为SOmin时的地面灌溉中 地表水流运动的模拟结果示意图；
[0062] 图4是本发明又一实施例提供的新疆建设兵团222团、灌溉时间为40min时的地面 灌溉中地表水流运动的模拟结果示意图；
[0063] 图5是本发明又一实施例提供的新疆建设兵团222团、灌溉时间为SOmin时的地面 灌溉中地表水流运动的模拟结果示意图；
[0064]图6是本发明又一实施例提供的山东簸箕李灌区、灌溉时间为40min时的地面灌溉 中地表水流运动的模拟结果示意图；
[0065]图7是本发明又一实施例提供的山东簸箕李灌区、灌溉时间为SOmin时的地面灌溉 中地表水流运动的模拟结果示意图。
[0066]附图标记分别表示：
[0067] h地表水深，
[0068] d入渗水深，
[0069] b地表高程。
【具体实施方式】
[0070] 为使本发明的目的、技术方案和优点更加清楚，下面将对本发明实施方式作进一 步地详细描述。
[0071] 本发明实施例提供了一种地面灌溉中地表水流运动的模拟方法，该模拟方法包括 以下步骤：
[0072]步骤101、选择向量形式的全水动力学方程作为模拟目标畦田地面灌溉中地表水 流运动的地面灌溉模型。
[0073]向量形式的全水动力学方程为：
[0074]
[0075]
[0076]
[0077]
[0078] 其中，F为物理通量，U为因变量向量，S为源项向量，t为目标畦田以灌溉起始时刻 作为原点的时间坐标，单位为s，x为目标畦田中任意测点与目标畦田的畦首之间的距离，单 位为m; u为水流沿目标畦田畦长方向的垂向均布流速，单位为m/s; h为地表水深，单位为m; q 为水流沿目标畦田畦长方向的单宽流量，单位为m3/(s · m)，且q = hu;g为重力加速度，单位 为m/s2; ic为实测得到的目标畦田地表水入渗率;b为实测得到的目标畦田地表相对高程，单 位为m; η为已知的目标畦田的畦面糙率系数，单位为s/m1/3。
[0079]步骤102、对向量形式的全水动力学方程进行时空离散化处理，得到全水动力学方 程的时空离散表达式。
[0080] 全水动力学方程的时空离散表达式为：
[0081]
[0082] 其中，Δ t为针对t得到的时间步长；Δ X为针对X得到的空间步长;nt为针对t得到 的时间离散节点；i为针对X得到的空间离散节点；i +1 / 2为空间离散节点i与i +1之间的中 点；i-1/2为空间离散节点i与i-Ι之间的中点；
[0083] 步骤103、提供利用q和h来表?
的第一表达式。
[0084] 第一表达式的为：
[0085]
〇
[0086] 步骤104、分别提供利用q和h来表i
^的第二表达式和第三表达式。
[0087] 第二表达式为：
[0088]
[0089]
[0090]
[0091]
[0092]
[0093] 步骤105、将苐一表达式、苐二表达式和苐三表达式代入全水动力学万程的时空离 散表达式，得到以q和h作为未知数的地表水流运动的模拟方程，通过对地表水流运动的模 拟方程进行求解，获取目标畦田中每一测点处的q和h，从而对地面灌溉中地表水流运动进 行模拟。
[0094] 本发明实施例提供的地面灌溉中地表水流运动的模拟方法，通过对向量形式的全 水动力学方程进行时空离散，获取该全水动力学方程的时空离散表达式，并将该表达式转 换成以单宽流量q和地表水深h作为未知数的地表水流运动的模拟方程。通过求解该模拟方 程，即可获取目标畦田中每一测点处的q和h，进而可获得目标畦田中所有测点处的q和h，进 而通过所求解得到的数据来对地面灌溉中地表水流运动进行模拟。本发明实施例所使用的 上述时空离散表达式简便可靠，能在保证计算精确度的前提下，有效缩短q和h的计算时间， 提高了地面灌溉中地表水流运动模拟数据的计算效率，而且，该方法能够适用于农田等流 速较缓的地面灌溉中地表水流运动的模拟，并可用于大规模地面灌溉中地表水流运动的模 拟。
[0095] 步骤101中，向量形式的全水动力学方程中的t为目标畦田以灌溉起始时刻作为原 点的时间坐标，目标畦田空间内所有测点采用统一的时间坐标，计算时需要对t进行时间离 散;X为目标畦田任意测点与目标畦田的畦首之间的距离，即将目标畦田的畦首位置定义为 〇m，则目标畦田中任意一个测点所在位置与畦首之间的距离则定义为xm，该距离可以通过 实测得到。根据向量形式的全水动力学方程，对于因变量向量U求其在t轴方向上的偏导数， 对于物理通量F求其在X轴方向上的偏导数，得到源项向量S，以对目标畦田中每个测点处的 q和h相对于t和X的变化情况进行模拟计算。
[0096] 可以理解的是，u为水流沿目标畦田畦长方向的垂向均布流速，其指的是与目标畦 田畦长方向垂直的水流速度;如附图1所示，h为地表水深，即目标畦田中某个测点处土壤表 面以上的灌溉水深度;q为水流沿目标畦田畦长方向的单宽流量，即目标畦田中单位面积和 单位时间内流过的水的体积，因此，可以得出q = hu。在向量形式的全水动力学方程中，上述 q、h、u三个量均为未知量，其中U可用q和h来表达，进而可得出仅仅关于q和h的方程，以对地 面灌溉中地表水流的运动进行模拟。
[0097]具体地，g为重力加速度，可以根据目标畦田所在的地理位置，确定重力加速度的 值，一般选取9.8m/s2。
[0098] ic为实测得到的目标畦田地表水入渗率，即目标畦田中地表水渗入土壤的速度， 通过地表水入渗率，即可求出某一水流推进时间内，目标畦田内某测点处水渗入土壤的深 度，即地表水的入渗水深d，其等于i。乘以t，从而对地面灌溉的情况进行模拟和预测。如附 图1所示，b为实测得到的目标畦田在某一测点处的地表高程b，即目标畦田相对于地表的高 度。地表高程b可以采用水准仪在间隔预设距离上进行测量，其中预设距离可以为5m或者其 他数值，例如4!11、6 111、8111、1〇111等，本领域技术人员可以根据实际情况进行设定。地表高程13可 以作为地表水深h以及入渗水深d的基准线，便于更直观地得到以地表高程b、地表水深h、入 渗水深d随畦田长度的变化关系来表示地面灌溉中地表水流运动的模拟曲线。
[0099] η为目标畦田的畦面糙率系数，在本发明实施例中，其作为已知值，通常情况下，畦 面糙率系数多在〇. 04-0.3之间选择，其具体值可通过如下方法来确定:对目标畦田畦首处 的地表水深h进行观测，获取目标畦田畦首处的地表水深h的观测数据，根据该观测数据和 Chgzy公式估算畦面糙率系数。本领域技术人员可以理解的是，目标畦田不同位置处的糙率 系数可能不相同，可以以不同位置处的糙率系数的平均值作为该目标畦田的糙率系数。通 过上述实测或估算得到的参数来辅助向量形式的全水动力学方程的求解，以实现地面灌溉 中对地表水流的运动的模拟。
[0100] 虽然在本发明实施例中，向量形式的全水动力学方程也是关于q和h的表达式，但 是该表达式为一个连续方程，无法得到目标畦田中每一测点处的q和h。所以在求解每一测 点处的q和h时，需要对向量形式的全水动力学方程进行离散，然后通过它的离散表达式进 行求解。
[0101] 向量形式的全水动力学方程的离散表达式为：
[0102]
[0103] 具体地，本发明实施例提供的地面灌溉中地表水流运动的模拟方法，由于其向量 形式的全水动力学方程是对时间和空间（即长度)进行的偏导，所以需要对向量形式的全水 动力学方程进行时空离散化处理。进行时空离散化处理时，t为目标畦田沿水流方向上的水 流推进时间，t被离散成0、1、2、3.......m、nt+1等有限个时间节点，时间节点之间的时间差 就是针对t得到的时间步长Δ t，Δ t可以选择I-IOs中的任意数值，例如选择Is、2s、5s、IOs 等;x为目标畦田水流方向上的长度，x也被离散成0、I、2、3.......i、i + Ι等有限个空间节 点，i + 1/2即表示空间离散节点i与i+Ι之间的中点；i-1/2为空间离散节点i与i-Ι之间的中 点，空间节点之间的距离就是针对X得到的空间步长A Χ，ΔΧ可以选择I-IOm中的任意数值， 例如选择1111、21]1、51]1、1〇1]1等。&1：和&1值是对时间和空间进行不同程度的离散，&1：和&1的 取值越小，说明离散程度越高，计算结果的精确度也越高。
[0104] 具体地，在求解上述全水动力学方程的时空离散表达式时，首先求解其中的
的计算公式中涉及到目标畦田地表水入渗率ic，ic是通过Kostiakov入渗公式计 算得到，Kostiakov入渗公式为：
[0105]
[0106] 其中，kin为土壤入渗参数，单位为m/sa;a为无量纲入渗指数参数;τ?η为入渗受水时 间，单位为s。具体地，通过双环渗水试验，可以测得目标畦田的土壤入渗参数以及无量纲入 渗指数参数。双环入渗水试验是一种常用的土壤入渗参数测量方法，其所用设备主要包括 内外两个同心设置的圆环以及马利奥特瓶。将内环和外环插入土层一定深度，利用马利奥 特瓶向内环注水，连续观测并记录马利奥特瓶内水位随时间的变化，即可获得目标畦田的 土壤入渗参数以及无量纲入渗指数参数。入渗受水时间即试验土壤从有水渗入直至土壤表 面没有水的时间，且目标畦田中每一个测点处的入渗受水时间是不同的。通过上述三个参 数，可以计算得到目标畦田的地表水入渗率。
[0107] 本领域技术人员可以理解的是，在一块目标畦田内，通常会设置多个土壤入渗参 数测点，每一个测点最终得到的土壤入渗参数kin以及无量纲入渗指数参数a可能会有所不 同。当多个测点得到的土壤入渗参数kin以及无量纲入渗指数参数α不相同时，在计算该目标 畦田的地表水入渗率时所用的土壤入渗参数kin以及无量纲入渗指数参数α的取值应当为各 测点得到的土壤入渗参数kin以及无量纲入渗指数参数α的平均值。
[0108] 具体地，求解全水动力学方程的时空离散表达式的关键步骤是求解其中的迎风格
式耗散项，即 和h表示出来，即用第二表达式和第三 表达式来表彡 通过经典Roe模拟方法)（
进行求解，经典Roe模拟方法一般适用于洪水等湍急的水流模式，求解过程十分复杂，计算 量庞大。而本发明实施例提供的地面灌溉中地表水流运动的模拟方法，对经典Roe模拟方法 中求得的一系列F值进行了算数平均，僅
的计算过程得到了简化。这种简化 使是建立在地面灌溉模式下的，且模拟出的数值曲线与经典Roe模拟方法十分接近，说明该 模拟方法准确度较高，其完全适用于农田地面灌溉中地表水流运动的模拟。采用本发明实 施例提供的地面灌溉中地表水流运动的模拟方法，缩短了模拟计算的时间，提高了模拟计 算的效率，使该模拟方法在进行大规模地面灌溉水流运动的模拟时，具有更突出的优势。
[0109] 具体地，通过将第一表达式、第二表达式和第三表达式代入全水动力学方程的时 空离散表达式，得到以q和h为未知数的地表水流运动的模拟方程，该模拟方程为：
[0112] 该模拟方程中b、n、kin、a、Tin、Δ X、Δ t为已知量，nt和i根据迭代的步数具体进行取 值，将q的初始条件设置为畦田内所有空间点均为无单宽流量值(即零值），边界条件为畦田 首部为实测获得的入畦单宽流量，该值采用超声波流量计进行测量。
[0113] 由于u = q/h，故通过联立求解上述两式，即可计算出每一测点处的q和h，并得到h 关于q变化的曲线，以此判断地面灌溉中地表水流沿畦长方向各点的入渗情况，从而选择合 适的灌溉水量，节约水资源。
[0114] 具体地，由于本发明实施例提供的模拟方法适用于农田的地面灌溉，所以沿目标 畦田畦长方向的单宽流量q的适用范围为0.002-0.01m 3/(s ·πι)，对具有适宜q值的地面灌 溉过程进行模拟，可以获得更加精准的模拟结果。
[0115] 以下将通过具体实施例进行详细阐述：
[0116] 实施例1
[0117] 本实施例提供了一种地面灌溉中地表水流运动的模拟方法，通过该模拟方法对河 北冶河灌区畦田中地面灌溉时地表水流的运动情况进行模拟，选取的河北冶河灌区畦田的 基本情况为：畦长l〇〇m，通过超声波流量计测得水流首次进入畦田时的单宽流量为3L/s · m，以作为求解q的初始条件，通过秒表测得畦田的受水时间为lOOmin，通过双环渗水试验测 得畦田的土壤入渗参数kin为0.00063m/s a，无量纲入渗指数参数α为〇. 23，通过畦首实测的 地表水深h估算得到畦面糙率系数η为0.08s/m1/3。
[0118] 具体模拟步骤如下：
[0119] 步骤1、选择向量形式的全水动力学方程作为模拟目标畦田地面灌溉中地表水流 运动的地面灌溉模型，其中，向量形式的全水动力学方程为：
[0120]
[0121]
[0122]
[0123]
[0124] 其中，F为物理通量，U为因变量向量，S为源项向量，t为目标畦田沿水流方向上的 水流推进时间，X为目标畦田水流方向上的长度;U为水流沿目标畦田畦长方向的垂向均布 流速，单位为m/s ;h为地表水深h，单位为m; q为水流沿目标畦田畦长方向的单宽流量，单位 为m3/(s ·πι)，且q = hu;g为重力加速度，选择9.8m/S2;ic为实测得到的目标畦田地表水入渗 率;b为实测得到的目标畦田地表高程b，单位为m;n为已知的目标畦田的畦面糙率系数，即 0.08s/m 1/3。
[0125] 其中，i。通过Kostiakov入渗公式计算得到，在本实施例中i。为：
[0126]
[0127] 步骤2、对向量形式的全水动力学方程进行时空离散化处理，得到全水动力学方程 的时空离散表达式：
[0128]
[0129] At为目标畦田受水时间变量所对应的时间步长，本实施例中At为Is; ΔΧ为目标 畦田的畦长变量所对应的空间步长，本实施例中A X为Im; nt为时间离散迭代步；i为空间离 散节点；i+1/2为空间离散节点i与i+Ι之间的中点；i-1/2为空间离散节点i与i-Ι之间的中 点。
[0130] 步骤3、提供利用q和h来表这
的第一表达式，具体地，第一表达式的为：
[0138]步骤5、将第一表达式、第二表达式和第三表达式代入全水动力学方程的时空离散 表达式，得到以q和h为未知数的地表水流运动的模拟方程，
[0141] 其中，i = 0，l，2,3,4,.....，1(1为空间节点总数），11* = 0，1，2,3,4,.....，N(N为 时间节点总数）。
[0142] 求解时，初始条件q = 0,边界条件为q〇 = 0.003。
[0143] 通过求解目标畦田中每一测点处的q和h，从而获取目标畦田中所有测点处的q和 h，以对地面灌溉中地表水流运动进行模拟。
[0144] 同时，本实施例还采用经典Roe模拟方法对该同一个目标畦田进行模拟，并将本实 施例提供的模拟方法与经典Roe模拟方法的模拟数值、模拟精度以及模拟计算所需时间进 行对比，以对本实施例提供的模拟方法的模拟效果进行检验。
[0145] 附图2和附图3分别为河北冶河灌区的畦田灌溉40min和80min时，采用本实施例提 供的地面灌溉中地表水流运动的模拟方法与采用经典Roe模拟方法得到的地表水深h的数 值模拟对比图，并通过畦田地表水入渗率与受水时间计算入渗水深d，绘制入渗水深d曲线。 而单宽流量q与地表水深h相对应，当地表水深h的模拟值较为准确时，则单宽流量q的值也 是较为准确的。因此，本领域技术人员在评价地面灌溉中地表水流运动的模拟方法的可行 性时，只需参考地表水深h与入渗水深d之间的管线图即可，即如附图2、3中所示。从附图2和 附图3中可以看出，通过上述两种方法模拟出的数值，无论是地表水深h还是入渗水深d都十 分接近，表示这两者的曲线基本重合在一起。由于经典Roe模拟方法的已经得到广泛认可和 应用，说明其模拟真实度较高，模拟效果较好，而通过本实施例提供的模拟方法得到的模拟 数值与经典Roe模拟方法得到的数值十分相近，说明本实施例提供的模拟方法也具有较好 的模拟效果，可以用来模拟地面灌溉中地表水流的运动情况。
[0146] 具体地，采用本实施例提供的地面灌溉中地表水流运动的模拟方法测得的地表水 深h和入渗水深d与采用经典Roe模拟方法得到的数值相比，在灌水时间为40min时，两种方 法所得到的地表水深h之间的平均相对误差为0.12%，入渗水深d之间的平均相对误差为 0.11 %。在灌水时间为SOmin时，两种方法所得到的地表水深h之间的平均相对误差为 0.09%，入渗水深d之间的平均相对误差为0.11 %。可见，利用本实施例提供的方法与经典 Roe模拟方法所模拟得到的数据之间的误差很小，模拟效果良好。
[0147] 然而，实验证明，利用本实施例提供的模拟方法进行上述各数值的模拟计算的时 间为2.92min，用经典Roe模拟方法进行模拟计算的时间则为9.42min，可见，采用本实施例 提供的模拟方法进行地面灌溉中地表水流运动的模拟计算，可使计算时间缩短69%，从而 提高了模拟计算的效率，更适用于大规模畦田的模拟计算。
[0148] 实施例2
[0149 ]本实施例提供了一种地面灌溉中地表水流运动的模拟方法，该方法的操作过程与 实施例1相同，所不同的是：
[0150] 进行试验的畦田为新疆建设兵团222团畦田，该畦田的基本情况为：畦长150m，通 过超声波流量计测得水流首次进入畦田时的单宽流量为6L/s · m( 即q的初始条件qo = 0.006 )，通过秒表测得畦田的受水时间为160min，通过双环渗水试验测得畦田的土壤入渗 参数kin为0.00016m/s a，无量纲入渗指数参数α为0.42，通过畦首实测的地表水深h估算得到 畦面糙率系数η为〇. 11 s/m1/3。
[0151]同时，采用经典Roe模拟方法对该同一畦田进行模拟，并将本实施例提供的模拟方 法与经典Roe模拟方法的模拟数值、模拟精度以及模拟计算所需时间进行对比，以对本实施 例提供的模拟方法的模拟效果进行检验。
[0152]附图4和附图5分别为新疆建设兵团222团的畦田灌溉40min和80min时，采用本实 施例提供的地面灌溉中地表水流运动的模拟方法与采用经典Roe模拟方法得到的地表水深 h的数值模拟对比图，并通过畦田地表水入渗率与受水时间计算入渗水深d，绘制入渗水深d 曲线。从附图4和附图5中可以看出，通过上述两种方法模拟出的数值，无论是地表水深h还 是入渗水深d都十分接近。由于经典Roe模拟方法的已经得到广泛认可和应用，说明其模拟 真实度较高，模拟效果较好，而通过本实施例提供的模拟方法得到的模拟数值与经典Roe模 拟方法得到的数值十分相近，说明本实施例提供的模拟方法也具有较好的模拟效果，可以 用来模拟地面灌溉中地表水流的运动情况。
[0153] 进一步地，采用本实施例提供的地面灌溉中地表水流运动的模拟方法测得的地表 水深h和入渗水深d与采用经典Roe模拟方法得到的数值相比，在灌水时间为40min时，地表 水深h之间的平均相对误差为0.13%，入渗水深d之间的平均相对误差为0.09%，在灌水时 间为SOmin时，地表水深h之间的平均相对误差为0.08%，入渗水深d之间的平均相对误差为 0.08%。可见，利用本实施例提供的方法与经典Roe模拟方法所模拟得到的数据之间的误差 很小，模拟效果良好。
[0154] 然而，试验证明，采用本实施例提供的模拟方法进行模拟计算的时间为4.25min， 用经典Roe模拟方法进行模拟计算的时间则为13.67min，可见，采用本实施例提供的模拟方 法进行地面灌溉中地表水流运动的模拟计算，可使计算时间缩短69%，从而提高了模拟计 算的效率，更适用于大规模畦田的模拟计算。
[0155] 实施例3
[0156] 本实施例提供了一种地面灌溉中地表水流运动的模拟方法，该方法的操作过程与 实施例1相同，所不同的是：
[0157] 进行试验的畦田为山东簸箕李灌区畦田，该畦田的基本情况为:畦长130m，通过超 声波流量计测得水流首次进入畦田时的单宽流量为4L/s · m(即q的初始条件qo = 0.004)， 通过秒表测得畦田的受水时间为120min，通过双环渗水试验测得畦田的土壤入渗参数kin为 0.00026m/s a，无量纲入渗指数参数α为0.38，通过畦首实测的地表水深h估算得到畦面糙率 系数 η 为 0.09s/m1/3。
[0158] 同时，采用经典Roe模拟方法对河北冶河灌区的同一畦田进行模拟，并将本实施例 提供的模拟方法与经典Roe模拟方法的模拟数值、模拟精度以及模拟计算所需时间进行对 比，以对本实施例提供的模拟方法的模拟效果进行检验。
[0159] 附图6和附图7分别为河北冶河灌区的畦田灌溉40min和80min时，采用本实施例提 供的地面灌溉中地表水流运动的模拟方法与采用经典Roe模拟方法得到的地表水深h的数 值模拟对比图，并通过畦田地表水入渗率与受水时间计算入渗水深d，绘制入渗水深d曲线。 从附图6和附图7中可以看出，通过上述两种方法模拟出的数值，无论是地表水深h还是入渗 水深d都十分接近。由于经典Roe模拟方法的已经得到广泛认可和应用，说明其模拟真实度 较高，模拟效果较好，而通过本实施例提供的模拟方法得到的模拟数值与经典Roe模拟方法 得到的数值十分相近，说明本实施例提供的模拟方法也具有较好的模拟效果，可以用来模 拟地面灌溉中地表水流的运动情况，且地面灌溉时间越长，模拟精度越高。
[0160] 进一步地，采用本实施例提供的地面灌溉中地表水流运动的模拟方法测得的地表 水深h和入渗水深d与采用经典Roe模拟方法得到的数值相比，在灌水时间为40min时，地表 水深h之间的平均相对误差为0.12%，入渗水深d之间的平均相对误差为0.11%，在灌水时 间为SOmin时，地表水深h之间的平均相对误差为0.10%，入渗水深d之间的平均相对误差为 0.11%。可见，利用本实施例提供的方法与经典Roe模拟方法所模拟得到的数据之间的误差 很小，模拟效果良好。
[0161] 然而，试验证明，采用本实施例提供的模拟方法进行模拟计算的时间为3.61min， 用经典Roe模拟方法进行模拟计算的时间则为11.58min，可见，采用本实施例提供的模拟方 法进行地面灌溉中地表水流运动的模拟计算，可使计算时间缩短69%，从而提高了模拟计 算的效率，更适用于大规模畦田的模拟计算。
[0162]以上所述仅为本发明的较佳实施例，并不用以限制本发明的保护范围，凡在本发 明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围 之内。
【主权项】
1. 一种地面灌溉中地表水流运动的模拟方法，其特征在于，所述模拟方法包括以下步 骤： 步骤a、选择向量形式的全水动力学方程作为模拟目标畦田地面灌溉中地表水流运动 的地面灌溉模型； 所述向量形式的全水动力学方程为：其中，F为物理通量，U为因变量向量，S为源项向量，t为所述目标畦田以灌溉起始时刻 作为原点的时间坐标，单位为s，x为所述目标畦田中任意测点与所述目标畦田的畦首之间 的距离，单位为m;u为水流沿所述目标畦田畦长方向的垂向均布流速，单位为m/s;h为地表 水深，单位为m;q为水流沿所述目标畦田畦长方向的单宽流量，单位为m 3/(s ? m)，且q = hu; g为重力加速度，单位为m/s2; i。为实测得到的所述目标畦田地表水入渗率;b为实测得到的 所述目标畦田地表相对高程，单位为m;n为已知的所述目标畦田的畦面糙率系数，单位为s/ m1/3； 步骤b、对所述向量形式的全水动力学方程进行时空离散化处理，得到所述全水动力学 方程的时空离散表达式； 所述全水动力学方程的时空离散表达式为：其中，△ t为针对所述t得到的时间步长；A x为针对所述x得到的空间步长;nt为针对所 述t得到的时间离散节点；i为针对所述x得到的空间离散节点；i+1/2为空间离散节点i与i+ 1之间的中点；i-1/2为空间离散节点i与i-1之间的中点； 步骤c、提供利用所述q和所述h来表达所述Sf4的第一表达式； 所述第一表达式的为：步骤d、分别提供利用所述q和所述h来表达所述的第二表达式和第三表 达式；步骤e、将所述第一表达式、所述第二表达式和所述第三表达式代入所述全水动力学方 程的时空离散表达式，得到以所述q和所述h作为未知数的地表水流运动的模拟方程，通过 对所述地表水流运动的模拟方程进行求解，获取所述目标畦田中每一测点处的所述q和所 述h，从而对地面灌溉中地表水流运动进行模拟。2. 根据权利要求1所述的模拟方法，其特征在于，所述目标畦田地表水入渗率i。通过 Kostiakov入渗公式计算得到，所述Kostiakov入渗公式为：其中，kin为土壤入渗参数，单位为m/sa; a为无量纲入渗指数参数;Tin为入渗受水时间， 单位为s。3. 根据权利要求2所述的模拟方法，其特征在于，所述畦面糙率系数n通过在所述目标 畦田的畦首实测畦首地表水深值，并根据所述畦首地表水深值估算得到。4. 根据权利要求1所述的模拟方法，其特征在于，所述时间步长为l-10s。5. 根据权利要求1所述的模拟方法，其特征在于，所述空间步长为l-10m。6. 根据权利要求1所述的模拟方法，其特征在于，在所述步骤e中，在求解所述地表水流 运动的模拟方程时，将水流首次进入畦田时的单宽流量作为求解所述q的初始条件。
【文档编号】G06F17/50GK106055774SQ201610363355
【公开日】2016年10月26日
【申请日】2016年5月26日
【发明人】章少辉, 许迪, 李益农, 白美健, 李福祥
【申请人】中国水利水电科学研究院