瓊東南盆地長昌凹陷火成巖侵入體對溫度場及烴源巖成熟度的影響

唐曉音，張功成，梁建設(shè)，楊樹春，饒 松，胡圣標(biāo)

1 中國科學(xué)院地質(zhì)與地球物理研究所 巖石圈演化國家重點實驗室，北京 100029

2 中海油研究總院，北京 100027

1 引 言

熱模擬可再現(xiàn)油氣盆地的熱演化歷史，在油氣生氣和儲集過程中，盆地的熱演化決定著烴類成熟度史.沉積盆地中巖漿侵入體是除斷裂活動帶、放射性元素富集區(qū)之外的最主要的局部熱源［1］，其對圍巖的熱影響一直以來在國內(nèi)外受到廣泛的重視［2－8］.王有孝、朱傳慶等［9－10］發(fā)現(xiàn)異常地?zé)釋N源巖的有機質(zhì)生烴過程和烴類成熟作用有明顯影響.Othman等［11］發(fā)現(xiàn)侵入體的存在會使烴源巖的生油窗深度變淺.馮喬等［12］認(rèn)為侵入體余熱可以使生油層早熟，并使其附近圍巖中的有機質(zhì)豐度降低.萬從禮等［13］認(rèn)為基性巖漿帶來大量的熱量和地幔物質(zhì)，并產(chǎn)生極高的壓力，加速了周圍烴源巖的成熟演化，且促使其排烴.侵入作用的熱效應(yīng)對油氣成藏既有有利的一面，即加速烴源巖的成熟和成烴作用［7，14－15］，也有不利的一面，如破壞形成的油藏或加速烴源巖進入過熟期［16－18］.

長昌凹陷位于南海北部大陸邊緣瓊東南盆地深水區(qū)，西鄰寶島凹陷，北接神狐隆起、順德凹陷（見圖1），其所在構(gòu)造帶是深水凹陷群中唯一的中央背斜構(gòu)造帶［19］.長昌凹陷內(nèi)發(fā)育始新統(tǒng)湖相烴源巖、下漸新統(tǒng)崖城組和上漸新統(tǒng)陵水組濱海沼澤相烴源巖，始新統(tǒng)烴源巖雖然尚未被鉆井所揭示，但據(jù)地震剖面解釋成果推測這套深淺湖相烴源巖應(yīng)該存在.崖城組和陵水組烴源巖總有機碳含量分別為0.81%和0.69%，H／C一般小于1.0，以Ⅲ干酪根為主［20］.王振峰等［21］分析認(rèn)為，長昌凹陷很可能是油氣兼生的凹陷，其潛在天然氣資源量約數(shù)千億方，為瓊東南盆地深水區(qū)東部首要鉆探有利區(qū)帶.根據(jù)中海油研究總院最新研究成果，長昌凹陷大體上由主洼、東洼、南洼組成，內(nèi)部分布著多個錐狀火成巖侵入體（見圖1），侵入年齡為17.32Ma.目前為止，關(guān)于長昌凹陷尤其是火成巖侵入體對長昌凹陷的影響的研究還處于空白.

本文擬基于最新的地震解釋成果和相關(guān)熱物性資料，通過有限元計算方法［22］評價了火成巖侵入體對區(qū)域溫度場的影響，并進一步結(jié)合地層冷卻史及Easy%Ro模型探討其對烴源巖成熟度的影響.這一結(jié)果對評價侵入體對研究區(qū)內(nèi)烴源巖生、排烴的影響，計算凹陷油氣資源量等有重要意義.

圖1 長昌凹陷構(gòu)造單元劃分（a）及地震測線位置圖（b）Fig.1 （a）Structural division of the Changchang Sag，（b）distribution of intrusions and location of the modelled seismic line

2 計算方法與原理

2.1 熱傳導(dǎo)基本方程

有限元法是求解熱傳導(dǎo)方程常用的數(shù)值方法之一，早已被運用于解決熱傳導(dǎo)問題［23］，評價區(qū)域構(gòu)造形態(tài)、高溫巖漿侵入和地形變化等對區(qū)域溫度場的影響［24－26］.其本質(zhì)上是將計算區(qū)域內(nèi)的連續(xù)溫度場離散成有限數(shù)量的溫度點，然后求出給定條件下這些溫度點上的溫度值來近似代表需要求解的溫度場.由于侵入體和圍巖的熱力學(xué)性質(zhì)及地質(zhì)的復(fù)雜性，我們假設(shè)侵入是瞬間的，散熱方式主要是熱傳導(dǎo)，侵入體和圍巖的熱擴散率相同.其基本方程為：

式中，T 為溫度（℃）；t為時間（s）；K 為熱導(dǎo)率（W／（m·K））；Q為巖石的熱產(chǎn)生值（μW／m3）；μ為熱容量，μ＝ρc，ρ為密度（kg／m3），c為常壓下的比熱值（J／（kg·K））.

本文所討論的邊界條件為：1）上邊界海底溫度為10℃；2）下邊界上的熱流值為45mW／m2（僅考慮侵入體的影響，屏蔽了形成侵入體時的高熱背景）［27］；3）左右邊界假定為絕熱邊界（即沒有熱傳遞）.

2.2 模型的建立

模型的建立基于凹陷內(nèi)一條穿越三個不同規(guī)模侵入體的地震測線，模擬剖面AB位置見圖1，長度為200km，高度為14km.侵入體（1）、（2）、（3）為錐狀體，在沉積基底面上的半徑分別為3km、6km、5km；高度分別為：9.2km、10km、3.7km.剖面特征見圖2.總體上具有瓊東南盆地雙層結(jié)構(gòu)，上層表現(xiàn)為披蓋式坳陷，下層則是斷陷和斷隆，具有典型的張性環(huán)境下斷陷盆地的構(gòu)造特征.根據(jù)地震反射特征及實際鉆進層位標(biāo)定可識別出T20、T30、T40、T50、T60、T70、T100反射界面，分別對應(yīng)于樂東組、鶯歌海組、黃流組、梅山組、三亞組、陵水組、崖城組、始新統(tǒng)的底部.研究區(qū)T60是破裂不整合面，代表了凹陷從斷陷向坳陷轉(zhuǎn)化的界面，界面之下發(fā)育了基底斷裂和漸新統(tǒng)斷裂體系，構(gòu)造活動強，且斷層控制沉積；界面之上斷層基本不活動，一般認(rèn)為是層序的“穩(wěn)定”充填期［28］.根據(jù)實測剖面情況，建立模型（模型形態(tài)見圖2右下角），并采用0.5km×0.5km網(wǎng)格，取時間步長為0.05Ma對其進行二維有限元熱模擬.

圖2 地震剖面及模擬模型簡圖Fig.2 Sketch showing the seismic profile and simulation model

2.3 模擬參數(shù)

巖漿的溫度隨組分的不同而不同［29］，目前，尚沒有文獻介紹瓊東南盆地新生代火山活動的特點及巖性組成.根據(jù)珠江口盆地鉆遇火山巖的情況來看，始新世早期以中酸性噴出巖和火山碎屑巖為主，中后期基性噴出巖逐漸增加，中新世以后，火山巖成分比較單一，全為玄武巖［30］.鄒和平等［31］認(rèn)為中酸性火山巖是由于古近紀(jì)的火成巖組合來自大陸地殼之下的富集型地幔的熔融，巖漿中的深源成分逐漸增多，說明巖石圈在古近系發(fā)生了明顯的減薄，軟流圈地幔抬升較高，因此，新近系巖漿中的深源成分增多.據(jù)張健等［32］分析，南海北部不同盆地由于所處的地質(zhì)構(gòu)造環(huán)境差異，地幔熱流不同，但總體上各坳陷地質(zhì)成因機制及其形成、演化和衰亡過程基本相似，所以深部地?zé)釥顩r具有相似的特點.結(jié)合鄒和平等［31］新近系巖漿中的深源成分增多的理論，推測瓊東南盆地可能具有類似的巖性組成和變化特征，因此區(qū)內(nèi)侵入體初始溫度取980℃.其它模擬參數(shù)見表1.

表1 模擬參數(shù)Table 1 Simulation parameters

3 侵入體對溫度場的影響

3.1 對圍巖溫度的影響

圖3 侵入體侵入后不同時刻的溫度場剖面（a）侵入0.2Ma后；（b）侵入0.5Ma后；（c）侵入1Ma后；（d）侵入2Ma后；（e）侵入5Ma后；（f）侵入10Ma后.Fig.3 Geothermal field along the modelled section for different times after intrusion（a）0.2Ma，（b）0.5Ma，（c）1Ma，（d）2Ma，（e）5Ma，（f）10Ma.

圖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 yca.7km）；（b）Bottom of intrusion 2（yca.2km）.

