＊張怡然 李?yuàn)W

（長(zhǎng)江大學(xué)地球物理與石油資源學(xué)院 湖北 430100）

引言

為了使油氣生產(chǎn)更加科學(xué)化、高效化，對(duì)測(cè)井方案的設(shè)計(jì)是進(jìn)行生產(chǎn)測(cè)井前的關(guān)鍵環(huán)節(jié)，設(shè)計(jì)環(huán)節(jié)中需要選取合適的測(cè)井儀器，估算產(chǎn)量，對(duì)井下流體進(jìn)行流型判斷等。其中，對(duì)井下流體流型的判斷是非常重要的步驟，流體流型判斷直接關(guān)系到后續(xù)對(duì)整個(gè)測(cè)井方案的制訂以及對(duì)測(cè)井儀器的選用，避免因?yàn)榉桨覆缓线m導(dǎo)致的浪費(fèi)。

近年來(lái)隨著計(jì)算機(jī)技術(shù)的不斷進(jìn)步，人們通過(guò)對(duì)已有的數(shù)據(jù)集進(jìn)行機(jī)器學(xué)習(xí)，能夠很好的對(duì)新的數(shù)據(jù)進(jìn)行預(yù)測(cè)，國(guó)內(nèi)外已經(jīng)有運(yùn)用機(jī)器學(xué)習(xí)神經(jīng)網(wǎng)絡(luò)等手段對(duì)流體物性參數(shù)（PVT）進(jìn)行處理預(yù)測(cè)的先例[1]，因而利用機(jī)器學(xué)習(xí)方案處理井下數(shù)據(jù)，進(jìn)而對(duì)井下流體流型進(jìn)行預(yù)測(cè)不失為一種可行的方案。同時(shí)通過(guò)將傳統(tǒng)測(cè)井與機(jī)器學(xué)習(xí)算法相結(jié)合，有助于整個(gè)測(cè)井行業(yè)進(jìn)入智能化、高效化的新發(fā)展階段，更加符合新時(shí)代下傳統(tǒng)產(chǎn)業(yè)數(shù)字化轉(zhuǎn)型，促進(jìn)新技術(shù)在傳統(tǒng)行業(yè)中發(fā)揮作用。

1.算法原理

(1)BP神經(jīng)網(wǎng)絡(luò)分類算法

其基本原理如下[2,3]：利用多層感知機(jī)（Multilayer Perceptron，MLP），根據(jù)現(xiàn)有的大量已知數(shù)據(jù)，將每條數(shù)據(jù)的特征匯總成數(shù)據(jù)集X={x1，x2……xn}，每條數(shù)據(jù)都對(duì)應(yīng)一個(gè)標(biāo)簽y，通過(guò)對(duì)數(shù)據(jù)集上訓(xùn)練來(lái)學(xué)習(xí)函數(shù)f(·):Rn→Ro，其中，n為數(shù)據(jù)集中數(shù)據(jù)的維數(shù)；o為輸出后的維數(shù)。和傳統(tǒng)線性回歸不同，在通過(guò)MLP進(jìn)行學(xué)習(xí)時(shí)，輸入層和輸出層之間，會(huì)存在一個(gè)或多個(gè)非線性層，稱為隱藏層，圖1為一個(gè)具有兩層隱藏層的MLP。

圖1 具有兩層隱藏層的BP神經(jīng)網(wǎng)絡(luò)

最左邊的層叫輸入層，將已有數(shù)據(jù)標(biāo)準(zhǔn)化后，作為代表輸入特征的神經(jīng)元；然后經(jīng)過(guò)一層或多層隱藏層中的神經(jīng)元將前一層的數(shù)據(jù)加權(quán)線性求和轉(zhuǎn)換w1x1+w2x2+……+wnxn，并通過(guò)非線性函數(shù)g(·)R→R，g(·)也被稱為激活函數(shù)，通常采用：

最后由輸出層接收最后一個(gè)隱藏層轉(zhuǎn)換后的值。

MLP采用了隨機(jī)梯度下降（Stochastic Gradient Descent）進(jìn)行訓(xùn)練，將需要調(diào)整的參數(shù)通過(guò)使用損失函數(shù)進(jìn)行梯度更新，以獲取更新后的參數(shù)：

其中，η為控制參數(shù)空間搜索步長(zhǎng)的學(xué)習(xí)速率；Loss為神經(jīng)網(wǎng)絡(luò)的損失函數(shù)。

MLP采用均方差損失函數(shù)：

