实施例:
如图1所示,为本发明一种优选风化层速度的确定方法的流程示意图,包括以下步骤:
步骤S101:预设代表风化层速度v0;所述预设代表风化层速度v0可选取750米/秒~850米/秒之间的任意一个值,例如:800米/秒。
步骤S102:根据所述预设代表风化层速度v0获取实际风化层速度变化区间,具体过程如下;
1)、根据剩余静校正值δt及截距时t的边界范围值,确定所述剩余静校正值δt、截距时t,在实际的地质环境下,为使剩余静校正量在后续的剩余静校正处理中可以被消除掉,通常剩余静校正值δt具体边界范围绝对值不超过20毫秒;同时对截距时t的取值范围,根据我国地质地貌的实际情况,可具体选取为截距时t的边界范围值不超过400毫秒,即风化层厚度不超过300米,因此,所述剩余静校正值δt的绝对值为20毫秒,所述截距时t的取值为400毫秒。
2)、将所述剩余静校正值δt及截距时t代入如下误差曲线公式中获取误差指数k;
当所述剩余静校正值δt为±20毫秒,及截距时t为400毫秒时,所述误差指数k等于±0.1。
3)、将所述误差指数k及分别将所述预设代表风化层速度v0代入如下误差曲线公式中分别获取剩余风化层速度dv0;如:将2)中获取的k等于±0.1,及在3)中通过预设选取的代表风化层速度预设值v0为800米/秒,代入如下误差曲线公式中,可获得对应于预设代表风化层速度v0(800米/秒)的剩余风化层速度dv0为±300米/秒。在以下公式中所示的v1为折射层速度、vr为替换速度,均可通过设定获取;
其中,v1为折射层速度、vr为替换速度;
4)、根据所述预设代表风化层速度v0及所述剩余风化层速度dv0对应获取实际风化层速度变化区间。如通过1)中当所述预设代表风化层速度v0为800米/秒,所述剩余风化层速度dv0为±300米/秒时,实际风化层速度变化区间为500米/秒~1100米/秒。
根据以上步骤S102,首先,可通过设定不同的剩余静校正值δt及截距时t的边界范围值,对误差指数k进行设定,之后,在确定了误差指数k的情况下,结合预设的代表风化层速度值v0,获得与其对应的不同实际风化层速度变化区间,如:预设代表风化层速度值v0为800米/秒、750米/秒、850米/秒,或750米/秒~850米/秒中间的任意一个值时,将获得不同的实际风化层速度变化区间,而获得的每个实际风化层速度变化区间与其预设的代表风化层速度值v0一一对应,如:当实际风化层速度变化区间为500米/秒~1100米/秒时,对应的预设代表风化层速度值v0为800米/秒;当实际风化层速度变化区间为450米/秒~1050米/秒时,对应的预设代表风化层速度值v0为750米/秒;当实际风化层速度变化区间为550米/秒~1150米/秒时,对应的预设代表风化层速度值v0为850米/秒。
步骤S103:确定为处理中使用的代表风化层速度。判断所述获得的实际风化层速度变化区间是否介于临界风化层速度区间中,若是,则将所述预设代表风化层速度v0确定为代表风化层速度,若否,则重新预设代表风化层速度v0,直到可获得代表风化层速度。如,当步骤S102中,由预设代表风化层速度为800米/秒获得的实际风化层速度变化区间为500米/秒~1100米/秒时,判断是否介于临界风化层速度区间300~1500米/秒(其中300米/秒是声波在空气中的传播速度,1500米/秒是声波在水中的传播速度)区间内,由于上述获得的实际风化层速度变化区间在临界风化层速度区间中,因此,确定预设代表风化层速度值v0为800米/秒为处理中使用的代表风化层速度。
需要说明的是,通过以上实施例中的方法可在临界风化层速度区间300~1500米/秒中,通过设立多个预设代表风化层速度v0,对应获得多个与之对应的多个区间,并可根据实际需要通过调整误差指数k的方式,对区间宽度进行调整,同时,可将获得值进行验证,因此在步骤S103步骤还可包括:当所述判断所述实际风化层速度变化区间是否介于临界风化层速度区间中,若是,则将所述预设代表风化层速度v0确定为处理中所用的代表风化层速度后;根据所述处理中所用的代表风化层速度进行静校正处理,判断所述剩余静校正值δt是否在剩余静校正值边界范围值取值内,若是,则保留所述处理中所用的代表风化层速度,若否,则重新预设代表风化层速度v0,对获取的处理中所用的代表风化层速度在进一步给予验证,从而获得满足要求的代表风化层速度。
为进一步说明实施例步骤S102中所述误差曲线公式的获得,以下将结合附图给予详细说明:
图2为本发明的一个近地表模型示意图。其中D为基准面高程,S为地表面高程,R为假定风化层正确时对应的高速顶界高程,R’为假定风化层错误时对应的高速顶界高程,h为风化层厚度,dh为风化层厚度误差。
根据折射理论,折射层截距时t与风化层厚度h间的关系由下式表示:
其中,v0是风化层速度,v1是折射层速度。
由任意物理点的实际风化层速度得到对应的风化层厚度为:
假设dv0是某物理点实际风化层速度与选择的代表风化层速度之差,由速度误差引起的风化层厚度误差为:
其中
在常规处理的基准面校正处理中,将数据用风化层速度剥离到高速顶界,然后用替换速度充填到统一基准面。
如果风化层速度是正确的,炮点或检波点静校正值将是:
假定风化层速度不正确,得到的错误炮点或检波点静校正值将是:
由不正确的风化层厚度引起的静校正差为:
把方程式(2)和(3)代入方程(6)得到误差曲线公式为:
其中, 为误差指数。
具体的验证过程如下,在代表风化层速度v0为800米/秒的情况下,采用不同的折射层速度和替换速度组合,得到图3里对应的不同误差指数曲线。图中,横坐标是实际风化层速度与处理中采用的代表速风化层速度之差。尽管静校正误差随风化层速度误差的变化而变化,并且速度误差越大,静校正误差也越大,但误差指数k仍然可以用来指导代表风化层速度的选择。在实际生产中,给定两个边界条件:其一,由风化层速度误差引起的剩余静校正量不大于20毫秒,这种静校正量在后续的剩余静校正处理中可以很容易被消除掉;其二,高速折射层对应的截距时不超过400毫秒,也就是说,风化层厚度不超过300米。根据(7)式知道,这时对应k值约等于±0.1。从图3很容易看出,k=±0.1对应的实际风化层速度变化范围在-300米/秒到300米/秒。这就是说,对于本研究工区,如果选800米/秒代表风化层速度作为折射波静校正求解时的风化层初始速度,当实际风化层速度在500米/秒到1100米/秒间变化时,其引起的剩余静校正误差的绝对值都在20毫秒以内,完全在后续的剩余静校正程序控制范围以内。
实际上,这一结论不仅适用于本研究工区,仍然适用于中国广大地区乃至世界各地。因为世界上任何一个地区,风化层的速度均大于300米/秒(声波在空气中的速度)且小于1500米/秒(声波在水中的速度)。另外,本研究结果适用的风化层厚度达到300米的巨厚沙漠,如果风化层厚度没有这么厚,k值的选择还可以再大一些,因此风化层速度的允许变化范围还可以再宽一些,这就是说,对于一个缺乏近地表风化层资料的地区,要想利用折射波法解决静校正问题,选800米/秒作为代表风化层速度是一个科学的选择。
根据以上选定风化层速度v0的获取方法的静校正的方法,具体包括以下步骤:
1)、根据地震炮集记录拾取各物理点初至波时间;
2)、根据如实施例一中所述方法获取选定风化层速度v0;
3)、根据所述物理点初至波时间获取截距时t、折射层速度v1,并通过如下公式获得风化层厚度h;
其中,为物理点对应的折射波入射角;
在以上步骤的实际实现过程中可以将拾取的初至波时间先写入数据道头,之后根据折射层速度v1、获取的风化层速度v0;再通过以上计算风化层厚度h的公式获取风化层厚度h值。
4)、根据所述风化层厚度h获取所述物理点的炮点和检波点的静校正值;
5)、对上述4)中的输出数据再用处理流里的其他程序通过互相关技术做剩余静校正,消除各物理点由于采用的风化层速度和实际值不吻合而残留的剩余静校正值(静校正的高频成分)。
图4是该工区利用上述研究结果得到的近地表模型。图中指示线41为实际地表高程曲线清楚的反映了该区高程变化,从左边的接***坦到右边的巨厚沙丘。相对厚度变化达到200米以上。在左边戈壁部分,计算到的地下折射面和地表形态不一致,但地质上是合理的,它展示了一条地下古河道的存在。尽管右边沙漠地段地表高程变化非常崎岖,折射面高程(指示线42所示)变化要圆滑得多。为了证实计算出的近地表模型是否和实际炮记录一致,选择了一些用2500米/秒速度做过线性校正后的近排列代表性单炮记录,然后把其所在的位置和近地表模型联系起来,炮集按位置顺序重新排队,为的是更清楚的看到风化层和降速层是怎样随测线横向变化而变化的,如图5所示。图6为原来处理的剖面和新处理过的剖面对比图。原剖面上由静校正引起的反射层不连续现象在重新处理的新剖面上不再显现,黑圈内的构造成像有了较大改进。
本发明实施例的技术方案具有以下优点,在进行静校正分解之前,将无法预先知道的实际风化层速度变化拆分成两个部分:一部分是全区不变的代表风化层速度;另一部分是各物理点实际风化层速度和代表风化层速度之差。设法找到该代表风化层速度的取值范围,使剩余风化层速度引起的静校正量在后续处理中能被控制。
通过以上的实施方式的描述,本领域的技术人员可以清楚地了解到本发明可借助软件加必需的通用硬件平台的方式来实现,当然也可以通过硬件,但很多情况下前者是更佳的实施方式。基于这样的理解,本发明的技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来,该计算机软件产品存储在一个存储介质中,包括若干指令用以使得一台终端设备(可以是手机,个人计算机,服务器,或者网络设备等)执行本发明各个实施例所述的方法。
以上所述仅是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以做出若干改进和润饰,这些改进和润饰也应视本发明的保护范围。