基底侵蚀作用下泥石流动力过程分析

华 曙 光

(中国人民武装警察部队黄金第五支队，陕西 西安 710100)



基底侵蚀作用下泥石流动力过程分析

华 曙 光

(中国人民武装警察部队黄金第五支队，陕西 西安 710100)

基底侵蚀是运动物体与浅层地表在沿非滑动接触面剪切的时候而产生的相互作用，该作用可以改变泥石流系统的动力作用机制，从而影响泥石流的危害程度，是泥石流研究中的一项重要议题。本研究提出侵蚀率概念，分析了初始运动物质与夹带物质系统的侵蚀关系；同时建立了相关的数学模型和流变本构的数值实现，构建了数值模拟试验进行参数化评估。详细阐述了基底侵蚀作用对泥石流运动行为和整体流动性的影响。选取的流变模型和参数会产生比较好的模拟结果，侵蚀速率概念以及无量纲形式和动态模型的提出会给具有强烈侵蚀泥石流分析提供一套科学的潜在应用，也为地质灾害防治提供一定的技术参考。

基底侵蚀；泥石流；动力过程

基底侵蚀是泥石流运动过程中最容易观察到的现象，它可以使得表层物质不断卷入泥石流中，使得运动系统的体积不断增大，同样物质也可以在尾部被分离。在泥石流的运动过程中，物质的改变影响着流动系统动力机制及其危害程度。关于基底侵蚀机理的研究是一项非常有意义却又有争议的课题。CANNON和SAVAGE(1988)发展了质量集中模型来研究物质卷入和分离过程。该模型引入了推动物质运动的作用与反作用力，直接作用于泥石流，当流动物质的质量增加或减少，物体动量便会自动的增加或减少。与此类似，VAN Gassen & CRUDEN (1989)揭示了泥石流沉积物质的传递，同时计算了泥石流运动距离。然而，HUNGR (1990a, 1990b)的相关研究表明，不存在这样的力来推动泥石流的运动，从流动物质中被排出的物质速度为零。否则典型泥石流在运动过程中具有的巨大推力需要额外的能量输入，但这是在重力驱动情况下是不存在的。ERLICHSON(1991)也指出，泥石流在运动过程中存在沉积物质与基底具有相互作用而不是物质之间具有相互作用，从而得出，运动物质的守恒方程并不直接适用于泥石流运动系统。泥石流在实际运动过程中基底侵蚀和浅表层物质的卷入是经常发生的，对泥石流、滑坡运动的影响十分关键，是一项重要的科学议题。鉴于以上分析，笔者采用泥石流概括模型，提出了侵蚀率概念，构建三维连续介质动态数值模型及求解方法，运用数值试验对泥石流运动过程、动力行为、基底侵蚀作用进行分析，以期为地质灾害防治提供一定的技术参考。

1 模型建立与参数选取

1.1 模型选取

泥石流在运动过程中属于非牛顿流体，由于速度梯度而引起摩擦，产生剪切力；而流变特性正是反映流体变形所受剪应力与剪切率的关系(CHEN et al., 2006; CROSTA et al., 2004)。通常采用摩擦模型、Voellmy模型和Bingham模型描述流变特性(CHEN and LI, 2000；2003；CHEN et al., 2006)，笔者选用摩擦模型(Friction model)和Voellmy model进行数值模拟试验研究(HUNGR，1990a, 1990b；HUNGR et al., 1984)。

1.1.1 摩擦模型 (Friction model)

在细颗粒较少而粗颗粒较多的条件下，底面剪切阻力主要表现出摩擦的行为，摩擦阻力与地面有效正应力成正比，如式(1)所示：

(1)

σZ为底面有效正应力，τ为摩擦阻力，孔隙水压力比率γu和动力底摩擦角：δ，作为流变参数引入摩擦模型中，孔隙水压力比率γu等于孔隙水压力u与底面应力σZ比值。根据太沙基理论，在底面总应力一定的条件下，孔隙水压力的存在将会减小有效正应力，从而减小底面摩擦阻力。孔隙水压力比率与动力摩擦角可以由单一变量体摩擦角φb统一表达，φb定义为：

(2)

摩擦阻力对于泥石流运动而言是一种非常重要的阻力类型，对于底面剪切阻力主要由摩擦产生的泥石流而言，底面孔隙水压力的存在常常用来解释其高速和长距离的运动原因。

1.1.2 Voellmy模型 (Voellmy model)

Voellmy剪切阻力模型主要考虑了摩擦和湍流行为，已经被广泛应用于颗粒流和岩石崩塌、流动的模型之中，用下式表示：

(3)

1.2 流变模型参数选取

在流变模型确定的基础上，选用流变参数。其中摩擦流变模型只与动力底摩擦角δ一个未知参数有关，其值与岩土体的含水量有关；Voellmy流变模型包括动力摩擦系数μ和湍流系数ξ2个未知参数(表1)，属于典型稀性泥石流范畴。

表1 流变模型参数选取表

1.3 侵蚀率计算与选取

侵蚀是泥石流动力过程中常见的现象，也是一个非常重要而又极其复杂的课题。由于异质性土壤理化性质和复杂的斜坡地形，总体等效的侵蚀速率只能通过已知事件的反分析进行确定。笔者根据CHEN and LI，2000，2003；CHEN et al., 2006所提出的侵蚀速率这一概念，进行数值模拟研究。

侵蚀速率E定义为法向沟床下降/流失速率与平均切向滑动速率之间的比率。根据诸多的研究成果(CHEN and LI, 2000，2003；CHEN et al., 2006)表明，E的数量级通常为10-3。尽管沉积体积可以达到初始体积的多倍，由于侵蚀物质是从全部泥石流冲痕面积中收集而来，沟床地形的变化很小。

由于异质性土壤理化性质和复杂的斜坡地形，总体/等效的侵蚀速率只能通过已知事件的后预报分析确定。然而，对于地质工程实践中初步估计来说，可以使用如下的简化准则：

(4)

式(4)中，Veroded是总侵蚀体积，Aeffect代表总体受影响的侵蚀面积，dcenter代表物质中心的移动距离，所有变量均可以通过野外调查大致确定。考虑到系统非线性，修正系数α选取需要考虑流变模型和底部地形的影响。

对于地质工程实践初步评估来说，可以使用式(4)计算。为了使侵蚀率估算更为容易，笔者通过数值模拟试验计算Aeffect和dcenter，获得相关系数，同时不断修正来真实反映地形和流变模型的可靠性。根据前人的研究成果(CHEN and LI, 2000，2003；CHEN et al., 2006)，对于摩擦模型，不论限制坡还是无限制坡均采用α=2.0；对于Voellmy模型，此时非线性更强，对于无限制坡面，采用α=2.2，对于限制坡面采用α=2.5。

1.4 模型几何结构

为了评估泥石流侵蚀行为，不同坡型的泥石流动力特征，本试验构建了不同的基础地形(A系列：无限制形坡，B系列：渠道坡，该坡可以在侧向限制物质运动)，采用2种流变模型Friction模型、Voellmy模型，3种侵蚀率(E=0，0.003，0.007)，进行数值模拟计算。模型几何结构设计成一倾斜的35°坡面，一个水平面，之间采用5°线性过渡区连接，x-y平面内的直角坐标系位于与初始物质质量中心投影一致的水平区域。x轴方向与斜坡倾向一致(图1)。所有模拟试验均是在斜坡上部同一位置，以相同数量的物源开始。

图1 数值模拟试验结构示意图Fig.1 Set-up of the numerical experiments

整个斜坡在y轴方向宽度为8 m (y:-4～+4)，在x轴方向坡面长11 m(x：-2～+9 m)，过渡区长2 m(x：+9～+11 m)，水平区域长9 m(x：+11～+20 m)。限制坡下凹部分宽4 m(y:-2～+2)，最深处位于x=0处，深度h均为0.25 m。物源区几何体为一球缺，在x-y平面投影的为圆，其投影半径R=2 m，球缺高H=0.25 m。因此，物源区体积为V0为0.4m3。材料在坡面上突然释放，经过过渡区域运动到水平面上。

2 结果与讨论

2.1 摩擦模型模拟结果

泥石流整个运动过程在x-y平面内的运动轮廓投影和运动物质深度情况如图2、图3所示。图中虚线代表运动物质不同时刻轮廓线，时间间隔0.5 s，等值线云图代表运动物质堆积区深度分布情况。

图2(a)展示了摩擦模型下(Friction model)无限制坡面在不同的侵蚀速率情况下，运动物质轮廓边界在x-y平面的投影情况。可以看出，球缺运动物质一旦释放，它会延斜坡迅速加速下滑，在运动过程中，由于没有边界限制，物质会向两侧扩散，轮廓从圆形变为前端较钝尾部较尖泪珠形状。当物质前缘越过过渡区进入水平区域时，前端迅速减速，而尾部依然向前推进，当物质最终静止时，物质轮廓呈现马蹄形。非侵蚀结果见图2(a)。计算结果非常典型，同KOCH et al. (1994)采用内摩擦角39°，底摩擦角29°的石英颗粒的试验结果相似。

图2 摩擦模型的泥石流运动轮廓模拟结果图Fig.2 Projection of debris profiles on horizontal plane of friction model

相同数量的物质在有限制的渠道中加速下滑，其运动情况见图2(b)。物质初始形状是圆形，物质轮廓沿着渠道拉长，逐渐形成蝌蚪状，横向扩展在整个运动过程中很小。尽管侵蚀速率不同，沿着倾斜渠道的运动模式没有明显的改变。当进入水平区域，尽管没有边界限制，纵向延伸依然明显，并且侵蚀速率越大，横向扩展越大。值得指出的是，非侵蚀结果图2(b)表明，其沉积物轮廓呈子弹状。其结果与GRAY et al. (1999) 采用内摩擦角40°，底摩擦角30°石英碎片(平均粒径2～4 mm)的渠道颗粒流试验结果极其相似。

从图2堆积区厚度等值线图中可以看出，在各种工况下，堆积区深度分布基本均呈现出中后部较厚而前端较薄，且侵蚀作用非常影响堆积区轮廓，即侵蚀作用越强，尾部堆积区范围越大。在运动过程中，如果大量的物质被侵蚀，无限制坡(系列A)较限制坡(系列B)而言，堆积区前端更钝；显然底部地形同样影响着堆积区轮廓。对于相同系列的运动情况，其侵蚀速率越大，运动距离越大。

从图2还可以看出，在运动初始阶段，运动物质尾部相对较厚、前端较薄。随着物体沿坡下降，尽管有额外的物质不断卷入，其平均深度仍然逐渐降低。当物质在堆积区不断沉积，随着累计体积不断增加，物质深度逐渐增加，物质最终形成了一个主体相对较厚，前端相对较薄的形态。

在运动过程中，其动力参数的变化过程见图4。可以看出，对于限制坡和无限制坡流动，平均速度、最大速度的变化几乎不受侵蚀速率的影响。也就是说采用摩擦模型时，即使有大量的物质被卷入，也并未明显改变整体下滑的运动性。同样可以看出，无限制坡流动的整体运动性只是轻微高于限制坡流动时的情况。体积随时间的变化展示了物质积累过程，在相同的地形条件下，在更大的侵蚀速率下，坡面流会卷走更多的物质，同样也具有更大的动能。在相同的侵蚀速率下，无限制坡坡面流与限制坡坡面流相比，会卷走更多的物质，这点表明底部地形条件在侵蚀过程中起到关键性作用(图3)。

图3 摩擦模型的泥石流动力参数模拟结果图Fig.3 Time records of some kinematic parameters of friction model

2.2 Voellmy模型模拟结果

图4(a)展示了Voellmy模型下，无限制坡坡面流边缘随时间的变化情况。受Voellmy流变特性中湍流的影响，物体边缘从最初的圆形变化成为泪珠状。侧面扩展十分显著，且更大的侵蚀速率会导致更平的后部。对于有限制的渠道流而言，纵向延伸依然是主要的，由于侧向限制，横向扩展相对较小，结果如图4(b)所示。由于受地形的部分限制，沉积区前端停留在水平区域内，呈凸型，而后端依然滞留于过渡区内，呈窄凹形。

图4 Voellmy模型的泥石流动力参数模拟结果图Fig.4 Projection of debris profiles on horizontal plane of Voellmy model

沉积物剖面和厚度分布情况也可在图4中体现出来。可以看出，无限制坡坡面流动其堆积区前端呈钝形，限制坡渠道流其堆积区前端呈舌形。泥石流的传播距离和堆积区的宽度与侵蚀速率的大小有关。对于相同的侵蚀速率，限制坡渠道流的传播距离比无限制流坡面流更远。还可以看出，尽管堆积物前端和尾部较薄，最大堆积物厚度大部分位于堆积区中段，或者轻微朝向前端；堆积区物质的中心基本在最大深处位置附近。移动物体轮廓从最初启动时的尾部较厚、前端较薄变化为堆积区内主体较厚、近尾端较薄的形状。侵蚀速率越大，其沉积层厚度越大。

比较图5(a)和图5(b)可以看出，侵蚀作用和冲出的几何形态强烈影响平均速度的变化情况。总的来说，由于侵蚀作用而导致的系统惯性增加，泥石流整体运动的越快，停止的也就越快。在相同侵蚀速率下，限制坡渠道流平均速度与无限制坡坡面流相比，有一个更高的峰值，但下降速度也更快。体积随时间的变化情况如图5所示，可以看出，对于相同的侵蚀速率，无限制坡坡面流与渠道流相比，在运动过程中可以卷走更多的沟床物质。前端速度(最大速度)的变化过程如图5所示，在各个侵蚀速率和地形条件下，它的变化趋势基本一致。在相同的地形条件下，较大的前端速度往往与较高的侵蚀速率相对应；在相同的侵蚀速率情况下，运动相同的距离，无限制坡坡面流相比渠道流而言，前端速度更慢。

2.3 侵蚀作用对泥石流动力过程的影响分析

摩擦模型(Friction model)和Voellmy模型都广泛应用于描述颗粒流、崩塌以及高速远程滑坡的运动特性。值得注意的是，采用摩擦模型(图2、图3)，底部侵蚀作用的存在只少量影响了2种坡型平均流速的变化模式。6种情况下的物质中心基本都在同一位置，这表明在摩擦流变模型下，尽管有大量的侵蚀物质卷入到运动系统中，但运动物质的纵向变形与整体下滑的运动性并未发生本质的变化。这项观察与HUNGR & EVANS(1997) 所做的分析实例完全一致，他猜测是由于离心力抵消了与夹带物质有关的最初潜在能量。而笔者的推断与此不同，认为造成这种现象的原因与主要下降方向的凹形地形有关。笔者通过监测在整个运动过程中的所有作用力的变化，表明作用于斜坡的剪切力T和在运动方向的重力G的数量级大约是净体积力的2倍，并且物质的改变对运动系统的影响是直接通过净体积力P实现的。沿下坡方向的变形比横向扩展更加明显，因此，当净体积力P项在技术上忽略的时候，沿斜坡方向的整体运动性不受夹带作用的明显影响。在水平区域内，重力G对系统运动没有贡献，此时通过净体积力P调整来调整物质的局部深度，使其很慢地达到土压力平衡状态。此时尽管P和T基本是同一数量级，但在堆积阶段没有发生明显的夹带作用。然而，值得注意的是，以上观测结果并不适用于VOELLMY流变模型，这是由于在VOELLMY流变模型中，随着夹带物质的增加使堆积物深度逐渐增加，能够明显减小剪切力T。因此，整体运动性会逐渐增加，并且推动物质中心移动更远。

在数值模拟试验中，在同一流变模型和同一侵蚀速率情况下，渠道流相比无限制坡坡面流前端可以传播更远(图2、图4)，且同时具有更高的平均速度和前端速度以及更大的总动能；这点与渠道的影响可以加强运动距离这一现象相一致，并且非限制坡与渠道坡相比可以会卷走更多的物质，这是因为限制坡强烈影响了运动物质的整体位移和堆积区厚度分布，并且和上部分的分析结果一致。在实际应用中，对于泥石流和高速远程滑坡的灾害减缓措施而言，如堤坝和线形渠道一些技术都会部署用以限制和转移物体流动，使得物体以可控制的方式运动。而且，通过图3和图5展示的体积积累过程，表明由于横向扩展作用，无限制坡坡面流与渠道流相比，可以侵蚀掉更多的沟道与坡面物质。野外观测到无限制坡坡面流动典型地呈现出宽范围的夹带和沉积体积，相反限制流会产生比较平稳的夹带物体积和较小的堆积体积 (FANNIN & WISE，2001)，这点与计算结果相一致。计算结果也表明下垫面地形条件是一项影响泥石流、滑坡运动和沉积过程的重要因子。

3 结论

(1)侵蚀是泥石流动力过程中常见的现象，是运动物体与浅层地表在沿接触面的剪切时候而产生的相互作用，这种动态行为影响着进入运动系统的体积。泥石流的危害规模更多的是受延程进入物质的体积所决定，而不是其初始物质的体积决定。通过对泥石流过去事件的调查，发现沟道侵蚀开始通常都是有一个相对小的表面侵蚀面积，但在长距离的运动过程中，都会出现较深的切沟；相反，开放的无限制坡面大体上都伴随大面积的浅层表面侵蚀出现。因此，在研究泥石流动力过程时，要考虑都基底侵蚀作用。

(2)下垫面地形条件以及侵蚀率是影响泥石流运动和沉积过程的重要因子，渠道流相比无限制坡坡面流，前端可以传播更远；由于横向扩展作用，无限制坡坡面流与渠道流相比，可以侵蚀掉更多的沟道与坡面物质。侵蚀速率概念和动态模型的提出会给有强烈侵蚀泥石流的分析提供一套科学的潜在应用，也为地质灾害防治提供一定的技术参考。

CANNON, S. H. & Savage, W. Z.A mass change model for the estimation of debris flow runout[J]. J. Geol.1988, 96, 221-227.

CHEN H, LEE, CF. Numerical simulation of debris flows[J]. Canadian Geotechnical Journal, 2000,37:146-160.

CHEN H, LEE CF. A dynamic model for rainfall-induced landslides on natural slopes[J]. Geomorphology, 2003, 51: 269-288.

CHEN H, DADSON S, CHI YG. Recent rainfall-induced landslides and debris flow in northern Taiwan[J]. Geomorphology, 2006, 77:112-125.

CROSTA GB, CHEN H, LEE CF. Replay of the 1987 Val Pola Landslide, Italian Alps[J].Geomorphology,2004,60, Nos 1-2：127-146.

ERLICHSON, H. A mass-change model for the estimation of debris-flow runout, a second discussion: Conditions for the application of the rocket equation[J]. J. Geol.1991, 99：633-634.

FANNIN, R. J. & WISE, M. P. An empirical-statistical modelfor debris flow travel distance. Can. Geotech[J]. 2001,38(5)：982-994.

GRAY, J. M. N. T., WIELAND, M. &Hutter, K. Gravity-drivenfree surface flow of granular avalanches over complex basaltopography. Proc. R. Soc. London, Ser.1999, A455：1841-1874.

HUNGR, O. & EVANS, S. G. A dynamic model for landslideswith changing mass.Proceedings of the IAEG internationalsymposium on engineering geology and the environment (eds P.G. Marinos, G. C. Koukis, G. C. Tsiambaos and G. C. Stournaras),1997,1： 719-724. Rotterdam: Balkema.

HUNGR, O. A mass-change model for the estimation of debris-flow runout: a discussion. J. Geol.1990a,98：791.

HUNGR, O. Momentum transfer and friction in the debris of rock avalanche: discussion. Can. Geotech.J.1990b,27：697.

HINGR O, MORGAN GC, KELLERHALS R. Quantitative analysis of debris torrent hazards for design of remedial measures[J]. Can GeotechJ,1984,21:663-677.

KOCH, T., GREVE, R. & HUTTER, K. Unconfined flow ofgranular avalanches along a partly curved surface. II. Experiments and numerical computations. Proc. R. Soc. London, Ser.1994, A445, 415-435.

