超低滲儲層孔隙演化特征及孔隙度計算
——以鄂爾多斯盆地吳倉堡油區長6段為例

張金金

(西安石油大學 地球科學與工程學院，陜西西安 710065)

超低滲儲層孔隙演化特征及孔隙度計算

——以鄂爾多斯盆地吳倉堡油區長6段為例

張金金

(西安石油大學 地球科學與工程學院，陜西西安 710065)

摘要：綜合應用巖心、電鏡、鑄體薄片、陰極發光等資料，分析鄂爾多斯盆地吳倉堡油區長6段儲層沉積作用特征和成巖作用特征。研究認為巖石成分成熟度低、不穩定組分含量高，以河口壩微相、水下分流河道微相為主的沉積砂層以及后期成巖作用強烈是研究區儲層物性差的主要因素。根據孔隙演化的特征，建立相應的孔隙度計算模型，結果表明：儲層的原始孔隙度平均值為38.6%，壓實作用、膠結作用強烈，致使孔隙損失嚴重，平均減孔率為52.1%;溶蝕作用、交代作用增加了孔隙空間，平均增孔率為8.6%;經過成巖作用過程，目前孔隙度平均值為8.2%。關鍵詞：超低滲儲層；沉積特征；成巖作用；孔隙演化；孔隙計算

儲層孔隙構成主要有原生孔隙和次生孔隙。原生孔隙主要受沉積物源、沉積環境的制約，在埋藏過程中經壓實、壓溶作用，損失嚴重；次生孔隙形成于成巖過程的溶蝕作用、交代作用，受埋藏深度、礦物成分成熟度、孔隙流體性質控制[1-4]。研究區超低滲儲層經歷漫長的沉積作用、成巖作用過程，期間伴隨著多種物理作用、化學作用。大量研究證實壓實作用、膠結作用使儲層孔隙減小，溶蝕、交代過程則會增加儲層孔隙，但大多是定性的研究。據實驗室測得鄂爾多斯盆地吳倉堡油區長6段儲層巖心平均滲透率為0.29×10-3μm2，孔隙度為8.2%，屬超低滲儲層。

研究區超低滲儲層孔隙演化過程的定量化研究理論欠缺，儲層物性受沉積特征、成巖特征的影響不明。因此，查明孔隙演化的特征，并建立孔隙度計算模型，有助于鉆前預測儲量，為油氣田勘探開發提供較準確的資料。

1孔隙演化特征分析

1.1 沉積作用特征

1.1.1 巖石學特征

礦物的成分是決定成巖過程中發生的主要反應的因素之一，不穩定組分含量高，溶蝕、交代反應較強，對次生孔隙發育的影響作用較大[5-6]。通過對區內12口井，128個薄片和巖心樣品分析認為，沉積顆粒以長石含量最高，體積分數為21%～63%，平均為47.6%，巖屑含量平均為28%，石英含量平均為24.4%。據電鏡掃描薄片、鑄體薄片及巖心分析可知，研究區礦物顆粒粒徑0.13～0.20 mm，以細粒為主；巖石分選性好，分選系數1.1～1.4，平均為1.26；顆粒之間以線接觸為主，局部凹凸接觸。

1.1.2 沉積微相特征

在研究區內長6段發育時期，湖盆進入萎縮階段，來自北東向的物源隨河流入湖后沉積在濱淺湖地帶。由于水動力減弱，懸浮顆粒下沉、堆積形成三角洲前緣亞相，是河湖共同作用地帶[7-9]。前人依據地區巖心資料和粒度分析資料，從顆粒粒度大小、分選度、層理發育、沉積韻律特征，判斷沉積水體環境，確定其沉積微相[10-12]。研究區巖心以淺灰色、灰色細砂、粉砂巖為主，粒度細，分選好。主要發育的層理有塊狀層理、波狀層理、水平層理、平行層理、透鏡狀層理，縱向上呈下粗上細的正韻律沉積，是水下分流河道微相沉積。

1.2成巖作用特征

細粒物質沉積之后，巖石即進入成巖演化階段。在巖石圍壓增大、溫度升高、孔隙流體的聯合作用下，巖石即可發生壓實、壓溶、溶蝕、膠結、交代作用[12-13]。

1.2.1 壓實作用特征

壓實作用發生于早期碎屑顆粒的埋藏，在上覆壓力不斷增大的情況下，剛性顆粒重新排列使孔隙減小，導致宏觀孔隙度減少[12-14]。壓實作用在本區的表現主要有：

(1)在高倍電鏡下觀察到的原生粒間孔隙數量少、空間小，多被巖屑、雜基充填。

