剩余静校正
地震数据处理第五章：静校正
总的低速带校正量为：
' j
hl )
静校正前
地面 V0
低速带底面 V
反射界面
第一步：井深校正后
V0 V
地面 低速带底面 反射界面
第二步：地形校正后
V0 V
基准面 低速带底面 反射界面
第三步：低速带校正后
基准面
反射界面
小结
1、符号约定：剥去地层时间为负，即减去静校正 量为负号；填充地层时间为正，即加上静校正量为 正号。 2、最终基准面校正量计算公式为
近地表沉积的介质相对深层而言，沉积年代
相对较短，长年的风化作用使近地表沉积的介质 疏松，无胶结或半胶结，地层中含水与不含水， 含水量的多少都会引起地球物理特征的变化。
近地表厚度和速度的各向异性、地表高程起
伏都会对地震波场造成不等量的延迟，延迟的大 小与近地表地层的物性有关，这种延迟时若不校 正，将会影响到叠加成像和构造形态的可靠性。
ESW—炮点处风化层高 程 EGW—检波点处风化层 高程 ESR—炮点处参考基准 面高程
EGR—检波点处参考基 准面高程
Hs —— 井深
炮点基准面静校正量为：
TS
TSW
- TSR
ES
- E SW vW
hs
ESW ESR VR
检波点基准面静校正量为：
TG
TGW
- TGR
EG - EGW v
（1）基准面校正；CMP叠加参考面校正；低降速 带底面校正；
（2）控制点数据线性内插法（微测井、小折射方 法等建立控制点数据）；
（3）沙丘曲线法（根据沙丘厚度在延迟时曲线上 找到对应的延迟时，计算静校正量）；
（4）相似系数法； （5）数据库法（建立导线成果、浮动基准面高程、 地表高程、小折射成果、高速层顶深度、潜水面深度 等数据库）。
第八节 地表一致性反射波剩余静校正
第八节 地表一致性反射波剩余静校正许多静校正方法均采用地表一致性模型。
检波器组在位置i 的延迟G i 和震源在位置j 的延迟S j ，对所有相应的地震道都是相同的。
如果i 和j 具有公共的坐标原点，其炮检距正比于（j －i ）。
如果沿测线有构造，我们以CMP （共中心点）号位置k 来表示延迟量L k ，其中k ＝（j - i ）/2，表示构造的深度比其它位置深L k 个单位，L k 是构造时移的均值。
对于平界面反射而言，它即指向中点位置。
如果倾角很缓，L k 对于共中心点来说几乎是常数。
如果动校正速度有误差，就会保留一些剩余时差M k ，它随炮检距的平方而变化。
如果不考虑炮点或检波点对测线的横向偏离，那么对于地表一致性模型，一个道总的时移量t i,j 为： t i,j =G i 十 S j ＋ L k 十 Mk （j – I ）2 (8.1)地表一致性模型不只限于确定静校正量的时移，在振幅调节、子波提取、反褶积及其它算法中，有时候都是基于地表一致性模型，均按上述相同的过程进行。
我们并不知道每一道的时移量，但可以利用互相关求一个道相对于另一个道的时移量（t i,j – t m,n ）：t i,j – t m,n =G i –G m +S j –S n +L i+j – L m+n +M i+j （j - i ）2 – M m+n （n –m ）2 （8.2）在此没有使用下标k ，而使用i+j ，是确保下标是整数，具体数值无关紧要。
对于一个CMP 道集，两两道互相关，就可以得出比未知数（G i ，S j ，L i+j ，M i+j ）个数还要多的方程式，是一个超定方程组。
但是，方程(8.2)在测量过程中也包含一些不确定的因素，例如等号左边就会出现误差。
一个“超定”，一个“不确定”，就使我们有办法来求解这个方程组，通常用最小平方法，有时也用迭代法。
最小平方问题是使误差e p 的平方和最小：()[]{}最小值∑∑=+-+-+---==++++2n m 2j i n m ,,2m) - (n M - i -j M L )(j i n j m i n m j i p L S S G G t t e E (8.3) 可用下式求解：,0,0,0,0=∂∂=∂∂=∂∂=∂∂++j i j i j i M E L E S E G E (8.4)通常，我们事先构成一个横型道，例如经过一般处理后的本道集的初步叠加道，或者是经过时移处理后的前一个道集的叠加道。
折射波剩余静校正方法
*北京市海淀区学院路中国地质大学地下信息探测技术与仪器教育部重点实验室,100083本文于2005年2月21日收到,修改稿于同年10月2日收到。
・处理方法・折射波剩余静校正方法段云卿*(中国地质大学地下信息探测技术与仪器教育部重点实验室)摘 要段云卿.折射波剩余静校正方法.石油地球物理勘探,2006,41(1):32～35山地、沙漠及其他复杂地表地区地震资料的线性散射噪声和随机噪声很强,有效反射信号弱,资料信噪比较低,静校正问题严重,使用常规剩余静校正方法难以见效。
本文利用折射波信噪比高的特点,将反射波剩余静校正方法应用于折射波资料处理,通过交互手段,逐段估算折射波的速度,用合适的速度对地震记录进行线性动校正,在共炮点或共中心点道集上,用相关方法计算各道与模型道时差,再用统计方法计算出炮点和检波点剩余静校正量。
将该方法应用于信噪比较低、反射波剩余静校正方法难以奏效的复杂地表区,获得良好处理效果。
关键词 剩余静校正 折射波法 共炮点道集 共中心点道集 复杂地表区 模型道1 引言静校正是地震资料处理中至关重要的一环。
我国西部地区地表条件极为复杂,静校正问题尤为严重。
如在沙漠、戈壁、黄土塬或山地等复杂地表区,地形起伏大,表层岩性变化非常剧烈,低降速带厚度变化大,激发和接收条件复杂,近地表条件纵横向千差万别。
近地表地形和低降速带的影响导致地震反射资料不能准确成像,也造成地下构造发生扭曲。
因此,研究复杂地表区静校正方法,对于提高地震勘探精度、降低勘探风险及节约勘探成本有着重要的意义。
本文基于反射波剩余静校正思路,提出一种实现折射波剩余静校正的方法,从而较好地解决了信噪比较低、反射波剩余静校正方法难以奏效地区的静校正问题。
2 方法实现思路静校正的常规步骤为:首先对地震资料进行野外静校正;随后进行折射波静校正;在动效正之后,再进行反射波剩余静校正。
通过这些处理,可初步解决长、中、短波长静校正问题。
但在山地、沙漠及其他复杂地表区,由于线性散射噪声和随机噪声强,有效反射信号弱,地震资料的信噪比往往较低,因此采用常规剩余静校正方法不能建立准确的模型道而达不到预期处理效果。
剩余静校正
（三）剩余静校正1．定义及分类（1）剩余静校正的定义由于技术上的原因或某些人为因素，例如低速带速度及厚度难以测准，使得野外实测资料往往不很准确，故进行了野外静校正后仍残存着剩余静校正量。
提取剩余静校正量并加以校正的过程叫剩余静校正．．．．．。
因为采用自动统计方法求取剩余静校正量，故也叫自动统计静校正．．．．．．．。
（2）剩余静校正的分类剩余静校正量‖‖短波长剩余静校正量高频剩余静校正量 ++长波长剩余静校正量低频剩余静校正量P125 图6.4-7①长波长或低频剩余静校正量长波长剩余静校正量是区域性的大范围内的异常，地形、低速带变化缓慢，自动统计方法无效。
②短波长或高频剩余静校正量短波长剩余静校正量是局部性的小范围内的异常，地形、低速带变化剧烈,自动统计方法有效。
2．自动统计法求取短波长剩余静校正量的假设和特点自动统计剩余静校正利用多次覆盖资料，所以假设和特点与多次覆盖有针对性。
(1)两点基本假设假设①：波在低速带内垂直于地面传播，即同一炮点或同一接收点的剩余静校正量相同。
假设②：在一个排列的长度上，各炮点（或接收点）的剩余静校正量是随机的，其均值为O 。
(2)三个特点特点①：某个记录道的相对剩余静校正量△τ包括炮点剩余静校正量△τO 和接收点剩余静校正量△τg ，即g O τττ∆+∆=∆特点②：在同一个共炮点道集中，因为各个道的炮点剩余静校正量相同，接收点剩余静校正量不同。
所以将各道的相对剩余静校正量求和再平均，就得到该炮点的剩余静校正量。
地面 V 0低速带底面 V反射界面 例如:一个共炮点道集有24个接收道：第O i 炮的第 1道的剩余静校正量：11g Oi τττ∆+∆=∆ 第O i 炮的第 2道的剩余静校正量：22g Oi τττ∆+∆=∆第O i 炮的第24道的剩余静校正量：2424g Oi τττ∆+∆=∆ 对以上24个式子相加再取平均值有：∑∑∑===∆+∆=∆241241241241241241j j gj Oi j j τττ τ∆ Oi τ∆ 0令左 边 项：∑=∆=∆241241j j ττ右边第一项：∑=∆=∆241241j Oi Oi ττ右边第二项：0241241=∆∑=j gj τ （根据假设②）上式成为：Oi ττ∆=∆上式物理意义：在同一个共炮点道集中，将各道的相对剩余静校正量求和再平均，就得到该炮点的剩余静校正量。
共炮(检)点剩余静校正方法
返回目录
４ ８
石 油 地 球 物 理 勘 探
２ ０ １ １年
［ ］ ８ 地震记录的总静校正量狋 ， 即 犻 犽 都可分解为 ４ 项 犼 ２ （ ） 狋 狉 犮 犺 １ 犻 犽 ＝狊 犻＋ 犽 犻 犽 犼 犼＋ 犼 ＋狔
根据式 （ ） 可知 ， 共 炮 点 道 集 经 过 叠 加 后， 生成的模 ４ 型道仅仅存在炮点静校正量的影响 。 使用 模 型 道 与 叠 前 道 进 行 相 关， 就可以计算出 对应检波点的剩余静校正量 。 在共炮点道集中求解 检波点剩余静校正量的过程可以作为一个最优化问 题， 即
： ６ １ ０ ５ ０ ０。Ｅ ｍ ａ ｉ ｌ ｌ ａ ｒ ｅ ９ ９ ９＠ｓ ｉ ｎ ａ ． ｃ ｏ ｍ 四川省成都市新都区西南石油大学资源与环境学院 ， ｇ 本文于 ２ 最终修改稿于 ２ ０ ０ ９年１ １ 月 ２ 日收到 ， ０ １ ０年１ １月２ ５ 日收到 。 基金项目 ： 教育部博士点基金项目 （ ） 及地球探测与信息技术四川省重点学科建设资助 。 ２ ０ ０ ８ ０ ６ １ ６ ０ ０ １ ５
２ 狋 狉 犮 犺 犻 犽 ＝狊 犻＋ 犽 犻 犽 烄 犼 犼＋ 犼 ＋狔 （ ） ２ 烅 ２ ， ， 狋 狉 犮 犺 犻 １， 犽 １ ＝狊 犻＋ １＋ 犽 １ 犻 １ ＋狔 犽 １ ＋ ＋ ＋ ＋ ＋ ＋ 烆 犼 犼 犼 从式 （ ） 可以看出 ， 在同 一 炮 中狊 如果将共炮 ２ 犻 相 同，
由于共炮点道集共检波点道集中存在相同的炮点检波点静校正量其静校正量的变化相对于cmp道集的静校正量变化小据此进行水平叠加可生成信噪比更高的模型道从而为解决低信噪比大剩余静校正量地震资料的静校正问题提供前提条件
０ １ １年２月 ２
第４ ６卷 第１期
山地地区地表一致性剩余静校正技术应用技巧
３ １ ２・
价 值 工 程
山地地 区地表一致性剩余静校正技术应用技巧
Ａｐ ｐ ｌ ｉ ｃ ａ ｔ ｉ ｏ ｎ Ｓ ｋ ｉ ｌ ｌ ｏ ｆ Ｓ ｕ ｒ ｆ ａ ｃ ｅ Ｃｏ ｎ ｓ ｉ ｓ ｔ ｅ ｎ ｔ Ｒｅ ｓ ｉ ｄ ｕ ａ ｌ Ｓ ｔ ａ ｔ ｉ ｃ Ｃｏ ｒ ｒ ｅ ｃ ｔ ｉ ｏ ｎ Ｔｅ ｃ ｈｎ ｏ ｌ ｏ ｇ ｙ ｉ ｎ Ｍｏ ｕ ｎ ｔ ａ ｉ ｎ ｏ ｕ ｓ Ａｒ ｅ ａ
Ａｂ ｓ ｔ ｒ ａ ｃ ｔ ：Ｄｕ ｅ ｔ ｏ ｃ ｏ ｍｐ ｌ ｅ ｘ ｇ ｅ ｏ ｌ ｏ ｇ ｉ ｃ ａ ｌ ｃ ｏ ｎ ｄ ｉ ｔ ｉ ｏ ｎ ｓ ｏ ｆ ｍｏ ｕ ｎ ｔ ａ ｉ ｎ ｏ ｕ ｓ ａ ｒ ｅ ａ ａ ｎ ｄ ｐ ｏ ｏ ｒ ｒ ｅ ｃ ｅ ｉ ｖ ｉ ｎ ｇ ｃ ｏ ｎ ｄ ｉ ｔ ｉ ｏ ｎ ， ｔ ｈ ｅ ｏ ｒ ｉ ｇ ｉ ｎ ａ ｌ ｍａ ｔ ｅ ｒ ｉ ａ ｌ Ｓ ＮＲ ｉ ｓ ｌ ｏ ｗ． Ｉ ｔ ｉ ｓ ｄ ｉ ｉｃ ｆ ｕ ｌ ｔ ｔ ｏ ａ ｃ ｃ ｕ ｒ ａ ｔ ｅ ｌ ｙ ｉ ｄ ｅ ｎ ｔ ｉ ｆ ｙ ｔ ｈ ｅ ｅ ｆ ｆ ｅ ｃ ｔ ｉ ｖ ｅ ｖ ｅ ｌ ｏ ｃ ｉ ｔ ｙ ｂ ｅ ｆ ｏ ｒ ｅ ｒ ｅ ｓ ｉ ｄ ｕ ａ ｌ ｓ ｔ ａ ｔ ｉ ｃ ｃ ｏ ｒ ｒ ｅ ｃ ｔ ｉ ｏ ｎ ．Ｔ ｈ ｅ ｅ ａ ｓ ｉ ｌ ｙ ｉ ｄ ｅ ｎ ｔ ｉ ｉ ｆ ｅ ｄ ｖ ｅ ｌ ｏ ｃ ｉ ｔ ｙ ｃ ａ ｎ ｂ ｅ ｇ ａ ｉ ｎ ｅ ｄ ａ ｆ ｔ ｅ ｒ ｓ ｅ ｖ ｅ ｒ ａ ｌ ｔ ｉ ｍｅ ｓ ｏ ｆ ｖ ｅ ｌ ｏ ｃ ｉ ｔ ｙ ａ ｎ ａ ｌ ｙ ｓ ｉ ｓ ａ ｎ ｄ ｒ ｅ ｓ ｉ ｄ ｕ ａ ｌ ｓ ｔ ａ ｔ ｉ ｃ ｃ ｏ ｒ ｅ ｃ ｔ ｉ ｏ ｎ ． Ｂ ｕ ｔ ｂ ｅ ｃ ａ ｕ ｓ ｅ ｏ ｆ ｔ ｈ ｅ ｅ ｒ ｒ ｏ ｒ ｏ ｆ ｉ ｎ ｉ ｔ ｉ ａ ｌ ｖ ｅ ｌ ｏ ｃ ｉ ｔ ｙ ａ ｎ ｄ ｉ ｎ ｉ ｔ ｉ ａ ｌ ｓ ｔ ａ ｔ ｉ ｃ ｓ ， ｔ ｈ ｅ ｉ ｔ ｅ ｒ ａ ｔ ｉ ｖ ｅ ｏ ｆ ｒ ｅ ｓ ｉ ｄ ｕ ｌ ａ ｓ ｔ ａ ｔ ｉ ｃ ｃ ｏ ｒ ｒ ｅ ｃ ｔ ｉ ｏ ｎ ｉ ｓ ａ ｌ ｗａ ｙ ｓ ｎ ｏ ｔ ｃ ｏ ｎ ｖ ｅ ｒ ｇ ｅ ｎ ｔ Ｓ Ｏ ｔ ｈ ｅ ｒ ｅ ｌ ａ ｕ ｎ ｄ ｅ ｒ ｇ ｒ ｏ ｕｎ ｄ ｓ ｔ ｒ ｕ ｃ ｔ ｕ ｒ ｌ ａ ｃ ｏ ｎ ｉ ｆ ｇ ｕ ｒ ａ ｔ ｉ ｏ ｎ ｃ ａ ｎ ｎ ｏ ｔ ｂ ｅ ｇ ａ ｉ ｎ ｅ ｄ ．Ｏｎ ｌ ｙ ｒ ｅ ｌ ａ ｔ ｉ ｖ ｅ ｌ ｙ ｒ ｅ ａ ｓ ｏ ｎ ａ ｂ ｌ ｅ ｉ ｎ ｉ ｔ ｉ ａ ｌ ｖ ｅ ｌ ｏ ｃ ｉ ｔ ｙ ｃ ａ ｌ ｃ ｕ ｌ ａ ｔ ｉ ｎ ｇ ｍｅ ｔ ｈ ｏ ｄ ｃ ａ ｎ ｅ ｎ ｓ ｕ ｒ ｅ ｔ ｈ ｅ ｒ ｅ ｌ ａ ａ ｎ ｄ ｒ ｅ ａ ｓ ｏ ｎ ａ ｂ ｌ ｅ ｉ ｍａ ｇ ｅ ｒ ｙ ．
一种新的剩余静校正方法研究及应用
一种新的剩余静校正方法研究及应用虞剑;李大卫【摘要】为了弥补自相关方法的缺陷,运用四阶累积量时延估计理论,同时利用高阶统计量对高斯噪声的盲目性和对非高斯噪声也可以起到压制作用的特点,对常规的最大能量法剩余静校正进行改进,得到一种基于四阶累积量时间延迟估计的剩余静校正处理新方法。
该方法在处理低信噪比、大剩余静校正量资料时更能体现出独到优势,更适用于山地等复杂地形地震资料的处理工作。
和传统方法相比,基于四阶累积量时间延迟估计的剩余静校正方法除在处理速度上略逊于自相关方法外,处理效果更好,鲁棒性更强。
%In order to fill a gap of autocorrelation method,this paper applies fourth-order cumulant to estimate time delay.At the same time,in light of the character of high order statistic＇s being insensitive to Gaussian noise,a new residual static correction method of estimating the time delay by the fourth-order cumulant is achieved through a certain improvement on the conventional residual static correction based on maximum-energy method.The method has the special predominance in processing the seismic data with low signal-to-noise ratio under the condition of complicated near surface or with relevant Gaussian noise and large residual pared with the traditional last static correct method,the model experiment and practical data process indicates that the method performance is excellent and it is more stable and robust excepted processing speed.【期刊名称】《重庆科技学院学报（自然科学版）》【年(卷),期】2011(013)006【总页数】4页(P108-110,123)【关键词】剩余静校正;高阶统计量;互相关;时间序列;时间延迟【作者】虞剑;李大卫【作者单位】西南石油大学,成都610500;西南石油大学,成都610500【正文语种】中文【中图分类】P631剩余静校正是地震资料处理流程中的一个重要环节。
GeoEast初至波剩余静校正
4
初至时间读取 FBTimeRead3D
计算 应用 FBDecom3D StApply
中国石油东方地球物理公司远程培训
概述
所谓的数据准备即为初至波拾取提供输入数据，此时用户需要在完成观 测系统定义的地震数据上进行野外校正量的直接应用（direct application of static correction选件），并得到完成野外校正量应用 之后的地震数据，然后在该数据上进行初至波拾取，作业例子如下：
输入数据：道头中含有初至时间，按inline依次排序的三维/二维CMP道集 输出数据：初至时间和相应道头信息的数据表
中国石油东方地球物理公司远程培训
25
FBTimeRead3D模块输出初至时间和相应道头信息的数据表，并保存在 对应项目、工区（线束）下的datatable目录中。模块参数取值如下图所 示，输出数据表文件名称：fbt－all。
中国石油东方地球物理公司远程培训
13
在静校正量窗口中移动光标，将激活 光标所在处炮的初至旅行时曲线（变 为红色）。
在控制点初至显示窗中，使用鼠标左键点击并拖动， 选择该控制点需要参与计算的偏移距范围。
中国石油东方地球物理公司远程培训
14
中国石油东方地球物理公司远程培训
15
用户可以根据自己需要进行下一个控制点的选择，在主窗口中滑动鼠标， 此时显示为黄色，如上图所示，在需要增加控制点的位臵点击鼠标左键， 弹出控制点初至显示窗，并选定范围即可完成另一个控制点的定义。此 时模块将根据多各控制点定义的范围自动进行插值和外推，形成全线计 算范围，并调整主窗口中初至时间曲线的显示颜色 。 对于已经定义好的控制点，用户可以将鼠标滑动到该控制点的初至时间 曲线，并点击，在弹出的控制点初至显示框进行修改，程序自动进行全 线计算范围的重新计算。 需要说明的时，第一次进入本模块时，无论需要选层与否，都要进行一 次插值计算，即在计算之前选择一次Calculate→ Interpolate，此时程 序按照缺省参数进行计算范围圈定。
剩余静校正
（三）剩余静校正1．定义及分类（1）剩余静校正的定义由于技术上的原因或某些人为因素，例如低速带速度及厚度难以测准，使得野外实测资料往往不很准确，故进行了野外静校正后仍残存着剩余静校正量。
提取剩余静校正量并加以校正的过程叫剩余静校正．．．．．。
因为采用自动统计方法求取剩余静校正量，故也叫自动统计静校正．．．．．．．。
（2）剩余静校正的分类剩余静校正量‖‖短波长剩余静校正量高频剩余静校正量 ++长波长剩余静校正量低频剩余静校正量P125 图6.4-7①长波长或低频剩余静校正量长波长剩余静校正量是区域性的大范围内的异常，地形、低速带变化缓慢，自动统计方法无效。
②短波长或高频剩余静校正量短波长剩余静校正量是局部性的小范围内的异常，地形、低速带变化剧烈,自动统计方法有效。
2．自动统计法求取短波长剩余静校正量的假设和特点自动统计剩余静校正利用多次覆盖资料，所以假设和特点与多次覆盖有针对性。
(1)两点基本假设假设①：波在低速带内垂直于地面传播，即同一炮点或同一接收点的剩余静校正量相同。
假设②：在一个排列的长度上，各炮点（或接收点）的剩余静校正量是随机的，其均值为O 。
(2)三个特点特点①：某个记录道的相对剩余静校正量△τ包括炮点剩余静校正量△τO 和接收点剩余静校正量△τg ，即g O τττ∆+∆=∆特点②：在同一个共炮点道集中，因为各个道的炮点剩余静校正量相同，接收点剩余静校正量不同。
所以将各道的相对剩余静校正量求和再平均，就得到该炮点的剩余静校正量。
地面 V 0低速带底面 V反射界面 例如:一个共炮点道集有24个接收道：第O i 炮的第 1道的剩余静校正量：11g Oi τττ∆+∆=∆ 第O i 炮的第 2道的剩余静校正量：22g Oi τττ∆+∆=∆第O i 炮的第24道的剩余静校正量：2424g Oi τττ∆+∆=∆ 对以上24个式子相加再取平均值有：∑∑∑===∆+∆=∆241241241241241241j j gj Oi j j τττ τ∆ Oi τ∆ 0令左 边 项：∑=∆=∆241241j j ττ右边第一项：∑=∆=∆241241j Oi Oi ττ右边第二项：0241241=∆∑=j gj τ （根据假设②）上式成为：Oi ττ∆=∆上式物理意义：在同一个共炮点道集中，将各道的相对剩余静校正量求和再平均，就得到该炮点的剩余静校正量。
地震数据处理重点整理
地震数据处理重点整理（个人观点）一、题型判断题20分/10个名词解释30分/5个简答题30分/3个计算题20分/2个二、名词解释1、地震剖面的“三高”：高信噪比、高分辨率和高保真度。
2、野外静校正：对陆上资料，把所有炮点和接收点位置均校正到一个公共基准面上以消除高程、低降速带和井深对旅行时的影响。
剩余静校正：野外静校正后，在地震数据中仍然残留有各种剩余静态时移，对这些的校正称为剩余静校正。
3、反褶积：沿时间坐标轴作用，通过压缩地震子波提高地震时间分辨率。
4、最小相位信号：是具有对相同振幅谱的物理可实现信号中相位最小的信号，或者说能量延迟最小的信号。
5、视波数：k=f/v,由于地震勘探是沿测线观测的，因此可以用视波长、视速度、视波数来描述地震波特征，可表示为k*=f/v*,其中k*为视波数。
6、预白化：为了解决带限问题，在地震信号的功率谱P（w）中，从低频到高频统一加一白噪。
7、子波整形反褶积：将不同相位的子波转变为最佳子波的反褶积。
8、速度分析：为叠加提供最佳叠加速度的方法。
9、静校正：存在地形起伏、低速带的厚度变化和速度的横向变化等，此时时距曲线发生畸变，对这些因素的校正，称为静校正。
10、动校正：在水平界面的情况下，从观测到的反射波旅行时中减去正常时差t，得到相当于自激自收的时间，这一过程叫做动校正。
11、正常时差：在界面水平时，对界面上某点以跑检距进行观测得到的反射波旅行时与自激自收观测的旅行时之差，称为正常时差。
12、拉伸畸变：动校正结果出现频率畸变，同相轴移向低频。
13、水平叠加：水平叠加是将CMP道集记录经NMO动校后叠加起来，目的是压制随机噪音，提高地震信噪比。
14、速度谱：把每一种速度所得的叠加结果并排显示在速度-双程零炮检距时间平面中，称此为速度谱。
15、速度扫描：应用一系列常速度在CMP道集进行动校正，并将结果并列显示，从中选出能使反射波同相轴拉平程度最高的速度作为NMO速度的速度分析方法称为速度扫描。
剩余静校正
• • （3-15）
• 式中，vstk为该排列范围内与CMP道集时间最接近的某条双曲线 的速度。
• 方程（3-15）描述的最佳叠加双曲线不一定是方程（3-15）给 出的小排列双曲线。
• 单层近地表模型的复杂性可归结为： • 风化层速度的横向变化； • 接收点和炮点位置高程的快速变化； • 折射层几何形态的横向变化，对折射静校正来说，它被定义为
基岩以上与风化层之间的分界面。
• 风化层底的折射能量经常包含共炮点道集 最先到达波，这些初至波的波前叫做初至。 初至的不同质量一定程度上依赖于震源
• (a) 从未进行处理的野外记录拾取并编辑初至；
• (b) 假设或从测井信息得到一个风化层速度值；
• (c) 对于井下震源应用测井校正；
• (d) 用折射静校正法在所有炮检点计算基岩速度和截距时间，如广义相遇法或最小平方法；
• (e) 用风化层速度、基岩速度和截距时间，在炮点和检波点位置计算基岩深度（方程3-53a,b）；
类型和近地表情况。图3.4-3中的共炮点道 集的初至有明显的起跳。线性初至时间的 偏离大多是由沿着测线高程变化引起的。
• •
• 静校正的处理流程
•
重新看一下图3.3-12的处理流程和图3.4-10的严格描述时差和静校正的近地表模型是很重要的。从未处理的野外记录开始，如图
3.3-12所示的详细的处理流程如下：
• 参看图3.1-16所示的旅行时可得以下结论： • (a) 方程（3-14）中，观察的双程零炮检距时
间OC=t（0），它可能与对应的最佳双曲线拟 合（方程3-15）的双程时间OB=tstk（0）不相 等。例如，在反射层以上存在某种速度非均匀 性时，会出现这种情况；
6.静校正
5.4地表一致性剩余静校正
基本概念
在地表一致性假设的前提下，经过野外静校正和动校正之后， 反射时差可以表示为4个分量之和 ' tij = si + g j + ek + M k x 2ij 式的4个分量中，后两个随反射时间(层位)的变化面变化，前两 个具有地表一致性特征，是我们要计算的炮点和检波点剩余静 校正量。 基于时差分解的剩余静校正方法一般分为三个步骤，首先 拾取每个地震道的时差；然后对时差进行分解，得到炮点和检 波点的剩余静校正量；最后在每个地震道上应用炮点和检波点 静校正量。
5.3初至折射静校正
A D G
加减法折射静校正
加减法折射静校正示意图
θc
B C H E F
图中有三个炮点检波点对，分别是 A → D 、 D → G 和 A → G ，现在定义 两个时间变量 t+ 和 t− 。
t+ = t ABCD + t DEFG − t ABFG
t+ = 2( zw z tan θ c CD CH ) = 2( ) − − w vw vb vw cos θ c vb v sin θ c = w vb
5.1静校正相关概念
基本概念
谢里夫(Sheriff)对静校正所做的定义为：用于补偿由于地 表高程变化、风化层的厚度和速度变化对地震资料的影响。其 目的是获得在一个平面上进行采集，且没有风化层或低速介质 存在时的反射波到达时间。 静校正通常称为地表一致性静校正。地表一致性指某一地 震道的静校正量只与炮点和检波点的地表位置有关，共炮点道 集有着相同的炮点静校正量，共检波点道集有着相同的检波点 静校正量。 静校正方法分类： （1）基准面静校正（野外静校正） （2）初至折射静校正 （3）反射波地表一致性剩余静校正
多道统计剩余静校正
多道统计剩余静校正X刘喜祥1,胡育波2(1.陕西延长石油集团研究院;2.东方地球物理公司长庆物探处,陕西西安 710021) 摘 要:针对短波长剩余静校正量在一共炮点或共深度点道集内是随机量的特点,本文从理论上提出了多道统计剩余静校正的原理,用互相关法求相对剩余静校正量,并用统计法分离,最后可求出记录道的剩余静校正量。
通过实际应用表明,该方法运算简单,工作效率高,有效的解决了短波长引起的误差,达到了精细勘探的目的。
关键词:多道统计;道集;短波长;剩余静校正 中图分类号:P631.4+43 文献标识码:A 文章编号:1006—7981(2012)04—0041—02 由于短波长剩余静校正量在一共炮点或共深度点道集内是随机的量,故采用剩余静校正量统计法分离出炮点和接收点的剩余静校正量。
求其二者的代数和,即得一道记录的剩余静校正量,然后把以基准面为标准换为以强反射层(即反射同相轴)为标准,用互相关法求相对剩余静校正量,并用统计法分离,最后可求出记录道的剩余静校正量,这就是该方法的基本思想。
1 基本原理本方法用多道同相叠加来获取模型道,当道数N足够大时,各道静校正值之和趋于零,该模型道的峰值时间位置必然就是零校正基准线。
所以多道模型道的实质就是用统计的方法求取静校正的零校正基准线,各道直接向此零校正线对齐,就达到了静校正的目的。
对于多次覆盖资料,这种校正必须完成两方面的任务:共深度点(CDP)道集内各道都应向模型道对齐,以保证叠加后反射波的能量增强。
为了保证CDP道集间光滑连续,必须使每一个CDP道集向模型道对齐后,其时间位置要正确的落在反射界面上。
2 方法步骤2.1 参考道的形成2.1.1 选定标准层。
静校正的假设条件表明了静校正量只与空间位置有关,而与时间无关,提取静校正量时,只要在记录中选择一个好的标准层就可以了。
用人工法选取时,可在初叠剖面上用人工对比解释的方法,选取波组连续性好、波形稳定、倾角小和能量强的反射波组作为标准层见(图1)。
GeoEast系统反射波剩余静校正
GeoEast系统反射波剩余静校正
中国石油东方地球物理公司远程培训
主要内容
概述 二维地表一致性剩余静校正 三维地表一致性剩余静校正 非地表一致性剩余静校正
2
中国石油东方地球物理公司远程培训
3
由于多种因素，一个CMP道集中的各道，经过野外静校正或者初至静 校后，仍然存在剩余静校正量，影响CMP叠加质量。因此在CMP叠加之 前，我们必须对它进行估算和校正，实现CMP道集的同相叠加。 GeoEast 处理系统中，对于反射波剩余静校正主要由： •剩余静校正量计算（RsStCalculat） •三维剩余静校正量调整（StHarmoniz3D） •模型道剩应用余时差拾取（SCRsCal3D） •三维地表一致性剩余时差分解（SCRsDecom3D） •叠前CMP道集时差微调（Trim3D） •静校正量应用（StApply）等模块来完成。
中国石油东方地球物理公司远程培训
13
一是控制点方式，如下图所示，在全线选择数个控制点，并填写相应的倾角 和时窗参数，控制点之间程序自动插值形成全线的计算时窗：
二是区域划分法，如下图所示，将整条二维测线划分成若干段，每段使用一 组时窗控制参数来定义计算时窗，段间保留一定的空间用于时窗插值。
中国石油东方地球物理公司远程培训
中国石油东方地球物理公司远程培训
5
输入数据描述
输入数据为动校正后的CMP道集数据。
输出数据描述 反射波剩余校正量，保存在数据库中或者数据表中供后续应用模块使用。
中国石油东方地球物理公司远程培训
主要内容
概述 二维地表一致性剩余静校正 三维地表一致性剩余静校正 非地表一致性剩余静校正
6
中国石油东方地球物理公司远程培训
maximum source number 处理区域中的最大炮号； minimum receiver number 处理区域中的最小接收点号；
模拟退火法剩余静校正在山地资料处理中的应用
低、 剩余静校 正量大的 山地资料处理 中效果 明显。 在 常规 地表一致性剩余静 校正方法失效 时依然有 效， 提高了剖面成像质量。但在使用该方法时 ， 需注
求 出的解不是最优解 ； 太大 ， 影响计算效率 ， 实践证 明， 迭代次数一般为６ 次就足够了。 ０
（ 下转第 ６ 页 ） ８
从 图 ２ 可 以看 出 ， ２ｂ 的处 理 效 果 明显 优 中 图 （）
在剩 余静 校正 里 ， 通过 给定 （ ， 的一个初 ｓ Ｒ）
值， 通过上述计算过程 ， 不断对初值进行扰动 ， 利用 概率分布 函数选取符合式 （） １的最优解 的扰动 ， 逼近 最优解 ， 求得全局收敛 的剩余静校正量 。
关键词 ： 拟退 火法； 模 剩余 静校正 ； 山地资 料处理
模 拟 退 火 法 （ｉｕａ ｄＡ ｎａｉｇＭｅｈｄ 简 称 Ｓｍ ｌｅ ｎ ｅｌ ｔｏ ， ｔ ｎ
出来的模型参数 ， 而对反演结果 的接受 可简化为最
优 化 问题 【 即 ｍｉＥＳ Ｒ 。 ３ ． ［， Ｊ Ｊ ｎ ，
．
一
泥 页岩 、 灰岩 为主 ， 泥 地表 地层 褶皱 扰动 ， 否则不接受该扰动。每 当一个参数 的扰动值 叠 系 的砂 岩 、 岩 、 地震地质条件差 ， 资料信噪比低 。该区深部地 被接受 ， 就修改该参数 的值 ， 所有后来 目 函数值的 剧烈 ， 标 震 地 质 条 件 复杂 ， 层 倾 角 大 ， 地 断裂 发 育 ， 在各 种 存 计算都根据该参数值 ， 直至它再次被修改为止 。
油 气 地 球 物 理
２ １年 ４ ０２ 月 Ｐ Ｔ ＯＬ ＵＭ ０Ｐ ＹＳＣ ＥＲ Ｅ ＧＥ Ｈ Ｉ Ｓ 第 １卷 第 ２ ０ 期
模拟退火法剩余静校正在山地资料处理中的应用
第八节 地表一致性反射波剩余静校正
第八节 地表一致性反射波剩余静校正许多静校正方法均采用地表一致性模型。
检波器组在位置i 的延迟G i 和震源在位置j 的延迟S j ，对所有相应的地震道都是相同的。
如果i 和j 具有公共的坐标原点，其炮检距正比于（j －i ）。
如果沿测线有构造，我们以CMP （共中心点）号位置k 来表示延迟量L k ，其中k ＝（j - i ）/2，表示构造的深度比其它位置深L k 个单位，L k 是构造时移的均值。
对于平界面反射而言，它即指向中点位置。
如果倾角很缓，L k 对于共中心点来说几乎是常数。
如果动校正速度有误差，就会保留一些剩余时差M k ，它随炮检距的平方而变化。
如果不考虑炮点或检波点对测线的横向偏离，那么对于地表一致性模型，一个道总的时移量t i,j 为： t i,j =G i 十 S j ＋ L k 十 Mk （j – I ）2 (8.1)地表一致性模型不只限于确定静校正量的时移，在振幅调节、子波提取、反褶积及其它算法中，有时候都是基于地表一致性模型，均按上述相同的过程进行。
我们并不知道每一道的时移量，但可以利用互相关求一个道相对于另一个道的时移量（t i,j – t m,n ）：t i,j – t m,n =G i –G m +S j –S n +L i+j – L m+n +M i+j （j - i ）2 – M m+n （n –m ）2 （8.2）在此没有使用下标k ，而使用i+j ，是确保下标是整数，具体数值无关紧要。
对于一个CMP 道集，两两道互相关，就可以得出比未知数（G i ，S j ，L i+j ，M i+j ）个数还要多的方程式，是一个超定方程组。
但是，方程(8.2)在测量过程中也包含一些不确定的因素，例如等号左边就会出现误差。
一个“超定”，一个“不确定”，就使我们有办法来求解这个方程组，通常用最小平方法，有时也用迭代法。
最小平方问题是使误差e p 的平方和最小：()[]{}最小值∑∑=+-+-+---==++++2n m 2j i n m ,,2m) - (n M - i -j M L )(j i n j m i n m j i p L S S G G t t e E (8.3) 可用下式求解：,0,0,0,0=∂∂=∂∂=∂∂=∂∂++j i j i j i M E L E S E G E (8.4)通常，我们事先构成一个横型道，例如经过一般处理后的本道集的初步叠加道，或者是经过时移处理后的前一个道集的叠加道。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
对近地表结构及观测点
高程变化所引起的地震 记录 时间延迟的补偿
（值大、长波长）
6
6 东方公司研究院处理中心
剩余静校正质控点
前言/野外．剩余
野外静校正定义
R
S
R´
地面
风化层底
S´ R´
S´
基准面
反射面
（a）
（b）
基准面校正示意图
（a）实际模型（炮点S、检波点R 在地面上） （b）设想模型（炮点S´、接收点R´ 在基准面上）
7
7
东方公司研究院处理中心
剩余静校正质控点
前言/野外．剩余
野外静校正定义
8
8
东方公司研究院处理中心
剩余静校正质控点
前言/野外．剩余
剩余静校正定义
剩余静校正是对基准面校正 误差的一种补充校正
其目的是保证一个CMP道集的反射同相 轴在动、静校正之后能很好的对齐，以 便提高道集和叠加数据的信噪比，最大 程度地改善叠加效果。 （值小、短波长）
剩余静校正质控点
IL21线第一次剩余静校正前的叠加
质控分析／实例－1
IL21线第一次剩余静校正后的叠加
54
54
东方公司研究院处理中心
剩余静校正质控点
IL21线第二次剩余静校正前的叠加
22
东方公司研究院处理中心
剩余静校正质控点 工程化质量控制系列讲座
模块举例
剩余静校正模块举例 a)旅行时分解法剩余静校正（MISER-Ω )
b) 叠加能量最大法剩余静校正(TDSAT-CGG)
C) 修饰静校正(TRIM- Ω) 23
23
东方公司研究院处理中心
剩余静校正质控点 工程化质量控制系列讲座
基本原理
E=∑ (Tij-T ij) i,j
,
2
E值最小，表明观测时间与求取时间的误差
最小，体现NMO的道集同相轴叠加能量最大
适应资料：覆盖次数低
校正量不大
34
34
S/N 高
东方公司研究院处理中心
剩余静校正质控点 模块举例/TDSAT 工程化质量控制系列讲座
相
关
法
2
基 本 原 理
P（△t）= Σ｛ F（t - △t）+ G（t）｝
参 数 选 择
44
44
计算时窗大一些 相关时移小一些（＜12ms) 数据准备好一些(NMO速度、地表一致性静校正）
效果检查细一些(构造及波组特征正常？）
东方公司研究院处理中心
剩余静校正质控点
目录
一 前言
二、剩余静校正模块举例 三、剩余静校正质控分析
45
45 东方公司研究院处理中心
剩余静校正质控点
东方公司研究院处理中心
剩余静校正质控点 模块举例/TDSAT 工程化质量控制系列讲座
基本原理
以定义的反射同相轴取得最佳叠加效果（振 幅强且连续性好）为标准来计算地表一致性 的炮点、检波点剩余静校正量（剩余动校正 量也可参与计算）
33
33 东方公司研究院处理中心
剩余静校正质控点
模块举例/TDSAT
28
28
东方公司研究院处理中心
模块举例/MISER 剩余静校正质控点 工程化质量控制系列讲座
构 造 项 平 滑
29
29
东方公司研究院处理中心
剩余静校正质控点 工程化质量控制系列讲座
模块举例
剩余静校正模块举例 a)旅行时分解法剩余静校正（MISER-Ω )
b) 叠加能量最大法剩余静校正(TDSAT-CGG)
• 最大校正值限制
• 校正值平滑因子长度
37
37 东方公司研究院处理中心
剩余静校正质控点 工程化质量控制系列讲座
模块举例
剩余静校正模块举例 a)旅行时分解法剩余静校正（MISER-Ω )
b) 叠加能量最大法剩余静校正(TDSAT-CGG)
C) 修饰静校正(TRIM- Ω) 38
38
东方公司研究院处理中心
48 东方公司研究院处理中心
剩余静校正质控点
质控分析／要求
1)用于NMO的速度场是否合适？
质 控 要 求
2)拾取时窗是否为高S/N数据段？
3)时移限制量选择是否有依据？
4)剩余静校正量有无局部异常？ 5)剩余静校正量是否收敛？ 6)剩余静校正对道集、叠加的 S/N是否有改善？
49Βιβλιοθήκη 49东方公司研究院处理中心
9
9 东方公司研究院处理中心
剩余静校正质控点
前言
• 野外／剩余静校正 • 长、短波长分量及SCMP值
• 地表一致性静校正
• 剩余静校正（与V分析）迭代
10
10 东方公司研究院处理中心
剩余静校正质控点
前言/长、短波长
λ ＞ L －长波长； λ <0.5 L –短波长
排列长度 L
短波长分量
长波长分量
除互相关权 不用于模型 建立或叠加 外，时差拾 取与MISER 拾取步相似
42
42
东方公司研究院处理中心
剩余静校正质控点 模块举例/TRIM 工程化质量控制系列讲座
TRIM- Ω
• 基本原理 • 参数选择
43
43
东方公司研究院处理中心
剩余静校正质控点 模块举例/TRIM 工程化质量控制系列讲座
模型范围小一些（3、5个面元亦可）
2
(旅行时方程)
a.方程数多于未知数 b.独立方程数少于未知数 c.方程组无唯一解，只有近似解 d.高斯－赛得尔迭代求解，分步实现 26
26 东方公司研究院处理中心
剩余静校正质控点 模块举例/MISER 工程化质量控制系列讲座
MISER-Ω
• • 基本原理 参数选择
27
27
东方公司研究院处理中心
剩余静校正质控点 模块举例/MISER 工程化质量控制系列讲座
剩余静校正质控点 模块举例/TRIM 工程化质量控制系列讲座
TRIM- Ω
• 基本原理 • 参数选择
39
39
东方公司研究院处理中心
剩余静校正质控点 模块举例/TRIM 工程化质量控制系列讲座
TRIM- Ω
• 基本原理 • 参数选择
40
40
东方公司研究院处理中心
剩余静校正质控点
模块举例/TRIM
剩余静校正质控点
质控分析
剩余静校正质控分析
. . .
50
50
质控要求 质控图件 质控实例分析
东方公司研究院处理中心
剩余静校正质控点
质控分析／图件
1)剩余静校正参数选择对比图
质 控 图 件
51
51
2)最终炮、检剩余静校正量（平面、剖面）图
3)（应用剩余静校正前后的）静校正迭代效果
比较图：
a) 叠加数据（含三维纵、横控制线）剖面比较显示
•长、短波长分量及SCMP值
•地表一致性静校正 •剩余静校正（与V分析）迭代
4
东方公司研究院处理中心
4
剩余静校正质控点
前言
• 野外／剩余静校正 • 长、短波长分量及SCMP值
• 地表一致性静校正
• 剩余静校正（与V分析）迭代
5
5 东方公司研究院处理中心
剩余静校正质控点
前言／野外．剩余
野外静校正定义
FSCS · -0.5 Scmp=炮点高频静校正值 FDCS-0.5 Scmp=检波点高频静校正值 · Scmp=炮点与检波点在该CDP点的低频静校正值之和 · 一个炮点（或检波点）各道的高、低频校正值通常不同 · 在 ·Scmp面上最适合进行速度分析和动校正
13
13
长、短 波长分 量及 SCMP值
基本原理
抛开地表一致性条件约束，利用NMO道集 和模型道数据的互相关来求取各道与面元对
齐的时移。其模型既可是本面元的叠加道，
也可是以它为中心多个面元叠加道的求和道
41
41 东方公司研究院处理中心
剩余静校正质控点 模块举例/TRIM 工程化质量控制系列讲座
● ● ● ● ●
纵向
横向
基 本 原 理
三维（3x3)模型道形成示意图
20
20
东方公司研究院处理中心
剩余静校正质控点
目录
一 前言
二、剩余静校正模块举例 三、剩余静校正质控分析
21
21 东方公司研究院处理中心
剩余静校正质控点 工程化质量控制系列讲座
剩余静校正模块举例 a)旅行时分解法剩余静校正（MISER-Ω )
b) 叠加能量最大法剩余静校正(TDSAT-CGG)
C) 修饰静校正(TRIM- Ω) 22
11
11
东方公司研究院处理中心
剩余静校正质控点
前言/长、短波长
Scmp值
RG－cor
Scmp =以k为中心局部范围内
k
(
s R )的平均值
i+ j
CMP－datcor
SCMP
1 ( Si R j ) N i, j
东方公司研究院处理中心
12
12
剩余静校正质控点
前言/长、短波长
• Scmp值是一个程双时间
模块举例
MISER-Ω
• • 基本原理 参数选择
24
24
东方公司研究院处理中心
剩余静校正质控点 模块举例/MISER 工程化质量控制系列讲座
MISER-Ω
• • 基本原理 参数选择
25
25
东方公司研究院处理中心
剩余静校正质控点
模块举例/MISER
基本原理
T
’
i j
=
Si