基于生态输水调整流域小气候和生态系统的方法

本发明涉及生态修复技术领域，具体涉及基于生态输水调整流域小气候和生态系统的方法。

背景技术：

塔里木河是我国最长的内陆河，从最长的源流叶尔羌河算起，全长2486公里，是中国乃至世界生态最为脆弱的流域之一。塔里木河流域的荒漠河岸胡杨林及各种树木和植物对稳定河道、维持流域生态平衡、防风固沙、调节绿洲气候和形成肥沃的森林土壤起到至关重要的作用，是保证绿洲发展的天然生态屏障。

但随着近年来经济社会对陆生生态系统的扰动加剧，导致流域胡杨林的景观格局及群落结构和功能发展发生一定的变化，各源流汇入的水量急剧减少，甚至断流，使得塔里木河流域的胡杨林及各种树木和植物的生态结构发生了萎缩，生态结构受到了严重的破坏，生态系统发生变化的同时使得周边的小气候(温度、湿度、光照、通风)也跟恶化。

技术实现要素：

针对现有技术中的上述不足，本发明提供的基于生态输水调整流域小气候和生态系统的方法能够通过生态输水调整研究流域的小气候和生态系统。

为了达到上述发明目的，本发明采用的技术方案为：

提供一种基于生态输水调整流域小气候和生态系统的方法，其包括：

s1、获取研究流域的植物品种及其覆盖面积，并计算所有植物品种的总覆盖面积；

s2、当总覆盖面积小于预设值时，计算研究流域上一年度的标准化降水蒸散指数，并判断标准化降水蒸散指数是否位于表征干旱的数值区间；若是，进入步骤s3，否则，进入步骤s6；

s3、获取所有植物处于最佳生长状态时的单位土壤截留量，采用最大单位土壤截留量和下一年度规划的植被覆盖总面积，计算土壤的总截流量；

s4、根据每种植物的覆盖面积，计算植物一年的生长需水量，计算采用下一年度规划新增面积种植草本植物时的生长需水量，累加所有生长需水量得到总生长需水量；

s5、累加总截流量和总生长需水量作为补给的包气带水量，并根据包气带水量及其在生态输水中的统计占比，计算得到研究流域下一年度的生态补水量；

s6、根据最大单位土壤截留量和下一年度规划新增面积，计算预计总截流量和种植草本植物时的预计生长需水量；

s7、累加预计总截流量与预计生长需水量作为包气带补水量，采用包气带补水量和统计占比，计算得到研究流域下一年度需新增的生态补水量。

进一步地，当总覆盖面积大于等于于预设值时，执行如下步骤：

s8、计算研究流域上一年度的土壤干旱指数，并判断土壤干旱指数是否大于预设干旱指数；若是，进入步骤s9，否则，按照上一年度的生态补水量进行补给；

s9、采用最大单位土壤截留量计算的土壤截留量及所有植物生长需水量及统计占比，计算得到研究流域下一年度的生态补水量。

本发明的有益效果为：本方案通过研究流域植被的总覆盖面积进行判断，首先可以确定研究流域的生态系统是否处于良性发展；在非良性发展的情况下，通过标准化降水蒸散指数分析植被总覆盖面积不足时是否干旱造成的，还是植被规划面积不足造成的，然后根据具体情况进行生态补水量，以对研究流域的生态系统进行改善，以进一步改善研究流域的小气候。

本方案结合土壤干旱指数可以在总覆盖面积达到规划面积时，对目前的生态补水量是否能够满足植被良性生长进行确认，以准确对生态补水量进行调整，以避免流域的生态系统发生退化，影响流域的生态平衡。

附图说明

图1为基于生态输水调整流域小气候和生态系统的方法的流程图。

具体实施方式

下面对本发明的具体实施方式进行描述，以便于本技术领域的技术人员理解本发明，但应该清楚，本发明不限于具体实施方式的范围，对本技术领域的普通技术人员来讲，只要各种变化在所附的权利要求限定和确定的本发明的精神和范围内，这些变化是显而易见的，一切利用本发明构思的发明创造均在保护之列。

参考图1，图1示出了基于生态输水调整流域小气候和生态系统的方法的流程图；如图1所示，该方法包括步骤s1至步骤s9。

在步骤s1中，获取研究流域的植物品种及其覆盖面积，并计算所有植物品种的总覆盖面积；

在步骤s2中，当总覆盖面积小于预设值时，计算研究流域上一年度的标准化降水蒸散指数，并判断标准化降水蒸散指数是否位于表征干旱的数值区间；若是，进入步骤s3，否则，进入步骤s6；

在本发明的一个实施例中，所述标准化降水蒸散指数的计算方法包括：

s21、采用thornthwaite方法计算逐月的潜在蒸散量：

其中，peti为第i个月潜在蒸散量；k为根据维度计算的修正系数；ti为第i个月的平均气温；i为年总加热指数；m为由i决定的系数；

s22、根据潜在蒸散量和月降水量，计算逐月降水与潜在蒸散的差值：

di＝pi-peti

其中，pi为第i个月降水量；di为第i个月的降水量与潜在蒸散量的差值；

s23、采用三参数的log-logistic分布对di进行拟合，并求出累计函数：

其中，α为尺度参数，β为形状参数，γ为origin参数，f(x)为概率密度函数，f(x)为概率分布函数；x为自变量；

s24、对概率分布函数进行标准化正态处理，得到标准化降水蒸散指数：

其中，spei为上一年度的标准化降水蒸散指数；p为与概率分布函数相关的参数，当p≤0.5时，p＝f(x)，当p>0.5时，p＝1-f(x)；w为中间参数；常数c0、c1、c2、d1、d2分别为c0＝2.515517，c1＝0.802853，c2＝0.010328，d1＝1.432788，d2＝0.001308；

s25、spei指数具有多时间尺度的特征，选取spei指数中12个月尺度的spei指数作为上一年度的标准化降水蒸散指数。

在步骤s3中，获取所有植物处于最佳生长状态时的单位土壤截留量，采用最大单位土壤截留量和下一年度规划的植被覆盖总面积，计算土壤的总截流量；其中，土壤截留量可以采用tdr土壤含水量监测器进行检测，之后再折算成单位面积的单位土壤截留量。

在步骤s4中，根据每种植物的覆盖面积，计算植物一年的生长需水量，计算采用下一年度规划新增面积种植草本植物时的生长需水量，累加所有生长需水量得到总生长需水量。

实施时，本法干优选根据植物的覆盖面积，计算植物一年的生长需水量的方法包括：

s41、在室内构建一块试验区域，将试验区域划分为若干试验田，并在每块试验田中分别种植上研究流域中存在的植物；

s42、调整试验区域的环境参数，使试验田中种植的植物处于最佳的生长状态，并根据试验区域的环境参数，采用彭曼公式计算每种植物的逐旬蒸散量：

其中，et0为植物的逐旬蒸散量；δ为饱和水汽压-温度曲线斜率；rn为植物冠层的净辐射；g为增热土壤消耗的能量；γ为湿度计常数；t为平均气温；u2为2m高风速；ea为饱和水汽压；ed为实际观测水汽压；

s43、累加同一植物在一年内的所有逐旬蒸散量，得到同一植物在试验田中的年度蒸散量，并将年度蒸散量折算为单位面积内的年度蒸散量；

s44、称取每块试验田中植物在年初和年末的质量，得到两者的质量差，并将其折算成单位面积内的年度植物储水量；

s45、根据研究流域每种植物的覆盖面积及单位面积内的年度蒸散量和年度植物储水量，计算得到植物一年的生长需水量。

