易 爽 夏潤(rùn)陶

（長(zhǎng)江大學(xué)，湖北 武漢 434023）

1 研究工區(qū)和數(shù)據(jù)集

1.1 工區(qū)概況和沉積環(huán)境

隨著地質(zhì)統(tǒng)計(jì)方法的發(fā)展，在80年代末期發(fā)展出指示克里金方法（IK），隨后A.Journel提出序貫指示模擬［1-4］。

應(yīng)用指示克里金（以下簡(jiǎn)稱IK）結(jié)合實(shí)際數(shù)據(jù)，分析環(huán)江油田環(huán)縣地區(qū)長(zhǎng)63的砂巖相展布特征，驗(yàn)證重力流沉積特征，對(duì)沉積相形態(tài)劃分做指導(dǎo)。環(huán)縣地區(qū)位于鄂爾多斯盆地西南部，長(zhǎng)6組3段處于湖盆衰退時(shí)期，由于構(gòu)造和地形因素，長(zhǎng)6、長(zhǎng)7深湖區(qū)域發(fā)育重力流沉積砂巖屬于三角洲前緣滑塌體沉積，實(shí)踐表明砂質(zhì)碎屑流沉積一般形成“下生上儲(chǔ)”或“上生下儲(chǔ)”的巖性油氣藏［5-8］，構(gòu)成了深水區(qū)域低孔低滲的儲(chǔ)集體。

研究工區(qū)處于環(huán)江油田環(huán)縣地區(qū)的東南部分，同時(shí)是眾多油井開(kāi)采區(qū)域。環(huán)縣地區(qū)位于鄂爾多斯盆地西緣，西邊靠近西緣沖斷帶——鄂爾多斯盆地一級(jí)構(gòu)造帶，處于天環(huán)坳陷帶內(nèi)，晚三疊世的構(gòu)造活動(dòng)造成盆地西陡東緩，東北高、西南低的古地理特點(diǎn)，沉降活動(dòng)主要發(fā)育在環(huán)縣以南［8-9］，晚三疊世長(zhǎng)6組處于湖盆衰退時(shí)期，由于坡折帶的存在，三角洲前緣沉積物滑塌在深湖泥之上，形成砂質(zhì)碎屑流舌狀體。這類砂巖底部具有印模構(gòu)造，發(fā)育在粉砂巖、細(xì)砂巖與泥巖接觸面之上，并且大多都是直接與泥巖接觸，構(gòu)成巖性突變。從測(cè)井曲線特征看，缺少河流沉積的旋回性，砂巖體成段出現(xiàn)，砂巖內(nèi)部為逆粒序——正粒序的變化規(guī)律?？梢?jiàn)，重力流沉積旋回中部密度較粗［6］，和頂、底兩端泥巖低孔低滲構(gòu)成一定圈閉性，形成巖性儲(chǔ)層。

本次研究將長(zhǎng)63層劃分為三個(gè)小層，每小層厚度在15m左右，主要參考了高振中［7］的沉積研究成果，按照重力流砂體內(nèi)部旋回性劃分。本次主力研究層段是長(zhǎng)632對(duì)應(yīng)砂體。

1.2 數(shù)據(jù)集

對(duì)應(yīng)工區(qū)內(nèi)分層結(jié)果，找到對(duì)應(yīng)長(zhǎng)632深度段以及對(duì)應(yīng)一套砂體，從取芯數(shù)據(jù)中篩選出24口井的孔隙度數(shù)據(jù)，對(duì)其求均值，代表各井對(duì)應(yīng)的孔隙度值，所有井的孔隙度值在9.35%和11.1%之間，屬于中低孔隙度，為含泥細(xì)砂巖至細(xì)砂巖。可以推測(cè)低值代表泥質(zhì)含量多的砂巖，高值代表較純凈的細(xì)砂巖。就沉積環(huán)境而言，低值屬于砂質(zhì)碎屑流邊緣或深湖泥區(qū)域，而高值可能對(duì)應(yīng)砂質(zhì)碎屑流中心區(qū)域，也是有利油氣儲(chǔ)集區(qū)域。

2 指示克里金