有機質(zhì)演化過程中化學(xué)反應(yīng)速率與溫度呈指數(shù)關(guān)系，而與時間只呈線性關(guān)系，所以有機質(zhì)的成熟度主要受溫度的控制［35］.本文侵入體對圍巖溫度的影響如圖3、圖4所示：侵入1Ma后，侵入體溫度迅速下降，溫度界線開始消失（見圖3c），2Ma后界線消失（見圖3d），5Ma以后對圍巖溫度的影響微弱（圖3e），10Ma以后侵入體與圍巖溫度完全一致（圖3f），與范桃園、張健的研究結(jié)論［1，36］相符；侵入體不同深度冷卻速率是不一樣的，以侵入體（2）（規(guī)模最大）為例，侵入體中上部短時間內(nèi)迅速冷卻，侵入0.2Ma后，侵入體中心的溫度降至750℃，降幅超過200℃（圖4a），而侵入體底部溫度在侵入0.2Ma后下降溫度小于100℃（圖4b）.由此可知，隨著深度的增加，侵入體對圍巖溫度的影響程度減小，但持續(xù)時間較長.

3.2 對熱流的影響

針對不同的地質(zhì)尺度，目前國內(nèi)外應(yīng)用于盆地?zé)崾返难芯糠椒ㄖ饕袃煞N：一是在巖石圈尺度上，根據(jù)盆地形成和演化的動力機制及地球內(nèi)部熱傳播方式，建立相應(yīng)的盆地數(shù)學(xué)模型，在已知或假定的初始和邊界條件下，通過調(diào)整模型參數(shù)，使得模擬計算結(jié)果擬合實際觀測的盆地構(gòu)造沉降史而確定盆地底部熱流史，進而結(jié)合盆地的埋藏史正演盆地的熱史，稱之為“構(gòu)造－熱演化法”.二是在盆地尺度上，利用盆地內(nèi)的沉積地層中記載古地溫信息的有機質(zhì)、礦物、流體等（稱作古溫標(biāo)或古地溫計）來反演地層的熱歷史；據(jù)“古溫計”直接反推出古地溫范圍，或據(jù)它們的熱動力學(xué)機理，構(gòu)造數(shù)學(xué)模型，定量反演地?zé)崾?，稱之為“古溫標(biāo)法”［37－38］，這些方法用于熱史重建各有優(yōu)劣.由于研究區(qū)內(nèi)目前為止還沒有相應(yīng)的“古溫標(biāo)”，所以熱史反演行不通.而基于 McKenzie［39］拉伸模型的構(gòu)造熱演化法只能計算得到基底熱流史，不能反映后期侵入體的熱流影響.

為了定量地評價火成巖侵入體對烴源巖成熟度的影響，我們分別選取了圍繞三個侵入體2km、5km的特殊點，分析了侵入體對這些人工井熱流史的影響.如圖5所示：規(guī)模最大的侵入體2對熱流的影響最大，短時間內(nèi)引起的熱流異?？沙^25mW／m2.個別人工井熱流異常呈現(xiàn)雙峰，可能是受到鄰近侵入體的影響.但總體來說，侵入體對熱流的影響與對圍巖溫度的影響一致，有顯著影響的時限不超過1Ma.

4 侵入體對烴源巖有機質(zhì)成熟度的影響

火成巖的異常熱作用對烴源巖的影響主要表現(xiàn)在使鏡質(zhì)體反射率出現(xiàn)異常，使有機質(zhì)成熟度增加，促使烴類的成熟、生成及排出.其影響程度與烴源巖距火成巖體的遠(yuǎn)近和火成巖體的規(guī)模有關(guān)［40－42］.

