抛物线Radon变换消除多次波应用要素分析.pdfVIP

分衄交流论文·解释篇 抛物线Radon变换消除多次波的应用要素分析 张军华王月英郭见乐 (石油大学华东资源信息学院 山东东营 257061) (中石化胜利油田有限公司物探研究院 山东东营 257022) 1．引言 南于在抛物线拉东域，多次波和一次波的分离程度比t—x、，一k和f—P域都要好，所以抛物线 Radon变换已成为最好的去多次波方法之一。但即便如此，国内实际应用还较少。究其原因，一个重要的 因素就是对PRT的物理意义、适用条件、算法和参数的选择等关键问题缺乏深入的研究。国内外大多数文 献还有一个共同的缺憾，就是只见公式和结果，中间关键环节讨论很少，难见“庐山真面目”，所以也制 约了此项技术的推广。 2．时间域PRT的基本公式和物理意义 2．1基本公式 对于二维剖面u阮一，其时间域抛物线拉东正变换定义为 v(q，f)=∑“(x，t=f+掣2) (1) 式中口是以零偏移距为中心的抛物线的曲率，r为截距时间，x为偏移距。抛物线拉东变换也被叫做 f—g变换。 2．2物理意义 抛物线抗东变换是曲率为玑截距为f的抛物线上的所有点，聚焦到一个点上。由于CSP或CMP道集 地震同相轴经过动校正或，2一拉伸后，多次波同相轴近似为抛物线形态，所以，所以可用抛物线拉东变换 来去除多次波。 2．3时间域PRT的特点与存在的问题 模型已做过动校正处理．有三组一次波。两组多次波，第一组反射系数为负，对应波谷，曲率 1 C／2———_ 1 2．4Ⅸ2’ 第二组曲率取半。对图】做时间域拉东变换，得到图2所示的拉东正变换结果。从图2可以清楚地看到一 次波集中在小道号，即零偏移距附近；而多次波集中在较大的q道。在拉东域我们还可以看到拉东变换共 同的典型特征：在聚焦点附近有蘸个呈“剪刀”叉似的尾巴，水平状的同相轴对应近偏移距能量发散，倾 斜的同相轴代表远偏移距能量发散。 虽然时间域PRT具有物理意义明确、直观的特点，但时间域拉东变换还存在不可回避的缺点：1)拉 东域近偏移距的能量发散很大．图2中无论是一次波还是多次波水平方向都有长长的尾巴，当多个波存在 时会与其它能量干涉，这会给多次波的切除带来麻烦；2)同时域褶积和快速傅立叶变换类似，时间域PRT 分扭交流论文·解释篇 一定在离散样点上，这就必须插值，这会产生误差。 用1理论模型 图2时间域拉东变换结果 3．PRT的频率域实现方法 3．1反变换最小平方解法 对(1)式两边同时取傅氏正变换，并利用时移特陛，得 V(q，∞)=∑u(五o)e…2 (2) 对每一个频率分量，将上式写成矩阵形式： 矿：R西 (3) JOX9-x； …P…。：z 8 I ! ’． ! (4) R=l eJOT．qxf…P”矿i L J 其中，R。々=e’”‘“，i=1,2，·一，Kg；七=1,2，’一，，z。。 (3)式就是频率域拉东正变换公式。从(4)式可以看出，月是一复系数矩阵，由于％并不一定等 于。，，所以南驴就得不到矿的唯一解(矿、驴已做过多次波切除)。这时可借助最小平方原理，来得 到最优解。用矩阵形式描述： 驴：f尺“R+22i)。R”V’