基于联合概率分布的干旱区地下水生态埋深分析方法与流程

1.本发明属于地球物理下生态水文分支技术领域，具体涉及一种基于联合概率分布的干旱区地下水生态埋深分析方法。


背景技术：

2.我国西北干旱区内陆河流域水资源稀缺、生态脆弱。流域降水集中在上游山区，中下游、盆地腹地降水稀少，而蒸发强烈。流域中下游自然植被生存主要依赖于浅层地下水，形成地下水依赖型陆生生态系统。为维系地下水依赖型陆生生态系统的生态服务功能正常发挥，需维持地下水生态埋深，以保障地下水依赖型植被生长的用水需求。在干旱区地下水生态埋深研究中，已有的分析方法大多基于地下水埋深和植被观测资料，其中植被观测资料多采用归一化植被指数(normalized difference vegetation index，ndvi)、植被覆盖度等表征植被生长状况的指标，利用回归分析构建植被指标与地下水埋深的响应关系曲线，进而基于植被指标
‑
地下水埋深响应关系曲线，通过拐点分析、突变点分析等方法确定地下水生态埋深。受土壤质地、养分、盐分等其它生境要素空间异质性的影响，植被指标与地下水埋深的相关关系复杂，植被指标对地下水埋深变化响应的规则性和随机性特征并存。回归拟合分析所得的植被指标
‑
地下水埋深响应关系曲线往往难以全面刻画植被指标与地下水埋深之间的复杂关系，影响地下水生态埋深分析成果的准确性和合理性，不利于干旱区内陆河流域中下游自然绿洲保护及流域可持续发展。为此，考虑采用copula联合概率分布函数构建ndvi与地下水埋深的相关关系，提出了基于联合概率分布的干旱区地下水生态埋深分析方法。


技术实现要素：