本文以崖城組烴源巖為例，分析不同規(guī)模侵入體對距其不同距離處有機質(zhì)成熟度的影響，基底熱流來自基于過長昌凹陷地震測線的構(gòu)造熱演化計算結(jié)果［43］，在此基礎(chǔ)上，根據(jù)前述侵入體對熱流的影響和持續(xù)時間，將其所產(chǎn)生的熱作用折算成等效附加大地?zé)崃?，定量分析侵入體對人工井烴源層有機質(zhì)成熟度Ro的影響.圖6、7、8分別表示不同熱史（考慮與不考慮侵入體對熱流的影響）正演得到的距三個侵入體2km、5km處人工井崖城組烴源巖有機質(zhì)成熟度史.可以明顯地看出，當(dāng)受熱時間持續(xù)1Ma時，侵入體1短時間內(nèi)對距其2km處崖城組烴源巖成熟度Ro的影響可達到1.1% （見圖6a），對距其5km處烴源巖成熟度Ro的影響很小，約為0.15% （見圖6b）；受熱時間持續(xù)1Ma時，侵入體2短時間內(nèi)對距其2km處崖城組烴源巖成熟度Ro的影響可達到1.6%（見圖7a）；對距其5km處烴源巖成熟度Ro的影響為0.4%（見圖7b）；受熱時間持續(xù)1Ma時，侵入體3短時間內(nèi)對距其2km處崖城組烴源巖成熟度Ro的影響為0.03%（見圖8a）；對距其5km處烴源巖成熟度Ro幾乎沒有影響（見圖8b）.總體說來，侵入體對烴源巖有機質(zhì)成熟度的影響隨侵入體規(guī)模、距侵入體的距離等因素不同而不同，侵入體3由于規(guī)模小，所以對人工井烴源巖有機質(zhì)成熟度Ro的影響微弱.

5 討論與結(jié)論

關(guān)于巖漿侵入體對盆地?zé)嵫莼吧鸁N影響的論述很多，但大多數(shù)只是從定性的角度討論，定量化較少［44－47］.張健等［1］根據(jù)中國北方沉積盆地 中巖漿 侵入的一般特點，利用二維有限單元法計算了巖漿侵入體對沉積盆地溫度結(jié)構(gòu)的影響和對烴類成熟度的影響，得出幾公里大小的侵入體熱影響可以在時間上持續(xù)數(shù)百萬年，在空間上擴散至十幾公里.范桃園等［36］將該方法用于模擬美國某盆地尺度約為3～4km的分侵入體對盆地?zé)嵫莼挠绊?，得出的結(jié)論為在巖漿侵入以后，最初約0.4百萬年內(nèi)，巖漿周圍一定范圍（大約5km內(nèi)）的巖石受到巖漿所釋放熱的顯著影響.隨著時間的演化，影響范圍逐漸擴大，同時影響程度卻逐漸降低.楊樹春等［27］也用此方法分析了潮汕坳陷火成巖侵入體的影響，并認(rèn)為研究區(qū)侵入體對烴源巖成熟度影響微弱.本文的模擬結(jié)果總體上符合前人的結(jié)論，只是根據(jù)研究區(qū)的情況不同，影響的時限、空間上有所不同.由于資料有限，本次模擬過程中，只考慮了巖漿的熱傳導(dǎo)過程，而未能考慮巖漿侵入過程中的熱流體活動的影響［48－49］.長昌凹陷的勘探開發(fā)還處于空白狀態(tài)，區(qū)內(nèi)實際鉆井資料缺乏，隨著今后實測烴源巖有機質(zhì)成熟度Ro等資料的豐富，本次的模擬結(jié)果將會得到更進一步的證明與完善.

圖6 距侵入體（1）2km處（a）和距侵入體（1）5km處（b）人工井的熱史及崖城組烴源巖有機質(zhì)成熟度史Ⅰ－背景熱史；Ⅱ－考慮火山影響的熱史；1－崖城組烴源巖基于區(qū)域背景熱史計算所得有機質(zhì)成熟度史；2－崖城組烴源巖基于受侵入體（1）影響熱史計算所得有機質(zhì)成熟度史.Fig.6 Thermal history and organic maturity for artificial wells（a）2km away from igneous intrusion No.1and（b）5km away from igneous intrusion No.1Ⅰ－Background thermal history；Ⅱ－Thermal history for artificial well 2km 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.

目前，幾乎沒有關(guān)于長昌凹陷侵入體的報道，本文根據(jù)最新的地震解釋成果、熱物性資料，利用二維有限元法計算了火成巖侵入體對凹陷區(qū)域溫度場與熱流史的影響，并進一步結(jié)合Easy%Ro模型探討了侵入體熱效應(yīng)對烴源巖成熟度史的影響（以崖城組為例），得出以下幾點認(rèn)識：

（1）侵入后1Ma以內(nèi)，侵入體溫度迅速下降；2Ma以后，侵入體對溫度場的影響微弱，5Ma以后幾乎沒有影響，10Ma以后，侵入體溫度與圍巖溫度一致.侵入體對熱流的影響與對溫度的影響一致，有明顯影響的時限不超過1Ma，距侵入體2km處人工井的熱流值增加5～25mW／m2.

（2）侵入體不同深度冷卻速率不相同，隨著深度的增加，侵入體對圍巖溫度的影響程度減小，但持續(xù)時間較長.

（3）研究區(qū)內(nèi)侵入體對烴源巖有機質(zhì)成熟度的影響隨侵入體的規(guī)模、距侵入體的距離不同而不同.當(dāng)受熱時間持續(xù)1Ma，侵入體1、2短時間內(nèi)對距其2km處崖城組烴源巖成熟度Ro的影響明顯，最高可使Ro值增加1.6%，對距其5km處崖城組烴源巖成熟度Ro的影響為0.4%.而侵入體3由于規(guī)模較小，對烴源巖成熟度的影響微弱.

本文在評價侵入體對溫度場、熱演化、烴源巖有機質(zhì)成熟度影響的時候，重點在于考慮單個火成巖侵入體的影響.如果區(qū)域內(nèi)侵入體分布密集，某些位置可能會受到多個侵入體的共同影響（見圖5中熱流雙峰），在評價侵入體對整個區(qū)域的影響時則需根據(jù)具體的地質(zhì)情況、侵入體的具體規(guī)模和位置等因素綜合考慮.

圖7 距侵入體（2）2km處（a）和距侵入體（2）5km處（b）人工井的熱史及崖城組烴源巖有機質(zhì)成熟度史Ⅰ－背景熱史；Ⅱ－考慮火山影響的熱史；1－崖城組烴源巖基于區(qū)域背景熱史計算所得有機質(zhì)成熟度史；2－崖城組烴源巖基于受侵入體（2）影響熱史計算所得有機質(zhì)成熟度史.Fig.7 Thermal history and organic maturity for artificial wells（a）2km away from igneous intrusion No.2and（b）5km away from igneous No.2Ⅰ－Background thermal history；Ⅱ－ Thermal history for artificial well 2km 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 距侵入體（3）2km處（a）和距侵入體（3）5km處（b）人工井的熱史及崖城組烴源巖有機質(zhì)成熟度史Ⅰ－背景熱史；Ⅱ－考慮火山影響的熱史；1－崖城組烴源巖基于區(qū)域背景熱史計算所得有機質(zhì)成熟度史；2－崖城組烴源巖基于受侵入體（3）影響熱史計算所得有機質(zhì)成熟度史.Fig.8 Thermal history and organic maturity for artificial wells（a）2km away from igneous intrusion No.3and（b）5km away from igneous intrusion No.3Ⅰ－Background thermal history；Ⅱ－ Thermal history for artificial well 2km 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.

（References）

［1］張健，石耀霖.沉積盆地巖漿侵入的熱模擬.地球物理學(xué)進展，1997，12（3）：55－64.Zhang J，Shi Y L.The thermal modeling of magma intrustion in sedimentary basins.Progress in Geophysics（in Chinese），1997，12（3）：55－64.

［2］Zhu C Q，Xu M，Yuan Y S，et al.Palaeogeothermal response and record of the effusing of Emeishan basalts in the Sichuan basin.Chinese Science Bulletin，2010，55（10）：949－956.

［3］Sun Y G，F(xiàn)u J M，Liu D H，et al.Effect of volcanism on maturation of sedimentary organic matter and its significance for hydrocarbon generation，a case：the East Sag of Liaohe Basin.Chinese Science Bulletin，1995，40（17）：1446－1450.

［4］Wang D Y，Song Y C，Liu Y，et al.The influence of igneous intrusions on the peak temperatures of host rocks：Finite－time emplacement，evaporation，dehydration，and decarbonation.Computers ＆Geosciences，2012，38（1）：99－106.

［5］Barker C E，Boney Y，Lewan M D.Fluid inclusion and vitrinite－reflectance geothermometry compared to heat－flow models of maximum paleotemperaure next to dikes，western onshore Gippsland Basin，Australia.International Journal of Coal Geology，1998，37（1－2）：73－111.

［6］Bishop A N，Abbott G D.Vitrinite reflectance and molecular geochemistry of Jurassic sediments－the influence of heating by Tertiary dykes（Northwest Scotland）.Organic Geochemistry，1995，22（1）：165－177.

［7］Fjeldskaar W，Helset H M，Johansen H，et al.Thermal modelling of magmatic intrusions in the Gjallar Ridge，Norwegian Sea：implications for vitrinite reflectance and hydrocarbon maturation.Basin Research，2008，20（1）：143－159.

