曾洪流，朱筱敏，朱如凱，張慶石

（1.Bureau of Economic Geology, Jackson School of Geosciences, The University of Texas at Austin；2.中國石油大學(xué)（北京）；3.中國石油勘探開發(fā)研究院；4.大慶油田有限責(zé)任公司勘探開發(fā)研究院）

0 引言

儲集層成巖相是影響油氣儲集層質(zhì)量的一個重要因素，多年來國內(nèi)外對成巖相的研究主要局限于巖心和鏡下薄片研究，很少有人將地震資料直接用于描述和預(yù)測成巖相。理論上，只要某特定成巖作用對成巖相單元巖體的聲學(xué)特性有顯著影響，就可利用地震信息預(yù)測該成巖相。例如，碎屑巖沉積盆地中砂巖儲集層普遍受到碳酸鹽（方解石、白云石等）膠結(jié)作用的改造，砂巖孔隙中碳酸鹽膠結(jié)物所占比例差異會導(dǎo)致巖石地震波速度、密度的差異，如果碳酸鹽膠結(jié)物與其他膠結(jié)物的地震波速度、密度差別較大，且碳酸鹽膠結(jié)物所占比例較大，其地震波速度、密度與圍巖相比可能會有比較明顯的差異。例如，方解石（基質(zhì)為主）膠結(jié)物的縱波速度（vp）為6.53 km/s，密度（ρb）為 2.71 g/cm3，波阻抗（Z）為 1.77×106g/（cm2·s），其與泥質(zhì)（各種黏土成巖礦物組成的膠結(jié)物，包括基質(zhì)和微孔隙，vp=4.00 km/s，ρb=2.50 g/cm3，Z=1×106g/（cm2·s））之間波阻抗相差較大（達(dá)44%），地震檢測應(yīng)當(dāng)相對容易。利用地震資料描述和預(yù)測成巖相的主要挑戰(zhàn)在于建立巖心成巖相與可檢測的地震屬性之間的聯(lián)系，及用地震資料劃分成巖相單元并識別成巖相。

Dutta等人首次試圖將砂巖中的鈣質(zhì)膠結(jié)物含量與砂巖聲波波阻抗聯(lián)系起來，并討論了膠結(jié)物礦物成分、膠結(jié)物含量、泥質(zhì)含量、粒度、膠結(jié)層位、層序地層界面等對波阻抗的影響，但作者未能將井和模型分析的結(jié)果與實際地震資料的解釋結(jié)果聯(lián)系起來并付諸實踐[1]。筆者基于松遼盆地齊家凹陷白堊系青山口組陸相砂巖大量取心及實驗室?guī)r礦分析資料，以及相對高質(zhì)量的三維地震數(shù)據(jù)，嘗試通過巖心分析確立其成巖相和成巖序列，用測井曲線建立主要成巖相與地震參數(shù)之間的聯(lián)系，劃分其地震沉積學(xué)和成巖相單元，最終利用地震資料預(yù)測成巖相。

受地震資料分辨率的限制，常規(guī)地震地層學(xué)作圖手段（如地震層序劃分和剖面地震相分析）無法達(dá)到地震成巖相分析要求的精度，目前唯一能從常規(guī)疊后地震資料提取相對高分辨率（單個同相軸成像，10 m以下薄層檢測效果）地震地質(zhì)圖像的解釋工具是地震沉積學(xué)。地震巖性學(xué)和地震地貌學(xué)是地震沉積學(xué)的兩個組成部分，用地震巖性學(xué)方法將三維地震數(shù)據(jù)體轉(zhuǎn)化為測井巖性數(shù)據(jù)體，在該巖性數(shù)據(jù)體中，使測井巖性曲線（如自然伽馬和自然電位）在一個小的允許誤差范圍內(nèi)盡量符合井旁地震道，以保證測井資料與地震數(shù)據(jù)在儲集層尺度的最佳擬合；運用地震地貌學(xué)可以以地層切片為載體將地震數(shù)據(jù)進(jìn)一步轉(zhuǎn)換為含有巖性標(biāo)記的沉積相和成巖相平面圖。因此，地震巖性體的建立和地層切片的制作是目前砂巖成巖相地震成像的基本條件[2-5]。

1 高精度層序格架與成巖相單元關(guān)系

松遼盆地白堊系泉頭組—嫩江組發(fā)育典型大型坳陷型陸相湖盆沉積[6-7]，齊家凹陷位于松遼盆地西北部（見圖1），西、北面與盆地斜坡連接，東部與大慶長垣毗鄰，本次研究目的層段青山口組厚400～500 m，地層層序的劃分沿用前期地震沉積學(xué)的成果[8]。遵循層序地層學(xué)原則，主要依據(jù)沉積層序中的最大洪泛面確立鉆井-地震地層格架。在三維地震體中追蹤了青山口組4個比較連續(xù)的地震層位（見圖2），這些層位大多對應(yīng)于最大洪泛面，具有等時地層意義。根據(jù)測井曲線層序地層分析和地震相特征分析，進(jìn)一步在青山口組劃分了3個三級層序和6個四級層序。依據(jù)自然伽馬、自然電位曲線旋回特征和地震同相軸趨勢可進(jìn)一步追蹤五級層序，但地震剖面上的五級層序邊界解釋受地震資料分辨率限制，不能保證等時，僅有象征意義，地層切片反映五級層序最具地質(zhì)意義。

圖1 松遼盆地齊家地區(qū)位置及三維地震測網(wǎng)分布圖[8]

圖2 松遼盆地齊家地區(qū)連井地震剖面A—A?的井-震對比和青山口組層序劃分[8]

青山口組三維地震資料的主頻約為50 Hz，對應(yīng)的調(diào)諧厚度（地震分辨率極限）為10 ms，對于縱波速度為4 000 m/s的砂巖，對應(yīng)厚度為20 m，這個厚度限定了本區(qū)地震地層研究的最小作圖單元。由于五級層序（中、短期基準(zhǔn)面旋回）的平均厚度約為10～50 m，與最小作圖單元厚度大致相當(dāng)，因此，地震同相軸與五級層序（G11—G42）對應(yīng)關(guān)系良好。每個五級層序平均包括兩個地震同相軸，圖 2下部紅色同相軸大致對應(yīng)基準(zhǔn)面上升半旋回，上部黑色同相軸大致對應(yīng)基準(zhǔn)面下降半旋回。

Taylor等人觀察野外露頭發(fā)現(xiàn)，砂巖中碳酸鹽膠結(jié)物的分布與沉積環(huán)境和地層結(jié)構(gòu)有關(guān)[9]。松遼盆地白堊系各種砂巖成巖相的母體均為湖相砂巖，由于湖相砂體的分布受湖盆沉積體系的控制，影響砂巖成巖作用的一些因素如砂質(zhì)顆粒成分、大小、分選及泥質(zhì)含量、孔隙水礦化度等，均與沉積環(huán)境有關(guān)。因此，成巖相單元與沉積相單元之間具有內(nèi)在聯(lián)系。為研究方便，本文假設(shè)主要成巖相單元和沉積相單元大致相同，則成巖相單元的地震劃分與地震沉積相劃分一致。由此推論，在利用地震資料研究成巖相時，20 m不僅是地層研究的最小作圖單元，也是成巖相作圖的最小單元。當(dāng)然，成巖相的實際厚度與檢測率極限有關(guān)，筆者甚至可在一定條件下檢測到厚度為3～4 m的成巖相。

2 鈣質(zhì)膠結(jié)和泥質(zhì)膠結(jié)砂巖成巖相

根據(jù)大量巖心、鑄體薄片、普通巖石薄片及掃描電鏡資料分析，齊家凹陷青山口組儲集層主要為方解石膠結(jié)或泥質(zhì)膠結(jié)細(xì)長石細(xì)砂巖（見圖 3、圖 4），少量為介形蟲屑灰?guī)r。砂巖中長石含量 36%～46%，石英含量36%～53%，巖屑含量12%～24%，巖屑類型主要為火山巖屑。巖石顆粒粒徑約0.02～0.16 mm，分選中等—好，磨圓以次圓狀為主；表現(xiàn)為較高成分成熟度、中等結(jié)構(gòu)成熟度的巖石學(xué)特征。部分巖石組分中含較多的介形蟲碎屑和碳酸鹽膠結(jié)物。

