青藏高原中南部地壳应力与变形特征：来自电性各向异性研究的启示

一、 研究概况

本研究的主要作者为Zhehan Liu（第一作者，现工作于义乌工商职业技术学院）、Hao Dong（通讯作者）、Sheng Jin（通讯作者）、Wenbo Wei、Gaofeng Ye和Letian Zhang。所有作者均来自中国地质大学（北京）地球物理与信息技术学院，部分作者同时隶属于地质过程与矿产资源国家重点实验室，Sheng Jin也同时任职于河北地质大学。该研究以题为“Stress and deformation characteristics of the central‐southern Tibetan crust: insights from electrical anisotropy studies”发表于*Journal of Geophysical Research: Solid Earth*期刊，于2026年发表（接收于2026年1月31日，在线发表于2024年10月21日）。该研究得到了国家自然科学基金等多个项目的资助。

二、 学术背景

本研究属于固体地球物理学领域，聚焦于利用地球物理方法探测大陆碰撞带深部结构与动力学过程。青藏高原作为印度板块与欧亚板块碰撞的产物，其隆升机制和内部变形过程是地球科学界长期争论的焦点。地壳流变学性质（即材料在应力下的变形行为）是理解高原变形模式的关键。以往的地震学研究揭示了高原广泛存在的地壳低速带和复杂的各向异性特征，而大地电磁法（Magnetotelluric, MT）研究则普遍发现了中下地壳存在大规模高导层，通常被解释为部分熔融体和/或流体的存在。然而，不同地球物理方法（如MT与地震学）估算的熔体分数存在显著差异（MT估算通常>10%，而地震学估算为2%-7%），且对高导层的成因及其与变形机制的联系缺乏系统性的电性各向异性约束。电性各向异性，即岩石电阻率随测量方向变化的特性，如同地震各向异性一样，能够反映深部结构的定向排列特征（如裂隙、矿物晶格或熔体囊的定向），从而为揭示地壳应力场和变形机制提供独特视角。因此，本研究旨在通过分析穿越青藏高原中南部（拉萨地块和喜马拉雅地块）的一条MT剖面数据，构建二维电性各向异性模型，以揭示该区域地壳的电性结构，特别是电性各向异性的分布特征，进而探讨其与地壳流变性质、应力状态和变形机制的关联，最终为理解高原隆升和地壳变形过程提供新的地球物理约束。

三、 详细研究流程

本研究流程严谨，主要包括数据采集与处理、维度分析、各向同性与各向异性反演建模、以及基于模型的熔体分数估算和流变学分析。

1. 数据采集与处理： 研究使用了来自“深部探测技术与实验研究专项”（Sinoprobe）和“中国地质大学（北京）3000线”（CUGb line-3000）项目的MT数据。剖面共包含22个宽频MT站点和4个长周期MT站点。数据采集使用Phoenix MTU-5A和LEMI-417仪器，分别记录了约24小时和7天的电磁时间序列。数据处理采用了稳健的统计算法和远程参考方法，估算了周期范围为0.01至10,000秒的MT和地磁传递函数。由于高原地区人文噪声低，数据质量普遍较高。

2. 维度与电性走向分析： 在反演前，研究团队对MT数据进行了系统的维度分析，以确定地下结构的复杂程度和最佳电性走向。使用了三种主要方法： * 相位张量（Phase Tensor）分析： 通过绘制相位张量椭圆及其偏斜角（|β|），评估结构的维度。结果显示，大部分站点在多数周期表现为一维或二维结构，但在缝合带和裂谷附近的长周期数据中观察到较大的偏斜角，暗示可能存在三维效应或各向异性。 * 感应矢量（Induction Vector）分析： 绘制了帕金森约定的实感应矢量图。短周期（0.1-10秒）矢量方向多变，反映局部异常；而长周期（100-1000秒）矢量变得短小且整体呈南北向排列，这可能指示存在一个电性各向异性层。 * Groom-Bailey张量分解： 用于估算区域电性走向，并分离 galvanic distortion（ galvanic distortion）。分析发现，随着周期增大，剖面可以一条虚线为界分为两部分：北侧站点优选走向约为N115°E，南侧站点优选走向近东西向。 综合这些分析，研究人员认为尽管存在局部三维效应，但整体可用二维框架进行近似。为了进行后续二维反演，将所有MT数据旋转到N90°E坐标系（即东西向），该方向大致平行或垂直于主要地表构造。

