原位溶浸采礦理論研究

梁衛(wèi)國(guó)，趙陽(yáng)升，徐素國(guó)，于艷梅

（太原理工大學(xué) 采礦研究所，太原 030024)

自從有人類(lèi)歷史以來(lái)，就有采石采礦活動(dòng)記錄。早在舊石器時(shí)代，人類(lèi)就有了采石與制作石器的技術(shù)，且主要為露天開(kāi)采；到新石器時(shí)代，采礦技術(shù)即有重大進(jìn)步，在露天開(kāi)采的基礎(chǔ)上，開(kāi)始了地下開(kāi)采。自此，以露天開(kāi)采與地下井工開(kāi)采為主的人類(lèi)采礦技術(shù)即初步形成，并一直延續(xù)至今。溶浸采礦是利用礦物自身的物理化學(xué)特性，將水、化學(xué)溶劑或微生物，注入礦床或礦堆，進(jìn)行有選擇性地溶解、浸出和回收其中有用礦物組份。鹽類(lèi)礦床水溶開(kāi)采即是溶浸采礦的典型代表，它利用鹽類(lèi)礦床易溶于水的特性，將淡水注入鹽巖礦床，溶解鹽巖后鹵水采汲至地表，再對(duì)其加熱蒸發(fā)結(jié)晶即獲得鹽巖晶體顆粒。利用水溶開(kāi)采的世界井礦鹽史，即源于中國(guó)且已有2000多年的歷史。在金屬礦（鈾、銅、金、銀等)開(kāi)采中，溶浸采礦是普遍應(yīng)用的一項(xiàng)技術(shù)，有地表堆浸、就地破碎溶浸和原地鉆孔溶浸采礦三種方法［1-3］。但其發(fā)展歷史遠(yuǎn)晚于鹽類(lèi)礦床的水溶開(kāi)采，從20世紀(jì)70年代才開(kāi)始在國(guó)內(nèi)外興起。在經(jīng)歷了數(shù)千年的露天開(kāi)采與地下井工開(kāi)采之后，人類(lèi)采礦技術(shù)的發(fā)展方向應(yīng)當(dāng)是原位溶浸開(kāi)采或固體礦物流體化開(kāi)采。

溶浸采礦技術(shù)在我國(guó)歷史悠久。20世紀(jì)60年代，安徽銅陵有色金屬公司松樹(shù)山銅礦率先應(yīng)用就地破碎浸礦法回收銅殘礦；90年代中后期，德興銅礦建成年產(chǎn)2000t電銅堆浸廠；2000年，中條山銅礦峪礦建成年產(chǎn)500t電解銅的地下溶浸提銅示范系統(tǒng)；2003年，云南官房銅礦建成處理含銅0.9%的原生硫化銅和次生硫化銅的生物堆浸廠。在國(guó)外，該技術(shù)已成功用于多種金屬的提取，目前已有幾十座銅、鈾、金礦山有大規(guī)模的工業(yè)應(yīng)用。早在16世紀(jì)，德國(guó)Harz地區(qū)和西班牙Rio Tinto礦山就采用溶浸方法處理硫化礦石。1947年Colmer和Hinkle又從酸性廢坑水中分離出氧化亞鐵硫桿菌，人們開(kāi)始逐漸認(rèn)識(shí)和利用細(xì)菌在硫化礦浸出中的作用。20世紀(jì)60年代，美國(guó)和蘇聯(lián)開(kāi)始試用原地溶浸法開(kāi)采鈾礦，并在懷俄明州建成世界第一座溶浸采鈾礦山，20世紀(jì)70年代，美國(guó)Kennecott礦業(yè)公司建成世界最大的微生物堆浸場(chǎng)，處理礦石量36億t，年產(chǎn)7.2萬(wàn)噸銅金屬。1986年南非Fairview礦建成世界第一個(gè)金礦生物浸出工廠，日處理金精礦10t。近20年，溶浸技術(shù)在美國(guó)、智利、澳大利亞等國(guó)得到推廣。

原位溶浸開(kāi)采是泛指通過(guò)地面鉆井，將溶浸、熱解流體（含高溫流體、微生物流體、水及空氣)注入固體礦層，通過(guò)注入流體與礦層有用組分發(fā)生物理、化學(xué)、熱、力耦合作用，將礦層中有用組分轉(zhuǎn)變?yōu)榱黧w或附著成分解吸游離，通過(guò)滲流方式傳輸至鉆井而后排采至地面的一類(lèi)采礦方法，如圖1所示。鹽類(lèi)礦床的水溶開(kāi)采，鈾礦、銅礦、金礦等的化學(xué)溶浸開(kāi)采，煤的地下原位熱解開(kāi)采，通過(guò)驅(qū)替技術(shù)而實(shí)施的煤層氣開(kāi)采，油頁(yè)巖油氣原位注熱開(kāi)采等等許多地下礦物的開(kāi)采，都屬于原位溶浸開(kāi)采之范疇。由于人員設(shè)備不必下入地下深部礦層，原位溶浸采礦具有安全、高效、環(huán)保的巨大技術(shù)優(yōu)勢(shì)，在地下礦產(chǎn)資源開(kāi)采、尤其是深部資源開(kāi)采中具有極為廣闊的發(fā)展及應(yīng)用前景。

圖1 原位溶浸采礦示意圖

自1999年以來(lái)，在對(duì)山西運(yùn)城鹽湖地下薄層芒硝礦、四川眉山難溶鈣芒硝礦開(kāi)采技術(shù)研究基礎(chǔ)上，太原理工大學(xué)對(duì)除金屬礦以外其他相關(guān)資源的原位溶浸開(kāi)采理論與技術(shù)進(jìn)行了深入細(xì)致研究，并取得一些初步研究成果。包括鹽類(lèi)礦床原位溶浸開(kāi)采、油頁(yè)巖原位熱解開(kāi)采油氣、煤的原位熱解、注氣驅(qū)替煤層瓦斯、高溫巖體地?zé)衢_(kāi)發(fā)等。借助先進(jìn)的顯微CT系統(tǒng)試驗(yàn)，揭示了溶浸開(kāi)采過(guò)程中這些礦物巖體的細(xì)觀孔裂隙結(jié)構(gòu)衍化規(guī)律特征，建立了礦物巖體溶解（熱解)——滲透互進(jìn)作用機(jī)理，發(fā)明了多項(xiàng)針對(duì)不同礦物的原位溶浸開(kāi)采新技術(shù)，極大推動(dòng)了原位溶浸采礦技術(shù)的進(jìn)步與發(fā)展［4-7］。本文將對(duì)上述相關(guān)礦物的原位溶浸開(kāi)采理論與技術(shù)進(jìn)行介紹分析，以期對(duì)原位溶浸采礦理論與技術(shù)的概念及內(nèi)涵進(jìn)一步拓展完善。

2 原位溶浸采礦基本原理

針對(duì)不同礦物特性及其地質(zhì)賦存特征，原位溶浸開(kāi)采技術(shù)方法存在一定差異。如鹽類(lèi)礦床的原位溶浸開(kāi)采，采用淡水作為基本溶劑，輔以一定壓力與溫度溶浸效果更佳；對(duì)油頁(yè)巖與煤體中油氣的原位開(kāi)采，則需要注入過(guò)熱高溫蒸汽，進(jìn)行能量轉(zhuǎn)化開(kāi)采其中油氣資源；對(duì)煤體中吸附的瓦斯氣體，則需要注入蒸汽或其他強(qiáng)吸附性氣體，對(duì)煤體中吸附瓦斯進(jìn)行解吸置換驅(qū)替，等等［8-11］。但是，所有這些技術(shù)都具有一個(gè)共同特征：注入地下一種介質(zhì)，即可在原位將賦存于礦物中的有用組份進(jìn)行物理或化學(xué)性質(zhì)改變，或固態(tài)轉(zhuǎn)變?yōu)橐簯B(tài)，或吸附態(tài)轉(zhuǎn)變?yōu)橛坞x態(tài)，并通過(guò)一定滲流擴(kuò)散運(yùn)移通道流出地表。顯然，影響原位溶浸開(kāi)采技術(shù)的主要因素有兩個(gè)方面：其一，注入介質(zhì)與礦體中有用組份的物理化學(xué)作用效果，其二為注入介質(zhì)及溶浸熱解產(chǎn)物進(jìn)出礦體的通道滲透性。為此，在實(shí)驗(yàn)室對(duì)不同礦物的溶浸熱解特性以及溶浸熱解前后礦體孔裂隙細(xì)觀結(jié)構(gòu)衍化及滲透性改變進(jìn)行了系統(tǒng)試驗(yàn)。

2.1 鹽類(lèi)礦物的溶解特性

