一种速度模型的平滑方法及叠前深度偏移速度建模方法.pdf

(19)中华人民共和国国家知识产权局(12)发明专利申请(10)申请公布号CN107390265A(43)申请公布日2017.11.24(21)申请号201710598584.0(22)申请日2017.07.21(71)申请人中国石油集团川庆钻探工程有限公司地球物理勘探公司地址610213四川省成都市双流县华阳镇华阳大道一段一号1-4(72)发明人王光银陈三平巫骏熊晶璇李德珍范昆刘鸿陈燕辉(74)专利代理机构成都中玺知识产权代理有限公司51233代理人张敏谭昌驰(51)Int.Cl.G01V1/28(2006.01)G01V1/34(2006.01)权利要求书1页说明书4页附图2页(54)发明名称一种速度模型的平滑方法及叠前深度偏移速度建模方法(57)摘要本发明提供了一种速度模型的平滑方法及叠前深度偏移速度建模方法。所述平滑方法采用随深度变化的平滑参数实现对速度模型的空变平滑。所述叠前深度偏移速度建模方法在速度建模过程中，采用本发明的平滑方法对速度模型进行平滑。本发明能够实现了速度模型的空变平滑，同时保留浅部的层析细节和平滑深部的成块更新量。CN107390265ACN107390265A权利要求书1/1页1.一种速度模型的平滑方法，其特征在于，所述平滑方法采用随深度变化的平滑参数对速度模型进行平滑。2.根据权利要求1中任意一项所述的速度模型的平滑方法，其特征在于，所述平滑方法包括对速度模型X方向和/或Z方向上数据进行平滑，其中，所述对速度模型Z方向的数据进行平滑的步骤包括：沿速度模型Z方向进行分段，每一段分别采用各自的平滑参数完成速度模型Z方向上数据的平滑；所述对速度模型X方向的数据进行平滑的步骤包括：采用随深度方向变化的平滑参数，完成速度模型X方向上数据的平滑。3.根据权利要求2所述的速度模型的平滑方法，其特征在于，所述分段步骤中分段数D＝[log2NZ]+1，每一段长度为相邻两段之间的重合长度不小于其中，NZ为速度模型在Z方向上的采样点数目，Z方向为深度方向。4.根据权利要求2所述的速度模型的平滑方法，其特征在于，每一段采用阶数随深度变化的Savitzky-Golay滤波器，完成速度模型Z方向上数据的平滑。5.根据权利要求4所述的速度模型的平滑方法，其特征在于，所述每一段应用的Savitzky-Golay滤波器的阶数通过以下方式获得：根据下式计算每个深度点z的阶数nz，式中，NZ为速度模型在Z方向的采样点数目，Z方向为深度方向，MAX为最大平滑参数，MIN为最小平滑参数；求取每一段内深度点的平均阶数i＝1、2、…、D且为正整数；将所述平均阶数作为相应段内Savitzky-Golay滤波器的阶数。6.根据权利要求2所述的速度模型的平滑方法，其特征在于，沿X方向采用随深度变化的Savitzky-Golay滤波器，完成对速度模型X方向上数据的平滑。7.根据权利要求6所述的速度模型的平滑方法，其特征在于，根据下式计算X方向随深度变化的阶数nz，式中，NZ为速度模型在Z方向的采样点数目，Z方向为深度方向，MAX为最大平滑参数，MIN为最小平滑参数，z代表当前深度点；将所述阶数nz作为Savitzky-Golay滤波器的阶数。8.根据权利要求1中任意所述的速度模型的平滑方法，其特征在于，所述平滑方法采用最小二乘法对速度模型进行平滑，最小二乘法的阶数随深度变化。9.根据权利要求1中所述的速度模型的平滑方法，其特征在于，所述平滑方法采用阶数随深度变化的Savitzky-Golay滤波器对速度模型进行平滑。10.一种叠前深度偏移速度建模方法，其特征在于，所述叠前深度偏移速度建模方法在速度建模过程中，采用如权利要求1至9中任意一项所述的平滑方法对速度模型进行平滑。2CN107390265A说明书1/4页一种速度模型的平滑方法及叠前深度偏移速度建模方法技术领域[0001]本发明属于石油天然气地震勘探领域，更具体地讲，涉及一种用于石油地震勘探中深度域速度建模中对层析结果的平滑方法。背景技术[0002]复杂构造区域，叠前深度偏移有助于提高地下构造的成像效果。但是叠前深度偏移需要较为精确的速度场，现阶段，采用反射层析估计深度域速度。[0003]由于介质初始速度不准确，导致深部同向轴在道集中不聚焦，而此时道集内同向轴的弯曲度不仅受当前点速度不准确的影响同时也受上层地层速度不准确的综合影响。因此，随着深度的逐渐增加，道集内同向轴的弯曲度越不能真实反映当前反射点的速度偏差。[0004]随着深度的逐渐增加，深部的射线密度较低，层析结果容易在深部出现空白更新区，有必要通过平滑对层析结果进行平滑，使层析结果更能反映速度变化的趋势。然而，常规平滑方式不能同时保留浅部的层析细节和平滑深部的成块更新量。发明内容[0005]针对现有技术中存在的问题，本发明的目的在于解