频率域位场处理和转换实验.docVIP

《重磁资料处理与解释》实验二 频率域位场处理和转换实验 :地测学院 :勘查技术与工程 :2018年1月9日 院称名号师期 名姓学老日 业生屯承六乂 学专学学指提 二0 —八年一月 目录 TOC \o 1-5 \h \z 1 穌腿 2 HYPERLINK \l bookmark3 \o Current Document \h 1.1位场的方程 2 HYPERLINK \l bookmark4 \o Current Document \h 1.2二维傅里叶变换及卷积性质 2 傅里叶变换 2 卷积性质 2 HYPERLINK \l bookmark5 \o Current Document \h 1.3频率域位场延拓原理 3 HYPERLINK \l bookmark6 \o Current Document \h 2输入/输出数据格式设计 3 HYPERLINK \l bookmark7 \o Current Document \h 2.1输入数裾格式设计 3 HYPERLINK \l bookmark8 \o Current Document \h 2.2输出数据格式设计 3 HYPERLINK \l bookmark9 \o Current Document \h 2.3参数文件数裾格式设计 3 3说、體i十 4 HYPERLINK \l bookmark11 \o Current Document \h 3.1频率域位场处理与转换的一般步骤 4 HYPERLINK \l bookmark12 \o Current Document \h 3.2软件总体设计结果流程图 4 4觀式關 5 HYPERLINK \l bookmark14 \o Current Document \h 4.1测试参数 5 向上延拓 5 向下延拓 5 HYPERLINK \l bookmark13 \o Current Document \h 4.2测试结果 6 HYPERLINK \l bookmark16 \o Current Document \h 5结论及建议 7 HYPERLINK \l bookmark17 \o Current Document \h 附录：源程序代码 8 1基本原理 1.1位场的方程 由场论知识可知，位场方程分为 W大类：有源的Possion方程以及无源的Laplace方程(V2f/=Q)。 手 ls=f:Laplace方程的第一边值问题((；|s=通常为Dirichlet 手 ls=f: 通常称为Nueman问题。若P点在S平面内称为内部问題，反之称 为外部问题。由唯一性定理可知，Dirichlet的内部和外部问题的解是唯一的，而 Nueman内部问题的解不是唯一的，有一常数差，但其外部问题解是唯一的。 外部问题的解的唯一性的原因：=0;I；?一 外部问题的解的唯一性的原因：=0; I；?一＞00 du dn r—oo 0。 无源区域位场可以表示为.? w（p） w（p） dn dn （1-1） -dedri（1-2） -dedri （1-2） ［（￡*-x）2+（77-y）2+（^- W（￡,Z7，^）*/?（x, y^-z） 1.2二维傅里叶变换及卷积性质 （1） 傅里叶变换 +00+00 TOC \o 1-5 \h \z G（uyv） = F［^（x,y）］ = J e，2^l，x+vy} dxdy （1-3） —00—00 g（x9y） = p~l［G（u,v）］= J Jg（w,v）^2/T＜+iv）^v （1-4） —00—00 （2） 卷积性质 Mgb，＞0* 厂（又，＞’）］=G（w,v）*P（w,v） （1-5） F1 ［G（n，v）* P（w, v）］ = g（x, y）^ p（x, y） （1-6） 1.3频率域位场延拓原理 当已知实测平面的异常时，换算场源以外的异常称之为延拓，分为向上延拓 和向下延拓。 半空间狄利克莱问题解析解： 傘（X，y,z）］=F［W（x,y^） W（义少乂- z）］ （1-7） =F［W（xyy^）］e^^^-z} 其中：为延拓因子，（为计算面Z坐标，Z轴向下为正方向, 为计算面频率域位场，为延拓面的频率域位场。 2输入/输出数据格式设计 2.1输入数据格式设计 观测面位场数据保存在filename_obser文件中，为.grd格式。 计算延拓误差时的精确场值文件保存在filename_obser2屮，为.grd格式。 2.2输出数据格式设计 实际计算得到的数据保存在filename_outputl和文件filename_output2中， 为.grd格式。 2.3参数文件数据格式设计 将所要读取的参数保存在一个文件中，该文件名变量为c