CN114201934A - 一种评价低频电场加热页岩油储层时温度场的方法 - Google Patents
一种评价低频电场加热页岩油储层时温度场的方法 Download PDFInfo
- Publication number
- CN114201934A CN114201934A CN202111529967.5A CN202111529967A CN114201934A CN 114201934 A CN114201934 A CN 114201934A CN 202111529967 A CN202111529967 A CN 202111529967A CN 114201934 A CN114201934 A CN 114201934A
- Authority
- CN
- China
- Prior art keywords
- reservoir
- unit
- grid
- time
- electric field
- 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
Links
- 239000003079 shale oil Substances 0.000 title claims abstract description 98
- 230000005684 electric field Effects 0.000 title claims abstract description 48
- 238000000034 method Methods 0.000 title claims abstract description 40
- 239000012530 fluid Substances 0.000 claims abstract description 74
- 239000011435 rock Substances 0.000 claims abstract description 55
- 238000010438 heat treatment Methods 0.000 claims abstract description 54
- 238000004134 energy conservation Methods 0.000 claims abstract description 11
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 claims description 84
- 239000003345 natural gas Substances 0.000 claims description 42
- 239000010779 crude oil Substances 0.000 claims description 40
- 230000008859 change Effects 0.000 claims description 33
- 239000011148 porous material Substances 0.000 claims description 26
- 239000008398 formation water Substances 0.000 claims description 23
- 239000003921 oil Substances 0.000 claims description 22
- 230000008569 process Effects 0.000 claims description 21
- 239000007789 gas Substances 0.000 claims description 20
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 17
- 239000000126 substance Substances 0.000 claims description 6
- 230000009466 transformation Effects 0.000 claims description 5
- 238000006243 chemical reaction Methods 0.000 claims description 3
- 230000008602 contraction Effects 0.000 claims description 3
- 239000002184 metal Substances 0.000 claims description 3
- 229910052751 metal Inorganic materials 0.000 claims description 3
- 238000011065 in-situ storage Methods 0.000 description 6
- 229930195733 hydrocarbon Natural products 0.000 description 4
- 150000002430 hydrocarbons Chemical class 0.000 description 4
- 230000004048 modification Effects 0.000 description 4
- 238000012986 modification Methods 0.000 description 4
- 239000004215 Carbon black (E152) Substances 0.000 description 3
- 230000015572 biosynthetic process Effects 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 238000003860 storage Methods 0.000 description 3
- 230000005672 electromagnetic field Effects 0.000 description 2
- 238000011084 recovery Methods 0.000 description 2
- 230000009467 reduction Effects 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000005336 cracking Methods 0.000 description 1
- 238000005553 drilling Methods 0.000 description 1
- 238000005485 electric heating Methods 0.000 description 1
- 230000005611 electricity Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 239000003027 oil sand Substances 0.000 description 1
- 230000035699 permeability Effects 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
- 238000010248 power generation Methods 0.000 description 1
- 230000005855 radiation Effects 0.000 description 1
- 239000002023 wood Substances 0.000 description 1
Images
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
- G06F30/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B43/00—Methods or apparatus for obtaining oil, gas, water, soluble or meltable materials or a slurry of minerals from wells
- E21B43/16—Enhanced recovery methods for obtaining hydrocarbons
- E21B43/24—Enhanced recovery methods for obtaining hydrocarbons using heat, e.g. steam injection
- E21B43/2401—Enhanced recovery methods for obtaining hydrocarbons using heat, e.g. steam injection by means of electricity
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/08—Fluids
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/08—Thermal analysis or thermal optimisation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
-
- 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
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/30—Assessment of water resources
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Theoretical Computer Science (AREA)
- Fluid Mechanics (AREA)
- Geology (AREA)
- Mining & Mineral Resources (AREA)
- General Physics & Mathematics (AREA)
- Computing Systems (AREA)
- Pure & Applied Mathematics (AREA)
- Geochemistry & Mineralogy (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Mathematical Physics (AREA)
- Algebra (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Environmental & Geological Engineering (AREA)
- Investigating Or Analyzing Materials By The Use Of Electric Means (AREA)
- Investigating Or Analyzing Materials Using Thermal Means (AREA)
Abstract
一种评价低频电场加热页岩油储层时温度场的方法,应用于石油开发领域。该方法步骤:①构建2口水平井井间页岩油储层三维网格模型;②依据泥页岩储层电阻率、加载的电场电压,建立低频电场加热页岩油储层的电场模型,获得电焦耳热模型;③建立一个单位时间内各个储层网格单元岩石骨架热能、流体热能、排出流体热能、热传导热能和电焦耳热能的变化量模型,构建它们之间的能量守恒模型;④确定在一个单位时间结束时各个储层网格单元的温度,以及其它参数值;⑤依据低频电场加热页岩油储层能量守恒模型,循环重复步骤②‑④计算页岩油储层温度场随时间变化特征。
Description
技术领域
本发明涉及石油勘探开发技术领域,一种评价低频电场加热页岩油储层时温度场的方法。
背景技术
针对页岩油开采的原位加热方式主要有热传导(电加热棒、蒸汽循环管、电介质加热等)、对流+传导(注入蒸汽、热有机气体和CO2)和辐射+传导(电场、电磁波和微波等)等三种。由于页岩油储层渗透率低,向页岩油储层注入热流体难度大、成本高。国内外众多学者开展低频电场(属于电介质加热)、微波和高频电磁场加热稠油、油砂、页岩油和页岩气等非常规油气储层等相关研究。壳牌等石油公司已经在若干区块利用加热管大规模加热页岩油储层方式开采页岩油进行了现场试验。低频电场、微波和高频电磁场具有体积式加热特点,能更有效、均匀地加热整个油藏,降低热量损失,提高能量利用率和采收率。原位加热改质技术通过提高储层温度,达到降低原油粘度、热膨胀孔隙流体和岩石骨架体积、转化原油中溶解气为游离气、甚至让泥页岩残留干酪根继续裂解生烃、增加储层流体压力、提高原油采收率的目的。
页岩油储层中具有较高含量的地层水,包括游离水、束缚水和结合水等,而且地层水盐度高达28%。可以将页岩油储层作为电阻元素,相邻两口井在目的层段的金属套管与地面供电设备相连,让低频电流通过页岩油储层并以热能形式损耗,达到加热页岩油储层的目的。利用低频电场加热页岩油储层可以降低页岩油粘度、增加储层流体压力、形成裂缝增加渗流能力、可能促使有机质少量生烃,增加页岩油井产量和降低页岩油开采成本,即对页岩油进行降粘、增压、增渗、增量、增产和降本开采。这样开采页岩油无需大规模水力压裂,也无需专门钻井放置加热管等,还能将产出烃类气体进行发电、继续加热储层。预计该方法能够大幅度降低成本、节约用水、减少对环境的影响。尤其是在西部鄂尔多斯、准噶尔、柴达木、三塘湖和吐哈等盆地,可以充分利用风能和太阳能,减少或杜绝风能、太阳能发电的“弃风、弃电”现象。
低频电场加热储层开采页岩油仍面临如下问题:首先,需要避免井筒周围因电流密度大、温度升温快,造成地层水转化为水蒸气而降低加热效率,需要让加载电场的井筒之间地层相对均匀地升温,且在合理时间段内加热页岩油储层。其次,被低频电场加热的页岩油储层温度能否让干酪根能否继续生烃。因此,温度场是低频电场加热页岩油储层开采页岩油的关键研究内容。
为此,本发明提出一种评价低频电场加热页岩油储层时温度场的方法。该方法可以定量评价低频电场加热页岩油储层过程中温度场随时间变化特征,为低频电场加热页岩油储层对页岩油进行原位改质开采提供重要的技术参数,为突破页岩油开采技术瓶颈提供支持。
发明内容
本发明的目的是:提供一种评价低频电场加热页岩油储层时温度场的方法,实现定量评价低频电场加热页岩油储层过程中温度场随时间变化特征,为低频电场加热页岩油储层对页岩油进行原位改质开采提供重要的技术参数,为突破页岩油开采技术瓶颈提供支持。
本发明采用的技术方案是:评价低频电场加热页岩油储层时温度场的方法,其特征在于:
步骤1:构建2口水平井井间页岩油储层三维网格模型,每个三维网格单元的长、宽、高均为1m,三维网格模型长度等于水平井水平段长度,宽度为2口水平井间距的2倍以上,高度等于页岩油储层厚度,两口井水平段金属套管连接地面电源作为电极,电极电压为U,并确定各网格单元的初始温度T0、初始孔隙度Ф0、初始含油饱和度So, 0、初始含气饱和度Sg, 0、初始含水饱和度Sw, 0、初始原油溶解天然气量Vdis-g, 0,初始横向电阻Rh, 0、初始垂向电阻Rv, 0, 初始岩石骨架比热容ζr, 0、初始各种流体比热容ζp, 0,初始岩石骨架热导率λr, 0、初始各种流体热导率λp, 0,初始岩石骨架热线膨胀系数βr, 0、初始各种流体热膨胀系数βp, 0,其中p=1,2,3为各种流体的编号,确定原油溶解天然气量、横向电阻、垂向电阻、比热容、热导率、热膨胀系数随温度的变化特征,电压的单位为V,温度的单位是℃,孔隙度的单位是%,含油饱和度、含气饱和度、含水饱和度的单位均为%,原油溶解天然气量m3/m3, 横向电阻率、垂向电阻的单位均为Ω,比热容的单位为J/kg/℃,热导率的单位为W/m/℃,热膨胀系数的单位为℃-1;
步骤2:将三维储层网格视为一系列叠置且垂直于水平段井轨迹在水平面投影的连续二维网格面,这些二维网格面厚度为1米而且电场特征相同,首先构建电阻模型:依据各网格单元横向电阻和垂向电阻以及它们随温度的变化特征,相应变化网格单元横向、垂向上的长度,并记录横向、垂向长度变化的倍数,确保伸缩变换之后长、宽、高均为1米的各个网格单元横向、垂向电阻相同,即构建新的均匀电介质网格单元,这样将非均匀电介质储层转换为均匀电介质模型;
然后搭建电势模型:两口相邻的井轨迹与该二维网格面的交点分别为正、负极,在该均匀电介质二维网格面中以正、负电极连线中点为原点、过正、负极直线为x轴构建直角坐标系,其中正电极A点坐标为(-a,0),负电极B点坐标为(a,0),依据均匀电介质中电场特征,该二维面中电场等势线除了y轴之外,其它等势线是以y轴对称的一系列圆,圆点分别为(-a(1+k2)/(1-k2),0)和(a(1+k2)/(1-k2),0),半径为2ka/(1-k2),其中k是等势线上点到正、负电极距离的比值,当k=1时等势线为y轴;
其次构建电流模型:在该均匀电介质二维网格面中每束电流由正极出发,电流线轨迹垂直于各电势线,均经过y轴上各网格单元中点并到达负极,这样一共就有n束电流,n为y轴上所分析的网格单元个数,假设点C(0,c)为y轴正向上任意一个网格单元中心点,每条电流轨迹经过点A(-a,0)、B(a,0)和C(0,c),而且垂直于各等势线,则该电流轨迹为以(0,(c 2-a 2)/(2c))为圆心、半径为(c 2+a 2)/2c的部分圆弧,其函数表达式如下公式,这样就确定了该均匀电介质二维网格面中各电流轨迹;
确定电焦耳热模型:确定了该均匀电介质二维网格面上n束电流轨迹之后,可以获得二维面中各网格单元上电流轨迹的条数和长度,根据各电流轨迹在均匀电介质网格单元内的长度,将每个网格单元贡献电阻分配给各电流轨迹,这样就能获得各束电流轨迹上累计电阻之和,结合电极电压确定各束电流轨迹上的电流强度,利用焦耳定律确定单位时间内各束电流在其轨迹所经过各均匀电介质网格单元上产生的电焦耳热能,单位时间可以是1小时、1分钟、1秒钟,结合电阻模型中储层网格单元伸缩变换特征,将原属于同一个储层网格单元的均匀电介质网格单元电焦耳热能累加,获得各储层网格单元电焦耳热能Qele,电焦耳热能的单位为J;
步骤3:建立一个单位时间内各个储层网格单元岩石骨架热能、流体热能、排出流体热能、热传导热能和电焦耳热能的变化量模型,构建它们之间的能量守恒模型,该能量守恒模型如下式,其中体系传导至上、下围岩的热能通过三维网格模型中位于边界处网格单元的热传导模型进行评价,在时间由ti-1至ti的单位时间内,每个储层网格单元岩石骨架热能变化量Qr等于岩石骨架体积V(1-Ф)、岩石骨架密度ρr、岩石骨架比热容ζr和温度变化量的乘积,每个储层网格单元流体热能变化量Qf等于每种流体体积VSpФ、流体密度ρp、流体比热容ζp和温度变化量的乘积之和,每个储层网格单元排出流体热能变化量Qff等于每种排出流体体积、排出流体密度ρp、排出流体比热容ζp和温度变化量的乘积之和,每个储层网格单元热传导至其四周单元格的热传导热能等于由体积百分比计算的储层热导率与相邻四个单元格距离、温度差的乘积之和,i=1,2,…,m,是大于零的整数,p=1, 2, … , np,是不同种类孔隙流体的编号,单位时间可以是1小时、1分钟、1秒钟,各种热能的单位为J;
Qr+ Qf+ Qff+ Qc= Qele
其中,
式中,Qr是单个储层网格单元从时间ti至ti+1过程中岩石骨架热能变化量,Qf是单个储层网格单元从时间ti至ti+1过程中流体热能变化量,Qff是单个储层网格单元从时间ti至ti+1过程中排出流体热能变化量,Qc是单个储层网格单元从时间ti至ti+1过程中热传导能量,热能变化量的单位均为J,Qele是单个储层网格单元从时间ti至ti+1过程中电焦耳热能,电焦耳热能的单位为J,V是储层网格单元体积,体积的单位为m3,Ф是页岩油储层孔隙度,孔隙度的单位为%,ρr和ρp分别是储层骨架和编号为p的孔隙流体密度,密度的单位为g/cm3,ti和ti+1分别是编号为i和i+1的时间点,i=1,2,…,m,是大于零的整数,t是时间,时间的单位是秒,T是温度,温度的单位是℃,p=1, 2, … , np,是不同种类孔隙流体的编号,ρ p 是编号为p孔隙流体密度,密度的单位为g/cm3,S p 是编号为p孔隙流体饱和度,饱和度的单位为%,ζ r 和ζ p 分别是页岩油储层和编号为p孔隙流体比热容,比热容的单位为J/K/kg,λ r 和λ p 分别是页岩油储层岩石骨架和编号为p孔隙流体导热率,导热率的单位为J/K/m, i=1,2,…,m,是大于零的整数,p=1, 2, … , np,是不同种类孔隙流体的编号;
步骤4:依据步骤3中能量守恒模型,确定在一个单位时间结束时各个储层网格单元的温度,并依据步骤1中各参数随温度的变化特征,确定各储层网格单元在温度变化后的原油溶解天然气量、横向电阻、垂向电阻、储层岩石骨架比热容、地层水比热容、原油比热容、天然气比热容、储层岩石骨架热导率、地层水热导率、原油热导率、天然气热导率、储层岩石骨架热膨胀系数、地层水热膨胀系数、原油热膨胀系数、天然气热膨胀系数,并且按照下列公式进一步确定各储层网格单元在温度变化后的孔隙度、含油饱和度、含气饱和度、含水饱和度;
Фi= Фi-1–3βr, i (Ti– Ti -1)
So, i = So, i -1 βo, i (Ti– Ti -1)/ (So, i -1 βo, i (Ti– T i -1) + Sg, i -1 βg, i (Ti–T i -1) + Sw, i -1 βw, i (Ti– T i -1) +ΔVdis-g, i )
Sg, i = (Sg, i -1 βg, i (Ti– T i -1) +ΔVdis-g, i)/ (So, i -1 βo, i (Ti– T i -1) +Sg, i -1 βg, i (Ti– T i -1) + Sw, i -1 βw, i (Ti - T i -1) +ΔVdis-g, i)
Sw, i = (Sw, i -1 βw, i (Ti– T i -1) +ΔVdis-g, i)/ (So, i -1 βo, i (Ti– T i -1) +Sg, i -1 βg, i (Ti– T i -1) + Sw, i -1 βw, i (Ti– T i -1) +ΔVdis-g, i)
其中,
ΔVdis-g, i = Vdis-g, i-1 – Vdis-g, i = 2.4ρrel-g(P﹒exp (1.77/ρoil– 0.001638Ti-1–1.67))1.205– 2.4ρrel-g(P﹒exp (1.77/ρoil– 0.001638Ti– 1.67))1.205
式中,Фi和Фi-1分别是低频电场加热页岩油储层过程中在时间ti和ti-1时的储层孔隙度,孔隙度的单位为%;βr, i、βo, i、βg, i和βw, i 分别是低频电场加热页岩油储层过程中在时间ti时储层骨架岩石、原油、天然气和地层水的热膨胀系数,热膨胀系数的单位为℃-1;Ti和Ti -1分别是低频电场加热页岩油储层过程中在时间ti和ti-1时的储层温度,储层温度的单位为℃;So, i和So, i -1分别是低频电场加热页岩油储层过程中在时间ti和ti-1时的储层含油饱和度,含油饱和度的单位为%;Sg, i和Sg, i -1分别是低频电场加热页岩油储层过程中在时间ti和ti-1时的储层含气饱和度,含气饱和度的单位为%;Sw, i和Sw, i -1分别是低频电场加热页岩油储层过程中在时间ti和ti-1时的储层含水饱和度,含水饱和度的单位为%; Vdis-g, i和Vdis-g, i-1分别是低频电场加热页岩油储层过程中在时间ti和ti-1时的储层原油溶解天然气量,原油溶解天然气量的单位为m3/m3;ΔVdis-g, i 是低频电场加热页岩油储层过程中在时间由ti-1至ti时储层原油溶解天然气量的变化量,原油溶解天然气量的变化量的单位为m3/m3,i=1,2,…,m,是大于零的整数;
步骤5:随着时间点的编号i每增加一次则循环重复步骤2、步骤3和步骤4,计算获得各页岩油储层网格单元的温度及其随时间变化特征,进而获得低频电场加热页岩油储层过程中储层温度场随时间变化特征。
本发明的有益效果:本发明评价低频电场加热页岩油储层时温度场的方法,该方法易于操作和实现。该方法可定量评价低频电场加热页岩油储层过程中温度场随时间变化特征,为低频电场加热页岩油储层对页岩油进行原位改质开采提供重要的技术参数,为突破页岩油开采技术瓶颈提供支持。
附图说明
图1是本发明评价低频电场加热页岩油储层时温度场的方法流程示意图。
图2是本发明利用低频电场加热页岩油储层原位改质开采页岩油的水平井布置方案图。
图3 是本发明利用低频电场加热页岩油储层时储层网格单元伸缩变换示意图。
图4 是本发明利用低频电场加热页岩油储层时均匀电介质二维剖面电场示意图。
图5是本发明利用低频电场加热页岩油储层46800秒后一个垂直于水平段井轨迹在水平面投影的二维网格面的温度等值线结果。
具体实施方式
实施例1:如图1所述,一种评价低频电场加热页岩油储层时温度场的方法,含有以下步骤:
步骤1:构建2口水平井井间页岩油储层三维网格模型(图2),这2口井水平段互相平行,2口井水平段距离为60米,两口井水平段金属套管连接地面电源作为电极,电压为1000V,每个三维网格单元的长、宽、高均为1m。三维网格模型长度等于水平井水平段长度,宽度为2口水平井间距的2倍以上,高度等于页岩油储层厚度。确定各网格单元的初始温度为90℃、初始孔隙度为2%、初始含油饱和度60%、初始含气饱和度0%、初始含水饱和度40%、初始原油溶解天然气量95.84m3/m3,初始横向电阻20Ω、初始垂向电阻40Ω, 初始岩石骨架比热容值为1127 J/kg/℃、初始地层水比热容为4198.13 J/kg/℃,初始原油比热容为2100J/kg/℃,初始天然气比热容为2087.78J/kg/℃,初始岩石骨架热导率为1.8093 W/m/℃,初始地层水热导率为0.6663 W/m/℃,初始原油热导率为0.13 W/m/℃,初始天然气热导率为0.0822W/m/℃,初始岩石骨架热膨胀系数为0.000007℃-1,初始地层水热膨胀系数为0.000512℃-1,初始原油热膨胀系数为0.000784℃-1,初始天然气热膨胀系数为0.002754℃-1。确定原油溶解天然气量、横向电阻、垂向电阻、比热容、热导率、热膨胀系数随温度T的变化特征。其中,原油比热容为定值2100J/kg/℃,原油热导率为定值0.13 W/m/℃,储层岩石骨架热膨胀系数为定值0.000007℃-1,原油热膨胀系数为定值0.000784℃-1。
原油溶解天然气量随温度T的变化特征为:
Vdis-g = 2.4ρrel-g(P﹒exp (1.77/ρoil– 0.001638T – 1.67))1.205
储层网格单元格横向电阻、垂向电阻随温度T的变化特征分别为:
Rh= 7.2299×10-9 T4– 9.4529×10-6 T3 + 4.5609 ×10-3 T2– 0.9659 T+76.9874
Rv= 1.446×10-8 T4– 1.8906 ×10-5 T3 + 9.1 218×10-3 T2– 1.9318 T+153.9748
储层岩石骨架、地层水和天然气的比热容随温度T的变化特征分别为:
ζr = – 0.7 T + 1190
ζw = 2.0238 ×10-6 T4– 1.056 ×10-3 T3 + 0.1784 T2 – 8.7463T4177.6075
ζg = 3.7731 T + 1748.2
储层岩石骨架、地层水和天然气的热导率随温度T的变化特征分别为:
λ r = –3×10-9 T3 + 5×10-5 T2– 0.003 T + 1.6765
λ w = –5×10-8 T3 + 2 ×10-5 T2– 0.0038T+ 0.8827
λ g = –1×10-9 T3 + 10-6 T2– 0.0004 T + 0.1108
地层水和天然气的热膨胀系数随温度T的变化特征分别为:
β w = 2 ×10-10 T3– 4×10-8 T2 + 8×10-6 T – 0.00003
β g = ((T + 273.15) / (T+ 272.15) – 1)
式中,Vdis-g是原油溶解天然气量,单位为m3/m3;ρrel-g是天然气相对密度,取值为0.6;P是储层孔隙流体压力,取值为25MPa;ρoil是原油的密度,取值为0.85g/cm3;T是储层的温度,单位为℃;Rh是储层网格单元格横向电阻,单位是Ω;Rv储层网格单元格垂向电阻,单位是Ω;ζr是储层岩石骨架的比热容,单位为J/kg/℃;ζw是地层水的比热容,单位为J/kg/℃;ζg是天然气的比热容,单位为J/kg/℃;λ r 是储层岩石骨架的热导率,单位为W/m/℃;λ w 是地层水的热导率,单位为W/m/℃;λ g 是天然气的热导率,单位为W/m/℃; β w 是地层水的热膨胀系数,单位为℃-1;β g 是天然气的热膨胀系数,单位为℃-1。
步骤2:将三维储层网格视为一系列叠置且垂直于水平段井轨迹在水平面的投影的连续二维网格面,这些二维网格面厚度为1米而且电场特征相同。以其中一个二维储层网格面为例开展电场模型的搭建。首先构建电阻模型:如图3所示,依据各网格单元横向电阻和垂向电阻以及它们随温度的变化特征,相应变化网格单元横向、垂向的长度,并记录横向、垂向长度变化的倍数,确保伸缩变换之后长、宽、高均为1米的各个网格单元横向、垂向电阻相同,即构建新的均匀电介质网格单元,这样将该二维非均匀电介质储层网格面转换为均匀电介质模型。
然后搭建电势模型:如图4所示,两口相邻的井轨迹与该二维网格面的交点分别为正、负极。在该均匀电介质二维网格面中以正、负电极连线中点为原点、过正、负极直线为x轴构建直角坐标系,其中正电极A点坐标为(-a,0),负电极B点坐标为(a,0)。依据均匀电介质中电场特征,该二维面中电场等势线除了y轴之外,其它等势线是以y轴对称的一系列圆,圆点分别为(-a(1+k2)/(1-k2),0)和(a(1+k2)/(1-k2),0),半径为2ka/(1-k2)。其中k是等势线上点到正、负电极距离的比值,当k=1时等势线为y轴。
其次构建电流模型:如图4所示,在该均匀电介质二维网格面中每束电流由正极出发,电流线轨迹垂直于各电势线,均经过y轴上各网格单元中点并到达负极。这样一共就有n束电流,n为y轴上所分析的网格单元个数,假设点C(0,c)为y轴正向上任意一个网格单元中心点,每条电流轨迹经过点A(-a,0)、B(a,0)和C(0,c),而且垂直于各等势线,则该电流轨迹为以(0, (c 2-a 2)/(2c))为圆心、半径为(c 2+a 2)/2c的部分圆弧,其函数表达式如下公式,这样就确定了该均匀电介质二维网格面中各电流轨迹。
确定电焦耳热模型:确定了该均匀电介质二维网格面上n束电流轨迹之后,可以获得二维面中各网格单元上电流轨迹的条数和长度,根据各电流轨迹在均匀电介质网格单元内的长度,将每个网格单元贡献电阻分配给各电流轨迹。这样就能获得各束电流轨迹上累计电阻之和,结合电极电压为1000V,确定各束电流轨迹上的电流强度,利用焦耳定律确定单位时间内各束电流在其轨迹所经过各均匀电介质网格单元上产生的电焦耳热能,单位时间是1秒钟。结合电阻模型中储层网格单元伸缩变换特征,将原属于同一个储层网格单元的均匀电介质网格单元电焦耳热能累加,获得各储层网格单元电焦耳热能Qele,电焦耳热能的单位为J。
步骤3:建立一个单位时间内各个储层网格单元岩石骨架热能、流体热能、排出流体热能、热传导热能和电焦耳热能的变化量模型,构建它们之间的能量守恒模型。该能量守恒模型如下式,其中体系传导至上、下围岩的热能通过三维网格模型中位于边界处网格单元的热传导模型进行评价。在时间由ti-1至ti的单位时间内,每个储层网格单元岩石骨架热能变化量Qr等于岩石骨架体积V(1-Ф)、岩石骨架密度ρr、岩石骨架比热容ζr和温度变化量的乘积;每个储层网格单元流体热能变化量Qf等于每种流体体积VSpФ、流体密度ρp、流体比热容ζp和温度变化量的乘积之和;每个储层网格单元排出流体热能变化量Qff等于每种排出流体体积、排出流体密度ρp、排出流体比热容ζp和温度变化量的乘积之和;每个储层网格单元热传导至其四周单元格的热传导热能等于由体积百分比计算的储层热导率与相邻四个单元格距离、温度差的乘积之和;i=1,2,…,m,是大于零的整数;p=1, 2, 3,分别是地层孔隙水、原油和天然气的编号,单位时间是1秒钟,各种热能的单位为J。
Qr+ Qf+ Qff+ Qc = Qele
其中,
式中,Qr是单个储层网格单元从时间ti至ti+1过程中岩石骨架热能变化量,Qf是单个储层网格单元从时间ti至ti+1过程中流体热能变化量,Qff是单个储层网格单元从时间ti至ti+1过程中排出流体热能变化量,Qc是单个储层网格单元从时间ti至ti+1过程中热传导能量,热能变化量的单位均为J,Qele是单个储层网格单元从时间ti至ti+1过程中电焦耳热能,电焦耳热能的单位为J,V是储层网格单元体积,体积的单位为m3,Ф是页岩油储层孔隙度,孔隙度的单位为%,ρr和ρp分别是储层骨架和编号为p的孔隙流体密度,密度的单位为g/cm3,ti和ti+1分别是编号为i和i+1的时间点,i=1,2,…,m,是大于零的整数,t是时间,时间的单位是秒,T是温度,温度的单位是℃,p=1, 2, … , np,是不同种类孔隙流体的编号,ρ p 是编号为p孔隙流体密度,密度的单位为g/cm3,S p 是编号为p孔隙流体饱和度,饱和度的单位为%,ζ r 和ζ p 分别是页岩油储层和编号为p孔隙流体比热容,比热容的单位为J/K/kg,λ r 和λ p 分别是页岩油储层岩石骨架和编号为p孔隙流体导热率,导热率的单位为J/K/m, i=1,2,…,m,是大于零的整数, p=1, 2, … , np,是不同种类孔隙流体的编号。
步骤4:依据步骤3中能量守恒模型,确定在一个单位时间结束时各个储层网格单元的温度,并依据步骤1中各参数随温度的变化特征,确定各储层网格单元在温度变化后的原油溶解天然气量、横向电阻、垂向电阻、储层岩石骨架比热容、地层水比热容、原油比热容、天然气比热容、储层岩石骨架热导率、地层水热导率、原油热导率、天然气热导率、储层岩石骨架热膨胀系数、地层水热膨胀系数、原油热膨胀系数、天然气热膨胀系数,并且按照下列公式进一步确定各储层网格单元在温度变化后的孔隙度、含油饱和度、含气饱和度、含水饱和度。
Фi= Фi-1–3βr, i (Ti– Ti -1)
So, i = So, i -1 βo, i (Ti– Ti -1)/ (So, i -1 βo, i (Ti– T i -1) + Sg, i -1 βg, i (Ti–T i -1) + Sw, i -1 βw, i (Ti– T i -1) +ΔVdis-g, i )
Sg, i = (Sg, i -1 βg, i (Ti– T i -1) +ΔVdis-g, i)/ (So, i -1 βo, i (Ti– T i -1) +Sg, i -1 βg, i (Ti– T i -1) + Sw, i -1 βw, i (Ti - T i -1) +ΔVdis-g, i)
Sw, i = (Sw, i -1 βw, i (Ti– T i -1) +ΔVdis-g, i)/ (So, i -1 βo, i (Ti– T i -1) +Sg, i -1 βg, i (Ti– T i -1) + Sw, i -1 βw, i (Ti– T i -1) +ΔVdis-g, i)
其中,
ΔVdis-g, i = Vdis-g, i-1 – Vdis-g, i = 2.4ρrel-g(P﹒exp (1.77/ρoil– 0.001638Ti-1–1.67))1.205– 2.4ρrel-g(P﹒exp (1.77/ρoil– 0.001638Ti– 1.67))1.205
式中,Фi和Фi-1分别是低频电场加热页岩油储层过程中在时间ti和ti-1时的储层孔隙度,孔隙度的单位为%;βr, i、βo, i、βg, i和βw, i 分别是低频电场加热页岩油储层过程中在时间ti时储层骨架岩石、原油、天然气和地层水的热膨胀系数,热膨胀系数的单位为℃-1;Ti和Ti -1分别是低频电场加热页岩油储层过程中在时间ti和ti-1时的储层温度,储层温度的单位为℃;So, i和So, i -1分别是低频电场加热页岩油储层过程中在时间ti和ti-1时的储层含油饱和度,含油饱和度的单位为%;Sg, i和Sg, i -1分别是低频电场加热页岩油储层过程中在时间ti和ti-1时的储层含气饱和度,含气饱和度的单位为%;Sw, i和Sw, i -1分别是低频电场加热页岩油储层过程中在时间ti和ti-1时的储层含水饱和度,含水饱和度的单位为%; Vdis-g, i和Vdis-g, i-1分别是低频电场加热页岩油储层过程中在时间ti和ti-1时的储层原油溶解天然气量,原油溶解天然气量的单位为m3/m3;ΔVdis-g, i 是低频电场加热页岩油储层过程中在时间由ti-1至ti时储层原油溶解天然气量的变化量,原油溶解天然气量的变化量的单位为m3/m3;i=1,2,…,m,是大于零的整数。
步骤5:随着时间点的编号i每增加一次则循环重复步骤2-步骤4,计算获得各页岩油储层网格单元的温度及其随时间变化特征,进而获得低频电场加热页岩油储层过程中储层温度场随时间变化特征。如图5所示,展示了低频电场加热页岩油储层46800秒,即13小时后的一个垂直于水平段井轨迹在水平面投影的二维网格面的温度等值线特征。
Claims (1)
1.一种评价低频电场加热页岩油储层时温度场的方法,其特征在于:
步骤1:构建2口水平井井间页岩油储层三维网格模型,每个三维网格单元的长、宽、高均为1m,三维网格模型长度等于水平井水平段长度,宽度为2口水平井间距的2倍以上,高度等于页岩油储层厚度,两口井水平段金属套管连接地面电源作为电极,电极电压为U,并确定各网格单元的初始温度T0、初始孔隙度Ф0、初始含油饱和度So, 0、初始含气饱和度Sg, 0、初始含水饱和度Sw, 0、初始原油溶解天然气量Vdis-g, 0,初始横向电阻Rh, 0、初始垂向电阻Rv, 0, 初始岩石骨架比热容ζr, 0、初始各种流体比热容ζp, 0,初始岩石骨架热导率λr, 0、初始各种流体热导率λp, 0,初始岩石骨架热线膨胀系数βr, 0、初始各种流体热膨胀系数βp, 0,其中p=1,2,3为各种流体的编号,确定原油溶解天然气量、横向电阻、垂向电阻、比热容、热导率、热膨胀系数随温度的变化特征,电压的单位为V,温度的单位是℃,孔隙度的单位是%,含油饱和度、含气饱和度、含水饱和度的单位均为%,原油溶解天然气量m3/m3, 横向电阻率、垂向电阻的单位均为Ω,比热容的单位为J/kg/℃,热导率的单位为W/m/℃,热膨胀系数的单位为℃-1;
步骤2:将三维储层网格视为一系列叠置且垂直于水平段井轨迹在水平面投影的连续二维网格面,这些二维网格面厚度为1米而且电场特征相同,首先构建电阻模型:依据各网格单元横向电阻和垂向电阻以及它们随温度的变化特征,相应变化网格单元横向、垂向上的长度,并记录横向、垂向长度变化的倍数,确保伸缩变换之后长、宽、高均为1米的各个网格单元横向、垂向电阻相同,即构建新的均匀电介质网格单元,这样将非均匀电介质储层转换为均匀电介质模型;
然后搭建电势模型:两口相邻的井轨迹与该二维网格面的交点分别为正、负极,在该均匀电介质二维网格面中以正、负电极连线中点为原点、过正、负极直线为x轴构建直角坐标系,其中正电极A点坐标为(-a,0),负电极B点坐标为(a,0),依据均匀电介质中电场特征,该二维面中电场等势线除了y轴之外,其它等势线是以y轴对称的一系列圆,圆点分别为(-a(1+k2)/(1-k2),0)和(a(1+k2)/(1-k2),0),半径为2ka/(1-k2),其中k是等势线上点到正、负电极距离的比值,当k=1时等势线为y轴;
其次构建电流模型:在该均匀电介质二维网格面中每束电流由正极出发,电流线轨迹垂直于各电势线,均经过y轴上各网格单元中点并到达负极,这样一共就有n束电流,n为y轴上所分析的网格单元个数,假设点C(0,c)为y轴正向上任意一个网格单元中心点,每条电流轨迹经过点A(-a,0)、B(a,0)和C(0,c),而且垂直于各等势线,则该电流轨迹为以(0, (c 2-a 2)/(2c))为圆心、半径为(c 2+a 2)/2c的部分圆弧,其函数表达式如下公式,这样就确定了该均匀电介质二维网格面中各电流轨迹;
确定电焦耳热模型:确定了该均匀电介质二维网格面上n束电流轨迹之后,可以获得二维面中各网格单元上电流轨迹的条数和长度,根据各电流轨迹在均匀电介质网格单元内的长度,将每个网格单元贡献电阻分配给各电流轨迹,这样就能获得各束电流轨迹上累计电阻之和,结合电极电压确定各束电流轨迹上的电流强度,利用焦耳定律确定单位时间内各束电流在其轨迹所经过各均匀电介质网格单元上产生的电焦耳热能,单位时间可以是1小时、1分钟、1秒钟,结合电阻模型中储层网格单元伸缩变换特征,将原属于同一个储层网格单元的均匀电介质网格单元电焦耳热能累加,获得各储层网格单元电焦耳热能Qele,电焦耳热能的单位为J;
步骤3:建立一个单位时间内各个储层网格单元岩石骨架热能、流体热能、排出流体热能、热传导热能和电焦耳热能的变化量模型,构建它们之间的能量守恒模型,该能量守恒模型如下式,其中体系传导至上、下围岩的热能通过三维网格模型中位于边界处网格单元的热传导模型进行评价,在时间由ti-1至ti的单位时间内,每个储层网格单元岩石骨架热能变化量Qr等于岩石骨架体积V(1-Ф)、岩石骨架密度ρr、岩石骨架比热容ζr和温度变化量的乘积,每个储层网格单元流体热能变化量Qf等于每种流体体积VSpФ、流体密度ρp、流体比热容ζp和温度变化量的乘积之和,每个储层网格单元排出流体热能变化量Qff等于每种排出流体体积、排出流体密度ρp、排出流体比热容ζp和温度变化量的乘积之和,每个储层网格单元热传导至其四周单元格的热传导热能等于由体积百分比计算的储层热导率与相邻四个单元格距离、温度差的乘积之和,i=1,2,…,m,是大于零的整数,p=1, 2, … , np,是不同种类孔隙流体的编号,单位时间可以是1小时、1分钟、1秒钟,各种热能的单位为J;
Qr+ Qf+ Qff+ Qc = Qele
其中,
式中,Qr是单个储层网格单元从时间ti至ti+1过程中岩石骨架热能变化量,Qf是单个储层网格单元从时间ti至ti+1过程中流体热能变化量,Qff是单个储层网格单元从时间ti至ti+1过程中排出流体热能变化量,Qc是单个储层网格单元从时间ti至ti+1过程中热传导能量,热能变化量的单位均为J,Qele是单个储层网格单元从时间ti至ti+1过程中电焦耳热能,电焦耳热能的单位为J,V是储层网格单元体积,体积的单位为m3,Ф是页岩油储层孔隙度,孔隙度的单位为%,ρr和ρp分别是储层骨架和编号为p的孔隙流体密度,密度的单位为g/cm3,ti和ti+1分别是编号为i和i+1的时间点,i=1,2,…,m,是大于零的整数,t是时间,时间的单位是秒,T是温度,温度的单位是℃,p=1, 2, … , np,是不同种类孔隙流体的编号,ρ p 是编号为p孔隙流体密度,密度的单位为g/cm3,S p 是编号为p孔隙流体饱和度,饱和度的单位为%,ζ r 和ζ p 分别是页岩油储层和编号为p孔隙流体比热容,比热容的单位为J/K/kg,λ r 和λ p 分别是页岩油储层岩石骨架和编号为p孔隙流体导热率,导热率的单位为J/K/m, i=1,2,…,m,是大于零的整数,p=1, 2, … , np,是不同种类孔隙流体的编号;
步骤4:依据步骤3中能量守恒模型,确定在一个单位时间结束时各个储层网格单元的温度,并依据步骤1中各参数随温度的变化特征,确定各储层网格单元在温度变化后的原油溶解天然气量、横向电阻、垂向电阻、储层岩石骨架比热容、地层水比热容、原油比热容、天然气比热容、储层岩石骨架热导率、地层水热导率、原油热导率、天然气热导率、储层岩石骨架热膨胀系数、地层水热膨胀系数、原油热膨胀系数、天然气热膨胀系数,并且按照下列公式进一步确定各储层网格单元在温度变化后的孔隙度、含油饱和度、含气饱和度、含水饱和度;
Фi= Фi-1–3βr, i (Ti – Ti -1)
So, i = So, i -1 βo, i (Ti – Ti -1)/ (So, i -1 βo, i (Ti – T i -1) + Sg, i -1 βg, i (Ti –T i -1) + Sw, i -1 βw, i (Ti – T i -1) +ΔVdis-g, i )
Sg, i = (Sg, i -1 βg, i (Ti – T i -1) +ΔVdis-g, i)/ (So, i -1 βo, i (Ti – T i -1) +Sg, i -1 βg, i (Ti – T i -1) + Sw, i -1 βw, i (Ti - T i -1) +ΔVdis-g, i)
Sw, i = (Sw, i -1 βw, i (Ti – T i -1) +ΔVdis-g, i)/ (So, i -1 βo, i (Ti – T i -1) +Sg, i -1 βg, i (Ti – T i -1) + Sw, i -1 βw, i (Ti – T i -1) +ΔVdis-g, i)
其中,
ΔVdis-g, i = Vdis-g, i-1 – Vdis-g, i = 2.4ρrel-g(P﹒exp (1.77/ρoil – 0.001638Ti-1 –1.67))1.205 – 2.4ρrel-g(P﹒exp (1.77/ρoil – 0.001638Ti – 1.67))1.205
式中,Фi和Фi-1分别是低频电场加热页岩油储层过程中在时间ti和ti-1时的储层孔隙度,孔隙度的单位为%;βr, i、βo, i、βg, i和βw, i 分别是低频电场加热页岩油储层过程中在时间ti时储层骨架岩石、原油、天然气和地层水的热膨胀系数,热膨胀系数的单位为℃-1;Ti和Ti -1分别是低频电场加热页岩油储层过程中在时间ti和ti-1时的储层温度,储层温度的单位为℃;So, i和So, i -1分别是低频电场加热页岩油储层过程中在时间ti和ti-1时的储层含油饱和度,含油饱和度的单位为%;Sg, i和Sg, i -1分别是低频电场加热页岩油储层过程中在时间ti和ti-1时的储层含气饱和度,含气饱和度的单位为%;Sw, i和Sw, i -1分别是低频电场加热页岩油储层过程中在时间ti和ti-1时的储层含水饱和度,含水饱和度的单位为%; Vdis-g, i和Vdis-g, i-1分别是低频电场加热页岩油储层过程中在时间ti和ti-1时的储层原油溶解天然气量,原油溶解天然气量的单位为m3/m3;ΔVdis-g, i 是低频电场加热页岩油储层过程中在时间由ti-1至ti时储层原油溶解天然气量的变化量,原油溶解天然气量的变化量的单位为m3/m3,i=1,2,…,m,是大于零的整数;
步骤5:随着时间点的编号i每增加一次则循环重复步骤2、步骤3和步骤4,计算获得各页岩油储层网格单元的温度及其随时间变化特征,进而获得低频电场加热页岩油储层过程中储层温度场随时间变化特征。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111529967.5A CN114201934B (zh) | 2021-12-15 | 2021-12-15 | 一种评价低频电场加热页岩油储层时温度场的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111529967.5A CN114201934B (zh) | 2021-12-15 | 2021-12-15 | 一种评价低频电场加热页岩油储层时温度场的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114201934A true CN114201934A (zh) | 2022-03-18 |
CN114201934B CN114201934B (zh) | 2022-08-26 |
Family
ID=80653688
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111529967.5A Active CN114201934B (zh) | 2021-12-15 | 2021-12-15 | 一种评价低频电场加热页岩油储层时温度场的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114201934B (zh) |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111794733A (zh) * | 2020-08-08 | 2020-10-20 | 东北石油大学 | 一种原位电加热开采页岩油储层温度场测量方法 |
CN112966422A (zh) * | 2021-03-18 | 2021-06-15 | 中国石油大学(华东) | 一种适用于页岩气藏热力开采的流动模拟方法 |
-
2021
- 2021-12-15 CN CN202111529967.5A patent/CN114201934B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111794733A (zh) * | 2020-08-08 | 2020-10-20 | 东北石油大学 | 一种原位电加热开采页岩油储层温度场测量方法 |
CN112966422A (zh) * | 2021-03-18 | 2021-06-15 | 中国石油大学(华东) | 一种适用于页岩气藏热力开采的流动模拟方法 |
Non-Patent Citations (1)
Title |
---|
王国营等: "原位注蒸汽开采横观各向同性油页岩储层的数值研究", 《太原理工大学学报》 * |
Also Published As
Publication number | Publication date |
---|---|
CN114201934B (zh) | 2022-08-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Wu et al. | Prospect of marine natural gas hydrate stimulation theory and technology system | |
Sresty et al. | Recovery of bitumen from tar sand deposits with the radio frequency process | |
US5042579A (en) | Method and apparatus for producing tar sand deposits containing conductive layers | |
Song et al. | Production performance of oil shale in-situ conversion with multilateral wells | |
Hascakir | Introduction to thermal Enhanced Oil Recovery (EOR) special issue | |
CN112084695A (zh) | 一种低频电场加热开采天然气水合物藏的数值模拟方法 | |
McGee et al. | The mechanisms of electrical heating for the recovery of bitumen from oil sands | |
Hao et al. | Temperature distribution simulation and optimization design of electric heater for in-situ oil shale heating | |
Wang et al. | Three-dimensional numerical simulation of enhancing shale gas desorption by electrical heating with horizontal wells | |
Bridges et al. | The IITRI in situ RF fuel recovery process | |
Maes et al. | Modelling in-situ upgrading of heavy oil using operator splitting method | |
Wang et al. | Numerical simulation of enhancing shale gas recovery using electrical resistance heating method | |
Ghannadi et al. | Overview of performance and analytical modeling techniques for electromagnetic heating and applications to steam-assisted-gravity-drainage process startup | |
Vinsome et al. | Electrical heating | |
CA2882069A1 (en) | Wide frequency range modeling of electromagnetic heating for heavy oil recovery | |
Zhong et al. | Feasibility study on produce heavy oil by gas and electrical heating assisted gravity drainage | |
CN114201934B (zh) | 一种评价低频电场加热页岩油储层时温度场的方法 | |
Tao et al. | Experimental Study on Well Placement Optimization for Steam‐Assisted Gravity Drainage to Enhance Recovery of Thin Layer Oil Sand Reservoirs | |
Jun et al. | Shale reservoir conductive mechanism simulation based on percolation network | |
CN205225594U (zh) | 地热中低温发电系统 | |
CN102834585A (zh) | 地下地层的低温感应加热 | |
Tian et al. | Modeling of geothermal power generation from abandoned oil wells using in-situ combustion technology | |
Pan et al. | Review of electric heating technology of oil shale and its simulation | |
Alajmi et al. | Investigating the performance of hot water injection in geostatistically generated permeable media | |
Wang et al. | Study on the production law of heavy oil gravity‐drainage‐assisted steam flooding |
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 |