董興朋

（中國地震局地震研究所，湖北武漢430071）

相似度－遺傳神經(jīng)網(wǎng)絡(luò)在儲層物性預(yù)測中的應(yīng)用

董興朋

（中國地震局地震研究所，湖北武漢430071）

傳統(tǒng)的測井解釋需要建立精確的數(shù)學模型，并常伴有嚴格的條件限制，因此很難得到真實反映儲層特性的結(jié)果。采用遺傳算法與BP神經(jīng)網(wǎng)絡(luò)相結(jié)合，利用遺傳算法的全局尋優(yōu)特點，優(yōu)化神經(jīng)網(wǎng)絡(luò)的連接權(quán)值和閥值，提高網(wǎng)絡(luò)的訓(xùn)練精度和預(yù)測精度，避免了BP算法易陷入局部極小的缺點，提高運算速度。將相似度的概念引入到測井中，定義相似度在測井中的計算公式，提出相似度與遺傳神經(jīng)網(wǎng)絡(luò)相結(jié)合的方法。實例研究表明，預(yù)測準確性較高，可以有效控制預(yù)測精度，避免因儲層差別大而造成的預(yù)測精度降低的現(xiàn)象。

地球物理測井；遺傳算法；相似度；神經(jīng)網(wǎng)絡(luò)；儲層物性；預(yù)測

0 引 言

通過巖心取樣進行室內(nèi)測試，獲取巖心的孔隙度精度很高，但取樣和測試的費用很高，致使獲取的儲層孔隙度數(shù)據(jù)有限，不能覆蓋整個工區(qū)［1］。常規(guī)測井解釋多通過經(jīng)驗公式或簡化地質(zhì)條件建立模型計算儲層參數(shù)。然而測井資料具有不確定性和非結(jié)構(gòu)化的特點，過去采用的統(tǒng)計方法未能擺脫線性方程的束縛，使求得的物性參數(shù)很難與實驗數(shù)據(jù)建立良好的對應(yīng)關(guān)系［2］。

BP神經(jīng)網(wǎng)絡(luò)［3］的學習基于梯度下降法，因此存在收斂速度慢及網(wǎng)絡(luò)參數(shù)和計算參數(shù)難以確定等缺點［4］。遺傳算法（Genetic Algorithm，GA）是一種借鑒生物界自然選擇和自然遺傳機制的搜索算法，它能在復(fù)雜空間中尋找最優(yōu)解，且有算法簡單、適用、魯棒性強等優(yōu)點?；谝陨?種算法各自的優(yōu)缺點，可以將遺傳算法和BP神經(jīng)網(wǎng)絡(luò)結(jié)合，提高網(wǎng)絡(luò)訓(xùn)練精度。利用GA－BP神經(jīng)網(wǎng)絡(luò)預(yù)測儲層物性參數(shù)，具有綜合利用多種測井參數(shù)，客觀反映輸入?yún)?shù)與儲層參數(shù)之間復(fù)雜規(guī)律和無需選擇測井解釋參數(shù)等特點。

本文將相似度的概念引入到測井中，提出了相似度與遺傳神經(jīng)網(wǎng)絡(luò)相結(jié)合，并在MATLAB中進行預(yù)測，可以有效控制預(yù)測精度，避免因儲層差別大而造成的神經(jīng)網(wǎng)絡(luò)預(yù)測精度降低的現(xiàn)象。

1 遺傳算法與神經(jīng)網(wǎng)絡(luò)

1.1 標準BP算法

BP（Back Propagation）網(wǎng)絡(luò)是一種按誤差逆?zhèn)鞑ニ惴ㄓ?xùn)練的多層前饋網(wǎng)絡(luò)，能學習和存儲大量的輸入－輸出模式映射關(guān)系，無需事前揭示描述這種映射關(guān)系的數(shù)學方程。它的學習規(guī)則是最速下降法，通過反向傳播不斷調(diào)整網(wǎng)絡(luò)的權(quán)值和閥值，使網(wǎng)絡(luò)的誤差平方和最小。理論證明，含有一個隱含層的BP網(wǎng)絡(luò)可以實現(xiàn)以任意精度近似任何連續(xù)非線性函數(shù)［5］。但在實際應(yīng)用過程中，存在易陷入局部極小和收斂速度慢等缺點。針對這些缺點，已有人對此提出了改進方案，如修正學習率、增加動量項、引入陡度因子及改進神經(jīng)元激勵函數(shù)等［3］。

