張憲國, 張育衡, 張 濤, 林承焰, 王宏偉

(1.中國石油大學(華東)深層油氣重點實驗室,山東青島 266580; 2.山東省油藏地質重點實驗室,山東青島 266580;3.中國石化河南油田勘探開發(fā)研究院，河南南陽 473132; 4.山東科技大學地球科學與工程學院,山東青島266590;5.中國石油大港油田采油六廠,河北滄州 061100)

辮狀河油氣儲層在新生代陸相盆地中廣泛發(fā)育,心灘作為其最重要的沉積單元成為油氣儲層研究的熱點[1-5]。針對心灘內部結構,國內外學者開展了大量研究,涵蓋露頭、現(xiàn)代沉積和地下儲層[6-11]。但主要是對心灘及其內部結構的靜態(tài)刻畫,對其形成過程的研究甚少[8-12],尤其是對地下心灘沉積的形成過程,缺乏深入認識。因此選取黃驊坳陷典型辮狀河沉積實例,分析心灘儲層構型并利用沉積數(shù)值模擬再現(xiàn)心灘的形成過程,認識心灘形成機制,指導地下儲層構型刻畫。

1 研究區(qū)地質概況

研究區(qū)位于黃驊坳陷羊二莊油田(圖1),是黃驊坳陷趙家堡斷層下降盤的一個穹隆背斜構造?；字献韵露习l(fā)育古近系孔店組、沙河街組和東營組、新近系館陶組和明化鎮(zhèn)組以及第四系平原組(圖2)。其中館陶組和明化鎮(zhèn)組是該油田的主要含油層系[13]。研究區(qū)館陶組發(fā)育辮狀河沉積,心灘是最主要的含油砂體類型[14]。經(jīng)過40余年的開發(fā),羊二莊油田進入高含水期。心灘砂體在開發(fā)中表現(xiàn)出較強的沉積非均質性,開發(fā)層內矛盾突出,故認識心灘內部結構成為預測剩余油分布和提高采收率的關鍵。

圖1 研究區(qū)位置圖Fig.1 Location of study area

圖2 黃驊坳陷地層綜合柱狀圖Fig.2 Comprehensive stratigraphic column of Huanghua Depression

2 解剖區(qū)心灘構型特征

辮狀河河道遷移頻繁,具有砂體連片發(fā)育的特點,砂體以心灘沉積為主,呈“疊覆泛砂體”的形態(tài)展布[15]。按照Miall的河流相儲層構型劃分方法[2,16],單一心灘沉積為四級構型單元,而心灘內部的加積單元為三級構型單元。

2.1 單一心灘識別與特征

地下儲層辮狀河沉積多是多期次疊加改造的殘留沉積,其心灘沉積多以心灘復合體的形式存在,識別單一心灘是刻畫心灘內部結構的前提和基礎。研究區(qū)為開發(fā)后期密井網(wǎng)區(qū),利用井點靜態(tài)和開發(fā)動態(tài)資料可以有效識別單一心灘。

利用同一心灘砂體在測井相上的相似性可以初步對其進行識別。此外心灘底部的相對高程差也是單一心灘識別的一個依據(jù)。通常情況下,在一定范圍內河床高程相對穩(wěn)定,因此單一心灘的底部與標準層的相對距離具有穩(wěn)定性,可以作為識別單一心灘的依據(jù)。單一心灘內部同期的落淤層沉積厚度橫向變化符合自然規(guī)律,這也是識別單一心灘的一個輔助依據(jù)。

根據(jù)上述識別方法和原則,綜合考慮心灘發(fā)育規(guī)模和平面形態(tài),對解剖區(qū)的單一心灘進行識別。以圖3所示的NgⅢ1-3單砂層典型心灘為例,在代表區(qū)域泛濫沉積的穩(wěn)定泥巖標準層拉平之后,莊5-15-6井心灘底面深度為1 885.29 m,與鄰近的莊5-15-2、莊5-15-3、莊5-15-5、莊5-15-10、莊5-15-3K以及莊5-15井同期次心灘砂體底面深度相近。而在相鄰的莊G3井,同一單砂層心灘底界面深度為1 912.34 m,與前述的幾口井深度差異大,因此莊G3井與前述各井鉆遇的心灘不是同一個單一心灘。同時從落淤層發(fā)育情況看,莊5-15-6井與莊G3井落淤層發(fā)育位置及指示的垂積體厚度差異顯著(圖3(a)、(c))。

圖3 解剖區(qū)心灘發(fā)育特征Fig.3 Sedimentary features of braided bars in detailed research area

研究區(qū)40余年的開發(fā)積累了大量動態(tài)資料,尤其是示蹤劑監(jiān)測數(shù)據(jù),可以作為判斷砂體連通性和單一心灘識別的輔助依據(jù)。以莊5-15-8井組為例,注水井莊5-15-8井注入示蹤劑后,在位于同一心灘的監(jiān)測井莊5-17K中有示蹤劑顯示,指示同一心灘內部砂體的良好連通性;而由于不同心灘之間四級構型界面的滲流屏障作用,位于鄰近心灘的監(jiān)測井莊5-15-2井未見到示蹤劑顯示。

2.2 心灘內部構型特征

以落淤層為界面的砂巖增生體疊置心灘模式是中國油氣儲層研究中使用最廣的心灘沉積模式,落淤層由辮狀河洪泛事件末期懸浮落淤沉積形成,具有指示事件性沉積的意義[4-5,11,15]。研究區(qū)心灘內落淤層以夾層的形式出現(xiàn),其巖性以粉砂質泥巖和泥巖為主,在測井曲線上特征顯著。以圖3中莊5-15-5井所在的心灘為例,測井曲線上可以識別出4期落淤層沉積,在自然電位(SP)和自然伽馬(GR)測井曲線上與圍巖幅度差明顯,孔隙度和滲透率均出現(xiàn)低值特征(圖4)。

在解剖區(qū)選取完整的單一心灘,在心灘的頭部、中部和尾部各選一口井建立連井剖面,在落淤層單井識別的基礎上,結合落淤層沉積特點進行不同井點的落淤層對比(圖3、5)。根據(jù)高程和單期垂向加積砂體的厚度變化趨勢,莊5-15-2井位于心灘主體部位,發(fā)育5期落淤層,而位于心灘頭部和尾部的莊5-15-5和莊5-15-6井均識別出4期落淤層沉積。位于心灘底部的一期落淤層規(guī)模較小,形成于心灘沉積早期,對應的其下部心灘垂積體僅覆蓋莊5-15-6井和莊5-15-2井;隨著心灘發(fā)育規(guī)模的增大,在順流加積作用下,心灘順流加積生長,向下游推進,在3口井上均鉆遇心灘砂體。心灘沉積過程中,在向下游順流加積的同時,心灘頭部迎水面受到水流侵蝕作用,導致心灘頭部增生砂體和落淤層難以沉積或保存,因此位于心灘頭部的莊5-15-6井頂部一期落淤層不發(fā)育。

圖4 莊5-15-5井心灘內落淤層及心灘增生體Fig.4 Slit layer and accretion sand body in braided bars of well Zhuang 5-15-5

圖5 心灘內部結構(井點平面位置見圖3)Fig.5 Reservoir architecture of braided bar (planar location of wells is in Fig.3)

3 心灘發(fā)育過程數(shù)值模擬

從研究區(qū)心灘構型來看,除了垂向加積外,在心灘形成過程中還存在側向和順流向的沉積作用。地下儲層中保留的心灘形態(tài)及其內部結構是上述作用綜合的結果。研究通過沉積過程的數(shù)值模擬,重現(xiàn)心灘形成過程,探究不同的沉積作用如何塑造心灘砂體沉積?；谀嗌乘畡恿W的沉積過程數(shù)值模擬能夠有效避免基于沉積物理模擬和現(xiàn)代沉積觀測的沉積過程研究中存在的不足[17-19],越來越多地被應用于沉積過程和沉積機制的研究[17-23]。

3.1 模型與模擬參數(shù)

根據(jù)實例解剖區(qū)辮狀河發(fā)育規(guī)模,建立3 km×20 km的模擬河段,縱向上設定3層網(wǎng)格;根據(jù)實例區(qū)構型單元尺度和研究精度要求,設定平面網(wǎng)格尺寸為50 m×50 m,確保單一心灘內有足夠數(shù)量的網(wǎng)格顯示其內部結構,以滿足對模擬結果的觀察和分析需要。根據(jù)現(xiàn)代沉積研究的觀測尺度以及氣候和降水量的周期性,同時考慮模擬運算量,將模擬時間設置為一年,模擬計算的時間步長設置為2 min,輸出觀測記錄的時間步長為6 h。依據(jù)研究區(qū)古地貌恢復結果,將模型的地形坡度設置為0.1‰。

