CN106053514A - 一种烃源岩有机质生烃转化率的评价方法 - Google Patents
一种烃源岩有机质生烃转化率的评价方法 Download PDFInfo
- Publication number
- CN106053514A CN106053514A CN201610388254.4A CN201610388254A CN106053514A CN 106053514 A CN106053514 A CN 106053514A CN 201610388254 A CN201610388254 A CN 201610388254A CN 106053514 A CN106053514 A CN 106053514A
- Authority
- CN
- China
- Prior art keywords
- hydrocarbon
- sample
- organic matter
- formula
- source rock
- 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
- 239000004215 Carbon black (E152) Substances 0.000 title claims abstract description 192
- 229930195733 hydrocarbon Natural products 0.000 title claims abstract description 189
- 150000002430 hydrocarbons Chemical class 0.000 title claims abstract description 189
- 239000011435 rock Substances 0.000 title claims abstract description 107
- 239000005416 organic matter Substances 0.000 title claims abstract description 73
- 238000006243 chemical reaction Methods 0.000 title claims abstract description 54
- 238000000034 method Methods 0.000 title claims abstract description 40
- 238000011156 evaluation Methods 0.000 claims abstract description 27
- 238000000197 pyrolysis Methods 0.000 claims abstract description 17
- 238000012937 correction Methods 0.000 claims abstract description 8
- 229910052739 hydrogen Inorganic materials 0.000 claims description 39
- 239000001257 hydrogen Substances 0.000 claims description 39
- UFHFLCQGNIYNRP-UHFFFAOYSA-N Hydrogen Chemical compound [H][H] UFHFLCQGNIYNRP-UHFFFAOYSA-N 0.000 claims description 38
- HEDRZPFGACZZDS-UHFFFAOYSA-N Chloroform Chemical compound ClC(Cl)Cl HEDRZPFGACZZDS-UHFFFAOYSA-N 0.000 claims description 30
- 238000005336 cracking Methods 0.000 claims description 25
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 claims description 23
- 229910052799 carbon Inorganic materials 0.000 claims description 23
- 239000007789 gas Substances 0.000 claims description 12
- 238000005070 sampling Methods 0.000 claims description 12
- 238000012360 testing method Methods 0.000 claims description 12
- 238000002474 experimental method Methods 0.000 claims description 10
- 241000209094 Oryza Species 0.000 claims description 9
- 235000007164 Oryza sativa Nutrition 0.000 claims description 9
- 235000009566 rice Nutrition 0.000 claims description 9
- 238000012216 screening Methods 0.000 claims description 6
- 239000000203 mixture Substances 0.000 claims description 4
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 claims description 3
- 230000008707 rearrangement Effects 0.000 claims description 3
- 238000009795 derivation Methods 0.000 claims description 2
- 239000002026 chloroform extract Substances 0.000 claims 2
- 230000009189 diving Effects 0.000 claims 1
- 238000010998 test method Methods 0.000 claims 1
- 230000000007 visual effect Effects 0.000 claims 1
- 238000004458 analytical method Methods 0.000 abstract description 7
- 239000000463 material Substances 0.000 abstract description 3
- 238000004364 calculation method Methods 0.000 abstract description 2
- 238000013178 mathematical model Methods 0.000 abstract description 2
- 239000003208 petroleum Substances 0.000 abstract 1
- 239000000523 sample Substances 0.000 description 70
- 239000004575 stone Substances 0.000 description 4
- 230000032683 aging Effects 0.000 description 3
- 230000035800 maturation Effects 0.000 description 3
- 239000007788 liquid Substances 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 235000008331 Pinus X rigitaeda Nutrition 0.000 description 1
- 235000011613 Pinus brutia Nutrition 0.000 description 1
- 241000018646 Pinus brutia Species 0.000 description 1
- 241000009298 Trigla lyra Species 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000005553 drilling Methods 0.000 description 1
- 150000002431 hydrogen Chemical class 0.000 description 1
- 238000011005 laboratory method Methods 0.000 description 1
- 239000006101 laboratory sample Substances 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000003908 quality control method Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 239000004079 vitrinite Substances 0.000 description 1
- 238000005303 weighing Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N25/00—Investigating or analyzing materials by the use of thermal means
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N1/00—Sampling; Preparing specimens for investigation
- G01N1/02—Devices for withdrawing samples
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N1/00—Sampling; Preparing specimens for investigation
- G01N1/28—Preparing specimens for investigation including physical details of (bio-)chemical methods covered elsewhere, e.g. G01N33/50, C12Q
Landscapes
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Production Of Liquid Hydrocarbon Mixture For Refining Petroleum (AREA)
- Analysing Materials By The Use Of Radiation (AREA)
Abstract
本发明涉及一种烃源岩有机质生烃转化率的评价方法,尤其是一种建立在物质平衡原理上的有机质在生烃转化率的评价方法,属于油气资源评价技术领域。本发明采用实测烃源岩热解数据和一个简便的数学模型进行计算,同时考虑了HI的校正,并将校正后的HI应用于有机质生烃转化率的分析评价,其结果的准确性以及可靠性得到显著的提升。本发明具有快速、实施成本低、可操作性强、所得结果精确的优势,便于大规模的工业化应用。
Description
技术领域
本发明涉及一种烃源岩有机质生烃转化率的评价方法,尤其是一种建立在物质平衡原理上的有机质在生烃转化率的评价方法,属于油气资源评价技术领域。
技术背景
有机质生烃转化率是一个盆地进行油气资源评价的基本参数,目前的获取方法有2种,一种是选取未熟烃源岩在实验室内开展热模拟实验(Lewan和Williams,1987;高岗和黄志龙,2000;卢双舫等,1994,2006),获得不同温度阶段的人工熟化样品,进而对人工熟化样品开展Rock‐Eval热解实验,得到产烃率与有机质热演化程度Tmax或镜质组反射率Ro的关系,再通过Tmax或Ro建立与深度的关系,转化为深度与生烃转化率的关系图版;第二种方法是采用生烃动力学理论模型对未熟‐低熟烃源岩样品热模拟实验产率进行拟合,得到烃源岩生烃动力学参数,最后结合埋藏史‐热史进行生烃转化率计算(王剑秋等,1984;金强等,1986;卢双舫,1996,1997,2001;王民等,2010,2011;李婷婷等,2011)。
第一种方法由于实验室内快速升温人工熟化实验样品获得的Tmax或Ro,与地质低温慢速升温条件下的成熟度数据不匹配使得这种方法得到的生烃转化率应用受到质疑。第二种方法中生烃动力学模型参数众多,多解性问题一致没有解决,使得这一方法的应用受到限制。同时,一个地区的埋藏史‐热史恢复本身也是世界级难题,存在极大不确定性。除此外,上述两种方法均没有采用实际自然演化样品的地球化学数据对生烃转化率结果进行验证、校正与对比,上述方法的应用争议很大(王飞宇等,2002;王民等,2011)。
众所周知,有机质生烃转化率是一个盆地进行油气资源评价的基本参数。而目前所常用的实验法和化学动力学模型法由于实验和地质条件的成烃转化率不匹配,且实验耗时费用高,另外无实测数据校正,使得生烃转化率结果获取困难,可信度较低。此外,国内外采用岩石热解地球化学数据评价烃源岩生烃潜力时均没有考虑泥岩孔隙中残留和有机质中吸附烃的影响,即尚未开展氢指数HI的校正并将其校正后的氢指数HI应用于机质生烃转化率的评价。
发明内容
本发明针对现有技术存在的不足提供一种便捷、节约且可信的烃源岩有机质生烃转化率评价评价的方法。
本发明一种烃源岩有机质生烃转化率的评价方法;包括下述步骤:
步骤A
步骤A1
在不同深度不同取样点取样,取样时对同一取样点采取的样本数大于等于2,从每一个取样点所取样本中取一份出来,通过热解实验测量其地球化学数据;所述地球化学数据包括所取样品的有机碳含量TOC、游离烃S1、裂解烃S2、裂解最大峰温Tmax;
步骤A2
剔除S2小于0.2mg烃/g岩石的样品,剔除TOC小于0.3wt.%的样品;剩余样品即为待测样品;
步骤A3
根据依据公式(1)计算每个待测样品的氢指数HI,
HI=S2/TOC×100 (1);
步骤A4
依据HI与Tmax判识有机质类型图版对待测样品的有机质类型进行划分,划分的型号为I型、II1型、II2型和III型;
步骤A5
数据抽稀,从地表往地下深处,以L米为刻度,每隔L米作为一个深度区间,选择每个深度区间内所取样品的HI最大值标记为HImax,并记录该样品的氢指数HImax以及该样品取样点距离地表的垂直距离h;
步骤A6
在步骤A5筛选抽稀所得样品的基础上,以深度h为纵坐标,不同类型有机质烃源岩样品HImax值为横坐标进行数据点的在直角坐标上投点,建立不同类型有机质烃源岩样品HImax与深度h的关系图;所述深度h为地表到取样点的垂直距离;
步骤A7
以步骤A6筛选得到样品的氢指数HImax以及该样品所处深度,建立起深度h和氢指数HImax的指数函数的拟合关系,
HImax=s+m×e(kh+t) (2)
其中s、m、k、t均为待求常数,通过已知的HImax和h求解;所述h为该样品所在深度;
步骤A8
令h=0米并将其带入根据步骤A7求解出s、m、k、t的公式(2)中求解,得到原始氢指数HIo;
步骤B
步骤B1
从步骤A1所取得每个取点所取样品中取出另一份出来,作为氯仿抽提备用
样;对氯仿抽提备用样进行氯仿抽提后,测量其有机碳含量TOC抽、游离
烃S1抽、裂解烃S2抽、裂解最大峰温Tmax抽;
步骤B2
将氯仿抽提备用样的所测得有机碳含量TOC抽、游离烃S1抽、裂解烃S2抽、裂解最大峰温Tmax抽代入公式(3)中求解该氯仿抽提备用样的HI';
HI'=S2抽/TOC抽×100 (3)
步骤C
将步骤A所得HIo、步骤B所得HI'代入公式(4)中计算得出每个取样点的烃源岩有机质生烃转化率TR;
公式(4)中,
TR为有机质成烃转化率,为0到1之间的小数;
HIo有机质原始氢指数,单位为mg烃/g岩石;其通过步骤A8得出;
HI'为校正后的现今实测的烃源岩有机质氢指数,单位为mg烃/g岩石;其通过公式(3)计算得出。
本发明一种烃源岩有机质生烃转化率的评价方法;步骤A1中,在不同深度不同取样点取样,取样时对同一取样点采取的样本数为2,从每一个取样点所取样本中取一份出来,通过热解实验测量其地球化学数据;所述地球化学数据包括所取样品的有机碳含量TOC、游离烃S1、裂解烃S2、裂解最大峰温Tmax。
本发明一种烃源岩有机质生烃转化率的评价方法;普遍认为S2小于0.2mg烃/g或TOC小于0.3wt.%的烃源岩为非烃源岩。
本发明一种烃源岩有机质生烃转化率的评价方法;步骤A4中,依据HI与Tmax判识有机质类型图版对待测样品的有机质类型进行划分时;所述有机质类型图版为Espitalie J的Rock‐Eval pyrolysis and its applications(Rock‐Eval热解仪及其应用)论文中755‐784页所记载的有机质类型图版。该论文刊登于Rev I Fpet(1985年40卷(第6期))。
本发明一种烃源岩有机质生烃转化率的评价方法;步骤A5中所述L选自5-10之间的一个定值。
本发明一种烃源岩有机质生烃转化率的评价方法;步骤A7中s、m、k、t均为待求常数,通过已知的HI和h采用最小二乘法拟合法进行求解;所述h为该样品所在深度;所述深度h为地表到取样点的垂直距离。
本发明一种烃源岩有机质生烃转化率的评价方法;步骤A8中由于有机质在浅层或地表时尚未成熟,属于未熟烃源岩,因此轮廓线与X轴的交点可视为原始生烃潜力,即当取h=0m时对应的拟合氢指数HI值为原始氢指数HIo。
本发明一种烃源岩有机质生烃转化率的评价方法;步骤B1中,从待测样品所在取样点的取出样品中取出另一份出来,作为氯仿抽提备用样;对氯仿抽提备用样进行氯仿抽提后,测量其有机碳含量TOC抽、游离烃S1抽、裂解烃S2抽、裂解最大峰温Tmax抽。
本发明一种烃源岩有机质生烃转化率的评价方法;步骤B2中,求得一系列样品的HI'后,根据一系列同一取样点样品的氢指数HI与HI',建立拟和方程(5)并求解拟和方程(5)
HI′=a×eb×HI (5)
公式(5)中a、b为常数,e为自然对数底数;
本发明建立拟和方程(5),大大节省了后续样品检测的步骤,通过计算样品的HI,即可通过公式(2)求出HIo、通过公式(5)求出HI';进而根据公式(4)直接求得所测样品的有机质生烃转化率TR。
本发明一种烃源岩有机质生烃转化率的评价方法;公式(4)的推导过程如下:
根据有机质生烃转化率的定义:有机质生烃转化率等于有机质已生成的烃量与原始生烃潜量的比值,可得出公式(6):
公式(6)中,HIo×CTOC o×Go为岩石的原始生烃潜量,HI'×CTOC×Gx为岩石的残余生烃潜量;
公式(6)中,CTOC o为原始有机碳在样品中的质量百分含量;Go烃源岩原始质量,单位为g;CTOC为现今岩石实测有机碳在样品中的质量百分含量;Gx为现今烃源岩的质量,即烃源岩原始质量Go经历生烃和排烃后的残留质量,单位为g;整理(6)式得:
随着烃源岩的逐渐埋深,有机质逐渐成熟,岩石中的干酪根开始向油气转化;通过各油田大量的数据统计发现,油气产物中碳的含量基本相当,平均含量为83.33%,即1200mg的油气(烃)中含有1g碳;于是公式(6)式可表达为:
整理(8)式,得到:
化简公式(9)得到公式(10):
(10)式两边乘以得到:
将(11)式带入(7)式,可得:
重新整理(12)式,可得式(4),即:
本发明的有益效果:
本发明具有如下明显的优势:①原理正确、计算简便、结果可靠。本发明所建立的模型依据物质平衡原理(烃源岩生烃前后总的生烃潜量保持不变),采用实测烃源岩热解数据和一个简单的数学模型进行计算,同时考虑了HI的校正。并且本方法采用了自然演化样品的实测数据,理论扎实,实际可靠。②计算结果连续、省时、节约。烃源岩Rock‐Eval热解数据(TOC、S1、S2、Tmax)是一个油田开展勘探的基本地球化学分析实验,在长期的勘探开发过程中各油田积累了大量的烃源岩热解实验数据,可以实现不同深度、不同区块的生烃转化率连续计算。同时烃源岩Rock‐Eval热解分析实验费用廉价(目前完成一块样品分析不足300元人民币),所用的分析时间也很短(一块样品仅耗时不到30分钟)。
附图说明
图1烃源岩有机质成熟演化模式图;
图2松辽盆地青山口组一段泥岩Tmax‐HI有机质类型判别图版;
图3松辽盆地青山口组一段筛选后的泥岩氢指数与深度关系图;
图4岩石样品抽提前后氢指数对比图。
具体实施方式
下面以松辽盆地青山口组一段烃源岩有机质成烃转化率的计算为实施例,对本发明作进一步的详细说明。
步骤一、烃源岩有机质生烃转化率计算公式的建立
随着烃源岩埋深的逐渐增加(图1),烃源岩有机质从未熟到成熟,未熟烃源岩有机质总生烃潜量为HIo×CTOC o×Go,成熟后烃源岩有机质残余生烃潜量为HI'×CTOC×Gx,其中转化掉的生烃潜量为HIo×CTOC o×Go‐HI'×CTOC×Gx,故此根据有机质成烃转化率的定义为有机质已生成的烃量与原始生烃潜量的比值,可知生烃转化率为:
式(1)中,TR为有机质成烃转化率,为0到1之间的小数;HIo有机质原始氢指数,mg烃/g岩石;HI'为校正后的现今实测的烃源岩有机质氢指数,mg烃/g岩石;CTOC o为原始有机碳含量,wt.%,即质量百分数;Go烃源岩原始质量,g;CTOC,为现今岩石实测有机碳含量,wt.%;Gx为烃源岩的现今质量,即Go的烃源岩生烃和排烃后质量,g。整理(1)式得:
随着烃源岩的逐渐埋深,有机质逐渐成熟,岩石中的干酪根开始向油气转化。一般认为油气产物中碳的含量为83.33%,即1200mg的油气(烃)中含有1g碳。于是(1)式可表达为:
整理(8)式,得到:
(10)式两边乘以得到:
将(11)式带入(7)式,可得:
重新整理(12)式,可得:
步骤二、有机质原始氢指数(HIo)的获取:
①收集整理油田钻井取芯岩石热解地球化学数据,包括样品的井位及取样深度,所取样品有机碳含量‐TOC、游离烃‐S1、裂解烃S2、裂解最大峰温Tmax;取样时要求样品从浅到深均有,即每间隔5米选一个取样点、每个取样点取2个样;
②地球化学数据的质量控制,剔除S2小于0.2mg烃/g岩石的样品,剔除TOC小于0.3wt.%的样品;
③依据氢指数HI计算公式(HI=S2/TOC×100)计算烃源岩样品氢指数;
④依据HI与Tmax判识有机质类型图版对挑选的烃源岩样品有机质类型进行划分,划分为I型、II1型、II2型和III型,详见图2(图版引自Espitalie et al,1985);
⑤数据抽稀,从地表往地下深处,以5m深度为刻度,每隔5m深度,选择该深度区间的HI最大值,并记录该氢指数HI为HImax以及其对应的深度h。
①在⑤基础山筛选抽稀后得到的样品,以深度为纵坐标,不同类型有机质烃源岩样品HImax值为横坐标进行数据点的在直角坐标上投点,建立不同类型有机质烃源岩样品HImax与深度关系图(图2);
⑦以⑥筛选得到的氢指数HImax及其对应的深度,建立起深度h和氢指数HImax的指数函数的拟合关系,即HImax=s+m×e(kh+t)。其中s、m、k、t均为未知待求参数(拟合方法采用最小二乘法拟合),求得
HImax=120+172×e(0.0054h+1.63) (I型)
HImax=65+155×e(0.0048h+1.24) (II1型)
HImax=20+134×e(0.0043h+0.92) (II2型)
HImax=20+72×e(0.0038h+0.58) (III型)
⑧求取原始氢指数HIo,即当取h=0米时对应的拟合氢指数HImax值,从图3中可以看出,I型有机质HIo为850mg烃/gTOC,II1有机质HIo为600mg烃/gTOC,II2有机质HIo为350mg烃/gTOC,III有机质HIo为150mg烃/gTOC。
步骤三、有机质残余氢指数(HI)的校正,即HI'的获取
HI'为现今岩石残留氢指数HI的校正结果。HI的定义为:
式中,S2为岩石残留裂解烃量,mg/g岩石,为300℃以后岩石热解挥发得到的烃量。已有的研究表明,岩石热解分析得到S2(裂解烃)中存在部分先前生成的液态烃,这部分液态烃分子量较大,Rock‐Eval分析中300℃之前尚未蒸发出(Jarvie,2012),在孔隙中残留和有机质中吸附。因此,岩石热解中得到的S2不全为残余裂解烃,因此,需要对其进行校正。
为此,设计开展岩石样品直接热解和进行氯仿抽提后再进行热解两种实验,分别获得抽提后岩石样品的TOC抽、S2抽和未抽提样品的TOC和S2,依据公式(9)可计算得到抽提前样品的氢指数(HI)和抽提后样品的氢指数HI’,作图4,通过指数函数拟合可以获得HI’的校正公式(10)。
HI'=59.192×e0.0037×HI (5)
最后,根据生烃转化率计算公式(8)及原始氢指数HIo和HI’校正公式(5)对松辽盆地青山口组一段泥岩生烃转化率进行计算,即烃源岩有机质生烃转化率
经过计算,各样品点的生烃转化率TR见表1。
表1 松辽盆地青山口组一段烃源岩有机质成烃转化率的计算成果表
Claims (9)
1.一种烃源岩有机质生烃转化率的评价方法;其特征在于包括下述步骤:
步骤A
步骤A1
在不同深度不同取样点取样,取样时对同一取样点采取的样本数大于等于2,从每一个取样点所取样本中取一份出来,通过热解实验测量其地球化学数据;所述地球化学数据包括所取样品的有机碳含量TOC、游离烃S1、裂解烃S2、裂解最大峰温Tmax;
步骤A2
剔除S2小于0.2mg烃/g岩石的样品,剔除TOC小于0.3wt.%的样品;剩余样品即为待测样品;
步骤A3
根据依据公式(1)计算每个待测样品的氢指数HI,
HI=S2/TOC×100 (1);
步骤A4
依据HI与Tmax判识有机质类型图版对待测样品的有机质类型进行划分,划分的型号为I型、II1型、II2型和III型;
步骤A5
数据抽稀,从地表往地下深处,以L米为刻度,每隔L米作为一个深度区间,选择每个深度区间内所取样品的HI最大值,并标记为HImax,并记录该样品的氢指数HImax以及该样品取样点距离地表的垂直距离h;
步骤A6
在步骤A5筛选抽稀所得样品的基础上,以深度h为纵坐标,不同类型有机质烃源岩样品HImax值为横坐标进行数据点的在直角坐标上投点,建立不同类型有机质烃源岩样品HImax与深度h的关系图;所述深度h为地表到取样点的垂直距离;
步骤A7
以步骤A6筛选得到样品的氢指数HImax以及该样品所处深度h,建立起深度h和氢指数HImax的指数函数的拟合关系,
HImax=s+m×e(kh+t) (2)
其中s、m、k、t均为待求常数,通过已知的HImax和h求解;所述h为该样品所在深度;
步骤A8
令h=0米并将其带入根据步骤A7求解出s、m、k、t的公式(2)中求解,得到原始氢指数HIo;
步骤B
步骤B1
从步骤A1所取得每个取点所取样品中取出另一份出来,作为氯仿抽提备用样;对氯仿抽提备用样进行氯仿抽提后,测量其有机碳含量TOC抽、游离烃S1抽、裂解烃S2抽、裂解最大峰温Tmax抽;
步骤B2
将氯仿抽提备用样的所测得有机碳含量TOC抽、游离烃S1抽、裂解烃S2抽、裂解最大峰温Tmax抽代入公式(3)中求解该氯仿抽提备用样的HI';
HI'=S2抽/TOC抽×100 (3)
步骤C
将步骤A所得HIo、步骤B所得HI'代入公式(4)中计算得出每个取样点的烃源岩有机质生烃转化率TR;
公式(4)中,
TR为有机质成烃转化率,为0到1之间的小数;
HIo有机质原始氢指数,单位为mg烃/g岩石;其通过步骤A8得出;
HI'为校正后的现今实测的烃源岩有机质氢指数,单位为mg烃/g岩石;其通过公式(3)计算得出。
2.根据权利要求1所述的一种烃源岩有机质生烃转化率的评价方法;其特征在于:步骤A1中,在不同深度不同取样点取样,取样时对同一取样点采取的样本数为2,从每一个取样点所取样本中取一份出来,通过热解实验测量其地球化学数据;所述地球化学数据包括所取样品的有机碳含量TOC、游离烃S1、裂解烃S2、裂解最大峰温Tmax。
3.根据权利要求1所述的一种烃源岩有机质生烃转化率的评价方法;其特征在于:步骤A4中,依据HI与Tmax判识有机质类型图版对待测样品的有机质类型进行划分时;所述有机质类型图版为Espitalie J等的Rock‐Eval pyrolysis and its applications论文中755‐784页所记载的有机质类型图版。
4.根据权利要求1所述的一种烃源岩有机质生烃转化率的评价方法;其特征在于:步骤A5中所述L选自5-10米之间的一个定值。
5.根据权利要求1所述的一种烃源岩有机质生烃转化率的评价方法;其特征在于:步骤A7中s、m、k、t均为待求常数,通过已知的HImax和h采用最小二 乘法拟合法进行求解;所述h为该样品所在深度;所述深度h为地表到取样点的垂直距离。
6.根据权利要求1所述的一种烃源岩有机质生烃转化率的测试方法;其特征在于:步骤A8中由于有机质在浅层或地表时尚未成熟,属于未熟烃源岩,因此轮廓线与X轴的交点可视为原始生烃潜力,即当取h=0米时对应的拟合氢指数HImax值为原始氢指数HIo。
7.根据权利要求1所述的一种烃源岩有机质生烃转化率的评价方法;其特征在于:步骤B1中,从待测样品所在取样点的取出样品中取出另一份出来,作为氯仿抽提备用样;对氯仿抽提备用样进行氯仿抽提后,测量其有机碳含量TOC抽、游离烃S1抽、裂解烃S2抽、裂解最大峰温Tmax抽。
8.根据权利要求1所述的一种烃源岩有机质生烃转化率的评价方法;其特征在于:步骤B2中,求得一系列样品的HI'后,根据一系列同一取样点样品的氢指数HI与HI',建立拟和方程(5)并求解拟和方程(5)
HI′=a×eb×HI (5)
公式(5)中a、b为常数,e为自然对数底数。
9.根据权利要求1所述的一种烃源岩有机质生烃转化率的评价方法;其特征在于:公式(4)的推导过程如下:
根据有机质生烃转化率的定义:有机质生烃转化率等于有机质已生成的烃量与原始生烃潜量的比值,可得出公式(6):
公式(6)中,HIo×CTOC o×Go为岩石的原始生烃潜量,HI'×CTOC×Gx为岩石的残余生烃潜量;
公式(6)中,CTOC o为原始有机碳在样品中的质量百分含量;Go烃源岩原始质量,单位为g;CTOC为现今岩石实测有机碳在样品中的质量百分含量;Gx为现今烃源岩的质量,即烃源岩原始质量Go经历生烃和排烃后的残留质量,单位为g;整理(6)式得:
随着烃源岩的逐渐埋深,有机质逐渐成熟,岩石中的干酪根开始向油气转化;通过各油田大量的数据统计发现,油气产物中碳的含量基本相当,平均含量为83.33%,即1200mg的油气烃中含有1g碳;于是公式(6)式可表达为:
整理(8)式,得到:
化简公式(9)得到公式(10):
(10)式两边乘以得到:
将(11)式带入(7)式,可得:
重新整理(12)式,可得式(4),即:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610388254.4A CN106053514B (zh) | 2016-06-02 | 2016-06-02 | 一种烃源岩有机质生烃转化率的评价方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610388254.4A CN106053514B (zh) | 2016-06-02 | 2016-06-02 | 一种烃源岩有机质生烃转化率的评价方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106053514A true CN106053514A (zh) | 2016-10-26 |
CN106053514B CN106053514B (zh) | 2019-04-09 |
Family
ID=57169442
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610388254.4A Expired - Fee Related CN106053514B (zh) | 2016-06-02 | 2016-06-02 | 一种烃源岩有机质生烃转化率的评价方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106053514B (zh) |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106908854A (zh) * | 2017-02-08 | 2017-06-30 | 中国石油天然气股份有限公司 | 目标区域的甜点区的确定方法和装置 |
CN109612794A (zh) * | 2018-12-06 | 2019-04-12 | 同济大学 | 一种烃源岩中不同赋存态有机质的分离与定量方法 |
CN110513106A (zh) * | 2019-08-08 | 2019-11-29 | 中国石油大学(北京) | 补偿轻烃蒸发量的烃源岩生排烃特征的表征方法及系统 |
CN110779831A (zh) * | 2019-11-19 | 2020-02-11 | 中国石油大学(华东) | 一种页岩储层有机质热裂解生烃转化率的计算方法及热解产物收集装置 |
CN112083034A (zh) * | 2020-09-14 | 2020-12-15 | 中国石油大学(华东) | 原位电加热过程中烃源岩生烃量确定方法及系统 |
CN112461780A (zh) * | 2020-10-20 | 2021-03-09 | 长江大学 | 一种烃源岩原始总有机碳含量的评价与计算方法 |
CN112632738A (zh) * | 2019-10-09 | 2021-04-09 | 中国石油天然气股份有限公司 | 泥页岩样品原位转化中生烃活化能的确定方法及装置 |
CN116741296A (zh) * | 2023-06-12 | 2023-09-12 | 中国石油大学(华东) | 一种生烃演化定量表征方法、系统、介质、设备及终端 |
CN116741296B (zh) * | 2023-06-12 | 2024-06-21 | 中国石油大学(华东) | 一种生烃演化定量表征方法、系统、介质、设备及终端 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103543470A (zh) * | 2013-10-18 | 2014-01-29 | 中国石油大学(华东) | 一种热解岩石中的游离烃/残留烃s1的轻、重烃校正方法 |
WO2014036077A1 (en) * | 2012-08-28 | 2014-03-06 | Saudi Arabian Oil Company | Method for reconstructing the total organic carbon content from compositional modeling analysis |
CN105572320A (zh) * | 2015-12-15 | 2016-05-11 | 中国石油天然气股份有限公司 | 一种烃源岩生、排、滞留烃量的定量评价方法 |
-
2016
- 2016-06-02 CN CN201610388254.4A patent/CN106053514B/zh not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2014036077A1 (en) * | 2012-08-28 | 2014-03-06 | Saudi Arabian Oil Company | Method for reconstructing the total organic carbon content from compositional modeling analysis |
CN103543470A (zh) * | 2013-10-18 | 2014-01-29 | 中国石油大学(华东) | 一种热解岩石中的游离烃/残留烃s1的轻、重烃校正方法 |
CN105572320A (zh) * | 2015-12-15 | 2016-05-11 | 中国石油天然气股份有限公司 | 一种烃源岩生、排、滞留烃量的定量评价方法 |
Non-Patent Citations (3)
Title |
---|
J.ESPITALIE ET AL.: "Rock-Eval Pyrolysis and Its Applications (Part two)", 《REV. INST. FR. PÉT.》 * |
LU SHUANGFANG ET AL.: "Chemical kinetics study of hydrocarbon regeneration from organic matter in carbonate source rocks and its significance", 《SCIENCE IN CHINA SERIES D: EARTH SCIENCES》 * |
王民: "有机质生烃动力学及火山作用的热效应研究与应用", 《中国博士学位论文全文数据库 基础科学辑》 * |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106908854A (zh) * | 2017-02-08 | 2017-06-30 | 中国石油天然气股份有限公司 | 目标区域的甜点区的确定方法和装置 |
CN106908854B (zh) * | 2017-02-08 | 2019-01-18 | 中国石油天然气股份有限公司 | 目标区域的甜点区的确定方法和装置 |
CN109612794A (zh) * | 2018-12-06 | 2019-04-12 | 同济大学 | 一种烃源岩中不同赋存态有机质的分离与定量方法 |
CN109612794B (zh) * | 2018-12-06 | 2021-05-11 | 同济大学 | 一种烃源岩中不同赋存态有机质的分离与定量方法 |
CN110513106A (zh) * | 2019-08-08 | 2019-11-29 | 中国石油大学(北京) | 补偿轻烃蒸发量的烃源岩生排烃特征的表征方法及系统 |
CN112632738A (zh) * | 2019-10-09 | 2021-04-09 | 中国石油天然气股份有限公司 | 泥页岩样品原位转化中生烃活化能的确定方法及装置 |
CN110779831A (zh) * | 2019-11-19 | 2020-02-11 | 中国石油大学(华东) | 一种页岩储层有机质热裂解生烃转化率的计算方法及热解产物收集装置 |
CN112083034A (zh) * | 2020-09-14 | 2020-12-15 | 中国石油大学(华东) | 原位电加热过程中烃源岩生烃量确定方法及系统 |
CN112461780A (zh) * | 2020-10-20 | 2021-03-09 | 长江大学 | 一种烃源岩原始总有机碳含量的评价与计算方法 |
CN112461780B (zh) * | 2020-10-20 | 2024-01-05 | 长江大学 | 一种烃源岩原始总有机碳含量的评价与计算方法 |
CN116741296A (zh) * | 2023-06-12 | 2023-09-12 | 中国石油大学(华东) | 一种生烃演化定量表征方法、系统、介质、设备及终端 |
CN116741296B (zh) * | 2023-06-12 | 2024-06-21 | 中国石油大学(华东) | 一种生烃演化定量表征方法、系统、介质、设备及终端 |
Also Published As
Publication number | Publication date |
---|---|
CN106053514B (zh) | 2019-04-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106053514B (zh) | 一种烃源岩有机质生烃转化率的评价方法 | |
CN104278991B (zh) | 盐湖相烃源岩有机碳和生烃潜量的多元测井计算方法 | |
Yu et al. | Porosity estimation in kerogen-bearing shale gas reservoirs | |
CN105221133A (zh) | 一种基于测井多参数确定烃源岩有机碳含量的方法和装置 | |
CN104298883A (zh) | 一种油气资源评价中的烃源岩产烃率图版的建立方法 | |
CN104656162A (zh) | 一种富铀烃源岩铀含量及有机碳含量确定方法 | |
CN103543470A (zh) | 一种热解岩石中的游离烃/残留烃s1的轻、重烃校正方法 | |
CN104897712A (zh) | 一种页岩含油量的测定方法与系统 | |
CN107966545B (zh) | 一种任意成熟度烃源岩的有机碳恢复方法及其装置 | |
CN112307601B (zh) | 一种复杂储层可压裂性评价方法 | |
CN104453874A (zh) | 一种基于核磁共振的砂砾岩储层含油饱和度的计算方法 | |
CN107780923B (zh) | 一种基于泥质校正的含水饱和度模型的建立、仿真方法 | |
CN108508182B (zh) | 快速确定笔石相热页岩中生物硅含量的测录井方法 | |
CN108561126B (zh) | 一种确定页岩气储层有机孔隙度的简易方法 | |
CN105138749A (zh) | 一种地质约束下高成熟烃源岩原始有机碳恢复方法 | |
CN105403585A (zh) | 烃源岩内分散液态烃丰度下限的确定方法 | |
CN107687986A (zh) | 一种评价沉积岩中总有机碳含量的方法 | |
CN103678778A (zh) | 一种放射性物化探信息综合方法 | |
CN105257284B (zh) | 一种利用元素俘获能谱测井确定凝灰质含量的方法及装置 | |
CN107807222A (zh) | 一种确定生排烃实验热成熟和生烃进程的方法和装置 | |
CN105205296A (zh) | 一种求取页岩气储层孔隙度的方法 | |
CN110320571A (zh) | 一种致密砂岩储层岩石脆性测井评价方法 | |
CN108333637A (zh) | 一种提高元素测井技术确定元素含量准确度的方法 | |
CN102914797B (zh) | 一种获得地层各向异性系数的方法及装置 | |
CN105003258A (zh) | 一种高温高压气层甲烷流体密度骨架参数的获取方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20190409 |