资料处理与地质解释

资料处理与解释成果~

首先运用近场校正软件对原始数据进行了近场校正；然后通过空间滤波的手段，可以平滑地表局部异常源的影响，从而减小静态效应。另外，还利用了平移法，将有明显静态效应的整条曲线平移到正常值位置。在做了上述处理后，再利用圆滑反演软件进行正式反演。
CSAMT法接收机所采集到的原始数据，只是对应一系列频率的没有考虑装置特性的平均视电阻率值，不能准确地反映地下地层的真实电性特征；必须结合收发距r的大小，AB偶极的长度，发射电流的大小，采样间距等装置的详细数据进行合理的反演计算，才能得到对应地下不同深度处的电性分布剖面。在反演时采用了整条剖面的圆滑反演技术。
如图2所示，8线电阻率（ρ）断面图及地质解释。1740～1860低阻异常（＜2000 Ω·m），异常平面中心在1820点。地表岩体出露，经地表探槽揭露有含矿矽卡岩脉，推测为中寒武统杨柳岗组“残留体”与岩浆岩脉相互穿插而破碎引起的。1940～2060低阻异常与南北向花岗斑岩脉对应，异常没有圈闭。推测为寒武系地层穿插于岩体中形成破碎带。1280～1440高阻异常区（＞6000 Ω·m）是燕山期花岗闪长岩引起的，而它的东侧中阻（2000～6000 Ω·m）花岗闪长岩属于燕山期不同次侵入，都为成矿母岩。绿线为推测地质界线。440～1680 地表为杨柳岗组大理岩化灰岩、白云岩，地表岩体接触带附近分布花岗斑岩脉及含矿矽卡岩脉。在地层内部也分布有矽卡岩矿化脉。

图2 CSAMT法勘查8线综合剖面图

如图3所示，0线电阻率（ρ）断面图及地质解释。1760～1900低阻异常（＜2000 Ω·m），空间分布范围增大。异常中心在1840点，北侧地表探槽揭露含矽卡岩矿化脉。南北向花岗斑岩脉带经过此段。1420～1600高阻区（＞6000 Ω·m）是杨柳岗组大理岩化灰岩、白云岩引起，此地层出露于地表，形成独立的异常中心，异常中心在1500 点。在地表与岩体接触带部位及地层内部分布大量的含矿矽卡岩带、花岗斑岩脉。所以，该高阻异常是重点找矿区段之一。绿线为推测地质界线。1220～1420高阻异常（＞6000 Ω· m）是花岗闪长岩引起的，与成矿关系密切，与东侧中阻（2000～6000 Ω·m）花岗闪长岩电性差异明显。西侧黄柏岭组地层中低阻（＜6000 Ω·m）异常变化范围亦大。接触带部位也具有较好的成矿条件。
如图4所示，7线电阻率（ρ）断面图及地质解释。1740～1960低阻异常（＜2000 Ω·m）范围收缩，埋深变浅。异常中心在1860点。1160～1500高阻（＞6000 Ω·m）花岗闪长岩与东侧中阻（2000～6000 Ω·m）花岗闪长岩接触带陡立。西侧黄柏岭组地层中低阻（＜6000 Ω· m）接触带也较陡。1500～1620有一弱的低阻异常，可能为晚期花岗斑岩脉破碎带引起的相对低阻，也可能为寒武系与岩体破碎带。周围分布多处矿化点。综上所述，3条电阻率剖面图对各地质体的反映基本对应，曲线形态吻合，证明CSAMT方法的有效性。从整体上看，初步了解矿区内的地层、构造与岩体的空间分布规律。绿线为推测地质界线。

图3 CSAMT法勘查0线综合剖面图


图4 CSAMT法勘查7线综合剖面图

图5反映平面异常整体分布。西端中低阻（＜6000 Ω·m）为黄柏岭组地层；中段高阻带为岩体和杨柳岗组大理岩化灰岩。图5 表明，异常在此深度表现最全。1 号低阻异常较清晰，有较好的地质成矿背景，平面上也有较大空间；2 号高阻（＞6000 Ω·m）异常反映完整，2 号高阻异常是已知矿赋存主体，其下部以及与岩体接触带附近还有很大的成矿空间。4号异常，地表揭露含矽卡岩矿化脉及大理岩；5号为高阻岩体中的低阻异常，可能存在隐伏地质体。

