专利名称：一种优选风化层速度的确定方法
技术领域：
本发明涉及石油勘探领域，应用于地震数据处理过程中，特别是涉及ー种优选风化层速度的确定方法。
背景技术：
现有技术的地震勘探解释的理论中，一般都假定激发点与接收点是在ー个水平面上，并且地层速度是均匀的。但实际上地面常常不平坦，各个激发点深度也可能不同，低速带中的波速与地层中的波速又相差悬殊，所以必将影响实测的时距曲线形状。为了消除这些影响，对原始地震数据要进行地形校正、激发深度校正、低速带校正等。由于这些校正对 同一观测点的不同地震界面都是不变的，因此统称静校正。而在静校正过程中，特别是长波长静校正问题是复杂近地表地区地震资料处理中提高资料处理质量、保证获得地下真实构造形态需要解决的主要问题。目前在地震资料处理领域中有很多关于折射波静校正的方法，如GLI，EGRM，加减法，TOMO等等。现有技术中折射波静校正的方法实现过程中涉及风化层速度填写方法通常为以下两种方式第一种方式是当研究エ区内能预先提供局部控制点风化层速度信息时，把这些信息分别填入程序各控制点位置对应的风化层速度參数栏。反演程序执行时，其他点的速度參数通过控制点信息插值获得；第二种方式是当研究エ区内不能预先提供任何风化层速度信息时。在程序的风化层速度參数栏填入任意ー个常数，用于全区反演运算，从而由于任意常数选择的盲目性，很容易造成选择的速度和实际速度有较大偏差，由此引起较大的模型反演错误和静校正错误，超过了后续剩余静校正程序所能修正的限制，不能获得理想的静校正結果。综上所述，发明人在实现本发明过程中发现现有技术中存在以下缺陷预获得可靠的折射波静校正結果，必须预先掌握近地表地层风化层速度的横向变化信息，因此需要预先进行诸如微测井等的近地表勘測，从而获取近地表真实数据，在以上的勘测过程中需要巨大的成本投入，同时由于其勘测点的分布及数据限制，其对不同勘测点的选择与数据处理，将直接影响到后期数据的精确度。同时，盲目的任意常数选择容易造成选择的速度和实际速度有较大偏差，最终无法达到预期效果。

发明内容
本发明实施例要解决的问题是提供一种优选风化层速度的确定方法，解决了因复杂地区近地表信息无法预先精确确定，由此无法准获取风化层速度，从而造成后续静校正精度偏低的问题。为达到上述目的，本发明实施例提出一种优选风化层速度的确定方法，具体包括预设代表风化层速度V0 ;根据所述预设代表风化层速度Vtl获取实际风化层速度变化区间；判断所述实际风化层速度变化区间是否介于临界风化层速度区间中，若是，则将所述预设代表风化层速度Vtl为处理中所用的代表风化层速度，若否，则重新预设代表风化层速度Vf!。根据本发明的实施例以上技术方案获得了以下技术效果在折射波静校正过程中，不再需要预先进行诸如微测井等的近地表调查，就可以使复杂地区为了获得地下勘探目的层精确构造成像而必需进行的短、长波长折射波静校正在精度方面得到提高。


为了更清楚地说明本发明实施例的技术方案，下面将对实施例描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。图I为本发明实施例中优选风化层速度的确定方法的流程示意图；图2为本发明的近地表模型示意图；图3为本发明的采用代表风化层速度Vtl为800米/秒、不同的折射层速度和替换速度组合得到的误差曲线。图4为本发明实施例中实际得到的近地表模型。图5为本发明实施例中的近地表模型和实际炮记录的吻合度检查图。图6为本发明实施例中原来处理的剖面和新处理过的剖面对比图。
具体实施例方式本发明实施例提供一种优选风化层速度的确定方法，为满足折射波静校正过程中校正精度，通过对剩余静校正值8 t数据及截距时的设定，通过误差曲线公式从而确定其风化层速度Vo值。下面结合附图和实施例，对本发明的具体实施方式
作进ー步详细描述实施例如图I所示，为本发明一种优选风化层速度的确定方法的流程示意图，包括以下步骤步骤SlOl :预设代表风化层速度Vtl ;所述预设代表风化层速度Vtl可选取750米/秒 850米/秒之间的任意一个值，例如800米/秒。步骤S102 :根据所述预设代表风化层速度Vtl获取实际风化层速度变化区间，具体过程如下；
I)、根据剩余静校正值S t及截距时t的边界范围值，确定所述剩余静校正值5 t、截距时t，在实际的地质环境下，为使剩余静校正量在后续的剩余静校正处理中可以被消除掉，通常剩余静校正值S t具体边界范围绝对值不超过20毫秒；同时对截距时t的取值范围，根据我国地质地貌的实际情况，可具体选取为截距时t的边界范围值不超过400毫秒，即风化层厚度不超过300米，因此，所述剩余静校正值S t的绝对值为20毫秒，所述截距时t的取值为400毫秒。2)、将所述剩余静校正值S t及截距时t代入如下误差曲线公式中获取误差指数^St —当所述剩余静校正值 S t为±20毫秒，及截距时t为400毫秒时，所述误差指数k等于±0. I。3)、将所述误差指数k及分别将所述预设代表风化层速度Vtl代入如下误差曲线公式中分别获取剰余风化层速度dV(l ;如将2)中获取的k等于±0. 1，及在3)中通过预设选取的代表风化层速度预设值Vtl为800米/秒，代入如下误差曲线公式中，可获得对应于预设代表风化层速度Vtl (800米/秒)的剰余风化层速度(1%为±300米/秒。在以下公式中所示的V1为折射层速度、Vr为替换速度，均可通过设定获取；k =I ノ(^77 "ノ + vrtos° ini = arcsin ^in = arcsin v°+dv°其中，V1为折射层速度、Vr为替换速度；4)、根据所述预设代表风化层速度Vtl及所述剩余风化层速度dV(l对应获取实际风化层速度变化区间。如通过I)中当所述预设代表风化层速度Vtl为800米/秒，所述剩余风化层速度dV(l为±300米/秒吋，实际风化层速度变化区间为500米/秒 1100米/秒。根据以上步骤S102，首先，可通过设定不同的剩余静校正值St及截距时t的边界范围值，对误差指数k进行设定，之后，在确定了误差指数k的情况下，结合预设的代表风化层速度值Vtl,获得与其对应的不同实际风化层速度变化区间，如预设代表风化层速度值Vtl为800米/秒、750米/秒、850米/秒，或750米/秒 850米/秒中间的任意ー个值时，将获得不同的实际风化层速度变化区间，而获得的每个实际风化层速度变化区间与其预设的代表风化层速度值Vtl—一对应，如当实际风化层速度变化区间为500米/秒 1100米/秒时，对应的预设代表风化层速度值Vtl为800米/秒；当实际风化层速度变化区间为450米/秒 1050米/秒时，对应的预设代表风化层速度值Vtl为750米/秒；当实际风化层速度变化区间为550米/秒 1150米/秒时，对应的预设代表风化层速度值Vtl为850米/秒。步骤S103 :确定为处理中使用的代表风化层速度。判断所述获得的实际风化层速度变化区间是否介于临界风化层速度区间中，若是，则将所述预设代表风化层速度Vtl确定为代表风化层速度，若否，则重新预设代表风化层速度Vtl，直到可获得代表风化层速度。如，当步骤S102中，由预设代表风化层速度为800米/秒获得的实际风化层速度变化区间为500米/秒 1100米/秒时，判断是否介于临界风化层速度区间300 1500米/秒(其中300米/秒是声波在空气中的传播速度，1500米/秒是声波在水中的传播速度)区间内，由于上述获得的实际风化层速度变化区间在临界风化层速度区间中，因此，确定预设代表风化层速度值\为800米/秒为处理中使用的代表风化层速度。需要说明的是，通过以上实施例中的方法可在临界风化层速度区间300 1500米/秒中，通过设立多个预设代表风化层速度Vtl，对应获得多个与之对应的多个区间，并可根据实际需要通过调整误差指数k的方式，对区间宽度进行调整，同时，可将获得值进行验证，因此在步骤S103步骤还可包括当所述判断所述实际风化层速度变化区间是否介于临界风化层速度区间中，若是，则将所述预设代表风化层速度\确定为处理中所用的代表风化层速度后；根据所述处理中所用的代表风化层速度进行静校正处理，判断所述剩余静校正值St是否在剩余静校正值边界范围值取值内，若是，则保留所述处理中所用的代表风化层速度，若否，则重新预设代表风化层速度Vtl，对获取的处理中所用的代表风化层速度在进ー步给予验证，从而获得满足要求的代表风化层速度。为进一歩说明实施例步骤S102中所述误差曲线公式的获得，以下将结合附图给予详细说明图2为本发明的ー个近地表模型示意图。其中D为基准面高程，S为地表面高程，R为假定风化层正确时对应的高速顶界高程，R’为假定风化层错误时对应的高速顶界高程，h为风化层厚度，Ah为风化层厚度误差。根据折射理论，折射层截距时t与风化层厚度h间的关系由下式表示
权利要求
1.一种优选风化层速度的确定方法，其特征在于，包括 预设代表风化层速度Vtl ； 根据所述预设代表风化层速度Vtl获取实际风化层速度变化区间； 判断所述实际风化层速度变化区间是否介于临界风化层速度区间中，若是，则将所述预设代表风化层速度Vtl确定为处理中所用的代表风化层速度，若否，则重新预设代表风化层速度Vf!。
2.如权利要求I所述的方法，其特征在于，所述根据所述预设代表风化层速度Vtl获取实际风化层速度变化区间步骤还包括 根据剩余静校正值S t及截距时t的边界范围值，确定所述剩余静校正值S t、截距时t ； 将所述剩余静校正值S t及截距时t代入如下误差曲线公式中获取误差指数k ；
3.如权利要求I或2所述的方法，其特征在干，所述临界风化层速度区间为300米/秒 1500米/秒。
4.如权利要求I或2所述的方法,其特征在于,所述剩余静校正值St的边界范围绝对值不超过20毫秒，所述截距时t的边界范围值不超过400毫秒。
5.如权利要求4所述的方法，其特征在于，所述将所述剩余静校正值St及截距时t代入如下误差曲线公式中获取误差指数k步骤还包括 当所述剩余静校正值S t为±20毫秒，及截距时t为400毫秒时，获取所述误差指数k等于±0. I。
6.如权利要求I中所述的方法，其特征在于，所述预设代表风化层速度Vtl包括800米/秒、750米/秒 850米/秒。
7.如权利要求6所述的方法，其特征在于，所述将所述误差指数k及所述预设代表风化层速度Vtl代入如下误差曲线公式中分别获取剰余风化层速度dV(l步骤还包括 当所述误差指数k为±0. I、所述预设代表风化层速度Vtl为800米/秒吋，获取剩余风化层速度dvQ为±300米/秒。
8.如权利要求7所述的方法，其特征在于，所述根据所述预设代表风化层速度Vtl及所述剩余风化层速度dV(l获取实际风化层速度变化区间步骤还包括 当所述预设代表风化层速度Vtl为800米/秒，所述剩余风化层速度dV(l为±300米/秒吋，获取所述实际风化层速度变化区间为500米/秒 1100米/秒；或当所述预设代表风化层速度\为750米/秒，所述剩余风化层速度dV(l为±300米/秒吋，获取所述实际风化层速度变化区间为450米/秒 1050米/秒； 或当所述预设代表风化层速度Vtl为850米/秒，所述剩余风化层速度dV(l为±300米/秒吋，获取所述实际风化层速度变化区间为550米/秒 1150米/秒。
9.如权利要求I所述的方法，其特征在于，所述判断所述实际风化层速度变化区间是否介于临界风化层速度区间中，若是，则将所述预设代表风化层速度Vtl确定为处理中所用的代表风化层速度步骤还包括 判断所述实际风化层速度变化区间是否介于临界风化层速度区间中，若是，则将所述预设代表风化层速度Vtl确定为处理中所用的代表风化层速度； 根据所述处理中所用的代表风化层速度进行静校正处理，判断所述剩余静校正值S t是否在剩余静校正值边界范围值取值内，若是，则保留所述处理中所用的代表风化层速度，若否，则重新预设代表风化层速度K。
全文摘要
本发明实施例公开了一种优选风化层速度的确定方法，具体包括预设一个代表风化层速度值，该值在风化层厚度不小于300米情况下，与实际风化层速度差引起的静校正量完全可限制在后续剩余静校正程序允许的范围以内，近地表模型的建立就可以获得稳定解。因此，在折射波静校正过程中，不再需要预先进行诸如微测井等的近地表调查，就可使复杂地区构造成像而必需进行的短、长波长折射波静校正在精度方面得到提高。
文档编号G01V1/36GK102645671SQ20121010369
公开日2012年8月22日 申请日期2012年4月10日 优先权日2012年4月10日
发明者吴媚, 孙庚文, 欧阳行艳, 王顺根, 谢桂生 申请人:恒泰艾普石油天然气技术服务股份有限公司