[发明专利]一种区域旱灾风险估计方法

权利要求书

1.一种区域旱灾风险估计方法，其特征在于，该方法包括：

从不同数据源获取感兴趣区域内的针对至少一个感兴趣点的至少一个孕灾因子、至少一个致灾因子、至少一个承灾体易损性评价因子和旱灾区划指数；

基于所述至少一个孕灾因子计算旱灾敏感性指数；

基于所述至少一个致灾因子计算旱灾危险性指数；

基于所述至少一个承灾体易损性评价因子计算旱灾脆弱性指数；

对所述旱灾区划指数、所述旱灾敏感性指数、所述旱灾危险性指数和所述旱灾脆弱性指数进行加权求和以计算旱灾风险指数；以及

根据计算出的旱灾风险指数来估计旱灾风险；

其中，所述方法还包括：

将旱灾区划指数根据分级标准分级，并对每个级赋以相应的值；

所述致灾因子包括：降水量距平百分率指数、土壤相对湿度指数、土壤热惯量指数和距平植被指数，

其中，降水量距平百分率指数为：

$P_l=\frac{P_r-\stackrel{\‾}{P}}{\stackrel{\‾}{P}}\×100\%,$

P_l是降水量距平百分率指数，P_r为最近时段降水量，为时段n年平均水量；

土壤相对湿度指数为：

$Q_i=\frac{Q_t-\stackrel{\‾}{Q_t}}{\stackrel{\‾}{Q_t}}\×100\%$

其中，Q_i为土壤相对湿度指数，Q_t为评价时段实测土壤相对湿度，为评价时段最适土壤相对湿度；

距平植被指数为：

$AVI={\mathrm{NDVI}}_i-\stackrel{\‾}{NDVI}$

AVI为距平植被指数，NDVI_i为某一年某一时期NDVI的值，为n年该时期NDVI的平均值，NDVI为归一化植被指数；

土壤热惯量指数按如下方式获得：

P＝B(1-A)/ΔTs

P为土壤热惯量，ΔTs为地表日较差，A为地表全波段反照率，A＝0.423ρ₁+0.577ρ₂，ρ₁，ρ₂分别是第一、第二波段的波谱反射率，B为经验系数；

根据土壤热惯量与土壤水分含量和土壤密度间如下关系建立土壤热惯量与土壤水分含量和土壤密度的查找表；

$P={\{2.1*{\mathrm{ds}}^{\[1.2-0.02*(ds/d)*w\]}*e^{\[-0.007*{(w*ds/d-20)}^2\]}+{\mathrm{ds}}^{\[0.8+0.02*(ds/d)*w\]}\}}^{1/2}*(0.2*w/d)*{\mathrm{ds}}^2/0.01$

其中，P为土壤热惯量，ds为土壤密度，d为水的密度，w为土壤水分含量的质量百分含量；

根据查找表确定土壤热惯量所对应的土壤水分含量，以土壤水分含量作为土壤热惯量指数。