通过部分暴露的I型-III型混合模式岩脉来估算岩浆的过压

【字体: 时间:2026年03月13日 来源:Journal of Structural Geology 2.9

编辑推荐:

  本研究扩展单侧椭圆法至混合模式I-III岩脉,通过自然与模拟数据验证,并应用于Dharwar克拉通的镁质岩脉,可靠估算岩浆过压(11-23 MPa)和源岩深度(3.8-7.8 km)。

  
阿努巴布·德伊(Anubhab Dey)| 西尔申杜·库马尔·比斯瓦斯(Sirshendu Kumar Biswas)| 苏拉夫·蒙达尔(Sourav Mondal)| 特里迪布·库马尔·蒙达尔(Tridib Kumar Mondal)
印度统计研究所地质研究部,203 B.T.路,加尔各答-700108,西孟加拉邦,印度

摘要

我们探讨了是否可以利用部分暴露的混合模式I-III岩脉来估算岩浆的过压。本研究扩展了之前用于从单端暴露岩脉中估算I型岩脉全尺寸的方法。这使得可以对混合模式I-III岩脉进行类似的重建,从而在遇到单端暴露的混合模式I-III岩脉时,能够一致地估算岩脉的尺寸、岩浆过压和起源深度。基于野外观察,我们认为岩脉最初的I型开口的残余椭圆性在随后可观察到的混合模式I-III几何结构中得以保留。测试数据集包括新墨西哥州Ship Rock的完全暴露的岩脉、印度东达尔瓦尔克拉通(Dharwar Craton)Closepet花岗岩中的长英质脉状体,以及使用COMSOL Multiphysics 6.2生成的两个数值模拟裂缝,以复制自然混合模式I-III特征的几何特性。结果表明,该方法能够为混合模式I-III岩脉的单端暴露提供稳健可靠的几何尺寸估算。

引言

岩浆从深处侵入地壳较浅层时所产生的压力可能会导致原本完整的围岩发生断裂(Anderson, 1939, 1951; Spence and Turcotte 1985; Lister and Kerr 1991)。这种由岩浆驱动的断裂形成了称为岩脉的平面岩浆体,它们是填充有岩浆的伸展裂缝。岩脉及其他由岩浆驱动形成的片状侵入体的几何形状对其形成机制提供了关键约束,同时记录了它们形成过程中的古应力状态以及局部和区域应力场的变化(Delaney et al., 1986; Glazner et al., 1999; Martinez-Poza et al., 2014; Stephens et al., 2018)。因此,它们的二维(2D)几何形状被用来理解远场应力,而它们的方向则通过三维(3D)莫尔圆(Mohr-circle)分析被广泛用于推断主应力的相对大小和方向(Jolly & Sanderson, 1997; Yamaji and Sato, 2011; Martinez-Poza & Druguet, 2016; Pascal, 2021)。因此,岩脉的几何形状在其基于实地的定量研究中起着重要作用。
当岩浆使原本完整的围岩发生断裂时,断裂形成的壁面会相互远离,沿着最小主应力轴移动。这导致岩脉以I型断裂的形式打开(Gudmundsson 2020; Pollard and Martel, 2020)。根据I型断裂形成的典型模型(Parker, 1981; Gudmundsson, 1983a, 1983b),以下假设成立:(i)围岩是均匀的、各向同性的且具有弹性;(ii)在断裂形成过程中最小主应力保持不变;(iii)最小主应力实际上是拉应力,并且垂直于断裂面。在岩浆驱动的断裂条件下,这些假设表明岩脉以椭圆形的I型断裂打开(见补充图表-1),其长宽比定义为沿走向方向的最大厚度与总长度的比值,该比值远小于1。这种椭圆形几何结构在大致垂直于岩脉体的剖面中观察最为明显。本研究中所有关于I型岩脉的几何描述和测量均基于此类剖面的观察结果。这种椭圆形几何结构的长宽比是进行关键定量分析的基础,包括估算岩浆过压和岩脉起源深度,前提是已知围岩的杨氏模量、泊松比和密度,以及侵入岩浆的密度(Gudmundsson, 1983a, 1983b; Kusumoto et al., 2013; Becerril et al., 2013; Biswas et al., 2023)。
岩脉研究中一个持续存在的实际挑战是部分野外暴露现象,即只有岩脉的一个尖端或根本没有尖端暴露出来。这限制了对其实际几何形状的准确重建。为了解决这一挑战,比斯瓦斯等人(Biswas et al. 2023)提出了一种几何方法,用于在仅暴露出一个尖端的情况下估算岩脉的完整尺寸。
该方法假设暴露的岩脉尖端与二维笛卡尔坐标系的原点重合,其中横坐标和纵坐标分别对应于椭圆岩脉露头的主轴和次轴(见补充图表-2)。然后根据沿岩脉暴露长度的系统厚度测量,在岩脉周边标出点。接着通过所有可能的点对构建椭圆。必须确保每个椭圆的主轴与岩脉的主轴平行,并且其中一个尖端与暴露的岩脉尖端重合。通过平均所有新生成的椭圆对应的参数,可以获得部分暴露岩脉的完整尺寸。这种方法此后被称为“单尖端椭圆法”,对于理想的椭圆形岩脉(I型)能够提供稳健的结果。然而,对于几何形状偏离这种理想形式的岩脉,其适用性仍需进一步研究。
野外研究表明,由于多种因素,岩脉的形状和几何结构往往与理论上的椭圆形有显著偏差,包括围岩的异质性(Gudmundsson, 2020; Pollard & Martel, 2020)、岩浆固化过程中岩脉尖端的楔形变形以及围岩的非弹性变形(Daniels et al., 2012)、主应力轴方向的变化(Pollard & Martel, 2020; Biswas & Mondal, 2024)、由于弹性自由表面(地球表面)的存在导致岩脉上端的前伸膨胀(Pollard et al., 1983; Rivalta & Dahm, 2006),以及存在层理或预先存在的断裂等结构各向异性(Jolly & Sanderson, 1995; Delaney et al., 1986)。Rivalta等人(2015)和Townsend等人(2017)的综述详细总结了在不同因素作用下岩脉传播的力学模型。虽然多种因素影响岩脉的最终几何形状,但本研究重点考察了由主应力轴重新定向形成的几何结构。
在岩脉侵入初始完整的围岩时,形成的断裂壁面会相互远离,沿着最小主应力轴移动。这导致岩脉以I型断裂的形式打开(Gudmundsson 2020; Pollard and Martel, 2020)。根据I型断裂形成的典型模型(Parker, 1981; Gudmundsson, 1983a, 1983b),以下假设成立:(i)围岩是均匀的、各向同性的且具有弹性;(ii)在断裂形成过程中最小主应力保持不变;(iii)最小主应力实际上是拉应力,并且垂直于断裂面。在岩浆驱动的断裂条件下,这些假设表明岩脉以椭圆形的I型断裂打开(见补充图表-1),其长宽比定义为沿走向方向的最大厚度与总长度的比值,该比值远小于1。这种椭圆形几何结构在大致垂直于岩脉体的剖面中观察最为明显。本研究中所有关于I型岩脉的几何描述和测量均基于此类剖面的观察结果。这种椭圆形几何结构的长宽比是进行关键定量分析的基础,包括估算岩浆过压和岩脉起源深度,前提是已知围岩的杨氏模量、泊松比和密度,以及侵入岩浆的密度(Gudmundsson, 1983a, 1983b; Kusumoto et al., 2013; Becerril et al., 2013; Biswas et al., 2023)。
岩脉研究中的一个持续实际挑战是部分野外暴露现象,即只有给定椭圆露头的一个尖端或没有尖端暴露出来。这限制了对实际几何形状的准确重建。为了解决这一挑战,比斯瓦斯等人(Biswas et al. 2023)提出了一种几何方法,用于在仅暴露出一个尖端的情况下估算岩脉的完整尺寸。
该方法假设暴露的岩脉尖端与二维笛卡尔坐标系的原点重合,其中横坐标和纵坐标分别对应于椭圆岩脉露头的主轴和次轴(见补充图表-2)。然后根据沿岩脉暴露长度的系统厚度测量,在岩脉周边标出点。接着通过所有可能的点对构建椭圆。必须确保每个椭圆的主轴与岩脉的主轴平行,并且其中一个尖端与暴露的岩脉尖端重合。通过平均所有新生成的椭圆对应的参数,可以获得部分暴露岩脉的完整尺寸。这种方法此后被称为“单尖端椭圆法”,对于理想的椭圆形岩脉(I型)能够提供稳健的结果。然而,对于几何形状偏离这种理想形式的岩脉,其适用性仍需进一步研究。
野外研究表明,由于多种因素,岩脉的形状和几何结构往往与理论上的椭圆形有显著偏差,包括围岩的异质性(Gudmundsson, 2020; Pollard & Martel, 2020)、岩浆固化过程中岩脉尖端的楔形变形以及围岩的非弹性变形(Daniels et al., 2012)、主应力轴方向的变化(Pollard & Martel, 2020; Biswas & Mondal, 2024)、由于存在弹性自由表面(地球表面)导致岩脉上端的前伸膨胀(Pollard et al., 1983; Rivalta & Dahm, 2006),以及存在层理或预先存在的断裂等结构各向异性(Jolly & Sanderson, 1995; Delaney et al., 1986)。Rivalta等人(2015)和Townsend等人(2017)的综述详细总结了在不同因素作用下岩脉传播的力学模型。虽然多种因素影响岩脉的最终几何形状,但本研究重点考察了由主应力轴重新定向形成的几何结构。
在岩脉传播过程中,局部应力场会发生变化,主应力轴的方向也会改变。这导致混合模式加载,包括初次打开(I型)和平面内剪切(II型)以及平面外剪切(III型)的组合(图1(a)(Delaney & Pollard, 1981; Nicholson & Pollard, 1985; Pollard & Martel, 2020)。本研究考察了由混合模式I-III加载形成的岩脉几何结构,这是由于III型加载叠加在最初以I型断裂打开的岩脉上所致。这些岩脉表现出分段的几何形状,各个段落以系统性的阶梯状排列,反映了剪切的方向。由于岩脉最初以I型断裂打开,因此与初次I型开口相关的残余椭圆性很可能在混合模式I-III岩脉中得到保留,正如本研究中从尖端到尖端的暴露测量结果所示。因此,我们首先使用单尖端椭圆法来测试从(i)新墨西哥州Ship Rock的东北岩脉(Delaney and Pollard 1981);(ii)印度卡纳塔克邦Pavagada的长英质脉状体(Biswas and Mondal 2024);(iii)使用基于有限元方法的软件生成的类似阶梯状裂缝的两种不同几何结构获得的混合模式I-III几何结构,以确定该方法对混合模式I-III岩脉单端暴露的适用性。在此过程中,我们开发了一种基于单端暴露的混合模式I-III岩脉几何尺寸量化的协议。最后,该方法应用于印度达尔瓦尔克拉通Closepet花岗岩中部分暴露的基性岩脉,以估算它们的岩浆过压和起源深度。在深入研究之前,有必要讨论混合模式I-III岩脉的渐进演化过程以及纯I型几何结构的残余椭圆性的保留情况,这将在以下部分详细阐述。

部分内容

混合模式I-III岩脉:形成机制和露头特征

岩脉最初的I型形成使其沿最大主压应力轴(σ1)传播,并沿最小主压应力轴(σ3)打开(Pollard et. al., 1982; Lister and Kerr 1991; Pollard and Martel, 2020, p. 380)。为了清晰起见,假设σ1轴是垂直的,而σ2轴和σ3轴位于水平面上(图1(b))。
当岩脉沿σ1轴传播时,周围应力场可能会发生变化,

测量概念

获取混合模式I-III岩脉的长度和厚度数据需要稍微修改比斯瓦斯等人(Biswas et al. 2023)提出的方法。测量首先确定岩脉露头的暴露尖端。一旦找到尖端,就确定一个适当的长度间隔,以便系统地记录其暴露轨迹上的岩脉厚度。间隔的选择至关重要:它既不能太大以至于只收集到少量数据点,也不能太小以至于

验证

单尖端椭圆法已针对具有混合模式I-III几何结构的自然和模拟数据集进行了测试。这些数据集包括:i) 新墨西哥州Ship Rock的记录详尽的东北岩脉(Delaney and Pollard, 1981);ii) 印度卡纳塔克邦Pavagada的长英质脉状体(Biswas and Mondal, 2024);以及iii) 使用COMSOL Multiphysics 6.2生成的两个数值模拟裂缝。

应用

单尖端椭圆法应用于在Closepet花岗岩转移带内花岗岩宿主岩中观察到的部分暴露的单尖端混合模式I-III基性岩脉露头,以估算它们的完整尺寸。研究区域在图6中用虚线框标出,位于Closepet花岗岩的转移带内。Moyen等人(2003)的磁化率(AMS)各向异性分析表明存在南北走向的磁层理

结果与讨论

单尖端椭圆法已在自然环境和数值模拟模型中观察到的混合模式I-III几何结构的数据集上进行了测试。该数据集涵盖了多种宿主岩性、力学性质和尺度,为验证该方法的适用性提供了全面的基础。
Ship Rock的东北岩脉的结果表明,单尖端椭圆法在预测次轴方面表现出出色的准确性和稳定性

结论

  • 该研究在验证了该方法对混合模式I-III岩脉的适用性后,将其扩展到了之前提出的单尖端椭圆法。
  • 通过对自然混合模式岩脉几何结构和数值模拟的混合模式几何结构进行测试,证实了即使在岩脉仅暴露了50%的情况下,该方法也能提供稳健可靠的岩脉尺寸估算。
  • 该方法应用于部分暴露的单尖端混合模式I-III基性岩脉

利益冲突声明

CRediT作者贡献声明

特里迪布·库马尔·蒙达尔(Tridib Kumar Mondal):撰写 – 审稿与编辑、可视化、监督、项目管理、资金获取、数据管理、概念化。苏拉夫·蒙达尔(Sourav Mondal):撰写 – 审稿与编辑、初稿撰写、软件开发、方法论、形式分析。西尔申杜·库马尔·比斯瓦斯(Sirshendu Kumar Biswas):撰写 – 审稿与编辑、初稿撰写、软件开发、方法论、数据管理、概念化。阿努巴布·德伊(Anubhab Dey):初稿撰写、软件开发、方法论、调查、形式分析,

未引用参考文献

Appledorn and Wright Jr, 1957; Das et al., 2025; Dey et al., 2015; Hakim and Karma, 2009; Martínez-Poza et al., 2014; Moyen et al., 2003; Ramakrishnan and Vaidyanadhan, 2008.

写作过程中生成式AI和AI辅助技术的声明

在准备本工作时,作者使用了ChatGPT-5来改进语言和可读性。使用该工具后,作者根据需要对内容进行了审查和编辑,并对出版物的内容负全责。

利益冲突声明

? 作者声明他们没有已知的可能会影响本文所述工作的财务利益或个人关系。

致谢

本研究是AD在印度统计研究所资助的博士研究的一部分。该研究得到了印度统计研究所的研究资助和MoES
相关新闻
生物通微信公众号
微信
新浪微博
  • 搜索
  • 国际
  • 国内
  • 人物
  • 产业
  • 热点
  • 科普

热点排行

    今日动态 | 人才市场 | 新技术专栏 | 中国科学人 | 云展台 | BioHot | 云讲堂直播 | 会展中心 | 特价专栏 | 技术快讯 | 免费试用

    版权所有 生物通

    Copyright© eBiotrade.com, All Rights Reserved

    联系信箱:

    粤ICP备09063491号