超分辨率地震速度模型

news2025/1/19 8:26:13

文献分享


1. Multitask Learning for Super-Resolution

原题目:Multitask Learning for Super-Resolution of Seismic Velocity Model

全波形反演(FWI)是估算地下速度模型的强大工具。与传统反演策略相比,FWI充分利用了地震波的运动学和动力特性,具有更高的精度和分辨率。近年来,低频和中频波数的FWI发展迅速,而高频FWI由于其计算成本巨大,亟待解决。由于深度学习技术在各种地球物理问题中取得了显著的性能,作者建议使用深度学习方法来提高FWI的效率和准确性。

Core parts:
1.super-resolution:目标是从给定的低分辨率图像中获得视觉上令人愉悦的高分辨率图像。作者专注于使用SR技术来执行FWI任务,即使用SR作为后处理操作来提高FWI的质量。
2.multitask learning:把多个相关(related)的任务放在一起学习,同时学习多个任务。通过在一定程度上提高辅助任务的性能和共享参数,可以更好地泛化原始任务。

1.1 方法

在这里插入图片描述
具有硬参数共享的MTL:边界恢复是当前SR处理的难点,并且这在地震速度模型的SR中尤为突出。在地球物理学中,地震速度模型的边缘信息是地质模型的重要组成部分,这证明了使用边缘图像重建作为辅助任务的相关性。并且由于地震速度模型与其边缘图像高度相关,因此采用硬参数共享策略更加高效。文章得到的边缘图像是原始图像经过sobel算子卷积的结果。

边缘图像:是对原始图像进行边缘提取后得到的图像。边缘是图像性区域和另一个属性区域的交接处,是区域属性发生突变的地方,是图像中不确定性最大的地方,也是图像信息最集中的地方,图像的边缘包含着丰富的信息。

混合损失函数:由MSE和TV组成,表达式如下:
loss ( Φ ) = mse ( Φ ) + TV ( Φ ) \text{loss}(\Phi)= \text{mse}(\Phi) +\text{TV}(\Phi) loss(Φ)=mse(Φ)+TV(Φ)
mse ( Φ ) = 1 N M ∑ k = 1 3 α k ∑ i = 1 N ∑ j = 1 M ( f ( I L , Φ ) i , j , k , − I H i , j , k ) 2 , \text{mse}(\Phi) =\frac{1}{NM}\sum^{3}_{k=1}\alpha_k\sum^{N}_{i=1}\sum^{M}_{j=1}(f(I_L,\Phi)^{i,j,k},-I_H^{i,j,k})^2, mse(Φ)=NM1k=13αki=1Nj=1M(f(IL,Φ)i,j,k,IHi,j,k)2, where ∑ k = 1 3 α k = 1 \sum^{3}_{k=1}\alpha_k=1 k=13αk=1 and a k l 1 ( I H − I L ) ∣ c = k ≡ C , k = 1 , 2 , 3 a_k l_1(I_H-I_L)|_{c=k}\equiv C, k=1,2,3 akl1(IHIL)c=kC,k=1,2,3

TV ( Φ ) = β 1 1 K 1 ∑ i = 1 N ∑ j = 2 M ∣ f ( I L , Φ ) i , j , 1 − f ( I L , Φ ) i , j − 1 , 1 ∣ +   β 2 1 K 2 ∑ i = 1 N ∑ j = 2 M ∣ f ( I L , Φ ) i , j , 1 − f ( I L , Φ ) i − 1 , j , 1 ∣ , \text{TV}(\Phi)=\beta_1\frac {1}{K_1}\sum^{N}_{i=1}\sum^{M}_{j=2}\lvert f(I_L,\Phi)^{i,j,1}-f(I_L,\Phi)^{i,j-1,1}\rvert+ \ \beta_2\frac {1}{K_2}\sum^{N}_{i=1}\sum^{M}_{j=2}\lvert f(I_L,\Phi)^{i,j,1}-f(I_L,\Phi)^{i-1,j,1}\rvert, TV(Φ)=β1K11i=1Nj=2Mf(IL,Φ)i,j,1f(IL,Φ)i,j1,1+ β2K21i=1Nj=2Mf(IL,Φ)i,j,1f(IL,Φ)i1,j,1,
where K 1 = N ( M − 1 ) , K 2 = ( N − 1 ) M K_1 = N(M-1), K_2=(N-1)M K1=N(M1),K2=(N1)M.

在这里插入图片描述
同时,在训练过程中,采用了学习率衰减方法。实际学习率表示为
η i = a i η 0 \eta_i = a^i\eta_0 ηi=aiη0
η 0 \eta_0 η0为初始学习率, a a a为衰减率。

1.2 实验结果

在这里插入图片描述
在这里插入图片描述
在这里插入图片描述
在这里插入图片描述

2. Super-resolution guided by seismic data

原题目:Super-Resolution of Seismic Velocity ModelGuided by Seismic Data
前面的提到的,M-RUDSR的多任务学习框架已经成功通过提高地震速度模型的分辨率,来提高FWI结果的精度。
然而,作者认为它并没有充分利用地震数据。但由于地震速度模型和地震数据处于不同的频段,单纯增加模型输入输出通道来实现地震数据的利用效果有限。在M-RUSDR的基础上,提出了M-RUSDRv2。
不同的是,M-RUSDR是仅将地震速度模型边缘图像的SR作为辅助任务以增强地震速度模型的分辨率,而M-RUSDRv2不仅考虑地震速度模型及其边缘图像,还考虑了地震数据及其边缘图像。并且M-RUSDRv2通过在特定数据的训练,海量数据上的微调,实现了较强的泛化能力。

2.1 方法

在这里插入图片描述
损失函数:与上文仅仅有一处不同,即 a k a_k ak
∑ k = 1 K a k = 1 , K = 6 \sum_{k=1}^K a_k=1,K=6 k=1Kak=1,K=6 and a k = m a k + 3 , m > 0 , k ∈ [ 1 , 3 ] a_k=ma_{k+3},m>0,k \in[1,3] ak=mak+3,m>0,k[1,3] and a k 1 ∥ I H k 1 − I L k 2 ∥ 1 = a k 2 ∥ I H k 2 − I L k 2 ∥ 1 , k 1 , k 2 ∈ [ 1 , 3 ] a_{k_1}\lVert I_H^{k_1}-I_{L}^{k_2}\rVert_1=a_{k_2}\lVert I_H^{k_2}-I_{L}^{k_2}\rVert_1, k_1,k_2 \in[1,3] ak1IHk1ILk21=ak2IHk2ILk21,k1,k2[1,3]

训练过程分为三个步骤:
1.对具有相同模糊级别的地震速度模型和地震数据进行训练,包括四种类型的训练数据(与上文类似)。上式参数 m m m设置为1。此时模型仅针对特定数据进行初步学习,泛化能力较差。
2.对地震速度模型和各种模糊级别的地震数据进行训练,包括16种训练数据。因此,该模型广泛学习各种数据,具有很强的泛化能力。 m m m值设置同上。
3.对地震速度模型和各种模糊级别的地震数据的训练与步骤2中的训练数据类似(图3中的所有部分)。此外,将参数 m m m调整为 m > 1 m > 1 m>1,这使得损失函数对地震速度模型更加敏感。

2.2 实验结果

在这里插入图片描述

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.coloradmin.cn/o/896015.html

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈,一经查实,立即删除!

相关文章

计算机丢失msvcp110.dll是什么意思?如何修复?

在日常使用电脑的时候,有时候会遇到一些使用问题。 比如,有一次遇到了这样一个问题。 那就是,因为“msvcp110.dll”这个文件丢失,有些软件安装不了。 计算机丢失msvcp110.dll是什么意思?该如何修复呢? ​…

05-微信小程序常用组件-表单组件

05-微信小程序常用组件-表单组件 文章目录 表单组件button 按钮案例代码 form 表单案例代码 image 图片支持长按识别的码案例代码 微信小程序包含了六大组件: 视图容器、 基础内容、 导航、 表单、 互动和 导航。这些组件可以通过WXML和WXSS进行布局和样式设…

LC-二叉树最大深度

LC-二叉树最大深度 链接:https://leetcode.cn/problems/maximum-depth-of-binary-tree/description/ 描述: 给定一个二叉树 root ,返回其最大深度。 二叉树的 最大深度 是指从根节点到最远叶子节点的最长路径上的节点数。 例1: …

一力破万法,Fiddler跟Charles抓包图文详解教程

一、Fiddler 1.1 Fiddler抓取HTTPS设置 1.1.1 配置证书 Tools菜单 —> Options —> HTTPS —> 勾选Decrypt HTTPS traffic选项。 说明: 勾选Decrypt HTTPS traffic选项:Decrypt HTTPS traffic:意思是解密HTTPS流量(…

[VS/C++]如何更好的配置DLL项目中的成品输出

注意,解决方案与项目不放在同一个文件夹中,即不选中图中选项 直入主题 首先右键项目选择属性,或者选中项目然后AltEnter 选择配置属性下的常规 分别在四种配置中编辑输出目录如下 注意,四种配置要分别配置,一个个来…

leetcode做题笔记86分隔链表

给你一个链表的头节点 head 和一个特定值 x ,请你对链表进行分隔,使得所有 小于 x 的节点都出现在 大于或等于 x 的节点之前。 你应当 保留 两个分区中每个节点的初始相对位置。 示例 1: 输入:head [1,4,3,2,5,2], x 3 输出&am…

用pytorch实现google net

GoogleNet(也称为Inception v1)是由Google在2014年提出的一个深度卷积神经网络架构。它在ImageNet Large Scale Visual Recognition Challenge (ILSVRC) 2014比赛中取得了优秀的成绩,并引起了广泛的关注。 GoogleNet的设计目标是构建一个更…

【学会动态规划】最长湍流子数组(23)

目录 动态规划怎么学? 1. 题目解析 2. 算法原理 1. 状态表示 2. 状态转移方程 3. 初始化 4. 填表顺序 5. 返回值 3. 代码编写 写在最后: 动态规划怎么学? 学习一个算法没有捷径,更何况是学习动态规划, 跟我…

LeetCode49.字母异味词分组

我一开始的思路就是用1个hashmap<Integer,List<String>>,Integer存的的是字符串所有字母ASCLL值的和&#xff0c;List里面放异位字符串&#xff0c;但是不是异位的字符串的ascll值也可能相同比如acd和abe&#xff0c;所以这个hashmap只能降低一点时间复杂度我还是要…

整理mongodb文档:索引

个人博客 整理mongodb文档:索引 个人博客&#xff0c;求关注&#xff0c;有问题的地方欢迎指出&#xff0c;觉得讲解的繁琐的也请指出 文章概叙 本文主要还是在shell下的操作。让大家了解下mongodb中index的用途&#xff0c;基本的索引的新增、查看、删除等&#xff0c;最后…

使用Dockker创建vwas容器时报错的解决方法

执行命令 docker run -it -d -p 13443:3443 --cap-add LINUX_IMMUTABLE secfa/docker-awvs没有详细看报错之前找了各种各样的解决办法&#xff0c;都无法解决。因此以后在看报错提示的时候耐心一点看关键词Error 后来才发现启动vwas时docker报了这个错&#xff1a; OSError: …

记忆正则表达式的基本元件

正则常见的三种功能&#xff0c;它们分别是&#xff1a;校验数据的有效性、查找符合要求的文本以及对文本进行切割和替换等操作。 正则表达式&#xff0c;简单地说就是描述字符串的规则。在正则中&#xff0c;普通字符表示的还是原来的意思&#xff0c;比如字符 a&#xff0c;…

【编织时空三:探究顺序表与链表的数据之旅】

本章重点 链表OJ题 1. 删除链表中等于给定值 val 的所有结点。 OJ链接 思路一&#xff1a;删除头结点时另做考虑&#xff08;由于头结点没有前一个结点&#xff09; struct ListNode* removeElements(struct ListNode* head, int val) {assert(head);struct ListNode* cur h…

利用Opencv实现人像迁移

前言&#xff1a; Hello大家好&#xff0c;我是Dream。 今天来学习一下如何使用Opencv实现人像迁移&#xff0c;欢迎大家一起参与探讨交流~ 本文目录&#xff1a; 一、实验要求二、实验环境三、实验原理及操作1.照片准备2.图像增强3.实现美颜功能4.背景虚化5.图像二值化处理6.人…

操作系统-笔记-第二章-进程调度

目录 二、第二章——【进程调度】 1、调度的概念 &#xff08;1&#xff09;五状态和七状态&#xff08;就绪挂起、阻塞挂起&#xff09; &#xff08;2&#xff09;三层调度 &#xff08;高级、中级、低级&#xff09; &#xff08;3&#xff09;总结 2、调度的切换 &a…

Vue--》打造个性化医疗服务的医院预约系统(六)

今天开始使用 vue3 + ts 搭建一个医院预约系统的前台页面,因为文章会将项目的每一个地方代码的书写都会讲解到,所以本项目会分成好几篇文章进行讲解,我会在最后一篇文章中会将项目代码开源到我的GithHub上,大家可以自行去进行下载运行,希望本文章对有帮助的朋友们能多多关…

RGOS日常管理操作

RGOS日常管理操作 一、前言二、RGOS平台概述2.1、锐捷设备的常用登陆方式2.2、使用Console登入2.3、Telnet远程管理2.4、SSH远程管理2.5、登陆软件&#xff1a;SecureCRT 三、CLI命令行操作3.1、CLI命令行基础3.2、CLI模式3.3、CLI模式互换3.4、命令行特性3.4.1、分屏显示3.4.2…

(六)、深度学习框架中的算子

1、深度学习框架算子的基本概念 深度学习框架中的算子&#xff08;operator&#xff09;是指用于执行各种数学运算和操作的函数或类。这些算子通常被用来构建神经网络的各个层和组件&#xff0c;实现数据的传递、转换和计算。 算子是深度学习模型的基本组成单元&#xff0c;它们…

07-微信小程序-注册页面-模块化

07-微信小程序-注册页面 文章目录 注册页面使用 Page 构造器注册页面参数Object初始数据案例代码 生命周期回调函数组件事件处理函数setData()案例代码 生命周期模块化 注册页面 对于小程序中的每个页面&#xff0c;都需要在页面对应的 js 文件中进行注册&#xff0c;指定页面…

Ribbon 源码分析

Ribbon 源码分析 Ribbon Debug 分析 断点 LoadBalancerInterceptor LoadBalancerInterceptor 实现了 ClientHttpRequestInterceptor 接口&#xff0c;重写了其中的 intercept 方法&#xff0c;用来拦截请求&#xff1b; 获取原始的 uri 和 服务名&#xff0c;调用 LoadBalanc…