CN108331565A - 一种改进型火驱数值模拟动力学模型的建模方法 - Google Patents

一种改进型火驱数值模拟动力学模型的建模方法 Download PDF

Info

Publication number
CN108331565A
CN108331565A CN201810062416.4A CN201810062416A CN108331565A CN 108331565 A CN108331565 A CN 108331565A CN 201810062416 A CN201810062416 A CN 201810062416A CN 108331565 A CN108331565 A CN 108331565A
Authority
CN
China
Prior art keywords
numerical simulation
fireflood
oil
modified
kinetic model
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.)
Pending
Application number
CN201810062416.4A
Other languages
English (en)
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 National Offshore Oil Corp CNOOC
CNOOC Research Institute Co Ltd
Original Assignee
China National Offshore Oil Corp CNOOC
CNOOC Research Institute Co Ltd
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 National Offshore Oil Corp CNOOC, CNOOC Research Institute Co Ltd filed Critical China National Offshore Oil Corp CNOOC
Priority to CN201810062416.4A priority Critical patent/CN108331565A/zh
Publication of CN108331565A publication Critical patent/CN108331565A/zh
Pending legal-status Critical Current

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
    • E21B43/00Methods or apparatus for obtaining oil, gas, water, soluble or meltable materials or a slurry of minerals from wells
    • E21B43/16Enhanced recovery methods for obtaining hydrocarbons
    • E21B43/24Enhanced recovery methods for obtaining hydrocarbons using heat, e.g. steam injection
    • E21B43/243Combustion in situ
    • 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
    • E21B47/00Survey of boreholes or wells
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION 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/00Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/02Agriculture; Fishing; Forestry; Mining

Landscapes

  • Engineering & Computer Science (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Mining & Mineral Resources (AREA)
  • Physics & Mathematics (AREA)
  • Geology (AREA)
  • Business, Economics & Management (AREA)
  • Environmental & Geological Engineering (AREA)
  • Fluid Mechanics (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Geochemistry & Mineralogy (AREA)
  • Marine Sciences & Fisheries (AREA)
  • Human Resources & Organizations (AREA)
  • Agronomy & Crop Science (AREA)
  • Geophysics (AREA)
  • Health & Medical Sciences (AREA)
  • Economics (AREA)
  • General Health & Medical Sciences (AREA)
  • Animal Husbandry (AREA)
  • Marketing (AREA)
  • Primary Health Care (AREA)
  • Strategic Management (AREA)
  • Tourism & Hospitality (AREA)
  • General Business, Economics & Management (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明涉及一种改进型火驱数值模拟动力学模型的建模方法,包括以下步骤:1)利用加速量热仪和热重分析仪求取原油的氧化动力学参数,同时确定原油中的重质组分以及轻质组分的分子量;2)利用热物性分析仪求取热物性参数;3)建立与实际油藏条件相同的孔隙度、渗透率以及含油饱和度条件,进行一维燃烧管实验,取得的一维燃烧管实验结果;4)利用CMG油藏数值模拟软件对步骤1)得到的氧化动力学参数以及步骤2)得到的热物性参数与步骤3)得到的一维燃烧管实验结果进行油藏数值模拟历史拟合,即可得到改进型火驱数值模拟动力学模型,该改进型火驱数值模拟动力学模型相比传统模型组分减少、方程式简化,不仅提高了运算精度,而且大大降低了运算时间。

Description

