历史与路径依赖14 分钟

同一断层为什么有时缓慢滑动,有时快速破裂?

断层是否快速破裂,取决于加载能提供多少驱动力、界面怎样逐步失去承载能力,以及滑动已经改变了什么。近年的透明材料实验把持续慢破裂、慢成核后加速和慢破裂尾部另起快破裂区分开来。初始局部扰动与背景应力的贡献可以此消彼长,而过程区、临界长度和界面宽度的关系决定哪些近似适用;这些机制线索尚不能直接变成天然地震预报。

“这条断层滑得慢”,究竟是哪一种慢

两块接触的固体被持续推挤,有时能一点一点错动,有时却突然大范围失稳。天然断层也有持续时间和辐射表现差别很大的滑动事件。若只给一条断层贴上“稳定”或“不稳定”的标签,很容易漏掉事件发生前后不断变化的条件。

先把三个速度分开。加载速度是外界怎样缓慢推动系统;滑动速度是界面两侧在某个位置怎样相对错动;破裂前沿速度则是“开始滑动或发生明显弱化”的边界沿界面推进多快。一处界面只错动很短距离,并不妨碍前沿在很短时间内扫过一大片区域。三者不能互换。1

“慢”还包含两类不同过程:一类前沿在较长距离上持续缓慢推进;另一类是快速事件之前的缓慢成核阶段。后者讨论从局部滑动区到大范围失稳的过渡,前者未必会成为快速事件的前奏。

本篇用透明材料摩擦界面的原始实验,把几个环节拆开:外界加载怎样提供驱动力,界面接触怎样损坏和恢复,已有的小滑动区怎样继续长大,以及什么时候必须考虑前沿附近并非一个无限小的点。最后再讨论这些线索能带到天然断层多远。实验室界面是看清过程的窗口,不能替代地壳中的复杂断层。

滑动总量相同,过程可以相差十万倍

设两个过程最后都相对错动100微米。一个用10秒完成,平均滑速为每秒10微米;另一个用0.1毫秒完成,平均滑速为每秒1米。累计滑动完全一样,时间却相差十万倍。

自选平滑滑速脉冲的积分,两条曲线最后都累积100微米。左右横轴分别为秒和毫秒,形状并非实验拟合;这里比较局部相对滑动,不是破裂前沿走过的距离。
自选平滑滑速脉冲的积分,两条曲线最后都累积100微米。左右横轴分别为秒和毫秒,形状并非实验拟合;这里比较局部相对滑动,不是破裂前沿走过的距离。

这不是为了证明慢滑动和地震一样,而是提醒:总共移动多少与移动得多突然,是两项不同的信息。 在同样的滑动面积和剪切模量下,二者可以具有同样的标量地震矩:

$$ M_0=\mu A\bar d. $$

这里 $\mu$ 是剪切模量,$A$ 是滑动面积,$\bar d$ 是平均滑移。这个关系没有给出辐射地震波的能量。辐射还取决于破裂过程、应力变化、速度时程等条件;同样的矩也不等于同样的地面摇晃。1

因此,“几乎感觉不到”不能直接解释成“没有发生显著滑动”;反过来,“释放了相当的滑动量”也不能直接翻译成“消除了某次未来地震”。后一个判断还需要知道周围应力被怎样重新分配。

一个前沿能不能前进,要看驱动力与代价一起怎样变

摩擦界面并不是两个处处完整贴合的理想平面。真实承载常集中在一批接触点上。前沿经过时,接触点逐步失效、重排或进入新的滑动状态;周围固体原先储存的弹性能随之释放。

断裂力学把其中一个问题写成:前沿每新扫过单位面积,能够供给的能量 $G$,能否支付前沿局部破坏和弱化的代价 $\Gamma_{\rm tip}$。在合适的线弹性、局部过程区和残余摩擦处理下,传播条件可写为:

$$ G(c_f,\ell,\text{加载与几何})=\Gamma_{\rm tip}. $$

$c_f$ 是前沿速度,$\ell$ 表示相关破裂尺度。这里的断裂能不是“整个界面所有摩擦发热”的别名。已经滑开的区域继续摩擦所耗散的能量,与前沿形成新滑动区域的局部代价,需要按照同一口径分别处理。1

