资源描述：
《长安大学重磁报告二》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。
1、重磁实验报告（二）位场边缘识别姓名:学号:专业：勘查技术与工程指导教师：鲁宝亮、王万银完成日期：2013年12月10日基本原理:1、重力异常计算公式.2、磁力异常计算公式四、五、目录输入/输出数据格式设计:1、输入数据文件名的格式设计2、输出数据文件名的格式设计3、重要变量的名称总体设计:测试结果:11233345结论及建议:错误!未定义书签。错误!未定义书签。附录：源程序代码Efff<-z）物^如今叫J松一彳埔一力2+（卜z）2f2式中：V为地质体的剩余质量对A点的单位质量产生的引力位；V为地质体的体积。我们还
2、可以推导出计算重力异常垂向梯度或重力垂向梯度异常的基本公式为：=v=r>Tfff2((—z)2—(f—x)2—(7—y)2dzJJJ[(^-%)24-(Z7-J)24-(^-z)2]5/2d^d/jdC计算重力异常水平梯度或重力水平梯度异常的基本公式为:(fz)(f—X)[0-又)2+("-y)2+((-z)2]5/2(^-z)(Z7-y)d^dTjdC[(^-x)2+(7-y)2+(^-z)2]5/2d^drjd^一阶导数类边缘识别计算公式：垂向导数：VDR（x,z）=3Ag（；V°?，Z）dz总水平导数：TH
3、DR(x,y,z)=.^S^y.z)}2+z3A^(x,y,z)dx解析信号振幅：ASM=JTHDR2+VDR2(2)(x,y，z)/at)/J“0，)’0，z0)Lz//zc直立长方体示意图均匀磁化直立长方体磁场各分量表达式:H(„(x,>',z)^{-Mvarctan(^-x)(r]-y)(^-x)2+r(^-z)+(^-z)2aIbIc(1)Mvlnk+(r_z)]+M2lnk+(7_y)l}广+!Z°12-V°~iZ°~2H“y(x，y,z)=^{Mx[r+(^-z)]-Mvarctan(7—y)2
4、+r(f—z)+(f—z)2aabX~2K-三
5、+MJn[r+(^-x)]}♦22—(2)Z“(x，y，z)=^{Mxln
6、r+(Z7—y)l+Mvln[r+(^-x)
7、4兀Mzarctan(£zM}abx.+一a2•v々abv—一。2^2e•+2e"2(3)把(1)、(2)、(3)式代入磁异常AT计算公式屮得：AT(%,>z)=^M{k^[r+^-x)]+女2•Ink+(77-y)]+&•ln[r+((-z)J+k4^arctan+众6•arctan(^-x)2+r(^-z)+(^-z)2>"I7I+k
8、5•arctan(7-y)2+r(^-z)+(^-z)2ab•、勹y^ClbX''-2门1厂((—z)其中r=[(
9、-x)2+(7-y)2+((-2)2]1/2,么77,(为异常体的坐标，Lo、Mo、为地磁场的方向余弦；贷为总磁化强度，M。、My、Ma为其在x、Y、Z坐标方向的分量。设M为H的模，I为磁化强度倾角，D为磁化强度磁偏角，且满足关系式Ma==McosIcosDMy=McosIsinD二、输入/输出数据格式设计:1、输入数据文件名格式设计：（1）输入数据及文件名对应如下：inputgravity_
10、fi1ename:anomaly,grdsquare_num:0i叩ut_gravity_filename对应anomaly,grd文件，存放重力异常数据；square_num赋值力0表不•将区域扩力256*256的方形区域。（2）观测面上的重力异常存放在“anomaly,grd”中，坐标单位为m，格式如下：2011000.0000001000.000000fieldmaxDSAA201-1000.000000-1000.000000fieldmin可知，网格化数据按GRD格式保存，有201条线、每条线201个点
11、，点（线）坐标从一1000m〜1000m，点（线）距为10m，fieldmin,fieldmax分别表示网格化数据的最小值与最大值。2、输出数据文件名格式设计：(1)输出数据及文件名对应如下：verticalderivativefilename:verticalderivative.grdtotal_horizontal_derivative_filename:total_horizontal_derivative.grdsignalamplitudefi1ename:signalamplitude.grd(2)
12、对给定平面上的重力异常，计算其垂向一阶导数存放在vertical_derivative.grd文件中；总水平导数存放在total_horizontal_derivative.grd文件中；解析信号振幅存放在signal_amplitude.grd文件中。坐标单位为m，格式如下：DSAA2011000.0000001000.000000fieldmax201-1000.00
