CN108595904B - 一种基于水位调控的灌区水田非点源污染负荷估算方法 - Google Patents
一种基于水位调控的灌区水田非点源污染负荷估算方法 Download PDFInfo
- Publication number
- CN108595904B CN108595904B CN201810172010.1A CN201810172010A CN108595904B CN 108595904 B CN108595904 B CN 108595904B CN 201810172010 A CN201810172010 A CN 201810172010A CN 108595904 B CN108595904 B CN 108595904B
- Authority
- CN
- China
- Prior art keywords
- paddy field
- water level
- model
- phosphorus
- nitrogen
- 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.)
- Active
Links
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 title claims abstract description 126
- 238000000034 method Methods 0.000 title claims abstract description 78
- 238000003973 irrigation Methods 0.000 title claims abstract description 29
- 230000002262 irrigation Effects 0.000 title claims abstract description 29
- 230000033228 biological regulation Effects 0.000 title claims abstract description 18
- IJGRMHOSHXDMSA-UHFFFAOYSA-N Atomic nitrogen Chemical compound N#N IJGRMHOSHXDMSA-UHFFFAOYSA-N 0.000 claims abstract description 88
- OAICVXFJPJFONN-UHFFFAOYSA-N Phosphorus Chemical compound [P] OAICVXFJPJFONN-UHFFFAOYSA-N 0.000 claims abstract description 49
- 229910052698 phosphorus Inorganic materials 0.000 claims abstract description 49
- 239000011574 phosphorus Substances 0.000 claims abstract description 49
- 229910052757 nitrogen Inorganic materials 0.000 claims abstract description 44
- 230000008595 infiltration Effects 0.000 claims abstract description 32
- 238000001764 infiltration Methods 0.000 claims abstract description 32
- 238000004364 calculation method Methods 0.000 claims abstract description 31
- 230000008569 process Effects 0.000 claims abstract description 29
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 17
- 238000004162 soil erosion Methods 0.000 claims abstract description 13
- 230000006872 improvement Effects 0.000 claims abstract description 8
- 239000002689 soil Substances 0.000 claims description 33
- 241000209094 Oryza Species 0.000 claims description 29
- 235000007164 Oryza sativa Nutrition 0.000 claims description 29
- 235000009566 rice Nutrition 0.000 claims description 29
- 230000008020 evaporation Effects 0.000 claims description 23
- 238000001704 evaporation Methods 0.000 claims description 23
- 238000003860 storage Methods 0.000 claims description 19
- 230000014759 maintenance of location Effects 0.000 claims description 9
- 230000003628 erosive effect Effects 0.000 claims description 7
- 230000004044 response Effects 0.000 claims description 7
- 239000002352 surface water Substances 0.000 claims description 7
- 230000006870 function Effects 0.000 claims description 6
- 239000002245 particle Substances 0.000 claims description 5
- 238000007654 immersion Methods 0.000 claims description 4
- 230000004048 modification Effects 0.000 claims description 4
- 238000012986 modification Methods 0.000 claims description 4
- 238000005457 optimization Methods 0.000 claims description 4
- 239000002344 surface layer Substances 0.000 claims description 4
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 claims description 3
- 238000010276 construction Methods 0.000 claims description 3
- 239000003621 irrigation water Substances 0.000 claims description 3
- 239000007788 liquid Substances 0.000 claims description 3
- 230000001960 triggered effect Effects 0.000 claims description 3
- MMDJDBSEMBIJBB-UHFFFAOYSA-N [O-][N+]([O-])=O.[O-][N+]([O-])=O.[O-][N+]([O-])=O.[NH6+3] Chemical compound [O-][N+]([O-])=O.[O-][N+]([O-])=O.[O-][N+]([O-])=O.[NH6+3] MMDJDBSEMBIJBB-UHFFFAOYSA-N 0.000 claims description 2
- 238000012821 model calculation Methods 0.000 claims description 2
- 125000001477 organic nitrogen group Chemical group 0.000 claims description 2
- 230000035515 penetration Effects 0.000 claims description 2
- 230000000717 retained effect Effects 0.000 claims description 2
- 238000001179 sorption measurement Methods 0.000 claims description 2
- 239000013589 supplement Substances 0.000 claims 1
- 230000002349 favourable effect Effects 0.000 abstract description 4
- 235000015097 nutrients Nutrition 0.000 abstract description 4
- 238000004088 simulation Methods 0.000 description 14
- 230000008859 change Effects 0.000 description 6
- 238000005516 engineering process Methods 0.000 description 6
- 238000012544 monitoring process Methods 0.000 description 5
- 238000011160 research Methods 0.000 description 3
- 238000005070 sampling Methods 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 2
- 230000008901 benefit Effects 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 238000012795 verification Methods 0.000 description 2
- 239000000853 adhesive Substances 0.000 description 1
- 230000001070 adhesive effect Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000002950 deficient Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 239000003344 environmental pollutant Substances 0.000 description 1
- 239000003337 fertilizer Substances 0.000 description 1
- 239000008187 granular material Substances 0.000 description 1
- 229910001385 heavy metal Inorganic materials 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 238000002844 melting Methods 0.000 description 1
- 230000008018 melting Effects 0.000 description 1
- 239000000575 pesticide Substances 0.000 description 1
- 231100000719 pollutant Toxicity 0.000 description 1
- 238000001556 precipitation Methods 0.000 description 1
- 238000003825 pressing Methods 0.000 description 1
- 238000004445 quantitative analysis Methods 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
- 238000002791 soaking Methods 0.000 description 1
- 231100000331 toxic Toxicity 0.000 description 1
- 230000002588 toxic effect Effects 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
- 238000003911 water pollution Methods 0.000 description 1
Images
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—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/02—Agriculture; Fishing; Forestry; 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
Landscapes
- Business, Economics & Management (AREA)
- Engineering & Computer Science (AREA)
- Human Resources & Organizations (AREA)
- Marketing (AREA)
- Marine Sciences & Fisheries (AREA)
- Mining & Mineral Resources (AREA)
- Agronomy & Crop Science (AREA)
- Health & Medical Sciences (AREA)
- Economics (AREA)
- General Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Animal Husbandry (AREA)
- Primary Health Care (AREA)
- Strategic Management (AREA)
- Tourism & Hospitality (AREA)
- Physics & Mathematics (AREA)
- General Business, Economics & Management (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Cultivation Of Plants (AREA)
- Processing Of Solid Wastes (AREA)
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 (5)
1.一种基于水位调控的灌区水田非点源污染负荷估算方法,其特征在于:利用FORTRAN语言对SWAT模型的壶穴模块进行改良,建立起能够表征水田水循环过程和非点源污染之间关联的水田模块,具体方法步骤如下:
步骤一:壶穴构型改良
利用FORTRAN 77语言对SWAT模型的壶穴模块的源代码进行编译,并据此构建符合水田灌溉和排水特征的水田模型;
步骤二:蒸散发算法的改良
利用编译器对壶穴模块的蒸散发算法进行修正,由原模型提出的基于叶面积指数与最大潜在蒸散发量之间关系的概化算法,改进为根据水田淹没或排空状态,分别计算叶片截留蒸发量、水稻蒸散发量、田面水蒸散发量、土壤蒸散发量的优化算法;
步骤三:下渗量计算方法的改良
利用编译器对壶穴模型的渗透算法进行编译,由原壶穴模型提出的基于土壤最大持水能力和土壤含水量的计算方法,改进为基于水田储水量与水田日稳定下渗量关系的优化算法;
步骤四:灌溉及排水控制模块的构建
基于水稻的实际设置规程,对水位控制单元进行编译;共设置最低三条水位线用于控制水田的排水及灌溉过程:最低水位、最适水位和最高水位;
步骤五:水田土壤侵蚀量计算方法的调整
为保持新水田模块与原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。
2.根据权利要求1所述的一种基于水位调控的灌区水田非点源污染负荷估算方法,其特征在于:所述步骤六的具体过程如下:颗粒态氮、磷的流失量依照Williams and Hann提出的公式(13)计算,溶解态氮磷流失过程的计算方法与原SWAT模型的方法一致;
SBP=0.001Con·(Sed/SA)·ε (13)
式中:SBP表示吸附在颗粒物上的氮、磷含量,kg/ha;Con表示表层10mm深度土壤中的氮磷含量,g/t;ε表示氮、磷在颗粒物上的吸附比率;Sed表示土壤侵蚀量,t。
3.根据权利要求1所述的一种基于水位调控的灌区水田非点源污染负荷估算方法,其特征在于:所述步骤七中,模型所需的基础输入数据包括:数字高程数据DEM、土地利用数据、土壤属性数据及气象数据;此外,壶穴面积的参数为POT_FR可在模型中直接调整,田埂占比和逐日水位数据需生成TXT文件录入。
4.根据权利要求1所述的一种基于水位调控的灌区水田非点源污染负荷估算方法,其特征在于:所述的田埂占比为1%~3%。
5.根据权利要求1所述的一种基于水位调控的灌区水田非点源污染负荷估算方法,其特征在于:所述步骤八中的水田总氮非点源污染负荷为地表流失的有机氮、地表流失的无机氮、土壤侧向流中的无机氮及由地下水进入河道中的无机氮的总负荷;水田的总磷非点源污染负荷为由地表流失的颗粒态无机磷、颗粒态有机磷及溶解态磷的总负荷。
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 CN108595904A (zh) | 2018-09-28 |
CN108595904B true 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) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113283718B (zh) * | 2021-05-12 | 2021-12-03 | 湖北省水利水电科学研究院 | 基于swat改进模型的虾稻共作区域水文模拟方法及系统 |
Citations (1)
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 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8303816B2 (en) * | 2011-08-22 | 2012-11-06 | Modular Wetland Systems, Inc. | Wetland biofilter chamber with peripheral catch basin and method of use thereof |
CN102867120B (zh) * | 2012-09-05 | 2016-02-24 | 环境保护部卫星环境应用中心 | 一种基于遥感像元的非点源污染计算方法 |
CN103020424B (zh) * | 2012-11-22 | 2016-04-27 | 北京师范大学 | 一种基于人工降雨实验的北方平原农田区非点源污染负荷估算方法 |
CN103294923B (zh) * | 2013-06-09 | 2016-03-30 | 北京师范大学 | 水稻田非点源污染流出负荷的观测方法及系统 |
-
2018
- 2018-03-01 CN CN201810172010.1A patent/CN108595904B/zh active Active
Patent Citations (1)
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 |
Also Published As
Publication number | Publication date |
---|---|
CN108595904A (zh) | 2018-09-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Xu et al. | Assessing the groundwater dynamics and impacts of water saving in the Hetao Irrigation District, Yellow River basin | |
Moriasi et al. | Evaluating hydrology of the Soil and Water Assessment Tool (SWAT) with new tile drain equations | |
Dripps et al. | A simple daily soil–water balance model for estimating the spatial and temporal distribution of groundwater recharge in temperate humid areas | |
Panagopoulos et al. | SWAT parameterization for the identification of critical diffuse pollution source areas under data limitations | |
Baffaut et al. | Hydrologic and water quality modeling: Spatial and temporal considerations | |
Dabney et al. | The application of the Revised Universal Soil Loss Equation, Version 2, to evaluate the impacts of alternative climate change scenarios on runoff and sediment yield | |
Knisel et al. | GLEAMS model application on a heavy clay soil in Finland | |
Pomeroy et al. | The cold regions hydrological modelling platform for hydrological diagnosis and prediction based on process understanding | |
Kersebaum et al. | Operational use of agro-meteorological data and GIS to derive site specific nitrogen fertilizer recommendations based on the simulation of soil and crop growth processes | |
Golmohammadi et al. | SWATDRAIN, a new model to simulate the hydrology of agricultural Lands, model development and evaluation | |
Van Dessel et al. | Predicting land cover changes and their impact on the sediment influx in the Lake Balaton catchment | |
Rad et al. | MOD $$ AT: A hydro-economic modeling framework for aquifer management in irrigated agricultural regions | |
Bowman et al. | Impacts of irrigation and drought on Illinois ground-water resources | |
CN108595904B (zh) | 一种基于水位调控的灌区水田非点源污染负荷估算方法 | |
Stoianov et al. | A framework linking urban park land use with pond water quality | |
Nelson et al. | Use of soft data for multicriteria calibration and validation of Agricultural Policy Environmental eXtender: Impact on model simulations | |
Srivastava et al. | GIS based integrated modelling framework for agricultural canal system simulation and management in Indo-Gangetic plains of India | |
De Vos et al. | Raising surface water levels in peat areas with dairy farming: Upscaling hydrological, agronomical and economic effects from farm-scale to local scale | |
Ndulue et al. | DRAINMOD simulation of drain spacing impact on canola yield in heavy clay soils in the Canadian prairies | |
Hack-ten Broeke | Irrigation management for optimizing crop production and nitrate leaching on grassland | |
Ghobadi et al. | Development and application of a seasonal furrow irrigation model (SFIM) | |
Ma et al. | Evapotranspiration predictions: a comparison among GLEAMS, Opus, PRZM-2, and RZWQM models in a humid and thermic climate | |
Dourte et al. | Agricultural management impacts on groundwater: simulations of existing and alternative management options in Peninsular India | |
Wu et al. | BASINMAN: a water balance model for farms with subsurface pipe drainage and an on-farm basin | |
Rajput et al. | Enhancement of field water use efficiency in the Indo‐Gangetic Plain of India |
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 |