一种潜山精细成像的地震数据处理方法及装置制造方法及图纸

本发明专利技术公开了一种潜山精细成像的地震数据处理方法及装置，该方法包括如下步骤：对潜山地震原始数据进行处理，获得处理后的地震数据；对处理后的地震数据采用精细化气泡处理，以提高提高潜山内幕特征识别；对精细化气泡处理后的地震数据基于三维稀疏taup变换和最优化鬼波参数反演，以实现粗糙海况下的最优化鬼波压制。针对潜山风化壳识别和潜山内幕低频不足的问题，本发明专利技术通过精细化气泡压制技术，有效压制了基底下延续相位；通过三维Taup域自适应鬼波压制技术，有效压制鬼波旁瓣，增强了潜山内幕低频信号和分辨率，更有利于潜山顶面识别和风化壳追踪解释。别和风化壳追踪解释。别和风化壳追踪解释。

【技术实现步骤摘要】
一种潜山精细成像的地震数据处理方法及装置


[0001]本专利技术涉及地震数据处理技术，具体涉及一种潜山精细成像的地震数据处理方法。

技术介绍

[0002]潜山是指现今盆地沉积前就已经存在的古代山体，潜山内幕是指潜山内部的构造形态。由于埋藏较深，构造形态复杂，导致地震数据资料信噪比低，内幕成像不清，潜山成为当前的一个地震勘探难点。

技术实现思路

[0003]本专利技术的目的在于克服上述现有技术的不足，提供一种潜山精细成像的地震数据处理方法及装置，以改善潜山内幕成像效果。
[0004]为实现上述目的，本专利技术的技术方案是：
[0005]第一方面，本专利技术提供一种潜山精细成像的地震数据处理方法，包括如下步骤：
[0006]对潜山地震原始数据进行处理，获得处理后的地震数据；
[0007]对处理后的地震数据采用精细化气泡处理，以提高提高潜山内幕特征识别；
[0008]对精细化气泡处理后的地震数据基于三维稀疏taup变换和最优化鬼波参数反演，以实现粗糙海况下的最优化鬼波压制。
[0009]进一步地，所述精细化气泡处理包括：
[0010]当得到精细气枪子波之后，利用预测反褶积或者直接截断的方法得到一个期望输出子波，然后通过的维纳滤波器，将精细气枪子波子波x(t)整形为期望输出的子波形状d(t)，以获得整形算子f(τ)，将f(τ)应用到地震数据数据压制气泡效应。
[0011]进一步地，所述基于三维稀疏taup变换和最优化鬼波参数反演包括：
[0012]鬼波视为一次波波场经海面反射下行所形成的下行波场，即
[0013]g(p
x
，p
y
，ω)＝R(p
x
，p
y
，ω)u
′
(p
x
，p
y
，ω)e
‑
iωΔt
ꢀꢀꢀ
(1)
[0014]式中，g表示下行的鬼波波场，u
′
表示一次波波场，R为鬼波入射到海平面的反射系数，ω为角频率，总波场u为
[0015]u(p
x
，p
y
，ω)＝u
′
(p
x
，p
y
，ω)+g(p
x
，p
y
，ω)＝
′
(p
x
，p
y
，ω)+R(P
x
,p
y
,ω)u
′
(P
x
,P
y
,ω)e
‑
iωΔt
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(2)
[0016][0017]利用式(2)求出一次波波场,即消除鬼波后的结果
[0018][0019]利用鬼波算子G来描述形成鬼波的波场延拓过程，即
[0020]G(p
x
，p
y
，ω)＝1+R(p
x
，p
y
，ω)e
‑
iωΔt
ꢀꢀꢀ
(4)
[0021]根据式(3)，得τ
‑
p域鬼波压制算法：
[0022][0023]其中u
′
为频率
‑
慢度域一次波波场，ε为保证分母不为0的稳定因子。
[0024]进一步地，所述方法还包括：
[0025]对实现粗糙海况下的最优化鬼波压制后的地震数据分别采用褶积法3DGSRME和模型法3DGSRME进行多次波建模，然后采用曲波域多模型自适应匹配减对多次波进行压制，以获得地震数据中的有效信号。
[0026]进一步地，所述方法还包括：对地震数据中的有效信号进行速度分析后来构建速度模型。
[0027]进一步地，所述速度模型的构建采用多信息约束
‑
多层
‑
多尺度网格层析+扫描层析的方式构建。
[0028]进一步地，所述多信息约束
‑
多层
‑
多尺度网格层析+扫描层析的方式包括如下步骤：
[0029]先基于“地质认识、构造层位、测井信息”建立一个具有背景趋势的大尺度的初始速度模型；
[0030]然后从浅到深逐步采用高精度网格层析进行速度的迭代更新；
[0031]针对深层及基低对速度不敏感的问题，采用速度扫描的方式求取相对准确的速度场。
[0032]进一步地，所述高精度网格层析包括：
[0033]首先进行目标线偏移；然后根据叠加剖面拾取小层、地层倾角、方位角及连续性信息作为辅助控制其次，自动拾取CIP道集剩余延迟，从没有拉平的共成像道集确定走时的错误；最后CIP TOMO利用在CIP道集上自动拾取出的剩余延迟、地层倾角、方位角及连续性信息作为输入和控制信息进行层析反演，成像道集剩余时差层析成像结合叠前深度偏移使得反演结果既受目标的控制，也受成像道集是否被拉平的限制；
[0034]如果速度模型正确，则CRP道集拉平；反之，则CRP道集存在一定的时差。
[0035]第二方面，本专利技术提供一种潜山精细成像的地震数据处理装置，包括存储器、处理器以及存储在所述存储器中并可在所述处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现如上任一所述方法的步骤。
[0036]第三方面，本专利技术提供一种计算机可读存储介质，所述计算机可读存储介质存储有计算机程序，所述计算机程序被处理器执行时实现如上任一所述方法的步骤。
[0037]本专利技术与现有技术相比，其有益效果在于：
[0038](1)针对潜山风化壳识别和潜山内幕低频不足的问题，本专利技术通过精细化气泡压制技术，有效压制了基底下延续相位；通过三维Taup域自适应鬼波压制技术，有效压制鬼波旁瓣，增强了潜山内幕低频信号和分辨率，更有利于潜山顶面识别和风化壳追踪解释；
[0039](2)针对基岩潜山内幕绕射多次波发育，本专利技术采用褶积法和模型法3DSRME联合建模，通过多方法联合建立更准确的多次波模型，继而通过稀疏域联合减有效衰减潜山内幕多次波，提高了潜山内幕信噪比；
[0040](3)针对基岩构造起伏大、内幕信噪比低，本文采用三维构造约束速度扫描和高精
度网格层析反演技术技术，建立了精细可靠的速度模型，有效提高了潜山内幕风化带和裂缝带的成像。
附图说明
[0041]图1为本专利技术实施例1提供的潜山精细成像的地震数据处理方法的整体流程图；
[0042]图2为采用曲波域多模型自适应匹配减对多次波进行压制流程图；
[0043]图3为本专利技术实施例2提供的潜山精细成像的地震数据处理装置的组成示意图。
具体实施方式
[0044]下面结合附图和实施例对本专利技术的技术方案做本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种潜山精细成像的地震数据处理方法，其特征在于，包括如下步骤：对潜山地震原始数据进行处理，获得处理后的地震数据；对处理后的地震数据采用精细化气泡处理，以提高提高潜山内幕特征识别；对精细化气泡处理后的地震数据基于三维稀疏taup变换和最优化鬼波参数反演，以实现粗糙海况下的最优化鬼波压制。2.如权利要求1所述的潜山精细成像的地震数据处理方法，其特征在于，所述精细化气泡处理包括：当得到精细气枪子波之后，利用预测反褶积或者直接截断的方法得到一个期望输出子波，然后通过的维纳滤波器，将精细气枪子波子波x(t)整形为期望输出的子波形状d(t)，以获得整形算子f(τ)，将f(τ)应用到地震数据数据压制气泡效应。3.如权利要求1或2所述的潜山精细成像的地震数据处理方法，其特征在于，所述基于三维稀疏taup变换和最优化鬼波参数反演包括：鬼波视为一次波波场经海面反射下行所形成的下行波场，即g(p
x
，p
y
，ω)＝R(p
x
，p
y
，ω)u
′
(p
x
，p
y
，ω)e
‑
iωΔt
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(1)式中，g表示下行的鬼波波场,u
′
表示一次波波场,R为鬼波入射到海平面的反射系数，ω为角频率，总波场u为u(p
x
,p
y
,ω)＝u
′
(p
x
,p
y
,ω)+g(p
x
,p
y
,ω)＝'(p
x
,p
y
,ω)+R(p
x
,p
y
,ω)u'(p
x
,p
y
,ω)e
‑
iωΔt
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(2)利用式(2)求出一次波波场,即消除鬼波后的结果利用鬼波算子G来描述形成鬼波的波场延拓过程，即G(p
x
，p
y
，ω)＝1+R(p
x
，p
y
，ω)e
‑
i...

【专利技术属性】
技术研发人员：张焱，邓盾，
申请(专利权)人：广州海洋地质调查局，
类型：发明
国别省市：