我国在国家战略层面上提出建设数字孪生流域以推动水利高质量发展
[1-3],构建数字孪生流域需要对水循环要素(径流量)进行精准化映射
[4-5]。但受环境、成本和人员的制约,单纯依靠地面监测站已无法满足径流量数据的获取,因此如何精确反演径流量是当下建设数字孪生流域的重点与难点。
近年来,卫星遥感数据的空间、时间分辨率不断提高,成为地面监测的重要补充
[6-7]。利用遥感技术获取河流水文、水力要素(如河宽、水位、水面面积、流速等)
[8]并反演径流量,已成为当前研究的前沿方向。目前,虽有学者就径流量反演的方法做出了定性研究
[9-10],但仍然缺少基于知识图谱等可视化计量工具进行的定量研究。
某一研究领域的发展态势和现状可以通过文献计量法进行分析
[11]。文献计量法是通过数学和统计学方法对文献进行量化分析,以揭示研究趋势、合作关系、影响力和研究热点,常通过CiteSpace软件进行分析,该软件能够以知识图谱的形式展现某一科学领域中的研究现状和前沿方向
[12-14]。此软件已广泛应用于诸多研究领域,如张馨等
[15]通过CiteSpace对萝芙木属植物在中医药领域的研究现状及前沿方向开展了可视化分析;郭亚丽等
[16]运用CiteSpace对温室气体与空气污染物协同减排领域的发展趋势与现状进行可了视化研究;焦红波等
[17]通过CiteSpce对水库移民后期扶持的研究脉络与现状进行了可视化分析。因此,基于文献计量学方法对径流量反演领域进行可视化分析具有重要意义
[18]。
基于此,本研究拟采用文献计量工具CiteSpace软件,对中国知网(CNKI)与Web of Science(WOS)数据库中2000—2023年径流量反演相关文献的发文量、发文机构、发文国家、共被引作者、共被引期刊和关键词6个方面进行量化分析,旨在揭示该领域研究趋势及研究现状,为进一步的理论和实践研究提供参考。
1 数据来源与研究方法
1.1 数据来源
本研究中的“径流量遥感反演”采用广义定义,泛指基于遥感技术与物理约束和统计推断进行估算流量的过程。中、英文文献分别基于CNKI和WOS数据库进行检索,检索时间范围均设置为2000—2023年,收集相关文献。在CNKI数据库中,以“遥感径流”“径流量反演”为单独词,“遥感数据”和“径流量”为组合词进行关键词检索,共得到886篇相关文章,使用EndNote软件剔除重复、相关度较低的文章后,最终得到343篇文献。在WOS数据库中,以“discharge retrieval”为单独词,“remote sensing”和“river discharge”为组合词进行关键词检索,共得到1 782篇文章,同样使用EndNote软件进行筛选、去重,最终得到1 383篇文献。
1.2 研究方法
基于筛选去重后两大数据库的文献数据,在CiteSpace 6.3.1软件中设置时间切片为1年,分别对发文量、发文机构、发文国家、作者被引、期刊被引和关键词6个方面进行可视化分析。其中,对发文量进行时间趋势分析并结合Mann-Kendall(MK)突变检验法计算具体突变年;对发文机构和发文国家基于普莱斯定律分析核心发文机构与国家;对共被引作者和共被引期刊进行中心性分析;对关键词运用Log-likelihood ratio(LLR)算法进行聚类与共现分析。其中,普莱斯定律相关计算公式如下
[19-21]:
核心作者和核心机构的最低发文量(n)计算公式为:
式中:max为最高产作者或机构的发文量。检验核心作者和核心机构是否形成核心群,计算公式为:
式中:
M为普莱斯指数结果占比,
M≥50%即表示已形成核心群
[22];
Nmax为高于普莱斯定律数值频数;
N为总词频数。
2 研究现状分析
2.1 发文量趋势
发文量是体现研究领域关注度的重要指标,一定程度上能够反映该领域的发展动态
[23]。CNKI和WOS数据库的发文量结果见
图1。由
图1(a
1)和(b
1)可知,CNKI数据库中径流量遥感领域文献数量近年呈增加趋势,但整体较WOS数据库波动更大、增速更低。利用MK突变检验法得到UF和UB曲线相交于2015年,此时纵坐标统计量位于99%置信区间外,说明2015年为突变年。按照文献发表量可分为2个时期:1)2000—2015年为低速稳定期,年发文量从2篇增至13篇,线性拟合斜率为0.97;2)2016—2023年为高速增长期,年发文量从23篇增至43篇,线性拟合斜率为3.37。2016年以来,CNKI数据库发文量出现激增,这与国产高分系列卫星(如高分一号、高分二号等)发射并投入使用以及《全国水文监测规划(2016—2020年)》的发布有关。高分卫星具有高空间分辨率、高时间分辨率和多光谱等特性,能够更精准地获取地表信息,为径流量监测提供了更丰富的数据支持;《全国水文监测规划(2016—2020年)》明确了“十三五”期间要通过技术手段创新提升监测能力,其中重要一项为利用遥感技术进行水文监测,进一步激发了国内学者对相关领域的探索热情,在政策因素与数据因素的共同影响下,使得2015年成为CNKI数据库遥感监测径流量领域发文量的突变年。以拟合线性趋势预测,未来国内径流量遥感反演研究文献数量将会持续升高。
由
图1(a
2)和(b
2)可知,WOS数据库文献量较CNKI数据库呈现更加明显增加趋势。利用MK突变检验法得到UF和UB曲线相交于2012年,此时纵坐标统计量位于99%置信区间外,说明2012年为突变年。按照文献发文量可分为2个时期:1)2000—2012年为低速稳定期,年发文量从7篇增至48篇,线性拟合斜率为3.36;2)2013—2023年为高速增长期,年发文量从62篇增至171篇,线性拟合斜率为11.66。2013年以来,文献激增与2013年国际水文计划(IHP)政府间理事会第二十一届会议召开有关,21世纪以来随着遥感技术的发展,水文领域也逐渐对遥感技术重视起来,此次会议明确了国际水文计划第八阶段(2014—2021年)的战略计划,其中一项设定为遥感技术在水文信息收集、水安全和水灾害等领域的应用,自此国际英文相关领域的发文量逐年迅速升高。以拟合线性趋势预测,未来国际径流量遥感反演研究文献数量将会保持高增长速率。
2.2 发文机构及国家
研究机构与国家的发文情况在一定程度上反映了研究团体分布,有助于学者了解径流量遥感领域合作关系。由于CNKI数据库相关文献量较少,合作关系不明显,无法进行可视化分析,因此,仅针对WOS数据库进行分析。2000—2023年WOS数据库径流遥感领域发文机构与发文国家知识图谱见
图2,发文量排名前10的机构与国家见
表1。
WOS数据库共得到593个机构,1 185次合作关系,网络密度为0.006 8,根据普莱斯定律,核心机构共有38个(
n=8),共发表文献712篇,占文献总数的51.48%(>50%),表明WOS数据库已形成核心机构群。发文量排名前3的机构依次为Chinese Academy of Sciences(中国科学院,121篇)、University of Chinese Academy of Sciences(中国科学院大学,46篇)、California Institute of Technology(加利福尼亚理工学院,35篇)。其中,中国科学院节点颜色从白色至红色均有分布,说明其在2000—2023年每年持续发文。由
表1可知,中国科学院发文数量远超其他机构,仅其中心性>0.1,表明其在研究领域内具有重要中心性地位,与其他机构联系紧密。总体来看,其他机构的发文数量相近,发文时间均衡,但中心性均未超过0.1,不足以构成中心重要节点。
发文量前3位的国家依次为USA(美国,421篇)、China(中国,311篇)和France(法国,112篇)。前2位国家的发文量总和(732篇)远超其他8位国家的总和(501篇),占前10位国家发文量总量(1 233篇)的59.37%,说明美国和中国的径流量反演研究发文量居于世界前列,为核心发文国家。根据中心性可知,前10位国家的中心性均>0.1,与其他国家联系最紧密的为美国(0.48),其次是中国(0.37),法国(0.31)位居第3。由
图2可知,前10位国家每个节点从内至外的颜色由白色渐变为红色,说明在每个时间段均有研究成果发表,较明显的是,以2014年的紫色圆环为分界,美国紫色圆环半径与其总圆环半径的比值要高于中国和法国,且中国和法国在2000—2023年的发文总量与美国接近,说明在2014年以前国际上发文主要以美国为主,2014年以来,中国和法国发文量占据更重要地位。
2.3 共被引作者及共被引期刊
通过分析2000—2023年径流量反演领域作者共被引和期刊共被引,可以识别研究领域高影响力作者和期刊,但由于CNKI数据库导出数据不包含参考文献等相关信息,因此,仅针对WOS数据库进行分析。
WOS数据库共951个共被引作者,5 082次合作关系,共被引量排名前10位的作者见
表2。中心性值显示,前10位作者中心性最高为0.08<0.1,无重要中心节点。Bjerklie D M和Smith L C共被引次数较高,2位作者曾合作发表多篇高质量论文
[24-26]。
WOS数据库共794个共被引期刊,6 540次合作关系。Remote Sensing of Environment共被引次数最多(856次),比第10位多470次,表明该期刊在径流量遥感领域发表的文章与其他期刊相比更具影响力。中心性结果显示,前10位期刊中心性均<0.1,其中Science中心性最高为0.06,具有较重要影响力。
3 研究前沿与热点
3.1 研究热点
分别对CNKI和WOS数据库文献进行关键词共现分析得到的重点关键词见
表3,合并中、英文文献,通过关键词共现分析得到以下6个主要研究方向:
1)利用多源遥感、深度学习等方法获取更高精度的水文、水力要素数据。径流量反演需要长时间序列的水文、水力要素数据。结合多源遥感技术,可以获取更全面的地表信息,利用深度学习(如卷积神经网络、生成对抗网络等)对遥感数据进行处理和分析,也可以提高水力要素的提取精度
[27]。
2)数据源与水文模型对不同流域的适配性研究。流域特性因地而异,而当下数据源与水文模型种类繁多,不同流域所适用的数据源与水文模型均有差异,因此当下有大量针对数据源和水文模型对流域适配性的研究
[28]。
3)以曼宁公式及其经验公式为理论基础,遥感数据为驱动的水动力学模型研究。曼宁公式广泛应用于水动力学模拟径流过程。利用遥感技术获取的数据(如水位、流速、河宽等)代入曼宁公式可以得到瞬时流量的相对真值。此外,还有学者结合实际观测数据,对曼宁公式和其他经验公式进行改进,提高结果精度
[26]。
4)通过水力几何法建立重要水文、水力要素与径流量关系模型。水文、水力要素与径流量之间存在直接或间接的定量关系,水力几何法通过建立相关数学模型,可以用于径流量的反演和预测
[29]。
5)以基于深度学习为代表的数理模型反演径流量。深度学习模型(如LSTM、CNN、Transformer等)可以捕捉复杂的非线性关系,部分研究者通过调整模型结构、优化超参数、引入新的训练算法等方法,不断提高深度学习模型的性能,对径流量实现高精度反演
[30]。
6)基于大环境变化下,研究径流量的变化趋势。气候变化和人类活动均对径流量产生显著影响,气候变化包括温度、降水和蒸发等气象要素,人类活动如农业灌溉和水利工程建设。当前也有不少学者利用长时间序列分析和统计等方法,分析径流量的长期变化趋势与原因
[31]。
方向1是方向2~5的基础,获取高精度水力要素数据是对径流量进行反演的前提;方向2~5是当前活跃在径流量遥感反演领域内的主流方法;方向6是方向2~5进行径流量反演后的延伸。对6个主要研究方向统计其发文量,分析其发展趋势,结果见
图3。6个方向发文量总体趋势均不断增加,且方向1~5均在2017年迅速增长,其中方向5增长最为明显,其在2017年之前发文量小于其他方向,但2017年后其发文量迅速增加,2019年发文量开始超过其他方向且持续高速增加,而方向1~4自2018年后发文量又趋于平稳。方向1~4于2017年发文量快速增加是因为该年前后遥感领域多个高质量卫星升空,如高分3号、吉林一号、ALOS-2和sentinel-2B相继升空,这为水力要素的高精度获取提供了更丰富、长时序的数据源。方向5于2017年后发文量持续高速增加涉及多个方面,如2016—2018年,GANs及其变种如条件GAN、InfoGAN、DCGAN等相继出现,解决了深度学习发展的一些难题,推动了无监督学习的发展;2015—2016年,Google和Facebook等相继推出TensorFlow、PyTorch和Keras等开源框架,降低了深度学习技术的使用门槛。虽然深度学习模型作为没有物理机制的“黑箱子”模型一直遭受质疑,但其热度不减的原因是该方法通过大量数据训练能够得到比传统考虑物理机制的水文模型更为准确的径流量反演结果。
3.2 关键词时间发展
对CNKI与WOS数据库进行关键词时间线分析和突现分析,结果分别见
表4和
图4。CNKI数据库分析得到Modularity
Q=0.746 1>0.3、Silhouette
S=0.928 5>0.8,聚类结构明显、效果显著。由
图4(a)可知,中文文献关键词最终得到10个大类,且不同关键词间存在着不同的时间演替,如“灾害综合风险评价”“洪涝灾害”最早出现在2000年,2008年被迭代演替,而“分类回归树”“gee”“偏差矫正” “机器学习”等关键词在2015年之后频繁出现。由
表4可知,关键词突现强度除较高的“遥感技术”(7.51)、“径流模拟”(5.81)外,其余突现强度均在2.5~4.5。结合
图4与
表4可将CNKI数据库的研究热点分为2个阶段:第一阶段为传统水文模型研究期。截至2015年,国内研究主要集中基于遥感数据获取水力要素,驱动或率定“SWAT模型”“SCS模型”和“降雨径流模型”等一系列水文模型,对不同流域径流量进行反演
[32-33]。第二阶段为技术融合快速发展期。2015年后,国内文章逐渐引进新的研究方法与思路,如“gee”“机器学习”“偏差矫正”等,这促使2015年后中文径流遥感领域文献数量呈爆发式增长,尤其在深度学习方向的文章更为突出
[34-35]。
WOS数据库分析得到Modularity
Q=0.378 3>0.3、Silhouette
S=0.734 2>0.7,聚类结构明显、效果显著。由
图4(b)可知,英文文献关键词最终得到7组聚类结果,与CNKI数据库文献聚类效果相比,WOS数据库文献聚类结果之间联系更为紧密,研究方向更全面,重要关键词出现时间更早。如20世纪初,英文文献中除了包含“model”“hydrological model”等表示水文模型的关键词外,还包括“synthetic aperture radar”“sensitivity analysis”“atmosphere correction”等代表不同数据源获取与分析方法的关键词。“atmosphere correction”在WOS数据库中2000年就已出现,而在CNKI数据库中“大气校正”直到2012年才出现;“google earth engine”在WOS数据库中2019年出现,但在CNKI数据库中“gee云平台”在2022年才出现。由
表4可知,排名前10位关键词突现强度集中在4~5。结合
图4与
表4可将WOS数据库研究热点分为2个阶段:第一阶段为多模型探索期。2000—2012年,遥感技术监测径流量处于探索阶段,主要依靠MODIS、Landsat等传感器获取数据,数据分辨率有限,研究人员尝试将遥感数据与水文模型、水动力模型和水力几何模型结合,估算流域参数来间接推算径流量,为后续发展奠定基础。第二阶段为融合发展期。2013—2023年,Sentinel和SWOT等高分辨率卫星投入使用,同时机器学习、深度学习的引入加速了领域发展,相关研究开始多方法融合,推动了径流量监测从单纯的数值计算向多元综合分析转变。
4 结论与展望
4.1 结论
本研究利用CiteSpace和Endnote软件对CNKI与WOS两大数据库中遥感反演径流量相关研究方向的共1 726篇文献进行数据处理与分析。主要得到如下结论:
1)2000—2023年,国内外径流量遥感反演相关文献在数量上呈现先低速增长再高速增长的趋势,其中英文发文量的突变点在2012年,中文发文量的突变点在2015年。
2)美国、中国和法国在该领域发文量稳居前列,中国科学院和加利福尼亚理工学院为高发文量机构,Bjerklie D M和Smith L C为高共被引作者,Remote Sensing of Environment和Journal of Hydrology为高共被引期刊。
3)关键词聚类可归纳为6个方向,方向1为方向2~5的基础,方向2~5为反演径流量具体方法,方向6是方向2~5的延伸。21世纪以来,6个方向发文量均不断提升,其中方向5(基于深度学习为代表的数理模型反演径流量)发文量于2017年开始迅速增加并于2019年远超其他方向,成为当下遥感反演径流量领域最热门的方向。
4)CNKI数据库研究热点可分为2个阶段:2000—2015年为传统水文模型探索期,研究主要围绕传统水文模型展开;2016—2023年为技术融合加速发展期,新技术的加入推动了该领域研究的快速发展。WOS数据库研究热点同样可分为2个阶段:2000—2012年为多模型探索阶段;2013—2023年为融合创新阶段,重点研究多学科方法相互交融。此外,英文文献在研究内容第一次出现的时间上普遍早于中文文献,并且在研究范围和深度上表现出更广泛的覆盖性与全面性。
4.2 展望
尽管遥感监测径流量研究领域已取得显著进展,但在一些方面仍存在诸多问题,如复杂地形处径流量反演困难、遥感数据密度无法支持长时间序列径流量反演、深度学习等方法存在模型可解释性问题、水文模型参数率定准确性问题等。针对目前面临的问题,对未来研究重点提出如下展望:
1)多源数据融合,拓展数据时空尺度。整合不同来源、类型与分辨率的数据资源,从时间尺度延长观测序列,从空间尺度拓宽覆盖范围,打破单一数据局限,提升动态变化捕捉精度。
2)改进优化模型,多模型耦合集成。深入模型机理,增强模型对复杂水文过程的模拟能力;考虑模型特点,耦合模型实现优势互补,如将深度学习与传统水文模型耦合,既保留了机理模型的物理可解释性,又发挥了数据驱动模型的性能优势。
3)依流域特性分类整理径流反演方法。依据流域的具体特性,对现有的径流量反演方法进行全面、系统地分类与整理,以实现方法与流域特性的精准匹配,筛选适合特定流域的径流反演方法,推动流域水文模拟、水资源管理与规划等相关工作。
4)气候变化与人类活动影响研究。在全球气候变化和人类活动影响日益显著的背景下,研究其对径流量的影响机制和趋势变化。利用遥感技术长期监测径流量的变化,并结合气候模型和人类活动数据,进行综合分析和预测。
国家重点研发计划(2023YFC3209302-03)
国家自然科学基金项目(52279031)
国家自然科学基金项目(52309040)
北京市自然科学基金项目(JQ21029)