频率域位场处理与转换(1225-).docxVIP

重磁实验报告 实验名称：频率域位场处理与转换 学生姓名：刘国强专业名称：勘查技术与工程学生学号：201226020126指导老师：巨鹏 王万银目录一．实验内容与要求2（1）实验内容：2（2）要求：2二．基本原理2（1）空间域位场延拓：2（2）频率域位场延拓：3三.输入输出数据格式设计31.输入数据：32.重要变量名3四．总体设计4五．测试结果5六．结论与建议6一．实验内容与要求实验内容：对一个网格化数据进行向上延拓，得到延拓后的结果，并计算延拓后的均方误差。要求：进行软件设计、代码编写、软件测试、编写实验报告。二．基本原理分别介绍空间域和频率域位场延拓的计算公式。（1）空间域位场延拓：延拓分为向上延拓和向下延拓。向上延拓是只将实测平面上的数据换到平面之上的平面上的计算。对于二度体（令z向下为正）： U(x,z)=- （z0）其中U（，0）为剖面上各点的实测值。Z为延拓的高度，即转换后的平面和观测平面的距离，由于z坐标向下为正，因此z0。空间域的延拓实际是积分的计算。 向下延拓是由实测磁场向磁源方向延拓。 （2）频率域位场延拓：设场源位于z=H平面以下（H0），则磁场为在z=H平面以上对x、y、z连续的函数，若z=0观测平面的磁场T（x，y，0）为已知，则由外部的狄利克莱问题： T（x，y，z）=dd 对其进行变量x，y进行傅里叶变换成 因此：= 在频率域中，延拓变成了对观测数据的傅里叶变换乘以延拓因子。三.输入输出数据格式设计1.输入数据：（1）建立文件：wave.par（2）观测面高度之差：height=3.3m（3）低高度网格化数据：A20_mag.grd（4）特征值：（6）输出文件名：expound_a20_2D.grd2.重要变量名（1）输入文件：input_filename(2) 输出文件：out_filename(3) 点数：mpoint(4)线数：nline(5)异常值：filed四．总体设计结束调用子程序，输出文件调用子程序，计算延拓因子调用子程序，对扩边数据进行快速傅里叶变换调用子程序，进行扩边调用程序，读入A20_mag.grd场值调用程序，计算扩边点位调用子程序，读入A20_mag.grd数据开始调用子程序，读入wave.par文件 五．测试结果A53_mag.filenameA20_mag_filenameExpound_A20_mag_filenameme分析：将底高度向上延拓3.3m后结果与高高度原场值图形近似，但是场值变小了，且异常部位相对不明显了。六．结论与建议结论：向上延拓可以压制浅部的重磁异常特征，突出深部异常。向下延拓是由实测磁源异常得到。建议：本次实验程序编写比较困难，以后要提高编程能力。最后得到了扩边后的数据，分析能力比较差，分析结果不够准确，多看资料，加强这方面能力。七．源程序及代码PROGRAM wavenumber_continuationINTEGER mpoint,line,m0,m1,m2,m3,n0,n1,n2,n3,m,nREAL height,eigval,Xmin,Xmax,Ymin,Ymax,dx,dyCHARACTER*80 input_filename,output_filename,cmd_filename,Expound_a20_2D_filenameREAL,ALLOCATABLE::Field(:,:),Freal(:,:),Fimage(:,:)REAL,ALLOCATABLE::U(:),v(:),w(:,:),Fcount_Real(:,:),Fcount_image(:,:)cmd_filename=wave.parCALL Read_CMD(cmd_filename,height,input_filename,Expound_a20_2D_filename,output_filename,eigval)CALL Read_mn(input_filename,mpoint,line,Xmin,Xmax,Ymin,Ymax)CALL Calculate_m1m2_dx(mpoint,m0,m1,m2,m3,m,Xmin,Xmax,dx)CALL Calculate_m1m2_dx(line,n0,n1,n2,n3,n,Ymin,Ymax,dy)ALLOCATE(U(m0:m3),V(n0:n3),W(m0:m3,n0:n3),Fcount_Real(m0:m3,n0:n3),Fcount_image(m0:m3,n0:n3))ALLOCATE(Field(m0:m3,n0:n3),Freal(m0:m3,n0:n3),Fimage(m0:m3,n0:n3)) CALL Input_GRD(input_filename,m0