王啟明，杜曉峰，加?xùn)|輝，徐春強(qiáng)

(中海石油(中國(guó))天津分公司渤海石油研究院，天津塘沽 300452)

遼中凹陷東二下亞段湖底扇數(shù)值沉積模擬研究

王啟明，杜曉峰，加?xùn)|輝，徐春強(qiáng)

(中海石油(中國(guó))天津分公司渤海石油研究院，天津塘沽 300452)

在無(wú)井的情況下,利用圍區(qū)古環(huán)境、三維地震等資料，應(yīng)用數(shù)值沉積模擬軟件，對(duì)渤海海域遼中凹陷東二下亞段湖底扇的形成過(guò)程及時(shí)空展布進(jìn)行了研究。研究區(qū)東二下亞段可以劃分為一個(gè)三級(jí)層序和三個(gè)體系域，湖平面曲線呈現(xiàn)脈沖式波動(dòng)變化，主物源在盆地東側(cè)和北側(cè)兩個(gè)方向。通過(guò)各個(gè)體系域的構(gòu)造圖，利用軟件恢復(fù)古地貌、古構(gòu)造，得到了相對(duì)真實(shí)的沉積背景，利用這些地質(zhì)參數(shù)進(jìn)行數(shù)值沉積模擬，進(jìn)一步得到了湖底扇分布的三維形態(tài)和演化規(guī)律。結(jié)果顯示，研究區(qū)湖底扇主要發(fā)育在在湖泛體系域內(nèi)，受先期不穩(wěn)定三角洲、湖泛域較強(qiáng)的水動(dòng)力條件、沉積坡折帶的耦合而形成，具有多期、雙物源、分布范圍廣的特點(diǎn)。

遼中凹陷；湖底扇；數(shù)值沉積模擬；湖泛體系域；儲(chǔ)層描述

斷陷湖盆凹陷中心廣泛發(fā)育湖底扇沉積，是陸相巖性地層油氣藏勘探重要的領(lǐng)域[1]。相關(guān)學(xué)者通過(guò)對(duì)湖底扇的形成機(jī)制、發(fā)育特征以及沉積微相構(gòu)成等方面的研究成果，近年來(lái)陸續(xù)發(fā)現(xiàn)了多個(gè)與湖底扇相關(guān)的巖性油氣藏[2-6]。隨著三維地震資料的應(yīng)用，在渤海海域遼中凹陷北部東二下亞段發(fā)現(xiàn)了多個(gè)在地震剖面上呈現(xiàn)“蠕蟲(chóng)狀”的地震反射特征，與渤海海域其他地區(qū)地震相對(duì)比[7]，證實(shí)為湖底扇沉積。由于研究區(qū)缺少鉆井條件，不能依賴像陸區(qū)大量的鉆測(cè)井、巖心資料，湖底扇鉆前儲(chǔ)層預(yù)測(cè)必須依賴地球物理技術(shù)，而地球物理信息多解性強(qiáng)，湖底扇富砂區(qū)精細(xì)預(yù)測(cè)存在一定的風(fēng)險(xiǎn)，難以滿足巖性油氣藏的勘探要求。數(shù)值沉積模擬是通過(guò)邏輯或數(shù)學(xué)描述方法來(lái)預(yù)測(cè)地質(zhì)系統(tǒng)和地質(zhì)沉積過(guò)程的綜合研究[8-10]，特別是在少井或者無(wú)井的地質(zhì)條件下，可以用定量方法進(jìn)行沉積機(jī)制的驗(yàn)證和輔助沉積體系時(shí)空演化規(guī)律分析[11]。本文在高精度層序地層格架的基礎(chǔ)上，對(duì)湖底扇發(fā)育層序位置、部位和形成機(jī)制進(jìn)行了研究，利用地震反射特征、古地貌、物源等資料，采用沉積模擬軟件，明確了湖底扇的形成過(guò)程及演化規(guī)律，以期在有限的條件下為海上勘探提供更可靠的信息，降低勘探風(fēng)險(xiǎn)。

1 區(qū)域地質(zhì)概況

