买专利,只认龙图腾
首页 专利交易 科技果 科技人才 科技服务 商标交易 会员权益 IP管家助手 需求市场 关于龙图腾
 /  免费注册
到顶部 到底部
清空 搜索

【发明授权】一种消除并行成像处理痕迹的方法、装置及系统_中国石油天然气股份有限公司_201810418162.5 

申请/专利权人:中国石油天然气股份有限公司

申请日:2018-05-04

公开(公告)日:2019-12-10

公开(公告)号:CN108983283B

主分类号:G01V1/28(20060101)

分类号:G01V1/28(20060101);G01V1/36(20060101)

优先权:

专利状态码:有效-授权

法律状态:2019.12.10#授权;2019.01.04#实质审查的生效;2018.12.11#公开

摘要:本申请实施方式公开了一种消除并行成像处理痕迹的方法、装置及系统,其中,消除并行成像处理痕迹的方法包括:分解地震叠前数据体,获得子数据体;对所述子数据体进行克希霍夫叠前深度偏移并行计算,获得子成像体;对每个子成像体经衰减滤波器处理;对衰减滤波器处理之后的子成像体进行叠加,获得最终成像体。本技术方案对每个子成像体应用一个衰减滤波器,解决了克希霍夫深度偏移并行成像在线照明次数正则化导致的并行痕迹问题,彻底消除了分段画弧和垂向强值等并行痕迹,增强了同相轴的连续性,提高了成像剖面的信噪比,成像的断层更加干脆。并且,本技术方案既发挥了并行处理高效率优势,又保证了并行成像的效果。

主权项:1.一种消除并行成像处理痕迹的方法,其特征在于,包括:分解地震叠前数据体,获得子数据体;对所述子数据体进行克希霍夫叠前深度偏移并行计算,获得子成像体;第一个子成像体经衰减滤波器的内部成像右侧区的权函数wrightin和衰减滤波器的外部成像右侧区的权函数wright处理;最后一个子成像体经衰减滤波器的内部成像左侧区的权函数wleftin和衰减滤波器的外部成像左侧区的权函数wleft处理;除第一个子成像体和最后一个子成像体之外的其他子成像体经衰减滤波器的权函数w处理;所述衰减滤波器的权函数w的表达式为: 其中,lmin表示深度偏移输入测线的最小线号,lmax表示深度偏移输入测线的最大线号,line为子成像体中线号;win为内部成像区的权函数,根据衰减滤波器的内部成像右侧区的权函数wrightin和衰减滤波器的内部成像左侧区的权函数wleftin确定;wout为外部成像区的权函数,根据衰减滤波器的外部成像右侧区的权函数wright和衰减滤波器的外部成像左侧区的权函数wleft确定;对衰减滤波器处理之后的子成像体进行叠加,获得最终成像体。