这也解释了为什么只看一个峰值摩擦系数不够。界面从承载到滑动,需要经历多长的弱化过程?残余阻力怎样随速度变化?之前滑过的地方还剩多少接触?外界此时继续加载,还是能量主要来自原先储存的弹性变形?这些条件都会改变传播。

常见的速度—状态摩擦描述,还用一个随接触停留和滑移更新的“状态变量”概括界面的记忆。停留时间长,接触可能逐渐增强;一旦滑动,旧接触被更新。状态变量是一个模型化的记账方式,不代表所有材料的记忆都来自同一种微观过程。4

速度增强有时能提供抑制加速的反馈,速度弱化有时能促进失稳,但它们也不是脱离加载刚度、尺度、状态演化和几何就能使用的全局标签。下面的实验恰好展示了这些条件的具体作用。

同一个界面,慢事件已经改变了后续快事件的条件

Shi和Fineberg在2025年研究PVC接触界面,记录到同一界面上的持续慢破裂与快速破裂。其补充材料用接触光强和数字图像相关互相对照,支持“接触面积明显下降”与局部滑动开始的对应。材料身份不变,传播方式却可以改变。2

这个结果的意义是:不能把两种事件的全部差异,都归因于两种永远不同的材料。它并没有证明界面的实际状态保持不变。

出版方图3至图5提供了更具体的区分:残余剪切强度在所测的低速段与高速段呈现不同速度依赖;接触面积也随停留时间和此前破裂而变化。作者据此建立速度相关的阻力与断裂能描述,比较相应传播解。这些曲线中的速度是论文所定义的破裂速度,不能直接当作任意位置的局部滑动速度。3

一个特别有用的观察是:快破裂可以在慢破裂经过后的尾部重新成核。 这时不是必须想象同一个慢前沿一路平滑加速;先前的慢过程已经改变了一段界面的接触与承载条件,新的快速事件在这片区域启动。

图4还比较了加载速度及加载暂停后恢复的情况。它们给出了历史和加载方式能够改变后续传播的实验入口,但不是把界面所有状态变量逐个固定后只改一个分子的干预。作者模型对两种传播方式的解释,应与观测事实分别看待。

本篇已核该研究的补充材料、出版方完整图注及作者讲演资料,未取得订阅主文全文和原始数据包。因此这里采用能回查的观测与模型解释,不重建其全部定量拟合。尤其是图5明确保留了中间一段速度范围缺少残余强度测量的空位,不能把两端测量之间画出的机制图当成整段已经测全。

一次小扰动,为什么能改变后续要等多久

另一个问题是:局部成核已经开始,为什么有的事件很快进入快速破裂,有的却要经过较长的慢阶段?

Fryer等人在2026年的PMMA界面实验中,比较了同一实验内三次事件。作者报告,事件前测得的应力分布较为接近,但成核时间分别约为12.7、3.6和0.6毫秒。这说明“当时测得的背景应力相近”不足以锁定同样的时间过程。4

他们追踪到的线索,是起始小事件留下的瞬态滑动:局部滑速先因扰动上升,再降低到一个较小值,随后可能重新加速。作者把这一中间阶段的局部瞬态滑速记作 $V_{\min}$。它和后续成核时间、空间尺度有关系。

这里有两项测量条件必须一起读:成核结束由前沿速度超过每秒10米来操作性定义,这不是所有材料的普适分界;很小的 $V_{\min}$ 则来自对滑移时程拟合后求导的估计,不能写成仪器直接读出的瞬时速度。

Fryer等2026年Fig4c公开表重绘。64条实验记录中,47条在所用三项量上均为正,可进入本双对数图;其余17条完整保留在下载数据中,未取绝对值或补值。点的颜色表示原表成核长度。三条线由论文近似关系独立计算,是模型参照,不是本库新拟合。
Fryer等2026年Fig4c公开表重绘。64条实验记录中,47条在所用三项量上均为正,可进入本双对数图;其余17条完整保留在下载数据中,未取绝对值或补值。点的颜色表示原表成核长度。三条线由论文近似关系独立计算,是模型参照,不是本库新拟合。

原表中另17条包含10条全零记录和7条负的滑速记录。数据表没有逐条说明这些值的物理意义,本篇不把它们擅自归为“没有慢阶段”,也不把反复加载同一装置的事件视为独立天然断层样本。5