VAN GASSEN, W. & CRUDEN, D. M. Momentum transfer and friction in the debris of rock avalanche. Can. Geotech. J.1989,26：623-628.

中国地质调查局资料馆西北分馆2016年下半年新增地质调查成果资料目录

续表

序号顺序号资 料 名 称形 成 单 位形成时间321592陕西省宁强县代家坝勉县艾叶口铜金矿调查地质报告西北有色地质勘查局地质勘查院2015．03331593西北黄土高原区地质灾害综合研究成果报告中国地质调查局西安地质调查中心2015．12341594北山祁连关键地区区域地质调查中止项目工作总结中国地质调查局西安地质调查中心2015．12351595青海赤布张错、温泉兵站、直根尕卡1/25万区域重力调查成果报告陕西省地质调查院2015．11361596新疆万滴水泉幅L45C0030041/25区域地质调查报告陕西省地质调查院2015．11371597新疆万北塔山牧场幅L46C0030011/25区域地质调查报告陕西省地质调查院2015．11381598新疆和静县哈不欺汗拉矿区铜多金属矿调查评价报告河南有色金属地质矿产局第五地质大队2014．12391599祁连成矿带西段1/5万航磁调查成果报告中国国土资源航空物探遥感中心2015．11401600青海祁漫塔格成矿带金属矿地质矿产调查报告南京地质调查中心2015．11411601甘肃省合水南部宁县北部煤炭资源调查评价报告甘肃煤炭地质勘查院2013．3421602青海省德令哈市雅沙图地区金锑矿调查评价报告青海省核工业地质局2015．12431603甘肃省重要地质遗迹调查报告甘肃省地质环境监测院2015．10441604黄河上游巨型滑坡形成机理研究成果报告中国地质环境监测院 青海省环境地质勘查局2015．5451605新疆西昆仑四连幅J44E020023、喀拉萨依幅J44E020024、苦牙克幅J44E021023、赛勒果勒幅J44E0210241/5区调报告新疆维吾尔自治区有色地质勘查局七○四队2015．12461606东昆仑成矿带地质矿产调查评价成果集成项目成果报告中国地质调查局西安地质调查中心2015．11471607新疆西昆仑阿恰勒幅J44E022003、塔木阿勒地幅J44E022004、康西瓦达坂幅J44E023003、康西瓦幅J44E0230041/5万区域地质调查报告新疆维吾尔自治区地质矿产开发局第二区域地质调查大队2015．08481608青海门源县银灿浪力克1/5万成矿预测成果报告中国地质调查局西安地质调查中心2015．12491609西北地区重要地质遗迹调查(新疆)报告新疆维吾尔自治区地质环境监测院2015．12501610陕西省凤县国安寺白杨沟地区矿产地质调查报告西北有色地质勘查局七一七队2015．2511611青海省化隆县高店昂思多地区矿产远景调查报告中国冶金地质总局山东正元地质勘查院2015．12521612青海盐池湾幅、德令哈真幅、诺木洪幅1/25万基础地质调查修测报告(区域重力)陕西省地质调查院2015．11531613阿尔泰关键地区区域地质调查中止项目工作总结中国地质调查局西安地质调查中心2016．3541614青海沱沱河地区I46E001005、I46E001006、I46E002005、I46E0020061/5万区域地质调查报告成都理工大学地质调查院2014．4．2551615新疆西准噶尔前进牧场地区L44E009024、L44E010024、L45E009001、L45E0100011/5万区域地质调查报告新疆地矿局第七地质大队2013．12561616青海省泽库县萨卡拉卡地区萨卡拉卡I47E004021、盖盖拉麻I47E004022、禾日寺I47E005021、欧和杂滩47E005022、泽库县I47E0060221/5万区域地质矿产调查报告西北有色地质勘查局地质勘查院2016．1571617青海龙亚拉地区温泉幅I46E017008、日龙玛亚楼幅I46E017009、姜梗日幅I46E018008、龙亚玛村幅I46E0180091/5区域地质调查成果报告四川省地质矿产勘查开发局川西北地质队2015．12581618新疆西昆仑克孜勒陶乡幅J43E009017、塔木村幅J43E010017、其木干村幅J43E010018、希力希力煤矿幅J43E011017、阿尔帕勒克村幅J43E0110181/5万区域地质调查报告江苏省地质调查研究院2015．7

续表

序号顺序号资 料 名 称形 成 单 位形成时间591619新疆尉犁县兴地地区矿产地质调查报告中化地质矿山总局地质研究院2015．11601620新疆觉罗塔格拜城一带铁磷矿资源远景调查成果报告中化地质矿山总局地质研究院2016．5611621新疆西准噶尔玛依塔巴克地区L45E011006、L45E011007、L45E012006、L45E0120071/5万区域地质调查报告长安大学地质调查研究院2014．12621622新疆巴里坤伊吾地区1/5万航磁调查综合研究成果报告中国国土资源航空物探遥感中心2015．12631623甘肃省矿山环境监测成果报告有色金属矿产地质调查院2015．10641624宁夏自治区矿山环境监测项目成果报告青海省地质调查院2015．12651625青海省矿山环境监测项目成果报告青海省地质调查院2015．12661626青海沱沱河巴斯湖地区多金属矿调查评价报告青海省第五地质矿产勘查院2015．5671627新疆西昆仑奥依亚依拉克幅J45E019001、其木来克幅J45E020001、巴什坑科勒幅J45E020002、阿克苏库勒J45E021002、阿力亚力克科勒幅J45E0210011/5万区域地质调查报告新疆维吾尔自治区地质调查院2015．4681628青海省格尔木市开木齐地区铁多金属矿评价报告青海省地质调查院2015．12691629沱沱河艰险地区I45E011024、I46E011001、I45E012024、I46E0120011/5万中止区调项目工作总结吉林省区域地质矿产调查院2014．5701630沱沱河艰险地区I46E011002、I46E011003、I46E012002、I46E0120031/5万中止区调项目工作总结吉林省区域地质矿产调查院2014．5711631新疆自治区矿山环境监测成果报告中国煤炭地质总局航测遥感局2015．11721632陕西省矿山环境监测成果报告中国煤炭地质总局航测遥感局2015．10731633青海省柴达木山地区哈腊昆德幅J46E012023、艾日温·查干哈达幅J46E012024、柯柯·额热格幅J46E013023、柯嫩·多勒格幅J46E0130241/5万矿产远景调查报告西北有色地质勘查局地质勘查院2015．12741634河西走廊石羊河流域国土资源综合监测成果集成与评价报告甘肃省环境监测院2016．2751635青藏高原中北部古特提斯岩浆演化与成矿地质背景中国地质大学(武汉)2015．10761636新疆东准噶尔地区隐伏矿地球化学勘查与深部矿体预测示范报告中国地质科学院地球物理地球化学勘查研究所2015．11771637新疆若羌县吐特喀维宝西北一带铁铜铅锌多金属矿产评价成果报告新疆地矿局物化探大队2015．12781638天水宝鸡地区1/25万基础地质调查修测工作报告中国地质调查局西安地质调查中心2015．12791639渭河盆地基底探测与选区地球物理攻关成果报告陕西省地质调查院2015．12801640东天山成矿带浅覆盖区1/25万基础地质调查修测成果报告中国地质科学院地球物理地球化学勘查研究所2016．5．16811641陕西省牟家坝私渡河地区矿产地质调查报告陕西地矿局汉中地质大队2015．12821642沱沱河艰险地区太阳湖幅I46E001003、小太阳湖幅I46E001004、马兰山幅I46E002003、乱沟山幅I46E0020041/5万中止区调项目工作总结成都理工大学2014．12831643青海省都兰县黑石山地区铜多金属矿调查评价报告青海省第一地质矿产勘查院2015．12841644新疆西昆仑天水湖I44E002002、土兹达万山I44E002003、土兹达万山口东I44E002004、大红柳滩南I44E0020051/5万区域地质调查报告新疆维吾尔自治区地质调查院2015．9．25

续表

序号顺序号资 料 名 称形 成 单 位形成时间851645北山地区斑岩型铜矿特征矿物测试技术研究成果报告中国地质调查局西安地质调查中心2015．12861646青海东台吉乃尔湖幅、大灶火幅(部分)、格尔木幅(部分)1/25万区域重力调查报告中国地质调查局西安地质调查中心2015．12871647青海矿产资源开发地质环境影响调查评价报告青海省水文地质工程地质环境地质调查院青海省环境地质调查局2015．12881648柴达木盆地循环经济实验区重点地区水文地质调查报告青海省水文地质工程环境地质调查院2016．3891649塔里木中北部地区航空重力调查成果报告中国地质调查局油气资源调查中心2016．6901650秦祁昆造山系地质构造演化与成矿规律研究成果报告中国地质调查局西安地质调查中心2016．4911651青海省格尔木市野马泉地区1/5万综合物探调查成果报告青海省第三地质矿产勘查院2015．12921652新疆重点地区航空高光谱调查与找矿预测技术研究成果报告有色金属矿产地质调查中心2015．11931653青海沱沱河地区扎木曲扎拉夏各涌铜多金属矿调查评价报告青海省第五地质矿产勘查院2016．1941654新疆西昆仑库拉木阿特达坂幅J43E016018、阿子岗沙勒幅J43E016019、莫莫克幅J43E016020、半德尔幅J43E017018、玉萨斯幅J43E017019、柯克亚煤矿幅J43E0170201/5万区域地质调查报告新疆维吾尔自治区地质调查院2015．6951655新疆呼图壁⁃沙湾地区1/25万多目标区域地球化学调查报告新疆维吾尔自治区地质调查院2015．5961656新疆青河县吉尔斯幅L46E011001、二台幅L46E012001、结勒得喀腊幅L46E013001、沙尔托海幅L46E0130021/5万区域地质调查报告陕西省地质矿产勘查开发局区域地质矿产研究院2015．12971657新疆西昆仑团结峰幅I44E005003、岔路口幅I44E005004、河尾滩幅I44E006003、兴山幅I44E006004、南大沟幅I44E0070031/5万区域地质调查报告四川省地质矿产勘查开发局区域地质调查队2016．2981658青海祁漫塔格成矿带金属矿地质矿产调查成果报告青海省第三地质矿产勘查院2015．12991659青海民和县幅、威武市幅1/25万基础地质调查修测结题报告(区域重力)青海省第三地质矿产勘查院2015．121001660青海沱沱河幅I46C002002、可可西里湖幅I46C001001、乌兰乌拉湖幅I46C002001、布喀达坂峰幅I46C0040011/25万区域重力调查成果报告青海省第三地质矿产勘查院2015．121011661青海丁字口幅、鱼卡幅、大柴旦镇幅1/25万区域重力调查成果报告青海省第三地质矿产勘查院2015．121021662青海省格尔木市卡而却卡地区1/5万区域地质综合调查成果报告青海省第三地质矿产勘查院2016．11031663青海省都兰县白兔滩J47E017001、巴日尔额格J47E018001、乌兰·哈布其勒J47E018002、额明尼克·乌拉南J47E0190031/5万区域地质矿产调查报告青海省第三地质矿产勘查院2015．121041664甘肃合作镇幅、临夏市幅1/25万基础地质调查修测成果报告(区域重力)中国地质科学院地球物理地球化学勘查研究所2015．121051665北山天山蛇绿岩成因类型、构造环境和SHRIMP年代学研究中国地质科学院地质研究所2014．121061666青藏高原中部古特提斯构造格架、演化过程与成矿报告中国地质科学院地质研究所、中国地质大学(北京)2015．61071667青海柴北缘J46E010018、J46E011018、J46E0110191/5万化探地球化学图说明书青海省第五地质矿产勘查院2016．1

续表