一种改进型火驱数值模拟动力学模型的建模方法
技术领域
本发明涉及一种建模方法,特别是关于一种改进型火驱数值模拟动力学模型的建模方法,属于油气田开发工程技术领域。
背景技术
火烧油层,亦称火驱,是指通过人工手段点燃油层,并向油层内源源不断地注入空气,利用稠油重质组分作为燃料,形成高温氧化前缘,受热的原油粘度不断降低并通过燃烧产生的烟道气不断驱动而产出的过程。火烧油层是稠油热采的关键技术之一,在罗马尼亚、北美等国家已经进行大规模的商业化开采。“十一五”以来,已在我国新疆油田、辽河油田等多个区块进行火烧油层矿场试验,在稠油注蒸汽开发采出程度有限的前提下,火烧油层作为其接替技术是极具潜力的一种热采方式。
火烧油层技术相比于其它热采方式机理更加复杂、工艺水平要求更高,开发生产规律需要重新进行认识。通常,在火烧油层开发过程中能够进行精确的数值模拟并进行方案设计以及动态跟踪调整是火驱矿场试验成败的关键。目前在国内外相关的火烧油层数值模拟研究工作较少,基于阿隆尼乌斯方程的传统动力学模型难以精确刻画火驱储层区带分布特征、预测并认识火驱生产动态规律。因此如何能建立一种精确描述火驱数值模拟动力学模型值得深入研究。
发明内容
针对上述问题,本发明的目的是提供一种改进型火驱数值模拟动力学模型的建模方法,该火驱数值模拟动力学模型能够降低运算时间并提高运算精度。
为实现上述目的,本发明采取以下技术方案:一种改进型火驱数值模拟动力学模型的建模方法,其包括以下步骤:
1)利用加速量热仪和热重分析仪求取原油的氧化动力学参数,包括:指前因子、反应活化能、放热量和焦炭摩尔数,同时确定原油中的重质组分以及轻质组分的分子量;
2)利用热物性分析仪求取热物性参数,包括:岩石体积热容、岩石热传导率和原油热传导率等;
3)进行原油的一维燃烧管实验,且燃烧管的渗透率、孔隙度、含油饱和度等油藏参数与实际油藏保持一致,取得的一维燃烧管实验结果包括:累产油量、产油速度、产出气体含量、燃烧前缘推进速度和不同测温点温度等;
4)利用CMG油藏数值模拟软件对步骤1)得到的氧化动力学参数以及步骤2)得到的热物性参数与步骤3)得到的一维燃烧管实验结果进行油藏数值模拟历史拟合,具体为:
①将原油组分划分为焦炭与可动油,参加化学反应的组分为焦炭,未参加化学反应的组分为可动油;
②根据步骤1)得到的不同原油中的重质组分以及轻质组分的分子量,建立拟合所依托的化学方程反应式(1)和(2):
1Oil1→(1-x)Oil2+xCoke (1)
Coke+O2→CO2+H2O (2)
上述式中,Oil1代表原始油;Oil2代表未参加反应的可动油;Coke代表焦炭;x代表焦炭摩尔数,0<x<1;
③将化学方程反应式(1)和(2)代入CMG油藏数值模拟软件的火驱数值模拟文件(即CMG油藏数值模拟软件读取的数值模拟体dat文件)中,并不断调试化学方程反应式(2)的配平系数、反应活化能以及放热量;
④当调试到与一维燃烧管实验结果:产油速度、累产油量、产出气体含量以及测温点温度等相吻合的过程时,则停止调试,此时经过调试化学方程反应式(1)和(2)即为改进型火驱数值模拟动力学模型。
本发明由于采取以上技术方案,其具有以下优点:本发明从火驱机理上重新认识了火烧油层过程,避免了传统火驱数值模拟动力学模型依赖阿隆尼乌斯方程导致的拟合精度不精确以及多组分导致的运算时间过长的不足;同时,本发明将燃烧过程划分为燃料沉积以及可动油流动的过程,组分数和方程数比传统火驱数值模拟动力学模型少2个,运算时间提高一倍,在相同网格数的前提下运算精度大幅提高。本发明可以提高油藏工程师的工作效率,对国家石油、能源安全具有重要的战略意义。
附图说明
图1是两种模型同一时刻一维燃烧管数值模拟结果对比图;
图2是两种模型不同网格数下瞬时产油量对比曲线图;
图3是线性交错井排两种模型同一时刻二维含油饱和度场分布图;
图4是两种模型火驱日产油量及产出端含氧量变化曲线图;
图5是1/4五点法井网两种模型同一时刻氧气浓度分布图;
图6是三维非均质模型纵向渗透率分布图;
图7是三维非均质条件下两种模型同一时刻温度分布图。
具体实施方式
下面结合附图和实施例对本发明进行详细的描述。然而应当理解,附图的提供仅为了更好地理解本发明,它们不应该理解成对本发明的限制。
实施例1:
建立多元热流体吞吐后转火驱开发的南堡35-2油田氧化动力学模型:
(1)首先利用加速量热仪和热重分析仪求取原油的氧化反应指前因子1500s-1Pa-1,反应活化能6×104J/mol,反应放热5.2×105J/mol,焦炭摩尔数50.9kg/m3
(2)利用热物性分析仪求取岩石体积热容2.575×106J/m3·℃,岩石热传导率为1.634×104J/m·day·℃,原油热传导率9.77×103J/m·day·℃;
(3)建立与实际油藏条件相同的孔隙度、渗透率以及含油饱和度条件,进行一维燃烧管实验;
(4)利用步骤(1)和(2)所得到的原油氧化动力学参数以及热物性参数与步骤(3)的一维燃烧管结果进行拟合;拟合过程中将原油定义为两种组分:焦炭与可动油,根据步骤(1)所得到的焦炭摩尔数,确定南堡35-2油田拟合方程中x数值为0.15,各组分基本参数如表1,最终的拟合方程为:
1Oil1→0.85Oil2+0.15Coke
Coke+1.18O2→0.8CO2+1.11H2O
表1 南堡35-2区块火驱数值模拟原油组分基本参数
下面对传统火驱数值模拟动力学模型(以下称传统模型)与改进型火驱数值模拟动力学模型(以下称改进模型)拟合的结果进行对比:
1)一维燃烧管模型
取长度为0.8m两端分别注气、采油的燃烧管模型对不同网格数的传统模型以及改进型模型的火驱模拟结果进行对比,两种模型同一时刻模拟结果如图1、图2所示。从图中可见,网格数为10的传统8组分、4方程模型燃烧前缘的传播严重滞后,模拟结果不准确,这是由于粗化的传统模型引入太多平均值,焦炭量增加,燃烧峰值温度增大导致。而10个网格的改进模型与细化后的传统模型(800网格)结果基本一致,由于组分以及方程式减少,其计算时间(如表2所示)大大降低。
表2 两种模型及火驱数值模拟运算时间
2)二维均质模型
建立线性交错井排(2注1采)模型,分别利用传统模型以及改进模型对不同网格数的火驱数值模拟结果进行对比,两种模型的同一时刻含油饱和度数值模拟结果如图3、图4所示。从图中可以看出,网格数为20×10的改进模型含油饱和度模拟结果与200×100的传统模型刻画的氧化前缘以及油墙位置基本一致,精确描述了火驱剩余油分布,而从表2来看,其运算时间大大降低。
实施例2:
建立旅大21-2油田火驱开发数值模拟模型:
(1)首先利用加速量热仪和热重分析仪求取原油的氧化反应指前因子33.4s-1Pa-1,活化能14.8×104J/mol,反应放热6.7×104J/mol,焦炭摩尔数15kg/m3
(2)利用热物性分析仪求取岩石体积热容2.35×106J/m3·℃,岩石热传导率3×105J/m·day·℃,原油热传导率1.2×104J/m·day·℃;
(3)建立与实际油藏条件相同的孔隙度、渗透率以及含油饱和度条件,进行一维燃烧管实验;
(4)利用步骤(1)和(2)所得到的原油氧化动力学参数以及岩石热物性参数与步骤(3)的一维燃烧管结果进行拟合;拟合过程中将原油定义为两种组分:焦炭与可动油,根据(1)所得到的焦炭摩尔数,确定南堡35-2油田拟合方程中x数值为0.20,各组分分子量如表3,最终的拟合方程为:
1Oil1→0.80Oil2+0.20Coke
1Coke+O2→0.5H2O+0.5CO2
表3 旅大21-2区块火驱数值模拟原油组分基本参数
下面对传统模型与改进模型拟合的结果进行对比:
1)均质二维模型
如图5所示,建立1/4五点法井网(1注1采)模型,对比不同网格数的传统模型与改进模型数值模拟结果。从图中可见,20×20网格数的改进模型与200×200的传统模型刻画的氧化前缘位置基本吻合,从表2中可以看出其运算时间大大降低。而粗化后的20×20的传统模型刻画的氧化前缘由于计算过程中引用了太多平均值而导致描述不精确。
2)非均质三维模型
如图6所示,截取旅大21-2油田非均质性较强的一部分油藏模型,从图中可以看出油藏纵向上存在着一定的渗透率差异。对不同网格数的传统模型以及改进模型进行数值模拟结果对比,数值模拟结果如图7所示,从图中可以看出,粗化后(网格数14×20)的改进模型与传统模型(网格数70×100)刻画的氧气含量分布图基本一致,同时精确地描述了由于纵向渗透率差异导致的空气超覆现象,同样的从表2中可以看出,使用粗化后的改进模型数值模拟运算时间大大降低。
上述各实施例仅用于说明本发明,其中各部件的结构、连接方式和制作工艺等都是可以有所变化的,凡是在本发明技术方案的基础上进行的等同变换和改进,均不应排除在本发明的保护范围之外。

Claims (5)

1.一种改进型火驱数值模拟动力学模型的建模方法,其特征在于,包括以下步骤:
1)利用加速量热仪和热重分析仪求取原油的氧化动力学参数,同时确定原油中的重质组分以及轻质组分的分子量;
2)利用热物性分析仪求取热物性参数;
3)建立与实际油藏条件相同的孔隙度、渗透率以及含油饱和度条件,进行一维燃烧管实验,取得的一维燃烧管实验结果;
4)利用CMG油藏数值模拟软件对步骤1)得到的氧化动力学参数以及步骤2)得到的热物性参数与步骤3)得到的一维燃烧管实验结果进行油藏数值模拟历史拟合,即可得到改进型火驱数值模拟动力学模型。
2.如权利要求1所述的一种改进型火驱数值模拟动力学模型的建模方法,其特征在于,上述步骤1)中,原油的氧化动力学参数包括:指前因子、反应活化能、放热量和焦炭摩尔数。
3.如权利要求2所述的一种改进型火驱数值模拟动力学模型的建模方法,其特征在于,上述步骤2)中,热物性参数包括:岩石体积热容、岩石热传导率和原油热传导率。
4.如权利要求3所述的一种改进型火驱数值模拟动力学模型的建模方法,其特征在于,上述步骤3)中,一维燃烧管实验结果包括:累产油量、产油速度、产出气体含量、燃烧前缘推进速度和不同测温点温度。
5.如权利要求4所述的一种改进型火驱数值模拟动力学模型的建模方法,其特征在于,上述步骤4)中,油藏数值模拟历史拟合包括以下步骤:
①将原油组分划分为焦炭与可动油,参加化学反应的组分为焦炭,未参加化学反应的组分为可动油;
②根据步骤1)得到的不同原油中的重质组分以及轻质组分的分子量,建立拟合所依托的化学方程反应式(1)和(2):
1Oil1→(1x)Oil2+xCoke (1)
Coke+O2→CO2+H2O (2)
上述式中,Oil1代表原始油;Oil2代表未参加反应的可动油;Coke代表焦炭;x代表焦炭摩尔数,0<x<1;
③将化学方程反应式(1)和(2)代入CMG油藏数值模拟软件的火驱数值模拟文件中,并不断调试化学方程反应式(2)的配平系数、反应活化能以及放热量;
④当调试到与一维燃烧管实验结果:产油速度、累产油量、产出气体含量以及测温点温度相吻合的过程时,则停止调试,此时经过调试化学方程反应式(1)和(2)即为改进型火驱数值模拟动力学模型。
CN201810062416.4A 2018-01-23 2018-01-23 一种改进型火驱数值模拟动力学模型的建模方法 Pending CN108331565A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810062416.4A CN108331565A (zh) 2018-01-23 2018-01-23 一种改进型火驱数值模拟动力学模型的建模方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810062416.4A CN108331565A (zh) 2018-01-23 2018-01-23 一种改进型火驱数值模拟动力学模型的建模方法

Publications (1)

Publication Number Publication Date
CN108331565A true CN108331565A (zh) 2018-07-27

Family

ID=62925652

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810062416.4A Pending CN108331565A (zh) 2018-01-23 2018-01-23 一种改进型火驱数值模拟动力学模型的建模方法

Country Status (1)

Country Link
CN (1) CN108331565A (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112196514A (zh) * 2020-09-30 2021-01-08 中国海洋石油集团有限公司 一种利用热重分析仪测量稠油注空气开发燃料沉积量的方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20020191004A1 (en) * 2000-08-09 2002-12-19 Ebersole John Franklin Method for visualization of hazards utilizing computer-generated three-dimensional representations
CN101798915A (zh) * 2010-04-26 2010-08-11 胡士清 一种利用分段火驱开采提高厚层块状稠油采收率的方法
CN101818637A (zh) * 2010-04-26 2010-09-01 胡士清 一种通过控制火烧注气速度提高厚层块状稠油油藏采收率的方法
CN104060975A (zh) * 2014-06-24 2014-09-24 中国石油大学(北京) 稠油燃烧过程中活化能的预测方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20020191004A1 (en) * 2000-08-09 2002-12-19 Ebersole John Franklin Method for visualization of hazards utilizing computer-generated three-dimensional representations
CN101798915A (zh) * 2010-04-26 2010-08-11 胡士清 一种利用分段火驱开采提高厚层块状稠油采收率的方法
CN101818637A (zh) * 2010-04-26 2010-09-01 胡士清 一种通过控制火烧注气速度提高厚层块状稠油油藏采收率的方法
CN104060975A (zh) * 2014-06-24 2014-09-24 中国石油大学(北京) 稠油燃烧过程中活化能的预测方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
刘新光: "海上稠油油田火烧油层数值模拟", 《特种油气藏》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112196514A (zh) * 2020-09-30 2021-01-08 中国海洋石油集团有限公司 一种利用热重分析仪测量稠油注空气开发燃料沉积量的方法
CN112196514B (zh) * 2020-09-30 2024-03-29 中国海洋石油集团有限公司 一种利用热重分析仪测量稠油注空气开发燃料沉积量的方法

Similar Documents

Publication Publication Date Title
CN111222271B (zh) 基于基质-裂缝非稳态窜流油藏裂缝数值模拟方法及系统
CN110334431A (zh) 一种低渗透致密气藏单井控制储量计算及剩余气分析方法
CN104989341B (zh) 一种确定低渗透油藏有效驱替注采井距的方法
CN103939066B (zh) 一种一注多采井组定注水量确定油井产液量的方法
CN106126899B (zh) 水驱油藏加密井井位确定方法
CN108959767A (zh) 一种窄河道型气藏不同井型凝析油伤害数值模拟方法
CN112360422B (zh) 一种页岩气藏压裂水平井产量预测方法及系统
CN103114830B (zh) 一种富气驱水/气交替注入方法
CN109611067B (zh) 深层石灰岩储层酸压酸液有效作用距离的数值计算方法
Wang et al. A new high-temperature gel for profile control in heavy oil reservoirs
CN107818188A (zh) 一种针对流体加热油页岩原位开采的模拟预测方法
CN101710021B (zh) 天然气超临界流动状态物性测量装置
CN112444470A (zh) 基于稳定流群孔抽水试验确定地下含水层渗透系数的方法
CN112943176A (zh) 模拟油藏型储气库注采的气油相对渗透率曲线测定方法
CN108331565A (zh) 一种改进型火驱数值模拟动力学模型的建模方法
CN115130274A (zh) 一种稠油油藏相似准则研究方法
CN106321057A (zh) 水平井压裂裂缝井网模型设计方法
CN203145918U (zh) 一种模拟稠油蒸汽吞吐采油的实验装置及其模型系统
CN104747154B (zh) 一种利用驱油效率比提高蒸汽驱剩余油研究精度的方法
Hu et al. Quick assessment to ascertain technical rational well spacing density in artificial water flooding oilfield
CN201714363U (zh) 一种通过控制火烧注气速度开采厚层块状稠油的装置
CN111852461B (zh) 一种模拟致密油藏基岩与裂缝间非稳态窜流的方法
Wang et al. Study on the production law of heavy oil gravity‐drainage‐assisted steam flooding
CN114017002A (zh) 一种测试油页岩自生热原位转化油收率的装置及方法
CN110162906B (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
RJ01 Rejection of invention patent application after publication

Application publication date: 20180727

RJ01 Rejection of invention patent application after publication