勘查技术与工程（物探）数据处理实习报告
图 2-5 数据反演后等值线图
2.1.3 实测数据分析解释
由反演结果图可知，电阻率为1000欧姆左右，基岩组要为灰岩，其电阻率在2000Ω左右。
2.2 最小二乘法反演
2.2.1 基本原理
一维电测深最小二乘法反演 建立目标函数：
（10-4）
其中 为观测电测深视电阻率数据， 为计算机模拟的视电阻率数据，都为列向量。 用泰勒公式略去高次项近似为：
（10-5）
其中 为初始模型模拟值， 为偏导数矩阵， 为模型修改量。这样：
（10-6）
令 为 ，令目标函数取极小值，则有 。因此得线性方程组：
6
勘查技术与工程（物探）数据处理实习报告
（10-7）
解这个线性方程组得模型修改量。修改模型得新的模型。对新模型进行正演，得模拟观测数据。代入目标函数，看是否满足精度要求。如果不满足则继续反演，继续计算模型修改量直到满足精度要求为止，此时的模型即为反演结果。
2.2.2 数据处理
2.2.2.1我们进行直接反演法，我们首先对我们的数据进行预处理，对于两个测深为MN/2的数据平均，得到反演原始数据。
图 2-6 数据预处理
2.2.2.2 然后对每个测点的数据画出Grapher图，知道每个图的曲线类型，从而建立模型文件，为最小二乘法的运行做准备工作。
图 2-7 原始数据Grapher图
图 2-8 反演模型
2.2.2.3 利用好模型文件，对每个测点做最小二乘法反演，得到反演数据，反演模型。
7
勘查技术与工程（物探）数据处理实习报告
图 2-9 数据处理过程
2.2.2.4 对得到的反演数据进行Grapher画图，然后画图模型图，得到模型电测深曲线类型平面。
图 2-10 最小二乘法分层图
2.2.3 实测数据分析解释
根据最小二乘法可以看出，表面为素填土，反演数据分层在0.3-3m的位置为素填土，在130-140号点的位置素填土很薄，因为这段地层表面为卵石，所以出现了薄的地层，160号到170号点都是出现薄地层。地层往下是粘土，一般为20m,根据图可以
8
勘查技术与工程（物探）数据处理实习报告
看出，在初始值大概为4m。，可以看出，粘土层大概为2m。接下来为含角砾粘土层，地层由在9m-18m的位置，地层大概为10m。往下为基岩，石灰岩，25m。
2.3 对称四极测深
2.3.1 基本原理
电测深法又称电阻率垂直测深法。该方法是对一个测点用一系列由小到大的级距进行视电阻率测量，反映由浅入深的地层 垂向变化情况。通过对现场实测曲线进行分析和解释，可对观测点处垂向各地电性层的厚度和电阻率的大小。在实际工作中，一般采用对称四极电测深装置。
通过图形可以看出，电阻率测深曲线会出现脱节现象，主要是由于在一些MN/2不同的情况下，由于测量的深度不一致，导致测量浅部物质不均一造成的。
2.3.2 数据处理
2.3.2.1 通过对各点的数据进行处理，画图对称四极的等值线图。
图 2-11 对称四极等值线图
2.3.3 实测数据分析解释
通过图可以看出，在表层都出现的都是低阻图，从而可以看出都是一下填土，出现的电阻率一般都是很小的。从而看出地层的的地层分布不均一。
2.4 高密度电阻率法
2.4.1 基本原理
9
勘查技术与工程（物探）数据处理实习报告
这个程序使用的反演程序是建立在抑制平滑度最小平方法的基础之上的。建 立在拟牛顿最优化技术基础之上的最小平方法的一个新执行也能被使用。对于大 型数据组，这种技术要比常规的最小平方法快10倍多，而且要求的内存比较小。 抑制平滑度最小平方法是建立在下面的方程的基础上的。 （J`J+μF）D=J`g
优点就是：阻尼系数和平直度滤波器能被调整来与数据的不同类型相配。
2.4.2 数据处理
2.4.2.1 数据准备：将原始数据处理成瑞典反演软件要求的格式（如图2-12），删除VP和IP两列，将滚动数改为电极距，在数据头和数据尾加上必要的一些代码（注：第一行：测线号；第二行：最小的电极间距；第三行：排列类型（温纳=1，二级=2，偶极-偶极=3，联剖=6，施伦贝尔=7）；第四行：测量数据点总数；第五行：测量数据点的X位置的类型，表示X 位置是排列的第一个电极的位置，输入0；如果是排列的中点 (也就是：在拟断面中的测量数据点的位置），最后将数据存为txt文件并改名为dat类型
图 2-12 数据预处理部分 图 2-13软件导入类型文件图 2.4.2.2 反演过程：打开瑞典高密度处理软件RES2DINV.EXE，点击“File”→“Read data file”,将准备好的文件读入程序，弹出对话框（如图2-14）点击“确定”按钮；完成数据读入后可点击“Edit”→“Exterminate bad datum points”删除异常点（如图2-15），完成后重新保存数据，并再次读入；然后点击“Inversion”→“Least-squares inversion”，弹出对话框，文件名填入“r.INV”文件类型选择“.inv”格式，然后点击
10