1. 引言

“高位远程滑坡碎屑流”是一种滑动距离长，滑后冲击力大的滑坡类型。高位远程滑坡碎屑流一般是由大型岩质滑坡或崩塌在源区失稳后经过强烈的碎屑化作用而形成的超高速远程运动的岩屑流。对于高位远程滑坡碎屑流的定义，国际上一般认为符合以下四种特征的滑坡都可定义为高位远程滑坡碎屑流：一是体积巨大(通常大于10⁶m³)；二是高速运动(运动速度一般大于20 m/s，多为40~80 m/s，最高可超过200 m/s)；三是岩屑流由斜坡基岩失稳后形成(即始于岩质崩塌和滑坡，运动过程中不断破碎解体，形成碎屑流，具有流态化特征)；四是具有超强的流动性(Φ ˂ 0.33，Φ即为滑坡源区最高点至堆积区最前缘连线的正切值)。由此可见，相较于其他的地质灾害事件，高位远程滑坡碎屑流虽然发生的概率更低，但是造成的危害却尤为突出。例如：2008年5·12汶川地震导致三溪村大字岩平台上部岩土体开裂形成贯通性裂缝，在强降雨的影响下，大量水体渗入到裂缝中，在动水压力和静水压力的作用下，山体结构破坏失稳，形成高位远程滑坡碎屑流，导致运动堆积区内11户居民建筑被毁、109人失踪。为此分析高位远程滑坡碎屑流的研究现状及未来发展趋势，对于高位远程滑坡碎屑流的研究及防范治理工作显得十分必要。

目前，国内相关学者在“高位远程滑坡碎屑流”领域发表了一系列的相关性文章，如林棋文，程谦恭等 [1] 系统地总结高位远程滑坡碎屑流碎屑化作用的沉积学特征，简要概述了碎屑化沉积结构对于滑坡动力学的指示，归纳并阐述了高位远程滑坡碎屑流动力破碎耗能和滑坡减阻这一争议性问题。陈果等 [2] 以三溪村滑坡为例，利用三维离散元软件PFC，结合无侧限单轴压缩试验进行参数标定，模拟三溪村滑坡启动、下滑、堆积和冲击山区房屋的整个过程，该实验成功模拟出了三溪村高位远程滑坡碎屑流的四个阶段，建立的冲击力–时程演化关系、峰值冲击力–运程分布关系和峰值冲击力–速度的关系，可为山区建筑防灾减灾治理工程提供参考。殷跃平等 [3] 详细地记录了贵州关岭大寨高位远程滑坡碎屑流发生的地质环境条件、滑坡–碎屑流堆积特征、滑坡–碎屑流运动特征以及滑坡–碎屑流形成条件，通过分析资料得出发生高位远程滑坡碎屑流的充要条件。程谦恭等 [4] 集百家之长，厘定了高位远程崩滑的定义，通过分析和总结国内外发生的重大高位远程滑坡碎屑流事件，详细地阐述了高位远程滑坡碎屑流灾害事件的几何学、运动学特征。由此概述了迄今为止所提出的高位远程崩滑形成的各种动力学机理，分析与讨论了高位远程崩滑动力学机理研究的发展趋势与动向，提出了今后进一步研究的总体思路、目标与应重点解决的关键问题。

自2006年CiteSpace引入国内后，大量的学者便利用该工具对其所研究的领域进行可视化分析。如张攀 [5] 基于CiteSpace工具对有关少儿科普的文献进行可视化分析，探析关于少儿科普研究领域的主要问题、研究趋势以及最新研究热点，以期为未来少儿科普的研究和发展提供参考。再比如贺义雄等 [6] 借助CiteSpace知识图谱工具，对关键词进行共线和聚类分析、研究作者分析、得到了相应的知识图谱，揭示了海洋碳汇研究进程。为针对性了解高位远程滑坡碎屑流领域的研究热点及发展趋势，本文基于CiteSpace软件进行可视化分析研究，对我国的高位远程滑坡碎屑流领域进行了细致的图谱分析。

2. 数据来源和分析方法

2.1. 数据来源

为了深度研究“高位远程滑坡碎屑流”的热点问题，文章对这一领域作者，研究机构和研究人员的合作关系进行了分析，并在此基础上对这一领域的发展趋势进行了研究。本文数据来源以中国知网数据库(CNKI)为基础，主题词以“高位远程滑坡碎屑流”进行检索，时间线为2006年~2022年，相关文献有241条。

2.2. 分析方法

2.2.1. 分析方向

通过CiteSpace对所得到的241条有关文献进行可视化分析，241条文献数据均来源于中国知网，通过主题词“高位远程滑坡”进行检索所得到。考虑到CNKI中推导的文献数据仅能分析作者与机构之间的合作网络关系以及完关键词共现图谱的绘制，所以本文主要是对上述知识图谱的展开分析来揭示学科内部作者与机构的关联，并且根据关键词共现程度，引出“高位远程滑坡碎屑流”这一领域发展趋势和其热点研究方向。

2.2.2. CiteSpace软件分析

CiteSpace是一款基于Java开发的文献计量软件，以科学知识为计量单位，用数学方程式方法来表达科学发展规律，基于共引分析理论和寻径算法网络对特定文域(集合)进行计量，分析学科领域演化的关键路径和学科拐点 [7]。目前，CiteSpace主要用于对研究领域的热点话题和发展方向进行可视化分析，以期对研究者跟踪这一领域的最新热点及今后的发展趋势有所帮助。本文基于CiteSpace软件对高位远程滑坡进行可视化分析，通过分析所得到的作者发文量知识图谱、发文机构知识图谱、关键词知识图谱、聚类知识图谱、时间线聚类图谱，进而得到近年来我国关于高位远程滑坡领域的研究热点和发展趋势。

2.2.3. LLR对数似然算法

LLR对数似然算法可以得出某个聚类的紧密程度。Ochiia相似系数能表现出文本之间的共现率

$\text{Cos}\left(A,B\right)=\frac{\left|A\cap B\right|}{\sqrt{\left|A\right|\left|B\right|}}\left(A\ge 0,B\ge 0\right)$

A，B分别代表关键词的出现频次， $A\cap B$代表关键词的共现频率，所以依据三角函数定理，当Cos(A, B) = 0时，A，B之间关联度为0，当Cos(A, B) = 1时，A，B之间关联度为最大 [8]。

3. 结果分析

3.1. 文献历年发文量分析

文章基于中国知网数据库(CNKI)，以“高位远程滑坡碎屑流”为主题词进行文献检索，检索的时间区间定为2006年1月到2022年8月。通过检索所得到的相关文献共计241篇，以此为对象展开本次的研究工作。图1给出了2016年~2022年每年有关高位远程滑坡碎屑流的发文数量。

由图1可知，从2008年开始，国内的发文数量开始增加，高位远程滑坡碎屑流问题逐步得到国内研究人员的重视。其中2016年发文量最多，达到了29篇。从2012年至2022年有关高位远程滑坡碎屑流的发文数量都在11~29篇以内。虽然最近五年的发文数量都较为稳定，但较之于2008年以前而言，我国在“高位远程滑坡碎屑流”领域的发文数量已有大幅增加。这说明在近些年以来，我国有关学者一直密切关注有关高位远程滑坡碎屑流方面的研究。

3.2. 文献作者群体分析

在由CiteSpace处理后得到的作者关系图谱中，作者文献发表量愈大，所代表节点也愈大，图中连线反映的是作者群体间的合作。在本次研究所得到的“高位远程滑坡碎屑流”作者关系图谱中，节点数有245个，连接438个，网络密度为0.0147。在高位远程滑坡碎屑流学科领域中，作者间的协作是围绕几大学者进行的，如殷跃平，程谦恭，王玉峰等学者，他们与其他作者之间具有较多的合作关系。但是也有一些学者发文量不多，发文量只有1~3份，彼此之间合作不多，说明研究高位远程滑坡碎屑流领域的作者较聚集，主要集中在少数学者中。在高位远程滑坡碎屑流研究领域中，核心作者分别构成了一个密切合作的联系网络，相互之间衔接紧凑，构成了两足鼎立之势，一部分是由殷跃平为核心形成的小合作网络，另一部分则是以程谦恭和王玉峰为核心而形成的合作关系网络。经图2、表1分析可知，在我国的“高位远程滑坡碎屑流”领域中，该领域的研究人员往往是以团队的形式存在，研究人员联系都较为密切。

3.3. 机构合作分析

对检索得到的241条文献进行可视化机构网络分析，CiteSpace中相关的设置如下：选取的节点类“NodeTypes = Institution”，检索区间“Time slicing”取2006~2022年，时间切片“Year Per Slice = 1”，阈值取Top = 50。由此可所得到的高位远程滑坡碎屑流发文机构图谱如图3所示，并以频次为先后顺序进行排名可得出表2。

经CiteSpace处理后所得到的发文机构知识图谱中，共有130个节点，151条连接，网络密度为0.018，说明了在高位远程滑坡碎屑流方面，机构间交流与合作关系较密切，这有助于高位远程滑坡碎屑流研究领域的持久开展。具体来看，图中出现了以中国地质调查局为核心，其中包含上海交通大学土木工程学院、长安大学地质工程与测绘学院等高校以及中国地质环境监测院和国土资源部新构造运动与地质灾害重点实验室等国家重点实验室。值得注意的是，在以中国地质调查局为核心的圈外，西南交通大学和成都理工大学这两所高校也凭借其独特的地理优势和国家重点实验室，为高位远程滑坡碎屑流领域做出了突出贡献。再者，这些研究所还与各大院校有合作，实现了研究资源共享，共同克服了关于高位远程滑坡碎屑流问题。相信在后续的发展研究过程中，各机构更能保持密切交流合作。

