地震波形指示反演在東營凹陷王家崗地區(qū)沙四上亞段灘壩砂的應(yīng)用

韓長城, 林承焰,2, 任麗華,2, 馬存飛, 魏 婷, 張憲國,2, 孫志峰

(1.中國石油大學(xué)球科學(xué)與技術(shù)學(xué)院,山東青島 266580;2.山東省油藏地質(zhì)重點(diǎn)實(shí)驗(yàn)室,山東青島 266580;3.中國石油新疆油田分公司石西油田作業(yè)區(qū),新疆克拉瑪依 834000)

地震波形指示反演在東營凹陷王家崗地區(qū)沙四上亞段灘壩砂的應(yīng)用

韓長城1, 林承焰1,2, 任麗華1,2, 馬存飛1, 魏 婷3, 張憲國1,2, 孫志峰1

(1.中國石油大學(xué)球科學(xué)與技術(shù)學(xué)院,山東青島 266580;2.山東省油藏地質(zhì)重點(diǎn)實(shí)驗(yàn)室,山東青島 266580;3.中國石油新疆油田分公司石西油田作業(yè)區(qū),新疆克拉瑪依 834000)

以東營凹陷王家崗沙四上亞段灘壩砂為例,在層序地層格架建立的基礎(chǔ)上,利用沉積學(xué)原理,基于貝葉斯理論將地震、測井和地質(zhì)統(tǒng)計學(xué)信息融合為地層模型參數(shù)的后驗(yàn)概率分布,采用Metropolis-Hastings抽樣算法對后驗(yàn)概率分布隨機(jī)抽樣,獲得反演解,從而預(yù)測有利灘壩相帶。結(jié)果表明:東營凹陷王家崗地區(qū)沙四上亞段低位域發(fā)育灘壩砂,其中準(zhǔn)層序組2沉積時期以壩砂為主,砂體厚且橫向分布穩(wěn)定,而準(zhǔn)層序組1和準(zhǔn)層序組3沉積時期灘砂和壩砂間互發(fā)育,砂體厚度相對較小;湖侵域和高位域主要發(fā)育碳酸鹽巖灘壩,灘壩砂零星狀分布。該方法體現(xiàn)了相控的思想,預(yù)測結(jié)果更符合地質(zhì)規(guī)律,在東營凹陷沙四上亞段灘壩儲層預(yù)測中取得了較好的效果,為同類型儲層預(yù)測提供了有利支持。

波形指示反演; 貝葉斯理論; 沙四上亞段; 灘壩砂

近年來,隨著中國油氣勘探開發(fā)進(jìn)入中—高勘探階段,油氣勘探技術(shù)逐漸成熟,油氣勘探逐步進(jìn)入隱蔽油氣藏勘探階段,灘壩砂巖儲層逐漸引起人們的關(guān)注[1]。國內(nèi)外許多學(xué)者對灘壩的地質(zhì)特征和識別預(yù)測技術(shù)進(jìn)行了大量研究,灘壩砂具有單層砂體厚度薄、橫向變化快、地震資料縱向分辨率低等特點(diǎn)。朱筱敏等[2]、陳世悅等[3]根據(jù)灘壩的主要物質(zhì)組成將其分為砂質(zhì)灘壩和生物碎屑灘壩,楊勇強(qiáng)等[4]建立了陸相斷陷湖盆的灘壩沉積模式,劉書會等[5]利用地震屬性對灘壩薄砂體進(jìn)行了預(yù)測,才巨宏[6]、韓宏偉[7]等利用波形分析技術(shù)和地震特征反演技術(shù)對灘壩儲層進(jìn)行了預(yù)測,黃捍東等[8]通過非線性隨機(jī)反演方法對陸相薄砂巖儲層進(jìn)行了預(yù)測。這些研究成果均較好地指導(dǎo)了灘壩勘探,但仍然存在較多問題:地震屬性多解性較強(qiáng);波阻抗反演分辨率太低,無法滿足薄層需求;測井約束反演高頻分量主要來源于初始阻抗模型,多解性強(qiáng);地質(zhì)統(tǒng)計學(xué)反演隨機(jī)性較強(qiáng),對井位分布要求較高。筆者基于相控模型的高精度反演,在層序地層格架約束下,運(yùn)用沉積學(xué)基本原理,充分利用地震波形的橫向變化反映儲層空間的相變規(guī)律,實(shí)現(xiàn)灘壩儲層和沉積相帶的預(yù)測。

1 反演方法原理

1.1 目標(biāo)函數(shù)建立

在地震勘探中,由聲波計算得到的波阻抗曲線通常表示為

Y=X+N.

(1)

式中,Y為由聲波計算得到的波阻抗值;X為待求解的地下地層實(shí)際波阻抗值;N為隨機(jī)噪音,一般情況下,假設(shè)N為服從高斯分布的平衡白噪音,其數(shù)學(xué)期望為0,協(xié)方差為σ。

(2)

波阻抗反演過程就是通過公式(2)研究X為何值時使J1達(dá)到最小。由于地震資料帶限寬度的限制,只能提供波阻抗反演的中頻成分,而高頻和低頻成分則主要通過測井?dāng)?shù)據(jù)、地質(zhì)信息獲得[9]。波阻抗反演存在不穩(wěn)定性和多解性,僅使用一種方法很難獲得較理想解,因此須在目標(biāo)函數(shù)中引入先驗(yàn)信息對其進(jìn)行約束,從而獲得較穩(wěn)定的反演解。在統(tǒng)計方法中,通常采用最大后驗(yàn)估計,則目標(biāo)函數(shù)為

J(Z)=J1(Z)+λ2J2(Z).

(3)

即可表示為測量數(shù)據(jù)估算可信度項與先驗(yàn)項的和。其中Z為波阻抗,先驗(yàn)項J2(Z)可以定義為勢函數(shù)的和,是與測井資料、地質(zhì)信息等先驗(yàn)信息有關(guān)的函數(shù)。J1(Z)為與某些后驗(yàn)信息有關(guān)的函數(shù)。λ為平滑參數(shù),用于協(xié)調(diào)J1(Z)和J2(Z)之間相互影響,為常數(shù)。公式(3)進(jìn)一步改寫為

(4)

式中,φ為勢函數(shù);λ和δ為正則參數(shù),λ用于平衡數(shù)據(jù)項和先驗(yàn)項之間的相互影響,δ是一種刻度參數(shù),調(diào)節(jié)不連續(xù)處梯度值[10]。公式(4)由兩部分組成,第一部分是關(guān)于測量模型和數(shù)據(jù)之間的某種一致性,當(dāng)其值達(dá)到最小時,一致性會變得很好,但解穩(wěn)定性變差;第二部分是關(guān)于解 “光滑性”的測量標(biāo)準(zhǔn),當(dāng)其值取最小時,獲得的近似解較穩(wěn)定。

1.2MCMC方法和Metropolis-Hastings抽樣準(zhǔn)則

馬爾科夫鏈蒙特卡羅(MCMC)方法是在貝葉斯框架下,用已有資料進(jìn)行約束,既可使最優(yōu)解滿足參數(shù)的統(tǒng)計特性,同時融入了先驗(yàn)信息,解的精度得到提高[11]。其基本思想是:通過構(gòu)建一個平穩(wěn)分布為π(x)的馬爾科夫鏈來獲得π(x)的樣本,被估計的模型參數(shù)值是鏈的狀態(tài)空間,被估計的后驗(yàn)分布為鏈的貝葉斯后驗(yàn)分布。

利用MCMC方法對貝葉斯理論下的后驗(yàn)概率分布進(jìn)行隨機(jī)抽樣,反演問題希望得到的是所估參數(shù)的后驗(yàn)概率分布,因此各個參數(shù)的馬爾科夫鏈應(yīng)收斂于所估計參數(shù)的后驗(yàn)概率分布[11]。

設(shè)M為某一空間,n為產(chǎn)生的總樣本數(shù),m為馬爾科夫鏈達(dá)到平穩(wěn)時的樣本數(shù),則MCMC的思路為[12]:

(1)構(gòu)建一條Markov鏈,使其收斂至平穩(wěn)分布π(x),常用的構(gòu)造轉(zhuǎn)移方法是Gibbs抽樣和Metropolis-Hastings抽樣。

(2)產(chǎn)生樣本。由M中的某一點(diǎn)x(0)出發(fā),用(1)中的Markov鏈進(jìn)行抽樣模擬,產(chǎn)生點(diǎn)序列:x1,…,xn。

(3)蒙特卡羅積分。任一函數(shù)f(x)的期望估計為

Metropolis-Hastings抽樣的基本思路[12]是:選擇一轉(zhuǎn)移函數(shù)q(x;x(i-1))和初始值x(0),若第i次迭代開始時的參數(shù)值為x(i-1),則第i次迭代過程為

(1)從q(x;x(i-1))中抽取一個備選值x′;

(3)以概率α(x(i-1),x′),置x(i)=x′,以概率1-α(x(i-1),x′),置x(i)=x(i-1);

(4)重復(fù)(1)～(3)步n次,則可以得到后驗(yàn)樣本x(1),x(2),x(3)…x(n),從而可計算后驗(yàn)分布的各階矩以進(jìn)行統(tǒng)計推理。

2 技術(shù)流程

圖1 波形指示反演流程Fig.1 Flow chart of waveform indication inversion

波形指示反演技術(shù)流程如圖1所示。首先對測井資料處理,包括測井曲線標(biāo)準(zhǔn)化、重構(gòu)和敏感曲線分析,建立層序地層格架,沉積特征及沉積相特征分析;然后分析測井?dāng)?shù)據(jù)和地震波形特征,參照樣本空間分布距離和曲線分布特征,優(yōu)選與待判別道波形特征相似程度較高的樣本井建立初始模型,統(tǒng)計其縱波阻抗建立先驗(yàn)概率函數(shù);將初始模型與地震波阻抗體進(jìn)行匹配濾波,求得似然函數(shù),對樣本進(jìn)行多尺度分解,逐步濾除高頻成分;最后基于貝葉斯理論,聯(lián)合似然分布與先驗(yàn)分布得到后驗(yàn)概率分布,并將其作為目標(biāo)函數(shù);基于波形指示優(yōu)選的樣本在空間上具有較好的相關(guān)性,采用Metropolis-Hastings抽樣算法對后驗(yàn)概率分布抽樣,選取目標(biāo)函數(shù)最大值時的解作為可行隨機(jī)實(shí)現(xiàn),求取多次可行實(shí)現(xiàn)的均值作為期望輸出。

3 應(yīng)用實(shí)例

東營凹陷王家崗地區(qū)沙四上亞段主要發(fā)育灘壩相沉積,灘壩砂巖油藏是其中一種重要的巖性油藏類型[13],分布面積廣、認(rèn)識程度較低,但勘探潛力巨大。實(shí)際勘探開發(fā)中往往會遇到一系列困難:巖性組合為砂泥巖薄互層,單層砂體多小于2 m,儲層物性差,橫向連續(xù)性差,分布規(guī)律難以掌握;同時地震資料主頻低,目的層段約為30 Hz,有效頻帶范圍10～70 Hz,目的層速度約為2 000 m/s,地震資料對儲層的分辨能力偏低,對灘壩砂巖儲層地震相特征認(rèn)識不清,導(dǎo)致灘壩砂體識別和預(yù)測難度大[14-16]。

3.1 區(qū)域地質(zhì)背景

渤海灣盆地是中國東部最重要的含油氣盆地,濟(jì)陽坳陷位于渤海灣盆地的東南部,其中東營凹陷位于濟(jì)陽坳陷的東南部,是渤海灣盆地中一個重要的三級構(gòu)造單元。王家崗地區(qū)位于東營凹陷南斜坡的東段,其西北部為牛莊洼陷,東北部為廣利洼陷,東部為八面河斷裂鼻狀構(gòu)造帶,南部為廣饒凸起。構(gòu)造位置處于近東西向陳官莊-王家港斷裂帶與東北向八面河斷裂鼻狀構(gòu)造帶交匯處,西南部為近北西向純化草橋斷裂鼻狀構(gòu)造帶(圖2)。

3.2 層序地層格架建立

建立層序地層格架為開展儲層預(yù)測和沉積相研究提供了地震地質(zhì)的等時格架,在地震沉積相研究中不可缺少。王家崗地區(qū)沙四上亞段為一個完整的三級層序[17],應(yīng)用初次湖泛面和最大湖泛面將其劃分為低水位體系域、湖侵體系域和高水位體系域三個體系域[18]。沙四上亞段頂界面在地震剖面上對應(yīng)于T6反射層,反射特征為強(qiáng)振幅、強(qiáng)連續(xù)反射,全區(qū)內(nèi)可以較大范圍追蹤,為沙四段與沙三段的分界面;沙四上亞段層序內(nèi)部以沉積作用轉(zhuǎn)換面T7為界,為低水位體系域向湖侵體系域過渡面,地震剖面上中強(qiáng)振幅、中連續(xù)反射特征;T7′為沙四上亞段和沙四下亞段的分界面。湖泊低水位體系域?qū)?yīng)準(zhǔn)層序組1、準(zhǔn)層序組2和準(zhǔn)層序組3,準(zhǔn)層序組1湖泊水體較淺,準(zhǔn)層序組2時期水體開始逐漸變淺至準(zhǔn)層序組3時最淺,主要發(fā)育砂泥間互的灘壩相沉積;湖侵體系域?qū)?yīng)準(zhǔn)層序組4,湖泊水體開始快速上升,發(fā)育一套分布穩(wěn)定高導(dǎo)泥巖;高水位體系域?qū)?yīng)于準(zhǔn)層序5、準(zhǔn)層序組6和準(zhǔn)層序組7,早期湖泊水體深且湖平面穩(wěn)定,沉積厚層泥巖,晚期水位逐漸下降,有一定量的沉積物供應(yīng),形成一定分布范圍的灘壩相沉積(圖3)。

圖2 東營凹陷區(qū)域構(gòu)造Fig.2 Structure units in Wangjiagang area, Dongying Depression

圖3 東營凹陷南北向?qū)有虻貙悠拭鍲ig.3 North-south profile of sequence stratigraphic in Wangjiagang area, Dongying Depression

3.3 主要技術(shù)環(huán)節(jié)

3.3.1 測井曲線處理

由于測井儀器誤差等因素導(dǎo)致測井?dāng)?shù)據(jù)出現(xiàn)偏差,采用頻率直方圖法進(jìn)行測井資料的標(biāo)準(zhǔn)化處理;統(tǒng)計分析尋找對巖性變化比較敏感的測井曲線,通過交會分析認(rèn)為波阻抗曲線對砂泥巖區(qū)分較差,砂巖和泥巖部分阻抗值重疊(圖4(a)),無法用現(xiàn)有波阻抗曲線進(jìn)行反演,但自然電位曲線對砂泥巖區(qū)分較好。通過對聲波曲線重構(gòu),提高地震反演的分辨率和精度,將聲波中體現(xiàn)地層背景速度的低頻信息與巖性敏感曲線自然電位曲線的高頻信息通過信息融合技術(shù)進(jìn)行融合,形成重構(gòu)聲波曲線。通過重構(gòu)聲波曲線得到的波阻抗對砂巖和泥巖的區(qū)分較好(圖4(b))。

圖4 王家崗地區(qū)灘壩砂巖和泥巖波阻抗概率直方圖Fig.4 Histogram of impedance of beach-bar sandstone and mudstone in Wangjiagang area, Dongying Depression

3.3.2 子波提取

子波的準(zhǔn)確提取是地震反演的基礎(chǔ)和關(guān)鍵。本文中根據(jù)地震資料的統(tǒng)計特性和測井資料的特點(diǎn),采用循環(huán)迭代法子波提取。其一般提取流程為:首先使用地震道的自相關(guān)性計算初始子波的振幅譜,結(jié)合測井曲線確定初始子波的相位譜(假定為線性相位);然后利用振幅譜和相位譜的信息合成一個理論雷克子波,用理論子波合成的記錄初步進(jìn)行標(biāo)定。在初步標(biāo)定的基礎(chǔ)上,提取井旁地震道實(shí)際地震子波,用該子波制作合成記錄并重新調(diào)整時深關(guān)系。如此反復(fù),直至得到相位振幅譜變化穩(wěn)定的精細(xì)子波。

3.3.3 頻率參數(shù)設(shè)置

頻率參數(shù)影響反演的分辨率,控制反演結(jié)果的有效頻帶范圍。地震資料有效信息主要集中在中頻;測井曲線屬于寬頻帶信息,能夠?yàn)榉囱萏峁┑皖l信息和高頻信息;在地震波形指示反演中,低頻成分主要從測井資料中獲取,中頻成分通過地震資料求取,高頻成分是在地震波形指示下進(jìn)行的隨機(jī)模擬結(jié)果,整個過程是由確定到隨機(jī)的逐步過渡,隨機(jī)成分較少。本次反演的低頻參數(shù)和高頻參數(shù)優(yōu)選如下:低頻參數(shù)中高通頻率取8 Hz,高截頻率取15 Hz;高頻參數(shù)中低截頻率取70 Hz,低通頻率取110 Hz,高通頻率取300 Hz,高截頻率取350 Hz。

3.4 沉積相類型及波形特征

3.4.1 碎屑巖灘壩

根據(jù)東營凹陷沙四段上亞段巖心、錄井、測井等資料將碎屑巖灘壩劃分為灘砂和壩砂,其具有不同的沉積特征(圖5)。

灘砂亞相發(fā)育于濱淺湖地區(qū),平行岸線分布[1],砂層多而厚度薄,平面上呈席狀或帶狀廣泛分布,物性較差,沉積構(gòu)造主要有波狀層理、浪成砂紋層理、沖洗交錯層理及壓扁層理等,可見植物根化石。根據(jù)沉積砂體厚度、粒度、分布特征等進(jìn)一步分為灘席和灘脊兩種微相[19]。灘席微相砂體粒度細(xì),巖性為灰色泥質(zhì)粉砂巖、粉砂巖,砂泥巖頻繁薄互層,砂體厚度小于2 m,自然電位曲線為中低幅鋸齒狀;灘脊微相巖性以粉砂巖、泥質(zhì)粉砂巖、粉砂質(zhì)泥巖為主,自然電位曲線為中高幅鋸齒狀。

