用扩散模型加速地震波场模拟,速度提升2.17倍且可跳过小步长限制。
Generative wave propagator

- 基于条件扩散模型,利用历史波场和速度模型预测下一时刻波场。
- 训练时采用时序加权损失,稳定长序列递推,支持十倍于传统方法的物理时间步长。
- 适合需要快速波场模拟的地震反演、全波形反演等场景。
地震波场模拟是地震学的基础,但传统有限差分(FD)方法受限于数值色散和稳定性约束,常需密集空间网格和小时间步长,严重制约迭代反演流程效率。本文提出一种基于条件扩散的波场传播器,可递归地从一个时间步推进到下一个。模型以近期波场快照、速度模型和时间步索引为条件,学习相邻物理状态间的条件转移。通过直接训练网络预测干净的下一波场快照,强物理约束使得反向扩散过程可被单次网络评估替代。为提升长时间递推稳定性,引入因果时间加权损失:自适应权重以快照训练误差的指数移动平均累积,强调与前向传播一致的方向,抑制单步预测误差放大。由于学习的传播器绑定于训练快照的时间间隔而非FD稳定性极限,其可用物理时间步长比底层求解器大十倍。在Overthrust、SEG/EAGE和Marmousi模型上的实验表明,该方法能准确复现波场快照和炮集,且在相同硬件条件下实现2.17倍于GPU加速的十阶交错网格FD实现的端到端加速。
原文摘要 · Abstract (English)
Seismic wavefield simulation is fundamental to seismology, but conventional finite-difference (FD) methods remain limited by numerical dispersion and stability constraints, which often require dense spatial grids and small time steps and thereby severely limit the effectiveness of iterative inversion workflows. We introduce a conditional diffusion-based wavefield propagator that advances seismic wavefields recursively from one time step to the next. Instead of learning an unconditional data distribution of wavefield evolution, the model is conditioned by a short history of recent wavefield time steps (snapshots), the velocity model, and the wavefield time step index, allowing it to represent the conditional transition between adjacent physical states. By training the network to directly predict the clean next wavefield snapshot, this strong physical conditioning makes it possible to replace the iterative reverse diffusion process with a single network evaluation for each predicted snapshot. To improve stability over long recursive rollouts, we further introduce a causal time-weighted loss, in which adaptive weights, accumulated as exponential moving averages of per-snapshot training errors, emphasize training directions that are consistent with the forward propagation sequence and reduce the amplification of one-step prediction errors. Because the learned propagator is tied to the temporal spacing of the training snapshots rather than to the FD stability limit, it can advance the wavefield using a physical time step ten times larger than that required by the underlying solver. Experiments on the Overthrust, SEG/EAGE, and Marmousi models show that the proposed method accurately reproduces wavefield snapshots and shot gathers and achieves an end-to-end speedup of 2.17 x over a GPU-accelerated tenth-order staggered-grid FD implementation under matched hardware conditions.
Thank you to arXiv for use of its open access interoperability. PaperDance 不是 arXiv 官方产品;中文卡片由大模型生成,请以原文为准。