自定义博客皮肤VIP专享

*博客头图:

格式为PNG、JPG,宽度*高度大于1920*100像素,不超过2MB,主视觉建议放在右侧,请参照线上博客头图

请上传大于1920*100像素的图片!

博客底图:

图片格式为PNG、JPG,不超过1MB,可上下左右平铺至整个背景

栏目图:

图片格式为PNG、JPG,图片宽度*高度为300*38像素,不超过0.5MB

主标题颜色:

RGB颜色,例如:#AFAFAF

Hover:

RGB颜色,例如:#AFAFAF

副标题颜色:

RGB颜色,例如:#AFAFAF

自定义博客皮肤

-+
  • 博客(22)
  • 收藏
  • 关注

原创 已知走时场和炮点检波点的位置做射线追踪的方法,注意各向同性和各向异性介质有所不同

而在各向异性介质中进行射线追踪时,需要沿着群速度方向,而不是简单地沿走时场的梯度方向。我们需要计算波的群速度,并沿群速度方向进行追踪。在各向同性介质中进行射线追踪时,可以简单地沿走时场的梯度方向既可以追踪到射线路径。

2024-07-25 21:56:25 124

原创 快速行进法(FMM)求解程函方程的基本原理,以及VTI介质和TTI介质中的走时计算

【代码】快速行进法(FMM)求解程函方程的基本原理。

2024-07-25 20:28:44 67

原创 基于程函方程(EIkonal Equation)走时计算的偏移成像

实际上就是绕射扫描叠加原理,或者说是用到Kirchhoff偏移的简化形式。该成像方式速度快,内存占用低。也可以用来做成像道集监控或者偏移速度分析。

2024-07-22 15:30:02 148

原创 龙格-库塔(Runge-kutta)法解Hamilton方程组(EIkonal方程)进行射线追踪(MATLAB代码)

