CN108984816B - 一种基于生油增压方程计算烃源岩超压的方法 - Google Patents
一种基于生油增压方程计算烃源岩超压的方法 Download PDFInfo
- Publication number
- CN108984816B CN108984816B CN201810374015.2A CN201810374015A CN108984816B CN 108984816 B CN108984816 B CN 108984816B CN 201810374015 A CN201810374015 A CN 201810374015A CN 108984816 B CN108984816 B CN 108984816B
- Authority
- CN
- China
- Prior art keywords
- oil
- volume
- kerogen
- overpressure
- crude oil
- 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
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Production Of Liquid Hydrocarbon Mixture For Refining Petroleum (AREA)
Abstract
一种基于生油增压方程计算烃源岩超压的方法,基于含油气盆地烃源岩生油,充分利用干酪根、原油、孔隙水三者质量、体积、压缩系数等参数的内在关系,构建基于生油增压计算模型,进而对超压进行计算,提高超压值计算精度的同时,将为提高油气地质评价提供技术支持;本发明在国内针对生油增压模型,提出了基于生油增压后的超压计算方法,能够有效地利用计算公式,提高了超压值的计算精度。
Description
技术领域
本发明属于油气勘探开发技术中的地质评价技术,特别涉及一种 基于生油增压方程计算烃源岩超压的方法。
背景技术
含油气盆地超压特征是目前研究难点,已经成为盆地分析中不可 缺少的组成部分,在油气成藏过程中起着越来越重要的作用。含油气 盆地超压与油气生成、运移和聚集具有密切关系。超压的形成与多种 因素有关,包括压实不均衡、孔隙流体热膨胀,黏土矿物脱水,构造 挤压和烃类生成等。
国内外通过烃源岩生油增压定量研究的目前还比较少,烃源岩生 油增压是因为高密度的干酪根转化成低密度的原油而使孔隙流体发 生膨胀的结果,因此干酪根生油作用被认为是可以使含油气盆地形成 大规模超压的主要成因机制,生油增压强度与有机质丰度、有机质类 型和成熟度以及封闭条件有关(osborne and swarbrick,1997)。Bredehoeft等(1994)认为由生油作用产生的孔隙流体超压可以在低渗 透性烃源岩中保持比较长的时期,其超压大小可以达到静岩压力。 Robert(1999)认为由生油作用产生的超压足以大到超过低渗透性烃源 岩毛细管压力,将油从烃源岩中驱动到储层。生油增压不像由压实不 均衡(欠压实)所形成的超压那样具有特定的判断依据,但已有的研 究证明生油增压沉积盆地存在一些与超压相关的特殊现象(Meissner, 1976)。郭小文于2011年提出有机质生油增压的计算模型,该模型依 据物质守恒和空间守恒两个原理,揭示了有机质在生油过程产生超压 的大小,但该模型存在2个主要不足,第一个不足是将地质参数带入 公式模型有时计算出现负值的情况,是因为该方程只考虑了石油的压 缩性,没有考虑到石油在实际地质情况下由于含有大量的溶解气,其 在地下的密度要小于地表,从而同质量的石油地下体积要大于地表体 积,因此郭小文公式计算的超压要小于实际产生的超压,第二个不足是公式中没有能够体现孔隙水残留情况的相关参数,烃源岩随压实的 进行,孔隙水不断排出,含水量逐渐减少,如果考虑到烃源岩生油时 同时还存在一部分孔隙水,那么实际产生的超压要大于郭小文公式计 算的超压。
发明内容
为了克服上述现有技术的不足,本发明的目的在于提供一种基于 生油增压方程计算烃源岩超压的方法,首次在模型中考虑了同质量石 油地下体积大于地表体积和孔隙水残留参数,充分利用孔隙水残留参 数与石油地下体积的内在关系,摒弃了前人计算模型存在的弊端,提 高生油增压计算精度的同时,将为油气地质开发评价精度提供技术支 持。
为了达到上述目的,本发明的技术方案为:
一种基于生油增压方程计算烃源岩超压的方法,包括以下步骤:
步骤一、生油增压的计算依据两个基本原理:
b.物质守恒:干酪根减少的质量(Mk)全部转化为油,因此生成 的油质量为:Mo=AF Mk-------------(2)
(1)、(2)式中:Vk1——原始干酪根体积,——孔隙度,Vw—— 原始孔隙水体积,Vo——原始油体积,Mo——生成的石油质量;Mk——转化为石油的干酪根质量,A——是为了换算单位定义的参数 A=IH/1000,IH——干酪根的氢指数,F——转化率。
步骤二、生成的液态油使孔隙流体膨胀将产生一定的超压,使孔 隙水和干酪根压缩更强烈,压缩后的孔隙水和干酪根的体积分别为:
Vw2=mVw1(1-CwΔP)-------------------------(3)
Vk2=(1-AF)(1-CkΔP)Vk1-------------(4)
生成的石油在地下占据的体积Vo2:为孔隙水减少的体积 (Vw1-Vw2)和干酪根减少的体积(Vk1-Vk2)之和:
Vo2=(Vw1-Vw2)+(Vk1-Vk2)-------------(5)
将公式(3)、(4)带入公式(5)得:
Vo2=(1-m)Vw1+Vw1mCwΔP+AFVk1+(1-AF)CkΔP Vk1;------(6)
(3)、(4)、(5)式中:Vw2——生油之后的孔隙水体积,m——孔隙 水残留系数,Vw1——未生油状态下孔隙水体积,Cw——地层水压缩 系数,ΔP——生油增加的压力,Vk2——生油之后的干酪根体积,Ck——干酪根压缩系数,Vk1——原始干酪根体积;
步骤三、生成石油的体积除了等于孔隙水减少的体积和干酪根减 少的体积之和,还应该等于生成油的质量与其密度比值,所以:
Vo2=Mo/ρo------(7)
将公式(2)带入公式(7)得:
Vo2=VoBo(1-ΔP Co)------(8)
结合公式(6)与(8)得:
VoBo(1-ΔP Co)=(1-m)Vw1+Vw1mCwΔP+AFVk1+(1-AF)CkΔP Vk1
整理得:
(7)、(8)、(9)、(10)式中:ρ-——原油密度,Co——石油压缩系 数,Bo——原油体积系数,m-——孔隙水残留系数,n——石油残留 系数。
步骤四、烃源岩超压确定:将实验获得的原油密度ρo,石油压缩 系数Co,原油体积系数Bo,孔隙水残留系数m和石油残留系数n输 入公式(10)得出△P,便可实现烃源岩超压的确定。
发明的效果
本发明基于生油增压公式计算的烃源岩超压确定方法,将干酪根 参数、原油地下体积与孔隙水残留参数有机结合在一起,所计算的烃 源岩超压与实验室模拟的超压基本重叠,其精度大大提高。
本发明抛弃了烃源岩生油物理模拟对超压模拟的高成本、长时间 等弊端,仅用到干酪根、原油、孔隙水一些参数,就可以精确获得超 压值,而且也不用依靠传统的物理模拟实验室,降低了成本。
本发明基于生油增压模型的一种计算方法,所计算的超压值与盆 地生油物理模拟的超压值相比,其精度不会降低只会升高、而成本大 大降低。
附图说明
图1为本发明中的生油增压计算方法流程图。
图2为本发明中的烃源岩生油增压模型。
具体实施方式
下面结合附图对本发明的技术方案做详细叙述。
参照图1,一种基于生油增压方程计算烃源岩超压的方法,包括 以下步骤:
步骤一、生油增压的计算依据两个基本原理,参照图2:
b.物质守恒:干酪根减少的质量(Mk)全部转化为油,因此生成 的油质量为:Mo=AF Mk-------------(2)
(1)、(2)式中:Vk1-原始干酪根体积,-孔隙度,Vw-原始孔隙水 体积,Vo-原始油体积,Mo-生成的石油质量;Mk-转化为石油的干酪 根质量,A是为了换算单位定义的参数A=IH/1000,IH-干酪根的氢指 数,F-转化率。
步骤二、生成的液态油使孔隙流体膨胀将产生一定的超压,使孔 隙水和干酪根压缩更强烈,压缩后的孔隙水和干酪根的体积分别为:
Vw2=mVw1(1-CwΔP)-------------------------(3)
Vk2=(1-AF)(1-CkΔP)Vk1-------------(4)
生成的石油在地下占据的体积Vo2:为孔隙水减少的体积 (Vw1-Vw2)和干酪根减少的体积(Vk1-Vk2)之和:
Vo2=(Vw1-Vw2)+(Vk1-Vk2)-------------(5)
将公式(3)、(4)带入公式(5)得:
Vo2=(1-m)Vw1+Vw1mCwΔP+AFVk1+(1-AF)CkΔP Vk1;------(6)
(3)、(4)、(5)式中:Vw2-生油之后的孔隙水体积,m-孔隙水残留 系数,Vw1-未生油状态下孔隙水体积,Cw-地层水压缩系数,ΔP-生油 增加的压力,Vk2-生油之后的干酪根体积,Ck-干酪根压缩系数,Vk1- 原始干酪根体积;
步骤三、生成石油的体积除了等于孔隙水减少的体积和干酪根减 少的体积之和,还应该等于生成油的质量与其密度比值,所以:
Vo2=Mo/ρo------(7)
将公式(2)带入公式(7)得:
Vo2=VoBo(1-ΔP Co)------(8)
结合公式(6)与(8)得:
VoBo(1-ΔP Co)=(1-m)Vw1+Vw1mCwΔP+AFVk1+(1-AF)CkΔP Vk1
整理得:
(7)、(8)、(9)、(10)式中:ρo-原油密度,Co-石油压缩系数,Bo- 原油体积系数,m-孔隙水残留系数,n-石油残留系数。
步骤四、烃源岩超压确定:将实验获得的原油密度ρo,石油压缩 系数Co,原油体积系数Bo,孔隙水残留系数m和石油残留系数n输 入公式(10)得出△P,便可实现烃源岩超压的确定。
为了验证建立的生油增压方程的准确性,采用生油增压方程计算 的超压值与中石化无锡石油地质研究所生油增压物理模拟结果进行 比较,选取的烃源岩样品为东濮凹陷某某井黑色泥岩,实验样品为刚 达到成熟门限,有机质类型为I型。采用生油增压方程计算时所用到 的烃源岩地球化学参数和模拟样品一致,岩石孔隙度取地表泥岩的孔 隙度62%,由于岩石孔隙度高,所以取密度为1700kg/m3,石油残留 系数n取1,因为生油增压物理模拟实验是在封闭条件下进行的,干 酪根的密度取1 200kg/m3,压缩系数为1.4×10- 3MPa-1,石油密度取 900kg/m3,压缩系数取2.2×10-3MPa-1,地层水的压缩系数为0.44×10-3MPa-1。由于模拟的样品孔隙流体是连通的,因此和单位质量的样品 产生的超压效果相同,在相同的孔隙度条件下,计算的生油增压量和 物理模拟结果具有可比性。计算得到烃源岩成熟度Ro分别在0.7%、0.75%、0.85%和1%时由于生油作用产生的累计增压量所示。采用建立的生油增压方程计算得到的生油增压量和物理模拟结果非常接 近,实测的超压与计算结果相差都在3MPa以内,相对误差都小于 5%,说明所建立的生油增压模型比较可靠,可以用于计算烃源岩生 油增压演化过程。
对比生油增压物理模拟方法与生油增压方程计算方法可知,新方 法使用后,能精确判断超压值大小。该方法在一定程度上降低了其他 参数的影响,比较快捷方便够满足含油气盆地超压特征评价的要求。 因此该法提高了超压值精度的同时,降低了时间、人力和财力等成本, 具有一定的推广应用价值。
本领域的技术人员应当理解,由于含油气盆地超压特征也受其他 地质因素的影响,为了保证该方法的有效可行性,必须保障干酪根、 原油、孔隙水一些参数(质量、密度、压缩系数等)的精确获得,基 于生油增压模型计算结果才具有较高的精度。
Claims (1)
1.一种基于生油增压方程计算烃源岩超压的方法,其特征在于,包括以下步骤:
步骤一、生油增压的计算依据两个基本原理:
b.物质守恒:干酪根减少的质量(Mk)全部转化为油,因此生成的油质量为:Mo=A FMk-------------(2)
(1)、(2)式中:Vk1——原始干酪根体积,——孔隙度,Vw——原始孔隙水体积,Vo——原始油体积,Mo——生成的石油质量;Mk——转化为石油的干酪根质量,A——是为了换算单位定义的参数A=IH/1000,IH——干酪根的氢指数,F——转化率;
步骤二、生成的液态油使孔隙流体膨胀将产生一定的超压,使孔隙水和干酪根压缩更强烈,压缩后的孔隙水和干酪根的体积分别为:
Vw2=mVw1(1-CwΔP)-------------------------(3)
Vk2=(1-AF)(1-CkΔP)Vk1-------------(4)
生成的石油在地下占据的体积Vo2:为孔隙水减少的体积(Vw1-Vw2)和干酪根减少的体积(Vk1-Vk2)之和:
Vo2=(Vw1-Vw2)+(Vk1-Vk2)-------------(5)
将公式(3)、(4)带入公式(5)得:
Vo2=(1-m)Vw1+Vw1mCwΔP+AFVk1+(1-AF)CkΔP Vk1;------(6)
(3)、(4)、(5)式中:Vw2——生油之后的孔隙水体积,m——孔隙水残留系数,Vw1——未生油状态下孔隙水体积,Cw——地层水压缩系数,ΔP——生油增加的压力,Vk2——生油之后的干酪根体积,Ck——干酪根压缩系数,Vk1——原始干酪根体积;
步骤三、生成石油的体积除了等于孔隙水减少的体积和干酪根减少的体积之和,还应该等于生成油的质量与其密度比值,所以:
Vo2=Mo/ρo------(7)
将公式(2)带入公式(7)得:
Vo2=VoBo(1-ΔP Co)------(8)
结合公式(6)与(8)得:
VoBo(1-ΔP Co)=(1-m)Vw1+Vw1mCwΔP+AFVk1+(1-AF)CkΔP Vk1整理得:
(7)、(8)、(9)、(10)式中:ρ-——原油密度,Co——石油压缩系数,Bo——原油体积系数,m-——孔隙水残留系数,n——石油残留系数;
步骤四、烃源岩超压确定:将实验获得的原油密度ρo,石油压缩系数Co,原油体积系数Bo,孔隙水残留系数m和石油残留系数n输入公式(10)得出△P,便可实现烃源岩超压的确定。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810374015.2A CN108984816B (zh) | 2018-04-24 | 2018-04-24 | 一种基于生油增压方程计算烃源岩超压的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810374015.2A CN108984816B (zh) | 2018-04-24 | 2018-04-24 | 一种基于生油增压方程计算烃源岩超压的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108984816A CN108984816A (zh) | 2018-12-11 |
CN108984816B true CN108984816B (zh) | 2023-03-31 |
Family
ID=64541895
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810374015.2A Active CN108984816B (zh) | 2018-04-24 | 2018-04-24 | 一种基于生油增压方程计算烃源岩超压的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108984816B (zh) |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104929624A (zh) * | 2015-04-22 | 2015-09-23 | 中国地质大学(武汉) | 一种超压驱动下原油二次运移速率的计算方法 |
CN106056459A (zh) * | 2016-05-31 | 2016-10-26 | 中国石油大学(华东) | 一种基于排烃效率的致密油源岩分级评价标准划分方法 |
CN106646659A (zh) * | 2017-01-19 | 2017-05-10 | 中国石油大学(华东) | 一种排烃门限以下地层超压的量化表征方法 |
CN106886058A (zh) * | 2017-04-17 | 2017-06-23 | 中国石油大学(北京) | 致密油藏有效储层物性下限值的确定方法和装置 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
FR2904654B1 (fr) * | 2006-08-04 | 2008-10-03 | Inst Francais Du Petrole | Methode pour quantifier la formation et la retention d'hydrocarbures dans une roche mere |
US8165817B2 (en) * | 2009-03-09 | 2012-04-24 | Schlumberger Technology Corporation | Method for integrating reservoir charge modeling and downhole fluid analysis |
-
2018
- 2018-04-24 CN CN201810374015.2A patent/CN108984816B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104929624A (zh) * | 2015-04-22 | 2015-09-23 | 中国地质大学(武汉) | 一种超压驱动下原油二次运移速率的计算方法 |
CN106056459A (zh) * | 2016-05-31 | 2016-10-26 | 中国石油大学(华东) | 一种基于排烃效率的致密油源岩分级评价标准划分方法 |
CN106646659A (zh) * | 2017-01-19 | 2017-05-10 | 中国石油大学(华东) | 一种排烃门限以下地层超压的量化表征方法 |
CN106886058A (zh) * | 2017-04-17 | 2017-06-23 | 中国石油大学(北京) | 致密油藏有效储层物性下限值的确定方法和装置 |
Non-Patent Citations (1)
Title |
---|
生油增压定量模型及影响因素;郭小文等;《石油学报》;20110715(第04期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN108984816A (zh) | 2018-12-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111441758B (zh) | 页岩油气甜点区的预测方法及装置 | |
CN104018829B (zh) | 一种利用煤层气井生产数据测量气水相渗曲线的方法 | |
CN110984970B (zh) | 一种利用地层测试确定启动压力梯度的方法 | |
CN105930932B (zh) | 基于含气指数的页岩气层标准化无阻流量的获取方法 | |
CN109991123B (zh) | 页岩油资源可动性的地球化学评价方法 | |
CN106056459A (zh) | 一种基于排烃效率的致密油源岩分级评价标准划分方法 | |
CN104007482B (zh) | 一种基于各向异性有效场的泥页岩岩石物理模型方法 | |
CN104948163A (zh) | 一种页岩气井产能测定方法 | |
CN106570262B (zh) | 一种储层构型结构的描述方法 | |
CN114372352B (zh) | 渗流-温度双场耦合数值模拟预测复杂断块油藏储气库调峰能力方法 | |
CN110095546B (zh) | 一种致密储层毫米级样品油源直接分析对比方法 | |
CN105138749A (zh) | 一种地质约束下高成熟烃源岩原始有机碳恢复方法 | |
CN112814669A (zh) | 一种页岩油藏全生命周期采收率预测方法和系统 | |
CN110108856A (zh) | 超压背景下烃源岩有机质成熟度的预测方法 | |
CN112487620B (zh) | 一种页岩油可动资源量的评价方法 | |
CN110952976A (zh) | 一种气藏开发模式下单井开采稳产潜力评价方法 | |
CN108984816B (zh) | 一种基于生油增压方程计算烃源岩超压的方法 | |
CN109188557B (zh) | 一种应用生烃增压模型计算排烃效率的方法 | |
CN112818501B (zh) | 基于动态监测数据校正碳酸盐岩油藏静态渗透率的方法 | |
CN110095584B (zh) | 一种储层油水饱和度校正方法 | |
CN106321076A (zh) | 一种注水井启动压力测试方法 | |
CN108072726B (zh) | 利用全油色谱研究油藏动用程度及预测剩余油分布的方法 | |
CN111577264A (zh) | 裂缝孔隙型油藏水平井产能预测方法及装置 | |
CN112746836B (zh) | 基于层间干扰的油井各层产量计算方法 | |
WO2022012368A1 (zh) | 一种页岩原位转化可采油气资源量的预测方法及装置 |
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 |