［8］Santos R V，Dantas E L，de Oliveira C G，et al.Geochemical and thermal effects of a basic sill on black shales and limestones of the Permian Irati Formation.Journal of South American Earth Sciences，2009，28（1）：14－24.

［9］王有孝，范璞，程學(xué)惠等.異常地?zé)釋Τ练e有機質(zhì)生烴過程的影響——以輝綠巖侵入體為例.石油與天然氣地質(zhì)，1990，11（1）：73－77.Wang Y X，F(xiàn)an P，Cheng X H，et al.Influence of diabase intrusion on hydrocarbon generation of depositional organic matter.Oil ＆Gas Geolog （in Chinese），1990，11（1）：73－77.

［10］朱傳慶，田云濤，徐明等.峨眉山超級地幔柱對四川盆地?zé)N源巖演化的影響.地球物理學(xué)報，2010，53（1）：119－127.Zhu C Q，Tian Y T，Xu M，et al.The effect on thermal evolution of source rock in Sichuan basin by Emeishan supper mantle plume.Chinese Journal of Geophysics（in Chinese），2010，53（1）：119－127.

［11］Othman R，Arouri K R，Ward C R，et al.Oil generation by igneous intrusions in the northern Gunnedah Basin，Australia.Organic Geochemistry，2001，32（10）：1219－1232.

［12］馮喬，湯錫元.巖漿活動對油氣藏形成條件的影響.地質(zhì)科技情報，1997，16（4）：59－65.Feng Q，Tang X Y.Magma activity′s influence of conditions forming oil and gas pools.Geological Science and Technology Information （in Chinese），1997，16（4）：59－65.

［13］萬從禮，金強.東營凹陷純西輝長巖對烴源巖異常生排烴作用研究.長安大學(xué)學(xué)報（地球科學(xué)版），2003，25（1）：20－25.Wan C L，Jin Q.Study on exceptional hydrocarbons generating and eliminating of gabbros to source rocks in Chunxi area of Dongying depression.Journal of Chang′an University （Earth Science Edition）（in Chinese），2003，25（1）：20－25.

［14］Chen R，He S，Wang Q.A preliminary discussion of magma activity on the maturation of organic matter：Taking Ggyucheng－Wenan area of Hebei Province as an example.Petroleun Expoloration and Development，1989，16（1）：29－38.

［15］Galushkin Y I.Thermal effects of igneous intrusions on maturity of organic matter：A possible mechanism of intrusion.Organic Geochemistry，1997，26（11－12）：645－658.

［16］Robinson N，Parnell J，Brassell S.Hydrocarbon compositions of bitumens from mineralized Devonian lavas and carboniferous sedimentary－rocks，central Scotland.Marine and Petroleum Geology，1989，6（4）：316－323.

［17］Finkelman R B，Bostick N H，Dulong F T，et al.Influence of an igneous intrusion on the inorganic geochemistry of a bituminous coal from Pitkin County，Colorado.International Journal of Coal Geology，1998，36（3－4）：223－241.

［18］Thorpe A N，Senftle F E，F(xiàn)inkelman R B，et al.Change in the magnetic properties of bituminous coal intruded by an igneous dike，Dutch Creek Mine，Pitkin County，Colorado.International Journal of Coal Geology，1998，36（3－4）：243－258.

［19］陶維祥，趙志剛，何仕斌等.南海北部深水西區(qū)石油地質(zhì)特征及勘探前景.地球?qū)W報，2005，26（4）：359－364.Tao W X，Zhao Z G，He S B，et al.Petroleum geological conditions and exploration prospects in deepwater area of northwestern South China Sea.Acta Geosicientia Sinica （in Chinese），2005，26（4）：359－364.

［20］龔再升，李思田，謝泰俊等.南海北部大陸邊緣盆地分析與油氣聚集.北京：科學(xué)出版社，1997.Gong Z S，Li S T，Xie T J，et al.Dynamic Research of Oil and Gas Accumulation in Northern Marginal Basions of South China Sea（in Chinese）.Beijing：Science Press，1997.

［21］王振峰，李緒深，孫志鵬等.瓊東南盆地深水區(qū)油氣成藏條件和勘探潛力.中國海上油氣，2011，23（1）：7－13.Wang Z F，Li X S，Sun Z P，et al.Hydrocarbon accumulation conditions and exploration potential in the deepwater region，Qiongdongnan basin.China Offshore Oil and Gas（in Chinese），2011，23（1）：7－13.