鹽類(lèi)物質(zhì)在水的溶解，是一個(gè)偶極水分子與鹽類(lèi)物質(zhì)晶體格架中的鹽離子相互吸引并破壞晶體格架的過(guò)程。溶解于水中的鹽離子以擴(kuò)散的方式在溶液中彌散運(yùn)移，而根據(jù)Fick擴(kuò)散定律，濃度差為擴(kuò)散的基本動(dòng)力。因此，溶液濃度、溫度以及流動(dòng)狀態(tài)決定著鹽類(lèi)物質(zhì)在水溶液中的溶解特性。

圖2 鹽巖溶解速度與溶液密度關(guān)系

通過(guò)大量試驗(yàn)研究發(fā)現(xiàn)：隨溶液濃度的增高，由于溶液中可溶鹽離子數(shù)量接近飽和，以及水分子對(duì)鹽離子吸引力的減弱，鹽類(lèi)物質(zhì)的溶解速度降低。圖2所示為芒硝鹽巖在兩個(gè)不同方向的溶解速度（V)與濃度（C)關(guān)系，兩個(gè)方向溶解速度與濃度均呈擬合度極高的線性關(guān)系，其中上向溶解V＝2.3205-0.072C（R2＝0.9987)；側(cè)向溶解V＝1.3105-0.0302C（R2＝0.9603)。當(dāng)溶液濃度接近飽和狀態(tài)時(shí)，上溶速度與側(cè)溶速度十分接近，溶解速度極低。

不同溫度淡水溶液中鈣芒硝溶解速率曲線如圖3所示。

圖3 鹽巖溶解速率與溶液溫度關(guān)系曲線

由圖3可以看出：鈣芒硝溶解速率由25℃的0.27g／cm2／h減小到0.23g／cm2／h附近，隨溫度升高，鈣芒硝溶解速率迅速增加，到40℃附近達(dá)到最大值0.36g／cm2／h左右。之后巖樣溶解速率迅速降低，到55℃，其值達(dá)到最小0.23g／cm2／h。溫度對(duì)鹽巖溶解促進(jìn)作用存在臨界點(diǎn)（～45℃)，當(dāng)溫度超過(guò)該臨界值后，鹽巖溶解速度不再隨溫度升高而增大。動(dòng)態(tài)溶解試驗(yàn)結(jié)果表明，鹽巖動(dòng)態(tài)溶解速率約為靜態(tài)條件下的1.2～1.5倍，且影響溶解速率的動(dòng)態(tài)流速也存在臨界點(diǎn)，當(dāng)超過(guò)該流速臨界值后，溶解速度不再隨流速的增大而升高。

上述溶解特性均為易溶鹽類(lèi)礦物試驗(yàn)結(jié)果。對(duì)鈣芒硝等難溶鹽類(lèi)礦物，其溶解機(jī)理與易溶鹽類(lèi)礦物相同，但由于其中難溶礦物與雜質(zhì)的存在，礦體的滲透特性影響決定著其溶解特性，且二者存在溶解滲透相互促進(jìn)的作用關(guān)系，2.3中細(xì)觀衍化結(jié)果對(duì)此進(jìn)行了很好說(shuō)明。影響鹽類(lèi)礦物溶解特性的主要因素有：鹽類(lèi)物質(zhì)自身物理化學(xué)特性、溶解環(huán)境（包括溶液溫度、濃度、壓力)以及溶解對(duì)礦物物理結(jié)構(gòu)（主要指孔隙率)等的影響。

2.2 煤與油頁(yè)巖的熱解特性

對(duì)煤、油頁(yè)巖這一類(lèi)有機(jī)變質(zhì)體或有機(jī)變質(zhì)賦存體，由于成煤、聚油過(guò)程中曾經(jīng)歷過(guò)溫度驟變作用，在高溫作用下，固體結(jié)構(gòu)發(fā)生改變的同時(shí)，其中飽含的油氣物質(zhì)的流動(dòng)性相應(yīng)增強(qiáng)，從而從煤及油頁(yè)巖中分離出來(lái)。對(duì)不同階煤加熱熱解過(guò)程中，還會(huì)有煤質(zhì)的變化作用，尤其低變質(zhì)煤會(huì)發(fā)生提質(zhì)改性，從而提高其經(jīng)濟(jì)價(jià)值。