序号顺序号资 料 名 称形 成 单 位形成时间1081668渭河盆地氦资源调查重磁测量与化探采样成果报告陕西省地质矿产勘查开发局第二综合物探大队2016．11091669新疆西昆仑三十里营房幅J44E022001、布孜阿日吐幅J44E022002、国庆桥幅J44E023001、哈巴克幅J44E0230021/5万区域地质调查报告四川省核工业地质调查院2015．121101670天山戈壁沙漠覆盖区成矿地质背景研究与成矿要素综合推断成果报告中国地质大学(武汉)2015．121111671新疆西天山斑岩铜矿找矿突破研究成果报告中国地质大学(北京)2016．31121672新疆乌恰县萨热克铜矿及外围铜矿调查评价报告有色金属矿产地质调查中心2016．4．11131673新疆阿克陶县克拉克尔地区铜矿调查评价报告有色金属矿产地质调查中心2015．71141674东天山成矿带浅覆盖区1/25万基础地质调查修测成果报告北京探矿工程研究所2015．111151675陕西省商洛市竹林关新庙一带脉石英(水晶)资源远景调查报告中国建筑材料工业地质勘查中心陕西总队2015．71161676陕北奥陶纪盐盆地钾盐资源调查评价工作项目结题报告中国地质科学院矿产资源研究所2015．121171677陕西紫阳岚皋地区铁铜矿成矿环境与勘查研究成果报告中国地质科学院矿产资源研究所2015．121181678青海省都兰县埃坑地区铜钼矿调查评价报告青海省第五地质矿产勘查院2016．11191679中央造山带西段(阿尔金祁连西秦岭)早古生代洋陆格局及造山机制中国地质科学院地质研究所2015．101201680青藏高原北缘主要大型断裂带的组成、结构与新生代活动性研究成果报告中国地质科学院地质研究所2015．101211681新疆自治区矿山环境监测成果报告新疆维吾尔自治区地质调查院2015．111221682青海省共和县加当根地区铜多金属矿调查评价报告青海省第一地质矿产勘查院2016．11231683渭河盆地地热井氦及伴生气调查与综合利用评价成果报告长安大学地质调查院、陕西燃气集团渭河能源开发有限责任公司2015．121241684西北地区国土资源遥感综合调查与监测成果报告新疆维吾尔自治区地质调查院新疆地质矿产勘查开发局第一区域地质调查大队2013．21251685青藏高原重点地区水文地质环境地质调查综合研究成果报告中国地质调查局水文地质环境地质调查中心2015．121261686宁夏中南部严重缺水地区地下水勘查与供水安全示范成果报告中国地质调查局水文地质环境地质调查中心2015．121271687关中城市群活动断裂地球物理调查成果报告中国地质科学院地球物理地球化学勘查研究所2016．11281688鄂尔多斯大型盆地地下水形成演化机理研究成果报告中国地质大学(北京)2016．11291689青海木乃地区I46E016010、I46E016011、I46E016012、I46E017010、I46E017011、I46E0170121/5万化探地球化学图说明书青海省第五地质矿产勘查院2014．121301690青海不冻泉卡巴纽尔多地区1/25万区域化探地球化学图说明书青海省第五地质矿产勘查院2015．111311691青海省德令哈市措纳合作哈幅J47E009006、瑙滚诺尔幅J47E010005、苏令郭勒幅J47E010006、敖果吞乌兰幅J47E0110051/5万区域地质调查报告青海省地质调查院2016．21321692甘肃省西和县崖湾一带金锑矿远景调查成果报告甘肃省地质调查院2015．121331693青海都兰县幅、茶卡镇幅、西宁市幅1/25基础地质调查修测成果报告(区域重力)青海省第三地质矿产勘查院2015．111341694青海茫崖采石沟地区阿卡腾能山幅J46E010002、采石场幅J46E010003、红沟子幅J46E009004、红沟子西幅J46E0090031/5万区域地质调查报告青海省有色地质矿产勘查院2012．9

续表

序号顺序号资 料 名 称形 成 单 位形成时间1351695青海当江地区索莫不久幅I46E015022、当江幅I46E015023、索莫不久下幅I46E016022、当江下幅I46E0160231/5万区域地质调查报告青海省有色地质矿产勘查院2012．21361696宁夏主要城市浅层地温能开发区1/5万水文地质调查报告宁夏回族自治区国土资源调查监测院2015．121371697新疆铁厂沟镇幅L45C002001、克拉玛依市幅L45C0030011/25万区域地质调查报告中国地质大学(武汉)2015．111381698西准噶尔克拉玛依后山地区三维地质调查成果报告中国地质大学(武汉)2015．91391699新疆西昆仑琼萨依幅J44E021022、红山顶幅J44E022022、牙赤克幅J44E022023、宿营地幅J44E0230221/5万区域地质调查报告中国地质大学(武汉)2015．121401700甘肃肃北县扫子山地区矿产地质调查成果报告甘肃省地质调查院2015．61411701新疆富蕴县巴勒克卓尼幅L45E011021、却木希拜幅L45E012021906、高地幅L45E013021、塔克尔巴斯陶幅L45E0130221/5万区域地质调查报告新疆维吾尔自治区地质调查院2014．61421702青海省矿产资源勘查与开发利用信息化建设项目成果报告青海省国土资源科技信息中心2015．121431703西北地区典型矿山地质环境治理关键技术示范成果报告中国地质调查局西安地质调查中心2015．111441704青海省冷湖镇J46E008015、J46E009016、J46E009017、J46E0090181/5万矿产远景调查报告青海省地质调查院2016．21451705青海阿拉克湖冬给措纳湖地区阿拉克湖幅I47E001001、冬给措纳湖幅I47E0010021/25万区域化探地球化学图说明书青海省第五地质矿产勘查院2015．31461706塔里木盆地孔雀河流域和什力克幅K45E014008、亚特玛斯幅K45E0150081/5万水文地质调查综合评价报告新疆地矿局第二水文工程地质大队2015．121471707新疆喀什经济开发区水文地质环境地质调查评价报告新疆地矿局第二水文工程地质大队2015．121481708青海省格尔木市忠阳山地区铜矿调查评价总结报告(20132014年)青海省第一地质矿产勘查院2015．121491709青海省天峻县苏里地区矿产远景调查报告青海省地质调查院2016．21501710柴达木东部盐湖构造中浅钻钻探工程成果报告中国石油青海油田勘探事业部2016．61511711宁夏中南部严重缺水地区地下水勘查与供水安全示范成果报告(2010年)宁夏回族自治区地质调查院2015．121521712宁夏自治区矿山环境监测成果报告宁夏回族自治区地质调查院2015．121531713塔里木盆地库尔勒市幅K45E014009、哈拉玉宫乡幅K45E0150091/5万水文地质调查综合评价报告新疆地矿局第一水文工程地质大队2015．121541714河西走廊石羊河流域国土资源综合监测示范报告甘肃省地地质环境监测院2016．21551715陕西省榆林市米脂县地质灾害详细调查报告中国人民武装警察部队黄金第五支队2016．11561716青海拉脊山J47E023023、J47E024023、J47E024024、J48E024001、J48E0240021/5万化探成果报告山东省物化探勘查院2015．121571717黄河上游共和段地质灾害发育机理及防治对策研究报告中国地质环境监测院2016．11581718青海省大柴旦行委大滩西依克湖南煤炭资源调查报告中国煤炭地质总局青海煤炭地质局2014．51591719新疆西昆仑鱼跃石幅I43E002024、英明羊滩幅I44E002001、喀喇昆仑山口幅I43E003024、天神达坂幅I44E0030011/5万区域地质调查报告四川省地质矿产勘查开发局四〇五地质队2016．81601720覆盖区矿产预测综合研究与成果汇总成果报告中国地质大学(武汉)2015．12

续表

序号顺序号资 料 名 称形 成 单 位形成时间1611721青海省杂多县宰里涌曲地区铅锌矿产资源调查评价工作总结青海崯田地球物理化学勘查股份合作公司2015．41621722青海东部严重缺水地区地下水勘查(平安、民和)报告青海省水文地质工程地质环境地质调查院2015．81631723青海东部严重缺水地区地下水勘查报告青海省水文地质工程地质环境地质调查院2015．81641724青海滩间山锡铁山铜金多金属成矿带找矿类型及靶区优选中国地质大学(武汉)2015．111651725新疆西昆仑八一八水库幅J44E021019、塔勒坎勒克幅J44E021020、乌鲁克萨依幅J44E022019、阿克苏塔格幅J44E022020、喀什塔什塔格幅J44E0230191/5万区域地质调查报告陕西省地质调查院2016．11661726新疆西昆仑I44E002006、I44E002007、I44E002008、I44E0020091/5万区域地质调查报告陕西省地质调查院2016．11671727新疆富蕴县塔尔结拜地区矿产地质调查报告新疆地矿局第四地质大队2016．11681728青海门源回族自治县幅、山丹县幅1/5万基础地质调查修测报告(区域重力)中国地质调查局西安地质调查中心2015．121691729陕西省镇安地质灾害调查成果报告中国地质环境监测院2016．11701730高原寒区典型地区地质环境调查评价成果报告中国地质调查局水文地质环境地质调查中心2015．121711731高原寒区地质环境调监测关键技术研究成果报告中国地质调查局水文地质环境地质调查中心2015．121721732青海天俊县尕何幅J47E007006、嘎巴玛尔当幅J47E007007、客尔浪幅J47E008007、玛尔洛休玛幅J47E008008、苏里幅J47E0080091/5万区域地质调查报告青海省地质调查院2016．21731733新疆巴里坤县纸房地区矿产地质调查成果报告中国冶金地质总局西北地质勘查院2016．11741734青海省泽库县瓦勒根地区金矿普查报告青海省第一地质矿产勘查院2015．121751735甘肃省照壁山马鬃山地区矿产地质调查成果报告甘肃省地质调查院2015．121761736青海雅西措地区铜多金属矿调查评价报告青海省柴达木综合地质矿产勘查院2015．121771737青海省格尔木市黑海北铜金矿调查评价报告青海省柴达木综合地质矿产勘查院2015．111781738柴达木南缘盐沼带钾盐矿产调查评价报告青海省柴达木综合地质矿产勘查院2015．121791739关中城市群城市建设地质环境适宜性评价成果报告长安大学地质调查院2015．121801740北山祁连成矿带矿产资源潜力调查总结报告中国地质调查局西安地质调查中心2016．11811741宁夏香山地区1/25万区域化探地球化学图说明宁夏地球物理地球化学勘查院2015．121821742柴达木盆地循环经济实验区重点地区地下水勘查报告青海省环境地质勘查局2015．121831743柴达木盆地都兰县幅、夏日哈幅1/5万水文地质调查报告青海省环境地质勘查局2016．51841744陕西省黄陇煤田煤中伴生金属资源调查评价成果报告中国煤炭地质总局特种技术勘探中心2015．121851745青海省平安乐都地区土地质量地球化学调查评价报告青海省第五地质矿产勘查院2016．11861746柴达木盆地铁奎滩幅J47E0220081/5万水文地质调查报告青海省柴达木综合地质矿产勘查院2016．71871747甘肃阳山金矿带成矿系统与找矿方向研究成果报告中国地质大学(北京) 中国人民武装警察部队黄金指挥部2016．11881748新疆北屯市拟建城区水文地质环境地质调查评价报告新疆生产建设兵团地质勘查中心 新疆兵团勘测设计院(集团)有限责任公司工程勘察院2015．111891749青海化隆循化地区1/25万土地环境地质调查报告青海省第五地质矿产勘查院2015．111901750青海诺木洪都兰地区1/25万基础地质调查修测(化探)2014年度工作总结青海省第五地质矿产勘查院2015．121911751青海省主要城市浅层地温能开发区1/5万水文地质调查报告青海省水文地质工程地质环境地质调查院2015．12

续表

序号顺序号资 料 名 称形 成 单 位形成时间1921752新疆富蕴县卡拉麦里地区矿产地质调查成果报告新疆维吾尔自治区有色地质勘查局2016．11931753新疆特克斯县伊什基里克山一带铁多金属矿远景调查成果报告新疆维吾尔自治区有色地质勘查局2014．121941754新疆西昆仑J44E021018、J44E022017、J44E022018、J44E023017、J44E0230181/5万区域地质调查报告新疆地质矿产勘查开发局第一区域地质调查大队2016．11951755新疆西准噶尔二四八公里道班L45E011008、阿腊德L45E011009、沥青矿L45E012008、夏四井L45E0120091/5万区域地质调查报告新疆地质矿产勘查开发局第一区域地质调查大队2015．121961756甘肃陇南白水江流域文县幅I48E0190111/5万地质灾害调查报告甘肃省地质环境监测院2015．121971757青藏高原冻土带天然气水合物调查评价成果报告中国地质科学院矿产资源研究所2016．31981758新疆乌恰县克河地区矿产地质调查成果报告有色金属矿产地质调查中心2015．111991759新疆阿什图市喀达塔木地区矿产地质调查工作总结报告有色金属矿产地质调查中心2015．112001760青海祁漫塔格成矿带金属矿地质矿产调查成果报告中国地质科学院地球物理地球化学勘查研究所2016．62011761青海省中南祁连地区1/5万航空磁测成果报告中国冶金地质总局地球物理勘查院2016．62021762青海省冷湖镇小赛什腾地区金铜多金属矿调查评价报告青海省第四地质矿产勘查院2015．72031763青海省祁连县走廊南山南坡一带铜多金属矿调查评价报告(2012—2015)青海省第四地质矿产勘查院2015．122041764青海省德令哈市地质灾害调查与区划报告青海省环境地质勘查局2015．102051765青海省玛多县地质灾害调查与区划报告青海省环境地质勘查局2007．112061766青海省久治县地质灾害调查与区划报告青海省环境地质勘查局2007．102071767青海省班玛县地质灾害调查与区划报告青海省环境地质勘查局2008．122081768西南三江成矿带北段地质矿产综合调查总结报告中国地质调查局西安地质调查中心2016．42091769新疆东天山雅勒伯克地区矿产地质调查报告新疆地质调查院2015．122101770新疆和田地区康西瓦大红柳滩地区卡皮达兰沟幅J44E024003、红旗桥幅J44E024004、四六二道班幅J44E024005、错鲁勒错幅I44E001003、大红柳滩西幅I44E001004、大红柳滩幅I44E0010051/5万矿产地质调查报告新疆维吾尔自治区地质调查院2015．112111771青海省同仁县曲库乎乡江龙牧业铅多金属矿调查评价报告青海省有色地质矿产勘查局七队2016．72121772青海省祁连县小沙龙南(C15、C26等异常区)铁矿预查报告青海省有色地质矿产勘查局七队2015．122131773西北大型盆地关键水文地质参数试验研究成果报告长安大学2016．62141774青海省赛多浦岗日地区矿产远景调查报告青海省第五地质矿产勘查院2015．112151775青海省曲麻莱县旁安地区金矿调查评价报告青海省第五地质矿产勘查院2016．12161776青海省曲麻莱县诺尔扎尕玛地区东温泉幅I46E002023、德德哈日乌拉幅I46E003023、诺尔扎尕玛幅I46E0040231/5万区域地质矿产调查报告青海省第五地质矿产勘查院2016．12171777青海省天峻县纳日综地区矿产远景调查报告中国冶金地质总局山东正元地质勘查院2015．122181778青海省刚察县牵陇幅J47E015015、观户贡玛幅J47E016015、科才贡巴幅J47E016016、陶勒曲果幅J47E017015、泉吉幅J47E0170161/5万区域地质矿产调查报告四川省地质矿产勘查开发局二零七地质队、四川省冶金地质勘查局六〇五大队2015．122191779新疆西天山阿吾拉勒铁矿整装勘查区成果集成及火山建造构造综合调查研究成果报告新疆维吾尔自治区地质调查院2014．122201780新疆塔什库尔干县马尔洋克其克一带铁矿调查评价报告新疆维吾尔自治区地质调查院2015．122211781新疆巴楚县瓦吉里塔格钒钛磁铁矿调查评价报告新疆维吾尔自治区地质调查院2014．2