［22］張菊明，熊亮萍.有限單元法在地?zé)嵫芯恐械膽?yīng)用.北京：科學(xué)出版社，1986.Zhang J M，Xiong L P.The Application of Finite Element Methods in Geothermal（in Chinese）.Beijing：Science Press，1986.

［23］Bonneville A，Capolsini P.THERMIC：a 2－D finite－element tool to solve conductive and advective heat transfer problems in Earth Sciences.Computers ＆Geosciences，1999，25（10）：1137－1148.

［24］張菊明，孫惠文，熊亮萍.區(qū)域地溫場的數(shù)學(xué)模擬及實例剖析.地質(zhì)科學(xué)，1982，（3）：315－321.Zhang J M，Sun H W，Xiong L P.Mathematical simulation of regional geothermal field and case analysis.Chinese Journal of Geology （in Chinese），1982，（3）：315－321.

［25］魏東平，周蕙蘭.渤海唐山地區(qū)的地下熱結(jié)構(gòu)及地表熱流分布的數(shù)值模擬.地震地質(zhì)，1991，13（3）：283－288.Wei D P，Zhou H L.Mathematical simulation of the thermal structure and heat flow distribution in the Bohai and Tangshan areas.Seismology and Geology （in Chinese），1991，13（3）：283－288.

［26］田華，張之立.渤海及其鄰區(qū)地殼熱結(jié)構(gòu)和熱應(yīng)力場的研究.地震學(xué)報，1992，14（1）：29－35.Tian H，Zhang Z L.Study on Bohai Sea and its adjacent crustal thermal structure and thermal stress field.Acta Seismologica Sinica （in Chinese），1992，14（1）：29－35.

［27］楊樹春，仝志剛，賀清等.潮汕坳陷中生界生烴歷史及火成巖侵入影響分析——以LF35－1－1井為例.中國海上油氣，2008，20（3）：152－156.Yang S C，Tong Z G，He Q，et al.Mesozoic hydrocarbon generation history and igneous intrusion impacts in Chaoshan depression，South China sea：a case of LF35－1－1well.China Offshore Oil and Gas（in Chinese），2008，20（3）：152－156.

［28］吳時國，孫啟良，吳拓宇等.瓊東南盆地深水區(qū)多邊形斷層的發(fā)現(xiàn)及其油氣意義.石油學(xué)報，2009，30（1）：22－26.Wu S G，Sun Q L，Wu T Y，et al.Polygonal fault and oilgas accumulation in deep－water area of Qiongdongnan Basin.Acta Petrolei Sinica （in Chinese），2009，30（1）：22－26.

［29］Touloukian Y S，Judd W R，Roy R F.Physical Properties of Rocks and Minerals.New York：McGraw－Hill，1981.

［30］陳長民，施和生，許仕策.珠江口盆地（東部）第三系油氣藏形成條件.北京：科學(xué)出版社，2003.Chen C M，Shi H S，Xu S C.Forming Condition of Tertiary Reservoirs in Pearl River Mouth Basin（East）（in Chinese）.Beijing：Science Press，2003.

［31］鄒和平，李平魯，饒春濤.珠江口盆地新生代火山巖地球化學(xué)特征及其動力學(xué)意義.地球化學(xué)，1995，24（增刊）：33－45.Zhou H P，Li P L，Rao C T.Geochemistry of Cenozoic volcanic rocks in Zhu Jiangkou Basin and its geodynamics significance.Geochimica （in Chinese），1995，24（Suppl.）：33－45.

［32］張健，熊亮萍，汪集旸.南海深部地球動力學(xué)特征及其演化機制.地球物理學(xué)報，2001，44（5）：602－610.Zhang J，Xiong L P，Wang J Y.Characteristics and mechanism of geodynamic evolution of the South China Sea.Chinese Journal of Geophysics （in Chinese），2001，44（5）：602－610.

［33］袁玉松.南海北部陸緣深水區(qū)構(gòu)造－熱演化與烴源巖成熟度史研究［博士論文］.北京：中國科學(xué)院地質(zhì)與地球物理研究所，2007.Yuan Y S.Tectonic－thermal evolution and source rock′s organic maturity history of the deepwater area in the Northern Margin of the South China Sea［Ph.D.thesis］（in Chinese）.Beijing：Institute of Geology and Geophysics，Chinese Academy of Sciences，2007.

［34］Shi X B，Zhou D，Qiu X L，et al.Thermal and rheological structures of the Xisha Trough，South China Sea.Tectonophysics，2002，351（4）：285－300.

