水盐模拟模型中暗管排水排盐边界的改进方法转让专利
申请号 : CN202010663696.1
文献号 : CN111896429B
文献日 : 2021-08-03
发明人 : 曾文治 , 刘义 , 敖畅 , 雷国庆 , 蒋永华 , 黄介生
申请人 : 武汉大学
摘要 :
权利要求 :
1.一种水盐模拟模型中暗管排水排盐边界的改进方法,其特征在于,包括如下步骤:S1:获取试验地的相关数据,根据试验地情况设计室内暗管排水试验以获得暗管排水排盐数据,其中,室内试验土壤情况与试验地一致;
S2:在水盐模拟模型的暗管原始边界周围建立虚拟土层边界,将试验地相关数据和获得的暗管排水排盐数据作为水盐模拟模型的输入参数,对虚拟土层的土壤水力参数进行率定,找出模拟的排水排盐量与实测的排水排盐量误差最小的一组土壤水力参数作为最优参数;
S3:以精确预测暗管排水排盐过程为目标,采用线性回归分析方法,建立不同模拟区域中虚拟土层的环宽RW与土壤水力参数中的饱和导水率Ks的函数关系,再将步骤S2中得到的最优参数中的饱和导水率Ks代入步骤S3得到的Ks函数中从而得到虚拟土层的环宽RW;
S4:将步骤S3得到的虚拟土层环宽RW作为改进后的水盐模拟模型的暗管边界参数,根据需要获得暗管排水过程的试验地的实际灌水情况对水盐模拟模型改进后的暗管边界进行验证和评价。
2.如权利要求1所述的水盐模拟模型中暗管排水排盐边界的改进方法,其特征在于,步骤S1中的室内暗管排水试验包括:
设计室内试验土槽,在土槽两个相邻的侧壁顶部设置有溢水孔,在土槽设置有溢水孔的一侧壁上还布置暗管埋设孔,在土槽底部设置排水管,排水管上设置有控制闸,此外土槽底部还设置有用于观测地下水的连通管,连通管上设置有刻度;
在试验地取土,保持土层的容重和层状分布与试验地原状土相同,在此过程中一并埋设暗管,暗管埋设在40‑60 cm的土层中,暗管外壁上包裹有一层可渗透的无纺布。
3.如权利要求1所述的水盐模拟模型中暗管排水排盐边界的改进方法,其特征在于,试验地相关数据包括原始土壤水力参数、各层土壤初始含水量与含盐量、气象数据以及地下水埋深数据,步骤S1中进行多次试验,每次施以不同的灌水量,暗管排水排盐数据包括暗管排水量和排盐量、水流速度、地下水位数据、土壤剖面各层最终的含水量与含盐量。
4.如权利要求1所述的水盐模拟模型中暗管排水排盐边界的改进方法,其特征在于,步骤S2具体还包括:
步骤1中的试验地相关数据包括原始土壤水力参数,先以原始土壤水力参数为基础对虚拟土层的土壤水力参数进行率定,然后将率定后的虚拟土层的土壤水力参数作为水盐模拟模型改进后的暗管边界的输入参数,最后分别运行原始边界和改进后的暗管边界模型,并评估改进后的暗管边界的模拟精度。
5.如权利要求4所述的水盐模拟模型中暗管排水排盐边界的改进方法,其特征在于,对虚拟土层的土壤水力参数进行率定的方法为:将率定时输入的土壤水力参数的变化范围设置为原始土壤水力参数的0.80‑1.20倍,其中,每层土壤的土壤水力参数变化相互独立,再通过水盐模拟模型计算得到多组模拟值,每组模拟值包括排水模拟值和排盐模拟值,然后通过数据处理软件将多组模拟值分别与对应的实测值进行比较,找出与实测值误差最小的一组模拟值对应的土壤水力参数作为最优参数。
6.如权利要求1所述的水盐模拟模型中暗管排水排盐边界的改进方法,其特征在于,步骤S3还包括:
设置不同的环宽RW,采用水盐模拟模型模拟计算不同环宽对应的各组排水量和排盐量,之后通过数据处理软件计算出与各组排水量和排盐量对应的饱和导水率Ks,得到多组不同的Ks值;
将不同的环宽RW与所对应的Ks值进行回归分析与线性拟合得到它们之间的函数关系。
说明书 :
水盐模拟模型中暗管排水排盐边界的改进方法
技术领域
背景技术
地下水位,以降低盐碱地的土壤盐分。但与沟渠排水相比,暗管排水提高了土地利用效率,
促进了农业机械化。此外,暗管排水增加了盐的横向排放,从而提高了盐的排放效率。然而,
地下管道的不合理布局往往需要更多的淋溶水,并导致更多的土壤养分流失和农业化学污
染。因此,有必要对土壤水盐运移过程进行精确量化,以优化暗管的布置和淋溶方案。
地的土壤质地和气候条件,缺乏普遍适用性的。此外,影响暗管布局的参数众多,无论是室
内试验还是野外试验都很难对这些影响因素进行系统的定量的研究。
参数容易获取,但对土壤中水和溶质运动的具体过程描述较少。虽然dual‑permeability模
型被广泛用于模拟田间暗管排水,但是模型参数的确定相对来说较为困难。相比之下,
HYDRUS‑2D模型具有灵活的边界条件,可以模拟二维饱和与非饱和土中的水盐运移规律。但
是现有的应用均是从土壤水、盐含量和其他化学物质迁移的角度来校准模型,忽略了排水
和排盐过程。其主要原因是HYDRUS‑2D中的渗透边界的排水能力远大于实际暗管,只运用渗
透边界很难模拟暗管的排水过程。
法的不足,有必要确定一种改进的暗管边界模型,进而能通过模型就能快速、系统地预测暗
管排水过程的方法。
发明内容
暗管边界,大大提高了模拟精度,解决现有技术中确定暗管排水过程费时费力和模型模拟
精度差的问题。
行率定,找出模拟的排水排盐量与实测的排水排盐量误差最小的一组土壤水力参数作为最
优参数;
数替换步骤S2中得到的最优参数中的饱和导水率Ks从而得到虚拟土层边界参数;
界进行验证和评价。
外土槽底部还设置有用于观测地下水的连通管,连通管上设置有刻度;
排水排盐数据包括暗管排水量和排盐量、水流速度、地下水位数据、土壤剖面各层最终的含
水量与含盐量。
盐模拟模型改进后的暗管边界的输入参数,最后分别运行原始边界和改进后的暗管边界模
型,并评估改进后的暗管边界的模拟精度。
参数变化相互独立,再通过水盐模拟模型计算得到多组模拟值,每组模拟值包括排水模拟
值和排盐模拟值,然后通过数据处理软件将多组模拟值分别与对应的实测值进行比较,找
出与实测值误差最小的一组模拟值对应的土壤水力参数作为最优参数。
组不同的Ks值;
从而刻画无纺布和小开孔率产生的阻力,从而使得重新建立的边界更接近实际管流情景,
即虚拟土层能够实现减缓暗管出流从而更能反映实际排水情况的目的,从而使得改进暗管
边界后的HYDRUS‑2D在不同模拟区域中实现了不仅对排水排盐量,而且对暗管排水排盐过
程都能准确预测;将改进暗管边界后的HYDRUS‑2D用于对复杂的田间试验进行系统的模拟,
能快速准确地模拟田间暗管排水过程,大大减少了时间成本和人工投入。
附图说明
具体实施方式
始含水量和含盐量可通过对土壤进行取样、测试及分析获得。其中,土壤水力参数包括残余
含水率(θr)、饱和含水率(θs)、影响水力函数形状的经验系数(α,n)、饱和导水率Ks,孔隙连
通性参数(l)。地下水埋深数据可通过布设水位观测井观测得到。
3.5cm处的正面侧壁上设置4个溢流孔11,在侧面侧壁上对应高度设置2个溢流孔11(每个溢
流孔2的直径 )。在土槽1正面侧壁上设置有一根暗管埋设孔12。在土槽1下方设
置地下水抽水孔。在土槽1的底部设置排水管,该排水管上设置有控制阀门。土槽1底部还设
置有用于观测地下水的连通管2,为了便于观察,该连通管2上设置有刻度。
质。为了模拟田间土壤分布,将土壤按照田间原状土形式分层填入土槽中。土壤分成4层,每
层厚度为20cm。每层的体积密度保持与田间相同。在土槽中填土时,一并将暗管埋设在土层
中,暗管为波纹管,其外径为54mm,内径为50mm,埋深50cm,其距土槽两侧均为0.75m。其中,
暗管开孔率为3%,暗管的坡度为2‰。在埋入之前,先在暗管外壁上包裹一层可渗透的无纺
布,以防止土壤颗粒堵塞暗管。
度、地下水位数据、土壤剖面各层最终的含水量与含盐量。具体的试验方法为:第一次试验
的总灌溉量为19.2cm,表层积水持续19.5h,入渗阶段地表积水均匀分布在表层土壤中。灌
溉4.2h后开始排水,排水历程持续了85.8h。第一次试验结束两周后进行第二次试验,其中,
第一次试验用于率定土壤水力参数,第二次试验用于验证所得的土壤水力参数。第二次试
验的总灌溉量为9.2cm,表面积水持续了9.6h。灌溉1.6h后开始排水,排水历程持续了
58.4h。用有机玻璃标尺测量表面积水的入渗流量,每30min记录一次;用翻斗式流量计测量
暗管排水流量,并且每30min收集一次排水样品,用电导率仪测量其电导率以分析暗管的排
盐量。
和75cm深度的土壤样品,每个土壤样品的体积为20cm 。用烘箱干燥法(105℃,8h)测定每个
样品的含水率从而获得土壤剖面的含水量。烘干后,每个土壤样品取15g,研磨,过筛,然后
放入含有75ml去离子水(水∶土壤质量比=1∶5)的150ml锥形瓶中。锥形烧瓶在振荡仪上充
分振荡。静置24小时后,用电导率仪测量上溶液过滤液的电导率以分析土壤剖面的含盐量。
总可溶性盐(TDS)含量由提取的EC1:5值计算,土壤盐含量和EC1:5之间的关系用干燥残渣法
校准如下:
行率定,找出模拟的排水排盐量与实测的排水排盐量误差最小的一组土壤水力参数作为最
优参数。
的暗管边界如图2所示。虚拟土层的主要作用是模拟无纺布和波纹管的结合界面从而刻画
无纺布层和小开孔率产生的阻力,图2中内环代表渗透边界,环内为虚拟土层。基于对实际
管流情景的模拟,虚拟土层能够实现减缓暗管出流从而更能反映实际排水情况的目的。
进后的暗管边界。
其中,虚拟土层的初始土壤水力参数设定为与暗管所在的40‑60cm土层对应的各项参数相
同。本发明先以原始土壤水力参数为基础对虚拟土层的环内土壤水力参数进行单独率定,
然后将率定后的环内土壤水力参数作为改进后的暗管边界的输入参数,最后分别运行原始
边界和改进后的暗管边界模型,并评估改进后的暗管边界的模拟精度。
的累积排水和累积排盐量。同时将模型中的观察点设在10、15、20、25、30、35、40、45、50、55、
60、65、70和75cm的深度来提取土壤剖面的盐分和水分含量。随后,通过PEST多目标校正方
法将测量值与模拟值进行比较并选出最优参数。以上率定过程的具体操作方法和过程如图
3所示,首先准备好PEST和HYDRUS‑2D所需的输入文件并进行初始化,即准备两者所需的输
入参数文件;然后读取并更改SELECTOR.IN文件中的土壤水力参数,更改方法为将该文件中
土壤水力参数的变化范围设置为原始土壤水力参数的0.8‑1.2倍,且每层土壤参数的变化
相互独立,即同层土壤水力参数中的各项参数用同一个比例,不同层土壤水力参数可以用
同一个比例也可用不同的比例。之后使用HYDRUS‑2D进行模拟计算并输出排水过程文件、排
盐过程文件、土壤含水率文件和土壤含盐率文件,从而获得各组模拟值,每组模拟值包括排
水模拟值和排盐模拟值。最后,使用PEST软件将以上输出的每组排水模拟值和排盐模拟值
与从试验中得到的对应的每组排水测量值和排盐测量值进行比较,选出与测量值误差最小
(RMSE和MAE最小)的模拟值对应的土壤水力参数,则认为此组土壤水力参数为最优参数,并
停止调参。否则,继续进行以上步骤,直到参数最优为止。
土壤饱和体积含水率(L·L );h为负压水头(L);K(h)为土壤非饱和导水率(L·T );Ks为
‑1
土壤饱和导水率(L·T );Se是有效饱和度;α是与土壤物理性质有关的经验参数;n是与多
孔介质有关的参数,其中m=1‑1/n;l为经验拟合参数,通常取值为0.5。
水动力弥散系数(L·T )。
言,很明显,在早期阶段,CWF被大大高估,曲线无法描述实际排水过程。对于改进后的暗管
边界来说,曲线形状与排水过程几乎一致。两种边界在率定和验证阶段的统计结果,见表1,
2
虚拟土层边界的RMSE和MAE分别为3.42L和2.91L,小于原始边界;且其决定系数(R)为
0.97,大于原始模型。结果表明,虚拟土层边界对于模拟排水过程的精度有了显著的提高。
与排水过程类似,对于排盐过程来说,虚拟土层边界在校准和验证阶段产生较小的RMSE和
2
MAE值以及较大的R值,这表明虚拟土层边界模型在CSF模拟中表现更好。
数替换步骤S2中得到的最优参数中的饱和导水率Ks参数从而得到虚拟土层边界参数。
水率对参数Ks和RW的敏感性,具体见表2。当参数Ks和RW变化50%分别导致暗管排水过程
(CWF)变化3.1%‑4.6%,排盐过程(CSF)变化3.4%‑4.5%;而土壤剖面盐分分布(SWC)和水
分分布(SSC)变化范围分别为0.6%‑7.14%和0.73%‑7.38%。这些结果表明暗管排水过程
和排盐过程对参数Ks有较高的敏感性,而土壤剖面盐分和水分对RW的敏感性较高。可见,Ks
和RW的变化对排水、排盐过程和土壤剖面水盐运移具有显著影响。为了定量表征它们之间
的关系,将虚拟土层环宽RW分别设置为0.5,1,1.5,2和2.5cm,采用HYDRUS‑2D模型计算不同
环宽对应的各组排水量和排盐量,然后通过数据处理软件PEST计算出与各组排水量和排盐
量对应的饱和导水率Ks,得到多组不同的Ks值。经过计算发现在不同的RW情景下,通过调整
Ks,该模型能够很好地模拟土壤剖面的盐分和水分以及暗管的排水和排盐过程。将所得的
RW与其对应的Ks进行线性拟合,拟合效果较好(图6),关系如下:
范围更广。
界进行验证和评价。
‑1
暗管半径为5cm,设置环宽RW为2cm,则对应的Ks取值为0.94cm·h 。HYDRUS‑2D模型输入数
据包括土壤数据、降雨灌溉数据、改进后的暗管边界条件数据。验证过程包括设置输入数
据、绘制单个暗管埋设下的土壤剖面图、确定初始条件、设置边界条件、模拟计算。经过模拟
2
和验证,模拟结果如图7所示,且模拟值与测量值之间的R 达0.98‑0.99,可见改进的暗管边
界能够准确预测田间暗管排水过程,可用于指导含黏土夹层地区的暗管埋设,并为研究暗
管排水条件下水盐运移理论进行科学指导。
以做出若干改进和变动,这些改进和变动也视为本发明的保护范围。