- 博客(67)
- 资源 (4)
- 问答 (4)
- 收藏
- 关注
RSS订阅转载 学习笔记：波动方程逆时偏移算法
分享下近期学习的逆时偏移相关算法，以下内容来源于，网站所提供的代码很全面（可直接下载），可以修改代码用于地震勘探数据处理。以下是计算结果。
2024-01-17 11:31:42
421
原创 静力触探数据智能预处理（4）
机器学习中有两个重要的任务，即分类和聚类。分类算法有逻辑回归、支持向量机、决策树、随机森林等。聚类算法有K均值聚类、层次聚类、密度聚类、谱聚类、EM聚类、模糊聚类等。K均值聚类是一种常用的无监督学习算法，用于将数据点划分为不同的簇，使得同一簇内的数据点彼此相似度较高，不同簇之间的数据点相似度较低。这个算法的目标是将数据点分为K个簇，其中K是用户定义的参数。K均值聚类的原理相对简单，主要思想是通过迭代寻找K个簇的中心点，将每个数据点分配给距离其最近的中心点，然后更新中心点的位置，直到满足停止条件为止。
2023-10-28 13:03:55
407
原创 不同组合地下管线的地质雷达响应特征分析
以混凝土管线为例，建立水平相邻管线电性模型，管径为60cm，中心埋深为70cm，管线长度为150cm，分别建立管线圆心相距150cm的两根相邻双管线和三管线模型，进行二维地质雷达正演模拟，研究地质雷达探测水平相邻管线电磁波的响应特征。
2023-09-28 17:13:10
282
2
原创 三维跨孔电磁波CT数据可视化框架搭建
利用matlab实现对跨孔电磁波CT实测数据反演，并搭建了三维CT数据可视化框架，可装填实测CT反演数据。
2023-09-04 18:08:52
632
5
原创 静力触探数据智能预处理（2）
数据处理方式已由手工1.0、计算机辅助2.0到人工智能3.0的趋势发展。现场采集的静力触探数据通常是由仪器厂家开发的数据采集软件保存，将原始数据导入Excel中，数据格式需要花费一定的时间整理，当数据量比较大的时候容易出错。利用matlab开发了静力触探数据快速读取程序，可以实现对原始数据的一键提取并快速分类，非专业编写，代码仅供参考。
2023-08-30 22:00:25
709
原创 层状介质一维大地电磁数值模拟
大地电磁测深法（ MT）是根据电磁感应原理研究天然场源在地下激励产生的交变电场或者磁场，通过地表观测到的电磁场或者通过电磁场计算出视电阻率或者相位等量来进行地下构造研究的一种电磁方法。由于它不需要建立人工场源，施工比较方便，另外它对良导体的分辨率能力强，受高阻地层的屏蔽影响较小，探测深度少则几十米，深可达几百公里，广泛应用于海上勘探、工程勘探、油气田勘探、地热资源调查、深部地质构造等地球物理领域。本文应用网上找到的开源MT代码，进行一维大地电磁数值模拟研究。
2023-08-08 17:11:47
541
2
原创 静力触探数据智能预处理（1）
目前，数据处理方式由手工1.0、计算机辅助2.0到人工智能3.0的趋势发展。现场采集的静力触探数据存在大量的异常数据，需要后期处理；但是目前还未见一个静力触探数据预处理的软件，数据预处理主要还是依靠在excel中人眼识别+手工剔除的方式。本博文尝试性的编写了静力触探数据预处理代码，非专业编写，代码仅供参考。
2023-06-28 21:15:02
593
原创 代码分享：gprMax钻孔地质雷达波场模拟
gprMax模拟地面地质雷达被广泛使用，但是在钻孔内进行地质雷达的模拟较少。本博文尝试利用gprMax进行钻孔地质雷达的模拟，代码仅供大家借鉴。为了展示电磁波在钻孔内传播的状态，利用paraview将模拟的结果动态展示。假设钻孔间距为20m，钻孔深度为30m，背景介质设置为花岗岩，设置两个方形充水空洞。模型图如下：雷达天线频率设置为32MHz，gprMax模拟的具体参数设置如下：本参数设置的不是很合理，但是不影响模拟出结果。运行模拟的主程序代码由python编辑，代码如下：注：程序运行出现bug，
2023-03-07 14:10:53
1716
1
原创 代码分享：面波数据快速成图
目前，物探数据主要用sufer软件成图，sufer软件具有强大的插值和绘图功能，成图比较美观。但是，生产过程中大量的物探数据，依靠excel和sufer来成图耗费人力时间成本。本博文在MATLAB平台上开发了一套用于面波数据快速成图的小程序，仅供大家借鉴。
2023-03-06 11:23:33
1285
原创 MATLAB代码实现钻孔地震波波场模拟
钻孔间地震波波场模拟对研究地震波在钻孔内传播非常重要，基于MATLAB实现了一种简单的钻孔地震波波场模拟，适合新手学习使用，对钻孔地震波的传播有个初步的认识。
2022-12-23 19:35:50
1053
8
原创 Pyqt学习笔记01
上面，我们通过一个简单的界面，认识了pyqt的基本结构，当然这非常基础，为了实现复杂的软件界面设计，还有许多内容需要学习。然后是创建一个空白的窗口，执行程序时窗口弹出来，显示在电脑桌面上，当关闭窗口时，窗口退出来。然后就可以在窗口上添加一些控件，直接在后面增加内容就可以了，比如我增加一个标签和一个按钮。一个基本的pyqt包含包和模块的导入、创建一个主界面、控件的操作、界面执行、退出。设计一个pyqt程序，首先需要导入对应的包和模块，为了方便初学，可以都导入进来。
2022-12-23 18:46:03
168
原创 一个简单的MATLAB脚本——快速行进算法（FMM））
介绍快速行进算法（FMM）的简单MATLAB脚本，不到20行代码实现快速行进算法的运算结果，而且计算速度非常快。给了两个实例模型来说明计算结果。
2022-12-14 15:04:12
1474
14
原创 交叉梯度函数的MATLAB实现及代码分享02
交叉梯度函数可用于反演成像中。作为一个连接不同物性参数的桥梁，交叉梯度函数可以实现不同物性参数的联合反演成像。本文给出了三维交叉梯度函数的定义及性质、计算公式、MATLAB代码，并用模型算例给出了计算结果
2022-12-14 10:56:19
955
原创 交叉梯度函数的MATLAB实现及代码分享01
交叉梯度函数可用于反演成像中。作为一个连接不同物性参数的桥梁，交叉梯度函数可以实现不同物性参数的联合反演成像。
2022-12-14 10:05:33
1558
1
原创 基于ttcrpy的跨孔CT高斯牛顿算法及python代码分享（2）
ttcrpy是加拿大学者伯纳德·吉鲁（Bernard Giroux）于2021年发布的开源python库，详见（https://github.com/groupeLIAMG），参考文献（Giroux B. 2021. ttcrpy: A Python package for traveltime computation and raytracing.）。ttcrpy库包含了三种射线追踪方法：快速扫描算法（FSM）、最短路径法（SPM）、动节点最短路径法（DSPM）。包含其二维与三维的实现。
2022-09-12 16:00:48
1139
18
原创 基于ttcrpy(三维射线追踪)的跨孔CT高斯牛顿算法及python代码分享（1）
ttcrpy是加拿大学者伯纳德·吉鲁（Bernard Giroux）于2021年发布的开源python库，详见（https://github.com/groupeLIAMG），参考文献（Giroux B. 2021. ttcrpy: A Python package for traveltime computation and raytracing.）。ttcrpy库包含了三种射线追踪方法：快速扫描算法（FSM）、最短路径法（SPM）、动节点最短路径法（DSPM）。包含其二维与三维的实现。
2022-09-11 21:54:46
1502
1
原创 基于python的三维射线追踪库-ttcrpy详解（5）
基于python的三维射线追踪库-ttcrpy详解（5）ttcrpy库中包含了矩形网格和三角形网格两种网格剖分形式，对于三角网格，要利用python中的vtk库和pygmsh库，本文研究ttcrpy中三角网格射线追踪。文章目录基于python的三维射线追踪库-ttcrpy详解（5）1、vtk库1.1、vtk库的安装1.2、vtk库的使用2、pygmsh库2.1、pygmsh库安装2.2、pygmsh库使用1、vtk库1.1、vtk库的安装安装：直接在cmd中pip一下即可。（正常情况下会报错，可以
2022-05-07 13:41:47
1650
原创 基于python的三维射线追踪库-ttcrpy详解（4）
基于python的三维射线追踪库-ttcrpy详解（4）python库ttcrpy中有三种射线追踪算法，他们是：FSM（快速扫描算法）、SPM（最短路径算法）、DSPM（动节点最短路径算法）。本内容为对比这三种算法的特点。文章目录基于python的三维射线追踪库-ttcrpy详解（4）1、模型一：低速异常体1、FSM（快速扫描算法）2、SPM（最短路径法）3、DSPM（动节点最短路径算法）4、三种方法追踪得到的走时对比5、模型二、高速异常体6、SPM与DSPM计算时间对比7、python代码搬砖不易，走
2022-05-02 14:48:48
2206
3
原创 基于python的三维射线追踪库-ttcrpy详解（3）
基于python的三维射线追踪库-ttcrpy详解（3）继续研究ttcrpy二维射线追踪，实现多发射点，多接收点的射线追踪。文章目录基于python的三维射线追踪库-ttcrpy详解（3）1、模型一2、模型二3、模型三4、python代码此代码在spyder上运行搬砖不易，走过路过，点个赞可好1、模型一2、模型二3、模型三4、python代码# -*- coding: utf-8 -*-"""Created on Wed Apr 27 15:37:28 2022@au
2022-04-29 16:16:39
941
原创 基于python的三维射线追踪库-ttcrpy详解（2）
基于python的三维射线追踪库-ttcrpy详解（1）二维射线追踪文章目录基于python的三维射线追踪库-ttcrpy详解（1）1、射线追踪效果2、python代码1、射线追踪效果2、python代码# -*- coding: utf-8 -*-"""Created on Wed Apr 27 15:37:28 2022@author: 86159"""import ttcrpy.rgrid as rgimport numpy as npimport matplotli
2022-04-28 23:04:03
1590
原创 gprMax中多相随机介质探地雷达三维建模与模拟
实际地下介质是非均匀介质，但数值模拟时常常把介质当做均匀介质，难以对实际介质产生准确认识。常规gprMax建模都是均匀介质建模、规则形状建模，本博文给出一种gprmax中多相随机介质建模的方法。
2022-04-25 20:43:49
7430
30
原创 基于MSFM算法与最速下降法的射线追踪技术
基于MSFM算法与最速下降法的射线追踪技术文章目录基于MSFM算法与最速下降法的射线追踪技术1、射线追踪效果2、技术原理3、代码上面的代码是是主程序代码，包含的两个子函数在网上可以找到，代码较长，我整理成文件夹上传到CSDN资源，欢迎感兴趣的读者下载学习交流。搬砖不易，走过路过的，点个赞支持一下，谢谢！1、射线追踪效果为了方便展示射线追踪的过程，用matlab将射线追踪过程做成了gif动图的形式，直观的展示一下跨孔方式的射线追踪效果，如下图所示，左图是波前走时等值线图（波前走时场），右图是射线路径图。
2022-04-22 17:32:43
1845
30
原创 gprMax中任意不规则形状建模与模拟
gprMax中任意不规则形状三维建模与模拟gprMax中给出了不规则几何形状的建模方法。常规的几何建模都是基于gprMax内置的几何命令建模，所建立的模型是规则的，如圆球、正方体、圆柱体、三角体等。不规则形状三维建模需要利用matlab或者python生成三维数组，然后将三维数组保存到HDF5格式文件中，因此，需要利用matlab编程生成任意不规则三维数组。本帖子给出一个生成任意不规则形状的代码。文章目录gprMax中任意不规则形状三维建模与模拟1、任意不规则几何形状建模2、matlab生成不规则形状的
2022-04-18 22:00:55
6770
43
原创 利用MSFM算法计算程函方程的走时CT胖射线反演法
利用FMM算法计算程函方程的走时CT胖射线反演法基于程函方程的弹性波走时CT胖射线反演方法，是在反演迭代过程中，利用程函方程计算走时场，并得到接收点处的走时数据，用这个数据代替正演数据。每次迭代完成后，速度模型发生了改变，所以FMM算法计算的走时场也会发生改变，正演数据的精度很高。但是层析成像的雅克比偏导数矩阵用胖射线的结果来近似代替，这样避免了反向求解射线路径造成过大的计算量。这种算法的计算速度非常快，比最短路径法快很多，而且收敛效果优于常规胖射线和直射线，反演效果较好。一、速度模型的建立及正演为了
2021-12-14 16:27:07
3641
40
转载 快速行进算法（fast_marching_kroon）的matlab代码
快速行进算法（fast_marching_kroon）的matlab代码快速行进算法用于求解程函方程得到走时场，源代码转载于link文章目录快速行进算法（fast_marching_kroon）的matlab代码模拟结果模拟结果function [T,Y]=msfm(F, SourcePoints, UseSecond, UseCross)% This function MSFM calculates the shortest distance from a list of% points
2021-11-28 20:57:36
1760
8
原创 Matlab 读取 gprmax 的 out 文件详细解释
Matlab 读取 gprmax 的 out 文件规律详细解释。out 文件是hdf5格式的文件，普通软件无法读取。out 文件有一定的规律，在matlab平台上简单操作一下，就可以深刻理解 out 文件含义了文章目录Matlab 读取 gprmax 的 out 文件规律详细解释。一、常用函数简介1、hdf5read2、h5read3、h5info4、h5disp二、读取out文件三、out文件数据格式一、常用函数简介1、hdf5read我用的 matlab 版本为2020a，帮助文档提示不推荐使
2021-11-15 12:55:48
7378
19
转载 MATLAB中文显示乱码的问题如何解决
MATLAB中文显示乱码的问题如何解决最近发现不同版本的 matlab 写的代码不能互用，比如在 matlab 2020 上写的代码，在matlab 2016 上显示，中文注释全部变成乱码，可能导致别人以为你分享的代码不想让对方看到注释。为此本人分析了一下这个问题。文章目录MATLAB中文显示乱码的问题如何解决一、问题一、问题matlab很多函数在读取中文后显示乱码，为了显示中文，应改为UTF-8方式或其他支持中文的编码方式，这在Matlab中的操作为：slCharacterEncoding(‘U
2021-11-05 16:35:32
59030
5
原创 MATLAB代码实现射线追踪的快速方法
10行MATLAB代码解决直射线追踪问题文章目录10行MATLAB代码解决直射线追踪问题% 直射线,射线水平的情况 if zFS(i1)==zJS(i2) if mod(zFS(i1),dz) == 0 if zFS(i1) == ymax dr(n,nx*floor(abs(zFS(i1)-1e-6)/dz)+1:nx*(ceil(abs(zFS(i1)-1e-6)/dz)))=delta/
2021-10-23 17:02:57
5237
11
转载 filedialog 函数用法
filedialog 函数用法文章目录filedialog 函数用法# -*- coding: utf-8 -*-"""Created on Sun Oct 17 16:21:07 2021@author: 86159"""import tkinter as tkfrom tkinter import filedialogroot = tk.Tk()root.withdraw()Folder_path = filedialog.askdirectory()file_path
2021-10-17 16:34:47
577
转载 varargin 函数用法
varargin 函数用法很多工具箱函数都用到 varargin 这个输入参数，当使用 varargin 作为函数的输入时，可以输入任意多的参数，且函数不会报错，使我们的函数更加的适用。文章目录varargin 函数用法测试 varargin 函数测试 varargin 函数matlab 帮助文档对这个函数说明的非常详细。close allclearclc% 此程序用来测试matlab的varargin函数c = test(1,2);function [c,d] = test(a,
2021-10-17 16:32:41
1234
原创 探地雷达（GPR）的 C-scan （三维图）的绘制
探地雷达（GPR）的 C-scan （三维图）的绘制作为一个初学者，看到教科书上探地雷达 C-scan 图绘制的比较高级，很好看，看起来也比二维图直观。因此，本文分享一下我绘制 C-scan 图的方法。文章目录探地雷达（GPR）的 C-scan （三维图）的绘制一、C-scan 图二、解析三、MATLAB代码分享一、C-scan 图首先展示一下我做出的结果（模仿教科书中绘制的，数据来源于 gprmax 软件正演）。上面的是灰度图，把它变为彩色图为换成红（高值）、蓝（低值）显示的图模型图大
2021-10-15 20:26:21
4680
21
原创 运用经验模态分解方法（EMD）对实测探地雷达数据进行处理
运用经验模态分解方法（EMD）对实测探地雷达数据进行处理经验模态分解（EMD）是一种新型自适应信号时频处理方法，特别适用于非线性非平稳信号的分析处理。探地雷达实测数据存在干扰，因此运用经验模态分解（EMD)对实测数据进行处理非常有必要。文章目录运用经验模态分解方法（EMD）对实测探地雷达数据进行处理一、问题的提出一、问题的提出运用上篇博客时间增益的方法对一个金属管线的GPR实测数据进行处理结果如下：...
2021-10-06 19:53:15
1326
2
原创 gprMax 正演模拟中Ex、Ey、Ez三个分量之间的关系分析
gprMax 正演模拟中Ex、Ey、Ez三个分量之间的关系分析在 GPR 应用中，电场分量通常是测得量。我们一般的正演模拟用哪个电场分量呢文章目录gprMax 正演模拟中Ex、Ey、Ez三个分量之间的关系分析一、电场E、磁场B和波传播方向K之间的关系二、模拟结果分析一、电场E、磁场B和波传播方向K之间的关系电场 E、磁场 B、波传播方向 K 是两两正交的，如下图所示一般的天线是电偶极子，激发产生的电磁波传播方向与电场方向垂直。二、模拟结果分析首先看 Ex 分量Ey 分量Ez 分量
2021-10-06 19:50:43
3539
7
原创 运用时间增益方法对 gprmax 模拟的探地雷达数据进行处理
运用时间增益方法对 gprmax 模拟的探地雷达数据进行处理GPR 对地下进行探测时，随着电磁信号向下传播，雷达信号的幅值往往衰减很快。与浅层目标的回波相比，深层回波信号的幅度很小。这些信号需要经过校正处理后才能同时显示出来，即对快速衰减的深层回波信号进行补偿处理。对时变幅度进行补偿在GPR应用中被称为时间增益或距离增益。文章目录运用时间增益方法对 gprmax 模拟的探地雷达数据进行处理一、时间增益效果1、原始的 GPR 图像2、时间增益后的 GPR 图像3、时间增益函数二、Matlab 代码分享一
2021-10-05 15:58:13
4255
8
原创 关于Gprmax正演模拟结果显示空白的原因分析
Gprmax正演模拟结果显示空白文章目录Gprmax正演模拟结果显示空白一、问题的提出二、问题分析三、问题解决一、问题的提出用Gprmax正演模拟时，偶然发现模拟出来的结果是空白的，结果如下：但是波形图却有结果：虽然波形不是很明显，但是可以看出部分结果。这是我用来模拟的In文件#title: 3Dkondon_100cm_20cm#material: 12.0 0.01 1.0 0.0 soil#material: 1.0 0 1.0 0.0 air#domain: 4.000
2021-09-30 20:35:21
2111
28
简单而简短的CPTU测试结果解释应用
2024-04-12
基于MATLAB的波动方程逆时偏移程序
2024-01-28
Traveltime tomography of crosshole radar data without ray tracin
2023-12-09
Structural joint inversion on irregular meshes
2023-12-09
Structural joint inversion of time‐lapse crosshole ERT and GPR
2023-12-09
跨孔电磁波与地震波CT交叉梯度联合反演算法研究及应用.pdf
2023-12-09
基于PYTHON的Marmousi模型
2022-12-21
最短路径法弯曲射线追踪的MATLAB代码
2022-12-19
一种基于gprmax的多相随机介质探地雷达三维建模与模拟
2022-04-25
基于MSFM算法与最速下降法射线追踪技术
2022-04-22
gprMax中不规则异常体建模及模拟
2022-04-20
探地雷达的 C-scan 堆积图是这样画的吗？
2021-10-13
gprMax正演模拟为什么会出现（如图）二次发射波？
2021-10-02
TA创建的收藏夹 TA关注的收藏夹
TA关注的人