郭 奇 莊天琳 何書梅 李 禎 魏朋朋 劉麗杰

(①中國(guó)石化勝利油田分公司勘探開發(fā)研究院，山東東營(yíng) 257015； ②中國(guó)石化勝利石油管理局有限公司博士后科研工作站，山東東營(yíng) 257000； ③中國(guó)石油大港油田公司勘探開發(fā)研究院，天津 300280)

0 引言

油藏進(jìn)入高含水開發(fā)期后，對(duì)儲(chǔ)層物性的精細(xì)識(shí)別是剩余油挖潛的關(guān)鍵。現(xiàn)場(chǎng)通常由測(cè)井資料解釋得到沿井的滲透率模型，然而不同學(xué)者的研究[1-5]證明，受注水沖刷影響，儲(chǔ)層滲透率在開發(fā)后期會(huì)發(fā)生很大變化，在高度水淹層投產(chǎn)的新井，其測(cè)井解釋結(jié)果已不能反映油藏開發(fā)初期的物性情況。

針對(duì)注水沖刷造成的儲(chǔ)層滲透率時(shí)變問題，很多學(xué)者從不同角度進(jìn)行了研究。孫國(guó)[6]、徐守余等[7]通過巖心測(cè)井資料研究了儲(chǔ)層宏觀參數(shù)變化規(guī)律，并利用神經(jīng)網(wǎng)絡(luò)方法建立儲(chǔ)層參數(shù)的四維地質(zhì)模型。彭仕宓等[8]通過建立注水開發(fā)前、后儲(chǔ)層參數(shù)數(shù)學(xué)演化模型，綜合巖心、測(cè)井和實(shí)驗(yàn)測(cè)試資料對(duì)注水過程中儲(chǔ)層參數(shù)的變化規(guī)律進(jìn)行了研究。束青林等[9]、李紅南等[10]充分利用巖心資料，通過神經(jīng)網(wǎng)絡(luò)方法建立了宏觀和微觀儲(chǔ)層參數(shù)的動(dòng)態(tài)模型，并對(duì)剩余油分布模式做出了評(píng)價(jià)。張楓等[11]利用不同開發(fā)時(shí)期的測(cè)井資料，建立了儲(chǔ)層參數(shù)隨開發(fā)時(shí)間變化的數(shù)學(xué)模型，形成了不同開發(fā)時(shí)間四維動(dòng)態(tài)演化模型，并對(duì)研究區(qū)剩余油分布進(jìn)行了預(yù)測(cè)。么忠文等[12]、郭敬民[13]利用不同含水時(shí)期的水淹測(cè)井資料建立了水淹四維地質(zhì)模型，以表征注水開發(fā)過程中儲(chǔ)層流體動(dòng)態(tài)變化。上述研究大多是通過實(shí)驗(yàn)或理論模型推導(dǎo)得到儲(chǔ)層參數(shù)時(shí)變情況，然而地下流—固耦合的變化是復(fù)雜的，單一巖心實(shí)驗(yàn)或理論模型并不能夠真實(shí)反映儲(chǔ)層滲透率的變化情況。

地震屬性反映地震波的幾何特征與物理特征，這些特征在不同程度上反映了儲(chǔ)層的巖性、滲透率以及孔隙流體的性質(zhì)。由于對(duì)采集、處理因素所導(dǎo)致的時(shí)間、振幅、頻率等方面的差異進(jìn)行了處理，因此時(shí)移地震資料地震響應(yīng)差異能夠連續(xù)反映開發(fā)過程中的動(dòng)態(tài)變化。本文充分利用三期時(shí)移地震資料預(yù)測(cè)研究區(qū)儲(chǔ)層參數(shù)，在降低非油氣藏因素影響的同時(shí)，提高油氣藏內(nèi)部滲透率參數(shù)預(yù)測(cè)精度，能夠更好地突顯油氣藏儲(chǔ)層滲透率參數(shù)在長(zhǎng)期注水開發(fā)過程中的變化情況。

1 研究區(qū)概況

研究區(qū)位于黃驊凹陷中區(qū)北大港二級(jí)斷裂構(gòu)造帶的東南部，是油田的主要開發(fā)區(qū)之一。其北以港東主斷層為界，南鄰歧口凹陷，西至聯(lián)盟地區(qū)，東到灘海區(qū)的港深64—港63-2井區(qū)，面積約50km2。研究區(qū)是河流相沉積的復(fù)合式砂巖油藏，在開發(fā)中油藏出現(xiàn)剩余油高度分散、注水沖刷儲(chǔ)層滲透率變化等問題。

目前研究區(qū)內(nèi)有1993年、2003年、2015年采集的三期全區(qū)覆蓋地震數(shù)據(jù)， 分別對(duì)應(yīng)中低、高和特高含水三個(gè)開發(fā)階段。針對(duì)資料因采集、處理因素導(dǎo)致的時(shí)間、振幅、頻率等方面的差異進(jìn)行了處理，提高了三期地震數(shù)據(jù)各項(xiàng)參數(shù)的一致性。通過偏最小二乘回歸方法對(duì)多期地震屬性進(jìn)行分析，評(píng)價(jià)建立模型的擬合度，預(yù)測(cè)得到三個(gè)含水階段的滲透率數(shù)據(jù)場(chǎng)，并通過示蹤劑資料對(duì)預(yù)測(cè)結(jié)果進(jìn)行驗(yàn)證，形成了一套儲(chǔ)層參數(shù)時(shí)變性研究計(jì)算方法。

2 時(shí)移地震屬性優(yōu)選

地震屬性優(yōu)化是地震屬性分析技術(shù)的關(guān)鍵步驟，屬性優(yōu)化的質(zhì)量直接影響預(yù)測(cè)結(jié)果的準(zhǔn)確度和可靠性。地震屬性攜帶有大量?jī)?chǔ)層信息，并且與預(yù)測(cè)對(duì)象之間的關(guān)系復(fù)雜，預(yù)測(cè)對(duì)象對(duì)不同類型的地震屬性具有不同的敏感度，因此需要優(yōu)選出對(duì)預(yù)測(cè)對(duì)象最為敏感的地震屬性。

優(yōu)選地震屬性的原則一般包括：①針對(duì)具體研究目標(biāo)，選擇對(duì)其相關(guān)性最高、最敏感、能夠反映其本質(zhì)特征的地震屬性或?qū)傩越M合； ②選擇相互獨(dú)立、盡可能低維的屬性集，減少冗余信息，達(dá)到最優(yōu)屬性集結(jié)構(gòu)； ③優(yōu)選屬性個(gè)數(shù)，使有效信息最大化，剔除干擾屬性，減少因?qū)傩苑N類增加帶來的不利影響。

通過文獻(xiàn)調(diào)研各類地震屬性物理意義發(fā)現(xiàn)，振幅類和頻率類地震屬性被較多應(yīng)用在預(yù)測(cè)儲(chǔ)層滲透率參數(shù)響應(yīng)特征。在此基礎(chǔ)上，對(duì)目的層段不受斷層影響且地震屬性數(shù)值正常的214口井?dāng)?shù)據(jù)提取了10種地震屬性，并對(duì)沿井地震屬性值與實(shí)際滲透率數(shù)值進(jìn)行了相關(guān)性分析。由表1統(tǒng)計(jì)數(shù)據(jù)可見，其中某些地震屬性與滲透率的相關(guān)性較差。單一屬性與滲透率的相關(guān)系數(shù)均小于0.50，說明單一屬性的預(yù)測(cè)效果相對(duì)較差，需要通過屬性組合更好地反映儲(chǔ)層參數(shù)。本次研究從10個(gè)屬性中優(yōu)選出5個(gè)相關(guān)性較強(qiáng)的地震屬性，分別是均方根振幅、反射強(qiáng)度、瞬時(shí)頻率、振幅包絡(luò)和主頻，為下一步開展時(shí)移地震多屬性預(yù)測(cè)滲透率變化奠定基礎(chǔ)。