由于膠結(jié)和壓實作用，青山口組砂巖固結(jié)致密，顆粒呈點—線接觸（見圖 3a、圖 4a），見石英次生加大，孔隙為原生粒間孔、粒間-粒內(nèi)孔的混合類型，碳酸鹽膠結(jié)物為含鐵方解石和白云石（見圖3a），黏土礦物以伊蒙混層為主，并發(fā)育伊利石和綠泥石膠結(jié)物（見圖4b），以上特征表明青山口組處于中成巖階段早期。砂巖儲集層經(jīng)歷了壓實、方解石膠結(jié)、長石和少量方解石膠結(jié)物溶蝕等一系列成巖過程，構(gòu)成的成巖序列為：早期壓實—方解石膠結(jié)—長石溶蝕—方解石交代石英顆粒，偶見壓溶作用等。

圖3 古933井2 145.8 m方解石膠結(jié)細(xì)粒巖屑長石砂巖照片

在各種成巖作用中，對儲集層質(zhì)量影響最大的成巖作用是方解石膠結(jié)，其發(fā)育程度（含量）直接影響孔隙度大小和滲透性。高方解石膠結(jié)（絕對含量大于10%）砂巖（見圖3）的粒間孔大多已被堵死，孔隙度很低，掃描電鏡下幾乎看不到任何粒間孔（見圖3b）；對比之下，低方解石膠結(jié)（絕對含量小于 10%）砂巖（見圖4）的粒間孔保留較好，孔隙度較高，掃描電鏡下粒間孔部分保留，有少量針狀伊利石和葉片狀綠泥石充填（見圖4b）。從成巖相角度出發(fā)，高方解石膠結(jié)砂巖屬鈣質(zhì)砂巖，含大量方解石膠結(jié)物，因此可將“鈣質(zhì)砂巖”命名為方解石膠結(jié)砂巖相；類似地，高孔隙、低方解石膠結(jié)砂巖一般有適量泥質(zhì)礦物存在，故將低方解石含量砂巖命名為泥質(zhì)膠結(jié)砂巖相。以下主要討論這兩種成巖相的地震定性判別問題。

3 方解石膠結(jié)砂巖相和泥質(zhì)膠結(jié)砂巖相的地震判別依據(jù)

圖4 金56井1 658.5 m泥質(zhì)膠結(jié)細(xì)粒巖屑長石砂巖照片

成巖相的地震判別依據(jù)來源于對連續(xù)取心井段巖心樣品的鏡下分析和實驗室測定資料的系統(tǒng)解剖。如在齊家凹陷金44井中，砂巖巖心薄片觀察測定的鈣質(zhì)膠結(jié)物（方解石）含量與巖心實測孔隙度（有效孔隙度）存在線性關(guān)系（見圖5）：當(dāng)孔隙膠結(jié)物中方解石相對含量從0增加到40%左右時，孔隙度從16%左右降低到小于4%。因此孔隙膠結(jié)物中方解石的比例直接影響砂巖儲集層質(zhì)量，方解石含量越高，儲集層質(zhì)量越差。孔隙膠結(jié)物中方解石相對含量與泥質(zhì)相對含量也有一定關(guān)系（見圖 5、圖 6）。巖心薄片統(tǒng)計的泥質(zhì)膠結(jié)物（伊利石、伊蒙混層、綠泥石等）相對含量為2%～20%，其對應(yīng)的方解石相對含量為0～45%。整體而言，砂巖膠結(jié)物中泥質(zhì)含量適度增加將導(dǎo)致方解石含量降低（見圖6），或儲集層質(zhì)量變好（孔隙度增高，見圖 5）。換言之，由于沉積水動力作用，砂粒之間可有少量泥質(zhì)存在的相帶（如河道充填砂體）后期方解石膠結(jié)弱，孔隙保存較好；泥質(zhì)含量較少的沉積相帶，如沿岸砂壩、三角洲前緣席狀砂等，可導(dǎo)致后期方解石強(qiáng)膠結(jié)，對砂巖成巖后孔隙保存不利。

圖5 巖心薄片測定的鈣質(zhì)膠結(jié)物（方解石）含量與巖心實測孔隙度之間關(guān)系

圖6 薄片測定的孔隙膠結(jié)物中方解石相對含量與泥質(zhì)相對含量的關(guān)系

