CN115542408B - 一种海洋瞬变电磁数据预处理与快速定区成像方法 - Google Patents
一种海洋瞬变电磁数据预处理与快速定区成像方法 Download PDFInfo
- Publication number
- CN115542408B CN115542408B CN202211545726.4A CN202211545726A CN115542408B CN 115542408 B CN115542408 B CN 115542408B CN 202211545726 A CN202211545726 A CN 202211545726A CN 115542408 B CN115542408 B CN 115542408B
- Authority
- CN
- China
- Prior art keywords
- data
- filtering
- delay
- time
- processing
- 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.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V3/00—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
- G01V3/38—Processing data, e.g. for analysis, for interpretation, for correction
-
- 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
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Image Processing (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
Description
技术领域
本发明涉及电磁数据处理技术领域,尤其是一种海洋瞬变电磁数据预处理与快速定区成像方法。
背景技术
海洋可控源电磁法(MCSEM)在海洋油气勘探、天然气水合物勘探等方面有诸多重要应用。由于MCSEM在油气勘探上的成功应用,石油工业界对相对成熟的频率域海洋可控源电磁法进行了大力推广,而时间域方法的研究相对较少。然而,随着海洋资源勘探以及浅海相关作业的需求,海洋电磁法有往浅海应用偏移的趋势,一维数值模拟研究表明,浅海时间域可控源电磁法相比频率域所受空气波影响更小且对高阻异常体分辨效果更佳,因此频率域的海洋可控源电磁法在浅海工程和浅部资源探测方面的应用方面失去了优势。
其中,海洋电磁快速成像对数据早期处理、粗处理和实时监控阶段有比较重要的意义。在频率域中,海洋电磁快速成像通常会处理数据形成MVO、PVO等曲线反应地下结构电阻率异常,但其不够直观,不能快速直接的反映出异常体的位置和大小。在时间域中,电磁法快速成像通常采用叠加抽道后的数据进行快速定区,在地面瞬变电磁,航空/半航空瞬变电磁中应用较多,在海洋瞬变电磁中还没有成熟的算法和软件系统。
因此,急需要提出一种逻辑简单、准确可靠的海洋瞬变电磁数据预处理与快速定区成像方法。
发明内容
针对上述问题,本发明的目的在于提供一种海洋瞬变电磁数据预处理与快速定区成像方法,本发明采用的技术方案如下:
一种海洋瞬变电磁数据预处理与快速定区成像方法,其包括以下步骤:
获取电磁信号叠加前测点的海洋瞬变电磁探测对应的感应电动势数据,并按道集提取分为数组第一数据集D;
采用均值圆滑滤波或样条插值滤波或高斯滤波对任一组第一数据集D进行滤波处理;
采用积分离散化将滤波处理后的第一数据集D转化成磁场数据;
预设延时道参数;
对延时断面数据进行均匀半空间的正演处理,并对延时断面数据进行归一化处理;
与现有技术相比,本发明具有以下有益效果:
(1)本发明巧妙地对感应电动势数据按道集提取拆分,以便于滤波和后续处理。另
外,本发明巧妙地使用均值圆滑滤波或样条插值滤波或高斯滤波对第一数据集D、电磁信号
叠加后的数据行滤波处理,其能够快速的去除高频噪声,衰减坏点的影响;
(2)本发明采用积分离散化将滤波处理后的第一数据集D转化成磁场数据,其目的在于,实现高精度的磁场转换过程,磁场数据更利于正反演计算与验证,保证数据处理的准确性;
(3)本发明巧妙地对延时断面数据进行均匀半空间的正演处理,并对延时断面数据进行归一化处理,其采用成像处理和正演分析,提高了数据的可视化,并能对地下典型结构进行快速有效的粗略定性分区,使实际资料处理更直观,更高效;
(4)本发明通过均匀半空间的正演数值的延时断面数据和预处理后的延时断面数据进行归一化计算,瞬变电磁正演数值模拟方法采用先频率域计算,再利用频时转换到时间域,在正演过程中的积分运算采用了汉克尔变换的滤波系数法以提高计算精度;
综上所述,本发明具有逻辑简单、准确可靠等优点,在电磁数据处理技术领域具有很高的实用价值和推广价值。
附图说明
为了更清楚地说明本发明实施例的技术方案,下面将对实施例中所需使用的附图作简单介绍,应当理解,以下附图仅示出了本发明的某些实施例,因此不应被看作是对保护范围的限定,对于本领域技术人员来说,在不付出创造性劳动的前提下,还可以根据这些附图获得其他相关的附图。
图1为本发明的逻辑流程图。
图2为本发明的数据预处理延时断面效果图。
图3为本发明的海洋瞬变电磁数据预处理阶段的测点滤波效果图。
图4为图3的局部放大图。
图5为本发明的海洋瞬变电磁数据预处理阶段的测道滤波前数据分布图。
图6为本发明的海洋瞬变电磁数据预处理阶段的测道滤波后数据分布图。
图7为本发明的积分离散化处理的第一效果图。
图8为本发明的积分离散化处理的第二效果图。
图9为本发明的正演数值模拟结果与延时成像结果展示图。
图10为本发明的快速定区成像效果图。
具体实施方式
为使本申请的目的、技术方案和优点更为清楚,下面结合附图和实施例对本发明作进一步说明,本发明的实施方式包括但不限于下列实施例。基于本申请中的实施例,本领域普通技术人员在没有做出创造性劳动的前提下所获得的所有其他实施例,都属于本申请保护的范围。
本本实施例中,术语“和/或”,仅仅是一种描述关联对象的关联关系,表示可以存在三种关系,例如,A和/或B,可以表示:单独存在A,同时存在A和B,单独存在B这三种情况。
本实施例的说明书和权利要求书中的术语“第一”和“第二”等是用于区别不同的对象,而不是用于描述对象的特定顺序。例如,第一目标对象和第二目标对象等是用于区别不同的目标对象,而不是用于描述目标对象的特定顺序。
在本申请实施例中,“示例性的”或者“例如”等词用于表示作例子、例证或说明。本申请实施例中被描述为“示例性的”或者“例如”的任何实施例或设计方案不应被解释为比其它实施例或设计方案更优选或更具优势。确切而言,使用“示例性的”或者“例如”等词旨在以具体方式呈现相关概念。
在本申请实施例的描述中,除非另有说明,“多个”的含义是指两个或两个以上。例如,多个处理单元是指两个或两个以上的处理单元;多个系统是指两个或两个以上的系统。
如图1至图10所示,本实施例提供了一种海洋瞬变电磁数据预处理与快速定区成像方法,其包括以下步骤:
第一步,获取电磁信号叠加前测点的海洋瞬变电磁探测对应的感应电动势数据,将原始二进制的感应电动势数据按道集提取分为数组第一数据集D,其表达式为:
在本实施例中,共计共M行6列,M为采样点数与采集站的乘积,6列数据分别为时间道数、接收机号和4个分量的数据。
第二步,本实施例采用均值圆滑滤波或样条插值滤波或高斯滤波对任一组第一数据集D进行滤波处理。其中,均值圆滑滤波直接利用前后数据直接求取均值。在样条插值滤波中,分别利用前后数据求取样条插值滤波的公式内的参数,再进行滤波处理。
第三步,采用积分离散化将滤波处理后的第一数据集D转化成磁场数据。具体来说,本实施例的积分离散化的表达式为:
第六步,快速定区成像算法需要对实际工况进行数值模拟,通过均匀半空间的正演数值的延时断面数据和预处理后的延时断面数据进行归一化计算,瞬变电磁正演数值模拟方法采用先频率域计算,再利用频时转换到时间域,在正演过程中的积分运算采用了汉克尔变换的滤波系数法以提高计算精度,其中,利用等对数间隔的方式计算延时时间道列表,其生成公式为:
第七步,采用归一化数据取对数,计算电阻率的变化率,再进行边缘检测,对归一化处理后的数据进行定区预测,并获得A和B两个区域。
在本实施例中,融合了正演数值模拟与快速定区成像,可对预处理之后的数据进行成像处理和正演分析,以提高数据的可视化。本发明能对地下典型结构进行快速有效的粗略定性分区,使实际资料处理更直观、更高效。
上述实施例仅为本发明的优选实施例,并非对本发明保护范围的限制,但凡采用本发明的设计原理,以及在此基础上进行非创造性劳动而作出的变化,均应属于本发明的保护范围之内。
Claims (6)
1.一种海洋瞬变电磁数据预处理与快速定区成像方法,其特征在于,包括以下步骤:
获取电磁信号叠加前测点的海洋瞬变电磁探测对应的感应电动势数据,并按道集提取分为数组第一数据集D;
采用均值圆滑滤波或样条插值滤波或高斯滤波对任一组第一数据集D进行滤波处理;
采用积分离散化将滤波处理后的第一数据集D转化成磁场数据;
预设延时道参数;
根据延时道参数和从滤波处理后的第一数据集D中获得电磁信号叠加后的数据dt;按延时抽道对应测点数据得到延时曲线数组T;
采用均值圆滑滤波或样条插值滤波或高斯滤波对电磁信号叠加后的数据dt行滤波处理,得到延时断面数据;
对延时断面数据进行均匀半空间的正演处理,并对延时断面数据进行归一化处理;
预设映射比例系数k,采用归一化的数据取对数,求得电阻率的变化率,并进行边缘检测和定区预测;采用均匀半空间的正演处理后的数据与延时断面数据通过映射比例系数k并取对数的归一化,得到用于快速成像的数据向量dappr,其表达式为:
dappr=klg(dpr/dfor)
其中,dpr表示预处理的数据向量,即延时断面数据;dfor表示对应时间道正演的数值向量,即采用均匀半空间的正演处理后的数据。
2.根据权利要求1所述的一种海洋瞬变电磁数据预处理与快速定区成像方法,其特征在于,所述第一数据集D的表达式为:
D=[tidx,ridx,ex,hx,hy,hz]M
其中,M表示第一数据集的行数;tidx表示延时道号;ridx表示接收机标号;ex表示x方向的电场数据;hx表示x方向的磁场数据;hy表示y方向的磁场数据,hz表示z方向磁场数据。
4.根据权利要求1所述的一种海洋瞬变电磁数据预处理与快速定区成像方法,其特征在于,所述延时道参数采用梯度间隔逐步增加的时间间隔。
6.根据权利要求1或4所述的一种海洋瞬变电磁数据预处理与快速定区成像方法,其特征在于,所述对延时断面数据进行均匀半空间的正演处理,包括:采用频率域对延时断面数据进行计算,并利用频时转换到时间域;正演处理过程中采用汉克尔变换的滤波系数法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211545726.4A CN115542408B (zh) | 2022-12-05 | 2022-12-05 | 一种海洋瞬变电磁数据预处理与快速定区成像方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211545726.4A CN115542408B (zh) | 2022-12-05 | 2022-12-05 | 一种海洋瞬变电磁数据预处理与快速定区成像方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN115542408A CN115542408A (zh) | 2022-12-30 |
CN115542408B true CN115542408B (zh) | 2023-03-28 |
Family
ID=84722070
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202211545726.4A Active CN115542408B (zh) | 2022-12-05 | 2022-12-05 | 一种海洋瞬变电磁数据预处理与快速定区成像方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115542408B (zh) |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105549101A (zh) * | 2016-01-28 | 2016-05-04 | 中国矿业大学 | 一种瞬变电磁数据微分电导解释方法 |
CN109828307A (zh) * | 2019-03-01 | 2019-05-31 | 山东大学 | 一种瞬变电磁多频率融合的探测方法及应用 |
CN113534270A (zh) * | 2021-07-20 | 2021-10-22 | 中铁二院工程集团有限责任公司 | 一种半航空瞬变电磁电导率-深度成像方法及设备 |
CN114265124A (zh) * | 2021-12-30 | 2022-04-01 | 成都理工大学 | 基于时间域瞬变电磁概率反演的不良地质体定位方法 |
CN115047530A (zh) * | 2022-06-15 | 2022-09-13 | 吉林建筑大学 | 基于一维反演的地空瞬变电磁数据三维频率域解释方法 |
Family Cites Families (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102419453A (zh) * | 2011-07-15 | 2012-04-18 | 中国科学院地质与地球物理研究所 | 长导线源瞬变电磁地空探测方法 |
CN106199742B (zh) * | 2016-06-29 | 2018-02-02 | 吉林大学 | 一种频率域航空电磁法2.5维带地形反演方法 |
US10920585B2 (en) * | 2017-12-26 | 2021-02-16 | Saudi Arabian Oil Company | Determining sand-dune velocity variations |
CN110717138B (zh) * | 2019-09-12 | 2023-06-30 | 吉林大学 | 一种拖曳式瞬变电磁数据预处理方法 |
CN112068212A (zh) * | 2020-07-30 | 2020-12-11 | 成都理工大学 | 一种无人机半航空时间域电磁探测数据分析解释方法 |
CN112114371A (zh) * | 2020-08-21 | 2020-12-22 | 广东电网有限责任公司 | 一种三线圈海缆探测的信号处理系统 |
CN112287545B (zh) * | 2020-10-29 | 2022-07-05 | 吉林大学 | 一种双相导电介质的时-空分数阶电导率建模及模拟方法 |
CN114460654B (zh) * | 2022-02-22 | 2022-10-14 | 成都理工大学 | 基于l1l2混合范数的半航空瞬变电磁数据反演方法及装置 |
CN114721059B (zh) * | 2022-03-22 | 2023-01-31 | 扬州蓝德森科技有限公司 | 一种基于瞬变电磁法的地下成像方法 |
-
2022
- 2022-12-05 CN CN202211545726.4A patent/CN115542408B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105549101A (zh) * | 2016-01-28 | 2016-05-04 | 中国矿业大学 | 一种瞬变电磁数据微分电导解释方法 |
CN109828307A (zh) * | 2019-03-01 | 2019-05-31 | 山东大学 | 一种瞬变电磁多频率融合的探测方法及应用 |
CN113534270A (zh) * | 2021-07-20 | 2021-10-22 | 中铁二院工程集团有限责任公司 | 一种半航空瞬变电磁电导率-深度成像方法及设备 |
CN114265124A (zh) * | 2021-12-30 | 2022-04-01 | 成都理工大学 | 基于时间域瞬变电磁概率反演的不良地质体定位方法 |
CN115047530A (zh) * | 2022-06-15 | 2022-09-13 | 吉林建筑大学 | 基于一维反演的地空瞬变电磁数据三维频率域解释方法 |
Also Published As
Publication number | Publication date |
---|---|
CN115542408A (zh) | 2022-12-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109035224B (zh) | 一种基于多波束点云的海底管道检测与三维重建方法 | |
CN104503432B (zh) | 一种基于小波能量的自主式水下机器人故障辨识方法 | |
CN110824429B (zh) | 深海环境下利用非同步垂直阵的宽带声源被动定位方法 | |
CN108444704B (zh) | 一种滚动轴承早期故障诊断方法 | |
CN107329137B (zh) | 基于box样条小波框架的探地雷达地下目标定位方法及系统 | |
CN109425906B (zh) | 一种磁异常探测矢量磁目标识别方法 | |
CN110542831A (zh) | 一种基于变分模态分解和s变换的故障行波检测方法 | |
CN109085477B (zh) | 用于电力电缆分布式局部放电监测系统的信号识别和定位方法 | |
CN113628227B (zh) | 一种基于深度学习的海岸线变化分析方法 | |
CN104655987A (zh) | 一种基于opgw光偏振态的输电线路雷击点定位新方法 | |
CN112130142A (zh) | 一种复杂运动目标微多普勒特征提取方法及系统 | |
CN111181634B (zh) | 一种分布式光纤振动信号快速定位方法 | |
CN107479037B (zh) | 一种pd雷达杂波区域判别方法 | |
CN107704835B (zh) | 一种利用光谱遥感图像识别海上人工设施的方法 | |
CN102141609A (zh) | 一种高频段一阶Bragg峰及其分裂谱峰的识别方法及系统 | |
CN115542408B (zh) | 一种海洋瞬变电磁数据预处理与快速定区成像方法 | |
CN115015911A (zh) | 一种基于雷达图像的导航地图制作和使用方法及系统 | |
CN113834547B (zh) | 一种河流虚拟站水位时序重建方法及系统 | |
Sørensen et al. | TDEM in groundwater mapping—a continuous approach | |
CN114152837A (zh) | 一种多尺度小波变换下的波头识别方法及装置 | |
CN110275150B (zh) | 基于经验模式分解和迭代端点拟合的变加速运动目标相参积累方法 | |
CN116973914A (zh) | 一种基于三维探地雷达的道路隐蔽病害三维重构方法 | |
CN109188504B (zh) | 一种地下管廊探测方法 | |
Di Martino et al. | An approach to detect lofar lines | |
CN111856400B (zh) | 一种水下目标声源定位方法及系统 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |