Research on Seepage Barriers of Massive and Superimposed Thick Oil Layers in the Low Oil Formation of No. 9 Block in Tahe Oilfield
-
摘要:
塔河油田9区三叠系阿克库勒组下油组储层为典型的块状叠置辫状河三角洲储层，夹层分布及连通性极其复杂，封闭断层和泥质屏障、钙质胶结帯等不同级次渗流屏障控制着油藏内部油水运动。为深化储层空间非均质表征，在储层构型分析基础上，笔者采用层次分析方法，动静态相结合，对渗流屏障类型及其级次进行划分，分析泥质屏障、钙质胶结帯形成机理及其受控因素，建立块状厚油层渗流屏障分布地质模型。研究表明：塔河油田9区下油组主要渗流屏障为封闭断层屏障和层间、层内发育的泥质及钙质隔夹层；依据渗流阻挡作用分为4个级次，东北向封闭性断层及复合河道顶部6级界面是1级渗流屏障，3～5级构型界面控制了2～4级泥质、钙质渗流屏障的发育，建立完全不遮挡型、部分遮挡型及完全遮挡型等3种渗流屏障模式。通过分级次定量表征渗流屏障的空间分布，深化储层空间非均质性研究，为基于流动单元的精细地质建模奠定基础。
Abstract:The low oil formation reservoir of No. 9 block in the Tahe Oilfield is a superimposedmassive braided river delta reservoir, moreover, distribution and connectivity of the interlayers as well as complexity produced the sealing faults, argillaceous barriers, calcareous cemented zones and other different–level seepage barriers to control the movement of oil and water within the reservoir is to deepen the heterogeneity of the reservoir space; and the analytic hierarchy process is adopted based on the analysis of the reservoir configuration, and combined with dynamic and static, to clarify the seepage barrier types and levels;the formation mechanism of calcareous cement zone and its controlling factors is analyzed by researcher, so that establishes a geological model for the distribution of seepage barriers in massive thick oil layers; Researches show that the main seepage barriers of the low oil formation in No. 9 Block of Tahe Oilfield are sealing faults barriers, mud barriers, and calcium cemented zones, what is more, it can be divided into four levels according to the seepage blocking effect, and northeast sealing faults and the 6–level interface at the top of the composite channel are the first–level seepage barrier, the 3～5 level configuration interface controls the 2～4 level mud and calcareous seepage barriers; to establish three types of seepage barrier modes composed by completely uncover, partial cover, and completely cover. The spatial distribution of seepage barriers is characterized by hierarchical and sub–quantitative methods, which deepening the study of reservoir spatial heterogeneity, and laying the foundation for fine geological modeling based on flow units.
-
青海东昆仑地区是中国重要的造山带之一（吴树宽等，2023），祁漫塔格成矿带位于该造山带西段（图1a）。近年来，随着该地区众多岩浆热液型多金属矿床（点）的发现，其已逐渐成为中国西部重要的、最具找矿潜力的成矿远景勘查区之一，引起了地质科研者们浓厚的研究兴趣，推动了该地区成矿规律研究工作的不断进行（丰成友等，2011；张爱奎，2012；高永宝等，2014；钟世华等，2017a，2017b）。
![]()
Figure 1. (a) Tectonic sketch map of east Kunlun orogenic beltand geology and (b) mineral deposit distribution in the Qimantage area牛苦头多金属矿床位于祁漫塔格成矿带中段野马泉–开木棋河东侧（图1b），前人对牛苦头矿区开展了一系列研究工作，认为牛苦头矿区矽卡岩底部隐伏的花岗闪长岩为成矿岩体（蒋成伍，2013；贾建团，2013），其成矿时代与区域大多数铅锌多金属矿床（卡尔却卡、虎头崖、肯德可克、野马泉、尕林格等铁多金属矿床）成矿时代一致，为晚三叠纪。随着牛苦头矿区勘查工作的进展以及研究工作的不断深入，前人逐渐意识到，华力西期（晚泥盆世）可能为矿区主要成矿岩体（李加多等，2019；王新雨等，2020，2021）。然而证据不够充分，尤其是矿石矿物年代学成矿研究工作较为空白。此外，矿区成矿物质来源问题一直没有很好的解决，前人研究仅仅从成矿岩体以及金属中S、Pb的物质来源着手（王新雨等，2021），金属成矿物质来源研究较为匮乏。
金属矿物Re–Os测年是近年来矽卡岩多金属矿床成矿时代研究较为较为有效的一种测年方法，其中，最为有效的属于辉钼矿Re–Os测年（李超等，2016）。然而，对于辉钼矿含量较低的多金属矿来说，辉钼矿测年具有较大的难度和局限性。而黄铁矿、黄铜矿中Re含量低，并且具有一定的初始Os，其Re–Os等时线年龄不仅能够限定金属矿物成矿时代，而且能够示踪金属矿物成矿物质来源（李超等，2009）。
笔者拟从成矿岩体锆石U–Pb年代学、矿石矿物Re–Os年代学入手，探讨成岩时代与成矿时代的耦合关系，为矿区铅锌矿床成矿时代提供新的有力证据。尝试采用黄铁矿的同位素初始比值对其成矿物质来源进行示踪，以解决牛苦头多金属矿床的成矿时代以及成矿来源等成矿相关问题，为研究野马泉–牛苦头–四角羊地区矽卡岩型铁锌多金属矿床成矿地质特征、矿床成因以及找矿预测工作提供相关理论依据。
1. 地质背景
1.1 区域地质
东昆仑造山带西侧的祁漫塔格成矿带位于青藏高原中北部，成矿条件较为优越（毛景文等，2012；Yu et al.，2017；刘渭等，2021）（图1a、图1b）。成矿带地区出露地层复杂，主要包括以下地层与岩性：新元古界金水口群主要为一套夹片岩、大理岩、片麻岩的混合岩；上奥陶统滩间山群岩组主要为一套夹火山岩组的碎屑岩地层；下石炭统大干沟组与缔敖苏组为一套沉积碳酸盐地层，岩性包括生物碎屑灰岩、大理岩、白云岩、鲕粒状灰岩和砂岩；下—中二叠统打柴沟组为一套灰岩夹炭质条带岩的岩石组合；第四系包括风积、冲积物和亚砂土。该地区构造活动发育较为复杂，构造受昆中、那陵郭勒断裂影响，产状主要呈NWW、NW向，局部近EW向，也常常发育近EW向的褶皱构造。区内岩浆岩被第四系覆盖严重，侵入岩发育强烈，时代分为华力西期和印支期2组，受NWW和N向构造断裂控制明显。
1.2 矿床地质
牛苦头矿区出露地层为上奥陶统滩间山群岩组与第四系（图2）。滩间山群在本矿区表现为一套浅海相沉积地层（碳酸盐岩、碎屑岩），岩性表现为灰岩、大理岩、条带状灰岩、碳质灰岩、粉砂岩、泥质粉砂岩。矿区内断裂分为NWW和NE向，以NWW向为主，为昆北断裂的次级分支断裂，但NE向也同样明显，NE向构造代表主要为牛苦头沟和四角羊沟断裂。矿区内岩浆岩种类复杂，包括辉长闪长岩、石英闪长岩、二长花岗岩、正长花岗岩、似斑状花岗岩以及富石英花岗岩类等，空间穿插关系复杂。
矿区M1矿段作为其主要开采地段，其矿体类型主要为层状或不规则脉状矿体，矿石成分主要为磁黄铁矿、闪锌矿、磁铁矿、方铅矿、黄铜矿、黄铁矿，可见少量毒砂与赤铁矿。矿石构造主要呈块状、浸染状、细脉浸染状、条带状以及班杂状。M4矿段矿体产状、矿石结构构造与M1基本类似，但其矿石成分与M1矿段有所区别，矿石中磁铁矿含量较大，其次为磁黄铁矿、闪锌矿、黄铁矿，方铅矿含量较少。
矿区主要矿体顺层产出，受地层和硅钙面控制明显（图3），局部受导矿构造影响，可切穿地层。脉石矿物与区域上其他矽卡岩矿床有所区别，局部矽卡岩Mn含量较高，与铅锌矿化关系密切。脉石矿物包括（锰）阳起石、（锰）黑柱石、（锰）钙铁辉石、石榴子石、硅辉石、石英、方解石、绿帘石、绿泥石等（图4）。花岗岩体位于矿体与矽卡岩之下，为二长花岗岩–花岗闪长岩组合，蚀变为绿泥石化、钾化及少量辉石化，局部花岗岩中可见星点状闪锌矿化与黄铁矿化。
![]()
图 3 牛苦头M1磁异常区10号地质勘探线剖面图Figure 3. Geological profile of No.10 exploration line n in M1 magnetic anomaly area of Niukutou deposit![]()
图 4 岩石、矿石手标本图a. 块状磁黄铁矿矿石；b. 斑杂状闪锌矿化，产于锰钙辉石矽卡岩中；c. 块状闪锌矿、磁黄铁矿石；d. 磁黄铁矿化石榴子石矽卡岩；e. 闪锌矿化含方解石黑柱石辉石矽卡岩；f. 闪锌矿化锰钙辉石矽卡岩；Q. 石英；Cal. 方解石；Ilv. 黑柱石； Jo. 锰钙辉石；Px. 辉石；Sp. 闪锌矿；Po. 磁黄铁矿；Grt. 石榴子石Figure 4. Specimens of rocks and ores2. 样品采集与测试方法
研究样品采自M1矿段10线钻孔底部，岩性分别为花岗闪长岩和二长花岗岩（图5a、图5b），取样位置见图3。
![]()
图 5 牛苦头矿区花岗岩类手标本及镜下照片a. M1矿段钻孔ZK1007底部强蚀变花岗闪长岩手标本；b. M1矿段钻孔ZK1009底部蚀变二长花岗岩发育绢云母化；c. M1矿段钻孔底部强蚀变花岗闪长岩中发育的绢云母化长石（正交偏光）；d. M1矿段钻孔底部蚀变二长花岗岩中发育的绢云母化长石与环带状中长石（正交偏光）Figure 5. Specimen and microscope photo of the granites from the Niukutou depositM1矿段钻孔底部花岗闪长岩，经鉴定岩性为中粒钾化黑云母花岗闪长岩（图5a），石英含量约为25%，自形–半自形结构，粒径为0.6～1.5 mm。斜长石，自形结构，含量约为55%，粒径为1～1.5 mm，斜长石表面多发生泥化、绢云母化。钾长石含量约为10%～15%，粒径为0.5～0.8 mm，多为正长石，局部表面发育泥化。黑云母含量约为3%～5%，片状，局部为泥化。花岗闪长岩中普遍发育钾化、硅化、绢云母化，表明M1底部岩体发生较强的热液蚀变（图5c）。
M1矿段钻孔底部二长花岗岩一般为中细粒结构（图5b），块状构造，其主要由石英（35%）、正长石（35%）、斜长石（25%）和黑云母（5%）组成。斜长石中可发现具有环带结构的中长石。蚀变主要为长石的钾化和绢云母化，同时普遍发育浸染状黄铁矿化（图5d）。
选取M1矿段10线钻孔ZK1007、ZK1009底部的2件花岗闪长岩（NZC01、NZC13）进行锆石U–Pb定年分析（图3）。Re–Os测年样品则选择M1矿段ZK1004钻孔底部黄铁矿闪锌矿矿石样品6件作为测试对象（图3）。
3. 测试方法
3.1 锆石LA–ICP–MS U–Pb测年
花岗岩锆石测年在合肥工业大学矿物原位分析实验室完成，分析仪器为Agilent7900电感耦合等离子质谱仪（美国），与之配套的为美国生产的Analyte气态准分子激光剥蚀系统。激光以He为载气，统一剥蚀半径为30 um。标样采用91500、PLE为外标的锆石进行校正，每测10个样品，用2个91500标样和1个PLE标样进行校正。随时观察仪器的信号，以保证数据的精确有效。数据处理及普通Pb校正采用CP–MS–Da–Ta–Cal（Liu，2010；李艳广等，2023）和EXCEL宏程序ComPbCorr#3-17（Andersen，2002），年龄谐和图解使用Isoplot3.0获取（王新雨等，2021）。
3.2 黄铁矿Re–Os测年
用于黄铁矿 Re–Os 同位素测试的样品采自牛苦头M1矿段ZK1004钻孔底部块状矿石，共 6 件，其中块状矿石中黄铁矿与磁黄铁矿、闪锌矿共生，为同阶段硫化物。将包含黄铁矿的闪锌矿矿石进行粉碎，剔除杂物，保证实验的黄铁矿样品纯度高于99%。
在中国地质科学院国家测试中心进行了黄铁矿 Re–Os 同位素组成测定。对于黄铁矿含量低的样品，采用逐级剥谱法扣除O同位素的干扰。其中Re同位素质量分馏校正采用的普通Re的185Re/187Re值为0.59738。而Os同位素质量分馏校正采用192Os/188Os值作为内标（杜安道等，2009）。详细测定方法及流程见2009，李超等（2016）。
4. 分析结果
4.1 锆石U–Pb定年
花岗闪长岩和二长花岗岩采样位置见图3，其锆石U–Pb同位素组成见表1。在CL图像下，花岗闪长岩中（NZC01）锆石呈无色透明的长–短柱状（图6）。长约为135～265 μm，宽为45～75 μm。大多数锆石具岩浆震荡环带，局部可见锆石形态残缺不完整，可能为晚阶段热液作用所致。锆石中Th含量为63×10−6～393×10−6，U含量为183×10−6～1085×10−6，Th/U值为0.3～0.56（均大于0.3），属于典型岩浆锆石成因。20个锆石LA–ICP–MS加权平均年龄为（362.2±2.7）Ma（MSWD=3.9，n=20）（图6），均为华力西期晚泥盆纪。
表 1 牛苦头成矿花岗岩锆石U–Pb数据统计表Table 1. U–Pb isotopic compositions of Niukutou granitoids测点号 U Th Th/U 同位素比值 年龄（Ma） （10−6） 206Pb/238U 1σ 207Pb/235U 1σ 207Pb/206Pb 1σ 206Pb/238U 1σ 207Pb/235U 1σ NZC01.1 1085 393 0.36 0.05738 0.0006 0.42570 0.0063 0.05385 0.0007 360 2 365 3 NZC01.2 225 74 0.33 0.05758 0.0006 0.41379 0.0106 0.05215 0.0013 361 2 292 5 NZC01.3 266 81 0.30 0.05725 0.0006 0.43111 0.0086 0.05461 0.0010 359 2 396 4 NZC01.4 240 91 0.38 0.06005 0.0006 0.44513 0.0114 0.05378 0.0013 376 2 362 5 NZC01.5 201 71 0.35 0.05615 0.0006 0.41780 0.0102 0.05403 0.0013 352 2 372 5 NZC01.6 358 173 0.48 0.05931 0.0006 0.44755 0.0094 0.05470 0.0010 371 2 400 4 NZC01.7 195 70 0.36 0.05930 0.0006 0.44163 0.0107 0.05402 0.0012 371 2 372 5 NZC01.8 270 114 0.42 0.05856 0.0006 0.42962 0.0087 0.05327 0.0011 367 2 340 4 NZC01.9 183 63 0.34 0.05754 0.0006 0.42627 0.0110 0.05371 0.0013 361 2 359 5 NZC01.10 271 82 0.30 0.05786 0.0006 0.43239 0.0090 0.05411 0.0010 363 2 376 4 NZC01.11 198 102 0.52 0.05783 0.0006 0.42305 0.0107 0.05312 0.0013 362 2 334 5 NZC01.12 259 119 0.46 0.05650 0.0006 0.41202 0.0088 0.05296 0.0011 354 2 327 4 NZC01.13 308 136 0.44 0.05815 0.0008 0.42186 0.0112 0.05264 0.0013 364 3 313 5 NZC01.14 284 115 0.41 0.05696 0.0007 0.41795 0.0097 0.05337 0.0012 357 2 344 5 NZC01.15 241 102 0.42 0.05745 0.0006 0.42972 0.0091 0.05423 0.0011 360 2 381 4 NZC01.16 190 91 0.48 0.05779 0.0007 0.42768 0.0137 0.05392 0.0017 362 3 368 6 NZC01.17 317 177 0.56 0.05796 0.0007 0.43976 0.0098 0.05501 0.0011 363 2 413 4 NZC01.18 387 195 0.50 0.05840 0.0006 0.43639 0.0080 0.05421 0.0009 366 2 380 3 NZC01.19 300 124 0.41 0.05734 0.0006 0.41698 0.0080 0.05282 0.0009 359 2 321 4 NZC01.20 332 162 0.49 0.05755 0.0006 0.41550 0.0079 0.05234 0.0009 361 2 300 3 NZC13.1 201 80 0.40 0.05828 0.0006 0.42108 0.0105 0.05242 0.0012 365 4 357 9 NZC13.2 280 103 0.37 0.05685 0.0006 0.41206 0.0091 0.05264 0.0010 356 4 350 8 NZC13.3 185 68 0.37 0.05879 0.0007 0.44450 0.0129 0.05476 0.0014 368 4 373 11 NZC13.4 208 119 0.57 0.05796 0.0006 0.43336 0.0114 0.05419 0.0013 363 4 366 10 NZC13.5 228 88 0.39 0.05721 0.0006 0.42424 0.0117 0.05394 0.0014 359 4 359 10 NZC13.6 565 178 0.32 0.06000 0.0007 0.44218 0.0084 0.05355 0.0008 376 5 372 7 NZC13.7 193 66 0.34 0.05623 0.0006 0.41640 0.0119 0.05383 0.0015 353 4 353 10 NZC13.8 182 85 0.47 0.05767 0.0006 0.41970 0.0119 0.05307 0.0014 361 4 356 10 NZC13.9 244 120 0.49 0.05763 0.0006 0.42948 0.0095 0.05409 0.0011 361 4 363 8 NZC13.10 185 66 0.36 0.05949 0.0007 0.43490 0.0117 0.05297 0.0013 373 4 367 10 NZC13.11 131 48 0.37 0.05665 0.0006 0.41935 0.0135 0.05376 0.0016 355 4 356 11 NZC13.12 247 91 0.37 0.05775 0.0006 0.41730 0.0090 0.05248 0.0010 362 4 354 8 NZC13.13 251 106 0.42 0.05773 0.0006 0.42749 0.0090 0.05371 0.0010 362 4 361 8 NZC13.14 284 107 0.38 0.05662 0.0006 0.41479 0.0081 0.05314 0.0009 355 4 352 7 NZC13.15 167 56 0.34 0.05839 0.0006 0.43116 0.0113 0.05354 0.0013 366 4 364 10 ![]()
图 6 牛苦头矿区成矿花岗闪长岩、二长花岗岩锆石年龄谐和及加权平均年龄图Figure 6. Zircon U–Pb Concordia and weighted average of granitoids in Niukutou ore district牛苦头M1矿段二长花岗岩中（NZC013）锆石CL图像晶形较好，呈无色透明短柱状，长约为150～250 μm，宽约为80～100 μm，少数锆石局部出现残缺，可能为晚期岩浆热液所致。锆石Th含量为48×10−6～178×10−6，U含量为131×10−6～565×10−6，Th/U值为0.34～0.57（均大于0.3），为典型岩浆锆石成因（王新雨等，2021）。 测得15个锆石LA–ICP–MS加权平均年龄为（361.8±3.4）Ma（MSWD=2.5，n=15）（图6），均属于华力西期晚泥盆纪。
4.2 Re–Os地球化学
黄铁矿Re–Os数据分析结果见表2。矿石样品中黄铁矿与闪锌矿共生关系见图7a。牛苦头铅锌矿床黄铁矿中普Re含量为0.031×10−9～7.887×10−9，187Os含量为0.20×10−12～29.96×10−12，普Os含量较低，相对于187Os可忽略不计。黄铁矿187Re/188Os值为61.7～32860，属于低含量、高放射成因硫化物（LLHR），黄铁矿Re–Os等时线年龄为（359.2±6.3）Ma（图7b），个别数据的等时线年龄为（352±15）Ma，初始187OS/188Os值为0.13±0.24（表2）。个别数据的等实线年龄与M4钻孔底部成矿岩体锆石U–Pb年龄一致（353.0±3.6 Ma）（王新雨等，2021）。
表 2 牛苦头铅锌矿床黄铁矿Re–Os同位素数据统计表Table 2. Re–Os data of pyrite from Niukutou Pb–Zn skarn deposit样品号 Re （10−9） 普Os （10−9） 187Re （10−9） 187Os（10−9） 187Re/188Os 187Os/188Os 模式年龄（Ma） 测定值 不确
定度测定值 不确
定度测定值 不确
定度测定值 不确
定度测定值 不确
定度测定值 不确
定度测定值 不确
定度NKC073 7.887 0.058 0.0012 0.0000 4.957 0.037 0.02996 0.00023 32860 381 198.0 1.3 361.7 3.7 NKC072 0.363 0.003 0.0009 0.0000 0.228 0.002 0.00136 0.00001 1963 20 11.69 0.02 356.0 3.7 NZC-115 0.887 0.007 0.0066 0.0001 0.557 0.004 0.00335 0.00003 648.9 6.7 3.880 0.010 NKC074 0.214 0.002 0.0028 0.0000 0.135 0.001 0.00079 0.00001 376.2 4.5 2.222 0.014 NKC075 0.254 0.002 0.0014 0.0000 0.160 0.001 0.00098 0.00001 852.6 16.9 5.191 0.078 NKC076 0.031 0.000 0.0024 0.0000 0.019 0.000 0.00020 0.00001 61.7 1.4 0.611 0.021 ![]()
图 7 闪锌矿与黄铁矿镜下共生的显微照片（a）和牛苦头M1矿段黄铁矿Re–Os同位素等时线图解（b）Figure 7. (a) micrograph of sphalerite intergrowth with pyriteand and (b) Re–Os isotope isochron of pyrite in M1 ore block of Niukutou ore district5. 讨论
5.1 牛苦头铅锌矿床成矿时代
黄铁矿Re–Os同位素定年方法，是解决金属矿床成矿时代问题一种有效的技术手段。牛苦头M1磁异常区铅锌矿与闪锌矿共生的黄铁矿Re–Os加权平均值年龄为（359.2±6.3）Ma，代表了牛苦头铅锌矿床的成矿时代为晚泥盆世。该年龄与矿区M1和M4华力西晚期钻孔底部成矿岩体年龄（363～362 Ma）基本一致，代表了祁漫塔格地区晚泥盆世存在一期重要的矽卡岩型铅锌成矿作用。综上所述，华力西期岩浆作用是祁漫塔格矿区一期重要的岩浆侵入作用，对应于晚泥盆世始特提斯洋后碰撞伸展作用。该期岩浆岩在整个牛苦头矿区普遍都存在，尤其在M3和M5地区。因此，在牛苦头地区M3和M5磁异常区具有寻找华力西期岩浆岩有关的铅锌矿床。
多数学者认为，中—晚三叠世是祁漫塔格造山带内非常重要的地质演化阶段（丰成友等，2012；Yu et al.，2017），是祁漫塔格地区主要的成矿时期。该阶段下，构造环境由挤压变为伸展，有利于岩浆混合作用的发生，并进行分异演化，为矿区内铁锌铜多金属矿化提供了较好的构造条件（高永宝等，2014）。
然而，也有部分学者对“印支期作为祁漫塔格地区唯一的矽卡岩成矿期”提出质疑，并通过事实和证据，提出该地区多金属成矿与泥盆纪岩浆作用有关 （高永宝等，2014；宋忠宝等，2014； 李加多等，2019）。但未深层次揭露岩体与矿体时空关系以及成矿作用过程。
测试结果表明，M1矿段钻孔底部成矿花岗闪长岩年龄为（362.2±2.7）～（361.8±3.4）Ma。 与与闪锌矿密切共生黄铁矿Re–Os等时线年龄为（359.2±6.3）Ma，年龄与矿区M1和M4华力西晚期钻孔底部成矿岩体年龄（363～362 Ma）基本一致，代表了祁漫塔格地区晚泥盆世存在一期重要的矽卡岩型铅锌成矿作用。结合前人发表的有关有关牛苦头矿区成矿岩体年龄（365～352 Ma）（李加多等，2019；王新雨等，2021），认为由此确定牛苦头铅锌矿床形成于362～352 Ma，属于华力西期晚泥盆世。
综合以上分析，认为牛苦头矿区M1、M4矿段大规模铅锌成矿作用时代为华力西期晚泥盆世。华力西期岩浆岩在整个牛苦头矿区普遍都存在，尤其在M3和M5地区。因此，在牛苦头地区M3和M5磁异常区具有寻找华力西期岩浆岩有关的铅锌矿床。
这个事实也进一步说明了祁漫塔格地区存在大规模“华里西期铅锌成矿作用”，为祁漫塔格地区区别于“印支期成矿作用的”的另一期中酸性岩浆岩找矿标志。
5.2 成矿物质来源
在岩浆演化及成矿流体运移过程中，由于Os亲铁、亲铜的特性，其可作为示踪成矿物质来源的元素，由于牛苦头矿床闪锌矿矿石中包含大量黄铁矿。因此，可以利用黄铁矿的Os同位素特征示踪牛苦头矿区闪锌矿的物质来源（丰成友等，2007；李超等，2009）。其中，地幔中187Os/188Os值约为0.12，而上地壳中的187Os/188Os值约为1，牛苦头黄铁矿–闪锌矿矿石的初始（187Os/188Os）i值为0.13±0.24，显示出幔源岩浆参与了成矿。矽卡岩矿床一般由成矿岩体与地层发生交代而形成，地层和成矿岩体可能同时贡献了成矿物质。牛苦头海西期岩浆岩为祁漫塔格地区晚泥盆世后碰撞造山环境下岩石圈减薄过程中壳幔混合作用的产物（莫宣学等，2007；高永宝等，2014），牛苦头铅锌多金属矿床成矿物质可能主要来源于壳幔作用的混合。
5.3 矿床成因
东昆仑造山带经历了4次岩浆构造旋回，形成的岩石构造组合基本上对应于4个时代，分别为：元古宙（前寒武纪）、早古生代、晚古生代—早中生代、晚中生代—新生代，其中以晚古生代—早中生代（华力西期—印支期）的岩浆岩较为发育（莫宣学等，2007；田龙等，2023）。各构造旋回相关的成矿作用主要集中在晚古生代—早中生代，其中东昆仑晚古生代成矿包括夏日哈木岩浆岩型铜镍矿床、白干湖矽卡岩–云英岩–石英脉型钨锡多金属矿床，矿床形成时代为430～422 Ma，属于志留世。该期成矿作用形成于志留纪碰撞造山后局部拉张环境（Zhong et al.，2018）。进入泥盆世（410～360 Ma），形成一系列花岗闪长岩–石英闪长岩–二长花岗岩类，与区内的矽卡岩型铁铅锌铜多金属成矿密切相关，这些侵入岩以 I 型花岗岩为主。目前，发现该时期矿床较少，近年来该时代矿床也逐渐被重视，其代表矿床包括牛苦头–四角羊铅锌多金属矿床、野马泉M13磁异常区铁锌矿床等（高永宝等，2014；王新雨等，2021），形成于区内早古生代—晚古生代早期构造–岩浆旋回的碰撞–后碰撞阶段，由古老陆壳重熔而成，加入部分地幔物质，壳幔岩浆混合作用可能是其大规模集中成矿的主要因素。
结合上述Re–Os同位素成矿物质来源分析，可以认为牛苦头矽卡岩型铅锌多金属矿床形成于晚泥盆世祁漫塔格洋俯冲于柴达木地块之下后的碰撞–后碰撞阶段，此时背景下，岩石圈减薄、软流圈上涌，强烈的壳幔相互作用诱发了大规模的岩浆活动，这也与牛苦头矿区成矿岩体中含有大量包体的事实相一致（李加多等，2019；王新雨等，2021）。中酸性晚泥盆世花岗闪长岩–二长花岗岩侵入到滩间山群发生了相互交代作用，晚泥盆世花岗岩类提供了成矿物质，岩浆与滩间山群大理岩接触带提供了赋矿空间，从而形成了牛苦头矽卡岩型铅锌多金属矿床。
6. 结论
（1）牛苦头矿区成矿二长花岗岩与花岗闪长岩时代为（362.2±2.7）～（361.8±3.4）Ma。 Re–Os等时线年龄为（359.2±6.3）Ma。成岩与成矿时代耦合，由此确定牛苦头铅锌矿床形成于362～359 Ma。
（2）根据黄铁矿初始值（187Os/188Os）i分析，可以认为牛苦头铅锌矿床金属成矿物质可能主要来源于壳幔相互作用的混合岩浆。
（3）结合已有资料，提出牛苦头铅锌矿床形成于碰撞–后碰撞阶段的拉伸背景之下。
致谢：野外工作得到了青海鸿鑫矿业有限公司技术中心工作人员郭天军、刘明、王燕的大力支持； Re–Os测试技术得到了中国地质科学院国家测试中心李超教授的悉心指导；审稿专家对论文提出了许多宝贵的意见和建议，对以上人员表示最衷心的谢意。
-
![]()
图 2 塔河油田9区下油组沉积柱状剖面图
Figure 2. Sedimentary column profile of lower oil formation in block 9 of Tahe oilfield
![]()
图 3 塔河油田9区S100井下油组构型单元及界面单井柱状图
Figure 3. Configuration unit and interface single well column diagram of S100 downhole oil formation in block 9 of Tahe oilfield
![]()
图 4 塔河油田9区三叠系下油组三角洲前缘储层构型模式图（据徐丽强等，2016修改）
Figure 4. The braided river delta reservoir architecture model of lower oil formation in block 9 of Tahe Oilfield
![]()
图 6 塔河油田9区三叠系下油组各类型夹层岩电特征
Figure 6. Electrical characteristics of various types of intercalations in the lower Triassic oil formation in block 9 of Tahe oilfield
![]()
图 7 塔河油田9区三叠系下油组1-1单层二级渗流屏障平面分布图
Figure 7. Plane distribution of 1-1 single layer secondary seepage barrier of Lower Triassic formation in block 9 of Tahe oilfield
![]()
图 8 塔河油田9区三叠系下油组三级渗流屏障遮挡剖面图
Figure 8. Seepage shielding profile of Lower Triassic oil formation in block 9 of Tahe Oilfield
![]()
图 9 塔河油田9区三叠系下油组1-1单层四级渗流屏障平面分布图
Figure 9. Plane distribution of 1-1 single layer four stage seepage barrier of Lower Triassic oil formation in block 9 of Tahe oilfield
![]()
图 10 基于构型控制的块状厚油层渗流屏障分布地质模型
①.一级渗流屏障：封闭断层；②.二级渗流屏障，5 级界面；③.三级渗流屏障，4 级界面（相界面）；④：四级渗流屏障：3 级界面（层内夹层）
Figure 10. Geological model of seepage barrier distribution in massive thick reservoir based on configuration control
![]()
图 11 塔河油田9区三叠系下油组1-1单层渗流屏障分布图
Figure 11. Distribution of single–layer seepage barrier of Lower Triassic formation 1-1 in block 9 of Tahe oilfield
表 1 渗流屏障级次及特征
Table 1 Seepage barrier grades and characteristics
渗流屏障
级别定义 对应构型界面级别 通体级别 一级 封闭性断层及垂向不同期次沉积体系间的非渗透性边界 6级（区域稳定沉积的泥岩隔层、不整合） 连通体 二级 同期水道与心滩坝沉积组合之间
非渗透性边界5级（沉积体系组合之间的泥岩或者泛滥平原沉积） 连通单元 三级 同期水道与心滩项之间的界面 4级（泥质河道、泥质半充填河道和洪水
漫流细粒沉积）连通单元 四级 心滩坝内部的非渗透性沉积 3级（心滩坝内部落淤层，坝上沟道 渗流单元 
下载: 导出CSV
表 2 不同构型单元非均质性参数统计表
Table 2 Statistical table of heterogeneity parameters of different configuration units
层号 辫状河道 心滩 水下/水上分流河道 水下/水上河道侧缘 平均值 变异
系数突进
系数级差 平均值 变异
系数突进
系数级差 平均值 变异
系数突进
系数级差 平均值 变异
系数突进
系数级差 T2a1-1 — — — — — — — — 103 0.6 4 49.3 93.7 0.7 6.7 563.5 T2a1-2 68.3 0.42 1.7 3.6 — — — — 120.3 0.52 4 69.4 140.9 0.62 3 14.9 T2a1-3 221.9 0.55 2.5 11.5 136.1 0.4 1.9 7.4 150.6 0.5 3 34.3 106.5 0.52 2.9 10.7 T2a1-4 207.6 0.58 3.6 32 160.4 0.56 3.5 51.3 — — — — — — — — 平均 165.9 0.52 2.6 15.7 148.3 0.48 2.7 29.3 125.4 0.54 3.7 51 113.7 0.61 4.2 196.4
下载: 导出CSV
-
陈欢庆, 赵应成, 舒治睿, 等. 储层构型研究进展[J]. 特种油气藏, 2013, 20(5): 7-12 doi: 10.3969/j.issn.1006-6535.2013.05.002 CHEN Huanqing, ZHAO Yingcheng, SHU Zhirui, et al. Advances in reservoir architecture research[J]. Special Oil & Gas Reservoirs, 2013, 20(5): 7-12. doi: 10.3969/j.issn.1006-6535.2013.05.002
段冬平, 侯加根, 郭素华, 等. 塔河油田九区三叠系油气藏隔夹层识别及其展布研究[J]. 科技导报, 2010, 28(19): 21-25 DUAN Dongping, HOU Jiagen, GUO Suhua, et al. Identification and distribution of interlayer in Triassic reservoir in block 9 of Tahe Oilfield [J]. Science and Technology Guide, 2010, 28 (19): 21-25.
付国民, 赵俊欣, 杨磊, 等. 塔河油田9区三叠系中上统高分辨率层序地层及沉积演化[J]. 兰州大学学报: 自然科学版, 2007, 43(3): 13-17. FU Guomin, ZHAO Junxin, YANG Lei, et al. High resolution sequence ofstratigraphy and sedimentary evolution of middle-upper Triassic in the9th area, Tahe oil field[J]. Journal of Lanzhou University ( NaturalSciences) , 2009, 43( 3) : 13 - 17.
付国民, 周丽梅, 刘蕊, 等. 塔河三叠系下油组河流相储层夹层成因类型及其分布规律[J]. 地球科学与环境学报, 2009, 31(3): 260-264 doi: 10.3969/j.issn.1672-6561.2009.03.006 FU Guomin, ZHOU Limei, LIU Rui, et al. Fluvial Facies Reservoir Interbed Genesis Category and Distribution Characteristic in Low Oil Group Triassic Tahe Oilfield[J]. Journal of Earth Sciences and Environment, 2009, 31(3): 260-264. doi: 10.3969/j.issn.1672-6561.2009.03.006
郭建华, 刘辰生, 朱锐. 阿克库勒地区三叠系层序地层学及储集砂体成因类型[J]. 沉积学报, 2007, 25( 2) : 169 - 177. . GUO Jianhua, LIU Chensheng, ZHU Rui. Sequence stratigraphy and sandbody genetic types of Triassic system in Akekule area[J]. Acta Sedmentologica Sinica, 2007, 25( 2) : 169 - 177.
何拓平, 李元昊, 陈朝兵, 等. 深水重力流储层宏观非均质性控制因素-以华庆地区长63为例[J]. 西北地质, 2020, 53(1): 177-188 HE Tuoping, LI Yuanhao, CHEN Zhaobin, et al. Macroscopic Heterogeneity Controlling Factors of Deepwater Gravity Flow Reservoirs[J]. Northwestern Geology, 2020, 53(1): 177-188.
贺婷婷, 段太忠, 赵磊, 等. 塔里木盆地塔河油田T区三叠纪沉积模式[J]. 石油与天然气地质, 2019, 40(4): 822-833 HE Tingting, DAUN Taizhong, ZHAO Lei, et al. Triassic sedimentary model in Block T of Tahe oilfield, Tarim Basin[J]. Oil & Gas Geology, 2019, 40(4): 822-833.
贺婷婷, 段太忠, 赵磊, 等. 塔河油田九区三叠纪下油组夹层识别及分布规律[J]. 东北石油大学学报, 2017, 41(6): 26-35 doi: 10.3969/j.issn.2095-4107.2017.06.004 HE Tingting, DUAN Taizhong, ZHAO Lei, et al. Identification and distribution of interlayer in Lower Triassic oil formation in block 9 of Tahe Oilfield [J]. Journal of Northeast Petroleum University, 2017, 41 (6): 26-35. doi: 10.3969/j.issn.2095-4107.2017.06.004
李顺明, 宋新民, 蒋有伟, 等. 高尚堡油田砂质辫状河储集层构型与剩余油分布[J]. 石油勘探与开发, 2011, 38(4): 474-482 LI Shunming, SONG Xinmin, JIANG Youwei, et al. Architecture and remaining oil distribution of the sandy braided river reservoir in the Gaoshangpu Oilfield[J]. Petroleum Exploration and Development, 2011, 38(4): 474-482.
孙天建, 穆龙新, 吴向红, 等. 砂质辫状河储层构型表征方法——以苏丹穆格莱特盆地 Hegli 油田为例[J]. 石油学报, 2014a, 35(4): 715–734 doi: 10.7623/syxb201404012 SUN Tianjian, MU Longxin, WU Xianghong, et al. A quantitative method for architectural characterization of sandy braided- river reservoirs: taking Hegli oilfield of Muglad Basin in Sudan as an example[J]. Acta Petrolei Sinica, 2014a, 35(4): 715–734. doi: 10.7623/syxb201404012
孙天建, 穆龙新, 赵国良. 砂质辫状河储集层隔夹层类型及其表征方法—以苏丹穆格莱特盆地Hegli油田为例[J]. 石油勘探与开发, 2014, 41(1): 112–120 SUN Tianjian, MU Longxin, ZHAO Guoliang. Classification and characterization of barrier-intercalation in sandy braided river reservoirs: Taking Hegli Oilfield of Muglad Basin in Sudan as an example[J]. Petroleum Exploration and Development,2014, 41(1): 112–120.
万琼华, 罗伟, 梁杰, 等. 基于储层构型的流动单元渗流屏障级次研究[J]. 西南石油大学学报(自然科学版), 2019, 41(1): 77-84 WAN Qionghua, LUO Wei, LIANG Jie, et al. Reservoir Architecture-based Classification of Seepafe Barriers of Flow Unit[J]Journal of Southwest Petroleum University (Science & Technology Edition), 2019, 41(1): 77-84.
王珂, 戴俊生, 贾开富, 等. 塔河油田1区三叠系储层流动单元研究[J]. 岩性油气藏, 2014, 26(3): 120-130 doi: 10.3969/j.issn.1673-8926.2014.03.020 WANG Ke, DAI Junsheng, JIA Kaifu, et al. Research on reservoir flow units of Triassic in block-1, Tahe Oilfield[J]. Lithologic Reservoirs, 2014, 26(3): 120-130. doi: 10.3969/j.issn.1673-8926.2014.03.020
王石, 万琼华, 陈玉琨, 等. 基于辫状河储层构型的流动单元划分及其分布规律[J]. 油气地质与采收率, 2015, 22(5): 47-51, 68 doi: 10.3969/j.issn.1009-9603.2015.05.008 WANG Shi, WAN Qionghua, CHEN Yukun, et al. Flow units division and their distribution law based on braided river reservoir architecture[J]. Petroleum Geology and Recovery Efficiency, 2015, 22(5): 47-51, 68. doi: 10.3969/j.issn.1009-9603.2015.05.008
吴胜和 , 岳大力 , 刘建民 , 等. 地下古河道储层构型的层次建模研究[J]. 中国科学: D 辑: 地球科学, 2008, 38(增刊Ⅰ): 111-121. WU Shenghe, YUE Dali, LIU Jianmin, et al. Hierarchy modeling of subsurface palaeochannel eservoir architecture[J]. Science in China: Series D: Earth Sciences, 2008, 51(Supp. Ⅱ): 126-137.
吴胜和. 储层表征与建模[M]. 北京: 石油工业出版社, 2010: 136–174 WU Shenghe. Reservoir characterization and modeling[M]. Beijing: Petroleum Industry Press, 2010: 136-174.
徐丽强, 李胜利, 于兴河, 等. 辫状河三角洲前缘储层构型分析——以彩南油田彩9井区三工河组为例[J]. 油气地质与采收率, 2016, 23(5): 50-57, 82 doi: 10.3969/j.issn.1009-9603.2016.05.008 XU Liqiang, LI Shengli, YU Xinghe, et al. Analysis of reservoir architecture in the braided river delta front: A case study of the Sangonghe Formation in Block Cai9 of Cainan oilfield[J]. Petroleum Geology and Recovery Efficiency, 2016, 23(5): 50-57, 82. doi: 10.3969/j.issn.1009-9603.2016.05.008
Kelly S. Scaling and hierarchy in braided rivers and their deposits: C[M]//Sambrook Smith G H, Best J L, Bristow C S, et al. Braided rivers: Process, deposits, ecology and management. Oxford: Blackwell Publishing, 2006: 75–106.
Leeder M R. Fluviatile fining upwards cycles and the magnitude of paleochann[J]. Geological Magazine, 1973, 110(3): 265-276. doi: 10.1017/S0016756800036098
Miall A D. Architectural elements analysis: A new method of facies analysis applied to fluvial deposits[J]. Earth Science Review, 1985, 22(4): 261-308. doi: 10.1016/0012-8252(85)90001-7
Miall A D. Architectural elements and bounding surfaces in fluvial deposits: anatomy of the Kayenta formation (lower jurassic), Southwest Colorado[J]. Elsevier, 1988, 55(3–4).
-
期刊类型引用(5)
1. 何佳乐，唐卫东，刘天航，高永宝，魏立勇，范堡程，于锴. 龙首山构造带罗城花岗闪长岩地质年代学、岩石成因和构造背景研究. 西北地质. 2025(01): 52-67 . 
本站查看
2. 王利，梁天意. 内蒙古东部钓鱼台地区火山岩年代学、地球化学特征及地质意义. 西北地质. 2025(01): 68-80 . 本站查看
3. 向浩予，刘松，康波，陈昌军，邓伟，邓修林，陈浩如. 班公湖-怒江成矿带西段白板地北部晚侏罗世花岗闪长岩锆石U-Pb年龄、微量元素组成及地质意义. 西北地质. 2025(01): 43-51 . 本站查看
4. 宋彦博，王继春，袁硕浦，王建平，翟德高，柳振江，许展，王东星，邱广东，杨永亮. 内蒙古下陶勒盖铁铜矿床成矿岩体地球化学特征及其成岩成矿启示. 西北地质. 2024(02): 90-106 . 本站查看
5. 程会，程旭，王新雨. 青海省祁漫塔格牛苦头铅锌矿床氧逸度特征与铅锌成矿作用关系研究. 地质与勘探. 2024(05): 863-881 . 百度学术
其他类型引用(0)
计量
- 文章访问数: 89
- HTML全文浏览量: 21
- PDF下载量: 26
- 被引次数: 5