该代码完全由ChatGPT生成,对于快速学习入门真的又好又快!除了各向同性的还有各向异性的以及三维的代码。(心得:万事不会求助ChatGPT!

2024-07-19 21:28:07 204

原创 各向同性介质二维声波波动方程有限差分法数值模拟(波场模拟动图)

弹性波这是二阶位移方程的波场模拟,没有交错网格。一阶速度应力-应力方程及其交错网格波场模拟见我前面的帖子。C/C++,MATLAB代码的帖子都可以找到。

2024-04-09 20:45:45 515

原创 各向同性与TTI介质各向异性走时场,射线路径,波场的差异

对称模型,走时场,射线路径,波场等波现象均有良好的对称性,容易现象。对称模型,但走时场,射线路径等波现象均无对称性,比各向同性复杂很多。

2024-04-08 19:09:40 271

原创 各向异性介质多参数化走时反演敏感核(辐射模式)分析

在阅读文章的基础上,很多朋友对此研究很感兴趣,想进一步深入学习,在此,提供了敏感核分析的代码,下面是Thomsen参数化的代码,其他参数化的代码一并打包上传。(好吧,发现不能打包上传,谁有办法,评论区告诉我!(需要代码的私信我就好))说明:做该项研究当时用好几种办法推导过 结果都很相似 验证了我的推导是对的。

2024-04-02 10:15:30 361 1

原创 快速傅里叶变换(FFT),时域<-->频域信号转换,以及频率域低通滤波、带通滤波、高通滤波

以下面的地震数据为例,原始地震数据包含低频和高频的噪音,经过4-40Hz带通滤波后,噪音明显被压制,有效信号得到加强。

2024-03-29 21:30:30 199 1

原创 二维和三维平滑滤波(Gauss平滑、SU Smooth、 Dave Hale‘s method)

Gauss平滑计算量大,在大模型平滑时非常耗时,在此基础上,提供另外两种高效的平滑程序,分别来自于SU软件和Dave Hale的递归滤波平滑的修改版本。

2024-03-29 21:26:27 286 1

原创 根据克里斯托弗方程(Christoffel equation)计算相速度和偏振(极化)方向,进而计算群速度

弹性波方程是通过应力(stress)应变(strain)关系(广义胡克定律,或本构关系)+运动微分方程(动量守恒线性方程 Linearized equation of momentum conservation Navier 方程或牛顿第二定律)导出,描述了波的所有动力学和运动学传播特性。而Christoffel 方程是各向异性弹性波动方程的平面波(单相波,对应相速度)解形式,可以用来独立描述三种类型的波的极化和相速度特征,即运动学(速度)和偏振(振幅以及方向)。3.几大方程的物理意义是什么?

2024-03-12 20:15:05 803 1

原创 三维各向同性快速扫描法3DFSM(附C++和MATLAB代码)

之前的快速扫描法的帖子中提供的二维代码稍显凌乱,为此追加提供一个新的三维的代码,更加清楚直观地展示三维FSM的过程,便于初学者学习int ny=101;i

2024-02-27 09:47:25 1148 1

原创 地震波形反演实用化道路:从理论全波形反演FWI ->分阶段,分尺度, 多信息联合反演

地震波形反演实用化道路,从理论全波形反演FWI -> 多信息,分阶段,分尺度联合反演

2023-11-23 15:58:59 457 1

原创 各向异性TTI介质快速扫描法(Fast Sweeping Method,FSM)走时计算

本人做了复现,虽然是复现,但是也花了很久的时间。其中的一些问题,想了很久才解决。当然实现了以后,就觉得所有的问题原来都不是什么特别的难题。2.沿x轴和z轴群速度的计算,可以采用Shooting方法(Qian,2012),或者Fomel(2004)近似公式。1.因果律的物理含义,代表群速度(射线方向)要和走时更新的方向一致,其实代表了射线的传播方向,即群速度方向。上面这俩式子很重要,又能用来推导沿x轴和z轴群速度,还能用来做因果律的判定。直接上结果了,代码有需要的私聊我。

2023-11-09 20:57:54 304 1

原创 程函方程伴随状态法走时层析,初至及反射,VTI介质多参数

实际资料例子:

2023-10-18 21:41:01 285 1

原创 散射积分法高斯-牛顿Hessian优化的射线层析

关于散射积分法的一些实验。

2023-10-18 21:34:23 97 1

原创 各向异性TTI介质地震波场模拟,速度-应力弹性波方程交错网格高阶有限差分法,PML吸收边界(C++和MATLAB两种语言的代码都有)

逐渐改进,从各向同性到VTI各向异性,再到TTI介质。模拟效果可以在走时计算的文章里看到,有对比。不懂的欢迎随时交流。

2023-10-18 21:14:10 699

原创 各向异性VTI介质地震波场模拟,速度-应力弹性波方程交错网格有限差分法,PML吸收边界(C++和MATLAB两种语言的代码都有)

参考各向同性的程序,进一步修改为VTI介质,不懂的随时交流。

2023-10-18 21:11:46 384 4

原创 各向同性介质地震波场模拟,速度-应力弹性波方程交错网格有限差分法

源码如下,不懂的随时交流。

2023-10-18 21:09:10 456 10

原创 各向异性TTI介质初至走时、反射走时计算及其射线路径追踪

涉及为未发表成果,代码暂时不公开,欢迎私下交流讨论,互相学习。

2023-10-18 20:47:40 149

原创 各向异性VTI介质快速扫描算法FSM走时计算及其改进算法

【代码】各向异性VTI介质快速扫描算法FSM走时计算及其改进算法。

2023-10-18 20:38:16 146

原创 地震初至波和反射波最短路径走时计算及射线追踪(射线路径回追)

涉及暂未发表的成果,算法思路和代码仅限私下交流,欢迎评论交流提问。

2023-10-18 20:33:57 232 5

原创 快速扫描法(Fast Sweeping Method,FSM)及最短路径辅助的快速扫描法SPAFSM走时计算

代码包含快速扫描法(FSM)走时计算的基本方法(Zhao 2005)以及高阶FSM方法和对FSM方法的一点改进:添加最短路径辅助解提高计算精度(SPAFSM,Zhang 2023),仅供初学者参考,如有疑问,欢迎随时交流讨论。//i是次数,j是具体下标。

2023-10-18 19:57:54 560

空空如也

空空如也

TA创建的收藏夹 TA关注的收藏夹

TA关注的人

提示
确定要删除当前文章?
取消 删除