渝東北地區(qū)WX2井頁巖氣賦存特征及其勘探指示意義

付常青，朱炎銘，陳尚斌，梁 峰

1.中國礦業(yè)大學煤層氣資源與成藏過程教育部重點實驗室，徐州221116；2.中國礦業(yè)大學資源與地球科學學院，徐州221008；3.中國石油勘探開發(fā)研究院廊坊分院，廊坊065007

渝東北地區(qū)WX2井頁巖氣賦存特征及其勘探指示意義

付常青1,2，朱炎銘1,2，陳尚斌1,2，梁 峰2,3

1.中國礦業(yè)大學煤層氣資源與成藏過程教育部重點實驗室，徐州221116；2.中國礦業(yè)大學資源與地球科學學院，徐州221008；3.中國石油勘探開發(fā)研究院廊坊分院，廊坊065007

為了深入研究渝東北地區(qū)龍馬溪組頁巖氣賦存特征，該文以WX2井頁巖高溫高壓等溫吸附及覆壓孔隙度實驗數據為基礎，通過誤差最小原則挑選了適合研究區(qū)的吸附模型，并基于孔隙度隨有效應力變化關系建立游離氣模型，綜合分析了吸附氣、游離氣及總含氣隨埋藏深度的變化特征。研究結果表明：WX2井頁巖不同溫度下過剩吸附量隨著壓力增大，均呈現先增大后減小的趨勢，隨著溫度的升高，最大吸附量逐漸減小，而校正后的絕對吸附量隨壓力增加，先迅速增大后增速放緩，且用D-A模型擬合絕對吸附量數據平均誤差最小，基本可以反應研究區(qū)頁巖真實吸附過程。頁巖樣品在加壓過程中孔隙及微裂隙會逐漸閉合，卸壓時絕大部分會重新打開，存在部分塑性變形造成的不可逆損傷，但不可逆損傷所占比重較輕。不同方向樣品孔隙度與有效應力之間具有負指數關系，富含層理頁巖平行樣品較垂直樣品具有更大的初始孔隙度以及更強的孔隙應力敏感性。頁巖氣賦存特征綜合受控于儲集層特征、吸附能力、溫度及壓力等因素，其中溫度對吸附氣和游離氣含量為負效應，儲層壓力為正效應；吸附氣、游離氣及總含氣量均遵循先增大后減小的總體趨勢，其中吸附氣及游離氣含量分別主要受控于溫度及儲層壓力。此外，臨界深度上下，頁巖吸附態(tài)與游離態(tài)相對含量發(fā)生變化，其對頁巖氣富集評價具有重要意義。

WX2井；賦存特征；含氣量；臨界深度；渝東北

中國南方地區(qū)下古生界五峰組—龍馬溪組海相頁巖因具有厚度大、分布穩(wěn)定、有機碳含量高、熱演化程度高以及含氣性良好等特點，成為頁巖氣勘探研究的重點目標層位（董大忠等，2012；郭旭升，2014；鄒才能等，2015）。長期以來，頁巖氣研究主要集中在以有利區(qū)優(yōu)選為目標的儲層特征及成藏條件方面（朱炎銘等，2010；聶海寬等，2012；紀文明等，2016），且近年來，長寧、昭通、威遠及焦石壩區(qū)塊頁巖氣成功實現商業(yè)開發(fā)，顯示我國頁巖氣勘探開發(fā)取得突破性進展（郭彤樓和張漢榮，2014；梁峰等，2016）。但除上述4個區(qū)塊外，其他區(qū)域目前尚未獲商業(yè)氣流也暴露出我國在正確認識頁巖氣賦存特征及揭示頁巖氣富集規(guī)律方面存在瓶頸。因此，下一階段研究重心將轉向以尋找頁巖氣富集區(qū)帶為目標的頁巖氣富集規(guī)律方面（金之鈞等，2016；何治亮等，2016）。沉積環(huán)境、儲層物性、頁巖氣賦存特征、氣體運移特征及頁巖氣保存條件等多種因素均會影響頁巖氣的富集（夏威等，2015；李艷芳等，2015；夏陽等，2015；劉樹根等，2016；劉軍等，2015）。因頁巖氣含量是頁巖氣富集程度最直觀的參數，并對頁巖氣地質儲量、可采資源量預測具有重要意義（郭旭升等，2016）；同時，四川盆地及其周緣五峰組—龍馬溪組頁巖氣勘探實踐揭示，厚度、有機碳含量、熱演化程度以及頂底板條件相似的頁巖，其含氣性存在很大差異（魏志紅，2015）；其次，我國南方古生界富有機質頁巖層一般均經歷了早期深埋及晚期抬升兩個階段，受后期構造作用的影響，使得現今頁巖氣層處于不同埋藏深度，但對不同埋藏深度下，頁巖氣各賦存相態(tài)含氣量變化及其對含量研究成果則極其匱乏，因此，研究頁巖氣賦存特征就顯得尤其重要。中國石油勘探開發(fā)研究院廊坊分院在四川盆地邊緣渝東北巫溪部署的WX2井，五峰組—龍馬溪組富有機質頁巖在有機質豐度、脆性礦物含量、富有機質頁巖連續(xù)厚度等多項靜態(tài)參數均優(yōu)于目前商業(yè)開發(fā)區(qū)塊的評價井，且頁巖最高實測含氣量超8m3/t（梁峰等，2016；鄒才能等，2016），因而正確認識WX2井頁巖氣賦存特征對揭示渝東北地區(qū)頁巖氣富集規(guī)律及指導頁巖氣勘探方面具有重要價值。

