【摘要】：快速精準的微震震源定位方法是微震監(jiān)測技術更好地發(fā)揮巖體穩(wěn)定性預測預警作用的基礎。精確的初至波走時計算是提高微震定位精度的關鍵。針對實際工程中帶有空洞和速度分區(qū)的復雜波速巖體,引入速度快、精度高和對復雜模型適應性強的多模板快速行進法multi-stencils fast marching methods(MSFM)用于初至波走時計算。該方法利用坐標旋轉的方式生成新的計算模板,使網格對角鄰點參與計算,提高了the Fast Marching Method(FMM)在對角方向的計算精度。首先,對比分析了一、二階FMM和MSFM計算精度和效率,表明二階MSFM具有更高的計算精度;其次,在分層速度巖體和帶空洞巖體中分別采用單一速度模型和二階MSFM計算初至波走時,通過與解析解對比發(fā)現,二階MSFM較單一速度模型的計算絕對誤差平均減小了97.65%和95.18%;然后,建立速度分層且?guī)в兴矶吹膸r體模型,驗證了二階MSFM算法對復雜速度模型適應性極強的特點;最后,提出了基于MSFM的復雜速度巖體微震定位算法,并將其應用到白鶴灘水電站左岸邊坡微震定位。通過對現場采集的4個爆破事件進行定位試算,得到使定位誤差平均值最小(9.6 m)的相對最優(yōu)層速度組合并用于考慮空洞的速度模型構建,爆破事件走時正演驗證了本速度模型的可靠性。對2015年5月的128個微震事件分別采用單一速度模型和提出的定位算法進行定位,通過對比定位效果,發(fā)現后者定位的事件較前者在巖體主要損傷區(qū)(層內錯動帶LS331和LS337)附近有更明顯呈條帶狀分布的聚集現象。研究表明MSFM算法作為走時正演方法具有很好的應用價值,提出的基于MSFM的微震定位方法能夠有效地提高震源定位精度。
【圖文】：

第36卷第2期郭亮等：基于MSFM的復雜速度巖體微震定位研究397MSFM算法通過坐標旋轉的方式使對角鄰點參與計算。如圖1[31]所示，點p1，p2，q1，q2為計算點x的4個對角鄰點，對角方向上的走時方向導數U=[U1U2]T表示為URT(x)(6a)圖1考慮對角鄰點的計算模板示意圖[31]Fig.1Thediagramofcalculationstencilconsideringthediagonalpoints[31]其中，T12Rrr，T()()()TTxyxxTx(6b)式中：r1，r2為對角線方向單位向量，r1=[r11r12]，r2=[r21r22]。經過變換有2TT121()()VTxURRU(7)由于有T111cos()sincos1RR(8)式中：為2個方向向量的夾角。那么：221122sinU2UUcosUV(9)式中：V為介質波速值。按照一階差分形式來近似方向導數，則有()max0，vvvTxTUxx(v=1，2)(10)式中：xv為已知節(jié)點空間坐標，Tv為其走時。通過求解方程式(9)可得到計算點x的走時T(x)。高階差分格式的MSFM求解公式本文不再贅述。將上述計算模板與FMM算法的基本模板聯(lián)合并在原FMM框架下實現[31-33]，便形成多模板快速行進法MSFM。2.2算法精度和效率分析本文建立模型為500m×500m的二維均勻模型。以模型某一角點為原點，方形的2個邊為X，Y軸建立坐標系，分別按照間距1，2，5m大小劃分3種網格，模型波速值為4000m/s。假定震源發(fā)生位置坐標為(0，0)，分別采用一階、二階的FMM及MSFM算法計算震源位置到其余各個節(jié)點的(初至)走時。采用幾何學中距離與速度之商作為走時解析解，將各節(jié)點走時計算值與解析解之差絕對值的平均值ξ作為精度評價指標，，ξ越大表示精度越

第36卷第2期郭亮等：基于MSFM的復雜速度巖體微震定位研究397MSFM算法通過坐標旋轉的方式使對角鄰點參與計算。如圖1[31]所示，點p1，p2，q1，q2為計算點x的4個對角鄰點，對角方向上的走時方向導數U=[U1U2]T表示為URT(x)(6a)圖1考慮對角鄰點的計算模板示意圖[31]Fig.1Thediagramofcalculationstencilconsideringthediagonalpoints[31]其中，T12Rrr，T()()()TTxyxxTx(6b)式中：r1，r2為對角線方向單位向量，r1=[r11r12]，r2=[r21r22]。經過變換有2TT121()()VTxURRU(7)由于有T111cos()sincos1RR(8)式中：為2個方向向量的夾角。那么：221122sinU2UUcosUV(9)式中：V為介質波速值。按照一階差分形式來近似方向導數，則有()max0，vvvTxTUxx(v=1，2)(10)式中：xv為已知節(jié)點空間坐標，Tv為其走時。通過求解方程式(9)可得到計算點x的走時T(x)。高階差分格式的MSFM求解公式本文不再贅述。將上述計算模板與FMM算法的基本模板聯(lián)合并在原FMM框架下實現[31-33]，便形成多模板快速行進法MSFM。2.2算法精度和效率分析本文建立模型為500m×500m的二維均勻模型。以模型某一角點為原點，方形的2個邊為X，Y軸建立坐標系，分別按照間距1，2，5m大小劃分3種網格，模型波速值為4000m/s。假定震源發(fā)生位置坐標為(0，0)，分別采用一階、二階的FMM及MSFM算法計算震源位置到其余各個節(jié)點的(初至)走時。采用幾何學中距離與速度之商作為走時解析解，將各節(jié)點走時計算值與解析解之差絕對值的平均值ξ作為精度評價指標，ξ越大表示精度越
【作者單位】： 四川大學水力學與山區(qū)河流開發(fā)保護國家重點實驗室;中國長江三峽集團公司;
【基金】：國家自然科學基金資助項目(51679158,51374149) 國家重點基礎研究發(fā)展計劃(973)項目(2015CB057903)~~
【分類號】：P631.4