图5 电阻率（Ω·m）平面等值线图

高密度电阻率法的测量系统在施工现场采集到大量关于地电断面结构特征的地质信息，并以数字的形式保存在随机存储器中。将其传入微机进行数据转换、处理与解释，然后生成供推断解释用的各类图件。图1.4.6为高密度电阻率法资料处理与解释系统框图，实际工作中可根据具体情况的需要选用其中某些过程。

图1.4.6 高密度电阻率法数据处理与解释流程框图

1.4.3.1 滤波处理方法
三电位电极系中，偶极和微分排列所测视电阻率曲线随着水平极距的加大，曲线由单峰变为双峰。绘成断面图时，除了和地质对象相对应的主异常外，一般还会出现强大的伴随异常。为消除或减弱三电位电极系视电阻率曲线中振荡成分的影响，简化异常形态，可以采用数字滤波方法，并将这种滤波称为扩展偏置滤波。
扩展偏置滤波器有4个非零的权系数：0.12，0.38，0.38，0.12。在滤波计算中，无论间隔系数为几的剖面测量结果，都应把滤波系数置于4个活动电极所对应的点上，在电极之间插入和电位相同的零系数。例如，n=2，滤波器长度为7，相应的权系数依次是：0.12，0，0.38，0，0.38，0，0.12。图1.4.7为二维地电模型正演模拟曲线的滤波处理结果。由图可见，未经滤波的剖面曲线随间隔系数的增大，曲线形态复杂；经滤波处理后的剖面曲线，形态大为简化，伴随异常的幅值减小并远离主极值。

图1.4.7 二维地电模型正演模拟曲线的滤波处理

点画线为滤波前原始曲线；实线为滤波后的结果
1.4.3.2 视比值参数
(1)Ts比值
高密度电阻率法的野外观测结果除了可以绘制相应装置的视电阻率断面图外，根据需要还可绘制两种比值参数图。考虑到三电位电极系中三种视电阻率异常的分布规律，选择了温纳β装置和温纳γ装置两种装置的测量结果为基础的一类比值参数。该比值参数的计算公式为

电法勘探

由于温纳β和温纳γ这两种装置在同一地电体上所获得的视电阻率总是具有相反的变化规律，因此用Ts参数绘制的比值断面图，在反映地电结构的分布形态方面，远比相应装置的视电阻率断面图清晰和明确的多。
图1.4.8 是对所谓地下石林模型的正演模拟结果。模型的电性分布如图1.4.8 c所示，其中温纳a装置的 拟断面图(图1.4.8a)几乎没有反映，而Ts比值断面图(图1.4.8b)则清楚地反映了上述模型的电性分布。

图1.4.8 模型上方视电阻率和比值参数Ts等值线断面图

(2)λs比值
另一类比值参数是利用联合三极装置的测量结果为基础组合而成的，其表达式为

电法勘探

式中：ρs(i)和ρs(i+1)分别表示剖面上相邻两点视电阻率值，计算结果示于i和i+1 点之间。
根据前面所学知识，有

电法勘探

则
lgλs(i，i+1)=lgFA(i)-lgFB(i+1)
而lgFA的差商为

电法勘探

令Δx=1，则lg λs为lg FA曲线的差商，或者说lg λs描述了岐离带曲线剖面水平方向的变化率。
图1.4.9表征比值参数λs在反映地电结构能力方面所作的模拟实验。视电阻率 断面图(图1.4.9a)只反映了基底的起伏变化，而λs比值断面图(图1.4.9b)却同时反映了基底起伏中的低阻构造。

图1.4.9 模型上方视电阻率和比值参数λs等值线断面图