表1 單一地震屬性與滲透率相關(guān)性統(tǒng)計(jì)表

3 不同含水階段時(shí)移地震多屬性預(yù)測(cè)方法研究

多屬性預(yù)測(cè)分析是將優(yōu)選出的敏感地震屬性集或經(jīng)降維處理的地震屬性集作為輸入，以預(yù)測(cè)測(cè)井?dāng)?shù)據(jù)為輸出，利用線性、非線性方法，將地震屬性與已知井的巖性、物性、含油氣性等儲(chǔ)層信息相結(jié)合，對(duì)儲(chǔ)層參數(shù)進(jìn)行定量預(yù)測(cè)[14]； 同時(shí)綜合地質(zhì)、測(cè)井、和實(shí)際生產(chǎn)等資料對(duì)預(yù)測(cè)結(jié)果的可靠性進(jìn)行評(píng)價(jià)和分析。目前常用的多屬性預(yù)測(cè)方法有多元逐步回歸方法[15-16]、神經(jīng)網(wǎng)絡(luò)方法[17-18]、協(xié)克里金方法、支持向量機(jī)方法[19]、非參數(shù)回歸分析方法等[20-22]。本次研究選用偏最小二乘回歸方法建立不同期次時(shí)移地震屬性與測(cè)井滲透率的相關(guān)性。偏最小二乘回歸算法集成了主成分分析、典型相關(guān)分析、非線性回歸分析的優(yōu)點(diǎn)，通過建立多自變量與因變量間的回歸關(guān)系對(duì)模型進(jìn)行預(yù)測(cè)。相較于普通的多元回歸方法，該方法可以更好地處理各自變量之間存在的多重相關(guān)性以及數(shù)據(jù)量小、樣本較少的問題[23]。

偏最小二乘基本原理如下：

(1)設(shè)有p個(gè)自變量的集合X和m個(gè)因變量的集合Y，且每個(gè)變量都有n個(gè)樣本數(shù)據(jù)：

(1)

(2)

對(duì)原始數(shù)據(jù)X進(jìn)行標(biāo)準(zhǔn)化處理

(3)

(2)求取主成分

從A0、B0中提取第一主成分t1和u1(t1、u1分別是A0與B0線性組合)

(4)

(5)

式中：α1和β1為權(quán)重系數(shù)，可通過拉格朗日乘數(shù)法求得； cov(·，·)表示協(xié)方差。

根據(jù)主成分分析原理要求，t1、u1盡可能多地?cái)y帶了各自數(shù)據(jù)集中的原始信息，并且自變量的主成分t1與因變量的主成分u1相關(guān)度達(dá)到最大，即求取t1、u1最大協(xié)方差。

(3)建立原自變量、因變量與主成分之間的回歸關(guān)系

(6)

用殘差矩陣A1、B1取代A0和B0后，求第二主成分t2和u2并建立回歸關(guān)系，如果殘差矩陣B1中元素絕對(duì)值接近0，可認(rèn)為回歸方程達(dá)到預(yù)期精度并中止算法，否則繼續(xù)重復(fù)步驟(2)～步驟(3)。

(4)建立因變量與自變量間回歸方程

若A0的秩為a，則有a個(gè)成分t1，t2，…,ta使

(7)

tk=αk1x11，αk2x12，…，αkpx1p(k=1,2,…,a)

(8)

將式(8)代入Y=t1r1+t2r2+…+tara則得到m個(gè)因變量的回歸方程

yil=αl1x11+αl2x12+…+αlpxmp

(l=1,2,…,m)

(9)

(5)交叉有效性檢驗(yàn)

偏最小二乘回歸法一般采用交叉有效性檢驗(yàn)方法確定主成分提取個(gè)數(shù)，即通過增加一個(gè)新的成分，判斷其能否提高回歸方程的預(yù)測(cè)精度。具體實(shí)現(xiàn)如下。

i=1,2,…,nl=1,2,…,m

(10)

則Y的預(yù)測(cè)誤差平方和為

(11)

(12)

定義交叉有效性

(13)

基于上述原理，利用偏最小二乘回歸模型建立研究區(qū)資料樣本集之間的聯(lián)系。分別對(duì)研究區(qū)的三期地震資料進(jìn)行地震屬性提取，并將優(yōu)選出的地震屬性體轉(zhuǎn)換為沿井地震屬性曲線； 以測(cè)井滲透率數(shù)據(jù)作為因變量，并將其對(duì)應(yīng)期次的沿井地震數(shù)據(jù)作為自變量，采用偏最小二乘回歸算法建立預(yù)測(cè)模型，通過預(yù)測(cè)模型對(duì)不同期次滲透率進(jìn)行全區(qū)預(yù)測(cè)，得到三個(gè)含水時(shí)期滲透率曲線； 最后利用驗(yàn)證數(shù)據(jù)集中的不同含水期樣本井的原始滲透率數(shù)據(jù)與預(yù)測(cè)滲透率數(shù)據(jù)進(jìn)行誤差分析，并結(jié)合研究區(qū)示蹤劑資料驗(yàn)證滲透率時(shí)變預(yù)測(cè)結(jié)果。

4 基于時(shí)移地震資料的儲(chǔ)層滲透率時(shí)變性研究

樣本數(shù)量的不對(duì)等會(huì)造成預(yù)測(cè)結(jié)果的偏差，經(jīng)統(tǒng)計(jì)，區(qū)內(nèi)上世紀(jì)90年代前、2000年左右和2010年左右的投產(chǎn)井?dāng)?shù)分別為71、96和78口，數(shù)目相差不大。測(cè)井解釋成果與三期時(shí)移地震資料相匹配，分別代表中低、高和特高含水階段，能夠?yàn)閮?chǔ)層參數(shù)預(yù)測(cè)提供較好的數(shù)據(jù)基礎(chǔ)。通過編制相應(yīng)的軟件，對(duì)研究區(qū)內(nèi)不同含水階段測(cè)井序列數(shù)據(jù)進(jìn)行回歸計(jì)算，以測(cè)井曲線中0.125m采樣點(diǎn)為精度，建立包含三期沿井時(shí)移地震屬性數(shù)據(jù)和各井投產(chǎn)時(shí)期測(cè)井解釋的滲透率數(shù)據(jù)共104983個(gè)數(shù)據(jù)集。隨機(jī)選取數(shù)據(jù)集中80%的數(shù)據(jù)構(gòu)建訓(xùn)練集，通過偏最小二乘回歸法建立當(dāng)期沿井時(shí)移地震屬性與當(dāng)期測(cè)井解釋滲透率之間的關(guān)系及預(yù)測(cè)模型； 選取20%的數(shù)據(jù)構(gòu)建驗(yàn)證集，將已知測(cè)井解釋的原始滲透率與通過預(yù)測(cè)模型得到的同期次預(yù)測(cè)滲透率進(jìn)行對(duì)比，用來驗(yàn)證結(jié)果的準(zhǔn)確性。在模型訓(xùn)練完成后，分別繪制驗(yàn)證集中不同含水階段通過測(cè)井解釋得到的原始滲透率與通過偏最小二乘法得到的對(duì)應(yīng)期次井的預(yù)測(cè)滲透率交會(huì)圖(圖1)。對(duì)于驗(yàn)證集中的數(shù)據(jù)，原始測(cè)井解釋滲透率和當(dāng)期次預(yù)測(cè)滲透率數(shù)據(jù)間的相關(guān)性越高代表預(yù)測(cè)模型越準(zhǔn)確，計(jì)算得到圖1中低、高、特高含水階段層段內(nèi)原始滲透率與預(yù)測(cè)滲透率交會(huì)的相關(guān)系數(shù)值分別為0.91、0.87、0.89，說明建立的偏最小二乘回歸模型能夠很好地進(jìn)行滲透率預(yù)測(cè)，其驗(yàn)證結(jié)果與實(shí)際值相符。

