热模拟可再现油气盆地的热演化历史，在油气生气和储集过程中，盆地的热演化决定着烃类成熟度史.沉积盆地中岩浆侵入体是除断裂活动带、放射性元素富集区之外的最主要的局部热源^[1]，其对围岩的热影响一直以来在国内外受到广泛的重视^[2-8].王有孝、朱传庆等^[9-10]发现异常地热对烃源岩的有机质生烃过程和烃类成熟作用有明显影响.Othman等^[11]发现侵入体的存在会使烃源岩的生油窗深度变浅.冯乔等^[12]认为侵入体余热可以使生油层早熟，并使其附近围岩中的有机质丰度降低.万从礼等^[13]认为基性岩浆带来大量的热量和地幔物质，并产生极高的压力，加速了周围烃源岩的成熟演化，且促使其排烃.侵入作用的热效应对油气成藏既有有利的一面，即加速烃源岩的成熟和成烃作用^{[7, 14-15]}，也有不利的一面，如破坏形成的油藏或加速烃源岩进入过熟期^[16-18].

长昌凹陷位于南海北部大陆边缘琼东南盆地深水区，西邻宝岛凹陷，北接神狐隆起、顺德凹陷（见图 1），其所在构造带是深水凹陷群中唯一的中央背斜构造带^[19].长昌凹陷内发育始新统湖相烃源岩、下渐新统崖城组和上渐新统陵水组滨海沼泽相烃源岩，始新统烃源岩虽然尚未被钻井所揭示，但据地震剖面解释成果推测这套深浅湖相烃源岩应该存在.崖城组和陵水组烃源岩总有机碳含量分别为0.81%和0.69%，H/C一般小于1.0，以Ⅲ干酪根为主^[20].王振峰等^[21]分析认为，长昌凹陷很可能是油气兼生的凹陷，其潜在天然气资源量约数千亿方，为琼东南盆地深水区东部首要钻探有利区带.根据中海油研究总院最新研究成果，长昌凹陷大体上由主洼、东洼、南洼组成，内部分布着多个锥状火成岩侵入体（见图 1），侵入年龄为17.32 Ma.目前为止，关于长昌凹陷尤其是火成岩侵入体对长昌凹陷的影响的研究还处于空白.

图 1

Fig. 1

图 1 长昌凹陷构造单元划分（a）及地震测线位置图（b） Fig. 1 (a)Structural division of the Changchang Sag, (b)distribution of intrusions and location of the modelled seismic line

本文拟基于最新的地震解释成果和相关热物性资料，通过有限元计算方法^[22]评价了火成岩侵入体对区域温度场的影响，并进一步结合地层冷却史及Easy%Ro模型探讨其对烃源岩成熟度的影响.这一结果对评价侵入体对研究区内烃源岩生、排烃的影响，计算凹陷油气资源量等有重要意义.

有限元法是求解热传导方程常用的数值方法之一，早已被运用于解决热传导问题^[23]，评价区域构造形态、高温岩浆侵入和地形变化等对区域温度场的影响^[24-26].其本质上是将计算区域内的连续温度场离散成有限数量的温度点，然后求出给定条件下这些温度点上的温度值来近似代表需要求解的温度场.由于侵入体和围岩的热力学性质及地质的复杂性，我们假设侵入是瞬间的，散热方式主要是热传导，侵入体和围岩的热扩散率相同.其基本方程为：

(1)

式中，T为温度（℃）；t为时间（s）；K为热导率（W/（m·K））；Q为岩石的热产生值（μW/m³）；μ为热容量，μ=ρc，ρ为密度（kg/m³），c为常压下的比热值（J/（kg·K））.

本文所讨论的边界条件为：1）上边界海底温度为10 ℃；2）下边界上的热流值为45 mW/m²（仅考虑侵入体的影响，屏蔽了形成侵入体时的高热背景）^[27]；3）左右边界假定为绝热边界（即没有热传递）.

模型的建立基于凹陷内一条穿越三个不同规模侵入体的地震测线，模拟剖面AB位置见图 1，长度为200km，高度为14km.侵入体（1）、（2）、（3）为锥状体，在沉积基底面上的半径分别为3km、6km、5km；高度分别为：9.2km、10km、3.7km.剖面特征见图 2.总体上具有琼东南盆地双层结构，上层表现为披盖式坳陷，下层则是断陷和断隆，具有典型的张性环境下断陷盆地的构造特征.根据地震反射特征及实际钻进层位标定可识别出T₂₀、T₃₀、T₄₀、T₅₀、T₆₀、T₇₀、T₁₀₀反射界面，分别对应于乐东组、莺歌海组、黄流组、梅山组、三亚组、陵水组、崖城组、始新统的底部.研究区T₆₀是破裂不整合面，代表了凹陷从断陷向坳陷转化的界面，界面之下发育了基底断裂和渐新统断裂体系，构造活动强，且断层控制沉积；界面之上断层基本不活动，一般认为是层序的“稳定”充填期^[28].根据实测剖面情况，建立模型（模型形态见图 2右下角），并采用0.5km×0.5km网格，取时间步长为0.05 Ma对其进行二维有限元热模拟.

图 2

Fig. 2

图 2 地震剖面及模拟模型简图 Fig. 2 Sketch showing the seismic profile and simulation model

岩浆的温度随组分的不同而不同^[29]，目前，尚没有文献介绍琼东南盆地新生代火山活动的特点及岩性组成.根据珠江口盆地钻遇火山岩的情况来看，始新世早期以中酸性喷出岩和火山碎屑岩为主，中后期基性喷出岩逐渐增加，中新世以后，火山岩成分比较单一，全为玄武岩^[30].邹和平等^[31]认为中酸性火山岩是由于古近纪的火成岩组合来自大陆地壳之下的富集型地幔的熔融，岩浆中的深源成分逐渐增多，说明岩石圈在古近系发生了明显的减薄，软流圈地幔抬升较高，因此，新近系岩浆中的深源成分增多.据张健等^[32]分析，南海北部不同盆地由于所处的地质构造环境差异，地幔热流不同，但总体上各坳陷地质成因机制及其形成、演化和衰亡过程基本相似，所以深部地热状况具有相似的特点.结合邹和平等^[31]新近系岩浆中的深源成分增多的理论，推测琼东南盆地可能具有类似的岩性组成和变化特征，因此区内侵入体初始温度取980 ℃.其它模拟参数见表 1.

表 1(Table 1)

表 1 模拟参数 Table 1 Simulation parameters 参数 单位 数值 热导率 W/(m·K) 2.5[33] 生热率 μW/m3 0.9~0.3[34] 热扩散系数 cm2/s 0.01~0.02[27] 比热容 J/(kg·K) 787.1[29] 侵入体密度 kg/m3 3010[29]

表 1 模拟参数 Table 1 Simulation parameters

有机质演化过程中化学反应速率与温度呈指数关系，而与时间只呈线性关系，所以有机质的成熟度主要受温度的控制^[35].本文侵入体对围岩温度的影响如图 3、图 4所示：侵入1 Ma后，侵入体温度迅速下降，温度界线开始消失（见图 3c），2 Ma后界线消失（见图 3d），5 Ma以后对围岩温度的影响微弱（图 3e），10 Ma以后侵入体与围岩温度完全一致（图 3f），与范桃园、张健的研究结论^{[1, 36]}相符；侵入体不同深度冷却速率是不一样的，以侵入体（2）（规模最大）为例，侵入体中上部短时间内迅速冷却，侵入0.2 Ma后，侵入体中心的温度降至750 ℃，降幅超过200 ℃（图 4a），而侵入体底部温度在侵入0.2 Ma后下降温度小于100 ℃（图 4b）.由此可知，随着深度的增加，侵入体对围岩温度的影响程度减小，但持续时间较长.

图 3

Fig. 3

图 3 侵入体侵入后不同时刻的温度场剖面（a）侵入0.2 Ma后；（b）侵入0.5 Ma后；（c）侵入1 Ma后；（d）侵入2 Ma后；（e）侵入5 Ma后；（f）侵入10 Ma后. Fig. 3 Geothermal field along the modelled section for different times after intrusion(a)0.2 Ma, (b)0.5 Ma, (c)1 Ma, (d)2 Ma, (e)5 Ma, (f)10 Ma.

图 4

Fig. 4

图 4 侵入体2侵入后温度随时间的变化（a）侵入体（2）中心（y约为7km）；（b）侵入体（2）底部（y约为2km）. Fig. 4 Temperature changing with time after intrusion for intrusion No.2(a)Intrusion center(at y ca.7km); (b)Bottom of intrusion 2(y ca.2km).

针对不同的地质尺度，目前国内外应用于盆地热史的研究方法主要有两种：一是在岩石圈尺度上，根据盆地形成和演化的动力机制及地球内部热传播方式，建立相应的盆地数学模型，在已知或假定的初始和边界条件下，通过调整模型参数，使得模拟计算结果拟合实际观测的盆地构造沉降史而确定盆地底部热流史，进而结合盆地的埋藏史正演盆地的热史，称之为“构造-热演化法”.二是在盆地尺度上，利用盆地内的沉积地层中记载古地温信息的有机质、矿物、流体等（称作古温标或古地温计）来反演地层的热历史；据“古温计”直接反推出古地温范围，或据它们的热动力学机理，构造数学模型，定量反演地热史，称之为“古温标法”^[37-38]，这些方法用于热史重建各有优劣.由于研究区内目前为止还没有相应的“古温标”，所以热史反演行不通.而基于McKenzie^[39]拉伸模型的构造热演化法只能计算得到基底热流史，不能反映后期侵入体的热流影响.

为了定量地评价火成岩侵入体对烃源岩成熟度的影响，我们分别选取了围绕三个侵入体2km、5km的特殊点，分析了侵入体对这些人工井热流史的影响.如图 5所示：规模最大的侵入体2对热流的影响最大，短时间内引起的热流异常可超过25mW/m².个别人工井热流异常呈现双峰，可能是受到邻近侵入体的影响.但总体来说，侵入体对热流的影响与对围岩温度的影响一致，有显著影响的时限不超过1Ma.

图 5

Fig. 5

图 5 距侵入体2km、5km处人工井热史侵入体（1）与侵入体（2）之间的人工井热流异常呈现双峰，可能是两个侵入体之间的相互影响. Fig. 5 Thermal history of artificial wells 2 km and 5 km away from the three igneous intrusions Bimodal heat flow anomalies for artificial wells may be caused by interaction between two igneous intrusions.

火成岩的异常热作用对烃源岩的影响主要表现在使镜质体反射率出现异常，使有机质成熟度增加，促使烃类的成熟、生成及排出.其影响程度与烃源岩距火成岩体的远近和火成岩体的规模有关^[40-42].

本文以崖城组烃源岩为例，分析不同规模侵入体对距其不同距离处有机质成熟度的影响，基底热流来自基于过长昌凹陷地震测线的构造热演化计算结果^[43]，在此基础上，根据前述侵入体对热流的影响和持续时间，将其所产生的热作用折算成等效附加大地热流，定量分析侵入体对人工井烃源层有机质成熟度Ro的影响.图 6、7、8分别表示不同热史（考虑与不考虑侵入体对热流的影响）正演得到的距三个侵入体2km、5km处人工井崖城组烃源岩有机质成熟度史.可以明显地看出，当受热时间持续1 Ma时，侵入体1短时间内对距其2km处崖城组烃源岩成熟度Ro的影响可达到1.1%（见图 6a），对距其5km处烃源岩成熟度Ro的影响很小，约为0.15%（见图 6b）；受热时间持续1 Ma时，侵入体2短时间内对距其2km处崖城组烃源岩成熟度Ro的影响可达到1.6%（见图 7a）；对距其5km处烃源岩成熟度Ro的影响为0.4%（见图 7b）；受热时间持续1 Ma时，侵入体3短时间内对距其2km处崖城组烃源岩成熟度Ro的影响为0.03%（见图 8a）；对距其5km处烃源岩成熟度Ro几乎没有影响（见图 8b）.总体说来，侵入体对烃源岩有机质成熟度的影响随侵入体规模、距侵入体的距离等因素不同而不同，侵入体3由于规模小，所以对人工井烃源岩有机质成熟度Ro的影响微弱.

图 6

Fig. 6

图 6 距侵入体（1）2km处（a）和距侵入体（1）5km处（b）人工井的热史及崖城组烃源岩有机质成熟度史Ⅰ-背景热史；Ⅱ-考虑火山影响的热史；1-崖城组烃源岩基于区域背景热史计算所得有机质成熟度史；2-崖城组烃源岩基于受侵入体（1）影响热史计算所得有机质成熟度史. Fig. 6 Thermal history and organic maturity for artificial wells(a)2 km away from igneous intrusion No.1 and (b) 5 km away from igneous intrusion No.1 Ⅰ-Background thermal history; Ⅱ-Thermal history for artificial well 2 km away from intrusion No.1;1-Organic maturity history of the Yacheng Formation based on the background thermal history; 2-Organic maturity history of the Yacheng Formation based on thermal history influenced by intrusion No.1.

图 7

Fig. 7

