非線性反演方法指通過觀測數(shù)據(jù)及一般原理或模型確定與之成非線性關(guān)系的模型參數(shù)估計(jì)值問題[1]。隨著當(dāng)前油氣勘探開發(fā)的不斷深入，非線性反演不僅包括微震數(shù)據(jù)處理、解釋，而且滲透到儲(chǔ)層預(yù)測、油藏監(jiān)測等油氣勘探開發(fā)各個(gè)領(lǐng)域，成為地球物理學(xué)家和油藏工程師共同關(guān)注的前沿課題。近年來,國內(nèi)外學(xué)者先后發(fā)表及實(shí)現(xiàn)了很多具有巨大潛力和應(yīng)用價(jià)值的非線性反演方法，如蒙特卡洛法(MCM) ,模擬退火法(SA)、人工神經(jīng)元法(ANN)及遺傳算法(GA)等[2]。非線性反演法是基于最優(yōu)化原理提出的模型,從大量已知模型正演結(jié)果中選出方差(或其他范數(shù)規(guī)則)最小的模型作為待求模型的解。

    本文針對非線性反演方法存在的極易陷入局部最優(yōu)值、分辨率低、迭代計(jì)算量大等問題，利用前沿交叉學(xué)科的最新進(jìn)展，基于BP-GA算法，提出了一種應(yīng)用于微震源定位的反演方法。該算法充分利用BP神經(jīng)網(wǎng)絡(luò)的擬合能力和遺傳算法的全局尋優(yōu)能力，為進(jìn)一步豐富和完善地球物理非線性全局智能優(yōu)化反演技術(shù)，提供了理論與技術(shù)參考依據(jù)。
1 微震源反演方法概述
1.1 微震源定位原理
    定位反演技術(shù)基于聲發(fā)射學(xué)和地震學(xué)?；驹頌椋旱叵聨r石由人為或自然因素發(fā)生破裂及移動(dòng)時(shí)，產(chǎn)生地震波向四周傳播。通過布置于周圍空間的檢波器組采集實(shí)時(shí)數(shù)據(jù)，經(jīng)相關(guān)算法處理并在三維空間顯示微震源位置。常規(guī)方法[2]有蒙特卡洛法(MCM)、模擬退火法(SA)、人工神經(jīng)元法(ANN)及雙重殘差法(DDA)等。假設(shè)地層模型為均勻速度模型，監(jiān)測網(wǎng)絡(luò)中設(shè)置若干個(gè)不同位置檢波器進(jìn)行P波拾取[3]，如圖1所示。
    以到時(shí)殘差最小為目標(biāo)求解模型，基于地震學(xué)聯(lián)合反演思想求震源解時(shí)，走時(shí)模型為：
  

2 BP神經(jīng)網(wǎng)絡(luò)與遺傳混合算法
    近年來,BP神經(jīng)網(wǎng)絡(luò)與遺傳算法是地球物理勘探領(lǐng)域應(yīng)用最為廣泛的完全非線性反演方法。遺傳算法根據(jù)自然界生物優(yōu)勝劣汰的進(jìn)化原則，有機(jī)結(jié)合定向和隨機(jī)搜索，不斷逼近最優(yōu)解,并廣泛應(yīng)用于非線性規(guī)劃問題、組合優(yōu)化問題、參數(shù)辨識(shí)、控制器優(yōu)化等領(lǐng)域[4]。BP算法是用大量計(jì)算單元（即神經(jīng)元）構(gòu)成非線性系統(tǒng)，模仿人腦神經(jīng)系統(tǒng)的處理信息、存儲(chǔ)信息及檢索功能，常用于函數(shù)擬合與樣本分類等領(lǐng)域。
    結(jié)合非線性全局優(yōu)化方法應(yīng)用時(shí)的局限性，諸如BP神經(jīng)網(wǎng)絡(luò)存在易陷入到局部最小值、訓(xùn)練時(shí)間過長及遺傳算法的早熟現(xiàn)象普遍等問題，本文研究了利用神經(jīng)網(wǎng)絡(luò)的非線性擬合能力和遺傳算法的非線性尋優(yōu)能力，求取目標(biāo)函數(shù)全局極值的BP-GA算法，如圖3所示。

3 微震源反演模型
　　結(jié)合微震監(jiān)測工程，本文采用多臺(tái)站井地聯(lián)合監(jiān)測進(jìn)行微震源定位反演。速度模型采用水平層狀均勻介質(zhì)模型,分別為V1=2 500 m/s,V2=3000 m/s,V3=3 500 m/s。假定微震源空間位置坐標(biāo)為（77，28，-98），地面檢波器位置范圍為（188，148，0）~（248，98，0），井下檢波器位置范圍為（248，58，-5）~（248，58，-96）。檢波器屬性信息如表1所示。

4 算法仿真與結(jié)果分析

    反演算法總體思路為：根據(jù)反演問題目標(biāo)函數(shù)的特點(diǎn)構(gòu)建自適應(yīng)BP神經(jīng)網(wǎng)絡(luò),以目標(biāo)函數(shù)的輸入輸出數(shù)據(jù)訓(xùn)練網(wǎng)絡(luò)后，可以在適宜的權(quán)值和閾值約束下預(yù)測非線性函數(shù)輸出。遺傳算法的極值尋優(yōu)將訓(xùn)練后的網(wǎng)絡(luò)預(yù)測結(jié)果作為個(gè)體適應(yīng)度值，經(jīng)過選擇、交叉和變異等操作，尋找目標(biāo)函數(shù)的全局最優(yōu)值及對應(yīng)輸入值。
4.1 BP神經(jīng)網(wǎng)絡(luò)算法實(shí)現(xiàn)
    BP網(wǎng)絡(luò)采用三層拓?fù)浣Y(jié)構(gòu)[5]，分別為輸入層n1、隱含層n2和輸出層n3。利用正演模型得出的目標(biāo)函數(shù)100組輸入輸出數(shù)據(jù)來訓(xùn)練網(wǎng)絡(luò)。樣本的3個(gè)輸入?yún)?shù)，為微震源點(diǎn)空間坐標(biāo)(X0、Y0、Z0)，輸出參數(shù)即為目標(biāo)函數(shù)輸出值。隱含層數(shù)結(jié)合網(wǎng)絡(luò)精度與模型復(fù)雜程度選擇為單隱含層。網(wǎng)絡(luò)結(jié)構(gòu)為3-7-1，共3×7+7×1=28個(gè)權(quán)值，7+1=8個(gè)閾值,遺傳算法優(yōu)化參數(shù)的個(gè)數(shù)為28+8=36。
4.2 遺傳算法實(shí)現(xiàn)
    個(gè)體采用實(shí)數(shù)編碼，由于目標(biāo)函數(shù)只有3個(gè)輸入?yún)?shù)，所以個(gè)體長度為3。個(gè)體適應(yīng)值為BP神經(jīng)網(wǎng)絡(luò)預(yù)測值，適應(yīng)度值越小，則個(gè)體越優(yōu)，交叉概率和變異概率分別設(shè)為0.4和0.2。
4.3 結(jié)果分析
    微震源定位反演以實(shí)現(xiàn)通過地面與井下檢波器組拾取的初至波走時(shí)差來確定微震源的空間坐標(biāo)為目的。圖4給出了混合算法的Matlab仿真結(jié)果。由BP神經(jīng)網(wǎng)絡(luò)預(yù)測誤差可以看出BP神經(jīng)網(wǎng)絡(luò)擬合性能較好，符合定位反演誤差精度要求。遺傳算法適應(yīng)度曲線中適應(yīng)度函數(shù)值的漸變體現(xiàn)了較快逼近目標(biāo)函數(shù)最優(yōu)值的趨勢。

    本文基于遺傳與BP神經(jīng)網(wǎng)絡(luò)兩種非線性方法的優(yōu)越性及局限性，建立了相互滲透、配合的BP-GA混合算法模型，并應(yīng)用于微震源定位反演,改善了傳統(tǒng)非線性全局尋優(yōu)方法的搜索性能，提高了微震參數(shù)反演中的運(yùn)算效率和反演精度[6]。非線性全局尋優(yōu)算法能夠解決不確定性、高度非線性、超大規(guī)模的復(fù)雜問題，所具備的通用性、魯棒性、自適應(yīng)等特征，不僅在地球物理反演領(lǐng)域具有良好的應(yīng)用效果，而且在組合優(yōu)化、自動(dòng)控制、圖像處理、人工生命、管理決策等領(lǐng)域都將得到廣泛的應(yīng)用。
參考文獻(xiàn)
[1] Wei Chao, Li Xiaogan, Zheng Xiaofong. The group searchbased parallel algorithm for the serial monte carlo inversion  method[J].Applied Geophysics,2010,7(2):127-134.
[2] 田玥，陳曉非.地震定位研究綜述[J]. 地球物理學(xué)進(jìn)展，2002，17(1):147-155.
[3] 鞏思園，竇林名，馬小平，等.煤礦礦震定位中異向波速模型的構(gòu)建與求解[J].地球物理學(xué)報(bào)，2012，55(5):1757-1763.
[4] 陳霄，王寧.基于混沌DNA遺傳算法的模糊遞歸神經(jīng)網(wǎng)絡(luò)建模[J].控制理論與應(yīng)用，2011，28(11):1589-1594.
[5] 楊立強(qiáng)，宋海斌，郝天姚.基于BP 神經(jīng)網(wǎng)絡(luò)的波阻抗反演及應(yīng)用[J].地球物理學(xué)進(jìn)展,2005,20(1):34-37.
[6] 張新兵，王家林，吳健生. 混合最優(yōu)化算法在地球物理學(xué)中的應(yīng)用現(xiàn)狀與前景[J].地球物理學(xué)進(jìn)展，2003，18(2):5-6.