Gravity and Magnetic Anomalies and Prospecting Prediction of Iron Polymetallic Deposits in the Downstream of Nalingguole River
-
摘要:
那陵郭勒河下游地区以第四系覆盖为主，已发现多个铁多金属矿床，同时受覆盖的影响，找矿工作难度增大。笔者以1∶5万航磁及1∶20万区域重力资料研究分析了那陵郭勒河下游重磁异常的展布特征，推断了研究区的断裂构造；根据南部基岩出露区地层、岩体与重磁场的关系，推断划分了覆盖区下地层及侵入岩的分布情况。结合地质特征与矿产的分布关系，认为中、小型铁多金属矿床分布主要与中部圈定的岩体密切相关，且位于推断岩体的边部，矿（化）点主要分布于推断断裂的附近，受断裂控制明显。圈定找矿靶区2处，为后期柴达木盆地周缘覆盖区（半覆盖区）的找矿及靶区优选提供有利的证据。
Abstract:The downstream area of the Nalingguole river is predominantly covered by Quaternary sediments. Several polymetallic ore deposits, primarily composed of iron, have been discovered in this region. However, the presence of sediment cover has increased the difficulty of prospecting activities. In this study, the author analysied of the distribution characteristics gravity and magnetic anomalies, and inferred the fault structure of the study area in the downstream structures of the Nalingguole river, focusing on 1∶50,000 aeromagnetic anomaly data and 1∶200,000 regional gravity data. Based on the relationship between strata, rock mass, gravity and magnetic field in the southern bedrock outcropping area, deduced and divided the distribution of strata and intrusive rocks under the overlying area. Moreover, the mineral occurrences are predominantly found near inferred fault lines, indicating a clear influence of fault control. Two prospective mining targets were identified, providing favorable evidence for future ore exploration and target selection in the (partially) covered areas surrounding the Qaidam basin.
-
地质灾害风险防控的前提是发现隐患，即地质灾害隐患早期识别。早期识别是指通过天、空、地等各种手段，对已经发生变形或者将来可能发生变形的地质体进行判识，识别可能发生地质灾害的位置和范围，分析地质灾害发生的可能性和风险（张茂省等，2021）。随着对地质灾害形成机理和诱发过程研究的不断深入和InSAR、LiDAR、无人机航测以及计算机信息等技术的快速发展，早期识别逐渐应用于地质灾害综合防控体系（刘传正，2019）。许强等（2019，2020）通过多学科交叉融合、多部门跨界合作、多种技术手段的综合应用，提出了重大地质灾害隐患早期识别“三查”（普查、详查、核查）体系。葛大庆等（2019）从光学遥感、InSAR、激光雷达等综合遥感测量的角度提出了以“形态、形变、形势”（三形）为观测内容的地质灾害隐患识别技术思路。综合遥感技术具有大范围连续跟踪地表形变的特性，现已成为地质灾害隐患早期识别的新技术（薛强等，2018；张茂省等，2019；张勤等，2022；张林梵，2023；贾俊等，2023；林明明等，2024）。
中国是世界上地质灾害最为严重的国家之一，经过多年努力，建成了以调查评价、监测预警、工程治理和应急处置为主要内容的地质灾害综合防治体系，防灾减灾成效显著（刘传正等，2020；殷跃平等，2022；唐辉明，2022；兰恒星等，2022；孙萍萍等，2022；杜鹏等，2024；唐然等，2024；吴铭洋等，2024）。近年来，随着极端降雨的频现和日趋强烈的人类工程扰动，地质灾害防控“两个80%”问题仍然突出，即：发生灾险情的80%不属于在册隐患点，80%的在册隐患点风险较低或从未发生灾险情。据统计，2015～2021年，黄土地区发生伤亡性地质灾害事件57起，共造成114人死亡失踪、45人受伤，其中，发生于在册隐患点的9起，发生于非在册隐患点的48起，发生于非在册隐患点的伤亡性地质灾害占比达到84.2%。暴露出地质灾害早期识别环节薄弱，早期识别还不能满足和指导地质灾害综合防控实际需要。
为精准掌握陕西省榆林地区地质灾害隐患现状，切实保护人民群众生命财产安全，笔者按照“不落一个梁峁，不丢一个沟壑，不漏一村一户”的要求，聚焦地质灾害“隐患在哪里”“风险有多高”等问题，开展了地质灾害隐患风险识别和核查，基本查明了榆林地区地质灾害隐患风险底数，构建了以早期识别—遥感核查—野外核查—风险防控为主线的地质灾害隐患风险识别、核查和防控技术方法体系，创建了“精准到坡、精准到户、精准管控”的地质灾害风险防控“榆林模式”。风险核查和防控模式在榆林地区汛期地质灾害防治中得到成功应用，增强了地质灾害防治的针对性和有效性，显著减轻了地质灾害造成的损失。
1. 风险识别内容与识别方法
1.1 风险识别和核查内容
地质灾害隐患风险识别和核查主要以县域为单元，坚持“以人为本、预防为主”的防灾原则，充分应用合成孔径雷达干涉测量技术（InSAR）、高精度数字高程模型（DEM）、航空遥感及无人机航测、地质灾害数字采集系统（APP）等新技术新方法，采用人力资源与现代科技相结合，航天航空与地面调查相结合，技术专家与干部群众相结合，以技术专家为主的风险核查模式，形成以群众报险报灾为基础，以干部为推动力，以技术专家为核心的地质灾害隐患识别和核查格局。对县域范围进行全面、系统的地质灾害隐患早期识别和风险核查，不留空白和死角，不放过一个隐患点，切实解决地质灾害调查精度低、早期识别环节薄弱、隐患点风险底数不清等一系列问题，通过掌握地质灾害隐患最新动态变化情况，增强地质灾害防控的针对性和有效性，实现防灾减灾3个转变，努力把地质灾害风险和损失降至最低，并建立地质灾害隐患识别和核查示范。具体识别和核查内容包括以下几个方面。
（1）高陡危险坡段识别和核查（DEM和遥感识别隐患点核查）：对通过DEM（2 m×2 m精度）和光学遥感（0.2 m分辨率）识别出来的危险坡段（坡度大于40°、坡高大于20 m且有威胁对象的斜坡段），逐一进行现场核查，研判其变形破坏模式，圈定可能的滑坡、崩塌范围，调查其威胁对象。主要核查隐患点类型、地质环境条件、主要特征、威胁对象，并进行危险性和风险评价，提出防治措施建议。本次榆林地区地质灾害隐患风险核查主要是基于DEM和遥感的地质灾害隐患点核查，核查过程包括DEM易发坡段识别→遥感危险坡段识别→风险野外核查。
（2）地表变形区段识别和核查（InSAR识别隐患点核查）：对通过InSAR技术识别出来的地表变形区段进行高精度遥感核查，识别是否存在威胁对象，对有威胁对象的疑似隐患点逐一开展现场核查，研判其变形破坏模式，圈定可能的滑坡、崩塌范围，调查其威胁对象。主要核查隐患点类型、地质环境条件、变形特征、威胁对象，并进行危险性和风险评价，提出防治措施建议。核查过程包括InSAR地面形变识别→遥感威胁对象核查→风险野外核查。
（3）群众报险报灾点核查（现场新发现隐患点核查）：群众报险报灾主要针对通过DEM、InSAR等手段未识别出，但出现变形迹象的隐患点。现场核查工作启动之后，充分发动群众，以村为单元，开展地质灾害群众报险报灾，上报坡体或房屋裂缝、渗水、鼓胀、滑塌、沉降、塌陷等变形信息，并对所在村域进行全面排查。主要核查隐患点类型、地质环境条件、变形破坏迹象、威胁对象，并进行危险性和风险评价，提出防治措施建议。
（4）在册地质灾害隐患点核查：现场核查已发生滑坡、崩塌、泥石流等隐患点及治理工程的变化情况，对其复活的可能性和危险性进行评估，重点核查滑坡崩塌裂缝、崩塌危岩体、泥石流沟物源有无新的发展，受威胁对象和危险区范围有无变化，监测仪器、警示标志和通讯设备是否完好，日常监测预警工作是否保持有效开展等。主要核查在册点地质环境条件、主要特征、威胁对象，以及防治工程、监测预警等风险减缓措施效果，并进行危险性和风险评价，提出防治措施建议。
1.2 风险识别和核查方法
地质灾害风险识别和核查总体按资料收集→早期识别→遥感核查→野外核查→风险管控的流程开展。充分应用和发挥新技术新方法的优势，研发地质灾害风险核查APP，建立基于DEM、InSAR的疑似隐患点识别阈值标准和基于高精度遥感影像的疑似隐患点识别方法。
具体技术路线为：在资料收集的基础上，基于高精度DEM，通过坡度和坡高识别高陡易发坡段，基于InSAR技术，通过形变监测识别地表变形区段；开展基于建筑物、基础设施等承灾体和基于高陡易发坡段、地表变形区段等致灾体的高精度遥感核查，圈定疑似滑坡、崩塌等地质灾害隐患；开展高陡危险坡段（DEM和遥感识别隐患点）、地表变形区段（InSAR识别隐患点）、群众报险报灾点（现场新发现隐患点）、在册地质灾害隐患点野外核查，现场填写地质灾害风险核查APP，建立地质灾害隐患点数据库；开展地质灾害风险管控对策研究，形成地质灾害风险识别和核查示范（图1）。
![]()
图 1 地质灾害风险识别和核查技术路线图Figure 1. Technical roadmap for geo-hazards risk identification and verification2. 研究区概况
榆林市位于陕西省最北部，地域东西长385 km，南北宽263 km，总面积为
43587 km2。地貌大体以长城为界，北部为风沙草滩区，地形较平坦，地貌类型以滩地和沙地为主；南部为黄土丘陵沟壑区，地形破碎，切割强烈，地貌以黄土梁峁和沟壑为主（图2）。榆林市多年平均降水量为402.29 mm，年内降水量分配不均，降水主要集中在6～9月，7、8月最大，多在60～100 mm。榆林市辖1市2区9县、155个乡镇、29个街道办事处、2967 个行政村，户籍人口385.59万人。榆林市地质环境条件脆弱，滑坡、崩塌等地质灾害多发频发，是陕西省地质灾害重灾区，滑坡、崩塌等地质灾害在空间上呈现“点多面广，东南集中”的分布规律。根据最新1∶5万地质灾害风险调查成果，榆林市地质灾害易发区总面积
32646 km2，占市域面积的74.9%，其中高易发区面积约7622 km2，中易发区面积约13267 km2，低易发区面积约11757 km2。榆林市现有在册地质灾害隐患点1589 处，其中滑坡541处、崩塌1005 处、泥石流18处、地面塌陷25处，滑坡、崩塌、泥石流地质灾害主要分布在榆林市东南部清涧县、绥德县、吴堡县、子洲县、米脂县、佳县等县黄土丘陵沟壑区，地面塌陷地质灾害主要分布在榆林市北部神木市、府谷县等市（县）采煤塌陷区。近年来，陕北黄土高原地区极端降雨频现，人工切坡活动加剧，近期发生的伤亡性地质灾害大多与降雨和人工切坡有关，降雨和人工切坡成为榆林市滑坡、崩塌等地质灾害的主要诱发因素。受地形限制，榆林市东南部清涧县、绥德县、吴堡县、子洲县、米脂县、佳县等县域大部分工程活动都涉及人工切坡，农村切坡建房（约占农户建房总数的70%）和沿山修建乡村道路也相当普遍，地质灾害防范压力巨大。
3. 地质灾害隐患早期识别
3.1 基于DEM的地质灾害隐患识别
榆林市地质灾害类型主要以黄土滑坡、崩塌为主，黄土滑坡、崩塌的发生主要受斜坡坡度、坡高等地形条件的控制和影响。根据不同坡度、坡高发生黄土滑坡、崩塌的比率（薛强等，2023），建立榆林市基于DEM的疑似隐患点识别阈值标准，开展地质灾害隐患易发坡段识别，为地质灾害隐患遥感核查提供有效靶区。榆林地区黄土滑坡、崩塌的发生主要集中在坡度大于40°、坡高大于20 m的高陡斜坡地带（张茂省等，2011），因此将坡度大于40°、坡高大于20 m作为榆林市黄土滑坡、崩塌隐患易发坡段识别的阈值标准。
基于ArcGIS平台利用2 m×2 m高精度DEM数据提取滑坡、崩塌隐患易发坡段，共提取坡度大于40°、坡高大于20 m的高陡易发坡段
493380 处，其中榆阳区38654 处、神木市36000 处、府谷县64748 处、横山区29528 处、定边县43396 处、靖边县36943 处、绥德县52210 处、吴堡县36724 处、佳县34382 处、子洲县38077 处、米脂县44716 处、清涧县38002 处（图3）。![]()
图 3 榆林市米脂县局部基于DEM的高陡易发坡段识别结果图Figure 3. Identification map of high and steep landslide-prone slopes based on DEM in local area of Mizhi County, Yulin City3.2 基于InSAR的地质灾害隐患识别
地质灾害孕育演化过程中最直观的表象就是地表的形变，合成孔径雷达干涉测量技术（InSAR）利用合成孔径雷达两次观测中雷达波相位差与空间距离之间的关系可有效获取区域地表三维形变信息，能够全天候、全天时获取数据，且不受天气影响。特别是InSAR所具有的大范围连续跟踪微小形变的特性，使其成为地质灾害隐患早期识别和监测的全新技术。InSAR技术的选择应根据观测对象、应用环境、观测精度、可观测的量程、所需数据量和观测频率、技术复杂程度等因素综合确定（郭华东，2000；廖明生等，2014；张茂省等，2021）。
本次采用2020年1月～2022年9月Sentinel-1升轨数据，首先采用Stacking InSAR方法初步确定形变危险区域，然后选择相干性高的干涉采用小基线集（SBAS）方法计算，获取区内最新地表变形速率。将形变结果与光学遥感影像、DEM等反映地物及地形特征的图层要素叠加，进行专业目视解译分析，可根据不同地物在光学影像中不同的纹理、色彩、亮度及地形地貌特征等判断滑坡、崩塌发育迹象，例如较明显、易解译的拉张裂缝、堆积平台等特征。对形变结果与遥感影像特征、地形特征较为吻合的斜坡进行圈定，获取基于InSAR技术识别的高形变区域（韩守富等，2020；苏晓军等，2021）。
基于榆林市历史灾害特征，结合高形变区域监测结果，圈定变形速率较大的不稳定斜坡，识别疑似地质灾害隐患（图4）。榆林市基于InSAR技术共识别出疑似地质灾害隐患区段635处。
![]()
图 4 榆林市米脂县典型区域基于InSAR技术识别的疑似地质灾害隐患Figure 4. Potential geo-hazards identified by InSAR technology in representative region in Mizhi County, Yulin City4. 地质灾害隐患精准核查
4.1 精准到坡式遥感核查
在DEM高陡易发坡段识别、InSAR地表变形区段识别的基础上，通过0.2 m分辨率遥感数据开展危险坡段识别和核查，为地质灾害风险实地核查提供有效靶区。遥感核查分为承灾体核查和致灾体核查两方面内容：承灾体核查主要核查研究区建筑物、基础设施等受威胁对象，并充分利用已有农房数据，对受到滑坡、崩塌隐患威胁的承灾体识别危险坡段；致灾体核查主要是核查高陡易发坡段和地表变形区段是否有威胁对象，对有威胁对象的高陡易发坡段和地表变形区段识别为危险坡段，对没有威胁对象的高陡易发坡段和地表变形区段不作为危险坡段。
榆林市基于高精度遥感和农房数据共识别和核查危险坡段（地质灾害疑似隐患）
31988 处，其中榆阳区1519 处、神木市2073 处、府谷县1416 处、横山区1176 处、定边县1186 处、靖边县1890 处、绥德县3529 处、吴堡县1519 处、佳县2580 处、子洲县6164 处、米脂县4198 处、清涧县4738 处（图5）。遥感识别和核查危险坡段数量占DEM、InSAR识别高陡易发坡段和地表变形区段数量的6.48%，缩小了地质灾害隐患识别范围，提高了识别精准性，为野外核查提供了有效靶区，实现地质灾害隐患识别“精准到坡”。![]()
图 5 榆林市米脂县局部基于遥感的危险坡段识别图Figure 5. Identification map of dangerous slope section based on remote sensing in local area of Mizhi County, Yulin City4.2 精准到户式野外核查
（1）核查方式：地质灾害隐患野外核查主要内容包括高陡危险坡段核查（DEM和遥感识别隐患点核查）、地表变形区段核查（InSAR识别隐患点核查）、群众报险报灾点核查（现场新发现隐患点核查）、在册地质灾害隐患点核查。野外核查工作采用技术专家、干部、群众相结合的工作模式，①充分发动群众报险报灾，上报坡体或房屋裂缝、渗水、鼓胀、沉降、塌陷等变形信息。②技术专家对DEM和遥感识别高陡危险坡段、InSAR识别地表变形区段、群众报险报灾点、在册地质灾害隐患点等逐一开展野外实地核查，根据坡体结构、变形迹象、威胁对象等进行风险分级，现场填写核查APP系统并实时上传提交至服务器平台，形成隐患点数据库（图6）。③野外实地核查中由乡村干部带队逐坡逐沟核查，每个隐患点以“村组+户主姓名”命名，精准摸清地质灾害隐患风险底数，形成以乡村干部和群众为基础、县级干部为助力、技术专家为核心的地质灾害隐患风险核查格局，实现地质灾害隐患核查“精准到户”。
（2）核查结果：野外根据坡体结构、变形迹象等判定斜坡失稳概率（即危险性），在此基础上根据威胁对象特征（即危害性）综合判定斜坡风险（薛强等，2023），将地质灾害隐患风险等级划分为极高风险、高风险、中风险、低风险和极低风险5级（图7）。榆林市共实地核查地质灾害隐患点
37523 处，其中极高风险隐患点82处，占总数的0.22%；高风险隐患点1378 处，占总数的3.67%；中风险隐患点2895 处，占总数的7.72%；低风险隐患点4946 处，占总数的13.18%；极低风险隐患点28222 处，占总数的75.21%（表1）。![]()
图 7 榆林市米脂县局部地质灾害隐患点分布图Figure 7. Distribution map of potential geo-hazards in local area of Mizhi County, Yulin City表 1 榆林市各县（市、区）核查地质灾害隐患点统计表Table 1. Statistical table of potential geo-hazards of verification in each county, Yulin City序号 县（市、区） 地质灾害隐患点（处） 合计 极高风险 高风险 中风险 低风险 极低风险 1 榆阳区 3 31 200 414 1037 1685 2 神木市 3 49 179 495 1326 2052 3 府谷县 6 100 109 236 956 1407 4 横山区 13 50 95 199 2299 2656 5 定边县 9 70 151 275 1650 2155 6 靖边县 1 176 142 288 1347 1954 7 绥德县 8 119 112 346 2546 3131 8 吴堡县 0 25 145 231 1356 1757 9 佳县 2 205 297 516 3203 4223 10 子州县 7 107 567 728 4199 5608 11 米脂县 7 183 391 596 3233 4410 12 清涧县 23 263 507 622 5070 6485 合计 82 1378 2895 4946 28222 37523 榆林市原有在册地质灾害隐患点
1589 处。通过对原有在册隐患点的实地核查，确定为极高风险的在册隐患点32处，占总数的2.02%；高风险在册隐患点137处，占总数的8.62%；中风险在册隐患点423处，占总数的26.62%；低风险在册隐患点396处，占总数的24.92%；极低风险在册隐患点601处，占总数的37.82%（图8）。可以看出，在册隐患点以中、低和极低风险为主，主要原因为部分在册隐患点采取了搬迁避让、工程治理、监测预警等防治措施，发生地质灾害的可能性较小，风险等级降低。同时，极高、高风险点中分别有50处、1241 处为非在册隐患点，分别占极高、高风险点总数的60.98%、90.06%（图9）。进一步说明了地质灾害防控面临的“两个80%”问题，即：发生灾险情的80%不属于在册隐患点，80%的在册隐患点风险较低或从未发生灾险情。![]()
图 8 榆林市在册地质灾害隐患点风险等级柱状图Figure 8. Histogram of risk grade of registered geo-hazards in Yulin City![]()
图 9 榆林市极高、高风险点中在册、非在册隐患点占比Figure 9. Proportion of registered and unregistered geo-hazards in extremely high and high risk points in Yulin City5. 地质灾害隐患风险精准管控
地质灾害隐患风险防控一般从两方面入手，一是通过坡面排水、坡脚加固、排危除险、工程治理等处置措施降低发生滑坡、崩塌的可能性；二是通过避险搬迁、监测预警、应急撤离等防控措施减少滑坡、崩塌等灾害发生之后造成的损失。笔者结合榆林市地质灾害防治的实际情况，从地质灾害防治适宜性和成本效益性角度出发，按照科学发展、人与环境协调的要求，提出地质灾害隐患防控措施建议。地质灾害隐患总体防控原则为：以人为本，生命至上；及时整改，确保安全；简单易行，注重实效；应急处置，统筹长远。
5.1 风险预警等级划分与预警响应
依据中国气象预警暴雨划分等级，并结合榆林市连续降雨情况，将地质灾害风险预警等级由低到高分别划分为蓝色预警、黄色预警、橙色预警、红色预警4级（表2），分别代表发生地质灾害的风险一般、风险较高、风险高、风险很高（董永路，2023）。
表 2 榆林市地质灾害预报预警等级划分标准Table 2. Classification standard of geo-hazards forecast and early warning in Yulin City预警等级 未来降雨量（气象预报） 前期降雨量 发生地质灾害风险 蓝色预警 气象预报12小时内降雨量将达50 mm以上 前24小时累计降雨量已达50 mm以上或
前3天累计降雨量已达100 mm以上风险一般 黄色预警 气象预报6小时内降雨量将达50 mm以上 前12小时累计降雨量已达50 mm以上或
前3天累计降雨量已达150 mm以上风险较高 橙色预警 气象预报3小时内降雨量将达50 mm以上 前6小时累计降雨量已达50 mm以上或
前3天累计降雨量已达200 mm以上风险高 红色预警 气象预报3小时内降雨量将达100 mm以上 前6小时累计降雨量已达100 mm以上或
前3天累计降雨量已达250 mm以上风险很高 注：未来降雨量和前期降雨量满足一项即可启动相应的预警等级。 蓝色预警：气象预报12小时内降雨量将达50 mm以上，或者前24小时累计降雨量已达50 mm以上或前期3天累计降雨量已达100 mm以上且气象预报降雨将持续。
黄色预警：气象预报6小时内降雨量将达50 mm以上，或者前12小时累计降雨量已达50 mm以上或前期3天累计降雨量已达150 mm以上且气象预报降雨将持续。
橙色预警：气象预报3小时内降雨量将达50 mm以上，或者前6小时累计降雨量已达50 mm以上或前期3天累计降雨量已达200 mm以上且气象预报降雨将持续。
红色预警：气象预报3小时内降雨量将达100 mm以上，或者前6小时累计降雨量已达100 mm以上或前期3天累计降雨量已达250 mm以上且气象预报降雨将持续。
根据风险核查结果和风险预警等级，制定不同风险级别地质灾害隐患预警响应方案（表3）。
表 3 不同风险隐患点预警防御方案应对简表Table 3. Early warning and defense scheme response table for potential geo-hazards of different risk levels风险等级 蓝色预警 黄色预警 橙色预警 红色预警 极高风险隐患点 应急搬迁避让 应急搬迁避让 应急搬迁避让 应急搬迁避让 高风险隐患点 加强监测 加强监测 应急搬迁避让 应急搬迁避让 中风险隐患点 加强监测 加强监测 应急搬迁避让 低风险隐患点 加强监测 加强监测 极低风险隐患点 加强监测 5.2 风险防控措施
按照不同风险级别地质灾害预警响应方案，需要重点管控的极高和高风险点共
1460 处，其中在册隐患点169处，仅占11.58%，进一步增强了地质灾害防控的有效性和精准性，减轻了地方政府汛期按风险区大范围临时撤离防控压力。联合地方政府部门，制定了地质灾害风险点管控方案，对82处极高风险点列入市安委会重大风险隐患督办台账。其中，对18处有临灾迹象的极高风险点全部转移撤离，对其余64处极高风险点要求按期实施综合治理；对
1378 处高风险点交由县级自然资源部门重点管控，并分期实施避险搬迁和综合治理；对2895 处中风险点和4946 处低风险点交由乡镇、街道办做好安全防范工作；对28222 处极低风险点交由村组进行日常监测管理。形成市、县、乡镇、村组分级防控格局，实现地质灾害隐患风险“精准管控”。6. 风险核查和防控模式应用
2022年，汛期榆林市出现6次极端强降雨过程，平均降水量592 mm，部分地区降雨量超过700 mm。榆林市根据风险核查结果和不同风险隐患点预警防御方案，及时发布风险预警，紧急避险转移
1058 处隐患点受威胁群众3703 户9861 人。2022～2023年，汛期榆林市上报地质灾害灾情或险情60起，新上报灾险情全部位于风险核查隐患点之列，其中极高风险隐患点3处、高风险隐患点40处、中风险隐患点12处、低风险隐患点5处，极高和高风险隐患点占比为71.67%（图10、图11）。由于临时撤离、应急处置等风险管控措施得当，各隐患点发生灾险情后均未造成人员伤亡，榆林市实现了丰水年汛期地质灾害“零伤亡”（薛强等，2023）。![]()
图 10 米脂县龙镇苏家沟村崩塌隐患风险核查结果与实际发生灾害对比（据薛强等，2023）a.风险核查结果遥感影像；b.灾害发生情况野外照片Figure 10. Comparison of risk verification result with actual disaster of Sujiagou collapse in Longzhen Town, Mizhi County![]()
图 11 米脂县杨家沟镇管家咀村滑坡隐患风险核查结果与实际发生灾害对比a.风险核查结果遥感影像；b.灾害发生情况野外照片Figure 11. Comparison of risk verification result with actual disaster of Guanjiazui landslide in Yangjiagou Town, Mizhi County2022～2023年，按照地质灾害风险点管控方案，榆林市对极高和高风险隐患点分期实施了避险搬迁和综合治理措施，有效控制了地质灾害隐患风险。在此基础上，对经过避险搬迁、综合治理的地质灾害隐患点核销或降低风险等级，对风险核查之后出现险情的地质灾害隐患点提高风险等级，实现了地质灾害隐患风险动态管控（图12～图16）。
![]()
图 12 米脂县城郊镇官庄村崩塌隐患风险动态管控过程a.野外核查为高风险隐患点；b.汛期发生2处滑塌；c.2023年封口搬迁Figure 12. Dynamic management and control process of Guanzhuang collapse risk in Chengjiao Town, Mizhi County![]()
图 13 米脂县沙家店镇牛圈塔村皇帝沟路旁崩塌隐患风险动态管控过程a.野外核查为高风险隐患点；b.汛期发生滑塌；c.2023年已完成工程治理Figure 13. Dynamic management and control process of Huangdigou road collapse risk in Niuquanta Village, Shajiadian Town, Mizhi County![]()
图 14 米脂县龙镇丁家沟村滑坡隐患风险动态管控过程a.野外核查为极高风险隐患点；b.2023年正在实施工程治理Figure 14. Dynamic management and control process of Dingjiagou landslide risk in Longzhen Town, Mizhi County![]()
图 15 米脂县印斗镇常家坪村崩塌隐患风险动态管控过程a.野外核查为高风险隐患点；b.2023年已完成工程治理Figure 15. Dynamic management and control process of Changjiaping collapse risk in Yindou Town, Mizhi County![]()
图 16 米脂县杨家沟镇何家岔村崩塌隐患风险动态管控过程a.野外核查为高风险隐患点；b.2023年已完成工程治理Figure 16. Dynamic management and control process of Hejiacha collapse risk in Yangjiagou Town, Mizhi County7. 结论
（1）根据榆林地区孕灾条件和黄土崩滑灾害发育分布特征，构建了以InSAR地表形变监测–高精度DEM高陡易发坡段识别–高精度遥感危险坡段核查–野外实地核查验证为主线的黄土地质灾害精细识别技术方法体系，创建了“精准到坡、精准到户、精准管控”的榆林地质灾害精细识别和风险防控模式。
（2）通过InSAR、DEM、光学遥感、野外实地核查等方法共识别榆林市地质灾害隐患点
37523 处，其中极高风险隐患点82处、高风险隐患点1378 处、中风险隐患点2895 处、低风险隐患点4946 处、极低风险隐患点28222 处，解决了榆林市地质灾害“隐患在哪里”“风险有多高”等关键问题。（3）按照风险等级，建立了榆林市地质灾害风险点预警防御和管控方案，风险识别、核查和防控模式得到及时、成功应用，有效控制了地质灾害隐患风险。
致谢：衷心感谢各县（市、区）核查单位和榆林市各级政府部门，在组织管理、系统研发、数据提供、野外核查等方面提供的帮助！
-
![]()
图 1 那陵郭勒河下游地质矿产图
1.第四系；2.下石炭统大干沟组；3.下石炭统石拐子组；4.下泥盆统契盖苏组；5.祁漫塔格群碳酸盐岩组；6.祁漫塔格群火山岩组；7.祁漫塔格群碎屑岩组；8.古元古代金水口岩群片麻岩组；9.晚三叠世石英闪长岩；10. 晚泥盆世石英二长岩；11.中泥盆世二长闪长岩；12.晚泥盆世石英闪长岩；13.中三叠世石英闪长岩；14.中泥盆世二长花岗岩；15. 早泥盆世二长花岗岩；16. 中泥盆世辉长岩；17.晚泥盆世花岗闪长岩；18.晚志留世花岗闪长岩；19.中三叠世花岗闪长岩；20.晚三叠世二长花岗岩；21.晚三叠世花岗闪长岩；22.中三叠世二长花岗闪；23.晚三叠世正长花岗岩；24.新元古代二长花岗片麻岩；25.铁多金属矿化点；26.铁多金属矿点；27.铁多金属小型矿床；28.铁多金属中型矿床及编号；29.地质界线；30.断裂；31.构造单元界线及编号；32.区域性大断裂；33山峰及名称；34.河流
Figure 1. Map of geology and mineral resources in the downstream of Nalingguole river
![]()
图 2 那陵郭勒河下游布格重力异常与断裂矿点套合图
1.铁多金属矿化点；2.铁多金属矿点；3.铁多金属小型矿床；4.铁多金属中型矿床及编号；5.推断主断裂及编号；6.推断次级断裂及编号
Figure 2. The combination map of bouguer gravity anomaly and fracture mine in the downstream of Nalingguole River
![]()
图 3 那陵郭勒河下游1∶5万航磁（△T）化极等值线及断裂矿点叠合图
1.铁多金属矿化点；2.铁多金属矿点；3.铁多金属小型矿床；4.铁多金属中型矿床及编号；5.推断主断裂及编号；6.推断次级断裂及编号；7.航磁（△T）化极等值线
Figure 3. The combination map of 1∶50,000 aeromagnetic (△T) polarizing contour and fracture mine in the downstream of Nalingguole river
![]()
图 4 那陵郭勒河下游1:5万航磁（△T）等值线与断裂矿点套合图
1.航磁（△T）等值线；2.铁多金属矿化点；3.铁多金属矿点；4.铁多金属小型矿床；5.铁多金属中型矿床及编号；6.推断主断裂及编号；7.推断次级断裂及编号
Figure 4. The combination map of 1:50,000 aeromagnetic (△T) contour and fracture mine in the downstream of Nalingguole river
![]()
图 5 布格重力水平总梯度模与断裂构造叠加图
1. 推断主断裂及编号；2. 推断次级断裂及编号
Figure 5. Map of bouguer gravity horizontal total gradient model and superposition diagram of fault structure
![]()
图 6 那陵郭勒河下游剩余重力异常和航磁（△T）化极叠合推断地质成果图
1.下石炭统大干沟组;2.下石炭统石拐子组；3.下泥盆统契盖苏组；4.祁漫塔格群碳酸盐岩组；5.祁漫塔格群火山岩组；6.祁漫塔格群碎屑岩组；7.中三叠世花岗闪长岩； 8.晚三叠世石英闪长岩； 9.中泥盆世二长闪长岩； 10.晚三叠世正长花岗岩；11.未分花岗闪长岩+二长花岗；12.铁多金属矿化点；13.铁多金属矿点；14.铁多金属小型矿床；15.铁多金属中型矿床及编号；16.推断地质界线；17.第四系界线；18.推断的找矿靶区及编号；19.航磁（△T）化极正等值线；20.航磁（△T）化极负等值线；21.航磁（△T）化极零等值线
Figure 6. The inferred geological results map of residual gravity anomaly and aeromagnetic (△T) polarization in the downstream of Nalingguole River
表 1 那陵郭勒河下游岩石物性统计表
Table 1 Statistical table of rock physical properties in the downstream in Nalingguole river
大类 代号 岩石名称 数量（块） 磁化率（10−6×4π·SI） 密度（g/cm3） 变化范围 平均值 变化范围 平均值 地层 Q 黏土、砂、砾石 30 1.43～2.24 1.83 C1d 灰岩 19 1.3～5.0 1.9 2.63～2.84 2.76 C1s 复成分岩屑砂岩 5 7.5～8.8 8.8 2.59～2.66 2.61 OQ3 灰岩 94 1.3～719.5 80.1 2.57～2.82 2.73 大理岩 74 18.2～1462.2 379.9 2.56～2.83 2.73 OQ2 大理岩 60 14.2～584.3 90.1 2.43～3.12 2.89 构造角砾岩 17 29.3～255.5 92.4 2.38～2.82 2.58 蚀变安山质玄武岩 50 16.8～888.9 88.2 2.84～3.14 3.02 OQ1 长石石英砂岩 60 37.9～1507.2 151.9 2.52～2.82 2.70 粉砂质板岩 120 2.5～3353.5 28.9 1.72～2.84 2.72 石英质板岩 5 2.5～23.9 15.1 2.69～2.82 2.71 Pt1J 片麻岩 82 3.8～45.2 22.6 2.59～3.07 2.73 侵入岩 正长花岗岩 12 2.5～2800.9 12.3 2.52～2.66 2.61 二长花岗岩 238 5.0～4797.9 348.8 2.54～2.87 2.62 花岗闪长岩 36 110.0～3007.8 783.8 弱磁铁矿化花岗岩 35 526.0～4489.1 2355.7 2.64～2.67 2.65 磁铁矿化花岗岩 40 290.4～6013.8 2556.8 2.65～2.68 2.66 矿石 磁铁矿 54 10209.2～467590.4 76112.0 3.24～4.72 4.17 
下载: 导出CSV
表 2 那陵郭勒河下游地层及侵入岩物性及重磁场特征统计表
Table 2 Statistical table of stratigraphic and intrusive petrophysical and gravity magnetic characteristics in the downstream in Nalingguole river
大类 代号 物性特征 可形成的重磁异常场 出露区重磁异常特征 地层 Q 低密度 重力低异常、负磁异常或无磁异常 东北角重力低异常、宽缓的负磁磁异常 C1d 高密度，低磁化率 重力高异常、负磁异常或无磁异常 重力高异常、正磁异常背景，局部有磁力高和磁力低伴生的磁异常 C1s 低密度，低磁化率 重力低异常、负磁异常或无磁异常 重力低、负磁异常带 OQ 高密度、中低磁化率 重力高异常、整体呈负磁异常或无磁异常，局部有弱磁异常 重力高异常、负磁异常为主，均布由正负伴生的磁异常 Pt1J 高密度、低磁化率 重力高异常、负磁异常或无磁异常 露头处形成重力高及低缓磁异常 侵入岩 低密度、中高磁化率 重力低异常、正磁异常 重力低异常、正磁异常区（西南角） 矿体 高密度、高磁化率 重力高异常、正磁异常 重力场梯级带，偏正异常，磁法正异常
下载: 导出CSV
-
白生龙,刘国燕,李建亮,等. 东昆仑沙丘低缓磁异常区深部铁多金属矿床找矿标志[J]. 矿产勘查,2019,10(1):118–124. doi: 10.3969/j.issn.1674-7801.2019.01.015 BAI Shenglong, LIU Guoyan, LI Jianliang, et al. Study on the prospecting signs of deep iron polymetallic deposits in the low-rate magnetic anomaly area of the east Kunlun Dune[J]. Mineral Exploration, 2019, 10(1): 118-124. doi: 10.3969/j.issn.1674-7801.2019.01.015
曹德智,袁桂林,郑振华,等. 它温查汉西铁多金属矿地球物理特征及磁法物探找矿方法优选研究[J]. 青海大学学报(自然科学版),2014,32(4):71–79. doi: 10.13901/j.cnki.qhwxxbzk.2014.04.014 CAO Dezhi, YUAN Guilin, ZHENG Zhenhua, et al. Study on optimization of magnetic--geophysical prospecting methods in thickly covered area in Tawenchahan in East Kunlun in Qinghai[J]. Journal of Qinghai University(Natural Science Edition), 2014, 32(4): 71-79. doi: 10.13901/j.cnki.qhwxxbzk.2014.04.014
丰成友,李东生,吴正寿,等. 东昆仑祁漫塔格成矿带矿床类型、时空分布及多金属成矿作用[J]. 西北地质,2010,43(4):10–17. doi: 10.3969/j.issn.1009-6248.2010.04.002 FENG Chengyou, LI Dongsheng, WU Zhengshou, et al. Major types, time-Space distribution and metallogeneses of polymetallic deposits in the qimantage metallogenic belt, Eastern Kunlun Area [J]. Northwestern Geology, 2010, 43(4): 10-17. doi: 10.3969/j.issn.1009-6248.2010.04.002
高鹏,耿涛,冀显坤,等. 东昆仑祁漫塔格地区激电测量中常见问题及解决方法[J]. 西北地质,2017,50(4):232–237. doi: 10.3969/j.issn.1009-6248.2017.04.025 GAO Peng, GENG Tao, JI Xiankun, et al. Common Problems and Countermeasures about the Induced Polarization Method Used in Qimantage area, East Kunlun [J]. Northwestern Geology, 2017, 50(4): 232-237. doi: 10.3969/j.issn.1009-6248.2017.04.025
郭广慧, 钟世华, 李三忠, 等. 运用机器学习和锆石微量元素构建花岗岩成矿潜力判别图解: 以东昆仑祁漫塔格为例[J/OL]. 西北地质, 2023: 1−14. doi: 10.12401/j.nwg.2023158 GUO Guanghui, ZHONG Shihua, LI Sanzhong, et al. Constructing Discrimination Diagrams for Granite Mineralization Potential by Using Machine Learning and Zircon Trace Elements: Example from the Qimantagh, East Kunlun[J/OL]. Northwestern Geology, 2023: 1−14. doi: 10.12401/j.nwg.2023158
李东生,张文权,田承盛,等. 青海祁漫塔格地区主要矿床类型找矿方法探讨[J]. 西北地质,2013,46(4):131–141. doi: 10.3969/j.issn.1009-6248.2013.04.012 LI Dongsheng, ZHANG Wenquan, TIAN Chengsheng, et al. Discussion on the metallogenic characteristics and ore-prospecting methods of qimantage Region, Qinghai Province [J]. Northwestern Geology, 2013, 46(4): 131-141. doi: 10.3969/j.issn.1009-6248.2013.04.012
李文渊. 祁漫塔格找矿远景区地质组成及勘查潜力[J]. 西北地质,2010,43(4):1–9. LI Wenyuan. The Geological Composition and Metallogenetic Prospect in the Qimantage Prospective Region, East Kunlun [J]. Northwestern Geology, 2010, 43(4): 1-9.
李玉春,张爱奎,张培青,等. 青海祁漫塔格沙丘地区侵入岩地质、地球化学特征及找矿意义[J]. 西北地质,2013,46(3):70–82. doi: 10.3969/j.issn.1009-6248.2013.03.005 LI Yuchun, ZHANG Aikui, ZHANG Peiqing, et al. The Geological and Geochemical Characteristics of the Intrusive Rock and Prospecting Significance in Shayiu Areas, Qimantage Metallogenic Belt, Qinghai Province [J], Northwestern Geology, 2013, 46(3): 70-82. doi: 10.3969/j.issn.1009-6248.2013.03.005
刘彦,严加永,吴明安. 基于重力异常分离方法寻找深部隐伏铁矿-以安徽泥河铁矿为例[J]. 地球物理学报,2012,55(12):4181–4193. LIU Yan, YAN Jiayong, WU Mingan. Exploring deep concealed ore bodies based on gravity separation methods: A case study of the Nihe iron deposit [J]. Chinese Journal of Geophysics, 2012, 55(12): 4181- 4193.
刘嘉情, 钟世华, 李三忠, 等. 基于机器学习和全岩成分识别东昆仑祁漫塔格斑岩-矽卡岩矿床成矿岩体和贫矿岩体[J/OL]. 西北地质, 2023: 1−16. doi: 10.12401/j.nwg.2023155 LIU Jiaqing, ZHONG Shihua, LI Sanzhong, et al. Identification of Mineralized and Barren Magmatic Rocks from the Qimantagh, East Kunlun Based on Machine Learning and Whole−Rock Compositions[J/OL]. Northwestern Geology, 2023: 1−16. doi: 10.12401/j.nwg.2023155
马忠元,刘光莲,汪周鑫,等. 东昆仑西段它温查汉西铁金多金属矿成矿模式探讨[J]. 新疆地质,2022,40(3):362–367. doi: 10.3969/j.issn.1000-8845.2022.03.007 MA Zhongyuan, LIU Guanglian, WANG Zhouxin, et al. Discussion on the Metallogenic Model of Tawenchahanxi Iron Polymetallic Deposit in the Western Part of East Kunlun Mountains[J]. Xinjiang Geology, 2022, 40(3): 362-367. doi: 10.3969/j.issn.1000-8845.2022.03.007
潘力, 何青林, 陈康, 等. 利用重磁资料解译川西地区深层断裂构造及预测火山岩顶界面深度[J]. 成都理工大学学报(自然科学版), 2023, 50(2): 240−248. PAN Li, HE Qinglin, CHEN Kang, et al. Application of gravity and magnetic data to the interpretation of deep fault structures and prediction of volcanic rock top interface depth in the western Sichuan Basin, China [J], Journal of Chengdu University of Technology (Science & Technology Edition), 2023, 50(2): 240−248.
乔耿彪,伍跃中. 东昆仑祁漫塔格地区花岗岩成因类型对成矿作用的控制[J]. 西北地质,2010,43(4):134–142. doi: 10.3969/j.issn.1009-6248.2010.04.016 QIAO Gengbiao, WU Yuezhong. Genetic Types of Granite Controlled the Mineralization in Qimantage Area, Eastern Kunlun [J], Northwestern Geology, 2010, 43(4): 134-142. doi: 10.3969/j.issn.1009-6248.2010.04.016
田承盛,丰成友,李军红,等. 青海它温查汉铁多金属矿床40Ar-39Ar年代学研究及意义[J]. 岩石矿物学杂志,2013,32(1):169–176. TIAN Chengsheng, FENG Chengyou, LI Junhong, et al. 40Ar-39Ar geochronology of Tawenchahan Fe-polymetallic deposit in qimantag Mountain of qinghai Province and its geological implications[J]. Mineral Deposits, 2013, 32(1): 169-176.
王谦身,滕吉文,张永谦,等. 鄂尔多斯-中秦岭-四川东部的重力异常场与深部地壳结构[J]. 地球物理学报,2015,58(2):532–541. doi: 10.6038/cjg20150216 WANG Qianshen, TENG Jiwen, ZHANG Yongqian, et al. Gravity anomalies and deep crustal structure of the Ordos basin-middle Qinling orogen-eastern Sichuan basin[J]. Chinese J. Geophys. (in Chinese), 2015, 58(2): 532-541. doi: 10.6038/cjg20150216
张文权,张爱奎,孟军海,等. 高精度磁测反演技术在沙丘地区找矿中的应用[J]. 西北地质,2012,45(1):277–282. doi: 10.3969/j.issn.1009-6248.2012.01.035 ZHANG Wenquan, ZHANG Aikui, MENG Junhai, et al. Application Inversion technology into Dune with High-Precision Magnetic Survey [J]. Northwestern Geology, 2012, 45(1): 277-282. doi: 10.3969/j.issn.1009-6248.2012.01.035
赵政璋, 李永铁, 叶和飞, 等. 青藏高原大地构造特征及盆地演化[M]. 北京: 科学出版社, 2001 ZHAO Zhengzhang, LI Yongtie, YE Hefei, et al. Tectonic characteristics and basin evolution of Qinghai-Tibet Plateau[M]. Beijing: Science Publishing House, 2001.
计量
- 文章访问数: 94
- HTML全文浏览量: 9
- PDF下载量: 52
