CN108318919A - 一种动态参数的可控震源非线性扫描信号设计系统及方法 - Google Patents
一种动态参数的可控震源非线性扫描信号设计系统及方法 Download PDFInfo
- Publication number
- CN108318919A CN108318919A CN201810114737.4A CN201810114737A CN108318919A CN 108318919 A CN108318919 A CN 108318919A CN 201810114737 A CN201810114737 A CN 201810114737A CN 108318919 A CN108318919 A CN 108318919A
- Authority
- CN
- China
- Prior art keywords
- dynamic parameter
- scanning signal
- frequency
- signal
- controlled source
- 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.)
- Granted
Links
- 238000013461 design Methods 0.000 title claims abstract description 44
- 238000000034 method Methods 0.000 title claims abstract description 15
- 238000004088 simulation Methods 0.000 claims abstract description 56
- 238000012360 testing method Methods 0.000 claims abstract description 30
- 238000004458 analytical method Methods 0.000 claims abstract description 27
- 230000000694 effects Effects 0.000 claims description 14
- 238000009795 derivation Methods 0.000 claims description 9
- 230000003247 decreasing effect Effects 0.000 claims description 3
- 239000004615 ingredient Substances 0.000 claims description 3
- 238000012545 processing Methods 0.000 claims description 2
- 238000001228 spectrum Methods 0.000 description 6
- 238000010586 diagram Methods 0.000 description 4
- 230000011514 reflex Effects 0.000 description 4
- 238000004364 calculation method Methods 0.000 description 3
- 238000006243 chemical reaction Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 239000000284 extract Substances 0.000 description 2
- 230000006870 function Effects 0.000 description 2
- 230000002787 reinforcement Effects 0.000 description 2
- 238000010408 sweeping Methods 0.000 description 2
- 238000012546 transfer Methods 0.000 description 2
- 238000007792 addition Methods 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 230000002708 enhancing effect Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
- 230000000644 propagated effect Effects 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
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/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/288—Event detection in seismic signals, e.g. microseismics
-
- 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/30—Analysis
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- Acoustics & Sound (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Business, Economics & Management (AREA)
- Emergency Management (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明公开了一种动态参数的可控震源非线性扫描信号设计的系统及方法,属于石油地球物理勘探技术领域;本发明设计了两种动态参数的可控震源非线性扫描信号,两种信号可以基于不同的地质目标,通过调节动态参数实现地震波优势频带的补偿。采用衰减介质模型对两种非线性扫描信号进行可控震源黏声正演模拟测试和分析,测试不同参数非线性扫描信号对谐波干扰的影响。通过测试发现,本发明基于不同的地质目标,通过调节动态参数实现优势频带补偿,在衰减介质中,采用本发明提出的两种非线性扫描信号通过调整动态参数可以较好地压制谐波干扰。
Description
技术领域
本发明属于石油勘探领域,具体涉及一种动态参数的可控震源非线性扫描信号设计系统及方法。
背景技术
可控震源采集信号能量受到可控震源车吨位的限制,为此在输出能量一定的情况下,如何设计高效且易于施工的扫描信号至关重要。目前,常见的可控震源信号大致分为:线性扫描信号、非线性扫描信号、伪随机扫描信号。其中,线性扫描信号利于控制且频带宽度较大,是野外普遍使用的扫描信号;非线性扫描信号能够整形振幅谱,达到补偿高频的目的,对高频地震勘探具有重要意义;伪随机扫描性信号的子波旁瓣最小,有利于提高采集数据的分辨率。实际勘探中,还应该通过选择合适可控震源扫描信号参数,以满足不同探区对采集资料品质的需求。
发明内容
针对现有技术中存在的上述技术问题,本发明提出了一种动态参数的可控震源非线性扫描信号设计系统及方法,设计合理,克服了现有技术的不足,具有良好的效果。
为了实现上述目的,本发明采用如下技术方案:
一种动态参数的可控震源非线性扫描信号设计系统,包括二次多项式动态参数扫描信号设计模块、指数动态参数扫描信号设计模块、可控震源黏声正演模块、二次多项式动态参数扫描信号正演模拟与分析模块、指数动态参数扫描信号正演模拟与分析模块以及动态参数扫描信号对谐波干扰影响模块;
二次多项式动态参数扫描信号设计模块,被配置为用于对二次多项式动态参数扫描信号进行设计;
指数动态参数扫描信号设计模块,被配置为用于对指数动态参数扫描信号进行设计;
可控震源黏声正演模块,被配置为用于对线性扫描信号进行可控震源黏声正演模拟;
二次多项式动态参数扫描信号正演模拟与分析模块,被配置为用于对二次多项式动态参数扫描信号的可控震源进行正演模拟测试和效果分析;
指数动态参数扫描信号正演模拟与分析模块,被配置为用于指数动态参数扫描信号的可控震源进行正演模拟测试和频带补偿效果分析;
动态参数扫描信号对谐波干扰影响模块,被配置为用于在指数动态参数扫描信号和二次多项式动态参数扫描信号正演模拟中加入谐波,测试指数动态参数扫描信号和二次多项式动态参数扫描信号对谐波干扰的影响。
此外,本发明还提到一种动态参数的可控震源非线性扫描信号设计方法,该方法采用如上所述的一种动态参数的可控震源非线性扫描信号设计系统,包括如下步骤:
步骤1:通过二次多项式动态参数扫描信号设计模块,设计二次多项式动态参数扫描信号,其表达式如公式(5)所示;
a≥1或a≤-1 (5);
其中,f1和f2分别表示起始频率和终止频率,T为扫描周期,a为动态参数;
步骤2:通过指数动态参数扫描信号设计模块,设计指数动态参数扫描信号,其表达式如公式(6)所示;
步骤3:通过可控震源黏声正演模块,对线性扫描信号进行可控震源黏声正演模拟,二维黏滞声波方程为:
其中,P表示为压力场,τσl和τεl分别表示为弛豫时间常数,GR表示弛豫模量;L代表标准线性体的个数,rl表示记忆变量,vx和vz分别为水平分量和垂直分量的速度;
步骤4:通过二次多项式动态参数扫描信号正演模拟与分析模块,测试二次多项式动态参数扫描信号对地震波高频或者低频的补偿效果,测试参数为:起始频率为5Hz,终止频率为80Hz,扫描长度为12s,动态参数a分别取a=1,a=2,a=-1,a=-2;
步骤5:通过指数动态参数扫描信号正演模拟与分析模块,测试指数动态参数扫描信号对地震波不同频率成分的补偿效果,测试参数为:起始频率为5Hz,终止频率为80Hz,扫描长度为12s,动态参数a分别取a=0,a=0.2,a=0.5,a=0.7,a=1;
步骤6:通过动态参数扫描信号对谐波干扰影响模块,在正演模拟中加入二阶谐波和三阶谐波,测试动态参数扫描信号对谐波干扰的影响。
优选地,在步骤1中,二次多项式动态参数扫描信号公式(5)的推导过程如下:
首先,设计瞬时频率的两个二次多项式:
分别对上式进行求导:
其中,F′1_1(t)的斜率为F1'_2(t)的斜率为由此可知,瞬时频率 F1_1(t)的增加速度越来越快,表示低频补偿;瞬时频率F1_2(t)的增加速度越来越慢,表示高频补偿;根据公式(1)和(2),引入动态参数a,实现二次多项式动态参数扫描信号设计,通过调节a,可实现对高频或者低频进行动态补偿;
二次多项式动态参数扫描信号F1(t)满足如下几个条件:
(5)F1(t)max=F1(T)=f2;
(6)F1(t)min=F1(0)=f1;
(7)当a>0时,F′1_1(t)单调递增;
(8)当a<0时,F′1_2(t)单调递减;
根据以上四个条件并结合公式(1)和(2),可推导出二次多项式动态参数扫描信号公式 (5)。
优选地,在步骤3中,二维黏滞声波方程中,为了方便求解,将二阶导数进行降阶处理,转化为一阶函数求导,从而提高了正演模拟的精度,且引入的辅助微分算子如下所示:
本发明所带来的有益技术效果:
本发明设计了两种动态参数的可控震源非线性扫描信号:二次多项式动态参数扫描信号和指数动态参数扫描信号;采用可控震源黏声正演模块进行测试与分析,通过调整参数,二次多项式动态参数扫描信号可实现高频或者低频补偿,而指数动态参数扫描信号不仅可以补偿高频或者低频,也可以基于地质目标实现中间频带的加强;当补偿高频成分时,能够在一定程度上补偿在地层衰减中的高频能量,可拓宽频带;当补偿低频成分时,分辨率降低,但因高频成分能量少,造成的衰减较少,因此深部反射能量相对更强,有利于在衰减比较严重的区域进行深层勘探,且对低频进行补偿后,地震记录的信噪比得到了一定的提升;低频补偿的非线性扫描信号可以很大程度上压制谐波噪音,提高地震资料的分辨率,在实际工区中,不同地区的地质构造和地质目标不同,造成了非线性扫描信号的需求不同,因此本论文提出的两种动态参数的非线性扫描信号具有广阔的应用前景。
附图说明
图1为本发明一种动态参数的可控震源非线性扫描信号设计方法的流程图;
图2为二次多项式动态参数扫描信号与线性升频扫描信号瞬时频率对比图;
图3为指数动态参数扫描信号与线性升频扫描信号瞬时频率对比图;
图4为线性升频扫描信号及频谱图,(a)为线性升频扫描信号,(b)为频谱图;
图5为正演模拟的速度模型与品质因子模型,(a)为速度模型,(b)为品质因子模型;
图6为线性扫描信号可控震源正演模拟炮记录;(a)为声波介质正演模拟得到的炮记录, (b)为黏声介质正演模拟得到的炮记录。
图7为从炮记录图(6)中抽取的第151道波形图,(a)全局图,(b)局部放大图;
图8为从炮记录图(6)中抽取的第151道频谱图;
图9为二次多项式动态参数扫描信号,(a)参数a=1的扫描信号,(b)参数a=-1的扫描信号,(c)参数a=2的扫描信号,(d)参数a=-2的扫描信号;
图10为二次多项式动态参数扫描信号的频谱图,(a)参数a=1的扫描信号频谱图,(b) 参数a=-1的扫描信号频谱图,(c)参数a=2的扫描信号频谱图,(d)参数a=-2的扫描信号频谱图;
图11为二次多项式动态参数扫描信号可控震源正演模拟炮记录,(a)参数a=1,(b)参数a=-1,(c)参数a=2,(d)参数a=-2;
图12为从图11炮记录中抽取的第120道单道记录,(a)参数a=1,(b)参数a=-1,(c) 参数a=2,(d)参数a=-2;
图13为指数动态参数扫描信号可控震源正演模拟炮记录,(a)参数a=0,(b)参数a=0.2, (c)参数a=0.5,(d)参数a=0.7,(e)参数a=1;
图14为从图13炮记录中抽取的第120道单道记录,(a)参数a=0,(b)参数a=0.2,(c) 参数a=0.5,(d)参数a=0.7,(e)参数a=1;
图15为指数动态参数扫描信号得到的包含谐波干扰的可控震源炮记录(a1-e1)与谐波干扰(a2-e2),(a1,a2)中参数a=0,(b1,b2)中参数a=0.2,(c1,c2)中参数a=0.5,(d1,d2) 中参数a=0.7,(e1,e2)中参数a=1;
图16为从图15(a2-e2)炮记录中抽取的第150道的单道记录,(a)参数a=0,(b)参数a=0.2,(c)参数a=0.5,(d)参数a=0.7,(e)参数a=1;
图17为本发明的实施方式中一种动态参数的可控震源非线性扫描信号设计系统的结构示意图。
具体实施方式
下面结合附图以及具体实施方式对本发明作进一步详细说明:
实施例1:
如图17所示,一种基于动态参数的可控震源非线性扫描信号设计系统,包括二次多项式动态参数扫描信号设计模块、指数动态参数扫描信号设计模块、可控震源黏声正演模块、二次多项式动态参数扫描信号正演模拟与分析模块、指数动态参数扫描信号正演模拟与分析模块以及动态参数扫描信号对谐波干扰影响模块;
二次多项式动态参数扫描信号设计模块,被配置为用于对二次多项式动态参数扫描信号进行设计;
指数动态参数扫描信号设计模块,被配置为用于对指数动态参数扫描信号进行设计;
可控震源黏声正演模块,被配置为用于对线性扫描信号进行可控震源黏声正演模拟;
二次多项式动态参数扫描信号正演模拟与分析模块,被配置为用于对二次多项式动态参数扫描信号的可控震源进行正演模拟测试和效果分析;
指数动态参数扫描信号正演模拟与分析模块,被配置为用于指数动态参数扫描信号的可控震源进行正演模拟测试和频带补偿效果分析;
动态参数扫描信号对谐波干扰影响模块,被配置为用于在指数动态参数扫描信号和二次多项式动态参数扫描信号正演模拟中加入谐波,测试指数动态参数扫描信号和二次多项式动态参数扫描信号对谐波干扰的影响。
实施例2
在上述实施例的基础上,本发明还提到一种动态参数的可控震源非线性扫描信号设计方法,其流程如图1所示,具体包括如下步骤:
步骤1:通过二次多项式动态参数扫描信号设计模块,设计含有动态参数的二次多项式非线性扫描信号,其表达式如公式(5)所示;
a≥1或a≤-1 (5);
其中,f1和f2分别表示起始频率和终止频率,T为扫描周期,a为动态参数;
二次多项式动态参数扫描信号公式(5)的推导过程如下:
首先,设计瞬时频率的两个二次多项式:
分别对上式进行求导:
其中,F1'_1(t)的斜率为F1'_2(t)的斜率为由此可知,瞬时频率 F1_1(t)的增加速度越来越快,表示低频补偿;瞬时频率F1_2(t)的增加速度越来越慢,表示高频补偿;根据公式(1)和(2),引入动态参数a,实现二次多项式动态参数扫描信号设计,通过调节a,可实现对可控震源产生的地震波的低频或高频进行动态补偿;
二次多项式动态参数扫描信号F1(t)满足如下几个条件:
(9)F1(t)max=F1(T)=f2;
(10)F1(t)min=F1(0)=f1;
(11)当a>0时,F′1_1(t)单调递增;
(12)当a<0时,F′1_2(t)单调递减;
根据以上四个条件并结合公式(1)和(2),可推导出二次多项式动态参数扫描信号公式 (5)。
步骤2:通过指数动态参数扫描信号设计模块,设计含有动态参数的指数非线性扫描信号,其表达式如公式(6)所示;
步骤3:通过可控震源黏声正演模块,对线性扫描信号进行可控震源黏声正演模拟,二维黏滞声波方程为:
其中,P表示为压力场,τσl和τεl分别表示为弛豫时间常数,GR表示弛豫模量;L代表标准线性体的个数,rl表示记忆变量,vx和vz分别为水平分量和垂直分量的速度;
在步骤3中,二维黏滞声波方程中,为了方便求解,将二阶导数进行降阶处理,转化为一阶函数求导,从而提高了正演模拟的精度,且引入的辅助微分算子如下所示:
步骤4:通过二次多项式动态参数扫描信号正演模拟与分析模块,测试二次多项式动态参数扫描信号对地震波高频或者低频的补偿效果,测试参数为:起始频率为5Hz,终止频率为80Hz,扫描长度为12s,动态参数a分别取a=1,a=2,a=-1,a=-2;
步骤5:通过指数动态参数扫描信号正演模拟与分析模块,测试指数动态参数扫描信号对地震波不同频率成分的补偿效果,测试参数为:起始频率为5Hz,终止频率为80Hz,扫描长度为12s,动态参数a分别取a=0,a=0.2,a=0.5,a=0.7,a=1;
步骤6:通过动态参数扫描信号对谐波干扰影响模块,在正演模拟中加入二阶谐波和三阶谐波,测试动态参数扫描信号对谐波干扰的影响。
本发明设计的二次多项式动态参数扫描信号,通过调整参数的值,可实现对地震波的低频或高频进行补偿。测试参数如下所示:起始频率为5Hz,终止频率为80Hz,扫描长度为 12s。动态参数分别取a=1;a=-1;a=-2;a=2,不同动态参数的扫描信号和线性信号的瞬时频率随时间变化的关系如图2所示。当a=1时,低频区域瞬时频率增加较快,高频区域频率增加较慢,对高频成分进行大量补偿;当a=2时,瞬时频率介于a=1与线性升频扫描信号之间,对高频成分进行少量补偿;当a=-1时,低频区域瞬时频率增加较慢,高频区域频率增加较快,对低频成分进行大量补偿;当a=-2时,瞬时频率介于a=-1与线性升频扫描信号之间,对低频成分进行少量补偿。为了保证信号为升频信号,需要满足a≥1或a≤-1。
本发明设计的指数动态参数的扫描信号,通过调整参数的值,可实现对地震波的低频、中频和高频的补偿。测试参数如下所示:起始频率为5Hz,终止频率为80Hz,扫描长度为12 s。动态参数分别取a=0,a=0.2,a=0.5,a=0.7,a=1,不同动态参数的扫描信号和线性信号的瞬时频率随时间变化的关系如图3所示;分别取a=0,a=0.2,a=0.5,a=0.7,a=1以及线性信号的瞬时频率随时间变化的关系如图3所示。当a=0和a=0.2时补偿高频,当a=1时补偿低频,当 a=0.5补偿25Hz-50Hz频率部分,当a=0.7补偿15Hz-35Hz频率部分。
本发明一种动态参数的可控震源非线性扫描信号设计方法,首先采用线性升频扫描信号对可控震源黏声正演模拟进行试算。起始频率为5Hz,终止频率为80Hz,扫描长度为12s。扫描信号及频谱如图4所示。采用如图5所示的某工区的模型进行试算,其中图5a为速度模型,图5b为品质因子Q模型,模型大小为3835m×3335m,网格间距为5m。垂直方向的定向力源在地表(1000m,0m)位置激发,记录时间为2s,时间采样间隔为0.5ms,767 个检波器均匀地分布于地表处。为了进行对比,分别采用声波介质正演模拟与黏声介质正演模拟进行试算,得到相关后的地震记录如图6所示,其中,图6a为声波介质正演模拟得到的炮记录,图6b为黏声介质正演模拟得到的炮记录。从炮记录可以看出,经过互相关之后,可控震源正演模拟得到了与炸药震源正演模拟相似的炮记录。对比图6a和图6b,采用黏声介质正演模拟得到的炮记录的能量存在明显的衰减,为了更清楚地进行对比,抽取了第151道的波形曲线如图7a所示,图中黑色实线表示声波正演模拟得到的波形,黑色虚线表示黏声正演模拟得到的波形。从单道记录中可以看出,声介质中地震波的能量比黏声介质中强很多。当可控震源产生的地震波在衰减介质中传播时,地震波会随着传播的进行而衰减,而在地震波无衰减介质中传播时只存在几何扩散衰减,0.3s~1.2s的局部放大图也能够证明这一点,而且可以看出,衰减介质的影响也会造成一定的相位差异。图8为从炮记录图(6)中抽取的第151道频谱。
本发明采用二次多项式动态参数扫描信号对图5所示的黏声介质模型进行可控震源正演模拟测试。动态参数a分别取a=1,a=2,a=-1,a=-2。四种参数的扫描信号如图9所示。图 10为对应的二次多项式动态参数扫描信号的频谱。通过频谱图可以看到,图10a和10c都是对地震波传播中衰减比较厉害的高频成分进行补偿,不同的是当a=1时,高频补偿明显(图 10a),而当a=2时,高频补偿没那么显著了。随着a的增大,高频补偿越来越小,当a趋近于正无穷时,该信号变为线性升频扫描信号,但实际上当a>10时,信号几乎就变成线性升频扫描信号了。同样,图10b和10d都是对低频能量进行加强,但是不同的是a=-1(10b)的加强幅度比a=-2(10d)大的很多,在10d中,频谱低频方向已经没有10b那么陡,随着a的减小,低频补偿越来越小,当a趋向负无穷小时,信号趋近线性信号。
采用上述四种动态参数扫描信号进行正演模拟,得到二次多项式动态参数扫描信号可控震源正演模拟炮记录如图11所示。从图11炮记录中抽取的第120道单道记录如图12所示。通过对不同振幅补偿方式的信号正演结果分析,可以看出,当能量从高频向低频转移时,主瓣分辨率有所降低,但是换来的是深层反射能量的增强。当能量从低频向高频转移时,主瓣分辨率提高,但是深层反射能量较弱,而且旁瓣能量较强,通过增加终止扫描处的斜坡长度可以减少旁瓣能量。对于某些特殊的勘探目的,比如浅层勘探,我们可以适当的补偿高频能量来提高分辨率;相反,当勘探目的层较深,补偿高平能量已经不能抵消越来越快的衰减作用,我们可以增加此时的优势频带,即衰减较少的低频成分的能量来实现我们的勘探目的。在选择参数时,根据地质目标同时兼顾深层反射能量、信噪比和分辨率,对于本工区模型,采用a=-2时得到的地震记录(图11d)在分辨率上与线性扫描信号得到的炮记录(图6b)相当,但深层能量明显更高,因此,选择a=-2时的二次多项式动态扫描信号更佳。
本发明对提出的指数动态参数扫描信号进行正演模拟与分析。同样采用图5所示的黏声模型进行测试,采用相同的扫描参数、正演模拟参数及观测系统。其中动态参数分别选用a=0; a=0.2;a=0.5;a=0.7;a=1得到的可控震源相关后的地震记录如图13所示。抽取的120道的单道记录如图14所示。从图中可以看出,当选用a=0;a=0.2的高频补偿参数时,得到的地震记录(图13a,b和图14a,b)主瓣分辨率更高,旁瓣干扰也更明显,但深层反射能量稍弱;而采用a=1的低频补偿参数时,信噪比更高,深部能量更强。当采用a=0.5补偿25Hz-50Hz 频率部分的参数和采用a=0.7补偿15Hz-35Hz频率部分的参数时,得到的炮记录(图13c,d) 更加接近于炸药震源得到的炮记录,波形曲线更加接近于雷克子波信号。当进行可控震源地震勘探时,可通过调整参数根据地质目标选择优势频带进行补偿。
本发明在正演模拟中加入谐波干扰,这里只加入二阶谐波和三阶谐波,采用两台震源车滑动扫描的采集方式,两台震源车之间的滑动时间为2s,通过滑动扫描谐波干扰特征可知,因为采用线性升频扫描信号,谐波出现在负时间轴上,因此只在第一炮上出现谐波干扰,而第二炮不受谐波干扰的污染,因此这里只展示了受到谐波干扰的第一炮的记录。图15(a1-e1) 所示的为指数动态扫描信号得到的包含谐波干扰的可控震源炮记录,图15(a2-e2)所示的为对应的谐波干扰。从图中,我们可以看出,随着参数a从0增大到1,频率从高频补偿到低频,谐波干扰的影响越来越弱。当参数a等于0,谐波干扰能量非常弱,几乎不对有效能量造成干扰。为了更清楚地展示,我们抽取了谐波干扰影响最大的零偏移距道的单道记录如图 16所示,从图中可以看出,随着补偿频率由高到低,谐波干扰的频率也逐渐变低,当采用参数a等于0时,几乎不存在谐波干扰。也可以在一定程度上证明,在黏弹性地质条件下,采用低频补偿的可控震源非线性扫描信号可以很好地压制谐波干扰。
本发明在常规线非线性扫描信号设计方法的基础上,引入了动态参数,实现了动态参数的可控震源非线性扫描信号设计。首先,基于瞬时频率的两个二次多项式,进行求导,引入动态参数a,实现了二次多项式动态参数的扫描信号设计,通过调节动态参数a,来实现动态高频或者低频补偿。针对可控震源地震勘探中,在某些特定目标区域,需要对中间频带区域进行补偿,设计了指数动态参数扫描信号;其次,分别采用线性升频扫描信号、二次多项式动态参数扫描信号和指数动态参数扫描信号进行可控震源正演模拟分析和测试,验证了本发明的有效性;最后在正演模拟中加入谐波干扰,通过测试不同参数非线性扫描信号对谐波干扰的影响,得出在衰减介质中,采用低频补偿的扫描信号可以较好地压制谐波干扰。因此,本实施方式提出的一种动态参数的可控震源非线性扫描信号设计系统及方法,相比于传统线性扫描信号,通过调整参数,二次多项式动态参数扫描信号可实现高频或者低频补偿,而指数动态参数扫描信号不仅可以补偿高频或者低频,也可以基于地质目标实现中间频带的加强。
当然,上述说明并非是对本发明的限制,本发明也并不仅限于上述举例,本技术领域的技术人员在本发明的实质范围内所做出的变化、改型、添加或替换,也应属于本发明的保护范围。
Claims (4)
1.一种动态参数的可控震源非线性扫描信号设计系统,其特征在于,包括二次多项式动态参数扫描信号设计模块、指数动态参数扫描信号设计模块、可控震源黏声正演模块、二次多项式动态参数扫描信号正演模拟与分析模块、指数动态参数扫描信号正演模拟与分析模块以及动态参数扫描信号对谐波干扰影响模块;
二次多项式动态参数扫描信号设计模块,被配置为用于对二次多项式动态参数扫描信号进行设计;
指数动态参数扫描信号设计模块,被配置为用于对指数动态参数扫描信号进行设计;
可控震源黏声正演模块,被配置为用于对线性扫描信号进行可控震源黏声正演模拟;
二次多项式动态参数扫描信号正演模拟与分析模块,被配置为用于对二次多项式动态参数扫描信号的可控震源进行正演模拟测试和效果分析;
指数动态参数扫描信号正演模拟与分析模块,被配置为用于指数动态参数扫描信号的可控震源进行正演模拟测试和频带补偿效果分析;
动态参数扫描信号对谐波干扰影响模块,被配置为用于在指数动态扫描信号和二次多项式动态扫描信号正演模拟中加入谐波,测试指数动态参数扫描信号和二次多项式动态参数扫描信号对谐波干扰的影响。
2.一种动态参数的可控震源非线性扫描信号设计方法,其特征在于,采用如权利要求1所述的动态参数的可控震源非线性扫描信号设计系统,包括如下步骤:
步骤1:通过二次多项式动态参数扫描信号设计模块,设计二次多项式动态参数扫描信号,其表达式如公式(5)所示;
其中,f1和f2分别表示起始频率和终止频率,T为扫描周期,a为动态参数;
步骤2:通过指数动态参数扫描信号设计模块,设计指数动态参数扫描信号,其表达式如公式(6)所示;
步骤3:通过可控震源黏声正演模块,对线性扫描信号进行可控震源黏声正演模拟,二维黏滞声波方程为:
其中,P表示为压力场,τσl和τεl分别表示为弛豫时间常数,GR表示弛豫模量;L代表标准线性体的个数,rl表示记忆变量,vx和vz分别为水平分量和垂直分量的速度;
步骤4:通过二次多项式动态参数扫描信号正演模拟与分析模块,测试二次多项式动态扫描信号对地震波高频或者低频的补偿效果,测试参数为:起始频率为5Hz,终止频率为80Hz,扫描长度为12s,动态参数a分别取a=1,a=2,a=-1,a=-2;
步骤5:通过指数动态参数扫描信号正演模拟与分析模块,测试指数动态扫描信号对地震波不同频率成分的补偿效果,测试参数为:起始频率为5Hz,终止频率为80Hz,扫描长度为12s,动态参数a分别取a=0,a=0.2,a=0.5,a=0.7,a=1;
步骤6:通过动态参数扫描信号对谐波干扰影响模块,在正演模拟中加入二阶谐波和三阶谐波,测试动态参数扫描信号对谐波干扰的影响。
3.根据权利要求2所述的动态参数的可控震源非线性扫描信号设计方法,其特征在于,在步骤1中,二次多项式动态参数扫描信号公式(5)的推导过程如下:
首先,设计瞬时频率的两个二次多项式:
分别对上式进行求导:
其中,F′1_1(t)的斜率为F′1_2(t)的斜率为由此可知,瞬时频率F1_1(t)的增加速度越来越快,表示低频补偿;瞬时频率F1_2(t)的增加速度越来越慢,表示高频补偿;根据公式(1)和(2),引入动态参数a,实现二次多项式动态参数扫描信号设计,通过调节a,可实现对高频或者低频进行动态补偿;
二次多项式动态参数扫描信号F1(t)满足如下几个条件:
(1)F1(t)max=F1(T)=f2;
(2)F1(t)min=F1(0)=f1;
(3)当a>0时,F′1_1(t)单调递增;
(4)当a<0时,F′1_2(t)单调递减;
根据以上四个条件并结合公式(1)和(2),可推导出二次多项式动态参数扫描信号公式(5)。
4.根据权利要求2所述的动态参数的可控震源非线性扫描信号设计方法,其特征在于,在步骤3中,二维黏滞声波方程中,为了方便求解,将二阶导数进行降阶处理,转化为一阶函数求导,且引入的辅助微分算子如下所示:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810114737.4A CN108318919B (zh) | 2018-02-06 | 2018-02-06 | 一种动态参数的可控震源非线性扫描信号设计系统及方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810114737.4A CN108318919B (zh) | 2018-02-06 | 2018-02-06 | 一种动态参数的可控震源非线性扫描信号设计系统及方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108318919A true CN108318919A (zh) | 2018-07-24 |
CN108318919B CN108318919B (zh) | 2020-05-15 |
Family
ID=62902004
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810114737.4A Active CN108318919B (zh) | 2018-02-06 | 2018-02-06 | 一种动态参数的可控震源非线性扫描信号设计系统及方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108318919B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111665536A (zh) * | 2019-03-05 | 2020-09-15 | 中石化石油工程技术服务有限公司 | 基于微测井子波定量化约束的井深设计方法 |
CN112505747A (zh) * | 2020-12-22 | 2021-03-16 | 吉林大学 | 基于多信号发生器协同可控震源振动畸变抑制系统及方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103777240A (zh) * | 2012-10-25 | 2014-05-07 | 中国石油化工股份有限公司 | 一种可控震源分频扫描频带补偿方法 |
CN104035128A (zh) * | 2013-03-06 | 2014-09-10 | 中国石油集团东方地球物理勘探有限责任公司 | 可控震源伪随机扫描信号生成方法 |
CN104950326A (zh) * | 2015-03-04 | 2015-09-30 | 中石化石油工程地球物理有限公司胜利分公司 | 基于目的层频谱的可控震源非线性扫描信号的设计方法 |
CN105445783A (zh) * | 2015-10-08 | 2016-03-30 | 吉林大学 | 一种适于复杂地表条件的电磁式可控震源广义预测控制方法 |
-
2018
- 2018-02-06 CN CN201810114737.4A patent/CN108318919B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103777240A (zh) * | 2012-10-25 | 2014-05-07 | 中国石油化工股份有限公司 | 一种可控震源分频扫描频带补偿方法 |
CN104035128A (zh) * | 2013-03-06 | 2014-09-10 | 中国石油集团东方地球物理勘探有限责任公司 | 可控震源伪随机扫描信号生成方法 |
CN104950326A (zh) * | 2015-03-04 | 2015-09-30 | 中石化石油工程地球物理有限公司胜利分公司 | 基于目的层频谱的可控震源非线性扫描信号的设计方法 |
CN105445783A (zh) * | 2015-10-08 | 2016-03-30 | 吉林大学 | 一种适于复杂地表条件的电磁式可控震源广义预测控制方法 |
Non-Patent Citations (2)
Title |
---|
李振春等: "可控震源两种谐波产生机理与特征研究", 《石油物探》 * |
王忠仁等: "可控震源非线性扫描地震响应的数值模拟", 《地球物理学进展》 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111665536A (zh) * | 2019-03-05 | 2020-09-15 | 中石化石油工程技术服务有限公司 | 基于微测井子波定量化约束的井深设计方法 |
CN111665536B (zh) * | 2019-03-05 | 2024-01-09 | 中石化石油工程技术服务有限公司 | 基于微测井子波定量化约束的井深设计方法 |
CN112505747A (zh) * | 2020-12-22 | 2021-03-16 | 吉林大学 | 基于多信号发生器协同可控震源振动畸变抑制系统及方法 |
CN112505747B (zh) * | 2020-12-22 | 2021-10-01 | 吉林大学 | 基于多信号发生器协同可控震源振动畸变抑制系统及方法 |
Also Published As
Publication number | Publication date |
---|---|
CN108318919B (zh) | 2020-05-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103376464B (zh) | 一种地层品质因子反演方法 | |
CN103885085B (zh) | 一种压制可控震源谐波干扰的方法 | |
CN102590869B (zh) | 一种人工场源频率域电法勘探方法及勘探系统 | |
CN102062872A (zh) | 可控震源多因子扫描方法 | |
CN103698807B (zh) | 梯状二维宽线观测系统设计方法 | |
CN109143374B (zh) | 一种井周散射体成像方法及系统 | |
CN104950326B (zh) | 基于目的层频谱的可控震源非线性扫描信号的设计方法 | |
CN109782335A (zh) | 低畸变宽频扫描信号设计方法 | |
CN108318919A (zh) | 一种动态参数的可控震源非线性扫描信号设计系统及方法 | |
CN106680874A (zh) | 基于波形形态特征稀疏化建模的谐波噪声压制方法 | |
CN104483704A (zh) | 基于avo异常类型约束的剩余相位校正方法 | |
CN109765626A (zh) | 一种基于最小二乘克希霍夫偏移的探月雷达数据处理方法 | |
CN103728661B (zh) | 一种高精度反q滤波地震资料处理方法 | |
CN102269822B (zh) | 一种混合的地层吸收补偿方法 | |
CN105277986A (zh) | 基于自适应匹配滤波算子的可控震源谐波压制方法 | |
CN102338888B (zh) | 提高地震分辨率的可控震源资料相关方法 | |
CN104570103B (zh) | 一种低信噪比地震资料的井约束速度谱拾取方法 | |
CN104793237A (zh) | 一种获得宽频可控震源扫描信号的方法和装置 | |
CN100412569C (zh) | 利用地震微测井进行地震信号高频补偿方法 | |
CN107238864A (zh) | 基于优势频带约束的地震波能量时频补偿方法 | |
CN109270574A (zh) | 一种基于多种震源采集数据的联合反褶积方法 | |
CN105277987A (zh) | 基于预测滤波法和纯相移法的可控震源谐波压制方法 | |
CN103777240B (zh) | 一种可控震源分频扫描频带补偿方法 | |
CN109975874B (zh) | 基于阻尼雷克子波的可控震源扫描信号设计方法 | |
CN107918152B (zh) | 一种地震相干层析成像方法 |
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 | ||
CB02 | Change of applicant information |
Address after: 065000 No. 84 Jinguang Road, Guangyang District, Langfang City, Hebei Province Applicant after: Institute of Geophysical and Geochemical Exploration under China Academy of Geos Applicant after: China Petroleum University (East China) Address before: 065000 No. 84 Jinguangdao, Guangyang District, Handan City, Hebei Province Applicant before: Institute of Geophysical and Geochemical Exploration under China Academy of Geos Applicant before: China Petroleum University (East China) |
|
CB02 | Change of applicant information | ||
GR01 | Patent grant | ||
GR01 | Patent grant |