生产井开发指标的分层分方向定量计算方法

张继成，任帅*，林立，赵帅，冯诗淼

（1.东北石油大学 石油工程学院，黑龙江 大庆 163318；2.大庆油田有限责任公司第三采油厂，黑龙江 大庆 163000）

生产井开发指标的分层分方向定量计算方法

张继成1，任帅1*，林立2，赵帅1，冯诗淼2

（1.东北石油大学 石油工程学院，黑龙江 大庆 163318；2.大庆油田有限责任公司第三采油厂，黑龙江 大庆 163000）

针对注水开发的多层砂岩油藏分层动态分析难度大等问题，在常规井层开发指标计算基础上，结合动、静态劈分方法，综合考虑渗透率、孔隙度、地层系数、含水饱和度、位置系数、措施系数及注水量系数，提出了一种既可将油、水井作为统一整体，又可对小层、方向流动分量开发指标进行定量计算的体现渗流力学本质的方法。用大庆油田N2-O1井组的产液剖面资料进行验证。结果表明，所提方法的计算结果与测量结果吻合度较高，精度平均值达75.11%。用该方法计算开发指标，适用性强，能较真实地反映各小层、各方向的产液情况，对现场应用具有指导意义。

注采关系；开发效果；数学模型；分层注水；流动分量

为更有效地制订分层调整措施和挖潜策略，提高油田的开发能力，须充分掌握油水井分层动用状况，目前常见的计算井层开发指标的方法有地层系数法［1-4］、数值模拟法［5-7］、动态劈分方程法［8-10］、数学统计与回归公式法［11-13］、物理实验法［14-16］等。地层系数法，劈分方法简单，参数单一，但因未考虑注采连通情况，只适合射孔层数少、储层非均质性强的砂岩油藏；数值模拟法，需要大量开发动态数据进行拟合，适合开发时间长的油田；动态劈分方程法，是基于渗流力学进行劈分计算的一种方法，在计算过程中因涉及在实际油田开发中难以确定的多个参数，实用性不高；数学统计与回归公式法，通常针对单个油田开发过程中的动静态参数建立经验公式，对其他油田不一定适用，不利于推广应用；物理实验法，成本高，且只能计算某时刻的开发指标，无法满足实际需要。此外，以上方法只能计算至小层开发指标，无法定量计算油水井连通方向的开发指标。

本文在以往研究的基础上，综合考虑渗透率、孔隙度、地层系数、含水饱和度、位置系数、措施系数及注水量系数等，提出了一种可定量计算小层、方向流动分量开发指标的方法。

1 采油井流动分量系数

1.1 位置系数

先根据王治国等［17］提出的速率关联度分析法确定井组的连通情况，再通过注采井之间的液量关联度定量描述井间的连通情况，进而计算位置系数。某井某层某方向的位置系数

其中，n为周围连通水井数，θ为某方向方位角。

如图1所示，在以油井O为中心，水井A，B，C为方向的注采井连通关系中，求水井A方向的位置系数，需先计算∠1和∠2。

图1 角度计算示意图Fig.1 Schematic diagram of angle calculation

假设水井A坐标为（x1，y1），水井B坐标为（x2，y2），水井C坐标为（x3，y3），油井O坐标为（x0，y0），则AB连线长为

同理，可求得|BO|，|AO|，|CO|和|AC|，所以有

其中，M为AB的中点。已知|AM|，|AO|和∠1，可得

所以有

同理，可求得∠NOA，最终可得水井A方向的波及角为。

设中心油井连通的水井数为n，则水井A方向的位置系数为

1.2 平均地层系数

平均地层系数=（水井地层系数+油井地层系数）/2，计算式为

其中，Ki为水井i的有效渗透率（μm2）；Hi为水井i的砂岩厚度（m）；Kk为油井k的有效渗透率（μm2）；Hk为油井k的砂岩厚度（m）；为水井i与油井k连通方向的平均地层系数（μm2·m）。

1.3 注水量系数

中心油井在各小层的产油量与产水量均可劈分到各方向，而各方向的产油量和产水量又受连通方向水井注水量的影响。因此，在进行各方向采出量劈分计算时，需引入注水量系数：

其中，Wik为水井i到油井k在j小层劈分的方向注水量（m3）。

1.4 连通状况系数

在以油井为中心、水井为方向的连通关系中，连通状况系数是衡量连通方向水井与中心油井砂体连通关系的一个指标，反映储层的宏观非均质性［18］，值越大，均质性越强。渗透率变异系数通常用于表征储层渗透率的非均质性，值越小，均质性越强，反之，非均质性越强。所以，可用渗透率变异系数的倒数表示连通状况系数：

其中，Z为水井i与油井k某小层连通方向的连通状况系数；V为渗透率变异系数；Ki为第i口水井在该小层的渗透率（μm2）；hi为第i口水井在该小层的折算有效厚度（m）；为渗透率平均值（μm2）；l为连通方向小层渗透率变化量。

1.5 措施系数

根据该地区的现场情况，确定连通油井k的措施系数Mk，见表1。

表1 措施系数Table 1 Measure coefficient

1.6 井 距

其中，xi为水井i的横坐标（m）；xk为油井k的横坐标（m）；yi为水井i的纵坐标（m）；yk为油井k的纵坐标（m）；为对数化后的井距（m）。

1.7 相 渗

根据含水饱和度和大庆油田相渗关系曲线求相渗。计算相渗前，先计算不同时刻含油饱和度。

（1）含油饱和度。

根据油水井所在最优微相选择相渗曲线，在相渗曲线上，水相相对渗透率零点所对应的含水饱和度为束缚水饱和度，1减去束缚水饱和度得到该方向的原始含油饱和度。

（2）各方向含油饱和度。

各方向原始地质储量等于该方向所在区域体积与孔隙度和含油饱和度的乘积：

其中，N为方向域内的地质储量（t）；H为小层折算厚度（m）；A为某四边形区域面积（m2）；为加权平均后的孔隙度；为加权平均后的原始含油饱和度；为地面原油密度（g·cm-3）；Boi为原油体积系数。

对于刚开发的油田，其油相相对渗透率Kro=1，水相相对渗透率Krw=0，设水井与油井连通方向第1年的产油量为N，则第1年后，该方向剩余地质储量为N1=N0-N，所以该方向剩余含油饱和度为

若受效油井数未发生变化，在计算第m年某方向波及区域的流动分量条件值时，该方向所有的油、水相相对渗透率均可取第m-1年剩余含油饱和度为o（m-1）时所对应的油、水相相对渗透率，得到第m年流动分量、剩余地质储量和剩余含油饱和度。

若受效油井数发生变化，当出现新井F时，水井W的受效油井数增加，如图2所示。其中WF方向波及区域为四边形WrFs。四边形WrFs的出现，导致四边形WpEq和WqAm减小，即WE和WA方向波及区域减小。

（a）对比图（b）出现新井图

图2 出现新井对比示意图

Fig.2 New well comparison diagram

其中，三角形Wyq范围内所有网格剩余含油饱和度与WE方向波及区域WpEr的剩余含油饱和度相同；三角形Wzq范围内所有网格剩余含油饱和度与WA方向波及区域WsAm的剩余含油饱和度相同；五边形Fszyr为原先未开发区域，其剩余含油饱和度仍为原始含油饱和度，见图3。

图3 出现新井含油饱和度计算示意Fig.3 The new well saturation calculation diagram appears

（3）相渗归一化处理。

由于油、水井连通方向上存在不同微相，同一含油饱和度对应不同相渗，因此需对油水井连通方向不同微相的油、水相渗进行归一化处理，使相渗更贴近油藏地质真实情况。

1.8 流动分量系数

（1）单层单方向产水流动分量条件值：

（2）单层单方向产油流动分量条件值：

