CN104793237A - 一种获得宽频可控震源扫描信号的方法和装置 - Google Patents

一种获得宽频可控震源扫描信号的方法和装置 Download PDF

Info

Publication number
CN104793237A
CN104793237A CN201510108954.9A CN201510108954A CN104793237A CN 104793237 A CN104793237 A CN 104793237A CN 201510108954 A CN201510108954 A CN 201510108954A CN 104793237 A CN104793237 A CN 104793237A
Authority
CN
China
Prior art keywords
spectrum
sweep signal
earth
coupling response
compensating
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.)
Pending
Application number
CN201510108954.9A
Other languages
English (en)
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 CN201510108954.9A priority Critical patent/CN104793237A/zh
Publication of CN104793237A publication Critical patent/CN104793237A/zh
Pending legal-status Critical Current

Links

Landscapes

  • Geophysics And Detection Of Objects (AREA)

Abstract

本申请实施例提供一种获得宽频可控震源扫描信号的方法和装置。该方法包括:获得通过预设扫描信号扫描得到的近道共偏移距道集;根据所述近道共偏移距道集进行统计子波计算,获得大地震动耦合响应谱;获得预设扫描信号谱;利用所述大地震动耦合响应谱对所述预设扫描信号谱进行能量补偿,并根据补偿后的频谱获得第一目标扫描信号。

Description

一种获得宽频可控震源扫描信号的方法和装置
技术领域
本发明涉及地球物理勘探技术领域,尤其涉及一种获得宽频可控震源扫描信号的方法和装置。
背景技术
一般地震勘探中常用的震源包括人工地面震源、炸药震源和可控震源。由于可控震源所产生的信号频谱和基本特性可以人为控制,可以在设计震源扫描信号时避开某些干扰频率,还能对地层对地震信号的吸收作用进行补偿,这是其它人工地面震源和炸药震源难于做到的,所以利用可控震源进行地震勘探可以得到反射能量足够,信噪比和信号分辨率能够满足地质勘探需要的资料。因此,可控震源技术在地震勘探中得到了广泛的应用。
在利用可控震源进行地震勘探中,要求采用震源机械-液压系统能够响应并可物理实现的信号,即信号频率宽度有限,其最低频率大于可控震源振动器所能激发信号的最低频率;最高频率不能超出振动器所能激发信号频率的上限,并且震源所激发信号的频带应在大地可传输信号通频带内,且信号具有良好分辨率的零相位相关子波。这种信号的振幅和频率都是时间的函数,称这种信号为扫描信号或扫频信号。
在地震勘探中,利用扫描信号进行扫描后能有效带回多少地下信息的决定性因素是地震波频宽。首先地震波反演技术的应用为利用地震属性直接确定多孔介质的物性参数、流体含量、边界以及孔隙度等提供了一种可靠的工具,而且只有包含低频(3Hz以下)的宽频地震波才能保证这一技术的可靠性。
但可控震源传到地下的信号通常不是完全与预先给定的扫描信号频率成分相同。可控震源传入大地的实际信号由震源系统的结构、以及震源系统与地表接触面的耦合条件等决定的。尤其对于地震波低频能量,震源系统与地表接触面的耦合条件会导致无法获得正确的地震波输出响应,也即无法获得均衡输出的扫描信号谱,因而无法得到包含低频(3Hz以下)的宽频地震波,也就不能准确的确定多孔介质的物性参数、流体含量、边界以及孔隙度等地质信息。
发明内容
本申请的目的是提供一种获得宽频可控震源扫描信号的方法和装置,以解决可控震源系统与大地震动耦合过程中对原设计信号的改造问题,使地震波的宽频均衡输出。
为了实现上述目的,本申请提供了一种获得宽频可控震源扫描信号的方法,该方法包括:
获得通过预设扫描信号扫描得到的近道共偏移距道集;
根据所述近道共偏移距道集进行统计子波计算,获得大地震动耦合响应谱;
获得预设扫描信号谱;
利用所述大地震动耦合响应谱对所述预设扫描信号谱进行能量补偿,并根据补偿后的频谱获得第一目标扫描信号。
在一个优选的实施例中,该方法还包括:获得震源出产输出特性谱;
利用所述震源出产输出特性谱和所述大地震动耦合响应谱对所述预设扫描信号谱进行能量补偿,并根据补偿后的频谱获得第二目标扫描信号。
在一个优选的实施例中,所述根据所述近道共偏移距道集进行统计子波计算,获得大地震动耦合响应谱包括:
获得所述近道共偏移距道集中每个接收点的地震子波,把所有地震子波进行最小二乘法拟合获得统计子波;
将所述统计子波进行傅里叶变换得到频谱,所述频谱作为大地震动耦合响应谱。
在一个优选的实施例中,所述利用所述大地震动耦合响应谱对所述预设扫描信号谱进行能量补偿,并根据补偿后的频谱获得第一目标扫描信号包括:
S11:根据所述大地震动耦合响应谱以第一预设迭代补偿公式对所述预设扫描信号谱进行能量补偿;
S12:获得当前频谱,将所述当前频谱进行傅里叶反变换获得当前扫描信号;
S13:判断根据所述当前扫描信号扫描大地后得到的地震波是否为宽频均衡输出;
S14:当所述判断为是时,所述当前扫描信号作为第一目标扫描信号;
S15:当所述判断为否时,根据所述大地震动耦合响应谱以第一预设迭代补偿公式对当前频谱进行能量补偿,并进入步骤S12。
在一个优选的实施例中,所述第一预设迭代补偿公式如下:
S n + 1 ( f ) = S n ( f ) E ( f )
其中,Sn(f)为第n次补偿迭代后的频谱;n为迭代的次数(n=0、1、2…),当n=0时,Sn(f)代表预设扫描信号谱;Sn+1(f)为第n+1次补偿迭代后的频谱;E(f)为大地震动耦合响应谱。
在一个优选的实施例中,所述利用所述震源出产输出特性谱和所述大地震动耦合响应谱对所述预设扫描信号谱进行能量补偿,并根据补偿后的频谱获得第二目标扫描信号包括:
S21:根据所述大地震动耦合响应谱以第二预设迭代补偿公式对所述预设扫描信号谱进行能量补偿;
S22:获得当前频谱,将所述当前频谱进行傅里叶反变换获得当前扫描信号;
S23:判断根据所述当前扫描信号扫描大地后得到的地震波是否为宽频均衡输出;
S24:当所述判断为是时,所述当前扫描信号作为第二目标扫描信号;
S25:当所述判断为否时,根据所述大地震动耦合响应谱以第二预设迭代补偿公式对当前频谱进行能量补偿,并进入步骤S22。
在一个优选的实施例中,所述第二预设迭代补偿公式如下:
S m + 1 ( f ) = S m ( f ) E ( f ) × G ( f )
其中,Sm(f)为第m次补偿迭代后的频谱;m为迭代的次数(m=0、1、2…),当m=0时,Sm(f)代表预设扫描信号谱;Sm+1(f)为第m+1次补偿迭代后的频谱;E(f)为大地震动耦合响应谱;G(f)震源出产输出特性谱。
本申请另一方面还提供一种获得宽频可控震源扫描信号的装置,该装置包括:
近道共偏移距道集获得单元,用于获得通过预设扫描信号扫描得到的近道共偏移距道集;
大地震动耦合响应谱获得单元,用于根据所述近道共偏移距道集进行统计子波计算,获得大地震动耦合响应谱;
预设扫描信号谱获得单元,用于获得所述预设扫描信号谱;
第一目标扫描信号获得单元,用于利用所述大地震动耦合响应谱对所述预设扫描信号谱进行能量补偿,并根据补偿后的频谱获得第一目标扫描信号。
在一个优选的实施例中,该装置还包括:
震源出产输出特性谱获得单元,用于获得震源出产输出特性谱;
第二目标扫描信号获得单元,用于利用所述震源出产输出特性谱和所述大地震动耦合响应谱对所述预设扫描信号谱进行能量补偿,并根据补偿后的频谱获得第二目标扫描信号。
在一个优选的实施例中,所述大地震动耦合响应谱获得单元包括:
统计子波获得模块,用于获得所述近道共偏移距道集中每个接收点的地震子波,把所有地震子波进行最小二乘法拟合获得统计子波;
傅里叶变换模块,用于将所述统计子波进行傅里叶变换得到频谱,所述频谱作为大地震动耦合响应谱。
由此可见,本申请实施例提供的技术方案在当前可控震源系统,利用预设扫描信号进行扫描,并获得相应的近道共偏移距道集,根据所述近道共偏移距道集进行统计子波计算,获得大地震动耦合响应谱。利用所述大地震动耦合响应谱对所述预设扫描信号进行近场输出能量补偿,并经过多次迭代获得第一目标扫描信号。后续可以利用近场输出能量补偿后得到的所述第一目标扫描信号进行扫描,即可得到宽频均衡输出的地震波。与现有技术相比,有效改善了震源输出响应的不确定性,使地震波的宽频均衡输出。
附图说明
为了更清楚地说明本申请实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本申请中记载的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。
图1是本申请一种获得宽频可控震源扫描信号的方法的第一实施例的流程图;
图2是利用预设扫描信号对某地表进行扫描得到的多个接收点近道子波频谱、各近道子波的拟合函数的频谱、以及该扫描信号谱的示意图;
图3是本申请一种获得宽频可控震源扫描信号的方法的第二实施例的流程图;
图4是利用本申请实施例提供技术方案获得的能量补偿后的第二目标扫描信号谱示意图;
图5是利用本申请实施例提供技术方案获得的能量补偿后的第二目标扫描信号对大地进行扫描的到的输出信号谱的示意图;
图6是本申请实施例提供的一种获得宽频可控震源扫描信号的装置的示意图;
图7是本申请实施例提供的大地震动耦合响应谱获得单元的模块示意图。
具体实施方式
为了使本技术领域的人员更好地理解本申请中的技术方案,下面将结合本申请实施例中的附图,对本申请实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本申请一部分实施例,而不是全部的实施例。基于本申请中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都应当属于本申请保护的范围。
下面以几个具体的例子详细说明本申请实施例的具体实现。
利用可控震源进行地震勘探时,一般可以在勘探区域布置二维或三维测线,使用可控震源激发地震波。可控震源所在点可以被称为炮点,也即激发点。沿着测线可以布置多个检波器来接收地震波信号,也即接收点。检波器或检波器组的数量可以根据需要进行设定,各检波器组之间排列方式可以不同,如中间放炮排列、端点放炮排列等,也可以根据需要进行设定。炮点和对应的检波器所在的检波点就可以形成一个地震道。检波器接收到地震波信号后,可以离散采样地震数据,并以数字形式记录在磁带上,作为该地震道的地震数据。
以下首先介绍本申请一种获得宽频可控震源扫描信号的方法的第一实施例。结合附图1,该实施例包括:
S110:获得通过预设扫描信号扫描得到的近道共偏移距道集。
这里的预设扫描信号是预先设置的扫描信号,且所述预设扫描信号在各频段是均衡输出的。但当可控震源系统利用所述预设扫描信号对大地进行扫描时,由于震源系统与大地耦合条件的影响,最终的输出并不像所述预设的扫描信号那样均衡输出,故需要对所述预设扫描信号进行近场输出能量补偿,也即对当前可控震源系统预设扫描信号进行因震源系统与大地耦合导致的输出能量损失的修正。
进一步的,一般共偏移距道集是指按照同一个偏移距,从不同共炮集或共道集数据抽取形成的道集,即该道集中的各接收点炮检距都相等;而这里的近道共偏移距道集是指所述偏移距为最小时的共偏移距道集。
进一步的,所述偏移距也即震源点与接收点之间的距离,所述最小的偏移距需要保证检波器在这一距离内能正常工作,且接收点接收到的地震信号为直达波。
S120:根据所述近道共偏移距道集进行统计子波计算,获得大地震动耦合响应谱。
所述根据所述近道共偏移距道集进行统计子波计算,获得大地震动耦合响应谱具体包括:
通过取所述道集中每个接收点与所述预设扫描信号的扫描长度相同的地震记录,所述地震数据作为该接收点的地震子波,把所有地震子波进行最小二乘法拟合获得统计子波。
将所述统计子波进行傅里叶变换得到频谱,所述频谱作为大地震动耦合响应谱。
如图2所示是用预设扫描信号对某地表进行扫描得到的多个接收点近道子波频谱、各近道子波的拟合函数的频谱、以及该扫描信号谱的示意图。所述近道子波的拟合函数的频谱也即大地震动耦合响应谱,图中包括预设扫描信号谱1101、多个接收点近道子波频谱(其中包括某一接收点近道子波频谱1102)以及大地震动耦合响应谱1103;其中横坐标代表频率,单位(HZ),纵坐标A代表能量。从图中可以看出,虽然采用了均衡的预设扫描进行扫描,但得到的大地震动耦合响应谱有了一定的能量损失,且各子波频谱输出不均衡。
S130:获得预设扫描信号谱。
这里的预设扫描信号谱是指利用预设扫描信号进行傅里叶变换得到的相应频谱。
S140:利用所述大地震动耦合响应谱对所述预设扫描信号谱进行能量补偿,并根据补偿后的频谱获得第一目标扫描信号。
利用所述预设扫描信号扫描下获得的大地震动耦合响应谱对所述预设扫描信号进行近场输出能量补偿,,并根据补偿后的频谱获得第一目标扫描信号具体的是:
S141:根据所述大地震动耦合响应谱以第一预设迭代补偿公式对所述预设扫描信号谱进行能量补偿;
S142:获得当前频谱,将所述当前频谱进行傅里叶反变换获得当前扫描信号;
S143:判断根据所述当前扫描信号扫描大地后得到的地震波是否为宽频均衡输出;
这里的地震波为宽频均衡输出具体包括:这里的宽频是大于5个倍频程(1.5HZ~3HZ为第1个倍频程,3HZ~6HZ为第2个倍频程,6HZ~12HZ为第3个倍频程,12HZ~24HZ为第4个倍频程,24HZ~48HZ为第5个倍频程),且当前接收到得到地震波在整个频谱上能量基本都是均衡分布的。
S144:当所述判断为是时,所述当前扫描信号作为第一目标扫描信号。
所述当前扫描信号作为第一目标扫描信号,所述当前频谱也即第一目标扫描信号谱。
S145:当所述判断为否时,根据所述大地震动耦合响应谱以第一预设迭代补偿公式对当前频谱进行能量补偿,并进入步骤S142。
所述第一预设迭代补偿公式如下:
S n + 1 ( f ) = S n ( f ) E ( f )
其中,Sn(f)为第n次补偿迭代后的频谱;n为迭代的次数(n=0、1、2…),当n=0时,Sn(f)代表预设扫描信号谱;Sn+1(f)为第n+1次补偿迭代后的频谱;E(f)为大地震动耦合响应谱。
由此可见,本申请第一方法实施例提供的技术方案在当前可控震源系统,利用预设扫描信号进行扫描,并获得相应的近道共偏移距道集,根据所述近道共偏移距道集进行统计子波计算,获得大地震动耦合响应谱。利用所述大地震动耦合响应谱对所述预设扫描信号进行近场输出能量补偿,并经过多次迭代获得第一目标扫描信号。后续可以利用近场输出能量补偿后得到的所述第一目标扫描信号进行扫描,即可得到宽频均衡输出的地震波。与现有技术相比,有效改善了震源输出响应的不确定性,使地震波的宽频均衡输出。
以下介绍本申请一种获得宽频可控震源扫描信号的方法的第二实施例。结合附图3,该实施例包括:
S210:获得通过预设扫描信号扫描得到的近道共偏移距道集。
S220:根据所述近道共偏移距道集进行统计子波计算,获得大地震动耦合响应谱。
S230:获得预设扫描信号谱。
S240:获得震源出产输出特性谱。
这里的震源出产输出谱是指当前可控震源系统可正常输出能量的特性曲线,进行傅里叶变换得到的相应频谱。
在实际勘探过程中,当震源系统出产后,其结构是固定,因此可以稳定均衡输出的频段是固定的,比如10HZ以上的频段都是均衡输出的,对于一些特殊的目标,可能需要其他频段的实现均衡输出,比如3HZ到10HZ频段也需要均衡输出时,震源的能量输出是受结构控制不能满振幅,并且是非线性的,那么就需要对预设扫描信号进行符合震源出产输出特性的振幅的补偿。
在本申请的具体实施例中,可以采用现有的地震数据采集系统获得所述震源出产输出特性谱。
S250:利用所述震源出产输出特性谱和所述大地震动耦合响应谱对所述预设扫描信号谱进行能量补偿,并根据补偿后的频谱获得第二目标扫描信号。
获得所述第二目标扫描信号可以包括:
S251:根据所述大地震动耦合响应谱以第二预设迭代补偿公式对所述预设扫描信号谱进行能量补偿;
S252:获得当前频谱,将所述当前频谱进行傅里叶反变换获得当前扫描信号;
S253:判断根据所述当前扫描信号扫描大地后得到的地震波是否为宽频均衡输出;
这里的地震波为宽频均衡输出具体包括:这里的宽频是大于5个倍频程(1.5HZ~3HZ为第1个倍频程,3HZ~6HZ为第2个倍频程,6HZ~12HZ为第3个倍频程,12HZ~24HZ为第4个倍频程,24HZ~48HZ为第5个倍频程),且当前接收到得到地震波在整个频谱上能量基本都是均衡分布的。
S254:当所述判断为是时,所述当前扫描信号作为第二目标扫描信号;
所述当前扫描信号作为第二目标扫描信号,所述当前频谱也即第二目标扫描信号谱。
如图4所示的利用本申请实施例提供技术方案获得的能量补偿后的第二目标扫描信号谱示意图。其中横坐标代表频率,单位(HZ),纵坐标A代表能量。
利用经过能量补偿后得到的所述第二目标扫描信号谱傅里叶反变换得到第二目标扫描信号进行扫描,即可得到宽频均衡输出的地震波。如图5所示,其中横坐标代表频率,单位(HZ),纵坐标A代表能量。
S255:当所述判断为否时,根据所述大地震动耦合响应谱以第二预设迭代补偿公式对当前频谱进行能量补偿,并进入步骤S252。
所述第二预设迭代补偿公式如下:
S m + 1 ( f ) = S m ( f ) E ( f ) × G ( f )
其中,Sm(f)为第m次补偿迭代后的频谱;m为迭代的次数(m=0、1、2…),当m=0时,Sm(f)代表预设扫描信号谱;Sm+1(f)为第m+1次补偿迭代后的频谱;E(f)为大地震动耦合响应谱;G(f)震源出产输出特性谱。
由此可见,本申请第二方法实施例提供的技术方案在当前可控震源系统,利用预设扫描信号进行扫描,并获得相应的近道共偏移距道集,根据所述近道共偏移距道集进行统计子波计算,获得大地震动耦合响应谱。利用所述大地震动耦合响应谱对所述预设扫描信号进行近场输出能量补偿、以及所述震源出产输出特性谱对预设扫描信号进行振幅的补偿,并经过多次迭代后获得第二目标扫描信号。后续可以所述第二目标扫描信号对大地进行扫描,即可得到宽频均衡输出的地震波,解决了由于震源系统结构、以及可控震源系统与大地震动耦合导致的能量损失,为准确的确定多孔介质的物性参数、流体含量、边界以及孔隙度等地质信息提供了保证。与现有技术相比,有效改善了震源输出响应的不确定性,使地震波的宽频均衡输出。
本申请另一方面还提供一种获得宽频可控震源扫描信号的装置,结合附图6,该装置300包括:近道共偏移距道集获得单310,大地震动耦合响应谱获得单元320,预设扫描信号谱获得单330,第一目标扫描信号谱获得单元340,第一目标扫描信号获得单元350。
其中,近道共偏移距道集获得单元310用于获得通过预设扫描信号扫描得到的近道共偏移距道集;大地震动耦合响应谱获得单元320用于根据所述近道共偏移距道集进行统计子波计算,获得大地震动耦合响应谱;预设扫描信号谱获得单元330用于获得预设扫描信号谱;第一目标扫描信号获得单元340用于利用所述大地震动耦合响应谱对所述预设扫描信号谱进行能量补偿,并根据补偿后的频谱获得第一目标扫描信号。
在一个优选的实施例中,该装置300还包括:
震源出产输出特性谱获得单元,用于获得震源出产输出特性谱;
第二目标扫描信号获得单元,用于利用所述震源出产输出特性谱和所述大地震动耦合响应谱对所述预设扫描信号谱进行能量补偿,并根据补偿后的频谱获得第二目标扫描信号。
在一个优选的实施例中,结合附图7,所述大地震动耦合响应谱获得单元320包括:
统计子波获得模块321,用于获得所述近道共偏移距道集中每个接收点的地震子波,把所有地震子波进行最小二乘法拟合获得统计子波;
傅里叶变换模块322,用于将所述统计子波进行傅里叶变换得到频谱,所述频谱作为大地震动耦合响应谱。
由此可见,本申请一种获得宽频可控震源扫描信号的方法和装置提供的技术方案在当前可控震源系统,利用预设扫描信号进行扫描,并获得相应的近道共偏移距道集,根据所述近道共偏移距道集进行统计子波计算,获得大地震动耦合响应谱。利用所述大地震动耦合响应谱对所述预设扫描信号进行近场输出能量补偿、以及所述震源出产输出特性谱对预设扫描信号进行振幅的补偿,并经过多次迭代后获得第二目标扫描信号。后续可以所述第二目标扫描信号对大地进行扫描,即可得到宽频均衡输出的地震波,解决了由于震源系统结构、以及可控震源系统与大地震动耦合导致的能量损失,为准确的确定多孔介质的物性参数、流体含量、边界以及孔隙度等地质信息提供了保证。与现有技术相比,有效改善了震源输出响应的不确定性,使地震波的宽频均衡输出。
本说明书中的各个实施例均采用递进的方式描述,每个实施例重点说明的都是与其他实施例的不同之处,各个实施例之间相同相似的部分互相参见即可。尤其,对于系统实施例而言,由于其基本相似于方法实施例,所以描述的比较简单,相关之处参见方法实施例的部分说明即可。
虽然通过实施例描绘了本申请,本领域普通技术人员知道,本申请有许多变形和变化而不脱离本申请的精神,希望所附的权利要求包括这些变形和变化而不脱离本申请的精神。