進(jìn)一步分析表明，青山口組砂巖方解石含量與其波阻抗和密度均呈線性正相關(guān)關(guān)系（見圖7，由于用聲波速度曲線代替實驗室測定巖樣速度，增加了誤差，圖上顯示的相關(guān)度可能與實際情況有偏差）。根據(jù)這個關(guān)系，10%的方解石含量變化將導(dǎo)致約0.02～0.03的反射系數(shù)變化，足以產(chǎn)生可檢測的地震信號。如果再加入巖性的影響，就可討論用波阻抗預(yù)測成巖相的可能性。圖 8為利用測井曲線計算的不同巖性波阻抗、有效孔隙度和泥質(zhì)含量之間的關(guān)系，圖中 4種主要巖性的物性分區(qū)比較清楚：泥巖和其他巖類的區(qū)分主要考慮泥質(zhì)含量；泥質(zhì)膠結(jié)砂巖與方解石膠結(jié)砂巖、灰?guī)r的區(qū)別表現(xiàn)在有效孔隙度；方解石膠結(jié)砂巖與灰?guī)r的差別是方解石含量。半深湖泥巖與方解石膠結(jié)砂巖用波阻抗很容易判別；泥質(zhì)膠結(jié)砂巖波阻抗整體比方解石膠結(jié)砂巖低，通常情況下應(yīng)當(dāng)可以分辨；但濱淺湖泥巖與泥質(zhì)膠結(jié)砂巖波阻抗相近，容易造成混淆，在解釋中需要綜合考慮其他因素（如地層切片上解釋的沉積特征的差別）。

圖7 青山口組砂巖方解石含量與密度和波阻抗關(guān)系

圖8 4種主要巖性波阻抗、有效孔隙度和泥質(zhì)相對含量分布特征

4 90°相位地震振幅的相對波阻抗含義

利用地震資料進(jìn)行高分辨率地震成巖相分析，最好首先將一個三維地震數(shù)據(jù)體轉(zhuǎn)化為一個地層巖性數(shù)據(jù)體。在這種地層巖性數(shù)據(jù)體上進(jìn)行地層切片，可以將物理意義上的地震屬性轉(zhuǎn)換為含有巖性標(biāo)記的高分辨率沉積相和成巖相平面圖。在目前技術(shù)條件下，實現(xiàn)常規(guī)地震資料巖性標(biāo)定的最經(jīng)濟(jì)有效的方法是地震道 90°相位化[10-11]。90°相位子波將地震響應(yīng)的主波瓣最大振幅點移至薄層中間點，此時的地震響應(yīng)對應(yīng)薄層中點，而不是對應(yīng)薄層的頂或底面，這使主要地震同相軸對應(yīng)于地質(zhì)上定義的儲集層單元，如砂巖層。這樣在0～1個波長范圍內(nèi)，地震極性可與巖性相對應(yīng)。雖然當(dāng)?shù)貙雍穸刃∮谒姆种徊ㄩL時并不非常準(zhǔn)確，但地層的頂?shù)酌婵梢员淮_定在振幅過零點上。當(dāng)將這些改進(jìn)應(yīng)用于實際資料時，地震同相軸和薄地層巖性單元之間將建立一一對應(yīng)關(guān)系，將使沉積巖性的地震解釋工作，如區(qū)分砂巖和泥巖，變得更容易[4]，而這些解釋上的優(yōu)點在零相位資料和其他相位資料中并不存在。另外，地震資料的巖性轉(zhuǎn)換也可通過地震反演技術(shù)、地震屬性分析和時頻分析技術(shù)實現(xiàn)。

當(dāng)?shù)貙又械膸r性類型比較復(fù)雜，尤其還要考慮成巖變化時，對地震屬性的巖性標(biāo)定將是一個復(fù)雜過程。從金44井自然電位和聲波速度曲線分析情況看（見圖9），青山口組的巖性-速度關(guān)系相當(dāng)復(fù)雜，如泥巖包括青一段低速半深湖泥巖及青二、三段的高速淺湖—濱淺湖泥巖，速度差別較大，砂巖包括相對低速的泥質(zhì)膠結(jié)砂巖及相對高速的方解石膠結(jié)砂巖，方解石膠結(jié)砂巖速度低于但接近灰?guī)r，而且，巖性-速度（波阻抗）分布還與深度有一定關(guān)系，并明顯受到層序發(fā)育的影響。在這種情況下，用剖面上90°相位地震振幅直接判斷巖性和成巖相會有較高的多解性，但至少能看到90°相位地震數(shù)據(jù)體與相對地震波阻抗的密切聯(lián)系。如在金44井連井地震剖面上（見圖10），雖然自然電位和自然伽馬曲線顯示砂、泥巖在剖面上的分布與地震同相軸關(guān)系不是很密切，但按同相軸平均的測井速度曲線卻與地震同相軸極性有很好的對比關(guān)系，表現(xiàn)為相對高波阻抗同相軸（紅色正極性）指示測井高速層，而相對低波阻抗同相軸（黑色負(fù)極性）指示測井低速層，這種對應(yīng)關(guān)系為用地震相對波阻抗檢測方解石膠結(jié)砂巖相和泥質(zhì)膠結(jié)砂巖相提供了基礎(chǔ)。

為減少成巖相地震預(yù)測的多解性，應(yīng)當(dāng)采用地層切片，而不是地震剖面作為解釋平臺[12]。平均速度與地震同相軸極性的吻合（見圖 10）說明地震反射極性和平均振幅主要反映沉積層序（湖水進(jìn)退或體系域遷移）的旋回性，而不是沉積體系內(nèi)部的巖性巖相變化。使用地層切片可以很方便地消除這種系統(tǒng)誤差，即忽略振幅在剖面上的絕對變化，只觀察其在平面（切片）上的相對變化。此外，使用地層切片可以充分利用振幅形態(tài)顯示的地震地貌學(xué)特征，從而在振幅（波阻抗）相近的情況下，根據(jù)相帶的差別解釋巖性和成巖相，以降低解釋的不確定性。

5 青山口組地震成巖相分析

成巖相的地震判別應(yīng)當(dāng)在地震沉積相解釋的基礎(chǔ)上進(jìn)行，因為：①假定主要成巖相單元與沉積相單元大致一致，成巖相的分布與沉積相的分布有內(nèi)在聯(lián)系，預(yù)測沉積相有助于預(yù)測成巖相；②某些成巖相容易與波阻抗相近的其他成巖相或巖相混淆，在沉積相背景下判斷成巖相應(yīng)當(dāng)能有效地減小誤差。

圖9 金44井主要巖性和砂巖成巖相的測井曲線解釋

圖10 金44井過井90°相位地震剖面[8]

圖11 是取自青山口組高頻層序G31下部的振幅地層切片。由于采用黑白色標(biāo)，大量河道狀低（負(fù)）振幅異常體清晰可見，但其他形態(tài)的地震地貌單元不易判別。同一地層切片采用光譜色標(biāo)（見圖 12）改善了多種非河道型地震地貌的識別，包括沿河道邊緣分布的不規(guī)則窄帶狀相對低振幅帶，河道前端的葉狀（朵狀）相對低振幅體，與河道大致呈90°交角的相對低振幅帶和相對高振幅帶，以及分布在河道區(qū)前端的席狀和不規(guī)則形狀的相對低振幅帶和相對高振幅帶。根據(jù)相對振幅與相對波阻抗的對應(yīng)關(guān)系（見圖10），可以解釋與各地貌單元對應(yīng)的巖性。河道內(nèi)最可能充填砂巖；河道邊緣、葉狀體和帶狀體可能為砂巖、粉砂巖或泥巖；席狀體則可能為泥巖（低振幅、低電阻率）或介形灰?guī)r（高振幅、高電阻率）。綜合上述分析，研究區(qū)地震沉積學(xué)作圖單元可解釋為湖相淺水三角洲沉積體系（見圖13），來自東北（大慶長垣）方向物源的大量砂泥沉積物經(jīng)分流河道搬運進(jìn)入淺水湖盆，砂質(zhì)沉積物主要沉積在河道，河口壩沉積物一部分堆積形成葉狀體，一部分則經(jīng)波浪和沿岸流作用改造，形成三角洲前緣席狀砂體。三角洲發(fā)育經(jīng)歷了多個前積期，岸線逐漸向前遷移，形成的三角洲前緣席狀砂體多呈與分流河道正交的條帶狀。

圖11 代表高頻層序G31下部低位體系域的振幅地層切片（黑白色標(biāo)）

圖12 代表高頻層序G31下部低位體系域的振幅地層切片（光譜色標(biāo)）

圖13 由圖12 解釋的G31下部低位期淺水三角洲沉積相圖

