CN108729912B - 适用于油藏数值模拟的产量劈分方法 - Google Patents

适用于油藏数值模拟的产量劈分方法 Download PDF

Info

Publication number
CN108729912B
CN108729912B CN201810356826.XA CN201810356826A CN108729912B CN 108729912 B CN108729912 B CN 108729912B CN 201810356826 A CN201810356826 A CN 201810356826A CN 108729912 B CN108729912 B CN 108729912B
Authority
CN
China
Prior art keywords
yield
well
coefficient
perforation
oil
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
CN201810356826.XA
Other languages
English (en)
Other versions
CN108729912A (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 Petroleum and Chemical Corp
Exploration and Development Research Institute of Sinopec Shengli Oilfield Co
Original Assignee
China Petroleum and Chemical Corp
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 Petroleum and Chemical Corp, Exploration and Development Research Institute of Sinopec Shengli Oilfield Co filed Critical China Petroleum and Chemical Corp
Priority to CN201810356826.XA priority Critical patent/CN108729912B/zh
Publication of CN108729912A publication Critical patent/CN108729912A/zh
Application granted granted Critical
Publication of CN108729912B publication Critical patent/CN108729912B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • EFIXED CONSTRUCTIONS
    • E21EARTH OR ROCK DRILLING; MINING
    • E21BEARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
    • E21B49/00Testing 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

Landscapes

  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Geology (AREA)
  • Mining & Mineral Resources (AREA)
  • Physics & Mathematics (AREA)
  • Environmental & Geological Engineering (AREA)
  • Fluid Mechanics (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Geochemistry & Mineralogy (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)
  • Filling Or Discharging Of Gas Storage Vessels (AREA)
  • Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)

Abstract

本发明提供一种适用于油藏数值模拟的产量劈分方法,该适用于油藏数值模拟的产量劈分方法包括:步骤1,根据井在各射孔单元的含油饱和度、油相流体密度和各层的有效孔隙体积,计算产层的产量系数;步骤2,计算多层油藏同时开发时的全井产量系数;步骤3,应用产层的产量系数和全井产量系数计算产层的劈分系数。该适用于油藏数值模拟的产量劈分方法建立了能够准确描述油藏开发及油藏数值模拟中的各产层的产出情况,对油气藏精细开发具有重要意义,为完善油藏数值模拟技术奠定基础。

Description

适用于油藏数值模拟的产量劈分方法
技术领域
本发明涉及油藏开发技术领域,特别是涉及到一种适用于油藏数值模拟的产量劈分方法。
背景技术
常规的KH产量劈分方法主要是考虑小层的渗透率和油层厚度的影响,将生产井的产量劈分到各小层。但是该方法不能识别储层性质。在常规的KH劈分方法的基础上引入相对渗透率的影响,称之为KHK劈分方法,该方法对地层系数的计算进行了改进,但是对于储层非均质较强的薄互层砂岩油藏适用性较差。KNK产量劈分方法是在KHK方法的基础上用油井小层控制地质储量对该模型进行修正得到的。上述方法均是在考虑小层流动能力的基础上建立的,没有反映近井地带流体的饱和程度。为此,我们发明了一种适用于油藏数值模拟的产量劈分方法,综合反映产层的储量规模、近井地带的生产状况和流体的饱和程度。为此我们发明了一种新的适用于油藏数值模拟的产量劈分方法,解决了以上技术问题。
发明内容
本发明的目的是提供一种综合反映产层的储量规模、近井地带的生产状况和流体的饱和程度,为完善油藏数值模拟技术提供技术支撑的适用于油藏数值模拟的产量劈分方法。
本发明的目的可通过如下技术措施来实现:适用于油藏数值模拟的产量劈分方法,该适用于油藏数值模拟的产量劈分方法包括:步骤1,根据井在各射孔单元的含油饱和度、油相流体密度和各层的有效孔隙体积,计算产层的产量系数;步骤2,计算多层油藏同时开发时的全井产量系数;步骤3,应用产层的产量系数和全井产量系数计算产层的劈分系数。
本发明的目的还可通过如下技术措施来实现:
在步骤1中,获取井所在的各射孔网格单元的含油饱和度、油相流体的密度和网格单元的有效孔隙体积,应用公式(1)计算各射孔网格单元的产量系数Mi,计算公式为:
Mi=ρo,i·so,i·Vi (1)
式中,Mi为射孔节点单元i的产量系数,kg;ρo,i为射孔节点单元i的油相流体的密度,kg/m3;so,i为射孔节点单元i的含油饱和度;Vi为射孔节点单元i所在产层的有效孔隙体积,m3
在步骤2中,根据生产井在所有射孔单元的产量系数Mi,应用公式(2)计算生产井的全井产量系数,计算公式为:
Figure BDA0001634249780000021
式中,M为生产井的全井产量系数,kg;N为生产井的射孔产层总数,整数。
在步骤3中,根据生产井在射孔节点单元的产量系数Mi和全井产量系数M,应用公式(3)计算生产井在各产层的劈分系数,计算公式为:
Figure BDA0001634249780000022
式中,ci为生产井在射孔节点单元i的劈分系数,小数。
在步骤3中,在油藏数值模拟中,设定一口井在某一给定时间步下的全井产量为Q,利用各个射孔节点单元的劈分系数ci,根据公式(4)获得各个模拟射孔节点单元的产量Qi,计算公式为:
Qi=ci×Q (4)。
本发明中的适用于油藏数值模拟的产量劈分方法,是油气藏开发技术、研究石油和天然气等地下流体在储层内渗流特征及油气藏数值模拟技术的重要方法。产量劈分是多层油藏精细开发、准确描述各产层储量动用情况及采出程度的重要方法,是油藏数值模拟技术是否可以准确描述各产层对产量贡献程度的重要方法。该方法针对油气藏数值模拟中产量劈分法的问题,在获得各产层平均油相流体密度、含油饱和度和有效体积的基础上,应用产量劈分模型,计算各产层的产量劈分系数。该方法建立了能够准确描述油藏开发及油藏数值模拟中的各产层的产出情况,对油气藏精细开发具有重要意义,为完善油藏数值模拟技术奠定基础。
附图说明
图1为本发明的适用于油藏数值模拟的产量劈分方法的一具体实施例的流程图。
具体实施方式
为使本发明的上述和其他目的、特征和优点能更明显易懂,下文特举出较佳实施例,并配合附图所示,作详细说明如下。
如图1所示,图1为本发明的适用于油藏数值模拟的产量劈分方法的流程图。
在步骤101中,获取井所在的各射孔网格单元的含油饱和度、油相流体的密度和网格单元的有效孔隙体积,应用公式(1)计算各射孔网格单元的产量系数Mi,计算公式为:
Mi=ρo,i·so,i·Vi (1)
式中,Mi为射孔节点单元i的产量系数,kg;ρo,i为射孔节点单元i的油相流体的密度,kg/m3;so,i为射孔节点单元i的含油饱和度;Vi为射孔节点单元i所在产层的有效孔隙体积,m3
在步骤102中,根据上述得到的生产井在所有射孔单元的产量系数Mi,应用公式(2)计算生产井的全井产量系数,计算公式为:
Figure BDA0001634249780000031
式中,M为生产井的全井产量系数,kg;N为生产井的射孔产层总数,整数。
在步骤103中,根据上述得到的生产井在射孔节点单元的产量系数和全井产量系数,应用公式(3)计算生产井在各产层的劈分系数,计算公式为:
Figure BDA0001634249780000032
式中,ci为生产井在射孔节点单元i的劈分系数,小数。
在油藏数值模拟中,设定一口井在某一给定时间步下的全井产量为Q,利用步骤103中计算得到的各个射孔节点单元的劈分系数ci,根据公式(4)获得各个模拟射孔节点单元的产量Qi,计算公式为:
Qi=ci×Q (4)
流程结束。
本发明中的适用于油藏数值模拟的产量劈分方法,综合反映了产层的储量规模、近井地带的生产状况和流体的饱和程度,对油气藏精细开发具有重要意义。同时,在油气藏数值模拟方面,合适的产量劈分方法,可以准确描述各产层的产出情况,准确的描述井的生产状况,提高油气藏数值模拟精度和可信度,为完善油气藏数值模拟技术奠定基础。

Claims (1)

1.适用于油藏数值模拟的产量劈分方法,其特征在于,该适用于油藏数值模拟的产量劈分方法包括:
步骤1,根据井在各射孔网格单元的含油饱和度、油相流体密度和各射孔网格单元的有效孔隙体积,计算射孔网格单元的产量系数;
步骤2,计算多层油藏同时开发时的全井产量系数;
步骤3,应用各射孔网格单元的产量系数和全井产量系数计算各射孔网格单元的劈分系数;
在步骤1中,获取井所在的各射孔网格单元的含油饱和度、油相流体的密度和各射孔网格单元的有效孔隙体积,应用公式(1)计算各射孔网格单元的产量系数Mi,计算公式为:
Mi=ρo,i·so,i·Vi (1)
式中,Mi为射孔网格单元i的产量系数,kg;ρo,i为射孔网格单元i的油相流体的密度,kg/m3;so,i为射孔网格单元i的含油饱和度;Vi为射孔网格单元i所在各射孔网格单元的有效孔隙体积,m3
在步骤2中,根据生产井在所有射孔网格单元的产量系数Mi,应用公式(2)计算生产井的全井产量系数,计算公式为:
Figure FDA0003472606180000011
式中,M为生产井的全井产量系数,kg;N为生产井的射孔产层总数,整数;
在步骤3中,根据生产井在射孔网格单元的产量系数Mi和全井产量系数M,应用公式(3)计算生产井在各射孔网格单元的劈分系数,计算公式为:
Figure FDA0003472606180000012
式中,ci为生产井在射孔网格单元i的劈分系数,小数;
在油藏数值模拟中,设定一口井在某一给定时间步下的全井产量为Q,利用各射孔网格单元的劈分系数ci,根据公式(4)获得各射孔网格单元的产量Qi,计算公式为:
Qi=ci×Q (4)。
CN201810356826.XA 2018-04-19 2018-04-19 适用于油藏数值模拟的产量劈分方法 Active CN108729912B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810356826.XA CN108729912B (zh) 2018-04-19 2018-04-19 适用于油藏数值模拟的产量劈分方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810356826.XA CN108729912B (zh) 2018-04-19 2018-04-19 适用于油藏数值模拟的产量劈分方法

Publications (2)

Publication Number Publication Date
CN108729912A CN108729912A (zh) 2018-11-02
CN108729912B true CN108729912B (zh) 2022-04-19

Family

ID=63939689

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810356826.XA Active CN108729912B (zh) 2018-04-19 2018-04-19 适用于油藏数值模拟的产量劈分方法

Country Status (1)

Country Link
CN (1) CN108729912B (zh)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109577944B (zh) * 2018-11-28 2022-10-21 王国亭 多层系低渗致密砂岩气井各层动态控制储量测定方法
CN109630102B (zh) * 2018-11-28 2022-11-11 郭建林 多层系低渗致密砂岩气井各层含气面积测定方法
CN110984976B (zh) * 2019-12-20 2023-02-28 陕西延长石油(集团)有限责任公司研究院 多层合采气井产量劈分方法
CN111382523B (zh) * 2020-03-17 2023-06-16 中国石油化工股份有限公司 合采油井产量劈分新方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104933246A (zh) * 2015-06-15 2015-09-23 中国海洋石油总公司 一种定量表征气田储量参数对储量变化贡献的方法
CN107339087A (zh) * 2017-08-10 2017-11-10 中国石油天然气股份有限公司 一种注水量劈分方法及装置
CN107795319A (zh) * 2016-08-29 2018-03-13 中国石油天然气股份有限公司 具有重力分异现象的油藏地质储量获取方法及装置
CN107869348A (zh) * 2017-10-27 2018-04-03 西北大学 一种厚层砂岩油藏生产井产量劈分的方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB2533847B (en) * 2014-11-06 2017-04-05 Logined Bv Local layer geometry engine with work zone generated from buffer defined relative to a wellbore trajectory

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104933246A (zh) * 2015-06-15 2015-09-23 中国海洋石油总公司 一种定量表征气田储量参数对储量变化贡献的方法
CN107795319A (zh) * 2016-08-29 2018-03-13 中国石油天然气股份有限公司 具有重力分异现象的油藏地质储量获取方法及装置
CN107339087A (zh) * 2017-08-10 2017-11-10 中国石油天然气股份有限公司 一种注水量劈分方法及装置
CN107869348A (zh) * 2017-10-27 2018-04-03 西北大学 一种厚层砂岩油藏生产井产量劈分的方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
薄互储层多层合采条件下的产量劈分新方法研究;曹肖萌 等;《科学技术与工程》;20140508;第14卷(第13期);166-171页 *

Also Published As

Publication number Publication date
CN108729912A (zh) 2018-11-02

Similar Documents

Publication Publication Date Title
CN108729912B (zh) 适用于油藏数值模拟的产量劈分方法
CN105354639B (zh) 致密油多重介质耦合渗流的全周期产能预测方法及其装置
CN105243182B (zh) 致密油压裂水平井的动态储量计算方法
CN104948163B (zh) 一种页岩气井产能测定方法
CN106481332B (zh) 用于确定页岩气多段压裂水平井内外区动态储量的方法
CN110334431A (zh) 一种低渗透致密气藏单井控制储量计算及剩余气分析方法
CN105386751B (zh) 一种基于油藏渗流模型的水平井测井产能预测方法
CN105404735B (zh) 裂缝与基质对超低渗透油藏单井产量贡献率的定量评价方法
CN104809277A (zh) 一种超低渗致密储层水平井地质建模方法
Kruger Determining areal permeability distribution by calculations
CN114427432B (zh) 一种气藏剩余气开发潜力确定方法
CN111353205B (zh) 用于致密气藏产水气井地层压力和动态产能的计算方法
CN109684685B (zh) 一种多分支井水合物降压开采条件下的产能及储层稳定性分析方法
US10359542B2 (en) Generating dynamically calibrated geo-models in green fields
CN103485770B (zh) 一种基于人工神经网络求取含油饱和度的方法及***
EP3094818A2 (en) Shared equation of state characterization of multiple fluids
CN108133086A (zh) 一种应力敏感储层中产水压裂气井裂缝半长反演方法
CN109958413A (zh) 一种特高含水期油藏动态流动单元划分方法
CN111476448B (zh) 一种多层合采直井动态产量劈分方法
CN113761780A (zh) 一种综合生产测井资料的产层剩余油饱和度最优化计算方法
CN109723435B (zh) 确定页岩气藏水平井重复压裂前现地应力场的方法及设备
CN112836438B (zh) 一种非均质边水气藏气井见水时间预测方法及装置
CN107288633A (zh) 具有重力分异现象的凝析气藏地质储量计算方法
CN111188613B (zh) 一种致密气藏气井井控半径确定方法及***
CN110863818A (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