粒子群優(yōu)化混合核極限學習機的構(gòu)造煤厚度預測方法

范 君，王 新，徐 慧

(中國礦業(yè)大學計算機科學與技術(shù)學院，江蘇徐州221116)

(*通信作者電子郵箱344298099@qq.com)

0 引言

構(gòu)造煤是指原生結(jié)構(gòu)煤在構(gòu)造應力作用下，煤層物理結(jié)構(gòu)甚至化學成分發(fā)生明顯變化的煤體。早在19世紀20年代就有學者研究構(gòu)造煤的產(chǎn)生和發(fā)育，至今已有研究表明煤礦瓦斯突出區(qū)域與構(gòu)造煤存在著必然的聯(lián)系［1－2］。變形強的構(gòu)造煤儲層的煤層結(jié)構(gòu)及其厚度變化是影響煤與瓦斯突出的主要因素，在受到強烈流變的構(gòu)造煤煤層中，煤層厚度值越大，煤層含氣量越大。理論上，煤層越破碎，煤粒表面能越增加，隨著煤層受到構(gòu)造應力的增強，煤中大分子結(jié)構(gòu)和空隙結(jié)構(gòu)發(fā)生變化，從而改變了甲烷等氣體的吸附能力，就會導致瓦斯壓力增高，增大突出條件。事實也證明，煤體結(jié)構(gòu)受到破壞的地區(qū)是瓦斯含量增高的地方，在高瓦斯礦區(qū)，順煤層滑動構(gòu)造發(fā)育的區(qū)域即是瓦斯突出區(qū)域。我國幅員遼闊，地質(zhì)活動頻繁，含煤地區(qū)大都經(jīng)歷了復雜的地質(zhì)構(gòu)造演化，這為我國煤層構(gòu)造煤提供了發(fā)育的環(huán)境，瓦斯爆炸事故時有發(fā)生。具不完全統(tǒng)計，自從1950年吉林省遼源礦務局富國二礦發(fā)生我國首次有記載的瓦斯突出事故以來，我國煤礦瓦斯突出事故時有發(fā)生，人員傷亡情況嚴重，經(jīng)濟損失慘重。通過定量預測構(gòu)造煤厚度來劃分瓦斯突出區(qū)域，將對煤礦安全管理和煤層氣的開發(fā)與利用起到至關(guān)重要的作用。傳統(tǒng)的構(gòu)造煤厚度的預測方法有煤壁觀測和鉆孔取芯等，其中煤壁觀測法只適用于對已揭露區(qū)的構(gòu)造煤進行編錄，而鉆孔取芯法由于其缺點明顯，例如取芯率低等，也無法保證煤礦的安全生產(chǎn)需求［3］。

當煤層含有構(gòu)造煤時，其彈性性質(zhì)明顯區(qū)別于正常煤層，含構(gòu)造煤煤層的地震屬性也明顯區(qū)別于正常煤層的地震屬性，因此運用地震勘探技術(shù)解決煤礦構(gòu)造煤的預測成為目前解決煤礦地質(zhì)問題的新方法。地震屬性與煤層構(gòu)造煤厚度的關(guān)系具有明顯的非線性特征，而機器學習算法很適合處理這種非線性特征的相關(guān)數(shù)據(jù)［4］。極限學習機(Extreme Learning Machine，ELM)是一種單隱含層前饋神經(jīng)網(wǎng)絡，與其他神經(jīng)網(wǎng)絡算法相比，ELM有明顯的優(yōu)點:具有更快的學習速度;輸入?yún)?shù)更少;可以避免局部極值問題;具有較好的泛化性能。但是由于傳統(tǒng)的ELM算法本身存在一些問題:例如隱含層節(jié)點數(shù)量無法確定、易產(chǎn)生奇點問題等。針對這些問題，將核函數(shù)與傳統(tǒng)的極限學習機相結(jié)合，提出一種核極限學習機模型，直接采用核函數(shù)代替隱含層節(jié)點的顯式映射，無需給定隱含層節(jié)點數(shù)，計算速度快，泛化性能好，得到廣泛應用。

本文在以上研究的基礎上，將多項式核和高斯核相結(jié)合，構(gòu)建混合核極限學習機(Hybrid Kernel Extreme Learning Machine，HKELM)預測模型，利用粒子群算法(Particle Swarm Optimization，PSO)較好的尋優(yōu)能力，優(yōu)化混合核極限學習機的核參數(shù)，并將預測模型應用至實際采區(qū)。在傳統(tǒng)的粒子群算法［5－6］中加入模擬退火的思想、隨進化次數(shù)逐漸減小的慣性權(quán)重以及基于反向?qū)W習變異操作，從而使得粒子群算法更加容易跳出局部最優(yōu)，將地震屬性數(shù)據(jù)運用主成分分析技術(shù)(Principal Component Analysis，PCA)降維消除相關(guān)性，然后輸入構(gòu)建的預測模型中，以提高采區(qū)構(gòu)造煤厚度預測的精度和可靠性。

1 基本原理

1.1 粒子群算法基本原理

由于混合核參數(shù)較多，構(gòu)造混合核極限學習機模型預測構(gòu)造煤厚度關(guān)鍵的一步為優(yōu)選核函數(shù)的參數(shù)。粒子群算法是一種群體智能的優(yōu)化算法。

PSO算法首先在可解空間中初始化一群粒子，算法中每個粒子表示最優(yōu)化問題的一個潛在解，利用位置、速度和適應度三種指標來表示一個已有粒子，其中適應度表示該粒子的優(yōu)劣，其值由適應度函數(shù)計算得到［7］。粒子在解空間中運動，通過個體極值和群體極值更新個體的位置，其中個體極值是指個體所經(jīng)歷位置中計算所得適應度值最優(yōu)位置，群體極值指的是所有粒子適應度最優(yōu)的位置。粒子每更新一次位置就計算一次適應度值，并選取最優(yōu)適應度的粒子更新個體極值和群體極值。

假設在一個D維的搜索空間中有個粒子組成的種群X=(X1，X2，…，Xn)，其中第 i個粒子表示一個 D維向量 Xi=［Xi1，Xi2，…，XiD］T，代表尋優(yōu)問題的一個潛在解。根據(jù)適應度函數(shù)初始化每一個粒子Xi對應的適應度值。設第i個粒子的速度為Vi= ［Vi1，Vi2，…，ViD］T，其個體極值為Pi= ［Pi1，Pi2，…，PiD］T，種群的全局極值為 Pg= ［Pg1，Pg2，…，PgD］T。粒子在每一次迭代過程中更新自身的速度和位置，更新公式如式(1)～(2):

其中:ω為慣性權(quán)重;d=1，2，…，D;i=1，2，…，n;k為當前迭代次數(shù);Vid為當前粒子速度;c1和c2為非負的常數(shù)，稱為加速度因子;r1和r2為分布于［0，1］的隨機數(shù)。為了防止粒子的盲目搜索，通常將粒子的位置和速度限制在［－Xmax，Xmax］、［－ Vmax，Vmax］。

1.2 核函數(shù)

核函數(shù)方法的原理是通過一個特征映射可以將低維輸入空間中的線性不可分數(shù)據(jù)映射到高維特征空間中的線性可分數(shù)據(jù)，其本質(zhì)是對應于高維空間中的內(nèi)積，從而與生成的高維空間的特征映射一一對應。這種對應關(guān)系不僅使待解決的問題變得線性可分，同時又避免了高維空間可能帶來的維數(shù)災難。

設x，z∈X，X屬于R(n)空間，非線性函數(shù)Φ實現(xiàn)從輸入空間X到特征空間F的映射，其中F屬于R(m)，n m。根據(jù)核函數(shù)理論，有:

其中:〈Φ(x)，Φ(z)〉為內(nèi)積;K(x，z)為核函數(shù)。從式(3)可以看出，核函數(shù)將m維高維空間的內(nèi)積運算轉(zhuǎn)化為n維低維輸入空間的核函數(shù)計算，從而解決了在高維特征空間中維數(shù)災難等問題。傳統(tǒng)構(gòu)造核函數(shù)是以Mercer定理為基礎，滿足Mercer條件的函數(shù)均可以作為核函數(shù)。目前使用較多的核函數(shù)主要有四類:

1)高斯核函數(shù):

2)多項式核函數(shù):

3)感知器核函數(shù):

4)線性核函數(shù):

此外，將簡單核函數(shù)構(gòu)造成復雜混合核函數(shù)仍然滿足Mercer定理對核函數(shù)的要求。

1.3 HKELM的基本原理

傳統(tǒng)的ELM是一個三層前饋神經(jīng)網(wǎng)絡，在這個網(wǎng)絡中輸入權(quán)值和隱含層閾值是隨機給的，影響模型性能的輸入?yún)?shù)是隱含層節(jié)點數(shù)，而整個模型需要訓練確定的是隱含層權(quán)值，其計算過程是通過計算一個矩陣的偽逆來實現(xiàn)的［8－9］。

若采用(a，b)表示輸入權(quán)值和隱含層閾值，樣本的訓練集用(x，t)表示，隱含層映射函數(shù)為g(x)，輸出權(quán)重為β，隱含層節(jié)點數(shù)為L，則極限學習機模型可表示為式(8):

令隱含層輸出矩陣為H，則式(8)可轉(zhuǎn)化為O=‖Hβ－T‖。極限學習機不斷學習的過程中，誤差O不斷減小，當無誤差學習時Hβ=T。隱含層節(jié)點數(shù)和樣本數(shù)相等時，H為非奇異矩陣，β可由H的逆矩陣求出，但一般隱含層節(jié)點數(shù)要遠小于樣本數(shù)，H不存在逆矩陣。此時廣義逆矩陣可以用于求解奇異矩陣的逆:

其中H+表示隱含層輸出矩陣的廣義逆矩陣［10］。

將核學習理論加入ELM模型中，隱含層將之前傳統(tǒng)激活函數(shù)顯式映射的任務交給核函數(shù)處理［11］。核函數(shù)的類型可以分為全局型核函數(shù)和局部型核函數(shù)，例如常用的高斯徑向基核函數(shù)為局部型核函數(shù)，在測試點附近數(shù)據(jù)點對核函數(shù)有影響，學習能力較強但泛化能力較弱;感知器核和多項式核為全局型核函數(shù)，遠離測試點的數(shù)據(jù)對核函數(shù)也有影響，泛化能力較強但學習能力較弱。

因此，本文將從兩類核函數(shù)中分別選取一種核函數(shù)組合成混合核，該混合核函數(shù)將作為混合核極限學習機模型的核函數(shù)［12－13］。

2 模型優(yōu)化

2.1 粒子群算法的改進

針對HKELM模型參數(shù)較多的問題，本文提出用粒子群算法進行參數(shù)尋優(yōu)［14］。為了提高粒子群算法尋優(yōu)能力，針對粒子群算法容易陷入局部最小的問題，對傳統(tǒng)的粒子群算法主要作出三點改進:首先，利用模擬退火算法具有較強全局收斂性的特點，將模擬退火算法的思想引入粒子群算法，在加溫、等溫、冷卻的過程中增強粒子的多樣性，使得粒子群算法更有機會跳出局部最優(yōu)的陷阱。其次，由于傳統(tǒng)的PSO算法中，慣性權(quán)重ω表示繼承先前粒子速度的能力，當ω較大時有利于算法的全局搜索能力，反之則有利于算法的局部搜索，ω值設置過大或過小容易造成PSO算法陷入局部最優(yōu)或是無法真正地找出最優(yōu)解。因此，加入隨迭代次數(shù)減小的慣性權(quán)重，慣性權(quán)重如式(10)所示:

其中:ωs為初始慣性權(quán)重;ωe為迭代至最大次數(shù)時的慣性權(quán)重;k為當前迭代代數(shù);Tmax為最大迭代代數(shù)。ω隨迭代次數(shù)的增加而減小，既保證了迭代前期較強的全局搜索能力，又保證了在迭代后期粒子的局部搜索能力。最后，在粒子群算法中引入基于隨機的反向?qū)W習變異操作。變異概率如式(11):

其中:P 為變異概率，且當P ＞ 0.3時，P=0.5;Ps和Pe為常數(shù)，實驗中取 Ps=0.35，Pe=0.1。變異概率隨迭代次數(shù)先增大后減小，從而在一定程度上保留優(yōu)秀粒子的同時增強粒子多樣性。變異公式為基于隨機的反向?qū)W習公式，如式(12)所示:

