[发明专利]一种基于局部峰值频率的陆相薄储层解释方法及系统

说明书

本发明公开了一种基于局部峰值频率的陆相薄储层解释方法及系统，属于地震勘探技术领域，解决了对薄层进行精细刻画较难的问题。一种基于局部峰值频率的陆相薄储层解释方法，包括以下步骤：将地震数据进行Q补偿，得到经Q补偿后的地震数据；对所述经Q补偿后的地震数据进行广义S变换，得到时频数据体；获取所述时频数据体的局部时频峰值所有对应的时刻，由局部时频峰值所有对应的时刻确定薄层位置。实现提高了薄层的解释精度。

技术领域

本发明涉及地震勘探技术领域，尤其是涉及一种基于局部峰值频率的陆相薄储层解释方法及系统。

背景技术

随油田开发的时间的变化，以及数据量增多，这些地震解释成果与地质认识不一致，以及地质认识与实际工程开发结果存在矛盾，很显然，要提供可靠的地震解释成果需要地震解释人员结合开发阶段的各类数据进行综合分析；

而陆相湖盆烃类储层以碎屑岩占绝对优势，储层非均质性强，尽管开发阶段井网密度很高，但仅靠井资料对井间关系进行推断仍非常不确定，而地震信息虽具备横向可追踪的特性，但常规地震资料的分辨率低，难以用于对薄储层(薄层)进行精细刻画。

发明内容

本发明的目的在于至少克服上述一种技术不足，提出一种基于局部峰值频率的陆相薄储层解释方法及系统。

一方面，本发明提供一种基于局部峰值频率的陆相薄储层解释方法，包括以下步骤：

采集地震数据，将地震数据进行Q补偿，得到经Q补偿后的地震数据；

对所述经Q补偿后的地震数据进行广义S变换，得到时频数据体；

获取所述时频数据体的局部时频峰值所有对应的时刻，由局部时频峰值所有对应的时刻确定薄层位置。

进一步地，所述将地震数据进行Q补偿，得到经Q补偿后的地震数据，具体包括：

利用公式对地震数据进行Q补偿，得到经Q补偿后的地震数据S_Q(t)，其中，r(τ)为反射系数序列，a_Q(τ，f)＝e^{-πfτQ+iH(πfτ/Q)}为Q的传递函数，H为在τ时刻对频率f的希尔伯特变换。

进一步地，对所述经Q补偿后的地震数据进行广义S变换，得到时频数据体，具体包括：

利用公式对所述经Q补偿后的地震数据进行广义S变换，得到时频数据体ST(τ,f)，其中，λ和q是控制广义S变换时频分辨率的两个参数，λ>1，q>1。

进一步地，获取所述时频数据体的局部时频峰值所有对应的时刻，具体包括：

求取时频数据体ST(τ，f)中频率的极大值，此时对应的时刻τ则为时频数据体的局部时频峰值所有对应的时刻。

另一方面，本发明还提供一种基于局部峰值频率的陆相薄储层解释系统，包括地震数据采集模块、补偿模块、时频数据体获取模块、薄层位置确定模块；

地震数据采集模块用于采集地震数据；

补偿模块用于将地震数据进行Q补偿，得到经Q补偿后的地震数据；

时频数据体获取模块用于对所述经Q补偿后的地震数据进行广义S变换，得到时频数据体；

薄层位置确定模块用于获取所述时频数据体的局部时频峰值所有对应的时刻，由局部时频峰值所有对应的时刻确定薄层位置。

进一步地，时频数据体获取模块利用公式对地震数据进行Q补偿，得到经Q补偿后的地震数据S_Q(t)，其中，r(τ)为反射系数序列，a_Q(τ，f)＝e^{-πfτQ+iH(πfτ/Q)}为Q的传递函数，H为在τ时刻对频率f的希尔伯特变换。