壩砂亞相分布在灘砂中,粒度相對較粗。砂體單層厚,呈條帶狀或透鏡狀,分布面積較小,物性較好,沉積構(gòu)造主要有平行層理、波狀層理、浪成砂紋層理和塊狀層理,見植物根化石、介殼化石等。壩亞相包括壩主體和壩邊緣兩種微相類型。壩主體微相是主體部分,粒度粗,巖性以厚層相對均質(zhì)的灰褐色粉砂巖、粉細(xì)砂巖為主,粒序呈反韻律或復(fù)合韻律,自然電位曲線為中高幅漏斗形、箱形組合型。壩邊緣微相巖性主要為灰色粉砂巖、泥質(zhì)粉砂巖、粉砂質(zhì)泥巖,粒序呈反韻律或復(fù)合韻律,自然電位曲線為中幅漏斗形或指形。

濱淺湖泥以灰色、深灰色泥巖、頁巖為主,砂體不發(fā)育,自然電位曲線為低幅或微幅齒形或平滑線形。

3.4.2 碳酸鹽巖灘壩

碳酸鹽灘壩多分布于靠近物源區(qū)且無大量陸源碎屑物質(zhì)供給的湖灣地區(qū)[20]。巖性主要為黃色頁狀灰?guī)r、深灰色碎屑石灰?guī)r和灰質(zhì)頁巖,灘單層厚度多小于2 m,壩單層厚度多大于2 m,感應(yīng)電導(dǎo)率、聲波時差為低值尖峰,電阻為高值,自然電位曲線值較高,低幅齒形。

3.4.3 波形特征

根據(jù)正演模擬及樣本井處沉積特征,分析不同巖性組合特征的波形特征:①壩砂亞相單砂層厚度較大,地震反射波形特征為半極值寬度大、幅度小、反射能量中等的中幅單峰波形,且壩砂反射振幅與壩砂厚度呈正相關(guān),壩砂發(fā)育處地震反射振幅明顯大于壩砂不發(fā)育處[21](圖5(a));②灘砂亞相砂泥巖薄互層,地震反射波形特征為中—強(qiáng)單峰波形和復(fù)合波波形。當(dāng)厚層泥巖夾薄層砂巖時,反射特征為半極值寬度小、幅度大、反射能量較強(qiáng)的中—強(qiáng)單峰波形(圖5(b));當(dāng)雙層砂巖夾薄層泥巖時,反射特征為復(fù)波波形,中等偏下振幅,隨上、下兩側(cè)砂巖的多少而呈不對狀,且振幅強(qiáng)處砂巖總厚度小而振幅弱處砂巖厚度反而大[7](圖5(c)、(d))。

圖5 東營凹陷沙四上亞段灘壩砂巖沉積及地震波形特征Fig.5 Sedimentary and seismic wave characteristics of beach-bar sand of E in Wangjiagang area, Dongying Depression

3.5 反演效果分析

采用約束稀疏脈沖反演和波形指示反演2種反演方法對目的層灘壩砂巖進(jìn)行預(yù)測,反演結(jié)果分別見圖6和圖7。低阻值代表泥巖,高阻值代表砂巖,約束稀疏脈沖反演和地震波形指示反演的結(jié)果整體趨勢相同,但約束稀疏脈沖反演分辨率較低,對薄層和橫向變化較快的儲層預(yù)測難度較大;而地震波形指示反演分辨率高,反演剖面中—高阻發(fā)育規(guī)律與連井剖面上砂體發(fā)育規(guī)律基本一致,反演阻抗尖滅點(diǎn)可以代表砂巖尖滅點(diǎn),砂巖尖滅自然,同時可以體現(xiàn)出砂巖的內(nèi)部結(jié)構(gòu)和沉積演化規(guī)律。地震波形指示反演具有寬頻阻抗信息,其中低頻成分通過井插值模型獲得,中頻成分來自地震資料的相對阻抗,而高頻成分是地震波形指示下的隨機(jī)模擬結(jié)果。地震波形是地下地質(zhì)體地震響應(yīng)參數(shù)的綜合反映,不同地層巖性、物性的變化都會引起反射特征的變化。該反演方法在提高反演分辨率的同時,充分利用地震波形的橫向變化來反映儲層空間的相變特征,體現(xiàn)了相控思想,是一種真正的井震結(jié)合反演方法。反演結(jié)果垂向分辨率和橫向分辨率明顯得到提高,符合地質(zhì)規(guī)律,反演效果好。該反演方法在沉積相穩(wěn)定且斷層發(fā)育較少,測井資料豐富,地震資料信噪比高,主頻越高越好,地震波形與沉積相變有明顯對應(yīng)關(guān)系的地區(qū)應(yīng)用效果較好,特別適用于成熟探區(qū)和開發(fā)區(qū)塊的薄層精細(xì)預(yù)測。

圖6 稀疏脈沖反演波阻抗剖面Fig.6 Parse spike inversion section of impedance

