CN113065705B - 海上薄互层砂岩油藏多层合采加密井可采储量预测方法 - Google Patents
海上薄互层砂岩油藏多层合采加密井可采储量预测方法 Download PDFInfo
- Publication number
- CN113065705B CN113065705B CN202110371858.9A CN202110371858A CN113065705B CN 113065705 B CN113065705 B CN 113065705B CN 202110371858 A CN202110371858 A CN 202110371858A CN 113065705 B CN113065705 B CN 113065705B
- Authority
- CN
- China
- Prior art keywords
- well
- production
- perforation
- encryption
- layer
- 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
Links
- 238000004519 manufacturing process Methods 0.000 title claims abstract description 225
- 238000000034 method Methods 0.000 title claims abstract description 49
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 182
- 238000002347 injection Methods 0.000 claims abstract description 181
- 239000007924 injection Substances 0.000 claims abstract description 181
- 238000004891 communication Methods 0.000 claims abstract description 63
- 239000011148 porous material Substances 0.000 claims abstract description 13
- 238000012360 testing method Methods 0.000 claims description 56
- 239000003921 oil Substances 0.000 claims description 26
- 230000001186 cumulative effect Effects 0.000 claims description 20
- 238000004458 analytical method Methods 0.000 claims description 13
- 230000015572 biosynthetic process Effects 0.000 claims description 10
- 239000012530 fluid Substances 0.000 claims description 8
- 230000003068 static effect Effects 0.000 claims description 6
- 238000010521 absorption reaction Methods 0.000 claims description 5
- 239000007788 liquid Substances 0.000 claims description 4
- 238000009530 blood pressure measurement Methods 0.000 claims description 3
- 238000013213 extrapolation Methods 0.000 claims description 3
- 239000008398 formation water Substances 0.000 claims description 3
- 238000013507 mapping Methods 0.000 claims description 3
- 230000035699 permeability Effects 0.000 claims description 3
- 239000010779 crude oil Substances 0.000 claims description 2
- 238000013461 design Methods 0.000 claims description 2
- 238000004364 calculation method Methods 0.000 abstract description 2
- 238000011161 development Methods 0.000 description 4
- 238000004088 simulation Methods 0.000 description 4
- 239000000243 solution Substances 0.000 description 4
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 238000005553 drilling Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000003345 natural gas Substances 0.000 description 1
- 239000002332 oil field water Substances 0.000 description 1
- 239000003209 petroleum derivative Substances 0.000 description 1
- 238000011084 recovery Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
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
- G06Q10/00—Administration; Management
- G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
-
- 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
-
- 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)
- Economics (AREA)
- Strategic Management (AREA)
- Physics & Mathematics (AREA)
- Tourism & Hospitality (AREA)
- Theoretical Computer Science (AREA)
- Marketing (AREA)
- General Physics & Mathematics (AREA)
- General Business, Economics & Management (AREA)
- Quality & Reliability (AREA)
- Animal Husbandry (AREA)
- Game Theory and Decision Science (AREA)
- Operations Research (AREA)
- Entrepreneurship & Innovation (AREA)
- Life Sciences & Earth Sciences (AREA)
- Agronomy & Crop Science (AREA)
- Development Economics (AREA)
- Marine Sciences & Fisheries (AREA)
- Mining & Mineral Resources (AREA)
- Health & Medical Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明涉及一种海上薄互层砂岩油藏多层合采加密井可采储量预测方法,包括如下步骤:1)统计动态注采连通系数,并计算生产井与注水井动态连通的射孔小层总流动系数;2)外推得到含水率达到98%时的已投产生产井单井可采储量;3)根据初期含水率分类建立不同初期含水率下单井可采储量与动态连通射孔小层总流动系数的关系图版;4)计算加密井区内各射孔小层原始孔隙体积;5)计算注水井在加密井区内各射孔小层总注水量;6)计算加密井区射孔层的注入PV数;7)利用油田注入PV数与油田含水率的关系图版插值计算加密井初期含水率;8)计算加密井与注水井动态连通的射孔小层总流动系数;9)计算加密井可采储量。
Description
技术领域
本发明涉及一种海上薄互层砂岩油藏多层合采加密井可采储量预测方法,属于石油天然气开采技术领域。
背景技术
对于处于开发中后期的海上薄互层砂岩油藏,由于纵向各层注采连通性差异、层间干扰以及投产初期见水等因素影响,其加密井可采储量预测难度较大。目前主要采用油藏数值模拟的方法对加密井的可采储量进行预测。然而对于薄互层砂岩油藏而言,即使花费大量时间和精力完成了长生产历史的拟合,其油藏数值模拟模型仍然有很大的不确定性。因此有必要在已投产生产井开发现状类比的基础上,通过油藏工程原理建立一套海上薄互层砂岩油藏多层合采加密井可采储量预测方法,用于指导加密井井位决策和油藏数值模拟研究。
发明内容
针对上述突出问题,本发明提供一种海上薄互层砂岩油藏多层合采加密井可采储量预测方法,该方法适用于采用定向井合采为主要开发方式、注水保压且处于开发中后期的海上薄互层砂岩油藏。
为实现上述目的,本发明采取以下技术方案:
一种海上薄互层砂岩油藏多层合采加密井可采储量预测方法,包括如下步骤:
1)根据测井水淹定量解释结果、MDT测压数据、PLT测试数据以及地质专业静态注采连通分析结果,分析生产井和注水井在各射孔小层的动态连通性,统计动态注采连通系数A,并计算生产井与注水井动态连通射孔小层总流动系数Msum;
2)按月统计已投产生产井单井的累积产油量Np和累积产水量Wp,对累积产油量与累积产水量的拟合关系式Np=NR-aNp 2/Wp q进行系数识别,并外推Np得到含水率达到98%时的已投产生产井单井可采储量Nplast,其中NR、a、q为拟合系数;
3)计算单井生产前3个月平均含水率作为初期含水率,根据单井的初期含水率分类建立不同初期含水率下所述步骤2)中的单井可采储量Nplast与所述步骤1)中动态连通射孔小层总流动系数Msum的关系图版;
4)根据加密井设计井位、射孔层位、地质油藏特征以及当前井网特征,从整个油藏中划分加密井区,计算加密井区内各射孔小层原始孔隙体积(Vpwell)n;
5)根据所述步骤1)中的动态注采连通分析结果,将与加密井区内的生产井存在动态连通关系的所有注水井的注水量劈分至各射孔小层内,计算注水井在加密井区内各射孔小层总注水量INJn;
6)根据所述步骤4)中的加密井区各射孔小层原始孔隙体积(Vpwell)n和所述步骤5)中注水井在加密井区内各射孔小层总注水量INJn,计算加密井区射孔层的注入PV数(PV数指孔隙体积倍数,即累积注入量与地下孔隙体积之比)PVall;
7)绘制油田注入PV数PVf与油田含水率fwf的关系图版,依据所述步骤6)中的加密井区射孔层的注入PV数PVall,利用图版插值计算加密井初期含水率Fw;
8)根据所述步骤1)计算加密井与注水井动态连通射孔小层总流动系数Msum_infill;
9)根据所述步骤7)中的加密井初期含水率Fw和所述步骤8)中加密井与注水井动态连通射孔小层总流动系数Msum_infill,利用所述步骤3)不同初期含水率下单井可采储量Nplast与动态连通射孔小层总流动系数Msum的关系图版计算加密井可采储量。
所述的可采储量预测方法,优选地,所述步骤1)中生产井与注水井动态连通射孔小层总流动系数Msum为:
式中,Abn为生产井在射孔小层n的动态注采连通系数,取值为0或1;G为生产井的射孔总层数;g为生产井射孔层编号,取值为1,2,3,…,G;kp为生产井射孔层的渗透率;hp为生产井射孔层的有效厚度;μo为生产井射孔层内的原油粘度。
所述的可采储量预测方法,优选地,所述步骤2)的具体外推方法如下:
①以生产井最后两个月累积产油量作为Np(0)和Np(1),等间距外推Np(2),Np(3),…,Np(j);
②根据所述步骤2)中拟合得到的系数NR、a、q,计算相应的累积产水量Wp(j):
式中,Np(j)为步骤①中外推生产井累积产油量;Wp(j)是与Np(j)对应的生产井累积产水量;NR、a、q为拟合系数;
③判断含水率达到98%时的单井可采储量Nplast,当Np(j)满足如下条件时:
此时Nplast=Np(j);
式中,Np(j)、Np(j-1)均是步骤①中外推生产井累积产油量;Wp(j)、Wp(j-1)分别是与Np(j)、Np(j-1)对应的生产井累积产水量;其中,j为步骤①中外推累积产油量和对应累积产水量的编号。
所述的可采储量预测方法,优选地,所述步骤3)中根据单井的初期含水率分类建立单井可采储量Nplast与动态连通射孔小层总流动系数Msum的关系图版,其分类建立原则为将初期含水率相同的生产井划分为同一类建立图版,其中初期含水率按四舍五入原则保留小数点后一位。
所述的可采储量预测方法,优选地,所述步骤4)中加密井区内各射孔小层原始孔隙体积(Vpwell)n为:
(Vpwell)n=(Sareahφ)n
式中,Sarea为加密井区内射孔小层n的面积;h为加密井区内射孔小层n的有效厚度;φ为加密井区内射孔小层n的地下孔隙度。
所述的可采储量预测方法,优选地,所述步骤5)中注水井在加密井区内各射孔小层总注水量INJn为:
式中,INJn为所有与加密井区内生产井动态连通的注水井在射孔小层n加密井区范围内的注水量;M为与加密井区内生产井动态连通的注水井数量;m为与加密井区内生产井动态连通的注水井编号,取值为1,2,3,…,M;(INJn)m为注水井m在射孔小层n加密井区范围内的注水量;(INJtotal_m)i为注水井m在第i次PLT测试和第i+1次PLT测试期间的累积注水量;(fPLT_m)i为注水井m第i次PLT测试中包含射孔小层n的测试段的吸水比例;NT为注水井PLT测试次数;i为PLT测试次数编号,取值为1,2,3,…,NT;(kwhw)n为注水井m在射孔小层n的地层系数;NN为包含射孔小层n的注水井PLT测试段内所有射孔小层总数;nn为PLT测试段内射孔小层编号,取值为1,2,3,…,NN;(kwhw)nn为注水井m在射孔小层nn的地层系数;[(Wpo)n_i]np为注水井m在第i次PLT测试和第i+1次PLT测试期间生产井np在射孔小层n的累积产液量;NP为加密井区内生产井数量;np为加密井区内生产井编号,取值为1,2,3,…,NP;D为与注水井m存在动态连通关系的生产井总井数;d为与注水井m存在动态连通关系的生产井编号,取值为1,2,3,…,D;[(Wpo)n_i]d为注水井m在第i次PLT测试和第i+1次PLT测试期间生产井d在射孔小层n的累积产液量;[Qtotal_i]np为生产井np在注水井m第i次PLT测试和第i+1次PLT测试期间的累积产液量;(knphnp)n为生产井np在射孔小层n的地层系数;Anp_n为生产井np在射孔小层n与加密井区内注水井的动态注采连通系数,取值为0或1,由所述步骤1)得到;G是生产井射孔层总数;g是生产井射孔层编号,取值为1,2,3,…,G;Anp_g是生产井np在射孔小层g与加密井区内注水井的动态注采连通系数,取值为0或1,由所述步骤1)得到;(knphnp)g为生产井np在射孔小层g的地层系数;[Qtotal_i]d是生产井d在注水井m第i次PLT测试和第i+1次PLT测试期间的累积产液量;(kdhd)n是生产井d在射孔小层n的地层系数;Ad_n是生产井d在射孔小层n与加密井区内注水井的动态注采连通系数,取值为0或1,由所述步骤1)得到;Ad_g是生产井d在射孔小层g与加密井区内注水井的动态注采连通系数,取值为0或1,由所述步骤1)得到;(kdhd)g是生产井d在射孔小层g的地层系数。
所述的可采储量预测方法,优选地,所述步骤6)中加密井区射孔层的注入PV数PVall为:
式中,NL为加密井射孔小层总数;n为加密井射孔小层编号,取值为1,2,3,…,NL;Bw为加密井区内射孔小层n的地层水体积系数。
本发明由于采取以上技术方案,其具有以下优点:
1、本发明提供了一种基于油藏工程原理和油田生产历史类比的海上薄互层砂岩油藏多层合采加密井可采储量预测方法,预测结果可靠。
2、本发明方法避免了油藏数值模拟方法成本高、周期长的问题,可以快速实现多层合采加密井的可采储量估算,为薄互层砂岩油藏的加密井井位决策和历史拟合研究提供指导。
3、本发明方法可以广泛应用于处于注水开发中后期的同类型薄互层砂岩油藏。
附图说明
图1为本发明一实施例提供的可采储量预测方法流程图;
图2为本发明该实施例提供的依据单井的初期含水率分类建立的不同初期含水率下单井可采储量和动态注采连通射孔小层的总流动系数关系图版。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚,下面对本发明中的技术方案进行清楚、完整地描述,显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
如图1所示,本发明提供一种海上薄互层砂岩油藏多层合采加密井的可采储量预测方法,包括如下步骤:
1)根据测井水淹定量解释结果、MDT测压数据、PLT测试数据以及地质专业静态注采连通分析结果,分析生产井和注水井在各射孔小层的动态连通性,统计动态注采连通系数A,并计算生产井与注水井动态连通的射孔小层总流动系数Msum;具体分析和计算过程如下:
①选择生产井B为分析井,根据B井射孔层位逐射孔小层分析B井的动态注采连通性;
②根据B井或其邻井钻后测井水淹定量解释结果分析动态注采连通性:若B井钻后测井解释在射孔小层n中已经水淹,且B井位于老井注采流线附近,则B井在射孔小层n中与注水井动态连通,动态注采连通系数Abn=1;若B井在射孔小层n中测井解释未水淹,但比B井投产晚的邻井在射孔小层n中水淹,且B井位于老井注采流线附近,则B井在射孔小层n中与注水井动态连通,动态注采连通系数Abn=1;若B井及其周围邻井均未水淹,则仅利用水淹解释结果无法判断动态注采连通性;
③根据B井的MDT(Modular Dynamics Tester,模块式地层动态测试器)测压资料分析动态注采连通性:对于B井及邻井均未水淹的情况,若MDT显示B井在射孔小层n中超压(目前地层压力高于原始地层压力)且地质分析结果为静态连通,则B井在射孔小层n中动态注采连通,动态注采连通系数Abn=1;若MDT测试显示B井在射孔小层n中欠压(目前地层压力低于原始地层压力)或无MDT测试资料,且地质分析结果为静态不连通或连通性差,则B井在射孔小层n动态不连通,动态注采连通系数Abn=0;若MDT测试显示B井在射孔小层n中欠压或无MDT测试资料,而地质分析结果为静态注采连通性好且对应注水井PLT(ProductionLogging Test,生产测井)测试显示射孔小层n明显吸水(吸水比例大于5%),则B井在射孔小层n动态注采连通,动态注采连通系数Abn=1;
④根据动态连通性分析结果,计算生产井B与注水井动态连通射孔小层总流动系数Msum:
式中,Abn为生产井在射孔小层n的动态注采连通系数,取值为0或1,由步骤②和③分析得到;G为生产井的射孔总层数,单位是层;g为生产井射孔层编号,无单位,取值为1,2,3,…,G;kp为生产井射孔层的渗透率,单位是mD;hp为生产井射孔层的有效厚度,单位是m;μo为生产井射孔层内的原油粘度,单位是mPa.s;
⑤对所有生产井按照①-④进行分析和计算。
2)按月统计已投产生产井单井的累积产油量Np和累积产水量Wp,对累积产油量与累积产水量的拟合关系式Np=NR-aNp 2/Wp q进行系数识别,并外推Np得到含水率达到98%时的已投产生产井单井可采储量Nplast,其中NR、a、q为拟合系数;具体外推方法如下:
①以生产井最后两个月累积产油量作为Np(0)和Np(1),等间距外推Np(2),Np(3),…,Np(j);
②根据所述步骤2)中拟合得到的系数NR、a、q,计算相应的累积产水量Wp(j):
式中,Np(j)为步骤①中外推生产井累积产油量,单位是m3;Wp(j)是与Np(j)对应的生产井累积产水量,单位是m3;NR、a、q为拟合系数;
③判断含水率达到98%时的单井可采储量Nplast,当Np(j)满足如下条件时:
此时Nplast=Np(j);
式中,Np(j)、Np(j-1)均是步骤①中外推生产井累积产油量,单位是m3;Wp(j)、Wp(j-1)分别是与Np(j)、Np(j-1)对应的生产井累积产水量,单位是m3;其中,j为步骤①中外推累积产油量和对应累积产水量的编号,无单位。
3)计算单井生产前3个月平均含水率作为初期含水率,根据单井的初期含水率分类建立不同初期含水率下所述步骤2)中单井可采储量Nplast与所述步骤1)中动态连通射孔小层总流动系数Msum的关系图版;(分类建立原则为将初期含水率相同的生产井划分为同一类建立图版,其中初期含水率按四舍五入原则保留小数点后一位。)
4)根据加密井(加密井是指为提高油藏的产量和采收率,在已有的开发井之间所钻的井)设计井位、射孔层位、地质油藏特征以及当前井网特征,从整个油藏中划分加密井区,计算加密井区内各射孔小层原始孔隙体积(Vpwell)n:
(Vpwell)n=(Sareahφ)n
式中,Sarea为加密井区内射孔小层n的面积,单位是m2;h为加密井区内射孔小层n的有效厚度,单位是m;φ为加密井区内射孔小层n的地下孔隙度,无量纲。其中加密井区划分采用本领域技术人员惯用方法,故不再赘述。
5)根据所述步骤1)中的动态注采连通分析结果,将与加密井区内的生产井存在动态连通关系的所有注水井的注水量劈分至各射孔小层内,计算注水井在加密井区内各射孔小层总注水量INJn;
式中,INJn为所有与加密井区内生产井动态连通的注水井在射孔小层n加密井区范围内的注水量,单位是m3;M为与加密井区内生产井动态连通的注水井数量,单位是口;m为与加密井区内生产井动态连通的注水井编号,无单位,取值为1,2,3,…,M;(INJn)m为注水井m在射孔小层n加密井区范围内的注水量,单位是m3;(INJtotal_m)i为注水井m在第i次PLT测试和第i+1次PLT测试期间的累积注水量,单位是m3;(fPLT_m)i为注水井m第i次PLT测试中包含射孔小层n的测试段的吸水比例,无量纲;NT为注水井PLT测试次数,单位是次;i为PLT测试次数编号,无单位,取值为1,2,3,…,NT;(kwhw)n为注水井m在射孔小层n的地层系数,单位是mD.m;NN为包含射孔小层n的注水井PLT测试段内所有射孔小层总数,单位是层;nn为PLT测试段内射孔小层编号,无单位,取值为1,2,3,…,NN;(kwhw)nn为注水井m在射孔小层nn的地层系数,单位是mD.m;[(Wpo)n_i]np为注水井m在第i次PLT测试和第i+1次PLT测试期间生产井np在射孔小层n的累积产液量,单位是m3;NP为加密井区内生产井数量,单位是口;np为加密井区内生产井编号,无单位,取值为1,2,3,…,NP;D为与注水井m存在动态连通关系的生产井总井数,单位是口;d为与注水井m存在动态连通关系的生产井编号,无单位,取值为1,2,3,…,D;[(Wpo)n_i]d为注水井m在第i次PLT测试和第i+1次PLT测试期间生产井d在射孔小层n的累积产液量,单位是m3;[Qtotal_i]np为生产井np在注水井m第i次PLT测试和第i+1次PLT测试期间的累积产液量,单位是m3;(knphnp)n为生产井np在射孔小层n的地层系数,单位是mD.m;Anp_n为生产井np在射孔小层n与加密井区内注水井的动态注采连通系数,取值为0或1,由所述步骤1)得到;G是生产井射孔层总数,单位是层;g是生产井射孔层编号,无单位,取值为1,2,3,…,G;Anp_g是生产井np在射孔小层g与加密井区内注水井的动态注采连通系数,取值为0或1,由所述步骤1)得到;(knphnp)g为生产井np在射孔小层g的地层系数,单位是mD.m;[Qtotal_i]d是生产井d在注水井m第i次PLT测试和第i+1次PLT测试期间的累积产液量,单位是m3;(kdhd)n是生产井d在射孔小层n的地层系数,单位是mD.m;Ad_n是生产井d在射孔小层n与加密井区内注水井的动态注采连通系数,取值为0或1,由所述步骤1)得到;Ad_g是生产井d在射孔小层g与加密井区内注水井的动态注采连通系数,取值为0或1,由所述步骤1)得到;(kdhd)g是生产井d在射孔小层g的地层系数,单位是mD.m。
6)根据所述步骤4)中的加密井区各射孔小层原始孔隙体积(Vpwell)n和所述步骤5)中注水井在加密井区内各射孔小层总注水量INJn,计算加密井区射孔层的注入PV数(PV数指孔隙体积倍数,即累积注入量与地下孔隙体积之比)PVall:
式中,NL为加密井射孔小层总数;n为加密井射孔小层编号,取值为1,2,3,…,NL;Bw为加密井区内射孔小层n的地层水体积系数,单位是m3/m3。
7)绘制油田注入PV数PVf与油田含水率fwf的关系图版,依据所述步骤6)中的加密井区射孔层的注入PV数PVall,利用图版插值计算加密井初期含水率Fw;其中油田注入PV数PVf与油田含水率fwf关系图版采用本领域技术人员惯用方法,故不再赘述。
8)根据所述步骤1)计算加密井与注水井动态连通射孔小层总流动系数Msum_infill。
9)根据所述步骤7)中的加密井初期含水率Fw和所述步骤8)中加密井与注水井动态连通射孔小层总流动系数Msum_infill,利用所述步骤3)不同初期含水率下单井可采储量Nplast与动态连通射孔小层总流动系数Msum的关系图版计算加密井可采储量。
最后应说明的是:以上实施例仅用以说明本发明的技术方案,而非对其限制;尽管参照前述实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。
Claims (6)
1.一种海上薄互层砂岩油藏多层合采加密井可采储量预测方法,其特征在于,包括如下步骤:
1)根据测井水淹定量解释结果、MDT测压数据、PLT测试数据以及地质专业静态注采连通分析结果,分析生产井和注水井在各射孔小层的动态连通性,统计动态注采连通系数A,并计算生产井与注水井动态连通射孔小层总流动系数Msum;
2)按月统计已投产生产井单井的累积产油量Np和累积产水量Wp,对累积产油量与累积产水量的拟合关系式Np=NR-aNp 2/Wp q进行系数识别,并外推Np得到含水率达到98%时的已投产生产井单井可采储量Nplast,其中NR、a、q为拟合系数;
3)计算单井生产前3个月平均含水率作为初期含水率,根据单井的初期含水率分类建立不同初期含水率下所述步骤2)中的单井可采储量Nplast与所述步骤1)中动态连通射孔小层总流动系数Msum的关系图版;
4)根据加密井设计井位、射孔层位、地质油藏特征以及当前井网特征,从整个油藏中划分加密井区,计算加密井区内各射孔小层原始孔隙体积(Vpwell)n;
5)根据所述步骤1)中的动态注采连通分析结果,将与加密井区内的生产井存在动态连通关系的所有注水井的注水量劈分至各射孔小层内,计算注水井在加密井区内各射孔小层总注水量INJn;
6)根据所述步骤4)中的加密井区各射孔小层原始孔隙体积(Vpwell)n和所述步骤5)中注水井在加密井区内各射孔小层总注水量INJn,计算加密井区射孔层的注入PV数PVall;
7)绘制油田注入PV数PVf与油田含水率fwf的关系图版,依据所述步骤6)中的加密井区射孔层的注入PV数PVall,利用图版插值计算加密井初期含水率Fw;
8)根据所述步骤1)计算加密井与注水井动态连通射孔小层总流动系数Msum_infill;
9)根据所述步骤7)中的加密井初期含水率Fw和所述步骤8)中加密井与注水井动态连通射孔小层总流动系数Msum_infill,利用所述步骤3)不同初期含水率下单井可采储量Nplast与动态连通射孔小层总流动系数Msum的关系图版计算加密井可采储量;
其中,所述步骤1)中生产井与注水井动态连通射孔小层总流动系数Msum为:
式中,Abn为生产井在射孔小层n的动态注采连通系数,取值为0或1;G为生产井的射孔总层数;g为生产井射孔层编号,取值为1,2,3,…,G;kp为生产井射孔层的渗透率;hp为生产井射孔层的有效厚度;μo为生产井射孔层内的原油粘度。
2.根据权利要求1所述的可采储量预测方法,其特征在于,所述步骤2)的具体外推方法如下:
①以生产井最后两个月累积产油量作为Np(0)和Np(1),等间距外推Np(2),Np(3),…,Np(j);
②根据所述步骤2)中拟合得到的系数NR、a、q,计算相应的累积产水量Wp(j):
式中,Np(j)为步骤①中外推生产井累积产油量;Wp(j)是与Np(j)对应的生产井累积产水量;NR、a、q为拟合系数;
③判断含水率达到98%时的单井可采储量Nplast,当Np(j)满足如下条件时:
此时Nplast=Np(j);
式中,Np(j)、Np(j-1)均是步骤①中外推生产井累积产油量;Wp(j)、Wp(j-1)分别是与Np(j)、Np(j-1)对应的生产井累积产水量;其中,j为步骤①中外推累积产油量和对应累积产水量的编号。
3.根据权利要求2所述的可采储量预测方法,其特征在于,所述步骤3)中计算单井生产前3个月平均含水率作为初期含水率,根据单井的初期含水率分类建立不同初期含水率下单井可采储量Nplast与动态连通射孔小层总流动系数Msum的关系图版,分类建立原则为将初期含水率相同的生产井划分为同一类建立图版,其中初期含水率按四舍五入原则保留小数点后一位。
4.根据权利要求3所述的可采储量预测方法,其特征在于,所述步骤4)中加密井区内各射孔小层原始孔隙体积(Vpwell)n为:
(Vpwell)n=(Sareahφ)n
式中,Sarea为加密井区内射孔小层n的面积;h为加密井区内射孔小层n的有效厚度;φ为加密井区内射孔小层n的地下孔隙度。
5.根据权利要求4所述的可采储量预测方法,其特征在于,所述步骤5)中注水井在加密井区内各射孔小层总注水量INJn为:
式中,INJn为所有与加密井区内生产井动态连通的注水井在射孔小层n加密井区范围内的注水量;M为与加密井区内生产井动态连通的注水井数量;m为与加密井区内生产井动态连通的注水井编号,取值为1,2,3,…,M;(INJn)m为注水井m在射孔小层n加密井区范围内的注水量;(INJtotal_m)i为注水井m在第i次PLT测试和第i+1次PLT测试期间的累积注水量;(fPLT_m)i为注水井m第i次PLT测试中包含射孔小层n的测试段的吸水比例;NT为注水井PLT测试次数;i为PLT测试次数编号,取值为1,2,3,…,NT;(kwhw)n为注水井m在射孔小层n的地层系数;NN为包含射孔小层n的注水井PLT测试段内所有射孔小层总数;nn为PLT测试段内射孔小层编号,取值为1,2,3,…,NN;(kwhw)nn为注水井m在射孔小层nn的地层系数;[(Wpo)n_i]np为注水井m在第i次PLT测试和第i+1次PLT测试期间生产井np在射孔小层n的累积产液量;NP为加密井区内生产井数量;np为加密井区内生产井编号,取值为1,2,3,…,NP;D为与注水井m存在动态连通关系的生产井总井数;d为与注水井m存在动态连通关系的生产井编号,取值为1,2,3,…,D;[(Wpo)n_i]d为注水井m在第i次PLT测试和第i+1次PLT测试期间生产井d在射孔小层n的累积产液量;[Qtotal_i]np为生产井np在注水井m第i次PLT测试和第i+1次PLT测试期间的累积产液量;(knphnp)n为生产井np在射孔小层n的地层系数;Anp_n为生产井np在射孔小层n与加密井区内注水井的动态注采连通系数,取值为0或1,由所述步骤1)得到;G是生产井射孔层总数;g是生产井射孔层编号,取值为1,2,3,…,G;Anp_g是生产井np在射孔小层g与加密井区内注水井的动态注采连通系数,取值为0或1,由所述步骤1)得到;(knphnp)g为生产井np在射孔小层g的地层系数;[Qtotal_i]d是生产井d在注水井m第i次PLT测试和第i+1次PLT测试期间的累积产液量;(kdhd)n是生产井d在射孔小层n的地层系数;Ad_n是生产井d在射孔小层n与加密井区内注水井的动态注采连通系数,取值为0或1,由所述步骤1)得到;Ad_g是生产井d在射孔小层g与加密井区内注水井的动态注采连通系数,取值为0或1,由所述步骤1)得到;(kdhd)g是生产井d在射孔小层g的地层系数。
6.根据权利要求5所述的可采储量预测方法,其特征在于,所述步骤6)中加密井区射孔层的注入PV数PVall为:
式中,NL为加密井射孔小层总数;n为加密井射孔小层编号,取值为1,2,3,…,NL;Bw为加密井区内射孔小层n的地层水体积系数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110371858.9A CN113065705B (zh) | 2021-04-07 | 2021-04-07 | 海上薄互层砂岩油藏多层合采加密井可采储量预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110371858.9A CN113065705B (zh) | 2021-04-07 | 2021-04-07 | 海上薄互层砂岩油藏多层合采加密井可采储量预测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113065705A CN113065705A (zh) | 2021-07-02 |
CN113065705B true CN113065705B (zh) | 2023-12-08 |
Family
ID=76565882
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110371858.9A Active CN113065705B (zh) | 2021-04-07 | 2021-04-07 | 海上薄互层砂岩油藏多层合采加密井可采储量预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113065705B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116291405B (zh) * | 2023-01-09 | 2024-08-23 | 中国海洋石油集团有限公司 | 一种海上多层合采油藏小层产液量劈分方法及装置 |
CN116044388B (zh) * | 2023-01-09 | 2024-08-06 | 中国海洋石油集团有限公司 | 一种海上水驱多层砂岩油藏小层注采平衡判断方法及装置 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102061910A (zh) * | 2010-06-30 | 2011-05-18 | 中国石油大学(北京) | 砂岩油藏注水开发可预测物理模型建立方法 |
CN102426390A (zh) * | 2011-10-21 | 2012-04-25 | 中国石油大学(北京) | 一种非均质泥砂岩储层储量确定方法 |
CN108446511A (zh) * | 2018-03-28 | 2018-08-24 | 中国石油化工股份有限公司 | 缝洞型油藏油、水的储量预测方法 |
CN111476448A (zh) * | 2019-01-24 | 2020-07-31 | 中国石油化工股份有限公司 | 一种多层合采直井动态产量劈分方法 |
CN112031757A (zh) * | 2020-09-08 | 2020-12-04 | 中国石油天然气股份有限公司 | 一种致密气藏井网加密潜力评价方法 |
CN112464136A (zh) * | 2020-11-25 | 2021-03-09 | 中海石油(中国)有限公司 | 海上薄互层砂岩油田定向井产能及开发效果预测方法 |
CN112541571A (zh) * | 2020-11-25 | 2021-03-23 | 中国石油大学(华东) | 基于双并联神经网络的机器学习的注采连通性确定方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6980940B1 (en) * | 2000-02-22 | 2005-12-27 | Schlumberger Technology Corp. | Intergrated reservoir optimization |
US11009620B2 (en) * | 2019-07-04 | 2021-05-18 | Chengdu University Of Technology | Method for determining favorable time window of infill well in unconventional oil and gas reservoir |
-
2021
- 2021-04-07 CN CN202110371858.9A patent/CN113065705B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102061910A (zh) * | 2010-06-30 | 2011-05-18 | 中国石油大学(北京) | 砂岩油藏注水开发可预测物理模型建立方法 |
CN102426390A (zh) * | 2011-10-21 | 2012-04-25 | 中国石油大学(北京) | 一种非均质泥砂岩储层储量确定方法 |
CN108446511A (zh) * | 2018-03-28 | 2018-08-24 | 中国石油化工股份有限公司 | 缝洞型油藏油、水的储量预测方法 |
CN111476448A (zh) * | 2019-01-24 | 2020-07-31 | 中国石油化工股份有限公司 | 一种多层合采直井动态产量劈分方法 |
CN112031757A (zh) * | 2020-09-08 | 2020-12-04 | 中国石油天然气股份有限公司 | 一种致密气藏井网加密潜力评价方法 |
CN112464136A (zh) * | 2020-11-25 | 2021-03-09 | 中海石油(中国)有限公司 | 海上薄互层砂岩油田定向井产能及开发效果预测方法 |
CN112541571A (zh) * | 2020-11-25 | 2021-03-23 | 中国石油大学(华东) | 基于双并联神经网络的机器学习的注采连通性确定方法 |
Non-Patent Citations (3)
Title |
---|
Combination of Horizontal Wells and Alkaline/Surfactant/Polymer Flooding: A Solution to Produce Hard-to-Recover Residual Oil in Polymer-Flooded Reservoirs;Shize Yin;《SPE Western Regional Meeting》;SPE-190121-MS * |
Horizontal Well Placement Optimization for Heavy-Oil Production in Complex Heterogeneous and Thin Reservoirs: A Case Study from the South China Sea;Donghong Luo et al.;《SPE Heavy and Extra Heavy Oil Conference: Latin America》;SPE-171042-MS * |
特低渗透砂岩油藏剩余油分布的主控因素及有利区块评价;高建等;《中国石油大学学报(自然科学版)》;第31卷(第1期);13-18 * |
Also Published As
Publication number | Publication date |
---|---|
CN113065705A (zh) | 2021-07-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Rogachev et al. | Improving the efficiency of using resource base of liquid hydrocarbons in jurassic deposits of Western Siberia | |
CN104747180B (zh) | 一种用于注水开发的缝洞型油藏分析方法及其应用 | |
CN103472484B (zh) | 基于rs三维敏感地震属性分析的水平井轨迹优化方法 | |
CN113065705B (zh) | 海上薄互层砂岩油藏多层合采加密井可采储量预测方法 | |
CA2707482A1 (en) | A method for performing oilfield production operations | |
CN115587674A (zh) | 油藏改建储气库扩容达产过程气井动态产能预测方法 | |
Mwakipunda et al. | Infill drilling optimization for enhanced oil recovery by waterflooding: A simulation study | |
CN111749688B (zh) | 一种优势渗流通道发育层位和方向的预测方法 | |
Araujo et al. | Production Experience of ICD/AICD for Heavy Oil at Peregrino Field | |
Rafiei | Improved oil production and waterflood performance by water allocation management | |
CN114427444B (zh) | 一种自喷期采油井的井底压力预测方法 | |
CN112464136B (zh) | 海上薄互层砂岩油田定向井产能及开发效果预测方法 | |
Birchenko | Analytical modelling of wells with inflow control devices. | |
Leemput et al. | Full-field reservoir modeling of Central Oman gas-condensate fields | |
CN115637972A (zh) | 一种页岩油水平井开发井排距合理性定量判识方法 | |
Ader et al. | Slaughter Estate Unit Tertiary Miscible Gas Pilot Reservoir Description | |
CN112465218A (zh) | 海上薄互层砂岩油田层系划分及射孔方案优化方法 | |
Aslam et al. | Capacitance Resistance Clustered Model for Mature Peripheral Waterflood Performance Prediction & Optimization | |
CA2671367A1 (en) | A method for performing oilfield production operations | |
Denney | Intelligent-well completions in Agbami: value added and execution performance | |
CN111382523A (zh) | 合采油井产量劈分新方法 | |
RU2301326C1 (ru) | Способ регулирования разработки нефтяного месторождения | |
Weijermars | Fast production and water-breakthrough analysis methods demonstrated using Volve Field data | |
Rafiei et al. | A new analytical approach for water allocation management by determining inter-well connectivity: North Sea Case Study | |
CN111550231B (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 |