發(fā)布時間：2017-12-28 22:20

  本文關(guān)鍵詞：隧道空間有限元地震波場模擬與逆時偏移成像　出處：《中國地質(zhì)大學(xué)》2017年博士論文　論文類型：學(xué)位論文

  更多相關(guān)文章： 隧道地震超前預(yù)報 三維有限元 波場模擬 邊界條件 彈性波逆時偏移

【摘要】：隧道工程作為交通建設(shè)的咽喉,其安全和高效地施工一直受到人們的重視。但隧道圍巖地質(zhì)情況復(fù)雜多變,施工過程中斷層破碎帶、巖溶等不良地質(zhì)體極易造成突水突泥、圍巖失穩(wěn)等地質(zhì)災(zāi)害。因此,開展隧道超前預(yù)報的工作,準確掌握隧道掌子面前方的地質(zhì)情況,并提前做好預(yù)防工作,是隧道工程安全和高效施工的重要保障。隧道超前預(yù)報技術(shù)主要是從20世紀70年代開始興起的,其中地震反射波法超前預(yù)報技術(shù)以其精度高,探測距離遠,而備受工程界青睞。但隧道復(fù)雜的圍巖地質(zhì)情況使得地震波場異常復(fù)雜,難以識別有效波場和獲取較可靠的成像結(jié)果。目前,針對隧道超前預(yù)報的研究主要是從波場正演模擬與成像兩方面進行的,但二者仍然在以下幾個方面值得進一步的研究和探討:1)對于隧道正演模擬來說,通常采用的方法是有限差分,有限差分在處理隧道這種特殊的弧形自由邊界時往往比較困難,不利于獲得符合實際的隧道空間地震波場響應(yīng)特征;2)對于隧道超前預(yù)報偏移成像方法的研究大多是采用Kirchhoff偏移或是繞射掃描偏移,而對目前公認的成像精度最高的逆時偏移也大多是基于二維空間,缺乏針對隧道地震超前預(yù)報三維逆時偏移成像的研究以及有利于成像的觀測系統(tǒng)討論;3)目前針對隧道地震超前預(yù)報中獲得的多分量地震資料,主要的處理流程是基于聲波偏移或是縱橫波分離式的標量偏移。由于隧道地震資料的復(fù)雜性,往往很難通過對地震記錄進行波場分離獲得純縱波或是純橫波資料,在偏移成像時,殘留的能量容易對成像造成干擾。針對目前隧道超前預(yù)報研究中值得進一步探討的問題,本文開展了基于有限元法的隧道三維空間彈性波波場模擬和基于彈性波逆時外推的隧道三維縱波逆時偏移成像研究。主要工作及成果如下:(1)基于二階位移彈性波波動方程,利用虛功原理推導(dǎo)并建立了三維彈性波波場模擬的有限元方程,并介紹了具有二階精度的顯式Newmark差分算法的實現(xiàn)過程;對比分析了行壓縮存儲(CSR)以及一維變帶寬存儲對稀疏矩陣的存儲效率;在給出了非結(jié)構(gòu)網(wǎng)格下有限元爆炸震源的加載方式下,實現(xiàn)了三維有限元彈性波全波場模擬,為后文開展逆時偏移成像研究提供了基礎(chǔ)。(2)基于完美匹配層理論,詳細推導(dǎo)并建立了基于二階位移波動方程的分裂完美匹配層(PML)邊界條件和非分裂卷積完美匹配層(CFS-UPML)邊界條件的有限元方程,通過與解析解對比驗證了本文有限元模擬面波的正確性,對比分析了兩種邊界條件對面波、掠射角入射波和近地表高泊松比介質(zhì)下的吸收效果,以及長時間計算的穩(wěn)定性情況;最后對有限元在三維傾斜自由表面模擬面波的正確性進行了驗證。研究結(jié)果表明:具有網(wǎng)格剖分靈活以及自動滿足自由邊界條件等優(yōu)勢的有限元法能夠獲得較精確的面波記錄,并且適應(yīng)于復(fù)雜環(huán)境下的波場模擬;CFS-UPML邊界條件在對掠射角入射波的吸收效果以及高泊松比介質(zhì)中的長時間穩(wěn)定性方面更具優(yōu)勢,但計算效率要低于PML邊界。(3)基于隧道施工過程中常見的不良地質(zhì)體(巖性分界面、溶洞)的結(jié)構(gòu)特點,構(gòu)建了三維隧道地質(zhì)地球物理模型,實現(xiàn)了隧道波場的三維彈性波全波場模擬,并進行了隧道波場特征分析,增強了對隧道波場傳播特征的系統(tǒng)認識,為隧道地震反射信息的識別提供了理論依據(jù)。分析結(jié)果顯示:炮點位于隧道邊墻時,巖性分界面產(chǎn)生的反射橫波能量相對來說要強于反射縱波;檢波器布置在靠近傾斜界面一側(cè)的隧道側(cè)壁上時,獲得的反射縱波能量強于另外一側(cè),更有利于反射縱波成像;溶洞充填空氣時,反射波的能量最強,其次是充填水,而充填軟弱土產(chǎn)生的反射波能量最弱。(4)在分析了彈性波逆時偏移成像基本理論的基礎(chǔ)上,對逆時偏移成像中的關(guān)鍵步驟(波場分離、成像條件以及成像噪音壓制)進行了詳細的分析和討論,并對其在隧道逆時偏移成像中的適用性進行了驗證。結(jié)果表明:采用基于Helmholtz分解的方法能夠有效的對隧道波場進行縱橫波解耦,以及邊界存儲策略的互相關(guān)成像條件的實施方式,能夠有效解決三維逆時偏移的內(nèi)存消耗問題;基于Poynting矢量和Laplace濾波的方法能夠有效消除隧道逆時偏移成像中的低頻噪音干擾。通過上述討論,最終建立了基于多分量地震資料-彈性波波場外推-縱橫波分離-彈性縱波成像的隧道地震超前預(yù)報三維彈性縱波逆時偏移成像的流程。(5)在隧道有限元三維彈性波波場模擬以及彈性波逆時偏移成像的研究基礎(chǔ)上,實現(xiàn)了基于彈性波逆時外推的隧道三維縱波逆時偏移成像。基于60°傾斜界面模型,對影響隧道成像結(jié)果的主要因素(觀測系統(tǒng)以及噪音)進行了討論,結(jié)果表明:基于彈性波逆時外推的三維彈性縱波成像能夠有效的進行隧道成像;雙邊觀測系統(tǒng)更有利于傾斜界面的成像,理論記錄成像結(jié)果表明采用一條排列即可滿足成像精度要求。噪聲討論結(jié)果表明:當(dāng)噪音比例小于等于50%時,獲得的成像剖面均較為可靠;隨著噪音比例的增加,成像結(jié)果中虛假同相軸的能量不斷增加,當(dāng)噪音比例大于70%時,成像結(jié)果才會受到虛假成像同相軸的干擾,可見,本文采用的隧道逆時偏移成像流程具有較好的抗噪性。(6)基于上述三維彈性縱波逆時偏移成像的流程,開展了隧道典型不良地質(zhì)體的理論模型(與隧道軸線成90°夾角)成像研究。結(jié)果表明:本文采用的基于多分量地震資料和彈性波波場外推的三維彈性縱波逆時偏移成像方法能夠有效的對掌子面正前方不良地質(zhì)體理論模型進行成像。通過與傾斜界面以及側(cè)面溶洞的成像結(jié)果對比可知,偏移剖面中與隧道軸線成90°夾角的地質(zhì)體成像能量相對較弱,因此,在高背景噪聲情況下,很難實現(xiàn)對與隧道軸線成90°夾角地質(zhì)體的精確成像。(7)基于本文討論的觀測系統(tǒng)以及三維彈性縱波成像方法進行了實際隧道地震超前預(yù)報資料的成像研究,通過對含有噪音記錄的成像,說明了本文采用方法具有一定的抗噪性;與二維偏移成像結(jié)果的對比,驗證了本文采用的成像方法的可行性和有效性。本文的創(chuàng)新性成果主要體現(xiàn)在:(1)基于二維分裂PML邊界條件,將其有效的擴展到了三維空間,并實現(xiàn)了基于二階位移方程的分裂PML邊界條件和CFS-UPML邊界條件的有限元隧道地震全波場模擬;(2)有效實現(xiàn)了三維隧道縱波逆時偏移成像,并對影響成像精度的關(guān)鍵參數(shù)進行了系統(tǒng)討論與分析,建立了滿足成像精度要求的觀測系統(tǒng),形成了基于多分量地震資料-彈性波波場外推-縱橫波分離-噪聲壓制-彈性縱波成像的隧道地震超前預(yù)報三維彈性縱波逆時偏移成像的流程。
[Abstract]:As the throat of traffic construction, the safety and efficient construction of tunnel engineering has been paid much attention by people. But the geological condition of tunnel surrounding rock is complex and changeable. During construction process, such unfavorable geological bodies as interruption zone, fracture zone and karst are very likely to cause geological disasters such as water bursting, mud outburst and surrounding rock instability. Therefore, to carry out the tunnel prediction work, accurately grasp the geological conditions of the tunnel face, and advance the prevention work, is an important guarantee for the safe and efficient construction of tunnel engineering. The advanced prediction technology of tunnel is mainly started in 1970s. The advanced prediction technology of seismic reflection wave is favored by engineering circles because of its high accuracy and long detection distance. However, the complex geological condition of the tunnel makes the seismic wave field very complicated, and it is difficult to identify the effective wave field and obtain more reliable imaging results. At present, the research of tunnel prediction is mainly from the wave field forward modeling and imaging of the two aspects, but the two still in the following aspects worthy of further research and discussion: 1) for the simulation of tunnel, is a commonly used method of finite difference, finite difference in the tunnel special arc free boundary is often difficult to get in line with the actual seismic wave field response characteristics of tunnel space; 2) for the prediction of migration method in tunnel mostly uses the Kirchhoff offset or diffraction scan migration, and the currently accepted imaging the highest accuracy of reverse time migration was also based on two-dimensional space in view of the lack of tunnel seismic prediction, 3D reverse time migration imaging research and observation system for imaging the discussion; 3) for the tunnel seismic prediction for multi points The main processing flow of seismic data is based on the wave offset or the scalar offset of the longitudinal and transverse wave separation. Because of the complexity of tunnel seismic data, it is often difficult to get pure longitudinal wave or pure shear wave data by separating wave field from seismic records. When migration imaging, residual energy can easily cause interference to imaging. Aiming at the problems that should be further explored in the study of tunnel advanced prediction, this paper develops three-dimensional elastic wave field simulation based on finite element method and 3D longitudinal wave reverse time migration imaging based on elastic wave inverse time extrapolation. The main work and results are as follows: (1) two order displacement based on the elastic wave equation is derived by using the principle of virtual work, and established a 3D elastic wave field simulation of finite element equations, and introduces explicit Newmark with two order accuracy difference algorithm implementation process; comparative analysis of the compressed storage (CSR) and one dimension variable bandwidth of sparse matrix storage efficiency; finite element is given under unstructured explosive source loading method, the three-dimensional finite element full elastic wave field simulation, for the development of reverse time migration imaging research provides a basis for. (2) based on the theory of perfectly matched layer, is derived and the matching layer two order displacement wave equation based on perfect splitting (PML) boundary conditions and non dividing convolutional perfectly matched layer (CFS-UPML) boundary conditions of finite element equations, by comparison with analytic solutions and verify the correctness of the finite element simulation of surface waves. Comparative analysis of two kinds of boundary conditions of surface wave, the grazing angle of incident wave and near surface absorption of medium high Poisson's ratio, and the stability of long time calculation; finally the finite element simulation of free surface in three-dimensional inclined surface wave are verified. The results show that: with the finite element method and flexible mesh automatically satisfy the free boundary conditions and other advantages can obtain more accurate records of surface wave, wave field simulation and adapt to the complex environment; CFS-UPML boundary conditions on the grazing angle of the incident wave absorption effect and high Poisson's ratio of long time stability of medium the advantages, but the computational efficiency is lower than the PML boundary. (3) common adverse geological tunnel construction process based on (lithologic interface, cave) structure characteristics, build the three-dimensional tunnel geological geophysical model, the 3D elastic wave tunnel wave field full wave field simulation, and analyzed the tunnel wave field characteristics, and enhance the understanding of wave propagation in tunnel the characteristics of the system, provides a theoretical basis for the identification of tunnel seismic reflection information. The analysis results show that the shot point is located in the tunnel wall, shear wave reflection energy generated by lithologic interface relatively stronger than the reflection wave; detector disposed near the tunnel sidewall tilt interface on one side, the other side is stronger than the P-wave reflection energy, more conducive to reflection wave imaging; cave filling the air, the strongest energy reflection wave, followed by water filling, filling and reflection of soft soil produced by wave energy. (4) in the analysis of the basic theory of elastic wave reverse time migration, the key steps in the reverse time migration imaging (wave field separation, imaging conditions and imaging noise suppression) are analyzed and discussed in detail, and the applicability of the offset in the tunnel in reverse time imaging is verified. The results show that the method based on Helmholtz decomposition of effective tunnel wave field wave decoupling, implementation method and boundary storage strategy of cross-correlation imaging condition, can effectively solve the problem of offset memory consumption inverse method; Poynting vector and Laplace filter can effectively eliminate the low-frequency noise migration in tunnel based on the inverse. Based on the above discussions, a three-dimensional elastic longitudinal wave reverse time migration imaging process based on multi-component seismic data elastic wave field extrapolation longitudinal and transverse wave separation elastic longitudinal wave imaging is finally established. (5) on the basis of the study of the three-dimensional elastic wave field simulation and the inverse time migration imaging of the elastic wave in the tunnel finite element, the tunnel based on the extrapolation of the elastic wave inverse time is realized.
【學(xué)位授予單位】：中國地質(zhì)大學(xué)
【學(xué)位級別】：博士
【學(xué)位授予年份】：2017
【分類號】：U452.11;P631.4

【相似文獻】

中國期刊全文數(shù)據(jù)庫 前10條

1 薛桂霞,鄧世坤,劉秀娟;逆時偏移在探地雷達信號處理中的應(yīng)用[J];煤田地質(zhì)與勘探;2004年01期

2 楊勤勇;段心標;;逆時偏移技術(shù)發(fā)展現(xiàn)狀與趨勢[J];石油物探;2010年01期

3 張春燕;李振春;孫小東;;逆時偏移方法技術(shù)進展綜述[J];勘探地球物理進展;2010年05期

4 滕厚華;許衛(wèi)華;于曉莉;尚新民;賀紅蕾;;無反射遞推算法疊后逆時偏移的研究與應(yīng)用[J];油氣地球物理;2011年04期

5 許璐;孟小紅;劉國峰;;逆時偏移去噪方法研究進展[J];地球物理學(xué)進展;2012年04期

6 王童奎;李瑩;郭愛華;石文武;楊曉利;孫明;潘歡;;逆時偏移技術(shù)在南堡1號構(gòu)造中的應(yīng)用研究[J];地球物理學(xué)進展;2012年06期

7 康智清;劉恩;司杰戈;陳康;;逆時偏移剖面假象產(chǎn)生機制及其消除方法[J];石油天然氣學(xué)報;2013年01期