Claims (10)

1.一种获得宽频可控震源扫描信号的方法,其特征在于,该方法包括:
获得通过预设扫描信号扫描得到的近道共偏移距道集;
根据所述近道共偏移距道集进行统计子波计算,获得大地震动耦合响应谱;
获得预设扫描信号谱;
利用所述大地震动耦合响应谱对所述预设扫描信号谱进行能量补偿,并根据补偿后的频谱获得第一目标扫描信号。
2.根据权利要求1所述的获得宽频可控震源扫描信号的方法,其特征在于,该方法还包括:
获得震源出产输出特性谱;
利用所述震源出产输出特性谱和所述大地震动耦合响应谱对所述预设扫描信号谱进行能量补偿,并根据补偿后的频谱获得第二目标扫描信号。
3.根据权利要求1或2所述的获得宽频可控震源扫描信号的方法,其特征在于,所述根据所述近道共偏移距道集进行统计子波计算,获得大地震动耦合响应谱包括:
获得所述近道共偏移距道集中每个接收点的地震子波,把所有地震子波进行最小二乘法拟合获得统计子波;
将所述统计子波进行傅里叶变换得到频谱,所述频谱作为大地震动耦合响应谱。
4.根据权利要求1所述的获得宽频可控震源扫描信号的方法,其特征在于,所述利用所述大地震动耦合响应谱对所述预设扫描信号谱进行能量补偿,并根据补偿后的频谱获得第一目标扫描信号包括:
S11:根据所述大地震动耦合响应谱以第一预设迭代补偿公式对所述预设扫描信号谱进行能量补偿;
S12:获得当前频谱,将所述当前频谱进行傅里叶反变换获得当前扫描信号;
S13:判断根据所述当前扫描信号扫描大地后得到的地震波是否为宽频均衡输出;
S14:当所述判断为是时,所述当前扫描信号作为第一目标扫描信号;
S15:当所述判断为否时,根据所述大地震动耦合响应谱以第一预设迭代补偿公式对当前频谱进行能量补偿,并进入步骤S12。
5.根据权利要求4所述的获得宽频可控震源扫描信号的方法,其特征在于,所述第一预设迭代补偿公式如下:
S n + 1 ( f ) = S n ( f ) E ( f )
其中,Sn(f)为第n次补偿迭代后的频谱;n为迭代的次数(n=0、1、2…),当n=0时,Sn(f)代表预设扫描信号谱;Sn+1(f)为第n+1次补偿迭代后的频谱;E(f)为大地震动耦合响应谱。
6.根据权利要求2所述的获得宽频可控震源扫描信号的方法,其特征在于,所述利用所述震源出产输出特性谱和所述大地震动耦合响应谱对所述预设扫描信号谱进行能量补偿,并根据补偿后的频谱获得第二目标扫描信号包括:
S21:根据所述大地震动耦合响应谱以第二预设迭代补偿公式对所述预设扫描信号谱进行能量补偿;
S22:获得当前频谱,将所述当前频谱进行傅里叶反变换获得当前扫描信号;
S23:判断根据所述当前扫描信号扫描大地后得到的地震波是否为宽频均衡输出;
S24:当所述判断为是时,所述当前扫描信号作为第二目标扫描信号;
S25:当所述判断为否时,根据所述大地震动耦合响应谱以第二预设迭代补偿公式对当前频谱进行能量补偿,并进入步骤S22。
7.根据权利要求6所述的获得宽频可控震源扫描信号的方法,其特征在于,所述第二预设迭代补偿公式如下:
S m + 1 ( f ) = S m ( f ) E ( f ) × G ( f )
其中,Sm(f)为第m次补偿迭代后的频谱;m为迭代的次数(m=0、1、2…),当m=0时,Sm(f)代表预设扫描信号谱;Sm+1(f)为第m+1次补偿迭代后的频谱;E(f)为大地震动耦合响应谱;G(f)震源出产输出特性谱。
8.一种获得宽频可控震源扫描信号的装置,其特征在于,该装置包括:
近道共偏移距道集获得单元,用于获得通过预设扫描信号扫描得到的近道共偏移距道集;
大地震动耦合响应谱获得单元,用于根据所述近道共偏移距道集进行统计子波计算,获得大地震动耦合响应谱;
预设扫描信号谱获得单元,用于获得预设扫描信号谱;
第一目标扫描信号获得单元,用于利用所述大地震动耦合响应谱对所述预设扫描信号谱进行能量补偿,并根据补偿后的频谱获得第一目标扫描信号。
9.根据权利要求8所述的获得宽频可控震源扫描信号的装置,其特征在于,该装置还包括:
震源出产输出特性谱获得单元,用于获得震源出产输出特性谱;
第二目标扫描信号获得单元,用于利用所述震源出产输出特性谱和所述大地震动耦合响应谱对所述预设扫描信号谱进行能量补偿,并根据补偿后的频谱获得第二目标扫描信号。
10.根据权利要求8或9所述的获得宽频可控震源扫描信号的装置,其特征在于,所述大地震动耦合响应谱获得单元包括:
统计子波获得模块,用于获得所述近道共偏移距道集中每个接收点的地震子波,把所有地震子波进行最小二乘法拟合获得统计子波;
傅里叶变换模块,用于将所述统计子波进行傅里叶变换得到频谱,所述频谱作为大地震动耦合响应谱。
CN201510108954.9A 2015-03-12 2015-03-12 一种获得宽频可控震源扫描信号的方法和装置 Pending CN104793237A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510108954.9A CN104793237A (zh) 2015-03-12 2015-03-12 一种获得宽频可控震源扫描信号的方法和装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510108954.9A CN104793237A (zh) 2015-03-12 2015-03-12 一种获得宽频可控震源扫描信号的方法和装置

Publications (1)

Publication Number Publication Date
CN104793237A true CN104793237A (zh) 2015-07-22

Family

ID=53558199

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510108954.9A Pending CN104793237A (zh) 2015-03-12 2015-03-12 一种获得宽频可控震源扫描信号的方法和装置

Country Status (1)

Country Link
CN (1) CN104793237A (zh)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109782335A (zh) * 2017-11-10 2019-05-21 中石化石油工程技术服务有限公司 低畸变宽频扫描信号设计方法
CN112526600A (zh) * 2019-09-19 2021-03-19 中国石油化工股份有限公司 一种扫频超声波激发信号的振幅补偿方法及系统
CN112666597A (zh) * 2019-10-16 2021-04-16 中国石油天然气集团有限公司 可控震源扫描方法和装置
CN113127799A (zh) * 2020-01-15 2021-07-16 中国石油天然气集团有限公司 一种扫描信号校正方法及装置
CN113138417A (zh) * 2020-01-20 2021-07-20 中国石油天然气集团有限公司 一种基于力信号反馈的扫描信号自动校正方法及装置

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20040044478A1 (en) * 2002-09-04 2004-03-04 Rainer Moerig Vibrator sweep shaping method
CN103777240A (zh) * 2012-10-25 2014-05-07 中国石油化工股份有限公司 一种可控震源分频扫描频带补偿方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20040044478A1 (en) * 2002-09-04 2004-03-04 Rainer Moerig Vibrator sweep shaping method
CN103777240A (zh) * 2012-10-25 2014-05-07 中国石油化工股份有限公司 一种可控震源分频扫描频带补偿方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
CLAUDIO BAGAINI ET AL.: "land seismic techniques for high-quality data", 《OILFIELD REVIEW》 *
LUO FEI ET AL.: "Design broadband sweep signal using the statistical near field wavelet", 《SOCIETY OF EXPLORATION GEOPHYSICISTS》 *
PETER MAXWELL ET AL.: "Extending low frequency bandwidth using pseudorandom sweeps", 《SEG DENVER 2010 ANNUAL MEETING》 *
陶知非,等: "可控震源低频信号激发技术的最新进展", 《物探装备》 *

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109782335A (zh) * 2017-11-10 2019-05-21 中石化石油工程技术服务有限公司 低畸变宽频扫描信号设计方法
CN112526600A (zh) * 2019-09-19 2021-03-19 中国石油化工股份有限公司 一种扫频超声波激发信号的振幅补偿方法及系统
CN112526600B (zh) * 2019-09-19 2023-12-12 中国石油化工股份有限公司 一种扫频超声波激发信号的振幅补偿方法及系统
CN112666597A (zh) * 2019-10-16 2021-04-16 中国石油天然气集团有限公司 可控震源扫描方法和装置
CN112666597B (zh) * 2019-10-16 2024-03-26 中国石油天然气集团有限公司 可控震源扫描方法和装置
CN113127799A (zh) * 2020-01-15 2021-07-16 中国石油天然气集团有限公司 一种扫描信号校正方法及装置
CN113127799B (zh) * 2020-01-15 2023-12-01 中国石油天然气集团有限公司 一种扫描信号校正方法及装置
CN113138417A (zh) * 2020-01-20 2021-07-20 中国石油天然气集团有限公司 一种基于力信号反馈的扫描信号自动校正方法及装置

Similar Documents

Publication Publication Date Title
Foti et al. Guidelines for the good practice of surface wave analysis: a product of the InterPACIFIC project
US8483009B2 (en) Characterizing spatial variability of surface waves in seismic processing
Cho et al. A new method to determine phase velocities of Rayleigh waves from microseisms
US9588241B2 (en) Frequency-varying filtering of simultaneous source seismic data
CN102305941B (zh) 由叠前时间偏移直接扫描确定地层叠加品质因子方法
KR20200014387A (ko) 지하 구조물의 검출
CN104793237A (zh) 一种获得宽频可控震源扫描信号的方法和装置
CN102109615B (zh) 地震资料处理真地表条件下的振幅补偿方法和装置
US20080170468A1 (en) Method for Processing at Least Two Sets of Seismic Data
CN106896409B (zh) 一种基于波动方程边值反演的变深度缆鬼波压制方法
CN104533396A (zh) 一种远探测声波的处理方法
WO2011034870A1 (en) Time-lapse seismic comparisons using pre-stack imaging and complex wave field comparisons to improve accuracy and detail
CN105116443A (zh) 一种低频信号的能量补偿方法及装置
CN101852864B (zh) 一种利用地表一致性统计频谱分析技术处理海量地震数据的方法
Gouveia et al. Deeper Vs profile constraining the dispersion curve with the ellipticity curve: A case study in Lower Tagus Valley, Portugal
CN108490488A (zh) 一种基于时窗能量比的槽波地震散射成像方法及装置
CN115877451A (zh) 一种多激发点瞬态面波勘探方法、设备及存储介质
CN115373023A (zh) 一种基于地震反射和车辆噪声的联合探测方法
CN108594301B (zh) 一种具有差异特征的地震数据融合的方法及处理终端
CN108121010B (zh) 基于孔巷联合的井下独头巷道槽波超前探测方法及系统
CN113514889B (zh) 一种提升海洋深反射地震数据中低频信号能量的处理方法
Rutenko et al. A method for estimating the characteristics of acoustic pulses recorded on the sakhalin shelf for multivariate analysis of their effect on the behavior of gray whales
US5991235A (en) Seismic data imaging system
CN111691876A (zh) 一种利用声波测井对邻井成像的方法、装置及存储介质
CN114721044B (zh) 一种多频率接收函数和振幅比联合反演地壳结构的方法及系统

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
EXSB Decision made by sipo to initiate substantive examination
SE01 Entry into force of request for substantive examination
WD01 Invention patent application deemed withdrawn after publication

Application publication date: 20150722

WD01 Invention patent application deemed withdrawn after publication