渤海海域遼中凹陷位于遼東灣地區(qū)中部，夾持在遼西低凸起和遼東凸起之間，為東斷西超的一個(gè)箕狀凹陷，沿北東方向呈現(xiàn)長(zhǎng)條狀分布。研究區(qū)位于遼中凹陷北部(圖1)，在古近系自下而上依次發(fā)育沙河街組和東營(yíng)組地層，其中，東營(yíng)組可以分為東三段、東二下亞段、東二上亞段[7]。東營(yíng)組主要發(fā)育大規(guī)模的辮狀河三角洲沉積，在東二下亞段凹陷區(qū)發(fā)現(xiàn)了多個(gè)在地震剖面上呈現(xiàn)“蠕蟲(chóng)狀”的地震反射特征，與圍區(qū)已經(jīng)鉆井證實(shí)的湖底扇地震特征較為類似。目前，研究區(qū)內(nèi)有6口鉆井，均未鉆遇該湖底扇砂體。湖底扇具有較好的成藏條件，是該區(qū)巖性油氣藏勘探的主要目標(biāo)之一。

2 地質(zhì)參數(shù)分析

2.1 層序劃分與對(duì)比

在研究層段內(nèi)進(jìn)行層序劃分與對(duì)比，建立等時(shí)的地層格架，是開(kāi)展湖底扇沉積模擬的首要任務(wù)。運(yùn)用經(jīng)典的層序地層分析方法，結(jié)合地層的接觸關(guān)系，在東二下亞段識(shí)別出一個(gè)完整的三級(jí)層序(圖2)，可以劃分為低位體系域、湖泛體系域和高位體系域，包括兩個(gè)層序界面以及初始湖泛面(ffs)和最大湖泛面(mfs)。

從層序的發(fā)育特征來(lái)看，T3m界面(28 Ma)為東二下層序的底界面，在層序界面之上，表現(xiàn)為明顯的底超和前積現(xiàn)象，界面之下表現(xiàn)為頂超現(xiàn)象，在全區(qū)較易識(shí)別與追蹤；初始湖泛面(27.5Ma)與最大湖泛面(26.5 Ma)表現(xiàn)為穩(wěn)定、一根或多根較強(qiáng)的連續(xù)反射軸，多為湖相泥巖沉積，易于識(shí)別和追蹤，為進(jìn)積與退積作用的轉(zhuǎn)換面；T3u界面(26 Ma)為東二下層序的頂界面，界面之下為頂超現(xiàn)象。從層序的構(gòu)成來(lái)看，在低位體系域內(nèi)主要發(fā)育厚層的砂泥巖互層，為大型的辮狀河三角洲沉積；湖泛體系域?yàn)楹嗄鄮r為主，局部發(fā)育湖底扇沉積，表現(xiàn)為蠕蟲(chóng)狀斷續(xù)反射特征。一般而言，大型復(fù)合湖底扇是由多期事件性沉積形成的疊置體，在不同期次的沉積記錄上有明顯的沉積間斷，這種沉積間斷表現(xiàn)為被不同尺度的湖泛面限定或分割。東二下層序復(fù)合湖底扇其頂面和底面分別為最大湖泛面和初始湖泛面限定，可以劃分為兩期(圖2)，中間被薄層泥巖所分割。

圖1 研究區(qū)構(gòu)造位置圖

圖2 東二下層序地層劃分

2.2 相對(duì)湖平面變化分析

在陸相斷陷湖盆中，相對(duì)湖平面變化受多種因素影響[12]，變化頻繁。由于研究區(qū)位于斜坡帶，在少井的條件下，可以利用海相的濱岸上超法，結(jié)合古生物孢粉資料，分析東二下亞段的湖平面相對(duì)變化趨勢(shì)。在三級(jí)層序內(nèi)識(shí)別層序界面上部的上超點(diǎn)向上遷移(LST上超點(diǎn)、TST上超點(diǎn)、早期HST上超點(diǎn))和界面下部的上超點(diǎn)向下遷移(晚期HST上超點(diǎn))，通過(guò)古生物、古地磁等方法，記錄對(duì)應(yīng)的沉積時(shí)間、橫向相對(duì)距離和縱向深度。其中，關(guān)鍵界面的內(nèi)插點(diǎn)沉積時(shí)間根據(jù)盆地內(nèi)部對(duì)應(yīng)等時(shí)點(diǎn)與關(guān)鍵界面間泥巖厚度占上下關(guān)鍵界面間泥巖總厚度的百分比來(lái)確定。

根據(jù)相對(duì)湖平面變化曲線可以看出(圖3)，在東二下層序低位時(shí)期，相對(duì)湖平面表現(xiàn)為快速下降，之后湖平面穩(wěn)定上升，到最大湖泛面附近水面達(dá)到最大，高位域時(shí)期湖平面又緩慢回落，變化幅度較小。從圍區(qū)已鉆井的巖性來(lái)看，低位域向上泥質(zhì)含量或泥巖夾層逐漸增多,至高位域則迅速變?yōu)榇筇啄鄮r,高位域泥巖夾層又逐漸減少，反映了沉積水深由快速下降到穩(wěn)定上升、最后緩慢下降的過(guò)程。

圖3 東二下層序相對(duì)湖平面變化分析

2.3 物源及發(fā)育機(jī)制分析

東二段沉積時(shí)期處于裂陷晚期，地勢(shì)趨于平緩，受盆外的長(zhǎng)軸方向物源的影響，表現(xiàn)為沉積作用控制下的湖盆遷移夷平過(guò)程，以大規(guī)模的辮狀河三角洲沉積為主；同時(shí)，東部斷層活動(dòng)性減弱，也發(fā)育了短軸方向的辮狀河三角洲沉積。因此，研究區(qū)東營(yíng)組時(shí)期具有很好的物源補(bǔ)給，大型沉積體能夠很容易地延伸至凹陷中心部位，為湖底扇的最終形成奠定了物質(zhì)基礎(chǔ)。

東營(yíng)組早期雖然斷裂作用較弱，但在凹陷邊緣的斜坡帶，廣泛發(fā)育繼承型同沉積構(gòu)造坡折帶。湖盆底形恢復(fù)表明，東側(cè)后期斷裂改造弱，沉積坡折繼承性發(fā)育，斜坡坡度突變帶特征清晰；而北部由于東三層序穩(wěn)定沉降后，東二下層序低位域?yàn)楦呓ㄔO(shè)三角洲的差異沉積，造就坡折坡度的積累，坡度在3°左右。坡折帶頂部和中上部的富砂部位，在湖平面相對(duì)上升期，很容易遭受湖浪沖擊破壞，在坡折帶的中下部形成規(guī)模較大的湖底扇。

與低位海底扇不同，斷陷湖盆湖底扇的類型更多，成因機(jī)制更為復(fù)雜[5]。斷陷湖盆除構(gòu)造運(yùn)動(dòng)、地震作用等觸發(fā)應(yīng)力外，湖平面的周期性變化對(duì)重力流的形成也十分有利。東營(yíng)組沉積時(shí)期整個(gè)湖泊水體開(kāi)闊，湖平面變化頻繁，濱岸線的遷移和水體的震蕩變化對(duì)先期沉積物改造作用非常強(qiáng)烈。從東營(yíng)組湖平面相對(duì)變化曲線看，東三段和東二下亞段在最大湖泛期湖平面大幅震蕩，東二下亞段甚至達(dá)到了60 m。正是這種變化造成了物源系統(tǒng)的動(dòng)態(tài)變化和沉積物在坡折帶不同部位的堆積。

3 模擬結(jié)果分析

用濱岸上超法，根據(jù)鉆井標(biāo)定的湖泛面來(lái)確定湖平面變化曲線。給定東部和北部為主要物源方向，根據(jù)地層厚度計(jì)算不同時(shí)間段內(nèi)物源供應(yīng)速率，通過(guò)古構(gòu)造、古地貌的恢復(fù)，進(jìn)行可容納空間的建立，利用Dionisos軟件對(duì)研究區(qū)東二下亞段湖泛體系域(27.4 Ma到26.5 Ma)進(jìn)行了沉積模擬。從研究區(qū)沉積模擬三維圖可見(jiàn)，由于盆地差異沉降和早期三角洲沉積，低位域沉積晚期在JZ17-2井區(qū)形成了較大坡折帶。27.4Ma到27Ma時(shí)，碎屑物質(zhì)從北側(cè)和東側(cè)向湖盆搬運(yùn)，遇到地形坡折處大量沉積下來(lái)，形成了第一期湖底扇沉積，地層在JZ17-2井區(qū)厚度最大，成“橢圓”狀展布，砂巖百分含量也最高(圖4中A到B)。湖平面繼續(xù)上升，形成了一套中間的薄層泥巖，在盆地中心區(qū)域厚度最大，向岸邊尖滅減薄。27.0 Ma到26.5 Ma時(shí)，湖平面持續(xù)上升，碎屑物質(zhì)同樣從北側(cè)和東側(cè)向湖盆推進(jìn)，在第一期湖底扇發(fā)育位置上再向前推進(jìn)，發(fā)育第二期湖底扇，模擬的兩個(gè)沉積中心特征更加明顯，在兩個(gè)水道前緣區(qū)域砂巖百分含量最高，呈“梨形”狀展布(圖4中C到D)。湖底扇的兩期砂體、中間泥巖層在縱向上主要呈現(xiàn)進(jìn)積的接觸關(guān)系，與地震資料觀察到的地質(zhì)現(xiàn)象和地震特征較為一致，基本上接近實(shí)際的沉積過(guò)程。由此，根據(jù)模擬結(jié)果，建立了研究區(qū)湖底扇的沉積模式。在東二下層序低位域時(shí)期,物源來(lái)自東部和北部方向，發(fā)育大型辮狀河三角洲前緣砂體；在湖泛體系域時(shí)期，隨著相對(duì)湖平面的上升，湖泛域有較強(qiáng)的水動(dòng)力條件，先期沉積的三角洲變得不穩(wěn)定，由于早期的同沉積坡折的影響，同時(shí)在斷裂活動(dòng)等突發(fā)事件的觸發(fā)下，在前三角洲或凹陷部位間歇性地發(fā)育兩期湖底扇沉積，中間被一套穩(wěn)定分布的薄層泥巖所分割。模擬結(jié)果驗(yàn)證了湖底扇在湖泛體系域內(nèi)的形成機(jī)制，并明確了其平面分布規(guī)律。

圖4 湖底扇發(fā)育過(guò)程演化圖

4 湖底扇儲(chǔ)層精細(xì)描述

上述沉積模擬結(jié)果為研究區(qū)平面沉積相分析提供了新的信息；結(jié)合宏觀地震相、地震屬性分析、地層厚度圖等信息，可以進(jìn)行湖底扇儲(chǔ)層的精細(xì)描述。從均方根振幅、瞬時(shí)頻率、瞬時(shí)相位和半能量做的屬性聚類分析圖來(lái)看，一期湖底扇受兩個(gè)物源供給方向的影響，發(fā)育兩個(gè)沉積主體，其中，東側(cè)物源沉積主體較北側(cè)規(guī)模大。屬性顯示由東北向西南逐漸由紅色區(qū)域過(guò)渡為藍(lán)色區(qū)域(圖5中A)，巖石地球物理相關(guān)性分析認(rèn)為，這反映了砂巖逐漸減少的變化趨勢(shì)，這與趨勢(shì)沉積模擬預(yù)測(cè)的砂巖百分比變化結(jié)果相一致。此外，在靠近東側(cè)和北側(cè)為綠色區(qū)域，為內(nèi)扇發(fā)育區(qū)，紅色區(qū)域?yàn)橹猩戎黧w部分，黃色和藍(lán)色區(qū)域?yàn)橹猩葌?cè)翼和遠(yuǎn)端。

從第二期湖底扇屬性聚類圖來(lái)看(圖5中C)，由于該沉積時(shí)相對(duì)湖平面進(jìn)一步擴(kuò)大，沉積中心往前推移，紅色區(qū)域?qū)τ诘膬蓚€(gè)中扇主體范圍變大，僅在邊緣地帶發(fā)育中扇遠(yuǎn)端。研究區(qū)地震預(yù)測(cè)與沉積模擬結(jié)果相結(jié)合而確定出的湖底扇內(nèi)部構(gòu)成(圖5中B和D)所示，可以看出，除了東側(cè)物源的供給，北側(cè)也為湖底扇的發(fā)育提供了物質(zhì)來(lái)源，兩期湖底扇平面分布范圍有差異，第二期湖底扇沉積主體較大。

圖5 湖底扇儲(chǔ)層精細(xì)描述

總的來(lái)看，東二下層序沉積時(shí)期湖泛體系域內(nèi)發(fā)育的湖底扇沉積，蓋層和烴源巖較好，儲(chǔ)蓋組合優(yōu)良，圈閉條件良好，具備有利的巖性油氣藏成藏條件[13-14]，兩期湖底扇分布范圍廣，最大面積達(dá)50 km2，其中扇主體富砂區(qū)面積達(dá)30 km2，勘探潛力巨大。

5 結(jié)論

(1)遼中凹陷區(qū)東二下亞段可以劃分為一個(gè)三級(jí)層序，根據(jù)地層的接觸關(guān)系，又可以識(shí)別出低位體系域、湖泛體系域、高位體系域,其中，湖底扇主要在湖泛體系域內(nèi)發(fā)育。

(2)模擬結(jié)果再現(xiàn)了湖底扇的形成過(guò)程。在東二下層序低位域時(shí)期,物源充沛，來(lái)自東部和北部方向，發(fā)育大型辮狀河三角洲前緣砂體；在湖泛體系域時(shí)期，隨著相對(duì)湖平面的上升，湖泛域有較強(qiáng)的水動(dòng)力條件，先期沉積的三角洲變得不穩(wěn)定，由于早期的同沉積坡折的影響，同時(shí)在斷裂活動(dòng)等突發(fā)事件的觸發(fā)下，在凹陷部位間歇性地發(fā)育兩期湖底扇沉積。

(3)通過(guò)模擬分析，研究區(qū)湖底扇具有兩期、雙物源、分布范圍廣的特點(diǎn)，根據(jù)模擬兩期湖底扇三維分布及形態(tài)，結(jié)合宏觀地震相和地震屬性等信息，可為湖底扇儲(chǔ)層的精細(xì)描述提供有力依據(jù)。

[1] 姜在興.沉積學(xué)[M].北京：石油工業(yè)出版社,2003：375-402.

[2] 鄧毅林,王天琦,曹正林,等.湖底扇沉積模式及地震響應(yīng)特征[J].大慶石油地質(zhì)與開(kāi)發(fā),2010,29(5):35-42.

[3] 楊有星,金振奎.黃驊坳陷港中地區(qū)古近系沙三段湖底扇沉積相特征[J].現(xiàn)代地質(zhì),2012,26(2):355-362.

[4] 王建功,王天琦,張順,等.松遼坳陷盆地水侵期湖底扇沉積特征及地球物理響應(yīng)[J].石油學(xué)報(bào),2009,30(3):361-366.

[5] 潘樹(shù)新,鄭榮才,衛(wèi)平生,等.陸相湖盆塊體搬運(yùn)體的沉積特征識(shí)別標(biāo)志與形成機(jī)制[J].巖性油氣藏，2013, 25(2):9-18.

[6] 龐軍剛,李文厚,石碩,等.鄂爾多斯盆地長(zhǎng)7段濁積巖沉積演化模式及石油地質(zhì)意義[J].巖性油氣藏, 2009, 21(4):73-77.

[7] 吳奎,吳俊剛,張中巧,等.遼中凹陷北部湖底扇沉積模式及地震響應(yīng)特征[J].東北石油大學(xué)學(xué)報(bào),2012, 36(5):33-37.

[8] 龐雄奇.地質(zhì)過(guò)程定量模擬[M].北京:石油工業(yè)出版社,2003：1-13.

[9] Burton R, Kendall C G C, Lerche I.On the impossibility of fathoming eustasy from the stratigraphic record[J].Earth Science Reviews, 1987，(24): 237-277.

[10] Cross T A. The scientific basis for stratigraphic inversion[J].Ciencia Tecnology Futuro , 2000， (1):824-845.

[11] 呂明,王穎,徐徽.沉積模擬方法在Bonaparte盆地的應(yīng)用[J].中國(guó)海上油氣,2010,22(2):83-90.

[12] 雷燕平,林暢松,劉景彥,等.海拉爾盆地貝爾凹陷下白堊統(tǒng)層序地層與沉積體系分析[J].石油地質(zhì)與工程,2007,21(5):11-15.

[13] 王寧,何幼斌,王思文,等.惠民凹陷商河油田商三區(qū)沙三上亞段濁積扇沉積特征[J].巖性油氣藏,2014,26(3):38-44.

[14] 孫曉明.遼河盆地西部凹陷近岸湖底扇油氣富集條件研究[J].石油地質(zhì)與工程,2007, 21(5):20-23.

編輯：吳官生

1673-8217(2015)04-0015-03

2014-09-01

王啟明，地質(zhì)工程師，碩士，1984年生，2010年畢業(yè)于中國(guó)地質(zhì)大學(xué)(北京)礦產(chǎn)普查與勘探專業(yè)，現(xiàn)從事石油地質(zhì)綜合研究工作。

國(guó)家十二五重大專項(xiàng)“大型油氣田及煤層氣開(kāi)發(fā)”(2011ZX05023-002)。

TE111.3

A