- 博客(17)
- 收藏
- 关注
RSS订阅原创 模拟各向同性弹性波波场-matlab程序
各向异性波场模拟是一种用于描述波在不同方向上传播特性的技术，广泛应用于地质勘探和地震学。在地质勘探中，各向异性波场模拟帮助确定油气储层的位置和特性，而在地震学中，它用于研究地壳中的应力分布。通过数值方法如有限差分法或伪谱法，模拟解决复杂介质中的波动方程，提高预测精度。在石油勘探中，它提升了资源定位的准确性，而在地震研究中，增强了解地壳结构和动力学的能力。模拟的高效率和处理复杂介质的能力，使其成为理解波传播的重要工具。- 均匀介质模型：vp=2000m/s, vs=1200m/s, ρ=2.60g/cm³。
2025-04-21 14:36:09
260
原创 地质三维建模案例
param layers: 地层参数列表，每个元素为(depth, permeability, porosity):param resolution: 网格分辨率(nx,ny,nz)# 2. 定义地层参数 (深度(m), 渗透率(m/s), 孔隙度):param size: 模型尺寸(x,y,z) 单位：米。(300, 1e-6, 0.15), # 黏土层。(100, 1e-4, 0.25), # 表层土。(200, 5e-5, 0.20), # 砂层。"""求解地下水流动方程"""
2025-04-20 23:18:19
479
原创 39-多道面波频散曲线自动提取
多道面波分析(MASW)是一种基于瑞利波频散特性的浅层地震勘探方法，其核心是通过提取频散曲线来反演地下横波速度结构。本次实现了从二维频散图像中提取基阶模态频散曲线及其上下边界的关键功能。
2025-04-20 23:12:43
454
原创 11-FWI全波形反演理论及程序剖析、运行效果展示
这段MATLAB代码实现了地震波正演模拟及数据处理流程，采用并行计算加速处理。将从功能、结构和技术细节三个方面进行详细分析。
2025-04-20 23:06:56
217
原创 9-MATLAB代码实现地震逆时偏移(RTM)成像算法
现代逆时偏移技术通常采用高阶有限差分法求解波动方程，配合PML吸收边界条件，并借助GPU并行计算提升效率。随着全波形反演（FWI）技术的发展，逆时偏移不仅作为独立的成像工具，更成为构建高精度速度模型的关键环节。尽管计算成本较高，但在深海勘探、非常规油气开发等前沿领域，逆时偏移已成为解决复杂地质问题的首选成像手段。这段MATLAB代码实现了地震逆时偏移(RTM)成像算法，用于地下速度结构成像。将从功能、结构和技术细节三个方面进行详细剖析。
2025-04-20 22:59:25
511
原创 7562i-对信号进行连续小波变换(CWT)分析
这段代码提供了完整的CWT分析流程，适合机械振动、生物信号等非平稳信号分析。通过调整参数可以适应不同应用场景的需求。
2025-04-20 22:51:57
422
原创 59276-交错网格有限差分法瑞利波模拟
本MATLAB代码实现了基于交错网格有限差分法的瑞利波场数值模拟，适用于地震波传播模拟和浅地表勘探研究。
2025-04-20 11:58:39
705
原创 4670-快速傅里叶变换时域信号图像
本代码实现完整的信号时频分析流程，包含数据读取、预处理、时域分析、FFT变换及可视化功能模块。适用于各类振动信号、声学信号的频谱特征分析。
2025-04-20 11:51:48
158
原创 31029-biot-pml边界条件波场模拟-双相介质波场模拟-带孔隙度参数-空间差分精度参数定义
本MATLAB代码实现了基于Biot理论的双相孔隙介质中弹性波传播的高阶有限差分数值模拟。代码采用速度-应力交错网格方法，通过分离固体骨架和孔隙流体的运动方程，模拟了快P波、慢P波和S波在双相介质中的传播特性。
2025-04-20 11:36:14
530
原创 20357-基于拉东变换的纵横波波场分离技术研究
本研究提出了一种基于拉东变换的弹性波场纵横波分离新方法。通过构建波数-频率域的拉东变换算子，实现了P波和S波的高精度分离。数值模拟结果表明，该方法在复杂介质条件下仍能保持90%以上的波场分离精度，且计算效率较传统方法提升约30%。理论分析和实际应用证实了拉东变换在解决弹性波场混叠问题上的独特优势。
2025-04-20 11:25:52
854
原创 103275-二维交错网格纵横波分离的弹性波模拟及应用
本研究提出了一种基于交错网格的二维弹性波数值模拟方案，着重阐述了纵横波分离算法的实现过程。采用高阶有限差分格式结合交错网格离散化技术，实现了弹性波场的高精度数值模拟，并通过波场分解技术有效分离出独立的P波和S波成分。数值实验证实了方法的可靠性和实用性，同时探讨了该方法在复杂地质建模和地震波场解析中的潜在应用价值。
2025-04-20 11:07:20
464
原创 seismic unix 制作数据 合成数据
这一部分是su的基础，在装配完su后，建议先熟悉环境，把linux基础操作，vim基础操作熟悉后，首先进行su脚本，语法规则等的基础学习。黑色的文件名，文件，在命令行下，more/vim+文件名可以查看内容。蓝色的文件名，文件夹，在命令行下cd+名字可以进到该文件夹内。绿色的文件名，可执行文件，在命令行下./+名字可以直接运行。ls -a 是列出当前文件夹下所有文件，包括隐藏文件/夹。ls是列出当前文件夹下文件，不包括隐藏文件/夹。cd .. 注意是两个点，返回上一级目录。
2023-10-10 13:00:13
265
原创 面波频散分析
利用地震数据进行近地表地质结构地震勘探，最为常见的方法有反射波法和 折射波法。近几年来，面波勘探法得到广泛应用，并 且面波的能量强，易于激发；同时，瑞雷 波携带大量的地下介质信息，所以通过面波的频散特性可以研究地层的结构和性质。通过分析面波的频散现象和地表结构之间的关系，使用标量传递矩阵进行频散曲线正演，确定敏感参数；2015年接触面波勘探，在**实习一年，做了面波频散分析-面波频散分析结果提取-面波频散曲线应用（反演），对面波的频散效果有了进一步认识。但是，为了毕业，没有把全部时间放在面波相关的创作上。
2023-10-10 12:37:53
1003
1
原创 seismic unix安装教程 linux红帽
2、修改环境变量，先在家目录下（终端打开就是，如果不在用cd命令后面不加参数就可以进去,用pwd就可以看路径，比如我的路径是 /home/user ，）用ls -a 查看文件。下载好安装包，目前版本是：cwp_su_all_40.gz，用tar -xvzf cwp_su_all_40.gz 解压出src文件夹,可以先在家目录下创建su文件夹，把src放在su内。末尾处有 /home/user/su/bin 就万事大吉了，关掉终端，重新打开一个，确保万无一失，如果不正确，则重新设置。
2023-10-10 12:22:21
654
原创 初识手写数字识别碰到的问题
if step % 5==0://这里每隔5个step打印报错，分析原因是因为step不是数值，不能和整数运算。db_test=db_test.map(preprocess).batch(batchsize)//不停的迭代。
2023-04-11 00:42:09
157
1
原创 损失函数的计算-LOSS（MSE、交叉熵）
Cross Entropy Loss（交叉熵）-专门应对分类问题，衡量不确定性。差的平方的累加，再平均。learningrate对数值比较大的loss起到调节作用。
2023-04-10 10:34:16
384
原创 tf.random.normal()和tf.random.uniform()
tf.random_normal()函数用于从“服从指定正态分布的序列”中随机取出指定个数的值。
2023-04-10 09:54:12
459
2
空空如也
空空如也
TA创建的收藏夹 TA关注的收藏夹
TA关注的人