CN105182429A - 一种海洋可控源电磁数据海水扰动噪声校正方法 - Google Patents
一种海洋可控源电磁数据海水扰动噪声校正方法 Download PDFInfo
- Publication number
- CN105182429A CN105182429A CN201510628332.9A CN201510628332A CN105182429A CN 105182429 A CN105182429 A CN 105182429A CN 201510628332 A CN201510628332 A CN 201510628332A CN 105182429 A CN105182429 A CN 105182429A
- Authority
- CN
- China
- Prior art keywords
- source electromagnetic
- electromagnetic data
- ocean
- seawater
- controllable source
- 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
Classifications
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/30—Assessment of water resources
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明涉及一种海洋可控源电磁数据海水扰动噪声校正方法。首先,根据海洋可控源电磁数据与海水扰动噪声特性选取适当的小波基,利用小波多分辨率分析原理,根据采样频率与采样时间决定分解级数,然后进行小波分解,利用最高分解级相对应的近似系数,重构出海水扰动噪声估计信号,最后将估计的海水扰动噪声从海洋可控源电磁原始数据中去除,即可以得到校正后的信号。本发明所述方法利用小波最高分解级相对应的近似系数,重构出海水扰动噪声的估计信号,达到校正海洋可控源电磁信号的目的,可以有效提升海洋可控源电磁数据的定性显示分析及定量反演解释效果。
Description
技术领域
本发明涉及一种基于小波多分辨率分析的海水扰动噪声校正方法,特别适用于海洋可控源电磁数据的海水扰动噪声校正方法。
背景技术
海洋可控源电磁法(MarineControlled-SourceElectectroMagneticMethod,MCSEM)是勘察海底油气资源和矿产资源的一种新兴海洋地球物理勘探技术。
海水能够屏蔽空间电磁噪声,有利于应用MCSEM探测海底高阻异常体。但是,MCSEM信号依然会受到随机噪声、空气波及海水扰动噪声的影响,导致信号失真,直接影响电磁数据的定性显示分析及定量反演解释效果。
目前,国内外的MCSEM数据处理方法主要有:基于最小二乘拟合与叠加算法的数据压缩及噪声估计法、基于波场分离理论的空气波去除技术、一阶差分预白及后黑法、Butterworth滤波器法、EMD消噪法、小波模极大值法、小波阈值法等。
但是,以上MCSEM数据处理方法主要针对随机噪声与空气波的去除,不利于消除海浪、暗流等运动产生的感应电磁场,即海水扰动噪声。
发明内容
本发明的目的是针对海水扰动噪声,提供一种基于小波多分辨率分析的海洋可控源电磁数据海水扰动噪声校正方法,其特征在于,该方法包括如下步骤:
(1)根据MCSEM数据与海水扰动噪声特性选取适当的小波基;
(2)利用小波多分辨率分析原理,根据MCSEM数据的采样频率fs与采样时间t决定小波分解级数L;
(3)对MCSEM数据进行L级小波分解,利用尺度L上的近似系数AL,重构出海水扰动噪声的估计信号X;
(4)从MCSEM数据中去除X,即可以得到校正后的信号。
所述步骤(1)具体为选取db8小波基。
所述步骤(2)具体为按照如下公式进行小波分解级数L的确定:
式中,t表示MCSEM数据的采样时间,fs表示MCSEM数据的采样频率,lw表示低通分解滤波器的长度,fix表示向零取整函数。
所述步骤(3)具体为采用db8小波,分解L层,得到L个细节系数D1~DL,与1个近似系数AL,利用AL重构出海水扰动噪声的估计信号X。
所述步骤(4)具体为采用如下公式得到校正后的信号sig1:
sig1=sig-X
式中,sig表示原信号,X表示利用AL重构出的海水扰动噪声估计信号。
本发明根据小波多分辨率分析原理,对MCSEM数据进行分解,利用小波最高分解级相对应的近似系数重构出海水扰动噪声的估计信号,达到校正MCSEM信号的目的。
附图说明:
图1是本发明一种海洋可控源电磁数据海水扰动噪声校正方法流程图;
图2是本发明的海底地层模型图;
图3是本发明海洋可控源电磁数据理论信号图;
图4是本发明海洋可控源电磁数据模拟噪声信号图;
图5是本发明重构出的海水扰动噪声估计信号图;
图6是本发明校正海水扰动噪声后的信号图。
具体实施方式:
本发明提出的一种海洋可控源电磁数据海水扰动噪声校正方法结合附图及实施例详细说明如下:
本发明方法的技术特点为,根据MCSEM数据特点与海水扰动噪声选取适当的小波基,根据采样频率与采样时间决定分解级数,然后进行小波分解,利用最高分解级相对应的近似系数,重构出海水扰动噪声的估计信号,最后将海水扰动噪声估计信号从MCSEM原始数据中去除,以达到校正MCSEM信号的目的,图1为校正方法流程图。
本发明具体实施例所采用的4层模型如图2所示,采用KerryKey's等人编写的一维海洋可控源电磁正演代码(Dipole1D)来计算理论响应信号,设定发射机发射频率为0.08Hz,发射机在海底上方50米,沿y轴方向移动,接收点位于发射机航迹正下方,接收机的采样频率fs为10Hz,采样时间t为5h(300min),共180000个数据点,截取其中第100分钟到第300分钟的数据,计算的理论信号如图3所示,在理论信号上叠加正弦低频信号模拟海水扰动噪声,模拟的噪声信号如图4所示,基于小波多分辨率分析的海水扰动噪声校正方法具体包括如下步骤:
1)根据MCSEM数据与海水扰动类噪声特点选取db8小波基;
2)利用小波多分辨率分析原理,根据公式:
式中,t表示MCSEM数据的采样时间,fs表示MCSEM数据的采样频率,lw表示低通分解滤波器的长度,fix表示向零取整函数,计算得到小波分解级数L为9;
3)对MCSEM数据采用db8小波基,进行9级小波分解,计算得到9尺度上的近似系数A9,利用A9重构出海水扰动噪声的估计信号X,X如图5所示;
4)利用如下公式:
sig1=sig-X
式中,sig表示原信号,X表示利用AL重构出的海水扰动噪声估计信号,从MCSEM数据中去除海水扰动噪声估计信号,即可以得到校正后的信号sig1,sig1如图6所示。
Claims (5)
1.一种海洋可控源电磁数据海水扰动噪声校正方法,其特征在于,该方法包括如下步骤:
(1)根据海洋可控源电磁数据与海水扰动噪声特性选取适当的小波基;
(2)利用小波多分辨率分析原理,根据海洋可控源电磁数据的采样频率fs与采样时间t决定小波分解级数L;
(3)对海洋可控源电磁数据进行L级小波分解,利用尺度L上的近似系数AL,重构出海水扰动噪声的估计信号X;
(4)从海洋可控源电磁数据中去除X,即可以得到校正后的信号。
2.按照权利要求1所述的海洋可控源电磁数据海水扰动噪声校正方法,其特征在于,所述步骤(1)选取db8小波基。
3.按照权利要求1所述的海洋可控源电磁数据海水扰动噪声校正方法,其特征在于,所述步骤(2)按照如下公式进行小波分解级数L的确定:
式中,t表示海洋可控源电磁数据的采样时间,fs表示海洋可控源电磁数据的采样频率,lw表示低通分解滤波器的长度,fix表示向零取整函数。
4.按照权利要求1所述的海洋可控源电磁数据海水扰动噪声校正方法,其特征在于,所述步骤(3)采用db8小波,分解L层,得到L个细节系数D1~DL,与1个近似系数AL,利用AL重构出海水扰动噪声的估计信号X。
5.按照权利要求1所述的海洋可控源电磁数据海水扰动噪声校正方法,其特征在于,所述步骤(4)采用如下公式得到校正后的信号sig1:
sig1=sig-X
式中,sig表示海洋可控源电磁数据原始信号,X表示利用AL重构出的海水扰动噪声估计信号。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510628332.9A CN105182429B (zh) | 2015-09-29 | 2015-09-29 | 一种海洋可控源电磁数据海水扰动噪声校正方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510628332.9A CN105182429B (zh) | 2015-09-29 | 2015-09-29 | 一种海洋可控源电磁数据海水扰动噪声校正方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105182429A true CN105182429A (zh) | 2015-12-23 |
CN105182429B CN105182429B (zh) | 2017-09-29 |
Family
ID=54904625
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510628332.9A Active CN105182429B (zh) | 2015-09-29 | 2015-09-29 | 一种海洋可控源电磁数据海水扰动噪声校正方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105182429B (zh) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108680966A (zh) * | 2018-03-21 | 2018-10-19 | 中国石油大学(华东) | 海洋可控源电磁勘探噪声降噪效果评估方法 |
CN108983299A (zh) * | 2018-05-30 | 2018-12-11 | 吉林大学 | 一种面向海洋可控源电磁数据消噪的新型小波基构造方法 |
CN109413982A (zh) * | 2018-12-11 | 2019-03-01 | 南京信息工程大学 | 一种海水电磁屏蔽体 |
CN110989005A (zh) * | 2019-12-23 | 2020-04-10 | 西北工业大学青岛研究院 | 基于标量磁力仪阵列的弱磁异常自适应实时探测方法 |
CN112198556A (zh) * | 2020-09-16 | 2021-01-08 | 南方科技大学 | 信号处理方法、装置、终端设备及存储介质 |
CN115903052A (zh) * | 2022-09-07 | 2023-04-04 | 航天恒星科技有限公司 | 基于扩展Prony算法的舰船尾迹电磁信号检测重构方法 |
CN116029131A (zh) * | 2023-01-09 | 2023-04-28 | 盛东如东海上风力发电有限责任公司 | 一种海上风电场水下金属结构磁信号模拟和处理方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20090133870A1 (en) * | 2005-11-01 | 2009-05-28 | Pavlov Dmitriy A | Method for Phase and Amplitude Correction in Controlled Source Electromagnetic Survey Data |
US20100109671A1 (en) * | 2008-11-03 | 2010-05-06 | Bruce Alan Hobbs | Method for acquiring controlled source electromagnetic survey data to assist in attenuating correlated noise |
CN101957458A (zh) * | 2009-07-21 | 2011-01-26 | 北京大学 | 一种处理超低频数据的方法 |
CN103543469A (zh) * | 2012-07-17 | 2014-01-29 | 中国石油化工股份有限公司 | 一种基于小波变换的小尺度阈值去噪方法 |
-
2015
- 2015-09-29 CN CN201510628332.9A patent/CN105182429B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20090133870A1 (en) * | 2005-11-01 | 2009-05-28 | Pavlov Dmitriy A | Method for Phase and Amplitude Correction in Controlled Source Electromagnetic Survey Data |
US20100109671A1 (en) * | 2008-11-03 | 2010-05-06 | Bruce Alan Hobbs | Method for acquiring controlled source electromagnetic survey data to assist in attenuating correlated noise |
CN101957458A (zh) * | 2009-07-21 | 2011-01-26 | 北京大学 | 一种处理超低频数据的方法 |
CN103543469A (zh) * | 2012-07-17 | 2014-01-29 | 中国石油化工股份有限公司 | 一种基于小波变换的小尺度阈值去噪方法 |
Non-Patent Citations (4)
Title |
---|
周潞 等: "小波收缩阈值消噪在海洋电磁中的应用", 《信息通信》 * |
李泽林 等: "海洋CESM信号去噪方法研究", 《中国地球科学联合学术年会2014》 * |
李肃义 等: "电性源时域地空电磁数据小波去噪方法研究", 《地球物理学报》 * |
马海舲: "海洋可控源电磁勘探的数据处理与解释", 《中国优秀硕士学位论文全文数据库·基础科学辑》 * |
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108680966A (zh) * | 2018-03-21 | 2018-10-19 | 中国石油大学(华东) | 海洋可控源电磁勘探噪声降噪效果评估方法 |
CN108983299A (zh) * | 2018-05-30 | 2018-12-11 | 吉林大学 | 一种面向海洋可控源电磁数据消噪的新型小波基构造方法 |
CN109413982A (zh) * | 2018-12-11 | 2019-03-01 | 南京信息工程大学 | 一种海水电磁屏蔽体 |
CN110989005A (zh) * | 2019-12-23 | 2020-04-10 | 西北工业大学青岛研究院 | 基于标量磁力仪阵列的弱磁异常自适应实时探测方法 |
CN110989005B (zh) * | 2019-12-23 | 2021-12-28 | 西北工业大学青岛研究院 | 基于标量磁力仪阵列的弱磁异常自适应实时探测方法 |
CN112198556A (zh) * | 2020-09-16 | 2021-01-08 | 南方科技大学 | 信号处理方法、装置、终端设备及存储介质 |
WO2022057305A1 (zh) * | 2020-09-16 | 2022-03-24 | 南方科技大学 | 信号处理方法、装置、终端设备及存储介质 |
CN115903052A (zh) * | 2022-09-07 | 2023-04-04 | 航天恒星科技有限公司 | 基于扩展Prony算法的舰船尾迹电磁信号检测重构方法 |
CN115903052B (zh) * | 2022-09-07 | 2023-12-19 | 航天恒星科技有限公司 | 基于扩展Prony算法的舰船尾迹电磁信号检测重构方法 |
CN116029131A (zh) * | 2023-01-09 | 2023-04-28 | 盛东如东海上风力发电有限责任公司 | 一种海上风电场水下金属结构磁信号模拟和处理方法 |
CN116029131B (zh) * | 2023-01-09 | 2024-02-23 | 盛东如东海上风力发电有限责任公司 | 一种海上风电场水下金属结构磁信号模拟和处理方法 |
Also Published As
Publication number | Publication date |
---|---|
CN105182429B (zh) | 2017-09-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105182429A (zh) | 一种海洋可控源电磁数据海水扰动噪声校正方法 | |
CN104793253B (zh) | 基于数学形态学的航空电磁数据去噪方法 | |
Sun et al. | Cross-correlation analysis and time delay estimation of a homologous micro-seismic signal based on the Hilbert–Huang transform | |
CN108490494B (zh) | 基于谱减法及小波分析的海洋磁测噪声抑制方法 | |
CN101598809A (zh) | 一种自适应消除线性规则噪声以及多次波干扰的方法 | |
CN109581516B (zh) | 曲波域统计量自适应阈值探地雷达数据去噪方法及系统 | |
CN108680966B (zh) | 海洋可控源电磁勘探噪声降噪效果评估方法 | |
Ma et al. | Deep residual encoder–decoder networks for desert seismic noise suppression | |
CN105116442A (zh) | 岩性油气藏弱反射地震信号的重构方法 | |
CN103713324B (zh) | 一种基于局部方差分析的自适应曲波阈值地震去噪方法 | |
Li et al. | Wavelet-based higher order correlative stacking for seismic data denoising in the curvelet domain | |
CN112198556A (zh) | 信号处理方法、装置、终端设备及存储介质 | |
Wei et al. | Comparative research on noise reduction of transient electromagnetic signals based on empirical mode decomposition and variational mode decomposition | |
Chen et al. | Seawater motion-induced electromagnetic noise reduction in marine magnetotelluric data using current meters | |
CN102236109B (zh) | 一种重磁干扰区的方差系数干扰校正方法 | |
CN109188542B (zh) | 一种波区相关性检测的远参考大地电磁阻抗计算方法 | |
CN107632326A (zh) | 地球物理信号去噪方法 | |
CN117574062A (zh) | 基于vmd-dnn模型的小回线瞬变电磁信号去噪方法 | |
CN108919345B (zh) | 一种海底电缆陆检噪声的衰减方法 | |
CN112713907B (zh) | 一种基于字典学习的海洋csem降噪方法及系统 | |
CN113397533B (zh) | 一种弱生命信号提取方法、装置、电子设备和存储介质 | |
Jeng et al. | A nonlinear method of removing harmonic noise in geophysical data | |
CN111856590B (zh) | 海洋大地电磁探测的海浪磁干扰压制方法 | |
CN109885906B (zh) | 一种基于粒子群优化的磁共振测深信号稀疏消噪方法 | |
CN108519619B (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 |