地震勘探资料数字处理

中国地质大学（北京） 课程名称：应用地震学 教师：段云卿 第25册

第四章：地震勘探资料数字处理

野外采集到的原始资料是以二进制的数字形式记录在磁带上，必须经过计算机的各种运算，才能输出供地震地质解释的各种资料，或直接输出某些解释成果，本章介绍如何进行数据处理。

§4.1校正和叠加处理

一、动校正

1．动校正的含义：（§3.5）

（1） 对于一次覆盖共炮点资料来说，把双曲线型或近似双曲线型反射波同

相轴拉直，也就是消去炮检距不为0对反射波旅行时的影响，使同相轴能直观地反映地下界面的构造形态。

（2） 对于共反射点道集来说，把各道均校正成共中心点M 处的自激自收道，

再叠加起来作为共中心点M 处的叠加道，使一次波同相叠加而加强，多次波等干扰波非同相叠加而减弱。

2．动校正公式（§3.5）

2

022V

t x

t =

∆ (6.2-26)

3．计算动校正量（使用共反射点道集） (1)公式

为了对共反射道集的每一道的整个道进行计算，将（6.2—26）改写为：

2

002

)

(2i i j

ij t V t x t =

∆ （j=1,2,……,n ； i=1,2,……,m ） (6.4-1)

j —— 道序号。 i —— 采样点序号。

x j —— 第j 道的炮检距。 n —— 覆盖次数。 M ——道长 t 0i ——为第i 个界面共中心点处自激自收时间。 (2)问题

不知什么地方有反射界面，就不知什么地方有反射波。 不知反射波的t 0时间。

中国地质大学（北京） 课程名称：应用地震学 教师：段云卿 第25册

（3）解决方法

地震道上有一个采样值就有一个反射波。

地震道上每一个采样点的时间i △，都看成一个t 0时间，记为t oi 。 （4）例子

①设采样间隔△=4ms

②长为0.5S -4.5S

的记录,就有1001个t 0值：

)(5.00,0s t = )(004.05.01,0s t += )(004.025.02,0s t ⨯+= )(004.05.0,0s i t i +=

)(004.010005.01000,0s t ⨯+=

③对任意一道就有1001个动校正量。例如炮检距为1000m 的第j 道，动校正量为：

)(207.0)

5.0(5.021000

2

2,0s V t j =⨯⨯=

∆

)

(205.0)

504.0()504.0(21000

2

2,1s V t j =⨯⨯=

∆

)

(204.0)

508.0()508.0(21000

2

2,2s V t j =⨯⨯=

∆

)

?()

004.05.0()004.05.0(21000

2

2,s i V i t j i =+⨯+⨯=

∆

)(000.0)

5.4()5.4(21000

2

2,1000s V t j =⨯⨯=

∆

st

同一个地震道上，每个采样点的动校正量不同，所以叫动校正。 同一个地震道上，x 是固定不变的，t 0↗,V ↗,Δt ↘。 ④手工动校正

在第j 道上，将ij i ij t t t ∆+=0处的振幅值搬到t 0i 处 上例：（0.5+0.207=0.707）s 处的振幅值搬到0.5s 处 （0.504+0.205=0.709）s 处的振幅值搬到0.504处

问题：手工动校正麻烦。如何计算机化？ 4．动校正的计算机化 ①∆∆/ij t

上例： j x (m) t 00/∆ ij t ∆(ms) ∆∆/ij t 1000 125 207 52 205 51 204 51 204 51 201 50 0 0

注意：有许多点有相同的动校正量

②简化

x j——第j道的炮检距。

t00/∆——起始时间是采样间隔的几倍。

t max/∆——最大动校正量是采样间隔的几倍。

k1——具有最大动校正量的个数。

k2——具有最大动校正量减1的个数。

k3——具有最大动校正量减2的个数。…………

③动校正量板表

1000，125，52，1，3，……156，32，1，3

5．动校正在计算机上的实现

①造动校正量板表例如：20，10，4，1，1，2，3，4

②动校正——逐点搬家法

10 10+4(t0+△t→t0: 10+4→10)

例如:

提醒：实际中动校正前、后用同一个一维数组的内存单元。

③程序要点

一道原始记录放在一维数组x中备注

N是道记录的长度一道原始记录放在一维数组x中从动校正量板中取来：k1 ,k2 ,k3 ,…动校正后仍然放在一维数组x中

m=t00/Δ, n=Δt max/Δ

DO I=1,m-1 且I≤N, m≥1 x（1）=0

x(I)=0 x(m-1)=0

DO I=m+n, k1+m+n 且I≤N x(m)=x(m+n)

x(I-n)=x(I) x(k1+m)=x(k1+m+n) 插值一个空格 x(k1+m+1)=插值

DO I=k1+m+n+1, k2+k1+m+n+1 且I≤N x(k1+m+2)=x(k1+m+n+1)

x(I-(n-1))=x(I) x(k2+k1+m+2)=x(k2+k1+m+n+1) 插值一个空格 x(k2+k1+m+3)=插值

DO I=k2+k1+m+n+2,k3+k2+k1+m+n+2且I≤N x(k2+k1+m+4)=x(k2+k1+m+n+2)

x(I-(n-2))=x(I) x(k3+k2+k1+m+4)=x(k3+k2+k1+m+n+2) 插值一个空格 x(k3+k2+k1+m+5)=插值

I=N 直到I=N时，停止

4．动校正对波形的畸变及处理方法

①动校正搬家后出现“空格”，波形被拉长，波形发生了畸变，浅层畸变大，

深层畸变小。

②对深层“空格”插值，一般用线性插值。对浅层畸变历害的进行“切除”。

二、静校正

（一）静校正的概念

1．低降速带

①低速带