［35］Tissot B P，Pelet R，Ungerer P.Thermal history of sedimentary basins，maturation indexes，and kinetics of oil and gas generation.AAPG Bulletin，1987，71（12）：1445－1466.

［36］范桃園，石耀霖，周炎如.沉積盆地?zé)嵫莼^程中的巖漿作用.中國科學(xué)院研究生院學(xué)報，1999，16（1）：63－69.Fan T Y，Shi Y L，Zhou Y R.Effect of magma intrusion on the thermal evolution in a sedimentary basin.Journal of Graduate School，Academia Sinica （in Chinese），1999，16（1）：63－69.

［37］胡圣標(biāo)，汪集旸.沉積盆地?zé)狍w制研究的基本原理和進展.地學(xué)前緣，1995，2（4）：171－180.Hu S B，Wang J Y.Principles and progresses on thermal regime of sedimentary basins－an overview.Earth Science Frontiers（in Chinese），1995，2（4）：171－180.

［38］何麗娟，熊亮萍，汪集旸等.鶯歌海盆地構(gòu)造熱演化模擬研究.中國科學(xué)（D輯：地球科學(xué)），2000，30（4）：415－419.He L J，Xiong L P，Wang J Y，et al.Tectono－thermal modeling of Yinggehai basin.Science in China （Series D）（in Chinese），2000，30（4）：415－419.

［39］McKenzie D.Some remarks on the development of sedimentary basins.Earth and Planetary Science Letters，1978，40（1）：25－32.

［40］陳榮書，王青玲，何生等.巖漿活動對有機質(zhì)成熟作用的影響初探——以冀中葛漁城—文安地區(qū)為例.石油勘探與開發(fā)，1989，16（1）：29－38.Chen R S，Wang Q L，He S，et al.Primary exploration of influence of the magmatic activity on organic maturation：Case study on the Geyucheng－Wen′an area of central Hebei.Petroleum Exploration and Development （in Chinese），1989，16（1）：29－38.

［41］孫永革，傅家謨，劉德漢等.火山活動對沉積有機質(zhì)演化的影響及其油氣地質(zhì)意義——以遼河盆地東部凹陷為例.科學(xué)通報，1995，40（11）：1019－1022.Sun Y G，F(xiàn)u J M，Liu D H，et al.Effect of volcanism on maturation of sedimentary organic matter and its significance for hydrocarbon generation，a case：The east sag of Liaohe Basin.Chinese Science Bulletin（in Chinese），1995，40（11）：1019－1022.

［42］Raymond A C，Murchison D G.Effect of volcanic activity on level of organic maturation in carboniferous rocks of east Fife.Midland Valley of Scotland Fuel，1988，67（8）：1164－1166.

［43］單競男.南海北部深水區(qū)熱體制及烴源巖成熟度史研究［博士論文］.北京：中國科學(xué)院地質(zhì)與地球物理研究所，2011.Shan J N.Thermal regime and source rock maturation history in deepwater area，the Northern Margin of the South China Sea［Ph.D.thesis］（in Chinese）.Beijing：Institute of Geology and Geophysics，Chinese Academy of Sciences，2011.

［44］Jaeger J C.The temperature in the neighborhood of a cooling intrusive sheet.American Journal of Science，1957，255（4）：306－318.

［45］Simoneit B R T，Brenner S，Peters K E，et al.Thermal alteration of Cretaceous black shale by basaltic intrusions in the Eastern Atlantic.Nature，1978，273（5663）：501－504.

［46］Mongelli F，Zito G，de Lorenzo S.Thermal effects of long intrusions.Geothermics，2000，29（3）：347－365.

［47］Jaeger J C.Thermal effects of intrusions.Reviews of Geophysics，1964，2（3）：443－466.

［48］張健，石耀霖.熱泉地區(qū)巖漿侵入的熱演化特征.中國科學(xué)院研究生院學(xué)報，1998，15（2）：171－176.Zhang J，Shi Y L.The characteristic of thermal evolution of magma intrusion in hot spring area.Journal of Graduate School，Academia Sinica （in Chinese），1998，15（2）：171－176.

［49］石耀霖，王其允.沉積載荷致密排水對沉積盆地溫度結(jié)構(gòu)和成油窗口的影響.石油學(xué)報，1992，13（3）：19－26.Shi Y L，Wang Q Y.Geotherms and oil generation window in sedimentary basins thermal effects of compaction and discharge.Acta Petrolei Sinica （in Chinese），1992，13（3）：19－26.