赞
踩
1、inversionnet_train_light 试运行——未成功
2、DL-FWI基础入门培训-1,2,以及作业1的完成——暂未完成作业
在补充完相关依赖后,运行仍有报错
产生原因:这个代码在当时师兄的库中存在相关申明,但现在没有了
解决方法:在 “InversionNet.py” 中补充以下代码
import torch.nn.functional as F
关于为什么解决完所有问题后仍然运行不成功:内存不够!
关于代码还需要说明的一点:目前只是将代码运行成功,实际上并没有逐句分析,在后面学习了相关知识后,建议开个篇章逐句理解。
学习目的:进一步认识 “通过全波形反演探测地震数据” 这一系统。
内容摘要:
(1)“全波形反演” 概述
A. 关于正演与反演的逻辑
B. FWI的数学表示
C. FWI的端到端深度学习解释
D. DL-FWI的优势
E. DL-FWI的基本发展和分类
F. 地震数据 (Seismic Data) 是什么样的?
G. 速度模型 (Velocity Model) 是什么?
H. 道 (Trace) 是什么?
笔记记录:
(一)何为全波形反演 (Full waveform inversion)?
从地表发出的震源弹性波会沿着不同方向传播, 并受到不同介质特性的影响而发生反射、折射、衍射等现象。
我们可以在地表放置多个接收器来检测这些反射信号, 从而来分析关于地下介质结构及物性参数等信息。当我们把每个接收器获取的弹性波绘制出来并拼接在一起就可以得到如下的地震数据。
而通过右图的结构推测左图构造的问题就称之为全波形反演, 后续我们简称为FWI.
A. 关于正演与反演的逻辑
正演 (Forward modeling) 指的是根据地下介质的物性参数, 运用物理方程和数值计算方法, 模拟真实地震波传播过程并得到合成地震记录的过程. 可以认为正演是反演的逆过程: 即通过地层图模拟地震波形.
反演问题是建立从地震数据到速度模型的映射. 一般认为是多解的, 即一个地震数据可能对应于多个速度模型.
正演问题是建立从速度模型到地震数据的映射.一般认为是单解的, 即一个速度模型只能生成一个地震数据.
B. FWI的数学表示
通过之前的表示, 常密度二维声波方程式:
若我们把上式简写为:
这个公式即将映射到
, 这里通常是非线性的, 这就是常规的正演思路.
在反演应用中, 求解上述方程的策略是开发一个相当快速且准确的逆算子.
现有的采用PDE (偏微分方程)的全波形反演方法将会使用迭代优化算法来最小化目标函数.
但是由于算子的非线性特性和测量
本身可能得不完善 (因为
本身也是通过数值化的正演模拟得到) 很难得到精确地下模型.
(吐槽一下这和电磁波有什么区别啊可恶,死去的电磁场与电磁波还在攻击我)
注:这部分的数学理论部分理解可参考 “数学物理方法” 和 “电磁场与电磁波” 等相关课程,计算靠机器。
C. FWI的端到端深度学习解释
而随着深度学习的壮大, 我们以深度网络的逼近定理为原则, 尝试用各种端到端结构去模拟目标函数.
由此FWI的端到端深度学习的方法可以总结为:
就是用已知的”地震数据”和”速度模型”来训练一个[模型].
然后通过将新的”地震数据”喂入这个[模型]中, 从而预测其”速度模型”, 从而推测地下构造.
D. DL-FWI的优势
传统的FWI的求解过程涉及到多次迭代, 时间开销是巨大的, 而且依赖于初始解, 存在病态性.
1. DL-FWI通常在训练阶段需要更多时间, 而实际预测花费的时间很少, 只需要几秒钟. 一旦建立了良好的通用网络, 那么就可以大大地减少FWI计算时间.
2. 神经网络可以自然地将空间相关性作为正则化引入生成速度模型, 这样可以抑制梯度中的噪声并减轻局部最小.
E. DL-FWI的基本发展和分类
DL-FWI最主要的发展是从2018年开始, 而2020左右是相关的重要架构产出的时代.
最开始人们采用的手段主要是”严格端到端”和”空间信息的嵌入”两种手段.
严格的端到端手法就是尝试通过设计构建完善的网络去拟合函数
Copyright © 2003-2013 www.wpsshop.cn 版权所有,并保留所有权利。