1.2 遺傳算法原理

遺傳算法是借助生物的進化規(guī)律，使所要解決的問題從初始解一步步逼近最優(yōu)解［6］。它是一種全局優(yōu)化搜索算法，可以有效克服BP神經(jīng)網(wǎng)絡(luò)易陷入局部最優(yōu)的不足。遺傳算法一般包括選擇、交叉和變異3個基本操作。適應(yīng)度函數(shù)（目標函數(shù)）被用來評價個體解的優(yōu)劣程度，從而對個體進行選擇操作?？刂艷A處理效果的主要參數(shù)包括群體規(guī)模與交叉概率、變異概率。其流程見圖1。

圖1 遺傳算法流程圖

1.3 權(quán)值優(yōu)化

用GA優(yōu)化神經(jīng)網(wǎng)絡(luò)可以使得神經(jīng)網(wǎng)絡(luò)具有自進化、自適應(yīng)能力，它主要包括網(wǎng)絡(luò)權(quán)值的進化、網(wǎng)絡(luò)拓撲結(jié)構(gòu)的進化和學習規(guī)則的進化，而最主要的是用來優(yōu)化網(wǎng)絡(luò)的權(quán)值［7］。遺傳算法優(yōu)化神經(jīng)網(wǎng)絡(luò)權(quán)值的一般步驟如下。

（1）確定BP神經(jīng)網(wǎng)絡(luò)的拓撲結(jié)構(gòu)和訓(xùn)練樣本。

（2）對神經(jīng)網(wǎng)絡(luò)權(quán)值編碼，遺傳控制參數(shù)初始化，產(chǎn)生含P個個體的初始種群（權(quán)值）。

（3）適應(yīng)度函數(shù)的確定。一般以誤差函數(shù)的倒數(shù)作為染色體評價函數(shù)。誤差越大適應(yīng)度值越小。

（4）計算每個個體的適應(yīng)度，選擇若干適應(yīng)度函數(shù)值大的個體直接進入小一代；適應(yīng)度小的個體被淘汰。

（5）利用交叉、變異等遺傳算子對當前新一代群體進行新一輪迭代，直到訓(xùn)練目標滿足終止條件為止。重復(fù)步驟（2）、（3）、（4）、（5），對新一代群體進行新一輪迭代，直到訓(xùn)練目標滿足終止條件為止。

2 相似度在測井中的應(yīng)用

2.1 儲層類型

影響儲層孔隙度的地質(zhì)因素有埋藏深度、構(gòu)造位置、沉積環(huán)境、巖性變化、成巖程度等。沉積環(huán)境對儲層物性的影響一直是定性說明，由于在不同的沉積環(huán)境中會形成不同的沉積序列，也就會有不同的砂地體積比值［8］。最常見的有3類：沖積扇沙礫巖儲層、河流砂體儲層和三角洲砂體儲層。在相同的儲層類型中，由于沉積環(huán)境相同，可以認為測井信息與儲層參數(shù)之間的非線性映射是相似的。

2.2 相似度在測井中的定義

在相同的沉積環(huán)境（即沖積扇沙礫巖儲層、河流砂體儲層和三角洲砂體儲層3類），取標準的儲層模型i，該模型指標特征值的向量μi＝（μi1，μi2，…，μik），其中，μi1，μi2，…，μik分別為儲層k種不同測井曲線值；取所要預(yù)測的儲層為樣本j，指標特征值向量為rj＝（rj1，rj2，…，rjk），k種不同的測井技術(shù)對儲層物性影響的權(quán)重分別為wj1，wj2，…，wjk則樣本j相對標準模型i的相似度為

有了相似度的定義，在進行儲層相似度預(yù)測時就可以先計算所要預(yù)測儲層與已知模型的相似度。若相似程度很高，就可以直接用經(jīng)過已知模型訓(xùn)練的神經(jīng)網(wǎng)絡(luò)對未知儲層的孔隙度進行預(yù)測。