其中，Kro为油相相对渗透率。

（3）单层单方向产水流动分量系数。

将某小层各方向产水流动分量条件值相加，该小层某方向产水流动分量条件值与该小层产水流动分量条件值之和的比值就是该方向在该小层的产水流动分量系数：

其中，Cw为某方向产水流动分量系数；Fwik为某方向产水流动分量条件值。

（4）单层单方向产油流动分量系数。

将某小层各方向产油流动分量条件值相加，该小层某方向产油流动分量条件值与该小层产油流动分量条件值之和的比值就是该方向在该小层的产油流动分量系数：

其中，Co为某方向产油流动分量系数；Foik为某方向产油流动分量条件值。

（5）单层产水流动分量条件值。

将某小层各方向流动分量条件值相加，得该小层的产水流动分量条件值：

其中，Fwj为某小层产水流动分量条件值。

（6）单层产油流动分量条件值。

将某小层各方向产油流动分量条件值相加，得该小层的产油流动分量条件值：

其中，Foj为某小层产油流动分量条件值。

（7）单层产水流动分量系数。

将某井各小层产水流动分量条件值相加，某井某小层产水流动分量条件值与该井产水流动分量条件值总和的比值就是该小层在该井的产水流动分量系数：

其中，Swj为单层产水流动分量系数。

（8）单层产油流动分量系数。

将某井各小层产油流动分量条件值相加，某井某小层产油流动分量条件值与该井产油流动分量条件值总和的比值就是该小层在该井的产油流动分量系数：

其中，Soj为单层产油流动分量系数。

2 注水井流动分量系数

2.1 流动分量条件值

（1）单层单方向流动分量条件值：

（2）单层流动分量条件值：

其中，Yj为某小层流动分量条件值。

2.2 流动分量系数

（1）单层单方向流动分量系数：

其中，Cik为某方向流动分量系数。

（2）单层流动分量系数：

其中，Sj为单层流动分量系数。

3 井层方向指标

3.1 注采比

其中，Rip为方向注采比，Wik为方向月注水量（m3），Qwik为方向月产水量（m3），Qoik为方向月产油量（m3），Bo为原油体积系数。

3.2 含水率

其中，fw为含水率。

4 应用实例

4.1 井组概况

N2-O1油井于1998年7月投产，共射开19个小层，主要沉积类型为河道砂、主体砂、非主体砂和表外砂，小层平均孔隙度15.1%，渗透率0.112～0.390 μm2，平均渗透率0.161 μm2，截至2021年9月，该井累计产油量14 443 m3，累计产液量221 807 m3，月产液量为955 m3，含水率为95.49%。周围连通的水井有N2-W1、N2-W2、N2-W3、N2-W4和N2-W5。N2-O1井组注采对应关系如图4所示。

图4 N2-O1井组注采关系示意Fig.4 Schematic of injection production relationship of well cluster N2-O1

4.2 注采指标

选取N2-O1井组的SIII3b小层作为研究对象。在SIII3b小层，N2-O1油井连通3个方向水井，注采关系如图5所示。

图5 N2-O1井组SIII3b小层注采关系示意Fig.5 Schematic of injection production relationship of SIII3b small layer of N2-O1 well cluster

整理N2-O1油井连通的N2-W4、N2-W3和N2-W2水井的物性参数，并计算SIII3b小层方向的流动分量系数，结果见表2～表4。

表2 N2-O1井组SIII3b小层物性参数Table 2 Physical property parameters of SIII3b sublayer of N2-O1 well cluster

表3 N2-O1油井SIII3b小层受效井各参数计算结果Table 3 Calculation results of parameters of affected wells in SIII3b sublayer of N2-O1 oil well

表4 N2-O1油井SIII3b小层受效井流动分量条件值计算结果Table 4 Calculation results of flow component condition parameters of effective well in SIII3b small layer of N2-O1 oil well

表5 N2-O1油井各小层开发指标计算结果Table 5 Calculation results of development indexes of each small layer of N2-O1 production well

当得到条件值中的各参数后，可计算小层各方向流动分量条件值，进而得到各小层流动分量条件值，最终得到各小层流动分量系数。将N2-O1井组各井资料代入相应计算式，与各井各小层流动分量系数相乘，得到各井各小层月产水量、月产油量、含水率等开发指标，见表5。

4.3 开发状况评价

4.3.1 结果验证

整理N2-O1油井各层段产液剖面测量数据，并与各层段流动分量计算结果进行对比，结果见图6和表6。

图6 产液剖面测量结果与计算结果对比Fig.6 Comparison between test and calculation results of liquid production profile of N2-O1 oil well

表6 N2-O1油井产液剖面测量值与计算结果对比Table 6 Comparison between test and calculation results of liquid production profile of N2-O1 oil well

可知，各层段产液占比测量结果与计算结果较为接近。层段1（SII5、SII6、SII13）与层段6（PII1a、PII9）的计算精度相对较低，其原因可能为实际测量时取的是单次测量结果，具有偶然性。整体看，各层段的产液占比计算值与测量值变化趋势基本一致，吻合度较高，平均精度达75.11%。

4.3.2 开发现状评价

由N2-O1油井各小层流动分量的计算结果和产液剖面实际产液量知，第4个层段SIII4b小层产液量最多，占全井产液量的30%左右，第2个层段SII15～SII16小层产液量最少，占全井产液量的3%左右。SII5小层的含水率最低，为71.28%，其次是SII15小层，含水率为78.58%。综合分析得到，N2-O1油井的SII5、SII15和SIII6小层潜力较大，应加强注水开发。

5 结论

5.1在以往只考虑小层开发指标计算的基础上，提出了方向开发指标计算方法，除考虑渗透率、厚度、位置系数、井距、含油饱和度、措施系数外，还加入了注水量系数，将油、水井作为整体考虑。

5.2用N2-O1井组的产液剖面现场资料验证本文方法，各层段产液占比计算结果与测量结果吻合度较高，二者变化趋势基本一致，平均精度达75.11%。对N2-O1井组进行综合评价，认为该井组的SII5、SII15和SIII6小层潜力较大，应加强注水开发。

［1］李源流，杨兆平，潘多寿，等. 考虑相对渗透率变化的砂岩油藏产量劈分方法［J］. 西安石油大学学报（自然科学版），2020，35（3）：72-76， 85. DOI：10.3969/j.issn.1673-064X.2020.03.010

LI Y L， YANG Z P，PAN D S，et al. A Production splitting method of sandstone reservoir considering the change of relative permeability［J］. Journal of Xi′an Shiyou University （Natural Science Edition），2020，35（3）：72-76， 85. DOI：10.3969/j.issn.1673-064X. 2020.03.010

［2］夏静，冯阳，陈新宇. 油井分层注水倍数计算方法［J］. 数学的实践与认识，2017，47（3）：114-119.

XIA J， FENG Y，CHEN X Y. Computing method of separating layer water injected volume about oil wells［J］. Mathematics in Practice and Theory， 2017，47（3）：114-119.

［3］邸井东. 结合含水饱和度场分布的水淹级别评价体系［J］. 石油化工高等学校学报，2019，32（3）：51-57. DOI：10.3969/j.issn.1006-396X.2019.03.009

DI J D. Evaluation system of water flooded grade combined with distribution of water saturation field［J］. Journal of Petrochemical Universities，2019，32（3）：51-57. DOI：10.3969/j.issn.1006-396X.2019.03.009

［4］张继成，王潇悦. 考虑含水饱和度的产量劈分方法及应用［J］. 浙江大学学报（理学版），2015，42（5）：626-630. DOI：10.3785/j.issn.1008-9497.2015.05.020

ZHANG J C， WANG X Y. The method of oil production splitting based on water saturation and its application［J］. Journal of Zhejiang University （Science Edition），2015，42（5）：626-630. DOI：10. 3785/j.issn.1008-9497.2015.05.020

［5］陈建华，晏庆辉，骆逸婷，等. 基于历史生产数据的多层合采井产量劈分新方法［J］.中国海上油气，2022，34（1）：110-116. DOI：10.11935/j.issn.1673-1506.2022.01.013

CHEN J H， YAN Q H，LUO Y T，et al. A historical production data based method for production splitting of multi-layer commingled gas wells［J］. China Offshore Oil and Gas， 2022， 34（1）： 110-116. DOI：10.11935/j.issn.1673-1506.2022.01.013

［6］LIU X，WANG P，LIANG X，et al. Research progress and surface/interfacial regulation methods for electrophotocatalytic hydrogen production from water splitting［J］. Materials Today Energy， 2020，18：100524. DOI：10.1016/j.mtener.2020.100524

［7］林孟雄，成育红，张林，等. 苏里格气田苏东区块多层产量劈分新方法［J］. 新疆地质，2019，37（3）：419-421. DOI：10.3969/j.issn.1000-8845.2019.03.024

LIN M X， CHENG Y H，ZHANG L，et al. A new method of production dividing for commingling production in the Sudong area［J］. Xinjiang Geology， 2019，37（3）：419-421. DOI：10.3969/j.issn.1000-8845.2019.03.024

［8］HEBA H F，SAYED Y S，NAGEH K A，et al. Emerging nanoporous anodized stainless steel for hydrogen production from solar water splitting［J］. Journal of Cleaner Production， 2020，274：122826. DOI：10.1016/j.jclepro.2020.122826

［9］周立国. 低渗砂岩油藏优势注采方向量化表征技术及应用［J］. 特种油气藏，2021，28（6）：98-104. DOI：10.3969/j.issn.1006-6535.2021.06.013

ZHOU L G. Quantitative characterization technology and application of dominant injection-production direction in low-permeability sandstone reservoirs［J］. Special Oil and Gas Reservoirs， 2021，28（6）：98-104. DOI：10.3969/j.issn.1006-6535.2021.06.013

［10］RODOLPHE D，ROBERT L P，ALANA G，et al. Reinventing residual reserves in the sea： Are we favouring ease of establishment over need for protection？［J］. Aquatic Conservation： Marine and Freshwater Ecosystems，2015，25（4）：480-504. DOI：10.1002/aqc.2445

［11］杨兆平，岳世俊，郑长龙，等. 薄互层砂岩油藏多因素综合约束的产量劈分方法［J］. 岩性油气藏，2018，30（6）：117-124. DOI：10.12108/yxyqc.20180614

YANG Z P， YUE S J，ZHENG C L，et al. Production split method restricted synthetically by multi-factors in thin interbed sandstone reservoirs［J］. Lithologic Reservoirs， 2018，30（6）：117-124. DOI：10.12108/yxyqc.20180614

［12］王立，喻高明，傅宣豪，等. 基于反距离加权插值法的产量劈分新方法［J］. 断块油气田，2018，25（5）：617-621. DOI：10.6056/dkyqt201805015

WANG L， YU G M，FU X H，et al. New method for production cleavage by inverse distance weighted interpolation［J］. Fault-block Oil and Gas Field， 2018，25（5）：617-621. DOI：10.6056/dkyqt201805015

［13］张继成，单五一，卢光夫，等. 一种预测注采关系调整效果的数学模型［J］. 数学的实践与认识，2020，50（4）：162-170.

ZHANG J C， SHAN W Y，LU G F，et al. A mathematical model to predict the effects of adjustment between injection and production［J］. Mathematics in Practice And Theory， 2020，50（4）：162-170.

［14］赵辉，刘邓，江厚顺，等. 基于数据物理驱动的油藏调剖动态预测新方法［J］. 长江大学学报（自然科学版），2021，18（1）：60-69. DOI：10.3969/j.issn.1673-1409.2021.01.007

ZHAO H， LIU D，JIANG H S，et al. A new method of reservoir profile control dynamic prediction based on the data-driven model［J］. Journal of Yangtze University （Natural Science Edition），2021，18（1）：60-69. DOI：10.3969/j.issn.1673-1409.2021.01.007

［15］HARIMURTI D. Indonesian milestone in production-sharing contract in perspective of government take， contractor take， cost recovery and production target［C］// Spe/Iatmi Asia Pacific Oil amp; Gas Conference amp; Exhibition. Jakarta：Silverchair，2017. DOI：10.2118/187008-MS

［16］谢书剑，周玉辉，丁艳雪，等. 水驱油藏分层注水效率评价新方法［J］. 中国科技论文，2021，16（5）：475-481. DOI：10.3969/j.issn.2095-2783.2021.05.004

XIE S J， ZHOU Y H，DING Y X，et al. A new method to evaluate the injection efficiency of different layers in waterflooding reservoir［J］. China Sciencepaper， 2021，16（5）：475-481. DOI：10. 3969/j.issn.2095-2783.2021.05.004

［17］王治国，张继成，宋考平，等. 应用速率关联度分析井间连通情况［J］. 特种油气藏，2011，18（3）：91-93， 140. DOI：10.3969/j.issn.1006-6535.2011.03.024

WANG Z G， ZHANG J C，SONG K P，et al. Analyze interwell communication by rate correlation［J］. Special Oil and Gas Reservoirs， 2011，18（3）：91-93， 140. DOI：10.3969/j.issn.1006-6535.2011.03.024

［18］娄瑞，刘宗堡，张云峰，等. 齐家地区高Ⅲ、高Ⅳ组砂岩输导层静态连通评价［J］. 断块油气田，2019，26（3）：285-289. DOI：10.6056/dkyqt201903003

LOU R， LIU Z B，ZHANG Y F，et al. Static connectivity evaluation of sandstone carrier bed in GⅢ-GⅣ groups in Qijia area［J］. Fault Block Oil and Gas Field， 2019，26 （3）：285-289. DOI：10.6056/dkyqt201903003

Quantitative calculation method of development indexes for layered and directional of production wells

ZHANG Jicheng1, REN Shuai1, LIN Li2, ZHAO Shuai1, FENG Shimiao2

（1.School of Petroleum Engineering，Northeast Petroleum University，Daqing163318，Heilongjiang Province，China；2. No.3Oil Production Plant，Daqing Oilfield Co.，Ltd.，Daqing163000，Heilongjiang Province，China）

On account of the difficulty of Layered Dynamic Analysis of multi-layer sandstone reservoir developed by water flooding, based on the calculation approach of conventional well layer development indexes, combined with dynamic and static splitting methods, and comprehensively considering the parameters such as permeability, porosity, formation coefficient, water saturation, position coefficient, measure coefficient and water injection coefficient, this paper propose a new method which can not only deal with oil and water wells as a unified whole, but also for small layers. The mathematical method for quantitative calculation of development index of directional flow component reflects the essence of seepage mechanics. The proposed method is verified by the liquid production profile data of N2-O1 well group in Daqing Oilfield, and shows high coincidence with the actual production, with average accuracy 75.11%. It has strong applicability to the calculation of development indicators and plays an important guiding role in field application.

injection-production relationship; development effect; mathematical model; layered water injection; flow component

P 548

A

1008⁃9497（2022）05⁃598⁃08

2021⁃11⁃03.

黑龙江省自然科学基金资助项目（E201407）.

张继成（1972—），ORCID：https：//orcid.org/0000-0002-6791-7767，男，博士，教授，主要从事石油与天然气工程研究，E-mail：zhangjc777@163.com.

通信作者，ORCID：https：//orcid.org/ 0000-0002-5274-1674，E-mail：769277251@qq.com.

10.3785/j.issn.1008-9497.2022.05.011