Geothermal field in Tongzhou District of Beijing and its influencing factors
-
摘要:
在收集整理大量地热井温度资料的基础上，系统分析了北京通州地区现今地温场特征，划分了地热异常区，并进一步估算了1000~3000 m埋深处的地层温度。研究表明，研究区现今地温梯度为1.34~3.12℃/100 m，平均值为1.88℃/100 m。大地热流值分布范围为47.1~75.9 mW/m2，平均值为63.2 mW/m2，其热流状态较好。在垂向上，研究区1000 m埋深处整体温度均超30℃，2000 m时高温异常区温度超60℃，3000 m时高温异常区温度超80℃。根据深部地层温度估算，在研究区的南部及东南部具有2处高温异常区，地温场分布形态与区域内主要断裂走向基本一致，主要受区域地质构造控制，岩石性质为地温场分布的重要影响因素。
Abstract:Based on a large number of geothermal well temperature data, the distribution of the current geothermal field in Tongzhou District of Beijing is systematically analyzed, the geothermal anomaly area is divided, and the formation temperature at the depth of 1000~3000 m is further estimated. The results show that the current geothermal gradient ranges from 1.34 ℃/100 m to 3.12 ℃/100 m with an average of 1.88 ℃/100 m. Terrestrial heat flow values range from 47.1 mW/m2 to 75.9 mW/m2 with an average of 63.2 mW/m2. Vertically, the estimated temperature in the study area is over 30 ℃ at depth of 1000 m, the temperature of the high temperature abnormal area is over 60 ℃ at 2000 m, and over 80 ℃ at 3000 m. According to the estimation of deep formation temperature, there exist two high temperature anomaly areas in its south and southeast. Controlled by regional geological structure, the distribution pattern of geothermal field is basically consistent with the strike of regional main faults. In addition, the lithology is also one of the important influencing factors.
-
川-滇-黔铅锌银多金属成矿区是中国主要的铅锌银矿产基地之一，共发现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同位素测年表明，噜鲁铅锌矿床成矿年龄为印支中晚期，印支运动对噜鲁铅锌成矿作用产生了重要影响。
硫、铅、铷-锶同位素数据研究成果表明，噜鲁铅锌矿床成矿流体中不同组分来源不同，但主要来自于基底地层，在印支运动强大的驱动力作用下，促使成矿元素活化-混合-迁移-聚集成矿。
-
![]()
图 1 研究区基岩地质构造图
F1—南苑-通县断裂；F2—燕郊断裂；F3—夏垫断裂；F4—牛堡屯断裂；F5—张家湾断裂
Figure 1. Geological map showing structure of the survey area
![]()
图 3 通州地区1000 m深度温度等值线图
Figure 3. Temperature contour at depth of 1000 m in Tongzhou District
![]()
图 4 通州地区2000 m深度温度等值线图
Figure 4. Temperature contour at depth of 2000 m in Tongzhou District
![]()
图 5 通州地区3000 m深度温度等值线图
Figure 5. Temperature contour map at depth of 3000 m in Tongzhou District
表 1 地热井地温梯度计算结果
Table 1 Geothermal gradients of geothermal well
地热井序号 成井深度/m 盖层深度/m 测温井段/m T/℃ G/(℃·100 m-1) G′/(℃·100 m-1) 1 811 0~445 100~800 29 1.95 1.81 2 642 0~473 100~600 28 2.46 3.16 3 1827 203~589 340~1800 43 1.64 2.18 4 2500 864~1418 1000~2417 62.5 2.03 2.16 5 2679 515~1065 1500~2670 56.8 1.62 - 6 2213 583~1108 1500~2210 53.2 1.80 - 7 2203 464~1049 1100~2200 53.8 1.83 - 8 2699 277~729 300~2670 57.1 1.63 1.84 9 2509 680~1225 100~2460 53.1 1.61 2.03 10 2480 439~979 1500~2400 55.34 1.74 4.09 11 2405 470~792 100~2350 45.2 1.34 1.10 12 2400 297~677 100~2400 51.5 1.58 1.24 13 2505 490~1073 100~2500 58.2 1.79 2.61 14 2106 926.5~1005 100~2100 49.89 1.73 2.32 15 2468 302~712 100~2400 57.1 1.82 3.38 16 3203 2404~2740 50~3200 84 2.21 1.66 17 3589 2817~3180 300~3425 89.5 2.22 2.17 18 2800 570~1086 50~2800 63.5 1.79 0.78 19 3001 255~618 100~3000 62.983 1.65 1.83 20 2800 0~478 100~2800 53.98 1.44 1.82 21 2800 - 300~2800 56.74 1.54 - 22 3000 - 286~3000 58.53 1.50 - 23 2204 504~1105 50~2190 53.7 1.84 2.91 24 2106 610~980 100~2100 45.94 1.54 0.46 25 2500 592~842 100~2500 51.4 1.51 0.76 26 3600 - 50~3600 116.4 2.87 - 27 3000 - 50~3000 97.61 2.82 - 28 3332 - 50~3330 116.9 3.12 - 注：T表示井底温度；G表示地热井地温梯度；G′表示盖层段地温梯度 
下载: 导出CSV
-
王均, 黄尚瑶, 黄歌山. 中国地温分布的基本特征[M]. 北京: 地震出版社, 1990, 26: 80-100. 蔺文静, 刘志明, 王婉丽, 等. 中国地热资源及其潜力评估[J]. 中国地质, 2013, 40(1): 312-321. doi: 10.3969/j.issn.1000-3657.2013.01.021 卫万顺, 郑桂森, 栾英波. 北京平原区浅层地温场特征及其影响因素研究[J]. 中国地质, 2010, 37(6): 1733-1739. doi: 10.3969/j.issn.1000-3657.2010.06.019 陈晓冬, 金旭, 管彦武, 等. 长春地区地表温度日变、年变对地温测量的影响[J]. 地球物理学进展, 2006, (3): 1008-1011. doi: 10.3969/j.issn.1004-2903.2006.03.047 栾英波, 郑桂森, 卫万顺. 北京平原区粉质粘土热导率影响因素实验研究[J]. 中国地质, 2013, 40(3): 981-988. doi: 10.3969/j.issn.1000-3657.2013.03.028 于仲坤, 赵洪, 刁慧, 等. 东海陆架盆地丽水凹陷热演化模拟及现今地温场特征[J]. 海洋地质与第四纪地质, 2020, 40(2): 124-134. https://www.cnki.com.cn/Article/CJFDTOTAL-HYDZ202002013.htm 周阳, 穆根胥, 张卉, 等. 关中盆地地温场划分及其地质影响因素[J]. 中国地质, 2017, 44(5): 1017-1026. https://www.cnki.com.cn/Article/CJFDTOTAL-DIZI201705014.htm 柏崇伟, 孙志国, 金洪文, 等. 陇东盆地地热系统地温场特征及平凉地热勘查评价[J]. 甘肃地质, 2017, 26(1): 74-80. https://www.cnki.com.cn/Article/CJFDTOTAL-GSDZ201701011.htm 李文, 孔祥军, 袁利娟, 等. 北京通州地区地热井酸化压裂增灌试验研究[J]. 城市地质, 2019, 14(4): 43-48. doi: 10.3969/j.issn.1007-1903.2019.04.008 高剑, 张进平, 孔祥军, 等. 北京通州地区深部地热井间的影响分析[J]. 城市地质, 2017, 12(4): 72-75. doi: 10.3969/j.issn.1007-1903.2017.04.014 蔡向民, 栾英波, 郭高轩, 等. 北京平原地区地质系统[J]. 城市地质, 2009, 4(3): 6-12. doi: 10.3969/j.issn.1007-1903.2009.03.002 宾德智, 刘久荣, 王小玲. 北京地热资源[C]//北京地热国际研讨会论文集, 2002: 169-177. 北京地质局. 北京市区域地质志[M]. 北京: 地质出版社, 1991: 213-260. 蔡向民, 郭高轩, 栾英波, 等. 北京山前平原区第四系三维结构调查方法研究[J]. 地质学报, 2009, 83(7): 1047-1057. doi: 10.3321/j.issn:0001-5717.2009.07.012 陈墨香. 华北地热[M]. 北京: 科学出版社, 1988: 1-100. 姜光政, 高堋, 饶松, 等. 中国大陆地区大地热流数据汇编(第四版)[J]. 地球物理学报, 2016, 59(8): 2892-2910. https://www.cnki.com.cn/Article/CJFDTOTAL-DQWX201608015.htm 徐世光, 郭远生. 地热学基础[M]. 北京: 科技出版社, 2009: 15. 雷晓东, 胡圣标, 李娟, 等. 北京平原区西北部大地热流与深部地温分布特征[J]. 地球物理学报, 2018, 61(9): 3735-3748. https://www.cnki.com.cn/Article/CJFDTOTAL-DQWX201809021.htm 汪集旸, 胡圣标, 庞忠和, 等. 中国大陆干热岩地热资源潜力评估[J]. 科技导报, 2012, 30(32): 25-31. doi: 10.3981/j.issn.1000-7857.2012.32.002 柯柏林, 林天懿, 李文, 等. 北京西山谷积山背斜地热系统成因模式及远景区预测[J]. 地质通报, 2019, 38(8): 1378-1385. https://www.cnki.com.cn/Article/CJFDTOTAL-ZQYD201908015.htm 阎敦实, 于英太. 京津冀油区地热资源评价与利用[M]. 武汉: 中国地质大学出版社, 2000: 19. 胡玉禄, 胡红文, 邱希青, 等. 华北地温场垂向变化[C]//山东省科学技术协会. 地质与可持续发展——华东六省一市地学科技论坛文集. 山东省科学技术协会, 2003: 3. 邱楠生. 中国西北部盆地岩石热导率和生热率特征[J]. 地质科学, 2002, (2): 196-206. https://www.cnki.com.cn/Article/CJFDTOTAL-DZKX200202006.htm -
期刊类型引用(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)
计量
- 文章访问数: 2876
- HTML全文浏览量: 325
- PDF下载量: 2630
- 被引次数: 12
