表层模型影响的高频成分中的时间校正量，利用（7)式计算表层 模型影响的高频成分中的高程改变量。
[0122] 图5、图6、图7和图8分别是计算得到的表层模型低频成分的顶界高程、底界高 程、平均速度和厚度，图9和图10分别是计算得到的表层模型影响的高频成分的时间校正 量和高程改变量。
[0123] 实施例2
[0124] 表层模型低频成分的平滑范围为以待计算测点为中心半径为R = 1000 m的圆，采 用一种随与待计算测点距离变化的权系数
[0125]
(8)
[0126] 式中I η-χ|是平面上Π点与X点间的距离，kfl/OR)、km=l/(8R)。
[0127] (1)计算表层模型低频成分的顶界高程
[0128] 采用（3)式计算表层模型低频成分的顶界高程，其中A(X)是落入以X点为中心半 径R= 1000 m的平滑范围内测点的集合，加权系数W(Ii-X)由⑶式计算。
[0129] (2)计算表层模型低频成分的底界高程
[0130] 采用（4)式计算表层模型低频成分的底界高程，其中A(X)是落入以X点为中心半 径R= 1000 m的平滑范围内测点的集合，加权系数W(Ii-X)由⑶式计算。
[0131] (3)计算表层模型低频成分的平均速度
[0132] 采用（5)式计算表层模型低频成分的平均速度，其中A(X)是落入以X点为中心半 径R= 1000 m的平滑范围内测点的集合，加权系数W(Ii-X)由⑶式计算。
[0133] (4)计算表层模型影响的高频成分
[0134] 采用（6)式计算表层模型影响的高频成分中的时间校正量，利用（7)式计算表层 模型影响的高频成分中的高程改变量。
[0135] 图11、图12、图13和图14分别是计算得到的表层模型低频成分的顶界高程、底界 高程、平均速度和厚度，图15和图16分别是计算得到的表层模型影响的高频成分的时间校 正量和高程改变量。
[0136] 实施例3
[0137] 表层模型低频成分的平滑范围是为以待计算测点为中心，半边长分别为dx = 900m、dy = 950m的矩形，矩形对角线长度的一冲
。采用一种随与 待计算测点距离变化的权系数
[0138]
(9)
[0139] 式中I η_χ|是平面上η点与X点间的距离，kfVORhkflASR)。
[0140] (1)计算表层模型低频成分的顶界高程
[0141] 采用（3)式计算表层模型低频成分的顶界高程，其中A(X)是落入以X点为中心半 边长分别为dx = 900m、dy = 950m的矩形内测点的集合，加权系数w( η-χ)由（9)式计算。
[0142] (2)计算表层模型低频成分的底界高程
[0143] 采用（4)式计算表层模型低频成分的底界高程，其中A(X)是落入以X点为中心半 边长分别为dx = 900m、dy = 950m的矩形内测点的集合，加权系数w( η-χ)由（9)式计算。
[0144] (3)计算表层模型低频成分的平均速度
[0145] 采用（5)式计算表层模型低频成分的平均速度，其中A(X)是落入以X点为中心半 边长分别为dx = 900m、dy = 950m的矩形内测点的集合，加权系数w( η-χ)由（9)式计算。
[0146] (4)计算表层模型影响的高频成分
[0147] 采用（6)式计算表层模型影响的高频成分中的时间校正量，利用（7)式计算表层 模型影响的高频成分中的高程改变量。
[0148] 图17、图18、图19和图20分别是计算得到的表层模型低频成分的顶界高程、底界 高程、平均速度和厚度，图21和图22分别是计算得到的表层模型影响的高频成分的时间校 正量和高程改变量。
[0149] 上述技术方案只是本发明的一种实施方式，对于本领域内的技术人员而言，在本 发明公开了应用方法和原理的基础上，很容易做出各种类型的改进或变形，而不仅限于本 发明上述【具体实施方式】所描述的方法，因此前面描述的方式只是优选的，而并不具有限制 性的意义。
【主权项】
1. 一种表层模型对地震数据影响的低频与高频成分的分解方法，其特征在于：所述方 法包括： (1) 输入表层模型数据、平滑范围和工区替换速度参数VR; (2) 计算表层模型低频成分的顶界高程，即平滑基准面高程； (3) 计算表层模型低频成分的底界高程； (4) 计算表层模型低频成分的平均速度； (5) 计算表层模型影响的高频成分，即时间校正量和高程改变量； (6) 重复步骤（2)到（5)完成所有测点的计算，然后输出所述低频成分的顶界高程、底 界高程、平均速度、表层模型影响的高频成分的时间校正量和高程改变量。2. 根据权利要求1所述的表层模型对地震数据影响的低频与高频成分的分解方法，其 特征在于：所述步骤（1)中，所述表层模型数据包括测点平面坐标、顶界高程、底界高程、平 均速度；所述平滑范围是计算表层模型低频成分的平滑范围； 实际表层模型由顶界高程、底界高程和平均速度3个参数描述，顶界高程即测点高程， 平均速度是顶界与底界之间介质的平均速度；表层模型低频成分同样由顶界高程、底界高 程和平均速度3个参数描述； 设平面位置x测点处实际表层模型的顶界高程、底界高程、平均速度分别为es(x)、eb(x)、v(x)，要计算的表层模型低频成分的顶界高程、底界高程、平均速度分别用Es(x)、 Eb (x)、V(x)表示； 给定的平滑范围是一个二维窗口，当地震勘探是一条二维测线时，平滑范围退化为一 维窗口； 相对于待计算测点，平滑范围包含待计算测点，并随着待计算测点位置的变化平移。3. 根据权利要求2所述的表层模型对地震数据影响的低频与高频成分的分解方法，其 特征在于：所述步骤（2)是这样实现的： 表层模型低频成分的顶界高程Es (x)是给定平滑范围内实际表层模型静校正量在给定 平滑范围内的加权平均值所对应的表层模型的顶界高程，即相对于x点的平滑范围内测点 实际表层模型顶界高程的加权平均值：(3) 式中A(x)是落入相对于x点的平滑范围内的测点的集合，w(n-x)是加权系数。4. 根据权利要求2所述的表层模型对地震数据影响的低频与高频成分的分解方法，其 特征在于：所述步骤（3)是这样实现的： 表层模型低频成分的底界高程Eb (x)是给定平滑范围内实际表层模型静校正量在给定 平滑范围内的加权平均值所对应的表层模型的底界高程，即相对于x点的平滑范围内测点 实际表层模型底界高程的加权平均值：(4) 式中A(x)是落入相对于x点的平滑范围内的测点的集合，w(n-x)是加权系数。5. 根据权利要求2所述的表层模型对地震数据影响的低频与高频成分的分解方法，其 特征在于：所述步骤（4)是这样实现的： 表层模型低频成分的平均速度V(x)是给定平滑范围内实际表层模型静校正量在给定 平滑范围内的加权平均值所对应的表层模型的速度，即相对于x点的平滑范围内测点实际 表层模型的平均速度：(5) 式中h(x) =es(X)-eb(X)是实际表层模型的厚度，A(x)是落入相对于x点的平滑范围 内的测点的集合，w(n-x)是加权系数。6. 根据权利要求3至5任一所述的表层模型对地震数据影响的低频与高频成分的分解 方法，其特征在于：所述加权系数W(n-x)是一个常数或者是跟n与X的相对关系有关的 函数； 当w(n-x)是跟n与X的相对关系有关的函数时，w(n-x)是随n与X之间的距离 增大而减小的函数，且w(n-x)大于等于〇。7. 根据权利要求2所述的表层模型对地震数据影响的低频与高频成分的分解方法，其 特征在于：所述步骤（5)是这样实现的： 表层模型影响的高频成分是实际表层模型影响与表层模型低频成分影响的差异在时 间校正量与高程改变量上的体现；表层模型影响高频成分中的时间校正量是实际表层模型 静校正量与表层模型低频成分静校正量的差： At(x) = - (es (x) -eb (x))/v(x) + (Es (x) -Eb (x))/V(x) + (Eb (x) -eb (x))/vR(6) 表层模型影响的高频成分中的高程改变量是表层模型低频成分顶高程与实际表层模 型顶高程之差： Az(x) =Es(x)-es(x) (7)。
【专利摘要】本发明提供了一种表层模型对地震数据影响的低频与高频成分的分解方法，属于地震勘探资料数据处理领域。本方法包括：(1)输入表层模型数据、平滑范围和工区替换速度参数vR；(2)计算表层模型低频成分的顶界高程，即平滑基准面高程；(3)计算表层模型低频成分的底界高程；(4)计算表层模型低频成分的平均速度；(5)计算表层模型影响的高频成分，即时间校正量和高程改变量；(6)重复步骤(2)到(5)完成所有测点的计算，然后输出所述低频成分的顶界高程、底界高程、平均速度、表层模型影响的高频成分的时间校正量和高程改变量。
【IPC分类】G01N1/36, G01V1/28, G01N1/30
【公开号】CN105093280
【申请号】CN201410207196
【发明人】林伯香, 袁联生, 徐颖, 朱海波
【申请人】中国石油化工股份有限公司, 中国石油化工股份有限公司石油物探技术研究院
【公开日】2015年11月25日
【申请日】2014年5月15日