楊 倉，鄒 荃，客偉利

(1.中國石油尼羅河公司，北京 100034；2.中國石油勘探開發(fā)研究院非洲所，北京 100083)

約束稀疏脈沖儲層反演技術(shù)在Sufyan凹陷的應(yīng)用

楊 倉1，鄒 荃2，客偉利2

(1.中國石油尼羅河公司，北京 100034；2.中國石油勘探開發(fā)研究院非洲所，北京 100083)

Sufyan凹陷深層AG組砂體存在厚度較薄、橫向變化快、儲層預(yù)測難度較大的問題，而約束稀疏脈沖反演作為儲層預(yù)測廣泛使用的方法，具有反演精度高、多解性少的特點(diǎn)，能夠直觀地刻畫巖性的空間分布規(guī)律。本文運(yùn)用約束稀疏脈沖儲層反演技術(shù)，通過地震資料分析、測井曲線標(biāo)準(zhǔn)化處理、精細(xì)井震標(biāo)定和子波提取等一系列工作，優(yōu)選反演參數(shù)，對Sufyan凹陷深層AG組進(jìn)行了儲層預(yù)測，為該區(qū)的儲層評價(jià)提供了有力支撐。實(shí)例應(yīng)用表明，該技術(shù)處理后的預(yù)測效果與鉆井資料區(qū)域沉積特征相吻合，能夠?yàn)楹罄m(xù)開發(fā)提供依據(jù)。

Sufyan凹陷；AG組；約束稀疏脈沖反演；儲層預(yù)測

隨著高分辨率地震技術(shù)的逐漸成熟，地震資料的應(yīng)用已經(jīng)不局限于基于地震同相軸追蹤進(jìn)行層位解釋的常規(guī)方法。借助于測井資料和地質(zhì)認(rèn)識，利用波阻抗反演提高地震資料薄儲層識別能力，已經(jīng)成為地震資料解釋和儲層預(yù)測的重要技術(shù)手段[1]。由于國外優(yōu)秀商業(yè)軟件的引入，反演程序的可操作性大大提高，然而反演的效果卻是因人而異，存在多解性。其原因是對待反演資料的研究程度不同，對反演軟件的基本原理和適用范圍掌握不夠，進(jìn)而影響到二者的有機(jī)結(jié)合，不能產(chǎn)生良好的反演效果[2]。因此本文運(yùn)用約束稀疏脈沖反演的方法，在充分掌握反演原理及實(shí)際工區(qū)資料的基礎(chǔ)上，對Sufyan凹陷三維區(qū)進(jìn)行了儲層預(yù)測研究，得到了很好的效果。

1 約束稀疏脈沖反演技術(shù)

約束稀疏脈沖反演(Constrained Sparse spike inversion, CSSI)技術(shù)是目前波阻抗反演技術(shù)中廣泛使用的一種方法。其原理是用地震道的振幅產(chǎn)生波阻抗模型，采用一個(gè)快速的趨勢約束脈沖反演算法，同地震解釋層位和測井約束控制波阻抗的趨勢和幅值范圍，脈沖算法產(chǎn)生了寬帶結(jié)果，恢復(fù)了缺失的部分低頻和高頻成分[3]。

約束稀疏脈沖反演流程如圖1所示，反演軟件為Jason軟件的Invertrace模塊。

圖1 約束稀疏脈沖反演流程Fig.1 CSSI flow diagram

反演的關(guān)鍵是低頻模型的建立、子波提取，約束條件，以及質(zhì)量控制參數(shù)λ的選取。除了要建立合理的低頻模型和提取高質(zhì)量的子波外，在約束稀疏脈沖反演處理中，還加入了多口井的波阻抗趨勢硬約束、地質(zhì)構(gòu)造框架模型控制及地震數(shù)據(jù)軟約束等，軟約束主要控制橫向變化，測井資料硬約束控制縱向變化。

