图 1.
试验区三维地形概览
Figure 1.
Three-dimensional terrain overview of the test area
| China Geological Survey Chinese Academy of Geological Sciences | Host |
| 科学出版社 | Publish |
| Citation: |
ZHANG Jinming, CHEN Guangting, CAI Hangjia, TIAN Chengxiu, LEI Xiaoqing. 2023. Geochemical characteristics and zircon U-Pb age of gabbros in the Zhaqiaohe ophiolite mélange, and its limitation on the ocean ridge environment[J]. Geology in China, 50(6): 1837-1847. doi: 10.12029/gc20200416002
|
Geochemical characteristics and zircon U-Pb age of gabbros in the Zhaqiaohe ophiolite mélange, and its limitation on the ocean ridge environment
-
Abstract
This paper is the result of geological survey engineering.
Objective Zhaqiaohe ophiolite is located at the junction of the Middle Qilian and South Qilian tectonic belt. The ophiolite in the belt is composed of serpentinite, serpentinized pyroxene peridotite, gabbro, diabase, block bedded altered basalt, a small amount of pillow basalt and lenticular silicite. It is the west extension of the Lachishan ophiolite melange belt in Qinghai Province.
Methods In this paper, the study on the chronology and petrogeochemistry of gabbro in Zhaqiaohe ophiolite has been carried out.
Results According to the characteristic judgment of the main chemical components of the gabbro, such as SiO2, Na2O, K2O, (Fe2O3 + FeO) and MgO, the gabbro is tholeiitic series, belonging to tholeiite. According to the REE distribution pattern, there is no Eu anomaly or weak negative anomaly, indicating that the rock experienced the separation and crystallization of plagioclase in the melting process, and the REE distribution curve is approximately flat; the trace elements have obvious Ta and Sr positive anomaly characteristics, indicating that gabbro in ophiolite in this area has similar characteristics with E-MORB. The results of LA-ICP-MS zircon U-Pb dating show that the gabbro crystallization age is (524.0±2.5) Ma.
Conclusions It is suggested that these gabbros originated from the initial slowness without large-scale melting. It is suggested that the Zhaqiaohe ophiolite may have formed in the oceanic ridge environment during the initial cracking stage of the paleo ocean basin.
-
Keywords:
- ophiolite belt /
- gabbro /
- zircon U-Pb dating /
- E-MORB /
- Central Qilian /
- geological survey engineering /
- Zhaqiaohe /
- Qinghai Province
-
0. 引言
高陡边坡上的危岩区域是水利水电工程建设中常见工程问题之一，其稳定状况事关水利枢纽工程建设和运行安全[1 − 2]。对边坡岩体进行地质勘察，获取边坡岩体的地质参数，是危岩体研究的首要工作。传统的地质勘察主要采用的是人工接触式方法，例如用皮尺和地质罗盘测量[3 − 5]。随着建设工程规模的不断增大，高陡边坡地质环境变得更加险峻和复杂，由于调查人员的可及范围有限，地质勘察工作的展开不仅效率低下，而且获取的地质参数不完整。传统的接触式调查已经难以满足实际工程的高效化和数字化要求。
在过去的20年中，许多研究人员采用了非接触式数据采集方法来解决这一问题[6]。Bitenc等[7]将优化后的激光雷达测量技术（light laser detection and ranging，LiDAR）应用于边坡地形测量数据采集，提高了边坡数据采集的精度；Salvini等[8]结合直升机摄影测量技术和地面激光扫描技术（terrestrial laser scanning，TLS）进行了边坡岩石块体和节理系统的详细地理信息绘制；Kumhálová等[9]，贾虎军等[10]通过机载激光扫描和倾斜摄影技术，获取精准的地质灾害隐患点、精细的地形地貌勘测数据，解决了设备难以进入的陡峭地形模型数据采集的问题。然而，激光扫描设备即便借助直升机等多地形工具，在进行高差过大，时限较短的数据采集任务时仍然显得笨重，难以满足需求。
近年来，为了满足快速获取复杂边坡表面特征的需求，成本更低、灵活方便、不受地形限制的无人机得到了一定的应用，它可以解决复杂边坡死角模型难以采集的问题。Niethammer等[11]，Riquelme等[12]，Zeybek等[13]采用无线电控制小型四旋翼无人机拍摄完整边坡滑坡的图像，进一步建立高分辨率的数字地形模型，实现了更快速获取坡面点云数据；赵婷婷等[14]，梁京涛等[15]，黄发明等[16]通过无人机倾斜摄影技术（unmanned aerial vehicle oblique photography technology，UAVOP）与多图像重建技术（structure from motion，SfM）采集复杂边坡表面点云模型，实现了对紧急，复杂的边坡数据采集任务的快速反应。基于无人机的数字摄影测量技术已经足以表现复杂边坡表面的详细形态，但这些研究对精确模型的利用率仍然较低，尚未充分地通过精确点云对其所反映的边坡隐蔽、细微的安全隐患（如危岩体）进行有效的勘察。
采用无人机开展高陡边坡的安全勘察研究时，由于数据采集方案单一，设备效率、模型精度较低或高精度模型的利用率较低[17 − 19]，尤其在高陡边坡中，由于畸变问题，会使得模型结果与实际情况产生一些偏差。为解决这一问题，本研究提出了一种基于无人机点云的超500 m高位边坡孤立单平面滑移式危岩体半自动提取的新方法，实现了以下三项关键的目标:
（1）基于免像控的无人机贴近摄影测量技术，设计了一套高陡边坡表面点云模型的测量方法，并构建了三维坐标精度为±10 cm、影像分辨率为1 cm的试验区精细边坡三维点云模型。
（2）基于高陡边坡点云模型提取了异型滑移式块体。利用异型滑移式块体突出于边坡表面的细微特征，手动分割确定了异型滑移式块体的位置和其具体形态。
（3）根据异型滑移式块体模型边界点云的密度特征拟合了后壁平面。此外，自动计算了异型滑移式块体模型的各类特征值，并进一步得出稳定性系数（K）以确定块体中的危岩体。最后，结合危岩体的方量和最大高差，对危险岩体进行了危害性评价。
1. 研究区概况
玉龙喀什水利枢纽工程（图1）是新疆和田玉龙喀什河山区河段的控制性水利枢纽工程，水库设计最大坝高230.5 m，总库容5.36×108 m3，工程主要有最大坝高233.5 m的面板堆石坝、“表、中、深”泄洪隧洞以及引水发电系统等组成，属大（2）型Ⅱ等工程。
试验区为工程左岸坝址区的高位边坡，最大高差超过500 m。边坡地形陡峭，坡面植被不发育，山体的表面覆盖有大量灰尘，边坡表面并未出露明显的结构面。考虑施工场地狭窄、施工周期长，坝基边坡开挖对上部高位边坡扰动大，而坝址区两侧自然高位边坡存在大量危岩体和不良地质块体，在施工期构成重大安全隐患。亟须对坝址区超高位边坡进行危岩调查和稳定性分析。
经过人工初步调查，勘察到可能的危岩体多为平面滑移式危岩体（见图2）。岩体突出于边坡表面，且与边坡斜面之间存在有明显的连接边缘。本研究对试验区表面此类孤立的单平面滑移式危岩体进行进一步的调查。
2. 技术方法
针对试验区高边坡广泛分布平面滑移式危岩体的情况，本文设计了一套从边坡点云模型数据采集到危岩体半自动识别及评价的边坡危岩体完整的地质勘察方法（图3）。程序主要由Python语言编写。该识别方法包括以下四个步骤：
（1）数据采集阶段：采用无人机贴近摄影测量结合精细的仿地飞行航路设计、智能建模等多种技术手段，获取高边坡精细地表LAS点云。
（2）孤立异型滑移式块体点云模型提取阶段：依据边坡突出块体边界特征将孤立的异型滑移式块体点云模型分别手动分割出来，并对模型进行点云下采样。
（3）滑移破坏后壁平面提取阶段：应用点云的密度特征自动提取突出块体与边坡连接的边界，此边界可近似为该块体的滑移破坏的平面的边界。
（4）孤立危岩体定性与危害性评估阶段：应用最小二乘法将边界点云拟合到后壁平面上，计算后壁倾角[20]。使用原块体点云和对后壁投影点之间的欧式距离计算最大高度差。利用Stokes公式和Alpha Shape算法确定投影点的平面边界，计算后壁平面的面积和块体点云体积。以上特征参数将用于定性孤立的危岩体（依据稳定性系数）。对于已经定性的危岩体，依据其体积与最大高差依次对所有危岩体的危害性进行进一步的评价。
2.1 数据采集及异型滑移式块体识别
本文使用的数据是新疆玉龙喀什水利工程左岸斜坡的有序点云数据，考虑到试验区山、谷、沟、梁交错的特殊地形，本文采用无人机贴近摄影测量[21 − 23]和全球导航卫星系统（global navigation satellite system, GNSS）实时动态差分定位技术[24]相结合的方法，综合利用免像控的无人机摄影测量技术、精细化地面模拟航路设计、智能建模等多种技术手段，获取500 m以上边坡的精细化地表信息。
本研究使用大疆M300 RTK无人机设备搭载Zenmuse P1全画幅云台相机（见图4）。无人机集成了双GNSS导航定位设备，通过集成的网络RTK模块，获取附近CORS站提供的实时动态差分信息，可为无人机提供实时厘米级定位数据，满足了免像控标准。
在对试验区模型采集的应用中，本研究设计的基于精细仿地飞行路线的贴近摄影测量技术方法主要采用三种综合航摄方案。
（1）初始飞行和航线规划（见图5）：采用航天飞机雷达地形测绘（shuttle radar topography mission, SRTM）采集30 m分辨率数字高程模型（digital elevation model, DEM），初始飞行航路较高。通过空中三角测量，可以获得调查区域的图像，并生成初始数字地表模型（digital surface model, DSM）[25]。
图 5. 第一张正射影照片路线规划（分辨率2.5 cm）Figure 5. Route planning of the first orthophoto photograph (resolution 2.5 cm)（2）二次精细仿地飞行规划与飞行（见图6）：基于初始模型DSM高程数据，规划精细地面模拟飞行路线，并为该路线建立多组使用精细仿地飞行的镜头拍摄角度。为了获得高坡度倾斜摄影测量图像，对近似直立的坡面使用了贴近摄影测量技术。通过无人机三角测量和建模，生成测量区倾斜摄影测量OSGB模型[26]。考虑到试验区边坡沟梁相间的地形，精细的仿地飞行路线沿地形梯度呈上下方向设置。
（3）三次手控飞行：检查OSGB模型存在的空洞和纹理拉伸的区域。如果出现这种情况，必须对关键区域进行额外的手动飞行拍摄。补拍完成后，再次进行空中三角测量和建模，生成测量区域完整的真实场景，即为最终的三维模型。针对研究区倒悬岩体等特殊场景，采用−10°等特殊视角和手控飞行的拍摄方案。
DJI Maps被用于无人机三维模型建模，采用基于计算统一设备架构（compute unified device architecture, CUDA）的智能重构算法[27]，并将场景图文件格式模型（open scene graph binary, OSGB）转换为二进制LAS点云模型（见图7）。
图 7. 无人机三维模型建模并生成点云模型注：（a）为构造三角网TIN；（b）为生成模型白膜；（c）为模型纹理映射；（d）LAS点云格式转换Figure 7. UAV three-dimensional modeling and point cloud model generation异型滑移式块体在边坡表面表现为一块完整的突出块体。对于边坡的点云模型，其与边坡接触处通常为一圈曲率较大的闭合点云。依据此特征在软件CloudCompare内使用点云切分工具将孤立的异型滑移式块体点云模型进行提取。此外，为了减小后续的计算量，使用滤波工具进一步进行点云下采样。
2.2 块体失稳平面自动提取
2.2.1 后壁平面边界点云
本研究认为，异型滑移式块体点云的边界近似处于滑移的后壁平面上，依据边界点云的密度与密度变化率特征进行后壁平面边界的提取。对于表面均匀的点云模型来说，处于模型边缘的点云一般密度（在某范围内的点云个数）相对较小[28 − 30]。此外，该处的点云在某个范围内密度会明显的下降。
对于单独的异型滑移式块体点云，构造点云的K邻近（K-nearest neighbor, KNN）搜索结构，对每个点查找其搜索范围内的点云密度。对于点P1，当该点点云密度较小，在R1范围内的点数为8（该点的密度）或者在一定的范围R2内，存在密度远超于它的点P2（密度为19）时，本研究认为其为一枚后壁平面的近似边界点（见图8）。
2.2.2 后壁平面拟合
异型滑移式块体边界点云可近似拟合于空间中的一块平面。考虑到试验区表面被大量的灰尘覆盖，边坡模型并未出露明显的结构面，则该平面被认为是当前异型滑移式块体可能的滑动结构面（后壁平面）。图9（a）中蓝色的点云为通过边界密度分析得到的后壁平面边界点，灰色点为该块体的其他点云。利用最小二乘法将边界点拟合为图9（b）中的黑色网格平面，该平面将为异型滑移式块体点云模型的各种特征计算提供依据。
2.3 孤立危岩体定性分析
2.3.1 三维特征估算
将异型滑移式块体点云整体向后壁平面进行投影，并对其进一步进行Alpha Shape边界搜索[31]，得到平面的边界点索引。该索引方向对应的边界点逆时针排列。对边界点利用Stokes的三维有向平面面积积分公式进行积分[32]，得到后壁平面的面积。后壁平面多边形
A 有N 个顶点P1,P2,⋯,Pi,⋯,Pn 呈逆时针方向排列，形成一个环L 顶点坐标表示为Pi=(Xi,Yi,Zi) ，如图10（a）所示。可求得其面积为SA=∬ 。
图 10. 异型滑移式块体的后壁面积与块体体积计算Figure 10. Calculation for the rear wall area and block volume ofprofiled block面积由Stokes公式求得：
\oint_LP\mathrm{d}_x+Q\mathrm{d}_y+R\mathrm{d}_{\text{z}}=\iint_A^{ }\left[\begin{gathered}\left(\frac{\partial R}{\partial y}-\frac{\partial Q}{\partial\text{z}}\right)\cos\left(n,x\right) \\ +\left(\frac{\partial P}{\partial\text{z}}-\frac{\partial R}{\partial x}\right)\cos\left(n,y\right) \\ +\left(\frac{\partial Q}{\partial x}-\frac{\partial P}{\partial y}\right)\cos\left(n,\text{z}\right) \\ \end{gathered}\right]\mathrm{d}_s (1) 式中：
L ——后壁平面多边形边界形成的环；P 、Q 、R ——代表空间中某一点依赖于x、y、z坐标 的向量场三个分量函数；\mathrm{d}_x 、\mathrm{d}_y 、\mathrm{d}_{\text{z}} ——表示线积分中的微分元素；A ——后壁平面多边形；\dfrac{\mathrm{\partial}R}{\mathrm{\partial}y}-\dfrac{\mathrm{\partial}Q}{\mathrm{\partial}\text{z}} ，\dfrac{{\partial P}}{{\partial {\textit{z}}}} - \dfrac{{\partial R}}{{\partial x}} ，\dfrac{{\partial Q}}{{\partial x}} - \dfrac{{\partial P}}{{\partial y}} ——旋度向量分别 指向x、y、z轴 的分量；\left( {n,x} \right) 、\left( {n,y} \right) 、\left( {n,{\textit{z}}} \right) ——多边形的单位法向量n与坐 标轴 x、y、z的夹角；\mathrm{d}_s ——平面上的微小面积元素。对完整的异型滑移式块体点云应用Alpha Shape算法得到异型滑移式块体点云的体积，如图10（b）所示。后壁平面与水平面之间的夹角为后壁倾角[33]，后壁倾角由公式（2）计算，可记为
\theta ，如图11（a）所示。后壁面与孤立岩体所有点之间的最大垂直距离为最大高差，见图11（b）。孤立岩体的最大高差可依据公式（3）进行计算。
图 11. 异型滑移式块体的后壁倾角与最大高差计算Figure 11. Calculation of the rear wall inclination angle and maximum height difference of profiled blocks\theta = {\cos ^{ - 1}}\left( {\frac{{\vec n \cdot \vec v }}{{\vec {\left| n \right|} \cdot \vec {\left| v \right|} }}} \right) (2) 式中：
\theta ——后壁面倾角；\vec n ——后壁平面的法向量；\vec v ——水平面的法向量。\Delta H = \max \left\{ {\left| {\overrightarrow {{A_1}{B_1}} } \right|,\left| {\overrightarrow {{A_2}{B_2}} } \right|,\cdots,\left| {\overrightarrow {{A_n}{B_n}} } \right|} \right\} (3) 式中：
\Delta H ——孤立岩体最大高差；{A_i} ——第i 个对应的投影点云；{B_i} ——孤立岩体的第i 个原点云。2.3.2 危岩体定性
得到了岩体点云模型的各种空间特征后，要判断提取的异型滑移式块体是否为危险岩体[34 − 35]，则应结合稳定系数（K）进行定性判断。稳定系数是定性判断沿后壁平面滑移破坏的危岩体的常用依据[36]。其公式已广泛应用于单滑移面危岩体稳定性试验领域。利用《中华人民共和国新能源行业标准》中《水电工程危险岩体工程地质调查与治理规范》（以下简称《规范》）对K进行计算和评价，如表1。
表 1. 单平面滑动岩体稳定性评价Table 1. Stability evaluation of single-plane sliding rock mass稳定性系数 稳定性分级 
稳定 
基本稳定 
欠稳定 
不稳定 | Show Table
DownLoad:
CSV
对于每一个可能的危险岩体对象，自动依次计算其对应的后壁面积、后壁倾角和块体方量，并对结果进行综合的稳定性分析。稳定性系数可依据式（4）计算。
K = \frac{{Q\cos \alpha \tan \varphi + cS}}{{Q\sin \alpha }} (4) 式中：
Q ——异型滑移式块体的重量/N；\alpha ——后壁的倾角/（°）；\varphi ——后壁滑动面的内摩擦角/（°）；c ——黏聚力/MPa；S ——异型滑移式块体的滑动面积/m2。2.4 危岩体危害性评估
应用安全系数判断异型滑移式块体是否为危岩体后，还需要综合评价危岩体方量和最大高差的特征[37]，以确定这些危岩体的危险等级和危害性，其中《规范》提供的危岩体规模分类依据如表2所示。
表 2. 水电工程危险岩体规模分级Table 2. Scale classification of dangerous rock mass in hydropower projects评价依据 小型 中型 大型 超大型 体积/m3 
| Show TableDownLoad:
CSV
3. 结果与讨论
3.1 数据采集与异型滑移式块体识别
通过精细化地面模拟航路设计、无人机贴近摄影测量、多图像三维智能建模等步骤重建获得玉龙喀什水利工程左岸500 m以上边坡的高精度三维点云模型信息（图12）。考虑到玉龙喀什水利枢纽孤立岩体多为异型滑移式块体，研究以左岸1块危岩体密集的区域进行危岩体半自动提取试验。
图 12. 左岸整体点云，左岸样本区区域Figure 12. Overall point cloud on the left bank and sample area on the left bank在软件CloudCompare内依据异型滑移式块体与边坡面的接触边界特征将该样本区可能的异型滑移式块体点云模型切分出来，如图13（a）。
共提取出5块异型滑移式块体（L1，L2，L3，L4，L5）见图13（b），所有的块体均具备明显的边界特征。此外，为了减少计算成本，对点云进行下采样。图14（b）为块体L1下采样后的点云（采样前点云个数:10486631，采样后点云个数：1171753）。
3.2 精确特征提取与危岩体定性
以块体L1为例，对于切分出的异型滑移式块体点云模型，使用点云密度与密度变化特征提取后壁平面边缘点云见图15（b），这类点可近似处地于后壁平面上。
图 15. 依据点云密度提取的异型滑移式块体后壁平面边缘点云Figure 15. Edge point cloud of the rear wall plane of profiled blocks extracted based on point cloud density图15（b）中蓝色点为通过点云密度分析所提取的后壁平面边缘点云，其表现为一圈闭合的边界。而灰色点为块体点云的其他部分。五块异型滑移式块体点云边界准确模型见图15（a）。利用最小二乘法将点云拟合为平面方程，此平面为岩体破坏时可能的滑动面（见图16）。
根据《规范》，应确定异形滑块的后壁倾斜度、后壁面积与方量三项特征，以判断和评价危岩体。以异型滑移式块体L1为例，对L1点云边界应用最小二乘法得到平面方程，由式（4）计算后壁面与水平面的法向量夹角，即后壁面倾角。将L1点云垂直投影到后壁上，计算原始点与投影点之间的欧氏距离。最大距离对应于L1的最大高差（见图17）。
图 17. 异型滑移式块体L1的后壁平面和点云投影Figure 17. Rear wall plane and point cloud projection of profiled block L1L1的面积和体积可以使用Alpha Shape算法和Stokes积分计算。图18（a）中的红色点云表示通过Alpha Shape算法提取的异形滑动块的后壁有向的精细边界。该边界将用于计算后壁面积。此外，图18（b）描绘了异形滑块外表面的三角形包络网络，将用于计算块体的体积。
图 18. Alpha Shape算法分别计算异型滑移式块体的后壁面积与体积Figure 18. Alpha Shape algorithm for separate calculation of the rear wall area and volume of profiled block如表3所示，采用提出的方法确定各异形滑块的后壁倾角、后壁面积、岩石体积和最大高差。
表 3. 异形滑移式块体特征数据统计Table 3. Statistical analysis of characteristic data for profiled blocks块体编号 后壁倾角/（°） 后壁面积/m2 块体体积/m3 最大高差/m L1 50.35 728.15 1712.30 7.65 L2 45.58 835.53 2398.60 10.60 L3 40.42 842.94 2299.47 9.48 L4 26.23 706.89 3294.37 14.53 L5 46.61 886.78 1690.18 7.71 | Show TableDownLoad:
CSV
通过工作人员的现场调查，可获得现场大部分岩质边坡危岩体以及后壁平面的物理力学参数（如表4）。进一步，对所有识别出的异形滑块进行稳定性K计算，从而判断其是否为危险岩体（如表5），可以看出，除了L5以外均为危岩体。
表 4. 异形滑移式块体物理力学参数Table 4. Physical and mechanical parameters of profiled block参数 块体重度/（kN·m−3） 结构面黏聚力/MPa 结构面内摩擦角/（°） 取值 25.8 0.100 35 | Show TableDownLoad:
CSV
表 5. 危险岩体稳定系数的计算与定性Table 5. Calculation and qualitative characterization of stability coefficient of dangerous rock mass块体编号 稳定性系数 稳定性分析 块体性质 L1 0.58 不稳定 是 L2 0.69 不稳定 是 L3 0.82 不稳定 是 L4 1.42 稳定 否 L5 0.66 不稳定 是 | Show TableDownLoad:
CSV
3.3 危岩体综合评价
在本研究中，在提取了试验区的危岩体后，还需要对危险岩体的危险性等级进行确定和评价。孤立岩体的稳定性和危险等级应结合稳定性系数、岩体体积和最大高差来确定。为了定量评价危险岩体的危险性，结合危险岩体的体积和高差特征对《规范》进行了评价。
如表6所示，本研究在试验区发现了一系列小山梁状危岩体，体积均在2000 m3左右，属大型的单平面滑动危岩体，危岩体最大高差在7～15 m，其中L4高差为14.53 m。经过现场调查，这些危险岩体已被爆破拆除。
表 6. 危险岩体与岩体体积评价Table 6. dangerous rock mass and rock massvolume evaluation危岩体编号 危岩体体积/m3 危岩体体积评价 L1 1712.30 大型 L2 2398.60 大型 L3 2299.47 大型 L5 1690.18 大型 | Show TableDownLoad:
CSV
3.4 适用性讨论
（1）方案没有考虑详细的地质参数。由于本研究的重点是快速识别危险岩体并提供初步评价，因此没有将岩体的具体地质特征作为依据。在未来的研究中，将自动识别过程与详细的地质参数相结合，实现对危险岩体更精细的评价。
（2）本研究只考虑了危岩体的单平面滑动的破坏模式，没有对其他滑动模型进行研究。这主要是由于研究区内普遍存在沿单一平面滑动的岩体。目前提出的半自动识别模型仅能提取沿后壁平面滑动的岩体。在未来的研究中，该模型将进一步扩展，以包含识别具有各种失稳和破坏模式的危岩体。
4. 结论
本研究基于无人机贴近摄影测量技术与点云模型数据处理算法，对玉龙喀什水利枢纽工程左岸边坡一块岩体密集区进行了初步的快速地质勘察，实现了试验区边坡点云模型重建、异型滑移式块体提取、稳定性分析等目标。该方法在快速获取高精度边坡数据、智能提取三维异型滑移式块体、准确计算三维特征、评估危险岩体等方面相对于传统方法具有一定的优势。本文提供了从数字边坡模型获取到危险岩体识别和初步评价的解决方法。为边坡工程数字化奠定一定的研究基础。研究得出了如下的主要结论:
（1）结合无人机贴近摄影测量、GNSS实时动态差分定位技术，以及自主设计的三次精细综合无人机飞行路线规划和手控校验，可为航空三角测量和建模提供分辨率优于1 cm的图像，解决了由于高差较大而产生模型畸变的问题。
（2）利用最小二乘法和Alpha Shape计算异形滑块的后壁倾角、后壁面积、方量等特征数据，确定异形滑块的稳定系数。从计算结果可以看出，试验区大部分块体为危险岩体（K<1.05），只有1个块体为稳定的异型滑移式块体（L4:K=1.42），解决了人工调查危岩精度较低的问题，与现场稳定性评价结果一致。
（3）对识别出的异形滑块进行特征计算，可快速初步确定危险岩体的危险等级。采用高差和体积相结合的方法对危险岩体进行评价。其中危害性最大的是一块大型危岩体（L2），高差为10.60 m，体积为2398.60 m3，对危岩体的提取与准确的评价为边坡工程的数字化奠定了一定的基础。
-
References
Chen Xuanhua, Shao Zhaogang, Xiong Xiaosong, Gao Rui, Liu Xuejun, Wang Caifu, Li Bing, Wang Zengzhen, Zhang Yiping. 2019. Fault system, deep structure and tectonic evolution of the Qilian Orogenic Belt, Northwest China[J]. Geology in China, 46(5): 995-1020 (in Chinese with English abstract). Coleman R G. 1977. Ophiolites: Ancient oceanic lithosphere?[J]. Mineralogical Magazine, 42(322): 308-308. Feng Yimin, Cao Xuanduo, Zhang Erpeng, Hu Yunxu, Yang Junlu, Pan Xiaoping, Jia Qunzi, Li Wenming. 2002. Structural Orogenic Process and Dynamics of the West Qinling Orogenic Belt[M]. Xi'an: Xi'an Map Publishing House, 263(in Chinese with English abstract). Feng Yimin, He Shiping. 1996. Geotectonics and Orogeny in Qilian Mountains[M]. Beijing: Geological Publishing House, 1-135(in Chinese with English abstract). Fu Changlei, Yan Zhen, Guo Xiaoqin, Niu Manlan. 2014. Geochemistry and SHRIMP zircon U-Pb age of diabases in the Lajishankou ophiolitic melange, South Qilian terrane[J]. Actor Petrologica Sinica, 30(6): 1695-1706 Hou Qingye, Zhang Hongfei, Zhang Benren, Zhao Zhidan, Zhu Yunhai. 2005. Paleogeochronicity of the Laji Mountain in the middle Qilian orogenic belt and its attribution: Geochemical evidence from basic volcanic rocks[J]. Earth Science, 30(1): 61-70(in Chinese with English abstract). Li Chunyu, Liu Yangwen, Zhu Baoqing, Feng Yimin, Wu Hanquan. 1978. Tectonic History of Qinling and Qilian Mountains[M]. Beijing: Geological Publishing House, 174-187(in Chinese with English abstract). Li Jianfeng, Zhang Zhicheng, Han Baofu. 2010. Geochronolgy and geochemistry of Early Paleozoic plutons from Subei and Shibaocheng areas, the western segment of Central Qilian and their geological implications[J]. Acta Petrologica Sinica, 26(8): 2431-2344(in Chinese with English abstract). Li Huaikun, Geng Jianzhen, Hao Shuang, Zhang Yongqing, Li Huiming. 2009. Determination of zircon U-Pb isotope ages by laser Ablation Multi-receiver Plasma Mass Spectrometer (LA-MC-ICPMS)[J]. Acta Mineralogica Sinica, 29(S1): 600-601 (in Chinese with English abstract). Middlemost E A K. 1994. Naming materials in the magma/igneous rock system[J]. Earth-science Reviews, 37(3-4): 215-224. doi: 10.1016/0012-8252(94)90029-9 Miyashiro A, Shido F. 1975. Tholeiitic and calc-alkalic series in relation to the behaviors of titanium, vanadium, chromium, and nickel[J]. American Journal of Science, 275 (3): 265-277. doi: 10.2475/ajs.275.3.265 Nicolas A. 1989. Structures of Ophiolites and Dynamics of Oceanic Lithosphere[M]. Dordrecht: Kluwer. 1-367. Pearce J A, Lippard S J, Roberts S. 1984. Characteristics and tectonic significance of super-subduction zone ophiolites[J]. Geological Society, London, Special Publications, 16(1): 74-94. Qiu Jiaxiang, Zeng Guangce, Wang Siyuan, Zhu Yunhai. 1997. Early Paleozoic Marine Volcanic Rocks and Mineralization in Laji Mountain[M]. Wuhan: China University of Geosciences Press, 1-118(in Chinese with English abstract). Sun S S, McDonough W F. 1989. Chemical and isotopic systematics of oceanic basalts: Implications for mantle composition and processes[J]. Geological Society, London, Special Publications, 42(1): 313-345 doi: 10.1144/GSL.SP.1989.042.01.19 Tseng C Y, Yang H J, Yang H Y, Liu D Y, Tsai C L, Wu H Q, Zuo G C. 2007. The Dongcaohe ophiolite from the North Qilian Mountains: A fossil oceanic crust of the Paleo-Qilian Ocean[J]. Chinese Science Bulletin, 52(7): 2390-2401 Wang Erqi, Zhang Qi, Burchfiel C B. 2000. Laji Mountain, Qinghai Province: A tectonic window of multi-stage uplift [J]. Chinese Journal of Geology, 35(4): 393-500(in Chinese with English abstract). Wang Jinrong, Zhang Xinhu, Zhang Xiang, Hui Xinwei, Zhang Benqi, Zhao Lei. 2021. Rock association, genesis and tectonic significance of Mustangshan structure in the western part of Middle Qilian[C]//Abstracts of the National Symposium on Petrology and Geodynamics, 466-468(in Chinese with English abstract). Wood D A. 1980. The application of a Th-Hf-Ta diagram to problems of tectonomagmatic classification and to establishing the nature of crustal contamination of basaltic lavas of the British Tertiary Volcanic Yrovince[J]. Earth and Planetary Science Letters, 50(1): 1-30. doi: 10.1016/0012-821X(80)90115-6 Xia Linqi, Xia Zuchun, Ren Youxiang, Peng Ligui, Zhang Cheng, Yang Jinghua, Wang Xingan, Li Zhipei, Han Song, Huang Zhongxiang. 1991. Marine Volcanic Kocks from Qilian and Qinling Mountains[M]. Wuhan: China University of Geosciences Press, 1-304 (in Chinese with English abstract) Xiao Xuechang, Chen Guoming, Zhu Zhizhi. 1978. Geological tectonic significance of the ancient ophiolite belt in Qilian Mountains[J]. Acta Geologica Sinica, 52(4): 281-295 (in Chinese with English abstract). Xiao W, Windley B F, Yong Y, Yan Z, Yuan C, Liu C, Li J. 2009. Early Paleozoic to Devonian multiple-accretionary model for the Qilian Shan, NW China[J]. Journal of Asian Earth Sciences, 35: 323-333. doi: 10.1016/j.jseaes.2008.10.001 Xu Zhiqin, Xu Huifen, Zhang Jianxin, Li Haibing, Zhu Zhizhi, Qu Jingchuan, Chen Daizhang, Chen Jinlu, Yang Kaichun. 1994. The proliferative terrane of Caledonian subduction complex in Nanshan of North Qilian Corridor and its dynamics[J]. Acta Geologica Sinica, 68 (1): 1-15 (in Chinese with English abstract). Yan Zhen, Wang Zongqi, Li Jiliang, Xu Zhiqin, Deng Jinfu. 2012. Tectonic properties and accretionary orogeny of West Qinling wedge[J]. Acta Petrologica Sinica, 28(6): 1808-1828(in Chinese with English abstract) Zhang Jianxin, Xu Zhiqin, Chen Wen, Xu Huifen. 1995. Discussion on the age of subduction-accretion complex/volcanic arc in the middle part of North Qilian [J]. Journal of Rock and Minerals, 16(2): 112-119 (in Chinese with English abstract) Zhang Xueting, Yang Shengde. 2007. Qinghai Province Plate Tectonics Research 1/1 Million Qinghai Province Geotectonic Map Description [M]. Beijing: Geological Publishing House, 1-221 (in Chinese with English abstract). Zhao Wenjin, Mechie J, Feng Mei, Shi Danian, Xue Guangqi, Shu Heping, Song Yang, Yang Hongwei, Liu Zhiwei. 2014. Orogeny in Qilian Mountains and the characteristic structure of the lithospheric mantle[J]. Geology in China, 41(5): 1411-1423(in Chinese with English abstract). Zheng Ying, Chen Guangting, Zhang Xiaoyong, Song Taizhong, Liang Kunxian. 2017. Geochemical characteristics and chronological significance of Ordovician quartz diorite in Suli area, Middle Qilian [J]. Northwest Geology, 50(4): 588-597 (in Chinese with English abstract). Zhong Mingjie. 1994. Discovery of Paleozoic in Laji Mountain, Qinghai Province[J]. Geological Review, 22(5): 347-347 (in Chinese with English abstract). Zuo Guochao, Zhang Shuling, Cheng Jiansheng, Gong Yanxue, Wang Yanbin, Wu Hanquan. 1996. Division of ophiolite belt in Qilian area and its tectonic significance [C]//Proceedings of the Symposium on Ophiolite and Geodynamics, 129-134(in Chinese with English abstract). 陈宣华, 邵兆刚, 熊小松, 高锐, 刘雪军, 王财富, 李冰, 王增振, 张义平. 2019. 祁连造山带断裂构造体系、深部结构与构造演化[J]. 中国地质, 46(5): 995-1020. 冯益民, 曹宣铎, 张二朋, 胡云绪, 杨军禄, 潘晓萍, 贾群子, 李文明. 2002. 西秦岭造山带结构造山过程及动力学[M]. 西安: 西安地图出版社, 263. 冯益民, 何世平. 1996. 祁连山大地构造与造山作用[M]. 北京: 地质出版社, 1-135. 付长垒, 闰臻, 郭现轻, 牛漫兰. 2014. 拉脊山口蛇绿混杂岩中辉绿岩的地球化学特征及SHRIMP锆石U-Pb年龄[J]. 岩石学报, 30(6): 1695-1706. 侯青叶, 张宏飞, 张本仁, 赵志丹, 朱云海. 2005. 祁连造山带中部拉脊山古地慢特征及其归属: 来自基性火山岩的地球化学证据[J]. 地球科学, 30(1): 61-70. 李春昱, 刘仰文, 朱宝清, 冯益民, 吴汉泉. 1978. 秦岭及祁连山构造发展史[M]. 北京: 地质出版社, 174-187. 李建锋, 张志诚, 韩宝福. 2010. 中祁连西段肃北、石包城地区早古生代花岗岩年代学、地球化学特征及其地质义[J]. 岩石学报, 26(8): 2431-2444. 李怀坤, 耿建珍, 郝爽, 张永清, 李惠民. 2009. 用激光烧蚀多接收器等离子体质谱仪(LA-MC-ICPMS)测定锆石U-Pb同位素年龄的研究[J]. 矿物学报, 29(增刊): 600-601. 丘家骧, 曾广策, 王思源, 朱云海. 1997. 拉脊山早古生代海相火山岩与成矿[M]. 武汉: 中国地质大学出版社, 1-118. 王二七, 张旗, Burchfiel C B. 2000. 青海拉鸡山: 一个多阶段抬升的构造窗[J]. 地质科学, 35(4): 393-500. 王金荣, 张新虎, 张翔, 翟新伟, 张本旗, 赵磊. 2012. 中祁连西段野马南山构造岩石组合、成因及其构造意义[C]//全国岩石学与地球动力学研讨会论文摘要集, 466-468. 夏林圻, 夏祖春, 任有祥, 彭礼贵, 张诚, 杨静华, 王兴安, 李智佩, 韩松, 黄忠祥. 1991. 祁连、秦岭山系海相火山岩[M]. 武汉: 中国地质大学出版社, 1-304. 肖序常, 陈国铭, 朱志直. 1978. 祁连山古蛇绿岩带的地质构造意义[J]. 地质学报, 52(4): 281-295. 许志琴, 徐惠芬, 张建新, 李海兵, 朱志直, 曲景川, 陈代璋, 陈金禄, 杨开春. 1994. 北祁连走廊南山加里东俯冲杂岩增生地体及其动力学[J]. 地质学报, 68(1): 1-15. 闫臻, 王宗起, 李继亮, 许志琴, 邓晋福. 2012. 西秦岭楔的构造属性及其增生造山过程[J]. 岩石学报, 28(6): 1808-1828. 张建新, 许志琴, 陈文, 徐惠芬. 1995. 北祁连中段俯冲-增生杂岩/火山弧的时代探讨[J]. 岩石矿物杂志, 16(2): 112-119. 张雪亭, 杨生德. 2007. 青海省板块构造研究1/100万青海省大地构造图说明书[M]. 北京: 地质出版社, 1-221. 赵文津, Mechie J, 冯梅, 史大年, 薛光琪, 宿和平, 宋洋, 杨宏伟, 刘志伟. 2014. 祁连山造山作用与岩石圈地幔的特型结构构造[J]. 中国地质, 41(5): 1411-1423. 郑英, 陈光庭, 张小永, 宋泰忠, 梁坤先. 2017. 中祁连苏里地区奥陶纪石英闪长岩地球化学特征及年代学意义[J]. 西北地质, 50(4): 588-597. 钟明杰. 1994. 青海拉脊山下古生界的发现[J]. 地质论评, 22(5): 347-347. 左国朝, 张淑玲, 程建生, 巩彦学, 王彦斌, 吴汉泉. 1996. 祁连地区蛇绿岩带划分及其构造意义[C]//蛇绿岩与地球动力学研讨会论文集, 129-134. -
Access History
Figures(7)
Tables(2)
Export File
Citation
ZHANG Jinming, CHEN Guangting, CAI Hangjia, TIAN Chengxiu, LEI Xiaoqing. 2023. Geochemical characteristics and zircon U-Pb age of gabbros in the Zhaqiaohe ophiolite mélange, and its limitation on the ocean ridge environment[J]. Geology in China, 50(6): 1837-1847. doi: 10.12029/gc20200416002
Format
Content
- 表 1. 单平面滑动岩体稳定性评价Table 1. Stability evaluation of single-plane sliding rock mass
稳定性系数 稳定性分级 稳定 基本稳定 欠稳定 不稳定 | Show TableDownLoad:
CSV
- 表 2. 水电工程危险岩体规模分级Table 2. Scale classification of dangerous rock mass in hydropower projects
评价依据 小型 中型 大型 超大型 体积/m3 | Show TableDownLoad:
CSV
- 表 3. 异形滑移式块体特征数据统计Table 3. Statistical analysis of characteristic data for profiled blocks
块体编号 后壁倾角/（°） 后壁面积/m2 块体体积/m3 最大高差/m L1 50.35 728.15 1712.30 7.65 L2 45.58 835.53 2398.60 10.60 L3 40.42 842.94 2299.47 9.48 L4 26.23 706.89 3294.37 14.53 L5 46.61 886.78 1690.18 7.71 | Show TableDownLoad:
CSV
- 表 4. 异形滑移式块体物理力学参数Table 4. Physical and mechanical parameters of profiled block
参数 块体重度/（kN·m−3） 结构面黏聚力/MPa 结构面内摩擦角/（°） 取值 25.8 0.100 35 | Show TableDownLoad:
CSV
- 表 5. 危险岩体稳定系数的计算与定性Table 5. Calculation and qualitative characterization of stability coefficient of dangerous rock mass
块体编号 稳定性系数 稳定性分析 块体性质 L1 0.58 不稳定 是 L2 0.69 不稳定 是 L3 0.82 不稳定 是 L4 1.42 稳定 否 L5 0.66 不稳定 是 | Show TableDownLoad:
CSV
- 表 6. 危险岩体与岩体体积评价Table 6. dangerous rock mass and rock massvolume evaluation
危岩体编号 危岩体体积/m3 危岩体体积评价 L1 1712.30 大型 L2 2398.60 大型 L3 2299.47 大型 L5 1690.18 大型 | Show TableDownLoad:
CSV
稳定性系数 稳定性分级 稳定 基本稳定 欠稳定 不稳定 评价依据 小型 中型 大型 超大型 体积/m3 块体编号 后壁倾角/（°） 后壁面积/m2 块体体积/m3 最大高差/m L1 50.35 728.15 1712.30 7.65 L2 45.58 835.53 2398.60 10.60 L3 40.42 842.94 2299.47 9.48 L4 26.23 706.89 3294.37 14.53 L5 46.61 886.78 1690.18 7.71 参数 块体重度/（kN·m−3） 结构面黏聚力/MPa 结构面内摩擦角/（°） 取值 25.8 0.100 35 块体编号 稳定性系数 稳定性分析 块体性质 L1 0.58 不稳定 是 L2 0.69 不稳定 是 L3 0.82 不稳定 是 L4 1.42 稳定 否 L5 0.66 不稳定 是 危岩体编号 危岩体体积/m3 危岩体体积评价 L1 1712.30 大型 L2 2398.60 大型 L3 2299.47 大型 L5 1690.18 大型
-
Figure 1.
Tectonic location map (a) and regional tectonic map of Zhaqiaohe area (b)
-
Figure 2.
Field and microscopic features of the Zhaqiaohe gabbro
-
Figure 3.
TAS diagram (a, after Middlemost, 1994) and SiO2-FeOT/MgO diagram (b, after Miyashiro and Shido, 1975) of altered gabbro from Zhaqiaohe area
-
Figure 4.
Chondrite normalized REE distribution patterns (a) and primitive mantle-normalized trace lements patterns (b) of the samples from Zhaqiaohe area (normalization values after Sun and McDonough, 1989)
-
Figure 5.
CL images of the zircons from Zhaqiaohe area
-
Figure 6.
Zircon U-Pb age harmonic diagram (a) and weighted diagram (b) of gabbros from Zhaqiaohe area
-
Figure 7.
Zr-Zr/Y diagram (a, after Pearce and Norry, 1979) and Hf/3-Th-Ta diagram (b, after Wood, 1980) of gabbros from Zhaqiaohe area