CN113139354A - 一种大流域尺度水氮迁移耦合模拟方法 - Google Patents

一种大流域尺度水氮迁移耦合模拟方法 Download PDF

Info

Publication number
CN113139354A
CN113139354A CN202110517237.7A CN202110517237A CN113139354A CN 113139354 A CN113139354 A CN 113139354A CN 202110517237 A CN202110517237 A CN 202110517237A CN 113139354 A CN113139354 A CN 113139354A
Authority
CN
China
Prior art keywords
nitrogen
soil
amount
river
water
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
CN202110517237.7A
Other languages
English (en)
Other versions
CN113139354B (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.)
China Institute of Water Resources and Hydropower Research
Original Assignee
China Institute of Water Resources and Hydropower Research
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 China Institute of Water Resources and Hydropower Research filed Critical China Institute of Water Resources and Hydropower Research
Priority to CN202110517237.7A priority Critical patent/CN113139354B/zh
Publication of CN113139354A publication Critical patent/CN113139354A/zh
Application granted granted Critical
Publication of CN113139354B publication Critical patent/CN113139354B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G06F30/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2113/00—Details relating to the application field
    • G06F2113/08—Fluids
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Mathematical Physics (AREA)
  • Fluid Mechanics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Computing Systems (AREA)
  • Pure & Applied Mathematics (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • Algebra (AREA)
  • Purification Treatments By Anaerobic Or Anaerobic And Aerobic Bacteria Or Animals (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明公开了一种大流域尺度水氮迁移耦合模拟方法，包括S1、根据地表水文过程，结合地表氮素总量，分别计算由于氮沉降和施肥沉降到地表的随地表径流流失、下渗到土壤以及储存在地表洼地的溶解态氮分量；S2、以土壤水文过程为基础，计算源汇项，并根据源汇项利用对流弥散方程计算土壤各层氮素浓度；S3、根据地下水水文过程，利用质量守恒方程计算地下水氮素浓度；S4、根据坡面沟水文过程，利用对流弥散方程计算坡面沟中氮素浓度；S5、根据河道水文过程，利用对流弥散方程计算河道断面氮素浓度。优点是：提高了水文过程模拟的准确性，并且可以很好的模拟冻土融化期河道断面流量和污染物氮素的浓度。

Description

一种大流域尺度水氮迁移耦合模拟方法
技术领域
本发明涉及寒区氮循环技术领域，尤其涉及一种大流域尺度水氮迁移耦合模拟方法。
背景技术
由于人口数量的增长以及工业化的加快发展，氮肥施用量、畜禽养殖量、城镇农村生活污水量等也在逐渐增长并且通过人类活动进入到生态系统，破坏了原有生态系统的平衡，造成土壤污染、水体富营养化等问题。在寒区土壤冻结过程中会有大量的面源氮污染物持留，在融冻过程中伴随着冻土和冰雪融化过程以及春季降水进入地表和地下水体，导致水质迅速下降，造成寒冷地区面源污染物流出具有累积性和突发性。
此外，国内主流的流域水文模型之一的大尺度流域水与能量转化过程模拟(Waterand Energy tranfer Processes in Large river basins，简称WEP-L)可以很好的模拟冻融过程中地表、土壤、地下水、坡面沟以及河道中的水分迁移；该模型为水文模型，未考虑土壤中污染物的迁移过程。
发明内容
本发明的目的在于提供一种大流域尺度水氮迁移耦合模拟方法，从而解决现有技术中存在的前述问题。
为了实现上述目的，本发明采用的技术方案如下：
一种大流域尺度水氮迁移耦合模拟方法，包括如下步骤，
S1、根据地表水文过程，结合地表氮素总量，分别计算由于氮沉降和施肥沉降到地表的随地表径流流失、下渗到土壤以及储存在地表洼地的溶解态氮分量；
S2、以土壤水文过程为基础，计算源汇项，并根据源汇项利用对流弥散方程计算土壤各层氮素浓度；
S3、根据地下水水文过程，利用质量守恒方程计算地下水氮素浓度；
S4、根据坡面沟水文过程，利用对流弥散方程计算坡面沟中氮素浓度；
S5、根据河道水文过程，利用对流弥散方程计算河道断面氮素浓度。
优选的，步骤S1具体为，根据地表水文过程，计算由于氮沉降和施肥沉积到地表的氮素随地表径流流失、下渗到土壤以及储存在地表洼地的比例，并利用比例乘以地表氮素总量分别计算随地表径流流失、下渗到土壤以及储存在地表洼地的溶解态氮量；计算公式分别如下，
Figure BDA0003062767340000021
Figure BDA0003062767340000022
Figure BDA0003062767340000023
Min＝Mfer+Mwet
Mrun＝raterun·Min
Minf＝rateinf·Min
Mstore＝ratestore·Min
其中，Rsuf_C为地表径流量；inf_c为地表下渗到土壤的水量；Rstore_C为洼地储留量；raterun为地表径流流失的水量占地表总水量的比例；rateinf为下渗到土壤的水量占地表总水量的比例；ratestore为储存在地表的水量占地表总水量的比例；Min为地表氮的输入总量；Mfer为地表施肥量；Mwet为地表氮沉降量；Mrun为随地表径流流失的溶解态氮量；Minf为下渗到土壤的溶解态氮量；Mstore为储存在地表的溶解态氮量。
优选的，步骤S2具体为，根据土壤水文过程，计算源汇相，源汇项包括由于土壤吸力引起的上下迁移量、氮素转化量以及通过壤中流流失量；进而利用源汇项使用对流弥散方程计算土壤各层氮素的浓度；
利用下式计算在dt时间内的源汇项
Figure BDA0003062767340000024
Figure BDA0003062767340000025
其中，V为土壤体积；Mfer为施氮量；Mwet为每日氮沉降量；Mupdown为由于土壤吸引力引起的氮素在相邻两层的迁移量；Mtran为氮素不同成分之间的转化量，Mtran＝-μiCi，μi为各种反应系数，Ci表示土壤中污染物的浓度；Muptake为植物根系吸收量；Mrunoff为壤中流损失量；
在冻融期，施氮量和植物根系吸收量为零；土壤分为11层，只在第一层土壤考虑氮沉降与施肥，层与层之间考虑由于重力势引起的下渗，由于吸引力引起的土壤水分在上下层的迁移量，进而引起氮素在上下层的迁移；
利用对流弥散方程计算土壤各层氮素浓度的计算公式如下，
Figure BDA0003062767340000031
其中，Cdis为土壤水中溶解态氮的浓度；Csorb为吸附态氮的浓度，Csorb＝Kd·Cdis，Kd表示固相、液相分配系数；·为土壤容重；v为土壤水流通量；z为表示各层土壤相对于基准面的深度；θ为t时刻位于深度z处的土壤体积含水率；D为水动力弥散系数。
优选的，步骤S3具体为，根据地下水水文过程，统计土壤下渗到地下水的氮量、该等高带地下水与周围等高带地下水氮素交换量、地下水与河道的交换量以及地下水氮素转化量，并利用质量守恒方程计算地下水氮素浓度；部分溶解态有机氮、氨氮和硝态氮会进一步淋溶到地下水中，其在地下水中的迁移转化过程采用质量守恒方程描述，
Figure BDA0003062767340000032
其中，P为饱和含水率；h为饱和带的厚度；qr为非饱和土壤带向饱和土壤带的渗漏速率；Cri为非饱和带渗漏水中溶解态氮的浓度；qriver为地下水与河道交换水量；qdgd为地下水与周围等高带地下水的交换水量。
优选的，步骤S4具体为，根据坡面沟水文过程，统计由地表径流以及壤中流进入坡面沟的氮量以及坡面沟中氮素转化量，利用对流弥散方程计算坡面沟中氮素浓度；除对流扩散过程之外还包括旁侧入流、一个等高带的坡面沟到下一个等高带的坡面沟以及不同组分之间的相互转化作用，具体计算公式为，
Figure BDA0003062767340000033
其中，C′dis为河道中溶解态氮的浓度；v为水流的流速；Ex为扩散系数；S为有机氮、氨氮或硝态氮的源汇项，用于表示旁侧入流以及不同组分之间的相互转化作用，具体指代哪一种氮素的源汇项，需要根据计算的氮素种类进行确定。
优选的，对于有机氮，矿化作用会使有机氮转化为氨氮从而减少，有机氮在坡面沟中或河道中的源汇项可表示为，
Figure BDA0003062767340000041
其中，M1为有机氮的旁侧入流；dx为河流的长度；A为河流的横截面积；K1为有机氮的矿化速率；C1为河流中有机氮的浓度；
对于氨氮，矿化作用会使有机氮转化为氨氮从而使河流中氨氮含量增加，硝化作用会使氨氮转化为硝态氮从而使河流中氨氮含量降低；则氨氮在坡面沟中或河道中的源汇项可表示为：
Figure BDA0003062767340000042
其中，M2表示氨氮的旁侧入流；K2表示氨氮的硝化速率；C2表示河流中氨氮的浓度；
对于硝态氮，硝化作用会使氨氮转化为硝态氮从而使河流中硝态氮含量增大，反硝化作用可使硝态氮转化为N2O，使河流中硝态氮含量降低；则硝态氮在坡面沟中或河道中的源汇项可表示为：
Figure BDA0003062767340000043
式中，M3表示硝态氮的旁侧入流；K3表示硝态氮的反硝化速率；C3表示河流中硝态氮的浓度。
优选的，步骤S5具体为，根据河道水文过程，统计坡面沟流入河道的量以及地下水与河道交换量，计算河道中氮素转化量，利用对流弥散方程计算河道断面氮素浓度；河道断面氮素浓度的计算方式与步骤S4相同；溶解态有机氮、氨氮和硝态氮在坡面沟与河道中的迁移转化过程均使用对流弥散方程进行描述，除对流扩散过程之外还包括旁侧入流、上一个子流域的河道到下一个子流域的河道以及不同组分之间的相互转化作用。
优选的，若坡面沟或河道窄而直，则Ex＝5.93·d·v，d表示河流深度；
若坡面沟或河道宽大，则Ex＝0.11·v·w2/d，w表示坡面沟或河道宽度。
优选的，河道中有机氮、氨氮以及硝态氮的旁侧入流包括地下水与河道的交换量以及该子流域各等高带的坡面沟汇入河道的氮量。
本发明的有益效果是：1、考虑了地表、土壤、地下水、坡面沟与河道中氮素的迁移转化过程，溶解态氮随水移动，考虑了随各个阶段氮素的迁移。2、使用对流弥散方程计算冻融过程中土壤、坡面沟、河道污染物的迁移过程，提高了水文过程模拟的准确性，并且可以很好的模拟冻土融化期河道断面流量和污染物氮素的浓度。
附图说明
图1是本发明实施例中模拟方法的流程图；
图2是本发明实施例中地表氮迁移转化模拟流程图；
图3是本发明实施例中土壤中氮迁移转化模拟流程图；
图4是本发明实施例中地下水中氮迁移转化模拟流程图；
图5是本发明实施例中坡面沟中氮迁移转化模拟流程图；
图6是本发明实施例中河道中氮迁移转化模拟流程图；
图7是本发明实施例中流域硝态氮模拟与实测对比图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白，以下结合附图，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施方式仅仅用以解释本发明，并不用于限定本发明。
如图1至图7所示，本实施例中，提供了一种大流域尺度水氮迁移耦合模拟方法，包括如下步骤，
S1、根据地表水文过程，结合地表氮素总量，分别计算由于氮沉降和施肥沉降到地表的随地表径流流失、下渗到土壤以及储存在地表洼地的溶解态氮分量；
S2、以土壤水文过程为基础，计算源汇项，并根据源汇项利用对流弥散方程计算土壤各层氮素的浓度；
S3、根据地下水水文过程，利用质量守恒方程计算地下水氮素浓度；
S4、根据坡面沟水文过程，利用对流弥散方程计算坡面沟中氮素浓度；
S5、根据河道水文过程，利用对流弥散方程计算河道断面氮素浓度。
本实施例中，模拟方法主要包括五部分内容，具体为，地表氮迁移转化模拟、土壤中氮迁移转化模拟、地下水中氮迁移转化模拟、坡面沟中氮迁移转化模拟和河道中氮迁移转化模拟，下面分别针对这五部分内容进行详细的解释说明。
一、地表氮迁移转化模拟
该部分内容主要对应步骤S1，如图2所示，步骤S1具体为，根据地表水文过程，计算由于氮沉降和施肥沉积到地表的氮素随地表径流流失、下渗到土壤以及储存在地表洼地的比例，并利用比例乘以地表氮素总量分别计算随地表径流流失、下渗到土壤以及储存在地表洼地的溶解态氮量；计算公式分别如下，
Figure BDA0003062767340000061
Figure BDA0003062767340000062
Figure BDA0003062767340000063
Min＝Mfer+Mwet
Mrun＝raterun·Min
Minf＝rateinf·Min
Mstore＝ratestore·Min
其中，Rsuf_C为地表径流量，m；inf_c为地表下渗到土壤的水量，m；Rstore_C为洼地储留量，m；raterun为地表径流流失的水量占地表总水量的比例；rateinf为下渗到土壤的水量占地表总水量的比例；ratestore为储存在地表的水量占地表总水量的比例；Min为地表氮的输入总量，kg；Mfer为地表施肥量，kg；Mwet为地表氮沉降量(该地表氮沉降量主要是氨氮和硝态氮的沉降量)，kg；Mrun为随地表径流流失的溶解态氮量，kg；Minf为下渗到土壤的溶解态氮量，kg；Mstore为储存在地表的溶解态氮量，kg。
二、土壤中氮迁移转化模拟
该部分内容对应步骤S2，如图3所示，步骤S2具体为，根据土壤水文过程，计算源汇相，源汇项包括由于土壤吸力引起的上下迁移量、使用一阶动力学反应方程计算氮素转化量以及通过壤中流流失量，在冻融期，施氮量以及植物根系吸收量为0；进而利用源汇项使用对流弥散方程计算土壤各层氮素的浓度；其中源汇项和土壤各层氮素的浓度分别计算如下：
1、利用下式计算在dt时间内的源汇项
Figure BDA0003062767340000071
Figure BDA0003062767340000072
其中，V为土壤体积，m3；Mfer为施氮量，kg；Mwet为每日氮沉降量，kg；Mupdown为由于土壤吸引力引起的氮素在相邻两层的迁移量，kg；Mtran为氮素不同成分之间的转化量(主要包括矿化作用、硝化作用和反硝化作用，可用一阶动力学反应方程进行计算)，kg，Mtran＝-μiCi，μi为各种反应系数，Ci表示土壤中污染物的浓度；Muptake为植物根系吸收量，kg；Mrunoff为壤中流损失量，kg；
土壤分为11层，只在第一层土壤考虑氮沉降与施肥，层与层之间考虑由于重力势引起的下渗，由于吸引力引起的土壤水分在上下层的迁移量，进而引起氮素在上下层的迁移。
2、利用对流弥散方程计算土壤各层氮素浓度(在非饱和带土壤，污染物的垂向运移主要采用对流-弥散方程进行计算)的计算公式如下，
Figure BDA0003062767340000073
其中，Cdis为土壤水中溶解态氮的浓度，kg/m3；Csorb为吸附态氮的浓度，kg/kg，Csorb＝Kd·Cdis，Kd表示固相、液相分配系数，m3/kg；ρ为土壤容重，kg/m3；v为土壤水流通量，m/d；z为表示各层土壤相对于基准面的深度，m；θ为t时刻位于深度z处的土壤体积含水率，m3/m3；D为水动力弥散系数，m2/d。
三、地下水中氮迁移转化模拟
该部分内容对应步骤S3，如图4所示，步骤S3具体为，根据地下水水文过程，统计土壤下渗到地下水的氮量、该等高带地下水与周围等高带地下水氮素交换量、地下水与河道的交换量以及地下水氮素转化量，并利用质量守恒方程计算地下水氮素浓度；
部分溶解态有机氮、氨氮和硝态氮会进一步淋溶到地下水中，其在地下水中的迁移转化过程采用质量守恒方程描述，
Figure BDA0003062767340000081
其中，P为饱和含水率，m3/m3；h为饱和带的厚度，m；qr为非饱和土壤带向饱和土壤带的渗漏速率，m/d；Cri为非饱和带渗漏水中溶解态氮的浓度，kg/m3；qriver为地下水与河道交换水量，m/d；qdgd为地下水与周围等高带地下水的交换水量，m/d。
四、坡面沟中氮迁移转化模拟
该部分对应步骤S4的内容，如图5所示，S4具体为，根据坡面沟水文过程，统计由地表径流以及壤中流进入坡面沟的氮量以及坡面沟中氮素转化量，利用对流弥散方程计算坡面沟中氮素浓度；溶解态有机氮、氨氮和硝态氮在坡面沟中的迁移转化过程使用对流弥散方程进行描述，除对流扩散过程之外还包括旁侧入流、一个等高带的坡面沟到下一个等高带的坡面沟以及不同组分之间的相互转化作用，具体计算公式为，
Figure BDA0003062767340000082
其中，C′dis为河道中溶解态氮的浓度，kg/m3；v为水流的流速，m/d；Ex为扩散系数，m2/d；S为有机氮、氨氮或硝态氮的源汇项，kg/m3/d，用于表示旁侧入流以及不同组分之间的相互转化作用，具体指代哪一种氮素的源汇项，需要根据计算的氮素种类进行确定。
若坡面沟窄而直，则Ex＝5.93·d·v，d表示河流深度，m；
若坡面沟宽大，则Ex＝0.11·v·w2/d，w表示坡面沟宽度，m。
对于有机氮，矿化作用会使有机氮转化为氨氮从而减少，有机氮在坡面沟中或河道中的源汇项可表示为，
Figure BDA0003062767340000083
其中，M1为有机氮的旁侧入流，kg；dx为河流的长度，m；A为河流的横截面积，m2；K1为有机氮的矿化速率，d-1；C1为河流中有机氮的浓度，kg/m3；
对于氨氮，矿化作用会使有机氮转化为氨氮从而使河流中氨氮含量增加，硝化作用会使氨氮转化为硝态氮从而使河流中氨氮含量降低；则氨氮在坡面沟中或河道中的源汇项可表示为：
Figure BDA0003062767340000091
其中，M2表示氨氮的旁侧入流，kg；K2表示氨氮的硝化速率，d-1；C2表示河流中氨氮的浓度，kg/m3；
对于硝态氮，硝化作用会使氨氮转化为硝态氮从而使河流中硝态氮含量增大，反硝化作用可使硝态氮转化为N2O，使河流中硝态氮含量降低；则硝态氮在坡面沟中或河道中的源汇项可表示为：
Figure BDA0003062767340000092
式中，M3表示硝态氮的旁侧入流，kg；K3表示硝态氮的反硝化速率，d-1；C3表示河流中硝态氮的浓度，kg/m3。
五、河道中氮迁移转化模拟
该部分内容对应步骤S5的内容，如图6所示，S5具体为，根据河道水文过程，统计坡面沟流入河道的量以及地下水与河道交换量，计算河道中氮素转化量，利用对流弥散方程计算河道断面氮素浓度；其计算方法与步骤S4的计算方式相同。溶解态有机氮、氨氮和硝态氮在坡面沟和河道中的迁移转化过程均使用对流弥散方程进行描述，除对流扩散过程之外还包括旁侧入流、上一个子流域的河道到下一个子流域的河道以及不同组分之间的相互转化作用。
若河道窄而直，则Ex＝5.93·d·v，d表示河流深度，m；
若河道宽大，则Ex＝0.11·v·w2/d，w表示河道宽度，m。
河道中有机氮、氨氮以及硝态氮的旁侧入流包括地下水与河道的交换量、该子流域各等高带的坡面沟汇入河道的氮量。
通过采用本发明公开的上述技术方案，得到了如下有益的效果：
本发明提供了一种大流域尺度水氮迁移耦合模拟方法，考虑了地表、土壤、地下水、坡面沟与河道中氮素的迁移转化过程，溶解态氮随水移动，考虑了随各个阶段氮素的迁移。使用对流弥散方程计算冻融过程中土壤、坡面沟、河道污染物的迁移过程，提高了水文过程模拟的准确性，并且可以很好的模拟冻土融化期河道断面流量和污染物氮素的浓度。
以上所述仅是本发明的优选实施方式，应当指出，对于本技术领域的普通技术人员来说，在不脱离本发明原理的前提下，还可以做出若干改进和润饰，这些改进和润饰也应视本发明的保护范围。

Claims (9)

1.一种大流域尺度水氮迁移耦合模拟方法，其特征在于：包括如下步骤，
S1、根据地表水文过程，结合地表氮素总量，分别计算由于氮沉降和施肥沉降到地表的随地表径流流失、下渗到土壤以及储存在地表洼地的溶解态氮分量；
S2、以土壤水文过程为基础，计算源汇项，并根据源汇项利用对流弥散方程计算土壤各层氮素浓度；
S3、根据地下水水文过程，利用质量守恒方程计算地下水氮素浓度；
S4、根据坡面沟水文过程，利用对流弥散方程计算坡面沟中氮素浓度；
S5、根据河道水文过程，利用对流弥散方程计算河道断面氮素浓度。
2.根据权利要求1所述的大流域尺度水氮迁移耦合模拟方法，其特征在于：步骤S1具体为，根据地表水文过程，计算由于氮沉降和施肥沉积到地表的氮素随地表径流流失、下渗到土壤以及储存在地表洼地的比例，并利用比例乘以地表氮素总量分别计算随地表径流流失、下渗到土壤以及储存在地表洼地的溶解态氮量；计算公式分别如下，
Figure FDA0003062767330000011
Figure FDA0003062767330000012
Figure FDA0003062767330000013
Min＝Mfer+Mwet
Mrun＝raterun·Min
Minf＝rateinf·Min
Mstore＝ratestore·Min
其中，Rsuf_C为地表径流量；inf_c为地表下渗到土壤的水量；Rstore_C为洼地储留量；raterun为地表径流流失的水量占地表总水量的比例；rateinf为下渗到土壤的水量占地表总水量的比例；ratestore为储存在地表的水量占地表总水量的比例；Min为地表氮的输入总量；Mfer为地表施肥量；Mwet为地表氮沉降量；Mrun为随地表径流流失的溶解态氮量；Minf为下渗到土壤的溶解态氮量；Mstore为储存在地表的溶解态氮量。
3.根据权利要求2所述的大流域尺度水氮迁移耦合模拟方法，其特征在于：步骤S2具体为，根据土壤水文过程，计算源汇相，源汇项包括由于土壤吸力引起的上下迁移量、氮素转化量以及通过壤中流流失量；进而利用源汇项使用对流弥散方程计算土壤各层氮素的浓度；
利用下式计算在dt时间内的源汇项
Figure FDA0003062767330000021
Figure FDA0003062767330000022
其中，V为土壤体积；Mfer为施氮量；Nwet为每日氮沉降量；Mupdown为由于土壤吸引力引起的氮素在相邻两层的迁移量；Mtran为氮素不同成分之间的转化量，Mtran＝-μiCi，μi为各种反应系数，Ci表示土壤中污染物的浓度；Muptake为植物根系吸收量；Mrunoff为壤中流损失量；
在冻融期，施氮量和植物根系吸收量为零；土壤分为11层，只在第一层土壤考虑氮沉降与施肥，层与层之间考虑由于重力势引起的下渗，由于吸引力引起的土壤水分在上下层的迁移量，进而引起氮素在上下层的迁移；
利用对流弥散方程计算土壤各层氮素浓度的计算公式如下，
Figure FDA0003062767330000023
其中，Cdis为土壤水中溶解态氮的浓度；Csorb为吸附态氮的浓度，Csorb＝Kd·Cdis，Kd表示固相、液相分配系数；ρ为土壤容重；v为土壤水流通量；z为表示各层土壤相对于基准面的深度；θ为t时刻位于深度z处的土壤体积含水率；D为水动力弥散系数。
4.根据权利要求3所述的大流域尺度水氮迁移耦合模拟方法，其特征在于：步骤S3具体为，根据地下水水文过程，统计土壤下渗到地下水的氮量、该等高带地下水与周围等高带地下水氮素交换量、地下水与河道的交换量以及地下水氮素转化量，并利用质量守恒方程计算地下水氮素浓度；部分溶解态有机氮、氨氮和硝态氮会进一步淋溶到地下水中，其在地下水中的迁移转化过程采用质量守恒方程描述，
Figure FDA0003062767330000031
其中，P为饱和含水率；h为饱和带的厚度；qr为非饱和土壤带向饱和土壤带的渗漏速率；Cri为非饱和带渗漏水中溶解态氮的浓度；qriver为地下水与河道交换水量；qdgd为地下水与周围等高带地下水的交换水量。
5.根据权利要求4所述的大流域尺度水氮迁移耦合模拟方法，其特征在于：步骤S4具体为，根据坡面沟水文过程，统计由地表径流以及壤中流进入坡面沟的氮量以及坡面沟中氮素转化量，利用对流弥散方程计算坡面沟中氮素浓度；除对流扩散过程之外还包括旁侧入流、一个等高带的坡面沟到下一个等高带的坡面沟以及不同组分之间的相互转化作用，具体计算公式为，
Figure FDA0003062767330000032
其中，C′dis为河道中溶解态氮的浓度；v为水流的流速；Ex为扩散系数；S为有机氮、氨氮或硝态氮的源汇项，用于表示旁侧入流以及不同组分之间的相互转化作用，具体指代哪一种氮素的源汇项，需要根据计算的氮素种类进行确定。
6.根据权利要求5所述的大流域尺度水氮迁移耦合模拟方法，其特征在于：
对于有机氮，矿化作用会使有机氮转化为氨氮从而减少，有机氮在坡面沟中或河道中的源汇项可表示为，
Figure FDA0003062767330000033
其中，M1为有机氮的旁侧入流；dx为河流的长度；A为河流的横截面积；K1为有机氮的矿化速率；C1为河流中有机氮的浓度；
对于氨氮，矿化作用会使有机氮转化为氨氮从而使河流中氨氮含量增加，硝化作用会使氨氮转化为硝态氮从而使河流中氨氮含量降低；则氨氮在坡面沟中或河道中的源汇项可表示为：
Figure FDA0003062767330000034
其中，M2表示氨氮的旁侧入流；K2表示氨氮的硝化速率；C2表示河流中氨氮的浓度；
对于硝态氮，硝化作用会使氨氮转化为硝态氮从而使河流中硝态氮含量增大，反硝化作用可使硝态氮转化为N2O，使河流中硝态氮含量降低；则硝态氮在坡面沟中或河道中的源 汇项可表示为：
Figure FDA0003062767330000041
式中，M3表示硝态氮的旁侧入流；K3表示硝态氮的反硝化速率；C3表示河流中硝态氮的浓度。
7.根据权利要求5所述的大流域尺度水氮迁移耦合模拟方法，其特征在于：步骤S5具体为，根据河道水文过程，统计坡面沟流入河道的量以及地下水与河道交换量，计算河道中氮素转化量，利用对流弥散方程计算河道断面氮素浓度；河道断面氮素浓度的计算方式与步骤S4相同；溶解态有机氮、氨氮和硝态氮在坡面沟与河道中的迁移转化过程均使用对流弥散方程进行描述，除对流扩散过程之外还包括旁侧入流、上一个子流域的河道到下一个子流域的河道以及不同组分之间的相互转化作用。
8.根据权利要求7所述的大流域尺度水氮迁移耦合模拟方法，其特征在于：
若坡面沟或河道窄而直，则Ex＝5.93·d·v，d表示河流深度；
若坡面沟或河道宽大，则Ex＝0.11·v·w2/d，w表示坡面沟或河道宽度。
9.根据权利要求7所述的大流域尺度水氮迁移耦合模拟方法，其特征在于：河道中有机氮、氨氮以及硝态氮的旁侧入流包括地下水与河道的交换量以及该子流域各等高带的坡面沟汇入河道的氮量。
CN202110517237.7A 2021-05-12 2021-05-12 一种大流域尺度水氮迁移耦合模拟方法 Active CN113139354B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110517237.7A CN113139354B (zh) 2021-05-12 2021-05-12 一种大流域尺度水氮迁移耦合模拟方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110517237.7A CN113139354B (zh) 2021-05-12 2021-05-12 一种大流域尺度水氮迁移耦合模拟方法

Publications (2)

Publication Number Publication Date
CN113139354A true CN113139354A (zh) 2021-07-20
CN113139354B CN113139354B (zh) 2021-12-07

Family

ID=76816927

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110517237.7A Active CN113139354B (zh) 2021-05-12 2021-05-12 一种大流域尺度水氮迁移耦合模拟方法

Country Status (1)

Country Link
CN (1) CN113139354B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114707371A (zh) * 2021-12-30 2022-07-05 水利部交通运输部国家能源局南京水利科学研究院 一种干湿交替带氮迁移转化耦合模拟方法及系统

Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20080142438A1 (en) * 2006-12-19 2008-06-19 Kent Greg B Wetland water treatment system
CN101908100A (zh) * 2010-07-26 2010-12-08 中国科学院生态环境研究中心 一种地下水环境的建模及数值模拟方法
CN102663267A (zh) * 2012-05-15 2012-09-12 南京大学 一种半湿润区流域面源污染负荷的确定方法
CN104143048A (zh) * 2014-07-23 2014-11-12 中国水利水电科学研究院 一种农田面源污染物入河量的计算方法及装置
CN104156570A (zh) * 2014-07-23 2014-11-19 中国水利水电科学研究院 一种土壤侵蚀类面源污染物入河量的计算方法及装置
CN105009768A (zh) * 2015-07-06 2015-11-04 中国农业科学院农业资源与农业区划研究所 一种流域尺度氮素化肥最大允许投入量的确定方法
CN106599471A (zh) * 2016-12-15 2017-04-26 中国水利水电科学研究院 一种分布式流域水沙过程耦合模拟方法
CN107590354A (zh) * 2017-10-26 2018-01-16 成都信息工程大学 基于稳定水流场的流域地表径流汇流模拟方法
CN110794111A (zh) * 2019-12-18 2020-02-14 济南大学 一种分析地下水中氯代烃污染迁移转化规律的方法

Patent Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20080142438A1 (en) * 2006-12-19 2008-06-19 Kent Greg B Wetland water treatment system
CN101908100A (zh) * 2010-07-26 2010-12-08 中国科学院生态环境研究中心 一种地下水环境的建模及数值模拟方法
CN102663267A (zh) * 2012-05-15 2012-09-12 南京大学 一种半湿润区流域面源污染负荷的确定方法
CN104143048A (zh) * 2014-07-23 2014-11-12 中国水利水电科学研究院 一种农田面源污染物入河量的计算方法及装置
CN104156570A (zh) * 2014-07-23 2014-11-19 中国水利水电科学研究院 一种土壤侵蚀类面源污染物入河量的计算方法及装置
CN105009768A (zh) * 2015-07-06 2015-11-04 中国农业科学院农业资源与农业区划研究所 一种流域尺度氮素化肥最大允许投入量的确定方法
CN106599471A (zh) * 2016-12-15 2017-04-26 中国水利水电科学研究院 一种分布式流域水沙过程耦合模拟方法
CN107590354A (zh) * 2017-10-26 2018-01-16 成都信息工程大学 基于稳定水流场的流域地表径流汇流模拟方法
CN110794111A (zh) * 2019-12-18 2020-02-14 济南大学 一种分析地下水中氯代烃污染迁移转化规律的方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
YANGWEN JIA 等: "Development of the WEP-L distributed hydrological model and dynamic assessment of water resources in the Yellow River Basin", 《JOURNAL OF HYDROLOGY》 *
刘扬: "基于流域水循环的微生物驱动氮素迁移转化机理研究——以青海湖流域为例", 《中国博士学位论文全文数据库 基础科学辑》 *
贺华翔 等: "基于二元水循环的流域分布式水质模型构建与应用", 《水利学报》 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114707371A (zh) * 2021-12-30 2022-07-05 水利部交通运输部国家能源局南京水利科学研究院 一种干湿交替带氮迁移转化耦合模拟方法及系统

Also Published As

Publication number Publication date
CN113139354B (zh) 2021-12-07

Similar Documents

Publication Publication Date Title
Jeon et al. Development of the HSPF-Paddy model to estimate watershed pollutant loads in paddy farming regions
Koch et al. Simulating unsteady flow, anabranching, and hyporheic dynamics in a glacial meltwater stream using a coupled surface water routing and groundwater flow model
Harlan et al. Rainfall-Runoff Modeling of citarum hulu river basin by using GR4j
Cebe et al. Monitoring and modeling land-based marine pollution
CN114021431A (zh) 一种基于地表-地下水耦合的tmdl模拟方法
CN113139354B (zh) 一种大流域尺度水氮迁移耦合模拟方法
Zhang et al. Field test of best management practice pollutant removal efficiencies in Shenzhen, China
CN116415456A (zh) 一种城市水循环模拟系统、方法、计算机设备及介质
Lodhi et al. Detention basins as best management practices for water quality control in an arid region
Noges Water and nutrient mass balance of the partly meromictic temperate Lake Verevi
De Baets Effects of hydroelectric dams on downstream oxygen and nitrogen in the Grand River and Conestogo River
Zorčič et al. Nitrate Concentration Changes In A River And Its Reservoir Within An Agriculturally-Influenced Watershed: The River Ledava (Se Austria And Ne Slovenia) Case Study
He et al. Unit process modeling of stormwater flow and pollutant sorption in a bioretention cell
Weintraub et al. WARMF: A watershed modeling tool for onsite wastewater systems
Koryto Evaluation of the Hydrologic and Water Quality Benefits of Regenerative Stormwater Conveyance.
Murphy Examining Nitrate Leaching Potential and Nitrogen Cycle Dynamics under Agricultural Managed Aquifer Recharge in the Central Valley of California
CN118536239A (zh) 基于一二维全水动力耦合模型的城市雨洪过程溯源方法
Wilby et al. Hydrological and hydrochemical dynamics
Eastman Field-scale nutrient transport monitoring and modeling of subsurface and naturally drained agricultural lands
Hellberg Modeling detention and pollutant fate in bioretention systems
Thompson Regenerative stormwater conveyances in Knoxville, TN: A performance analysis and investigation of design feature alternatives
Zhang et al. Shifts of the pond area ratio for lowland polders: Implication for nutrient control
Elci Experimental and theoretical study of subsurface flow and transport in a riparian wetland
Fonder Hydraulic and Removal Efficiencies of Horizontal Flow Treatemnt Wetlands
Kluge et al. Transport processes between Lake Belau and its drainage basin

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