CN111522064B - 一种高阶伪随机电磁勘探信号生成方法及系统 - Google Patents

一种高阶伪随机电磁勘探信号生成方法及系统 Download PDF

Info

Publication number
CN111522064B
CN111522064B CN202010344733.2A CN202010344733A CN111522064B CN 111522064 B CN111522064 B CN 111522064B CN 202010344733 A CN202010344733 A CN 202010344733A CN 111522064 B CN111522064 B CN 111522064B
Authority
CN
China
Prior art keywords
signals
basic unit
phase
signal
frequency
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
Application number
CN202010344733.2A
Other languages
English (en)
Other versions
CN111522064A (zh
Inventor
何继善
杨洋
李帝铨
瓮晶波
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Hunan Geosun Hi Technology Co ltd
Shandong University
Central South University
Original Assignee
Hunan Geosun Hi Technology Co ltd
Shandong University
Central South University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Hunan Geosun Hi Technology Co ltd, Shandong University, Central South University filed Critical Hunan Geosun Hi Technology Co ltd
Priority to CN202010344733.2A priority Critical patent/CN111522064B/zh
Publication of CN111522064A publication Critical patent/CN111522064A/zh
Application granted granted Critical
Publication of CN111522064B publication Critical patent/CN111522064B/zh
Priority to US17/601,972 priority patent/US11983511B2/en
Priority to PCT/CN2021/083727 priority patent/WO2021218524A1/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

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
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F7/00—Methods or arrangements for processing data by operating upon the order or content of the data handled
    • G06F7/58—Random or pseudo-random number generators
    • G06F7/582—Pseudo-random number generators
    • 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/12—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation operating with electromagnetic waves
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F7/00—Methods or arrangements for processing data by operating upon the order or content of the data handled
    • G06F7/60—Methods or arrangements for performing computations using a digital non-denominational number representation, i.e. number representation without radix; Computing devices using combinations of denominational and non-denominational quantity representations, e.g. using difunction pulse trains, STEELE computers, phase computers
    • G06F7/68—Methods or arrangements for performing computations using a digital non-denominational number representation, i.e. number representation without radix; Computing devices using combinations of denominational and non-denominational quantity representations, e.g. using difunction pulse trains, STEELE computers, phase computers using pulse rate multipliers or dividers pulse rate multipliers or dividers per se

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Theoretical Computer Science (AREA)
  • Remote Sensing (AREA)
  • Environmental & Geological Engineering (AREA)
  • Geology (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Geophysics (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • General Engineering & Computer Science (AREA)
  • Electromagnetism (AREA)
  • Computing Systems (AREA)
  • Mathematical Physics (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

本发明公开了一种高阶伪随机电磁勘探信号生成方法及系统，所述方法包括：根据勘探需求构造两个或两个以上基本单元信号，其中，所述基本单元信号是多个同相位周期方波信号叠加得到的阶梯信号，所述多个方波信号相邻之间频率之比为2；将所述两个或两个以上基本单元信号叠加，得到叠加阶梯信号，并修正幅值使其与方波信号幅值一致，得到高阶2n序列伪随机信号。本发明能够在有限频率区间内构造不同阶数的2n序列阶梯信号，改善在电磁勘探时的分辨率。

Description

一种高阶伪随机电磁勘探信号生成方法及系统
技术领域
本发明涉及地球物理勘探领域，尤其涉及一种高阶伪随机电磁勘探信号生成方法及系统。
背景技术
本部分的陈述仅仅是提供了与本公开相关的背景技术信息，不必然构成在先技术。
目前，有很多种电磁勘探类方法，大多采用不同伪随机编码方式进行信号编码，尤其是在频率域勘探领域。但是过去很多伪随机信号的有效频率频带范围往往很窄，在勘探过程中需要不断改变勘探信号的频率成分，以获得不同频率区间的勘探信息，再依次发送，大大浪费时间，且提高了勘探成本。
发明内容
为克服上述现有技术的不足，本发明提供一种高阶伪随机电磁勘探信号生成方法及系统，能够在有限频率区间内构造不同阶数的2n序列伪随机信号，相比传统2n序列伪随机信号具有更高的分辨率。
为实现上述目的，本发明的一个或多个实施例提供了如下技术方案：
一种高阶伪随机电磁勘探信号生成方法，包括以下步骤：
根据勘探需求构造两个或两个以上基本单元信号，其中，所述基本单元信号是多个同相位周期方波信号叠加得到的阶梯信号，所述多个方波信号相邻之间频率之比为2；
将所述两个或两个以上基本单元信号叠加，得到叠加阶梯信号，并修正幅值使其与方波信号幅值一致，得到高阶2n序列伪随机信号。
一个或多个实施例提供了一种终端设备，包括处理器和计算机可读存储介质，处理器用于实现各指令；计算机可读存储介质用于存储多条指令，其特征在于，所述指令适于由处理器加载并执行所述的高阶伪随机电磁勘探信号生成方法。
一个或多个实施例提供了一种计算机可读存储介质，其中存储有多条指令，其特征在于，所述指令适于由终端设备的处理器加载并执行所述的高阶伪随机电磁勘探信号生成方法。
一个或多个实施例提供了一种电法勘探信号发送机，采用所述的高阶伪随机电磁勘探信号生成方法生成高阶2n序列伪随机信号。
一个或多个实施例提供了一种电法勘探信号发送机，采用所述的高阶伪随机电磁勘探信号生成方法生成高阶2n序列伪随机信号。
以上一个或多个技术方案存在以下有益效果：
采用以上技术方案，能够根据勘探需求，在有限频段范围内增加主频个数，在有限频率区间内构造不同阶数的2n序列阶梯信号，以提高主频之间的密度，有助于改善在电磁勘探时的分辨率。
附图说明
构成本发明的一部分的说明书附图用来提供对本发明的进一步理解，本发明的示意性实施例及其说明用于解释本发明，并不构成对本发明的不当限定。
图1为本发明实施例中高阶伪随机电磁勘探信号生成方法流程图；
图2为本发明实施例中1Hz周期正弦信号与1Hz周期方波信号(相位φ＝0)对应关系示意图；
图3为本发明实施例中7个频率叠加后台阶信号示意图；
图4为本发明实施例中7个频率叠加后台阶信号(a)及其频谱示意图(b)；
图5为本发明实施例中6个频率叠加后台阶信号(a)及其频谱示意图(b)；
图6为本发明实施例中包含13频波叠加信号及对应生成的伪随机信号。
图7为本发明实施例中包含13频波伪随机信号时域波形(a)及频谱示意图(b)；
图8为本发明实施例中第二个构造单元选取不同相位情况下的主频频谱相对均方误差曲线；
图9为本发明实施例中3Hz-6频-2π/5相位构造单元对应时域波形(a)及频谱特征(b)；
图10为本发明实施例中1Hz-7频-0相位、3Hz-6频-2π/5相位时对应伪随机信号时域波形(a)及频谱特征(b)；
图11为本发明实施例中第二个构造单元不同相位情况下的主频频谱相对均方误差曲线；
图12为本发明实施例中1Hz-11频-0相位与3Hz-10频-67π/180相位构造单元生成伪随机信号L2-F21—1Hz-1536Hz的时域波形(a)及频谱特征(b)；
图13为本发明实施例中3阶信号构造单元不同相位变化时主频频谱幅值相对均方误差变化曲线；
图14为本发明实施例中3阶2n序列信号(L3-F29—1Hz-1536Hz)波形(a)及频谱(b)；
图15为本发明实施例中4阶信号相位变化曲线；
图16为本发明实施例中4阶2n序列信号(L3-F29—1Hz-1536Hz)波形(a)及频谱(b)；
图17为本发明实施例中5阶信号随相位变化主频频谱相对均方误差变化曲线；
图18为本发明实施例中5阶2n序列信号(L5-F45—1Hz-1536Hz)波形(a)及频谱(b)；
图19为本发明实施例中2阶信号随相位变化主频频谱相对均方误差变化曲线；
图20为本发明实施例中1Hz-7频-π/9相位、3Hz-6频-π/9相位时对应伪随机信号时域波形(a)及频谱(b)；
图21为本发明实施例中适用于浅层勘探的2阶信号(L2-F27-0.25Hz-3072Hz)时域波形(a)及频谱(b)；
图22为本发明实施例中适用于浅层勘探的3阶信号(L3-F39-0.25Hz-3072Hz)时域波形(a)及频谱(b)。
以上附图涉及的频谱图中，仅折线的拐点存在频谱值，为更清楚的展示频谱幅值的大小和变化趋势，本发明实施例中均采用折线图表示。
具体实施方式
应该指出，以下详细说明都是示例性的，旨在对本发明提供进一步的说明。除非另有指明，本文使用的所有技术和科学术语具有与本发明所属技术领域的普通技术人员通常理解的相同含义。
需要注意的是，这里所使用的术语仅是为了描述具体实施方式，而非意图限制根据本发明的示例性实施方式。如在这里所使用的，除非上下文另外明确指出，否则单数形式也意图包括复数形式，此外，还应当理解的是，当在本说明书中使用术语“包含”和/或“包括”时，其指明存在特征、步骤、操作、器件、组件和/或它们的组合。
在不冲突的情况下，本发明中的实施例及实施例中的特征可以相互组合。
为了演示，实施例中将以计算机仿真信号为例进行说用，所采用采样频率远高于伪随机信号中最高主频频率，计算机引入后对信号离散化所引起计算误差可以忽略。
实施例一
本实施例公开了一种高阶伪随机电磁勘探信号生成方法，如图1所示，所述方法包括：
步骤1：根据勘探需求构造两个或两个以上基本单元信号，其中，所述基本单元信号是多个同相位周期方波信号叠加得到的阶梯信号，所述多个方波信号相邻之间频率之比为2。
其中，所述周期方波的生成方法为：
生成与目标周期方波频率一致的正弦信号，并将这些正弦信号中，大于0的值置为A，小于0的值置为-A，等于0的位置按照索引奇偶分别置A和-A，具体地，索引识别信号中0的所有位置，按照奇偶项设置等于A和-A，如第1、3、5…个0位置数值置为A，第2、4、6…个0位置数值置为-A，得到一系列相位相同的周期方波信号，其中，A≠0，相位值可根据需求调整，本实施例中相位为0。
其中，周期方波的主频频率范围与勘探深度相关，主频频率的最低与最高频率可以根据待探测的深度需求(例如深层、浅层或其他具体探测深度)，根据经验公式辅助选择。具体地，获取用户需求的待探测深度范围；根据待探测深度范围确定最高频率和最低频率。然后结合获取的需求的频率增加倍数，得到目标周期方波相应的频率。
以1Hz为例说明，首先生成一个幅值为100、相位为0的1Hz正弦信号：S(t)＝A sin(2πft+φ)，其中，S(t)为正弦信号，A为幅值，φ为相位，f为信号频率，幅值为100时，A＝100相位为0时φ＝0，频率为1Hz时，f＝1。
对这些正弦信号中大于0的值置100，小于0的值置-100，等于0的位置按照位置索引奇偶分别置100和-100，在这样的方式下，则可以生成1Hz的0相位方波信号，如图2所示，依此方式则可以生成所有频率幅值为100、相位为0的频率。
设R(t)为周期方波信号，当正弦信号在t位置S(t)为正值时，将信号数值置为A，即对应时间位置t位置R(t)＝A；当正弦信号在t位置S(t)为负值时，将信号数值置为-A，即对应时间位置t位置R(t)＝-A，本实施例中将正弦信号的相位都设为0，取A＝100，对这些正弦信号中大于0的值置100，小于0的值置-100，等于0的位置按照位置索引奇偶分别置100和-100。
本实施例中，为表述方便称周期方波最低频率为基频，图3生成信号的基频为1Hz。注意，在此叠加信号需为2n序列频率，即相邻方波频率之间需以连续2的倍数增加，依据此方法能够生成以1Hz为基频的其他频率如2Hz、4Hz、8Hz…64Hz相位为0时的周期方波信号，而通过叠加1Hz、2Hz、4Hz、8Hz、16Hz、32Hz、64Hz的0相位周期方波后，则可以获得含有多个频率的台阶信号，叠加后波形及对应频谱如图4所示。在频率域叠加后信号，1Hz、2Hz、4Hz、8Hz、16Hz、32Hz、64Hz对应频谱值是相同的，为100。这是因为每个方波只存在奇次谐波，不存在偶次谐波，所以1Hz的奇次谐波成分没有影响2Hz及其谐波频谱成分，同理2Hz及其奇次谐波成分并未对4Hz及其奇次谐波成分造成影响，同理依次类推，所以所有主频及其对应谐波的幅值是相同的，具有非常好的特征，此台阶信号作为第一个基本单元信号，该基本单元含有7个频率，最低频率为1Hz，相位φ＝0，则称之为1Hz-7频-0相位2n序列基本单元。
设最低频率为3Hz,基本单元中的频率数为6，其构造信号及频谱图为图5所示，6个周期方波频率为3Hz、6Hz、12Hz、24Hz、48Hz、96Hz，相位同样为0，该叠加后的阶梯信号作为第二个基本单元信号，称为3Hz-6频-0相位2n序列基本单元。
基本单元的最低频率可以改变，频率数不一定为奇数，可以为偶数。本实施例中，第一个基本单元的最低频率称之为基频，第二个基本单元的最低频率为基频的3倍，第三个基本单元的最低频率为基频的5倍，第四个基本单元的最低频率为基频的7倍，第n个基本单元的最低频率为基频的2n+1倍。
本实施例中，要求对不同基本单元中主频数目进行求和后为奇数，不应为偶数，数学上，能够证明对奇数个周期方波(信号中数值只存在100与-100)相互叠加时，其幅值不可能为0，只可能是100或者-100的整数倍，所以经过叠加后，只需将大于100位置置为数值100，小于-100位置置为数值-100。假设第一个基本单元中主频数目为a1、第二个基本单元中主频数目a2、第三个基本单元中主频个数a3，Numf为所有基本单元主频数目的和，Numf＝a1+a2+a3。不难发现，当构造n阶2n序列伪随机信号时，Numf＝a1+a2+a3+…+an，其中an为第n个基本单元中的主频数目。本发明中，在构造高阶2n序列伪随机信号时，Numf须为奇数，且对每个基本单元中的数目是否为奇数并不要求。同时，一般情况下，a1数目应大于其他基本单元的主频数，以此来约束高阶2n序列伪随机信号中最高频率。
以上构造要求简化了后续叠加信号的处理难度，能够高效的得到勘探所需高阶2n序列伪随机信号。
本实施例中，基本单元中频率数目是可以改变的，可以将基本单元1Hz-11频-0相位和基本单元3Hz-10频-67π/180相位叠加构造2阶2n序列伪随机信号，图11和图12所示，其中复杂2n序列伪随机信号的编码方式为L2-F21--1Hz-1536Hz，其中L为Level的缩写，意为阶数，F为Frequency的缩写，意为信号包含主频个数，1Hz表示信号中最低频率主频，1536Hz表示信号中的最高频率主频。不同阶数代表主频之间的疏密程度，2阶比1阶密，3阶要比2阶密，以此类推。
本实施例根据勘探需求，在有限频段范围内增加主频个数，在有限频率区间内构造不同阶数的2n序列阶梯信号，以提高主频之间的密度，有助于改善在电磁勘探时的分辨率。
步骤2：将所述两个或两个以上基本单元信号叠加，得到叠加阶梯信号，并修正幅值使其与方波信号幅值一致，得到高阶2n序列伪随机信号。
具体地，将以上两组基本单元叠加后获得一组新的叠加阶梯信号，并对叠加阶梯信号中幅值大于A的位置幅值置为A，幅值小于-A的位置幅值置为-A，获得2n序列伪随机信号。如图6所示，红色为叠加阶梯信号，蓝色为伪随机信号。该伪随机信号波形及频谱如图7所示。在时间域，其信号中只存在100与-100数值，而在频率域含有13个主频频率，即两个基本单元中7个主频与6个主频的相加数,且由于基本单元中每个方波只有100和-100,数学上可以证明，13个100或-100叠加后一定不等于0，一定是大于等于100，或者小于等于-100的数值，因此通过将叠加后的数值进行大于100置100和小于-100置-100后，其信号中只会存在100和-100，不会存在等于0的位置。
但是，直接相加两个基本单元所得到的伪随机信号，主频频谱能量差异较大，并不适合应用于频率域电磁勘探，因此，还对叠加方式进行寻优。
步骤3：对于所述两个或多个基本单元信号中的至少一个进行相位调整，寻找使得叠加后高阶2n序列伪随机信号中各主频对应频谱幅值的均方差最小的相位，作为相应基本单元信号的最优相位；根据最优相位将所述两个或两个以上基本单元信号叠加，并修正幅值使其与方波信号幅值一致，得到高阶2n序列伪随机信号。
本实施例通过计算主频对应频谱幅值的均方差大小来判断其幅值之间的均匀性，选择使均方差数据最小的相位作为最优相位，以使各主要频率(主频)对应频谱幅值均匀，即主频能量分布更加均匀。具体地，通过不断改变信号构造基本单元的相位，即正弦信号生成公式中φ，进而生成对应相位的周期方波叠加后的阶梯信号，利用这种组合方式生成高阶2n序列伪随机信号，以使各主要频率(主频)对应频谱幅值均匀，即主频能量分布更加均匀。
具体地，对于进行相位调整的基本单元信号，首先将半个周期的相位弧度π拆分成N份(大于半个周期相位弧度π后，信号反向，但是信号整体形态一致，所以只需要考虑半个周期相位即可)，即相位改变基本单元为π/N；每次调整一个相位改变基本单元，记录叠加后高阶2n序列伪随机信号中各主频对应频谱值的均方差；绘制均方差大小随相位的变化曲线，寻找使得均方差最小的相位作为该基本单元信号的最佳相位。
在本例中，假设将伪随机信号生成过程中将半个周期的相位弧度π，即拆分成π/180，每次相位改变π/180，依次增加，即每次改变公式中φ的值，其可能的值依次为0、π/180、2π/180、3π/180...179π/180。第一次构造高阶2n序列伪随机信号时，1Hz-7频-0相位基本单元不变，3Hz-6频-0相位作为第二个基本单元，两个基本单元进行叠加后，进行大于100位置置为数值100和小于-100位置置为数值-100后，计算并记录此时13个主频对应频谱值的相对均方误差值；改变第二个构造单位的相位至π/180，变为3Hz-6频-π/180，计算并记录此时13个主频对应频谱值的相对均方误差值，以此类推获得2π/180、3π/180...179π/180，180个不同相位下对应伪随机信号的相对均方误差曲线，通过比较获得相对均方误差最小值对应第二个基本单元的相位，如图8所示，相对均方误差最小点位于横坐标数值73，即第二个基本单元对应相位为2π/5，以此作为此13个频率下伪随机信号的最优构造方式，即第一个基本单元0相位，即为1Hz-7频-0相位，第二个基本单元相位2π/5时，即3Hz-6频-2π/5相位基本单元，如图9所示，获得L2-F13-1Hz-96Hz伪随机信号的最优构造方式，如图10所示。
第一基本单元对应相位φ允许改变。如针对1Hz-7频基本单元和3Hz-6频基本单元进行叠加合成2n序列信号时，两个基本单元相位同时改变，当N＝36时对应相对均方误差变化曲线如图19所示，其最小值点两个基本单元对应相位分别为π/9与2π/9，此时对应2阶2n序列伪随机信号时域波形和频谱如图20所示。
若要通过三个基本单元构造3阶2n序列伪随机信号时，在构造此高阶2n序列伪随机信号时，将半个周期的相位弧度π拆分成60个基本单元(即N＝60)，固定第一个构造单位相位为0，则另外两个构造基本单元的相位变化组合共有3600种，计算并记录每一种组合情况下的伪随机信号组合方式，比较获得主频频谱幅值最小相对均方误差位置，并获取在此位置对应不同基本单元的相位信息，如图13和图14所示。
通过不断改变构造基本单元的相位获取主频频谱幅值相对均方误差最小值点，即频谱成分最优分配2n序列伪随机信号。
如图14所示，其伪随机信号编号为L3-F29—1Hz-1536Hz，即为3阶29频、最低主频1Hz，最高主频1536Hz的2n序列伪随机信号，其三个基本单元分别为1Hz-11频-0相位、3Hz-10频-π/2相位及5Hz-8频-9π/20相位，三个基本单元在对应相位条件下获得此频段主频能量最优分配的2阶2n序列伪随机信号，如图14b频谱图所示。
本发明所提出构造方法，在固定频率区间内，可以构造更高阶信号，如图15和图16中所示4阶2n序列伪随机信号L4-F37—1Hz-1536Hz，图14和图15所示5阶2n序列伪随机信号L5-F45—1Hz-1536Hz。
本发明中生成高阶2n序列伪随机信号频率区间，是能够根据勘探需要改变的，一般情况下频率区间的最低频率不低于0.015625Hz，最高频率不高于8192Hz。
采用本实施例的方法，能够构造适用于浅层高精度勘探的信号：L2-F27-0.25Hz-3072Hz，其中，主频为0.25Hz构造单元相位为10度，主频为0.75Hz构造单元相位为120度；以及，L3-F39-0.25Hz-3072Hz，其中，主频为0.25Hz构造单元相位为0度，主频为0.75Hz构造单元相位为90度，主频1.25Hz构造单元相位为90度。
实施例二
本实施例的目的是提供一种高阶伪随机电磁勘探信号生成系统，包括：
基本单元信号构造模块，根据勘探需求构造两个或两个以上基本单元信号，其中，所述基本单元信号是多个同相位周期方波信号叠加得到的阶梯信号，所述多个方波信号相邻之间频率之比为2；
高阶伪随机信号生成模块，将所述两个或两个以上基本单元信号叠加，得到叠加阶梯信号，并修正幅值使其与方波信号幅值一致，得到高阶2n序列伪随机信号。
实施例三
本实施例的目的是提供一种计算机可读存储介质，其中存储有多条指令，所述指令适于由终端设备的处理器加载并执行：
根据勘探需求构造两个或两个以上基本单元信号，其中，所述基本单元信号是多个同相位周期方波信号叠加得到的阶梯信号，所述多个方波信号相邻之间频率之比为2；
将所述两个或两个以上基本单元信号叠加，得到叠加阶梯信号，并修正幅值使其与方波信号幅值一致，得到高阶2n序列伪随机信号。
实施例四
本实施例的目的是提供一种终端设备，包括处理器和计算机可读存储介质，处理器用于实现各指令；计算机可读存储介质用于存储多条指令，所述指令适于由处理器加载并执行：
根据勘探需求构造两个或两个以上基本单元信号，其中，所述基本单元信号是多个同相位周期方波信号叠加得到的阶梯信号，所述多个方波信号相邻之间频率之比为2；
将所述两个或两个以上基本单元信号叠加，得到叠加阶梯信号，并修正幅值使其与方波信号幅值一致，得到高阶2n序列伪随机信号。
实施例五
本实施例的目的是提供一种电法勘探信号发送机，采用实施例一种所述信号生成方法，得到符合勘探需求的高阶2n序列伪随机信号。
以上所述仅为本发明的优选实施例而已，并不用于限制本发明，对于本领域的技术人员来说，本发明可以有各种更改和变化。凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。
上述虽然结合附图对本发明的具体实施方式进行了描述，但并非对本发明保护范围的限制，所属领域技术人员应该明白，在本发明的技术方案的基础上，本领域技术人员不需要付出创造性劳动即可做出的各种修改或变形仍在本发明的保护范围以内。

Claims (7)

1.一种高阶伪随机电磁勘探信号生成方法，其特征在于，包括以下步骤：
根据勘探需求构造两个或两个以上基本单元信号，其中，所述基本单元信号是多个同相位周期方波信号叠加得到的阶梯信号，所述多个方波信号相邻之间频率之比为2；
将所述两个或两个以上基本单元信号叠加，得到叠加阶梯信号，并修正幅值使其与方波信号幅值一致，得到高阶2n序列伪随机信号；
所述方法还包括：
对于所述两个或多个基本单元信号中的至少一个进行相位调整，寻找使得叠加后高阶2n序列伪随机信号中各主频对应频谱幅值的均方差最小的相位，作为相应基本单元信号的最优相位；
根据最优相位将所述两个或两个以上基本单元信号叠加，并修正幅值使其与方波信号幅值一致，得到最终的高阶2n序列伪随机信号；
最优相位寻找方法为：
对于进行相位调整的基本单元信号，首先将半个周期的相位弧度π拆分成N份，得到相位改变基本单元π/N；
每次调整一个相位改变基本单元，记录叠加后高阶2n序列伪随机信号中各主频对应频谱值的均方差；
绘制均方差大小随相位的变化曲线，寻找使得均方差最小的相位作为该基本单元信号的最佳相位；
根据所述方法构造得到的适用于浅层勘探的信号如下：
二阶2n序列伪随机信号，其中，两个基本单元信号的主频和相位分别为：0.25Hz和10度，0.75Hz和120度；以及，三阶2n序列伪随机信号，其中，三个基本单元信号的主频和相位分别为：0.25Hz和0度，0.75Hz和90度，1.25Hz和90度。
2.如权利要求1所述的一种高阶伪随机电磁勘探信号生成方法，其特征在于，所述两个或两个以上基本单元信号满足以下要求：
第一个基本单元信号中主频的最低频率记为基频，第n个基本单元的最低频率为基频的2n+1倍，n为大于1的自然数；
所有基本单元信号的主频数目之和为奇数。
3.如权利要求1所述的一种高阶伪随机电磁勘探信号生成方法，其特征在于，周期方波的生成方法为：
首先生成与目标周期方波频率一致的正弦信号，将这些正弦信号中，大于0的值置为A，小于0的值置为-A，得到一系列相位相同的周期方波信号，其中，A≠0；
对于周期方波中值为0的信号，按位置赋值A或-A，得到不含0值的周期方波。
4.一种高阶伪随机电磁勘探信号生成系统，其特征在于，包括：
基本单元信号构造模块，根据勘探需求构造两个或两个以上基本单元信号，其中，所述基本单元信号是多个同相位周期方波信号叠加得到的阶梯信号，所述多个方波信号相邻之间频率之比为2；
高阶伪随机信号生成模块，将所述两个或两个以上基本单元信号叠加，得到叠加阶梯信号，并修正幅值使其与方波信号幅值一致，得到高阶2n序列伪随机信号；
对于所述两个或多个基本单元信号中的至少一个进行相位调整，寻找使得叠加后高阶2n序列伪随机信号中各主频对应频谱幅值的均方差最小的相位，作为相应基本单元信号的最优相位；
根据最优相位将所述两个或两个以上基本单元信号叠加，并修正幅值使其与方波信号幅值一致，得到最终的高阶2n序列伪随机信号；
最优相位寻找方法为：
对于进行相位调整的基本单元信号，首先将半个周期的相位弧度π拆分成N份，得到相位改变基本单元π/N；
每次调整一个相位改变基本单元，记录叠加后高阶2n序列伪随机信号中各主频对应频谱值的均方差；
绘制均方差大小随相位的变化曲线，寻找使得均方差最小的相位作为该基本单元信号的最佳相位；
根据所述方法构造得到的适用于浅层勘探的信号如下：
二阶2n序列伪随机信号，其中，两个基本单元信号的主频和相位分别为：0.25Hz和10度，0.75Hz和120度；以及，三阶2n序列伪随机信号，其中，三个基本单元信号的主频和相位分别为：0.25Hz和0度，0.75Hz和90度，1.25Hz和90度。
5.一种终端设备，包括处理器和计算机可读存储介质，处理器用于实现各指令；计算机可读存储介质用于存储多条指令，其特征在于，所述指令适于由处理器加载并执行权利要求1-3中任一项所述的高阶伪随机电磁勘探信号生成方法。
6.一种计算机可读存储介质，其中存储有多条指令，其特征在于，所述指令适于由终端设备的处理器加载并执行权利要求1-3中任一项所述的高阶伪随机电磁勘探信号生成方法。
7.一种电法勘探信号发送机，其特征在于，采用如权利要求1-3任一项所述的高阶伪随机电磁勘探信号生成方法生成高阶2n序列伪随机信号。
CN202010344733.2A 2020-04-27 2020-04-27 一种高阶伪随机电磁勘探信号生成方法及系统 Active CN111522064B (zh)

Priority Applications (3)

Application Number Priority Date Filing Date Title
CN202010344733.2A CN111522064B (zh) 2020-04-27 2020-04-27 一种高阶伪随机电磁勘探信号生成方法及系统
US17/601,972 US11983511B2 (en) 2020-04-27 2021-03-30 Method and system for generating high-order pseudo-random electromagnetic exploration signal
PCT/CN2021/083727 WO2021218524A1 (zh) 2020-04-27 2021-03-30 一种高阶伪随机电磁勘探信号生成方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010344733.2A CN111522064B (zh) 2020-04-27 2020-04-27 一种高阶伪随机电磁勘探信号生成方法及系统

Publications (2)

Publication Number Publication Date
CN111522064A CN111522064A (zh) 2020-08-11
CN111522064B true CN111522064B (zh) 2021-01-26

Family

ID=71910885

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010344733.2A Active CN111522064B (zh) 2020-04-27 2020-04-27 一种高阶伪随机电磁勘探信号生成方法及系统

Country Status (3)

Country Link
US (1) US11983511B2 (zh)
CN (1) CN111522064B (zh)
WO (1) WO2021218524A1 (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111522064B (zh) * 2020-04-27 2021-01-26 山东大学 一种高阶伪随机电磁勘探信号生成方法及系统
CN115236748B (zh) * 2022-06-30 2023-10-27 深圳市吉奥地球科技有限公司 一种频率域组合波归一化表达方法以及装置、存储介质
CN115794028B (zh) * 2023-02-13 2023-04-25 山东大学 幅值靶向增强的伪随机信号构建方法、系统、介质及设备

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101672928A (zh) * 2009-09-29 2010-03-17 中南大学 伪随机多频仪接收机自动测试系统
CN102098073A (zh) * 2011-01-27 2011-06-15 王红星 一种基于重构滤波的时限带通椭圆球面波函数信号产生方法
CN111505723A (zh) * 2020-04-27 2020-08-07 山东大学 一种2n序列伪随机电磁勘探信号生成方法及系统
CN111505722A (zh) * 2020-04-27 2020-08-07 山东大学 一种对数非均匀伪随机电磁勘探信号生成方法及系统

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1163764C (zh) * 2000-05-19 2004-08-25 何继善 一种主动源频率域电法勘探方法
EP3103474A1 (en) 2015-06-12 2016-12-14 Institut Pasteur Live recombinant measles-m2 virus and its use in eliciting immunity against influenza viruses
JP7051050B2 (ja) 2018-04-02 2022-04-11 株式会社Ihi 物理探査方法及び音源
CN109257102B (zh) * 2018-09-30 2021-06-25 西南交通大学 一种基于光子技术的多阶微波跳频信号发生器
CN109005005B (zh) 2018-11-05 2019-04-05 湖南继善高科技有限公司 一种伪随机信号混合编码方法及系统
CN111522064B (zh) 2020-04-27 2021-01-26 山东大学 一种高阶伪随机电磁勘探信号生成方法及系统

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101672928A (zh) * 2009-09-29 2010-03-17 中南大学 伪随机多频仪接收机自动测试系统
CN102098073A (zh) * 2011-01-27 2011-06-15 王红星 一种基于重构滤波的时限带通椭圆球面波函数信号产生方法
CN111505723A (zh) * 2020-04-27 2020-08-07 山东大学 一种2n序列伪随机电磁勘探信号生成方法及系统
CN111505722A (zh) * 2020-04-27 2020-08-07 山东大学 一种对数非均匀伪随机电磁勘探信号生成方法及系统

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
E-Hz广域电磁方法研究;佟铁钢;《中国博士学位论文全文数据库 基础科学辑》;20120215(第02期);14、16、20、22 *

Also Published As

Publication number Publication date
CN111522064A (zh) 2020-08-11
US11983511B2 (en) 2024-05-14
WO2021218524A1 (zh) 2021-11-04
US20220308838A1 (en) 2022-09-29

Similar Documents

Publication Publication Date Title
CN111522064B (zh) 一种高阶伪随机电磁勘探信号生成方法及系统
CN111505722B (zh) 一种对数非均匀伪随机电磁勘探信号生成方法及系统
CN111505723B (zh) 一种2n序列伪随机电磁勘探信号生成方法及系统
Porter et al. Finite-element ray tracing
Wang Multichannel matching pursuit for seismic trace decomposition
MXPA01010291A (es) Metodo de vigilancia sismica de una zona subterranea por utilizacion simultanea de varias fuentes vibro-sismicas.
Chong et al. Joint inversion of body‐wave receiver function and Rayleigh‐wave ellipticity
Varma et al. Comparison of binary black hole initial data sets
CN102062872A (zh) 可控震源多因子扫描方法
CN109061724A (zh) 一种基于自适应变分模态分解的地震数据降噪方法
CN110161311B (zh) 一种谐波与间谐波的检测方法
CN105675126A (zh) 一种用于检测多频多源复杂稳定声场声压的新方法
RU2434250C1 (ru) Способ регистрации сейсмических сигналов на акватории моря при поиске подводных залежей углеводородов
Kumar et al. Resolving phase wrapping by using sliding transform for generation of dispersion curves
Mostafanejad et al. Velocity structure of the northern Mississippi Embayment sediments, Part II: Inversion of teleseismic P‐wave transfer functions
Wu et al. Mesh-free radial-basis-function-generated finite differences and their application in reverse time migration
Porter et al. Anomalous parabolic equation results for propagation in leaky surface ducts
Hornbostel et al. Waveform design for electroseismic exploration
Islam Study of eccentric binary black hole mergers using numerical relativity and an inspiral-merger-ringdown model
CN107850997A (zh) Cs本振序列生成方法、装置、发射机及接收机
JP2001311781A (ja) 地盤速度構造の推定方法
Stiassnie et al. Prediction of long forcing waves for harbor agitation studies
Shi et al. Calculation of normal and leaky modes in planar waveguides based on a semi-analytical spectral element method
CN110349574A (zh) 基于声控设备的账户绑定方法及装置
Scholtz Pseudo-random sweeps for built-up area seismic surveys

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