CN107422382A - 自适应二维多项式拟合的滤波方法和装置 - Google Patents
自适应二维多项式拟合的滤波方法和装置 Download PDFInfo
- Publication number
- CN107422382A CN107422382A CN201610350006.0A CN201610350006A CN107422382A CN 107422382 A CN107422382 A CN 107422382A CN 201610350006 A CN201610350006 A CN 201610350006A CN 107422382 A CN107422382 A CN 107422382A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- picture point
- breakpoint
- window
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 48
- 238000001914 filtration Methods 0.000 title claims abstract description 14
- 230000005012 migration Effects 0.000 claims description 12
- 238000013508 migration Methods 0.000 claims description 12
- 238000005070 sampling Methods 0.000 claims description 9
- 238000003384 imaging method Methods 0.000 claims description 7
- 230000008569 process Effects 0.000 abstract description 10
- 238000003325 tomography Methods 0.000 abstract description 3
- 238000012545 processing Methods 0.000 description 13
- 238000010586 diagram Methods 0.000 description 12
- 230000000694 effects Effects 0.000 description 8
- 230000006870 function Effects 0.000 description 7
- 238000005516 engineering process Methods 0.000 description 5
- 238000004590 computer program Methods 0.000 description 4
- 230000009466 transformation Effects 0.000 description 4
- 230000009471 action Effects 0.000 description 3
- 241001269238 Data Species 0.000 description 2
- 238000013459 approach Methods 0.000 description 2
- 230000005540 biological transmission Effects 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 2
- RYGMFSIKBFXOCR-UHFFFAOYSA-N Copper Chemical compound [Cu] RYGMFSIKBFXOCR-UHFFFAOYSA-N 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 229910052802 copper Inorganic materials 0.000 description 1
- 239000010949 copper Substances 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 230000008030 elimination Effects 0.000 description 1
- 238000003379 elimination reaction Methods 0.000 description 1
- 230000002349 favourable effect Effects 0.000 description 1
- 239000000835 fiber Substances 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 239000013307 optical fiber Substances 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 239000004576 sand Substances 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 230000005945 translocation Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/36—Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
- G01V1/364—Seismic filtering
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/20—Trace signal pre-filtering to select, remove or transform specific events or signal components, i.e. trace-in/trace-out
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/30—Noise handling
- G01V2210/32—Noise reduction
- G01V2210/322—Trace stacking
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Image Processing (AREA)
Abstract
本发明公开了自适应二维多项式拟合的滤波方法和装置。本发明依据多项式拟合原理,对二维多项式拟合公式进行了离散化和具体化,此处的二维分别指叠前道集的偏移距和CDP方向,以更好地保持地震信号振幅的连续性。同时,本发明中沿局部地层倾角方向进行二维多项式拟合来消除随机噪音,能够实现自动处理断层边界,以保证断层断点的清晰度。
Description
技术领域
本发明涉及地震勘探中反射地震资料处理领域,更具体地,涉及一种自适应二维多项式拟合的滤波方法和一种自适应二维多项式拟合的滤波方法。
背景技术
在油气勘探中对小断层及薄砂体的识别,是认识油气疏导体系、识别有利储层的重要环节。因此,在地震勘探资料处理过程中,提高地震成像的分辨率一直是一个重要的努力目标。而由于地震数据在激发与接收过程存在干扰,并且在后续处理过程中也会引入人为因素的噪音,因此如何准确有效的去除噪音干扰影响一直以来都是高分辨率地震资料处理流程中的重要步骤。它不仅能够提高最终叠加剖面的信噪比和分辨率,而且对叠前道集的进一步应用,如岩性反演,AVO/AVA分析等起到至关重要的作用。
叠前地震数据往往由于进行了反Q叠前偏移成像,在补偿地层粘性吸收、拓宽地震数据频带和提高高频有效信息的同时引入了大量的高频噪音,因此如何准确有效地去除高频随机噪音就成为了一个难题。而在前人的滤波方法研究中,诸如常规的预测滤波及小波变换滤波等方法在去除噪音的同时均不可避免的对地震波的有效频带尤其是高频信息有一定的损伤(利用沃尔什变换提高地震资料信噪比,徐云霞等,岩性油气藏,2009;扩散滤波方法在地震资料处理中的应用研究,陈可洋等,岩性油气藏,2014)。
多项式拟合去噪技术在国内最早由俞寿朋等人(用地震信号多项式拟合提高叠加剖面信噪比,石油地球物理勘探,1988,23(2):131-139.)提出并取得了较好的应用效果,该方法基于叠前道集中同相轴沿偏移距和CDP方向上能够保持一定的连续性的假设,首先通过某数据点周围的局部离散点来拟合一个二维高次多项式(一般为二次或者三次),然后通过该多项式在该数据点处的函数值来取代原始数据以达到消除随机噪音的目的。之后也有很多研究者提出了多项式拟合的改进方法并取得了很好的效果,杨云飞等人(正交多项式逼近去白噪,石油地球物理勘探,2001,36(6):704-706.)介绍了一种多项式拟合去白噪方法,利用切比雪夫多项式或勒让德多项式拟合逼近数据列去除白噪声,并使用该方法对测井数据和野外地震数据进行了处理,取得的良好的效果。但是,由于在多项式阶数上还不好准确把握,有效信号与白噪之间还存在混叠问题,而且,实际资料中的白噪往往也并非完全意义上的随机干扰。这样一来,一方面阶数偏差会对有效信号造成一定损失;另一方面,混叠与失真也影响了对白噪的去除效果。钟伟等人(多项式拟合技术在强噪声地震资料中的应用研究,地球物理学进展,2006,21(1):184-189.)对多项式拟合去噪技术进行了研究,通过对实际数据的进行处理,结果表明多项式拟合去噪方法能够很好地消除随机噪声,其缺点是多项式拟合去噪方法中由于波形的估计对振幅有一定程度的改造(即振幅保真性差的问题),使处理后的剖面看起来有点“假”。李向云等人(改进的正交多项式变换在地震资料去噪中的应用,地球物理学进展,2013,28(4):2085-2093.)提出了一种基于SVD(奇异值分解)技术和小波变换的正交多项式联合去噪的方法,利用SVD技术估计的有效信号能量能自适应地确定各个时间多项式的阶数,有效地压制了高阶上的随机噪声,同时利用小波阀值处理能在一定程度上减弱系数谱上混叠影响。该方法适用于同相轴近似水平的地震资料,而且当地震资料的信噪比很低时,由于严重的低阶混叠现象,无法取得很好的去噪效果。
发明内容
本发明提出了一种能够去除叠前道集的噪音同时保持叠前道集中有效信号的相对振幅及频带的滤波方法。本发明还提出了相应的装置。
根据本发明的一方面,提出了一种自适应二维多项式拟合的滤波方法,该方法包括:得到叠前偏移距道集;将叠前偏移距道集整理为维度分别为三维数据体G(xi,offj,tl),xi表示像点的水平位置,offj表示偏移距长度,tl表示像点的垂直旅行时;在三维数据体G(xi,offj,tl)中,按照垂直旅行时对各个偏移距对应的数据进行叠加,得到二维叠加剖面P(xi,tl);基于二维叠加剖面P(xi,tl)得到成像空间内反射面的地层倾角并识别三维数据体G(xi,offj,tl)中的断点和非断点;针对每一被识别为非断点的像点,沿所述地层倾角方向设置多项式拟合窗口,进行二维多项式拟合以更新该像点。
根据本发明的另一方面,提出了一种自适应二维多项式拟合的滤波装置,该装置包括:偏移距道集获取单元,用于得到叠前偏移距道集;三维数据体获取单元,用于将叠前偏移距道集整理为维度分别为三维数据体G(xi,offj,tl),xi表示像点的水平位置,offj表示偏移距长度,tl表示像点的垂直旅行时;二维剖面获取单元,用于在三维数据体G(xi,offj,tl)中,按照垂直旅行时对各个偏移距对应的数据进行叠加,得到二维叠加剖面P(xi,tl);倾角和非断点识别单元,用于基于二维叠加剖面P(xi,tl)得到成像空间内反射面的地层倾角并识别三维数据体G(xi,offj,tl)中的断点和非断点;非断点滤波单元,用于针对每一被识别为非断点的像点、沿所述地层倾角方向设置多项式拟合窗口,进行二维多项式拟合以更新该像点。
本发明的各方面通过对二维(即叠前道集的偏移距和CDP方向两个维度)多项式拟合离散化和具体化,从而有利于保持地震信号振幅的连续性,同时,本发明中沿局部地层倾角方向进行二维多项式拟合来消除随机噪音,能够实现自动处理断层边界,保证断层断点的清晰度。
附图说明
通过结合附图对本发明示例性实施方式进行更详细的描述,本发明的上述以及其它目的、特征和优势将变得更加明显,其中,在本发明示例性实施方式中,相同的参考标号通常代表相同部件。
图1示出了根据本发明的一个实施例的自适应二维多项式拟合的滤波方法的流程图。
图2示出了某原始CDP道集。
图3示出了应用本发明对图2所示原始CDP道集进行滤波后得到的CDP道集。
图4示出了某实际工区的原始二维叠加剖面。
图5示出了应用本发明对图4所示工区的叠前地震资料去噪后的叠加剖面。
具体实施方式
下面将参照附图更详细地描述本发明的优选实施方式。虽然附图中显示了本发明的优选实施方式,然而应该理解,可以以各种形式实现本发明而不应被这里阐述的实施方式所限制。相反,提供这些实施方式是为了使本发明更加透彻和完整,并且能够将本发明的范围完整地传达给本领域的技术人员。
实施例1
图1示出了根据本发明的一个实施例的自适应二维多项式拟合的滤波方法的流程图。在本实施例中,该方法包括:
步骤101,得到叠前偏移距道集;
步骤102,将叠前偏移距道集整理为维度分别为三维数据体G(xi,offj,tl),xi表示像点的水平位置(即在设定的成像坐标下,像点的水平坐标值),offj表示偏移距长度,tl表示像点的垂直旅行时;
步骤103,在三维数据体G(xi,offj,tl)中,按照垂直旅行时对各个偏移距对应的数据进行叠加,得到二维叠加剖面P(xi,tl);
步骤104,基于二维叠加剖面P(xi,tl)得到成像空间内反射面的地层倾角并识别三维数据体G(xi,offj,tl)中的断点和非断点;
步骤105,针对每一被识别为非断点的像点,沿所述地层倾角方向设置多项式拟合窗口,进行二维多项式拟合以更新该像点。
本实施例中,通过对二维多项式拟合离散化和具体化从而有利于保持地震信号振幅的连续性,并沿局部地层倾角方向进行二维多项式拟合来消除随机噪音,实现自动处理断层边界以保证断层断点的清晰度。
在一个示例中,在上述步骤104中,可以采用扫描法得到所述地层倾角,其具体可以包括:
可以在二维叠加剖面P(xi,tl)内设置扫描窗口,所述扫描窗口可以大于所述多项式拟合窗口,这有利于保证数据的稳定性;
可以在所述扫描窗口内的每一像点p(xi,tl)处、在预设角度范围(例如0°~180°的范围)内以预设角度间隔(例如1°)进行扫描,以计算对应于相应角度α的度量值r(α):
其中,可以设扫描窗口在像点水平坐标方向的窗口长度为2M+1,其对应的像点水平坐标可以依次编号为x-M,...,x0,...,xM,可以为在像点p(xi,tl)处的角度为α的扫描线与像点水平坐标为xk的数据线的交点,dx可以为在像点水平坐标方向的采样间隔,V(xi,tl)可以为对应于像点p(xi,tl)的速度参数;
如果在像点p(xi,tl)处扫描得到的度量值r(α)中的最大度量值不小于预设阈值ε,则可以确定该最大度量值对应的角度α为像点p(xi,tl)处的地层倾角。ε可以根据需要设置,例如,可以设置为0~0.5之间的某个值。
需要注意的是,考虑到二维叠加剖面P(xi,tl)在时间方向为离散采样,可能出现落在离散点外的情况,在这种情况下,可以根据周围的像点插值得到例如,可以取沿时间轴方向的上下两个像点的像素值的算术平均。
进一步地,在上述步骤104中,识别断点和非断点可以包括:
如果在像点p(xi,tl)处扫描得到的度量值r(α)中的最大度量值小于预设阈值ε,则可以标定像点p(xi,tl)为参考断点;
针对三维数据体G(xi,offj,tl)中的每个像点g(xi,offj,tl),如果该像点g(xi,offj,tl)与边界的距离和与在沿该像点的倾角方向上标定的所有参考断点的距离中的最小距离值dmin<dthres_hoid,则可以识别该像点g(xi,offj,tl)为断点;否则,可以识别该像点g(xi,offj,tl)为非断点。dmin可以用该段距离上的采样点的数量来衡量,其中,dthres_hoid可以是设定的边缘断点距离阈值。因为在后续处理中,所用因子是B×B二维个点,例如5×5个点,即目标断点范围为前后各2个点,此时,dthres_hoid=2,当某离散点与边界的距离小于2个点或者与某参考断点间的距离小于2个点时,可以将该离散点视为断点。。
进一步地,可以根据非断点的最小距离值dmin来调整多项式拟合窗口的窗口长度。可以设当前多项式拟合窗口的窗口长度为(2N+1)×(2N+1),针对被识别为非断点的像点g(xi,offj,tl),如果其对应的最小距离值dmin有2≤dmin<N,则可以更新N=dmin;如果其对应的dmin≥N,则可以N保持不变。从而可自适应地调整多项式拟合窗口的大小,有利于进一步提高处理性能。
在一个示例中,在上述步骤105中,可以通过下式得到被识别为非断点的当前像点g(xi,offj,tl)的更新后的像点g_new(xi,offj,tl):
其中,可以设所述多项式拟合窗口的窗口长度为(2N+1)×(2N+1), dx可以为在像点水平坐标方向的采样间隔,V(xk1,tl)可以为对应于像点p(xk1,tl)的速度参数,α可以为像点p(xk1,tl)处的地层倾角。
实施例2
根据本实施例,公开了一种自适应二维多项式拟合的滤波装置。该装置包括:
偏移距道集获取单元,用于得到叠前偏移距道集;
三维数据体获取单元,用于将叠前偏移距道集整理为维度分别为三维数据体G(xi,offj,tl),xi表示像点的水平位置,offj表示偏移距长度,tl表示像点的垂直旅行时;
二维剖面获取单元,用于在三维数据体G(xi,offj,tl)中,按照垂直旅行时对各个偏移距对应的数据进行叠加,得到二维叠加剖面P(xi,tl);
倾角和非断点识别单元,用于基于二维叠加剖面P(xi,tl)得到成像空间内反射面的地层倾角并识别三维数据体G(xi,offj,tl)中的断点和非断点;
非断点滤波单元,用于针对每一被识别为非断点的像点、沿所述地层倾角方向设置多项式拟合窗口,进行二维多项式拟合以更新该像点。
在一个示例中,在倾角和非断点识别单元中,可以采用扫描法得到所述地层倾角,具体为:
可以在二维叠加剖面P(xi,tl)内设置扫描窗口,所述扫描窗口可以大于所述多项式拟合窗口;
可以在所述扫描窗口内的每一像点p(xi,tl)处、在预设角度范围内以预设角度间隔进行扫描,以计算对应于相应角度α的度量值r(α):
其中,可以设扫描窗口在像点水平坐标方向的窗口长度为2M+1,其对应的像点水平坐标可以依次编号为x-M,...,x0,...,xM,可以为在像点p(xi,tl)处的角度为α的扫描线与像点水平坐标为xk的数据线的交点,dx可以为在像点水平坐标方向的采样间隔,V(xi,tl)可以为对应于像点p(xi,tl)的速度参数;
如果在像点p(xi,tl)处扫描得到的度量值r(α)中的最大度量值不小于预设阈值ε,则可以确定该最大度量值对应的角度α为像点p(xi,tl)处的地层倾角。
进一步地,在倾角和非断点识别单元中,识别断点和非断点可以包括:
如果在像点p(xi,tl)处扫描得到的度量值r(α)中的最大度量值小于预设阈值ε,则可以标定像点p(xi,tl)为参考断点;
针对三维数据体G(xi,offj,tl)中的每个像点g(xi,offj,tl),如果该像点g(xi,offj,tl)与边界的距离和与在沿该像点的倾角方向上标定的所有参考断点的距离中的最小距离值dmin<dthres_hold,则可以识别该像点g(xi,offj,tl)为断点;否则,可以识别该像点g(xi,offj,tl)为非断点,其中,dthres_hold可以是设定的断点边缘距离阈值。
进一步地,可以设当前多项式拟合窗口的窗口长度为(2N+1)×(2N+1),针对被识别为非断点的像点g(xi,offj,tl),如果其对应的2≤dmin<N,则可以更新N=dmin;如果其对应的dmin≥N,则N可以保持不变。
在一个示例中,在非断点滤波单元中,针对被识别为非断点的像点g(xi,offj,tl),可以更新得到像点g_new(xi,offj,tl):
其中,可以设所述多项式拟合窗口的窗口长度为(2N+1)×(2N+1),
dx可以为在像点水平坐标方向的采样间隔,V(xk1,tl)可以为对应于像点p(xk1,tl)的速度参数,α可以为像点p(xk1,tl)处的地层倾角。
应用示例
为便于理解本发明实施例的方案及其效果,以下给出一个具体应用示例。本领域技术人员应理解,该示例仅为了便于理解本发明,其任何具体细节并非意在以任何方式限制本发明。
图2示出了某原始CDP道集,图3示出了应用本发明对图2所示原始CDP道集进行滤波后得到的CDP道集,其横坐标均表示偏移距(单位:千米),其纵坐标均表示时间(单位:秒)。对比图2和图3,可以明显看出,应用本发明处理后的叠前CDP道集,能够在保持有效信号的同时显著提高信噪比。
图4示出了某实际工区的原始二维叠加剖面。图5示出了应用本发明对图4所示工区的叠前地震资料去噪后的叠加剖面。对比图4和图5,可以明显看出,应用本发明能够很好地消除随机噪声,并且本发明能够很好地适用于复杂构造区。
本发明可以是系统、方法和/或计算机程序产品。计算机程序产品可以包括计算机可读存储介质,其上载有用于使处理器实现本发明的各个方面的计算机可读程序指令。
计算机可读存储介质可以是可以保持和存储由指令执行设备使用的指令的有形设备。计算机可读存储介质例如可以是――但不限于――电存储设备、磁存储设备、光存储设备、电磁存储设备、半导体存储设备或者上述的任意合适的组合。计算机可读存储介质的更具体的例子(非穷举的列表)包括:便携式计算机盘、硬盘、随机存取存储器(RAM)、只读存储器(ROM)、可擦式可编程只读存储器(EPROM或闪存)、静态随机存取存储器(SRAM)、便携式压缩盘只读存储器(CD-ROM)、数字多功能盘(DVD)、记忆棒、软盘、机械编码设备、例如其上存储有指令的打孔卡或凹槽内凸起结构、以及上述的任意合适的组合。这里所使用的计算机可读存储介质不被解释为瞬时信号本身,诸如无线电波或者其他自由传播的电磁波、通过波导或其他传输媒介传播的电磁波(例如,通过光纤电缆的光脉冲)、或者通过电线传输的电信号。
这里所描述的计算机可读程序指令可以从计算机可读存储介质下载到各个计算/处理设备,或者通过网络、例如因特网、局域网、广域网和/或无线网下载到外部计算机或外部存储设备。网络可以包括铜传输电缆、光纤传输、无线传输、路由器、防火墙、交换机、网关计算机和/或边缘服务器。每个计算/处理设备中的网络适配卡或者网络接口从网络接收计算机可读程序指令,并转发该计算机可读程序指令,以供存储在各个计算/处理设备中的计算机可读存储介质中。
用于执行本发明操作的计算机程序指令可以是汇编指令、指令集架构(ISA)指令、机器指令、机器相关指令、微代码、固件指令、状态设置数据、或者以一种或多种编程语言的任意组合编写的源代码或目标代码,所述编程语言包括面向对象的编程语言—诸如Smalltalk、C++等,以及常规的过程式编程语言—诸如“C”语言或类似的编程语言。计算机可读程序指令可以完全地在用户计算机上执行、部分地在用户计算机上执行、作为一个独立的软件包执行、部分在用户计算机上部分在远程计算机上执行、或者完全在远程计算机或服务器上执行。在涉及远程计算机的情形中,远程计算机可以通过任意种类的网络—包括局域网(LAN)或广域网(WAN)—连接到用户计算机,或者,可以连接到外部计算机(例如利用因特网服务提供商来通过因特网连接)。在一些实施例中,通过利用计算机可读程序指令的状态信息来个性化定制电子电路,例如可编程逻辑电路、现场可编程门阵列(FPGA)或可编程逻辑阵列(PLA),该电子电路可以执行计算机可读程序指令,从而实现本发明的各个方面。
这里参照根据本发明实施例的方法、装置(系统)和计算机程序产品的流程图和/或框图描述了本发明的各个方面。应当理解,流程图和/或框图的每个方框以及流程图和/或框图中各方框的组合,都可以由计算机可读程序指令实现。
这些计算机可读程序指令可以提供给通用计算机、专用计算机或其它可编程数据处理装置的处理器,从而生产出一种机器,使得这些指令在通过计算机或其它可编程数据处理装置的处理器执行时,产生了实现流程图和/或框图中的一个或多个方框中规定的功能/动作的装置。也可以把这些计算机可读程序指令存储在计算机可读存储介质中,这些指令使得计算机、可编程数据处理装置和/或其他设备以特定方式工作,从而,存储有指令的计算机可读介质则包括一个制造品,其包括实现流程图和/或框图中的一个或多个方框中规定的功能/动作的各个方面的指令。
也可以把计算机可读程序指令加载到计算机、其它可编程数据处理装置、或其它设备上,使得在计算机、其它可编程数据处理装置或其它设备上执行一系列操作步骤,以产生计算机实现的过程,从而使得在计算机、其它可编程数据处理装置、或其它设备上执行的指令实现流程图和/或框图中的一个或多个方框中规定的功能/动作。
附图中的流程图和框图显示了根据本发明的多个实施例的系统、方法和计算机程序产品的可能实现的体系架构、功能和操作。在这点上,流程图或框图中的每个方框可以代表一个模块、程序段或指令的一部分,所述模块、程序段或指令的一部分包含一个或多个用于实现规定的逻辑功能的可执行指令。在有些作为替换的实现中,方框中所标注的功能也可以以不同于附图中所标注的顺序发生。例如,两个连续的方框实际上可以基本并行地执行,它们有时也可以按相反的顺序执行,这依所涉及的功能而定。也要注意的是,框图和/或流程图中的每个方框、以及框图和/或流程图中的方框的组合,可以用执行规定的功能或动作的专用的基于硬件的系统来实现,或者可以用专用硬件与计算机指令的组合来实现。
以上已经描述了本发明的各实施例,上述说明是示例性的,并非穷尽性的,并且也不限于所披露的各实施例。在不偏离所说明的各实施例的范围和精神的情况下,对于本技术领域的普通技术人员来说许多修改和变更都是显而易见的。本文中所用术语的选择,旨在最好地解释各实施例的原理、实际应用或对市场中的技术的改进,或者使本技术领域的其它普通技术人员能理解本文披露的各实施例。
Claims (10)
1.一种自适应二维多项式拟合的滤波方法,该方法包括:
得到叠前偏移距道集;
将叠前偏移距道集整理为维度分别为三维数据体G(xi,offj,tl),xi表示像点的水平位置,offj表示偏移距长度,tl表示像点的垂直旅行时;
在三维数据体G(xi,offj,tl)中,按照垂直旅行时对各个偏移距对应的数据进行叠加,得到二维叠加剖面P(xi,tl);
基于二维叠加剖面P(xi,tl)得到成像空间内反射面的地层倾角并识别三维数据体G(xi,offj,tl)中的断点和非断点;
针对每一被识别为非断点的像点,沿所述地层倾角方向设置多项式拟合窗口,进行二维多项式拟合以更新该像点。
2.根据权利要求1所述的滤波方法,其中,采用扫描法得到所述地层倾角,具体为:
在二维叠加剖面P(xi,tl)内设置扫描窗口,所述扫描窗口大于所述多项式拟合窗口;
在所述扫描窗口内的每一像点p(xi,tl)处、在预设角度范围内以预设角度间隔进行扫描,以计算对应于相应角度α的度量值r(α):
<mrow>
<mi>r</mi>
<mrow>
<mo>(</mo>
<mi>&alpha;</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mo>|</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>k</mi>
<mo>=</mo>
<mo>-</mo>
<mi>M</mi>
</mrow>
<mi>M</mi>
</munderover>
<mi>p</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mi>k</mi>
</msub>
<mo>,</mo>
<msub>
<mi>t</mi>
<mi>l</mi>
</msub>
<mo>-</mo>
<mfrac>
<mrow>
<mi>d</mi>
<mi>x</mi>
<mo>*</mo>
<mi>k</mi>
<mo>*</mo>
<mi>t</mi>
<mi>a</mi>
<mi>n</mi>
<mi>&alpha;</mi>
</mrow>
<mrow>
<mi>V</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mi>i</mi>
</msub>
<mo>,</mo>
<msub>
<mi>t</mi>
<mi>l</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<mo>|</mo>
<mo>/</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>k</mi>
<mo>=</mo>
<mo>-</mo>
<mi>M</mi>
</mrow>
<mi>M</mi>
</munderover>
<mo>|</mo>
<mi>p</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mi>k</mi>
</msub>
<mo>,</mo>
<msub>
<mi>t</mi>
<mi>l</mi>
</msub>
<mo>-</mo>
<mfrac>
<mrow>
<mi>d</mi>
<mi>x</mi>
<mo>*</mo>
<mi>k</mi>
<mo>*</mo>
<mi>t</mi>
<mi>a</mi>
<mi>n</mi>
<mi>&alpha;</mi>
</mrow>
<mrow>
<mi>V</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mi>i</mi>
</msub>
<mo>,</mo>
<msub>
<mi>t</mi>
<mi>l</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<mo>|</mo>
<mo>,</mo>
</mrow>
其中,设扫描窗口在像点水平坐标方向的窗口长度为2M+1,其对应的像点水平坐标依次编号为x-M,...,x0,...,xM,为在像点p(xi,tl)处的角度为α的扫描线与像点水平坐标为xk的数据线的交点,dx为在像点水平坐标方向的采样间隔,V(xi,tl)为对应于像点p(xi,tl)的速度参数;
如果在像点p(xi,tl)处扫描得到的度量值r(α)中的最大度量值不小于预设阈值ε,则确定该最大度量值对应的角度α为像点p(xi,tl)处的地层倾角。
3.根据权利要求2所述的方法,其中,识别断点和非断点包括:
如果在像点p(xi,tl)处扫描得到的度量值r(α)中的最大度量值小于预设阈值ε,则标定像点p(xi,tl)为参考断点;
针对三维数据体G(xi,offj,tl)中的每个像点g(xi,offj,tl),如果该像点g(xi,offj,tl)与边界的距离和与在沿该像点的倾角方向上标定的所有参考断点的距离中的最小距离值dmin<dthres_hold,则识别该像点g(xi,offj,tl)为断点;否则,识别该像点g(xi,offj,tl)为非断点,其中,dthres_hold是设定的断点边缘距离阈值。
4.根据权利要求3所述的方法,其中,设当前多项式拟合窗口的窗口长度为(2N+1)×(2N+1),针对被识别为非断点的像点g(xi,offj,tl),如果其对应的dthres_hold≤dmin<N,则更新N=dmin;如果其对应的dmin≥N,则N保持不变。
5.根据权利要求1所述的方法,其中,针对被识别为非断点的像点g(xi,offj,tl),更新得到像点g_new(xi,offj,tl):
<mrow>
<msub>
<mi>g</mi>
<mrow>
<mo>_</mo>
<mi>n</mi>
<mi>e</mi>
<mi>w</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mi>i</mi>
</msub>
<mo>,</mo>
<msub>
<mi>off</mi>
<mi>j</mi>
</msub>
<mo>,</mo>
<msub>
<mi>t</mi>
<mi>l</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>m</mi>
<mn>0</mn>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>n</mi>
<mn>2</mn>
</msub>
<mo>+</mo>
<msub>
<mi>n</mi>
<mn>3</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msub>
<mi>n</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>m</mi>
<mrow>
<mn>2</mn>
<mi>i</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>m</mi>
<mrow>
<mn>2</mn>
<mi>j</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>n</mi>
<mn>0</mn>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>n</mi>
<mn>2</mn>
</msub>
<mo>+</mo>
<msub>
<mi>n</mi>
<mn>3</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mn>2</mn>
<msubsup>
<mi>n</mi>
<mn>1</mn>
<mn>2</mn>
</msubsup>
</mrow>
</mfrac>
<mo>,</mo>
</mrow>
1
其中,设所述多项式拟合窗口的窗口长度为(2N+1)×(2N+1), dx为在像点水平坐标方向的采样间隔,V(xk1,tl)为对应于像点p(xk1,tl)的速度参数,α为像点p(xk1,tl)处的地层倾角。
6.一种自适应二维多项式拟合的滤波装置,该装置包括:
偏移距道集获取单元,用于得到叠前偏移距道集;
三维数据体获取单元,用于将叠前偏移距道集整理为维度分别为三维数据体G(xi,offj,tl),xi表示像点的水平位置,offj表示偏移距长度,tl表示像点的垂直旅行时;
二维剖面获取单元,用于在三维数据体G(xi,offj,tl)中,按照垂直旅行时对各个偏移距对应的数据进行叠加,得到二维叠加剖面P(xi,tl);
倾角和非断点识别单元,用于基于二维叠加剖面P(xi,tl)得到成像空间内反射面的地层倾角并识别三维数据体G(xi,offj,tl)中的断点和非断点;
非断点滤波单元,用于针对每一被识别为非断点的像点、沿所述地层倾角方向设置多项式拟合窗口,进行二维多项式拟合以更新该像点。
7.根据权利要求6所述的滤波装置,其中,在倾角和非断点识别单元中,采用扫描法得到所述地层倾角,具体为:
在二维叠加剖面P(xi,tl)内设置扫描窗口,所述扫描窗口大于所述多项式拟合窗口;
在所述扫描窗口内的每一像点p(xi,tl)处、在预设角度范围内以预设角度间隔进行扫描,以计算对应于相应角度α的度量值r(α):
<mrow>
<mi>r</mi>
<mrow>
<mo>(</mo>
<mi>&alpha;</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mo>|</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>k</mi>
<mo>=</mo>
<mo>-</mo>
<mi>M</mi>
</mrow>
<mi>M</mi>
</munderover>
<mi>p</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mi>k</mi>
</msub>
<mo>,</mo>
<msub>
<mi>t</mi>
<mi>l</mi>
</msub>
<mo>-</mo>
<mfrac>
<mrow>
<mi>d</mi>
<mi>x</mi>
<mo>*</mo>
<mi>k</mi>
<mo>*</mo>
<mi>t</mi>
<mi>a</mi>
<mi>n</mi>
<mi>&alpha;</mi>
</mrow>
<mrow>
<mi>V</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mi>i</mi>
</msub>
<mo>,</mo>
<msub>
<mi>t</mi>
<mi>l</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<mo>|</mo>
<mo>/</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>k</mi>
<mo>=</mo>
<mo>-</mo>
<mi>M</mi>
</mrow>
<mi>M</mi>
</munderover>
<mo>|</mo>
<mi>p</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mi>k</mi>
</msub>
<mo>,</mo>
<msub>
<mi>t</mi>
<mi>l</mi>
</msub>
<mo>-</mo>
<mfrac>
<mrow>
<mi>d</mi>
<mi>x</mi>
<mo>*</mo>
<mi>k</mi>
<mo>*</mo>
<mi>t</mi>
<mi>a</mi>
<mi>n</mi>
<mi>&alpha;</mi>
</mrow>
<mrow>
<mi>V</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mi>i</mi>
</msub>
<mo>,</mo>
<msub>
<mi>t</mi>
<mi>l</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<mo>|</mo>
<mo>,</mo>
</mrow>
其中,设扫描窗口在像点水平坐标方向的窗口长度为2M+1,其对应的像点水平坐标依次编号为x-M,...,x0,...,xM,为在像点p(xi,tl)处的角度为α的扫描线与像点水平坐标为xk的数据线的交点,dx为在像点水平坐标方向的采样间隔,V(xi,tl)为对应于像点p(xi,tl)的速度参数;
如果在像点p(xi,tl)处扫描得到的度量值r(α)中的最大度量值不小于预设阈值ε,则确定该最大度量值对应的角度α为像点p(xi,tl)处的地层倾角。
8.根据权利要求7所述的装置,其中,在倾角和非断点识别单元中,识别断点和非断点包括:
如果在像点p(xi,tl)处扫描得到的度量值r(α)中的最大度量值小于预设阈值ε,则标定像点p(xi,tl)为参考断点;
针对三维数据体G(xi,offj,tl)中的每个像点g(xi,offj,tl),如果该像点g(xi,offj,tl)与边界的距离和与在沿该像点的倾角方向上标定的所有参考断点的距离中的最小距离值dmin<dthres_hold,则识别该像点g(xi,offj,tl)为断点;否则,识别该像点g(xi,offj,tl)为非断点,其中,dthres_hold是设定的断点边缘距离阈值。
9.根据权利要求8所述的装置,其中,设当前多项式拟合窗口的窗口长度为(2N+1)×(2N+1),针对被识别为非断点的像点g(xi,offj,tl),如果其对应的dthres_hold≤dmin<N,则更新N=dmin;如果其对应的dmin≥N,则N保持不变。
10.根据权利要求5所述的装置,其中,在非断点滤波单元中,针对被识别为非断点的像点g(xi,offj,tl),更新得到像点g_new(xi,offj,tl):
<mrow>
<msub>
<mi>g</mi>
<mrow>
<mo>_</mo>
<mi>n</mi>
<mi>e</mi>
<mi>w</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mi>i</mi>
</msub>
<mo>,</mo>
<msub>
<mi>off</mi>
<mi>j</mi>
</msub>
<mo>,</mo>
<msub>
<mi>t</mi>
<mi>l</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>m</mi>
<mn>0</mn>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>n</mi>
<mn>2</mn>
</msub>
<mo>+</mo>
<msub>
<mi>n</mi>
<mn>3</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msub>
<mi>n</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>m</mi>
<mrow>
<mn>2</mn>
<mi>i</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>m</mi>
<mrow>
<mn>2</mn>
<mi>j</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>n</mi>
<mn>0</mn>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>n</mi>
<mn>2</mn>
</msub>
<mo>+</mo>
<msub>
<mi>n</mi>
<mn>3</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mn>2</mn>
<msubsup>
<mi>n</mi>
<mn>1</mn>
<mn>2</mn>
</msubsup>
</mrow>
</mfrac>
<mo>,</mo>
</mrow>
其中,设所述多项式拟合窗口的窗口长度为(2N+1)×(2N+1), dx为在像点水平坐标方向的采样间隔,V(xk1,tl)为对应于像点p(xk1,tl)的速度参数,α为像点p(xk1,tl)处的地层倾角。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610350006.0A CN107422382A (zh) | 2016-05-24 | 2016-05-24 | 自适应二维多项式拟合的滤波方法和装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610350006.0A CN107422382A (zh) | 2016-05-24 | 2016-05-24 | 自适应二维多项式拟合的滤波方法和装置 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN107422382A true CN107422382A (zh) | 2017-12-01 |
Family
ID=60422712
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610350006.0A Pending CN107422382A (zh) | 2016-05-24 | 2016-05-24 | 自适应二维多项式拟合的滤波方法和装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107422382A (zh) |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2006111543A1 (en) * | 2005-04-19 | 2006-10-26 | Cgg Services | Method of processing seismic data for avo or avoa characterisation |
CN102681014A (zh) * | 2012-05-23 | 2012-09-19 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 基于多项式拟合的规则线性干扰压制方法 |
CN102866421A (zh) * | 2012-09-04 | 2013-01-09 | 中国科学院地质与地球物理研究所 | 识别小断距断点的散射波叠前成像方法 |
CN105242309A (zh) * | 2015-10-19 | 2016-01-13 | 中国石油天然气集团公司 | 一种地震数据中规则干扰的压制方法及装置 |
-
2016
- 2016-05-24 CN CN201610350006.0A patent/CN107422382A/zh active Pending
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2006111543A1 (en) * | 2005-04-19 | 2006-10-26 | Cgg Services | Method of processing seismic data for avo or avoa characterisation |
CN102681014A (zh) * | 2012-05-23 | 2012-09-19 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 基于多项式拟合的规则线性干扰压制方法 |
CN102866421A (zh) * | 2012-09-04 | 2013-01-09 | 中国科学院地质与地球物理研究所 | 识别小断距断点的散射波叠前成像方法 |
CN105242309A (zh) * | 2015-10-19 | 2016-01-13 | 中国石油天然气集团公司 | 一种地震数据中规则干扰的压制方法及装置 |
Non-Patent Citations (1)
Title |
---|
秦培锐: "自适应二维多项式拟合在叠前道集去噪中的应用", 《长江大学学报(自科版)》 * |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111164462B (zh) | 一种人工源面波勘探方法、面波勘探装置及终端设备 | |
Liu et al. | Seismic data interpolation beyond aliasing using regularized nonstationary autoregression | |
US10705238B2 (en) | Method and apparatus for processing seismic data | |
CN103582825B (zh) | 同步源地震数据的变频滤波 | |
CN102854533B (zh) | 一种基于波场分离原理提高地震资料信噪比的去噪方法 | |
CN107272062B (zh) | 一种数据驱动的地下介质q场估计方法 | |
CN106094024A (zh) | 一种可控震源滑动扫描数据采集处理方法及装置 | |
CN102831588B (zh) | 一种三维地震图像的降噪处理方法 | |
AU2018225300A1 (en) | Generating geophysical images using directional oriented wavefield imaging | |
CN108693555A (zh) | 智能化时变盲反褶积宽频处理方法及装置 | |
WO2018026875A1 (en) | Surface consistent statics solution and amplification correction | |
CN103926623A (zh) | 一种压制逆时偏移低频噪音的方法 | |
CN113156515A (zh) | 一种声波远探测成像降噪处理方法和装置 | |
CN102736108B (zh) | 基于样条拟合的真三维地震数据噪声压制方法 | |
Abedi et al. | Nonhyperbolic stretch-free normal moveout correction | |
CN103901469B (zh) | 地震数据的恢复方法 | |
CN111624658A (zh) | 深度域成像模拟方法及系统 | |
CN105510975A (zh) | 提高地震数据信噪比的方法及装置 | |
CN107255833B (zh) | 地震叠加剖面的确定方法和装置 | |
Shang et al. | Seismic random noise suppression using an adaptive nonlocal means algorithm | |
CN114428343A (zh) | 基于归一化互相关的Marchenko成像方法及系统 | |
CN108037533B (zh) | 一种基于块重组期望对数似然的地震勘探噪声压制方法 | |
CN116068644B (zh) | 一种利用生成对抗网络提升地震数据分辨率和降噪的方法 | |
CN104459792B (zh) | 一种构造约束下的边缘保护滤波方法 | |
CN107422382A (zh) | 自适应二维多项式拟合的滤波方法和装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20171201 |