岩石学报  2014, Vol. 30 Issue (2): 341-349   PDF    
岩浆热场:它的基本特征及其与地热场的区别
张旗1, 金惟俊1, 李承东2, 焦守涛1     
1. 中国科学院地质与地球物理研究所 北京 100029;
2. 中国地质调查局天津地质矿产研究所, 天津 300170
摘要:“岩浆热场”指的是由岩浆引发的瞬间热场。热场的热主要来自未固结的岩浆,岩浆加热了围岩,使下地壳、中地壳和上地壳的下部在一个短暂的时间内保持一种高热状态。岩浆热场与地热场有许多不同:(1)热的来源不同。地热场的热主要来自地壳物质放射性生成的热;岩浆热场的热来自岩浆。(2)热的分布不同。地热场的等温面总体上呈水平分布,温度随深度增加而增加;岩浆热场的等温面则围绕岩体分布,靠近岩体温度高,远离岩体温度低,故岩浆热场的等温面是大体垂直于地热场等温面分布的。(3)热场的规模不同。地热场是全球性的,岩浆热场是局部性的,只在有岩浆的地方才出现。岩体小则规模小(热场宽度仅几米或几十米),岩体大则规模大(宽约几千米);如果存在大规模岩浆活动,岩浆热场的长宽均可达几百或上千千米,如在中国东部中生代大规模岩浆活动期间。(4)热持续的时间不同。地热场可以持续很长的时间(几十、几百或几千个百万年);岩浆热场是瞬间的突发性事件,持续的时间从几年到几个百万年。岩浆热场最重要的意义是,它是热液赖以上升的通道,它有利于来自下地壳底部和壳幔过渡带的流体(热液)的活动,使含矿热液得以顺利上升,并在热场范围内进行充分的活动、对流循环、萃取围岩中的成矿金属元素,并在地壳浅部岩浆热场之上合适的部位沉淀富集成矿。“岩浆热场”的概念依赖于对岩浆物理性质和过程的深入了解,由于我们这方面的知识相对贫乏,所以目前对岩浆热场的了解还是很肤浅的。
关键词地热场     岩浆热场     等温线     流体     区别     规模     瞬间    
Magma-thermal field:Its basic characteristics, and differences with geothermal field
ZHANG Qi1, JIN WeiJun1, LI ChengDong2, JIAO ShouTao1    
1. Institute of Geology and Geophysics, Chinese Academy of Sciences, Beijing 100029, China;
2. Tianjin Institute of Geology and Mineral Resources, China Geological Survey, Tianjin 300170, China
Abstract: The magma-thermal field indicates the instantaneous thermal field caused by magma. The heat in the thermal field mainly from unconsolidated magma has heated the surrounding rock. The lower crust, middle crust and the lower part of the upper crust maintains a high temperature state in a short time. There are many differences between a magma-thermal field and a geothermal field: (1) The different heat source. The heat in geothermal field is generated by the crustal radioactive material of the overlay rocks; otherwise the heat in a magma-thermal field comes from the magma. (2) The different distribution. The isotherm of the geothermal field was of horizontal, and it increases with increasing depth; the temperature contour of magma thermal fields is of surrounding about the intrusions, it is of high temperature close to intrusion, and lower one away from the intrusions, the isotherm of the magma thermal fields is perpendicular to the geothermal field isotherm. (3) The different scales of thermal field. The geothermal field is global; the magma thermal field is localized only where the magma exists, small scale with smaller intrusion (only a few meters or tens of meters in the width), large scale with larger intrusion (a few kilometers in the width). The magma thermal field is of up to hundreds or thousands of kilometers in the length and width if there is large-scale magmatic activity, such as large-scale magmatic activity in eastern China during the Mesozoic era. (4) The different duration of heat maintain. The geothermal field can last for a very long time (tens, hundreds or thousands of millions of years); the magma thermal field is instantaneous event, last from a few years to a few million years. It is the most important implication that the magma thermal field supplies the hydrothermal rise channel on which it is fever to fluid (hydrothermal) activities under the crust bottom and crust-mantle transition zone, so that the ore-bearing hydrothermal rise smoothly, active fully, convect, and extract mineralization metal elements in the surrounding rock, and precipitation enrichment and mineralization in a suitable site on the shallow crust within the thermal field scope. The concept 'Magma-thermal field' relies on the deep understanding of the physical properties and processes of magma, and is still very superficial due to our relative lack of knowledge.
Key words: Geothermal field;Magma-thermal field;Isotherms;Fluid;Difference;Scale;Instantly    
1 引言

“岩浆热场”指的是由岩浆引发的瞬间热场,这不是一个新概念。岩浆热场实际上指的就是岩浆对围岩的热效应,这个问题早在20世纪初就引起了国外学者的广泛重视。尤其是20世纪60年代Jaeger的成果发表以来(Jaeger, 1961, 1964),岩浆的热效应问题已成为地学界重要的科学问题。国外很早就开始尝试用数学计算对侵入体周围的热传导和热液活动过程进行模拟计算(Carslaw and Jaeger, 1959; Feoktistov, 1972; Fedotov, 1976; Delaney, 1987; Duffield and Ruiz, 1992; Rubin, 1995; Webber et al., 1999; Gunson et al., 2000; Ascencio et al., 2006; Eldursi et al., 2009; Kol’tsov, 2010)。常识告诉我们,炽热的岩浆侵位必定在周围形成一个热场(切列缅斯基,1977罗文积和陈家清,1997)。“岩浆热场”的术语也很早就出现在国内文献中了(魏洲龄,1994吕新彪等,1997罗文积和陈家清,1997高贵荣和张勉斌,1998饶超,2007邵飞,2007金旭东等,2010),许多人还讨论了现代和古代岩浆活动与地热场的关系(周江羽等,1997万志军等,2005杨兴科等,2005郭进京和周安朝,2008周安朝等,2010朱传庆,2010a, b赵慈平等,2012)。但是,把岩浆热场作为一个科学概念加以阐述却并不多见。其中,以罗文积和陈家清(1997)的研究令人刮目相看。他们早在20世纪90年代就明确表述了对“岩浆热场”的认识。他们指出,花岗岩浆侵入围岩必然要形成一个以岩浆熔体为中心的热场。花岗岩浆缓慢冷却的时间可达数百万年,花岗岩浆的温度大约在700~800℃以上,过热的岩浆还可更高些(切列缅斯基,1977Petford et al., 2000; 罗文积和陈家清,1997周安朝等,2010)。当这些高温融浆大规模侵入围岩时,它将有足够的时间和热容量使围岩增温形成热场,并引起围岩热力变质。与侵入体直接接触部分温度较高,随远离侵入体围岩温度逐渐降低,形成一个以岩浆熔体为中心的热场,并出现一个明显的由内向外递降的温度梯度。侵入体对围岩的热力影响有一个较大的范围。岩浆热场范围的大小和形状还与侵入体的温度、成分、形态、大小、侵入深度及围岩性质有关(罗文积和陈家清1997)。

众所周知,中国东部中生代存在大规模岩浆活动和大规模成矿作用(陈衍景等,2004毛景文,20052007),但是,为什么大规模岩浆活动与大规模成矿作用有密切的关系?按照流行的见解这似乎不成为问题,因为流行的理论认为,岩浆与成矿在时间、空间和成因上有关,热液和成矿金属皆来源于岩浆。但是,这种解释存在几个困难:(1)这种解释的前提是花岗岩质岩浆必须富水,而且能够分异和演化,在岩浆固结以后温度压力降低能够分出大量的流体,花岗岩体越大分出的流体越多。上述理论与花岗岩的实际不符。实际情况是:花岗岩绝大多数是在缺水条件下部分熔融的(Rutter and Wyllie, 1988; Thompson et al., 1995; Harris et al., 1995; 刘勇胜和高山,1998张旗,2012),侵位的花岗岩是干的岩浆,没有水分出,尤其越是大的花岗岩越没有水被分出的现象(花岗岩边部出现的是烘烤现象和接触变质作用),水仅出现在少数小岩体及其边部。(2)野外存在岩浆活动与成矿作用不匹配的现象,如胶东存在3期岩浆活动,而成矿只与后两期岩浆活动有关,最早的一期无矿化。广西存在加里东期、海西期、印支期和燕山期花岗岩,而钨锡成矿主要与最后一期岩浆活动有关,虽然上述各期花岗岩的地球化学性质接近(均为有利于钨锡成矿的南岭型花岗岩,张旗和李承东,2012)。(3)研究表明,在下地壳底部,流体的出现早于岩浆(张旗,2012),而大量的矿床实例却告诉我们,成矿往往在成岩之后。

我们知道,流体是极其活泼的,流体从下地壳底部上升到地壳浅部的速度可能比岩浆快得多,但是,需要能够保持一个高温状态。如果温度很快下降了,流体可能迅速冷却就无力继续上升了。含矿热液相当于水,但不是水;水在0℃结冰,不同成分的含矿热液大约在500~150℃范围内沉淀。因此,为了流体的上升,前提是需要在地壳内存在一个高温热场,而且需要保持一个相当长的时间(如几年、几万年或几个百万年),以使热液有足够的时间上升、对流、循环、萃取赋存于围岩中成矿元素并最终在地壳浅部沉淀成矿。构造、岩浆都有可能造就一个高温热场。但是,相对而言,岩浆的侵入是形成高温热场的最佳可能。因为岩浆(尤其是酸性岩浆)上升的速度慢、降温慢、向围岩散热范围广、持续时间长。热场持续的时间取决于岩浆冷却的速度。一般来说,热场存续的时间很短暂,或许只有几年,几万年或几百万年,上述时间对于地质演化来说仅仅是一瞬间。因此,岩浆热场属于突发性热事件产生的瞬间热场,是一个由岩浆引发的规模很小,时间很短的热场。为了区别于“地热场(Geothermal field)”,本文称其为“岩浆热场(Magma-thermal field)”。

“岩浆热场”的确立很主要的一个方面是依赖于对岩浆物理性质和岩浆动力学的了解。例如,岩浆是怎么形成的?怎么侵位的?上升和冷却的速度是多少?向围岩散发热量有多少?最远距离多少?保持的时间多少?新叠加的热场对围岩有什么作用?如何使围岩中的水进一步被释放出来?释放出来的水是个什么状态?在热场中水是怎么循环对流的?对于成矿有什么作用?等等。由于我们对岩浆物理性质和岩浆动力学研究的长期忽视,只有少数学者在这方面作了一些探索(马昌前,1986; 马昌前等, 1994吴传荣,1992周珣若等,1993周涛发等,1995冯乔和汤锡元,1997杨坤光和杨巍然,1997张健和石耀霖,1997王涛,1999a, b万天丰等,2000冯佐海,2003罗照华等,2007章邦桐等,2007李德东等,2011Zhao et al., 2003);而国外对这方面的研究比我们早50~100年,而且取得了重要的研究成果(Barner, 1979; Turcott and Schubert, 1982; McBirney and Murase, 1984; Sparks et al., 1984; Marsh, 1982; Castro, 1987; Paterson and Tobisch, 1992; Pitcher, 1987, 1993; Evans et al., 1994; Speer et al., 1994; Petford et al., 1993, 2000; Kay and Kay, 1993; Weinberg and Podladchikov, 1994; Roman et al., 1997, Roman-Berdiel et al., 1995; Ferré et al., 1995; Hutton et al., 1990)。例如,许多人对火成岩侵入自身的热性质和热效应进行了大量的研究(Jaeger, 1961, 1964; Hall, 1971; Hildreth, 1981; Brandeis and Jaupart, 1986; Webber et al., 1999; Gerya and Burg, 2007);对不同形状侵入体(包括岩墙、岩床、岩株、岩筒、岩盖以及一些不规则的侵入体)的热效应进行了二维、三维模拟(Jaeger, 1964; Simmons, 1967; Guerrero-Martínez et al., 2013);计算了瞬时侵入导致的热影响的规模(Carslaw and Jaeger, 1959; Feoktistov, 1972; Fedotov, 1976; Delaney, 1987; Duffield and Ruiz, 1992);解释了岩墙和岩床边部缺少变质晕的现象(Kontorovich et al., 1981; Delaney and Pollard, 1982; Raymond and Murchison, 1988; Thrasher, 1992);对围绕侵入体的热液过程、热传导、流体循环以及成矿作用进行了研究(Rubin, 1995; Webber et al., 1999; Gunson et al., 2000; Ascencio et al., 2006; Eldursi et al., 2009; Kol’tsov, 2010);准确评估了热场对围岩有机质热成熟度的影响(Galushkin, 1997)等。国内对于这方面的研究起步很晚,至今也很少受到关注(马昌前,1986; 马昌前等,1994吴传荣,1992周珣若等,1993周涛发等,1995冯乔和汤锡元,1997杨坤光等,1997张健和石耀霖,1997王涛,1999a, b万天丰等,2000冯佐海,2003罗照华,20072009,2011章邦桐等,2007李德东等,2011Zhao et al., 2003)。

本文不属于综述性或评论性文章,而是笔者研究花岗岩与成矿关系的一点体会、感想或猜想。因为,岩浆热场理论是建立在岩浆物理学、岩浆动力学和流体动力学基础上的,它需要专门的知识,而笔者于这些方面是门外汉。由于它可能对成矿有特殊的意义,笔者才关注它。本文只是抛砖引玉,希望能够引起有关专家和有兴趣同行的关注和批评。

2 岩浆热场的基本特征

据我们的初步认识,岩浆热场可能具备如下几个特征:

(1)它是由岩浆引起的,是岩浆侵位时带来的热导致的。罗文积和陈家清(1997)指出,花岗岩浆侵入围岩必然要形成一个以岩浆熔体为中心的热场。花岗岩浆的温度大约在700~800℃以上,过热的岩浆还可更高些。当这些高温融浆大规模侵入围岩时,它将有足够的时间和热容量使围岩增温形成热场,并引起围岩热力变质。与侵入体直接接触部分温度较高,随远离侵入体围岩温度逐渐降低,形成一个以岩浆熔体为中心的热场,并出现一个明显的由内向外递降的温度梯度。侵入体对围岩的热力影响有一个较大的范围。岩浆热场范围的大小和形状还与侵入体的温度、成分、形态、大小、侵入深度及围岩性质有关。一般认为,地球内热有以下3种传送方式:热传导、热辐射和热对流。热传导不是地球内热传导的主要方式,因为岩石是极不良的热导体,热导率很小。由于在近红外和可见光区域内硅酸盐矿物不易发生辐射,故热辐射也不可能是地球内热传送的主要方式。于是,热对流即成为地球内热传导的主要方式(Roy et al., 1972; Villas and Norton, 1977; Torcotte and Schubert, 1982; Chapman and Rybach, 1985; Best and Christiansen, 2001; 切列缅斯基,1977马东升,1998万志军等,2005章邦桐等,2007王满等,2012张旗等,2013)。

(2)它是瞬间发生的,只代表岩浆侵入及其冷却过程中发生的事件,不包括岩浆冷却固结后由于放射性元素蜕变产生的热。由于岩浆本身携带的热量不同,岩浆性质的不同,岩浆冷却时间的不同,这个瞬间可以短到几千年,也可能短到几个百万年。唐晓音等(2013)最近对琼东南盆地长昌凹陷烃源岩成熟度的研究表明,一个面积为300km2厚约10km的侵入体,对围岩温度场有显著影响的时限不超过1Myr,对烃源岩有机质成熟度影响的最大距离不超过2km。Petford et al. (2000)认为,岩浆上升侵位的时间可能不到0.1Myr。笔者认为这个时间可能不够,因为,酸性岩浆的黏性很高,而且岩浆并不是统统沿断裂侵入的,许多岩体呈浑圆状产出,如房山岩体、姑婆山岩体等,表明它们是底辟侵位的。因此,岩浆侵位及其冷却时间在几个百万年以内是可以接受的。

(3)热异常和等温面大体是垂直分布的,就如一个盐丘底辟造成的等温面的分布一样(图 1,引自切列缅斯基,1977)。由于许多岩体的侵入,使热场的形状更加复杂一些而已。这一点明显不同于大致呈水平分布的地热场,这是由于岩浆的垂直上升决定的,如图 2所示。

图 1 盐丘的地热剖面
灰色示盐丘,实线示等温线
Fig. 1 Geothermal profile of salt domes

图 2 岩浆热场示意图
图中呈水平分布的为地热场,岩浆热场等温线大体为垂直分布,切割地热场.花岗岩来自下地壳底部部分熔融区
Fig. 2 Schematic diagram of magma-thermal field

(4)热场的规模很小,通常只离岩体几米或几千米。与岩体的规模、性质,围岩的传热性以及流体的有无有关。如果考虑温度可以降低至300~100℃,足够引发煤质的改变以及油气藏的变化,其影响的范围可能达到几、十几至几十千米。瞬间侵入的岩体对围岩温度的影响,在有孔隙水气化的情况下估计的岩席-围岩的接触温度Tc最高,达到852℃,无孔隙水也能达到818℃,而缓慢侵入的热传导模型计算的接触温度最低,仅为706℃(王大勇等,2011)。对比模拟结果和镜质体反射率测量数据表明,瞬间岩浆的热影响的范围不超过75m,大约为侵入体厚度的5倍。看来,小岩墙、岩株、岩脉影响的范围很小,顶多不超过岩体接触带几米或十几米;而岩基的影响范围就大一些,有些接触变质带的宽度可达1km,热场则可超过岩体边界几千米。因此,岩浆热场的规模变化很大。如果是一个辉绿岩墙,由于它的黏性低,侵位的速度快,所形成的热场的宽度就大大不如花岗岩岩墙,大约只有几米或几十米,继续的时间也非常短暂(图 2a)。50m宽的侵入岩体的有效影响范围约200m,对于100m宽的侵入岩体的有效影响范围约500m(王满等,2012)。规模较大的热场可能由一个大岩体或岩基引起,热场的宽度比岩墙大多了,持续的时间也长多了(图 2b)。图 2c为规模巨大的岩浆热场,由大规模的岩浆活动引起,大规模岩浆活动必定有多期岩浆的侵入,图中分为4期,分别以红色、蓝色、绿色和黄色表示,持续的时间约几个或十几个百万年。

(5)热场可大致可分为高温热场和低温热场,笔者估计高温热场的温度至少>600~500℃(图 2b和c表示了温度为600℃和400℃的热场范围),在这个温度以上,可能有利于流体的运动,有利于流体从围岩中萃取出金属组分;<400℃为低温热场。低温热场的温度虽然低,但是,能够明显改变煤和石油的品质,对于沉积热液矿床的富集可能也是有意义的。

总之,“岩浆热场”指的是:在一个很短的时间内,在一个地区出现的岩浆活动,可能会造成该区域地热梯度相对周边地区明显的上升,使之形成一个局部的热场。热场的热主要来自未固结的岩浆,岩浆加热了围岩,使下地壳、中地壳和上地壳的下部在一个短暂的时间内保持一种高热状态。

岩浆热场最重要的意义是,它是热液赖以上升的场所。这样的热场有利于来自下地壳底部和壳幔过渡带的流体(热液)的活动(张旗,2012),使含矿热液得以顺利上升,并在热场范围内进行充分的活动、对流循环、萃取围岩中的成矿金属元素,并在地壳浅部岩浆热场之上合适的部位沉淀富集成矿。

3 岩浆热场与地热场的区别

岩浆热场不同于地热场,它们的区别是:

(1)热的来源不同。地热场以温度、地热梯度、地热流密度以及热源放热(吸热)的强度及其分布来表示。地壳及其地表的热力场取决于宇宙和地球内部的热源以及地球内部的动力作用,受多种因素制约。地壳的热力场是由放热和吸热热源的强度和分布、地层的热力学性质以及地表温度决定的(切列缅斯基,1977)。多数人认为地壳和地幔所含放射性元素的衰变是地球内部热源最主要的来源(Roy et al., 1972; 切列缅斯基,1977从柏林,1978)。而岩浆热场的热主要来自岩浆的侵入,包括基性和酸性的岩浆,但主要是中酸性的岩浆。

(2)热的分布不同。地热场的温度是随地壳深度递增的,上部温度低,深部温度高,因而是近水平分布的。地热场的分布是大体均匀的,在一个地方测定的地热增温率可以代表一个较大区域的热状况。而岩浆热场的热的分布是极不均匀的,主要取决于岩浆的性质、岩浆的规模、岩体的形状,并且可能受构造活动的影响。岩浆热场的热是围绕岩体分布的,靠近岩体温度高,远离岩体温度低,而不是受地壳深度控制的,故岩浆热场的等温线大体是垂直分布的,不同于地热场的水平分布。岩浆热场的顶部形态因岩体侵入的位置不同而变。此外,岩浆热场的分布还与构造有关,如果有开放性的构造,可以把热传递到远离岩体的部位。但对于岩浆热场来说,这种现象是局部的。

