張 勇，鐘 薇

（1.中國石化華東分公司物探研究院，江蘇 南京 210007；2.中國石化華東分公司非常規(guī)資源處，江蘇 南京 210007）

隨著油氣田勘探和開發(fā)的不斷深入，勘探目標由構(gòu)造油氣藏向巖性油氣藏和隱蔽油氣藏轉(zhuǎn)變，儲層預測技術(shù)逐漸成為巖性油氣藏勘探中不可或缺的重要技術(shù)。地震反演是儲層預測的核心[1]，主要包括疊后反演和疊前反演技術(shù)。疊后反演以褶積模型為基礎，假設地下的強反射系數(shù)界面（不整合面或主要的巖性分界面）是稀疏的、不連續(xù)分布的，目的是在最少的反射系數(shù)條件下達到合成記錄與地震記錄的最小誤差，從而得到縱波阻抗數(shù)據(jù)[2]。疊前反演[3-7]能反映地震反射振幅隨偏移距不同（AV0）或入射角不同而變化（AVA）的更多的特點，能得到縱波阻抗、橫波阻抗、泊松比等彈性參數(shù)，但受限于疊前道集處理以及計算量等原因，現(xiàn)階段勘探開發(fā)工作中主要以疊后反演為主。疊后反演中，以約束稀疏脈沖反演為代表的反演只能反映大套砂體展布情況，對較薄的砂巖層及砂泥巖互層識別能力不夠，難以滿足非構(gòu)造油氣藏勘探開發(fā)需要。而建立在約束稀疏脈沖反演基礎上的疊后地質(zhì)統(tǒng)計學反演[8-11]，引入了隨機建模的思想，兼顧了地震數(shù)據(jù)的橫向分辨率高、地質(zhì)信息豐富的特點，將測井數(shù)據(jù)在單點垂向分辨率上的優(yōu)勢發(fā)揮出來，能更加準確預測單砂體變化規(guī)律。

研究區(qū)位于海安凹陷曲塘深凹北斜坡帶，面積82 km2，全三維覆蓋，擁有6口探井完井資料。目的層為下第三系阜寧組三段，物源主要來自西北，為三角洲前緣亞相，阜寧組三段分為上砂組和下砂組，中間一套較穩(wěn)定泥巖層。上砂組為油氣順斷層運移而在高部位物性較好、蓋層保存條件好的區(qū)域成藏，以構(gòu)造控油為主。上砂組巖性為灰白色、褐灰色粉—細砂巖夾灰綠色、灰黑色泥巖，砂巖儲層厚度較薄，層數(shù)少，砂體變化快,單砂體厚度3～5 m。下砂組為垂向近源運移到物性較好的區(qū)域成藏，主要為物性控藏。下砂組巖性同為灰白色、褐灰色粉—細砂巖夾灰綠色、灰黑色泥巖。砂巖儲層較厚，層數(shù)較多、砂體連通情況復雜，單砂體厚度4～8 m。試油結(jié)果顯示阜寧組三段油層具有較好的勘探前景，為了深入挖掘研究區(qū)潛力，需對阜寧組三段油層進行測井小層對比、精細構(gòu)造解釋以加強對阜寧組三段油層的認識。在阜寧組三段進行約束稀疏脈沖反演基礎上開展地質(zhì)統(tǒng)計學反演，研究砂體展布特征，預測有利成藏區(qū)帶。

1 方法原理

地質(zhì)統(tǒng)計學反演是將隨機模擬理論與地震反演相結(jié)合的反演方法，其過程主要由隨機模擬和在對研究區(qū)地質(zhì)認識的基礎上，模擬結(jié)果進行優(yōu)選兩部分組成。地質(zhì)統(tǒng)計學反演認可輸入信息（測量數(shù)據(jù)、地質(zhì)認識等）都包含不確定性，得到的多個實現(xiàn)都滿足所有輸入信息，能反映多個數(shù)據(jù)的不確定性。反演主要算法包括：統(tǒng)計學建模、貝葉斯推斷、馬爾科夫鏈—蒙特卡羅[12-13]算法。具體反演過程中通過概率密度函數(shù)（PDF）和變差函數(shù)對離散屬性（巖性、孔隙度）和彈性屬性（波阻抗）的分析。概率密度函數(shù)主要描述某一屬性在空間的概率分布情況，可以通過高斯、對數(shù)高斯、Binned等子區(qū)間、常數(shù)分布等一個和多個組合類型來描述。離散屬性的概率密度函數(shù)一般從巖相曲線統(tǒng)計的百分比或波阻抗反演的巖性概率體中統(tǒng)計獲得。彈性屬性的概率密度函數(shù)可以從測井曲線數(shù)據(jù)擬合獲得。變差函數(shù)主要描述某屬性的空間展布特征隨距離的變化，是距離的一個函數(shù)。在求取變差函數(shù)過程中，首先是分析縱橫向?qū)嶒炞儾詈瘮?shù)，縱向?qū)嶒炞儾詈瘮?shù)是根據(jù)測井曲線統(tǒng)計獲得，橫向?qū)嶒炞儾詈瘮?shù)從稀疏脈沖波阻抗反演所提取的平面屬性獲得，實驗變差函數(shù)學表達式為：

