Deformation Due to a Pair of Infinite and a Pair of infinite-finite planar faults in Fractional Standard Linear Solid

技术摘要：分数阶标准线性固体模型中一对无限大与一对无限-有限平面断层的变形研究

问题陈述
地震活跃区域经常包含多个相互作用的断层，其中一个断层的滑动会显著影响相邻断层的应力积累与分布。虽然以往的研究利用整数阶导数在粘弹性介质（如 Maxwell 或标准线性固体模型）中对断层相互作用进行了建模，但这些传统的处理方法往往忽略了地质材料固有的“记忆”效应。本文通过模拟两种特定的相互作用蠕动走滑断层场景，研究了岩石圈-软流圈系统在无震期的变形情况：

1. 情形 I： 两条无限大埋藏断层之间的相互作用。
2. 情形 II： 一条无限大埋藏断层与一条有限长地表破裂断层之间的相互作用。

本研究假定地球介质表现为一种由分数阶标准线性固体 (FSLS) 模型表征的粘弹性半空间，与整数阶模型相比，该模型能更好地捕捉地壳岩石随时间变化的流变行为和记忆效应。

研究方法
作者开发了理论数学模型，以推导位移、应力和应变分量的解析解。研究方法包含以下关键步骤：

• 数学公式化： 使用包含 Caputo 分数阶导数算子 (${}_{0}D^{\alpha}_{t}$) 的分数阶标准线性固体本构方程来定义介质，其中阶数 $\alpha \in (0, 1]$。应力-应变关系结合了有效刚度 ($\mu$) 和有效粘度 ($\eta$)。
• 控制方程： 分析假设为准静态变形（忽略惯性力）。针对无限大和有限断层几何形状，建立了应力运动方程和边界条件。
• 求解技术： 作者采用 拉普拉斯变换 (Laplace transformation) 来处理随时间变化的分数阶微分方程。解的推导利用了：
  • 对应原理 (Correspondence Principle) (Rybicki, 1971)。
  • 修正格林函数法 (Modified Green's Function technique) (Maruyama, 1966)，以解释断层面两侧的位错条件。
  • 米洛塔函数 (Mittag-Leffler function) ($E_{\alpha}$)，该函数源于分数阶算子的拉普拉斯逆变换。
• 分析场景：
  • 推导了断层蠕动之前（构造应力积累）和断层蠕动之后（由于滑动导致的位错）的状态解。
  • 应用蠕动速度 ($V_1, V_2$) 和随深度变化的蠕动函数来模拟断层 $F_1$ 和 $F_2$ 的运动。
  • 使用基于地质数据（例如，刚度 $\mu = 3.5 \times 10^{10}$ N/m²，粘度 $\eta = 5 \times 10^{19}$ Pa·s）的参数进行数值计算，以生成针对不同断层倾角 ($\theta$)、分数阶 ($\alpha$) 和蠕动速度的图形结果。

主要贡献

1. 断层相互作用的分数阶建模： 本文通过引入分数阶微积分，扩展了现有的断层相互作用模型，从而能够更准确地表征岩石圈具有记忆依赖性的粘弹性行为。
2. 复杂几何结构的解析解： 研究为两种不同的配置——一对无限大埋藏断层以及无限大与有限长地表破裂断层的混合系统——提供了位移、应力和应变场的闭式解析表达式。
3. 量化记忆效应： 研究明确展示了分数阶 $\alpha$ 如何影响应力的积累与释放，将数学参数与地质介质的物理“记忆”联系起来。

结果与数值发现
支持数值分析的图形化表示得出以下具体观察结果：

• 位移量级： 断层附近的表面位移变化量级为 $10^{-3}$ m/年。这与观测到的地球物理数据（例如 USGS 报告的 2–6 mm/年）高度吻合，而之前的整数阶模型预测值在 $10^{-5}$ 范围内，显示出与现实世界的差异。
• 断层几何结构的影响：
  • 在情形 I（两条无限大断层）中，第二条断层附近的位移量大于第一条断层附近。
  • 在情形 II（无限大与有限长断层）中，有限长断层附近的位移量低于无限大断层附近。
  • 位移在断层迹线附近达到峰值，并随着距离增加而减小，最终趋于零。
• 分数阶 ($\alpha$) 的影响： 研究发现，应力积累和释放的速率高度依赖于 $\alpha$。较低的 $\alpha$ 值（表示较强的记忆效应）会导致更高的剪切应力积累量和更慢的松弛过程。随着 $\alpha$ 趋向于 1（接近整数阶行为），应力释放速率增加。
• 驱动力与应力分布： 跨越断层的总驱动力在蠕动开始后迅速释放。该释放量的大小随断层倾角和蠕动速度而变化。等值线图显示，应力积累在自由表面和断层位错带附近最为显著，而应力降低在对应于断层位置的深度处最为明显。
• 应变分析： 剪切应变 ($e_{12}$) 的变化范围在 $10^{-4}$ 量级，这与 2011 年东日本大地震（Tohoku-oki）等真实事件的数据一致。

意义与主张
作者声称，其分数阶模型提供了比传统整数阶模型更通用且更准确的断层相互作用动力学表征。具体而言：

• 分数阶方法成功捕捉了粘弹性介质中的记忆效应和非线性，而这些特性在标准模型中往往被忽视。
• 该模型的位移速率预测（$10^{-3}$）与观测到的地球物理数据具有强一致性，表明分数阶微积分是预测地震行为和地下变形的优越工具。
• 研究结果表明，断层间的相互作用对应力场的影响比单个断层更为广泛，且分数阶 $\alpha$ 是决定应力释放速率的关键参数。

论文得出结论，这些理论模型增强了对地下变形及其在影响断层动力学和地震发生中作用的理解，为分析地震活跃区的无震期提供了一个稳健的框架。