利用鏡質(zhì)體反射率方法尋找隱伏巖體
——巖漿熱場應(yīng)用的一個實(shí)例

張　旗， 金維浚， 李承東， 王金榮， 焦守濤，， 陳萬峰

（1.中國科學(xué)院 地質(zhì)與地球物理研究所， 巖石圈演化國家重點(diǎn)實(shí)驗(yàn)室， 北京 100029； 2.中國地質(zhì)調(diào)查局天津地質(zhì)調(diào)查中心， 天津 300170； 3.蘭州大學(xué) 地質(zhì)科學(xué)與礦產(chǎn)資源學(xué)院， 甘肅 蘭州 730000； 4.北京礦產(chǎn)地質(zhì)研究所， 北京 100012）

利用鏡質(zhì)體反射率方法尋找隱伏巖體
——巖漿熱場應(yīng)用的一個實(shí)例

張旗1， 金維浚1， 李承東2， 王金榮3， 焦守濤1，4， 陳萬峰3

（1.中國科學(xué)院 地質(zhì)與地球物理研究所， 巖石圈演化國家重點(diǎn)實(shí)驗(yàn)室， 北京 100029； 2.中國地質(zhì)調(diào)查局天津地質(zhì)調(diào)查中心， 天津 300170； 3.蘭州大學(xué) 地質(zhì)科學(xué)與礦產(chǎn)資源學(xué)院， 甘肅 蘭州 730000； 4.北京礦產(chǎn)地質(zhì)研究所， 北京 100012）

尋找隱伏巖體是熱液礦床找礦中一項(xiàng)重要的任務(wù)。識別隱伏巖體的方法很多， 如磁法、電法、重力法等， 很多很有效。但是， 缺少了一個既直接又簡便的方法： 鏡質(zhì)體反射率方法。巖漿侵位必然在周圍形成一個熱場， 且隨著距離巖漿的遠(yuǎn)近形成一個溫度梯度場， 疊加在地?zé)釄鲋?。因此?依據(jù)巖漿熱場的變化即可發(fā)現(xiàn)隱伏巖體。如果我們能夠知道隱伏巖漿的性質(zhì)、圍巖的熱傳導(dǎo)參數(shù)， 甚至可以定量地查明隱伏巖體的深度、位置和規(guī)模。本文介紹煤田和石油地質(zhì)部門最常用的鏡質(zhì)體反射率（vitrinite reflectance， 以Ro表示）方法， 具有簡單、經(jīng)濟(jì)、有效、實(shí)用、快捷的優(yōu)點(diǎn)。利用該方法識別隱伏巖體的有效距離根據(jù)巖體規(guī)模、成分、深度不同大約在幾米至幾公里范圍。無論新老礦區(qū)， 只要有志留紀(jì)以后的含泥質(zhì)的巖石， 都可以嘗試采用這種方法。該方法既適用于地表也適用于鉆孔。對于老礦區(qū)， 有條件開展三維立體層面上的研究； 新礦區(qū)鉆孔少， 應(yīng)注意了解鉆孔不同深度Ro變化的情況， 尤其注意鉆孔中某些深度Ro突變的情況以及隨鉆進(jìn)深度增加Ro變化的情況， 以便及時指導(dǎo)找礦。筆者認(rèn)為， 在礦區(qū)和非礦區(qū)， 隱伏巖體的意義是不一樣的。在礦區(qū)（無論新礦區(qū)老礦區(qū)）， 巖體與流體有相伴出現(xiàn)的可能性。如果在礦區(qū)發(fā)現(xiàn)Ro出現(xiàn)異常（Ro＞2%）， 則應(yīng)注意考察是否有與隱伏巖體相伴的礦化的可能性。筆者希望將該方法引入到礦床學(xué)研究中來， 并結(jié)合熱液金屬礦床成礦作用的特征， 開拓出新的應(yīng)用前景， 推進(jìn)深部找礦取得新的進(jìn)展。

鏡質(zhì)體反射率； 隱伏巖體； 巖漿熱場； 方法； 成礦； 花崗巖

卷（Volume）39， 期（Number）6， 總（SUM）149

頁（Pages）1094～1107， 2015， 12（December， 2015）

識別隱伏巖體是熱液礦床找礦中一項(xiàng)重要的任務(wù)。識別隱伏巖體的方法很多， 如磁法、電法、重力法等， 很多很有效， 王維和湯靜如（2013）對此有很好的歸納。但是， 缺少了一個既直接又簡便的方法： 巖漿熱場法， 即巖漿產(chǎn)生的熱場或熱效應(yīng)的方法。巖漿是熱的， 巖漿侵位， 必然在周圍形成一個熱場， 且隨著距離巖漿的遠(yuǎn)近， 溫度逐漸下降， 形成一個溫度梯度場， 疊加在地?zé)釄鲋稀Ｈ藗兗纯筛鶕?jù)這種熱場溫度的分布去發(fā)現(xiàn)隱伏的巖漿侵入體。如果我們能夠知道隱伏巖漿的性質(zhì)， 圍巖的熱傳導(dǎo)參數(shù)， 甚至可以定量地查明隱伏巖體的深度、位置和規(guī)模。這并不是什么新鮮的見解， 國外從20世紀(jì)50～60年代開始， 即非常注意對這個問題的研究（Carslaw and Jaeger， 1959； Jaeger， 1961， 1964；Feoktistov， 1972； Fedotov， 1976； Duffield and Ruiz，1992； Rubin， 1995； Webber et al.， 1999； Gunson et al.，2000； Ascencio et al.， 2006； Fjeldskaar et al.， 2008；Eldursi et al.， 2009； Kol'tsov， 2010）， 國內(nèi)這項(xiàng)研究開展得也比較早， 尤以煤田和石油部門研究極為出色（如： 楊起和任德貽， 1981； 楊文寬， 1982； 楊起等，1987a， b， 1996； 陳榮書等， 1989； 楊起， 1989， 1999；陳儒慶和袁奎榮， 1990； 吳傳榮， 1992； 吳江和李思田， 1993； 馬昌前等， 1994； 張正階等， 1996； 羅文積和陳家清， 1997； 馮喬和湯錫元， 1997； 張健和石耀霖， 1997； 單業(yè)華等， 1998； 范桃園等， 1999； 楊起和湯達(dá)禎， 2000； 張映紅和顧家裕， 2003； 萬志軍等，2005； 劉洪林等， 2005； 楊興科等， 2005； 章邦桐等，2007， 2012； 周安朝等， 2010； 朱傳慶等， 2010； 王大勇等， 2011； 王民等， 2011； 王滿等， 2012； 齊天，2012； 趙慈平等， 2012； 唐曉音等， 2013； 馬野牧等，2013）， 只是專注于固體地球科學(xué)的巖漿巖巖石學(xué)家和礦床學(xué)家沒有給予適當(dāng)?shù)淖⒁舛选?/p>

測定巖漿熱場的方法很多， 在微觀尺度上有鏡質(zhì)體反射率、牙形石色變指數(shù)、孢粉體的半透明度、孢子體的熒光性、干酪根中碳的重量、黏土礦物的轉(zhuǎn)化序列、流體包裹體測溫、磷灰石裂變徑跡、牙形石色變指數(shù)（CAI）、自生礦物組合、氧同位素和40Ar/39Ar比值等； 大尺度上有居里面及熱泉分布等（吳江和李思田， 1993）。此外， 還有接觸變質(zhì)暈方法、蝕變分帶、礦化分帶方法等。根據(jù)煤田和石油地質(zhì)部門多年實(shí)踐的經(jīng)驗(yàn)， 鏡質(zhì)體（或鏡質(zhì)組）反射率（vitrinite reflectance）是一個非常簡便、有效、實(shí)用、精確的方法。需要提出的是， 楊瑞琰和鮑征宇（2005）對江西銀山礦床隱伏巖體的研究采用了流體動力學(xué)的熱質(zhì)輸運(yùn)-反應(yīng)體系模擬了該區(qū)的溫度場， 即包含了巖漿熱場的概念。那么， 什么是巖漿熱場呢？簡略地說， 它指的是： 在一個很短的時間內(nèi)， 在一個局部的地區(qū)出現(xiàn)的巖漿活動， 使該區(qū)域地?zé)崽荻让黠@上升， 形成一個局部區(qū)域的瞬間熱場。熱場的規(guī)模通常很小， 離巖體約幾米至幾公里。熱異常和等溫線疊加在地?zé)釄鲋希?大體垂直于地?zé)釄龇植?。巖漿熱場范圍的大小和形狀與侵入體的溫度、成分、形態(tài)、大小、侵入深度以及流體、構(gòu)造、圍巖性質(zhì)等有關(guān)。巖漿熱場與地?zé)釄霾煌?包括它們各自熱的來源不同， 熱的分布不同， 地溫梯度不同， 熱場規(guī)模不同， 持續(xù)的時間不同， 熱與流體的關(guān)系不同以及研究方法不同等（見張旗等， 2013， 2014a， 2014b）。

1　什么是鏡質(zhì)體反射率？

鏡質(zhì)體是高等植物木質(zhì)素經(jīng)過生物化學(xué)降解作用和凝膠化作用而形成的膠狀體， 在煤和炭質(zhì)泥頁巖中含量最高， 在海相碳酸鹽巖中含量最低。鏡質(zhì)體在受熱時， 化學(xué)組成和結(jié)構(gòu)會發(fā)生相應(yīng)的變化，其中最具特征的是反射率的變化。鏡質(zhì)體反射率即鏡質(zhì)體表面反射光與入射光之間的比率， 通常用油浸物鏡下測得的反射率Ro（%）來表示（蔣國豪等，2001； 王強(qiáng)， 2007）。

鏡質(zhì)組反射率是煤化程度的函數(shù)， 取決于溫度、壓力和時間等因素。平均反射率值主要受溫度和溫度作用的時間所控制， 而反射率各向異性（即雙反射率）則主要受煤化過程中應(yīng)力條件的影響（曹代勇和王文俠， 1990）。Barker （1986）專門研究過鏡質(zhì)體反射率與古地溫之間的關(guān)系， 通過世界上35個地區(qū)600多個腐殖型有機(jī)質(zhì)的平均鏡質(zhì)體反射率Ro值與對應(yīng)的最高溫度Tpeak統(tǒng)計(jì)分析得出：

該方法剔除了時間因素的影響， 尤其適用于研究存在地層剝蝕的多旋回沉積盆地（疊合盆地）， 但是該方法在Ro值小于0.9%時計(jì)算的古地溫偏低， 而在較高成熟度階段計(jì)算的古地溫偏高， 因此必然帶來計(jì)算出的古地溫梯度也會是較高的（王強(qiáng)， 2007）。

鏡質(zhì)體反射率的測定方法類似于光薄片研究方法， 需將巖石樣品制成光片， 用顯微光度計(jì)測定。由于鏡質(zhì)體反射率與溫度之間存在一定的函數(shù)關(guān)系，還由于它具有對熱的敏感性和變化的不可逆性， 且易于精確測定， 被視為靈敏的地質(zhì)溫度計(jì)（表1）， 被石油和煤田地質(zhì)部門作為最佳古地溫方法普遍加以應(yīng)用（楊起等， 1987a， 1987b； 王強(qiáng)， 2007）。鏡質(zhì)體反射率與煤質(zhì)、煤變質(zhì)的溫度及其埋深深度有密切的關(guān)系， 可用來反推其形成的溫度和埋深深度（表2）。

表1　鏡質(zhì)體最大反射率（Ro，max）與其形成溫度對應(yīng)值（楊起等， 1987b）Table 1　Maximum vitrinite reflectances （Ro，max） and their corresponding temperatures

鏡質(zhì)體反射率受影響因素很多， 包括溫度、構(gòu)造、壓力等， 不同實(shí)驗(yàn)室測定的數(shù)據(jù)不同也有影響（曹代勇和王文俠， 1990； 蔣國豪等， 2001； 馬安來和張大江， 2002）。Ro主要與溫度有關(guān)， 斷裂、構(gòu)造等也能造成局部Ro的升高， 有人甚至認(rèn)為構(gòu)造帶的影響可能比溫度更加明顯。筆者認(rèn)為， 這可能更適合Ro比較小的情況。對于尋找隱伏巖體來說， 首先， 很低的Ro可能是地?zé)嵩鰷芈室鸬模?我們只關(guān)注Ro數(shù)值高的部分（至少Ro＞2%）。其次， 對Ro的精確度要求并不高， 即使Ro誤差達(dá)到±0.5%， 也不會影響我們的判斷。我們的研究認(rèn)為， 如果Ro＞3%， 應(yīng)當(dāng)主要與巖漿熱場有關(guān)， 無煙煤分布區(qū)的下部大概就有隱伏巖體出現(xiàn)。第三， 斷裂的作用導(dǎo)致鏡質(zhì)體反射率等值線主要為拉長的條帶狀分布， 如果鏡質(zhì)體反射率等值線為渾圓狀分布， 即可能是巖體的影響而非構(gòu)造的影響。再說， 即使是斷裂的影響， 斷裂導(dǎo)致Ro數(shù)值升高也是與斷裂帶來的熱和熱水有關(guān)， 根源還是下部有巖漿。

表2　煤系烴源巖有機(jī)質(zhì)演化階段及煤化作用參數(shù)表（轉(zhuǎn)引自程喆等， 2011）Table 2　The evolution stages of organic matters in hydrocarbon source rocks and coalification parameters

如何判斷Ro的變化是否與巖漿有關(guān)并非難事。國外根據(jù)25個煤田， 時代自二疊紀(jì)至第三紀(jì)的22540個鏡質(zhì)組反射率的統(tǒng)計(jì)分析， 在深成變質(zhì)作用下， 與埋深1000 m、2000 m、3000 m、4000 m所對應(yīng)的Ro眾數(shù)值分別為0.344%、0.415%、0.530%、0.630%（楊起等， 1988）。如桂西北地區(qū)， 上二疊統(tǒng)煤層的最大埋深小于6000 m， 按照上述標(biāo)準(zhǔn)， 對應(yīng)的Ro應(yīng)小于1%。但實(shí)際上區(qū)內(nèi)的Ro值遠(yuǎn)大于1%（2%～5%）， 顯然是一種異常的情況， 表明有來自其他熱源的強(qiáng)烈影響（吳江和李思田， 1993）。這種情況可以用圖1予以說明。圖1是一個虛擬的巖漿熱場圖。在圖1a中， 地?zé)釄龅牡葴鼐€大體是水平方向分布的， 巖漿熱場等溫線圍繞巖漿侵入體， 大體是垂直分布的， 由侵入巖向外溫度逐漸降低。如果在相應(yīng)的深度得到的Ro值與全球平均值接近（例如，當(dāng)深度達(dá)到5000 m時， Ro仍然＜1%， 圖1b）， 其演化線（圖1b中的綠色粗實(shí)線）指示該區(qū)為正常的地?zé)釄?。圖1c鉆孔2的上部顯示為正常地?zé)釄龅奶卣?，下部Ro逐漸增高， 在同樣的深度上， Ro明顯大于鉆孔1得到的Ro值（＞3%， 圖1c）。下部Ro值的明顯增高解釋為巖漿熱場的影響（圖1c中的紅色粗實(shí)線）。鉆孔3開孔在巖漿熱場范圍， Ro值很高， 大多在2%～5%之間。Ro值與深度關(guān)系的曲線變緩（紅色粗實(shí)線）， 指示受巖漿熱場的控制（圖1d）。鉆孔4同鉆孔1一樣， Ro值變化屬于正常地?zé)釄龅那闆r（以綠色粗實(shí)線代表）。但在局部深度Ro出現(xiàn)異常高值（Ro＞3%，圖1e的紅色粗實(shí)線）， 過了這個深度又恢復(fù)到正常地?zé)釄龅那闆r（Ro值＜1%， 綠色粗實(shí)線）。暗示鉆孔旁側(cè)有巖床或巖株出現(xiàn)， 雖然在鉆孔中并未見到侵入體。圖1展示的是一個巖漿熱場在剖面上變化的情況， 該熱場在平面圖上通常會出現(xiàn)一個Ro高的等值線范圍。

總的說來， 目前鏡質(zhì)體反射率是古地溫研究中應(yīng)用最廣泛、研究最深入也最為有效的指標(biāo)。雖然利用它來研究古地溫也存在著一定的局限性和誤差，但總體說來這是一種費(fèi)用少、準(zhǔn)確性較高的方法（李善鵬和邱楠生， 2003）。

2　應(yīng)用鏡質(zhì)體反射率方法尋找隱伏巖體的幾個實(shí)例

應(yīng)用鏡質(zhì)體反射率方法尋找隱伏巖體已經(jīng)取得許多成果， 下面舉幾個實(shí)例。下述實(shí)例中， 古溫度測量有些采用了鏡質(zhì)體反射率方法， 有些不是。部分作者還根據(jù)古溫度場的分布對隱伏巖體進(jìn)行了預(yù)測，是很有預(yù)見性和超前性的。

2.1北京房山巖體

宋鴻林和朱寧（1998）對北京房山巖體圍巖的變形作用和熱場進(jìn)行了研究， 得出了圍繞房山巖體的古地溫分布圖（圖2a）。他們采用4種方法估計(jì)古溫度： （1）變質(zhì)礦物組合所反映的變質(zhì)相方法； （2）礦物地質(zhì)溫壓計(jì)方法（包括角閃石-斜長石、石榴石-黑云母、白云石-方解石溫度計(jì)等）； （3）方解石脈中鈣鎂含量方法； （4）包裹體測溫法。房山巖體熱場分布范圍很寬， 以溫度降低至200 ℃為界， 巖漿熱場的寬度約達(dá)30 km， 顯然不是房山巖體地表出露的部分能夠解釋的， 推測可能深部有隱伏的大巖體（圖2b）。

圖1　巖漿熱場與鏡質(zhì)體反射率的關(guān)系示意圖Fig.1　The relationship between magma thermal field and vitrinite reflectance

早先的研究大多認(rèn)為房山巖體（又稱為周口店巖體）為小巖株， 是底辟侵位的（馬昌前， 1988； 王人鏡和馬昌前， 1989； 燕濱等， 2008）。而根據(jù)巖漿熱場的分布它應(yīng)當(dāng)為一個規(guī)模較大的巖基， 地表出露的部分僅是冰山之一角。圖2a系根據(jù)宋鴻林和朱寧（1998）的資料簡化， 圖2b為筆者解釋溫度場分布試擬的， 僅供參考。

2.2湖南騎田嶺巖體

煤田地質(zhì)研究表明， 圍繞中生代的騎田嶺花崗巖煤質(zhì)存在有規(guī)律的變化： 靠近巖體為高變質(zhì)的石墨和無煙煤， 遠(yuǎn)離巖體為貧煤。圖3即是根據(jù)騎田嶺巖體周圍煤的揮發(fā)分變化和一些熱變質(zhì)礦物的分布結(jié)合X-衍射等對煤結(jié)構(gòu)的測定得出來的。圖3中紅色虛線示揮發(fā)分等值線（1%～10%）； 藍(lán)色鎖線示煤變質(zhì)帶界線； I、II、III、IV分別代表變質(zhì)帶級別的編號。自巖體向外依次可分為4個煤變質(zhì)帶（畢汝泉等， 1981）： 第I變質(zhì)帶為石墨帶和超無煙煤帶， 寬約幾百米， 煤和圍巖中出現(xiàn)大量電氣石、紅柱石等變質(zhì)礦物， 為近巖體的熱接觸變質(zhì)帶； 第II變質(zhì)帶為高-中變質(zhì)無煙煤帶， 寬度2～5 km， 溫度大約在250～410 ℃之間； 第III帶為低變質(zhì)無煙煤帶， 寬約9 km， 溫度大約在250 ℃以下； 第IV帶為貧煤帶（包括部分低變質(zhì)無煙煤）， 最寬約17 km。研究表明，騎田嶺花崗巖接觸面向外傾（畢汝泉等， 1981）， 因此導(dǎo)致熱場范圍很寬。熱場實(shí)際寬度估計(jì)不超過幾公里。

2.3山西太原西山煤田

太原西山煤田是位于沁水盆地西北側(cè)的一個石炭系-二疊系煤田， 在燕山期發(fā)生的巖漿侵入事件在該煤田形成了異常的古地?zé)釄觯?使煤級在短暫地質(zhì)時限內(nèi)迅速增高。劉洪林等（2005）對太原西山煤田的研究發(fā)現(xiàn)， 鏡質(zhì)組反射率圍繞狐偃山巖體分布， 據(jù)此還發(fā)現(xiàn)在西山煤田的東南部還存在一個更大的隱伏侵入巖體， 該侵入體已被重磁探測和煤田鉆探證實(shí)， 是一個位于石炭系-二疊系之下的燕山期石英二長巖， 命名為祁縣花崗巖（圖4）。圖4中實(shí)線為Ro等值線， 虛線為煤層露頭范圍； 紅色示狐堰山巖體，虛線圈定的范圍為祁縣隱伏巖體。

圖2　北京房山變質(zhì)核雜巖中生代早期地?zé)岙惓５葴鼐€略圖（據(jù)宋鴻林和朱寧， 1998修改）Fig.2　Isothermal map of the Early Mesozoic paleogeothermal anomaly around the Fangshan metamorphic core complex， Beijing

圖3　湖南梅田礦區(qū)龍?zhí)督M12號煤層變質(zhì)帶示意圖（據(jù)畢汝全等， 1981簡化）Fig.3　Schematic map of the metamorphic belt of No. 12 coalbed of the Longtan Formation in the Meitian mine， Hunan province

圖4　太原西山煤田巖漿活動與Ro等值線圖（據(jù)劉洪林等，2005）Fig.4　Map showing the magmatism and Rocontour in the Taiyuan Xishan coal field

2.4湖南邵陽煤田

邵陽地區(qū)中生代花崗巖分布廣， 對煤的變質(zhì)作用有很大的影響， 使該區(qū)晚二疊世煤的煤級普遍達(dá)氣煤或肥煤階段。在靠近侵入巖的部位， 煤的變質(zhì)程度提高到焦煤、瘦煤、貧煤、無煙煤甚至石墨階段。如鄧家鋪、羅城、牛頭寨-高桂山地區(qū)（圖5）， 有中生代侵入體及隱伏巖體分布， 導(dǎo)致形成了以鄧家鋪-牛頭寨-高桂山為中心的煤變質(zhì)環(huán)帶。從圖5看，邵東、祁陽地區(qū)Ro，max＜1.0%， 而在有中生代侵入巖分布的地區(qū)， Ro，max通常＞3.0%， 在鄧家鋪、牛頭寨-高桂山一帶甚至＞5.0%。在Ro，max＞3.0%和熱液礦點(diǎn)密集分布的地方， 周春光等（1996）推測可能存在3個隱伏巖體（圖5虛線圈定的棕色范圍）。

圖5　邵陽含煤區(qū)龍?zhí)督M2煤Ro，max（%）等值線圖（據(jù)周春光等， 1996）Fig.5　Isogram for Ro，max（%） of coal No.2 in the Longtan Formation in the Shaoyang coal-bearing region

2.5湖南漣源煤田

湘中漣源凹陷周邊有印支期（淺紫色）和燕山期（紅色）花崗巖分布。石炭系測水組煤系Ro在靠近巖體處最大值為10.58%， 向盆地中心降為3.00%（畢華等， 1997）， 導(dǎo)致從盆地邊緣向中心烴源巖熱演化程度遞減的現(xiàn)象。盆地東北為印支期的溈山花崗巖，但是， 并不是Ro最高的部位（圖6）。最高的Ro值在南部， 畢華等（1997）推測可能是由龍山2個隱伏的花崗巖（黃色）引起的。

3　根據(jù)Ro推測隱伏巖體的幾個實(shí)例

下面根據(jù)某些地區(qū)已經(jīng)得到的鏡質(zhì)體反射率，嘗試推測可能的隱伏巖體。

3.1西藏比如盆地

王先美等（2011）報(bào)道了青藏高原比如盆地中-上侏羅統(tǒng)烴源巖研究的成果， 他們對該區(qū)鏡質(zhì)體反射率研究的結(jié)果展示在圖7中。從圖7看， 比如盆地鏡質(zhì)體反射率普遍較高（Ro＞2.0%）， 其中Ro＞4.0%的數(shù)據(jù)分布在盆地中部和東部。該區(qū)燕山期和喜山期巖漿巖發(fā)育， 但是， 該區(qū)地表出露的巖漿巖似乎還不足以解釋該區(qū)Ro的分布， 尤其中部Ro＞4.0%的部位。因此， 筆者推測在比如和那曲之間Ro＞4.0%的部位之下可能有隱伏巖體（圖7中用白色粗虛線表示，紅色示侵入體， 數(shù)字示Ro（%）。

圖6　漣源煤盆地測水組煤變質(zhì)分帶（據(jù)畢華等， 1997）Fig.6　The coal metamorphic belt of the Ceshui Formation in the Lianyuan basin

圖7　比如盆地拉貢塘組Ro分布圖（據(jù)王先美等， 2011修改）Fig.7　Ro（%） distribution of the Lagongtang Formation in the Biru Basin

3.2松潘-阿壩地區(qū)

松潘阿壩與四川盆地僅一山（龍門山）之隔， 地質(zhì)情況卻完全不同： 四川盆地沉積地層發(fā)育， 中生代地層盛產(chǎn)油氣， 而松潘-阿壩地區(qū)中生代僅有三疊紀(jì)地層出露， 且多印支期花崗巖侵入體（金維浚等，2005）。四川盆地三疊紀(jì)為前陸盆地， 后逐漸向山間盆地演化， 侏羅紀(jì)以來則主要為河流湖沼相沉積。而松潘阿壩地區(qū)三疊紀(jì)后整體抬升， 侏羅系、白堊系及第三系僅零星分布（范明等， 2006）。因此， 四川盆地和松潘-阿壩地區(qū)的Ro值也明顯不同， 松潘阿壩地區(qū)Ro明顯地高， 普遍＞2%， 局部＞4%， 主要受巖漿活動的影響（范明等， 2006）。從圖8看， Ro＞4%的范圍主要分布在圖的南部和西北部（松潘-紅原-阿壩以南、久治以北和瑪曲以西）。據(jù)筆者初步考察， 該區(qū)三疊紀(jì)地層出露廣泛， 厚度巨大， 褶皺強(qiáng)烈， 部分地區(qū)經(jīng)歷了明顯的低級-甚低級變質(zhì)作用，為板巖和千枚巖。上述地區(qū)三疊紀(jì)巖漿活動的規(guī)模很小， 僅紅原之南的達(dá)古冰川花崗巖規(guī)模較大， 不足以解釋該區(qū)Ro明顯普遍高的現(xiàn)象。筆者推測， 該區(qū)Ro值普遍較高（＞2%）可能與區(qū)域強(qiáng)烈構(gòu)造變動導(dǎo)致的區(qū)域熱變質(zhì)作用有關(guān)， 而該區(qū)幾個Ro特別高的地區(qū)（Ro＞4%～8%）， 可能與巖漿侵入活動有關(guān)。但目前地表出露的零星小巖體無法解釋Ro值的分布狀況， 估計(jì)深部可能有較大的隱伏巖體存在（圖8虛線圈定的黃色部分）。松潘以南Ro明顯呈圓形分布， Ro值最高， 類似變質(zhì)核雜巖的熱分布狀態(tài)， 顯然可能是巖漿侵入造成的。毛爾蓋之南的Ro高值與達(dá)古冰川花崗巖出露的位置大體吻合， 根據(jù)Ro的分布， 巖體深部可能規(guī)模更大。而紅原之南的Ro高值可能暗示在馬爾康地區(qū)（已出圖）應(yīng)該有一個隱伏的大巖體。此外， 在久治以北和瑪曲以西也可能有一個隱伏巖體。

圖9　南華北地區(qū)上古生界山西組Ro等值線（據(jù)趙俊峰等， 2010， 2011簡化、修改）Fig.9　Contour map of the Rovalues from the Shanxi Formation of the Upper Paleozoic in the southern North China Craton

圖8　松潘-阿壩地區(qū)Ro分布圖（據(jù)范明等， 2006修改）Fig.8　Rodistribution of the Songpan-Aba area

3.3南華北盆地

煤田和石油地質(zhì)部門對南華北盆地煤田資源進(jìn)行了許多研究， 他們得出了一份很有價值的鏡質(zhì)體反射率分布圖（圖9）。該圖據(jù)趙俊峰等（2010， 2011）的資料修改。Ro等值線僅選用了Ro＞2.0%的數(shù)據(jù)，Ro＜2.0%的等值線沒有予以表示。筆者推測該區(qū)Ro＞3.0%的資料可能與巖漿熱場有關(guān)， 尤其濟(jì)源地區(qū)。圖中淮北地區(qū)的資料據(jù)趙俊峰等（2010， 2011）、李雷（2011）和汪青松（2010）。

據(jù)研究， 南華北盆地沉積巖厚度很大， 如果按照正常的地?zé)嵩鰷芈剩?深度達(dá)到3000～4000 m， Ro大約在1.0%左右， 而該區(qū)不少地方Ro超過3%， 甚至達(dá)到6%， 這顯然不是正常的情況。許多人認(rèn)為， 南華北盆地異常高的Ro可能是區(qū)域變質(zhì)作用和局部巖漿作用兩個因素導(dǎo)致的結(jié)果。但是， 喆程等（2011）認(rèn)為， 該區(qū)中生代以來的巖漿作用是造成局部地區(qū)熱演化程度較高的根本原因。筆者贊同這一見解，認(rèn)為導(dǎo)致南華北Ro高的原因可能主要是中生代巖漿活動造成的， 可能的依據(jù)如下：

（1） 石油和煤田地質(zhì)勘探在若干地區(qū)已經(jīng)鉆遇了中生代巖漿巖， 包括中酸性火山巖和侵入巖。趙俊峰等（2011）指出， 南華北已發(fā)現(xiàn)的巖漿巖主要來自燕山期和喜山期兩個時代， 本文只關(guān)注燕山期的巖漿活動。周口坳陷燕山期巖漿活動主要分布在周口中南部凹陷帶， 在沈丘凹陷中部周參10井、周參11井已經(jīng)鉆遇燕山期中酸性火山巖及花崗巖， 汝南凹陷周參17井見1700 m安山巖， 同位素年齡為135 Ma。太康隆起太參2井、新太參1井也鉆遇了燕山期花崗巖及中酸性噴出巖（趙俊峰等， 2010， 2011， 圖9）。

（2） 據(jù)趙俊峰等（2010）報(bào)道， 南華北盆地北部的濟(jì)源、滎鞏及焦作等煤田（圖9的西北角）大多經(jīng)歷了異常高、但又不均勻的差異熱演化過程。如濟(jì)源地區(qū)的克井礦區(qū)， 山西組下部煤層Ro值高達(dá)6.0%～7.0%， 已為超級無煙煤； 向西不到20 km的下冶礦區(qū)Ro減小為1.89%～2.5%， 達(dá)貧、瘦及無煙煤階段。偃龍煤田10個地點(diǎn)的反射率都在5.0%以上， 到諸葛礦Ro又降到2.91%。由東向西逐漸變低的煤變質(zhì)趨勢十分明顯。滎鞏煤田Ro在4.2%～5.7%之間。趙俊峰等（2010）指出， 上述煤田雖熱演化程度很高，但少見巖漿侵入煤系的直接證據(jù)， 故不認(rèn)為是巖漿活動造成的。但是， 筆者認(rèn)為， 上述地區(qū)短距離內(nèi)Ro的明顯變化， 極有可能是巖漿侵入的結(jié)果。

（3） 據(jù)趙俊峰等（2011）報(bào)道， 在南華北及周邊某些煤田中， 發(fā)育高變質(zhì)含煤區(qū)， 在這些煤田內(nèi)及其周邊已發(fā)現(xiàn)多處燕山期巖漿活動侵入煤系的現(xiàn)象。如永城-商丘一帶有閃長巖和輝綠巖以巖株、巖墻或巖脈形式侵入煤層。淮北朔里礦高變質(zhì)煤與北鄰丁里巖體有關(guān)。陜澠礦區(qū)的支建煤礦、平頂山石龍礦區(qū)、確山吳桂橋等地也發(fā)現(xiàn)了巖漿侵入煤系的現(xiàn)象（吳桂橋英安巖的U-Pb鋯石年齡為123 Ma）。據(jù)汪青松（2010）報(bào)道， 在徐州、永城、淮北宿遷地區(qū)， 燕山期巖漿活動以閃長玢巖、石英閃長玢巖、石英二長閃長玢巖為主， 分布面積已達(dá)數(shù)百余平方千米， 與鐵、銅、金等礦產(chǎn)關(guān)系密切。

（4） 中國東部巖漿活動研究表明， 華北在中生代存在大規(guī)模巖漿活動。只要是有基巖出露的地區(qū)，中生代巖漿巖總是占有相當(dāng)?shù)臄?shù)量， 如冀北、燕山、膠東、遼東、豫西等地的實(shí)例。因此， 圖9從濟(jì)源經(jīng)鄭州、開封到太康的區(qū)域內(nèi)， 鏡質(zhì)體反射率（Ro）大于3%的部分， 很可能深部有較大的隱伏巖體。

4　尋找隱伏巖體的意義

尋找隱伏巖體對找礦是很有意義的， 尤其在礦區(qū)范圍內(nèi)外， 如果含礦熱液發(fā)育， 且與巖漿活動有關(guān)， 這時， 尋找隱伏巖體即是深部找礦的重要手段。下面舉幾個實(shí)例：

4.1滇黔桂金三角卡林型金礦利用熱作用參數(shù)找礦的嘗試

吳江和李思田（1993）研究了滇黔桂“金三角”地區(qū)熱作用參數(shù)（熱作用參數(shù)即巖漿熱場， 筆者注）與卡林型金礦的關(guān)系。他們的研究表明， 熱作用參數(shù)既是古地?zé)峄顒恿粝碌暮圹E， 也是恢復(fù)古地?zé)峄顒拥墓ぞ?。?yīng)用熱作用參數(shù)一是建立古地?zé)釄觯?通過成礦期古地?zé)釄龅慕ⅲ?結(jié)合成礦溫度， 推斷成礦作用在空間上的分布情況； 二是判斷隱伏巖體， 揭示出那些貌似與巖漿活動無關(guān)的成礦區(qū)的真實(shí)情況，并推斷可能由巖漿熱液活動所帶來的礦質(zhì)。

4.2美國內(nèi)華達(dá)州卡林型金礦研究的啟示

美國內(nèi)華達(dá)州是卡林型金礦的發(fā)源地， 那里的許多卡林型金礦區(qū)除常見脈巖外， 還有侵入巖和火山巖的發(fā)育（Ressel and Henry， 2006）。研究表明， 該區(qū)卡林型金礦與始新世（42～36 Ma）火山巖和淺成侵入巖有關(guān)。有些礦區(qū)無淺成侵入巖出露， 但是， 根據(jù)巖脈的分布以及航磁異常推測， 深部可能存在隱伏巖體（圖10a）。圖10a是據(jù)Ressel and Henry （2006）的資料修改的， 不同顏色的虛線代表不同溫度的等溫線， 系筆者試擬， 只表示等溫線分布的變化， 沒有定量的概念。圖10b也是筆者試擬的， 代表可能的鏡質(zhì)體反射率Ro變化的趨勢。筆者認(rèn)為， 如果該剖面根據(jù)航磁資料推測的隱伏巖體的深度是可能的，則可形成一個圍繞巖體分布的熱場。由于隱伏巖體南部淺北部深（圖10a）， 故巖漿熱場的溫度是南高北低。如果沿剖面做Ro測量， 其值也應(yīng)當(dāng)是南高北低（圖10b）。如果上述情況是合理的， 即可推測在北部的深部還有找到卡林型金礦的可能。

4.3甘肅南部卡林型金礦擴(kuò)大找礦的可能性

甘肅南部位于陜甘川金三角的西北部， 筆者曾經(jīng)對該區(qū)印支期巖漿巖進(jìn)行過研究（金維浚等， 2005；張旗等， 2009）， 發(fā)現(xiàn)該區(qū)北部巖漿巖分布較多， 而南部較少。北部成礦有卡林型、矽卡巖型和斑巖型等， 南部主要是卡林型（如陽山、大水、馬腦殼、拉爾瑪、鹿兒壩、早子溝等）。鑒于西秦嶺印支期金銅尤其是金礦分布多， 儲量大， 可能是中國金礦最具潛力的地區(qū)之一。我們設(shè)想： 卡林型金礦出露的地區(qū)雖然很少有侵入體的顯示， 或僅有一些小巖脈、巖株、巖床， 但是， 深部可能隱伏有較大的巖體?？中徒鸬V最可能出現(xiàn)在距大巖體一定距離的一定的巖漿熱場范圍內(nèi)， 如圖10a所示。為此， 可以嘗試從查明隱伏巖體入手： 首先利用鏡質(zhì)體反射率和其他方法查明該區(qū)巖漿熱場的分布， 識別可能的隱伏巖體的信息（如Ro高值區(qū)）； 其次在若干Ro值高的區(qū)域進(jìn)行物化探研究， 查明是否有成礦作用的顯示； 第三， 在上述研究的基礎(chǔ)上選擇若干有成礦潛力的區(qū)域開展鉆探研究。該區(qū)范圍很大， 長寬約500×300 km2， 似乎無從下手。鑒于該區(qū)已經(jīng)積累了很多資料，可以先從航磁、重力、化探、遙感等方法入手， 選擇可能有隱伏巖體顯示的區(qū)域， 開展小比例尺鏡質(zhì)體反射率測量。然后， 逐漸縮小范圍， 進(jìn)行大比例尺鏡質(zhì)體反射率測量， 確定進(jìn)一步研究的區(qū)域， 實(shí)現(xiàn)找礦的突破。

圖10　美國內(nèi)華達(dá)州卡林礦帶北-中部剖面示意圖（Ressel and Henry， 2006）Fig.10　Schematic cross section of the northern and central Carlin trend in Nevada， America

4.4銅陵礦集區(qū)深部找礦問題

銅陵礦區(qū)深部巖體及找礦問題已經(jīng)有了許多的研究， 例如嚴(yán)加永等（2009）利用航磁3D反演成像技術(shù)識別隱伏巖體， 取得了很好的效果。采用該方法可以獲得巖體的3D形態(tài)， 如頂部埋深、界面起伏、空間形態(tài)等（嚴(yán)加永等， 2009）。此外， 許多人認(rèn)為銅陵深部可能存在一個大巖體， 銅陵的銅和金， 包括許多小巖體， 即來源于這個深部的未知大巖體。王慶飛等（2004）的研究支持這種見解。為此， 建議采用鏡質(zhì)體反射率方法對上述研究進(jìn)行驗(yàn)證： （1）對地表早于侵入巖時代的地層進(jìn)行采樣（主要采泥質(zhì)巖、含炭質(zhì)的泥質(zhì)巖）， 圈定Ro等值線， 查明深部是否有隱伏巖體出現(xiàn)的可能。（2）沿鉆孔取樣， 考察鉆孔中是否有Ro突變的情況（如圖1e鉆孔4的情況）， 如果有Ro突變， 則暗示鉆孔旁可能有隱伏巖體。當(dāng)然， 一孔之見不足以確定隱伏巖體的位置， 應(yīng)當(dāng)利用多個鉆孔， 編制Ro等值線三維圖予以確定。（3）選擇鉆進(jìn)深度大的若干鉆孔， 查明鉆孔底部Ro值的變化情況，如果該值超過正常深度Ro的允許值且沒有降低的趨勢（如Ro＞3%）， 則暗示在終孔位置以下的深部有隱伏巖體存在的可能性。銅陵礦區(qū)碳酸鹽地層出露廣泛， 是一個不利于該方法使用的因素， 但是， 仍然有一些含泥質(zhì)巖的地層， 尤其是二疊系的龍?zhí)督M，泥頁巖豐富， 夾薄煤層， 是最佳測定對象。

4.5與巖漿作用有關(guān)的熱液礦床、矽卡巖型礦床等的找礦問題

這類礦床很多， 如金銅鎢錫鉛鋅鉬銀汞銻等。

首先， 對一些老礦區(qū)、研究程度較高的礦區(qū)， 可以考慮利用鏡質(zhì)體反射率方法從尋找隱伏巖體入手去擴(kuò)大找礦。老礦區(qū)資料多、鉆孔多， 有條件開展鏡質(zhì)體反射率的三維立體研究。

其次， 對一些新礦區(qū)開展鏡質(zhì)體反射率測量也是有益的， 雖然新區(qū)鉆孔少是一個不足， 但是， 對于了解鉆孔中沒有見到而鉆孔旁側(cè)可能出現(xiàn)的隱伏巖體可能是有啟示的。因此， 對新區(qū)地表和鉆孔兩方面的鏡質(zhì)體反射率研究都不應(yīng)忽視。

第三， 礦區(qū)隱伏巖體有一個特點(diǎn)， 即巖體與流體（熱液蝕變）活動可能有關(guān)。因此， 在尋找隱伏巖體的同時不可忽略對流體的追蹤。找隱伏巖體不是我們的目的， 目的是找礦， 找礦離不開流體。如果一旦發(fā)現(xiàn)了隱伏巖體的跡象， 應(yīng)立即開展地質(zhì)及相應(yīng)的物化探研究， 查明是否有含礦流體出現(xiàn)及其與隱伏巖體的關(guān)系。

第四， 礦區(qū)如果小巖體多， 礦體多， 構(gòu)造必定多。因?yàn)椋?巖體和礦體可能主要受構(gòu)造制約， 沿一定的構(gòu)造發(fā)育。這時， 不僅巖體是熱的（形成巖漿熱場），流體是熱的（形成流體熱場）， 而且構(gòu)造部位也可能是熱的（形成構(gòu)造熱場）。因此， 一旦發(fā)現(xiàn)鏡質(zhì)體反射率出現(xiàn)異常的高值， 即可著手找礦方面的研究。

5　結(jié)　語

（1） 尋找隱伏巖體是有意義的： 一是對了解大片覆蓋區(qū)巖漿活動期間的地殼厚度和地殼熱狀態(tài)有用（張旗， 2014）； 二是對熱液礦床找礦有用。因此，伴有礦化的隱伏巖體是值得特別關(guān)注的， 尋找隱伏巖體， 找隱伏的含礦巖體， 是找礦的重要目標(biāo)。

（2） 尋找隱伏巖體的方法很多， 如常用的磁、重力、電法等， 但是， 忽略了一個熱法。相比之下， 利用巖漿熱場理論， 采用鏡質(zhì)體反射率方法尋找隱伏巖體， 可能是一個最簡單、精確、實(shí)用的方法：

簡單是因?yàn)樗?jīng)濟(jì)省錢， 設(shè)備簡單， 只需要一臺光學(xué)顯微鏡即可。我們甚至可以不用添加設(shè)備，不用建立實(shí)驗(yàn)室， 野外采樣直接送煤田或石油部門實(shí)驗(yàn)室進(jìn)行測試即可。

精確是由于鏡質(zhì)體屬于有機(jī)質(zhì)， 對溫度十分敏感， 可以感知十幾度范圍內(nèi)溫度的變化； 因此， Ro的靈敏度很高， 數(shù)據(jù)精確可靠（與地質(zhì)方法、變質(zhì)暈方法比較）。盡管鏡質(zhì)體反射率具有非均質(zhì)性， 但是，如果取多次測量的平均值可避免這種缺陷。鏡質(zhì)體反射率一個最大的特點(diǎn)是其具有的不可逆性， 盡管我們所研究的對象巖漿巖早已冷卻了， 但是， Ro所記錄的卻是該樣品所達(dá)到的最大溫度值。

