专利名称：校正声波测井曲线进行储层预测的方法及设备的制作方法
技术领域：
本发明关于地质勘探技术领域，特别是关于石油天然气地球物理勘探技术领域中基于模型反演进行储层预测的技术，具体的讲是一种校正声波测井曲线进行储层预测的方法及设备。
背景技术：
在石油天然气的勘探过程中，整合地震资料时需要建立储层模型，地震反演是建立储层模型最重要的一种方式，其主要分为叠前地震反演和叠后地震反演两大类。从地震反演的实现思路来看，其大致可分为如下三大类相对波阻抗反演(道积分)、递推反演以及基于模型的反演。其中基于模型的反演优于其它两类方法，它突破了传统意义上的地震分辨率的限制，理论上可得到与测井资料相同的分辨率，是油田开发阶段储层精细描述的关键技术。多解性是测井约束地震反演的固有特性，建立尽可能接近实际地层情况的初始波阻抗模型，是减少反演结果多解性的根本途径。而声波测井曲线是建立初始模型的基础资料和地质解释的基本依据。通常情况下，声波测井受到井孔环境的影响而产生误差，同一口井的不同层段，不同井的同一层段误差大小亦不相同。因此，用于制作初始波阻抗模型的测井资料必须经过环境校正。由于水基泥浆的侵蚀，井筒中的地层常发生泥岩蚀变及井眼垮塌。声波测井是在井筒中作业的，泥浆浸泡导致声波测井曲线失真，影响程度随浸泡时间增加而增强。而长期以来人们一般只对储集层的声波时差质量有兴趣，并研制了相应的校正技术，对泥岩段的声波时差质量关注较少。这对储集层评价影响较小，但在储层预测中，声波是联系测井和地震的桥梁，如曲线畸变严重，将导致储层预测的失败。为了揭示测井的环境影响，研究有效和实用的校正方法，胜利油田曾部署了一批油基泥浆取心井。在油基泥浆条件下，井眼十分规则，无明显扩径现象，此时测得的声波时差值基本代表原状地层的声波时差。更换为水基泥浆后，泥岩段就出现地层明显蚀变、扩径等现象。通过多口井的统计对比研究，总结出了声波、密度校正的地区经验公式。但是油基泥浆钻井成本昂贵，井数极少，另外得到的地区经验公式也难以推广。
发明内容
为了克服现有技术存在的上述问题，本发明提供了一种校正声波测井曲线进行储层预测的方法及设备，利用模拟退火算法，通过对声波测井曲线进行校正，得到了与井旁道更加吻合的合成记录，从而为后续的模型反演降低多解性、提高反演精度提供了重要的基础数据，进一步提高了储层预测的准确性。本发明的目的之一是，提供一种校正声波测井曲线进行储层预测的方法，包括采集井旁道资料、声波测井曲线以及原始密度曲线；根据所述的声波测井曲线得到原始纵波速度曲线；根据所述的井旁道资料、原始纵波速度曲线以及原始密度曲线确定目标函数；采集原始自然电位曲线或原始自然伽玛曲线；根据所述的原始自然电位曲线或原始自然伽玛曲线、原始纵波速度曲线、原始密度曲线以及井旁道资料确定新的目标函数；设定退火温度；根据模拟退火算法、所述的退火温度、所述的目标函数以及新的目标函数确定校正后的目标函数；根据所述校正后的目标函数确定校正后的声波测井曲线；根据所述校正后的声波测井曲线进行储层预测。本发明的目的之一是，提供了一种校正声波测井曲线进行储层预测的设备，包括第一采集装置，用于采集井旁道资料、声波测井曲线以及原始密度曲线；原始纵波速度确定装置，用于根据所述的声波测井曲线得到原始纵波速度曲线；目标函数确定装置，用于根据所述的井旁道资料、原始纵波速度曲线以及原始密度曲线确定目标函数；第二采集装置，用于采集原始自然电位曲线或原始自然伽玛曲线；新目标函数确定装置，用于根据所述的原始自然电位曲线或原始自然伽玛曲线、原始纵波速度曲线、原始密度曲线以及井旁道资料确定新的目标函数；退火温度设定装置，用于设定退火温度；目标函数校正装置，用于根据模拟退火算法、所述的退火温度、所述的目标函数以及新的目标函数确定校正后的目标函数；声波测井曲线校正装置，用于根据所述校正后的目标函数确定校正后的声波测井曲线；储层预测装置，用于根据所述校正后的声波测井曲线进行储层预测。本发明的有益效果在于，利用模拟退火算法，通过对声波测井曲线进行校正，得到了与井旁道更加吻合的合成记录，从而为后续的模型反演降低多解性、提高反演精度提供了重要的基础数据，进一步提高了储层预测的准确性。为让本发明的上述和其他目的、特征和优点能更明显易懂，下文特举较佳实施例，并配合所附图式，作详细说明如下。
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。图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，重新对测井曲线进行校正，直至达到具体实施方式
中要求的校正精度。本具体的实施例中，可按照双曲线下降型冷却进度表更新温度，即按照如下公式更新温度tk = a k*t0 ；其中α取为常数，t0为初始温度，k为叠代次数。在具体的实施方式中，还可采集泥岩段的其他声波测井曲线，对每个测井曲线均按照上述的步骤进行校正，直至校正出精度最高效果最好的测井曲线。本发明提供的校正声波测井曲线进行储层预测的方法，利用模拟退火算法，通过对声波测井进行校正，得到了与井旁道更加吻合的合成记录，从而为后续的模型反演降低多解性、提高反演精度提供了重要的基础数据，进一步提高了储层预测的准确性。
图2为图1中的步骤S103的具体流程图，由图2可知，步骤S103具体包括S201 :根据所述的井旁道资料确定子波。在具体的实施例中，子波为16Hz雷克子波。S202 :根据所述的原始纵波速度曲线以及原始密度曲线确定反射系数。S203 :将所述的子波与所述的反射系数进行褶积得到合成记录，在具体的实施方式中，合成记录可通过诸如syn(i)来表示；S204 :根据所述的合成记录、所述的井旁道资料确定目标函数。在具体的实施方式中，井旁道可用诸如seismic (i)来表示，目标函数可为如下公式
权利要求
1.一种校正声波测井曲线进行储层预测的方法，其特征是，所述的方法包括 采集井旁道资料、声波测井曲线以及原始密度曲线； 根据所述的声波测井曲线得到原始纵波速度曲线； 根据所述的井旁道资料、原始纵波速度曲线以及原始密度曲线确定目标函数； 采集原始自然电位曲线或原始自然伽玛曲线； 根据所述的原始自然电位曲线或原始自然伽玛曲线、原始纵波速度曲线、原始密度曲线以及井旁道资料确定新的目标函数； 设定退火温度； 根据模拟退火算法、所述的退火温度、所述的目标函数以及新的目标函数确定校正后的目标函数； 根据所述校正后的目标函数确定校正后的声波测井曲线； 根据所述校正后的声波测井曲线进行储层预测。
2.根据权利要求1所述的方法，其特征是，根据所述的井旁道资料、原始纵波速度曲线以及原始密度曲线确定目标函数具体包括 根据所述的井旁道资料确定子波； 根据所述的原始纵波速度曲线以及原始密度曲线确定反射系数； 将所述的子波与所述的反射系数进行褶积得到合成记录； 根据所述的合成记录、所述的井旁道资料确定目标函数。
3.根据权利要求2所述的方法，其特征是，所述的目标函数为
4.根据权利要求1所述的方法，其特征是，根据所述的原始自然电位曲线或原始自然伽玛曲线、原始纵波速度曲线、原始密度曲线以及井旁道资料确定新的目标函数具体包括 根据所述的原始自然电位曲线或原始自然伽玛曲线、原始纵波速度曲线以及原始密度曲线确定最大纵波速度以及最大密度； 根据所述的原始纵波速度曲线确定初始纵波速度； 根据所述的原始密度曲线确定初始密度； 根据所述的原始纵波速度曲线、初始纵波速度以及最大纵波速度生成新的纵波速度曲线. 根据所述的原始密度曲线、初始密度以及最大密度生成新的密度曲线； 根据所述的井旁道资料、新的纵波速度曲线以及新的密度曲线生成新的目标函数。
5.根据权利要求4所述的方法，其特征是，根据所述的原始自然电位曲线或原始自然伽玛曲线、原始纵波速度曲线以及原始密度曲线确定最大纵波速度以及最大密度具体包括 根据所述的原始自然电位曲线或原始自然伽玛曲线确定解释井段纯泥岩层的测井值以及纯砂岩层的测井值； 根据所述的解释井段纯泥岩层的测井值、纯砂岩层的测井值确定泥岩段的泥质含量；根据所述的泥岩段的泥质含量确定泥岩段的深度范围； 根据所述的深度范围从所述的原始纵波速度曲线以及原始密度曲线中确定出泥岩段深度范围内的声波测井曲线以及密度测井曲线； 从所述泥岩段深度范围内的声波测井曲线以及密度测井曲线中搜索最大纵波速度以及最大密度。
6.根据权利要求5所述的方法，其特征是，所述的泥岩段的泥质含量通过如下公式进行
7.根据权利要求4所述的方法，其特征是，所述的新的纵波速度曲线通过如下公式得到 v"m = 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.一种校正声波测井曲线进行储层预测的设备，其特征是，所述的设备包括 第一采集装置，用于采集井旁道资料、声波测井曲线以及原始密度曲线； 原始纵波速度确定装置，用于根据所述的声波测井曲线得到原始纵波速度曲线； 目标函数确定装置，用于根据所述的井旁道资料、原始纵波速度曲线以及原始密度曲线确定目标函数； 第二采集装置，用于采集原始自然电位曲线或原始自然伽玛曲线； 新目标函数确定装置，用于根据所述的原始自然电位曲线或原始自然伽玛曲线、原始纵波速度曲线、原始密度曲线以及井旁道资料确定新的目标函数； 退火温度设定装置，用于设定退火温度； 目标函数校正装置，用于根据模拟退火算法、所述的退火温度、所述的目标函数以及新的目标函数确定校正后的目标函数； 声波测井曲线校正装置，用于根据所述校正后的目标函数确定校正后的声波测井曲线. 储层预测装置，用于根据所述校正后的声波测井曲线进行储层预测。
15.根据权利要求14所述的设备，其特征是，所述的目标函数确定装置具体包括 子波确定单元，用于根据所述的井旁道资料确定子波； 反射系数确定单元，用于根据所述的原始纵波速度曲线以及原始密度曲线确定反射系数； 褶积单元，用于将所述的子波与所述的反射系数进行褶积得到合成记录； 目标函数确定单元，用于根据所述的合成记录、所述的井旁道资料确定目标函数。
16.根据权利要求15所述的设备，其特征是，所述的目标函数为
17.根据权利要求14所述的设备，其特征是，新目标函数确定装置具体包括 最大值确定单元，用于根据所述的原始自然电位曲线或原始自然伽玛曲线、原始纵波速度曲线以及原始密度曲线确定最大纵波速度以及最大密度； 初始纵波确定单元，用于根据所述的原始纵波速度曲线确定初始纵波速度； 初始密度确定单元，用于根据所述的原始密度曲线确定初始密度； 新纵波曲线生成单元，用于根据所述的原始纵波速度曲线、初始纵波速度以及最大纵波速度生成新的纵波速度曲线；新密度曲线生成单元，用于根据所述的原始密度曲线、初始密度以及最大密度生成新的密度曲线； 新目标函数生成单元，用于根据所述的井旁道资料、新的纵波速度曲线以及新的密度曲线生成新的目标函数。
18.根据权利要求17所述的设备，其特征是，所述的最大值确定单元具体包括 测井值确定单元，用于根据所述的原始自然电位曲线或原始自然伽玛曲线确定解释井段纯泥岩层的测井值以及纯砂岩层的测井值； 泥质含量确定单元，用于根据所述的解释井段纯泥岩层的测井值、纯砂岩层的测井值确定泥岩段的泥质含量； 深度范围确定单元，用于根据所述的泥岩段的泥质含量确定泥岩段的深度范围；泥岩段曲线确定单元，用于根据所述的深度范围从所述的原始纵波速度曲线以及原始密度曲线中确定出泥岩段深度范围内的声波测井曲线以及密度测井曲线； 搜索单元，用于从所述泥岩段深度范围内的声波测井曲线以及密度测井曲线中搜索最大纵波速度以及最大密度。
19.根据权利要求18所述的设备，其特征是，所述的泥岩段的泥质含量通过如下公式进行
20.根据权利要求17所述的设备，其特征是，所述的新纵波速度曲线生成单元通过如下公式进行vr 二 V,:+弋 其中，Vi为初始纵波速度，e i为一随机数且0〈 e i< Avi, A Vi=Vmax-Vi, Vmax为最大速度，W为新的纵波速度。
21.根据权利要求17所述的设备，其特征是，所述的新密度曲线生成单元通过以下公式进行 Pr =Pi+^l1 其中，P i为初始密度，Hi为一随机数且0〈 rii< A Pi, Api= Pniax-P i; Pniax为最大密度，P广为新的密度。
22.根据权利要求17所述的设备，其特征是，所述的新目标函数生成单元具体包括 子波确定单元，用于根据所述的井旁道资料确定子波； 反射系数确定单元，用于根据所述的新的纵波速度曲线以及新的密度曲线确定新的反射系数； 褶积单元，用于将所述的子波与所述的反射系数进行褶积得到新的合成记录； 新目标函数确定单元，用于根据所述的新的合成记录、所述的井旁道资料确定新的目标函数。
23.根据权利要求14所述的设备，其特征是，所述的目标函数校正装置具体包括概率密度确定单元，用于根据
24.根据权利要求14所述的设备，其特征是，所述的声波测井曲线校正装置具体包括 合成记录确定单元，用于根据所述的校正后的目标函数确定其对应的合成记录； 参数确定单元，用于根据所述的合成记录确定与校正后的目标函数对应的纵波速度曲线以及密度曲线； 测井曲线确定单元，用于根据所述的与校正后的目标函数对应的纵波速度曲线以及密度曲线确定校正后的声波测井曲线。
25.根据权利要求14所述的设备，其特征是，所述的设备还包括退火温度更新单元，用于按照冷却进度表更新退火温度。
26.根据权利要求25所述的设备，其特征是，所述的退火温度更新单元通过如下公式进行 tk = a k*t0 ； 其中ct取为常数，h为初始温度，k为叠代次数。
全文摘要
本发明提供一种校正声波测井曲线进行储层预测的方法及设备，所述的方法包括采集井旁道资料、声波测井曲线以及原始密度曲线；根据井旁道资料、原始纵波速度曲线以及原始密度曲线确定目标函数；采集原始自然电位曲线或原始自然伽玛曲线；根据原始自然电位曲线或原始自然伽玛曲线、原始纵波速度曲线、原始密度曲线以及井旁道资料确定新的目标函数；设定退火温度；根据校正后的目标函数确定校正后的声波测井曲线；根据校正后的声波测井曲线进行储层预测。利用模拟退火算法，通过对声波测井曲线进行校正，得到了与井旁道更加吻合的合成记录，从而为后续的模型反演降低多解性、提高反演精度提供了重要的基础数据，进一步提高了储层预测的准确性。
文档编号G01V1/40GK103018778SQ20121053742
公开日2013年4月3日 申请日期2012年12月12日 优先权日2012年12月12日
发明者苏明军, 倪长宽, 陈启林, 刘化清 申请人:中国石油天然气股份有限公司