STUDY ON ORE-CONTROLLING FACTORS AND MINE EXPLORATION PROPOSAL OF THE GOLD ORE FIELD IN NORTH HANYIN, SOUTH QINLING
-
摘要: 选取南秦岭汉阴北部金矿田内黄龙、长沟和金斗坡典型矿区为研究区,通过野外大比例尺构造-岩相填图、剖面实测、综合编图等工作,结合前人研究成果,查明了金矿田的控矿要素,并指出了有利的找矿勘探方向。综合研究表明,汉阴北部金矿田的矿化严格受控于志留系梅子垭组的不同岩性段;脆-韧性剪切带控制了矿床的空间分布,次级断层控制矿体产出的具体部位;S2面理对金矿化有显著的制约作用;促使石榴子石和黑云母变斑晶发育的热变质作用能够引起原始岩层中金元素发生活化,且热变质作用的晚期为退变质期,退变质期的适宜温度有利于金元素富集成矿。矿田内的金矿勘探应重点选在脆-韧性剪切带内,尤其针对不同岩性的接触面、次级断层面附近、石榴子石与黑云母变斑晶密集发育区周缘等区域,矿体追索方向应依照S2面理的空间优势展布方向。Abstract: Taking Huanglong, Changgou and Jindoupo gold ore areas located in Northern Hanyin ore field, Southern Qinling as the typical study area, through the large-scale tectonic lithofacies mapping, measuring geologic section and complex map compilation, the ore-controlling factors of this ore-field has eventually been acquired which would provide the directions of gold exploration based on previous research results. Comprehensive analysis shows that gold mineralization of this field is bound up with Silurian Meiziya Group stratum, and the brittle-plastic shear zone restricts the spatial distribution of gold deposits, and the locations of the ore bodies are decided by the secondary faults. S2 structural foliation obviously controls the gold mineralization, and the thermometamorphism which develops biotite and garnet porphyroblasts can cause the gold element to reactivate. And in the late stage it is the retrograde metamorphism period during which the temperature is very available for the form of gold deposits. Therefore the gold exploration in this area should be focus on the brittle-plastic shear zones in the future, especially the positions of those connection of different rock layer, the region near the secondary faults, and the scopes where biotite and garnet porphyroblasts massively appear. Meanwhile the direction of gold seeking should follow the main spatial distribution direction of S2 foliation.
-
0. 引言
边坡稳定性一直都是岩土地质研究中经久不衰的课题,地震及降雨引发边坡坍塌、滑坡以及泥石流等边坡失稳问题尤其如此。目前,单工况下边坡稳定性的研究层出不穷,多工况耦合作用下边坡稳定性的研究则报道较少。蒋中明等[1]通过自编FISH函数实现了边坡三维饱和渗流计算分析;刘春玲等[2]利用FLAC3D对某边坡地震稳定性进行了研究;齐信等[3]基于GIS技术对汶川地震强震区地震诱发滑坡与后期降雨诱发滑坡控制因子进行了耦合分析;龚文俊等[4]采用GEOSTUDIO软件以甘肃西和Ⅲ号滑坡为例研究了降雨和地震耦合作用对滑坡稳定性的影响。边坡稳定性与人类生产生活息息相关,因此,对边坡稳定性的研究显得尤为重要。
本文以成兰铁路岷江双线特大桥高边坡为例,结合多种数值分析方法的特点,开展天然工况、地震工况、降雨工况以及降雨-地震耦合工况下边坡稳定性研究和工程安全性评价,为单工况及耦合工况下边坡稳定性研究提供参考数据及技术依据。
1. 研究区概况
成兰铁路岷江双线特大桥高边坡位于四川省九寨沟县大录乡八郎沟玉瓦寨,边坡所在的2号岩堆大致呈扇形展布于斜坡上,岩堆体纵坡较陡,坡度约在30°-38°之间,坡面长满杂树,树林茂密。该区域处于我国地貌阶梯起伏最剧烈的深切陡峻高山峡谷。铁路线路沿岷江两岸而上,在该区域穿越了4条活动断裂带,共涉及规模不等的十几条活动断裂,平均滑动速率均大于0.5 mm/a以上[5]。边坡区域地层岩性纷杂,不良地质分布广,在天然条件下会引起坡内岩土侧移导致轻微滑坡[6]。岷江地震带属高烈度地震强震频发区,地震基本烈度为Ⅷ度,地震峰值加速度为0.2 g[7]。该区域内气候多变,降水量较多,易形成暴雨。
2. 模型建立
2.1 数值模拟理论基础
本文采用分阶段弹塑性求解法、容重增加法与强度折减法的结合以及非线性动力反应分析法对边坡天然工况、降雨工况、地震工况及降雨-地震耦合工况下的应力应变特征进行研究。
2.1.1 分阶段弹塑性求解法
设置摩尔-库伦模型为材料的本构模型,先设置粘聚力和抗拉强度为大值,计算至平衡后,再设置粘聚力和抗拉强度为分析所采用的值求解,生成初始的应力场[8]。
2.1.2 强度折减法
首先将边坡强度参数粘聚力c和内摩擦角ϕ同时除以折减系数F,得到一组新的粘聚力和内摩擦角值(c′和ϕ′),然后,将c′和ϕ′值作为一组新的材料参数输入计算;不断增加F值,直到边坡达到极限状态,发生失稳破坏,这时所对应的F值即为边坡的稳定性系数[9]。用公式表示如下:
tanϕ′=(tanϕ)F (1) c′=cF (2) 式中:F为折减系数;c为岩土边坡最初的粘聚力,kPa;ϕ为岩土边坡初始内摩擦角,(°);c′、ϕ′分别为经过折减后的粘聚力和内摩擦角。
2.1.3 容重增加法
保持岩土的粘聚力及内摩擦角不变,逐步增加重力加速度G反复进行计算(逐步增加重力加速度等同于增加岩土体的容重),直至边坡破坏,此时采用的重力加速度Glimit与实际重力加速度(G0,通常取值为9.8 m/s2)之比即为该边坡的安全系数[10],即:
Fs+Glimit/G0 (3) 2.2 模型尺寸及结构面划分
研究所采用的边坡模型(见图 1)长180 m,宽80 m,高110 m。平均坡度降约0.48。坡体中上部易滑坡区域的坡度约36°。为便于分析,沿易滑坡区域对称轴取剖面1-1,易滑坡区对称轴与正北方向呈51°角。模型共由4个结构面组成。结构面Ⅰ是砂岩、板岩夹灰岩与上层的块石土以及碎石土的分界面,倾角40°;结构面Ⅱ是碎石土与块石土的分界面,倾角16°;结构面Ⅲ是块石土与粉质黏土的分界面,倾角34°;结构面Ⅳ是粉质黏土与上层覆土的分界面,倾角38°。边坡岩土物理力学参数(见表 1)根据现场勘测所获取资料确定。
表 1 边坡岩石力学参数Table 1. Natural and rainfall condition parameters工况 岩土体类型 D/(kg·m-3) c/kPa ϕ/(°) E/GPa μ σt/MPa 安全系数 天然 砂岩、板岩夹灰岩 2400 3.30E+07 55.0 2.30E+10 0.26 5.00E+05 1.31875 块石土 2200 1.05E+04 40.0 8.20E+08 0.28 0.00E+00 碎石土 2100 1.11E+04 40.0 8.10E+08 0.23 0.00E+00 粉质黏土 1900 2.20E+04 20.0 1.00E+08 0.35 0.00E+00 25 mm/d(降雨量) 砂岩、板岩夹灰岩 2450 2.81E+07 46.7 2.30E+10 0.26 4.25E+05 1.00125 块石土 2320 8.93E+03 34.0 8.20E+08 0.28 0.00E+00 碎石土 2210 9.44E+03 34.0 8.10E+08 0.23 0.00E+00 粉质黏土 2030 1.87E+04 17.0 1.00E+08 0.35 0.00E+00 35 mm/d(降雨量) 砂岩、板岩夹灰岩 2470 2.60E+07 43.8 2.30E+10 0.26 3.93E+05 0.93245 块石土 2340 8.29E+03 31.4 8.20E+08 0.28 0.00E+00 碎石土 2230 8.62E+03 31.4 8.10E+08 0.23 0.00E+00 粉质黏土 2100 1.73E+04 15.7 1.00E+08 0.35 0.00E+00 45 mm/d(降雨量) 砂岩、板岩夹灰岩 2485 2.39E+07 39.9 2.30E+10 0.26 3.73E+05 0.91225 块石土 2357 7.58E+03 29.2 8.20E+08 0.28 0.00E+00 碎石土 2247 8.11E+03 29.2 8.10E+08 0.23 0.00E+00 粉质黏土 2148 1.62E+04 15.1 1.00E+08 0.35 0.00E+00 注:D-容重;c-粘聚力;ϕ-内摩擦角;E-弹性模量;μ-泊松比;σt-抗拉强度 2.3 控制条件设置
2.3.1 天然工况
天然工况即自重工况。赋予摩尔-库仑破坏准则,采用分阶段弹塑性求解法生成初始地应力场,利用自编强度折减法计算得出天然工况下边坡的安全系数。
2.3.2 降雨工况
降雨为该地区的主要工况。在天然工况的基础上,采用强度折减法的思想对岩土体的抗剪强度参数进行折减,其后利用自编强度折减法计算得出不同降雨条件下的边坡安全系数。
2.3.3 地震工况
研究边坡位于岷江地震带,地震频发。采用弹塑性模型以及摩尔-库仑破坏准则。为减少变形过大对结果的影响,体系最大不平衡力设定为100,时间步长设置为1~5,采用最小临界阻尼比0.005、最小中心频率12.8的瑞利阻力。在模型底部导入峰值加速度为3 m/s2的汶川地震波,其达到峰值的时间为9.97 s,持续时间为20 s。
2.3.4 降雨-地震耦合工况
考虑研究区特殊的环境条件,综合强度折减法和容重增加法,在降雨工况的基础上施加地震动力荷载,比较分析在不同降雨量和地震共同作用下边坡的应力应变特征及其稳定性。各耦合工况设置情况见表 2。
表 2 耦合工况设置Table 2. Multi load-case set工况 降雨条件/(mm·d-1) 地震条件 耦合工况一 25 汶川地震波 耦合工况二 35 汶川地震波 耦合工况三 45 汶川地震波 3. 各工况计算及结果分析
3.1 天然工况及降雨工况
模拟过程中,考虑不同降雨情况对边坡稳定性的影响,设置25 mm/d、35 mm/d、45 mm/d等3种降雨情况,对摩擦角、粘聚力、体积模量、剪切模量等参数分别进行0.85、0.78、0.72倍折减。同时,考虑不同降雨情况下雨水入渗的加载作用,结合容重增加法,对土体容重进行不同情况的扩大。各降雨工况岩土体参数与其安全系数见表 1。
由表 1可知,自然工况下边坡的安全系数为1.318,有较高的安全储备,边坡处于稳定状态。降雨量为25 mm/d时,安全系数为1.001,较天然工况降低24.23%,安全储备所剩无几,边坡处于临界滑坡状态,微小的不利影响即可启动滑坡。降雨量为35 mm/d时,安全系数为0.932,较天然降雨工况降低29.29%,安全储备为0,此时边坡的抗滑力小于下滑力,边坡发生破坏。降雨量为45 mm/d时,安全系数为0.912,较天然降雨工况降低30.80%,安全储备为0,较35 mm/d降雨工况,滑体产生更大的位移。
图 2显示了降雨工况下安全系数随降雨强度的变化,由图 2可知,当降雨量为0时,边坡处于稳定状态;随着降雨量的增加,边坡稳定性逐渐降低,且当降雨量小于25 mm/d时降低幅度较大,降雨量大于25 mm/d时降低幅度较小。该曲线的拟合方程为:
y=1.87664E−4x2−0.01753x+1.31933 (4) 图 3为各降雨工况下边坡的位移分布图。从图 3可知,不同降雨工况下的位移高值区均分布于坡腰中上部,但范围随降雨量增加逐步扩大。降雨量为25 mm/d时,最大位移值为39.3 mm(见图 3a),位移高值区分布于高程2683~2714 m范围内。降雨量为35 mm/d时,最大位移值为341.7 mm(见图 3b),位移高值区分布于高程2682~2719 m范围内。降雨量为45 mm/d时,最大位移值为528.1 mm(见图 3c),位移高值区分布于高程2681~2733 m范围内。
分析可知,降雨工况与天然工况相比,位移高值区分布大致相同,随降雨量增加,雨水对土体的软化和加载作用增强,位移高值区范围沿坡面逐步扩大,且边坡稳定性逐渐降低。当降雨量达到35 mm/d时,滑体位移由39.3 mm突变至341.7 mm,表征其在雨水入渗和加载的综合作用下发生局部滑坡,这与实际工程情况较吻合。
3.2 地震工况
数值模拟模型采用弹塑性模型及摩尔-库仑准则。模型阻尼采用瑞利阻尼,其最小临界阻尼比取0.005,最小中心频率取12.8。施加汶川地震动力作用,汶川地震波的峰值加速度为3 m/s2,持时为20 s。汶川地震波加速度时程曲线见图 4。
在模型的关键位置布置4个监测点,具体位置(见图 5)为:A点位于高程2740 m的坡肩处,B点位于高程2723 m的上坡腰处,C点位于高程2670 m的下坡腰处,D点位于高程2627 m的坡脚处,记录各监测点的加速度、位移(见图 6、图 7)。
由图 6可知,边坡坡面各监测点加速度响应差异较大。坡脚D点加速度峰值最小,几乎为0;坡肩A点其次,加速度在±0.6 m/s2范围内变动;坡腰C点的加速度在±1.23 m/s2范围内变动;坡腰B点的加速度峰值最大,加速度在-3.3 m/s2~+1.4 m/s2范围内变动。
由图 7可知,坡脚D点位移几乎为0;坡肩A点在地震波持时3.5 s时开始移动,11 s时位移发生突变,14 s时位移达到最大值7.3 cm;坡腰C点在3.5 s时开始移动,12 s时位移发生突变,17 s时位移达到最大值-11.5 cm;坡腰B点在1 s时开始移动,12 s时位移发生突变,18 s时达到最大值-57 cm。
分析认为,在地震波持时18 s时各监测点位移均已达到稳定状态,表征各监测点在汶川地震作用下产生塑性变形。坡腰C点所在土层(粉质黏土)与坡腰B点所在土层(砂岩、板岩夹灰岩)在地震动力作用下的位移有极大差异,表征发生较大相对位移,即C点处土层发生局部破坏,边坡处于不安全状态。
3.3 降雨与地震耦合工况
由各耦合工况坡腰B点加速度图(见图 8)可知,从0 s开始各工况坡腰B点就有较大的加速度,其中耦合工况一最小,耦合工况二其次,耦合工况三最大。耦合工况一坡腰B点处加速度从0 s时的-1.27 m/s2锐减到0,然后在-1.1~1.08 m/s2范围内变动;耦合工况二坡腰B点处加速度从0 s时的4.71 m/s2锐减到0,然后在-2.01~3.15 m/s2范围内变动;耦合工况三坡腰B点处加速度从0 s时的-7.23 m/s2锐减到0,然后在-4.02~3.15 m/s2范围内变动。
图 9显示了各耦合工况下坡腰B点的位移。由图 9可知,耦合工况一、二、三的坡腰B点位移在0~17 s的范围内大致呈线性增加。耦合工况一坡腰B点处的位移在17.8 s时达到最大值1.9 m,之后趋于稳定;耦合工况二坡腰B点处的位移在18.3 s时达到最大值5.2 m,之后趋于稳定;耦合工况三坡腰B点处的位移在18.5 s时达到最大值13.6 m,之后趋于稳定。
由降雨与耦合工况坡腰B点位移对比图(见图 10)可知,降雨-地震耦合工况下坡腰的位移明显比降雨工况大。降雨量为25 mm/d时,坡腰位移在降雨工况下为0.039 m,在耦合工况下为1.9 m;降雨量35 mm/d时,坡腰位移在降雨工况下为0.34 m,耦合工况下为5.2 m;降雨量45 mm/d时,坡腰位移在降雨工况下为0.53 m,在耦合工况下为13.6 m。
分析认为,降雨工况使得边坡土体的容重和含水率增加,粘聚力和内摩擦角降低,地震动力作用效应得到了放大,导致在降雨-地震耦合工况下,边坡发生破坏,其位移量达到单一降雨工况和单一地震工况时的数十倍。
4. 结论
结合工程实例,运用FLAC3D有限差分软件详细地对成兰铁路岷江双线特大桥高边坡进行建模,并模拟、分析了天然、降雨、地震以及降雨-地震耦合等4种工况下边坡的位移、加速度等变化特征,对坡体的稳定性进行安全性评价。
通过天然工况和降雨工况坡体稳定性分析,得到天然工况和不同降雨工况的强度折减安全系数表以及降雨量与安全系数关系式,表明随降雨量增加,边坡稳定性逐步降低,为量化分析降雨量与滑坡的关系提供基础数据。
通过对比不同降雨工况与不同耦合工况,得到降雨工况和耦合工况的关系曲线,表明降雨和地震的共同作用,会大大削弱边坡的稳定性,使边坡稳定性较降雨和地震单独作用时要低,极其容易造成滑坡。故须对类似特殊环境的边坡采取合适的工程措施,增加其稳定性。
较详细地对成兰铁路岷江双线特大桥高边坡进行数值模拟分析,为单工况以及耦合工况下边坡的稳定性研究提供理论依据。同时,也对大巴山、龙门山、龙泉山等其他滑坡频发山区的研究和预警提供借鉴。
-
表 1 典型金矿区矿体特征
Table 1. Characteristics of ore-bodies in the typical gold mine areas
矿区 矿体编号 矿体形态 矿体长度/m 矿体厚度/m 平均品位/10-6 矿体产状 主要矿石类型 金斗坡 K1 透镜状 500 1.00~2.97 1.30 280°—330°∠29°—43° 黑云母变斑晶绢云石英片岩含炭绢云石英片岩 K2 似层状 250 1.90~2.74 1.18 K3 似层状 250 1.90~3.95 1.13 黄龙 Ⅲ 似层状 150 1.00~6.31 5.60 5°—22°∠20°—55° 黑云母变斑晶绢云石英片岩含炭绢云石英片岩 Ⅳ 似层状 500 0.80~7.14 3.15 Ⅷ-1 似层状 200 0.97~1.84 4.01 长沟 K1 似层状 460 0.85~3.80 1.54 190°—220°∠35°—55° 黑云母变斑晶绢云石英片岩含炭绢云石英片岩夹变石英砂岩含金石英脉 K2 似层状 190 0.95~4.45 1.52 K4 似层状 1100 0.85~4.35 2.60 15°—35°∠40°—70° 注:部分资料据文献[10] 表 2 不同类型变质作用阶段划分及矿物生长序列[7]
Table 2. Dividing of different phases on metamorphism and the order of mineral development
-
[1] 吕古贤, 李秀章, 张迎春, 等.矿田地质学的研究和发展问题[J].地质与勘探, 2012, 48(6):1143~1150. http://www.cnki.com.cn/Article/CJFDTOTAL-DZKT201206012.htmLV Gu-xian, LI Xiu-zhang, ZHANG Ying-chun, et al. Research and development of ore field geology[J]. Geology and Exploration, 2011, 48(6): 1143~1150. http://www.cnki.com.cn/Article/CJFDTOTAL-DZKT201206012.htm [2] 吴淦国, 吕承训.矿田构造的研究历史, 现状与发展[J].地质通报, 2011, 31(4):461~468. http://www.cnki.com.cn/Article/CJFDTOTAL-ZQYD201104003.htmWU Gan-guo, LV Cheng-xun. Research history, present situation and development of ore field structure study[J]. Geological Bulletin of China, 2011, 30(4): 461~468. http://www.cnki.com.cn/Article/CJFDTOTAL-ZQYD201104003.htm [3] 刘林, 吴景平, 杨译, 等.陕西石泉—汉阴北部下志留统黑色岩系中的金矿床[J].地质与资源, 2011, 20(4):304~309. http://www.cnki.com.cn/Article/CJFDTOTAL-GJSD201104020.htmLIU Lin, WU Jing-ping, YANG Yi, et al. Gold deposit in the Lower Silurian black series in northern Shiquan-Hanyin, Shaanxi Province[J]. Geology and Resources, 2011, 20(4): 304~309. http://www.cnki.com.cn/Article/CJFDTOTAL-GJSD201104020.htm [4] 李福让, 侯俊富, 杨弘, 等.陕西石泉县羊坪湾金矿床地质-地球化学特征及金的富集规律[J].中国地质, 2009, 36(1):174~182. http://www.cnki.com.cn/Article/CJFDTOTAL-DIZI200901019.htmLI Fu-rang, HOU Jun-fu, YANG Hong, et al. Geological and geochemical characteristics and gold enrichment regularity of the Yangpingwan gold deposit in Shiquan, Shaanxi Province[J]. Geology in China, 2009, 36(1): 174~182. http://www.cnki.com.cn/Article/CJFDTOTAL-DIZI200901019.htm [5] 唐永忠, 杨兴科, 张宝荣, 等.南秦岭造山带安康石梯—旬阳神河早古生代热水沉积盆地构造-沉积相与热水聚矿特征[J].中国地质, 2012, 39(5):1261~1270. http://www.cnki.com.cn/Article/CJFDTOTAL-DIZI201205012.htmTANG Yong-zhong, YANG Xing-ke, ZHANG Bao-rong, et al. Structure-sedimentary facies and hydrothermal mineralization characteristics of Shiti (in Ankang)-Shenhe (in Xunyang) Early Paleozoic hydrothermal basin in Qinling orogenic belt[J]. Chinese Geology, 2012, 39(5): 1261~1270. http://www.cnki.com.cn/Article/CJFDTOTAL-DIZI201205012.htm [6] 杜坚.安康北部志留系金矿地质特征及找矿前景探讨[J].中国矿山工程, 2014, 43(5):39~43. http://www.cnki.com.cn/Article/CJFDTOTAL-YSKS201405012.htmDU Jian. Discussion on the geological characteristics and prospecting potential of Silurian gold-bearing in north Ankang[J]. China Mine Engineering, 2014, 43(5): 39~43. http://www.cnki.com.cn/Article/CJFDTOTAL-YSKS201405012.htm [7] 张康. 南秦岭汉阴北部志留系金矿田构造控矿规律及找矿标志研究[D]. 西安: 长安大学, 2012.ZHANG Kang. Ore-controlling rules and prospecting indicator of Silurian gold deposit in North Hanyin, South Qinling[D]. Xi'an: Chang'an University, 2012. [8] 张国伟, 张宗清, 董云鹏, 等.秦岭造山带主要构造岩石地层单元的构造性质及其大地构造意义[J].岩石学报, 2009, 11(2):101~114. http://www.cnki.com.cn/Article/CJFDTOTAL-YSXB199502000.htmZHANG Guo-wei, ZHANG Zong-qing, DONG Yun-peng, et al. Nature of main tectono-lithostratigraphic units of the Qinling orogen: Implications for the tectonic evolution[J]. Acta Petrologica Sinica, 1995, 11(2): 101~114. http://www.cnki.com.cn/Article/CJFDTOTAL-YSXB199502000.htm [9] 杨兴科, 何虎军, 张康, 等. 陕西安康汉阴县铁佛寺—汉滨区流芳志留系金矿田构造及快速评价勘查方法[R]. 西安: 长安大学, 2012.YANG Xing-ke, HE Hu-jun, ZHANG Kang, et al. Geologic structure of Silurian gold ore-field and instant method of assessment in Tiefosi, Hanyin-Liufang, Hanbin, in Ankang, Shaanxi Province[R]. Xi'an: Chang'an University, 2012. [10] 吴闻人, 王北颖. 陕西省南秦岭造山带中部韧性剪切带的形成、演化及其与金矿成矿关系研究[R]. 西安: 陕西省地矿局, 1991.WU Wen-ren, WANG Bei-ying. Formation and development of ductile shear belt in the middle of South Qinling orogenic belt, Shaanxi Province and its relation to gold mineralization[R]. Xi'an: Geology and Mineral Resources Bureau of Shaanxi, 1991. [11] 李会民.石泉—汉阴北部金矿带地质特征[J].陕西地质, 1997, 15(2):48~57. http://www.cnki.com.cn/Article/CJFDTOTAL-SXDY199702005.htmLI Hui-min. The geological characteristics of gold deposit belt in the northern Shiquan-Hanyin, Shaanxi[J]. Geology of Shaanxi, 15(2): 48~57. http://www.cnki.com.cn/Article/CJFDTOTAL-SXDY199702005.htm [12] 徐开礼, 朱志澄.构造地质学[M].北京:地质出版社, 1987:1~200.XU Kai-li, ZHU Zhi-cheng. Structural geology[M]. Beijing: Geological Publishing House, 1987: 1~200. [13] 周喜文, 魏春景, 卢良兆.高温变泥质岩石中石榴石-黑云母地质温度计的应用——以胶北荆山群富铝岩石为例[J].地学前缘, 2004, 10(4):353~363. http://www.cnki.com.cn/Article/CJFDTOTAL-DXQY200304005.htmZHONG Xi-wen, WEI Chun-jing, LU Liang-zhao, et al. Application of the garnet-biotite geothermometer in high grade metapelite: Al-rich rock from the Jingshan Group in north Jiaodong, China[J]. Earth Science Frontiers, 2003, 10(4): 353~363. http://www.cnki.com.cn/Article/CJFDTOTAL-DXQY200304005.htm [14] 吴春明, 陈泓旭.变质作用温度与压力极限值的估算方法[J].岩石学报, 2013, 29(5):1499~1510. http://www.cnki.com.cn/Article/CJFDTOTAL-YSXB201305006.htmWU Chun-ming, CHEN Hong-xu. Estimation of minimum or maximum pressure or temperature conditions in metamorphic rocks[J]. Acta Petrologica Sinica, 2013, 29(5): 1499~1510. http://www.cnki.com.cn/Article/CJFDTOTAL-YSXB201305006.htm [15] 何德锋, 钟宏, 朱维光.石榴石-黑云母地质温度计在四川拉拉铜矿床的应用[J].矿物学报, 2008, 28(2):127~134. http://www.cnki.com.cn/Article/CJFDTOTAL-KWXB200802004.htmHE De-feng, ZHONG Hong, ZHU Wei-guang. Appication of the garnet-biotite geothermometer in the Lala copper deposit, Sichuan Province[J]. Acta Mineralogica Sinica, 2008, 28(2): 127~134. http://www.cnki.com.cn/Article/CJFDTOTAL-KWXB200802004.htm [16] Perchuk L L. Thermodynamical control of metamorphic processes[C]//Saxena S K, Bhattacharji S. Energetics of geological processes. New York: Springer-Verlag, 1977: 285~352. [17] JIAO Jian-gang, YANG Xing-ke, QIU Dong-dong, et al. Report on Silurian structure-lithofacies mapping and metallogenic regularity of gold deposit in Shiquan County, Shaanxi[J]. Xi'an: Chang'an University, 2014. [18] 王友文, 李兴柱, 阎竹斌, 等.中国金矿地质及普查勘探方法[M].西安:陕西地质矿产局, 1987:56~59.WANG You-wen, LI Xing-zhu, YAN Zhu-bin, et al. Method of investigation and exploration on Chinese gold mines[M]. Xi'an: Shaanxi Bureau of Geology and Mineral Resources, 1987: 56~59. [19] 隗合明. 安康北部志留系岩金矿品位指标研究论证报告[R]. 西安: 长安大学, 2006.WEI He-ming. Research report of Silurian gold grade index in north Ankang[R]. Xi'an: Chang'an University, 2006. [20] 任小华, 高航, 李福让, 等.羊坪湾金矿地质特征及成矿模式[J].西北地质, 2015, 48(1):127~136. http://www.cnki.com.cn/Article/CJFDTOTAL-XBDI201501017.htmREN Xiao-hua, GAO Hang, LI Fu-rang, et al. Geological features and metallogenic model of Yangpingwan gold deposit[J]. Northwestern Geology, 2015, 48(1): 127~136. http://www.cnki.com.cn/Article/CJFDTOTAL-XBDI201501017.htm [21] 高怀雄, 李卫波, 许锋, 等. 陕西石泉—旬阳金矿成矿规律研究与矿产资源预测[R]. 西安: 陕西省地质调查中心, 2016.GAO Huai-xiong, LI Wei-bo, XU Feng, et al. Research on gold metallogy and prediction of mineral resources in Shiquan-Xunyang area, Shaanxi[R]. Xi'an: Geological Survey Center of Shaanxi Province, 2016. [22] 赵进江, 刘洪祺, 杨鹏, 等.陕西省汉阴县长沟金矿地质和地球化学特征[J].科技信息, 2013, 23:399~401. doi: 10.3969/j.issn.1001-9960.2013.21.320ZHAO Jin-jiang, LIU Hong-qi, YANG Peng, et al. Alteration and geochemical features of Changgou gold ore in Hanyin County, Shaanxi Province[J]. Science and Technology Information, 2013, 23: 399~401. doi: 10.3969/j.issn.1001-9960.2013.21.320 -