渤中凹陷西南部烴源巖TOC含量預(yù)測(cè)

王 祥 馬勁風(fēng)* 王德英 王震亮 張新濤 王飛龍

(①西北大學(xué)地質(zhì)學(xué)系，陜西西安 710069； ②西北大學(xué)地質(zhì)學(xué)系二氧化碳捕集與封存技術(shù)國(guó)家地方聯(lián)合工程研究中心，陜西西安 710069； ③西北大學(xué)大陸動(dòng)力學(xué)國(guó)家重點(diǎn)實(shí)驗(yàn)室，陜西西安 710069；④中國(guó)海油天津分公司勘探開(kāi)發(fā)研究院，天津 300452)

0 引言

烴源巖的總有機(jī)碳(TOC)含量是有機(jī)質(zhì)豐度評(píng)價(jià)的重要指標(biāo)。巖心測(cè)試僅能獲得離散的TOC含量數(shù)據(jù)，尤其是井位稀少或者取心較少的勘探區(qū)塊，僅靠地球化學(xué)分析手段難以對(duì)烴源巖的生烴潛力和資源量進(jìn)行三維評(píng)價(jià)[1-2]。受構(gòu)造演化和沉積作用等地質(zhì)因素的影響，烴源巖在縱向和橫向上都可能具有很強(qiáng)的非均質(zhì)性[3-4]，所以如何識(shí)別優(yōu)質(zhì)的烴源巖，研究烴源巖的空間分布規(guī)律，降低勘探風(fēng)險(xiǎn)，是一個(gè)長(zhǎng)期的研究重點(diǎn)。

地球物理測(cè)井方法因其較高的縱向分辨率和較好的連續(xù)性，在油氣勘探和開(kāi)發(fā)階段都發(fā)揮著至關(guān)重要的作用。自1945年，學(xué)者們開(kāi)始根據(jù)測(cè)井曲線研究烴源巖，即分別使用了自然伽馬(GR)、聲波時(shí)差(DT)、電阻率等測(cè)井曲線及交會(huì)圖分析測(cè)井參數(shù)與烴源巖有機(jī)質(zhì)豐度的對(duì)應(yīng)關(guān)系[5-10]。Passey等[11]基于Archie公式推導(dǎo)了基于電阻率曲線和孔隙度曲線的定量計(jì)算TOC含量曲線的方法——ΔlogR法，烴源巖的TOC含量預(yù)測(cè)自此進(jìn)入了定量化的時(shí)代。在中國(guó)針對(duì)陸相地層對(duì)該方法做了諸多分析和改進(jìn)，相繼提出了簡(jiǎn)化ΔlogR法[12]、廣義ΔlogR法[13]、變系數(shù)ΔlogR法[14]等。多元回歸法[15-17]也是一種應(yīng)用廣泛的預(yù)測(cè)TOC的測(cè)井方法，該方法通過(guò)交會(huì)圖優(yōu)選與TOC相關(guān)度高的測(cè)井參數(shù)，進(jìn)行多元擬合，進(jìn)而實(shí)現(xiàn)TOC含量實(shí)測(cè)數(shù)據(jù)的曲線化，但公式具有很強(qiáng)的地區(qū)特殊性，且ΔlogR方法的物理意義并不明確。近年來(lái)，隨著人工智能技術(shù)的興起，也產(chǎn)生了很多測(cè)井資料與人工智能相結(jié)合計(jì)算TOC含量的新思路和新方法[18-20]。

雖然測(cè)井拓展預(yù)測(cè)的TOC含量曲線可以直接反映井位置已鉆遇烴源巖質(zhì)量的優(yōu)劣，但是在勘探初期階段，由于鉆井?dāng)?shù)量較少或未鉆遇目的烴源巖層，因此如何分析井間和深層未鉆遇的烴源巖分布，是備受關(guān)注的問(wèn)題。由于地震資料具有較高的橫向分辨率，在識(shí)別烴源巖的工作中發(fā)揮著重要的作用。在研究工作的早期，人們更加關(guān)注地震相和地震速度在烴源巖層段的響應(yīng)特征，忽視了對(duì)巖石的彈性參數(shù)、地震資料屬性等與烴源巖有機(jī)質(zhì)的關(guān)系的研究[1-2,21-22]。隨著研究的深入，逐步利用地震屬性精細(xì)地研究烴源巖，由此衍生出了使用地震多屬性分析技術(shù)進(jìn)行TOC含量預(yù)測(cè)的方法。這些方法首先使用測(cè)井方法計(jì)算TOC含量曲線；再提取井旁道的地震屬性，以TOC含量曲線為因變量，地震屬性為自變量，優(yōu)選與TOC含量相關(guān)度高的地震屬性進(jìn)行多元擬合，實(shí)現(xiàn)TOC含量的地震預(yù)測(cè)。針對(duì)TOC含量的多屬性預(yù)測(cè)方法，研究熱點(diǎn)主要集中在更好地建立多屬性預(yù)測(cè)TOC含量的井—震預(yù)測(cè)模型[23-30]。近年還出現(xiàn)了使用BP神經(jīng)網(wǎng)絡(luò)技術(shù)的多屬性預(yù)測(cè)TOC含量的方法[29]。此外，地震反演方法也在預(yù)測(cè)TOC含量方面具有廣泛的應(yīng)用。L?seth等[23]分析了波阻抗等彈性參數(shù)與TOC含量的交會(huì)圖，為后續(xù)尋找彈性參數(shù)與TOC含量的關(guān)系提供了理論依據(jù)，從而間接證明了地震反演TOC含量的可行性。反演方法較地震屬性方法更為豐富，主要方法有三種： ①根據(jù)密度與TOC含量的高相關(guān)性，使用疊前反演的密度體計(jì)算TOC含量體[31-35]； ②使用分頻反演技術(shù)預(yù)測(cè)TOC含量數(shù)據(jù)，由于此過(guò)程考慮到了不同厚度下振幅與頻率的關(guān)系(AVF)信息，因此反演結(jié)果的可信度更高[36]； ③結(jié)合巖石物理交會(huì)圖分析，使用疊后波阻抗反演方法預(yù)測(cè)TOC含量[37]。

雖然地震預(yù)測(cè)方法難以得到電阻率、自然伽馬等與TOC含量密切相關(guān)的測(cè)井參數(shù)，其縱向分辨率不如測(cè)井方法，但是根據(jù)優(yōu)選的地震屬性、彈性參數(shù)反演體可預(yù)測(cè)TOC含量的空間分布。這發(fā)揮了地震資料橫向分辨率較高的優(yōu)勢(shì)，也充分利用了井間的地震信息，相比于測(cè)井預(yù)測(cè)TOC含量曲線，進(jìn)行井間插值獲得TOC含量數(shù)據(jù)體的方法更加可靠。

近年來(lái)，伴隨著巖石物理技術(shù)的快速發(fā)展，開(kāi)始考慮TOC含量對(duì)巖石彈性參數(shù)和泥頁(yè)巖各向異性參數(shù)的影響，并且建立了富含有機(jī)質(zhì)的烴源巖巖石物理模型[38-43]。精確的正演模型是獲得良好反演結(jié)果的基礎(chǔ)，通過(guò)研究富含有機(jī)質(zhì)烴源巖的巖石物理性質(zhì)，尋找TOC含量與烴源巖彈性參數(shù)的本質(zhì)聯(lián)系，是未來(lái)利用地震資料定量識(shí)別烴源巖的基礎(chǔ)[44]。

本文以渤海灣盆地渤中凹陷西南部三維疊前地震資料和測(cè)井資料為基礎(chǔ)，結(jié)合構(gòu)造、沉積、巖石物理、實(shí)測(cè)TOC含量數(shù)據(jù)等資料，優(yōu)選了本研究區(qū)的TOC含量測(cè)井預(yù)測(cè)方法(拓展ΔlogR方法)；并通過(guò)巖石物理交會(huì)分析，建立考慮基于縱波速度(VP)、密度(DEN)的地震預(yù)測(cè)TOC含量模型，使用疊前反演方法分層(東營(yíng)組、沙河街組)計(jì)算研究區(qū)TOC含量反演體；最終利用沉積相驗(yàn)證反演結(jié)果，并通過(guò)反演井、檢驗(yàn)井井旁道TOC含量、測(cè)井預(yù)測(cè)TOC含量與實(shí)測(cè)含量TOC進(jìn)行對(duì)比，驗(yàn)證了該TOC含量空間評(píng)價(jià)方法的可靠性。

1 研究區(qū)地質(zhì)背景

渤中凹陷面積近1×104km2，是渤海灣盆地中最大的凹陷之一，也是中國(guó)近海地區(qū)最重要的一個(gè)油氣勘探凹陷。研究區(qū)位于渤海灣盆地渤中凹陷西南部，北部為沙壘田凸起，西部為埕北低凸起，東南部為渤南低凸起，渤中19-6構(gòu)造脊位于研究區(qū)中部(圖1)[45]。前人通過(guò)構(gòu)造演化特征、沉積作用，結(jié)合鉆、測(cè)井資料，認(rèn)為該區(qū)主要發(fā)育沙三段、沙一二段、東營(yíng)組下段三套湖相烴源巖。伴隨多期的構(gòu)造運(yùn)動(dòng)，不同沉積時(shí)期的沉積中心發(fā)生了較為明顯的遷移。凹陷自沙三段沉積時(shí)期至東二下段沉積時(shí)期，總體經(jīng)歷了沙三段(快速斷陷沉降時(shí)期，凹陷形成了多個(gè)沉降區(qū))——沙一二段(沙三組沉積末期抬升剝蝕后，重新接受沉積，湖盆面積變大，水體變淺)——東營(yíng)組沉積末期(快速斷陷沉降時(shí)期，沉積了渤中凹陷特有的一套湖相泥巖)的三個(gè)階段。研究區(qū)沙三段主要發(fā)育半深湖—深湖沉積，沉積中心位于研究區(qū)中南部；沙一二段主要發(fā)育濱淺湖沉積，沉積中心位于研究區(qū)中部；東三段沉積時(shí)期，研究區(qū)北部主要發(fā)育半深湖—深湖沉積，南部主要發(fā)育濱淺湖沉積；東二下段沉積時(shí)期，研究區(qū)主要發(fā)育辮狀河三角洲沉積[45-46]。

研究區(qū)內(nèi)渤中19-6構(gòu)造脊及其周?chē)l(fā)育了一系列低潛山群，多套湖相烴源巖披覆在潛山之上，鮮見(jiàn)發(fā)育優(yōu)質(zhì)烴源巖的凹陷深部鉆井和測(cè)試資料[47]。

圖1 渤中凹陷地質(zhì)圖(據(jù)文獻(xiàn)[46]修改)及工區(qū)內(nèi)鉆井分布圖

2 研究區(qū)烴源巖TOC含量預(yù)測(cè)方法

綜合研究區(qū)地質(zhì)背景，在鉆、測(cè)井資料缺乏的情況下預(yù)測(cè)深層良好烴源巖成為了亟需解決的問(wèn)題。研究區(qū)內(nèi)覆蓋了高精度的三維地震資料，筆者結(jié)合西南部較為豐富的鉆、測(cè)井和巖心實(shí)測(cè)數(shù)據(jù)，優(yōu)選了測(cè)井方法和疊前地震反演方法進(jìn)行烴源巖TOC含量的三維預(yù)測(cè)，為凹陷深部的烴源巖研究提供理論依據(jù)。具體流程見(jiàn)圖2。

圖2 研究區(qū)烴源巖TOC含量地球物理預(yù)測(cè)方法流程

2.1 巖石物性交會(huì)分析

根據(jù)前人研究，良好的陸相烴源巖具有“三高一低”的特征，即： 高GR、高DT、高RD和較低的DEN[48]。通過(guò)建立巖石物性交會(huì)圖板，分析對(duì)TOC含量敏感的彈性參數(shù)。本次研究使用渤中凹陷西南部8口井(圖1)的TOC含量測(cè)試數(shù)據(jù)與測(cè)井?dāng)?shù)據(jù)進(jìn)行分析。首先將DT轉(zhuǎn)化為VP，再分別繪制VP、DEN、RD、GR與TOC含量實(shí)測(cè)數(shù)據(jù)建立交會(huì)圖。經(jīng)過(guò)對(duì)比，發(fā)現(xiàn)VP與TOC含量相關(guān)性最高，其次為DEN，RD、GR與TOC含量相關(guān)性較低(圖3)。因此，本研究區(qū)烴源巖的彈性參數(shù)(Vp、DEN)與TOC含量存在較高的相關(guān)性，可以通過(guò)地震反演預(yù)測(cè)TOC含量的空間分布。

2.2 地球物理測(cè)井方法

2.2.1 拓展ΔlogR方法

調(diào)研發(fā)現(xiàn)，通過(guò)ΔlogR系列方法在縱向上將TOC含量實(shí)測(cè)值擴(kuò)展為T(mén)OC含量曲線比較可靠，可實(shí)現(xiàn)井—震聯(lián)合預(yù)測(cè)TOC含量的目的。考慮到研究區(qū)復(fù)雜的地質(zhì)背景和強(qiáng)烈的壓實(shí)作用，本文采用一種拓展的ΔlogR方法，即以Passey等[11]提出的ΔlogR方法為基礎(chǔ)，引入基于GR和DEN曲線預(yù)測(cè)TOC含量

(1)

TOC=ΔlogR(a+b×DEN+c×lgGR)+d

(2)

式中：ΔlogR為孔隙度曲線(本文選用DT計(jì)算)與電阻率曲線RD反向疊合后的幅度差；RD、R基線、Δt、Δt基線分別為RD值、RD基線值、DT值、DT基線值；a、b、c、d為模型擬合系數(shù)。分析式(1)和式(2)可以看出，當(dāng)不考慮GR和DEN的影響(即b、c為0)時(shí)，該公式可化簡(jiǎn)為Passey等提出的初始ΔlogR公式，即此公式引入了更多與TOC敏感的測(cè)井參數(shù)GR和DEN，拓展了測(cè)井曲線預(yù)測(cè)TOC含量的參數(shù)范圍[49]。

2.2.2 多元回歸法

根據(jù)測(cè)井參數(shù)交會(huì)圖分析(圖3)，優(yōu)選與TOC含量相關(guān)度高的測(cè)井參數(shù)，以實(shí)測(cè)TOC含量為因變量、優(yōu)選的測(cè)井參數(shù)為自變量進(jìn)行多元擬合

TOC=e+f×DT+g×RD+h×DEN

(3)

式中e、f、g、h為模型擬合系數(shù)。

2.3 地震疊前反演方法

地震疊前反演僅能獲得IS(橫波阻抗)、IP(縱波阻抗)、VP、VS(橫波速度)、VP/VS(縱橫波速度比)等彈性參數(shù)，不能直接反演得到DEN和GR等測(cè)井參數(shù)。經(jīng)過(guò)巖石物理交會(huì)分析和對(duì)區(qū)域地質(zhì)背景的研究，發(fā)現(xiàn)VP和DEN與TOC含量的相關(guān)性最高，因此使用TOC含量實(shí)測(cè)值與彈性參數(shù)(VP、DEN)擬合反演公式

TOC地震反演=k+m×VP+n×DEN

(4)

圖3 研究區(qū)實(shí)測(cè)TOC含量與測(cè)井參數(shù)VP(a)、DEN(b)、RD(c)、GR(d)的交會(huì)圖

式中k、m、n為擬合系數(shù)。使用研究區(qū)高精度三維疊前地震資料進(jìn)行疊前同時(shí)反演，得到VP反演體和DEN反演體，再代入式(4)，獲得TOC反演體。

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

3.1 測(cè)井預(yù)測(cè)TOC含量

本文采用拓展ΔlogR方法和多元回歸法進(jìn)行TOC含量測(cè)井預(yù)測(cè)，并進(jìn)行對(duì)比優(yōu)選。研究區(qū)發(fā)育多套湖相烴源巖，東營(yíng)組與沙河街組的烴源巖在沉積背景上存在較大差異，直接影響了不同層位烴源巖的巖石物性的差異，因此使用研究區(qū)的TOC含量實(shí)測(cè)數(shù)據(jù)與測(cè)井?dāng)?shù)據(jù)進(jìn)行分層擬合計(jì)算，并使用區(qū)內(nèi)兩口未參與模型建立的井(B井、E井)進(jìn)行驗(yàn)證。

東營(yíng)組和沙河街組(圖4a和圖4c)拓展ΔlogR方法擬合公式分別為

(5)

(6)

式中R表示相關(guān)系數(shù)。

通過(guò)B井和E井檢驗(yàn)，該方法計(jì)算的B井TOC含量平均相對(duì)誤差為13.58%，E井平均相對(duì)誤差為14.89%。

東營(yíng)組和沙河街組(圖4b和圖4d)多元回歸法擬合公式分別為

(7)

(8)

通過(guò)B井和E井的驗(yàn)證，該方法計(jì)算的B井TOC含量平均相對(duì)誤差為24.62%，E井為27.67%。

通過(guò)對(duì)比以上數(shù)據(jù)，本次研究選擇拓展ΔlogR方法預(yù)測(cè)TOC含量曲線。

圖4 拓展ΔlogR方法與多元回歸法應(yīng)用誤差對(duì)比圖

3.2 地震預(yù)測(cè)TOC含量

采用疊前同時(shí)反演預(yù)測(cè)研究區(qū)烴源巖TOC含量。根據(jù)研究區(qū)的6口井(A、C、D、F、G、H井)數(shù)據(jù)進(jìn)行反演，以B井、E井作為檢驗(yàn)井，檢驗(yàn)反演結(jié)果。為了驗(yàn)證反演結(jié)果的可靠性，抽取TOC含量實(shí)測(cè)數(shù)據(jù)較豐富的反演井(A井、C井、D井)和驗(yàn)證井(B井、E井)井旁TOC含量剖面圖與拓展ΔlogR方法預(yù)測(cè)的TOC含量曲線、TOC含量實(shí)測(cè)數(shù)據(jù)進(jìn)行對(duì)比。

首先根據(jù)巖石物理交會(huì)圖分析結(jié)果，得到東營(yíng)組和沙河街組地震分層反演TOC含量的計(jì)算公式分別為

(9)

(10)

使用HRS反演軟件導(dǎo)入測(cè)井曲線、連片的三維疊前地震數(shù)據(jù)、層位以及TOC預(yù)測(cè)曲線；制作合成地震記錄，進(jìn)行井震標(biāo)定，保證在井位置獲得可靠的時(shí)深關(guān)系；再建立疊前反演模型，進(jìn)行反演分析； 進(jìn)而反演得到VP反演體和DEN反演體。對(duì)比反演的VP、VS和DEN與E井剖面可以看出，二者對(duì)應(yīng)關(guān)系較好，證實(shí)了反演的一手資料——彈性參數(shù)的可靠性(圖5)。最后使用式(9)、式(10)計(jì)算得到東營(yíng)組和沙河街組的TOC含量反演體。并將反演結(jié)果與測(cè)井預(yù)測(cè)TOC含量曲線、實(shí)測(cè)TOC含量進(jìn)行對(duì)比(圖6～圖8)。對(duì)比可見(jiàn)，反演的TOC含量、測(cè)井方法預(yù)測(cè)的TOC含量與實(shí)測(cè)的TOC含量吻合度較高。

圖5 過(guò)E井(檢驗(yàn)井)疊前VP(a)、VS(b)、DEN(c)反演剖面

圖7 過(guò)B井(檢驗(yàn)井)TOC含量反演剖面

圖8 過(guò)E井(檢驗(yàn)井)TOC含量反演剖面

3.3 井震預(yù)測(cè)結(jié)果對(duì)比

由于地震數(shù)據(jù)的縱向分辨率不如測(cè)井，因此在對(duì)比反演結(jié)果與巖心實(shí)測(cè)TOC含量時(shí)，將井旁道TOC含量曲線、測(cè)井預(yù)測(cè)TOC含量曲線分為東二段、東三段、沙一二段、沙三段，分別計(jì)算各層的平均值，并與實(shí)測(cè)巖心TOC含量平均值進(jìn)行對(duì)比(圖9)。

圖9 井位置巖心實(shí)測(cè)、地震預(yù)測(cè)、測(cè)井預(yù)測(cè)TOC含量統(tǒng)計(jì)

以下幾點(diǎn)需要特別說(shuō)明：

(1)東三段誤差較大的原因：此段相比于其他層位比較薄，TOC含量實(shí)測(cè)樣本較少，因此測(cè)量結(jié)果的平均值誤差較大，故東三段的TOC含量評(píng)價(jià)應(yīng)盡量以沉積相圖和TOC含量平面圖的對(duì)比為主；

(2)B井、D井僅有東二下段的實(shí)測(cè)TOC含量數(shù)據(jù)，因此圖9中除東二下段外，未對(duì)B井、D井進(jìn)行誤差對(duì)比。

通過(guò)對(duì)比地震井旁反演結(jié)果、測(cè)井預(yù)測(cè)結(jié)果與TOC含量實(shí)測(cè)數(shù)據(jù)，總體上可以看出：

(1)測(cè)井方法在縱向上預(yù)測(cè)的TOC含量相對(duì)誤差較小，準(zhǔn)確性更高；

(2)沙河街組烴源巖TOC含量高于東營(yíng)組烴源巖。

3.4 沉積相驗(yàn)證

基于VP反演體和DEN反演體計(jì)算得到TOC含量反演體后，再根據(jù)地震層位分別求取東二下段、東三段、沙一二段、沙三段的TOC含量平均值，得到烴源巖層的TOC含量平面圖，并與各層的沉積相圖進(jìn)行對(duì)比(圖10)。

3.4.1 東二下段(E3d2L)