圖1 不同含水時(shí)期原始滲透率和預(yù)測(cè)滲透率交會(huì)圖(a)中低含水時(shí)期； (b)高含水時(shí)期； (c)特高含水時(shí)期

為了細(xì)化研究區(qū)內(nèi)滲透率變化規(guī)律，選取G2-62-1井區(qū)進(jìn)行對(duì)比分析。該井區(qū)范圍內(nèi)主力開發(fā)小層包括NgⅠ1-1和NmⅣ3-1，兩主力層均無明顯斷層分布，且構(gòu)造較為平緩。從物性上看，NgⅠ1-1小層為中孔、中滲儲(chǔ)層，井區(qū)范圍內(nèi)原始滲透率平均值為132mD； 而NmⅣ3-1小層在該井控制區(qū)域內(nèi)平均原始滲透率為7293mD，屬于高孔、高滲儲(chǔ)層。G2-62-1井區(qū)范圍內(nèi)兩個(gè)主力層滲透率分布圖顯示(圖2)，在同一小層內(nèi)滲透率差異較小，NgⅠ1-1和NmⅣ3-1層的滲透率變異系數(shù)分別為0.14和0.23，按照滲透率變異系數(shù)分類[24]，兩小層在該井區(qū)范圍內(nèi)均屬于弱非均質(zhì)儲(chǔ)層。

圖2 G2-62-1井區(qū)范圍內(nèi)NmⅣ3-1(左)、NgⅠ1-1(右)小層原始滲透率分布

將G2-62-1井兩期預(yù)測(cè)滲透率曲線和一期測(cè)井原始滲透率曲線成果進(jìn)行粗化，構(gòu)建形成中低、高和特高含水期滲透率剖面圖(圖3)，對(duì)比不同含水階段滲透率的變化。由圖可見，兩套儲(chǔ)層受注水沖刷后，其滲透率變化并非簡(jiǎn)單的單調(diào)增加或減小，而是與巖石膠結(jié)程度和泥質(zhì)含量等因素有關(guān)[25-28]。從圖3預(yù)測(cè)結(jié)果判斷，該井控制范圍內(nèi)明化鎮(zhèn)組和館陶組的儲(chǔ)層滲透率變化存在差異，其中NgⅠ1-1小層經(jīng)注水沖刷后，滲透率變化表現(xiàn)為緩慢減小的趨勢(shì)，由中低含水期的154mD變?yōu)樘馗吆诘?20mD； 而NmⅣ3-1小層經(jīng)注水沖刷后滲透率變化表現(xiàn)為增加的趨勢(shì)，由中低含水期的7211mD變?yōu)樘馗吆诘?903mD。

圖3 G2-62-1井三期滲透率剖面圖

