地球物理重磁勘探实验实验名称：直立六面体正演程序设计实验 学院名称：专业名称：姓名：学生学号：指导老师：目录一，实验内容：2二，基本原理3(1) 模型示意及有关变量的意义，并分别给出重力异常和磁力异常的计算公式，化极。3(2) 方向描述及磁化强度分量计算4(3) 计算点输入坐标格式设计4(4) 计算点输入坐标格式设计5三，输入/输出数据格式设计5(1) 参数5(2) 类型5四，总体设计6五，测试结果6(1) 命名文件：filename.txt7(2) 场源参数文件：source.txt7(3) 平面地形数据文件：bxyz53.txt8(4) 输出的重力异常文件：gravity.dat8(5) 输出的磁力异常文件： magnetic.dat9(6) 化极磁力异常文件：huaji.dat9(7) 重力异常，磁力异常，化极磁力异常由surfer得到的图形10六，结论及建议12附录：源程序代码(本程序按照自由格式编写，对于粘贴中难免有一些问题)12一，实验内容:对给定的密度体和磁性体计算平面规则网上的重力异常、磁力异常以及化极 磁力异常。(1) 曲面规则网：用网格数据(.grd)的格式输入(2) 平面规则网：也用网 格数据( grd)的格式输入DSAADSAAAffmpffint), N(Une)M,Nmax二，基本原理(1)模型示意及有关变量的意义，并分别给出重力异常和磁力异常 的计算公式，化极。直立点元(直立六面体法三度体正演Green基本函数垂向一阶导数表达式r（x9yfz） = f ftr兀(3=J J J宜$呵可叱爲Jrur-亠= -fe-x）ln（7-y）+r+（7-j）lQfe-xhr-（-j）arctanJA空二竺竺+竺竺+空空% dx at0 dY at0 dz at04GcrV - JJJ dxdydzAg = GoVZPTP 二 Z ,a丄(2)方向描述及磁化强度分量计算Y吩0倾角应该满足二-90c 7 =JsinJu2)平面规则网：也用网格数据( grd)的格式输入(3) 计算点输入坐标格式设计(1)曲面规则网：用网格数据(.grd)的格式输入DSAAM(mpo iut) N(line)M,Nv a V min max(4)计算点输入坐标格式设计曲面非规则网(4)平面非规则网耳EH1耳Fx2.y2,zx29Y2y y 7三，输入/输出数据格式设计(1)参数场 源 系 数 设 计 ai_source,ax_source,ay_source,xl_source,x2_source,yl_source,y2_source ,zl_source,z2_source场源参数文件：sourcefilename 输入平面规则网文件：Inputfilename 重力异常输出文件：gravityfilename 磁力异常输出文件：magnaticfilename 化极磁力异常：RTPfilename 地磁的倾角：ai, 点数线数：mpoint，nline 平面规则网数据：xmin,(2)类型ax, aydensity ,Jmay ，xmax， ymin， ymax， zmin， zmax曲面规则网：输出用GRD格式平面规则网：输出也用GRD格式曲面非规则网：输出采用DAT格式，形式如下:X , Y , A g1 iiX , Y , A g2 2 2平面非规则网:X , Y , Z , A g111 1X , Y , Z , A g2 2 2 2X , Y , A gmmm，总体设计开始五，测试结果测试参数如下:地磁场倾角为50，与x轴的夹角为85，与y轴夹角为5 平面坐标范围与曲面坐标范围相同，但其z坐标值为-5.3km。(1)命名文件：filenametxt(2)场源参数文件：sourcetxt UltraEdit - C:UsersLHDesktopdiergecxlsource.txtJ8文件(F)编辑(E)搜索插入(n)项目(P)视圏(V) tetCD列(L)宏(M)脚本(i)高级(A)窗口(W)帮助(H) 口叵 S零。曰己目I昌釦圍I吕I釦釦dl冷-IdilIIsource.txt打开资源管理器丽匕卜El0 I I I I I I I I .1 卩I I I I I I1 0.7,20000,50,85,52 0.8,34000,50,85,53 0.9,17000,50,85,5i i 2卩i i i i i i i i 3卩i i i i i i i i 4卩i i i i i i i i 5卩i i i i i i i i .6 卩i r-3,-3,4,8,2,7F0F7p0r5r3r7t一4 8f-8r-3f 2f 7s-:s-:l-.:l-.s-.s-二二G:FTP帐号M a |uj 5 5= Al J % 1 1 EbF a ra ISP F録敢为？行1列1 coDOS修改：2O12/5/11 13:13:04 文件大小：98播平面地形数据文件：bxyz53txt防 Ultra Ed it - C:UsersLHDesktopdiergecxlbxyz53.txt卫文件(F)碍(E)搜索插入(n)项目(P)视囹(V)格式E列(L)宏(M)脚本(i)高级(A)窗口(W)耕助(H) 回富曰己冒|白皿團|圏|寻丨曲|訓右謹|bxyz53.txt资源管理器 环FF 筛选： * 反c D巳巳GH 二因豊 S-.S-.S-.S-.S-.S-0 I I I I I I I I 丄卩I I I I I I I I 2U I I I I I I I 3Pi1 DSAA2 27 273 -26 264 -26 265 5.3 5.34,0,5,0,6.0,勺遂塞逵卷旦一 阪赫曲F7 Z1 $録翠蚩？亀盅回隘與寿？ p w离：m “ ？行1,列1, coDOS修改:2012/5/11 13:13:56文件大小:38插(4) 输出的重力异常文件：gravity.dat筛选:-26.0000026.000003C：D:E:F:G：H：日曰m 9出电 :-;.ll-.ll-:.l-:.-:a-10111214151617192021-26.00000-8.508221-4.3975830E-02-8.1668086E-0201370148-0.1534061-0.1125061-6.1373446E-02-4.8933409E-02-9.6258543E-02-0.1740448-0.2067329-0.133803970291348E-02-5.4489136E-02-0.1139618-0.2251480-0.276757826.00000-4.9746703E-02-9.1928095E-02-0.1466003-0.1533073-0.1007938-5.3923462E-02一55937193E-02-0.1100632-0.1884720-0.1989502-0.1223658-6.1082266E-02-6.2802121E-02-0.1323868-0.2474976一02644031-5.6352232E-02-6.3836671E-02-7.2-0.1030273-0.1145248-0.1一01540558-0.1454430-8.9584157E-02-6.3987732E-02-0.1252075-0.1870106-0.1070859-0.1533508-0.2656265一02457157-0.1587837-0.135511873320009E-02-0.1413605-0.2072212-9.3219638E-02-8.4281921E-02-0.1764053-0.2773788-0.2228092-0.1-0.1-8.4-0.1-0.2-0.1-8.0-9.7-0.2-0.2-0.1-0.1731281-0.1497066-0.1285880-0.1100235-8.0490015E-0224-6.0366821E-02一70277400E-02-8.2310483E-02-0.12501350540-0.1892227-0.2226540-0.2-0.3824959-0.2973068-0.3817779-0.3614151-0.2943157-0.2|靄龟盘回隘讯錨？ W w昌I:翹* Q ?行1,列1, coDOS惨改:2012/5/12 1:07:10 文件大小：11387（5）输出的磁力异常文件：magnetic. dat（6）化极磁力异常文件：huaji dat重力异常，磁力异常，化极磁力异常由surfer得到的图形重力异常:-25-20-15-105101 520250-0.5-1-1.5-2-2.5-3-3.5-4-4.5-5-5.5-6-6.5-7-7.5-8-8.5-9磁力异常:420