在步骤s5中，累加总截流量和总生长需水量作为补给的包气带水量，并根据包气带水量及其在生态输水中的统计占比，计算得到研究流域下一年度的生态补水量；

在步骤s6中，根据最大单位土壤截留量和下一年度规划新增面积，计算预计总截流量和种植草本植物时的预计生长需水量。

一般情况下，在流域中若是供水充足，首先生长的是一年生草本植物，若是连续几年存在供水较充足时，才会有较大的灌木生长，为了保证生态补水量的合理使用，本方案规划在新增面积上种植草本植物，先通过草本植物对水土进行改善后，再慢慢栽种灌木。

在步骤s7中，累加预计总截流量与预计生长需水量作为包气带补水量，采用包气带补水量和统计占比，计算得到研究流域下一年度需新增的生态补水量。

实施时，本方案优选所述统计占比的计算方法包括：

s71、获取研究流域临近下一年度的若干次输水过程中的流域来水量、流域出水量和包气带补给量。

河水在流动的过程中通过蒸发和渗流转化为大汽水、土壤水、地下水。若以河道流域为研究对象，可将流量的沿程损失称之为河道的耗水过程。河道耗水量取决于流量、河床岩性、河道形状及河道附近的植被状况等自然属性。通过对河道上各监测断面的数据分析，利用水量平衡的关系，则可以得出不同河段的耗水量特征，用以下关系式表达：

δw＝w1-w2＝w包+w蒸+w地，dw＝δw/l，k＝dw/w1

其中，w1为河道流域上断面的流域来水量；w2为河道流域下断面的流域流出量；δw为断面间的区间耗水量；w包为包气带水量；w蒸为河道流域的蒸发量；w地为补入地下水的水量；l为两断面间的河道长度；dw为单位河长耗水量；k为为单位河长耗水率。

s72、计算每次输水时，包气带补给量相对流域来水量和流域出水量间差值(相当于断面间的区间耗水量)的占比；并采用所有的占比的平均值作为统计占比。

当总覆盖面积大于等于于预设值时，执行如下步骤：

在步骤s8中，计算研究流域上一年度的土壤干旱指数，并判断土壤干旱指数是否大于预设干旱指数；若是，进入步骤s9，否则，按照上一年度的生态补水量进行补给；

其中，土壤干旱指数的计算方法包括：

s81、获取植被归一化指数ndvi、湿边拟合方程系数a1和b1及干边拟合方程系数a2和b2，计算陆面温度极低值tsmin和陆面温度极高值tsmax：

tsmin＝a1+b1ndvi

tsmax＝a2+b2ndvi

其中，ρ为反射率修正系数；

s82、根据陆面温度极低值tsmin和陆面温度极高值tsmax，计算土壤干旱指数：

其中，tvdi为土壤干旱指数，其值介于0～1之间，tvdi越接近0，表明土壤水分越接近土壤持水量，越接近1，表明土壤水分含量越接近萎蔫点；ts表示陆面温度。

在步骤s9中，采用最大单位土壤截留量计算的土壤截留量及所有植物生长需水量及统计占比，计算得到研究流域下一年度的生态补水量。

步骤s9可以进一细化为：根据最大单位土壤截留量和总覆盖面积，分别计算土壤截留量和植物生长需水量，之后累加土壤截留量和植物生长需水量作为包气带水量，采用包气带水量及统计占比，计算得到研究流域下一年度的生态补水量(生态补水量＝包气带水量/统计占比)。

本方案的流域出水量通过在研究流域与下游交汇面修建大坝进行调控，通过大坝可以对新增的生态补水量进行拦截存储，以保证新增的生态补水量慢慢浸润或者漫溢至研究流域的生态系统处的土地中。

综上所述，本方案提供的基于生态输水调整流域小气候和生态系统的方法能够通过标准化降水蒸散指数和土壤干旱指数对研究流域的生态补水量进一调整，以使研究流域的的生态系统和小气候得以调整。