CN108595904A - 一种基于水位调控的灌区水田非点源污染负荷估算方法 - Google Patents

一种基于水位调控的灌区水田非点源污染负荷估算方法 Download PDF

Info

Publication number
CN108595904A
CN108595904A CN201810172010.1A CN201810172010A CN108595904A CN 108595904 A CN108595904 A CN 108595904A CN 201810172010 A CN201810172010 A CN 201810172010A CN 108595904 A CN108595904 A CN 108595904A
Authority
CN
China
Prior art keywords
paddy field
water level
indicate
formula
day
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201810172010.1A
Other languages
English (en)
Other versions
CN108595904B (zh
Inventor
欧阳威
魏鹏
谢先红
郝增超
张璇
刘莲华
郝新
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Beijing Normal University
Original Assignee
Beijing Normal University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Beijing Normal University filed Critical Beijing Normal University
Priority to CN201810172010.1A priority Critical patent/CN108595904B/zh
Publication of CN108595904A publication Critical patent/CN108595904A/zh
Application granted granted Critical
Publication of CN108595904B publication Critical patent/CN108595904B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q50/00—Systems or methods specially adapted for specific business sectors, e.g. utilities or tourism
    • G06Q50/02—Agriculture; Fishing; Mining
    • G—PHYSICS
    • G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16Z—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS, NOT OTHERWISE PROVIDED FOR
    • G16Z99/00—Subject matter not provided for in other main groups of this subclass
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2219/00—Indexing scheme relating to application aspects of data processing equipment or methods
    • G06F2219/10—Environmental application, e.g. waste reduction, pollution control, compliance with environmental legislation

Abstract

本发明涉及一种基于水位调控的灌区水田非点源污染负荷估算方法，步骤如下：步骤一：壶穴构型改良；步骤二：蒸散发算法的改良；步骤三：下渗量计算方法的改良；步骤四：灌溉及排水控制模块的构建；步骤五：水田土壤侵蚀量计算方法的调整；步骤六：水田氮、磷流失量的计算方法；步骤七：改进模型的运行；步骤八：计算水田非点源污染负荷。本发明方法，分别考虑了水田在浸没状态下和排空期内的不同水文特征，同时也对各种算法进行了优化改良，有助于得到更为精确的稻田水文循环过程和养分流失过程。

Description

一种基于水位调控的灌区水田非点源污染负荷估算方法
【技术领域】
本发明涉及一种灌区水田非点源污染负荷估算方法。尤其涉及一种基于水位调控的灌区水田非点源污染负荷估算方法，属于农业环境保护技术。
【背景技术】
流域农业非点源污染是指在降水、融雪及农田退水等过程中，部分土壤颗粒物、营养元素及有毒害作用的农药和重金属等污染物随径流从农田进入水体，进而形成的水体污染。而径流作为非点源污染形成的重要驱动因子，在非点源污染生成、运移及转化过程中起着举足轻重的作用。水田作为人工营造的湿地系统，其表层的溢流和排水及地下渗透过程贯穿水稻整个生长过程，同时也势必加剧氮磷等营养元素的流失，威胁脆弱的流域水环境。
传统的田间采样技术虽然能够用于分析田间尺度的水田非点源污染负荷，但该技术成本高，且监测周期较长，在采样区选取方面也具有一定的技术难度。此外，在流域尺度下，非点源污染的输出过程具有普遍的时空异质性规律，因此也限制田间采样法的应用和推广。相较而言，模型技术在流域非点源污染时空序列绘制及风险源区识别方面有着突出的优势。在模型率定和模型验证技术的支持下，模型技术能够直观的对流域非点源污染进行量化分析。但目前主流的流域非点源预测模型(如SWAT，HSPF和AnnAGNPS)在灌区的应用过程中也存在着明显的缺陷：(1)主流流域非点源预测模型并未内置水田模块，仅通过径流曲线法、Green-Ampt入渗模型等将水田简化为旱田处理；(2)未考虑水田的淹没状态，在灌溉发生时，除部分灌溉水按下渗量计算，其余均视为地表径流量；(3)不能模拟水田的水位的年季波动特性，无法在非点源污染与水田水位调控之间建立起联系。
目前的非点源污染的预测方法正在不断发展，但考虑水田实际灌溉特征的预测方法却十分匮乏，例如申请号为CN201410619286.1，名称为一种基于GIS的分布式水田氮素径流流失负荷估算方法，公开号CN104361523A，仅利用GIS系统、土地利用情景和肥料施用情况构建了简单的输出系模型，用于评价水田非点源污染输出负荷，缺乏足够的模型机理支持，也无助与建立水位调控与非点源污染的对应关系。又如申请号为CN201310232100.2，名称为水田非点源污染流出负荷的观测方法及系统，公开号CN103294923A，构建了一种田间尺度的水田非点源污染分析及估算方法，但因成本问题难以在流域尺度展开。再如申请号为CN201210326701.5，名称为一种基于遥感像元的非点源污染计算方法，公开号为CN102867120A，该方法根据土地利用、植被盖度、土壤属性及氮磷平衡规律估算非点源污染释放规律，但该方法难以表征浸没状态下水田的氮磷释放规律。鉴于此种背景，有必要对现有的流域非点源污染预测模型进行修正，使其能够高效、准确的模拟水位变化条件下的水田非点源污染输出负荷。
【发明内容】
1、目的：本发明的目的是提供一种基于水位调控的灌区水田非点源污染负荷估算方法，进而可以模拟水田水位变化条件下水田的非点源污染输出规律，有助于精准、有效的灌区的流域非点源污染进行分析。
2、技术方案：本发明一种基于水位调控的灌区水田非点源污染负荷估算方法，其特征在于利用FORTRAN语言对SWAT模型的壶穴模块进行改良，建立起能够表征水田水循环过程和非点源污染之间关联的水田模块。其具体方法步骤如下：
步骤一：壶穴构型改良
利用FORTRAN 77语言对SWAT模型的壶穴模块的源代码进行编译，并据此构建符合水田灌溉和排水特征的水田模型。本专利对壶穴模型的构型进行了必要改良，将其由原SWAT模型中的圆锥体设定改良为截面积恒定的多边形结构，并基于改进方法增设了田埂面积的补充算法。田埂面积和改进的水田有效储水面积的计算方法如公式1和2所示。
SR＝SHD·FR (1)
SA＝SHD·(1-FR) (2)
式中：SR表示田埂面积，ha；FR表示田埂面积与水文响应单元面积之比，％；SHD表示水田所处的水文响应单元的总面积，ha；SA表示水田有效储水面积，ha。
步骤二：蒸散发算法的改良
利用编译器对壶穴模块的蒸散发算法进行修正，由原模型提出的基于叶面积指数与最大潜在蒸散发量之间关系的概化算法，改进为根据水田淹没或排空状态，分别计算叶片截留蒸发量、水稻蒸散发量、田面水蒸散发量、土壤蒸散发量的优化算法。当叶片截留量大于水田日潜在蒸散发量时，水田模块将先计算叶片截留蒸发量，见公式(3)。当叶片截留量小于水田日潜在蒸散发量时，其差值由水稻蒸散发量弥补，见公式(4)。当水稻蒸散发量及叶片截留量之和仍小于水田日潜在蒸散发量时，水田模块将根据水位条件，分别进行计算：在浸没状态下，蒸发量的差额由田面水蒸散发量补足，见公式(5)；在排空状态下，蒸发量的差额由土壤蒸散发量补足，见公式(6)。
其中，水田日潜在蒸散发量由Penman-Monteith方程计算。水稻蒸散发量可根据稻田的叶面积指数，分别由公式(7)和公式(8)计算。当水田处于排空状态下，其土壤蒸散发规律与旱田相近，因此仍沿用原SWAT模型的计算方法，不进行改动。
Ea＝Ecan Ecan>PET (3)
Ea＝Ep+Ecan Ecan<PET<Ecan+Ep (4)
Ea＝Ep+Ecan+EW PET>Ecan+Ep (5)
Ea＝Ep+Ecan+Es PET>Ecan+Ep (6)
Ep＝(PET-Ecan)·L/3 0≤L≤3 (7)
Ep＝PET-Ecan L>3 (8)
式中：PET表示水田日潜在蒸散发量，mm；Ea表示水田实际蒸发量，mm；Ecan表示叶片截留蒸发量，mm；Ep表示水稻蒸散发量，mm；Ew表示田面水蒸发量，mm；Es表示土壤蒸散发量，mm；L表示叶面积指数。
步骤三：下渗量计算方法的改良
利用编译器对壶穴模型的渗透算法进行编译，由原壶穴模型提出的基于土壤最大持水能力和土壤含水量的计算方法，改进为基于水田储水量与水田日稳定下渗量关系的优化算法。当水田处于浸没状态时，若水田储水量大于水田日稳定下渗量时，水田实际下渗量与水田日稳定下渗量在数值上相等，如公式(9)所示；若水田储水量小于水田日稳定下渗量时，则水田储水量即为水田实际下渗量，如公式(10)。当水田处于排空状态时，其水文循环过程与旱田类似，下渗量的计算方法与原SWAT模型旱田的下渗量计算方式一致，不做修改。
Vinf＝Ks·SA＝Vmax Vs≥Vmax (9)
Vinf＝Vs Vs<Vmax (10)
式中：Vinf表示水田实际下渗量，m3；Ks表示土壤稳定入渗速率，cm/day；SA表示土壤有效储水面积，ha；Vs表示水田储水量，m3；Vmax表示水田日稳定下渗量。
步骤四：灌溉及排水控制模块的构建
基于水稻的实际设置规程，对水位控制单元进行编译。共设置最低三条水位线用于控制水田的排水及灌溉过程：最低水位、最适水位和最高水位。当稻田水位小于最低水位时，水田模块将触发灌溉功能，直至水位达到最适水位，如公式(11)。在降雨事件中，当水田的水位上升至最高水位时，将触发水田的排水功能，多余的田面水将以溢流的形式排入河网，直至水位回落至最高水位以下，如公式(12)。三条水位线均可设定逐日数据，并以txt文本格式保存，供模型读取。
Ri＝(H-Pd,i)·SA Pd,i<Hmin (11)
Oi＝(Pd,i-Hmax)·SA Pd,i>Hmax (12)
式中：Ri表示灌溉用水量，m3；Pd,i表示水田的浸没深度，mm；Oi表示水田的溢流量，m3；Hmin表示最低水位，mm；H表示最适水位，mm；Hmax表示最高水位，mm。
步骤五：水田土壤侵蚀量计算方法的调整
为保持新水田模块与原SWAT模型的参数一致性，仍采用改进通用流失方程作为水田土壤侵蚀量的计算公式；并根据水田的水位条件，对计算方法进行了优化调整。当水田处于浸没状态时，若未发生溢流或排水，则不计算土壤侵蚀过程。若有溢流或排水发生，则仅计算田埂侵蚀量，水田表土由于受到自由液面保护，不计算其侵蚀量。当水田处于排空状态下，水田田埂及田面均按可侵蚀面计算。
步骤六：水田氮、磷流失量的计算方法
氮、磷的流失量按其流失形态分别计算，氮的流失形态包括颗粒态氮和硝态氮，磷的流失形态包括颗粒态有机磷、颗粒态无机磷和溶解态磷。其中颗粒态氮、磷的流失量依照Williams and Hann提出的公式(13)计算，溶解态氮磷流失过程的计算方法与原SWAT模型的方法一致。
SBP＝0.001Con·(Sed/SA)·ε (13)
式中：SBP表示吸附在颗粒物上的氮、磷含量，kg/ha；Con表示表层10mm深度土壤中的氮磷含量，g/t；ε表示氮、磷在颗粒物上的吸附比率；Sed表示土壤侵蚀量，t。
步骤七：改进模型的运行
录入模型所需的基础数据，根据模型率定结果对参数进行调整；设置与改进模型相关的参数：包括壶穴面积、田埂占比、水田逐日水位数据及运行模型。其中，模型所需的基础输入数据包括：数字高程数据DEM、土地利用数据、土壤属性数据及气象数据。此外，壶穴面积的参数为POT_FR可在模型中直接调整，田埂占比和逐日水位数据需生成TXT文件录入。
其中田埂占比以1％～3％为宜。
步骤八：计算水田非点源污染负荷
在模型输出文件中，读取水田的氮磷流失数据，分别计算流域中各水文响应单元的总氮及总磷的非点源污染负荷。
其中，步骤八中的水田总氮非点源污染负荷为地表流失的有机氮、地表流失的无机氮、土壤侧向流中的无机氮及由地下水进入河道中的无机氮的总负荷；水田的总磷非点源污染负荷为由地表流失的颗粒态无机磷、颗粒态有机磷及溶解态磷的总负荷。
3、优点及功效：
本发明所述的一种基于水位调控的灌区水田非点源污染负荷估算方法，以改进SWAT模型为基础，能够实现流域尺度的水田非点源污染负荷的估算，并有助于反映非点源污染的时空异质性特征，也有助于实现稻区的流域水环境的优化管理。
本发明所述的一种基于水位调控的灌区水田非点源污染负荷估算方法，充分考虑了水稻生长周期内的水位变化过程，有助于建立稻田水位控制与非点源污染之间的对应关系，同时也可以为水田的科学管理提供指导和帮助。
本发明所述的一种基于水位调控的灌区水田非点源污染负荷估算方法，分别考虑了水田在浸没状态下和排空期内的不同水文特征，同时也对各种算法进行了优化改良，有助于得到更为精确的稻田水文循环过程和养分流失过程。
【附图说明】
图1为改良水田模型运行流程的示意图。
图2为改良水田模型河道径流率定结果图。
图3为0cm～15cm原模型与改良模型的水田土壤水率定结果比较图。
图4为15cm～30cm原模型与改良模型的水田土壤水率定结果比较图。
图5为改良水田模型年均氮流失负荷模拟值与前人监测结果比较图。
图6为改良水田模型年均磷流失负荷模拟值与前人监测结果比较图。
图7为改良模型水田月均氮流失负荷模拟值的结果图。
图8为改良模型的水田月均磷流失负荷模拟值的结果图。
【具体实施方式】
见图1～图8，下面结合附图和东北三江平原的阿布胶灌区的实际案例，对本发明做进一步说明。
本发明一种基于水位调控的灌区水田非点源污染负荷的计算方法，该方法具体步骤如下：
步骤一：改进壶穴形状设定
基于图1中所示的改进模型运行框架，利用Intel Visual FORTRAN中的FORTRAN77语言对SWAT模型的壶穴模块的源代码进行编译，建立能够符合水田的水文循环及污染输出特征的改进模型单元，即水田模块。对原模型中壶穴构型进行修正，由原SWAT模型中的圆锥体设定修正为截面积恒定的多边形结构，保证水田面积为恒定值，且不随水位变化而出现波动。此外,在原模型的基础上，增加了田埂面积的计算方法。其具体计算方法见公式(1)、公式(2)。
步骤二：改进水田蒸散发计算方法
在保证水田面积恒定的基础上，进一步修改壶穴模型源代码，对其蒸发量的计算方法进行调整。根据稻田的水位变化特征，分别编译排空状态和浸没状态下的两种水田蒸发量计算方法。先利用Penman-Monteith方程计算水田日潜在蒸散发量，此后根据各项蒸发量与水田日潜在蒸散发量的差值，对计算方法不断进行调整。当叶片截留量大于水田日潜在蒸散发量时，水田模块将只计算叶片截留蒸发量，如公式(3)；当叶片截留量小于水田日潜在蒸散发量时，其差值由水稻蒸散发量弥补，如公式(4)，其中水稻蒸散发量的具体计算方法如公式(7)、公式(8)所示。当水稻蒸散发量及叶片截留蒸发量之和仍小于水田日潜在蒸散发量时，水田模块将根据水位条件，分别进行计算：在浸没状态下，蒸发量的差额由田面水蒸发量补足，如公式(5)；在排空状态下，蒸发量的差额由土壤蒸散发量补足，如公式(6)，其中土壤蒸散发的计算方法与原SWAT模型一致。
步骤三：改进水田下渗量的计算方法
在步骤一和步骤二的基础上，进一步改进壶穴模块的代码，构建基于水位变化的水田下渗量计算方法。如公式(9)所示，当水田处于浸没状态时，若水田储水量大于水田日稳定下渗量，则水田实际下渗量等于水田日稳定下渗量；若水田储水量小于水田日稳定下渗量，水田储水量均视为水田实际下渗量，如公式(10)。当水田处于排空状态时，其水田实际下渗量的计算方法与原模型旱田下渗量的计算方式相一致。
步骤四：灌溉及排水控制单元的构建
基于水田的逐日水位控制需求，结合原SWAT模型的灌溉模块，构建三条水位控制线。最低水位为灌溉触发水位，低于此水位时，触发模型灌溉功能，如公式(11)；最适水位为水稻的最适宜生长水位，灌溉操作在触及此水位时停止；最高水位为排水触发水位，在降雨过程中，当水田的水位上升至此高度时，模型将触发排水功能，直至水位回落至最高水位线以下，如公式(12)。三条水位线均需设置逐日的水位值，并分别将对应水位输入至rice_paddy_depth.txt文件中。稻田水位在设置过程中应充分参考水稻的品种及各生长阶段的水位需求，本研究所模拟的水稻品种为三江1号，并根据泡田期、排水期、种植期、复青期、分蘖期、抽穗期、灌浆期、成熟期、晒田期、休耕期的水稻生产技术规程，对水稻生长周期内的水位进行了分别设定。
步骤五：水田土壤侵蚀量计算方法的调整
在完成步骤一至步骤四中水文循环模块的改良的基础上，对水田污染负荷的计算方法进行改良。为保持改进水田模块与原SWAT模型参数的一致性，仍采用改进通用流失方程MUSLE对水田的土壤侵蚀量的进行计算。此处所用方法与原模型一致，只在计算方式上有所调整。当水田处于浸没状态时，田面由于受到液面的保护，不会发生土壤侵蚀，此时，仅考虑田埂部分的侵蚀计算；当水田处于排空状态时，即最高水位设置值为0，其侵蚀量为田面侵蚀量与田埂侵蚀量之和。其中田埂占比在rice_paddy_depth.txt文件中设定，并以水田面积的1％～3％为宜，本实施例中设置为1％。
步骤六：水田氮、磷流失量的计算方法
基于土壤侵蚀量的计算结果，通过设置氮、磷在土壤颗粒上的分配系数，计算水田颗粒态氮磷的流失量，如公式(13)；田面水中溶解态氮、磷的含量通过表层10mm土壤与田面水之间的交换量计算。完成上述所有编译工作后，利用编译器生成可执行的exe文件，并将其命名为swatuser.exe。
步骤七：改进水田模型的运行
运行SWAT模型，并输入DEM数据、河网数据、土地利用数据、土壤属性数据及气象数据，并将土地利用类型为水田的水文响应单元的POT_FR参数值设置为1。完成所有模型基础设置后，将步骤五和步骤七中准备好的rice_paddy_depth.txt文件和swatuser.exe文件复制至swat工程文件下的Textinout文件夹下。完成设置工作后，利用SWAT CUP软件对模型参数进行率定，得到符合流域特征的基本参数。如图2所示，本案例中河道径流的率定和验证结果的决定系数和纳什系数均大于0.6，表明模型具有较好的适用性。在水田土壤水模拟结果方面，在0cm～15cm和15cm～30cm深度土壤范围内，模拟值与监测值的均方根误差分别为14.56mm和24.55mm(图3和图4)，显著优于原模型的模拟结果。在模型的可靠性得到验证之后，用SWATCUP所得的模型参数替换原模型中的参数，并再次运行模型，输出模拟结果。
步骤八：计算水田非点源污染负荷
根据输出结果，分别计算水田的氮、磷年均流失负荷，并结合其他研究中的水田非点源负荷的场地监测值，对模拟结果的合理性进行验证。如图5及图6所示，改进水田模型的氮、磷年均流失负荷的模拟结果均处于置信区间范围内。且从流域水田月均氮、磷流失负荷的模拟结果分析(图7和图8)，改进水田模型的模拟结果也符合案例研究区水田的氮、磷流失规律。其模拟结果显示，氮、磷流失的关键时期为排水和溢流频繁的7月和8月。因此，基于水位调控的灌区水田非点源污染负荷估算方法将有助于准确、快速的对流域水田的非点源污染进行计算。

Claims (9)

1.一种基于水位调控的灌区水田非点源污染负荷估算方法，其特征在于：利用FORTRAN语言对SWAT模型的壶穴模块进行改良，建立起能够表征水田水循环过程和非点源污染之间关联的水田模块，具体方法步骤如下：
步骤一：壶穴构型改良
利用FORTRAN 77语言对SWAT模型的壶穴模块的源代码进行编译，并据此构建符合水田灌溉和排水特征的水田模型；
步骤二：蒸散发算法的改良
利用编译器对壶穴模块的蒸散发算法进行修正，由原模型提出的基于叶面积指数与最大潜在蒸散发量之间关系的概化算法，改进为根据水田淹没或排空状态，分别计算叶片截留蒸发量、水稻蒸散发量、田面水蒸散发量、土壤蒸散发量的优化算法；
步骤三：下渗量计算方法的改良
利用编译器对壶穴模型的渗透算法进行编译，由原壶穴模型提出的基于土壤最大持水能力和土壤含水量的计算方法，改进为基于水田储水量与水田日稳定下渗量关系的优化算法；
步骤四：灌溉及排水控制模块的构建
基于水稻的实际设置规程，对水位控制单元进行编译；共设置最低三条水位线用于控制水田的排水及灌溉过程：最低水位、最适水位和最高水位；
步骤五：水田土壤侵蚀量计算方法的调整
为保持新水田模块与原SWAT模型的参数一致性，仍采用改进通用流失方程作为水田土壤侵蚀量的计算公式；并根据水田的水位条件，对计算方法进行了优化调整；当水田处于浸没状态时，若未发生溢流或排水，则不计算土壤侵蚀过程；若有溢流或排水发生，则仅计算田埂侵蚀量，水田表土由于受到自由液面保护，不计算其侵蚀量；当水田处于排空状态下，水田田埂及田面均按可侵蚀面计算；
步骤六：水田氮、磷流失量的计算方法
氮、磷的流失量按其流失形态分别计算，氮的流失形态包括颗粒态氮和硝态氮，磷的流失形态包括颗粒态有机磷、颗粒态无机磷和溶解态磷；
步骤七：改进模型的运行
录入模型所需的基础数据，根据模型率定结果对参数进行调整；设置与改进模型相关的参数：包括壶穴面积、田埂占比、水田逐日水位数据及运行模型；
步骤八：计算水田非点源污染负荷
在模型输出文件中，读取水田的氮磷流失数据，分别计算流域中各水文响应单元的总氮及总磷的非点源污染负荷。
2.根据权利要求1所述的一种基于水位调控的灌区水田非点源污染负荷估算方法，其特征在于：所述步骤一具体过程如下：对壶穴模型的构型进行改良，将其由原SWAT模型中的圆锥体设定改良为截面积恒定的多边形结构，并基于改进方法增设了田埂面积的补充算法；田埂面积和改进的水田有效储水面积的计算方法如公式(1)和(2)所示：
SR＝SHD·FR (1)
SA＝SHD·(1-FR) (2)
式中：SR表示田埂面积，ha；FR表示田埂面积与水文响应单元面积之比，％；SHD表示水田所处的水文响应单元的总面积，ha；SA表示水田有效储水面积，ha。
3.根据权利要求1所述的一种基于水位调控的灌区水田非点源污染负荷估算方法，其特征在于：所述步骤二具体过程如下：当叶片截留量大于水田日潜在蒸散发量时，水田模块将先计算叶片截留蒸发量，见公式(3)；当叶片截留量小于水田日潜在蒸散发量时，其差值由水稻蒸散发量弥补，见公式(4)；当水稻蒸散发量及叶片截留量之和仍小于水田日潜在蒸散发量时，水田模块将根据水位条件，分别进行计算：在浸没状态下，蒸发量的差额由田面水蒸散发量补足，见公式(5)；在排空状态下，蒸发量的差额由土壤蒸散发量补足，见公式(6)；
其中，水田日潜在蒸散发量由Penman-Monteith方程计算，水稻蒸散发量可根据稻田的叶面积指数，分别由公式(7)和公式(8)计算；当水田处于排空状态下，其土壤蒸散发规律与旱田相近，因此仍沿用原SWAT模型的计算方法，不进行改动；
Ea＝Ecan Ecan>PET (3)
Ea＝Ep+Ecan Ecan<PET<Ecan+Ep (4)
Ea＝Ep+Ecan+EW PET>Ecan+Ep (5)
Ea＝Ep+Ecan+Es PET>Ecan+Ep (6)
Ep＝(PET-Ecan)·L/3 0≤L≤3 (7)
Ep＝PET-Ecan L>3 (8)
式中：PET表示水田日潜在蒸散发量，mm；Ea表示水田实际蒸发量，mm；Ecan表示叶片截留蒸发量，mm；Ep表示水稻蒸散发量，mm；Ew表示田面水蒸发量，mm；Es表示土壤蒸散发量，mm；L表示叶面积指数。
4.根据权利要求1所述的一种基于水位调控的灌区水田非点源污染负荷估算方法，其特征在于：所述步骤三具体过程如下：当水田处于浸没状态时，若水田储水量大于水田日稳定下渗量时，水田实际下渗量与水田日稳定下渗量在数值上相等，如公式(9)所示；若水田储水量小于水田日稳定下渗量时，则水田储水量即为水田实际下渗量，如公式(10)；当水田处于排空状态时，其水文循环过程与旱田类似，下渗量的计算方法与原SWAT模型旱田的下渗量计算方式一致，不做修改；
Vinf＝Ks·SA＝Vmax Vs≥Vmax (9)
Vinf＝Vs Vs<Vmax (10)
式中：Vinf表示水田实际下渗量，m3；Ks表示土壤稳定入渗速率，cm/day；SA表示土壤有效储水面积，ha；Vs表示水田储水量，m3；Vmax表示水田日稳定下渗量。
5.根据权利要求1所述的一种基于水位调控的灌区水田非点源污染负荷估算方法，其特征在于：所述步骤四具体过程如下：当稻田水位小于最低水位时，水田模块将触发灌溉功能，直至水位达到最适水位，如公式(11)；在降雨事件中，当水田的水位上升至最高水位时，将触发水田的排水功能，多余的田面水将以溢流的形式排入河网，直至水位回落至最高水位以下，如公式(12)；三条水位线均可设定逐日数据，并以txt文本格式保存，供模型读取；
Ri＝(H-Pd,i)·SA Pd,i<Hmin (11)
Oi＝(Pd,i-Hmax)·SA Pd,i>Hmax (12)
式中：Ri表示灌溉用水量，m3；Pd,i表示水田的浸没深度，mm；Oi表示水田的溢流量，m3；Hmin表示最低水位，mm；H表示最适水位，mm；Hmax表示最高水位，mm。
6.根据权利要求1所述的一种基于水位调控的灌区水田非点源污染负荷估算方法，其特征在于：所述步骤六的具体过程如下：颗粒态氮、磷的流失量依照Williams and Hann提出的公式(13)计算，溶解态氮磷流失过程的计算方法与原SWAT模型的方法一致；
SBP＝0.001Con·(Sed/SA)·ε (13)
式中：SBP表示吸附在颗粒物上的氮、磷含量，kg/ha；Con表示表层10mm深度土壤中的氮磷含量，g/t；ε表示氮、磷在颗粒物上的吸附比率；Sed表示土壤侵蚀量，t。
7.根据权利要求1所述的一种基于水位调控的灌区水田非点源污染负荷估算方法，其特征在于：所述步骤七中，模型所需的基础输入数据包括：数字高程数据DEM、土地利用数据、土壤属性数据及气象数据；此外，壶穴面积的参数为POT_FR可在模型中直接调整，田埂占比和逐日水位数据需生成TXT文件录入。
8.根据权利要求7所述的一种基于水位调控的灌区水田非点源污染负荷估算方法，其特征在于：所述的田埂占比以1％～3％为宜。
9.根据权利要求1所述的一种基于水位调控的灌区水田非点源污染负荷估算方法，其特征在于：所述步骤八中的水田总氮非点源污染负荷为地表流失的有机氮、地表流失的无机氮、土壤侧向流中的无机氮及由地下水进入河道中的无机氮的总负荷；水田的总磷非点源污染负荷为由地表流失的颗粒态无机磷、颗粒态有机磷及溶解态磷的总负荷。
CN201810172010.1A 2018-03-01 2018-03-01 一种基于水位调控的灌区水田非点源污染负荷估算方法 Active CN108595904B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810172010.1A CN108595904B (zh) 2018-03-01 2018-03-01 一种基于水位调控的灌区水田非点源污染负荷估算方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810172010.1A CN108595904B (zh) 2018-03-01 2018-03-01 一种基于水位调控的灌区水田非点源污染负荷估算方法

Publications (2)

Publication Number Publication Date
CN108595904A true CN108595904A (zh) 2018-09-28
CN108595904B CN108595904B (zh) 2020-05-05

Family

ID=63625535

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810172010.1A Active CN108595904B (zh) 2018-03-01 2018-03-01 一种基于水位调控的灌区水田非点源污染负荷估算方法

Country Status (1)

Country Link
CN (1) CN108595904B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113283718A (zh) * 2021-05-12 2021-08-20 湖北省水利水电科学研究院 基于swat改进模型的虾稻共作区域水文模拟方法及系统

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5928403A (en) * 1993-09-30 1999-07-27 Trustees Of University Of Arkansas Methods for decreasing non-point source pollution from poultry manure
CN102867120A (zh) * 2012-09-05 2013-01-09 环境保护部卫星环境应用中心 一种基于遥感像元的非点源污染计算方法
CN103020424A (zh) * 2012-11-22 2013-04-03 北京师范大学 一种基于人工降雨实验的北方平原农田区非点源污染负荷估算方法
CN103294923A (zh) * 2013-06-09 2013-09-11 北京师范大学 水稻田非点源污染流出负荷的观测方法及系统
US20140251928A1 (en) * 2011-08-22 2014-09-11 Modular Wetland Systems, Inc. Horizontal flow biofilter system and method of use thereof

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5928403A (en) * 1993-09-30 1999-07-27 Trustees Of University Of Arkansas Methods for decreasing non-point source pollution from poultry manure
US20140251928A1 (en) * 2011-08-22 2014-09-11 Modular Wetland Systems, Inc. Horizontal flow biofilter system and method of use thereof
CN102867120A (zh) * 2012-09-05 2013-01-09 环境保护部卫星环境应用中心 一种基于遥感像元的非点源污染计算方法
CN103020424A (zh) * 2012-11-22 2013-04-03 北京师范大学 一种基于人工降雨实验的北方平原农田区非点源污染负荷估算方法
CN103294923A (zh) * 2013-06-09 2013-09-11 北京师范大学 水稻田非点源污染流出负荷的观测方法及系统

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
GANGSHENG WANG ET AL: "Improvement of SWAT2000 modelling to assess the impact of dams and sluices on streamflow in the Huai River basin of China", 《HYDROLOGICAL PROCESSES》 *
GREY R. EVENSON ET AL: "Geographically isolated wetlands and watershed hydrology: A modified model analysis", 《JOURNAL OF HYDROLOGY》 *
J.G.ARNOLD ET AL: "SWAT2000:current capabilities and research opportunities in applied watershed modelling", 《HYDROLOGICAL PROCESSES》 *
欧阳威等: "基于SWAT模型的平原灌区水量平衡模拟研究", 《灌溉排水学报》 *
胡媛: "基于SWAT模型的东圳库区非点源污染模拟研究", 《中国优秀硕士学位论文全文数据库 农业科技辑》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113283718A (zh) * 2021-05-12 2021-08-20 湖北省水利水电科学研究院 基于swat改进模型的虾稻共作区域水文模拟方法及系统
CN113283718B (zh) * 2021-05-12 2021-12-03 湖北省水利水电科学研究院 基于swat改进模型的虾稻共作区域水文模拟方法及系统

Also Published As

Publication number Publication date
CN108595904B (zh) 2020-05-05

Similar Documents

Publication Publication Date Title
Sakaguchi et al. Development and evaluation of a paddy module for improving hydrological simulation in SWAT
Sahu et al. Modeling the effects of riparian buffer zone and contour strips on stream water quality
Zipper et al. Quantifying indirect groundwater-mediated effects of urbanization on agroecosystem productivity using MODFLOW-AgroIBIS (MAGI), a complete critical zone model
Jung et al. Water management practices and SCS curve numbers of paddy fields equipped with surface drainage pipes
Mathur et al. Modeling water uptake by plant roots
Wu et al. Improvement and testing of SWAT for multi-source irrigation systems with paddy rice
Tu et al. Water budget components estimation for a mature citrus orchard of southern China based on HYDRUS-1D model
Ren et al. Hydrological complexities in irrigated agro-ecosystems with fragmented land cover types and shallow groundwater: Insights from a distributed hydrological modeling method
Pereira et al. Surface irrigation
Zhang et al. Experimental analysis of soil moisture response to rainfall in a typical grassland hillslope under different vegetation treatments
Wayne Skaggs Drainage
Holt et al. Improved water and economic sustainability with low-input compact bed plasticulture and precision irrigation
Rahbeh et al. Numerical and conceptual evaluation of preferential flow in Zarqa River Basin, Jordan
CN108595904A (zh) 一种基于水位调控的灌区水田非点源污染负荷估算方法
Jovanovic et al. Total evaporation estimates from a Renosterveld and dryland wheat/fallow surface at the Voëlvlei Nature Reserve (South Africa)
Awe et al. Efficiency of small pan evaporimeter in monitoring evapotranspiration under poly-covered house and open-field conditions in a hot, tropical region of Nigeria
Theiveyanathan et al. An irrigation-scheduling model for application of saline water to tree plantations
Mao et al. Assessment and water saving issues for Ningxia paddies, upper Yellow River Basin
Campos et al. Water saving in the Yellow River basin, China. 1. Irrigation demand scheduling
Napoli et al. Modeling soil and nutrient runoff yields from an Italian vineyard using SWAT
Hajare et al. New technique for evaluation of crop water requirement
Bakhsh et al. Simulating tillage effects on nonpoint source pollution from agricultural lands using GLEAMS
Dabhi et al. Application of Dual Crop Coefficient Approach for estimation of crop water requirement for summer sesame using SIMDualKc Model
Carbon et al. Simulation of the water balance for plants growing on coarse textured soils
Rojano-Cruz et al. Impacts of a hydroinfiltrator rainwater harvesting system on soil moisture regime and groundwater distribution for olive groves in semi-arid Mediterranean regions

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant