CN107422381A - 一种基于eemd‑ica的地震低频信息流体预测方法 - Google Patents
一种基于eemd‑ica的地震低频信息流体预测方法 Download PDFInfo
- Publication number
- CN107422381A CN107422381A CN201710842580.2A CN201710842580A CN107422381A CN 107422381 A CN107422381 A CN 107422381A CN 201710842580 A CN201710842580 A CN 201710842580A CN 107422381 A CN107422381 A CN 107422381A
- Authority
- CN
- China
- Prior art keywords
- mrow
- signal
- earthquake
- imf
- eemd
- 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 29
- 239000012530 fluid Substances 0.000 title claims abstract description 15
- 238000004458 analytical method Methods 0.000 claims abstract description 35
- 238000000354 decomposition reaction Methods 0.000 claims abstract description 15
- 238000013519 translation Methods 0.000 claims abstract description 9
- 238000002360 preparation method Methods 0.000 claims abstract description 7
- 230000035939 shock Effects 0.000 claims abstract description 7
- 238000010183 spectrum analysis Methods 0.000 claims abstract description 7
- 238000000513 principal component analysis Methods 0.000 claims abstract description 5
- 230000009467 reduction Effects 0.000 claims abstract description 4
- 239000002131 composite material Substances 0.000 claims description 4
- 238000002156 mixing Methods 0.000 claims description 4
- 230000008901 benefit Effects 0.000 abstract description 5
- 238000005516 engineering process Methods 0.000 description 11
- 230000009466 transformation Effects 0.000 description 11
- 239000003921 oil Substances 0.000 description 6
- 238000000605 extraction Methods 0.000 description 5
- 238000012545 processing Methods 0.000 description 5
- 230000007547 defect Effects 0.000 description 4
- 238000001514 detection method Methods 0.000 description 4
- 230000003044 adaptive effect Effects 0.000 description 3
- 230000008859 change Effects 0.000 description 3
- 238000009826 distribution Methods 0.000 description 3
- 238000003384 imaging method Methods 0.000 description 3
- 238000001228 spectrum Methods 0.000 description 3
- 230000007812 deficiency Effects 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 230000004807 localization Effects 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 239000004215 Carbon black (E152) Substances 0.000 description 1
- 206010044565 Tremor Diseases 0.000 description 1
- 230000002159 abnormal effect Effects 0.000 description 1
- 238000010521 absorption reaction Methods 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 238000007405 data analysis Methods 0.000 description 1
- 238000005553 drilling Methods 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 229930195733 hydrocarbon Natural products 0.000 description 1
- 150000002430 hydrocarbons Chemical class 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 239000010742 number 1 fuel oil Substances 0.000 description 1
- 230000010355 oscillation Effects 0.000 description 1
- 230000003595 spectral 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/30—Analysis
- G01V1/306—Analysis for determining physical properties of the subsurface, e.g. impedance, porosity or attenuation profiles
-
- 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/62—Physical property of subsurface
- G01V2210/624—Reservoir parameters
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
本发明公开了一种基于EEMD‑ICA的地震低频信息流体预测方法,包括:输入地震叠后单道CDP信号,进行EEMD经验模态分解,得到IMF分量;对IMF分量整体平均;识别出存在模态混叠的IMF分量信号;得到模态混叠分量所有的平移信号;对平移信号利用主成分分析进行降维;对PCA降维处理后得到的信号,通过Fast‑ICA计算得到解混后的分量信号;对得到信号分别进行频谱分析,并与原信号进行相似系数计算;对地震CDP数据按照道号,依次进行循环运算;对地震数据低频信息进行强振幅异常分析。本发明的优点在于:有效解决了EEMD的随机性、不稳定性、含噪性的问题,能够得到更加精确的提取地震信号的频率域信息,从而使得低频信息的获取更为准确,提高流体预测的精度。
Description
技术领域
本发明涉及油气地震勘探领域领域,特别涉及一种基于EEMD-ICA的地震低频信息流体预测方法。
背景技术
Chapman(2006)等指出在低频范围内,油气储层显示出高异常。Goloshubin利用实际地震数据证实富含油储层使得低频范围内的反射系数显著的增加,拓展了有效地震带宽到更低频端,并表明利用地震波的低频信息预测流体的潜力。
近年来,以瞬时谱分析为基础,利用低频地震信号的振幅(或能量)变化检测储层及识别油气方面出现了一些新的技术进展和成功的应用成果,比较典型的有利用低频阴影、低频振幅异常、瞬时谱能量吸收分析技术、成像属性分析技术等直接检测油气;另外,利用叠前资料的低频谱能量分析进行油气检测已成为一个新的发展趋势。地震波低频信息的潜在意义和应用价值已经越来越受到地球物理勘探界的重视。
低频振幅异常(Low frequency anomalies,LFA)是利用地震波中的低频成分信息直接指示油气储层的分布,也是HCI中的一种油气识别标志。Goloshubin通过实际资料处理证实,反射地震波中的低频成分包含了与油气储层有关的极其重要的信息,对于油气储层检测和指示显示了惊人的成像能力,与钻井结果极其吻合,并指出了常规资料处理中忽视地震信号低频信息的局限性。其它学者也展示了低频振幅异常检测油气储层的成功实例。
挖掘地震资料低频信息的方法有低通滤波、带通滤波、具有时频分辨率的时频分析方法,后者具有瞬时定位分析能力,用于地震资料分析和解释时,被从事地震资料处理和解释的地球物理学者们广泛接受。
地震数据时频分析技术是一种基于频谱分解的地震成像解释方法。时频分析技术通过相关数学变换,将时域信号转换成时频域的信息。时频分析技术可以直接提取地震数据低频信息的信号处理技术。
时频分析技术是提取地震低频信息最直接有效的方法。前人在利用时频分析技术提取地震低频信息方面已经开展了大量工作。
Gabor于1946年引进窗口傅里叶变换改进了傅里叶变换在时-频局部分析方面的不足,但是由于其时-频局部化窗口大小固定不变,这与时-频局部化分析中要求:“反映信号高频成分需要窄的时间窗,反映信号低频成分需要宽的时间窗”不一致。因此国外内学者相继发展了小波变换分析方法用于时-频局部分析。小波分析作为傅里叶变换分析的发展,既保留了傅里叶变换分析的优点,又弥补了傅里叶变换分析的不足。但是小波分析只对信号的低频端信息成分做进一步的分解,而对高频成分,即信号的细节部分不再继续分解,所以小波分析能够很好地表征一大类以低频信息为主要成分的信号,但是不能很好地分解和表示大量细节的信息。
为了克服小波变换分析的不足,Stockwell(1996)提出S变换分析方法用于局部时-频精细分析。S变换的窗函数不仅能随着频率的增加而自适应地减小分析时窗,并且时窗具有一定的自适应性、不需要满足小波容许性条件,还能与傅里叶变换计算建立直接联系。
针对S变换的小波函数是固定的,不能随具体的应用而调整,缺乏灵活性,高静怀、贺振华、陈学华等相继提出了小波函数变化的广义S变换。广义S变换不仅具有小波分析的特点,并且对信号的细节,及信号的高频部分也能进行很好的时-频分析。从时频分析技术的发展及其各自时频分析技术的优点和时频聚集性能认识到提取地震波高频信息的最佳方法是改进的S变换分析,基于时频聚集性度量准则的一种自适应S变换进行地震数据时频分析。
近年来,经验模态分解EMD方法被逐渐应用于地震信号分析,并用于地震低频信息提取方面。经验模态分解法(EMD)是由Huang等提出的一种非平稳信号分析方法,该方法通过提取复杂信号在每一个时刻局部的振荡模式,按由高到低的自适应频率分解模式寻找信号内蕴的高频信息,进而分解得到若干个平稳信号分量,即模态函数分量(IMF)。EMD分解法其应用领域已经遍及地震、雷达和语音信号处理及图像分析等各个面。
因为EMD的自适应性,缺乏约束条件,使其不可避免的存在缺陷,这种缺陷称为"模态混叠"。任何信号我们都可以看成其由若干个固有模态函数(IMF)组成的,一个模态描述一个单一的震动状态,而如果IMF之间相互重叠,则形成复合信号。在经验模态分析过程中,期望将这些单一的模态干净的分离出来,传统的经验模态分解方法(EMD)由于算法本身的局限性,在分离出来的平稳信号中会包含多个模态,从而造成模态混叠,其结果会造成频谱分析的错误。针对EMD的这一缺陷,Wu和Huang提出了EEMD方法,利用高斯白噪声迭代求取各阶模态函数,EEMD方法在每次迭代时引入噪声,能够部分解决模态混叠现象。
目前,现有的低频信息获取方法主要有傅里叶变换、小波变换、S变换,其中S变换是介于STFT和小波变换之间的一种时频分析方法,它吸收了短时傅立叶变换的思想,引进小波的多分辨分析,改进了短时傅立叶变换中时频分辨率固定的问题,其窗函数与频率(即尺度)自适应地成反比变化,低频部分具有较高的频率分辨率而高频部分具有较高的时间分辨率,满足了地震信号的特征。这种方法中,时间窗口长度取决于频率,S变换的小波基函数是固定的。在实际地震信号处理中,时频分布的特性与地震信号本身及地震子波都有关系,S变换中小波基函数的固定不能满足实际地震数据处理的需求。
对数据分析的最终目标不是找到数据的数学特性,而是挖掘隐藏在数据中的物理表现及其意义。因此,对于非平稳和非线性数据,尤其是地震信号,自适应性是非常重要的。因此,近年来,引入EMD(经验模态分解)方法来进行信号分解,EMD方法不需要预先设定任何基函数,它利用数据自身的时间尺度特征来进行信号的分解。与短时傅立叶变换、小波分解等方法相比,这种方法是直观的、直接的、后验的和自适应的。但是,EMD分解过程中,最大的缺点且是不可避免的就是会产生模态混叠效应,即某个分解的IMF分量信号中包含其他IMF分量的信号。因此,前人将该方法进行了改进,形成了EEMD方法,该方法能够较好的解决模态混叠效应,但是在分解过程中加入随机高斯白噪,造成分解得到的IMF分量存在随机性,且信号中包含了白噪,并不能真实提取地震的低频信息。
发明内容
本发明针对现有技术的缺陷,提供了一种基于EEMD-ICA的地震低频信息流体预测方法,能有效的解决上述现有技术存在的问题。
为了实现以上发明目的,本发明采取的技术方案如下:
一种基于EEMD-ICA的地震低频信息流体预测方法,包括以下步骤:
步骤一:输入地震叠后单道CDP信号,进行EEMD经验模态分解,得到IMF分量;将待处理地震单道CDP信号中加入等长度不等幅高斯白噪声,对复合信号EEMD分解,重复操作k次,得到IMF分量cik和余项rik;
步骤二:对IMF分量整体平均;
步骤三:令步骤二对原始地震单道CDP信号s(t)进行EEMD分解得到IMF分量结果,所有的IMF分量为f1(t)、f2(t)至fn(t),对所有IMF分量信号进行频谱分析,识别出存在模态混叠的IMF分量信号fk(t);
步骤四:对模态混叠的IMF分量fk(t)分别平移τm个位置,得到模态混叠分量fk(t)所有的平移信号分别为fk1(t)、fk2(t)至fkm(t);
步骤五:对步骤四得到的fk1(t)、fk2(t)至fkm(t),共计m个平移信号;对这m个平移信号利用主成分分析进行降维;
步骤六:对PCA降维处理后得到的信号,通过Fast-ICA计算得到gk1(t)、gk2(t)至gkn(t),共计n个信号,得到的信号即为解混后的分量信号;
步骤七:对步骤六得到n个信号分别进行频谱分析,并与原信号进行相似系数计算,去除噪声信号与非相似性冗余信号,保留相似系数最大的单个低频信号分量;
执行公式:
式中:s(t)为源地震CDP信号;g(t)为ICA解混后得到的分量信号;t为CDP信号时间采样点个数;
步骤七:按照步骤一到步骤六,对地震CDP数据按照道号,依次进行循环运算;
步骤八:对得到地震数据低频信息进行强振幅异常分析,强振幅异常区即为流体分布区域。
进一步地,所述步骤二对IMF分量整体平均执行公式:
式中,n为EMD分解的次数,i为信号的采样点个数。
与现有技术相比本发明的优点在于:在EMD、EEMD的算法基础上,利用PCA-ICA的方法进行了IMF信号分量的解混,去掉白噪信号与冗余信号。有效解决了EEMD的随机性、不稳定性、含噪性的问题,能够得到更加精确的提取地震信号的频率域信息,从而使得低频信息的获取更为准确,提高流体预测的精度。
具体实施方式
为使本发明的目的、技术方案及优点更加清楚明白,以下举实施例,对本发明做进一步详细说明。
一种基于EEMD-ICA的地震低频信息流体预测方法,包括以下步骤:
步骤一:输入地震叠后单道CDP信号,进行EEMD经验模态分解,得到IMF分量;将待处理地震单道CDP信号中加入等长度不等幅高斯白噪声,对复合信号EEMD分解,重复操作k次,得到IMF分量cik和余项rik;
步骤二:对IMF分量整体平均,执行公式:
式中,n为EMD分解的次数,i为信号的采样点个数;
步骤三:令步骤二对原始地震单道CDP信号s(t)进行EEMD分解得到IMF分量结果,所有的IMF分量为f1(t)、f2(t)至fn(t),对所有IMF分量信号进行频谱分析,识别出存在模态混叠(即多个频率信号混合)的IMF分量信号fk(t)。
步骤四:由于单个模态混叠分量fk(t)无法利用ICA算法分离出混叠的多个分量。同时,由于时间延迟不影响信号的周期性,所以对模态混叠的IMF分量fk(t)分别平移τm个位置,得到模态混叠分量fk(t)所有的平移信号分别为fk1(t)、fk2(t)、…fkm(t)。
步骤五:对步骤四得到的fk1(t)、fk2(t)、…fkm(t),共计m个平移信号。对这m个平移信号利用主成分分析(Principal Component Analysis,PCA)进行降维。
步骤六:对PCA降维处理后得到的信号,通过Fast-ICA计算得到gk1(t)、gk2(t)、…gkn(t),共计n个信号,得到的信号即为解混后的分量信号;
步骤七:对步骤六得到n个信号分别进行频谱分析,并与原信号进行相似系数计算,去除噪声信号与非相似性冗余信号,保留相似系数最大的单个低频信号分量;
定义信号相似系数如下:s(t)为源地震CDP信号,g(t)为ICA解混后得到的分量信号,为步骤六得到的gk1(t)、gk2(t)、…gkn(t),t为CDP信号时间采样点个数:
执行公式:
步骤七:按照步骤一到步骤六,对地震CDP数据按照道号,依次进行循环运算;
步骤八:对得到地震数据低频信息进行强振幅异常分析,强振幅异常区即为流体分布区域。
本领域的普通技术人员将会意识到,这里所述的实施例是为了帮助读者理解本发明的实施方法,应被理解为本发明的保护范围并不局限于这样的特别陈述和实施例。本领域的普通技术人员可以根据本发明公开的这些技术启示做出各种不脱离本发明实质的其它各种具体变形和组合,这些变形和组合仍然在本发明的保护范围内。
Claims (2)
1.一种基于EEMD-ICA的地震低频信息流体预测方法,其特征在于包括以下步骤:
步骤一:输入地震叠后单道CDP信号,进行EEMD经验模态分解,得到IMF分量;将待处理地震单道CDP信号中加入等长度不等幅高斯白噪声,对复合信号EEMD分解,重复操作k次,得到IMF分量cik和余项rik;
步骤二:对IMF分量整体平均;
步骤三:令步骤二对原始地震单道CDP信号s(t)进行EEMD分解得到IMF分量结果,所有的IMF分量为f1(t)、f2(t)至fn(t),对所有IMF分量信号进行频谱分析,识别出存在模态混叠的IMF分量信号fk(t);
步骤四:对模态混叠的IMF分量fk(t)分别平移τm个位置,得到模态混叠分量fk(t)所有的平移信号分别为fk1(t)、fk2(t)至fkm(t);
步骤五:对步骤四得到的fk1(t)、fk2(t)至fkm(t),共计m个平移信号;对这m个平移信号利用主成分分析进行降维;
步骤六:对PCA降维处理后得到的信号,通过Fast-ICA计算得到gk1(t)、gk2(t)至gkn(t),共计n个信号,得到的信号即为解混后的分量信号;
步骤七:对步骤六得到n个信号分别进行频谱分析,并与原信号进行相似系数计算,去除噪声信号与非相似性冗余信号,保留相似系数最大的单个低频信号分量;
执行公式:
<mrow>
<mi>&sigma;</mi>
<mo>=</mo>
<mi>&sigma;</mi>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>,</mo>
<mi>g</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</munderover>
<mi>s</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mi>g</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
<msqrt>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</munderover>
<msup>
<mi>s</mi>
<mn>2</mn>
</msup>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</munderover>
<msup>
<mi>g</mi>
<mn>2</mn>
</msup>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
</msqrt>
</mfrac>
<mo>;</mo>
</mrow>
式中:s(t)为源地震CDP信号;g(t)为ICA解混后得到的分量信号;t为CDP信号时间采样点个数;
步骤七:按照步骤一到步骤六,对地震CDP数据按照道号,依次进行循环运算;
步骤八:对得到地震数据低频信息进行强振幅异常分析,强振幅异常区即为流体分布区域。
2.根据权利要求1所述的一种基于EEMD-ICA的地震低频信息流体预测方法,其特征在于:所述步骤二对IMF分量整体平均执行公式:
<mrow>
<msub>
<mi>IMF</mi>
<mrow>
<mi>i</mi>
<mi>k</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mi>n</mi>
</mfrac>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>k</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</munderover>
<msub>
<mi>IMF</mi>
<mrow>
<mi>i</mi>
<mi>k</mi>
</mrow>
</msub>
<mo>;</mo>
</mrow>
<mrow>
<msub>
<mi>r</mi>
<mi>i</mi>
</msub>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mi>n</mi>
</mfrac>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>k</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</munderover>
<msub>
<mi>r</mi>
<mrow>
<mi>i</mi>
<mi>k</mi>
</mrow>
</msub>
<mo>;</mo>
</mrow>
式中,n为EMD分解的次数,i为信号的采样点个数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710842580.2A CN107422381B (zh) | 2017-09-18 | 2017-09-18 | 一种基于eemd-ica的地震低频信息流体预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710842580.2A CN107422381B (zh) | 2017-09-18 | 2017-09-18 | 一种基于eemd-ica的地震低频信息流体预测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107422381A true CN107422381A (zh) | 2017-12-01 |
CN107422381B CN107422381B (zh) | 2019-07-02 |
Family
ID=60432604
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710842580.2A Active CN107422381B (zh) | 2017-09-18 | 2017-09-18 | 一种基于eemd-ica的地震低频信息流体预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107422381B (zh) |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108152855A (zh) * | 2017-12-14 | 2018-06-12 | 西南石油大学 | 一种基于eemd-svd的地震流体识别方法 |
CN109116424A (zh) * | 2018-10-11 | 2019-01-01 | 中国石油天然气股份有限公司 | 地震波数据的低波数噪音分离方法、装置及存储介质 |
CN109258509A (zh) * | 2018-11-16 | 2019-01-25 | 太原理工大学 | 一种生猪异常声音智能监测系统与方法 |
CN109740453A (zh) * | 2018-12-19 | 2019-05-10 | 吉林大学 | 一种基于小波变换的卫星磁场数据地震前兆异常提取方法 |
CN110688981A (zh) * | 2019-10-14 | 2020-01-14 | 江西理工大学 | 一种振动信号去噪的模态混叠消除方法 |
CN111859274A (zh) * | 2020-07-03 | 2020-10-30 | 自然资源部第一海洋研究所 | 一种用于度量复杂动力系统本征因果关系的方法 |
CN112464855A (zh) * | 2020-12-09 | 2021-03-09 | 中国石油天然气集团有限公司 | 基于eemd的随钻泥浆正脉冲信号处理方法及装置 |
CN112711070A (zh) * | 2019-10-24 | 2021-04-27 | 中国石油化工股份有限公司 | 一种基于地震信号分解的油气检测方法及装置 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102855408A (zh) * | 2012-09-18 | 2013-01-02 | 福州大学 | 基于ica的改进emd过程中imf判定方法 |
US20130182537A1 (en) * | 2012-01-12 | 2013-07-18 | Cggveritas Services Sa | Device and method for estimating time-shifts |
CN105044777A (zh) * | 2015-07-01 | 2015-11-11 | 中国石油天然气股份有限公司 | 基于经验模态分解检测地震标志层强反射振幅消除的方法 |
CN106898361A (zh) * | 2017-03-16 | 2017-06-27 | 杭州电子科技大学 | 基于反馈变分模式分解的单通道盲源分离方法 |
CN107024718A (zh) * | 2017-05-31 | 2017-08-08 | 西南石油大学 | 基于ceemd‑spwvd时频谱分析的叠后地震流体预测方法 |
CN107102356A (zh) * | 2017-06-02 | 2017-08-29 | 成都理工大学 | 基于ceemd的地震信号高分辨率处理方法 |
-
2017
- 2017-09-18 CN CN201710842580.2A patent/CN107422381B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20130182537A1 (en) * | 2012-01-12 | 2013-07-18 | Cggveritas Services Sa | Device and method for estimating time-shifts |
CN102855408A (zh) * | 2012-09-18 | 2013-01-02 | 福州大学 | 基于ica的改进emd过程中imf判定方法 |
CN105044777A (zh) * | 2015-07-01 | 2015-11-11 | 中国石油天然气股份有限公司 | 基于经验模态分解检测地震标志层强反射振幅消除的方法 |
CN106898361A (zh) * | 2017-03-16 | 2017-06-27 | 杭州电子科技大学 | 基于反馈变分模式分解的单通道盲源分离方法 |
CN107024718A (zh) * | 2017-05-31 | 2017-08-08 | 西南石油大学 | 基于ceemd‑spwvd时频谱分析的叠后地震流体预测方法 |
CN107102356A (zh) * | 2017-06-02 | 2017-08-29 | 成都理工大学 | 基于ceemd的地震信号高分辨率处理方法 |
Non-Patent Citations (2)
Title |
---|
赵迎等: "《CEEMD与小波变换联合去噪方法研究》", 《地球物理学进展》 * |
黄艳波: "《单通道盲源分离算法的研究》", 《杭州电子科技大学硕士学位论文》 * |
Cited By (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108152855A (zh) * | 2017-12-14 | 2018-06-12 | 西南石油大学 | 一种基于eemd-svd的地震流体识别方法 |
CN109116424B (zh) * | 2018-10-11 | 2020-06-09 | 中国石油天然气股份有限公司 | 地震波数据的低波数噪音分离方法、装置及存储介质 |
CN109116424A (zh) * | 2018-10-11 | 2019-01-01 | 中国石油天然气股份有限公司 | 地震波数据的低波数噪音分离方法、装置及存储介质 |
CN109258509A (zh) * | 2018-11-16 | 2019-01-25 | 太原理工大学 | 一种生猪异常声音智能监测系统与方法 |
CN109258509B (zh) * | 2018-11-16 | 2023-05-02 | 太原理工大学 | 一种生猪异常声音智能监测系统与方法 |
CN109740453B (zh) * | 2018-12-19 | 2022-03-29 | 吉林大学 | 一种基于小波变换的卫星磁场数据地震前兆异常提取方法 |
CN109740453A (zh) * | 2018-12-19 | 2019-05-10 | 吉林大学 | 一种基于小波变换的卫星磁场数据地震前兆异常提取方法 |
CN110688981A (zh) * | 2019-10-14 | 2020-01-14 | 江西理工大学 | 一种振动信号去噪的模态混叠消除方法 |
CN110688981B (zh) * | 2019-10-14 | 2023-06-09 | 江西理工大学 | 一种振动信号去噪的模态混叠消除方法 |
CN112711070A (zh) * | 2019-10-24 | 2021-04-27 | 中国石油化工股份有限公司 | 一种基于地震信号分解的油气检测方法及装置 |
CN112711070B (zh) * | 2019-10-24 | 2024-02-20 | 中国石油化工股份有限公司 | 一种基于地震信号分解的油气检测方法及装置 |
CN111859274A (zh) * | 2020-07-03 | 2020-10-30 | 自然资源部第一海洋研究所 | 一种用于度量复杂动力系统本征因果关系的方法 |
CN112464855A (zh) * | 2020-12-09 | 2021-03-09 | 中国石油天然气集团有限公司 | 基于eemd的随钻泥浆正脉冲信号处理方法及装置 |
Also Published As
Publication number | Publication date |
---|---|
CN107422381B (zh) | 2019-07-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107422381B (zh) | 一种基于eemd-ica的地震低频信息流体预测方法 | |
Herrera et al. | Applications of the synchrosqueezing transform in seismic time-frequency analysis | |
Lu et al. | Seismic spectral decomposition using deconvolutive short-time Fourier transform spectrogram | |
Pinnegar et al. | The S-transform with windows of arbitrary and varying shape | |
Jian et al. | On the denoising method of prestack seismic data in wavelet domain | |
Smith et al. | Extending seismic bandwidth using the continuous wavelet transform | |
CN102798891B (zh) | 基于短时分数阶傅里叶变换的地震信号时频分解方法 | |
Pan et al. | Identification of stratigraphic formation interfaces using wavelet and Fourier transforms | |
CN103995289A (zh) | 基于时频谱模拟的时变混合相位地震子波提取方法 | |
Lapins et al. | An examination of the continuous wavelet transform for volcano-seismic spectral analysis | |
CN102692647B (zh) | 一种高时间分辨率的地层含油气性预测方法 | |
CN103728663B (zh) | 一种时频分析方法 | |
CN106707334B (zh) | 一种提高地震资料分辨率的方法 | |
CN101923176B (zh) | 一种利用地震数据瞬时频率属性进行油气检测的方法 | |
Zoukaneri et al. | A combined Wigner-Ville and maximum entropy method for high-resolution time-frequency analysis of seismic data | |
CN103364832A (zh) | 一种基于自适应最优核时频分布的地震衰减定性估计方法 | |
CN104360382A (zh) | 一种利用叠后地震数据进行油气检测的方法 | |
Liner | An overview of wavelet transform concepts and applications | |
CN106054250A (zh) | 基于变频分量扩散滤波融合的地震资料噪声消减方法 | |
Li et al. | Waveform characteristics of earthquakes induced by hydraulic fracturing and mining activities: Comparison with those of natural earthquakes | |
CN109557578A (zh) | 一种储层含气性检测方法及装置 | |
Wang et al. | Seismic thin interbeds analysis based on high-order synchrosqueezing transform | |
CN107918146A (zh) | 一种基于非线性挤压s时频变换的弱信号检测方法 | |
CN104422956A (zh) | 一种基于稀疏脉冲反演的高精度地震谱分解方法 | |
Biao et al. | Rotating machine fault diagnosis by a novel fast sparsity-enabled feature-energy-ratio method |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |