CN106772629B - 一种基于Gold编码的多发多收电磁探测方法 - Google Patents
一种基于Gold编码的多发多收电磁探测方法 Download PDFInfo
- Publication number
- CN106772629B CN106772629B CN201710026992.9A CN201710026992A CN106772629B CN 106772629 B CN106772629 B CN 106772629B CN 201710026992 A CN201710026992 A CN 201710026992A CN 106772629 B CN106772629 B CN 106772629B
- Authority
- CN
- China
- Prior art keywords
- gold
- gold code
- survey line
- coding
- impulse response
- 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
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/08—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation operating with magnetic or electric fields produced or modified by objects or geological structures or by detecting devices
- G01V3/083—Controlled source electromagnetic [CSEM] surveying
Landscapes
- Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Life Sciences & Earth Sciences (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
本发明公开了一种基于Gold编码的多发多收电磁探测方法,属于地球物理勘探技术领域。目的是提供一种能够实现多条测线的同步探测,实现提高电磁探测的工作效率的基于Gold编码的多发多收电磁探测方法,方法具体包括:步骤一生成Gold码;步骤二将Gold码分配给各条测线作为激励源;步骤三大地脉冲响应的辨识。本发明通过在不同测线使用独立的Gold编码作为激励源,利用Gold码优良的互相关特性与较佳的正交性,降低不同测线之间信号源的相互影响,利用Gold码尖锐的自相关特性,实现大地脉冲响应的辨识。本发明所述***能够实现多条测线的同步勘察,大大提高了电磁探测的工作效率。
Description
技术领域
本发明具体涉及一种基于Gold编码的多发多收电磁探测方法,属于地球物理勘探技术领域。
背景技术
在地球物理勘探领域中,传统的电磁探测常采用双极性方波作为激励源,波形简单、频谱单一、易受干扰,若采用多测线同步勘察时不同激励源的波场信号会在接收区混叠,而在后期处理中无法辨别不同激励源对应的响应信号,因此在进行矿区勘察时往往以单发多收的方式以测线为单位进行探测,工时长、效率低。
随着编码电流探测的兴起,以m序列作为激励源的电磁探测***逐渐成为了研究热点,并已得到初步应用。虽然m序列具备尖锐的自相关特性,能够提高大地脉冲响应曲线的辨识精度,但不同编码的m序列互相关函数波动较大,即若使用不同编码的m序列进行多发多收探测,会导致与发射方波情况下相同的问题。
发明内容
因此,为解决上述问题,本发明提供一种能够实现多条测线的同步探测,实现提高电磁探测的工作效率的基于Gold编码的多发多收电磁探测方法,
本发明的方法具体包括:
步骤一生成Gold码;
步骤二将Gold码分配给各条测线作为激励源;
步骤三大地脉冲响应的辨识。
进一步的,根据需求以及工区情况确定测线数与测线位置,由于高阶Gold码具备更低的互相关特性,构造第一m序列发生器,通过更改本原多项式以匹配公式构造第二m序列发生器,在同一时钟脉冲的控制下,产生的两个m序列通过模二加法后得到Gold码,通过配置第二m序列发生器的不同初始状态,得到不同的Gold序列。
其中,所述匹配公式为:
6.进一步的,所述步骤二具体为:
使用Gold码对发射电流进行编码,产生双极性发射波形,若发射电流为I,则编码映射关系为:1→I、0→-I;在进行发射时,以周期为单位进行数轮发射,每条测线的收发相互独立,互不干扰。
2.进一步的,所述步骤三具体为:
将大地等效为LTI***,其脉冲响应为g1(t),经Gold码编码的发射信号为i1(t),接收信号r(t)表示为:r(t)=i1(t)*g1(t),若考虑其它测线发射源的影响,r(t)表示为:r(t)=i1(t)*g1(t)+i2(t)*g2(t)+…+in(t)*gn(t),
其中i2(t)…in(t)为其它测线的Gold码激励源信号,g1(t)…gn(t)为其它测线到本测线间的大地脉冲响应,等号两边对i 1(t)求相关,可写作:
R(r,i1)≈g1(t)*R(i1,i1),
R(r,i1)≈g1(t)*R(i1,i1),简写为a(t)=g(t)*b(t),对该式进行频域转换,可得:A(w)=G(w)·B(w),时间域中的卷积变为频域乘积,解得G(w):
对上式变换到时间域即可解得g(t),即完成大地脉冲响应的辨识。
本发明的有益效果在于:本发明的基于Gold编码的多发多收电磁探测方法,通过在不同测线使用独立的Gold编码作为激励源,利用Gold码优良的互相关特性与较佳的正交性,降低不同测线之间信号源的相互影响,利用Gold码尖锐的自相关特性,实现大地脉冲响应的辨识。本发明所述***能够实现多条测线的同步勘察,大大提高了电磁探测的工作效率。
附图说明
图1为本发明基于Gold编码的多发多收电磁探测的收发布置图。
图2为本发明Gold码产生示意图。
具体实施方式
下面结合附图对本发明的具体实施方式进行说明:
本发明的技术方案的基本原理为:Gold码有着较低的互相关值的同时兼具类似于m序列的尖锐自相关特性。能够在降低不同发射源相互影响的同时,实现大地脉冲响应曲线的辨识,能够实现多发多收的多测线同步勘察,提高了电测探测的工作效率。
多发多收电磁探测***的收发布置如图1所示。
Gold码原理为:Gold码是两个m序列的模二和,作为产生Gold码的两个m序列,它们的互相关是有界的,需要满足以下关系:其中n为m序列的阶数:
因此,Gold序列并非m序列,它只是m序列的组合码;但它保留了m序列的优良特性,其零延迟相关值R(0)远大于延迟相关值R(k)(k不等于0),即具有尖锐的自相关特性;同时其互相关函数相较于m序列而言要小得多,其最大的互相关值不会超过同阶数m序列间的最大互相关值;另外Gold码具有远多于同阶m序列的独立码组。
由于Gold的低互相关特性,不同的Gold激励源在接收区带来的混叠影响在后期的数据处理中可通过相关辨识进行压制,因此可实现在一块测区进行多测线同步探测。首先根据需求以及工区情况确定测线数与测线位置,由于高阶Gold码具备更低的互相关特性,因此若为提高工作效率而选择较低的测线间隔时,应选择不同的高阶Gold码作为各测线激励源。
基于Gold编码的多发多收电磁探测方法的具体实施步骤为:
首先构造第一m序列发生器,通过更改本原多项式以匹配公式1构造第二m序列发生器,在同一时钟脉冲的控制下,产生的两个m序列通过模二加后得到Gold码,通过配置第二m序列发生器的不同初始状态,可以得到不同的Gold序列,如图2所示。
产生不同的Gold码后,即可分配给各条测线作为激励源,具体为使用Gold码对发射电流进行编码,产生双极性发射波形,若发射电流为I,则编码映射关系为:1→I、0→-I;在进行发射时,为获取完备的收发数据,往往以周期为单位进行数轮发射,而不同Gold序列经过循环移位后,其互相关特性不变,因此每条测线无需统一的发射授时,即每条测线的收发相互独立,互不干扰。
各条测线的大地脉冲响应相关辨识算法具体为:将大地等效为LTI(线性时不变)***,其脉冲响应为g1(t),本测线经Gold码编码的发射信号为i1(t),若不考虑噪声,且仅考虑本测线发射源,则接收信号r(t)可表示为:r(t)=i1(t)*g1(t),若考虑其它测线发射源的影响,r(t)可表示为:r(t)=i1(t)*g1(t)+i2(t)*g2(t)+…+in(t)*gn(t).
上式中i2(t)…in(t)为其它测线的Gold码激励源信号,g1(t)…gn(t)为其它测线到本测线间的大地脉冲响应,等号两边对i1(t)求相关,可写作:
R(r,i1)=g1(t)*R(i1,i1)+g2(t)*R(i1,i2)+…+gn(t)*R(i1,in)
由于不同测线使用的Gold码激励源之间具备较低的互相关值,即R(i1,i2)…R(i1,in)均近似为零,因此上式可写作R(r,i1)≈g1(t)*R(i1,i1),即,在其它测线均有激励源的情况下,本测线的数据处理不受其它测线激励源的影响。
考虑R(r,i1)≈g1(t)*R(i1,i1),简写为a(t)=g(t)*b(t),对该式进行频域转换,可得:A(w)=G(w)·B(w),时间域中的卷积变为频域乘积,因此可解得G(w):
对上式变换到时间域即可解得g(t),即完成大地脉冲响应的辨识。
以上所述是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明所述原理的前提下,还可以作出若干改进和润饰,这些改进和润饰也应视为本发明的保护范围。
Claims (1)
1.一种基于Gold编码的多发多收电磁探测方法,其特征在于,所述方法包括以下步骤:
步骤一生成Gold码;
步骤二将Gold码分配给各条测线作为激励源;
步骤三大地脉冲响应的辨识;
所述步骤一具体为:
根据需求以及工区情况确定测线数与测线位置,由于高阶Gold码具备更低的互相关特性,构造第一m序列发生器,通过更改本原多项式以匹配公式构造第二m序列发生器,在同一时钟脉冲的控制下,产生的两个m序列通过模二加法后得到Gold码,通过配置第二m序列发生器的不同初始状态,得到不同的Gold序列;
所述匹配公式为:
所述步骤二具体为:
使用Gold码对发射电流进行编码,产生双极性发射波形,若发射电流为I,则编码映射关系为:1→I、0→-I;在进行发射时,以周期为单位进行数轮发射,每条测线的收发相互独立,互不干扰;
所述步骤三具体为:
将大地等效为LTI***,其脉冲响应为g1(t),经Gold码编码的发射信号为i1(t),接收信号r(t)表示为:r(t)=i1(t)*g1(t),若考虑其它测线发射源的影响,r(t)表示为:r(t)=i1(t)*g1(t)+i2(t)*g2(t)+…+in(t)*gn(t),
其中i2(t)…in(t)为其它测线的Gold码激励源信号,g2(t)…gn(t)为其它测线到本测线间的大地脉冲响应,等号两边对i1(t)求相关,可写作:
R(r,i1)≈g1(t)*R(i1,i1),
R(r,i1)≈g1(t)*R(i1,i1),简写为a(t)=g(t)*b(t),对该式进行频域转换,可得:A(w)=G(w)·B(w),时间域中的卷积变为频域乘积,解得G(w):
对上式变换到时间域即可解得g(t),即完成大地脉冲响应的辨识。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710026992.9A CN106772629B (zh) | 2017-01-15 | 2017-01-15 | 一种基于Gold编码的多发多收电磁探测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710026992.9A CN106772629B (zh) | 2017-01-15 | 2017-01-15 | 一种基于Gold编码的多发多收电磁探测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106772629A CN106772629A (zh) | 2017-05-31 |
CN106772629B true CN106772629B (zh) | 2019-03-12 |
Family
ID=58946679
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710026992.9A Active CN106772629B (zh) | 2017-01-15 | 2017-01-15 | 一种基于Gold编码的多发多收电磁探测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106772629B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107589461B (zh) * | 2017-09-04 | 2019-11-05 | 长江大学 | 一种基于双向编码的尖脉冲时域电磁深层探测方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102508292A (zh) * | 2011-10-27 | 2012-06-20 | 吉林大学 | 可控震源匹配扫描方法 |
CN102721982A (zh) * | 2012-06-26 | 2012-10-10 | 中国科学院电工研究所 | 基于扩频编码技术的地面电磁勘探方法及其探测*** |
CN104035128A (zh) * | 2013-03-06 | 2014-09-10 | 中国石油集团东方地球物理勘探有限责任公司 | 可控震源伪随机扫描信号生成方法 |
CN104217099A (zh) * | 2014-07-18 | 2014-12-17 | 中国石油化工股份有限公司 | 一种用于产生并处理可控震源的扫描信号的方法 |
CN104502980A (zh) * | 2014-12-08 | 2015-04-08 | 中国科学院电子学研究所 | 一种电磁大地冲激响应的辨识方法 |
CN105938205A (zh) * | 2016-03-04 | 2016-09-14 | 中国科学院地质与地球物理研究所 | 一种多道瞬变电磁法接收波形记录合成方法与装置 |
-
2017
- 2017-01-15 CN CN201710026992.9A patent/CN106772629B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102508292A (zh) * | 2011-10-27 | 2012-06-20 | 吉林大学 | 可控震源匹配扫描方法 |
CN102721982A (zh) * | 2012-06-26 | 2012-10-10 | 中国科学院电工研究所 | 基于扩频编码技术的地面电磁勘探方法及其探测*** |
CN104035128A (zh) * | 2013-03-06 | 2014-09-10 | 中国石油集团东方地球物理勘探有限责任公司 | 可控震源伪随机扫描信号生成方法 |
CN104217099A (zh) * | 2014-07-18 | 2014-12-17 | 中国石油化工股份有限公司 | 一种用于产生并处理可控震源的扫描信号的方法 |
CN104502980A (zh) * | 2014-12-08 | 2015-04-08 | 中国科学院电子学研究所 | 一种电磁大地冲激响应的辨识方法 |
CN105938205A (zh) * | 2016-03-04 | 2016-09-14 | 中国科学院地质与地球物理研究所 | 一种多道瞬变电磁法接收波形记录合成方法与装置 |
Non-Patent Citations (2)
Title |
---|
伪随机编码源激发下的时域电磁信号合成;王若 等;《地球物理学报》;20161215;第59卷(第12期);4414-4423 * |
基于m伪随机序列的电磁法抗噪能力分析;王显祥 等;《地球物理学报》;20160531;第59卷(第5期);1861-1874 * |
Also Published As
Publication number | Publication date |
---|---|
CN106772629A (zh) | 2017-05-31 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP2705901B2 (ja) | カオスコードを利用するレーダシステム | |
RU2251715C2 (ru) | Сигналы источника для электросейсмической разведки | |
AU2017248407B2 (en) | Coded signals for marine vibrators | |
CN104181509A (zh) | 基于跳频和多相交替码的非相干散射雷达信号处理方法 | |
WO1992001957A1 (fr) | Procede et appareil pour stratigraphie du sous-sol a l'aide d'un radar | |
CN112285786B (zh) | 一种基于超高密度电法设备提高勘探深度的方法 | |
CN106772629B (zh) | 一种基于Gold编码的多发多收电磁探测方法 | |
Sharma et al. | Four-phase orthogonal code design for MIMO radar systems | |
JP5454475B2 (ja) | 位置検出システム、送信装置、受信装置、位置検出方法、位置検出プログラム | |
Homier et al. | Hybrid fixed-dwell-time search techniques for rapid acquisition of ultra-wideband signals | |
US8098712B2 (en) | Optical correlation apparatus and method | |
CN104459802B (zh) | 一种用于电法勘探***中的发射装置 | |
EP1072092B1 (en) | Monopulse generator | |
CN106549683B (zh) | 一种扩频编码探测信号发射*** | |
Sharma et al. | MIMO radar ambiguity analysis of frequency hopping pulse waveforms | |
García et al. | Multilevel complementary sets of sequences and their application in UWB | |
CN107733547A (zh) | 一种电磁探测的扩频编码信号及其发生*** | |
RU7211U1 (ru) | Устройство измерения задержки сигнала | |
CA2280386A1 (en) | Electromagnetic geophysical mapping instruments | |
RU2012019C1 (ru) | Способ акустического каротажа | |
Wong | Spread spectrum techniques for seismic data acquisition | |
Cai et al. | Orthogonal waveform set design for netted radar | |
RU2012018C1 (ru) | Способ акустического каротажа | |
RU2642430C1 (ru) | Способ измерения дальности | |
Donato et al. | Uncorrelated binary sequences of lengths 2a3b4c5d7e11f13g based on nested Barker codes and complementary sequences |
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 |