[论文解读] 3rd-order Spectral Representation Method: Part I -- Multi-dimensional random fields with fast Fourier transform implementation

Lohit Vandanapu, Michael D. Shields|arXiv (Cornell University)|Oct 14, 2019

Soil Geostatistics and Mapping参考文献 38被引用 5

一句话总结

本文提出了一种广义的三阶谱表示法（3rd-order SRM），通过同时引入功率谱和三谱约束，用于模拟具有非对称非线性的多维非高斯随机场。该方法利用快速傅里叶变换（FFT）实现高效计算，能够在标准硬件上实现可扩展的、高保真度的2D与3D场模拟，计算量随维度增加呈指数级减少。

ABSTRACT

This paper introduces a generalised 3rd-order Spectral Representation Method for the simulation of multi-dimensional stochastic fields with asymmetric non-linearities. The simulated random fields satisfy a prescribed Power Spectrum and Bispectrum. The general d-dimensional simulation equations are presented and the method is applied to simulate 2D and 3D random fields. The differences between samples generated by the proposed methodology and the existing classical Spectral Representation Method are analysed. An important feature of this methodology is that the formula can be implemented efficiently with the Fast Fourier Transform, details of which are presented. Computational savings are shown to grow exponentially with dimensionality as a testament of the scalability of the simulation methodology.

研究动机与目标

• 解决经典二阶方法在捕捉多维随机场中非高斯、非对称非线性特性方面的局限性。
• 实现具有指定功率谱与三谱的随机场的精确模拟，尤其适用于具有偏斜、非高斯行为的场。
• 通过快速傅里叶变换（FFT）实现计算高效的算法实现，使高维模拟在标准硬件上成为可能。
• 通过2D与3D域中的数值示例，展示该方法的可扩展性与保真度。
• 验证模拟场与目标的二阶与三阶统计特性（包括偏度与三谱特征）的一致性。

提出的方法

• 构建一个广义的d维随机展开，同时包含二阶（功率谱）与三阶（三谱）谱信息。
• 通过复值傅里叶模态与由功率谱和三谱导出的振幅的随机变量调制，定义三阶SRM。
• 利用快速傅里叶变换（FFT）实现模拟算法，以加速计算，将高维情况下的复杂度从O(N³)降低至O(N log N)。
• 通过使用相同的随机相位实现，确保模拟中相位的一致性，以便与二阶方法进行比较。
• 通过逆FFT对谱分量进行处理，构建随机场实现，同时正确处理共轭对称性，以确保输出为实值。
• 通过在波数空间中使用均匀网格对谱域进行离散化（Δκ₁=Δκ₂=Δκ₃=0.314），并采用N₁=N₂=N₃=16，M₁=M₂=M₃=32，将该方法应用于2D与3D场。

实验结果

研究问题

• RQ1 能否为d维随机场构建广义的三阶谱表示法，以捕捉非对称的非线性相互作用？
• RQ2 如何利用快速傅里叶变换（FFT）高效实现三阶SRM，以降低多维场景下的计算成本？
• RQ3 模拟场在多大程度上表现出由三谱决定的预期非高斯特征，如偏度与方向非对称性？
• RQ4 与经典二阶SRM生成的场相比，模拟场在结构与统计特征上存在哪些差异？
• RQ5 三阶SRM在维度上的可扩展性如何？它是否能够在标准台式计算机上实现3D随机场的实际模拟？

主要发现

• 三阶SRM成功生成了匹配目标功率谱与三谱的2D与3D随机场，三谱B₁、B₂与B₃对应的模拟偏度值分别为0.00680、0.00682与0.00661，与目标偏度0.00580非常接近。
• 三阶场的方差估计值（178.96–178.98）与目标方差（179.0812）的偏差在0.1%以内，证实了二阶矩的准确匹配。
• 三阶与二阶模拟场之间的差异揭示了与特定角度（如±73°与45°）对齐的非对称、拉长的特征，具体取决于各平面中三谱的各向异性。
• 场中非对称特征的尺度直接由三谱中的相对各向异性决定——例如，由于ω₁₁²的权重更高，B₁在x₁方向上的尺度更小。
• FFT实现带来的计算节省随维度增加呈指数增长，使得在台式计算机上实现16×16×16离散化的高效3D场生成成为可能。
• 尽管无法直观显示完整的三谱，但方法保持了正确的三谱结构，这由一致的偏度与矩统计量所证实。