1. 引言

随着我国城市化进程的加快，越来越多的地区为了提升区域竞争力，创造更多的社会经济价值，加大了对区域自然资源的开发和利用力度；随即区域资源短缺、生态环境破坏等现象频繁发生，进一步影响到区域整体的可持续发展。水资源是区域生态环境内部的基本自然资源 [1]，对其合理的利用有助于区域社会经济、生态环境以及人类群体本身健康良性的发展 [2] [3] [4] [5]；但目前国内许多快速城市化地区(诸如西安市)水资源同样面临日益短缺的问题。20世纪90年代，生态足迹(Ecological Footprint)模型，由加拿大生态经济学家William和Wackernagel提出、发展，并被广泛应用于区域资源可持续利用发展的综合评估研究当中 [6] [7]。同时该模型方法也为水资源(单项资源)的可持续利用定量评价提供了新的思路 [8] [9]。国内学者在有关区域性水资源生态足迹方面也做了诸多的研究：卢艳等 [10] 通过分析计算发现河南省大多城市水资源出现生态赤字；张养安等 [11] 对陕西省汉中市水域生态承载力进行了研究，结果表明：汉中市水域生态承载力整体处于安全状态。张义等 [12] 提出了改进型水资源生态足迹模型，并计算了广西的水资源生态足迹。王文国等 [13] 通过分析计算发现：四川省人均水资源生态足迹呈整体上升趋势但有一定的生态盈余。刘子刚等 [14] 通过将水生态足迹划分为水产品生态足迹、水资源生态足迹与水污染生态足迹三部分，建立了相应的水生态足迹模型。西安市作为典型的内陆快速城市化地区，而且也是中国西部地区重要的中心城市，它的可持续发展状态将影响和辐射至整个西北地区。因此，分析和评估该区域的水资源可持续发展状况，具有重要的现实意义。

本文以陕西省西安市为研究区域，评估计算了2010~2020年西安市水资源生态足迹和生态承载力，并对其空间分布特征进行分析，同时应用LMDI法分析了影响西安市水资源生态足迹的因素。文中研究结果可为西安市水资源科学的管理、有效的利用以及区域可持续发展的评估提供基础依据，同时也可为同类型城市水资源的可持续评价提供参考。

2. 材料与方法

2.1. 研究区概况

陕西省西安市位于我国西北东部、关中平原中部，北濒渭河、南依秦岭，位于33˚39'~34˚44' N、107˚41'~109˚49'E之间，地势西北低东南高，辖境东西长约204 km，南北宽约116 km，土地总面积约为10,108 km²，西安市自2017年代管西咸新区，故2017年及以后的土地总面积视为10,752 km²。暖温带半湿润大陆性季风气候在市内的平原地区分布，四季分明(见图1)。年降水量522.4~719.5 mm，由北向南递增，全年中7月、9月为两个明显降水高峰月。据陕西省2020年水资源公报显示，西安市降水量位于陕西省所有地级市的第六位，水资源总量位于第八位，水资源总量较为匮乏。

2.2. 数据来源

本文研究区人口、GDP等数据来源于《西安市统计年鉴(2010~2020)》和《陕西省统计年鉴(2010~2020)》；水资源总量等相关数据来源于《陕西省水资源公报(2020)》和《西安市水资源公报(2010~2020)》。

2.3. 研究方法

2.3.1. 水资源生态足迹

水资源生态足迹反映区域人口、社会发展对水资源的占用情况。具体涵义是在一定的人口规模和经济发展水平条件下，将人类生存和社会发展等所必需的水资源量折算为所必要的生态生产性面积，再通过水资源产量因子和均衡因子等参数的校正，得到区域内可以相互比较的均衡值 [1]：

$WF=T\times a_{wf}=T\times \epsilon \times \left(W/a\right)$(1)

式中：WF为区域水资源总生态足迹(hm²)；T为区域的人口数量； $a_{wf}$为区域的人均水资源生态足迹(hm²/人)； $\epsilon$为区域水资源均衡因子；W为区域人均消耗的水资源量(m³)；a为世界水资源的平均生产能力(m³/hm²)。本文水资源均衡因子 $\epsilon$取值5.19 (黄林楠等 [15])；水资源全球平均生产能力a取值3140 m³/hm²(孟丽红等 [16])。

2.3.2. 水资源生态承载力

水资源生态承载力反映区域水资源的供给状况。具体涵义是在一定的历史发展阶段，水资源最大供给能够支持维系有限发展目标的资源–环境–社会–经济可持续发展最大规模 [17]，计算公式如下 [18]：

$WC=T\times a_{wc}=0.4\times \theta \times \epsilon \times \left(Q/a\right)$(2)

