馬嶺油田長7致密油儲層裂縫特征及分布規(guī)律

趙國璽　劉浩宇　張瀚丹　何小娟

（1．中國石油長慶油田分公司勘探開發(fā)研究院　陜西　西安　710018；2.低滲透油氣田勘探開發(fā)國家工程實驗室　陜西　西安　710018；3.長慶油田第四采油廠　陜西　靖邊　718500）

?

馬嶺油田長7致密油儲層裂縫特征及分布規(guī)律

趙國璽1,2劉浩宇3張瀚丹1,2何小娟1,2

（1．中國石油長慶油田分公司勘探開發(fā)研究院陜西西安710018；2.低滲透油氣田勘探開發(fā)國家工程實驗室陜西西安710018；3.長慶油田第四采油廠陜西靖邊718500）

摘要本文通過地質(zhì)露頭及巖心觀察、巖石聲發(fā)射實驗、古地磁裂縫定向，巖石壓縮實驗、差應(yīng)變實驗、波速各向異性實驗等方法,對馬嶺油田長7儲層裂縫開展定性和定量研究。確定馬嶺長7儲層主要發(fā)育層理縫、張性縫、壓扭縫、張扭縫。裂縫發(fā)育優(yōu)勢方位為NEE向和NNW向。應(yīng)用數(shù)值模擬方法開展裂縫定量預(yù)測，劃分了裂縫相對發(fā)育區(qū)、次發(fā)育區(qū)和不發(fā)育區(qū),為今后馬嶺油田長7致密油制定規(guī)模開發(fā)技術(shù)政策提供了地質(zhì)依據(jù)。

關(guān)鍵詞馬嶺油田；長7致密油；裂縫特征；裂縫分布規(guī)律

1　引言

近年來致密油逐漸成為全世界非常規(guī)油氣勘探開發(fā)的新熱點[1～4]。美國、加拿大、澳大利亞等國家都已經(jīng)對致密油進行了大規(guī)模商業(yè)開發(fā)[5～7]。中國的致密油氣資源很豐富，目前還處于堪探開發(fā)試驗階段[8～9]。根據(jù)勘探開發(fā)實際，長慶油田將空氣滲透率小于0.3mD的儲層稱為致密儲層，賦存于其中的石油稱為致密油，主要分布于鄂爾多斯盆地長7儲層[10]。

馬嶺油田長7儲層空氣滲透率0.19mD，屬于典型的致密油儲層，室內(nèi)試驗和現(xiàn)場實踐證實馬嶺油田長7儲層孔喉細微，物性差，但天然裂縫普遍較發(fā)育。天然裂縫發(fā)育特征及分布規(guī)律是影響單井產(chǎn)量的主要因素之一。天然裂縫特征及分布規(guī)律的研究對于認識油氣成藏機理、制定開發(fā)技術(shù)政策及儲層壓裂工藝等方面均具有重要意義。

2　長7致密油儲層裂縫特征

2.1裂縫發(fā)育特征參數(shù)

根據(jù)野外露頭及室內(nèi)巖心觀察，馬嶺油田長7儲層裂縫根據(jù)產(chǎn)狀和形成機理可劃分為4種類型。包括高角度張扭縫、壓扭縫，低角度張性縫、層理縫和滑脫縫。

研究區(qū)46口井長7目的層段開展巖心觀察，其中41口井觀察到不同程度的宏觀裂縫發(fā)育，裂縫鉆遇率達89.1%。觀測到的各類裂縫總計1080條。其中，高角度壓扭縫鉆遇率最高為60.87%，張扭縫鉆遇率為34.78%，低角度張性縫鉆遇率為43.48%，層理縫鉆遇率為36.96%（圖1）。

圖1　馬嶺油田長7儲層裂縫類型分布圖

在高角度裂縫中有75%被方解石充填。所有46口井巖心裂縫中，有39口井發(fā)現(xiàn)油跡，所占比例為84.8%。所有類型的裂縫均有可能成為油氣運移通道。

2.2儲層構(gòu)造裂縫優(yōu)勢方位

圖2　馬嶺油田長7儲層巖石密度等值線

鄂爾多斯盆地中生界構(gòu)造裂縫主要形成于晚中生代燕山期構(gòu)造運動[11]。根據(jù)野外露頭觀察，延長組長7段砂巖構(gòu)造裂縫大多成組出現(xiàn)，主要為未穿層裂縫，共軛裂縫近于垂直地層層面分布，主要發(fā)育NEE向和NNW向兩組裂縫，次發(fā)育NE向裂縫，裂縫產(chǎn)狀穩(wěn)定，延伸較遠，為區(qū)域性裂縫。

表1　馬嶺油田長7儲層力學(xué)模型巖石力學(xué)參數(shù)選擇

圖3　格里非斯張破裂擴展方向示意圖

圖4　巖石剪切破裂分析圖（a-莫爾圓；b-剪切破裂）

圖5　馬嶺油田長7裂縫有利區(qū)預(yù)測圖

根據(jù)馬嶺油田22口井長7儲層巖芯裂縫古地磁定向結(jié)果，壓扭縫走向分為NEE (60.6°～79.9°)和NWW(96.7°～123°)兩組，具有較為明顯的共軛特征；張扭縫主要為近EW向（81°～88.5°），近似為共軛壓扭縫的平分線。與相似露頭區(qū)裂縫系統(tǒng)優(yōu)勢方位一致。

3　構(gòu)造裂縫分布定量預(yù)測

致密油儲層裂縫定量預(yù)測的關(guān)鍵是地質(zhì)模型、力學(xué)模型和數(shù)學(xué)模型的建立[12]。

3.1地質(zhì)模型建立

儲層巖石密度能較好的反映儲層的非均質(zhì)性。因此，模擬對象的地質(zhì)體用不同的巖石密度區(qū)帶來表示，并對不同的巖石密度區(qū)帶賦予不同的巖石力學(xué)參數(shù)。根據(jù)馬嶺油田130余口井密度解釋成果，建立了長7儲層構(gòu)造裂縫預(yù)測地質(zhì)模型（圖2）。

3.2力學(xué)模型建立

力學(xué)模型是將地質(zhì)模型轉(zhuǎn)化為數(shù)值計算的關(guān)鍵，包括巖石力學(xué)參數(shù)的綜合選取，邊界力的作用方式、方向及邊界約束條件等[13]。

巖石力學(xué)參數(shù)的選取是根據(jù)馬嶺油田不同巖性巖石力學(xué)性質(zhì)測量數(shù)據(jù)，萃取不同地區(qū)泥巖、砂巖端元參數(shù)，包括彈性模量、泊松比、巖石抗拉強度、內(nèi)聚力和內(nèi)摩擦角等參數(shù)。依據(jù)研究區(qū)地質(zhì)模型巖石密度，參考實測數(shù)據(jù)，賦予不同地區(qū)不同巖石物理力學(xué)屬性(表1)。

研究區(qū)主要形變期燕山期應(yīng)力場最大主壓方位由地層共軛配套節(jié)理和區(qū)域地層形變規(guī)律而得來。馬嶺油田燕山期水平投影方位為84°；相應(yīng)最小主應(yīng)力傾伏角近似于水平，其方位為174°，此次模擬σ1方向取值為84°。

根據(jù)研究區(qū)巖石差應(yīng)變實驗結(jié)果，現(xiàn)今水平最大有效主應(yīng)力方向約為76°，水平最小有效主應(yīng)力方向約為166°。

根據(jù)聲發(fā)射法應(yīng)力測試結(jié)果，裂縫發(fā)育期的水平最大有效主應(yīng)力約44MPa，水平最小有效主應(yīng)力約24MPa?，F(xiàn)今水平最大有效主應(yīng)力約27MPa，水平最小有效主應(yīng)力約15Mpa。

在選取巖石力學(xué)參數(shù)確定應(yīng)力場邊界條件的基礎(chǔ)上，建立了馬嶺油田長7致密油儲層的力學(xué)模型。

3.3數(shù)學(xué)模型建立

根據(jù)研究區(qū)長7儲層巖心裂縫觀察的力學(xué)性質(zhì)既有張性，又有剪性，故有必要對這兩種裂縫發(fā)育和分布分別進行預(yù)測研究。在力學(xué)模型的基礎(chǔ)上，利用巖石破裂法[8]建立數(shù)學(xué)模型，定量預(yù)測構(gòu)造裂縫分布規(guī)律。按照庫倫-摩爾準(zhǔn)則[14]和格里菲斯準(zhǔn)則[15]分別計算出剪破裂系數(shù)（R）、張破裂系數(shù)(η)。

（1）格里菲斯張破裂準(zhǔn)則

根據(jù)格里菲斯張破裂準(zhǔn)則,脆性物體的破壞是由內(nèi)部存在裂隙決定的。巖石內(nèi)存在很多微小裂隙，在裂隙的尖端存在應(yīng)力集中，使裂隙擴展，以至破壞形成張裂縫（圖3）。

當(dāng)σ1+3σ3≥0(壓應(yīng)力為正,張應(yīng)力為負) 時,判斷張破裂的格里菲斯準(zhǔn)則的二維修正表達式如下：

式中:σ1為最大主應(yīng)力，σ3為最小主應(yīng)力，σT為二維等效張應(yīng)力。張破裂系數(shù)(η)表述為:

式中:σTC為巖石的張破裂強度,通過實驗確定。當(dāng)η≥1時,發(fā)生張破裂。η值愈大,張裂縫愈發(fā)育。張裂縫的破裂方位可根據(jù)破裂面與σ1之間的夾角β來確定,表示為:

當(dāng)（σ1+σ3）＜0時,σT=-σ3。這時張破裂方向與最大主應(yīng)力σ1的方向一致。

（2）庫侖—莫爾剪破裂準(zhǔn)則

根據(jù)庫侖—莫爾剪破裂準(zhǔn)則，巖石的破裂主要是某一個面上發(fā)生的剪切破壞，這種剪切破壞與作用在該面上的正應(yīng)力σN和剪應(yīng)力τN的相對狀態(tài)有關(guān)。表達式為:

式中:σN為巖石破裂面的剪應(yīng)力;C為巖石粘聚力,MPa;φ是巖石的內(nèi)摩擦角,°,均由實驗測定。當(dāng)某一面上的正應(yīng)力σN和剪應(yīng)力τN滿足關(guān)系式(5)時，該面上達到極限剪切平衡，開始發(fā)生剪破裂（圖4），形成共軛剪切面。

其方向可用破裂面的法線N與最大主應(yīng)力σ1的夾角α來表示，兩組破裂面空間方位不同，發(fā)生在垂直于σ1-σ3組成的平面內(nèi)，且最大主應(yīng)力平分兩組剪破裂的夾角。破裂面與最大主應(yīng)力的夾角α= （45°-φ/2）。式（5）可用主應(yīng)力表示為：

R為剪破裂系數(shù)，反映剪破裂發(fā)育程度。R＜1時地層破裂不會發(fā)生，越接近1，巖石破裂可能性增大，當(dāng)R＝1時，即巖石發(fā)生破裂；R＞1時，地層破裂可能性愈大，也就是說地層破裂的幾率愈大。同時，可以預(yù)測共軛剪破裂方位。

3.4構(gòu)造裂縫分布定量預(yù)測

根據(jù)巖心裂縫統(tǒng)計和力學(xué)性質(zhì)分析,馬嶺油田長7儲層中有效剪裂縫與有效張裂縫分別約占92.2%和7.8%。定義裂縫綜合評價指標(biāo)(Fy)表達式：

然后根據(jù)裂縫綜合評價指標(biāo)(Fy)對長7儲層裂縫發(fā)育情況進行綜合評價。一般而言,Fy愈大,裂縫發(fā)育幾率或發(fā)育程度愈大。

結(jié)合裂縫發(fā)育強度統(tǒng)計、裂縫發(fā)育方向和裂縫綜合評價指標(biāo)，對馬嶺油田長7裂縫發(fā)育有利區(qū)進行預(yù)測（圖5）。

根據(jù)馬嶺油田長7的裂縫預(yù)測結(jié)果，大致以研究區(qū)的NW- SE向?qū)蔷€為界，高角度裂縫發(fā)育區(qū)主要分布在NE半?yún)^(qū)塊，低角度裂縫發(fā)育區(qū)主要分布在SW半?yún)^(qū)塊；裂縫不發(fā)育區(qū)僅分布在西南角,局部地區(qū)也具有穿插分布的關(guān)系。

3　結(jié)論

（1）馬嶺油田長7致密油儲層露頭地層和巖心裂縫定量觀測表明,研究區(qū)長7段儲層主要發(fā)育NEE向和NNW向兩組裂縫，次發(fā)育NE向裂縫，裂縫產(chǎn)狀穩(wěn)定，延伸較遠，為區(qū)域性裂縫。研究區(qū)裂縫類型包括高角度張扭縫、壓扭縫，低角度張性縫和層理縫。以高角度裂縫發(fā)育為主，鉆遇率84.6%。所有類型的裂縫均有可能成為油氣運移通道。

（2）以統(tǒng)計學(xué)研究為基礎(chǔ),綜合巖石破裂法和能量法,以構(gòu)造應(yīng)力場數(shù)值模擬為橋梁,對馬嶺油田長7儲層構(gòu)造裂縫特征和分布規(guī)律進行了系統(tǒng)研究,建立了長7儲層構(gòu)造裂縫分布定量預(yù)測數(shù)學(xué)模型,合理預(yù)測了目的層裂縫相對發(fā)育區(qū)、次發(fā)育區(qū)和不發(fā)育區(qū)。陜西水利

參考文獻

[1]孫贊東,賈承造,李相方,等.非常規(guī)油氣勘探與開發(fā):上冊[M].石油工業(yè)出版社, 2011:1.150.

[2]鄒才能,陶士振,侯連華,等.非常規(guī)油氣地質(zhì)[M].地質(zhì)出版社,2011:1.310.

[3]梁狄剛,冉隆輝,戴彈申,等.四川盆地中北部侏羅系大面積非常規(guī)石油勘探潛力的再認識[J].石油學(xué)報,2011,32(1):8.17.

[4]賈承造,鄭民,張永峰.中國非常規(guī)油氣資源與勘探開發(fā)前景[J].石油勘探與開發(fā), 2012,39(2):129.136.

[5]《頁巖氣地質(zhì)與勘探開發(fā)實踐叢書》編委會.北美地區(qū)頁巖氣勘探開發(fā)新進展[M].石油工業(yè)出版社,2009.

[6]林森虎,鄒才能,袁選俊,等.美國致密油開發(fā)現(xiàn)狀及啟示[J].巖性油氣藏,2011,23(4): 25.30.

[7]盧雪梅.美國致密油成開發(fā)新熱點[N].中國石化報,2011.12.30(5).

[8]關(guān)德師,牛嘉玉,郭麗娜.中國非常規(guī)油氣地質(zhì)[M].石油工業(yè)出版社,1995.

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

[10]楊華,李士祥,劉顯陽,鄂爾多斯盆地致密油、頁巖油特征及資源潛力[J].石油學(xué)報, 2013.34.1:1.11.

[11]張泓.鄂爾多斯盆地中新生代構(gòu)造應(yīng)力場[J].華北地質(zhì)礦產(chǎn)雜志,1996,11(1):87.92.

[12]周新桂,操成杰,袁嘉音.儲層構(gòu)造裂縫定量預(yù)測與油氣滲流規(guī)律研究現(xiàn)狀和進展[J].地球科學(xué)進展,2003,18(3):398.404.

[13]周新桂,張林炎,屈雪峰.沿河灣探區(qū)低滲透儲層構(gòu)造裂縫特征及分布規(guī)律定量預(yù)測[J].石油學(xué)報,2009,30(2):195.200.

[14]周新桂,張利容.遼河歡喜嶺油田中潛山帶成藏期應(yīng)力場及裂縫分布規(guī)律研究[J].新疆石油地質(zhì),2003,24(1):50.54.

[15] Narr W, Lerche I. A method for estimating sub surface fracture density in core. AAPG Bull, 1984,68(5):637.648.

（責(zé)任編輯：暢妮）

中圖分類號：TE122.2

文獻標(biāo)識碼：A