在火山岩油气藏研究中,储层岩性、岩相分类、孔渗变化、火山作用对烃源岩成熟度影响以及无机成因烃的形成机理等方面研究较多,对于火山作用的有机质成烃热效应几乎没有研究。研究火山作用的有机质成烃热效应可以为评价其对烃源岩生烃量的增减和油藏的热破坏程度提供技术手段,还有助于正确认识火山作用在油气成藏中的贡献和完善火山岩油气藏理论。这方面研究不仅需要建立火山岩体的热容模型、热传导模型,还需建立有机质成烃动力学模型。为此,建立了多个化学动力学模型及岩浆侵入体热传导模型,并建立了各模型参数优化求取方法。进一步在不同类型有机质的生烃(油、气)热模拟实验,获得产物产率与温度的关系数据基础上,对不同模型、不同实验体系结果、不同类型有机质生烃动力学特征的进行对比研究,并开展了动力学模型参数稳定性、烃源岩非均质性对动力学参数的影响以及岩浆侵入体规模、性质、与烃源岩时空匹配关系不同时对围岩成熟度、生烃热效应的影响研究。最后结合徐家围子、英台断陷烃源岩参数进行天然气资源量评价。本次研究对今后应用生烃动力学方法评价资源潜力以及定量评价火山热作用对烃源岩生烃量的增减有重要指导意义。取得以下成果和认识:(1)密闭体系下有机质、原油热裂解过程中总烃气质量产率下降的温度与重烃气质量产率的拐点温度并不相同,表明高温阶段甲烷的来源仍然存在有机质/原油的初次裂解贡献,不完全是重烃气(C2-5)裂解贡献,但对于组成相对简单的液态烃或纯化合物这两个温度比较接近。(2)TG-MS实验下煤岩生甲烷终止温度约为850℃,对应的Ro约为5.3%(10℃/min升温速率),密闭体系下煤岩在650℃时(Ro约为4.9%,升温速率2℃/hour)煤岩生气能力尚未结束,气态烃质量产率一直呈增长趋势,表明煤岩在高温阶段仍具有生甲烷能力。可能是密闭体系下低温阶段热解的正构烷烃(n-alkyl)产物通过环化和芳香化作用与沥青或者干酪根发生缩聚/再结合作用形成了具有较高热稳定性的产物,这一产物在高温阶段可以再次生成甲烷。而在开放体系下正构烷烃被载气携走,不具备这类反应发生条件。(3)热传导模型结合EasyRo%模型可以较好的拟合实测Ro数据,还可以方便地模拟计算二维、三维空间内温度场及成熟度史演化。岩浆侵入体的热作用范围是有限的,地质情况不同,影响范围广度也不同,对于相同厚度的侵入体,初始温度越高,作用范围越广,但一般X/D<3(X/D指接触面的距离与侵入体厚度的比值)。对于不同厚度侵入体而言,侵入体热作用影响范围X/D<2。(4)总包反应动力学方程简单,但不适合描述复杂有机质成烃过程,有机质生烃也不完全符合连串反应机理,可能结合连串反应和平行反应两种模型描述有机质生烃过程比较合适,但是这类模型参数标定比较复杂,实验室工作比较繁重。串联反应模型无论从动力学理论上还是标定动力学参数的结果上都不适合。从各模型在烃源岩潜力评价中应用情况及各模型参数的优化效果来看,活化能服从离散分布的平行一级反应模型描述有机质生烃过程最为合适。(5)随指前因子的增加,表观活化能逐渐增加,且不同类型干酪根的表观活化能按如下顺序排列E II2-III >E I >E II1。(6)以SFF模型描述有机质生烃特征,优化的动力学参数—指前因子由小到大变化,其优化效果也存在一个极值。随着指前因子的增加,活化能分布形态基本上不发生变化,但整体向活化能增高方向偏移,指前因子每增加10倍,平均活化能增加约12kJ/mol,以3.3℃/Ma的地质升温速率进行外推,TR0.5(成烃转化率为50%)对应的地质温度逐渐增加,但增加值逐渐减少。(7)采用SFF模型研究不同类型有机质成烃动力学参数特征表明湖相I型有机质样品化合物结构及化学键类型相对单一,以单个活化能分布形态为主,湖相II2及陆相III型有机质样品化合物结构的非均质性较强,活化能分布范围宽,且“优势活化能”所占反应的分数随其氢指数降低而降低。同时不同活化能分布间隔动力学模型地质外推结果显示,活化能分布间隔对成烃转化率的影响不大,可以采用适当个数的平行反应模型描述有机质成烃过程。地质外推表明,I型有机质生烃速率范围分布较窄,且曲线圆滑。II1型、II2—III型有机质生烃速率分布范围较宽,其中II2—III型有机质生烃速率存在多个局部高点,波动频繁。(8)采用MFF(Multiple Frequence Factor)模型可以避免采用SFF(Single Frequence Factors)模型出现的低估具有较低和较高活化能有机质的生烃潜力(反应分数)这一问题。(9)烃源岩非均质性对采用生烃动力学方法评价资源量有一定的影响,可以通过采用MFF模型和S2加权平均法获得的动力学参数进行资源评价。同时建议结合地质上其他数据(如S1、Tmax、S1+S2等)对动力学参数进行优选。(10)相同条件的侵入体对于具有不同初始Ro的有机质热成熟度和生烃的热效应并不相同,表现在Ro<0.9%(对应大量生烃前)时,热效应随着初始Ro的增加而逐渐增加,当初始Ro>0.9%(对应大量生烃后)之后,热效应随着初始Ro的增加而逐渐降低。因此,具有相同初始条件的侵入体可以导致不同的热影响范围。在烃源岩大量生烃前,岩浆侵位时间越晚对生烃的热效应也越大,反之,则越小。(11)对龙深1井侵入体附近泥岩有机质成熟度及地球化学参数进行了分析,表明随着与接触面距离的减小,围岩中有机质成熟度逐渐升高,由原来的1.6变化到2.1左右,岩石热解参数Tmax则逐渐增大,H/C和O/C值逐渐降低,围岩中TOC逐渐降低,氯仿沥青"A"先增加后降低。在两套侵入体之间,Tmax值变化呈现”V”字型特征,围岩中TOC出现先增大后降低的趋势。同样,围岩氢指数(HI)、热解烃S2也出现类似变化。(12)采用MFF模型结合龙深1井区埋藏史、热史对侵入体的生烃热效应进行了研究,表明在侵入体影响范围内,成气转化率迅速增长。侵入体不同距离处的生气史则表明在很短的时间内,生气转化率从0增加到80%以上。可见,将化学动力学模型、埋藏史模型、侵入体热传导模型和正常热史模型结合起来可以很好的研究有机质复杂的成烃过程。(13)评价结果显示英台断陷深层天然气生成总量为5.1×1012m3,资源量为1072×108~1608×108m3(运聚系数取1.6%~2.4%)。徐家围子断陷深层天然气总生成量为33.75×1012m3,资源量为5020×108~7530×108m3(运聚系数取1.6%~2.4%)。