實(shí)用是該方法應(yīng)用領(lǐng)域廣， 只要是志留紀(jì)及其以后的含泥質(zhì)比較豐富的地層即可滿足測定Ro的需要。鏡質(zhì)體反射率主要與溫度有關(guān)， 如果Ro值超過了由地?zé)釄鰧?dǎo)致的溫度變化， 則可判斷是由巖漿巖引起的。盡管熱的流體、構(gòu)造帶也可以導(dǎo)致溫度升高， 形成熱場， 但其規(guī)模很小， 形態(tài)變化大， 易與巖漿熱場區(qū)分開。說到底， 熱流、斷裂所攜帶的熱的來源仍然是巖漿， 只是巖漿熱場的另外一種表現(xiàn)而已。有人用莫霍面抬升來解釋局部性的熱場， 在地殼淺部， 這種解釋的可能性是值得商榷的。

（3） 利用鏡質(zhì)體反射率方法找礦不受比例尺的限制， 可以大到宏觀的區(qū)域上的規(guī)模， 如全國、一省、數(shù)省的范圍； 可以是一個礦集區(qū)的范圍（如銅陵、德興、騎田嶺）； 也可以小到一個礦體、一個鉆孔、一個坑道的規(guī)模。如楊起等（1987b）對全國煤質(zhì)分布與中生代巖漿活動關(guān)系的研究， 如馬東升（1999）對大華南礦床分布的研究， 如吳江和李思田（1993）對滇黔桂金礦的研究等。

（4） 該方法既適合老礦區(qū)也適合新礦區(qū)。老礦區(qū)資料多、鉆孔多， 使用該方法的空間大， 范圍廣， 有條件開展三維立體研究； 對于新礦區(qū)來說， 鉆孔少是一個嚴(yán)重的限制， 應(yīng)注意了解鉆孔鉆進(jìn)過程中Ro的變化。尤其注意鉆孔中某些深度Ro突變的情況以及隨鉆進(jìn)深度增加Ro增加的情況， 以便及時指導(dǎo)找礦。前者暗示在Ro突變處附近有巖體的顯示（如圖1e）； 后者暗示深部可能有巖體存在（圖1c和1d）。

（5） 礦區(qū)和非礦區(qū)隱伏巖體的意義是不一樣的。在礦區(qū)， 無論新老礦區(qū)， 巖體有與流體相伴的可能性， 當(dāng)然不是指每個巖體。因此， 如果在礦區(qū)發(fā)現(xiàn)Ro出現(xiàn)異常， 應(yīng)相應(yīng)地開展礦化方面的研究， 考察隱伏巖體與成礦的關(guān)系。

（6） 由于鏡質(zhì)體是一種有機(jī)物殘留體， 只出現(xiàn)在志留紀(jì)及其以后的沉積巖中。它在煤和泥質(zhì)巖中的含量很高， 在碳酸鹽巖中含量很低。因此， 該方法的應(yīng)用也受一定的限制， 它不適合老于志留紀(jì)的地層， 不適合碳酸鹽巖地層。幸好我們隨后的研究發(fā)現(xiàn)還有一類海相鏡質(zhì)體， 它在早古生代的泥質(zhì)巖和碳酸鹽巖中保存較多， 恰好可以彌補(bǔ)鏡質(zhì)體反射率方法的不足（詳見陳萬峰等， 待刊）。

（7） 地球是一個大熱場， 熱場溫度高不一定都是巖漿引起的。如果熱場溫度很低（例如＜100 ℃），變化很小， 則很難區(qū)分是地?zé)釄鲞€是巖漿熱場。巖漿熱場是局部性的， 如果在一個小范圍內(nèi)， 溫度的變化明顯， 等溫線是渾圓狀或橢圓狀的， 則等溫線環(huán)繞的中心即可能是隱伏巖體的表現(xiàn)。鏡質(zhì)體反射率與溫度有關(guān)， 溫度低， 鏡質(zhì)體反射率Ro即低。因此， 對尋找隱伏巖體來說， 我們主要關(guān)注Ro數(shù)值高的數(shù)據(jù)。在中國， 如果Ro＞3%， 基本上可以確定其下有隱伏巖體出現(xiàn)； Ro＜1%的數(shù)據(jù)基本上是沒有意義的； 對于Ro在2%～3%的數(shù)據(jù)應(yīng)當(dāng)具體分析， 在某些情況下， 它可能是隱伏巖體的反映（如四川盆地上三疊統(tǒng)須家河組從上到下鏡質(zhì)體反射率的變化遠(yuǎn)遠(yuǎn)超出了地?zé)嵩鰷芈首兓那闆r）， 在某些情況下， 可能與隱伏巖體無關(guān)。

鏡質(zhì)體反射率是煤田和石油地質(zhì)部門非常熟悉、幾乎須臾不可離開的研究方法。幾十年來， 已經(jīng)積累了數(shù)以幾十萬計(jì)、甚至上百萬計(jì)的資料； 所涉及的地區(qū)既有廣大的無巖漿巖出露的中新生代沉積盆地區(qū)， 也有巖漿活動非常發(fā)育的地區(qū)， 其中有些對于礦床學(xué)家來說可能是有用的， 可以予以借鑒。但是， 重要的仍然是將該方法引入到礦床學(xué)研究中來， 結(jié)合熱液金屬礦床成礦的特征， 開拓出新的應(yīng)用前景， 推進(jìn)深部找礦工作更上一層樓。

后記： 鏡質(zhì)體反射率方法是筆者學(xué)習(xí)煤田和石油部門研究成果時了解的， 這是一個好方法， 對于尋找隱伏巖體非常有益， 可以嘗試應(yīng)用到巖石學(xué)和礦床學(xué)的研究中來。當(dāng)然， 應(yīng)用時需要做一些調(diào)整和改變。感謝中國地質(zhì)大學(xué)（武漢）馬昌前教授和另外兩位匿名審稿人對本文的評論與批評， 特致謝忱。

（References）：

畢華， 彭格林， 趙志忠， 鐘建華. 1997. 漣源煤盆地測水組煤的區(qū)域巖漿熱變質(zhì)作用. 煤炭學(xué)報(bào)， 22（4）：349- 354.

畢汝全， 曲星武， 王金城. 1981. 對梅田礦區(qū)煤變質(zhì)作用的探討. 煤炭科學(xué)技術(shù)， （3）： 22- 55.

曹代勇， 王文俠. 1990. 鏡質(zhì)組反射率各向異性分析技術(shù)及其在構(gòu)造研究中的應(yīng)用. 中國煤田地質(zhì)， 2（1）： 1- 8.

陳榮書， 何生， 王青玲， 蘭廷澤. 1989. 巖漿活動對有機(jī)質(zhì)成熟作用的影響初探——以冀中葛漁城-文安地區(qū)為例. 石油勘探與開發(fā)， （1）： 29- 38.

陳儒慶， 袁奎榮. 1990. 巖漿作用與煤變質(zhì)——煤的巖漿熱液接觸變質(zhì)作用. 桂林冶金地質(zhì)學(xué)院學(xué)報(bào)， 10（2）：191- 200.

陳萬峰, 張旗, 王金榮, 康文凱, 焦守濤. 海相鏡質(zhì)體:一個利用鏡質(zhì)體反射率尋找隱伏巖體的補(bǔ)充方法.大地構(gòu)造與成礦學(xué), 待刊.

程喆， 徐旭輝， 王榮新， 武明輝， 辛藜莉. 2011. 南華北地區(qū)上古生界烴源巖異常熱演化因素探討. 石油實(shí)驗(yàn)地質(zhì)， 33（2）： 142- 147.

范明， 秦建中， 張渠. 2006. 松潘阿壩地區(qū)烴源巖有機(jī)質(zhì)熱演化特征. 沉積學(xué)報(bào)， 24（3）： 440- 445.

范桃園， 石耀霖， 周炎如. 1999. 沉積盆地?zé)嵫莼^程中的巖漿作用. 中國科學(xué)院研究生院學(xué)報(bào)， 16（1）： 63- 69.

馮喬， 湯錫元. 1997. 巖漿活動對油氣藏形成條件的影響.地質(zhì)科技情報(bào)， 16（4）： 59- 65.

蔣國豪， 胡瑞忠， 方維萱. 2001. 鏡質(zhì)體反射率（Ro）推算古地溫研究進(jìn)展. 地質(zhì)地球化學(xué)， 29（1）： 40- 45.

金維浚， 張旗， 何登發(fā)， 賈秀琴. 2005. 西秦嶺埃達(dá)克巖的SHRIMP定年及其構(gòu)造意義. 巖石學(xué)報(bào)， 21（3）： 959- 966.

李雷. 2011. 巖漿侵入作用影響下北辰礦瓦斯賦存規(guī)律研究及區(qū)域預(yù)測. 焦作： 河南理工大學(xué)碩士學(xué)位論文.

李善鵬， 邱楠生. 2003. 應(yīng)用鏡質(zhì)體反射率方法研究東營凹陷古地溫. 西安石油學(xué)院學(xué)報(bào)（自然科學(xué)版）， 18（6）：9- 11.

劉洪林， 王紅巖， 趙國良， 李貴中， 楊帆， 劉洪建. 2005.燕山期構(gòu)造熱事件對太原西山煤層氣高產(chǎn)富集影響.天然氣工業(yè)， 25（1）： 28- 32.

羅文積， 陳家清. 1997. 雙向匯聚熱液成礦. 甘肅地質(zhì)學(xué)報(bào)， 6（增刊）： 44- 49.

馬安來， 張大江. 2002. 壓力對鏡質(zhì)組反射率與烴類生成的影響. 地質(zhì)地球化學(xué)， 30（1）： 85- 90.

馬昌前. 1988. 北京周口店巖株侵位和成分分帶的巖漿動力學(xué)機(jī)理. 地質(zhì)學(xué)報(bào)， （4）： 329- 341.

馬昌前， 楊坤光， 唐仲華等. 1994. 花崗巖類巖漿動力學(xué)——理論方法及鄂東花崗巖類例析. 武漢： 中國地質(zhì)大學(xué)出版社.

馬東升. 1999. 華南中、低溫成礦帶元素組合和流體性質(zhì)的區(qū)域分布規(guī)律——兼論華南燕山期熱液礦床的巨型分帶現(xiàn)象和大規(guī)模成礦作用. 礦床地質(zhì)， 18（4）：347- 358.

馬野牧， 陸現(xiàn)彩， 張雪芬， 李曉昭， 胡文瑄， 汪愷. 2013.花崗巖侵入體-泥質(zhì)圍巖熱傳輸過程的數(shù)值模擬及其地質(zhì)意義： 以粵東典型接觸帶剖面為例. 高校地質(zhì)學(xué)報(bào)， 19（2）： 307- 315.

齊天. 2012. 熱流模型在地質(zhì)資源勘探上的應(yīng)用研究.大連： 大連理工大學(xué)碩士學(xué)位論文.

單業(yè)華， 袁鼎， 李志安. 1998. 巖漿熱變質(zhì)作用數(shù)值模擬.中國煤田地質(zhì)， 10（2）： 19- 26.

宋鴻林， 朱寧. 1998. 北京西山南部中生代早期的構(gòu)造變形相和古地?zé)岙惓? 現(xiàn)代地質(zhì)， 12（3）： 302- 310.

唐曉音， 張功成， 梁建設(shè)， 楊樹春， 饒松， 胡圣標(biāo). 2013.瓊東南盆地長昌凹陷火成巖侵入體對溫度場及烴源巖成熟度的影響. 地球物理學(xué)報(bào)， 56（1）： 159- 169.

萬志軍， 趙陽升， 康建榮. 2005. 高溫巖體地?zé)豳Y源模擬與預(yù)測方法. 巖石力學(xué)與工程學(xué)報(bào)， 24（6）： 945- 949.

汪青松. 2010. 淮北地區(qū)矽卡巖型鐵銅礦床控礦條件分析與成礦模式. 資源調(diào)查與環(huán)境， 31（2）： 103- 111.

王大勇， 陸現(xiàn)彩， 徐士進(jìn)， 胡文瑄， 齊天. 2011. 沉積盆地內(nèi)侵入巖席對富含有機(jī)質(zhì)圍巖熱影響的熱傳輸模型研究. 南京大學(xué)學(xué)報(bào)（自然科學(xué)版）， 47（1）： 45- 50.

王滿， 王英偉， 薛林福. 2012. 巖漿侵入對圍巖中有機(jī)質(zhì)成熟度影響的有限元模擬. 斷塊油氣田， 19（2）： 172- 176.

王民， 盧雙舫， 劉大為， 劉楊， 武靜. 2011. 巖漿侵入體熱傳導(dǎo)模型優(yōu)選及應(yīng)用. 吉林大學(xué)學(xué)報(bào)（地球科學(xué)版），41（1）： 71- 78.

王強(qiáng). 2007. 黏土礦物地溫計(jì)研究與應(yīng)用. 青島： 中國石油大學(xué)碩士學(xué)位論文.

王慶飛， 鄧軍， 黃定華， 張強(qiáng). 2004. 銅陵礦集區(qū)淺層隱伏巖體預(yù)測及其形態(tài)分析. 礦床地質(zhì)， 23（3）： 405- 409.

王人鏡， 馬昌前. 1989. 北京周口店侵入體特征及其侵位機(jī)制. 地球科學(xué)——中國地質(zhì)大學(xué)學(xué)報(bào)， 14（4）： 399- 406.

王維， 湯靜如. 2013. 隱伏巖體預(yù)測的理論方法研究進(jìn)展.中國科技信息， （6）： 36- 38.

王先美， 伍新和， 鄭和榮， 王毅. 2011. 青藏高原比如盆地中-上侏羅統(tǒng)烴源巖. 石油學(xué)報(bào)， 32（1）： 41- 48.

吳傳榮. 1992. 煤的巖漿熱變質(zhì)作用. 煤田地質(zhì)與勘探，20（4）： 24- 29.

吳江， 李思田. 1993. 熱作用參數(shù)在成金作用分析中的應(yīng)用——以滇黔桂“金三角區(qū)”為例. 地球科學(xué)——中國地質(zhì)大學(xué)學(xué)報(bào)， 18（6）： 725- 734.

嚴(yán)加永， 呂慶田， 孟貴祥， 趙金花. 2009. 銅陵礦集區(qū)中酸性巖體航磁3D成像及對深部找礦方向的指示. 礦床地質(zhì)， 28（6）： 838- 849.

燕濱， 何斌， 徐義剛， Parterson S， 張濤. 2008. 北京西山房山巖體西北部強(qiáng)變形帶的成因新解. 大地構(gòu)造與成礦學(xué)， 32（4）： 521- 527.

楊起. 1989. 中國煤變質(zhì)研究. 地球科學(xué)——中國地質(zhì)大學(xué)學(xué)報(bào)， 14（4）： 341- 345.

楊起. 1999. 中國煤的疊加變質(zhì)作用. 地學(xué)前緣， 6（增刊）：1- 8.

楊起等. 1987a. 煤的光性研究， 煤化作用 // 楊起等. 煤地質(zhì)學(xué)進(jìn)展. 北京： 科學(xué)出版社.

楊起等. 1988. 華北石炭、二疊紀(jì)煤變質(zhì)特征與地質(zhì)因素探討. 北京： 地質(zhì)出版社.

楊起， 潘治貴， 翁成敏， 蘇玉春， 汪正平. 1987b. 區(qū)域巖漿熱變質(zhì)作用及其對我國煤質(zhì)的影響. 現(xiàn)代地質(zhì)，1（1）： 123- 130.

楊起， 任德貽. 1981. 中國煤變質(zhì)問題的探討. 煤田地質(zhì)與勘探， （2）： 1- 10.

楊起， 湯達(dá)禎. 2000. 華北煤變質(zhì)作用對煤層含氣量和滲透率的影響. 地球科學(xué)——中國地質(zhì)大學(xué)學(xué)報(bào)， 25（3）：273- 278.

楊起， 吳沖龍， 湯達(dá)禎， 康西棟， 劉大錳. 1996. 中國煤變質(zhì)作用. 地球科學(xué)——中國地質(zhì)大學(xué)學(xué)報(bào)， 21（3）：311- 319.

楊瑞琰， 鮑征宇. 2005. 應(yīng)用熱質(zhì)輸運(yùn)-反應(yīng)體系對銀山隱伏巖體預(yù)測. 地質(zhì)找礦論叢， 20（2）： 87- 92.

楊文寬. 1982. 巖漿余熱對地溫和有機(jī)質(zhì)影響的定量估計(jì).石油實(shí)驗(yàn)地質(zhì)， 4（3）： 191- 199.

楊興科， 劉池洋， 楊永恒， 季麗丹， 趙亮， 魏振權(quán)， 徐曉尹. 2005. 熱力構(gòu)造的概念分類特征及其研究進(jìn)展.地學(xué)前緣， 12（4）： 385- 396.

張健， 石耀霖. 1997. 沉積盆地巖漿侵入的熱模擬. 地球物理學(xué)進(jìn)展， 12（3）： 55- 64.

張旗. 2014. 大陸花崗巖的地球動力學(xué)意義. 巖石礦物學(xué)雜志， 33（4）： 785- 795.