郭莉等[24]對(duì)油田注水開發(fā)過程中儲(chǔ)層參數(shù)的變化規(guī)律按類別進(jìn)行了分析，研究表明對(duì)于高孔、高滲儲(chǔ)層，其黏土礦物和膠結(jié)物含量隨注水沖刷而減少，溶蝕孔隙增多，整體表現(xiàn)為儲(chǔ)層孔喉增大、物性變好的趨勢(shì)。對(duì)于中孔、中滲油藏則出現(xiàn)注水沖刷導(dǎo)致喉道堵塞、孔喉減小、物性變差的趨勢(shì)，注入水進(jìn)入油層后，由于水對(duì)黏土礦物的聚散和水化膨脹作用，大孔道中的黏土礦物隨水流帶出被沖散、遷移，使孔道變得暢通，擴(kuò)大了喉道直徑； 另一方面，一些被剝落或沖散的黏土在小孔隙中重新聚集，使小孔隙變得更小，增強(qiáng)了儲(chǔ)層微觀非均質(zhì)性。采用常規(guī)巖心驅(qū)替實(shí)驗(yàn)分析單井三期滲透率變化差異的原因，得到的儲(chǔ)層參數(shù)變化受到巖心尺度、實(shí)驗(yàn)條件、實(shí)驗(yàn)樣本數(shù)等因素的影響，只能反映單一巖心儲(chǔ)層參數(shù)變化規(guī)律。然而實(shí)際油藏儲(chǔ)層物性分布復(fù)雜，時(shí)移地震作為一種能夠反映全區(qū)參數(shù)場(chǎng)分布的資料，在儲(chǔ)層參數(shù)預(yù)測(cè)方面能夠發(fā)揮重要作用。圖3結(jié)果與高孔、高滲及中孔、中滲油藏驅(qū)替實(shí)驗(yàn)結(jié)果相符，進(jìn)一步說明通過時(shí)移地震資料進(jìn)行儲(chǔ)層參數(shù)預(yù)測(cè)比單一巖心驅(qū)替實(shí)驗(yàn)結(jié)果更具可靠性。

5 成果驗(yàn)證

利用示蹤劑資料對(duì)比一個(gè)井組內(nèi)不同方向注水竄流速度，進(jìn)一步驗(yàn)證滲透率預(yù)測(cè)精度和準(zhǔn)確性。在G2-62-1井明化鎮(zhèn)組注入示蹤劑檢測(cè)資料，表2記錄了三口監(jiān)測(cè)井水驅(qū)速度數(shù)據(jù)，G1-62井水驅(qū)速度為1.76m/d，速度較慢； G2-61-1井未見示蹤劑； 而監(jiān)測(cè)井G1-61-2在短時(shí)間內(nèi)檢測(cè)出示蹤劑，計(jì)算水驅(qū)速度達(dá)到7.49m/d。

表2 G2-62-1注采井組水驅(qū)速度監(jiān)測(cè)

選取G2-62-1注采井組范圍內(nèi)各井1993年和2015年滲透率曲線，用NmⅣ3-1小層2015年單井滲透率數(shù)值除以1993年滲透率數(shù)值，得到各井點(diǎn)處滲透率變化倍數(shù)。以砂體分布作為約束條件對(duì)井點(diǎn)數(shù)據(jù)進(jìn)行平面插值，得到NmⅣ3-1小層中低含水期—特高含水期的滲透率變化倍數(shù)分布圖(圖4)。其中數(shù)值小于1表示經(jīng)過注水沖刷后儲(chǔ)層物性變差，從中低含水期到特高含水期滲透率數(shù)值變?。?數(shù)值大于1表示從中低含水期到特高含水期儲(chǔ)層滲透率數(shù)值變大。

圖4 G2-62-1井區(qū)NmⅣ3-1層滲透率變化倍數(shù)分布

由圖2已知NmⅣ3-1小層整體原始滲透率分布差異較小，砂體內(nèi)部無大的斷層分布，滲透率變化與注水沖刷量有直接關(guān)系，而物性的變化又影響了儲(chǔ)層滲流特征。由于G2-61-1井處于NmⅣ3-1小層有效砂體外部，與G2-62-1注水井無直接連通關(guān)系，因此未見示蹤劑響應(yīng)。G1-62井有示蹤劑顯示，該井與G2-62-1注水井間水驅(qū)速度為1.76m/d，滿足達(dá)西定律，同時(shí)從圖4可以看出，G1-62井位于滲透率變化明顯區(qū)域的外部，未與注水井形成直接連通的滲流通道。G1-61-2井特高含水時(shí)期滲透率較中低含水時(shí)期增大近三倍，受G2-62-1注水井影響，G2-62-1與G1-61-2兩井之間的儲(chǔ)層滲透率發(fā)生變化，形成高滲通道，導(dǎo)致水驅(qū)速度為7.49m/d，遠(yuǎn)高于達(dá)西流的滲流速度，兩井間滲透率變化已呈現(xiàn)高速非達(dá)西流的特征。

