一种流体替换的低渗透油藏水淹信息的计算方法
【技术领域】
[0001] 本发明属于石油测井技术领域，设及一种流体替换的低渗透油藏水淹信息的计算 方法。
【背景技术】
[0002] 研究区目的储层物性差，储集空间小，孔隙度分布范围为4~16%，渗透率主要分 布在0.1~2Χ10-3μπι 2之间，属低孔隙度低渗透率油藏。历经20多年的注水开发，导致储层的 流体性质、孔隙结构和油水分布等发生变化。伴随着含水率的上升，水淹状况复杂，电性变 化越来越大，水淹级别判识困难。水淹级别的准确划分是注水开发油田调整工作中必须解 决的问题，在油田老区提高采收率和稳定产量等方面起着至关重要的作用。目前，水淹层水 淹级别评价方法较多，如测井曲线定性识别法、交会图版法和可动流体分析法等，但国内各 油田由于其地质状况和开发程度的不同，目前的评价方法都具有区域局限性，没有一种通 用的快速定性的低渗透油藏水淹信息的计算方法，难W推广应用。重构水淹层原始地层电 阻率快速定性的划分水淹级别是有效方法之一，但由于运方面资料较少，且重构方法多基 于大量的岩屯、实验数据建立的，受昂贵取屯、费用的限制和岩屯、实验分析资料的匿乏，应用 范围不大，尤其是油田老区二次开发中，难W满足实际生产需求。
【发明内容】
[0003] 本发明的目的在于克服上述现有技术的缺点，提供了一种流体替换的低渗透油藏 水淹信息的计算方法，该方法准确得到低渗透油藏水淹信息。
[0004] 为达到上述目的，本发明所述的流体替换的低渗透油藏水淹信息的计算方法包括 W下步骤：
[0005] 1)获取区域中储层水淹前后测井的响应特征，并从所述区域中储层水淹前后测井 的响应特征中选取水淹前后岩性和物性无明显变化的自然伽马GR及声波时差AC测井曲线；
[0006] 2)选取区域中未水淹的模型井的测井信息，再根据未水淹的模型井的测井信息建 立储层孔隙中所含原始流体的电阻率与自然伽马及声波时差的函数关系；
[0007] 3)利用步骤2)得到的储层孔隙中所含原始流体的电阻率与自然伽马及声波时差 的函数关系计算对区域中已水淹的模型井，确定已水淹的模型井的储层孔隙中所含原始流 体的电阻率ti，并将已水淹的模型井的储层孔隙中所含原始流体的电阻率ti与实际测量得 到的已水淹的模型井的电阻率t2进行对比，得已水淹的模型井的储层孔隙中所含原始流体 的电阻率ti与实际测量得到的已水淹的模型井的电阻率t2之间的电阻率变化值ARt;
[000引4)由区域中已水淹的模型井的试油信息确定含水率Fw，根据所述含水率FwW及步 骤2)得到的电阻率变化值Δ Rt建立电阻率变化值Δ Rt识别水淹级别的划分标准；
[0009] 5)利用步骤2)得到的储层孔隙中所含原始流体的电阻率与自然伽马及声波时差 的函数关系对区域中待评价井进行计算，得区域中待评价井的储层孔隙中所含原始流体的 电阻率，然后将区域中待评价井的储层孔隙中所含原始流体的电阻率与实际测量得到的待 评价井的电阻率进行对比，确定区域中待评价井的储层孔隙中所含原始流体的电阻率与实 际测量得到的待评价井的电阻率之间的电阻率变化值A Rs，然后根据所述电阻率变化值Δ RsW及步骤4)得到的水淹级别的划分标准确定区域中待评价井的水淹信息。
[0010] 根据未水淹的模型井的测井信息建立储层孔隙中所含原始流体的电阻率与自然 伽马及声波时差的函数关系的具体操作为：
[0011] 根据未水淹的模型井的测井信息基于模糊聚类数学方法建立储层孔隙中所含原 始流体的电阻率与自然伽马及声波时差的函数关系。
[0012] 根据未水淹的模型井的测井信息建立储层孔隙中所含原始流体的电阻率与自然 伽马及声波时差的函数关系的具体操作为：
[0013] 根据所述模型井的测井信息确定待测区域的电阻率样本Rt，再对待测区域的电阻 率样本化进行分类，并设置待测区域的电阻率样本Rt的模糊识别矩阵U和中屯、指标向量S， 然后通过循环迭代求解最优模糊识别矩阵r、最优模糊聚类中屯、指标S'及变量权重W'，然 后根据最优模糊识别矩阵U'、最优模糊聚类中屯、指标s'及变量权重W'确定储层孔隙中所含 原始流体的电阻率与自然伽马和声波时差的函数关系，其中，
[0014] Rt.=aXH+b;H=f(GR AC U· S· r P)
[0015] 式中，Rt'为储层孔隙中所含原始流体的电阻率曲线数据，a和b为模型系数，f为映 射的函数关系，GR和AC为归一化后的自然伽马和声波时差测井数据，0.9 < P。
[0016] 采用模糊聚类数学方法确定模型系数a和b，得a = 0.1783，b = 1.003。
[0017]采用模糊聚类数学方法确定最优模糊聚类中屯、指标S '，得 矿_(0.4625 0.5150 0.5857 0.6779 ~'0.5884 0.6175 0.6293 0.7399
[001引采用模糊聚类数学方法确定变量权重W'，得W' = (0.6425 0.3575);
[0019] p = 0.9017。
[0020] 本发明具有W下有益效果：
[0021] 本发明所述的流体替换的低渗透油藏水淹信息的计算方法在计算低渗透油藏水 淹信息的过程中，通过建立未水淹的模型井的测井信息建立储层孔隙中所含原始流体的电 阻率与自然伽马及声波时差的函数关系，然后计算得到待测区域中储层孔隙中所含原始流 体的电阻率，再根据实测的与计算得到的电阻率的差值得到待测区域的低渗透油藏水淹信 息，本发明能够满足实际生产的需要，经油田老区二次开发的120口实际井资料处理，获取 的低渗透油藏水淹信息符合率达到83% W上，具有较好的应用价值。
【附图说明】
[0022] 图1为本发明的实施例一中模型井电阻率样本Rt的累计频率分布图；
[0023] 图2为本发明的实施例一中含水率Fw与电阻率变化值Δ Rt的变化规律图；
[0024] 图3(a)为本发明的实施例一中60-63号层的储层电阻率变化值ARt示意图；
[0025] 图3(b)为本发明的实施例一中64-68号层的储层电阻率变化值Δ Rt示意图；
[00%]图3(C)为本发明的实施例一中56-59号层的储层电阻率变化值Δ Rt示意图。
【具体实施方式】
[0027] 下面结合附图对本发明做进一步详细描述：
[0028] 本发明所述的流体替换的低渗透油藏水淹信息的计算方法包括W下步骤：
[0029] 1)获取区域中储层水淹前后测井的响应特征，并从所述区域中储层水淹前后测井 的响应特征中选取水淹前后岩性和物性无明显变化的自然伽马GR及声波时差AC测井曲线；
[0030] 2)选取区域中未水淹的模型井的测井信息，再根据未水淹的模型井的测井信息建 立储层孔隙中所含原始流体的电阻率与自然伽马及声波时差的函数关系；
[0031] 3)利用步骤2)得到的储层孔隙中所含原始流体的电阻率与自然伽马及声波时差 的函数关系计算对区域中已水淹的模型井，确定已水淹的模型井的储层孔隙中所含原始流 体的电阻率ti，并将已水淹的模型井的储层孔隙中所含原始流体的电阻率ti与实际测量得 到的已水淹的模型井的电阻率t2进行对比，得已水淹的模型井的储层孔隙中所含原始流体 的电阻率ti与实际测量得到的已水淹的模型井的电阻率t2之间的电阻率变化值ARt;
[0032] 4)由区域中已水淹的模型井的试油信息确定含水率Fw，根据所述含水率FwW及步 骤2)得到的电阻率变化值Δ Rt建立电阻率变化值Δ Rt识别水淹级别的划分标准；
[0033] 5)利用步骤2)得到的储层孔隙中所含原始流体的电阻率与自然伽马及声波时差 的函数关系对区域中待评价井进行计算，得区域中待评价井的储层孔隙中所含原始流体的 电阻率，然后将区域中待评价井的储层孔隙中所含原始流体的电阻率与实际测量得到的待 评价井的电阻率进行对比，确定区域中待评价井的储层孔隙中所含原始流体的电阻率与实 际测量得到的待评价井的电阻率之间的电阻率变化值A Rs，然后根据所述电阻率变化值Δ RsW及步骤4)得到的水淹级别的划分标准确定区域中待评价井的水淹信息。
[0034] 根据未水淹的模型井的测井信息建立储层孔隙中所含原始流体的电阻率与自然 伽马及声波时差的函数关系的具体操作为:根据未水淹的模型井的测井信息基于模糊聚类 数学方法建立储层孔隙中所含原始流体的电阻率与自然伽马及声波时差的函数关系。
[0035] 根据未水淹的模型井的测井信息建立储层孔隙中所含原始流体的电阻率与自然 伽马及声波时差的函数关系的具体操作为：
[0036] 根据所述模型井的测井信息确定待测区域的电阻率样本Rt，再对待测区域的电阻 率样本Rt进行分类，并设置待测区域的电阻率样本Rt的模糊识别矩阵U和中屯、指标向量S，然 后通过循环迭代求解最优模糊识别矩阵r、最优模糊聚类中屯、指标S'及变量权重W'，然后 根据最优模糊识别矩阵U'、最优模糊聚类中屯、指标S'及变量权重W'确定储层孔隙中所含原 始流体的电阻率与自然伽马和声波时差的函数关系，其中，
[0037] Rt.=aXH+b;H=f(GR AC U· S· W· P)
[0038] 式中，Rt'为储层孔隙中所含原始流体的电阻率曲线数据，a和b为模型系数，f为映 射的函数关系，GR和AC为归一化后的自然伽马和声波时差测井数据，0.9