申请/专利权人：中国水利水电科学研究院

申请日：2023-09-08

公开（公告）日：2024-04-02

公开（公告）号：CN117454062B

主分类号：G06F17/11

分类号：G06F17/11;G06F18/10;G06F17/18;G06F18/24

优先权：

专利状态码：有效-授权

法律状态：2024.04.02#授权;2024.02.13#实质审查的生效;2024.01.26#公开

摘要：本发明公开了一种植被恢复潜力的计算方法，所述方法包括以下步骤：步骤1、收集数据：收集研究区的气象数据、植被数据以及土壤数据；步骤2、数据预处理：对收集的数据进行设置或计算得到降水参数、植被参数、土壤参数以及潜在蒸散发参数；步骤3、将数据匹配到空间数据计算网格；步骤4、植被恢复潜力计算。本发明所述方法基于植被生长对水分和光热条件的需求，结合植被生长过程，能够科学计算植被恢复潜力，机理明确、精确度高、可靠性强、适用性广，能够在不同地区推广应用，解决了区域植被修复合理规划的问题。在海河流域山区实际应用中，精确度能够达到87％，大大提高了计算结果的准确性。

主权项：1.一种植被恢复潜力的计算方法，其特征在于，所述方法包括以下步骤：步骤1、收集数据：收集研究区的气象数据、植被数据以及土壤数据；所述气象数据包括研究区气象站点的逐日降水数据、最低气温、最高气温、平均气温、风速、相对湿度、日照时数、大气压；所述植被数据包括植被类型数据；所述土壤数据包括土壤透水系数、土壤含水量、土壤热通量；步骤2、数据预处理：对收集的气象数据、植被数据以及土壤数据进行设置或计算得到降水参数、植被参数、土壤参数以及潜在蒸散发参数；步骤3、将数据匹配到空间数据计算网格：首先划分空间数据计算网格，空间分辨率为5km，将气象数据利用插值法插补到预先设置好的格点上；统计各网格内各植被类型的比例；将土壤数据映射到对应的网格上；数据序列保证时间和空间上的完整性，对缺失的数据采用临近数据进行插补延长；步骤4、植被恢复潜力计算：建立水量平衡方程，计算区域水分和光热条件支撑的植被最优覆盖度，进而计算植被恢复潜力；步骤2所述降水参数根据气象站点的逐日降水数据计算得到，所述植被参数根据植被类型数据设置得到，所述土壤参数根据土壤数据计算得到，所述潜在蒸散发参数包括非生长季饱和裸土潜在蒸散发参数与生长季饱和裸土潜在蒸散发参数，所述生长季为4-10月；步骤4所述植被恢复潜力计算具体包括以下步骤：步骤41、建立水量平衡方程：建立生长季系列的水量平衡方程为： 式中，I为植被冠层蒸腾量；II为冠层降水截留量；III为土壤水蓄变量；VI为裸土蒸发量；V为地表径流量；VI为深层渗漏量；VII为毛细管上升量；mh为生长季平均单次降水深，mm；I中，M为生长季植被覆盖度，无量纲；kv为植被实际冠层导度，无量纲；βv为裸土蒸发效率系数，无量纲；m′tb为平均降水间隙中用于蒸腾的平均时间，day；Eps为生长季饱和裸土潜在蒸散发参数，mm·day-1；II中，为冠层降水截留量，mm； 包括植被截留量和裸土截留量，由下式计算得到: 式中，h0为平均截留深度，mm；η0为气孔面积与叶面积之比，无量纲；β为叶倾角的余弦值，无量纲；Lt为叶面积指数，无量纲；III中，mv为生长季降水次数；ΔS为土壤需水变化量，mm；所述土壤需水变化量的计算公式为：ΔS＝Pd-1-MdEpdmd-Yd3式中，Pd为非生长季降水，mm；Md为非生长季植被盖度，day；Epd为非生长季饱和裸土潜在蒸散发参数，mm·day-1；md为非生长季时长，day；Yd为非生长季径流，mm；VI中，βs为裸土蒸发系数；m″tb为平均降水间隙中用于裸土蒸发的时间，day；V中，产流方式设为超渗产流，参数G和σ由下式计算得到 κ0＝λ0mh7 式中，ω为平均降水强度mi的倒数，day·mm-1；K1为土壤饱和导水率，mm·day-1；s0为研究时段内平均土壤含水量，无量纲；c为土壤透水系数，无量纲；w为毛细管上升速率，mm·day-1；δ为概率密度参数，无量纲；m为土壤孔隙大小分布指数，无量纲；k0为形状参数，无量纲；ne为有效土壤孔隙度，无量纲；λ0为降水泊松分布的尺度参数，无量纲；为mh的方差；ψ1为土壤饱和基质势，mm；φc,s0为土壤吸附扩散度，无量纲；VI中，Pt为生长季平均降水，mm；mt为生长期日数，day；VII中，zw为地下水埋深，mm；基于水量平衡方程，生长季植被覆盖度和植被潜在冠层导度由以下公式计算获得： 式中，Ve为可供植被利用的水分，mm；kv*为植被潜在冠层导度，无量纲；步骤42、计算植被实际冠层导度：所述生长季饱和裸土潜在蒸散发参数的计算公式为： 式中，λ为蒸发潜热，Kcal·g-1；Rn为地表净辐射量，MJ·m-2·day-1；ρ为空气质量密度，kg·m-3；cp为大气的定压比热，MJ·kg-1·℃-1；es和ea分别为饱和水汽压和实际水汽压，kPa；ra为空气动力阻抗，s·m-1；Δ为温度-饱和水汽压曲线斜率，kPa·℃-1；γ0为干湿表常数，kPa·℃-1；生长季植被潜在蒸腾的计算公式为： 式中，Ev为植被潜在蒸腾，mm·day-1；rc为冠层阻抗，s·m-1；rcra为阻抗比；计算植被实际冠层导度，计算公式为： 式中，kv为植被实际冠层导度，无量纲；植被冠层的阻抗比通过M＝0和M＝1对应的阻抗比之间进行线性插值得到： 式中，和分别为M＝0和M＝1对应的阻抗比；对于M→0的稀疏冠层，根据以下公式计算： 式中，hsh为树干分数，即树干的高度和植被单株高度之比；对于M→1的闭合冠层，植被冠层阻抗比计算公式为： 步骤43、计算区域水分和光热条件支撑的植被最优覆盖度：联立水量平衡方程和植被潜在冠层导度，将kv*～M曲线绘制在同一坐标系里，得到最优植被覆盖度的交点；当0M1时，曲线相交，即： 两条曲线的交点即为最优植被覆盖度M*；若M达到1时，曲线仍不相交，即： 则最优植被覆盖度M*＝1；步骤44、计算植被恢复潜力：对植被恢复潜力进行计算，植被恢复潜力即为在研究区区域水分和光热条件支撑的植被最优覆盖度与当前现状植被覆盖度对比，其差值即为植被恢复潜力。

全文数据：

权利要求：

百度查询： 中国水利水电科学研究院 一种植被恢复潜力的计算方法

免责声明
1、本报告根据公开、合法渠道获得相关数据和信息，力求客观、公正，但并不保证数据的最终完整性和准确性。
2、报告中的分析和结论仅反映本公司于发布本报告当日的职业理解，仅供参考使用，不能作为本公司承担任何法律责任的依据或者凭证。

阅读全文 双屏查看 官方信息 专利公告 收藏专利 下载PDF 下载WORD