1.4.3.3 统计处理
统计处理包括以下内容：
(1)利用滑动平均计算视电阻率的有效值，例如，三点平均
ρx(i)=[(ρs(i-1)+ρs(i)+ρs(i+1)]/3
式中：ρx(i)为i点的视电阻率有效值，i=1，2，3，…。
(2)计算统计参数
平均值： = · ρx(i)(N 为断面上的总测点数)
标准差：σA=∑
(3)计算电极调整系数

电法勘探

式中： (L)为电极距为L时全部视电阻率观测数据平均值。
(4)计算相对电阻率

电法勘探

通过计算相对电阻率，可以在一定程度上消除背景影响，相对视电阻率等值线断面图更能明显地反映地电体沿剖面的横向变化。
(5)对视参数分级
为了对视参数进行分级，首先必须按平均值和标准差关系确定视参数的分级间隔。间隔太小，等级过密，间隔太大，等级过稀，这都不利于反映地电体的分布。一般情况下，以采用五级制为宜，即根据平均值和标准差的关系划分四个界限

电法勘探

利用上述视参数的分级间隔，可将断面上各点的ρs(i)或ρy(i)划分成不同的等级，用不同的符号或灰阶表示时，便得到视参数异常灰度图。例如，ρs(i)＜D1，低阻；ρs(i)=D1～D2，较低阻；ρs(i)=D2～D3，中等；ρs(i)=D3～D4，较高阻；ρs(i)＞D4，高阻。
视参数的等级断面图在一定条件下能比较直观和形象反映地下介质的电性分布特征。当然，我们也可用不同的颜色来表示不同的等级，得到视参数异常色谱图。
统计处理原则上适应于三电位电极系中各种电极排列的测量结果，只是在考虑视电阻率参数图示时，由于偶极和微分两种排列的异常和地电体之间具有复杂得对应关系，因此一般只对温纳a装置的测量结果进行统计处理。温纳-斯伦贝尔格热装置的测量结果也可进行统计处理。
1.4.3.4 高密度电阻率法二维反演
在地球物理学中，地球物理反演是利用在地球表面观测到的物理数据推测地球内部介质物理状态的空间变化及物性结构。如果把地球物理问题分为资料采集、数据处理和反演解释三个阶段的话，那么，资料采集是基础，数据处理是手段，反演解释才是地球物理工作最终目的。
在高密度电阻率法中，仅根据高密度电阻率法的视参数等值线断面图(ρs、Ts、λs)或视参数分级灰度图来进行定性解释显然是很不够的。为了获得地下介质电性分布更为精确的电性结构与空间分布，目前最常用的方法是进行视电阻率断面二维反演。
(1)二维反演的基本原理
二维反演同一维反演的数学实质相同，都是要寻找一个地电模型，使其对应的理论计算值与实测视电阻率在一定法则下重合最好。数学上将求解这类问题的方法称为最优化方法。
用计算机对物探异常进行定量解释的最优化算法，实际上是求解多元函数极值的一种方法。最优化算法种类很多，如最速下降法(高斯法)、共轭梯度法、牛顿法、变尺度法(拟牛顿法)、最小二乘法、阻尼最小二乘法(马奎特法)、单纯性法等。
在电法资料定量解释中，最优化算法的基本步骤可归结为：
1)给出实测视电阻率离散值。
2)根据已知物性资料、地质资料和定性解释成果，确定地电模型，即给出地电模型初值。
3)通过正演计算得到地电断面的理论值。
4)评定理论计算值和实测视电阻率的拟合程度，即通过计算拟合差来判定拟合程度。
5)若拟合差不符合要求，则修改模型参数值，并根据修改后的地电模型参数重新计算理论值。
6)再次评定拟合程度，反复修改地电模型参数，直至拟合差达到事先给定精度为止。
7)此时地电模型参数即为解释结果。
完成上述过程的关键是用何种方法来评定拟合程度，以及如何修改模型参数。
评定拟合程度时，由于精确的地球模型我们并不知道，只能用方差来衡量它的质量，通常采用理论计算值和实测视电阻率两者偏差平方来作为衡量标准，即

电法勘探

式中：ε(x)为目标函数；δk(x)为偏差函数；Ek为实测数据；fk(x)为理论计算值；k为离散点之顺序；m为采样点数；x为理论模型参量。
式(1.4.4)的意思是，建立一个目标函数，使总误差恰好是向量δ(x)的欧几里得长度的平方。在数学上，范数经常作为长度或大小的某种度量，并且用一组双竖直线段来表示。如 即为向量e的范数。以向量元素的n次幂之和为基础的范数是最常用的，并且称之为Ln，如
L1范数：