(3)地温梯度不同。地热场的温度随深度增加,地温梯度总体上却有一个随深度增加减小的趋势(图 3a)。也就是说,在地壳浅部,地温梯度可能比较大,随着深度增加,地温梯度变小了(万志军等,2005)。由岩石本身放射性元素获得的地温曲线直达下地壳底部,温度可能也不会超过800℃,依靠地温梯度的增加可能不会导致花岗岩的部分熔融。在某些地区,地球物理探测如果深部某个部位存在一个或几个低速带(由未固结的岩浆引起),地温梯度将在低速带处明显升高,过了低速带,温度下降,地温梯度曲线为瘤状(图 3b)。岩浆热场不同,地温曲线很复杂:在地表附近,地温梯度即异常的高(图 3C,假定为50℃/km),进入岩浆热场范围地温梯度进一步升高(达到约90℃/km),到达侵入体顶面温度有一个跳跃式的升高,即深度不增加,温度可增加100~200℃或更高。这是因为侵入体接触部分与围岩之间存在一个明显的温度差所致。进入未固结的岩浆后,地温梯度即急剧下降(<10℃/km),一直到下地壳底部,温度仅升高100~200℃(图 3c)。

图 3 示意的地热场与岩浆热场地温梯度对比
(a)-正常地热场地温梯度分布;(b)-地热场中有低速带存在的分布;(c)-岩浆热场地温梯度分布.详细说明见正文
Fig. 3 Schematic comparing to geothermal gradient between geothermal field and magma thermal field

(4)热场的规模不同。地热场是全球性的,在全球任何地方都存在,只是温度梯度不同而已。岩浆热场是局部性的,只在有岩浆活动的地方才出现,且规模很小,无法与地热场相匹敌。岩浆热场的规模小到几米或几十米,由岩墙的侵入所引起;大的与大规模岩浆活动有关,长宽均可达几百至上千千米,如中国东部中生代大规模岩浆活动期间(示意见图 2)。

(5)热持续的时间不同。地热场可以持续很长的时间,如果没有构造变动,可以持续几十、几百或几千个百万年。朱传庆等(2010a, b)利用石油钻井的Ro(镜质体反射率)资料采取古热流恢复方法得出的四川盆地的热流史表明,在加里东期之前,该区热状态是稳定的,热流值较低;至海西期热流开始增大,在259Ma左右热流值达到最高,多数钻井的最高古热流值在60~80MW/m2之间,少数钻井经历的最高古热流超过了100MW/m2,此后热流持续降低直到现今(图 4)。图 4中从加里东期至今,地热场大体不变或变化不大,原因是维持地热场的放射性热的来源没有变。在259Ma出现的高热流属于岩浆热场,其余时间段属于地热场。岩浆热场是由于峨眉山玄武岩的喷发造成的,反映了东吴运动期间玄武岩喷发时岩浆活动导致的(瞬时)热效应朱传庆等(2010a, b)。

图 4 四川盆地代表性钻井的古热流(据朱传庆等,2010a) Fig. 4 The drilling representation paleo-heat flow of Sichuan basin (after Zhu et al., 2010a)

(6)热与水的关系不同。地热场提供的热是使地层中的水释放出来,随着深度的增加和热的增加,赋存在地层中的水逐渐减少。地热是将水向上驱赶,温度越高水含量越低。岩浆热场不同,通常伴有明显的热液活动。热液一部分来自深部(下地壳底部、壳幔过渡带),一部分来自温度升高对地层中原有水的萃取。如果构造条件合适,岩层流通性好,热液可以在热场中进行充分的对流循环,汲取围岩中的金属元素,成为金属矿床的来源之一(Norton and Taylor, 1975; Ferry, 1984; 罗文积和陈家清1997马东升,1998张映红和顾家裕,2003)。在火山岩发育地区和热泉活动的地区的热虽然归入地热场的范畴,实际上是来自岩浆热场的(与侵入和喷出的岩浆岩有关)。

(7)研究方法不同。地热场是地球物理学研究的内容之一,现代地热场可以通过地表和钻孔测量以及地球物理方法得到。而岩浆活动时期的岩浆热场则根据岩浆活动带给围岩的一系列变化来估计。包括岩浆岩的规模、岩浆侵入的深度、热液蚀变的宽度、构造活动的强度、围岩孔隙开放的程度等。

4 初步的认识

由于岩浆热场的存在,可造成地壳局部地区温度的急剧升高,高温(>700~800℃)可从地壳浅部开始,维持达几十千米的深度,这种现象是极具吸引力的。但是,关于岩浆热场我们的认识还很肤浅。例如,岩浆热场的边界在哪里?以多高的温度为标志比较合适?岩浆热场的范围多大?受什么因素制约?在诸多因素中哪些是主要的,哪些是次要的,哪些因素在哪些情况下是变化的?

瞬间热场的出现把局部地区的地温梯度急剧提高,它对地热场会产生什么影响?流体和热液在岩浆热场中是怎样运动和循环的?对温度的限制是多少?哪些因素有利于流体的对流循环?哪些情况下它可以从围岩中萃取出金属元素来?次生金属的加入会使热液发生什么变化?包裹体测温测到的700~800℃的高温是发生在什么地方?热场与成矿是什么关系?成矿在热场之内还是之外?在什么情况下沉淀成矿?这些,已经有过一些研究,但是,许多问题仍然是需要继续探索的。

令人惊叹的是,岩浆热场说对于热液成矿可能具有非常重要的意义,这也是笔者为什么努力去探索岩浆热场说的原因之一。例如,与花岗岩有关的金铜与钨锡成矿有些是相悖的,有些是相伴的,热场说可能会给出一个新的解释。热场说还解释了为什么大多数矿床是多金属成矿的问题。关于多金属成矿学术界已经有很多很好的讨论,笔者相信,很多矿床中的低品位的伴生金属元素可能来自于热场范围内的围岩,是被在热场中活动的流体萃取出来的。热场说还解决了成矿为什么总是晚于成岩的问题,为什么大规模岩浆活动与大规模成矿作用密切相关的问题,它也从一个新的角度解释了“就矿找矿”的理论为什么仍然适用的问题。此外,热场说还可以解释一些沉积热液矿床的成矿问题、如某些卡林型金矿问题、某些层控SEDEX铅锌矿床再富集问题、某些远离侵入岩的夕卡岩成因问题。不仅如此,热场说对煤层热演化和煤层气生成也有显著的控制作用,对油气的生成、运移、聚集和油气藏的形成与保存也有明显影响,还可能解释某些地区的地热异常问题。

总之,岩浆热场是一个极具创新性的课题,其意义究竟怎样,还有许多是我们不知道的。本文只是揭开了岩浆热场面貌的一角,更加精彩的部分正等待人们进一步地去挖掘。

后记