3 GA－BP網(wǎng)絡(luò)模型及MATLAB實現(xiàn)

用神經(jīng)網(wǎng)絡(luò)預(yù)測儲層物性參數(shù)，就是尋求測井信息與儲層物性參數(shù)之間的某種非線性映射。BP網(wǎng)絡(luò)是一種按誤差逆?zhèn)鞑ニ惴ㄓ?xùn)練的多層前饋網(wǎng)絡(luò)，很好地解決了輸入和輸出之間的關(guān)系和復(fù)雜的非線性問題，為儲層參數(shù)的預(yù)測提供了一種有效的途徑［2］。

MATLAB中的神經(jīng)網(wǎng)絡(luò)工具箱以神經(jīng)網(wǎng)絡(luò)為基礎(chǔ)，利用MATLAB腳本語言構(gòu)造出典型的神經(jīng)網(wǎng)絡(luò)激活函數(shù)，使設(shè)計者對所選網(wǎng)絡(luò)的計算變成對激活函數(shù)的調(diào)用，為用戶提供了極大方便［2］?；贛ATLAB的GA－BP網(wǎng)絡(luò)實現(xiàn)儲層物性預(yù)測的一般步驟如下。

（1）原始資料的預(yù)處理。巖心實驗數(shù)據(jù)及測井數(shù)據(jù)的預(yù)處理包括巖心數(shù)據(jù)和測井數(shù)據(jù)的深度歸位，使得巖心分析所得的物性參數(shù)和測井信息來源于同一深度。

（2）選取學習樣本。

（3）數(shù)據(jù)的歸一化處理。使用premnmx函數(shù)進行歸一化處理，使數(shù)據(jù)分布在［－1～1］之間。

（4）構(gòu)建預(yù)測網(wǎng)絡(luò)模型。根據(jù)神經(jīng)網(wǎng)絡(luò)的連續(xù)函數(shù)映射定理，具有3層（1個輸入層、1個隱含層和1個輸出層）的網(wǎng)絡(luò)結(jié)構(gòu)，即可實現(xiàn)任意連續(xù)函數(shù)從輸入空間向輸出空間的映射［3］。本文采用3層的神經(jīng)網(wǎng)絡(luò)結(jié)構(gòu)，輸入神經(jīng)元數(shù)為測井參數(shù)的個數(shù)，輸出神經(jīng)元為1個（孔隙度或滲透率），隱含層節(jié)點數(shù)較難確定，節(jié)點數(shù)太多，網(wǎng)絡(luò)訓(xùn)練時間長；節(jié)點太少，誤差精度又達不到要求。用學習樣本集訓(xùn)練網(wǎng)絡(luò)，直到精度滿足要求為止，保存網(wǎng)絡(luò)，建立儲層參數(shù)預(yù)測網(wǎng)絡(luò)模型。

（5）應(yīng)用上述建立的預(yù)測網(wǎng)絡(luò)模型預(yù)測未知儲層層段的物性參數(shù)。

4 孔隙度的預(yù)測

4.1 基于GA－BP神經(jīng)網(wǎng)絡(luò)的孔隙度預(yù)測

以寧東5井延長組目的層段的7個樣本為訓(xùn)練樣本作為神經(jīng)網(wǎng)絡(luò)的輸入，巖心孔隙度作為訓(xùn)練樣本輸出值，選取反映儲層孔隙度參數(shù)敏感的聲波（AC）、密度（DEN）、中子測井（CNL）和自然伽馬（GR）等4種測井曲線作為輸入元，隱含層的神經(jīng)元數(shù)為25，輸出層神經(jīng)元為1，建立GA－BP神經(jīng)網(wǎng)絡(luò)模型（樣本訓(xùn)練數(shù)據(jù)見表1）。

表1 樣本訓(xùn)練數(shù)據(jù)

