丽江黑龙潭地下河系统的水均衡分析与预测评价.pdf

丽江黑龙潭地下河系统的水均衡分析与预测评价任世川1,2,3，杨晓艳1,2,3，杨颖彬1,2,3，刘海峰1,2,3，杨帆1,2,3（1.自然资源部高原山地地质灾害预报预警与生态修复重点实验室,云南昆明650216；2.云南省高原山地地质灾害预报预警与生态修复重点实验室,云南昆明650216；3.云南省地质环境监测院,云南昆明650216）摘要：滇西北丽江黑龙潭流量大幅减少、间歇式断流问题日趋频繁，已成为地方经济社会发展的重要制约因素。在对丽江盆地区域水文地质条件和黑龙潭地下水系统动态调查的基础上，考虑社会经济发展过程中人类工程活动的影响，以近 70 年的数据为基础，采用水均衡法定量分析黑龙潭地下水系统流量的动态变化及其影响因素，得出不同时段黑龙潭断流的主要原因，并根据近年黑龙潭断流的影响因素及变化趋势，预测未来黑龙潭流量的变化趋势，20212030 年，黑龙潭将成为季节性泉水，2030 年之后，黑龙潭泉群将长期处于断流状态。该认识可为黑龙潭的地下水动态恢复和丽江城区地下水水资源开发及环境保护提供依据。关键词：地下水系统；水均衡分析；黑龙潭断流；丽江市中图分类号：P641.134文献标识码：A文章编号:10014810（2023）06118310开放科学（资源服务）标识码（OSID）：0引言岩溶水系统是受基础地质条件和岩溶发育程度控制的，随着气候条件、人类活动等诸多因素变化而不断演化的动态系统，泉水的动态变化是岩溶水系统演化过程的反映，其影响因素具有复杂多样、时空变化大、相互影响强及存在不同滞后期等特点，特别是人类活动作为环境变化的主要驱动力之一，短时间内造成环境改变的速率及强度远超一般的自然演替，其产生的环境效应及诱发水文地质问题极难辨别和量化，一直是国内外学者和地方政府关注的焦点16。丽江盆地位于晚新生代地壳强烈抬升区，发育着典型的隐伏、半隐伏岩溶水系统，孕育了黑龙潭、西关泉等岩溶大泉（泉群），这些泉水既是丽江市重要的生活水源，也是丽江的亮点和旅游名片。自2000 年开始，黑龙潭断流现象不断发生，并有日趋加剧趋势。为查找泉水断流原因，众多学者从不同角度开展了相关研究，普遍认为降水量减少是黑龙潭断流的主要原因713，而人类活动对岩溶水系统的改造和影响研究相对较少。因此，研究人类活动对岩溶水循环过程及其资源环境演化的驱动机理及响应变化，查明黑龙潭断流的主要原因，并预测变化趋势，有利于科学推进大中型城市岩溶水资源的开发利用、生态环境修复和防灾减灾等活动6。1地下河系统的边界及结构1.1系统概况和边界条件黑龙潭地下水系统位于丽江盆地北东侧，呈倒三角形由北东向南西斜插入丽江盆地，是丽江盆地地下水系统最为重要的岩溶水子系统，为三面封闭一面开放的地下水系统，汇水总面积 172.69km2。系第一作者简介：任世川（1979），男，工程硕士，正高级工程师，主要从事水工环地质方面的调查研究工作。E-mail：。收稿日期：20230101第42卷第6期中国岩溶Vol.42No.62023年12月CARSOLOGICASINICADec.2023任世川，杨晓艳，杨颖彬，等.丽江黑龙潭地下河系统的水均衡分析与预测评价J.中国岩溶，2023，42（6）：1183-1192.DOI：10.11932/karst20230604统西侧边界沿山盆界线分布，边界以西为清溪龙潭地下水子系统。根据收集资料78、1112、14和本次开展的钻孔、示踪试验、水化学分析、同位素测定等技术方法验证二者之水力联系微弱，属隔水边界；北侧边界为地表和地下分水岭，本项目示踪试验证明九子海一带落水洞仅补给黑龙潭；东侧边界为丽江金棉断裂，由隔水断层和隔水层构成隔水边界12；南西侧以丽江古城外围为界，为透水边界，地下水侧向补给丽江城地下水子系统，多次示踪试验证明黑龙潭与古城内泉水水力联系密切（图 1）。系统北部九子海至腊日光一带属地下水补给区，面积约 122.9km2，以低矮缓丘、漏斗、洼地、落水洞等地表岩溶组合地形为主，含水层为中三叠统北衙组灰岩、白云质灰岩，入渗系数 0.411，枯季径流模数 6.07Ls1km2，富水性强，但不均匀。中南部腊日光至象山一带属地下水径流区，面积约 47.3km2，地形平缓，坡度 515，季节性溪沟发育，沟水由顶部向东西两侧径流，出露地层以古近系始新统丽江组砂砾岩夹泥岩为主，入渗系数 0.367，枯季径流模数4.63Ls1km2，富水性较强，地下水多补给下伏北衙组岩溶含水层，少部分以泉水形式排泄。象山至古城一带为地下水排泄区，面积约 2.5km2，地势平缓，出露第四系松散堆积层，以冰水砾石层为主1416。1.2系统结构黑龙潭地下水系统是由多种岩溶空间组合而成的结构复杂的导蓄水岩溶空间系统。浅表部以各种溶蚀裂隙、风化裂隙网络、表层岩溶带等岩溶空间为主，局部有稀少的溶洞、岩溶管道或较大的溶蚀裂隙（缝）及第四系风化残余土中的孔隙等含水介质空间分布，具有较好的蓄水功能，并具有一定的输水作用；中深部则以分布范围狭窄但发育规模较大的岩溶管道或深大岩溶裂隙（缝）为主，具有较强的输水功能，较差的蓄水作用。2地下水径流、排泄过程及动态特征2.1径流、排泄过程大气降水入渗补给是黑龙潭地下水系统唯一补给来源，受两侧隔水层控制，地下水通过节理、裂隙、溶隙等空间向中部九子海至腊日光一带汇集，由以裂隙流为主逐渐转变为以管道流为主，通过深部岩溶管道系统快速向丽江盆地边缘的象山一带径流，九子海（2802m）至黑龙潭（2420m）距离约 14.5km，高差 382m，水力坡度较大，水动力条件较好，循环交替迅速，雨季，管道流到达时间一般 1015 天，裂隙流到达时间一般 24 个月。排泄区南侧为古近系碎屑岩和少量河湖相沉积物，东侧为碎屑岩和玄武岩，均属隔水层，具有较好的地下水储存条件，大量的地下水在象山一带储存，形成承压水。西侧为第四系冰水堆积物，透水性极好，但均被后期的第四系河湖积沉积物覆盖，排泄不畅，受其控制黑龙潭地下水系统部分地下水沿断裂带上升以泉水形式排泄，形成黑龙潭泉群；部分地下水通过第四系冰水堆积物孔隙侧向补给丽江城地下水系统，以蒸发排泄为主。因此，黑龙潭地下水系统属典型的“自产二流型”岩溶水系统，以泉水流出和侧向排泄是该地下水系统的主要排泄方式7（图 1）。2.2泉水动态和断流特征。根据黑龙潭系列监测资料（19902011 年，2012 年 2 月断流）13，黑龙潭流量年际变化剧烈，多年平均年径流量为4040 万m3，最大径流量为8038 万m3（2001 年），最小年径流量约 1277 万 m3（2007 年），相差达 6 倍。年内流量动态变化也较剧烈，呈现明显的年周期特征，在一个周期内，动态表现为单峰单谷，属波态型岩溶泉水，变幅一般可达 35 倍。据丽江县志和相关资料记载，历史上黑龙潭共出现 27 次不同程度的间歇式断流713、17，自 17291960 年的 200 多年间共断流 9 次，平均约 25.7 年断流一次，每次断流时间一般 23 个月，无跨年断流；19611990 年断流 5 次，平均 6 年断流一次，其中 3次跨年断流，断流时间最长约 29 个月，断流情况呈现加剧趋势；19912020 年断流 13 次，平均约 2.3 年断流一次，其中 2 次跨年断流，断流时间最长 43 个月，断流情况明显加剧（图 2）。3地下河水平衡方程建立与参数处理3.1水平衡方程建立3.1.1黑龙潭地下水系统水均衡方程式：KP=w+Qh+Qg+Qs式中：入渗系数；K补给面积（m2）；P降水量1184中国岩溶2023年岩溶水裂隙水孔隙水碳酸盐岩灰质角砾岩长石砂岩第四系松散堆积丽江盆地边界一级地下水系统界线二级地下水系统界线富水块段边界断层推测断层地层界线推测地层底部边界剖面位置及编号泉或泉群左泉名,右流量 Ls1钻孔及名称落水洞地下水流向024 km-水文地质剖面图1 8002 1002 4002 7003 0003 3001 8002 1002 4002 7003 0003 300210160标高/m标高/m蒸发排泄大气降水地下河管道大气降水西关泉群漾弓江漾弓江黑龙潭FFFFFFT2T2T2T2T2T2T3EEQ腊日光巴忙山二级排泄面一级排泄面三级排泄面丽江市九子海玉龙县鹤庆盆地北端图1丽江盆地水文地质及地下水系统关系图Fig.1RelationshipbetweenhydrogeologyandgroundwatersystemintheLijiangbasin第42卷第6期任世川等：丽江黑龙潭地下河系统的水均衡分析与预测评价1185（m）；地下水储存变化量（m3a1），是蓄水深度H 的函数，即 f（H）；Qh黑龙潭泉群流量（m3a1）；Qg系统内其它泉水排泄量（m3a1）；Qs侧向排泄量（m3a1）（补给丽江城地下水系统）。其中，Qs受控于丽江城地下水系统环境要素变化，因此，要定量分析黑龙潭流量变化的影响因素，必须分析丽江城地下水系统水平衡。3.1.2丽江城区地下水系统水平衡方程式：Qnc=w+Qy+Qr+E+QqQjQdQqcQg式中：Qnc黑龙潭侧向补给量（m3a1），等于上式的Qs；为地下水储存量变化（m3a1），为蓄水深度H 的函数，即 f（H）；Qy越流排泄量（m3a1）；Qr人工排泄量（m3a1）；E蒸发排泄量（m3a1）；Qq其它变化量（m3a1），城市化带来的其它影响（难以细分和量化）；Qj降水入渗补给量（m3a1）；Qd地表河系渗漏补给量（m3a1）；Qqc其它区域侧向补给量（m3a1）；Qg灌溉水入渗补给量（m3a1）。其中，降水入渗补给量均采用以下计算公式Q=KP式中参数同前。蒸发强度计算采用阿维扬诺夫公式计算：=0(1h/l)n式中：蒸发强度(m)；0水面蒸发强度(m)；h水位埋深(m)，取值 1.402m；l极限蒸发深度，取经验值 5.3m；n蒸发指数，取 3。再据下式计算蒸发排泄量：E=K式中：E蒸发排泄量（m3a1）；蒸发强度(m)；K排泄区面积（m2）。地表河系渗漏补给量采用目前我国广泛应用的考斯加科夫经验公式计算：S=0.01AQ1mn式中：S渗漏损失，m3(skm)1；A透水性系统，取经验值 0.4；Qn通过渠道的净流量，m3s1，取年平均值；m透水性指数，取经验值 0.3。再根据下式计算地表河系渗漏补给量Qd=31536000SL式中：L丽江城地下水系统内地表河系长度（km），其它参数同前。3.2水均衡要素及变化特征根据水均衡方程式，影响黑龙潭流量变化的要素主要包括降水量、蒸发量、农田灌溉水量、与地表河系交换量、盆地入渗面积、人工开采量、丽江城市横向发展对蒸发排泄量的影响和城市纵向发展对地下水环境的影响。丽江盆地在唐代就已经有了原始的村落，但至解放前一直以农耕为主，社会经济发展缓慢，当时区内人口稀少，村落稀疏，人类活动对环境影响微弱，地下水环境长期处于自然状态，后期才有以现有河、泉的简单利用活动，对地下水环境改造几乎为零17。20 世纪 5090 年代，盆地内社会经济发展仍较为缓慢，多属满足本地生产和居民生活需求的基础建设18，如水利、道路等基础设施有了明显地改善；同时，该阶段境内木材经济快速发展，森林植被破坏严重，森林面积减少了 42.1%19，且原始森林几乎消耗殆尽，各种树木平均年龄从 450 年下降至 69 年17。20 世纪 90 年代后期，随着旅游业的快速兴起，盆地内社会经济进入快速发展时期，地下水环境变化显著，尤以城乡聚落的扩展、铁（公）路、水利等基础建设用地的迅速增加，耕地、冰川和积雪面积的大幅度减少为特征。同期，生态环境治理力度不断增强，丽江市森林覆盖率迅速增长，面积增加了 157.7%18、2026。全境各环境要素变化详见表 1。3.3数据筛选与处理此次水平衡计算数据涉及的年限为 1950 至2017 年。新中国成立后，丽江市各部门根据业务范围陆续展开相关的监测、统计、整理工作，地下水水平衡方程式所涉及的参数数据才逐渐被提取出来，但各参数出现的时间和间隔跨度并不相同，少数参010203040501729190019011930193119601961199019912020(年)年断流概率/%最长断流时间/月图2黑龙潭历史断流情况变化趋势曲线图Fig.2Trendcurveofthehistoricaldrying-upoftheHeilongtanspring1186中国岩溶2023年数从未被关注，相关数据严重不足。若按年度分析黑龙潭流量影响因素及权重，必须对数据进行处理和补充。因而，本次计算主要采用以下四种方法补齐相关数据。3.3.1离散处理法对工作区仅有数年连续数据，对年份不全的参数采用离散函数法进行补充，如九子海站降水量仅有 20012017 年的数据，通过对比丽江站与九子海站 17 年数据，发现两站降水量呈线形关系，仅个别年份(2014 年)不符合该规律（图 3），故根据两站降水量相关关系，将九子海气象站降水量资料延长至1950 年。3.3.2插值法对仅有数年断续数据的参数，根据丽江盆地地下水环境发展变化以插值法补充完善，如灌溉用水量、地下水开采量、地表河系入渗补给量等。其中，丽江盆地灌溉水用水量仅有 2 年数据，分别为 1980年的 3100 万 m3和 2010 年的 2160 万 m39，1990年以前丽江盆地农业用水量变化较小，始终延用1980 年的 3100 万 m3，1990 年之后，随着城市迅速发展，农业用水量随之下降，以 2010 年的 2160 万 m3为基础，用插值法补全其它年份农业用水量，入渗系数取经验值为 0.3。3.3.3类比法对一般的泉水仅有本项目的长观数据或偶测数据，少量泉水有近几年云南省水文水资源局丽江分局长观资料，这些泉水变化虽然直接影响黑龙潭排泄量，但其地下水环境未有明显改变，主要受降水量控制，因此，采用数模处理方法进行类比处理，得到其系列数据。根据丽江盆地地下水补给、径流、排泄特征可知，表 1丽江盆地社会经济发展和地下水环境变化情况统计表Table1StatisticsofsocioeconomicdevelopmentandgroundwaterenvironmentchangeintheLijiangbasin变化因素1940s1950s1960s1970s1980s1990s2000s2010s备注丽江县（古城区和玉龙县）国民生产总值/亿元0.30.60.71.23.342.5194.2收集数据极值丽江城人口/万人3.906.9114.6222.70收集数据极值城市建成区面积/km21.702.103.227.5530.5848.90收集数据极值公路里程/km0153351335.95内耕地面积/km2129.0989.0577.44水利工程蓄水量/万m30300011483城市供水规模/万m30071.480.58763504收集数据极值地下水开采量/万m303737373751.05589.831412.36丽江市植被覆盖率/%47.2030.4027.3043.661.1070.34收集数据极值九子海站年平均降水量/mm1289.7 1199.4 1256.5 1218.5 1215.1 1319.9 1274.31108.8极值差/mm507229446441518523489777年平均蒸发量/mm2003.6 2118.1 2178.5 2240.3 2079.8 1551.61528.4年平均气温/12.512.612.612.712.813.313.7年平均日照数/h2453.3 2572.8 2518.8 2450.5 2420.9 2366.72403.6玉龙雪山冰雪面积/km2155.6397.4085.34113.8417.97收集数据极值冰川面积/km211.618.505.303.98收集数据极值02 0004 0006 0008 00010 00012 00014 00016 0002000 年2001 年2002 年2003 年2004 年2005 年2006 年2007 年2008 年2009 年2010 年2011 年2012 年2013 年2014 年2015 年2016 年2017 年0.1 mm丽江站(mm)九子海站(mm)图3丽江站与九子海站 2001 年至 2017 年降水关系曲线图Fig.3PrecipitationrelationshipbetweenLijiangstationandJiuzihaistationfrom2001to2017第42卷第6期任世川等：丽江黑龙潭地下河系统的水均衡分析与预测评价1187西关龙潭地下水补给来源主要包括 2 种类型，一是补给区降水入渗补给，根据丽江站降水量系列数据求得其降水入渗补给量。二是盆地越流补给，1990年以前，主要受降水控制，以数模处理的方法进行类比处理，得到其系列数据；1990 年以后，受人类工程活动影响，越流补给量逐年减少，以 1978 年和 2017年的数据为基础，采用插值法获得系列数据。3.3.4推算法对工作区少部分无相关数据的计算参数，以现有与之有关的资料进行推算补充，如近年来人工地下水开采量、丽江城地下水系统耕地面积数据、地表河系渗漏补给公式的部分参数等。根据丽江市人民政府网公开信息丽江：全面落实最严格水资源管理制度和丽江市水井排查结果，截至 2019 年 12 月 19 日，丽江盆地共有水井 18169口，其中民井 17490 口，按 4 人计，每人每天用水定额为 100L，每 年 用 水 量 为 255.35 万 m3；深 水 井679 口，全部关停，根据古城区 2017 年关停开采井节约用水量推算，深水井取水量为 1150 万 m3，根据丽江自来水公司 2019 年和 2015 年供水量之差推算深水井取水量为 1155.4 万 m3。因此，2018 年，深井地下水开采量约 1412.36 万 m3。其中丽江城地下水系统内为 1300.7 万 m3。1990 年以前丽江城地下水系统内耕地约占盆地的耕地的 50%，1991 至 2018 年，随着社会经济发展，系统内大量耕地被建设用地替代，因此，耕地减少面积以系统内城市建成区面积增加量进行计算，最终得出丽江城地下水系统内灌溉水入渗补给量。地表河系渗漏补给公式中 A 值 0.7，雨季丽江城地下水系统内地下水水位上升，受其影响，河流、水渠渗漏补给地下水量减小，综合考虑，将 A 值调整为 0.4；1991 年至 2000 年，随着经济发展，大约 70%的河、渠逐渐修建了三面光工程，以插值法计算；2000 年之后，以 30%的交换量进行计算。4水均衡结果和可信度分析4.1计算结果4.1.1黑龙潭地下水系统水资源变化20 世纪 5080 年代，黑龙潭泉群排泄量呈缓慢减少的趋势，由占补给量的 54.72%减少到 44.24%，90 年代因降水丰富，黑龙潭排泄量略有回升，占补给量的 48.17%，20002010 年，黑龙潭排泄量呈断崖式下降，由占补给量的 44.74%减少到 3.42%，侧向排泄量变化趋势与黑龙潭排泄量呈相反的增长趋势（表 2）。4.1.2黑龙潭断流主要原因20 世纪 5080 年代，因人类活动造成区域植被更替和减少，导致自然蒸发能力增强，地下水蒸发排泄量持续增加 1082.74 万 m3，占系统多年平均补给量的 13.07%，占资源减少总量的 87.00%，是黑龙潭断流概率和跨度增加的主要原因。20 世纪 90 年代21 世纪 10 年代，随着植被的恢复，蒸发排泄量相比20 世纪 80 年代持续减少 1830.34 万 m3，有利于黑龙潭涵养水源恢复，但人工开采、地表河系渗漏、城市开发等因素造成地下水资源持续减少 5058.69 万 m3，表 219532017 年黑龙潭地下水系统补给量和排泄量 10 年平均值变化统计表Table2Statisticsofthe10-yearaveragevaluechangeofreplenishmentanddischargeofthegroundwaterintheHeilongtanspringfrom1952to2017年代20世纪21世纪多年平均5360617071808190910001101117降水入渗补给量/万m38129.738494.858238.228214.888923.728615.186998.368230.71黑龙潭泉群水量/万m34448.274544.824081.523633.944298.413854.67239.223585.84占补给量/%54.7253.5049.5444.2448.1744.743.4242.62侧向排泄量水量/万m33412.473668.963884.134309.134330.064475.466533.484373.38占补给量/%41.9843.1947.1552.4648.5251.9593.3654.09其它排泄水量/万m3268.99281.07272.57271.81295.25285.05225.66271.49占补给量/%3.313.313.313.313.313.313.223.301188中国岩溶2023年总体表现为不利于黑龙潭地下水的恢复，是造成黑龙潭断流并迅速恶化的主要原因，其中，人工开采地下水和城市开发造成资源量减少 4580.52 万 m3，占地下水系统多年平均补给量的 55.30%，占减少总量的 90.54%（表 3）。表 3黑龙潭及丽江城地下水系统的水资源量变化统计表Table3StatisticsofthequantitychangeingroundwaterresourceoftheHeilongtanspringandLijiangCity年代变化变化因素降水（浮动）蒸发灌溉水人工开采地表河系渗漏盆地入渗面积变化城市化其它影响合计19511990变化量/万m35991082.7450.0537.00074.7501244.54等价降水量/mm1607.55.5011018419912018变化量/万m3777+1830.34243.631300.70234.54279.823000.003228.35等价降水量/mm+2703618736414594894.2可信度分析根据丽江井水位监测资料和丽江城地下水系统水平衡方程式计算结果，整体看人类活动对地下水干扰相对较轻的 19922003 年间，年平均水位上升时补给量大于排泄量，年平均水位下降时，补给量小于排泄量；年平均水位变化较小时，补给量和排泄量基本平衡。如：1994 年 1 月和 2004 年 1 月，地下水水位分别为正 2.7m 和正 2.8m，水位上升 0.1m。水均衡计算结果显示，黑龙潭地下水系统和丽江城区地下水系统的补给量平均值比排泄量多 78.38 万 m3，仅占平均补给量的 1.38%，基本平衡。根据水平衡方程式，19902004 年，九子海多年平均降水量 1318mm，类比计算黑龙潭多年平均排泄量应为 4505.06 万 m3a1，略小于实际监测结果4896.00 万 m3a1，约为 389.94 万 m3a1。根据计算结果反推，1950 年之前黑龙潭断流的降水量阈值约为 960mm，19512018 年间仅 3 年降水量小于该值，平均 22.7 年出现一次，略高于 17291960 年间平均约 25.7 年断流一次的概率，符合自然断流规律。综上所述，丽江城和黑龙潭地下水系统水平衡方程式符合实际情况，计算因子和参数选择合理，计算结果可信。5地下河水均衡的变化趋势预测评价19891992 年和20062010 年年均降水量均在1300mm 左右，属丰水年，黑龙潭泉群由不断流变为频繁断流，地下河深部排泄量由年均 3438.51万 m3增加到 6438.89 万 m3，年均增加 166.69万 m3（表 4）。表 4不同年份黑龙潭地下水系统排泄量变化一览表Table4DischargechangesoftheHeilongtangroundwatersystemindifferentyears时间节点/年份水位/m(井口为0m)降水量/mma1入渗补给量/万m3a1黑龙潭排泄量/万m3a1其它泉水排泄量/万m3a1侧向排泄量/万m3a1断流次数/次19891992无1318.48913.215179.79294.913438.510199319971.040.941237.98369.372657.01276.925435.442200620101.751.571300.38790.892061.14290.866438.894以此推算，盆地内地下水环境若无明显改善，20212030 年，丰水年黑龙潭将成为季节性泉水，平水年和枯水年黑龙潭跨年断流；2030 年之后，系统侧向排泄量将超过多年平均降水入渗补给量，黑龙潭泉群将长期处于断流状态，再现的可能性小。6结论（1）黑龙潭地下水系统属典型的“自产二流型”岩溶水系统，泉水和侧向排泄是该系统的主要排泄第42卷第6期任世川等：丽江黑龙潭地下河系统的水均衡分析与预测评价1189方式，二者之和约占总排泄量的 96.7%，黑龙潭泉群流量受侧向排泄量影响显著。1949 年后，黑龙潭泉群排泄量由缓慢下降逐渐转变成迅速下降，20 世纪5080 年代，黑龙潭泉群排泄量持续减少了 19.2%，90 年代因降水丰富回升了 7.2%，21 世纪 0010 年代，黑龙潭排泄量与 90 年代相比减少了 92.9%；侧向排泄量变化趋势与黑龙潭排泄量呈相反的趋势。（2）降水量的年际变化是黑龙潭断流的主要控制因素，19511990 年间人类破坏植被引起蒸发量的增加是其断流加剧的主要因素，19912018 年人类工程活动（如地下水人工开采、城市化进程的影响等）迅速增强，逐渐成为黑龙潭断流频繁发生并不断加剧的主要因素。（3）盆地内地下水环境若无明显改善，20212030 年，丰水年黑龙潭将成为季节性泉水，平水年和枯水年黑龙潭跨年断流；2030 年之后，系统侧向排泄量将超过多年平均降水入渗补给量，黑龙潭泉群将长期处于断流状态，再现的可能性小。（4）科学规避人类工程活动对地下水环境的改造和影响，最大限度地减小地下水人工开采量和降低城市化进程的不利影响，是避免黑龙潭断流进一步加剧的必然之路。参考文献Xing Liting,Huang Linxian,Chi Guangyao,Yang Lizhi,LiChangsuo,HouXinyu.Adynamicstudyofakarstspringbasedonwaveletanalysisandthemann-kendalltrendtestJ.Water,2018,10（6）：698.1王宇,何绕生,刘海峰,王梓溦,晏祥省,双灵,彭淑惠.昆明翠湖九龙池泉群断流原因及恢复措施J.中国岩溶,2014,33（3）：263-271.WANGYu,HERaosheng,LIUHaifeng,WANGZiwei,YANXiangsheng,SHUANGLing,PENGShuhui.DroughtcausesandrestorationmeasuresforJiulongchispringgroupwithinCuihulakeJ.CarsologicaSinica,2014,33（3）：263-271.2梁永平,王维泰,赵春红,王玮,唐春雷.中国北方岩溶水变化特征及其环境问题J.中国岩溶,2013,32（1）：34-42.LIANGYongping,WANGWeitai,ZHAOChunhong,WANGWei,TANGChunlei.Variationsofkarstwaterandenvironmen-talproblemsinNorthChinaJ.CarsologicaSinica,2013,32（1）：34-42.3于苗,邢立亭,吴吉春,侯玉松,朱恒华,彭玉明,李罡.基于时间序列分形的济南岩溶大泉动态研究J.地质学报,2020,94（8）：2509-2519.YU Miao,XING Liting,WU Jichun,HOU Yusong,ZHUHenghua,PENGYuming,LIGang.Studyoflargekarstsprings4usingthetimeseriesfractalmethodinJinanJ.ActaGeologicaSinica,2020,94（8）：2509-2519.韩宝平.国外岩溶水文地质学进展J.中国岩溶,1993,12（4）：400-406.HANBaoping.AdvancesinforeignkarsthydrogeologyJ.Car-sologicaSinica,1993,12（4）：400-406.5张军以,王腊春,苏维词,曾春芬.岩溶地区人类活动的水文效应研究现状及展望J.地理科学进展,2014,33（8）：1125-1135.ZHANGJunyi,WANGLachun,SUWeici,ZENGChunfen.Sta-tusandprospectofthehydrologicaleffectsofhumanactivitiesinthekarstareaJ.ProgressinGeography,2014,33（8）：1125-1135.6曾成,杨睿,杨明明,胡君春,武贵华,樊宇红.丽江市黑龙潭泉群断流的人工神经网络模拟J.中国岩溶,2013,32（4）：391-397.ZENGCheng,YANGRui,YANGMingming,HUJunchun,WUGuihua,FANYuhong.ArtificialneuralnetworksimulationtozeroflowoftheHeilongtanspringgroupsinLijiangJ.Carso-logicaSinica,2013,32（4）：391-397.7康晓波,王宇,张华,曹瑾.丽江黑龙潭泉群水文地质特征及断流的影响因素分析J.中国岩溶,2013,32（4）：398-403.KANGXiaobo,WANGYu,ZHANGHua,CAOJin.Hydrogeo-logicfeaturesandinfluencefactorsofzeroflowoftheHeilong-tanspringgroupinLijiangJ.CarsologicaSinica,2013,32（4）：398-403.8周友妹.丽江古城黑龙潭泉水断流原因分析J.人民珠江,2015,36（1）：54-58.9和菊芳,方金鑫.丽江黑龙潭泉水动态变化与降水关系初探J.中国农村水利水电,2019（8）：26-27,35.10李豫馨.基于时间序列分析的丽江黑龙潭泉域动态研究D.成都：成都理工大学,2016.LI Yuxin.Dynamic research of the Lijiang Heilongtan springsystembasedonthetimeseriesanalysisD.Chengdu:ChengduUniversityofTechnology,2016.11高伟.云南省丽江市黑龙潭泉域地下水系统分析D.成都：成都理工大学,2016.GAO Wei.The systematic analysis of Heilongtan spring areakarst-water system in Lijiang City of Yunnan ProvinceD.Chengdu:ChengduUniversityofTechnology,2016.12唐克旺,谭乃元,谭大书.丽江黑龙潭泉群保护理论与实践M.北京：中国水利水电出版社,2015.13中国人民解放军00931部队.区域水文地质普查报告(丽江幅)R.1979.14四川省地质调查院.石鼓、丽江县等5幅区域地质调查报告R.2016.15丽江纳西族自治县县志编纂委员会.丽江纳西族自治县志M.昆明：云南人民出版社,2001.16王瑞雪,史茂,余其芬.基于遥感影像解读丽江古城的历史变迁J.昆明理工大学学报,2000,35（3）：6-11.WANG Ruixue,SHI Mao,YU Qifen.Historical changes of171190中国岩溶2023年ancientLijiangtownbasedonremotesensingimagesJ.Jour-nalofKunmingUniversityofScienceandTechnology,2000,35（3）：6-11.方宝初.丽江地区森林历史变迁及其对环境的影响J.林业调查规划,2002,27（4）：42-45.FANGBaochu.HistoricaltransitionoftheforestsofLijiangpre-fectureanditseffectonenvironmentJ.ForestInventoryandPlanning,2002,27（4）：42-45.18蒋雪中.滇西北漾弓江流域近期景观格局演变分析J.华东师范大学学报(自然科学版),2005（3）：84-91.JIANG Xuezhong.The evolution analysis of the spatial land-scapepatternintheYanggongriverwatershed,northwestYun-nanProvinceduringlast15yearsJ.JournalofEastChinaNor-malUniversity(NaturalSciences),2005（3）：84-91.19杨海潮.丽江古城及形象的形成与变迁J.西南民族大学学报(人文社会科学版),2009,219（11）：44-51.YANGHaichao.FormationandtransformationoftheLijiangoldtownanditsimagesJ.JournalofSouthwestMinzuUniversity(HumanitiesandSocialScience),2009,219（11）：44-51.20杜建括.基于遥感与实测的玉龙雪山冰川变化D.兰州：兰州大学,2008.DUJiankuo.TheresearchofglacierchangeintheYulongSnowMountains based on remote sensing and observation dataD.Lanzhou:LanzhouUniversity,2008.21董智杰,施昆,王瑞雪,隋玉成.基于遥感影像解读玉龙雪山的22冰雪动态变化J.科学技术与工程,2009,11（9）：2075-2093.DONGZhijie,SHIKun,WANGRuixue,SUIYucheng.Inter-pretedthedynamicchangesoftheiceonYulongJokulbasedonremotesensingimagesJ.ScienceTechnologyandEngineering,2009,11（9）：2075-2093.李敏敏,张增祥,赵晓丽,汪潇,刘芳.丽江市空间扩展与外部形态演变J.遥感信息,2015,30（6）：102-110.LIMinmin,ZHANGZengxiang,ZHAOXiaoli,WANGXiao,LIU Fang.Urban spatial expansion and external morphologyevolutionofLijiangCityJ.RemoteSensingInformation,2015,30（6）：102-110.23郑香.玉龙雪山冰川变化研究J.社会研究,2018：279-281.ZHENG Xiang.A study of glacier changes in Yulong SnowMountainJ.SocialResearch,2018：279-281.24李芸.浅析丽江盆地地下水开发利用现状与保护对策J.水土保持,2016,4（1）：7-11.LIYun.Briefanalys