李冰,申春生,李林,張君博,胡治華
中海石油(中國)有限公司天津分公司渤海石油研究院,天津 300459
薄互層是指在縱向上儲層(砂巖)與非儲層(泥)相互交替出現且厚度均較小的一種沉積類型的巖體,垂向上呈現頻繁的砂泥交替,在地震上一個反射相位是一套砂泥交替出現的儲層綜合反映[1-4]。針對薄互層內部的砂體分布特征,很難進一步通過地震資料來預測,如何提高薄互層及內部砂體的精細研究,就成為油田開發迫切需要研究的問題。
自1973 年Widess 首次運用楔狀模型分析薄層厚度與地震響應的關系以來,國內外學者在薄互層內部定量預測薄層做了大量研究性工作。孫魯平等[5]通過從地震資料中提取薄層頻率屬性(如峰值頻率、中心頻率等)建立了地震峰值頻率與薄層厚度的關系。黃文鋒等[6]構建了楔狀薄互層模型,建立了薄互層地震響應與薄互層參數之間的聯系。李雪英等[7]利用深度域相移法對等厚薄互層進行正演模擬,實現了等厚薄互層的定量預測。上述薄互層描述方法一般以Widess 模型為基礎,利用地震屬性識別厚度小于1/4λ 的孤立薄層。然而,在實際油田中經常遇到一組由多個厚度都小于1/4λ的薄層和之間的隔夾層組成的薄互層組。其地震響應不僅與組內的單一薄層厚度有關,而且與薄層層數、物性、流體性質差異以及反射子波頻率和衰減系數有關。
利用上述時頻域方法不能有效地描述薄互層內部的砂體分布特征,如何提高薄互層及內部砂體的精細研究,就成為油田開發迫切需要研究的問題。受沉積作用的影響,薄互層內部砂體之間有密切的成因聯系,將沉積過程[8-13]融入到薄互層儲層中,可以分析薄互層儲層的合理性,重塑薄互層儲層內部結構,實現薄互層儲層內部單層砂體的形成機理,進一步加深薄互層儲層的認識。為此,筆者以渤海灣P19 油田為例,基于沉積環境分析完成了三角洲前緣薄互層儲層的精細刻畫,為類似薄互層的展布特征及剩余油分布提供參考和依據。
P19 油田位于渤海中南部海域渤南低凸起中段的東北端(圖1)。為斷裂背斜翼部構造,斷裂系統發育,地層傾角7°~18°,主要含油層段位于新近系館陶組,以辮狀河三角洲前緣為主,縱向上單油層多且厚度較薄,約70%層數的厚度小于3 m,表現為砂泥巖薄互層的特征。根據巖性及旋回特征,館陶組可劃出多套砂層組,單一砂層組主要由砂、泥巖不等厚疊置構成,為地震可識別的最小單元,具有明顯的砂—泥—砂結構(圖2)。

圖1 研究區構造位置Fig.1 Structural location in the study area
研究區主要為灰色、灰綠色細砂巖、中砂巖,巖芯上粒度呈反旋回特征(圖3a),為河口壩沉積物,碎屑顆粒呈次棱角狀(圖3b,c),接觸類型為孔隙式—接觸式膠結,分選中等—差,成熟度低,黏土礦物主要為伊利石、高嶺石,可見少量黃鐵礦(圖3d),具有辮狀河三角洲前緣沉積的特點[14-16]。
P19油田2016年實施開發井,分層系開發,井距100~300 m,鉆井資料較為豐富,含油面積3 km2,具有很大的調整空間。研究區豐富的地震資料和動態資料為本次研究奠定了堅實的基礎。

圖2 研究區館陶組地層綜合柱狀圖Fig.2 Comprehensive stratigraphic column of Guantao Formation in the study area
本次研究在層序地層學及標志層約束基礎上,地震資料約束和分級控制原則相結合,對薄互層儲層分為砂層組和單層砂體兩個級次進行精細分析。本文以L70-III砂層組為例,利用研究區地震、鉆井及測井資料,在精細時深標定基礎上進行古地形分析及鉆井揭示的沉積環境分析,再結合地震最小振幅屬性,得到砂層組的展布范圍;然后根據沉積過程分析[17-20],得到薄互層內部單層砂體級次的沉積模式及演化特征,以此確定薄互層內部單層砂體的分布范圍。從而實現了薄互層由砂層組級次到單層砂體級次的精細解剖。

