馬朗凹陷蘆草溝組泥頁巖儲層微觀孔隙特征及地質(zhì)意義

郭小波，黃志龍，柳 波，陳 旋，張文丹，陳 晨，王麗麗，孫美麗

(1.中國石油大學油氣資源與探測國家重點實驗室，北京 102249;2.東北石油大學地球科學學院，黑龍江大慶 163318;3.中國石油吐哈油田分公司，新疆哈密 839009)

隨著中國對油氣資源需求的不斷加大，蘊藏在致密儲層中的非常規(guī)油氣開始受到高度重視，繼頁巖氣之后，頁巖油、致密油、致密砂巖氣又成為油氣勘探中又一戰(zhàn)略性突破領(lǐng)域[1－3]。新疆三塘湖盆地馬朗凹陷二疊系蘆草溝組二段泥頁巖的地球化學分析顯示，其有機質(zhì)豐度高，有機質(zhì)類型以Ⅰ型和Ⅱ1型為主，成熟度(Ro/%)介于0.5%～0.9%，已經(jīng)開始大量生烴，整體上為一套優(yōu)質(zhì)烴源巖，同時其也是致密儲集層。因此，該含油泥頁巖層系可統(tǒng)稱為頁巖油或源內(nèi)致密油。在該區(qū)現(xiàn)已有5口井獲得了工業(yè)油流[4]。泥頁巖作為致密儲層，具有典型的“自生自儲”特征，其納米—微米級的儲集空間表征，是深入研究富烴機理、孔滲特征等的重要基礎(chǔ)。為此，本文以馬朗凹陷蘆草溝組泥頁巖為例，首次采用微米—納米CT、比表面及孔徑分析等數(shù)字巖心技術(shù)，對蘆草溝組泥頁巖儲層微觀地質(zhì)特征進行研究，以期為頁巖油富集機理的研究與勘探開發(fā)提供科學依據(jù)。

1 樣品與實驗

本次研究主要選取馬朗凹陷蘆草溝組二段泥頁巖樣品。脈沖滲透率測試顯示，樣品基質(zhì)滲透率極低，介于 0.158 4×10－6～0.964 0×10－6μm2，平均為 0.467 1 ×10－6μm2，屬于典型的低滲透儲層。國際純化學和應用化學聯(lián)合會(IUAPC)將微觀孔隙分為3類，即微孔(直徑＜2 nm)，介孔(直徑 2 ～50 nm)和宏孔(直徑 ＞50 nm)[5]，這3類微觀孔隙在研究區(qū)泥頁巖儲層中均有發(fā)育。泥頁巖儲層微觀孔隙定性研究技術(shù)有掃描電鏡、原子力顯微鏡、CT掃描等，前兩種技術(shù)目前應用廣泛，但存在對泥頁巖有損傷，易造成一些人為儲集空間，不能呈現(xiàn)微觀孔隙三維空間分布特征等缺陷。CT掃描技術(shù)即計算機斷層成像技術(shù)(Computed Tomography)，是近年才發(fā)展起來的數(shù)字巖心技術(shù)，是基于X射線的全方位、無損傷巖石掃描成像功能，可對同一樣品進行微米—納米CT的多尺度掃描成像，獲取不同尺度下致密儲層微觀孔隙結(jié)構(gòu)的特征切片，再通過數(shù)字合成技術(shù)得到三維立體圖像。其具有能清晰展示所測樣品的內(nèi)部結(jié)構(gòu)關(guān)系及物質(zhì)組成，揭示巖石納米級三維孔隙結(jié)構(gòu)特征等功能[6－7]，是致密儲層微觀孔隙結(jié)構(gòu)研究的理想技術(shù)。本次微觀孔隙CT掃描研究，采用臺灣生產(chǎn)的 Xradia－Micro－NanoCT－200掃描儀，在溫度差不超過2℃、無塵環(huán)境下，對樣品進行二維切片掃描，進而合成三維立體圖像和動畫。

