超越Hurst指数：谱前因子在断层粗糙度表征中的作用及其合成验证

《Journal of Structural Geology》：Beyond the Hurst Exponent: The Role of the Spectral Pre-factor in Fault Roughness Characterization and Its Synthetic Validation

　　

功率谱密度分析（power spectral density analysis, PSD）是表征断层表面粗糙度的标准方法。Hurst指数（H）常作为自仿射标度行为的主要描述符，即使具有相同H值的表面在粗糙度振幅和表观形态上表现出显著差异。利用合成表面，研究人员将定性观察（例如粗糙和光滑）与定量参数协调起来，并量化结果对三种谱平均变体（spectral averaging variants）和仪器噪声（instrumental noise）的敏感性。例如，滑移平行剖面（slip-parallel profiles）尽管具有较低H值却更光滑的观察，可通过同时考虑H和前因子（pre-factor, α）加以解释。恢复的参数处于周期图估计器（periodogram estimator）预期的统计变异性内，这是该方法本身的特征，而非数据质量的反映。带限仪器噪声（band-limited instrumental noise）导致H系统性低估，其严重程度在信号衰减更快（即H更高）的方向上最大，并在更低频率处沉降至噪声基底（noise floor）以下。因此，测得的H各向异性（anisotropy）被降低甚至反转。多尺度验证（multi-scale validation）表明，匹配α对于不同尺寸表面之间的谱对齐（spectral alignment）是必要的。在合成测试的指导下，研究人员将该分析应用于意大利中亚平宁（central Apennines）Roccapreturo断层的一个露头尺度斑块（outcrop-scale patch），该斑块使用SfM-MVS摄影测量（Structure from Motion Multi-View Stereo photogrammetry）获取，采用两个不同分辨率的相机系统和两条处理流程（processing pipelines）。两次采集在低空间频率处收敛，但在阈值以上发散，为一次采集的可用频率范围提供了内部判据（internal criterion）。

**研究背景与问题**
断层粗糙度表征滑面上凹凸体的尺寸和分布，这些凹凸体由不同尺度过程形成。在断层带尺度，不规则性可源自断层结构演化、非平面或弯曲传播路径、岩性-力学控制以及断层交汇；在单个凹凸体尺度，粗糙度主要归因于滑动岩体间摩擦阻力相关的磨蚀磨损，同震摩擦生热也可通过纳米尺度熔融和粒度减小产生镜面状滑移平滑。断层粗糙度已在微观到大陆尺度被描述，并具有尺度依赖形态。自仿射（self-affine）意味着沿平行和垂直滑移方向的形态属性不同。Hurst指数（H）常被作为主要描述符，但相同H值的表面可在粗糙度振幅和表观形态上显著不同。功率谱密度（PSD, power spectral density）方法可同时获得H和谱前因子（pre-factor, α），但以往研究常侧重H而忽视α；谱平均策略、拟合频带和仪器噪声也会影响PSD估计。长期观察显示，滑移平行剖面H较低却更光滑，这与单独用H解释相冲突。因此，研究人员利用已知参数的合成各向异性自仿射表面和自然断层SfM-MVS（Structure from Motion Multi-View Stereo，运动恢复结构-多视角立体摄影测量）模型，检验H与α联合表征、三种谱平均与拟合变体的敏感性、仪器噪声影响，以及跨尺度谱对齐和可用频段判据。

**主要技术方法**
研究人员生成已知H_⊥、H_∥和目标Rq或P2V的各向异性自仿射合成表面：在频率域对高斯白噪声施加可分离各向异性振幅函数，再经反傅里叶变换得到零均值高程场。随后注入带限高频噪声，控制截止波长λ_c、噪声RMS Rq_noise和信噪比（SNR, signal-to-noise ratio）。PSD分析对合成和自然表面沿X_⊥、Y_∥方向提取剖面，进行3% taper后FFT，计算功率谱密度，并用allmean、progmean、userbounded progmean三种平均/拟合变体估计H和α，计算理论RMS粗糙度Rq_th。自然样本为意大利中亚平宁Roccapreturo断层露头尺度斑块，2024年7月13日于42°11'46.75"N、13°42'7.47"E采集，使用SfM-MVS摄影测量，两台相机系统（iPhone 13 Pro和Nikon D5300）及两条Agisoft Metashape Pro处理流程。MATLAB流水线公开于GitHub。

**研究结果**
4.1 单尺度PSD验证（H_x=0.8, H_y=0.6）：30个0.5 m×0.5 m、点间距0.001 m、目标Rq=0.001 m的合成表面显示，三种变体给出的H和α存在系统性差异；userbounded progmean未显示统计显著偏差，而allmean和progmean在多数方向/变体组合中偏差显著。progmean的Rq_th接近目标Rq，验证了Parseval关系。无噪声条件下每频log₁₀P(k)标准差均值约0.489–0.505个数量级，接近周期图估计器理论预期。

4.2 跨组合频率范围的多尺度PSD验证：将30个小表面与50 m×50 m、dx=0.1 m、目标Rq=3.98×10^-2 m的大表面配对，并按Rq(L)∝L^H使α一致。合并谱约覆盖四个数量级频率。α在小网格和大网格拟合间无显著差异（p=0.984），而合并拟合未显著降低H偏差；H_large无偏，H_small与H_combined存在小但显著偏差。说明α一致性而非H保证跨尺度谱对齐。

4.3 仪器噪声对谱估计的影响：对30个表面加入λ_c=0.1 m、Rq_noise=0.0001 m、SNR=20 dB的带限噪声。三种变体均低估H，allmean最严重，userbounded progmean最稳健；H较高方向（X）低估更大，可缩小甚至反转表观各向异性。谱在高频平坦化，谱离散仍接近无噪声条件但略高。

4.4 H和α在控制表面粗糙度中的独立作用：四个各向同性表面组合H=0.90/0.50和Rq=0.001/0.010 m。相同H但不同Rq的表面内部结构相似而振幅差异大；相同Rq但不同H的表面振幅相近而小尺度不规则性不同，说明H和α独立控制表面粗糙度。

4.5 Roccapreturo断层的谱表征：高质量Nikon模型H_⊥=0.803、H_∥=0.693，α_⊥=3.58×10^-5 m³、α_∥=7.79×10^-6 m³；低质量iPhone模型H_⊥=0.589、H_∥=0.438，α_⊥=1.06×10^-5 m³、α_∥=1.67×10^-6 m³。两模型低频一致，超过阈值后低质量模型谱平坦化。高质量模型谱离散接近无噪声合成，低质量模型显著更高。两模型均保留α_⊥>α_∥和H_⊥>H_∥层级。

**讨论总结与结论翻译**
讨论部分指出，PSD图中谱估计的垂直离散是周期图估计器的固有特性，而非单纯数据质量差或地质复杂性标志；无噪声合成表面即约有0.48–0.51个数量级离散，理论期望约0.557个数量级，带噪声合成约0.54–0.55，高质量Nikon模型约0.561–0.566，低质量iPhone模型约0.687–0.730，后者指示较强仪器噪声。H和α共同指定平稳高斯自仿射场的二阶统计；仅报告H会混淆尺度与振幅信息，而滑移方向H和α均较低可解释其更光滑。H接近或超过1可能反映大尺度确定性结构，不宜直接解释为物理H值。三种拟合变体引入系统差异，可能被误认为真实地质变化；带限噪声使高频谱平坦，H低估在较高H方向更严重，各向异性可缩小或反转，α变化幅度大于H。多尺度合并需α一致。合成高斯场与自然断层存在差异，但PSD管道估计H和α的可靠性可验证；周期图统计性质可转移，偏差幅度不一定转移。大尺度波状起伏可抬高X方向低频功率，可能使H_⊥偏高。Nikon/iPhone低频收敛显示采集内一致性，阈值以上发散给出可用频段判据，但不构成几何精度外部验证。局限包括合成实现、参数组合和自然表面有限，频带固定；未来可自适应选频带、二维谱分析、更高阶统计和空间相关函数，并将α或Rq_th与接触面积、孔隙度、渗透率等联系。

结论翻译：1. 仅H不足：来自同一回归的谱前因子α必须与H一同报告。具有相同H的两个表面可因α不同而在振幅上相差一个数量级；一旦同时考虑α，长期观察到的断层表面沿滑移方向尽管H较低却看似更光滑的现象即可得到解释，因为沿滑移方向H和α均趋于较低。2. 合成验证表明，H的恢复处于周期图估计器预期的统计离散范围内（无噪声条件下log₁₀P(k)约0.48–0.51个数量级），这是估计器本身的性质，而非数据质量差的标志。3. 带限仪器噪声引入了方向依赖的H低估，在真实H较高的方向上更严重，可使表面的表观各向异性缩小甚至反转；在所测试的三种拟合变体中，userbounded progmean对该偏差最稳健，allmean最敏感。4. 合并物理尺寸差异很大的表面证实，允许其谱合并的是α的一致性，而非H。5. 将该流程应用于Roccapreturo断层时，相同流程在真实数据上再现了这些模式。同一斑块的高质量和低质量SfM-MVS模型尽管裁剪范围不同，但在低空间频率处一致；它们在由低质量模型噪声基底设定的阈值以上发散。该收敛与发散点提供了一种直接从PSD判断SfM-MVS采集可靠频率范围的实用方法。两个模型均保留了H_⊥>H_∥的层级。对于使用此类分析的人员，有三点结论：将H和α一同报告，而非仅报告H；说明使用了哪种谱平均与拟合策略，因为它可测量地改变两个参数；并将每个频率处PSD估计值的离散视为方法属性，而非丢弃剖面或频率范围的理由。这些结论仍基于有限的合成实现、单一自然案例研究以及一个固定而非自适应选择的频带。相同框架应扩展到其他断层、仪器和尺度。