重力异常数据处理方法

第四章 重力异常的数据处理

布格重力异常反映了地壳内部物质密度的不均匀性，即从地表到地下

几十公里的地壳深部，只要物质密度横向发生变化，在地下不同的空间和范田内形成剩余质量，就可以引起地表的重力异常。 定性解释侧重于判断引起异常的地质原因，并粗略估计产生异常的地

质体的形状、产状及埋深等。

定量解释则是通过理论计算．对地质体的规模、形状、产状及埋深等

作出具体解答。

重力异常的推断解释的步骤：

① 阐明引起异常的地质因素 具体地说，就是确定异常是浅部因素还是深部因

素引起，是矿体还是构造或其它密度不均匀体(岩性变化、侵入体等)的反映。——定性解释

② 划分和处理实测异常 重力异常图往往是地表到地球深处所有密度不均匀体

产生的异常的叠加图象。为了获取探测对象产生的异常，需要将它们进行划分。不同的研究目的提取的异常信息不同，例如，矿产调查要提取队是矿体或没部构造产生的局部异常；而深部重力研究的目标正好相反，需要划分出的是反映地壳深部及上地幔的区域异常。

③ 确定地质体或地质构造的赋存形态 一是根据已知地质体或地质构造的

形状、产状及埋深等．研究它们引起的异常的特征，包括异常的形状、幅度、梯度及变化规律等。二是根据异常的形态及变化规律等，确定地质体或地

实验数据记录与处理

石灰石的密度2.7g/cm3 质量m1=240g 体积V1=88.89ml

水的密度1g/ cm3 质量m2=(12000-88.89)*1=11911.91 体积V2=11911.91 石灰水的密度=（11911.91+240）/(11911.91+88.89)=1.013 g/cm3 过滤面积A=0.0216m2

Δq=ΔV/A Δt=10s (1)当过滤压力P=40KPa时： 序号 1 2 3 4 5 6 7 8 9 10 滤液量M/g 11.8 9.2 7.4 6.2 5.4 4.8 4.2 3.9 3.5 3.2 滤液量V/ml Δq(m3/m2) 11.64856861 9.081934847 7.305034551 6.120434353 5.330700888 4.738400790 4.146100691 3.849950642 3.455083909 3.158933860 8000070000Δτ/Δq 18543.0508 23783.4783 29568.6486 35291.6129 40520.0000 45585.0000 52097.1429 56104.6154 6251

合成一个面

如果操作过程中提示你操作会失去原来的面或者线的时候，不妨把面或者线先copy，操作了之后再paste就好。

Solid 用来生成体。

第一栏用来直接生成一些规则的体。Sweep是通过旋转、拉伸面模型得到体。

第二栏是对体进行一些布尔操作，如加减等。Split是将一个体沿一个面（xy、yz、xz）劈开成两部分，可以选择要保留的部分。在减操作时，如有必要，还是先copy一下被减模型。

第三栏cover surface是通过闭合的曲面生成体。

Arrange 选取模型组件后，对模型组件进行移动、旋转、镜像（不保存原模型）、缩放等操作。 Options 用来进行一些基本的设置。单位的转换，检查两个体是否有重叠（保存的时候会自动检查）、设置background大小、定义公式以及设置颜色。

二、 材料设置

相对比较简单，Maxwell材料库自带了一些常用的材料，如果没有可以自己新建一个材料。Material—〉Add，输入文件名，及相关的参数即 可。如果BH曲线是非线性的，就，在

B-H Nonlinear Material 前面打勾，就会有自己输入BH曲线的选项，自己输入就好。但是要注意 BH曲线是单调递增的。

新建的材料还可以设置为理想导体和各向异性的材

6级5级异常数据核心基地

初始怪 1430 30 长明监狱长 第一波 1988 78 长明监狱长 第二波 3084 60 长明看护者 第三波 1580 96 核心堡垒

初始怪 2054 63 新生看护者 第一波 2944 48 苏醒的典狱官 第二波 3518 162 长明守卫 第三波 3252 144 明镜

