CN117169880B - 一种伪随机信号发射方法、接收方法及系统 - Google Patents
一种伪随机信号发射方法、接收方法及系统 Download PDFInfo
- Publication number
- CN117169880B CN117169880B CN202311450878.0A CN202311450878A CN117169880B CN 117169880 B CN117169880 B CN 117169880B CN 202311450878 A CN202311450878 A CN 202311450878A CN 117169880 B CN117169880 B CN 117169880B
- Authority
- CN
- China
- Prior art keywords
- frequency
- pseudo
- electromagnetic wave
- random signal
- low
- 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
- 238000000034 method Methods 0.000 title claims abstract description 37
- 238000005070 sampling Methods 0.000 claims description 37
- 230000008054 signal transmission Effects 0.000 claims description 6
- 230000007274 generation of a signal involved in cell-cell signaling Effects 0.000 claims description 3
- 230000005540 biological transmission Effects 0.000 abstract description 22
- 238000001514 detection method Methods 0.000 abstract description 18
- 238000010276 construction Methods 0.000 abstract description 8
- 238000006243 chemical reaction Methods 0.000 abstract description 6
- 238000005096 rolling process Methods 0.000 abstract description 4
- 230000000694 effects Effects 0.000 description 6
- 238000010586 diagram Methods 0.000 description 5
- 238000005516 engineering process Methods 0.000 description 5
- 238000001228 spectrum Methods 0.000 description 4
- 230000004927 fusion Effects 0.000 description 3
- 238000011161 development Methods 0.000 description 2
- 230000001939 inductive effect Effects 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 230000002035 prolonged effect Effects 0.000 description 2
- 241001522296 Erithacus rubecula Species 0.000 description 1
- 238000009933 burial Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000005672 electromagnetic field Effects 0.000 description 1
- 229910052500 inorganic mineral Inorganic materials 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 239000011707 mineral Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000000737 periodic effect Effects 0.000 description 1
- 230000035699 permeability Effects 0.000 description 1
- 230000010363 phase shift Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 239000011435 rock Substances 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明涉及一种伪随机信号发射方法、接收方法及系统,属于电磁探测技术领域,包括:根据待探测深度范围确定电磁波的最高频率和最低频率;根据待探测目标的径深比确定纵向分辨率;由所述最高频率、所述最低频率和所述纵向分辨率确定电磁波频率组;将每个电磁波频率组均生成一个伪随机信号;按照所述电磁波频率组的基频从高到低或低到高发射所有所述伪随机信号。本发明将每个电磁波频率组的多个频率电磁波融合生成相应的伪随机信号进行发射,并且采用变频方式,按照基频从高到低或低到高,轮巡地发射所有伪随机信号,极大地提高了勘探施工效率,并且可实现分布式滚动式勘探。
Description
技术领域
本发明涉及电磁探测技术领域,特别涉及一种伪随机信号发射方法、接收方法及系统。
背景技术
早期电磁法勘探采用人工换频方式,一个频率电磁波一个频率电磁波地发射,随着GPS同步技术和控制技术的发展,实现了自动换频率发射,后续发展了谐波勘探技术,但其奇次谐波在对数上不等间距,且随着谐波次数的增加,信号衰减的越快。为了解决方波信号的谐波不均匀、快速衰减的问题,提出了2n伪随机信号方案,最典型的是伪随机7频率波信号,一次可发射7个频率电磁波,且在对数上等间隔、幅值接近,但是其频带宽度依然有限,在施工时需要人工变换频率组以满足业务需求,施工效率低。
发明内容
为了解决现有技术的缺陷,本发明提供了一种伪随机信号发射方法、接收方法及系统。
为了解决以上技术问题,本发明采用以下技术方案:
第一方面,本发明实施例提供一种伪随机信号发射方法,包括如下步骤:
根据待探测深度范围确定电磁波的最高频率和最低频率;
根据待探测目标的径深比确定纵向分辨率;
由所述最高频率、所述最低频率和所述纵向分辨率确定电磁波频率组;
将每个电磁波频率组均生成一个伪随机信号;
按照所述电磁波频率组的基频从高到低或低到高发射所有所述伪随机信号。
在一些实施例中,所述由所述最高频率、所述最低频率和所述纵向分辨率确定电磁波频率组包括:
根据所述最低频率确定频率阶数;
由所述最高频率、所述最低频率和所述纵向分辨率确定电磁波序列;
当所述频率阶数大于等于预设阶数时,将每个电磁波序列均分成高频电磁波序列和低频电磁波序列;
将每个所述电磁波序列均生成一个电磁波频率组。
在一些实施例中,所述按照所述电磁波频率组的基频从高到低或低到高发射所有所述伪随机信号包括:
计算当前时间距离预设时间点的整秒数,将所述整秒数除以发射周期获得余数,然后定位该余数在一个发射周期内的波形位置,再从该波形位置按照所述电磁波频率组的基频从高到低或低到高发射所有所述伪随机信号。
第二方面,本发明实施例提供一种伪随机信号接收方法,包括如下步骤:
接收伪随机信号,并计算当前时间距离预设时间点的整秒数,将所述整秒数除以接收周期获得余数,然后定位该余数在一个接收周期内的波形位置,再判断该波形位置的伪随机信号的基频类型;
当伪随机信号的基频是高频时,使用高采样率进行伪随机信号数据采集;
当伪随机信号的基频是低频时,使用低采样率进行伪随机信号数据采集。
在一些实施例中,所述当伪随机信号的基频是高频时,使用高采样率进行伪随机信号数据采集包括:
当伪随机信号的基频是高频时,使用高采样率进行伪随机信号数据采集,将抗混叠滤波器调整到高的拐角频率;
所述当伪随机信号的基频是低频时,使用低采样率进行伪随机信号数据采集包括:
当伪随机信号的基频是低频时,使用低采样率进行伪随机信号数据采集,将抗混叠滤波器调整到低的拐角频率。
第三方面,本发明实施例提供一种伪随机信号发射系统,包括:
最高低频率确定模块,用于根据待探测深度范围确定电磁波的最高频率和最低频率;
纵向分辨率确定模块,用于根据待探测目标的径深比确定纵向分辨率;
电磁波频率组确定模块,用于由所述最高频率、所述最低频率和所述纵向分辨率确定电磁波频率组;
伪随机信号生成模块,用于将每个电磁波频率组均生成一个伪随机信号;
伪随机信号发射模块,用于按照所述电磁波频率组的基频从高到低或低到高发射所有所述伪随机信号。
在一些实施例中,所述由所述最高频率、所述最低频率和所述纵向分辨率确定电磁波频率组包括:
根据所述最低频率确定频率阶数;
由所述最高频率、所述最低频率和所述纵向分辨率确定电磁波序列;
当所述频率阶数大于等于预设阶数时,将每个电磁波序列均分成高频电磁波序列和低频电磁波序列;
将每个所述电磁波序列均生成一个电磁波频率组。
在一些实施例中,所述按照所述电磁波频率组的基频从高到低或低到高发射所有所述伪随机信号包括:
计算当前时间距离预设时间点的整秒数,将所述整秒数除以发射周期获得余数,然后定位该余数在一个发射周期内的波形位置,再从该波形位置按照所述电磁波频率组的基频从高到低或低到高发射所有所述伪随机信号。
第四方面,本发明实施例提供一种伪随机信号接收系统,包括:
伪随机信号接收模块,用于接收伪随机信号,并计算当前时间距离预设时间点的整秒数,将所述整秒数除以接收周期获得余数,然后定位该余数在一个接收周期内的波形位置,再判断该波形位置的伪随机信号的基频类型;
伪随机信号数据采集模块,用于当伪随机信号的基频是高频时,使用高采样率进行伪随机信号数据采集;
当伪随机信号的基频是低频时,使用低采样率进行伪随机信号数据采集。
在一些实施例中,所述当伪随机信号的基频是高频时,使用高采样率进行伪随机信号数据采集包括:
当伪随机信号的基频是高频时,使用高采样率进行伪随机信号数据采集,将抗混叠滤波器调整到高的拐角频率;
所述当伪随机信号的基频是低频时,使用低采样率进行伪随机信号数据采集包括:
当伪随机信号的基频是低频时,使用低采样率进行伪随机信号数据采集,将抗混叠滤波器调整到低的拐角频率。
本发明提供的一种伪随机信号发射方法、接收方法及系统,与现有技术相比,本发明取得的技术效果包括:
1.本发明将每个电磁波频率组的多个频率电磁波融合生成相应的伪随机信号进行发射,并且采用变频方式,按照基频从高到低或低到高,轮巡地发射所有伪随机信号,极大地调高了勘探施工效率,并且可实现分布式滚动式勘探。与此同时每个伪随机信号由7、9、11或13个频率电磁波融合生成,即采用一阶伪随机信号,具有高分频系数特性,取保了信号能量强度,有利于压制干扰获取信噪比更高的数据,提高了勘探精度。
2.用高采样率采集高频伪随机信号,用低采用率采集低频伪随机信号,即不同频率的伪随机信号使用不同的采用率,在保证高质量采集数据的同时,有效减少了采集数据文件量。
3.用高拐角频率采集高频伪随机信号,用低拐角频率采集低频伪随机信号,即不同频率的伪随机信号采集使用不同的拐角频率或采用不同数据链,有效减少了混叠效应影响。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅涉及本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本发明实施例一种伪随机信号发射方法流程图;
图2为步骤S3的方法流程图;
图3为轮巡发射伪随机信号示意图;
图4为本发明实施例一种伪随机信号接收方法流程图;
图5为本发明实施例一种伪随机信号发射系统结构示意图;
图6为本发明实施例一种伪随机信号接收系统结构示意图。
具体实施方式
为了使本领域技术人员更好地理解本发明的方案,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然所描述的实施例仅仅是本发明的一部分实施例,而不是全部的实施例。基于本发明的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的其他方案,都应当属于本发明的保护范围。
为了便于本领域技术人员理解本发明的实施例,以下对本发明涉及的技术术语做出解释。
电磁法勘探是基于地下岩石的导电性和导磁性,通过观测由人工制造的电磁场或天然电磁场诱导的电流分布,可以确定地下岩石的导电性和导磁性,进而推断地下矿藏、油气资源或地层结构等。
伪随机信号,是由周期性数字序列经过调幅、相移等处理后得出的,它具有类似于随机噪声的某些统计特性,同时又能够重复产生。
实施例1
图1示出了本发明实施例一种伪随机信号发射方法流程图,由图1可知,所述伪随机信号发射方法包括如下步骤:
S1、根据待探测深度范围确定电磁波的最高频率和最低频率。
S2、根据待探测目标的径深比确定纵向分辨率。
目前广域电磁法的最大探测深度约为7km,由于电磁法测深部分反演算法采用的是由表及里的剥离式方式,因此浅层信息也很重要,因此探测时,最高频率一般定义为仪器所能发射接收的最高频率,一般为9600Hz,当然电磁波的最高频率和最低频率的确定最终是由实际的探测需求决定。
优选地,首先确定待探测深度范围,然后将该待探测深度范围的两端值代入电磁法勘探探测深度H公式,得到电磁波的最高频率和最低频率。电磁法勘探探测深度H公式为:
;
其中H为探测深度,为勘探区域的背景电阻率,f为勘探电磁波频率。
可通过收集勘探区域以往物探、地质或测井资料了解背景电阻率。通过上述公式,根据待探测深度范围计算出所需电磁波的最高频率和最低频率,为了确保勘探深度,实际作业时的最低频率应比理论计算的最低频率还要低几个频点,优选地,选取最低频率理论值的一半。
实际勘探中,可根据勘探目标体规模确定纵向分辨率,优选地,根据待探测目标的径深比确定纵向分辨率,其中径深比是待探测目标体的纵向厚度与埋藏深度的比值。优选地,由于电磁法勘探的纵向分辨率极限为10%左右,因此对于计算出的小于10%的纵向分辨率按照10%计算。
S3、由所述最高频率、所述最低频率和所述纵向分辨率确定电磁波频率组。
图2示出了步骤S3的方法流程图,由图2可知,该步骤S3包括如下子步骤:
S31、根据所述最低频率确定频率阶数。
将所述最低频率代入如下公式中,计算得到频率阶数。
;
其中n为频率阶数,ceil为向上取整,f为勘探电磁波频率。
S32、由所述最高频率、所述最低频率和所述纵向分辨率确定电磁波序列。
本实施例电磁波序列参考广域电磁法序列或参考GDP32频率表,由所述最高频率和所述最低频率确定,其中电磁波序列中各个频率值成等比数列或近似等比数列,公比可以是2、3或4等等,公比根据勘探需求设定,优选地公比为2,此时频率密度更大。
例如,当最高频率为8192,最低频率为0.015625时,参考广域电磁法序列或GDP32频率表,确定电磁波序列为如下表1中的序列1。仅仅由所述最高频率和所述最低频率确定电磁波序列依然无法满足所有实际勘探需求,此时纵向分辨率有点低,因此再增加电磁波序列,即增加表1中的序列2、3和4,当然还可以增加更多的序列。从表1中计算出,同时使用4个序列即序列1-4时,理论纵向分辨率约为8%;使用序列1和3时,理论纵向分辨率约为16%,使用序列1时,理论纵向分辨率约为29%,综上,可根据实际勘探需求得到的纵向分辨率选择相应的序列。也就是说,由所述最高频率、所述最低频率和所述纵向分辨率确定电磁波序列。
表1-电磁波序列表
S33、当所述频率阶数大于等于预设阶数时,将每个电磁波序列均分成高频电磁波序列和低频电磁波序列。
S34、将每个所述电磁波序列均生成一个电磁波频率组。
考虑到伪随机信号的频谱分布特性,本实施例预设阶数k,k的取值范围为(10,18),优选地k取值为14,此时生成的伪随机信号具有良好的频谱分布特性。
当所述频率阶数大于等于预设阶数时,将每个电磁波序列均分成高频电磁波序列和低频电磁波序列;当所述频率阶数小于预设阶数时,每个电磁波序列都不需要拆分。
例如,当选取表1中的电磁波序列1-4,预设阶数k为14时,频率阶数大于等于预设阶数,表1中的每个电磁波序列将拆分为10个高频电磁波序列和10个低频电磁波序列,由于伪随机信号构造时,偶数频率电磁波将会出现过零的信号,此时发射端需要假负载来平衡发电机的输出,不便于施工,为不降低频率宽影响探测深度,因此将每个电磁波序列拆分成11个高频电磁波序列和11个低频电磁波序列,如表1-2所示,将电磁波序列1拆分成高频电磁波序列1和低频电磁波序列1;将电磁波序列2拆分成高频电磁波序列2和低频电磁波序列2;将电磁波序列3拆分成高频电磁波序列3和低频电磁波序列3;将电磁波序列4拆分成高频电磁波序列4和低频电磁波序列4;优选地,为了既有足够宽的频率又不出现重叠频率,使发出去的伪随机信号能量更高些,以及减少探测时间,将每个电磁波序列拆分成9个高频电磁波序列和11个低频电磁波序列,如表1-2所示,将电磁波序列1拆分成由高频电磁波序列1的前9个频率构成的高频电磁波序列(8192、4096、2048、1024、512、256、128、64、32)和低频电磁波序列1,以此类推,其他电磁波序列拆分不赘述,然后将每个所述电磁波序列均生成一个对应的电磁波频率组。
表2-高低频电磁波序列表
S4、将每个电磁波频率组均生成一个伪随机信号。
S5、按照所述电磁波频率组的基频从高到低或低到高发射所有所述伪随机信号。
根据每个电磁波频率组构造生成对应的伪随机信号,每个伪随机信号的周期由对应电磁波频率组的最低频率决定,将每个电磁波频率组的最低频率定义为基频,即高频电磁波序列1对应的电磁波频率组的基频为8,低频电磁波序列1对应的电磁波频率组的基频为1/64等等,以此类推。按照所述电磁波频率组的基频从高到低或低到高发射所有所述伪随机信号,以减少感抗的影响,从而减少发射机工作负载因不同频率的伪随机信号而产生的波动,进而保护发射端系统,延长发射端的使用寿命。
优选地首先计算当前时间距离预设时间点如UTC时间2000/01/01 00:00:00的整秒数,将所述整秒数除以发射周期获得余数,然后定位该余数在一个发射周期内的波形位置,再从该波形位置起进行码流输出控制发射机进行发射伪随机信号工作,按照所述电磁波频率组的基频从高到低或低到高发射所有所述伪随机信号。
图3示出了伪随机信号发射顺序,从图3可知,按照所述电磁波频率组的基频从高到低发射所有伪随机信号,同理也可以从低到高发射所有伪随机信号(图中未示出),并且可根据勘探需求,设定在一个发射周期内,每个基频的伪随机信号发射次数,不同伪随机信号发射次数可相同,也可不同,其中发射周期等于每个伪随机信号的周期与相应的发射次数相乘再求和。例如,在一个发射周期内发射顺序可以是基频为32、32、28、24、20、20、20、1/16、7/128、3/64、5/128、5/128、5/128、5/128的系列伪随机信号。
本实施例将一个电磁波频率组的多个频率电磁波融合生成一个伪随机信号进行发射,并且采用变频方式,按照基频从高到低或低到高轮巡地发射所有伪随机信号,极大地提高了勘探施工效率,并且可实现分布式滚动式勘探。与此同时每个伪随机信号由7、9、11或13个频率电磁波融合生成,即采用低阶伪随机信号,具有高分频系数特性,取保了信号能量强度,有利于压制干扰获取信噪比更高的数据,提高了勘探精度。
实施例2
图4示出了本发明实施例一种伪随机信号接收方法流程图,由图4可知,所述伪随机信号接收方法包括如下步骤:
S21、接收伪随机信号,并计算当前时间距离预设时间点的整秒数,将所述整秒数除以接收周期获得余数,然后定位该余数在一个接收周期内的波形位置,再判断该波形位置的伪随机信号的基频类型。
首先计算当前时间距离预设时间点如UTC时间2000/01/01 00:00:00的整秒数,将所述整秒数除以接收周期获得余数,然后定位该余数在一个接收周期内的波形位置,再判断该波形位置的伪随机信号的基频类型,其中接收周期等于每个伪随机信号的周期与相应的接收次数相乘再求和,也就是说,接收周期与发射周期相等。
S22、当伪随机信号的基频是高频时,使用高采样率进行伪随机信号数据采集。
S23、当伪随机信号的基频是低频时,使用低采样率进行伪随机信号数据采集。
当接收到伪随机信号时,判断每个伪随机信号的基频类型,如果伪随机信号的基频是高频时,然后采用高采样率进行伪随机信号数据采集,同时通过接收端的电子开关,如mos开关,将抗混叠滤波器调整到高的拐角频率;
如果伪随机信号的基频是低频时,然后采用低采样率进行伪随机信号数据采集,同时通过接收端的电子开关,如mos开关,将抗混叠滤波器调整到低的拐角频率。当然,也可以采用两个独立的信号链路,高采样率和低采样率分别走不同的数据链,以减少混叠效应影响。
本实施例用高采样率采集高频伪随机信号,用低采用率采集低频伪随机信号,即不同频率的伪随机信号使用不同的采用率,在保证高质量采集数据的同时,有效减少了采集数据文件量。
本实施例用高拐角频率采集高频伪随机信号,用低拐角频率采集低频伪随机信号,即不同频率的伪随机信号采集使用不同的拐角频率或采用不同数据链,有效减少了混叠效应影响。
实施例3
图5示出了本发明实施例一种伪随机信号发射系统结构示意图,由图5可知,所述伪随机信号发射系统包括:最高低频率确定模块51、纵向分辨率确定模块52、电磁波频率组确定模块53、伪随机信号生成模块54和伪随机信号发射模块55;
最高低频率确定模块51,用于根据待探测深度范围确定电磁波的最高频率和最低频率。
纵向分辨率确定模块52,用于根据待探测目标的径深比确定纵向分辨率。
目前广域电磁法的最大探测深度约为7km,由于电磁法测深部分反演算法采用的是由表及里的剥离式方式,因此浅层信息也很重要,因此探测时,最高频率一般定义为仪器所能发射接收的最高频率,一般为9600Hz,当然电磁波的最高频率和最低频率的确定最终是由实际的探测需求决定。
优选地,首先确定待探测深度范围,然后将该待探测深度范围的两端值代入电磁法勘探探测深度H公式,得到电磁波的最高频率和最低频率。电磁法勘探探测深度H公式为:
;
其中H为探测深度,为勘探区域的背景电阻率,f为勘探电磁波频率。
可通过收集勘探区域以往物探、地质或测井资料了解背景电阻率。通过上述公式,根据待探测深度范围计算出所需电磁波的最高频率和最低频率,为了确保勘探深度,实际作业时的最低频率应比理论计算的最低频率还要低几个频点,优选地,选取最低频率理论值的一半。
电磁波频率组确定模块53,用于由所述最高频率、所述最低频率和所述纵向分辨率确定电磁波频率组。
根据所述最低频率确定频率阶数。
将所述最低频率代入如下公式中,计算得到频率阶数。
;
其中n为频率阶数,ceil为向上取整,f为勘探电磁波频率。
由所述最高频率、所述最低频率和所述纵向分辨率确定电磁波序列。
本实施例电磁波序列参考广域电磁法序列或参考GDP32频率表,由所述最高频率和所述最低频率确定,其中电磁波序列中各个频率值成等比数列或近似等比数列,公比可以是2、3或4等等,公比根据勘探需求设定,优选地公比为2,此时频率密度更大。
例如,当最高频率为8192,最低频率为0.015625时,参考广域电磁法序列或GDP32频率表,确定电磁波序列为如下表1中的序列1。仅仅由所述最高频率和所述最低频率确定电磁波序列依然无法满足所有实际勘探需求,有时纵向分辨率有点低,因此再增加电磁波序列,即增加表1中的序列2、3和4,当然还可以增加更多的序列。从表1中计算出,同时使用4个序列即序列1-4时,理论纵向分辨率约为8%;使用序列1和3时,理论纵向分辨率约为16%,使用序列1时,理论纵向分辨率约为29%,综上,可根据实际勘探需求得到的纵向分辨率选择相应的序列。也就是说,由所述最高频率、所述最低频率和所述纵向分辨率确定电磁波序列。
当所述频率阶数大于等于预设阶数时,将每个电磁波序列均分成高频电磁波序列和低频电磁波序列。
将每个所述电磁波序列均生成一个电磁波频率组。
考虑到伪随机信号的频谱分布特性,本实施例预设阶数k,k的取值范围为(10,18),优选地k取值为14,此时生成的伪随机信号具有良好的频谱分布特性。
当所述频率阶数大于等于预设阶数时,将每个电磁波序列均分成高频电磁波序列和低频电磁波序列;当所述频率阶数小于预设阶数时,每个电磁波序列都不需要拆分。
例如,当选取表1中的电磁波序列1-4,预设阶数k为14时,频率阶数大于预设阶数,表1中的每个电磁波序列将拆分为10个高频电磁波序列和10个低频电磁波序列,由于伪随机信号构造时,偶数频率电磁波将会出现过零的信号,此时发射端需要假负载来平衡发电机的输出,不便于施工,为不降低频率宽影响探测深度,因此将每个电磁波序列拆分成11个高频电磁波序列和11个低频电磁波序列,如表1-2所示,将电磁波序列1拆分成高频电磁波序列1和低频电磁波序列1;将电磁波序列2拆分成高频电磁波序列2和低频电磁波序列2;将电磁波序列3拆分成高频电磁波序列3和低频电磁波序列3;将电磁波序列4拆分成高频电磁波序列4和低频电磁波序列4;优选地,为了既有足够宽的频率又不出现重叠频率,使发出去的伪随机信号能量更高些,以及减少探测时间,将每个电磁波序列拆分成9个高频电磁波序列和11个低频电磁波序列,如表1-2所示,将电磁波序列1拆分成由高频电磁波序列1的前9个频率构成的高频电磁波序列(8192、4096、2048、1024、512、256、128、64、32)和低频电磁波序列1,以此类推,其他电磁波序列拆分不赘述,然后将每个所述电磁波序列均生成一个对应的电磁波频率组。
伪随机信号生成模块54,用于将每个电磁波频率组均生成一个伪随机信号。
伪随机信号发射模块55,用于按照所述电磁波频率组的基频从高到低或低到高发射所有所述伪随机信号。
根据每个电磁波频率组构造生成对应的伪随机信号,每个伪随机信号的周期由对应电磁波频率组的最低频率决定,将每个电磁波频率组的最低频率定义为基频,即高频电磁波序列1对应的电磁波频率组的基频为8,低频电磁波序列1对应的电磁波频率组的基频为1/64等等,以此类推。按照所述电磁波频率组的基频从高到低或低到高发射所有所述伪随机信号,以减少感抗的影响,从而减少发射机工作负载因不同频率的伪随机信号而产生的波动,进而保护发射端系统,延长发射端的使用寿命。
优选地首先计算当前时间距离预设时间点如UTC时间2000/01/01 00:00:00的整秒数,将所述整秒数除以发射周期获得余数,然后定位该余数在一个发射周期内的波形位置,再从该波形位置起进行码流输出控制发射机进行发射伪随机信号工作,按照所述电磁波频率组的基频从高到低或低到高发射所有所述伪随机信号。
图3示出了伪随机信号发射顺序,从图3可知,按照所述电磁波频率组的基频从高到低发射所有伪随机信号,同理也可以从低到高发射所有伪随机信号(图中未示出),并且可根据勘探需求,设定在一个发射周期内,每个基频的伪随机信号发射次数,不同伪随机信号发射次数可相同,也可不同,其中发射周期等于每个伪随机信号的周期与相应的发射次数相乘再求和。例如,在一个发射周期内发射顺序可以是基频为32、32、28、24、20、20、20、1/16、7/128、3/64、5/128、5/128、5/128、5/128的系列伪随机信号。
本实施例将一个电磁波频率组的多个频率电磁波融合生成一个伪随机信号进行发射,并且采用变频方式,按照基频从高到低或低到高轮巡地发射所有伪随机信号,极大地调高了勘探施工效率,并且可实现分布式滚动式勘探。与此同时每个伪随机信号由7、9、11或13个频率电磁波融合生成,即采用低阶伪随机信号,具有高分频系数特性,取保了信号能量强度,有利于压制干扰获取信噪比更高的数据,提高了勘探精度。
实施例4
图6示出了本发明实施例一种伪随机信号接收系统结构示意图,由图6可知,所述伪随机信号接收系统包括:伪随机信号接收模块61和伪随机信号数据采集模块62;
伪随机信号接收模块61,用于接收伪随机信号,并计算当前时间距离预设时间点的整秒数,将所述整秒数除以接收周期获得余数,然后定位该余数在一个接收周期内的波形位置,再判断该波形位置的伪随机信号的基频类型。
首先计算当前时间距离预设时间点如UTC时间2000/01/01 00:00:00的整秒数,将所述整秒数除以接收周期获得余数,然后定位该余数在一个接收周期内的波形位置,再判断该波形位置的伪随机信号的基频类型,其中接收周期等于每个伪随机信号的周期与相应的接收次数相乘再求和,也就是说,接收周期与发射周期相等。
伪随机信号数据采集模块62,用于当伪随机信号的基频是高频时,使用高采样率进行伪随机信号数据采集;当伪随机信号的基频是低频时,使用低采样率进行伪随机信号数据采集。
当接收到伪随机信号时,判断每个伪随机信号的基频类型,如果伪随机信号的基频是高频时,然后采用高采样率进行伪随机信号数据采集,同时通过接收端的电子开关,如mos开关,将抗混叠滤波器调整到高的拐角频率;
如果伪随机信号的基频是低频时,然后采用低采样率进行伪随机信号数据采集,同时通过接收端的电子开关,如mos开关,将抗混叠滤波器调整到低的拐角频率。当然,也可以采用两个独立的信号链路,高采样率和低采样率分别走不同的数据链,以减少混叠效应影响。
本实施例用高采样率采集高频伪随机信号,用低采用率采集低频伪随机信号,即不同频率的伪随机信号使用不同的采用率,在保证高质量采集数据的同时,有效减少了采集数据文件量。
本实施例用高拐角频率采集高频伪随机信号,用低拐角频率采集低频伪随机信号,即不同频率的伪随机信号采集使用不同的拐角频率或采用不同数据链,有效减少了混叠效应影响。
上述实施例为本发明较佳的实施方式,但本发明的实施方式并不受上述实施例的限制,其他的任何未背离本发明的精神实质与原理下所作的改变、修饰、替代、组合、简化,均应为等效的置换方式,都包含在本发明的保护范围之内。
Claims (6)
1.一种伪随机信号发射方法,其特征在于,包括如下步骤:
根据待探测深度范围确定电磁波的最高频率和最低频率;
根据待探测目标的径深比确定纵向分辨率;
由所述最高频率、所述最低频率和所述纵向分辨率确定电磁波频率组;
将每个电磁波频率组均生成一个伪随机信号;
按照所述电磁波频率组的基频从高到低或低到高发射所有所述伪随机信号;
所述由所述最高频率、所述最低频率和所述纵向分辨率确定电磁波频率组包括:
根据所述最低频率确定频率阶数;
由所述最高频率、所述最低频率和所述纵向分辨率确定电磁波序列;
当所述频率阶数大于等于预设阶数时,将每个电磁波序列均分成高频电磁波序列和低频电磁波序列;
将每个所述电磁波序列均生成一个电磁波频率组。
2.根据权利要求1所述的伪随机信号发射方法,其特征在于,所述按照所述电磁波频率组的基频从高到低或低到高发射所有所述伪随机信号包括:
计算当前时间距离预设时间点的整秒数,将所述整秒数除以发射周期获得余数,然后定位该余数在一个发射周期内的波形位置,再从该波形位置按照所述电磁波频率组的基频从高到低或低到高发射所有所述伪随机信号。
3.一种伪随机信号接收方法,其特征在于,包括如下步骤:
接收权利要求1或2所述的伪随机信号发射方法发射的伪随机信号,并计算当前时间距离预设时间点的整秒数,将所述整秒数除以接收周期获得余数,然后定位该余数在一个接收周期内的波形位置,再判断该波形位置的伪随机信号的基频类型;
当伪随机信号的基频是高频时,使用高采样率进行伪随机信号数据采集;
当伪随机信号的基频是低频时,使用低采样率进行伪随机信号数据采集;
所述当伪随机信号的基频是高频时,使用高采样率进行伪随机信号数据采集包括:
当伪随机信号的基频是高频时,使用高采样率进行伪随机信号数据采集,将抗混叠滤波器调整到高的拐角频率;
所述当伪随机信号的基频是低频时,使用低采样率进行伪随机信号数据采集包括:
当伪随机信号的基频是低频时,使用低采样率进行伪随机信号数据采集,将抗混叠滤波器调整到低的拐角频率。
4.一种伪随机信号发射系统,其特征在于,包括:
最高低频率确定模块,用于根据待探测深度范围确定电磁波的最高频率和最低频率;
纵向分辨率确定模块,用于根据待探测目标的径深比确定纵向分辨率;
电磁波频率组确定模块,用于由所述最高频率、所述最低频率和所述纵向分辨率确定电磁波频率组;
伪随机信号生成模块,用于将每个电磁波频率组均生成一个伪随机信号;
伪随机信号发射模块,用于按照所述电磁波频率组的基频从高到低或低到高发射所有所述伪随机信号;
所述由所述最高频率、所述最低频率和所述纵向分辨率确定电磁波频率组包括:
根据所述最低频率确定频率阶数;
由所述最高频率、所述最低频率和所述纵向分辨率确定电磁波序列;
当所述频率阶数大于等于预设阶数时,将每个电磁波序列均分成高频电磁波序列和低频电磁波序列;
将每个所述电磁波序列均生成一个电磁波频率组。
5.根据权利要求4所述的伪随机信号发射系统,其特征在于,所述按照所述电磁波频率组的基频从高到低或低到高发射所有所述伪随机信号包括:
计算当前时间距离预设时间点的整秒数,将所述整秒数除以发射周期获得余数,然后定位该余数在一个发射周期内的波形位置,再从该波形位置按照所述电磁波频率组的基频从高到低或低到高发射所有所述伪随机信号。
6.一种伪随机信号接收系统,其特征在于,包括:
伪随机信号接收模块,用于接收权利要求4或5所述的伪随机信号发射系统发射的伪随机信号,并计算当前时间距离预设时间点的整秒数,将所述整秒数除以接收周期获得余数,然后定位该余数在一个接收周期内的波形位置,再判断该波形位置的伪随机信号的基频类型;
伪随机信号数据采集模块,用于当伪随机信号的基频是高频时,使用高采样率进行伪随机信号数据采集;
当伪随机信号的基频是低频时,使用低采样率进行伪随机信号数据采集;
所述当伪随机信号的基频是高频时,使用高采样率进行伪随机信号数据采集包括:
当伪随机信号的基频是高频时,使用高采样率进行伪随机信号数据采集,将抗混叠滤波器调整到高的拐角频率;
所述当伪随机信号的基频是低频时,使用低采样率进行伪随机信号数据采集包括:
当伪随机信号的基频是低频时,使用低采样率进行伪随机信号数据采集,将抗混叠滤波器调整到低的拐角频率。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202311450878.0A CN117169880B (zh) | 2023-11-03 | 2023-11-03 | 一种伪随机信号发射方法、接收方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202311450878.0A CN117169880B (zh) | 2023-11-03 | 2023-11-03 | 一种伪随机信号发射方法、接收方法及系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN117169880A CN117169880A (zh) | 2023-12-05 |
CN117169880B true CN117169880B (zh) | 2024-02-09 |
Family
ID=88947284
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202311450878.0A Active CN117169880B (zh) | 2023-11-03 | 2023-11-03 | 一种伪随机信号发射方法、接收方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN117169880B (zh) |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1896769A (zh) * | 2006-05-17 | 2007-01-17 | 中国地质大学(北京) | 海底平面波电磁场探测装置及测量方法 |
JP2014029663A (ja) * | 2012-07-05 | 2014-02-13 | Sumitomo Electric Ind Ltd | 監視システムおよび監視方法 |
CN105607131A (zh) * | 2016-01-11 | 2016-05-25 | 甘肃省有色地质调查院 | 编码源电磁测深法获取大地全区视电阻率谱的方法及装置 |
CN108427145A (zh) * | 2018-01-26 | 2018-08-21 | 吉林大学 | 地空频率域电磁法可控频率源探测信号脉宽调制方法 |
CN111505722A (zh) * | 2020-04-27 | 2020-08-07 | 山东大学 | 一种对数非均匀伪随机电磁勘探信号生成方法及系统 |
CN111796330A (zh) * | 2020-07-13 | 2020-10-20 | 中国科学院地质与地球物理研究所 | 一种时频联合探测波的合成方法和装置以及探测方法 |
CN114646820A (zh) * | 2020-12-18 | 2022-06-21 | 施耐德电气美国股份有限公司 | 改进与检测异常状况相关联的问题的识别的系统和方法 |
CN115236748A (zh) * | 2022-06-30 | 2022-10-25 | 深圳市吉奥地球科技有限公司 | 一种频率域组合波归一化表达方法以及装置、存储介质 |
CN116203639A (zh) * | 2023-03-01 | 2023-06-02 | 吉林大学 | 一种拖曳式频域电磁探测深度聚焦发射系统及方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
AU2009262987B2 (en) * | 2008-06-24 | 2013-09-12 | Exxonmobil Upstream Research Company | Method for electromagnetic prospecting waveform design |
-
2023
- 2023-11-03 CN CN202311450878.0A patent/CN117169880B/zh active Active
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1896769A (zh) * | 2006-05-17 | 2007-01-17 | 中国地质大学(北京) | 海底平面波电磁场探测装置及测量方法 |
JP2014029663A (ja) * | 2012-07-05 | 2014-02-13 | Sumitomo Electric Ind Ltd | 監視システムおよび監視方法 |
CN105607131A (zh) * | 2016-01-11 | 2016-05-25 | 甘肃省有色地质调查院 | 编码源电磁测深法获取大地全区视电阻率谱的方法及装置 |
CN108427145A (zh) * | 2018-01-26 | 2018-08-21 | 吉林大学 | 地空频率域电磁法可控频率源探测信号脉宽调制方法 |
CN111505722A (zh) * | 2020-04-27 | 2020-08-07 | 山东大学 | 一种对数非均匀伪随机电磁勘探信号生成方法及系统 |
CN111796330A (zh) * | 2020-07-13 | 2020-10-20 | 中国科学院地质与地球物理研究所 | 一种时频联合探测波的合成方法和装置以及探测方法 |
CN114646820A (zh) * | 2020-12-18 | 2022-06-21 | 施耐德电气美国股份有限公司 | 改进与检测异常状况相关联的问题的识别的系统和方法 |
CN115236748A (zh) * | 2022-06-30 | 2022-10-25 | 深圳市吉奥地球科技有限公司 | 一种频率域组合波归一化表达方法以及装置、存储介质 |
CN116203639A (zh) * | 2023-03-01 | 2023-06-02 | 吉林大学 | 一种拖曳式频域电磁探测深度聚焦发射系统及方法 |
Also Published As
Publication number | Publication date |
---|---|
CN117169880A (zh) | 2023-12-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Anderson et al. | An integrated approach to marine electromagnetic surveying using a towed streamer and source | |
CN102183792B (zh) | 人工源频率域电磁测深装置及测量方法 | |
CN102721982B (zh) | 基于扩频编码技术的地面电磁勘探方法及其探测系统 | |
CN102590869B (zh) | 一种人工场源频率域电法勘探方法及勘探系统 | |
Trichtchenko et al. | Modeling geomagnetically induced currents using geomagnetic indices and data | |
US20220308249A1 (en) | Method and system for generating logarithmic non-uniform pseudo-random electromagnetic exploration signal | |
CN111239839B (zh) | 一种针对频率域电磁法的频谱密度计算方法及装置 | |
CN110187394B (zh) | 双场源电磁测深法获取地层电阻率各向异性的方法及装置 | |
CN103955003A (zh) | 一种在超导瞬变电磁应用中的噪声抑制方法 | |
CN102053278A (zh) | 电法勘探方法及测量装置 | |
CN110879418B (zh) | 一种频率域电性源电磁场视电导率计算方法 | |
Samanes et al. | Estimating the VLF modal interference distance using the South America VLF Network (SAVNET) | |
CN110031901A (zh) | 一种变深电法勘测系统和一种勘测方法 | |
Möllhoff et al. | Seismic noise characterization in proximity to strong microseism sources in the Northeast Atlantic | |
Chand et al. | VLF modal interference distance and nighttime D region VLF reflection height for west-east and east-west propagation paths to Fiji | |
Maxwell et al. | Natural noise fields from 1 cps to 100 kc | |
CN108761540B (zh) | 一种频率域天然电场三维勘探方法 | |
CN117169880B (zh) | 一种伪随机信号发射方法、接收方法及系统 | |
CN108594301B (zh) | 一种具有差异特征的地震数据融合的方法及处理终端 | |
CN104793237A (zh) | 一种获得宽频可控震源扫描信号的方法和装置 | |
Kolobov et al. | The KVVN-7 multifunction digital measuring station for electromagnetic monitoring of seismoactive zones | |
CN115755200A (zh) | 一种井下直流电法监测的多频数据抗干扰方法 | |
Hornbostel et al. | Waveform design for electroseismic exploration | |
Allen | The influence of attenuation in earthquake ground-motion and magnitude estimation: implications for Australian earthquake hazard | |
Bourdillon et al. | Effects of geomagnetic pulsations on the Doppler shift of HF backscatter radar echoes |
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 |