其中z為全局最佳粒子。

2.2 混合核極限學習機

極限學習機的核函數(shù)按照對數(shù)據(jù)的影響可以分為兩大類，局部型核函數(shù)和全局型核函數(shù)。局部型核函數(shù)學習能力強，但泛化性能相對弱;全局型核函數(shù)學習能力一般，但泛化性能強［15］。將局部核和全局核混合成為的混合核函數(shù)具備較強的學習能力，同時又有良好的泛化能力。

通常將局部的高斯徑向基核函數(shù)與全局的多項式核函數(shù)混合生成混合核構(gòu)建的HKELM模型既具有較強的學習能力，又具有不錯的泛化性能，因此本文在此基礎上構(gòu)建預測構(gòu)造煤厚度的預測模型。高斯徑向基核函數(shù)具有較強的學習能力，能夠很容易地將訓練樣本在特征空間中線性分開，并能夠?qū)ο嗑嘁欢ǚ秶鷥?nèi)的樣本準確預測，但對超過一定范圍的樣本則無法準確預測。

如圖1 所示，以X=0.4為測試點，分別取σ為0.1、0.2、0.3、0.4。徑向基核函數(shù)具有局部特性，體現(xiàn)在它只對測試點附近的值產(chǎn)生影響，且離測試點越近的點影響越大，由此可知，局部性核函數(shù)學習能力較強，泛化能力較弱。

圖1 徑向基核函數(shù)曲線Fig.1 Radial basis kernel function curve

多項式核函數(shù)具有較強的泛化能力，對即使相距很遠的樣本點也能產(chǎn)生影響，但對局部的樣本能力較弱［16］。

如圖2所示，以X=0.4為測試點，m=1，n=1，分別取d為1、2、3、4。多項式核函數(shù)的全局性特性體現(xiàn)在與測試點相距較遠的樣本點也能對函數(shù)值產(chǎn)生一定的影響，具有較強的泛化能力，學習能力較弱［16］。

圖2 多項式核函數(shù)曲線Fig.2 Polynomial kernel function curve

將多項式核函數(shù)與徑向基核函數(shù)結(jié)合，得到混合核函數(shù):

如圖3 所示，以 X=0.4 為測試點，分別取 λ 為0.2、0.3、0.4、0.5，σ取值為0.1，d為2。由圖3可以看出，混合核函數(shù)不僅對測試點附近的樣本具有一定的影響，并且當樣本點遠離測試點時，仍能產(chǎn)生影響。因此，混合核函數(shù)結(jié)合了多項式核函數(shù)和高斯徑向基核函數(shù)的特點，有效地彌補了單核核函數(shù)的缺點，從而提高了擬合能力。

此外，為了防止模型復雜度過高引起過擬合現(xiàn)象，在核函數(shù)的基礎上加入L2正則項，有效地避免了噪聲和異常點對模型泛化性能的影響。

圖3 混合核函數(shù)曲線Fig.3 Hybrid kernel function curve

2.3 算法描述

改進PSO優(yōu)化混合核極限學習機參數(shù)的算法具體過程如下:1)數(shù)據(jù)標準化，進行主成分分析處理，保存結(jié)果集備用;2)將σ、m、n、λ設為粒子(d=2)，隨機初始化粒子的位置和速度;

3)根據(jù)式(13)計算訓練集隱含層節(jié)點的輸出，并加入L2正則項;

4)計算隱含層輸出權(quán)值矩陣;

5)計算驗證集隱含層節(jié)點的輸出，并根據(jù)隱含層節(jié)點輸出計算得到測試集的預測值;

6)將均方誤差mse作為粒子適應度，計算個體極值和種群極值;

7)根據(jù)更新公式更新粒子位置和速度;

8)根據(jù)變異概率進行粒子變異;

9)計算更新后粒子的適應度，更新個體極值和種群極值，若未滿足最大迭代次數(shù)，則返回7)，否則繼續(xù)10);

10)保存種群最優(yōu)適應度對應的粒子，即最佳混合核函數(shù)參數(shù)，將參數(shù)代入HKELM模型，獲得預測構(gòu)造煤厚度的預測模型，并用測試集進行測試。