全文数据:一种消除并行成像处理痕迹的方法、装置及系统技术领域[0001]本申请涉及地震数据处理技术领域,特别涉及一种消除并行成像处理痕迹的方法、装置及系统。背景技术[0002]地震勘探是一种利用地下介质弹性和密度的差异,通过观测和分析大地对人工激发地震波的响应,推断地下构造形态和岩层性质的方法。地震勘探是钻探前勘测石油与天然气资源的重要手段,在煤田和工程地质勘察、区域地质研究和地壳研究等方面,地震勘探也得到广泛应用。地震勘探在分层的详细程度和勘查的精度上,都优于其他地球物理勘探方法。[0003]地震勘探通常分为折射波地震勘探和反射波地震勘探,目前在石油和天然气资源勘查领域,主要采用反射波地震勘探。地震勘探过程由地震数据采集、地震数据处理和地震资料解释三个阶段组成。第一阶段是地震数据采集,其主要任务是在油气勘探有利区带,布置二维或三维测线,使用炸药震源或可控震源激发地震波,在测线上等间距布置多个检波器来接收地震波信号,以等时间间隔离散采样地震数据,并以数字形式记录在磁盘上。第二阶段是地震数据处理,其主要任务是利用计算机和软件,加工处理野外采集的地震数据,将地震数据变成能够反映地下构造的地震剖面图和能够反映地下岩性变化的地震波振幅、频率和传播速度等信息图。地震数据处理主要包括噪声压制、反褶积、水平叠加和偏移成像等关键环节。第三个阶段是地震资料解释,其主要任务是分析地震数据处理所提供的资料,确定地下岩层的产状和构造关系,找出有利的含油气目标,提供钻井位置。地震资料解释主要包括构造解释、储层预测和流体识别等。[0004]地震资料偏移成像是地震资料处理的关键环节之一,也是提交处理成果的最后环节。目前常用的偏移成像技术主要分为克希霍夫叠前深度偏移成像和波动方程叠前深度偏移成像。尽管近年来以逆时偏移和全波形反演为代表的波动方程叠前深度偏移成像方法成为研究的热点,但工业上最常用的方法仍然是克希霍夫叠前深度偏移成像。目前工业上都采用大型并行计算机来进行地震资料叠前深度偏移成像处理,经过并行成像处理后的子成像体需要经过合并来形成最终的成像体。实践表明,这种最终成像体中隐含了或多或少的并行痕迹,也就是说在最终成像体中存在并行切块留下的痕迹,这些痕迹表现为分段画弧或垂向强值。发明内容[0005]本申请实施方式的目的是提供一种消除并行成像处理痕迹的方法、装置及系统,对每个并行处理的子成像体应用一个衰减滤波器,解决了克希霍夫叠前深度偏移并行成像在线照明次数正则化导致的并行痕迹问题。[0006]为实现上述目的,本申请实施方式提供一种消除并行成像处理痕迹的方法,包括:[0007]分解地震叠前数据体,获得子数据体;[0008]对所述子数据体进行克希霍夫叠前深度偏移并行计算,获得子成像体;[0009]对每个子成像体经衰减滤波器处理;[0010]对衰减滤波器处理之后的子成像体进行叠加,获得最终成像体。[0011]优选地,对每个子成像体经衰减滤波器处理的步骤包括:[0012]第一个子成像体经衰减滤波器的内部成像右侧区的权函数Wrlghtln和衰减滤波器的外部成像右侧区的权函数Wrlght处理;[0013]最后一个子成像体经衰减滤波器的内部成像左侧区的权函数Wleftin和衰减滤波器的外部成像右侧区的权函数Wleft处理;[0014]除第一个子成像体和最后一个子成像体之外的其他子成像体经衰减滤波器的权函数w处理。[0015]优选地,所述衰减滤波器的权函数w的表达式为:[0016][0017]其中,Imin表示深度偏移输入测线的最小线号,Imax表示深度偏移输入测线的最大线号,line为子成像体中线号;Win为内部成像区的权函数,根据衰减滤波器的内部成像右侧区的权函数Wrightin和衰减滤波器的内部成像左侧区的权函数Wlrftin确定;Wmjt为外部成像区的权函数,根据衰减滤波器的外部成像右侧区的权函数Wright和衰减滤波器的外部成像右侧区的权函数Wleft确定。[0018]优选地,获得子数据体的步骤包括:[0019]按照线号对所述地震叠前数据体进行分解,获得多个子数据体;其中,所述子数据体的个数由所用并行计算机台数确定。[0020]优选地,获得最终成像体的步骤包括:[0021]将每个子成像体在同一成像空间位置上的成像值对应相加,求和结果除以子成像体个数,结果值为最终成像体的成像值。[0022]为实现上述目的,本申请实施方式提供一种消除并行成像处理痕迹的装置,包括:[0023]子数据体模块,用于分解地震叠前数据体,获得子数据体;[0024]子成像体模块,用于对所述子数据体进行克希霍夫叠前深度偏移并行计算,获得子成像体;[0025]衰减滤波器处理模块,用于对每个子成像体经衰减滤波器处理;[0026]叠加模块,用于对衰减滤波器处理之后的子成像体进行叠加,获得最终成像体。[0027]优选地,所述衰减滤波器处理模块包括:[0028]第一个子成像体处理单元,用于第一个子成像体经衰减滤波器的内部成像右侧区的权函数Wrightin和衰减滤波器的外部成像右侧区的权函数Wright处理;[0029]最后一个子成像体处理单元,用于最后一个子成像体经衰减滤波器的内部成像左侦怄的权函数Wleftin和衰减滤波器的外部成像右侧区的权函数Wleft处理;[0030]中间子成像体处理单元,用于除第一个子成像体和最后一个子成像体之外的其他子成像体经衰减滤波器的权函数w处理。[0031]优选地,还包括:衰减滤波器设计模块;其中,所述衰减滤波器设计模块获得的衰减滤波器的权函数W的表达式为:[0032][0033]其中,Imin表示深度偏移输入测线的最小线号,Imax表示深度偏移输入测线的最大线号,line为子成像体中线号;Win为内部成像区的权函数,根据衰减滤波器的内部成像右侧区的权函数Wrightin和衰减滤波器的内部成像左侧区的权函数Wlrftin确定;Wmjt为外部成像区的权函数,根据衰减滤波器的外部成像右侧区的权函数Wright和衰减滤波器的外部成像右侧区的权函数Wleft确定。[0034]优选地,所述子数据体模块按照线号对所述地震叠前数据体进行分解,获得多个子数据体;其中,所述子数据体的个数由所用并行计算机台数确定。[0035]优选地,所述叠加模块将每个子成像体在同一成像空间位置上的成像值对应相加,求和结果除以子成像体个数,结果值为最终成像体的成像值。[0036]为实现上述目的,本申请实施方式提供一种消除并行成像处理痕迹的系统,所述系统包括:存储器和处理器,所述存储器中存储计算机程序,所述计算机程序被所述处理器执行时,实现以下功能:[0037]分解地震叠前数据体,获得子数据体;[0038]对所述子数据体进行克希霍夫叠前深度偏移并行计算,获得子成像体;[0039]对每个子成像体经衰减滤波器处理;[0040]对衰减滤波器处理之后的子成像体进行叠加,获得最终成像体。[0041]优选地,对每个子成像体经衰减滤波器处理,所述计算机程序被所述处理器执行时,实现以下功能:[0042]第一个子成像体经衰减滤波器的内部成像右侧区的权函数Wrlghtln和衰减滤波器的外部成像右侧区的权函数Wrlght处理;[0043]最后一个子成像体经衰减滤波器的内部成像左侧区的权函数Wleftin和衰减滤波器的外部成像右侧区的权函数Wleft处理;[0044]除第一个子成像体和最后一个子成像体之外的其他子成像体经衰减滤波器的权函数W处理。[0045]优选地,获得子数据体,所述计算机程序被所述处理器执行时,实现以下功能:[0046]按照线号对所述地震叠前数据体进行分解,获得多个子数据体;其中,所述子数据体的个数由所用并行计算机台数确定。[0047]优选地,对衰减滤波器处理之后的子成像体进行叠加,所述计算机程序被所述处理器执行时,实现以下功能:[0048]将每个子成像体在同一成像空间位置上的成像值对应相加,求和结果除以子成像体个数,结果值为最终成像体的成像值。[0049]为实现上述目的,本申请实施方式再提供一种电子设备,包括存储器、处理器以及存储在所述存储器上并可在所述处理器上运行的计算机程序,所述处理器执行所述计算机程序时实现上述所述的消除并行成像处理痕迹的方法。[0050]为实现上述目的,本申请实施方式再提供一种可读存储介质,其上存储有计算机程序,所述计算机程序被执行时实现上述所述的消除并行成像处理痕迹的方法的步骤。[0051]由上可见,与现有技术相比较,本技术方案对每个子成像体应用一个衰减滤波器,解决了克希霍夫深度偏移并行成像在线照明次数正则化导致的并行痕迹问题,彻底消除了分段画弧和垂向强值等并行痕迹,增强了同相轴的连续性,提高了成像剖面的信噪比,成像的断层更加干脆。并且,本技术方案既发挥了并行处理高效率优势,又保证了并行成像的效果。附图说明[0052]为了更清楚地说明本申请实施方式或现有技术中的技术方案,下面将对实施方式或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本申请中记载的一些实施方式,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。[0053]图1为本申请实施例提出一种消除并行成像处理痕迹的方法流程图;[0054]图2为本申请实施例的子成像体衰减滤波器及其作用区域的示意图;[0055]图3为本申请实施例的衰减滤波器的内部成像区权函数参数示意图;[0056]图4为本申请实施例的衰减滤波器的外部成像区权函数参数示意图;[0057]图5为本发明实施例中衰减滤波器应用之前的单个成像体横剖面示意图;[0058]图6为本发明实施例中衰减滤波器应用之后的单个成像体横剖面示意图;[0059]图7为本发明实施例中衰减滤波器应用之前的最终成像体横剖面示意图;[0060]图8为本发明实施例中衰减滤波器应用之后的最终成像体横剖面示意图;[0061]图9为本发明实施例中衰减滤波器应用之前的最终成像体深度1700米切片示意图;[0062]图10为本发明实施例中衰减滤波器应用之后的最终成像体深度1700米切片示意图;[0063]图11为本发明实施例中衰减滤波器应用之前的第三孔径成像体深度1700米切片示意图;[0064]图12为本发明实施例中衰减滤波器应用之后的第三孔径成像体深度1700米切片示意图;[0065]图13为本发明实施例中衰减滤波器应用之前的第四孔径成像体深度1700米切片示意图;[0066]图14为本发明实施例中衰减滤波器应用之后的第四孔径成像体深度1700米切片示意图;[0067]图15为本申请实施例提出的一种消除并行成像处理痕迹的装置功能框图;[0068]图16为本申请实施例提出的一种消除并行成像处理痕迹的系统示意图。具体实施方式[0069]为了使本技术领域的人员更好地理解本申请中的技术方案,下面将结合本申请实施方式中的附图,对本申请实施方式中的技术方案进行清楚、完整地描述,显然,所描述的实施方式仅仅是本申请一部分实施方式,而不是全部的实施方式。基于本申请中的实施方式,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施方式,都应当属于本申请保护的范围。[0070]目前常用的偏移成像技术主要分为克希霍夫叠前深度偏移成像和波动方程叠前深度偏移成像。尽管近年来以逆时偏移和全波形反演为代表的波动方程叠前深度偏移成像方法成为研究的热点,但工业上最常用的方法仍然是克希霍夫叠前深度偏移成像。目前工业上都采用大型并行计算机来进行地震资料叠前深度偏移成像处理,经过并行成像处理后的子成像体需要经过合并来形成最终的成像体。实践表明,这种最终成像体中隐含了或多或少的并行痕迹,也就是说在最终成像体中存在并行切块留下的痕迹,这些痕迹表现为分段画弧或垂向强值。分析其原因,这是因为在偏移成像过程中一般都要进行在线照明次数正则化处理,由此在相邻子成像体的重叠区就会保留并行痕迹。而采用离线照明次数正则化处理,即在并行成像处理过程中不进行照明次数正则化处理,而是在合并之后再进行照明次数正则化处理,这样势必需要多一倍的内存和输入输出,以便保存照明次数,这对于海量地震数据而言,现有计算机难以满足此要求。[0071]为了满足海量地震数据克希霍夫叠前深度偏移并行成像处理的计算效率和成像效果,提出一种消除并行成像处理痕迹的技术方案。该技术方案的特点是对每个并行处理的子成像体应用一个衰减滤波器,解决了克希霍夫叠前深度偏移并行成像在线照明次数正则化导致的并行痕迹问题,这样既发挥了并行处理高效率优势,又保证了并行成像的效果。[0072]基于上述描述,本申请实施例提出一种消除并行成像处理痕迹的方法,如图1所示。包括:[0073]步骤101:分解地震叠前数据体,获得子数据体。[0074]在本实施例中,按线号分解地震叠前数据体,形成多个子数据体。本实施例中,整体成像区线号范围为7069-7369,CDP范围为3100-3600,成像区深度为0-3000米,深度样点数nz=601。成像点间距ΔX=25米、Δy=25米、Δz=5米,分别为成像点在x、y和z方向间距。而输入的地震叠前数据体是经过面元规则化的共中心点CDP地震数据,共中心点道集的范围是:线号7069-7369,⑶P3100-3600。将此地震叠前数据体分解为12个子数据体,其线号范围分别是(7069,7098、(7099,7122、(7123,7146、(7147,7170、(7171,7194、7195,7218、(7219,7242、(7243,7266、(7267,7290、(7291,7314、(7315,7338和7339,7369。这样分解地震叠前数据体,确保最终成像体的纵测线内不会出现并行成像处理痕迹,而并行成像处理痕迹只会出现在最终成像体的横测线内。对输入的地震记录有如下要求:经过了常规地震资料预处理,包括振幅均衡处理、频带均衡处理和静校正等。输入的地震道要包含线号道头字。[0075]步骤102:对所述子数据体进行克希霍夫叠前深度偏移并行计算,获得子成像体。[0076]在本实施例中,要求克希霍夫深度偏移输出的子成像体道头字tracr中设定子孔径编号,如果偏移中只用一个孔径时,就设置为1。[0077]步骤103:对每个子成像体经衰减滤波器处理。[0078]在本实施例中,图2显示了子成像体衰减滤波器及其作用区域的示意图。横坐标为线号,纵坐标为地下成像深度。这里显示了3个子成像体。其中,第一个实框表示首个成像体分布区域,第一个虚框表不相应的输入数据线号分布区;第二个实框表不中间成像体分布区域,第二个虚框表不相应的输入数据线号分布区;第三个实框表不最后成像体分布区域,第三个虚框表示相应的输入数据线号分布区。在子成像体重叠区(如图中阴影区所示需要应用衰减滤波器。图中,曲线1表示首个成像体所用的衰减滤波器权函数,从左到右权函数从1降为0;对于中间成像体,左边的权函数从O升为1、右边的权函数从1降为0,如图中曲线2所示;曲线3表示最后成像体所用的衰减滤波器权函数,从左到右权函数从0升为1。图中所示的子成像体重叠区是由三角阴影区和矩形阴影区组成,三角阴影区决定了分段画弧痕迹分布范围,而矩形阴影区决定了垂向强值痕迹分布范围。[0079]由图2可知,首个成像体的衰减滤波器作用区域在右侧、最后成像体的衰减滤波器作用区域在左侧、中间成像体的衰减滤波器作用区域分布在左右两侧。下面我们以中间成像体为例,计算子成像体衰减滤波器的权函数,具体步骤如下:[0080]㈧确定三角阴影区和矩形阴影区的分界点深度Z0。[0081][0082]其中,Zbeg为成像区起始深度米),zref为参考深度米),radius为参考深度上最大孔径米),tabry为走时表y方向范围的一半米),ncycle为孔径总个数,k为子孔径编号k=1,2,…,ncycle。当只用一个孔径偏移时,ncycle=l,k=l。[0083]为了满足不对称情况,提供一个参数ml作为左侧成像区相对减少的道数,以便区分左右两侦扮界点深度。右侦扮界点深度为ZOr=ZO,左侧分界点深度为ZOl:[0084][0085]其中,dyo为y方向成像点网格间距米)。[0086]本实施例中,成像区起始深度Zbeg=0米,参考深度Zref=3000米,参考深度上最大孔径:radius=3500米,走时表y方向范围一半tabry=2200米,孔径总个数ncycle=10,左侧成像区减少的道数ml=7,y方向成像点网格间距dyo=25米,子孔径编号k可以从子成像体的道头tracr中提取。如图5所示,为本发明实施例中衰减滤波器应用之前的单个成像体横剖面示意图。如图6所示,为本发明实施例中衰减滤波器应用之后的单个成像体横剖面示意图。图5和图6显示的是前五个子孔径k=l,2,3,4,5的成像体之和,横测线号为3350。[0087]⑻衰减滤波器设计[0088]将成像区分为内部成像区和外部成像区,其中内部成像区是指输入测线范围内的成像区,外部成像区为输入测线范围之外的成像区域。设Imin和Imax分别为深度偏移输入测线的最小线号和最大线号,line为子成像体中线号。则衰减滤波器的权函数w为:[0089][0090]其中Win为内部成像区的权函数,Wciut为外部成像区的权函数。[0091]本实施例中图5和图6所用的lmin=7195、lmax=7218,这是由步骤101中分解的子数据体线号范围确定的。[0092]a内部成像区的权函数Win[0093]图3显示了内部成像区权函数参数示意图。由于内部成像区中左侧三角区和矩形区是由前一个子成像体的右侧三角区和矩形区决定的,所以图3中左侧三角区和矩形区的分界点深度为Zor。同理,图3中右侧三角区和矩形区的分界点深度为Zc1U图3中参数计算如下:[0094][0095][0096][0097][0098]其中,[0099][0100]为满足互补性权函数要求(即两个重叠窗口中权值之和为1,并且相互对称),我们采用球函数作为权函数。于是内部成像区的权函数Win可表示为:[0101]Win—P〇WX0·5+0·5XWleftinX0·5+0·5XWrightin,1ίηΘImax[0102]其中,pow是在0.3到I之间的一个常数用户参数)。且[0105]其中,z为成像点深度米)。且[0103][0104][0106][0107]这里d=ni_line,a=0.5Xwidth,width为衰减区道数(用户参数)。[0108][0109][0110][0111]这里d=line_n2,a=0.5Xwidth。[0112][0113]本实施例中,常数pow=0.6,衰减区道数width=88道。[0114]⑹外部成像区的权函数Wout[0115]图4显示了外部成像区权函数参数示意图,图中Slmin和Slmax分别为总成像区的最小线号和最大线号,其他参数计算如下:[0116]Ii=Imin-ty+ml[0117]l2=lmax+ty[0118]I3=li+width[0119]U=b-width[0120]外部成像区的权函数win可表示为:[0121][0122][0123][0124][0125]其中,[0126][0127][0128][0129][0130][0131][0132][0133][0134][0135][0136][0137][0138]这里d=line_m4,a=a2rz。[0139]本实施例中,总成像区的最小线号slmin=7069、最大线号slmax=7369。[0140]⑹衰减滤波器应用[0141]以子成像体中线号line作为变量,将每个子成像体乘以这个衰减滤波器w,就可获得衰减滤波之后的子成像体。由于内部成像区左右两侧衰减区可能会重叠,为此可进行如·cl~3Wleftin—Wrightin—I、Wleft—I、£1¾Wright—1日寸,—I;b^Wleftin—I、或者Wrightin=I时,取W=0.5+0.5XwieftinX0.5+0.5XWrightin。对于首个成像体,应用右侧衰减滤波器Wrightin和Wright;对于最后成像体,应用左侧衰减滤波器Wliftin和Wlift;对于中间成像体,应用左右两侧衰减滤波器。[0142]步骤104:对衰减滤波器处理之后的子成像体进行叠加,获得最终成像体。[0143]在本实施例中,将所有衰减滤波之后的子成像体进行叠加,就可以输出无并行痕迹的最终成像体。该过程表示为[0144][0145]其中,Iix,y,z为第i个子成像体该成像空间位置x,y,z上的成像值,N为子成像体个数,IX,y,z为输出的最终成像值。[0M6]图5是本发明实施例中衰减滤波器应用之前的单个成像体横剖面示意图。图6是本发明实施例中衰减滤波器应用之后的单个成像体横剖面示意图。对比图5和图6,可以明显看出衰减滤波器压制的效果。图7是本发明实施例中衰减滤波器应用之前的最终成像体横剖面示意图,图8是本发明实施例中衰减滤波器应用之后的最终成像体横剖面示意图。本实施例中,N=12。对比图7和图8,可以看出,图8显示了横测线号为3350,明显消除了图7中存在的分段画弧和垂向强值等并行成像处理痕迹,增强了同相轴的连续性,提高了成像剖面的信噪比,成像的断层更加干脆。图9是本发明实施例中衰减滤波器应用之前的最终成像体深度1700米切片示意图,图10是本发明实施例中衰减滤波器应用之后的最终成像体深度1700米切片示意图。对比图9和图10,可以看出图10中成像更加清晰,明显减少了横向痕迹。图11是本发明实施例中衰减滤波器应用之前的第三孔径成像体深度1700米切片示意图,图12是本发明实施例中衰减滤波器应用之后的第三孔径成像体深度1700米切片示意图。对比图11和图12,可以看出图12中明显消除了并行成像处理痕迹,成像分辨率和清晰度得到显著提高。图13是本发明实施例中衰减滤波器应用之前的第四孔径成像体深度1700米切片示意图,图14是本发明实施例中衰减滤波器应用之后的第四孔径成像体深度1700米切片示意图。对比图13和图14,可以看出图14中成像非常清晰,并行成像处理痕迹完全消失了。由此实际资料应用效果,可以更加清楚本发明的效果。所以,图10、图12和图14分别显示了衰减滤波器应用之后的深度1700米切片,它们分别对应于最终成像体、第三孔径成像体和第四孔径成像体,明显消除了图9、图11和图13中并行成像处理的痕迹,提高了成像的分辨率和清晰度。对比可以看出,消除并行成像处理痕迹的方法既发挥了并行处理高效率优势,又保证了并行成像的效果。[0M7]如图15所示,为本申请实施例提出的一种消除并行成像处理痕迹的装置功能框图。包括:[0148]子数据体模块201,用于分解地震叠前数据体,获得子数据体;[0149]子成像体模块202,用于对所述子数据体进行克希霍夫叠前深度偏移并行计算,获得子成像体;[0150]衰减滤波器处理模块203,用于对每个子成像体经衰减滤波器处理;[0151]叠加模块204,用于对衰减滤波器处理之后的子成像体进行叠加,获得最终成像体。[0152]在本实施例中,所述衰减滤波器处理模块203包括:[0153]第一个子成像体处理单元,用于第一个子成像体经衰减滤波器的内部成像右侧区的权函数Wrightin和衰减滤波器的外部成像右侧区的权函数Wright处理;[0154]最后一个子成像体处理单元,用于最后一个子成像体经衰减滤波器的内部成像左侦怄的权函数Wleftin和衰减滤波器的外部成像右侧区的权函数Wleft处理;[0155]中间子成像体处理单元,用于除第一个子成像体和最后一个子成像体之外的其他子成像体经衰减滤波器的权函数w处理。[0156]在本实施例中,消除并行成像处理痕迹的装置还包括:衰减滤波器设计模块;其中,所述衰减滤波器设计模块获得的衰减滤波器的权函数w的表达式为:[0157][0158]其中,Imin表示深度偏移输入测线的最小线号,1^^表示深度偏移输入测线的最大线号,line为子成像体中线号;Win为内部成像区的权函数,Wciut为外部成像区的权函数。[0159]在本实施例中,所述子数据体模块按照线号对所述地震叠前数据体进行分解,获得多个子数据体;其中,所述子数据体的个数由所用并行计算机台数确定。[0160]在本实施例中,所述叠加模块将每个子成像体在同一成像空间位置上的成像值对应相加,求和结果除以子成像体个数,结果值为最终成像体的成像值。[0161]如图16所示,为本申请实施例提出的一种消除并行成像处理痕迹的系统示意图。所述系统包括:存储器和处理器,所述存储器中存储计算机程序,所述计算机程序被所述处理器执行时,实现以下功能:[0162]分解地震叠前数据体,获得子数据体;[0163]对所述子数据体进行克希霍夫叠前深度偏移并行计算,获得子成像体;[0164]对每个子成像体经衰减滤波器处理;[0165]对衰减滤波器处理之后的子成像体进行叠加,获得最终成像体。[0166]在本实施例中,对每个子成像体经衰减滤波器处理,所述计算机程序被所述处理器执行时,实现以下功能:[0167]第一个子成像体经衰减滤波器的内部成像右侧区的权函数Wrightin和衰减滤波器的外部成像右侧区的权函数Wrlght处理;[0168]最后一个子成像体经衰减滤波器的内部成像左侧区的权函数Wleftin和衰减滤波器的外部成像右侧区的权函数Wleft处理;[0169]除第一个子成像体和最后一个子成像体之外的其他子成像体经衰减滤波器的权函数w处理。[0170]在本实施例中,获得子数据体,所述计算机程序被所述处理器执行时,实现以下功能:[0171]按照线号对所述地震叠前数据体进行分解,获得多个子数据体;其中,所述子数据体的个数由所用并行计算机台数确定。[0172]在本实施例中,对衰减滤波器处理之后的子成像体进行叠加,所述计算机程序被所述处理器执行时,实现以下功能:[0173]将每个子成像体在同一成像空间位置上的成像值对应相加,求和结果除以子成像体个数,结果值为最终成像体的成像值。[0174]在本实施方式中,所述存储器包括但不限于随机存取存储器(RandomAccessMemory,RAM、只读存储器(Read-OnlyMemory,ROM、缓存(Cache、硬盘(HardDiskDrive,HDD或者存储卡MemoryCard。[0175]在本实施方式中,所述处理器可以按任何适当的方式实现。例如,所述处理器可以采取例如微处理器或处理器以及存储可由该微处理器执行的计算机可读程序代码例如软件或固件)的计算机可读介质、逻辑门、开关、专用集成电路ApplicationSpecificIntegratedCircuit,ASIC、可编程逻辑控制器和嵌入微控制器的形式等等。[0176]本说明书实施方式提供的消除并行成像处理痕迹的系统,其存储器和处理器实现的具体功能,可以与本说明书中的前述实施方式相对照解释,并能够达到前述实施方式的技术效果,这里便不再赘述。[0177]综上,由上述描述可知,本技术方案可以快速消除并行成像处理产生的分段画弧和垂向强值等并行痕迹。本技术方案可用于经过预处理的叠前地震数据,包括共炮集数据、共炮检距道集数据、共中心点道集资料或规则化之后的道集数据。本发明的效果是彻底消除了分段画弧和垂向强值等并行痕迹,增强了同相轴的连续性,提高了成像剖面的信噪比,成像的断层更加干脆。[0178]在20世纪90年代,对于一个技术的改进可以很明显地区分是硬件上的改进例如,对二极管、晶体管、开关等电路结构的改进还是软件上的改进对于方法流程的改进)。然而,随着技术的发展,当今的很多方法流程的改进已经可以视为硬件电路结构的直接改进。设计人员几乎都通过将改进的方法流程编程到硬件电路中来得到相应的硬件电路结构。因此,不能说一个方法流程的改进就不能用硬件实体模块来实现。例如,可编程逻辑器件ProgrammableLogicDevice,PLD例如现场可编程门阵列(FieldProgrammableGateArray,FPGA就是这样一种集成电路,其逻辑功能由用户对器件编程来确定。由设计人员自行编程来把一个数字系统“集成”在一片PLD上,而不需要请芯片制造厂商来设计和制作专用的集成电路芯片。而且,如今,取代手工地制作集成电路芯片,这种编程也多半改用“逻辑编译器logiccompiler”软件来实现,它与程序开发撰写时所用的软件编译器相类似,而要编译之前的原始代码也得用特定的编程语言来撰写,此称之为硬件描述语言HardwareDescriptionLanguage,HDL,而HDL也并非仅有一种,而是有许多种,如ABELAdvancedBooleanExpressionLanguage、AHDLAlteraHardwareDescriptionLanguage、Confluence、CUPLCornellUniversityProgrammingLanguage、HDCaI、JHDLJavaHardwareDescriptionLanguage、Lava、Lola、MyHDL、PALASM、RHDLRubyHardwareDescriptionLanguage等,目前最普遍使用的是VHDLVery-High-SpeedIntegratedCircuitHardwareDescriptionLanguage与Verilog2〇本令页域技术人员也应该清楚,只需要将方法流程用上述几种硬件描述语言稍作逻辑编程并编程到集成电路中,就可以很容易得到实现该逻辑方法流程的硬件电路。[0179]本领域技术人员也知道,除了以纯计算机可读程序代码方式实现客户端、服务器以外,完全可以通过将方法步骤进行逻辑编程来使得客户端、服务器以逻辑门、开关、专用集成电路、可编程逻辑控制器和嵌入微控制器等的形式来实现相同功能。因此这种客户端、服务器可以被认为是一种硬件部件,而对其内部包括的用于实现各种功能的装置也可以视为硬件部件内的结构。或者甚至,可以将用于实现各种功能的装置视为既可以是实现方法的软件模块又可以是硬件部件内的结构。[0180]通过以上的实施方式的描述可知,本领域的技术人员可以清楚地了解到本申请可借助软件加必需的通用硬件平台的方式来实现。基于这样的理解,本申请的技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来,该计算机软件产品可以存储在存储介质中,如R0MRAM、磁碟、光盘等,包括若干指令用以使得一台计算机设备可以是个人计算机,服务器,或者网络设备等执行本申请各个实施方式或者实施方式的某些部分所述的方法。[0181]本说明书中的各个实施方式均采用递进的方式描述,各个实施方式之间相同相似的部分互相参见即可,每个实施方式重点说明的都是与其他实施方式的不同之处。尤其,针对客户端的实施方式来说,均可以参照前述方法的实施方式的介绍对照解释。[0182]本申请可以在由计算机执行的计算机可执行指令的一般上下文中描述,例如程序模块。一般地,程序模块包括执行特定任务或实现特定抽象数据类型的例程、程序、对象、组件、数据结构等等。也可以在分布式计算环境中实践本申请,在这些分布式计算环境中,由通过通信网络而被连接的远程处理设备来执行任务。在分布式计算环境中,程序模块可以位于包括存储设备在内的本地和远程计算机存储介质中。[0183]虽然通过实施方式描绘了本申请,本领域普通技术人员知道,本申请有许多变形和变化而不脱离本申请的精神,希望所附的权利要求包括这些变形和变化而不脱离本申请的精神。

权利要求:1.一种消除并行成像处理痕迹的方法,其特征在于,包括:分解地震叠前数据体,获得子数据体;对所述子数据体进行克希霍夫叠前深度偏移并行计算,获得子成像体;对每个子成像体经衰减滤波器处理;对衰减滤波器处理之后的子成像体进行叠加,获得最终成像体。2.如权利要求1所述的方法,其特征在于,对每个子成像体经衰减滤波器处理的步骤包括:第一个子成像体经衰减滤波器的内部成像右侧区的权函数Wrlghtln和衰减滤波器的外部成像右侧区的权函数Wright处理;最后一个子成像体经衰减滤波器的内部成像左侧区的权函数Wleftin和衰减滤波器的外部成像右侧区的权函数们也处理;除第一个子成像体和最后一个子成像体之外的其他子成像体经衰减滤波器的权函数W处理。3.如权利要求2所述的方法,其特征在于,所述衰减滤波器的权函数w的表达式为:其中,1_表示深度偏移输入测线的最小线号,Imax表示深度偏移输入测线的最大线号,line为子成像体中线号;Win为内部成像区的权函数,根据衰减滤波器的内部成像右侧区的权函数Wrightin和衰减滤波器的内部成像左侧区的权函数Wleftin确定;Wciut为外部成像区的权函数,根据衰减滤波器的外部成像右侧区的权函数Wright和衰减滤波器的外部成像右侧区的权函数Wleft确定。4.如权利要求1或2所述的方法,其特征在于,获得子数据体的步骤包括:按照线号对所述地震叠前数据体进行分解,获得多个子数据体;其中,所述子数据体的个数由所用并行计算机台数确定。5.如权利要求1或2所述的方法,其特征在于,获得最终成像体的步骤包括:将每个子成像体在同一成像空间位置上的成像值对应相加,求和结果除以子成像体个数,结果值为最终成像体的成像值。6.—种消除并行成像处理痕迹的装置,其特征在于,包括:子数据体模块,用于分解地震叠前数据体,获得子数据体;子成像体模块,用于对所述子数据体进行克希霍夫叠前深度偏移并行计算,获得子成像体;衰减滤波器处理模块,用于对每个子成像体经衰减滤波器处理;叠加模块,用于对衰减滤波器处理之后的子成像体进行叠加,获得最终成像体。7.如权利要求6所述的装置,其特征在于,所述衰减滤波器处理模块包括:第一个子成像体处理单元,用于第一个子成像体经衰减滤波器的内部成像右侧区的权函数Wrightin和衰减滤波器的外部成像右侧区的权函数Wright处理;最后一个子成像体处理单元,用于最后一个子成像体经衰减滤波器的内部成像左侧区的权函数们也^和衰减滤波器的外部成像右侧区的权函数Wleft处理;中间子成像体处理单元,用于除第一个子成像体和最后一个子成像体之外的其他子成像体经衰减滤波器的权函数W处理。8.如权利要求7所述的装置,其特征在于,还包括:衰减滤波器设计模块;其中,所述衰减滤波器设计模块获得的衰减滤波器的权函数w的表达式为:其中,1_表示深度偏移输入测线的最小线号,Imax表示深度偏移输入测线的最大线号,line为子成像体中线号;Win为内部成像区的权函数,根据衰减滤波器的内部成像右侧区的权函数Wrightin和衰减滤波器的内部成像左侧区的权函数Wleftin确定;Wciut为外部成像区的权函数,根据衰减滤波器的外部成像右侧区的权函数Wright和衰减滤波器的外部成像右侧区的权函数Wleft确定。9.如权利要求6或7所述的装置,其特征在于,所述子数据体模块按照线号对所述地震叠前数据体进行分解,获得多个子数据体;其中,所述子数据体的个数由所用并行计算机台数确定。10.如权利要求6或7所述的装置,其特征在于,所述叠加模块将每个子成像体在同一成像空间位置上的成像值对应相加,求和结果除以子成像体个数,结果值为最终成像体的成像值。11.一种消除并行成像处理痕迹的系统,其特征在于,所述系统包括:存储器和处理器,所述存储器中存储计算机程序,所述计算机程序被所述处理器执行时,实现以下功能:分解地震叠前数据体,获得子数据体;对所述子数据体进行克希霍夫叠前深度偏移并行计算,获得子成像体;对每个子成像体经衰减滤波器处理;对衰减滤波器处理之后的子成像体进行叠加,获得最终成像体。12.如权利要求11所述的系统,其特征在于,对每个子成像体经衰减滤波器处理,所述计算机程序被所述处理器执行时,实现以下功能:第一个子成像体经衰减滤波器的内部成像右侧区的权函数Wrlghtln和衰减滤波器的外部成像右侧区的权函数Wright处理;最后一个子成像体经衰减滤波器的内部成像左侧区的权函数Wleftin和衰减滤波器的外部成像右侧区的权函数们也处理;除第一个子成像体和最后一个子成像体之外的其他子成像体经衰减滤波器的权函数W处理。13.如权利要求11或12所述的系统,其特征在于,获得子数据体,所述计算机程序被所述处理器执行时,实现以下功能:按照线号对所述地震叠前数据体进行分解,获得多个子数据体;其中,所述子数据体的个数由所用并行计算机台数确定。14.如权利要求11或12所述的系统,其特征在于,对衰减滤波器处理之后的子成像体进行叠加,所述计算机程序被所述处理器执行时,实现以下功能:将每个子成像体在同一成像空间位置上的成像值对应相加,求和结果除以子成像体个数,结果值为最终成像体的成像值。15.—种电子设备,包括存储器、处理器以及存储在所述存储器上并可在所述处理器上运行的计算机程序,其特征在于,所述处理器执行所述计算机程序时实现权利要求1〜5任意一项权利要求所述的消除并行成像处理痕迹的方法。16.—种可读存储介质,其上存储有计算机程序,其特征在于,所述计算机程序被执行时实现权利要求1〜5任意一项权利要求所述的消除并行成像处理痕迹的方法的步骤。

百度查询: 中国石油天然气股份有限公司 一种消除并行成像处理痕迹的方法、装置及系统

免责声明
1、本报告根据公开、合法渠道获得相关数据和信息,力求客观、公正,但并不保证数据的最终完整性和准确性。
2、报告中的分析和结论仅反映本公司于发布本报告当日的职业理解,仅供参考使用,不能作为本公司承担任何法律责任的依据或者凭证。

相关技术
相关技术
相关技术
相关技术