电法勘探

L2范数：

电法勘探

随着幂次的逐步提高，e的最大元素的权也逐步增大。最小二乘法系用L2范数来定量表示长度。这里有一个问题，为什么采用L2范数而不是采用其他范数呢?这要看数据的特性。依据前述可知，幂次不同，对不同大小的误差所加的权也不同。换句话说，低阶范数对不同尺度给予的权相差不大，而高阶范数对大误差给予优惠。当数据服从高斯分布时，此时应用L2范数解是最合理的。然而，虽然高斯分布的数据在地球物理测量中是很常见的，但是大多数据集中总有那么几个“坏”数据，它们的存在说明数据并不完全服从高斯分布的统计性质。在反演之前不曾把这些坏数据去掉的话，那么最好用使L1范数最小化的反演方法，如线性规划法等，因为这时最小二乘法给予坏数据的权仍然过大。
如果给定初值x(0)，修正量为Δx，修正后的参量为x(1)，则有
x(1)=x(0)+Δx
参量修正后使目标函数ε(x)为最小值，即

电法勘探

此时拟合程度最佳。这就是最小二乘法意义下的最优化方法，即最小二乘法。最小二乘法实际上是一种反问题的长度解法，由于其简单易行，在电法资料解释中得到广泛应用。
如何求取修正参量Δx很重要。求取方法有多种，但不管那种算法都是通过不同途径修正Δx大小和方向，以达到拟合最佳的目的。
(2)二维断面反演的阻尼最小二乘法
设在被研究的地电断面Ω域上，待求的地电模型参数为M个(即M个待求的电阻率ρ1，ρ2，…，ρM)。通常，我们是将Ω域按一定法则剖分，例如，可按照图1.4.10进行剖分(变网格)，每个网格单元对应一个电阻率。这些单元中的模型参量，可用向量ρ=(ρ1，ρ2，…，ρM)表示。实测视电阻率拟断面上，取N个视电阻率采样值(即 ，i=1，2，…，N)。

图1.4.10 二维断面反演网格剖分示意图

二维拟断面反演是不断修改地电模型的电阻率参数，使理论模型拟断面向实测拟断面逼近。在理论计算值( ，i=1，2，…，N)向实测视电阻率 逼近过程中，通过不断改变电阻率值ρ，使理论计算的 与实测视电阻率 之间的误差尽可能小(一般＜5%)，以此作为衡量实测视电阻率 和理论计算 间拟合程度。通常采用对数型拟合方差F作为拟合视电阻率的目标函数(注意，以下的对数运算均系对电阻率的数值{ρi}Ω·m进行)，但简化记为

电法勘探

式(1.4.6)中 是根据初始模型参数正演计算的结果，它是地电断面参数ρ和电极距的函数，即 = (ρ，di)(di为与电极距有关的量)。目标函数F反映了实测拟断面数据与理论拟断面数据间的拟合程度，是模型参数的函数。二维视电阻率拟断面反演的目的就是要找到一组模型参数ρ=(ρ1，ρ2，…，ρM)，使目标函数取得最小值，即

电法勘探

由于理论计算的 (ρ，di)是模型函数ρ的非线性函数，故式(1.4.7)被称为非线性最小二乘问题，求取模型参数ρ拟合过程相当于数学上求多元函数极小值问题。对于非线性函数F直接求出ρ是很困难的，为此需要对非线性函数进行线性化近似处理。对假定的地电断面，给出一组模型参数初值ρ0=
( ， ，…， )，将ln (ρ，di)在初值ρ0附近做泰勒级数展开，将二阶和二阶以上的偏导数项略去，展开式的结果如下

电法勘探

令

电法勘探


电法勘探

则式(1.4.8)可写为

电法勘探

将式(1.4.11)代入式(1.4.6)，可得到目标函数F的近似表达式

电法勘探

将式(1.4.12)右端记为 ，则非线性最小二乘问题式(1.4.7)即可转换为线性最小二乘问题

电法勘探

