楊立中　劉金輝　孫占學(xué), 　王安東　萬建軍　周　毅

(1.東華理工大學(xué)水資源與環(huán)境工程學(xué)院；2.省部共建核資源與環(huán)境國家重點(diǎn)實(shí)驗(yàn)室培育基地)

?

漳州巖體放射性生熱率特征及干熱巖資源潛力*

楊立中1劉金輝1孫占學(xué)1,2王安東2萬建軍2周毅1

(1.東華理工大學(xué)水資源與環(huán)境工程學(xué)院；2.省部共建核資源與環(huán)境國家重點(diǎn)實(shí)驗(yàn)室培育基地)

摘要分析了漳州巖體地表鉆孔巖石樣品的密度以及鈾、釷、鉀含量，并進(jìn)行了放射性生熱率計(jì)算及生熱率特征分析。在此基礎(chǔ)上，對(duì)漳州巖體放射性生熱率進(jìn)行比較、干熱巖地?zé)岬刭|(zhì)特征分析和干熱巖資源量估算，結(jié)果表明：漳州巖體具有較高的放射性生熱率、良好的干熱巖地?zé)岬刭|(zhì)背景、較大的干熱巖資源儲(chǔ)量，總熱能為3.244×10(21)J，折合標(biāo)準(zhǔn)煤約110.71億t，提取2%折合2.21億t標(biāo)準(zhǔn)煤，提取20%折合22.14億t標(biāo)準(zhǔn)煤，提取40%折合44.28億t標(biāo)準(zhǔn)煤，且主要以高溫形式存在，占總資源量的91.37%，干熱巖資源潛力較大。上述研究成果可為福建地區(qū)干熱巖資源勘查評(píng)價(jià)及開發(fā)利用提供依據(jù)。

關(guān)鍵詞放射性生熱率干熱巖資源潛力儲(chǔ)量估算資源勘查評(píng)價(jià)

在地殼巖石中含有多種放射性元素，該類放射性元素衰變釋放出的熱能是地球內(nèi)熱的主要熱源之一[1]。對(duì)地球內(nèi)熱有明顯貢獻(xiàn)的放射性元素須具備足夠豐度、生熱量大、半衰期長等3個(gè)條件，U、Th、K等放射性元素具備上述條件，屬于放射性生熱元素。我國東南沿海地區(qū)花崗巖類分布廣泛，大面積出露，并與多種金屬礦床(如鈾礦床等)和地?zé)岙惓?如溫泉等)關(guān)系密切[2-3]，表明該區(qū)花崗巖體有重要的放射性地?zé)岬厍蚧瘜W(xué)研究價(jià)值。我國東南沿海地區(qū)花崗巖放射性生熱率分布以粵中—贛南為核心，向外擴(kuò)展呈下降趨勢，其中福建、江西花崗巖放射性生熱率均值為3.3 μW/m3，浙江為3.4 μW/m3，安徽為2.1 μW/m3，廣東高達(dá)5.7 μW/m3，可見該區(qū)花崗巖的放射性生熱率普遍較高[4]?？傮w來講，我國東南沿海地區(qū)花崗巖放射性生熱率研究程度較高，但就福建省而言，僅有89個(gè)巖石放射性生熱率數(shù)據(jù)(其中花崗巖55個(gè)、非花崗巖34個(gè)，花崗巖放射性生熱率數(shù)據(jù)中鉆孔數(shù)據(jù)38個(gè)，地表露頭數(shù)17個(gè)[5-6])，可見福建巖石放射性生熱率數(shù)據(jù)資料較少，無法滿足干熱巖勘查工作要求。為進(jìn)一步完善福建地區(qū)巖石放射性生熱率資料，本研究對(duì)漳州巖體采集16塊具代表性的地表露頭樣品和38塊鉆孔巖芯樣品進(jìn)行放射性生熱率特征研究，為該區(qū)干熱巖勘查提供參考。

1漳州巖體地質(zhì)概況與樣品采集

1.1漳州巖體地質(zhì)概況

漳州巖體位于福建東南沿海漳州—長泰一帶，出露面積達(dá)900 km2，為由燕山晚期輝長蘇長巖、石英二長閃長巖、花崗閃長巖、更長環(huán)斑二長花崗巖、細(xì)?；◢弾r、晶洞堿長花崗巖和堿長花斑巖等7次侵入體組成的復(fù)式巖基，巖體位于福安—南靖和仙游—云霄NE向斷裂與廈門—漳州EW向斷裂交匯處北側(cè)，同時(shí)受湯川—鐘魏NE向隱伏背斜控制，為一等軸狀巖基(圖1)，圍巖主要由上侏羅統(tǒng)南園組英安質(zhì)-流紋質(zhì)晶屑凝灰熔巖和晶屑凝灰?guī)r、下侏羅統(tǒng)梨山組和上三疊統(tǒng)文賓山組含長石石英砂巖、粉砂巖等組成[7]。

1.2漳州巖體地?zé)岙惓?/p>

巖體地?zé)岙惓Ｖ饕詿崛?、溫泉形式顯示，漳州地?zé)崽锸悄壳拔覈鴸|南沿海地區(qū)(除臺(tái)灣地區(qū)以外)已發(fā)現(xiàn)的中—低溫地?zé)崽镏袦囟茸罡叩?，鉆孔揭露的熱水溫度：孔口106 ℃，孔底90 m 121.5 ℃[8]?？梢姡闹輲r體具有良好的地?zé)岜尘啊?/p>

1.3巖體樣品采集

圖1　漳州巖體構(gòu)造位置[7]

圖2　漳州巖體取樣點(diǎn)布置

圖3　漳州巖體

2漳州巖體放射性生熱率特征

2.1密度

將采集的樣品進(jìn)行適當(dāng)處理后，用電子天平蠟封法測量所有樣品的密度(誤差±5%以內(nèi))，結(jié)果表明：地表樣品密度2.36～2.84 g/cm3，鉆孔樣品密度2.18～2.78 g/cm3，平均2.60 g/cm3。

2.2U、Th、K含量

U、Th、K含量均在廣州澳實(shí)測試中心完成，U、 Th含量測定采用ICP-MS溶液法，數(shù)據(jù)精度在 ±10%以內(nèi)，檢測下限4×10-6，上限1 000×10-6；K2O含量采用X射線熒光光譜法(XRF)測定，檢測下限0.01%，上限50%，測試精度誤差在±5%以內(nèi)。測試結(jié)果均在檢查范圍內(nèi)，無異常值，見表1。由表1可知，w(U)平均7.57×10-6，w(Th)平均30.41×10-6，w(K)平均4.16%，漳州巖體U、Th、K含量分布具不均一性。

表1　漳州巖體U、Th、K含量

2.3放射性生熱率計(jì)算

放射性生熱率是指單位體積巖石中所含放射性元素在單位時(shí)間內(nèi)經(jīng)放射性衰變所產(chǎn)生的能量。本研究采用Rybach 根據(jù)修正后的天然放射性核素計(jì)算公式進(jìn)行計(jì)算：

(1)

式中，A為巖石放射性生熱率，μW/m3；ρ為巖石密度，g/cm3；CU、CTh分別為巖石中U、Th含量，(×10-6)；CK為巖石中K含量，%。

表2　漳州巖體放射性生熱率

2.4放射性生熱率分布特征

放射性生熱率在花崗巖體中的分布主要由水平分布和垂向分布組成。區(qū)內(nèi)花崗巖體放射性生熱率水平分布呈多樣性，不存在統(tǒng)一規(guī)律(圖4)[9]。在垂向上，因巖石放射性生熱率以酸性巖類最高，超基性巖類最低，隨深度的增加巖石越來越偏基性，理論上講，放射性生熱率應(yīng)隨深度的增加而遞減[10]，但地殼淺部巖石的放射性生熱率并未隨著深度的增加而成指數(shù)衰減，多呈鋸齒狀變化(圖5、圖6)[9]。