3. 建模与反演： 研究使用了MARE2DEM软件进行二维各向同性和各向异性反演。该软件采用非结构化三角形网格和有限元算法进行正演模拟，并利用Occam反演方法迭代求解最光滑的电阻率模型。 * 各向同性反演： 作为基线模型。反演联合使用了TE模式、TM模式和Tipper数据。最终模型显示，在雅鲁藏布江缝合带（YZS）以北的上地壳出现了一系列电阻率和导电性异常交替出现的现象，中下地壳也显示出类似的交替序列。作者指出，这种现象可能与各向异性数据被不充分地用各向同性模型拟合所产生的假象有关。 * 各向异性反演： 鉴于维度分析显示各向异性迹象以及各向同性模型出现的疑似假象，研究进行了二维各向异性反演。在MARE2DEM中，电导率用轴向各向异性模型表示，允许电阻率在三个主轴方向（平行于剖面走向的ρ_xx、垂直于走向的ρ_yy和垂直方向的ρ_zz）变化。本研究主要关注方位角各向异性（ρ_yy ≠ ρ_xx）。通过合成测试，选择了合适的各向异性惩罚权重因子（α=0.3）。各向异性反演得到的归一化均方根误差（RMS misfit）从8.67降至1.46，优于各向同性模型（1.74）。各向异性模型清晰地揭示了一个关键特征：在中下地壳存在一个双层电性各向异性结构。上层的各向异性轴（即电流更易流动的方向）近南北向（ρ_yy < ρ_xx），深度约20-35公里（中地壳）；下层的各向异性轴近东西向（ρ_xx < ρ_yy），深度从约35公里延伸至莫霍面（下地壳）。 * 引入撕裂不连续面的反演： 为了解决MT反演中高导体下边界可能向下“涂抹”（smearing）到地幔的问题，并更清晰地约束各向异性层的边界，研究在模型中引入了两个“撕裂”不连续面（tear discontinuities）。上不连续面（D1）设置在约35公里深度，对应于中下地壳电性各向异性转变的边界；下不连续面（D2）设置在根据接收函数研究得到的莫霍面深度（约55-60公里）。引入不连续面后，反演结果的基本特征保持不变，但高导体的下边界更加清晰，模型更为合理。

4. 基于模型的分析： * 有效电导率与熔体分数估算： 研究计算了中地壳（20-30公里）和下地壳（35-55公里）导电层的有效体积电导率（σ_eff）。结果显示，中地壳表现为南北向各向异性（σ_yy > σ_xx），而下地壳表现为东西向各向异性（σ_xx > σ_yy），各向异性强度可达3-5倍。为了解释高导层的成因并估算熔体分数，研究首先采用了修正的阿尔奇定律（Modified Archie‘s Law）进行各向同性估算。然而，该模型要求过高的熔体分数（例如，电导率0.2 S/m对应约10.9%的熔体），与地震学估算存在矛盾，且无法解释观测到的电性各向异性。 * 各向异性导电模型： 为了解决上述矛盾，研究引入了一个各向异性导电模型（基于Berryman and Hoversten， 2013的自洽经验公式）。该模型假设熔体/流体以椭球状包裹体形式存在，并在应力作用下定向排列（即具有特定的纵横比α）。计算表明，当椭球体被拉长（α变小，如α=0.1，即长短轴比10：1）时，更低的熔体分数（如3.4%-7.3%）就能产生与观测相符的有效电导率（如0.2 S/m）和强烈的各向异性。这一模型为调和MT与地震学估算的熔体分数差异提供了新思路，并强有力地支持了“应力导致熔体囊定向排列”的假设。 * 流变学分析： 利用经验公式将电导率转换为有效粘度。计算表明，在下地壳的软弱区，粘度在10^18–10^20 Pa·s之间，这与数值地球动力学研究的估计一致。粘度低于约10^19 Pa·s时，地壳具备发生塑性流动（地壳流）的能力。这种低粘度环境使得部分熔融物质在应力作用下易于发生定向变形。

四、 主要研究结果

1. 揭示双层电性各向异性结构： 通过二维各向异性反演，本研究首次在青藏高原中南部地壳清晰地揭示了一个双层电性各向异性结构。中地壳（~20-35公里）表现为近南北向的电性各向异性（电流更易南北向流动），而下地壳（~35公里至莫霍面）则表现为近东西向的电性各向异性。这一特征是本研究最核心的发现。
2. 各向异性模型优于各向同性模型： 各向异性反演不仅显著降低了数据拟合误差，而且消除了各向同性模型中出现的、难以用地质构造解释的电阻率与导电体交替出现的假象。模型还更好地刻画了可能代表俯冲印度板块的高阻体（R1）。
3. 各向异性起源与变形机制关联： 研究结合岩石物理实验和地质背景，对双层各向异性的成因提出了解释：
   • 下地壳（东西向各向异性）： 对应于韧性变形域。该深度温度压力条件（~800-845°C， ~1.2 GPa）可能引发角闪岩的脱水熔融，产生部分熔体。在区域南北向挤压应力背景下，熔体囊被拉长、定向排列，其长轴趋向于东西向（垂直于主压应力方向或平行于伸展方向）。电流更易沿熔体囊的长轴方向（东西向）流动，从而形成观测到的东西向电性各向异性。各向异性导电模型表明，仅需较低（~7%）且更合理的熔体分数即可解释观测到的高导性和强各向异性。
   • 中地壳（南北向各向异性）： 对应于脆-韧性过渡域（半脆性变形）。该深度温度可能不足以引发大规模熔融，但富角闪石的中地壳在东西向伸展应力作用下，可能通过半脆性蠕变形成南北向排列的微裂隙网络。富含离子的流体充填于这些定向微裂隙中，导致电流优先沿南北向（平行于微裂隙走向）传导，形成南北向电性各向异性。
   • 上地壳： 主要表现为脆性变形，其各向异性可能与宏观断裂、岩性不均一性或地堑的倾角变化有关。