南方下古生界五峰組—龍馬溪組頁巖處于高過成熟度階段（鄒才能等，2010），頁巖氣除少部分呈溶解狀態(tài)賦存于有機質和結構水以外，絕大部分以游離狀態(tài)賦存于孔隙和裂隙中，或以吸附狀態(tài)賦存于礦物顆粒和有機質的表面（張金川等，2004；Curtis et al.，2002），頁巖氣總含量基本是游離態(tài)頁巖氣和吸附態(tài)頁巖氣之和。對頁巖氣的賦存特征研究對容積法評價頁巖氣地質儲量有重要意義（張雪芬等，2010）。等溫吸附是研究頁巖吸附氣含量的重要途徑，而吸附模型是對等溫吸附曲線進行非線性擬合得到具有物理意義的參數，從而對吸附過程進行描述并對吸附氣含量進行預測。頁巖吸附模型主要有Lanmuir模型（L模型）、Freundlich模型（F模型）、Langmuir-Freundlich模型（L-F模型）、Dubibin-Astakhov模型（D-A模型）及Dubibin-Radushkevich模型（D-R模型）（Harpalani et al.,2006;Ji et al.,2015;Singh and Javadpour，2016;Wang etal.,2016），但各模型在研究區(qū)的適用性未有研究成果提及。游離氣含量多采用PVT方程進行計算，其含量受外界環(huán)境及儲集層孔隙特征的控制，而現階段針對儲集層特征的研究主要集中在孔隙表征及分類，缺少應力條件下頁巖孔隙變化研究，造成對不同溫壓條件下游離氣含量認識不清。鑒于此，本文通過WX2井龍馬溪組鉆井巖心實驗分析，旨在揭示以埋藏深度為基礎的儲層溫度及壓力控制下的頁巖氣賦存特征，為渝東北地區(qū)海相頁巖氣富集規(guī)律研究及勘探開發(fā)提供依據。

圖1 WX2井位置及采樣點信息Fig.1 The location ofWX2 welland sample information

1 地質概況

WX2井地處重慶東北部巫溪縣東（圖1），構造上位于南大巴山弧形褶皺帶南部，鐵溪—巫溪隱伏斷裂帶以北，構造位置位于田壩背斜北翼，背斜地表構造形態(tài)呈現箱狀背斜或斷背斜形態(tài)，背斜兩翼為斷層，地表出露地層為上二疊統，背斜北翼地層傾角10°～50°，深部地層平緩（梁峰等，2016；鄒才能等，2016）。該井從頂至底鉆遇地層分別為上二疊統長興組、吳家坪組，下二疊統茅口組、棲霞組、銅礦溪組，中志留統徐家壩組、下志留統龍馬溪組，上奧陶統五峰組、臨湘組等，其中缺失上志留統、泥盆系、石炭系。WX2井所在地區(qū)地溫梯度約19.5℃/km（鄧賓等，2013；湯慶艷等，2013；曹環(huán)宇等，2016），龍馬溪組底部頁巖儲層壓力梯度約1.4MPa/hm（劉若冰，2015；郭彤樓，2016；劉洪林等，2016），為超壓儲層。目標層為五峰組—龍馬溪組，其中TOC大于2.0%層段位于1543～1636m井段（梁峰等，2016）。五峰組—龍馬溪組富有機質頁巖段巖性以硅質頁巖、炭質泥頁巖、黑色泥頁巖及黑色頁巖為主，其中黑色頁巖層理及黃鐵礦晶體十分發(fā)育，筆石化石豐富，缺乏底棲生物，但含有低等菌藻類及硅質海綿骨針和放射蟲，指示深水、低能、缺氧的深水陸棚相沉積環(huán)境（熊小輝等，2015）。

2 樣品與實驗

2.1 樣品

本次取樣點位于WX2井實測高含氣層段，取樣深度1623.20m，實測地溫46.65℃，儲層壓力22.72MPa，現場解吸含氣量約4.90mL/g（圖1）。樣品中含有大量筆石化石，頁巖層理極為發(fā)育，有機質類型為Ⅰ-Ⅱ1型，Ro值為2.76%，處于過成熟演化階段（梁峰等，2016）。

2.2 實驗方法（1）高壓甲烷等溫吸附

將巖芯樣品破碎至60～80目（0.18～0.25mm），取100～150 g空氣干燥基樣品，在貴州省煤田地質局采用GAI-100型高壓氣體等溫吸附儀進行實驗。儀器最高壓力69 MPa，壓力傳感器精度達到0.05%，實驗溫度最高177℃，采用先進的數字式油浴，控溫精度為0.1℃。為研究不同埋藏深度對應儲層溫度下頁巖吸附能力，實驗中最高吸附壓力為29MPa，吸附溫度為30℃、40℃、50℃、60℃、70℃及80℃。每個壓力點吸附平衡時間一般大于12 h，甲烷吸附介質純度為99.99%。因現階段頁巖氣吸附實驗尚未發(fā)布國家標準，參照《煤的高壓等溫吸附試驗方法》進行實驗。

（2）覆壓孔隙度

將實驗樣品按照《工程巖體試驗方法標準》要求，分別沿層理方向和垂直層理方向加工成直徑25mm高度50mm的標準圓柱試樣（后文中分別稱為：平行樣及垂直樣），加工后的試樣兩端面平整度偏差在0.05mm內，試樣側面光滑筆直，軸向垂直度不超過0.001 rad。實驗在貴州省煤田地質局采用運用PDP-200覆壓孔滲測試儀進行。每個壓力點加壓時間大于2 h，實驗氣體為氦氣，介質純度為99.99%，實驗中保持孔隙壓力不變，通過改變圍壓來改變有效應力從而完成受力過程中孔隙度的測試。

表1 龍馬溪組頁巖等溫吸附實驗測試數據Table1 Isotherm experimentaldataof the LongmaxiFormation

3 結果與討論

3.1 結果

3.1.1 等溫吸附實驗結果

WX2井頁巖等溫吸附測試結果見表1。不同溫度下，吸附量隨壓力的增大均呈現先增大后減小的趨勢，其中30℃時吸附量介于0～2.59 m L/g，40℃時吸附量介于0～2.14mL/g，50℃時吸附量介于0～1.85mL/g，60℃時吸附量介于0～1.68mL/g， 70℃時吸附量介于0～1.54mL/g，80℃時吸附量介于0～1.47 mL/g。各溫度條件下，最大吸附量對應壓力均在6MPa附近，且伴隨著溫度的升高，最大吸附量逐漸減小。

3.1.2 覆壓孔隙度實驗結果

受控于實驗儀器限制，實驗過程中，設置初始圍壓2MPa，步長4MPa，最大圍壓22MPa，實驗過程中保持孔隙壓力0.5MPa不變，WX2井頁巖覆壓孔隙度測試結果見表2。圍壓從2～22MPa，平行樣品孔隙度從4.69%降至2.63%，垂直樣品孔隙度從4.10%降至2.68%。兩樣品隨著圍壓的增大孔隙度均逐漸減小，并表現出初期迅速衰減，后期衰減幅度放緩的整體趨勢。

表2 樣品覆壓孔隙度測試數據Table2 Sample porosity data in confining stress

3.2 討論

3.2.1 吸附模型適用性

現有研究表明，體積法所測試的吸附數據是過剩吸附量，也稱作吉布斯表面過剩吸附量（Harpalanietal.，2006；Wang etal.，2016），在高壓吸附實驗中，過剩吸附量明顯小于相應的絕對吸附量，絕對吸附量的計算公式（郇璇等，2015）：

式中：nabs（mol/g）為絕對吸附摩爾量，nGibbs（mol/g）為壓力p（MPa）時過剩吸附摩爾量，ρsorbed（mol/mL）為吸附相的平均密度摩爾密度，ρgas（mol/m L）為壓力與溫度的函數。

將吸附摩爾量轉換成吸附體積，則絕對吸附量計算公式為：

式中：Vabs（mL/g）為絕對吸附量；VGibbs（mL/g）為過剩吸附量；ρsorbed多取定值0.421g/mL（Wang etal.，2016），ρgas按照標準氣態(tài)方程轉換后的計算公式進行計算：

式中：M（g/mol）為摩爾質量，甲烷為16，R（cm·MPa/mol·K）為理想氣體常數，取值8.314，T（K）為熱力學溫度。

將表1實測值帶入公式（2），通過計算，龍馬溪組頁巖過剩及絕對等溫吸附曲線如圖2所示。校正后的絕對吸附量曲線均呈現隨測試壓力的增加，先迅速增大，后增大放緩的變化趨勢；此外，隨著測試溫度的升高，絕對吸附量呈逐漸減小的整體趨勢，與過剩等溫吸附曲線保持一致。

利用現階段常用吸附理論模型L模型、F模型、L-F模型、D-A模型及D-R模型（Harpalaniet al.，2006；Ji et al.，2015；Singh and Javadpour，2016；Wangetal.，2016）對頁巖絕對等溫吸附數據進行擬合研究。

（1）L模型

Langmuir從動力學觀點出發(fā)，提出了單分子層吸附的狀態(tài)方程。L模型的表達式為：

式中：V（mL/g）為吸附平衡時的吸附量；VL（m L/g）為Langmuir體積，即飽和吸附量；p（MPa）為平衡壓力；pL（MPa）為Langmuir壓力。

（2）F模型

Freundlich通過熱力學方法推導出來的等溫吸附方程可以描述多層吸附，F模型的表達式為：

式中：a（m3/(t·MPan)）為Freundlich系數；n為Freundlich指數。

（3）L-F模型

L-F模型方程同時考慮了吸附劑表面的非均勻性以及被吸附分子之間的作用等復雜因素，模型表達式為：

式中：n為與溫度及頁巖孔隙分布有關的模型參數，用來校正吸附位與吸附分子，當n=1時，即為L模型。

（4）D-A模型

基于微孔充填及吸附勢理論，Dubinin和Astakhov提出的D-A吸附模型為：

式中：V0（mL/g）為吸附劑的微孔體積；D在特定吸附質和吸附劑為常數，其表達式為(RT/E)n，其中E（L·MPa/mol）為吸附特征能，n為與吸附劑非均勻性相關的參數。

