CN107729588A - 考虑氮气注入的数值模拟方法 - Google Patents
考虑氮气注入的数值模拟方法 Download PDFInfo
- Publication number
- CN107729588A CN107729588A CN201610666182.5A CN201610666182A CN107729588A CN 107729588 A CN107729588 A CN 107729588A CN 201610666182 A CN201610666182 A CN 201610666182A CN 107729588 A CN107729588 A CN 107729588A
- Authority
- CN
- China
- Prior art keywords
- nitrogen
- numerical simulation
- nitrogen injection
- data
- 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.)
- Withdrawn
Links
Classifications
-
- 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
Abstract
本发明提供一种考虑氮气注入的数值模拟方法,该考虑氮气注入的数值模拟方法包括:步骤1,建立三相三组分模型;步骤2,确定地下粘温关系和不同相平衡常数;步骤3,确定导热系数与氮气浓度关系,并建立油藏数值模拟模型;步骤4,利用生产历史和数值模拟模型计算的结果进行拟合;步骤5,定义不同的氮气注入量,对计算结果进行优化,确定合理的氮气注入量。该考虑氮气注入的数值模拟方法技术思路清楚、应用简单,具有创新性和实用性,利于推广,为指导现场实施并实现浅薄层经济有效开发具有重要的指导意义。
Description
技术领域
本发明涉及油藏数值模拟领域,特别涉及利用蒸汽+氮气采油提高稠油油藏开发效果的油藏数值模拟方法。
背景技术
胜利油田目前探明稠油地质储量6.4×108t,其中东部油区探明稠油地质储量5.89×108t,目前已动用4.69×108t;西部油区相继发现春风、春晖油田,探明地质储量8029×104t,目前已动用2471×104t。随着新理论与新技术的推广应用,目前,胜利稠油的年产量持续增长,2012年产量达到了475×104t,是胜利油田稳产增产的主力阵地之一。
目前春风油田主要采用HDNS、VDNS、HNS、HDS、VNS等方法进行开发,取得了较好的开发效果。但是不同的区块油藏地质特征不尽相同,原油性质也不完全一样,采用何种开发方式,助剂的注入量的合理范围,由于以上各种开发方式应用时间短,如何在现场应用中通过选择合理的注入参数而是油田实现经济高效开发是面临的最大难题。为此我们发明了一种新的考虑氮气注入的数值模拟方法,解决了以上技术问题。
发明内容
本发明的目的是提供一种可以有效的指导现场实施,实现油藏的高效开发的考虑氮气注入的数值模拟方法。
本发明的目的可通过如下技术措施来实现:考虑氮气注入的数值模拟方法,该考虑氮气注入的数值模拟方法包括:步骤1,建立三相三组分模型;步骤2,确定地下粘温关系和不同相平衡常数;步骤3,确定导热系数与氮气浓度关系,并建立油藏数值模拟模型;步骤4,利用生产历史和数值模拟模型计算的结果进行拟合;步骤5,定义不同的氮气注入量,对计算结果进行优化,确定合理的氮气注入量。
本发明的目的还可通过如下技术措施来实现:
在步骤1中,在模型中定义三相三组分模型,所述的三相三组分是油、水、氮气三个组分和油、水、气三个相态。
在步骤2中,在确定地下粘温关系时,收集三组分的粘度随温度变化的实验数据,该数据中原油粘度的数据为地面脱气原油粘度,收集目标区块的溶解油气比,利用公式对原油粘度数据进行修正,转化为地下原油粘度随温度的变化数据,
A=10.715(5.615R+100)-0.515
B=5.44(5.615R+150)-0.338
其中,μ为含气原油粘度,mPa.s;μ0为脱气原油粘度,mPa.s;R为气油比。
在步骤2中,在确定不同相平衡常数时,利用Wilson’s方程计算各组分在不同比例在不同温度、压力下的平衡常数;
LnKi=5.37(1+ωi)(1-Tci/T)+ln(pci/p)
Ki—组分i的平衡常数;
ωi—组分i的偏心因子;
Tci—组分i的临界温度;
T—温度;
pci—组分i的临界压力;
p—压力。
该考虑氮气注入的数值模拟方法还包括,在步骤2之后,收集油藏的生产和注入历史数据,生成生产动态数据和生产历史文件。
在步骤3中,收集岩石导热系数随氮气的注入浓度变化的实验数据,给出不同氮气浓度下岩石的导热系数。
在步骤3中,根据已知的油藏地质静态资料、粘温资料、相渗资料、相平衡常数以及动态数据,建立油藏数值模拟模型。
在步骤4中,利用生产历史和数值模拟模型计算的结果进行拟合,精度达到10%以内即可进行合理注入量预测,若误差较大,要返回步骤3修改岩石导热系数与氮气浓度关系数据,再次进行拟合计算,直到达到精度要求。
本发明的考虑氮气注入的数值模拟方法,包括1、合理的注氮量;2、注入氮气前后对地层温度的变化;3、注入氮气前后地层压力的变化。该方法在利用蒸汽+氮气开发过程中,利用油藏数值模拟方法确定合理的注汽强度和注氮量,本发明技术思路清楚、应用简单,开发了稠油氮气辅助蒸汽吞吐的数值模拟方法,具有创新性和实用性,利于推广。为指导现场实施并实现浅薄层经济有效开发具有重要的指导意义。
附图说明
图1为本发明的考虑氮气注入的数值模拟方法的一具体实施例的流程图;
图2为本发明的一具体实施例中注入地层后氮气分布图;
图3为本发明的一具体实施例中是否注氮气地层内注入焓的对比曲线;
图4为本发明的一具体实施例中常规蒸汽吞吐与HDNS方法地层压力对比曲线。
具体实施方式
为使本发明的上述和其他目的、特征和优点能更明显易懂,下文特举出较佳实施例,并配合所附图式,作详细说明如下。
如图1所示,图1为本发明的考虑氮气注入的数值模拟方法的流程图。
步骤一:在模型中定义两相三组分模型。所述的三相三组分是油、水、氮气三个组分和油、水、气三个相态。
步骤二:收集三组分的粘度随温度变化的实验数据。该数据中原油粘度的数据为地面脱气原油粘度。收集目标区块的溶解油气比。利用公式对原油粘度数据进行修正,转化为地下原油粘度随温度的变化数据。
A=10.715(5.615R+100)-0.515
B=5.44(5.615R+150)-0.338
μ含气原油粘度,mPa.s;
μ0脱气原油粘度,mPa.s;
R气油比。
步骤三:利用Wilson’s方程计算各组分在不同比例在不同温度、压力下的平衡常数;
LnKi=5.37(1+ωi)(1-Tci/T)+ln(pci/p)
Ki—组分i的平衡常数;
ωi—组分i的偏心因子;
Tci—组分i的临界温度;
T—温度;
pci—组分i的临界压力;
p—压力。
步骤四,收集油藏的生产和注入历史数据,生成生产动态数据和生产历史文件,步骤四中的数据为生产井和注入井的日度或月度数据。
步骤五:收集岩石导热系数随氮气的注入浓度变化的实验数据,给出不同氮气浓度下岩石的导热系数;在一具体实例中由于气体较轻,注入后氮气分布在油层的顶部,改变岩石的导热系数,阻止热量向顶层的散失,图2为注入地层后氮气分布图,由于气体较轻,注入后氮气分布在油层的顶部,阻止热量向顶层的散失;注入氮气后,由于氮气的导热系数低,氮气分布在油层顶部形成隔热“被”,使注入地层的热量较好的保存,图3为是否注氮气地层内注入焓的对比曲线。曲线显示出,注入氮气的模型较未注氮气的热焓高,主要是由于氮气的导热系数低,氮气分布在油层顶部形成隔热“被”,阻止向地层热量散失,而使注入地层的热量较好的保存。
步骤六:根据已知的油藏地质静态资料、粘温资料、相渗资料、相平衡常数以及动态数据,建立油藏数值模拟模型。
步骤七,利用生产历史和数值模拟模型计算的结果进行拟合,精度达到10%以内即可进行合理注入量预测,若误差较大,要返回步骤五修改岩石导热系数与氮气浓度关系数据,再次进行拟合计算,直到达到精度要求。
步骤八,在步骤七中的油藏地质参数基础上,定义不同的氮气注入量,对计算结果进行优化,确定合理的氮气注入量。
在应用本发明的一具体实例中,图2、3、4为本发明的一具体实施例中常规蒸汽吞吐与HDNS方法气体分布、热焓以及地层压力对比曲线,对比了是否注入氮气的结果,从结果看,注入氮气之后地层压力提高0.5~1MPa,同时由于氮气注入地层后主要集中在油藏的顶部,氮气的导热系数较岩石、油、水的导热系数低1个数量级,减少了注汽过程中向顶部盖层的热损失,提高了油藏的热利用率,地层热焓提高7.3%。说明注入的氮气确实起到了提高地层压力和提高地层热利用率的作用,在模型中有较好的体现。
Claims (8)
1.考虑氮气注入的数值模拟方法,该考虑氮气注入的数值模拟方法包括:
步骤1,建立三相三组分模型;
步骤2,确定地下粘温关系和不同相平衡常数;
步骤3,确定导热系数与氮气浓度关系,并建立油藏数值模拟模型;
步骤4,利用生产历史和数值模拟模型计算的结果进行拟合;
步骤5,定义不同的氮气注入量,对计算结果进行优化,确定合理的氮气注入量。
2.根据权利要求1所述的考虑氮气注入的数值模拟方法,其特征在于,在步骤1中,在模型中定义三相三组分模型,所述的三相三组分是油、水、氮气三个组分和油、水、气三个相态。
3.根据权利要求1所述的考虑氮气注入的数值模拟方法,其特征在于,在步骤2中,在确定地下粘温关系时,收集三组分的粘度随温度变化的实验数据,该数据中原油粘度的数据为地面脱气原油粘度,收集目标区块的溶解油气比,利用公式对原油粘度数据进行修正,转化为地下原油粘度随温度的变化数据,
<mrow>
<mi>&mu;</mi>
<mo>=</mo>
<msubsup>
<mi>A&mu;</mi>
<mn>0</mn>
<mi>B</mi>
</msubsup>
</mrow>
A=10.715(5.615R+100)(0.515
B=5.44(5.615R+150)-0.338
其中,μ为含气原油粘度,mPa.s;μ0为脱气原油粘度,mPa.s;R为气油比。
4.根据权利要求1所述的考虑氮气注入的数值模拟方法,其特征在于,在步骤2中,在确定不同相平衡常数时,利用Wilson’s方程计算各组分在不同比例在不同温度、压力下的平衡常数;
LnKi=5.37(1+ωi)(1-Tci/T)+ln(pci/p)
Ki—组分i的平衡常数;
ωi—组分i的偏心因子;
Tci—组分i的临界温度;
T—温度;
pci—组分i的临界压力;
p—压力。
5.根据权利要求1所述的考虑氮气注入的数值模拟方法,其特征在于,该考虑氮气注入的数值模拟方法还包括,在步骤2之后,收集油藏的生产和注入历史数据,生成生产动态数据和生产历史文件。
6.根据权利要求1所述的考虑氮气注入的数值模拟方法,其特征在于,在步骤3中,收集岩石导热系数随氮气的注入浓度变化的实验数据,给出不同氮气浓度下岩石的导热系数。
7.根据权利要求1所述的考虑氮气注入的数值模拟方法,其特征在于,在步骤3中,根据已知的油藏地质静态资料、粘温资料、相渗资料、相平衡常数以及动态数据,建立油藏数值模拟模型。
8.根据权利要求1所述的考虑氮气注入的数值模拟方法,其特征在于,在步骤4中,利用生产历史和数值模拟模型计算的结果进行拟合,精度达到10%以内即可进行合理注入量预测,若误差较大,要返回步骤3修改岩石导热系数与氮气浓度关系数据,再次进行拟合计算,直到达到精度要求。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610666182.5A CN107729588A (zh) | 2016-08-12 | 2016-08-12 | 考虑氮气注入的数值模拟方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610666182.5A CN107729588A (zh) | 2016-08-12 | 2016-08-12 | 考虑氮气注入的数值模拟方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN107729588A true CN107729588A (zh) | 2018-02-23 |
Family
ID=61201455
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610666182.5A Withdrawn CN107729588A (zh) | 2016-08-12 | 2016-08-12 | 考虑氮气注入的数值模拟方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107729588A (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110965992A (zh) * | 2018-09-27 | 2020-04-07 | 中国石油化工股份有限公司 | 一种确定地层含气原油粘度的方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103902752A (zh) * | 2012-12-27 | 2014-07-02 | 中国地质大学(北京) | 一种新型黑油、组分一体化混合数值模拟方法及装置 |
CN104389568A (zh) * | 2014-09-29 | 2015-03-04 | 中国石油大学(北京) | 蒸汽辅助重力泄油过程中气体辅助用量的获取方法及装置 |
CN104899438A (zh) * | 2015-06-02 | 2015-09-09 | 中国地质大学(北京) | 一种基于凝胶泡沫的数值模拟方法 |
US20160122628A1 (en) * | 2014-10-22 | 2016-05-05 | John A. BABCOCK | Y-grade ngl stimulation fluids |
-
2016
- 2016-08-12 CN CN201610666182.5A patent/CN107729588A/zh not_active Withdrawn
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103902752A (zh) * | 2012-12-27 | 2014-07-02 | 中国地质大学(北京) | 一种新型黑油、组分一体化混合数值模拟方法及装置 |
CN104389568A (zh) * | 2014-09-29 | 2015-03-04 | 中国石油大学(北京) | 蒸汽辅助重力泄油过程中气体辅助用量的获取方法及装置 |
US20160122628A1 (en) * | 2014-10-22 | 2016-05-05 | John A. BABCOCK | Y-grade ngl stimulation fluids |
CN104899438A (zh) * | 2015-06-02 | 2015-09-09 | 中国地质大学(北京) | 一种基于凝胶泡沫的数值模拟方法 |
Non-Patent Citations (4)
Title |
---|
CUPLONG: "WinProp相态程序用户指南", 《豆丁》 * |
冯保华: "浅薄稠油油藏HDNS技术数值模拟研究", 《中国优秀硕士学位论文全文数据库_工程科技Ⅰ辑》 * |
刘文章: "高粘原油机械开采工艺的研究", 《石油钻采工艺》 * |
李星: "河南油田超稠油油藏热化学辅助蒸汽吞吐技术研究", 《油气藏评价与开发》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110965992A (zh) * | 2018-09-27 | 2020-04-07 | 中国石油化工股份有限公司 | 一种确定地层含气原油粘度的方法 |
CN110965992B (zh) * | 2018-09-27 | 2023-04-07 | 中国石油化工股份有限公司 | 一种确定地层含气原油粘度的方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Qu et al. | Influence of different fracture morphology on heat mining performance of enhanced geothermal systems based on COMSOL | |
Yuan et al. | Residual oil distribution characteristic of fractured-cavity carbonate reservoir after water flooding and enhanced oil recovery by N2 flooding of fractured-cavity carbonate reservoir | |
CN103590818B (zh) | 一种井筒内温度分布半解析确定方法及装置 | |
CN100419208C (zh) | 非均质油层多段塞等流度聚能平行同步驱油方法 | |
CN103670350B (zh) | 一种碳酸盐岩缝洞型油藏变强度注水开采方法 | |
Davarpanah et al. | A simulation study to control the oil production rate of oil-rim reservoir under different injectivity scenarios | |
CN103590799B (zh) | 一种致密砂岩储层水平井三维压裂设计方法 | |
CN104389568B (zh) | 蒸汽辅助重力泄油过程中气体辅助用量的获取方法及装置 | |
Wang et al. | Dynamic fractures are an emerging new development geological attribute in water-flooding development of ultra-low permeability reservoirs | |
CN106437674A (zh) | 仿水平井注水开发井网适配方法 | |
Huang et al. | Heat extraction and power production forecast of a prospective Enhanced Geothermal System site in Songliao Basin, China | |
CN103590798B (zh) | 一种超稠油注蒸汽采油焖井时间的确定方法及计算装置 | |
CN109543206A (zh) | 稠油热采水平井储层经济需热量优化方法 | |
CN107729587A (zh) | 考虑降粘剂注入的数值模拟方法 | |
CN104727790A (zh) | 水驱油藏剩余油再富集成藏周期预测方法 | |
Tan et al. | Flow model of a multi-stage hydraulic fractured horizontal well based on tree-shaped fractal fracture networks | |
CN109763805A (zh) | 一种深层页岩气螺旋式变参数压裂方法 | |
CN106295235B (zh) | 一种页岩地层天然裂缝剪切滑移量的计算方法 | |
CN107729588A (zh) | 考虑氮气注入的数值模拟方法 | |
Altunina et al. | Promising physical-chemical IOR technologies for Arctic oilfields | |
CN106468160B (zh) | 一种确定co2驱泡沫流油组分的方法以及co2驱的模拟方法 | |
Guo et al. | Oil production rate predictions for steam assisted gravity drainage based on high-pressure experiments | |
Ziabakhsh-Ganji et al. | Investigation of the synergy potential of oil and geothermal energy from a fluvial oil reservoir | |
CN106050204A (zh) | 一种考虑岩石孔隙压缩系数温敏效应的稠油热采模拟方法 | |
CN104675372A (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 | ||
WW01 | Invention patent application withdrawn after publication | ||
WW01 | Invention patent application withdrawn after publication |
Application publication date: 20180223 |