在構(gòu)造研究的基礎(chǔ)上，針對Sufyan凹陷的儲層特征，結(jié)合勘探實(shí)際需要，對Sufyan三維區(qū)Abu Gabra組二段開展了地震反演及儲層預(yù)測相關(guān)工作。實(shí)際工作中，利用地質(zhì)、測井物探資料，開展儲層地質(zhì)分析、儲層地震響應(yīng)特征分析，在綜合分析的基礎(chǔ)上，進(jìn)行反演適應(yīng)方法試驗(yàn)，優(yōu)選適合本區(qū)地質(zhì)特點(diǎn)的反演方法——約束稀疏脈沖反演，開展全區(qū)地震反演。

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

2.1 工區(qū)基本概況及需解決的問題

Muglad盆地位于非洲板塊中部(圖2)，是在穩(wěn)定的前寒武系基底之上發(fā)育起來的中、新生代裂谷盆地[4-5]。其是在中非剪切帶的右旋剪切應(yīng)力場背景下拉張形成的，屬于差異應(yīng)力引起的被動(dòng)裂谷盆地[6-7]。

Sufyan凹陷是Muglad盆地西北部的一個(gè)二級構(gòu)造單元，位于蘇丹共和國的南部，呈菱形、東西向展布。自白堊紀(jì)起，受控于中非斷裂帶剪切運(yùn)動(dòng)形成的張剪性構(gòu)造應(yīng)力背景，Sufyan凹陷開始斷陷，形成以白堊系為主體的斷陷盆地[8-9]。

圖2 Sufyan凹陷區(qū)域構(gòu)造位置Fig.2 Sufyan sag regional structural location

目前，隨著Sufyan凹陷勘探程度的不斷提高，2013—2016年圍繞Sufyan凹陷中央構(gòu)造帶部署的預(yù)探井和評價(jià)井共14口，其中有2口井失利，其余均獲得商業(yè)油流，但是主要目的層AG2組單砂體厚度僅有3～10 m，儲層標(biāo)定及地震識別難度大。本次選用約束稀疏脈沖反演對儲層展布特征進(jìn)行定量預(yù)測，大大地提高了地質(zhì)人員對儲層乃至油氣分布特征的認(rèn)識程度，降低了投資風(fēng)險(xiǎn)。

2.2 約束稀疏脈沖反演技術(shù)的應(yīng)用

(1)對測井曲線進(jìn)行標(biāo)準(zhǔn)化處理。

高質(zhì)量的測井資料是定量化油藏描述項(xiàng)目研究的基本條件。所加載的井資料包括P波聲波、密度、伽馬曲線，以及巖石物理解釋的結(jié)果，如孔隙度、泥質(zhì)含量曲線等[10]。聲波和密度曲線能生成聲波阻抗信息。利用泥質(zhì)含量曲線，以50%為門檻值區(qū)分巖性。為了保持反演結(jié)果的橫向一致性，需要對所有測井曲線進(jìn)行標(biāo)準(zhǔn)化處理。

標(biāo)準(zhǔn)化處理的主要目標(biāo)是聲波、密度和伽馬曲線，目的是消除井間由于施工年度、儀器設(shè)備刻度精度、鉆井液、環(huán)境因素等環(huán)節(jié)引起的系統(tǒng)誤差，使該區(qū)所有井曲線的整體特點(diǎn)與地質(zhì)特征更加吻合，同時(shí)標(biāo)準(zhǔn)化后也使反演基礎(chǔ)數(shù)據(jù)建立在統(tǒng)一的基準(zhǔn)之上[11]。

采用的方法是利用交匯單井波阻抗曲線的方法[12]，對測井曲線進(jìn)行標(biāo)準(zhǔn)化處理。通過直方圖對每條曲線確定極大值、極小值，通過線性回歸方程對測井資料進(jìn)行重新刻度。在獲得了縱向質(zhì)量可靠、橫向具有一致性的測井資料后，才可以開展標(biāo)定、反演等工作，所得到的結(jié)果才是可靠、真實(shí)的。

