石油软件概论-实验报告.docxVIP

本科生实验报告实验课程 石油地球物理软件技术概论 学院名称 地球物理学院 专业名称 勘查技术与工程（石油物探） 学生姓名 学生学号 指导教师 赵宪生 实验地点 地球物理学院5417 实验成绩 二〇一五年4月-二〇一五年5月目 录摘 要1第一章seg-y数据滤波21.1 方法21.2 源程序21.3 结果分析81.3.1 滤波因子81.3.2 原始数据及结果显示与分析9第二章 自动统计剩余时差校正132.1方法132.2源程序142.3结果分析19第三章 交互集成软件213.1方法213.2源程序213.3结果显示22参考文献23心得及评语24资料处理软件与交互软件设计摘 要地震勘探采集的数据需要经过一定的处理，才能被解释使用，地震勘探采集数据多为seg_y格式，石油勘探采集处理越来越依赖于计算机技术的发展，使用计算机技术实现地球物理处理软件设计、交互软件设计和图形软件设计及可视化和集成化基本技术。本实验报告中，主要针对seg_y格式数据进行滤波处理，比较带通滤波和低通滤波的结果差异；对给定信号进行剩余时差校正，并写出相应的C语言源程序；并在MFC框架程序下，设计出相应的交互集成软件。关键字：seg_y，带通滤波，低通滤波，剩余时差校正，交互集成软件第一章seg-y数据滤波1.1 方法分析已知seg-y数据格式的程序，补充该程序中计算滤波因子的子函数段，选择参数对数据文件366.sgy完成滤波处理，以及对某一道数据做谱分析，比较低通滤波和带通滤波及参数对结果的影响。用Fimage.exe程序显示原数据和滤波处理数据，比较分析结果。1.2 源程序#include stdio.h#include math.h#include string.h#include ctype.h#define LEN 3200#define NQ2 30000#define PI 3.1415926//==========低通滤波因子=========//float filter1(float h[],int n){int i;float fc;double t,dt=0.002;printf(输入fc: \n);scanf(%f,fc);for(i=1;i=n/2;i++){t=i*dt;h[i]=sin(2*PI*fc*t)/(PI*t);h[i+n/2]=h[i];}h[25]=2*fc;for(i=1;i=n/2;i++){h[n/2-i]=h[n/2+i];}}//============带通滤波因子=========// float filter2(float h[],int n){int i;float f1,f2,df,f0;double t，dt=0.002;printf(输入f1、f2: \n);scanf(%f%f,f1,f2);f0=(f1+f2)/2;df=abs(f2-f1)/2;for(i=0;i=n/2;i++){t=i*dt;h[i]=2*cos(2*PI*f0*t)*sin(2*PI*df*t)/(PI*t);h[i+n/2]=h[i];}h[25]=4*df;for(i=0;i=n/2;i++){h[n/2-i]=h[n/2+i];}}//==========求褶积函数==========// int cvbx0(float x[], float h[], float r[],int kc, int nt){int m,np,i,j;float sum;m=kc/2;for(np=0;npnt;np++)r[np]=0.0;for(i=0;int;i++){sum=0.0;for(j=0;jkc;j++){if((i-m+j)0||(i-m+j)=nt) continue; sum=sum+x[i-m+j]*h[j];}r[i]=sum;}}//================傅里叶函数==============//int IBIT(int j, int m){int i, it, j1, j2;it=0;j1=j;for(i=1; i=m; i++){j2=j1/2;it=it*2+(j1-2*j2);j1=j2;}return it;}void fft(float *xr, float *xi, int nr, float T){int i, j, k, l, n, n2, nr1, i1, j1, k2, k1n2;float c, s, p, tr, ti, trc, tic, ars;n=(int)(pow(2,nr));n2=n/2;nr1=nr-1;k=0;for(l=1; l=nr; l++){loop1: for(j