CN107621532A - 基于突变分析确定浅水湖泊底泥起悬临界切应力的方法 - Google Patents
基于突变分析确定浅水湖泊底泥起悬临界切应力的方法 Download PDFInfo
- Publication number
- CN107621532A CN107621532A CN201710865987.7A CN201710865987A CN107621532A CN 107621532 A CN107621532 A CN 107621532A CN 201710865987 A CN201710865987 A CN 201710865987A CN 107621532 A CN107621532 A CN 107621532A
- Authority
- CN
- China
- Prior art keywords
- mrow
- shear stress
- curve
- critical shear
- outstanding
- 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.)
- Pending
Links
Landscapes
- Measuring Or Testing Involving Enzymes Or Micro-Organisms (AREA)
Abstract
本发明公开了基于突变分析确定浅水湖泊底泥起悬临界切应力的方法,包括如下步骤:获取数据,包括水‑土界面处悬浮物浓度与底层切应力;通过曲线拟合,根据上述数据建立底层切应力与悬浮物浓度的拟合函数,根据拟合曲线划分初步判定临界切应力的范围;应用非参数Mann‑Kendall检验法对上述数据进行突变分析,确定底泥起悬的临界切应力。本发明简单易操作,图像直观清晰,客观、科学、定量地确定了湖泊底泥起悬临界切应力,为湖泊沉积物再悬浮与内源释放理论研究、湖泊富营养化治理提供可靠依据。具有很强的实用性和广泛的适用性。
Description
技术领域
本发明涉及一种湖泊底泥起悬临界切应力的方法,具体涉及一种基于突变分析确定浅水湖泊底泥起悬临界切应力的方法。
背景技术
当前湖泊面临着诸多问题,如富营养化、有机污染、西部湖泊咸化、生态环境以及湖泊萎缩、水量减少等,尤其是富营养化,不但制约了湖泊资源的可利用性,而且直接影响到人类的健康生存与社会经济的持续发展。
湖泊富营养化是指由于水体中氮、磷营养物质的富集,引起藻类及其他浮游生物迅速繁殖、水体溶解氧量下降鱼类及其他生物大量死亡、水质恶化的现象。从富营养污染控制分析,可将污染源分为外源和内源。外源污染包括工业废水、城镇生活污水、固废处置场、城镇地表径流、农牧区地表径流、大气降水等,内源污染指的是入湖的各种营养物、污染物的主要蓄积而成的污染底泥等。湖泊外源污染的控制及治理目前已取得了一定实效。外源得到控制后,内源的影响会变得明显,底泥的释放速率加快,使水体的污染状况持续下去。
沉积物再悬浮是底边界层内频繁发生的动力过程,是实现沉积物与水体物质交换的重要途径,在沉积物与水体的物质交换过程中有着重要作用。水动力引起的沉积物再悬浮会影响到水下光场分布、浊度以及内源释放,进而影响到蓝、绿藻的生长。当水-土界面的底泥受到风浪的频繁扰动时,会发生再悬浮,使得底泥间隙水中高浓度的营养释放出来,进一步加剧了内源的负荷。
众多研究表明动力扰动引起的底泥起悬及内部营养物的释放量远大于浓度梯度的静释放量,在水体生态环境的演变中扮演了重要角色。底泥悬浮深度的大小受控于风浪的扰动强度,扰动越剧烈,水体底部产生的波浪切应力越大,底泥悬浮量越大。
因此,湖泊底泥起悬临界切应力的确定,可为湖泊沉积物再悬浮与内源释放理论研究、湖泊富营养化治理提供可靠依据。
发明内容
为解决现有技术的不足,本发明的目的在于提供一种简单易操作,图像直观清晰的基于突变分析确定浅水湖泊底泥起悬临界切应力的方法。
为了实现上述目标,本发明采用如下的技术方案:
基于突变分析确定浅水湖泊底泥起悬临界切应力的方法,包括如下步骤:
S1、获取水-土界面处悬浮物浓度与底层切应力的数据;
S2、通过曲线拟合,根据上述数据建立底层切应力与悬浮物浓度的拟合函数,根据拟合曲线划分初步判定临界切应力的范围;
S3、应用非参数Mann-Kendall检验法进行突变分析,确定底泥起悬的临界切应力。
上述步骤S1中水-土界面处悬浮物浓度与底层切应力的数据通过野外监测实验同步获取。测量时,数据应同步测量,且测量采样数据频率保持高频率,以确保数据量充足,且保证分析的精确度。
上述步骤S2中曲线拟合通过SPSS、EXCEL、MATLAB数值分析软件完成。
上述步骤S2中建立底层切应力与悬浮物浓度的拟合函数,具体步骤为:选择回归模型进行数据拟合,通过曲线的拟合优度检验,确定拟合度最佳的函数。
进一步的,上述回归模型包括Linear线性模型、Quadratic二次模型、Compound复合模型、Growth生长模型、Logarithmic对数模型、S型模型、Cubic抛物线模型、Exponential指数模型、Inverse倒数模型、Power幂函数模型、Logistic模型。
进一步的,上述拟合优度检验的方法如下:
定义指数R2,表示误差平方和与均方差的比值,其计算公式为:
式中,Y表示实际值,Y1表示实际平均值,Y2表示拟合曲线所求得的理论值。
若误差平方和与均方差的比值越小,则说明实际观察值与估计值越接近,曲线拟合的越好。
上述步骤S2中根据拟合曲线划分初步判定临界切应力的范围,具体方法为:根据曲线斜率k的变化,将曲线划分为悬浮物浓度基本不变、悬浮物浓度随切应力变大而逐渐变大、悬浮物浓度随切应力变大迅速变大三部分,分别对应底泥未起悬、小规模起悬大规模起悬三阶段。
此方法可确定临界切应力的范围,下一步分析精确的临界切应力值。
上述步骤S3中应用非参数Mann-Kendall检验法进行突变分析,最终确定底泥起悬的临界切应力,具体步骤为:
A1:将获取的悬浮物浓度与底层切应力的同步数据,根据底层切应力大小升序排列,构造一组悬浮物浓度序列xi(i=1,2,…,n);
A2:构造一秩序列aij,表示xi>xj(1≤j≤i)的样本累积数,如下:
定义Sk:
其中k=2,3,…,n (3)
Sk均值E(Sk)以及方差var(Sk)定义如下:
定义统计量:
其中k=2,3,…,n (6)
其中UF1=0;
A3:将序列xi逆序排列,即xn,xn-1,…,x1,重复步骤A2,得出统计量UFk *令UBk=UFk *,给出显著水平α,取α=0.05,查正态分布表,得临界值u0.05=±1.96。
A4:将UFk和UBk两个统计量序列曲线和信度线(±1.96两条直线,上限称为LU,下限称为LD)绘在同一张图上,分析绘出的UFk和UBk曲线图,当UFk大于0时,表明序列呈上升趋势;若大于1.96,则表明序列上升趋势显著;当UFk小于0时,表明序列呈下降趋势;若小于-1.96,则表明序列下降趋势显著;当UFk曲线与UBk曲线出现交点时,且交点处于±1.96之间,则表明存在突变点,且交点即为突变点;
将UFk曲线与UBk曲线交点所对应的切应力定义为τ1,为底泥由未起悬向小规模起悬状态过渡的临界切应力,UFk曲线与LU直线交点所对应的切应力定义为τ2,为底泥由小规模起悬状态向大规模起悬状态过渡的临界切应力,并结合拟合曲线结果分析临界切应力取值的合理性。
本发明的有益之处在于:
本发明的基于突变分析确定浅水湖泊底泥起悬临界切应力的方法,简单易操作,图像直观清晰,客观、科学、定量地确定了湖泊底泥起悬临界切应力,为湖泊沉积物再悬浮与内源释放理论研究、湖泊富营养化治理提供可靠依据,具有很强的实用性和广泛的适用性。
附图说明
图1是本发明实施例中底层悬浮物浓度随波浪切应力变化的拟合曲线。
图2是本发明实施例中底层悬浮物浓度随切应力变化的M-K统计量曲线图。
具体实施方式
以下结合附图和具体实施例对本发明作具体的介绍。
基于突变分析确定浅水湖泊底泥起悬临界切应力的方法,包括如下步骤:
S1、同步获取水-土界面处悬浮物浓度与底层切应力的数据,采样频率为10Hz;
S2、通过SPSS、EXCEL、MATLAB数值分析软件完成曲线拟合,根据上述数据建立底层切应力与悬浮物浓度的拟合函数:选择回归模型,例如Linear线性模型、Quadratic二次模型、Compound复合模型、Growth生长模型、Logarithmic对数模型、S型模型、Cubic抛物线模型、Exponential指数模型、Inverse倒数模型、Power幂函数模型、Logistic模型进行拟合,通过曲线的拟合优度检验,确定拟合度最佳的函数;再根据拟合曲线划分初步判定临界切应力的范围;
拟合优度检验方法为:
定义相关指数R2,表示误差平方和与均方差的比值,其计算公式为:
式中,Y表示实际值,Y1表示实际平均值,Y2表示拟合曲线所求得的理论值。
根据拟合曲线划分初步判定临界切应力的范围的方法为:根据曲线斜率k的变化,划分为悬浮物浓度基本不变、悬浮物浓度随切应力变大而逐渐变大、悬浮物浓度随切应力变大迅速变大三部分,分别对应底泥未起悬、小规模起悬、中等规模及大规模起悬三阶段。
S3、应用非参数Mann-Kendall检验法进行突变分析,确定底泥起悬的临界切应力,具体步骤为:
A1:已获取悬浮物浓度与底层切应力的同步数据,根据底层切应力大小升序排列,构造一组悬浮物浓度序列xi(i=1,2,…,n)
A2:构造一秩序列aij,表示xi>xj(1≤j≤i)的样本累积数,如下:
定义Sk:
其中k=2,3,…,n (3)
Sk均值E(Sk)以及方差var(Sk)定义如下:
定义统计量:
其中k=2,3,…,n(6)
其中UF1=0;
A3:将序列xi逆序排列,即xn,xn-1,…,x1,重复步骤A2,得出统计量UFk *,令UBk=-UFk *,给出显著水平α,一般取α=0.05,查正态分布表,得临界值u0.05=±1.96。
A4:将UFk和UBk两个统计量序列曲线和信度线(±1.96两条直线,上限称为LU,下限称为LD)绘在同一张图上,分析绘出的UFk和UBk曲线图,当UFk大于0时,表明序列呈上升趋势;若大于1.96,则表明序列上升趋势显著;当UFk小于0时,表明序列呈下降趋势;若小于-1.96,则表明序列下降趋势显著;当UFk曲线与UBk曲线出现交点时,且交点处于±1.96之间,则表明存在突变点,且交点即为突变点。
将UFk曲线与UBk曲线交点所对应的切应力定义为τ1,为底泥由未起悬向小规模起悬状态过渡的临界切应力,UFk曲线与LU直线交点所对应的切应力定义为τ2,为底泥由小规模起悬状态向中等规模及大规模起悬状态过渡的临界切应力;并结合拟合曲线结果分析临界切应力取值的合理性。
下面以取自太湖的野外同步监测数据为例,说明本发明实施例方法的处理过程。
1、通过ADV(Son Tek;10Hz)、OBS(每3min输出一个值)获取水-土界面处悬浮物浓度与底层切应力的同步高分辨率数据,五分钟一组,共364组数据。
2、通过曲线拟合,建立底层切应力与悬浮物浓度的最佳拟合函数,根据拟合曲线划分初步判定临界切应力的范围。
底层悬浮物浓度随波浪切应力变化的拟合曲线如附图1所示,利用SPSS软件进行曲线拟合,最佳拟合函数为二次多项式,R2=0.7956:
根据曲线斜率k的变化,划分为悬浮物浓度基本不变、悬浮物浓度随切应力变大而逐渐变大、悬浮物浓度随切应力变大迅速变大三部分,分别对应底泥未起悬、小规模起悬、中等规模及大规模起悬三阶段。可确定底泥由未起悬向小规模起悬状态过渡的临界切应力范围为0.01~0.1N/m2,底泥由小规模起悬状态向大规模起悬状态过渡的临界切应力范围为0.05~0.5N/m2,精确的临界切应力值需要下一步分析。
3、应用非参数Mann-Kendall检验法进行突变分析,最终确定底泥起悬的临界切应力。
底层悬浮物浓度随切应力变化的M-K统计量曲线图如附图2所示,在±1.96置信范围内UF曲线与UB曲线存在交点,说明底层悬浮物浓度在不断增大的过程中存在突变情况,突变点对应的切应力的值为0.062N/m2。
此外单独对UF曲线的变化趋势进行分析,由附图2可知:
当波浪切应力小于0.062N/m2,UF曲线上下波动,其数值在绝大多数情况下主要集中在-1.96~0之间,表明底层悬浮物浓度呈现下降及显著下降趋势;
当波浪切应力为0.062~0.105N/m2时,UF值为正但小于1.96,表明底层底层悬浮物浓度呈现上升趋势;
当波浪切应力大于0.105N/m2时,UF值大于1.96,表明底层底层悬浮物浓度具有显著上升趋势。
结合底泥起悬经历的不同状态及切应力的取值范围,本次监测过程中底泥起悬可分为以下几个阶段:
底泥未起悬,τw<0.062N/m2,底层悬浮物浓度平均值为40.79±11.22mg/l;
底泥小规模起悬,0.062≤τw<0.105N/m2,底层悬浮物浓度平均值为53.78±10.46mg/l;
底泥中等规模及大规模起悬,τw≥0.105N/m2,底层悬浮物浓度平均值为148.33±63.76mg/l。
结合拟合曲线结果,临界切应力取值合理。
以上显示和描述了本发明的基本原理、主要特征和优点。本行业的技术人员应该了解,上述实施例不以任何形式限制本发明,凡采用等同替换或等效变换的方式所获得的技术方案,均落在本发明的保护范围内。
Claims (8)
1.基于突变分析确定浅水湖泊底泥起悬临界切应力的方法,其特征在于,包括如下步骤:
S1、获取数据,包括水-土界面处悬浮物浓度与底层切应力;
S2、通过曲线拟合,根据上述数据建立底层切应力与悬浮物浓度的拟合函数,根据拟合曲线划分初步判定临界切应力的范围;
S3、应用非参数Mann-Kendall检验法对上述数据进行突变分析,确定底泥起悬的临界切应力。
2.根据权利要求1所述的基于突变分析确定浅水湖泊底泥起悬临界切应力的方法,其特征在于,所述步骤S1中水-土界面处悬浮物浓度与底层切应力的数据通过野外监测实验同步获取。
3.根据权利要求1所述的基于突变分析确定浅水湖泊底泥起悬临界切应力的方法,其特征在于,所述步骤S2中曲线拟合通过SPSS、EXCEL、MATLAB数值分析软件完成。
4.根据权利要求1所述的基于突变分析确定浅水湖泊底泥起悬临界切应力的方法,其特征在于,所述步骤S2中建立底层切应力与悬浮物浓度的拟合函数,具体步骤为:选择回归模型进行数据拟合,通过曲线的拟合优度检验,确定拟合度最佳的函数。
5.根据权利要求4所述的基于突变分析确定浅水湖泊底泥起悬临界切应力的方法,其特征在于,所述回归模型包括Linear线性模型、Quadratic二次模型、Compound复合模型、Growth生长模型、Logarithmic对数模型、S型模型、Cubic抛物线模型、Exponential指数模型、Inverse倒数模型、Power幂函数模型、Logistic模型。
6.根据权利要求4所述的基于突变分析确定浅水湖泊底泥起悬临界切应力的方法,其特征在于,所述拟合优度检验的方法如下:
定义指数R2,表示误差平方和与均方差的比值,其计算公式为:
<mrow>
<msup>
<mi>R</mi>
<mn>2</mn>
</msup>
<mo>=</mo>
<mn>1</mn>
<mo>-</mo>
<mfrac>
<mrow>
<mi>&Sigma;</mi>
<msup>
<mrow>
<mo>(</mo>
<mi>Y</mi>
<mo>-</mo>
<msub>
<mi>Y</mi>
<mn>2</mn>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
<mrow>
<mi>&Sigma;</mi>
<msup>
<mrow>
<mo>(</mo>
<mi>Y</mi>
<mo>-</mo>
<msub>
<mi>Y</mi>
<mn>1</mn>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
式中,Y表示实际值,Y1表示实际平均值,Y2表示拟合曲线所求得的理论值。
7.根据权利要求1所述的基于突变分析确定浅水湖泊底泥起悬临界切应力的方法,其特征在于,所述步骤S2中根据拟合曲线划分初步判定临界切应力的范围,具体方法为:根据曲线斜率k的变化,将曲线划分为悬浮物浓度基本不变、悬浮物浓度随切应力变大而逐渐变大、悬浮物浓度随切应力变大迅速变大三部分,分别对应底泥未起悬、小规模起悬及大规模起悬三阶段。
8.根据权利要求1所述的基于突变分析确定浅水湖泊底泥起悬临界切应力的方法,其特征在于,所述步骤S3中应用非参数Mann-Kendall检验法进行突变分析,最终确定底泥起悬的临界切应力,具体步骤为:
A1:将获取的悬浮物浓度与底层切应力的同步数据,根据底层切应力大小升序排列,构造一组悬浮物浓度序列xi(i=1,2,…,n);
A2:构造一秩序列aij,表示xi>xj(1≤j≤i)的样本累积数,如下:
<mrow>
<msub>
<mi>a</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>=</mo>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mn>1</mn>
</mtd>
<mtd>
<mrow>
<msub>
<mi>x</mi>
<mi>i</mi>
</msub>
<mo>></mo>
<msub>
<mi>x</mi>
<mi>j</mi>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mrow>
<msub>
<mi>x</mi>
<mi>i</mi>
</msub>
<mo>&le;</mo>
<msub>
<mi>x</mi>
<mi>j</mi>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>,</mo>
<mn>1</mn>
<mo>&le;</mo>
<mi>j</mi>
<mo>&le;</mo>
<mi>i</mi>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</mrow>
定义Sk:
其中k=2,3,…,n (3)
Sk均值E(Sk)以及方差var(Sk)定义如下:
<mrow>
<mi>E</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>S</mi>
<mi>k</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mrow>
<mi>k</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
<mn>4</mn>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>4</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mi>var</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>S</mi>
<mi>k</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mrow>
<mi>k</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mi>k</mi>
<mo>+</mo>
<mn>5</mn>
<mo>)</mo>
</mrow>
</mrow>
<mn>72</mn>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>5</mn>
<mo>)</mo>
</mrow>
</mrow>
定义统计量:
其中k=2,3,…,n (6)
其中UF1=0;
A3:将序列xi逆序排列,即xn,xn-1,…,x1,重复步骤A2,得出统计量UFk *,令UBk=-UFk *,给出显著水平α,取α=0.05,查正态分布表,得临界值u0.05=±1.96。
A4:将UFk和UBk两个统计量序列曲线和信度线(±1.96两条直线,上限称为LU,下限称为LD)绘在同一张图上,分析绘出的UFk和UBk曲线图,当UFk大于0时,表明序列呈上升趋势;若大于1.96,则表明序列上升趋势显著;当UFk小于0时,表明序列呈下降趋势;若小于-1.96,则表明序列下降趋势显著;当UFk曲线与UBk曲线出现交点时,且交点处于±1.96之间,则表明存在突变点,且交点即为突变点;
将UFk曲线与UBk曲线交点所对应的切应力定义为τ1,为底泥由未起悬向小规模起悬状态过渡的临界切应力,UFk曲线与LU直线交点所对应的切应力定义为τ2,为底泥由小规模起悬状态向大规模起悬状态过渡的临界切应力,并结合拟合曲线结果分析临界切应力取值的合理性。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710865987.7A CN107621532A (zh) | 2017-09-22 | 2017-09-22 | 基于突变分析确定浅水湖泊底泥起悬临界切应力的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710865987.7A CN107621532A (zh) | 2017-09-22 | 2017-09-22 | 基于突变分析确定浅水湖泊底泥起悬临界切应力的方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN107621532A true CN107621532A (zh) | 2018-01-23 |
Family
ID=61090194
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710865987.7A Pending CN107621532A (zh) | 2017-09-22 | 2017-09-22 | 基于突变分析确定浅水湖泊底泥起悬临界切应力的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107621532A (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108827871A (zh) * | 2018-08-17 | 2018-11-16 | 河海大学 | 一种管式泥沙侵蚀试验装置中泥沙表面切应力确定方法 |
CN108984771A (zh) * | 2018-07-23 | 2018-12-11 | 河海大学 | 基于河道断面坡度值的Mann-Kendall突变检验的河道宽度提取方法 |
CN110162905A (zh) * | 2019-05-29 | 2019-08-23 | 河海大学 | 一种基于内源释放的大型浅水湖泊模型及其构建方法 |
CN112464169A (zh) * | 2020-12-03 | 2021-03-09 | 上海卫星工程研究所 | 基于统计分布的大数据单参数阈值挖掘与修正方法及*** |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101789055A (zh) * | 2010-03-23 | 2010-07-28 | 北京师范大学 | 一种人为干扰河流生态需水阈值计算方法 |
CN101944161A (zh) * | 2010-09-03 | 2011-01-12 | 北京师范大学 | 一种基于对生态***扰动程度评价的湿地生态需水量计算方法 |
CN105046051A (zh) * | 2015-06-19 | 2015-11-11 | 长安大学 | 一种实际环境冻融作用次数预估方法 |
CN105354415A (zh) * | 2015-10-26 | 2016-02-24 | 江西省水利科学研究院 | 一种湖泊生态***稳态转换突变的分析方法 |
CN106205060A (zh) * | 2016-08-19 | 2016-12-07 | 临沂大学 | 用于露天矿坑尾矿库边坡滑坡预警预报方法 |
-
2017
- 2017-09-22 CN CN201710865987.7A patent/CN107621532A/zh active Pending
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101789055A (zh) * | 2010-03-23 | 2010-07-28 | 北京师范大学 | 一种人为干扰河流生态需水阈值计算方法 |
CN101944161A (zh) * | 2010-09-03 | 2011-01-12 | 北京师范大学 | 一种基于对生态***扰动程度评价的湿地生态需水量计算方法 |
CN105046051A (zh) * | 2015-06-19 | 2015-11-11 | 长安大学 | 一种实际环境冻融作用次数预估方法 |
CN105354415A (zh) * | 2015-10-26 | 2016-02-24 | 江西省水利科学研究院 | 一种湖泊生态***稳态转换突变的分析方法 |
CN106205060A (zh) * | 2016-08-19 | 2016-12-07 | 临沂大学 | 用于露天矿坑尾矿库边坡滑坡预警预报方法 |
Non-Patent Citations (2)
Title |
---|
王胜艳等: "风浪模拟扰动下梅梁湾沉积物的实验研究", 《环境科学与技术》 * |
鲁春霞等: "安康水库下游径流演变及水库运行对径流影响分析", 《长江流域资源与环境》 * |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108984771A (zh) * | 2018-07-23 | 2018-12-11 | 河海大学 | 基于河道断面坡度值的Mann-Kendall突变检验的河道宽度提取方法 |
CN108827871A (zh) * | 2018-08-17 | 2018-11-16 | 河海大学 | 一种管式泥沙侵蚀试验装置中泥沙表面切应力确定方法 |
CN108827871B (zh) * | 2018-08-17 | 2020-11-10 | 河海大学 | 一种管式泥沙侵蚀试验装置中泥沙表面切应力确定方法 |
CN110162905A (zh) * | 2019-05-29 | 2019-08-23 | 河海大学 | 一种基于内源释放的大型浅水湖泊模型及其构建方法 |
CN110162905B (zh) * | 2019-05-29 | 2022-09-16 | 河海大学 | 一种基于内源释放的大型浅水湖泊模型及其构建方法 |
CN112464169A (zh) * | 2020-12-03 | 2021-03-09 | 上海卫星工程研究所 | 基于统计分布的大数据单参数阈值挖掘与修正方法及*** |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107621532A (zh) | 基于突变分析确定浅水湖泊底泥起悬临界切应力的方法 | |
CN108108911B (zh) | 基于水生态***健康的大型湖泊分区水质目标制定方法 | |
Ge et al. | Analysis of the impact on ecosystem and environment of marine reclamation--A case study in Jiaozhou Bay | |
CN102937637A (zh) | 水体实验装置及方法 | |
CN110580327A (zh) | 一种河流生态流量计算方法 | |
CN110162905B (zh) | 一种基于内源释放的大型浅水湖泊模型及其构建方法 | |
Hu | A review of the models for Lake Taihu and their application in lake environmental management | |
CN104914227B (zh) | 基于多高斯核自优化相关向量机的污水水质软测量方法 | |
CN111461503B (zh) | 基于微生物p/r值评估河流减水河段环境流量的方法 | |
CN103471991A (zh) | 一种基于菌胶团显微图像分析的活性污泥特性确定方法 | |
Voleti | Experimental studies of vertical mixing in an open channel raceway for algae biofuel production | |
CN101440613A (zh) | 悬沙模型试验分布加沙***及方法 | |
CN106404639A (zh) | 一种基于计数框的浮游藻类定量方法 | |
Garcia et al. | A model for predicting the temporal evolution of dissolved oxygen concentration in shallow estuaries | |
CN107818220B (zh) | 基于生态***动力学综合模型对海湾环境容量的估算方法 | |
CN110031592A (zh) | 一种快速测定植物抗盐能力的方法 | |
CN112183123B (zh) | 用于水华防控的河道内生态环境需水量分析方法 | |
Mukherjee et al. | An ex situ and in vitro approach to delineate pennate diatom species with bioindicator potentials in a well mixed tropical estuarine ecosystem | |
Bazgirkhoob et al. | A Numerical Tool for Dissolved Oxygen Simulation in the Western Mississippi Sound | |
CN110208458B (zh) | 用于研究水动力学条件下有机物降解系数变化的试验装置及方法 | |
CN109580514B (zh) | 一种测量非生物颗粒态磷和生物颗粒态磷含量的方法 | |
CN211856261U (zh) | 温度和紊流同时可控的泥沙沉降速率观测实验装置 | |
CN102332110A (zh) | 一种“基于gis网格”的沼泽化程度评价技术 | |
CN111475910B (zh) | 基于分形理论的河口盐淡水混合程度的检测与统计方法 | |
CN213580015U (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 | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20180123 |
|
RJ01 | Rejection of invention patent application after publication |