黄土丘陵区沟道土壤水分估算方法,申请号CN201310165487.4-传众专利搜索

黄土丘陵区沟道土壤水分估算方法

申请公布号：CN103246809A

申请号：CN201310165487.4

申请日期：2013.05.07

申请公布日期：2013.08.14

申请人：

西北农林科技大学

发明人：吴普特;赵西宁;高晓东;史银光;王荚文

分类号：G06F19/00(2011.01)I

主分类号：G06F19/00(2011.01)I

代理机构：

上海汉声知识产权代理有限公司 31236

代理人：牛山;郭国中

地址：712100 陕西省西安市杨凌示范区西农路22号

摘要：本发明提供了一种黄土丘陵区沟道土壤水分估算方法，包括：合理布设沟道和坡面土壤水分监测样点，覆盖研究区不同土地类型和空间位置；计算沟道和坡面土壤水分每次采样的空间算术均值，获取所需土壤水分时间序列数据；分别对沟道和坡面土壤水分时间序列数据进行秩排序，并计算对应位置土壤水分差值；然后采用多项式拟合所计算差值与坡面秩序列土壤水分关系，并利用最小二乘法确定多项式最佳阶数，进而得到观测算子；最后基于观测算子和坡面秩序列土壤水分数据，得到沟道秩序列土壤水分，然后将其转换成时间序列数据，并计算误差。本发明只需通过坡面土壤水分数据便可获得沟道对应土层土壤水分，无需再进行沟道土壤水分采样，估算精度高。

主权项：1.一种黄土丘陵区沟道土壤水分估算方法，其特征在于，包括如下步骤：步骤1：在研究区沟道和坡面分别布设土壤水分监测样点，其中，样点覆盖沟道和坡面不同土地类型和地形位置，以便能够代表沟道和坡面土壤水分实际状况；步骤2：通过连续测定方法，获取小流域沟道和坡面的土壤水分数据，连续测定方法中的测定时间涵盖研究区土壤干湿状况，以便能够代表所述研究区生长季土壤水分变化状况，其中，所述土壤水分数据包括样点土壤含水量；步骤3：假设沟道和坡面各自所有样点土壤含水量的算术平均值分别代表沟道和坡面的土壤平均含水量，在剔除错误数值的基础上，将沟道所有样点土壤含水量算术平均值的时间序列作为实测沟道土壤水分时间序列数据ηj，将坡面所有样点土壤含水量算术平均值的时间序列作为实测坡面土壤水分时间序列数据μj；<maths num="0001"><![CDATA[<math><mrow><msub><mi>η</mi><mi>j</mi></msub><mo>=</mo><munderover><mi>Σ</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><msub><mi>η</mi><mi>ij</mi></msub><mo>;</mo><msub><mi>μ</mi><mi>j</mi></msub><mo>=</mo><munderover><mi>Σ</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></munderover><msub><mi>μ</mi><mi>ij</mi></msub></mrow></math>]]></maths>式中ηij和μij分别为沟道和坡面第i个样点第j次测得的土壤含水量；N和M分别是沟道和坡面所布设土壤水分监测样点数目；步骤4：分别将实测沟道土壤水分时间序列数据和实测坡面土壤水分时间序列数据进行秩排序分别得到实测沟道土壤水分秩序列数据ηk和实测坡面土壤水分秩序列数据μk，进而构建估算沟道土壤水分所需的观测算子<img file="FDA00003153131400012.GIF" wi="89" he="94" />步骤5：通过观测算子<img file="FDA00003153131400013.GIF" wi="56" he="98" />和实测坡面土壤水分秩序列数据μk得到估算沟道土壤水分秩序列数据<img file="FDA00003153131400014.GIF" wi="90" he="83" /><maths num="0002"><![CDATA[<math><mrow><msub><mover><mi>η</mi><mo>~</mo></mover><mi>k</mi></msub><mo>=</mo><msub><mi>μ</mi><mi>k</mi></msub><mo>+</mo><msub><mover><mi>δ</mi><mo>~</mo></mover><mi>k</mi></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow></math>]]></maths>式（1）中，<img file="FDA00003153131400016.GIF" wi="70" he="85" />为估算沟道土壤水分秩序列数据，μk为实测坡面土壤水分秩序列数据，<img file="FDA00003153131400017.GIF" wi="64" he="101" />为观测算子；然后将估算沟道土壤水分秩序列数据<img file="FDA00003153131400018.GIF" wi="60" he="85" />和实测坡面土壤水分秩序列数据μk分别转换为估算沟道土壤水分时间序列数据和实测坡面土壤水分时间序列数据，并计算误差；最后采用不同水文年土壤水分时间序列数据对观测算子进行验证，若验证误差在允许范围内，则表明所建观测算子能够适用于该研究区不同水文年；若验证误差在允许范围外，则需构建适用于不同水文年的观测算子；其中，所述验证是指，首先基于验证期内不同水文年的实测坡面土壤水分时间序列数据μj和所构建观测算子<img file="FDA00003153131400021.GIF" wi="92" he="101" />通过式（1）计算得到估算沟道土壤水分秩序列数据<img file="FDA00003153131400022.GIF" wi="96" he="82" />再转换为估算沟道土壤水分时间序列数据，然后将估算沟道土壤水分时间序列数据与验证期实测沟道土壤水分时间序列数据进行对比分析，计算误差。