3.4. 关键词分析

3.4.1. 关键词共线图谱分析

一般一篇文献的关键词之间常有联系，对关键词进行共现分析，能够反应出学科领域中的重要研究领域，及逐年学科演化和发展，也能直观地反映出不同时序下热点领域、分析视角和研究方法发生了转变 [9]。首先，在Citespace的设置中将时间切片定为“Year Per Slice = 1”，阈值定为Top N = 50，得出关键词知识网络图谱，在该图谱中网络线的颜色体现了关键词首次被引的时间，因此从整体上分析网络线的颜色变化就可以清楚地知道该研究领域的新旧情况。如图4所示，为2006~2022年高位远程滑坡碎屑流关键词图谱图。

结合图4可知，2006~2022年高位远程滑坡碎屑流领域研究，出现次数最多的是数值模拟，一共出现了36次。其次碎屑流，汶川地震、动力破碎出现的频率也较为频繁。据此，我们可以分析出，过去十多年来，我国“高位远程滑坡碎屑流”研究的许多方向是围绕图4所示4个方向展开的。通过量化关键词可以使我们的分析更为准确，由此进行统计分析可得到表3。

根据表3可知，“数值模拟”出现频次最高，达36次，其次“碎屑流”(28次)，“汶川地震”(19次)，“动力破碎”(18次)，“滑坡”(16次)，“高速远程”(16次)等，以上的关键词都表明了高位远程滑坡碎屑流领域的重要研究热点和发展趋势。可以发现，最近几年来我国的研究人员多从有限元数值模拟的角度去对高位远程滑坡领域进行研究探索，例如张艳玲团队 [10]，以哀牢山地区芭蕉树滑坡为研究对象，根据芭蕉树滑坡的基本特征，利用DAN-W软件建立出了芭蕉树滑坡的计算模型。该模型能够较为准确地呈现芭蕉树滑坡的整个动力模拟过程，并根据滑坡特征，结合滑坡分区和滑坡运动过程中的铲刮效应，将滑坡滑动路径分为了滑源段、铲刮段、堆积段3段。该项目的数值模拟所得的结果与调查与所得的滑坡滑移结果吻合度高。由此可见，对于“高位远程滑坡碎屑流”领域数值模拟这个环节是必不可少的。

3.4.2. 关键词聚类LLR算法分析

通过基于LLR对数似然算法(Log-Likelihood Ratio)的聚类分析，可以推测出高位远程滑坡碎屑流领域的研究热点。本文利用CiteSpace的快速聚类方法，提取出关键词，再通过阈值调整，得到图5所示的知识图谱。该图中模块值(Modularity)为Q = 0.6914 > 0.3，说明聚类结果显著，平均轮廓值(Mean)为S = 0.8807 > 0.7，说明聚类结果的可信度是非常高的。由此可得，2006~2022年高位远程滑坡碎屑流领域的研究主要是从图5所示的9大聚类模块进行展开。从平均年份可以看出该领域的发展初期在2011年左右，起步较晚，主要原因是2008年以前我国鲜有发生这类大规模且速度极快的滑坡现象，高位远程滑坡碎屑流在当时并未得到国内研究人员的重视。但是自从2008年汶川地震之后，高位远程滑坡碎屑流地质灾害事件频繁出现，国内学者开始对高位远程滑坡碎屑流逐渐加大研究力度。

表4的紧密程度体现了各聚类中的所有关键词同质性，若紧密程度值越大，那么，代表这个聚类成员相似性就越大。在给出的表4中每个聚类板块紧密程度都超过了0.8，这说明本次聚类效果很好，关键词之间联系紧密。如在“地质灾害”聚类中，“滑坡体”、“减灾防灾”和“地址选线”是最紧密的三个关键词，在“颗粒流”聚类中，“模型试验”、“细观流态”、“科里奥利效应”联系最为最为紧密。

