CN107621350A - 一种模拟超临界co2驱替天然气流动的方法 - Google Patents
一种模拟超临界co2驱替天然气流动的方法 Download PDFInfo
- Publication number
- CN107621350A CN107621350A CN201610562387.9A CN201610562387A CN107621350A CN 107621350 A CN107621350 A CN 107621350A CN 201610562387 A CN201610562387 A CN 201610562387A CN 107621350 A CN107621350 A CN 107621350A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- msubsup
- mfrac
- mtd
- 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
Landscapes
- Organic Low-Molecular-Weight Compounds And Preparation Thereof (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明提出了一种模拟超临界CO2驱替天然气流动的方法,所述方法包括如下步骤:根据气藏地质特征,建立三维地质模型;通过改进状态方程建立流体模型确定气藏流体的基础参数;确定CO2的扩散系数;建立考虑扩散的超临界CO2驱天然气多组分数值模拟模型;应用模型计算不同CO2注入量和天然气采出量的对应关系。本发明提出了新的数学模型及处理方法模拟超临界CO2驱替天然气流动特征,与常规数值模型相比能更好的评价CO2驱气流动特点及驱替效果。
Description
技术领域
本发明属于气藏开发领域,具体地,涉及一种模拟超临界CO2驱替天然气流动的方法。
背景技术
由于人类对石化燃料(煤、石油、天然气)的过度依赖,工业和人类生活过程中产生的温室气体排放量日益增加,由此导致的空气污染和温室效应正在严重地威胁着人类赖以生存的环境。温室气体二氧化碳(CO2)使全球气候变暖,对人类的生存和社会经济的发展构成了严重的威胁。CO2有效利用已成为了全球广泛关注的课题。通过不断的发展及完善,在油藏中注入CO2成为较为成熟的技术。
但在气藏中注入CO2相态和驱替机理与油藏完全不同,注CO2提高气藏采收率成为一项全新的课题。国外在上世纪九十年代中后期开始CO2EGR研究,国内近两年初步开展室内理论研究。但目前国内外仅开展实验及机理模拟,仍处于探索阶段,未进入现场应用,此项研究仍处于起步阶段。
在CO2与天然气驱替模拟方面,国内外均应用数值模拟软件开展了CO2提高采收率机理研究。Oldenburg等人应用TOUGH2软件开展等温多相多组分数值模拟,模拟结果表明底部注气有利于减少CO2向上锥进和延缓CO2与CH4混合,提高平面及垂向驱替效率。TorstenClemens以实例气田为例模拟CO2提高采收率驱替效率。模拟结果表明,气藏开采到废弃阶段注入CO2能提高采收率9.4%。A.AL-Hashami应用机理模型模拟底水气藏CO2驱替效率,通过数模研究表明CO2有助于提高气藏采收率,能增产8-11%。目前,普遍采用商业数值模拟软件分析CO2提高采收率驱替封存效率,模型考虑因素较为简单。未考虑CO2扩散对驱替效果的影响。CO2驱气机理不明确,驱替机理论证不系统,不充分。
发明内容
为了解决现有技术的不足,本发明探索提出体积改进状态方程方法,解决低温时高含CO2气体实验数据拟合难的问题;并且考虑了不同压力条件下CO2扩散对驱替天然气的影响,建立了考虑扩散的超临界CO2驱天然气多组分数值模拟模型,模拟超临界CO2驱替天然气驱替机理及驱替效果。
根据本发明的一个方面,提供一种模拟超临界CO2驱替天然气流动的方法,所述方法包括如下步骤:
根据气藏地质特征,建立三维地质模型;
通过改进状态方程建立流体模型确定气藏流体的基础参数;
确定CO2的扩散系数;
建立考虑扩散的超临界CO2驱天然气多组分数值模拟模型;
应用模型计算不同CO2注入量和天然气采出量的对应关系。
进一步地,确定气藏流体的基础参数,包括:实验测试不同CO2天然气混合物偏差因子、密度及粘度。
进一步地,改进状态方程,获得描述整个压力温度范围内临界性质的热力学参数。
进一步地,引入体积修正项对状态方程进行修正,对临界点体积参数进行改进,改进项由两部分组成:①临界点改进;②温度改进:
c(T)=ccf(Tr) (1)
其中,c(T)体积修正项,m3/mol;cc临界点改进体积修正项,m3/mol;f(Tr)温度改进系数,小数。
进一步地,引入临界点改进项:
其中,临界压缩因子Zc为定值,R气体常数,8.314J/mol·K;Tc临界温度,K;Pc临界压力,MPa;vEOS状态方程计算摩尔体积,m3/mol;vExp实验计算摩尔体积,m3/mol;cc体积修正项,m3/mol。
进一步地,引入温度改进项进行修正,温度改进通过相平衡实验确定Zc与T的对应关系,通过实验拟合得方程式(3)-(5):
f(Tr)=β+(1-β)exp(0.5β+8.3) (3)
β=-2.8431exp([-64.2184(0.3074-Zc)]+0.1735) (4)
β=T-Tr (5)
其中,T气体所处的温度,K;Tr=T/Tc对比温度,K;β温度校正项,K。
进一步地,将临界点改进项c(T)=ccf(Tr)带入PR状态方程:
状态方程中,ab由临界参数和偏心因子计算。式中,
a=α(T)·ac;
m=0.37464+1.54226ω-0.26992ω2;
Tr=T/Tc;
b=0.77796·R·Tc/Pc
其中,p压力,MPa;v摩尔体积,m3/mol;ω偏心因子,无因次。
进一步地,确定CO2的扩散系数包括:实验测量不同压力条件下CO2在天然气中的扩散浓度变化规律,应用插值计算方法描述多孔介质中CO2的扩散系数。
进一步地,建立考虑扩散的超临界CO2驱天然气多组分数值模拟模型包括:应用改进状态方程劈分计算各组分热力学参数,质量守恒方程中考虑扩散项对流动的影响。
进一步地,考虑扩散的超临界CO2驱天然气多组分数值模拟模型包括,如下:
其中,k渗透率,D;krg气相渗透率,小数;μg气相粘度mPa·s;ρg气相密度kg/m3;yig组分i在气相中的摩尔分数,小数;气相的势;φ孔隙度,小数;Sg气相饱和度,小数;qik组分i在k层中的产量或流量,m2/d;Dig为组分i在气相中的扩散系数,m2/d;nf为井的射孔层段;
扩散张量:
式中,
其中,ρi组分i密度kg/m3;Sj j相饱和度,小数;Dig为组分i在气相中的扩散系数,m2/d;为组分i在j相中的分子扩散系数,m2/d;uj为速度相量,m2/d;ujx为x方向速度相量,m2/d;ujy为y方向速度相量,m2/d;ujz为z方向速度相量,m2/d;τj迂曲率,小数;j相的沿程弥散系数,小数;αjt j相的横向弥散系数,小数。
进一步地,应用模型计算不同CO2注入量和天然气采出量的对应关系包括,根据长岩心驱替实验的物性参数,建立实验尺度数值模型,导入束缚水饱和度、渗透率、孔隙度参数,计算不同CO2注入量和天然气采出量的对应关系。
本发明解决了超临界CO2天然气混合物近临界区高压物性拟合难的问题。
常规气体偏差系数计算模型主要有:Dranchuk–Purvis-Robinson(DPR)法、Hall&Yarborough(HY)法、Sarem方法、Dranchuk-Abu-Kassem(DAK)法、李相方(LXF)法等。采用经验公式对比实验数据,发现计算误差较大;而应用状态方程,高温区域热力学参数计算误差小,但在低温(近临界温度)区且处于近临界区(近临界压力)时,热力学参数误差较大。
应用改进的状态方程,计算偏差因子与实验值误差小于3%、计算不同CO2含量天然气混合物高压物性参数(密度、粘度与实验值)误差小于3%。
本发明考虑不同压力条件下CO2扩散对驱替天然气的影响,描述超临界CO2天然气驱替前沿及驱替效率。
通过长岩心驱替实验,获得CO2注入量与出口端CO2摩尔含量及天然气采收率对应关系数据。根据实验数据,建立长岩心驱替模型,在相态拟合的基础上,分别应用常规组分模型和本文考虑吸附的组分模型计算在长岩心中CO2驱替天然气,使用本发明方法建立的组分模型计算的采收率与实验值更接近。
附图说明
通过结合附图对本公开示例性实施方式进行更详细的描述,本公开的上述以及其它目的、特征和优势将变得更加明显,其中,在本公开示例性实施方式中,相同的参考标号通常代表相同部件。
图1为根据本发明实施例的方法流程图。
图2显示了Zc与β实验测试拟合图。
图3显示了CO2在甲烷中的扩散系数。
图4显示了状态方程热力学参数计算结果。
图5显示了不同CO2含量天然气混合物密度拟合图。
图6显示了不同CO2含量天然气混合物粘度拟合图。
图7显示了注入CO2孔隙体积与出口端CO2摩尔含量计算数据。
图8显示了注入CO2孔隙体积与采收率计算数据。
具体实施方式
下面将参照附图更详细地描述本公开的优选实施方式。虽然附图中显示了本公开的优选实施方式,然而应该理解,可以以各种形式实现本公开而不应被这里阐述的实施方式所限制。相反,提供这些实施方式是为了使本公开更加透彻和完整,并且能够将本公开的范围完整地传达给本领域的技术人员。
本发明提供了一种模拟超临界CO2驱替天然气流动的方法,属于注气提高采收率模拟方法领域。本方法以组分模型为基础,基于超临界CO2气体的相行为特殊,考虑近临界及温度的影响,改进状态方程,探索描述超临界CO2注气过程中与天然气混合的组分表征方法。在流动表征方面,在流动项表征时考虑不同压力CO2扩散对驱替天然气的影响,建立了考虑扩散的超临界CO2驱天然气多组分数值模拟模型,实现了超临界CO2驱气前缘的模拟,本发明提出了新的数学模型及处理方法模拟超临界CO2驱替天然气流动特征,与常规数值模型相比能更好的评价CO2驱气流动特点及驱替效果。
目前,超临界CO2驱天然气多组分数值模拟模型存在以下问题:
在组分表征方面,主要为经验公式法及状态方程法。经验公式法参变量少,易于计算,但是计算误差较大;而状态方程以热力学理论为基础,描述分子间引力斥力作用,模型通用性很好,计算热力学参数为组分模拟器中各组分劈分的基础,但是,在近临界区高压物性计算精度差。CO2气体的相行为特殊,在近临界区严重偏移理想气体,CO2-天然气混合物也具有相似的行为。因此,本发明探索提出描述超临界CO2注气过程中与天然气混合的组分表征方法。
在流动表征方面,超临界CO2驱替天然气用达西定律表征,但是,实验表明CO2扩散对驱替前沿存在影响,流动项表征时需考虑扩散对CO2驱替的影响。因此,需要在实验的基础上,考虑不同压力条件下CO2扩散对驱替天然气的影响。
针对以上存在的问题,本发明提出了新的数学模型及处理方法模拟超临界CO2驱替天然气流动特征,比常规数值模型相比能更好的评价CO2驱气流动特点及驱替效果。
参照图1,本发明提出了一种模拟超临界CO2驱替天然气流动的方法,所述方法包括如下过程:
根据气藏地质特征,建立三维地质模型;
通过改进状态方程建立流体模型确定气藏流体的基础参数;
确定CO2的扩散系数;
建立考虑扩散的超临界CO2驱天然气多组分数值模拟模型;
应用模型计算不同CO2注入量和天然气采出量的对应关系。
进一步地,确定气藏流体的基础参数,包括:实验测试不同CO2天然气混合物偏差因子、密度及粘度;通过引入体积修正项改进状态方程,获得描述整个压力温度范围内临界性质的热力学参数。
优选地,确定扩散系数包括:实验测量不同压力条件下CO2在天然气中的扩散浓度变化规律,应用插值计算方法描述多孔介质中CO2的扩散系数。
优选地,建立考虑扩散的超临界CO2驱天然气多组分数值模拟模型,包括:应用改进状态方程劈分计算各组分热力学参数,质量守恒方程中考虑扩散项对流动的影响。
优选地,根据长岩心驱替实验的物性参数,建立实验尺度数值模型,导入束缚水饱和度、渗透率、孔隙度等基本参数,计算不同CO2注入量和天然气采出量的对应关系,验证模型的实用性。
优选地,根据气藏地质特征,建立三维地质模型。导入气田生产井产气量、注入CO2量及压力数据,进行流体流动预测分析,计算不同方案生产井CO2含量及天然气采出程度。
优选地,超临界CO2-天然气状态方程改进方法,如下:
CO2相行为特殊,当温度大于31.04C,压力大于7.38MPa,CO2具有超临界特性,超临界CO2密度与液体接近,大于天然气密度。在近临界区CO2物性严重偏移理想气体,CO2-天然气混合物也具有相似的行为。应用常规状态方程,CO2-天然气混合物高压物性在近临界区难以拟合,另外,当温度较低,在近临界温度附近时,计算误差较大。
因此引入体积修正项对状态方程进行修正,反映CO2-天然气混合物相变特征。改进项首先考虑立方体状态方程自身存在的缺陷,立方体状态方程临界压缩因子Zc为定值,最为接近真实气体的PR状态方程Zc为0.307,而CO2Zc 0.274真实气体Zc0.23-0.29,存在差异,因此需要对临界点体积参数进行改进。改进项由两部分组成:①临界点改进;②温度改进:
c(T)=ccf(Tr) (1)
其中,c(T)体积修正项,m3/mol;cc临界点改进体积修正项,m3/mol;f(Tr)温度改进系数,小数。
引入临界点改进项:
其中,临界压缩因子Zc为定值,R气体常数,8.314J/mol·K;Tc临界温度,K;Pc临界压力,MPa;vEOS状态方程计算摩尔体积,m3/mol;vExp实验计算摩尔体积,m3/mol;cc体积修正项,m3/mol。
引入温度改进项进行修正,温度改进通过相平衡实验确定Zc与T的对应关系,通过实验拟合得方程式(3)-(5):
f(Tr)=β+(1-β)exp(0.5β+8.3) (3)
β=-2.8431exp([-64.2184(0.3074-Zc)]+0.1735) (4)
β=T-Tr (5)
其中,T气体所处的温度,K;Tr=T/Tc对比温度,K;β温度校正项,K。
进一步地,将临界点改进项c(T)=ccf(Tr)带入PR状态方程:
状态方程中,ab由临界参数和偏心因子计算。式中,
a=α(T)·ac;
m=0.37464+1.54226ω-0.26992ω2;
Tr=T/Tc;
b=0.77796·R·Tc/Pc
其中,p压力,MPa;v摩尔体积,m3/mol;ω偏心因子,无因次。
应用改进的PR状态方程(6)劈分各组分物性参数,获得超临界CO2驱替天然气过程中整个压力温度范围内的热力学参数。
进一步地,考虑扩散的超临界CO2驱天然气多组分数值模拟模型包括,如下:
其中,k渗透率,D;krg气相渗透率,小数;μg气相粘度mPa·s;ρg气相密度kg/m3;yig组分i在气相中的摩尔分数,小数;气相的势,MPa;φ孔隙度,小数;Sg气相饱和度,小数;qik组分i在k层中的产量或流量,m2/d;Dig为组分i在气相中的扩散系数,m2/d;nf为井的射孔层段。
扩散张量:
式中,
其中,ρi组分i密度kg/m3;Sj j相饱和度,小数;Dig为组分i在气相中的扩散系数,m2/d;为组分i在j相中的分子扩散系数,m2/d;uj为速度相量,m2/d;ujx为x方向速度相量,m2/d;ujy为y方向速度相量,m2/d;ujz为z方向速度相量,m2/d;τj迂曲率,小数;j相的沿程弥散系数,小数;αjt j相的横向弥散系数,小数。
优选地,CO2扩散系数的求取,如下:
应用高压PVT仪器测量容器中不同压力条件下不同部位超临界CO2与天然气混合浓度变化规律,计算CO2分子扩散系数。扩散系数与压力有关,绘制扩散系数与压力的关系图,见图3。
与压力有关,通过插值确定不同压力时CO2在甲烷中的分子扩散系数。代入考虑扩散的超临界CO2驱天然气多组分数值模拟模型计算。
为便于理解本发明实施例的方案及其效果,以下给出一个具体应用示例。本领域技术人员应理解,该示例仅为了便于理解本发明,其任何具体细节并非意在以任何方式限制本发明。
下面参照图4至图8,描述本发明的一个应用例,将本发明应用于CO2驱替长岩心实验模拟。
首先,根据气藏地质特征,建立三维地质模型,开展网格划分并粗化处理。应用CMG软件,以长岩心实验为例,建立实验尺度长岩心地质模型,网格确定为横向上20个网格,纵向上1层,网格数量为20×1×1=20个。横向网格长度为5cm,纵向上为2.5cm。平均孔隙度9.9%,渗透率0.652mD,束缚水饱和度为37.3%。横向上第1个网格处设为注入端,第20个网格处设为产出端。每个网格原始地层压力设为25MPa。
接下来,确定流体参数。实验测试不同CO2天然气混合物偏差因子、密度及粘度;应用改进状态方程,计算不同CO2天然气混合物偏差因子、密度及粘度,验证模型的准确性。应用改进状态方程,劈分各组分,获得描述整个压力温度范围内临界性质的热力学参数。将获得的各组分热力学参数录入到CMG软件中。
接下来,确定CO2扩散系数。实验测量CO2在天然气中的不同压力条件下不同部位CO2浓度变化规律。通过实验测试,在扩散作用下,随着压力的不断降低,CO2与天然气的混合时间越来越短,25MPa时CO2和天然气混合的时间为210小时左右,15MPa时两者混合的时间为120小时左右,而5MPa下为25小时左右。进一步定量评价二者的扩散速度,计算多孔介质条件下CO2在甲烷中的扩散系数。将菲克定律推广到测定气体通过岩石的扩散系数。假定岩心长度为Z(而不是孔径长度),面积为S(而不是有效面积),则求得的扩散系数D应是有效扩散系数。公式为
其中
式中:C0‘、C0‘’——定压取样装置气体反应区上部取样点CO2和下部取样点CO2的初始浓度;
Ct‘、Ct‘’——t时间定压取样装置气体反应区上部取样点CO2和下部取样点CO2的浓度;
V‘、V‘’——CO2与天然气的初始体积;
S——定压取样装置内部截面积;
Z——定压取样装置上部和下部取样点的距离。
通过计算获得定压取样装置内25MPa时CO2在天然气中的扩散系数为2.74×10- 4cm2/s,15MPa时的扩散系数为6.55×10-4cm2/s,5MPa时的扩散系数为1.57×10-3cm2/s。
接下来,建立考虑扩散的超临界CO2驱天然气多组分数值模拟模型。与压力有关,通过插值确定不同压力时CO2在甲烷中的分子扩散系数,应用于考虑扩散的超临界CO2驱天然气多组分数值模拟模型(公式6)。在模型应用中,对不同压力步录入扩散系数。组分模型求解方式采用自适应隐式方法。
接下来,模拟多级降压开采,第20个网格处设为产出端,定压力求产。降压间隔均为3MPa,最终压力降至8MPa。第1个网格设为注入端,注入端压力为11MPa,出口端回压为8MPa,驱替压差3MPa,定压力求产。
模型结果能显示不同网格压力,CO2分布特征,密度粘度的变化规律,及不同CO2注入量天然气的产出量。通过数模计算,当出口端CO2为10%时,比衰竭方案提高采收率9.46%。数值模型论证了超临界CO2提高气藏采收率的可行性。
表1长岩心CO2驱替实验结果与数值模型计算结果对比表
应用上述方程,计算实验数据,与实验值拟合误差小于5%注入倍数在0.5左右时出口端气体中CO2摩尔含量超过10%。
以上已经描述了本公开的各实施例,上述说明是示例性的,并非穷尽性的,并且也不限于所披露的各实施例。在不偏离所说明的各实施例的范围和精神的情况下,对于本技术领域的普通技术人员来说许多修改和变更都是显而易见的。本文中所用术语的选择,旨在最好地解释各实施例的原理、实际应用或对市场中的技术改进,或者使本技术领域的其它普通技术人员能理解本文披露的各实施例。
Claims (10)
1.一种模拟超临界CO2驱替天然气流动的方法,其特征在于,所述方法包括如下步骤:
根据气藏地质特征,建立三维地质模型;
通过改进状态方程建立流体模型确定气藏流体的基础参数;
确定CO2的扩散系数;
建立考虑扩散的超临界CO2驱天然气多组分数值模拟模型;
应用模型计算不同CO2注入量和天然气采出量的对应关系。
2.根据权利要求1所述的方法,其特征在于,确定气藏流体的基础参数,包括:实验测试不同CO2天然气混合物偏差因子、密度及粘度;
通过改进状态方程获得描述整个压力温度范围内临界性质的热力学参数。
3.根据权利要求1所述的方法,其特征在于,确定CO2的扩散系数包括:实验测量不同压力条件下CO2在天然气中的扩散浓度变化规律,应用插值计算方法描述多孔介质中CO2的扩散系数。
4.根据权利要求1所述的方法,其特征在于,建立考虑扩散的超临界CO2驱天然气多组分数值模拟模型包括:应用改进状态方程劈分计算各组分热力学参数,质量守恒方程中考虑扩散项对流动的影响。
5.根据权利要求1所述的方法,其特征在于,应用模型计算不同CO2注入量和天然气采出量的对应关系包括,根据长岩心驱替实验的物性参数,建立实验尺度数值模型,导入束缚水饱和度、渗透率、孔隙度参数,计算不同CO2注入量和天然气采出量的对应关系。
6.根据权利要求2所述的方法,其特征在于,引入体积修正项对状态方程进行修正,对临界点体积参数进行改进:
c(T)=ccf(Tr) (1)
其中,c(T)体积修正项,m3/mol;cc临界点改进体积修正项,m3/mol;f(Tr)温度改进系数。
7.根据权利要求6所述的方法,其特征在于,引入临界点改进项:
<mrow>
<msub>
<mi>c</mi>
<mi>c</mi>
</msub>
<mo>=</mo>
<msub>
<mi>v</mi>
<mrow>
<mi>E</mi>
<mi>O</mi>
<mi>S</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>v</mi>
<mrow>
<mi>E</mi>
<mi>x</mi>
<mi>p</mi>
</mrow>
</msub>
<mo>=</mo>
<mrow>
<mo>(</mo>
<mn>0.3074</mn>
<mo>-</mo>
<msub>
<mi>Z</mi>
<mi>c</mi>
</msub>
<mo>)</mo>
</mrow>
<mfrac>
<mrow>
<msub>
<mi>RT</mi>
<mi>c</mi>
</msub>
</mrow>
<msub>
<mi>P</mi>
<mi>c</mi>
</msub>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</mrow>
其中,临界压缩因子Zc为定值,R气体常数;Tc临界温度,K;Pc临界压力,MPa;vEOS状态方程计算摩尔体积,m3/mol;vExp实验计算摩尔体积,m3/mol;cc体积修正项,m3/mol。
8.根据权利要求6所述的方法,其特征在于,引入温度改进项进行修正,温度改进通过相平衡实验确定Zc与T的对应关系,通过实验拟合得方程式(3)-(5):
f(Tr)=β+(1-β)exp(0.5β+8.3) (3)
β=-2.8431exp([-64.2184(0.3074-Zc)]+0.1735) (4)
β=T-Tr (5)
其中,T气体所处的温度,K;Tr=T/Tc对比温度,K;β温度校正项,K。
9.根据权利要求7所述的方法,其特征在于,将临界点改进项和温度改进项带入PR状态方程:
<mrow>
<mi>P</mi>
<mo>=</mo>
<mfrac>
<mrow>
<mi>R</mi>
<mi>T</mi>
</mrow>
<mrow>
<mi>v</mi>
<mo>+</mo>
<mi>c</mi>
<mo>-</mo>
<mi>b</mi>
</mrow>
</mfrac>
<mo>-</mo>
<mfrac>
<mi>a</mi>
<mrow>
<msup>
<mrow>
<mo>(</mo>
<mi>v</mi>
<mo>+</mo>
<mi>c</mi>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>+</mo>
<mn>2</mn>
<mi>b</mi>
<mrow>
<mo>(</mo>
<mi>v</mi>
<mo>+</mo>
<mi>c</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msup>
<mi>b</mi>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>6</mn>
<mo>)</mo>
</mrow>
</mrow>
状态方程中,ab由临界参数和偏心因子计算,式中,
a=α(T)·ac;
<mrow>
<mi>&alpha;</mi>
<mrow>
<mo>(</mo>
<mi>T</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msup>
<mrow>
<mo>&lsqb;</mo>
<mn>1</mn>
<mo>+</mo>
<mi>m</mi>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<msubsup>
<mi>T</mi>
<mi>r</mi>
<mn>0.5</mn>
</msubsup>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>;</mo>
</mrow>
m=0.37464+1.54226ω-0.26992ω2;
Tr=T/Tc;
<mrow>
<msub>
<mi>a</mi>
<mi>c</mi>
</msub>
<mo>=</mo>
<mn>0.457235</mn>
<mo>&CenterDot;</mo>
<msup>
<mi>R</mi>
<mn>2</mn>
</msup>
<mo>&CenterDot;</mo>
<msubsup>
<mi>T</mi>
<mi>c</mi>
<mn>2</mn>
</msubsup>
<mo>/</mo>
<msub>
<mi>P</mi>
<mi>c</mi>
</msub>
<mo>;</mo>
</mrow>
b=0.77796·R·Tc/Pc
其中,p压力,MPa;v摩尔体积,m3/mol;ω偏心因子,无因次。
10.根据权利要求8所述的方法,其特征在于,考虑扩散的超临界CO2驱天然气多组分数值模拟模型包括,如下:
<mrow>
<mo>&dtri;</mo>
<mo>&lsqb;</mo>
<munder>
<mo>&Sigma;</mo>
<mi>g</mi>
</munder>
<mfrac>
<mrow>
<msub>
<mi>kk</mi>
<mrow>
<mi>r</mi>
<mi>g</mi>
</mrow>
</msub>
</mrow>
<msub>
<mi>&mu;</mi>
<mi>g</mi>
</msub>
</mfrac>
<msub>
<mi>&rho;</mi>
<mi>g</mi>
</msub>
<msub>
<mi>y</mi>
<mrow>
<mi>i</mi>
<mi>g</mi>
</mrow>
</msub>
<mo>&dtri;</mo>
<msubsup>
<mi>&Phi;</mi>
<mi>g</mi>
<mo>*</mo>
</msubsup>
<mo>&rsqb;</mo>
<mo>+</mo>
<mo>&dtri;</mo>
<mo>&lsqb;</mo>
<munder>
<mo>&Sigma;</mo>
<mi>g</mi>
</munder>
<msub>
<mi>&rho;</mi>
<mi>g</mi>
</msub>
<msub>
<mi>&phi;</mi>
<mi>g</mi>
</msub>
<msub>
<mi>S</mi>
<mi>g</mi>
</msub>
<msub>
<mi>D</mi>
<mrow>
<mi>i</mi>
<mi>g</mi>
</mrow>
</msub>
<mo>&dtri;</mo>
<msub>
<mi>y</mi>
<mrow>
<mi>i</mi>
<mi>g</mi>
</mrow>
</msub>
<mo>&rsqb;</mo>
<mo>+</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>k</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<msub>
<mi>n</mi>
<mi>f</mi>
</msub>
</munderover>
<munder>
<mo>&Sigma;</mo>
<mi>g</mi>
</munder>
<msub>
<mi>&rho;</mi>
<mi>i</mi>
</msub>
<msub>
<mi>q</mi>
<mrow>
<mi>i</mi>
<mi>k</mi>
</mrow>
</msub>
<msub>
<mi>y</mi>
<mrow>
<mi>i</mi>
<mi>g</mi>
</mrow>
</msub>
<mo>-</mo>
<mfrac>
<mrow>
<mo>&part;</mo>
<mo>&lsqb;</mo>
<munder>
<mo>&Sigma;</mo>
<mi>g</mi>
</munder>
<msub>
<mi>&phi;&rho;</mi>
<mi>g</mi>
</msub>
<msub>
<mi>y</mi>
<mrow>
<mi>i</mi>
<mi>g</mi>
</mrow>
</msub>
<msub>
<mi>S</mi>
<mi>g</mi>
</msub>
<mo>&rsqb;</mo>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>t</mi>
</mrow>
</mfrac>
<mo>=</mo>
<mn>0</mn>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>7</mn>
<mo>)</mo>
</mrow>
</mrow>
i=1,...,n,k=1,...,nf
其中,k渗透率,D;krg气相渗透率;μg气相粘度mPa·s;ρg气相密度kg/m3;yig组分i在气相中的摩尔分数;气相的势,MPa;φ孔隙度;Sg气相饱和度;qik组分i在k层中的产量或流量,m2/d;Dig为组分i在气相中的扩散系数,m2/d;nf为井的射孔层段;
扩散张量:
<mrow>
<msub>
<mi>D</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msubsup>
<mi>D</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
<mrow>
<mi>x</mi>
<mi>x</mi>
</mrow>
</msubsup>
</mtd>
<mtd>
<msubsup>
<mi>D</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
<mrow>
<mi>x</mi>
<mi>y</mi>
</mrow>
</msubsup>
</mtd>
<mtd>
<msubsup>
<mi>D</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
<mrow>
<mi>x</mi>
<mi>z</mi>
</mrow>
</msubsup>
</mtd>
</mtr>
<mtr>
<mtd>
<msubsup>
<mi>D</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
<mrow>
<mi>y</mi>
<mi>x</mi>
</mrow>
</msubsup>
</mtd>
<mtd>
<msubsup>
<mi>D</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
<mrow>
<mi>y</mi>
<mi>y</mi>
</mrow>
</msubsup>
</mtd>
<mtd>
<msubsup>
<mi>D</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
<mrow>
<mi>y</mi>
<mi>z</mi>
</mrow>
</msubsup>
</mtd>
</mtr>
<mtr>
<mtd>
<msubsup>
<mi>D</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
<mrow>
<mi>z</mi>
<mi>x</mi>
</mrow>
</msubsup>
</mtd>
<mtd>
<msubsup>
<mi>D</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
<mrow>
<mi>z</mi>
<mi>y</mi>
</mrow>
</msubsup>
</mtd>
<mtd>
<msubsup>
<mi>D</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
<mrow>
<mi>z</mi>
<mi>z</mi>
</mrow>
</msubsup>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>8</mn>
<mo>)</mo>
</mrow>
</mrow>
式中,
<mrow>
<msubsup>
<mi>D</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
<mrow>
<mi>y</mi>
<mi>y</mi>
</mrow>
</msubsup>
<mo>=</mo>
<mfrac>
<msubsup>
<mi>D</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
<mo>*</mo>
</msubsup>
<msub>
<mi>&tau;</mi>
<mi>j</mi>
</msub>
</mfrac>
<mo>+</mo>
<mfrac>
<mrow>
<mo>(</mo>
<msub>
<mi>&alpha;</mi>
<mrow>
<mi>j</mi>
<mi>l</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>&alpha;</mi>
<mrow>
<mi>j</mi>
<mi>t</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mrow>
<msub>
<mi>&phi;S</mi>
<mi>j</mi>
</msub>
<mo>|</mo>
<msub>
<mi>u</mi>
<mi>j</mi>
</msub>
<mo>|</mo>
</mrow>
</mfrac>
<msubsup>
<mi>u</mi>
<mrow>
<mi>j</mi>
<mi>y</mi>
</mrow>
<mn>2</mn>
</msubsup>
<mo>+</mo>
<mfrac>
<mrow>
<msub>
<mi>&alpha;</mi>
<mrow>
<mi>j</mi>
<mi>t</mi>
</mrow>
</msub>
<mo>|</mo>
<msub>
<mi>u</mi>
<mi>j</mi>
</msub>
<mo>|</mo>
</mrow>
<mrow>
<msub>
<mi>&phi;S</mi>
<mi>j</mi>
</msub>
</mrow>
</mfrac>
</mrow>
<mrow>
<msubsup>
<mi>D</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
<mrow>
<mi>z</mi>
<mi>z</mi>
</mrow>
</msubsup>
<mo>=</mo>
<mfrac>
<msubsup>
<mi>D</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
<mo>*</mo>
</msubsup>
<msub>
<mi>&tau;</mi>
<mi>j</mi>
</msub>
</mfrac>
<mo>+</mo>
<mfrac>
<mrow>
<mo>(</mo>
<msub>
<mi>&alpha;</mi>
<mrow>
<mi>j</mi>
<mi>l</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>&alpha;</mi>
<mrow>
<mi>j</mi>
<mi>t</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mrow>
<msub>
<mi>&phi;S</mi>
<mi>j</mi>
</msub>
<mo>|</mo>
<msub>
<mi>u</mi>
<mi>j</mi>
</msub>
<mo>|</mo>
</mrow>
</mfrac>
<msubsup>
<mi>u</mi>
<mrow>
<mi>j</mi>
<mi>z</mi>
</mrow>
<mn>2</mn>
</msubsup>
<mo>+</mo>
<mfrac>
<mrow>
<msub>
<mi>&alpha;</mi>
<mrow>
<mi>j</mi>
<mi>t</mi>
</mrow>
</msub>
<mo>|</mo>
<msub>
<mi>u</mi>
<mi>j</mi>
</msub>
<mo>|</mo>
</mrow>
<mrow>
<msub>
<mi>&phi;S</mi>
<mi>j</mi>
</msub>
</mrow>
</mfrac>
</mrow>
<mrow>
<msubsup>
<mi>D</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
<mrow>
<mi>x</mi>
<mi>y</mi>
</mrow>
</msubsup>
<mo>=</mo>
<msubsup>
<mi>D</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
<mrow>
<mi>y</mi>
<mi>x</mi>
</mrow>
</msubsup>
<mo>=</mo>
<mfrac>
<mrow>
<mo>(</mo>
<msub>
<mi>&alpha;</mi>
<mrow>
<mi>j</mi>
<mi>l</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>&alpha;</mi>
<mrow>
<mi>j</mi>
<mi>t</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mrow>
<msub>
<mi>&phi;S</mi>
<mi>j</mi>
</msub>
<mo>|</mo>
<msub>
<mi>u</mi>
<mi>j</mi>
</msub>
<mo>|</mo>
</mrow>
</mfrac>
<msub>
<mi>u</mi>
<mrow>
<mi>j</mi>
<mi>x</mi>
</mrow>
</msub>
<msub>
<mi>u</mi>
<mrow>
<mi>j</mi>
<mi>y</mi>
</mrow>
</msub>
</mrow>
<mrow>
<msubsup>
<mi>D</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
<mrow>
<mi>x</mi>
<mi>z</mi>
</mrow>
</msubsup>
<mo>=</mo>
<msubsup>
<mi>D</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
<mrow>
<mi>z</mi>
<mi>x</mi>
</mrow>
</msubsup>
<mo>=</mo>
<mfrac>
<mrow>
<mo>(</mo>
<msub>
<mi>&alpha;</mi>
<mrow>
<mi>j</mi>
<mi>l</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>&alpha;</mi>
<mrow>
<mi>j</mi>
<mi>t</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mrow>
<msub>
<mi>&phi;S</mi>
<mi>j</mi>
</msub>
<mo>|</mo>
<msub>
<mi>u</mi>
<mi>j</mi>
</msub>
<mo>|</mo>
</mrow>
</mfrac>
<msub>
<mi>u</mi>
<mrow>
<mi>j</mi>
<mi>x</mi>
</mrow>
</msub>
<msub>
<mi>u</mi>
<mrow>
<mi>j</mi>
<mi>z</mi>
</mrow>
</msub>
</mrow>
<mrow>
<msubsup>
<mi>D</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
<mrow>
<mi>z</mi>
<mi>y</mi>
</mrow>
</msubsup>
<mo>=</mo>
<msubsup>
<mi>D</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
<mrow>
<mi>y</mi>
<mi>z</mi>
</mrow>
</msubsup>
<mo>=</mo>
<mfrac>
<mrow>
<mo>(</mo>
<msub>
<mi>&alpha;</mi>
<mrow>
<mi>j</mi>
<mi>l</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>&alpha;</mi>
<mrow>
<mi>j</mi>
<mi>t</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mrow>
<msub>
<mi>&phi;S</mi>
<mi>j</mi>
</msub>
<mo>|</mo>
<msub>
<mi>u</mi>
<mi>j</mi>
</msub>
<mo>|</mo>
</mrow>
</mfrac>
<msub>
<mi>u</mi>
<mrow>
<mi>j</mi>
<mi>z</mi>
</mrow>
</msub>
<msub>
<mi>u</mi>
<mrow>
<mi>j</mi>
<mi>y</mi>
</mrow>
</msub>
</mrow>
<mrow>
<mo>|</mo>
<msub>
<mi>u</mi>
<mi>j</mi>
</msub>
<mo>|</mo>
<mo>=</mo>
<msup>
<mrow>
<mo>(</mo>
<msubsup>
<mi>u</mi>
<mrow>
<mi>j</mi>
<mi>x</mi>
</mrow>
<mn>2</mn>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>u</mi>
<mrow>
<mi>j</mi>
<mi>z</mi>
</mrow>
<mn>2</mn>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>u</mi>
<mrow>
<mi>j</mi>
<mi>y</mi>
</mrow>
<mn>2</mn>
</msubsup>
<mo>)</mo>
</mrow>
<mrow>
<mn>1</mn>
<mo>/</mo>
<mn>2</mn>
</mrow>
</msup>
</mrow>
其中,ρi组分i密度kg/m3;Sj j相饱和度;Dig为组分i在气相中的扩散系数,m2/d;为组分i在j相中的分子扩散系数,m2/d;uj为速度相量,m2/d;ujx为x方向速度相量,m2/d;ujy为y方向速度相量,m2/d;ujz为z方向速度相量,m2/d;τj迂曲率;αjl j相的沿程弥散系数;αjt j相的横向弥散系数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610562387.9A CN107621350B (zh) | 2016-07-15 | 2016-07-15 | 一种模拟超临界co2驱替天然气流动的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610562387.9A CN107621350B (zh) | 2016-07-15 | 2016-07-15 | 一种模拟超临界co2驱替天然气流动的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107621350A true CN107621350A (zh) | 2018-01-23 |
CN107621350B CN107621350B (zh) | 2019-11-05 |
Family
ID=61087989
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610562387.9A Active CN107621350B (zh) | 2016-07-15 | 2016-07-15 | 一种模拟超临界co2驱替天然气流动的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107621350B (zh) |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108460203A (zh) * | 2018-02-26 | 2018-08-28 | 中国地质大学(北京) | 页岩油藏循环注气油井产量分析方法及其应用 |
CN110298048A (zh) * | 2018-03-22 | 2019-10-01 | 中国石油化工股份有限公司 | 一种考虑界面相的超临界co2-凝析气多相渗流模拟方法 |
CN110318741A (zh) * | 2018-03-28 | 2019-10-11 | 中国石油化工股份有限公司 | 一种异常高压、高含水致密低渗气藏剩余气分布描述方法 |
CN111027211A (zh) * | 2019-12-10 | 2020-04-17 | 中国石油大学(华东) | 一种考虑孔喉受限机理的致密油藏流体模拟方法及在注气开发模拟的应用 |
CN111579463A (zh) * | 2020-06-29 | 2020-08-25 | 西南石油大学 | 有水气藏封存二氧化碳物理模拟装置及其模拟方法 |
CN113063703A (zh) * | 2021-03-29 | 2021-07-02 | 西南石油大学 | 含束缚水多孔介质中天然气与co2扩散系数的测试方法及装置 |
CN115081217A (zh) * | 2022-06-25 | 2022-09-20 | 西南石油大学 | 一种超临界条件下干气-凝析气动态扩散系数计算方法 |
CN116411896A (zh) * | 2023-06-07 | 2023-07-11 | 太原理工大学 | 超临界二氧化碳驱替开采煤层气的相态识别装置及方法 |
CN117828942A (zh) * | 2024-01-04 | 2024-04-05 | 成都理工大学 | 一种准平衡边界条件下co2在液相中扩散系数计算方法 |
CN118008225A (zh) * | 2024-03-25 | 2024-05-10 | 中国石油大学(北京) | 基于非完全混相特征的注二氧化碳开发优化方法及系统 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20050165593A1 (en) * | 2002-03-20 | 2005-07-28 | Roland Lenormand | Method for modelling the production of hydrocarbons by a subsurface deposit which are subject to depletion |
CN102943654A (zh) * | 2012-11-01 | 2013-02-27 | 中国海洋石油总公司 | 利用高含co2气藏提高高含烃气藏采收率的方法 |
CN104975827A (zh) * | 2014-04-03 | 2015-10-14 | 中国石油化工股份有限公司 | 预测二氧化碳驱油藏指标的物质平衡方法 |
-
2016
- 2016-07-15 CN CN201610562387.9A patent/CN107621350B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20050165593A1 (en) * | 2002-03-20 | 2005-07-28 | Roland Lenormand | Method for modelling the production of hydrocarbons by a subsurface deposit which are subject to depletion |
CN102943654A (zh) * | 2012-11-01 | 2013-02-27 | 中国海洋石油总公司 | 利用高含co2气藏提高高含烃气藏采收率的方法 |
CN104975827A (zh) * | 2014-04-03 | 2015-10-14 | 中国石油化工股份有限公司 | 预测二氧化碳驱油藏指标的物质平衡方法 |
Non-Patent Citations (2)
Title |
---|
许兰兵 等: "注入二氧化碳提高天然气采收率技术研究进展", 《内蒙古石油化式》 * |
郭平 等: "天然气藏CO2驱及地质埋存技术研究进展", 《岩性油气藏》 * |
Cited By (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108460203A (zh) * | 2018-02-26 | 2018-08-28 | 中国地质大学(北京) | 页岩油藏循环注气油井产量分析方法及其应用 |
CN110298048B (zh) * | 2018-03-22 | 2022-06-17 | 中国石油化工股份有限公司 | 一种考虑界面相的超临界co2-凝析气多相渗流模拟方法 |
CN110298048A (zh) * | 2018-03-22 | 2019-10-01 | 中国石油化工股份有限公司 | 一种考虑界面相的超临界co2-凝析气多相渗流模拟方法 |
CN110318741A (zh) * | 2018-03-28 | 2019-10-11 | 中国石油化工股份有限公司 | 一种异常高压、高含水致密低渗气藏剩余气分布描述方法 |
CN111027211A (zh) * | 2019-12-10 | 2020-04-17 | 中国石油大学(华东) | 一种考虑孔喉受限机理的致密油藏流体模拟方法及在注气开发模拟的应用 |
CN111027211B (zh) * | 2019-12-10 | 2020-09-08 | 中国石油大学(华东) | 一种考虑孔喉受限机理的致密油藏流体模拟方法及在注气开发模拟的应用 |
CN111579463A (zh) * | 2020-06-29 | 2020-08-25 | 西南石油大学 | 有水气藏封存二氧化碳物理模拟装置及其模拟方法 |
CN113063703A (zh) * | 2021-03-29 | 2021-07-02 | 西南石油大学 | 含束缚水多孔介质中天然气与co2扩散系数的测试方法及装置 |
CN115081217A (zh) * | 2022-06-25 | 2022-09-20 | 西南石油大学 | 一种超临界条件下干气-凝析气动态扩散系数计算方法 |
CN115081217B (zh) * | 2022-06-25 | 2024-08-16 | 西南石油大学 | 一种超临界条件下干气-凝析气动态扩散系数计算方法 |
CN116411896A (zh) * | 2023-06-07 | 2023-07-11 | 太原理工大学 | 超临界二氧化碳驱替开采煤层气的相态识别装置及方法 |
CN116411896B (zh) * | 2023-06-07 | 2023-08-18 | 太原理工大学 | 超临界二氧化碳驱替开采煤层气的相态识别装置及方法 |
CN117828942A (zh) * | 2024-01-04 | 2024-04-05 | 成都理工大学 | 一种准平衡边界条件下co2在液相中扩散系数计算方法 |
CN117828942B (zh) * | 2024-01-04 | 2024-09-24 | 成都理工大学 | 一种准平衡边界条件下co2在液相中扩散系数计算方法 |
CN118008225A (zh) * | 2024-03-25 | 2024-05-10 | 中国石油大学(北京) | 基于非完全混相特征的注二氧化碳开发优化方法及系统 |
Also Published As
Publication number | Publication date |
---|---|
CN107621350B (zh) | 2019-11-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107621350A (zh) | 一种模拟超临界co2驱替天然气流动的方法 | |
Sun et al. | Numerical study of CO2 enhanced natural gas recovery and sequestration in shale gas reservoirs | |
Khasanov et al. | Mathematical model of injection of liquid carbon dioxide in a reservoir saturated with methane and its hydrate | |
Gou et al. | Numerical simulation of carbon dioxide injection for enhanced gas recovery (CO 2-EGR) in Altmark natural gas field | |
CN107144889A (zh) | 一种基于等效孔隙理论的砂岩岩石物理建模方法 | |
Yanze et al. | The role of diffusion for nonequilibrium gas injection into a fractured reservoir | |
CN110298048B (zh) | 一种考虑界面相的超临界co2-凝析气多相渗流模拟方法 | |
Singh et al. | Gas-condensate pseudopressure in layered reservoirs | |
Wang et al. | Numerical simulation of CH4 recovery from gas hydrate using gaseous CO2 injected into porous media | |
Gao et al. | Pore-scale mechanisms and hysteresis effect during multi-cycle injection and production process in underground hydrogen storage reservoir | |
Bejestani et al. | Effect of petrophysical matrix properties on bypassed oil recovery from a matrix-fracture system during CO2 near-miscible injection: Experimental investigation | |
CN110318741A (zh) | 一种异常高压、高含水致密低渗气藏剩余气分布描述方法 | |
Zhang et al. | Study of gas-liquid two-phase flow characteristics in hydrate-bearing sediments | |
Male | Application of a one dimensional nonlinear model to flow in hydrofractured shale gas wells using scaling solutions | |
Huseby et al. | Use of natural geochemical tracers to improve reservoir simulation models | |
CN111287741B (zh) | 一种致密油藏体积压裂改造区渗透率的快速计算方法 | |
Liu et al. | Study on the two-component gas apparent permeability model in shale nanopores considering Knudsen number correction | |
Salimi et al. | Compositional flow simulations of mixed CO2-water injection into geothermal reservoirs: Geothermal energy combined with CO2 storage | |
Zhang et al. | The development of a new reservoir simulator to model polymer flooding and advanced wells for enhanced heavy oil recovery in Bohai bay | |
Glavnov et al. | EOR miscible gas project in oil-gas condensate field | |
Wang et al. | A two-phase, multi-component model for efficient CO2 storage and enhanced gas recovery in low permeability reservoirs | |
Kuk et al. | Estimation of carbon dioxide injection into reservoirs potential to enhance oil production and reduce CO2 emissions to the atmosphere | |
Obidi et al. | Analytical solution for compositional profile driven by gravitational segregation and diffusion | |
Syzdykov | Evaluating gas injection performance in very low permeability, thick carbonate gas-condensate reservoirs to improve ultimate liquid yield | |
Fu et al. | Sensitivity Analysis of Diffusion-Based Hydrocarbon Gas Huff-n-Puff Injection in Eagle Ford Shale |
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 |