由于混合核極限學習機參數(shù)較多，而用傳統(tǒng)的粒子群算法優(yōu)化混合核極限學習機無法找出最優(yōu)參數(shù)下的預測模型，加入隨迭代次數(shù)逐漸減小的慣性權(quán)重、模擬退火思想以及基于反向?qū)W習的變異公式使得PSO在尋優(yōu)能力和收斂速度上都有了改善。相比傳統(tǒng)ELM激活函數(shù)的顯示映射，混合核函數(shù)在較優(yōu)的參數(shù)下，兼顧了數(shù)據(jù)的局部性和全局性，既具有較強的學習能力，又具有較好的泛化性能。因此，通過改進PSO尋找到最優(yōu)的核函數(shù)參數(shù)，獲得最優(yōu)核參數(shù)下的HKELM模型，理論上可以作為預測構(gòu)造煤厚度的預測模型。

2.4 預測模型建立

2.4.1 測試數(shù)據(jù)

為了建立適合構(gòu)造煤厚度預測的HKELM模型，首先建立了構(gòu)造煤模型如圖4所示，煤層的頂?shù)装鍏?shù)為最薄處為0 m，最厚處為10 m;原生煤最薄處為0 m，最厚處也為10 m;直接頂/底是厚度為2 m的泥巖，老頂/底為砂巖。根據(jù)褶積模型，利用主頻為50 Hz的Ricker子波進行正演模擬，獲得其對應的正演地震剖面。利用圖4所示的構(gòu)造煤模型正演地震剖面如圖5所示，分別提取出剖面中的正負相位的地震屬性。

利用圖4、圖5建立的模型，以所提取出的負相位地震屬性為例，分別進行有噪聲與無噪聲數(shù)據(jù)進行測試。地震屬性維數(shù)較大，為去除屬性中的相關(guān)屬性，實驗使用PCA對地震屬性進行降維去噪，取累計貢獻率大于95%的主成分并求其主成分得分作為實驗數(shù)據(jù)進行后續(xù)的實驗分析。

圖4 構(gòu)造煤模型Fig.4 Tectonic coal model

2.4.2 參數(shù)設置

在運用改進的粒子群算法優(yōu)化混合核極限學習機構(gòu)建的構(gòu)造煤厚度預測模型中，學習因子c1、c2均設置為2，最大迭代數(shù)MAXGEN=300，種群粒子數(shù)SIZEPOP=20，模擬退火算法初始溫度T=100，退火常數(shù)L=0.25，初始慣性權(quán)重ωs為1.2，迭代至最大次數(shù)時的慣性權(quán)重ωe為0.4，核極限學習機中核函數(shù)的參數(shù)由改進粒子群算法尋優(yōu)得到。此外，為了探究本文提出改進的模型預測性能，將改進后的模型與傳統(tǒng)的BP神經(jīng)網(wǎng)絡和用K折交叉驗證方法優(yōu)化參數(shù)的支持向量機(Support Vector Machines，SVM)進行預測結(jié)果對比。其中BP神經(jīng)網(wǎng)絡采用Matlab自帶函數(shù)newff()建立模型，輸入層神經(jīng)元為9，輸出層神經(jīng)元設為1，隱含層神經(jīng)元數(shù)設為10，輸入層和隱含層、隱含層和輸出層之間的傳遞函數(shù)分別選擇tansig和logsig函數(shù)，最大訓練次數(shù)設為500，學習速率為0.1，訓練算法設置為 LM算法，動量因子為0.9，期望誤差設置為10－3，性能函數(shù)采用均方誤差函數(shù)(Mean Squared Error，MSE)和決定系數(shù)(Coefficient of Correlation)記為R2。利用Libsvm工具箱建立SVM預測模型，交叉驗證數(shù)的參數(shù)v設置為5，即參數(shù)c和參數(shù)g由5折交叉驗證得出，其中搜索空間為［－30，30］，步長設為0.5，核函數(shù)類型參數(shù)t設置為2，即采用徑向基函數(shù)作為核函數(shù)，支持向量機的類型參數(shù)s設置為3，即使用εSVR類型，即ε-支持向量機回歸算法，εSVR中損失值 ε 設置為0.1。

