实验二：三度直立六面体正演程序设计实验报告教程.docxVIP

《重磁资料处理与解释》实验二 直立六面体正演程序设计 专业名称： 地球物理学 学生姓名： 学生学号： 指导老师：王万银、纪新林、纪晓琳、邱世灿 提交日期：2016-11-30 目 录  TOC \o 1-3 \h \z \u  HYPERLINK \l _Toc465871370 1 基本原理  PAGEREF _Toc465871370 \h 1  HYPERLINK \l _Toc465871371 2 输入/输出数据格式设计  PAGEREF _Toc465871371 \h 1  HYPERLINK \l _Toc465871372 2.1 场源参数数据格式  PAGEREF _Toc465871372 \h 1  HYPERLINK \l _Toc465871373 2.2 计算点坐标数据格式设计  PAGEREF _Toc465871373 \h 2  HYPERLINK \l _Toc465871374 2.3 计算结果输出数据格式设计  PAGEREF _Toc465871374 \h 2  HYPERLINK \l _Toc465871375 2.4 参数文件数据格式设计  PAGEREF _Toc465871375 \h 2  HYPERLINK \l _Toc465871376 3 总体设计  PAGEREF _Toc465871376 \h 2  HYPERLINK \l _Toc465871377 4 测试结果  PAGEREF _Toc465871377 \h 3  HYPERLINK \l _Toc465871378 4.1 测试参数  PAGEREF _Toc465871378 \h 3  HYPERLINK \l _Toc465871379 4.2 测试结果  PAGEREF _Toc465871379 \h 4  HYPERLINK \l _Toc465871380 5 结论及建议 6  HYPERLINK \l _Toc465871381 附录：源程序代码 7  PAGE \* MERGEFORMAT29 1 基本原理 在空间直角坐标系o-xyz中，形体（直立六面体）模型如图1所示。设该直立六面体x方向的坐标范围为，y方向坐标为，z方向（铅垂向下为正）坐标为；又设该直立六面体剩余密度为， 根据正演理论得知，其在空间任意一点处产生的重力异常为 （1-1） 式中，为万有引力常数，在国际单位制中其值为；为计算点到场源点的距离，可表示为 （1-2） 2 输入/输出数据格式设计 2.1 场源参数数据格式设计 场源参数按照一个直立六面体为一个记录进行设计，在数据文件中占一行。第一列为剩余密度density_source(g/cm3)；第二列为磁化强度mag_source(nT)；第三列为磁化方向倾角Az_source (DEG)；第四列为磁化方向与x轴的夹角Ax_source (DEG)；第五列为磁化方向与y轴的夹角Ay_source (DEG)；第六列～第七列为x坐标的起点X1和终点X2(km)；第八列～第九列为y坐标的起点Y 1和终点Y2 (km)；第十列～第十一列为z坐标的起点Z1和终点Z2)(km，向下为正)。以上各量均为实型变量，各量的意义见图1所示。 2.2 计算点坐标数据格式设计 计算点坐标数据格式设计为1：曲面非规则网；2：曲面规则网；3：平面非规则网；4：平面规则网；共四种形式。其中非规则网采用一个计算点为一个记录的方式设计。第1列保存计算点x坐标x_coord，第2列保存计算点y坐标y_coord，第3列保存计算点z坐标z_coord。以上各量均为实型变量。规则网采用grd文件形式进行设计。 2.3 计算结果输出数据格式设计 计算结果输出数据格式与输入格式对应。非规则网时，采用一个计算点为一个记录的方式设计。第1列保存计算点x坐标x_coord，第2列保存计算点y坐标y_coord，第3列保存计算点计算结果field。以上各量均为实型变量。规则网时采用grd文件形式进行输出。 2.4 参数文件数据格式设计 将以上部分量保存在一个文件中，该文件名变量为cmdfile，字符串变量，长度不超过80，全路径名。在该文件中保存的参数如下： 场源参数文件名filename_source，字符串变量，长度不超过80； 计算点坐标文件名filename_obser，字符串变量，长度不超过80； 计算结果输出文件名filename_field，字符串变量，长