GA－BP神經(jīng)網(wǎng)絡(luò)染色體編碼采用實數(shù)編碼，適應(yīng)度函數(shù)為誤差平方和的倒數(shù)。種群規(guī)模為50，遺傳代數(shù)為100，網(wǎng)絡(luò)最大訓(xùn)練次數(shù)為1 000，訓(xùn)練誤差精度為0.000 1。用訓(xùn)練樣本訓(xùn)練網(wǎng)絡(luò)，經(jīng)過100代遺傳操作，得到遺傳算法優(yōu)化的BP網(wǎng)絡(luò)初始權(quán)值和閥值，適應(yīng)度函數(shù)值為2 300（見圖2）；再用共軛梯度BP算法訓(xùn)練網(wǎng)絡(luò)，經(jīng)過35次迭代，均方差誤差達到訓(xùn)練誤差精度0.000 1（見圖3）。

取3組數(shù)據(jù)作為測試樣本進行預(yù)測。取訓(xùn)練樣本與測井輸入曲線的平均值作為標準儲層模型的指標特征值，根據(jù)式（1），取權(quán)重均為0.25，計算3組測試樣本的相似度。預(yù)測結(jié)果見表2，誤差隨相似度變化曲線見圖4。

由表2和圖4知，當所預(yù)測樣本與訓(xùn)練樣本的相似度高時，預(yù)測精度很高；反之，預(yù)測精度不高，存在較大的誤差。隨著相似度的增大，GA－BP神經(jīng)網(wǎng)絡(luò)預(yù)測的誤差逐漸減小。因此，在孔隙度預(yù)測之前，須先計算預(yù)測樣本的相似度，以判斷是否適合，這樣，可以有效控制預(yù)測精度，避免因儲層差別很大而造成的神經(jīng)網(wǎng)絡(luò)預(yù)測精度降低的現(xiàn)象。

圖4 誤差隨相似度的變化曲線

表2 GA－BP神經(jīng)網(wǎng)絡(luò)孔隙度預(yù)測數(shù)據(jù)表

4.2 GA－BP算法與標準BP算法預(yù)測結(jié)果對比

同樣，以表2中3組樣本為測試樣本，采用標準BP算法對樣本孔隙度值進行預(yù)測（預(yù)測結(jié)果見表3、圖5）。由表3可知，采用標準BP算法得到的預(yù)測值的誤差高于采用GA－BP算法的誤差。樣本1由于相似度較低，2種算法得到的誤差都較高；樣本2和3采用GA－BP算法得到的預(yù)測值的誤差遠小于標準BP算法，證明了GA－BP算法相對于傳統(tǒng)的標準BP算法的優(yōu)勢；圖5中，GA－BP算法和標準BP算法預(yù)測的誤差都隨著相似度的增加而減小，證明了將相似度引入到儲層孔隙度預(yù)測中的必要性。

表3 GA－BP算法同標準BP算法預(yù)測結(jié)果對比

5 滲透率的預(yù)測

巖石滲透率受孔隙度、粒度、分選、磨圓、礦物組成和黏土分布等諸多因素直接或間接影響［2］。選取孔隙度以及自然伽馬（GR）、電阻率（Rt）、聲波（AC）、密度（DEN）、中子測井（CNL）等作為網(wǎng)絡(luò)的輸入，巖心滲透率為網(wǎng)絡(luò)的輸出。以寧東5井為例，選取7個樣本為訓(xùn)練樣本作為神經(jīng)網(wǎng)絡(luò)的輸入，預(yù)測模型采用3層的神經(jīng)網(wǎng)絡(luò)，其中輸入神經(jīng)元6個，隱含層神經(jīng)元25個，輸出神經(jīng)元1個。GA－BP神經(jīng)網(wǎng)絡(luò)設(shè)置與上文中孔隙度預(yù)測設(shè)置相同，求相似度時，自然伽馬（GR）、電阻率（Rt）、聲波（AC）、密度（DEN）、中子測井（CNL）權(quán)重均取0.15，孔隙度權(quán)重取0.25。訓(xùn)練樣本數(shù)據(jù)和預(yù)測結(jié)果見表4和表5?？梢钥闯?，同孔隙度預(yù)測一樣，隨相似度的減小，GA－BP神經(jīng)網(wǎng)絡(luò)預(yù)測精度降低。

圖5 GA－BP算法同標準BP算法預(yù)測結(jié)果對比圖

表4 訓(xùn)練樣本數(shù)據(jù)