续表

Numerical Study about the Influence of the Basal Erosional and Block Effectson Dynamic Process of Debris Dlow

HUA Shuguang

(No.5 Gold Geological Party of Chinese Armed Police Force, Xi’an 710100, Shaanxi, China)

The mechanism of mass changes and block in debris transportation process is an important topic in the study of fast landslides, debris flows and avalanches. Basal erosion is recognized as a dynamic interaction between the original moving material and the entrained basal topsoil sheared along their non-slip contact surface.The dam is used to control mobility behavior of debris flows, to decrease the runout speed of debris flows, to reduce erosion capability and discharge, and tomake the flow direction diverge.In this paper, a new concept of yield rate has been proposed, and the erosional relationship between original moving material and the entrained basal topsoil has been analyzed. Apertinent mathematical model and numerical implementation have been formulated. Parametric numerical experiments have been conducted to compare the erosional effects. The influence of involved erosive material and block effects on the runout behavior and global mobility of moving material has been elucidated. The simulated results are consistent with available experimental and field observations. The excellent match to the field data suggests that it is a plausible potential application to analyze the gravity-driven flow with significant erosion. The defined dimensionless form and the proposed yield rate can be estimated conveniently in general geotechnical practices.

basal erosion; block dam;debris flow; rheological model; numerical modelling

序号顺序号资 料 名 称形 成 单 位形成时间2221782陕西省主要城市浅层地温能开发区1/5万水文地质调查报告陕西省地质调查院2015．122231783新疆西南天山皮羌地区矿产地质调查成果报告陕西省地质调查院2015．122241784陕西镇巴岚皋地区矿产地质调查成果报告陕西省地质调查院2015．122251785区域化探方法技术研究与成果集成(西北地区)成果报告中国地质调查局西安地质调查中心2015．112261786甘肃莫坝尔幅J46E005014、阿尔金山幅J46E005015、阿克塞沟幅J46E005016、黑山头幅J46E006014、托拉依格达坂幅J46E006015、大苏干湖北幅J46E0060161/5万区域地质调查报告甘肃省地质调查院2013．52271787青海万托来牧场幅、阳康幅1/25基础地质调查修测成果报告(区域重力)北京勘查技术工程有限公司2015．122281788青海沱沱河地区巴音莽鄂阿幅I46E006009、诺日加玉幅I46E006010、巴音查乌马幅I46E007009、查香结德幅I46E0070101/5万中止区调项目工作总结中国地质大学(武汉)2015．122291789新疆焉耆县哈都虎拉山地区矿产地质调查报告中国冶金地质总局第一地质勘查院2015．122301790甘肃白山泉、大一明镇、明水、白头山幅1/5万化探成果报告甘肃省地质调查院2015．52311791新疆苏立杰幅K44E009022、查汗沙拉幅K44E009023、布鲁坦哈尔沙拉幅K44E009024、雄呼尔音代勒幅K45E0090011/5万区域地质调查报告陕西省地质调查院2016．42321792青海省德令哈市查斯台·郭勒幅J47E012007、阿尔格·陶勒盖幅J47E012008、三岔东幅J47E013007、查干莫孙幅J47E013008、拜兴图幅J47E014007、乌兰·陶勒盖幅J47E0140081/5万区域地质矿产调查报告甘肃省地质矿产勘查开发局第一地质矿产勘查院2015．112331793青海省茫崖镇大通沟南山地区J46E007007、J46E007008、J46E008006、J46E008007、J46E0080081/5万矿产远景调查报告西北有色地质勘查局地质勘查院2015．122341794青海省唐古拉山地区I46E019018、I46E019019、I46E020018、I46E0200191/5万区域地质调查项目中止阶段性总结报告江苏华东地质调查集团有限公司 华东有色地质矿产勘查开发院2016．92351795东昆仑早古生代构造热事件及其造山作用成果报告中国地质科学院地质研究所2015．62361796秦祁昆结合部古生代中生代构造体制转换与成矿作用成果报告中国地质科学院地质研究所2015．52371797新疆西昆仑乔尔天山岔路口地区矿产地质调查报告新疆地矿局第十一地质大队2015．32381798青海省格尔木市雪水沟地区铅锌矿调查评价报告青海省地质调查院2016．42391799青海省红石沟地区铜多金属矿调查评价报告青海省地质调查院2016．42401800甘肃省红泉幅K46E014023、鲤鱼梁幅K46E014024、红岩井幅K47E014001、勒巴泉幅K47E014002、平头山幅K47E0150021/5万区域地质调查报告甘肃省地质调查院2015．122411801新疆西昆仑库孜幅J43E014017、瓦尔果勒幅J43E015017、布仑木莎幅J43E016017、皮勒幅J43E0170171/5万区域地质调查报告山东省第四地质矿产勘查院2015．122421802阿尔泰准噶尔北缘成矿带矿产资源调查总结报告中国地质调查局西安地质调查中心2015．122431803阿尔泰准噶尔北缘成矿带地质矿产调查综合研究报告中国地质调查局西安地质调查中心2015．122441804青海省天峻县快尔玛幅J47E016012、哈德木幅J47E016013、江河幅J47E016014、关角幅J47E017012、天峻幅J47E017013、天棚幅J47E0170141/5万地球化学普查成果报告河北省区域地质矿产调查研究所2016．22451805青海省兴海县长虫山地区查库贡巴幅I47E006013、下大武幅I47E006014、怡仓幅I47E006015、尕科合幅I47E0060161/5万区域地质调查报告四川省冶金地质勘查院2016．22461806青海东昆仑清水河—树维门科地区铜多金属矿调查评价报告青海省地质调查院2016．4

2016-05-29；

2016-08-08

国家自然科学基金项目(41302224)及陕西省科学技术研究发展计划项目(2014KJXX-20)共同资助

华曙光(1970-), 男, 安徽桐城人，博士，高级工程师，从事矿山地质环境保护与恢复治理、地质灾害等方面的研究工作。E-mail:huashuguang88@sina.com

TD824.7

A

1009-6248(2016)04-0219-09