圖7 波形指示反演剖面Fig.7 Waveform indication inversion section of impedance

3.6 沉積體系展布規(guī)律

沙四上純下沉積時期東營湖盆處于斷陷初期,氣候由干旱向潮濕轉(zhuǎn)變,凹陷南斜坡相對平緩的古地貌背景,決定了廣闊的湖區(qū)以濱淺湖灘壩沉積為主[22]。基于層序地層格架下對反演數(shù)據(jù)體提取地層切片,從平面分布上預(yù)測砂體的展布范圍和演化規(guī)律(圖8)。低位域時期湖盆面積小,沉積水體整體較淺。準(zhǔn)層序組1沉積期剛進(jìn)入濱淺湖沉積初期,物源主要來自東南部,物源供應(yīng)不足,以灘砂沉積為主;準(zhǔn)層序組2沉積期水體稍有加深,之后水體維持相對穩(wěn)定,沉積物由東南向西北推進(jìn),在水下高地處沉積灘壩砂體,壩砂非常發(fā)育,且砂體分布穩(wěn)定,由北部向南砂體厚度逐漸減薄;準(zhǔn)層序組3沉積期基準(zhǔn)面持續(xù)上升,基準(zhǔn)面上升緩慢,沉積物以加積為主,灘砂和壩砂間互發(fā)育,灘壩砂體相對發(fā)育,砂體分布面積增大而泥巖面積縮小,西北部遠(yuǎn)離物源處的水下高地沉積碳酸鹽巖灘壩。湖侵域時期基準(zhǔn)面迅速上升并相對穩(wěn)定,只零星發(fā)育碳酸鹽巖灘壩和灘砂,以碳酸鹽巖灘壩為主。高位域時期湖盆面積較湖侵域變化不大,水體相對變淺。準(zhǔn)層序組5和準(zhǔn)層序組6沉積期,灘壩砂呈孤立橢圓狀分布,遠(yuǎn)離物源處的水下高地沉積碳酸鹽巖灘壩;準(zhǔn)層序組7沉積期基準(zhǔn)面略有下降,沉積物整體具進(jìn)積特征,水下高地處沉積灘壩砂體,西北部遠(yuǎn)離物源處沉積碳酸鹽巖灘壩,碳酸鹽巖灘壩分布面積較大(圖9)。

圖8 東營凹陷沙四上亞段灘壩反演體切片F(xiàn)ig. 8 Inversion-slice of beach bar sandstone of E in Wangjiagang area, Dongying Depression

圖9 東營凹陷沙四上亞段灘壩沉積微相Fig.9 Sedimentary microfacies of beach-bar sandstone of E in Wangjiagang area, Dongying Depression

4 結(jié) 論

(1)本文中反演方法是在層序地層格架下利用沉積學(xué)原理,基于貝葉斯理論下的反演過程,充分利用地震波形的橫向變化反映儲層空間的相變規(guī)律,更好地體現(xiàn)了相控的思想,分辨率明顯提高,預(yù)測結(jié)果更符合地質(zhì)規(guī)律,是一種新的井震結(jié)合反演方法。

(2)東營凹陷沙四上亞段在低位域時期準(zhǔn)層序組2時期以發(fā)育壩砂為主,砂體厚且橫向分布穩(wěn)定,準(zhǔn)層序組1沉積期和準(zhǔn)層序組3沉積期以灘砂和壩砂間互發(fā)育,砂體厚度相對較小;湖侵域和高位域時期主要發(fā)育碳酸鹽巖灘壩,砂質(zhì)灘壩呈零星狀分布。

[1] 商曉飛,侯加根,董越,等.板橋凹陷同沉積斷層主控的灘壩砂體沉積機(jī)制與分布樣式[J].中國石油大學(xué)學(xué)報(自然科學(xué)版),2014,38(6):32-39. SHANG Xiaofei, HOU Jiagen, DONG Yue, et al. Sedimentary mechanism and distribution pattern of beach-bar sandbodies mainly dominated by contemporaneous faults in Banqiao Sag[J]. Journal of China University of Petroleum (Edition of Natural Science), 2014,38(6):32-39.

[2] 朱筱敏,信荃麟,張晉仁,等.斷陷湖盆灘壩儲集體沉積特征及沉積模式[J].沉積學(xué)報,1994,12(2):20-28. ZHU Xiaomin, XIN Quanlin, ZHANG Jinren, et al. Sedimentary characteristics and models of beach-bar reservoirs in faulted down lacustrine basin[J]. Acta Sedimentologica Sinca, 1994,12(2):20-28.

[3] 陳世悅,楊劍萍,操應(yīng)長,等.惠民凹陷西部下第三系沙河街組兩種灘壩沉積特征[J].煤田地質(zhì)與勘探,2000,28(3):1-4. CHEN Shiyue, YANG Jianping, CAO Yingchang, et al. Sedimentary characteristics of two kinds of beach-bar of Oligocene Shahejie Formation in western Huimin Depression, China[J]. Coal Geology & Exploration, 2000,28(3):1-4.

[4] 楊勇強(qiáng),邱隆偉,姜在興,等.陸相斷陷湖盆灘壩沉積模式:以東營凹陷古近系沙四上亞段為例[J].石油學(xué)報,2011,32(3):417-423. YANG Yongqiang, QIU Longwei, JIANG Zaixing, et al. A depositional pattern if beach bar in continental rift lake basins: a case study on the upper part of the fourth member of Shahejie Formation in the Dongying Sag[J]. Acta Petrolei Sinica, 2011,32(3):417-423.

[5] 劉書會.薄層屬性分析中存在的問題及解決方法:以東營凹陷梁108地區(qū)灘壩砂巖為例[J].油氣地質(zhì)與采收率,2006,13(2):56-58. LIU Shuhui. Problems in thin-layer attribute analysis and solved methods: a case study in beach bar sandstones in Liang 108 area of Dongying Sag[J]. Petroleum Geology and Recovery Efficienfy, 2006,13(2):56-58.

[6] 才巨宏.綜合應(yīng)用波形分析及地震特征反演技術(shù)預(yù)測灘壩砂巖:以博興洼陷梁108地區(qū)為例[J].油氣地質(zhì)與采收率,2005,12(3):42-44. CAI Juhong. Comprehensive application of technology with wave form analysis and seismic characteristic inversion to beach-bar sandstones prediction-taking the area of Liang108 in Boxing Subsag as an example[J]. Petroleum Geology and Recovery Efficiency, 2005,12(3):42-44.

[7] 韓宏偉.薄互層地震波形特征研究:以博興洼陷沙四段灘壩砂為例[J].地學(xué)前緣,2009,16(3):349-355. HAN Hongwei. Research on the characteristics of thin-alternating-bed seismic waveform: a case study of the beach bar sandstones of Es4in Boxing Sag[J]. Earth Science Frontiers, 2009,16(3):349-355.

[8] 黃捍東,張如偉,魏世平.地震非線性隨機(jī)反演方法在陸相薄砂巖儲層預(yù)測中的應(yīng)用[J].石油學(xué)報,2009,30(3):386-390. HUANG Handong, ZHANG Ruwei, WEI Shiping. Research on application of seismic nonlinear random inversion to reservoirs prediction in the thin sandstone of continental deposits[J]. Acta Petrolei Sinica, 2009,30(3):386-390.

[9] 張宏兵,楊長春.正則參數(shù)控制下的波阻抗約束反演[J].地球物理學(xué)報,2003,46(6):827-834. ZHANG Hongbing, YANG Changchun. A constrained impedance inversion method controlled by regularized[J]. Chinese Journal of Geophysics, 2003,46(6):827-834.

[10] 張宏兵,尚作萍,楊長春,等.波阻抗反演正則參數(shù)估計[J].地球物理學(xué)報,2005,48(1):181-188. ZHANG Hongbing, SHANG Zuoping, YANG Changchun, et al. Estimation of regular parameters for the impedance inversion[J]. Chinese Journal of Geophysics, 2005,48(1):181-188.

[11] 張廣智,王丹陽,印興耀.利用MCMC方法估算地震參數(shù)[J].石油地球物理勘探,2011,46(4):605-609. ZHANG Guangzhi, WANG Danyang, YIN Xingyao. Seismic parameter estimation using Markov Chain Monte Carlo Method[J]. OGP, 2011,46(4):605-609.

[12] 朱新玲.馬爾科夫鏈蒙特卡羅方法研究綜述[J].統(tǒng)計與決策,2009(21):151-153. ZHU Xinling. Review on the research of Markov Monte Carlo method[J]. Statistics & Decision,2009 (21):151-153.

[13] 袁紅軍.東營凹陷博興洼陷濱淺湖相灘壩砂巖儲層預(yù)測[J].石油與天然氣地質(zhì),2007,28(4):497-503. YUAN Hongjun. Prediction of beach bar sand reservoirs of shore-shallow lake facies in Boxing subsag of Dongying sag[J]. Oil & Gas Geology, 2007,28(4):497-503.

[14] 郭建卿,林承焰.頻譜分解及地質(zhì)模型反演新方法在灘壩砂沉積特征及發(fā)育模式研究中的應(yīng)用[J].中國石油大學(xué)學(xué)報(自然科學(xué)版),2013,37(3):37-43. GUO Jianqing, LIN Chengyan. Application of new method of spectrum decomposition and seismic inversion in research of beach-bar sand sedimentary characteristics and development model[J]. Journal of China University of Petroleum(Edition of Natural Science), 2013,37(3):37-43.

[15] 商曉飛,侯加根,孫福亭,等.砂質(zhì)灘壩儲集層內(nèi)部結(jié)構(gòu)特征及構(gòu)型模式:以黃驊坳陷板橋油田古近系沙河街組為例[J].石油學(xué)報,2014,35(6):1160-1171. SHANG Xiaofei, HOU Jiagen, SUN Futing, et al. Architectural characteristics and sedimentary models of beach-bar sandstone reservoirs: a case study of Paleogene Shahejie Formation in Banqiao oilfield, Huanghua depression[J]. Acta Petrolei Sinica, 2014,35(6):1160-1171.

[16] 張軍華,劉培金,朱博華,等.灘壩砂儲層地震解釋存在的問題及對策[J].石油地球物理勘探,2014,49(1):167-175. ZHANG Junhua, LIU Peijin, ZHU Bohua, et al. Problems and countermeasures in seismic interpretation of beach bar sandstone reservoirs[J]. OGP, 2014,49(1):167-175.

[17] 蘇永進(jìn),蔣有錄,房新娜,等.博興洼陷古近系沙四上和沙三段油氣成藏的差異性研究[J].中國石油大學(xué)學(xué)報(自然科學(xué)版),2006,30(2):11-15. SU Yongjin, JIANG Youlu, FANG Xinna, et al. Otherness on pool-forming in the upper E2s4and E2s3of Palaeogene in Boxing subsag[J]. Journal of China University of Petroleum (Edition of Natural Science), 2006,30(2):11-15.

[18] 田繼軍,姜在興.東營凹陷沙河街組四段上亞段層序地層特征與沉積體系演化[J].地質(zhì)學(xué)報,2009,83(6):836-846. TIAN Jijun, JIANG Zaixing. Sequence stratigraphy characteristics and sedimentary system evolution of upper Es4in the Dongying depression[J]. Acta Geologica Sinica, 2009,83(6):836-846.

[19] 操應(yīng)長,王健,劉惠民,等.東營凹陷南坡沙四上亞段灘壩砂體的沉積特征及模式[J].中國石油大學(xué)學(xué)報(自然科學(xué)版),2009,33(6):5-10. CAO Yingchang, WANG Jian, LIU Huimin, et al. Sedimentary characteristics and models of beach-bar sandbodies in theupper part of the fourth member of Paleogene in the south slope of Dongying depression[J]. Journal of China University of Petroleum(Edition of Natural Science), 2009,33(6):5-10.

[20] 李國斌,姜在興,陳詩望,等.利津洼陷沙四上亞段灘壩沉積特征及控制因素分析[J].中國地質(zhì),2008,35(5):911-921. LI Guobin, JIANG Zaixing, CHEN Shiwang, et al. Sedimentary characteristics and controlling factors of beach bars in the upper submember of the fourth member of the Shahejie formation in the Lijin subasin[J]. Geology in China, 2008,35(5):911-921.

[22] 楊勇強(qiáng),邱隆偉,姜在興,等.東營凹陷沙四上亞段灘壩物源體系[J].吉林大學(xué)學(xué)報(地球科學(xué)版),2011,41(1):46-53. YANG Yongqiang, QIU longwei, JIANG Zaixing, et al. Beach bar-provenance system on the upper part of fourth member of Shahejie Formation, in Dongying sag[J]. Journal of Jilin University(Earth Science Edition), 2011,41(1):46-53.

(編輯 修榮榮)

HAN Changcheng1, LIN Chengyan1,2, REN Lihua1,2, MA Cunfei1,WEI Ting3, ZHANG Xianguo1,2, SUN Zhifeng1

(1.SchoolofGeosciencesinChinaUniversityofPetroleum,Qingdao266580,China;2.KeyLaboratoryofReservoirGeologyinShandongProvince,Qingdao266580,China;3.ShixiOilProductionPlant,XinjiangOilfieldCompany,PetroChina,Karamay834000,China)

2016-05-13

國家重大科技專項(2017ZX05009-001);國家自然科學(xué)基金項目(41672129)

韓長城(1984-)，男，博士研究生，研究方向?yàn)橛蜌獠亻_發(fā)地質(zhì)。E-mail:517892849@qq.com。

1673-5005(2017)02-0060-10

10.3969/j.issn.1673-5005.2017.02.007

TE 122.1

A

韓長城,林承焰,任麗華,等. 地震波形指示反演在東營凹陷王家崗地區(qū)沙四上亞段灘壩砂的應(yīng)用[J]. 中國石油大學(xué)學(xué)報(自然科學(xué)版),2017,41(2):60-69.