摘要：维地质建模工作的最终目标是建立尽可能精确的储层参数模型来展示三维空间的地质情况，并为油气藏数值模拟提供地质基础。震控储层建模技术是以地质统计学和变差函数分析等常规数学建模方法为基础，通过将合理的地震参数数据作为约束，实现井间测井数据更合理内插的一种储层建模方法。该技术可以将含有丰富地质、岩性、物性信息的地震数据和测井数据结合在一起，既体现了测井数据的垂向分辨率，又综合了地震数据体反映的储层横向变化特征，特别适用于井网密度不高、储层非均质性强的油气田。为此，基于该方法在普光气田的实际应用效果，结合地震反演数据在储层建模中的应用成果，把地震数据作为宏观约束条件，有效地解决了井网密度不高时储层属性模型的精度问题，并指出了储层建模中整合地震数据方法的若干发展趋势。
关键词：普光气田；震控地质建模；井震结合地震反演；多级建模储集层
1 震控建模研究方法
所谓震控建模，就是在利用随机模拟方法对测井数据进行插值计算的同时，通过协克里金方法，将地震数据作为第二变量，对模拟计算进行加权和条件约束，使测井数值的插值与地震数据的分布特征相近似。利用地震数据进行约束，可实现确定性建模和随机性建模相结合，充分发挥地震储层横向预测的优势，使所建模型更符合地质实际[1]。
具体地说，震控建模首先根据对比和标定，分析地震属性参数与储层物性参数之间的关系，寻找两者之间对应性最好的参数对。通过选择合理的计算方法，将地震数据引入地质模型的模拟计算中。由于地震数据是采集得到的原始数据，比较忠实地反映了地下地层的各种地质特征，因而可以使地质模型的模拟计算更为接近地下的实际地质情况。但地震数据也存在着分辨率低、多解性强的问题。因此这种方法适应于钻井密度低、地震数据质量好、处于开发阶段早期的油气藏。
震控建模的核心是寻找出一种与储层物性参数有相关性的地震属性参数，例如在普光气田，波阻抗与储层的孔隙度之间具有明显的负相关性，即低阻抗区往往与高孔隙度区相对应，而高阻抗区往往与低孔隙度区相对应。然后，在测井曲线插值计算时借鉴前期储层反演预测得到的地震波阻抗数据体的空间分布规律，对测井曲线的插值进行加权或约束[2]。
需要说明的是，震控并不是将地震数据直接转换为属性模型，模型的计算还是以测井曲线的插值计算为基础，只是在插值过程中参考地震属性在空间上的变化规律。因此在井点处属性模型还是与测井资料完全一致，而在井间或无井控制区，使测井数据插值计算出的数据分布规律(例如等值线的形态)与地震属性的数据分布规律(或等值线的形态)具有相似性。
2 震控建模技术流程及步骤
目前将地震数据引入到测井数据插值计算中进行条件约束的数学方法主要是协克里金法。在协克里金方法中将插值计算的数值(例如孔隙度)作为主变量，将参考地震数据(例如地震波阻抗)作为次级变量。首先分别统计两个变量的数据分布函数，然后建立起二者之间的关联，再利用次级变量的分布函数对主变量的插值计算进行加权和约束。在计算中还可以通过设定相关系数来控制地震数据对测井数据插值的影响程度。利用地震数据约束地质模型的插值计算进行震控地质建模，主要有以下几个重要的步骤。
1) 对地震数据进行多种特殊处理(如反演计算)，得到各种属性参数。
2) 进行各种地震属性与各种储层参数之间的相关性分析。通过标定、回归、统计等方法，寻找出与储层地质参数有相关性的地震属性参数。例如在碎屑岩地区波阻抗与泥质含量、砂体厚度的关系；碳酸盐岩地区波阻抗与储层孔隙度的关系；波阻抗、地层倾角与储层裂缝的关系；地震多参数聚类与沉积相带的关系等。这一步骤是震控建模的关键，是必不可少的一步。其主要作用是防止地震属性的“错用”或“滥用”。
3) 根据三维地质模型的网格对地震数据体进行重采样，充填到三维地质模型网格空间内，形成地震属性模型，为以后的模拟计算提供基础。
4) 确定第一个属性模型。计算的第一个模型必须是和地震属性关系最为密切的储层参数，再以这一模型为基础计算其他模型。在碎屑岩地区，地震属性往往与砂体分布、泥质含量等参数有关。因此可以先建立砂体分布模型(岩相模型)或泥质含量模型，然后再通过相控的方法利用这些模型对孔隙度、渗透率等物性参数的模型进行控制。在普光地区波阻抗体直接与孔隙度有较好的相关性，而且孔隙发育区呈连片分布。另外相关性分析也没有发现哪种地震属性与沉积相、岩性等参数具有相关性，因而可以跳过碎屑岩建模常用的相控建模的方法，直接利用波阻抗体作为约束条件建立孔隙度模型，不需首先建立相模型。
5) 进行震控模型的计算。利用测井数据进行模型的插值计算，同时通过协克里金的方法利用地震属性进行条件约束，即在插值计算井间或无井区某一点或区域的数值时，首先分析同一区域地震数据的分布特征和变化趋势，然后参考地震数据的变化趋势决定测井数据在井间的插值趋势，从而达到利用地震数据约束(或控制)测井数据插值(震控)的目标。
6) 计算其他地质模型。在建立起第一个属性模型后，再利用第一个模型为约束条件(也可以同时再加上地震属性)计算其他地质模型。例如在普光地区首先建立孔隙度模型，然后再利用孔隙度模型作为约束条件，计算渗透率模型[3]。
若地震数据与岩相、沉积相之间的相关性不够理想，则以地震数据作为约束条件对能够反映相变化的地球物理特征曲线进行随机模拟，建立地球物理特征曲线模型，再以地球物理特征曲线模型为约束进行相建模，并以有限的相模型和地球物理特征曲线模型为约束进行储层参数建模。这种应用多来源信息、逐级约束、多级建模的方法比较适用于多砂体叠置、单砂体成藏、相变快、规律性差的油气田[4～5]。
3 震控建模在普光气田的应用实例
普光地区储层属于碳酸盐岩礁滩相沉积，针对储层储集空间的类型多样，储层厚度大、非均质性较强，存在多套气水关系等地质特点，综合利用地震、测井及地质等资料，在对储层主控因素分析的基础上，根据确定性建模与随机建模相结合的原则、等时建模原则、地震数据约束建模原则，利用序贯高斯模拟法、克里金模拟法，对储层物性空间分布的非均质性进行了细致的描述，建立普光气田主体构造储层模型。
3.1 构造模型建立
构造模型是三维地质建模工作的基础，构造模型的可靠性和细致程度直接决定了三维地质模型对储集单元描述的准确性。在建模工作中通过精细i维地震解释、地震数据与构造模型的三维可视化交互编辑等技术努力提高构造模型的准确性，可以使构造模型比较真实地反映地下地质体的整体构造形态和特征。
建立构造模型，首先要建立用以后面时深转化的三维速度场。采用了波形对比与平面趋势对比相结合的方法，先制作普光气田主体地区50口井合成记录，搭建了钻井数据和地震数据之间最重要的桥梁。在Petrel软件提供的Layereake变速转换方法，根据合成地震记录提供的时深关系，计算出T1f4、T1f3、T1f1-2和P2ch共4个层系的平均速度平面图，再将各套地层的平均速度平面图进行叠合，得到三维速度模型。在叠合过程中，进一步利用钻井分层数据进行校正。通过三维速度场可以将三维地震数据体转换到深度域，为下一步的属性建模工作打基础。
普光主体工区共发育了3条大型断层，即普光7、普光3、老君庙断层。先利用地震解释成果在时间域内建立了3条断层的断层模型，再利用上面建立的速度场进行时间域到时深度域的转换。由于地震数据的分辨率和精度较低，还需要利用钻井数据对断面进行校正，这样就建立了工区的断层模型。
从地震解释层位中可以得到T1f4顶、T1f3底、T1f1-2底、P2ch底4个界面。因此，在时间域内先利用这些层面直接网格化计算地层界面。将断层与层面模型进行组合，建立地层的空间格架，并进行网格化，最终建立了地层构造模型(图1)。
3.2 地震属性参数优选
地震资料的优势在于能够提供远离井点的储集层信息.但该数据体是否能用于约束储集层参数建模，还需要进一步分析其与储集层物性参数的相关性，只有两者相关性好，方可用地震资料反演约束建模。选择何种地震属性作为约束条件，以及利用地震数据体对何种储层参数做约束是这种建模方法中很重要的一个环节。通过分析认为，孔隙度曲线与地震反演波阻抗之间存在良好的负相关(图2)，即高孔隙度层段通常都对应于低波阻抗反射波，而低孔隙度层对应于高波阻抗。
3.3 储层属性模型建立
在建立孔隙度模型时，利用上面的分析结果，充分发挥地震储集层横向预测的优势，利用地震波阻抗预测成果作为约束条件，实现确定性建模与随机性建模相结合，可以使所建属性模型更符合地质实际。
在孔隙度模型的模拟计算中采用了随机模拟中最常用的序贯高斯模拟方法。模拟计算的同时用地震波阻抗体进行约束，即将地震波阻抗体通过速度模型转换到深度域，按照三维地质模型的网格进行重新采样，利用波阻抗体作为第二变量通过协克里金方法对孔隙度模型的计算进行体对体的条件约束，使波阻抗反演体对模拟计算起到有效的趋势控制作用，定量、直观地表达不同层系储集层孔隙度的空间变化(图3)。
普光气田储集空间类型多样，以孔隙为主，局部发育裂缝。孔隙度和渗透率的关系受多种因素的控制，相关性不好。为此，在建立渗透率模型时分两步进行：①由于储层孔隙类型不同导致孔隙度与渗透率关系多样化，在计算渗透率模型时，首先通过单井测井曲线声波与密度的相对幅度差判断孔隙结构类型，按照不同孔隙结构类型的孔渗关系式，由各井测井孔隙度计算出单井各层面不同孔隙类型的渗透率；②由于裂缝的存在，孔隙度变化不大，而渗透率却成倍增大，因而在对渗透率模型计算时引入裂缝指数数据体对渗透率进行裂缝校正。
裂缝指数模型是通过裂缝研究成果(单井裂缝指数)这一指标，利用序贯指数计算方法建立的裂缝分布模型(图4)。
当渗透率模型用孔隙度计算出来以后，再用裂缝指数模型作为第2个变量加权校正，使渗透率模型比较好地反映裂缝的影响。从最终渗透率模型(图5)来看，其整体的数据分布特征与孔隙度模型有一定的相似性，同时又反映了各层面不同孔隙结构类型的渗透率变化特征，在预测有裂缝发育的部位表现出高渗区沿裂缝呈条带状延伸的特点。
4 结论
1) 随着地震资料处理、解释精度的提高，时间、空间、频率各种域内数据综合分析处理，根据差异互补原则，利用各种域的相对优势，处理其他域内难以解决的描述和预测问题。地震数据逐渐成为储层模拟的重要约束条件之一。
2) 地震数据与其他多种数据(如地质、测井、试井以及动态数据等)进一步融合，多来源、不同精度数据互补优势越来越明显，使传统静态地质模型表征逐渐转向动态、跟踪建模。
3) 利用非线性、随机模拟技术整合地震资料必将成为主流技术。
4) 震控、相控多级建模，既可以发挥地震数据横向分辨率高的优势，也可以用相控来减少地震数据的多解性，理论上比单一的相控、震控更有优势。
参考文献
[1] 张旺青，刘阳平，程超，等.储层相控随机建模研究[J].断块油气田，2008，15(5)：44-46.
[2] 于兴河，李胜利，赵舒，等.河流相油气储层的井震结合相控随机建模约束方法[J].地学前缘，2008，15(4)：33-41.
[3] 朱怡翔，张明禄，王根久，等.苏里格气田相控建模及有利储层预测[J].中国石油勘探，2004，9(1)：52-58.
[4] 潘少伟，杨少春，杨柏，等.相控建模技术在江苏油田庄2断块中的应用[J].天然气地球科学，2009，20(6)：935-940.
[5] 严申斌，李少华，邓恒.相控储层建模在胜南油田的应用[J].断块油气田，2008，15(1)：16-18.
(本文作者：姜贻伟1 刘红磊2 杨福涛1 张纪喜2 张雪松2 1.中国石化中原油田分公司；2.中国石化中原油田勘探开发科学研究院)
您可以选择一种方式赞助本站
支付宝转账赞助
微信转账赞助
- 注解：本资料由会员及群友提供仅供阅读交流学习，不得用于商业用途！