模型中,状态更新需要累积一定的滑移距离。若滑速更小,走过这段距离通常更费时间;但空间扩展和摩擦演化会修正简单的“距离除以速度”。原文Methods式17—18给出近似关系,包含速度相关修正,并用有效平均滑速连接局部估计值。本库按图4c参数重新计算了三条参照线,与所附3,000个表值的最大相对差约为 $1.19\times10^{-5}$。

这里的0.03、0.1和0.3微米,是状态演化的特征滑移距离,不是成核区的长度;后者在图中以厘米计。复算参照式也不等于重跑完整的摩擦运动方程或重新验证实验拟合。

走完一个最小模型:停下来、慢一阵、一直较快

怎样理解“同样的背景,起始扰动不同,后面就可能很不一样”?原研究将背景加载和起始局部滑动对前沿的作用分开:随着滑动区扩大,前者的贡献可以增加,后者的贡献则逐渐减弱。4

我们保留这两个相反趋势,另选一个简单模型实际计算。它从已经存在的小破裂区出发,不模拟接触从完好到形成尖端的全部过程。

令无量纲破裂半长为 $\ell$,前沿驱动的应力强度因子为:

$$ K_{\rm drive}=q\sqrt{\ell}+\frac{I}{\sqrt{\ell}}. $$

$q$ 表示固定背景项,$I$ 表示初始局部扰动留下的贡献。应力强度因子概括尖端附近应力集中的强弱,它不是能量本身。第一项随着区域变大而上升,第二项则下降,两项之和会经历一个低谷。

再选一个随前沿速度增加而提高的有效韧度,用简化的亚瑞利传播近似连接供给与代价:

$$ K_{\rm drive}\sqrt{1-u} =1+\eta\ln(1+u/u_0),\qquad 0\le u<1. $$

其中 $u=c_f/c_R$,$c_R$ 为材料的瑞利波速;取 $\eta=0.05$、$u_0=0.001$。这些都是解释机制的自选量,不是从PVC或天然断层拟合出的常数。若 $K_{\rm drive}\le1$,模型令前沿停止;否则求出速度,并积分 $d\ell/dt=u$。长度和时间均为自选无量纲单位。

自选模型,非原论文完整运动方程。左图放大驱动低谷附近,初期较高值超出纵轴显示范围;虚线为零速门槛。右图为从同一初始半长0.01出发的实际轨迹。灰色驱动曲线在更大长度处虽然再次上升,前沿已在第一次过门槛处停住,不能自动跨过中间区域。
自选模型,非原论文完整运动方程。左图放大驱动低谷附近,初期较高值超出纵轴显示范围;虚线为零速门槛。右图为从同一初始半长0.01出发的实际轨迹。灰色驱动曲线在更大长度处虽然再次上升,前沿已在第一次过门槛处停住,不能自动跨过中间区域。

固定 $q=0.6$,只改 $I$,就得到三种路径:

起始扰动贡献 计算过程 到同一终点的结果
$I=0.2$ 局部扰动衰减后,背景项还接不上,驱动落到门槛 趋近半长0.054处的停滞位置,达不到终点
$I=0.45$ 驱动低谷仍略高于门槛,前沿大幅减速后恢复 到半长3约需650个时间单位
$I=1$ 低谷离门槛较远,整个过程仍较快 到同一半长约需8.13个时间单位

为什么低谷如此关键?对这个自选式求极小值得:

$$ \ell_{\min}=\frac{I}{q},\qquad K_{\min}=2\sqrt{qI}. $$

中等扰动的低谷约为1.039,只比门槛1高一点,对应的最低前沿速度约为 $0.00117c_R$;较大扰动对应的最低速度约为 $0.310c_R$。于是,背景条件没变,从初始扰动过渡到背景驱动的那段过程,就能拉开很大的持续时间差。

停滞例子的曲线只画到距解析停滞位置约百万分之一处;图上截断时间不代表系统在那个有限时刻严格达到零速。数值积分已与独立求积交叉检查。

这个例子解释了一个可能的路径分叉,没有证明“真实断层只要知道一个扰动强度就能预报”。它未纳入有限宽度、过程区生成、接触恢复、流体或离断层损伤,速度修正也不适用于超剪切破裂。较大扰动能否越过其他障碍,要回到具体系统。

小破裂区还没有“长成熟”时,几何不能单独作决定

前面用了一个近似清楚的尖端。但真实接触不是在一个数学点突然全部失效:尖端附近有一段区域仍在逐步弱化,这个区域的长度称为过程区或黏聚长度 $\xi$。它若与界面宽度差不多,甚至更大,尖端近似就要谨慎使用。

Shi和Fineberg在2026年的聚碳酸酯PC实验中,直接比较了界面宽度、过程区和临界破裂尺度。这里材料是PC,不是前面的PVC或PMMA。6

他们看到,局部滑动区先在界面内二维扩大;当它跨满有限宽度、进入近似一维扩展后,仍可保持慢传播,随后才发生快速加速。跨满宽度这个几何变化,不总是快速破裂立刻开始的时刻。

解释这一差别,需要同时比较三把尺:逐步弱化需要多大的过程区,当前加载下多大的破裂才有足够驱动力,以及装置提供了多宽的界面。该研究把它与先前较脆材料的实验区间对照,说明几何主导和过程区主导可以出现在不同尺度关系中,不能把某个装置里的顺序推广为唯一顺序。

在特定二维中心裂纹、线弹性和均匀应力降等假设下,Griffith尺度可写成:

$$ \ell_G=\frac{E\Gamma}{\pi(\Delta\tau)^2}. $$

这里 $\ell_G$ 是半长,$E$ 为相应弹性参数,$\Delta\tau$ 为剪切应力降。改变裂纹几何或弹性约束时,需要核对系数与模量口径。

即便这个式子适用,“更难破坏,所以成核区必然更大”仍需加条件。固定其他量时,提高断裂能会增大临界长度;但停留更久或法向载荷变化,可能同时改变断裂能和应力降。若两者都增加一倍,按这个式子临界长度反而减半。一个材料处理不是只转动一枚互不牵连的旋钮。

带回天然断层:能解释哪些变化,还缺哪些证据

这些研究共同推进的,是把“同一断层”从一块永远不变的材料,改写为一个有历史的力学系统:接触会恢复和损坏,先前滑动会改写残余阻力,局部扰动会改变早期扩展过程,而几何和过程区决定哪些前沿近似可用。

但几组实验并不是相互独立地证明同一个完整模型。2025年PVC研究重点在同界面持续慢破裂与后续快破裂;2026年PMMA研究重点在瞬态滑动与成核时间;2026年PC研究重点在过程区、临界尺度和有限宽度的先后关系。它们各自钉住不同环节。

要进一步区分机制,值得看的是具体条件变化:暂停加载是否改变接触恢复和下一次传播;同样背景下起始局部滑动的差别能否说明时间差;改变宽度后快速加速是否始终紧随跨宽发生;慢过程结束后,后续快速事件是在原前沿继续推进,还是另起一个核。这些比事后只比较“都出现过慢和快”更有判别力。

天然断层还多了流体压力、粗糙度、温度、非均匀应力与界外变形;地表仪器往往只给出经过空间平均和模型反演的滑动线索。实验室透明界面中可直接分开的环节,在地下可能混在同一个观测量里。2026年研究对若干天然前震—慢滑动序列作了比较,但个案关联、反演量和事后筛选都不能自动变成前瞻预报能力。4

因此,更可靠的理解是:慢与快是加载、阻力演化和空间扩展共同产生的过程状态。 慢滑动可以持续,也可以停滞;它可能改变后续快速事件的条件,但不保证每次都会引出快事件。快速事件也不保证在现有仪器可分辨的时空范围内,先出现足够长、足够清楚的慢阶段。

本文图表和计算可回查:计算结果、全部64条原表记录及绘图标记、原式与参照表逐点比较、最小模型轨迹。正文中的合成例子与原始实验始终分别标明。