表3對(duì)比了井組內(nèi)各單井控制范圍內(nèi)的水驅(qū)速度與滲透率變化倍數(shù)之間的關(guān)系，直接對(duì)預(yù)測(cè)結(jié)果進(jìn)行驗(yàn)證。通過時(shí)移地震資料預(yù)測(cè)，G1-62井控制區(qū)域內(nèi)滲透率并未見明顯變化，其水驅(qū)速度符合達(dá)西定律； G1-61-2井與G2-62-1井控區(qū)域內(nèi)預(yù)測(cè)滲透率變化倍數(shù)分別為2.7和3.4，井間儲(chǔ)層呈現(xiàn)孔喉增大、物性變好的趨勢(shì)，其水驅(qū)速度明顯高于G1-62井，呈現(xiàn)高速非達(dá)西流的趨勢(shì)，說明基于多期時(shí)移地震屬性預(yù)測(cè)的滲透率變化結(jié)果與水驅(qū)開發(fā)特征相一致。

表3 G2-62-1注采井組預(yù)測(cè)滲透率變化倍數(shù)與水驅(qū)速度對(duì)比

6 結(jié)論

(1)地震屬性參數(shù)的合理優(yōu)選直接影響預(yù)測(cè)結(jié)果的準(zhǔn)確度和可靠性。本文提取研究區(qū)時(shí)移地震資料的10種地震屬性進(jìn)行了相關(guān)性分析，優(yōu)選其中均方根振幅、反射強(qiáng)度、瞬時(shí)頻率、振幅包絡(luò)和主頻5種與滲透率相關(guān)性較高的屬性，為下一步時(shí)移地震多屬性預(yù)測(cè)滲透率時(shí)變規(guī)律提供了數(shù)據(jù)基礎(chǔ)。

(2)研究選用的偏最小二乘回歸算法較普通多元回歸法能夠更好地處理自變量之間存在的多重相關(guān)性以及數(shù)據(jù)量小、樣本較少的問題，應(yīng)用回歸模型建立多自變量和因變量間的回歸關(guān)系并對(duì)模型進(jìn)行預(yù)測(cè)，確定了研究區(qū)低、高、特高含水階段原始滲透率和預(yù)測(cè)滲透率的相關(guān)系數(shù)值分別達(dá)到0.91、0.87、0.89，建立的回歸模型能夠很好地對(duì)滲透率進(jìn)行預(yù)測(cè)。

(3)不同類型儲(chǔ)層滲透率時(shí)變規(guī)律存在差異，將建立的偏最小二乘回歸預(yù)測(cè)模型應(yīng)用于研究區(qū)NgⅠ1-1中孔、中滲儲(chǔ)層，其滲透率變化表現(xiàn)為緩慢減小的趨勢(shì)，由中低含水期的154mD變?yōu)樘馗吆诘?20mD； 而NmⅣ3-1高孔、高滲儲(chǔ)層滲透率變化表現(xiàn)為增加的趨勢(shì)，由中低含水期的7211mD變?yōu)樘馗吆诘?903mD。通過示蹤劑資料進(jìn)一步證明，建立的多期沿井滲透率預(yù)測(cè)模型在高含水油藏儲(chǔ)層精細(xì)描述過程中具有可借鑒作用，該方法為指導(dǎo)高含水開發(fā)后期剩余油定量挖潛提供了依據(jù)。