圖3 研究區館陶組巖石特征Fig.3 Petrology characteristics of Guantao Formation in the study area
2.1.1 精細時深標定
聲波測井曲線是地震分析過程中必不可少的資料,隨鉆縱、橫波測井曲線的準確預測對該油田剩余油分布預測及井位優化研究工作具有很強的現實意義[21-27]。由于研究區大多數井都沒有聲波測井曲線,在沒有聲波測井資料的情況下,必須對其進行預測,目前常用的橫波測井曲線的模擬預測是基于Han 公式、Greengerg-Castagna(GC)公式和Xu-White 模型完成的。本次研究基于油田實際地震資料,實際改進Han 巖石物理公式(公式1)得到可靠時深關系,得到回歸聲波數據曲線,制作合成地震記錄(圖4a),將井揭示儲層和地震資料標定,再通過人工尋找標志層的辦法,進行進一步的精細標定,從而保證砂層組級次的高精度控制,完成砂層組級次等時地層格架的建立。

式中:DTC 為擬合縱波時差,μs/ft;TVDSS 為海拔深度,m;VCLAY 為泥質含量,%;PORT 為總孔隙度,%;PORE 為有效孔隙度,%;SW 為含水飽和度,%;a,b,c,d,e,f,g表示常數系數。
2.1.2 地震屬性分析
研究區館陶組地震主頻約為32 Hz,地震可識別砂層組的厚度約為9.7 m(波長/8)。在此基礎上,以純波地震資料為基礎,利用道積分算法將地震資料的零相位與砂泥巖阻抗界面相對應,從而實現利用振幅包絡直接刻畫目的薄互層頂底界面(圖4b),并根據該頂底界面分別提取一定時間厚度內的均方根振幅、最大振幅、最小振幅及平均振幅等地震屬性,通過與鉆井砂巖對比及分析,對比認為最小振幅屬性值與井點厚度值吻合程度較高(圖4c,d),儲層平面預測鉆井吻合率達到85%以上,呈現出砂巖越發育振幅越強的特征(圖4c,e),為研究砂層組的平面展布和連通性提供了基礎。

圖4 研究區館陶組砂層組地震剖面特征及最小振幅屬性平面圖Fig.4 Seismic profile map and minimum amplitude isoline map of Guantao Formation in the study area
2.1.3 沉積環境特征分析
古地形對沉積體系有一定的控制作用,波折帶決定了沉積的卸載場所,沉積物一般會沿著低勢地貌堆積,古地形的恢復有助于研究儲層發育與分布特點、判斷古地理環境及沉積演化特征等[28-31]。本研究首先根據實施井軌跡與地層的空間關系,消除視地層厚度的影響,結合上述等時精細地層格架,然后基于層拉平古地貌恢復技術,依據沉積環境認識,將目的層頂面拉平,視為古沉積時期的湖平面,此時底面的形態就可以認為是地層沉積前的相對古地貌。再通過編制研究區地層等厚圖,對目的層的古地貌形態加以驗證,這樣就基本恢復了古地形形態。經計算,主要沉積區域坡度一般在0.5°~2°,呈現緩坡沉積特點(圖5a),儲層展布規律受控于分流河道,沿近東西向分流,主河道寬度350~500 m。主要包括河谷、高地、斜坡三種古地貌單元。
通過地震反射特征分析,結合砂體疊置樣式,得出不同地震反射條件下的砂體疊置特點,結合巖芯觀察與測井信息建立辮狀河三角洲前緣下的主要微相單元模式(圖5b),確定了砂層組級次下的沉積微相分布圖,明確了不同成因砂體接觸關系,展現了砂體的平面連續性,最終得到砂層組平面上的砂體沉積、遷移方向、垂向砂層組內砂體發育結構等模式認知(圖5c)。河谷:為地形較低的古地貌單元,有一定的坡度,主要的水下分流河道卸載區,形成以河口壩沉積為主的微相單元,河谷內的砂體疊合厚度大于10 m,薄互層砂體多以“厚—薄—厚”式疊置;斜坡:為高地與河谷的連接區,主要發育水下分流河道邊緣沉積,當河流作用加劇時,可在斜坡形成連續砂體,砂體疊合厚度5~10 m,砂體以多個砂泥等厚疊置為主;高地:為地形最高的地貌單元,在北部、南部各發育一個小高地,為小型分流河道或溢岸的卸載區域,當水量增多時,可形成多條小型分支河道,砂體疊合厚度0~5 m,砂體以多個薄層疊置或孤立薄層為主,為薄層發育區;平原:為地形最低的古地貌單元,坡度較小,主要的水體發育區,以泥質沉積為主,分布在東側一帶。

圖5 研究區古地形及沉積環境分析Fig.5 Ancient landform and sedimentary environment in the study area
受古環境的影響,河道砂體間垂向切割、側向疊置頻繁,深化沉積模式認識對于研究砂體的組合、疊置模式有著重要意義[32-35]。研究區薄互層的形成是由于辮狀河順著水流方向的間歇性沉積引起的,當湖水較發育時,主體區以泥質沉積為主,局部會有條帶狀薄層砂體沉積,且存在擺動,當河流作用增大時,大量沉積物沿著河谷方向卸載沉積下來,形成一套優質砂體。
通過對井剖面的精細對比,反應了沉積環境水體向下變深,砂巖厚度增大、泥巖減薄、砂泥比值加大的特征,它代表了可容納空間減小的過程,辮狀河三角洲沉積中心逐漸向湖中心推進。隨著沉積時期推移,單層砂體厚度規模和展布范圍逐漸擴大,連續性變好,呈現向湖盆增長的趨勢(圖6),可在單一方向上確定單層砂體的演化關系,共經歷三個階段:1)受古地形的控制影響,氣候條件干旱,水流變小,湖泊推進,周邊物源供給減少,碎屑沉積受湖水波浪作用較大,主要形成一些小型河道沉積;2)匯入新的主干河道,水下分流河道規模有所增大,且向北遷移;3)大量辮狀河砂體順著緩坡河谷區匯入研究區,在河谷沉積形成大量砂體。這些砂體縱向上相互疊置,平面上疊加連片,向東南部逐漸尖滅,形成三角洲前緣沉積,其分布范圍取決于溝谷的寬度。

圖6 研究區三角洲前緣沉積模式Fig.6 Depositional model of braided delta front in the study area

圖7 研究區館陶組薄互層砂體疊置樣式Fig.7 Superposition pattern of thin interbed layers of Guantao Formation in the study area
辮狀河三角洲前緣薄互層內部單層砂體以多種疊置樣式存在[36-40],對于砂體的連通關系及注水效果具有很重要的影響,對于這種薄互層縱向、平面沉積變化,其精細研究可以將該砂層組的沉積模式與成因相結合,更能體現其本質。
2.3.1 縱向期次劃分
不同階段形成的砂體其規模與薄互層的沉積過程有著較好的耦合性:初期由于水下分流河道較多,沉積物較分散,隨著演化沉積物逐漸集中,多以側向沉積為主,至末期河道流量達到高峰,水下分流河道水動力較強,多以垂向加積為主,薄互層砂體呈現多種疊置樣式。
本文通過分析認為在研究區主要有三種砂體疊置關系,同位疊覆型(圖7a):主要分布在研究區高部位,其水動力較強,大量砂體優先卸載,單層砂體持續性堆積,在同一位置多期疊置;側向切疊型(圖7b):主要在河谷的邊部,受古地形與水流動力的影響,沉積不穩定,容易在側向上多期疊置,夾層數量增多;側翼接觸型(圖7c):處于研究區的邊部,為沉積遠端,為小型水下分流河道發育區,不同的河道砂體相對獨立,當水流量增多時,部分河道會相遇,導致砂體有所接觸。
對于疊置復合河道,在形成過程中新河道對老河道存在不同程度的侵蝕沖刷,因此沖刷面即為兩期河道的沉積界面。本文通過巖芯觀察發現界面回返幅度值在0.8 之下的界面在測井曲線上響應是明顯的,可以作為未取芯井中河道間界面識別和期次劃分依據。因此,在砂層組精細地層格架下,利用測井資料將薄互層細分到單層砂體級別。劃分結果:目的砂層組內部主要發育三個單層砂體,砂體厚度2~8 m,泥巖隔層厚度1~4 m(圖8a)。
2.3.2 平面砂體演化
受沉積環境的影響,相似的沉積單元具有相似的砂體分布信息。類比現代河流及周邊相似油田,辮狀河三角洲前緣水下分流河道砂體平面上多呈條帶狀,剖面具有中間厚、兩側均勻減薄或尖滅的特點[41-45],統計厚度2~6 m 的單層砂體的寬厚比為80~100。在上述砂層組沉積展布范圍基礎上,以沉積過程為指導,以相似砂體分布參數作約束,結合井點實鉆砂體厚度,最終實現單層砂體的邊界刻畫。

圖8 研究區單層砂體疊置及平面分布樣式Fig.8 Superposition and plane distribution pattern of thin interbed layers in the study area

圖9 研究區應用實例分析Fig.9 Application case analysis in the study area
不同階段形成的砂體其規模與薄互層的沉積過程有著較好的耦合性,初期小型水下分流河道眾多,以多條沉積帶為主,后期沉積物卸載集中,形成主要沉積帶,沿其方向砂體呈連續沉積,導致平面上砂體的展布空間有所不同。平面上研究區主要有兩種砂體接觸關系,軸向疊置:水動力較強,單層砂體持續性向前沉積形成,不會發生側向遷移擺動;側向遷移:單層砂體在推進過程中,填洼補平,水道多有側向遷移,然后鄰近地區再沉積形成的(圖8b)。
從圖8c 中可以看出,砂層組內部砂體從下到上逐漸變寬變厚,展現了一個湖水縮進的演化過程。通過統計,研究區館陶組單層砂體呈北西—南東向窄條狀展布,寬400~600 m,寬厚比70~150 m,北西—南東向具有一定的連通性,河道邊緣為薄層發育區。
對研究區館陶組進行了薄互層儲層的發育狀況及內部單層砂體的精細分析,取得了良好的應用效果。在實際生產中驗證了薄互層劃分結果的合理性及單層砂體的平面連續性,以J39 井組為例,B42ST02、J39井完鉆后砂體出現不同程度的虧壓,且J39、J44井完鉆后出現分段水淹(圖9a),即證實了砂體連續性,也為剩余油分布規律及下一步優化調整提供了依據。在綜合調整項目中指導井位部署和優化40 井次,以J38 井為例,根據薄互層精細刻畫結果,原設計井位位于砂體邊緣,向西南砂體變厚,據此將該井位置向西南優化了60 m,該井實施后,生產油層厚度比鉆前設計增加15%,產能增加48 m3/d(圖9b)。滾動評價新增儲量500×104m3,產能超設計1.4倍。通過對單層砂體級次的研究,提高注采對應率23%,保證了開發井的高產與注水井的顯著效果。
(1)針對地震資料難以進一步刻畫薄互層儲層的難點,本文在沉積過程分析基礎上,對三角洲前緣薄互層儲層進行了分級精細描述,為類似油田的研究和開發提供了指導。
(2)根據實際改進Han 巖石物理公式取得了較好的時深標定結果,結合古地形、地震屬性及測井分析,恢復了砂層組級次的沉積環境特征。
(3)闡明了研究區L70-III 砂層組的形成過程,共經歷了三個演化階段,不同階段單層砂體的規模及疊置樣式與其沉積過程具有較好的相關性。