泥頁巖微觀孔隙特征定量分析技術(shù)主要有壓汞技術(shù)、氣體吸附技術(shù)。由于測試原理的差異，氣體吸附法對微孔和介孔的分析精度較高，對宏孔的分析精度較低，而壓汞法只對宏孔的分析精度較高，二者結(jié)合，可對致密儲層微觀孔隙結(jié)構(gòu)進行較全面的定量描述[8－9]?？紤]到微孔對頁巖油的富集意義不大，本次采用的氮氣吸附孔徑分析儀的理論檢測精度為0.35 nm，分析數(shù)據(jù)主要體現(xiàn)介孔特征。氮氣吸附孔徑分析試驗采用美國生產(chǎn)的全自動比表面及孔徑分析儀(型號:ASAP2020M)，檢測精度為比表面積0.000 1 m2/g，比孔容0.000 1 mL·g－1。實驗在恒溫條件下進行，最終得到等溫吸附線和脫附線，選用多點BET模型，通過線性回歸得到樣品的比表面積、比孔容等數(shù)據(jù)[10]。

2 泥頁巖儲層微觀特征

2.1 泥頁巖微觀儲集空間類型

馬朗凹陷蘆草溝組泥頁巖中微觀儲集空間類型多樣，成因也各有差異，大多為納米—微米級的微孔隙和微裂縫[11]。蘆草溝組泥頁巖有機質(zhì)豐度高，有機質(zhì)孔隙發(fā)育，一般為親油性質(zhì)，是泥頁巖中重要的儲集空間類型[11－12]。前期研究發(fā)現(xiàn)的有機質(zhì)孔多為有機質(zhì)生烴殘留孔隙，呈凹坑狀分布(見圖1A)，本次還發(fā)現(xiàn)有機質(zhì)中存在原生的有機質(zhì)孔，是有機質(zhì)自身結(jié)構(gòu)造成的，可呈管束狀，長度可達500 μm 以上，孔隙寬度約10 μm，能譜分析元素主要為碳(見圖1B，C)。本研究采用CT多尺度掃描，從三維立體角度展現(xiàn)富有機質(zhì)泥頁巖儲層中各類孔隙的結(jié)構(gòu)與空間展布。微米CT研究發(fā)現(xiàn)，有機質(zhì)在發(fā)育大量分散、孤立微孔隙的同時，也發(fā)育呈管束狀的微米級大孔隙，且連通性好(見圖1D)，與圖1B所示原生型有機質(zhì)孔相似。馬朗凹陷蘆草溝組發(fā)育富有機質(zhì)紋層巖，主要由淺色白云質(zhì)泥巖與富有機質(zhì)暗色泥巖組成，儲層物性較好，含油級別高，是富集頁巖油的主要巖石類型之一。由實驗樣品的微米CT掃描三維圖像中可見，有機質(zhì)與無機礦物呈層狀分布，在二者的接觸面處，微孔隙密集且連片分布(見圖1E)。這些界面處連片發(fā)育的孔隙，在微觀尺度上，可為烴源巖向臨近的無機礦物儲層以非幕式排運液態(tài)烴提供通道，同時也可使層面間具備儲烴和流體滲流能力。這是紋層巖作為頁巖油有利儲層的微觀機制。這些微孔隙的大小多在納米級別，也有孔隙達到微米級別。研究認為，烴源巖中烴類分子、瀝青質(zhì)、鏈烷烴和甲烷的形成演化呈現(xiàn)一種連續(xù)譜的特征，分子直徑大小自瀝青質(zhì)的10 nm變化到甲烷的0.38 nm[6]。由此可見，泥頁巖儲層中大多數(shù)孔隙對烴類的賦存是有效的。

圖1 馬朗凹陷蘆草溝組泥頁巖中微孔、微縫特征Fig.1 Micro－pore and micro－crack in Lucaogou Formation shale in Malang Sag

納米CT顆粒分割圖像中可見微裂縫和粒間微孔隙，粒間微孔隙屬于殘留的原生孔隙(見圖1F)。泥頁巖儲層的微觀孔隙分布具有明顯的非均質(zhì)性，不同位置的微孔隙大小、形態(tài)、豐度存在差異，這種微觀差異會影響致密儲集空間的連通性和儲層的滲透性。納米CT孔隙連通性分析表明，微裂縫的連通性好，具有溝通微孔隙、改善滲透性的能力;其他呈孤立狀的納米孔隙連通性較差，但處于烴源巖內(nèi)部，具有一定的儲烴能力;集中分布的微孔隙可相互疊置，表現(xiàn)出較好的連通性(見圖1G)。在無損傷條件下，CT掃描展現(xiàn)了泥頁巖儲層內(nèi)部微孔隙的結(jié)構(gòu)與分布，可見有相當部分的微孔隙是封閉的、不連通的，這會使一般技術(shù)下測得的孔隙度值偏低。此外，蘆草溝組泥頁巖中還發(fā)育碳酸鹽巖溶蝕孔、長石溶蝕孔、白云石晶間孔、黃鐵礦晶間孔等[11，13]，連同微裂縫、有機質(zhì)孔共同組成了泥頁巖的儲集空間。

2.2 泥頁巖儲集空間的結(jié)構(gòu)特征

2.2.1 宏孔的結(jié)構(gòu)特征 壓汞數(shù)據(jù)顯示，蘆草溝組泥頁巖孔喉直徑平均值最大為911 nm，宏孔的體積百分數(shù)為68.1%(見圖2)。由此可見，蘆草溝組泥頁巖中，宏孔所占的比例很大。同時，孔隙結(jié)構(gòu)以微孔、細喉、細歪度、分選差為主要特征。

圖2 馬朗凹陷蘆草溝組泥頁巖孔隙直徑分布Fig.2 Pore diameter distribution of Lucaogou Formation shale in Malang Sag

2.2.2 介孔結(jié)構(gòu)特征 對蘆草溝組13件泥頁巖樣品進行氮氣等溫吸附法實驗，得到吸附－脫附曲線。實驗中，隨著相對壓力的增加，氮氣逐漸吸附于泥頁巖孔隙內(nèi)壁，之后壓力降低，出現(xiàn)解吸。不同微孔隙的形態(tài)和結(jié)構(gòu)不同，吸附－脫附過程不同，在吸附－脫附曲線上形成的形態(tài)也不同。因此，可以通過吸附－脫附曲線的形態(tài)對泥頁巖樣品的孔隙結(jié)構(gòu)進行研究。若吸附線沒出現(xiàn)吸附飽和狀態(tài)，表明樣品中有較大的孔隙存在;形成吸附回線，表明孔隙呈開放狀態(tài);回線越明顯，吸附線上升速率越大，表明孔隙的開放程度也就越高[10，14]。

本次分析的樣品都沒出現(xiàn)吸附飽和狀態(tài)，表明樣品中均有較大孔隙的存在。根據(jù)曲線特征，可將實驗樣品分為3類，典型曲線形態(tài)見圖3。Ⅰ類主要為蘆草溝組二段中部的樣品，吸附回線明顯，上升速率整體較高，表明孔隙開放程度較高;Ⅱ類主要為蘆草溝組二段上部的樣品，吸附回線出現(xiàn)于相對壓力(P/P0)高值段，吸附曲線上升速率較大，表明介孔中直徑較大的孔隙開放程度相對比較高;Ⅲ類主要為蘆草溝組二段下部的樣品，吸附回線不明顯，雖然開始出現(xiàn)回線的相對分壓較Ⅰ類低，但吸附曲線平緩且上升速率低，最大吸附量也較低，表明樣品中缺少大孔隙，且開放程度低。脈沖滲透率平均值顯示，Ⅰ類樣品為0.564 ×10－6μm2，Ⅱ類樣品為 0.378 ×10－6μm2，Ⅲ類樣品為0.312×10－6μm2。由此可見，隨埋深增加，孔隙的開放程度先變好，后變差;埋深較大時，較大孔隙存在的概率降低，甚至消失。同時，介孔開放程度與泥頁巖滲透性具有相同的變化規(guī)律，說明介孔的形態(tài)與分布應是影響泥頁巖儲層滲透性的重要微觀因素。所以，以Ⅰ類樣品為代表的蘆草溝組二段中部泥頁巖儲層，微孔隙發(fā)育條件最好，巖石基質(zhì)滲透性也相對最好。

3 礦物成分對儲層微觀性質(zhì)的影響

3.1 儲層物性與礦物組成

對90件蘆草溝組泥頁巖樣品進行X衍射分析，其礦物有石英、鉀長石、斜長石、方解石、白云石、黏土、黃鐵礦和方沸石。其中，碳酸鹽巖量最高，質(zhì)量百分含量分布在8.0% ～74.0%，平均為34.1%(白云石平均質(zhì)量百分含量為25.4%，方解石平均質(zhì)量百分含量為15.7%);石英質(zhì)量百分含量次之，在10.0% ～62.0%，平均為31.1%;黏土礦物質(zhì)量百分含量在4.0% ～38.0%，平均為14.6%;多數(shù)樣品含有斜長石(質(zhì)量百分含量2.0% ～29.0%)、鉀長石(質(zhì)量百分含量1.0% ～25.0%)，少量樣品檢測到黃鐵礦和方沸石。黏土礦物組分中，沒有蒙脫石，伊利石或伊/蒙混層占主導，伊利石質(zhì)量百分含量最高為96%，平均為34.5%，伊/蒙混層質(zhì)量百分含量為4% ～100%，平均為52.3%。由此可見，蘆草溝組二段泥頁巖脆性礦物含量高，塑性礦物含量低，儲層的脆性和可改造性強。

圖3 馬朗凹陷蘆草溝組泥頁巖典型樣品等溫吸附曲線Fig.3 The typical isothermal adsorption curve of Lucaogou Formation shale samples in Malang Sag

3.2 礦物成分對比表面積和比孔容的影響

根據(jù)等溫吸附－脫附曲線，通過多點BET模型計算表明，蘆草溝組泥頁巖比表面積介于0.26～2.68 m2·g－1，平均為 1.34 m2·g－1;總比孔容介于2.665 ×10－3～ 9.878 ×10－3mL·g－1，平均為6.608×10－3mL·g－1;介孔直徑介于 3.98 ～83.41 nm，平均為27.85 nm。比表面積是指單位質(zhì)量物料所具有的總面積，比表面積越大，表明物料的吸附能力越強。在泥頁巖儲層中，其表示礦物對烴類的吸附能力，比表面積越大，越不利于烴類的滲流，使頁巖油的可采性降低。比孔容是指單位質(zhì)量物料中孔隙的容積，與比表面積結(jié)合可反映微觀孔隙的大小與分布。比表面積和微觀孔隙分布受泥頁巖中礦物組成及含量的影響。

分析顯示，黏土礦物與比表面積呈正相關(guān)，與比孔容關(guān)系不明顯(見圖4A)，但整體關(guān)系表明，黏土礦物主要是有利于分散狀、小直徑孔隙的形成，從而增大比表面積，增強儲層的吸附能力。碳酸鹽巖和長石的含量均與巖樣比表面積呈負相關(guān)，而比孔容與碳酸鹽巖正相關(guān)(見圖4B，C)，與長石關(guān)系不明確(見圖4C)。這表明，碳酸鹽巖和長石的存在有利于分布集中、直徑較大的孔隙的形成(比如次生溶蝕孔隙)，這與掃描電鏡下觀察到的結(jié)果一致。而石英含量與比表面積呈正相關(guān)，與比孔容關(guān)系不確定(見圖4D)。單一礦物對比孔容的影響均不明顯，這正是泥頁巖儲集空間類型多樣、影響因素復雜的體現(xiàn)。但石英含量會增加巖石比表面積，表現(xiàn)為增強儲層對烴類的吸附能力，降低其流動性，這與傳統(tǒng)認識存在一定差異。初步分析，這可能與蘆草溝組泥頁巖中自生微晶石英的形成有關(guān)。

通過掃描電鏡可觀察到，蘆草溝組泥頁巖中存在一定數(shù)量的自生微晶石英，晶形一般較好，大小分布在5～30 μm，充填于孔隙之中或者鑲嵌于其他礦物之上，可孤立分布，也可集中分布(見圖5)。從自生石英的分布形態(tài)分析，自生微晶石英充填于儲層的大孔隙，形成微小孔隙，或自生石英以孤立狀分布，這應是造成隨石英含量增加泥頁巖比表面積增大的主要原因。

Thybegr(2010)和Metwally(2012)分別通過實驗研究表明，泥頁巖成巖過程中黏土礦物轉(zhuǎn)化，尤其是蒙脫石向伊利石轉(zhuǎn)化會釋放大量的二氧化硅，并可在泥頁巖內(nèi)部就地沉淀，形成自生微晶石英，增強泥頁巖的脆性[15－16]。蘆草溝組泥頁巖中黏土礦物含量相對較低，其形成的自生石英數(shù)量也有限。馬朗凹陷蘆草溝組作為兩套火山巖系(上覆條湖組與下伏卡拉岡組)之間的一套近海陸相湖盆沉積[17－18]，其沉積物不免會受到火山活動、地幔熱流體的影響。以此推測，大量出現(xiàn)的自生微晶石英還應有深部熱液成因或其他成因。雖然自生石英可充填于泥頁巖孔隙，阻塞喉道，但同時具有增強泥頁巖脆性、有助于形成微裂縫等作用。分析表明，石英含量與泥頁巖總孔隙度具有良好的正相關(guān)關(guān)系(見圖6)。泥頁巖中自生石英的定量分析還存在較大困難，但據(jù)此綜合分析，石英或者自生石英對泥頁巖有利儲層形成的建設(shè)性作用強于破壞性作用。