表5 GA－BP神經(jīng)網(wǎng)絡(luò)滲透率預(yù)測數(shù)據(jù)表

6 結(jié) 論

（1）利用BP神經(jīng)網(wǎng)絡(luò)可以建立起測井響應(yīng)和儲層物性參數(shù)之間的非線性映射關(guān)系，將遺傳算法與BP神經(jīng)網(wǎng)絡(luò)相結(jié)合，優(yōu)化了神經(jīng)網(wǎng)絡(luò)的權(quán)值和閥值，避免了標準BP算法易陷入局部極小的缺點，為儲層物性預(yù)測提供了一條很好的途徑。

（2）通過MATLAB工具箱庫函數(shù)的調(diào)用就可以很容易實現(xiàn)GA－BP神經(jīng)網(wǎng)絡(luò)的建模，免去了繁雜的編程，提高了研究效率。

（3）將相似度的概念引入到測井儲層物性預(yù)測中。實例研究表明，相似度越高，GA－BP神經(jīng)網(wǎng)絡(luò)預(yù)測誤差越小。因此，在儲層物性參數(shù)預(yù)測之前，須先計算預(yù)測樣本的相似度，以判斷是否適合，這樣可以有效地控制預(yù)測精度，避免因儲層差別大而造成的神經(jīng)網(wǎng)絡(luò)預(yù)測精度降低的現(xiàn)象。

［1］ 連承波，李漢林，渠芳，等.基于測井資料的BP神經(jīng)網(wǎng)絡(luò)模型在孔隙度定量預(yù)測中的應(yīng)用［J］.天然氣地球科學，2006，17（3）：382－384.

［2］ 陳蓉，王峰.基于MATLAB的BP神經(jīng)網(wǎng)絡(luò)在儲層物性預(yù)測中的應(yīng)用［J］.測井技術(shù)，2009，33（1）：75－78.

［3］ 陳鋼花，董維武.遺傳神經(jīng)網(wǎng)絡(luò)在煤質(zhì)測井評價中的應(yīng)用［J］.測井技術(shù)，2011，35（2）：171－175.

［4］ 林香，姜青山，熊騰科.一種基于遺傳BP神經(jīng)網(wǎng)絡(luò)的預(yù)測模型［J］.計算機研究與發(fā)展，2006，43（suppl.）：338－343.

［5］ 謝立春.BP神經(jīng)網(wǎng)絡(luò)算法的改進及收斂性分析［J］.計算技術(shù)與自動化，2007，26（3）：52－56.

［6］ 趙軍，祁興中，宋帆，等.遺傳算法在測井識別凝析氣藏中的應(yīng)用［J］.測井技術(shù)，2006，30（4）：313－316.

［7］ 洪露，馬長山，謝宗安.基于遺傳算法的神經(jīng)網(wǎng)絡(luò)權(quán)值優(yōu)化［J］.貴州工業(yè)大學學報：自然科學版，2003，32（6）：48－51.

［8］ 王朋巖.用神經(jīng)網(wǎng)絡(luò)預(yù)測儲層的孔隙度［J］.大慶石油學院學報，2003，27（2）：5－7.

Application of Similarity－Genetic Neural Network to Reservoir Parameters Prediction

DONG Xingpeng
（Institute of Seismology，China Earthquake Administration，Wuhan，Hubei 430071，China）

Due to the anisotropy of reservoir，using the linear method it is difficult to get actual reservoir characteristics.On the basis of a combine of genetic algorithm and BP neural network，the global random hunting function of the genetic algorithm is used to optimize neural network connection weights and threshold，which enhances the network tranining precision and parameters prediction accuracy，and as well，increases computing speed by avoiding its disadvantages that standard BP algorithm is apt to trap in local minimal solution.At the same time，we define the similarity in well logging and its calculation formula，and propose a predicting method to combine similarity and Genetic Neural Network.Compared with real samples，the predictive accuracy is higher and effectively controlled，enhance the neural network prediction accuracy.

geophysics logging，genetic algorithm，similarity，neural network，reservoir physical property，prediction

P631

A

2011－12－05 本文編輯 李總南）

1004－1338（2012）03－0267－05

董興朋，男，1988年生，在讀研究生，研究方向為固體地球物理。