岩浆热场问题是笔者从思考花岗岩与成矿关系中悟出来的,这只是一个猜想,它是否可信,是允许怀疑和讨论的。笔者觉得,在自然界可能存在一个瞬间的岩浆热场,它与岩浆活动有关,与成矿的关系可能也很密切。岩浆热场说建立在地热学、岩浆物理学、岩浆动力学和流体动力学理论的基础上,笔者于上述几方面均是外行,虽然临时抱佛脚学了一点点,只是皮毛而已,欢迎大家评头论足。科学研究的本质是探索,是向未知领域探求真理。岩浆热场说是否有价值,希望今后的实践予以检验。论文撰写过程中向赵平研究员讨教过意见,两位审稿人提出了很好的意见,罗照华教授与笔者交流、提供了许多重要的参考文献,使笔者受益匪浅,在此一并表示衷心的感谢。

参考文献
[1] 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(1): 70-78
[2] Barner HL. 1979. Geochemistry of Hydrothermal Ore Deposits. 2nd Edition. A Wiley-Interscience Publication, John Wiley and Sons
[3] Best MG and Christiansen EH. 2001. Igneous Petrology. Blackwell Science Ltd., 1-458
[4] Brandeis G and Jaupart C. 1986. On the interaction between convection and crystallization in cooling magma chambers. Earth and Planetary Science Letters, 77(3-4): 345-361
[5] Carslaw HS and Jaeger JC. 1959. Conduction of Heat in Solids. New York: Oxford University Press, 1-386
[6] Castro A. 1987. On granitoid emplacement and related structures: A review. Geol. Rundsch., 76(1): 101-124
[7] Chapman DS and Rybach L. 1985. Heat flow anomalies and their interpretations. J. Geodynamics, (4): 3-37
[8] Chen YJ, Pirajno F, Lai Y et al. 2004. Metallogenetic time and tectonic setting of the Jiaodong gold province, eastern China. Acta Petrologica Sinica, 20(4): 907-922 (in Chinese with English abstract)
[9] Cong BL. 1978. Magmatism and Igneous Associations. Beijing: Geological Publishing House (in Chinese)
[10] Delaney PT and Pollard DD. 1982. Solidification of basaltic magma during flow in a dike. Amer. J. Sci., 282(6): 856-885
[11] Delaney PT. 1987. Heat transfer theory applied to mafic dike intrusions. In: Halls HC and Fahrig WF (eds.). Mafic Dike Swarms. Geol. Assoc. Canada, Spec. Paper., 34: 31-46
[12] Duffield WA and Ruiz J. 1992. Evidence for the reversal of gradients in the uppermost parts of silicic magma reservoirs. Geology, 20(12): 1115-1118
[13] Eldursi K, Branquet Y, Guillou-Frottier L and Marcoux E. 2009. Numerical investigation of transient hydrothermal processes around intrusions: Heat-transfer and fluid-circulation controlled mineralization patterns. Earth and Planetary Science Letters, 288(1-2): 70-83
[14] Evans DJ, Rowley WJ, Chadwick RA et al. 1994. Seismic reflection data and the internal structure of the Lake District batholith, Cumbria, northern England. Proc. Yorkshire Geol. Soc., 50: 11-24
[15] Fedotov SA. 1976. Uplift of mafic magma in Earth crust and mechanism of fracture basalt eruption. lzvestiya AN SSSR, Ser. Geol., 10: 5-23 (in Russian)
[16] Feng Q and Tang XY. 1997. The magma activity's influence on conditions forming oil and gas accumulation. Northwest Geoscience, 18(1): 56-62(in Chinese with English abstract)
[17] Feng ZH. 2003. Emplacement process and structural analysis of Guposhan-Huashan granitic pluton, Guangxi. Ph. D. Dissertation. Changsha: Central South University (in Chinese with English summary)
[18] Feoktistov GD. 1972. Metamorphism of Clay-sandy Rocks near Contact Zone with Intrusion. Moscow: Nauka, 1-100 (in Russian)
[19] Ferré E, Gleize G, Bouchez J et al. 1995. Intenal fabric and strike-slip emplacement of the Pan-African granite of Solli Hills, northern Nigeria. Tectonics, 14(5): 1205-1219
[20] Ferry M. 1984. Reaction progress: A monitor of fluid-rock interaction during metamorphic and hydrothermal events. In: Walther JV and Wood BJ (eds.). Fluid-Rock Interactions during Metamorphism. New York: Springer-Verlag, 60-89
[21] Galushkin Y. 1997. Thermal effects of igneous intrusions on maturity of organic matter: A possible mechanism of intrusion. Org. Geochem., 26(11-12): 645-658
[22] Gao GR and Zhang MB. 1998. Mineralization characters ore-controlling and ore-searching prospect in Yinkeng multimetal orefield, Jiangxi Province. Volcanology & Mineral Resources, 19(4): 347-356(in Chinese with English abstract)
[23] Gerya TV and Burg JP. 2007. Intrusion of ultramac magmatic bodies into the continental crust: Numerical simulation. Physics of the Earth and Planetary Interiors, 160(2): 124-142
[24] Guerrero-Martínez FJ and Verma SP. 2013. Three dimensional temperature simulation from cooling of two magma chambers in the Las Tres Vírgenes geothermal eld, Baja California Sur, Mexico. Energy, 52: 110-118
[25] 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(2): 271-282
[26] Guo JJ and Zhou AC. 2008. Discussion on the determination methods and strategy for hot dry rocks. Energy Engineering, (6): 1-4(in Chinese with English abstract)
[27] Hall A. 1971. The Relationship between geothermal gradient and the composition of granitic magmas in orogenic belts. Contr. Mineral. Petrol., 32: 186-192
[28] Harris N, Ayres M and Massey J. 1995. Geochemistry of granitic melts produced during the incongruent melting of muscovite: Implications for the extraction of Himalayan leucogranite magmas. J. Geophys. Res., 100(B8): 15766-15777
[29] Hildreth W. 1981. Gradients in silicic magma chambers: Implications for lithospheric magmatism. J. Geophys. Res., 86(BⅡ): 10153-10192
[30] Hutton DHW, Dempster TJ, Brown PE et al. 1990. A new mechanism of granite emplacement: Intrusion in active extensional shear zones. Nature, 343: 452-455
[31] Jaeger JC. 1961. The effect of the drilling fluid on temperature measured in bore holes. J. Geophys. Res., 66(2): 563-569
[32] Jaeger JC. 1964. Thermal effects of intrusions. Reviews in Geophysics, 2(3): 433-466
[33] Jin XD, Zhang DH and Wan TF. 2010. Upper part and deeper area ore-prospecting for hidden rock mass. Geological Bulletin of China, 29(2-3): 392-400(in Chinese with English abstract)
[34] Kay RW and Kay SM. 1993. Delamination and delamination magmatism. Tectonophysics, 219(1-3): 177-189
[35] Kol'tsov AB. 2010. Hydrothermal mineralization in the fields of temperature and pressure gradients. Geochemistry International, 48(11): 1097-1111
[36] Kontorovich AE, Surkov VC, Troftmuk AA et al. 1981. Oil and Gas Geology of West Siberian Platform. Moscow: Nedra, 1-550 (in Russian)
[37] Li DD, Luo ZH, Zhou JL, Yang ZF and Liu C. 2011. Constraints of dike thicknesses on the metallogenesis and its application to the Shihu gold deposit. Earth Science Frontiers, 18(1): 166-178 (in Chinese with English abstract)
[38] Liu YS and Gao S. 1998. Anatexis of the continental crust and granites as tracers of lower crustal composions. Geological Science and Technology Information, 17(3): 31-38 (in Chinese with English abstract)
[39] Luo WJ and Chen JQ. 1997. Two-way convergence and hydrothermal mineralization. Acta Geologica Gansu, 6(Suppl.): 44-49 (in Chinese with English abstract)
[40] Luo ZH, Huang ZM and Ke S. 2007. An overview of granitoid. Geological Review, 53(Suppl.): 180-226 (in Chinese with English abstract)
[41] Luo ZH, Lu XX and Chen BH. 2009. Introduction to the Metallogenic Theory on the Transmagmatic Fluids. Beijing: Geological Publishing House, 1-168 (in Chinese with English abstract)
[42] Luo ZH, Liu JQ, Zhao CP, Guo ZF, Cheng LL, Li XH and Li DP. 2011. Deep fluids and magmatism: The deep processes beneath the Tengchong volcano group. Acta Petrologica Sinica, 27(10): 2855-2862(in Chinese with English abstract)
[43] Lv XB, Yao SZ and Zhou ZG. 1997. Magmatic rock system and metallogenesis in Jiurui region, Jiangxi Province. Geology and Mineral Resources of South China, (1): 27-36(in Chinese with English abstract)
[44] Ma CQ, Yang KG, Tang ZH et al. 1994. Magma Dynamics of Granitoids: Theory, Method and Case Study of the Eastern Hubei Granitoids. Wuhan: Press of China University of Geosciences, 1-260 (in Chinese with English abstract)
[45] Ma CQ. 1986. Some kinetic parameters estimation methods in magmatic activity: The application of the principles of fluid dynamics. Geological Science and Technology Information, 5(3): 47-54 (in Chinese with English abstract)
[46] Ma DS. 1998. Large-scale fluid flow systems in the crust and their implications for metallogenesis of hydrothermal ore deposits. Geological Journal of China Universities, 4(3): 250-261 (in Chinese with English abstract)
[47] Mao JW, Xie GQ, Zhang ZH et al. 2005. Mesozoic large-scale metallogenic pulses in North China and corresponding geodynamic settings. Acta Petrologica Sinica, 21(1): 169-188 (in Chinese with English abstract)
[48] Mao JW, Xie GQ, Guo CL et al. 2007. Large-scale tungsten-tin mineralization in the Nanling region, South China: Metallgenic ages and corresponding geodynamic processes. Acta Petrologica Sinica, 23(10): 2329-2338 (in Chinese with English abstract)
[49] Marsh BD. 1982. On the mechanics of igneous diapirism, stoping and zone melting. Am. J. Sci., 282(6): 808-855
[50] McBirney AR and Murase T. 1984. Rheological properties of magmas. Ann. Rev. Earth Planet. Sci., 12: 337-357
[51] Norton D and Taylor HP. 1975. Quantitatives simulation of the hydrothermal systems of crystallizing magmas on the basic of transport theory and oxygen isotope data: An analysis of the Skaergaard intrusion. Journal of Petrology, 20: 421-486
[52] Paterson SR and Tobisch OT. 1992. Rates of processes in magmatic arcs: Implications for the timing and nature of pluton emplacement and wall rock deformation. Journal of Structural Geology, 14(3): 291-300
[53] Petford N, Kerr RC and Lister JR. 1993. Dyke transport of granitoid magma. Geology, 21(9): 845-848
[54] Petford N, Cruden AR, McCaffrey KJW et al. 2000. Granite magma formation, transport and emplacement in the Earth's crust. Nature, 408: 669-673
[55] Pitcher WS. 1987. Granites and yet more granites forty years on. Geol. Rund., 76(1): 51-59
[56] Pitcher WS. 1993. The Nature and Origin of Granite. Glasgow and London: Blackie Academic, 182-190
[57] Qieliemianski KA. 1977. Applied Geothermal Science. In: Zhao Y and Chen M (Trans.). 1982. Beijing: Geological Publishing House (in Chinese)
[58] Rao C. 2007. Numerical simulation of heat and liquid transfer in the hydrothermal ore-forming progresses. Master Degree Thesis. Chengdu: Chengdu University of Tecnology (in Chinese with English summary)
[59] Raymond AC and Murchison DG. 1988. Development of organic maturation in the thermal aureoles of sills and its relation to sediment compaction. Fuel, 67(12): 1599-1608
[60] Roman-Berdiel T, Pueyo-Morer EL and Casas-Sainz AM. 1995. Granite emplacement during contemporary shortening and normal faulting: Structure and magnetic study of the Veiga Massif (NW Spain). J. Struct. Geol., 17(12): 1689-1706
[61] Roman-Berdiel T, Gapais D and Brun JP. 1997. Granite intrusion along strike-slip zones in experiment and nature. Am. J. Sci., 297(6): 651-678
[62] Roy RF, Klackwell DD and Decker ER. 1972. Continental heat flow. In: Robertson EC (ed.). The Nature of the Solid Earth. New York: McGraw-Hill, 506-543
[63] Rubin AM. 1995. Propagation of magma-filled cracks. Annual Review of Earth and Planetary Sciences, 23: 287-336
[64] Rutter MJ and Wyllie PJ. 1988. Melting of vapour-absent tonalite at 10kbar to simulate dehydration-melting in the deep crust. Nature, 331(6152): 159-160
[65] Shao F. 2007. Study on water-rock interaction and its relation with uranium metallogenesis. Ph. D. Dissertation. Wuhan: China University of Geosciences (in Chinese with English summary)
[66] Simmons G. 1967. Interpretation of heat flow anomalies 2. flux due to initial temperature of intrusives. Reviews of Geophysics, 5(2): 109-120
[67] Sparks RSJ, Huppert HE, Turner JS et al. 1984. The fluid dynamics of evolving magma chambers. Phil. Trans. R. Soe. Lond., A310(1514): 511-534
[68] Speer JA, McSween HY Jr and Gates AE. 1994. Generation, segregation, ascent, and emplacement of Alleghanian pluton in the southern Applachians. The Journal of Geology, 102(3): 249-267
[69] Tang XY, Zhang GC, Liang JS, Yang SC, Rao S and Hu SB. 2013. Influence of igneous intrusions on the temperature field and organic maturaty of the Changchang sag, Qiongdongnan basin, South China Sea. Chinese J. Geophys., 56(1): 159-169 (in Chinese with English abstract)
[70] Thompson AB, James A and Connolly D. 1995. Melting of the continental crust: Some thermal and petrological constraints on anatexis in continental collision zones andother tectonic settings. J. Geophys. Res., 100(B8): 15565-15579
[71] Thrasher J. 1992. Thermal effect of the Tertiary Cuillins intrusive complex in the Jurassic of the Hebrides: An organic geochemical study. In: Parnell J (ed.). Basins on the Atlantic Seaboard: Petroleum Geology, Sedimentology, and Basin Evolution. Geol. Soc. Spec. Publ., 62: 35-49
[72] Turcotte DL and Schubert G. 1982. Geodynamics: Application of Continuum Physics to Geological Problems. New York: John Wiley and Sons. Inc.
[73] Villas RN and Norton DL. 1977. Irreversible mass transfer between circulating hydrothermal fluids and the Mayflower Stock. Econ. Geol., 72(8): 1471-1540
[74] Wan TF, Teyssier C, Zhen HL et al. 2001. Emplacement mechanism of Linglong granitoid complex, Shandong Peninsula, China. Science in China (Series D), 44(6): 535-544
[75] Wan ZJ, Zhao YS and Kang JR. 2005. Simulation and forecast method of geothermal resources in hot dry rock. Chinese Journal of Rock Mechanics and Engineering, 24(6): 945-949 (in Chinese with English abstract)
[76] Wang DY, Lu XC, Xu SJ et al. 2011. Heat-transfer-model analysis of thermal effect of intrusive sills on organic-rich host rocks in sedimentary basins. Journal of Nanjing University (Natural Sciences), 47(1): 45-50 (in Chinese with English abstract)
[77] Wang M, Wang YW and Xue FL. 2012. Finite element simulation of effect of magmatic intrusion on maturity of organic matter in surrounding rock. Fault-Block Oil and Gas Field, 19(2): 172-176(in Chinese with English abstract)
[78] Wang T, Wang XX and Li WP. 1999a. Multiple emplacement mechanism and space of granitoid plutons. Geological Review, 45(2): 142-150 (in Chinese with English abstract)
[79] Wang T, Zhang GW, Wang XX et al. 1999. Growth partterns of granitoid plutons and their implications for tectonics, kinematics and dynamics: Examples from granitoid plutons in the core of the Qinling orogenic belt, China. Chinese Journal of Geology, 34(3): 326-335 (in Chinese with English abstract)
[80] Webber KL, Simmons WB, Falster AU and Foord EE. 1999. Cooling rates and crystallization dynamics of shallow level pegmatite-aplite dikes, San Diego County, California. American Mineralogist, 84(5-6): 708-717
[81] Wei ZL. 1994. Oil-gas pools in transition spans and analysis of their formation conditions. Geotectonica et Metallogenia, 18(1): 26-32(in Chinese with English abstract)
[82] Weinberg RF and Podladchikov YY. 1994. Diapiric ascent of magmas through power law crust and mantle. J. Geophys. Res., 99(B5): 9543-9559
[83] Wu CR. 1992. Magmatic thermo-metamorphism of coal. Coal Geology and Exploration, 20(4): 24-29(in Chinese with English abstract)
[84] Yang KG and Yang WR. 1997. Post-collision orogeny process and origin of huge quantity granites in orogenic belt. Geological Science and Technology Information, 16(4): 16-22(in Chinese with English abstract)
[85] Yang XK, Liu CY, Yang YH et al. 2005. The concept, classification and research progress on thermal structure. Earth Science Frontiers, 12(4): 385-396(in Chinese with English abstract)
[86] Zhang BT, Wu JQ, Ling HF et al. 2007. The calculation of U-Th-K radiogenic heat and effects on granitic cooling-crystallization process, and its geological significance. Science in China (Series D), 37(2): 155-159 (in Chinese)
[87] Zhang J and Shi YL. 1997. The thermal modeling of magma intrusion in sedimentary basins. Progress in Geophysics, 12(3): 55-64(in Chinese with English abstract)
[88] Zhang Q. 2012. Granite and metallogenic relationship discussion (1): Comment on the popular magmatic hydrothermal mineralization theory. Gansu Geology, 21(4): 1-14(in Chinese with English abstract)
[89] Zhang Q and Li CD. 2012. Granites: Implications for Continental Geodynamics. Beijing: Ocean Press, 1-276 (in Chinese with English abstract)
[90] Zhang Q, Jin WJ, Li CD and Jiao ST. 2013. Identification and implication of magma thermal field in the geothermal field. Progress in Geophys., 28(5): 2495-2507 (in Chinese with English abstract)
[91] Zhang YH and Gu JY. 2003. Hydrothermal circular flow: Significant influence on reservoirs in intrusion rock-exomorphic zone. Special Oil and Gas Reservoirs, 10(1): 86-89(in Chinese with English abstract)
[92] Zhao CB, Lin G, Hobbs BE, Ord A, Wang YJ and Mühlhaus HB. 2003. Effects of hot intrusions on pore-uid ow and heat transfer in uid-saturated rocks. Comput. Methods Appl. Mech. Engrg., 192(16-18): 2007-2030
[93] Zhao CP, Ran H and Wang Y. 2012. Present-day mantle-derived helium release in the Tengchong volcanic field, Southwest China: Implications for tectonics and magmatism. Acta Petrologica Sinica, 28(4): 1189-1204(in Chinese with English abstract)
[94] Zhou AC, Zhao YS, Guo JJ et al. 2010. Study of geothermal extraction scheme of hot dry rock in Tibetan Yangbajing region. Chinese Journal of Rock Mechanics and Engineering, 29(Suppl. 2): 89-95(in Chinese with English abstract)
[95] Zhou JY, Wu ZL and Zhuang XG. 1997. Geothermal field study and its significance in eastern Zhejiang, Fujian and Guangdong provinces. Geological Science and Technology Information, 16(2): 7-12(in Chinese with English abstract)
[96] Zhou TF, Yue SC and Lan TY. 1995. Magmatic dynamics and its relations to metallization of diorites in Yueshan, Anhui Province. Mineral Deposits, 14(4): 303-313(in Chinese with English abstract)
[97] Zhou XR, Wu CL, Huang XC et al. 1993. Characteristics of cognate inclusions in intermediate-acid intrusive rocks of Tongling area and their magmatic dynamcs. Acta Petrologica et Mineralogica, 12(1): 20-31(in Chinese with English abstract)
[98] Zhu CQ, Tian YT, Xu M et al. 2010a. The effect of Emeishan supper mantle plume to the thermal evolution of source rocks in the Sichuan basin. Chinese J. Geophys., 53(1): 119-127(in Chinese with English abstract)
[99] Zhu CQ, Xu M, Yuan YS et al. 2010b. Palaeo-geothermal response and record of the effusing of Emeishan basalts in Sichuan basin. Chinese Sci. Bull., 55(10): 949-956
[100] 陈衍景, Pirajno F, 赖勇等. 2004. 胶东矿集区大规模成矿时间和构造环境. 岩石学报, 20(4): 907-922
[101] 从柏林. 1978. 岩浆活动与火成岩组合. 北京: 地质出版社
[102] 冯乔, 汤锡元. 1997. 岩浆活动与油气成藏地质条件的关系. 西北地质科学, 18(1): 56-62
[103] 冯佐海. 2003. 广西姑婆山-花山花岗岩体侵位过程及构造解析. 博士学位论文.长沙: 中南大学
[104] 高贵荣, 张勉斌. 1998. 江西省于都县银坑贵多金属矿田矿化特征、成矿控制及找矿方向. 火山地质与矿产, 19(4): 347-356
[105] 郭进京, 周安朝. 2008. 高温岩体圈定的思路与方法探讨. 能源工程, (6): 1-4
[106] 金旭东, 张德会, 万天丰. 2010. 隐伏岩体顶上带与深部成矿预测. 地质通报, 29(2-3): 392-400
[107] 李德东, 罗照华, 周久龙, 杨宗锋, 刘翠. 2011. 岩墙厚度对成矿作用的约束: 以石湖金矿为例. 地学前缘, 18(1): 166-178
[108] 刘勇胜, 高山. 1998. 地壳深熔(anatexis)与花岗岩对下地壳的示踪作用. 地质科技情报, 17(3): 31-38
[109] 罗文积, 陈家清. 1997. 双向汇聚热液成矿. 甘肃地质学报, 6(增刊): 44-49
[110] 罗照华, 黄忠敏, 柯珊. 2007. 花岗质岩石的基本问题. 地质论评, 53(增刊): 180-226
[111] 罗照华, 卢欣祥, 陈必河. 2009. 透岩浆流体成矿作用导论. 北京: 地质出版社, 1-168
[112] 罗照华, 刘嘉麒, 赵慈平, 郭正府, 程黎鹿, 李晓惠, 李大鹏. 2011. 深部流体与岩浆活动: 兼论腾冲火山群的深部过程. 岩石学报, 27(10): 2855-2862
[113] 吕新彪, 姚书振, 周宗桂. 1997. 江西九瑞地区岩浆岩系统与成矿. 华南地质与矿产, (1): 27-36
[114] 马昌前. 1986. 岩浆活动中某些动力学参数的估算方法: 流体动力学原理的应用. 地质科技情报, 5(3): 47-54
[115] 马昌前, 杨坤光, 唐仲华等. 1994. 花岗岩类岩浆动力学--理论方法及鄂东花岗岩类例析. 武汉: 中国地质大学出版社
[116] 马东升. 1998. 地壳中流体的大规模流动系统及其成矿意义. 高校地质学报, 4(3): 250-261
[117] 毛景文, 谢桂青, 张作衡等. 2005. 中国北方中生代大规模成矿作用的期次及其地球动力学背景. 岩石学报, 21: 169-188
[118] 毛景文, 谢桂青, 郭春丽等. 2007. 南岭地区大规模钨锡多金属成矿作用: 成矿时限及地球动力学背景. 岩石学报, 23(10): 2329-2338
[119] 切列缅斯基ГА. 1977. 实用地热学.见:赵羿, 陈明译. 1982.北京: 地质出版社
[120] 饶超. 2007. 热液成矿过程中热及流体传输的计算机数值模拟. 硕士学位论文, 成都: 成都理工大学
[121] 邵飞. 2007. 水-岩相互作用及其与铀矿成矿关系研究--以相山铀矿田为例. 博学位论文, 武汉: 中国地质大学
[122] 唐晓音, 张功成, 梁建设等. 2013. 琼东南盆地长昌凹陷火成岩侵入体对温度场及烃源岩成熟度的影响. 地球物理学报, 56(1): 159-169
[123] 万天丰, Teyssier C, 曾华霖等. 2000. 山东玲珑花岗质岩体侵位机制. 中国科学(D辑), 30(4): 337-340
[124] 万志军, 赵阳升, 康建荣. 2005. 高温岩体地热资源模拟与预测方法. 岩石力学与工程学报, 24(6): 945-949
[125] 王大勇, 陆现彩, 徐士进等. 2011. 沉积盆地内侵入岩席对富含有机质围岩热影响的热传输模型研究. 南京大学学报(自然科学版), 47(1): 45-50
[126] 王满, 王英伟, 薛林福. 2012. 岩浆侵入对围岩中有机质成熟度影响的有限元模拟. 断块油气田, 19(2): 172-176
[127] 王涛, 王晓霞, 李伍平. 1999a. 试论花岗质深成岩体的复合定位机制及定位空间问题. 地质论评, 45(2): 142-150
[128] 王涛, 张国伟, 王晓霞等. 1999b. 花岗岩体生长方式及构造运动学、动力学意义. 地质科学, 34(3): 326-335
[129] 魏洲龄. 1994. 过渡际油气藏及其形成条件分析. 大地构造与成矿学, 18(1): 26-32
[130] 吴传荣. 1992. 煤的岩浆热变质作用. 煤田地质与勘探, 20(4): 24-29
[131] 杨坤光, 杨巍然. 1997. 碰撞后的造山过程及造山带巨量花岗岩的成因. 地质科技情报, 16(4): 16-22
[132] 杨兴科, 刘池洋, 杨永恒等. 2005. 热力构造的概念分类特征及其研究进展. 地学前缘, 12(4): 385-396
[133] 章邦桐, 吴俊奇, 凌洪飞等. 2007. U-Th-K放射成因热对花岗岩冷却-结晶过程影响的计算及地质意义. 中国科学(D辑), 37(2): 155-159
[134] 张健, 石耀霖. 1997. 沉积盆地岩浆侵入的热模拟. 地球物理学进展, 12(3): 55-64
[135] 张旗. 2012. 花岗岩与成矿关系的讨论之一: 评流行的岩浆热液成矿理论. 甘肃地质, 21(4): 1-14
[136] 张旗, 李承东. 2012. 花岗岩: 地球动力学意义. 北京: 海洋出版社, 1-276
[137] 张旗, 金惟俊, 李承东, 焦守涛. 2013. 地热场中"岩浆热场"的识别及其意义. 地球物理学进展, 28(5): 2495-2507
[138] 张映红, 顾家裕. 2003. 热液环流: 侵入岩-外变质带储层发育的重要影响因素. 特种油气藏, 10(1): 86-89
[139] 赵慈平, 冉华, 王云. 2012. 腾冲火山区的现代幔源氦释放: 构造和岩浆活动意义. 岩石学报, 28(4): 1189-1204
[140] 周安朝, 赵阳升, 郭进京等. 2010. 西藏羊八井地区高温岩体地热开采方案研究. 岩石力学与工程学报, 29(增刊2): 89-95
[141] 周江羽, 吴冲龙, 庄新国. 1997. 浙闽粤东部地热场研究及其意义. 地质科技情报, 16(2): 7-12
[142] 周涛发, 岳书仓, 兰天佑. 1995. 安徽月山地区闪长岩类岩浆动力学及其与成矿作用的联系. 矿床地质, 14(4): 303-313
[143] 周珣若, 吴才来, 黄许陈等. 1993. 铜陵中酸性侵入岩同源包体特征及岩浆动力学. 岩石矿物学杂志, 12(1): 20-31
[144] 朱传庆, 田云涛, 徐明等. 2010a. 峨眉山超级地幔柱对四川盆地烃源岩热演化的影响. 地球物理学报, 53(1): 119-127
[145] 朱传庆, 徐明, 袁玉松等. 2010b. 峨眉山玄武岩喷发在四川盆地的地热学响应. 科学通报, 55(6): 474-482