太原理工大學(xué)在實(shí)驗(yàn)室內(nèi)對(duì)褐煤、長(zhǎng)焰煤、貧瘦煤、氣煤等煤種進(jìn)行了不同程度的熱解試驗(yàn)，發(fā)現(xiàn)所有煤種均存在影響煤體熱解及產(chǎn)氣速率的臨界溫度值，貧瘦煤加熱超過(guò)310℃時(shí)氣體產(chǎn)量達(dá)到最大，而長(zhǎng)焰煤在450～500℃為產(chǎn)氣最佳溫度段。對(duì)油頁(yè)巖進(jìn)行的熱解試驗(yàn)發(fā)現(xiàn)，與煤熱解相同，存在熱解臨界溫度值，在400～600℃之間，油頁(yè)巖熱解失重占總失重的80%。與鹽類(lèi)礦物溶解相類(lèi)似，影響煤與油頁(yè)巖熱解的因素包括：煤或油頁(yè)巖自身特性、熱解環(huán)境（溫度、壓力、升溫速率)以及熱解對(duì)礦體物理結(jié)構(gòu)（孔隙率)等的影響。

2.3 溶解、熱解前后礦體細(xì)觀結(jié)構(gòu)演化與滲透性

原位溶浸采礦過(guò)程中，礦體結(jié)構(gòu)及孔隙率衍化非常重要，它影響決定著原位溶浸采礦的效率與效益。對(duì)以鈣芒硝為代表的難溶鹽類(lèi)礦物溶浸試驗(yàn)研究發(fā)現(xiàn)，礦物溶解與礦體孔隙結(jié)構(gòu)衍化存在互進(jìn)作用。裂隙首先在固體軟弱膠結(jié)面處產(chǎn)生，溶液沿著生成裂隙進(jìn)入，由于裂隙中水壓的影響使得裂隙寬度發(fā)生改變，從而導(dǎo)致裂隙滲透規(guī)律的改變。隨著裂隙的擴(kuò)展溶液繼續(xù)沿著裂隙進(jìn)入，從而促進(jìn)礦物的溶解提高采礦的效率。對(duì)煤及油頁(yè)巖這類(lèi)富含雜質(zhì)的特殊資源而言，在低溫段由于內(nèi)部吸附氣體及水蒸氣的溢出以及固體顆粒間的熱脹冷縮作用，直接導(dǎo)致固體內(nèi)部產(chǎn)生豐富的孔、裂隙，提高了固體的滲透性，為高溫?zé)峤猱a(chǎn)生的有用氣體及油氣的傳輸提供通道。

如圖4、5分別為鈣芒硝在淡水溶浸不同時(shí)間細(xì)觀結(jié)構(gòu)衍化圖及孔隙率變化曲線。

圖4 淡水中鈣芒硝溶浸不同時(shí)間作用下細(xì)觀結(jié)構(gòu)衍化圖

圖5 淡水中溶浸12h鈣芒硝孔隙率變化曲線

由圖5可以看出，隨著溶蝕時(shí)間的增大，鈣芒硝孔隙率不斷增加。

圖6為瘦煤熱解過(guò)程中孔隙率隨溫度變化曲線。由圖6可以看出，在常溫到200℃，瘦煤孔隙率變化不大。隨著溫度的升高，孔隙率極速增加，到400℃達(dá)到最大。從400℃到600℃，瘦煤孔隙率逐漸減小。

圖6 瘦煤熱解過(guò)程中孔隙率隨溫度變化曲線

顯然隨溶浸熱解的進(jìn)行，礦物內(nèi)部孔裂隙結(jié)構(gòu)與孔隙率都在發(fā)生不斷的變化，這一變化直接影響著礦物溶解熱解的深入進(jìn)行。當(dāng)孔隙率增大時(shí)，有利于流體或氣體的滲入與流出，促進(jìn)礦物的溶解或熱解反應(yīng)。相反，則不利于原位溶解熱解的進(jìn)行。因此，在原位溶浸采礦技術(shù)中，礦體結(jié)構(gòu)衍化對(duì)溶浸采礦效果起著十分重要的作用。如何提高礦體的滲透性、增強(qiáng)溶浸液滲入與生產(chǎn)物流出是原位溶浸采礦的核心與關(guān)鍵。

3 原位溶浸采礦理論架構(gòu)

由文章引言部分對(duì)原位溶浸采礦過(guò)程描述可知，原位溶浸采礦是一個(gè)涉及流體運(yùn)移、礦物溶解（熱解)、溶質(zhì)擴(kuò)散、熱量交換、孔裂隙與固體骨架變形的復(fù)雜過(guò)程，其基本理論涉及流體力學(xué)、物理化學(xué)、傳熱傳質(zhì)學(xué)、固體力學(xué)等基礎(chǔ)學(xué)科，并與眾多工程學(xué)科相互交叉。

由于在礦體中礦物溶解與溶浸液滲透相互促進(jìn)作用，原位溶浸采礦過(guò)程實(shí)質(zhì)為一多孔介質(zhì)多場(chǎng)耦合作用過(guò)程，其理論架構(gòu)包括四個(gè)方面：

1)耦合作用本構(gòu)規(guī)律，研究某一物理場(chǎng)本構(gòu)規(guī)律和控制方程形式受其他物理場(chǎng)作用而發(fā)生變化；

2)多場(chǎng)耦合作用控制方程組或數(shù)學(xué)模型，研究某一物理場(chǎng)方程中因變量或源匯項(xiàng)受其他物理場(chǎng)作用而變化的數(shù)學(xué)描述，同時(shí)包括本構(gòu)規(guī)律影響在控制方程中的反映；

3)多孔介質(zhì)多場(chǎng)耦合作用控制方程組的求解方法及數(shù)值模擬理論與技術(shù)；

4)依據(jù)耦合理論求解方法而實(shí)施的工程方案制定及復(fù)雜規(guī)律的研究。原位溶浸采礦理論架構(gòu)中，各物理場(chǎng)影響因素及其相互耦合作用機(jī)制如圖7—10所示。

圖7 固體結(jié)構(gòu)及其力學(xué)性態(tài)影響因素

4 原位溶浸采礦本構(gòu)規(guī)律與數(shù)學(xué)模型

根據(jù)上述理論架構(gòu)，原位溶浸采礦本構(gòu)規(guī)律應(yīng)當(dāng)包括固體變形、流體運(yùn)移、傳熱傳質(zhì)，以及各物理場(chǎng)之間耦合因子變化規(guī)律。每一物理場(chǎng)的本構(gòu)規(guī)律中又存在多個(gè)影響因素，需要進(jìn)行大量的實(shí)驗(yàn)與理論研究。圖11為框圖式的本構(gòu)規(guī)律組成。在各個(gè)物理場(chǎng)本構(gòu)規(guī)律建立的基礎(chǔ)上，即可建立用于數(shù)值模擬求解的耦合數(shù)學(xué)模型。如圖12所示為原位溶浸采礦多物理場(chǎng)耦合數(shù)學(xué)模型，巖鹽水溶開(kāi)采耦合數(shù)學(xué)模型可表示為：

圖8 流體滲流及流場(chǎng)性態(tài)影響因素

圖9 溫度場(chǎng)影響因素

圖10 濃度場(chǎng)傳質(zhì)影響因素

圖11 原位溶浸采礦本構(gòu)規(guī)律組成

圖12 原位溶浸采礦多場(chǎng)耦合數(shù)學(xué)模型

原位注蒸汽開(kāi)采的熱流固耦合控制方程可表示為：

針對(duì)不同工程及開(kāi)采技術(shù)問(wèn)題，依據(jù)上述數(shù)學(xué)模型進(jìn)行精細(xì)建模，輔以一定邊界約束條件，采用一定的解耦及模擬求解方法，可獲得相應(yīng)的工程模擬結(jié)果，為相應(yīng)工程實(shí)踐提供重要指導(dǎo)。

5 結(jié)束語(yǔ)

原位溶浸開(kāi)采是地下固體礦床流體化開(kāi)采的一種方法。將溶浸、熱解流體通過(guò)鉆孔注入固體礦層，與礦層有用組分發(fā)生物理、化學(xué)、熱、力耦合作用，將礦層中有用組分轉(zhuǎn)變?yōu)榱黧w或附著成分解吸游離，通過(guò)滲流方式傳輸至鉆井而后排采至地面，具有安全、高效、環(huán)保的巨大技術(shù)優(yōu)勢(shì)，在地下礦產(chǎn)資源開(kāi)采、尤其是深部資源開(kāi)采中具有極為廣闊的發(fā)展及應(yīng)用前景。鹽類(lèi)礦床的水溶開(kāi)采，鈾礦、銅礦、金礦等的化學(xué)溶浸開(kāi)采，煤的地下原位熱解開(kāi)采，通過(guò)驅(qū)替技術(shù)而實(shí)施的煤層氣開(kāi)采，油頁(yè)巖油氣原位注熱開(kāi)采等等許多地下礦物的開(kāi)采，都屬于原位溶浸開(kāi)采之范疇。

