样条曲线拟合初至波剩余静校正方法材料.docVIP

样条曲线拟合初至波剩余静校正方法 周强1,2，曹琳昱1，曹中林1,2，陈爱萍1,2 （1中国石油集团川庆钻探工程有限公司地球物理勘探公司，四川成都 610213 2 中国石油天然气集团公司山地地震技术试验基地，四川成都 610213） 摘要：针对复杂探区因模型道品质差导致反射波剩余静校正效果不好的难题，本文提出了一种基于样条曲线拟合的初至波剩余静校正方法。该方法相对直线拟合法，可以用的初至信息更多，统计的短波长静校正量更为准确；相对折线拟合法，它不用拾取折射拐点，操作方便快捷，更具有可操作性。它采用三次样条函数将应用了长波长静校正量的初至拟合成一条光滑曲线，根据地表一致性原则将延迟时分解为炮点、检波点短波长静校正量。实际资料应用结果表明，基于样条曲线拟合的单炮初至变得光滑，叠加剖面同相轴连续性得到增强，可为反射波剩余静校正提供较好品质的模型道，对低信噪比资料剩余静校正问题的解决具有实际意义。 关键字：静校正、剩余、短波长、初至波、样条曲线 收稿日期：2013.05.26 修回日期：2014.08.04 作者简介：周强（1985.04-）男，四川资中人，工程师，硕士研究生，主要从事复杂地区静校正方法研究.(E-mail: zhouq_wt@) 1 引言 静校正问题是制约我国复杂地区油气勘探发展的一个瓶颈问题，通常采用折射或层析方法，计算长波长静校正量，恢复地下构造形态，但是折射、层析方法求取的近地表结构，都是对速度、厚度的一个近似解，不可避免存在剩余静校正问题，影响成像效果[1-4]。根据处理对象不同，剩余静校正分为反射波、初至波剩余静校正。反射波剩余静校正应用最广、最常规、最可靠，但在复杂探区，由于资料信噪比过低而不能建立较好品质的模型道导致难以取得较好的效果[5-9]。物探工作者探索进行初至波剩余静校正，期望增强同相轴连续性，为反射波剩余静校正提供较好的模型道。 技术方法 研究表明，采用数学方法将初至拟合光滑的同时，可使反射波同相轴双曲线形态更加明显，使得能提高叠加剖面的成像精度[5-9]。初至波剩余静校正应该在在进行反射波剩余静校正之前做，它的一个关键是将离散的初至时间拟合成一条趋势线，有学者应用最小二乘法等数学方法将初至时间拟合成一条直线，这就需要根据偏移距对初至时间进行限制，那么可用的初至时间数量会减少，统计信息的可靠性会降低；在此基础上，有学者提出将初至按偏移距分段拟合成一条折线，统计信息有所增加，但不足之处在于计算机无法自动识别直达波、折射波、反射波，无法判断哪些初至时间来自同一层而影响计算结果。样条曲线在数学建模领域已被广泛使用，不仅能离散点拟合光滑，而且也能保持原始点的趋势变化，固将其引入地震勘资料处理领域，用于拟合初至时间。 1.1 样条曲线拟合初至时间 拟合的具体实施方法如下： 将坐标点及初至时间按如下表示 其中是应用了长波长静校正量后的初至时间， 为偏移距，为初至个数。 要拟合函数得到拟合后的函数须满足 （1） 且带有约束 （2） 其中和为预先提供的参数。为控制拟合函数光滑度的参数；为控制拟合函数偏差度的参数，即越大，允许拟合函数与原节点的偏差度越大。一般取范围为 将上述带约束的求极小问题转化为求解泛函极小问题[6]，即 （3） 其中为拉格朗日系数，是为了将（2）式化为等式所添加的实数型参数，（3）式以三次样条为目标函数，假设具有以下形式 （4） 根据的定义，其系数具有关系式： 其中， 。 计算中实际上只有2个未知量参与计算。 将带入（3）式，要求解（3）的泛函最小所对应的即要求解 （5） 其中 T是n-1阶的三对角矩阵，其对角元为，次对角元为。 Q是n+1行n-1列的三对角矩阵，其矩阵元素为 二维数据分正负偏移距建立方程组，三维分四象限建立方程组，利用牛顿迭代法求解方程组（5），即可获得拟合函数，将样点代入函数即得到拟合后的初至时间。 1.2 分解延迟时 本文方法在炮域、检波点域分别分解时差，得到炮点、检波点静校正量[10]。 延迟时是拟合前后的初至时间之差，表达式如下 （6） 其中，、分别为某一炮中偏移距为的这一道拟合前后的初至时间，为该道的延迟时间，其值等于检波点j的剩余静校量与炮点k的剩余静校量之和。 如果求取炮点桩号为k的炮点静校量，首先提取所对应的共炮点道集，假设有m道，就对应m个延迟时间，可建立含有m个方程的方程组，将其整理后