CiteSpace所得到的时间线图(Time-line)将聚类分析和时间线结合到一起，在时间线图中同一聚类的节点按照时间顺序会被排布在同一水平线上，同时不同聚类之间的联系也可以清晰展现在图中，通过时间线图可以清晰直观地了解到聚类的时间跨度以及历史进程。从经过CiteSpace可视化处理后得到的图6中我们可以看出，以上的九个节点分别对应了各聚类的历史发展历程，通过这些时间节点我们可以做以下分析：1) 聚类#0碎屑流：研究人员先是从碎屑流成因的原理出发进行研究，随后采用能量法分析，接着建立各种关于碎屑流的模型。例如郜颖超，位伟，姜清辉 [11] 采用数值流形方法，对牛圈沟强震诱发滑坡所形成的高速碎屑流运动过程进行模拟，分析了其堆积状态及运动特性。在建立碎屑流模型的过程中，研究人员也从碎屑流的各个方面进行研究。2) 聚类#6动力破碎：将理论和实验两者结合，以此来研究动力破碎对高位远程滑坡碎屑流的影响。在这个过程中研究人员进行了“模型试验”与“斜槽实验”中。如郑光 [12] 基于现场典型灾害实例调查和已发生碎屑流数据统计，确定了影响远程运动距离的主要因素，从而建立了一系列碎屑流体的斜槽堆积实验。

3.4.3. 主题演进与前沿分析

利用Burst detection功能我们可以获得某点时间内突然大量使用的关键词。在进行关键词突现分析时，在其余设置均为不变的情况下，将突现词的最小持续时间定为2年，由此得到突现词2个。表5对突现关键词的强度、出现年份、结束年份、持续时间进行整理，从表5可以看出关于关键词“汶川地震”的强度较高。

强度最高的突现词为“汶川地震”，起始于2008年，于2011年强度开始降低。2008年，汶川地震引发了多起高位远程滑坡碎屑流灾害事件，分析其诱发原因、运动特征以及防范措施在当时是迫在眉睫的任务。因此从汶川地震发生以后，国内便掀起了对高位远程滑坡碎屑流的研究的热潮。例如殷跃平 [13] 通过解剖牛圈沟滑坡–碎屑流、北川城西滑坡和青川东河口滑坡–碎屑流3个典型实例，分析出汶川八级地震远程高速滑坡具有的五个特征——岩性条件、抛掷效应、碰撞效应、铲刮效应、气垫效应，此外也指出了汶川地震高速滑坡初始滑动主要来自于强震的抛掷，从而确定了高速远程的初始速度要求是来源于强震的抛掷而并非平常所认为的来自于滑带强度的突然丧失。吴鑫培 [14] 则是在系统研究高位远程滑坡碎屑流运动机理的基础上，进一步探讨研究区内形成的地质灾害链的成灾模式及形成条件。通过建立PFC3D模型，研究了灾害链链生机理的方式，将整条文家沟灾害链分成两个阶段来进行分析。这为工程建设和防灾减灾提供科学的技术支持，具有重要的理论和实践意义。以上案例都是基于汶川地震所引发的灾害事件而开展的有关高位远程滑坡碎屑流的研究。

除此之外，“数值模拟”也是高位远程滑坡碎屑流研究领域绕不开的话题。关于数值模拟软件的使用，苏生瑞 [15] 等采用2D-Block软件，建立了二维离散元数值模拟模型，研究了相应地质体在不同阶段下的运动特征。齐超等 [16] 采用了滑坡动力分析软件DAN-W分别建立了摩擦模型、Voellmy模型和F-V等3种不同的滑坡数值模型，对东河口滑坡–碎屑流的运动距离、不同时刻运动速度特性、堆积物分布规律及滑坡堆积体积进行了模拟，同时对滑坡运动时间进行了估算。

4. 结论与展望

本次研究以CNKI数据库为平台，运用CiteSpace分析2006~2022年间高位远程滑坡碎屑流研究领域相关文献知识图谱结构，从文献的历年发文量，文献的作者群体，机构合作和关键词等作了详细说明，具体说来，可以得到如下三点结论。

1) 目前，我国在高位远程滑坡碎屑流研究领域较2008年前取得了长足的进步，不仅是发文数量上有上升，学科有关的研究人员的人数也呈现上升趋势。在此基础上，与早期单方面的研究相比，2008年后聚类方向较多，高位远程滑坡碎屑流学科进一步发展和提高。从发文机构方面得知，西南交通大学作为该领域的核心研究机构其发文量远高于其他机构。从机构分布地区区域来看，在西南片区，西南交通大学和成都理工大学双方联系最为密切。

2) 根据论文发表数量趋势来看，接下来的几年关于高位远程滑坡碎屑流的发文数量仍会持续增加，将会有越来越多的学者关注这个领域。但是在该领域中研究结构的联系合作较少，望在后续的发展研究过程中，各机构能保持密切交流合作。希望作者之间能加大交流合作力度，携手攻克高位远程滑坡碎屑流的问题。

3) 通过所得到的关键词共现图谱来看，“数值模拟”、“碎屑流”、“动力破碎”等为近几年的研究热点，就此说明了目前高位远程滑坡碎屑流领域的前沿工作是有“数值模拟”和“碎屑流”，因此我们在研究该领域的过程中，应当重点把握这两个方面的学习。

基金项目

中国水利水电第七工程局有限公司2021年项目赞助。

NOTES

^*通讯作者。

参考文献