在未取芯地區(qū)，經(jīng)常要將孔隙度取一系列閾值（thresholds）來(lái)識(shí)別巖相。如果一種巖相能通過(guò)一個(gè)閾值來(lái)識(shí)別，則這種巖相在位置uα存在用1表示，不存在用0表示：

其中zk是進(jìn)行指示變換的一系列閾值；

這樣，原始數(shù)據(jù)轉(zhuǎn)變成了k類巖相新數(shù)據(jù)，每一類巖相只用0和1表示。這樣做的目的就是為了得到井上的數(shù)據(jù)小于或大于閾值的概率分布。區(qū)域A內(nèi)孔隙度小于一個(gè)閾值的概率為：

相對(duì)的，可以得到該處孔隙度大于閾值的概率：

如果樣本數(shù)據(jù)足夠多，即n足夠大，則有：

這里wi（i=1，2，…n）代表克里金體系里殘差值計(jì)算的權(quán)重值，對(duì)應(yīng)各個(gè)測(cè)量點(diǎn)。指示變換是直接對(duì)未測(cè)量區(qū)域不確定性分布的估計(jì)。輸入數(shù)據(jù)集的全局概率分布也同樣通過(guò)各個(gè)閾值來(lái)估計(jì)。

3 完善數(shù)據(jù)集

原有數(shù)據(jù)集來(lái)自24口井的巖心取樣孔隙度，對(duì)應(yīng)長(zhǎng)632的一套砂巖體，由于采樣密度較大，每口井做了求平均處理。由于缺少數(shù)據(jù)，從井位圖上可見(jiàn)，在東北和西南分布比較稀疏，因此，新添加6口“模擬井”到工區(qū)內(nèi)。6口井井位數(shù)據(jù)是工區(qū)內(nèi)已有井（但是缺少孔隙度值），最終孔隙度數(shù)據(jù)根據(jù)線性克里金插值模型獲得。由圖1可以看見(jiàn)具體井位和孔隙度值。選擇已存在井也是為了方便今后得到數(shù)據(jù)可以對(duì)孔隙度分布結(jié)果進(jìn)行驗(yàn)證。最終得到的30口井?dāng)?shù)據(jù)。

本文前提假設(shè)是：不同孔隙度值代表砂體泥質(zhì)含量不一樣，也代表不同的沉積環(huán)境和位置。理論上講，閾值劃分在5～10個(gè)是比較合適的，本文劃分6個(gè)閾值：9.3、9.73、10.02、10.2、10.32和11.08%。然后針對(duì)每一個(gè)閾值，通過(guò)SGeMS軟件計(jì)算其變差函數(shù)。

4 IK模擬結(jié)果

對(duì)于連續(xù)型變量，指示克里金算法為各個(gè)門檻值提供了獨(dú)立離散的模型，指示克里金允許不同概率值有不同的變程，相類之間區(qū)域的插值通過(guò)指定變差函數(shù)模型（如本次使用Exponential指數(shù)模型）完成過(guò)度。這樣模型具有很強(qiáng)的實(shí)用性，可以綜合軟數(shù)據(jù)。在網(wǎng)格設(shè)置里，X和Y均設(shè)置200，模型共有40 000個(gè)網(wǎng)格。圖2顯示了井位圖，顏色代表含不同孔隙度值。紅線是概率密度值，藍(lán)線是孔隙度概率累積分布。

圖1 工區(qū)井位

選定好6個(gè)閾值后，每個(gè)閾值做一個(gè)概率插值圖，這樣可以觀察孔隙度在相鄰兩個(gè)閾值之間的砂巖分布狀況變化，而且通過(guò)比較不同閾值的概率圖，可以發(fā)現(xiàn)有些砂體在不同閾值范圍內(nèi)消失，而有些高值部分砂體顯現(xiàn)更加清晰。從砂體展布形態(tài)來(lái)看，主要是高值部分，更加獨(dú)立，沒(méi)有很好的繼承性關(guān)系，孔隙度10.3%對(duì)應(yīng)的概率圖可見(jiàn)，孔隙度最高值分布于工區(qū)南北兩端，而中間部分傳統(tǒng)的人為劃分沉積相圖是很難將砂體具體形態(tài)展布顯示出來(lái)，通過(guò)IK模擬，工區(qū)內(nèi)不同特征的砂體細(xì)節(jié)關(guān)和形態(tài)系可以很好展現(xiàn)，從而可以得到較為細(xì)致的沉積相分布形狀和范圍。

圖2 對(duì)應(yīng)6個(gè)閾值的孔隙度概率圖，概率1（紅色）表示該處孔隙度值確定小于閾值孔隙度，概率0（藍(lán)色）表示該處孔隙度值大于閾值孔隙度。

5 結(jié)論

IK十分適合展示平面的不確定性，與其他克里金方法是線性的不同，IK是非線性的差值方法［10］。應(yīng)用IK得到了不同的截?cái)嘀岛蛯?duì)應(yīng)孔隙度分布，可以據(jù)此很好地揭示工區(qū)內(nèi)長(zhǎng)632砂體展布特征，即：砂質(zhì)碎屑流延伸的前端發(fā)育于三角洲前緣坡折帶以下的湖盆區(qū)域，沿物源方向（本次主要是西北）上延伸不連續(xù)，但垂直物源方向疊置連片［11-13］，是湖盆斜坡與湖盆中心有利的目標(biāo)區(qū)域。盡管砂質(zhì)碎屑流引入湖盆相解釋深水重力流理論還比較新，但是新的理論可以幫助找到新的有利建產(chǎn)區(qū)，借助IK可以圈定有利區(qū)的大致理論范圍，對(duì)后續(xù)的勘探部署有指導(dǎo)意義。

［1］傅強(qiáng)，呂苗苗，劉永斗.鄂爾多斯盆地晚三疊世湖盆濁積巖發(fā)育特征及地質(zhì)意義［J］.沉積學(xué)報(bào)，2008，26（2）：186-191.

［2］陳飛，胡光義，等.鄂爾多斯盆地富縣地區(qū)上三疊統(tǒng)延長(zhǎng)組砂質(zhì)碎屑流沉積特征及其油氣勘探意義［J］.沉積學(xué)報(bào)，2012，30（6）：1042-1052.

［3］Journel， A.G..Non-parametric estimation of spatial distribution.Mathematical Geology，1983，15（3）：445-468.

［4］Journel， A.G..Fundamentals of Geostatistics in Five Lessons［J］.ShortCourse in Geology.V.8.AGU.1989.

［5］楊俊杰.鄂爾多斯盆地構(gòu)造演化與油氣分布規(guī)律［M］.北京：石油工業(yè)出版社，2002：50-107.

［6］高振中.鄂爾多斯盆地環(huán)縣地區(qū)上三疊統(tǒng)延長(zhǎng)組長(zhǎng)63砂層組沉積微相特征及新認(rèn)識(shí)［J］.古地理學(xué)報(bào)，2012（6）：37-40.

［7］Sen，Z..Spatial Modeling Principles in Earth Sciences.Springer-Verlag，2009.

［8］鄒才能，趙政璋，楊華，等.陸相湖盆深水砂質(zhì)碎屑流成因機(jī)制與分布特征—以鄂爾多斯盆地為例［J］.沉積學(xué)報(bào)，2009，27（6）：1065—1075.

［9］吳崇筠，李純菊，劉國(guó)華，等.斷陷湖盆中的濁積巖［C］.北京：石油工業(yè)出版社，2011.

［10］趙俊興，李鳳杰，申曉莉，等.鄂爾多斯盆地南部長(zhǎng)6長(zhǎng)7油層濁流事件的沉積特征及發(fā)育模式［J］.石油學(xué)報(bào)，2008，29（3）：186-191.

［11］劉昭君.湖泊水下扇沉積特征及影響因素—以伊通盆地莫里青斷陷雙陽(yáng)組為例［J］.沉積學(xué)報(bào)，2003，21（1）：148-154.

［12］張金亮，沈鳳.烏爾遜凹陷大磨拐河組近岸水下扇儲(chǔ)層特征［J］，石油學(xué)報(bào)，1991，12（3）：25-35.

［13］鄒才能，陶士振，袁選俊，等.連續(xù)型油氣藏形成條件與分布特征.［J］，石油學(xué)報(bào)，2009，30（3）：324-331.