拟合初至波获取短波长静校正量的方法
【专利摘要】本发明提供了一种拟合初至波获取短波长静校正量的方法。所述方法包括：A)选取测线上的接收点，获取和所述接收点对应的炮点集合；B)拾取对应于所述接收点和炮点集合的共接收点道集的第一初至时间集；C)对所拾取的共接收点道集的第一初至时间集进行样条约束拟合以得到样条约束拟合后的第二初至时间集；D)基于所述共接收点道集的第一初至时间集和样条约束拟合后的第二初至时间集获取所述接收点的剩余静校正量。通过本方法可使地震剖面成像效果好并且计算效率高。
【专利说明】拟合初至波获取短波长静校正量的方法
【技术领域】
[0001]本发明涉及地震勘探资料处理领域。更具体地讲，涉及一种获取复杂地区二维资料的静校正量的方法。
【背景技术】
[0002]静校正问题是制约我国复杂地区油气勘探发展的一个瓶颈问题。静校正问题分为长波长、短波长校正问题，而短波长静校正问题解决的好坏直接影响地震剖面成像效果，所以大批的地震勘探工作者对剩余静校正问题进行了大量研究，形成了初至波、折射波、反射波剩余静校正方法。
[0003]基于反射波的剩余静校正方法与速度拾取相结合多次迭代计算，在资料品质较好的地区可得到较好效果，但在诸如山地、沙漠、厚黄土塬等地区，地形起伏大，近地表速度、厚度变化剧烈，激发、接收条件差，资料信噪比低，导致反射波剩余静校正方法因构建的模型道难以见到有效反射而失效，因此在做反射波剩余静校正之前有必要做初至波、折射波剩余静校正。
[0004]折射波剩余静校正方法在一些地区取得了一定效果，但是这种方法也需要对时窗内的地震数据做互相关，计算量较大，计算耗时长。另外，折射波尽管能量较强、易识别，但有的地震道因为接收条件差，折射波出现奇异值，影响计算的精准性。而初至波剩余静校正方法，可人机联合拾取初至时，去除不可靠的奇异点，保证计算结果准确、可靠。
[0005]基于初至波的剩余静校正方法拾取延迟时的前提是将初至时间拟合成一条趋势线，有学者应用最小二乘法等数学方法将初至时间拟合成一条直线，或分步长拟合成一条折线。前者需要根据偏移距对初至时间进行限制，那么可用的初至时间数量会减少，统计信息的可靠性会降低；后者的不足之处在于计算机无法自动识别直达波、折射波、反射波，无法判断哪些初至时间来自同一层，导致计算结果不是很精确。
【发明内容】
[0006]本发明的目的在于提供一种地震剖面成像效果好并且计算效率高的获取初至波剩余静校正量的方法。
[0007]为了实现上述目的，根据本发明的一方面，提供了一种拟合初至波获取短波长静校正量的方法，包括:A)选取测线上的接收点，获取和所述接收点对应的炮点集合；B)拾取对应于所述接收点和炮点集合的共接收点道集的第一初至时间集；C)对所拾取的共接收点道集的第一初至时间集进行样条约束拟合以得到样条约束拟合后的第二初至时间集；D)基于所述共接收点道集的第一初至时间集和样条约束拟合后的第二初至时间集获取所述接收点的剩余静校正量。
[0008]可选地，步骤D)可包括:D1)将所述共接收点道集的第一初至时间集包括的每个初至时间和拟合后的第二初至时间集中相应的初至时间做差，并将每个做差后的差值累加；D2)基于与所述接收点对应的炮点数对Dl)中所累加的差值做平均，获取所述接收点的剩余静校正量。
[0009]可选地，所述的方法还可包括:E)检测所述接收点是否是测线上最后的点，当所述接收点不是测线上最后的点时，改变接收点，返回执行步骤A)以对测线上的接收点进行遍历，并获取测线上所有接收点的剩余静校正量。
[0010]根据本发明的另一方面，提供了一种拟合初至波获取短波长静校正量的方法，包括:A)选取测线上的炮点，获取和所述炮点对应的接收点集合；B)拾取对应于所述炮点和接收点集合的共炮点道集的第一初至时间集；C)对所拾取的共炮点道集的第一初至时间集进行样条约束拟合以得到样条约束拟合后的第二初至时间集；D)基于所述共炮点道集的第一初至时间集和样条约束拟合后的第二初至时间集获取所述炮点的剩余静校正量。
[0011]可选地，步骤D)可包括:D1)将所述共炮点道集的第一初至时间集包括的每个初至时间和拟合后的第二初至时间集中相应的初至时间做差，并将每个做差后的差值累加；D2)基于与所述炮点对应的接收点数量对Dl)中所累加的差值做平均，获取所述炮点的剩余静校正量。
[0012]可选地，所述的方法还可包括:E)检测所述炮点是否是测线上最后的点，当所述炮点不是测线上最后的点时，改变炮点，返回执行步骤A)以对测线上的炮点进行遍历，并获取测线上所有炮点的剩余静校正量。
[0013]将在接下来的描述中部分阐述本发明另外的方面和/或优点，还有一部分通过描述将是清楚的，或者可以经过本发明的实施而得知。
【专利附图】
【附图说明】
[0014]通过下面结合附图进行的详细描述，本发明的上述和其它目的、特点和优点将会变得更加清楚，其中:
[0015]图1示出了根据本发明示例性实施例的拟合初至波获取接收点短波长静校正量的方法的流程图；
[0016]图2示出了根据本发明示例性实施例的拟合初至波获取炮点短波长静校正量的方法的流程图。
【具体实施方式】
[0017]以下，将参照附图更充分地描述本发明的示例性实施例，示例性实施例在附图中示出。然而，可以以许多不同的形式实施示例性实施例，并且不应被解释为局限于在此阐述的示例性实施例。相反，提供这些实施例从而本公开将会彻底和完整，并将完全地将示例性实施例的范围传达给本领域的技术人员。
[0018]图1示出了拟合初至波获取接收点短波长静校正量的方法的流程图。
[0019]参照图1，在步骤S101，选取测线上的接收点，获取和所述接收点对应的炮点集
八
口 ο
[0020]例如，可通过在地震资料的炮点集合信息中检索具有所述接收点号的炮点的方式来获取与所述接收点对应的炮点集合。
[0021]在步骤S102，拾取对应于所述接收点和炮点集合的共接收点道集的第一初至时间集。[0022]在初至时间信息中检索同时具有所述接收点号和在步骤SlOl中获取的炮点集合包括的炮点号的初至时间，以获取第一初至时间集。
[0023]在步骤S103，对所拾取的共接收点道集的第一初至时间集进行样条约束拟合以得到样条约束拟合后的第二初至时间集。
[0024]例如，可以使用炮点坐标作为变量建立初至时间和炮点坐标的函数，对初至时间和炮点坐标的函数进行样条约束拟合。
[0025]在步骤S104，基于所述共接收点道集的第一初至时间集和样条约束拟合后的第二初至时间集获取所述接收点的剩余静校正量。
[0026]例如，可将所述共接收点道集的第一初至时间集包括的每个初至时间和拟合后的第二初至时间集中相应的初至时间做差，并将每个做差后的差值累加；可基于与所述接收点对应的炮点数对所累加的差值做平均，获取所述接收点的剩余静校正量。[0027]作为一个示例，第一初至时间集和第二初至时间集包括一一对应的元素，可将第二初至时间集中的每个元素减第一初至时间集中相应的元素得到差值，将每个元素相减的差值累加后除以炮点数获得平均值，所获得的平均值是所述接收点的剩余静校正量。
[0028]此外，所述方法还可包括以下步骤(未示出):可检测所述接收点是否是测线上最后的点，当所述接收点不是测线上最后的点时，改变接收点，返回执行步骤S102至步骤S104以对测线上的接收点进行遍历，并获取测线上所有接收点的剩余静校正量。
[0029]下面将示出拟合初至波获取短波长静校正量的方法的详细实施例以对步骤SlOl至步骤S104进行说明。
[0030]作为一个示例，以一个接收点&为例，首先获取与该接收点对应的所有炮点S”i=l, 2，…，N。
[0031]拾取接收点&与各炮点所对应的共接收点道集的初至时间Tijt5
[0032]提取所有所对应的共炮点道集，对道集中的初值时间进行拟合处理，得到接收点Rj与各炮点所对应的共接收点道集新的初至时间T’ ij。
[0033]拟合的具体实施方法如下:
[0034]将坐标点及初至时间按如下表示
[0035]f (Xi)=Yi, i=0, I,..., η
[0036]其中，Xi表示炮点坐标，Yi表示对应的初至时间，且令
[0037]x0<x1<...<χη
[0038]要拟合函数f (Xi)=yi得到拟合后的函数g(x)=y须满足
[0039]min Jg"(x)2dx(I)
[0040]且带有约束
[0041]
【权利要求】
1.一种拟合初至波获取短波长静校正量的方法，包括: A)选取测线上的接收点，获取和所述接收点对应的炮点集合； B)拾取对应于所述接收点和炮点集合的共接收点道集的第一初至时间集； C)对所拾取的共接收点道集的第一初至时间集进行样条约束拟合以得到样条约束拟合后的第二初至时间集； D)基于所述共接收点道集的第一初至时间集和样条约束拟合后的第二初至时间集获取所述接收点的剩余静校正量。
2.根据权利要求1所述的方法，其中，步骤D)包括: Dl)将所述共接收点道集的第一初至时间集包括的每个初至时间和拟合后的第二初至时间集中相应的初至时间做差，并将每个做差后的差值累加； D2)基于与所述接收点对应的炮点数对Dl)中所累加的差值做平均，获取所述接收点的剩余静校正量。
3.根据权利要求1所述的方法，还包括: E)检测所述接收点是否是测线上最后的点，当所述接收点不是测线上最后的点时，改变接收点，返回执行步骤A)以对测线上的接收点进行遍历，并获取测线上所有接收点的剩余静校正量。
4.一种拟合初至波获取短波长静校正量的方法，包括: A)选取测线上的炮点，获取和所述炮点对应的接收点集合； B)拾取对应于所述炮点和接收点集合的共炮点道集的第一初至时间集； C)对所拾取的共炮点道集的第一初至时间集进行样条约束拟合以得到样条约束拟合后的第二初至时间集； D)基于所述共炮点道集的第一初至时间集和样条约束拟合后的第二初至时间集获取所述炮点的剩余静校正量。
5.根据权利要求4所述的方法，其中，步骤D)包括: Dl)将所述共炮点道集的第一初至时间集包括的每个初至时间和拟合后的第二初至时间集中相应的初至时间做差，并将每个做差后的差值累加； D2)基于与所述炮点对应的接收点数量对Dl)中所累加的差值做平均，获取所述炮点的剩余静校正量。
6.根据权利要求4所述的方法，还包括: E)检测所述炮点是否是测线上最后的点，当所述炮点不是测线上最后的点时，改变炮点，返回执行步骤A)以对测线上的炮点进行遍历，并获取测线上所有炮点的剩余静校正量。
【文档编号】G01V1/36GK103837895SQ201410086448
【公开日】2014年6月4日 申请日期:2014年3月10日 优先权日:2014年3月10日
【发明者】周强, 陈爱萍, 何光明, 曹琳昱, 杨海涛 申请人:中国石油集团川庆钻探工程有限公司地球物理勘探公司