[发明专利]地面灌溉中地表水流运动的模拟方法

摘要

本发明公开了一种地面灌溉中地表水流运动的模拟方法，属于地面灌溉领域。该模拟方法包括：首先，选择向量形式的全水动力学方程作为模拟目标畦田地面灌溉中地表水流运动的地面灌溉模型。然后，对向量形式的全水动力学方程进行时空离散化处理，得到全水动力学方程的时空离散表达式；并将时空离散表达式中的和分别用q和h来表达。最后，将全水动力学方程的时空离散表达式整理成以q和h为未知数的地表水流运动模拟方程，通过求解该方程，获取目标畦田中每一测点处的q和h，从而对目标畦田地面灌溉中地表水流运动进行模拟。该模拟方法计算简便，有效缩短了模拟计算的时间，且计算精度高，适用于大规模地面灌溉中地表水流运动的模拟。

主权项

1.一种地面灌溉中地表水流运动的模拟方法，其特征在于，所述模拟方法包括以下步骤：步骤a、选择向量形式的全水动力学方程作为模拟目标畦田地面灌溉中地表水流运动的地面灌溉模型；所述向量形式的全水动力学方程为：其中，F为物理通量，U为因变量向量，S为源项向量，t为所述目标畦田以灌溉起始时刻作为原点的时间坐标，单位为s，x为所述目标畦田中任意测点与所述目标畦田的畦首之间的距离，单位为m；u为水流沿所述目标畦田畦长方向的垂向均布流速，单位为m/s；h为地表水深，单位为m；q为水流沿所述目标畦田畦长方向的单宽流量，单位为L/(s·m)，且q＝hu；g为重力加速度，单位为m/s²；i_c为实测得到的所述目标畦田地表水入渗率；b为实测得到的所述目标畦田地表相对高程，单位为m；n为已知的所述目标畦田的畦面糙率系数，单位为s/m^1/3；步骤b、对所述向量形式的全水动力学方程进行时空离散化处理，得到所述全水动力学方程的时空离散表达式；所述全水动力学方程的时空离散表达式为：其中，Δt为针对所述t得到的时间步长；Δx为针对所述x得到的空间步长；n_t为针对所述t得到的时间离散节点；i为针对所述x得到的空间离散节点；i+1/2为空间离散节点i与i+1之间的中点；i‑1/2为空间离散节点i与i‑1之间的中点；步骤c、提供利用所述q和所述h来表达所述的第一表达式；所述第一表达式的为：步骤d、分别提供利用所述q和所述h来表达所述和的第二表达式和第三表达式；所述第二表达式为：其中，所述第三表达式为：步骤e、将所述第一表达式、所述第二表达式和所述第三表达式代入所述全水动力学方程的时空离散表达式，得到以所述q和所述h作为未知数的地表水流运动的模拟方程，通过对所述地表水流运动的模拟方程进行求解，获取所述目标畦田中每一测点处的所述q和所述h，从而对地面灌溉中地表水流运动进行模拟。