来源与依据

  1. S1 领域Perspective与理论/观测边界综合

    Earthquake energy dissipation in a fracture mechanics framework

    David S. Kammer, Gregory C. McLaskey, Rachel E. Abercrombie, Jean-Paul Ampuero et al. · 2024

    Nature Communications15;PMC11148115公开全文 · 相关全文已核

    Fundamentals及Box1–3;tip/tail能量、LEFM条件、观测量与自然应用限制

    核G与Gamma口径、残余摩擦、breakdown work与局部断裂能不自动相等,矩M0=µAd不决定辐射能。速度修正为指定几何、加载和亚瑞利近似,不作普适公式;本库K模型为自选说明。

  2. S2 原始PVC摩擦实验的补充材料;主文只取摘要

    Concurrent slow and fast frictional ruptures in laboratory earthquakes — Supplementary Information

    Songlin Shi, Jay Fineberg · 2025

    Nature Physics21,1019–1025;2025-04-14在线;原始SI五注已读 · 相关全文已核

    Contact imaging与DIC、Strain/contact/Slip velocity、aging/running-in及不同速度下材料响应等五注;主文摘要另存shi2025.html

    主文全文与Figshare数据未得;SI核接触光强与局部滑移联系、近稳态形状下Vp/VL约1/Δε、慢/快累计滑移例子。本篇不使用SI关于辐射能与矩成比例的泛化句;SI的µ=E/(1+ν)与标准各向同性关系相差2,不据此重算应力,也无证据判定该排印问题是否进入作者原分析。

  3. S3 s2同一原始研究的完整出版图注

    Concurrent slow and fast frictional ruptures — published figure captions 3–5

    Songlin Shi, Jay Fineberg · 2025

    Nature Physics正式图3、4、5;同一论文不算三个独立实验 · 相关全文已核

    Fig4a–c加载速度/暂停恢复和慢尾新快核;另存shi2025-figure3.html的接触/残余强度/aging、figure5.html的模型解和缺测中间速度区

    图4c明确新快裂纹在慢破裂尾部成核,不能写成同一前沿必然连续加速。图5模型采用先前滑动后的Gamma变化,非直接独立操纵Gamma;低高速之间部分tau_res(v)未测。作者讲演fineberg2025-slides.pdf提供同源解释,不当独立验证。

  4. S4 原始PMMA实验、速度状态Griffith型模型及天然个案比较

    Foreshock-induced slip transients set mainshock nucleation timing

    Barnaby Fryer, Dmitry Garagash, Mathias Lebihain, François Passelègue · 2026

    Nature653,752–757,2026-05-06在线,05-21卷期;正式全文PDF含Methods/Extended Data · 相关全文已核

    主文Two decades/Weakening/Foreshocks/EoM/Transients;Fig1–4;Methods measurement、Vmin拟合求导、EoM与式17–18参数及自然外推

    十实验94动态事件不等于Fig4c64行;压力表bar不当局部界面法向应力。12.7/3.6/.6ms取同一次实验;10m/s为操作性结束门槛。Vmin来自拟合时程导数,模型有效平均约取其2倍。式17的L是状态滑移距离而非裂纹半长。Fig4c a/b=.55、V0=1e-6与全模型拟合a=.016,b=.029,V0=1.35e-6并非同组参数。Methodsplane-stress标签与所列Hooke矩阵约束口径有不一致,未用该矩阵独立计算绝对应力;不据此否定整项研究。自然量经反演/代理并有个案关联争议,不证明前瞻可预测性。

  5. S5 s4作者公开表格汇总与参照曲线,非独立实验

    Source Data for Fryer et al. Fig. 4

    Barnaby Fryer et al. · 2026

    出版附件MOESM5;Fig4c表完整读取和复算 · 相关全文已核

    Fig4c1001×7;64行实验记录、1000×3模型参照;Fig4a/b/d只读结构与局部字段

    47行Vmin/时间/长度均正可双对数图,17行原样保存,缺每行零负值物理分类及实验块ID。原始传感器信号未重分析;fryer-fig3.xlsx及fryer-extended6.xlsx也取得,前者核结构,后者只读表头,不声称逐曲线复算完整EoM。

  6. S6 原始PC界面成核与接触/滑移成像实验、尺度比较

    Interplay of Cohesive, Griffith, and Geometric Scales in the Nucleation of Friction

    Songlin Shi, Jay Fineberg · 2026

    Physical Review Letters136,076203,2026-02-20正式版;APS harvest公开PDF · 相关全文已核

    主文pp1–5、Fig1–4;W/xi/Griffith尺度、中心裂纹半长与1D全长、cohesive建立及PMMA尺度对照

    PC不是PVC或PMMA;跨W后仍可慢,不推广所有自然地震必有可探测前兆。Gamma通过接触面积下降校准推断,非本库独立测量;另附SI未得。xi为渐进弱化的有效长度,具体阈值定义不能写成所有接触完全断开。Gamma与Δtau共变,固定参数偏导不冒充实际干预的总效应。

最后更新 2026-09-09

来源