李瑞友，張淮清，吳昭

(重慶大學(xué) 輸配電裝備及系統(tǒng)安全與新技術(shù)國家重點實驗室,重慶 400044)

0 引言

瞬變電磁法(transient electromagnetic method, TEM)是地球物理勘探中最常用的方法之一，廣泛應(yīng)用于礦產(chǎn)、地下水勘探及淺層地質(zhì)調(diào)查等領(lǐng)域[1-3]。瞬變電磁法是一種利用快速關(guān)斷脈沖電流產(chǎn)生二次接收電壓的電磁探測方法，而TEM反演是根據(jù)其得到的探測數(shù)據(jù)推導(dǎo)出地電結(jié)構(gòu)參數(shù)，后者是一個典型的具有非線性特征和不適定參數(shù)估計的反演問題，它一直是地球物理研究中的重要問題。在過去的幾十年里，各種解釋地電阻率數(shù)據(jù)的方法相繼發(fā)表[4-6]。其中，人工神經(jīng)網(wǎng)絡(luò)在TEM數(shù)據(jù)解釋中展現(xiàn)了最為廣闊的潛力。

人工神經(jīng)網(wǎng)絡(luò)反演算法能夠模擬人腦神經(jīng)元信息學(xué)習(xí)和處理，直接從樣本數(shù)據(jù)中構(gòu)造數(shù)學(xué)模型，以較強的學(xué)習(xí)能力和泛化能力在地球物理反演解釋中得到廣泛應(yīng)用[7-10]。其中，Srinivas等[11]比較了一維垂直電測深數(shù)據(jù)中BPNN (back propagation neural network)、RBFNN(radial basis function neural network)和GRNN(general regression neural network)的反演性能表現(xiàn)；Maiti等[12]研究了一種以一維電測深數(shù)據(jù)為基礎(chǔ)的貝葉斯神經(jīng)網(wǎng)絡(luò)訓(xùn)練方法；Johnson等[13]利用人工神經(jīng)網(wǎng)絡(luò)方法成功地估算了層狀大地模型電阻率參數(shù)。綜上所述，在地球物理反演解釋中神經(jīng)網(wǎng)絡(luò)得到了廣泛的應(yīng)用(具有良好的非線性擬合能力)，但在瞬變電磁反演應(yīng)用中仍然存在收斂速度慢、精度低以及訓(xùn)練中過擬合等局限性[14]。因此，復(fù)雜的電磁響應(yīng)數(shù)據(jù)和模型參數(shù)(電阻率和層厚)反演需要更高效的神經(jīng)網(wǎng)絡(luò)模型和更完善的優(yōu)化程序來獲得更好的結(jié)果和解釋。

極限學(xué)習(xí)機是由Huang等[15]提出的一種新的單隱含層前饋神經(jīng)網(wǎng)絡(luò)(single hidden layer feedforward neural network, SLFN)模型，其輸入層權(quán)值和隱含層偏差均是隨機產(chǎn)生，輸出層權(quán)值通過解析計算可獲得唯一的全局最優(yōu)解，很好地克服了傳統(tǒng)神經(jīng)網(wǎng)絡(luò)(如BPNN、RBFNN等)存在的訓(xùn)練時間長、易陷入局部極值及過擬合等問題。該算法具有結(jié)構(gòu)簡單、學(xué)習(xí)速度快和泛化性能強等特點，已廣泛應(yīng)用于回歸與分類[16-17]、圖像識別[18-19]和決策支持[20]等領(lǐng)域。但是，傳統(tǒng)的ELM基本上是批量學(xué)習(xí)算法，即一次將樣本輸入訓(xùn)練，得到輸出，然而在大多數(shù)情況下，數(shù)據(jù)樣本都不能一次性得到，而是不斷地有新的樣本加入，因此限制了ELM在實時或非平穩(wěn)情況下的適用性。最近，Liang等[21]提出了一種在線慣序極限學(xué)習(xí)機(OSELM)，它可以通過固定或不同長度的逐個或逐塊學(xué)習(xí)訓(xùn)練數(shù)據(jù)。與傳統(tǒng)極限學(xué)習(xí)機算法相比，OSELM可以提供更好的泛化性能，以更快的速度學(xué)習(xí)?；谶@些優(yōu)點，OSELM在預(yù)測問題領(lǐng)域得到了廣泛的應(yīng)用[22-24]。

本文針對TEM反演的高維、非凸及非線性特征，提出了一種在線慣序極限學(xué)習(xí)機反演算法，該方法通過隨機設(shè)定隱層參數(shù)來簡化反演模型的學(xué)習(xí)過程，提高了反演算法的計算效率；同時，OSELM的Moore-Penrose廣義逆計算過程和逐塊式學(xué)習(xí)訓(xùn)練數(shù)據(jù)保證了反演算法的全局最優(yōu)和泛化能力，較好地解決了TEM的局部收斂問題。文中還給出了TEM正演方法、樣本采集方法及非線性反演模型的流程，并通過兩個經(jīng)典的層狀地電模型和一個異常體模型反演驗證本文算法的有效性。

1 層狀地電模型

1.1 TEM正演

假設(shè)n層水平大地各層介質(zhì)的電阻率和層厚表示為ρ1,h1;ρ2,h2;…;ρn,hn;hn→∞。在地表(h=0)的中心回線發(fā)射線圈Tx向地下發(fā)射斜階躍波形式的激勵信號，同時利用在地表的接收線圈Rx采集具有地下信息的電磁響應(yīng)數(shù)據(jù)Hz，如圖1所示。

圖1 層狀地電模型和TEM方法示意Fig.1 Schematic diagram of the layered geoelectric model and TEM method

在準靜態(tài)(忽略位移電流)下柱坐標系中，層狀大地中心回線源地表(z=0)處的垂直磁場頻域響應(yīng)表達式為[25]：

采用Hankel變換求解上式中的貝塞爾函數(shù)，再通過Gaver-Stehfest變換實現(xiàn)頻—時域的轉(zhuǎn)換，可得到垂直磁場時域響應(yīng)：

式中：Kn為系數(shù)；sn=ln2·t-1·n,n由計算機位決定，一般n=12。

1.2 TEM數(shù)據(jù)采集方法

針對層狀大地電阻率測深進行非線性反演，反演算法訓(xùn)練數(shù)據(jù)集來源于瞬變電磁正演理論計算方法。其中，中心回線瞬變電磁法的各項理論參數(shù)設(shè)置如下：發(fā)射線圈半徑r=100 m，發(fā)射電流為斜階躍響應(yīng)，關(guān)斷時間為1 μs，幅值為1 A。由于在野外數(shù)據(jù)采集中往往含有不可避免的噪聲，于是添加5%高斯噪聲到正演響應(yīng)計算的代碼中。

通常采用隨機改變電阻率和層厚度值構(gòu)建訓(xùn)練樣本的地電模型，但是，由于實際構(gòu)造不可能發(fā)生劇烈的地電參數(shù)突變[26]，所以TEM反演一般采用限定約束多層模型，可以降低模型的復(fù)雜性?；诖?，提出了一種去除冗余特征的樣本分布極限策略，每層地電模型參數(shù)構(gòu)造方式如下:

i=0,1,…,N-1。

式中：fix為取整函數(shù)，pmin、pmax分別是地電模型參數(shù)(電阻率和厚度)的最小值和最大值，N為樣本數(shù)目。為了使本文的算法能夠有效適用于常規(guī)的觀測數(shù)據(jù)，應(yīng)使模型參數(shù)取值范圍盡可能大，因此設(shè)置pmin和pmax如表1所示(以5層和9層地電模型為例)。文中算法的輸入層節(jié)點數(shù)為31個(時間采樣點Hz)，輸出層節(jié)點數(shù)為地電結(jié)構(gòu)參數(shù)數(shù)目(電阻率和層厚)。通過以上模型構(gòu)造方式，設(shè)計了2.5萬組用于訓(xùn)練和測試的樣本模型，其中5 000組作為測試，其他組作為訓(xùn)練。

表1 5層、9層地電模型各層參數(shù)的最大值與最小值

1.3 在線慣序極限學(xué)習(xí)機

ELM是由Huang等[15]提出的一種非迭代學(xué)習(xí)的神經(jīng)網(wǎng)絡(luò)算法。ELM的顯著優(yōu)點是所有隱藏節(jié)點參數(shù)(輸入權(quán)值和偏差)都是隨機生成的，并利用Moore-Penerose廣義逆計算輸出權(quán)值，克服了基于梯度下降學(xué)習(xí)理論的機器學(xué)習(xí)的固有缺陷，具有學(xué)習(xí)速度快、不易陷入局部極值、泛化能力強等優(yōu)點。其三層的網(wǎng)絡(luò)拓撲結(jié)構(gòu)如圖2所示。

圖2 ELM神經(jīng)網(wǎng)絡(luò)結(jié)構(gòu)Fig.2 ELM neural network structure

在圖2中，x為ELM的輸入值,f(x)為ELM的預(yù)測輸出值,G(w,b,x)為ELM的激活函數(shù)(本文均采用Sigmoid激活函數(shù))，β為ELM的輸出權(quán)值。隱含層用于擬合從垂直磁場到地電模型的本構(gòu)關(guān)系。

OSELM是在ELM的基礎(chǔ)上由Liang等[21]提出的一種增量學(xué)習(xí)算法，它的主要目的是解決分批次到達數(shù)據(jù)的學(xué)習(xí)問題(ELM本質(zhì)上是批處理/批量或離線學(xué)習(xí)方法)。OSELM過程可分為兩部分(圖3)：第一，在初始化階段，OSELM由給定的少量訓(xùn)練數(shù)據(jù)估計初始輸出權(quán)重矩陣β(0)和初始矩陣H0；第二，在線學(xué)習(xí)階段，每當接收到新的數(shù)據(jù)塊時，將再次運行ELM并獲得新的輸出權(quán)重。新、舊輸出權(quán)重合并，產(chǎn)生輸出權(quán)重矩陣β0以完成神經(jīng)網(wǎng)絡(luò)的更新。算法1給出了OSELM偽代碼，其中I為恒等矩陣，詳見文獻[21] 。

圖3 OSELM算法計算流程Fig.3 Block diagram of the calculation flow of the OSELM algorithm

得到初始輸出權(quán)值矩陣：

在線慣序?qū)W習(xí)階段：當接收到新的數(shù)據(jù)樣本時，假設(shè)有N1樣本進入模型，可以根據(jù)ELM的思想得到

其中，β(1)由β(0)表示，即

因此，可以得到在線學(xué)習(xí)的遞推關(guān)系

最后，設(shè)置迭代次數(shù)k=k+1，然后重復(fù)在線順序?qū)W習(xí)階段。

2 層狀地電模型反演與分析

2.1 評價指標

本文將展示OSELM算法進行層狀地電模型反演示例，算法仿真在Matlab R2016a中進行，PC采用Intel(R) Core(TM) i5-7500處理器，主頻為3.40 GHz，內(nèi)存為8.0 GB。為了客觀地評價所提方法在瞬變電磁層狀反演中的性能，采用絕對百分比誤差(absolute percentage error, APE)、相對均方根誤差(relative root mean square error, RRMSE)和相關(guān)系數(shù)(R2)進行評估，具體公式如下。

絕對百分比誤差(APE)：

相對均方根誤差(RRMSE)：

相關(guān)系數(shù)(R2)：

2.2 反演結(jié)果與比較

為了驗證OSELM反演算法的性能，進行了2個實驗,實驗1為2種學(xué)習(xí)機的反演性能比較，實驗2為綜合實例分析。

在實驗1中，對OSELM與ELM收集的TEM數(shù)據(jù)集進行反演分析。經(jīng)過反復(fù)試驗驗證，采用相同數(shù)量的隱藏節(jié)點(隱藏節(jié)點為1 000)，兩種極限學(xué)習(xí)機方法反演性能的比較結(jié)果見表2。

表2 兩種極限學(xué)習(xí)機方法反演性能的比較

從表2中不難看出，整體而言，本文提出的OSELM反演精度要高于ELM。OSELM方法有最低的APE、RRMSE和最高R2，而且所需反演時間更短,因此，該方法具有更大的優(yōu)越性，可以有效地減少非線性特性引起的反演誤差。OSELM反演質(zhì)量高的原因在于引入逐塊批量式學(xué)習(xí)策略，增強極限學(xué)習(xí)機的泛化能力，提高了反演方法的性能。在實驗2中，用一個復(fù)雜的綜合數(shù)據(jù)集來證明該方法的一致性。

2.3 層狀地電模型反演

為了驗證OSELM反演算法的性能，選取了兩種典型的地電模型(表3),合成數(shù)據(jù)集使用上述相同的正演方法和參數(shù)生成,且合成數(shù)據(jù)(垂直磁場數(shù)據(jù))受到5%高斯噪聲的干擾。圖4是OSELM與ELM方法反演一維地電模型構(gòu)造圖，圖5是一維反演模型的正演響應(yīng)數(shù)據(jù)與理論數(shù)據(jù)的擬合曲線。

表3 5層、9層地電模型反演理論值

圖4表明這兩種方法的反演結(jié)果與理論模型基本一致，但OSELM方法更精確，與理論模型的吻合程度更好。同時，這兩種方法淺層地電參數(shù)反演精度不錯，但隨著模型深度增加，地電參數(shù)略有突變，反演精度降低。從圖5也可以看出，正演響應(yīng)曲線與理論響應(yīng)曲線后期的不擬合程度更高(通常響應(yīng)后期反映著深層地電模型)，但是OSELM比ELM始終具有更低的反演偏差。因此，以上結(jié)果表明，OSELM具有最佳的泛化能力，也證明了該方法對TEM數(shù)據(jù)反演的有效性。

圖4 5層地電模型(a)與9層地電模型(b)的不同算法反演結(jié)果Fig.4 Inversion results of different algorithms for of 5-layer(a) and 9-layer(b) geoelectric models

圖5 5層地電模型(a)與9層地電模型(b)的正演響應(yīng)曲線Fig.5 Forward response curves of 5-layer(a) and 9-layer(b) geoelectric models

2.4 擬二維模型反演驗證

利用5層地電模型驗證OSELM算法對瞬變電磁法反演的可行性，理論模型及測點位置如圖6所示。圖7為OSELM和ELM反演擬剖面，從中可以看出這兩種方法反演結(jié)果與實際模型較吻合，ELM算法反演低電阻率曲線分布不夠均勻，中間模型邊界區(qū)域突變明顯,而OSELM算法反演的異常體位置稍微均勻點，形狀較清晰，反演數(shù)據(jù)更接近實際模型。顯然，與ELM相比，OSELM對層狀模型反演效果更好，更適合于地電模型反演。因此，OSELM算法在瞬變電磁反演問題上具有很高的適用性。

圖6 擬二維模型示意圖及測量位置Fig.6 Schematic diagram of quasi two dimensional model and measurement position

圖7 擬二維地電模型OSELM方法(a)和ELM方法(b)的反演結(jié)果Fig.7 Inversion results of OSELM method (a) and ELM method (b) of quasi-two-dimensional geoelectric model

3 結(jié)論

本文提出了一種在線慣序極限學(xué)習(xí)機(OSELM)的瞬變電磁層狀大地電阻率測深非線性快速反演的方法。層狀地電理論模型反演結(jié)果表明:①本文提出的去除冗余特征的樣本采集極限策略，能夠更好地反映層狀地電模型非線性反演的輸入輸出特性；選取了五層和九層地電模型的最佳OSELM網(wǎng)絡(luò)結(jié)構(gòu)，取得了較好的反演效果;②證明了OSELM算法能夠減少對初始模型的依賴，不易陷入局部極值，可以有效地應(yīng)用于層狀大地TEM數(shù)據(jù)反演中，并且同ELM算法相比，該算法反演精度更高，穩(wěn)定性更好，即泛化能力和學(xué)習(xí)能力更強，具備更好的適應(yīng)性, 能夠有效運用于一維電阻測深的神經(jīng)網(wǎng)絡(luò)反演。

總之，通過瞬變電磁層狀地電模型驗證了OSELM方法在TEM反演中的可行性和有效性。這是一種新穎的反演技術(shù)，下一步，我們將研究如何提高其實際應(yīng)用效果，并完善反演算法的性能評價指標。