在最初的隨機(jī)化權(quán)重時(shí)，多層感知機(jī)通過(guò)重復(fù)更新權(quán)重，使得損失的值最小，在計(jì)算損失之后，反向傳遞至輸出層的前一層，為每一個(gè)權(quán)重參數(shù)提供一個(gè)更新值，從而減少誤差。

隨著梯度下降，計(jì)算得到對(duì)每一個(gè)權(quán)重所損失的梯度▽Lossw，并從總損失中減去，即：

由于在對(duì)井下流體流型的預(yù)測(cè)中，井下流體參數(shù)并不是唯一的，因而需要解決多個(gè)類同時(shí)存在的情況，故在隱藏層中學(xué)習(xí)函數(shù)f(x)=W2g(W1Tx+b1)+b2本身為一個(gè)n維向量，通過(guò)Softmax函數(shù)：

其中，zi為Softmax中輸入的第i個(gè)元素，其對(duì)應(yīng)與相應(yīng)的第i類；K是所有類別的總數(shù)量，輸出的值包含了樣本中屬于每一個(gè)類的概率的向量，最終輸出的結(jié)果為概率最高的類。

(2)樹(shù)模型分類算法

其基本原理如下[4]：通過(guò)創(chuàng)造一個(gè)模型，使程序?qū)υ撃Ｐ瓦M(jìn)行學(xué)習(xí)，最終從模型數(shù)據(jù)特征中推斷出的簡(jiǎn)單決策規(guī)則來(lái)預(yù)測(cè)新的輸入變量的值。

將已有的數(shù)據(jù)標(biāo)準(zhǔn)化后分條作為訓(xùn)練向量xi∈Rn，i=1,……,l和類向量y∈Rl，利用決策樹(shù)遞歸劃分空間，將相同類別的樣本劃分至同一個(gè)組，節(jié)點(diǎn)m處的數(shù)據(jù)集用Q表示，對(duì)于一個(gè)由特征j和閾值tm組成的階段劃分?jǐn)?shù)據(jù)集θ=(j,tm)，將數(shù)據(jù)分別作為兩個(gè)分支：

節(jié)點(diǎn)處的不存度用不存度函數(shù)H(·)計(jì)算：

其中，Pmk表示在具有Nm個(gè)觀測(cè)值的區(qū)域Rm中，節(jié)點(diǎn)m在k類中所占的比例，故Pmk公式為：

Xm為節(jié)點(diǎn)m中所包含的已知數(shù)據(jù)。

井下流體包含多種流體種類，因而得出的結(jié)論也是多種類的，故需要將類向量集定義為一個(gè)二維的數(shù)組；對(duì)每一種輸出的結(jié)果需要建立一個(gè)估計(jì)器，依據(jù)每一種類別的輸出平均減少量來(lái)作為每一棵分支的產(chǎn)生標(biāo)準(zhǔn)。

2.方法應(yīng)用

由于實(shí)際生產(chǎn)情況具有多種不同井況條件，因而分別采用了在井斜角度0°、60°、85°、90°下，含水率20%、40%、60%、80%、90%，流量100m3/d，300m3/d，600m3/d的實(shí)驗(yàn)數(shù)據(jù)，各個(gè)數(shù)據(jù)間差異較大，如果直接使用會(huì)造成部分?jǐn)?shù)據(jù)所占權(quán)重過(guò)大影響最終結(jié)果的情況，因而需要對(duì)數(shù)據(jù)進(jìn)行標(biāo)準(zhǔn)化處理，采用離差標(biāo)準(zhǔn)化處理，使得所有數(shù)據(jù)結(jié)果值映射到[0,1]區(qū)間內(nèi)，所用函數(shù)如下：

其中，xmax為數(shù)據(jù)中的最大值；xmin為數(shù)據(jù)中最小值。

然后將標(biāo)準(zhǔn)化處理的數(shù)據(jù)集中每個(gè)數(shù)據(jù)分別進(jìn)行加權(quán)線性求和，利用損失函數(shù)計(jì)算損失后反向傳播至前一層以減小誤差，最終達(dá)到設(shè)定的最大迭代數(shù)時(shí)結(jié)束計(jì)算，BP分類算法結(jié)束。

再對(duì)標(biāo)準(zhǔn)化處理的數(shù)據(jù)進(jìn)行評(píng)估，選出其中最具有決定性的特征，然后將原始數(shù)據(jù)集劃分為幾個(gè)子集，這些子集分布在第一個(gè)決策點(diǎn)的所有分支上，如果一個(gè)分支上所有數(shù)據(jù)都是同一類別，則不需要繼續(xù)分類，若包含有不同類別，則需要繼續(xù)選出最具有決定性特征來(lái)繼續(xù)劃分，直至所有葉子節(jié)點(diǎn)上的類同屬于一類特征，圖2為決策樹(shù)模型運(yùn)算結(jié)果，決策樹(shù)算法結(jié)束[5]。

