地震地面运动信号的标度分析

《Entropy》：Scaling Analysis of Seismic Ground Motion Signals

　　

地震在目录层级表现出公认的统计规律性，例如Gutenberg–Richter震级-频度关系与Omori–Utsu余震衰减，常被解释为地震活动作为远离热力学平衡的外驱动耗散系统的特征。然而，在单个地面运动记录中是否可直接检测到尺度不变性与反常扩散等可比特征仍属开放问题。本研究探讨2009年Mw 6.1 L'Aquila地震记录的加速度、速度与位移信号是否具备与非平衡复杂系统一致的统计性质，以及不同地震相是否携带可区分且可复现的统计特征。研究人员采用AR-AIC（自回归-赤池信息准则）估计P波与S波到时，以CRUST1.0速度模型理论到时为中心设自适应搜索窗；尾波（coda）起始采用三种互补判据的中位集成确定，将记录分割为至多五个时间窗。在各时间窗与信号类型下，于强反常扩散（strong anomalous diffusion）框架内分析位移矩标度，得到标度指数。研究人员针对尾波起始方法、尾波起止经验阈值、矩标度拟合所用τ子区间及地面运动信号滤波频带系统评估鲁棒性。反常扩散证据呈现细微结构：其符号与量级依赖地震相，仅少数配置在所有分割方案下保持稳定。结果表明，反常标度特征若存在并非普适，系统鲁棒性分析对区分真实物理效应与分割伪影至关重要。

研究背景方面，地震作为地壳非线性动力学的典型表现，在目录层级已揭示Gutenberg–Richter关系与Omori–Utsu衰减等尺度不变规律，被归入远离平衡的外驱动耗散系统。此类系统可出现重尾分布、矩的反常标度及强反常扩散（即第q阶矩标度指数非线性依赖于q）。尽管目录层统计已被充分研究，单个地面运动记录（加速度、速度、位移）是否携带与失衡动力学一致的反常标度特征，以及P波、S波、尾波、事件前噪声等相态是否具有可复现统计签名，仍缺乏系统回答。已有少数研究将多重分形技术用于地震波形，但未在相态分割前提下采用矩标度框架。为此，研究人员以2009年L'Aquila M_w 6.1地震为案例，利用意大利加速度归档库（ITACA）提供的192条三分量记录（经筛选后62台站186条），开展相态识别、信号分割与矩标度分析，并在《Entropy》发表，旨在检验地面运动信号中反常扩散签名的存在性及其对方法与相态的依赖性。

关键技术方法上，样本为ITACA归档的2009 L'Aquila地震三分量地面运动（加速度、速度、位移），采样率200 Hz，基线校正后用于增量矩计算。相态检测采用ObsPy的ar_pick实现AR-AIC，以CRUST1.0分层速度模型算理论P/S到时并设距离分箱自适应窗；尾波起始并行采用Rautian（2倍S走时）、Arias强度（95%累积能）、包络衰减（Hilbert变换后平滑包络降至S波峰30%）三法及其中位集成，尾波终止配对应准则（包络10%稳持2 s或Arias D99.5）。每记录分割为事件前、P波、S波、尾波、事件后至多五窗。矩标度分析按Vollmer等框架，对位移增量Δx计算q阶矩?|Δx(τ)|^q?∝τ^ζ(q)，q取-3至3步长0.15，τ对数取样，OLS回归得ζ(q)，以空间集网平均降低方差；并对尾波阈值、P波τ子区间、滤波频带做敏感性z-score与符号检验。

研究结果部分，3.1节相态检测与分割显示AR-AIC的P波残差标准差1.92–2.32 s，S波离散更大，位移SNR P违例率达26.5%，MonoP/MonoS违例率超82%源于真实速度横向不均而非纯拾取错；SPO台站图示五窗可用。3.2节矩标度与反常扩散中，P波窗三信号ζ(q)随q近似线性增长且明显高于正态扩散参考ζ=q/2，加速度ζ(3)≈2.5、位移ζ(3)≈3.0，但P波窗仅数秒致τ动态范围窄；S波窗加速度与速度ζ近零并转负，位移弱正ζ(1)但ζ(3)负，呈非线性；尾波窗加速度ζ负、速度ζ负，位移ζ(1)正但ζ(3)依尾波法变号（Rautian/Arias负，Envelope/Median正）。3.3节敏感性分析表明：尾波起止阈值变动除Envelope法20%起始在尾波窗z≈2.1外均|z|<2；P波τ上下子区间符号检验全信号拒绝随机（p<0.01），速度z达9说明单幂律不充分；滤波频带对P波窗影响最大（近全带z_max≈4.8），S/尾波窗稳健。

讨论部分指出，AR-AIC单调违例反映地壳横向非均质与射线路径效应，非拾取失效；位移因二次积分丢失高频时序信息致SNR P违例高。P波窗高ζ(q)虽指示偏离正态扩散，但受短窗与滤波敏感双重制约，应读作“存在偏离”而非精确定量；其ζ(q)随q增源于脉冲式到达与稀有大幅值波动，与AR-AIC能拾取P波同源。S波窗多子相叠加近正态扩散且对方法与滤波稳健。尾波窗散射能衰减近背景噪声，位移Arias法得最约束拟合（R²>0.9），Envelope法Peak违例最高说明其窗界与PGD位置最不一致。结论部分翻译如下：本研究探讨2009 M_w 6.1 L'Aquila地震加速度、速度、位移信号是否显示非平衡复杂系统统计签名及跨信号类型、地震相与方法选择的复现性。各记录经AR-AIC与四种尾波起始定义分割，矩标度分析并对尾波阈值、τ子区间、滤波带评估鲁棒性。反常扩散证据呈细微性：可靠性依地震相及尾波定义而异；P波窗指数最一致但最不受尾波相关选择影响却最受τ子区间与滤波敏感；S波窗负指数弱约束但对尾波法与滤波稳健；尾波窗最弱且方法依赖，位移常跨法变号。条件性是核心发现：地面运动反常标度主张不能视为方法无关，跨尾波定义系统比较可辨真伪。结果兼容地震活动为远平衡外驱动耗散系统解释，但未识别临界指数或普适类，不作相变断言。局限含尾波法依赖、AR-AIC残差不可分拾取误与速度模型误、单事件限制；未来需扩事件级数与台网几何。该研究提供了一套以鲁棒性框架鉴别地面运动标度伪影的方法论范例。