CN111856552B - 可控震源扫描信号的生成方法及装置 - Google Patents

可控震源扫描信号的生成方法及装置 Download PDF

Info

Publication number
CN111856552B
CN111856552B CN201910362233.9A CN201910362233A CN111856552B CN 111856552 B CN111856552 B CN 111856552B CN 201910362233 A CN201910362233 A CN 201910362233A CN 111856552 B CN111856552 B CN 111856552B
Authority
CN
China
Prior art keywords
pseudo
amplitude spectrum
random
scanning signal
scanning
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
CN201910362233.9A
Other languages
English (en)
Other versions
CN111856552A (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.)
China National Petroleum Corp
BGP Inc
Original Assignee
China National Petroleum Corp
BGP Inc
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 China National Petroleum Corp, BGP Inc filed Critical China National Petroleum Corp
Priority to CN201910362233.9A priority Critical patent/CN111856552B/zh
Publication of CN111856552A publication Critical patent/CN111856552A/zh
Application granted granted Critical
Publication of CN111856552B publication Critical patent/CN111856552B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/02—Generating seismic energy
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/36—Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
    • 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
    • Y02D—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN INFORMATION AND COMMUNICATION TECHNOLOGIES [ICT], I.E. INFORMATION AND COMMUNICATION TECHNOLOGIES AIMING AT THE REDUCTION OF THEIR OWN ENERGY USE
    • Y02D30/00—Reducing energy consumption in communication networks
    • Y02D30/70—Reducing energy consumption in communication networks in wireless communication networks