(2)在掃描電鏡及鑄體薄片下觀察，顆粒接觸類型以線接觸為主，局部凹凸接觸，塑性巖屑多被擠壓變形為長條狀充填于剛性顆粒之間。

(3)部分剛性顆粒(石英、長石顆粒)發育擠壓裂紋，在裂紋發育的地方伴隨溶蝕孔的出現。

綜合鏡下薄片觀察，研究區壓實作用減孔嚴重，是孔隙演化的主要控制因素。

1.2.2 交代、溶蝕作用特征

溶蝕作用與交代作用發生于中成巖期，巖石骨架基本形成，溶蝕、交代的基本特征有：①均可以增加孔隙空間；②二者均是發生在長石、巖屑、云母表面，且都發生化學反應；③反應過程持續的時間長。綜合對比分析，巖石物性較好的層段，孔隙流體受構造運動、地溫變化、流體重力、油氣侵入的影響，進入地層，溶蝕交代較多顆粒。在地下高壓的環境下，一切反應應向巖石體積減小(即孔隙增大)的方向發展，即交代作用、溶蝕作用可以增加孔隙。在鏡下觀察，研究區內的長石主要是以鉀長石、鈣長石為主，溶蝕、交代沿解理縫、顆粒表面發生，使儲集空間增大，喉道變寬，從而大大改善了儲層物性。

1.2.3 膠結作用特征

區內膠結物總量大，體積分數為3%～27%，平均13%，對孔隙的破壞性強，是造成儲層致密的重要因素。區內膠結物主要因溶液濃度過飽和或孔隙流體酸堿度或電位平衡打破析出晶體。圖1左圖表明膠結物與石英基本無關；圖1右圖可知膠結物與長石、巖屑含量負相關，即當前長石、巖屑碎屑含量低，表示此類顆粒被溶蝕、交代而形成交代產物或膠結物，表明膠結物的主要物質來源是長石與巖屑。

碳酸鹽巖膠結物：由電鏡掃描、陰極發光資料顯示，碳酸鹽巖膠結物多充填于粒間孔隙，極大降低孔隙空間。研究區內的方解石膠結物含量為2%～18%，平均為6%。目前的研究表明鄂爾多斯盆地碳酸鹽巖主要來源有以下幾種：①早成巖階段埋藏較淺時，細菌的新陳代謝及甲烷菌產生的CO2與硫酸鹽經還原作用形成碳酸鹽膠結物，一般認為此類膠結物多被溶解；②對區內方解石膠結物的包裹體含烴鹽水測溫，以及δ13C值含量的分析認為，方解石膠結物主要生成于中成巖A期，有機質發生大量脫羥基作用[15-16]，孔隙流體呈中性至弱堿性，提供充足的碳氧化合物與孔隙流體中的鈣離子結晶析出方解石[17]；③長石云母類礦物發生的交代反應生成。

圖1　膠結物含量與長石、巖屑及石英含量對比圖

黏土膠結物：本區的黏土膠結物主要是以綠泥石、伊利石為主，分布廣，含量豐富，體積分數最高為12%，平均為3.7%。

綠泥石膠結物在研究區的主要產狀有包膜、襯里、絨球狀、蜂窩狀綠泥石，呈葉片狀綠泥石包繞碎屑顆粒，降低了孔隙空間，縮小喉道半徑，甚至堵塞喉道，使儲層物性變差。當鐵鎂礦物充足、地溫較高時利于長石的溶蝕、交代生成綠泥石。另外也有研究表明綠泥石存在的地方，較少有石英次生加大，因而可改善儲層物性[18]。伊利石的膠結物主要生成于成巖階段中-晚期。早前研究表明，早成巖時期生成的伊利石包膜對巖石骨架有支撐作用，減少壓實作用對孔隙的破壞作用，對儲層孔隙有建設性作用。黏土膠結物增加喉道彎曲度，也可堵塞喉道，導致滲透性變差[19-20]。

硅質膠結物：硅質膠結物主要以自生石英為主，本區分布廣泛，含量最高達到10%，平均含量3%，多以微晶石英集合體產出，充填于粒間孔隙當中，對孔隙的破壞作用較大；石英的次生加大在粒度較粗、黑云母含量較低、埋藏較深的砂巖中發育，主要由鉀長石、鈣長石、云母交代反應的生成物。

超低滲儲層由沉積物顆粒的堆積，形成原始孔隙空間，隨后因埋深等上覆壓力增大，機械壓實作用導致孔隙空間減少，在地下溫度升高以后，與孔隙流體作用，一些不穩定的碎屑組分如長石、云母等顆粒表面就被溶蝕或者交代了，溶蝕掉的離子與孔隙流體中原有離子重組生成膠結物，這些膠結物與交代產物占據了溶蝕孔以及原生粒間孔，致使儲層致密。

