![](https://img-blog.csdnimg.cn/20201014180756928.png?x-oss-process=image/resize,m_fixed,h_64,w_64)
地球物理
张建明0704
地球物理学在读,研究兴趣主要但不限于地球物理正反演算法,包括波场模拟,走时计算,层析成像,全波形反演,偏移成像,弹性和各向异性介质理论方法。
展开
-
各向同性介质二维声波波动方程有限差分法数值模拟(波场模拟动图)
弹性波这是二阶位移方程的波场模拟,没有交错网格。一阶速度应力-应力方程及其交错网格波场模拟见我前面的帖子。C/C++,MATLAB代码的帖子都可以找到。原创 2024-04-09 20:45:45 · 483 阅读 · 0 评论 -
各向同性与TTI介质各向异性走时场,射线路径,波场的差异
对称模型,走时场,射线路径,波场等波现象均有良好的对称性,容易现象。对称模型,但走时场,射线路径等波现象均无对称性,比各向同性复杂很多。原创 2024-04-08 19:09:40 · 264 阅读 · 0 评论 -
各向异性介质多参数化走时反演敏感核(辐射模式)分析
在阅读文章的基础上,很多朋友对此研究很感兴趣,想进一步深入学习,在此,提供了敏感核分析的代码,下面是Thomsen参数化的代码,其他参数化的代码一并打包上传。(好吧,发现不能打包上传,谁有办法,评论区告诉我!(需要代码的私信我就好))说明:做该项研究当时用好几种办法推导过 结果都很相似 验证了我的推导是对的。原创 2024-04-02 10:15:30 · 355 阅读 · 1 评论 -
二维和三维平滑滤波(Gauss平滑、SU Smooth、 Dave Hale‘s method)
Gauss平滑计算量大,在大模型平滑时非常耗时,在此基础上,提供另外两种高效的平滑程序,分别来自于SU软件和Dave Hale的递归滤波平滑的修改版本。原创 2024-03-29 21:26:27 · 276 阅读 · 1 评论 -
快速傅里叶变换(FFT),时域<-->频域信号转换,以及频率域低通滤波、带通滤波、高通滤波
以下面的地震数据为例,原始地震数据包含低频和高频的噪音,经过4-40Hz带通滤波后,噪音明显被压制,有效信号得到加强。原创 2024-03-29 21:30:30 · 190 阅读 · 1 评论 -
根据克里斯托弗方程(Christoffel equation)计算相速度和偏振(极化)方向,进而计算群速度
弹性波方程是通过应力(stress)应变(strain)关系(广义胡克定律,或本构关系)+运动微分方程(动量守恒线性方程 Linearized equation of momentum conservation Navier 方程或牛顿第二定律)导出,描述了波的所有动力学和运动学传播特性。而Christoffel 方程是各向异性弹性波动方程的平面波(单相波,对应相速度)解形式,可以用来独立描述三种类型的波的极化和相速度特征,即运动学(速度)和偏振(振幅以及方向)。3.几大方程的物理意义是什么?原创 2024-03-12 20:15:05 · 766 阅读 · 1 评论 -
三维各向同性快速扫描法3DFSM
之前的快速扫描法的帖子中提供的二维代码稍显凌乱,为此追加提供一个新的三维的代码,更加清楚直观地展示三维FSM的过程,便于初学者学习int ny=101;i原创 2024-02-27 09:47:25 · 1141 阅读 · 1 评论 -
地震波形反演实用化道路:从理论全波形反演FWI ->分阶段,分尺度, 多信息联合反演
地震波形反演实用化道路,从理论全波形反演FWI -> 多信息,分阶段,分尺度联合反演原创 2023-11-23 15:58:59 · 454 阅读 · 1 评论 -
各向异性TTI介质快速扫描法(Fast Sweeping Method,FSM)走时计算
本人做了复现,虽然是复现,但是也花了很久的时间。其中的一些问题,想了很久才解决。当然实现了以后,就觉得所有的问题原来都不是什么特别的难题。2.沿x轴和z轴群速度的计算,可以采用Shooting方法(Qian,2012),或者Fomel(2004)近似公式。1.因果律的物理含义,代表群速度(射线方向)要和走时更新的方向一致,其实代表了射线的传播方向,即群速度方向。上面这俩式子很重要,又能用来推导沿x轴和z轴群速度,还能用来做因果律的判定。直接上结果了,代码有需要的私聊我。原创 2023-11-09 20:57:54 · 296 阅读 · 1 评论 -
程函方程伴随状态法走时层析,初至及反射,VTI介质多参数
实际资料例子:原创 2023-10-18 21:41:01 · 277 阅读 · 1 评论 -
散射积分法高斯-牛顿Hessian优化的射线层析
关于散射积分法的一些实验。原创 2023-10-18 21:34:23 · 94 阅读 · 1 评论 -
各向异性TTI介质地震波场模拟,速度-应力弹性波方程交错网格高阶有限差分法,PML边界
逐渐改进,从各向同性到VTI各向异性,再到TTI介质。模拟效果可以在走时计算的文章里看到,有对比。不懂的欢迎随时交流。原创 2023-10-18 21:14:10 · 676 阅读 · 0 评论 -
各向异性VTI介质地震波场模拟,速度-应力弹性波方程交错网格有限差分法,PML边界(C++和MATLAB两种语言的代码都有)
参考各向同性的程序,进一步修改为VTI介质,不懂的随时交流。原创 2023-10-18 21:11:46 · 369 阅读 · 4 评论 -
各向同性介质地震波场模拟,速度-应力弹性波方程交错网格有限差分法
源码如下,不懂的随时交流。原创 2023-10-18 21:09:10 · 448 阅读 · 10 评论 -
各向异性VTI介质快速扫描算法FSM走时计算及其改进算法
【代码】各向异性VTI介质快速扫描算法FSM走时计算及其改进算法。原创 2023-10-18 20:38:16 · 141 阅读 · 0 评论 -
地震初至波和反射波最短路径走时计算及射线追踪(射线路径回追)
涉及暂未发表的成果,算法思路和代码仅限私下交流,欢迎评论交流提问。原创 2023-10-18 20:33:57 · 228 阅读 · 5 评论 -
快速扫描法(Fast Sweeping Method,FSM)及最短路径辅助的快速扫描法SPAFSM走时计算
代码包含快速扫描法(FSM)走时计算的基本方法(Zhao 2005)以及高阶FSM方法和对FSM方法的一点改进:添加最短路径辅助解提高计算精度(SPAFSM,Zhang 2023),仅供初学者参考,如有疑问,欢迎随时交流讨论。//i是次数,j是具体下标。原创 2023-10-18 19:57:54 · 550 阅读 · 0 评论 -
各向异性TTI介质初至走时、反射走时计算及其射线路径追踪
涉及为未发表成果,代码暂时不公开,欢迎私下交流讨论,互相学习。原创 2023-10-18 20:47:40 · 148 阅读 · 0 评论