- 博客(25)
- 收藏
- 关注
RSS订阅原创 各种环境下安装软件的命令对比与总结
环境安装命令特点Linux APTDebian/Ubuntu 系统专用Linux YUMCentOS/RHEL 系统专用Python pipPython 包专用conda跨语言，适合科学计算环境npmNode.js 包管理源码编译最通用但最复杂跨发行版沙盒应用选择哪种安装方式取决于：你的操作系统类型要安装的软件类型是否需要特定版本是否需要环境隔离通常建议优先使用系统包管理器，其次是语言专用包管理器，最后考虑通用安装方式。
2025-03-31 20:40:18
928
原创 全球地震数据和各个国家的地震数据网络介绍
这些网络和机构提供的地震数据可以通过其官方网站访问，支持科学研究和工程应用。根据您的需求，可以选择合适的网络和机构获取数据。
2024-12-25 19:06:56
3001
原创 已知走时场和炮点检波点的位置做射线追踪的方法，注意各向同性和各向异性介质有所不同
而在各向异性介质中进行射线追踪时，需要沿着群速度方向，而不是简单地沿走时场的梯度方向。我们需要计算波的群速度，并沿群速度方向进行追踪。在各向同性介质中进行射线追踪时，可以简单地沿走时场的梯度方向既可以追踪到射线路径。
2024-07-25 21:56:25
599
1
原创 基于程函方程（Eikonal Equation）走时计算的偏移成像
实际上就是绕射扫描叠加原理，或者说是用到Kirchhoff偏移的简化形式。该成像方式速度快，内存占用低。也可以用来做成像道集监控或者偏移速度分析。
2024-07-22 15:30:02
911
原创 龙格-库塔（Runge-kutta）法解Hamilton方程组（Eikonal方程）进行射线追踪（MATLAB代码）
该代码完全由ChatGPT生成，对于快速学习入门真的又好又快！除了各向同性的还有各向异性的以及三维的代码。（心得：万事不会求助ChatGPT！
2024-07-19 21:28:07
900
原创 各向同性介质二维声波波动方程有限差分法数值模拟（波场模拟动图）
弹性波这是二阶位移方程的波场模拟，没有交错网格。一阶速度应力-应力方程及其交错网格波场模拟见我前面的帖子。C/C++，MATLAB代码的帖子都可以找到。
2024-04-09 20:45:45
1918
原创 各向同性与TTI介质各向异性走时场，射线路径，波场的差异
摘要：研究对比了各向同性介质与TTI各向异性介质（Eps=0.4,Del=0.2，倾角45°）的波传播特性。各向同性介质具有对称的走时场和射线路径，而TTI各向异性介质则表现出明显的不对称性和复杂性。通过起伏界面模型和BP模型的对比发现，在复杂介质中，是否考虑各向异性（特别是倾角）对走时、射线路径及波前匹配均有显著影响。VTI与TTI的处理方式差异会导致不同的模拟结果，凸显了各向异性因素在地震波研究中的重要性。
2024-04-08 19:09:40
703
1
原创 各向异性介质多参数化走时反演敏感核（辐射模式）分析
在阅读文章的基础上，很多朋友对此研究很感兴趣，想进一步深入学习，在此，提供了敏感核分析的代码，下面是Thomsen参数化的代码，其他参数化的代码一并打包上传。（好吧，发现不能打包上传，谁有办法，评论区告诉我！（需要代码的私信我就好））说明：做该项研究当时用好几种办法推导过 结果都很相似 验证了我的推导是对的。
2024-04-02 10:15:30
594
1
原创 快速傅里叶变换（FFT），时域＜--＞频域信号转换，以及频率域低通滤波、带通滤波、高通滤波
以下面的地震数据为例，原始地震数据包含低频和高频的噪音，经过4-40Hz带通滤波后，噪音明显被压制，有效信号得到加强。
2024-03-29 21:30:30
652
1
原创 二维和三维平滑滤波（Gauss平滑、SU Smooth、 Dave Hale‘s method）
Gauss平滑计算量大，在大模型平滑时非常耗时，在此基础上，提供另外两种高效的平滑程序，分别来自于SU软件和Dave Hale的递归滤波平滑的修改版本。
2024-03-29 21:26:27
661
1
原创 根据克里斯托弗方程（Christoffel equation）计算相速度和偏振(极化)方向，进而计算群速度
弹性波方程是通过应力（stress）应变（strain）关系（广义胡克定律，或本构关系）+运动微分方程（动量守恒线性方程 Linearized equation of momentum conservation Navier 方程或牛顿第二定律）导出，描述了波的所有动力学和运动学传播特性。而Christoffel 方程是各向异性弹性波动方程的平面波（单相波，对应相速度）解形式，可以用来独立描述三种类型的波的极化和相速度特征，即运动学（速度）和偏振（振幅以及方向）。3.几大方程的物理意义是什么？
2024-03-12 20:15:05
2600
1
原创 三维各向同性快速扫描法3DFSM（附C++和MATLAB代码）
之前的快速扫描法的帖子中提供的二维代码稍显凌乱，为此追加提供一个新的三维的代码，更加清楚直观地展示三维FSM的过程，便于初学者学习int ny=101;i
2024-02-27 09:47:25
1657
1
原创 地震波形反演实用化道路：从理论全波形反演FWI -＞分阶段，分尺度， 多信息联合反演
本文提出了一种基于波动方程的初至波多信息联合反演方法（JTEW），用于提高近地表速度建模精度。该方法创新性地在统一波动理论框架下，分阶段利用初至波的走时、包络和波形信息进行联合反演。相比传统方法，JTEW具有三大优势：（1）只需一次波动方程正演即可同时计算三种信息，无需额外射线追踪；（2）突破射线理论高频假设限制，提升复杂地质条件下的适用性；（3）通过权重因子动态调节不同信息比例，有效缓解目标函数漂移问题，降低对初始模型的依赖。数值实验和实际资料处理表明，该方法能显著提高反演精度和效率，为初至波波形反演的实
2023-11-23 15:58:59
903
1
原创 二维和三维TTI各向异性介质快速扫描法（Fast Sweeping Method，FSM）走时计算，以及因式分解程函方程解法（Factored Eikonal Equation）
本人做了复现，虽然是复现，但是也花了很久的时间。其中的一些问题，想了很久才解决。当然实现了以后，就觉得所有的问题原来都不是什么特别的难题。2.沿x轴和z轴群速度的计算，可以采用Shooting方法（Qian，2012），或者Fomel（2004）近似公式。1.因果律的物理含义，代表群速度（射线方向）要和走时更新的方向一致，其实代表了射线的传播方向，即群速度方向。上面这俩式子很重要，又能用来推导沿x轴和z轴群速度，还能用来做因果律的判定。直接上结果了，代码有需要的私聊我。
2023-11-09 20:57:54
1134
5
原创 程函方程伴随状态法走时层析，初至及反射，VTI介质多参数、双差走时层析，及其在勘探与固体地球物理中的应用
本文系统梳理了伴随状态法走时层析技术的发展历程与研究进展。该方法由Sei和Symes（1994）首创，经国内外学者持续改进：Leung和Qian（2006）、Taillandier（2009）完善理论框架；我国学者谢春、刘玉柱、董良国等（2014起）引入并发展出预条件伴随状态法、VTI各向异性反演及反射波走时反演技术。针对反射层析的关键问题，作者提出ASCRT方法解决界面定位难题（2023-2024），并拓展至多分量地震数据反演（2025）。创新性工作包括：程函方程特征反射层析、伴随双差走时层析（ASDDT
2023-10-18 21:41:01
1526
13
原创 各向异性TTI介质地震波场模拟，速度-应力弹性波方程交错网格高阶有限差分法，PML吸收边界（C++和MATLAB两种语言的代码都有）
逐渐改进，从各向同性到VTI各向异性，再到TTI介质。模拟效果可以在走时计算的文章里看到，有对比。不懂的欢迎随时交流。
2023-10-18 21:14:10
2145
1
原创 各向异性VTI介质地震波场模拟，速度-应力弹性波方程交错网格有限差分法，PML吸收边界（C++和MATLAB两种语言的代码都有）
参考各向同性的程序，进一步修改为VTI介质，不懂的随时交流。
2023-10-18 21:11:46
1111
6
原创 地震初至波和反射波最短路径走时计算及射线追踪（射线路径回追）
摘要：本文提出了一种融合最短路径法(SPM)和快速扫描法(FSM)优势的新方法。通过交替扫描顺序替代传统SPM中的堆排序，显著提高了计算效率；通过设置次级网格增加射线传播方向，提升了计算精度。该方法采用四种扫描顺序获取全局最小走时，并基于费马原理实现射线路径回溯。数值实验表明，该方法在计算效率和精度之间实现了良好平衡，且更易于实现。
2023-10-18 20:33:57
810
5
原创 快速扫描法（Fast Sweeping Method，FSM）及最短路径辅助的快速扫描法SPAFSM走时计算
代码包含快速扫描法（FSM）走时计算的基本方法（Zhao 2005）以及高阶FSM方法和对FSM方法的一点改进：添加最短路径辅助解提高计算精度（SPAFSM，Zhang 2023），仅供初学者参考，如有疑问，欢迎随时交流讨论。//i是次数，j是具体下标。
2023-10-18 19:57:54
1488
3
空空如也
空空如也
TA创建的收藏夹 TA关注的收藏夹
TA关注的人