圖2 龍馬溪組頁巖過剩等溫吸附曲線（左）及絕對等溫吸附曲線（右）Fig.2 Excessadsorption isotherm(left)and absoluteadsorption isotherms(right)of the LongmaxiFormation

（5）D-R模型

Dubinin和Radushkevitvh研究表明，D-A模型中，當n=2的值可能適合某些情況，由此建立的D-R模型為：

本次利用誤差最小原則挑選適合研究區(qū)的等溫吸附模型。基于計算結果，對各模型的擬合結果與實測結果的誤差進行計算，方程為（Harpalanietal.，2006；Wangetal.，2016）：

式中：Vexp（mL/g）為實驗吸附量，Veq（mL/g）為擬合吸附量，r為誤差。

將每個溫度點下的吸附點誤差取平均值，作為此溫度下的平均誤差：

式中：ARE為每個溫度點下的平均誤差值，j為每個溫度點下的測量數，abs(ri)為每個測量點誤差的絕對值。

從各模型擬合結果的誤差來看（表3），F模型平均誤差最大，D-R模型次之，L-F模型及D-A模型擬合效果較好，誤差值均小于1%，且以D-A模型誤差最小。此外，采樣測試點儲層條件下（地溫46.65℃，儲層壓力22.72MPa），D-A模型誤差僅0.191%（測試溫度50℃），并基于研究區(qū)地溫梯度，五峰組—龍馬溪組目標層溫度在50～80℃之間，而此溫度范圍內D-A模型誤差值同樣最小。鑒于此，D-A模型最適合用來擬合研究區(qū)等溫吸附實驗數據，其基本可以反應頁巖真實吸附過程。

圖3 不同埋深條件下頁巖的吸附特征Fig.3 Adsorption characteristicsof shalesunder different burialdepths

表3 各模型擬合參數及誤差Table3 Fitting resultsand fitting errorsof differentmodels

3.2.2 吸附氣含量計算

吸附氣含量嚴格受控于頁巖的吸附能力。根據前人研究成果，研究區(qū)地溫場介于19～25℃/km，儲層壓力系數介于0.80～1.50（換算為儲層壓力梯度介于0.8～1.5MPa/hm）（鄧賓等，2013；湯慶艷等，2013；劉若冰，2015；曹環(huán)宇等，2016；郭彤樓，2016；劉洪林等，2016）。為全面考慮溫壓控制下頁巖氣含量，研究中分別選取不同地溫梯度（20℃/km和25℃/km）及儲層壓力梯度（0.8MPa/hm和1.5MPa/hm）進行對比研究。假設地表溫度為恒溫帶溫度15℃，儲層壓力為0MPa，利用地溫梯度和壓力梯度，可計算不同埋深的溫度及儲層壓力，再根據優(yōu)選的D-A等溫吸附模型，即可求得各溫壓條件下頁巖吸附氣含量，以此反演不同深度（溫壓）環(huán)境中頁巖吸附氣量特征（圖3）。

預測0到4000m埋深條件下，頁巖吸附氣量呈現先增大后減小的整體趨勢，最大吸附氣含量均出現在埋深750m附近。相同埋深條件下，儲層壓力梯度越高，地溫梯度越小，則吸附氣含量越大。反之，儲層壓力梯度越低，地溫梯度越大，則吸附氣含量越小。由此可見，溫度對頁巖吸附的影響是消極的，在其他條件不變的情況下，吸附能力隨溫度增大而降低，而儲層壓力對頁巖吸附的影響卻是積極的，在其他條件不變的情況下，吸附能力隨儲層壓力的增大而升高。隨著埋藏深度的增大，相同地溫條件下吸附氣含量趨于相同，且現今全球獲工業(yè)氣流頁巖氣開采深度基本大于1000m（鄒才能等，2015），對比可知，在此埋深條件下，吸附氣含量受儲層壓力影響微弱，而主要受控于儲層溫度。

3.2.3 覆壓孔隙模型求取

眾多學者研究認為孔隙度隨有效應力增加而減小，且總結了孔隙度的變化關系式（Athy,1930; Gilesetal.,1998），較具有代表性的是：

式中：?t（%）為t時刻的孔隙度，?0（%）為初始狀態(tài)的孔隙度，Cp（MPa-1）為孔隙壓縮系數，pt（MPa）為t時刻的有效應力。

本次研究中，有效應力根據Terzaghi的定義（耶格等，1981），計算公式為：

式中：σ（MPa）為外部應力，即為本實驗中施加的圍壓；pp（MPa）為孔隙壓力。

對WX2垂直樣品先加壓而后卸壓，當最大有效應力為21.5 MPa，回復后孔隙度損失量為10.976%（圖4），表明頁巖樣品在加壓過程中孔隙及微裂隙會逐漸關閉，卸壓時絕大部分已經關閉的孔隙和微裂隙會重新打開，仍存在部分塑性變形造成的不可逆損傷。但由于不可逆損傷所占比重較輕，因此，利用頁巖覆壓孔隙度實驗，可以表征頁巖在經歷構造抬升，有效應力變化時孔隙的特征。

圖4 垂直樣品孔隙度隨有效應力相應關系Fig.4 The relationship between verticalsample porosity and effective stress

圖5 有效應力與孔隙度關系Fig.5 The relationship betw een porosity and effective stress

利用表2實驗結果，進行有效應力與孔隙度的擬合（圖5）。

不同方向樣品孔隙度隨應力擬合曲線相關系數分別為0.9803與0.9473，表明孔隙度與有效應力之間符合負指數關系。相比之下，平行樣品較垂直樣品具有更大的初始孔隙度，但隨著有效應力的增加，平行樣品孔隙度衰減幅度較大，孔隙度迅速減小并小于垂直樣品孔隙度。究其原因為：WX2井龍馬溪組富有機質頁巖頁理發(fā)育，而頁理作為薄弱面具有更多的孔隙空間（李傳亮等, 2003），使其得平行樣品較垂直樣品具有更大的初始孔隙度，但同時，因樣品具有較大的孔隙度，使具有更強的孔隙應力敏感性（Milliken et al., 2013），最終表現為隨應力增加，孔隙度降低趨勢更明顯。

對于富含層理頁巖而言，要想較為準確獲取樣品孔隙度，避免頁巖結構各向異性的影響，需綜合平行于層面與垂直于層面柱狀樣品孔隙度，取其平均值。因此，利用表2中平均孔隙度數據，建立WX2頁巖覆壓孔隙度綜合模型，見式13。

式中：?t（%）為平均孔隙度。

3.2.4 游離氣含量計算

頁巖氣中的游離氣主要儲集于頁巖孔隙中（游聲剛等，2015），因此準確獲得不同溫壓條件下頁巖孔隙度的變化，即可對頁巖中游離氣含量進行研究。由于實驗未獲得溫度對孔隙度的影響數據，本次僅研究不同儲層壓力梯度對不同埋藏深度頁巖游離氣的影響，但包含溫度對氣體體積膨脹產生的影響。

對游離氣含量的計算采用理想氣態(tài)方程：

式中：p（MPa）為某深度時儲層壓力，V（m L/g）為某深度時對應孔隙體積；T（K）為某深度時儲層開氏溫度，式中p0、V0、T0為折算到0℃、101.325 kPa時所對應的壓力、游離氣含量及開氏溫度。

由溫度與深度的關系可知，溫度與埋深成正比，建立溫度與深度的關系：

式中：h（m）為頁巖埋藏深度，dT（℃/km）為溫度梯度。

儲層壓力的計算公式為：

式中：g（N/kg）為重力加速度，值取10.0，η為儲層系數；ρW（kg/m3）為水體的密度，值取1.0× 103。

儲層有效應力計算公式為：

式中：pt（MPa）為儲層有效應力，ρR（kg/m3）為巖石平均的密度，取平均值2.6×103。

利用上述公式，分別計算不同埋深的儲層孔隙體積、溫度和壓力（與吸附氣含量計算中地溫梯度及儲層壓力梯度取值相同），再根據理想氣態(tài)方程，即可反演不同深度（溫壓）環(huán)境中頁巖游離氣量特征（圖6）。

圖6 不同埋深條件下頁巖的游離氣特征Fig.6 Characteristicsof free gas in shale under different burialdepths

預測0到4000 m埋深范圍內，游離氣量同樣呈現先增大后減小的整體趨勢，最大游離氣含量在不同儲層壓力梯度下差異較大，分別出現在1900 m及2900m附近，而基本不受地溫梯度的影響。相同埋深條件下，儲層壓力梯度越高，地溫梯度越小，則游離氣含量越大；反之，儲層壓力梯度越低，地溫梯度越大，則游離氣含量越小。由此可見，溫度對頁巖游離氣的影響同樣是消極的，而儲層壓力梯度對頁巖游離氣的影響同樣是積極地，且對比可知，儲層壓力較儲層溫度的影響更大。因此，對于不同埋藏深度頁巖，其游離氣含量受溫度影響微弱，而主要受控于儲層壓力，且隨著埋藏深度的增大，不同溫壓條件下游離氣含量差異逐漸增大，表明儲層壓力對頁巖游離氣含量的影響隨埋藏深度的增大而逐漸增大。

3.2.5 總含氣量及吸附氣與游離氣比例

頁巖氣除少部分呈溶解狀態(tài)，因此總含氣量可近似于吸附氣與游離氣含量之和。分別將相同地溫及儲層壓力下的吸附氣及游離氣相加，即可反演不同埋深（溫壓）下頁巖總含氣量特征（圖7）。隨埋藏深度的增大，總含氣量呈現先增大后減小的整體趨勢。相同地溫梯度下，儲層壓力梯度越高，相同埋深下總含氣量越大，且隨埋深增大，儲層壓力梯度越高，總含氣量減緩趨勢越??；相同壓力梯度下，隨埋深增大，總含氣量趨勢基本相近，且儲層壓力梯度越高，總含氣量衰減幅度越小。因此，對于不同埋藏深度頁巖，其總含氣量受地溫影響小，而主要受控于儲層壓力，儲層壓力是控制總含氣量的主要因素。

利用吸附氣及游離氣與埋深的相互關系，計算吸附氣與游離氣比例見圖8。吸附氣與游離氣之比隨著埋藏深度的增大逐漸減小，表明吸附氣占總含氣量比重逐漸減小，游離氣比重逐漸增大。當儲層壓力梯度為0.8MPa/hm時，比值均大于1，表明總含氣量中吸附氣占主導，當儲層壓力梯度為1.5MPa/hm時，在埋深1700m附近，比值為1，頁巖吸附氣與游離氣量相同，可將此深度定義為頁巖氣賦存相態(tài)轉換的臨界深度。即表示，在1.5MPa/hm條件下，當頁巖埋藏深度小于1700m，吸附氣含量大于游離氣含量，頁巖氣以吸附氣占主導，當埋深大于1700m，游離氣含量逐漸超過吸附氣含量，頁巖氣以游離氣占主導。隨著儲層壓力梯度增大，臨界深度逐漸變淺，相同儲層壓力梯度條件下，地溫梯度越高，臨界深度越淺，但比較之下，溫度對比值影響較小，儲層壓力則影響較大。

圖7 不同埋深條件下頁巖總含氣量特征Fig.7 Characteristicsof totalgas contentunder different burialdepths

圖8 不同埋深條件下吸附氣與游離氣比值Fig.8 Ratiosofadsorbed gasand freegasunder differentburialdepths

圖9 WX2井龍馬溪組底部含氣量特征Fig.9 Characteristicsofgas-bearing strataat thebaseof the Longmaxi Formation inWellWX 2

3.2.6 勘探意義

為確定本次吸附氣及游離氣預測的準確性，采用WX2井鉆井巖芯實測含氣量與預測氣量進行對比。結合實測等溫吸附及覆壓孔隙度數據，模擬龍馬溪組底部吸附氣、游離氣及含氣量隨頁巖埋藏深度特征如圖9所示。模型顯示，WX2井頁巖樣品吸附氣含量介于0到3.087mL/g，最大吸附氣含量埋深約800m，游離氣含量介于0到2.25mL/g，最大游離氣含量對應埋深約2700m，總含氣量介于0到4.46mL/g，最大含氣量出現在1500m附近。吸附氣、游離氣及總含氣量隨埋深均呈現先增大后減小的趨勢。WX2井頁巖采樣深度為1623.20m，對應預測含氣量總含氣量約為4.46mL/g，其值與現場實測值誤差為8.98%（梁峰等，2016），證明本次頁巖氣賦存研究結果較為可靠。對WX2井頁巖吸附氣與游離氣比值分析可知，在埋深為2100m附近吸附氣與游離氣相等，即2100m為WX2井采樣點頁巖氣賦存相態(tài)轉換臨界深度（圖10）。推至研究區(qū)，可用已知地溫梯度和壓力梯度，計算不同埋深的溫度及儲層壓力，結合實驗測試結果，即可初步估算不同深度（溫壓）環(huán)境中各相態(tài)頁巖氣賦存特征。

WX2井頁巖樣品吸附氣含量在較小埋深即趨近于最大值，且隨著埋藏深度的增大，吸附氣含量逐漸減小，而在很小埋深時，游離氣含量較小，由此可知，頁巖氣富集需具有一定的埋藏深度。此外，由于吸附氣主要受頁巖吸附能力的控制，而游離氣受控于頁巖氣保存能力（唐相路等，2015；肖賢明等，2015），吸附氣與游離氣富集的主控因素存在差異，因此頁巖氣各賦存狀態(tài)含氣量及比例對頁巖氣富集評價具有重要意義。

綜上所述，頁巖氣賦存特征綜合受控于儲集層特征、吸附能力、地溫梯度及儲層壓力等因素。對某一研究區(qū)而言，需綜合研究頁巖吸附能力與儲集空間特征，明確含氣量與埋藏深度相關關系，并獲得最大含氣量埋深范圍，以此來確定研究區(qū)頁巖氣最宜開采深度，使頁巖氣井產能最大化。

圖10 WX2井龍馬溪組底部吸附氣與游離氣比值Fig.10 Adsorbed gasand freegas ratio at thebaseof the Longmaxi Formation inWellWX2

4 結論

（1）WX2井頁巖過剩吸附量隨壓力的增大呈先增大后減小的趨勢，最大吸附量對應壓力均在6MPa附近，伴隨著溫度的升高，最大吸附量逐漸減小，而校正后的絕對吸附量呈隨壓力的增加，先迅速增大，后增大放緩的變化趨勢，且用D-A模型來擬合絕對吸附量數據，平均誤差僅為0.604%，基本可以反頁巖真實吸附過程。

（2）WX2井頁巖樣品在加壓過程中孔隙及微裂隙會逐漸關閉，卸壓時絕大部分會重新打開，存在部分塑性變形造成的不可逆損傷，但不可逆損傷所占比重較輕。不同方向樣品孔隙度隨應力擬合曲線相關系數分別為0.9803與0.9473，指示孔隙度與有效應力之間具有負指數關系，且富含層理頁巖平行樣品較垂直樣品具有更大的初始孔隙度以及更強的孔隙應力敏感性。

（3）頁巖氣賦存特征綜合受控于儲層特征、吸附能力、儲層溫度及儲層壓力等因素，溫度對吸附氣和游離氣含量為負效應，儲層壓力為正效應，吸附氣、游離氣及總含氣量均遵循先增大后減小的總體趨勢，其中吸附氣及游離氣含量分別主要受控于溫度及儲層壓力。此外，臨界深度上下，頁巖吸附態(tài)與游離態(tài)相對含量發(fā)生變化，其對頁巖氣富集評價具有重要意義。

（

）：

曹環(huán)宇，朱傳慶，邱楠生.2016.川東地區(qū)古生界主要泥頁巖最高古溫度特征[J].地球物理學報，59(3):1017-1029.

鄧賓，劉樹根，王國芝，等.2013.四川盆地南部地區(qū)新生代隆升剝露研究——低溫熱年代學證據[J].地球物理學報，56(6):1958-1973.

董大忠，鄒才能，楊樺，等.2012.中國頁巖氣勘探開發(fā)進展與發(fā)展前景[J].石油學報，（增刊1）:107-114.

郭彤樓，張漢榮.2014.四川盆地焦石壩頁巖氣田形成與富集高產模式[J].石油勘探與開發(fā)，41(1):28-36.

郭彤樓.2016.涪陵頁巖氣田發(fā)現的啟示與思考[J].地學前緣，23(1): 29-43.

郭旭升，胡東風，李宇平，等.2016.海相和湖相頁巖氣富集機理分析與思考：以四川盆地龍馬溪組和自流井組大安寨段為例[J].地學前緣，23(2):18-28.

郭旭升.2014.南方海相頁巖氣“二元富集”規(guī)律——四川盆地及周緣龍馬溪組頁巖氣勘探實踐認識[J].地質學報，88(7):1209-1218.

何治亮，聶海寬，張鈺瑩.2016.四川盆地及其周緣奧陶系五峰組-志留系龍馬溪組頁巖氣富集主控因素分析[J].地學前緣，23(2):8-17.

紀文明，宋巖，姜振學，等.2016.四川盆地東南部龍馬溪組頁巖微—納米孔隙結構特征及控制因素[J].石油學報，37(2):182-195.

金之鈞，胡宗全，高波，等.2016.川東南地區(qū)五峰組-龍馬溪組頁巖氣富集與高產控制因素[J].地學前緣，23(1):1-10.

李傳亮，孔祥言，徐獻芝，等.2003.多孔介質的流變模型研究[J].力學學報，35(2):230-234.

李艷芳，邵德勇，呂海剛，等.2015.四川盆地五峰組—龍馬溪組海相頁巖元素地球化學特征與有機質富集的關系[J].石油學報，36(12): 1470-1483.

梁峰，拜文華，鄒才能，等.2016.渝東北地區(qū)巫溪2井頁巖氣富集模式及勘探意義[J].石油勘探與開發(fā)，43(3):1-9.

劉洪林,王紅巖,方朝合,等.2016.中國南方海相頁巖氣超壓機制及選區(qū)指標研究[J].地學前緣,23(2):48-54.

劉軍,姚艷斌,朱正杰,等.2015.渝東南地區(qū)黑色巖系儲集層特征及含氣性評價[J].礦物巖石地球化學通報,34(5):1021-1026.

劉若冰.2015.超壓對川東南地區(qū)五峰組—龍馬溪組頁巖儲層影響分析[J].沉積學報,33(4):817-827.

劉樹根,鄧賓,鐘勇,等.2016.四川盆地及周緣下古生界頁巖氣深埋藏-強改造獨特地質作用[J].地學前緣,23(1):11-28.

聶海寬,包書景,高波,等.2012.四川盆地及其周緣上奧陶統—下志留統頁巖氣成藏體系研究[J].石油實驗地質,34(2):115-119.

湯慶艷,張銘杰,余明,等.2013.頁巖氣形成機制的生烴熱模擬研究[J].煤炭學報,38(5):742-747.

唐相路,姜振學,張鶯鶯,等.2015.渝東南地區(qū)頁巖氣富集區(qū)差異性分布成因[J].西安石油大學學報(自然科學版),30(3):24-30.

魏志紅.2015.四川盆地及其周緣五峰組-龍馬溪組頁巖氣的晚期逸散[J].石油與天然氣地質,36(4):659-665.

夏威,于炳松,孫夢迪.2015.渝東南YK1井下寒武統牛蹄塘組底部黑色頁巖沉積環(huán)境及有機質富集機制[J].礦物巖石,35(2):70-80.

夏陽,金衍,陳勉,等.2015.頁巖氣滲流數學模型[J].科學通報,60(24): 2259-2271.

肖賢明,王茂林,魏強,等.2015.中國南方下古生界頁巖氣遠景區(qū)評價[J].天然氣地球科學,26(8):1433-1445.

熊小輝,王劍,余謙,等.2015.富有機質黑色頁巖形成環(huán)境及背景的元素地球化學反演——以渝東北地區(qū)田壩剖面五峰組—龍馬溪組頁巖為例[J].天然氣工業(yè),35(4):25-32.

郇璇,張小兵,韋歡文.2015.基于不同類型煤吸附甲烷的吸附勢重要參數探討[J].煤炭學報,40(8):1859-1864.

游聲剛,郭茜,耿小燼,等.2015.頁巖含氣量的影響因素分析及含氣量測試方法[J].中國礦業(yè),24(12):80-85.

張金川,金之鈞,袁明生.2004.頁巖氣成藏機理和分布[J].天然氣工業(yè), 24(7):15-18.

張雪芬,陸現彩,張林曄,等.2010.頁巖氣的賦存形式研究及其石油地質意義[J].地球科學進展,25(6):597-604.

朱炎銘，陳尚斌，方俊華，等.2010.四川地區(qū)志留系頁巖氣成藏的地質背景[J].煤炭學報，35(7):1160-1164.

鄒才能，董大忠，王社教，等.2010.中國頁巖氣形成機理、地質特征及資源潛力[J].石油勘探與開發(fā)，37(6):641-653.

鄒才能，董大忠，王玉滿，等.2015.中國頁巖氣特征、挑戰(zhàn)及前景(一)[J].石油勘探與開發(fā)，42(6):689-701.

鄒才能，董大忠，王玉滿，等.2016.中國頁巖氣特征、挑戰(zhàn)及前景(二)[J].石油勘探與開發(fā)，43(2):166-178.

耶格JC和庫克NGW.1981.巖石力學基礎[M].北京:科學出版社：268-272.

Athy L F.1930.Density,porosity,and compaction of sedimentary rocks[J]. AAPG Bulletin,14(1):1-24.

Curtis JB.2002.Fractured shale-gas systems[J].AAPG bulletin,86(11): 1921-1938.

Giles M R,Indrelid S L and James D M D.1998.Compaction the great unknown in basin modeling[J].London:Geological Society,Special Publications,141(1):15-43.

HarpalaniS,Prusty BK and Dutta P.2006.Methane/CO2sorptionmodeling for coalbed methane production and CO2sequestration[J].Energy and Fuels,20(4):1591-1599.

JIW M,Song Y,Jiang ZX,et al.2015.Estimation ofmarine shalemethane adsorption capacity based on experimental investigations of Lower Silurian Longmaxi formation in the Upper Yangtze Platform,south China[J].Marineand Petroleum Geology,68:94-106.

Milliken K L,RudnickiM,Awwiller DN,etal.2013.Organicmatterehosted pore system,Marcellus Formation(Devonian),Pennsylvania[J].AAPG Bull,97:177-200.

Singh H and Javadpour F.2016.Langmuir slip-Langmuir sorption permeabilitymodelofshale[J].Fuel,164:28-37.

Wang Y,Zhu Y M,Liu SM,et al.2016.Methane adsorption measurements and modeling for organic-rich marine shale samples[J].Fuel,172: 301-309.

ShaleGasOccurrenceCharacteristicsand Exploration Significanceof WX2Wellin theNortheastChongqing

FU Changqing1,2,ZHU Yanm ing1,2,CHEN Shangbin1,2,LIANG Feng2,3
1.Key Laboratory ofCoalbed Methane Resources&Reservoir Formation Process,Ministry ofEducation, China University ofMining&Technology,Xuzhou 221116,China; 2.SchoolofResources&Geosciences,China University ofMining&Technology,Xuzhou 221008,China; 3.PetroChina Research Instituteof Petroleum Exploration and Development-Langfang Branch,Langfang065007,China

To study the shale gas occurrence characteristics of the Longmaxi Formation in northeast Chongqing,the isothermal adsorption data under high temperature and high pressure and porosity data under burden pressure of the well WX2 shale were analyzed.According to the principle ofminimum error,we have chosen a suitable adsorptionmodel for the study area.Based on the relationship between porosity and effective stress,the free gas model is established.In addition,comprehensive analysis of the characteristics of adsorption gas,free gas,and changes in total gaswith the buried depth was conducted.The results show thatexcess adsorption quantity ofwellWX2 shale increased firstand then decreased under different temperatureswith an increasing pressure.Themaximum adsorption capacity decreased gradually as the temperature increased.After correction,the absolute adsorption quantity shows a rapid increase firstand then a slow increasewith increasing pressures.As the absoluteadsorption data can be fitted with D-A modelwith theminimum average error,itcan represent the realadsorption processofshale in the study area.Pores andmicro-cracksof shale samples will gradually close during pressurization,and reopen during depressurization.This process will cause irreversible damage by plastic deformation,but irreversible damage accounts for small proportion.Porosities of samples from different directions have a negative exponent relationship with effective stress.Parallel samples of lamina-rich shales have greater initial porosity and stronger pore stress sensitivity than the vertical samples.Shale gas occurrence characteristic is comprehensively controlled by reservoir characteristics,adsorption capacity,reservoir temperature,and reservoir pressure.Among them,temperature has the negativeeffectson adsorption gas content and free gas content,and the reservoir pressure has the positive effect.The adsorption gas content,free gas content and total content increases first then decreases.Adsorption gas content and free gas content are mainly controlled by temperature and reservoir pressure,respectively.In addition,relative abundance of adsorption gas and free gaswill change near the criticaldepth,and thiswillbeofgreat importance forevaluating shalegasenrichment.

WellWX2;occurrence characteristics;gas content;criticaldepth;northeastChongqing

FUChangqing,Ph.D.Candidate;E-mail:spirit1208@163.com

P618.1

A文獻標識碼：1006-7493（2016）04-0679-11

10.16108/j.issn1006-7493.2016076

2016-05-31；

2016-08-20

煤層氣資源與成藏過程教育部重點實驗室（中國礦業(yè)大學）開放基金資助項目（2016-005）；國家自然科學基金項目（41272155）；國家基礎研究發(fā)展計劃（973）項目（2012CB214702）聯合資助

付常青，男，1990年生，博士研究生，從事非常規(guī)油氣地質研究；E-mail：spirit1208@163.com