圖2 良好決策節(jié)點(diǎn)進(jìn)行決策分支的決策樹(shù)

3.實(shí)驗(yàn)

(1)設(shè)計(jì)實(shí)驗(yàn)

實(shí)驗(yàn)所用的數(shù)據(jù)包括實(shí)驗(yàn)數(shù)據(jù)和實(shí)際數(shù)據(jù)，將實(shí)驗(yàn)數(shù)據(jù)作為數(shù)據(jù)集以供機(jī)器學(xué)習(xí)算法學(xué)習(xí)，然后將實(shí)際數(shù)據(jù)輸入學(xué)習(xí)后的算法進(jìn)行預(yù)測(cè)，將預(yù)測(cè)后的結(jié)果和實(shí)際流體流型進(jìn)行對(duì)比，以檢驗(yàn)算法的可行性。

實(shí)驗(yàn)數(shù)據(jù)包含井斜角度、流體流量（m3/d）、含水率、溫度、管徑，將實(shí)際數(shù)據(jù)帶入至通過(guò)實(shí)驗(yàn)數(shù)據(jù)訓(xùn)練學(xué)習(xí)后的模型，得到對(duì)實(shí)際數(shù)據(jù)預(yù)測(cè)的流體流型，與實(shí)際結(jié)果進(jìn)行對(duì)照；通過(guò)分析不同流量、井斜、含水率下的準(zhǔn)確度，判斷算法的可行性。

(2)預(yù)測(cè)結(jié)果分析

在學(xué)習(xí)算法結(jié)束之后，隨機(jī)選取不同井斜、含水率、流量的12組已知數(shù)據(jù)，對(duì)數(shù)據(jù)進(jìn)行預(yù)處理后，然后使用訓(xùn)練過(guò)的神經(jīng)網(wǎng)絡(luò)進(jìn)行預(yù)測(cè)，將預(yù)測(cè)結(jié)果與實(shí)際數(shù)據(jù)進(jìn)行對(duì)比，以檢驗(yàn)BP分類算法和樹(shù)模型分類算法對(duì)流型預(yù)測(cè)的準(zhǔn)確率。

BP分類算法和樹(shù)模型預(yù)測(cè)結(jié)果如表1所示。

表1 BP分類算法與樹(shù)模型預(yù)測(cè)結(jié)果

從表1可知，在井斜角度較低時(shí)，BP分類算法得到的預(yù)測(cè)結(jié)果較井斜角度較大時(shí)的準(zhǔn)確度有大幅度提升，在井斜角度較大時(shí)預(yù)測(cè)準(zhǔn)確率較低，總的預(yù)測(cè)準(zhǔn)確度在75%左右。而從圖3來(lái)看，決策樹(shù)模型在決策點(diǎn)選擇較好的情況下，在各流量、各井斜、各含水率上均有較為準(zhǔn)確的預(yù)測(cè)結(jié)果，總體準(zhǔn)確度能達(dá)到90%以上。

圖3 兩種算法與實(shí)際流型結(jié)果比較

4.結(jié)論

將實(shí)驗(yàn)數(shù)據(jù)分別按照BP分類算法和樹(shù)模型算法進(jìn)行預(yù)測(cè)后與實(shí)際流體流型對(duì)比后可以得到如下結(jié)論：

（1）通過(guò)機(jī)器學(xué)習(xí)BP分類算法和樹(shù)模型算法進(jìn)行井筒流體流型預(yù)測(cè)，對(duì)井下流體流型預(yù)測(cè)提供了新的解決方案。其中BP分類算法在井斜角度較大的情況下能取得95%以上的準(zhǔn)確識(shí)別率，但是在其他條件下效果較差；樹(shù)模型算法在各種井斜、流量、含水率下均有較為準(zhǔn)確的預(yù)測(cè)，綜合準(zhǔn)確度在91.7%左右。

（2）通過(guò)較為合適的數(shù)據(jù)標(biāo)準(zhǔn)化處理，可以降低某些數(shù)據(jù)在預(yù)測(cè)中占有過(guò)高權(quán)重，減少了數(shù)據(jù)處理后的誤差，提高了數(shù)據(jù)處理的效率。

（3）通過(guò)不斷完善訓(xùn)練數(shù)據(jù)集的數(shù)據(jù)，可以不斷提升預(yù)測(cè)結(jié)果的準(zhǔn)確度，為后續(xù)該算法在實(shí)際應(yīng)用中打下良好基礎(chǔ)。