Landscapes

  • Engineering & Computer Science (AREA)
  • Remote Sensing (AREA)
  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Acoustics & Sound (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

本发明公开了一种可控震源扫描信号的生成方法及装置，该方法包括：获取可控震源的扫描参数，其中，扫描参数至少包括：扫描长度和采样间隔；根据扫描长度和采样间隔，计算采样点数N；生成预设区间范围且均匀分布的N个伪随机数；根据N个伪随机数，确定一个伪随机序列；基于线性扫描信号的振幅谱，对伪随机序列的振幅谱进行整形处理，得到可控震源的伪随机扫描信号，其中，线性扫描信号是根据所述可控震源的扫描参数生成的线性信号。本发明能够为可控震源提供一种振幅谱平坦的伪随机扫描信号，使得可控震源能够用于建筑设施敏感地区的地震资料采集。

Description

可控震源扫描信号的生成方法及装置
技术领域
本发明涉及石油勘探领域，尤其涉及一种可控震源扫描信号的生成方法及装置。
背景技术
本部分旨在为权利要求书中陈述的本发明实施例提供背景或上下文。此处的描述不因为包括在本部分中就承认是现有技术。
对于陆上石油地震勘探，可控震源是一种广泛应用的安全、环保的激发装备。可控震源在振动时需要预先设计一个参考扫描信号，震源的电控系统根据该参考扫描信号驱动震源进行振动。
一般预先设计的参考扫描信号有如下三种类型：线性扫描信号、非线性扫描信号和伪随机扫描信号，其中，线性扫描信号和非线性扫描是指扫描速率是线性的或非线性的，而伪随机扫描信号的扫描速率是伪随机变化的。
线性扫描信号具有能量强、振幅谱平坦无抖动等优点，是应用最广泛的扫描信号；非线性扫描信号主要通过补偿某些频段的能量，应用于低频扩展或高频补偿的特殊项目中；伪随机扫描信号由于扫描速率伪随机变化，不会与建筑设施发生共振，可以减小对建筑设计的破坏，主要应用于建筑设施敏感地区的施工。
目前，现有伪随机扫描信号设计方法生成的伪随机扫描信号，存在振幅谱抖动剧烈，以及相关子波旁瓣能量强等缺点，使得伪随机扫描的应用受到限制。
发明内容
本发明实施例提供一种可控震源扫描信号的生成方法，用以解决现有技术中，用于可控震源的伪随机扫描信号振幅谱抖动剧烈的技术问题，该方法包括：获取可控震源的扫描参数，其中，扫描参数至少包括：扫描长度和采样间隔；根据扫描长度和采样间隔，计算采样点数N；生成预设区间范围且均匀分布的N个伪随机数；根据N个伪随机数，确定一个伪随机序列；基于线性扫描信号的振幅谱，对伪随机序列的振幅谱进行整形处理，得到可控震源的伪随机扫描信号，其中，线性扫描信号是根据所述可控震源的扫描参数生成的线性信号；
通过如下公式对伪随机序列的振幅谱进行整形处理，得到整形后的振幅谱：
Pshape(f)＝PRj(f)×Ptarget(f)/(PRj(f)+ε)
其中，Pshape(f)表示整形后的振幅谱；PRj(f)表示伪随机序列的振幅谱；Ptarget(f)表示目标振幅谱；ε表示白噪声系数；目标振幅谱由线性扫描信号的振幅谱进行归一化处理后得到；
根据扫描参数确定线性扫描信号通过如下公式来实现：
其中，
其中，S(t)表示线性扫描信号；A(t)表示Blackman窗函数或Cosine函数定义的振幅包络；T表示扫描长度；f0表示扫描的起始频率；fm表示扫描的终止频率；f(t)表示瞬时频率变换函数。
本发明实施例还提供一种可控震源扫描信号的生成装置，用以解决现有技术中，用于可控震源的伪随机扫描信号振幅谱抖动剧烈的技术问题，该装置包括：扫描参数获取模块，用于获取可控震源的扫描参数，其中，扫描参数至少包括：扫描长度和采样间隔；采样点数计算模块，用于根据扫描长度和采样间隔，计算采样点数N；伪随机数生成模块，用于生成预设区间范围且均匀分布的N个伪随机数；伪随机序列确定模块，用于根据N个伪随机数，确定一个伪随机序列；伪随机扫描信号生成模块，用于基于线性扫描信号的振幅谱，对伪随机序列的振幅谱进行整形处理，得到可控震源的伪随机扫描信号，其中，线性扫描信号是根据所述可控震源的扫描参数生成的线性信号；
通过如下公式对伪随机序列的振幅谱进行整形处理，得到整形后的振幅谱：
Pshape(f)＝PRj(f)×Ptarget(f)/(PRj(f)+ε)
其中，Pshape(f)表示整形后的振幅谱；PRj(f)表示伪随机序列的振幅谱；Ptarget(f)表示目标振幅谱；ε表示白噪声系数；目标振幅谱由线性扫描信号的振幅谱进行归一化处理后得到；
根据扫描参数确定线性扫描信号通过如下公式来实现：
其中，
其中，S(t)表示线性扫描信号；A(t)表示Blackman窗函数或Cosine函数定义的振幅包络；T表示扫描长度；f0表示扫描的起始频率；fm表示扫描的终止频率；f(t)表示瞬时频率变换函数。
本发明实施例还提供一种计算机设备，用以解决现有技术中，用于可控震源的伪随机扫描信号振幅谱抖动剧烈的技术问题，该计算机设备包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，处理器执行计算机程序时实现上述可控震源扫描信号的生成方法。
本发明实施例还提供一种计算机可读存储介质，用以解决现有技术中，用于可控震源的伪随机扫描信号振幅谱抖动剧烈的技术问题，该计算机可读存储介质存储有执行上述可控震源扫描信号的生成方法的计算机程序。
本发明实施例中，根据可控震源的扫描参数生成一个伪随机序列，并根据该可控震源的扫描参数生成一个线性扫描信号，然后基于该线性扫描信号的振幅谱，对伪随机序列的振幅谱进行整形处理，从而得到与线性扫描信号振幅谱接近的伪随机扫描信号，以作为可控震源振动时的参考扫描信号。
通过本发明实施例，能够为可控震源提供一种振幅谱平坦的伪随机扫描信号，更加接近相同频宽的线性扫描信号的振幅谱，可以应用于建筑设施敏感地区的可控震源施工。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。在附图中：
图1为本发明实施例中提供的一种可控震源扫描信号的生成方法流程图；
图2为本发明实施例中提供的一种线性扫描信号示意图；
图3为本发明实施例中提供的一种目标振幅谱示意图；
图4为本发明实施例中提供的一种伪随机扫描信号示意图；
图5为本发明实施例中提供的一种线性扫描信号与伪随机扫描信号子波频谱对比示意图；
图6为本发明实施例中提供的一种采用6Vibpro箱体生成的伪随机扫描信号子波频谱示意图；
图7为本发明实施例中提供的一种采用7VE464箱体生成的伪随机扫描信号子波频谱示意图；
图8为本发明实施例中提供的一种可控震源扫描信号的生成装置示意图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚明白，下面结合附图对本发明实施例做进一步详细说明。在此，本发明的示意性实施例及其说明用于解释本发明，但并不作为对本发明的限定。
在本说明书的描述中，所使用的“包含”、“包括”、“具有”、“含有”等，均为开放性的用语，即意指包含但不限于。参考术语“一个实施例”、“一个具体实施例”、“一些实施例”、“例如”等的描述意指结合该实施例或示例描述的具体特征、结构或者特点包含于本申请的至少一个实施例或示例中。在本说明书中，对上述术语的示意性表述不一定指的是相同的实施例或示例。而且，描述的具体特征、结构或者特点可以在任何的一个或多个实施例或示例中以合适的方式结合。各实施例中涉及的步骤顺序用于示意性说明本申请的实施，其中的步骤顺序不作限定，可根据需要作适当调整。
本发明实施例中提供了一种可控震源扫描信号的生成方法，图1为本发明实施例中提供的一种可控震源扫描信号的生成方法流程图，如图1所示，该方法包括如下步骤：
S101，获取可控震源的扫描参数，其中，扫描参数至少包括：扫描长度和采样间隔。
需要说明的是，通过S101获取的可控震源的扫描参数，可以是用户预先设置的扫描参数，需要注意的是，预先设置的扫描参数除了上述扫描长度和采样间隔外，还可以包括但不限于：起始频率、终止频率(起始频率小于终止频率)、斜坡类型、起始斜坡长度、终止斜坡长度。可选地，设置的斜坡类型可以是但不限于Blackman窗函数或Cosine函数，用来定义扫描信号振幅包络。
S102，根据扫描长度和采样间隔，计算采样点数N。
需要说明的是，采样点数N为待生成伪随机扫描信号的采样点数。作为一种可选的实施方式，采样点数的计算方法为：
N＝T/dt (1)
其中，T表示扫描长度；dt表示采样间隔。
S103，生成预设区间范围且均匀分布的N个伪随机数。
作为一种可选的实施方式，上述S103可以生成均匀分布且范围在[-π,π]之间的N个伪随机数Rj：
Rj＝Random(j)其中，-π≤Rj≤π (2)
S104，根据N个伪随机数，确定一个伪随机序列。
以均匀分布且范围在[-π,π]之间的N个伪随机数Rj为例，作为一种可选的实施方式，上述可以具体包括：计算N个伪随机数的正弦值或余弦值；将N个伪随机数的正弦值或余弦值，确定为伪随机序列。
如果将N个伪随机数的正弦值作为伪随机序列，则可以通过下面的公式确定一个伪随机序列：
SRj＝sin(Rj) (3)
其中，SRj表示正弦值伪随机序列。
S105，基于线性扫描信号的振幅谱，对伪随机序列的振幅谱进行整形处理，得到可控震源的伪随机扫描信号，其中，线性扫描信号是根据所述可控震源的扫描参数生成的线性信号。
作为一种可选的实施例，上述S105可以具体包括如下步骤：计算线性扫描信号的振幅谱；对线性扫描信号的振幅谱进行归一化处理，得到目标振幅谱；基于目标振幅谱，对伪随机序列的振幅谱进行整形处理，得到可控震源的伪随机扫描信号。
进一步地，可以通过如下步骤来计算线性扫描信号的振幅谱：根据扫描参数确定线性扫描信号；对线性扫描信号进行傅里叶变换运算，得到线性扫描信号对应的复数序列；根据线性扫描信号对应的复数序列，计算线性扫描信号的振幅谱。
其中，根据扫描参数确定线性扫描信号可以通过如下公式来实现：
其中，
其中，S(t)表示线性扫描信号；A(t)表示Blackman窗函数或Cosine函数定义的振幅包络；T表示扫描长度；f0表示扫描的起始频率；fm表示扫描的终止频率；f(t)表示瞬时频率变换函数。
作为第一种可选的实施方案，当斜坡类型为Blackman窗函数的情况下，振幅包络A(t)可以表示为：
其中，
作为第二种可选的实施方案，当斜坡类型为Cosine函数的情况下，振幅包络A(t)可以表示为：
其中，
其中，T1表示起始斜坡长度；T2表示终止斜坡长度。
在根据扫描参数确定线性扫描信号之后，可以对线性扫描信号S(t)进行傅里叶变换运算，得到如下复数序列：
Cj＝aj+ibj (6)
其中，j＝1,2,…,N；i表示虚数单位。
进一步地，可以通过如下公式计算线性扫描信号的振幅谱PLj：
通过下式对线性扫描信号的振幅谱PLj进行归一化处理，得到目标振幅谱：
在一种可选的实施例中，基于目标振幅谱，对伪随机序列的振幅谱进行整形处理，得到可控震源的伪随机扫描信号，具体可以包括如下步骤：对伪随机序列进行傅里叶变换运算，得到伪随机序列对应的复数序列；根据伪随机序列对应的复数序列，计算伪随机序列的振幅谱；基于目标振幅谱，对伪随机序列的振幅谱进行整形处理，得到整形后的振幅谱；对整形后的振幅谱进行傅里叶逆变换运算，得到整形信号；对整形信号进行归一化处理，得到可控震源的伪随机扫描信号。
以正弦值伪随机序列SRj为例，对其进行傅里叶变换运算，可得到如下复数序列：
CRj＝cj+idj (9)
其中，j＝1,2,…,N；i表示虚数单位。
进一步地，可以通过如下公式计算伪随机序列SRj的振幅谱PRj：
在基于目标振幅谱，对伪随机序列的振幅谱进行整形处理，得到整形后的振幅谱的时候，作为一种可选的实施方式，通过如下公式对伪随机序列的振幅谱进行整形处理，得到整形后的振幅谱：
Pshape(f)＝PRj(f)×Ptarget(f)/(PRj(f)+ε) (11)
其中，Pshape(f)表示整形后的振幅谱；PRj(f)表示伪随机序列的振幅谱；Ptarget(f)表示目标振幅谱；ε表示白噪声系数。
进一步地，对整形后的振幅谱Pshape进行傅里叶逆变换运算，可得到整形信号shape(t)；对整形信号shape(t)进行归一化处理，可得到可控震源的伪随机扫描信号R(t)：
由上可知，本发明实施例中，根据可控震源的扫描参数生成一个伪随机序列，并根据该可控震源的扫描参数生成一个线性扫描信号，然后基于该线性扫描信号的振幅谱，对伪随机序列的振幅谱进行整形处理，从而得到与线性扫描信号振幅谱接近的伪随机扫描信号，以作为可控震源振动时的参考扫描信号。
通过本发明实施例提供的可控震源扫描信号的生成方法，能够为可控震源提供一种振幅谱平坦的伪随机扫描信号，更加接近相同频宽的线性扫描信号的振幅谱，可以应用于建筑设施敏感地区的可控震源施工。
作为一种优选的实施方式，本发明实施例提供的可控震源伪随机扫描信号的生成方法可以通过如下步骤来实现：
①设置可控震源扫描参数；
②计算伪随机扫描信号采样点数N；
③生成均匀分布的范围在[-π，π]之间的N个伪随机数；
④计算[-π，π]之间的N个伪随机数的正弦值或余弦值，得到伪随机序列；
⑤计算线性扫描信号的振幅谱，作为目标谱；
⑥计算伪随机序列的振幅谱，并对其进行谱整形，得到整形信号；
⑦对整形信号进行归一化处理，得到伪随机扫描信号。
其中，可控震源扫描参数包括起始频率f0、终止频率fm、扫描长度T、斜坡类型、起始斜坡长度T1、终止斜坡长度T2和采样间隔dt。其中，扫描的起始频率小于终止频率；斜坡类型一般为Blackman窗函数或Cosine函数，用来定义扫描信号振幅包络。
假设可控震源的扫描参数为：起始频率f0＝5Hz、终止频率fm＝96Hz、扫描长度T＝16s、斜坡类型采用Blackman窗函数、起始斜坡长度T1＝250ms、终止斜坡长度T2＝250ms和采样间隔dt＝0.5ms，则根据扫描长度和采样间隔计算得到待生成伪随机扫描信号的采样点数为32000。
图2所示为根据可控震源的扫描参数生成的线性扫描信号；对图2所示的线性扫描信号进行归一化处理，可得到目标振幅谱，如图3所示。图4示出了基于图3所示的目标振幅谱对根据可控震源的扫描参数生成的伪随机序列进行整形处理后的伪随机扫描信号。图5是线性扫描信号和伪随机扫描信号的频谱对比示意图。
图6是常见VibPro震源箱体采用相同扫描参数生成伪随机扫描信号的频谱示意图，图7是常见VE464震源箱体采用相同扫描参数生成伪随机扫描信号的频谱示意图；由图6和图7可以看出，VibPro箱体和VE464箱体生成的伪随机扫描信号的频谱抖动剧烈，而采用本发明实施例生产伪随机扫描信号的子波频谱抖动很小，与线性扫描信号的子波频谱很接近。
本发明实施例中还提供了一种可控震源扫描信号的生成装置，如下面的实施例所述。由于该装置实施例解决问题的原理与可控震源扫描信号的生成方法相似，因此该装置实施例的实施可以参见方法的实施，重复之处不再赘述。
图8为本发明实施例中提供的一种可控震源扫描信号的生成装置示意图，如图8所示，该装置包括：扫描参数获取模块81、采样点数计算模块82、伪随机数生成模块83、伪随机序列确定模块84和伪随机扫描信号生成模块85。
其中，扫描参数获取模块81，用于获取可控震源的扫描参数，其中，扫描参数至少包括：扫描长度和采样间隔；采样点数计算模块82，用于根据扫描长度和采样间隔，计算采样点数N；伪随机数生成模块83，用于生成预设区间范围且均匀分布的N个伪随机数；伪随机序列确定模块84，用于根据N个伪随机数，确定一个伪随机序列；伪随机扫描信号生成模块85，用于基于线性扫描信号的振幅谱，对伪随机序列的振幅谱进行整形处理，得到可控震源的伪随机扫描信号，其中，线性扫描信号是根据所述可控震源的扫描参数生成的线性信号。
由上可知，本发明实施例提供的可控震源扫描信号的生成装置，根据可控震源的扫描参数生成一个伪随机序列，并根据该可控震源的扫描参数生成一个线性扫描信号，然后基于该线性扫描信号的振幅谱，对伪随机序列的振幅谱进行整形处理，从而得到与线性扫描信号振幅谱接近的伪随机扫描信号，以作为可控震源振动时的参考扫描信号。
通过本发明实施例提供的可控震源扫描信号的生成装置，能够为可控震源提供一种振幅谱平坦的伪随机扫描信号，更加接近相同频宽的线性扫描信号的振幅谱，可以应用于建筑设施敏感地区的可控震源施工。
在一种可选的实施例中，本发明实施例生成的伪随机数的区间范围为[-π,π]，则上述伪随机序列确定模块84还用于计算N个伪随机数的正弦值或余弦值；将N个伪随机数的正弦值或余弦值，确定为伪随机序列。
在一种可选的实施例中，本发明实施例提供的可控震源扫描信号的生成装置中，伪随机扫描信号生成模块85还包括：线性扫描信号振幅计算子模块，用于计算线性扫描信号的振幅谱；目标振幅谱计算子模块，用于对线性扫描信号的振幅谱进行归一化处理，得到目标振幅谱；伪随机扫描信号生成子模块，用于基于目标振幅谱，对伪随机序列的振幅谱进行整形处理，得到可控震源的伪随机扫描信号。
基于上述实施例，作为一种可选的实施方式，线性扫描信号振幅计算子模块还用于根据扫描参数确定线性扫描信号；对线性扫描信号进行傅里叶变换运算，得到线性扫描信号对应的复数序列；根据线性扫描信号对应的复数序列，计算线性扫描信号的振幅谱。
需要说明的是，本发明实施例提供的可控震源扫描信号的生成装置中，扫描参数还可以包括：起始频率、终止频率、斜坡类型、起始斜坡长度、终止斜坡长度，其中，斜坡类型为Blackman窗函数或Cosine函数，线性扫描信号振幅计算子模块还用于通过如下公式根据扫描参数确定线性扫描信号：
其中，
其中，
其中，S(t)表示线性扫描信号；A(t)表示Blackman窗函数的振幅包络；T1表示起始斜坡长度；T2表示终止斜坡长度；T表示扫描长度；f0表示起始频率；fm表示终止频率；f(t)表示瞬时频率变换函数。
基于上述实施例，作为一种可选的实施方式，伪随机扫描信号生成子模块还用于对伪随机序列进行傅里叶变换运算，得到伪随机序列对应的复数序列；根据伪随机序列对应的复数序列，计算伪随机序列的振幅谱；基于目标振幅谱，对伪随机序列的振幅谱进行整形处理，得到整形后的振幅谱；对整形后的振幅谱进行傅里叶逆变换运算，得到整形信号；对整形信号进行归一化处理，得到可控震源的伪随机扫描信号。
可选地，伪随机扫描信号生成子模块还用于通过如下公式对伪随机序列的振幅谱进行整形处理，得到整形后的振幅谱：
Pshape(f)＝PRj(f)×Ptarget(f)/(PRj(f)+ε)；
其中，Pshape(f)表示整形后的振幅谱；PRj(f)表示伪随机序列的振幅谱；Ptarget(f)表示目标振幅谱；ε表示白噪声系数。
本发明实施例中还提供了一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，处理器执行计算机程序时实现上述方法实施例中任意一种可选的或优选的可控震源扫描信号的生成方法。
本发明实施例中还提供了一种计算机可读存储介质，计算机可读存储介质存储有执行上述方法实施例中任意一种可选的或优选的可控震源扫描信号的生成方法的计算机程序。
综上所述，本发明实施例提供了一种可控震源伪随机扫描信号的生成方法，根据可控震源的扫描参数生成线性扫描信号，并计算线性扫描信号的振幅谱，将该振幅谱作为待生成的伪随机扫描信号的目标振幅谱；对伪随机序列按照目标谱进行谱整形，得到伪随机扫描信号。将本发明实施例提供的可控震源伪随机扫描信号用于可控震源激发的地震勘探，特别是在建筑设施附近进行施工时，可以减小安全距离，缩小采集资料缺口，与直接采用线性扫描信号的可控震源相比，采用本发明实施例提供的伪随机扫描信号的可控震源，可以减轻对建筑设施的破坏。
本领域内的技术人员应明白，本发明的实施例可提供为方法、系统、或计算机程序产品。因此，本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本发明是参照根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
以上所述的具体实施例，对本发明的目的、技术方案和有益效果进行了进一步详细说明，所应理解的是，以上所述仅为本发明的具体实施例而已，并不用于限定本发明的保护范围，凡在本发明的精神和原则之内，所做的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。

Claims (8)

1.一种可控震源扫描信号的生成方法，其特征在于，包括：
获取可控震源的扫描参数，其中，所述扫描参数至少包括：扫描长度和采样间隔；
根据扫描长度和采样间隔，计算采样点数N；
生成预设区间范围且均匀分布的N个伪随机数；
根据所述N个伪随机数，确定一个伪随机序列；
基于线性扫描信号的振幅谱，对所述伪随机序列的振幅谱进行整形处理，得到所述可控震源的伪随机扫描信号，其中，所述线性扫描信号是根据所述可控震源的扫描参数生成的线性信号；
通过如下公式对伪随机序列的振幅谱进行整形处理，得到整形后的振幅谱：
Pshape(f)＝PRj(f)×Ptarget(f)/(PRj(f)+ε)
其中，Pshape(f)表示整形后的振幅谱；PRj(f)表示伪随机序列的振幅谱；Ptarget(f)表示目标振幅谱；ε表示白噪声系数；目标振幅谱由线性扫描信号的振幅谱进行归一化处理后得到；
根据扫描参数确定线性扫描信号通过如下公式来实现：
其中，
其中，S(t)表示线性扫描信号；A(t)表示Blackman窗函数或Cosine函数定义的振幅包络；T表示扫描长度；f0表示扫描的起始频率；fm表示扫描的终止频率；f(t)表示瞬时频率变换函数。
2.如权利要求1所述的方法，其特征在于，所述预设区间范围为[-π,π]，其中，根据所述N个伪随机数，确定一个伪随机序列，包括：
计算所述N个伪随机数的正弦值或余弦值；
将所述N个伪随机数的正弦值或余弦值，确定为所述伪随机序列。
3.如权利要求1所述的方法，其特征在于，基于线性扫描信号的振幅谱，对所述伪随机序列的振幅谱进行整形处理，得到所述可控震源的伪随机扫描信号，包括：
计算线性扫描信号的振幅谱；
对所述线性扫描信号的振幅谱进行归一化处理，得到目标振幅谱；
基于所述目标振幅谱，对所述伪随机序列的振幅谱进行整形处理，得到所述可控震源的伪随机扫描信号。
4.如权利要求3所述的方法，其特征在于，计算线性扫描信号的振幅谱，包括：
根据所述扫描参数确定线性扫描信号；
对所述线性扫描信号进行傅里叶变换运算，得到所述线性扫描信号对应的复数序列；
根据所述线性扫描信号对应的复数序列，计算所述线性扫描信号的振幅谱。
5.如权利要求3所述的方法，其特征在于，基于所述目标振幅谱，对所述伪随机序列的振幅谱进行整形处理，得到所述可控震源的伪随机扫描信号，包括：
对所述伪随机序列进行傅里叶变换运算，得到所述伪随机序列对应的复数序列；
根据所述伪随机序列对应的复数序列，计算所述伪随机序列的振幅谱；
基于所述目标振幅谱，对所述伪随机序列的振幅谱进行整形处理，得到整形后的振幅谱；
对所述整形后的振幅谱进行傅里叶逆变换运算，得到整形信号；
对所述整形信号进行归一化处理，得到所述可控震源的伪随机扫描信号。
6.一种可控震源扫描信号的生成装置，其特征在于，包括：
扫描参数获取模块，用于获取可控震源的扫描参数，其中，所述扫描参数至少包括：扫描长度和采样间隔；
采样点数计算模块，用于根据扫描长度和采样间隔，计算采样点数N；
伪随机数生成模块，用于生成预设区间范围且均匀分布的N个伪随机数；
伪随机序列确定模块，用于根据所述N个伪随机数，确定一个伪随机序列；
伪随机扫描信号生成模块，用于基于线性扫描信号的振幅谱，对所述伪随机序列的振幅谱进行整形处理，得到所述可控震源的伪随机扫描信号，其中，所述线性扫描信号是根据所述可控震源的扫描参数生成的线性信号；
通过如下公式对伪随机序列的振幅谱进行整形处理，得到整形后的振幅谱：
Pshape(f)＝PRj(f)×Ptarget(f)/(PRj(f)+ε)
其中，Pshape(f)表示整形后的振幅谱；PRj(f)表示伪随机序列的振幅谱；Ptarget(f)表示目标振幅谱；ε表示白噪声系数；目标振幅谱由线性扫描信号的振幅谱进行归一化处理后得到；
根据扫描参数确定线性扫描信号通过如下公式来实现：
其中，
其中，S(t)表示线性扫描信号；A(t)表示Blackman窗函数或Cosine函数定义的振幅包络；T表示扫描长度；f0表示扫描的起始频率；fm表示扫描的终止频率；f(t)表示瞬时频率变换函数。
7.一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，其特征在于，所述处理器执行所述计算机程序时实现权利要求1至5任一项所述可控震源扫描信号的生成方法。
8.一种计算机可读存储介质，其特征在于，所述计算机可读存储介质存储有执行权利要求1至5任一项所述可控震源扫描信号的生成方法的计算机程序。
CN201910362233.9A 2019-04-30 2019-04-30 可控震源扫描信号的生成方法及装置 Active CN111856552B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910362233.9A CN111856552B (zh) 2019-04-30 2019-04-30 可控震源扫描信号的生成方法及装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910362233.9A CN111856552B (zh) 2019-04-30 2019-04-30 可控震源扫描信号的生成方法及装置

Publications (2)

Publication Number Publication Date
CN111856552A CN111856552A (zh) 2020-10-30
CN111856552B true CN111856552B (zh) 2023-07-25

Family

ID=72966654

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910362233.9A Active CN111856552B (zh) 2019-04-30 2019-04-30 可控震源扫描信号的生成方法及装置

Country Status (1)

Country Link
CN (1) CN111856552B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114384580B (zh) * 2021-12-31 2023-05-02 同济大学 一种基于可控震源的理想子波定制方法

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4093988A (en) * 1976-11-08 1978-06-06 General Electric Company High speed frequency response measurement
CN102062872A (zh) * 2009-11-11 2011-05-18 中国石油天然气集团公司 可控震源多因子扫描方法
CN102426386A (zh) * 2011-10-27 2012-04-25 吉林大学 脉冲编码可控震源的多维匹配冲击方法
CN104035128A (zh) * 2013-03-06 2014-09-10 中国石油集团东方地球物理勘探有限责任公司 可控震源伪随机扫描信号生成方法
CN104217099A (zh) * 2014-07-18 2014-12-17 中国石油化工股份有限公司 一种用于产生并处理可控震源的扫描信号的方法
CN104950326A (zh) * 2015-03-04 2015-09-30 中石化石油工程地球物理有限公司胜利分公司 基于目的层频谱的可控震源非线性扫描信号的设计方法
CN105467434A (zh) * 2016-01-15 2016-04-06 吉林大学 冲击夯可控震源变频伪随机编码扫描的控制方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9568552B2 (en) * 2013-06-06 2017-02-14 Mentor Graphics Corporation Logic built-in self-test with high test coverage and low switching activity

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4093988A (en) * 1976-11-08 1978-06-06 General Electric Company High speed frequency response measurement
CN102062872A (zh) * 2009-11-11 2011-05-18 中国石油天然气集团公司 可控震源多因子扫描方法
CN102426386A (zh) * 2011-10-27 2012-04-25 吉林大学 脉冲编码可控震源的多维匹配冲击方法
CN104035128A (zh) * 2013-03-06 2014-09-10 中国石油集团东方地球物理勘探有限责任公司 可控震源伪随机扫描信号生成方法
CN104217099A (zh) * 2014-07-18 2014-12-17 中国石油化工股份有限公司 一种用于产生并处理可控震源的扫描信号的方法
CN104950326A (zh) * 2015-03-04 2015-09-30 中石化石油工程地球物理有限公司胜利分公司 基于目的层频谱的可控震源非线性扫描信号的设计方法
CN105467434A (zh) * 2016-01-15 2016-04-06 吉林大学 冲击夯可控震源变频伪随机编码扫描的控制方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
伪随机扫描在可控震源上的应用探讨;马铁荣;《物探装备》;20161031;第26卷(第5期);319-323 *
伪随机扫描在可控震源探测中的应用;T.Dean;《世界地震译丛》;20160131;第47卷(第1期);67-89 *

Also Published As

Publication number Publication date
CN111856552A (zh) 2020-10-30

Similar Documents

Publication Publication Date Title
CN107389478B (zh) 一种基于小波包分析的材料疲劳寿命的预测方法
CN109782335B (zh) 低畸变宽频扫描信号设计方法
US9160448B2 (en) Apparatus and method for estimating intra-channel nonlinear damage
JP2017537318A (ja) フーリエ変換質量分析法
CN111856552B (zh) 可控震源扫描信号的生成方法及装置
US20180252741A1 (en) Method and device for determining machine speeds
CN110596758B (zh) 一种地震信号低频能量补偿方法
Catherall et al. High resolution spectrograms using a component optimized short-term fractional Fourier transform
CN105675126A (zh) 一种用于检测多频多源复杂稳定声场声压的新方法
CN113359207B (zh) 基于太赫兹雷达的声致水面微动特征提取方法及装置
Kupriyanova et al. Long-period pulsations of the thermal microwave emission of the solar flare of June 2, 2007 from data with high spatial resolution
CN107367758B (zh) 一种提高信噪比的可控震源加权相关的方法
Sun et al. Wavelet transform based de-noising method for self mixing interferometry signals
CN114019236A (zh) 一种电网谐波单通道混叠目标信号检测方法和装置
CN107219551A (zh) 拓宽地震数据频带的方法及装置
Baidakov et al. Estimation of broadening of the spectra of short surface waves in the presence of long waves
Sergeev et al. Temporal evolution of HF-excited plasma waves, measured at different pump frequencies by stimulated electromagnetic emission (SEE)
JP6182836B2 (ja) 信号処理装置及び方法、並びにプログラム及び記録媒体
Sergeev et al. Study of the plasma turbulence dynamics by measurements of diagnostic stimulated electromagnetic emission. I. Experimental results
Ochiai et al. Impulse response measurement with constant signal-to-noise ratio over a wide frequency range
CN113126166B (zh) 一种基于分段线性调频的声波测井偶极声源自动调节方法
Scholtz Optimum pseudo-random sweeps for vibratory surveys to reduce risk of damaging buildings
CN116992254B (zh) 变频电机的壳振信号的重构方法、装置、设备及存储介质
Scholtz Pseudo-random sweeps for built-up area seismic surveys
Bork Exponential sweep check using hilbert transform

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