CN113625337A - 一种极浅水高精度地震资料快速成像方法 - Google Patents
一种极浅水高精度地震资料快速成像方法 Download PDFInfo
- Publication number
- CN113625337A CN113625337A CN202110870298.1A CN202110870298A CN113625337A CN 113625337 A CN113625337 A CN 113625337A CN 202110870298 A CN202110870298 A CN 202110870298A CN 113625337 A CN113625337 A CN 113625337A
- Authority
- CN
- China
- Prior art keywords
- frequency
- amplitude
- data
- time
- correction
- 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
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 title claims abstract description 26
- 238000003384 imaging method Methods 0.000 title claims abstract description 14
- 238000012937 correction Methods 0.000 claims abstract description 60
- 230000003068 static effect Effects 0.000 claims abstract description 36
- 238000000034 method Methods 0.000 claims abstract description 30
- 230000001629 suppression Effects 0.000 claims abstract description 18
- 230000002159 abnormal effect Effects 0.000 claims abstract description 17
- 230000008569 process Effects 0.000 claims abstract description 8
- 238000012545 processing Methods 0.000 claims abstract description 8
- 238000006243 chemical reaction Methods 0.000 claims abstract description 5
- 238000010521 absorption reaction Methods 0.000 claims description 21
- 238000001514 detection method Methods 0.000 claims description 18
- 238000001914 filtration Methods 0.000 claims description 15
- 230000003044 adaptive effect Effects 0.000 claims description 9
- 239000000523 sample Substances 0.000 claims description 9
- 238000000354 decomposition reaction Methods 0.000 claims description 6
- 230000008859 change Effects 0.000 claims description 5
- 238000011410 subtraction method Methods 0.000 claims description 5
- 238000010586 diagram Methods 0.000 claims description 4
- 230000002708 enhancing effect Effects 0.000 claims description 4
- 238000009792 diffusion process Methods 0.000 claims description 3
- 238000007619 statistical method Methods 0.000 claims description 3
- 238000013459 approach Methods 0.000 claims description 2
- 238000005259 measurement Methods 0.000 claims description 2
- 238000005516 engineering process Methods 0.000 abstract description 9
- 238000011835 investigation Methods 0.000 description 3
- 238000003825 pressing Methods 0.000 description 3
- 238000010276 construction Methods 0.000 description 2
- 230000007547 defect Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000003672 processing method Methods 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 1
- 230000007123 defense Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 239000007789 gas Substances 0.000 description 1
- 230000004807 localization Effects 0.000 description 1
- NMJORVOYSJLJGU-UHFFFAOYSA-N methane clathrate Chemical compound C.C.C.C.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O NMJORVOYSJLJGU-UHFFFAOYSA-N 0.000 description 1
- 230000008092 positive effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 239000004576 sand Substances 0.000 description 1
Images
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/282—Application of seismic models, synthetic seismograms
-
- 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/32—Transforming one recording into another or one representation into another
- G01V1/325—Transforming one representation into another
-
- 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/34—Displaying seismic recordings or visualisation of seismic data or attributes
- G01V1/345—Visualisation of seismic data or attributes, e.g. in 3D cubes
-
- 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
-
- 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
- G01V1/362—Effecting static or dynamic corrections; Stacking
-
- 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
- G01V1/364—Seismic filtering
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/20—Trace signal pre-filtering to select, remove or transform specific events or signal components, i.e. trace-in/trace-out
-
- 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
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/50—Corrections or adjustments related to wave propagation
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/50—Corrections or adjustments related to wave propagation
- G01V2210/53—Statics correction, e.g. weathering layer or transformation to a datum
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/70—Other details related to processing
- G01V2210/74—Visualisation of seismic data
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
本发明一种极浅水高精度地震资料快速成像方法,设计了适合于该地震资料数据处理的技术流程,通过基于域转换的异常振幅噪音压制、分频自适应相减的多次波压制、多波束与地震联合的剩余静校正技术、潮汐校正技术、空间振幅补偿与分频信号增强等技术,同时解决了极浅水高分辨率地震数据多次波压制与静校正两个难点问题,获得海底反射真实形态,增强了有效同向轴的连续性,浅部断层、河道等有效地质信息得到精细刻画,最终实现该地震资料的快速成像。
Description
技术领域
本发明属于地震勘探数据处理领域,具体涉及一种极浅水高精度地震资料快速成像方法。
背景技术
极浅水高分辨率地震资料采集是一种基于水声学原理的连续走航式探测浅部地层结构和构造的地球物理方法,因其具有低耗、高效、直观的特点,在探查海底浅部地层结构,了解断裂构造的分布、埋藏古河道、浅层气、海底塌陷和滑坡等地质灾害情况,在航道建设、管道检测等近海工程以及海沙资源调查、天然气水合物调查、冷泉探测等资源调查中得到了广泛的应用。为国民经济可持续发展、海域划界和国防建设提供基础地质资料。
极浅水高分辨率地震资料采集在实际作业中容易受到周围环境的干扰,如涌浪、机械振动、螺旋桨转动以及电力干扰等次生干扰源的影响,测区的地震地质条件及水深条件等也是影响浅地层剖面质量的重要因素。导致很大程度上降低了采集资料的信噪比和分辨率,进而影响对地质结构的解释,因此,必须对原始资料进行有效处理才能获得高质量的数据剖面。
由于极浅水高分辨率地震资料无法获取速度信息,无法利用叠加等技术来压制噪音,在数据处理中可以使用的技术方法有限。目前,对于该资料的处理主要是参照常规地震资料的处理方法,而单独针对极浅水高分辨率地震资料数据处理方法的研究不多。
发明内容
本发明针对现有技术中存在的缺陷,提出一种极浅水高分辨率地震资料快速成像方法,通过异常振幅噪音压制、多次波压制、剩余静校正、潮汐校正、空间振幅补偿、分频信号增强等技术,解决了该地震资料中多次波压制与静校正等问题,最终实现该地震资料的快速成像。
本发明是采用以下的技术方案实现的:一种极浅水高精度地震资料快速成像方法,首先对获得的地震资料进行数据解编和检查,处理时具体包括以下步骤:
步骤A、基于域转换的异常振幅噪音压制:将地震数据从时间域转换到频率域,基于中值滤波压制异常振幅;
步骤B、分频自适应相减的多次波压制:进行海底反射时间拾取,确定海底相关多次波周期,基于自适应相减法来压制海底相关多次波;
步骤C、剩余静校正:用多道同相叠加获取模型道,基于统计方法求取静校正的零校正基准线,各道直接向此零校正基准线对齐,达到静校正目的;
步骤D、潮汐校正:通过计算潮汐值,将GPS坐标校正到震源处坐标,实现潮汐时差校正;
步骤E、空间振幅补偿:在保持相对振幅不变的前提下,通过统计学的方式校正各种因素对振幅的影响,使整个剖面上的能量达到均衡;
步骤F、分频信号增强:对地震信号进行分频段增强处理,压制随机噪音,并增强同相轴的连续性,最后输出成果剖面。
进一步的,所述步骤B具体包括以下步骤:
步骤B1、海底反射时间拾取:在地震剖面上进行海底反射时间拾取,拾取海底反射同相轴的波峰或者波谷,确定海底相关多次波周期;
步骤B2、自适应相减:主通道输入含有希望提取的信号S和噪声n1,参考输入只有噪声n2,通过自适应滤波器的权矢量ω调整其输出y=ωt·x,使输出y在最小均方误差意义下接近主通道噪声n1,通过相减器,将主通道的噪声分量n1抵消掉。
进一步的,所述步骤C中剩余静校正具体采用以下方式:
步骤C1、建立模型道:搜索给定炮检点范围内的初至时间,建立初至时间-炮检距散点图,通过给定炮检距范围和滑动步长逐步拟合此炮检距范围内的初至线性关系,将该曲线分解得对应炮检点范围内的模型道,进而通过给定搜索半径和滑动步长得到全工区范围内的模型道;
步骤C2、地表一致性分解:
(1)统计拟合得到每个共炮点道集的模型道后,第Si炮的第Rj接收点道的模型道与实际道之间的旅行时间差表示为
(2)设每个共炮点道集中所有检波点剩余静校正量之和为零,每个共检波点道集中炮点剩余静校正量之和为零;ΔTi,j还表示为该道的炮点和检波点剩余静校正量之和,即
(3)采用多波束数据得到海底反射时间,用多道同相叠加来获取模型道,当道数N足够大时,各道静校正值之和趋于零,该模型道的峰值时间位置就是零校正基准线,各道直接向此零校正基准线对齐,达到静校正的目的。
进一步的,所述步骤D进行潮汐校正时,包括以下两种方式:
(1)利用SkyFix XP高程数据:当采集数据通过SkyFix XP定位***采集潮汐值时,先对SkyFix XP定位***记录的高程数据进行异常值编辑,并对其进行平滑滤波;然后采用高程拟合法将大地高程转换为海拔高程,计算出地震采集过程中的实时潮汐值,然后进行校正;
(2)利用实测水深资料:当采集过程没有高程数据,只有每炮记录的水深数据,海底为水平界面,计算实测水深的变化量为潮汐变化量,然后进行校正。
进一步的,所述步骤E中进行空间振幅补偿步骤如下:
步骤E1、基于时间域窄带通滤波实现时频分解:设第j道地震记录为xj(t),t为时间,把地震数据在最小频率和高截止频率之间等分多个频带,设频带个数为M,即得到M个窄带通滤波门,其带通滤波分频表达式为:
式中h(t,fk)表示第k个频带的带通滤波算子;xj(t,fk)表示第j到第k频带的分频数据;
步骤E2、根据模型每一个频带从浅至深各时窗平均均方根振幅能量中值的变化,采用最小二乘法,导出大地吸收衰减系数的线性方程,拟合大地吸收衰减曲线;
若振幅的球面扩散与吸收衰减满足时间和频率的指数函数的衰减规律,则第j道第k频带振幅为
式中,A(0)为初始振幅,a(t,fk)为n次多项式,将第k频带大地吸收衰减补偿曲线为
a0(fk),a1(fk)t,...,an(fk)tn为第k个频带不同阶数的吸收衰减系数,n值取决于不同频段球面发散的复杂程度;计算各频带每一时刻相应的大地吸收补偿系数;
步骤E3、将分频数据用大地吸收补偿系数做振幅补偿,即
xj(t,fk)=xj(t,fk)×y(t,fk)
然后对整道数据xj(t,fk)统计,拟合得到空间域振幅补偿系数y(t,fk),并对数据做空间振幅补偿。
与现有技术相比,本发明的优点和积极效果在于:
本发明根据极浅水高分辨率地震资料的特点,设计了适合于该地震资料数据处理的技术流程,通过基于域转换的异常振幅噪音压制、分频自适应相减的多次波压制、多波束与地震联合的剩余静校正技术、潮汐校正技术、空间振幅补偿与分频信号增强等技术,同时解决了极浅水高分辨率地震数据多次波压制与静校正两个难点问题,获得海底反射真实形态,增强了有效同向轴的连续性,浅部断层、河道等有效地质信息得到精细刻画。
附图说明
图1为本发明实施例快速成像方法流程示意图;
图2为本发明实施例异常振幅压制前后剖面示意图,(a)异常振幅压制前剖面;(b)异常振幅压制后剖面;
图3为本发明实施例多次波压制前后剖面示意图;(a)多次波压制前剖面;(b)多次波压制后剖面;
图4为本发明实施例剩余校正前后剖面示意图,(a)剩余静校正前剖面;(b)剩余静校正后剖面;
图5为本发明实施例潮汐校正前后剖面示意图,(a)潮汐校正前剖面;(b)潮汐校正后剖面;
图6为本发明实施例空间振幅补偿前后剖面示意图,(a)空间振幅补偿前剖面;(b)空间振幅补偿后剖面;
图7为本发明实施例信号能量增强前后剖面示意图,(a)信号能量增强前剖面;(b)信号能量增强后剖面。
具体实施方式
为了能够更加清楚地理解本发明的上述目的、特征和优点,下面结合附图及实施例对本发明做进一步说明。在下面的描述中阐述了很多具体细节以便于充分理解本发明,但是,本发明还可以采用不同于在此描述的其他方式来实施,因此,本发明并不限于下面公开的具体实施例。
步骤A:基于域转换的异常振幅噪音压制,该方法是将地震数据从时间域转换到频率域,应用中值滤波来压制异常振幅,设定门槛值,在一定的时窗内与振幅中值差异较大的振幅,对其衰减或用相邻道插值。
中值滤波原理:首先在不同频带上,时窗范围内计算每道振幅的平均能量:
式(1)中,Eftk为在时窗t的第k道的f频率段的振幅能量;Aiftk为在时窗t的第k道的f频率段的第i采样点的振幅能量;nftk为在时窗t的第k道的f频率段的样点数。
设定中值滤波的宽度,计算振幅中值偏差值:
MADftk=medk[|Eftk-medk(Eftk)] (2)
式(2)中,MADftk为在时窗t的第k道的f频率段的中值绝对偏差;medk为在中值滤波宽度内,时窗t频带f的所有k道的振幅中值;
设定门槛值,中值滤波通常有两种方法,即加减法和乘除法,加减法表达方式为:
Tftk=Mftk±TFt·MADftk (3)
式(3)中,Tftk为在时窗t的第k道的f频率段设定的门槛振幅值;Mftk为在时窗t的第k道的f频率段的中值振幅;TFt为时窗t的门槛系数,乘除法表达方式为:
Tftk=Mftk·TFt (4)
Tftk=Mftk/TFt (5)
分析地震资料的异常振幅衰减过程中,首先分析噪声分布特征,频率分布范围,振幅能量与有效波的差异;然后设定时窗长度,起始时间尽量避开强振幅的海底反射波,避免产生假频,时窗之间要有一定的重合;对于设定中值滤波时窗的宽度,确保噪声宽度不超过时窗宽度的1/3;针对噪声存在的频率范围,将数据分解成多个频率段;设定门槛值,对异常振幅进行衰减或用相邻道插值代替。
步骤B:分频自适应相减的多次波压制;该方法首先在进行多次压制前进行海底反射时间拾取;其次,确定海底相关多次波周期;然后采用自适应相减法来压制海底相关多次波。
(1)海底反射时间拾取
海底反射时间拾取主要在地震剖面上进行拾取工作,拾取海底反射同相轴的波峰或者波谷。
(2)自适应相减
自适应相减法主要原理:主通道输入含有希望提取的信号S和噪声n1,参考输入只有噪声n2,并通过自适应滤波器的权矢量ω调整其输出y=ωt·x,使y在最小均方误差意义下最接近主通道噪声n1,这样通过相减器,将主通道的噪声分量n1抵消掉。
步骤C:剩余静校正技术:采用模型曲线法,假设应用基准面静校正量之后的初至时间等于正常的旅行时与剩余校正量之和,模型道初至时间与实际初至时间之间的差为该道的炮点和检波点剩余静校正量之和。假设每个共炮点道集中所有检波点剩余静校正量之和为零,每个共检波点道集中炮点剩余静校正量之和为零,可通过给定区域的初至时间-炮检距散点图拟合相应的模型道,并通过地表一致性分解得到每道的炮点和检波点剩余静校正量。
(1)模型道的建立
初至时间模型道的准确性对模型法初至波剩余静校正量的计算影响至关重要(条件允许的话,可采用多波束数据得到精确的初至时间模型道)。为此,在模型道建立过程中采用多道统计方式做线性拟合。即搜索给定炮检点范围内的初至时间,建立初至时间-炮检距散点图,通过给定炮检距范围和滑动步长逐步拟合此炮检距范围内的初至线性关系(初至时间炮检距),将该曲线分解即可得对应炮检点范围内的模型道;通过给定搜索半径和滑动步长得到全工区范围内的模型道。
(2)地表一致性分解
统计拟合得到每个共炮点道集的模型道后,第Si炮的第Rj接收点道的模型道与实际道之间的旅行时间差即可表示为
按照假设前提,ΔTi,j还可表示为该道的炮点和检波点剩余静校正量之和,即
本方法采用多波束数据可精确得到海底反射时间,然后在用多道同相叠加来获取模型道,当道数N足够大时,各道静校正值之和趋于零,该模型道的峰值时间位置必然就是零校正基准线。所以多道模型道的实质就是用统计的方法求取静校正的零校正基准线,各道直接向此零校正线对齐,就达到了静校正的目的。
步骤D:潮汐校正技术;(1)利用SkyFix XP高程数据。当采集数据通过SkyFix XP定位***采集潮汐值时,先对SkyFix XP定位***记录的高程数据,进行异常值编辑,并对其进行平滑滤波,然后采用高程拟合法将大地高程转换为海拔高程,计算出地震采集过程中的实时潮汐值,然后进行校正。(2)利用实测水深资料。当采集过程没有高程数据,只有每炮记录的水深数据,海底为水平界面,可计算实测水深的变化量为潮汐变化量,然后进行校正。
步骤E:空间振幅补偿;该技术是在保持相对振幅不变的前提下,通过统计学的方式校正各种因素对振幅的影响,使整个剖面上的能量达到均衡。空间域振幅补偿步骤如下:
(1)应用时间域窄带通滤波实现时频分解。窄带通滤波分频的关键是如何选择相等带宽的带通滤波门,使其具有较小的吉普斯响应和频带间的频率泄露。设第j道地震记录为xj(t),t为时间,把地震数据在最小频率和高截止频率之间等分10-20多个频带,设频带个数为M,即可得到M个窄带通滤波门。每个门都是一个等腰梯形,前一个梯形门的右边与后一个梯形门的左边在各自中点相交,滤波门陡度通常为4-8Hz。其带通滤波分频表达式为
式中h(t,fk)表示第k个频带的带通滤波算子;xj(t,fk)表示第j到第k频带的分频数据,带通滤波算子可通过下式求得,即
式中f1,f2,f3,f4表示第1到M频带中某频带的带通滤波的滤波门;Δ表示采样间隔;n表示滤波算子长度;
(2)根据模型每一个频带从浅至深各时窗平均均方根振幅能量中值的变化,用最小二乘法,导出大地吸收衰减系数的线性方程,拟合大地吸收衰减曲线。
若振幅的球面扩散与吸收衰减满足时间和频率的指数函数的衰减规律,则第j道第k频带振幅为
式中,A(0)为初始振幅,a(t,fk)为n次多项式,可以表达为
a(t,fk)=a0(fk)+a1(fk)t+...+an(fk)tn
a0(fk),a1(fk)t,...,an(fk)tn为第k个频带不同阶数的吸收衰减系数,n值取决于不同频段球面发散的复杂程度。从而第k频带大地吸收衰减补偿曲线为
可以计算各频带每一时刻相应的大地吸收补偿系数。
(3)将分频数据用大地吸收补偿系数做振幅补偿,即
xj(t,fk)=xj(t,fk)×y(t,fk)
然后对整道数据xj(t,fk)统计,拟合得到空间域振幅补偿系数y(t,fk),并对数据做空间振幅补偿。
步骤F:分频信号增强:
为了使地震波组特征更清晰,同相轴连续性更好,需要对地震信号进行增强处理,压制随机噪音,并增强同相轴的连续性。
分频信号增强采用的是Hilbert-huang变换方法,其关键是求取固有模态函数的方法-经验模态分解法(Empirical Mode Decomposition,EMD)。EMD方法可以将原始的地震信号分解为多个固有模态函数相加的形式,该方法较好地克服了希尔伯特变换中信号局部化的缺点,并且对于含噪声的地震信号的处理效果显著。
以上所述,仅是本发明的较佳实施例而已,并非是对本发明作其它形式的限制,任何熟悉本专业的技术人员可能利用上述揭示的技术内容加以变更或改型为等同变化的等效实施例应用于其它领域,但是凡是未脱离本发明技术方案内容,依据本发明的技术实质对以上实施例所作的任何简单修改、等同变化与改型,仍属于本发明技术方案的保护范围。
Claims (5)
1.一种极浅水高精度地震资料快速成像方法,首先对获得的地震资料进行数据解编和检查,其特征在于,处理时具体包括以下步骤:
步骤A、基于域转换的异常振幅噪音压制:将地震数据从时间域转换到频率域,基于中值滤波压制异常振幅;
步骤B、分频自适应相减的多次波压制:进行海底反射时间拾取,确定海底相关多次波周期,基于自适应相减法来压制海底相关多次波;
步骤C、剩余静校正:用多道同相叠加获取模型道,基于统计方法求取静校正的零校正基准线,各道直接向此零校正基准线对齐,达到静校正目的;
步骤D、潮汐校正:通过计算潮汐值,将GPS坐标校正到震源处坐标,实现潮汐时差校正;
步骤E、空间振幅补偿:在保持相对振幅不变的前提下,通过统计学的方式校正各种因素对振幅的影响,使整个剖面上的能量达到均衡;
步骤F、分频信号增强:对地震信号进行分频段增强处理,压制随机噪音,并增强同相轴的连续性,最后输出成果剖面。
2.根据权利要求1所述的极浅水高精度地震资料快速成像方法,其特征在于:所述步骤B具体包括以下步骤:
步骤B1、海底反射时间拾取:在地震剖面上进行海底反射时间拾取,拾取海底反射同相轴的波峰或者波谷,确定海底相关多次波周期;
步骤B2、自适应相减:主通道输入含有希望提取的信号S和噪声n1,参考输入只有噪声n2,通过自适应滤波器的权矢量ω调整其输出y=ωt·x,使输出y在最小均方误差意义下接近主通道噪声n1,通过相减器,将主通道的噪声分量n1抵消掉。
3.根据权利要求1所述的极浅水高精度地震资料快速成像方法,其特征在于:所述步骤C中剩余静校正具体采用以下方式:
步骤C1、建立模型道:搜索给定炮检点范围内的初至时间,建立初至时间-炮检距散点图,通过给定炮检距范围和滑动步长逐步拟合此炮检距范围内的初至线性关系,将该曲线分解得对应炮检点范围内的模型道,进而通过给定搜索半径和滑动步长得到全工区范围内的模型道;
步骤C2、地表一致性分解:
(1)统计拟合得到每个共炮点道集的模型道后,第Si炮的第Rj接收点道的模型道与实际道之间的旅行时间差表示为
(2)设每个共炮点道集中所有检波点剩余静校正量之和为零,每个共检波点道集中炮点剩余静校正量之和为零;ΔTi,j还表示为该道的炮点和检波点剩余静校正量之和,即
(3)采用多波束数据得到海底反射时间,用多道同相叠加来获取模型道,当道数N足够大时,各道静校正值之和趋于零,该模型道的峰值时间位置就是零校正基准线,各道直接向此零校正基准线对齐,达到静校正的目的。
4.根据权利要求1所述的极浅水高精度地震资料快速成像方法,其特征在于:所述步骤D进行潮汐校正时,包括以下两种方式:
(1)利用SkyFix XP高程数据:当采集数据通过SkyFix XP定位***采集潮汐值时,先对SkyFix XP定位***记录的高程数据进行异常值编辑,并对其进行平滑滤波;然后采用高程拟合法将大地高程转换为海拔高程,计算出地震采集过程中的实时潮汐值,然后进行校正;
(2)利用实测水深资料:当采集过程没有高程数据,只有每炮记录的水深数据,海底为水平界面,计算实测水深的变化量为潮汐变化量,然后进行校正。
5.根据权利要求1所述的极浅水高精度地震资料快速成像方法,其特征在于:所述步骤E中进行空间振幅补偿步骤如下:
步骤E1、基于时间域窄带通滤波实现时频分解:设第j道地震记录为xj(t),t为时间,把地震数据在最小频率和高截止频率之间等分多个频带,设频带个数为M,即得到M个窄带通滤波门,其带通滤波分频表达式为:
式中h(t,fk)表示第k个频带的带通滤波算子;xj(t,fk)表示第j到第k频带的分频数据;
步骤E2、根据模型每一个频带从浅至深各时窗平均均方根振幅能量中值的变化,采用最小二乘法,导出大地吸收衰减系数的线性方程,拟合大地吸收衰减曲线;
若振幅的球面扩散与吸收衰减满足时间和频率的指数函数的衰减规律,则第j道第k频带振幅为
式中,A(0)为初始振幅,a(t,fk)为n次多项式,将第k频带大地吸收衰减补偿曲线为
a0(fk),a1(fk)t,...,an(fk)tn为第k个频带不同阶数的吸收衰减系数,n值取决于不同频段球面发散的复杂程度;计算各频带每一时刻相应的大地吸收补偿系数;
步骤E3、将分频数据用大地吸收补偿系数做振幅补偿,即
xj(t,fk)=xj(t,fk)×y(t,fk)
然后对整道数据xj(t,fk)统计,拟合得到空间域振幅补偿系数y(t,fk),并对数据做空间振幅补偿。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110870298.1A CN113625337B (zh) | 2021-07-30 | 2021-07-30 | 一种极浅水高精度地震资料快速成像方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110870298.1A CN113625337B (zh) | 2021-07-30 | 2021-07-30 | 一种极浅水高精度地震资料快速成像方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113625337A true CN113625337A (zh) | 2021-11-09 |
CN113625337B CN113625337B (zh) | 2022-10-28 |
Family
ID=78381679
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110870298.1A Active CN113625337B (zh) | 2021-07-30 | 2021-07-30 | 一种极浅水高精度地震资料快速成像方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113625337B (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114167494A (zh) * | 2021-11-29 | 2022-03-11 | 哈尔滨工程大学 | 一种基于集合经验模态分解的海底地震波降噪方法 |
CN114636637A (zh) * | 2022-05-07 | 2022-06-17 | 青岛海洋地质研究所 | 一种悬浮物浓度的原位测定装置及工作方法 |
CN115144899A (zh) * | 2022-06-24 | 2022-10-04 | 中国地质大学(北京) | 崎岖海底obn弹性波联合偏移成像方法和装置 |
CN117233839A (zh) * | 2023-11-10 | 2023-12-15 | 山东科技大学 | 地震数据大地吸收衰减三维空间质控方法、***以及设备 |
CN117908108A (zh) * | 2024-03-20 | 2024-04-19 | 山东省地质矿产勘查开发局第二水文地质工程地质大队(山东省鲁北地质工程勘察院) | 一种海洋地震实时监测*** |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5920828A (en) * | 1997-06-02 | 1999-07-06 | Baker Hughes Incorporated | Quality control seismic data processing system |
US5971095A (en) * | 1996-01-09 | 1999-10-26 | Schlumberger Technology Corporation | Noise filtering method for seismic data |
US20160299244A1 (en) * | 2015-04-07 | 2016-10-13 | Korea Institute Of Geoscience And Mineral Resource | METHOD FOR SWELL EFFECT AND MIS-TIE CORRECTION IN HIGH-RESOLUTION SEISMIC DATA USING MULTI-BEAM Echo SOUNDER DATA |
CN109856680A (zh) * | 2019-03-27 | 2019-06-07 | 中国地质科学院地球物理地球化学勘查研究所 | 一种沿海滩涂区拖曳式浅层地震探测方法 |
CN113093280A (zh) * | 2021-04-07 | 2021-07-09 | 青岛海洋地质研究所 | 基于相干函数控制的虚反射走时电缆等浮校正方法 |
-
2021
- 2021-07-30 CN CN202110870298.1A patent/CN113625337B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5971095A (en) * | 1996-01-09 | 1999-10-26 | Schlumberger Technology Corporation | Noise filtering method for seismic data |
US5920828A (en) * | 1997-06-02 | 1999-07-06 | Baker Hughes Incorporated | Quality control seismic data processing system |
US20160299244A1 (en) * | 2015-04-07 | 2016-10-13 | Korea Institute Of Geoscience And Mineral Resource | METHOD FOR SWELL EFFECT AND MIS-TIE CORRECTION IN HIGH-RESOLUTION SEISMIC DATA USING MULTI-BEAM Echo SOUNDER DATA |
CN109856680A (zh) * | 2019-03-27 | 2019-06-07 | 中国地质科学院地球物理地球化学勘查研究所 | 一种沿海滩涂区拖曳式浅层地震探测方法 |
CN113093280A (zh) * | 2021-04-07 | 2021-07-09 | 青岛海洋地质研究所 | 基于相干函数控制的虚反射走时电缆等浮校正方法 |
Non-Patent Citations (9)
Title |
---|
刘俊,等: "南黄海地震资料叠前去噪技术应用", 《海洋地质与第四纪地质》 * |
刘喜祥,等: "多道统计剩余静校正", 《内蒙古石油化工》 * |
干大勇,等: "提高浅层地震成像品质的处理技术-以川中地区沙溪庙组为例", 《石油地球物理勘探》 * |
潘龙,等: "初至波剩余静校正技术在南安集海地区的应用", 《石油地球物理勘探》 * |
王小杰,等: "南黄海中部浅地层剖面数据处理新进展", 《海洋地质前沿》 * |
王振力,等: "一种基于双通道自适应噪声对消的语音增强法", 《信号处理》 * |
胡斐,等: "利用 SkyFixXP高程对高分辨率三维地震资料进行潮汐校正的方法", 《中国海上油气》 * |
范旭,等: "时频空间域振幅补偿方法及其应用", 《新疆石油地质》 * |
赵维娜等: "海底地震仪浅海广角探测的数据特征与噪声组合压制――以南黄海OBS2016测线为例", 《地球物理学报》 * |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114167494A (zh) * | 2021-11-29 | 2022-03-11 | 哈尔滨工程大学 | 一种基于集合经验模态分解的海底地震波降噪方法 |
CN114636637A (zh) * | 2022-05-07 | 2022-06-17 | 青岛海洋地质研究所 | 一种悬浮物浓度的原位测定装置及工作方法 |
CN114636637B (zh) * | 2022-05-07 | 2023-09-01 | 青岛海洋地质研究所 | 一种悬浮物浓度的原位测定装置及工作方法 |
CN115144899A (zh) * | 2022-06-24 | 2022-10-04 | 中国地质大学(北京) | 崎岖海底obn弹性波联合偏移成像方法和装置 |
CN115144899B (zh) * | 2022-06-24 | 2023-02-17 | 中国地质大学(北京) | 崎岖海底obn弹性波联合偏移成像方法和装置 |
CN117233839A (zh) * | 2023-11-10 | 2023-12-15 | 山东科技大学 | 地震数据大地吸收衰减三维空间质控方法、***以及设备 |
CN117233839B (zh) * | 2023-11-10 | 2024-01-26 | 山东科技大学 | 地震数据大地吸收衰减三维空间质控方法、***以及设备 |
CN117908108A (zh) * | 2024-03-20 | 2024-04-19 | 山东省地质矿产勘查开发局第二水文地质工程地质大队(山东省鲁北地质工程勘察院) | 一种海洋地震实时监测*** |
CN117908108B (zh) * | 2024-03-20 | 2024-05-28 | 山东省地质矿产勘查开发局第二水文地质工程地质大队(山东省鲁北地质工程勘察院) | 一种海洋地震实时监测*** |
Also Published As
Publication number | Publication date |
---|---|
CN113625337B (zh) | 2022-10-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN113625337B (zh) | 一种极浅水高精度地震资料快速成像方法 | |
CN106405651B (zh) | 一种基于测井匹配的全波形反演初始速度模型构建方法 | |
CN109669212B (zh) | 地震数据处理方法、地层品质因子估算方法与装置 | |
CN108983284B (zh) | 一种适用于海上斜缆数据的f-p域鬼波压制方法 | |
CN108196305B (zh) | 一种山地静校正方法 | |
CN112327362B (zh) | 速度域的海底多次波预测与追踪衰减方法 | |
CN104820242B (zh) | 一种面向叠前反演的道集振幅分频补偿方法 | |
US9952341B2 (en) | Systems and methods for aligning a monitor seismic survey with a baseline seismic survey | |
CN104570116A (zh) | 基于地质标志层的时差分析校正方法 | |
CN102323618B (zh) | 基于分数阶傅里叶变换的相干噪声抑制方法 | |
CN106950600A (zh) | 一种近地表散射面波的去除方法 | |
CN111239814B (zh) | 基于同相轴分频追踪平滑的浅剖数据机械干扰压制方法 | |
CN110261899B (zh) | 地震数据z字形干扰波去除方法 | |
CN108919345B (zh) | 一种海底电缆陆检噪声的衰减方法 | |
CN113075732B (zh) | 一种消除高分辨率小多道地震地层异常起伏的方法 | |
CN112200069B (zh) | 时频域谱减法和经验模态分解联合的隧道滤波方法及*** | |
CN111538088B (zh) | 一种海上斜缆波场校正方法 | |
CN109884701B (zh) | 地质体散射角度导向深度成像方法 | |
CN112099090B (zh) | 地震资料视速度域非一致性长波长静校正方法 | |
Yao et al. | Microseismic signal denoising using simple bandpass filtering based on normal time–frequency transform | |
CN105572742A (zh) | 一种确定海水深度的方法和装置 | |
CN116125535B (zh) | 三维vsp成像的方法及装置 | |
CN112394393B (zh) | 一种crp道集数据体重构的方法 | |
CN112327361B (zh) | 基于线性同相轴迭代追踪衰减的倾斜干扰剔除方法 | |
CN117148443B (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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |