杜潤(rùn)林，劉 展

1.中國(guó)石油大學(xué)(華東)地球科學(xué)與技術(shù)學(xué)院， 山東 青島 266555 2.青島海洋地質(zhì)研究所國(guó)土資源部海洋油氣資源與環(huán)境地質(zhì)重點(diǎn)實(shí)驗(yàn)室， 山東 青島 266071

?

基于粒子群優(yōu)化的細(xì)胞神經(jīng)網(wǎng)絡(luò)油氣重力異常信息提取

杜潤(rùn)林1，2，劉 展1

1.中國(guó)石油大學(xué)(華東)地球科學(xué)與技術(shù)學(xué)院， 山東 青島 266555 2.青島海洋地質(zhì)研究所國(guó)土資源部海洋油氣資源與環(huán)境地質(zhì)重點(diǎn)實(shí)驗(yàn)室， 山東 青島 266071

細(xì)胞神經(jīng)網(wǎng)絡(luò)方法是一種有效的重力異常提取方法，提取出的重力異常信息具有較強(qiáng)的橫向分辨能力，但如何訓(xùn)練優(yōu)化模板是影響該方法提取精度的關(guān)鍵所在。筆者引入粒子群優(yōu)化算法進(jìn)行參數(shù)優(yōu)化，結(jié)合細(xì)胞神經(jīng)網(wǎng)絡(luò)和粒子群方法各自的特點(diǎn)，通過(guò)細(xì)胞神經(jīng)網(wǎng)絡(luò)動(dòng)態(tài)性能分析，對(duì)模板加入約束條件，并在粒子群算法中引入收縮因子，提出一種基于改進(jìn)的粒子群優(yōu)化的細(xì)胞神經(jīng)網(wǎng)絡(luò)算法。使用該方法進(jìn)行了模型驗(yàn)證和實(shí)際資料處理，處理結(jié)果表明，該方法提高了重力異常提取的穩(wěn)定性，能夠更準(zhǔn)確地提取油氣藏重力異常信息。

細(xì)胞神經(jīng)網(wǎng)絡(luò)；粒子群優(yōu)化；重力異常；收縮因子;油氣

0 引言

油氣儲(chǔ)層往往比圍巖密度小，會(huì)在其上方產(chǎn)生重力低異常值，但油氣產(chǎn)生的重力異常普遍很微弱，而且被背景構(gòu)造產(chǎn)生的重力異常掩蓋，難以識(shí)別。國(guó)內(nèi)外學(xué)者利用不同方法從不同角度進(jìn)行了油氣重力異常的提取研究，如歸一化總梯度方法[1]、插值切割法[2]、小波提取法[3-4]等，可是這些方法對(duì)于微弱的油氣異常提取效果并不好。

細(xì)胞神經(jīng)網(wǎng)絡(luò)(cellular neural network，CNN)是L. O. Chua在1988年提出的[5]。Muhittin等[6-7]首次將CNN用于重力異常的分離，說(shuō)明該方法能夠突出淺層異常、壓制區(qū)域異常。劉展等[8]使用擬BP算法改進(jìn)了CNN方法，對(duì)花溝地區(qū)火成巖重力異常進(jìn)行提取，并且總結(jié)了CNN方法的適用條件。但以上方法存在參數(shù)設(shè)置比較復(fù)雜，容易陷入局部極小解等缺點(diǎn)。針對(duì)以上不足，筆者將粒子群算法應(yīng)用到CNN模板設(shè)計(jì)中。粒子群優(yōu)化(particle swarm optimization, PSO)算法是由Jim Kennedy 于1995年提出的隨機(jī)全局優(yōu)化方法[9]，同其他算法相比較，其優(yōu)勢(shì)為沒(méi)有過(guò)多參數(shù)需要調(diào)整，能比較快速地找到最優(yōu)解[10-15]。在本文中，筆者結(jié)合CNN動(dòng)態(tài)性能分析，使用粒子群算法優(yōu)化設(shè)計(jì)CNN的模板參數(shù)，并對(duì)油氣重力異常信息進(jìn)行提取，用粒子群算法對(duì)CNN模板參數(shù)優(yōu)化并進(jìn)行油氣重力異常信息提取，通過(guò)對(duì)比分析說(shuō)明得到參數(shù)的合理性和有效性。

1 細(xì)胞神經(jīng)網(wǎng)絡(luò)原理

細(xì)胞神經(jīng)網(wǎng)絡(luò)中的一個(gè)細(xì)胞僅與它的近鄰細(xì)胞相連接，近鄰細(xì)胞可以彼此直接相互作用，因?yàn)樯窠?jīng)網(wǎng)絡(luò)連續(xù)時(shí)間動(dòng)力學(xué)的傳播效應(yīng)，非近鄰細(xì)胞不直接連接[16]。細(xì)胞神經(jīng)網(wǎng)絡(luò)中細(xì)胞的鄰域C(i,j)定義為

(1)

式中：r為半徑；M，N分別為橫坐標(biāo)和縱坐標(biāo)方向的搜索最大值。則狀態(tài)方程為

(2)

其中：xij為狀態(tài)變量；yij為輸出變量；uij為輸入變量；I為偏量；A為反饋系數(shù)矩陣；B為控制吸收系數(shù)矩陣；t為模型參數(shù)。輸出方程為

(3)

其中，f(x)為分段線性函數(shù)，通過(guò)對(duì)A、B、I的設(shè)計(jì)優(yōu)化實(shí)現(xiàn)油氣重力異常信息提取。

2 粒子群算法原理

PSO中，每個(gè)優(yōu)化問(wèn)題的解都是搜索空間的一個(gè)“粒子”的狀態(tài)，所有粒子都由一個(gè)適應(yīng)函數(shù)(fitness function)決定適應(yīng)值(fitness value)，每個(gè)粒子還有一個(gè)速度直接影響它們的飛翔方向和距離，粒子根據(jù)當(dāng)前自身情況和粒子群情況在解空間中搜索[17]。假設(shè)在D維搜索空間中，有m個(gè)粒子，其中第i個(gè)粒子的位置和速度分別是xi=(xi1,xi2,…,xiD)，vi=(vi1,vi2,…,viD)，i=1，2，…，m，將vi帶入目標(biāo)函數(shù)可計(jì)算其適應(yīng)值。記第i個(gè)粒子搜索到的最優(yōu)位置為piBest=(pi1,pi2,…,piD)，整個(gè)粒子群搜索到的最優(yōu)位置為pgBest=(pg1,pg2,…,pgD)，則根據(jù)公式(4)和(5)可更新粒子的速度和位置：

(4)

(5)

其中：w是非負(fù)常數(shù)，稱為慣性因子，w也可以隨著迭代線性地減少；學(xué)習(xí)因子c1和c2是非負(fù)常數(shù)；r1和r2是[0，1]之間的隨機(jī)數(shù)；vi∈[vmin,vmax]，vmin和vmax由人為設(shè)定。迭代中止條件為最大迭代次數(shù)或粒子群迄今為止搜索到的最優(yōu)位置滿足閾值。

3 基于粒子群算法的細(xì)胞神經(jīng)網(wǎng)絡(luò)模板學(xué)習(xí)方法

3.1 種群參數(shù)的確定

由于細(xì)胞神經(jīng)網(wǎng)絡(luò)具有空間不變性，因此進(jìn)行油氣重力異常信息提取時(shí)，采用如下模板格式：

(6)

由模板可知，需要用PSO算法對(duì)x1，x2，x3，x4，x5五個(gè)參數(shù)進(jìn)行優(yōu)化。

3.2 適應(yīng)度函數(shù)

將細(xì)胞神經(jīng)網(wǎng)絡(luò)輸出端誤差E(ω)的平方根作為PSO算法的適應(yīng)度函數(shù)[18-20]：

(7)

式中：yij是指神經(jīng)網(wǎng)絡(luò)的輸出；ωij是期望輸出。

3.3 帶約束條件的粒子群算法(C-PSO)優(yōu)化CNN模板參數(shù)

通過(guò)對(duì)CNN的動(dòng)態(tài)性能分析[21-24]，得到用于油氣重力異常信息提取的五個(gè)模板參數(shù)有如下關(guān)系：

(8)

將式(8)作為一種約束對(duì)CNN模板進(jìn)行優(yōu)化，可以使種群較快地收斂到最優(yōu)值。此外，為了提高收斂速度，對(duì)式(4)引入收縮因子λ：

(9)

3.4C-PSO算法的實(shí)現(xiàn)

基于粒子群優(yōu)化細(xì)胞神經(jīng)網(wǎng)絡(luò)模板的流程圖如圖1所示。

圖1 基于粒子群優(yōu)化細(xì)胞神經(jīng)網(wǎng)絡(luò)模板的流程圖Fig.1 Flow chart of templates of CNN based on PSO

4 模型試算和實(shí)際資料處理

4.1 模板訓(xùn)練

為了驗(yàn)證方法的有效性，首先建立三球體模型：用一個(gè)放置在深部的球體(球1)模擬區(qū)域異常，在淺層放置兩個(gè)較小的球體(球2、球3)模擬局部異常，參數(shù)如表1所示。此三球體模型的疊加重力異常如圖2a所示，圖2b為兩個(gè)淺層球體(球2、球3)的疊加重力異常。將淺部小球球2和球3的疊加重力異常作為細(xì)胞神經(jīng)網(wǎng)絡(luò)的期望輸出，使用擬BP、PSO、C-PSO三種算法分別對(duì)模板參數(shù)進(jìn)行優(yōu)化。表2為三種算法得到的提取重力異常信息的模板，表3為三種算法的迭代次數(shù)和時(shí)間。從表3可以看出，C-PSO算法尋找最優(yōu)值的迭代次數(shù)和迭代時(shí)間都明顯小于擬BP和PSO算法。

4.2 疊加模型試驗(yàn)

建立四球體模型，參數(shù)如表4所示。使用表1中的三個(gè)模板對(duì)模型進(jìn)行重力異常信息提取，結(jié)果如圖3所示。從圖3可以看出：C-PSO算法優(yōu)化模板提取的重力異常信息明顯要比擬BP算法的準(zhǔn)確，抗噪性更好；而與PSO算法對(duì)比，兩者得到的結(jié)果幾乎一致。但由4.1節(jié)知C-PSO算法要比PSO算法快。綜上所述，C-PSO算法在優(yōu)化提取重力異常信息CNN模板時(shí)要優(yōu)于其他兩種算法。

表1 三球體模型參數(shù)

4.3 油氣藏模型試算

建立背斜油氣藏模型，縱向和橫向長(zhǎng)度均為100 km，油氣儲(chǔ)層位于第三層，剩余密度為-0.10 g/cm3，從上向下各層剩余密度依次為0.50、0.75、1.00、1.50 g/cm3。圖4為在y=50 km處、平行于橫軸的模型橫切面圖。

圖5a為油氣藏模型產(chǎn)生的重力異常。由于區(qū)域異常的掩蓋，在圖5a中看不出油氣異常信息。使用表2中C-PSO算法得到的模板進(jìn)行油氣異常信息提取，結(jié)果如圖5b所示。從圖5b中可以看出，背斜構(gòu)造引起的區(qū)域異常被消除，提取的相對(duì)重力負(fù)異常對(duì)應(yīng)儲(chǔ)層位置。說(shuō)明該方法能夠提取油氣異常信息，消除區(qū)域異常干擾。

圖2 三個(gè)(a)和兩個(gè)(b)淺層球體重力疊加異常Fig.2 Gravity anomaly of three spheres (a) and two shallow spheres (b)

擬BP算法A=00．078200．07820．48390．078200．07820?è???，B=-0．0893-0．0893-0．0893-0．08930．8438-0．0893-0．0893-0．0893-0．0893?è???，I=0．0325PSO算法A=00．025700．02570．41710．025700．02570?è???，B=-0．0638-0．0638-0．0638-0．06380．6152-0．0638-0．0638-0．0638-0．0638?è???，I=0．0472CPSO算法A=00．085700．08570．52710．085700．08570?è???，B=-0．0738-0．0738-0．0738-0．07380．7152-0．0738-0．0738-0．0738-0．0738?è???，I=0．0472

表3 模板訓(xùn)練的迭代次數(shù)和迭代時(shí)間

Table 3 Number of iterations and iteration time of template training

迭代次數(shù)迭代時(shí)間/s擬BP算法7333PSO算法3211CPSO算法177

表4 四球體模型參數(shù)

a. 四球模型疊加重力異常；b. 擬BP算法處理結(jié)果；c. PSO算法處理結(jié)果；d. C-PSO算法處理結(jié)果。圖3 四球體模型原始重力異常及三種方法提取的異常Fig.3 Original gravity anomaly of four spheres and the anomaly extracted by three methods

圖4 油氣藏模型橫切面圖Fig.4 Cross-sectional diagram of reservoir model

4.4 實(shí)際資料處理

筆者利用C-PSO算法優(yōu)化的CNN模板對(duì)柴達(dá)木盆地三湖地區(qū)澀北實(shí)測(cè)重力資料進(jìn)行油氣重力異常信息提取。該地區(qū)天然氣主要產(chǎn)自第四系湖相沉積地層，氣田區(qū)地層密度比周圍同一套地層密度低0.10 g/cm3左右，是典型的背斜油氣藏。因此，使用本文4.3節(jié)中的模型進(jìn)行模板訓(xùn)練，得到的模板為

(10)

圖6a為澀北的實(shí)測(cè)重力異常，氣田區(qū)地層密度比周圍地層密度低，因此在其上方產(chǎn)生相對(duì)重力低；但由于構(gòu)造等因素的影響，并不能確定氣田范圍。使用模板(10)對(duì)澀北勘探區(qū)進(jìn)行油氣重力異常信息提取，結(jié)果如圖6b所示。從圖6b中可以看出提取的油氣重力異常信息與已知的工業(yè)氣井有很好的對(duì)應(yīng)關(guān)系；說(shuō)明該方法能夠壓制區(qū)域異常，提取油氣異常信息，圈定油氣邊界，并且有較強(qiáng)的抗干擾能力。

圖5 原始重力異常(a)及C-PSO算法提取的異常(b)Fig.5 Original gravity anomaly (a) and the anomaly extracted C-PSO algorithm (b)

Δg為重力異常。圖6 澀北原始重力異常(a)及C-PSO算法處理結(jié)果(b)Fig.6 Gravity anomaly of Sebei (a) and reservoir gravity anomaly extracted by C-PSO (b)

5 結(jié)語(yǔ)

本文提出的C-PSO算法優(yōu)化細(xì)胞神經(jīng)網(wǎng)絡(luò)模板在提取重力異常信息時(shí)，無(wú)論是在迭代步數(shù)、時(shí)間、搜索精度上均優(yōu)于擬BP和PSO算法，該算法能夠根據(jù)自己的期望訓(xùn)練模板，更精準(zhǔn)地提取油氣重力異常信息，消除或減弱其他干擾。

需要指出的是，在提取油氣重力異常信息時(shí)，該方法需要大量初始模型訓(xùn)練模板，因此要足夠多的先驗(yàn)信息才能夠保證結(jié)果的準(zhǔn)確性。后續(xù)工作中，我們將重點(diǎn)研究如何在保證精度的前提下，減小方法對(duì)先驗(yàn)信息的依賴。

[1] 張鳳旭，孟令順，張鳳琴，等.利用Hilbert變換計(jì)算重力歸一化總梯度[J].地球物理學(xué)報(bào)，2005，48(3):704-709. Zhang Fengxu, Meng Lingshun, Zhang Fengqin, et al. Calculating Normalized Full Gradient of Gravity Anomaly Using Hilbert Ttransform[J]. Chinese Journal of Geophysics, 2005,48(3)：704-709.

[2] 張剛. 油氣重磁異常識(shí)別及提取方法研究[D].青島: 中國(guó)石油大學(xué), 2011. Zhang Gang. The Method Rresearch of Identification and Extraction of Oil Gravity and Magnetic Anomalies[D]. Qingdao: China University of Petroleum, 2011.

[3] 耿喜哲，張春華. 綜合小波分析與分形方法進(jìn)行弱磁異常識(shí)別[J].物探與化探, 2002, 26(2): 126-130. Geng Xizhe, Zhang Chunhua. The Application of Wavelet Analysis and Fractal Technique to the Recoganion of Micro-Magnetic Anomaly[J]. Geophysical & Geochemical Exploration, 2002, 26(2): 126-130.

[4] 宋雙, 劉天佑, 蘇俊青. 重磁小波多尺度分解及在火山巖識(shí)別中的應(yīng)用[J]. 工程地球物理學(xué)報(bào), 2008, 5(6): 675-679. Song Shuang, Liu Tianyou, Su Junqing. The Wavelet Multi-Scale Decomposition of Gravity and Magnet and Its Application in the Identification of Volcanic Rock[J]. Chinese Journal of Engineering Geophysics, 2008, 5(6): 675-679.

[5] Chua L O, Yang L. Cellular Neural Networks：Theory[J]. IEEE Trans on Circuits, 1988,35(5): 1273-1290.

[6] Albora A M，Ucan O N，Ozmen A，et al. Septaration of Bouguer Anomaly Map Using Cellular Neural Network[J]. Journal Applied Geophysics, 2001, 46:129-142.

[7] Albora A M，Ozmen A. Wavelet Based Design of Dynamical Neural Networks for Magnetic Map Processing[R].Istanbul:University of Istanbul Press,2002:133-140.

[8] 劉展，劉茂成，魏巍，等．重力異常分離的細(xì)胞神經(jīng)網(wǎng)絡(luò)方法研究[J]．中國(guó)石油大學(xué)學(xué)報(bào):自然科學(xué)版，2010，31(1)：31-34. Liu Zhan, Liu Maocheng, Wei Wei, et al. A Gravity Anomaly Separation Method Based on Cellular Neural Network[J]. Journal of China University of Petroleum: Edition of Natural Science, 2010, 31(1): 31-34.

[9] Kennddy J, Eberhart R. Particle Swarm Optimization[C]// Proceedings of IEEE International Conference on Neutral Networks. Perth: IEEE, 1995:1942-1948.

[10] Eberhart R C, Shi Y. Particle Swarm Optimization: Developments, Applicationgs and Resources[C]// Pro Congress on Evolutonary Computation 2001. Piscataway: IEEE, 2001: 81-86.

[11] 盧珊萍,于盛林.基于粒子群算法的細(xì)胞神經(jīng)網(wǎng)絡(luò)模板參數(shù)設(shè)計(jì)[J].計(jì)算機(jī)技術(shù)與發(fā)展,2009,19(4):83-86. Lu Shanping, Yu Shenglin. A Template Design Method for Cellular Neural Network Based Particle Swarm Optimizer Algorithm[J]. Computer Technology and Development, 2009,19(4):83-86.

[12] Karakuzu C. Parameter Tuning of Fuzzy Sliding Mode Controller Using Particle Swarm Optimization[J]. International Journal of Innovative Computing Information and Control, 2010, 6(10): 4755-4770.

[13] 馬國(guó)慶,黃大年,杜曉娟,等.Hartley變換在位場(chǎng)(重、磁)異常導(dǎo)數(shù)計(jì)算中的應(yīng)用[J].吉林大學(xué)學(xué)報(bào):地球科學(xué)版,2014,44(1):328-335. Ma Guoqing, Huang Danian, Du Xiaojuan, et al. Hartley Transform in the Application of the Derivatives Date[J]. Journal of Jilin University: Earth Science Edition, 2014, 44(1):328-335.

[14] 管志寧，侯俊勝.重磁異常反演的擬BP神經(jīng)網(wǎng)路方法及其應(yīng)用[J].地球物理學(xué)報(bào)，1998，41(2)：242-251. Guan Zhining, Hou Junsheng. Inversion of Gravity and Magnetic Anomalies Using Pseudo-BP Neural Network Method and Its Application[J]. Chinese Journal of Geophysics, 1998，41(2)：242-251.

[15] 曾琴琴,王永華,吳文賢.二維磁異常的粒子群快速成像方法及其應(yīng)用[J].吉林大學(xué)學(xué)報(bào):地球科學(xué)版,2013,43(2):616-622. Zeng Qinqin, Wang Yonghua, Wu Wenxian. Fast Imaging of 2D Magnetic Anomaly by Particle Swarm Optimization and Its Application[J]. Journal of Jilin University: Earth Science Edition, 2013, 43(2): 616-622.

[16] Chua L, Thiran P. An Analytic Method for Designing Simple Cellular Neural Networks[J]. IEEE Trans on Circuits and Systems, 1991, 38(11):1332-1341.

[17] 范正國(guó),劉前坤,黃旭釗,等.自適應(yīng)重磁對(duì)應(yīng)分析[J].吉林大學(xué)學(xué)報(bào):地球科學(xué)版,2012,42(6):1904-1919. Fan Zhengguo, Liu Qiankun, Huang Xuzhao, et al. Adaptive Correlative Analysis Between Aeromagnetic and Gravity Anomaly[J]. Journal of Jilin University: Earth Science Edition, 2012, 42(6): 1904-1919.

[18] Zarandy A, Stoeffets A K, Roska T,et al. Implementation of Binary and Gray-Scale Mathematical Morphology on the CNN Universal Machine[J]. IEEE Trans on CAS, 1998,45(2): 163-168.

[19] 侯遵澤，楊文采.中國(guó)重力異常的小波變換與多尺度分析[J].地球物理學(xué)報(bào)，1997，40(1):85-95. Hou Zunze, Yang Wencai. Wavelet Transform and Multi-Scale Analysis on Gravity Anomalies of China[J]. Chinese Journal of Sinica, 1997,40(1):85-95.

[20] 楊文采，施志群，侯遵澤，等.離散小波變換與重力異常多重分解[J].地球物理學(xué)報(bào), 2001，44(4)：534-541. Yang Wencai, Shi Zhiqun, Hou Zunze, et al. Discrete Wavelet Transform for Multiple Decomposition of Gravity Anomalies[J]. Chinese Journal of Geophysics, 2001，44(4)：534-541.

[21] Hinze W J, Aiken C V,Brozena J M,et al. New Standard for Reducing Gravity Data: The North American Gravity Database[J]. Geophysics, 2005,70(4): 346-365.

[22] Wu Jianshe, Jiao Licheng. Synchronization in Dynamic Networks with Nonsymmetrical Time-Delay Coupling Based on Linear Feedback Controllers[J]. Physica, 2008, 387(3): 2111-2119.

[23] Zhong W S, Stefanovski J D,Dimirovski G M,et al. Decentralized Control and Synchronization of Time-Varying Complex Dynamical Network[J]. Physica, 2009, 142(7): 151-167.

[24] Krink T, Vesterstrom J S, Riget J. Particle Swarm Optimization with Spatial Particle Extension[C]//Proceeding of International Conference on Evolutionary Computation.[S.l.].IEEE, 2002: 1474-1497.

Gravity Anomaly Extraction for Hydrocarbon Based on Particle Swarm Optimization and Cellular Neural Network

Du Runlin1, 2, Liu Zhan1

1.SchoolofGeosciences,ChinaUniversityofPetroleum(EastChina),Qingdao266555，Shandong,China2.TheKeyLaboratoryofMarineHydrocarbonResourcesandEnvironmentGeology,QingdaoInstituteofMarineGeology,Qingdao266071,Shandong，China

Cellular neural network is an effective method for extracting useful anomalies with a strong capability of horizontal resolution. Its key point is to optimize the template when using cellular neural network to extract hydrocarbon gravity anomaly information. We presented an improved algorithm based on the characteristics of the cellular neural network and the particle swarm optimization, joined constraint condition in template, and introduced the shrinkage factors in particle swarm optimization individual characteristics through the analysis of the dynamic performance of cellular neural network. Through the model testing and the actual data processing, the algorithm proved to be more stable to optimize the template for extracting reservoir gravity anomaly information.

cellular neural network; particle swarm optimization; gravity anomaly；shrinkage factors;petroleum gas

10.13278/j.cnki.jjuese.201503301.

2014-08-15

國(guó)家“863”計(jì)劃項(xiàng)目(2012AA09A20107)；國(guó)土資源部地質(zhì)調(diào)查項(xiàng)目(GZH200900501)；國(guó)家自然科學(xué)基金青年基金項(xiàng)目(41306062)

杜潤(rùn)林(1987--)，男，博士研究生，主要從事地球物理資料反演研究，E-mail:durunlin123@163.com

劉展(1957--)，男，教授，博士生導(dǎo)師，主要從事區(qū)域地球物理資料處理與解釋研究，E-mail:liuzhan5791@sina.com。

10.13278/j.cnki.jjuese.201503301

P631.1

A

杜潤(rùn)林，劉展. 基于粒子群優(yōu)化的細(xì)胞神經(jīng)網(wǎng)絡(luò)油氣重力異常信息提取.吉林大學(xué)學(xué)報(bào):地球科學(xué)版,2015,45(3):926-933.

Du Runlin, Liu Zhan. Gravity Anomaly Extraction for Hydrocarbon Based on Particle Swarm Optimization and Cellular Neural Networks.Journal of Jilin University:Earth Science Edition,2015,45(3):926-933.doi:10.13278/j.cnki.jjuese.201503301.