在求取實驗變差函數(shù)后，并不能直接用于屬性模擬，實驗變差函數(shù)不能保證被預測點的變差是正值，以及實驗變差函數(shù)中按某些滯后距可能搜索不到數(shù)據(jù)，算不出對應的實驗變差函數(shù)。但模擬需要對每個整數(shù)倍滯后距計算對應的變差函數(shù)，不允許缺失數(shù)據(jù)。模型變差函數(shù)能滿足這兩點要求，在反演中都只能使用實驗變差函數(shù)擬合后的模型變差函數(shù)參與計算。模型變差函數(shù)主要分為指數(shù)型和高斯型，或兩者按照一定比例分配權(quán)重，指數(shù)型變差函數(shù)兩點間的相關性快速降低，常用于模擬變化快的數(shù)據(jù)，比如彈性屬性模擬結(jié)果變化很快，看起來不平滑；高斯型變差函數(shù)兩點間的相關性降低很慢，常用于模擬非常連續(xù)的數(shù)據(jù)，比如離散屬性，模擬結(jié)果變化很慢，看起來很平滑。模型變差函數(shù)數(shù)學表達式為：

在分析縱橫向變程時，要根據(jù)不同研究區(qū)情況選擇不同縱橫向變程。

2 應用效果分析

在地質(zhì)統(tǒng)計學反演之前需要完成一個高質(zhì)量的稀疏脈沖波阻抗反演，高質(zhì)量的稀疏脈沖波阻抗反演需要在以下各環(huán)節(jié)都要高質(zhì)量完成，包括正確的井震標定、合適的子波提取、符合地質(zhì)認識的低頻模型以及合適的質(zhì)控參數(shù)。而正確的井震標定（圖1）是完成后續(xù)工作的基礎，它將時間域的地震資料與深度域的測井資料結(jié)合起來，將直接影響稀疏脈沖波阻抗反演結(jié)果。在完成了稀疏脈沖波阻抗反演的基礎上，才能繼續(xù)開展地質(zhì)統(tǒng)計學反演。地質(zhì)統(tǒng)計學反演時，利用了研究區(qū)內(nèi)5口探井參與反演，1口探井（well3）作為反演結(jié)果預測性的驗證井。

2.1 剖面識別效果

地質(zhì)統(tǒng)計學反演結(jié)果是否正確，需要與稀疏脈沖波阻抗反演結(jié)果進行對比，從過well1井波阻抗反演結(jié)果（圖2）和地質(zhì)統(tǒng)計學反演結(jié)果（圖3）上可以看出，研究區(qū)高阻抗儲層在稀疏脈沖波阻抗反演和地質(zhì)統(tǒng)計學反演結(jié)果中展布情況較吻合，上下兩套高阻抗層，中間一套較穩(wěn)定的低波阻抗層。地質(zhì)統(tǒng)計學反演波阻抗較稀疏脈沖波阻抗反演分辨率明顯提高，薄層儲細節(jié)特征清楚。well3井實際鉆探結(jié)果顯示，在阜寧組三段含有砂巖儲層十三層，從過well3井反演結(jié)果巖性剖面（圖4）看，在well3井沒有作為約束井的情況下，地質(zhì)統(tǒng)計學反演巖性剖面能與well3井巖性解釋結(jié)果較好吻合，地質(zhì)統(tǒng)計學反演結(jié)果具有較好的預測性。

2.2 平面識別效果

研究區(qū)處于三角洲前緣亞相，砂體展布應該呈現(xiàn)出從西北部到東南部發(fā)育逐漸減少。從上砂組砂體平面展布特征（圖5）看，上砂組砂體在西北部較發(fā)育，東南部砂體發(fā)育程度差，與單井揭示砂體發(fā)育吻合；下砂組砂體（圖6）在整個研究區(qū)內(nèi)砂體較發(fā)育，有從西北部往東南部逐漸減少的趨勢，所反映的砂體展布規(guī)律與該地區(qū)三角洲前緣亞相沉積環(huán)境及井數(shù)據(jù)符合。

