< P。
[0039] 采用模糊聚类数学方法确定模型系数a和b，得a = 0.1783，b = 1.003。
[0040] 采用模糊聚类数学方法确定最优模糊聚类中屯、指标S '，得 及·（0.4625 0.5150 0.5857 0.6779 >。 0.5884 0.6175 0.6293 0.7399 U
[0041] 采用模糊聚类数学方法确定变量权重W'，得W' = (0.6425 0.3575);
[0042] p = 0.9017。
[00创实施例一
[0044]根据区域油藏开发背景，通过对研究区储层水淹前后测井响应特征的对比分析， 优选了水淹前后岩性和物性无明显变化的自然伽马和声波时差测井曲线；按照《水淹层测 井资料处理与解释规范(SY/T6178-2011)》中含水率的划分界限，结合试油资料，选取有代 表性未水淹的模型井(Fw含10%为未水淹)的测井信息，采用模糊聚类数学方法建立储层孔 隙中所含原始流体的电阻率与自然伽马和声波时差的函数关系，图1为模型井电阻率Rt样 本集的累计频率分布图，通过电阻率样本集可划分为4类，其函数关系为:Rt'=aXH+b;H = f(GR AC U' S' W' P)，式中，Rt'为储层孔隙中所含原始流体的电阻率曲线数据;a和b为模 型的系数，a = 0.1783，b= 1.003;f为映射的函数关系;GR和AC为归一化后的自然伽马和声 波时差测井数据;S '为最优模糊聚类中屯、矩阵， , 1 。， 0/4625 0.5150 0"巧57 0.6779 .、，女 ^ 玉革 0045 矿二（ ）;W为变量权重， 0.5884 0.6175 0.6293 0,7399
[0046] W· = (0.6425 0.3575);P为相关系数，P = 0.9017，利用上述建立的函数关系，对所 述研究区选取有代表性已水淹的模型井进行计算，计算出储层孔隙中所含原始流体的电阻 率，并与实测电阻率对比，实测电阻率与计算电阻率之差为电阻率变化值A Rt，电阻率变化 值A Rt由储层孔隙中流体性质变化引起的;对电阻率变化值Δ Rt与含水率Fw进行分析，两者 呈现正相关，且满足一元二次函数关系，然后结合含水率Fw的划分界限，利用电阻率变化值 Δ Rt得低渗透油藏水淹信息。在实际应用时，水淹级别的划分标准为A Rt。. 83 Ω .m为弱水 淹，3.83Ω .m<A化巧1.34Ω .m为中水淹，ARt>5l.34Q .m为强水淹，图2为电阻率变化值 ARt与含水率Fw的关系图。
[0047] 实际资料处理过程中，反演储层孔隙中所含原始流体的电阻率曲线通过编写程序 实现。图3为基于流体替换的水淹级别测井解释成果图：参考图3(a)，60号层基于流体替换 法反演的储层孔隙中所含原始流体的电阻率接近于实测电阻率，电阻率变化值A Rt为1.46 Ω -m，小于3.83Ω -m，依据水淹级别的划分标准，测井解释为弱水淹。经生产动态数据验 证，该井在1435.0~1441. Om井段射孔，试油：油7.5t，水4.05m3，含水率35.06 % ;投产：油 5.84t，水2.29m3,含水率14.5%，试油结论为弱水淹，测井解释与生产实际相符。参考图3 (b)，65号层基于流体替换法计算的储层孔隙中所含原始流体电阻率比实测电阻率小，电阻 率变化值A Rt为16.44 Ω . m，介于3.83和51.34 Ω . m，依据水淹级别的划分标准，测井解释 为中水淹。经生产动态数据验证，该井在1325~1328m井段射孔，试油：油10.08t，水2.74m 3， 含水率21.3 %，投产：油1.59t，水2.56m3，含水率61.6 %，试油结论为中水淹，测井解释与生 产实际相符。参考图3(c)，58号层基于流体替换法计算的储层孔隙中所含原始流体电阻率 远小于实测电阻率，电阻率变化值A Rt为116.74 Ω . m，大于51.34 Ω . m，依据水淹级别的 划分标准，测井解释为为强水淹。经生产动态数据验证，该井在1642~1644m井段射孔，试 油：油花，水12.6m 3，投产:油0.79t，水4.95m3，含水率86.2 %，试油结论为强水淹，测井解释 符合生产实际。经120口实际井资料的处理，并与试油测试结果对比，解释符合率达到83%， 验证了基于流体替换的低渗透油藏水淹信息计算方法的有效性，具有较好的应用效果。
[0〇4引 W上所述仅为本发明的较佳实施图，并不用W限制本发明，对于本领域的技术人 员，凡在本发明的精神和原则之内，所作的任何修改与改进等，均应包含本发明的保护范围 内。
【主权项】
1. 一种流体替换的低渗透油藏水淹信息的计算方法，其特征在于，包括以下步骤： 1) 获取区域中储层水淹前后测井的响应特征，并从所述区域中储层水淹前后测井的响 应特征中选取水淹前后岩性和物性无明显变化的自然伽马GR及声波时差AC测井曲线； 2) 选取区域中未水淹的模型井的测井信息，再根据未水淹的模型井的测井信息建立储 层孔隙中所含原始流体的电阻率与自然伽马及声波时差的函数关系； 3) 利用步骤2)得到的储层孔隙中所含原始流体的电阻率与自然伽马及声波时差的函 数关系计算对区域中已水淹的模型井，确定已水淹的模型井的储层孔隙中所含原始流体的 电阻率 tl，并将已水淹的模型井的储层孔隙中所含原始流体的电阻率t与实际测量得到的 已水淹的模型井的电阻率^进行对比，得已水淹的模型井的储层孔隙中所含原始流体的电 阻率 tl与实际测量得到的已水淹的模型井的电阻率^之间的电阻率变化值ARt; 4) 由区域中已水淹的模型井的试油信息确定含水率Fw，根据所述含水率Fw以及步骤2) 得到的电阻率变化值△ Rt建立电阻率变化值A Rt识别水淹级别的划分标准； 5) 利用步骤2)得到的储层孔隙中所含原始流体的电阻率与自然伽马及声波时差的函 数关系对区域中待评价井进行计算，得区域中待评价井的储层孔隙中所含原始流体的电阻 率，然后将区域中待评价井的储层孔隙中所含原始流体的电阻率与实际测量得到的待评价 井的电阻率进行对比，确定区域中待评价井的储层孔隙中所含原始流体的电阻率与实际测 量得到的待评价井的电阻率之间的电阻率变化值A Rs，然后根据所述电阻率变化值△匕以 及步骤4)得到的水淹级别的划分标准确定区域中待评价井的水淹信息。2. 根据权利要求1所述的流体替换的低渗透油藏水淹信息的计算方法，其特征在于，根 据未水淹的模型井的测井信息建立储层孔隙中所含原始流体的电阻率与自然伽马及声波 时差的函数关系的具体操作为： 根据未水淹的模型井的测井信息基于模糊聚类数学方法建立储层孔隙中所含原始流 体的电阻率与自然伽马及声波时差的函数关系。3. 根据权利要求2所述的流体替换的低渗透油藏水淹信息的计算方法，其特征在于，根 据未水淹的模型井的测井信息建立储层孔隙中所含原始流体的电阻率与自然伽马及声波 时差的函数关系的具体操作为： 根据所述模型井的测井信息确定待测区域的电阻率样本Rt，再对待测区域的电阻率样 本Rt进行分类，并设置待测区域的电阻率样本Rt的模糊识别矩阵U和中心指标向量S，然后通 过循环迭代求解最优模糊识别矩阵IT、最优模糊聚类中心指标S'及变量权重W'，然后根据 最优模糊识别矩阵IT、最优模糊聚类中心指标S'及变量权重Γ确定储层孔隙中所含原始流 体的电阻率与自然伽马和声波时差的函数关系，其中， Rt.=aXH+b;H=f(GR AC U· S· W· P) 式中，Rt'为储层孔隙中所含原始流体的电阻率曲线数据，a和b为模型系数，f为映射的 函数关系，GR和AC为归一化后的自然伽马和声波时差测井数据，0.9 < P。4. 根据权利要求3所述的流体替换的低渗透油藏水淹信息的计算方法，其特征在于，采 用模糊聚类数学方法确定模型系数a和b，得a = 0.1783，b = 1.003。5. 根据权利要求3所述的流体替换的低渗透油藏水淹信息的计算方法，其特征在于，采 用模糊聚类数学方法确定最优模糊聚类中心指标S '，得6. 根据权利要求3所述的流体替换的低渗透油藏水淹信息的计算方法，其特征在于，采 用模糊聚类数学方法确定变量权重W·，得W· = (0.6425 0.3575)。7. 根据权利要求3所述的流体替换的低渗透油藏水淹信息的计算方法，其特征在于，P = 0.9017。
【专利摘要】本发明公开了一种流体替换的低渗透油藏水淹信息的计算方法，包括：1)获取区域中储层水淹前后测井的响应特征，并从其中选取水淹前后岩性和物性无明显变化的自然伽马GR及声波时差AC测井曲线；2)选取未水淹的模型井的测井信息，建立储层孔隙中所含原始流体的电阻率与自然伽马及声波时差的函数关系；3)所述函数关系对区域中已水淹的模型井进行计算，得已水淹的模型井中储层孔隙中所含原始流体的电阻率t1，并与实际测量得到的电阻率t2对比，得电阻率变化值；4)建立利用电阻率变化值ΔRt识别水淹级别的划分标准；5)根据步骤4)建立的利用电阻率变化值ΔRt识别水淹级别的划分标准准确有效的确定区域中待评价井的水淹信息。
【IPC分类】G06Q50/02
【公开号】CN105447762
【申请号】CN201510901137
【发明人】成志刚, 石玉江, 罗少成, 陈玉林, 李素娟, 张海涛, 郑小敏, 肖飞, 吴有彬, 唐冰娥
【申请人】中国石油天然气集团公司, 中国石油集团测井有限公司
【公开日】2016年3月30日
【申请日】2015年12月8日