改实验一：二度直立柱体正演程序设计实验报告(202,24-李龙坤,2016-11-18)教程.docx

《重磁资料处理与解释》实验一 二度直立柱体正演程序设计 专业名称：勘查技术与工程 学生姓名：李龙坤 学生学号：201326020224 指导老师：王万银、纪新林、纪晓琳、邱世灿 提交日期：2016-11-18 目 录  TOC \o 1-3 \h \z \u  HYPERLINK \l _Toc468715707 1 基本原理  PAGEREF _Toc468715707 \h 1  HYPERLINK \l _Toc468715708 2 输入/输出数据格式设计  PAGEREF _Toc468715708 \h 1  HYPERLINK \l _Toc468715709 2.1 场源参数数据格式设计  PAGEREF _Toc468715709 \h 1  HYPERLINK \l _Toc468715710 2.2 计算点坐标数据格式设计  PAGEREF _Toc468715710 \h 2  HYPERLINK \l _Toc468715711 2.3 计算结果输出数据格式设计  PAGEREF _Toc468715711 \h 2  HYPERLINK \l _Toc468715712 2.4 参数文件数据格式设计  PAGEREF _Toc468715712 \h 2  HYPERLINK \l _Toc468715713 3 总体设计  PAGEREF _Toc468715713 \h 2  HYPERLINK \l _Toc468715714 4 测试结果  PAGEREF _Toc468715714 \h 3  HYPERLINK \l _Toc468715715 4.1 测试参数  PAGEREF _Toc468715715 \h 3  HYPERLINK \l _Toc468715716 4.2 测试结果  PAGEREF _Toc468715716 \h 4  HYPERLINK \l _Toc468715717 五 结论及建议  PAGEREF _Toc468715717 \h 4  HYPERLINK \l _Toc468715718 附录：源程序代码  PAGEREF _Toc468715718 \h 5 1 基本原理 1 2 输入/输出数据格式设计 1 2.1 场源参数数据格式 1 2.2 计算点坐标数据格式设计 2 2.3 计算结果输出数据格式设计 2 2.4 参数文件数据格式设计 2 3 总体设计 2 4 测试结果 3 4.1 测试参数 3 4.2 测试结果 4 5 结论及建议 5 附录：源程序代码 6  PAGE \* MERGEFORMAT13 图1 二度直立柱体示意图 1 基本原理 在二度直立柱体示意图，但不太确定。 直角坐标系o-xz中，形体（二度直立柱体）模型 如图1所示。设该二度直立柱体x方向的坐标范围为， z方向（铅垂向下为正）坐标为；又设该二度直立柱体剩余密度为。 根据正演理论得知，其在空间任意一点处产生的重力异常为 （1-1） （1-1） 式中，为万有引力常数，在国际单位制中其值为； 。 若有n个二度直立柱体，其引起的重力异常符合叠加原理，由此得到 （1-2） 2 输入/输出数据格式设计 2.1 场源参数数据格式设计 场源参数按照一个二度直立柱体为一个记录进行设计，在数据文件中占一行。第一列为剩余密度density_source(g/cm3)；第二列～第三列为x坐标的起点和终点起点X1和终点X2(m)；第四列～第五列为z坐标的起点和终点Z1和终点Z2)(m，向下为正)。以上各量均为实型变量，各量的意义见图1所示。例如： 0.2 -100 -50 50 200 0.3 -50 50 50 200 2.2 计算点坐标数据格式设计 计算点坐标数据格式设计为规则网测线，采用一个计算点为一个记录的方式设计。第1列保存计算点x坐标x_coord（m），第z列保存计算点z坐标z_coord（m）。以上各量均为实型变量。保存在Input_file_source文件名变量中。例如：用一个模型数据在此示意一下会更好。以下数据格式设计相同。 已经知道了 -400 -20 -398 -20 2.3 计算结果输出数据格式设计 计算结果输出数据格式与输入格式对应，设计为非规则网规则测线，采用一个计算点为一个记录的方式设计。