成巖相地震判別的主要依據(jù)是砂巖中方解石膠結(jié)物含量與波阻抗的關(guān)系（見圖7）。泥質(zhì)膠結(jié)砂巖相（方解石絕對含量小于 10%）波阻抗較?。ㄒ话阈∮?.15×106g/（cm2·s），見圖8），在地層切片上振幅偏低；相反，方解石膠結(jié)砂巖相（方解石絕對含量大于10%）波阻抗較大（一般大于 1.15×106g/（cm2·s），見圖 8），在地層切片上振幅偏高。其他巖性的存在對兩種成巖相的地震判別也會有一定干擾，高波阻抗灰?guī)r的反射振幅偏強(qiáng)，與方解石膠結(jié)砂巖相類似，不易區(qū)分，但它們同屬非儲集層相，在應(yīng)用上沒有區(qū)分必要。濱淺湖泥巖與泥質(zhì)膠結(jié)砂巖相波阻抗類似，解釋時容易混淆，只有在沉積相帶清楚時才有可能分開。按照這些標(biāo)準(zhǔn)，可將地震沉積相圖（見圖 13）轉(zhuǎn)換成成巖相圖（見圖14）：分流河道發(fā)育泥質(zhì)膠結(jié)砂巖相；分流河道邊緣、三角洲前緣朵狀體及低幅/低阻三角洲前緣席狀砂體均發(fā)育泥質(zhì)膠結(jié)砂巖相（也可能外緣混有泥巖），而高幅/高阻三角洲前緣席狀砂體應(yīng)對應(yīng)方解石膠結(jié)砂巖相。最后，根據(jù)相帶位置可劃分席狀相，再根據(jù)振幅和相對波阻抗大小將它們細(xì)分為淺湖泥巖沉積區(qū)和介形灘灰?guī)r沉積區(qū)。現(xiàn)有井的測井曲線可用于標(biāo)定成巖相圖的解釋成果，如三角洲前緣席狀砂體中的方解石膠結(jié)砂巖相在井中對應(yīng)于高聲波速度層，而在所有鉆遇泥質(zhì)膠結(jié)砂巖相的井中，聲波速度都明顯降低，筆者總結(jié)的高頻層序G31下部各相帶的地震地貌學(xué)特征、振幅和波阻抗特征、巖性和沉積相解釋以及對應(yīng)的地震成巖相分析結(jié)論如表1所示。

圖14 G31下部低位體系域地震成巖相圖

從地震成巖相分析可以看出，齊家凹陷青山口組的主要有利油氣儲集層是分流河道泥質(zhì)膠結(jié)砂巖。河道砂體沉積時未經(jīng)過充分淘洗，顆粒之間含有一定量泥質(zhì)。這些泥質(zhì)對埋藏后成巖作用過程中阻礙方解石膠結(jié)、保存原生粒間孔隙起到了重要作用。當(dāng)然，也有一部分河道砂體最終還是被方解石膠結(jié)，波阻抗增加，可能已經(jīng)失去了河道地震振幅特征（比較圖13沉積相圖中解釋的河道與圖12在地層切片上實際識別的河道）。分流河道邊緣、三角洲前緣朵狀體及低幅/低阻三角洲前緣席狀砂體的泥質(zhì)膠結(jié)砂巖、粉砂巖也有一定勘探潛力。分流河道相泥質(zhì)膠結(jié)砂巖應(yīng)當(dāng)是今后齊家凹陷巖性地層圈閉鉆探的主要目標(biāo)。

表1 G31下部低位體系域地震沉積學(xué)和地震成巖相特征

砂巖成巖相地震檢測研究剛剛開始，還有許多問題需要探討，如：判斷波阻抗相近的不同成巖相時存在多解性，需根據(jù)不同情況尋找對成巖相更加敏感的地震屬性；對地層單元和成巖單元、沉積相和成巖相關(guān)系的了解還很膚淺；地震相和成巖相的標(biāo)定是定性的，還有一定隨意性等等，僅憑一個研究實例還很難解決這些問題。

6 結(jié)論

