haoNote

原创 用matplotlib绘制成像剖面

【代码】用matplotlib绘制成像剖面。

原创 如何建立一支小而精的研究团队

在管理学和心理学视角下，建立一支 “小而精” 的研究团队需兼顾结构设计、成员匹配、动力激发和文化塑造。

原创 地震勘探成像技术在声学、医学影像与无损检测领域的应用

地震勘探成像技术基于波场传播理论、信号处理及反演算法，在声学、医学影像和无损检测（NDT）领域展现了显著的交叉应用潜力。以下分领域详述其技术移植与创新应用：

原创 2025年2月读文献：基于快速扫描法与伴随状态技术的透射与反射联合走时层析成像

基于快速扫描法与伴随状态技术的透射与反射联合走时层析成像

原创 程函方程的详细推导

以下是基于非均匀介质弹性波方程（无纵波假设）推导程函方程的详细过程，完整考虑纵波（P 波）和横波（S 波）的耦合效应：

原创 速度分析是怎么做的？叠加速度是波传播速度吗？

速度分析的方法与叠加速度的本质解析...

原创 什么是均方根速度（Root Mean Square Velocity）

均方根速度（RMS 速度）是地震勘探中用于描述 层状介质 中地震波传播速度的等效参数，是叠加速度的理论基础，也是连接叠加速度与真实层速度的桥梁。以下从定义、计算方法和实际应用三方面详细解析...

原创 如何从共炮点数据构建速度场，为射线层析或全波形反演提供初始模型（DeepSeek R1回答）

为了从原始的共炮点数据构建初始速度场，为射线层析或全波形反演提供初始模型，需按照以下步骤进行...

原创 用python生成动图

【代码】用python生成透明背景动图。

转载 矢量算符运算法则

来源：

原创 2023中国科学院文献情报中心期刊分区表（GEOCHEMISTRY & GEOPHYSICS 地球化学与地球物理）

原创 MATLAB计算一元二次方程ax^2+bx+c=d并保留5位小数

【代码】MATLAB计算一元二次方程ax^2+bx+c=d并保留5位小数。

原创 zip分卷压缩与解压

【代码】zip分卷压缩与解压。

原创 Madagascar v3.2 绘图: shell脚本示例

【代码】Madagascar v3.2 绘图shell脚本示例。

原创 MATLAB绘制波数谱

【代码】MATLAB绘制波数谱。

原创 Linux解压rar文件

【代码】Linux解压rar文件。

原创 使用SU画波数谱

suaddhead < pp.dat n1=234 n2=663 |sushw key=dt a=1000| suspeck1k2 | suximage cmap=rgb1 legend=1 perc=99

原创 使用SU做波数域滤波

# 滤出低波数suaddhead < pp_lap_pow.dat n1=234 n2=663 | sukfilter k=0.,0.3,0.7,1 amps=1,1,0,0 verbose=1 |sushw key=dt a=1000 | sustrip> pp_lap_pow_lk.dat# 滤出中波数suaddhead < pp_lap_pow.dat n1=234 n2=663 | sukfilter k=0.,0.3,0.7,1 amps=0,1,1,0 verbose=

原创 用SU保存矢量图

用SU保存矢量图img=$1psimage < ${img}.dat n1=240 \label1="Depth(km)" label2="Trace" labelfont=Times-Roman labelsize=24 \f1=0 d1=0.055 f2=0 d2=1 perc=99 \legend=0 lstyle=vertright lbeg=-0.1 lend=0.1 ldnum=0.05 lgrid=dot \height=6.0 width=10.0 \verbose=

原创 SU应用场景3：地震记录抽稀

# 从in.su每隔20道抽出一道存到out.susuwind < in.su key=tracl j=20 >out.su

原创 SU应用场景2：地震记录重采样

SU应用场景2：地震记录插值问题描述命令问题描述现有地震记录640道(无道头)，每道256个采样点，采样间隔5ms.现在需要将采样点插值成1280个采样点命令suaddhead < input.dat ns=256 | sushw key=dt a=5000 >output1.susuresamp <output1.su nt=1280 dt=1000 >output2.susustrip < output2.su >output2.dat...

原创 SU应用场景1：从多个二进制炮记录各抽一道拼在一起最后画频谱

suaddhead < file1.dat ns=100 | suwind dt=0.001 key=tracl min=1 max=1 >f1.susuaddhead < file2.dat ns=100 | suwind dt=0.001 key=tracl min=2 max=2 >f2.susuaddhead < file3.dat ns=100 | suwind dt=0.001 key=tracl min=3 max=3 >f3.sucat f1.su

原创 使用Madagascar实现三维地震数据体的动态展示

#!/bin/shplot(){ segyFile=3D.segy segyDir=. imgDir=. title=3D ny= number Of Inline nx= number Of Crossline crl0= the first Crossline to show inl0= the first Inline to show gain=98 inlshow= the default inline to show

原创 使用Madagascar绘制剖面图

目录1.查看要绘制的segy文件的总道数，以确定绘制范围2.快速预览从第a道起的b道剖面（注意:a+b1.查看要绘制的segy文件的总道数，以确定绘制范围sfsegyread < 文件名.segy endian=1 | sfinn1 : 每一道的采样点数n2 : 总道数2.快速预览从第a道起的b道剖面（注意:a+b<n2）sfsegyread < 文件名.segy endian=1 | sfput d1=0.01 unit1=km label1=Depth title= 标题

原创 基于Ubuntu16.04的Anconda3安装及虚拟环境配置（python3.7 + cuda10.1 + cuDNN7.6 + tensorflow2.1.0）

一、安装Anconda在清华大学开源镜像站或者Anconda官网选择适合你机器的版本下载。我这里下载的安装脚本文件名为Anaconda3-2020.11-Linux-x86_64.sh运行安装脚本sh Anaconda3-2020.11-Linux-x86_64.sh输入配置环境变量sudo gedit /etc/profile打开文本五后在最后一行加上：export PATH=~/anaconda3/bin:$PATH一、查看显卡信息lspci | grep -i vga我手上

原创 使用python3读取segy文件道头信息

# code for read headers of each traceimport matplotlib.pyplot as mpimport numpy as np import sysimport structimport binasciifileName = "BAOFU_amp.segy"nTrace = 4nSample = 3001fSegy = open(fileName,"rb")fSegyHeader = {'tracr':[5,8] ,'ep':[17,2

原创 使用python3读IBM格式segy文件并绘制变密度剖面图

# code for read IBM segy format seismic data files import matplotlib.pyplot as mpimport numpy as np import sysimport structimport binasciidef ibm2ieee(ibm): if ibm == 0: return 0.0 sign = ibm >> 31 & 0x01 exponent = ibm

原创 使用python3读IEEE格式segy文件并绘制变密度剖面图

# code for read IEEE segy format seismic data files import matplotlib.pyplot as mpimport numpy as np import sysimport structimport binasciifileName = "2267x1401.segy"nTrace = 2207nSample = 1401fSegy = open(fileName,"rb")data = np.zeros((nSample,

原创 判断segy文件的编码格式是IBM还是IEEE

# code for judge a segy file IEEE or IBM format = 'nothing'fileName = "vp_marmousi-ii.segy"fSegy = open(fileName,"rb")fSegy.seek(3225,0) tempValue = fSegy.read(1)hexValue = tempValue.hex()decValue = int(hexValue,16)if decValue == 1: format = 'I

原创 如何基于ubuntu16.04完整安装Madagascar3.1.1

基于ubbuntu16.04安装Madagascar3.1.11.准备库2.下载安装包3.解压安装包3.安装4.测试1.准备库在你的终端执行以下命令：sudo apt-get install libxaw7-devsudo apt-get install libnetpbm10-devsudo apt-get install libtiff5-devsudo apt-get install libgd-devsudo apt-get install libplplot-devsudo apt

转载 Linux操作系统：如何实现vim的跨文件复制粘贴

用vim打开一个文件，例如：a.cpp在普通模式下，输入：":sp"（不含引号）横向切分一个窗口，或者":vsp"纵向切分一个窗口，敲入命令后，你将看到两个窗口打开的是同一个文件在普通模式下，输入：":e b.cpp"，在其中一个窗口里打开另一个文件切换到含有源文件（a.cpp）的窗口，在普通模式下，把光标移到你需要复制内容的起始行，然后输入你想复制的行的数量（从光标所在行往下计算），在行数后面接着输入yy，这样就将内容复制到临时寄存器里了（在 普通模式下ctrl+w，再按一下w，可以在两个窗口之间.

原创 Linux操作系统：如何解决cannot open shared object file: No such file or directory一类问题