根据极值存在的必要条件，使F达到最小的Δρj(j=1，2，…，M)，应满足下列方程组

电法勘探

整理后得

电法勘探

将式(1.4.15)写成矩阵形式
ATAΔρ=ATΔG (1.4.16)
式(1.4.16)即为目标函数(式(1.4.13))的法方程。其中：N×M矩阵A称为雅可比矩阵，其元素由式(1.4.9)来确定；ΔG为列向量，其元素为

电法勘探

求解方程式(1.4.16)，可得出模型参数的修正量Δρ，取ρ1=ρ0+Δρ作为新的模型近似值，若
F(ρ1)＜F(ρ0)
且F(ρ1)＜ε(ε为给定精度)，则ρ1作为二维地电断面的反演解释结果。若达不到精度，则以ρ1取代ρ0重复以上过程，直至求出符合精度要求的模型参数为止。
以上求解过程的特点是将非线性最小二乘问题(式(1.4.7))转化为求解一系列最小二乘问题——式(1.4.13)，虽然每一步求得的Δρ只是 =min的极小元，还不能使F(ρ1+Δρ)达到极小，但只要模型参数初值ρ0选取得当，这种逐步线性化的过程是收敛的。
法方程组的系数矩阵ATA，一般病态十分严重，甚至奇异。为保证反演过程收敛，增强法方程线的数值稳定性，可采用改进的阻尼最小二乘法(马奎特法)
(ATA+λS)Δρ=ATΔG (1.4.18)
式中：λ为阻尼因子；S为对角矩阵。
总结以上阻尼最小二乘法的反演思路我们可将其归纳为以下四个主要步骤：
1)给出初值；
2)计算理论拟断面(用2.5维数值模拟方法作正演计算)；
3)解法方程：(ATA+λS)Δρ=ATΔG；
4)让ρ1=ρ0+Δρ作为新模型参数，重新迭代反演。
可以看出，阻尼最小二乘法的主要计算工作量是解法方程，而要解法方程组，关键在于求出系数矩阵A(雅可比矩阵，Jacobian Matrix)，A中元素用式(1.4.9)来计算。可见雅可比矩阵的计算是反演成像计算中很重要的一步。
计算雅克比矩阵的方法有很多种。均匀介质模型是最为理想化的初始模型，利用这个理想化的模型，就可用解析的方法来计算偏导数值雅可比矩阵，这样就可以节约大量的机时，这里我们不再进行讨论。
1.4.3.5 高密度电阻率法的应用
(1)野外工作技术
A.测网布置
地球物理工作的测区一般是由地质任务确定的，总的原则与剖面法和测深法相同。对主要应用于工程及环境地质调查中的高密度电法而言，按工程地质任务所给出的测区往往是非常有限的，我们只能在需要解决工程问题的有限范围内布设测线、测网，可供选择的余地往往有限，这是一般工程物探经常遇到的情况。测网布设除了建立测区的坐标系统外，还包含了技术人员试图以多大的网度和怎样的工作模式去解决所遇到的工程地质问题。在这里，经验和技巧非常重要。特殊情况下，高密度电阻率法可布设不规则的测线和测网，尽可能在有限的测区内获得更多的测量数据。
B.装置选择
通常使用的装置还如上述四种类型(温纳、偶极-偶极、三极和温纳-斯伦贝尔格热装置)。不同厂家的测量系统基本上以这几种装置为主，但也各有特点。有的高密度电阻率仪提供了十多种装置以供选择。不同装置可联合使用，也可根据需要单独使用某一种装置。在高密度电阻率法中，合理地选择工作装置或其组合装置，可以提高采集数据对目标体的敏感度，放大异常，从而提高分辨率。选择一个合适的工作装置应考虑以下方面因素：探测目标的特性、探测深度、有效探测范围、信号强度、装置对地下电阻率水平或垂向变化分辨能力、场地噪声本底水平以及仪器灵敏度等。
当一个排列布好后，应尽量多地采集几个装置数据，以供后续处理和解释之用；同时，在一个地区工作时，除特殊情况外，应尽可能保持装置类型的一致，这样做的目的也是为了剖面的连接和解释的方便。
C.最小电极间距和排列长度的选择
最小电极间距和排列长度的选择取决于地质对象的大小和埋藏深度。要保证有足够的横向分辨率，探测目标体横向上至少要有2～3根电极通过。同时，由于高密度电阻率法实际上是一种二维测深剖面方法，所以在保证最大极距能够探测到主要地质对象的前提下，还要考虑围岩背景也能在二维断面图中得到充分的反映。如对小而深的探测目标体，要求较小的电极间距和较多的电极数。
对于长剖面，可以通过电极的移动来获得连续的断面数据。图1.4.11是温纳-斯伦贝尔格热装置通过两次移动来获得18a剖面长度的例子。一般地，在剖面对接时要重叠3～4个点，重叠点的数据取两次测量的平均值。

图1.4.11 温纳-斯伦贝谢装置移动测量示意图

(2)高密度电阻率法在工程与环境地质中的应用
近年来，高密度电阻率法在场地勘察、公路及铁路隧道勘查、坝基及桥墩选址、采空区和溶洞探测、地裂缝调查以及水库渗漏研究等领域得到广泛应用，取得了明显的地质效果和显著的经济效益。下面用几个高密度电阻率法的实际例子来说明该方法的应用效果。
A.在煤气管道勘查中的应用
场地地形有较小起伏，测线左侧是水田，右侧是砂石路。煤气管道包裹石棉网并涂有沥青，为高阻反映。实测时，最小电极间隔Δx=0.3 m，电极数为N=30，nmax=9。图1.4.12是经过地形改正后的ρs等值线断面图，图1.4.13是经过地形改正后并进行比值计算得到的视参数Ts等值线断面图，两图的高值等值线圈位置与管道位置相符，但Ts等值线断面图比ρs等值线等面图异常明显，这一结果与探地雷达探测结果一致，并与开挖后实际情况相符。

图1.4.12 煤气管道ρs等值线断面图

B.在水库坝体探测中的应用
中坊水库位于江西省广昌县境内，该水库建于20世纪50年代初期，坝体主要用素填土筑成。1989年初发现左坝肩背水面漏水，漏水位置位于土坝体与左坝肩接触面不远。当时发现渗水时渗水量很小，看不到明水。后来渗水量越来越大，到1998年渗水量达2 373.8 m3/d。如不采取措施，将影响坝体的安全。为了确定漏水通道及为续后灌浆加固处理提供资料，我们采用了高密度电阻率法和自然电场法进行了综合勘查。

图1.4.13 煤气管道Ts等值线断面图

以坝体与左坝肩接触面为中心，沿着坝体走向布置了一条高密度电阻率法剖面。最小电极间隔Δx=1 m，电极总数Psum=60。设第一根电极处x=0，接触面位置在x=29.5 m处。图1.4.14是经二维反演后得到的温纳a装置电阻率等值线断面图，图中清晰地反映了坝体与左坝肩接触面位置。

图1.4.14 横穿坝体与左坝肩接触面高密度剖面电阻率等值线断面图

通过本次工作(综合高密度电法、自然电场法以及水文地质调查结果)，基本查清水库漏水地段及漏水通道位置。观察结果表明，土坝体是完整的，从土坝与山体的接触面(接触面向坝体一侧倾斜)起一直延伸到山体内的范围内是渗水区。山体为花岗斑岩、风化程度比较高，裂隙发育，土坝体与山体的接触面及山体是漏水的通道，为后续打钻灌浆加固处理提供了依据。
近年来，在重大场地的工程地质调查、坝基及桥墩选址、采空区及地裂缝探测等方面，高密度电阻率法均取得较好的地质效果，实例很多。

浅层地震勘查资料处理和解释基本上都已实现了微机自动化。折射波和反射波法各有不同，就各自的流程框图作些简要介绍。

6.2.2.1 浅层折射波的资料处理与地质解释

由图6.2.2可见，首先对取得的折射波资料，进行观察判别，是否要进行预处理。如果干扰较强，要作滤波处理，先进行切除、均衡等抑制干扰措施。第二步是提取有效波的初至时间和绘制相应的时距曲线。第三步是进行地质解释。如果是水平地层或近似单斜地层，用常规方法，就可以取得好的效果。如果地层界面复杂，需采取相应的解释方法。常用的有t₀法和非时距曲线法等。此外还有哈里斯法、时间场法、射线追踪法等。