初始怪 3334 54 第一波 3406 132 第二波 3168 84 奇怪的能源读数

初始怪 1742 39 第一波 2406 144 第二波 3366 78 核心兵站

初始怪 802 0 第一波 1542 72 第二波 1920 72 第三波 1402 33 核心要塞

初始怪 692 0 第一波 1928 24 第二波 756 18 第三波 2718 24 傲怒子欧

6级5级异常数据核心基地

初始怪 1430 30 长明监狱长 第一波 1988 78 长明监狱长 第二波 3084 60 长明看护者 第三波 1580 96 核心堡垒

初始怪 2054 63 新生看护者 第一波 2944 48 苏醒的典狱官 第二波 3518 162 长明守卫 第三波 3252 144 明镜

初始怪 3334 54 第一波 3406 132 第二波 3168 84 奇怪的能源读数

初始怪 1742 39 第一波 2406 144 第二波 3366 78 核心兵站

初始怪 802 0 第一波 1542 72 第二波 1920 72 第三波 1402 33 核心要塞

初始怪 692 0 第一波 1928 24 第二波 756 18 第三波 2718 24 傲怒子欧

Bernese5.0 PPP数据结算步骤

BERN的数据处理方法。基本上可以分为下面几个步骤： 一、处理数据的准备

此步骤包括准备观测文件、星历文件，以及更新数据处理所需的表文件，然后把RINEX格式的数据转化成Bernese二进制格式文件，目的为加速数 据读取速率。

RINEX格式的文件分别为观测文件(ssssdddf．yyo)、导航文件(ssssdddf．yyN，ssssdddf．yyG)和气象文件 (ssssdddf．yyM)。观测文件转换成BERNESE格式有如下四种格式，它们分别为： *.PZH(相位非差头文件) *.PZO(相位非差观测文件) *.CZH(码非差头文件) *.CZO(码非差观测文件)

导航文件和气象文件的转换与此相似。在原始文件由RINEX格式转换成BERNESE格式过程中，有时会出错，认为该接收机类型与Phas_igs．O1 文件不匹配，造成转换不成功。其主要问题是RINEX格式的原始文件中可能存在非法字符，该问题通常可以通过检查原始文件是否有非法字符或用数据管理软件 teqc使其标准化，该软件可以从www.unavco.ucar.edu下 载。 在运行该软件前，首先要准备好所必须的文件，具体包括：

(1)原始文

实验数据处理方法的合理选择

中国科技论文统计源期刊

（中国科技核心期刊）

实验技术与管理

ＥＸＰＥＲＩＭＥＮＴＡＬＴＥＣＨＮＯＬＯＣＹＡＮＤＭＡＮＡＧＥＭＥＮＴ

Ｖｄ．２３№．４

ＡＰＲ．２００６

实验数据处理方法的合理选择

李潮锐

（中山大学物理系．广州５１０２７５）

摘要：实验技术与教学资源的优化促进了教学内容的拓展和对实验结果的定量分析，因此对传统教学中所采用的定性或半定量分析方法的合理性产生质疑。通过具有教学普遍意义的“夫兰克一赫兹实验”中原子定态能级分析，指出常用的平均法或“类平均法”的局限性。选择数据拟台方法，不仅定量地获得实验测量结果及其误差，同时也为更深人了解实验的物理实质提供了合理的分析手段。

关键词：物理实验；数据处理；教学方法中圈分类号：Ｇ４２４．３１

文献标识码：Ａ

文章编号：１００２４９５６（２００６）０４一呻２３－０２

ＴｈｅＲｅａｓｏｎａｂｌｅＰｒｏｃｅｓｓｉｎｇｆｏｒＤａｔａＡｎａｌｙｓｉｓｉｎＰｈｙｓｉｃｓＥｘｐｅｒｉｍｅｎｔｓ

ＬＩＣｈａｏ．ｍｉ

（Ｐｈｙ西ｃｓ

Ｄ８ｐ８ｒｔⅡ１ｅｎｔ，ｓｕｎ

Ｙａｔ—ｓ即（ｚｂｏｎ础舶）ｕｎｉｖ嘲ｉｔｙ，Ｇ啪ｇｚｈｏｕ５１０２７５，ｃｈｉ∞）

ｑｕａｌｉ协ｔｉｖｅ粕ｄ∞ｍｉ—ｑＩｌ枷ｔ撕ｖｅｄ血ｐｒｏｃｅ鹤ｉ

用excel处理大学物理实验数据：

大学物理试验特点：好做，数据难处理，尤其是不确定度以及方差之类的东西。用计算器算太困难。下面本人介绍一种简单实际的计算机处理方法。（excel高手就不必看了） 1、 新建excel表格，根据自己的表格建立excel表格，此处不详细介绍了，不会可以上网找，

或者去借本书自己看。 2、 将原始数据填入表格。 3、 用公式和函数处理数据。

公式部分：根据你的原始数据输入公式求出要求的值，如下图：

图中“=B1*9.8”即为公式，B1是数字2的坐标，然后回车就OK了，接下来就到了excel的独特魅力的地方了，批量处理数据

当结果出来以后，把鼠标放在19.6那一格的右下角处，当出现“＋”符号时，点住鼠标左键向右拖，拖到第五个格时所有的数据都有了。

OK，会了吗？

介绍一下符号 乘：“*” 除：“/ ” 平方：“^2”即“5^2”就是求5的平方，同样3次方就把“^”后面的2改成3，开平方就改成0.5， 顺序先后就加括号就行了，比如：

处理复杂的批量运算很爽的!!

函数功能，更能体现excel的优势了， 这里介绍几个常用的函数把，

sum函数：求和，可以求一列的和或一行的和 average函数：求平均值。

st

题目：数据的预处理问题

摘要

数据处理贯穿于社会生产和社会生活的各个领域。数据处理技术的发展及其应用的广度和深度，极大地影响着人类社会发展的进程。数据补充，异常数据的鉴别及修正，在各个领域也起到了重要作用。 对于第一问，我们采用了多元线性回归的方法对缺失数据进行补充，我们将1960-2015.xls（见附表一）中的数据导入matlab。首先作出散点图，设定y(X59287)与x1(X54511)、x2(X57494)的关系为二元线性回归模型，即y=b0+b1x1+b2x2。之后作多元回归，求出系数b0=18.014，b1=0.051，b2=0.354，所以多元线性回归多项式为:Y=18.014+0.051*x1+0.354*x2。再作出残差分析图验证拟合效果，残差较小，说明回归多项式与源数据吻合得较好。若x1=30.4，x2=28.6时，y的数据缺失，则将x1，x2带入回归多项式，算出缺失值y=29.6888。类似地，若x1=40.6,x2=30.4时，y的数据缺失，则将x1，x2带入回归多项式，算出缺失值y=30.8462，即可补充缺失数据。 对于第二问，我们使用了异常值检验中标准差未知的t检验法。将除可疑测定值xd以外的其余测定值当做一个

1、打开EXCEL,找到XPS谱图的表。 2、用C1s的峰值减去284.6eV(等于 1.11)。 3、把需要分析的峰的数据拷贝到Origin里 面。 4、点击A(X)→点击鼠标右键 →点击Set Column Values →col(A)-1.11 → OK。

5、去除多余的数值(剩下453.04~ 462.09)。 6、作图。 7.点击Analysis→smoothing→FFT Filter→5→OK。

8、 双击左上角的1→移走data1_c→OK。 9、 双击smoothed把出现的数据拷贝到新 建的文本文档中并保存 。 10、打开XPS分峰软件。 11、点击data→import(ASCII)→新建 文本文档→打开。

12、点击background,选背景(一般为 linear),背景范围455~462。 13、add peak,选择Peak type:p。然后 添加合理的峰位、FWHM和L-G值。 14、调节峰位、峰面积、FWHM和L-G值。 使虚拟的峰位和实际峰位基本重合。

15、分别点击Save XPS、 Export(spectrum)和 Export to clipboard保存图。 16、打开画图板复制入图形可以得到想要 的数