校正声波测井曲线进行储层预测的方法及设备技术

本发明专利技术提供一种校正声波测井曲线进行储层预测的方法及设备，所述的方法包括：采集井旁道资料、声波测井曲线以及原始密度曲线；根据井旁道资料、原始纵波速度曲线以及原始密度曲线确定目标函数；采集原始自然电位曲线或原始自然伽玛曲线；根据原始自然电位曲线或原始自然伽玛曲线、原始纵波速度曲线、原始密度曲线以及井旁道资料确定新的目标函数；设定退火温度；根据校正后的目标函数确定校正后的声波测井曲线；根据校正后的声波测井曲线进行储层预测。利用模拟退火算法，通过对声波测井曲线进行校正，得到了与井旁道更加吻合的合成记录，从而为后续的模型反演降低多解性、提高反演精度提供了重要的基础数据，进一步提高了储层预测的准确性。

【技术实现步骤摘要】

本专利技术关于地质勘探
，特别是关于石油天然气地球物理勘探
中基于模型反演进行储层预测的技术，具体的讲是一种校正声波测井曲线进行储层预测的方法及设备。
技术介绍
在石油天然气的勘探过程中，整合地震资料时需要建立储层模型，地震反演是建立储层模型最重要的一种方式，其主要分为叠前地震反演和叠后地震反演两大类。从地震反演的实现思路来看，其大致可分为如下三大类相对波阻抗反演(道积分)、递推反演以及基于模型的反演。其中基于模型的反演优于其它两类方法，它突破了传统意义上的地震分辨率的限制，理论上可得到与测井资料相同的分辨率，是油田开发阶段储层精细描述的关键技术。多解性是测井约束地震反演的固有特性，建立尽可能接近实际地层情况的初始波阻抗模型，是减少反演结果多解性的根本途径。而声波测井曲线是建立初始模型的基础资料和地质解释的基本依据。通常情况下，声波测井受到井孔环境的影响而产生误差，同一口井的不同层段，不同井的同一层段误差大小亦不相同。因此，用于制作初始波阻抗模型的测井资料必须经过环境校正。由于水基泥浆的侵蚀，井筒中的地层常发生泥岩蚀变及井眼垮塌。声波测井是在井筒中作业的，泥浆浸泡导致声波测井曲线失真，影响程度随浸泡时间增加而增强。而长期以来人们一般只对储集层的声波时差质量有兴趣，并研制了相应的校正技术，对泥岩段的声波时差质量关注较少。这对储集层评价影响较小，但在储层预测中，声波是联系测井和地震的桥梁，如曲线畸变严重，将导致储层预测的失败。为了揭示测井的环境影响，研究有效和实用的校正方法，胜利油田曾部署了一批油基泥浆取心井。在油基泥浆条件下，井眼十分规则，无明显扩径现象，此时测得的声波时差值基本代表原状地层的声波时差。更换为水基泥浆后，泥岩段就出现地层明显蚀变、扩径等现象。通过多口井的统计对比研究，总结出了声波、密度校正的地区经验公式。但是油基泥浆钻井成本昂贵，井数极少，另外得到的地区经验公式也难以推广。
技术实现思路
为了克服现有技术存在的上述问题，本专利技术提供了一种校正声波测井曲线进行储层预测的方法及设备，利用模拟退火算法，通过对声波测井曲线进行校正，得到了与井旁道更加吻合的合成记录，从而为后续的模型反演降低多解性、提高反演精度提供了重要的基础数据，进一步提高了储层预测的准确性。本专利技术的目的之一是，提供一种校正声波测井曲线进行储层预测的方法，包括采集井旁道资料、声波测井曲线以及原始密度曲线；根据所述的声波测井曲线得到原始纵波速度曲线；根据所述的井旁道资料、原始纵波速度曲线以及原始密度曲线确定目标函数；采集原始自然电位曲线或原始自然伽玛曲线；根据所述的原始自然电位曲线或原始自然伽玛曲线、原始纵波速度曲线、原始密度曲线以及井旁道资料确定新的目标函数；设定退火温度；根据模拟退火算法、所述的退火温度、所述的目标函数以及新的目标函数确定校正后的目标函数；根据所述校正后的目标函数确定校正后的声波测井曲线；根据所述校正后的声波测井曲线进行储层预测。本专利技术的目的之一是，提供了一种校正声波测井曲线进行储层预测的设备，包括第一采集装置，用于采集井旁道资料、声波测井曲线以及原始密度曲线；原始纵波速度确定装置，用于根据所述的声波测井曲线得到原始纵波速度曲线；目标函数确定装置，用于根据所述的井旁道资料、原始纵波速度曲线以及原始密度曲线确定目标函数；第二采集装置，用于采集原始自然电位曲线或原始自然伽玛曲线；新目标函数确定装置，用于根据所述的原始自然电位曲线或原始自然伽玛曲线、原始纵波速度曲线、原始密度曲线以及井旁道资料确定新的目标函数；退火温度设定装置，用于设定退火温度；目标函数校正装置，用于根据模拟退火算法、所述的退火温度、所述的目标函数以及新的目标函数确定校正后的目标函数；声波测井曲线校正装置，用于根据所述校正后的目标函数确定校正后的声波测井曲线；储层预测装置，用于根据所述校正后的声波测井曲线进行储层预测。本专利技术的有益效果在于，利用模拟退火算法，通过对声波测井曲线进行校正，得到了与井旁道更加吻合的合成记录，从而为后续的模型反演降低多解性、提高反演精度提供了重要的基础数据，进一步提高了储层预测的准确性。为让本专利技术的上述和其他目的、特征和优点能更明显易懂，下文特举较佳实施例，并配合所附图式，作详细说明如下。附图说明为了更清楚地说明本专利技术实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本专利技术的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。图1为本专利技术实施例提供的一种校正声波测井曲线进行储层预测的方法的实施方式一的流程图；图2为图1中的步骤S103的具体流程图；图3为图1中的步骤S105的具体流程图；图4为图3中的步骤S301的具体流程图；图5为图3中的步骤S306的具体流程图；图6为本专利技术实施例提供的一种校正声波测井曲线进行储层预测的设备的实施方式一的结构框图；图7为本专利技术实施例提供的一种校正声波测井曲线进行储层预测的设备中目标函数确定装置的结构框图；图8为本专利技术实施例提供的一种基校正声波测井曲线进行储层预测的设备中的新目标函数确定装置的结构框图；图9为本专利技术实施例提供的一种基校正声波测井曲线进行储层预测的设备中的最大值确定单元的机构框图10为本专利技术实施例提供的一种校正声波测井曲线进行储层预测的设备中的新目标函数生成单元的结构框图；图11为具体实施例中的泥质含量曲线及自然伽玛测井曲线；图12为具体实施例中的原声波曲线与校正后声波曲线；图13为具体实施例中的井旁地震道、原合成记录与校正后合成记录；图14为具体实施例中的原误差曲线与校正后的误差曲线；图15为具体实施例中的校正声波测井曲线进行储层预测的流程图。具体实施方式 下面将结合本专利技术实施例中的附图，对本专利技术实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本专利技术一部分实施例，而不是全部的实施例。基于本专利技术中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本专利技术保护的范围。图1为本专利技术实施例提供的一种校正声波测井曲线进行储层预测的方法的实施方式一的流程图，由图1可知，所述方法包括SlOl :采集井旁道资料、声波测井曲线以及原始密度曲线；S102 :根据所述的声波测井曲线得到原始纵波速度曲线；S103:根据所述的井旁道资料、原始纵波速度曲线以及原始密度曲线确定目标函数；S104 :采集原始自然电位曲线或原始自然伽玛曲线；S105:根据所述的原始自然电位曲线或原始自然伽玛曲线、原始纵波速度曲线、原始密度曲线以及井旁道资料确定新的目标函数；S106 :设定退火温度；S107 :根据模拟退火算法、所述的退火温度、所述的目标函数以及新的目标函数确定校正后的目标函数；S108 :根据所述校正后的目标函数确定校正后的声波测井曲线；S109 :根据所述校正后的声波测井曲线进行储层预测。本专利技术提供的预测方法除了上述步骤外，在具体的实施方式中，根据对测井曲线的校正精度的不同需求，还可包括如下步骤按照冷却进度表更新退火温度。退火温度改变后，可返回执行步骤SlOl至步骤S108，重新对测井曲线进行校正，直至达到具体实施本文档来自技高网...