在實(shí)驗(yàn)室內(nèi)，對(duì)鈣芒硝鹽礦的溶解滲透試驗(yàn)、對(duì)油頁(yè)巖及不同煤種的熱解試驗(yàn)結(jié)果表明：在礦體溶解或熱解過(guò)程中，礦體中固體骨架及孔隙結(jié)構(gòu)隨溶解時(shí)間或熱解溫度不斷衍化，這一孔隙結(jié)構(gòu)及滲透性的改變，影響著溶浸液的滲流及固體固件的變形，二者相互作用并影響著原位溶浸開(kāi)采的效率及效益。如何提高礦體的滲透性、增強(qiáng)溶浸液滲入與生產(chǎn)物流出是原位溶浸采礦的核心與關(guān)鍵。

根據(jù)試驗(yàn)結(jié)果，結(jié)合理論分析，初步建立了原位溶浸采礦的理論架構(gòu)、本構(gòu)規(guī)律及耦合數(shù)學(xué)模型，為不同礦體原位溶浸開(kāi)采理論建立及工程模擬，提供了重要理論指導(dǎo)。原位溶浸采礦是地下礦產(chǎn)資源開(kāi)采的重要方向，在理論架構(gòu)及數(shù)學(xué)模型建立的基礎(chǔ)上，對(duì)不同條件下的本構(gòu)規(guī)律、耦合模型求解等問(wèn)題還需要進(jìn)行大量深入細(xì)致的研究工作。與此同時(shí)，在原位溶浸開(kāi)采過(guò)程中，在一定范圍內(nèi)流體運(yùn)移及污染物控制也是未來(lái)需要研究的重要課題。

［1］吳愛(ài)祥，尹升華，王洪江，等.堆浸過(guò)程溶質(zhì)運(yùn)移機(jī)理與模型［J］.中南大學(xué)學(xué)報(bào)，2006，37（2)：385-389.

［2］李壯闊，桂斌旺，段希祥.德興銅礦堆浸廠的生產(chǎn)實(shí)踐及技術(shù)研究［J］.礦冶工程，2002，22（1)：46-48.

［3］水浩東.廢酸堆浸氧化銅鋅礦工藝［J］.有色礦冶，2003，19（1)：29-30.

［4］梁衛(wèi)國(guó).鹽類(lèi)礦床水壓致裂水溶開(kāi)采的多場(chǎng)耦合理論及應(yīng)用研究［J］.巖石力學(xué)與工程學(xué)報(bào)，2005，24（6).

［5］康志勤，呂兆興，楊棟，等.油頁(yè)巖原位注蒸汽開(kāi)發(fā)的固-流-熱-化學(xué)耦合數(shù)學(xué)模型研究［J］.西安石油大學(xué)學(xué)報(bào)（自然科學(xué)版)，2008，23（4)：30-34.

［6］于艷梅，胡耀青，梁衛(wèi)國(guó)，等.應(yīng)用CT技術(shù)研究瘦煤在不同溫度下孔隙變化特征［J］.地球物理學(xué)報(bào)，2012，55（2)：637-644.

［7］梁衛(wèi)國(guó).鹽類(lèi)礦床水壓致裂水溶開(kāi)采的多場(chǎng)耦合理論及應(yīng)用研究［D］.太原理工大學(xué)，2004.

［8］梁衛(wèi)國(guó)，趙陽(yáng)升，徐素國(guó).鹽礦群井致裂控制水溶開(kāi)采技術(shù)及應(yīng)用［J］.礦業(yè)研究與開(kāi)發(fā)，2005，25（4)：7-10.

［9］孫可明，趙陽(yáng)升，楊棟.非均質(zhì)熱彈塑性損傷模型及其在油頁(yè)巖地下開(kāi)發(fā)熱破裂分析中的應(yīng)用［J］.巖石力學(xué)與工程學(xué)報(bào)，2008，27（1)：42-52.

［10］梁衛(wèi)國(guó)，吳迪，趙陽(yáng)升.CO2驅(qū)替煤層CH4試驗(yàn)研究［J］.巖石力學(xué)與工程學(xué)報(bào)，2010，29（4)：665-673.

［11］趙陽(yáng)升，胡耀青，楊棟，等.三維應(yīng)力下吸附作用對(duì)煤巖體氣體滲流規(guī)律影響的試驗(yàn)研究［J］.巖石力學(xué)與工程學(xué)報(bào)，1999，18（6)：651-653.

［12］趙陽(yáng)升.多孔介質(zhì)多場(chǎng)耦合作用及其工程響應(yīng)［M］.北京：科學(xué)出版社，2010.