4. 提出“双层电性各向异性模型”： 基于上述发现和解释，研究提出了一个概念模型（图9，图10），将观测到的电性各向异性模式与地壳不同层次的变形机制和应力状态直接联系起来。中上地壳的近南北向各向异性反映了近脆性（包括脆性和半脆性）变形，与南北向挤压应力下形成的定向裂隙/微裂隙相关；而下地壳的近东西向各向异性则反映了韧性变形，与应力作用下定向排列的熔体囊相关。两者之间的边界（D1，~35公里）可能代表了从富水流体主导到部分熔融主导的流变学转换界面。

五、 研究结论与意义

本研究通过对青藏高原中南部MT数据的各向异性反演和分析，得出了以下核心结论： 1. 青藏高原中南部地壳存在一个清晰的双层电性各向异性结构，其中地壳和下地壳的各向异性优势方向正交。 2. 这种各向异性结构很可能是不同深度地壳物质在不同变形机制（上中地壳的近脆性变形与下地壳的韧性变形）和区域应力场（南北向挤压与东西向伸展）共同作用下的综合体现。 3. 引入各向异性导电模型后，可以用更低的、与地震学估算更为一致的熔体分数来解释观测到的高导异常，为调和不同地球物理方法在估算熔体含量上的分歧提供了新途径。 4. 电性各向异性可以作为探测地壳应力、变形和流变学分层的有效诊断指标，为理解大陆碰撞带的深部动力学过程提供了新的视角和方法。

本研究的科学价值在于： * 提供了新的深部探测约束： 首次系统揭示了高原中南部地壳的电性各向异性精细结构，为高原隆升和地壳变形模型提供了关键的地球物理证据。 * 建立了电性各向异性与变形机制的桥梁： 将观测到的电性各向异性模式与具体的岩石变形机制（脆性破裂、半脆性蠕变、韧性流动）和熔体/流体分布联系起来，深化了对电性异常成因的理解。 * 提出了创新的解释框架： “双层电性各向异性模型”及其与应力-变形分层的关联，为解释其他构造活动区的类似观测提供了一个可检验的模板。 * 推动了方法学发展： 展示了各向异性反演在解决复杂构造区电性结构中的重要性，以及结合各向异性导电模型来更合理估算熔体分数的潜力。

六、 研究亮点

1. 重要发现： 明确揭示了青藏高原中南部地壳存在“中地壳南北向、下地壳东西向”的双层电性各向异性结构，这是该研究最突出的原创性发现。
2. 方法新颖： 系统应用并展示了二维电性各向异性反演在高原复杂构造区解析深部结构的优势；创新性地引入“各向异性导电模型”来解释高导层的成因和估算熔体分数，有效调和了地球物理观测之间的矛盾。
3. 多学科交叉融合： 研究深度整合了大地电磁学、岩石物理学、流变学、构造地质学等多学科知识，不仅给出了地球物理图像，还深入探讨了其背后的物理机制和地质含义，形成了从观测到解释的完整逻辑链条。
4. 明确的构造指示意义： 将电性各向异性特征直接与地壳的应力状态和变形类型（脆性/半脆性/韧性）相关联，使地球物理观测成为解读构造动力学过程的直接工具。

七、 其他有价值内容

研究也讨论了模型的局限性，例如二维反演对可能存在的三维结构分辨率有限，以及数据周期范围对约束地幔各向异性的不足。作者建议未来应利用区域阵列MT数据进行三维各向异性反演，并布设更多长周期站点以探测岩石圈地幔的应力与变形机制。同时，指出目前高温高压实验对下地壳导电性的研究多集中于各向同性样品，未来需要开展更多针对各向异性样品的实验，以定量化应力与电性各向异性之间的关系。这些都为后续研究指明了方向。