2孔隙定量計算

2.1原始孔隙

原始孔隙計算依據Beard[21]等人給出的公式：

φ1=20.91+22.90/Sd

(1)

(2)

式中：φ1為原始孔隙度；Sd為 Trask 分選系數；φ75為粒度概率累積75分位的截距。

φ2=φ等大球體排列模型/Sd

(3)

國內學者張創經分析論證認為，當分選系數小于1.5時的未固結砂巖原始孔隙度φ1與砂巖的分選系數Sd應用的效果較好。通過對研究區128個巖石粒度分析結果計算出Sd值，代入式(1)得出原始孔隙度為38.5%～41.2%，平均為39.6%。數據分析結果表明：假設顆粒是等大球體模型，此時顆粒未壓實，理論孔隙度是47.64%，將其與Sd值之比得到的值作為原始孔隙計算模型(公式3)，這個值與式(1)的值對比分析，得到圖2的校正圖版，得到研究區原始孔隙校正公式。

圖2　φ1與φ2校正圖版

2.2壓實后孔隙

根據公式(3)建立的模型，假設等大球體顆粒被完全壓實，及顆粒菱形排列，其孔隙度為25.95%，將此值代入式(3)模擬出壓實后孔隙度，將計算值代入圖2中的校正方程，即可得到壓實后孔隙度值。計算得到研究區壓實之后的孔隙度為18%～24%，平均孔隙度值21%，平均減孔率為52.1%(注：文中采用的平均增孔率或減孔率值為增加或損失孔隙平均值與原始孔隙平均值之比)。

2.3交代、溶蝕增孔計算

交代、溶蝕作用的母質主要是長石、云母，據前人的研究表明[22-24]，主要發生如下反應：

3KAlSi3O8(鉀長石)+2H+→KAl3Si3O10(OH)2(伊利石)+6SiO2(石英)+2K+

(4)

CaAl2Si2O8(鈣長石)+SiO2+5Ca2++Mg2++4H2O →MgAl2Si3O10(OH)8(綠泥石)+6CaCO3

(5)

2KAlSi3O8(鉀長石) +5Ca2++Mg2++4H2O→MgAl2Si3O10(OH)8(綠泥石)+5CaCO3

(6)

KAlSi3O8(鉀長石)+4H+→3SiO2+K++Al3+

(7)

KAl2(AlSi3O10)(OH)2(云母)→KAl3Si3O10(OH)2(伊利石)

(8)

前人研究表明：鈣長石的摩爾體積為93.48 cm3/mol，鉀長石的摩爾體積為109.1 cm3/mol，石英的摩爾體積為22.64 cm3/mol，伊利石的摩爾體積為140.7 cm3/mol，綠泥石的摩爾體積為205.85 cm3/mol，方解石的摩爾體積為36.83 cm3/mol。將數據代入上述方程得：反應式(4)增加孔隙15.51%，反應式(5)增加孔隙2.25%，反應式(6)增加孔隙27.5%，反應式(7)增加孔隙37.75%，反應式(8)對孔隙沒有影響。將此系數與相應膠結物含量相乘，即可計算出溶蝕交代作用的增孔率，計算結果為0.25%～9.5%，平均增孔率為8.6%。

2.4膠結物減孔計算

從掃描電鏡、鑄體薄片、陰極發光照片觀察，膠結、交代產物充填于粒間孔隙，占據孔隙體積較大，最小值為6%，最大值為27%，平均減孔率為36.7%。

3計算值與實測值驗證

目前孔隙是由殘余粒間孔、溶蝕孔構成，原始孔隙空間先經壓實作用，孔隙減少，后經溶蝕作用增加了孔隙空間，同時由于生成的膠結物占據了孔隙空間，剩余的孔隙體積比就是目前的孔隙度。即有：

孔隙度=原始孔隙度-壓實后孔隙度+溶蝕交代作用增孔率-膠結物減孔率

將數據代入計算得到的孔隙度與實驗巖心測出的孔隙度絕對誤差±1.1，相對誤差小于18.5%，取得較好的應用效果，見圖3。

圖3　計算孔隙度與實測孔隙度關系圖

4結論

(1)鄂爾多斯盆地吳倉堡油區長6儲層物性差，砂巖類型主要以長石砂巖為主，巖屑含量較高，成分成熟度低；沉積顆粒以細砂為主，顆粒分選好，原始孔隙與分選相關，分選越好，原始孔隙越大。

(2)研究區內壓實作用強，使孔隙損失嚴重，平均減孔率為52.1%；交代溶蝕作用較強，平均增孔率為8.6%；膠結物含量高，導致孔隙大量損失，平均減孔率36.7%。

(3)根據儲層成巖演化過程得出了超低滲儲層孔隙度的計算方法，理論值與實測值絕對誤差小于±1.1，適合研究區孔隙度計算。若將該模型應用于其它地區，則需要結合當地實際沉積特征與成巖特征，綜合巖石礦物組合、成巖埋藏史及其他條件，對該方法校正，以達到理想的應用效果。

參考文獻

[1]吳小斌,侯加根,孫衛.特低滲砂巖儲層微觀結構及孔隙演化定量分析[J].中南大學學報,2011,42(11):3439-3445.

[2]劉國勇,金之鈞,張劉平.碎屑巖成巖壓實作用模擬實驗研究[J].沉積學報,2006,24(3):408-412.

[3]于紅巖,李洪奇,黨娜,等.鄂爾多斯盆地超低滲透儲層形成的宏觀地質因素分析[J].科技導報,2012,30(9):44-47.

[4]王嵐,鄒才能,林潼,等.鄂爾多斯盆地白豹-華池地區延長組長6油層組沉積環境及成巖相分析[J].天然氣地球科學,2011, 22(5):797-805.

[5]朱筱敏.沉積巖石學[M].北京,石油工業出版社,2008:241-248.

[6]姜在興.沉積學[M].北京,石油工業出版社,2010:180-181.

[7]李書恒,方國慶,楊金龍,等.鄂爾多斯盆地超低滲透儲層成因研究及意義[J].巖性油氣藏,2012,24(6):32-36.

[8]王 嵐,鄒才能,林 潼,等.鄂爾多斯盆地白豹-華池地區延長組長6油層組沉積環境及成巖相分析[J].天然氣地球科學,2011,22(5):797-805.

[9]鄒才能,趙政璋,楊華,等.陸相湖盆深水砂質碎屑流成因機制與分布特征-以鄂爾多斯盆地為例[J].沉積學報,2009,27(6):1067-1074.

[10]蘭葉芳,鄧秀芹,程黨性,等.鄂爾多斯盆地三疊系延長組次生孔隙形成機制[J].地質科技情報,2014,33(6):128-135.

[11]徐蕾,師永民,徐常勝,等.長石族礦物對致密油儲滲條件的影響-以鄂爾多斯盆地長6油層組為例[J].石油勘探與開發,2013,40(4):448-453.

[12]黃可可,黃思靜,佟宏鵬,等.長石溶解過程的熱力學計算及其在碎屑巖儲層研究中的意義[J].地質通報,2009,28(4): 474-482.

[13]張瑞,王琪,姚涇利,等.鄂爾多斯盆地延長世湖盆中部長6段儲層成巖特征[J].天然氣地球科學,2010,21(6):890-896.

[14]曹斌風,孫衛.吳旗地區薛岔區塊延長組長6砂巖儲層成巖作用研究[J].天然氣地球科學,2011,22(6):952-959.

[15]覃建雄.流體包裹體在沉積盆地分析中的應用[J].地質科技情報,1994,13( 2):39-44．

[16]羅靜蘭,劉小洪,林潼,等.成巖作用與油氣侵位對鄂爾多斯盆地延長組砂巖儲層物性的影響[J].地質學報,2006,80(5):664-674.

[17]孟萬斌,呂正祥,馮明石,等.致密砂巖自生伊利石的成因及其對相對優質儲層發育的影響[J].石油學報,2011,32(5):783-790.

[18]黃思靜,謝連文,張萌,等.中國三疊系陸相砂巖中自生綠泥石的形成機制及其與儲層孔隙保存的關系[J].成都理工大學學報,2004,31(3):274-280.

[19]田建鋒,高永利,張蓬勃,等.鄂爾多斯盆地合水地區長7致密油儲層伊利石成因[J].石油天然氣地質,2013,34(5):701-706.

[20]王 芳,馮勝斌,何 濤,等.鄂爾多斯盆地西南部延長組長7致密砂巖伊利石成因初探[J].西安石油大學學報,2012,27(4):20-22.

[21]Beard D C, Weyl P K.Influence of texture on porosity and permeability of unconsolidated Sand[J].AAPG Bulletin, 1973,57(2): 349-369

[22]潘兆櫓.結晶學及礦物學[M].北京.地質出版社,1994：35-46.

[23]陳世悅.礦物巖石學[M].山東東營,中國石油大學出版社,2002:80-120.

[24]羅明高.定量儲層地質學[M].北京,地質出版社,1998:19-34.

編輯：王金旗

文章編號：1673-8217(2016)02-0056-04

收稿日期：2015-10-26

作者簡介：張金金,1989年生，西安石油大學地質工程專業在讀碩士研究生，從事地球物理測井技術方向的研究。

中圖分類號：TE348

文獻標識碼：A