圖4 馬朗凹陷蘆草溝組泥頁巖比孔容、比表面積與礦物的關(guān)系Fig.4 Relationships of specific pore volume，specific surface and minerals of Lucaogou Formation shale in Malang Sag

4 儲層微觀特征研究的地質(zhì)意義

蘆草溝組泥頁巖油藏，不同于傳統(tǒng)的泥巖裂縫性油藏，其儲層宏觀裂縫不發(fā)育，儲集空間以微孔隙和微裂縫為主;其也不同于含氣泥頁巖儲層，頁巖氣可以擴散方式在納米級微孔隙中發(fā)生滲流[19]，由于油的高密度、低流動性等特征，對儲層微觀性質(zhì)的要求較高。泥頁巖儲層中，頁巖油可以游離態(tài)、吸附態(tài)和溶解態(tài)存在。由于賦存形式的多樣性，使得馬朗凹陷蘆草溝組泥頁巖錄井含油級別較高，油跡、油斑普遍存在，但勘探實際表明，錄井顯示好，并不都預示著是頁巖油的有利儲層發(fā)育段。小孔隙直徑、大比表面積的儲層，容易對烴類形成分子篩效應，不利于烴類的滲流。以當前的開發(fā)技術(shù)，賦存于連通性較好、較大的微孔隙或微裂縫中的游離態(tài)頁巖油才最有意義。通過以上分析可知，馬朗凹陷蘆草溝組泥頁巖儲層可歸納為以下特征:①黏土礦物含量低。作為增加泥頁巖比表面積和增強儲層塑性的主要礦物，其含量較低，使得泥頁巖吸附能力降低，儲層脆性加強，有利于游離態(tài)頁巖油的賦存。②泥頁巖儲層中孔隙類型多樣，宏孔占較大比例，介孔也可具有相對較好的連通性。同時，蘆草溝組泥頁巖主生烴段的泥頁巖具有較強的生烴能力，會使連通性較差的微孔隙也成為有效的儲油空間，在自然裂縫或人工裂縫的溝通下可具備滲流能力，成為可采資源。③碳酸鹽巖、長石等有利于直徑較大、分布集中的孔隙形成，石英(包含自生石英)對提高儲層物性、可改造性具有明顯的建設(shè)性作用。這3種礦物在蘆草溝組泥頁巖中含量高，有效地改善了泥頁巖的儲集性能?；谝陨夏囗搸r儲層的微觀地質(zhì)條件，馬朗凹陷蘆草溝組有利于游離態(tài)頁巖油的賦存與滲流，這是建成頁巖油產(chǎn)能的重要因素。因此，在頁巖油勘探開發(fā)中，需充分考慮泥頁巖儲層的微觀地質(zhì)條件，這也是影響頁巖油富集與可采能力的重要因素。

圖5 馬朗凹陷蘆草溝組泥頁巖中自生微晶石英Fig.5 Authigenic microcrystalline quartz in in Lucaogou Formation shale in Malang Sag

圖6 馬朗凹陷蘆草溝組泥頁巖中石英含量與孔隙度關(guān)系Fig.6 Relationships of porosity and quartz content of Lucaogou Formation shale

5 結(jié)論

1)采用CT多尺度掃描技術(shù)對馬朗凹陷蘆草溝組泥頁巖進行三維刻畫，表明其儲層中微孔隙分布具有較強的空間非均質(zhì)性;無機礦物和有機質(zhì)中均可發(fā)育較大孔隙;納米孔隙通過連片分布、相互疊置或被微裂縫溝通，可具有相對較好的連通性。

2)蘆草溝組泥頁巖儲層以發(fā)育宏孔為主;隨埋深增加，介孔的開放程度先變好，后變差，與巖石滲透性變化規(guī)律一致;介孔既是泥頁巖儲集空間的重要組成部分，其結(jié)構(gòu)形態(tài)也是影響儲層滲透性的重要微觀因素。

3)蘆草溝組泥頁巖中碳酸鹽巖、石英、長石含量高，碳酸鹽巖、長石有利于集中分布的大孔隙形成，石英(包含自生石英)可增強巖石脆性形成微裂縫、抵抗壓實作用等，對儲層物性具有積極作用;蘆草溝組泥頁巖儲層具有較低的比表面積;較大的微孔隙發(fā)育且連通性較好，使其有利于游離態(tài)頁巖油的富集，并建成產(chǎn)能。

[1]楊華，李士祥，劉顯陽.鄂爾多斯盆地致密油、頁巖油特征及資源潛力[J].石油學報，2013，34(1):1－11.

[2]魏虎，孫衛(wèi)，李達，等.鄂爾多斯盆地致密砂巖氣藏微觀孔隙結(jié)構(gòu)對氣井生產(chǎn)的影響[J].西北大學學報(自然科學版)，2011，41(5):869－875.

[3]鄒才能，楊智，崔景偉，等.頁巖油形成機制、地質(zhì)特征及發(fā)展對策[J].石油勘探與開發(fā)，2013，40(1):14－26.

[4]柳波，呂延防，趙榮，等.三塘湖盆地馬朗凹陷蘆草溝組泥頁巖系統(tǒng)地質(zhì)超壓與頁巖油富集機理[J].石油勘探與開發(fā)，2012，39(6):699－705.

[5]ROUQUEROL J，AVNIR D，F(xiàn)AIRBRIDGE C W，et al.Recommendations for the characterization of porous solids[J].Pure ＆ Applied Chemistry，1994，66(8):1739－1785.

[6]鄒才能，朱如凱，白斌，等.中國油氣儲層中納米孔首次發(fā)現(xiàn)及其科學價值[J].巖石學報，2011，27(6):1857－1864.

[7]白斌，朱如凱，吳松濤，等.利用多尺度 CT成像表征致密砂巖微觀孔喉結(jié)構(gòu)[J].石油勘探與開發(fā)，2013，40(3):329－333.

[8]田華，張水昌，柳少波，等.壓汞法和氣體吸附法研究富有機質(zhì)頁巖孔隙特征[J].石油學報，2012，33(3):419－427.

[9]MAYKA S，CELSO P F，JOSé A B，et al.Characterization of pore systems in seal rocks using Nitrogen Gas Adsorption combined with Mercury Injection Capillary Pressure techniques[J].Marine and Petroleum Geology，2013，39:138－149.

[10]陳尚斌，朱炎銘，王紅巖，等.川南龍馬溪組頁巖氣儲層納米孔隙結(jié)構(gòu)特征及其成藏意義[J].煤炭學報，2012，37(2):439－444.

[11]黃志龍，郭小波，柳波，等.馬朗凹陷蘆草溝組源巖油儲集空間特征及其成因[J].沉積學報，2012，30(6):1115－1122.

[12]朱日房，張林曄，李鉅源.渤海灣盆地東營凹陷泥頁巖有機儲集空間研究[J].石油實驗地質(zhì)，2012，34(4):352－356.

[13]吳林鋼，李秀生，郭小波，等.馬朗凹陷蘆草溝組頁巖油儲層成巖演化與溶蝕孔隙形成機制[J].中國石油大學學報(自然科學版)，2012，36(3):38－43.

[14]張紅日，劉常洪.吸附回線與煤的孔結(jié)構(gòu)分析[J].煤炭工程師，1993(2):23－27.

[15]THYBERG B，JAHREN J，WINJE T，et al.Quartz cementation in late Cretaceous mudstones，northern North Sea:Changes in rock properties due to dissolution of smectite and precipitation of micro－quartz crystals[J].Marine and Petroleum Geology，2010，27:1752－1764.

[16]METWALLY Y M，CHESNOKOV E M.Clay mineral transformation as a major source for antigenic quartz in thermo－mature gas shale [J].Applied Clay Science，2012，55:138－150.

[17]李紅，柳益群，梁浩，等.三塘湖盆地二疊系陸相熱水沉積方沸石巖特征及成因分析[J].沉積學報，2012，30(2):205－218.

[18]柳益群，焦鑫，李紅，等，新疆三塘湖躍進溝二疊系地幔熱液噴流型原生白云巖[J].中國科學(地球科學)，2011，41(12):1862－1871.

[19]王瑞，張寧生，劉曉娟，等.頁巖氣擴散系數(shù)和視滲透率的計算與分析[J].西北大學學報(自然科學版)，2013，43(1):75－80.