本次研究對該區(qū)13口井的測井曲線進(jìn)行處理，消除了該區(qū)測井曲線因不同年限、不同儀器、不同采集處理參數(shù)所引起的差異，使所有井的同種曲線數(shù)值在量綱上一致起來(圖3)。

圖3 測井曲線標(biāo)準(zhǔn)化處理對比圖Fig.3 Standardization of well logging curve processing contrast

(2)儲層地球物理特征分析。

為了更好地利用測井資料對儲層進(jìn)行分析，需對測井曲線進(jìn)行儲層地球物理特征分析，確定反映儲層地球物理響應(yīng)特征的曲線，并做各種交會分析，確定儲層特征識別參數(shù)；對儲層的物性做定量分析，建立測井屬性與儲層巖性、物性、電性和含油性的關(guān)系。儲層地球物理特征分析主要采用交會圖分析法，通過曲線交會來分析地下儲層及油氣藏的敏感參數(shù)。

眾所周知，阻抗是聯(lián)系地震數(shù)據(jù)和儲層屬性的橋梁，因此直方圖和相關(guān)性分析是地震反演的關(guān)鍵步驟。本次研究主要以直方圖的方式分析了阻抗與巖性以及阻抗與聲波(DT)、密度(DEN)的相關(guān)性。從AG1、AG2和AG3儲層物性直方圖可以看出，縱波阻抗可以用來直接區(qū)分砂巖和泥巖(圖4)。

圖4 AG組儲層阻抗直方圖Fig.4 Impedance histogram of AG formation reservoir wave

(3)儲層精細(xì)標(biāo)定與子波提取。

層位標(biāo)定及子波反演是聯(lián)系地震和測井?dāng)?shù)據(jù)的橋梁，在地震反演中占有重要地位。對基于模型的反演，地震子波更是聯(lián)系地震記錄與初始模型的橋梁，模型反演結(jié)果與地震記錄、初始模型、地震子波密切相關(guān)。子波是模型反演的關(guān)鍵因素，子波與模型反射系數(shù)褶積產(chǎn)生合成地震記錄，合成地震記錄與實(shí)際地震資料的誤差最小是終止迭代的約束條件[13]。

層位標(biāo)定的好壞直接影響子波反演的結(jié)果，而子波的正確性又對層位的準(zhǔn)確標(biāo)定產(chǎn)生制約，只有通過子波反演和層位標(biāo)定交互迭代才能獲取最佳標(biāo)定和最佳子波。在子波反演過程中，保持聲波曲線標(biāo)定好的時(shí)深關(guān)系不變，以目的層的合成記錄與井旁道地震波組特征有較好的對應(yīng)關(guān)系為原則，采用最小平方法提取一個(gè)初始子波，然后利用廣義線性反演方法，迭代求取最佳子波。

結(jié)合地震解釋層位和測井地質(zhì)分層，采用處理解釋兩套系統(tǒng)相互驗(yàn)證，對工區(qū)內(nèi)10口井進(jìn)行了精細(xì)的合成記錄標(biāo)定。標(biāo)定工作依據(jù)波形、波組和能量最大相關(guān)為原則，本次先用ricker子波標(biāo)定制作合成記錄，然后提取井旁道子波精細(xì)標(biāo)定。正確的合成記錄在獲得與地震記錄較好相關(guān)性的同時(shí)，還必須保證井的時(shí)深關(guān)系的合理性。從標(biāo)定情況來看，合成記錄與井旁地震道有良好相關(guān)性，地震解釋層位和地質(zhì)分層有很好的匹配(圖5)。

(4)子波估算。

子波估算是所有地震反演的核心，估算地震子波波形的準(zhǔn)確與否將會嚴(yán)重影響儲層評價(jià)的質(zhì)量。子波估算時(shí)，首先用地震數(shù)據(jù)體(振幅譜)估算子波，然后用常相位譜技術(shù)確定地震相位。最后，將井資料與地震資料相結(jié)合估算每口井子波的振幅和相位譜。在得到每口井的估算子波之后，借助井編輯和嵌入子波工具，從本工區(qū)13口井資料中，選取相關(guān)的井，并應(yīng)用多口井子波平均技術(shù)估算最終的平均子波。