图 7 距侵入体（2）2km处（a）和距侵入体（2）5km处（b）人工井的热史及崖城组烃源岩有机质成熟度史Ⅰ-背景热史；Ⅱ-考虑火山影响的热史；1-崖城组烃源岩基于区域背景热史计算所得有机质成熟度史；2-崖城组烃源岩基于受侵入体（2）影响热史计算所得有机质成熟度史. Fig. 7 Thermal history and organic maturity for artificial wells(a)2 km away from igneous intrusion No.2 and(b)5 km away from igneous No.2 Ⅰ-Background thermal history; Ⅱ-Thermal history for artificial well 2 km away from intrusion No.2;1-Organic maturity history of the Yacheng Formation based on the background thermal history; 2-Organic maturity history of the Yacheng Formation based on thermal history influenced by intrusion No.2.

图 8

Fig. 8

图 8 距侵入体（3）2km处（a）和距侵入体（3）5km处（b）人工井的热史及崖城组烃源岩有机质成熟度史Ⅰ-背景热史；Ⅱ-考虑火山影响的热史；1-崖城组烃源岩基于区域背景热史计算所得有机质成熟度史；2-崖城组烃源岩基于受侵入体（3）影响热史计算所得有机质成熟度史. Fig. 8 Thermal history and organic maturity for artificial wells(a)2 km away from igneous intrusion No.3 and(b)5 km away from igneous intrusion No.3 Ⅰ-Background thermal history; Ⅱ-Thermal history for artificial well 2 km away from intrusion No.2;1-Organic maturity history of the Yacheng Formation based on the background thermal history; 2-Organic maturity history of the Yacheng Formation based on thermal history influenced by intrusion No.3.

关于岩浆侵入体对盆地热演化及生烃影响的论述很多，但大多数只是从定性的角度讨论，定量化较少^[44-47].张健等^[1]根据中国北方沉积盆地中岩浆侵入的一般特点，利用二维有限单元法计算了岩浆侵入体对沉积盆地温度结构的影响和对烃类成熟度的影响，得出几公里大小的侵入体热影响可以在时间上持续数百万年，在空间上扩散至十几公里.范桃园等^[36]将该方法用于模拟美国某盆地尺度约为3~4km的分侵入体对盆地热演化的影响，得出的结论为在岩浆侵入以后，最初约0.4百万年内，岩浆周围一定范围（大约5km内）的岩石受到岩浆所释放热的显著影响.随着时间的演化，影响范围逐渐扩大，同时影响程度却逐渐降低.杨树春等^[27]也用此方法分析了潮汕坳陷火成岩侵入体的影响，并认为研究区侵入体对烃源岩成熟度影响微弱.本文的模拟结果总体上符合前人的结论，只是根据研究区的情况不同，影响的时限、空间上有所不同.由于资料有限，本次模拟过程中，只考虑了岩浆的热传导过程，而未能考虑岩浆侵入过程中的热流体活动的影响^[48-49].长昌凹陷的勘探开发还处于空白状态，区内实际钻井资料缺乏，随着今后实测烃源岩有机质成熟度Ro等资料的丰富，本次的模拟结果将会得到更进一步的证明与完善.

目前，几乎没有关于长昌凹陷侵入体的报道，本文根据最新的地震解释成果、热物性资料，利用二维有限元法计算了火成岩侵入体对凹陷区域温度场与热流史的影响，并进一步结合Easy%Ro模型探讨了侵入体热效应对烃源岩成熟度史的影响（以崖城组为例），得出以下几点认识：

（1）侵入后1 Ma以内，侵入体温度迅速下降；2 Ma以后，侵入体对温度场的影响微弱，5 Ma以后几乎没有影响，10 Ma以后，侵入体温度与围岩温度一致.侵入体对热流的影响与对温度的影响一致，有明显影响的时限不超过1 Ma，距侵入体2km处人工井的热流值增加5~25mW/m².

（2）侵入体不同深度冷却速率不相同，随着深度的增加，侵入体对围岩温度的影响程度减小，但持续时间较长.

（3）研究区内侵入体对烃源岩有机质成熟度的影响随侵入体的规模、距侵入体的距离不同而不同.当受热时间持续1 Ma，侵入体1、2短时间内对距其2km处崖城组烃源岩成熟度Ro的影响明显，最高可使Ro值增加1.6 %，对距其5km处崖城组烃源岩成熟度Ro的影响为0.4%.而侵入体3由于规模较小，对烃源岩成熟度的影响微弱.

本文在评价侵入体对温度场、热演化、烃源岩有机质成熟度影响的时候，重点在于考虑单个火成岩侵入体的影响.如果区域内侵入体分布密集，某些位置可能会受到多个侵入体的共同影响（见图 5中热流双峰），在评价侵入体对整个区域的影响时则需根据具体的地质情况、侵入体的具体规模和位置等因素综合考虑.

本文研究使用地震测线及解释成果由中海石油（中国）有限公司北京研究中心提供.