3.发明目的：针对以上问题，本发明提出一种基于联合概率分布的干旱区地下水生态埋深分析方法。
4.技术方案：为实现本发明的目的，本发明所采用的技术方案是：一种基于联合概率分布的干旱区地下水生态埋深分析方法，具体包括如下步骤：
5.(1)收集整理研究区的生态水文基础数据资料，包括地下水埋深、ndvi等栅格文件资料；具体地：
6.①
地下水埋深记为h，单位m；植被记为v，采用ndvi表征植被的生长状况。
7.②
地下水埋深栅格文件与ndvi栅格文件在空间上相互匹配对应，每个栅格即为一个植被样点；每个植被样点数据包括一个地下水埋深数据和一个ndvi数据，形成一个数据组对，记为(h
j
,v
j
)，其中j＝1,2,...,n，n为植被样点总数。
8.③
基于地下水埋深资料，以1m为地下水埋深间距，等间距划分地下水埋深，得到若干个地下水埋深区间，划分结果为0
‑
1m,1
‑
2m,2
‑
3m,3
‑
4m,4
‑
5m,5
‑
6m,6
‑
7m等；第i个地下水埋深区间的下边界记为h
i,l
，h
i,l
＝0,1,2,...,m
‑
1，第i个地下水埋深区间的上边界记为h
i,u
，h
i,u
＝1,2,3,...,m，i＝1,2,3,...,m，m为地下水埋深区间总数。
9.(2)对于不同的地下水埋深区间，统计各地下水埋深区间范围内的植被出现频次，拟合植被出现频次—地下水埋深响应关系曲线；具体地：
10.①
分析各个地下水埋深区间的植被出现频次，植被出现频次的统计公式为：
[0011][0012]
式中，p
i
为第i个地下水埋深区间的植被出现频次，n
i
为第i个地下水埋深区间内的植被样点数目；p
i
定量反映了某一地下水埋深区间范围内生长植被的可能性，p
i
数值越大，说明该地下水埋深区间的生境更适合植被生长，在该地下水埋深区间范围内植被生长且长势较好(反映为ndvi数值较大)的可能性越大。
[0013]
②
将地下水埋深区间的植被出现频次p
i
作为地下水埋深区间的中点的植被出现频次，以地下水埋深区间的中点值为自变量，地下水埋深区间的中点值即0.5m,1.5m,2.5m,3.5m,4.5m,5.5m,6.5m等地下水埋深区间的中点的植被出现频次为因变量，采用高斯曲线方程拟合植被出现频次—地下水埋深响应关系函数，记为f
h
：
[0014]
p
h
＝f
h
(h)
[0015]
式中，p
h
为地下水埋深h＝h时的植被出现频次，植被出现频次—地下水埋深响应关系函数f
h
的作用为：计算不同地下水埋深对应的植被出现频次。
[0016]
(3)基于植被出现频次—地下水埋深响应关系函数，计算各个植被样点的地下水埋深所对应的植被出现频次；对植被出现频次进行频率分析，拟合植被出现频次的边缘分布函数；对ndvi进行频率分析，拟合ndvi的边缘分布函数；在边缘分布函数分析的基础上，构造地下水埋深和ndvi的联合概率分布函数；具体地：
[0017]
①
对于每个植被样点的地下水埋深h
j
，采用植被出现频次—地下水埋深响应关系函数f
h
，计算该地下水埋深h
j
所对应的植被出现频次对于所有植被样点的植被出现频次，进行频率分析，采用核分布函数(kernel distribution function)拟合植被出现频次的边缘分布函数，记为
[0018]
②
对于所有植被样点的ndvi v
j
，进行频率分析，采用广义极值分布(generalized extreme value distribution)拟合ndvi的边缘分布函数，记为f
v
。
[0019]
③
对于所有植被样点的地下水埋深和ndvi数据组对(h
j
,v
j
)，采用植被出现频次—地下水埋深响应关系函数f
h
，计算各组对中地下水埋深h
j
所对应的植被出现频次生成植被出现频次和ndvi数据组对对于所有植被样点的植被出现频次和ndvi数据组对进行双变量频率分析，基于植被出现频次的边缘分布函数f
ph
和ndvi的边缘分布函数f
v
，采用gumbel copula函数构建植被出现频次—ndvi联合概率分布函数f
c0
。
[0020]
④
对于不同的地下水埋深，采用植被出现频次—地下水埋深响应关系函数f
h
计算该地下水埋深对应的植被出现频次，通过植被出现频次的边缘分布函数f
ph
计算该地下水埋深对应的概率值；对于不同的ndvi，采用ndvi的边缘分布函数f
v
计算该ndvi对应的概率值；基于不同地下水埋深对应的概率值和不同ndvi对应的概率值，通过植被出现频次—ndvi联合分布函数f
c0
计算联合概率值，生成地下水埋深—ndvi联合概率分布，记为f
c
。
[0021]
(4)设置合理的地下水依赖型植被保护目标；基于地下水埋深—ndvi联合概率分
布，分析不同地下水埋深区间实现植被保护目标的条件概率，以条件概率最大值对应的地下水埋深作为地下水生态埋深；具体地：
[0022]
①
设置合理的ndvi阈值v
t
，作为地下水依赖型植被的保护目标，满足v＞v
t
，v表示不同ndvi值；
[0023]
②
基于地下水埋深—ndvi联合概率分布f
c
，分别计算第i个地下水埋深区间实现植被保护目标的条件概率p
i
(v＞v
t
|h
i,l
<h≤h
i,u
)：
[0024][0025]
式中，采用植被出现频次—地下水埋深响应关系函数f
h
，计算地下水埋深h
i,u
对应的植被出现频次再采用植被出现频次的边缘分布函数计算对应的概率值
[0026]
采用植被出现频次—地下水埋深响应关系函数f
h
，计算地下水埋深h
i,l
对应的植被出现频次再采用植被出现频次的边缘分布函数计算对应的概率值采用地下水埋深—ndvi联合概率分布f
c
，计算f
c
(h
i,u
,v
t
)；
[0027]
采用地下水埋深—ndvi联合概率分布f
c
，计算f
c
(h
i,l
,v
t
)。
[0028]
③
不同地下水埋深所形成的生境条件不同，因此对植被生长的适宜性不同，p
i
(v＞v
t
|h
i,l
<h≤h
i,u
)随地下水埋深变化而变化；当p
i
(v＞v
t
|h
i,l
<h≤h
i,u
)在某一地下水埋深条件下达到最大值时，表明该地下水埋深条件较为适宜植被生长；取p
i
(v＞v
t
|h
i,l
<h≤h
i,u
)最大值对应的地下水埋深，作为地下水生态埋深。
[0029]
有益效果：与现有技术相比，本发明技术方案具有以下有益技术效果：
[0030]
本发明提出一种基于联合概率分布的干旱区地下水生态埋深分析方法，克服现有方法难以全面刻画植被指标与地下水埋深相关关系的不足，提高ndvi、地下水埋深等生态水文观测资料的利用效率，提升地下水生态埋深分析成果的准确性和合理性，完善干旱区内陆河流域生态需水分析方法体系，为保障干旱区内陆河流域地下水依赖型陆生生态系统结构与功能稳定、优化自然生态系统与社会经济系统水资源配置格局提供技术支撑和科学依据，以促进人与自然和谐共生、助力生态文明建设、实现绿色可持续发展。
附图说明
[0031]
图1是一种实施例下本发明的基于联合概率分布的干旱区地下水生态埋深分析方法的流程图；
[0032]
图2是一种实施例下吐鲁番盆地地下水依赖型植被样点位置图；
[0033]
图3是一种实施例下吐鲁番盆地地下水依赖型植被出现频次与地下水埋深关系图；
[0034]
图4是一种实施例下吐鲁番盆地地下水依赖型植被ndvi与地下水埋深联合概率分布图；
[0035]
图5是一种实施例下吐鲁番盆地实现地下水依赖型植被保护目标的条件概率图。
具体实施方式
[0036]
下面结合附图和实施例对本发明的技术方案作进一步的说明。
[0037]
本发明实施例的应用研究区域为吐鲁番盆地。吐鲁番盆地位于我国新疆维吾尔自治区中部，面积约5万km2，三面环山，是一个典型的干旱区内陆河流域。火焰山、盐山沿东西走向贯穿吐鲁番盆地中部，将盆地分为北盆地和南盆地；南盆地有我国内陆最低点—艾丁湖，海拔为
‑
154m。吐鲁番盆地的气候类型为典型的暖温带干旱荒漠气候，盆地地下水主要受周围山区的降水和冰雪融水补给，地下水依赖型植被主要分布在盆地中心腹地。
[0038]
参考图1，为本发明实施例的基于联合概率分布的干旱区地下水生态埋深分析方法的流程图，具体包括如下步骤：
[0039]
第一步：收集整理吐鲁番盆地地下水依赖型植被分布区的基础数据资料，包括地下水埋深栅格资料和ndvi栅格资料，空间分辨率1km
×
1km。以1m为地下水埋深间距，等间距划分地下水埋深区间，即0
‑
1m,1
‑
2m,
……
,29
‑
30m。
[0040]
参考图2，为本发明实施例的吐鲁番盆地地下水依赖型植被样点位置图。
[0041]
第二步：对于不同的地下水埋深区间，统计各地下水埋深区间范围内的植被出现频次，采用高斯曲线(gaussian curve)方程拟合植被出现频次—地下水埋深响应关系曲线。参考图3，为本发明实施例的吐鲁番盆地地下水依赖型植被出现频次与地下水埋深关系图。
[0042]
第三步：基于植被出现频次—地下水埋深响应关系函数，计算各个植被样点的地下水埋深所对应的植被出现频次；对各个植被样点的植被出现频次进行频率分析，采用核分布函数(kernel distribution function)拟合植被出现频次的边缘分布函数；对各个植被样点的ndvi进行频率分析，采用广义极值分布(generalized extreme value distribution)拟合ndvi的边缘分布函数；在边缘分布函数分析的基础上，采用gumbel copula函数构建植被出现频次—ndvi联合概率分布函数，生成地下水埋深—ndvi联合概率分布。参考图4，为本发明实施例的吐鲁番盆地地下水依赖型植被ndvi与地下水埋深联合概率分布图。
[0043]
第四步：设置吐鲁番盆地合理的地下水依赖型植被保护目标，v
t
＝0.15，满足v＞0.15；基于地下水埋深—ndvi联合概率分布，分析不同地下水埋深区间实现植被保护目标v＞0.15的条件概率。参考图5，为本发明实施例的吐鲁番盆地实现地下水依赖型植被保护目标的条件概率图。以实现地下水依赖型植被保护目标条件概率最大值对应的地下水埋深作为地下水生态埋深。不同地下水埋深区间实现植被保护目标的条件概率的最大值为0.47，对应地下水埋深区间6
‑
7m。进而得到吐鲁番盆地地下水依赖型植被的地下水生态埋深为6
‑
7m。
[0044]
本发明所述一种基于联合概率分布的干旱区地下水生态埋深分析方法的更多应用是：可以用于分析干旱区内陆河流域地下水依赖型植被的地下水生态埋深，为干旱区内陆河流域生态需水评估提供一个有效工具，支撑地下水依赖型陆生生态系统的保护和管理，为保障干旱区内陆河流域自然绿洲植被生态需水、提高水资源利用效益、优化水资源配置格局提供技术支撑和科学依据。
[0045]
需要说明的是，以上说明仅是本发明的优选实施方式，应当理解，对于本领域技术人员来说，在不脱离本发明技术构思的前提下还可以做出若干改变和改进，这些都包括在
本发明的保护范围内。