【技术保护点】
一种校正声波测井曲线进行储层预测的方法，其特征是，所述的方法包括：采集井旁道资料、声波测井曲线以及原始密度曲线；根据所述的声波测井曲线得到原始纵波速度曲线；根据所述的井旁道资料、原始纵波速度曲线以及原始密度曲线确定目标函数；采集原始自然电位曲线或原始自然伽玛曲线；根据所述的原始自然电位曲线或原始自然伽玛曲线、原始纵波速度曲线、原始密度曲线以及井旁道资料确定新的目标函数；设定退火温度；根据模拟退火算法、所述的退火温度、所述的目标函数以及新的目标函数确定校正后的目标函数；根据所述校正后的目标函数确定校正后的声波测井曲线；根据所述校正后的声波测井曲线进行储层预测。

【技术特征摘要】
1.一种校正声波测井曲线进行储层预测的方法，其特征是，所述的方法包括 采集井旁道资料、声波测井曲线以及原始密度曲线； 根据所述的声波测井曲线得到原始纵波速度曲线； 根据所述的井旁道资料、原始纵波速度曲线以及原始密度曲线确定目标函数； 采集原始自然电位曲线或原始自然伽玛曲线； 根据所述的原始自然电位曲线或原始自然伽玛曲线、原始纵波速度曲线、原始密度曲线以及井旁道资料确定新的目标函数； 设定退火温度； 根据模拟退火算法、所述的退火温度、所述的目标函数以及新的目标函数确定校正后的目标函数； 根据所述校正后的目标函数确定校正后的声波测井曲线； 根据所述校正后的声波测井曲线进行储层预测。2.根据权利要求1所述的方法，其特征是，根据所述的井旁道资料、原始纵波速度曲线以及原始密度曲线确定目标函数具体包括 根据所述的井旁道资料确定子波； 根据所述的原始纵波速度曲线以及原始密度曲线确定反射系数； 将所述的子波与所述的反射系数进行褶积得到合成记录； 根据所述的合成记录、所述的井旁道资料确定目标函数。3.根据权利要求2所述的方法，其特征是，所述的目标函数为4.根据权利要求1所述的方法，其特征是，根据所述的原始自然电位曲线或原始自然伽玛曲线、原始纵波速度曲线、原始密度曲线以及井旁道资料确定新的目标函数具体包括 根据所述的原始自然电位曲线或原始自然伽玛曲线、原始纵波速度曲线以及原始密度曲线确定最大纵波速度以及最大密度； 根据所述的原始纵波速度曲线确定初始纵波速度； 根据所述的原始密度曲线确定初始密度； 根据所述的原始纵波速度曲线、初始纵波速度以及最大纵波速度生成新的纵波速度曲线. 根据所述的原始密度曲线、初始密度以及最大密度生成新的密度曲线； 根据所述的井旁道资料、新的纵波速度曲线以及新的密度曲线生成新的目标函数。5.根据权利要求4所述的方法，其特征是，根据所述的原始自然电位曲线或原始自然伽玛曲线、原始纵波速度曲线以及原始密度曲线确定最大纵波速度以及最大密度具体包括 根据所述的原始自然电位曲线或原始自然伽玛曲线确定解释井段纯泥岩层的测井值以及纯砂岩层的测井值； 根据所述的解释井段纯泥岩层的测井值、纯砂岩层的测井值确定泥岩段的泥质含量；根据所述的泥岩段的泥质含量确定泥岩段的深度范围； 根据所述的深度范围从所述的原始纵波速度曲线以及原始密度曲线中确定出泥岩段深度范围内的声波测井曲线以及密度测井曲线； 从所述泥岩段深度范围内的声波测井曲线以及密度测井曲线中搜索最大纵波速度以及最大密度。6.根据权利要求5所述的方法，其特征是，所述的泥岩段的泥质含量通过如下公式进行7.根据权利要求4所述的方法，其特征是，所述的新的纵波速度曲线通过如下公式得到 vm = Vf + Si 其中，Vi为初始纵波速度，e i为一随机数且0〈 e i< Avi, A Vi=Vmax-Vi, Vmax为最大速度，Vrj为新的纵波速度。8.根据权利要求4所述的方法，其特征是，所述的新的密度曲线通过以下公式得到 P 二 P +>1, 其中，Pi为初始密度，Hi为一随机数且0〈 ni〈 A p i，APi=Pniax-Pi, PniaxS最大密度P力新的密度。9.根据权利要求4所述的方法，其特征是，根据所述的井旁道资料、新的纵波速度曲线以及新的密度曲线生成新的目标函数具体包括 根据所述的井旁道资料确定子波； 根据所述的新的纵波速度曲线以及新的密度曲线确定新的反射系数； 将所述的子波与所述的反射系数进行褶积得到新的合成记录； 根据所述的新的合成记录、所述的井旁道资料确定新的目标函数。10.根据权利要求1所述的方法，其特征是，根据模拟退火算法、所述的退火温度、所述的目标函数以及新的目标函数确定校正后的目标函数具体包括 根据11.根据权利要求1所述的方法，其特征是，根据所述校正后的目标函数确定校正后的声波测井曲线具体包括 根据所述的校正后的目标函数确定其对应的合成记录； 根据所述的合成记录确定与校正后的目标函数对应的纵波速度曲线以及密度曲线；根据所述的与校正后的目标函数对应的纵波速度曲线以及密度曲线确定校正后的声波测井曲线。12.根据权利要求1所述的方法，其特征是，所述的方法还包括按照冷却进度表更新退火温度。13.根据权利要求12所述的方法，其特征是，按照冷却进度表更新退火温度通过如下公式进行 tk= a k*t0 ； 其中ct取为常数，h为初始温度，k为叠代次数。14.一种校正声波测井曲线进行储层预测的设备，其特征是，所述的设备包括 第一采集装置，用于采集井旁道资料、声波测井曲线以及原始密度曲线； 原始纵波速度确定装置，用于根据所述的声波测井曲线得到...

【专利技术属性】
技术研发人员：苏明军，倪长宽，陈启林，刘化清，
申请(专利权)人：中国石油天然气股份有限公司，
类型：发明
国别省市：