文章编号：1001—1749(2008)01—0017—06地震剖面匹配闭合差校正高荣涛，张军华，周振晓，钟磊，田连玉(中国石油大学地球资源与信息学院，山东东营257061)摘要：当用不同施工时期采集的地震剖面对同一地区进行综合对比解释时，可能会遇到相交剖面产生闭合差情况。通常这些闭合差是由于激发因素、接收条件的不同，导致交叉点上振幅、时间及相位存在差异造成的，这将给地震资料的综合解释带来不便。因此，随着四维地震技术的推广，做好地震剖面的匹配校正是一件很有意义的工作。这里给出了2D地震测线上各个交点的振幅、时间和相位差匹配校正的原理公式，利用迭代最小平方算法，通过模型试算，详细阐述匹配流程并进行了效果分析。关键词：闭合差；校正；振幅；相位旋转；时移;中图分类号：P631.4文献标识码：A0前言在分析测线网络中的数据剖面时，常会看到交叉点处剖面之间的地震道不匹配，通常这种交叉点的不匹配称为数据闭合差。这些闭合差有些在视觉上是可以捕捉到的如明显整体幅值差异和纵向时差。它的存在给我们的解释评价工作造成了很大的不便，特别是随着勘探难度的增加和四维地震技术的发展，这要求我们利用一切可以获得的资料信息，进行综合储层预测，而要想有效地进行综合解释，必须做好新、老资料之间的匹配整合。通过匹配整合使老资料具有利用价值,这比重新勘探在经济上具有绝对的优势，同时避免了风险投资,因此这种整合处理在叠后十分有意义。一般方法是针对产生闭合差因素的各个影响值，采取相应的改正方法逐一校正。但目前的研究趋势与此相反都趋向于在可靠的理论分析和因素分析基础上,对地震道的差值结果进行一次性校正，将基于地震道的校正结果作为对所有野外采集、处理参数等因素的修正,其剩余校正值作为偏移剖面固有的差异。基于Bishop⑴提出的利用最小平方校正来消除闭合差的思路和方法，作者对模型数据和实际地震资料，进行了匹配整合处理,达到了不错的效果。该方法简单易行，同时能够解决测线网络中多道交叉匹配校正的情况。1迭代匹配法原理假设在某个测线网络中有m条测线，n个交点，由)和秫)分别是第i个交点处第k、l条测线的某个标度(时间、振幅和相位)，通常令1(订＞k(i)，我们定义b(i)为第i个交点处某个标度的闭合差，即b(i)=xi(i)-xl(i).在最小平方约束下，对于剖面j(j=l,2,…，m)，假设其第「次迭代过程产生的估算值为可,这里「=0、1、…表示迭代次数。在每次迭代中，对估算校正值⑵进行校正是通过加进该线的平均观测闭合差与平均计算闭合差的误差来实现的(在这里适当地考虑其符号)。在一个交叉点处计算的闭合差，是从各相交线的当前估计校正最中产生的所以,对于第"1次迭代估计值可“与第「次迭代估计值之间，存在下面的关系W"=可+土£（气［6-（%；（,）-x；（0）］）令xj=0（1）式中•1,若k（i）=j.ay=•-1»若心）=j»%=I。，其它为统计出测线，上的交点数。该过程中用观测值和计算值的均方根误差来监控该方法的收敛性。£=J+£［6-（x*（i）-%心））］2（2）收敛准则为心菩W&3是一个很小的数（通常e取为0.00001）。这种方法可以简单有效地计算出闭合差校正值，且Bishop等已经证明了该法所得结果与最小平方解等价。因此，对于交点个数大于剖面条数的超定问题就简化为，只要知道所有相交剖面之间的某个标度的闭合差距阵B=|6Ifb29,・•,bit…，以及闭合差所对应的相交剖面的序号，就可以方便地计算出每张剖面需要校正的某个标度的最佳值。交点塞金项目：中国石油大学研究生创新项目(S200624)收稿日期：2006-12-21改回日期：2007-06-192闭合差的计算前面给出了利用闭合差来计算校正值,现在来讨论求取闭合差的原理公式。通常闭合差有振幅、时差和相位三个标度，在实际的地震资料中这三个标度是结合在一起的需要我们分开来计算闭合差。2.1时间闭合差的求取为了求取相交地震剖面k,的时间闭合差△》,在交点位置一定时窗长度内取地震数据。时窗长度可以不相同，如先在剖面，上开时窗十,则在剖面k上以a所在位置为中心，上、下各多取100个点，此时时窗Tk=T,+200o然后利用互相关计算公式（3）,做滑动互相关处理⑶。互相关公式为j=成（3）项2（力项2其中为、"表'示...