東二下段沉積時(shí)期，研究區(qū)中部主要發(fā)育辮狀河三角洲沉積，北部為半深湖—深湖沉積和湖底扇沉積，南部為濱淺湖沉積，水體動(dòng)蕩的沉積背景不利于形成大范圍的優(yōu)質(zhì)烴源巖，研究區(qū)北部TOC含量約為0.6%～1.5%，南部約為0.3%～0.6%，研究區(qū)北部TOC含量較高，南部較低，總體呈現(xiàn)由北向南TOC含量逐漸降低的趨勢(shì)(圖10a)。

3.4.2 東三段(E3d3)

東三段沉積時(shí)期，研究區(qū)北部為半深湖—深湖沉積，南部為濱淺湖沉積，該沉積時(shí)期從北到南水體逐漸變淺，研究區(qū)北部TOC含量約為1.5%～2.2%，南部約為0.5%～1.5%，研究區(qū)北部TOC含量較高，南部較低，總體也呈現(xiàn)由北向南TOC含量逐漸降低的趨勢(shì)(圖10b)。

3.4.3 沙一二段(E2s1+2)

沙一二段沉積時(shí)期，研究區(qū)主要為濱淺湖沉積，該沉積時(shí)期水體較淺，因此TOC含量相較沙三段沉積時(shí)期偏低，A、C、E井區(qū)及其附近區(qū)域?yàn)樵摃r(shí)期的沉積中心，TOC含量約為2.0%～2.5%，其余濱淺湖沉積地區(qū)TOC含量約為0.5%～2.0%，總體呈現(xiàn)以A、B、C、E井區(qū)為中心、TOC含量向外逐漸降低的趨勢(shì)(圖10c)。

圖10 工區(qū)TOC含量平面分布(左)與沉積相(右)對(duì)比

3.4.4 沙三段(E2s3)

沙三段沉積時(shí)期，研究區(qū)主要為半深湖—深湖沉積，沉積中心在A、C井區(qū)附近及其南部地區(qū)，該沉積時(shí)期水體較深，有利于優(yōu)質(zhì)烴源巖的形成，因此TOC含量平面圖中A、C井區(qū)及南部部分地區(qū)的TOC含量較高，約為2.5%～2.8%，南部部分無(wú)井地區(qū)的TOC含量可達(dá)3.0%～3.5%；其余半深湖—深湖地區(qū)約為1.5%～2.5%，北部和西南部的濱淺湖地區(qū)約為0.5%～1.5%，總體也呈現(xiàn)以A、C井區(qū)及其西部、南部地區(qū)為中心，TOC含量向外逐漸降低的趨勢(shì)(圖10d)。

4 結(jié)論

本文使用拓展ΔlogR方法、疊前同時(shí)反演方法對(duì)渤海灣盆地渤中凹陷西南部某三維地震工區(qū)的烴源巖TOC含量進(jìn)行了預(yù)測(cè)，得到以下結(jié)論。

(1)渤中凹陷西南部烴源巖使用基于GR和DEN的一種拓展ΔlogR方法，預(yù)測(cè)的TOC含量誤差較小，優(yōu)于多元回歸法；該地區(qū)地質(zhì)背景復(fù)雜，因此應(yīng)采取多種測(cè)井方法預(yù)測(cè)，進(jìn)行橫向?qū)Ρ葍?yōu)選。

(2)研究區(qū)烴源巖TOC含量與VP、DEN相關(guān)度最高，因此通過(guò)VP、DEN建立TOC含量反演公式。以研究區(qū)測(cè)井?dāng)?shù)據(jù)、TOC含量實(shí)測(cè)數(shù)據(jù)、三維疊前地震數(shù)據(jù)為基礎(chǔ)，進(jìn)行疊前同時(shí)反演得到VP、DEN反演體，再計(jì)算得到TOC含量反演體。井震預(yù)測(cè)誤差分析和沉積相檢驗(yàn)證明了該三維預(yù)測(cè)方法是可靠的。

(3)總體來(lái)說(shuō)，研究區(qū)發(fā)育上下疊置的多套烴源巖，且橫向和縱向非均質(zhì)性明顯。東營(yíng)組TOC含量均呈現(xiàn)北高南低的趨勢(shì)，東二下段TOC含量整體偏低，東三段整體偏高；沙河街組TOC含量整體較高，為研究區(qū)的主力烴源巖層，TOC含量呈現(xiàn)自沉積中心向外逐漸降低的趨勢(shì)，沙一二段A、B、C、E井區(qū)及其附近區(qū)域、沙三段A、C井區(qū)及其西部、南部部分地區(qū)TOC含量最高，具有良好的生烴潛力。