最近在linux操作系统运行程序时出现了libiomp5.so: cannot open shared object file: No such file or directory这种报错，比较简单的解决步骤是：定位所需库文件的位置loacte libiomp5.so2.拷贝该库文件如果能在本地找到该库文件，则将其拷贝到运行程序的路径如果本地没有该库文件，则从其他计算机拷贝...

原创 用SeismicUnix（SU）读segy文件画地震剖面

1、读.segy生成.susegyread tape=文件名.segy verbose=1 ns=采样点数 conv=1>文件名.su2、用.su画变密度剖面图（红白蓝）suximage < 文件名.su legend=1 d2=tr.d2 d1=tr.d1 title=“图名” cmap=rgb1 &3.自定义时窗时：suximage < 文件名.su legend=1 d2=tr.d2 d1=tr.d1 title=“图名” cmap=rgb1 x1beg=纵坐标起点

原创 c++ 输出保留两位小数怎样实现

题目描述：班上有学生若干名，给出每名学生的年龄（整数），求班上所有学生的平均年龄，保留到小数点后两位。输入描述:第一行有一个整数n（1<= n <= 100），表示学生的人数。其后n行每行有1个整数，取值为15到25。输出描述:可能有多组测试数据，对于每组数据，输出一行，该行包含一个浮点数，为要求的平均年龄，保留到小数点后两位。示例1输入21817输出17....

原创 STL—map的使用

参考资料题目描述：给出一个01字符串（长度不超过100），求其每一个子串出现的次数。输入描述:输入包含多行，每行一个字符串。输出描述:对每个字符串，输出它所有出现次数在1次以上的子串和这个子串出现的次数，输出按字典序排序。示例1：输入10101输出0 201 21 310 2101 2#include <iostream>#include <...

原创 吃糖果

题目描述:名名的妈妈从外地出差回来，带了一盒好吃又精美的巧克力给名名（盒内共有 N 块巧克力，20 > N >0）。 妈妈告诉名名每天可以吃一块或者两块巧克力。 假设名名每天都吃巧克力，问名名共有多少种不同的吃完巧克力的方案。 例如： 如果N=1，则名名第1天就吃掉它，共有1种方案； 如果N=2，则名名可以第1天吃1块，第2天吃1块，也可以第1天吃2块，共有2种方案； 如果N=3，则...

原创 c++ 自学笔记：多态

c++ 自学笔记：多态1.什么是多态2.动态多态的原理1.什么是多态C++面向对象三大特征：封装、继承、多态多态分两类：静态多态：包括 函数重载、运算符重载，复用函数名动态多态：包括 派生类、虚函数实现运行时多态静态多态与动态多态区别：静态多态的函数地址早绑定——编译阶段确定函数地址动态多态的函数地址晚绑定——运行阶段确定函数地址动态多态满足条件：有继承关系子类重写...

原创 子数组异或查询

有一个正整数数组 arr，现给你一个对应的查询数组 queries，其中 queries[i] = [Li, Ri]。对于每个查询 i，请你计算从 Li 到 Ri 的 XOR 值（即 arr[Li] xor arr[Li+1] xor … xor arr[Ri]）作为本次查询的结果。并返回一个包含给定查询 queries 所有结果的数组。示例 1：输入：arr = [1,3,4,8], qu...

原创 解码字母到整数映射

给你一个字符串 s，它由数字（‘0’ - ‘9’）和 ‘#’ 组成。我们希望按下述规则将 s 映射为一些小写英文字符：字符（‘a’ - ‘i’）分别用（‘1’ - ‘9’）表示。字符（‘j’ - ‘z’）分别用（‘10#’ - ‘26#’）表示。返回映射之后形成的新字符串。题目数据保证映射始终唯一。示例 1：输入：10#11#12输出：jkab示例 2：输入：1326#输出：acz...

原创 内存分区模型

内存分区模型1.总览2.程序运行前3.程序运行后4.new操作符1.总览代码区：存放函数体的二进制代码，由操作系统进行管理的区域全局区：存放全局变量和静态变量以及常量的区域栈区 ：由编译器自动分配释放，存放函数的参数值、局部变量等堆区：由程序员分配和释放，若程序员不释放，程序结束时由操作系统回收不同区域存放的数据，生命周期不同2.程序运行前在程序编译后，生成了可执行程序（如Wi...