IPC分类号 : G06F30/20,G06F17/11,G01N5/00,G06F113/08
专利摘要
本发明公开一种天然气解析过程中碳同位素分馏计算方法及应用,属于非常规天然气开发领域,包括:计算初始时刻游离气中12CH4和13CH4质量;计算第一时间间隔流出气体中12CH4和13CH4质量;算出解吸气12CH4和13CH4质量;结合初始时刻游离气12CH4和13CH4质量、流出气体及解吸气12CH4和13CH4质量,计算末尾时刻游离气中12CH4和13CH4质量及碳同位素值;经历下一个时间间隔,以第一时间间隔末尾时刻各物理量作为第二时间间隔的初始量,重复上述过程,迭代计算至岩石内压力等于外界压力。该方法能够合理表征天然气开发中碳同位素分馏全过程。
权利要求
1.一种天然气解析过程中碳同位素分馏计算方法,其特征在于,包括:
(1)获取初始时刻初始条件下的第一游离气质量和吸附气质量,计算所述第一游离气和所述吸附气的碳同位素初始值,并计算所述第一游离气中
(2)待开采开始时,根据质量守恒方程计算第一时间间隔内流出地层的流出气体质量和所述第一时间间隔末尾时刻的地层压力,再根据所述流出气体质量与所述初始时刻时所述第一游离气质量之间的关系计算所述第一时间间隔内所述流出气体的碳同位素值,并计算所述流出气体中
(3)计算所述吸附气含量、解吸气质量以及所述吸附气和所述解吸气的碳同位素值,并确定所述解吸气中
(4)根据所述初始时刻时所述第一游离气中
(5)待开采过程进入相邻的第二时间间隔时,以所述第一时间间隔末尾时刻时所述第一游离气中
2.根据权利 要求1所述的天然气解析过程中碳同位素分馏计算方法,其特征在于,所述初始时刻初始条件下的初始碳同位素值、所述第一游离气的碳同位素值和所述吸附气的碳同位素值之间满足的关系式为:
其中,m
3.根据权利 要求1所述的天然气解析过程中碳同位素分馏计算方法,其特征在于,所述初始时刻初始条件下的所述第一游离气中
其中,
4.根据权利 要求1所述的天然气解析过程中碳同位素分馏计算方法,其特征在于,所述计算所述流出气体中
获取所述初始时刻初始条件下所述第一游离气的质量、碳同位素初始值以及质量守恒方程;
根据所述质量守恒方程计算所述第一时间间隔内所述第一游离气的流出质量和地层压力;
根据所述第一游离气的含量和碳同位素初始值以及所述第一时间间隔内所述第一游离气的流出质量计算所述第一时间间隔内所述流出气体的碳同位素值;
根据所述第一时间间隔内所述第一游离气的流出质量和所述流出气体的碳同位素值计算所述第一时间间隔内所述流出气体中
5.根据权利 要求4所述的天然气解析过程中碳同位素分馏计算方法,其特征在于,所述质量守恒方程为:
其中,m
6.根据权利 要求4所述的天然气解析过程中碳同位素分馏计算方法,其特征在于,所述第一时间间隔内所述流出气体的碳同位素值的计算公式为:
其中,δ
7.根据权利 要求4所述的天然气解析过程中碳同位素分馏计算方法,其特征在于,所述第一时间间隔内所述流出气体中
其中,m
8.根据权利 要求1所述的天然气解析过程中碳同位素分馏计算方法,其特征在于,所述解吸气中
待开采开始时,地层压力从初始压力p
根据所述第一时间间隔末尾时刻所述解吸气的碳同位素值
其中,所述第一时间间隔内所述解吸气质量m
其中,ρ
所述第一时间间隔末尾时刻所述解吸气的碳同位素值
所述第一时间间隔内所述解吸气中
其中,m
9.根据权利 要求1所述的天然气解析过程中碳同位素分馏计算方法,其特征在于,所述第一时间间隔末尾时刻所述第一游离气中
其中,
10.如权利 要求1-9中任一项所述的计算方法在构建页岩气生产过程气体碳同位素演化模型中的应用。
说明书
技术领域
本发明属于非常规天然气开发技术领域,涉及一种天然气解析过程中碳同位素分馏计算方法及应用。
背景技术
随着非常规天然气藏的商业开发,生产过程中的气体碳同位素变化引起了越来越多的关注,相对于常规气藏稳定的气体碳同位素值,非常规气藏生产过程中的碳同位素值出现了不同的演化特征,比如页岩气藏碳同位素值具有先变轻后变重的特点。
前人研究证实,复杂孔隙结构中气体流动产生的碳同位素分馏效应存在扩散作用和吸附-解吸作用两种机理。Xia和Tang(Xia Xinyu&Tang Yongchun.Isotopefractionation of methane dur ing natural gas flow with coupled diffusion andadsorption/desorption.Geochimica et Cosmochimica Acta,2012,77,489-503)提出同时考虑吸附-解吸作用和扩散作用的连续流动模型定量评价天然气中甲烷从烃源岩中流出的同位素分馏效应。在此基础上,Cao等(Cao Gaohui,Zhang Huanxu,Jiang Wenbin, WuSheng,Zhu Di,Lin Mian.A new gas-content-evaluation method for organic-richshale using the fractionation of carbon isotopes of methane.SPE Journal,2019,2574-2589)同时考虑有机质和无机质孔隙中气体流动及赋存形式,建立了甲烷同位素分馏模型模拟气体从页岩岩屑中流出的碳同位素分馏效应。李文镖等(李文镖,卢双舫,李俊乾,张鹏飞,王思远,冯文俊,魏永波.页岩气运移过程中的碳同位素分馏:机理、表征及意义.中国科学:地球科学,2020,50,553-569)认为吸附-解吸作用是引起同位素分馏的主要机理,假设完整的解析过程符合一系列平行反应组成的动力学模型,该模型适用于表征吸附气的解吸过程,但无法描述游离气的渗流过程,因此尚不能用于表征同位素分馏全过程。
并且,以上模型都有一个共同的特点:将甲烷视为由
发明内容
为了解决现有模型无法合理表征同位素分馏全过程的技术问题,本发明提供一种天然气解析过程中碳同位素分馏计算方法。该方法能够表征天然气开发中碳同位素分馏全过程,可用于模拟页岩气、煤层气和致密气等以甲烷为主体的非常规气藏在生产过程中气体碳同位素的演化。
本发明还提供一种天然气解析过程中碳同位素分馏计算方法在构建页岩气或煤层气或致密气生产过程气体碳同位素演化模型中的应用。
本发明通过以下技术方案实现:
一种天然气解析过程中碳同位素分馏计算方法,包括:
(1)获取初始时刻初始条件下的第一游离气质量和吸附气质量,计算所述第一游离气和所述吸附气的碳同位素初始值,并计算所述第一游离气中
(2)待开采开始时,根据质量守恒方程计算第一时间间隔内流出地层的流出气体质量和所述第一时间间隔末尾时刻的地层压力,再根据所述流出气体质量与所述初始时刻时所述第一游离气质量之间的关系计算所述第一时间间隔内所述流出气体的碳同位素值,并计算所述流出气体中
(3)计算所述吸附气含量、解吸气质量以及所述吸附气和所述解吸气的碳同位素值,并确定所述解吸气中
(4)根据所述初始时刻时所述第一游离气中
(5)待开采过程进入相邻的第二时间间隔时,以所述第一时间间隔末尾时刻时所述第一游离气中
进一步的,所述初始时刻初始条件下的初始碳同位素值、所述第一游离气的碳同位素值和所述吸附气的碳同位素值之间满足的关系式为:
其中,m
进一步的,所述初始时刻初始条件下的所述第一游离气中
其中, 为所述初始时刻初始条件下所述第一游离气的含量, 为所述初始时刻初始条件下所述第一游离气中
进一步的,所述计算所述流出气体中
获取所述初始时刻初始条件下所述第一游离气的质量、碳同位素初始值以及质量守恒方程;
根据所述质量守恒方程计算所述第一时间间隔内所述第一游离气的流出质量和地层压力;
根据所述第一游离气的含量和碳同位素初始值以及所述第一时间间隔内所述第一游离气的流出质量计算所述第一时间间隔内所述流出气体的碳同位素值;
根据所述第一时间间隔内所述第一游离气的流出质量和所述流出气体的碳同位素值计算所述第一时间间隔内所述流出气体中
进一步的,所述质量守恒方程为:
其中,mf表示所述第一游离气质量,t表示时间,r表示岩石半径,m表示岩石形状因子,ρ表示气体密度,μ表示气体粘度,k表示气体渗透率,p表示地层压力。
进一步的,所述第一时间间隔内所述流出气体的碳同位素值的计算公式为:
其中,δ
进一步的,所述第一时间间隔内所述流出气体中
其中,mins(12)为所述第一时间间隔内所述流出气体中
进一步的,所述解吸气中
待开采开始时,地层压力从初始压力p0下降为pt,计算所述第一时间间隔内所述解吸气质量md;
根据所述第一时间间隔末尾时刻所述解吸气的碳同位素值 所述第一时间间隔末尾时刻所述吸附气的碳同位素值 和所述初始时刻初始条件下所述吸附气的碳同位素值 之间关系计算所述第一时间间隔末尾时刻所述吸附气的碳同位素值 和所述第一时间间隔末尾时刻所述解吸气的碳同位素值
其中,所述第一时间间隔内所述解吸气质量md计算公式为:
其中,ρstd为甲烷标况下气体密度,VL和PL分别为岩石的Langmuir体积和压力,mr为岩石质量;
所述第一时间间隔末尾时刻所述解吸气的碳同位素值 所述第一时间间隔末尾时刻所述吸附气的碳同位素值 和所述初始时刻初始条件下所述吸附气的碳同位素值 之间满足公式为:
所述第一时间间隔内所述解吸气中
其中,md(12)为所述第一时间间隔内所述解吸气中
进一步的,所述第一时间间隔末尾时刻所述第一游离气中
其中, 为所述第一时间间隔末尾时刻所述游离气中
一种天然气解析过程中碳同位素分馏计算方法在构建页岩气生产过程气体碳同位素演化模型中的应用。
本发明中,初始时刻初始条件下的原地含气量和游离气与吸附气比值可根据CN110687006A专利所公开的计算方法求得,具体过程在此不再详述。
本发明计算过程的方程中,a为吸附气与游离气之间的同位素分馏系数,与地层温度有关, T为地层温度,通过钻井数据获得。
本发明实施例中的一个或多个技术方案,至少具有如下技术效果或优点:
本发明一种天然气解析过程中碳同位素分馏计算方法,该方法先计算出初始时刻初始压力条件下游离气
附图说明
为了更清楚地说明本发明实施例中的技术方案,下面将对实施例描述中所需要使用的附图作一简单地介绍,显而易见地,下面描述中的附图是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其它的附图。
图1是本发明碳同位素分馏计算方法的流程图。
图2是本发明实施例1吸附气和游离气随时间的变化过程图。
图3是本发明实施例1岩芯流出甲烷气体的碳同位素值变化过程图。
具体实施方式
下文将结合具体实施方式和实施例,具体阐述本发明,本发明的优点和各种效果将由此更加清楚地呈现。本领域技术人员应理解,这些具体实施方式和实施例是用于说明本发明,而非限制本发明。
在整个说明书中,除非另有特别说明,本文使用的术语应理解为如本领域中通常所使用的含义。因此,除非另有定义,本文使用的所有技术和科学术语具有与本发明所属领域技术人员的一般理解相同的含义。若存在矛盾,本说明书优先。
除非另有特别说明,本发明中用到的各种原材料、试剂、仪器和设备等,均可通过市场购买得到或者可通过现有方法制备得到。
本发明实施例提供的技术方案为解决上述技术问题,总体思路如下:
现有天然气生产碳同位素分馏模型的设计通常忽略了
本发明首先将甲烷气体作为单一气体,分别计算基质内游离气和吸附气在生产(压力下降)过程中的减少,游离气的减少符合连续流动方程,碳同位素分馏效应遵循扩散机理,而吸附气的减少符合朗格谬尔方程,且碳同位素分馏效应遵循吸附-解吸机理。然后,根据碳同位素分馏效应和游离气/吸附气的减少质量,逐步计算吸附气、解吸气、游离气和流出气的碳同位素值变化。
下面将结合实施例对本申请一种天然气解析过程中碳同位素分馏计算方法进行详细说明。
实施例1
一种天然气解析过程中碳同位素分馏计算方法,包括:
(1)获取初始时刻初始条件下的第一游离气质量和吸附气质量,计算所述第一游离气和所述吸附气的碳同位素初始值,并计算所述第一游离气中
(2)待开采开始时,根据质量守恒方程计算第一时间间隔内流出地层的流出气体质量和所述第一时间间隔末尾时刻的地层压力,再根据所述流出气体质量与所述初始时刻时所述第一游离气质量之间的关系计算所述第一时间间隔内所述流出气体的碳同位素值,并计算所述流出气体中
(3)计算所述吸附气含量、解吸气质量以及所述吸附气和所述解吸气的碳同位素值,并确定所述解吸气中
(4)根据所述初始时刻时所述第一游离气中
(5)待开采过程进入相邻的第二时间间隔时,以所述第一时间间隔末尾时刻时所述第一游离气中
进一步的,所述初始时刻初始条件下的初始碳同位素值、所述第一游离气的碳同位素值和所述吸附气的碳同位素值之间满足的关系式为:
其中,m
进一步的,所述初始时刻初始条件下的所述第一游离气中
其中, 为所述初始时刻初始条件下所述第一游离气的含量, 为所述初始时刻初始条件下所述第一游离气中
进一步的,所述计算所述流出气体中
获取所述初始时刻初始条件下所述第一游离气的质量、碳同位素初始值以及质量守恒方程;
根据所述质量守恒方程计算所述第一时间间隔内所述第一游离气的流出质量和地层压力;
根据所述第一游离气的含量和碳同位素初始值以及所述第一时间间隔内所述第一游离气的流出质量计算所述第一时间间隔内所述流出气体的碳同位素值;
根据所述第一时间间隔内所述第一游离气的流出质量和所述流出气体的碳同位素值计算所述第一时间间隔内所述流出气体中
进一步的,所述质量守恒方程为:
其中,mf表示所述第一游离气质量,t表示时间,r表示岩石半径,m表示岩石形状因子,ρ表示气体密度,μ表示气体粘度,k表示气体渗透率,p表示地层压力。
进一步的,所述第一时间间隔内所述流出气体的碳同位素值的计算公式为:
其中,δ
进一步的,所述第一时间间隔内所述流出气体中
其中,mins(12)为所述第一时间间隔内所述流出气体中
进一步的,所述解吸气中
待开采开始时,地层压力从初始压力p0下降为pt,计算所述第一时间间隔内所述解吸气质量md;
根据所述第一时间间隔末尾时刻所述解吸气的碳同位素值 所述第一时间间隔末尾时刻所述吸附气的碳同位素值 和所述初始时刻初始条件下所述吸附气的碳同位素值 之间关系计算所述第一时间间隔末尾时刻所述吸附气的碳同位素值 和所述第一时间间隔末尾时刻所述解吸气的碳同位素值
其中,所述第一时间间隔内所述解吸气质量md计算公式为:
其中,ρstd为甲烷标况下气体密度,VL和PL分别为岩石的Langmuir体积和压力,mr为岩石质量;
所述第一时间间隔末尾时刻所述解吸气的碳同位素值 所述第一时间间隔末尾时刻所述吸附气的碳同位素值 和所述初始时刻初始条件下所述吸附气的碳同位素值 之间满足公式为:
所述第一时间间隔内所述解吸气中
其中,md(12)为所述第一时间间隔内所述解吸气中
进一步的,所述第一时间间隔末尾时刻所述第一游离气中
其中, 为所述第一时间间隔末尾时刻所述游离气中
一种天然气解析过程中碳同位素分馏计算方法及应用专利购买费用说明
Q:办理专利转让的流程及所需资料
A:专利权人变更需要办理著录项目变更手续,有代理机构的,变更手续应当由代理机构办理。
1:专利变更应当使用专利局统一制作的“著录项目变更申报书”提出。
2:按规定缴纳著录项目变更手续费。
3:同时提交相关证明文件原件。
4:专利权转移的,变更后的专利权人委托新专利代理机构的,应当提交变更后的全体专利申请人签字或者盖章的委托书。
Q:专利著录项目变更费用如何缴交
A:(1)直接到国家知识产权局受理大厅收费窗口缴纳,(2)通过代办处缴纳,(3)通过邮局或者银行汇款,更多缴纳方式
Q:专利转让变更,多久能出结果
A:著录项目变更请求书递交后,一般1-2个月左右就会收到通知,国家知识产权局会下达《转让手续合格通知书》。
动态评分
0.0