CN108279295B - 一种综合性确定碳酸盐岩层系烃源岩供烃底限的方法 - Google Patents

一种综合性确定碳酸盐岩层系烃源岩供烃底限的方法 Download PDF

Info

Publication number
CN108279295B
CN108279295B CN201711475358.XA CN201711475358A CN108279295B CN 108279295 B CN108279295 B CN 108279295B CN 201711475358 A CN201711475358 A CN 201711475358A CN 108279295 B CN108279295 B CN 108279295B
Authority
CN
China
Prior art keywords
hydrocarbon
source rock
vitrinite
lowest limit
reflectance
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
Application number
CN201711475358.XA
Other languages
English (en)
Other versions
CN108279295A (zh
Inventor
庞雄奇
郑天昱
汪文洋
王恩泽
郑定业
赵毅
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China University of Petroleum Beijing
Original Assignee
China University of Petroleum Beijing
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by China University of Petroleum Beijing filed Critical China University of Petroleum Beijing
Priority to CN201711475358.XA priority Critical patent/CN108279295B/zh
Publication of CN108279295A publication Critical patent/CN108279295A/zh
Application granted granted Critical
Publication of CN108279295B publication Critical patent/CN108279295B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N33/00—Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
    • G01N33/24—Earth materials

Landscapes

  • Life Sciences & Earth Sciences (AREA)
  • Health & Medical Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Chemical & Material Sciences (AREA)
  • Food Science & Technology (AREA)
  • Analytical Chemistry (AREA)
  • Geology (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Environmental & Geological Engineering (AREA)
  • Medicinal Chemistry (AREA)
  • Physics & Mathematics (AREA)
  • Remote Sensing (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Investigating Or Analysing Materials By Optical Means (AREA)
  • Investigating Or Analyzing Non-Biological Materials By The Use Of Chemical Means (AREA)

Abstract

本发明提供了一种综合性确定碳酸盐岩层系烃源岩供烃底限的方法，其包括以下步骤：根据氢碳原子比H/C与Ro的变化关系，拟合H/C值演化曲线，以H/C＝0.01对应的Ro为供烃底限Ro1；根据氧碳原子比O/C与Ro的变化关系，拟合O/C值演化曲线，以O/C＝0.02对应的Ro作为供烃底限Ro2；根据生烃潜力指数“(S1+S2)/TOC”与Ro的变化关系，拟合排烃曲线，并计算出生烃量，以生烃量的最后1％所对应的Ro作为供烃底限Ro3；根据Ro1、Ro2和Ro3中的最大值和最小值，确定供烃底限范围；然后根据不同地区Ro和深度H的变化关系，对应供烃底限范围，得到不同地区的供烃底限深度范围。

Description

一种综合性确定碳酸盐岩层系烃源岩供烃底限的方法
技术领域
本发明属于储层评价技术领域，涉及一种综合性确定碳酸盐岩层系烃源岩供烃底限的方法。
背景技术
碳酸盐岩层系烃源岩在全球广泛分布但是成熟程度普遍高，这些高-过成熟烃源岩在晚期还能够生成多少油气，这关乎到我国碳酸盐岩油气潜力和勘探前景。
塔里木盆地位于中国新疆境内，是中国最大的含油气沉积盆地，面积达56×104m2。塔里木盆地台盆区系指古生代海相地层的地区，主要包括塔北隆起、北部坳陷区(包括阿瓦提凹陷和满加尔凹陷)、中央隆起区(包括巴楚隆起、塔中隆起、塔东低凸起)等构造单元，如图1所示。截止2010年，台盆区累计探明26个油气田，发现了36个含油气构造，探明石油地质储量为12.06×108t，约占整个盆地探明储量的95％，而天然气探明储量3940.1×1012，约占整个盆地探明储量的35％，是海相油气富集的主要区域。
塔里木盆地台盆区自下而上发育古生界、中生界和新生界。研究表明，下古生界寒武系和奥陶系是台盆区烃源岩主要发育层位，烃源岩包括泥质和碳酸盐岩两种，碳酸盐岩烃源岩可分为寒武系-下奥陶统和中上奥陶统两套，其中寒武系-下奥陶统以中下寒武统烃源岩为主，而上寒武统-下奥陶统烃源岩较差。中下寒武统酸盐岩层系烃源岩主要分布在塔里木盆地东部欠补偿盆地相和西部蒸发泻湖相，有机质类型为Ⅰ型，TOC为0.5％-5.52％，Ro为1.64％-3.61％，处于高-过成熟的生气阶段。中-上奥陶统碳酸盐岩层系烃源岩主要分布在塔里木盆地中部和满加尔凹陷西部的台缘斜坡灰泥丘相，有机质类型以Ⅰ型为主并混有Ⅱ2-Ⅲ型，TOC为0.5％-5.4％，Ro为0.81％-1.30％，处于成熟-高成熟生油阶段。
因此，探究塔里木盆地台盆区烃源岩还能够生成多少油气，确定该碳酸盐岩层系烃源岩供烃底限的探讨便具有重要理论和现实意义。目前已有一些科学家对供烃底限的影响因素、获取做过一些研究，并取得了一些重要成果，但也存在一些问题。有的研究方法单一而局限[1-3]，主要是通过生烃热模拟实验或者生烃动力学方法；有的对供烃底限的确定只是定性的分析，而基本没有定量的计算[4-7]；要么多是针对泥页岩和煤，而对碳酸盐岩烃源岩供烃底限的研究很少。而且现有关于碳酸盐岩层系烃源岩供烃底限的确定方法并不准确[1,8-9]。
[1]Tissot B P,Welte D H,Petroleum Formation and occurrence[M].Berlin,Heidelberg,New York:Spring-Vevlag,1978；
[2]Cramer B,Faber E,Gerling P,et al.Reaction kinetics of stablecarbon isotope in natural gas insights from dry,open system pyrolysisexperiments[J].Energy Fuels,2001,15:517-532；
[3]张有生，秦勇，刘焕杰，等.沉积有机质二次生烃热模拟实验研究[J].地球化学,2002,31(3):273-282；
[4]Su Long,Zheng Jianjing,Chen Guojun,et al.The upper limit ofmaturity of natural gas generation and its implication for the YachengFormation in the Qiongdongnan Basin,China[J].Journal of Asian Earth Sciences,2012,(54/55):203-213；
[5]秦勇，张有生，朱炎铭，等.煤中有机质二次生烃迟滞性及其反应动力学机制[J].地球科学：中国地质大学学报，2000,25(3):278-282；
[6]谢增业，蒋助生，张英，等.全岩热模拟新方法及其在气源岩评价中的应用[J].沉积学报，2002,20(3):510-514；
[7]王云鹏，赵长毅，王兆利，等.用生烃动力学方法确定海相有机质的主生气期及其初步应用[J].石油勘探与开发，2005,32(4):153-158；
[8]陈建平，赵文智，王招明，等.海相干酪根天然气生成成熟度上限与生气潜力极限探讨——以塔里木盆地研究为例[J].科学通报，2007,52(增刊1)：95-100；
[9]胡国艺，魏志平，肖志华，等.腐殖煤干酪根裂解气主生气成熟度上限探讨[J].沉积学报，2006,24(4):585-589。
发明内容
鉴于上述现有技术的缺点，本发明的目的在于提供一种综合性确定碳酸盐岩层系烃源岩供烃底限的方法，该方法能从烃源岩的生烃、残留烃、排烃的机理出发，对碳酸盐岩烃源岩的供烃底限进行综合系统的研究和定量的确定。
为了达到前述的发明目的，本发明提供一种综合性确定碳酸盐岩层系烃源岩供烃底限的方法，其包括以下步骤：
步骤一：在待测区域采集不同深度目标层段的烃源岩样品；
步骤二：测试每个烃源岩样品所含有机质的镜质体反射率Ro；
步骤三：测试每个烃源岩样品中干酪根的氢碳原子比H/C和氧碳原子比O/C；
步骤四：对每个烃源岩进行热解定量评价，获得所述烃源岩的可溶烃含量S1、裂解烃含量S2，以及有机质丰度TOC含量；
步骤五：根据各烃源岩样品测得的镜质体反射率Ro及干酪根的氢碳原子比H/C，拟合形成氢碳原子比H/C值随镜质体反射率Ro的演化曲线，并以氢碳原子比H/C＝0.1所对应的镜质体反射率Ro为供烃底限Ro1；
步骤六：根据各烃源岩样品测得的镜质体反射率Ro及干酪根的氧碳原子比O/C，拟合形成氧碳原子比O/C值随镜质体反射率Ro的演化曲线，并以氧碳原子比O/C＝0.02所对应的镜质体反射率Ro为供烃底限Ro2；
步骤七：根据各烃源岩样品测得的镜质体反射率Ro及生烃潜力指数P，拟合形成残余生烃潜力曲线，根据所述残余生烃潜力曲线计算生烃量，并根据生烃总量Q的99％所对应的镜质体反射率Ro为供烃底限Ro3，其中，所述生烃潜力指数P由以下公式(1)计算获得：
P＝(S1+S2)/TOC (1)；
步骤八：取所述Ro1、Ro2和Ro3中的最大值和最小值，作为该待测区域烃源岩的供烃底限范围。
上述综合性确定碳酸盐岩层系烃源岩供烃底限的方法中，S1为热解试验中岩石中残留烃的含量(即可溶烃含量，单位为mg(烃)/g(岩石))；S2为岩石中干酪根在热解过程中新生成的烃类的含量(即裂解烃含量，单位为mg(烃)/g(岩石))；TOC为有机质丰度，以百分比％表示；Ro为镜质体反射率，以百分比％，Q为生烃总量，单位为mg(烃)/g(岩石))。
根据本发明的具体实施例，优选地，所述方法还包括确定待测区域供烃底限深度范围的步骤，具体包括：
根据各烃源岩样品测得的镜质体反射率Ro及其所在的深度，拟合形成镜质体反射率Ro与深度的函数关系式；
将所述Ro1、Ro2和Ro3中的最大值和最小值分别代入所述函数关系式中计算出供烃底限深度范围的两个端值，得到所述供烃底限深度范围。
根据本发明的具体实施例，优选地，采用SPSS专业统计软件进行拟合，形成镜质体反射率Ro与深度的函数关系式。
根据本发明的具体实施例，优选地，采用SPSS专业统计软件进行拟合，形成氢碳原子比H/C值随镜质体反射率Ro的演化曲线，氧碳原子比O/C值随镜质体反射率Ro的演化曲线以及残余生烃潜力曲线。
根据本发明的具体实施例，优选地，所述残余生烃潜力曲线随Ro先增大后减小，且所述残余生烃潜力曲线转折点所对应的Ro为排烃作用产生的门限值Roe。
根据本发明的具体实施例，优选地，所述生烃总量Q通过以下公式(2)计算：
所述公式(2)中，f(Ro)表示所述残余生烃潜力曲线中生烃潜力指数P关于镜质体反射率Ro的函数表达式，所述Ro max表示镜质体反射率Ro的最大值，所述Romin表示镜质体反射率Ro的最小值。所述Ro max，Ro min可以为本领域常规认定的Ro的最大、最小值。
根据本发明的具体实施例，优选地，所述残余生烃潜力曲线包括两部分曲线，第一部分曲线生烃潜力指数P随Ro的增大而增大，第二部分曲线烃潜力指数P随Ro的增大而减小，所述残余生烃潜力曲线的转折点处的烃潜力指数P最大，该转折点对应的Ro值为排烃作用产生的门限值Roe。
根据本发明的具体实施例，优选地，所述第一部分曲线对应的生烃量Q1通过以下公式(3)计算，所述第二部分曲线对应的生烃量Q2通过以下公式(4)，
所述公式(3)中，f1(Ro)表示所述第一部分残余生烃潜力曲线中生烃潜力指数P关于镜质体反射率Ro的函数表达式，f2(Ro)表示所述第二部分残余生烃潜力曲线中生烃潜力指数P关于镜质体反射率Ro的函数表达式；
所述生烃总量Q为所述生烃量Q1与生烃量Q2的总和。
根据本发明的具体实施例，优选地，所述生烃总量的99％所对应的镜质体反射率Ro为供烃底限Ro3通过以下公式(5)或(6)计算：
根据本发明的具体实施例，优选地，所述步骤一中，在同一深度目标层段的烃源岩样品个数为2-10个。
与现有技术相比，本发明的有益效果在于：本发明提供的综合性确定碳酸盐岩层系烃源岩供烃底限的方法从烃源岩的生烃、残留烃、排烃的机理出发，利用多种方法对碳酸盐岩烃源岩的供烃底限进行系统的研究和定量的确定，获得在实际地质条件下的供烃底限。
附图说明
图1是塔里木盆地地理位置分布图；
图2A是实施例1采集的不同深度塔里木盆地台盆区寒武系-奥陶系烃源岩样品中干酪根的氢碳原子比H/C随镜质体反射率Ro变化的分布图；
图2B是利用SPSS专业统计软件对图2A中的氢碳原子比H/C与镜质体反射率Ro变化关系初步拟合形成的曲线图；
图2C是利用SPSS专业统计软件对图2B中的氢碳原子比H/C与镜质体反射率Ro变化关系拟合形成的函数曲线图；
图2D是对图2C拟合形成的函数曲线的分析图；
图3A是实施例1采集的不同深度塔里木盆地台盆区寒武系-奥陶系烃源岩样品中干酪根的氧碳原子比O/C随镜质体反射率Ro变化的分布图；
图3B是利用SPSS专业统计软件对图3A中的氧碳原子比O/C与镜质体反射率Ro变化关系初步拟合形成的曲线图；
图3C是利用SPSS专业统计软件对图2B中的氧碳原子比O/C与镜质体反射率Ro变化关系拟合形成的函数曲线图；
图4A是实施例1采集的不同深度塔里木盆地台盆区寒武系-奥陶系烃源岩样品的生烃潜力指数(S1+S2)/TOC随镜质体反射率Ro变化的分布图；
图4B是利用SPSS专业统计软件对图4A中的生烃潜力指数(S1+S2)/TOC与镜质体反射率Ro变化关系初步拟合形成的曲线图；
图4C是利用SPSS专业统计软件对图4B中的生烃潜力指数(S1+S2)/TOC与镜质体反射率Ro变化关系在门限值Roe前拟合形成的函数曲线图；
图4D是利用SPSS专业统计软件对图4B中的生烃潜力指数(S1+S2)/TOC与镜质体反射率Ro变化关系在门限值Roe后拟合形成的函数曲线图；
图5A是实施例1确定门限值Roe前生烃量的原理图；
图5B是实施例1确定门限值Roe后生烃量的原理图；
图6A塔里木盆地台盆区中塔中地区单井镜质体反射率Ro与深度关系图；
图6B塔里木盆地台盆区中英东2地区单井镜质体反射率Ro与深度关系图；
图6C塔里木盆地台盆区中巴楚地区单井镜质体反射率Ro与深度关系图；
图6D塔里木盆地台盆区中轮南46地区单井镜质体反射率Ro与深度关系图。
具体实施方式
为了对本发明的技术特征、目的和有益效果有更加清楚的理解，现对本发明的技术方案进行以下详细说明，但不能理解为对本发明可实施范围的限定。
实施例1
本实施例提供了一种综合性确定塔里木盆地台盆区寒武系-奥陶系供烃底限的方法，其包括以下步骤：
步骤一：采集塔里木盆地台盆区寒武系-奥陶系在不同深度目标层段中的烃源岩样品，每一深度采集5个样品；
步骤二：测试每个烃源岩样品所含有机质的镜质体反射率Ro；
步骤三：测试每个烃源岩样品中干酪根的氢碳原子比H/C和氧碳原子比O/C；
步骤四：对每个烃源岩进行热解定量评价，获得所述烃源岩的可溶烃含量S1、裂解烃含量S2，以及有机质丰度TOC含量；
步骤五：对各烃源岩样品进行氢碳原子比H/C有机元素法分析，具体步骤为：
统计各烃源岩样品测得的镜质体反射率Ro及干酪根的氢碳原子比H/C，制成图表，如图2A所示；
将图2A中每个Ro值对应的氢碳原子比H/C值取平均值，即以作为分析点，然后在的分析点集合中以Ro为坐标每间隔0.1取点，形成新的分析点集群；利用SPSS专业统计软件对该分析点集群的氢碳原子比H/C与镜质体反射率Ro变化关系初步拟合，形成氢碳原子比H/C值随镜质体反射率Ro的演化曲线，如图2B所示；
进一步利用SPSS专业统计软件对图2B中的氢碳原子比H/C与镜质体反射率Ro变化关系拟合形成的函数曲线图，如图2C所示，在图2C中，氢碳原子比H/C(设为x)与镜质体反射率Ro(设为y)的函数关系满足关系式y＝-1.428In(x)+1.0908；相关性R2＝0.9966；
将图2C拟合形成的函数曲线进行翻转，即将图2B中的氢碳原子比H/C(设为y)与镜质体反射率Ro(设为x)，然后利用SPSS专业统计软件对图2B中的氢碳原子比H/C与镜质体反射率Ro变化关系重新拟合形成新的函数曲线图，如图2D所示，在图2D中，氢碳原子比H/C(y)与镜质体反射率Ro(x)的函数关系满足关系式y＝2.1364e-0.698x，相关性R2＝0.9966，与图2C拟合获得的函数关系相比，图2D拟合获得的函数关系更贴近实际分布；也更便于分析，以实际烃源岩样品测得的Ro为参照范围，结合常规的Ro范围，暂取Ro的最大值Ro max为6％，最小值Ro min为0.5％，由图2D拟合的函数曲线可知，当H/C为0.1时，镜质体反射率Ro的最小值为4.35％，因此，采用氢碳原子比H/C有机元素法分析获得的供烃底限Ro1为4.35％；
步骤六：对各烃源岩样品进行氧碳原子比O/C有机元素法分析，具体步骤为：
统计各烃源岩样品测得的镜质体反射率Ro及干酪根的氧碳原子比O/C，制成图表，如图3A所示；
将图3A中每个Ro值对应的氧碳原子比O/C值取平均值，即以作为分析点，然后在的分析点集合中以Ro为坐标每间隔0.1取点，形成新的分析点集群；利用SPSS专业统计软件对该分析点集群的氧碳原子比O/C与镜质体反射率Ro变化关系初步拟合，形成氧碳原子比O/C值随镜质体反射率Ro的演化曲线，如图3B所示；
进一步利用SPSS专业统计软件对图3B中的氧碳原子比O/C与镜质体反射率Ro变化关系拟合形成的函数曲线图，如图3C所示，在图3C中，氧碳原子比O/C(设为x)与镜质体反射率Ro(设为y)的函数关系满足关系式y＝-2.038In(x)-3.775；相关性R2＝0.9982；
以实际烃源岩样品测得的Ro为参照范围，结合常规的Ro范围，暂取Ro的最大值Romax为6％，最小值Ro min为0.5％；由图3C拟合的函数曲线可知，当O/C为0.02时，镜质体反射率Ro的最小值为4.26％，因此，采用氧碳原子比O/C有机元素法分析获得的供烃底限Ro2为4.26％。
步骤七：对各烃源岩样品进行生烃潜力法分析，生烃潜力法的原理是基于在源岩热解定量评价中，通常用可溶烃(S1)与裂解烃(S2)的和表示源岩的生烃潜力。此处采用一个综合热解参数，即生烃潜力指数((S1+S2)/TOC)来表征源岩的生烃潜力。当源岩的生烃潜力指数在地史演化过程中开始减小时，则表明有烃类开始排出，其相应的埋深则代表了源岩的排烃门限。该生烃潜力法分析的具体步骤为：
统计各烃源岩样品测得的镜质体反射率Ro及生烃潜力指数(S1+S2)/TOC，制成图表，如图4A所示；
将图4A中每个Ro值对应的生烃潜力指数(S1+S2)/TOC取平均值，即以作为分析点，然后在的分析点集合中以Ro为坐标每间隔0.1取点，形成新的分析点集群；利用SPSS专业统计软件对该分析点集群的生烃潜力指数(S1+S2)/TOC与镜质体反射率Ro变化关系初步拟合，形成生烃潜力指数(S1+S2)/TOC值随镜质体反射率Ro的演化曲线，即残余生烃潜力曲线，如图4B所示；该残余生烃潜力曲线包括两部分曲线，第一部分曲线生烃潜力指数(S1+S2)/TOC随Ro的增大而增大，第二部分曲线烃潜力指数(S1+S2)/TOC随Ro的增大而减小，因此，该残余生烃潜力曲线的转折点处的烃潜力指数(S1+S2)/TOC最大，该转折点对应的Ro值为排烃作用产生的门限值Roe，图4B中残余生烃潜力曲线的门限值Roe＝0.9；
利用SPSS专业统计软件对图4B中门限值Roe＝0.9之前的残余生烃潜力曲线进行拟合，形成的函数曲线图，如图4C所示，在图4C中，生烃潜力指数(S1+S2)/TOC(设为x)与镜质体反射率Ro(设为y)的函数关系f1满足关系式y＝0.3957In(x)-1.4742；相关性R2＝0.9966；
利用SPSS专业统计软件对图4B中门限值Roe＝0.9之后的残余生烃潜力曲线进行拟合，形成的函数曲线图，如图4D所示，在图4D中，生烃潜力指数(S1+S2)/TOC(设为x)与镜质体反射率Ro(设为y)的函数关系f2满足关系式y＝-0.957In(x)+6.9951；相关性R2＝0.9911；
以实际烃源岩样品测得的Ro为参照范围，结合常规的Ro范围，暂取Ro的最大值Romax为6％，最小值Ro min为0.5％；
将图4C拟合形成的函数曲线进行翻转，即将图4B中的生烃潜力指数(S1+S2)/TOC(设为y)与镜质体反射率Ro(设为x)，然后利用SPSS专业统计软件对图4B中的生烃潜力指数(S1+S2)/TOC与镜质体反射率Ro变化关系重新拟合形成新的函数曲线图，如图5A所示，在图5A中，生烃潜力指数(S1+S2)/TOC与镜质体反射率Ro(x)的函数关系F1满足关系式y＝41.96e2.5184x，相关性R2＝0.9966，与图4C拟合获得的函数关系相比，图5A拟合获得的函数关系更贴近实际分布；也更便于分析，以实际烃源岩样品测得的Ro为参照范围，结合常规的Ro范围，暂取Ro的最小值Ro min为0.5％，由如4B可知，门限值Roe＝0.9；对该F1在Ro∈[0.5，0.9]范围积分，计算门限值前的生烃量Q1，即
将图4D拟合形成的函数曲线进行翻转，即将图4B中的生烃潜力指数(S1+S2)/TOC(设为y)与镜质体反射率Ro(设为x)，然后利用SPSS专业统计软件对图4B中的生烃潜力指数(S1+S2)/TOC与镜质体反射率Ro变化关系重新拟合形成新的函数曲线图，如图5B所示，在图5B中，生烃潜力指数(S1+S2)/TOC与镜质体反射率Ro(x)的函数关系F2满足关系式y＝1459.8e-1.035x，相关性R2＝0.9911，与图4D拟合获得的函数关系相比，图5B拟合获得的函数关系更贴近实际分布；也更便于分析，以实际烃源岩样品测得的Ro为参照范围，结合常规的Ro范围，暂取Ro的最大值Ro min为6％，由如4B可知，门限值Roe＝0.9；对该F2在Ro∈[0.9，6.0]范围积分，计算门限值后的生烃量Q2，即
令计算出Ro3的值为4.83％；
因此，采用生烃潜力法分析获得的供烃底限Ro3为4.83％。
步骤八：由上述步骤五至步骤七计算获得Ro1、Ro2和Ro3统计结果如表1所示，表1为不同方法确定的供烃底限Ro对比表。由表1比较可知，该塔里木盆地台盆区寒武系-奥陶系烃源岩的供烃底限范围为4.26％-4.83％，比常规指标更贴近实际地质条件下的生油气特征。
表1
步骤九：对各烃源岩样品进行供烃底限深度分析，具体包括：
分别根据区域各烃源岩样品测得的镜质体反射率Ro及其所在的深度，统计各烃源岩样品测得的镜质体反射率Ro及其所在的深度H，制成图表，如图6A-6D所示；
将图6A中每个Ro值对应的深度H取平均值，即h，以(Ro，h)作为分析点，然后在(Ro，h)的分析点集合中以Ro为坐标每间隔0.1取点，形成新的分析点集群；利用SPSS专业统计软件对该分析点集群的深度h与镜质体反射率Ro变化关系初步拟合，形成深度h值随镜质体反射率Ro的演化曲线，利用SPSS专业统计软件对图6A中的曲线进行拟合，形成的函数曲线图，如图6A所示，在图6A中，深度h(设为x)与镜质体反射率Ro(设为y)的函数关系满足关系式：y＝2936.4In(x)+4770.1；相关性R2＝0.8349；
将上述计算得到的Ro1、Ro2和Ro3中的最大值和最小值分别代入上述函数关系式中计算出供烃底限深度范围的两个端值，得到供烃底限深度范围。
采用相同的方法对图6B-6D中的镜质体反射率Ro及深度H关系进行拟合，形成的函数曲线图，如图6B所示，在图6B中，深度h(设为x)与镜质体反射率Ro(设为y)的函数关系满足关系式y＝1872.6In(x)+2804.7；相关性R2＝0.7179；
如图6C所示，在图6C中，深度h(设为x)与镜质体反射率Ro(设为y)的函数关系满足关系式y＝5030.8In(x)+1706.4；相关性R2＝0.9034；
如图6D所示，在图6D中，深度h(设为x)与镜质体反射率Ro(设为y)的函数关系满足关系式y＝801.64In(x)+6014.8；相关性R2＝0.9805；
根据不同地区Ro与深度H的关系式，将上述计算得到的Ro1、Ro2和Ro3中的最大值和最小值分别代入上述各函数关系式中计算出不同地区供烃底限深度范围的两个端值，得到不同地区供烃底限深度范围，然后根据图6A-6D得到的函数关系式判断获得塔里木盆地不同地区碳酸盐岩层系烃源岩供烃底限Ro对应的深度H，如表2所示，表2为塔里木盆地台盆区不同地区碳酸盐岩烃源岩供烃底限深度表。
表2
由表2可知，塔里木盆地台盆区寒武系-奥陶系碳酸盐岩层系烃源岩主要分布在四个地区：塔东地区(包括满加尔凹陷)、塔中隆起、塔北轮南和巴楚隆起。而塔里木盆地不同地区古地温条件存在差异，这导致了不同地区烃源岩热演化程度的不同，表现在不同地区深度H与Ro关系不同(如图6A-6D所示)。
本实施例通过对塔里木盆地碳酸盐岩层系烃源岩供烃底限的确定，解决了目前对碳酸盐岩层系烃源岩基本没有定量的计算的问题，为该地区甚至我国碳酸盐岩油气潜力和勘探前景提供了一个评价的依据。
由上述实施例可知，本发明把供烃底限定义为沉积盆地油气生成的死亡线对应的有机质成熟度Ro或深度H，并对塔里木盆地台盆区海相碳酸盐岩层系烃源岩供烃底限Ro进行的研究，操作更简单可行；有些和某些学者研究的结果较为相近，也有些结果偏大，但都更符合实际地质条件下的生油气特征，因此，本发明提供的求取的供烃底限的方法是相对准确和可信的。
值得注意的是，本文应用的方法或多或少需要利用数据统计的结果，而统计数据的多少、差异和可靠性可能对求取结果有影响，因此需要对获取的数据进行分析筛选，选择层位相同、岩性和沉积相类似的数据，并删除为储层样品点的数据。
以上所述的具体实施方式，对本发明的目的、技术方案和有益效果进行了进一步详细说明，所应理解的是，以上所述仅为本发明的具体实施方式而已，并不用于限制本发明，凡在本发明的精神和原则之内，所做的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。

Claims (10)

1.一种综合性确定碳酸盐岩层系烃源岩供烃底限的方法，包括以下步骤：
步骤一：在待测区域采集不同深度目标层段的烃源岩样品；
步骤二：测试每个烃源岩样品所含有机质的镜质体反射率Ro；
步骤三：测试每个烃源岩样品中干酪根的氢碳原子比H/C和氧碳原子比O/C；
步骤四：对每个烃源岩进行热解定量评价，获得所述烃源岩的可溶烃含量S1、裂解烃含量S2，以及有机质丰度TOC含量；
步骤五：根据各烃源岩样品测得的镜质体反射率Ro及干酪根的氢碳原子比H/C，拟合形成氢碳原子比H/C值随镜质体反射率Ro的演化曲线，并以氢碳原子比H/C＝0.1所对应的镜质体反射率Ro为供烃底限Ro1；
步骤六：根据各烃源岩样品测得的镜质体反射率Ro及干酪根的氧碳原子比O/C，拟合形成氧碳原子比O/C值随镜质体反射率Ro的演化曲线，并以氧碳原子比O/C＝0.02所对应的镜质体反射率Ro为供烃底限Ro2；
步骤七：根据各烃源岩样品测得的镜质体反射率Ro及生烃潜力指数P，拟合形成残余生烃潜力曲线，根据所述残余生烃潜力曲线计算生烃量，并根据生烃总量Q的99％所对应的镜质体反射率Ro为供烃底限Ro3，其中，所述生烃潜力指数P由以下公式(1)计算获得：
P＝(S1+S2)/TOC (1)；
步骤八：取所述Ro1、Ro2和Ro3中的最大值和最小值，作为该待测区域烃源岩的供烃底限范围。
2.根据权利要求1所述的综合性确定碳酸盐岩层系烃源岩供烃底限的方法，其特征在于：所述方法还包括确定待测区域供烃底限深度范围的步骤，具体包括：
根据各烃源岩样品测得的镜质体反射率Ro及其所在的深度，拟合形成镜质体反射率Ro与深度的函数关系式；
将所述Ro1、Ro2和Ro3中的最大值和最小值分别代入所述函数关系式中计算出供烃底限深度范围的两个端值，得到所述供烃底限深度范围。
3.根据权利要求2所述的综合性确定碳酸盐岩层系烃源岩供烃底限的方法，其特征在于：采用SPSS专业统计软件进行拟合，形成镜质体反射率Ro与深度的函数关系式。
4.根据权利要求1所述的综合性确定碳酸盐岩层系烃源岩供烃底限的方法，其特征在于：采用SPSS专业统计软件进行拟合，形成氢碳原子比H/C值随镜质体反射率Ro的演化曲线，氧碳原子比O/C值随镜质体反射率Ro的演化曲线以及残余生烃潜力曲线。
5.根据权利要求1所述的综合性确定碳酸盐岩层系烃源岩供烃底限的方法，其特征在于：所述残余生烃潜力曲线随Ro先增大后减小，且所述残余生烃潜力曲线转折点所对应的Ro为排烃作用产生的门限值Roe。
6.根据权利要求1所述的综合性确定碳酸盐岩层系烃源岩供烃底限的方法，其特征在于：所述生烃总量Q通过以下公式(2)计算：
所述公式(2)中，f(Ro)表示所述残余生烃潜力曲线中生烃潜力指数P关于镜质体反射率Ro的函数表达式，所述Ro max表示镜质体反射率Ro的最大值，所述Ro min表示镜质体反射率Ro的最小值。
7.根据权利要求6所述的综合性确定碳酸盐岩层系烃源岩供烃底限的方法，其特征在于：所述残余生烃潜力曲线包括两部分曲线，第一部分曲线生烃潜力指数P随Ro的增大而增大，第二部分曲线烃潜力指数P随Ro的增大而减小，所述残余生烃潜力曲线的转折点处的烃潜力指数P最大，该转折点对应的Ro值为排烃作用产生的门限值Roe。
8.根据权利要求7所述的综合性确定碳酸盐岩层系烃源岩供烃底限的方法，其特征在于：所述第一部分曲线对应的生烃量Q1通过以下公式(3)计算，所述第二部分曲线对应的生烃量Q2通过以下公式(4)，
所述公式(3)中，f1(Ro)表示所述第一部分残余生烃潜力曲线中生烃潜力指数P关于镜质体反射率Ro的函数表达式，f2(Ro)表示所述第二部分残余生烃潜力曲线中生烃潜力指数P关于镜质体反射率Ro的函数表达式；
所述生烃总量Q为所述生烃量Q1与生烃量Q2的总和。
9.根据权利要求8所述的综合性确定碳酸盐岩层系烃源岩供烃底限的方法，其特征在于：所述生烃总量的99％所对应的镜质体反射率Ro为供烃底限Ro3通过以下公式(5)或(6)计算：
10.根据权利要求1所述的综合性确定碳酸盐岩层系烃源岩供烃底限的方法，其特征在于：所述步骤一中，在同一深度目标层段的烃源岩样品个数为2-10个。
CN201711475358.XA 2017-12-29 2017-12-29 一种综合性确定碳酸盐岩层系烃源岩供烃底限的方法 Active CN108279295B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201711475358.XA CN108279295B (zh) 2017-12-29 2017-12-29 一种综合性确定碳酸盐岩层系烃源岩供烃底限的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201711475358.XA CN108279295B (zh) 2017-12-29 2017-12-29 一种综合性确定碳酸盐岩层系烃源岩供烃底限的方法

Publications (2)

Publication Number Publication Date
CN108279295A CN108279295A (zh) 2018-07-13
CN108279295B true CN108279295B (zh) 2019-10-11

Family

ID=62802613

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201711475358.XA Active CN108279295B (zh) 2017-12-29 2017-12-29 一种综合性确定碳酸盐岩层系烃源岩供烃底限的方法

Country Status (1)

Country Link
CN (1) CN108279295B (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111077586B (zh) * 2018-10-18 2023-08-01 中国石油化工股份有限公司 一种复杂断陷盆地油气主成藏期的确定方法及装置
CN110513106A (zh) * 2019-08-08 2019-11-29 中国石油大学(北京) 补偿轻烃蒸发量的烃源岩生排烃特征的表征方法及系统
CN113407899B (zh) * 2021-08-20 2021-10-26 中国科学院地质与地球物理研究所 海相高过成熟烃源岩排烃量评价方法

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101929960A (zh) * 2009-06-22 2010-12-29 中国石油天然气股份有限公司 定量测定烃源岩中干酪根结构组成和成熟度的方法
CN103760307A (zh) * 2014-01-16 2014-04-30 中国石油大学(华东) 一种评价干酪根类型的方法
CN103926388A (zh) * 2013-01-10 2014-07-16 中国石油大学(北京) 一种判别低丰度有效烃源岩的方法
CN104298883A (zh) * 2014-10-16 2015-01-21 中国石油大学(华东) 一种油气资源评价中的烃源岩产烃率图版的建立方法
CN104700316A (zh) * 2015-04-02 2015-06-10 中国石油大学(华东) 一种非均质性评价泥页岩生油气量的方法
CN106803021A (zh) * 2017-01-23 2017-06-06 中国石油大学(北京) 一种常规与非常规储层的油气资源量的评价方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101929960A (zh) * 2009-06-22 2010-12-29 中国石油天然气股份有限公司 定量测定烃源岩中干酪根结构组成和成熟度的方法
CN103926388A (zh) * 2013-01-10 2014-07-16 中国石油大学(北京) 一种判别低丰度有效烃源岩的方法
CN103760307A (zh) * 2014-01-16 2014-04-30 中国石油大学(华东) 一种评价干酪根类型的方法
CN104298883A (zh) * 2014-10-16 2015-01-21 中国石油大学(华东) 一种油气资源评价中的烃源岩产烃率图版的建立方法
CN104700316A (zh) * 2015-04-02 2015-06-10 中国石油大学(华东) 一种非均质性评价泥页岩生油气量的方法
CN106803021A (zh) * 2017-01-23 2017-06-06 中国石油大学(北京) 一种常规与非常规储层的油气资源量的评价方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
四川盆地上三叠系须家河组油气资源评价;郑定业等;《特种油气藏》;20170831;第24卷(第4期);第67-72页 *

Also Published As

Publication number Publication date
CN108279295A (zh) 2018-07-13

Similar Documents

Publication Publication Date Title
Hu et al. Oil content evaluation of lacustrine organic-rich shale with strong heterogeneity: A case study of the Middle Permian Lucaogou Formation in Jimusaer Sag, Junggar Basin, NW China
Xu et al. Characteristics of source rocks and genetic origins of natural gas in deep formations, Gudian Depression, Songliao Basin, NE China
Tang et al. Upper Paleozoic coal measures and unconventional natural gas systems of the Ordos Basin, China
Peters et al. Recognition of an Infracambrian source rock based on biomarkers in the Baghewala-1 oil, India
Pang et al. The dead line for oil and gas and implication for fossil resource prediction
Katz et al. A new look at the components of the petroleum system of the South Caspian Basin
CN108279295B (zh) 一种综合性确定碳酸盐岩层系烃源岩供烃底限的方法
Xulong et al. Carboniferous source rocks in the Ludong-Wucaiwan area, Junggar basin, NW China
Zheng et al. Hydrocarbon generation and expulsion characteristics of the source rocks in the third member of the Upper Triassic Xujiahe Formation and its effect on conventional and unconventional hydrocarbon resource potential in the Sichuan Basin
Li et al. Organic geochemistry, distribution and hydrocarbon potential of source rocks in the Paleocene, Lishui Sag, East China Sea Shelf Basin
Hakimi et al. Petroleum source rock characterisation and hydrocarbon generation modeling of the Cretaceous sediments in the Jiza sub-basin, eastern Yemen
Hall et al. Cooper Basin source rock geochemistry
Sylwan Source rock properties of Vaca Muerta formation, Neuquina basin, Argentina
Gong et al. Geochemical characteristics of Carboniferous coaly source rocks and natural gases in the Southeastern Junggar Basin, NW China: Implications for new hydrocarbon explorations
Wang et al. Stratigraphic thermohistory and its implications for regional geoevolution in the Tarim Basin, NW China
Zhao et al. Hydrocarbon-generating potential of the upper Paleozoic section of the Huanghua Depression, Bohai Bay Basin, China
Wang et al. Organic geochemical and petrographic characteristics of the coal measure source rocks of Pinghu Formation in the Xihu Sag of the East China Sea Shelf Basin: Implications for coal measure gas potential
Zhu et al. Identification of effective source rocks in different sedimentary environments and evaluation of hydrocarbon resources potential: A case study of paleogene source rocks in the Dongpu Depression, Bohai Bay Basin
Goodwin et al. Diamondoids and basin modelling reveal one of the world's deepest petroleum systems, South Caspian Basin, Azerbaijan
Botterell et al. Geochemical advances in Arctic Alaska oil typing–North Slope oil correlation and charge history
Mazumdar Stable isotopic study of late Neoproterozoic-early Cambrian (?) sediments from Nagaur-Ganganagar basin, western India: Possible signatures of global and regional C-isotopic events
Liu et al. Evaluation of the tight oil “sweet spot” in the Middle Permian Lucaogou Formation (Jimusaer Sag, Junggar Basin, NW China): Insights from organic petrology and geochemistry
Cao et al. Astronomical cycles calibrated the sea-level sequence durations of Late Miocene to Pliocene in Qiongdongnan Basin, south China sea
Farouk et al. Source and thermal maturity assessment of the Paleozoic-Mesozoic organic matter in the Risha gas field, Jordan
Li et al. Discrimination of effective source rocks and evaluation of the hydrocarbon resource potential in Marsel, Kazakhstan

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