本次反演分別從13口井資料中提取各自的子波，并用這些提取的子波計(jì)算平均子波。圖6是Sufyan三維區(qū)制作合成地震記錄所提取的統(tǒng)計(jì)子波及其振幅譜、相位譜與頻譜?？梢钥闯鏊袇⑴c反演的單井的子波形態(tài)、能量、相位、振幅及頻寬基本一致，相位穩(wěn)定，從而說明統(tǒng)計(jì)子波可靠合理。所得到的平均子波，形態(tài)近似rick子波，主瓣能量集中，旁瓣較小，相位穩(wěn)定，是比較理想的平均子波。

(5)低頻模型的建立。

由于地震采集系統(tǒng)的限制，地震直接反演結(jié)果中不包含10 Hz 以下的低頻成分，須從其他資料提取予以補(bǔ)償。從地震資料出發(fā)，以測井資料和鉆井?dāng)?shù)據(jù)為基礎(chǔ)，建立基本反映沉積體地質(zhì)特征的低頻初始模型。在地震反演中，初始地質(zhì)模型的合理建立是很重要的，特別是對模型反演來說，反演結(jié)果的好壞很大程度上由初始模型即先驗(yàn)地質(zhì)認(rèn)識決定，因此建立初始模型是做好基于模型反演的關(guān)鍵。

圖5 儲層精細(xì)標(biāo)定結(jié)果Fig.5 Results of detailed reservoir calibration

圖6 最終的子波Fig.6 Final results of average wavelet

建立盡可能接近實(shí)際地層條件的初始波阻抗模型，是減少其最終結(jié)果多解性的根本途徑。測井資料在縱向上詳細(xì)揭示了巖層的變化細(xì)節(jié)，地震資料則連續(xù)記錄了界面的橫向變化，二者的結(jié)合為精確建立空間波阻抗模型提供了必要的條件。

低頻模型的建立包括兩部分內(nèi)容：建立地質(zhì)模型及生成低頻模型，主要受地震解釋層位的控制。地質(zhì)模型的建立產(chǎn)生所有建模所需的參數(shù)(儲層框架和井比重等)。地質(zhì)模型建好后，基于地質(zhì)模型對井曲線進(jìn)行內(nèi)插，并在層位的控制下生成低頻模型，再用這個(gè)低頻模型進(jìn)行約束稀疏脈沖反演。

在建模過程中，地震解釋層位AG1、AG2 和 AG3用來確保內(nèi)部層序的一致性，采用了克里金插值法，該方法在井網(wǎng)分布不規(guī)則時(shí)，能有效反映儲層參數(shù)宏觀變化趨勢，精度相對高。圖7是內(nèi)插的低頻模型。

圖7 連井波阻抗初始模型剖面Fig.7 Initial cross section model of impedance

(6)砂體厚度預(yù)測。

根據(jù)前面的巖石物理分析可知，縱波阻抗可以較好區(qū)分AG組砂泥巖，同時(shí)根據(jù)交會結(jié)果分析可確定區(qū)分砂巖的門檻值為8000 g/cm3*m/s。在求取砂巖厚度的過程中，首先利用砂巖屬性門檻值提取砂巖時(shí)間厚度，然后根據(jù)地震速度初步計(jì)算出砂巖厚度，最后利用單井砂巖厚度值對預(yù)測砂巖厚度進(jìn)行趨勢校正，得到最終砂巖厚度預(yù)測等厚圖以及含砂率預(yù)測等厚圖[14](圖8、圖9)。

從Sufyan地區(qū)砂體厚度預(yù)測圖和含砂率預(yù)測圖可以看出，在AG3沉積時(shí)期，砂體主要分布在西部斜坡帶，最厚達(dá)到65 m(圖8d)，凹陷中部以泥質(zhì)沉積為主，含砂率分布與砂巖厚度分布相似，在西斜坡及Suf構(gòu)造帶含量高，向東南減少(圖9d)。在AG2下段沉積時(shí)期，整個(gè)凹陷內(nèi)部砂體厚度明顯變薄，主要分布在中央構(gòu)造帶等構(gòu)造較高的部位，砂體厚度小且分布范圍小(圖8c)，含砂率較低，約15%以下(圖9c)。在AG2上段沉積時(shí)期，砂體厚度較AG2下段時(shí)期明顯增厚且分布范圍變大(圖8b)。在AG1沉積時(shí)期，砂體向凹陷內(nèi)推進(jìn)，到達(dá)中

央構(gòu)造帶內(nèi)部(圖8a)，AG1段砂體發(fā)育，含砂率高，最高可達(dá)85%(圖9a)；由于埋藏相對較淺，分選較好，泥質(zhì)含量低，儲層物性好。

總體而言，Sufyan凹陷在AG期以半深湖、三角洲前緣為主，砂體含砂量和厚度在不同構(gòu)造區(qū)內(nèi)變化較大，且分布不均勻，儲層預(yù)測對油氣勘探部署非常重要。

2.3 反演結(jié)果分析

根據(jù)實(shí)際地質(zhì)認(rèn)識以及圖4的分析結(jié)果可以認(rèn)為，利用8000 g/cm3*m/s為波阻抗的門檻值來區(qū)分砂泥巖，砂巖在反演剖面上對應(yīng)高波阻抗同相軸，而泥巖則對應(yīng)低波阻抗同相軸，波阻抗反演成果能有效識別砂體展布特點(diǎn)。

通過合成記錄的標(biāo)定，將13口井資料與波阻抗剖面疊合。通過將地震反演結(jié)果與井上的波阻抗曲線進(jìn)行對比，其目的層段砂泥巖的分布與井上實(shí)際井曲線基本吻合，說明了地震反演結(jié)果的可靠性(圖10)，可認(rèn)為本次稀疏脈沖反演很好地揭示了Sufyan地區(qū)AG組儲集砂體內(nèi)幕和側(cè)向變化特征。

圖8 Sufyan地區(qū)主要目的層砂體厚度預(yù)測圖Fig.8 Main formations sand thickness prediction of Sufyan sag

圖9 Sufyan地區(qū)主要目的層含砂率預(yù)測圖Fig.9 Main formations sand ratios prediction of Sufyan sag

3 結(jié)論與建議

(1)約束稀疏脈沖反演技術(shù)應(yīng)用過程中，需要注意以下幾點(diǎn)：①測井曲線標(biāo)準(zhǔn)化處理能夠消除非地質(zhì)因素造成的井間一致性問題，明確儲層地球物理特征，使低頻模型準(zhǔn)確、儲層預(yù)測可靠。②井震標(biāo)定過程中，要考慮井震資料之間的速度差異，根據(jù)多個(gè)標(biāo)志層來確定時(shí)深漂移量，保障合理的時(shí)深關(guān)系、子波的準(zhǔn)確估算。

圖10 CSSI反演結(jié)果與井上波阻抗曲線對比圖Fig.10 Contrast of CSSI inversion results and impedance curve

(2)利用稀疏脈沖反演技術(shù)在Sufyan凹陷開展儲層預(yù)測，取得了較好的效果，主要表現(xiàn)在：①反演結(jié)果與已鉆井吻合程度高。②目的層段薄砂層橫向分布規(guī)律明確，可以指導(dǎo)工區(qū)沉積微相分析工作；結(jié)合鉆井，能為單砂體刻畫提供依據(jù)。

(3)反演人員對地震工區(qū)內(nèi)各種勘探資料以及前人研究成果的認(rèn)識和掌握程度也對反演結(jié)果有著重大影響，這體現(xiàn)在波阻抗約束條件的制定等反演過程的控制以及對地震資料所缺失的低頻分量的補(bǔ)充等多個(gè)方面。

[1] 李國發(fā),王艷倉,熊金良,等.地震波阻抗反演實(shí)驗(yàn)分析[J].石油地球物理勘探,2010,45(6):868-872.

