资源描述:
《石油软件概论报告-》由会员上传分享,免费在线阅读,更多相关内容在教育资源-天天文库。
1、本科生实验报告实验课程石油地球物理软件技术概论学院名称地球物理学院专业名称勘查技术与工程学生姓名学生学号201105060416指导教师赵宪生实验地点地球物理学院5419实验成绩二〇一四年十一月二〇一四年十二月资料处理软件与交互软件设计摘要地球物理软件技术发展表现出的先进性,开放性一体化,自由化,网络化,普遍化等特征。利用地震资料的横向相关性的特点,从分频处理的思想出发,把小波分析原理应用于提高地震资料信噪比和分辨率的处理上.通过对实际资料的处理表明,使用小波变换方法对地震资料进行处理,具有更为精细的特点,处
2、理效果有明显的改善.该文的重点在于小波变换地震资料处理软件包的设计与实现.采用自顶向下逐步求精的结构化方法,按照功能特点逐步细化,得出整个软件包的系统流程图、数据流图、层次结构图,最终将软件结构转化为软件过程化描述,实现每个模块的算法,使它们有机地结合成为一个高效系统伴随石油地球物理不断地向前发展,我们对于数据的收集量和处理量都急剧攀升,不论是科研还是实际勘探都面临着数据资料处理效率带来的压力。目前市面上的数据资料处理软件类目繁杂、价格昂贵,且软件与软件之间几乎无法实现数据的共享。关键词:石油地球物理软件;地
3、震资料;数据整理;数据流图第1章segy数据滤波1.1方法segy是指“segy格式地震数据”。标准SEG-Y文件一般包括三部分,第一部分是EBCDIC文件头(3200字节),由40个卡组成(例如:每行80个字符*40行),用来保存一些对地震数据体进行描述的信息;第二部分是二进制文件头(400字节)用来存储描述SEG-Y文件的一些关键信息,包括SEG-Y文件的数据格式、采样点数、采样间隔、测量单位等一些信息,这些信息一般存储在二进制文件头的固定位置上;第三部分是实际的地震道,每条地震道都包含240字节的道头信
4、息和地震道数据。道头数据中一般保存该地震道对应的线号、道号、采样点数、大地坐标等信息,但一些关键的参数位置(如线号、道号在道头中的位置)并不固定。地震道数据是对地震信号的波形按一定时间间隔Δt进行取样,再把这一系列的离散振幅值以某种方式记录下来[1]。地震数据格式可以是IBM浮点型、IEEE浮点型、整型、长整型等,一个三维地震工区同一次处理的地震数据格式是唯一的。地震道采样点数由该地震道道头中采样点数决定,大部分SEG-Y文件的所有地震道采样点数是一致的,但也存在不同地震道采样点数不同得情况,一般称这种SEG
5、-Y文件为变道长格式的SEG-Y文件。1.2程序//++++++++++++++++++++++++++++++++++++++++++++++++++++////读segy数据方法////++++++++++++++++++++++++++++++++++++++++++++++++++++//#include#include#include#include#defineLEN3200#defineNQ230000voidfilter(
6、float*h,intn){inti;floatt,f=40.,dt=0.001,pi=3.14;for(i=1;i<=n/2;i++){t=i*dt;h[i]=sin(2*pi*f*t)/(t*pi);//算出滤波函数的右半边h[i+n/2]=h[i];//将右半边向右移动25}h[25]=2*f;//定义在零点的值for(i=0;i<=n/2;i++){h[n/2-i]=h[n/2+i];//左半边由右半边翻转定义}}voidcvbx0(float*ki,float*h,float*ko,intn,int
7、len1){inti,j,k;floats;for(i=0;i8、
9、k>=len1)continue;//当j为0到24和len1-24到len1时,褶积不完全当(k<0
10、
11、k>=len1)时ki没有数与h褶积s=s+h[i]*ki[k];//褶积}ko[j]=s;//i每循环完一次赋给ko}}voidcvbx0(float*ki
12、,float*h,float*ko,intn,intlen1){inti,j,k;floats;for(i=0;i13、
14、k>=len1)continue;//当j为0到24和len1-24到len1时,褶积不完全。这时当(k<0
15、
16、k>=l