何梅興，裴發(fā)根,2，方慧，張鵬輝，何大雙，張小博，呂琴音，仇根根

(1.中國地質(zhì)科學(xué)院地球物理地球化學(xué)勘查研究所，河北廊坊，065000；2.成都理工大學(xué)地球物理學(xué)院，四川成都，610059)

我國祁連山凍土區(qū)與極地地區(qū)的天然氣水合物在成藏條件和分布特征方面存在很大差異，如加拿大馬更些三角洲地區(qū)(Mackenzie Delta)地溫梯度小，凍土層厚度較大，天然氣水合物主要賦存在古近系碎屑沉積層內(nèi)，氣體來源為熱成因氣，與常規(guī)油氣資源共生，水合物飽和度高[1-3]。而祁連山地區(qū)屬中緯度高海拔凍土區(qū)，地溫梯度大，凍土層相對(duì)薄，水合物主要賦存在中侏羅統(tǒng)江倉組，地層產(chǎn)狀陡、斷裂發(fā)育，水合物埋深淺、飽和度低[4]。相比極地地區(qū)，祁連山凍土區(qū)的天然氣水合物成藏條件較差，礦體規(guī)模小，成藏規(guī)律更加復(fù)雜。在木里聚乎更煤礦區(qū)三露天井田內(nèi)26 口天然氣水合物鉆探井中，其中11 口鉆井鉆獲的天然氣水合物呈不連續(xù)分布，規(guī)律不明顯[5-7]，這導(dǎo)致了水合物資源勘探難度大，亟待加強(qiáng)研究該區(qū)有關(guān)天然氣水合物不均勻性分布的地質(zhì)控制因素等基礎(chǔ)性問題。

前人從地層巖性、氣源類型、斷層破碎帶和凍土層屬性等角度探討分析了與天然氣水合物分布的規(guī)律，認(rèn)為斷層是最主要控制因素，氣源成因類型、凍土層的構(gòu)成及巖性對(duì)水合物均有一定影響[5,8]。但這些分析主要基于有限鉆井巖心的測試結(jié)果，木里地區(qū)地質(zhì)情況復(fù)雜，鉆井資料有限，還缺乏天然氣水合物與其成藏控制因素的整體認(rèn)識(shí)，而且涉及的地球物理資料有限，因此，對(duì)天然氣水合物空間分布與主要控制因素的對(duì)應(yīng)關(guān)系認(rèn)識(shí)還不足。

針對(duì)天然氣水合物分布規(guī)律不明顯與地質(zhì)控制因素不清等問題，在以往工作基礎(chǔ)上，從木里三露天的電性結(jié)構(gòu)與天然氣水合物分布關(guān)系的角度，綜合分析音頻大地電磁測深結(jié)果，探討該區(qū)天然氣水合物運(yùn)移通道和凍土層分布等對(duì)天然氣水合物分布的控制作用，為進(jìn)一步研究青藏高原凍土區(qū)天然氣水合物成藏機(jī)制和下一步勘探開發(fā)提供借鑒與參考。

1 地質(zhì)概況及地球物理場特征

圖1所示為祁連造山帶構(gòu)造單元分布。由圖1可見：研究區(qū)處于中祁連構(gòu)造帶中西部，石炭紀(jì)—二疊紀(jì)，區(qū)域接受了一套海相和海陸交互相沉積，三疊紀(jì)受印支運(yùn)動(dòng)影響，地層隆升剝蝕，侏羅紀(jì)的燕山運(yùn)動(dòng)導(dǎo)致局部拉伸，沉積了一套煤系地層。晚中生代的燕山運(yùn)動(dòng)構(gòu)造抬升，沉積面積萎縮，白堊系被剝蝕而保留了部分侏羅統(tǒng)。自新近紀(jì)以來，受喜馬拉雅運(yùn)動(dòng)的影響，祁連山地區(qū)處于強(qiáng)烈擠壓構(gòu)造應(yīng)力區(qū)，發(fā)生了大幅度的地殼縮短、抬升和旋扭。天然氣水合物發(fā)現(xiàn)區(qū)木里坳陷聚乎更礦區(qū)主要由2個(gè)向斜和1個(gè)背斜構(gòu)成的復(fù)向斜構(gòu)造組成[9-11]。

圖1 祁連造山帶構(gòu)造單元分布Fig.1 Distribution of tectonic units in Qilian orogenic belt

根據(jù)祁連山重力場特征，研究區(qū)可劃分為北祁連重力高值梯度帶、中祁連重力場過渡帶和南祁連重力場平緩區(qū)，前2個(gè)區(qū)總體呈北西—南東走向重力梯度帶，祁連山西南部重力低反映了地殼增厚，東北至西南莫霍面逐漸加深，最深處位于祁連山西南部[12]。航磁異常分布特征顯示祁連山東、西兩段基底巖系出露區(qū)磁異常強(qiáng)度為50～250 nT，中部晚古生代—早中生代分布區(qū)磁異常強(qiáng)度為-25～-75 nT，磁異常強(qiáng)度差異性表明中部基底埋深大，可能存在北東方向的基底隱伏斷裂[13]。

2 工作方法

在天然氣水合物地球物理調(diào)查中，電磁法勘探是一種重要的方法手段，該方法根據(jù)水合物儲(chǔ)層與圍巖之間的電阻率差異，推斷天然氣水合物賦存位置[14]。在青藏高原地球物理方法有效性對(duì)比實(shí)驗(yàn)中，該方法對(duì)于天然氣水合物形成、運(yùn)移和封堵所需要的斷裂構(gòu)造和凍土蓋層也有較好的探測效果[15]。

2.1 數(shù)據(jù)采集

測區(qū)位于木里坳陷聚乎更礦區(qū)三露天井田，測線穿過多個(gè)天然氣水合物鉆井，完全覆蓋了已發(fā)現(xiàn)水合物區(qū)域。測線南端跨越了大通山北緣主控?cái)嗔?，測線方向垂直于北西-南東向的背斜構(gòu)造，如圖2所示。音頻大地電磁測深(AMT)采用五分量張量測量方式，即同時(shí)測量天然電磁場5個(gè)分量(電場分量Ex和Ey；磁場分量Hx，Hy和Hz；其中，x，y和z為觀測坐標(biāo)系統(tǒng)方向)，布設(shè)采用正“十”形測量裝置。

圖2 木里坳陷三露天音頻大地電磁測點(diǎn)與地層分布圖Fig.2 AMT measuring points and stratigraphic distribution

2.2 數(shù)據(jù)處理

對(duì)采集的音頻大地電磁測深電場和磁場時(shí)間序列進(jìn)行頻譜分析，采用Robust 方法估算了實(shí)測坐標(biāo)系的阻抗張量元素，為了更好地對(duì)比觀測結(jié)果，進(jìn)行了二維和三維反演。由于大地電磁信號(hào)信噪比低和信號(hào)干擾，導(dǎo)致某些頻點(diǎn)的數(shù)據(jù)誤差較大，甚至一些頻點(diǎn)或頻段發(fā)生了畸變，在數(shù)據(jù)反演之前先進(jìn)行主軸旋轉(zhuǎn)、測線圓滑和靜態(tài)校正等預(yù)處理。結(jié)合地質(zhì)情況，選定TM模式二維反演結(jié)果作為剖面解釋依據(jù)。

三維反演采用SIRIPUNVARAPORN 教授在Occam反演算法基礎(chǔ)上改進(jìn)的數(shù)據(jù)空間算法[16]，該算法將模型空間轉(zhuǎn)換到數(shù)據(jù)空間，數(shù)據(jù)遠(yuǎn)小于模型參數(shù)數(shù)據(jù)，極大地減小了反演的計(jì)算量。三維反演采用測區(qū)完成的AMT 測量189 個(gè)測點(diǎn)，點(diǎn)距為100 m，面積約為6 km2，三維模型網(wǎng)格x，y和z軸分別剖分為45，51和25個(gè)單元格，邊界網(wǎng)格數(shù)為5 個(gè)，水平方向單元網(wǎng)格長×寬(x×y)為100 m×100 m，深度(z)單元網(wǎng)格首層為20 m，以1.05倍數(shù)遞增，有效反演深度為1 000 m，反演迭代10 次，反演擬合誤差為1.53%。

3 研究區(qū)電性特征

3.1 地層電性特征

祁連山木里地區(qū)中侏羅統(tǒng)江倉組下段的下部及其下伏的木里組均為相對(duì)高阻的含煤地層，砂巖、粉砂巖為次高阻，中侏羅統(tǒng)江倉組下段的上部泥巖地層為低阻特征，上三疊系砂巖、含煤層為中高阻特征。含天然氣水合物地層段均表現(xiàn)出較明顯高電阻率特征，凍土層和水合物賦存帶與相同巖性非水合物賦存帶相比表現(xiàn)出相對(duì)高阻特征，各地層的電阻率差異為電磁法試驗(yàn)提供了物性條件。

3.2 研究區(qū)整體電性結(jié)構(gòu)

圖3所示為木里三露天三維電性結(jié)構(gòu)模型。由圖3可見，高阻異常區(qū)對(duì)應(yīng)為三疊系和木里組含煤地層，主要分布在水合物區(qū)西南部和東北部，地層巖性相對(duì)完整。低阻異常主要對(duì)應(yīng)侏羅系泥巖、油頁巖，范圍主要分布在中東部，含泥巖、頁巖的侏羅系和斷裂構(gòu)造較為發(fā)育。整體電性特征與地質(zhì)構(gòu)造特征具有較好對(duì)應(yīng)，表現(xiàn)為“東西向分塊、南北向分帶”的特征。

圖3 木里三露天三維電性結(jié)構(gòu)模型Fig.3 Three-dimensional electrical structure model in Sanlutian of Muli area

4 天然氣水合物地質(zhì)控制因素

4.1 氣源

木里凍土區(qū)天然氣水合物以熱解成因氣為主，微生物成因氣較少[17]，但對(duì)熱解成因中是油型氣還是煤型氣還存在分歧。從天然氣水合物氣體組成和同位素基本特征分析，烴類氣體主要是油型氣[18]。但如果根據(jù)天然氣水合物在煤系地層的賦存狀態(tài)、碳同位素和氣體組成，則指示為煤成氣[19]。氣源對(duì)比表明中侏羅紀(jì)泥巖暗色泥頁巖和晚三疊世尕勒得寺組暗色泥巖是主要的烴源巖，中侏羅統(tǒng)與上三疊統(tǒng)烴源巖是今后的主要目標(biāo)層[20]。在氣源分布差異方面，三露天的鉆井中巖心顯示，烴類氣體有效供應(yīng)量從西至東逐次變差[5]，氣源的有效供應(yīng)狀況可能決定了天然氣水合物分布的差異性。

4.2 斷裂構(gòu)造

燕山期木里地區(qū)受區(qū)域擠壓應(yīng)力作用，主要斷裂性質(zhì)由張性轉(zhuǎn)變?yōu)閴盒?，南?cè)發(fā)育規(guī)模較大的逆沖推覆斷裂F1和F2[21-23]，控制著坳陷南邊界。在電阻率剖面上反映了逆沖推覆構(gòu)造的電性特征，表現(xiàn)為低阻異常和高阻異常傾向一致向西南傾斜，反映了斷層向北逆沖，傾向?yàn)槲髂?，根?jù)低阻異常特征反映的逆沖斷裂，可以劃分出主斷裂F1 和F2，目前發(fā)現(xiàn)的水合物主要分布主斷裂附近，構(gòu)造控礦明顯。

圖4所示為L01 與L02 測線電阻率剖面的低阻異常帶電性特征。由圖4可見，從西至東F1 和F2推覆斷裂產(chǎn)狀及規(guī)模不同，總體特征為發(fā)育規(guī)模變小，傾角變陡，推覆距離變小。推覆斷裂的發(fā)育規(guī)模、傾角和推覆距離對(duì)水合物成藏都有控制作用，推覆斷裂弧形前鋒附近斷裂最為發(fā)育，F(xiàn)1和F2 組成的斷距最大，該處還發(fā)育北東向平移斷裂，形成了該區(qū)斷裂構(gòu)造發(fā)育區(qū)，傾角較東部小，推覆距離較遠(yuǎn)，鉆孔發(fā)現(xiàn)的水合物最多集中在該區(qū)，說明該處相對(duì)低角度的推覆斷裂起著封堵作用，為烴源巖地層提供了保存條件，烴類氣體最富集，同時(shí)提供了運(yùn)移和成藏空間，在溫壓條件下，形成了天然氣水合物分布集中區(qū)。因此，破碎帶發(fā)育的規(guī)模大、傾角小和推覆距離遠(yuǎn)的推覆斷裂前鋒區(qū)域?qū)λ衔锍刹馗欣?/p>

圖4 木里三露天AMT測線電性斷面及綜合地質(zhì)解釋圖Fig.4 Resistivity profile and comprehensive interpretation of AMT survey line in Sanlutian of Muli area

圖5所示為三露天測區(qū)電阻率深度切片圖。由圖5可見：測區(qū)區(qū)西南部存在1 條明顯的北西-南東走向的低阻異常帶，與大通山北緣斷裂構(gòu)造方向一致，根據(jù)北西-南東方向的低阻帶可推測斷裂發(fā)育情況。從淺部(埋深為-187 m)到深部(埋深為-853 m)，低阻異常發(fā)育范圍逐漸擴(kuò)大，表明了深部斷裂發(fā)育規(guī)模逐漸變大，推測斷裂構(gòu)造發(fā)育深度大于1 000 m。從深部(埋深為-583 m)到淺部斷裂走向上，存在局部高阻阻隔區(qū)，破碎帶連通性逐漸變差，表明了木里坳陷南部主斷裂在走向與垂向上分布的差異性。

由圖5可見，淺部方向約為北偏西45°，深部方向約為北偏西60°，說明主斷裂深淺存在角度差，表明在推覆作用下，主斷裂從深部北西西向旋轉(zhuǎn)至淺部北西向。從地表地質(zhì)情況可知，DK3處為弧形推覆構(gòu)造的前鋒，測區(qū)處于弧形推覆構(gòu)造的東側(cè)，結(jié)合電性特征和推覆作用差異性，主斷裂構(gòu)造附近存在2 處局部低阻異常區(qū)，推斷為2個(gè)斷裂破碎通道，分別在DK3和DK8-19附近，其中DK3 附近范圍內(nèi)發(fā)育規(guī)模較大。斷裂通道與水合物發(fā)現(xiàn)區(qū)有較好的對(duì)應(yīng)性，說明斷裂通道對(duì)水合物分布起了重要控制作用。在測區(qū)中部(DK10-16和DK10-18附近)區(qū)域，中淺部對(duì)應(yīng)為中高阻地層，缺乏貫穿深部的斷裂通道，應(yīng)是未鉆遇水合物的原因之一。

圖5 三露天測區(qū)電阻率深度切片圖Fig.5 Resistivity depth slice diagram in Sanlutian of Muli area

為分析斷裂破碎帶的連通情況，建立了電性結(jié)構(gòu)體小于30 Ω·m的三維等值面圖，如圖6所示。由圖6可見：測區(qū)存在破碎帶構(gòu)成的2個(gè)“氣煙囪”通道，通道位置及鉆遇水合物井的分布有較好對(duì)應(yīng)性。測區(qū)也發(fā)育泥頁巖低阻層，根據(jù)“氣煙囪”通道形態(tài)可推測低阻通道應(yīng)該是低阻泥巖層和斷裂破碎帶兩者的綜合反映?！皻鉄焽琛蓖ǖ莱氏聦捝险螒B(tài)，通道延伸至深部，從深及淺發(fā)育2 個(gè)“氣煙囪”式的斷裂通道，該通道與已發(fā)現(xiàn)的水合物鉆孔對(duì)應(yīng)性較好，說明了烴類氣源沿著“氣煙囪”通道運(yùn)移并成藏。

圖6 三露天測區(qū)三維電阻率等值面與氣源通道Fig.6 Three-dimensional electrical structure model and source gas channel in Sanlutian of Muli area

4.3 凍土層

凍土層具有致密性高和通透性差的特點(diǎn)，可阻止烴類氣體向上遷移至空氣中，穩(wěn)定分布在地表的凍土層。木里地區(qū)海拔一般為4 000～4 300 m，年平均地表地溫最低為-2.4 ℃，地表凍土發(fā)育。祁連山木里地區(qū)具有豐富烴類氣體，氣源條件良好，凍土條件已成為天然氣水合物成藏的關(guān)鍵因素之一[24]。

由于凍土層含冰量高，具有地表上連續(xù)或間斷分布的高電阻層特征，應(yīng)用音頻大地電磁測深探查凍土層應(yīng)用效果良好。依據(jù)音頻大地電磁測深獲取的天然氣水合物底板埋深及厚度，與鉆井實(shí)測的凍土層厚度和底板埋深基本一致，結(jié)合音頻大地電磁測深和高密度電法結(jié)果表明該區(qū)凍土較發(fā)育，平均厚度約82 m。

圖7所示為三露天實(shí)驗(yàn)區(qū)永凍土底界厚度分布圖。由圖7可見：大通山北坡水合物分布區(qū)DK13-11 至DK10-17 范圍凍土層相對(duì)發(fā)育，DK10-16 至DK10-18的東部凍土層發(fā)育差?？傮w上三露天凍土以島狀形態(tài)分布，西南部發(fā)育，東部發(fā)育差。受地形影響，地表水流向具有方向性，大通山北坡向三露天的地表水流向大致向DK11-14 至DK-10區(qū)域，地表水聚集區(qū)凍土層發(fā)育，可能是造成該區(qū)凍土層發(fā)育主要因素。除此之外，凍土層發(fā)育還受高程、地形和巖性等因素影響，三露天凍土層發(fā)育的不均勻性，在一定程度上限制了水合物成藏的有利空間。

圖7 三露天實(shí)驗(yàn)區(qū)永凍土底界厚度分布圖Fig.7 Distribution of thickness of permafrost in Sanlutian of Muli area

5 成藏模式探討

天然氣水合物主要在氣源、構(gòu)造和凍土層之間耦合機(jī)制體系下成藏[25-27]，水合物鉆探顯示天然氣水合物產(chǎn)出層段多為零星、空間不連續(xù)分布，天然氣水合物成藏機(jī)理復(fù)雜。侏羅系下伏的烴源巖的生排烴時(shí)間在距今90～100 Ma，而距今1.2 Ma才開始發(fā)育高原凍土，天然氣水合物形成于距今1.2 Ma之后[20]，這期間喜山運(yùn)動(dòng)導(dǎo)致祁連山地區(qū)遭受了強(qiáng)烈的地殼壓縮、隆起和剝蝕等，地層改造強(qiáng)烈，油氣保存條件相對(duì)差，烴源巖的完整性和圈閉構(gòu)造是水合物成藏的關(guān)鍵因素。

依據(jù)低阻異常范圍的深度，大通山北緣斷層發(fā)育深度大于1 km，側(cè)面說明了氣源來源于中深部，推測來自于約1 km 深度或更深處。在三露天天然氣水合物測區(qū)斷層連通性方面，斷層已連通深部至淺部，存在2個(gè)“氣煙囪”通道，水合物發(fā)現(xiàn)區(qū)域與斷層通道分布吻合，遠(yuǎn)離“氣煙囪”通道的水合物鉆井未鉆遇水合物。如DK6 和DK9 鉆井對(duì)比(圖8)，地層巖性都含泥巖、油頁巖等烴源巖層，應(yīng)該都具備了生成天然氣水合物的氣源條件，DK9 鉆井位于“氣煙囪”通道上方，凍土層厚度約100 m，而DK6遠(yuǎn)離“氣煙囪”通道，凍土層厚度只有65 m，DK9 在氣源通道和凍土層封閉條件都比DK6 有利，DK9 鉆井鉆遇了3 層含水合物層，表明水合物成藏與“氣煙囪”通道及凍土層分布密切有關(guān)，缺乏氣源通道運(yùn)移和凍土封閉條件是造成水合物不均勻性分布的主要因素。

圖8 水合物鉆井測井曲線與水合物分布圖Fig.8 Hydrate well logging and hydrate distribution

圖9所示為祁連山木里凍土區(qū)天然氣水合物成藏模式。木里地區(qū)斷陷接受侏羅系沉積后，受燕山后期區(qū)域壓應(yīng)力場作用下，中生代地層開始褶皺變形，隨著喜山期擠壓作用的加劇，坳陷南北兩側(cè)發(fā)育對(duì)沖的逆沖斷層，形成了褶—斷組合形態(tài)，在三露天的南部大通山北緣發(fā)育了北西西走向的疊瓦狀由南向北逆沖斷層，傾向南西。斷層逆推過程中，因侏羅系中發(fā)育低強(qiáng)度、高塑性的煤和泥巖層等軟弱層位，三疊系沿著軟弱滑脫層發(fā)生逆沖推覆。從三露天地層分布可知，受局部應(yīng)力和軟弱層位分布差異，逆沖斷層平面分布形態(tài)呈弧形，弧形前鋒推覆距離至少1 km，天然氣水合物主要位于應(yīng)力集中的逆沖斷層前鋒帶下盤，擠壓作用下形成大量逆沖斷層，易形成巖性圈閉，三露天天然氣水合物控礦構(gòu)造樣式以逆沖前鋒型為主。

圖9 祁連山木里凍土區(qū)天然氣水合物成藏模式Fig.9 Pattern for gas hydrate accumulation Muli area in Qilian mountain

6 結(jié)論

1)木里三露天測區(qū)發(fā)育2 個(gè)“氣煙囪”通道，主要由斷裂破碎帶構(gòu)成，該通道呈“下寬上窄”形態(tài)，淺部斷裂通道與深部斷裂連通，該通道與已發(fā)現(xiàn)的天然氣水合物鉆孔對(duì)應(yīng)性較好，指示烴類氣源沿著“氣煙囪”式的斷裂通道運(yùn)移并成藏。

2)據(jù)音頻大地電磁法反演結(jié)果得出木里三露天測區(qū)整體三維電性結(jié)構(gòu)模型，測區(qū)整體電性特征以“東西向分塊、南北向分帶”為主，反映了地層及地質(zhì)構(gòu)造特征。

3)木里三露天天然氣水合物主要位于逆沖斷層前鋒帶下盤，控礦構(gòu)造樣式以逆沖前鋒型為主，木里三露天天然氣水合物空間分布不連續(xù)，斷裂通道及凍土層共同作用下可形成水合物成藏有利條件，天然氣水合物分布具有“逆沖推覆前鋒帶、斷裂破碎帶和凍土層發(fā)育區(qū)”的規(guī)律。