圖4　漳州巖體放射性生熱率水平分布(單位：μW/m3)

圖孔放射性生熱率隨深度變化

圖孔放射性生熱率隨深度變化

2.5漳州巖體放射性熱流貢獻(xiàn)率

地殼巖石放射性熱流是影響地表熱流大小的因素之一，尤其是花崗巖體大面積分布的地區(qū)，地表熱流中的放射性熱流占相當(dāng)大的比例。中國東南地區(qū)放射性熱流占地表熱流的40%以上[11]，若漳州地區(qū)背景熱流(73.2 mW/m2)的40%由放射性生熱提供,則需厚5 km厚的花崗巖體，若全部由放射性熱流提供，則需13 km厚的花崗巖體[5]，福建地區(qū)由地殼放射性生熱產(chǎn)生的平均熱流值為71.81 mW/m2，占地表熱流的51.17%[12]。

漳州地區(qū)深成巖體一般厚7±2 km，最厚不超過12.0 km，區(qū)內(nèi)地殼分為上地殼和下地殼，上地殼厚16.5～18.8 km，下地殼厚12.0～13.0 km，上地殼分為上下2部分，在上地殼下部有1個(gè)低速層，速度約6.0 km/s，低速層頂部深約12.0 km，厚約5.0 km[13]，總體上，漳州巖體厚度在12 km以內(nèi)。本研究漳州巖體實(shí)測生熱率為4.22 μW/m3，對(duì)于厚度不超過5 km的巖體，放射性生熱率取值為4.22 μW/m3；對(duì)于厚5～12 km的巖體，放射性生熱率取值為3.11 μW/m3，背景異常熱流值取95.0 mW/m2，漳州巖體放射性熱流貢獻(xiàn)率計(jì)算結(jié)果見表3。由表3可知：漳州巖體放射性熱流和放射性熱流貢獻(xiàn)率隨厚度的增大而增大，放射性熱流平均值為26.71 mW/m2，放射性熱流貢獻(xiàn)率平均值為28.12%，可見，漳州巖體放射性熱流是該區(qū)地表熱流異常值的重要組成部分。

表3　漳州巖體放射性熱流及熱流貢獻(xiàn)率

3漳州巖體干熱巖資源潛力評(píng)價(jià)

干熱巖是溫度一般大于200 ℃、埋深數(shù)千米、內(nèi)部不存在流體或僅有少量地下流體的高溫巖體。該定義是基于商業(yè)開發(fā)目的得出的，即埋深在鉆探范圍內(nèi)溫度大于200 ℃的高溫巖體具有商業(yè)開發(fā)價(jià)值。我國干熱巖資源開發(fā)利用仍處于起步階段，國內(nèi)尚無干熱巖試驗(yàn)電站，從我國實(shí)際情況出發(fā)，可將4 km深度內(nèi)且?guī)r體溫度不低于150 ℃的地區(qū)作為干熱巖發(fā)電試驗(yàn)的重點(diǎn)靶區(qū)和干熱巖資源遠(yuǎn)景區(qū)。我國大陸3～10 km深度段干熱巖地?zé)豳Y源總量相當(dāng)于715～860萬億t標(biāo)準(zhǔn)煤，按2%的可開采資源量計(jì)算，為我國大陸2010年能源消耗總量的 4 400～5 200倍[14-16]，其中福建是我國干熱巖資源較豐富的省份之一，據(jù)初步估算，10 km以內(nèi)深度的資源儲(chǔ)量超過50萬億t標(biāo)準(zhǔn)煤，而漳州市干熱巖儲(chǔ)量占福建省的40%，相當(dāng)于20萬億t標(biāo)準(zhǔn)煤[17]。表明我國干熱巖資源總體豐富，福建省漳州市干熱巖開發(fā)前景較好。

3.1巖體放射性產(chǎn)熱能力評(píng)價(jià)

為綜合評(píng)價(jià)漳州巖體產(chǎn)熱能力，將其放射性生熱率與國內(nèi)外其他花崗巖體進(jìn)行對(duì)比，部分花崗巖體放射性生熱率見表4。由表4可知：巖體平均放射性生熱率為4.13 μW/m3。本研究得出的漳州巖體放射性生熱率為4.22 μW/m3，盡管目前尚無具體根據(jù)花崗巖體放射性生熱率判斷巖體產(chǎn)熱高低的具體標(biāo)準(zhǔn)，但從比較結(jié)果看，漳州巖體應(yīng)屬中等偏高的放射性產(chǎn)熱巖體。

表4部分花崗巖體放射性生熱率[6]

μW/m3

3.2巖體干熱巖地?zé)岬刭|(zhì)特征

通常將1個(gè)潛在的適合干熱巖開發(fā)的區(qū)域劃分為2個(gè)部分，即埋藏于地下一定深度的變質(zhì)巖或結(jié)晶巖體(基底)以及上覆于基底之上的蓋層?？梢姡挥谢缀蜕w層具備一定條件時(shí)才適合作為干熱巖體進(jìn)行開發(fā)。目前，世界上關(guān)于干熱巖體的評(píng)價(jià)尚無具體的地?zé)岬刭|(zhì)學(xué)標(biāo)準(zhǔn)，但從理論上講，具備干熱巖體開發(fā)的巖體應(yīng)具有一些必要的地?zé)岬刭|(zhì)學(xué)特征。為此，分別從地質(zhì)、地球物理、地?zé)岬刭|(zhì)等角度總結(jié)了漳州巖體的地?zé)岬刭|(zhì)學(xué)特征，見表5。由表5可知：漳州巖體基本具備了良好的干熱巖地?zé)岬刭|(zhì)特征。

3.3漳州巖體干熱巖資源量估算

干熱巖資源量(Q)主要采取體積法進(jìn)行估算，公式[15-16]如下：

(2)

式中，ρ為巖石密度，g/cm3；Cp為巖石比熱容，J/(kg·℃)；V為巖體體積，m3；T(z)為計(jì)算深度的巖體溫度，℃；z為巖體深度，m；T0為地表溫度，℃。

在穩(wěn)態(tài)熱流狀態(tài)下的深部溫度計(jì)算公式[15-16]為：

(3)

式中，q0為地表熱流值，mW/m2；K為巖石熱導(dǎo)率，w/(m·K)。

本研究漳州巖體各參數(shù)取值見表6。

表5　漳州巖體干熱巖地?zé)岬刭|(zhì)學(xué)指標(biāo)[5,7,18-21]

表6　漳州巖體干熱巖資源量估算參數(shù)[6-7,10,14-15]

注：D為放射性元素富集層厚度，S為巖體面積，T0取漳州地區(qū)多年平均溫度。

本研究干熱巖的評(píng)價(jià)深度為3～10 km[15]，據(jù)表6及式(2)、式(3)，分別計(jì)算了漳州巖體3～10 km深度內(nèi)各段的干熱巖資源潛力，結(jié)果見表7。由表7可知：漳州巖體3～10 km深處干熱巖熱能總量為3.244×1021J，折合標(biāo)準(zhǔn)煤約110.71億t，干熱巖能量提取上限為40%，按下限2%提取折合2.21億t 標(biāo)準(zhǔn)煤，按中值提取20%折合22.14億t標(biāo)準(zhǔn)煤，按上限提取40%折合44.28億t標(biāo)準(zhǔn)煤，其中≥150 ℃的干熱巖熱能為2.964×1021J，約占熱能總量的91.37%，可見，漳州巖體干熱巖資源以高溫資源為主，且干熱巖資源總量超過100.0億t標(biāo)準(zhǔn)煤。為此，該區(qū)不僅有利于干熱巖開發(fā)，而且應(yīng)成為我國深層地?zé)豳Y源開發(fā)的首選地區(qū)[22-23]。從放射性產(chǎn)熱能力、放射性熱流貢獻(xiàn)率、干熱巖地?zé)岬刭|(zhì)背景和干熱巖資源總量等方面看，漳州巖體干熱巖開發(fā)潛力較大。

