CN107807390B - 地震数据的处理方法及系统 - Google Patents
地震数据的处理方法及系统 Download PDFInfo
- Publication number
- CN107807390B CN107807390B CN201610812740.4A CN201610812740A CN107807390B CN 107807390 B CN107807390 B CN 107807390B CN 201610812740 A CN201610812740 A CN 201610812740A CN 107807390 B CN107807390 B CN 107807390B
- Authority
- CN
- China
- Prior art keywords
- result
- seismic data
- frequency
- function
- analytical function
- 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.)
- Active
Links
- 238000003672 processing method Methods 0.000 title claims abstract description 21
- 238000012545 processing Methods 0.000 claims abstract description 112
- 230000009466 transformation Effects 0.000 claims abstract description 50
- 238000000034 method Methods 0.000 claims abstract description 24
- 238000010276 construction Methods 0.000 claims description 15
- 235000013399 edible fruits Nutrition 0.000 claims description 6
- 238000004364 calculation method Methods 0.000 claims description 4
- 238000010586 diagram Methods 0.000 description 38
- 238000001228 spectrum Methods 0.000 description 18
- 230000008859 change Effects 0.000 description 7
- 230000008901 benefit Effects 0.000 description 5
- 230000000694 effects Effects 0.000 description 5
- 238000005516 engineering process Methods 0.000 description 5
- 238000009795 derivation Methods 0.000 description 3
- 238000003384 imaging method Methods 0.000 description 3
- 238000005070 sampling Methods 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 206010034719 Personality change Diseases 0.000 description 1
- 238000010521 absorption reaction Methods 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000006854 communication Effects 0.000 description 1
- 230000000052 comparative effect Effects 0.000 description 1
- 238000000354 decomposition reaction Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 239000011435 rock Substances 0.000 description 1
- 238000010183 spectrum analysis Methods 0.000 description 1
- 230000002123 temporal effect Effects 0.000 description 1
- 230000002087 whitening effect 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/32—Transforming one recording into another or one representation into another
-
- 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
-
- 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
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/282—Application of seismic models, synthetic seismograms
-
- 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/30—Analysis
- G01V1/307—Analysis for determining seismic attributes, e.g. amplitude, instantaneous phase or frequency, reflection strength or polarity
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/14—Fourier, Walsh or analogous domain transformations, e.g. Laplace, Hilbert, Karhunen-Loeve, transforms
- G06F17/141—Discrete Fourier transforms
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/40—Transforming data representation
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/40—Transforming data representation
- G01V2210/43—Spectral
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geophysics (AREA)
- Mathematical Physics (AREA)
- Computational Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Algebra (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Discrete Mathematics (AREA)
- Geophysics And Detection Of Objects (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
- Radar Systems Or Details Thereof (AREA)
- Complex Calculations (AREA)
Abstract
本发明公开了一种地震数据的处理方法及系统。其中方法包括:获取原始单道地震数据;对原始单道地震数据分别进行傅里叶变换和希尔伯特变换,得到傅里叶变换结果和希尔伯特变换结果;依次对预设的分频范围内的每个待处理频率,根据原始单道地震数据、傅里叶变换结果和希尔伯特变换结果,得到该待处理频率对应的处理结果;根据分频范围内所有待处理频率对应的处理结果,得到输出结果。本发明可有效补偿因常规采集仪器的带宽原因造成的低频端缺失,使得地震资料的频带得到明显的均匀展宽,从而大幅提高了地震资料分辨率,并且避免了由于频带的拓展造成的波形畸变并且保持子波运动学特征不变。
Description
技术领域
本发明涉及数字信号处理技术领域,尤其涉及一种地震数据的处理方法,还涉及一种地震数据的处理系统。
背景技术
油气勘探的目标已经从简单型转向复杂型,从构造油气藏转向隐蔽岩性油气藏。可以看出,地震数据处理技术越来越重要,其中,如何提高地震资料的分辨率是一项至关重要的环节。开展高分辨率地震资料采集的成本较高,而常规采集仪器的频带较窄,特别是5Hz以下的低频端严重缺失。现有的处理成果很难满足岩性油气藏勘探的需要,高分辨率处理技术迫在眉睫。
地震信号分辨率分为纵向分辨率和横向分辨率。一般我们通常所说的分辨率都是针对纵向分辨率而言的。影响分辨率的因素很多,主要有以下几个方面。(1)岩石吸收衰减作用对分辨率的影响。地震波在地下介质传播过程中,其振幅会出现一定的衰减,其衰减与传播距离、频率、Q值的倒数都是呈指数关系。(2)采样率对分辨率的影响。在资料处理中,地震记录是一系列的离散数据,时间采样率的大小直接决定了数据的最高频率。比如1ms采样时,最高频率可达500Hz,就目前所采集的资料来说,即使我们考虑仪器上的去假频滤波器,我们所使用的采样也刚能基本满足目前的分辨率要求。(3)子波频带宽度的影响。地震勘探的分辨率由地震子波的持续时间或脉冲宽度确定。给定一个脉冲的频带宽度,就有一个确定的最小脉冲宽度,即确定了一个潜在的最大分辨率。因此分辨率取决于子波的频带宽度,要提高分辨率,必须有效地展宽子波的有效频带,压缩子波,这也是高分辨率处理要解决的主要问题。(4)子波相位的影响。在振幅谱相同的情况下,零相位子波具有最好的分辨率,因为零相位子波小于其它子波的长度,边沿振幅值小,且反射时间出现在子波的峰值处。因此,理想的情况下最终的剖面上子波应为零相位,目前的技术还不能准确地知道子波的相位,无法提取准确的相位信息,而只能对子波进行统计估计,尽量向零相位靠近,以提高分辨率,目前绝大多数反褶积都是基于这一点。
公开发表的技术文档和参考文献表明,现有采用的主要高分辨率处理技术是各种改进的反褶积算法,如确定性子波反褶积方法、时变谱白化方法、特征值提高分辨率法、独立分量法(ICA)和盲反褶积等。比起常规反褶积,这些方法技术实现的分辨率会得到明显提高,但大多真假难辩。如何在不产生假同相轴的前提下,把分辨率提高到高于常规反褶积的程度,对地震资料处理和解释非常必要。此外,现有技术难以实现地震数据频带的高、低频端的均匀拓展,一味追求高频端的拓展,可能导致地震子波的波形畸变或旅行时即运动学特征发生改变,由此影响速度建模和成像精度。
发明内容
本发明的目的在于:在地震子波不发生波形畸变且不改变其运动学特征的前提下,实现地震数据频带的高频端和低频端均匀展宽,从而有效补偿由常规采集仪器带宽原因造成的数据低频端缺失,使得地震数据的分辨率得到明显提高,为后续地震数据反演及成像的处理和解释提供技术支撑。
为了实现上述目的,本发明提供提供了一种基于分频迭代的地震数据的处理方法及系统。
根据本发明的一个方面,提供了一种地震数据的处理方法,其包括:
获取原始单道地震数据;
对所述原始单道地震数据分别进行傅里叶变换和希尔伯特变换,得到傅里叶变换结果和希尔伯特变换结果;
依次对预设的分频范围内的每个待处理频率,根据所述原始单道地震数据、傅里叶变换结果和希尔伯特变换结果,得到所述待处理频率对应的处理结果;
根据所述分频范围内所有待处理频率对应的处理结果,得到输出结果。
优选的是,根据所述原始单道地震数据、傅里叶变换结果和希尔伯特变换结果,得到所述待处理频率对应的处理结果,包括:
根据得到所述待处理频率f对应的处理结果Kr(t);
其中,处理结果Kr(t)为构建的输出函数K(t)的实部,x(t)为所述原始单道地震数据,Xr(t)为所述傅里叶变换结果的实部,h(t)为所述希尔伯特变换结果。
优选的是,上述处理方法还包括:构建所述输出函数K(t)。
优选的是,构建所述输出函数K(t),包括:
构建第一解析函数E(t),并使所述第一解析函数E(t)满足:E(t)=x(t)+jh(t);
构建第二解析函数Y(t),并使所述第二解析函数Y(t)满足:Y(t)=Xr(t)+j2πft;
根据所述第二解析函数Y(t)构建第三解析函数Z(t),并使所述第三解析函数Z(t)满足:
对所述第一解析函数E(t)和所述第三解析函数Z(t)进行乘积运算,得到所述输出函数K(t),所述输出函数K(t)满足:
优选的是,根据所述分频范围内所有待处理频率对应的处理结果,得到输出结果,包括:
对所述分频范围内所有待处理频率对应的处理结果进行求和运算,得到所述输出结果。
根据本发明的另一个方面,提供了一种地震数据的处理系统,其包括:
数据获取模块,设置为获取原始单道地震数据;
傅里叶变换模块,设置为对所述原始单道地震数据进行傅里叶变换,得到傅里叶变换结果;
希尔伯特变换模块,设置为对所述原始单道地震数据进行希尔伯特变换,得到希尔伯特变换结果;
处理结果确定模块,设置为依次对预设的分频范围内的每个待处理频率,根据所述原始单道地震数据、傅里叶变换结果和希尔伯特变换结果,得到所述待处理频率对应的处理结果;
输出结果确定模块,设置为根据所述分频范围内所有待处理频率对应的处理结果,得到输出结果。
优选的是,所述处理结果确定模块具体设置为:
根据得到所述待处理频率f对应的处理结果Kr(t);
其中,处理结果Kr(t)为构建的输出函数K(t)的实部,x(t)为所述原始单道地震数据,Xr(t)为所述傅里叶变换结果的实部,h(t)为所述希尔伯特变换结果。
优选的是,上述处理系统还包括构建模块,所述构建模块设置为构建所述输出函数K(t)。
优选的是,所述构建模块包括:
第一构建单元,设置为构建第一解析函数E(t),并使所述第一解析函数E(t)满足:E(t)=x(t)+jh(t);
第二构建单元,设置为构建第二解析函数Y(t),并使所述第二解析函数Y(t)满足:Y(t)=Xr(t)+j2πft;
第三构建单元,设置为根据所述第二解析函数Y(t)构建第三解析函数Z(t),并使所述第三解析函数Z(t)满足:
输出函数确定单元,对所述第一解析函数E(t)和所述第三解析函数Z(t)进行乘积运算,得到所述输出函数K(t),所述输出函数K(t)满足:
优选的是,所述输出结果确定模块具体设置为:
对所述分频范围内所有待处理频率对应的处理结果进行求和运算,得到所述输出结果。
与现有技术相比,上述方案中的一个或多个实施例可以具有如下优点或有益效果:
本发明可有效补偿因常规采集仪器的带宽原因造成的低频端缺失,使得地震资料的频带得到明显的均匀展宽,从而大幅提高了地震资料分辨率,并且避免了由于频带的拓展造成的波形畸变并且保持子波运动学特征不变。
本发明的其它特征和优点将在随后的说明书中阐述,并且部分地从说明书中变得显而易见,或者通过实施本发明而了解。本发明的目的和其他优点可通过在说明书、权利要求书以及附图中所特别指出的结构来实现和获得。
附图说明
附图用来提供对本发明的进一步理解,并且构成说明书的一部分,与本发明的实施例共同用于解释本发明,并不构成对本发明的限制。在附图中:
图1示出了本发明实施例地震数据的处理方法的一种流程示意图;
图2示出了本发明实施例地震数据的处理方法的另一种流程示意图;
图3示出了本发明实施例中构建输出函数的方法的流程示意图;
图4示出了本发明实施例地震数据的处理系统的一种结构示意图;
图5示出了本发明实施例地震数据的处理系统的另一种结构示意图;
图6示出了本发明实施例中构建模块的结构示意图;
图7a示出了理论子波示意图;
图7b示出了经本发明实施例(分频范围为(0,10))处理后的子波示意图;
图7c示出了原始子波对应的频谱示意图;
图7d示出了经处理的原始子波的频谱示意图;
图8a示出了应用本发明实施例之前靶区CMP道集记录及其速度谱示意图;
图8b示出了应用本发明实施例之后靶区CMP道集记录及其速度谱示意图;
图9a示出了应用本发明实施例之前靶区CMP道集记录的原频谱示意图;
图9b示出了应用本发明实施例(分频范围为(0,10)Hz)之后的频谱示意图;
图10a示出了靶区的第444线原叠加剖面示意图;
图10b示出了应用本发明实施例对图10a所示的第444线进行叠前处理后的叠加剖面示意图;
图11a示出了靶区的第452线原叠加剖面示意图;
图11b示出了应用本发明实施例对图11a所示的第452线进行叠前处理后的叠加剖面示意图;
图12a示出了靶区的第460线原叠加剖面示意图;
图12b示出了应用本发明实施例对图12a所示的第460线进行叠前处理后的叠加剖面示意图。
具体实施方式
以下将结合附图及实施例来详细说明本发明的实施方式,借此对本发明如何应用技术手段来解决技术问题,并达成技术效果的实现过程能充分理解并据以实施。需要说明的是,只要不构成冲突,本发明中的各个实施例以及各实施例中的各个特征可以相互结合,所形成的技术方案均在本发明的保护范围之内。
现有采用的主要高分辨率处理技术是各种改进的反褶积算法。比起常规反褶积,这些方法技术实现的分辨率会得到明显提高,但大多真假难辩。如何在不产生假同相轴的前提下,把分辨率提高到高于常规反褶积的程度,对地震资料处理和解释非常必要。此外,现有技术难以实现地震数据频带的高、低频端的均匀拓展,一味追求高频端的拓展,可能导致地震子波的波形畸变或旅行时即运动学特征发生改变,由此影响速度建模和成像精度。
为了解决上述技术问题,本发明实施例提供了一种地震数据的处理方法。
实施例一
图1示出了本发明实施例地震数据的处理方法的一种流程示意图。如图1所示,本发明实施例主要包括步骤101至步骤106。
在步骤101中,获取原始单道地震数据。这里,原始单道地震数据用x(t)表示。
在步骤102中,对原始单道地震数据进行傅里叶变换(Fourier Transform),得到傅里叶变换结果。
具体地,傅里叶正变换和傅里叶反变换可以实现信号在时间域和频率域的相互转变。通常,利用式(1)对原始单道地震数据x(t)进行傅里叶正变换,得到傅里叶正变换结果X(w)。利用式(2)对傅里叶正变换结果X(w)进行傅里叶反变换,得到原始单道地震数据x(t)。
在式(1)中,Xr(w)为傅里叶变换后的实部,对Xr(w)进行傅里叶反变换,得到傅里叶反变换结果的实部Xr(t)。Xi(w)为傅里叶变换后的虚部,其振幅函数满足其相位函数满足
在频率域,傅里叶变换的典型用途是将信号分解成幅值谱显示,以便对信号进行频谱分析。当原始地震信号低频化时,在频率域表现为子波主频小、频带较窄,因而具有较低的分辨率,这将不利于信号的分析和进一步的地震解释工作。分辨率取决于子波的频带宽度,要提高分辨率,必须有效地展宽子波的有效频带,压缩子波,这也是本发明要解决的主要问题。
在步骤103中,对原始单道地震数据进行希尔伯特变换,得到希尔伯特变换结果。
具体地,希尔伯特变换(HT,Hilbert Transform)是信号分析中的重要工具。利用式(3)对原始单道地震数据x(t)进行希尔伯特变换,得到希尔伯特变换结果h(t)。
在步骤104中,依次对预设的分频范围内的每个待处理频率,根据原始单道地震数据、傅里叶变换结果和希尔伯特变换结果,得到待处理频率对应的处理结果。
具体地,根据式(4)得到待处理频率f对应的处理结果Kr(t)。
在式(4)中,处理结果Kr(t)为构建的输出函数K(t)的实部。输出函数K(t)可以在线构建也可以离线构建,具体的构建方法将在下文中结合图2进行详细地阐述。x(t)为原始单道地震数据,Xr(t)为傅里叶变换结果的实部,h(t)为希尔伯特变换结果。
在步骤105中,判断是否得到分频范围内的所有待处理频率对应的处理结果。
在步骤106中,在判断出得到分频范围内的所有待处理频率对应的处理结果时,根据分频范围内所有待处理频率对应的处理结果,得到输出结果。在判断出未得到分频范围内的所有待处理频率对应的处理结果时,返回步骤104。
具体地,分频范围为离线预设的。这里,利用(fmin,fmax)表示分频范围。其中,fmin表示分频范围的下限,fmax表示分频范围的上限。可以看出,本文所述的分频范围可以看成是待处理频率f的集合。
当判断出未全部得到分频范围内的所有待处理频率对应的处理结果时,返回步骤104继续进行处理。
当判断出全部得到分频范围内的所有待处理频率对应的处理结果时,根据得到的所有处理结果得到输出结果。在本发明一优选的实施例中,对分频范围内所有待处理频率对应的处理结果进行求和运算,得到输出结果。即,根据式(5)得到处理后的实际输出y(t)。
应用本实施例所述的地震数据的处理方法,依次得到分频范围内的每个待处理频率对应的处理结果,然后根据所有的处理结果得到最终的输出结果。此外,在计算每个待处理频率对应处理结果时,引入了原始单道地震数据、傅里叶变换结果和希尔伯特变换结果,并且希尔伯特变换结果用于对三瞬属性(瞬时振幅、瞬时频率和瞬时相位)进行约束。可见,本实施例加入了对原始单道地震数据的约束以及对三瞬属性的约束,从而能够避免信号在频率域变换时可能发生的畸变现象。
总体来讲,本实施例将地震数据进行基于单道的分频迭代高分辨率处理,可以有效提高地震数据的分辨率并且保持子波运动学特征不变,在拓宽地震数据的高频端的同时也完成了低频端的拓宽处理,使得有效信号的频带得到明显的展宽。具体地,单道的分频迭代高分辨率处理是基于傅里叶变换和希尔伯特变换进行解析函数的构建,可以在不同尺度对数据进行单道单频率的高分辨率处理,避免了由于频谱的拓展造成的波形畸变的产生。因此,本实施例大幅提高了地震资料分辨率,为后续的处理和解释提供技术支撑。
实施例二
图2示出了本发明实施例地震数据的处理方法的另一种流程示意图。如图2所示,本实施例在实施例一的基础上增加了步骤201。
在步骤201中,构建输出函数K(t)。这里,可以在线构建输出函数K(t),也可以离线构建输出函数K(t)。
图3示出了本发明实施例中构建输出函数K(t)的方法的流程示意图。如图3所示,本实施例构建输出函数K(t)的方法主要包括步骤301至步骤304。
在步骤301中,构建第一解析函数E(t),并使第一解析函数E(t)满足式(6)。
E(t)=x(t)+jh(t) (6)
具体地,第一解析函数E(t)由原始单道地震数据x(t)和希尔伯特变换结果h(t)构建,从而加入了对输入信号(即原始单道地震数据)和信号三瞬属性的约束,这里,信号三瞬属性指的是瞬时振幅属性、瞬时频率属性和瞬时相位属性。
在步骤302中,构建第二解析函数Y(t),并使第二解析函数Y(t)满足式(7)。
Y(t)=Xr(t)+j2πft (7)
在步骤303中,根据第二解析函数Y(t)构建第三解析函数Z(t),并使第三解析函数Z(t)满足式(8)。
具体地,第二解析函数Y(t)的构建主要是为了构建第三解析函数Z(t),这样逐层的构建方式便于公式的推导。
第三解析函数Z(t)由原始单道地震数据x(t)的傅里叶变换结果的实部Xr(t)和待处理频率f的三角函数(sin(2πft)、cos(2πft))构建。在第三解析函数Z(t)的表达式(8)中,指数函数用于振幅约束,三角函数用于相位约束。
在步骤304中,对第一解析函数E(t)和第三解析函数Z(t)进行乘积运算,得到输出函数K(t),输出函数K(t)满足式(9)。
具体地,将第一解析函数E(t)和第三解析函数Z(t)的乘积作为输出函数K(t),可见输出函数K(t)是在上述几种约束下(具体包括:对输入信号的约束,对信号三瞬属性的约束、对振幅的约束和对相位的约束)构建的关于单道单频率的函数,约束的加入可以避免信号在频率域变换时可能发生的畸变现象。
在本实施例中,第一、第二和第三解析函数的构建主要是为了逐层推导出输出函数而构建的,这样可以在加入约束的同时使得输出函数的推导思路清晰。此外,单道单频率的输出函数K(t)是本发明实施例的核心函数,该输出函数通过对输入信号、三瞬属性、振幅和相位进行约束,来避免信号在频率域变换时可能发生的畸变现象。
综上所述,本实施例所述的地震数据的处理方法,将地震数据进行基于单道的分频迭代高分辨率处理,可以有效提高地震数据的分辨率并且保持子波运动学特征不变,在拓宽地震数据的高频端的同时也完成了低频端的拓宽处理,使得有效信号的频带得到明显的展宽。此外,单道的分频迭代高分辨率处理是基于傅里叶变换和希尔伯特变换进行解析函数的构建,可以在不同尺度对数据进行单道单频率的高分辨率处理,从而避免了由于频谱的拓展造成的波形畸变的产生。由此,本发明实施例大幅提高了地震资料分辨率,为后续的处理和解释提供有力的技术支撑。
为了更好地验证由本发明实施例带来的上述有益效果,下面结合图7a至图12b进行详细地阐述。具体地,通过理论数据和实际数据的处理验证,来验证本发明实施例方法的正确性和有效性。
具体地,图7a示出了一个雷克子波,图7b示出了经本发明实施例所述的处理方法(分频范围为(0,10)Hz)处理后的高分辨率处理结果。由图7a与7b的对比可知,经本发明实施例的方法处理后的子波分辨率得到明显提高,并且子波主瓣对应的时间位置没有发生位移。图7c与图7d是处理前后对应的频谱示意图。对比图7c与图7d可知,处理后的子波主频得到提高,并且频带也得到明显的拓宽。具体地,频带不仅在高频端得到拓展,而且在低频端也得到一定的拓展。这使得波形的频率成分更加丰富。
以下为实际资料应用本发明实施例的处理方法处理前后的结果对比。靶区为中国西部某工区三维地震资料。虽然本发明实施例是单道处理方法,主要针对叠前处理,但是这并不意味着不能针对叠后处理。
图8a示出了应用本发明实施例之前靶区CMP道集记录及其速度谱示意图。图8b示出了应用本发明实施例之后该靶区CMP道集记录及其速度谱示意图。对比图8a和图8b可知,应用本实施例方法处理前后的速度谱能量团的物理位置是吻合的,这表明本发明实施例的处理方法没有改变子波的运动学特征。
图9a示出了应用本发明实施例之前靶区CMP道集记录的原频谱示意图。参照图9a,处理之前频带较窄,分辨率较低,特别是5Hz以下的低频端因地震数据采集仪器的带宽原因而严重缺失。图9b示出了应用本发明实施例(分频范围为(0,10)Hz)之后的频谱示意图。对比图9a和图9b,频带得到有效拓展,尤其是低频端得到了有效补偿。
图10a示出了靶区的第444线(CMP1380-1520,时间2-3秒)原叠加剖面示意图。图10b示出了应用本发明实施例对图10a所示的第444线进行叠前处理后的叠加剖面示意图。图中的白色粗线表示对比提示线。这里,对比提示线并不限于此。
图11a示出了靶区的第452线(CMP760-900,时间2-3秒)原叠加剖面示意图。图11b示出了应用本发明实施例对图11a所示的第452线进行叠前处理后的叠加剖面示意图。图中的白色粗线表示对比提示线。这里,对比提示线并不限于此。
图12a示出了靶区的第460线(CMP1560-1700,时间2-3秒)原叠加剖面示意图。图12b示出了应用本发明实施例对图12a所示的第460线进行叠前处理后的叠加剖面示意图。图中的白色粗线表示对比提示线。这里,对比提示线并不限于此。
对图10a和图10b进行对比,对图11a和图11b进行对比,并且对图12a和图12b进行对比。对比结果表明,应用本发明实施例处理方法进行叠前处理后的分辨率有明显的提高,且未产生假同相轴,不存在真假难辩的现象。可见,应用本实施例处理方法处理的数据,为后续的处理和解释提供了技术支撑。
综上所述,本发明实施例可有效补偿因常规采集仪器的带宽原因造成的低频端缺失,使得地震资料的频带得到明显的均匀展宽,从而大幅提高了地震资料分辨率,并且避免了由于频带的拓展造成的波形畸变并且保持子波运动学特征不变。
采用宽频采集仪器无疑会大幅增加地震数据的采集成本,在相当长的时期内低成本的常规采集仪器仍将发挥重要作用。而本发明实施例的方法优势非常适用于由常规采集仪器采集到的地震资料的高分辨率处理。
实施例三
对应于上述实施例一和实施例二,本发明实施例还提供了一种地震数据的处理系统。
图4示出了本发明实施例地震数据的处理系统的一种结构示意图。如图4所示,本实施例地震数据的处理系统主要包括数据获取模块401、傅里叶变换模块402、希尔伯特变换模块403、处理结果确定模块404和输出结果确定模块405。其中,数据获取模块401分别与傅里叶变换模块402和希尔伯特变换模块403连接。傅里叶变换模块402和希尔伯特变换模块403均与处理结果确定模块404连接。处理结果确定模块404与输出结果确定模块405连接。
具体地,数据获取模块401,设置为获取原始单道地震数据。
傅里叶变换模块402,设置为对原始单道地震数据进行傅里叶变换,得到傅里叶变换结果。
希尔伯特变换模块403,设置为对原始单道地震数据进行希尔伯特变换,得到希尔伯特变换结果。
处理结果确定模块404,设置为依次对预设的分频范围内的每个待处理频率,根据原始单道地震数据、傅里叶变换结果和希尔伯特变换结果,得到待处理频率对应的处理结果。
特别地,处理结果确定模块404具体设置为:根据得到待处理频率f对应的处理结果Kr(t)。其中,处理结果Kr(t)为构建的输出函数K(t)的实部,x(t)为原始单道地震数据,Xr(t)为傅里叶变换结果的实部,h(t)为希尔伯特变换结果。
输出结果确定模块405,设置为根据分频范围内所有待处理频率对应的处理结果,得到输出结果。特别地,输出结果确定模块405具体设置为:对分频范围内所有待处理频率对应的处理结果进行求和运算,得到输出结果。
应用本实施例所述的地震数据的处理系统,依次得到分频范围内的每个待处理频率对应的处理结果,然后根据所有的处理结果得到最终的输出结果。此外,在计算每个待处理频率对应处理结果时,引入了原始单道地震数据、傅里叶变换结果和希尔伯特变换结果,并且希尔伯特变换结果用于对三瞬属性(瞬时振幅、瞬时频率和瞬时相位)进行约束。可见,本实施例加入了对原始单道地震数据的约束以及对三瞬属性的约束,从而能够避免信号在频率域变换时可能发生的畸变现象。
总体来讲,本实施例将地震数据进行基于单道的分频迭代高分辨率处理,可以有效提高地震数据的分辨率并且保持子波运动学特征不变,在拓宽地震数据的高频端的同时也完成了低频端的拓宽处理,使得有效信号的频带得到明显的展宽。具体地,单道的分频迭代高分辨率处理是基于傅里叶变换和希尔伯特变换进行解析函数的构建,可以在不同尺度对数据进行单道单频率的高分辨率处理,避免了由于频谱的拓展造成的波形畸变的产生。因此,本实施例大幅提高了地震资料分辨率,为后续的处理和解释提供技术支撑。
实施例四
参照图5,本实施例在实施例三的基础上增加了构建模块501。构建模块501与处理结果确定模块404连接。该构建模块501设置为构建输出函数K(t)。
图6示出了本发明实施例中构建模块501的结构示意图。参照图6,本实施例构建模块501主要包括第一构建单元601、第二构建单元602、第三构建单元603和输出函数确定单元604。其中,第一构建单元601与输出函数确定单元604连接。第二构建单元602通过第三构建单元603与输出函数确定单元604连接。
具体地,第一构建单元601,设置为构建第一解析函数E(t),并使第一解析函数E(t)满足:E(t)=x(t)+jh(t)。
第二构建单元602,设置为构建第二解析函数Y(t),并使第二解析函数Y(t)满足:Y(t)=Xr(t)+j2πft。
第三构建单元603,设置为根据第二解析函数Y(t)构建第三解析函数Z(t),并使第三解析函数Z(t)满足:
输出函数确定单元604,对第一解析函数E(t)和第三解析函数Z(t)进行乘积运算,得到输出函数K(t),输出函数K(t)满足:
在本实施例中,第一、第二和第三解析函数的构建主要是为了逐层推导出输出函数而构建的,这样可以在加入约束的同时使得输出函数的推导思路清晰。此外,单道单频率的输出函数K(t)是本发明实施例的核心函数,该输出函数通过对输入信号、三瞬属性、振幅和相位进行约束,来避免信号在频率域变换时可能发生的畸变现象。
需要说明的是,实施例三和实施例四中各模块或各单元中的操作的具体细化,可参见上面结合图1至图3、图7a至图12b对本发明方法的说明,在此不再详细赘述。
综上所述,本发明实施例可有效补偿因常规采集仪器的带宽原因造成的低频端缺失,使得地震资料的频带得到明显的均匀展宽,从而大幅提高了地震资料分辨率,并且避免了由于频带的拓展造成的波形畸变并且保持子波运动学特征不变。采用宽频采集仪器无疑会大幅增加地震数据的采集成本,在相当长的时期内低成本的常规采集仪器仍将发挥重要作用。而本发明实施例的处理系统优势非常适用于由常规采集仪器采集到的地震资料的高分辨率处理。
本领域的技术人员应该明白,上述的本发明的各模块或各步骤可以用通用的计算装置来实现,它们可以集中在单个的计算装置上,或者分布在多个计算装置所组成的网络上,可选地,它们可以用计算装置可执行的程序代码来实现,从而,可以将它们存储在存储装置中由计算装置来执行,或者将它们分别制作成各个集成电路模块,或者将它们中的多个模块或步骤制作成单个集成电路模块来实现。这样,本发明不限制于任何特定的硬件和软件结合。
虽然本发明所公开的实施方式如上,但所述的内容只是为了便于理解本发明而采用的实施方式,并非用以限定本发明。任何本发明所属技术领域内的技术人员,在不脱离本发明所公开的精神和范围的前提下,可以在实施的形式上及细节上作任何的修改与变化,但本发明的保护范围,仍须以所附的权利要求书所界定的范围为准。
Claims (10)
1.一种地震数据的处理方法,其特征在于,包括:
获取原始单道地震数据;
对所述原始单道地震数据分别进行傅里叶变换和希尔伯特变换,得到傅里叶变换结果和希尔伯特变换结果;
依次对预设的分频范围内的每个待处理频率,根据所述原始单道地震数据、傅里叶变换结果和希尔伯特变换结果,得到所述待处理频率对应的处理结果;以及
根据所述分频范围内所有待处理频率对应的处理结果,得到输出结果,
其中,根据得到所述待处理频率f对应的处理结果,式中,x(t)为所述原始单道地震数据,Xr(t)为所述傅里叶变换结果的实部,h(t)为所述希尔伯特变换结果。
2.根据权利要求1所述的方法,其特征在于,所述待处理频率f对应的处理结果Kr(t)为构建的输出函数K(t)的实部。
3.根据权利要求2所述的方法,其特征在于,还包括:构建所述输出函数K(t)。
4.根据权利要求3所述的方法,其特征在于,构建所述输出函数K(t),包括:
构建第一解析函数E(t),并使所述第一解析函数E(t)满足:E(t)=x(t)+jh(t);
构建第二解析函数Y(t),并使所述第二解析函数Y(t)满足:Y(t)=Xr(t)+j2πft;
根据所述第二解析函数Y(t)构建第三解析函数Z(t),并使所述第三解析函数Z(t)满足:
对所述第一解析函数E(t)和所述第三解析函数Z(t)进行乘积运算,得到所述输出函数K(t),所述输出函数K(t)满足:
5.根据权利要求1至4中任一项所述的方法,其特征在于,根据所述分频范围内所有待处理频率对应的处理结果,得到输出结果,包括:
对所述分频范围内所有待处理频率对应的处理结果进行求和运算,得到所述输出结果。
6.一种地震数据的处理系统,其特征在于,包括:
数据获取模块,设置为获取原始单道地震数据;
傅里叶变换模块,设置为对所述原始单道地震数据进行傅里叶变换,得到傅里叶变换结果;
希尔伯特变换模块,设置为对所述原始单道地震数据进行希尔伯特变换,得到希尔伯特变换结果;
处理结果确定模块,设置为依次对预设的分频范围内的每个待处理频率,根据所述原始单道地震数据、傅里叶变换结果和希尔伯特变换结果,得到所述待处理频率对应的处理结果;以及
输出结果确定模块,设置为根据所述分频范围内所有待处理频率对应的处理结果,得到输出结果,
其中,所述处理结果确定模块具体设置为根据:得到所述待处理频率f对应的处理结果,式中,x(t)为所述原始单道地震数据,Xr(t)为所述傅里叶变换结果的实部,h(t)为所述希尔伯特变换结果。
7.根据权利要求6所述的系统,其特征在于,所述待处理频率f对应的处理结果Kr(t)为构建的输出函数K(t)的实部。
8.根据权利要求7所述的系统,其特征在于,还包括构建模块,所述构建模块设置为构建所述输出函数K(t)。
9.根据权利要求8所述的系统,其特征在于,所述构建模块包括:
第一构建单元,设置为构建第一解析函数E(t),并使所述第一解析函数E(t)满足:E(t)=x(t)+jh(t);
第二构建单元,设置为构建第二解析函数Y(t),并使所述第二解析函数Y(t)满足:Y(t)=Xr(t)+j2πft;
第三构建单元,设置为根据所述第二解析函数Y(t)构建第三解析函数Z(t),并使所述第三解析函数Z(t)满足:
输出函数确定单元,对所述第一解析函数E(t)和所述第三解析函数Z(t)进行乘积运算,得到所述输出函数K(t),所述输出函数K(t)满足:
10.根据权利要求6至9中任一项所述的系统,其特征在于,所述输出结果确定模块具体设置为:
对所述分频范围内所有待处理频率对应的处理结果进行求和运算,得到所述输出结果。
Priority Applications (5)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610812740.4A CN107807390B (zh) | 2016-09-09 | 2016-09-09 | 地震数据的处理方法及系统 |
CA2974134A CA2974134A1 (en) | 2016-09-09 | 2017-07-19 | Method and system for seismic data processing |
AU2017210665A AU2017210665B2 (en) | 2016-09-09 | 2017-08-07 | Method and system for seismic data processing |
FR1757929A FR3055974B1 (fr) | 2016-09-09 | 2017-08-28 | Procede et systeme de traitement de donnees sismiques |
RU2017131554A RU2751088C2 (ru) | 2016-09-09 | 2017-09-08 | Способ и система обработки сейсмических данных |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610812740.4A CN107807390B (zh) | 2016-09-09 | 2016-09-09 | 地震数据的处理方法及系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107807390A CN107807390A (zh) | 2018-03-16 |
CN107807390B true CN107807390B (zh) | 2019-08-23 |
Family
ID=61387318
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610812740.4A Active CN107807390B (zh) | 2016-09-09 | 2016-09-09 | 地震数据的处理方法及系统 |
Country Status (5)
Country | Link |
---|---|
CN (1) | CN107807390B (zh) |
AU (1) | AU2017210665B2 (zh) |
CA (1) | CA2974134A1 (zh) |
FR (1) | FR3055974B1 (zh) |
RU (1) | RU2751088C2 (zh) |
Families Citing this family (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108802820B (zh) * | 2018-05-28 | 2019-10-11 | 中国石油天然气股份有限公司 | 一种深度域反假频方法、装置及系统 |
CN111142161B (zh) * | 2019-12-31 | 2020-10-09 | 北京东方盈科伟业科技有限公司 | 基于地震数据的复数域地质影像方法及电子设备 |
CN111242198B (zh) * | 2020-01-06 | 2023-05-09 | 中国石油化工股份有限公司 | 基于稀疏表示的道集调谐优化增维方法 |
CN113960671B (zh) * | 2020-07-20 | 2024-04-16 | 中国石油化工股份有限公司 | 频变子波压缩处理方法、装置、计算机设备和存储介质 |
CN111856562B (zh) * | 2020-07-30 | 2022-07-26 | 成都理工大学 | 一种广义高阶同步挤压地震信号时频分解与重构方法 |
CN114114418B (zh) * | 2020-08-28 | 2024-06-18 | 中国石油化工股份有限公司 | 地震信号增益方法、装置、计算机设备和存储介质 |
CN113466936B (zh) * | 2021-07-01 | 2023-07-18 | 中海石油(中国)有限公司上海分公司 | 一种断层阴影区crp道集获取方法、装置、设备及存储介质 |
CN113805225B (zh) * | 2021-08-12 | 2023-09-15 | 中国石油大学(华东) | 一种相位约束的高分辨率地震反演方法及计算机设备 |
CN117233839B (zh) * | 2023-11-10 | 2024-01-26 | 山东科技大学 | 地震数据大地吸收衰减三维空间质控方法、系统以及设备 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102854526A (zh) * | 2011-06-30 | 2013-01-02 | 中国石油化工股份有限公司 | 一种多分量地震资料的处理方法 |
CN103217713A (zh) * | 2012-01-19 | 2013-07-24 | 中国石油化工股份有限公司 | 油气勘探地震速度分析数据优化方法 |
CN103869357A (zh) * | 2012-12-13 | 2014-06-18 | 中国石油集团东方地球物理勘探有限责任公司 | 一种地震数据频谱奇异性的检测方法与装置 |
EP2757391A2 (en) * | 2013-01-22 | 2014-07-23 | Schlumberger Technology B.V. | Automatic processing of ultrasonic data |
CN104459772A (zh) * | 2013-09-25 | 2015-03-25 | 中国石油化工股份有限公司 | 一种地震数字信号的提频方法及装置 |
CN105044769A (zh) * | 2015-06-10 | 2015-11-11 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 提高地震信号的分辨率的方法 |
Family Cites Families (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5852262A (en) * | 1995-09-28 | 1998-12-22 | Magnetic Pulse, Inc. | Acoustic formation logging tool with improved transmitter |
US6131071A (en) * | 1996-12-06 | 2000-10-10 | Bp Amoco Corporation | Spectral decomposition for seismic interpretation |
US5870691A (en) * | 1996-12-06 | 1999-02-09 | Amoco Corporation | Spectral decomposition for seismic interpretation |
RU2119679C1 (ru) * | 1997-12-05 | 1998-09-27 | Станислав Васильевич Васильев | Способ поиска и разведки нефтяных и газовых залежей |
US5995907A (en) * | 1998-02-05 | 1999-11-30 | Geoquest | Seismic signal processing method and apparatus for generating time slice or horizon maps in response to seismic traces and quadrature traces to determine geologic features |
US6931324B2 (en) * | 2003-10-16 | 2005-08-16 | Rdspi, L.P. | Method for determining formation quality factor from seismic data |
US7230879B2 (en) * | 2005-02-12 | 2007-06-12 | Chevron U.S.A. Inc. | Method and apparatus for true relative amplitude correction of seismic data for normal moveout stretch effects |
US7376517B2 (en) * | 2005-05-13 | 2008-05-20 | Chevron U.S.A. Inc. | Method for estimation of interval seismic quality factor |
CN103376464B (zh) * | 2012-04-13 | 2016-04-06 | 中国石油天然气集团公司 | 一种地层品质因子反演方法 |
US10564306B2 (en) * | 2014-08-29 | 2020-02-18 | Pgs Geophysical As | Methods and systems to evaluate noise content in seismic data |
-
2016
- 2016-09-09 CN CN201610812740.4A patent/CN107807390B/zh active Active
-
2017
- 2017-07-19 CA CA2974134A patent/CA2974134A1/en active Pending
- 2017-08-07 AU AU2017210665A patent/AU2017210665B2/en active Active
- 2017-08-28 FR FR1757929A patent/FR3055974B1/fr active Active
- 2017-09-08 RU RU2017131554A patent/RU2751088C2/ru active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102854526A (zh) * | 2011-06-30 | 2013-01-02 | 中国石油化工股份有限公司 | 一种多分量地震资料的处理方法 |
CN103217713A (zh) * | 2012-01-19 | 2013-07-24 | 中国石油化工股份有限公司 | 油气勘探地震速度分析数据优化方法 |
CN103869357A (zh) * | 2012-12-13 | 2014-06-18 | 中国石油集团东方地球物理勘探有限责任公司 | 一种地震数据频谱奇异性的检测方法与装置 |
EP2757391A2 (en) * | 2013-01-22 | 2014-07-23 | Schlumberger Technology B.V. | Automatic processing of ultrasonic data |
CN104459772A (zh) * | 2013-09-25 | 2015-03-25 | 中国石油化工股份有限公司 | 一种地震数字信号的提频方法及装置 |
CN105044769A (zh) * | 2015-06-10 | 2015-11-11 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 提高地震信号的分辨率的方法 |
Non-Patent Citations (1)
Title |
---|
基于快速傅里叶变换的复地震道分析及应用;禾苗 等;《成都航空职业技术学院学报》;20100331(第82期);第61-65页 |
Also Published As
Publication number | Publication date |
---|---|
FR3055974B1 (fr) | 2021-09-03 |
CN107807390A (zh) | 2018-03-16 |
RU2017131554A3 (zh) | 2020-10-15 |
AU2017210665B2 (en) | 2023-04-20 |
AU2017210665A1 (en) | 2018-03-29 |
RU2017131554A (ru) | 2019-03-11 |
RU2751088C2 (ru) | 2021-07-08 |
CA2974134A1 (en) | 2018-03-09 |
FR3055974A1 (fr) | 2018-03-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107807390B (zh) | 地震数据的处理方法及系统 | |
CN102109612B (zh) | 一种地震波吸收衰减补偿方法 | |
CN107942375B (zh) | 一种基于声波方程的非线性优化隐式时空域有限差分数值模拟方法 | |
CN104570110A (zh) | 一种基于纵横波匹配的多分量资料联合速度分析方法 | |
US11209563B2 (en) | Data optimization method and integral prestack depth migration method | |
CN106934183B (zh) | 频散曲线确定方法和装置,及纵横波速度确定方法和装置 | |
CN109633741B (zh) | 基于双凸优化稀疏约束的混合震源数据一次波分离方法 | |
CN105425289A (zh) | 确定低频波阻抗的方法和装置 | |
CN106569262B (zh) | 低频地震数据缺失下的背景速度模型重构方法 | |
CN110058303A (zh) | 声波各向异性逆时偏移混合方法 | |
CN108845357A (zh) | 一种基于同步挤压小波变换估计地层等效品质因子的方法 | |
CN104391324A (zh) | 依赖频率的avo反演前的地震道集动校拉伸校正预处理技术 | |
CN103364826A (zh) | 基于独立分量分析的地震盲源反褶积方法 | |
CA3090760A1 (en) | Systems and methods to enhance 3-d prestack seismic data based on non-linear beamforming in the cross-spread domain | |
CN107229075B (zh) | 深度域地震子波的确定方法和装置 | |
CN104422956B (zh) | 一种基于稀疏脉冲反演的高精度地震谱分解方法 | |
CN109856672B (zh) | 基于深度波数谱的瞬变波包提取方法、存储介质与终端 | |
CN104765063B (zh) | 基于频谱计算吸收衰减属性的油气检测方法及装置 | |
CN106547021B (zh) | 基于单井褶积算法建立初始模型的方法和装置 | |
Zheng et al. | Research on application of Hilbert transform in radar signal simulation | |
CN109270576A (zh) | 多通道三维地震物理模拟数据采集方法及装置 | |
CN114114421B (zh) | 基于深度学习的导向自学习地震数据去噪方法及装置 | |
CN102353991B (zh) | 基于匹配地震子波的物理小波的地震瞬时频率分析方法 | |
CN113917533B (zh) | Ti介质双联动全方位成像的系统性实现方法 | |
Margrave et al. | Improving explicit seismic depth migration with a stabilizing Wiener filter and spatial resampling |
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 |