【相似文獻】

相關期刊論文 前10條

1 魏利娜;薛世鵬;周春艷;丁慧哲;李靜;;微震信號現場監(jiān)測試驗及特征研究[J];河北工程大學學報(自然科學版);2014年02期

2 ;文摘[J];世界地震工程;1985年03期

3 周立功;微震臺網在地質工作中的應用[J];中國地質;1983年01期

4 魯振華,張連城;門頭溝礦微震的近場監(jiān)測效能評估[J];地震;1989年05期

5 唐禮忠;張君;李夕兵;汪令輝;周建雄;劉濤;;基于定量地震學的礦山微震活動對開采速率的響應特性研究[J];巖石力學與工程學報;2012年07期

6 王鵬;常旭;王一博;王璐琛;翟鴻宇;;基于時頻稀疏性分析法的低信噪比微震事件識別與恢復[J];地球物理學報;2014年08期

7 王鵬;桂志先;謝宋雷;王朝鋒;張建萍;;壓裂微震數據的快速讀取及可視化[J];物探化探計算技術;2010年02期

8 肖和平;湖南邵陽市郊微震活動特征及成因[J];華南地震;1989年03期

9 韓伯文;關于地微震若干問題的看法與認識[J];勘察科學技術;2001年06期

10 王昌銀;董曉偉;;淺論微震系統(tǒng)與大孔徑卸壓中的相互應用[J];科技風;2012年16期

相關會議論文 前7條

1 李國營;劉虎;;唐口礦井受采動影響下微震活動規(guī)律研究[A];綜采放頂煤技術理論與實踐的創(chuàng)新發(fā)展——綜放開采30周年科技論文集[C];2012年

2 王永仁;歐陽振華;孔令海;趙善坤;李曉璐;;不同地質條件下微震信號衰減規(guī)律研究[A];第3屆全國工程安全與防護學術會議論文集[C];2012年

3 高峰;呂進國;;基于HHT方法的微震信號時頻分析[A];北京力學會第18屆學術年會論文集[C];2012年

4 王鵬;常旭;王一博;王璐琛;;水力壓裂誘發(fā)微震事件的頻譜特征統(tǒng)計分析[A];中國地球物理2013——第二十三專題論文集[C];2013年

5 張傳玖;毛學健;;強沖擊煤層綜采工作面微震活動分析研究[A];煤炭開采新理論與新技術——中國煤炭學會開采專業(yè)委員會2012年學術年會論文集[C];2012年

6 周勝健;劉超;;采動下新莊孜礦62113工作面地質預警微震實踐[A];2010年度淮南礦業(yè)集團煤炭學會學術交流會論文匯編（一）·地質專業(yè)[C];2010年

7 李學政;唐云凱;;微震工程檢測與應用技術[A];國家安全地球物理叢書（七）——地球物理與核探測[C];2011年

相關重要報紙文章 前1條

1 本報記者　張靜;三峽微震[N];21世紀經濟報道;2003年

相關博士學位論文 前7條

1 黃志平;深埋隧洞開挖卸荷巖爆孕育過程及微震預警分析[D];東北大學;2015年

2 于洋;深埋隧洞即時型巖爆孕育過程的微震信息特征分析及分形研究[D];東北大學;2014年

3 楊威;煤巖變形破裂電磁和微震信號關聯(lián)響應機理及特征研究[D];中國礦業(yè)大學（北京）;2014年

4 劉建坡;深井礦山地壓活動與微震時空演化關系研究[D];東北大學;2011年

5 胡永泉;地面微震資料去噪方法研究[D];西南石油大學;2013年

6 劉超;采動煤巖瓦斯動力災害致災機理及微震預警方法研究[D];大連理工大學;2011年

7 葉慶東;大別蘇魯地區(qū)背景噪聲成像與汶川地震科學鉆探井孔附近微震定位[D];中國地震局地球物理研究所;2014年

相關碩士學位論文 前10條

1 周莉;基于礦山巖石破裂的微震信號的檢測與研究[D];昆明理工大學;2015年

2 商鵬;爆破掘進煤巖破裂微震信號特征研究[D];河南理工大學;2014年

3 趙中流;基于Fermat原理和加速度估計的微震定位研究[D];吉林大學;2016年

4 李保林;煤礦微震與爆破信號特征提取及識別研究[D];中國礦業(yè)大學;2016年

5 王強;巖層破裂微震動同步數據采集網絡系統(tǒng)研究[D];江蘇科技大學;2016年

6 樊手;地下洞室圍巖與支護損傷演化的微震釋放研究[D];西南石油大學;2016年

7 朱超;微震實時在線監(jiān)測系統(tǒng)的研究與實現[D];武漢科技大學;2012年

8 景林波;煤礦微震信號特征及傳播規(guī)律研究[D];中國礦業(yè)大學;2014年

9 劉佳;地下礦山實時信號的微震震源分析[D];武漢科技大學;2013年

10 閆憲磊;東灘煤礦深井綜放工作面微震發(fā)生規(guī)律研究[D];遼寧工程技術大學;2012年

本文編號：2551844