![](https://img-blog.csdnimg.cn/20200613153042322.jpeg?x-oss-process=image/resize,m_fixed,h_224,w_224)
地震勘探
地震勘探技术
吃最辣的锅
这个作者很懒,什么都没留下…
展开
-
Madagascar v3.2 绘图: shell脚本示例
【代码】Madagascar v3.2 绘图shell脚本示例。原创 2023-11-13 10:32:16 · 63 阅读 · 0 评论 -
使用SU画波数谱
suaddhead < pp.dat n1=234 n2=663 |sushw key=dt a=1000| suspeck1k2 | suximage cmap=rgb1 legend=1 perc=99原创 2021-11-29 15:38:24 · 410 阅读 · 1 评论 -
使用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=原创 2021-11-26 21:18:11 · 409 阅读 · 0 评论 -
用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=原创 2021-11-09 16:32:41 · 323 阅读 · 0 评论 -
SU应用场景3:地震记录抽稀
# 从in.su每隔20道抽出一道存到out.susuwind < in.su key=tracl j=20 >out.su原创 2021-07-19 19:40:23 · 284 阅读 · 0 评论 -
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...原创 2021-07-18 23:05:24 · 590 阅读 · 0 评论 -
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原创 2021-05-17 11:32:40 · 277 阅读 · 0 评论 -
使用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原创 2021-04-09 14:31:17 · 1150 阅读 · 0 评论 -
使用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= 标题原创 2021-04-09 14:19:27 · 1622 阅读 · 0 评论 -
基于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我手上原创 2021-04-07 11:17:19 · 608 阅读 · 3 评论 -
使用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原创 2021-03-15 20:33:42 · 1965 阅读 · 1 评论 -
使用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原创 2021-03-15 20:32:35 · 602 阅读 · 2 评论 -
使用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,原创 2021-03-15 20:31:37 · 570 阅读 · 1 评论 -
判断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原创 2021-03-15 20:29:18 · 676 阅读 · 0 评论 -
如何基于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原创 2020-11-26 10:36:34 · 597 阅读 · 0 评论 -
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.拷贝该库文件如果能在本地找到该库文件,则将其拷贝到运行程序的路径如果本地没有该库文件,则从其他计算机拷贝...原创 2020-06-16 17:03:35 · 1720 阅读 · 0 评论 -
用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=纵坐标起点原创 2020-06-15 16:21:24 · 3487 阅读 · 2 评论