CN115874997A - 一种特高含水后期低效无效循环条带定量判识方法 - Google Patents
一种特高含水后期低效无效循环条带定量判识方法 Download PDFInfo
- Publication number
- CN115874997A CN115874997A CN202111156399.9A CN202111156399A CN115874997A CN 115874997 A CN115874997 A CN 115874997A CN 202111156399 A CN202111156399 A CN 202111156399A CN 115874997 A CN115874997 A CN 115874997A
- Authority
- CN
- China
- Prior art keywords
- well
- water injection
- water
- seepage
- evaluation
- 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.)
- Pending
Links
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 title claims abstract description 254
- 238000000034 method Methods 0.000 title claims abstract description 54
- 238000002347 injection Methods 0.000 claims abstract description 168
- 239000007924 injection Substances 0.000 claims abstract description 168
- 238000004519 manufacturing process Methods 0.000 claims abstract description 37
- 238000011161 development Methods 0.000 claims abstract description 11
- 238000011156 evaluation Methods 0.000 claims description 60
- 239000011159 matrix material Substances 0.000 claims description 34
- 239000003129 oil well Substances 0.000 claims description 32
- 238000010521 absorption reaction Methods 0.000 claims description 28
- 230000035699 permeability Effects 0.000 claims description 14
- 230000003068 static effect Effects 0.000 claims description 10
- 239000004576 sand Substances 0.000 claims description 9
- 238000004891 communication Methods 0.000 claims description 6
- 238000004422 calculation algorithm Methods 0.000 claims description 5
- 238000005406 washing Methods 0.000 claims description 5
- 208000035126 Facies Diseases 0.000 claims description 4
- 238000012067 mathematical method Methods 0.000 claims description 4
- 238000004364 calculation method Methods 0.000 claims description 3
- 238000012216 screening Methods 0.000 claims description 3
- 238000007689 inspection Methods 0.000 claims description 2
- 238000010606 normalization Methods 0.000 claims description 2
- LFQSCWFLJHTTHZ-UHFFFAOYSA-N Ethanol Chemical compound CCO LFQSCWFLJHTTHZ-UHFFFAOYSA-N 0.000 claims 1
- 230000001186 cumulative effect Effects 0.000 claims 1
- 238000011084 recovery Methods 0.000 abstract description 2
- 239000010410 layer Substances 0.000 description 77
- 230000018109 developmental process Effects 0.000 description 9
- 238000012360 testing method Methods 0.000 description 7
- 239000002356 single layer Substances 0.000 description 5
- 230000015572 biosynthetic process Effects 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 238000012407 engineering method Methods 0.000 description 3
- 239000012530 fluid Substances 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 238000012512 characterization method Methods 0.000 description 2
- 238000006073 displacement reaction Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012544 monitoring process Methods 0.000 description 2
- 238000004445 quantitative analysis Methods 0.000 description 2
- 239000000243 solution Substances 0.000 description 2
- 239000000126 substance Substances 0.000 description 2
- 238000006467 substitution reaction Methods 0.000 description 2
- 238000009412 basement excavation Methods 0.000 description 1
- 238000012790 confirmation Methods 0.000 description 1
- 230000005611 electricity Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000002452 interceptive effect Effects 0.000 description 1
- 238000011835 investigation Methods 0.000 description 1
- 238000005342 ion exchange Methods 0.000 description 1
- 230000035515 penetration Effects 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
- 239000011148 porous material Substances 0.000 description 1
- 230000033764 rhythmic process Effects 0.000 description 1
- 239000011435 rock Substances 0.000 description 1
- 230000011218 segmentation Effects 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
- 230000010148 water-pollination Effects 0.000 description 1
Images
Classifications
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A10/00—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE at coastal zones; at river basins
- Y02A10/40—Controlling or monitoring, e.g. of flood or hurricane; Forecasting, e.g. risk assessment or mapping
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明属于油田勘探开发技术领域，具体涉及一种特高含水后期低效无效循环条带定量判识方法。所述方法包括以下步骤：定量识别强注水井点，确定低效无效循环起源；识别强注水井点的优势渗流层段，确定低效无效循环位置；明确优势渗流层段平面注采流动方向，确定优势渗流方向判别参数及其界限，识别低效无效循环条带。本发明方法大幅提高了低效无效循环条带识别的准确性，易于推广应用，对提高油藏采收率具有重要的指导意义。
Description
技术领域
本发明属于油田勘探开发技术领域，具体涉及一种特高含水后期低效无效循环条带定量判识方法。
背景技术
注水开发是胜利油田开发的主体，进入特高含水期的注水开发油田，仍有大量的地质储量滞留于地下，准确识别低效无效循环条带是指导油藏有效注采调整的重要依据。
CN106593409B公开了一种油田注采关系无效循环的识别方法，包括以下步骤：无效循环静态识别的步骤，测井曲线识别的步骤，无效循环动态识别的步骤。所述无效循环静态识别的步骤，进一步包括：注采井组选择的子步骤；沉积相类型选择的子步骤。所述测井曲线识别的步骤，进一步包括：厚薄油层划分的子步骤，水淹解释结果选择的子步骤，测井曲线识别的子步骤。该发明方法综合运用精细地质研究成果、剖面测试、测井解释等多项资料，明确了无效循环场的分布状况，并开展现场试验，取得了较好效果，为后续水驱开发提供了可靠依据。
付艳波、刘仁强在《运用“四参数法”判别低效无效注水循环体》中结合高精度注入，产出剖面测试资料解释结果进行分析，判断低效无效循环层；同时利用静态资料，对可疑的低效无效循环层位加以确认，利用油水井动态资料判断循环体形成的时间(付艳波,刘仁强.运用"四参数法"判别低效无效注水循环体[J].大庆石油地质与开发,2009,28(3):74-78.)。
姜红芹根据选出的低效循环井层，总结出了油田高含水后期无效和低效循环的几种类型；以"动态普查,静态确认，测试验证"为指导思想,根据高精度注，产剖面等测试资料识别层位及方向，以能大量,方便录取，信息来源充足为原则，选取油水井日常生产过程中录取的动,静态资料，通过分析筛选,选取具有独立性，能分别代表低效循环条带的不同特征及与形成低效循环条带有直接的相关性的参数，建立四参数决策系数法判定低效循环条带，圈定低效无效注水循环的可疑井组，初步应用动态分析手段判断低效无效循环形成的大体时间，最后利用静态资料加以确认(姜红芹.特高含水期低效无效循环识别研究[D].东北石油大学,2010.)。
目前低效无效循环判识方法，有的采用定性或半定量方式，参数权重人为赋值，主观性太强；有的只制定了简单的低效无效循环井层识别标准，但未对小层低效无效循环条带进行识别；有的方法操作难度大，需进行模型拟合，其拟合吻合率决定识别准确度。
特高含水后期低效无效循环条带判识方法的准确性仍需进一步提高。
发明内容
本发明主要目的在于提供特高含水后期低效无效循环条带定量判识方法，所述方法将模糊数学、测试分析、油藏工程、渗流力学等多技术手段相结合，结合静态、动态、监测、数模结果等资料，层层递进，形成一套系统识别小层各注采方向低效无效循环条带的定量识别方法。本发明方法能更加准确地判断出低效无效循环条带，为剩余油挖潜提供指导依据。
为实现上述目的，本发明采用以下技术方案：
本发明提供一种特高含水后期低效无效循环条带定量判识方法，其包括以下步骤：
定量识别强注水井点，确定低效无效循环起源；识别强注水井点的优势渗流层段，确定低效无效循环位置；明确优势渗流层段平面注采流动方向，确定优势渗流方向判别参数及其界限，识别低效无效循环条带。
进一步地，应用模糊数学方法定量识别强注水井点，具体包括以下步骤：
筛选评价注水状况的表征参数，作为评价参数；设定评价参数的评价等级，选取模糊隶属度函数计算评价参数对各评价等级的隶属度，得到单注水井点模糊评价矩阵；应用模糊层次分析法构建评价参数判断矩阵，并一致性检验；应用特征根法求取判断矩阵特征向量，归一化得到评价因素的权重向量；权重向量与模糊评价矩阵进行模糊运算，并根据最大隶属度原则，确定各水井注水能力强弱，识别出强注水井点。
更进一步地，通过分析油藏开发各动静态参数之间的相互制约关系，筛选能反映水井注水状况的渗透率、有效厚度、视吸水指数和单位厚度累注入量作为评价参数。
更进一步地，评价等级设定为注水能力强和注水能力不强两级；根据各评价参数与水井注水强弱的影响关系，选取半梯形隶属度函数，计算评价参数对各评价等级的隶属度，得到单注水井点模糊评价矩阵。
进一步地，识别强注水井点的优势渗流层段，确定低效无效循环位置的方法包括以下步骤：依据岩心水洗特征、测井曲线特征、水淹特征及沉积相资料分析解释定性识别优势渗流井层；计算强注水井点各注水层位的吸水强度和吸水强度突进系数这两个判别参数，定量识别出优势渗流井层；低效无效循环位于识别出的优势渗流井层。
进一步地，优势渗流井层平面注采流动方向的确定同时满足以下判别条件：
1)在同一砂体；
2)在同一时间生产；
3)存在连通路径；
4)流动路线上不存在其它井正在生产，为一线连通；
5)油水流动方向与其它油水流向一致，压力传导需合理。
进一步地，将瞬时注水量和注水强度作为优势渗流方向判别参数，根据这两个参数分布曲线拐点值，确定判别参数界限，同时大于或等于这两个参数判别界限的为优势渗流方向，将各优势渗流方向的含水率与区块低效无效循环含水率界限比较，识别出低效无效循环条带。
更进一步地，根据分流量方程计算各优势渗流方向的含水率
更进一步地，针对优势渗流井层，以注水井为中心，根据动静态参数，采用渗流阻力系数劈分注水量算法计算出各受效油井方向的的渗流阻力系数，确定劈分注水比例系数，计算出各受效油井方向的瞬时注水量和注水强度。
更进一步地，各受效油井方向的瞬时注水量和注水强度的计算具体包括以下步骤：
(1)计算第i口水井与第j口受效油井小层渗流阻力系数
ko＝k·kro(Sw)
kw＝k·krw(Sw)
(2)计算第i口注水井小层平面总渗流阻力
(3)计算第i口注水井第j口受效油井方向的小层注水比例系数
(4)根据第i口注水井小层注入量，计算第j口受效油井方向的瞬时注水量
Wij＝λ平面ij×Wi
(5)根据瞬时注水量，得到第j口受效油井方向的注水强度
Iij＝Wij/Hj
式中：Rij为ij两点间小层渗流阻力系数；Lij为ij间井距，m；Hj为第j口油井的小层有效厚度，m；μw为水粘度，μo为油粘度，mPa·s；koj为第j口油井方向小层平均油相渗透率，kwj为第j口油井方向小层平均水相渗透率，mD；Wi为第i口水井小层日注水量。
本发明还提供一种特高含水后期低效无效循环条带定量判识系统，其具有用于执行以上所述识别方法的功能模块。
与现有技术相比，本发明具有以下优势：
本发明方法首先应用模糊数学方法对油藏强注水能力井形成的影响参数开展综合评价，根据评价结果识别出强注水井点，明确低效无效循环起源；进一步，根据识别出的强注水井点的岩心、测井、监测等资料，应用定性解释和油藏工程方法，确定低效无效循环存在层位；再建立考虑动静态参数的渗流阻力系数劈分注水量算法，算出优势渗流井层各注采流动方向的瞬时注水量和注水强度，根据瞬时注水量、注水强度及单元含水率，识别出低效无效循环条带。
本发明方法大幅提高了低效无效循环条带识别的准确性，对提高油藏采收率具有重要的指导意义；该方法步骤简单，易于推广应用。
附图说明
构成本发明的一部分的说明书附图用来提供对本发明的进一步理解，本发明的示意性实施例及其说明用于解释本发明，并不构成对本发明的不当限定。
图1为本发明特高含水后期低效无效循环条带定量判识方法的流程图；
图2为本发明的一具体实施例所述强注水井点定量识别流程图。
具体实施方式
应该指出，以下详细说明都是示例性的，旨在对本发明提供进一步的说明。除非另有指明，本文使用的所有技术和科学术语具有与本发明所属技术领域的普通技术人员通常理解的相同含义。
需要注意的是，这里所使用的术语仅是为了描述具体实施方式，而非意图限制根据本发明的示例性实施方式。如在这里所使用的，除非上下文另外明确指出，否则单数形式也意图包括复数形式，此外，还应当理解的是，当在本说明书中使用术语“包含”和/或“包括”时，其指明存在特征、步骤、操作和/或它们的组合。
为了使得本领域技术人员能够更加清楚地了解本发明的技术方案，以下将结合具体的实施例详细说明本发明的技术方案。
实施例1
如图1所示，所述特高含水后期低效无效循环条带定量判识方法，包括以下步骤：
在步骤101，低效无效循环起源模糊判别主要应用模糊数学方法，主要包括模糊层次分析法、模糊隶属度函数及模糊最大隶属度原则。
如图2所示，首先全面分析各个参数之间的相互制约关系，优选出能科学、客观评价水井注水强弱的关键表征参数渗透率、有效厚度、视吸水指数和单位厚度累注入量作为评价因素，设定注水能力强(Ⅰ级)和注水能力不强(Ⅱ级)两个评价等级作为评价集。再选用隶属度函数确定各注水井点注水强弱评价矩阵，根据各评价参数与水井注水强弱的影响关系，选取半梯形隶属度函数，实现各注水井点各项参数的标准化，得到单注水井点模糊评价矩阵。
另外，根据Saaty理论确定各评价因素的标度，进而得到各评价因素的判断矩阵，为了检验判断矩阵是否具有满意的一致性，需比较平均随机一致性指标R.I.和一致性指标C.I.，当C.R.＝C.I./R.I.<0.10时，判断矩阵有可接受的一致性；否则，需重新修正判断矩阵。
应用特征根法求取判断矩阵的值，首先求取判断矩阵的最大特征值，根据矩阵特征值与特征向量的关系公式，求取判断矩阵最大特征值对应的特征向量，再将特征向量归一化，即可得到评价因素的权重值。最后，权重向量和评价矩阵相乘并根据最大隶属度原则，判别出区块强注水井点，即低效无效循环起源。流程进入到步骤102。
在步骤102，应用定性解释和油藏工程方法进行低效无效循环位置综合判定。
定性方面主要是依据岩心水洗特征、测井曲线特征、水淹特征及沉积相资料分析解释判定优势渗流井层；油藏工程方法是应用近几年正常吸水剖面或分层测试等监测资料，计算步骤101中判别出的强注水井点各注水层位的吸水强度和吸水强度突进系数这两个判别参数，吸水强度突进系数为单层吸水强度与全井平均单层吸水强度的比值，与判别参数界限对比，判定出优势渗流井层。低效无效循环位于判定出的优势渗流井层。进入到步骤103。
在步骤103，低效无效循环条带定量识别。
应用达西公式和水电相似性原理等知识创建渗流阻力系数劈分注水量算法，计算步骤102判定出的优势渗流井层各注采流动方向的瞬时注水量和注水强度，并结合单元含水率界限，定量识别出低效无效循环条带。
该方法首先定义出五个判别条件：注水井和采油井在某层必须在同一砂体、必须存在连通路径、必须在同一时间生产、流动路线上不存在其它生产井正在生产、油水流动方向与其它油水流向一致。同时满足这五个条件时，优势渗流井层注采井间存在流动关系。
确定出优势渗流井层存在流动关系的油水井后，以注水井为中心，计算各流动方向的渗流阻力系数，根据水电相似原理，计算对应注水井小层平面总渗流阻力，进而得到各受效油井方向的注水比例系数，根据注水比例系数，求得各流动方向的瞬时注水量和注水强度。根据单元瞬时注水量和注水强度分布规律，确定瞬时注水量和注水强度界限，将瞬时注水量和注水强度同时大于等于界限值的识别为优势渗流方向。
各注采流动方向的瞬时注水量和注水强度的计算具体包括以下步骤：
(1)计算第i口水井与第j口受效油井小层渗流阻力系数
ko＝k·kro(Sw)
kw＝k·krw(Sw)
(2)计算第i口注水井小层平面总渗流阻力
(3)计算第i口注水井第j口受效油井方向的小层注水比例系数
(4)根据第i口注水井小层注入量，计算第j口受效油井方向的瞬时注水量
Wij＝λ平面ij×Wi
(5)根据瞬时注水量，得到第j口受效油井方向的注水强度
Iij＝Wij/Hj
式中：Rij为ij两点间小层渗流阻力系数；Lij为ij间井距，m；Hj为第j口油井的小层有效厚度，m；μw为水粘度，μo为油粘度，mPa·s；koj为第j口油井方向小层平均油相渗透率，kwj为第j口油井方向小层平均水相渗透率，mD；Wi为第i口水井小层日注水量。
识别出的优势渗流方向又分为低效、无效和有效循环，根据分流量方程计算出各优势渗流方向的含水率，与单元低效、无效循环含水率界限比较，识别出低效、无效循环条带。
实施例2
所述特高含水后期低效无效循环条带定量判识方法，以胜利油区孤岛油田某单元为例，作详细说明如下：
该单元综合含水达到97.4％，已进入特高含水后期；平均单井日油1.5t/d，开发效果较差。目前该单元层间、层内、平面驱替不均衡，低效无效循环问题突出，为了控制低效无效循环注水，改善开发效果，需开展低效无效循环条带识别方法研究，比较准确地识别出低效无效循环条带，有针对性地指导注采调整，使驱替均衡。
所述方法包括以下步骤：
步骤1、低效无效循环起源判别：
如图2所示，首先全面分析各个参数之间的相互制约关系，优选出能科学、客观评价水井注水强弱的关键表征参数：渗透率、有效厚度、视吸水指数和单位厚度累注入量，并作为评价因素，设定注水能力强(Ⅰ级)和注水能力不强(Ⅱ级)两个评价等级。
隶属度函数可使原本不可比指标实际值，转化为同一规定范围内、具有可比性的数值，根据各评价参数与水井注水强弱的影响关系，选取半梯形隶属度函数计算出注水井点各评价因素对两个评价等级的隶属度，实现各注水井点各项参数的标准化，从而得到单井点的模糊评价矩阵。
半梯形隶属度函数：
rij＝[yij-min(yi)]/[max(yi)-min(yi)]
其中：rij为第i个参数第j口水井隶属度，yij为第i个参数第j口水井参数实际值，min(yi)为该单元第i个参数实际值的最小值，max(yi)为该单元第i个参数实际值的最大值。
模糊评价矩阵：
其中：rij表示第i个参数对评价集中第j个评价等级的隶属度。
对评价参数两两进行对比，根据Saaty理论，按重要性程度评价等级确定各评价参数的标度，判断矩阵元素Aij(i＝1,2，……4；j＝1,2，……4)表示参数Ai对Aj的重要性标度值，A1为视吸水指数，A2为渗透率，A3为有效厚度，A4为单位厚度累注，进而得到各评价因素的判断矩阵。
表1 Satty理论层次分析法标度表
判断矩阵：
对判断矩阵进行一致性检验(防止A1比A2重要、A2比A3重要，而A3又比A1重要的情况发生)，比较平均随机一致性指标R.I.和一致性指标C.I.，当C.R.＝C.I./R.I.<0.10时，判断矩阵有可接受的一致性，否则，需重新修正判断矩阵。其中，一致性指标为C.I.＝(λmax-n)/(n-1)，λmax为最大特征根。
表2 平均随机一致性指标R.I.数值表
| 阶数n | 1 | 2 | 3 | 4 | 5 | 6 | 7 | 8 | 9 |
| R.I. | 0.0 | 0.0 | 0.58 | 0.9 | 1.12 | 1.24 | 1.32 | 1.41 | 1.45 |
应用特征根法求取判断矩阵的值，首先求取判断矩阵的最大特征值，根据矩阵特征值与特征向量的关系公式，求取判断矩阵最大特征值对应的特征向量，再将特征向量归一化，即可得到评价因素视吸水指数、渗透率、有效厚度和单位厚度累注入量的权重向量为：
(A1,A2,A3,A4)＝(0.38,0.31,0.23,0.08)
权重向量与模糊评价矩阵相乘得到注水井点的综合评判结果：
运用最大隶属原则优选结果，模糊综合评价结果向量B＝(b1，b2)中bi最大，则判别结果为i级，最终评价结果I级的为强注水井点，即为低效无效循环起源。
表3 单元强注水井点定量识别结果
步骤2、低效无效循环位置综合判定
(1)定性解释方法
依据岩心水洗特征、测井曲线特征、水淹特征及沉积相资料分析解释判定优势渗流井层。
岩心水洗特征及其韵律性判断优势渗流层段分布于：河道砂体多段叠合正韵律分段底部。
优势渗流层段的测井曲线具备两个特征：①泥质与油质含量下降，离子交换阻力减小，自然电位幅度上升；②孔隙流体以水为主，岩性向亲水转变，电阻率幅度差下降。
纵向和平面上的水淹特征可以识别出优势渗流层段的分布位置及规律：垂向上强、中水淹层段为优势渗流层段，平面上强、中水淹主要分布在主注采线上。
按照一类河道＞二类河道＞三类河道＞废弃河道＞决口河道＞溢岸砂这个优先级顺序形成优势渗流层段。
(2)油藏工程定量判定方法
应用近几年正常吸水剖面或分层测试资料，计算各注水层位的吸水强度。引入概念吸水强度突进系数表征单层突进现象，即单层吸水强度与全井平均单层吸水强度的比值：
式中：λ为吸水强度突进系数，d为单个油层吸水强度，n为全井吸水层数。
将吸水强度和吸水强度突进系数作为评价是否是优势渗流层段的判别参数。根据计算得到的单元吸水强度和吸水强度突进系数分布规律，确定判别参数界限，将同时大于或等于这两个判别参数界限的识别为优势渗流井层，判识结果“I”级的为优势渗流井层。低效无效循环位于优势渗流井层中。
表4 单元优势渗流井层定量识别结果
步骤3、低效无效循环条带定量识别
对识别出的优势渗流井层，首先应用分层注采流动关系识别技术，判断优势渗流层段对应油水井存在连通流动关系，需同时满足五个判别条件。
判别条件一：注水井与采油井在某小层必须在同一砂体内。判别条件二：注水井与采油井间在某小层必须存在连通路径，流体仅能在砂体内流动，由于砂体的形态、砂体内可能有泥岩和断层等因素的影响，两口井之间可能没有连通路径。判别条件三：注水井与采油井在某小层必须同时生产。注水井与采油井必须在同一小层均射孔，采油井正在采油、注水井正在注水，才可能有流体在两井间流动。判别条件四：注水井与采油井在某小层必须为一线连通。在注水井与采油井之间的油水流动路径周围，如果存在一口其它注水井或采油井并同时在该小层生产，则表明这两口井间存在干扰井，即不是一线连通。干扰井的存在导致这两口井间的油水流动性大大降低，甚至没有流动，本次不考虑这种极弱的流动关系。判别条件五：注水井与采油井间在某小层压力传导必须合理。当两条注水井与采油井间的油水流动路线相交时，如果相交水流方向一致，则两条油水流动路线同时存在，否则只有一条油水流动路线存在，注水井的水流将会朝着阻力较小的采油井方向流动。
确定出优势渗流层段存在流动关系的油水井后，创建渗流阻力系数劈分注水量算法，计算出优势渗流井层各流动方向的渗流阻力系数，进而确定劈分注水比例系数，算得各注采方向的瞬时注水量和注水强度，判断出平面优势水流分布。
步骤1.以注水井为中心，计算第i口水井与第j口受效油井小层渗流阻力系数
ko＝k·kro(Sw)
kw＝k·krw(Sw)
式中：Rij为ij两点间小层渗流阻力系数；Lij为ij间井距，m；Hj为第j口油井的小层有效厚度，m；μw为水粘度，μo为油粘度，mPa·s；koj为第j口油井方向小层平均油相渗透率，kwj为第j口油井方向小层平均水相渗透率，mD。
步骤2.计算第i口注水井小层平面总渗流阻力
步骤3.计算第i口注水井第j口受效油井方向的小层注水比例系数
步骤4.根据第i口注水井小层注入量，计算第j口受效油井方向的瞬时注水量
Wij＝λ平面ij×Wi
式中：Wi为第i口水井小层日注水量。
步骤5.根据瞬时注水量，得到第j口受效油井方向的注水强度
Iij＝Wij/Hj
将瞬时注水量和注水强度作为优势渗流方向判别参数，依据分布规律，确定判别参数界限。将这两个判别参数同时大于或等于其界限的识别为优势渗流方向，识别结果“I”级的为优势渗流方向。
表5 单元优势渗流方向定量识别结果
识别出的优势渗流方向又分为低效、无效和有效循环，首先根据分流量方程计算各优势渗流方向的含水率：
根据区块低效和无效循环含水率界限，将各优势渗流方向的含水率与界限比较，识别出低效无效循环条带。
表6 单元低效无效循环含水率界限
表7 单元低效无效循环识别结果
上述实施例为本发明较佳的实施方式，但本发明的实施方式并不受上述实施例的限制，其他的任何未背离本发明的精神实质与原理下所作的改变、修饰、替代、组合、简化，均应为等效的置换方式，都包含在本发明的保护范围之内。
Claims (10)
1.一种特高含水后期低效无效循环条带定量判识方法，其特征在于，其包括以下步骤：
定量识别强注水井点，确定低效无效循环起源；
识别强注水井点的优势渗流层段，确定低效无效循环位置；
明确优势渗流层段平面注采流动方向，确定优势渗流方向判别参数及其界限，识别低效无效循环条带。
2.根据权利要求1所述的方法，其特征在于，应用模糊数学方法定量识别强注水井点，具体包括以下步骤：
筛选评价注水状况的表征参数，作为评价参数；设定评价参数的评价等级，选取模糊隶属度函数计算评价参数对各评价等级的隶属度，得到单注水井点模糊评价矩阵；应用模糊层次分析法构建评价参数判断矩阵，并一致性检验；应用特征根法求取判断矩阵特征向量，归一化得到评价因素的权重向量；权重向量与模糊评价矩阵进行模糊运算，并根据最大隶属度原则，确定各水井注水能力强弱，识别出强注水井点。
3.根据权利要求2所述的方法，其特征在于，通过分析油藏开发各动静态参数之间的相互制约关系，筛选能反映水井注水状况的渗透率、有效厚度、视吸水指数和单位厚度累注入量作为评价参数。
4.根据权利要求2所述的方法，其特征在于，评价等级设定为注水能力强和注水能力不强两级；根据各评价参数与水井注水强弱的影响关系，选取半梯形隶属度函数，计算评价参数对各评价等级的隶属度，得到单注水井点模糊评价矩阵。
5.根据权利要求1所述的方法，其特征在于，识别强注水井点的优势渗流层段，确定低效无效循环位置的方法包括以下步骤：依据岩心水洗特征、测井曲线特征、水淹特征及沉积相资料分析解释识别优势渗流井层；计算强注水井点各注水层位的吸水强度和吸水强度突进系数这两个判别参数，定量识别出优势渗流井层；低效无效循环位于识别出的优势渗流井层。
6.根据权利要求1所述的方法，其特征在于，优势渗流层段平面注采流动方向的确定需同时满足以下判别条件：
1)在同一砂体；
2)在同一时间生产；
3)存在连通路径；
4)流动路线上不存在其它井正在生产，为一线连通；
5)油水流动方向与其它油水流向一致，压力传导需合理。
7.根据权利要求1所述的方法，其特征在于，将瞬时注水量和注水强度作为优势渗流方向判别参数，根据这两个参数分布曲线拐点值，确定判别参数界限，同时大于或等于这两个参数判别界限的为优势渗流方向，将各优势渗流方向的含水率与区块低效无效循环含水率界限比较，识别出低效无效循环条带。
8.根据权利要求7所述的方法，其特征在于，针对优势渗流层段，以注水井为中心，根据动静态参数，采用渗流阻力系数劈分注水量算法计算出各受效油井方向的的渗流阻力系数，确定劈分注水比例系数，计算出各注采方向的瞬时注水量和注水强度。
9.根据权利要求8所述的方法，其特征在于，各注采方向的瞬时注水量和注水强度的计算具体包括以下步骤：
(1)计算第i口水井与第j口受效油井小层渗流阻力系数
ko＝k·kro(Sw)
kw＝k·krw(Sw)
(2)计算第i口注水井小层平面总渗流阻力
(3)计算第i口注水井第j口受效油井方向的小层注水比例系数
(4)根据第i口注水井小层注入量，计算第j口受效油井方向的瞬时注水量
Wij＝λ平面ij×Wi
(5)根据瞬时注水量，得到第j口受效油井方向的注水强度
Iij＝Wij/Hj
式中：Rij为jj两点间小层渗流阻力系数；Lij为jj间井距，m；Hj为第j口油井的小层有效厚度，m；μw为水粘度，μo为油粘度，mPa·s；koj为第j口油井方向小层平均油相渗透率，kwj为第j口油井方向小层平均水相渗透率，mD；Wi为第i口水井小层日注水量。
10.一种特高含水后期低效无效循环条带定量判识系统，其特征在于，具有用于执行根据权利要求1-9任一项所述识别方法的功能模块。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111156399.9A CN115874997A (zh) | 2021-09-29 | 2021-09-29 | 一种特高含水后期低效无效循环条带定量判识方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111156399.9A CN115874997A (zh) | 2021-09-29 | 2021-09-29 | 一种特高含水后期低效无效循环条带定量判识方法 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN115874997A true CN115874997A (zh) | 2023-03-31 |
Family
ID=85756554
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202111156399.9A Pending CN115874997A (zh) | 2021-09-29 | 2021-09-29 | 一种特高含水后期低效无效循环条带定量判识方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN115874997A (zh) |
-
2021
- 2021-09-29 CN CN202111156399.9A patent/CN115874997A/zh active Pending
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN112343587B (zh) | 一种特低渗透油藏优势渗流通道识别表征方法 | |
| CN104747185B (zh) | 非均质油藏储层综合分类评价方法 | |
| CN104790926B (zh) | 一种缝洞型油藏注水开发效果评价方法 | |
| CN110288258A (zh) | 一种高含水油藏剩余油挖潜方法 | |
| CN104832166A (zh) | 一种页岩气水平井初期产能预测方法 | |
| CN110717301B (zh) | 一种基于支持向量机算法的流动单元信息分类识别方法 | |
| CN105005712B (zh) | 灰岩含水层富水性评价方法 | |
| CN112132454B (zh) | 一种煤层顶板或底板含水层富水性综合评价方法 | |
| CN112100727A (zh) | 一种断层破碎带影响下的富水隧道突水的预警防控方法 | |
| US20230160304A1 (en) | Method and system for predicting relative permeability curve based on machine learning | |
| CN111444621A (zh) | 基于动静态渗流界面的高含水油藏流动单元划分方法 | |
| CN112647930A (zh) | 一种水平井油藏工程找水方法 | |
| CN117131971A (zh) | 基于Xgboost算法的储层优势渗流通道预测方法 | |
| CN115874997A (zh) | 一种特高含水后期低效无效循环条带定量判识方法 | |
| CN111335871A (zh) | 一种基于分层产能评价的查层补孔技术方法 | |
| CN112561279B (zh) | 一种高耗水层带的识别方法和系统 | |
| CN113589371B (zh) | 碳酸盐岩内断层的封闭性评价方法 | |
| CN111155980A (zh) | 一种水流优势通道识别方法及装置 | |
| CN113032996B (zh) | 一种用于低渗裂缝型油藏水平井的水窜通道识别方法 | |
| CN109973068B (zh) | 油藏注水诱导裂缝的识别方法及装置 | |
| CN115705452A (zh) | 整装砂岩油藏开发中后期的新型采收率预测方法 | |
| CN115146427A (zh) | 一种计算缝洞型油藏调流道堵剂用量的方法及系统 | |
| CN112632864B (zh) | 一种动静结合识别高耗水层带的矿场方法 | |
| CN117449846B (zh) | 一种致密砂岩气藏产水来源综合诊断方法 | |
| Sha et al. | Research on Accurate Identification Method for Water Flooding Channel in Low Permeability Fractured Reservoirs |
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 |