基于條件對抗增強的Transformer 煤礦微震定位方法?

丁琳琳 胡永亮 李昱達 王凱璐 王慧穎

（1.遼寧大學(xué)信息學(xué)院 沈陽 110036）（2.國網(wǎng)遼寧省電力有限公司信息通信分公司 沈陽 110065）

1 引言

近年來，由于煤礦生產(chǎn)采掘深度的逐漸增大和煤礦企業(yè)開采掘進速度的逐漸加快，頻繁發(fā)生一系列因采礦引發(fā)的煤礦動力災(zāi)害，比如煤礦沖擊地壓災(zāi)害和瓦斯突出災(zāi)害等［1］。為了減少事故災(zāi)害造成的人員傷亡，及時組織救援，同時監(jiān)測礦區(qū)的微震能量變化，目前大多數(shù)礦區(qū)均采用微震監(jiān)測系統(tǒng)作為煤礦動力災(zāi)害的主要監(jiān)測手段，如波蘭的SOS系統(tǒng)、南非的ISSI系統(tǒng)和加拿大的ESG系統(tǒng)以及國內(nèi)的部分系統(tǒng)等［2］。

微震監(jiān)測系統(tǒng)通過組網(wǎng)臺站接收的震源時空數(shù)據(jù)來反演震源位置，反演出微震事件發(fā)生的位置并對定位結(jié)果做出適當?shù)脑u價。但是定位結(jié)果受到監(jiān)測臺組網(wǎng)布置、定位算法選取、到時拾取、可用定位震相、地殼速度模型等諸多方面的影響［3］，導(dǎo)致震源定位結(jié)果與實際震源位置相差較大，影響了在礦山的實際應(yīng)用效果。為了避免上述諸多因素帶來的影響，文獻［4］介紹了一種將全卷積神經(jīng)網(wǎng)絡(luò)與三維高斯分布結(jié)合的方式，將定位問題轉(zhuǎn)化為分類問題求解震源坐標；文獻［5］介紹了一種將卷積神經(jīng)網(wǎng)絡(luò)CNN 與Transformer 模型結(jié)合并引入了臺站位置信息的方式，進一步提高了震源定位精準度，保證了定位精度優(yōu)于深度學(xué)習(xí)基線。

針對上述問題，本文提出一種基于條件對抗增強的Transformer 煤礦微震定位方法CGAN-Transformer，該方法首先通過條件生成對抗網(wǎng)絡(luò)CGAN［6］對現(xiàn)有的微震信號（波形數(shù)據(jù)）進行學(xué)習(xí)和模仿，利用其生成器生成與指定標簽的微震數(shù)據(jù)具有相同數(shù)據(jù)分布的偽微震數(shù)據(jù)；其次通過Transformer編碼器層將微震信號轉(zhuǎn)換成特征向量后再通過以α的倍率放大首個P 波臺站在微震信號數(shù)據(jù)中所占比重的方式抵消了不同地質(zhì)條件對微震波傳播的消極影響后，同時也利用其注意力機制進一步學(xué)習(xí)微震波形數(shù)據(jù)深層次特征和復(fù)雜的站間依賴關(guān)系；最后通過混合密度輸出層特征向量與震源位置之間的映射依賴關(guān)系，并利用高斯分布參數(shù)計算最優(yōu)的震源位置。在本文的最后，通過在兩個波形數(shù)據(jù)集上的對比實驗驗證了模型的實際效果。

2 相關(guān)工作

為了提高定位精度，有許多研究人員在此震源定位領(lǐng)域上進行了大量的研究和優(yōu)化，大致可分為三大類，即傳統(tǒng)定位算法、機器學(xué)習(xí)定位算法、深度學(xué)習(xí)定位算法。

在傳統(tǒng)定位算法中，文獻［7］介紹了一種求解到時函數(shù)的偏導(dǎo)數(shù)與逆矩陣，經(jīng)過迭代公式求解出到時誤差的極小值，即震源參數(shù)；文獻［8］介紹了一種避免了求解偏導(dǎo)數(shù)或者逆矩陣，通過選取最小到時臺站空間位置即可；文獻［9］介紹了一種搜索迭代的方式尋找符合微震事件條件震源點。以上傳統(tǒng)算法對算法的初始值的依賴性較強，如果采用靜態(tài)波速模型，在迭代求解的過程中會造成誤差累加效應(yīng)，使算法陷入局部最優(yōu)。文獻［10］雖然在一定程度上減小了定位誤差，但是影響微震事件定位精度的相關(guān)參量并不是相互獨立的，所以在求解過程中隱形地增加了困難。

在機器學(xué)習(xí)定位算法中，文獻［11］介紹了一種利用高密度臺陣的和粒子群算法分析了定位誤差與傳感器數(shù)量、位置之間的關(guān)系；文獻［12］介紹了一種利用方位角優(yōu)化的方法結(jié)合模擬退火算法與粒子群算法提出了改進全局搜索算法；文獻［13］介紹了一種差分進化遺傳算法進行尋優(yōu)的定位方法，這些方法提高定位精度的同時實現(xiàn)了從傳統(tǒng)的“數(shù)學(xué)模型”到“仿生模型”的演變。

在深度學(xué)習(xí)定位算法中，文獻［14］介紹了一種利用卷積神經(jīng)網(wǎng)絡(luò)（CNN）的方法從單臺站三分量波形數(shù)據(jù)中檢測地震事件；文獻［15］介紹了一種考慮到單臺站波形數(shù)據(jù)不足以代表整個地震事件，使用深度卷積對多臺站波形數(shù)據(jù)進行特征提取。

然而，上述的定位模型需要一定規(guī)模的標簽事件，才能在訓(xùn)練階段充分學(xué)習(xí)數(shù)據(jù)特征，并對無標簽的事件進行定位。由于數(shù)據(jù)的限制，如訓(xùn)練數(shù)據(jù)中高質(zhì)量事件的數(shù)量不足、標簽缺失、規(guī)模較小，最終導(dǎo)致模型定位精度不高、魯棒性差等問題。文獻［16］介紹了一種利用CGAN 生成數(shù)據(jù)樣本的能力，解決了在故障診斷領(lǐng)域由于數(shù)據(jù)量小導(dǎo)致的數(shù)據(jù)不平衡問題；文獻［17］介紹了一種利用遷移學(xué)習(xí)將大樣本源域地震數(shù)據(jù)的特征遷移到小樣本源域地震數(shù)據(jù)的模型中，從而改善了數(shù)據(jù)規(guī)模不足、特征不足等問題。

3 CGAN-Transformer模型結(jié)構(gòu)

模型的整體結(jié)構(gòu)可以分為三層，條件對抗增強層、Transformer編碼器層以及混合密度輸出層。條件對抗增強層：先利用判別器D來學(xué)習(xí)真實微震信號所具有的數(shù)據(jù)特征分布后，再利用生成器G生成符合真實微震分布的偽微震數(shù)據(jù)。Transformer 編碼器層：先將微震信號提取成特征向量后，再利用其注意力機制對數(shù)據(jù)特征進一步混合提取，尋找深層次微震特征和復(fù)雜的站間依賴關(guān)系，同時引入了監(jiān)測臺站的地理位置編碼并通過放大首個P波臺站的數(shù)據(jù)比重的方式，抵消由地質(zhì)結(jié)構(gòu)不確定引起的定位結(jié)果遠離煤礦區(qū)域的問題。最后通過混合密度輸出層獲取震源的高斯分布參數(shù)，根據(jù)分布參數(shù)計算最優(yōu)的震源位置。完整的模型結(jié)構(gòu)如圖1所示。

圖1 CGAN-Transformer模型結(jié)構(gòu)圖

3.1 條件對抗增強層

雖然煤礦微震監(jiān)測系統(tǒng)的監(jiān)測臺站每時每刻都在記錄著傳感器所受到的原始應(yīng)力數(shù)據(jù)，但對于本項工作來說，并不是所有的應(yīng)力數(shù)據(jù)都是有用的。只有在發(fā)生煤礦微震事件時，處于[t0,t1]區(qū)間的波形數(shù)據(jù)才是有效數(shù)據(jù)。其中，t0=tp-5s；t1=tp+25s。tp為P波到時。從微震監(jiān)測數(shù)據(jù)中發(fā)現(xiàn)的數(shù)據(jù)不平衡問題，是影響震源定位精度不高的主要原因；再加上微震事件的震源位置與微震信號是一對多的關(guān)系（即不同的微震信號可引起同震源位置的微震事件）降低了模型的震源定位精度。因此本文采用一種基于CGAN 架構(gòu)的數(shù)據(jù)增強方法對現(xiàn)有的微震數(shù)據(jù)集進行擴充增強。條件對抗增強層主要由判別器D、生成器G和條件Cond三部分組成。具體結(jié)構(gòu)如圖1（a）所示。生成器G 將輸入的隨機噪聲z 映射到輸出信號y 中；判別器D 將其輸入信號y映射成概率值；將條件Cond分別輸入到生成器G和判別器D中用控制生成器G的輸出。

