CN105182417B - 一种基于形态成分分析的面波分离方法及系统 - Google Patents
一种基于形态成分分析的面波分离方法及系统 Download PDFInfo
- Publication number
- CN105182417B CN105182417B CN201510578757.3A CN201510578757A CN105182417B CN 105182417 B CN105182417 B CN 105182417B CN 201510578757 A CN201510578757 A CN 201510578757A CN 105182417 B CN105182417 B CN 105182417B
- Authority
- CN
- China
- Prior art keywords
- wave signal
- surface wave
- signal component
- dictionary
- dimensional
- 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.)
- Expired - Fee Related
Links
- 238000004458 analytical method Methods 0.000 title claims abstract description 53
- 238000000926 separation method Methods 0.000 title claims abstract description 43
- 230000000877 morphologic effect Effects 0.000 claims description 48
- 238000005070 sampling Methods 0.000 claims description 29
- 238000012545 processing Methods 0.000 claims description 25
- 238000000034 method Methods 0.000 claims description 24
- NRNCYVBFPDDJNE-UHFFFAOYSA-N pemoline Chemical compound O1C(N)=NC(=O)C1C1=CC=CC=C1 NRNCYVBFPDDJNE-UHFFFAOYSA-N 0.000 claims description 12
- 230000014759 maintenance of location Effects 0.000 claims description 6
- 230000000694 effects Effects 0.000 abstract description 12
- 230000010365 information processing Effects 0.000 abstract description 2
- 238000010586 diagram Methods 0.000 description 9
- 238000001914 filtration Methods 0.000 description 4
- 238000013519 translation Methods 0.000 description 4
- 239000006185 dispersion Substances 0.000 description 3
- 238000002474 experimental method Methods 0.000 description 3
- 238000012986 modification Methods 0.000 description 3
- 230000004048 modification Effects 0.000 description 3
- 230000009466 transformation Effects 0.000 description 3
- 238000005056 compaction Methods 0.000 description 2
- 230000007547 defect Effects 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 230000002194 synthesizing effect Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 239000002131 composite material Substances 0.000 description 1
- 229910052500 inorganic mineral Inorganic materials 0.000 description 1
- 239000011707 mineral Substances 0.000 description 1
- 238000003825 pressing Methods 0.000 description 1
- 230000001629 suppression Effects 0.000 description 1
- 238000011426 transformation method Methods 0.000 description 1
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明属于地震波信息处理技术领域，尤其涉及一种基于形态成分分析的面波分离方法及系统。所述基于形态成分分析的面波分离方法包括：步骤a：制作合成地震记录，并构造基于二维字典形态成分分析的面波分离模型；步骤b：选取二维非抽样小波字典稀疏表示面波信号分量，并选取二维局部离散余弦字典稀疏表示反射波信号分量；步骤c：采用二维字典下形态成分分析对面波信号分量与反射波信号分量进行分离。本发明既考虑面波信号特性的同时也考虑到反射波信号的特性，能更加有效地分离出面波与反射波信号。本发明充分利用了道与道之间的相关性，同时考虑在二维非抽样小波字典下，面波信号集中在低频高波数域中，使得面波信号的分离效果更好。
Description
技术领域
本发明属于地震波信息处理技术领域，尤其涉及一种基于形态成分分析的面波分离方法及系统。
背景技术
地震勘探是指通过人工地震方法探测地质结构以及地下矿藏资源。目前，该方法已经成为地球物理勘探最广泛采用的一种方法。通常野外采集到的地震资料中包含大量噪声干扰，因此需要对野外采集的数据进行处理，噪声压制成为地震资料处理中必不可少的环节。面波是反射波地震勘探中主要的规则干扰之一，是对地震资料信噪比影响较大的干扰波，因此对面波的压制有利于提高地震资料的信噪比，便于后续的地震资料处理以及解释工作。
面波主要分布在炮点附近，呈扇形形状分布，具有低频率、强能量及频散等特点。现有技术中，常用滤除面波的方法包括一维傅里叶变换滤波法、f-k变换滤波法、τ-p变换法、小波变换法及曲波变换法等。现有的滤除面波方法通常仅采用单一变换的方法滤除面波信号，单一变换只能利用面波的低频、强能量及低视速度这些特点中的部分特性对面波信号进行压制，不能完全地消除面波信号。
近年来，稀疏表示理论成为地震资料处理领域内较为热门的方法，基于稀疏表示理论的形态成分分析方法可以根据信号不同成分的形态结构特征差异选择合适的字典稀疏表示信号的不同成分，目前该方法已经成功用于单道地震记录的面波与反射波的分离处理中。但是，现有的基于形态成分分析的面波分离方法，采用的地震信号表示字典为一维字典，逐道对地震数据进行处理，不能有效利用地震资料中道与道的相关性；同时，采用的一维处理算法在处理二维地震数据时处理速度较慢。
发明内容
本发明提供了一种基于形态成分分析的面波分离方法及系统，旨在解决现有的基于形态成分分析的面波分离方法采用的地震信号表示字典为一维字典，逐道对地震数据进行处理，不能有效利用地震资料中道与道的相关性，且在处理二维地震数据时处理速度较慢的技术问题。
本发明是这样实现的，一种基于形态成分分析的面波分离方法，包括：
步骤a：制作合成地震记录，并构造基于二维字典形态成分分析的面波分离模型；
步骤b：选取二维非抽样小波字典稀疏表示面波信号分量，并选取二维局部离散余弦字典稀疏表示反射波信号分量；
步骤c：采用二维字典下形态成分分析对面波信号分量与反射波信号分量进行分离。
本发明实施例采取的技术方案还包括：在所述步骤a中，所述合成地震记录为：
x＝xg+xr+n (1)
式(1)中，xg为面波信号分量，xr为反射波信号分量，n为随机噪声。
本发明实施例采取的技术方案还包括：在所述步骤a中，所述面波分离模型为：
式(2)中，λ为拉格朗日乘子，控制着面波信号分量与反射波信号分量之和拟合地震记录的失真度；Tg、Tr分别稀疏表示面波信号分量与反射波信号分量。
本发明实施例采取的技术方案还包括：在所述步骤c中，所述采用二维字典下形态成分分析对面波信号分量与反射波信号分量进行分离还包括：将地震记录分为面波信号分量与反射波信号分量，将二维非抽样离散小波字典与二维局部离散余弦字典作为地震记录的联合稀疏表示字典，构建稀疏表示模型；所述稀疏表示模型通过对式(2)进行转换得到：
本发明实施例采取的技术方案还包括：在所述步骤c中，所述采用二维字典下形态成分分析对面波信号分量与反射波信号分量进行分离还包括：通过块协调松弛算法对稀疏表示模型进行求解，得到分离出来的面波信号分量以及反射波信号分量。
本发明实施例采取的技术方案还包括：所述通过块协调松弛算法对稀疏表示模型进行求解具体包括：
(1)初始化拉格朗日乘子λ，迭代次数N，阈值δ＝Lmax，截止阈值δmin，其中Lmax为地震记录x在两种字典下表示系数的最大值中较小的一个；
(2)进行N次迭代：
A部分：假设xr不变，更新xg
①计算残差：r＝x-xg-xr；
②计算xg+r经过二维非抽样离散小波变换后的系数，
③保留系数αg低频高波数域中系数不变，对系数αg的其它域小波系数进行软阈值处理，阈值为δ，得到系数
④通过重建xg；
B部分：假设xg不变，更新xr
⑤计算残差r＝x-xg-xr；
⑥计算xr+r经过二维局部离散余弦变换后的系数，
⑦对系数αr进行软阈值处理，阈值为δ，得到系数
⑧通过重建xr；
(3)更新阈值δ＝δ*λ；
(4)如果δ>δmin转至步骤(2)，否则结束。
本发明实施例采取的另一技术方案为：一种基于形态成分分析的面波分离系统，包括地震记录制作模块、二维字典选取模块和面波分离模块；
所述地震记录制作模块用于制作合成地震记录，并构造基于二维字典形态成分分析的面波分离模型；
所述二维字典选取模块用于选取二维非抽样小波字典稀疏表示面波信号分量，并选取二维局部离散余弦字典稀疏表示反射波信号分量；
所述面波分离模块用于采用二维字典下形态成分分析对面波信号分量与反射波信号分量进行分离。
本发明实施例采取的技术方案还包括：所述合成地震记录为：
x＝xg+xr+n (1)
式(1)中，xg为面波信号分量，xr为反射波信号分量，n为随机噪声；
所述面波分离模型为：
式(2)中，λ为拉格朗日乘子，控制着面波信号分量与反射波信号分量之和拟合地震记录的失真度；Tg、Tr分别稀疏表示面波信号分量与反射波信号分量。
本发明实施例采取的技术方案还包括：所述面波分离模块采用二维字典下形态成分分析对面波信号分量与反射波信号分量进行分离还包括：将地震记录分为面波信号分量与反射波信号分量，将二维非抽样离散小波字典与二维局部离散余弦字典作为地震记录的联合稀疏表示字典，构建稀疏表示模型；所述稀疏表示模型通过对式(2)进行转换得到：
本发明实施例采取的技术方案还包括：所述面波分离模块采用二维字典下形态成分分析对面波信号分量与反射波信号分量进行分离还包括：通过块协调松弛算法对稀疏表示模型进行求解，得到分离出来的面波信号分量以及反射波信号分量；所述通过块协调松弛算法对稀疏表示模型进行求解具体包括：
(1)初始化拉格朗日乘子λ，迭代次数N，阈值δ＝Lmax，截止阈值δmin，其中Lmax为地震记录x在两种字典下表示系数的最大值中较小的一个；
(2)进行N次迭代：
A部分：假设xr不变，更新xg
①计算残差：r＝x-xg-xr；
②计算xg+r经过二维非抽样离散小波变换后的系数，
③保留系数αg低频高波数域中系数不变，对系数αg的其它域小波系数进行软阈值处理，阈值为δ，得到系数
④通过重建xg；
B部分：假设xg不变，更新xr
⑤计算残差r＝x-xg-xr；
⑥计算xr+r经过二维局部离散余弦变换后的系数，
⑦对系数αr进行软阈值处理，阈值为δ，得到系数
⑧通过重建xr；
(3)更新阈值δ＝δ*λ；
(4)如果δ>δmin转至步骤(2)，否则结束。
本发明实施例的基于形态成分分析的面波分离方法及系统根据地震记录中面波信号与反射波信号的结构特征差异，分别选取二维非抽样小波字典稀疏表示面波信号，选取二维局部离散余弦字典稀疏表示反射波信号，相对于单一变换，在分离面波与反射波的过程中，既考虑面波信号特性的同时也考虑到反射波信号的特性，从而能更加有效地分离出面波与反射波信号。同时，本发明利用二维字典对二维地震数据进行整体稀疏表示与处理，充分利用了道与道之间的相关性，得到较好的分离效果。并且面波信号与反射波信号在二维字典下具有更好的辨识度，同时考虑在二维非抽样小波字典下，面波信号集中在低频高波数域中，使得面波信号的分离效果更好。
附图说明
图1是本发明实施例的基于形态成分分析的面波分离方法的流程图；
图2为实际地震记录切片效果图；
图3为图2所示的地震记录在二维非抽样离散小波字典与二维局部离散余弦字典表示中，稀疏性与重构相对误差的关系比较图；
图4为含有强能量面波干扰的实际地震记录图；
图5为图4所示的地震记录在二维非抽样离散小波字典与二维局部离散余弦字典表示中，稀疏性与重构相对误差的关系比较图；
图6为合成的反射波地震记录；
图7为合成的含有面波与反射波的地震记录；
图8为对图7所示的合成地震记录利用二维字典的形态成分分析压制面波后的结果；
图9为对合成地震记录利用二维字典的形态成分分析得到的面波；
图10为对图4所示的实际地震记录利用二维字典的形态成分分析压制面波后的结果；
图11为对实际地震记录利用二维字典的形态成分分析得到的面波；
图12是本发明实施例的基于形态成分分析的面波分离系统的结构示意图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本发明，并不用于限定本发明。
请参阅图1，是本发明实施例的基于形态成分分析的面波分离方法的流程图。本发明实施例的基于形态成分分析的面波分离方法包括以下步骤：
步骤100：制作合成地震记录，并构造基于二维字典形态成分分析的面波分离模型；
在步骤100中，本发明实施例制作合成地震记录具体为：合成三层反射波地震记录，由面波与反射波两种信号分量组合而成。反射波信号分量由主频为35Hz的雷克子波与反射系数褶积而成，最大振幅为1；面波信号分量由正弦函数构成的扫描信号合成，频带范围为5-15Hz，最大振幅为3。假设合成地震记录大小为N*N，该地震记录可表示为：
x＝xg+xr+n (1)
式中，xg为面波信号分量，xr为反射波信号分量，n为随机噪声。
构造出两种二维字典Tg、Tr分别稀疏表示面波信号分量与反射波信号分量，这样地震记录x在Tg和Tr组成的联合字典下的稀疏表示问题可以表示为求解下式：
式(2)即为面波分离模型，其中λ为拉格朗日乘子，控制着面波信号分量与反射波信号分量之和拟合地震记录的失真度。通过求解该式不仅可以得到面波信号分量与反射波信号分量的表示系数，而且可以有效抑制地震资料中存在的随机噪声干扰。
步骤200：选取二维非抽样小波字典稀疏表示面波信号分量，并选取二维局部离散余弦字典稀疏表示反射波信号分量；
在步骤200中，根据面波信号低频、窄带及频散等特性，本发明选取二维非抽样小波字典稀疏表示面波信号分量，二维非抽样小波变换可以对面波信号分量进行多尺度的分析，经过二维非抽样小波变换后的面波信号分量主要集中在低频高波数子带内，因此对于面波信号分量具有较好的辨识能力。本发明选取的二维非抽样离散小波变换是一种超完备的小波变换，它克服了正交小波变换对平移敏感的缺点，具有平移不变性的特点。
根据反射波信号局部表现为波动信号并且相关性强等特性，本发明选取二维局部离散余弦字典稀疏表示反射波信号，二维离散余弦变换对相关性较强的信号具有能量集中的作用，它是加窗型的离散余弦变换，能够有效表示信号的局部特征，非常适合用来稀疏表示局部具有波动信号特征的反射波信号。本发明通过实验验证所选取的字典合理性，具体请一并参阅图2、图3、图4和图5，图2为实际地震记录切片效果图，从图2中可以明显看出反射波信号成分。图3为图2所示的地震记录在二维非抽样离散小波字典(UDWT2)与二维局部离散余弦字典(LDCT2)表示中，稀疏性与重构相对误差的关系比较图；从图3中可以看出二维非抽样小波字典对于反射波信号的稀疏表示能力较弱，二维离散余弦字典对于反射波信号的稀疏表示能力较强。图4为含有强能量面波干扰的实际地震记录图，图5为图4所示的地震记录在二维非抽样离散小波字典与二维局部离散余弦字典表示中，稀疏性与重构相对误差的关系比较图。从图5中可以明显看出扇形状分布的面波信号，可以看出二维非抽样小波字典对面波信号的稀疏表示能力较强，二维离散余弦字典对面波信号的稀疏表示能力较弱。
步骤300：采用二维字典下形态成分分析对面波信号分量与反射波信号分量进行分离；
在步骤300中，采用二维字典下形态成分分析对面波信号分量与反射波信号分量进行分离具体包括：将地震记录分为面波信号分量与反射波信号分量，将二维非抽样离散小波字典与二维局部离散余弦字典作为地震记录的联合稀疏表示字典，构建稀疏表示模型；并通过块协调松弛算法对稀疏表示模型进行求解，得到分离出来的面波信号分量xg以及反射波信号分量xr；其中，稀疏表示模型通过对式(2)进行转换得到：
在使用块协调松弛算法分离面波信号分量xg以及反射波信号分量xr时，由于面波信号分量主要集中在信号小波变换的低频高波数域内，并且该域内反射波信号能量较低，做出以下修改：在对面波信号分量进行阈值处理阶段，对面波信号分量集中的低频高波数域内系数不进行阈值处理，将其保留下来；这样使得分离面波与反射波的效果更加明显。通过块协调松弛算法对稀疏表示模型进行求解的具体步骤如下：
(1)初始化拉格朗日乘子λ，迭代次数N，阈值δ＝Lmax，截止阈值δmin，其中Lmax为地震记录x在两种字典下表示系数的最大值中较小的一个；
(2)进行N次迭代：
A部分：假设xr不变，更新xg
①计算残差：r＝x-xg-xr；
②计算xg+r经过二维非抽样离散小波变换后的系数，
③保留系数αg低频高波数域中系数不变，对系数αg的其它域小波系数进行软阈值处理，阈值为δ，得到系数
④通过重建xg。
B部分：假设xg不变，更新xr
⑤计算残差r＝x-xg-xr；
⑥计算xr+r经过二维局部离散余弦变换后的系数，
⑦对系数αr进行软阈值处理，阈值为δ，得到系数
⑧通过重建xr。
(3)更新阈值δ＝δ*λ；
(4)如果δ>δmin转至步骤(2)，否则结束。
为了进一步说明本发明的有效性和准确性，本发明对合成地震记录以及实际地震记录进行面波与反射波信号的分离实验，具体请一并参阅图6、图7、图8、图9、图10和图11，图6为合成的反射波地震记录，图7为合成的含有面波与反射波的地震记录，图8为对图7所示的合成地震记录利用二维字典的形态成分分析压制面波后的结果，从图8可以看出，面波大部分能够得到压制；图9为对合成地震记录利用二维字典的形态成分分析得到的面波；图10为对图4所示的实际地震记录利用二维字典的形态成分分析压制面波后的结果，面波成分得到较好的压制，反射波信号同相轴清晰地显示出来。图11为对实际地震记录利用二维字典的形态成分分析得到的面波。
请参阅图12，是本发明实施例的基于形态成分分析的面波分离系统的结构示意图。本发明实施例的基于形态成分分析的面波分离系统包括地震记录制作模块、二维字典选取模块和面波分离模块；其中，
地震记录制作模块用于制作合成地震记录，并构造基于二维字典形态成分分析的面波分离模型；本发明实施例制作合成地震记录具体为：合成三层反射波地震记录，由面波与反射波两种信号分量组合而成。反射波信号分量由主频为35Hz的雷克子波与反射系数褶积而成，最大振幅为1；面波信号分量由正弦函数构成的扫描信号合成，频带范围为5-15Hz，最大振幅为3。假设合成地震记录大小为N*N，该地震记录可表示为：
x＝xg+xr+n (1)
式中，xg为面波信号分量，xr为反射波信号分量，n为随机噪声。
构造出两种二维字典Tg、Tr分别稀疏表示面波信号分量与反射波信号分量，这样地震记录x在Tg和Tr组成的联合字典下的稀疏表示问题可以表示为求解下式：
式(2)中λ为拉格朗日乘子，控制着面波信号分量与反射波信号分量之和拟合地震记录的失真度。通过求解该式不仅可以得到面波信号分量与反射波信号分量的表示系数，而且可以有效抑制地震资料中存在的随机噪声干扰。
二维字典选取模块用于选取二维非抽样小波字典稀疏表示面波信号分量，并选取二维局部离散余弦字典稀疏表示反射波信号分量；在本发明实施例中，根据面波信号低频、窄带及频散等特性，选取二维非抽样小波字典稀疏表示面波信号分量，二维非抽样小波变换可以对面波信号分量进行多尺度的分析，经过二维非抽样小波变换后的面波信号分量主要集中在低频高波数子带内，因此对于面波信号分量具有较好的辨识能力。本发明选取的二维非抽样离散小波变换是一种超完备的小波变换，它克服了正交小波变换对平移敏感的缺点，具有平移不变性的特点。
根据反射波信号局部表现为波动信号并且相关性强等特性，选取二维局部离散余弦字典稀疏表示反射波信号，二维离散余弦变换对相关性较强的信号具有能量集中的作用，它是加窗型的离散余弦变换，能够有效表示信号的局部特征，非常适合用来稀疏表示局部具有波动信号特征的反射波信号。本发明通过实验验证所选取的字典合理性，具体请一并参阅图2、图3、图4和图5，图2为实际地震记录切片效果图，从图2中可以明显看出反射波信号成分。图3为图2所示的地震记录在二维非抽样离散小波字典(UDWT2)与二维局部离散余弦字典(LDCT2)表示中，稀疏性与重构相对误差的关系比较图；从图3中可以看出二维非抽样小波字典对于反射波信号的稀疏表示能力较弱，二维离散余弦字典对于反射波信号的稀疏表示能力较强。图4为含有强能量面波干扰的实际地震记录图，图5为图4所示的地震记录在二维非抽样离散小波字典与二维局部离散余弦字典表示中，稀疏性与重构相对误差的关系比较图；从图5中可以明显看出扇形状分布的面波信号，可以看出二维非抽样小波字典对面波信号的稀疏表示能力较强，二维离散余弦字典对面波信号的稀疏表示能力较弱。
面波分离模块用于采用二维字典下形态成分分析对面波信号分量与反射波信号分量进行分离；其中，采用二维字典下形态成分分析对面波信号分量与反射波信号分量进行分离具体包括：将地震记录分为面波信号分量与反射波信号分量，将二维非抽样离散小波字典与二维局部离散余弦字典作为地震记录的联合稀疏表示字典，构建稀疏表示模型；并通过块协调松弛算法对稀疏表示模型进行求解，得到分离出来的面波信号分量xg以及反射波信号分量xr；其中，稀疏表示模型通过对式(2)进行转换得到：
在使用块协调松弛算法分离面波信号分量xg以及反射波信号分量xr时，由于面波信号分量主要集中在信号小波变换的低频高波数域内，并且该域内反射波信号能量较低，做出以下修改：在对面波信号分量进行阈值处理阶段，对面波信号分量集中的低频高波数域内系数不进行阈值处理，将其保留下来；这样使得分离面波与反射波的效果更加明显。通过块协调松弛算法对稀疏表示模型进行求解的具体步骤如下：
(1)初始化拉格朗日乘子λ，迭代次数N，阈值δ＝Lmax，截止阈值δmin，其中Lmax为地震记录x在两种字典下表示系数的最大值中较小的一个；
(2)进行N次迭代：
A部分：假设xr不变，更新xg
①计算残差：r＝x-xg-xr；
②计算xg+r经过二维非抽样离散小波变换后的系数，
③保留系数αg低频高波数域中系数不变，对系数αg的其它域小波系数进行软阈值处理，阈值为δ，得到系数
④通过重建xg。
B部分：假设xg不变，更新xr
⑤计算残差r＝x-xg-xr；
⑥计算xr+r经过二维局部离散余弦变换后的系数，
⑦对系数αr进行软阈值处理，阈值为δ，得到系数
⑧通过重建xr。
(3)更新阈值δ＝δ*λ；
(4)如果δ>δmin转至步骤(2)，否则结束。
本发明实施例的基于形态成分分析的面波分离方法及系统根据地震记录中面波信号与反射波信号的结构特征差异，分别选取二维非抽样小波字典稀疏表示面波信号，选取二维局部离散余弦字典稀疏表示反射波信号，相对于单一变换，在分离面波与反射波的过程中，既考虑面波信号特性的同时也考虑到反射波信号的特性，从而能更加有效地分离出面波与反射波信号。同时，本发明利用二维字典对二维地震数据进行整体稀疏表示与处理，充分利用了道与道之间的相关性，得到较好的分离效果。并且面波信号与反射波信号在二维字典下具有更好的辨识度，同时考虑在二维非抽样小波字典下，面波信号集中在低频高波数域中，使得面波信号的分离效果更好。
以上所述仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等，均应包含在本发明的保护范围之内。
Claims (4)
1.一种基于形态成分分析的面波分离方法，包括：
步骤a：制作合成地震记录，并构造基于二维字典形态成分分析的面波分离模型；
步骤b：选取二维非抽样小波字典稀疏表示面波信号分量，并选取二维局部离散余弦字典稀疏表示反射波信号分量；
步骤c：采用二维字典下形态成分分析对面波信号分量与反射波信号分量进行分离；
在所述步骤a中，所述合成地震记录为：
x＝xg+xr+n (1)
式(1)中，xg为面波信号分量，xr为反射波信号分量，n为随机噪声；
在所述步骤a中，所述面波分离模型为：
式(2)中，λ为拉格朗日乘子，控制着面波信号分量与反射波信号分量之和拟合地震记录的失真度；Tg、Tr分别稀疏表示面波信号分量与反射波信号分量；
通过块协调松弛算法对稀疏表示模型进行求解，得到分离出来的面波信号分量以及反射波信号分量；具体包括：
(1)初始化拉格朗日乘子λ，迭代次数N，阈值δ＝Lmax，截止阈值δmin，其中Lmax为地震记录x在两种字典下表示系数的最大值中较小的一个；
(2)进行N次迭代：
A部分：假设xr不变，更新xg
①计算残差：r＝x-xg-xr；
②计算xg+r经过二维非抽样离散小波变换后的系数，
③保留系数αg低频高波数域中系数不变，对系数αg的其它域小波系数进行软阈值处理，阈值为δ，得到系数
④通过重建xg；
B部分：假设xg不变，更新xr
⑤计算残差r＝x-xg-xr；
⑥计算xr+r经过二维局部离散余弦变换后的系数，
⑦对系数αr进行软阈值处理，阈值为δ，得到系数
⑧通过重建xr；
(3)更新阈值δ＝δ*λ；
(4)如果δ＞δmin转至步骤(2)，否则结束。
2.根据权利要求1所述的基于形态成分分析的面波分离方法，其特征在于，在所述步骤c中，所述采用二维字典下形态成分分析对面波信号分量与反射波信号分量进行分离还包括：将地震记录分为面波信号分量与反射波信号分量，将二维非抽样离散小波字典与二维局部离散余弦字典作为地震记录的联合稀疏表示字典，构建稀疏表示模型；所述稀疏表示模型通过对式(2)进行转换得到：
3.一种基于形态成分分析的面波分离系统，其特征在于，包括地震记录制作模块、二维字典选取模块和面波分离模块；
所述地震记录制作模块用于制作合成地震记录，并构造基于二维字典形态成分分析的面波分离模型；
所述二维字典选取模块用于选取二维非抽样小波字典稀疏表示面波信号分量，并选取二维局部离散余弦字典稀疏表示反射波信号分量；
所述合成地震记录为：
x＝xg+xr+n (1)
式(1)中，xg为面波信号分量，xr为反射波信号分量，n为随机噪声；
所述面波分离模型为：
式(2)中，λ为拉格朗日乘子，控制着面波信号分量与反射波信号分量之和拟合地震记录的失真度；Tg、Tr分别稀疏表示面波信号分量与反射波信号分量；
所述面波分离模块用于采用二维字典下形态成分分析对面波信号分量与反射波信号分量进行分离；通过块协调松弛算法对稀疏表示模型进行求解，得到分离出来的面波信号分量以及反射波信号分量；所述通过块协调松弛算法对稀疏表示模型进行求解具体包括：
(1)初始化拉格朗日乘子λ，迭代次数N，阈值δ＝Lmax，截止阈值δmin，其中Lmax为地震记录x在两种字典下表示系数的最大值中较小的一个；
(2)进行N次迭代：
A部分：假设xr不变，更新xg
①计算残差：r＝x-xg-xr；
②计算xg+r经过二维非抽样离散小波变换后的系数，
③保留系数αg低频高波数域中系数不变，对系数αg的其它域小波系数进行软阈值处理，阈值为δ，得到系数
④通过重建xg；
B部分：假设xg不变，更新xr
⑤计算残差r＝x-xg-xr；
⑥计算xr+r经过二维局部离散余弦变换后的系数，
⑦对系数αr进行软阈值处理，阈值为δ，得到系数
⑧通过重建xr；
(3)更新阈值δ＝δ*λ；
(4)如果δ＞δmin转至步骤(2)，否则结束。
4.根据权利要求3所述的基于形态成分分析的面波分离系统，其特征在于，所述面波分离模块采用二维字典下形态成分分析对面波信号分量与反射波信号分量进行分离还包括：将地震记录分为面波信号分量与反射波信号分量，将二维非抽样离散小波字典与二维局部离散余弦字典作为地震记录的联合稀疏表示字典，构建稀疏表示模型；所述稀疏表示模型通过对式(2)进行转换得到：
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510578757.3A CN105182417B (zh) | 2015-09-11 | 2015-09-11 | 一种基于形态成分分析的面波分离方法及系统 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510578757.3A CN105182417B (zh) | 2015-09-11 | 2015-09-11 | 一种基于形态成分分析的面波分离方法及系统 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105182417A CN105182417A (zh) | 2015-12-23 |
| CN105182417B true CN105182417B (zh) | 2018-06-22 |
Family
ID=54904613
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510578757.3A Expired - Fee Related CN105182417B (zh) | 2015-09-11 | 2015-09-11 | 一种基于形态成分分析的面波分离方法及系统 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105182417B (zh) |
Families Citing this family (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105700015B (zh) * | 2016-02-02 | 2017-07-28 | 中国矿业大学(北京) | 一种小尺度不连续地质体检测方法和装置 |
| CN107065006B (zh) * | 2017-01-23 | 2019-06-11 | 武汉大学 | 一种基于在线字典更新的地震信号编码方法 |
| CN107083957B (zh) * | 2017-05-25 | 2020-05-26 | 电子科技大学 | 钻井液随钻信号的泵冲干扰消除方法及系统 |
| CN109782346B (zh) * | 2019-01-14 | 2020-07-28 | 西安交通大学 | 一种基于形态成分分析的采集脚印压制方法 |
| CN111308555B (zh) * | 2020-03-11 | 2021-04-13 | 清华大学 | 一种基于频间相似性低频重构的面波干扰消除方法 |
| CN113077386A (zh) * | 2021-04-06 | 2021-07-06 | 电子科技大学 | 基于字典学习和稀疏表征的地震资料高分辨率处理方法 |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4935905A (en) * | 1989-08-04 | 1990-06-19 | Terra Linda Canada, Inc. | Method and apparatus for seismic exploration |
| CN102103215A (zh) * | 2009-12-16 | 2011-06-22 | 中国石油天然气集团公司 | 三维高密度地震勘探记录叠前面波压制方法 |
| CN102944895A (zh) * | 2012-10-26 | 2013-02-27 | 中国石油大学(北京) | 一种地震记录面波干扰压制方法及装置 |
| CN104570115A (zh) * | 2013-10-16 | 2015-04-29 | 中国石油天然气集团公司 | 一种面波衰减方法及装置 |
| CN104614769A (zh) * | 2015-02-05 | 2015-05-13 | 中铁隧道集团有限公司 | 一种压制地震面波的聚束滤波方法 |
-
2015
- 2015-09-11 CN CN201510578757.3A patent/CN105182417B/zh not_active Expired - Fee Related
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4935905A (en) * | 1989-08-04 | 1990-06-19 | Terra Linda Canada, Inc. | Method and apparatus for seismic exploration |
| CN102103215A (zh) * | 2009-12-16 | 2011-06-22 | 中国石油天然气集团公司 | 三维高密度地震勘探记录叠前面波压制方法 |
| CN102944895A (zh) * | 2012-10-26 | 2013-02-27 | 中国石油大学(北京) | 一种地震记录面波干扰压制方法及装置 |
| CN104570115A (zh) * | 2013-10-16 | 2015-04-29 | 中国石油天然气集团公司 | 一种面波衰减方法及装置 |
| CN104614769A (zh) * | 2015-02-05 | 2015-05-13 | 中铁隧道集团有限公司 | 一种压制地震面波的聚束滤波方法 |
Non-Patent Citations (1)
| Title |
|---|
| 基于稀疏表示理论的地震信号处理方法研究;李海山;《中国博士学位论文全文数据库 基础科学辑》;20150615(第06期);第55-57、66-77页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105182417A (zh) | 2015-12-23 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105182417B (zh) | 一种基于形态成分分析的面波分离方法及系统 | |
| Chen et al. | Empirical low-rank approximation for seismic noise attenuation | |
| Liu et al. | Random noise suppression in seismic data: What can deep learning do? | |
| CN108845352B (zh) | 基于vmd近似熵与多层感知机的沙漠地震信号去噪方法 | |
| CN107144880B (zh) | 一种地震波波场分离方法 | |
| Chen et al. | Random noise reduction using a hybrid method based on ensemble empirical mode decomposition | |
| CN105974468B (zh) | 一种能够同时进行五维地震数据重建和噪声压制的方法 | |
| Zhou et al. | Robust noise attenuation based on nuclear norm minimization and a trace prediction strategy | |
| Wang et al. | An iterative curvelet thresholding algorithm for seismic random noise attenuation | |
| WO2007041720A2 (en) | Coherent wave energy removal from seismic data | |
| Li et al. | Wavelet-based higher order correlative stacking for seismic data denoising in the curvelet domain | |
| US20150117151A1 (en) | Automatic Tracking of Faults by Slope Decomposition | |
| CN108828670B (zh) | 一种地震数据降噪方法 | |
| Sun et al. | A noise attenuation method for weak seismic signals based on compressed sensing and CEEMD | |
| CN106199532A (zh) | 基于混合傅立叶‑小波分析的探地雷达信号降噪方法 | |
| Zhang et al. | 3D simultaneous seismic data reconstruction and noise suppression based on the curvelet transform | |
| CN110031899A (zh) | 基于压缩感知的弱信号提取算法 | |
| Zheng et al. | The surface wave suppression using the second generation curvelet transform | |
| Yi et al. | A least-squares correlation-based full traveltime inversion for shallow subsurface velocity reconstruction | |
| Bai et al. | Seismic signal enhancement based on the low‐rank methods | |
| CN105676292B (zh) | 一种基于二维曲波变换的三维地震数据去噪方法 | |
| Sun et al. | Multiple attenuation using λ-f domain high-order and high-resolution Radon transform based on SL0 norm | |
| CN102854530B (zh) | 基于对数时频域双曲平滑的动态反褶积方法 | |
| CN103645504A (zh) | 基于广义瞬时相位及p范数负模的地震弱信号处理方法 | |
| CN109782346B (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 | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20180622 Termination date: 20190911 |
|
| CF01 | Termination of patent right due to non-payment of annual fee |