圖1 well4井井震標定Fig.1 Calibrating downhole seismic of well4

圖2 過well1井稀疏脈沖反演波阻抗剖面圖Fig.2 Rarefaction pulse inversion wave impedance profile through well1

圖3 過well1井地質(zhì)統(tǒng)計學反演波阻抗剖面圖Fig.3 Geostatistical inversion wave impedance profile through well1

圖4 過well3井地質(zhì)統(tǒng)計學反演巖性概率剖面圖Fig.4 Geostatistical inversion lithology probability profile through well3

圖5 阜寧組三段上砂組砂體平面展布圖Fig.5 Sandbody plane distribution of the upper reservoir of the 3rdmember in Funing group

圖6 阜寧組三段下砂組砂體平面展布圖Fig.6 Sandbody plane distribution of the lower reservoir of the 3rdmember in Funing group

3 結(jié)論和認識

以海安凹陷曲塘深凹為例，綜合利用地質(zhì)、測井、地震資料，在稀疏脈沖波阻抗反演基礎上，應用地質(zhì)統(tǒng)計學反演進行儲層預測。充分利用了測井資料垂向分辨率和地震數(shù)據(jù)橫向分辨率的優(yōu)勢，較好反映儲層砂體縱橫向變化，取得了良好的預測效果。預測結(jié)果顯示：研究區(qū)阜寧組三段上砂組和下砂組砂體在西北部砂體較發(fā)育，東南部砂體欠發(fā)育的特征，與三角洲前緣沉積特征一致，符合研究區(qū)的地質(zhì)認識以及實際鉆探情況，表明地質(zhì)統(tǒng)計學反演能有效預測薄儲層發(fā)育，能清楚反映單砂體在井間的變化，為研究區(qū)后續(xù)的勘探和開發(fā)提供有效依據(jù)。

[1]王雅春，王璐.地質(zhì)統(tǒng)計學反演在杏北西斜坡區(qū)儲層預測中的應用[J].地球物理學進展，2013，28（5）：2554-2560.

[2]原福堂，董紅，朱再文，等.地質(zhì)統(tǒng)計隨機反演與稀疏脈沖反演在儲層預測中的應用對比[J].石油天然氣學報，2014，36（7）：55-59.

[3]姜巖，徐立恒，張秀麗，等.疊前地質(zhì)統(tǒng)計學反演方法在長垣油田儲層預測中的應用[J].地球物理學進展，2013，28（5）：2579-2586.

[4]郎曉玲，彭仕宓，康洪全，等.疊前同時反演方法在流體識別中的應用[J].石油物探，2010，49（2）：164-169.

[5]李愛山，印興耀，張繁昌，等.疊前AVA多參數(shù)同步反演技術(shù)在含氣儲層預測中的應用[J].石油物探，2007，46（1）：64-68.

[6]崔成軍，龔姚進，申大媛.波阻抗反演在儲層預測研究中的應用[J].地球物理學進展，2010，25（1）：9-15.

[7]安鴻偉，李正文，李仁甫，等.稀疏脈沖波阻抗反演在YX油田開發(fā)中的應用[J].石油物探，2002，41（1）：56-60.

[8]潘昱潔，李大衛(wèi)，楊鍇.確定性反演和隨機反演對井約束條件的需求分析[J].石油物探，2011，50（4）：345-360.

[9]劉百紅，李建華，魏小東，等.隨機反演在儲層預測中的應用[J].地球物理學進展，2009，24（2）：581-589.

[10]王香文，劉紅，滕彬彬，等.地質(zhì)統(tǒng)計學反演技術(shù)在薄儲層預測中的應用[J].石油與天然氣地質(zhì)，2012，33（5）：730-735.

[11]葉云飛，劉春成，劉志斌，等：地質(zhì)統(tǒng)計學反演技術(shù)研究與應用[J].物探化探計算技術(shù)，2014，36（4）：446-451.

[12]趙林.馬爾可夫鏈蒙特卡羅模擬在儲層反演中的應用[J].石油天然氣學報，2009，31（2）：249-252.

[13]張小艷，譚勇.變差函數(shù)球狀模型的自動擬合與實現(xiàn)[J].物探與化探，2010，34（2）：253-257.