8 胡明順;潘冬明;董守華;金紅娣;趙新;;基于散射理論的逆時偏移隨機邊界構(gòu)建策略優(yōu)化分析[J];地球物理學(xué)進展;2013年04期

9 張叔倫;陳靜;張關(guān)泉;;逆時偏移的方向?qū)?shù)算法[J];石油地球物理勘探;1991年06期

10 Ruey-Chyuan Shih;A.R.Levander;陳炳文;;層剝離逆時偏移[J];國外油氣勘探;1995年04期

中國重要會議論文全文數(shù)據(jù)庫 前10條

1 杜啟振;秦童;;逆時偏移最佳匹配層邊界實現(xiàn)[A];中國地球物理·2009[C];2009年

2 陳康;;逆時偏移剖面的波數(shù)域分析及高階拉普拉斯濾波算法的實現(xiàn)[A];“地球·資源”全國博士生學(xué)術(shù)論壇會議論文摘要集[C];2011年

3 龍云;韓立國;韓利;鄧武斌;;單程波逆時偏移[A];中國地球物理學(xué)會第二十七屆年會論文集[C];2011年

4 李慶洋;李振春;黃建平;;基于雙變算法的精細逆時偏移成像[A];中國地球物理2013——第二十二專題論文集[C];2013年

5 曹曉莉;黃建平;李振春;;最小二乘逆時偏移高精度成像方法研究[A];中國地球物理2013——第二十二專題論文集[C];2013年

