A discussion on the calculation method of instability probability of landslide due to rainfall
-
摘要:
为了解决滑坡风险评价中的滑坡失稳概率计算问题，利用前人在降雨阈值的研究成果，结合气象学中降雨概率分布理论，以云南省盐津县庙坝滑坡为例进行计算，建立降雨型滑坡失稳概率计算模型。结果表明，盐津县降雨型滑坡的降雨阈值类型为累积降雨量-历时关系阈值，即为单日降雨阈值，降雨阈值为29.7 mm；盐津县在当日降雨量达到或超过阈值水平时可能诱发滑坡，对滑坡影响的滞后天数最大为5天；庙坝滑坡在8月20—25日6天内单日降雨达到或超过29.7 mm的降雨概率为46.49%；庙坝滑坡在8月25日因前5天或当天单日降雨量超过29.7 mm而失稳的概率为0.2853%。
Abstract:In order to solve the landslide risk evaluation of landslide failure instability probability calculation problem, this paper summarizes the research achievements of previous researchers in rainfall threshold and, in combination with the previous research results in a threshold rainfall, uses the theory of statistical to perform coupling analysis of the historical record of the rainfall landslide and rainfall data, so as to create a rainfall landslide failure instability probability calculation model, with Miaoba landslide in Yanjin County of Yunnan Province as an example for demonstration.The results show that the rainfall threshold type of Yanjin County rainfall landslide is cumulative rainfall duration threshold, which means the one-day rainfall threshold, with the rainfall threshold being 29.7 mm; most of rainfall landslide in Yanjin County is caused by the daily rainfall which reaches rainfall threshold within 5 days; the possibility that the rain reaches or exceeds 29.7 mm from August 20th to 25th is 46.49%;the possibility of the failure of Miaoba landslide in August 25th is 0.2853%.
-
川-滇-黔铅锌银多金属成矿区是中国主要的铅锌银矿产基地之一，共发现400多个铅-锌矿床和矿化点。前人主要针对震旦系—二叠系沉积岩中的典型铅锌矿床进行了较系统的矿床地质、矿床地球化学及成矿年代学研究，取得了诸多重要进展[1-4]：① 矿床分布明显受构造控制，属于后生热液矿床；② 成矿物质具有多来源特征，其中基底岩石是重要的物源岩；③ 尽管成矿年龄差别较大，但主要集中于192~226Ma，暗示矿床形成与印支运动有关。而滇中铅锌矿集区赋矿层位多，从昆阳群到第四系各时代地层中均有铅锌矿床赋存，几乎涵盖了所有川-滇-黔成矿区赋矿地层。云南禄劝噜鲁铅锌矿床滇中矿集区赋存于下寒武统梅树村组顶部白云岩地层中。前人仅对滇中矿集区做过报道[5-10]，对云南禄劝噜鲁铅锌矿床研究较薄弱，地球化学资料未见报道，制约了该矿床成矿过程和成矿预测研究。
本次在详实的矿床地质研究基础上，开展系统的同位素地球化学和年代学研究，以期查明噜鲁铅锌矿床的成矿特征、物质来源和形成时代，分析矿床形成机制及成矿动力学背景。
1. 区域及矿床地质概况
滇中铅锌矿集区位于扬子地块西南缘，上扬子铅锌成矿带上，其大地构造位置大致北起东川东西向的宝台厂-洪门厂断裂，南到红河断裂西至绿汁江断裂，东达小江断裂。滇中铅锌矿集区地壳结构复杂，具有多层结构的特征。滇中铅锌矿集区地层从古元古界大红山群（苴林群）到新生界第四系均有分布。沉积岩岩性以碳酸盐岩为主。滇中铅锌矿集区经历了多期的构造运动，断裂构造主要有元谋-绿汁江断裂、普渡河-滇池断裂、西昌-易门断裂、小江断裂、红河断裂、弥勒-师宗断裂（图 1）。
![]()
图 1 噜鲁铅锌矿床矿区地质图P2β—上二叠统峨眉山玄武岩组；P1q+m—下二叠统茅口组+栖霞组；∈1d+s—中寒武统双龙潭+陡坡寺组；∈1l—下寒武统龙王庙组；∈1c—下寒武统仓浪铺组；∈1q—下寒武统筇竹寺组；∈1m2—下寒武统梅树村组上段；∈1m1—下寒武统梅树村组下段；Zbdn—上震旦统灯影组；1—城市名；2—断裂（①—元谋-绿汁江断裂；②—西昌-易门断裂；③—普渡河-滇池断裂；④—小江断裂；⑤—红河断裂；⑥—弥勒-师宗断裂）；3—湖泊; 4—铅锌矿体；5—断层；6—地层产状；7—地理点；8—地质界线；9-勘探线Figure 1. Geological map of the Lulu Pb-Zn ore district滇中矿集区受扬子板块与印度板块碰撞及板内攀西裂谷作用的影响，岩浆活动频繁，从古—中生代均有岩浆活动，最重要的是二叠纪基性火山活动，其次为燕山期中酸性岩浆侵入活动。
噜鲁铅锌矿床位于昆明市禄劝县，大地构造位置地处扬子地块西南边缘，小江深大断裂与普渡河-滇池深大断裂之间，目前该矿床仍处于勘查阶段，资源量已达小型矿床规模。
矿区出露地层有新元古界震旦系灯影组（Zbdn）、下寒武统梅树组（∈1m）、筇竹寺组（∈1q）、沧浪铺组（∈1c）、龙王庙组（∈1l）、陡坡寺+双龙潭组（∈1d+s）、下二叠统茅口+栖霞组（P1q+m）、上二叠统峨眉山玄武岩组（P2β）。铅锌矿体赋存于下寒武统梅树组下段（∈1m1）（图 1）。
矿区内构造简单，为宏宽背斜东翼的单斜岩层构造，走向北北西，倾向约70°，倾角多为12°~20°，局部出现小的次级褶曲。在矿区范围内断层不发育，仅南部有1条，且离铅锌矿体较远，对矿体没有明显的破坏作用。
矿区内岩浆活动较强，主要为中上二叠统（P2β），形成上部、下部2个旋回。第一旋回岩性主要以灰绿色斜斑玄武岩为主夹致密玄武岩或斜斑玄武岩。第二旋回岩性主要为底部火山角砾岩、火山角砾玄武质凝灰岩、凝灰岩夹灰岩；向上为熔岩，上部为杏仁状玄武岩；中部致密状玄武岩为主夹杏仁状玄武岩；下部以斜长玄武岩为主。
矿区目前圈定铅锌矿体1个（图 2），矿体赋存于下寒武统梅树村组下段（∈1m1）顶部白云岩中，矿体产状与地层基本一致，呈似层状、局部地段有脉状矿体。矿体地表出露长约150m，倾向南，倾角15°~ 20°，矿体厚度一般为1.70~10.71m，平均厚约3.26m，随深度增大，矿体厚度有逐步变厚的趋势。矿体倾向延伸控制约350m。随着矿体的延深见少量的白云岩夹石呈透镜状分布。矿体Pb品位为0.44%~ 20.00%，平均品位为1.78%；Zn品位为0.56%~9.06%，平均品位为2.01%。
![]()
图 2 噜鲁铅锌矿床3勘探线剖面∈1c—下寒武统仓浪铺组；∈1q—下寒武统筇竹寺组；∈1m2—下寒武统梅树村组上段；∈1m1—下寒武统梅树村组下段；1—白云岩；2—砂质白云岩；3—泥质粉砂岩；4—粉砂岩；5—页岩；6—铅锌矿体；7—产状；8—钻孔Figure 2. Geological section along No. 3 exploration line of the Lulu Pb-Zn deposit矿石中金属矿物主要为方铅矿、闪锌矿、黄铁矿、白铁矿、毒砂等；非金属矿物主要为重晶石、石英、方解石。矿石中常见菊花状、放射状重晶石（图 3），也可见少量乳滴状沥青。矿石结构主要为自形-他形中细晶状结构、压碎结构、纤维状结构、交代溶蚀结构等。矿石构造主要为块状构造、网脉状构造、浸染状构造、脉状构造、条带状-条纹状构造。围岩蚀变主要有黄铁矿化、硅化、重晶石化、碳酸盐化等。其中黄铁矿化常见黄铁矿呈团块状、星散状、脉状等分布，脉宽最大可达1m，在坑道及老硐见矿位置均普遍分布。重晶石化常见重晶石呈团块状、放射状及细脉状分布于铅锌矿体及黄铁矿周围，较富集。碳酸盐化以白云石化为主，广泛发育于含矿层及顶板、底板，并贯穿于成矿过程的始终。
![]()
图 3 噜鲁铅锌矿床矿物手标本与镜下特征a—块状矿石，矿石矿物组合为方铅矿(Gn)、黄铁矿(Py)、重晶石(Brt)；b—方铅矿、闪锌矿(Sp)、重晶石等矿物呈现星点状分布；c—黄铁矿沿闪锌矿粒间充填交代，呈不规则状分布；d—方铅矿穿插交代闪锌矿，少量黄铁矿包含于闪锌矿中Figure 3. Mineral hand specimens and microscopic features of Lulu Pb-Zn deposit矿区经历了多期成矿作用，根据矿脉穿插关系、矿物组合及矿石结构，结合区域构造演化过程，将噜鲁铅锌的成矿作用划分为沉积成岩期、热液成矿期和表生氧化期3期，其中热液成矿期可进一步划分为3个成矿阶段，即重结晶白云岩-石英-黄铁矿阶段、黄铁矿-方铅矿-闪锌矿阶段（主成矿阶段）和黄铁矿阶段。
2. 样品及分析方法
样品采自噜鲁铅锌矿床矿区坑道内，新鲜。用于Rb-Sr同位素定年的黄铁矿与方铅矿为矿物共生组合。
Rb-Sr化学分离与质谱测试在中国科学院地质与地球物理研究所固体同位素地球化学实验室完成。同位素比值测试在高精度固体热电离质谱计（IsoProbe-T）上完成。采用单W和单Ta灯丝载入Sr和Rb样品，以纯化的TaF5作为发射剂，以Faraday接收器静态测量锶同位素组成，以Daly接收器监测Rb对Sr的同质异素干扰。测定方法和仪器参数见参考文献[9-10]，锶同位素比值采用保守外精度误差0.05‰。
硫同位素组成分析在中国科学院贵阳地球化学研究所环境地球化学国家重点实验室，采用EAIRMS法在连续流质谱仪上完成。用国标GBW04415和04414，Ag2S做内标，以CDT为标准，分析精度约为0.2‰（2σ）。
铅同位素样品测试在核工业北京地质研究院分析测试研究中心进行，先称取适量样品放入聚四氟乙烯坩埚中，加入氢氟酸中、高氯酸溶样。样品分解后，将其蒸干，再加入盐酸溶解蒸干，加入0.5N氢溴酸（HBr）溶液溶解样品进行铅的分离；将样品溶解倒入预先处理好的强碱性阴离子交换树脂中进行铅的分离，用0.5N氢溴酸（HBr）溶液淋洗树脂，再用2N盐酸（HCl）溶液淋洗树脂，最后用6N盐酸（HCl）溶液解脱，将解脱溶液蒸干备质谱测定；用热表面电离质谱法进行铅同位素测量，仪器型号为JY/T004-1996，对1μg的铅，208Pb/206Pb测量精度优于0.005%，NBS981标准值（2δ）：208Pb/206Pb= 2.1681 ± 0.0008，207Pb/206Pb=0.91464 ± 0.00033，204Pb/ 206Pb=0.059042±0.000037。
3. 测试结果
3.1 Rb-Sr同位素组成
本次对噜鲁铅锌矿床中的5件硫化物样品进行测试。根据测试结果可知，噜鲁铅锌矿床方铅矿单矿物Rb、Sr含量存在变化。由表 1可见，全部硫化物样品的Rb含量较低（0.1476×10-6~2.013×10-6），Sr含量不高（0.2279 × 10-6~2.957 × 10-6），Rb/Sr值为0.1142~4.0158；87Rb/86Sr值为0.3205~12.17，87Sr/86Sr值为0.712331~0.746457。
表 1 噜鲁铅锌矿床硫化物Rb-Sr同位素组成Table 1. Rb-Sr isotopic composition of sulfides from the Lulu Pb-Zn deposit编号 矿物 Rb/10-6 Sr/10-6 87Rb/86Sr 87Sr/86Sr LL-0-02 黄铁矿 0.1476 1.293 0.3205 0.712331±10 LL-0-04 黄铁矿 2.013 2.957 2.036 0.717113±8 LL-13 方铅矿 0.3405 0.9726 1.037 0.714388±9 LL-34 方铅矿 1.594 1.231 3.814 0.722339±11 LL-36 方铅矿 0.9152 0.2279 12.17 0.746457±8 3.2 硫同位素组成
本次共测试噜鲁铅锌矿床硫化物和硫酸盐矿物样品40件，硫同位素组成见表 2。其中13件方铅矿的δ34S值介于10.32‰~12.82‰之间，平均值为11.32‰；6件闪锌矿的δ34S值介于9.72‰~18.76‰之间，平均值为12.87‰；12件黄铁矿的δ34S值介于12.52‰~22.44‰之间，平均值为17.29‰；9件重晶石的δ34S值介于26.08‰ ~29.34‰之间，平均值为27.29‰；矿石中硫化物δ34S的变化规律为：方铅矿δ34S＜闪锌矿δ34S＜黄铁矿δ34S＜重晶石δ34S，说明热液成矿作用过程中硫同位素分馏达到平衡。硫同位素直方图浓度中心为9‰~15‰（图 4）。
表 2 噜鲁铅锌矿床硫同位素组成Table 2. Sulfur isotopic composition of the Lulu Pb-Zn deposit编号 测试矿物 δ34S/‰ 编号 测试矿物 δ34S/‰ LL-0-01 方铅矿 11.09 LL-27 黄铁矿 16.80 LL-0-02 方铅矿 10.86 LL-29 黄铁矿 13.78 LL-0-04 方铅矿 10.77 LL-32 黄铁矿 16.86 LL-12 方铅矿 11.63 LL-42 黄铁矿 22.44 LL-13 方铅矿 11.66 LL-45 黄铁矿 13.30 LL-16 方铅矿 12.10 LL-0-07 闪锌矿 14.22 LL-28 方铅矿 12.82 LL-32 闪锌矿 10.97 LL-32 方铅矿 11.08 LL-38 闪锌矿 18.76 LL-34 方铅矿 12.27 LL-0-10 闪锌矿 10.38 LL-36 方铅矿 11.01 LL-33 闪锌矿 13.18 LL-41 方铅矿 10.87 LL-0-08 闪锌矿 9.72 LL-42 方铅矿 10.32 LL-0-07 重晶石 27.02 LL-45 方铅矿 10.67 LL-34 重晶石 27.12 LL-01 黄铁矿 22.30 LL-0-08 重晶石 26.36 LL-12 黄铁矿 11.63 LL-0-04 重晶石 27.21 LL-13 黄铁矿 21.57 LL-12 重晶石 26.08 LL-14 黄铁矿 14.35 LL-0-10 重晶石 29.34 LL-14-1 黄铁矿 17.33 LL-13 重晶石 27.92 LL-16 黄铁矿 14.53 LL-16 重晶石 28.47 LL-17 黄铁矿 21.67 LL-14 重晶石 26.09 3.3 铅同位素组成
本次共测试噜鲁铅锌矿床硫化物矿物样品12件，铅同位素组成见表 3。样品206Pb/204Pb值为18.259~18.342，平均值为18.307；207Pb/204Pb值为15.608~15.639，平均值为15.623；208Pb/204Pb值为38.46~38.821，平均值为38.540；铅同位素组成相对均一。样品ω值为36.070~40.860，平均值为37.740；Th/U值为3.690~4.160，平均值为3.838；μ值在9.470~9.560之间，平均值为9.516。样品μ值均小于9.81（现代海洋沉积物的μ值）。
表 3 噜鲁铅锌矿床铅同位素组成Table 3. The lead isotope composition of the Lulu Pb-Zn ore deposits样号 矿物 206Pb/204Pb 207Pb/204Pb 208Pb/204Pb 206Pb/207Pb μ ω Th/U LL-42 黄铁矿 18.342 15.629 38.543 1.174 9.530 37.600 3.820 LL-45 黄铁矿 18.339 15.626 38.532 1.174 9.520 37.540 3.820 LL-0-01 方铅矿 18.283 15.608 38.564 1.171 9.490 37.820 3.860 LL-0-02 方铅矿 18.298 15.615 38.557 1.172 9.510 37.770 3.840 LL-0-04 方铅矿 18.315 15.624 38.561 1.172 9.520 37.780 3.840 LL-13 方铅矿 18.329 15.617 38.545 1.174 9.510 37.570 3.820 LL-32 方铅矿 18.316 15.612 38.46 1.173 9.500 37.240 3.790 LL-34 方铅矿 18.291 15.627 38.558 1.171 9.530 37.930 3.850 LL-36 方铅矿 18.331 15.618 38.537 1.174 9.510 37.530 3.820 LL-41 方铅矿 18.326 15.632 38.535 1.172 9.540 37.690 3.820 LL-42 方铅矿 18.259 15.639 38.539 1.168 9.560 38.150 3.860 LL-45 方铅矿 18.262 15.636 38.548 1.168 9.550 38.140 3.870 4. 讨论
4.1 成矿年代学及成矿构造背景
热液矿物中Rb-Sr主要赋存在矿物晶格中，或赋存在固态微包体、流体包裹体中。随着测试技术的发展，以及Rb-Sr等时线法成功地运用于MVT铅锌矿床定年，众多学者对Rb-Sr定年法进行了研究，证明了Rb-Sr等时线法对铅锌矿床定年的有效性[11-12]。同理，许多学者也运用Rb-Sr法对研究区周边铅锌矿床进行了大量的定年研究[13-15]。前人对铅锌矿床闪锌矿的Rb-Sr同位素定年研究较多，这是因为闪锌矿较其他热液矿物Rb、Sr含量较高，且Rb/Sr值较高、变化范围较大[16]。随着测试技术的发展，方铅矿和黄铁矿的Rb-Sr同位素定年也成为可能[17-19]。田世洪等[18]对玉树地区东莫扎抓铅锌矿床方铅矿与黄铁矿矿物共生组合Rb-Sr等时线年龄进行了研究，并与闪锌矿单矿物Rb-Sr等时线年龄进行比较，两者的等时线年龄在误差范围内一致。黄华等[19]对云南保山金厂河铁铜铅锌多金属矿床方铅矿与闪锌矿矿物共生组合Rb-Sr等时线年龄进行了研究，并与闪锌矿单矿物Rb-Sr等时线年龄进行比较，两者的等时线年龄在误差范围内一致。刘建明等[20]认为，由于不同矿物相具有不同的化学势，化学性质不同的Rb和Sr将发生化学分异，从而使同一成矿母液中沉淀出的一组共生矿物具有不同的Rb/Sr值，用热液矿物组合Rb-Sr等时线测定热液矿床的成矿时代比较理想。因此，笔者运用噜鲁铅锌矿床方铅矿和黄铁矿矿物组合进行Rb-Sr等时线定年。经过Isoplot软件计算，5件硫化Rb-Sr同位素初始值（87Sr/86Sr）i为0.71135± 0.00012，得到202.8±1.4Ma等时线年龄（MSWD= 0.58）（图 5），表明该矿床形成于印支晚期—燕山早期。
![]()
图 5 噜鲁铅锌矿床方铅矿与黄铁矿矿物组合Rb-Sr同位素等时线图Figure 5. Rb-Sr isochron diagram of intergrowth pyrite and galena mineral association from the Lulu Pb-Zn deposit前人对区域上其他铅锌矿床做了大量的定年工作，黄智龙[13]等利用闪锌矿单矿物颗粒Rb-Sr法测定的会泽铅锌矿床成矿年龄为226±1Ma、225±1Ma、226±7Ma；蔺志永等[14]利用闪锌矿单矿物颗粒Rb-Sr法测定跑马铅锌矿床成矿年龄为200.1±4.0Ma；包广萍等[15]利用方解石Sm-Nd法测定的茂租铅锌矿成矿年龄为196±13Ma，前人获得年代学数据与本次取得的年代学数据在误差范围内一致。
前人在研究87Sr/86Sr值的二元混合体系时提出，一组参加等时线拟合的样品，通过对其87Sr/86Sr-Sr-1、87Rb/86Sr-Rb-1作图，可以判断是否为混合线。当二者为正相关时，为混合线；二者非正相关时，获得的等时线年龄代表成矿阶段的成矿年龄[18-19]。本次将噜鲁铅锌矿床硫化物Rb-Sr同位素的87Sr/86Sr-Sr-1、87Rb/86Sr-Rb-1作图（图 6）。由图 6可知，典型铅锌矿床各要素之间不存在相关性。因此认为，硫化物Rb-Sr样品拟合的直线具有等时线意义。
![]()
图 6 噜鲁铅锌矿床方铅矿与黄铁矿矿物组合87Sr/86Sr-Sr-1、87Rb/86Sr-Rb-1关系图Figure 6. 87Sr/86Sr-Sr-1, 87Rb/86Sr-Rb-1 relationship diagrams of intergrowth pyrite and galena mineral association in the Lulu Pb-Zn deposit川滇黔成矿域C-O-S-Pb-Sr同位素研究表明，峨眉山玄武岩与铅锌矿床成矿存在联系[13]。此外，在川滇黔成矿域，幔源的碱性岩浆岩体形成时代为204~225Ma[21]。晚二叠世—三叠纪，扬子地块和与特提斯洋闭合有关的邻近板块发生碰撞，即印支运动，前人研究认为，印支晚期古特提斯洋向北不断俯冲消减，使华南板块与印支板块在早中三叠世沿松马缝合带碰撞拼合，区内随后向陆内构造体制转化[22]。年代学数据表明，川滇黔成矿域的大规模成矿主要发生在印支晚期—燕山早期，成矿作用和华南板块与印支板块拼合造山作用的时限相吻合。印支运动可能与以上岩浆活动存在联系。韩润生等[23]认为，印支期造山作用对川-滇-黔成矿区铅锌多金属成矿作用产生重要影响。刘肇昌等[24]认为，该区毗邻的龙门山于三叠纪早期开始隆升；高振敏等[25]认为，南盘江-右江在三叠纪几乎同时封闭造山。因此，龙门山造山带、南盘江-右江增生弧型冲褶带在印支期演化过程中必然对该区铅锌成矿作用产生重要影响。前人也认为，川滇黔成矿区天桥铅锌矿床逆冲断层与印支运动存在联系[26]。而噜鲁铅锌矿床成矿时代为印支运动晚期—燕山早期，因此，印支运动可能与噜鲁铅锌矿床成矿存在联系。
4.2 成矿物质来源
4.2.1 硫来源及形成机制
噜鲁铅锌矿床存在硫化物和硫酸盐共同沉淀析出现象。根据矿物组合规律推测，噜鲁铅锌矿床具有较高的氧逸度fO2，成矿流体的总硫同位素组成需综合考虑硫化物中的硫同位素及重晶石的硫同位素[27]，即δ34SΣS＞δ34S硫化物，δ34SΣS﹤δ34S重晶石。因此，下文将详细讨论该矿床的硫化物和重晶石的硫同位素组成，从而得出成矿流体中硫的来源。噜鲁铅锌矿床矿石矿物黄铁矿、闪锌矿和方铅矿相对富集重硫，其δ34S值浓度中心为9‰~15‰。明显不同于δ34S值在0附近的幔源硫。柳贺昌等[28]认为，区域上震旦系灯影组、寒武系龙王庙组及石炭系大塘组、摆佐组、马平组、黄龙组等中均有石膏、重品石等硫酸盐矿物，其δ34S值约为15%。与噜鲁铅锌矿床硫同位素组成相近，因而噜鲁铅锌矿床硫化物中的硫可能主要来自多个时代地层，为海相硫酸盐的还原产物。
噜鲁铅锌矿床重晶石的δ34S值介于26.08‰~ 29.34‰之间，平均值为27.29‰。滇中矿集区缺少下寒武统硫酸盐硫同位素数据，但Claypool等[29]报道世界其他地区寒武纪海相硫酸盐的δ34S值为30‰左右；于津生等[30]在总结前人研究成果的基础上，认为晚前寒武纪海相硫酸盐岩的δ34S值约为18‰，进入显生宙后，海相硫酸盐岩δ34S值急剧增加到30‰左右，尔后逐渐减小。笔者认为，该矿床中重晶石的硫主要来源于下寒武统海相硫酸盐。因此，硫化物与重晶石中的硫具有同源性。
海相硫酸盐的还原过程主要存在3种观点：热化学还原、细菌还原和有机质热降解。热化学还原较细菌还原发生在相对高温条件且能产生许多还原态硫，形成还原态硫的δ34S值变化范围较小，细菌还原作用发生条件与热化学还原相反[31]。热水塘铅锌矿床δ34S值集中在8‰~13‰之间。笔者测得噜鲁铅锌矿床方解石流体包裹体的均一温度主要分布在12~180℃之间（未发表数据）。前文已述，矿石沉淀时成矿流体中硫已达平衡，可利用矿物对δ34S差值计算成矿温度，利用统一手标本黄铁矿和方铅矿的δ34S差值（表 1）计算的成矿温度为122~ 227℃[32]。综合上述观点，成矿温度超过了细菌可以存活的温度[33]。因此，细菌还原作用对噜鲁锌矿床可能影响较小。有机质热降解通常发生在100~ 150℃, 虽然矿石中常发育有沥青，证明有机质在一定程度上参与成矿，但噜鲁铅锌矿床还没有由碳酸盐岩中有机质热降解作用产生硫的报道，有机质热降解作用贡献大小不能确定。这些特征均表明, 噜鲁铅锌矿床成矿流体中的硫可能为各时代碳酸盐岩地层的海相硫酸盐热化学还原的产物，有机质可能参与了成矿作用。
4.2.2 铅来源
川滇黔铅-锌-银多金属成矿域绝大部分铅锌矿床赋存于不同时代的碳酸盐岩地层的粗晶白云岩、白云质灰岩中，区域分布大面积峨眉山玄武岩。柳贺昌等[28]认为，川滇黔成矿域成矿物质由碳酸盐岩地层和玄武岩提供；李连举等[34]认为，上震旦统、下寒武统、中上泥盆统和石炭系是区域重要的矿源层；胡耀国[35]则认为，成矿物质主要来源于区域基底岩石。
前已述及，206Pb/204Pb、207Pb/204Pb、208Pb/204Pb组成相对均一。据Zartman等[36]的铅同位素构造模式图解，可将数据投在206Pb/204Pb-207Pb/204Pb图（图 7）上，数据点落入上地壳和造山带演化曲线之间，表明在成矿作用过程中，铅是多源混合的产物。该区域在基底岩石的铅同位素组成范围内（图 7），靠近不同时代碳酸盐地层和峨眉山玄武岩，与震旦系灯影组白云岩分布范围相差较远（表 4），与会泽及天桥铅锌矿床也存在差异。这些铅同位素组成特征表明，本区成矿物质主要由基底岩石提供，区域上不同时代碳酸盐岩地层和峨眉山玄武岩可能提供了部分成矿物质。
![]()
表 4 噜鲁铅锌矿床铅同位素组成统计对比Table 4. Statistics of Pb isotopic compositions of various rock units in the area统计对象 样品个数 206Pb/204Pb 207Pb/204Pb 208Pb/204Pb 会泽铅锌矿 95 18.251~18.530 15.439~15.855 38.487~39.433 天桥铅锌矿 33 18.378~18.601 15.519~15.811 38.666~39.571 栖霞~茅口组 2 18.189~18.759 15.609~16.522 38.493~38.542 黄龙组 2 15.656~16.675 18.136~18.167 38.204~38.236 摆佐组 8 18.120~18.673 15.500~16.091 38.235~39.685 大埔组 5 18.397~18.828 15.537~16.499 38.463~39.245 灯影组 10 18.198~18.517 15.699~15.987 38.547~39.271 昆阳群 27 17.781~20.993 15.582~15.985 37.178~40.483 会理群 6 18.094~18.615 15.630~15.827 38.274~38.932 峨眉山玄武岩 16 18.175~19.019 15.528~15.662 38.380~39.928 噜鲁铅锌矿床铅同位素组成相对稳定，可能是由于成矿物质来源铅同位素组成相对稳定的地质体或成矿物质来源于铅同位素组成相对不稳定的地质体，成矿流体在沉淀成矿之前存在均一化过程[13]。从表 2可见，区域上基底岩石、玄武岩及碳酸盐岩地层铅同位素组成均有较宽的变化范围，即本区不存在铅同位素组成相对稳定的地质体。因此，笔者认为，本矿床成矿流体在成矿之前可能存在均一化过程。
4.2.3 （87Sr/86Sr）i组成及成矿物质来源
87Sr/86Sr值是判断成岩成矿物质来源的重要指标，在矿床地质研究中，常利用其示踪成矿物质来源。但应用锶同位素组成示踪成矿物质来源时，需要对潜在源区的锶同位素组成进行成矿年龄校正。噜鲁铅锌矿床成矿年龄为202.8±1.4Ma，且近年来，川滇黔成矿域成矿年代学研究取得了重要进展，成矿作用时间集中于约200Ma[28]。经年龄（200Ma）校正后，印支运动产生的基性玄武岩在成矿时的87Sr/86Sr值为0.7039~0.7078[38]；扬子地块西南缘各时代沉积的碳酸盐岩在成矿时的87Sr/86Sr值为0.7073~0.7111[23, 41-42]；昆阳群基底地层在成矿时的87Sr/86Sr值为0.7243~0.7288 [23]。地幔的87Sr/86Sr值为0.704±0.002[26]。而噜鲁铅锌矿床硫化物在成矿时的87Sr/86Sr值为0.7112~0.7115。通过分析，矿床中硫化物较地幔和峨眉山玄武岩有更多放射性成因锶，成矿物质可能来自于基底地层，或成矿流体流经这些富放射性成因岩体。
4.3 成矿机制
噜鲁铅锌矿床成矿时代主要发生在印支晚期，成矿作用和华南板块与印支板块拼合造山作用的时限相吻合。硫化物和硫酸盐矿物的硫同位素分析结果说明，硫主要为海相硫酸盐热化学还原的产物；硫化物铅同位素组成和Rb-Sr同位素组成证实，成矿物质主要来自基底岩石或成矿流体流经这些富放射性成因岩石。因此，噜鲁铅锌矿床成矿流体中组分来源不同，主要来自于基底岩石，或成矿流体流经这些富放射性成因岩体。基底岩石为噜鲁铅锌矿床提供了成矿流体及重要的成矿物质。噜鲁铅锌矿床围岩和黄铁矿稀土元素揭示，噜鲁铅锌矿床赋矿地层可能为成矿提供了成矿流体，成矿流体可能为富Cl流体，暗示成矿温度应该为中低温，成矿流体可能遭受了外来热液的混入（另文探讨）。
综合上述认识，基底岩石是矿床的成矿流体与成矿物质的重要来源。将该矿床成矿过程描述为，印支期构造旋回过程中，伴随着西面地区大面积玄武岩喷发与沉积建造，产生的热、气与地下热水，驱动来自基底岩石以氯化物络合物形式携带大量成矿金属的含矿热卤水溶液，含矿热卤水溶液沿途不断萃取、溶解高背景地层中的铅-锌-银等矿质，并还原不同时代地层沉积岩中海相硫酸盐的大气降水或层间水发生流体混合，伴随着物理、化学等成矿条件的改变，在有利的空间部分聚集成矿，金属硫化物沉淀形成具有工业价值的矿体。
5. 结论
噜鲁铅锌矿床成矿流体中不同组分来源不同。原生矿石硫化物硫同位素组成δ34S在9.72‰~ 22.44‰之间，多数为9‰~15‰。重晶石δ34S值介于26.08‰~29.34‰之间，平均为27.29‰。矿床中的硫同位素组成变化较大，δ34S值在6.33‰~9.75‰之间。噜鲁铅锌矿床主要是碳酸盐地层的海相硫酸盐热化学还原的产物。铅同位素组成特征表明，本区成矿物质可能由基底地层提供，区域上不同时代碳酸盐岩地层可能提供了部分成矿物质。由锶同位素可知，噜鲁铅锌矿床成矿物质可能来自基底地层，或成矿流体流经这些富放射性成因岩体。Rb-Sr同位素测年表明，噜鲁铅锌矿床成矿年龄为印支中晚期，印支运动对噜鲁铅锌成矿作用产生了重要影响。
硫、铅、铷-锶同位素数据研究成果表明，噜鲁铅锌矿床成矿流体中不同组分来源不同，但主要来自于基底地层，在印支运动强大的驱动力作用下，促使成矿元素活化-混合-迁移-聚集成矿。
致谢: 感谢云南省昭通气象局提供气象数据，感谢盐津县提供地质灾害记录。审稿专家提出的建议和意见使得本文得以完善。 -
![]()
图 2 庙坝滑坡剖面示意图
T1y—三叠纪永宁镇组；T1f—三叠纪飞仙关组
Figure 2. The schematic diagram of Miaoba landslide
![]()
图 4 降雨型滑坡破坏前单日降雨量最小值及发生次数
Figure 4. The minimum of daily rainfall and the occurrences before landslide damage
表 1 三种常用降雨概率分布函数及统计参数
Table 1 Three commonly used rainfall probability distribution function and statistical parameter
概率曲线及统计参数 公式 参数 P-Ⅲ型 P(X≥x)=βαΓ(α)+∞∫x(x−a0)α−1e−β(x−a0)dx a0=E(x)(1−2CVCS) α=4C2S β=2E(x)CVCS 指数分布 P(X≥x)=+∞∫xf(x)dx=+∞∫xαe−α(x−β)dx=e−α(x−β) \alpha = \frac{1}{b} = \frac{1}{{E(x)}}, {\rm{ \mathit{ β} = (}}{{\rm{X}}_{\rm{i}}}{{\rm{)}}_{\min }} 耿贝尔分布 P(X \ge x) = 1 - {e^{[ - {e^{ - a(x - b)}}]}} a = \frac{{1.2825}}{\sigma }, b = E(x) - 0.45005\sigma 统计参数 E(x) = x = \frac{1}{n}, \sum\limits_{i = 1}^n {{x_i}, {K_i} = \frac{{{x_i}}}{x}} {C_v} = \sqrt {\frac{1}{{n - 1}}\sum\limits_{i = 1}^n {{{({K_i} - 1)}^2}} } , {C_S} = \frac{{\sum\limits_{i = 1}^n {{{({K_i} - 1)}^3}} }}{{\left( {n - 3} \right)C_v^3}}, \sigma = \sqrt {\frac{{\sum\limits_{i = 1}^n {{{({x_i} - x)}^2}} }}{{n - 1}}} 
下载: 导出CSV
表 2 8月20—25日降雨概率分布曲线拟合参数结果
Table 2 Fitting parameter results of rainfall probability distribution curve from August 20th to 25th
分布类型及参数 8月20日 8月21日 8月22日 8月23日 8月24日 8月25日 皮尔逊Ⅲ型分布 E(x) 6.3 3.6 9.4 8.6 7.8 7.6 Cv 0.6 0.6 0.6 0.6 0.6 0.6 Cs 1.2 1.2 1.2 1.2 1.2 1.2 s 12.6 8.9 15.8 17.9 19.4 10.8 指数分布 a 0 0 0 0 0 0 b 6.4 3.6 9.5 8.5 7.9 7.7 s 10.8 7.8 7.8 15.2 9.9 8.3 耿贝尔分布 a 0.08366 0.12396 0.06349 0.05908 0.07936 0.0885 b -0.4682 -1.04356 0.44242 -1.00348 0.65429 1.1432 s 9.9 7.1 11 13.3 8.4 6.9 统计参数 E(Si) 6.3 3.6 9.6 8.6 7.8 7.6 D(Si) 235.0 107.0 408.0 471.2 261.2 210.2 Max(Si) 97.2 60.7 103.8 103.9 84.0 52.4 Min(Si) 0.0 0.0 0.0 0.0 0.0 0.0
下载: 导出CSV
表 3 8月20—25日降雨概率分布曲线及超过阈值概率
Table 3 The probability distribution curve of rainfall and the threshold exceeding the probability from August 20th to 25th
日期 超过概率分布公式 P(X≥x=29.7) 8月20日 P(X \ge x) = 1 - {e^{[ - {e^{ - 0.08366(x + 0.46820)}}]}} 0.07702 8月21日 P(X \ge x) = 1 - {e^{[ - {e^{ - 0.12396(x + 1.04356)}}]}} 0.02188 8月22日 P(X \ge x) = 1 - {e^{ - 0.10526x}} 0.04388 8月23日 P(X \ge x) = 1 - {e^{[ - {e^{ - 0.05908(x + 1.00348)}}]}} 0.15041 8月24日 P(X \ge x) = 1 - {e^{[ - {e^{ - 0.07936(x + 0.65429)}}]}} 0.09494 8月25日 P(X \ge x) = 1 - {e^{[ - {e^{ - 0.08850(x + 1014320)}}]}} 0.07677
下载: 导出CSV
-
王涛, 吴树仁, 石菊松.国际滑坡风险评估与管理指南研究综述[J].地质通报, 2009, 28(8):1006-1019. doi: 10.3969/j.issn.1671-2552.2009.08.002 石菊松, 石玲, 吴树仁, 等.滑坡风险评估实践中的难点与对策[J].地质通报, 2009, 28(8):1020-1030. doi: 10.3969/j.issn.1671-2552.2009.08.003 吴树仁, 石菊松, 张春山, 等.滑坡风险评估理论与技术[M].北京:科学出版社, 2012. 李媛, 孟晖, 董颖, 等.中国地质灾害类型及其特征——基于全国县市地质灾害调查成果分析[J].中国地质灾害与防治学报, 2004, 15(2):29-34. doi: 10.3969/j.issn.1003-8035.2004.02.005 陈洪凯, 魏来, 谭玲.降雨型滑坡经验性降雨阈值研究综述[J].重庆交通大学学报(自然科学版), 2012, 31(5):990-996. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=cqjtxyxb201205018 谢剑明, 刘礼领, 殷坤龙, 等.浙江省滑坡灾害预警预报的降雨阀值研究[J].地质科技情报, 2004, 22(4):101-105. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=dzkjqb200304019 陈洪凯, 魏来, 谭玲.降雨型滑坡经验性降雨阈值研究综述[J].重庆交通大学学报(自然科学版), 2012, 31(5):990-996. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=cqjtxyxb201205018 Maitir M.Assessing the triggering rainfall-induced landslip events in the ShivKhola Watershed of DARJILING Himalaya, west Bengal[J].European Journal of Geography, 2013, 4(3):21-37. https://www.researchgate.net/publication/264903246_ASSESSING_THE_TRIGGERING_RAINFALL-INDUCED_LANDSLIP_EVENTS_IN_THE_SHIVKHOLAWATERSHED_OF_DARJILING_HIMALAYA_WEST_BENGAL
Galeandro A, Doglioni A, Simeone V, et al.Analysis of infiltration processes into fractured and swelling soils as triggering factors of landslides[J].Environmental earth sciences, 2014, 71(6):2911-2923. doi: 10.1007/s12665-013-2666-7
Capparelli G, Versace P.Analysis of landslide triggering conditions in the Sarno area using a physically based model[J].Hydrology and Earth System Sciences, 2014, 18(8):3225-3237. doi: 10.5194/hess-18-3225-2014
Chien Yuan C, Tien Chien C, Fan Chieh Y, et al.Analysis of time-varying rainfall infiltration induced landslide[J].Environmental Geology, 2005, 48(4/5):466-479. doi: 10.1007/s00254-005-1289-z
Muntohar A S, Liao H J.Rainfall infiltration:infinite slope model for landslides triggering by rainstorm[J].Natural hazards, 2010, 54(3):967-984. doi: 10.1007/s11069-010-9518-5
魏来.降雨诱发滑坡预测模型研究[D].重庆交通大学硕士学位论文, 2013. 高华喜, 殷坤龙.降雨与滑坡灾害相关性分析及预警预报阀值之探讨[J].岩土力学, 2007, 28(5):1055-1060. doi: 10.3969/j.issn.1000-7598.2007.05.039 张珍, 李世海, 马力.重庆地区滑坡与降雨关系的概率分析[J].岩石力学与工程学报, 2005, 24(17):3185-3191. doi: 10.3321/j.issn:1000-6915.2005.17.029 陈剑, 杨志法, 李晓.三峡库区滑坡发生概率与降水条件的关系[J].岩石力学与工程学报, 2005, 24(17):3052-3056. doi: 10.3321/j.issn:1000-6915.2005.17.008 陈正洪, 孟斌.湖北省降雨型滑坡泥石流及其降雨因子的时空分布, 相关性浅析[J].岩土力学, 1995, 16(3):62-69. http://www.cnki.com.cn/Article/CJFDTotal-YTLX503.009.htm 陈丽霞, 殷坤龙, 刘长春.降雨重现期及其用于滑坡概率分析的探讨[J].工程地质学报, 2012, 20(5):745-750. doi: 10.3969/j.issn.1004-9665.2012.05.013 陈丽霞, 殷坤龙, 刘礼领, 等.江西省滑坡与降雨的关系研究[J].岩土力学, 2008, 29(4):1114-1120. doi: 10.3969/j.issn.1000-7598.2008.04.049 张桂荣, 殷坤龙, 陈丽霞, 等.滑坡预测预报的地质-气象耦合模型研究[J].水利水电技术, 2005, 36(3):5-18. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=slsdjs200503005 殷坤龙, 陈丽霞, 张桂荣.区域滑坡灾害预测预警与风险评价[J].地学前缘, 2008, 14(6):85-97. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=dxqy200706011 杜娟.单体滑坡灾害风险评价研究[D].中国地质大学(北京)博士学位论文, 2012. 黄振平.水文统计学[M].南京:河海大学出版社, 2003. Caine N.The rainfall intensity-duration control of shallow landslides and debris flows[J].Geografiska Annaler:Series A, Physical Geography, 1980, 62(1/2):23-27. doi: 10.2307/520449
Guzzetti F, Peruccacci S, Rossi M, et al.The rainfall intensity-duration control of shallow landslides and debris flows:an update[J].Landslides, 2008, 5(1):3-17. doi: 10.1007/s10346-007-0112-1
Huang Q X, Xu X T, Kulatilake P H S W, et al.Formation mechanism of a rainfall triggered complex landslide in southwest China[J].Journal of Mountain Science, 2020, 17(5):1128-1142. doi: 10.1007/s11629-019-5736-9
周伟杰, 徐卫亚, 王如宾, 等.暴雨及久雨作用下东岭信滑坡堆积体的渗流特性及稳定性分析[J].三峡大学学报(自然科学版), 2020, 42(2):28-33. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=whsldldxxb-yc202002006 张勇, 温智, 程英建.四川巴中市滑坡灾害与降雨雨型关系探讨[J].水文地质工程地质, 2020, 47(2):178-182. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=swdzgcdz202002023 甘建军, 储小东.中低山台风暴雨型泥石流形成机制和动力特征——以江西德安杜家沟泥石流为例[J].灾害学, 2020, 35(1):150-155. doi: 10.3969/j.issn.1000-811X.2020.01.027 向小龙, 孙炜锋, 李国伟, 等.云南盐津地区地质灾害影响因素分析[J].地质力学学报, 2015, 21(1):97-105. doi: 10.3969/j.issn.1006-6616.2015.01.011 -
期刊类型引用(9)
1. 刘玲，张创，孙明. 志丹油田纸坊北油区三叠系延长组长6—长9储层致密史与油藏成藏史研究. 地质与资源. 2023(03): 327-334 . 
百度学术
2. 韩晓洁，范昌育，高潮，张丽霞，尹锦涛，王成达，王宁. 构造抬升区欠压实超压恢复方法——以鄂尔多斯盆地下寺湾地区延长组为例. 天然气地球科学. 2023(07): 1163-1172 . 百度学术
3. 徐泽阳，刘震，赵靖舟，党佳城，李军，唐子怡. 基于测井响应特征的烃源岩古超压成因分析——以鄂尔多斯盆地三叠系延长组7段烃源岩为例. 天然气地球科学. 2023(11): 1950-1960 . 百度学术
4. 王苗，吴柏林，李艳青，刘池阳，郝欣，刘明义，张婉莹，李琪，姚璐航，张效瑞. 鄂尔多斯盆地深部富铀烃源岩提供铀源可能性的实验研究. 地球科学. 2022(01): 224-239 . 百度学术
5. 邹敏，云金表，王濡岳，于岚，伍岳，吉圆圆，成立，何维领，赵刚. 镇泾地区中侏罗统延安组低幅度圈闭成因类型与评价. 断块油气田. 2022(03): 313-318+330 . 百度学术
6. 胡艳飞，孔庆莹. 鄂尔多斯盆地西南部长8油层储层主控因素及分布规律. 吉林大学学报(地球科学版). 2022(04): 1078-1090 . 百度学术
7. 马立元，邱桂强，胡才志，徐士林，陈纯芳，李松. 鄂尔多斯盆地红河油田延长组致密砂岩成藏机制分析. 沉积学报. 2022(06): 1762-1773 . 百度学术
8. 郑登艳，王震亮，王联国，黄昊. 低幅度构造特征及其对油气成藏的控制作用：以鄂尔多斯盆地彭阳地区延安组为例. 现代地质. 2021(04): 1114-1123 . 百度学术
9. 庞宏，吴松，胡英杰，王海朋，丁旭光，刘兴周，惠沙沙，陈昌，郭军敏. 辽河拗陷牛心坨沙四段致密砂岩油形成主控因素及有利区预测. 石油科学通报. 2020(04): 467-482 . 百度学术
其他类型引用(3)
计量
- 文章访问数: 3242
- HTML全文浏览量: 285
- PDF下载量: 2620
- 被引次数: 12