首先，對判別器D 進行訓(xùn)練。該判別器由7 層卷積層和一個多層感知機（Multilayer Perceptron，MLP）組成，本文將真實微震事件的震源坐標作為條件Cond 添加到判別器D 的輸入端，判別器D 既要判斷輸入數(shù)據(jù)的真?zhèn)?，也要判斷輸入?shù)據(jù)與作為條件信息的震源坐標是否匹配。設(shè)輸入判別器D的數(shù)據(jù)X是一個三維的張量Z0n,t,c，代表一個臺站截取的有效波形數(shù)據(jù)，其中n?{1,2,3}代表微震波三分量，c=1 代表通道數(shù)。則每層卷積的計算過程如下公式所示：

其中，σ(x)代表的是LeakyReLU 函數(shù)；c代表的是第i層的第c通道；t代表的是第i層的第c通道的第t元素；同理可以得，c'和t'代表的是第i-1 層的通道和數(shù)據(jù)元素；Ci-1代表第i-1 層的通道數(shù)；s代表卷積核的步幅；Ti-1代表的是第i-1 層的第Ci-1通道的數(shù)據(jù)元素總長度。在經(jīng)過7層卷積層的內(nèi)容學(xué)習(xí)和特征提取后，通過一層MLP 將特征數(shù)據(jù)進行展平并映射到一維的評價標量上，用于指導(dǎo)生成器G的訓(xùn)練。

其次，對生成器G 進行訓(xùn)練。該生成器由7 層反卷積層和一層MLP 組成，即為判別器D 的逆過程。不同于判別器D 的是除最后一層σ(x)為Tanh函數(shù)外，其余層均為ReLU函數(shù)。

其中，Pdata(x)為真實微震事件特征分布；Pg(z)為噪聲分布；D(x|cond)表示判別器D 判斷為真實微震事件為真的概率；D(G(z|cond)|cond)表示隨機噪聲生成的偽微震事件特征分布；D（G（z））為判別器D 判斷生成器G 輸出為真的概率；c表示該真實微震事件對應(yīng)的震源標簽。為了避免模式坍塌，保證生成數(shù)據(jù)的多樣性，本文采用Wasserstein距離［19］作為Pdata(x)分布與Pg(z)分布之間距離的衡量標準。

3.2 帶有位置嵌入的Transformer編碼器層

Transformer 模型是一種基于注意力機制的序列到序列的經(jīng)典模型，GPT 以及Bert等大型的語言模型都是基于Transformer 模型衍變而來的。該模型的主要貢獻有以下兩點：第一點則是通過引入注意力機制和位置編碼的方式，解決了傳統(tǒng)網(wǎng)絡(luò)記憶長度過短，提取出的信息無法保留較多上下文信息，以至于難以處理長序列數(shù)據(jù)等的問題；第二點則是通過引入多頭注意力機制，解決傳統(tǒng)NLP模型中并行計算能力不足的問題。

本文通過注意力機制來量化在一次微震事件中每個激發(fā)臺站（波形特征和位置）和其他臺站之間的相關(guān)度或依賴程度，從而在確定震源位置時知道更偏向哪幾個臺站。Transformer 編碼器層的具體結(jié)構(gòu)如圖1（b）所示。

首先，由于條件對抗增強層是直接對微震波形數(shù)據(jù)進行增強，屬于線下增強。因此，在進行自注意力運算之前需要重新對波形數(shù)據(jù)進行特征提取。該卷積提取與式（1）的不同之處在于：σ(x)代表的是ReLU 函數(shù)；Ci-1={1,8,32,64,128,64,32,8}。在經(jīng)過卷積提取之后，先將特征向量展平，再通過多層感知機MLP 將特征向量微震映射到dmodel維。