表7　漳州巖體3～10 km深處干熱巖資源估算結(jié)果

4結(jié)語

通過對(duì)漳州地區(qū)巖體進(jìn)行系統(tǒng)的地表和鉆孔取樣，并對(duì)樣品進(jìn)行密度以及U、Th、K含量測試，經(jīng)相關(guān)計(jì)算獲得了漳州巖體樣品放射性生熱率值，據(jù)此對(duì)巖體分別從放射性生熱率特征、放射性熱流貢獻(xiàn)率、放射性產(chǎn)熱能力、干熱巖資源潛力等方面進(jìn)行了詳細(xì)分析。研究成果可在一定程度上完善福建地區(qū)巖石放射性生熱率數(shù)據(jù)，對(duì)于該區(qū)干熱巖資源勘查及開發(fā)有一定的參考價(jià)值。

參考文獻(xiàn)

[1]汪集旸，孫占學(xué).神奇的地?zé)醄M].北京：清華大學(xué)出版社，2001.

[2]張祖還，章邦桐.華南產(chǎn)鈾花崗巖及有關(guān)鈾礦床研究[M].北京：原子能出版社，1991.

[3]胡郁樂，張惠.深部地?zé)徙@井與成井技術(shù)[M].北京：中國地質(zhì)大學(xué)出版社，2013.

[4]趙平，汪集旸，汪緝安，等.中國東南地區(qū)巖石生熱率分布特征[J].巖石學(xué)報(bào)，1995，11(3)：292-305.

[5]龐忠和.漳州盆地地?zé)嵯到y(tǒng)成因模式、熱能潛力與熱水分布規(guī)律[D].北京：中國科學(xué)院地質(zhì)與地球物理研究所，1987.

[6]趙平.中國東南地區(qū)巖石生熱率研究[D].北京：中國科學(xué)院地質(zhì)與地球物理研究所，1993.

[7]周珣若，吳克隆.漳州I-A型花崗巖[M].北京：科學(xué)出版社，1994.

[8]熊亮萍，汪集旸，龐忠和.漳州熱田的對(duì)流熱流和傳導(dǎo)熱流的研究[J].地球物理學(xué)報(bào)，1990(6)：702-711.

[9]趙平，汪集旸，汪緝安.熱流和巖石生熱率關(guān)系的研究[J].地質(zhì)科學(xué)，1996，31(3)：297-305.

[10]卡普邁耶，海涅爾.地?zé)釋W(xué)及其應(yīng)用[M].北京：科學(xué)出版社，1981.

[11]胡圣標(biāo)，汪集旸.中國東南地區(qū)地殼生熱率與地慢熱流[J].中國科學(xué)，1994，24(2)：185-193.

[12]周江羽，吳沖龍，莊新國.浙閩粵東部地?zé)釄鲅芯考捌湟饬x[J].地質(zhì)科技情報(bào)，1997，16(2)：7-12.

[13]朱金芳，方盛明，張先康，等.漳州盆地及其鄰區(qū)地殼深部結(jié)構(gòu)的探測與研究[J].中國地震，2006，22(4):405-417.

[14]袁玉松，馬永生，胡圣標(biāo)，等.中國南方現(xiàn)今地?zé)崽卣鱗J].地球物理學(xué)報(bào)，2006，49(4)：1118-1126.

[15]汪集旸，胡圣標(biāo)，龐忠和，等.中國大陸干熱巖地?zé)豳Y源潛力評(píng)估[J].科技導(dǎo)報(bào)，2012，30(32)：25-31.

