CN108761534B - 陆上地震加速度信号应用新方法 - Google Patents
陆上地震加速度信号应用新方法 Download PDFInfo
- Publication number
- CN108761534B CN108761534B CN201810486300.3A CN201810486300A CN108761534B CN 108761534 B CN108761534 B CN 108761534B CN 201810486300 A CN201810486300 A CN 201810486300A CN 108761534 B CN108761534 B CN 108761534B
- Authority
- CN
- China
- Prior art keywords
- acceleration
- signal
- seismic
- noise
- velocity
- 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
- 230000001133 acceleration Effects 0.000 title claims abstract description 150
- 238000000034 method Methods 0.000 title claims abstract description 40
- 238000012545 processing Methods 0.000 claims abstract description 18
- 238000003384 imaging method Methods 0.000 claims abstract description 17
- 238000005516 engineering process Methods 0.000 claims abstract description 10
- 230000001629 suppression Effects 0.000 claims abstract description 5
- 230000007613 environmental effect Effects 0.000 claims description 7
- 238000006243 chemical reaction Methods 0.000 claims description 6
- 230000009466 transformation Effects 0.000 claims description 6
- 238000004088 simulation Methods 0.000 claims description 5
- 230000015572 biosynthetic process Effects 0.000 claims description 3
- 238000006073 displacement reaction Methods 0.000 claims description 3
- 230000008569 process Effects 0.000 claims description 3
- 229910052704 radon Inorganic materials 0.000 claims description 2
- SYUHGPGVQRZVTB-UHFFFAOYSA-N radon atom Chemical compound [Rn] SYUHGPGVQRZVTB-UHFFFAOYSA-N 0.000 claims description 2
- 230000008859 change Effects 0.000 abstract description 16
- 230000004044 response Effects 0.000 abstract description 12
- 230000035945 sensitivity Effects 0.000 abstract description 5
- 230000000704 physical effect Effects 0.000 abstract description 4
- 238000012360 testing method Methods 0.000 description 12
- 230000000694 effects Effects 0.000 description 10
- 238000011161 development Methods 0.000 description 5
- 230000006872 improvement Effects 0.000 description 5
- 238000004458 analytical method Methods 0.000 description 4
- 238000001914 filtration Methods 0.000 description 4
- 239000011435 rock Substances 0.000 description 4
- 238000010586 diagram Methods 0.000 description 3
- 230000001965 increasing effect Effects 0.000 description 3
- 239000003921 oil Substances 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 238000001228 spectrum Methods 0.000 description 3
- 230000008901 benefit Effects 0.000 description 2
- 230000008878 coupling Effects 0.000 description 2
- 238000010168 coupling process Methods 0.000 description 2
- 238000005859 coupling reaction Methods 0.000 description 2
- 238000005755 formation reaction Methods 0.000 description 2
- 230000006870 function Effects 0.000 description 2
- 238000009533 lab test Methods 0.000 description 2
- 230000007246 mechanism Effects 0.000 description 2
- 239000013307 optical fiber Substances 0.000 description 2
- 235000019738 Limestone Nutrition 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 239000010620 bay oil Substances 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000005553 drilling Methods 0.000 description 1
- 230000002708 enhancing effect Effects 0.000 description 1
- 230000005284 excitation Effects 0.000 description 1
- 239000006028 limestone Substances 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
- 239000004576 sand Substances 0.000 description 1
- 230000003595 spectral effect Effects 0.000 description 1
- 230000008961 swelling Effects 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/30—Analysis
- G01V1/306—Analysis for determining physical properties of the subsurface, e.g. impedance, porosity or attenuation profiles
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/30—Noise handling
- G01V2210/32—Noise reduction
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/40—Transforming data representation
- G01V2210/46—Radon transform
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/40—Transforming data representation
- G01V2210/48—Other transforms
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/62—Physical property of subsurface
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
本发明提供一种陆上地震加速度信号应用新方法,该陆上地震加速度信号应用新方法包括:步骤1,获取加速度正演模拟地震记录,估算加速度信号采集覆盖次数,采集地震加速度信号;步骤2,压制与消除高频噪声;步骤3,进行加速度地震信号成像,将经过高频噪声压制处理后的数据应用地震资料处理技术来完成加速度地震资料成像处理。该陆上地震加速度信号应用新方法有效提高了地震勘探的成像精度,主要表现在地震信号的保真度高、分辨率高、频带宽、对介质物性变化响应灵敏度高,是陆上提高地震勘探精度的有效方法。
Description
技术领域
本发明涉及勘探地球物理信号应用技术领域,特别是涉及到一种陆上地震加速度信号应用新方法。
背景技术
油气勘探中地质勘探目标越来越复杂,对地震勘探精度的要求也越来越高,东部渤海湾油田对3-5km深度地层分辨要求由10-15m减小到5-10m,勘探精度要求提高并不意味着勘探技术已经得到了相应提高,事实上,目前陆上地震勘探所用的检波器主要是常规动圈检波器,接收的地震波是速度型信号,这种地震信号频率最高只能到100-120Hz,相对较窄的频带限制了分辨地质目标的能力,要分辨更薄的地层存在较大困难。为了寻找提高地震分辨率的方法1993年李庆忠院士在专著中提出了高分辨率地震勘探系统工程的概念,论述了开展高分辨率地震勘探的相关技术,俞寿朋对地震波的激发、接收及资料处理等方面提出了改进方案,高分辨率地震勘探接收技术主要提高了检波器的性能和检波器的耦合能力,促使了检波器类型的增多和耦合工艺的提高。1993年荷兰SENSOR公司首先推出了SM-4SH允差更小的超级检波器,接着美国WESTERN GEOPHASICAL和MARK公司、日本OYO公司、国内西安石油仪器厂也相继推出了超级检波器,但这些改进与提高都没有摆脱对速度信号的依赖。可喜的是1999年以来MEMS、陆用压电、光纤等新型加速度检波器的研发推动了地震波接收技术的发展(Maxwell P.M.et al.,1999;Mougenot,D.,2004;Lv,G.H.,2005;Wang,J.Y.et al.,2016),2007年Hons,M等人进行了MEMS检波器与常规动圈检波器对比试验(Hons,Met al.,2007),2011-2016年Stotter,C等人经过试验取得了良好效果(Stotter,C.Et al.,2011;Bai,X;Liu,Z.D.Et al.,2012.M.Et al.,2014,Criss,J,2016),信号频带得到拓宽,信号保真度、相位一致性提高,并得出加速度检波器接收地震波频带宽、分辨率高的结论。然而由于对加速度信号认识的不足加速度检波器的应用效果并没有得到广泛认可,2013年Wei,J D等人对MEMS、陆用压电加速度检波器提出了质疑,认为MEMS、陆用压电、光纤检波器是加速度型检波器,常规动圈检波器是速度型检波器,两种检波器要转到同一域进行比较,Hons,M.S等也有类似的观点(Hons,M.S.et al.,2008;),并得出MEMS与常规动圈检波器接收效果相似,加速度检波器不能拓宽频带的结论。
分辨率和成像精度是陆上地震勘探发展要解决的关键问题,特别是垂向分辨率一直无法达到复杂地质目标对地震勘探精度的要求,这与采用的地震速度信号有较大关系,因为速度信号具有带通滤波效应,对低频和高频都有压制作用,且受子波拉伸畸变、频带窄、分辨率低的制约,难以对复杂介质进行5m以下高精度分辨,使地震勘探精度无法再上新台阶。速度信号波形畸变大、保真度低,难以完成高保真地震勘探的任务,且速度信号对介质物性变化反应灵敏度低,这都不利于对复杂介质的高精度成像。
由于以上原因,地震速度信号难以满足高精度地震勘探的需要,需要寻求更好的高精度地震勘探方法。为此我们从弹性波理论、岩石物理学、正演模拟、实验室测试、野外试验等方面开展了大量的研究工作,查阅了大量工程地震、矿山监测、大洋探测等领域对探测信号的应用文献,经过试验、分析和研究总结形成了一套地震加速度信号应用新方法,通过多个项目的应用获得了良好的效果。
发明内容
本发明的目的是提供一种通过加速度信号的应用进一步提高地震勘探成像精度,为复杂介质提供更精确的勘探地球物理技术的陆上地震加速度信号应用新方法。
本发明的目的可通过如下技术措施来实现:陆上地震加速度信号应用新方法,该陆上地震加速度信号应用新方法包括:步骤1,获取加速度正演模拟地震记录,估算加速度信号采集覆盖次数,采集地震加速度信号;步骤2,压制与消除高频噪声;步骤3,进行加速度地震信号成像,将经过高频噪声压制处理后的数据应用地震资料处理技术来完成加速度地震资料成像处理。
本发明的目的还可通过如下技术措施来实现:
在步骤1中,在获取加速度正演模拟地震记录时,首先,根据地质模型将地层密度ρ、速度v参数代入波动方程进行正演模拟得到目的层纵波速度信号反射炮集记录,模拟采用完全弹性介质模型,纵波波动方程为:
式中:为位移场标量位,vp为纵波速度;
其次,将速度地震模拟记录转化为加速度地震记录,速度信号v与加速度信号a之间存在关系式:
a=v' (2)。
在步骤1中,在估算加速度信号采集覆盖次数时,根据加速度正演模拟结果确定覆盖次数,提取速度信号地震记录环境噪声,将噪声振幅Av按(3)式转化为加速度信号噪声振幅Aa:
Aa=2πfAv (3)
式中,f为噪声频率。
在步骤1中,在估算加速度信号采集覆盖次数时,将转化得到的噪声数据加入到模拟加速度记录中,对带有噪声的速度和加速度模拟资料分别做不同覆盖次数叠加处理,将加速度信号叠加剖面信噪比与速度信号剖面信噪比相当的覆盖次数确定为采集加速度信号的覆盖次数。
在步骤1中,在估算加速度信号采集覆盖次数时,根据高信噪比速度信号地震资料覆盖次数估算采集加速度信号的覆盖次数,选择信噪比较高的地震剖面反射层位,计算其有效覆盖次数,将速度信号覆盖次数的2-3倍确定为采集加速度信号的覆盖次数。
在步骤1中,在采集地震加速度信号时,用加速度检波器直接加速度信号,将加速度检波器的谐振频率点做到800Hz以上,将地震波频带内的有效信号无畸变地记录下来。
在步骤1中,在采集地震加速度信号时,先用速度检波器采集速度信号,然后根据速度信号v与加速度信号a之间存在的关系式,在室内将速度信号转换为加速度信号。
步骤2包括:
步骤a,在频率域对高频噪声进行压制;
步骤b,将频域加速度信号转换为时域信号;
步骤c,进行噪声去除处理。
在步骤a中,对长度为M的离散信号x(n)按式(4)进行N点离散傅里叶变换:
式中,
然后,将X(k)振幅除以加速度信号噪声降为速度信号水平,提高资料信噪比。
在步骤b中,将步骤a压制高频噪声后得到的频域信号按式(5)进行离散傅里叶反变换:
此时,得到压制高频噪声后的加速度信号。
在步骤c中,用拉东变换、傅立叶变换、小波变换去除面波、线性干扰这些规则干扰,用曲波、压缩感知、字典学习方法去除环境干扰,进一步提高加速度信号信噪比。
本发明中的陆上地震加速度信号应用新方法,根据加速度信号产生机理及其与传播介质之间的关系,通过理论分析、正演模拟、实验室测试及野外试验总结形成了一套陆上地震加速度信号应用方法,应用加速度信号的低通效应、高保真度、宽频带、高分辨率、对介质物性反应灵敏等独特的信号属性提高地震勘探精度。
附图说明
图1为本发明的步骤101中速度、加速度信号正演模拟地震记录图;
图2为本发明的步骤101中速度、加速度信号正演模拟地震剖面图;
图3为本发明的一具体实施例中速度、加速度信号保真度测试对比图;
图4为本发明的一具体实施例中信号拉伸畸变率与频率关系曲线图;
图5为本发明的一具体实施例中速度、加速度信号实验室脉冲测试信号频谱图;
图6为本发明的一具体实施例中速度、加速度信号与介质弹性参数关系曲线图;
图7为本发明的一具体实施例中东风港三维单炮对比图;
图8为本发明的一具体实施例中东风港三维单炮频谱对比图;
图9为本发明的一具体实施例中东风港三维剖面对比图;
图10为本发明的一具体实施例中罗家三维单炮对比图;
图11为本发明的一具体实施例中罗家三维单炮频谱对比图;
图12为本发明的一具体实施例中罗家三维剖面对比图;
图13为本发明的一具体实施例中罗家三维时间切片效果对比图;
图14为本发明的陆上地震加速度信号应用新方法的一具体实施例的流程图。
具体实施方式
为使本发明的上述和其他目的、特征和优点能更明显易懂,下文特举出较佳实施例,并配合附图所示,作详细说明如下。
如图14所示,图14为本发明的陆上地震加速度信号应用新方法的流程图。
步骤101、地震加速度信号采集
1)获取加速度正演模拟地震记录
首先,根据地质模型将地层密度ρ、速度v参数代入波动方程进行正演模拟得到目的层纵波速度信号反射炮集记录(见图1)。模拟采用完全弹性介质模型,纵波波动方程为:
式中:为位移场标量位,vp为纵波速度。
其次,将速度地震模拟记录转化为加速度地震记录,速度信号v与加速度信号a之间存在关系式:
a=v' (2)
2)加速度信号采集覆盖次数估算
加速度信号具有增强高频的作用,频率越高增强倍数越大,对地震波而言,高频部分多为环境噪声,这使得加速度信号频噪声较强,信噪比降低,在采集加速度信号时必须采用更高的覆盖次数才能保证信号有足够高的信噪比。
覆盖次数通过两种方法综合确定:
i).根据高信噪比速度信号地震资料覆盖次数估算采集加速度信号的覆盖次数,选择信噪比较高的地震剖面反射层位,计算其有效覆盖次数,将速度信号覆盖次数的2-3倍确定为采集加速度信号的覆盖次数;
ii).根据加速度正演模拟结果确定覆盖次数,提取速度信号地震记录环境噪声,将噪声振幅Av按(3)式转化为加速度信号噪声振幅Aa:
Aa=2πfAv (3)
式中,f为噪声频率。
然后将转化得到的噪声数据加入到模拟加速度记录中,对带有噪声的速度和加速度模拟资料分别做不同覆盖次数叠加处理(见图2),将加速度信号叠加剖面信噪比与速度信号剖面信噪比相当的覆盖次数确定为采集加速度信号的覆盖次数。
最后,综合考虑①和②两步分析最终确定加速度信号采集覆盖次数。
3)加速度信号采集
可用两种方法采集:
i).用加速度检波器直接加速度信号,这种采集方法最好,因为加速度检波器具有低通滤波效应,将检波器的谐振频率点做到800Hz以上,就可将地震波频带内的有效信号无畸变地记录下来;而速度检波器因为具有带通滤波作用,所以地震波低频和高频都受到一定程度压制,从而使地震波波形发生畸变。
ii).先用速度检波器采集速度信号,然后根据(2)式在室内将速度信号转换为加速度信号,但这种方式采集到的加速度信号存在两种损失:第一,速度检波器的带通滤波作用使地震波频带变窄,波形发生一定程度畸变;第二,信号转化过程中数值计算会引入一部分噪声。这种采集方式适合于采集精度要求相对较低的情况。
步骤102、高频噪声压制与消除
加速度信号与速度信号相比高频噪声得到增强,可分三步对噪声进行消除处理:
第一步,在频率域对高频噪声进行压制
首先,对长度为M的离散信号x(n)按式(4)进行N点离散傅里叶变换:
式中,
然后,将X(k)振幅除以加速度信号噪声降为速度信号水平,提高资料信噪比。
第二步,频域加速度信号转换为时域信号
将第一步压制高频噪声后得到的频域信号按式(5)进行离散傅里叶反变换:
此时,得到压制高频噪声后的加速度信号。
第三步,噪声去除处理
用拉东Radon变换、傅里叶Fourier变换、小波变换去除面波、线性干扰等规则干扰,用曲波、压缩感知、字典学习等方法去除环境干扰,进一步提高加速度信号信噪比。
步骤103、加速度地震信号成像
经过高频噪声压制处理后的数据即可应用常规地震资料处理技术来完成加速度地震资料成像处理。
在应用本发明的具体实施例中,通过在东风港和罗家两个三维开发地震项目中应用地震加速度信号,地震资料品质得到明显提高,获得了良好的地质勘探效果,优质完成了油气资源勘探任务。
1、东风港开发三维地震的应用
在东风港地震采集项目中做了加速度信号保真度、频带宽度、介质物性响应灵敏度实验室测试与分析。图3是脉冲信号测试,看到DSU3和BK-Z两种加速度检波器信号波形与测试波形非常相似,信号尾部旁瓣小,拉伸畸变小(畸变率低于10%);PS-5R、20DX速度检波器波形则发生了明显畸变,信号脉宽由5ms拉伸到14ms(拉伸率180%),对比不同频率脉冲测试波形,畸变率与频率近于线性增大关系(图2),且在略高于检波器自然频率(20DX检波器10Hz)处存在一个拐点,低于自然频率的信号拉伸率更大,信号畸变更严重。速度检波器信号尾部旁瓣幅值高,保真度和分辨率较低。
图4可以看到脉冲测试信号的脉宽由5ms逐渐增大到20ms,即频率由高到低测试,分析频率变化与拉伸畸变的关系,看到检波器拉伸率与检波器自然频率成反比,频率越高拉伸率越低,低频检波器信号畸变最大。
观察图5脉冲测试记录频谱,DSU3和BK-Z加速度检波器58Hz以上高频信号强,且频带宽;PS-5Hz、20DX速度型检波器低频好,但频带窄;PS-28Hz检波器低频受到压制与加速度检波器相似,高频与速度检波器相似。
图6是速度、加速度信号对介质物性响应图,看出加速度响应密度变化的灵敏度明显高于速度的响应,密度越小加速度值变化越大;而速度响应密度的变化比较平缓。保持密度值不变,分析纵波速度、加速度随拉梅常数的变化规律,看出纵波加速度对拉梅常数的变化响应比速度大的多,变化近似呈线性关系,速度响应的变化量相对较小。对于剪切模量而言,加速度对剪切模量较小的疏松岩层响应变化较大,对剪切模量大于1GPa的火成岩、灰岩等致密岩石响应变化变缓,这说明对剪切模量较大的岩层纵波加速度响应涨缩变化的能力比响应剪切变化强。分析横波速度、加速度随密度的变化,当密度值较小时加速度变化较明显,随着密度的增大变化趋于平缓;速度对密度变化的反应则小得多。对剪切模量的响应,加速度的响应幅度呈线性变化,相比速度变化有明显的优势。
东风港开发地震主要任务是提高地震资料分辨率,老的地震资料上在车西洼陷南斜坡的浊积岩信噪比、分辨率低,波组连续性差,难于准确追踪描述砂体,小断块的断点不清、低序级断层、低幅构造刻画难度大,构造落实程度差,使地质人员确定钻井增加了风险。本次采集选择了LHKJ-1A陆用压电加速度型单个检波器接收,如图7-9所示获得的地震资料提高了反射波的保真度,拓宽了单炮频带25Hz,偏移剖面目的层断点、小断块得到精确成像,同相轴连续性改善,地层分辨能力显著提高,圆满完成了预期目标。
2、罗家三维开发地震的应用
如图10-13所示,罗家三维开发地震老资料Es3、Es4、Mz地层在信噪比和成像方面都比较差,难以满足精细构造研究及储层预测的要求,本次采集要提高地震资料分辨率和信噪比,选择LHKJ-1A陆用压电加速度型单个检波器接收,通过加速度信号的应用单炮频带拓宽40-50Hz,地震剖面成像精度得到显著提高,在水平时间切片上,加速度陆用压电检波器资料不但河道、断层清晰可见,而且小尺度纵横交错的断裂也能有效分辨,相比老资料有了质的提升。
从时间切片看,加速度信号资料信噪比、分辨率显著提高,地质构造清晰,浅层小断裂、河道变化趋势明显。
本发明中的陆上地震加速度信号应用新方法是在加速度信号产生机理基础上分析了信号与地下传播介质之间的关系和变化规律,根据加速度信号保真度高、频带宽、子波分辨率高,但高频噪声强的特点设计了加速度信号应用方法,这种方法提高了地震信号的保真度,拓宽了地震波频带,压缩了地震子波,提高了地震资料的分辨率和成像精度,在多个工区的应用中取得了良好的效果,说明了加速度信号对提高复杂构造区勘探精度非常有效,在勘探难度越来越大,复杂程度越来越高的油气勘探领域具有广阔的应用前景。
Claims (8)
1.陆上地震加速度信号应用新方法,其特征在于,该陆上地震加速度信号应用新方法包括:
步骤1,获取加速度正演模拟地震记录,估算加速度信号采集覆盖次数,采集地震加速度信号;
步骤2,压制与消除高频噪声;
步骤3,进行加速度地震信号成像,将经过高频噪声压制处理后的数据应用地震资料处理技术来完成加速度地震资料成像处理;
步骤2包括:步骤a,在频率域对高频噪声进行压制;步骤b,将频域加速度信号转换为时域信号;步骤c,进行噪声去除处理;
在步骤a中,对长度为M的离散信号x(n)按式(4)进行N点离散傅里叶变换:
式中,
然后,将X(k)振幅除以加速度信号噪声降为速度信号水平,提高资料信噪比;
在步骤b中,将步骤a压制高频噪声后得到的频域信号按式(5)进行离散傅里叶反变换:
此时,得到压制高频噪声后的加速度信号。
2.根据权利要求1所述的陆上地震加速度信号应用新方法,其特征在于,在步骤1中,在获取加速度正演模拟地震记录时,首先,根据地质模型将地层密度ρ、速度vp参数代入波动方程进行正演模拟得到目的层纵波速度信号反射炮集记录,模拟采用完全弹性介质模型,纵波波动方程为:
式中:为位移场标量位,vp为纵波速度;
其次,将速度地震模拟记录转化为加速度地震记录,速度信号v与加速度信号a之间存在关系式:
a=v' (2)。
3.根据权利要求1所述的陆上地震加速度信号应用新方法,其特征在于,在步骤1中,在估算加速度信号采集覆盖次数时,根据加速度正演模拟结果确定覆盖次数,提取速度信号地震记录环境噪声,将噪声振幅Av按(3)式转化为加速度信号噪声振幅Aa:
Aa=2πfAv (3)
式中,f为噪声频率。
4.根据权利要求3所述的陆上地震加速度信号应用新方法,其特征在于,在步骤1中,在估算加速度信号采集覆盖次数时,将转化得到的噪声数据加入到模拟加速度记录中,对带有噪声的速度和加速度模拟资料分别做不同覆盖次数叠加处理,将加速度信号叠加剖面信噪比与速度信号剖面信噪比相当的覆盖次数确定为采集加速度信号的覆盖次数。
5.根据权利要求1所述的陆上地震加速度信号应用新方法,其特征在于,在步骤1中,在估算加速度信号采集覆盖次数时,根据高信噪比速度信号地震资料覆盖次数估算采集加速度信号的覆盖次数,选择信噪比较高的地震剖面反射层位,计算其有效覆盖次数,将速度信号覆盖次数的2-3倍确定为采集加速度信号的覆盖次数。
6.根据权利要求1所述的陆上地震加速度信号应用新方法,其特征在于,在步骤1中,在采集地震加速度信号时,用加速度检波器直接加速度信号,将加速度检波器的谐振频率点做到800Hz以上,将地震波频带内的有效信号无畸变地记录下来。
7.根据权利要求1所述的陆上地震加速度信号应用新方法,其特征在于,在步骤1中,在采集地震加速度信号时,先用速度检波器采集速度信号,然后根据速度信号v与加速度信号a之间存在的关系式,在室内将速度信号转换为加速度信号。
8.根据权利要求1所述的陆上地震加速度信号应用新方法,其特征在于,在步骤c中,用拉东变换、傅立叶变换、小波变换去除面波、线性干扰,用曲波、压缩感知、字典学习方法去除环境干扰,进一步提高加速度信号信噪比。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810486300.3A CN108761534B (zh) | 2018-05-18 | 2018-05-18 | 陆上地震加速度信号应用新方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810486300.3A CN108761534B (zh) | 2018-05-18 | 2018-05-18 | 陆上地震加速度信号应用新方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108761534A CN108761534A (zh) | 2018-11-06 |
CN108761534B true CN108761534B (zh) | 2024-03-29 |
Family
ID=64007189
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810486300.3A Active CN108761534B (zh) | 2018-05-18 | 2018-05-18 | 陆上地震加速度信号应用新方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108761534B (zh) |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102520441A (zh) * | 2011-11-30 | 2012-06-27 | 朱德兵 | 一种井巷定向探测地震勘探数据采集方法 |
CN103412336A (zh) * | 2013-07-22 | 2013-11-27 | 中国石油天然气股份有限公司 | 一种非均质油藏中岩石系统的纵波速度预测方法 |
WO2015196779A1 (zh) * | 2014-06-24 | 2015-12-30 | 夏正元 | 利用叠加速度求取高精度地震波速度的方法 |
CN107065013A (zh) * | 2016-12-27 | 2017-08-18 | 中国石油大学(华东) | 一种地震尺度下的层速度确定方法及装置 |
CN107132578A (zh) * | 2017-04-06 | 2017-09-05 | 吉林大学 | 一种微地震地面监测速度模型校正算法 |
CN107479092A (zh) * | 2017-08-17 | 2017-12-15 | 电子科技大学 | 一种基于方向导数的频率域高阶声波方程正演模拟方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8700372B2 (en) * | 2011-03-10 | 2014-04-15 | Schlumberger Technology Corporation | Method for 3-D gravity forward modeling and inversion in the wavenumber domain |
-
2018
- 2018-05-18 CN CN201810486300.3A patent/CN108761534B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102520441A (zh) * | 2011-11-30 | 2012-06-27 | 朱德兵 | 一种井巷定向探测地震勘探数据采集方法 |
CN103412336A (zh) * | 2013-07-22 | 2013-11-27 | 中国石油天然气股份有限公司 | 一种非均质油藏中岩石系统的纵波速度预测方法 |
WO2015196779A1 (zh) * | 2014-06-24 | 2015-12-30 | 夏正元 | 利用叠加速度求取高精度地震波速度的方法 |
CN107065013A (zh) * | 2016-12-27 | 2017-08-18 | 中国石油大学(华东) | 一种地震尺度下的层速度确定方法及装置 |
CN107132578A (zh) * | 2017-04-06 | 2017-09-05 | 吉林大学 | 一种微地震地面监测速度模型校正算法 |
CN107479092A (zh) * | 2017-08-17 | 2017-12-15 | 电子科技大学 | 一种基于方向导数的频率域高阶声波方程正演模拟方法 |
Also Published As
Publication number | Publication date |
---|---|
CN108761534A (zh) | 2018-11-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CA2684737A1 (en) | Methods of hydrocarbon detection using spectral energy analysis | |
CN104237945A (zh) | 一种地震资料自适应高分辨处理方法 | |
CN109884709B (zh) | 一种基于面波旅行时层析的转换波静校正方法 | |
CN213813970U (zh) | 基于分布式光纤声波传感的时移vsp数据采集系统 | |
RU2412454C2 (ru) | Способ обработки сейсмических данных с использованием дискретного вейвлет-преобразования | |
CN109541689B (zh) | 一种基于反射波能量特征的介质密实度评价方法 | |
CN102323618B (zh) | 基于分数阶傅里叶变换的相干噪声抑制方法 | |
Zeng et al. | High‐resolution shallow structure at brady hot springs using ambient noise tomography (ANT) on a trenched distributed acoustic sensing (DAS) array | |
CN110737018B (zh) | Vsp地震资料各向异性建模方法 | |
CN109188528B (zh) | 井间弹性波层析成像系统和方法 | |
CN112230274B (zh) | 面向随钻导向的声波方程频率域逆时偏移快速成像方法 | |
CN110244383B (zh) | 基于近地表数据的地质岩性综合模型创建方法 | |
CN108761534B (zh) | 陆上地震加速度信号应用新方法 | |
CN110261899B (zh) | 地震数据z字形干扰波去除方法 | |
CN109324344A (zh) | 基于纯纵波和拟声波反演的页岩厚度预测方法及系统 | |
Jiang et al. | Distributed acoustic sensing for shallow structure imaging using mechanical noise: A case study in Guangzhou, China | |
CN107918152B (zh) | 一种地震相干层析成像方法 | |
CN115236738A (zh) | 一种基于汉克尔变换的勒夫波频散提取方法 | |
Zou et al. | Assessing the reliability of low frequencies in geophone records | |
CN106291675A (zh) | 一种基于基追踪技术的地震数据重构方法 | |
CN108646296A (zh) | 基于自适应谱峭度滤波器的沙漠地震信号噪声消减方法 | |
CN114415234A (zh) | 基于主动源面波频散和h/v确定浅地表横波速度的方法 | |
CN113296149A (zh) | 一种基于微动探测技术的滑坡体地层分布情况探测方法 | |
CN112526611A (zh) | 表层地震波品质因子的提取方法及装置 | |
Wills et al. | Time-lapse measurements of shallow horizontal wave velocity over a compacting field |
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 |