其次，目前將臺站的位置信息嵌入到特征向量中有二種方式：第一種方式為ConCatenate，即直接在特征向量的尾部串聯(lián)上對應(yīng)臺站的經(jīng)緯度和深度信息，則嵌入之后的特征向量維度為；第二種方式為PositionEmbedding，即通過三角函數(shù)將臺站位置信息映射到特征向量中，特征向量維度保持不變。本文采用基于不同波長的正弦函數(shù)對臺站位置信息進行編碼，具體編碼計算方式如下公式所示：

其中，PE 為二維矩陣[numstation,dmodel]，矩陣的列數(shù)和上游提取出來的波形特征向量列數(shù)保持一致；矩陣PE 的每一行代表微震事件觸發(fā)臺站的位置向量，列表示位置向量的維度。pos1,2,3代表的是臺站的經(jīng)度、緯度和深度；經(jīng)度即該位置向量的第一個維度。經(jīng)度是地理位置的一個重要指標，用來表示一個地點在東西方向上的位置。緯度即該位置向量的第二個維度。緯度是地理位置的另一個重要指標，用來表示一個地點在南北方向上的位置。深度即該位置向量的第三個維度。深度是用來表示一個地點在垂直方向上的位置，通常用距離地表的距離來表示。dmodel表示位置向量維度，i表示位置向量具體的維度，i的取值范圍是[0,dmodel/2-1]；PE(pos1,2,3,2i)代表在臺站位置向量中的偶數(shù)位置的值，PE(pos1,2,3,2i+1)表示奇數(shù)位置的值。對經(jīng)度、緯度、深度全部按照該編碼函數(shù)映射成dmodel維的向量后再逐元素相加，組合成一個位置向量。

理論上來說，當一個微震事件發(fā)生時距離震源最近的臺站是第一個接收到微震P 波的；但由于震源到臺站之間的地質(zhì)結(jié)構(gòu)復(fù)雜多變，吸收波、反射波等情況的存在，所以實際上不一定是離震源最近的臺站最先接收到P 波?？紤]到上述情況的存在，本文以概率α?[0,1]來模擬這種情況的發(fā)生：當α=0 時，則認為存在上述情況，不做任何處理；當α=1 時，則認為不存在上述情況，通過放大特征向量的方式表達該特征向量與震源的關(guān)聯(lián)度最大；當α=(0,1)時，則認為上述情況不同程度的存在；隨機變量α的取值符合高斯分布。對于微震事件激發(fā)的任意臺站，其位置向量展開如下公式所示：

3.3 混合密度輸出層

本文提出的煤礦微震震源定位模型中，Transformer 編碼器層的輸出只能獲取包含上下文及臺站位置信息的特征向量，還需要考慮到輸入的微震波形數(shù)據(jù)與最終震源位置標簽之間的依賴關(guān)系。因此，模型中采用混合密度輸出層建立特征向量到震源位置的依賴關(guān)系來獲取全局最優(yōu)的震源位置?；旌厦芏染W(wǎng)絡(luò)是一種特殊的用于密度估計的高斯混合模型，該模型常常被應(yīng)用于多元預(yù)測任務(wù)上［20］?；旌厦芏容敵鰧拥木唧w結(jié)構(gòu)如圖1（c）所示。

混合密度網(wǎng)絡(luò)（Mixture Density Network，MDN）是一種神經(jīng)網(wǎng)絡(luò)模型，也是一種生成模型，用于對多峰分布的條件概率密度函數(shù)進行建模。MDN 的思想是通過利用混合高斯分布來建立對數(shù)據(jù)的概率密度函數(shù)。即：每個高斯分布函數(shù)都由一個均值、方差和加權(quán)系數(shù)組成，其中加權(quán)系數(shù)代表的是每個高斯分布函數(shù)所占的比重。

混合密度網(wǎng)絡(luò)先假設(shè)某一區(qū)域所有微震事件為總體樣本（包括未發(fā)生的事件）X=x1,x2,…,xn，并且該樣本集的總體分布是由K 個高斯分布加權(quán)混合而成，則對于微震事件來說屬于該樣本集的概率就是分別屬于k個高斯分布概率的加權(quán)和。即

4 實驗結(jié)果與分析

為了驗證本文提出的基于Transformer 條件對抗增強的煤礦震源定位方法的有效性，同時與其他基于深度學(xué)習(xí)的震源定位模型進行對比，本文采用遼寧某礦真實數(shù)據(jù)集和智利地震數(shù)據(jù)集（M<1.5）進行實驗驗證。