6 許璐;孟小紅;劉國峰;;成像后壓制逆時偏移噪聲[A];中國地球物理學(xué)會第二十七屆年會論文集[C];2011年

7 薛霆哠;張書杰;;逆時偏移中坡印廷向量濾波的改進方法[A];中國地球物理2013——第二十二專題論文集[C];2013年

8 楊仁虎;毛偉健;凌云;;逆時偏移儲層成像波形校正研究[A];中國地球物理2013——第二十專題論文集[C];2013年

9 曾同生;李振春;;基于疊前逆時偏移的速度分析方法研究[A];中國地球物理學(xué)會第二十七屆年會論文集[C];2011年

10 胡昊;劉伊克;王一博;;地震疊前逆時偏移不同邊界處理方式研究[A];中國地球物理學(xué)會第二十七屆年會論文集[C];2011年

中國重要報紙全文數(shù)據(jù)庫 前2條

1 王躍;物探院逆時偏移成像技術(shù)躋身國際前列[N];中國石化報;2013年

2 本報記者 劉洪宇;石油HPC走向“平衡計算”[N];中國計算機報;2011年

中國博士學(xué)位論文全文數(shù)據(jù)庫 前10條

1 王京;隧道空間有限元地震波場模擬與逆時偏移成像[D];中國地質(zhì)大學(xué);2017年