針對數(shù)值模擬中河流流量、含沙量等水文地質參數(shù),采用參考沉積規(guī)模相似的現(xiàn)代河流相關參數(shù)的方法來確定。根據(jù)研究區(qū)辮狀河構型刻畫結果,以河道寬度、地貌(坡度)、沉積物類型和粒度、心灘規(guī)模和幾何形態(tài)等為對比依據(jù),選取與實例區(qū)具有很高相似性的黃河龍門水文站河段和潼關水文站河段,作為相似現(xiàn)代沉積參考河段。根據(jù)近年來龍門水文站和潼關水文站的水文資料,設置模擬的入口邊界流量為40 m3/s(單位網(wǎng)格),沉積物總質量濃度為9 kg/m3。

3.2 心灘沉積演化分析

沉積過程模擬結果顯示,在辮狀河沉積開始階段,沉積物在向下游搬運的過程中逐漸沉積下來,形成心灘的早期雛形。之后在垂向加積、順流加積和側向加積的作用下,心灘規(guī)模逐漸增大,同時心灘的形態(tài)也在不斷的加積和侵蝕共同作用中產生變化和遷移(圖6)。

3.2.1 順流加積與垂向加積

在心灘發(fā)育過程中,心灘的頭部受到水流的沖刷侵蝕,而灘尾處則以沉積作用為主,出現(xiàn)順流逐漸生長的特征,在宏觀上表現(xiàn)為心灘的整體順流遷移。同時心灘頂部也在垂向加積作用下逐漸加厚(圖7)。

以圖7中標示的心灘為例,從模擬的948步到1 168步,根據(jù)平面上沉積物厚度隨時間的變化可以看出,該心灘尾部順流加積作用明顯,在沉積物不斷堆積的過程中形態(tài)從分叉狀逐漸合并成一個完整的灘尾。心灘迎水面的侵蝕作用對心灘形態(tài)的塑造要弱于灘尾的加積作用,這一點在順流向的縱向剖面圖上更加顯著。在模擬時間由948步到1 168步的過程中,灘頭侵蝕長度為141 m,遷移速度為0.64 m/步,但灘尾加積范圍達到298 m,遷移速度達到1.35 m/步(圖8縱軸正值代表沉積,負值代表侵蝕;剖面位置見圖7剖面線AB;圖中紅色虛線框為觀測的心灘位置)。灘頭侵蝕的速度和規(guī)模都要弱于灘尾的加積。

圖6 辮狀河心灘發(fā)育過程的模擬平面結果Fig.6 Simulation result of braided river deposition process on plane

心灘的垂向加積在剖面上最為顯著。以初始沉積厚度為基準,在上述沉積過程中,心灘最大沉積厚度增加了0.3 m,從剖面上不同時間的砂體頂面高度可以看到(圖8),垂向加積作用在心灘形成過程中具有普遍性,是主要的砂體沉積作用類型,但在不同階段垂向加積的規(guī)模和速度存在差異。

3.2.2 側向侵蝕與加積

在心灘順流加積和垂向加積的過程中,由于心灘兩側水道的流速和水流能量差異,流速較高的一側水道對心灘沖刷,常出現(xiàn)以侵蝕為主;流速較小的一側水道對心灘的侵蝕速度慢或者發(fā)生砂體的加積,兩側水道對心灘作用的這種差異在宏觀上表現(xiàn)為心灘側向遷移。以圖7中標示的心灘為例,從水流速度模擬結果的橫剖面上可以看到(圖9),該心灘兩側水道流速分別為0.87～1.09和0.64～0.87 m/s。左側辮狀水道流速大,對照該剖面上砂體厚度隨時間的變化發(fā)現(xiàn)(圖10),心灘左側以侵蝕作用為主,砂體厚度隨時間減小,靠近心灘出現(xiàn)水道下切,而右側水道的侵蝕作用較弱,從剖面上的砂體厚度變化來看,右側砂體加積。

圖7 心灘順流加積過程中的平面形態(tài)變化Fig.7 Planar evolution of braided bar with downstream accretion

圖8 順流向縱剖面上心灘砂體頂面與河床底部距離隨時間的變化Fig.8 Variation of distance between braided bar top and channel bottom with time at downstream section

圖9 河道橫截面水流速度分布(剖面位置見圖7剖面線DC)Fig.9 Distribution of water velocity in cross profile of channel (location of the profile is marked in Fig.7 as line DC)

圖10 辮狀河橫剖面上心灘頂面與河床底部距離隨時間的變化Fig.10 Variation of distance between braided bar top and channel bottom with time at crosssection

從圖10的心灘橫剖面砂體厚度變化可以看到,在模擬時間從1 020步到1 140步的過程中,心灘左側侵蝕寬度約為26 m,右側的側向加積寬度約為6 m,侵蝕速度大于加積速度,在該剖面位置心灘橫向寬度變小,宏觀上呈現(xiàn)向右側的遷移。在側向遷移的同時,綜合圖7和圖8的砂體厚度變化可以看到,心灘也出現(xiàn)整體向下游方向的遷移。

4 解剖區(qū)心灘發(fā)育過程

通過研究心灘內部結構的刻畫結果和心灘沉積過程的數(shù)值模擬,對解剖區(qū)心灘發(fā)育過程進行探討。解剖區(qū)心灘在形成初期規(guī)模較小,心灘生長過程中經(jīng)歷多次洪水期,心灘規(guī)模逐漸擴大。從保留下來的落淤層發(fā)育情況看,有明顯地層記錄的大規(guī)模洪水沉積事件有5期。在第一個大規(guī)模洪水期過后,心灘的發(fā)育雛形已經(jīng)明確,心灘在順流加積和頭部侵蝕作用下,宏觀上表現(xiàn)為整體的不斷向下游遷移。在遷移過程中,由于頭部侵蝕速率小于尾部的加積速率,因此心灘整體在順流方向上尺度變大(圖5、11)。在沉積數(shù)值模擬中也得到了相同的規(guī)律和定量化的依據(jù)。

心灘頂部的一期落淤層在侵蝕作用下發(fā)育范圍局限,在莊5-15-6井處未能完整保留,在測井曲線上無法明顯識別(圖5);從縱剖面上看,第一期心灘規(guī)模小,這是在莊5-15-5井上沒有看到第一期增生體頂部落淤層沉積的原因。其后在順流加積作用下心灘順流生長到莊5-15-5井位置,在該井處發(fā)育之后的幾期增生砂體。從橫剖面上看,心灘單側加積作用明顯,心灘頭部和中部在右側呈現(xiàn)出以侵蝕作用為主的特點,形成心灘兩側的不對稱加積,在宏觀上表現(xiàn)為心灘橫向遷移(圖11)。

圖11 研究區(qū)心灘演化模式(井平面位置見圖3)Fig.11 Evolution model of braided bar in study area (well position is as in Fig.3)

5 結 論

(1)將動靜態(tài)資料相結合,刻畫研究區(qū)單一心灘及其內部構型,發(fā)現(xiàn)垂向加積、順流加積、側向加積和侵蝕作用共同塑造心灘的形態(tài)和內部結構。

(2)心灘的形成是一個動態(tài)的過程,灘頭的侵蝕與灘尾的順流加積是一個同時的過程,但是二者的作用強度會存在差異,侵蝕作用會造成灘頭部位落淤層發(fā)育差。心灘兩側河道的流速差異會造成心灘一側加積而另一側侵蝕,存在不對稱的側向加積。垂向加積伴隨著心灘的整個生長過程,但是在不同階段以及心灘的不同部位,垂向加積速率存在差異。

(3)辮狀河水動力特征復雜,心灘的成因機制和影響因素多樣,研究是對地下儲層心灘發(fā)育的探討,兩岸地貌差異、物源變化、氣候等因素對心灘發(fā)育的影響沒有單獨進行詳細探討,這是下一步需要關注的問題。此外基于泥沙水動力學的數(shù)值模擬多用在水利研究中,設計模擬時間尺度小于地質時間尺度,兩個不同時間尺度上模擬結果的可對比性是沉積學家所關注的重點?？紤]到季節(jié)周期對河流沉積的控制,年度內的沉積模擬對認識辮狀河沉積演化具有參考意義。