2.4.3 測試結(jié)果

無噪聲數(shù)據(jù)提取出26個地震屬性，通過PCA屬性降維獲得7個線性不相關(guān)屬性。含噪聲數(shù)據(jù)獲得11個地震屬性，進行屬性降維獲得9個線性不相關(guān)屬性。兩種數(shù)據(jù)集分別取98個樣本進行實驗并隨機抽取其70%作為訓練集，剩余30%作為測試集測試改進PSO-HKELM模型。預測結(jié)果如圖6所示。

圖6預測結(jié)果顯示，在無噪聲數(shù)據(jù)集下，改進的粒子群算法優(yōu)化混合核極限學習機構(gòu)建的構(gòu)造煤厚度預測模型有很好的擬合效果，擬合度可以達到0.989，均方誤差為0.0041971;而含噪聲數(shù)據(jù)集的預測結(jié)果比無噪聲數(shù)據(jù)集的預測結(jié)果要差，擬合度為0.9293，均方誤差為0.0195，但在噪聲數(shù)據(jù)的影響下，同樣體現(xiàn)了較好的擬合效果。

圖6 模擬不同數(shù)據(jù)預測結(jié)果Fig.6 Prediction results of simulating different data

目前國內(nèi)外已有利用智能算法預測構(gòu)造煤厚度的研究［17］，為進一步驗證模型的優(yōu)劣，在當前學者研究的基礎上，將改進的粒子群算法優(yōu)化混合核極限學習機構(gòu)建的構(gòu)造煤厚度預測模型分別與用K折交叉驗證方法優(yōu)化參數(shù)的SVM模型以及BP神經(jīng)網(wǎng)絡進行實驗比較［18］。取98個含噪聲數(shù)據(jù)集樣本進行實驗構(gòu)建模型，每個模型分別運行10次，且每次均隨機抽取其70%作為訓練集，剩余30%作為測試集，預測結(jié)果如圖7所示。

圖7 不同算法預測結(jié)果比較Fig.7 Comparison of prediction results of different algorithms

由圖7可以看出，采用的改進PSO-HKELM模型的預測準確率明顯高于另兩個模型預測準確率。在10次實驗中，改進PSO-HKELM模型的平均預測準確度達到90%以上。此外，將運行時間作為判斷預測模型時間復雜度的標準，模擬含噪聲數(shù)據(jù)在相同環(huán)境和參數(shù)下運行5次，運行時間如表1所示(保留兩位小數(shù))。

表1 模擬含噪聲數(shù)據(jù)的運行時間 sTab.1 Running time of simulated data with noise s

由表1可以發(fā)現(xiàn)，改進模型運行時間較長，算法復雜度較高，但考慮到其預測精度提高較為明顯，預測結(jié)果較為穩(wěn)定，且在實際開采中，預測構(gòu)造煤厚度對實時要求不高，因此，本文提出的一種改進PSO-HKELM預測模型在預測構(gòu)造煤厚度上的運用是可行的，且預測效果較好。

3 實例預測

3.1 研究區(qū)概況

本文選取陽煤集團新景煤礦蘆南二采區(qū)中部的15#煤層作為研究區(qū)域，結(jié)合實際地震屬性，對構(gòu)造煤厚度進行預測。陽煤集團新景煤礦底層層滑構(gòu)造特征明顯，瓦斯含量較高，極易發(fā)生煤與瓦斯突出，因此陽煤集團新景煤礦成為研究煤與瓦斯突出課題的理想?yún)^(qū)域。新景煤礦二采區(qū)面積為3.42 km2，其中15#煤層埋深 575.9 m ～640.93 m，平均深埋623.8 m，煤層厚度為6.25 m ～10.4 m，平均厚度 7.925 m;構(gòu)造煤厚度為0 m～4.3 m，平均厚度為1.57 m。鉆孔共有10口，詳細情況見表2。

表2 15#煤層鉆孔揭露煤層信息 mTab.2 Exposing coal information of drills in 15#coal seam m

根據(jù)鉆孔揭露的煤層厚度和構(gòu)造煤厚度，插值生成煤層厚度分布圖，如圖8所示。

