CN109958413B - 一种特高含水期油藏动态流动单元划分方法 - Google Patents

一种特高含水期油藏动态流动单元划分方法 Download PDF

Info

Publication number
CN109958413B
CN109958413B CN201910154632.6A CN201910154632A CN109958413B CN 109958413 B CN109958413 B CN 109958413B CN 201910154632 A CN201910154632 A CN 201910154632A CN 109958413 B CN109958413 B CN 109958413B
Authority
CN
China
Prior art keywords
permeability
water phase
water
phase
flow unit
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201910154632.6A
Other languages
English (en)
Other versions
CN109958413A (zh
Inventor
崔传智
杨勇
孙业恒
杜玉山
刘志宏
黄迎松
隋迎飞
唐从见
郭长春
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China University of Petroleum East China
Exploration and Development Research Institute of Sinopec Shengli Oilfield Co
Original Assignee
China University of Petroleum East China
Exploration and Development Research Institute of Sinopec Shengli Oilfield Co
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by China University of Petroleum East China, Exploration and Development Research Institute of Sinopec Shengli Oilfield Co filed Critical China University of Petroleum East China
Priority to CN201910154632.6A priority Critical patent/CN109958413B/zh
Publication of CN109958413A publication Critical patent/CN109958413A/zh
Application granted granted Critical
Publication of CN109958413B publication Critical patent/CN109958413B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • E—FIXED CONSTRUCTIONS
    • E21—EARTH DRILLING; MINING
    • E21B—EARTH DRILLING, e.g. DEEP DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
    • E21B43/00—Methods or apparatus for obtaining oil, gas, water, soluble or meltable materials or a slurry of minerals from wells
    • E21B43/16—Enhanced recovery methods for obtaining hydrocarbons
    • E21B43/20—Displacing by water
    • E—FIXED CONSTRUCTIONS
    • E21—EARTH DRILLING; MINING
    • E21B—EARTH DRILLING, e.g. DEEP DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
    • E21B49/00—Testing the nature of borehole walls; Formation testing; Methods or apparatus for obtaining samples of soil or well fluids, specially adapted to earth drilling or wells
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation

Abstract

本发明涉及针对注水开发油藏在特高含水期实现动态流动单元划分的方法，该方法包括：步骤(1)：根据油田已有的相渗曲线，回归得到不同绝对渗透率下的相渗曲线表达式，进而计算得到各网格点的水相相对渗透；步骤(2)：计算水相渗流系数，由步骤(1)得到各网格点的水相相对渗透，根据渗流系数定义，计算得到模型中每个网格点的水相渗流系数；步骤(3)：确定划分界限，绘制水相渗流系数半对数累积曲线图，找出两个平行段(斜率较小，几乎为0)以及两个平行段间斜率较大的一段，作为动态流动单元的划分界限；步骤(4)：网格所属分区判定，在步骤(3)划分界限的基础上，对各个网格的水相渗流系数进行分区判定，进而划分流动单元。

Description

一种特高含水期油藏动态流动单元划分方法
技术领域
本发明涉及注水开发油藏综合调整技术领域，特别是涉及一种油藏进入特高含水期后，在分析地质参数的基础上，基于油藏动态参数进行流动单元划分的方法。
背景技术
我国东部油田目前大都已进入高含水、特高含水阶段，在寻找新的后备储量难度越来越大、储采比严重失调的形势下，开发老油田剩余油已成为我国东部老油田持续发展的重大战略。揭示剩余油的数量和分布是设计和优化注水开发方案和三次采油方案的重要依据，而油藏内剩余油的数量大小、空间分布受控于储层宏观和微观非均质性质。流动单元概念的提出和及其研究方法的形成与发展为认识油藏的非均质性提供了有效手段。流动单元定义为在空间上连续分布的具有相似岩石物理特征和渗流特征的储集体，流动单元的研究对分析地下油水分布规律，预测剩余油分布具有很大的实际意义，然而，作为一种地质体描述定量化手段，现有的流动单元划分方法并没有很好地将地质与油藏结合起来，无法描述与解决开发后期无效驱动的问题，为油田开发调整提供依据。本专利提出将渗流系数作为流动单元划分的依据指标，能够将地质参数和开发动态参数很好的结合起来，提出了特高含水期流动单元划分依据指标的政策界限，建立了一套油藏特高含水期流动单元划分方法，从而为开发后期制定针对性的高效调控技术提供支撑。
发明内容
本发明的目的是提供一种适用于油藏特高含水期的动态流动单元划分方法。针对特高含水期由于油藏非均质引起的不均衡驱替，以渗流系数作为依据指标，在地质参数划分的基础上，以数值模拟计算结果为基础，实现动态流动单元的划分。
本发明可以通过如下技术措施来实现：
步骤(1)：计算不同网格点的水相相对渗透率。不同绝对渗透率对应的相对渗透率曲线不同，根据不同绝对渗透率对应的相对渗透率曲线表达式和各网格点含水饱和度，计算得到各网格点的水相相对渗透率；
步骤(2)：计算水相渗流系数。由步骤(1)得到各网格点的水相相对渗透率，根据渗流系数定义，计算得到地层中每个网格点的水相渗流系数；
步骤(3)：确定划分界限。绘制水相渗流系数半对数累积曲线图，找出两个平行段(斜率较小，几乎为0)以及两个平行段间斜率较大的一段，作为动态流动单元的划分界限；
步骤(4)：网格所属分区判定。在步骤(3)划分界限的基础上，对各个网格的水相渗流系数进行分区判定，进而划分流动单元。
关键技术要点包括：
1.水相相对渗透率的计算
不同绝对渗透率储层对应的相渗曲线不同。水相相对渗透率(krw)为关于渗透率(k)和含水饱和度(Sw)的函数，采用的表达式为：
Figure BDA0001982522570000021
其中：
krw为水相相对渗透率，无因次；
krwi为水相端点渗透率，即残余油饱和度下的水的相对渗透率；
Sw为油藏含水饱和度；
Swc为束缚水饱和度；
Soi为原始含油饱和度；
Sor为残余油饱和度；
n为水相幂指数；
通过回归方法，建立相渗曲线特征参数值与绝对渗透率的关系，特征参数包括：原始含油饱和度(Soi)、残余油饱和度(Sor)、水相幂指数(n)以及残余油下水相渗透率(krwi)。方法如下：
首先，对上式变形，方便下面计算，变形如下：
Figure BDA0001982522570000031
然后，对不同渗透率的实际相对渗透率曲线(已知)进行处理，利用上式，分别得到不同实际相对渗透率曲线中，各特征参数值与渗透率的对应关系，并分别绘制各特征参数值与渗透率的散点图；
最后，根据各特征参数与渗透率散点图，分别回归出各特征参数与渗透率的关系式，并给出渗透率的范围(即实际相对渗透率曲线的渗透率范围)。
根据各特征参数与渗透率的关系式和水相相对渗透率表达式，由每个网格的绝对渗透率和含水饱和度，求得对应的水相相对渗透率。
2.水相渗流系数的计算
流动单元定义为具有相似渗流能力的储集体，因此提出渗流系数为流动能力的表征指标。渗流系数的表达式为：
Figure BDA0001982522570000032
其中：
k为储层绝对渗透率，μm2；
μw为水相粘度，mPa·s。
3.划分界限的确定
计算得到每个网格点的水相渗流系数(Zw)，并绘制其半对数累积分布图，将出现的第一个平行段(斜率较小，几乎为0)作为第一类流动单元的划分标准，此平行段之前(包括平行段)为第一类流动单元；将出现的第二个平行段(斜率较小，几乎为0)作为第三类流动单元的划分标准，此平行段之后(包括平行段)为第三类流动单元；将两个平行段之间，斜率较大的上升段作为第二类流动单元的划分标准，划分第二类流动单元。
4.流动单元的划分
依据步骤3的划分标准，以及步骤2计算的水相渗流系数，对每个网格的水相渗流系数对数值进行判别，进而划分为相应的流动单元类型；对所有网格判别完成后，最终在平面图中显示出来，完成流动单元的划分(如图3所示)。
附图说明
图1为实例中所用油藏概念模型示意图；
图2为本发明的特高含水期动态流动单元划分方法——水相渗流系数半对数累积分布图；
图3为本发明的特高含水期流动单元划分结果平面图。
具体实施方式
为了进一步描述本发明，采用实例作详细说明。
考虑油藏的非均质性，根据所需条件，建立五点法数值模拟模型，并对其特高含水期进行流动单元划分，模型参数设置如下：
模型注采井距均为300m；模型有一条高渗带，渗透率为1000mD，高渗带两侧渗透率为500mD；边界处渗透率较低，最低为50mD，且向两侧渗透率渐变；注入井(P1、P2、P3、P4)分布在模型四周，生产井(I1)在模型中心。模型示意图如图1所示。
步骤1，以胜利油田100条相对渗透率曲线回归整理，建立的相渗曲线特征参数与绝对渗透率的关系式如下：
原始含油饱和度：Soi＝0.0327ln(k)+0.6764；
残余油饱和度：Sor＝-0.0131ln(k)+0.1946；
水相幂指数：n＝0.0056k+1.9146；
残余油下水相相对渗透率：krwi＝-0.038ln(k)+0.2594；
渗透率k的变化范围：0.06～6.3μm2，可计算该范围内任意渗透率对应的相渗曲线。
根据相对渗透率表达式，以及上述各式，得到所需的渗透率为500mD和1000mD时的相对渗透率表达式。
渗透率为500mD时，水相相对渗透率表达式为
Figure BDA0001982522570000051
渗透率为1000mD时，水相相对渗透率表达式为
Figure BDA0001982522570000052
当油藏模型含水率达到90％时，得到每个网格的含水饱和度Sw，根据500mD和1000mD的水相对渗透率表达式，分别得到每个网格的水相相对渗透率(krw)。
步骤2，根据公式
Figure BDA0001982522570000053
以及每个网格水相相对渗透率(krw)值，得到每个网格的水相渗流系数。
步骤3，根据水相渗流系数(Zw)的数据，作出其半对数累积分布图(如图2所示)，根据划分标准，找出两个平行段(斜率较小，几乎为0)以及两个平行段间斜率较大的一段，将斜率突变处水相渗流系数值作为划分界限，如图1中虚线所示(突变处Zw值分别为1.803×10-4和206.3)，并依据此划分为三类流动单元。
步骤4，依据步骤3的界限，根据每个网格的水相渗流系数进行判别，从而完成流动单元的划分(如图3所示)。

Claims (3)

1.特高含水期油藏动态流动单元划分方法，包括以下步骤：
步骤(1)：计算水相相对渗透率；根据不同绝对渗透率对应的相渗曲线表达式和各网格点含水饱和度，计算得到各网格点的水相相对渗透率；
步骤(2)：计算水相渗流系数；由步骤(1)得到各网格点的水相相对渗透率，根据渗流系数计算公式
Figure FDA0002922170010000011
计算得到地层中每个网格点的水相渗流系数，其中，Zw为渗流系数，krw为水相相对渗透率，无因次；Sw为油藏含水饱和度；k为储层绝对渗透率，μm2；μw为水相粘度，mPa·s；
步骤(3)：确定划分界限；绘制水相渗流系数半对数累积曲线图，找出两个平行段以及两个平行段间斜率较大的一段，将斜率突变处水相渗流系数值的作为动态流动单元的划分界限；
步骤(4)：网格所属分区判定；在步骤(3)划分界限的基础上，对各个网格的水相渗流系数进行分区判定，进而划分流动单元。
2.根据权利要求1所述的特高含水期油藏动态流动单元划分方法，其特征在于，在步骤(1)中，根据相渗曲线的表达公式，通过对多条相渗曲线回归，得到油藏不同渗透率下的相渗曲线表达式，进而得到每个网格点的水相相对渗透率。
3.根据权利要求1-2中任一项所述的特高含水期油藏动态流动单元划分方法，其中将水相渗流系数Zw作为流动单元划分的政策界限，在Zw半对数累积曲线图中，找出两个平行段以及两个平行段间斜率较大的一段，将斜率突变处水相渗流系数值作为划分界限，最终完成特高含水期流动单元的划分。
CN201910154632.6A 2019-03-01 2019-03-01 一种特高含水期油藏动态流动单元划分方法 Active CN109958413B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910154632.6A CN109958413B (zh) 2019-03-01 2019-03-01 一种特高含水期油藏动态流动单元划分方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910154632.6A CN109958413B (zh) 2019-03-01 2019-03-01 一种特高含水期油藏动态流动单元划分方法

Publications (2)

Publication Number Publication Date
CN109958413A CN109958413A (zh) 2019-07-02
CN109958413B true CN109958413B (zh) 2021-03-30

Family

ID=67023916

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910154632.6A Active CN109958413B (zh) 2019-03-01 2019-03-01 一种特高含水期油藏动态流动单元划分方法

Country Status (1)

Country Link
CN (1) CN109958413B (zh)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112282727B (zh) * 2019-07-11 2023-12-26 中国石油天然气股份有限公司 油藏渗流场分区的方法和装置
CN110424944B (zh) * 2019-07-22 2022-08-23 中国海洋石油集团有限公司 巨厚碳酸盐岩油藏油水拟相渗曲线的建立方法
CN110702579B (zh) * 2019-10-21 2022-05-20 长江大学 基于渗流实验的深层非均质滩相气藏流动单元划分方法
CN111535787B (zh) * 2020-04-09 2022-02-22 中国石油化工股份有限公司 高含水油藏动态渗流界面的判识模型及判识界限构建方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CA2749284A1 (en) * 2009-01-19 2010-07-22 Schlumberger Canada Limited Estimating petrophysical parameters and invasion profile using joint induction and pressure data inversion approach
CN105298477A (zh) * 2014-07-07 2016-02-03 中国石油化工股份有限公司 一种基于流动单元的地层孔隙结构解释方法
CN107358374A (zh) * 2017-08-07 2017-11-17 中国石油集团渤海钻探工程有限公司 一种基于测试资料和流动单元的压裂决策方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CA2749284A1 (en) * 2009-01-19 2010-07-22 Schlumberger Canada Limited Estimating petrophysical parameters and invasion profile using joint induction and pressure data inversion approach
CN105298477A (zh) * 2014-07-07 2016-02-03 中国石油化工股份有限公司 一种基于流动单元的地层孔隙结构解释方法
CN107358374A (zh) * 2017-08-07 2017-11-17 中国石油集团渤海钻探工程有限公司 一种基于测试资料和流动单元的压裂决策方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
Developing a workflow to represent fractured carbonate reservoirs for simulation models under uncertainties based on flow unit concept;Mahjour, SK 等;《OIL & GAS SCIENCE AND TECHNOLOGY-REVUE D IFP ENERGIES NOUVELLES》;20190227;第74卷(第15期);第2-3页 *
基于辫状河储层构型的流动单元划分及其分布规律;王石 等;《油气地质与采收率》;20150930;第22卷(第5期);第47-51页 *

Also Published As

Publication number Publication date
CN109958413A (zh) 2019-07-02

Similar Documents

Publication Publication Date Title
CN109958413B (zh) 一种特高含水期油藏动态流动单元划分方法
CN106651610B (zh) 一种浅层超低渗砂岩油藏注水开发动态分析方法
CN110334431A (zh) 一种低渗透致密气藏单井控制储量计算及剩余气分析方法
CN103745022B (zh) 聚合物驱后变流线调整方法
CN103939066B (zh) 一种一注多采井组定注水量确定油井产液量的方法
CN110617042B (zh) 高耗水层带发育油藏分层注水调控方法
CN105756654A (zh) 水驱油藏井网优化设计方法
CN105386751B (zh) 一种基于油藏渗流模型的水平井测井产能预测方法
CN108301813B (zh) 零散砂体发育的多层砂岩油藏开发调整方法
CN103244087B (zh) 一种低渗透油藏调剖堵水选井决策方法
CN102777157A (zh) 一种co2驱油气水异井注入油藏混驱开发方法
CN105822271A (zh) 厚层油藏变流线矢量调整方法
CN105626010A (zh) 一种分段注水井中注水层段合理划分方法
CN106150463B (zh) 一种砾岩油藏聚合物驱注入压力升幅确定方法
CN109359332A (zh) 一种浅薄层油藏数值模拟模型建立方法及转蒸汽驱的方法
CN109858177A (zh) 一种基于拟流线法的底水油藏水平井水驱数值模拟方法
CN104033137A (zh) 利用油田污水提高断块油藏采收率方法
CN116108572A (zh) 页岩凝析气井体积压裂外区产能贡献分析方法
CN108561129A (zh) 小断块油藏油水过渡带剩余油快速评价方法
CN111350485B (zh) 井网调整方法及装置
CN114676652A (zh) 一种欠饱和储层煤层气井生产边界定量预测的计算方法
CN105275457B (zh) 整装老油田河流相储层水平井部署最佳模式研究方法
CN114818229A (zh) 一种二氧化碳驱开发模式下提高采收率评价方法
CN113468476A (zh) 一种聚合物驱增油量预测方法及装置
CN113653474B (zh) 一种实现高含水老油田绿色开发的方法

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant