CN112102101A - 一种水驱油藏井组调堵决策方法 - Google Patents
一种水驱油藏井组调堵决策方法 Download PDFInfo
- Publication number
- CN112102101A CN112102101A CN202010937743.7A CN202010937743A CN112102101A CN 112102101 A CN112102101 A CN 112102101A CN 202010937743 A CN202010937743 A CN 202010937743A CN 112102101 A CN112102101 A CN 112102101A
- Authority
- CN
- China
- Prior art keywords
- water
- seepage channel
- dominant
- water content
- dominant seepage
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 44
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 225
- 230000035699 permeability Effects 0.000 claims abstract description 64
- 238000002347 injection Methods 0.000 claims abstract description 63
- 239000007924 injection Substances 0.000 claims abstract description 63
- 238000004519 manufacturing process Methods 0.000 claims abstract description 62
- 238000011161 development Methods 0.000 claims abstract description 43
- 238000004364 calculation method Methods 0.000 claims abstract description 33
- 239000011148 porous material Substances 0.000 claims abstract description 26
- 239000003795 chemical substances by application Substances 0.000 claims abstract description 24
- 238000013461 design Methods 0.000 claims abstract description 7
- 238000011160 research Methods 0.000 claims abstract description 6
- 238000012216 screening Methods 0.000 claims abstract description 5
- 239000012071 phase Substances 0.000 claims description 22
- 238000005325 percolation Methods 0.000 claims description 15
- 230000008569 process Effects 0.000 claims description 10
- 230000000694 effects Effects 0.000 claims description 9
- 206010011416 Croup infectious Diseases 0.000 claims description 8
- 201000010549 croup Diseases 0.000 claims description 8
- 230000033228 biological regulation Effects 0.000 claims description 6
- 239000012530 fluid Substances 0.000 claims description 6
- 239000002245 particle Substances 0.000 claims description 6
- 230000000704 physical effect Effects 0.000 claims description 6
- 230000001105 regulatory effect Effects 0.000 claims description 6
- 238000012407 engineering method Methods 0.000 claims description 5
- 238000004422 calculation algorithm Methods 0.000 claims description 4
- 238000009795 derivation Methods 0.000 claims description 4
- 238000009499 grossing Methods 0.000 claims description 4
- 239000008346 aqueous phase Substances 0.000 claims description 3
- 150000001875 compounds Chemical class 0.000 claims description 3
- 230000001186 cumulative effect Effects 0.000 claims description 3
- 230000007423 decrease Effects 0.000 claims description 3
- 238000005065 mining Methods 0.000 claims description 3
- 238000005457 optimization Methods 0.000 claims description 3
- 230000037361 pathway Effects 0.000 claims description 3
- 238000004088 simulation Methods 0.000 claims description 3
- 239000000243 solution Substances 0.000 claims description 3
- 239000000126 substance Substances 0.000 claims description 3
- 235000020681 well water Nutrition 0.000 claims description 3
- 239000002349 well water Substances 0.000 claims description 3
- 230000002902 bimodal effect Effects 0.000 claims description 2
- 238000012360 testing method Methods 0.000 description 6
- 230000008901 benefit Effects 0.000 description 4
- 238000002474 experimental method Methods 0.000 description 2
- 238000012544 monitoring process Methods 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 230000001276 controlling effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 239000003814 drug Substances 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 239000004005 microsphere Substances 0.000 description 1
- 238000013508 migration Methods 0.000 description 1
- 230000005012 migration Effects 0.000 description 1
- 238000011158 quantitative evaluation Methods 0.000 description 1
- 239000000700 radioactive tracer Substances 0.000 description 1
- 230000000087 stabilizing effect Effects 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 239000002699 waste material Substances 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/02—Agriculture; Fishing; Forestry; Mining
-
- 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0639—Performance analysis of employees; Performance analysis of enterprise or organisation operations
- G06Q10/06393—Score-carding, benchmarking or key performance indicator [KPI] analysis
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T11/00—2D [Two Dimensional] image generation
- G06T11/20—Drawing from basic elements, e.g. lines or circles
- G06T11/203—Drawing of straight lines or curves
-
- 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
- Engineering & Computer Science (AREA)
- Business, Economics & Management (AREA)
- Human Resources & Organizations (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Strategic Management (AREA)
- Economics (AREA)
- Development Economics (AREA)
- Tourism & Hospitality (AREA)
- General Business, Economics & Management (AREA)
- Marketing (AREA)
- Educational Administration (AREA)
- Entrepreneurship & Innovation (AREA)
- General Engineering & Computer Science (AREA)
- Agronomy & Crop Science (AREA)
- General Health & Medical Sciences (AREA)
- Mining & Mineral Resources (AREA)
- Marine Sciences & Fisheries (AREA)
- Primary Health Care (AREA)
- Animal Husbandry (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Geometry (AREA)
- Evolutionary Computation (AREA)
- Computer Hardware Design (AREA)
- Game Theory and Decision Science (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种水驱油藏井组调堵决策方法：对于研究区块井组的生产井生产动态数据消除噪点计算含水率导数，绘制含水率、含水率导数与时间的关系曲线，判断是否存在优势渗流通道发育，确定优势渗流通道及正常储层的见水时间；求取优势渗流通道与正常储层渗透率之比，得到优势渗流通道渗透率；求取优势渗流通道的平均孔吼半径；求取优势渗流通道体积；利用优势渗流通道与正常储层渗透率之比筛选调堵井组，利用优势渗流通道的平均孔吼半径和优势渗流通道体积，指导调堵剂类型和注入量的设计。本发明实现优势渗流通道是否发育的快速判定，以及优势渗流通道孔吼半径和发育规模的快速计算，可以进行水驱油藏井组调堵决策。
Description
技术领域
本发明涉及油气田开发研究领域，更具体地说，是涉及一种水驱油藏井组调堵决策方法。
背景技术
油水井堵水调剖是优势渗流通道发育的水驱油稳油控水、提高水驱开发效果的重要手段。但其成功的关键在于准确选取发育优势渗流通道的井组，不适当的措施则会抑制生产井产能、注水井的注入能力造成开发效果变差。进行调堵措施前，需要先获取两方面认识：一是判断是否发育优势渗流通道，它是决策是否有必要进行调堵工作的依据；二是定量描述优势渗流通道的参数，这是决定油田现场堵剂用量的关键参数。此处优势渗流通道的关键参数主要包括两个：一是优势渗流通道的孔吼尺寸，孔吼尺寸决定后续调剖调驱等措施注入段塞的浓度或者颗粒大小，注入段塞浓度过低、颗粒过小无法形成有效封堵，浓度过高、颗粒过大会造成注入困难，近井堵塞影响堵剂运移范围和措施效果；二是优势渗流通道的发育规模，发育规模决定了堵剂用量，用量过少难以实现深部封堵，造成堵水有效期较短或难以实现控水效果，用量过多注入困难且造成浪费，影响项目经济性。
目前业界常用优势渗流通道识别方法主要有：测井方法、井间动态监测方法、试井方法、取心方法和示踪剂方法等。但测井、取心等方法仅仅反应测试点处，无法反应优势渗流通道整体发育规模和井间情况，更无法准确计算优势渗流通道的发育规模及井间尺寸。示踪剂及试井依赖于专门的测试，周期长、费用高。现有的井间动态监测方法大多聚焦于对是否发育优势渗流通道做定性判断，对其参数的定量描述欠缺，并且上述现有技术往往受矿场测试少、费用高、周期长、占用海上平台作业空间等限制，应用受限，难以满足海上油田开发生产的实际需求，需要寻求更为便捷、准确的油藏工程方法实现优势渗流通道的定量描述，并建立一套水驱油藏井组调堵决策方法。
发明内容
为了解决现有技术在老油田优势渗流通道定量描述、井组调堵决策方面的不足，本发明基于油藏工程方法，矿场基础静态参数、生产动态参数和岩心实验分析结果，提出一种水驱油藏井组调堵决策方法，实现优势渗流通道是否发育的快速判定，以及优势渗流通道孔吼半径和发育规模的快速计算，可以进行水驱油藏井组调堵决策；所用参数简单易获取，可以有效指导矿场实践。
本发明的目的是通过以下技术方案实现的。
本发明水驱油藏井组调堵决策方法，包括以下过程：
第一步：对于研究区块井组的生产井生产动态数据，采用平滑算法消除噪点计算含水率导数，绘制含水率、含水率导数与时间的关系曲线，根据含水率导数曲线是否呈现双峰特征，判断是否存在优势渗流通道发育，确定优势渗流通道及正常储层的见水时间；
第二步：根据以B-L方程为基础推导出的见水时间和物性关系式，采用比值的方式消除参数，最终得到含水率导数特征值与物性关计算式，求取优势渗流通道与正常储层渗透率之比，进而得到优势渗流通道渗透率；
第三步：利用第二步所求得的优势渗流通道渗透率，基于所在矿区岩心实验结果，求取优势渗流通道的平均孔吼半径；
第四步：利用第二步所求得的优势渗流通道渗透率，结合已知开发参数求取优势渗流通道体积；
第五步：利用第二步中优势渗流通道与正常储层渗透率之比筛选调堵井组，利用第三步求得的优势渗流通道的平均孔吼半径和第四步求得的优势渗流通道体积，指导调堵剂类型和注入量的设计。
第一步中计算含水率导数，绘制含水率、含水率导数与时间的关系曲线，根据含水率导数曲线是否呈现双峰特征，判断是否存在优势渗流通道发育，确定优势渗流通道及正常储层的见水时间，具体过程如下:
为方便叙述计算过程区分优势渗流通道和正常储层参数，将参数附加下标i，i＝1时代表优势渗流通道，i＝2时代表正常储层；
采用如下公式对研究区块井组的生产井生产动态数据降噪：
式中：fw、f′w分别为生产井综合含水率、含水率导数，f；fwj-1、fwj、fwj+1分别为第j-1、j、j+1时刻的含水率，f；tj-1、tj、tj+1分别为第j-1、j、j+1时刻的生产时间，d；下标j代表时间步长；
存在优势渗流通道的井组含水率实际为正常储层与优势渗流通道储层含水率的叠加反映；
式中：Qi为优势渗流通道或正常储层的产液量，m3/s；fwi为优势渗流通道或正常储层的含水率，f；
正常储层见水后含水率上升平滑，含水率导数曲线在见水时刻出现峰值，随后单调递减；发育优势渗流通道的井：注入水首先沿高渗优势渗流通道到达生产井，即第一见水时刻，这也是含水率导数曲线的第一个峰值对应的时间；随后注入水沿渗透率较低正常储层到达生产井，即第二见水时刻，由于水驱前缘含水率是一个跃变，这也是对应含水率导数曲线的第二个峰值；发育优势渗流通道的井含水率曲线特征为迅速上升后有一个平缓段，然后再上升，含水率导数曲线的特征是出现双峰；通过含水率和含水率导数曲线，确定优势渗流通道是否存发育，并确定优势渗流通道或正常储层的见水时间Ti值。
第二步中优势渗流通道渗透率的计算方法及理论推导包括：
(ⅰ)建立定压条件下的双层非活塞式水驱油模型
优势渗流通道或正常储层的渗流符合非活塞式水驱油，水驱前缘位置为：
式中：xfi为优势渗流通道或正常储层的水驱前缘位置，m；swf为水驱前缘含水饱和度，f；f′w(swf)为水驱前缘饱和度下的含水率导数，f；
为优势渗流通道或正常储层的孔隙度，f；
为从0到t时刻优势渗流通道或正常储层的累积注入量，m3；Ai为优势渗流通道或正常储层的渗流截面积，m2；
依据B-L理论，两相区某一点a与注入端压差：
式中：Krw为水相相对渗透率，10-3μm2；Δpai为优势渗流通道或正常储层的某点a与注入端之间压差，MPa；xa为a位置与原点距离，m；pai为优势渗流通道或正常储层x＝a位置处压力，MPa；pinj为注入井底压力，MPa；v为总渗流速度，m/s；Ki为优势渗流通道或正常储层的储层渗透率，10-3μm2；fw(Sw)为含水饱和度Sw下的含水率，f；μw为水相黏度，mPa·s；
结合式(3)、(4)，注入井至水驱前缘任意处压差与含水饱和度关系：
式中，
式中：f′w(sw)为含水饱和度Sw下的含水率导数，f；swa为x＝a位置的含水饱和度，当swa＝swf，则得到F(Swf)＝C(常数)；联立上述两个方程积分整理后，得到0～t时刻两相区的累积注入量为：
式中：Δp为注入端与产出端之间压差，MPa；t为生产时间，s；L为注采井距，m；μo为油相黏度，mPa·s；
上式与等饱和度面移动方程结合：
水驱前缘到达采油井，xfi＝L，见水时间由下式求得：
(ⅱ)计算优势渗流通道渗透率
采用比值方法约去相关参数，可以得到优势渗流通道与正常储层的见水时间之比为：
令
进行简化，则优势渗流通道的渗透率计算式：
第三步中优势渗流通道的平均孔吼半径计算方法：根据油田已有的孔喉半径与渗透率关系，将第二步中获得的优势渗流通道渗透率带入其中，获取优势渗流通道的平均孔吼半径。
第四步中优势渗流通道发育体积计算方法：
在T2时刻，由于优势渗流通道早已见水，完全处于两相渗流区，而正常储层的水驱前缘刚刚抵达生产井，若将该时刻含水率取跃变前一刻的值，此时fw2(T2)＝0，因此得到全井含水率满足：
式中：fwT2为T2时刻全井综合含水率值，f；fw1(T2)与fw2(T2)分别为T2时刻优势渗流通道和正常储层的含水率，f；Q1(T2)与Q2(T2)分别为T2时刻优势渗流通道和正常储层的产液量，f；
计算前先对产量方程进行处理变形：
由等饱和度面移动方程取微分得：
含水率方程：
式中：Kro、Krw分别为油相、水相相对渗透率，f；
则产量方程：
其中，
式中：fw(swf)为水驱前缘饱和度下的含水率，f；F′w(swei)为优势渗流通道或正常储层的出口端含水饱和度下的含水率导数，f；
本计算方法的具体步骤为：
a)首先已知注采井间渗流截面A，假设其中优势渗流通道的渗流截面A1，则其中正常储层渗流截面A2＝A-A1；
b)分别计算初始时刻各层的注入量，设置时间步长，采用迭代计算至T2时刻，若此时含水率等于FwT2，则保留A1，否则修改A1值重新计算；
c)利用计算得到的A1值计算优势渗流通道的发育体积：
式中：V1为优势渗流通道发育体积，m3；V为注采井间波及体积，m3。
第五步中筛选调堵井组、指导调堵剂类型和注入量的设计：
当优势渗流通道与正常储层渗透率比值
视为优势渗流通道弱发育，优势渗流通道与正常储层渗透率比值
视为优势渗流通道强发育，建议进行调堵措施；调堵剂地下颗粒半径应达到优势渗流通道平均孔吼半径的1/10～1/3，以达到有效封堵优势渗流通道的目的；调堵剂注入量根据计算得到的优势渗流通道体积进行设计：
式中：V调堵为调堵段塞用量，m3；V为注采井间波及体积，m3；δ为注入孔隙体积倍数，f；
注入孔隙体积倍数δ，可以根据具体调堵剂，采用油藏工程方法进行模拟优化，通过模拟不同注入孔隙体积倍数条件下堵剂用量对应的措施效果，优化出最优的堵剂用量。
与现有技术相比，本发明的技术方案所带来的有益效果是：
(1)本发明针存在优势渗流通道的注采井见水后含水率的叠加跃变现象，利用求导方法放大差异，实现了优势渗流通道发育情况的判别。
(2)本发明对注水开发的油田提供了一套优势渗流通道参数定量计算的方法，包括优势渗流通道的孔吼尺寸和发育体积，该方法依据常规的动态参数和矿场基本数据，以非活塞水驱油理论为基础，建立优势渗流通道与正常储层各项参数之间的关系式，进而联立求取优势渗流通道参数。
(3)本发明基于得到的优势渗流通道参数可以进行油水井调堵决策，同时根据优势渗流通道孔吼半径尺寸和发育体积指导调堵措施的设计。
(4)基于以上三点，本发明进一步提高了优势渗流通道的计算精度、油水井调堵措施效果，并且无需额外测试工作，节省了费用及时间，实现了优势渗流通道的快速定量评价、调堵措施定量决策，为油田开发提供有力保障。
附图说明
图1是本发明水驱油藏井组调堵决策方法的流程图；
图2为本发明实施例提供的含水率导数特征值确定方法图；
图3为本发明实施例提供的Q油田岩心平均孔喉半径与渗透率关系图；
图4为本发明实施例提供的优势渗流通道发育体积计算程序图；
图5为本发明实施例提供的X井生产曲线图；
图6为本发明实施例提供的X井含水率及导数曲线图；
图7为本发明实施例提供的X井组调驱后的效果曲线。
具体实施方式
下面结合附图对本发明作进一步的描述。
本发明水驱油藏井组调堵决策方法，如图1所示，具体实现过程如下：
第一步：对于研究区块井组的生产井生产动态数据，采用平滑算法消除噪点计算含水率导数，绘制含水率、含水率导数与时间的关系曲线，根据含水率导数曲线是否呈现双峰特征，判断是否存在优势渗流通道发育，确定优势渗流通道及正常储层的见水时间。具体过程如下：
为方便叙述计算过程区分优势渗流通道和正常储层参数，将见水时间、渗透率等参数附加下标i，i＝1时代表优势渗流通道，i＝2时代表正常储层。
为了保障含水率及其导数曲线具有较好的光滑性和比进度，提高准确性和可靠性，采用如下公式对研究区块井组的生产井生产动态数据降噪：
式中：fw、f′w分别为生产井综合含水率、含水率导数，f；fwj-1、fwj、fwj+1分别代表第j-1、j、j+1时刻的含水率，f；tj-1、tj、tj+1分别代表第j-1、j、j+1时刻的生产时间，d。下标j代表时间步长。
存在优势渗流通道的井组含水率实际为正常储层与优势渗流通道储层含水率的叠加反映。
式中：Qi为优势渗流通道或正常储层的产液量，m3/s；fwi为优势渗流通道或正常储层的含水率，f。
正常储层见水后含水率上升平滑，含水率导数曲线在见水时刻出现峰值，随后单调递减。发育优势渗流通道的井：注入水首先沿高渗优势渗流通道到达生产井，即第一见水时刻，这也是含水率导数曲线的第一个峰值对应的时间；随后注入水沿渗透率较低正常储层到达生产井，即第二见水时刻，由于水驱前缘含水率是一个跃变，这也是对应含水率导数曲线的第二个峰值。发育优势渗流通道的井含水率曲线特征为迅速上升后有一个平缓段，然后再上升，含水率导数曲线的特征是出现双峰。如图2所示，通过含水率和含水率导数曲线，可确定优势渗流通道是否存发育，并确定Ti值，Ti为优势渗流通道或正常储层的见水时间，s。
第二步：根据以B-L方程为基础推导出的见水时间和物性关系式，采用比值的方式消除压差、井距等参数，最终得到含水率导数特征值与物性关计算式，求取优势渗流通道与正常储层渗透率之比，进而得到优势渗流通道渗透率。
其中，优势渗流通道渗透率的计算方法及理论推导包括：
(ⅰ)建立定压条件下的双层非活塞式水驱油模型
优势渗流通道或正常储层的渗流符合非活塞式水驱油，水驱前缘位置为：
式中：xfi为优势渗流通道或正常储层的水驱前缘位置，m；swf为水驱前缘含水饱和度，f；f′w(swf)为水驱前缘饱和度下的含水率导数，f；
为优势渗流通道或正常储层的孔隙度，f；
为从0到t时刻优势渗流通道或正常储层的累积注入量，m3；Ai为优势渗流通道或正常储层的渗流截面积，m2；
依据B-L理论，两相区某一点a与注入端压差：
式中：Krw为水相相对渗透率，10-3μm2；Δpai为优势渗流通道或正常储层的某点a与注入端之间压差，MPa；xa为a位置与原点距离，m；pai为优势渗流通道或正常储层x＝a位置处压力，MPa；pinj为注入井底压力，MPa；v为总渗流速度，m/s；Ki为优势渗流通道或正常储层的储层渗透率，10-3μm2；fw(Sw)为含水饱和度Sw下的含水率，f；μw为水相黏度，mPa·s；
结合式(3)、(4)，注入井至水驱前缘任意处压差与含水饱和度关系：
式中，
式中：f′w(sw)为含水饱和度Sw下的含水率导数，f；swa为x＝a位置的含水饱和度，当swa＝swf，则得到F(Swf)＝C(常数)。联立上述两个方程积分整理后，可以得到0～t时刻两相区的累积注入量为：
式中：Δp为注入端与产出端之间压差，MPa；t为生产时间，s；L为注采井距，m；μo为油相黏度，mPa·s；
上式与等饱和度面移动方程结合：
水驱前缘到达采油井，xfi＝L，见水时间可以由下式求得：
(ⅱ)计算优势渗流通道渗透率
为了消除等参数波动及计量带来的误差并进一步简化计算，采用比值方法约去相关参数，可以得到优势渗流通道与正常储层的见水时间之比为：
发育优势渗流通道后孔隙度并未发生变化，主要是孔吼结构导致的渗透率变化，因此可以令
进行简化，则优势渗流通道的渗透率计算式：
第三步：利用第二步所求得的优势渗流通道渗透率，基于所在矿区岩心实验结果(见图3)，求取优势渗流通道的平均孔吼半径。
其中，优势渗流通道的平均孔吼半径计算方法：根据油田已有的孔喉半径与渗透率关系，将第二步中获得的优势渗流通道渗透率带入其中，获取优势渗流通道的平均孔吼半径。
第四步：利用第二步所求得的优势渗流通道渗透率，结合相渗、井距、注采井间波及面积等已知开发参数求取优势渗流通道体积。
其中，优势渗流通道发育体积计算方法：
在T2时刻，由于优势渗流通道早已见水，完全处于两相渗流区，而正常储层的水驱前缘刚刚抵达生产井，若将该时刻含水率取跃变前一刻的值，此时fw2(T2)＝0，因此可以得到全井含水率满足：
式中：fwT2为T2时刻全井综合含水率值，f；fw1(T2)与fw2(T2)分别为T2时刻优势渗流通道和正常储层的含水率，f；Q1(T2)与Q2(T2)分别为T2时刻优势渗流通道和正常储层的产液量，f。
上述方程中，Qi在每个时刻都是变值，两类储层的渗流截面之和即为井组渗流截面，因此优势渗流通道的渗流截面采用迭代求解。由于优势渗流通道一般发育在正常储层底部，注采井间波及体积已知，进而可以得出优势渗流通道的发育体积。
计算前先对产量方程进行处理变形：
由等饱和度面移动方程取微分得：
含水率方程：
式中：Kro、Krw分别为油相、水相相对渗透率，f，均与含水饱和度有关。
则产量方程：
其中，
式中：fw(swf)为水驱前缘饱和度下的含水率，f；f′w(swei)为优势渗流通道或正常储层的出口端含水饱和度下的含水率导数，f。
本计算方法的计算程序见图4，具体步骤为：
a)首先已知注采井间渗流截面A，假设其中优势渗流通道的渗流截面A1，则其中正常储层渗流截面A2＝A-A1；
b)分别计算初始时刻各层的注入量，设置时间步长，采用迭代计算至T2时刻，若此时含水率等于fwT2，则保留A1，否则修改A1值重新计算；
c)利用计算得到的A1值计算优势渗流通道的发育体积：
式中：V1为优势渗流通道发育体积，m3；V为注采井间波及体积，m3。
第五步：利用第二步中优势渗流通道与正常储层渗透率之比筛选调堵井组，利用第三步求得的优势渗流通道的平均孔吼半径和第四步求得的优势渗流通道体积，指导调堵剂类型和注入量的设计。具体过程如下：
当优势渗流通道与正常储层渗透率比值
视为优势渗流通道弱发育，优势渗流通道与正常储层渗透率比值
视为优势渗流通道强发育，建议进行调堵措施。调堵剂地下颗粒半径应达到优势渗流通道平均孔吼半径的1/10～1/3，以达到有效封堵优势渗流通道的目的。调堵剂注入量根据计算得到的优势渗流通道体积进行设计：
式中：V调堵为调堵段塞用量，m3；V为注采井间波及体积，m3；δ为注入孔隙体积倍数，m3。
注入孔隙体积倍数δ，可以根据具体调堵剂，采用油藏工程方法进行模拟优化，通过模拟不同注入孔隙体积倍数条件下堵剂用量对应的措施效果，优化出最优的堵剂用量。
如图5所示为本发明实施例提供的X井生产曲线图，利用该井曲线进行平滑处理绘制含水率及导数曲线图。
图6为本发明实施例提供的X井含水率及导数曲线图。含水率导数曲线出现明显偏折，说明该井发育优势渗流通道，依据本发明算法及区块基本参数计算结果该井组的优势渗流通道参数如表1所示。解释成果指导对应注水井采用适用于6000～8000mD的纳米微球药剂进行调驱作业，如图7所示调驱后X井含水率大幅下降，产油量大幅增加，证实对于优势渗流通道的预测和计算较为准确，可以有效指导现场实践。
表1 X井优势渗流通道参数解释成果表
尽管上面结合附图对本发明进行了描述，但本发明并不局限于上述，上述的具体实施方式仅仅是示意性的，而不是限制性的，本领域的普通技术人员在本发明的启示下，还可以做出很多形式，这些均属于本发明的保护之内。
Claims (6)
1.一种水驱油藏井组调堵决策方法，其特征在于，包括以下过程：
第一步：对于研究区块井组的生产井生产动态数据，采用平滑算法消除噪点计算含水率导数，绘制含水率、含水率导数与时间的关系曲线，根据含水率导数曲线是否呈现双峰特征，判断是否存在优势渗流通道发育，确定优势渗流通道及正常储层的见水时间；
第二步：根据以B-L方程为基础推导出的见水时间和物性关系式，采用比值的方式消除参数，最终得到含水率导数特征值与物性关计算式，求取优势渗流通道与正常储层渗透率之比，进而得到优势渗流通道渗透率；
第三步：利用第二步所求得的优势渗流通道渗透率，基于所在矿区岩心实验结果，求取优势渗流通道的平均孔吼半径；
第四步：利用第二步所求得的优势渗流通道渗透率，结合已知开发参数求取优势渗流通道体积；
第五步：利用第二步中优势渗流通道与正常储层渗透率之比筛选调堵井组，利用第三步求得的优势渗流通道的平均孔吼半径和第四步求得的优势渗流通道体积，指导调堵剂类型和注入量的设计。
2.根据权利要求1所述的水驱油藏井组调堵决策方法，其特征在于，第一步中计算含水率导数，绘制含水率、含水率导数与时间的关系曲线，根据含水率导数曲线是否呈现双峰特征，判断是否存在优势渗流通道发育，确定优势渗流通道及正常储层的见水时间，具体过程如下:
为方便叙述计算过程区分优势渗流通道和正常储层参数，将参数附加下标i，i＝1时代表优势渗流通道，i＝2时代表正常储层；
采用如下公式对研究区块井组的生产井生产动态数据降噪：
式中：fw、f′w分别为生产井综合含水率、含水率导数，f；fwj-1、fwj、fwj+1分别为第j-1、j、j+1时刻的含水率，f；tj-1、tj、tj+1分别代表第j-1、j、j+1时刻的生产时间，d；下标j代表时间步长；
存在优势渗流通道的井组含水率实际为正常储层与优势渗流通道储层含水率的叠加反映；
式中：Qi为优势渗流通道或正常储层的产液量，m3/s；fwi为优势渗流通道或正常储层的含水率，f；
正常储层见水后含水率上升平滑，含水率导数曲线在见水时刻出现峰值，随后单调递减；发育优势渗流通道的井：注入水首先沿高渗优势渗流通道到达生产井，即第一见水时刻，这也是含水率导数曲线的第一个峰值对应的时间；随后注入水沿渗透率较低正常储层到达生产井，即第二见水时刻，由于水驱前缘含水率是一个跃变，这也是对应含水率导数曲线的第二个峰值；发育优势渗流通道的井含水率曲线特征为迅速上升后有一个平缓段，然后再上升，含水率导数曲线的特征是出现双峰；通过含水率和含水率导数曲线，确定优势渗流通道是否存发育，并确定优势渗流通道或正常储层的见水时间Ti值。
3.根据权利要求1所述的水驱油藏井组调堵决策方法，其特征在于，第二步中优势渗流通道渗透率的计算方法及理论推导包括：
(ⅰ)建立定压条件下的双层非活塞式水驱油模型
优势渗流通道或正常储层的渗流符合非活塞式水驱油，水驱前缘位置为：
式中：xfi为优势渗流通道或正常储层的水驱前缘位置，m；swf为水驱前缘含水饱和度，f；f′w(swf)为水驱前缘饱和度下的含水率导数，f；
为优势渗流通道或正常储层的孔隙度，f；
为从0到t时刻优势渗流通道或正常储层的累积注入量，m3；Ai为优势渗流通道或正常储层的渗流截面积，m2；
依据B-L理论，两相区某一点a与注入端压差：
式中：Krw为水相相对渗透率，10-3μm2；Δpai为优势渗流通道或正常储层的某点a与注入端之间压差，MPa；xa为a位置与原点距离，m；pai为优势渗流通道或正常储层x＝a位置处压力，MPa；pinj为注入井底压力，MPa；v为总渗流速度，m/s；Ki为优势渗流通道或正常储层的储层渗透率，10-3μm2；fw(Sw)为含水饱和度Sw下的含水率，f；μw为水相黏度，mPa·s；
结合式(3)、(4)，注入井至水驱前缘任意处压差与含水饱和度关系：
式中，
式中：f′w(sw)为含水饱和度Sw下的含水率导数，f；swa为x＝a位置的含水饱和度，当swa＝swf，则得到F(Swf)＝C(常数)；联立上述两个方程积分整理后，得到0～t时刻两相区的累积注入量为：
式中：Δp为注入端与产出端之间压差，MPa；t为生产时间，s；L为注采井距，m；μo为油相黏度，mPa·s；
上式与等饱和度面移动方程结合：
水驱前缘到达采油井，xfi＝L，见水时间由下式求得：
(ⅱ)计算优势渗流通道渗透率
采用比值方法约去相关参数，可以得到优势渗流通道与正常储层的见水时间之比为：
令
进行简化，则优势渗流通道的渗透率计算式：
4.根据权利要求1所述的水驱油藏井组调堵决策方法，其特征在于，第三步中优势渗流通道的平均孔吼半径计算方法：根据油田已有的孔喉半径与渗透率关系，将第二步中获得的优势渗流通道渗透率带入其中，获取优势渗流通道的平均孔吼半径。
5.根据权利要求1所述的水驱油藏井组调堵决策方法，其特征在于，第四步中优势渗流通道发育体积计算方法：
在T2时刻，由于优势渗流通道早已见水，完全处于两相渗流区，而正常储层的水驱前缘刚刚抵达生产井，若将该时刻含水率取跃变前一刻的值，此时fw2(T2)＝0，因此得到全井含水率满足：
式中：fwT2为T2时刻全井综合含水率值，f；fw1(T2)与fw2(T2)分别为T2时刻优势渗流通道和正常储层的含水率，f；Q1(T2)与Q2(T2)分别为T2时刻优势渗流通道和正常储层的产液量，f；
计算前先对产量方程进行处理变形：
由等饱和度面移动方程取微分得：
含水率方程：
式中：Kro、Krw分别为油相、水相相对渗透率，f；
则产量方程：
其中，
式中：fw(swf)为水驱前缘饱和度下的含水率，f；fw′(swei)为优势渗流通道或正常储层的出口端含水饱和度下的含水率导数，f；
本计算方法的具体步骤为：
a)首先已知注采井间渗流截面A，假设其中优势渗流通道的渗流截面A1，则其中正常储层渗流截面A2＝A-A1；
b)分别计算初始时刻各层的注入量，设置时间步长，采用迭代计算至T2时刻，若此时含水率等于fwT2，则保留A1，否则修改A1值重新计算；
c)利用计算得到的A1值计算优势渗流通道的发育体积：
式中：V1为优势渗流通道发育体积，m3；V为注采井间波及体积，m3。
6.根据权利要求1所述的水驱油藏井组调堵决策方法，其特征在于，第五步中筛选调堵井组、指导调堵剂类型和注入量的设计：
当优势渗流通道与正常储层渗透率比值
视为优势渗流通道弱发育，优势渗流通道与正常储层渗透率比值
视为优势渗流通道强发育，建议进行调堵措施；调堵剂地下颗粒半径应达到优势渗流通道平均孔吼半径的1/10～1/3，以达到有效封堵优势渗流通道的目的；调堵剂注入量根据计算得到的优势渗流通道体积进行设计：
式中：V调堵为调堵段塞用量，m3；V为注采井间波及体积，m3；δ为注入孔隙体积倍数，f；
注入孔隙体积倍数δ，可以根据具体调堵剂，采用油藏工程方法进行模拟优化，通过模拟不同注入孔隙体积倍数条件下堵剂用量对应的措施效果，优化出最优的堵剂用量。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010937743.7A CN112102101B (zh) | 2020-09-08 | 2020-09-08 | 一种水驱油藏井组调堵决策方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010937743.7A CN112102101B (zh) | 2020-09-08 | 2020-09-08 | 一种水驱油藏井组调堵决策方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112102101A true CN112102101A (zh) | 2020-12-18 |
| CN112102101B CN112102101B (zh) | 2023-01-24 |
Family
ID=73751219
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010937743.7A Active CN112102101B (zh) | 2020-09-08 | 2020-09-08 | 一种水驱油藏井组调堵决策方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112102101B (zh) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113586046A (zh) * | 2021-06-23 | 2021-11-02 | 中海油能源发展股份有限公司 | 一种跟踪评价水平井调驱效果的方法 |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106022626A (zh) * | 2016-05-27 | 2016-10-12 | 中国海洋石油总公司 | 一种水驱油藏定位调堵水驱拟前缘的计算方法 |
| CN109948272A (zh) * | 2019-03-27 | 2019-06-28 | 长江大学 | 基于井间连通性的调堵动态预测方法和系统 |
-
2020
- 2020-09-08 CN CN202010937743.7A patent/CN112102101B/zh active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106022626A (zh) * | 2016-05-27 | 2016-10-12 | 中国海洋石油总公司 | 一种水驱油藏定位调堵水驱拟前缘的计算方法 |
| CN109948272A (zh) * | 2019-03-27 | 2019-06-28 | 长江大学 | 基于井间连通性的调堵动态预测方法和系统 |
Non-Patent Citations (1)
| Title |
|---|
| 黄斌等: "基于模糊综合评判法的优势渗流通道定量识别方法研究", 《数学的实践与认识》 * |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113586046A (zh) * | 2021-06-23 | 2021-11-02 | 中海油能源发展股份有限公司 | 一种跟踪评价水平井调驱效果的方法 |
| CN113586046B (zh) * | 2021-06-23 | 2023-07-28 | 中海油能源发展股份有限公司 | 一种跟踪评价水平井调驱效果的方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112102101B (zh) | 2023-01-24 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Jing et al. | Mechanism simulation of oil displacement by imbibition in fractured reservoirs | |
| Furui et al. | A rigorous formation damage skin factor and reservoir inflow model for a horizontal well | |
| CN111236908A (zh) | 一种适用在低渗透致密气藏中的多段压裂水平井产能预测模型及产能敏感性分析的方法 | |
| US11692426B2 (en) | Method and apparatus for determining integrated exploitation approach of shale and adjacent oil reservoirs | |
| CN106522928B (zh) | 一种酸化压裂后停泵测井口压降不稳定试井方法 | |
| CN104060985A (zh) | 一种层状油藏调剖堵水堵剂进入深度测试方法及系统 | |
| CN112102101B (zh) | 一种水驱油藏井组调堵决策方法 | |
| CN114880962B (zh) | 基于地层渗流-井筒管流耦合的致密气藏储气库单井注采能力一体化分析方法 | |
| Wang et al. | The dynamic convection-diffusion effect of CO2 Huff-n-Puff on Pre-Darcy flow behavior | |
| Moradi et al. | Near-well simulation of oil production from a horizontal well with ICD and AICD completions in the Johan Sverdrup field using OLGA/ROCX | |
| CN107355200B (zh) | 一种纳微米颗粒分散体系改善水驱选井方法 | |
| Dietrich | Relative permeability during cyclic steam stimulation of heavy-oil reservoirs | |
| Tremblay | Cold flow: a multi-well cold production (CHOPS) model | |
| CN106321076B (zh) | 一种注水井启动压力测试方法 | |
| CN115221666A (zh) | 一种排水采气工艺优选方法及系统 | |
| CN110714755B (zh) | 水驱油藏剩余油二次富集速度快速预测方法 | |
| CN108984952B (zh) | 一种油田注入能力变化预测方法 | |
| CN113971528A (zh) | 一种优势渗流通道的识别方法 | |
| CN112647930B (zh) | 一种水平井油藏工程找水方法 | |
| Zhang et al. | The optimal model of water control completion based on source function and network model | |
| CN114417687A (zh) | 一种软硬双层结构粘弹性颗粒在油藏中的调驱模拟方法 | |
| CN114439460A (zh) | 一种用于海上疏松砂岩油藏大孔道的试井分析方法 | |
| CN113177363A (zh) | 一种储层大孔道参数定量表征的油藏工程方法 | |
| CN111022007B (zh) | 一种新型油藏流场优势通道识别方法 | |
| Chen et al. | Modeling of counter-current spontaneous imbibition in independent capillaries with unequal diameters |
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 |