4.1 實驗數(shù)據(jù)

實驗數(shù)據(jù)來源于遼寧某礦的真實微震數(shù)據(jù)和智利公開的地震數(shù)據(jù)集。其中遼寧某礦的所有微震數(shù)據(jù)均為微震監(jiān)測系統(tǒng)在非人為干預(yù)情況下的自動拾取，計算過程中涉及的運算參數(shù)與微震信號所在礦區(qū)實際匹配。具體包括了從2019年到2021年礦區(qū)所發(fā)生的全部微震事件，共計240 個微震事件；智利數(shù)據(jù)集選取震級小于1.5，觸發(fā)臺站個數(shù)等于20的微地震事件，共計2665個微地震事件；兩個數(shù)據(jù)集的詳細統(tǒng)計如表1所示。

表1 數(shù)據(jù)集信息統(tǒng)計表

4.2 實驗環(huán)境

實驗?zāi)Ｐ筒捎糜蒄acebook 人工智能研究院（FAIR）搭建的Pytorch 框架。其中，條件對抗增強模型訓(xùn)練了大致3.5 個小時，Transformer 震源定位模型訓(xùn)練了大致0.8 個小時，實驗中使用了Dropout層防止模型過擬合。具體的實驗環(huán)境如表2 所示，實驗中采用的超參數(shù)如表3所示。

表2 實驗環(huán)境統(tǒng)計表

表3 實驗超參數(shù)設(shè)置表

其中，在條件對抗增強模型中1e-4為生成器G的學(xué)習(xí)率，3e-4為判別器D的學(xué)習(xí)率。并且該學(xué)習(xí)率隨著迭代次數(shù)的增加而線性衰減，幫助模型在后期快速收斂。在Batch大小的選擇上通過采用觸發(fā)臺站個數(shù)的倍數(shù)，進而控制輸入模型批量訓(xùn)練個數(shù)為微震事件的整數(shù)倍。

4.3 實驗結(jié)果

為了驗證本文提出的基于Transformer 條件對抗增強的煤礦微震定位方法的性能，同時比較其他的微震震源定位方法，本節(jié)引入基于卷積神經(jīng)網(wǎng)絡(luò)的CNN-EL 模型［15］和基于圖卷積神經(jīng)網(wǎng)絡(luò)的DGNN 模型［21］作為對比模型。在本文提出的基于Transformer 條件對抗增強的煤礦微震定位模型中不僅增加了條件對抗增強層，還對Transformer編碼器層的位置向量嵌入進行了改進，其中條件對抗增強層可以將隨機噪聲轉(zhuǎn)換為指定震源對應(yīng)的波形數(shù)據(jù)，Transformer編碼器層引入了隨機變量。為了驗證這些改動的有效性，本文引入原生的Transformer 震源定位模型在智利數(shù)據(jù)集和遼寧某礦數(shù)據(jù)集上作進一步對比實驗，其中Transformer為原生模型，Concatenate 與PositionEmbedding 均為本文提出的基于條件對抗增強的Transformer 微震定位模型，區(qū)別在于監(jiān)測臺站位置信息的嵌入方式不同。

五個模型分別在智利數(shù)據(jù)集和遼寧某礦數(shù)據(jù)集上進行對比實驗，采用震中誤差（Epicenter Error）和震源誤差（Hypocentral Error）作為評價指標，通過最終的實驗結(jié)果驗證該微震震源定位方法的有效性。其中震中誤差EE 和震源誤差HE 分別表示為

其中，(x,y,z)為震源的預(yù)測坐標；(x',y',z')為震源的真實坐標，五個模型在智利數(shù)據(jù)集和遼寧某礦數(shù)據(jù)集上的震中誤差變化曲線如圖2 所示、震源誤差變化曲線如圖3所示。

圖2 五種模型在兩個數(shù)據(jù)集上的震中誤差變化曲線

圖3 五種模型在兩個數(shù)據(jù)集上的震源誤差變化曲線

從圖2 中可以看出，本文提出的基于Transformer 條件對抗增強的微震定位方法在智利數(shù)據(jù)集和遼寧某礦數(shù)據(jù)集上均取得了最優(yōu)成績，在智利數(shù)據(jù)集上震中誤差降低到6.29km，比CNN-EL模型和DGNN 模型分別降低72%和54%；改進后的CGAN-Transformer 模 型 比 原 生Transformer 模 型 震中誤差降低34%，PositionEmbedding方式比Concatenate 方法震中誤差降低了8%。在遼寧某礦數(shù)據(jù)集上，CGAN-Transformer 模型震中誤差比CNN-EL模型和DGNN 模型分別降低46%和34%；改進后的CGAN-Transformer 模 型 比 原 生Transformer 模 型 震中誤差降低11%，PositionEmbedding方式比Concatenate 方法震中誤差降低了2%。這說明不管是在智利數(shù)據(jù)集上還是遼寧某礦數(shù)據(jù)集上，本文提出的基于條件對抗增強的Transformer 微震定位模型的精準度是優(yōu)于CNN-EL 模型和DGNN 模型以及原生Transformer模型的，而監(jiān)測臺站位置的嵌入方式在智利數(shù)據(jù)集上對定位精準度的影響較大，在遼寧某礦數(shù)據(jù)集上的對定位精度的影響較小。

從圖3 中可以看出，本文提出的基于Transformer 條件對抗增強的微震定位方法在智利數(shù)據(jù)集和遼寧某礦數(shù)據(jù)集上均取得了最優(yōu)成績，在智利數(shù)據(jù)集上震中誤差降低到6.52km，比CNN-EL模型和DGNN 模型分別降低75%和57%；改進后的CGAN-Transformer 模 型 比 原 生Transformer 模 型 震中誤差降低42%，PositionEmbedding方式比Concatenate 方法震中誤差降低了21%。在遼寧某礦數(shù)據(jù)集上，CGAN-Transformer 模型震中誤差比CNN-EL模型和DGNN 模型分別降低46%和36%；改進后的CGAN-Transformer 模 型 比 原 生Transformer 模 型 震中誤差降低13%，PositionEmbedding方式比Concatenate方法震中誤差降低了5%。

從圖4 中可以看出，五個模型在兩個數(shù)據(jù)集上的誤差大致呈右偏分布，本文提出的基于條件對抗增強的Transformer 微震定位模型在智利數(shù)據(jù)集上表現(xiàn)較好，誤差分布較為集中，其震中誤差峰值和震源誤差峰值都低于10km，其中位誤差也低于其他四個模型，定位效果較好。但是本文提出的基于條件對抗增強的Transformer 微震定位模型在遼寧某礦數(shù)據(jù)集上的誤差分布較為離散、數(shù)據(jù)波動較大、定位效果不穩(wěn)定，這可能是因為遼寧某礦數(shù)據(jù)集中可能存在一些數(shù)據(jù)質(zhì)量問題，例如數(shù)據(jù)噪聲、不準確的測量值等，并且數(shù)據(jù)集的數(shù)據(jù)規(guī)?？赡茌^小，相比智利數(shù)據(jù)集而言可能數(shù)據(jù)量有限。這些問題會導(dǎo)致模型在處理數(shù)據(jù)時產(chǎn)生不穩(wěn)定性，從而導(dǎo)致誤差分布較為離散。但其中位誤差依舊低于其他四個模型。因此，可以得出條件對抗增強層的引入以及首個P 波臺站隨機增強的引入提升了CGAN-Transformer模型的最終定位效果。

圖4 誤差分布式小提琴圖

5 結(jié)語

本文針對煤礦微震定位領(lǐng)域中出現(xiàn)的微震事件單一導(dǎo)致的定位模型泛化弱、定位精度不高，提出了一種基于條件對抗增強的Transformer 煤礦微震定位方法，該方法通過對抗增強的方式自主學(xué)習(xí)已有微震事件的特征并由生成器G 生成符合微震特征的偽微震事件對現(xiàn)有微震事件進行多樣變種，在Transformer 編碼器層通過高斯隨機變量模擬在不同地質(zhì)條件下首個P 波臺站情況。實驗表明本文提出的CGAN-Transformer 模型在智利數(shù)據(jù)集上的震中誤差和震源誤差分別降低到6.29km 和6.52km，比改進前分別降低了34%和42%（即誤差降低了38%）；其在遼寧某礦數(shù)據(jù)集上的震中誤差和震源誤差比改進前分別降低了11%和13%（即誤差降低了12%），說明了改進的有效性。在后續(xù)工作中，將在更多的數(shù)據(jù)集上與更多其他的基于深度學(xué)習(xí)的模型進行對比，同時進一步優(yōu)化模型。