CN111191389B - 一种多元热流体热采油藏数值模拟方法 - Google Patents
一种多元热流体热采油藏数值模拟方法 Download PDFInfo
- Publication number
- CN111191389B CN111191389B CN202010176052.XA CN202010176052A CN111191389B CN 111191389 B CN111191389 B CN 111191389B CN 202010176052 A CN202010176052 A CN 202010176052A CN 111191389 B CN111191389 B CN 111191389B
- Authority
- CN
- China
- Prior art keywords
- equation
- thermal fluid
- fluid
- reservoir
- temperature
- 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
Links
Images
Classifications
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH DRILLING; MINING
- E21B—EARTH DRILLING, e.g. DEEP 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/20—Displacing by water
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH DRILLING; MINING
- E21B—EARTH DRILLING, e.g. DEEP 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
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)
Abstract
本发明属于稠油油藏热采技术领域,具体公开了一种多元热流体热采油藏数值模拟方法,步骤一:进行PVT实验,利用测试得到的数据得到多元热流体PR‑EOS状态方程;基于该状态方程求出多元热流体‑稠油体系的相图;步骤二:进行储层伤害评价实验,建立多元热流体储层伤害数学模型;步骤三:利用步骤一、步骤二得到的新型多元热流体PR‑EOS状态方程和储层伤害数学模型,建立多元热流体数值模拟方法,并进行编程实现,得到多元热流体数值模拟器;通过数值模拟器得到多个在不同条件下多元热流体热采的施工方案,进行比较,得到最优的方案。使用本方案,能有效解决现有多元热流体热采油藏数值模拟方法无法准确描述多元热流体相态特征和多元热流体对储层伤害的问题。
Description
技术领域
本发明属于稠油油藏热采技术领域,尤其涉及一种多元热流体热采油藏数值模拟方法。
背景技术
渤海特殊稠油资源量大,三级地质储量近7.8亿吨,但对于粘度大于350mPa.s的特殊稠油,冷采产能低,甚至不流动。采用多元热流体吞吐或蒸汽吞吐开发,流体粘度大幅度降低,产能大幅度提高,开发效果明显得到改善,因此热采开发势在必行。多元热流体吞吐主要通过注入蒸汽和烟道气体(CO2和N2),然后关井焖井,数日后开井排液采油的方式,来达到降低稠油粘度、提高稠油流动性及产能的目的,该方法目前已成为渤海稠油油田主要热采试验方式之一。
热采数值模拟可以研究合理的开发方案,选择最佳的开采参数,稠油热采模拟是海上稠油热采开发方案编制的必要手段。我国2012年才正式将多元热流体吞吐的开采方式应用在稠油开采中,稠油热采数值模拟的实操经验还相对比较欠缺。除了实操经验不足外,本次开采的稠油对象跟以前还存在很大区别。我们以前开采的稠油主要针对的是浅层稠油(开采深度在1000米以内),而且稠油是低温、低压的状态(相态是气态),而渤海的稠油是深层稠油,除了埋深深外,稠油还是高温、高压的状态(相态是液态),因此开采本次的深层稠油的热采数值模拟不能直接借鉴现有的热采数值模拟方式,需要重新调整。
除了上述提到的问题外,现有的数值模拟考虑的因素不多、不够全面,导致模拟精度较低。如目前针对单一组分的增油机理研究成果较多,但是多元热流体以固定配比混合后,各组分之间的耦合作用以及协同作用仍不清楚;N2在稠油中的溶解度与温度呈正相关,CO2在稠油中的溶解度与温度呈负相关关系,因此导致CO2+N2混合气体在稠油中的溶解度呈现非单调性,但目前没有研究针对CO2和N2的混合气体与稠油作用系数BIP进行调整,分析其混合气体对稠油高压物性的影响;向储层中注入多元热流体时,多元热流体会与储层岩石发生复杂的物理化学反应,孔隙中流体压力会大幅度增加,油藏温度会发生剧烈的变化,岩石颗粒也会随之发生位移,但现有的数值模拟方法没有将这些由于多元热流体的注入引起变化的参数考虑在内。
虽然采用多元热流体吞吐的方式能够极大的提高产能,但是将该方式用于渤海稠油的开采还存在很多有待解决的难题。
发明内容
本发明的目的在于提供一种多元热流体热采油藏数值模拟方法,以解决现有多元热流体热采油藏数值模拟方法模拟精度较低的问题。
为了达到上述目的,本发明的基础方案为:一种多元热流体热采油藏数值模拟方法,包括以下步骤:
步骤一:进行PVT实验,利用测试得到的数据得到新型多元热流体PR-EOS状态方程;基于新型多元热流体PR-EOS状态方程求出多元热流体-稠油体系的相图;
步骤二:进行储层伤害评价实验,描述多元热流体对储层渗透率的伤害程度,建立多元热流体储层伤害数学模型;
步骤三:利用步骤一、步骤二得到的新型多元热流体PR-EOS状态方程和储层伤害数学模型,建立多元热流体数值模拟方法,并进行编程实现,得到了多元热流体数值模拟器;通过数值模拟器得到多个在不同条件下多元热流体热采的施工方案,比较不同的施工方案,得到最优的多元热流体热采施工方案。
本基础方案的有益效果在于:本方案中的多元热流体热采油藏数值模拟方法,充分考虑了多元热流体的相态特征以及多元热流体对储层性质的影响,修正了现有的多元热流体热采油藏数值模拟模型,能够极大的提高现有多元热流体热采油藏数值模拟方法的精度。
进一步,步骤一中,PVT实验的具体操作步骤为:
(1)准备PVT筒,将PVT筒清洗干净并抽真空,将一定量的多元热流体和稠油样品以80℃恒定温度转入PVT腔中;
(2)保持混合物在压力大于200bar的条件下注入PVT筒,且初始状态时为液态单相;
(3)将混合物在80℃的条件下搅拌均匀并保持12h以上;
(4)降低压力,记录在200bar,150bar,100bar,50bar四个压力点的体积,保证每一个体积记录点相态平衡;
(5)绘制PV曲线,拐点为饱和压力;
(6)在150℃和280℃的条件下,重复实验。
进一步,步骤一中PR-EOS模型为:
其中,
其中P为压力,T为温度,V为摩尔体积,Pc为临界压力,Tc为临界温度,R为通用气体常数,Tr为对比温度;
Soave函数为,
如果偏心因子小于0.49,
m=0.37464+1.54226ω-0.26992ω2
如果偏心因子大于0.49,
m=0.3796+1.485ω-0.1644ω2+0.01667ω3
对于混合物,应用范德瓦尔斯混合规则计算参数a和b如下:
其中,kij是两个分量之间的BIP。
进一步,BIP影响PR-EOS模型的准确性,采用临界体积法计算BIP:
其中kij是第i个分量和第j个分量之间的BIP,vci和vcj分别为第i分量和第j分量的临界体积,θ为该临界体积法的指数常数。
进一步,BIP影响PR-EOS模型的准确性,采用临界温度法计算BIP:
其中kij是第i个分量和第j个分量之间的BIP,Tci和Tcj分别为第i分量和第j分量的临界温度,σ为指数常数。
进一步,步骤二中,储层伤害评价实验的具体操作步骤为:
(1)选用高温高压岩心驱替装置进行实验,高温高压岩心驱替装置包括蒸汽发生器和恒温箱;
(2)利用地层水作为驱替介质,测试岩心的水测渗透率;
(3)设置多元热流体注入速度是1mL/min;调整蒸汽发生器的温度与方案设计温度相同,调整恒温箱的温度与地层温度相同,当蒸汽发生器温度和恒温箱温度达到设计要求后,注入24PV多元热流体,记录注入多元热流体过程中的压力数值;
(4)当注入方案设计的多元热流体的PV数后,注入地层水,第二次测试岩心的水测渗透率;
(5)变化多元热流体的温度,重复上述实验步骤,多元热流体的温度分别是150℃,200℃,250℃和300℃。
进一步,步骤二中,建立的多元热流体储层伤害数学模型为:
进一步,步骤三中,得到数值模拟器的过程为:
数值模拟模型由4大类方程构成,分别是质量守恒方程、能量守恒方程、约束方程和边界条件方程;其中,质量守恒方程的基础是质量守恒原理,对于任意组分,质量守恒原理表述为:
ΔM=M1-M2+M3
根据质量守恒原理,水组分的质量守恒方程是:
气组分的质量守恒方程是:
轻质油组分的质量守恒方程是:
重质油组分的质量守恒方程是:
采用显式的时变渗透率处理方法,根据第t-1时间步的绝对渗透率,PV数和温度,计算第t时间步绝对渗透率,具体计算方法如下:
K(t)=(beta1×PV5+beta2×PV4+beta3×PV3+beta4×PV2+beta5×PV+beta6×T3+beta7×T2+beta8×T+beta9)×K(t-1)
考虑地层温度的变化,由能量守恒定律可知:
ΔQ=Q1+Q2+Q3-Q4
根据能量守恒定律,得到多元热流体的能量守恒方程:
约束方程主要分为4类,分别是饱和度约束方程,摩尔分数约束方程,相平衡方程,毛管力方程;其中,饱和度约束方程是:
Sw+So+Sg=1
摩尔分数约束方程:
x1+x2+x3=1
y1+y2+y3=1
相平衡方程:
K1(p,T)=y2/x1
K2(p,T)=y3/x2
y1=Psat/Pg
毛管力方程:
Pw=Po-Pcow(Sw)
Pg=Po-Pcog(Sg)
边界条件方程:
当给定产量和注入量生产时,内边界条件表述为:
qI=const
当给定井底流动压力生产时,内边界条件表述为:
pwf=const
初始条件:
P=Pi(t=0)
T=Ti(t=0)
多元热流体热采数值模拟理论的差分离散:
为了离散控制方程,将数值模拟模型划分为Nx×Ny×Nz个均匀网格,设网格(i,j,k)的空间步长为Δx,Δy,Δz,时间步长为Δt,上标n表示t=n·Δt时刻,
流体的导热系数,取值为两个相邻网格的调和平均,即:
网格中流体的速度采用下式表示:
上式中,绝对渗透率与粘度之比等于两个相邻网格的调和平均,即:
相对渗透率取流体源流动方向上网格的值,即:
综上,对控制方程进行差分离散,水组分的质量守恒方程的离散形式为:
溶解气组分的质量守恒方程的离散形式为:
轻质油组分的质量守恒方程的差分离散形式为:
重质油组分的质量守恒方程的离散形式为:
能量守恒方程的离散形式为:
具体求解过程:
采用Newton-Raphson迭代法对非线性方程进行线性化求解,4个组分的离散方程写为:
F(X)=0
式中,向量F代表4个组分的离散方程,向量X代表待求的变量;利用Newton-Raphson迭代方法,对第m个方程Fm进行泰勒展开,得到:
上式表示为:
Fm(X+δX)=Fm(X)+J·δX+o(δX2)
其中J是雅可比行列式,采用下式计算
令Fm(X+δX)=0,并且忽略o(δX2),则得到:
J·δX=-F
采用共轭梯度法进行求解,在求解得到了δX后,判断是否收敛,若不收敛,则把δX作为上一个迭代初值的增量,得到下一个新迭代步的数值,
Xnew=Xold+δX
不断的重复迭代,直到满足下列收敛条件:
通过以上步骤,求得离散方程组的解。
附图说明
图1是PVT实验采用的装置示意图;
图2是储层伤害评价实验采用的装置示意图;
图3是岩心的原始渗透率和经过多元热流体伤害后的岩心渗透率测试结果图;
图4是渗透率与多元热流体注入PV数关系曲线;
图5是用于优化BIP的流程图。
具体实施方式
下面通过具体实施方式进一步详细说明:
本实施例公开了一种多元热流体热采油藏数值模拟方法,包括以下步骤:
步骤一:进行PVT实验,PVT实验采用的装置如图1所示;利用测试得到的数据得到新型多元热流体PR-EOS状态方程;基于新型多元热流体PR-EOS状态方程求出多元热流体-稠油体系的相图。
其中PVT实验的具体操作步骤为:
(1)准备PVT筒,将PVT筒清洗干净并抽真空,将一定量的多元热流体和稠油样品以80℃恒定温度转入PVT腔中;
(2)保持混合物在压力大于200bar的条件下注入PVT筒,且初始状态时为液态单相;
(3)将混合物在80℃的条件下搅拌均匀并保持12h以上;
(4)降低压力,记录在200bar,150bar,100bar,50bar四个压力点的体积,保证每一个体积记录点相态平衡;
(5)绘制PV曲线,拐点为饱和压力;
(6)在150℃和280℃的条件下,重复实验。
其中PR-EOS模型为:
其中,
其中P为压力,T为温度,V为摩尔体积,Pc为临界压力,Tc为临界温度,R为通用气体常数,Tr为对比温度;
Soave函数为,
α(Tr)=[1+m(1-Tr 0.5)]2
如果偏心因子小于0.49,
m=0.37464+1.54226ω-0.26992ω2
如果偏心因子大于0.49,
m=0.3796+1.485ω-0.1644ω2+0.01667ω3
对于混合物,应用范德瓦尔斯混合规则计算参数a和b如下:
其中,kij是两个分量之间的BIP(二元相互作用参数)。
BIP影响PR-EOS模型的准确性,可采用临界体积法或临界温度法计算BIP,
采用临界体积法计算BIP:
其中kij是第i个分量和第j个分量之间的BIP(二元相互作用参数),vci和vcj分别为第i分量和第j分量的临界体积,θ为该临界体积法的指数常数。
采用临界温度法计算BIP:
其中kij是第i个分量和第j个分量之间的BIP(二元相互作用参数),Tci和Tcj分别为第i分量和第j分量的临界温度,σ为指数常数。
在本实施例中采用临界体积法进行计算。BIP与温度有关,对BIP在不同温度下的指数进行优化,图5显示了用于优化BIP的流程图,其中图5中的绝对平均相对偏差(AARD)计算如下:
步骤二:进行储层伤害评价实验,描述多元热流体对储层渗透率的伤害程度,建立多元热流体储层伤害数学模型。
其中土层伤害评价实验的具体操作步骤为:
(1)选用高温高压岩心驱替装置进行实验,高温高压岩心驱替装置包括蒸汽发生器和恒温箱,该装置的具体结构如图1所示;实验用岩心是新疆某油田实际岩心,岩心经过洗油烘干,测量岩心的基础参数,岩心基础参数表如表1所示:
表1
(2)利用岩心饱和地层水作为驱替介质,测试岩心的水测渗透率;
(3)设置多元热流体注入速度是1mL/min;调整蒸汽发生器的温度与方案设计温度相同,调整恒温箱的温度与地层温度相同,当蒸汽发生器温度和恒温箱温度达到设计要求后,注入24PV多元热流体,记录注入多元热流体过程中的压力数值;
(4)当注入方案设计的多元热流体的PV数后,注入岩心饱和地层水,第二次测试岩心饱和地层水的测渗透率;
(5)变化多元热流体的温度,重复上述实验步骤,多元热流体的温度分别是150℃,200℃,250℃和300℃。
利用高温高压岩心驱替装置,按照实验步骤分别测试了不同温度条件下储层伤害后的储层渗透率。岩心的原始渗透率和经过多元热流体伤害后的岩心渗透率测试结果见图3。从图3的实验结果可以看到,多元热流体对渗透率500×10-3μm2的岩心和渗透率1200×10-3μm2的岩心均会造成伤害,多元热流体温度越高,对储层渗透率的伤害程度越严重。
将储层伤害后的渗透率数值除以原始储层渗透率数值定义为储层渗透率伤害系数(k/ki),统计储层渗透率伤害系数结果与多元热流体注入PV数结果绘制于图4中(散点是实验测试数据,曲线是方程拟合数据)。通过图4可以看出,多元热流体注入岩心2PV后,岩心的渗透迅速下降,多元热流体温度是573K时,渗透率下降为原始渗透率的62%,多元热流体温度时413K时,渗透率下降为原始渗透率的88%;多元热流体注入岩心6PV后,岩心渗透率下降速度明显变缓,表明多元热流体并不会对储层造成持续性的伤害。
通过实验结果可以看到,多元热流体储层渗透率下降系数会受到多元热流体温度和多元热流体注入PV数的影响。为了进一步描述多元热流体对储层渗透率的伤害程度,本实施例建立了多元热流体储层伤害数学模型,见公式A,利用多元非线性回归方法,计算得到了多元热流体储层伤害数学模型的待定系数,见表2。将公式A的计算结果与实验实际测试结果进行对比(图3),可以发现本实施例建立的多元热流体储层伤害数学模型准确可靠。
参数 | 数值 | 参数 | 数值 | 参数 | 数值 |
beta1 | -7.199799E-07 | beta2 | 5.652760E-05 | beta3 | -1.709992E-03 |
beta4 | 2.478296E-02 | beta5 | -1.711991E-01 | beta6 | -7.771773E-08 |
beta7 | 5.550085E-05 | beta8 | -1.430857E-02 | beta9 | 2.29120371 |
表2
步骤三:利用步骤一、步骤二得到的新型多元热流体PR-EOS状态方程和储层伤害数学模型,建立多元热流体数值模拟方法,并进行编程实现,得到了多元热流体数值模拟器;通过数值模拟器得到多个在不同条件下多元热流体热采的施工方案,比较不同的施工方案,得到最优的多元热流体热采施工方案。
得到数值模拟器的过程为:
建立多元热流体数值模拟方法,必须要考虑多元热流体的温度变化,相态变化,多元热流体与储层岩石之间的相互作用,多元热流体与储层流体的相互作用。多元热流体渗流规律极其复杂,本实施例在热力采油数值模拟方法的基础上,综合考虑多元热流体的特点,对模型进行了如下的基本假设:
(1)流体中存在油、水、气三相流体;
(2)流体流动符合达西定律;
(3)流体流动过程中没有化学反应;
(4)流体由3相(油相、水相、气相)4组分(轻质油组分、重质油组分、水组分、气组分)组成,组分的相态变化服从相平衡原理;
(5)流体渗流是不等温渗流过程;
(6)热量传递的方式是对流和传热;
数值模拟模型由4大类方程构成,分别是质量守恒方程、能量守恒方程、约束方程和边界条件方程。其中,质量守恒方程的基础是质量守恒原理。对于任意组分,质量守恒原理可以表述为
ΔM=M1-M2+M3 (1)
根据质量守恒原理,水组分的质量守恒方程是:
气组分的质量守恒方程是:
轻质油组分的质量守恒方程是:
重质油组分的质量守恒方程是:
考虑到多元热流体对储层渗透率的伤害,多元热流体数值模拟模型中,绝对渗透率不再是始终不变的,本实施例采用显式的时变渗透率处理方法,根据第t-1时间步的绝对渗透率,PV数和温度,计算第t时间步绝对渗透率,具体计算方法如下:
K(t)=(beta1×PV5+beta2×PV4+beta3×PV3+beta4×PV2+beta5×PV+beta6×T3+beta7×T2+beta8×T+beta9)×K(t-1) (6)
多元热流体数值模拟过程中,必须要考虑地层温度的变化。由能量守恒定律可知:
ΔQ=Q1+Q2+Q3-Q4 (7)
根据能量守恒定律,可以得到多元热流体的能量守恒方程:
约束方程主要分为4类,分别是饱和度约束方程,摩尔分数约束方程,相平衡方程,毛管力方程。其中,饱和度约束方程是:
Sw+So+Sg=1 (9)
摩尔分数约束方程:
x1+x2+x3=1 (10)
y1+y2+y3=1 (11)
相平衡方程:
K1(p,T)=y2/x1 (12)
K2(p,T)=y3/x2 (13)
y1=Psat/Pg (14)
毛管力方程:
Pw=Po-Pcow(Sw) (15)
Pg=Po-Pcog(Sg) (16)
边界条件方程:
当给定产量和注入量生产时,内边界条件表述为:
qI=const (17)
当给定井底流动压力生产时,内边界条件表述为:
pwf=const (18)
初始条件:
P=Pi(t=0) (19)
T=Ti(t=0) (20)
多元热流体热采数值模拟理论的差分离散:
为了离散控制方程,将数值模拟模型划分为Nx×Ny×Nz个均匀网格,设网格(i,j,k)的空间步长为Δx,Δy,Δz,时间步长为Δt。上标n表示t=n·Δt时刻。
流体的导热系数,取值为两个相邻网格的调和平均,即:
网格中流体的速度采用下式表示
上式中,绝对渗透率与粘度之比等于两个相邻网格的调和平均,即:
相对渗透率取流体源流动方向上网格的值,即:
综上所述,可以对控制方程进行差分离散。水组分的质量守恒方程的离散形式为:
溶解气组分的质量守恒方程的离散形式为:
轻质油组分的质量守恒方程的差分离散形式为:
重质油组分的质量守恒方程的离散形式为:
能量守恒方程的离散形式为:
具体求解过程:
每个网格节点包含4个离散控制方程,每个离散控制方程具有强非线性,为了提高计算稳定性和加快收敛速度,本实施例采用Newton-Raphson迭代法对非线性方程进行线性化求解。4个组分的离散方程可以写为:
F(X)=0
式中,向量F代表4个组分的离散方程,向量X代表待求的变量。利用Newton-Raphson迭代方法,对第m个方程Fm进行泰勒展开,可以得到:
上式可以表示为:
Fm(X+δX)=Fm(X)+J·δX+o(δX2)
其中J是雅可比行列式,可以采用下式计算:
令Fm(X+δX)=0,并且忽略o(δX2),则可以得到:
J·δX=-F
行列式J是一个大型的稀疏矩阵,本实施例采用共轭梯度法进行求解。在求解得到了δX后,判断是否收敛,若不收敛,则把δX作为上一个迭代初值的增量,得到下一个新迭代步的数值。
Xnew=Xold+δX
不断的重复迭代,直到满足下列收敛条件:
通过以上步骤,可以得到离散方程组。
Claims (7)
1.一种多元热流体热采油藏数值模拟方法,其特征在于:包括以下步骤:
步骤一:进行PVT实验,利用测试得到的数据得到新型多元热流体PR-EOS状态方程;基于新型多元热流体PR-EOS状态方程求出多元热流体-稠油体系的相图;
步骤二:进行储层伤害评价实验,描述多元热流体对储层渗透率的伤害程度,建立多元热流体储层伤害数学模型;
步骤三:利用步骤一、步骤二得到的新型多元热流体PR-EOS状态方程和储层伤害数学模型,建立多元热流体数值模拟方法,并进行编程实现,得到多元热流体数值模拟器;通过数值模拟器得到多个在不同条件下多元热流体热采的施工方案,比较不同的施工方案,得到最优的多元热流体热采施工方案。
5.根据权利要求1所述的一种多元热流体热采油藏数值模拟方法,其特征在于:步骤二中,储层伤害评价实验的具体操作步骤为:
(1)选用高温高压岩心驱替装置进行实验,高温高压岩心驱替装置包括蒸汽发生器和恒温箱;
(2)利用地层水作为驱替介质,测试岩心的水测渗透率;
(3)设置多元热流体注入速度是1mL/min;调整蒸汽发生器的温度与方案设计温度相同,调整恒温箱的温度与地层温度相同,当蒸汽发生器温度和恒温箱温度达到设计要求后,注入24PV多元热流体,记录注入多元热流体过程中的压力数值;
(4)当注入方案设计的多元热流体的PV数后,注入地层水,第二次测试岩心的水测渗透率;
(5)变化多元热流体的温度,重复上述实验步骤,多元热流体的温度分别是150℃,200℃,250℃和300℃。
7.根据权利要求1所述的一种多元热流体热采油藏数值模拟方法,其特征在于:步骤三中,得到数值模拟器的过程为:
数值模拟模型由4大类方程构成,分别是质量守恒方程、能量守恒方程、约束方程和边界条件方程;其中,质量守恒方程的基础是质量守恒原理,对于任意组分,质量守恒原理表述为:
ΔM=M1-M2+M3
根据质量守恒原理,水组分的质量守恒方程是:
气组分的质量守恒方程是:
轻质油组分的质量守恒方程是:
重质油组分的质量守恒方程是:
采用显式的时变渗透率处理方法,根据第t-1时间步的绝对渗透率,PV数和温度,计算第t时间步绝对渗透率,具体计算方法如下:
K(t)=(beta1×PV5+beta2×PV4+beta3×PV3+beta4×PV2+beta5×PV+beta6×T3+beta7×T2+beta8×T+beta9)×K(t-1)
考虑地层温度的变化,由能量守恒定律可知:
ΔQ=Q1+Q2+Q3-Q4
根据能量守恒定律,得到多元热流体的能量守恒方程:
约束方程主要分为4类,分别是饱和度约束方程,摩尔分数约束方程,相平衡方程,毛管力方程;其中,饱和度约束方程是:
Sw+So+Sg=1
摩尔分数约束方程:
x1+x2+x3=1
y1+y2+y3=1
相平衡方程:
K1(p,T)=y2/x1
K2(p,T)=y3/x2
y1=Psat/Pg
毛管力方程:
Pw=Po-Pcow(Sw)
Pg=Po-Pcog(Sg)
边界条件方程:
当给定产量和注入量生产时,内边界条件表述为:
qI=const
当给定井底流动压力生产时,内边界条件表述为:
pwf=const
初始条件:
P=Pi(t=0)
T=Ti(t=0)
多元热流体热采数值模拟理论的差分离散:
为了离散控制方程,将数值模拟模型划分为Nx×Ny×Nz个均匀网格,设网格(i,j,k)的空间步长为Δx,Δy,Δz,时间步长为Δt,上标n表示t=n·Δt时刻,
流体的导热系数,取值为两个相邻网格的调和平均,即:
网格中流体的速度采用下式表示:
上式中,绝对渗透率与粘度之比等于两个相邻网格的调和平均,即:
相对渗透率取流体源流动方向上网格的值,即:
综上,对控制方程进行差分离散,水组分的质量守恒方程的离散形式为:
溶解气组分的质量守恒方程的离散形式为:
轻质油组分的质量守恒方程的差分离散形式为:
重质油组分的质量守恒方程的离散形式为:
能量守恒方程的离散形式为:
具体求解过程:
采用Newton-Raphson迭代法对非线性方程进行线性化求解,4个组分的离散方程写为:
F(X)=0
式中,向量F代表4个组分的离散方程,向量X代表待求的变量;利用Newton-Raphson迭代方法,对第m个方程Fm进行泰勒展开,得到:
上式表示为:
Fm(X+δX)=Fm(X)+J·δX+o(δX2)
其中J是雅可比行列式,采用下式计算
令Fm(X+δX)=0,并且忽略o(δX2),则得到:
J·δX=-F
采用共轭梯度法进行求解,在求解得到了δX后,判断是否收敛,若不收敛,则把δX作为上一个迭代初值的增量,得到下一个新迭代步的数值,
Xnew=Xold+δX
不断的重复迭代,直到满足下列收敛条件:
通过以上步骤,求得离散方程组的解。
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010176052.XA CN111191389B (zh) | 2020-03-13 | 2020-03-13 | 一种多元热流体热采油藏数值模拟方法 |
PCT/CN2021/080333 WO2021180189A1 (zh) | 2020-03-13 | 2021-03-12 | 一种多元热流体热采油藏数值模拟方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010176052.XA CN111191389B (zh) | 2020-03-13 | 2020-03-13 | 一种多元热流体热采油藏数值模拟方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111191389A CN111191389A (zh) | 2020-05-22 |
CN111191389B true CN111191389B (zh) | 2022-03-01 |
Family
ID=70708664
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010176052.XA Active CN111191389B (zh) | 2020-03-13 | 2020-03-13 | 一种多元热流体热采油藏数值模拟方法 |
Country Status (2)
Country | Link |
---|---|
CN (1) | CN111191389B (zh) |
WO (1) | WO2021180189A1 (zh) |
Families Citing this family (17)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111191389B (zh) * | 2020-03-13 | 2022-03-01 | 重庆科技学院 | 一种多元热流体热采油藏数值模拟方法 |
CN113327647B (zh) * | 2021-04-13 | 2022-12-30 | 中国石油大学(华东) | 一种模拟水合物藏中不同来源甲烷运移轨迹的方法 |
CN114136838B (zh) * | 2021-11-19 | 2023-11-17 | 中国海洋石油集团有限公司 | 确定海上稠油不同含水阶段注热水驱原油粘度界限的方法 |
CN114352248B (zh) * | 2021-12-17 | 2023-09-26 | 常州大学 | 一种稠油热采二维物理模拟实验装置及其使用方法 |
CN114295791A (zh) * | 2022-03-09 | 2022-04-08 | 新疆广陆能源科技股份有限公司 | 热流体发生器尾气检测系统和多元介质热流体发生系统 |
CN114357810B (zh) * | 2022-03-18 | 2022-05-17 | 中国石油大学(华东) | 一种页岩气藏有机质孔隙内流体临界参数及相图计算方法 |
CN114861396A (zh) * | 2022-03-30 | 2022-08-05 | 西北核技术研究所 | 考虑细砂冲击压缩特性和吸热相变的数学模型及建模方法 |
CN114893172B (zh) * | 2022-05-26 | 2023-08-08 | 常州大学 | 一种模拟稠油热流体驱替的方法及系统 |
CN115030699A (zh) * | 2022-06-30 | 2022-09-09 | 常州大学 | 一种用于油砂sagd二维可视化模拟实验追踪加热装置 |
CN116130015B (zh) * | 2022-12-22 | 2024-04-30 | 大连理工大学 | 孔隙尺度下二氧化碳驱油封存数值模拟方法 |
CN116050293B (zh) * | 2022-12-30 | 2023-08-08 | 常州大学 | 一种确定稠油泡沫调驱技术参数最优数值指示图版的方法 |
CN116044355B (zh) * | 2023-03-30 | 2023-07-04 | 中国石油大学(华东) | 一种模拟非均一粘度场的可视化驱替装置及其工作方法 |
CN116307280B (zh) * | 2023-05-18 | 2023-08-01 | 成都理工大学 | 一种酸性气藏气井储层硫堵损伤量化评估方法 |
CN116738779B (zh) * | 2023-05-24 | 2024-01-26 | 成都理工大学 | 超临界地热流体电导率的计算方法 |
CN117113873B (zh) * | 2023-08-15 | 2024-04-09 | 西南石油大学 | 一种多相流体地层渗流的数值模拟方法及应用 |
CN117217393B (zh) * | 2023-11-08 | 2024-01-26 | 新疆智能港环保科技有限公司 | 一种通过渗析扩容提高油气井产量检测修正系统 |
CN117688874B (zh) * | 2024-02-04 | 2024-04-16 | 北京理工大学 | 一种超临界流体热力学状态判别与计算方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN201705320U (zh) * | 2010-05-24 | 2011-01-12 | 中国石油天然气股份有限公司 | 一种注多元热流体采油三维模拟试验装置 |
CN106886649A (zh) * | 2017-03-01 | 2017-06-23 | 中国海洋石油总公司 | 一种多元热流体吞吐注入参数优化方法 |
CN110097254A (zh) * | 2019-04-01 | 2019-08-06 | 中国石油化工股份有限公司 | 多元热流体驱油藏潜力的动态评价方法 |
Family Cites Families (16)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5980695A (en) * | 1996-11-14 | 1999-11-09 | Transglobal Industrial Services, Inc. | Apparatus for vacuum distillation of solvents |
KR100326599B1 (ko) * | 1998-06-27 | 2002-11-11 | 문상우 | 소각기의 폐열이용장치 |
US9020793B2 (en) * | 2005-12-22 | 2015-04-28 | Chevron U.S.A. Inc. | Method, system and program storage device for reservoir simulation utilizing heavy oil solution gas drive |
CN101476458B (zh) * | 2008-12-03 | 2010-12-08 | 刘其成 | 一种油藏开发模拟系统、油藏模型本体及其数据处理方法 |
CN101818636B (zh) * | 2010-05-24 | 2013-04-24 | 中国石油天然气股份有限公司 | 注多元热流体采油三维模拟试验装置 |
CN102230372A (zh) * | 2011-06-24 | 2011-11-02 | 中国海洋石油总公司 | 一种稠油井多元热流体热采工艺 |
CN102606121B (zh) * | 2012-03-15 | 2015-07-22 | 中国海洋石油总公司 | 一种用于稠油油藏的多元热流体生产工艺以及热采工艺 |
CN102678096B (zh) * | 2012-05-08 | 2015-07-08 | 中国石油天然气股份有限公司 | 一种热水辅助重力泄油开采高凝油油藏的方法 |
CN104265255B (zh) * | 2014-09-26 | 2017-01-25 | 中国石油天然气股份有限公司 | 一种稠油二维微观可视化驱替模拟实验系统及其使用方法 |
KR101715513B1 (ko) * | 2015-10-06 | 2017-03-13 | 충북대학교 산학협력단 | 열매유를 이용한 케로신의 열전달, 열분해 실험장치 |
CN105422068B (zh) * | 2015-11-12 | 2017-10-17 | 中国石油天然气股份有限公司 | 水平井分段体积压裂和压裂充填组合开发稠油油藏的方法 |
CN105909214A (zh) * | 2016-04-14 | 2016-08-31 | 中国石油大学(华东) | 一种利用长水平井自循环开采致密干热岩地热能的方法 |
CN106294282B (zh) * | 2016-08-01 | 2019-04-12 | 中国石油天然气股份有限公司 | 黑油油藏模拟方法及装置 |
CN106761733A (zh) * | 2017-01-06 | 2017-05-31 | 中国海洋石油总公司 | 一种稠油油藏水平井蒸汽吞吐初期产能预测方法 |
CN110685653A (zh) * | 2019-10-11 | 2020-01-14 | 中海石油(中国)有限公司 | 考虑变启动压力梯度的水驱稠油数值模拟方法 |
CN111191389B (zh) * | 2020-03-13 | 2022-03-01 | 重庆科技学院 | 一种多元热流体热采油藏数值模拟方法 |
-
2020
- 2020-03-13 CN CN202010176052.XA patent/CN111191389B/zh active Active
-
2021
- 2021-03-12 WO PCT/CN2021/080333 patent/WO2021180189A1/zh active Application Filing
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN201705320U (zh) * | 2010-05-24 | 2011-01-12 | 中国石油天然气股份有限公司 | 一种注多元热流体采油三维模拟试验装置 |
CN106886649A (zh) * | 2017-03-01 | 2017-06-23 | 中国海洋石油总公司 | 一种多元热流体吞吐注入参数优化方法 |
CN110097254A (zh) * | 2019-04-01 | 2019-08-06 | 中国石油化工股份有限公司 | 多元热流体驱油藏潜力的动态评价方法 |
Also Published As
Publication number | Publication date |
---|---|
CN111191389A (zh) | 2020-05-22 |
WO2021180189A1 (zh) | 2021-09-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111191389B (zh) | 一种多元热流体热采油藏数值模拟方法 | |
Settari | Simulation of hydraulic fracturing processes | |
CN107818188B (zh) | 一种针对流体加热油页岩原位开采的模拟预测方法 | |
Li et al. | Modeling tracer flowback in tight oil reservoirs with complex fracture networks | |
Pahari et al. | Optimal pumping schedule with high-viscosity gel for uniform distribution of proppant in unconventional reservoirs | |
Cheng et al. | Mechanics of hydraulic fracturing | |
CN105089649A (zh) | 一种井间连通性模型建立方法 | |
Yu et al. | Modified zipper fracturing in enhanced geothermal system reservoir and heat extraction optimization via orthogonal design | |
CN111075441A (zh) | 一种边底水稠油油藏热采后冷采三维物理模拟实验装置及方法 | |
Xiong et al. | A novel model to investigate the effects of injector-producer pressure difference on SAGD for bitumen recovery | |
Ito | The introduction of the microchanneling phenomenon to cyclic steam stimulation and its application to the numerical simulator (sand deformation concept) | |
Guo et al. | A new SAGD comprehensive multi-stage model for oil production using a concave parabola geometry | |
Ma et al. | Embedded discrete fracture modeling: Flow diagnostics, non-Darcy flow, and well placement optimization | |
Hemmati et al. | Laboratory validation of steady-state-transient test to determine relative permeability and capillary pressure | |
Shi et al. | Injection energy utilization efficiency and production performance of oil shale in-situ exploitation | |
Sharafisafa et al. | Coupled thermo-hydro-mechanical simulation of hydraulic fracturing in deep reservoirs using finite-discrete element method | |
Liu et al. | A heat and mass transfer coupling model for transition interface of expanding solvent steam-assisted gravity drainage | |
Gonzalez Abad | Development of a compositional reservoir simulator for asphaltene precipitation based on a thermodynamically consistent model | |
Takbiri-Borujeni et al. | Upscaling the steam-assisted-gravity-drainage model for heterogeneous reservoirs | |
Du et al. | Expansion velocity model of steam-assisted gravity drainage considering thermal convection | |
Chen et al. | The application of Stefan problem in calculating the lateral movement of steam chamber in SAGD | |
Zhang et al. | Prediction of the thermal conductivity of rocks using group method of data handling (GMDH) | |
Yu et al. | A multicomponent thermal fluid numerical simulation method considering formation damage | |
Zhang et al. | Hydrocarbon gas huff-n-puff optimization of multiple horizontal wells with complex fracture networks in the M unconventional reservoir | |
Wang et al. | Experimental Study on the Variations of Oil Sand Permeability during SAGD Processes |
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 |