張旗， 金惟俊， 李承東， 焦守濤. 2013. 地?zé)釄鲋小皫r漿熱場”的識別及其意義. 地球物理學(xué)進(jìn)展， 28（5）：2495- 2507.

張旗， 金惟俊， 李承東， 焦守濤. 2014a. 巖漿熱場： 它的基本特征及其與地?zé)釄龅膮^(qū)別. 巖石學(xué)報(bào)， 30（2）：341- 349.

張旗， 金惟俊， 李承東， 焦守濤. 2014b. “巖漿熱場”說及其成礦意義（上） . 甘肅地質(zhì)， 23（1）： 1- 18.

張旗， 殷先明， 殷勇， 金維浚， 王元龍， 趙彥慶. 2009. 西秦嶺與埃達(dá)克巖和喜馬拉雅型花崗巖有關(guān)的金銅成礦及找礦問題. 巖石學(xué)報(bào)， 25（12）： 3103- 3122.

張映紅， 顧家裕. 2003. 熱液環(huán)流： 侵入巖-外變質(zhì)帶儲層發(fā)育的重要影響因素. 特種油氣藏， 10（1）： 86- 89.

張正階， 宋謝炎， 王玉蘭， 馬潤則. 1996. 席狀巖漿房中的巖漿冷卻及動力穩(wěn)定性. 巖石學(xué)報(bào)， 12（1）： 1- 16.

章邦桐， 吳俊奇， 凌洪飛， 陳培榮. 2007. U-Th-K放射成因熱對花崗巖冷卻-結(jié)晶過程影響的計(jì)算及地質(zhì)意義.中國科學(xué)（D輯）， 37（2）： 155- 159.

章邦桐， 吳俊奇， 凌洪飛， 陳培榮. 2012. 金雞嶺產(chǎn)鈾花崗巖體印支期侵位的巖漿動力學(xué)證據(jù)及其構(gòu)造意義.鈾礦地質(zhì)， 28（1）： 11- 20.

趙慈平， 冉華， 王云. 2012. 騰沖火山區(qū)的現(xiàn)代幔源氦釋放：構(gòu)造和巖漿活動意義. 巖石學(xué)報(bào)， 28（4）： 1189- 1204.

趙俊峰， 劉池洋， 何爭光， 劉永濤. 2010. 南華北地區(qū)主要層系熱演化特征及其油氣地質(zhì)意義. 石油實(shí)驗(yàn)地質(zhì)， 32（2）： 101- 107

趙俊峰， 劉池洋， 劉永濤， 何爭光， 毛偉， 朱斌. 2011. 南華北地區(qū)上古生界熱演化史恢復(fù). 石油與天然氣地質(zhì)， 32（1）： 64- 74.

周安朝， 趙陽升， 郭進(jìn)京， 張寧. 2010. 西藏羊八井地區(qū)高溫巖體地?zé)衢_采方案研究. 巖石力學(xué)與工程學(xué)報(bào)，29（增刊2）： 89- 95.

周春光， 龔玉紅， 張惠良. 1996. 湖南中生代巖漿活動與晚二疊世煤變質(zhì)特征. 現(xiàn)代地質(zhì)， 10（4）： 470- 477.

朱傳慶， 徐明， 袁玉松， 趙永慶， 單竟男， 何志國， 田云濤， 胡圣標(biāo). 2010. 峨眉山玄武巖噴發(fā)在四川盆地的地?zé)釋W(xué)響應(yīng). 科學(xué)通報(bào)， 55（6）： 474- 482.

Ascencio F， Samaniego F and Rivera J. 2006. Application of a spherical-radial heat transfer model to calculate geothermal gradients from measurements in deep boreholes. Geothermics， 35： 70-78.

Barker C E. 1986. The correlation of vitrinite reflectance with maximum temperature in Humic organic matter // Lecture Notes in Earth Science. Berlin： Springer-Verlag， （5）.

Carslaw H S and Jaeger J C. 1959. Conduction of Heat in Solids. Oxford University Press， New York： 1- 386.

Duffield W A and Ruiz J. 1992. Evidence for the reversal of gradients in the uppermost parts of silicic magma reservoirs. Geology， 20： 1115- 1118.

Eldursi K， Branquet Y， Guillou-Frottier L and Marcoux E. 2009. Numerical investigation of transient hydrothermal processes around intrusions： Heat-transfer and fluidcirculation controlled mineralization patterns. Earth and Planetary Science Letters， 288： 70-83.

Fedotov S A. 1976. Uplift of mafic magma in Earth crust and mechanism of fracture basalt eruption. lzvestiya AN SSSR， Series Geology， 10： 5- 23 （in Russian）.

Feoktistov G D. 1972. Metamorphism of Clay-sandy Rocks near Contact Zone with Intrusion. Nauka， Moscow：1- 100 （in Ru ssian）.

Fjeldskaar W， Helset H M and Johansen H， et al. 2008. Thermal modeling of magmatic intrusions in the Gjallar Ridge， Norwegian Sea： Implications for vitrinite reflectance and hydrocarbon maturation. Basin Research，20： 143- 159.

Gunson M， Hall G and Johnston M. 2000. Foraminiferal coloration index as a guide to hydrothermal gradients around the Porgera intrusive complex， Papua New Guinea. Economic Geology， 95： 271-282.

Jaeger J C. 1961. The effect of the drilling fluid on temperature measured in bore holes. Journal Geophyscis Resources， 66（2）： 563-569.

Jaeger J C. 1964. Thermal effects of intrusions. Reviews in Geophysics， 2： 433-466.

Kol'tsov A B. 2010. Hydrothermal mineralization in the fields of temperature and pressure gradients. Geochemistry International， 48（11）： 1097-1111.

Ressel M W and Henry C D. 2006. Igneous geology of the Carlin Trend， Nevada： Development of the Eocene plutonic complex and significance for Carlin-type gold deposits. Economic Geology， 101： 347- 383.

Rubin A M. 1995. Propagation of magma-filled cracks. Annual Review of Earth and Planetary Sciences， 23： 287- 336.

Webber K L， Simmons W B， Falster A U and Foord E E. 1999. Cooling rates and crystallization dynamics of shallow level pegmatite-aplite dikes， San Diego County，California. American Mineralogist， 84： 708- 717.

Vitrinite Reflectance Method to Identify Concealed Intrusion： An Example of Application of the Magmatic Thermal Field Theory

ZHANG Qi1， JIN Weijun1， LI Chengdong2， WANG Jinrong3，JIAO Shoutao1，4and CHEN Wanfeng3
（1. State Key Laboratory of Lithospheric Evolution， Institute of Geology and Geophysics， Chinese Academy of Sciences， Beijing 100029， China； 2. Tianjin Center， China Geological Survey， Tianjin 300170， China； 3. School of Earth Sciences， Lanzhou University， Lanzhou 730000， Gansu， China； 4. Beijing Institute of Geology and Mineral Resources， Beijing 100012， China）

It is of crucial important to identify the undiscovered intrusions in hydrothermal deposit prospecting. There are lots of effective methods to identify undiscovered intrusions， such as magnetic， electric， and gravitational methods. However， the vitrinite reflectance method for the geothermal indicator is a direct and simple method. Magmatic emplacement is bound to form a thermal field around， and form a temperature gradient field with the distance of magma and superimposed on the geothermal field. Therefore， we can find undiscovered intrusions through analysis of the variation of thermal field. If we know the heat transfer parameters of magma and surrounding rocks， we even can quantitatively determine the depth， the location and the scale of the undiscovered intrusion. The vitrinite reflectance（expressed as a Ro） method commonly used in coal and petroleum geology is a simple， economic， and effective method to identify undiscovered intrusion. The effective distance of the method is around a few meters to a few kilometers，depending on the size， composition， and depth of different intrusion. Regardless of the new and old mining area， you can try this method so long as there is post-Silurian mud rock. This method applies to the surface can also be applied to drilling. For old mining area， there are the research data for three dimension images， for new ones， the depth of Rovariation should be noted in order to guide the prospecting work in time. The authors argue that the intrusion is commonly associated with the fluid flow in the mining area. If Roanomalies （Ro＞2%） are found in the mining area，attention should be paid to the possibility of undiscovered intrusion and the accompanying mineralization in geological expectations. The authors hope that the introduction of this method， combining with the characteristics of hydrothermal ore deposit to the study of ore mineralization， will carve out a new application and make new progress in the deep prospecting.

vitrinite reflectance； undiscovered intrusive； magma thermal field； method； mineralization； granite

P622

A

1001-1552（2015）06-1094-014

10.16539/j.ddgzyckx.2015.06.011

2014-10-20； 改回日期： 2015-05-27

項(xiàng)目資助： 國家自然科學(xué)基金重大研究計(jì)劃項(xiàng)目（91014001）資助。

張旗（1937-）， 男， 研究員。巖石學(xué)和地球化學(xué)專業(yè)。Email： zq1937@126.com