[16]藺文靜，劉志明，馬峰，等.我國陸區(qū)干熱巖資源潛力估算[J].地質(zhì)學(xué)報(bào)，2012(5)：807-811.

[17]莊慶祥.福建省干熱巖研究第I期工程項(xiàng)目中地?zé)岚l(fā)電內(nèi)容的探討[J].能源與環(huán)境，2012(2)：53-54.

[18]熊紹柏，金東敏，孫克忠，等.福建漳州地?zé)崽锛捌溧徑貐^(qū)的地殼深部構(gòu)造特征[J].地球物理學(xué)報(bào)，1991，34(1)：55-63.

[19]林錦華.漳州斷陷盆地及發(fā)震構(gòu)造[J].地殼形變與地震，2001，21(1)：85-89.

[20]陳宗基.福建地?zé)岬厍蛭锢硌芯縖J].北京：中國科學(xué)技術(shù)出版社，1992.

[21]萬天豐，褚明記，陳明佑.福建省巖石圈的熱狀態(tài)與地?zé)豳Y源的遠(yuǎn)景評(píng)價(jià)[J].地質(zhì)學(xué)報(bào)，1988(2)：178-198.

[22]廖志杰，萬天豐，張振國.增強(qiáng)型地?zé)嵯到y(tǒng)：潛力大、開發(fā)難[J].地學(xué)前緣，2015，22(1)：335-343.

[23]黃少鵬.中國地?zé)崮茉撮_發(fā)的機(jī)遇與挑戰(zhàn)[J].地?zé)崮埽?015(1)：11-16.

Study of the Characteristics of Radioactive Heat Production Rate and Hot Dry Rock Resources Potential in Zhangzhou City

Yang Lizhong1Liu Jinhui1Sun Zhanxue1,2Wang Andong2Wan Jianjun2Zhou Yi1

(1.School of Water Resources & Environmental Engineering,East China Institute of Technology;2.State Key Laboratory Breeding Based of Nuclear Resources and Environment)

AbstractThe contents of U, Th, K of the rock samples collected in Zhangzhou city are tested and analyzed, besides that, the radioactive heat production rate is computed and the radioactive heat production rate characteristics is discussed.Based on the above analysis results, the radioactive heat production rate of the rock samples is conducted contrastive analysis, the geological characteristics of the hot dry rock resources is analyzed and the hot dry rock resources is estimated. The results show that the rock mass in Zhangzhou city with high radioactive heat production rate, good geothermal geological background and large reserves of hot dry rock resources,the total heat energy of the hot dry rock resources in Zhangzhou city is 3.244×10(21)J, which is equivalent to the heat energy provided by 11.071 billion ton standard coal, the 2% of the total heat energy of the hot dry rock resources in Zhanghzou city is equivalent to the heat energy provided by 221 million ton standard coal, the 20% of the total heat energy of the hot dry rock resources in Zhanghzou city is equivalent to the heat energy provided by 2.214 billion ton standard coal, the 40% of the total heat energy of the hot dry rock resources in Zhanghzou city is equivalent to the heat energy provided by 4.428 billion ton standard coal,the heat energy of the hot dry rock resources in Zhangzhou city is existed mainly in the form of high temperature, which is the 91.37% of the total resources reserves, the hot dry rock resources potential of the area is large. The above research result can provide reference for the exploration,evaluation, development and utilization of the hot dry rock resources in this area.

KeywordsRadioactive heat production rate, Hot dry rock, Resources potential, Reserves estimation

(收稿日期2016-01-03)

*中國地質(zhì)調(diào)查局計(jì)劃項(xiàng)目(編號(hào)：[2011]01-17-31)。

楊立中(1989—)，男，碩士研究生，330013 江西省南昌市昌北經(jīng)濟(jì)開發(fā)區(qū)廣蘭大道418號(hào)。