CN107273686A - 雨洪径流氮磷输出负荷估算方法 - Google Patents
雨洪径流氮磷输出负荷估算方法 Download PDFInfo
- Publication number
- CN107273686A CN107273686A CN201710454167.9A CN201710454167A CN107273686A CN 107273686 A CN107273686 A CN 107273686A CN 201710454167 A CN201710454167 A CN 201710454167A CN 107273686 A CN107273686 A CN 107273686A
- Authority
- CN
- China
- Prior art keywords
- mrow
- data
- msub
- output load
- history
- 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
Links
Classifications
-
- 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
Abstract
本发明涉及雨洪径流氮磷输出负荷估算方法，包括：选定目标流域，并选定该流域出口处的监测断面；获取历史数据；根据相应数据拟合数学模型；验证数学模型的模拟效果；选择模拟效果更好的数学模型作为预测模型用于估算。本发明针对性强，且能根据实际情况选择更加合适的估算模型，科学合理。
Description
技术领域
本发明涉及一种雨洪径流氮磷输出负荷估算方法，属于水文监测技术领域。
背景技术
降雨是流域面源污染发生的主要驱动因子，面源污染物一般随流域雨洪径流过程输出。我国东部位于东亚季风气候区，雨洪过程主要集中在6月至9月。研究氮磷随雨洪径流过程输出特征对于流域氮磷污染控制具有重要意义。由于面源污染的发生和降雨径流过程有关，在开展面源污染输出负荷估算时需要同时监测水量和水质，这需要投入大量的人力物力和财力。若能在已有监测资料基础上，利用水文监测要素与氮磷输出负荷之间的关系，建立基于水文要素的次洪过程氮磷输出负荷预测模型，估算次洪过程氮磷输出负荷，可以节省面源污染负荷估算成本。
经检索发现，专利号CN201110109508.1、授权公告号CN102201034B、名称《一种临河村落非点源污染负荷的估算方法》的中国发明专利，其方法包括：1)基础资料收集：获取某地小雨、中雨、大到暴雨等三种情况的多年累积平均雨量R1、R2、R3，以及村庄面积A等基本情况；2)野外实地的监测：对有效降雨场次的降雨历时、降雨量和降雨径流的次径流系数、径流历时、产流面积进行监测，收集的水样进行化学需氧量、总氮、总磷的监测；3)非点源污染负荷的模型估算：采用公式计算临河村落非点源污染物化学需氧量、总氮、总磷的年污染负荷。
然而，以上述技术方案为代表的现有技术并没有提供专门针对雨洪径流的估算方法，亟待加以研究。
发明内容
本发明所要解决的技术问题是：克服现有技术存在的问题，提供一种雨洪径流氮磷输出负荷估算方法，针对性强，且能根据实际情况选择更加合适的估算模型。
本发明解决其技术问题的技术方案如下：
一种雨洪径流氮磷输出负荷估算方法，其特征是，包括以下步骤：
第一步、选定目标流域，并选定该流域出口处的监测断面；选定历史时期，获取目标流域在该历史时期内的历史降水量数据，并获取该监测断面在该历史时期内的历史流量数据、总氮浓度数据以及总磷浓度数据；
第二步、根据历史降水量数据，筛选出次降水事件；所述次降水事件指单次且连续的降水事件；所述次降水事件中，从降水开始的时刻至降水结束的时刻，为次降水事件期间；
第三步、在第二步所得的所有次降水事件期间：
根据相应的历史流量数据、总氮浓度数据以及总磷浓度数据，计算出总氮输出负荷数据以及总磷输出负荷数据；根据历史流量数据按预定算法计算出历史直接径流量数据；将历史降水量数据与总氮输出负荷数据、历史降水量数据与总磷输出负荷数据、历史直接径流量数据与总氮输出负荷数据、历史直接径流量数据与总磷输出负荷数据分别一一对应；
第四步、根据第三步得出的所有对应数据，拟合数学模型：
将历史降水量数据与总氮输出负荷数据、历史降水量数据与总磷输出负荷数据分别拟合为总氮输出负荷、总磷输出负荷与降水量相关的二次多项式模型；
将历史直接径流量数据与总氮输出负荷数据、历史直接径流量数据与总磷输出负荷数据分别拟合为总氮输出负荷、总磷输出负荷与直接径流量相关的线性模型；
其中，与总氮输出负荷数据相关的二次多项式模型和线性模型属于总氮负荷预测模型，与总磷输出负荷数据相关的二次多项式模型和线性模型属于总磷负荷预测模型；
第五步、根据历史降水量数据，筛选出不同于第二步的次降水事件；
在这些次降水事件期间：
根据相应的历史流量数据、总氮浓度数据以及总磷浓度数据，计算出总氮输出负荷数据以及总磷输出负荷数据，作为实际数据；根据历史流量数据按预定算法计算出历史直接径流量数据；将历史降水量数据、历史直接径流量数据分别代入到第四步所得的数学模型中，并计算得出相应的总氮输出负荷数据以及总磷输出负荷数据，作为各数学模型的预测数据；
第六步、针对各数学模型，将其预测数据与实际数据进行比较，并计算得出评价因子，所述评价因子包括相对平均误差、决定系数和模拟效率系数；根据评价因子，分别在总氮负荷预测模型、总磷负荷预测模型中选择模拟效果更好的数学模型作为最终的总氮负荷预测模型和总磷负荷预测模型；
第七步、将第六步所得的总氮负荷预测模型和总磷负荷预测模型用于估算。
发明人经深入实践研究发现，不同的流域，估算其雨洪径流氮磷输出负荷往往需要不同的预测模型；在此发现的基础上，发明人经进一步地实践研究，终于得出了上述方法，先分别拟合与降水量相关的二次多项式模型和与直接径流量相关的线性模型，再从中选择模拟效果更好的模型用于估算，即能专门针对雨洪径流，并实现良好的估算效果。
本发明进一步完善的技术方案如下：
优选地，第一步中，所述历史时期为连续N年，N为大于或等于3的整数，即从第一年至第N年；第二步中，从第一年起至少两年中筛选出次降水事件；第五步中，从其余历史时期中筛选出次降水事件。
具体而言，第一步中，所述历史时期为连续三年，即第一年、第二年以及第三年；第二步中，从第一年和第二年中筛选出次降水事件；第五步中，从第三年中筛选出次降水事件。
采用该优选方案，能进一步确保数学模型的模拟效果。
优选地，第三步及第五步中，所述预定算法为按以下公式进行计算：
qf(t)＝q(t)-qb(t)
且
其中，qb(t)为t时刻的基流；qb(t-1)为t-1时刻的基流；q(t)为t时刻实测径流量，且q(t)来源于历史流量数据；t为时间，单位为d；qf(t)为t时刻的直接径流量，且qf(t)构成历史直接径流量数据；BFImax为最大基流指数；α为退水常数；
当t＝1时，qb(t-1)＝所选历史时期内历史流量数据的最小值×调整系数，1.1≤调整系数≤1.5。
具体而言，0<BFImαx<1，且0<α<1。
在具体应用时，BFImax＝0.97，α＝0.8。
采用该优选方案，可更加简单、准确地将各时刻流量分割为直接径流量和基流量，进而利于确保与直接径流量相关数学模型的模拟效果。
优选地，第四步中，二次多项式模型为：L＝A+B·P+C·P2，其中，L为总氮或总磷输出负荷，P为降水量，A、B、C为常数；
线性模型为：L＝A+B·DR，其中，L为总氮或总磷输出负荷，DR为直接径流量，A、B为常数。
采用该优选方案，可进一步简化数学模型，利于提高整个估算方法的效率。
优选地，第六步中，所述评价因子由相对平均误差RME、决定系数R2和模拟效率系数Nash-Sutcliffe组成。
相对平均误差RME按下式计算：
决定系数R2按下式计算：
模拟效率系数Nash-Sutcliffe按下式计算：
其中，Poi为第i次观测得到的负荷值，为全部观测负荷值的平均值，Psi为由模型计算得到与第i次观测对应的负荷值，为由模型计算所得全部负荷值的平均值，n为总观测次数。
此外，第六步中，模拟效果的评价过程为：
S1、评价原则为：相对平均误差的绝对值接近零者为佳；决定系数的数值接近1者为佳；模拟效率系数的数值接近1者为佳；
S2、评价同属于总氮或总磷负荷预测模型的二次多项式模型和线性模型：比较两模型的相对平均误差、决定系数以及模拟效率系数；若其中一个模型有两个以上评价因子均更佳，则该模型的模拟效果更好。
采用该优选方案，可快速得出评价结果，利于提高整个估算方法的效率。
与现有技术相比，本发明针对性强，且能根据实际情况选择更加合适的估算模型，科学合理。
附图说明
图1至图4为试验案例中的比较结果图。
具体实施方式
下面参照附图并结合实施例对本发明作进一步详细描述。但是本发明不限于所给出的例子。
本发明具体实施的雨洪径流氮磷输出负荷估算方法，包括：
第一步、选定目标流域，并选定该流域出口处的监测断面；选定历史时期，获取目标流域在该历史时期内的历史降水量数据，并获取该监测断面在该历史时期内的历史流量数据、总氮浓度数据以及总磷浓度数据。
第二步、根据历史降水量数据，筛选出次降水事件；次降水事件指单次且连续的降水事件；次降水事件中，从降水开始的时刻至降水结束的时刻，为次降水事件期间。
第三步、在第二步所得的所有次降水事件期间：
根据相应的历史流量数据、总氮浓度数据以及总磷浓度数据，计算出总氮输出负荷数据以及总磷输出负荷数据；根据历史流量数据按预定算法计算出历史直接径流量数据；将历史降水量数据与总氮输出负荷数据、历史降水量数据与总磷输出负荷数据、历史直接径流量数据与总氮输出负荷数据、历史直接径流量数据与总磷输出负荷数据分别一一对应。
第四步、根据第三步得出的所有对应数据，拟合数学模型：
将历史降水量数据与总氮输出负荷数据、历史降水量数据与总磷输出负荷数据分别拟合为总氮输出负荷、总磷输出负荷与降水量相关的二次多项式模型；
将历史直接径流量数据与总氮输出负荷数据、历史直接径流量数据与总磷输出负荷数据分别拟合为总氮输出负荷、总磷输出负荷与直接径流量相关的线性模型；
其中，与总氮输出负荷数据相关的二次多项式模型和线性模型属于总氮负荷预测模型，与总磷输出负荷数据相关的二次多项式模型和线性模型属于总磷负荷预测模型。
具体而言，二次多项式模型为：L＝A+B·P+C·P2，其中，L为总氮或总磷输出负荷，P为降水量，A、B、C为常数；
线性模型为：L＝A+B·DR，其中，L为总氮或总磷输出负荷，DR为直接径流量，A、B为常数。
第五步、根据历史降水量数据，筛选出不同于第二步的次降水事件；
在这些次降水事件期间：
根据相应的历史流量数据、总氮浓度数据以及总磷浓度数据，计算出总氮输出负荷数据以及总磷输出负荷数据，作为实际数据；根据历史流量数据按预定算法计算出历史直接径流量数据；将历史降水量数据、历史直接径流量数据分别代入到第四步所得的数学模型中，并计算得出相应的总氮输出负荷数据以及总磷输出负荷数据，作为各数学模型的预测数据；
第六步、针对各数学模型，将其预测数据与实际数据进行比较，并计算得出评价因子，评价因子包括相对平均误差、决定系数和模拟效率系数；根据评价因子，分别在总氮负荷预测模型、总磷负荷预测模型中选择模拟效果更好的数学模型作为最终的总氮负荷预测模型和总磷负荷预测模型。
具体而言，评价因子由相对平均误差RME、决定系数R2和模拟效率系数Nash-Sutcliffe组成。
相对平均误差RME按下式计算：
决定系数R2按下式计算：
模拟效率系数Nash-Sutcliffe按下式计算：
其中，Poi为第i次观测得到的负荷值，为全部观测负荷值的平均值，Psi为由模型计算得到与第i次观测对应的负荷值，为由模型计算所得全部负荷值的平均值，n为总观测次数。
模拟效果的评价过程为：
S1、评价原则为：相对平均误差的绝对值接近零者为佳；决定系数的数值接近1者为佳；模拟效率系数的数值接近1者为佳；
S2、评价同属于总氮或总磷负荷预测模型的二次多项式模型和线性模型：比较两模型的相对平均误差、决定系数以及模拟效率系数；若其中一个模型有两个以上评价因子均更佳，则该模型的模拟效果更好。
第七步、将第六步所得的总氮负荷预测模型和总磷负荷预测模型用于估算。
其中，第一步中，历史时期为连续N年，N为大于或等于3的整数，即从第一年至第N年；第二步中，从第一年起至少两年中筛选出次降水事件；第五步中，从其余历史时期中筛选出次降水事件。
比如，第一步中，历史时期为连续三年，即第一年、第二年以及第三年；第二步中，从第一年和第二年中筛选出次降水事件；第五步中，从第三年中筛选出次降水事件。
第三步及第五步中，预定算法为按以下公式进行计算：
qf(t)＝q(t)-qb(t)
且
其中，qb(t)为t时刻的基流；qb(t-1)为t-1时刻的基流；q(t)为t时刻实测径流量，且q(t)来源于历史流量数据；t为时间，单位为d；qf(t)为t时刻的直接径流量，且qf(t)构成历史直接径流量数据；BFImax为最大基流指数；α为退水常数。当t＝1时，qb(t-1)＝所选历史时期内历史流量数据的最小值×调整系数，1.1≤调整系数≤1.5。
0<BFImax<1，且0<α<1。
注：根据水文学基本原理，将总径流分割为基流和直接径流，其中，基流包括深层壤中流和地下径流，直接径流流包括表层径流和浅层壤中流。
以下为具体试验案例。
以花某山流域为目标流域，并选定其出口处的某监测断面。选定历史时期为2012-2014年，通过当地水文局获取：该历史时期内的历史降水量数据，该监测断面在该历史时期内的历史流量数据、总氮浓度数据以及总磷浓度数据。
筛选出2012年、2013年的全部次降水事件，在这些次降水事件期间，根据相应的历史流量数据、总氮浓度数据以及总磷浓度数据，计算出总氮输出负荷数据以及总磷输出负荷数据；根据历史流量数据按预定算法计算出历史直接径流量数据(BFImax＝0.97，α＝0.8)；将历史降水量数据与总氮输出负荷数据、历史降水量数据与总磷输出负荷数据、历史直接径流量数据与总氮输出负荷数据、历史直接径流量数据与总磷输出负荷数据分别一一对应。
然后根据这些对应数据，拟合数学模型如下：
L总氮＝-8.84+0.25×P+0.0015×P2
L总磷＝8.61-0.1697×P+0.00125×P2
L总氮＝-2.47+1.78×DR
L总磷＝0.34+0.27×DR。
之后，筛选出2014年的全部次降水事件，在这些次降水事件期间：
根据相应的历史流量数据、总氮浓度数据以及总磷浓度数据，计算出总氮输出负荷数据以及总磷输出负荷数据，作为实际数据；
根据历史流量数据按预定算法计算出历史直接径流量数据；将历史降水量数据、历史直接径流量数据分别代入到上述数学模型中，并计算得出相应的总氮输出负荷数据以及总磷输出负荷数据，作为各数学模型的预测数据。
针对各数学模型，将其预测数据与实际数据进行比较，并计算得出评价因子，评价因子由相对平均误差RME、决定系数R2和模拟效率系数Nash-Sutcliffe组成。结果如图1至图4所示。
采用前文提到的模拟效果评价办法进行评价，可知，总氮负荷预测模型、总磷负荷预测模型均应采用与直接径流量相关的线性模型用于估算。
除上述实施例外，本发明还可以有其他实施方式。凡采用等同替换或等效变换形成的技术方案，均落在本发明要求的保护范围。
Claims (10)
1.一种雨洪径流氮磷输出负荷估算方法，其特征是，包括以下步骤：
第一步、选定目标流域，并选定该流域出口处的监测断面；选定历史时期，获取目标流域在该历史时期内的历史降水量数据，并获取该监测断面在该历史时期内的历史流量数据、总氮浓度数据以及总磷浓度数据；
第二步、根据历史降水量数据，筛选出次降水事件；所述次降水事件指单次且连续的降水事件；所述次降水事件中，从降水开始的时刻至降水结束的时刻，为次降水事件期间；
第三步、在第二步所得的所有次降水事件期间：
根据相应的历史流量数据、总氮浓度数据以及总磷浓度数据，计算出总氮输出负荷数据以及总磷输出负荷数据；根据历史流量数据按预定算法计算出历史直接径流量数据；将历史降水量数据与总氮输出负荷数据、历史降水量数据与总磷输出负荷数据、历史直接径流量数据与总氮输出负荷数据、历史直接径流量数据与总磷输出负荷数据分别一一对应；
第四步、根据第三步得出的所有对应数据，拟合数学模型：
将历史降水量数据与总氮输出负荷数据、历史降水量数据与总磷输出负荷数据分别拟合为总氮输出负荷、总磷输出负荷与降水量相关的二次多项式模型；
将历史直接径流量数据与总氮输出负荷数据、历史直接径流量数据与总磷输出负荷数据分别拟合为总氮输出负荷、总磷输出负荷与直接径流量相关的线性模型；
其中，与总氮输出负荷数据相关的二次多项式模型和线性模型属于总氮负荷预测模型，与总磷输出负荷数据相关的二次多项式模型和线性模型属于总磷负荷预测模型；
第五步、根据历史降水量数据，筛选出不同于第二步的次降水事件；
在这些次降水事件期间：
根据相应的历史流量数据、总氮浓度数据以及总磷浓度数据，计算出总氮输出负荷数据以及总磷输出负荷数据，作为实际数据；根据历史流量数据按预定算法计算出历史直接径流量数据；将历史降水量数据、历史直接径流量数据分别代入到第四步所得的数学模型中，并计算得出相应的总氮输出负荷数据以及总磷输出负荷数据，作为各数学模型的预测数据；
第六步、针对各数学模型，将其预测数据与实际数据进行比较，并计算得出评价因子，所述评价因子包括相对平均误差、决定系数和模拟效率系数；根据评价因子，分别在总氮负荷预测模型、总磷负荷预测模型中选择模拟效果更好的数学模型作为最终的总氮负荷预测模型和总磷负荷预测模型；
第七步、将第六步所得的总氮负荷预测模型和总磷负荷预测模型用于估算。
2.根据权利要求1述及的雨洪径流氮磷输出负荷估算方法，其特征是，第一步中，所述历史时期为连续N年，N为大于或等于3的整数，即从第一年至第N年；第二步中，从第一年起至少两年中筛选出次降水事件；第五步中，从其余历史时期中筛选出次降水事件。
3.根据权利要求2述及的雨洪径流氮磷输出负荷估算方法，其特征是，第一步中，所述历史时期为连续三年，即第一年、第二年以及第三年；第二步中，从第一年和第二年中筛选出次降水事件；第五步中，从第三年中筛选出次降水事件。
4.根据权利要求1述及的雨洪径流氮磷输出负荷估算方法，其特征是，第三步及第五步中，所述预定算法为按以下公式进行计算：
qf(t)＝q(t)-qb(t)
且
<mrow>
<msub>
<mi>q</mi>
<mrow>
<mi>b</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<msub>
<mi>BFI</mi>
<mrow>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
</mrow>
</msub>
<mo>)</mo>
<mo>&CenterDot;</mo>
<mi>&alpha;</mi>
<mo>&CenterDot;</mo>
<msub>
<mi>q</mi>
<mrow>
<mi>b</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
</msub>
<mo>+</mo>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<mi>&alpha;</mi>
<mo>)</mo>
<mo>&CenterDot;</mo>
<msub>
<mi>BFI</mi>
<mrow>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
</mrow>
</msub>
<mo>&CenterDot;</mo>
<msub>
<mi>q</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</msub>
</mrow>
<mrow>
<mn>1</mn>
<mo>-</mo>
<mi>&alpha;</mi>
<mo>&CenterDot;</mo>
<msub>
<mi>BFI</mi>
<mrow>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
</mrow>
</msub>
</mrow>
</mfrac>
</mrow>
其中，qb(t)为t时刻的基流；qb(t-1)为t-1时刻的基流；q(t)为t时刻实测径流量，且q(t)来源于历史流量数据；t为时间，单位为d；qf(t)为t时刻的直接径流量，且qf(t)构成历史直接径流量数据；BFImax为最大基流指数；α为退水常数；
当t＝1时，qb(t-1)＝所选历史时期内历史流量数据的最小值×调整系数，1.1≤调整系数≤1.5。
5.根据权利要求4述及的雨洪径流氮磷输出负荷估算方法，其特征是，0<BFImax<1，且0<α<1。
6.根据权利要求5述及的雨洪径流氮磷输出负荷估算方法，其特征是，BFImax＝0.97，α＝0.8。
7.根据权利要求1述及的雨洪径流氮磷输出负荷估算方法，其特征是，第四步中，二次多项式模型为：L＝A+B·P+C·P2，其中，L为总氮或总磷输出负荷，P为降水量，A、B、C为常数；
线性模型为：L＝A+B·DR，其中，L为总氮或总磷输出负荷，DR为直接径流量，A、B为常数。
8.根据权利要求1述及的雨洪径流氮磷输出负荷估算方法，其特征是，第六步中，所述评价因子由相对平均误差RME、决定系数R2和模拟效率系数Nash-Sutcliffe组成。
9.根据权利要求1述及的雨洪径流氮磷输出负荷估算方法，其特征是，相对平均误差RME按下式计算：
<mrow>
<mi>R</mi>
<mi>M</mi>
<mi>E</mi>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mi>n</mi>
</mfrac>
<msubsup>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</msubsup>
<mfrac>
<mrow>
<msub>
<mi>P</mi>
<mrow>
<mi>o</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>P</mi>
<mrow>
<mi>s</mi>
<mi>i</mi>
</mrow>
</msub>
</mrow>
<msub>
<mi>P</mi>
<mrow>
<mi>o</mi>
<mi>i</mi>
</mrow>
</msub>
</mfrac>
</mrow>
决定系数R2按下式计算：
<mrow>
<msup>
<mi>R</mi>
<mn>2</mn>
</msup>
<mo>=</mo>
<mfrac>
<msup>
<mrow>
<mo>&lsqb;</mo>
<msubsup>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</msubsup>
<mrow>
<mo>(</mo>
<msub>
<mi>P</mi>
<mrow>
<mi>o</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mover>
<mi>P</mi>
<mo>&OverBar;</mo>
</mover>
<mi>o</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>&times;</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>P</mi>
<mrow>
<mi>s</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mover>
<mi>P</mi>
<mo>&OverBar;</mo>
</mover>
<mi>s</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
</mrow>
<mn>2</mn>
</msup>
<mrow>
<msubsup>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</msubsup>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>P</mi>
<mrow>
<mi>o</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mover>
<mi>P</mi>
<mo>&OverBar;</mo>
</mover>
<mi>o</mi>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>&times;</mo>
<msubsup>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</msubsup>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>P</mi>
<mrow>
<mi>s</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mover>
<mi>P</mi>
<mo>&OverBar;</mo>
</mover>
<mi>s</mi>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
</mrow>
模拟效率系数Nash-Sutcliffe按下式计算：
<mrow>
<mi>N</mi>
<mi>a</mi>
<mi>s</mi>
<mi>h</mi>
<mi> </mi>
<mi>S</mi>
<mi>u</mi>
<mi>t</mi>
<mi>c</mi>
<mi>l</mi>
<mi>i</mi>
<mi>f</mi>
<mi>f</mi>
<mi>e</mi>
<mo>=</mo>
<mn>1</mn>
<mo>-</mo>
<mfrac>
<mrow>
<msubsup>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</msubsup>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>P</mi>
<mrow>
<mi>o</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>P</mi>
<mrow>
<mi>s</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
<mrow>
<msubsup>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</msubsup>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>P</mi>
<mrow>
<mi>o</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mover>
<mi>P</mi>
<mo>&OverBar;</mo>
</mover>
<mi>o</mi>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
</mrow>
其中，Poi为第i次观测得到的负荷值，为全部观测负荷值的平均值，Psi为由模型计算得到与第i次观测对应的负荷值，为由模型计算所得全部负荷值的平均值，n为总观测次数。
10.根据权利要求1述及的雨洪径流氮磷输出负荷估算方法，其特征是，第六步中，模拟效果的评价过程为：
S1、评价原则为：相对平均误差的绝对值接近零者为佳；决定系数的数值接近1者为佳；模拟效率系数的数值接近1者为佳；
S2、评价同属于总氮或总磷负荷预测模型的二次多项式模型和线性模型：比较两模型的相对平均误差、决定系数以及模拟效率系数；若其中一个模型有两个以上评价因子均更佳，则该模型的模拟效果更好。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710454167.9A CN107273686B (zh) | 2017-06-15 | 2017-06-15 | 雨洪径流氮磷输出负荷预测方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710454167.9A CN107273686B (zh) | 2017-06-15 | 2017-06-15 | 雨洪径流氮磷输出负荷预测方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107273686A true CN107273686A (zh) | 2017-10-20 |
| CN107273686B CN107273686B (zh) | 2018-08-17 |
Family
ID=60066795
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710454167.9A Active CN107273686B (zh) | 2017-06-15 | 2017-06-15 | 雨洪径流氮磷输出负荷预测方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107273686B (zh) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109117984B (zh) * | 2018-07-10 | 2020-06-19 | 上海交通大学 | 稻田径流量预测及氮磷流失估测方法 |
| CN112348241A (zh) * | 2020-10-28 | 2021-02-09 | 中国农业科学院农业资源与农业区划研究所 | 一种保护地农用地面源污染排放量预测方法及系统 |
| CN112541611A (zh) * | 2020-10-28 | 2021-03-23 | 中国农业科学院农业资源与农业区划研究所 | 一种雨养农用地面源污染排放量预测方法及系统 |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102201034A (zh) * | 2011-04-26 | 2011-09-28 | 北京师范大学 | 一种临河村落非点源污染负荷的估算方法 |
| CN102663267A (zh) * | 2012-05-15 | 2012-09-12 | 南京大学 | 一种半湿润区流域面源污染负荷的确定方法 |
| CN102677626A (zh) * | 2012-05-15 | 2012-09-19 | 临沂大学 | 一种构建小流域面源污染综合防治体系的方法 |
| CN103020424A (zh) * | 2012-11-22 | 2013-04-03 | 北京师范大学 | 一种基于人工降雨实验的北方平原农田区非点源污染负荷估算方法 |
| CN106202857A (zh) * | 2016-06-22 | 2016-12-07 | 中国水利水电科学研究院 | 一种地表水质对旱涝事件的响应分析方法 |
| CN106250666A (zh) * | 2016-06-22 | 2016-12-21 | 中国水利水电科学研究院 | 一种灌溉系统发达区湖泊氮磷污染负荷模拟方法 |
-
2017
- 2017-06-15 CN CN201710454167.9A patent/CN107273686B/zh active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102201034A (zh) * | 2011-04-26 | 2011-09-28 | 北京师范大学 | 一种临河村落非点源污染负荷的估算方法 |
| CN102663267A (zh) * | 2012-05-15 | 2012-09-12 | 南京大学 | 一种半湿润区流域面源污染负荷的确定方法 |
| CN102677626A (zh) * | 2012-05-15 | 2012-09-19 | 临沂大学 | 一种构建小流域面源污染综合防治体系的方法 |
| CN103020424A (zh) * | 2012-11-22 | 2013-04-03 | 北京师范大学 | 一种基于人工降雨实验的北方平原农田区非点源污染负荷估算方法 |
| CN106202857A (zh) * | 2016-06-22 | 2016-12-07 | 中国水利水电科学研究院 | 一种地表水质对旱涝事件的响应分析方法 |
| CN106250666A (zh) * | 2016-06-22 | 2016-12-21 | 中国水利水电科学研究院 | 一种灌溉系统发达区湖泊氮磷污染负荷模拟方法 |
Non-Patent Citations (6)
| Title |
|---|
| 任玮，等: "基于改进输出系数模型的云南宝象河流域非点源污染负荷估算", 《中国环境科学》 * |
| 余进祥，等: "鄱阳湖流域泥沙流失及吸附态氮磷输出负荷评估", 《生态学报》 * |
| 张丽，等: "流域非点源磷素输出负荷估算方法研究", 《水文》 * |
| 徐立红，等: "基于流域降雨强度的氮磷输出系数模型改进及应用", 《农业工程学报》 * |
| 马永胜，等著: "《水资源保护理论与实践》", 30 June 2009, 中国水利水电出版社 * |
| 龙天渝，等: "基于SLURP模型和输出系数法的三峡库区非点源氮磷负荷预测_基于SLURP模型和输出系数法的三峡库区非点源氮磷负荷预测", 《环境科学学报》 * |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109117984B (zh) * | 2018-07-10 | 2020-06-19 | 上海交通大学 | 稻田径流量预测及氮磷流失估测方法 |
| CN112348241A (zh) * | 2020-10-28 | 2021-02-09 | 中国农业科学院农业资源与农业区划研究所 | 一种保护地农用地面源污染排放量预测方法及系统 |
| CN112541611A (zh) * | 2020-10-28 | 2021-03-23 | 中国农业科学院农业资源与农业区划研究所 | 一种雨养农用地面源污染排放量预测方法及系统 |
| CN112348241B (zh) * | 2020-10-28 | 2024-01-05 | 中国农业科学院农业资源与农业区划研究所 | 一种保护地农用地面源污染排放量预测方法及系统 |
| CN112541611B (zh) * | 2020-10-28 | 2024-03-29 | 中国农业科学院农业资源与农业区划研究所 | 一种雨养农用地面源污染排放量预测方法及系统 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107273686B (zh) | 2018-08-17 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Chen et al. | Urban flood risk warning under rapid urbanization | |
| Trepel et al. | Quantifying nitrogen retention in surface flow wetlands for environmental planning at the landscape-scale | |
| JP5537883B2 (ja) | 災害予測システム及び災害予測方法 | |
| CN107273686B (zh) | 雨洪径流氮磷输出负荷预测方法 | |
| CN102201034A (zh) | 一种临河村落非点源污染负荷的估算方法 | |
| CN103810532B (zh) | 优化城市排水系统运行状况的方法 | |
| CN106706033A (zh) | 一种海绵城市性能监视系统和方法 | |
| CN112163346B (zh) | 一种基于透水-不透水格局的面源污染迁移能力估算方法 | |
| CN111815955B (zh) | 一种基于交通流基本图的城市道路积水智能识别方法 | |
| CN115186858A (zh) | 基于不同影响类型的变电站洪水淹没风险预警方法及系统 | |
| CN114511990A (zh) | 泥石流概率测算方法、泥石流多要素协同监测预警方法 | |
| CN107169589A (zh) | 一种基于knn和轮盘赌算法的低压网格精细化负荷预测方法 | |
| CN109033589B (zh) | 一种基于swmm模型考虑lid净化作用的污染物去除率计算方法 | |
| CN112905721B (zh) | 一种海绵城市建设项目的数据库系统 | |
| Jayasooriya et al. | Development of a framework for the valuation of Eco-System Services of Green Infrastructure | |
| Eric et al. | Modeling low impact development potential with hydrological response units | |
| CN108364090B (zh) | 一种预测城市道路表面苯系物累积负荷的方法 | |
| Liu et al. | Study on LID facilities comprehensive effect evaluation: a case in campus | |
| CN111047213B (zh) | 基于多要素联合判别的中长期水资源量临界预警指标方法 | |
| CN113792367B (zh) | 基于PySWMM的排水系统多来源入流入渗和出渗量动态估算方法 | |
| CN112528563B (zh) | 一种基于svm算法的城市内涝预警方法 | |
| CN112163347B (zh) | 一种城市地块尺度面源污染动态输出特征评估的方法 | |
| CN112330150A (zh) | 大区域泥石流灾害可发性的早期判识方法 | |
| KR20050112700A (ko) | 비점오염원 모니터링 방법 및 시스템 | |
| McGarity | Optimization for urban watershed management: Stormwater runoff and nonpoint pollution control |
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 |