CN109188558B - 一种应用生烃增压模型计算致密油充注距离的方法 - Google Patents
一种应用生烃增压模型计算致密油充注距离的方法 Download PDFInfo
- Publication number
- CN109188558B CN109188558B CN201811155522.3A CN201811155522A CN109188558B CN 109188558 B CN109188558 B CN 109188558B CN 201811155522 A CN201811155522 A CN 201811155522A CN 109188558 B CN109188558 B CN 109188558B
- Authority
- CN
- China
- Prior art keywords
- hydrocarbon
- volume
- rock
- kerogen
- hydrocarbon generation
- 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
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V9/00—Prospecting or detecting by methods not provided for in groups G01V1/00 - G01V8/00
Abstract
本发明提供了一种应用生烃增压模型计算致密油充注距离的方法。所述方法包括如下步骤:(1)构建基于干酪根与生成油气密度比的生烃增压模型,生烃过程所产生的源储压差作为油气充注的主要驱动力;(2)以步骤(1)构建的生烃增压模型为基础,建立充注距离与剩余源储压差之间的计算公式,并根据建立的充注距离计算公式计算得到所述致密油充注距离。本发明方法相对简便,参数易于获取。
Description
技术领域
本发明涉及石油和天然气地质勘探技术领域,具体的说,本发明涉及一种应用生烃增压模型计算致密油充注距离的方法。
背景技术
致密油是指以吸附或游离状态赋存于生油岩中,或与生油岩互层、紧邻的致密砂岩、致密碳酸盐岩等储集岩中,未经过大规模长距离运移的非常规石油聚集。致密油从生油层排出之后经过充注过程进入储集层中,致密储层中缺乏可动水,油气运移的主要驱动力不是浮力,而是烃源岩中生烃过程压力增加(即生烃增压)形成的源储压差。
郭小文等(2011)研究了生油增压定量模型,讨论了生油增压受有机质丰度(TOC)、类型(原始氢指数,HI)、石油残留系数(α)的影响因素。但在郭小文等(2011)的模型中,没有建立生烃增压与油气运移的对应关系,不能应用于运移和充注距离计算。目前还没有关于生烃增压与致密油充注距离方面的研究报道。张焕旭等(2018)介绍了国内外生烃增压研究现状,强调了生烃增压产生的前提条件,同时也提到了生烃增压可作为,生烃增压是推动油气在致密储层中运移的主要动力。但在张焕旭等(2018)的综述中,并未开展生烃增压与油气充注方面的论述。
发明内容
本发明的目的在于提供一种应用生烃增压模型计算致密油充注距离的方法,以确定致密油的充注范围,为致密油分布预测及甜点区选择提供可靠依据。
为达上述目的,本发明提供了一种应用生烃增压模型计算致密油充注距离的方法,其中,所述方法包括如下步骤:
(1)构建基于干酪根与生成油气密度比的生烃增压模型,生烃过程所产生的源储压差作为油气充注的主要驱动力;
(2)以步骤(1)构建的生烃增压模型为基础,建立充注距离与剩余源储压差之间的计算公式,并根据建立的充注距离计算公式计算得到所述致密油充注距离。
根据本发明一些具体实施方案,其中,所述方法包括如下步骤:
(1)在设定条件为没有发生油气排出和充注前提下构建生烃增压模型,建立生烃增压与孔隙水初始体积Vw1、干酪根初始体积Vk1、生烃转化率F、生成油气密度ρo之间的数学关系式;
(2)以步骤(1)构建的生烃增压模型为基础,在设定条件为发生油气充注的主要动力为生烃增压前提下,建立油气充注之后的剩余源储压差ΔP、烃源岩孔隙水初始体积Vw1、干酪根初始体积Vk1、生烃转化率F、生成油气密度ρo、储层孔隙度之间的关系式,从而计算得到所述致密油充注距离。
根据本发明一些具体实施方案,其中,步骤(1)中设定条件还包括生烃前后体积相等:Vr1+Vw1+Vk1=Vr2+Vw2+Vk2+Vo2;其中Vr1,Vw1,Vk1分别为生烃过程发生前,烃源岩中岩石骨架体积、地层水体积和干酪根体积,单位均为cm3;Vr2,Vw2,Vk2,Vo2分别为生烃过程发生后,烃源岩中岩石骨架体积、地层水体积、剩余的干酪根体积以及新生成的油气体积,单位均为cm3。
根据本发明一些具体实施方案,其中,Vk1=ρr÷ρk×TOC×HI÷100000;Vt1=1-Vk1-Vw1。其中,ρr和ρk分别为岩石密度和干酪根密度,单位g/cm3;TOC为岩石总有机碳含量,单位%;HI为烃源岩原始氢指数,单位‰;为烃源岩的初始孔隙度,单位%。
根据本发明一些具体实施方案,其中,Vr2=Vr1×(1-Cr×ΔP);Vw2=Vw1×(1-Cw×ΔP);Vk2=(1-F)×Vk1×(1-Ck×ΔP);Vo2=(ρk÷ρo)×F×Vk1×(1-Co×ΔP);
其中Cr、Cw、Ck和Co分别为岩石、地层水、干酪根、石油的压缩系数,单位为10-3MPa-1;F为生烃转化率,取值范围为0-1;ΔP为油气充注之后的剩余源储压差(即生烃后与生烃前的压力差)。
根据本发明一些具体实施方案,其中,Cr、Cw、Ck和Co分别取值为0.2、0.44、1.4和2.2。
根据本发明一些具体实施方案,其中,ΔP=P2-P1;
其中,P1为生烃过程发生前烃源岩中的流体压力,P2为生烃过程发生后烃源岩中的流体压力,单位MPa。
根据本发明一些具体实施方案,其中,步骤(1)构建的生烃增压模型由如下式(1)表示:
根据本发明一些具体实施方案,其中,步骤(2)建立的充注距离计算公式由如下式(2)表示:
根据本发明一些具体实施方案,其中,本发明具体包括如下步骤:
(1)在设定条件为没有发生油气充注前提下构建生烃增压模型,建立生烃增压与孔隙水初始体积Vw1、干酪根初始体积Vk1、生烃转化率F以及生成油气密度ρo之间的数学关系式;
(2)以步骤(1)重新构建的生烃增压模型为基础,建立充注距离与剩余源储压差ΔP、烃源岩孔隙水初始体积Vw1、干酪根初始体积Vk1、生烃转化率F、生成油气密度ρo、储层孔隙体积度之间的数学关系式,从而计算得到所述致密油充注距离。
在一个实施例中,生烃过程发生前,烃源岩中包括岩石骨架、地层水和干酪根三部分,三者体积分别为Vt1,Vw1,Vk1,烃源岩中的流体压力为P1;干酪根生烃过程发生后,烃源岩中包括岩石骨架、地层水、剩余的干酪根以及新生成的油气四部分,体积分别为Vr2,Vw2,Vk2,Vo2,烃源岩中的流体压力为P2。生烃后与生烃前的压力差ΔP=P2-P1。干酪根向油气的转化率为F,取值范围是0~1;干酪根的密度为ρk,取值1.4g/cm3;生成的油气密度为ρo,可通过实际样品测定;岩石、地层水、干酪根、石油的压缩系数分别为Cr、Cw、Ck和Co,单位是10- 3MPa-1,分别取值为0.2、0.44、1.4和2.2。
在一个实施例中,首先假定没有发生油气充注,即生成的油气全部滞留在烃源岩中(如图1所示)。生烃前后体积相等,即Vr1+Vw1+Vk1=Vr2+Vw2+Vk2+Vo2。据此可通过以下公式构建生烃增压数学模型:
此压差即为烃源岩和邻近储集层之间的压力差,简称源储压差。
在一个实施例中,生成的油气并没有完全滞留在烃源岩中,而是一部分充注进入了储集层中,充注距离为L2(如图2所示)。由于油气的充注进入储集层中,烃源岩和接受油气充注的储集层段连通,即油气保存的空间增大。据此,可建立充注距离的计算公式:
在一个实施例中,随着烃源岩埋藏深度的增加,干酪根转化率(F)增加,生成的油气密度(ρo)降低,相应的充注距离增加。
综上所述,本发明提供了一种应用生烃增压模型计算充注距离的方法。本发明的方法具有如下优点:
一是本发明提供了一种致密油开采甜点的预测方法,通过致密油充注距离的计算,为有利勘探目标区优选提供了科学依据。同时,本发明还可进行烃源岩和邻近储集层地层压力预测,为致密油开发提供有效参数。
二是本发明相对简便,参数易于获取,只需要获得烃源岩岩石总有机碳含量(TOC)、烃源岩原始氢指数(HI)、生成的油气密度(ρo)、烃源岩生烃转化率(F)、烃源岩中的地层压力(P2)、烃源岩孔隙度储层孔隙度等,这些都是在地质研究中易于获取的参数。
三是本发明可以计算任意转化率条件下的致密油的充注距离。烃源岩在地质条件下的生烃转化是一个动态过程,只要调整烃源岩转化率(F)、生成油气密度(ρo)两项参数,就可以及时不同成熟度烃源岩的充注距离。
附图说明
图1为没有油气充注的生烃增压模型;
图2为油气充注距离示意图。
具体实施方式
以下通过具体实施例详细说明本发明的实施过程和产生的有益效果,旨在帮助阅读者更好地理解本发明的实质和特点,不作为对本案可实施范围的限定。
实施例1
生烃过程发生前,烃源岩中包括岩石骨架、地层水和干酪根三部分,三者体积分别为Vr1,Vw1,Vk1,烃源岩中的流体压力为P1;干酪根生烃过程发生后,烃源岩中包括岩石骨架、地层水、剩余的干酪根以及新生成的油气四部分,体积分别为Vr2,Vw2,Vk2,Vo2,烃源岩中的流体压力为P2。生烃后与生烃前的压力差ΔP=P2-P1,此压差亦即油气充注之后的剩余源储压差。干酪根向油气的转化率为F,取值范围是0~1;烃源岩的密度为ρr,取值2.5g/cm3;干酪根的密度为ρk,取值1.4g/cm3;生成的烃类密度为ρo,可通过实际样品测定;岩石、地层水、干酪根、石油的压缩系数分别为Cr、Cw、Ck和Co,单位是10-3MPa-1,可分别取值为0.2、0.44、1.4和2.2。
其中,Vk1=ρr÷ρk×TOC×HI÷100000,TOC为总有机碳含量(单位:%),HI为烃源岩原始氢指数(单位:‰),这两项参数通过烃源岩岩石热解和碳硫分析仪测试获取,分析方法见标准(GB/T 18602-2012和GBT 19145-2003);为烃源岩的初始孔隙度(单位:%);Vr1=1-Vk1-Vw1。
Vr2=Vr1×(1-Cr×ΔP),Vw2=Vw1×(1-Cw×ΔP),Vk2=(1-F)×Vk1×(1-Ck×ΔP),Vo2=D×F×Vk1×(1-Co×ΔP)。
第一步,假定烃源岩是一个封闭体系,即没有油气排出,生成的油气全部滞留在烃源岩中(如图1所示)。生烃前后体积相等,即Vr1+Vw1+Vk1=Vr2+Vw2+Vk2+Vo2。据此可通过以下公式构建生烃增压模型:
在该模型中,烃源岩总有机碳含量(TOC)、原始氢指数(HI)、烃源岩初始孔隙度干酪根生烃转化率(F)、生成的油气密度(ρo)是输入参数;干酪根密度、岩石、地层水、干酪根、石油的压缩系数是常数;生烃增压量ΔP是输出结果。
通过敏感性分析表明,生烃增压与干酪根初始体积(Vk1,用TOC、HI来表示,Vk1=ρr÷ρk×TOC×HI÷100000)正相关,干酪根初始体积越大,生烃增压越高;生烃增压与干酪根转化率(F)正相关,干酪根转化率越高,生烃增压越高;生烃增压与生成油气密度(ρo)负相关,油气密度越小,生烃增压越高。生烃增压与烃源岩孔隙水初始体积(Vw1)负相关,水体积越多,生烃增压也高。
此时,由于没有烃类充注,源储压差达到最大值。
第二步,生成的油气并没有完全滞留在烃源岩中,而是一部分充注进入了储集层中,充注距离为L2。由于油气的充注进入储集层中,烃源岩和接受油气充注的储集层段连通,即油气保存的空间增大。据此,可建立充注距离的计算公式:
在该模型中,烃源岩总有机碳含量(TOC)、原始氢指数(HI)、烃源岩初始孔隙度干酪根生烃转化率(F)、生成的油气密度(ρo)、储层孔隙度剩余源储压差(ΔP)、烃源岩厚度(L1)是输入参数;干酪根密度、岩石、地层水、干酪根、石油的压缩系数是常数;油气充注距离L2是输出结果。
通过敏感性分析表明,充注距离与干酪根初始体积(Vk1)正相关,干酪根初始体积越大,充注距离越大;与干酪根转化率(F)正相关,干酪根转化率越高,充注距离越大;与烃源岩厚度(L1)正相关,厚度越大,充注距离越大;与生成油气密度(ρo)负相关,油气密度越小,充注距离越大;与剩余源储压差(ΔP)负相关,剩余源储压差越小,充注距离越大;与储层孔隙度负相关,储层孔隙度越小,充注距离越大。
第三步,建立烃源岩热演化过程中的生烃增压与充注距离之间的关系。生烃增压与充注距离符合第二步建立的数学模型。输入参数中,烃源岩总有机碳含量(TOC)、原始氢指数(HI)、烃源岩初始孔隙度储层孔隙度为静态参数;而干酪根生烃转化率(F)、生成的油气密度(ρo)为动态参数,与有机质成熟度(Ro)和生成的烃类成分有关。
随着烃源岩埋藏深度的增加,烃源岩成熟度(Ro)增加,生成的油气成分变轻,生成的油气密度(ρo)降低,充注距离(L2)也会逐渐增加。具体计算过程中,需根据实测Ro值,确定烃源岩转化率、生成油气密度(ρo)等参数的变化情况,同时结合实际地质情况选择相应的源储压差(ΔP),即可计算油气的充注距离(L2)。
Claims (1)
1.一种应用生烃增压模型计算致密油充注距离的方法,其中,所述方法包括如下步骤:
(1)在设定条件为没有发生油气排出和充注前提下构建生烃增压模型,建立生烃增压与孔隙水初始体积Vw1、干酪根初始体积Vk1、生烃转化率F、生成油气密度ρo之间的数学关系式,生烃过程所产生的源储压差作为油气充注的主要驱动力;生烃前后体积相等:Vr1+Vw1+Vk1=Vr2+Vw2+Vk2+Vo2;其中Vr1,Vw1,Vk1分别为生烃过程发生前,烃源岩中岩石骨架体积、地层水体积和干酪根体积,单位均为cm3;Vr2,Vw2,Vk2,Vo2分别为生烃过程发生后,烃源岩中岩石骨架体积、地层水体积、剩余的干酪根体积以及新生成的油气体积,单位均为cm3;Vk1=ρr÷ρk×TOC×HI÷100000;Vr1=1-Vk1-Vw1;其中,ρr和ρk分别为岩石密度和干酪根密度,单位g/cm3;TOC为岩石总有机碳含量,单位%;HI为烃源岩原始氢指数,单位‰;为烃源岩的初始孔隙度,单位%;Vr2=Vr1×(1-Cr×ΔP);Vw2=Vw1×(1-Cw×ΔP);Vk2=(1-F)×Vk1×(1-Ck×ΔP);Vo2=(ρk÷ρo)×F×Vk1×(1-Co×ΔP);其中,Cr、Cw、Ck和Co分别为岩石、地层水、干酪根、石油的压缩系数,单位为10-3MPa-1;F为生烃转化率,取值范围为0-1;ΔP为油气充注之后的剩余源储压差;ΔP=P2-P1;其中,P1为生烃过程发生前烃源岩中的流体压力,P2为生烃过程发生后烃源岩中的流体压力,单位MPa;
构建的基于干酪根和生成油气比的生烃增压模型由如下式(1)表示:
(2)以步骤(1)构建的生烃增压模型为基础,在设定条件为发生油气充注的主要动力为生烃增压前提下,建立油气充注之后的剩余源储压差ΔP、烃源岩孔隙水初始体积Vw1、干酪根初始体积Vk1、生烃转化率F、生成油气密度ρo、储层孔隙度之间的关系式,从而计算得到所述致密油充注距离;
步骤(2)建立的充注距离计算公式由如下式(2)表示:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811155522.3A CN109188558B (zh) | 2018-09-29 | 2018-09-29 | 一种应用生烃增压模型计算致密油充注距离的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811155522.3A CN109188558B (zh) | 2018-09-29 | 2018-09-29 | 一种应用生烃增压模型计算致密油充注距离的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109188558A CN109188558A (zh) | 2019-01-11 |
CN109188558B true CN109188558B (zh) | 2020-08-11 |
Family
ID=64907087
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811155522.3A Active CN109188558B (zh) | 2018-09-29 | 2018-09-29 | 一种应用生烃增压模型计算致密油充注距离的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109188558B (zh) |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102681027A (zh) * | 2012-02-28 | 2012-09-19 | 北京师范大学 | 一种广布型复杂岩性油气藏成藏模式评价系统和方法 |
CN106285652A (zh) * | 2015-05-29 | 2017-01-04 | 中国石油化工股份有限公司 | 确定页岩游离气体饱和度的方法 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8967271B2 (en) * | 2012-06-07 | 2015-03-03 | Kellogg Brown & Root Llc | Subsea overpressure relief device |
CN106646659B (zh) * | 2017-01-19 | 2018-09-07 | 中国石油大学(华东) | 一种排烃门限以下地层超压的量化表征方法 |
CN109188557B (zh) * | 2018-09-29 | 2020-07-10 | 中国石油天然气股份有限公司 | 一种应用生烃增压模型计算排烃效率的方法 |
CN109633754B (zh) * | 2018-12-29 | 2021-05-28 | 中国地质调查局油气资源调查中心 | 一种天然气水合物开发模拟实验装置的模拟方法 |
-
2018
- 2018-09-29 CN CN201811155522.3A patent/CN109188558B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102681027A (zh) * | 2012-02-28 | 2012-09-19 | 北京师范大学 | 一种广布型复杂岩性油气藏成藏模式评价系统和方法 |
CN106285652A (zh) * | 2015-05-29 | 2017-01-04 | 中国石油化工股份有限公司 | 确定页岩游离气体饱和度的方法 |
Non-Patent Citations (1)
Title |
---|
中国陆相湖盆致密油成藏主控因素综述;马洪,等;《石油实验地质》;20141130;第36卷(第6期);668-677 * |
Also Published As
Publication number | Publication date |
---|---|
CN109188558A (zh) | 2019-01-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106056459B (zh) | 一种基于排烃效率的致密油源岩分级标准划分方法 | |
Guo et al. | Quantitative estimation of overpressure caused by oil generation in petroliferous basins | |
Laurent et al. | Pore-pressure influence in the poroelastic behavior of rocks: experimental studies and results | |
CN111638317B (zh) | 一种大型淡水湖盆陆相泥页岩油形成与演化评价方法 | |
CN104697914B (zh) | 一种泥页岩不同类型有机孔的预测方法 | |
CN103983746A (zh) | 一种富火山物质烃源岩有机碳恢复方法 | |
CN105181909A (zh) | 一种基于油气生排烃机理的有机碳恢复系数方法 | |
CN111694068A (zh) | 一种大型淡水湖盆陆相泥页岩油形成与富集评价方法 | |
Huang et al. | Quantitative modelling of hydrocarbon expulsion and quality grading of tight oil lacustrine source rocks: Case study of Qingshankou 1 member, central depression, Southern Songliao Basin, China | |
CN106646659B (zh) | 一种排烃门限以下地层超压的量化表征方法 | |
Gu et al. | Numerical analysis on gas production from silty hydrate reservoirs in the South China sea by depressurizing: The effect of permeability reduction caused by pore compression | |
Hinz et al. | A model for effective permeability in an unconsolidated hydrate reservoir | |
CN109188557B (zh) | 一种应用生烃增压模型计算排烃效率的方法 | |
CN103454408B (zh) | 一种测算致密砂岩油聚集量的方法与装置 | |
CN109188558B (zh) | 一种应用生烃增压模型计算致密油充注距离的方法 | |
CN108663722B (zh) | 一种油气资源评价中外部氢指数的确立方法及其应用 | |
Schneider et al. | Ceres2D: A numerical prototype for HC potential evaluation in complex area | |
CN116950654A (zh) | 一种低渗致密砂岩气藏开发效果评价方法 | |
Yan et al. | Prediction model of coal reservoir pressure and its implication for the law of coal reservoir depressurization | |
CN113944459B (zh) | 一种页岩原位转化可采油气资源量的预测方法及装置 | |
Adeyanju et al. | Experimental studies of sand production from unconsolidated sandstone petroleum reservoirs in Niger-Delta | |
Cook et al. | Simulation of a North Sea field experiencing significant compaction drive | |
CN107478768A (zh) | 一种泥页岩古抗剪强度确定方法 | |
Zeng et al. | Impact of swelling area expansion from the fracture wall into the matrix on the evolution of coal permeability | |
Huo et al. | Evolution and recovery of original total organic carbon for muddy source rocks with different total organic carbon in the Western Sichuan Depression, Sichuan Basin, China |
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 |