高级检索

    川西孝泉—丰谷地区须二段深层致密砂岩储层地应力大小评价方法及其工程应用

    黄滔, 刘岩, 何建华, 叶泰然, 邓虎成, 李瑞雪, 李可赛, 张家维

    黄滔,刘岩,何建华,叶泰然,邓虎成,李瑞雪,李可赛,张家维. 2024. 川西孝泉—丰谷地区须二段深层致密砂岩储层地应力大小评价方法及其工程应用[J]. 中国地质, 51(1): 89−104. DOI: 10.12029/gc20230306005
    引用本文: 黄滔,刘岩,何建华,叶泰然,邓虎成,李瑞雪,李可赛,张家维. 2024. 川西孝泉—丰谷地区须二段深层致密砂岩储层地应力大小评价方法及其工程应用[J]. 中国地质, 51(1): 89−104. DOI: 10.12029/gc20230306005
    Huang Tao, Liu Yan, He Jianhua, Ye Tairan, Deng Hucheng, Li Ruixue, Li Kesai, Zhang Jiawei. 2024. Evaluation method and engineering application of in−situ stress of deep tight sandstone reservoir in the second member of Xujiahe Formation in Xiaoquan−Fenggu area, western Sichuan[J]. Geology in China, 51(1): 89−104. DOI: 10.12029/gc20230306005
    Citation: Huang Tao, Liu Yan, He Jianhua, Ye Tairan, Deng Hucheng, Li Ruixue, Li Kesai, Zhang Jiawei. 2024. Evaluation method and engineering application of in−situ stress of deep tight sandstone reservoir in the second member of Xujiahe Formation in Xiaoquan−Fenggu area, western Sichuan[J]. Geology in China, 51(1): 89−104. DOI: 10.12029/gc20230306005

    川西孝泉—丰谷地区须二段深层致密砂岩储层地应力大小评价方法及其工程应用

    基金项目: 深地科学与工程教育部重点实验室开放基金项目(DESEYU202102)与四川省科技厅重点苗子项目(22MZGC0159)联合资助。
    详细信息
      作者简介:

      黄滔,男,1998年生,硕士生,主要从事非常规油气储层地应力场精细描述研究;E-mail: 1136830431@qq.com

      通讯作者:

      刘岩,女,1983年生,副教授,主要从事储层特征及储层评价技术研究;E-mail: liuyan08@cdut.cn

    • 中图分类号: P618.13

    Evaluation method and engineering application of in−situ stress of deep tight sandstone reservoir in the second member of Xujiahe Formation in Xiaoquan−Fenggu area, western Sichuan

    Funds: Supported by the Open Found from Key Lab of Deep Earth Science and Engineering, Ministry of Education (No.DESEYU202002), Sichuan Provincial Science and Technology Department Key Seedling Project (No.22MZGC0159).
    More Information
      Author Bio:

      HUANG Tao, male, born in 1998, master candidate, engaged in the study of fine description of the in−situ stress field; E-mail: 1136830431@qq.com

      Corresponding author:

      LIU Yan, female, born in 1983, associate professor, mainly engaged in reservoir characteristics and reservoir evaluation technology; E-mail: liuyan08@cdut.cn.

    • 摘要:
      研究目的 

      川西坳陷孝泉—丰谷地区须二段砂岩气藏的勘探开发潜力巨大,但该地区埋藏较深且构造复杂、断缝系统多期叠加,使得地应力频繁变化,制约了该区井位轨迹设计与压裂改造的有效实施,故需对该区地应力大小进行精细评价,为工程开发提供建议从而提高产能。

      研究方法 

      基于岩石力学、声发射实验及差应变分析等实验测试方法,并结合常规测井、特殊测井及水力压裂等资料分析,优选了适应于深层块状均质致密砂岩储层的地应力大小实验测试方法,并在单点地应力大小准确评价的基础之上,构建了研究区分构造变形单元分层的单井地应力大小连续测井解释模型,查明了纵向上地应力大小变化结构类型及分布规律。

      研究结果 

      研究表明差应变分析法计算的地应力大小精确度最高,为更能够准确表征深层均质块状致密砂岩地应力大小的实验方法。测试结果显示须二段属于Ⅲ类地应力类型,处于走滑应力状态,存在部分逆冲挤压应力状态;在井点测试的基础上,形成了基于井壁影像反演的地应力大小评价技术;地应力大小结构变化在纵向上分为5种类型,其中南北向(SN)三级以上断层和南北向(SN)褶皱或北东东向(NEE)褶皱变形的高部位以低低高(LLH)型和低应力(LC)型为主,在小规模断层或平缓构造区以高低高(HLH)型或高低低(HLL)型为主。

      结论 

      建议选择低低高(LLH)型地应力剖面进行工程开发,其纵向上可穿透更多含气层,同时避开底层底水,预防生产早期快速见水,故应选择二—三级南北向断层和南北或北东向纵弯褶皱区须二2中上段进行水力压裂改造。

      创新点:

      形成了基于井壁影像反演的地应力大小评价技术;结合不同地应力剖面类型和水力压裂缝形态对须二段进行储层改造。

      Abstract:

      This paper is the result of oil and gas exploration survey engineering.

      Objective 

      Xiaoquan−Fenggu area in the western Sichuan Basin has huge potential for exploration and development of Xujiahe Formation gas reservoirs. However, due to the region's deep burial, complex structure, and multiple superimposed fault systems that cause frequent variation in stress orientation, effective well placement design and hydraulic fracturing practices have been limited. Therefore, it is necessary to evaluate the precise magnitude of in−situ stress in this area to provide recommendations for engineering development and increase production capacity.

      Methods 

      Experimental methods such as rock mechanics, acoustic emission testing, and differential strain analysis combined with conventional logging, special logging, and hydraulic fracturing data were used to experimentally test the in−situ stress magnitude in deep heterogeneous blocky tight sandstone reservoirs. Based on accurate evaluation of single−point in−situ stress magnitude, a logging interpretation model was established for the subdivision of tectonic units in the study area, examining structural variation of in−situ stress magnitude and distribution along a single well.

      Results 

      Our study showed that differential strain analysis provided the most accurate measurement of stress in heterogeneous tight sandstone reservoirs. Test results indicated that the Xujiahe Formation belongs to the type III in−situ stress category and exists in a strike−slip stress state with partial compression and thrust stress states. Based on single−point test data, we developed a technique to evaluate in−situ stress magnitude by utilizing borehole image inversion. Structural changes in in−situ stress magnitude were divided vertically into five types, whereby high positions of north−south (SN) faults with grades above three and folds in SN or northeast−trending (NEE) resulted predominantly in low−low−high (LLH) or low stress (LC) profiles. Meanwhile, small−scale faults or gentle deformation areas had high−low−high (HLH) or high−low−low (HLL) profiles.

      Conclusion 

      A low−low−high (LLH) stress profile was suggested for engineering development to penetrate more gas layers vertically, while avoiding bottom water and preventing rapid water breakthrough during early production. Therefore, it is recommended to select the second to third order north−south(SN) trending faults and north−south(SN) or northeast trending(NE) longitudinal flexure zones located in the middle−upper part of the second layer of the second member of Xujiahe Formation.

      Highlights:

      The in−situ stress evaluation technology based on borehole wall image inversion is formed. Combined with different in−situ stress profile types and hydraulic fracture morphology, the reservoir reconstruction of the second member of Xujiahe Formation was carried out.

    • 依据寄主围岩的不同,天然金刚石矿化类型包括金伯利岩型、钾镁煌斑岩(又称橄榄金云火山岩)型、蛇绿岩型、碱性基性岩型、超高压变质岩型、陨石相关型和冲积型等(张培元, 1998Shirey et al., 2013Yang et al., 2014),前四种统归为幔源型,其中金伯利岩型是产出宝石级金刚石(钻石)最重要的原生矿床。金伯利岩是一种自然界分布稀少、岩浆起源最深的浅成-喷出岩(Heaman et al., 2003),属于由幔源物质和含挥发分熔体/流体组成、富集相容和不相容元素的钾质偏碱性超基性岩,呈蓝、绿、灰黑色等色调,具斑状和火山碎屑结构,块状、角砾状和凝灰状等构造,常以岩筒、岩脉和岩床产出(郑建平, 1989Heaman et al., 2003;Becker et al., 2006;路凤香, 2008Brook, 2017雷雪英, 2017)。幔源物质包括岩石捕虏体和晶屑两部分,前者主要以石榴石或尖晶石二辉橄榄岩、尖晶石方辉橄榄岩和纯橄岩为主,还有各类辉石岩、榴辉岩和交代型岩石(Stachel et al., 2004路凤香, 2008);后者主要有橄榄石、镁钛铁矿、石榴石、铬尖晶石、铬透辉石、钙钛矿、金云母等(池际尚等, 1996a路凤香, 2008赵欣等, 2015de Wit et al., 2016许康康等, 2016)。非洲南部金伯利岩被分为玄武岩型和金云母型(由岩相学划分),或克拉通型(on-craton)和克拉通边缘(off-craton)型(由产出位置划分,de Wit et al., 1992;Becker et al., 2006),或GroupⅠ(80~114 Ma)和GroupⅡ(114~150 Ma)型(由侵位时代和地球化学特征划分,Smith, 1983)等不同类型。根据金刚石包裹体的来源和成因,金刚石被分为岩石圈(lithospheric)型和超深(superdeep)型(Shirey et al., 2013)。前者约占含包裹体金刚石的99%,一般被分为地幔橄榄岩(或超基性岩,P)型、榴辉岩(E)型和二辉石岩(W)型(Stachel et al., 2008),其中P型金刚石根据石榴石包裹体中Ca和Cr含量的不同分为方辉橄榄岩、二辉橄榄岩和异剥橄榄岩型(Stachel et al., 2008Koornneef et al., 2017);超深型可以来源于软流圈、地幔过渡带、下地幔甚至核幔边界(Stachel, 2001Kaminsky, 2012陆琦等, 2012)。

      全球金伯利岩的分布与克拉通(或古陆核)密切相关。世界上首粒金伯利岩型金刚石发现于卡拉哈里(Kalahari)克拉通内,于1869年南非金伯利附近的道斯方丹(Dorstfontein)和布特丰坦(Bultfontein)农场(Field et al., 2008)。1887年卡维尔刘易斯(H. Carvill Lewis)教授在给曼彻斯特英国协会的报告中首次用当时英国高级长官金伯利(Kimberley)的名字将含金刚石的岩石命名为“金伯利岩”(Wilson et al., 2007)。中国分别于1965年和1971年在华北克拉通内发现山东蒙阴和辽宁瓦房店(旧称复县)两处含金刚石金伯利岩矿带(郑建平, 1989路凤香, 2008)。相比北美、西伯利亚、非洲南部卡拉哈里等克拉通内的原生矿床,中国金刚石矿产资源比较贫乏,已探明金刚石优质资源已近枯竭,目前可供开采的产地严重不足,金刚石找矿亟待突破(彭艳菊等, 2013刘陟娜等, 2016)。蒙阴和瓦房店含矿金伯利岩的时代为中奥陶世(池际尚, 1996b张宏福等, 2007Yang et al., 2009),而此时的华北克拉通并未发生岩石圈减薄(吴福元等, 1999郑建平, 2009朱日祥等, 2011),其岩石圈厚度和结构与非洲南部卡拉哈里克拉通类似。那么,中国金伯利岩型金刚石资源的“匮乏”的原因是什么,是否有进一步找矿的可能?

      笔者多次考察了蒙阴、瓦房店和应县金伯利岩带,2016年8月参加第35届国际地质大会,会前考察了南非的奥兰治—瓦尔(Orange-Vaal)河流域的金刚石砂矿,2017年9月参加了在哈博罗内举行的第11届国际金伯利岩会议,并野外考察了朱瓦能(Jwaneng)、奥拉帕(Orapa)、莱特拉卡内(Letlhakane)、丹姆沙(Damtshaa)和卡罗韦(Karowe)金伯利岩矿区。本文在前人研究的基础上,归纳了中国不同类型金刚石的分布,系统对比了博茨瓦纳和中国金伯利岩型金刚石矿床的区域地质、矿床特征及成因,总结了找矿标志,以期为探寻中国类似金刚石矿产资源提供素材。

      博茨瓦纳人口约233万,面积581730 km2,2/3国土面积被卡拉哈日(Kalahari)沙漠覆盖(Brook, 2017)。从英国殖民统治结束的1966年之前,博茨瓦纳是非洲第三穷的国家,而自1967年该国首次发现金刚石以来,尤其1980年至今较稳定政局的投资环境,使得金刚石产业发展为该国三大经济支柱之一(其他为畜牧业和旅游)(Brook, 2017De Wit, 2018)。

      博茨瓦纳目前金刚石产量居世界第二,总价值居世界第一(de Wit et al., 2016)。黛比瓦纳(Debswana)钻石有限公司是博茨瓦纳最有影响力的钻石开采和销售公司,由政府和黛比尔斯(De Beers)公司各占50%股份,控制了最有价值的朱瓦能和奥拉帕两大含矿金伯利岩区。博茨瓦纳金刚石勘察工作最早始于1955年,目前共发现12个金伯利岩区,根据发现的先后顺序,分别为卡特伦(Kgatleng)、奥拉帕、朱瓦能、Kikau-Khutse、科空(Kokong)、察邦(Tshabong)、莱科都(Lekgodu)、Ghaghoo(2013年之前称为Gope)、马布阿塞胡贝(Mabuasehube)、勒雷拉(Lerala,又称Martin`s Drift)、奥卡瓦(Okwa)和Nxsunxau金伯利岩区(Brook, 2017)(图 1)。目前正在开采的金刚石矿共8个分布在4个金伯利岩区中,按照开采时间包括奥拉帕(1971年,A/K1)、莱特拉卡内(1975年,D/K1和D/K2)、朱瓦能(1982年,DK2)、丹姆沙(2003,B /K1、B/K9、B/K12、B/K15)、勒雷拉(2008,K002~K006)、B/K11(2010)、卡罗韦(2012,A/K6),Ghaghoo(2013),前七个为露天开采,后一个为地下开采(Brook, 2017De Wit, 2018)(图 1)。

      图  1  博茨瓦纳12个金伯利岩群在行政区图和构造地质简图中的分布(据Brook, 2017修改)
      Figure  1.  Administrative and geological maps of Botswana depicting locations of 12 kimberlite clusters

      博茨瓦纳金伯利岩型金刚石矿床的分布与非洲南部克拉通密切相关,该克拉通主要由三部分组成,即卡普瓦尔(Kaapvaal,3.7~2.7 Ga)克拉通和津巴布韦克拉通(3.5~2.6 Ga,Kusky, 1998),以及古元古代的奥卡瓦杂岩体(Okwa Inlier),前两者被2.7~2.6 Ga林波波(Limpopo)构造带(De Wit et al., 1992)连结形成卡拉哈里(Kalahari)古陆(Miensopust et al., 2011)。其周边拼贴一些古、中元古代增生带(Jacobs et al., 2008),包括迈冈迪—开斯—奥卡瓦带(Makondi-Kheis-Okawa belt,MKO)、雷霍博特区(Rehoboth Province)、纳马夸—纳塔尔带(Namaqua-Natal belt)和达马拉(Damaran)造山带等(图 1)。

      MKO位于卡拉哈里古陆西北缘,其中迈冈迪带具有津巴布韦克拉通岩石圈根(Luguet et al., 2015),分布于津巴布韦克拉通西北部,为北东走向的古元古代活动带,由一系列2.06~1.96 Ga绿片岩相至麻粒岩相变质岩组成,其源岩包括(1)约2.26 Ga Deweras群,为津巴布韦克拉通被动陆缘裂解相关的沉积火山岩组合;(2)2.20 Ga Lomagundi群,为一套弧后初始盆地相浅海沉积岩;(3)2.16 Ga Pirowiri群,为陆缘深水大陆架和大陆坡沉积岩。迈冈迪带于2.0~1.85 Ga拼贴到津巴布韦克拉通,随后进入被动陆缘沉积阶段直到中元古代ChomaKalomo地块增生,强烈伸展并伴随大量火山活动发生在1.1 Ga(Jacobs et al., 2008)。开斯带位于卡普瓦尔克拉通西南缘,时代(2.1~2.0 Ga)与迈冈迪带类似,代表被动陆缘初始裂解至大洋形成的火山沉积过程。奥卡瓦(Okwa)杂岩体位于卡普瓦尔克拉通西北缘,上部被新元古代的杭济—乔贝带(GhanziChobe)不整合覆盖。奥卡瓦杂岩体主要由变流纹岩、糜棱化眼球状花岗岩、富钾长石二长花岗岩和细晶花岗岩组成,侵位时代约为2.06 Ga (Mapeo et al., 2006),其上部被低级变质碎屑沉积岩和灰岩构造覆盖。虽然3个带的变质变形时代均主要介于2.1~1.96 Ga,但部分学者认为它们可能不属于一个带,如Moen et al.(2008)认为迈冈迪带与津巴布韦克拉通外缘的林波波带相连,后者发育2.65 Ga和2.0 Ga两期高级变质事件。

      雷霍博特区(Rehoboth Province)位于MKO带西部,最长边约700 km,由太古宙地核(2.98~2.70 Ga)和古元古代花岗岩类组成,与卡普瓦尔克拉通拼合时间在2.2~2.05 Ga(Van Schijndel et al., 2011)。上部发育1.93~1.85 Ga的陆源碎屑岩盆地沉积,之后经历了长达450 Ma的岩浆平静期,直到1.35 Ga发生纳马夸(Namaqua)被动陆缘裂解,1.30~1.20 Ga形成岛弧和弧后盆地岩浆事件,约1.20 Ga西侧岛弧和地体与雷霍博特碰撞增生到雷霍博特—卡普瓦尔克拉通(Van Schijndel et al., 2011)。

      纳马夸—纳塔尔带(Namaqua-Natal belt),在博茨瓦纳西北部该带又称为杭济—乔贝带(GhanziChobe belt) (Modie, 2000),是陆陆碰撞或弧陆碰撞的高压变质带,时间发生在1.2~1.0 Ga之间(Jacobs et al., 2008),出露大量花岗片麻岩和变沉积岩。纳马夸带主要为1.29~1.24 Ga具有岛弧性质的变火山沉积岩,内部侵入大量1.2~1.08 Ga花岗岩类,其上部被Koras群不整合覆盖,后者为陆内伸展环境下未变形和低级变质的双峰式火山岩和不成熟的碎屑岩夹少量碳酸盐岩等沉积岩组成,形成时代为1.17~1.10 Ga。纳塔尔带于1.14 Ga逆冲到卡普瓦尔克拉通东北缘,它包括1.25~1.10 Ga俯冲相关的岛弧岩浆事件(Mccourt et al., 2006),1.07~1.03 Ga弧陆碰撞相关的奥长花岗岩和紫苏花岗岩,伴生高温麻粒岩相变质作用(Eglington, 2006)。

      达马拉(Damaran)造山带是刚果克拉通与卡拉哈日克拉通于550 Ma碰撞后的界限(Kinabo et al., 2008),主要以角闪岩相的花岗片麻岩穹隆以及530~510 Ma碰撞后高钾超铝质S型花岗岩为特征,上部被新元古代变沉积岩覆盖。

      除了上述前寒武系结晶基底外,博茨瓦纳金伯利岩重要的上部围岩为石炭—侏罗系卡鲁超群(Karoo Supergroup,300~183 Ma)页岩、粉砂岩和砂岩等陆源沉积岩和侵入其中的基性岩墙群。卡鲁超群分为5个地层单元,包括Tlapana, Thlabala, Mosolotsane, Ntane, Stormberg组,覆盖了非洲南部大部分区域,记录了潘吉亚超大陆聚合至早期裂解过程,形成石炭—二叠纪弧后前陆盆地和183~174 Ma地幔柱岩浆事件,后者形成非洲南部著名的卡鲁大陆溢流玄武岩省(Hastie et al., 2014)。

      奥拉帕金伯利岩筒群(Orapa kimberlite cluster, OKC)位于古元古代迈冈迪活动带内,接近于迈冈迪和林波波带的交汇部位(图 1),金刚石的形成受到古、中元古代增生过程的影响(Aulbach et al., 2017)。OKC岩筒群由88个岩筒组成,整体呈北西西方向展布(Mccourt et al., 2004)(图 2)。该区域航空磁测资料显示强烈的北西西向异常,结合地表分布北西向线性构造和岩墙群,推测OCK金伯利岩的侵入可能受北西西向断裂控制(Brook, 2017)。

      图  2  奥拉帕金伯利岩群中金刚石矿的分布图(Kruger et al., 2017)
      Figure  2.  Positions of the diamondiferous mines in the Orapa kimberlite cluster

      奥拉帕镇是奥拉帕、莱特拉卡内和丹姆沙矿山的行政中心,是一个人口超过1.2万人的封闭城镇(图 2)。OKC目前有7个矿区在开采,即黛比瓦纳公司的奥拉帕A/K1、莱特拉卡内的D/K1和D/K2、丹姆沙的B/K9和B/K12,以及卢卡拉(Lucara)钻石公司的卡罗韦A/K6和火石(Firestone)钻石公司的B/K11矿区。奥拉帕A/K1矿区位于首都哈博罗内以北约550 km处,位于弗朗西斯敦(Francistown)以西240 km处;莱特拉卡内D/K1矿区距离A/K1东南50 km,距弗朗西斯敦约190 km,D/K1和D/K2相距约800 m;卡罗韦A/K6距离D/K1约25 km,距A/K1约23 km,以产大颗粒优质Ⅱ a型钻石而闻名(Wudrick et al., 2017)。

      (1) 奥帕拉A/K1

      A/K1矿区是OKC最有名的矿区,于1967年由Manfred Marx领导的戴比尔斯地质团队发现,以产量大和产出大粒优质金刚石著称,年产量平均约1200万克拉,每年贡献了约15%的GDP。A/K1由南、北2个岩筒组成,两岩筒地表覆盖面积约118英亩(约477529 m2),南北长约2 km,东西宽约1 km,采深约250 m(图 3ab),2026年预计深达450 m,是世界上最大的露天开采矿区。岩筒侵入太古代英云闪长质和花岗质片麻岩结晶基底和上部的卡鲁火山沉积岩超群,岩筒最上部为一系列火山口湖泊相的表层碎屑金伯利岩,包括含有白垩纪动植物化石的页岩、砂岩、砂砾岩和混杂碎石层,地表大部分被卡拉哈日沙漠覆盖(Rayner et al., 1991)。

      图  3  奥拉帕A/K1金刚石矿的3D模型及其野外地质特征
      a—A/K1的露头采坑;b—A/K1的3D模型显示由两个金伯利岩筒组成(缩写见正文);c—北筒NPK−普遍可见花岗岩和玄武岩岩屑,岩屑粒径多小于1cm,橄榄石晶屑普遍蛇纹石化和方解石化;d—北筒MVK2−BBX,普遍可见棱角状玄武岩角砾;e—南筒黑色块状火山碎屑金伯利岩(SDVK),比重较大,含大量新鲜橄榄石和少量玄武岩捕虏体;f—南筒南部火山碎屑金伯利岩(SVK)上层,可见大量岩屑和橄榄石斑晶,偶见单斜辉石和钛铁矿
      Figure  3.  Field photos and a 3D view of the A/K1 kimberlite in which the lithologic units are shown in the text in details
      a-Orapa A/K1 open pit containing North and South pipes; b-3D model of the A/K1 kimberlitic pipes; c-Northern Proclastic Kimberlite (NPK) in North pipe showing pervasively granite and basalt debris with size less than 1 cm. Olivine is commonly altered by serpentinization and calcitization; d-Massive Volcaniclastic Kimberlite (MVK)− Basalt breccia xenolith (BBX) in North pipe; e-Southern Dark Volcaniclastic Kimberlite (SDVK) including an amount of olivines in the South Pipe; f-Basalt breccia, olivine, rare clinopyroxene and ilmenite are commonly observed in the upper unit of the Southern Volcaniclastic kimberlite (SVK)

      北筒岩石组合主要包括北部玻基斑状火山碎屑(Northern Pyroclastic)金伯利岩(NPK,绿灰色,厚层状构造,含粒径<10 mm的花岗岩基底岩屑(图 3c),橄榄石晶屑普遍蛇纹石化和少量碳酸盐化,局部可见石榴石)、黑色北部玻基斑状火山碎屑金伯利岩(NPKD,为NPK中的薄层,富含黑绿-黑灰色橄榄石,相比其它类型岩石,其岩屑极少)、富含花岗质片麻岩角砾的玻基斑状火山碎屑金伯利岩(NPK-GG,暗绿色,层状构造,富含大量玄武岩和基底角砾)、块状火山碎屑金伯利岩(MVK2,暗灰绿色金伯利岩,含少量蚀变基底角砾)和角砾块状火山碎屑金伯利岩(MVK2-BBX,可见大量棱角明显的玄武质岩屑,图 3d),外缘含有垮塌成因的玄武质角砾岩透镜体。其中NPK类似于MVK2,也类似于火山通道相凝灰质金伯利质角砾岩,指示岩筒总体剥蚀程度低。

      南筒地表直径约800 m,岩石组合相对复杂,包括近水平层状火山口相和近直立火山通道相的岩石单元(图 3b),主要有黑色异相岩屑角砾岩(Dark Heterolithic Breccias (DHB),是南筒中特有的一种岩石类型,含大量砂岩和泥岩岩屑及少量玄武岩和金伯利岩角砾,基质主要由细粒砂/泥和碳酸盐组成),南部黑色块状火山碎屑岩金伯利岩(SDVK,黑绿色-灰绿色,呈致密块状,含有大量橄榄石斑晶,少量岩屑为玄武岩,图 3e),南部火山碎屑金伯利岩(SVK,分为上中下部分:上层为灰绿色层状金伯利岩,分选差至分选显著,含橄榄石和岩屑捕虏体,橄榄石斑晶含量约11%,粒径多在5 mm以上,可见金属光泽的钛铁矿(图 3f);中层为块状绿色金伯利岩,中粗粒基质胶结,岩屑呈半棱角状,伴随碳酸盐化,碳酸盐脉发育;下层黑绿色块状金伯利岩,相比中层,花岗质片麻岩等基底岩屑含量明显减少),玄武角砾岩(BBR,分选较差的棱角状玄武岩,多由碳酸盐充填)。

      南北岩筒的金伯利岩具有许多相似之处,如普遍蛇纹石化和碳酸盐化,局部可见新鲜的橄榄石斑晶,围岩捕虏体主要为卡鲁超群的玄武岩、泥岩和少量基底花岗岩(图 3c);地表被层状岩屑碎石堆(Bedded Talus Cone,A3T,图 3b)覆盖。然而根据南筒层状火山碎屑岩不整合穿切北筒块状火山碎屑金伯利岩指示北筒相比南筒喷发稍早(Gernon et al., 2009),北筒块状火山碎屑金伯利岩基质中含有透辉石微晶,而南筒未见,表明南筒金伯利岩具有较高的侵位温度,玄武岩岩屑的剩磁研究指示北筒的北筒温度为500~600℃,而南筒在<380℃(Gernon et al., 2009)。南筒的稀释率(dilution,即围岩捕虏体体积含量,42%)明显高于北筒(10%)(Brook, 2017)。

      金伯利岩中幔源捕虏体、金刚石及其包裹体是研究上地幔最直观的样品(Stiefenhofer et al., 1997)。奥拉帕A/K1金伯利岩中捕虏体主要为榴辉岩和少量辉石岩(Aulbach et al., 2017)。两组榴辉岩的单斜辉石和石榴石成分均记录了津巴布韦克拉通边缘多次增生过程中未被改造的洋壳印记:一种为单斜辉石87Sr/86Sr初始比值较低(0.7026~0.7046),单斜辉石87Sr/86Sr与火成岩低压熔融参数(如Eu/ Eu*、Mg#、Sr/Y)相关,可能与海底低温蚀变导致的较大范围δ18O值(单斜辉石:4.83‰~8.36‰;石榴石:5.05‰~8.18‰),这些证据指示了古洋壳岩石(辉长岩)经历了分离结晶和海底低压蚀变;另一种为单斜辉石87Sr/86Sr>0.7048~0.7091,而单斜辉石和石榴石均具有地幔型的δ18O值(5.1‰~5.9‰),指示受到少量来自古老的金云母富集地幔流体交代作用。该交代作用具体表现在不同类型榴辉岩中单斜辉石的87Sr/86Sr值上,即含金刚石榴辉岩87Sr/86Sr初始比值(0.7039)接近于N-MORB(0.703),但明显低于受流体交代的含石墨和不含石墨的榴辉岩中单斜辉石(分别为0.7064、0.7051),指示少量流体/熔体的氧化-交代作用可能不仅导致了金刚石向石墨的转化,还催化了碳的消失。该推理与含金刚石和含石墨榴辉岩和辉石岩捕虏体已移出金刚石稳定域(150~200 km(Shirey et al., 2013))的计算结果(单斜辉石-石榴石压力计:3.2~4.9 GPa(Aulbach et al., 2017))相吻合。

      奥拉帕A/K1金伯利岩主要产出E型金刚石,且至少存在两种年龄不同的IaAB型金刚石,以及两种Ⅱ型金刚石(Chinn et al., 2018)。老金刚石(形成于1.7Ga,T=1175℃)N的含量、N的聚合方式以及生长环带样式较年轻金刚石(形成于1.1 Ga,T=1150℃)多变,老金刚石表现为N含量较高,δ13C =-28.6‰~-2.1‰,δ15N=-10.4‰至+23.0‰,而年轻的金刚石δ13C =-9.6‰~-3.3‰,δ15N=-10.2‰~-4.4‰,该特征与环带明显的金刚石从核部到边部δ13C值逐渐增加一致(Chinn et al., 2018)。具有高N和复杂N聚合形式、低δ13C老金刚石的形成与俯冲沉积物流体有关,而年轻金刚石的形成可能与来自软流圈的流体有关(Chinn et al., 2018)。与含金刚石榴辉岩相比,金刚石的榴辉质包裹体的FeO含量高,SiO2含量低,CaO和Na2O含量低,它们记录了高平衡温度和石榴石具有地幔型δ18O同位素特征,这一特征是由一种来自富含单斜辉石岩陆下岩石圈地幔熔体造成的,可能与布什维尔德2.0 Ga的岩浆事件有关(Aulbach et al., 2017)。

      (2) 莱特拉卡内(Letlhakane)

      莱特拉卡内金刚石矿发现于1968年,包括D/ K1和D/K2两个金伯利岩筒,分别11.6英亩和3.6英亩(分别46943 m2和14569 km2),为不同喷发期次的岩筒。2003年该矿生产了106万克拉的金刚石,可采矿石品位约为0.26克拉/t;2015年矿山年采矿石210多万t,处理矿石230多万t,分选了58.3万克拉金刚石。莱特拉卡内矿山的露天作业于2017年停止,尾矿加工厂已建成投产,投产后将使用寿命延长至2043年,年产量将达80万克拉。莱特拉卡内金伯利岩的围岩与奥拉帕A/K1类似,下部为太古代花岗质片麻岩结晶基底,上部为卡鲁超群,顶部被厚4~10 m不等的卡拉哈日沙漠覆盖。

      根据金伯利岩的构造和蚀变程度、岩屑大小和含量、橄榄石颗粒大小和含量,D/K1共分为6种岩相:(1)火山碎屑金伯利岩(VK1)厚度超过600 m,主要由块状和似层状火山碎屑岩及围岩角砾岩组成,围岩粒径颗粒较大分层明显,角砾磨圆稍差;(2)黑色火山碎屑金伯利岩(DVK),为含有大量蛇纹石化橄榄石斑晶(体积含量>50%)、岩屑(体积含量<25%),基质呈黑灰至黑色的坚硬火山碎屑岩;(3)灰白色火山碎屑金伯利岩(PVK),呈灰白色,富集围岩角砾岩;(4)斑点状火山碎屑金伯利岩(SVK),由含粒径为5~15 mm的角砾状花岗质岩屑和黄棕色蚀变橄榄石斑晶(体积含量约30%)和浅灰色至黑棕色基质组成;(5)玄武角砾岩(BBX),由围岩重力不稳引发崩塌和泥石流等过程产生的沉积物,呈浅棕灰色,橄榄石斑晶含量10%~20%,完全蛇纹石化,岩屑含量60%~90%以上,粒径10~20 cm,最大可达2 m;(6)相关金伯利岩(Coherent Kimberlite, CK),含量较少,呈黑灰色至黑色,块状构造,橄榄石含量高达60%以上,整体新鲜,局部颗粒具有富钙边缘,岩屑含量较少(约7%),普遍发生流体交代作用(Kruger et al., 2017)。

      D/K2仅于1985—1986年期间短暂开采,D/K2岩相主要包括巨晶金伯利岩(MK,呈黑灰色,含相当多新鲜橄榄石,岩屑主要为玄武岩和花岗片麻岩)、火山碎屑金伯利岩(VK,呈灰绿色,橄榄石普遍蛇纹石化呈黑绿色,玄武岩、泥岩等岩屑呈拉长状蚀变强烈)、围岩角砾岩(WRB,弱角砾状碎屑支撑的金伯利岩,可见大的岩块)、异相角砾岩(HB,基质至碎屑支撑的外来岩屑,由卡鲁地壳岩石和金伯利岩混合而成)(Kruger et al., 2017)。

      莱特拉卡内D/K1岩筒金伯利岩中常见地幔橄榄岩、辉石岩、榴辉岩、云母岩、MARID(由云母、角闪石、金红石、钛铁矿、透辉石组成)地幔岩等捕虏体,其中地幔橄榄岩普遍较新鲜,粒径一般在5~10 cm,主要由含石榴石二辉橄榄岩、不含石榴石二辉橄榄岩、尖晶石方辉橄榄岩组成,与金伯利岩有明显的反应边,指示普遍经历了交代作用。地幔橄榄岩中石榴石和斜方辉石常被金云母、单斜辉石、钛铁矿和尖晶石取代形成石榴石金云母地幔橄榄岩和金云母地幔橄榄岩(Luguet et al., 2015)。尖晶石地幔橄榄岩的Re-Os同位素TRD年龄(3.75±0.54)Ga (Wainwright et al., 2015)和2.5~2.8 Ga(Luguet et al., 2015),代表了迈冈迪造山带下部克拉通根的稳定时代,指示该克拉通根可能是津巴布韦克拉通的西延部分(Luguet et al., 2015Wainwright et al., 2015)。

      金刚石中硫化物包裹体的Re-Os同位素的TMA年龄为1.06~2.38 Ga,与相邻奥拉帕金刚石的年龄类似(1.0~2.9 Ga)。金刚石N含量变化较大,从很低含量(21×10-6~43×10-6)至较高的聚合N(53%~90% ⅠaB)。碳同位素δ13C值为-19.3‰~-22.7‰。阴极发光图像中没有明显环带构造的两粒金刚石的硫化物等时线年龄分别为(0.92 ± 0.23)Ga(初始187Os/188Os=1.31 ± 0.24)和(0.93 ± 0.36)Ga(初始187Os/188Os=0.69±0.44),指示含有循环Os元素的金刚石结晶环境极度不均一。硫化物包裹体球粒陨石标准化的PGE元素图解富集Ir(3.4~33)和Pt (2.3~28.1),指示一种富集地幔过程。硅酸盐包裹体的Sm-Nd等时线年龄为(2.3 ± 0.02)Ga、(1.0 ± 0.14)Ga和(0.25±0.04)Ga,以上特征指示D/K1金刚石经历了多阶段生长过程(Gress et al., 2017)。

      (3) 丹姆沙(Damtshaa)

      丹姆沙矿床距离奥拉帕A/K1矿点约15 km,包括B/K1、B/K9、B/K12、B/K15等金伯利岩筒,发现于1967—1972年间,开采于2003年(Deines et al., 2009),其中B/K9和B/K12是两个最大矿点,也是正在开采的矿点,两者分别占地15、8英亩(分别为60703 m2、32375 m2Field et al., 2008)。丹姆沙矿床预计能开采31年,总产量约500万克拉。金伯利岩围岩与奥拉帕和莱特拉卡内相似,下部为太古宙花岗质片麻岩,上部为二叠—侏罗系卡鲁超群的泥岩、砂岩和德拉肯斯堡(Drakensberg)组的拉斑玄武岩。

      B/K9矿点在地下250 m处分为北、中、南岩筒(图 4),3个岩筒在地表合并在一起,其中南、北2个岩筒岩岩相类似,主要由浅成-半深成相黑色火山碎屑金伯利岩(DVKns,ns指南北岩筒)组成,边部夹少量玄武质角砾岩(BB)透镜体及玄武岩和火山碎屑金伯利岩巨大块体;中部岩筒由浅层火山口相火山沉积物和BB,及深部DVKc(c指中部岩筒)组成。DVKns呈特征的黑绿-灰色,岩屑和晶屑散布在黑色细粒基质中,岩筒内部分层不明显呈块状,而边部层状构造明显。晶屑主要为橄榄石和少量钛铁矿和金云母,大部分蛇纹石化、水榴石化和碳酸盐化等;岩屑包括大量地幔捕虏体和其他次要岩石,岩屑粒径从0.5 cm至10~26 cm,岩屑比例在1.5%~10%,最大可达18%。BBL边部最厚可达20 m,夹少量DVK,岩层倾向岩筒内部15~20°,向中心逐渐尖灭。DVKc呈棕色至蓝绿色,以含大量橄榄石为特征(40%~45%)。玄武岩岩屑和地幔捕虏体约占7%。晶屑普遍具有核边结构,核部大多为新鲜的橄榄石,少量为金云母和钛铁矿,边部为蛇纹石、钛磁铁矿、水榴石和钙钛矿等蚀变矿物。总之,B/K9南、北岩筒的两种DVK呈块状构造,零星分布的岩屑由基质支撑,岩石组成随深度几乎不发生变化,说明DVK的形成经历了充分的混合,但不完全均匀化,岩屑粒径多为几厘米,最大岩屑不超过10~20 cm(Buse et al., 2011)。

      图  4  丹姆沙B/K9矿点岩筒地质简图(据Buse et al., 2011修改)
      Figure  4.  Schematic geological map illustrating the features of the Damtshaa B/K9(modified from Buse et al., 2011)

      B/K12矿点包括3个火山碎屑金伯利岩单元和一个富含玄武质角砾岩的中部带,整体研究程度较低,捕虏体、金刚石特征等方面还没有报道(Field et al., 2008)。

      Deines et al.(2009)报道了丹姆沙95粒含包裹体金刚石中55粒为P型、38粒为E型和2粒为W型。其中P型包裹体主要有橄榄石、斜方辉石、铬镁铝榴石,E型有镁铝榴石-铁铝榴石和绿辉石;W型有斜方辉石、单斜辉石和石榴石包裹体。55个P型δ13C值为-1.44‰~-9.20‰,38个E型δ13C值为-3.78‰~-17.60‰,2个W型的δ13C值为-4.91‰~-10.63‰。所有类型的金刚石具有一个-4.75‰的δ13C峰(-5.0‰~-4.5‰, Deines et al., 2009)。

      丹姆沙E型金刚石的δ13C值和N含量与莱特拉卡内没有明显差异,而与奥拉帕具有较大差异。相比奥拉帕,丹姆沙E型金刚石非常富集δ13C,含有较高的N含量,最高超过1000×10-6δ13C为-5‰的金刚石比更低δ13C金刚石具有更高的N含量。相比奥拉帕和莱特拉卡内,丹姆沙含有橄榄石的P型金刚石具有更高的δ13C值,橄榄石的MgO含量也稍高于莱特拉卡内,但是总体来说3个矿区的Mg#值变化不大。丹姆沙橄榄石的CaO含量稍高于奥拉帕。金刚石的铬尖晶石成分不同,丹姆沙铬尖晶石的MgO含量稍低于奥拉帕,而后者的MnO含量很高(0.77%)。相比奥拉帕,丹姆沙石榴石包裹体Cr2O3、MgO含量和Mg#较高,而Al2O3、CaO含量较低(Deines et al., 2009)。

      W型金刚石中的石榴石的Cr含量和Mg#值位于P和E型之间。δ13C值介于-21.3‰~-3.9‰。W型与P型石榴石成分存在过渡,δ13C值明显不同;W与E型石榴石成分存在差异,但两者的δ13C值重合。如果低δ13C值是壳源成分引起的,那么来自地壳的石榴石的Al2O3/Cr2O3比值大于地幔。然而低Al2O3/Cr2O3比值的金刚石,δ13C值也低(<-7‰),而Al2O3/Cr2O3比值高的金刚石,δ13C值约为-5‰。该推理与传统认为E型金刚石的低δ13C值是有机碳(δ13C平均为-25‰)俯冲的结果相反,说明地幔中可能存在低Al2O3/Cr2O3、低δ13C值的地幔区域。以上数据说明,相比奥拉帕,丹姆沙和莱特拉卡内P型石榴石包裹体源自非常亏损的纯橄岩-方辉橄榄岩区域(CaO<1.8%),较高Cr含量,指示源自更深地幔。奥拉帕金刚石来自160 km深度,丹姆沙和拉特拉卡内金刚石可能分别来自190 km和200 km (Deines et al., 2009)。

      (4) 卡罗韦(Karowe)A/K6

      卡罗韦A/K6金刚石矿(25° 28'13"E,21° 30' 35"S)是近年来世界上产出大颗粒优质钻石的重要矿区之一,尤其以出产大颗粒IIa型(N<10×10-6)钻石闻名(Moore, 2014)。Lucara金刚石公司于2015年11月在该矿发现世界第二大钻石(1109克拉的Lesedi La Rona,或Our Light in Setswana),仅次于1905年南非发现的3106克拉“库里南”,此外813克拉的“星座”钻石成交价达到创纪录的6310万美元。该矿最早由De Beers公司1969年发现,由于较小的地表出露面积(3.3英亩,约13355 m2,图 5a~c)和较低的工业品位(3.5 ct/100 t),于1998年被放弃。A/K6矿于2010年由Lucara公司完全控股,开采于2012年,计划露天开采至2026年。2012年6月至2016年12月,该矿已产出180万克拉金刚石,售价超过10.2亿美元,2017年产出29~31万克拉,优质大颗粒金刚石估计能开采15年(Campbell et al., 2017),2018年的开采32.5~35万克拉,税收约1.9亿美元,计划开采深度324 m,矿石3620万t,储量630万克拉;如果采深达400 m,将有5100万t矿石和820万克拉的资源量,400~750 m之间的2100万t矿石含有400万克拉资源量

      卡罗韦矿海拔1022 m,大致呈南北走向,地下120 m处的面积约7英亩(约28328 m2),地表多被卡拉哈日沙漠覆盖,其围岩与奥拉帕类似。A/K6包括北、中、南3个金伯利筒,然而由于地表强烈蚀变,无法识别3个岩筒的界限,故将岩筒顶部区域划为一个独立单元,即钙质和硅质胶结砂砾岩(图 5a),其厚度一般为10 m,局部可达20 m。岩筒上部10~50 m为强烈风化金伯利岩(WK),主要矿物组成和结构构造已难鉴别。风化程度随深度增加而减小,新鲜金伯利岩一般出露在70~90 m以下。玄武质角砾岩(BBX),混合有少量(通常<10%)的金伯利岩,呈不连续分布于三个岩筒中,粒径从几米至金伯利岩基质中的毫米级不等,普遍具有不同程度的破裂和碳酸盐脉充填(Lynn et al., 2014)。

      图  5  卡罗韦A/K6金刚石矿床岩筒模型和野外地质特征
      a—卡罗韦A/K6南、中、北岩筒模型剖面图(南筒模拟深度约300 m(Chinn et al., 2010));b—AK6的3个金伯利岩矿筒中产出50~100克拉金刚石的分布图;c—AK6露头采坑;d—采坑的玄武质围岩;e—南筒的玻基斑状火山碎屑金伯利岩(VK)和块状初始岩浆金伯利岩(MPK)
      Figure  5.  Geological model and field photographs of diamondiferous kimberlite features in the Karowe mine
      a-Geological section of A/K6 mine showing the features of north, central and south lobes; b-Distribution map of diamonds with 50 to 100 carats from north, central and south lobes; c-Open pin of the Karowe (A/K6); d-Country rock of basalt from south lobe; e-Fragmental kimberlite (FK) sample from the central lobe and magmatic/pyroclastic kimberlite (MPK) sample from the south lobe

      北筒主要为破碎金伯利岩(FK),淡绿灰色,火山角砾粒径主要以4~32 mm不等,基质支撑,分选较差的块状,基质中普遍蛇纹石化和方解石化。玄武岩是主要的围岩捕虏体,另含少量卡鲁沉积岩和花岗片麻岩结晶基底。中筒与北筒大体相似,均含有非碎屑的岩浆金伯利岩和火山碎屑金伯利岩,两者结构和颜色变化较大交替出现,但彼此间为渐变过渡关系。围岩捕虏体均以卡鲁玄武岩为主。然而相比北筒,中筒金伯利岩碳酸盐含量较高(即富含碳酸盐的碎屑金伯利岩,CFK,图 5a)。CFK通常比FK含有更高的橄榄岩大晶体和更低的围岩捕虏体,基质中不透明矿物的含量也略高。BBX常出现于FK而非CFK中。南筒在岩相和地球化学等方面明显不同于北中筒。南筒上部、西部以风化玄武质角砾岩(WBBX)为主,下部为未蚀变玄武质角砾岩(BBX)和大块玄武岩体(图 5d)。南筒主要由中粗粒(4~32 mm)块状、基质支撑、分选差的粗晶岩浆金伯利岩或玻基斑状火山碎屑金伯利岩(M/PK)充填。M/PK宏观上呈灰色,含5%~10%热交代/蚀变的围岩捕虏体。橄榄石相对新鲜,含有丰富的不透明矿物。新鲜的钙镁橄榄石(monticellite)随着深度的增加而丰富。围岩捕虏体以玄武岩为主,基底和卡鲁沉积岩较少,但整体地壳稀释率很低(一般<10%),很少达到25%。局部可见少量厚层状PK,厚0.16~1.5 m不等,由橄榄岩等晶体和近水平定向排列的地壳捕虏体聚集而成(Lynn et al., 2014)。

      南筒与北中筒具有不同的岩相(图 5a)和喷发期次,北中筒的内部结构构造复杂岩相,而南部则呈均一的块状,多产出大颗粒优质钻石,10.8克拉以上的钻石在北筒的比例为1.8%,而南筒达到5%(图 5b)。南筒的地表直径相对较大(约280 m×180 m,图 5, 图 c),岩相学特征与北、中筒明显不同,主要岩相为岩浆/玻基斑状金伯利岩(M/PK)为主(图 5e),橄榄石、钛铁矿等斑晶均较新鲜,且随深度越深越新鲜,围岩捕虏体含量小于10%。

      A/K6出露大量地幔橄榄岩捕虏体而榴辉岩较少,该特征与奥拉帕相反。Wudrick et al.(2017)报道了卡罗韦矿区19个地幔橄榄岩捕虏体,其中14个为尖晶石二辉橄榄岩,5个为石榴石二辉橄榄岩,石榴石二辉橄榄岩中含有韭闪石指示经历地幔交代作用。橄榄石的Fo值=91~93,与卡普瓦尔克拉通的橄榄石Fo平均值(92.6)一致,放射性187Os/188Os = 0.1260~0.1295,Re-Os模式年龄(TDR)=2.6~2.0 Ga,与南非同样产Ⅱa型Premier矿床年龄(2.0 Ga)类似,可能指示卡普瓦尔克拉通中晚太古代岩石圈地幔为亏损型地幔,并叠加了布什维尔德侵入体式岩浆交代或约1.8 Ga迈冈迪带形成时的交代作用(Wudrick et al., 2017)。

      金刚石中矿物包体在限定金刚石和地幔性质及其关系具有重要意义。Motsamai et al.(2018)从A/K6矿区107颗金刚石中发现15粒氧化物、22粒硫化物和134粒硅酸盐,其中53%为E型,44%为P型,1%为W型。P型石榴石的CaO含量较低(<1.8 wt%)而Cr2O3含量较高(14.7 wt%),橄榄石的Fo平均值为93.1,指示金刚石形成时克拉通岩石圈非常亏损且深度大于220 km。3个E型石榴石具有镁铁榴石成分,指示这种金刚石来自330~420 km的深度,位于软流圈和过渡带边界部位(Motsamai et al.(2018))。

      朱瓦能矿是目前世界上价值最高的钻石矿床,发现于1972年的纳勒迪河谷(Naledi River Valley,又称“星星谷”)中,于1982年8月全面投产,贡献了黛比瓦纳公司的60%~70%的年收益。朱瓦能岩筒群位于哈博罗内以西约160 km,包括约28个金伯利岩筒,海拔1192 m。黛比尔斯公司于1971年发现DK1,1972年发现DK2、DK7,目前DK2是主要的开矿区,DK7还未开采。DK2采坑地表面积约54英亩(约218530 m2图 6a),岩筒壁较陡(80°以上),包含11个岩筒(de Wit et al., 2016),主要由北、中、南3个岩管和隐伏4号筒组成(Brown et al., 2008)。

      图  6  朱瓦能DK2矿区的4个金伯利岩管分布图及野外特征
      a—DK2露头采坑;b—DK2采坑平面地质简图(de Wit et al, 2016);c—南筒火山碎屑金伯利岩和玄武岩围岩露头
      Figure  6.  Jwaneng D/K2 pit composed of north, central, south and forth lobes (a); Schematic geological map showing the DK2 features (b); Photographs of the volcanoclastic kimberlite and wall rock basalt from the south lobe (c)

      DK2金伯利岩的围岩为新太古—古元古代(2.65~2.05 Ga)德兰士瓦(Transvaal)超群白云岩、炭质和石英质页岩等硅酸盐-碳酸盐岩组合,厚度>500 m(Creus et al., 2018),其不整合覆盖于太古宙坎耶(Kanye)火山岩、花岗岩及绿岩带上,与布什维尔德杂岩体同期(2.05 Ga,Creus et al., 2018),地表被卡拉哈日沙漠沉积物覆盖(de Wit et al., 2016)。火山口相火山碎屑金伯利岩厚度在600 m以上,上部被沉积层理发育的再造火山碎屑物覆盖,下部为块状火山碎屑金伯利岩组成了岩筒的主体(图 6b),火山碎屑岩发育陡峭的斜层理和剪切构造,边部可见围岩角砾。

      DK2岩筒的侵位先后顺序为南、中、北筒。南、中筒主体为块状火山碎屑金伯利岩,上部被再造火山碎屑金伯利岩覆盖。边部主要为含有围岩角砾金伯利岩和剪切金伯利岩,还包括卡鲁超群下部的泥岩残余,这表明金伯利岩曾侵入卡鲁下部沉积层,但后来遭到了侵蚀。南筒主要由致密块状基本不发育层理构造的VK组成,其边部可见石英富集的火山碎屑金伯利岩(图 6c)。中筒的火山碎屑金伯利岩整体蚀变强烈,橄榄石多已蚀变为蛇纹石和富Mg黏土,颜色由浅灰绿色至深绿色,岩石碎屑含量较低(<5%~15%),主要为中细粒的石英质页岩、泥岩和白云岩,石榴石和钛铁矿捕虏晶常见,地幔橄榄岩和榴辉岩捕虏体较少,在5%以下。中筒西南角含有石英质页岩和白云岩角砾,与火山碎屑金伯利岩相接触截然(Mmualefe, 2017)。中筒含有一块宽楔状的围岩角砾岩,是金伯利岩侵位时深部爆炸形成的(de Wit et al., 2016)。

      北筒不同于南、中筒,其核部主要由厚层火山碎屑金伯利岩组成,夹较细火山碎屑金伯利岩薄层,上部为薄层再造火山碎屑(RVK)和氧化金伯利岩(OK),下部块状玻基斑状(PK)和火山碎屑金伯利岩(VK)组成。其中PK主要黑色致密块状碎屑岩石和蛇纹石化橄榄石斑晶、不同含量细小岩屑组成,形成致密块状结构。

      博茨瓦纳金刚石交易公司对朱瓦能2013— 2017年间生产的近13万粒0.1~0.8克拉金刚石统计显示,大多为八面体和菱形十二面体,仅5%的金刚石含有包裹体(粒径大于10 μm),其中>88%为E型,6%为P型。包裹体颜色和结构具有以下特征:榴辉质镁铝榴石-铁铝榴石呈橙色、绿辉石质单斜辉石为灰绿色、蓝晶石为蓝色,地幔橄榄岩型镁铝榴石呈紫色,铬尖晶石为棕色至黑色,硫化物包体周围普遍可见特征的玫瑰状破裂(Gress et al., 2018)。CL和FTIR研究显示,E型和P型金刚石至少分别经历了4次和3次生长过程。该过程分别对应了卡普瓦尔克拉通大约1.0 Ga、1.5 Ga和2.9 Ga演化过程(年龄来自硫化物包裹体Re-Os等时线和P型石榴石和单斜辉石包体的Sm-Nd等时线),后者与南非金伯利岩筒中俯冲相关的P型硫化物年龄2.89 Ga相吻合,记录了东、西卡普瓦尔克拉通在Colesberg磁异常构造带聚合过程(Richardson et al., 2004Gress et al., 2018)。

      中国金刚石矿化类型丰富,不仅产出大量幔源型金刚石,还出露超高压变质型和冲积型矿床。其中产于铬铁矿和地幔橄榄岩中幔源型的蛇绿岩型金刚石目前已在新特提斯雅鲁藏布江缝合带(东波、普兰、当穷、休古嘎布、日喀则、泽当、罗布莎)、中特提斯班公湖—怒江带(丁青、东巧)、古亚洲洋贺根山蛇绿岩带和萨尔托海蛇绿岩带中发现(Yang et al., 2015);幔源型的碱性基性岩型金刚石目前仅出露于安徽栏杆(蔡逸涛等, 2018),超高压变质型主要分布于柴北缘、南阿尔金、北秦岭、大别—苏鲁、松多、西南天山、喜马拉雅西构造结等构造带中(杨经绥等, 2009, 及内部文献)。冲积型矿床主要出露于郯庐断裂区沂水河流域(产出158.786克拉的常林钻石)和复州河—岚崮河流域,以及湖南常德一带沅水流域,在大别山区、新疆喀拉喀什河及西藏安多等也有分布(涂怀奎, 2001)。

      中国金伯利岩型和钾镁煌斑岩型金刚石矿床(化)分布于华北克拉通、华南克拉通和塔里木陆块中。华北克拉通与金刚石相关的岩体主要为金伯利岩,后两者主要为钾镁煌斑岩。华北克拉通目前共发现12个金伯利岩岩群和2个钾镁煌斑岩岩群,分属于3个金伯利岩带:(1)郯庐金伯利岩带,沿郯城—庐江深大断裂系分布,包括辽宁桓仁、铁岭、瓦房店(旧称复县)、葫芦岛(齐玉兴等, 1998宋瑞祥, 2013),吉林集安(张复顺, 1992),山东蒙阴(刘继太, 2002)金伯利岩岩群以及吉林通化金厂(张复顺, 1992)钾镁煌斑岩岩群;(2)华北中央金伯利岩带,分布于东、西华北地块1.85 Ga聚合形成的华北中央造山带(Zhao et al., 2005)内,包括山西应县水沟门(杨建民, 1995)、大同采凉山(杨建民, 1995)、柳林尖家沟(魏云峰等, 2015),河北涉县(河北省地质矿产局综合地质大队, 1988),河南鹤壁(侯广顺等, 2016)金伯利岩岩群以及山西大同饮牛沟钾镁煌斑岩岩群(杨建民, 1995);(3)华北北缘金伯利岩带,主要沿华北克拉通北缘的柳河—铁岭—赤峰—白云鄂博断裂带分布,包括内蒙古四子王旗龙头山(中国地质科学院地质研究所, 1970)金伯利岩岩群(图 7)。其中有经济价值的金伯利岩型金刚石矿床为郯庐带中的蒙阴和瓦房店。

      图  7  中国金伯利岩和钾镁煌斑岩分布简图(板块边界和主断裂带据Zheng et al., 2013许文良等, 2013
      1—山东蒙阴;2—辽宁瓦房店(旧称复县);3—辽宁桓仁;4—辽宁葫芦岛;5—辽宁铁岭;6—吉林集安;7—吉林通化钾镁煌斑岩;8—河南鹤壁;9—河北涉县;10—山西大同饮牛沟钾镁煌斑岩;11—山西大同采凉山;12—山西应县;13—山西柳林;14—内蒙古四子王旗;15—湖北大洪山彭家垮金伯利岩;16—湖北大洪山王关−徐家冲钾镁煌斑岩;17—湖南宁乡;18—广西融水;19—广西都安;20—广西大化;21—贵州镇远;22—贵州施秉;23—贵州麻江;24—扬子陆块西缘钾质煌斑岩;25—浙江龙游;26—江西安远;27—新疆巴楚;28—新疆皮山
      Figure  7.  Schematic map of kimberlites and lamproites distribution in China
      1-Mengyin, Shandong; 2-Wafangdian (formerly known as Fuxian), Liaoning; 3-Huanren, Liaoning; 4-Huludao, Liaoning; 5-Tieling, Liaoning; 6-Ji'an, Jilin; 7-Potassium-magnesium lamprophyre in Tonghua, Jilin; 8-Hebi, Henan; 9-Shexian, Hebei; 10-Yinniugou potassiummagnesium lamprophyre in Datong, Shanxi; 11-Cailiangshan, in Datong, Shanxi; 12-Yingxian, Shanxi; 13-Liulin, Shanxi; 14-Siziwang Banner, Inner Mongolia 15-Pengjia collapse kimberlite in Dahongshan, Hubei; 16-Wangguan-Xujiachong K-Mg lamprophyre in Dahongshan, Hubei; 17-Ningxiang, Hunan; 18-Guangxi meltwater; 19-Du'an, Guangxi; 20-Dahua, Guangxi; 21-Zhenyuan, Guizhou; 22-Shibing, Guizhou; 23-Majiang, Guizhou; 24-K-lamprophyre on the western margin of Yangtze block; 25-Longyou, Zhejiang; 26-Anyuan, Jiangxi; 27-Bachu, Xinjiang 28-Pishan, Xinjiang

      华南克拉通钾镁煌斑岩主要分布于四个带中:(1)江南钾镁煌斑岩带,分布于扬子陆块和华夏陆块新元古代聚合形成的江南造山带内,包括贵州镇远—麻江—施秉钾镁煌斑岩和金伯利岩区(罗会文等, 1989任怀翔等, 1993张锡贵等, 2015;)、湖南宁乡钾镁煌斑岩区(林玮鹏, 2007)、广西融水—都安—大化钾镁煌斑岩区(王新宇等, 2013郑翔等, 2016);(2)华南北缘钾镁煌斑岩带,分布于华南克拉通北缘,出露湖北大洪山彭家垮金伯利岩和王关-徐家冲钾镁煌斑岩区(叶德隆等, 1991);(3)扬子陆块西缘钾质煌斑岩区(舒小辛, 1994),分布于华南克拉通西缘;(4)江绍似金伯利质煌斑岩带,出露浙江游龙—江西安远路径似金伯利质煌斑岩(秦正永等, 2001叶松等, 2007),大致沿江绍断裂带分布。其中只有江南带的镇远和宁乡发现一定量的金刚石,但不具经济价值(梅厚钧等, 1998)。

      塔里木陆块被认为是中国金刚石原生矿最有前景的找矿靶区之一(梁日暄等, 1992)。新疆和田地区喀拉喀什河流域的墨玉县已发现11粒金刚石,但自1945年发现塔里木首粒金刚石以来原位金刚石找矿工作一直未取得进展,通过金刚石指示矿物(如铬尖晶石)推测金刚石的寄主岩石可能为钾镁煌斑岩(徐向珍, 2006),该推测与皮山县发现的两条钾镁煌斑岩脉吻合(赵磊等, 1998)。呈岩脉和岩管产出的289 Ma巴楚角砾状超镁铁岩含有大量橄榄石、金云母斑晶及单斜辉石、钙铁榴石、钙钛矿等微晶与金伯利岩类似,但镁铝榴石、铬透辉石含量极低,全岩地球化学高CaO、TFe2O3,低MgO、LREE和(La/Yb)n比值,不同于典型的金伯利岩,被中国学者命名为金伯利质隐爆角砾岩(鲍佩声等, 2009李弦, 2012)。

      中国目前具工业价值的含矿金伯利岩产地仅有山东蒙阴和辽宁瓦房店,两者分布在郯庐断裂带的两侧(图 7),南北距离约550 km(乔秀夫等, 2002田洪水等, 2017)。郯庐断裂带是中国东部一条规模巨大的断裂构造带,向南至少延至苏鲁造山带,向北延至俄罗斯境内,虽然其发育时代、是否发生平移及平移距离、南北段的延伸及演化机制等还存在争议(乔秀夫等, 2002罗志立等, 2005吴根耀等, 2007林宗满, 2011),然而其产状陡倾,切割深达上地幔,为地幔上涌和岩浆贯入的通道,控制(含矿)金伯利岩的侵位等基本达成共识。其活动时代普遍认为包括新元古代—古生代古郯庐断裂带活动阶段和中—新生代活动阶段(乔秀夫等, 2002吴根耀等, 2007田洪水等, 2017)。郯庐断裂带在山东境内又称沂沐断裂带(罗志立等, 2005),长约330 km,宽20~80 km(平均约50 km),主要由鄌郚—葛沟、沂水—汤头、安丘—莒县、昌邑—大店等4条断裂组成,将山东分为鲁西、鲁东次级陆块。鲁西东部(聊城—兰考断裂以东至安丘—莒县以西)发育一系列北西向相间分布的断块隆起和凹陷盆地。蒙阴金伯利岩带位于该区隆断带的中部,走向北北东,长55 km,宽15~20 km,由常马庄、西峪、坡里岩带组成,仅前两个含有原生金刚石(池际尚, 1996b),三带分别受蒙山、新泰—垛庄、铜冶店—蔡庄等北西向断裂限制(罗声宣等, 1999乔秀夫等, 2002)(图 8)。在已有的58个岩体中,含矿岩体45个,20个达到工业品位(宋瑞祥, 2013)。常马庄带位于蒙山地幔凸起区,断续长约14km,宽约8km,由8组岩脉和一对岩管(胜利Ⅰ号)组成。该带不仅分布中国最早发现的金伯利岩体(红旗Ⅰ号岩脉),而且金刚石品位高(最高可达6.35ct/m3),盛产大颗粒金刚石(王久华, 2012)。两个岩管地表相距20 m,地下约250 m合并为一个岩管,其中小管是蒙阴岩带中含矿最富的岩体。大小管岩性主要由粗晶斑状金伯利岩和金伯利角砾岩组成,含少量蛇纹石化碎裂岩(宋瑞祥, 2013)。西峪带位于新甫山地幔凸起和地幔凹陷区的过渡带,分布在高都镇山头、西峪和薛家峪一带,与常马庄带相距14 km,沿5~20°方向展布,长约15 km,宽0.5~1.0 km,由15条岩脉和12个岩管组成,其中红旗5号岩脉具工业意义,主要岩性为斑状金伯利岩、凝灰质金伯利岩。坡里带位于金星山地幔凹陷内(罗声宣等, 1999),分布于岱崮镇至野店镇一带,岩带长约18 km,宽约0.5 km,由25条岩脉组成,未达到工业品位,主要岩性为斑状富金云母金伯利岩(池际尚, 1996b)。

      图  8  山东蒙阴金伯利岩地质简图(据宋瑞祥, 2013修改)
      F1—蒙山断裂;F2—新泰—垛庄断裂;F3—铜冶店—蔡庄断裂;F4—上五井断裂;F5—鄌郚−葛沟断裂;F6—沂水−汤头断裂;F7—安丘—莒县断裂;F8—昌邑—大店断裂
      Figure  8.  Geological map of Mengyin kimberlites showing the locations of three mines of Changmazhuang, Xiyu and Poli and main faults(modified from Song, 2013)
      F1−Mengshan fault; F2−Xintai-Duozhuang fault; F3−TongyedianCaizhuang fault; F4−Shangwujing fault; F5−Zanzenggou fault; F6−Yishui−Tangtou fault; F7−Anqiu−Juxian fault; F8−Changyi−Dadian fault

      蒙阴金伯利岩岩石类型主要有斑状金伯利岩、细粒金伯利岩、金伯利角砾岩、凝灰质金伯利岩、蚀变金伯利岩等。其中金伯利角砾岩主要出露于岩管和局部宽脉的部位,分为灰岩质和花岗质金伯利角砾岩;斑状金伯利岩是构成岩脉和岩管的主要岩石类型,包括斑状金伯利岩、斑状金云母金伯利岩、斑状含镁铝榴石金伯利岩、斑状含岩球金伯利岩、含花岗质围岩金伯利岩等,呈灰绿、墨绿、灰蓝色及暗棕色,圆斑或卵斑结构,块状构造、球状构造,局部为流动构造,斑晶主为橄榄石、次为金云母、镁铝榴石、铬透辉石等。基质为显微斑状结构,斑晶为蛇纹石、金云母、磁铁矿、钙钛矿、磷灰石等。蚀变类型主要有蛇纹石化、碳酸盐化和硅化等。金云母含量从常马庄、西峪、坡里岩带含量逐渐增加(分别为<10%、20%、30%~40%)。3个岩带的金刚石品位和成矿条件不同,含矿性和岩石组合亦有较大差异,从南向北岩带含矿性依次变贫。最南端的常马庄含矿性较好,其中最富的为“胜利Ⅰ号”小管,以产大颗粒的金刚石著称,产出至今最大的原生矿金刚石是“蒙山Ⅰ号”(119.01克拉);岩性以斑状镁铝榴石金伯利岩为主,次为斑状金云母金伯利岩;西峪岩带以斑状金云母金伯利岩为主,次为斑状金伯利岩和斑状含岩球金伯利岩;坡里以斑状富金云母金伯利岩为主。常马庄和西峪岩带中可见少量石榴石纯橄岩、尖晶石方辉橄榄岩、尖晶石二辉橄榄岩、榴辉岩等捕虏体,大多见于岩管中,粒径一般5~10 cm(宋瑞祥, 2013)。

      蒙阴金刚石颜色整体以无色、淡黄、浅棕色为主,次有不同浓度的绿、灰黄、橙等色;形态有单晶、双晶、聚合晶等,以阶梯状八面体和曲面菱形十二面体最常见;包裹体矿物以石墨为主,达到80%~90%,少量可见镁铝榴石、铬尖晶石、橄榄石、铬透辉石、金刚石等。3个岩带中金刚石的颜色、形态和包裹体含量和种类上均有明显差异。常马庄金刚石有色金刚石较多,无色者占36.43%,以菱形十二面体为主,八面体次之,含包裹体金刚石较少(约占32.6%),Ⅰ型金刚石多,Ⅱ型少。相比之下,西峪无色金刚石达75%,坡里金刚石形态以八面体为主,菱形十二面体次之,西峪金刚石介于两者之间(宋瑞祥, 2013)。

      辽宁省金刚石储量全国第一(占52%),目前已在瓦房店、铁岭、桓仁、葫芦岛等4个地区发现了含金刚石金伯利岩和3个冲积型砂矿(宋瑞祥, 2013)。含矿金伯利岩区分属瓦房店—铁岭、庄河—桓仁、绥中—北镇3条金刚石成矿带(赵春强等, 2018),均位于凌源—北票—沙河断裂以北的华北克拉通东北部(赵光慧等, 2011)(图 9)。辽宁金伯利岩型金刚石矿带分布受郯庐深大断裂和金州等数条次级断裂带控制。郯庐断裂带在辽宁境内涉及沈阳—渤海段,又称营潍断裂带(张鹏等, 2010),走向NE 20~30°,南起大连湾,由金州西经瓦房店至沈阳以东和铁岭,分割东侧长期处于隆起区的胶辽地盾和西侧复州凹陷区(池际尚等, 1996a赵建军等, 2011)。

      图  9  辽宁省含金刚石金伯利岩和主要断裂构造分布图(据赵光慧等, 2011宋瑞祥, 2013修改)
      Figure  9.  Distribution map of diamondiferous kimberlites and major faults in Liaoning Province (modified from Zhao et al., 2011; Song, 2013)

      瓦房店金伯利岩型金刚石矿是辽宁省目前唯一达到工业开采要求的原生矿,位于华北克拉通冀辽陆核南部辽东台隆-普兰店地幔隆起区,距郯庐断裂带东30~40 km处。NNE向复州河断裂、松木岛—松树镇断裂和金州断裂系与近EW向岚崮山弧形断裂联合控制了金伯利岩侵入(赵建军等, 2011宋瑞祥, 2013)。瓦房店金伯利岩群分为脉状、管状和隐伏体3种产出类型,近EW向断裂控制着区内大部分金伯利岩管和岩脉的长轴方向,脉状金伯利岩体均赋存于NEE向至近EW向构造带内;而管状金伯利岩体主要出现在NNE断裂与NEE向断裂的交接部位(赵建军等, 2011)。基底由太古宙深变质鞍山群与浅变质辽河群组成,元古宙辉绿岩脉使地块刚性增强,盖层主要为新元古代青白口纪至震旦纪地层。金伯利岩体主要侵入于盖层复式背斜轴部,即新元古代钓鱼台、南芬、桥头组的石英岩、粉砂岩、泥灰岩、页岩和寒武纪毛庄组和徐庄组灰岩、页岩等中(徐俊, 2013)。

      瓦房店金伯利岩群南北长32 km,东西宽约30 km,主要分布于头道沟、老虎屯、李店、二道沟一带,由24个岩管和88条岩脉组成,其中达到工业品位的有6个岩体(30、42、50、51、68、74)。岩群划分为5个金伯利岩带,包括:(1)大马圈子—太阳沟带,长约28 km,宽2 km,包括14个岩管和67个岩脉,中国目前最大的42号岩管和瓦房店第二大岩管(30号)及9、10、11岩脉分布于该矿带中;(2)头道沟—岚崮店带,长约15 km,宽1~2 km;包括8个岩管和17条岩脉,整体品位好价值高,含矿最好的50号岩管分布其中;(3)邓屯—大高家屯带,长约6 km,宽0.25~1 km,含2个岩管和4条岩脉;(4)石灰窑子—李家店带;(5)大王沟—西洼带(图 10),以前3个带为主(宋瑞祥, 2013赵春强等, 2018)。瓦房店金伯利岩岩性复杂,主要包括斑状金伯利岩、斑状金云母金伯利岩(金云母含量5%~15%)、斑状富金云母金伯利岩(金云母含量>15%)、细粒金伯利岩、含岩球斑状金伯利岩、含围岩角砾状金伯利岩、金伯利(凝灰)角砾岩、金伯利凝灰岩和含金伯利物质角砾岩等。热液蚀变主要包括金云母化、蛇纹石化、碳酸盐化、滑石化、硅化、绿泥石化和褐铁矿化,反映了岩浆期后含挥发分热液活动强烈。50号岩管深远捕虏体主要为地幔橄榄岩(尖晶石方辉橄榄岩、石榴石方辉橄榄岩)和金云母岩及少量云母橄长岩。

      图  10  辽宁瓦房店金伯利岩群构造地质简图(据(宋瑞祥, 2013赵春强等, 2018)修改)
      Figure  10.  Sketch geological map of Wafangdian kimberlite cluster in Liaoning Province including five kimberlite zones (modified from Zhao et al., 2011; Song, 2013; Zhao et al., 2018)

      除瓦房店含矿金伯利岩外,辽宁其他3个金伯利岩群金刚石含量低颗粒小,目前不具有工业价值(宋瑞祥, 2013)。桓仁金伯利岩群共含60多个金伯利岩体,以岩脉为主,含金刚石极少,分布于太子河—浑河断裂和庄恒断裂带交汇的桓仁凸起内(宋瑞祥, 2013)(图 9)。金伯利岩脉集中产出在桓仁四平街的湾沟、马圈子、苇塘沟一带,走向30~40°,脉长30~300 m,宽0.5~2.5 m不等,脉岩以斑状金云母金伯利岩为主,有不同程度的蛇纹石化和碳酸盐化,围岩为古元古代混合片麻岩(宋瑞祥, 2013)。铁岭金伯利岩群共含28个金伯利岩体,呈岩脉和岩管产出,走向NE和NNE,分布在铁岭八宝岭—帽园子一带,位于郯庐断裂带东侧30 km,辽东台隆的铁岭—靖宇隆起的边部。铁矿金伯利岩脉一般长100~150 m,最长1260 m,宽约30 m;岩管长30~150 m,宽20~40 m,岩性以斑状金云母金伯利岩为主,围岩为古元古代混合片麻岩、斜长角闪岩等。葫芦岛金伯利岩位于山海关台隆绥中凸起内,基底为太古宙建平岩群大营子组混合花岗岩,位于向要路沟断裂带、北大山—锦西断裂带和药王庙—班吉塔断裂交汇处。金伯利岩主要分布在苇子沟东平台和石灰窑荀家沟,两者相距约25 km,以NEE向岩脉产出。其中苇子沟岩脉侵入寒武系馒头组页岩和灰岩中,脉长40~80 m,宽0.2~2.4 m,岩性以斑状富金云母金伯利岩为主(宋瑞祥, 2013)。

      限定金伯利岩的侵位时代是合理解释地球深部动力学和建立金伯利岩成因的时空模型的关键和难点之一。金伯利岩定年早期方法有金云母和全岩Rb-Sr等时线年龄、金云母K-Ar、40Ar-39Ar定年等,然而这些结果常由于样品后期蚀变、低封闭温度、潜在过剩Ar含量等影响而出现不确定性(Wu et al., 2010Stanley et al., 2016)。由于锆石普遍以捕虏晶产出且常与钛铁矿连生,锆石U-Pb测年结果普遍被认为稍早于金伯利岩实际喷发年龄(Moore et al., 2005),而部分高温环境可使锆石U-Pb系统处于开放状态,其年龄也可能代表金伯利岩喷发的时间(Griffin et al., 2014)。然而中温条件(150~400℃)锆石裂变径迹年龄、金红石和钙钛矿U-Pb年龄可以准确限定金伯利岩的喷发时代,尤其钙钛矿测年方法应用最广,因为钙钛矿是金伯利岩浆结晶的产物,它含有高的U、Th且不受蚀变影响(Griffin et al., 2014Stanley et al., 2016)。

      辽宁瓦房店金伯利岩还没有钙钛矿U-Pb年龄的报道,山东蒙阴胜利1号金伯利岩中钙钛矿U-Pb年龄为(456±8)Ma(Dobbs et al., 1994)、457 Ma(池际尚, 1996b)、(470±4)Ma(Yang et al., 2009),低于金云母Rb-Sr等时线年龄(493 ± 10)Ma(Lu et al., 1998),(475±3)Ma(Dobbs et al., 1994)),分别与蒙阴和辽宁瓦房店金伯利岩中金云母巨晶的40Ar-39Ar等时线年龄(464.9±2.7)Ma和(463.9±6.3)Ma相吻合(张宏福等, 2007),也与瓦房店金伯利岩中金云母的Rb-Sr模式等时线年龄(461.7±4.8)Ma接近,说明两地金伯利岩群的侵位时间可能均为中奥陶世(456~470 Ma)。

      Griffin et al.(2014)对非洲南部135个金伯利岩筒开展了钙钛矿U-Pb定年,其中朱瓦能金伯利岩的钙钛矿U-Pb年龄为238 Ma,与锆石U-Pb年龄(235±7)Ma(Kinny et al., 1989)一致,稍低于金刚石中单斜辉石的40Ar-39Ar年龄(244 Ma, Burgess et al., 1992)。奥拉帕岩群中AK1金伯利岩中钙钛矿的U-Pb年龄为95~92 Ma;丹姆沙矿区的BK9为94 Ma,BK12为90 Ma,外围BK1为97 Ma;莱特拉卡内DK1金伯利岩中钙钛矿的U-Pb年龄为95~94 Ma (Griffin et al., 2014),以上年龄与奥拉帕金伯利岩锆石裂变径迹年龄((92.4 ± 6.1) Ma, (87.4 ± 5.7)Ma, Haggerty et al., 1983)一致。Gope岩群中有经济价值的Gope 25矿区的钙钛矿U-Pb为86 Ma,Gope 04为94.3 Ma,Gope136为89 Ma,Gope137为93 Ma。

      金伯利岩一般以岩筒(或岩管)、岩脉、岩席(或岩床)等形态产出(董振信, 1991Field et al., 2008Brook, 2017)。博茨瓦纳8个在采的含矿金伯利岩区均以岩筒产出。中国蒙阴和瓦房店含矿金伯利岩以岩脉为主,岩筒次之。根据岩筒垂向产状,金伯利岩岩相分为火山口相、火山通道相和根部相(Field et al., 2008路凤香, 2008),三者横截面逐渐变小,金刚石品位一般随深度增加而下降(王雪木等, 2015)。火山口相大致呈圆形盆地相,倾角一般在25~70°,直径一般50~1500 m不等,深度在150~300 m,岩石一般分为熔岩、火山碎屑金伯利岩和表层金伯利质碎屑岩,火山口相以发育明显的层状构造、弱固结、混乱岩屑流和火山碎屑物组成的沉积物为特征,这些沉积物常被湖盆沉积物覆盖。火山通道相呈圆柱状或圆锥状,筒壁垂直陡立,倾角一般在75~85°,岩壁常呈光滑状、条纹状或穿插接触。火山通道相岩石包括凝灰质金伯利岩和凝灰质金伯利角砾岩(>4 mm岩屑超过15%),常含球状火山砾、棱角状围岩和浑圆状地幔橄榄岩和榴辉岩等地幔岩捕虏体及地幔深部矿物分布于细粒基质中。根部相一般位于火山口相2~3 km以下,厚约0.5 km,沿垂向和走向上呈不规则脉状、席状或穹隆状产出,所有的岩管根部均由不同岩脉或岩脉交叉构成(Venter, 1998)。根部相岩石最接近初始金伯利岩,它由挥发分富集岩浆直接结晶而成,一般由透辉石、钙镁橄榄石、金云母、方解石和蛇纹石等5种主要矿物组成,具有斑状构造和岩浆分异和多期次侵入特征(Venter, 1998)。

      火山口相和通道相在非洲南部金伯利矿区被称为“黄土带(yellow ground)”,较容易开采,根部相又称为“蓝土带(blue ground)”,包括岩浆质或半深成矿物,常见橄榄石、金云母巨晶,以及石榴石、钛铁矿、尖晶石和铬透辉石捕掳晶,开采相对困难(路凤香, 2008)。博茨瓦纳奥拉帕和朱瓦能金伯利岩群普遍保留火山口相,而南非的金伯利、亚格斯方丹(Jagersfontein)、咖啡方丹(Koffiefontein)为火山通道相(Field et al., 2008Brook, 2017),而整个华北克拉通如辽宁瓦房店、山东蒙阴和河北涉县等金伯利岩群呈“脉多筒少”的根部相特点(池际尚, 1996b路凤香, 2008)。

      中国含矿金伯利岩筒地表出露面积普遍小于博茨瓦纳。虽然中国最大的岩筒(瓦房店42号岩筒,由3个岩管组成)地表出露面积(412500 m2,约102英亩,宋瑞祥, 2013)与博茨瓦纳最大的奥拉帕A/K1岩筒地表面积(118英亩,约477529 m2)类似,但中国金伯利岩筒地表出露面积普遍较小,如瓦房店第二大岩筒(33号)的地表出露面积约为14000 m2(3.46英亩),由两管组成的蒙阴胜利Ⅰ号岩筒地表出露面积为3988 m2(约0.99英亩)。这种较小的地表产出形态普遍与岩筒剥蚀至根部相有关,相对富集金刚石的火山口相和火山通道相剥蚀殆尽,指示中国具有寻找金刚石砂矿的巨大潜力。

      太古代克拉通陆下岩石圈地幔(SCLM)是上地幔顶部坚硬(非对流)相对较冷的亏损的岩层,具有足够低的氧化态,具有导热性,厚度一般在150~250 km,密度相对年轻地幔较小(Jacob, 2004O'Reilly et al., 2010),其下部为岩石圈-软流圈边界(LAB),与对流的温度较高相对富集的软流圈相隔。大部分SCLM在3 Ga之前就已克拉通化,高度亏损玄武质成分(如高Mg和低Ca、Al),随后多期次熔体和流体淋滤导致地球化学成分富集(O'Reilly et al., 2010)。由于太古宙克拉通SCLM底部是目前人类获得工业价值金刚石最重要的部位,前人对其进行了较为深入的研究,然而对于地幔岩石圈的成因、年龄、组成和厚度还知之甚少。金伯利岩中捕虏体能够用于示踪SCLM性质,计算周围对流地幔热、氧化还原条件、温度压力和化学状态(O'Reilly et al., 2010Aulbach et al., 2017),是研究克拉通地幔的组成、形成演化及地球深部动力学过程难得的直接研究样品之一(郑建平, 2009O'Reilly et al., 2010Shirey et al., 2013)。

      SCLM主要由地幔橄榄岩和少量榴辉岩和辉石岩等组成(Agashev et al., 2013),对应的金伯利岩中捕虏体也以地幔橄榄岩为主。其中榴辉岩捕虏体被普遍解释为俯冲洋壳的变质产物,是大洋板片俯冲至深部地幔过程的重要证据(Jacob, 2004Aulbach et al., 2017)。博茨瓦纳含金刚石金伯利岩中捕虏体类型复杂,即便同一个岩筒群,金伯利岩中捕虏体成分差异显著,如奥拉帕A/K1和朱瓦能DK2金伯利岩的捕虏体主要为榴辉岩和少量辉石岩(Aulbach et al., 2017),而距奥拉帕A/K1仅23 km的卡罗韦A/K6金伯利岩出露大量地幔橄榄岩捕虏体而榴辉岩较少,距A/K1约50 km的莱特拉卡内D/ K1捕虏体类型也不同,常见地幔橄榄岩、辉石岩、榴辉岩、云母岩、MARID地幔岩等捕虏体(Luguet et al., 2015Stanley et al., 2016Aulbach et al., 2017)。位于同一个郯庐金伯利岩带的中国蒙阴和瓦房店金伯利岩捕虏体成分亦存在差异,蒙阴常马庄和西峪金伯利岩带中常见石榴石纯橄岩、石榴石方辉橄榄岩、石榴石/尖晶石二辉橄榄岩、榴辉岩等捕虏体(宋瑞祥, 2013);辽宁瓦房店榴辉岩捕虏体极少,主要为二辉橄榄岩、纯橄岩、尖晶石方辉橄榄岩、石榴石方辉橄榄岩、金云母岩和少量麻粒岩、云母橄长岩(路凤香等, 1991郑建平, 2009宋瑞祥, 2013)。

      捕虏体中橄榄石Mg#值(Mg/(Mg+Fe))大于92代表难熔地幔,小于90代表富集地幔,在92~90之间为过渡型地幔(Zheng et al., 2001郑建平, 2009)。蒙阴和瓦房店地幔橄榄岩捕虏体/捕虏晶的橄榄石Mg#值大部分>92,仅极个别在92~90之间(郑建平等, 1999Zheng et al., 2006),类似于莱特拉卡内低Ca方辉橄榄岩和二辉橄榄岩中橄榄石的值(分别为92.9、91.5),均代表了古老的亏损难熔地幔。奥拉帕榴辉岩中单斜辉石87Sr/86Sr比值较低(0.7026~0.7046),单斜辉石87Sr/86Sr与Eu/Eu*、Mg#、Sr/Y等比值相关,可能与海底低温蚀变导致的较大范围δ18O值(单斜辉石:4.83‰~8.36‰;石榴石:5.05‰~8.18‰),这些证据指示了原岩为古洋壳岩石的榴辉岩来源于亏损地幔,并经历了分离结晶和海底低压蚀变(Aulbach et al., 2017)。郑建平等(1999)报道了蒙阴和瓦房店地幔橄榄岩全岩的86Sr/87Sr和δ18O值,其中蒙阴86Sr/87Sr比值变化较大(0.70393~0.71109),但整体低于瓦房店橄榄岩的比值(0.71216~0.71398);蒙阴地幔橄榄岩的δ18O值(4.21‰~4.45‰),明显低于瓦房店橄榄岩的值(6.68‰~7.26‰)(郑建平, 2009),也比地幔型δ18O值(5.1‰~5.9‰)稍低,这些数据反映了古生代华北克拉通SCLM可能经历了从亏损部分熔融向熔体/流体交代富集的过程。

      金伯利岩原始岩浆及其地幔源区一直存在较大争议。由于金伯利岩岩浆从产生到升至地表过程中,受到多期次的混染交代作用,使得金伯利岩中混有多源的地幔和地壳捕虏体、多源的斑晶和捕虏晶等,以至于很难确定金伯利岩原始岩浆的组成。雷雪英总结了鉴别金伯利岩是否受到混染的参数,认为未受混染的金伯利岩Al2O3<5%,SiO2= 25%~35%,Na2O/K2O<0.5,混染指数(SiO2+Al2O3+ Na2O)/(MgO + 2K2O)>1.5;而受混染的金伯利岩Al2O3>5%,SiO2>35%,混染指数则为1~1.5。并推测金伯利岩原始岩浆主量成分为SiO2=25%~32%,MgO=22%~41%,CaO=8%~20%,FeOT(全铁)=5.6%~13 %,Al2O3=1.2~3.3%,K2O=0.4%~1.9%,TiO2=0.5% ~3.4 %,CO2=4%~4%,H2O =6%~12.3%。微量元素中相容元素丰度与超基性岩相似,不相容元素丰度与碱性岩相似。

      金伯利岩的源区争议涉及上地幔(SCLM、软流圈地幔)(Becker et al., 2006;张宏福等, 2007Tappe et al., 2012)、下地幔或核幔边界(Wu et al., 2010)、地幔柱(Chalapathi Rao et al., .2016)等观点。普遍认为放射性同位素是示踪金伯利岩地幔源区的重要方法(Wu et al., 2010)。传统研究大多基于全岩同位素成分,例如Smith(1983)根据全岩Pb-Sr-Nd同位素和矿物地球化学特征,结合金伯利岩的侵位年龄,将非洲南部金伯利岩分为Ⅰ型(大部分114~80 Ma,低Sr、高Nd)和Ⅱ型(150~114 Ma,高Sr、低Nd)金伯利岩,分别对应于玄武岩型和云母型。然而由于金伯利岩包含不同类型的壳源和幔源捕虏体以及受到后期蚀变影响,该分类中全岩同位素值可能代表非真实的金伯利岩浆同位素值(Yang et al., 2009),但部分作者认为Nd同位素可用于追溯地幔源区(Tappe et al., 2012),而钙钛矿作为初始岩浆结晶矿物,其Sr-Nd-Hf放射性同位素是示踪金伯利岩源区另一个有效方法(Wu et al., 2010Tappe et al., 2012Griffin et al., 2014)。

      Griffin et al.(2014)对非洲南部110个地区的135个金伯利岩筒开展了钙钛矿的原位U-Pb定年和Sr、Nd同位素分析,获得了1800~1600 Ma、1200~1000 Ma、800~500 Ma、180~170 Ma和130~50 Ma五个阶段年龄峰,指示岩石圈交代岩富集事件早在1.8 Ga已经开始,而100~80 Ma是最重要的岩浆事件。非洲南部包括奥拉帕岩群在内的(90±10) Ma的Ⅰ型金伯利岩喷发之前至少经历了大约180 Ma和130~120 Ma的Ⅱ型金伯利岩的局部缓慢的岩浆活动,高度富集的SCLM是Ⅱ型金伯利岩的地幔源区(高87Sr/86Sr,低εNd值,高Rb和LREE/HREE,源自MARID和金云母-钾-锰闪石交代的地幔橄榄岩),这种富集SCLM在1.8 Ga既已存在于岩石圈深部,在热的软流圈物质或地幔柱上升过程中岩石圈发生减薄,导致富集部分在110 Ma之前被完全消融,以至于年轻的Ⅰ型金伯利岩具有低87Sr/86Sr和高εNd (0~+5)值,源区来自更浅的亏损SCLM,且几乎未与老的富集SCLM发生交代反应(Griffin et al., 2014)。Wu et al.(2010)也报道了非洲南部Ⅰ型金伯利岩中的钙钛矿类似的εNd值(+1.2~+3.1),其远远低于亏损软流圈地幔值(+10, Salters et al., 2004)),结合中等的87Sr/86Sr比值(0.704~0.705),指示未分异或弱亏损的地幔源区,暗示Ⅰ型金伯利岩浆可能源自下地幔(Wu et al., 2010)。

      山东蒙阴金伯利岩全岩87Sr/86Sr比值普遍较高(0.70548~0.70668,池际尚,1996b),0.704489~0.709783,张宏福等, 2007),被解释为源自富集地幔,可能为华北SCLM受到大洋板片俯冲流体交代的产物(张宏福等, 2007),然而上述全岩87Sr/86Sr明显高于钙钛矿值(0.70371±12,Yang et al., 2009)说明受到蚀变混染的影响。蒙阴金伯利岩中钙钛矿的初始87Sr/86Sr比值(0.70371±12)较低,并具有较窄范围的εNd值(0.13±0.22)(Yang et al., 2009),21个全岩147Sm/144Nd比值为0.512072~0.512408(池际尚, 1996b),小于球粒陨石平均值(0.512638),以上特征类似于非洲南部Ⅰ型金伯利岩。结合早古生代华北陆块东部从辽宁至山东无大洋俯冲记录,金伯利岩中发现碳化钨、自然铁等还原矿物(池际尚等, 1996b),指示蒙阴金伯利岩浆源自弱亏损地幔或初始地幔源区,可能与深部地幔物质上涌或地幔柱上侵至华北岩石圈底部有关(Yang et al., 2009)。复县金伯利岩类似于Ⅰ型,目前还没有钙钛矿同位素数据的报道,然而根据其与蒙阴金伯利岩分布于郯庐带,具有类似的根部相、岩石组合及侵位时代,笔者认为两者具有一致的地幔源区特征。

      华北克拉通面积约200万km2,保存3.8 Ga古老残余陆壳,其基底主要由太古宙—古元古代TTG片麻岩、麻粒岩和斜长角闪岩等高级变质组成,中元古代(约1.8 Ga)东、西陆块碰撞拼合,完成终极克拉通化并进入稳定地台发展阶段(Zhao et al., 2005),上部覆盖巨厚中、新元古代至早古生代沉积盖层。然而自465~470 Ma至中石炭世,发生明显的沉积间断,该间断不仅与金伯利岩的侵位构造事件一致,也与蒙古洋向南俯冲到华北克拉通之下,形成一系列俯冲增生和碰撞造山作用,形成华北克拉通北部巨大的中亚造山带等构造事件相吻合。中生代早期,华北克拉通岩石圈地幔被俯冲流体/熔体交代,导致克拉通北部岩石圈地幔减薄破坏,并在T3发生广泛伸展活动(Niu et al., 2016)。

      华北克拉通岩石圈破坏的时间为中生代,峰期在125 Ma左右(朱日祥等, 2011),而中国东部含金刚石金伯利岩的喷发时间在456~470 Ma,说明在中生代华北岩石圈破坏之前喷发的“蒙阴—瓦房店”型且源自初始地幔的金伯利岩均可能含矿。考虑到华北克拉通自古生代以来受到古亚洲洋(或古蒙古—鄂霍次洋)、古/现代太平洋的俯冲作用,并经历了470~320 Ma长期剥蚀及随后晚石炭至早三叠世的海陆交互相稳定沉积(朱日祥等, 2011),暗示中生代之前喷发的“蒙阴—瓦房店”型含矿金伯利岩可能被晚石炭至早三叠世沉积岩覆盖或者部分被剥蚀。由此推测,中国金伯利岩型和冲积型金刚石矿床均有找矿潜力。

      非洲南部白垩纪金伯利岩侵位的主构造背景为大西洋打开,其中(90±10) Ma金伯利岩的喷发与非洲和南美陆块的板块的裂离有关,地幔柱活动和古大洋岩石圈俯冲是大陆裂解可能的驱动力。金伯利岩一般沿线性和弧形破碎带和深大断裂分布,卡拉哈日克拉通内基底断裂和破碎带走向南西-北东,长约1600 km,宽50~90 km,成因与大西洋洋中脊右旋转换有关。Griffin et al.(2014)统计了非洲南部金伯利岩中钙钛矿的εSr值和εNd值与Vs波在150 km深度的关系发现金伯利岩的地球化学富集特征与低Vs波带密切重合,指示延伸到岩石圈底部的长期构造活动带,如主要区域断裂、克拉通边界剪切带和破碎带,控制金伯利岩的交代作用和侵位。该认识与中国学者认为金伯利岩的产出与深切地幔的断裂系密切相关一致,它是金伯利岩型金刚石重要的控岩要素(黄蕴慧等, 1992池际尚等, 1996b)。例如南非金伯利岩群沿北东向岩石圈尺度的芬特斯多普(Ventersdorp)深大断裂排列,中国蒙阴和瓦房店含矿金伯利岩群均沿郯庐深大断裂分布。

      金刚石能够抵御蚀变和再结晶对包裹体的干扰,能够保存壳源、幔源甚至核幔边界的物质组分,记录了金刚石的形成环境、形成时代和形成机制等重要信息,是人类认识地球深部过程,探讨壳-幔之间碳等物质循环最直接的样品(连东洋等, 2018)。金伯利岩型金刚石中包裹体丰富,不仅有源自岩石圈地幔的包裹体,还有少量可能来自于软流圈、地幔过渡带、下地幔矿物和流体(Harte, 2010Stachel et al., 2015Nestola et al., 2018)。因为克拉通结构控制石墨-金刚石形成的压力温度条件,金刚石一般形成在克拉通150 km以下相对冷却、还原和低密度岩石圈根域,只有在太古宙—早元古代克拉通相对低的温度(<1200℃)下,才能保存较大范围的金刚石稳定“生存”区。统计显示,含包裹体的金伯利岩型金刚石中约99%金刚石形成于岩石圈地幔,该类包裹体主要为橄榄石、辉石、石榴石、尖晶石以及硫化物等矿物(Shirey et al., 2013连东洋等, 2018)。约1%含包裹体的金刚石来自于软流圈、地幔过渡带、下地幔甚至核幔边界,如博茨瓦纳卡罗韦金伯利岩型金刚石中发现超硅石榴石(majorite)(Motsamai et al., 2018),它是鉴别来自软流圈和地幔过渡带金刚石的主要指示矿物(Stachel, 2001Tappert et al., 2005);南非、巴西、加拿大金伯利岩型金刚石中有来源于下地幔标志性矿物组合:瓦硅钙钡石(walstromite)和斜硅钙石(larnite)(Smith et al., 2018)、铁方镁石+镁硅钙钛矿,以及似榴辉岩型矿物组合、斯石英、碳酸岩型矿物组合等(Stachel, 2001Bulanova et al., 2010Kaminsky, 2012);南非Monastery和委内瑞拉Guaniamo金伯利岩型金刚石发现了可能来自核幔边界D”层的矿物组合:(1)方铁矿(wüstite)+方镁石(periclase)和(2)铁的碳化物(ironcarbide)+自然铁矿物组合(Kaminsky, 2012连东洋等, 2018)。中国山东蒙阴和辽宁复县的金刚石具有相似的岩石圈矿物包裹体组合,主要包含橄榄石、铬尖晶石、透辉石、石榴石、石盐、自然铁、硫化物等矿物,主体为P型,仅有少数矿物显示E型包裹体的化学特征(Wang et al., 1998殷莉, 2008)。除P、E型包裹体外,陆琦等(2011, 2012)在瓦房店金刚石中发现有碳化钛(TiC)、超硅石榴石和钙钛矿等高压-超高压矿物,指示一部分金刚石可能来源于华北克拉通的软流圈甚至地幔过渡带。

      硫化物包裹体常用于示踪金刚石生长、地幔熔融、俯冲物质加入和地球内部物质循环等过程(Kemppinen et al., 2018)。金刚石中贱金属硫化物(Base metal sulphide,BMS)包裹体是认识地幔中亲铜和亲铁元素组成的独特样品,它比地幔捕虏体、造山带和蛇绿岩地幔橄榄岩能更好的认识地幔的微量元素、Re-Os和S同位素组成。P型和E型金刚石的BMS包裹体成分明显不同(如前者Ni含量明显高于后者),P型金刚石具有较高的PGE含量(尤其Os, Ir, Ru)和Te、As元素含量(Mcdonald et al., 2017)。金刚石包裹体如硫化物Re-Os同位素年龄和硅酸盐(如石榴石和单斜辉石)Sm-Nd等时线年龄可以确定金刚石形成的时代,结果显示同一个矿区普遍存在多期年龄的金刚石(Gress et al., 2018),例如朱瓦能E型金刚石硅酸盐Sm-Nd等时线上下交代年龄分别为(1540±20) Ma和1 Ga,硫化物的ReOs等时线上下交点年龄为2.9 Ga和1.5 Ga (Thomassot et al., 2009)。此外,石榴石和单斜辉石包裹体的Rb-Sr、Sm-Nd、氧同位素组成也能够鉴别俯冲沉积物和亏损或富集地幔源区(Timmerman et al., 2017),比如博茨瓦纳奥拉帕和莱特拉卡内具有多期生长环带的单晶金刚石(δ13C值为-5.0‰~-33.6‰)与地幔中多源含碳流体有关,其辉石和石榴石包裹体的87Sr/86Sr比值(0.7033~0.7097)亦反映了俯冲的生物沉积物、俯冲的富集碳沉积物和亏损地幔三个源区端元(Timmerman et al., 2017)。

      SCLM内的钻石生长与大陆裂解和聚合及地幔柱活动等大规模岩浆活动之间存在潜在耦合关系(Richardson et al., 2004),这些构造热事件过程中结晶的金刚石常含有C-H-O-S富集的超临界流体和硅酸盐熔体以及少量记录构造环境变化的矿物(Koornneef et al., 2017)。如南非威尼西亚(Venetia) 2.95 Ga方辉橄榄岩型金刚石具有低CaO、低HREE和“正弦型”REE配分模式的石榴石包裹体,其成因与低温很低Ca、Fe、Ti和Al流体为主的交代作用有关,指示一种异常亏损的地幔残留特征。此外该类金刚石的δ13C值(-7.1±2.7‰)与石榴石CaO含量正相关,而含钙石榴石的金刚石一般具有较低的碳同位素含量,反之亦然。这种石榴石没有典型俯冲特征(Nb/La>1),因此该金刚石不可能与再循环有机碳有关,而形成于相对少量的超临界还原流体(CH4)发生Raleigh型分离结晶导致同位素分离至更低的δ13C值,这种P型金刚石记录了津巴布韦克拉通南缘裂解过程中以温度相对较低的流体促使太古宙金刚石结晶。(1.15±0.11) Ga金刚石的石榴石包裹体的CaO、FeO和REE含量较高,指示中等富集C-H-O-LREE的流体交代了亏损的SCLM。这种类型的交代作用通常与源自软流圈高温(T>1190℃)硅酸盐熔体交代作用形成二辉橄榄岩型石榴石有关(Koornneef et al., 2017)。结合金刚石的δ13C平均值为(-4.6‰)类似于软流圈值,该类金刚石生长记录了约1.1 Ga乌肯多(Umkondo)大火成岩省有关的流体交代津巴布韦克拉通的亏损地幔(Koornneef et al., 2017)。类似的实例还有朱瓦能3.2 Ga(硫化物的Re-Os等时线年龄,下同)方辉橄榄岩型金刚石的形成与西部卡普瓦尔陆核的形成有关,2.9 Ga和1.5 Ga E型金刚石生长分别与东、西卡普瓦尔克拉通聚合以及克拉通的后期改造相一致(Richardson et al., 2004)。

      金刚石的碳、氮同位素和氮含量可以为金刚石物质来源、生长过程和成因提供重要的信息(Shirey et al., 2013Cartigny et al., 2014连东洋等, 2018)。如奥拉帕A/K1矿区的黄色具有环带构造的E型金刚石,核部为黄色,边部为无色。单斜辉石-石榴石矿物对包裹体获得金刚石形成的压力和温度分别为5 Gpa和1080~1160℃(Deines et al., 2004Timmerman et al., 2018)。氮含量较低(普遍<50×10-6),为未聚合的单原子氮;核部δ13C值介于-2.9‰~-10.6‰之间,边部δ13C值位于-6.6‰~-10.5‰之间(Timmerman et al., 2018),核边不同的碳同位素值指示金刚石经历了多阶段生长,边部较窄的δ13C值可能是亏损地幔和俯冲沉积物混合的结果(Timmerman et al., 2018)。

      Shirey et al.(2013)Cartigny et al.(2014)连东洋等(2018)等国内外学者系统综述了不同类型金刚石的碳、氮同位素及氮含量,发现全球金刚石的δ13C值范围极广(-41‰和+5‰),接近于沉积岩的范围,其中超过72%的金刚石的δ13C值介于-8‰~-2‰,与源自下地幔金刚石(-8.5‰~-0.5‰)、纤维型/复合型(fibrous/coated, -8.1‰和-4.9‰)的δ13C值基本重合,并在(-5±1)‰的位置出现峰值,与洋中脊玄武岩、洋岛、碳酸岩、金伯利岩等幔源岩石δ13C值类似(Cartigny et al., 2014)。他们还发现不同生长环境下金刚石的δ13C值差异显著:P型金刚石的δ13C值介于-26.4‰~+0.2‰,而E型金刚石的δ13C值介于-41.3‰~+2.7‰,过渡带和下地幔的E型金刚石的δ13C值<-10‰,形成于超高压俯冲变质岩中金刚石的δ13C值介于-30‰~-3‰,发现于中非和巴西沉积物中的黑金刚石(carbonado,由金刚石、石墨和无定形碳组成的多晶金刚石)的δ13C值介于-32‰和-5‰(Shirey et al., 2013),不同类型陨石中的金刚石δ13C值介于-38.8‰~-32.5‰(Russell et al., 1996),蛇绿岩型金刚石的δ13C值介于-28.4‰~-18.0‰ (Lian et al., 2018连东洋等, 2018)。俯冲带沉积物和地表碳酸盐矿物的δ13C值平均值分别为-25‰和0 (Cartigny et al., 2014)。此外相同矿体中相同成因的金刚石,含硫化物包裹体的与含硅酸盐的金刚石碳同位素亦不同(Cartigny et al., 2014)。一般认为P型、地幔过渡带和下地幔金刚石的碳来自地幔,而E型、超高压变质型和蛇绿岩型金刚石的碳源自俯冲沉积物。但不同类型金刚石的碳源鉴别还存在着争论,比如Cartigny et al.(2001)认为在开放系统瑞利(Rayleigh)条件下或铁碳化物还原条件下,碳同位素在地幔中会发生显著分馏,那么仅根据金刚石的碳同位素不一定能确定其碳的来源(Cartigny et al., 2001)。此外,碳作为太阳系中含量第四高的元素,碳同位素的分布,特别是在地核中的分布,在很大程度上仍是未知的(Satish-Kumar et al., 2011)。除铁、镍元素外,地核还含有地球上大部分的铬和钒、硫、磷和碳元素(Mcdonough, 2013),那么在地幔柱活动异常活跃的地球早期,地核可能是形成金刚石的重要碳源之一。

      俯冲带是壳幔碳循环的关键部位,地表有机碳的俯冲不仅影响了地球深部金刚石的生长,还与古元古代地球大气中氧含量升高密切相关(Duncan et al., 2017)。许多古老金刚石及其榴辉质包裹体的碳同位素组分显示板片俯冲作用最早开始于3.0~3.5 Ga,有机碳俯冲至深部地幔是导致2.2~2.5 Ga大气氧含量上升的主要原因(Duncan et al., 2017)。Duncan et al.(2017)以碳酸盐、地幔橄榄岩、有机碳(δ13C平均值分别为0、-5‰和-25‰)作为碳源的三个端元,并假设初始源区只有少量δ13C分离结晶,模拟得出δ13C值小于-11‰的金刚石可能来自循环的有机物,认为碳能够在含水流体、碎屑沉积物的部分熔体中,或在板片的玄武岩/榴辉岩中从俯冲板片中释放出来。CO2含量随着压力的增加而增加,而与温度没有明显相关性。即使在非常高的温度下,由有机碳或海底碳酸盐还原而产生的石墨碳大部分仍保留在俯冲板块中,进而有效的赋存于地幔中,从而提供了一种从地表储集层长期去除和储存还原轻质碳的机制,即通过CO2+H2O=CH2O+O2途径,有效地从大气中吸收碳的机制,建立全球系统的俯冲循环和地幔中还原性碳的长期封存。地球早期(>2.5 Ga)的俯冲一定导致了石墨化碳优先于碳酸盐的俯冲,太古宙的热俯冲过程可能没有氧化能力,甚至可能在局部降低了地幔氧含量(Duncan et al., 2017)。这些过程不仅为金刚石结晶提供了碳源,还改变了金刚石生长的氧化还原条件,进而影响金刚石中氮元素的聚集。

      氮是天然金刚石中最常见的杂质元素,不同类型金刚石的氮含量差异明显(连东洋等, 2018):超高压变质型金刚石氮含量范围最大,介于0~9600×10-6 (连东洋等, 2018);纤维型/复合型金刚石具有较高的氮含量(600×10-6~1800×10-6)峰值在900×10-6~1200×10-6;蛇绿岩型金刚石的氮含量数据较少,介于7×10-6~541×10-6 (Lian et al., 2018);P型和E型金刚石的氮含量类似,介于0~2100×10-6,峰值在0~300×10-6。不同类型金刚石的氮同位素也存在差异:纤维型/复合型金刚石的δ15N值在-12‰~0;而P型和E型金刚石分别在-24‰~+12‰和-12‰~+18‰;超高压变质型金刚石δ15N值均为正值,与变质沉积物具有相似的同位素组成;蛇绿岩型金刚石δ15N值为-5.6‰~+ 28.6‰。总的来说,幔源型金刚石具有类似的氮同位素组成,明显区别于超高压变质型。

      Mikhail et al.(2016)认为虽然地幔是氮的最大储库,但人类对氮在地幔的地球化学性质的认识是有偏差的,因为>90%氮的数据来自金刚石,而金刚石能在多大程度上记录流体/熔体氮的特征还不清楚。氮在金刚石中的相容性一直存在着争论,最初作为单原子氮进入金刚石晶格中的氮是高度不稳定的,以分子形式出现的氮(如N2或NH4+)在金刚石晶格中也是不相容的。那么不同种类氮在钻石形成过程中如何分解产生瞬时单原子氮的,是决定金刚石中氮丰度的因素,即氮在金刚石中的丰度主要受介质中氮参与金刚石生长的程度决定,意味着介质中氮的浓度对金刚石的氮含量不是主要影响因素,例如当富氮介质并没有发生释放单原子氮的反应时,无氮(Ⅱ型)金刚石可以从富氮的介质中生长,相反含氮(Ⅰ型)金刚石可以从含氮量较低的流体中生长(Mikhail et al., 2016)。总的来说,促使介质发生释放氮的反应条件主要包括氮的摩尔丰度、压力、温度、生长速率和氧化还原、酸碱度等。金刚石的平均氮丰度普遍随深度的增加而降低,这可能反映了随着深度(温度、压力)的增加,氮在其他矿物相中随着深度的增加而更加相容。快速不平衡生长过程中常形成富氮金刚石,其N/C比接近于流体或熔体;相比之下,在接近平衡缓慢生长条件下,无论生长介质中氮的含量如何,只可能形成Ⅱ型金刚石。金刚石中氮含量最大的影响因素是氧逸度和酸碱度。在还原条件下硫化物普遍存在可能会抑制氮化物的形成。当氧化/还原和酸/碱条件发生变化时,地幔中N2、NH4+、CO2或CH4可能分解成单原子氮进入金刚石晶格中。

      寻找金伯利岩的指示矿物一般为抗物理化学风化强、比重大的矿物,主要包括锆石、铬尖晶石、钛铁矿、镁铝榴石、铬透辉石以及少量橄榄石、金云母、顽火辉石/古铜辉石、钾锰闪石/镁红闪石等,其中石榴石、钛铁矿和铬尖晶石的二次物理磨损或化学风化形态特征能够进一步评估距离初始源区的距离,是最有效的指示矿物(Venter, 1998)。董振信(1991)包淑华(2013)等学者系统总结了金伯利岩产状、结构构造和矿物组成与含矿性的关系,认为岩管状金伯利岩金刚石含矿性较好,而岩脉状岩体很少含金刚石;具斑状结构、块状构造的金伯利岩含金刚石较富,而微斑状及碎屑结构的金伯利岩含金刚石较贫;粗斑或多斑状结构的金伯利岩和富含高温、高压矿物的岩体,一般金刚石含矿物性强,否则含矿性较差;在碎屑结构的金伯利岩中,含同源碎屑多的金伯利岩金刚石较富,含异源碎屑较多的含金刚石较贫;斑晶中橄榄石(或假象)颗粒越粗、总含量越高,金刚石就较富,反之就较贫;含铬尖晶石、铬镁铝榴石、硅灰石越多,金刚石就越富;金云母不是影响金刚石含量的重要因素;锐钛矿的含量则与金刚石含量呈负相关,这与金刚石中元素钛和金刚石含量关系呈反比的规律是一致的(包淑华, 2013)。

      金红石能够直接确定金伯利岩的侵位年龄和物源研究,揭示幔源碱性岩浆岩的存在和含矿性(Malkovets et al., 2016)。金红石是金伯利岩中常见的重矿物,其寄主岩石和赋存矿物包括广泛的地壳岩石、地幔捕虏晶、捕虏体(榴辉岩、辉石岩、MARID组合、交代地幔橄榄岩、金红石巨晶、金刚石包裹体、金刚石共生)等。金红石中Pb封闭温度600~640℃,其U-Pb年龄可以指示金伯利岩的侵位时代。克拉通地幔交代岩石中金红石的Cr2O3含量变化较大(0.64%~9.75%),地壳岩石中金红石的Cr2O3含量最低(<0.4%),克拉通边缘(off-craton)碱性玄武岩中交代尖晶石地幔橄榄岩捕虏体中金红石的Cr2O3含量位于两者之间。幔源碱性岩浆岩(金伯利岩、碳酸岩、钾镁煌斑岩、煌斑岩等)中富铬金红石通常与交代的幔源原生岩有关。通常,Cr含量<1000×10-6指示源自变泥质源区,而Cr含量>1000× 10-6指示源自基性源区,而源自交代地幔捕虏体和捕虏晶的金红石的Cr2O3可达9.75%。高Cr含量的金红石被解释为继承自原本亏损的克拉通地幔,而Ti、Fe、Nb等不相容元素则被交代流体/熔体带入。金红石中Cr2O3含量高于1.7%被归为地幔克拉通,而小于1.7%可能源自克拉通边缘和克拉通地壳和地幔源区。因此,富Cr金红石可作为寻找含金刚石金伯利岩一种潜在的指示矿物(Malkovets et al., 2016)。

      地球物理异常是覆盖区寻找金伯利岩的重要标志,由于金伯利岩体所处地质环境的特殊性,各物化探方法的效果是不一致的,重力、磁法、常规电法、大地电磁法等综合地球物理方法对寻找金刚石原生矿有互补作用(张建太, 2017)。博茨瓦纳早期金伯利岩型金刚石矿的发现绝大多数是通过土壤取样分析石榴石、镁钛铁矿、单斜辉石等指示矿物获得的(De Wit, 2018)。最近60年该国使用航空磁测等地球物理方法寻找新的金伯利岩体显著增加,然而在过去的15年里,博茨瓦纳发现金伯利岩的几率亦显著下降,这是由于对土壤取样的勘探工作显著减少的结果。因此,地球物理方法应与传统找矿方法结合,尤其在弱磁性和无磁性地区寻找金伯利岩体时,传统的土壤取样分析指示矿物仍是重要方法。

      (1) 非洲南部和中国有经济价值的金刚石原生岩石均分布在克拉通及其周缘古元古代构造带中。博茨瓦纳含矿原生岩石仅为金伯利岩,与卡普瓦尔和津巴布韦克拉通密切相关。该国目前共发现12个金伯利岩区,其中4个岩区的8个矿在开采,分别是奥拉帕(A/K1)、莱特拉卡内(D/K1)、朱瓦能(DK2)、丹姆沙(B/K9、B/K12)、勒雷拉(K002~K006)、B/K11、卡罗韦(A/K6)和Ghaghoo。中国含矿原生岩石分布复杂,金伯利岩主要分布在华北克拉通的郯庐、华北中央和华北北缘金伯利岩带,仅蒙阴和瓦房店具有工业价值,分布于郯庐金伯利岩带中;钾镁煌斑岩主要出露在华南克拉通,重点分布于江南和华南北缘钾镁煌斑岩带中;塔里木陆块研究程度较低,仅皮山县发现少量钾镁煌斑岩脉。郯庐金伯利岩带、江南钾镁煌斑岩带和塔里木陆块是中国含金刚石原生岩石重要的找矿区。

      (2) 根据来源深度,金伯利岩型金刚石被分为岩石圈型和超深型,其中99%的含包裹体金刚石为岩石圈型,分布于克拉通下150~250 km范围内。超深型可以来自于软流圈、地幔过渡带、下地幔甚至核幔边界。奥拉帕A/K1和朱瓦能DK2金伯利岩体是世界上为数不多的主要产出榴辉岩捕虏体和E型金刚石的岩筒之一,而同位于奥拉帕岩群的莱特拉卡内、丹姆沙和卡罗韦岩体与中国郯庐带的金伯利岩体类似,均主要产出地幔橄榄岩捕虏体以及P型和E型金刚石。

      (3) 钙钛矿是金伯利岩浆结晶的产物,其U-Pb定年和Sr、Nd和Hf同位素可以有效限定金伯利岩的时代和成因。238 Ma朱瓦能、86~97 Ma奥拉帕岩筒群的金伯利岩与456~470 Ma蒙阴和瓦房店金伯利岩均属于Ⅰ型金伯利岩,具有低87Sr/86Sr(0.703~0.705)和中等εNd(t)值(-0.09~+5),指示金伯利岩浆源自弱亏损深部地幔或初始地幔源区。博茨瓦纳金伯利岩形成于三叠—白垩纪,相比之下,中国含矿金伯利岩形成较早,被中奥陶世以来地层覆盖,因此隐伏矿床是中国重要的找矿方向。

      (4) 金伯利岩侵位受克拉通相关的深大断裂控制,因此克拉通内部和周缘深大断裂带是寻找含矿金伯利岩重要的控岩构造要素。石榴石、镁钛铁矿、铬透辉石、铬尖晶石和铬金红石等是寻找含金刚石金伯利岩重要的指示矿物。航磁等地球物理找矿方法需与土壤取样分析指示矿物方法相结合才能取得更好的找矿效果。

      (5) 一般来说从金伯利岩筒的火山口相、火山通道相至根部相金刚石的品位逐渐降低。博茨瓦纳金伯利岩筒绝大部分为火山口相,而中国蒙阴和瓦房店等金伯利岩筒均为根部相,表明后者经历了强烈的剥蚀,这可能是中国金伯利岩型金刚石资源“匮乏”的重要原因之一。金伯利岩筒的强烈剥蚀作用对冲积型矿床的形成有利,暗示中国具有冲积型金刚石矿床的巨大找矿潜力。

      (6) 金伯利岩和金刚石是人类认识地球深部(地幔和地核)性质和组成、壳幔和核幔物质循环以及古环境变化等科学问题的重要研究对象。博茨瓦纳和中国金刚石均记录了至少3~4次生长过程,然而金伯利岩初始岩浆的组成、来源和演化等,金刚石碳源和成因(流体结晶还是固态相变),氮硼等元素进入金刚石的机制,多源包裹体对金刚石生长的影响,超深型金刚石中核幔物质的来源和演化等问题存在争议或没有解决。中国在上述方面的研究起步较晚,目前还处于学步和跟跑阶段,亟需开展进一步研究。

      注释

      http://www.debswana.com/

      https://en.wikipedia.org/wiki/Letlhakane_diamond_minewww.debswana.com/Operations/Pages/Letlhakane-Mine.aspx

      http://www.debswana.com/Operations/Pages/Damtshaa-Mine.aspx

      https://www.lucaradiamond.com/

      https://www.lucaradiamond.com/mine-operations-and-exploration/karowe-mine-summary

      http://www.debswana.com/operations/pages/jwaneng-mine.aspx

      1➊刘向君. 2008. 川西地区工程地质特征研究[R]. 四川: 中石化西南油气分公司工程技术研究院.
    • 图  1   新场—合兴场—丰谷地区TX22断层平面分布图及地层特征

      a—断层平面分布图;b—须二段砂组划分剖面图(Xinshen 1井)

      Figure  1.   TX22 fault plane distribution and stratigraphic characteristics in Xinchang–Hexingchang–Fenggu area

      a–Plane distribution diagram of faults; b–Section diagram of sand group division of the Xu 2 member (Xinshen 1 well)

      图  2   新场—合兴场—丰谷地区南西—北东向地震剖面解释

      Figure  2.   Interpretation of southwest−northeast seismic profile in Xinchang–Hexingchang–Fenggu area

      图  3   声发射岩心取样示意图及XS1井声发射实验试件

      Figure  3.   Schematic diagram of acoustic emission test core sampling and XS1 well test specimen

      图  4   差应变实验的试件及应变片贴法示意图

      Figure  4.   Schematic diagram of specimen and strain gauge sticking method for differential strain test

      图  5   XS1井4487 m,1−S4号样品累计能量对数、振铃计数与到达时间关系图

      Figure  5.   The relationship between cumulative energy logarithm, ringing count and arrival time of No.1−S4 sample at 4487 m in Well XS1

      图  6   CX560井差应变结果图

      Figure  6.   CX560 Well differential strain results

      图  7   多种方法计算单点最小水平主应力值对比图

      Figure  7.   Comparison of various methods to calculate the minimum horizontal principal stress value of single point

      图  8   诱导缝形态与差应力系数之间的关系图

      Figure  8.   Diagram of the relationship between induced fracture morphology and differential stress coefficient

      图  9   鞍形诱导缝岩心示意图及X10井岩样照片

      Figure  9.   Saddle induced fracture core diagram and X10 well rock sample photos

      图  10   不同崩落宽度下水平应力量值示意图

      a—CG561井,4810.00 m,崩落宽度20°; b—X10井,4905.20 m,崩落宽度40°

      Figure  10.   Schematic diagram of horizontal stress value under different caving widths

      a–CG561 well, 4810.00 m, caving width of 20°; b–X10 well, 4905.20 m, caving width of 40°

      图  11   井壁崩落宽度与水平两向差应力关系图

      Figure  11.   Relationship between the width of the wellbore caving and the horizontal two-dimension differential stress

      图  12   X11井须二段各小层岩石力学与地应力综合剖面

      Figure  12.   Comprehensive section of rock mechanics and in-situ stress of each sublayer of the second member of Xujiahe Formation in Well X11

      图  13   研究区须二段典型地应力剖面类型

      Figure  13.   Typical type of geostress profile in the second segment of Xujiang Formation in the study area

      图  14   研究区TX22小层井上地应力剖面类型主要分布图

      Figure  14.   Main distribution of surface stress profile types of TX22 sublayer in the study area

      图  15   不同储隔层应力差模拟水力裂缝扩展形态图

      a—隔层应力为97 MPa,储层应力为96 MPa对裂缝形态影响;b—隔层应力为97 MPa,储层应力为90 MPa对裂缝形态影响

      Figure  15.   Simulated hydraulic fracture propagation pattern of different reservoir and interlayer stress difference

      a–The effect of interlayer stress of 97 MPa and reservoir stress of 96 MPa on fracture morphology; b–The effect of interlayer stress of 97 MPa and reservoir stress of 90 MPa on fracture morphology

      表  1   水力压裂法计算的地应力大小的数据结果

      Table  1   Results of the magnitude of in-situ stress calculated using hydraulic fracturing

      井号层位测试段中深/m最大主应力/MPa最小主应力/MPa垂向主应力/MPa井所在构造位置
      CH127须二4581.50141.35101.35114.61F6上盘
      GM2须二14713.00146.99122.99116.51F23上盘
      GM2须二24776.00150.31122.31116.19F23上盘
      GM2须二34845.00149.87109.87118.82F23上盘
      GM3须二34920.00157.31125.31121.15H7上盘
      GM4须二44882.00153.06113.06113.98F44上盘
      X10须二24717.00138.68103.00117.99F24上盘
      X10须二65042.00153.00112.40125.92F24上盘
      X10-2须二24732.00143.9683.96111.41新场构造七郎庙高点
      X10-2须二44880.00154.1894.18114.89新场构造七郎庙高点
      X10-2须二44851.50154.2294.22114.22新场构造七郎庙高点
      X11须二44917.00145.00113.00121.45新场构造七郎庙高点北翼
      X209须二44868.00153.66104.46124.16新场构造七郎庙高点
      X5须二34880.00155.05107.05119.68F4上盘
      X8-1H须二24880.00141.92101.52117.29F3下盘
      XC12须二85250.13162.31102.31128.76孝泉构造
      XS1须二14496.00127.3787.37112.03F8上盘
      XS101须二14598.00144.79112.79114.12F9-1上盘
      XS101须二24807.00149.07109.07120.25F9-1上盘
      XS101须二24786.00150.43110.43118.79F9-1上盘
      XS101须二24709.00144.63104.63115.49F9-1上盘
      XS101须二24634.00143.50103.50112.29F9-1上盘
      下载: 导出CSV

      表  2   孝泉—丰谷地区须二段声发射实验测试的三向应力值的测试结果数据

      Table  2   Test results of the three−direction stress value of the acoustic emission test of the second member of the Xuxujiahe Formation in the Xiaoquan–Fenggu area

      井号层位井深/mKaiser点应力值/MPa最大主应力/MPa最小主应力/MPa垂向主应力/MPa
      45°90°垂直
      XS1须二14487.00101.3382.1088.62101.29118.7690.04110.71
      X11须二24762.24116.1692.7499.78119.56135.26100.68129.56
      CG561须二44943.80112.7489.5197.28109.74132.7198.07120.12
      XC12须二44766.83111.6795.32106.53114.13133.13105.09124.14
      X501须二65270.70120.8799.48110.79117.80144.01109.79128.87
      下载: 导出CSV

      表  3   孝泉—丰谷地区须二段差应变实验测试的三向应力值的测试结果数据

      Table  3   Test results of triaxial stress value of the second member of Xujiahe Formation in Xiaoquan–Fenggu area

      井号 层位 深度/m 三向主应力/MPa
      最大主应力 最小主应力 垂向主应力
      CX560 须二2 4810.97 141.00 102.00 125.00
      CH127 须二1 4566.14 135.09 100.65 111.87
      CH127 须二2 4639.29 137.60 104.72 113.66
      X10 须二4 4882.06 148.86 107.93 119.61
      X10 须二4 4884.04 140.75 104.76 119.66
      下载: 导出CSV

      表  4   研究区样品地应力大小计算实例

      Table  4   Calculation example of in-situ stress of samples in the study area

      样品编号 深度/m 诱导缝形态 差应力系数 端面尺寸/mm Δd/μm 弹性模量/GPa 泊松比 最大水平主应力/MPa 最小水平主应力/MPa
      dmax dmin
      X10井 4928.92 马鞍形 0.25~0.30 69.207 69.157 50 45.359 0.182 120.23~138.72 92.48~110.98
      下载: 导出CSV

      表  5   不同构造变形单元内不同层位构造应力系数反算结果

      Table  5   Back−calculation results of tectonic stress coefficient of different layers in different tectonic deformation units

      层位构造应力系数构造变形单元
      新场合兴场丰谷
      TX21–TX23A0.9821.1021.060
      B0.5100.6500.610
      TX24–TX210A1.0651.1101.105
      B0.6240.6600.648
      下载: 导出CSV
    • [1]

      Cai Meifeng. 1993. Review of principles and methods for rock stress measurement[J]. Chinese Journal of Rock Mechanics and Engineering, 12(3): 275−283 (in Chinese with English abstract).

      [2]

      Chen Nian, Wang Chenghu, Gao Guiyun, Wang Pu. 2021. Characteristics of in–situ stress field in the powerhouse area on the right bank of Baihetan based on stress polygon and borehole breakout method[J]. Rock and Soil Mechanics, 42(12): 3376−3384 (in Chinese with English abstract).

      [3]

      Cheng Yuanfang, Shen Haichao, Zhao Yizhong. 2008. A simplified differential strain analysis stress method for in–situ stress measurement[J]. Oil Drilling and Production Technology, 30(2): 61−64 (in Chinese with English abstract).

      [4]

      Dai Jinxing, Ni Yunyan, Wu Xiaoqi. 2012. Tight gas in China and its significance in exploration and exploitation[J]. Petroleum Exploration and Development, 39(3): 257−264 (in Chinese with English abstract).

      [5]

      Fu Xiaomin, Wang Xudong. 2007. The research on data processing about in–Situ stress measurement with AE[J]. Research and Exploration in Laboratory, (11): 296−299 (in Chinese with English abstract).

      [6]

      He Jianhua, Cao Feng, Deng Hucheng, Wang Yuanyuan, Li Yong, Xu Qinglong. 2023. Evaluation of in–situ stress in dense sandstone reservoirs in the second member of Xujiahe formation of the HC area of the Sichuan Basin and its application to dense sandstone gas development[J]. Geology in China: 50(4): 1107–1121 (in Chinese with English abstract).

      [7]

      Huang Rongzun, Zhuang Jinjiang. 1986. A new method for predicting formation fracture pressure[J]. Oil Drilling and Production Technology, (3): 1−14 (in Chinese).

      [8]

      Jia Chengzao, Zheng Min, Zhang Yongfeng. 2012. Unconventional hydrocarbon resources in China and the prospect of exploration and development[J]. Petroleum Exploration and Development, 39(2): 129−136 (in Chinese with English abstract).

      [9]

      Jiang Yongdong, Xian Xuefu, Xu Jiang. 2005. Research on application of Kaiser effect of acoustic emission to measuring initial stress in rock mass[J]. Rock and Soil Mechanics, 26(6): 946−950 (in Chinese with English abstract).

      [10]

      Li Guoxin, Lei Zhengdong, Dong Weihong, Wang Hongyan, Zheng Xingfan, Tan Jian. 2022. Progress, challenges and prospects of unconventional oil and gas development of CNPC[J]. China Petroleum Exploration, 27(1): 1−11 (in Chinese with English abstract).

      [11]

      Liu Yaqun, Li Haibo, Jing Feng, Luo Chaowen, Chen Bingrui, Li Junru, Zhou Qingchun. 2007. Determination of in–situ stress by hydraulic fracturing tests on preexisting fractures considering stress gradient and its engineering application[J]. Chinese Journal of Rock Mechanics and Engineering, 26(6): 1145−1149 (in Chinese with English abstract).

      [12]

      Mao H Y, Luo T T, Lai F Q, Zhang G T, Zhong L L. 2019. Experimental analysis and logging evaluation of in–situ stress of mud shale reservoir: Taking the deep shale gas reservoir of Longmaxi Formation in western Chongqing as an example[J]. IOP Conference Series:Earth and Environmental Science, 384(1): 012129.

      [13]

      Qiu Zhen, Zou Caineng, Li Jianzhong, Guo Qiulin, Wu Xiaozhi, Hou Lianhua. 2013. Unconventional petroleum resources assessment: Progress and future prospects[J]. Natural Gas Geoscience, 24(2): 238−246 (in Chinese with English abstract).

      [14]

      Schmitt D R, Currie C A, Zhang L. 2012. Crustal stress determination from boreholes and rock cores: Fundamental principles[J]. Tectonophysics, 580: 1−26. doi: 10.1016/j.tecto.2012.08.029

      [15]

      Shen Haichao, Cheng Yuanyuan, Wang Jingyin, Zhao Yizhong, Zhang Jianguo. 2008. Principal direction differential strain method for in–situ stress measurement and application[J]. Xinjiang Petroleum Geology, (2): 250−252 (in Chinese with English abstract).

      [16]

      Shi Lin, Zhang Xudong, Jin Yan, Chen Mian. 2004. New method for measurement of in–situ stresses at great depth[J]. Chinese Journal of Rock Mechanics and Engineering, (14): 2355−2358 (in Chinese with English abstract).

      [17]

      Wang Chenghu, Gao Guiyun, Wang Hong, Wang Pu. 2020. Integrated determination of principal stress and tensile strength of rock based on the laboratory and field hydraulic fracturing tests[J]. Journal of Geomechanics, 26(2): 167−174 (in Chinese with English abstract).

      [18]

      Wang Pu, Wang Chenghu, Yang Ruhua, Hou Zhengyang, Wang Hong. 2019. Preliminary investigation on the deep rock stresses prediction method based on stress polygon and focal mechanism solution[J]. Rock and Soil Mechanics, 40(11): 4486−4496 (in Chinese with English abstract).

      [19]

      Wang Xiaoqiong, Ge Hongkui, Song Lili, Hetai Ming, Xinwei. 2011. Experimental study of two types of rock sample acoustic emission events and Kaiser effect point recognition approach[J]. Chinese Journal of Rock Mechanics and Engineering, 30(3): 580−588 (in Chinese with English abstract).

      [20]

      Yin Shuai, Ding Wenlong, Lin Lifei, Liu Hanlin, Li Airong. 2023. Characteristics and the controlling effect on hydrocarbon accumulation of fractures in Yanchang Formation in Zhidan–Wuqi area, western Ordos Basin[J]. Earth Science, 48(7): 2614−2629 (in Chinese with English abstract).

      [21]

      Yin Shuai, Sun Xiaoguang, Wu Zhonghu, Wang Yingbin, Zhao Jingzhou, Sun Weihao, Yan Haolin. 2022. Coupling control of tectonic evolution and fractures on the Upper Paleozoic gas reservoirs in the northeastern margin of the Ordos Basin[J]. Journal of Central South University (Science and Technology), 53(9): 3724−3737 (in Chinese with English abstract).

      [22]

      Yin Xingyao, Ma Ni, Ma Zhengqian, Zong Zhaoyun. 2018. Review of in–situ stress prediction technology[J]. Geophysical Prospecting for Petroleum, 57(4): 488−504 (in Chinese with English abstract).

      [23]

      Zhang Chongyuan, Wu Manlu, Chen Qunze, Liao Chunting, Feng Chengjun. 2012. Reviews of in–situ stress measurement methods[J]. Journal of Henan Polytechnic University (Natural Science), 31(3): 305−310 (in Chinese with English abstract).

      [24]

      Zhang Jincai, Fan Xin, Huang Zhiwen, Liu Zhongqun, Qi Yuanchang. 2021. Assessment of anisotropic in–situ stressses in the upper Triassic Xujiahe Formation reservoirs in Western Sichuan Depression of the Sichuan Basin[J]. Oil and Gas Geology, 42(4): 963−972 (in Chinese with English abstract).

      [25]

      Zhang Kang, Zhang Liqin, Liu Dongmei. 2022. Situation of China’s oil and gas exploration and development in recent years and relevant suggestions[J]. Acta Petrolei Sinica, 43(1): 15−28,111 (in Chinese with English abstract). doi: 10.1038/s41401-021-00637-0

      [26]

      Zhao Gang, Dong Shier. 2009. The theory of the measurement of ground stress by hydraulic fracturing method and its application[J]. Shanxi Architecture, 35(36): 77−78 (in Chinese with English abstract).

      [27]

      Zheng Herong, Liu Zhongqun, Xu Shilin, Liu Zhenfeng, Liu Junlong, Huang Zhiwen, Huang Yanqing, Shi Zhiliang, Wu Qingzhao, Fan Lingxiao, Gao Jinhui. 2021. Progress and key research directions of tight gas exploration and development in Xujiahe Formation, Sinopec exploration areas, Sichuan Basin[J]. Oil and Gas Geology, 42(4): 765−783 (in Chinese with English abstract).

      [28]

      Zhou Wen. 2006. The Characteristics of In–situ Earth Stress and its Application Research in Engineering Geology of Petroleum on Compact Reservoir in Western Sichuan Depression[D]. Chengdu: Chengdu University of Technology, 1–142 (in Chinese with English abstract).

      [29]

      Zhu Haiyan, Song Yu jia, Xu Yun, Li Kuidong, Tang Xuanhe. 2021. Four–dimensional in–situ stress evolution of shale gas reservoirs and its impact on infill well complex fractures propagation[J]. Acta Petrolei Sinica, 42(9): 1224−1236 (in Chinese with English abstract).

      [30]

      Ziegler M, Valley B. 2021. Evaluation of the diametrical core deformation and discing analyses for in–situ stress estimation and application to the 4.9 km deep rock core from the Basel geothermal borehole, Switzerland[J]. Rock Mechanics and Rock Engineering, 54(12): 6511−6532. doi: 10.1007/s00603-021-02631-8

      [31]

      Zoback M D, Barton C A, Brudy M, Castillo D A, Finkbeiner T, Grollimund B R, Moos D B, Peska P, Ward C D, Wiprut D J. 2003. Determination of stress orientation and magnitude in deep wells[J]. International Journal of Rock Mechanics and Mining Sciences, 40(7): 1049−1076.

      [32]

      Zou Caineng, Tao Shizhen, Yang Zhi, Yuan Xuanjun, Zhu Rukai, Hou Lianhua, Jia Jinhua, Wang Lan, Wu Songtao, Bai Bin, Gao Xiaohui, Yang Chun. 2012. New advance in unconventional petroleum exploration and research in China[J]. Bulletin of Mineralogy, Petrology and Geochemistry, 31(4): 312–322 (in Chinese with English abstract).

      [33]

      Zou Caineng, Yang Zhi, He Dongbo, Wei Yunsheng, Li Jian, Jia Ailin, Chen Jianjun, Zhao Qun, Li Yilong, Li Jun, Yang Shen. 2018. Theory, technology and prospects of conventional and unconventional natural gas[J]. Petroleum Exploration and Development, 45(4): 575−587 (in Chinese with English abstract).

      [34] 蔡美峰. 1993. 地应力测量原理和方法的评述[J]. 岩石力学与工程学报, 12(3): 275−283.
      [35] 陈念, 王成虎, 高桂云, 王璞. 2021. 基于应力多边形与钻孔崩落的白鹤滩右岸厂房区地应力场特征研究[J]. 岩土力学, 42(12): 3376−3384.
      [36] 程远方, 沈海超, 赵益忠. 2008. 一种简化的差应变地应力测量技术[J]. 石油钻采工艺, 30(2): 61−64.
      [37] 戴金星, 倪云燕, 吴小奇. 2012. 中国致密砂岩气及在勘探开发上的重要意义[J]. 石油勘探与开发, 39(3): 257−264.
      [38] 付小敏, 王旭东. 2007. 利用岩石声发射测试地应力数据处理方法的研究[J]. 实验室研究与探索, (11): 296−299.
      [39] 何建华, 曹峰, 邓虎成, 王园园, 李勇, 徐庆龙. 2023. 四川盆地HC地区须二段致密砂岩储层地应力评价及其在致密气开发中的应用[J]. 中国地质, 50(4): 1107−1121.
      [40] 黄荣樽, 庄锦江. 1986. 一种新的地层破裂压力预测方法[J]. 石油钻采工艺, (3): 1−14.
      [41] 贾承造, 郑民, 张永峰. 2012. 中国非常规油气资源与勘探开发前景[J]. 石油勘探与开发, 39(2): 129−136.
      [42] 姜永东, 鲜学福, 许江. 2005. 岩石声发射Kaiser效应应用于地应力测试的研究[J]. 岩土力学, 26(6): 946−950.
      [43] 李国欣, 雷征东, 董伟宏, 王红岩, 郑兴范, 谭健. 2022. 中国石油非常规油气开发进展、挑战与展望[J]. 中国石油勘探, 27(1): 1−11.
      [44] 刘亚群, 李海波, 景锋, 罗超文, 陈炳瑞, 李俊如, 周青春. 2007. 考虑应力梯度的原生裂隙水压致裂法地应力测量的原理及工程应用[J]. 岩石力学与工程学报, 26(6): 1145−1149.
      [45] 邱振, 邹才能, 李建忠, 郭秋麟, 吴晓智, 侯连华. 2013. 非常规油气资源评价进展与未来展望[J]. 天然气地球科学, 24(2): 238−246.
      [46] 沈海超, 程远方, 王京印, 赵益忠, 张建国. 2008. 主方向差应变地应力测量方法[J]. 新疆石油地质, (2): 250−252.
      [47] 石林, 张旭东, 金衍, 陈勉. 2004. 深层地应力测量新方法[J]. 岩石力学与工程学报, (14): 2355−2358.
      [48] 王成虎, 高桂云, 王洪, 王璞. 2020. 利用室内和现场水压致裂试验联合确定地应力与岩石抗拉强度[J]. 地质力学学报, 26(2): 167−174.
      [49] 王璞, 王成虎, 杨汝华, 侯正阳, 王洪. 2019. 基于应力多边形与震源机制解的深部岩体应力状态预测方法初探[J]. 岩土力学, 40(11): 4486−4496.
      [50] 王小琼, 葛洪魁, 宋丽莉, 和泰名, 辛维. 2011. 两类岩石声发射事件与Kaiser效应点识别方法的试验研究[J]. 岩石力学与工程学报, 30(3): 580−588.
      [51] 尹帅, 丁文龙, 林利飞, 刘翰林, 李爱荣. 2023. 鄂尔多斯盆地西部志丹–吴起地区延长组裂缝特征及其控藏作用[J]. 地球科学, 48(7): 2614−2629.
      [52] 尹帅, 孙晓光, 邬忠虎, 王应斌, 赵靖舟, 孙伟豪, 闫浩霖. 2022. 鄂尔多斯盆地东北缘上古生界构造演化及裂缝耦合控气作用[J]. 中南大学学报(自然科学版), 53(9): 3724−3737.
      [53] 印兴耀, 马妮, 马正乾, 宗兆云. 2018. 地应力预测技术的研究现状与进展[J]. 石油物探, 57(4): 488−504.
      [54] 张重远, 吴满路, 陈群策, 廖椿庭, 丰成君. 2012. 地应力测量方法综述[J]. 河南理工大学学报(自然科学版), 31(3): 305−310.
      [55] 张金才, 范鑫, 黄志文, 刘忠群, 亓原昌. 2021. 四川盆地川西坳陷上三叠统须家河组储层各向异性地应力评价方法[J]. 石油与天然气地质, 42(4): 963−972.
      [56] 张抗, 张立勤, 刘冬梅. 2022. 近年中国油气勘探开发形势及发展建议[J]. 石油学报, 43(1): 15−28,111.
      [57] 赵刚, 董事尔. 2009. 水压致裂法测量地应力理论与应用[J]. 山西建筑, 35(36): 77−78. doi: 10.3969/j.issn.1009-6825.2009.36.048
      [58] 郑和荣, 刘忠群, 徐士林, 刘振峰, 刘君龙, 黄志文, 黄彦庆, 石志良, 武清钊, 范凌霄, 高金慧. 2021. 四川盆地中国石化探区须家河组致密砂岩气勘探开发进展与攻关方向[J]. 石油与天然气地质, 42(4): 765−783.
      [59] 周文. 2006. 川西致密储层现今地应力场特征及石油工程地质应用研究[D]: 成都: 成都理工大学, 1–142.
      [60] 朱海燕, 宋宇家, 胥云, 李奎东, 唐煊赫. 2021. 页岩气储层四维地应力演化及加密井复杂裂缝扩展规律[J]. 石油学报, 42(9): 1224−1236.
      [61] 邹才能, 陶士振, 杨智, 袁选俊, 朱如凯, 侯连华, 贾进华, 王岚, 吴松涛, 白斌, 高晓辉, 杨春. 2012. 中国非常规油气勘探与研究新进展[J]. 矿物岩石地球化学通报, 31(4): 312−322.
      [62] 邹才能, 杨智, 何东博, 位云生, 李剑, 贾爱林, 陈建军, 赵群, 李易隆, 李君, 杨慎. 2018. 常规–非常规天然气理论、技术及前景[J]. 石油勘探与开发, 45(04): 575−587.
    • 期刊类型引用(15)

      1. 刘礼广,熊志强,李和禄,刘福海,张慧龙,吴大天. 辽宁省瓦房店金刚石原生矿找矿方法组合模型的建立. 中国地质调查. 2024(01): 16-26 . 百度学术
      2. 贾志业,安战锋,赵廷严. 辽南瓦房店地区金伯利岩航磁垂向梯度异常特征及找矿意义. 地质与勘探. 2023(03): 617-626 . 百度学术
      3. 李兆丽,刘飞,郑永春,余晓艳,肖丙建,杨经绥. 金刚石的撞击成因与找矿新思考. 地质学报. 2023(10): 3493-3496 . 百度学术
      4. 吕青,刘飞,褚志远,葛跃进,刘潇,焦永鑫. 山东蒙阴三个金伯利岩带金刚石的矿物学特征与对比. 地质学报. 2022(04): 1225-1238 . 百度学术
      5. 刘礼广,朱婉婷,王涛,王旭. 辽南42号金伯利岩管地质特征及找矿预测. 地质找矿论丛. 2022(02): 174-190 . 百度学术
      6. 刘飞,杨经绥,冯光英,牛晓露,李观龙,张超凡. 古太平洋板片晚二叠-早三叠世俯冲后撤:来自海南岛弧岩浆作用的制约. 岩石学报. 2022(11): 3455-3483 . 百度学术
      7. 刘礼广,吴大天,韩双,孙海涛,李和禄,熊志强. 辽南瓦房店金刚石矿田110号金伯利岩管的地质特征及成因探讨. 矿产勘查. 2021(06): 1339-1354 . 百度学术
      8. 刘礼广,吴大天,熊志强. 辽宁省金刚石原生矿成矿规律与成矿模式研究. 山东国土资源. 2021(10): 29-36 . 百度学术
      9. 张鑫刚,周佐民,袁启玉. 博茨瓦纳重要矿产资源潜力与投资前景分析. 地质与勘探. 2021(06): 1297-1303 . 百度学术
      10. 张国仁,仲米山,潘玉啟,付海涛,李显东,高永钊,高福亮,吴子杰,王海鹏. 辽宁瓦房店金刚石成矿带发现逆冲推覆构造. 地质与资源. 2020(03): 294-298 . 百度学术
      11. 付海涛. 辽宁瓦房店金刚石矿田金伯利岩侵位机制分析. 地质学报. 2020(09): 2640-2649 . 百度学术
      12. 程光锁,刘卫东,褚志远,周军,赵德庆,胡戈,梁太涛. 山东蒙阴县西峪金伯利岩带深部综合地球物理特征及其意义. 地质学报. 2020(09): 2772-2782 . 百度学术
      13. 刘飞,杨经绥,连东洋,熊发挥,吴魏伟. 蛇绿岩中金刚石和异常矿物的全球记录. 地质学报. 2020(09): 2588-2605 . 百度学术
      14. 刘礼广,吴大天,韩双,万方来,林少全,杨献忠. 金伯利岩型金刚石矿床研究及其成矿模式探讨——以辽宁瓦房店地区金刚石原生矿矿床为例. 地质学报. 2020(09): 2650-2665 . 百度学术
      15. 向璐,郑建平. 中国原生金刚石形成条件及找矿思考. 地质学报. 2020(09): 2574-2587 . 百度学术

      其他类型引用(2)

    图(15)  /  表(5)
    计量
    • 文章访问数:  1945
    • HTML全文浏览量:  1035
    • PDF下载量:  396
    • 被引次数: 17
    出版历程
    • 收稿日期:  2023-03-05
    • 修回日期:  2023-05-22
    • 网络出版日期:  2024-02-03
    • 刊出日期:  2024-01-24

    目录

    /

    返回文章
    返回
    x 关闭 永久关闭