Geochemical Evidences of Hydrocarbon Expulsion Threshold and Its Application
-
摘要: 排烃门限概念的提出为有效源岩的判识奠定了理论基础, 它在油气勘探与资源评价中具有重要意义.采用地球化学的方法论证了排烃门限的客观存在.研究表明, 源岩进入排烃门限后, 其反映生烃潜能的指标(S1+S2)/TOC、氢指数以及反映烃源岩残留烃量的指标氯仿沥青“A”都开始明显降低; 烃源岩中可溶有机质的化学成分、生物标志物的含量与分布特征也都发生明显的变化, 没有进入排烃门限的源岩可溶有机质与原油差别大, 进入门限后的源岩可溶有机质开始与原油组分较为一致.渤海湾盆地东营凹陷沙四段泥岩和页岩的排烃门限分别为2 200 m和2 000 m左右, 排烃强度分别为0~1 100 kg/m2和0~150 kg/m2, 累积排出烃量为7.365×108 t和1.435×108 t.Abstract: The concept of "Hydrocarbon expulsion threshold" is taken as an important theoretical base in the identification of effective source rock, and therefore, plays a significant role in petroleum exploration and resource assessment. Occurrence of hydrocarbon expulsion threshold is discussed by means of geochemical approaches. The results show that the parameters of (S1+S2)/TOC, HI showing source rock potential and the abundance of resident hydrocarbons-"A" decrease obviously with increasing depth once source rocks reach the expulsion threshold, and obvious variation of chemical composition as well as (abundance) and relative distribution of biomarkers can also be observed. Similar chemical compositions between crude oils and soluble organic matters are observed after the source rocks entered the expulsion threshold, however, the discrepancies are present within the source rocks below the expulsion threshold. Taking as an example of the Dongying depression of the Bohaiwan basin, the expulsion thresholds of Es4 mudstone and shale in Bamianhe area are about 2 200 and 2 000 m respectively, with expulsion intensity of 0-1 100 and 0-150 kg/m2 and expulsed accumulative hydrocarbons of 7.365×108 and 1.435×108 t.
-
众所周知, Tissot and Welte (1984)提出的干酪根晚期降解成烃理论已被几十年来的油气勘探实践所证实, 其在指导全球油气勘探实践中发挥了巨大的作用.近年来, 油气成烃理论飞速发展, 其中, 较为引人注目的有未熟-低熟油理论, 被认为是对干酪根晚期降解理论的补充与完善(张林晔和张春荣, 1999).自世界各地未熟-低熟油的发现与报道以来(Davis and Bery, 1969; Snowdon, 1980; Snowdon and Powell, 1982), 该理论的学术以及应用价值一直为广大石油地质工作者所重视.本文研究排烃门限存在的地球化学证据, 无论是未熟油还是成熟油, 它们生成后只有排出源岩外才能对油气藏形成作贡献.排烃门限概念的提出和确定为生烃理论走向实际应用架设了桥梁(庞雄奇等, 1993, 2001; 庞雄奇, 1995), 已有的研究从地质、物理模拟、数值模拟3方面论证了排烃门限的存在, 探索了它们在油气勘探中的应用; 并探索应用排烃门限理论确定有效源岩分布范围的方法和途径.
1. 排烃门限的基本概念与研究意义
排烃门限指烃源岩在埋藏演化过程中, 由于生烃量满足了自身吸附、孔隙水溶、油溶(气)和毛细管饱和等多种形式的残留需要并开始以游离相大量排出的临界点(图 1), 该点亦是烃源岩在演化过程中从欠饱和烃到过饱和烃, 从只能以水溶、扩散相排烃到能以游离相等多种形式排烃, 从少量排烃到大量排烃的转折点(庞雄奇等, 1993, 2001; 庞雄奇, 1995), 油气只有排出源岩后才能对油气成藏作贡献.排烃门限并非一个概念性名词, 而是可以加以定性与定量的地质指标, 已从物理模拟实验和数值模拟诸方面得到证据(图 2), 还应从微观的地球化学角度寻找直接的证据.研究排烃门限及其控油作用对于判别源岩、计算排烃量和确定有利勘探区带都具重要意义.
图 1 排烃门限地质概念模型(据庞雄奇等, 2001)Qp.累积生烃量; Qrm.残留烃临界饱和量; Qr.实际残烃量; Qe.源岩排烃总量; Qew, Qed, Qes.源岩水溶相、扩散相、游离相排烃量Fig. 1. Geological model for hydrocarbon expulsion threshold2. 排烃门限存在的地球化学证据
按照“干酪根晚期热降解作用生烃”的理论模式, 生烃门限的确定一般依据镜质体反射率Ro值是否达到0.5%为临界点.类似地, 排烃门限同样可用某些确定的指标来加以定性与定量.传统的排烃门限的识别方法是依据模拟实验结果、物质平衡计算结果、生烃门限和粘土矿物大量转化脱水门限等(庞雄奇等, 1993), 较偏重于地质因素.事实上, 排烃门限还有显著的地球化学证据可循, 该点以往未曾引起人们足够的重视.
2.1 源岩生烃潜力变化特征反映排烃门限的存在
在热解参数中, S1和S2分别代表游离烃、热解烃, “S1+S2”反映烃源岩生烃潜能, 包含残存于烃源岩中的烃以及尚未生成的烃.TOC可分为3部分, 一部分为“死碳”; 一部分以烃的形式已排出(若烃源岩已越过排烃门限); 第三部分仍存在于烃源岩中并以S1和S2形式存在. (S1+S2) /TOC为当前烃源岩的生烃潜力指数.假设有机碳是最初的固定值, 则(S1+S2) /TOC必存在一最大值; 如果TOC为未恢复值, 则相对于“S1+S2”, 其变化率也是相对较小的, 不论有机碳恢复与否, 实例研究表明, (S1+S2) /TOC值都表现出随埋深增大先增加, 然后达到某一极大值后再相对降低的趋势.统计表明, (S1+S2) /TOC曲线一般为近似抛物线的“大肚子”曲线, 特定的(S1+S2) /TOC曲线形状可因研究区甚至岩性而异.本研究所采集的八面河地区(东营凹陷牛庄洼陷) 烃源岩样(少部分为收集样) 如图 3a, 3b所示.类似地, 笔者对苏北盆地金湖凹陷进行研究后也得到类似的曲线(图 3c).生烃潜力“S1+S2”、TOC或(S1+S2)/TOC值的降低都意味着烃源岩已达排烃门限并开始排出烃类.除总有机碳(TOC) 与生烃潜量指标外, 氢指数(HI) 的变化也程度不等地与排烃作用有关(图 3d).依据烃源岩生烃潜能指标, 八面河地区沙四段页岩的排烃门限在2 000 m左右; 泥岩为2 200 m左右.
2.2 源岩残留烃量变化特征反映排烃门限的存在
氯仿沥青“A”反映烃源岩中的残存有机质, 理论上, 在热演化过程中, 氯仿沥青“A”值有逐渐增大然后减少的趋势.在原始烃源岩有机质丰度和类型保持一致的前提下, 氯仿沥青“A”含量的降低表明烃类的排出(图 4).在早期演化阶段, 烃转化率“A”/TOC可能变化不大(张林晔和张春荣, 1999), 并且很多未熟-低熟烃源岩表现为富含氯仿沥青“A”. 此时, 可能恰恰反映烃源岩尚未到达排烃门限, 如果氯仿沥青“A”的化学成分与该区原油相距太远, 这种含可溶有机质的烃源岩也不一定是有效烃源岩.研究表明, 未熟-低熟烃源岩族组分富含NSO的非烃, 后者需进一步演化才转化为石油烃类.
2.3 油-源对比特征反映排烃门限存在
传统排烃门限的确认方法比较侧重于从宏观地质角度讨论烃源岩的排烃门限、排烃过程与特征(庞雄奇等, 1993, 2001; 庞雄奇, 1995), 事实上, 从分子级水平判断排烃门限更为精确.对于成熟度较低的烃源岩, 抽提物常表现为富杂原子的醇、酸及酯类, 与石油烃类相距甚远.只有当烃源岩中可溶有机质的化学组成与原油成分较为接近时, 应当才有较为充分的理由说明已达到烃源岩的生、排烃门限.可溶有机质族组成特征、总烃转化率、烃类化合物分布规律(丰度、特征性生物标志物异构化程度) 都应当成为确认烃源岩是否达到生、排烃门限的重要识别指标.抽提物族组分中高含量的非烃、总烃转化率过低、特征性生物标志物异构化程度过低、某些生物标志物浓度过高或过低, 都可能标志着烃源岩并未开始大量生烃, 此时很难进入大量排烃门限.
八面河地区烃源岩可溶物中的特征性烃类生物标志物的分布特征可清楚地反映烃源岩演化、生排烃门限与生排烃过程.甾烷异构化参数C29甾烷ααα20S/ (S+R) (图 5a)、αα/βα-三环萜烷(图 5b)、C31藿烷22S/(S+R) (图 5c) 以及标样定量结果(图 5d) 反映出烃源岩从未熟至成熟的演化过程(Li et al., 2003; Pang et al., 2003).烃源岩的排烃门限应该在其可溶有机质的化学组成与原油组成接近的深度点附近, 对八面河油田而言, 油气主要来自2 500~2 700 m以下, 尽管上文的生烃潜能等指标反映埋深超过2 000~2 200 m的烃源岩开始进入排烃门限, 但可能不足以形成一定规模的油藏, 因为进入排烃门限后不久排出的烃类不能即刻成藏, 它们因吸附、水溶、扩散等作用而损耗在运移途中.进一步分析表明, 其化学组成远远没有达到原油成分级别(Li et al., 2003; Pang et al., 2003).因此, 进入了排烃门限的源岩需要排出一定的油气后才能进入成藏门限(庞雄奇等, 2000).
3. 排烃门限理论实际应用
3.1 确认源岩分布范围
在判定出源岩的排烃门限后可以圈定出有效源岩的分布范围, 计算它们的排烃强度和排烃量, 从而为有利区带优选和资源潜力评价提供理论指导.依据排烃门限可有效地确认排烃烃源岩的空间展布, 并进一步寻找可能聚集油气的圈闭, 提高勘探精度.依据烃类化学成分对比, 确认八面河油田的油气与沙四段2 700 m以下层段有较好的相关性(图 5) (Li et al., 2003; Pang et al., 2003), 即2 700 m以下的烃源岩始为主力排烃烃源岩, 这一临界点确定以后, 就可避免在缺少与有效烃源岩连通的浅部地层寻找工业规模的油气藏, 并提高油气资源评价的精度.
3.2 计算排烃强度和排烃量
目前比较流行的一种观点是, 在正常生油门限之上, 烃源岩中的可溶有机质可以直接形成烃类, 此为未熟—低熟油成因的重要机理之一.但以这种生烃机理为主的烃源岩在早期未熟—低熟阶段能形成多少油气, 特别是排出量有多大才是关系到油气勘探方向的实质性问题.按照排烃门限理论, 只有当烃源岩生成的烃类满足了吸附、水溶、扩散与毛细管饱和等需要后, 烃类才能有效排出.因此, 适用于不同生烃机理的有效的排烃量计算方法是解决问题的关键.在现有多种定量评价排烃量的方法(庞雄奇, 1995; Lafargue et al., 1990; Cooles et al., 1986; Burnham and Braun, 1990; Ungerer, 1990) 中, 最为常用的是物质平衡法, 首先估测出生成的烃与残留的烃, 排出的烃即Q排=Q生-Q残.但该方法所得结果相当程度上取决于生烃量计算方法, 残烃量可直接测定, 而生烃量的计算常存在诸多问题(Lewan et al., 2002).
本研究试图采用一种简单可行的新方法, 计算模型如图 6所示.正如前文所述, S1代表烃源岩中的游离烃, S2代表热解烃, (S1+S2) /TOC (生烃潜率指数) 值反映总剩余烃的潜能, 该值的降低意味着烃类的排出以及排烃门限的客观存在.分析具有统计意义的、有代表性的处于不同演化阶段的同层段样品, 可获得标准的样品包络线(图 6中曲线1).最大值即为源岩在地史过程中的最大生烃潜力, 最大生烃潜力与源岩埋藏过程中的生烃潜力之差即为排出烃量, 计算模型为:
(1) (2) 上列式中, Qrm、Qpm、Qe分别表示单位重量有机母质(TOC) 残留生烃潜量(mg/g)、最大生烃潜量(mg/g) 及排烃强度(mg/m2); TOC、ρr、H、S分别表示源岩的母质丰度含量(小数)、密度(g/cm3)、厚度(m)和面积(m2); z、z0分别表示源岩埋深(m)和对应于排烃门限的深度(m); Qe表示源岩层在S范围内进入排烃门限后的排烃总量.
这里强调说明, 源岩层厚度(H)、母质丰度(TOC)、密度(ρr) 以及残留生烃潜量(Qrm)等均是埋藏深度(z) 的函数, 获得它们在单位面上的等值线图后, 即可利用式(1)和(2) 计算出排烃强度(Qe)和排烃总量(Qe).
该方法在八面河地区的应用中获得了较好的效果.图 7指示八面河油田所在的牛庄洼陷南斜坡沙四段页岩相对较为发育, 但其排烃强度远不及同层段泥岩; 沙四段泥岩、页岩的排烃中心分别在牛庄洼陷和广利洼陷.各层段排烃强度等值线图还显示, 沙三段油页岩的排烃中心在牛庄洼陷西侧, 沙三段泥岩的排烃中心在牛庄洼陷, 但沙三中、上段泥岩的排烃强度远小于沙三下段.沙四段泥岩的排烃强度为0~1 100 kg/m2, 页岩为0~150 kg/m2.
八面河地区各层段烃源岩的排烃量计算结果如表 1, 该区总排烃量约为17.745×108 t, 如以2 700 m (Ro=0.5%) 为未熟与成熟阶段的划分界限, 则未熟烃排量约为0.735×108 t, 成熟烃17.01×108 t.总排烃量最高的是沙四段泥岩, 为7.365×108 t, 沙三泥岩稍低(6.031×108 t), 沙三、四段页岩分别为2.914×108 t, 1.435×108 t.沙三下、中、上亚段泥岩的排烃量分别为4.238×108 t, 1.758×108 t及0.035×108 t.
表 1 不同层位烃源岩排烃量统计Table Supplementary Table Total hydrocarbons expulsed from diverse strata of Bamianhe area3.3 各类烃源岩相对贡献量评估
在多套烃源岩发育的地区, 需要对各类烃源岩的相对贡献量作出客观的评估.依据上述方法, 不同层段烃源岩在不同阶段的排烃特征是, 沙四段页岩在未熟阶段的相对排烃量最高, 占17.352%, 其次是沙四段泥岩(5.974%); 沙三段烃源岩在未熟阶段的相对贡献量极低(约1%) (表 1).显然, 未熟烃基本由沙四段烃源岩提供, 沙四段页岩和泥岩在未熟阶段的排烃量各占未熟阶段总排烃量的33.878%和59.864% (表 1).以往研究中, 发育于牛庄洼陷南斜坡的沙四段页岩被认为是该区主力烃源岩, 然而, 本次研究结果表明, 沙四段2 700 m以上页岩的排烃量仅占沙四段总排烃量的2.83%, 油源对比表明, 后者为八面河油田的主力烃源岩.总之, 八面河地区各类烃源岩的相对贡献为, 从岩性而言, 以泥岩为主, 排烃量占75.494%, 页岩占24.507%;从区块而言, 以牛庄洼陷为主要排烃区, 其次是广利洼陷, 相对贡献分别为78.62%, 21.4%;从层位而言, 沙三、沙四排烃量相近, 分别为49.589%, 50.411%;从成熟度而言, 成熟烃源岩的贡献占95.862%, 未熟烃源岩的贡献占4.138% (表 2).
表 2 八面河地区源岩排烃量相对贡献大小评价Table Supplementary Table Relative contribution of hydrocarbons expulsed from Bamianhe area值得提出的是, 排烃理论在油气资源评价中的应用仍然需要结合有机地球化学工作, 以地质背景为前提.尽管计算结果表明, 八面河地区沙三段泥岩有较高的贡献量, 但其形成的烃类可能主要局限于生油中心的透镜体油藏中(如牛庄油田), 并未有效地运移至八面河油田所在的南斜坡带.
4. 结论
(1) 地球化学研究结果表明, 源岩排烃门限客观存在.进入排烃门限前, 源岩生烃潜力大、残留烃量随生烃量增大而增加、圈闭中聚集的油与源岩中的烃有显著差异.进入排烃门限后, 源岩生烃潜力迅速减小、残留烃量降低、油源对比显示亲缘关系. (2) 八面河油田周边沙四段泥岩和页岩的排烃门限分别为2 200 m和2 000 m; 八面河油田的油与周边2 700 m以下的源岩对比有亲缘关系.这说明, 源岩排烃门限早于油气成藏门限. (3) 八面河地区沙三段和沙四段泥岩和页岩的排烃强度约为0~1 100 kg/m2; 它们的累积排烃总量约为17.745×108 t, 未熟与成熟阶段排烃比率为1∶23.155, 页岩和泥岩比率为1∶3.082, 沙三段与沙四段比率为1∶0.984.
-
图 1 排烃门限地质概念模型(据庞雄奇等, 2001)
Qp.累积生烃量; Qrm.残留烃临界饱和量; Qr.实际残烃量; Qe.源岩排烃总量; Qew, Qed, Qes.源岩水溶相、扩散相、游离相排烃量
Fig. 1. Geological model for hydrocarbon expulsion threshold
表 1 不同层位烃源岩排烃量统计
Table 1. Total hydrocarbons expulsed from diverse strata of Bamianhe area
表 2 八面河地区源岩排烃量相对贡献大小评价
Table 2. Relative contribution of hydrocarbons expulsed from Bamianhe area
-
[1] Burnham, A.K., Braun, R.L., 1990. Development of a detailed model of petroleum formation, destruction, and expulsion from lacustrine and marine source rocks. Organic Geochemistry, 16(1-3): 27-39. doi: 10.1016/0146-6380(90)90023-S [2] Cooles, G.P., Mackenzie, A.S., Quigley, T.M., 1986. Calculation of petroleum masses generated and expelled from source rocks. Organic Geochemistry, 10(1-3): 235-245. doi: 10.1016/0146-6380(86)90026-4 [3] Davis, B.P., Bery, E.E., 1969. Analysis of oil and cap rock from Chanllenger(Sigsbee) Knoll. In: Ewing, W.M., ed., Initial Reports of the Deep Sea Drilling Project, 1: 415-499. [4] Lafargue, E., Espitalie, J., Jacobsen, T., et al., 1990. Experimental simulation of hydrocarbon expulsion. Organic Geochemistry, 16(1-3): 121-131. doi: 10.1016/0146-6380(90)90032-U [5] Lewan, M.D., Henry, M.E., Higley, D.K., et al., 2002. Material-balance assessment of the new Albany-Chesterian petroleum system of the Illinois basin. AAPG Bulletin, 86 (5): 754-777. https://pubs.er.usgs.gov/publication/70024630 [6] Li, S.M., Li, M.W., Pang, X.Q., et al., 2003. Geochemistry of petroleum systems in the Niuzhuang south slope of Bohai Bay basin-Part 1: Source rock characterization. Organic Geochemistry, 34(3): 389-412. doi: 10.1016/S0146-6380(02)00210-3 [7] Pang, X.Q., 1995. Theory of expulsion threshold controlling oil distribution and its application. Petroleum Industry Press, Beijing, 15-133(in Chinese). [8] Pang, X.Q., Chen, Z.M., Chen, F.J., 1993. Digital simulation of tectonic, thermal, hydrocarbon generation, residence and expulsion history, and source rock quantitative evaluation. Geological Publishing House, Beijing, 10-85(in Chinese). [9] Pang, X.Q., Ian, L., Wang, Y.C., et al., 2001. Theoretical study and application of hydrocarbon expulsion threshold of coal measures. Petroleum Industry Press, Beijing, 270 (in Chinese). [10] Pang, X.Q., Jiang, Z.X., Li, J.Q., et al., 2000. Geologic thresholds in the process of forming oil and gas reservoir and their functions of controlling petroleum. Journal of Universityof Petroleum, China, 24(4): 53-58(in Chinese with English abstract). [11] Pang, X.Q., Li, M.W., Li, S.M., et al., 2003. Geochemistry of petroleum systems in the Niuzhuang south slope of Bohai Bay basin-Part 2: Evidence for significant contribution of mature source rocks to" immature oils" in the Bamianhe field. Organic Geochemistry, 34(7): 931-950. doi: 10.1016/S0146-6380(03)00032-9 [12] Snowdon, L.R., 1980. Resinite-A potential petroleum source in the Upper Cretaceous/Tertiary of the BeaufortMachenzie basin. In: Miall, A.D., ed., Facts and principles of world petroleum occurrence. Canadian Soc. of Petroleum Geologists Mem., 6: 509-521. [13] Snowdon, L.R., Powell, T.G., 1982. Immature oil and condensate-Modification of hydrocarbon generation model for terrestrial organic matter. AAPG Bulletin, 66(6): 775-788. https://pubs.geoscienceworld.org/aapgbull/article/66/6/775/37520/Immature-Oil-and-Condensate-Modification-of [14] Tissot, B.P., Welte, D.H., 1984. Petroleum formation and Occurrence. Springer-Verlag, New York, 23-78. [15] Ungerer, P., 1990. State of the art of research in kinetic modeling of oil formation and expulsion. Organic Geochemistry, 16(1-3): 1-25. doi: 10.1016/0146-6380(90)90022-R [16] Zhang, L.Y., Zhang, C.R., 1999. Origins of immature oils and related petroleum systems: Case studies from the southern slope of Niuzhuang sag, Jiyang depression. Petroleum Industry Press, Beijing, 35-110(in Chinese). [17] 庞雄奇, 1995. 排烃门限控油气理论与应用. 北京: 石油工业出版社, 15-133. [18] 庞雄奇, 陈章明, 陈发景, 1993. 含油气盆地地史、热史、生留排烃史数值模拟研究与烃源岩定量评价. 北京: 地质出版社, 10-85. [19] 庞雄奇, Ian, L., 王雅春, 等, 2001. 煤系源岩排烃门限理论研究与应用. 北京: 石油工业出版社, 270. [20] 庞雄奇, 姜振学, 李建青, 等, 2000. 油气成藏过程中的地质门限及其控油气作用. 石油大学学报, 24(4): 53-58. https://www.cnki.com.cn/Article/CJFDTOTAL-SYDX200004012.htm [21] 张林晔, 张春荣, 1999. 低熟油生成机理及成油体系-以济阳坳陷牛庄洼陷南部斜坡为例. 北京: 石油工业出版社, 35-110. -