高频探地雷达在古城墙修复效果评价中的应用.pdf

文章编号：1007-2993（2024）02-0196-07高频探地雷达在古城墙修复效果评价中的应用高频探地雷达在古城墙修复效果评价中的应用孟田华1黄丽1张海蛟2王浩航1卢玉和1姜子超1（1.山西大同大学，山西大同037009；2.山西省考古研究院，山西太原030001）【摘要】在长期环境污染、自然因素和人为因素的共同影响和作用下，明长城得胜堡墙体结构出现了裂隙、空鼓开裂、墙体基础掏蚀、酥碱、植物根系等病害。研究利用探地雷达技术分别对得胜堡三个代表性城墙区域进行检测，探地雷达剖面图显示修缮后的城墙体状况有明显的改观，特别是经过传统逐层夯筑法修缮后的古城墙体，每层夯筑层连续性好，无明显异常，清晰可获得均匀的夯筑层厚度为 0.02 m，可有效防治裂隙、空鼓、毛细现象及动植物等病害。采用土坯砖垒砌法修缮的城墙虽然在一定程度上起到了保护的作用，但由于土坯砖间连接的不密实，容易为毛细水、自然降水、动植物等提供产生裂隙、空鼓、生物等病害的条件，对于城墙的长久保存效果仍不够理想。采用探地雷达可快速了解古城墙的保存现状和修复性能，适宜开展大面积古城墙调查。【关键词】探地雷达；无损检测；古城墙修缮；夯土城墙【中图分类号】TU195 【文献标识码】Adoi:10.3969/j.issn.1007-2993.2024.02.012Application of High-frequency Ground Penetrating Radar in the Evaluation ofAncient City Wall Restoration EffectMeng Tianhua1Huang Li1Zhang Haijiao2Wang Haohang1Lu Yuhe1Jiang Zichao1(1.Shanxi Datong University,Datong 037009,Shanxi,China；2.Archaeology Academy of Shanxi Province,Taiyuan 030001,Shanxi,China)【Abstract】Under the combined influence and action of long-term environmental pollution,the wall structure of Deshengbao onthe Ming Great Wall has developed cracks,hollowing and cracking,erosion of the wall foundation,alkalinity,and plant root systemdiseases.In this research,the three representative wall areas of Deshengbao were detected by ground-penetrating radar technology.Theground penetrating radar profiles showed that the condition of the repaired city wall has improved,especially the repaired ancient citywall with the traditional technique of layer-by-layer ramming process.The continuity of each layer of the ramming layer is perfect,andthe thickness of the ramming layer of 0.02 m can be obtained easily,which can effectively prevent cracks,hollowing,capillarity,plantroots,and other diseases.Although the adobe brick method may have helped preserve the city wall,the long-term preservation effect ofthe wall is not the best due to the uncompacted adobe brick,which is easy to provide conditions for capillary water,natural precipita-tion,animals,and plants to produce cracks,hollowing,plant roots,and other diseases.So the method of ground penetrating radar canquickly detect the preservation status and repair performance of the ancient city wall,which is suitable for large-scale investigation ofthe ancient city wall.【Key words】ground-penetrating radar；non-destructive testing；repair of ancient city walls；rammed earth walls 0 引言得胜堡位于山西省大同市新荣区堡子湾乡，是明代长城大同镇重要关隘。得胜堡修筑于嘉靖二十八年（1548 年），先为土筑万历二年（1574 年）包砖。经历长期环境作用和人为破坏，部分墙面和墙体已出现裂缝、分层、空鼓及墙皮脱落等病害。当地文保局分别使用土坯砖垒砌法和传统的逐层夯筑法对典型区域试点进行了修缮。对文物的修复首先要满足干预最小化、修复可逆性、修复材料兼容性和修复可辨识性等基本原则，此外，材料的选择与应用将直接影 基金项目：国家自然科学基金（62071312）；省筹资金资助回国留学人员科研项目(2020-135）；山西省应用基础研究项目（202203021221212）；山西大同大学云冈专项（2020YGZX003,2020YGZX004,2020YGZX007）；山西大同大学专项（2021XSZX1）作者简介：孟田华，女，1980 年生，汉族，山西大同人，博士，副教授，主要从事文物保护技术研究。E-mail： 第 38 卷第 2 期岩土工程技术Vol.38 No.22024 年4 月Geotechnical Engineering TechniqueApr，2024响到修复效果1-2。因此，寻求一种操作简单、结果可靠、且对得胜堡遗址无损伤的修复性能检测，即墙体内部密实程度和稳定性评估的技术显得尤为重要。探地雷达（ground penetrating radar，GPR）作为一种常规的无损检测方法，具有快速、无损、操作简便、在深度方向上分辨率较高、不易受外部环境干扰等优点，已在考古发掘、矿产资源勘查、公路质量检测、管线探测等领域得到了广泛应用3。1904 年，Christian Hlsmeyer 基于无线电波回声探测原理发明了电动镜（telemobiloscope），对能见度较差的水域进行交通监控，进而可实现海上船舶相撞预防等潜在商业应用4。1975 年，Vicker 研究组首次将探地雷达应用于考古探测，分别对墨西哥州查科峡谷的三个测区进行了探测5。其他文物研究科学家陆续使用探地雷达进行考古勘察，如 Sternberg 等对美国亚利桑那州南部进行了考古研究6，Whit-ing 等对巴巴多斯岛史前遗址进行了勘察7，Poreba等对波兰 Karmanowice 地区新石器时代区域进行了调查8，Thomas 等对以色列下加利利的铜石时代遗址 Marj Rabba 进行了探地雷达调查9。国内探地雷达的研究始于 19 世纪 70 年代初，自 19 世纪 80 年代开始引进国外探地雷达技术及设备，19 世纪 90 年代以来，探地雷达得到了广泛的应用和研究10。2000 年后，探地雷达也被用于文物领域的无损检测研究，如姚萌等对古墓遗址的探地雷达探测11，祝炜平等对杭州雷峰塔遗址地下遥感考古研究12，梅宝等利用探地雷达技术在云冈石窟修复中的监控应用13，王亮等将探地雷达用于金沙遗址的考古探测14，林金鑫等对邺城遗址探地雷达考古调查研究15，赵文柯等应用探地雷达对北庭故城古建筑基址、良渚古遗址和南诏城墙古遗址分别进行了探测并得到了较好的结果16-17。方云等对北庭故城古建筑基址和奉先寺进行了勘查18。由于探地雷达探测速度快、探测过程连续、分辨率高、操作方便灵活，因此也可被应用于古城墙墙体内部病害的检测。本研究利用高频探地雷达探测了不同材料和方法修缮后城墙的内部结构，并与未修缮的古城墙进行对比，分析了三个测区的探地雷达图像差异，进而对其修缮效果进行客观评价，为古城墙的修缮工作提供资料和技术参考。1 得胜堡城墙的探地雷达检测原理探地雷达是利用天线发射和接收高频电磁波来探测介质内部物质特性和分布规律的一种地球物理方法19。探地雷达通过发射天线向被测目标体发射4Z2+X2高频脉冲电磁波，接收天线接收反射电磁波，电磁波在墙内介质中传播，遇到存在电磁差异的界面时发生反射，脉冲波行程时间 t=/v，可根据接收到的电磁波波形、振幅强度和时间的变化特征，推断出墙内介质的空间位置，结构形态以及深度20（见图 1）。城墙距离/mO双程旅行时/ns异常体天线天线天线 图 1 探地雷达原理示意图 对墙体采用高频电磁波进行测量。根据电磁波传播理论，高频电磁波在介质中的传播服从麦克斯韦方程组，即E=Bt（1）H=J+Dt（2）B=0（3）D=（4）式中：为电荷密度，C/m3；J 为电流密度，A/m2；E 为电场强度，V/m；D 为电位移，C/m2；B 为磁感应强度，T；H 为磁场强度，A/m。Dt式（1）为微分形式的法拉第电磁感应定律，式（2）为安培电流环路定律，其中由麦克斯韦引入的一项为位移电流密度 Jd，即Jd=Dt（5）探地雷达所发射的脉冲电磁波，传播中遇到存在电磁差异界面时发生反射，其脉冲波行程时间为：t=4Z2+X2/v（6）式中：t 为电磁波的双程行程时间；z 为目的体的埋深；x 为天线距；v 为电磁波在介质中的传播速度。墙体内的介质不同，对应的电阻率、介电常数、磁导率等有较大的差异。对于城墙而言，岩石、土是主要组成部分，相对介电常数在 47 的范围内。土的成分、结构等更加复杂，土的含水量与介电常数关系很大，其表达式见式（7）。孟田华等：高频探地雷达在古城墙修复效果评价中的应用197=(1)m+SwSw+(1Sw)a（7）Sw式中：下标 m、w 和 a 分别为土骨质、水和空气；为含水饱和度；为除土壤固体以外的成分含量。探地雷达中天线的基本功能是能量转换和辐射。其中，天线的效率由式（8）确定。A=PPin=PP+Pd（8）PdPPin式中：为损耗功率；为天线辐射的功率；为进入天线的功率。UinIin天线输入阻抗是指天线馈电点呈现的阻抗值，直接决定了和馈电系统之间的匹配状态，从而影响反馈到天线上的功率及馈电系统效率。输入阻抗 Zin与输入电压和电流的关系为Zin=UinIin=Rin+jXin（9）RinXinIinI0RinR式中：为输入电阻；为输入电抗。注意，输入阻抗应等于输入电流=的辐射阻抗，其实部为=。所以，通过探地雷达数据中的振幅类异常反映地下介质横向变化特征，而频率衰减变化来确定目标体与周围介质的物性差异，相位差异可确定介质分层以及异常体的接触边缘变化。基于探地雷达这些优良的特性21，本研究利用较高分辨率为 10 ps，探测天线频率为 1600 MHz 的探地雷达，采用剖面法对得胜堡古城墙进行测量，获得城墙体内部的异常情况，进而对修复效果进行表征。其他主要技术参数见表 1。表 1 探地雷达技术参数及指标技术参数指标扫描速度128 道/秒采样模式距离模式时间窗口15 ns增益方式自动探测深度75 cm 2 得胜堡城墙的探地雷达检测区域得胜堡长城的城墙体主要由黄土夯筑而成，由于夯土城墙的坚固程度较差，很容易受到风、雨等自然侵害。为了更好地评估对古城墙的修复效果，选择了具有代表性的三段城墙体进行探地雷达检测，测线大致沿古城墙的走向东西方向布设。具体测线布置情况如下：第一测区如图 2（a）所示，选定宽 1.7 m、长 5 m的长方形区域，自下而上横向布置十条测线，各测线间距均为 17 cm。该测区的修缮方法是土坯砖垒砌法，即利用当地精土制作土坯砖，对表面进行夯土填补后用，再使用当地的精土泥进行粘接。第二测区如图 2（b）所示，选定宽 1.36 m、长 4.5 m 的长方形区域，自下而上横向布置八条测线，各测线间距均为 17 cm。该测区的修缮方法是仿照传统工艺的逐层夯筑法。即在原夯土层上，采用灰土进行逐层补夯，其中补夯所用灰土是由当地的块状生石灰磨细得到的白灰与精泥土以 2:8 的比例拌和而成。第三测区如图 2（c）所示，选定宽 1.36 m，长 2.5 m 的长方形区域，自下而上横向布置八条测线，各测线间距均为 17 cm；该测区为未经修缮但保存良好的一段古城墙墙体。(a)第一测区现场图(a)第二测区现场图(a)第三测区现场图1.7 m1.36 m4.5 m5 m102182.5 m1.36 m12821 图 2 得胜堡三个代表性测区现场图进行城墙体雷达无损检测过程为：利用小滑轮将收发一体的高频天线沿城墙体所布测线往返多次缓慢移动，覆盖测区的所有测点，以获得城墙体内部198岩土工程技术2024 年第 2 期分布的高分辨率雷达图像。2.1 第一测区探地雷达图像分析解译图 3 为得胜堡古城墙第一测区其中 4 条测线的探地雷达图像，图 3（a）3（d）分别对应第 1、第 4、第 7 和第 10 条测线的探地雷达图像。总体看第一测区各测线对应浅层部分都有三层清晰的强反射雷达异常，且与周围土存在明显的性质差异，这对应土坯砖垒砌法中精土制作的三层土坯砖。每一层是表面夯土填补后的土坯砖，由于压制密实多块土坯砖逐层排列垒砌的修复工艺，在雷达图像上表现为紧密的近乎各向同性介质，振幅几乎无变化，出现明显的分层现象。各层振幅较强的异常间多组连续的反射同向轴，则对应精土泥进行粘接土坯砖的修复工艺。但随着深度的加深，土坯砖的连续性越来越差，且界面也变得越来越模糊。表现为雷达图像剖面中出现的由湿砂土和水饱和砂土介电常数差异引起的纵向连续强反射轴，研判为主要由毛细水引起的裂隙、酥碱等病害。001.53.04.500.20.40.62.55.07.510.012.515.0时间标尺/ns深度标尺/m距离/m001.53.04.500.20.40.62.55.07.510.012.515.0时间标尺/ns深度标尺/m距离/m001.53.04.56.000.20.40.62.55.07.510.012.515.0时间标尺/ns深度标尺/m距离/m001.53.04.56.000.20.40.62.55.07.510.012.515.0时间标尺/ns深度标尺/m距离/m(a)第 1 条测线雷达图像(b)第 4 条测线雷达图像(c)第 7 条测线雷达图像(d)第 10 条测线雷达图像 图 3 第一测区选取 4 条测线的雷达图像 如图 3（a）所示，为该测区离地面最高的第 1 条测线雷达图像。在深度较浅范围反射层整体波形连续性都比较好，表明了土坯砖垒砌层的完整性，即使在振幅特征不明显的地方，依然有着一致的波形连续性，这一连续性位置可表征砖层之间的间隔，即说明虽然这一测区的土坯砖完整，但波阻抗界面在前后两层砖之间存在差异，这对应于层层垒砌前对土坯砖表面进行夯土填补工艺的界面，而且可清晰看到每层砖的间隔。可由雷达电磁波信号得到在 0.02 m、0.06m、0.09 m、0.12 m、0.16 m、0.18 m、0.21 m、0.24 m、0.27 m 处产生分界，据此可划分 00.02 m、0.020.09 m、0.090.12 m、0.120.16 m、0.160.18 m、0.180.21 m、0.210.24 m 和 0.240.27 m 共 8 个土坯砖垒砌层，土坯砖厚度约为 0.03 m。但随着深度的加深，土坯砖的连续性变差，界面也变得越来越模糊。表现为雷达图像剖面中 0.090.27 m 处出现的由湿砂土和水饱和砂土介电常数差异所引起的纵向连续强反射轴，研判为主要由毛细水引起的裂隙、酥碱等病害。由图 3（b）3（d）还可见，与第 1 测线的图 3（a）雷达剖面图像类似，均可看到浅层较清晰的土坯砖层及净土泥层，但较深部分，由于毛细水等的作用使得城墙体内的土坯砖砖体松散，电导率发生变化，与精土泥间的界面变得越来越模糊，且随着高度的降低，即由图 3（a）到图 3（d）劣化情况越来越严重，特别是接近地面的第 7 和第 10 条测线的雷达图像在深度0.090.27 m 范围间已无法辨别出土坯砖层。为了更好地分析城墙浅层随着古城墙体高度变化的修复效果，选取从低到高所测 10 条测线的探地雷达图像中分别截取了深度范围为 00.4 m 的图像进行对比分析，如图 4 所示。随着高度的降低，雷达图像上对应表面夯土层下 0.090.27 m 范围内呈现强弱相间、不连续、平行的条带状反射轴逐渐增多，这是由毛细水、自然降水渗入所引起的细小裂隙、土坯砖裂缝等所致。而且，相比浅层的三层连续土坯砖层，可识别的较深部位的第四、五、六层土坯砖层间每块土坯砖间已变得不连续，出现了间隙，清晰可见变松散的各块土坯砖。2.2 第二测区探地雷达图像分析解译图 5 为得胜堡古城墙第二测区其中 4 条测线的探地雷达图像，图 5（a）5（d）分别对应第 2、第 4、第 6 和第 8 条测线的探地雷达图像。由四幅图像清晰可见分层界面，且层与层之间大部分表现均匀，偶尔出现双曲线特征。这些特点对应逐层夯筑法中，在原夯土层上逐层补夯压制密实的良好性能表现为雷达图像的均匀性，夯筑结束后，对外侧的拍土表现为孟田华等：高频探地雷达在古城墙修复效果评价中的应用199间隔基本相同的夯筑层厚度。可由雷达反射波信号读出 0.02 m、0.04 m、0.06 m、0.08 m、0.10 m 等分界值，进而推断出约 0.02 m 的夯筑层厚度。此外，夯土原材料中可能混入的小石子表现为雷达图像中的双曲线特征。00.2深度标尺/m距离/m00.200.200.200.200.200.200.200.200123456789100.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 4.5 5.0 5.50.2 图 4 第一测区不同高度的浅层雷达图像对比图 同样，为了更好地分析城墙浅层夯筑随城墙体高度变化的修复效果，选取从低到高所测 8 条测线的探地雷达图像，并从中分别截取了深度范围为00.4 m 的图像进行对比分析，如图 6 所示。随着高度的降低，雷达图像的变化趋势基本相同，小石子引起的尖锐异常并无规律可循，且无明显的雨水侵入而产生裂隙的异常现象，说明该密实的夯筑修补方法可有效加固城墙体，对于城墙的长久保存具有明显的作用。2.3 第三测区探地雷达图像分析解译图 7 为得胜堡古城墙第三测区随机选取 4 条测线的探地雷达图像，图 7（a）7（d）分别对应第 2、第4、第 6 和第 8 条测线的探地雷达图像。比较四幅图像，发现这段古城墙体内部保持得非常好，逐层夯筑的界面清晰，且内部基本无明显异常，清晰可见分层 02.55.07.510.000.20.40.612.515.0时间标尺/ns深度标尺/m01243距离/m02.55.07.510.000.20.40.612.515.0时间标尺/ns深度标尺/m01243距离/m02.55.07.510.000.20.40.612.515.0时间标尺/ns深度标尺/m01243距离/m(a)第 2 条测线雷达图像(b)第 4 条测线雷达图像(c)第 6 条测线雷达图像(d)第 8 条测线雷达图像02.55.07.510.000.20.40.612.515.0时间标尺/ns深度标尺/m01243距离/m 图 5 第二测区选取 4 条测线的雷达图像 深度标尺/m00.200.200.200.200.200.200.200.212345678距离/m00.51.01.52.02.53.03.54.04.5 图 6 第二测区不同高度的浅层雷达图像对比图200岩土工程技术2024 年第 2 期夯筑技术所对应的分层界面，整个剖面基本无明显异常。而且由分界面深度数值 0.02 m、0.04 m、0.06 m、0.08 m 等可断定当时夯筑层的厚度为 0.02 m。此外，由图 8 所示的深度范围 00.4 m 间第三测区不同高度的浅层雷达图像对比图可见，深度为 0.080.28m 间出现了夯层的不连续，但基本无明显其他异常现象。02.55.07.510.000.20.40.612.515.0时间标尺/ns深度标尺/m距离/m00.51.02.52.01.502.55.07.510.000.20.40.612.515.0时间标尺/ns深度标尺/m距离/m00.51.02.52.01.502.55.07.510.000.20.40.612.515.0时间标尺/ns深度标尺/m距离/m(a)第 2 条测线雷达图像(b)第 4 条测线雷达图像(c)第 6 条测线雷达图像(d)第 8 条测线雷达图像00.51.02.01.502.55.07.510.000.20.40.612.515.0时间标尺/ns深度标尺/m距离/m00.51.02.01.5 图 7 第三测区选取 4 条测线的雷达图像 总体比较城墙体三个测区的雷达图像可见，第一测区土坯砖垒砌法相较第二测区和第三测区的逐层夯筑工艺更容易产生病害，如雷达图像上所呈现的强弱相间、不连续、平行的条带状反射轴所对应的每块土坯砖间的间隙为裂隙、空鼓、动植物病害等提供场所；雷达图像剖面中纵向连续强反射轴所对应的与地面较近的低层墙体所受毛细水等的侵袭及酥碱等病害。而逐层夯筑压制密实的工艺确保每层夯筑层连续性很好，无明显病害，对应雷达图像均匀，基本无异常现象，该方法可有效防治裂隙、空鼓、毛细现象及动植物等病害。所以，传统逐层夯筑工艺更适宜于古城墙的修复。3 结论利用探地雷达无损检测技术分别对得胜堡长城中使用不同修复方式的墙体及未修复墙体进行了检测，获得以下结论：（1）由雷达剖面图可有效识别修复工艺，获得土坯砖厚度、夯筑层厚度等主要参数。两种修复方式对成墙体的致密性、完整性等均有明显的改善，其中传统工艺逐层夯筑城墙方式较好，在防治病害方面表现更为突出。（2）采用土坯砖垒砌法修缮的城墙虽然在一定程度上起到了保护的作用，但由于土坯砖间连接的不密实而埋下隐患，容易为毛细水、自然降水、植物和动物等提供产生裂隙、空鼓、生物等病害的条件，对城墙的长久保存效果不是最好。（3）探地雷达方法对城墙的病害检测和修复性能评估是可行和高效的。探地雷达这一快速、无损、高精度、操作简便的技术在古城墙检测方面的应用将会得到进一步发展，成为古城墙无损探测的重要手段，为古城墙的修复和养护工程提供科学依据。参考文献 马涛.文物修复的理论、原则与程序C/文物保护科学论文集.西安:西安文物保护与修复中心,2004.1 黄克忠.走向二十一世纪的中国文物科技保护J.敦煌研究，2000，（1）：5-9.2 深度标尺/m00.200.200.200.200.200.20.20.20012345678距离/m00.51.01.52.0 图 8 第三测区不同高度的浅层雷达图像对比图孟田华等：高频探地雷达在古城墙修复效果评价中的应用201李嘉，郭成超，王复明，等.探地雷达应用概述J.地球物理学进展，2007，（2）：629-637.3 GRIFFITHS H，KNOTT P，KOCH W.ChristianHlsmeyer:Invention and Demonstration of Radar,1904J.IEEE Aerospace and Electronic SystemsMagazine，2019，34（9）：56-60.4 TANKERSLEY K B，DUNNING N P，THRESS J，et al.Evaluating soil salinity and water management in ChacoCanyon,New MexicoJ.Journal of ArchaeologicalScience Reports，2016，9：94-104.5 STERNBERG B K，MCGILL J W.Archaeology studiesin southern Arizona using ground penetrating radarJ.Journal of Applied Geophysics，1995，33（1）：209-225.6 WHITING B M，MCFARLAND D P，HACKENBERGER S.Three-dimensional GPR study ofa prehistoric site in Barbados,West IndiesJ.Journal ofApplied Geophysics，2001，47（3）：217-226.7 POREBA A，WACAW M Z，JOLANTA N C，et al.Archaeological objects in loesses recognized by GPRresearch at the site Karmanowice,PolandJ.ActaGeophysica，2007，5（4）：640-651.8 THOMAS M U，YORKE M R，MORAG M K.Ground-penetrating radar investigations at Marj Rabba,aChalcolithic site in the lower Galilee of IsraelJ.Journalof Archaeological Science，2014，46：96-106.9 陈义群，肖柏勋.论探地雷达现状与发展J.工程地球物理学报，2005，2（2）：149-155.10姚萌，刘树人，杨燕.基于探地雷达的古墓遗址探测及数据后处理方法J.遥感学报，2001，5（4）：317-11320.祝炜平，黄世强，李江林.杭州雷峰塔遗址地下遥感考古研究J.地球信息科学，2002，4（2）：104-107.12梅宝，邓世坤，胡朝彬.探地雷达技术在云冈石窟维护中的应用J.工程地球物理学报，2006，3（6）：448-451.13王亮，王绪本，李正文.探地雷达在金沙遗址考古探测中的应用研究J.物探与化探，2008，32（4）：401-403.14林金鑫,田钢,石战结,等.邺城遗址探地雷达考古调查研究C/中国地球物理学会.中国地球物理.2012.15赵文轲,田钢.应用探地雷达多属性叠合分析对北庭故城古建筑基址的勘查C/中国地球物理学会.中国地球物理.2012.16赵文轲，田钢.地下古城墙遗址的地球物理探测J.科学，2013，65（2）：4,38-41.17方云，翟国林，乔梁.探地雷达探测技术在奉先寺保护工程中的应用J.物探与化探，2014，38（4）：815-819.18罗回国，张清利.不同中心频率天线探地雷达在古城墙病害体探测中的应用研究J.岩土工程技术，2018，32（5）：226-230,265.19周伟.探地雷达法在古城墙内部病害检测上的应用J.智能城市，2018，4（11）：25-26.20宗鑫，王心源，刘传胜，等.探地雷达在地下考古遗存探测中的实验与应用J.地球信息科学学报，2016，18（2）：272-281.21收稿日期：2023-06-06202岩土工程技术2024 年第 2 期