- 博客(22)
- 收藏
- 关注
原创 已知走时场和炮点检波点的位置做射线追踪的方法,注意各向同性和各向异性介质有所不同
而在各向异性介质中进行射线追踪时,需要沿着群速度方向,而不是简单地沿走时场的梯度方向。我们需要计算波的群速度,并沿群速度方向进行追踪。在各向同性介质中进行射线追踪时,可以简单地沿走时场的梯度方向既可以追踪到射线路径。
2024-07-25 21:56:25 124
原创 基于程函方程(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
原创 各向异性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
原创 快速扫描法(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关注的人