3.2 構(gòu)造煤厚度的預測

利用第2章建立的預測模型，對研究目標區(qū)域的構(gòu)造煤厚度進行預測。如果只采取已知的10口鉆孔數(shù)據(jù)作為訓練樣本，則訓練樣本太少，預測模型的可靠性難以保證。因此，將鉆孔附近較小范圍內(nèi)(15 m×15 m)的地震屬性提取出作為訓練集(共2232個樣本)，以提高訓練集的數(shù)量。15#煤層收集地震屬性共14維，首先對14維數(shù)據(jù)進行主成分分析，降維后的數(shù)據(jù)變?yōu)?維。然后采用2232×6的數(shù)據(jù)訓練構(gòu)造煤厚度預測模型。將BP、SVM和改進PSO-HKELM模型分別運用到整個15#煤層的構(gòu)造煤厚度預測中，將三種方法預測出的10口鉆孔處的構(gòu)造煤厚度提取出來，與實際的鉆孔數(shù)據(jù)作對比。10口鉆孔的值與預測誤差如表3所示。

圖8 15#煤層厚度分布Fig.8 15#coal seam thickness distribution

表3 不同方法鉆孔處預測值對比 mTab.3 Prediction value comparison of drilling holes of different methods m

由表3可以得出，其中改進PSO-HKELM模型預測出的相對誤差最小，誤差均值為0.05(其中最大誤差為0.06，最小誤差為0.00)，而且預測結(jié)果均比較穩(wěn)定;K折交叉驗證方法優(yōu)化參數(shù)的SVM的預測精度其次，誤差均值為0.21(其中最大誤差為0.97，最小誤差為0.00);BP效果最弱，誤差均值為0.27(其中最大誤差為 1.02，最小誤差為 0.01)。因此，利用本文提出的改進PSO-HKELM模型在構(gòu)造煤厚度的預測上效果更好。

數(shù)據(jù)經(jīng)過PCA處理，運用改進PSO-HKELM模型預測出的構(gòu)造煤厚度分布圖如圖9所示。

圖9 預測目標區(qū)構(gòu)造煤厚度分布Fig.9 Distribution of tectonic coal thickness in prediction target area

根據(jù)圖9總體上可以看出，預測分布圖與煤層整體的構(gòu)造煤分布具有很高的一致性。構(gòu)造煤主要分布在井3-1381、井3-146和井3-137為中心的三個獨立區(qū)域。

4 結(jié)語

本文利用地震屬性預測構(gòu)造煤厚度，首先運用主成分分析的方法對數(shù)據(jù)進行預處理，然后針對傳統(tǒng)預測方法的不足，建立了改進的PSO-HKELM的構(gòu)造煤厚度預測的模型。采用較為常用的高斯核和多項式核組合成為混合核極限學習機的核函數(shù)，并用改進的粒子群算法對混合核參數(shù)進行優(yōu)化。在改進的粒子群算法中加入模擬退火的思想、隨迭代次數(shù)減小的慣性權(quán)重以及基于隨機反向?qū)W習的變異操作，能有效地彌補粒子群算法較易陷入局部最優(yōu)的缺點。最后，將建立預測模型應用于陽煤集團新景煤礦#15采區(qū)的實際煤層。通過模擬數(shù)據(jù)和新景煤礦#15區(qū)煤層數(shù)據(jù)的實驗分析，采用改進的PSO-HKELM算法對構(gòu)造煤厚度建立預測模型，其性能將優(yōu)于傳統(tǒng)的BP神經(jīng)網(wǎng)絡和用K折交叉驗證優(yōu)化參數(shù)SVM預測模型，且其精度較高，為預測煤礦采區(qū)煤與瓦斯突出形勢提供了一種可行的新途徑。構(gòu)建的混合核極限學習機的構(gòu)造煤厚度預測模型可以較為準確地預測實際采區(qū)構(gòu)造煤厚度，但預測結(jié)果受地震屬性信噪比的影響較為嚴重，當訓練數(shù)據(jù)中含有噪聲的樣本較多時，預測結(jié)果不理想，因此尋找一種降低數(shù)據(jù)信噪比的處理方法需要進一步的研究。