CN105700014A - 一种基于频域显著性检测的地震属性分析方法 - Google Patents
一种基于频域显著性检测的地震属性分析方法 Download PDFInfo
- Publication number
- CN105700014A CN105700014A CN201610051252.6A CN201610051252A CN105700014A CN 105700014 A CN105700014 A CN 105700014A CN 201610051252 A CN201610051252 A CN 201610051252A CN 105700014 A CN105700014 A CN 105700014A
- Authority
- CN
- China
- Prior art keywords
- frequency domain
- seismic
- significance
- section
- phase spectrum
- 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
- 238000004458 analytical method Methods 0.000 title claims abstract description 18
- 238000001228 spectrum Methods 0.000 claims abstract description 27
- 238000001914 filtration Methods 0.000 claims abstract description 3
- 238000001514 detection method Methods 0.000 claims description 26
- 230000009466 transformation Effects 0.000 claims description 21
- 238000000034 method Methods 0.000 description 14
- 238000004422 calculation algorithm Methods 0.000 description 6
- 238000011160 research Methods 0.000 description 6
- 230000015572 biosynthetic process Effects 0.000 description 3
- 238000011161 development Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 230000000007 visual effect Effects 0.000 description 3
- 230000006835 compression Effects 0.000 description 2
- 238000007906 compression Methods 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 208000035126 Facies Diseases 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000007621 cluster analysis Methods 0.000 description 1
- 230000001427 coherent effect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 239000012530 fluid Substances 0.000 description 1
- 238000012067 mathematical method Methods 0.000 description 1
- 238000003786 synthesis reaction Methods 0.000 description 1
- 238000004454 trace mineral analysis Methods 0.000 description 1
- 230000003945 visual behavior 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/30—Analysis
- G01V1/301—Analysis for determining seismic cross-sections or geostructures
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
一种基于频域显著性检测的地震属性分析方法,旨在将显著性检测应用于地震领域,获得新的且有效的地震属性,其方案为读入地震剖面x(t,q),q为地震剖面的道号;对每一道地震数据做傅里叶变换F,获得地震剖面的频域表示X(f,q);计算地震剖面的频域表示X(f,q)的相位谱P(f,q);利用相位谱信息P(f,q)进行傅里叶反变换F-1,初步获得显著性剖面S(t,q);用低通滤波器g对显著性剖面S(t,q)进行滤波,作为最终的地震显著性属性剖面。
Description
技术领域
本发明属于地震信号处理领域,具体涉及一种基于频域显著性检测的地震属性分析方法。
背景技术
“地震属性”最早可追溯到20世纪60年代,当时国内有很多的译名,如地震特征、地震参数、地震标志等等,直到90年代初才基本统一称为地震属性。从应用地球物理学的角度看,地震属性是地震数据中反映不同地质信息的子集,是刻画、描述地层结构、岩性以及物性等地质信息的地震特征量。
地震属性的发展大致经历了以下三个阶段:
第一阶段是起步阶段,20世纪60年代到70年代末,人们通过各种观察提出了“亮点”、“暗点”和“平点”技术,利用这个技术可以直接进行油气检测,随着数学方法的引入,又提出了瞬时属性和复数道分析技术。然而,这时的属性研究既没有考虑地震资料的运动学、动力学特征,也没有特定的地质含义,仅仅是对地震剖面特征的一种定性描述与分析;
第二阶段是迅速发展阶段,20世纪80年代初期,一方面利用振幅随炮检距变化的规律,另一方面出现了大量的属性定量提取方法,提取出来的属性多达几十种,但是这样提取出来的地震属性没有明确的地质意义,在应用过程中导致了人们对它的不信任;
第三阶段是基本成熟阶段,20世纪90年代初,以相干、倾角、方位角等为代表的一批多维属性分析技术开始出现,地震属性有了更明确的地质意义。这类属性的一个显著特点是能够直观地反映地层的结构性信息,因此,它们的开发与应用,很快得到有关各方的广泛认同。
目前,地震属性分析技术在构造解释、地层岩性解释、储层评价、油藏描述以及油藏流体动态检测等领域得到了广泛应用,在模型正演、相干体技术、聚类分析、地震相分析,多属性综合分析等方面也有了较大的发展,地震属性分析技术在油气勘探开发中发挥着越来越重要的作用。
人类视觉系统是人类认识世界的一个重要生理系统,也被认为是人类已知最强大最复杂的生理系统之一。视觉显著性检测的概念就是在研究人类视觉系统时提出,很早就在生物学上开始研究。基于计算机的显著性检测研究开始于20世纪90年代,处理的对象是图像或者视频。最开始,研究主要集中在如何利用生物学上关于视觉显著性的研究成果,在计算机上建立相似的模型并且采用近似的处理方式来模拟人类的视觉行为。这些模型比较复杂,用计算机实现起来效率不高,检测的效果也比较一般。研究人员逐渐意识到生物学上对视觉显著性的研究是不彻底的,很多的结论具有猜测性。在这样不可靠的理论上建立的模拟算法很难有可靠的效果。逐渐,很多方法突破了严格的生物学模型,模型中的很多部分被简化,采用更为简单直接的计算实现。这些算法在检测效果有了较大幅度的提高。
2007年,Hou等人首次提出基于变换域的显著性检测,提出一种基于傅里叶剩余谱的显著性检测模型(spectralresidual,SR)。SR方法不需要提取特征,也不计算对比度,而是直接对图像做傅里叶变换,并对傅里叶变换的幅度谱进行分析。根据高频抑制原理,人的视觉系统能抑制那些高频出现的特征,同时能对不常出现的特征反应敏感。因此该方法通过计算图像的傅里叶变换,并用对数表示其幅度谱,然后用对数幅度谱减去经过高斯低通滤波器后的对数幅度谱,经过高斯低通滤波器滤波后能抑制图像的高频成分,因此能保留不常出现的显著部分。SR算法的计算十分简单,只需要用到图像处理常用的傅里叶变换和反变换,但是其显著性检测效果一般,只能处理简单背景图像的显著性检测。
正是由于基于频域分析的显著性检测模型计算简单等特点,越来越多的基于变换域的显著性检测模型被提出。Guo等人在SR的基础上,用实验证明了图像的相位谱特征是计算显著性的关键特征,并提出了基于傅里叶变换相位谱的显著性检测模型(phasespectrumofFouriertransform,PFT),该方法对图像的亮度特征进行傅里叶变换,并用其相位谱的反变换结果与合适大小的二维高斯滤波器滤波后的灰度图作为其最终的显著图。在此基础上,Guo等人又提出一种基于四元数傅里叶变换相位谱的显著性检测模型(phasespectrumofquaternionFouriertransform,PQFT),该模型不仅仅使用图像的亮度特征,还计算了图像的两个对立颜色特征和视频序列中的运动特征。用视频图像的四个近似独立的特征通道构成图像四元数表示,并用图像四元傅里叶变换来代替传统的傅里叶变换。此外,Guo等人还将PQFT的多分辨率显著性检测算法应用于视频图像压缩,并能获得比传统编码算法更高的压缩率。该算法由于引入了视频的运动特征,因此其不仅能够检测出图像的显著性区域,还能检测到视频序列的显著性区域。该模型对图像分辨率的鲁棒性较高,它不仅可以检测粗糙图像显著性,还能检测高分辨率的显著性区域。
发明内容
本发明提供了一种基于频域显著性检测的地震属性分析方法,旨在将显著性检测应用于地震领域,获得新的且有效的地震属性。
为了解决上述技术问题,达到上述目的,本发明采用如下技术方案:
一种基于频域显著性检测的地震属性分析方法,其特征在于,包括以下步骤:
步骤1:读入地震剖面x(t,q),q为地震剖面的道号;
步骤2:对每一道地震数据做傅里叶变换F,获得地震剖面的频域表示X(f,q);
步骤3:计算地震剖面的频域表示X(f,q)的相位谱P(f,q);
步骤4:利用相位谱信息P(f,q)进行傅里叶反变换F-1,初步获得显著性剖面S(t,q);
步骤5:用低通滤波器g对显著性剖面S(t,q)进行滤波,作为最终的地震显著性属性剖面。
上述技术方案中,其中所述步骤2涉及傅里叶变换如下:
设傅里叶变换算子为F,对地震剖面x(t,q)的每一道数据做傅里叶变换,获得频域表示X(f,q)为
式中,q表示地震剖面的道号,f是频率域变量,t是时间变量,j是虚数单位,e是自然常数。
上述技术方案中,其中所述步骤3涉及频域相位谱如下:
对于频域表示X(f,q),其相位谱P(f,q)为
式中,Im{·}表示虚部,Re{·}表示实部,arctan表示反正切函数。
上述技术方案中,其中所述步骤4涉及利用相位谱信息进行傅里叶反变换如下:
已知相位谱P(f,q),显著性剖面S(t,q)表示为
S(t,q)=|F-1{ejP(f,q)}|(8)
式中,q表示地震剖面的道号,t是时间变量,f是频率域变量,j是虚数单位,e是自然常数,F-1表示傅里叶反变换,即
上述技术方案中,其中所述步骤5涉及滤波过程,如下:
对二维显著性剖面S(t,q),q表示地震剖面的道号,采用低通滤波器g,进行滤波后地震显著性剖面FS(t,q)表示为
FS(t,q)=g*|S(t,q)|2(10)。
与现有技术相比,本发明具有如下优势:
本发明将利用傅里叶变换相位谱提取地震显著性属性,有效地将显著性检测应用于地震领域,为地震解释提供了新的地震属性分析方法,可以有效地进行储层预测和油气解释。
附图说明
图1为方法流程图;
图2为地震剖面;
图3为地震剖面的频域表示;
图4为地震剖面的频域表示的相位谱;
图5为显著性剖面;
图6为最终的地震显著性属性剖面。
具体实施方式
下面结合附图对本发明进行具体说明:
步骤1:读入地震剖面x(t,q)(如图2所示),q为地震剖面的道号;
步骤2:对每一道地震数据做傅里叶变换F,获得地震剖面的频域表示X(f,q)(如图3所示)
步骤3:计算地震剖面的频域表示X(f,q)的相位谱P(f,q)(如图4所示)
步骤4:利用相位谱信息P(f,q)进行傅里叶反变换F-1,初步获得显著性剖面S(t,q)(如图5所示)
步骤5:用低通滤波器g对显著性剖面S(t,q)进行滤波,滤波结果为
FS(t,q)=g*|S(t,q)|2(14)
FS(t,q)即最终的地震显著性属性剖面(如图6所示,该显著性剖面可以展示地下储层分布信息,专业地震解释人员可以利用该剖面确定储层区域)。
Claims (5)
1.一种基于频域显著性检测的地震属性分析方法,其特征在于,包括以下步骤:
步骤1:读入地震剖面x(t,q),q为地震剖面的道号;
步骤2:对每一道地震数据做傅里叶变换F,获得地震剖面的频域表示X(f,q);
步骤3:计算地震剖面的频域表示X(f,q)的相位谱P(f,q);
步骤4:利用相位谱信息P(f,q)进行傅里叶反变换F-1,获得显著性剖面S(t,q);
步骤5:用低通滤波器g对显著性剖面S(t,q)进行滤波,作为最终的地震显著性属性剖面。
2.根据权利要求1所述的一种基于频域显著性检测的地震属性分析方法,其特征在于,其中所述步骤2涉及傅里叶变换如下:
傅里叶变换算子为F,对地震剖面x(t,q)的每一道数据做傅里叶变换,获得频域表示X(f,q)
式中,q表示地震剖面的道号,f是频率域变量,t是时间变量,j是虚数单位,e是自然常数。
3.根据权利要求1所述的一种基于频域显著性检测的地震属性分析方法,其特征在于,其中所述步骤3涉及频域相位谱如下:
对于频域表示X(f,q),其相位谱P(f,q)为
式中,Im{·}表示虚部,Re{·}表示实部,arctan表示反正切函数。
4.根据权利要求1所述的一种基于频域显著性检测的地震属性分析方法,其特征在于,其中所述步骤4涉及利用相位谱信息进行傅里叶反变换如下:
已知相位谱P(f,q),显著性剖面S(t,q)表示为
S(t,q)=|F-1{ejP(f,q)}|(3)
式中,q表示地震剖面的道号,t是时间变量,f是频率域变量,j是虚数单位,e是自然常数,F-1表示傅里叶反变换,即
5.根据权利要求1所述的一种基于频域显著性检测的地震属性分析方法,其特征在于,步骤5涉及滤波过程,如下:
对显著性剖面S(t,q),q表示地震剖面的道号,采用低通滤波器g,进行滤波后地震显著性剖面FS(t,q)表示为
FS(t,q)=g*|S(t,q)|2(5)。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610051252.6A CN105700014B (zh) | 2016-01-26 | 2016-01-26 | 一种基于频域显著性检测的地震属性分析方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610051252.6A CN105700014B (zh) | 2016-01-26 | 2016-01-26 | 一种基于频域显著性检测的地震属性分析方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105700014A true CN105700014A (zh) | 2016-06-22 |
CN105700014B CN105700014B (zh) | 2018-05-15 |
Family
ID=56228590
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610051252.6A Active CN105700014B (zh) | 2016-01-26 | 2016-01-26 | 一种基于频域显著性检测的地震属性分析方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105700014B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109164483A (zh) * | 2018-08-29 | 2019-01-08 | 中国科学院地球化学研究所 | 多分量地震数据矢量去噪方法及多分量地震数据矢量去噪装置 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102279416A (zh) * | 2010-12-30 | 2011-12-14 | 北京诺克斯达石油科技有限公司 | 远近偏移距叠加剖面相位频率的匹配方法 |
CN102305940A (zh) * | 2011-05-24 | 2012-01-04 | 中国石油集团川庆钻探工程有限公司 | 流体因子提取方法 |
WO2012170090A1 (en) * | 2011-06-08 | 2012-12-13 | Chevron U.S.A. Inc. | System and method for data inversion with phase unwrapping |
CN104966286A (zh) * | 2015-06-04 | 2015-10-07 | 电子科技大学 | 一种3d视频显著性检测方法 |
CN105093325A (zh) * | 2014-05-22 | 2015-11-25 | 中国石油化工股份有限公司 | 一种定量的提频方法 |
CN105116443A (zh) * | 2015-07-27 | 2015-12-02 | 中国石油天然气集团公司 | 一种低频信号的能量补偿方法及装置 |
-
2016
- 2016-01-26 CN CN201610051252.6A patent/CN105700014B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102279416A (zh) * | 2010-12-30 | 2011-12-14 | 北京诺克斯达石油科技有限公司 | 远近偏移距叠加剖面相位频率的匹配方法 |
CN102305940A (zh) * | 2011-05-24 | 2012-01-04 | 中国石油集团川庆钻探工程有限公司 | 流体因子提取方法 |
WO2012170090A1 (en) * | 2011-06-08 | 2012-12-13 | Chevron U.S.A. Inc. | System and method for data inversion with phase unwrapping |
CN105093325A (zh) * | 2014-05-22 | 2015-11-25 | 中国石油化工股份有限公司 | 一种定量的提频方法 |
CN104966286A (zh) * | 2015-06-04 | 2015-10-07 | 电子科技大学 | 一种3d视频显著性检测方法 |
CN105116443A (zh) * | 2015-07-27 | 2015-12-02 | 中国石油天然气集团公司 | 一种低频信号的能量补偿方法及装置 |
Non-Patent Citations (1)
Title |
---|
黄虹: "自然图像统计显著图模型与质量盲评估测度", 《中国优秀硕士学位论文全文数据库 信息科技辑》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109164483A (zh) * | 2018-08-29 | 2019-01-08 | 中国科学院地球化学研究所 | 多分量地震数据矢量去噪方法及多分量地震数据矢量去噪装置 |
US11467298B2 (en) | 2018-08-29 | 2022-10-11 | Institute Of Geochemistry, Chinese Academy Of Sciences | Vector denoising method and device for multicomponent seismic data |
Also Published As
Publication number | Publication date |
---|---|
CN105700014B (zh) | 2018-05-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Lu et al. | Seismic spectral decomposition using deconvolutive short-time Fourier transform spectrogram | |
CN107817527B (zh) | 基于块稀疏压缩感知的沙漠地震勘探随机噪声压制方法 | |
Pan et al. | Identification of stratigraphic formation interfaces using wavelet and Fourier transforms | |
Li et al. | Time-synchroextracting general chirplet transform for seismic time–frequency analysis | |
CN103995289A (zh) | 基于时频谱模拟的时变混合相位地震子波提取方法 | |
CN103728663B (zh) | 一种时频分析方法 | |
CN102692647B (zh) | 一种高时间分辨率的地层含油气性预测方法 | |
CN101632587A (zh) | 一种基于小波变换的步态信号提取方法 | |
CN104280765A (zh) | 基于变子波反射系数反演的地震高分辨处理方法 | |
CN106680874A (zh) | 基于波形形态特征稀疏化建模的谐波噪声压制方法 | |
CN107255831A (zh) | 一种叠前频散属性的提取方法 | |
CN107356964A (zh) | S变换域基于变分原理的q值估计与补偿方法 | |
Soliman et al. | Application of wavelet transform to the analysis of pressure-transient data | |
Pan et al. | Multiscale frequency-domain seismic inversion for fracture weakness | |
CN106707334A (zh) | 一种提高地震资料分辨率的方法 | |
CN107632326B (zh) | 地球物理信号去噪方法 | |
CN105334532A (zh) | 一种地震子波估计方法 | |
CN109655883A (zh) | 一种针对目标的地震分频方法及系统 | |
CN104991280B (zh) | 一种基于微观孔喉结构指数的射孔层段选择方法及装置 | |
CN105700014B (zh) | 一种基于频域显著性检测的地震属性分析方法 | |
CN102928872B (zh) | 一种提高滩坝砂储层地震识别性能与描述精度的方法 | |
CN105572737B (zh) | 一种基于分数域显著性检测的地震属性分析方法 | |
Shang et al. | Automatic microseismic events detection using morphological multiscale top-hat transformation | |
Xudong et al. | Pre-stack gather optimization technology based on an improved bidimensional empirical mode decomposition method | |
CN106291675A (zh) | 一种基于基追踪技术的地震数据重构方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |