基于地质模型的映射采样地震数据体时深转换方法及应用

[收稿日期]20190320

[基金项目]中国石油天然气集团有限公司 十三五 重大科技专项(2016E -0213)㊂ [作者简介]白振华(1982),男,硕士,工程师,现主要从事地球物理与油气地质方面的综合研究工作,u p

c b z h @163.c o m ㊂[引著格式]白振华.基于地质模型的映射采样地震数据体时深转换方法及应用[J ].长江大学学报(自然科学版),2019,16(10

):23~28.

基于地质模型的映射采样地震数据体时深转换方法及应用

白振华

(中石油大庆油田有限责任公司勘探开发研究院,黑龙江大庆163712

)[摘要]在油田开发中后期的精细油藏描述阶段,地震数据体对油藏储层模拟的约束和指导有至关重要的作用,但基于传统速度模型的时深转换存在较大不确定性,因此,提出了基于地质模型的映射采样法时深转换,规避了速度场的累计误差和不确定性㊂首先,在区域标志层面上实现井震精确标定,追踪解释成图时深转换后,分时间域和深度域同步建立高精度的地质模型;然后,通过地震重采样方法把地震数据体采样到时间域模型;最后,通过映射采样实现时间域模型到深度域模型的转换㊂与传统方法相比,该方法完全保证了转换的精度和保真度,提高了油田开发地质与油藏描述的精度,应用效果较好㊂[关键词]速度场;时间域;深度域;时深转换;映射采样[中图分类号]P 631.44

[文献标志码]A [文章编号]16731409(2019)10002306

随着地质统计学建模技术的发展,储层模拟技术在油田开发阶段越来越受到重视,但受井网密度㊁井位空间分布和储层非均质性的综合影响,单靠井上硬数据的地质统计学规律很难实现储层空间描述的准确性㊂地震数据平面上采集密度大,在横向上具有较好的连续性和预测性,可以对储层横向分布起到

较好的趋势约束作用[

1,2

],但是地震信息在垂向上分辨率低,面临从时间域到深度域的转换精度问题㊂在层面构造进行时深转换已形成了多种成熟的方法,在二维层面上的转换可以实现高精度,然而这些方

法并不适用于三维数据体的精细时深转换[

3~5]

㊂因此,三维精细速度场的建立和高精度的时深转换是研究人员需要解决的一个重要且具有挑战性的问题[

6~8]

㊂一个可靠的速度场需要结合所有可用的㊁最佳的来自地震和井的速度信息,根据适当的分层方案,

考虑岩性㊁地质条件变化和地层内各向异性的影响,以地震速度分析为基础㊁井点速度为标尺,建立一个包含不同层速度的速度模型㊂对于地震速度的精确分析求取,研究人员进行了长期大量的研究工作,

初期主要集中于传统的叠加速度分析,应用D i x 公式方法求取层速度[9,1

0]

,但该方法仅适用于水平层状㊁横向无变速介质㊂为解决速度横向变化问题,学者们研究了基于贝叶斯理论的速度反演分析,在一定程度上解决了该问题,但仍是基于水平层状介质的假设,倾斜地层对速度分析精度影响较大㊂20世纪中期发展起来的偏移地震成像技术,促进了多种以地震成像最优的偏移速度分析方法的形成,在存在横向速度变化的情况下,偏移成像相当于求解描述成像射线运动学和几何扩散的偏微分方程,其过程就是找到一个精确的区间速度㊂偏移速度具有更高的分析精度,但依然不能满足横向速度变化剧烈㊁构造复杂的成像精度㊂继而提出了基于地震波运动学和动力学分析的层析成像方法,进一步提高了成像速度

分析精度[

11,12]

,但受地震资料品质等多条件影响,速度分析始终存在较大的不确定性㊂无论哪种速度场建立方法,其空间不确定性都始终存在㊂为此,笔者在伊拉克哈法亚油田M 层油

藏开展了基于地质模型的映射采样地震数据体时深转换方法研究,以规避速度场空间变化的不确定性,提高时深转换精度,充分发挥地震数据在油田开发中后期对储层精细描述的约束和指导作用㊂

㊃

32㊃长江大学学报(自然科学版) 2019年第16卷第10期

J o u r n a l o fY a n g t z eU n i v e r s i t y (N

a t u r a l S c i e n c eE d i t i o n ) 2019,V o l .16N o .10

㊃42㊃地质资源与地质工程2019年10月1传统速度场的不确定性

建立速度场的数据源主要有3个:地震数据的速度㊁声波测井资料得到的速度和V S P(垂直地震剖面)井间地震得到的速度㊂

由地震数据得到的速度,不论处理解释得到的叠加速度谱㊁偏移速度或者反演得到的层速度及其他类型的速度资料,都需要经过修饰校正处理,即以地震解释层位为约束,利用D i x公式将叠加速度转换为层速度和平均速度㊂由于D i x公式是在水平层状介质中射线垂直入射条件下建立起来的速度关系式,前提条件非常苛刻,当地层产状复杂时,用其求取的层速度或平均速度会产生较大误差[13,14],横向连续性好但纵向分辨率不够㊂

由声波测井资料得到的速度一般通过拟合公式㊁合成记录曲线回归㊁克里金插值法㊁随机模拟法等[15~17]建立速度场,其核心是依据井上速度按照趋势插值计算,建立的速度场在近井区精度较高,但在远井区速度不确定性增加㊂

V S P井间地震得到的平均速度和层速度是比较可靠的地层真实速度,但是V S P资料成本高,信息量少,一般达不到单独建立速度场的要求㊂

任何一种方法求取㊁建立速度场都受到一定地质条件的限制,建立的速度场无法兼顾全区不同地质条件变化的影响,不能避免速度场在井间的不确定性和由上至下的误差累计,难以满足三维数据体精细时深转换的要求㊂

2映射采样时深转换方法

要实现映射采样时深转换有2个关键点:一是区域稳定标志层的精确井震标定,得到可靠的时间域构造图,并实现层面的高精度时深转换;二是时间域和深度域地质模型网格建立的完全一致,保证时间和深度域实现一一对应的映射关系㊂

2.1层位标定解释与时深转换

综合录井㊁测井和区域地质研究,找到稳定的区域标志层,通过多井地震合成记录标定,建立地震记录与地质分层的时深联系,保证单井的合成记录与地震反射的相关性;在该基础上进行层面解释,生成时间域层位构造图,根据全区合成记录和V S P资料校正建立的时深关系生成深度域层位构造图㊂2.2网格模型与时深转换

2.2.1深度域网格模型

基于P e t r e l建模软件,根据时间与深度构造图分别建立地质网格模型㊂首先选取建模工区范围,根据地震采集面元㊁井密度和构造方向设置平面网格面元大小和I㊁J方向(见图1);再根据地震解释深度域构造图建立地层框架模型(H o r i z o n);然后在单井地质分层与地震面趋势约束下,建立各层内的小层组地层网格(Z o n e);最后根据小层组地层厚度和垂向分辨率需要等比例设置小层组的细分网格数(L a y e r)㊂深度域模型完成后对网格质量进行质量检查,保证网格没有负体积和扭曲㊂

2.2.2时间域网格模型

时间域地质模型建立步骤与深度域模型保持一致,使用同一平面网格系统,根据时间域构造图建立地层框架模型,在建立时间域小层组地层网格时,采用I s o c h o r e厚度约束方式由深度域模型的小层组厚度做约束,保证小层组在时间和深度域分布一致;然后对小层组进行等比例细分网格,细分数目保持与深度域一致㊂由此建立的深度域和时间域模型在平面和垂向有一致的映射关系,搭建了时间域到深度域转换的物理指数通道(见图1)㊂

2.2.3映射采样时深转换

通过S e i s m i c r e s a m p l i n g地震重采样模块把地震资料采样到时间域的网格中,该过程为时间域地震数据到时间域模型网格采样,数据完全保真㊂因为映射采样只能在同一域下进行,所以在地震采样完成