式中：WC为区域的水资源承载力(hm²)； $a_{wc}$为区域的人均水资源承载力(hm²/人)； $\theta$为区域水资源产量因子(该区域水资源平均生产能力与世界水资源平均生产能力之比)，文中取值1.5 (全球生态足迹网于2017年发布的《Working Guidebook to the National Footprint Accounts》 [19])；Q为区域的水资源总量(m³)； $\epsilon$和a同式(1)。已有研究表明：区域水资源开发利用率超过30%~40%将会引发生态环境的恶化 [20]，因此水资源承载力中需预留60%用于维护生态环境稳定 [21]。

2.3.3. 水资源生态盈亏

水资源生态盈亏反映区域水资源的供需平衡状态，用于判别和衡量区域水资源的可持续利用情况 [18]。具体涵义是生态承载力与水资源生态足迹的差值，计算公式如下：

$D=WC-WF$(3)

式中：D为水资源的生态盈亏。D大于0，表明区域内水资源的供大于需，处于盈余状态，水资源可持续利用；D等于0，表明区域水资源处于供给与需求的临界状态；D小于0，表明区域内水资源的供不应求，处于赤字状态，水资源不可持续利用。

2.3.4. 万元GDP水资源生态足迹

万元GDP水资源生态足迹用来判断与衡量区域的用水效率 [18]，具体涵义为区域水资源生态足迹与GDP的比，计算公式如下：

$R=WF/GDP$(4)

式中：R为区域的万元GDP生态足迹。R值越大，区域水资源利用效率越高，反之亦然。

2.3.5. 水资源生态压力

生态压力指数反映人类的活动对生态系统的干扰程度 [22]，具体涵义为区域水资源生态足迹与水资源承载力之比比，计算公式如下：

$EPI=WF/WC$(5)

式中：EPI为区域水资源生态压力。当EPI_W小于1时，表明人类活动对生态系统的干扰强度尚未超过区域生态系统自身的反馈阈值，生态系统仍处于平衡状态。而EPI大于1，表明人类活动对生态系统的干扰超过了生态系统自身的反馈阈值；并且其值越大，干扰越强，区域就越有可能面临生态环境安全和生态系统崩溃的问题 [23]。

2.3.6. 对数平均迪氏指数法

近年来在资源环境领域，因素分解法逐渐成为一种重要的研究方法，该方法能够提供一个完整的框架来分析各变量的影响程度 [24]。LMDI模型是由Ang等 [25] 提出的将研究对象划分为若干影响因子的分解法 [26]，不产生无法分解的残差，适用范围较广，用于测度不同因子对研究对象的影响及其驱动力分析 [27]。

基于LMDI模型首先将水资源生态足迹模型进行指数分解，如下：

$W{F}_t={\displaystyle \underset{i}{\sum }W{F}_{it}}={\displaystyle \underset{i}{\sum }\frac{W{F}_{it}}{W{F}_t}\cdot \frac{W{F}_t}{g_t}\cdot \frac{g_t}{p_t}\times p_t}={\displaystyle \underset{i}{\sum }{s}_{it}{t}_t{e}_t{p}_t}$(6)

$\Delta W{F}_t=W{F}_t-W{F}_0={\displaystyle \underset{i}{\sum }{s}_{it}{t}_t{e}_t{p}_t}-{\displaystyle \underset{i}{\sum }{s}_{i0}{t}_0{e}_0{p}_0}=\Delta W{F}_s+\Delta W{F}_t+\Delta W{F}_e+\Delta W{F}_p$(7)

其中， $s_{it}=W{F}_{it}/W{F}_t$， $t_t=W{F}_t/g_t$， $e_t=g_t/p_t$，将区域水资源生态足迹，按照生产、生活和生态等3类生产账户进行分解。式中： $W{F}_{it}$为第t时期第i类水资源的生态足迹(hm²)； $W{F}_t$为第t时期水资源总生态足迹(hm²)； $g_t$为第t时期区域国内生产总值GDP (亿元)； $p_t$为第t时期区域常住人口数(万人)； $W{F}_t,W{F}_0$分别是第t年和基准年的水资源生态足迹量(hm²)。

将西安市水足迹变化的驱动因素分解为结构因素、人口因素、经济因素和技术因素四方面： $s_{it}$为第t时期第i类水资源利用占水资源总生态足迹比重，代表水足迹结构因素； $t_t$为第t时期单位GDP水资源生态足迹，代表水足迹技术因素； $e_t$为第t时期人均GDP，代表水足迹经济因素； $p_t$代表水足迹人口因素； $\Delta W{F}_t$为水资源总生态足迹变化量(hm²/a)。接下来对这四个因素分别引起的水足迹变化量进行计算，即利用LMDI加和分解对(7)式进行分解，如下：

$\Delta W{F}_s=\ln \frac{s_{it}}{s_{i0}}\underset{i}{{\displaystyle \sum }}\frac{W{F}_{it}-W{F}_{i0}}{\ln W{F}_{it}-\ln W{F}_{i0}}$(8)

$\Delta W{F}_t=\ln \frac{t_t}{t_0}\underset{i}{{\displaystyle \sum }}\frac{W{F}_{it}-W{F}_{i0}}{\ln W{F}_{it}-\ln W{F}_{i0}}$(9)

$\Delta W{F}_e=\ln \frac{e_t}{e_0}\underset{i}{{\displaystyle \sum }}\frac{W{F}_{it}-W{F}_{i0}}{\ln W{F}_{it}-\ln W{F}_{i0}}$(10)

$\Delta W{F}_p=\ln \frac{p_t}{p_0}\underset{i}{{\displaystyle \sum }}\frac{W{F}_{it}-W{F}_{i0}}{\ln W{F}_{it}-\ln W{F}_{i0}}$(11)

分解结果若大于0，则表明该因素起到推动作用，否则，表明该因素具有抑制效果 [28]。

3. 结果与分析

3.1. 西安市2010~2020年间水资源生态足迹各指标分析

西安市地处内陆地区，水资源并不丰富，但随着经济的发展，产业的扩张，水资源生态足迹呈逐年平稳上升又缓慢下降(自2019年起)趋势；如图2所示：本文研究时段内，其最小值为2.94 × 10⁶hm²(2010年)，最大值为3.70 × 10⁶hm²(2018年)，涨幅达25.87%，年均涨幅2.92%。

资料来源：《西安市水资源公报》《西安市统计年鉴》。

西安市是典型的内陆城市，其水资源总量，强烈依赖于降水量。如图3所示，本文研究时段内，西安市水资源生态承载力最大值为3.05 × 10⁶hm²(2011年)，最小值为1.87 × 10⁶hm²(2012年)，其他年份整体平稳波动。如表1所示：2011年西安市降水总量为83.42 × 10⁸m³，水资源总量为30.74 × 10⁸m³，同比上年增加27.57%和27.55%；该年分水资源生态承载力也同比上年增加27.62%，属于丰水年。将这11年间区域水资源生态承载力与区域降水量之间进行回归分析，如图4所示：两者之间的相关系数达0.6557，呈显著正相关。数据分析表明：西安市水资源的承载力同样受区域降水量(自然因素)的强烈影响。

如表2所示：2010~2020年11年间，西安市的水资源生态盈亏状态，大多数年份呈现生态赤字，只有2011年有少量的盈余。在水资源生态足迹逐年增加和降水量年际波动的交互影响下，区域生态亏损加剧。生态赤字最大值出现在2018年，其值为−1.53 × 10⁶hm²，生态盈余唯一出现在2011年，约为0.02 × 10⁶hm²。西安市水资源的供给，对自然的依赖比较严重，随着经济的不断发展，对水资源的需求进一步加剧，区域水资源的供需失衡，甚至区域水资源的短缺，都是未来该区域可持续发展将要面临和亟待解决的重要问题。

资料来源：《西安市水资源公报》。

如图5所示，西安市万元GDP水资源生态足迹最大值出现在2010年，为9.06 × 10⁻⁶hm²，最小值出现在2020年，为3.37 × 10⁻⁶hm²。从各年数据来看，万元GDP水资源生态足迹总体上呈缓慢下降趋势，下降幅度最大的年份在2011年、2017年。究其原因，一方面西安市在发展经济的同时，民众节水意识增强，水资源利用效率提高；另一方面政府在各类污水的处理上取得了较好的成效。

如图6所示，西安市2010~2018年生态压力指数呈波动上升趋势，在2019年后有显著下降的趋势。最小值出现在2011年，约为0.99；最大值出现在2018年，约为1.70。上述情形表明：随着西安市社会经济活动的不断发展与加强，对水资源生态系统的干扰越来越大，面临水资源生态安全和水资源生态系统崩溃的风险和可能性也越来越大，但是在2019年之后，生态压力指数也有逐渐降低的趋势，这表明政府和社会对区域生态以及水资源可持续发展越来越重视。

3.2. 西安市2010~2020年间各生产账户生态足迹分析

将区域水资源生态足迹，按照生产、生活和生态等3类生产账户进行分解，并分析各账户内的水资源生态足迹随时间的变化趋势。如图7所示，各账户中生产账户占比最大，约占总账户的60%~75%。其次为生活账户，约占20%~25%，最后是生态账户，约占总账户的3%~15%。但随着可持续发展理念的不断落实与深化，生态文明建设的不断推进，西安市生态账户的比重逐年上升，且增速较快，生产账户比重有所下降，生活账户基本不变。

3.3. 西安市2010~2020年间三大产业水资源生态足迹分析

将区域水资源生态足迹中生产用水生态足迹，按照三大产业进一步细化区分，并分析各产业人均生产用水

生态足迹随时间的变化趋势。如图8所示：西安市第一产业人均水资源足迹整体上呈现下降的趋势(2019~2020年下降幅度最大，其余年份间变化平稳)；第二产业呈上升后陡降的态势(2020年受新冠肺炎疫情影响第二产业水资源生态足迹大幅下降)；第三产业的前期平稳，2016年略有下降，随后逐渐上升，2019~2020年增幅最大，整体变动趋势与第二产业相反。

3.4. 西安市2010~2020年间水资源生态足迹变化驱动因素分析

为了进一步探讨西安市水资源生态足迹变化的驱动因素，本文以2010年为基准年，采用LMDI模型的加和形式对其进行因素分解，计算2010~2020年西安市水资源生态足迹的分解因素，并分析其趋势如下。

图9可知，水足迹经济因素、人口因素和结构因素对西安市水资源生态足迹的增长起到正效应，水足迹技术因素起到负效应。经济因素是西安市水资源生态足迹增长的主要因素，结构和人口因素次之。上述结果表明，程度十分有限，且在2019年后有所下降；对于技术因素来说，整体呈下降趋势，对水资源生态足迹过快增长的抑制作用明显，反映了西安市加大污水治理和提高水资源重复回收利用等方面的技术投入。这就说明通过技术手段提高单位GDP的用水量，可以在一定程度上减少水资源生态足迹。技术因素虽然很大程度上抑制西安市水资源生态足迹的增长，但是其抑制作用远远比不上经济因素与人口因素引起的拉动作用，从而导致西安市水资源生态足迹整体上呈现增长趋势。

4. 讨论

总体而言，西安市是一个典型水资源不足的内陆城市，为了缓解区域水资源短缺，改善水生态赤字的情况，首先：应该做好节水的相关规划，采取切实可行的措施，实现节水的目标 [29] [30]；其次：应该加强环境保护与治理，尤其需要进一步加强渭河、黑河等现有水资源及其重要供水河流的环境保护与污染治理，改善水环境，提高水功能，控制地表水污染、减少地下水超采、增大污水和雨水的回用量，使有效供水量增加，提高城市水资源和水环境的承载力 [31] [32]；再次：由本研究对西安市水资源生态足迹在各种账户、各个产业之间变化趋势的对比分析，也可以看出产业调整对水资源需求的影响。今后西安市也应该在充分考虑社会经济现实情况的基础上，进行产业结构的调整和优化；促成水资源在各产业之间的合理分配，提高其利用效率 [14] [33]，从而促进西安市整体社会的可持续发展。

从水资源生态足迹驱动因素分析方面来看，经济发展因素是导致水资源生态足迹增长的主要原因，技术因素是降低水资源生态足迹的关键因素，西安市应加大清洁生产、节能降耗、回收利用等技术方面的投入。此外，在上述各项措施仍无法满足西安市整体用水需求时，亦可适当考虑实施开辟当地新水源或者采取区外调水的措施 [31]。从而统筹兼顾，多种渠道谋取解决西安市区域水资源短缺的问题。

5. 结论

1) 2010~2020年西安市水资源生态足迹呈逐年上升又下降(自2019年起)趋势；总水资源生态足迹中生产用水生态足迹占比最大，生态用水生态足迹以较快增速逐年上升。三类产业人均水资源生态足迹中，一产呈现缓慢下降趋势，2019~2020年下降幅度最大，二产呈上升后陡降的态势(2020年受新冠肺炎疫情影响大幅下降)；三产前期平稳，随后逐渐上升，2019~2020年增幅最大，整体变动趋势与二产相反。万元GDP水资源生态足迹总体呈缓慢下降趋势。

2) 西安市水资源生态承载力受区域降水影响较大，二者呈明显的正相关关系；除2011年外，西安市的水资源生态盈亏大多年份处于生态赤字状态，而且赤字在2019年前有上升的趋势，在2019年后赤字情况得到一定的缓解。在研究时段内西安市水资源生态压力指数除个别年份外，大多数年份呈上升的趋势，但整体有下降趋势，水资源生态压力有所放缓。

3) 影响西安市水资源生态足迹的最主要驱动因子为经济因素。技术因素对水资源生态足迹的增加起到抑制作用，通过提高单位GDP的用水量，可以减少水资源生态足迹。

各种迹象表明：西安市随经济的发展，用水需求增加，水资源生态足迹逐年扩大；但水资源的供给对自然的依赖强，供需矛盾突出，长期处于生态赤字状态，未来水资源短缺或将成为区域可持续发展重要的制约因素。

基金项目

国家自然科学基金：多源观测信息与陆表水碳通量过程模型的数据同化研究(41301451)；中央高校基本科研业务费：基于路径和空间异质运算的遥感影像地表信息提取分析方法研究(2452018144)。

参考文献

参考文献