CN102768366A - 基于高阶统计量的线性与非线性融合的地震子波提取方法 - Google Patents
基于高阶统计量的线性与非线性融合的地震子波提取方法 Download PDFInfo
- Publication number
- CN102768366A CN102768366A CN2011101148500A CN201110114850A CN102768366A CN 102768366 A CN102768366 A CN 102768366A CN 2011101148500 A CN2011101148500 A CN 2011101148500A CN 201110114850 A CN201110114850 A CN 201110114850A CN 102768366 A CN102768366 A CN 102768366A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- model
- theta
- cumulant
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 89
- 239000011159 matrix material Substances 0.000 claims abstract description 40
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 29
- 238000005457 optimization Methods 0.000 claims abstract description 19
- 230000001364 causal effect Effects 0.000 claims description 22
- 238000000605 extraction Methods 0.000 claims description 18
- 241001123248 Arma Species 0.000 claims description 17
- 239000013598 vector Substances 0.000 claims description 15
- 230000002068 genetic effect Effects 0.000 claims description 9
- 238000002474 experimental method Methods 0.000 claims description 7
- 238000004458 analytical method Methods 0.000 claims description 6
- 238000001914 filtration Methods 0.000 claims description 6
- 230000004927 fusion Effects 0.000 claims description 4
- 230000004044 response Effects 0.000 claims description 4
- 230000003247 decreasing effect Effects 0.000 claims description 3
- 230000002194 synthesizing effect Effects 0.000 claims description 3
- 238000011156 evaluation Methods 0.000 claims description 2
- 238000012545 processing Methods 0.000 abstract description 9
- 230000036039 immunity Effects 0.000 abstract 1
- 230000001373 regressive effect Effects 0.000 abstract 1
- 230000006870 function Effects 0.000 description 28
- 238000004088 simulation Methods 0.000 description 6
- 238000001228 spectrum Methods 0.000 description 6
- 230000008569 process Effects 0.000 description 5
- 230000000694 effects Effects 0.000 description 4
- 241000282326 Felis catus Species 0.000 description 3
- 239000000654 additive Substances 0.000 description 3
- 230000000996 additive effect Effects 0.000 description 3
- 230000000739 chaotic effect Effects 0.000 description 3
- 230000007547 defect Effects 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 238000007500 overflow downdraw method Methods 0.000 description 3
- 238000009825 accumulation Methods 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 238000013507 mapping Methods 0.000 description 2
- 238000013459 approach Methods 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 230000005284 excitation Effects 0.000 description 1
- 230000008450 motivation Effects 0.000 description 1
- 230000008439 repair process Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
Images
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明公开一种基于高阶统计量的线性与非线性融合的地震子波提取方法,其特点在于:对ARMA模型子波建立拟合目标函数后,利用线性与非线性相结合的方法求解地震子波参数。首先采用累积量矩阵方程法对子波模型进行初步估计,然后将得到的子波预估计值用于拟合优化算法初始参数搜索空间的确定,并在此基础上用累积量拟合误差调整矩阵方程法的阈值,寻找精确寻优区间,进而由累积量拟合方法得到模型准确的阶数及参数值。本发明提出的子波提取方法抗噪性能强,在数据较短时也能提取出较高精度的地震子波,在实际地震数据处理中有很好的应用价值。
Description
技术领域:
本发明属于地震信号处理领域。
背景技术:
地震子波估计是地震信号处理中的关键问题之一,是地震正演、地震反演和反褶积等地震处理技术的基础。高精度的地震子波估计具有重要的实用价值,是油气地震勘探领域高分辨率处理技术亟待解决的重要研究课题。
在信号处理和系统理论等领域使用高阶统计量的主要动机与出发点可以归结为:(1)抑制加性有色噪声(其功率谱未知)的影响;(2)辨识非因果、非最小相位系统或重构非最小相位信号;(3)抽取由于高斯性偏离引起的各种信息;(4)检测和表征信号中的非线性以及辨识非线性系统;(5)检测和表征信号中的循环平稳性以及分析和处理循环平稳信号。高阶统计量能够保留系统的相位信息、并且理论上可以完全抑制高斯噪声,在解决非线性、非高斯和非最小相位系统的分析和处理方面已逐渐显示出独特优势,为高精度的地震子波提取提供了可能。由高阶累积量缺陷的分析可知,对于高阶统计的运算来说,数据量越大越有利于提高估计质量,该累积量的缺陷本可用于增加地震记录数据(如以多道联合求解)的方式进行缩减或消除。然而不实用的地方在于:一般来说,对于实际的地震记录,地震子波是时空变的,即使同一道记录,浅层和深层的地震子波也是不相同的,这主要是因为大地的滤波作用,使地震子波频率随着地层深度增加而降低。为了能够准确估计出地震子波,在地震子波为时不变模型的假设下,要求参与运算的地震记录长度不能太长,甚至应将浅层和深层数据分开来计算,导致了所求得的地震数据累积量有很大的估计方差。
累积量矩阵方程法及累积量拟合法均可从理论上提取出非高斯随机过程的模型参数。然而,这两种方法又各有自己的缺陷:相对于非线性参数估计方法(累积量拟合法),线性化参数估计方法(累积量矩阵方程法)运算速度较快,但由于仅利用了地震记录累积量的特殊切片信息,在对高噪信比、短数据地震记录进行子波提取时,该方法可能存在较大的估计误差;累积量拟合优化方法充分利用了地震记录的累积量信息,子波提取具有更高的精确度和数值稳定性,但用累积量拟合法确定模型初始参数范围较为困难,严重影响了算法的寻优效率。虽然比较稳定,但该方法对模型参数的求取异常困难。这是由于累积量拟合目标函数多为一多维多峰值的非线性函数,通常的优化算法对该目标函数的求解无能为力。在对目标函数的仿真中会发现,该目标函数的拟合误差对模型参数的变化极为敏感,使得线性化优化算法的参数估计结果严重依赖于初始估值,尤其在ARMA描述下,这一特点将更为突出。
针对累计量矩阵方程法和累积量拟合法的缺陷,将两者结合,提出一种基于高阶累积量的线性与非线性融合的地震子波提取方法,首先利用矩阵方程法确定子波模型的阶数及参数范围,然后根据累积量拟合公式精确估计模型参数,以弥补两种方法单独使用时的缺陷。
考虑利用累积量拟合误差反馈矩阵方程法所得初始估计,以拟合误差的大小对累积量矩阵方程法的估计参数及阈值进行调整,并根据调整过的参数重新确定参数寻优空间,进而求得模型的全局最优解。
发明内容:
本发明的目的在于提出一种基于高阶累积量的线性与非线性方法融合的高精度子波提取方法。
本发明的特征在于,对ARMA模型子波建立拟合目标函数后,利用线性与非线性相结合的方法求解子波参数。首先采用矩阵方程法对子波模型进行初步估计,然后将得到的子波预估计值用于拟合优化算法初始参数搜索空间的确定,并在此基础上用累积量拟合误差调整矩阵方程法的阈值,寻找精确寻优区间,进而由累积量拟合方法得到模型准确的阶数及参数值。该方法依次含有以下步骤:
步骤1数据生成:用待估计ARMA子波模型与满足反射系数序列假设的随机序列合成地震数据记录y(n);
步骤2构建拟合目标函数E(θ):基于地震褶积模型的假设,构建多峰多极值目标函数E(θ): 其中 C4y(τ1,τ2,τ3|θ)分别为地震记录四阶累积量和估计模型的四阶矩归一化后的值,令待估计子波的参数向量为θ′=(a1,…,ap,b1,…,bq),子波模型的四阶累积量可表示为: 其中,h(i,θ)为模型参数为θ的ARMA模型的冲激响应;n≥max(lc,lnc),lc,lnc分别为参数向量θ′对应子波模型因果、非因果部分的长度,时延τi的取值区间仅需包含:{(0≤τ1≤q+3p),max(0,q-p)≤τ2≤min(q+3p,τ1),0≤τ3≤min(2p,τ2)}。
步骤3阶数预估计:用累积量矩阵方程法初步确定待估计模型的阶数p,q;
步骤4确定参数搜索空间:在给定的模型阶数值p,q下,用矩阵方程法初步估计模型参数,并以此确定模型参数向量θ的搜索范围;
步骤5拟合优化:在给定参数向量的搜索范围内用优化算法对累积量拟合目标函数E(θ)进行参数精确估计,寻找最优解θe;
步骤6评价函数分析:将拟合误差作为目标适应度函数,并计算最优解θe的目标适应度函数值,当其有明显增加或降低时,转至步骤8,否则,转至步骤7;
步骤7阶数调整:根据累积量拟合误差对模型的阶数进行调整,生成一组新的p,q值,并根据拟合误差和生成的新阶数值调整累积量矩阵方程法中PODE或SVD的参数阈值,然后转至步骤4;
步骤8结束:将当次阶数调整前所得阶数p,q及对应精确参数估计值θe视为模型最优解。
本发明的特征之二在于基于地震褶积模型假设构建多峰多极值目标函数E(θ),采用累积量矩阵方程法对子波模型进行初步估计,即利用观测信号累积量组成Hankel矩阵的扩展矩阵来估计AR反因果部分的秩和参数,用估计出的AR反因果部分参数滤波原观测信号,用基于样本高阶累积量估计因果系统AR参数的SVD-TLS算法来估计模型AR因果部分参数,利用AR模型预测残差矩阵有效秩和PODE(对角元素乘积)实验相结合确定MA阶次q,然后将得到的子波预估计值作为改进的遗传优化算法初始值并用于初始参数搜索空间的确定,在此基础上用累积量拟合误差即E(θ)调整矩阵方程法的阈值,寻找精确寻优区间,进而由累积量拟合方法得到最佳模型参数向量θ,准确的阶数p,q和模型参数a1,…,ap,b1,…bq。 实验仿真说明:
本发明通过合成地震记录的子波估计实验来检验线性与非线性结合方法提取地震子波算法的效果。首先合成一长度为800(采样间隔为1ms)的短数据地震记录。反射系数序列服从伯努利分布,子波为非因果、非最小相位,其参数如表1所示,长度取100ms,并分别加入5%(噪声与信号的能量比)的高斯有色噪声进行实验。分别采用累积量矩阵方程法和线性与非线性相结合方法提取地震子波,所得子波参数如表1所示,两种方法所提取出的子波波形如图2所示。
从表1和图2中可以看出,使用线性与非线性融合方法提地震子波可以在短数据情况下有效地估计出非因果、非最小相位子波。在这个例子中,由于合成地震记录的长度仅为800ms,若仅使用基于累积量的矩阵方程法所得的子波于原始波形相差较大,而采用基于线性和非线性的混合算法的估计误差并不十分明显,相对具有更高的估计精度。
为了分析在不同数据长度下,线性与非线性方法提取的子波的稳定性,进行如下实验仿真:在25次试验中每次均随机生成长度分别为20000ms,5000ms,2000ms,500ms的合成 地震记录,采用线性与非线性融合方法提取子波,结果如图3所示。由图3可以看出,对于不同数据长度的地震记录,提取的地震子波与真实子波十分接近,因此该方法可提取高精度的地震子波,特别是在地震数据短的情况。
为了分析在不同噪声强度和成份下,线性与非线性融合方法提取的子波的稳定性,进行如下实验仿真:为考察加性噪声对本发明方法的影响,每次实验均任意生成一段长度为20000ms的合成记录,对该记录分别加入不同强度的高斯、非高斯色噪声后应用本发明方法进行子波提取。
由图4(a)看出,本发明方法在高斯色噪声成份较高的情况下仍能较为准确的辨识出实际的波形。由图4(b)展示了非高斯色噪声的加入后对子波提取的效果影响并不明显。
为验证本发明方法的实际应用价值,对一实际地震剖面进行地震子波提取试验。如图5(a)所示,该地震剖面数据为201道,采样率为1ms,数据记录长度为2s。地震子波提取方式为每隔20道均随机选取该地震剖面中的一道应用本发明方法进行模型参数估计,其结果如图5(b)所示。应用本发明方法所提取出的10道地震子波,其波形在相位上保持了一定的连续性,从另一个方面表明了本发明方法的一致性。在实际地震资料的子波提取过程中发现,由于随机性地选择20道中的一道地震记录,对剖面中某些道数据进行子波提取时,子波提取结果并非完全一致,由图5(b)也可看出,各地震波波形均有一定的偏差,分析认为是实际地震资料中微曲多次反射以及反射系数序列分布偏离白噪假设时对地震子波波形的滤波效应所致。
综上所述,本发明提出的方法可以有效的提高地震子波的估计精度,抗噪性能强,特别是在地震数据短的情况下,也能得到精度比较高的地震子波。
附图说明:
图1、本发明所属方法的计算机流程图
图2、使用矩阵方程法和线性与非线性融合方法提取的子波比较
图3、不同数据长度下提取的子波模型参数
图4、不同强度噪声成份下提取的子波波形
图5、实际地震剖面(a)和真实地震资料中提取的子波波形(b)
具体实施方式:
本发明提出了基于高阶累积量的线性与非线性优化算法相结合的ARMA模型参数估计方法。这一算法充分利用了线性和非线性模型参数辨识方法的优点,且具有以下特点:
1.相对于MA模型,基于高阶累积量的ARMA模型地震子波提取方法对地震子波的描述还具有参数更为节俭的特点,理论仿真分析和实际数据处理结果表明了该方法的可行性 和高效性。
2.运用累积量矩阵方程法求解ARMA模型初始值和初始搜索范围,提高了子波提取的运算效率。
3.运用累积量拟合法构建拟合目标函数,进而求得基于高阶累积量的ARMA模型描述下的地震子波提取方法,充分地利用了地震记录的累积量信息,使方法有更好的稳定性,参数估计结果受数据量大小以及噪声强度的影响较小。
4.利用累积量拟合误差反馈矩阵方程法所得初始估计,以拟合误差的大小对初始估计进行参数调整,并根据调整过的参数重新确定参数寻优空间,进而求得模型的全局最优解,可稳定而准确地从短数据高噪声污染的地震记录中提取出混合相位的地震子波。
5.本发明所提出的技术也可以应用于系统辨识领域。
本发明按以下步骤实施:
1.用待估计ARMA子波模型与满足反射系数序列假设的随机序列合成地震数据记录y(n)。
2.基于地震褶积模型假设,构建多峰多极值目标函数E(θ)。
3.用累积量矩阵方程法初步确定待估计模型的阶数p,q,确定模型参数向量θ的搜索范围。
4.将累积量矩阵方程法得到的模型阶数和参数作为预估计值,在给定参数向量的搜索范围内用优化算法对累积量拟合目标函数E(θ)进行参数精确估计,寻找最优解θe。
5.将拟合误差作为目标适应度函数,当其有明显增加或降低时,转至6,否则,根据累积量拟合误差对模型的阶数进行调整,生成一组新的p,q值,转至4。
6.将当次阶数调整前所得阶数p,q及对应精确参数估计值θe视为模型最优解,输出地震子波。 具体原理如下:
本发明利用累积量拟合误差反馈矩阵方程法所得初始估计,以拟合误差的大小对累积量矩阵方程法的估计参数及阈值进行调整,并根据调整过的参数重新确定参数寻优空间,进而求得模型的全局最优解。
1.应用累积量矩阵方程法对模型参数进行初步估计
通常地震记录可以假设为一零均值平稳随机过程y(n),构建一随机ARMA过程如下: y(n)=x(n)+v(n) (1) 其中x(n)为满足下列差分方程的非高斯信号,是ARMA(p,q)模型在输入为反射系数序列r(n) 时的响应: v(n)是满足下式的有色高斯观测噪声:
对上述信号作如下假设:
假设2:v(n)为环境噪声、仪器噪声及激发产生的多次反射噪声等的合成噪声信号,与r(n)统计独立,一般为加性色噪声,且其高斯成分远大于非高斯部分。
假设3:与式(2)对应的ARMA(p,q)模型为非因果、非最小相位系统,其中a(0)=1。其传递函数为: 其中x(n)为稳定信号,所以当|z|=1时,A(z)≠0,并且H(z)中无零极点对消。
假设4:H(z)有p1个圆内极点,有p2=p-p1个圆外极点。令αj为H(z)的极点,则有 A(z)=A1(z)A2(z)g(z) (5) 式中 称之为H(z)的因果AR部分;而 称 为H(z)的反因果AR部分,显然 (常数);
假设5:H(z)有q1个圆内零点,有q2=q-q1个圆外零点。则与极点情况类似有 B(z)=B1(z)B2(z)s(z)其中B1(z)为最小相位部分,B2(z)s(z)为非最小相位部分。
应用累积量矩阵方程法求解子波初始解,步骤如下:
(1)求功率谱等价意义下的因果最小相位模型的AR参数
由式(1)~(3)可以得到:
其中Y(z)、R(z)、C(z)、D(z)、E(z)分别是y(n)、r(n)、c(i)、d(j)、e(n)的Z变换,根据式(5)和式(7)易知观测信号y(n)的功率谱可以表示为:
由于A2(z)的根均在单位圆外,将A2(z)所对应的ARMA模型的极点取共轭倒数,而反演到单位圆内得到 则 的根均在单位圆内。因 与A2(z)的功率谱等价,现构造 则A1(z)A2(z)与 功率谱等价。由功率谱等价原理,并利用Cadzow的自相关法估计AR参数的SVD-TLS(奇异值分解-总体最小二乘)算法可以估计出 这种基于观测信号自相关来估计因果信号AR参数的算法已经比较成熟,且仿真结果证明它有很高的分辨率和精度。
(2)求H(z)的反因果部分AR参数
(3)求H(z)的因果部分AR参数及MA参数
利用估计出的A2(z)g(z)滤波原观测信号y(n)得到y2(n),则其Z变换为:
由上式可见y2(n)是一因果信号和一高斯有色噪声的叠加。利用基于样本高阶累积量估计因果系统AR参数的SVD-TLS算法可估计出因果部分的AR参数a(1),a(2),...,a(p1)。然后利用AR模型预测残差矩阵有效秩和PODE(对角元素乘积)实验相结合可确定MA阶次q。用求得的因果部分AR参数a(0),a(1),...,a(p1)对y2(n)进行滤波可得到一纯MA过程y3(n),其 Z变换为: 最后采用只用高阶累积量确定MA参数b(0),b(1),...,b(q)的方法即可得到MA参数。
2.累积量拟合法优化参数
基于地震褶积模型的假设,构建目标函数: 其中, C4y(τ1,τ2,τ3|θ)分别为地震记录四阶累积量和估计模型的四阶矩归一化后的值。修改参数向量θ,令待估计子波的参数向量为θ′=(a1,...,ap,b1,...,bq),子波模型的四阶累积量C4y(τ1,τ2,τ3|θ)可表示为: 其中,h(i,θ)为模型参数为θ,且描述符合式(3)的ARMA模型的冲激响应;n≥max(lc,lnc),lc,lnc分别为参数向量θ′对应子波模型因果、非因果部分的长度。由累积量的对称性,时延τi的取值区间仅需包含:{(0≤τ1≤q+3p),max(0,q-p)≤τ2≤min(q+3p,τ1),0≤τ3≤min(2p,τ2)}即可。
目标函数E(θ)多峰多极值函数,采用改进的遗传算法求解参数精确解。不仅各参数间的依赖性较大,且函数值模型对某些特定参数尤其是AR部分参数的变化较为敏感。因此,求解该目标函数所用寻优算法须兼有对参数向量整体的全局随机搜索能力和对单个参数的深度搜索能力,遗传算法在搜索空间的深度搜索和广度搜索之间维持了很好的平衡性,结合本课题目标函数的特点,本发明采用改进的遗传算法对其进行寻优,将遍历均匀性好、不容易陷入小循环或者不动点的猫映射(cat map)引入到混沌遗传算法中并进行了改进,提出了较遗传算法性能更为稳定有效的、基于猫映射的混沌遗传算法。计算最优解θe的目标函数适应度即拟合误差E(θ)(确定在3.5左右),当其有明显增加或降低时,继续直到接近目标函数适应度,否则,结束。根据累积量拟合误差对模型的阶数进行调整,生成一组新的p,q值,并根据拟合误差和生成的新阶数值调整累积量矩阵方程法中PODE或SVD的参数阈值。
Claims (2)
1.基于高阶累积量的线性与非线性融合的地震子波提取方法的思想,其特征在于,该方依次含有以下步骤:
步骤(1)初始数据生成y(n):用待估计ARMA子波模型与满足反射系数序列假设的随机序列合成地震记录y(n);
步骤(2)构建拟合目标函数E(θ):基于地震褶积模型的假设,构建多峰多极值目标函数:
其中,h(i,θ)为模型参数为θ的ARMA模型的冲激响应;n≥max(lc,lnc),lc,lnc分别为参数向量θ′对应子波模型因果、非因果部分的长度,时延τi的取值区间仅需包含:{(0≤τ1≤q+3p),max(0,q-p)≤τ2≤min(q+3p,τ1),0≤τ3≤min(2p,τ2)};
步骤(3)阶数预估计:用累积量矩阵方程法初步确定待估计模型的阶数p,q;
步骤(4)确定参数搜索空间:对给定阶数p,q值,用矩阵方程法对该阶数进行参数预估计,并由此确定解空间(a,b),确定模型参数向量θ的搜索范围;
步骤(5)拟合优化:在给定解空间内用拟合优化算法对累积量拟合目标函数E(θ)进行参数精确估计,寻找最优解θe;
步骤(6)评价函数分析:将拟合误差作为目标适应度函数,并计算最优解θe的目标适应度函数值,当其有明显增加或降低时,转至(8),否则,转至步骤(7);
步骤(7)阶数调整:根据累积量拟合误差对模型的阶数进行调整,生成一组新的p,q值,并根据拟合误差和生成的新阶数值调整累积量矩阵方程法中PODE或SVD的参数阈值,然后转至步骤(4);
步骤(8)结束:将当次扰动前所得阶数p,q及参数估计值θe视为模型最优解。
2.基于高阶累积量的线性与非线性相结合的地震子波提取方法,其特征在于基于地震褶积模型假设构建多峰多极值目标函数E(θ),采用累积量矩阵方程法(线性化方法)对子波模型进行初步估计,即利用观测信号累积量组成Hankel矩阵的扩展矩阵来估计AR反因果部分的秩和参数,用估计出的AR反因果部分参数滤波原观测信号,用基于样本高阶累积量估计因果系统AR参数的SVD-TLS算法来估计模型AR因果部分参数,利用AR模型预测残差矩阵有效秩和PODE(对角元素乘积)实验相结合确定MA阶次q,然后将得到的子波预估计值作为改进的遗传优化算法初始值并用于初始参数搜索空间的确定,在此基础上用累积量拟合误差即E(θ)作为适应度函数调整矩阵方程法的阈值,寻找精确寻优区间,进而由累积量拟合方法得到最佳模型参数向量θe、准确的阶数p,q和模型参数(a1,…,ap,b1, …,bq)。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2011101148500A CN102768366A (zh) | 2011-05-04 | 2011-05-04 | 基于高阶统计量的线性与非线性融合的地震子波提取方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2011101148500A CN102768366A (zh) | 2011-05-04 | 2011-05-04 | 基于高阶统计量的线性与非线性融合的地震子波提取方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN102768366A true CN102768366A (zh) | 2012-11-07 |
Family
ID=47095829
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2011101148500A Pending CN102768366A (zh) | 2011-05-04 | 2011-05-04 | 基于高阶统计量的线性与非线性融合的地震子波提取方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102768366A (zh) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104516771A (zh) * | 2015-01-22 | 2015-04-15 | 黄国庆 | 一种非平稳随机过程高效模拟方法 |
CN104849758A (zh) * | 2015-05-05 | 2015-08-19 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 针对地震数据中的规则干扰的压制方法 |
CN108646291A (zh) * | 2018-05-04 | 2018-10-12 | 北京信息科技大学 | 基于果蝇神经网络算法的子波整形反褶积处理方法及装置 |
CN108983292A (zh) * | 2018-04-28 | 2018-12-11 | 西安交通大学 | 基于三角滤波的局部倾斜叠加波束形成偏移假频压制方法 |
CN110727030A (zh) * | 2019-10-22 | 2020-01-24 | 中国石油化工股份有限公司 | 一种基于子波多极值特征联合分析的储层预测方法 |
CN111427080A (zh) * | 2020-03-13 | 2020-07-17 | 王仰华 | 地震数据空变广义子波的提取方法 |
CN114002740A (zh) * | 2020-07-28 | 2022-02-01 | 北京大学 | 一种低信噪比地震初至信号自动增强方法及系统 |
-
2011
- 2011-05-04 CN CN2011101148500A patent/CN102768366A/zh active Pending
Non-Patent Citations (1)
Title |
---|
戴永寿等: "基于高阶累积量ARMA模型线性非线性结合的地震子波提取方法研究", 《地球物理学报》 * |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104516771A (zh) * | 2015-01-22 | 2015-04-15 | 黄国庆 | 一种非平稳随机过程高效模拟方法 |
CN104849758A (zh) * | 2015-05-05 | 2015-08-19 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 针对地震数据中的规则干扰的压制方法 |
CN104849758B (zh) * | 2015-05-05 | 2017-10-03 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 针对地震数据中的规则干扰的压制方法 |
CN108983292A (zh) * | 2018-04-28 | 2018-12-11 | 西安交通大学 | 基于三角滤波的局部倾斜叠加波束形成偏移假频压制方法 |
CN108983292B (zh) * | 2018-04-28 | 2019-10-11 | 西安交通大学 | 基于三角滤波的局部倾斜叠加波束形成偏移假频压制方法 |
CN108646291A (zh) * | 2018-05-04 | 2018-10-12 | 北京信息科技大学 | 基于果蝇神经网络算法的子波整形反褶积处理方法及装置 |
CN110727030A (zh) * | 2019-10-22 | 2020-01-24 | 中国石油化工股份有限公司 | 一种基于子波多极值特征联合分析的储层预测方法 |
CN111427080A (zh) * | 2020-03-13 | 2020-07-17 | 王仰华 | 地震数据空变广义子波的提取方法 |
CN111427080B (zh) * | 2020-03-13 | 2021-08-13 | 王仰华 | 地震数据空变广义子波的提取方法 |
CN114002740A (zh) * | 2020-07-28 | 2022-02-01 | 北京大学 | 一种低信噪比地震初至信号自动增强方法及系统 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102768366A (zh) | 基于高阶统计量的线性与非线性融合的地震子波提取方法 | |
CN111596366B (zh) | 一种基于地震信号优化处理的波阻抗反演方法 | |
CN103389094B (zh) | 一种改进的粒子滤波方法 | |
CN113240170A (zh) | 一种基于季节性循环神经网络的空气质量预测方法 | |
CN114662414B (zh) | 一种基于图小波神经网络模型的油藏生产预测方法 | |
CN103645500B (zh) | 一种频率域的混合相位地震子波估算方法 | |
CN116029459B (zh) | 结合图傅里叶变换tmgcn交通流预测模型的提取方法 | |
CN102184408A (zh) | 基于自回归模型的雷达高分辨距离像目标识别方法 | |
CN105467446A (zh) | 基于径向高斯核的自适应最优核时频分析方法 | |
CN102096101B (zh) | 提高地震数据处理精度的方法及装置 | |
CN113486574A (zh) | 基于历史数据以及机器学习的声速剖面补全方法及装置 | |
CN115240418B (zh) | 基于因果门控-低通图卷积网络的短时交通流量预测方法 | |
CN111159888A (zh) | 一种基于互相关函数的协方差矩阵稀疏迭代时延估计方法 | |
CN103500455A (zh) | 一种基于无偏有限冲击响应滤波器(ufir)的改进机动目标跟踪方法 | |
CN102768365A (zh) | 基于高阶统计量和arma模型的高分辨率地震子波提取方法 | |
CN105043390A (zh) | 基于泛克里金法的重力场插值方法 | |
CN102928875A (zh) | 基于分数阶傅里叶域的子波提取方法 | |
CN111142134B (zh) | 一种坐标时间序列处理方法及装置 | |
CN109582915B (zh) | 应用于纯方位跟踪的改进非线性可观测度自适应滤波方法 | |
CN108957529B (zh) | 一种基于属性的无井区子波估计方法 | |
CN112578438A (zh) | 一种地震子波提取方法及系统 | |
CN108646291B (zh) | 基于果蝇神经网络算法的子波整形反褶积处理方法及装置 | |
CN104469795A (zh) | 基于领域搜索的FARIMA模型Hurst参数估计方法 | |
CN107346300B (zh) | 一种基于绝对传递率函数的传递路径分析方法 | |
CN114428291A (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 | ||
C02 | Deemed withdrawal of patent application after publication (patent law 2001) | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20121107 |