一种构造约束的速度平滑方法和系统技术方案

本发明专利技术提供一种构造约束的速度平滑方法和系统。该方法包括：拾取构造层位；确定加权平滑因子；利用加权平滑因子对速度模型进行构造约束的平滑处理。本发明专利技术首先在地震剖面上拾取构造层位，然后沿着层位在一个矩形窗口内对速度数据加权平均滤波，通过纵向和横向的加权因子参数控制速度模型纵向和横向的平滑程度。由于构造层位的约束，经过平滑后的速度模型基本反映了地质构造的变化趋势。

【技术实现步骤摘要】
一种构造约束的速度平滑方法和系统
本专利技术涉及地球物理勘探领域，更具体地，涉及偏移速度建模研究领域，通过地质构造层位的约束，使平滑后的速度模型更接近地质构造形态，因而偏移成像的结果更合理。
技术介绍
当前，关于速度平滑的方法有很多，最常用的方法有：中值滤波方法、高斯平滑方法和滑动平均方法等，这些方法只考虑了数据本身的特性，采用统计的方法，去除数据中的噪音，没有考虑地质构造的约束。不管是叠前深度偏移还是叠前反演，都需要准确的速度模型，而且速度模型要反映地质构造的变化趋势。但在一些地震资料信噪比低的探区，或者是构造复杂的探区，由层析反演建立的速度模型往往不能反映地质构造的变化趋势，因而影响叠前深度偏移的效果。因此，本领域急需一种能够反映地质构造的变化趋势的速度平滑方法。
技术实现思路
本专利技术的目的是提出一种构造约束的速度平滑方法。本专利技术的构造约束的速度平滑方法是采用加权平均的方法，在一个矩形窗口内，根据加权因子，沿着构造层位进行平滑，使平滑后的速度模型具有地质构造的变化形态。本专利技术提出的构造约束的速度平滑方法，根据地质层位信息，沿着层位在一个矩形窗口内对数据加权平均滤波，即消除了速度的异常值，又使速度模型能够反映地质构造的形态。根据本专利技术的一个方面，提供一种构造约束的速度平滑方法，该方法包括：拾取构造层位；确定加权平滑因子；利用加权平滑因子对速度模型进行构造约束的平滑处理。进一步地，确定ir1、ir2和alfa为加权平滑因子，ir1和ir2为大于1的整数，alfa为大于零小于1的浮点数。进一步地，利用以下公式进行平滑处理：vnew(j,k)＝vave*alfa+(1-alfa)*ss，其中mz＝[(k-iup)*horsur(i+1,mj)+(idown-k)*horsur(i,mj)]/(idown-iup)ir＝ir1+(ir2-ir1)/n式中，vnew(j,k)表示平滑后的速度，vave为层内平均速度，n表示构造层位数，i＝1,2,3,…,n，第i层和i+1层在CDP＝j位置的深度为iup和idown，i层和i+1层之间深度为k位置的原始速度记为vold(j,k)，horsur(i,mj)为第i层在CDP＝mj位置的层位深度。进一步地，ir1和ir2用于控制横向平滑程度，alfa用于控制纵向平滑程度。进一步地，通过试验选取ir1、ir2和alfa合适的数值。根据本专利技术的另一个方面，提供一种构造约束的速度平滑系统，该系统包括：用于拾取构造层位的装置；用于确定加权平滑因子的装置；利用加权平滑因子对速度模型进行构造约束的平滑处理的装置。针对中国南方某探区低信噪比的地震资料，由层析成像反演得到的速度模型的构造形态与偏移剖面不一致的问题，专利技术人首先在偏移剖面上拾取构造层位，采用本专利技术的平滑方法，对层析反演的速度模型进行了处理，平滑后的速度模型的构造形态基本与偏移剖面一致，说明本方法是有效的。附图说明通过结合附图对本专利技术示例性实施方式进行更详细的描述，本专利技术的上述以及其它目的、特征和优势将变得更加明显，其中，在本专利技术示例性实施方式中，相同的参考标号通常代表相同部件。图1示出了本专利技术构造约束的速度平滑方法的流程图。图2示出了构造约束的速度平滑示意图。附图3示出了中国南方某探区的构造约束速度平滑计算实例，其中(a)示出了地震偏移剖面；(b)示出了由速度谱拾取得到的均方根速度转换到深度域的速度模型剖面，可以看出速度剖面没有反应出图(a)所示的构造形态；(c)是对图(b)所示的速度模型，由构造约束速度平滑得到的速度剖面，基本反映了地震剖面中的构造形态。具体实施方式下面将参照附图更详细地描述本专利技术的优选实施方式。虽然附图中显示了本专利技术的优选实施方式，然而应该理解，可以以各种形式实现本专利技术而不应被这里阐述的实施方式所限制。本专利技术针对构造复杂、低信噪比地震资料速度建模精度低，速度模型往往不能反映地质构造的变化趋势等问题，提出构造约束的速度平滑技术。本专利技术提供一种构造约束的速度平滑方法，该方法包括：拾取构造层位；确定加权平滑因子；利用加权平滑因子对速度模型进行构造约束的平滑处理。具体地，首先在地震剖面上拾取构造层位，然后沿着层位在一个矩形窗口内对速度数据加权平均滤波，通过纵向和横向的加权因子参数控制速度模型纵向和横向的平滑程度。由于构造层位的约束，经过平滑后的速度模型基本反映了地质构造的变化趋势。优选地，确定ir1、ir2和alfa为加权平滑因子，ir1和ir2为大于1的整数，alfa为大于零小于1的浮点数。ir1和ir2用于控制横向平滑程度，alfa用于控制纵向平滑程度。优选地，利用以下公式进行平滑处理：vnew(j,k)＝vave*alfa+(1-alfa)*ss，其中mz＝[(k-iup)*horsur(i+1,mj)+(idown-k)*horsur(i,mj)]/(idown-iup)ir＝ir1+(ir2-ir1)/n式中，vnew(j,k)表示平滑后的速度，vave为层内平均速度，n表示构造层位数，i＝1,2,3,…,n，第i层和i+1层在CDP＝j位置的深度为iup和idown，i层和i+1层之间深度为k位置的原始速度记为vold(j,k)，horsur(i,mj)为第i层在CDP＝mj位置的层位深度。优选地，通过试验选取ir1、ir2和alfa合适的数值。在实际资料的处理中，需要先选择小范围数据做实验处理，确定合适的加权因子参数，再进行大数据体的处理。本专利技术在中国南方某探区取得了较好的效果。此外，本专利技术还提供一种构造约束的速度平滑系统，其特征在于，该系统包括：用于拾取构造层位的装置；用于确定加权平滑因子的装置；利用加权平滑因子对速度模型进行构造约束的平滑处理的装置。接下来参照附图详细描述本专利技术的实施例。如图1所示，在本实施例中，构造约束的速度平滑方法包括：拾取构造层位；确定加权平滑因子；利用加权平滑因子对速度模型进行构造约束的平滑处理；最后输出结果。接下来参照图2和图3描述本实施例的具体过程。如图2所示，假设有n个构造层位，i＝1,2,3,…,n，第i层和i+1层在CDP＝j位置的深度为iup和idown，i层和i+1层之间深度为k位置的原始速度记为vold(j,k)，平滑后的速度为vnew(j,k)，则：vnew(j,k)＝vave*alfa+(1-alfa)*ss(1)mz＝[(k-iup)*horsur(i+1,mj)+(idown-k)*horsur(i,mj)]/(idown-iup)ir＝ir1+(ir2-ir1)/nir1(>＝2)、ir2(>＝2)和alfa(0-1)为控制参数，称为加权平滑因子，ir1和ir2控制横向平滑程度，数值越大越平滑；alfa控制纵向平滑程度，数值越大构造特征越明显。vave为层内平均速度，horsur(i,mj)为第i层在CDP＝mj位置的层位深度(采样点数)。首先，拾取构造层位。在偏移剖面上，根据构造特征，拾取具有代表性的构造层位，主要是具有地质意义的层位，如：岩性分界面、风化剥蚀面等。接下来，确定参数ir1、ir2、alfa。ir1和ir2为大于1的整型数，alfa为大于零小于1的浮点数。具体地，通过试验选取合适的数值，即选本文档来自技高网...

【技术保护点】
1.一种构造约束的速度平滑方法，其特征在于，该方法包括：拾取构造层位；确定加权平滑因子；利用加权平滑因子对速度模型进行构造约束的平滑处理。

【技术特征摘要】
1.一种构造约束的速度平滑方法，其特征在于，该方法包括：拾取构造层位；确定加权平滑因子；利用加权平滑因子对速度模型进行构造约束的平滑处理。2.如权利要求1所述的构造约束的速度平滑方法，其特征在于，确定ir1、ir2和alfa为加权平滑因子，ir1和ir2为大于1的整数，alfa为大于零小于1的浮点数。3.如权利要求2所述的构造约束的速度平滑方法，其特征在于，利用以下公式进行平滑处理：vnew(j,k)＝vave*alfa+(1-alfa)*ss，其中mz＝[(k-iup)*horsur(i+1,mj)+(idown-k)*horsur(i,mj)]/(idown-iup)ir＝ir1+(ir2-ir1)/n式中，vnew(j,k)表示平滑后的速度，vave为层内平均速度，n表示构造层位数，i＝1,2,3,…,n，第i层和i+1层在CDP＝j位置的深度为iup和idown，i层和i+1层之间深度为k位置的原始速度记为vold(j,k)，horsur(i,mj)为第i层在CDP＝mj位置的层位深度。4.如权利要求2所述的构造约束的速度平滑方法，其特征在于，ir1和ir2用于控制横向平滑程度，alfa用于控制纵向平滑程度。5.如权利要求2所述的构造约束的速度平滑方法，其特征在于，通过试验选取ir1、ir2和alfa合适的数值。6.一种构造...

【专利技术属性】
技术研发人员：周巍，胡光辉，张兵，
申请(专利权)人：中国石油化工股份有限公司，中国石油化工股份有限公司石油物探技术研究院，
类型：发明
国别省市：北京,11