CN113552305B - 一种多孔纤维材料埋设对产流水质影响的机理识别方法 - Google Patents

一种多孔纤维材料埋设对产流水质影响的机理识别方法 Download PDF

Info

Publication number
CN113552305B
CN113552305B CN202110827396.7A CN202110827396A CN113552305B CN 113552305 B CN113552305 B CN 113552305B CN 202110827396 A CN202110827396 A CN 202110827396A CN 113552305 B CN113552305 B CN 113552305B
Authority
CN
China
Prior art keywords
rainfall
runoff
experimental
nitrogen
phosphorus
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
Application number
CN202110827396.7A
Other languages
English (en)
Other versions
CN113552305A (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 CN202110827396.7A priority Critical patent/CN113552305B/zh
Publication of CN113552305A publication Critical patent/CN113552305A/zh
Application granted granted Critical
Publication of CN113552305B publication Critical patent/CN113552305B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N33/00—Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
    • G01N33/18—Water

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Chemical & Material Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Food Science & Technology (AREA)
  • Medicinal Chemistry (AREA)
  • Physics & Mathematics (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Cultivation Of Plants (AREA)

Abstract

本发明提供了一种多孔纤维材料埋设对产流水质影响的机理识别方法，属于农业环保技术领域，以多孔纤维材料体积和布局方式为控制因子设置正交试验；根据当地施肥经验确定肥料情况；确定作物生长期间的人工降雨次数、降雨强度和降雨历时；计算得到氮磷流失量；计算得到氮磷累计流失量，以及氮磷浓度的均值和变异系数，对氮磷浓度和累计流失量进行显著性检验，通过对照组和实验组之间的氮磷浓度和累计流失量相对差异值确定多孔纤维材料用量对于产流水质的影响，以及通过计算不同布局的实验组之间的氮磷浓度和累计流失量相对差异值确定多孔纤维材料布局对于氮磷浓度和累计流失量的影响。本发明为减小农田面源污染和河湖富营养化的目标做出探索。

Description

一种多孔纤维材料埋设对产流水质影响的机理识别方法
技术领域
本发明属于农业环保技术领域，尤其涉及一种多孔纤维材料埋设对产流水质影响的机理识别方法。
背景技术
农业面源污染已成为我国重要环保问题，源于化肥和农药滥用是导致河流湖泊富营养化等污染的主要原因。现有技术控制农田氮磷流失主要依靠减少化肥投入及增大利用效率等方面，但与农作物增产的目标存在相斥的方案使得控制农田面源污染效果和推广进程较差。
目前关于多孔纤维材料(岩棉)的应用集中在无土栽培领域，将其作物生长基质来促进作物对于营养液的吸收，但很少有人将岩棉的研究应用在改善农田地表产流水质上来，岩棉以自身所具有的高孔隙性、吸水性、稳定性等特性能够增强土壤最大持水能力，增大入渗并减小地表产流，而水体是氮磷运移的载体，减小地表径流就意味这减小氮磷流失风险，进而改善河流湖泊水质。然而截至到目前，多孔纤维材料改善产流水质的机理并未完全弄清。
发明内容
针对现有技术中的上述不足，本发明提供的一种多孔纤维材料埋设对产流水质影响的机理识别方法，通过埋设不同体积和不同布局多孔纤维材料，识别多孔纤维材料对于产流水质的影响机理，为减小农田面源污染和河湖富营养化的目标做出探索。
为了达到以上目的，本发明采用的技术方案为：
本方案提供一种多孔纤维材料埋设对产流水质影响的机理识别方法，包括以下步骤：
S1、将实验小区分为未埋设多孔纤维材料的空白对照组以及按梯度设置多孔纤维材料埋设体积的实验组，并确定各实验小区的规格、内容以及多孔纤维材料的布设方式；
S2、在作物种植过程中，根据当地种植经验在播种前进行施肥，同时在返青期进行追肥；
S3、基于当地多年平均降雨量及灌溉水量，结合作物不同生长期的需水规律，分别设定降雨时长、降雨次数、降雨强度和降雨总量进行降雨产流实验；
S4、设置第一时间间隔，并按第一时间间隔记录降雨产流实验开始后单位时间内的产流总量，并根据所述产流总量计算得到某一时刻的产流流量，绘制产流流量过程曲线；
S5、设置第二时间间隔，在降雨产流实验开始后按第二时间间隔取实验小区出口断面的产流水体，测量产流水体的氮磷浓度以及绘制氮磷浓度变化曲线，并通过产流流量过程曲线和氮磷浓度变化曲线对实验小区产流过程中各时段氮磷流失量进行累加，得到氮磷累计流失量；
S6、计算空白对照组和实验组产流水体中氮磷浓度的均值和变差系数，并根据所述均值和变差系数对产流水体中氮磷浓度及氮磷累计流失量进行检验，判断多孔纤维材料对于产流水质是否有影响；
S7、根据判断结果，通过空白对照组和实验组之间产流水体中氮磷浓度和氮磷累计流失量的相对差异值确定多孔纤维材料用量对于产流水质的影响，并通过实验组不同布局之间的氮磷浓度和氮磷产流累计流失量相对差异值确定多孔纤维材料布局对于氮磷浓度和氮磷产流累计流失量的影响，量化多孔纤维材料对产流过程中水质的影响程度大小，完成多孔纤维材料埋设对产流水质影响机理的识别。
本发明的有益效果是：本发明通过埋设不同体积和不同布局多孔纤维材料，识别多孔纤维材料对于产流水质的影响机理，为减小农田面源污染和河湖富营养化做出探索，通过对比分析在不同体积和布局的多孔纤维材料下土体单元产流水质的变化曲线，计算多孔纤维材料埋设下产流水体氮磷浓度和累计流失量的相对差异值，定量化识别多孔纤维材料埋设对产流水质的影响机理。本发明探寻多孔材料对于产流水质的影响机理，并给出详细的计算方法，为减小农田地表产流所导致的面源污染提供思路。
进一步地，所述步骤S1包括以下步骤：
S101、设所述实验小区为8个，其包括2个未埋设多孔纤维材料的空白对照组以及6个按梯度设置多孔纤维材料埋设体积的实验组，其中，所述6个实验组分为3组多孔纤维材料体积埋设方式，且每两组实验组之间多孔纤维材料埋设体积相同但埋设布局方式不同；
S102、在埋设等体积多孔纤维材料的两个实验小区之间设置A和B两种布局方式，将等体积的多孔纤维材料均分为m和n块，并将两个实验小区之间的埋块长边分别设置为东西方向和南北方向，即编号为2、3、4的实验小区设置为南北向，多孔材料埋设体积依次为V1、2V1、3V1，设为A布局，编号为6、7、8的实验小区设置为东西向，多孔材料埋设体积依次为V1、2V1、3V1，设为B布局，1号和5号实验小区为空白对照组；
S103、设多孔纤维材料埋块的中心位置高度与实验小区有效深度内的中心位置高度相等，并均匀埋设于实验小区的同等深度；
S104、在各实验小区边缘位置布置隔离板进行分隔，在实验小区边缘位置设置防风帷幕挡风，根据当地地形情况在实验小区内设置一定的倾斜角度，在实验小区坡脚处设置出口断面，并在出口断面外设置三角堰以测量流量及取水样。
上述进一步方案的有益效果是：本发明通过设计实验小区的各个组成部分，最大限度地减小环境因素的因素，并确定多孔材料的埋设和布局方案，为接下来的实验进行做铺垫。
再进一步地，所述步骤S3包括以下步骤：
S301、通过当地实验站获取历年降雨数据，计算得到作物生长期间的多年平均降雨量，并根据多年平均降雨量得到整个作物生长期间的降雨总量；
S302、根据所述降雨总量，结合多年平均月降雨分布规律以及作物不同生长期的需水规律设置降雨实验，将灌溉水量按天数均分，在需水强度低的作物生长期按多年平均降雨量与平均灌溉水量之和进行人工降雨补水，以供植物正常生长，在需水强度高的作物生长期进行降雨产流实验，以验证多孔纤维材料对于产流水质的影响，且整个作物生长期内的人工降雨量为多年平均降雨量与灌溉水量之和，得到不同生长期的单次降雨量以及降雨次数；
S303、根据当地土壤类型及获取的历年降雨数据，确定降雨强度以及降雨时长，并根据单次降雨量、灌溉水量、降雨强度以及降雨次数进行降雨产流实验，其中，所述单次降雨量为作物不同生长期的设计降雨量，其包括：需水强度低的作物生长期人工降雨量低，仅供正常生长；需水强度高的作物生长期人工降雨量高，进行降雨产流实验。
上述进一步方案的有益效果是：本发明根据设计降雨量，契合作物生长期内降雨规律，最大限度地模拟当地环境下的作物生长条件，进而确定在作物生长过程中，多孔材料埋设对产流水质的影响。
再进一步地，所述步骤S4包括以下步骤：
S401、设置第一时间间隔，针对降雨产流实验开始后，利用实验小区出口断面的三角堰测量单位时间内的产流总量；
S402、根据所述单位时间内的产流总量，计算得到不同实验小区、不同降雨场次和不同时间段的产流流量，并绘制产流流量过程曲线。
上述进一步方案的有益效果是：根据设定的实验，确定的时间间隔进行产流流量测定以及水样氮磷监测，可以看出多孔材料对于不同生长期、不同时间段的产流流量影响，判断多孔材料对于产流流量的影响，为后续氮磷相关计算提供支持。
再进一步地，所述步骤S402中产流流量的表达式如下：
Figure BDA0003174107950000051
其中，Qijk表示在第i个实验小区、第j场降雨和第k个时间段的产流流量，单位为cm3/s，Wijk表示在第i个实验小区、第j场降雨和第k个时间段的单位时间内产流总量，单位为cm3，tijk在第i个实验小区、第j场降雨和第k个时间段时取得Wijk体积水量所用的时间，单位为s，i表示实验小区编码，取值为1-8，其中1号和5号实验小区为空白实验组，多孔材料埋设用量为0，2-4实验小区为A布局，多孔纤维材料埋设用量依次为V1、2V1和3V1，6-8实验小区为B布局，多孔纤维材料埋设用量依次为0、V1、2V1和3V1，j表示第j次人工降雨场次，取值为1-G，k表示第k个时间段，取值为1-N，G表示人工降雨总次数，N表示总的时间段。
上述进一步方案的有益效果是：本发明通过计算得出i个实验小区、第j场降雨和第k个时间段的产流流量，判断不同埋设条件下的多孔材料对产流流量的影响，为后续氮磷浓度及累计流失量的计算做铺垫。
再进一步地，所述步骤S5包括以下步骤：
S501、设置第二时间间隔，针对降雨产流实验开始后，在测量产流流量的同时按第二时间间隔从实验小区出口断面的三角堰上取200ml的产流水样，并保持样品温度；
S502、测量产流水样的氮磷浓度，并绘制氮磷浓度曲线；
S503、根据所述产流流量曲线和氮磷浓度变化曲线，通过产流流量过程曲线和氮磷浓度变化曲线对实验小区产流过程中各时段氮磷流失量进行累加，得到对应实验小区、对应降雨场次和对应时间段内的氮磷累计流失量。
上述进一步方案的有益效果是：本发明通过以上设计，以确定多孔材料对于不同时间产流水体中氮磷浓度的影响，以及根据上述流量数据计算出累计流失量，确定多孔材料对于不同生长期的单次降雨0-Tk时间段内氮磷累计流失量的变化，确定不同实验小区氮磷流失量随时间的变化特征。
再进一步地，所述步骤S503中氮磷累计流失量的表达式如下：
Figure BDA0003174107950000061
Figure BDA0003174107950000062
其中，WNijk表示在第i个实验小区、第j场降雨和第T(k-1)至Tk时间段内的总氮累计流失量，单位为mg，T表示产流流量测定或水样采集的时间间隔，单位为s,k表示第k个时间段，取值为1-N,Tk表示从产流开始第0分钟至第Tk分钟期间的时长，单位为min，WPijk表示在第i个实验小区、第j场降雨和第T(k-1)至Tk时间段内的总磷累计流失量，单位为mg，Qijk表示在第i个实验小区、第j场降雨和第k个时间段的产流流量，单位为cm3/s，Qij(k-1)表示在第i个实验小区、第j场降雨和第k-1个时间段的产流流量，Pij(k-1)表示第i个实验小区、第j场降雨和第k-1个时间段的产流水样总磷浓度，单位为mg/L，Pijk表示在第i个实验小区、第j场降雨和第k个时间段的产流水样总磷浓度，单位为mg/L，Nij(k-1)表示在第i个实验小区、第j场降雨和第k-1个时间段的产流水样总氮浓度，Nijk表示在第i个实验小区、第j场降雨和第k个时间段的产流水样总氮浓度，单位为mg/L。
上述进一步方案的有益效果是：本发明通过以上设计计算得出第i个实验小区、第j场降雨在产流开始后第0秒至Tk秒时间段内的总氮累计流失量，通过绘图可以看出在多孔材料的埋设体积和布局差异下，降雨产流期间不同实验小区、不同生长期氮磷累计流失量随时间变化的特征。
再进一步地，所述步骤S6中氮磷浓度均值的表达式如下：
Figure BDA0003174107950000071
Figure BDA0003174107950000072
其中，
Figure BDA0003174107950000073
表示第i个实验小区、第j场降雨的总氮平均浓度，单位为mg/L，
Figure BDA0003174107950000074
表示第i个实验小区、第j场降雨的总磷平均浓度，单位为mg/L；
所述氮磷浓度变异系数的表达式如下：
Figure BDA0003174107950000075
Figure BDA0003174107950000076
其中，CV(Nij)表示第i个实验小区、第j场降雨下总氮平均浓度的变差系数，CV(Pij)表示第i个实验小区、第j场降雨下总磷平均浓度的变差系数。
上述进一步方案的有益效果是：通过氮磷浓度的均值判断多孔材料对于产流过程中水体氮磷浓度的影响，通过氮磷浓度的变异系数则是进一步说明在这个过程中水体氮磷浓度随时间变化过程中的离散程度。
再进一步地，所述步骤S7中多孔纤维材料用量对于产流水质的影响包括：不同多孔材料用量对产流水体中氮磷浓度的影响以及不同多孔材料用量对产流实验中氮磷累计流失量的影响；
所述不同多孔材料用量对产流水体中氮磷浓度的影响的计算表达式如下：
Figure BDA0003174107950000081
Figure BDA0003174107950000082
所述不同多孔材料用量对产流实验中氮磷累计流失量的影响的计算表达式如下：
Figure BDA0003174107950000083
Figure BDA0003174107950000084
其中，
Figure BDA0003174107950000085
表示在不同的多孔纤维材料用量对产流水体氮浓度的影响大小，即实验组和空白对照组产流水体总氮浓度的相对差异值，单位为mg/L，
Figure BDA0003174107950000086
表示在不同的多孔纤维材料用量对产流水体磷浓度的影响大小，即实验组和空白对照组产流水体总磷浓度的相对差异值，单位为mg/L，ΔWNi表示在不同的多孔纤维材料用量对产流水体氮累计流失量的影响大小，即实验组和空白对照组产流水体总氮的累积流失量的相对差异值，单位为mg，ΔWPi表示在不同的多孔纤维材料用量对产流水体磷累计流失量的影响大小，即实验组和空白对照组产流水体总磷的累积流失量的相对差异值，单位为mg，i表示实验小区编号，此处i取值为2、3、4，
Figure BDA0003174107950000091
表示第i号实验小区在第j场降雨中的总氮平均浓度，单位为mg/L，
Figure BDA0003174107950000092
表示第i+4号实验小区在第j场降雨中的总氮平均浓度，单位为mg/L，
Figure BDA0003174107950000093
表示第1号实验小区在第j场降雨中的产流水体总氮平均浓度，单位为mg/L，
Figure BDA0003174107950000094
表示第5号实验小区在第j场降雨中的产流水体总氮平均浓度，单位为mg/L，
Figure BDA0003174107950000095
表示第i号实验小区在第j场降雨中的产流水体总磷平均浓度，单位为mg/L，
Figure BDA0003174107950000096
表示第i+4号实验小区在第j场降雨中的产流水体总磷平均浓度，单位为mg/L，
Figure BDA0003174107950000097
表示第1号实验小区在第j场降雨中的产流水体总磷平均浓度，单位为mg/L，
Figure BDA0003174107950000098
表示第5号实验小区在第j场降雨中的产流水体总磷平均浓度，单位为mg/L，WNij表示第i号实验小区在第j场降雨中的产流水体总氮流失量，单位为mg，WN(i+4)j表示第i+4号实验小区在第j场降雨中的产流水体总氮流失量，单位为mg，WN1j表示第1号实验小区在第j场降雨中的产流水体总氮流失量，单位为mg，WN5j表示第5号实验小区在第j场降雨中的产流水体总氮流失量，单位为mg，WPij表示第i号实验小区在第j场降雨中的产流水体总磷流失量，单位为mg，WP(i+4)j表示第i+4号实验小区在第j场降雨中的产流水体总磷流失量，单位为mg，WP1j表示第1号实验小区在第j场降雨中的产流水体总磷流失量，单位为mg，WP5j表示第5号实验小区在第j场降雨中的产流水体总磷流失量，单位为mg。
上述进一步方案的有益效果是：通过实验组和对照组氮磷浓度的相对差异值来判断多孔材料对于产流水体中氮磷浓度的影响，而通过氮磷流失总量的相对差异值来判断产流水体中氮磷累计流失量的影响，进而确定多孔材料埋设用量对于产流水质的影响机理。
再进一步地，所述步骤S7中实验组不同布局之间的氮磷浓度和产流累计流失量相对差异值的表达式如下：
Figure BDA0003174107950000101
Figure BDA0003174107950000102
Figure BDA0003174107950000103
Figure BDA0003174107950000104
其中，ΔN表示在不同的多孔纤维材料布局对产流水体氮浓度的影响大小，即A布局和B布局之间的相对差异值，单位为mg/L，ΔP表示在不同的多孔纤维材料布局对产流水体磷浓度的影响大小，即A布局和B布局之间的相对差异值，mg/L，ΔWN表示在不同的多孔纤维材料布局对产流水体氮累计流失量的影响大小，即A布局和B布局之间的相对差异值，mg，ΔWP表示在不同的多孔纤维材料布局对产流水体磷累计流失量的影响大小，即A布局和B布局之间的相对差异值，单位为mg，i表示埋设多孔纤维材料的实验小区变化，取值为2、3、4。
上述进一步方案的有益效果是：通过实验组之间A和B布局下氮磷浓度的相对差异值来判断多孔材料对于产流水体中氮磷浓度的影响，而通过氮磷流失总量的相对差异值来判断产流水体中氮磷累计流失量的影响，进而确定多孔材料布局对于产流水质的影响机理，最终通过多孔材料埋设用量和布局确定其对于产流水质的影响机理。
附图说明
图1为本发明的方法流程图。
图2为本实例中设置的多孔材料埋设示意图。
具体实施方式
下面对本发明的具体实施方式进行描述，以便于本技术领域的技术人员理解本发明，但应该清楚，本发明不限于具体实施方式的范围，对本技术领域的普通技术人员来讲，只要各种变化在所附的权利要求限定和确定的本发明的精神和范围内，这些变化是显而易见的，一切利用本发明构思的发明创造均在保护之列。
实施例
如图1所示，本发明提供了一种多孔纤维材料埋设对产流水质影响的机理识别方法，其实现方法如下：
S1、将实验小区分为未埋设多孔纤维材料的空白对照组以及按梯度设置多孔纤维材料埋设体积的实验组，并确定各实验小区的规格、内容以及多孔纤维材料的布设方式，其实现方法如下：
S101、设所述实验小区为8个，其包括2个未埋设多孔纤维材料的空白对照组以及6个按梯度设置多孔纤维材料埋设体积的实验组，其中，所述6个实验组分为3组多孔纤维材料体积埋设方式，且每两组实验组之间多孔纤维材料埋设体积相同但埋设布局方式不同；
S102、在埋设等体积多孔纤维材料的两个实验小区之间设置A和B两种布局方式，将等体积的多孔纤维材料均分为m和n块，并将两个实验小区之间的埋块长边分别设置为东西方向和南北方向，即编号为2、3、4的实验小区设置为南北向，多孔材料埋设体积依次为V1、2V1、3V1，设为A布局，编号为6、7、8的实验小区设置为东西向，多孔材料埋设体积依次为V1、2V1、3V1，设为B布局，1号和5号实验小区为空白对照组，如图2所示，图2中，编号为2.3.4的实验小区设置为南北向，多孔纤维材料埋设用量为1.08、2.16、3.24m3，设为A布局，编号为6.7.8的实验小区设置为东西向，多孔纤维材料埋设用量为1.08、2.16、3.24m3，设为B布局。
本实施例中，如实验小区的尺寸为3.8*5.3m，有效深度为1m，设置多孔纤维材料埋设水平为0m3、1.08m3、2.16m3、3.24m3。
本实施例中，在实验地块具有代表性的地方利用环刀取出原状土壤样品，然后通过室内环刀法得到当地土壤和多孔材料的最大吸水率分别为A1和A2，测量实验小区内部土壤最大吸水率为A1，并测定多孔材料饱和吸水率为A2，以提高土壤最大持水能力5％、10％、15％为目标，结合实验小区的尺寸和有效深度，计算得出单位面积上每埋设V1体积的多孔纤维材料可以使土壤最大持水能力提高K，以提高土壤最大持水能力为目标，从而确定各实验小区材料埋设体积分别为0、V1、2V1、3V1。
本实施例中，如实验小区的尺寸为3.8*5.3m，有效深度为1m，以提高土壤最大理论持水能力5％、10％和15％为目标，设置多孔纤维材料埋设水平为0m3、1.08m3、2.16m3、3.24m3。
S103、设所述实验小区为8个，其包括2个未埋设多孔纤维材料的空白对照组以及6个按梯度设置多孔纤维材料埋设体积的实验组，其中，所述6个实验组分为3组多孔纤维材料体积埋设方式，且每两组实验组之间多孔纤维材料埋设体积相同但埋设布局方式不同；
本实施例中，实验小区为8个，包括2个空白对照组和6个实验组，空白实验组未埋设材料，另外6个实验组分为3组材料体积埋设水平将实验小区编号为1-8，1和5号实验小区为空白对照组，其余均为埋设不同体积和布局多孔纤维材料的实验组。
本实施例中，以多孔纤维材料体积和布局方式为控制因子设置正交试验，多孔纤维材料体积设置3个水平，布局方式设置2个水平，设置6个埋设材料的实验组，同时设置2个空白对照组，空白实验组未埋设材料，另外6个实验组分为3组材料体积埋设水平，每两组实验组之间材料体积埋设相同但具体布置方向和尺寸存在差异，然后将多孔纤维材料均匀地埋设于实验小区中。
S104、在埋设等体积多孔纤维材料的两个实验小区之间设置A和B两种布局方式，将等体积的多孔纤维材料均分为m和n块，并将两个实验小区之间的埋块长边分别设置为东西方向和南北方向，即编号为2.3.4的实验小区设置为东西向，设为A布局，编号为6.7.8的实验小区设置为南北向，设为B布局；
本实施例中，在6个埋设有多孔纤维材料的实验组内，埋设等体积材料的两个实验小区之间，设置两种布局方式，每两组实验组之间材料体积埋设相同但具体布置方向存在差异，等体积且不同布局方式的两个实验小区之间的埋块长边分别设置为东西向和南北向。将等体积的多孔纤维材料均分为m和n块，即各实验小区的埋块尺寸按式计算：
Vi＝mi×ai×bi×h
Vq＝mq×aq×bq×h
其中，i表示2-4号实验小区编号，取值为2.3.4；表示6-8号实验小区编号，取值为6、7、8；Vi和Vq表示i、q实验小区埋块总体积，取值为V1、2V1、3V1，m3；mi和mq表示i、q实验小区埋块块数，块；ai和aq表示i、q实验小区长边长度，m；bi和bq表示i、q实验小区短边长度，m，h表示埋块高度，取值为0.5m。
S105、设多孔纤维材料埋块的中心位置高度与实验小区有效深度内的中心位置高度相等，并均匀埋设于实验小区的同等深度；
S106、在各实验小区边缘位置布置隔离板进行分隔，在实验小区边缘位置设置防风帷幕挡风，根据当地地形情况在实验小区内设置一定的倾斜角度，在实验小区坡脚处设置出口断面，并在出口断面外设置三角堰以测量流量及取水样。
本实施例中，如实验小区尺寸为3.8*5.3m，经过计算得出材料埋设体积分别为0、1.08m3、2.16m3、3.24m3三种水平，2、3、4、6、7、8号实验小区的埋块数*长*宽的分别为8*0.45*0.75m、10*0.45*1.2m、6*0.45*3m、6*0.45*1m、8*0.45*1.5m、4*0.45*4.5m，高均为40cm,埋设在土壤30-70cm的空间上。在6个埋设有多孔纤维材料的实验组内，埋设等体积材料的两个实验小区之间，设置两种布局方式。本实验中，多孔材料埋设的具体布置见附图2，将(2.3.4)和(6.7.8)的长边分别设置为南北向和东西向，均匀地埋设于实验小区内。在所有小区边缘位置设置2mm厚和1m深的隔离钢板，实验小区坡度设置为3°，同时在小区边缘设置3.8m宽*3m高的聚乙烯塑料布充当防风帷幕，在实验小区坡脚处设置出口断面，并在出口断面设置三角堰，以便测量流量及取水样。
S2、在作物种植过程中，根据当地种植经验在播种前进行施肥，同时在返青期进行追肥；
本实施例中，在作物生长过程中，依据当地种植施肥经验，在播种前按每亩地15kg施加尿素和每亩地50kg施加复合肥，并在返青期进行追肥，追肥标准为每亩地施加5kg尿素，经过面积换算后将对应质量的化肥均匀的撒在每个实验小区内部，同时每次施肥过后马上进行灌溉或降雨实验，使化肥充分渗入土壤。
S3、基于当地多年平均降雨量及当地农业部门建议的灌溉水量，结合作物不同生长期的需水规律，分别设定降雨时长、降雨次数(根据不同生长期需水规律设计)、降雨强度和降雨总量进行降雨产流实验，其实现方法如下：
S301、通过当地实验站获取历年降雨数据，计算得到作物生长期间的多年平均降雨量，并根据多年平均降雨量得到整个作物生长期间的降雨总量；
本实施例中，通过当地实验站获取历年降雨数据，计算得出作物生长期间(10月-次年6月)的多年平均降雨量，然后查询当地省份农业部门发布的作物建议灌溉水量，得出在整个作物生长期间总降雨量。
S302、根据所述降雨总量，结合多年平均月降雨分布规律以及作物不同生长期的需水规律设置降雨实验，将灌溉水量按天数均分，在需水强度低的作物生长期按多年平均降雨量与平均灌溉水量之和进行人工降雨补水，以供植物正常生长，在需水强度高的作物生长期进行降雨产流实验，以验证多孔纤维材料对于产流水质的影响，且整个作物生长期内的人工降雨量为多年平均降雨量与灌溉水量之和，得到不同生长期的单次降雨量以及降雨次数；
本实施例中，结合多年平均月降雨分布规律以及作物不同生长期需水规律，拔节期之前按多年平均降雨量灌溉保证作物正常生长，然后在拔节期及之后设计5次典型降雨实验，计算得出典型降雨实验单次降雨量为W。
S303、根据当地土壤类型及获取的历年降雨数据，确定降雨强度以及降雨时长，并根据单次降雨量、灌溉水量、降雨强度以及降雨次数进行降雨产流实验，其中，所述单次降雨量为作物不同生长期的设计降雨量，其包括：需水强度低的作物生长期人工降雨量低，仅供正常生长；需水强度高的作物生长期人工降雨量高，进行降雨产流实验。
本实施例中，根据作物各生长期需水规律来设置典型降雨实验，越冬期之前按当地年平均降雨量保证作物正常生长，然后分别在作物不同生长期和裸土下进行人工降雨实验。结合当地土壤入渗特性及地表产流要求，每个实验小区内进行降雨实验时的降雨强度均设置为P，降雨强度参照当地暴雨级别，降雨总量均设置为W，分别在作物的不同生育期下进行降雨产流实验，整个作物生长期间人工降雨量等于当地多年平均降雨量与灌溉水量之和。
本实施例中，结合淮北平原五道沟实验站实验数据，自1986-2014年间(当年10月-次年6月)平均降雨量为466.1mm，加之安徽省灌溉量为50mm，将人工降雨量设计为516.1mm，同时根据实验站所测得的冬小麦不同生长期需水规律，在越冬期前仅按照多年平均降雨量进行正常灌溉，在冬小麦的拔节期、孕穗期、扬花期和灌浆期进行降雨产流实验，每次降雨量为100mm，同时设置降雨强度为50mm/h,对应当地降雨强度等级的暴雨级别。同时在冬小麦收割后进行同等条件下的降雨实验进行对比。
S4、设置第一时间间隔，并按第一时间间隔记录降雨产流实验开始后单位时间内的产流总量，并根据所述产流总量计算得到某一时刻的产流流量，绘制产流流量过程曲线，其实现方法如下：
S401、设置第一时间间隔，针对降雨产流实验开始后，利用实验小区出口断面的三角堰测量单位时间内的产流总量；
S402、根据所述单位时间内的产流总量，计算得到不同实验小区、不同降雨场次和不同时间段的产流流量，并绘制产流流量过程曲线。
本实施例中，进行降雨实验时，自产流开始后，设定固定间隔时间T，利用实验小区出口断面的三角堰测量单位时间内的产流水量，然后计算不同实验小区、不同降雨场次和不同时间段的产流流量，并绘制流量过程曲线。产流流量计算公式如下：
Figure BDA0003174107950000161
其中，Qijk表示在第i个实验小区、第j场降雨和第k个时间段的产流流量，单位为cm3/s，Wijk表示在第i个实验小区、第j场降雨和第k个时间段的单位时间内产流总量，单位为cm3，tijk在第i个实验小区、第j场降雨和第k个时间段时取得Wijk体积水量所用的时间，单位为s，i表示实验小区编码，取值为1-8，其中1-4实验小区为A布局，多孔纤维材料埋设用量依次为0、V1、2V1和3V1，5-8实验小区为B布局，多孔纤维材料埋设用量依次为0、V1、2V1和3V1，j表示第j次人工降雨，取值为1-G，代表该次降雨实验对应的场数，k表示第k个时间段，取值为1-N，G表示人工降雨总次数，N表示总的时间段。
本实施例中，每次进行降雨实验时，需保证在无风或微风条件下进行，同时使用并固定防风帷幕进行挡风。自产流实验开始后，设定固定间隔时间5分钟，利用实验小区出口断面的三角堰使用容器采集水体，并记录各实验小区和降雨场次的单位时间内采集到的产流水体质量，所有的产流水体密度均取1g/ml，然后将产流水体的质量换算成产流水体体积Wijk，然后计算不同实验小区、不同降雨场次和不同时间段的产流流量，并绘制流量过程曲线。
S5、设置第二时间间隔，在降雨产流实验开始后按第二时间间隔取实验小区出口断面的产流水体，测量产流水体的氮磷浓度以及绘制氮磷浓度变化曲线，并通过产流流量过程曲线和氮磷浓度变化曲线对实验小区产流过程中各时段氮磷流失量进行累加，得到氮磷累计流失量，其实现方法如下：
S501、设置第二时间间隔，针对降雨产流实验开始后，在测量产流流量的同时按第二时间间隔从实验小区出口断面的三角堰上取200ml的产流水样，并保持样品温度；
S502、测量产流水样的氮磷浓度，并绘制氮磷浓度曲线；
本实施例中，降雨实验开始后，设定时间间隔为T，通过实验小区出口断面的三角堰取200ml的水样，保温然后送至实验室，分别通过碱性过硫酸钾消解紫外分光光度法和钼酸铵分光光度法分别测量水体氮磷浓度，并绘制氮磷浓度变化曲线。
S503、根据所述产流流量曲线和氮磷浓度变化曲线，通过产流流量过程曲线和氮磷浓度变化曲线对实验小区产流过程中各时段氮磷流失量进行累加，得到对应实验小区、对应降雨场次和对应时间段内的氮磷累计流失量：氮磷累计流失量计算公式如下：
Figure BDA0003174107950000181
Figure BDA0003174107950000182
其中，WNijk表示在第i个实验小区、第j场降雨和第T(k-1)至Tk时间段内的总氮累计流失量，单位为mg，T表示产流流量测定或水样采集的时间间隔，单位为s,k表示第k个时间段，取值为1-N,Tk表示从产流开始第0分钟至第Tk分钟期间的时长，单位为min，WPijk表示在第i个实验小区、第j场降雨和第T(k-1)至Tk时间段内的总磷累计流失量，单位为mg，Qijk表示在第i个实验小区、第j场降雨和第k个时间段的产流流量，单位为cm3/s，Qij(k-1)表示在第i个实验小区、第j场降雨和第k-1个时间段的产流流量，Pij(k-1)表示第i个实验小区、第j场降雨和第k-1个时间段的产流水样总磷浓度，单位为mg/L，Pijk表示在第i个实验小区、第j场降雨和第k个时间段的产流水样总磷浓度，单位为mg/L，Nij(k-1)表示在第i个实验小区、第j场降雨和第k-1个时间段的产流水样总氮浓度，Nijk表示在第i个实验小区、第j场降雨和第k个时间段的产流水样总氮浓度，单位为mg/L。
本实施例中，设定第二时间间隔为T，自产流实验开始后，在测定产流流量的同时从实验小区出口断面的三角堰上取出200ml的产流水样，保持样品温度，然后送到实验室测量水体氮磷浓度，并绘制氮磷浓度变化曲线。根据产流流量曲线和相应氮磷浓度变化曲线，计算得出对应实验小区、对应场次和对应时间段的氮磷累计流失量。
S6、计算空白对照组和实验组产流水体中氮磷浓度的均值和变差系数，并根据所述均值和变差系数对产流水体中氮磷浓度及氮磷累计流失量进行检验，判断多孔纤维材料对于产流水质是否有影响；
本实施例中，根据产流流量变化曲线和氮磷浓度变化曲线，对实验小区产流过程中的各时段氮磷流失量进行累加，绘制出氮磷累计流失量变化曲线，计算得出对照组和实验组产流水体氮磷浓度的均值和变差系数，然后使用F检验对氮磷浓度和氮磷累计流失量进行显著性检验，根据对照组和实验组之间氮磷浓度变化的均值和变差系数，以及氮磷浓度变化和累计流失量的显著性检验结果判断多孔纤维材料用量和布局对于产流水质是否有影响。
Figure BDA0003174107950000191
Figure BDA0003174107950000192
其中，
Figure BDA0003174107950000193
表示第i个实验小区、第j场降雨的总氮平均浓度，单位为mg/L，
Figure BDA0003174107950000194
表示第i个实验小区、第j场降雨的总磷平均浓度，单位为mg/L；
所述氮磷浓度变异系数的表达式如下：
Figure BDA0003174107950000195
Figure BDA0003174107950000196
其中，CV(Nij)表示第i个实验小区、第j场降雨下总氮平均浓度的变差系数，CV(Pij)表示第i个实验小区、第j场降雨下总磷平均浓度的变差系数。
S7、根据判断结果，通过空白对照组和实验组之间产流水体中氮磷浓度和氮磷累计流失量的相对差异值确定多孔纤维材料用量对于产流水质的影响，并通过实验组不同布局之间的氮磷浓度和产流累计流失量相对差异值确定多孔纤维材料布局对于氮磷浓度和产流累计流失量的影响，量化多孔纤维材料对产流过程中水质的影响程度大小，完成多孔纤维材料埋设对产流水质影响机理的识别。
本实施例中，通过对照组和实验组之间的氮磷浓度和累计流失量相对差异值确定多孔纤维材料用量对于产流水质的影响程度大小，通过(2-4)和(6-8)实验组不同布局之间的氮磷浓度和累计流失量相对差异值确定多孔纤维材料布局对于氮磷浓度和累计流失量的影响，进而明确多孔纤维材料布设对产流水质的影响。
本实施例中，所述步骤S7中空白对照组和实验组之间产流水体中氮磷浓度和氮磷累计流失量的相对差异的表达式如下：
Figure BDA0003174107950000201
Figure BDA0003174107950000202
Figure BDA0003174107950000203
Figure BDA0003174107950000204
其中，
Figure BDA0003174107950000211
表示在不同的多孔纤维材料用量对产流水体氮浓度的影响大小，单位为mg/L，
Figure BDA0003174107950000212
表示在不同的多孔纤维材料用量对产流水体磷浓度的影响大小，单位为mg/L，ΔWNi表示在不同的多孔纤维材料用量对产流水体氮累计流失量的影响大小，单位为mg，ΔWPi表示在不同的多孔纤维材料用量对产流水体磷累计流失量的影响大小，i表示实验小区编号，此处i取值为2、3、4，
Figure BDA0003174107950000213
表示第i号实验小区在第j场降雨中的总氮平均浓度，单位为mg/L，
Figure BDA0003174107950000214
表示第i+4号实验小区在第j场降雨中的总氮平均浓度，单位为mg/L，
Figure BDA0003174107950000215
表示第1号实验小区在第j场降雨中的产流水体总氮平均浓度，单位为mg/L，
Figure BDA0003174107950000216
表示第5号实验小区在第j场降雨中的产流水体总氮平均浓度，单位为mg/L，
Figure BDA0003174107950000217
表示第i号实验小区在第j场降雨中的产流水体总磷平均浓度，单位为mg/L，
Figure BDA0003174107950000218
表示示第i+4号实验小区在第j场降雨中的产流水体总磷平均浓度，单位为mg/L，
Figure BDA0003174107950000219
表示第1号实验小区在第j场降雨中的产流水体总磷平均浓度，单位为mg/L，
Figure BDA00031741079500002110
表示第5号实验小区在第j场降雨中的产流水体总磷平均浓度，单位为mg/L，WNij表示第i号实验小区在第j场降雨中的产流水体总氮流失量，单位为mg，WN(i+4)j表示第i+4号实验小区在第j场降雨中的产流水体总氮流失量，单位为mg，WN1j表示第1号实验小区在第j场降雨中的产流水体总氮流失量，单位为mg，WN5j表示第5号实验小区在第j场降雨中的产流水体总氮流失量，单位为mg，WPij表示第i号实验小区在第j场降雨中的产流水体总磷流失量，单位为mg，WP(i+4)j表示第i+4号实验小区在第j场降雨中的产流水体总磷流失量，单位为mg，WP1j表示第1号实验小区在第j场降雨中的产流水体总磷流失量，单位为mg，WP5j表示第5号实验小区在第j场降雨中的产流水体总磷流失量，单位为mg。
本实施例中，所述步骤S7中实验组不同布局之间的氮磷浓度和产流累计流失量相对差异值的表达式如下：
Figure BDA0003174107950000221
Figure BDA0003174107950000222
Figure BDA0003174107950000223
Figure BDA0003174107950000224
其中，ΔN表示在不同的多孔纤维材料布局对产流水体氮浓度的影响大小，即A布局(2、3、4号实验小区)和B布局(6、7、8号实验小区)之间的相对差异值，单位为mg/L，ΔP表示在不同的多孔纤维材料布局对产流水体磷浓度的影响大小，即A布局(2、3、4号实验小区)和B布局(6、7、8号实验小区)之间的相对差异值，mg/L，ΔWN表示在不同的多孔纤维材料布局对产流水体氮累计流失量的影响大小，即A布局(2、3、4号实验小区)和B布局(6、7、8号实验小区)之间的相对差异值，mg，ΔWP表示在不同的多孔纤维材料布局对产流水体磷累计流失量的影响大小，即A布局(2、3、4号实验小区)和B布局(6、7、8号实验小区)之间的相对差异值，单位为mg，i表示埋设多孔纤维材料的实验小区变化，取值为2、3、4。

Claims (9)

1.一种多孔纤维材料埋设对产流水质影响的机理识别方法，其特征在于，包括以下步骤：
S1、将实验小区分为未埋设多孔纤维材料的空白对照组以及按梯度设置多孔纤维材料埋设体积的实验组，并确定各实验小区的规格、内容以及多孔纤维材料的布设方式；
S2、在作物种植过程中，根据当地种植经验在播种前进行施肥，同时在返青期进行追肥；
S3、基于当地多年平均降雨量及灌溉水量，结合作物不同生长期的需水规律，分别设定降雨时长、降雨次数、降雨强度和降雨总量进行降雨产流实验；
S4、设置第一时间间隔，并按第一时间间隔记录降雨产流实验开始后单位时间内的产流总量，并根据所述产流总量计算得到某一时刻的产流流量，绘制产流流量过程曲线；
S5、设置第二时间间隔，在降雨产流实验开始后按第二时间间隔取实验小区出口断面的产流水体，测量产流水体的氮磷浓度以及绘制氮磷浓度变化曲线，并通过产流流量过程曲线和氮磷浓度变化曲线对实验小区产流过程中各时段氮磷流失量进行累加，得到氮磷累计流失量；
所述步骤S5包括以下步骤：
S501、设置第二时间间隔，针对降雨产流实验开始后，在测量产流流量的同时按第二时间间隔从实验小区出口断面的三角堰上取200ml的产流水样，并保持样品温度；
S502、测量产流水样的氮磷浓度，并绘制氮磷浓度曲线；
S503、根据所述产流流量曲线和氮磷浓度变化曲线，通过产流流量过程曲线和氮磷浓度变化曲线对实验小区产流过程中各时段氮磷流失量进行累加，得到对应实验小区、对应降雨场次和对应时间段内的氮磷累计流失量；
S6、计算空白对照组和实验组产流水体中氮磷浓度的均值和变差系数，并根据所述均值和变差系数对产流水体中氮磷浓度及氮磷累计流失量进行检验，判断多孔纤维材料对于产流水质是否有影响；
S7、根据判断结果，通过空白对照组和实验组之间产流水体中氮磷浓度和氮磷累计流失量的相对差异值确定多孔纤维材料用量对于产流水质的影响，并通过实验组不同布局之间的氮磷浓度和氮磷产流累计流失量相对差异值确定多孔纤维材料布局对于氮磷浓度和氮磷产流累计流失量的影响，量化多孔纤维材料对产流过程中水质的影响程度大小，完成多孔纤维材料埋设对产流水质影响机理的识别。
2.根据权利要求1所述的多孔纤维材料埋设对产流水质影响的机理识别方法，其特征在于，所述步骤S1包括以下步骤：
S101、设所述实验小区为8个，其包括2个未埋设多孔纤维材料的空白对照组以及6个按梯度设置多孔纤维材料埋设体积的实验组，其中，所述6个实验组分为3组多孔纤维材料体积埋设方式，且每两组实验组之间多孔纤维材料埋设体积相同但埋设布局方式不同；
S102、在埋设等体积多孔纤维材料的两个实验小区之间设置A和B两种布局方式，将等体积的多孔纤维材料均分为m和n块，并将两个实验小区之间的埋块长边分别设置为东西方向和南北方向，即编号为2、3、4的实验小区设置为南北向，多孔材料埋设体积依次为V1、2V1、3V1，设为A布局，编号为6、7、8的实验小区设置为东西向，多孔材料埋设体积依次为V1、2V1、3V1，设为B布局，1号和5号实验小区为空白对照组；
S103、设多孔纤维材料埋块的中心位置高度与实验小区有效深度内的中心位置高度相等，并均匀埋设于实验小区的同等深度；
S104、在各实验小区边缘位置布置隔离板进行分隔，在实验小区边缘位置设置防风帷幕挡风，根据当地地形情况在实验小区内设置一定的倾斜角度，在实验小区坡脚处设置出口断面，并在出口断面外设置三角堰以测量流量及取水样。
3.根据权利要求2所述的多孔纤维材料埋设对产流水质影响的机理识别方法，其特征在于，所述步骤S3包括以下步骤：
S301、通过当地实验站获取历年降雨数据，计算得到作物生长期间的多年平均降雨量，并根据多年平均降雨量得到整个作物生长期间的降雨总量；
S302、根据所述降雨总量，结合多年平均月降雨分布规律以及作物不同生长期的需水规律设置降雨实验，将灌溉水量按天数均分，在需水强度低的作物生长期按多年平均降雨量与平均灌溉水量之和进行人工降雨补水，以供植物正常生长，在需水强度高的作物生长期进行降雨产流实验，以验证多孔纤维材料对于产流水质的影响，且整个作物生长期内的人工降雨量为多年平均降雨量与灌溉水量之和，得到不同生长期的单次降雨量以及降雨次数；
S303、根据当地土壤类型及获取的历年降雨数据，确定降雨强度以及降雨时长，并根据单次降雨量、灌溉水量、降雨强度以及降雨次数进行降雨产流实验，其中，所述单次降雨量为作物不同生长期的设计降雨量，其包括：需水强度低的作物生长期人工降雨量低，仅供正常生长；需水强度高的作物生长期人工降雨量高，进行降雨产流实验。
4.根据权利要求3所述的多孔纤维材料埋设对产流水质影响的机理识别方法，其特征在于，所述步骤S4包括以下步骤：
S401、设置第一时间间隔，针对降雨产流实验开始后，利用实验小区出口断面的三角堰测量单位时间内的产流总量；
S402、根据所述单位时间内的产流总量，计算得到不同实验小区、不同降雨场次和不同时间段的产流流量，并绘制产流流量过程曲线。
5.根据权利要求4所述的多孔纤维材料埋设对产流水质影响的机理识别方法，其特征在于，所述步骤S402中产流流量的表达式如下：
Figure FDA0004024783240000041
其中，Qijk表示在第i个实验小区、第j场降雨和第k个时间段的产流流量，单位为cm3/s，Wijk表示在第i个实验小区、第j场降雨和第k个时间段的单位时间内产流总量，单位为cm3，tijk在第i个实验小区、第j场降雨和第k个时间段时取得Wijk体积水量所用的时间，单位为s，i表示实验小区编码，取值为1-8，其中1号和5号实验小区为空白实验组，多孔材料埋设用量为0，2-4实验小区为A布局，多孔纤维材料埋设用量依次为V1、2V1和3V1，6-8实验小区为B布局，多孔纤维材料埋设用量依次为0、V1、2V1和3V1，j表示第j次人工降雨场次，取值为1-G，k表示第k个时间段，取值为1-N，G表示人工降雨总次数，N表示总的时间段。
6.根据权利要求5所述的多孔纤维材料埋设对产流水质影响的机理识别方法，其特征在于，所述步骤S503中氮磷累计流失量的表达式如下：
Figure FDA0004024783240000042
Figure FDA0004024783240000043
其中，WNijk表示在第i个实验小区、第j场降雨和第T(k-1)至Tk时间段内的总氮累计流失量，单位为mg，T表示产流流量测定或水样采集的时间间隔，单位为s,k表示第k个时间段，取值为1-N,Tk表示从产流开始第0分钟至第Tk分钟期间的时长，单位为min，WPijk表示在第i个实验小区、第j场降雨和第T(k-1)至Tk时间段内的总磷累计流失量，单位为mg，Qijk表示在第i个实验小区、第j场降雨和第k个时间段的产流流量，单位为cm3/s，Qij(k-1)表示在第i个实验小区、第j场降雨和第k-1个时间段的产流流量，Pij(k-1)表示第i个实验小区、第j场降雨和第k-1个时间段的产流水样总磷浓度，单位为mg/L，Pijk表示在第i个实验小区、第j场降雨和第k个时间段的产流水样总磷浓度，单位为mg/L，Nij(k-1)表示在第i个实验小区、第j场降雨和第k-1个时间段的产流水样总氮浓度，Nijk表示在第i个实验小区、第j场降雨和第k个时间段的产流水样总氮浓度，单位为mg/L。
7.根据权利要求6所述的多孔纤维材料埋设对产流水质影响的机理识别方法，其特征在于，所述步骤S6中氮磷浓度均值的表达式如下：
Figure FDA0004024783240000051
Figure FDA0004024783240000052
其中，
Figure FDA0004024783240000053
表示第i个实验小区、第j场降雨的总氮平均浓度，单位为mg/L，
Figure FDA0004024783240000054
表示第i个实验小区、第j场降雨的总磷平均浓度，单位为mg/L；
所述氮磷浓度变差系数的表达式如下：
Figure FDA0004024783240000055
Figure FDA0004024783240000056
其中，CV(Nij)表示第i个实验小区、第j场降雨下总氮平均浓度的变差系数，CV(Pij)表示第i个实验小区、第j场降雨下总磷平均浓度的变差系数。
8.根据权利要求7所述的多孔纤维材料埋设对产流水质影响的机理识别方法，其特征在于，所述步骤S7中多孔纤维材料用量对于产流水质的影响包括：不同多孔材料用量对产流水体中氮磷浓度的影响以及不同多孔材料用量对产流实验中氮磷累计流失量的影响；
所述不同多孔材料用量对产流水体中氮磷浓度的影响的计算表达式如下：
Figure FDA0004024783240000061
Figure FDA0004024783240000062
所述不同多孔材料用量对产流实验中氮磷累计流失量的影响的计算表达式如下：
Figure FDA0004024783240000063
Figure FDA0004024783240000064
其中，
Figure FDA0004024783240000065
表示在不同的多孔纤维材料用量对产流水体氮浓度的影响大小，即实验组和空白对照组产流水体总氮浓度的相对差异值，单位为mg/L，
Figure FDA0004024783240000066
表示在不同的多孔纤维材料用量对产流水体磷浓度的影响大小，即实验组和空白对照组产流水体总磷浓度的相对差异值，单位为mg/L，ΔWNi表示在不同的多孔纤维材料用量对产流水体氮累计流失量的影响大小，即实验组和空白对照组产流水体总氮的累积流失量的相对差异值，单位为mg，ΔWPi表示在不同的多孔纤维材料用量对产流水体磷累计流失量的影响大小，即实验组和空白对照组产流水体总磷的累积流失量的相对差异值，单位为mg，i表示实验小区编号，此处i取值为2、3、4，
Figure FDA0004024783240000071
表示第i号实验小区在第j场降雨中的总氮平均浓度，单位为mg/L，
Figure FDA0004024783240000072
表示第i+4号实验小区在第j场降雨中的总氮平均浓度，单位为mg/L，
Figure FDA0004024783240000073
表示第1号实验小区在第j场降雨中的产流水体总氮平均浓度，单位为mg/L，
Figure FDA0004024783240000074
表示第5号实验小区在第j场降雨中的产流水体总氮平均浓度，单位为mg/L，
Figure FDA0004024783240000075
表示第i号实验小区在第j场降雨中的产流水体总磷平均浓度，单位为mg/L，
Figure FDA0004024783240000076
表示第i+4号实验小区在第j场降雨中的产流水体总磷平均浓度，单位为mg/L，
Figure FDA0004024783240000077
表示第1号实验小区在第j场降雨中的产流水体总磷平均浓度，单位为mg/L，
Figure FDA0004024783240000078
表示第5号实验小区在第j场降雨中的产流水体总磷平均浓度，单位为mg/L，WNij表示第i号实验小区在第j场降雨中的产流水体总氮流失量，单位为mg，WN(i+4)j表示第i+4号实验小区在第j场降雨中的产流水体总氮流失量，单位为mg，WN1j表示第1号实验小区在第j场降雨中的产流水体总氮流失量，单位为mg，WN5j表示第5号实验小区在第j场降雨中的产流水体总氮流失量，单位为mg，WPij表示第i号实验小区在第j场降雨中的产流水体总磷流失量，单位为mg，WP(i+4)j表示第i+4号实验小区在第j场降雨中的产流水体总磷流失量，单位为mg，WP1j表示第1号实验小区在第j场降雨中的产流水体总磷流失量，单位为mg，WP5j表示第5号实验小区在第j场降雨中的产流水体总磷流失量，单位为mg。
9.根据权利要求8所述的多孔纤维材料埋设对产流水质影响的机理识别方法，其特征在于，所述步骤S7中实验组不同布局之间的氮磷浓度和产流累计流失量相对差异值的表达式如下：
Figure FDA0004024783240000079
Figure FDA0004024783240000081
Figure FDA0004024783240000082
Figure FDA0004024783240000083
其中，ΔN表示在不同的多孔纤维材料布局对产流水体氮浓度的影响大小，即A布局和B布局之间的相对差异值，单位为mg/L，ΔP表示在不同的多孔纤维材料布局对产流水体磷浓度的影响大小，即A布局和B布局之间的相对差异值，mg/L，ΔWN表示在不同的多孔纤维材料布局对产流水体氮累计流失量的影响大小，即A布局和B布局之间的相对差异值，mg，ΔWP表示在不同的多孔纤维材料布局对产流水体磷累计流失量的影响大小，即A布局和B布局之间的相对差异值，单位为mg，i表示埋设多孔纤维材料的实验小区变化，取值为2、3、4。
CN202110827396.7A 2021-07-21 2021-07-21 一种多孔纤维材料埋设对产流水质影响的机理识别方法 Active CN113552305B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110827396.7A CN113552305B (zh) 2021-07-21 2021-07-21 一种多孔纤维材料埋设对产流水质影响的机理识别方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110827396.7A CN113552305B (zh) 2021-07-21 2021-07-21 一种多孔纤维材料埋设对产流水质影响的机理识别方法

Publications (2)

Publication Number Publication Date
CN113552305A CN113552305A (zh) 2021-10-26
CN113552305B true CN113552305B (zh) 2023-02-24

Family

ID=78103950

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110827396.7A Active CN113552305B (zh) 2021-07-21 2021-07-21 一种多孔纤维材料埋设对产流水质影响的机理识别方法

Country Status (1)

Country Link
CN (1) CN113552305B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114996911A (zh) * 2022-05-09 2022-09-02 中国水利水电科学研究院 一种农田产流的识别分析方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6332287B1 (en) * 1997-03-06 2001-12-25 Carroll M. Geraldson Plant cultivation apparatus and method
CA2843086A1 (en) * 2013-02-19 2014-08-19 Furbish Company, Llc Vegetative roofing systems
CN108918827A (zh) * 2018-07-02 2018-11-30 中国水利水电科学研究院 一种识别多孔材料对土壤持水性影响的方法
CN108956948A (zh) * 2018-07-02 2018-12-07 中国水利水电科学研究院 一种多孔材料对平原区产流影响的识别方法
CN108982807A (zh) * 2018-07-02 2018-12-11 中国水利水电科学研究院 一种识别多孔材料对土壤持水性影响的动态监测系统

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP5522626B2 (ja) * 2009-09-08 2014-06-18 国立大学法人島根大学 土壌改質方法
EP3777920A1 (en) * 2015-11-30 2021-02-17 Toray Industries, Inc. Phosphorus adsorbent, porous fiber and phosphorus adsorption column
CN108717453B (zh) * 2018-05-21 2020-05-05 河海大学 一种基于gis平台的平原河网区污染负荷计算方法
CN109190827B (zh) * 2018-09-07 2021-02-02 中国水利水电科学研究院 一种垄沟布局对降雨产流的影响机理识别方法
US20200229361A1 (en) * 2019-01-21 2020-07-23 James Canyon Method and System of Determining Soil-Water Properties
CN110264052B (zh) * 2019-06-06 2023-05-30 郑州市污水净化有限公司 城市河流水质提升系统

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6332287B1 (en) * 1997-03-06 2001-12-25 Carroll M. Geraldson Plant cultivation apparatus and method
CA2843086A1 (en) * 2013-02-19 2014-08-19 Furbish Company, Llc Vegetative roofing systems
CN108918827A (zh) * 2018-07-02 2018-11-30 中国水利水电科学研究院 一种识别多孔材料对土壤持水性影响的方法
CN108956948A (zh) * 2018-07-02 2018-12-07 中国水利水电科学研究院 一种多孔材料对平原区产流影响的识别方法
CN108982807A (zh) * 2018-07-02 2018-12-11 中国水利水电科学研究院 一种识别多孔材料对土壤持水性影响的动态监测系统

Non-Patent Citations (6)

* Cited by examiner, † Cited by third party
Title
Factors affecting the ability of extensive green roofs to reduce nutrient pollutants in rainfall runoff;Yongwei Gong等;《Science of the Total Environment》;20200825;第732卷;第1-14页 *
Porous-fiber module increases infiltration and reduces runoff;Zhenyu Lv等;《Agronomy Journal》;20201031;第112卷(第5期);第4420-4436页 *
基于土柱实验的绿色屋顶种植基质氮磷淋溶特征;何松等;《中国给水排水》;20210608;第1-11页 *
应用岩棉材料提高丘陵区经济林土壤水分保蓄能力;顾佳悦等;《干旱地区农业研究》;20200531;第38卷(第3期);第10-18页 *
沂河上游典型小流域磷元素输出特征分析;李振炜;《中国优秀博硕士学位论文全文数据库(硕士)工程科技Ⅰ辑》;20120815(第8期);B027-76 *
沂蒙山区典型小流域特殊降雨的磷素输出特征;李振炜等;《环境科学》;20120430;第33卷(第4期);第1152-1158页 *

Also Published As

Publication number Publication date
CN113552305A (zh) 2021-10-26

Similar Documents

Publication Publication Date Title
Van Driel et al. Denitrification of agricultural drainage using wood-based reactors
Zhao et al. Nitrogen fate and environmental consequence in paddy soil under rice-wheat rotation in the Taihu lake region, China
Li et al. Increased nitrogen use efficiencies as a key mitigation alternative to reduce nitrate leaching in north china plain
Salazar et al. Evaluation of the DRAINMOD–N II model for predicting nitrogen losses in a loamy sand under cultivation in south-east Sweden
Liu et al. Use of a new controlled-loss-fertilizer to reduce nitrogen losses during winter wheat cultivation in the Danjiangkou reservoir area of China
Gao et al. Effects of vegetation cover on phosphorus loss from a hillslope cropland of purple soil under simulated rainfall: a case study in China
Zhang et al. Advances and prospects of soil erosion research in the black soil region of Northeast China
Chen et al. Effects of terracing on soil properties in three key mountainous regions of China
CN104239729B (zh) 一种基于水循环的农田面源污染测算方法及装置
Xia et al. Is indirect N2O emission a significant contributor to the agricultural greenhouse gas budget? A case study of a rice paddy-dominated agricultural watershed in eastern China
Chen et al. Comparison of nitrogen loss weight in ammonia volatilization, runoff, and leaching between common and slow-release fertilizer in paddy field
Nicholls et al. Nutrients in subsurface and runoff waters of the Holland Marsh, Ontario
CN113552305B (zh) 一种多孔纤维材料埋设对产流水质影响的机理识别方法
Helwig et al. Modelling nitrate losses in drainage water using DRAINMOD 5.0
CN114936339B (zh) 一种提升农业面源污染防治效果的模拟方法
Guo et al. Performance analysis and experimental study on rainfall water purification with an extensive green roof matrix layer in Shanghai, China
Ale et al. Comparison of the performances of DRAINMOD-NII and ADAPT models in simulating nitrate losses from subsurface drainage systems
CN110231451B (zh) 一种旱生植被土壤适生水分的调查方法
CN104782301A (zh) 一种兼顾区域产量与环境风险的北京地区夏玉米分区施氮的方法
CN112352494A (zh) 一种利用生物炭、膨润土和有机肥改良风沙土的方法
CN104969681A (zh) 一种利用引黄泥沙快速改良粘质盐土的方法
CN104813789A (zh) 一种兼顾区域产量与环境风险的北京地区冬小麦分区施氮的方法
CN113673160B (zh) 一种面向作物水分高效利用的多孔材料优化埋设方法
WANG et al. Review of soil nitrate formation, leaching transport and their control measures in upland farming systems
CN113273370B (zh) 基于相关性分析和数据修正的精准施肥方法

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