t₀法又称t₀差数时距曲线法，是常用于相遇时距曲线解释的方法。当折射界面的曲率半径大于深度很多时，效果最好。

根据射线的几何光学原理，如图6.2.3所示，有地下折射曲面R激发点O₁、O₂，取得相遇时距曲线S₁、S₂，在接收点D两波对应的走时分别为

环境地球物理学概论

在O₁、O₂点，两者走时相等（见图6.2.3），称互换时（T）。

环境地球物理学概论

如果界面R曲率半径不很大，则DBC可视为等腰三角形，由几何学原理可以导出深度

图6.2.2 折射波法资料处理流程图

环境地球物理学概论

若以t₀=t₁+t₂-T；K=v₁/2·cosi，则（6.2.3）式可改写为

环境地球物理学概论

根据斯奈尔定律，可得

环境地球物理学概论

式中v₁可根据直达波求得，而v₂一般根据数差时距曲线微分方程求得

环境地球物理学概论

图6.2.3 t₀法原理图

式中dt₁/dx和dt₂/dx分别为时距曲线（见图6.2.3）S₁和S₂的斜率。经过变换之后，可得

环境地球物理学概论

若界面倾角不大（小于15°时），可近似为

环境地球物理学概论

根据（6.2.6）式可求得dx/dθ（x），计算v₂；从而根据（6.2.5）求得K值，可以计算每个测点的深度h。

6.2.2.2 反射波法资料处理与解释

反射波的资料处理比折射波复杂，因为反射波采集的资料主要是多次覆盖得到的共激发点记录，伴随的干扰波较多。首先要对资料进行切除、校正和滤波处理，然后才能求取界面深度，进行地质解释。

其流程图，如图6.2.4所示。首先将采集的数据逐点输入计算机，并格式一致，切除干扰严重的纪录段，进行地形起伏和地表速度变化引起的时差校正。再将其反射点的纪录道从共激发点的地震纪录中抽出来，并按一定顺序构成新的共反射点道集（即CDP道集称抽道集）。不同类型的地震波具有不同的频谱成分，因此可以用频率滤波方法抑制干扰，这就是频谱分析和数值滤波。大多数还要进行速度谱分析，进行视速度滤波，在数值滤波中往往两者结合起来同时进行。通常称之为二维滤波。

图6.2.4 浅层地震反射波法资料处理流程图

对非水平地层，界面产状变化较大时，如倾斜、凹面、断层等，还要进行偏归位处理。这时按水平界面道集就不是真正的共反射点道集，偏移处理使反射波归位，反映真实界面形状。

流程（图6.2.4）中的纵坐标是以时间表示的时间剖面。为此，必须输入相应的速度参数，将时间剖面转即变为深度剖面需计算出反射界面的深度。确定反射面的构造形态和断层分布。

图6.2.5 古河道上高分辨率地震剖面

应用实例（探测古河道）。图6.2.25是一横切古河道的地震剖面。由图可见在530 ms处，反射波在剖面中部缺失而浅部和深部的反射波连续性都很好。这是古河道的特征，图中缺失部分即为古河道。

---

舒兰市13279489091： 地震勘探的勘探方法 - ？
司泪津博： 包括反射法、折射法和地震测井(见钻孔地球物理勘探).三种方法在陆地和海洋均可应用. 研究很浅或很深的界面、寻找特殊的高速地层时,折射法比反射法有效.但应用折射法必须满足下层波速大于上层波速的特定要求,故折射法的应用范...

舒兰市13279489091： 地震剖面的解释问题 - ？
司泪津博： 首先利用时间剖面进行解释 , 出时间域层面构造图 其次利用深度剖面进行解释 , 出深度域层面构造图 对照同一地质层面的时间域和深度域图件, 寻找速度陷阱. 由此判断是否使用深度域数据.在我国,大部分使用深度域的构造图,可以由深度偏移地震数据得来 ; 也可以利用时间偏移地震数据经过钻井标定后时深转换得来 . 如果处理得当 ,两种数据体的效果应该是一样的.

舒兰市13279489091： 地震资料解释员 是干什么的 - ？
司泪津博： 这项工作主要应用是能源勘探开发行业,也就是石油天然气和煤炭勘探开发行业,其他固体矿产也有少量人员从事.具体工作主要是用人工激发的地震波资料来解释地下的地质构造,主要是解释地层、圈闭和断层,现在基本都是在微机或工作站上人机交互解释,很少有在纸剖面上解释的了,所以最主要的问题是在电脑前呆久了的电磁幅射问题. 另外地震勘探是个很宽的行业,包括野外采集、地震资料数据处理、地震资料解释、地震资料储层预测和开发地震等方面,不知道你想问的是狭意的地震资料解释还是其它的. 从就业前景和待遇上来说,这个行业还不错,技术含量也可以,想出国留学也很容易,主要问题前面说了,一是在办公室坐着和电磁辐射,另一个是经常加班.

舒兰市13279489091： 地震资料解释员工作好吗 - ？
司泪津博： 比较累,成天坐在电脑前对着两个屏幕,做解释的最好再往地质方向发展一下,不然光在物探本专业的话从对自己解释的角度也不全面.但是看着图出来比较有成就感,当然看别人用你的图定井位的时候也会想自己应该学点综合地质.就收入方面做石油的都差不多,做解释相对而言自己从单位出来找工作或是干私活比较容易一些.

舒兰市13279489091： 在地震勘探观测系统的设计中应考虑哪些主要因素 - ？
司泪津博： 地震勘探的三大环节: 野外采集, 室内处理, 资料解释. (1) 野外采集: 按照预先设计的观测系统,炮点激发、检波器接收、仪器记录,得到原始地震资料(按时分道).数据通常记成SEGB或SEGD格式,班报有电子格式的和手写格式的.这一部分工作由物探地震小队完成. (2)室内处理: 将野外采集的原始地震资料转化为可用于地质解释的地震剖面 包括:预处理、常规处理和特殊处理三块内容. 这部分工作由资料处理中心完成. (3)资料解释: 结合地质、测井、录井、油藏工程等,进行综合解释.多由物探研究院、物探公司、地质研究院、采油厂地质所等完成.

舒兰市13279489091： 如何避免或降低测井资料解释应用的多解性 - ？
司泪津博： 地震资料的地质解释,指根据地震信息确定地质构造形态和空间位置,推测地层的岩性、厚度及层间接触关系,确定地层含油气的可能性,直接为钻探提供井位. 地震勘探的地质成效,在很大程度上取决于地震资料的正确与否.而要正确地解释...

舒兰市13279489091： 地球探测信息与技术这门学科是中国地质大学的更好还是中国石油大学的更好,适合女生吗 - ？
司泪津博： 石油大学主要是石油物探,以地震和测井为主,地质大学除了油气的地震和测井还有重磁电放,应用于固体矿产更多些,作为女生最好的话还是学地震解释或者测井解释,不需要出野外

舒兰市13279489091： 资源勘查和勘察技术与工程 - ？
司泪津博： 资源勘查和勘察技术与工程:本专业学生在学习数学、物理、化学、外语、计算机的基础上,主要学习普通地质学、机械设计与制造基础、工程力学、流体力学、电子技术、钻井工程、工程管理等方面的基本理论和基本知识,受到工程师的基本...

舒兰市13279489091： 磁测数据的处理与解释 - ？
司泪津博： (一)磁测数据的处理 在环境与工程测量中获得的磁测数据的处理与解释方法与矿产勘查中数据处理与解释方法基本相同.数据处理大体上可分为滤除干扰的一般处理和提取信息的专项处理两类.一般处理的目的在于滤除干扰,得到能客观反映...

舒兰市13279489091： 应用地球物理是一个怎样的专业? - ？
司泪津博： 业务培养目标:培养从事运用地球物理理论和方法研究地球内部结构,进行区域地质调查,金属与非金属矿产、油气资源勘查,水文地质与工程地质调查等方面工作的高级工程技术人才. 业务培养要求:本专业学生主要学习各类地球物理方法的...