CN111736221A - 振幅保真度确定方法及系统 - Google Patents
振幅保真度确定方法及系统 Download PDFInfo
- Publication number
- CN111736221A CN111736221A CN202010410797.8A CN202010410797A CN111736221A CN 111736221 A CN111736221 A CN 111736221A CN 202010410797 A CN202010410797 A CN 202010410797A CN 111736221 A CN111736221 A CN 111736221A
- Authority
- CN
- China
- Prior art keywords
- gather
- reflection energy
- seismic data
- average reflection
- determining
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 55
- 238000004590 computer program Methods 0.000 claims description 15
- 230000015572 biosynthetic process Effects 0.000 claims description 7
- 238000003786 synthesis reaction Methods 0.000 claims description 6
- 238000012545 processing Methods 0.000 abstract description 24
- 230000008569 process Effects 0.000 abstract description 17
- 230000008859 change Effects 0.000 abstract description 10
- 238000012544 monitoring process Methods 0.000 abstract description 9
- 239000002131 composite material Substances 0.000 description 6
- 238000010586 diagram Methods 0.000 description 6
- 238000004458 analytical method Methods 0.000 description 5
- 230000005540 biological transmission Effects 0.000 description 4
- 238000011161 development Methods 0.000 description 3
- 230000006870 function Effects 0.000 description 3
- 238000013461 design Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 238000003672 processing method Methods 0.000 description 2
- 238000004891 communication Methods 0.000 description 1
- 238000007405 data analysis Methods 0.000 description 1
- 238000005553 drilling Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 239000000835 fiber Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000004321 preservation Methods 0.000 description 1
- 238000004451 qualitative analysis Methods 0.000 description 1
- 238000003908 quality control method Methods 0.000 description 1
- 238000004445 quantitative analysis Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
Images
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
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/63—Seismic attributes, e.g. amplitude, polarity, instant phase
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)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明提供一种振幅保真度确定方法及系统。该振幅保真度确定方法包括:确定叠前合成记录道集的平均反射能量、原始地震数据道集的平均反射能量和去噪后地震数据道集的平均反射能量;根据叠前合成记录道集的平均反射能量和原始地震数据道集的平均反射能量确定第一相似度;根据叠前合成记录道集的平均反射能量和去噪后地震数据道集的平均反射能量确定第二相似度;根据第一相似度和第二相似度确定该去噪后地震数据道集的振幅保真度。本发明可以确定去噪处理前后数据振幅保真度的变化,进而对地震数据处理过程进行量化的质量监控。
Description
技术领域
本发明涉及地震数据技术领域,具体地,涉及一种振幅保真度确定方法及系统。
背景技术
地震数据数字处理的目的就是改造野外地震数据并从中提取有用的地质信息,为地震勘探的地质解释提供可靠信息。要达到这个目的,地震数据处理就要满足三高的要求,即高信噪比、高分辨率和高保真度。其中,提高信噪比,突出有效信号是地震数据的首要任务,也是最为重要的任务。但为了后续储层分析和预测能正确、顺利地进行,需要在振幅保真的前提下进行目标性处理,并以此作为使用的处理手段是否可行的一个衡量工具。
“保真度”是指输出结果对研究对象的描述所能保持的真实性程度。保持真实性的程度高,即称其为保真度高,否则为保真度低。反射地震信号保真分析就是确定信号处理前后信号强度保持相对不变,即处理后的地震信号对某一反射点的记录在各道都是一致的,没有被“人工修饰”。当然,就目前的技术条件看,在一个相长的时期内达到完全意义上的保真处理是非常困难的,而保真度的提高会使储层预测、油藏描述精度及钻井成功率等有一个较大提升。
近年来,开发地震在老油田的勘探开发中的作用越来越大,与之相应的岩性目标处理技术也得到了很大的发展,人们更加重视地震数据的相对保持振幅、频率、相位和波形的处理,采用叠前数据沿层振幅随炮检距变化情况分析和反演来获得更为丰富、更为精细的地下信息,这对叠前地震数据处理能否保真提出了更高的要求。为此,也出现了各种相对保持振幅的处理方法。
保真度可用来表示地震记录的样点值代表反射系数的程度,只有在相邻反射系数的间距超过分辨率极限的情况下才能讨论。子波旁瓣幅度越大,保真度越差。子波旁瓣数目多、延续长,一个子波的主峰可能与前后多个子波的旁瓣重叠,保真度也可能变差。其中一个保真度确定方法为利用调谐曲线来观察保真度。由于调谐效应,反射系数的大小会与合成记录的反射振幅有差异。所以这种方法一般用于在地层厚度超过调谐厚度的情况下进行叠后数据的分析。另外一种方法主要涉及逐炮的能量对比和时间方向反射能量的变化曲线对比。但这种方案方法主要是对处理后的叠加结果反射能量进行定性的分析,没有考虑叠前数据沿层振幅变化趋势,也没有进行定量的分析对比。
叠前反演或沿层振幅随炮检距变化情况分析的理论基础就是反射波的振幅比接近于相应的反射界面的反射系数比,这就要求前期的地震数据处理能够保持反射波的相对振幅关系,特别是保持角道集数据中不同入射角振幅的变化规律与反射点不同入射角反射系数的变化趋势一致。要满足这个要求,就要对每一个处理步骤都进行相对振幅保持分析,确定所使用的处理方法是否是“保真的”。
发明内容
本发明实施例的主要目的在于提供一种振幅保真度确定方法及系统,以确定去噪处理前后数据振幅保真度的变化,进而对地震数据处理过程进行量化的质量监控。
为了实现上述目的,本发明实施例提供一种振幅保真度确定方法,包括:
确定叠前合成记录道集的平均反射能量、原始地震数据道集的平均反射能量和去噪后地震数据道集的平均反射能量;
根据叠前合成记录道集的平均反射能量和原始地震数据道集的平均反射能量确定第一相似度;
根据叠前合成记录道集的平均反射能量和去噪后地震数据道集的平均反射能量确定第二相似度;
根据第一相似度和第二相似度确定该去噪后地震数据道集的振幅保真度。
本发明实施例还提供一种振幅保真度确定系统,包括:
平均反射能量确定单元,用于确定叠前合成记录道集的平均反射能量、原始地震数据道集的平均反射能量和去噪后地震数据道集的平均反射能量;
第一相似度确定单元,用于根据叠前合成记录道集的平均反射能量和原始地震数据道集的平均反射能量确定第一相似度;
第二相似度确定单元,用于根据叠前合成记录道集的平均反射能量和去噪后地震数据道集的平均反射能量确定第二相似度;
振幅保真度单元,用于根据第一相似度和第二相似度确定该去噪后地震数据道集的振幅保真度。
本发明实施例还提供一种计算机设备,包括存储器、处理器及存储在存储器上并在处理器上运行的计算机程序,处理器执行计算机程序时实现所述的振幅保真度确定方法的步骤。
本发明实施例还提供一种计算机可读存储介质,其上存储有计算机程序,计算机程序被处理器执行时实现所述的振幅保真度确定方法的步骤。
本发明实施例的振幅保真度确定方法及系统先根据叠前合成记录道集和原始地震数据道集的平均反射能量确定第一相似度,根据叠前合成记录道集和去噪后地震数据道集的平均反射能量确定第二相似度,再根据第一相似度和第二相似度确定该去噪后地震数据道集的振幅保真度,可以确定去噪处理前后数据振幅保真度的变化,进而对地震数据处理过程进行量化的质量监控。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1是本发明实施例中振幅保真度确定方法的流程图;
图2是本发明实施例中S102的流程图;
图3是本发明实施例中S103的流程图;
图4是本发明实施例中测井数据的示意图;
图5是本发明实施例中叠前合成记录道集、原始地震数据道集和去噪后地震数据道集的对比示意图;
图6是本发明实施例中平均反射能量随入射角变化的曲线对比示意图;
图7是本发明实施例中振幅保真度确定系统的结构框图;
图8是本发明实施例中计算机设备的结构框图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
本领域技术人员知道,本发明的实施方式可以实现为一种系统、装置、设备、方法或计算机程序产品。因此,本公开可以具体实现为以下形式,即:完全的硬件、完全的软件(包括固件、驻留软件、微代码等),或者硬件和软件结合的形式。
由于噪音污染,叠前地震数据的信噪比通常比较低,尤其是在进行叠前沿层振幅分析用于识别储层含油气性时,相对保持振幅的去噪和信号增强是关键环节。常规的去噪方法和分析信号增强处理后评价信号振幅的保真度是必要的质控手段。
鉴于现有技术无法确定准确的保真度,本发明实施例提供一种振幅保真度确定方法,以确定去噪处理前后数据振幅保真度的变化,进而对地震数据处理过程进行量化的质量监控。以下结合附图对本发明进行详细说明。
图1是本发明实施例中振幅保真度确定方法的流程图。如图1所示,振幅保真度确定方法包括:
S101:确定叠前合成记录道集的平均反射能量、原始地震数据道集的平均反射能量和去噪后地震数据道集的平均反射能量。
一实施例中,确定叠前合成记录道集的平均反射能量包括:
1、根据测井数据确定不同入射角度的地震波的纵波反射系数。
对单一反射界面,假设地震波从上覆介质入射,可以通过如下公式确定地震波的纵波反射系数:
其中,测井数据为已经完成时深标定后的测井数据,包括:入射角θ、横波反射角φ1、上覆介质的纵波速度α1、上覆介质的横波速度β1、上覆介质的密度ρ1、下覆介质的纵波速度α2、下覆介质的横波速度β2、下覆介质的密度ρ2、纵波透射角θ2和横波透射角φ2。RPP(Rpp(t,θ))为地震波的纵波反射系数,RPS为地震波的转换波反射系数,TPP为地震波的纵波反射透射系数,TPS为地震波的转换波透射系数。
2、根据纵波反射系数和预先获得的子波确定叠前合成记录道集。
具体实施时,可以通过如下公式确定叠前合成记录道集:
ssyn(t,θ)=w(t)×Rpp(t,θ);
其中,ssyn(t,θ)为t时刻入射角为θ的叠前合成记录道集,w(t)为t时刻的子波(理论子波或统计方法估算的地震子波),Rpp(t,θ)为t时刻入射角为θ的纵波反射系数。
图4是本发明实施例中测井数据的示意图。如图4所示,图4的横坐标分别为上覆介质的密度,上覆介质的纵波速度、上覆介质的横波速度、地震波的纵波反射系数和叠前合成记录道集,纵坐标为时间。
3、根据叠前合成记录道集确定叠前合成记录道集的平均反射能量。
具体实施时,可以通过如下公式确定叠前合成记录道集在第k个时窗[(k-1)l,kl]的平均反射能量:
其中,Esyn(k,θ)为第k个时窗入射角为θ的叠前合成记录道集的平均反射能量,ssyn(t,θ)为t时刻入射角为θ的叠前合成记录道集(地震道反射振幅值),时窗长度为l。k的取值范围是[1,N],N=L/l,L为地震道长,N为时窗总个数;时窗长度l不应大于子波长度。
S102:根据叠前合成记录道集的平均反射能量和原始地震数据道集的平均反射能量确定第一相似度。
图2是本发明实施例中S102的流程图。如图2所示,S102包括:
S201:根据叠前合成记录道集的平均反射能量确定叠前合成记录道集的平均反射能量的标准方差。
具体实施时,可以通过如下公式确定叠前合成记录道集的平均反射能量的标准方差:
S202:根据原始地震数据道集的平均反射能量确定原始地震数据道集的平均反射能量的标准方差。
在执行S202之前,可以通过如下公式确定原始地震数据道集的平均反射能量:
其中,Eraw(k,θ)为第k个时窗入射角为θ的原始地震数据道集的平均反射能量,sraw(t,θ)为t时刻入射角为θ的原始地震数据道集。
具体实施时,可以通过如下公式确定原始地震数据道集的平均反射能量的标准方差:
S203:根据叠前合成记录道集的平均反射能量和原始地震数据道集的平均反射能量确定叠前合成记录道集的平均反射能量与原始地震数据道集的平均反射能量的协方差。
具体实施时,可以通过如下公式确定叠前合成记录道集的平均反射能量与原始地震数据道集的平均反射能量的协方差:
Cov(Esyn(k,θ),Eraw(k,θ))=(Eraw(k,θ)-σraw(k))(Esyn(k,θ)-σsyn(k));
其中,Cov(Esyn(k,θ),Eraw(k,θ))为叠前合成记录道集的平均反射能量Esyn(k,θ)与原始地震数据道集的平均反射能量Eraw(k,θ)的协方差。
S204:根据叠前合成记录道集的平均反射能量的标准方差、原始地震数据道集的平均反射能量的标准方差和叠前合成记录道集的平均反射能量与原始地震数据道集的平均反射能量的协方差确定第一相似度。
具体实施时,可以通过如下公式确定第一相似度:
其中,rraw(k)为第k个时窗的第一相似度。
S103:根据叠前合成记录道集的平均反射能量和去噪后地震数据道集的平均反射能量确定第二相似度。
图5是本发明实施例中叠前合成记录道集、原始地震数据道集和去噪后地震数据道集的对比示意图。如图5所示,图5的横坐标从左到右依次为叠前合成记录道集、原始地震数据道集和去噪后地震数据道集,纵坐标为时间。图6是本发明实施例中平均反射能量随入射角变化的曲线对比示意图。如图6所示,图6的横坐标为入射角,纵坐标为1100ms处的时窗[-25,25]的平均反射能量。图6中的点划线为叠前合成记录道集的平均反射能量随入射角变化的曲线,实线为原始地震数据道集的平均反射能量随入射角变化的曲线,虚线为去噪后地震数据道集的平均反射能量随入射角变化的曲线。
图3是本发明实施例中S103的流程图。如图3所示,S103包括:
S301:根据叠前合成记录道集的平均反射能量确定叠前合成记录道集的平均反射能量的标准方差。
S302:根据去噪后地震数据道集的平均反射能量确定去噪后地震数据道集的平均反射能量的标准方差。
在执行S302之前,可以通过如下公式确定去噪后地震数据道集的平均反射能量:
其中,Enew(k,θ)为第k个时窗入射角为θ的去噪后地震数据道集的平均反射能量,snew(t,θ)为t时刻入射角为θ的去噪后地震数据道集。
具体实施时,可以通过如下公式确定去噪后地震数据道集的平均反射能量的标准方差:
S303:根据叠前合成记录道集的平均反射能量和去噪后地震数据道集的平均反射能量确定叠前合成记录道集的平均反射能量与去噪后地震数据道集的平均反射能量的协方差。
具体实施时,可以通过如下公式确定叠前合成记录道集的平均反射能量与去噪后地震数据道集的平均反射能量的协方差:
Cov(Esyn(k,θ),Enew(k,θ))=(Enew(k,θ)-σnew(k))(Esyn(k,θ)-σsyn(k));
其中,Cov(Esyn(k,θ),Enew(k,θ))为叠前合成记录道集的平均反射能量Esyn(k,θ)与去噪后地震数据道集的平均反射能量Enew(k,θ)的协方差。
S304:根据叠前合成记录道集的平均反射能量的标准方差、去噪后地震数据道集的平均反射能量的标准方差和叠前合成记录道集的平均反射能量与去噪后地震数据道集的平均反射能量的协方差确定第二相似度。
具体实施时,可以通过如下公式确定第二相似度:
其中,rnew(k)为第k个时窗的第二相似度。
S104:根据第一相似度和第二相似度确定该去噪后地震数据道集的振幅保真度。
具体实施时,可以通过如下公式确定振幅保真度:
其中,p(k)为第k个时窗的振幅保真度。
图1所示的振幅保真度确定方法的执行主体可以为计算机。由图1所示的流程可知,本发明实施例的振幅保真度确定方法先根据叠前合成记录道集和原始地震数据道集的平均反射能量确定第一相似度,根据叠前合成记录道集和去噪后地震数据道集的平均反射能量确定第二相似度,再根据第一相似度述第二相似度确定该去噪后地震数据道集的振幅保真度,可以确定去噪处理前后数据振幅保真度的变化,进而对地震数据处理过程进行量化的质量监控。
本发明实施例的具体流程如下:
1、根据测井数据确定不同入射角度的地震波的纵波反射系数,根据纵波反射系数和预先获得的子波确定叠前合成记录道集。
2、根据叠前合成记录道集确定叠前合成记录道集的平均反射能量。
3、根据叠前合成记录道集的平均反射能量确定叠前合成记录道集的平均反射能量的标准方差,根据叠前合成记录道集的平均反射能量确定叠前合成记录道集的平均反射能量的标准方差。
4、根据叠前合成记录道集的平均反射能量和原始地震数据道集的平均反射能量确定叠前合成记录道集的平均反射能量与原始地震数据道集的平均反射能量的协方差。
5、根据叠前合成记录道集的平均反射能量的标准方差、原始地震数据道集的平均反射能量的标准方差和叠前合成记录道集的平均反射能量与原始地震数据道集的平均反射能量的协方差确定第一相似度。
6、根据叠前合成记录道集的平均反射能量确定叠前合成记录道集的平均反射能量的标准方差,根据去噪后地震数据道集的平均反射能量确定去噪后地震数据道集的平均反射能量的标准方差。
7、根据叠前合成记录道集的平均反射能量和去噪后地震数据道集的平均反射能量确定叠前合成记录道集的平均反射能量与去噪后地震数据道集的平均反射能量的协方差。
8、根据叠前合成记录道集的平均反射能量的标准方差、去噪后地震数据道集的平均反射能量的标准方差和叠前合成记录道集的平均反射能量与去噪后地震数据道集的平均反射能量的协方差确定第二相似度。
9、根据第一相似度和第二相似度确定该去噪后地震数据道集的振幅保真度。
综上,本发明实施例的振幅保真度确定方法先根据叠前合成记录道集和原始地震数据道集的平均反射能量确定第一相似度,根据叠前合成记录道集和去噪后地震数据道集的平均反射能量确定第二相似度,再根据第一相似度和第二相似度确定该去噪后地震数据道集的振幅保真度,可以确定去噪处理前后数据振幅保真度的变化,进而对地震数据处理过程进行量化的质量监控。
基于同一发明构思,本发明实施例还提供了一种振幅保真度确定系统,由于该系统解决问题的原理与振幅保真度确定方法相似,因此该系统的实施可以参见方法的实施,重复之处不再赘述。
图7是本发明实施例中振幅保真度确定系统的结构框图。如图7所示,振幅保真度确定系统包括:
平均反射能量确定单元,用于确定叠前合成记录道集的平均反射能量、原始地震数据道集的平均反射能量和去噪后地震数据道集的平均反射能量;
第一相似度确定单元,用于根据叠前合成记录道集的平均反射能量和原始地震数据道集的平均反射能量确定第一相似度;
第二相似度确定单元,用于根据叠前合成记录道集的平均反射能量和去噪后地震数据道集的平均反射能量确定第二相似度;
振幅保真度单元,用于根据第一相似度和第二相似度确定该去噪后地震数据道集的振幅保真度。
在其中一种实施例中,第一相似度确定单元具体用于:
根据叠前合成记录道集的平均反射能量确定叠前合成记录道集的平均反射能量的标准方差;
根据原始地震数据道集的平均反射能量确定原始地震数据道集的平均反射能量的标准方差;
根据叠前合成记录道集的平均反射能量和原始地震数据道集的平均反射能量确定叠前合成记录道集的平均反射能量与原始地震数据道集的平均反射能量的协方差;
根据叠前合成记录道集的平均反射能量的标准方差、原始地震数据道集的平均反射能量的标准方差和叠前合成记录道集的平均反射能量与原始地震数据道集的平均反射能量的协方差确定第一相似度。
在其中一种实施例中,第二相似度确定单元具体用于:
根据叠前合成记录道集的平均反射能量确定叠前合成记录道集的平均反射能量的标准方差;
根据去噪后地震数据道集的平均反射能量确定去噪后地震数据道集的平均反射能量的标准方差;
根据叠前合成记录道集的平均反射能量和去噪后地震数据道集的平均反射能量确定叠前合成记录道集的平均反射能量与去噪后地震数据道集的平均反射能量的协方差;
根据叠前合成记录道集的平均反射能量的标准方差、去噪后地震数据道集的平均反射能量的标准方差和叠前合成记录道集的平均反射能量与去噪后地震数据道集的平均反射能量的协方差确定第二相似度。
在其中一种实施例中,平均反射能量确定单元具体用于:
根据测井数据确定地震波的纵波反射系数;
根据纵波反射系数和预先获得的子波确定叠前合成记录道集;
根据叠前合成记录道集确定叠前合成记录道集的平均反射能量。
综上,本发明实施例的振幅保真度确定系统先根据叠前合成记录道集和原始地震数据道集的平均反射能量确定第一相似度,根据叠前合成记录道集和去噪后地震数据道集的平均反射能量确定第二相似度,再根据第一相似度和第二相似度确定该去噪后地震数据道集的振幅保真度,可以确定去噪处理前后数据振幅保真度的变化,进而对地震数据处理过程进行量化的质量监控。
本发明实施例还提供能够实现上述实施例中的振幅保真度确定方法中全部步骤的一种计算机设备的具体实施方式。图8是本发明实施例中计算机设备的结构框图,参见图8,所述计算机设备具体包括如下内容:
处理器(processor)801和存储器(memory)802。
所述处理器801用于调用所述存储器802中的计算机程序,所述处理器执行所述计算机程序时实现上述实施例中的振幅保真度确定方法中的全部步骤,例如,所述处理器执行所述计算机程序时实现下述步骤:
确定叠前合成记录道集的平均反射能量、原始地震数据道集的平均反射能量和去噪后地震数据道集的平均反射能量;
根据叠前合成记录道集的平均反射能量和原始地震数据道集的平均反射能量确定第一相似度;
根据叠前合成记录道集的平均反射能量和去噪后地震数据道集的平均反射能量确定第二相似度;
根据第一相似度和第二相似度确定该去噪后地震数据道集的振幅保真度。
综上,本发明实施例的计算机设备先根据叠前合成记录道集和原始地震数据道集的平均反射能量确定第一相似度,根据叠前合成记录道集和去噪后地震数据道集的平均反射能量确定第二相似度,再根据第一相似度和第二相似度确定该去噪后地震数据道集的振幅保真度,可以确定去噪处理前后数据振幅保真度的变化,进而对地震数据处理过程进行量化的质量监控。
本发明实施例还提供能够实现上述实施例中的振幅保真度确定方法中全部步骤的一种计算机可读存储介质,所述计算机可读存储介质上存储有计算机程序,该计算机程序被处理器执行时实现上述实施例中的振幅保真度确定方法的全部步骤,例如,所述处理器执行所述计算机程序时实现下述步骤:
确定叠前合成记录道集的平均反射能量、原始地震数据道集的平均反射能量和去噪后地震数据道集的平均反射能量;
根据叠前合成记录道集的平均反射能量和原始地震数据道集的平均反射能量确定第一相似度;
根据叠前合成记录道集的平均反射能量和去噪后地震数据道集的平均反射能量确定第二相似度;
根据第一相似度和第二相似度确定该去噪后地震数据道集的振幅保真度。
综上,本发明实施例的计算机可读存储介质先根据叠前合成记录道集和原始地震数据道集的平均反射能量确定第一相似度,根据叠前合成记录道集和去噪后地震数据道集的平均反射能量确定第二相似度,再根据第一相似度和第二相似度确定该去噪后地震数据道集的振幅保真度,可以确定去噪处理前后数据振幅保真度的变化,进而对地震数据处理过程进行量化的质量监控。
以上所述的具体实施例,对本发明的目的、技术方案和有益效果进行了进一步详细说明,所应理解的是,以上所述仅为本发明的具体实施例而已,并不用于限定本发明的保护范围,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
本领域技术人员还可以了解到本发明实施例列出的各种说明性逻辑块(illustrative logical block),单元,和步骤可以通过电子硬件、电脑软件,或两者的结合进行实现。为清楚展示硬件和软件的可替换性(interchangeability),上述的各种说明性部件(illustrative components),单元和步骤已经通用地描述了它们的功能。这样的功能是通过硬件还是软件来实现取决于特定的应用和整个系统的设计要求。本领域技术人员可以对于每种特定的应用,可以使用各种方法实现所述的功能,但这种实现不应被理解为超出本发明实施例保护的范围。
本发明实施例中所描述的各种说明性的逻辑块,或单元,或装置都可以通过通用处理器,数字信号处理器,专用集成电路(ASIC),现场可编程门阵列或其它可编程逻辑装置,离散门或晶体管逻辑,离散硬件部件,或上述任何组合的设计来实现或操作所描述的功能。通用处理器可以为微处理器,可选地,该通用处理器也可以为任何传统的处理器、控制器、微控制器或状态机。处理器也可以通过计算装置的组合来实现,例如数字信号处理器和微处理器,多个微处理器,一个或多个微处理器联合一个数字信号处理器核,或任何其它类似的配置来实现。
本发明实施例中所描述的方法或算法的步骤可以直接嵌入硬件、处理器执行的软件模块、或者这两者的结合。软件模块可以存储于RAM存储器、闪存、ROM存储器、EPROM存储器、EEPROM存储器、寄存器、硬盘、可移动磁盘、CD-ROM或本领域中其它任意形式的存储媒介中。示例性地,存储媒介可以与处理器连接,以使得处理器可以从存储媒介中读取信息,并可以向存储媒介存写信息。可选地,存储媒介还可以集成到处理器中。处理器和存储媒介可以设置于ASIC中,ASIC可以设置于用户终端中。可选地,处理器和存储媒介也可以设置于用户终端中的不同的部件中。
在一个或多个示例性的设计中,本发明实施例所描述的上述功能可以在硬件、软件、固件或这三者的任意组合来实现。如果在软件中实现,这些功能可以存储与电脑可读的媒介上,或以一个或多个指令或代码形式传输于电脑可读的媒介上。电脑可读媒介包括电脑存储媒介和便于使得让电脑程序从一个地方转移到其它地方的通信媒介。存储媒介可以是任何通用或特殊电脑可以接入访问的可用媒体。例如,这样的电脑可读媒体可以包括但不限于RAM、ROM、EEPROM、CD-ROM或其它光盘存储、磁盘存储或其它磁性存储装置,或其它任何可以用于承载或存储以指令或数据结构和其它可被通用或特殊电脑、或通用或特殊处理器读取形式的程序代码的媒介。此外,任何连接都可以被适当地定义为电脑可读媒介,例如,如果软件是从一个网站站点、服务器或其它远程资源通过一个同轴电缆、光纤电缆、双绞线、数字用户线(DSL)或以例如红外、无线和微波等无线方式传输的也被包含在所定义的电脑可读媒介中。所述的碟片(disk)和磁盘(disc)包括压缩磁盘、镭射盘、光盘、DVD、软盘和蓝光光盘,磁盘通常以磁性复制数据,而碟片通常以激光进行光学复制数据。上述的组合也可以包含在电脑可读媒介中。
Claims (10)
1.一种振幅保真度确定方法,其特征在于,包括:
确定叠前合成记录道集的平均反射能量、原始地震数据道集的平均反射能量和去噪后地震数据道集的平均反射能量;
根据所述叠前合成记录道集的平均反射能量和所述原始地震数据道集的平均反射能量确定第一相似度;
根据所述叠前合成记录道集的平均反射能量和所述去噪后地震数据道集的平均反射能量确定第二相似度;
根据所述第一相似度和所述第二相似度确定该去噪后地震数据道集的振幅保真度。
2.根据权利要求1所述的振幅保真度确定方法,其特征在于,确定第一相似度包括:
根据所述叠前合成记录道集的平均反射能量确定所述叠前合成记录道集的平均反射能量的标准方差;
根据所述原始地震数据道集的平均反射能量确定所述原始地震数据道集的平均反射能量的标准方差;
根据所述叠前合成记录道集的平均反射能量和所述原始地震数据道集的平均反射能量确定所述叠前合成记录道集的平均反射能量与所述原始地震数据道集的平均反射能量的协方差;
根据所述叠前合成记录道集的平均反射能量的标准方差、所述原始地震数据道集的平均反射能量的标准方差和所述叠前合成记录道集的平均反射能量与所述原始地震数据道集的平均反射能量的协方差确定第一相似度。
3.根据权利要求1所述的振幅保真度确定方法,其特征在于,确定第二相似度包括:
根据所述叠前合成记录道集的平均反射能量确定所述叠前合成记录道集的平均反射能量的标准方差;
根据所述去噪后地震数据道集的平均反射能量确定所述去噪后地震数据道集的平均反射能量的标准方差;
根据所述叠前合成记录道集的平均反射能量和所述去噪后地震数据道集的平均反射能量确定所述叠前合成记录道集的平均反射能量与所述去噪后地震数据道集的平均反射能量的协方差;
根据所述叠前合成记录道集的平均反射能量的标准方差、所述去噪后地震数据道集的平均反射能量的标准方差和所述叠前合成记录道集的平均反射能量与所述去噪后地震数据道集的平均反射能量的协方差确定第二相似度。
4.根据权利要求1所述的振幅保真度确定方法,其特征在于,确定叠前合成记录道集的平均反射能量包括:
根据测井数据确定地震波的纵波反射系数;
根据所述纵波反射系数和预先获得的子波确定叠前合成记录道集;
根据所述叠前合成记录道集确定所述叠前合成记录道集的平均反射能量。
5.一种振幅保真度确定系统,其特征在于,包括:
平均反射能量确定单元,用于确定叠前合成记录道集的平均反射能量、原始地震数据道集的平均反射能量和去噪后地震数据道集的平均反射能量;
第一相似度确定单元,用于根据所述叠前合成记录道集的平均反射能量和所述原始地震数据道集的平均反射能量确定第一相似度;
第二相似度确定单元,用于根据所述叠前合成记录道集的平均反射能量和所述去噪后地震数据道集的平均反射能量确定第二相似度;
振幅保真度单元,用于根据所述第一相似度和所述第二相似度确定该去噪后地震数据道集的振幅保真度。
6.根据权利要求5所述的振幅保真度确定系统,其特征在于,所述第一相似度确定单元具体用于:
根据所述叠前合成记录道集的平均反射能量确定所述叠前合成记录道集的平均反射能量的标准方差;
根据所述原始地震数据道集的平均反射能量确定所述原始地震数据道集的平均反射能量的标准方差;
根据所述叠前合成记录道集的平均反射能量和所述原始地震数据道集的平均反射能量确定所述叠前合成记录道集的平均反射能量与所述原始地震数据道集的平均反射能量的协方差;
根据所述叠前合成记录道集的平均反射能量的标准方差、所述原始地震数据道集的平均反射能量的标准方差和所述叠前合成记录道集的平均反射能量与所述原始地震数据道集的平均反射能量的协方差确定第一相似度。
7.根据权利要求5所述的振幅保真度确定系统,其特征在于,所述第二相似度确定单元具体用于:
根据所述叠前合成记录道集的平均反射能量确定所述叠前合成记录道集的平均反射能量的标准方差;
根据所述去噪后地震数据道集的平均反射能量确定所述去噪后地震数据道集的平均反射能量的标准方差;
根据所述叠前合成记录道集的平均反射能量和所述去噪后地震数据道集的平均反射能量确定所述叠前合成记录道集的平均反射能量与所述去噪后地震数据道集的平均反射能量的协方差;
根据所述叠前合成记录道集的平均反射能量的标准方差、所述去噪后地震数据道集的平均反射能量的标准方差和所述叠前合成记录道集的平均反射能量与所述去噪后地震数据道集的平均反射能量的协方差确定第二相似度。
8.根据权利要求5所述的振幅保真度确定系统,其特征在于,所述平均反射能量确定单元具体用于:
根据测井数据确定地震波的纵波反射系数;
根据所述纵波反射系数和预先获得的子波确定叠前合成记录道集;
根据所述叠前合成记录道集确定所述叠前合成记录道集的平均反射能量。
9.一种计算机设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,其特征在于,所述处理器执行所述计算机程序时实现权利要求1至4任一项所述的振幅保真度确定方法的步骤。
10.一种计算机可读存储介质,其上存储有计算机程序,其特征在于,所述计算机程序被处理器执行时实现权利要求1至4任一项所述的振幅保真度确定方法的步骤。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010410797.8A CN111736221B (zh) | 2020-05-15 | 2020-05-15 | 振幅保真度确定方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010410797.8A CN111736221B (zh) | 2020-05-15 | 2020-05-15 | 振幅保真度确定方法及系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111736221A true CN111736221A (zh) | 2020-10-02 |
CN111736221B CN111736221B (zh) | 2023-08-22 |
Family
ID=72647295
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010410797.8A Active CN111736221B (zh) | 2020-05-15 | 2020-05-15 | 振幅保真度确定方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111736221B (zh) |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2011109839A2 (en) * | 2010-03-05 | 2011-09-09 | Vialogy Llc | Active noise injection computations for improved predictability in oil and gas reservoir discovery and characterization |
US20120253681A1 (en) * | 2011-03-31 | 2012-10-04 | Chevron U.S.A. Inc. | System and method for processing seismic data |
CN102778693A (zh) * | 2011-05-13 | 2012-11-14 | 中国石油化工股份有限公司 | 一种基于反射波层拉平提取并消除的绕射波分离处理方法 |
CN103293552A (zh) * | 2013-05-24 | 2013-09-11 | 中国石油天然气集团公司 | 一种叠前地震资料的反演方法及系统 |
CN104808245A (zh) * | 2015-05-19 | 2015-07-29 | 中国石油化工股份有限公司江汉油田分公司勘探开发研究院 | 道集优化处理方法及其装置 |
CN104820242A (zh) * | 2015-03-27 | 2015-08-05 | 中国石油化工股份有限公司 | 一种面向叠前反演的道集振幅分频补偿方法 |
CN105182419A (zh) * | 2015-09-18 | 2015-12-23 | 电子科技大学 | 一种基于全局优化的叠前地震信号同相轴拉平方法 |
US20150369938A1 (en) * | 2011-03-31 | 2015-12-24 | Chevron U.S.A. Inc. | System and method for processing seismic data |
CN108107470A (zh) * | 2017-12-04 | 2018-06-01 | 中国石油天然气集团公司 | 一种地震数据处理方法及装置 |
-
2020
- 2020-05-15 CN CN202010410797.8A patent/CN111736221B/zh active Active
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2011109839A2 (en) * | 2010-03-05 | 2011-09-09 | Vialogy Llc | Active noise injection computations for improved predictability in oil and gas reservoir discovery and characterization |
US20120253681A1 (en) * | 2011-03-31 | 2012-10-04 | Chevron U.S.A. Inc. | System and method for processing seismic data |
CA2816341A1 (en) * | 2011-03-31 | 2012-10-04 | Chevron U.S.A. Inc. | System and method for processing seismic data |
US20150369938A1 (en) * | 2011-03-31 | 2015-12-24 | Chevron U.S.A. Inc. | System and method for processing seismic data |
CN102778693A (zh) * | 2011-05-13 | 2012-11-14 | 中国石油化工股份有限公司 | 一种基于反射波层拉平提取并消除的绕射波分离处理方法 |
CN103293552A (zh) * | 2013-05-24 | 2013-09-11 | 中国石油天然气集团公司 | 一种叠前地震资料的反演方法及系统 |
CN104820242A (zh) * | 2015-03-27 | 2015-08-05 | 中国石油化工股份有限公司 | 一种面向叠前反演的道集振幅分频补偿方法 |
CN104808245A (zh) * | 2015-05-19 | 2015-07-29 | 中国石油化工股份有限公司江汉油田分公司勘探开发研究院 | 道集优化处理方法及其装置 |
CN105182419A (zh) * | 2015-09-18 | 2015-12-23 | 电子科技大学 | 一种基于全局优化的叠前地震信号同相轴拉平方法 |
CN108107470A (zh) * | 2017-12-04 | 2018-06-01 | 中国石油天然气集团公司 | 一种地震数据处理方法及装置 |
Non-Patent Citations (4)
Title |
---|
刘玉金; 李振春; 周卿: "基于局部相似度的叠前非稳态相位校正方法", 石油地球物理勘探 * |
印兴耀,张繁昌,孙成禹著: "《叠前地震反演》", 31 December 2010, pages: 151 - 154 * |
王正军,苏卫民: "有限覆盖密度下采集面元的选择原则", 《中国石油勘探》 * |
王正军,苏卫民: "有限覆盖密度下采集面元的选择原则", 《中国石油勘探》, 30 November 2015 (2015-11-30), pages 46 - 49 * |
Also Published As
Publication number | Publication date |
---|---|
CN111736221B (zh) | 2023-08-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US10705238B2 (en) | Method and apparatus for processing seismic data | |
CN102109612B (zh) | 一种地震波吸收衰减补偿方法 | |
Wu et al. | Noise attenuation for 2-D seismic data by radial-trace time-frequency peak filtering | |
US9435906B2 (en) | Simultaneous wavelet extraction and deconvolution processing in the time domain | |
US20160091623A1 (en) | Method and Device for Estimating Quality Factor Based on Zero Offset Vertical Seismic Profile Data | |
Xiantai et al. | Adaptive picking of microseismic event arrival using a power spectrum envelope | |
CN110687592B (zh) | 一种微地震震相识别初至拾取方法、设备和存储介质 | |
US11630228B1 (en) | Physical embedded deep learning formation pressure prediction method, device, medium and equipment | |
Ji et al. | Effects of nonlinear soil behavior on Kappa (κ): Observations from the KiK‐Net database | |
CN112183407B (zh) | 一种基于时频域谱减法的隧道地震波数据去噪方法及系统 | |
US10429531B2 (en) | Advanced noise reduction in acoustic well logging | |
Zheng et al. | Microseismic event denoising via adaptive directional vector median filters | |
Rao et al. | Detecting karst voids based on dominant frequencies of seismic profiles | |
CN111736221A (zh) | 振幅保真度确定方法及系统 | |
Hu et al. | A first arrival detection method for low SNR microseismic signal | |
CN112782763B (zh) | 一种地震品质因子估算方法、装置、设备及存储介质 | |
Torii et al. | Shear slowness evaluation by parametric inversion of dipole waveforms | |
Feng et al. | Transform learning in the synchrosqueezing frequency domain—A novel denoising strategy for optical fiber seismic records | |
Dong et al. | A multi-scale dense-connection denoising network for DAS-VSP records | |
Li et al. | Adaptive linear TFPF for seismic random noise attenuation | |
Yao et al. | Microseismic signal denoising using simple bandpass filtering based on normal time–frequency transform | |
Zhang et al. | Q estimation via the discriminant method based on error modeling | |
CN112130209A (zh) | 岩溶储层预测方法及装置 | |
Bakulin et al. | Seismic speckle as multiplicative noise explaining land reflections distorted by near-surface scattering | |
Li et al. | A Novel Polarization Estimation Method for Seismic Recordings |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |