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

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

Info

Publication number
CN108729912A
CN108729912A CN201810356826.XA CN201810356826A CN108729912A CN 108729912 A CN108729912 A CN 108729912A CN 201810356826 A CN201810356826 A CN 201810356826A CN 108729912 A CN108729912 A CN 108729912A
Authority
CN
China
Prior art keywords
coefficient
production
numerical simulation
perforation
payzone
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
Application number
CN201810356826.XA
Other languages
English (en)
Other versions
CN108729912B (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 Henan Oilfield Branch Co
Original Assignee
China Petroleum and Chemical Corp
Exploration and Development Research Institute of Sinopec Henan Oilfield Branch 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 Henan Oilfield Branch 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

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)计算生产井的全井产量系数,计算公式为:
式中,M为生产井的全井产量系数,kg;N为生产井的射孔产层总数,整数。
在步骤3中,根据生产井在射孔节点单元的产量系数Mi和全井产量系数M,应用公式(3)计算生产井在各产层的劈分系数,计算公式为:
式中,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)计算生产井的全井产量系数,计算公式为:
式中,M为生产井的全井产量系数,kg;N为生产井的射孔产层总数,整数。
在步骤103中,根据上述得到的生产井在射孔节点单元的产量系数和全井产量系数,应用公式(3)计算生产井在各产层的劈分系数,计算公式为:
式中,ci为生产井在射孔节点单元i的劈分系数,小数。
在油藏数值模拟中,设定一口井在某一给定时间步下的全井产量为Q,利用步骤103中计算得到的各个射孔节点单元的劈分系数ci,根据公式(4)获得各个模拟射孔节点单元的产量Qi,计算公式为:
Qi=ci×Q (4)
流程结束。
本发明中的适用于油藏数值模拟的产量劈分方法,综合反映了产层的储量规模、近井地带的生产状况和流体的饱和程度,对油气藏精细开发具有重要意义。同时,在油气藏数值模拟方面,合适的产量劈分方法,可以准确描述各产层的产出情况,准确的描述井的生产状况,提高油气藏数值模拟精度和可信度,为完善油气藏数值模拟技术奠定基础。

Claims (5)

1.适用于油藏数值模拟的产量劈分方法,其特征在于,该适用于油藏数值模拟的产量劈分方法包括:
步骤1,根据井在各射孔单元的含油饱和度、油相流体密度和各层的有效孔隙体积,计算产层的产量系数;
步骤2,计算多层油藏同时开发时的全井产量系数;
步骤3,应用产层的产量系数和全井产量系数计算产层的劈分系数。
2.根据权利要求1所述的适用于油藏数值模拟的产量劈分方法,其特征在于,在步骤1中,获取井所在的各射孔网格单元的含油饱和度、油相流体的密度和网格单元的有效孔隙体积,应用公式(1)计算各射孔网格单元的产量系数Mi,计算公式为:
Mi=ρo,i·so,i·Vi (1)
式中,Mi为射孔节点单元i的产量系数,kg;ρo,i为射孔节点单元i的油相流体的密度,kg/m3;so,i为射孔节点单元i的含油饱和度;Vi为射孔节点单元i所在产层的有效孔隙体积,m3
3.根据权利要求1所述的适用于油藏数值模拟的产量劈分方法,其特征在于,在步骤2中,根据生产井在所有射孔单元的产量系数Mi,应用公式(2)计算生产井的全井产量系数,计算公式为:
式中,M为生产井的全井产量系数,kg;N为生产井的射孔产层总数,整数。
4.根据权利要求1所述的适用于油藏数值模拟的产量劈分方法,其特征在于,在步骤3中,根据生产井在射孔节点单元的产量系数Mi和全井产量系数M,应用公式(3)计算生产井在各产层的劈分系数,计算公式为:
式中,ci为生产井在射孔节点单元i的劈分系数,小数。
5.根据权利要求4所述的适用于油藏数值模拟的产量劈分方法,其特征在于,在步骤3中,在油藏数值模拟中,设定一口井在某一给定时间步下的全井产量为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 true CN108729912A (zh) 2018-11-02
CN108729912B 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)

Cited By (4)

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

Citations (5)

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

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20160130916A1 (en) * 2014-11-06 2016-05-12 Schlumberger Technology Corporation Local layer geometry engine with work zone generated from buffer defined relative to a wellbore trajectory
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
曹肖萌 等: "薄互储层多层合采条件下的产量劈分新方法研究", 《科学技术与工程》 *

Cited By (6)

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

Also Published As

Publication number Publication date
CN108729912B (zh) 2022-04-19

Similar Documents

Publication Publication Date Title
CN108729912A (zh) 适用于油藏数值模拟的产量劈分方法
CN107829718B (zh) 基于均衡水驱理念的油藏井网及注采方案优化设计方法
CN105224704B (zh) 一种海底不稳定性的评价方法
CN102041995B (zh) 复杂油藏水淹状况监测***
CN104632200B (zh) 河流相多层及厚层油藏精细建模方法
CN103485770B (zh) 一种基于人工神经网络求取含油饱和度的方法及***
CN108089227B (zh) 新的基于三维地震资料的地层孔隙压力预测方法
CN109670220A (zh) 一种基于非结构网格的水平井气水两相数值模拟方法
CN105808311B (zh) 一种基于降维策略的油藏模拟快速拟合方法
CN110056346B (zh) 一种基于趋势变化函数的油藏三维原始含水饱和度模拟方法
CN110838175B (zh) 注气开发油藏的地质模型建立方法
CN107991714A (zh) 基于湖盆古地形恢复的定量方法
CN106150477A (zh) 一种确定单井控制储量的方法
CN110439515A (zh) 注采参数优化方法及装置
CN109684685B (zh) 一种多分支井水合物降压开采条件下的产能及储层稳定性分析方法
CN105095642B (zh) 油藏数值模拟中历史拟合调整参数确定方法及装置
CN105205273B (zh) 致密储层多重介质中油气接力流动的模拟方法及其装置
CN110263434A (zh) 一种基于多尺度混合有限元的流动单元数值模拟方法
CN108518219A (zh) 一种基于天然能量开发油田规模提液选井的评价方法
CN109958413A (zh) 一种特高含水期油藏动态流动单元划分方法
CN106321057B (zh) 水平井压裂裂缝井网模型设计方法
CN108457630A (zh) 利用多配点小压差注水提高油藏采收率的方法
CN104156617B (zh) 用于多层砂岩气藏气层品质分类表征的六阶段建模方法
CN104183018B (zh) 用于带水碳酸盐岩气藏气水分布表征的六阶段建模方法
Zhao et al. INSIM-FPT-3D: A data-driven model for history matching, water-breakthrough prediction and well-connectivity characterization in three-dimensional 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
GR01 Patent grant
GR01 Patent grant