发明名称 一种土壤含水率间接获取及快速评价方法
摘要 本发明公开了一种土壤含水率间接获取及快速评价方法,该方法依次包括如下步骤:1、基于研究区域,设计取样方案;2、使用环刀法对每个取样点取样,测量土壤容重,土壤含水率和土壤紧实度;3、使用模型对数据进行拟核,获得不同模型的待确定系数和相关系数,确定本研究区域的应用模型;4、对每个取样点进行土壤紧实度测试,将所得到的数据代到应用模型中,得到土壤含水率的结果;5、使用面积加权平均值的方法,对土壤含水率进行整体计算,以快速评价研究区域的墒情。该发明的实施方案系统条理,获取速度快,样本数量多,可有效的降低土壤空间异质性导致的不确定性,可应用于不同作物覆盖,土壤质地和地形地貌的地区。
申请公布号 CN104331639A 申请公布日期 2015.02.04
申请号 CN201410719050.5 申请日期 2014.12.01
申请人 中国矿业大学 发明人 杨德军;卞正富;雷少刚;熊集兵
分类号 G06F19/00(2011.01)I;G01N33/24(2006.01)I 主分类号 G06F19/00(2011.01)I
代理机构 南京瑞弘专利商标事务所(普通合伙) 32249 代理人 杨晓玲
主权项 一种土壤含水率间接获取及快速评价方法,其特征在于包括如下步骤:步骤1、基于研究区域,综合考虑土壤质地、作物类型和地形地貌,采取等距取样法,设计取样方案,布置取样点;每个取样点的面积S<sub>i</sub>(i=1,2,...N)为3~5m<sup>2</sup>,取样数N为15~25个;步骤2、使用环刀法对每个取样点取样,使用烘干法对所取样品分别测量土壤容重BD<sub>i</sub>(i=1,2,...N)和土壤含水率θ<sub>i</sub>(i=1,2,...N);使用土壤紧实度仪对每个取样点土层的土壤紧实度CI<sub>i</sub>(i=1,2,...N)多次测试取平均值;步骤3、基于取样点的面积S<sub>i</sub>、土壤容重BD<sub>i</sub>、土壤含水率θ<sub>i</sub>和土壤紧实度CI<sub>i</sub>(i=1,2,...N)的数据,利用模型:模型1:Ayers模型<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mi>CI</mi><mo>=</mo><mrow><mo>(</mo><msub><mi>A</mi><mn>1</mn></msub><mo>&times;</mo><msup><mi>BD</mi><msub><mi>A</mi><mn>4</mn></msub></msup><mo>)</mo></mrow><mo>/</mo><mrow><mo>(</mo><msub><mi>A</mi><mn>2</mn></msub><mo>+</mo><msup><mrow><mo>(</mo><mi>&theta;</mi><mo>-</mo><msub><mi>A</mi><mn>3</mn></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000621710130000011.GIF" wi="954" he="77" /></maths>模型2:Upadhyaya模型<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mi>CI</mi><mo>=</mo><msub><mi>B</mi><mn>1</mn></msub><msup><mi>BD</mi><msub><mi>B</mi><mn>2</mn></msub></msup><msup><mi>e</mi><mrow><msub><mrow><mo>-</mo><mi>B</mi></mrow><mn>3</mn></msub><mi>&theta;</mi></mrow></msup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000621710130000012.GIF" wi="953" he="74" /></maths>模型3:Busscher模型<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><mi>CI</mi><mo>=</mo><msub><mi>C</mi><mn>1</mn></msub><msup><mi>&theta;</mi><msub><mi>C</mi><mn>2</mn></msub></msup><msup><mi>BD</mi><msub><mi>C</mi><mn>3</mn></msub></msup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000621710130000013.GIF" wi="954" he="74" /></maths>模型4:Hernanz模型<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><mi>CI</mi><mo>=</mo><msub><mi>D</mi><mn>1</mn></msub><msup><mi>&theta;</mi><msub><mi>D</mi><mn>2</mn></msub></msup><msup><mi>BD</mi><msub><mi>D</mi><mn>3</mn></msub></msup><msup><mi>d</mi><msub><mi>D</mi><mn>4</mn></msub></msup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000621710130000014.GIF" wi="957" he="76" /></maths>进行拟核,得到所选取模型的待确定系数和相关系数,相关系数最大的模型确定为本研究区域的应用模型;式中:CI为土壤紧实度,kPa;BD为土壤容重,g/cm<sup>3</sup>;θ为土壤质量含水率;d为深度,mm;A,B,C,D为依赖于土壤类型的拟合常数,且恒为正;步骤4、使用土壤紧实度仪对每个取样点土层的土壤紧实度CI<sub>j</sub>(j=1,2,...M)多次测试取平均值,将土壤紧实度CI<sub>j</sub>代入应用模型,得到土壤含水率θ<sub>j</sub>(j=1,2,...M);步骤5、利用面积加权平均值的方法,通过公式;θ=Σθ<sub>j</sub>×S<sub>j</sub>/ΣS<sub>j</sub>  (5)对研究区域的土壤含水率θ进行整体计算,以此评价研究区域的墒情。
地址 221116 江苏省徐州市大学路1号中国矿业大学科研院