重磁实验二.docxVIP

《重、磁资料处理与解释》上机实验报告实验二：频率域位场处理和转换实验姓 名：学 号：专 业：地球物理学指导教师：王万银、纪晓琳完成时间：2017.1.101基本原理当已知实测平面的异常时，换算场源以外的异常称之为延拓，分为向上延拓和向下延拓。?半空间狄利克莱问题解析解：?其中：称为延拓因子，为计算面Z坐标，Z轴向下为正方向，为计算面频率域位场，为延拓面的频率域位场。2 输入/输出数据格式设计2.1 输入/输出数据文件名输入数据和输出数据文件名均保存在“parameter.txt”中。第一行为输入的低高度观测面数据文件；第二行为输出的高高度观测面数据文件；第三行～第四行依次为输入的扩边比例因子和延拓高度。A20_mag.grdA53_mag.grd1.53.32.2重要变量名filename_Field:低高度观测面数据文件filename_Conti:高高度观测面数据文件Field(m,n): 低高度观测面数据Conti(m,n): 高高度观测面数据error: 延拓后的均方误差factor_x: 扩边比例因子（1.0）height: 延拓高度（0：向上延拓，0：向下延拓）Factor_Conti: 延拓因子point: 点数line: 线数m: 扩边以后总点数（满足2的幂次方）n: 扩边以后总线数（满足2的幂次方）3 总体设计I输入有关参数：filename_Field,filename_Conti,factor_x,height从文件输入有关参数：point,line计算有关参数：m,m1,m2,n,n1,n2从文件输入有关数据：xmin,xmax,ymin,ymax,dx,dy,FieldP进行扩边处理傅里叶正变换计算延拓因子进行乘积运算傅里叶反变换O输出计算结果：Conti,error4 测试结果5 结论及建议由高高度观测面位场等值线图（图2）可以看出，向上延拓的效果较好，能清晰地反映磁异常位置，等值线光滑，均方误差仅为1.660418。与低高度观测（图1）相比，向上延拓3.3m的高高度观测（图2）得到的磁异常幅度有所降低，位场等值线图明显对1、2号这两个局部异常体进行了压制，突出了3号区域异常体。本次程序设计过程中的一些认识：扩边参数m,n计算完成后再开辟数组，否则会导致在扩边过程中数组空间不够用；二维扩边可利用一维扩边子程序，将大大减轻二维的繁琐程度；对老师提供的FFT和计算圆频率的子程序要理解输入/输出参数的含义并灵活运用；理解延拓因子的计算公式中各个参量的含义才能进行程序实现；最大的感受就是在一开始虽然对整体的程序流程有了大致的把握，但落实到具体的程序编写时却发现了对各种细节理解的偏颇之处，正所谓纸上谈兵就是如此吧，只有在实践中才能真正地理解和提升。附录：源程序代码!******************************************************************c!! 功能：频率域位场延拓!! 文件参数说明：! filename_Field:低高度观测面数据文件! filename_Conti:高高度观测面数据文件!! 输入参数说明：! factor_x: 扩边比例因子（1.0）! height: 延拓高度（0：向上延拓，0：向下延拓）! point: 点数! line: 线数! xmin,xmax: 点坐标的最大值，最小值! ymin,ymax: 线坐标的最大值，最小值! dx: x 方向点距! dy: y 方向线距!! 扩边参数说明：! m1,m2: 实际数据起点位置和终点位置(m2-m1+1=point)! m: 扩边以后总点数（满足2的幂次方）! n1,n2: 实际数据起点位置和终点位置(n2-n1+1=line)! n: 扩边以后总线数（满足2的幂次方）!! 输入数据说明：! Field(m,n): 低高度观测面数据! Freal(m,n): 低高度观测面数据数据实部! Fimage(m,n): 低高度观测面数据数据虚部!! 延拓参数说明：! U(m): x