用地震資料劃分地層和成巖相單元、建立巖心成巖相與地震屬性的聯(lián)系以及尋找有效的地震成巖相成圖方法，是地震成巖相研究的 3個關(guān)鍵步驟。在松遼盆地齊家凹陷青山口組砂泥巖層段通過高分辨率井-震層序?qū)Ρ群偷貙忧衅闹谱?，建立了儲集層尺度?0 m，檢測率極限3～4 m）的沉積層序和成巖相單元格架。基于巖心的砂巖成巖相分析揭示了壓實、方解石膠結(jié)、長石顆粒和方解石膠結(jié)物溶蝕等一系列成巖現(xiàn)象并證實方解石膠結(jié)作用是影響儲集層質(zhì)量的最重要的成巖作用，方解石含量較低的泥質(zhì)膠結(jié)砂巖是主要儲集層。通過對孔隙度、泥質(zhì)含量、方解石含量和波阻抗等儲集層參數(shù)和巖石物性的分析，得出了方解石含量和儲集層孔隙度與波阻抗的關(guān)系，并認(rèn)定泥質(zhì)膠結(jié)砂巖相和方解石膠結(jié)砂巖相是適用于地震檢測的主要成巖相單元。通過90°相位化處理將地震數(shù)據(jù)體轉(zhuǎn)換成地震巖性體，用作成巖相的振幅（波阻抗）檢測信號。在用地震巖性體制作的地層切片上，識別了湖相淺水三角洲沉積體系。在沉積相圖上分析沉積相、波阻抗和成巖相的關(guān)系，最終作出了高分辨率層序的成巖相圖，可作為今后在齊家凹陷尋找?guī)r性地層圈閉的參考依據(jù)。

致謝：參加本課題研究的還有孫玉、王瑞、周川閩、白斌等，在此一并表示感謝。

[1]Dutta T, Mukerji T, Mavko G.How does carbonate cementation in sandstones affect seismic response?[R].Tulsa: SEG Las Vegas 2008 Annual Meeting, 2008.

[2]Zeng Hongliu, Backus M M, Barrow K T, et al.Stratal slicing: Part I:Realistic 3-D seismic model[J].Geophysics, 1998, 63(2): 502-513.

[3]Zeng Hongliu, Henry S C, Riola J P.Stratal slicing: Part Ⅱ: Real 3-D seismic data[J].Geophysics, 1998, 63(2): 514-522.

[4]Zeng Hongliu, Hentz T F.High-frequency sequence stratigraphy from seismic sedimentology: Applied to Miocene, Vermilion Block 50, Tiger Shoal area, offshore Louisiana[J].AAPG Bulletin, 2004,88(2): 153-174.

[5]曾洪流.地震沉積學(xué)在中國: 回顧和展望[J].沉積學(xué)報, 2011,29(3): 61-70.Zeng Hongliu.Seismic sedimentology in China: A review[J].Acta Sedimentologica Sinica, 2011, 29(3): 61-70.

[6]Yang W L.Daqing Oilfield, People’s Republic of China, a giant field with non-marine origin[J].AAPG Bulletin, 1985, 69: 1101-1111.

[7]Feng Z Q, Jia C Z, Xie X N, et al.Tectonostratigraphic units and stratigraphic sequences of the nonmarine Songliao Basin, northeast China[J].Basin Research, 2010, 22: 79-95.

[8]曾洪流, 朱筱敏, 朱如凱, 等.陸相坳陷型盆地地震沉積學(xué)研究規(guī)范[J].石油勘探與開發(fā), 2012, 39(3): 275-284.Zeng Hongliu, Zhu Xiaomin, Zhu Rukai, et al.Guidelines for seismic sedimentologic study in non-marine postrift basins[J].Petroleum Exploration and Development, 2012, 39(3): 275-284.

[9]Taylor K G, Gawthorpe R L, Curtis C D, et al.Carbonate cementation in a sequence-stratigraphic framework: Upper Cretaceous sandstones,Book Cliffs, Utah-Colorado[J].Journal of Sedimentary Research,2000, 70(2): 360-372.

[10]Zeng Hongliu, Backus M M.Interpretive advantages of 90°-phase wavelets: Part 1: Modeling[J].Geophysics, 2005, 70(3): 7-15.

[11]Zeng Hongliu, Backus M M.Interpretive advantages of 90°-phase wavelets: Part 2: Seismic applications[J]: Geophysics, 2005, 70(3):17-24.

[12]張憲國, 林承焰, 張濤, 等.大港灘海地區(qū)地震沉積學(xué)研究[J].石油勘探與開發(fā), 2011, 38(1): 40-46.Zhang Xianguo, Lin Chengyan, Zhang Tao, et al.Seismic sedimentologic research in shallow sea areas, Dagang[J].Petroleum Exploration and Development, 2011, 38(1): 40-46.