25
receiver谱 Receiver 谱= receiver谱 Cmp 谱= 4590 Continue 4595 Continue
保留本点叠代后谱，用于下次迭代 保留本点叠代后谱， 4600 continue 4900 continue 5000 continue
东方公司物探技术研究中心 保留最终炮点，接收点谱，CMP，偏移距域的谱
原始谱与地表一致性反褶积
15
原始
用4项（炮，检，trace,CMP)求解算子
东方公司物探技术研究中心
地表一致性反褶积，用3项与4项差不多
16
用3项（炮，检，trace）求解算子
用4项（炮，检，trace,CMP)求解算子
东方公司物探技术研究中心
地表一致性反褶积（一步法与三步法）对比
17
没有起始频率，但做了 起始频率10 高通滤波，（8，12）
12
1.FFT 子程序不统一。 2.GOS_ABORT()调用后没有return，多处。 3.输入道在反褶积之前平均振幅级别规格化到5000. （有一个理由可以解释：输入数据没做振幅补偿或 做振幅补偿没有规格化到某个级别） 4.程序中考虑炮点项，检波点项，CMP项，偏移距项， 程序中偏移距项用道号代替，在2D可能还可以， 3D是不是有误差？ 5.输出之前做了高通滤波H(8,12），也许是损失低频的原因之一。 6.三步法中有一个起始，终止频率，起始缺省为10，低频保留要比 加高通好。
东方公司物探技术研究中心
3D数据，高频噪声谱与正常道谱
5
存在大量野值，不去 野值，做地表一致性 反褶积，不论单道还 是多道，其谱似乎损 失低频，其实是野值 造成。此谱是多道统 计，有效信号受压。
东方公司物探技术研究中心
3D数据，去野值后单炮及频谱
6
原始数据及谱
东方公司物探技术研究中心
一步法与三步法地表一致性反褶积对比 （3D数据参数相同）
7
内部加高通 滤波 反褶积后依 然存在
一步法
三步法
东方公司物探技术研究中心
一步法与三步法地表一致性反褶积谱对比 （3D数据参数相同）
8
低频面 少低频？ 波
高频台的 高
一步法
三步法
东方公司物探技术研究中心
SCDC前后的叠加剖面及频谱对比 前后的叠加剖面及频谱对比 Omega
9
原始
似乎少低 频？
PD02
东方公司物探技术研究中心
三步法与一步法地表一致性反褶积对比 （2D数据参数相同）
2
三步法
一步法
东方公司物探技术研究中心
3D数据高频噪声谱与正常道谱
3
东方公司物探技术研究中心
SCDC前后的单炮及频谱对比 前后的单炮及频谱对比(1) 前后的单炮及频谱对比 Omega
4
原始
PD02
高频台的 相对高
GeoEast
19
三步法（有起始频率10）
不加高通滤波，保留低频，但是出现低频噪声
东方公司物探技术研究中心
原始与三步法
20
东方公司物探技术研究中研究中心
一步法（加滤波（8，12），与不加滤波）
22
东方公司物探技术研究中心
一步法（加滤波H（1，6）），谱要好一些
23
NO filter
三步法（S-R）
一步法（S-R-T-CMP)
如果说损失低频，怀疑加高通滤波造成。
东方公司物探技术研究中心
地表一致性反褶积（一步法与三步法）对比
18
没有起始频率，但做了 高通滤波，（8，12）
起始频率10
三步法（S-R）
一步法（S-R)
三步法两项与一步法两项
东方公司物探技术研究中心
加高通（8，12）太高，可加（1，5）适当压制一下
东方公司物探技术研究中心
三步法地表一致性反褶积存在问题
13
1.低频面波有放大可能性，起始频率用6HZ与用12HZ面波仍有。 2.分解时可以四项，但应用时只用两项。 3.自相关一致性不好与低频有关。
东方公司物探技术研究中心
原始谱与地表一致性反褶积
14
原始
只用两项（炮，检）求解算子
东方公司物探技术研究中心
H(1,6）
东方公司物探技术研究中心
某道分解的各域的谱
24
某道下边五项之和的 log谱 某道对应的各域log谱 rrr ：receiver sss：shot ooo：offset ccc: cmp 全测线平均log谱
东方公司物探技术研究中心
Do 5000 kf=1, 迭代次数 Do 4900 kw =1,4 (四个域 offset,shot,receiver,cmp) (四个域 同一点的最大道数：如对于shot 1，最大道数420 Do 4600 I=1,param(60+ku) 同一点的最大道数：如对于shot 1，最大道数420 J=1,rmax从此范围内找同一道： 420道 10炮 Do 4595 J=1,rmax-rmin+1 （从此范围内找同一道：1炮420道，共10炮， 找第1 找第1道，需从1-3781道顺序内即可 需从1 3781道顺序内即可 找第2 找第2道，需从2-3781 道顺序内即可 需从2 。。。。。。。。。。 其它域意思一样。 其它域意思一样。 Do 4590 ii =1,4 （1） 找到本道在其它域中的谱 （2） 本道的求和谱 – 本道在其它域中的谱 迭代： 迭代： Offset Shot receiver谱 shot谱 Cmp谱 谱= offset 谱 –receiver谱 - shot谱 - Cmp谱 receiver 谱= receiver谱 offset谱 Cmp谱 shot 谱 –receiver谱 - offset谱 - Cmp谱 receiver shot谱 offset谱 Cmp谱 shot谱 - offset谱 - Cmp谱 shot谱 offset offset谱 receiver谱 shot谱 –offset谱 - receiver谱 cmp谱 cmp谱 -
一步法与三步法地表一致性反褶积对比
1
本测试通过与一步法地表一致性反褶积对比，试图回答三个问题： （如下问题出现在三步法地表一致性反褶积中，属于常用模块） 1. 地表一致性反褶积：炮域统计自相关结果横向稳定性较差 。 2. 地表一致性反褶积处理当预测步长较小时， 出现低频信号损失的现象，在叠加剖面上尤为明显 。 3. 高频信号的能量出现极不均匀情况， 随着预测步长的增大，现象逐渐消失。 没有解释清楚为啥三步法地表一致性反褶积后频谱抖动大。 本测试可供参与此模块完善时用。
GeoEast
东方公司物探技术研究中心
一步法与三步法地表一致性反褶积谱对比 （3D数据参数相同）
10
原始
高频台的 高 三步法
一步法
东方公司物探技术研究中心
一步法与三步法地表一致性反褶积自相关对比 （3D数据参数相同）
11
一步法整道自相关 没有对 齐
三步法整道自相关
东方公司物探技术研究中心
一步法地表一致性反褶积存在问题