2 李靜爽;基于近似解析離散化算子的逆時偏移方法及其應(yīng)用[D];清華大學(xué);2014年

3 程飛;多井井間三維波場數(shù)值模擬與逆時偏移成像研究[D];中國地質(zhì)大學(xué);2015年

4 馮蘭蘭;保辛幾何結(jié)構(gòu)不變的逆時偏移方法研究[D];清華大學(xué);2015年

5 楊佳佳;多分量地震波逆時偏移的關(guān)鍵技術(shù)研究[D];中國海洋大學(xué);2015年

6 高新成;基于云計算的逆時偏移數(shù)據(jù)處理方法研究[D];東北石油大學(xué);2016年

7 Naveed Akram;基于彈性波模擬的逆時偏移以及地震數(shù)據(jù)解釋[D];中國科學(xué)技術(shù)大學(xué);2017年

8 王娟;TTI介質(zhì)準縱波逆時偏移方法研究[D];中國石油大學(xué)(華東);2015年

9 趙旭;時域聲波方程逆時偏移成像及其角道集計算[D];吉林大學(xué);2015年

10 蔡志成;基于Q補償?shù)酿椊橘|(zhì)逆時偏移成像研究[D];中國地質(zhì)大學(xué);2017年

中國碩士學(xué)位論文全文數(shù)據(jù)庫 前10條

1 陳洪杰;基于聲波方程的數(shù)值模擬與逆時偏移方法研究[D];吉林大學(xué);2010年

2 魏巍;基于疊前逆時偏移方法的圓管損傷識別[D];哈爾濱工業(yè)大學(xué);2015年

3 張敬秋;交錯網(wǎng)格有限差分法地震資料逆時偏移[D];成都理工大學(xué);2015年

4 劉詩竹;VSP數(shù)據(jù)成像方法研究[D];東北石油大學(xué);2015年

5 柯璇;復(fù)雜構(gòu)造疊前逆時偏移成像方法研究[D];東北石油大學(xué);2015年

6 侯維麗;逆時偏移成像[D];東北石油大學(xué);2015年

7 郭雪豹;VSP地震資料逆時偏移方法研究[D];東北石油大學(xué);2015年

8 張千祥;純P波各向異性介質(zhì)波場數(shù)值模擬及逆時偏移研究[D];吉林大學(xué);2016年

9 于江龍;聲波逆時偏移及噪聲壓制[D];吉林大學(xué);2016年

10 方修政;各向異性介質(zhì)逆時偏移方法研究[D];東北石油大學(xué);2016年

，

本文編號：1347664