[2] 齊金成,厲玉樂,韋學(xué)銳,等.約束稀疏脈沖波阻抗反演及其應(yīng)用[J].內(nèi)蒙古石油化工,2004,30(1):120-122.

[3] 何麗娟,熊冉,雷海飛,等.約束稀疏脈沖反演在大慶茂興地區(qū)泉四段儲層研究中的應(yīng)用[J].沉積與特提斯地質(zhì),2009,29(1):92-95.

[4] GENIK G J. Petroleum geology of Cretaceous_Tertiary rift basins in Niger, Chad, and Central African Republic[J].AAPG,1993,77(8):1405-1434.

[5] CHULL T J. Rift basins of interior Sudan: petroleum exploration and discovery[J]. AAPG Bulletin, 1988,72:1128-1142.

[6] 童曉光,竇立榮,田作基,等,蘇丹穆格萊特盆地的地質(zhì)模式和成藏模式[J].石油學(xué)報(bào),2004,25(1):19-24．

[7] 竇立榮,潘校華,田作基,等.蘇丹裂谷盆地油氣藏的形成與分布——兼與中國東部裂谷盆地對比分析[J].石油勘探與開發(fā),2006,33(3):255-261.

[8] 張亞敏,陳發(fā)景.穆格萊德盆地形成特點(diǎn)與勘探潛力[J].石油與天然氣地質(zhì),2002,23(3):236-240.

[9] 張亞敏.蘇丹國穆格萊德盆地構(gòu)造特征及演化[J].西安石油大學(xué)學(xué)報(bào),2008,23(3):38-42.

[10] 于秀英,古正富,賈俊杰,等.儲層孔隙結(jié)構(gòu)測井表征的新方法[J].非常規(guī)油氣,2016,3(1):14-20.

[11] 邵才瑞,宋慧瑩,張福明,等.利用標(biāo)準(zhǔn)砂層組采用趨勢面分析進(jìn)行多井資料標(biāo)準(zhǔn)化[J].地球物理學(xué)進(jìn)展,2011,26(3):1031-1038.

[12] 徐德龍,李濤,黃寶華,等.利用交匯圖法識別國外M油田巖性與流體類型的研究[J].地球物理學(xué)進(jìn)展,2012,27(3):1123-1132.

[13] 陸基孟.地震勘探原理[M].北京:石油大學(xué)出版社,1983:50-60.

[14] 許翠霞,馬朋善,賈志坤,等.致密火山碎屑巖含氣敏感屬性研究[J].非常規(guī)油氣,2016,3(5):8-13.

TheApplicationofConstrainedSparseSpikeInversioninSufyanSag

Yang Cang1, Zou Quan2, Ke Weili2

(1.CNPCInternational(Nile)Ltd.,Beijing100034,China; 2.PetroChinaResearchInstituteofPetroleumExploration&Development,Beijing100083,China)

Due to the sand bodies are thinner and the lateral reservoir change seriously of deep AG member in Sufyan sag, the reservoir prediction is very difficult to do. However, the constrained sparse spike inversion is widely used for reservoir prediction in many oil fields of the world. This method has high precision and has less unknown various, and can preferably overcome the multi-answers problem. It is able to visually depict the characteristics of the spatial distribution of lithology. We used the constrained sparse spike inversion technique, through a series of work, such as the analysis of seismic data, well logging curve standardized processing, fine well-to-seismic calibration and wavelet extraction and so on, optimized inversion parameters, deep AG member in Sufyan sag has carried on the reservoir prediction, provided a strong support for the reservoir evaluation in this area. Application showed that the forecast effects after technical processing are consistent with the drilling data and the sedimentary characteristics in this area, it can provide the basis for future development.

Sufyan sag; AG formation; constrained sparse spike inversion; reservoir prediction

楊倉(1983—)，男，北京人，漢族，在職碩士研究生，中國石油尼羅河公司工程師，主要從事物探工作。郵箱：kamangta@126.com.

P31

:B