【摘要】：偏移成像是地震勘探中關(guān)鍵環(huán)節(jié),逆時偏移成像技術(shù)能夠適應(yīng)復(fù)雜速度模型,成像結(jié)果更為精確,成為當前地震勘探研究熱點之一。地震勘探中得到的地震記錄常存在介質(zhì)黏彈性導(dǎo)致的地震波衰減影響,受強烈衰減的記錄中能夠精細刻畫深層地區(qū)的信息模糊不清,使得勘探分辨率降低,更接近實際介質(zhì)假設(shè)的偏移方法越來越受到關(guān)注。黏彈性模型能夠更準確表征實際地下介質(zhì)地震波傳播過程,考慮黏彈介質(zhì)吸收衰減規(guī)律和波場傳播特征的逆時偏移方法能夠指導(dǎo)地震勘探資料處理解釋,幫助解決油氣勘探面臨的復(fù)雜問題,以提高勘探精度�；诖苏撐拈_展了基于Q補償?shù)酿椊橘|(zhì)逆時偏移成像研究,在逆時偏移波場延拓過程中對黏彈介質(zhì)的衰減影響進行處理,準確的補償衰減和校正相位畸變效應(yīng),得到更高分辨率的偏移成像結(jié)果。主要研究內(nèi)容和成果如下:1)論文從地震勘探頻帶內(nèi)較為準確的常Q黏彈模型出發(fā),推導(dǎo)了分數(shù)階時間偏導(dǎo)黏彈介質(zhì)波動方程,利用拉普拉斯算子將分數(shù)階時間偏導(dǎo)轉(zhuǎn)換至空間偏導(dǎo),得到了偽譜法常Q模型黏彈介質(zhì)拉普拉斯分數(shù)階波場延拓算子,實現(xiàn)了偽譜法黏彈介質(zhì)波動方程正演數(shù)值模擬。基于常Q模型的黏彈波動方程能夠較好模擬地震波黏彈介質(zhì)中傳播相位延遲振幅減弱的衰減現(xiàn)象,受強烈黏彈性衰減影響的地震記錄偏移時存在多分量波場混合和有效反射波受衰減信號較弱的問題。2)論文研究了彈性波逆時偏移成像方法以應(yīng)對黏彈地震記錄數(shù)據(jù)多分量波場混和問題。利用波場分離,得到了縱橫波分離的彈性波逆時偏移成像結(jié)果。為壓制逆時偏移過程中低頻噪音,發(fā)展了一階速度應(yīng)力能量范數(shù)成像條件的彈性波逆時偏移方法,引入了中間應(yīng)變分量,發(fā)展了一階速度應(yīng)力方程能量范數(shù)成像條件逆時偏移波場分量更新方式�；谀芰糠稊�(shù)的彈性波逆時偏移能夠有效壓制低頻干擾得到較高分辨率彈性波逆時偏移結(jié)果。3)論文發(fā)展了補償黏彈逆時偏移波場延拓算子構(gòu)建思路,改造了相位與振幅衰減解耦的黏彈介質(zhì)波動方程,重新構(gòu)建了相位振幅變化類波場延拓算子。重構(gòu)的相位校正類黏彈算子波場延拓過程能夠獨立校正相位不衰減振幅,振幅補償類黏彈算子模擬波傳播過程中能夠獨立補償振幅不改變相位。相位延遲校正和振幅衰減補償算子可以相互組合,構(gòu)成了靈活變化黏彈波場延拓技術(shù)序列。4)論文給出了基于Q補償黏彈逆時偏移框架,以解決黏彈地震記錄偏移需要校正相位振幅的問題。發(fā)展了補償?shù)酿椖鏁r偏移方法,提出了互相成像條件類和震源歸一化類成像條件補償黏彈逆時偏移策略�；ハ喑上駰l件類選用相位衰減振幅波場黏彈算子,震源歸一化類成像條件采用僅相位延遲和相位延遲振幅衰減黏彈算子,實現(xiàn)了補償黏彈逆時偏移成像算法;發(fā)展了縱橫波波場分離的補償黏彈逆時偏移技術(shù)。補償黏彈逆時偏移成像算法在正反傳波場延拓中進行相位校正和振幅補償,能夠得到更高分辨率偏移結(jié)果。論文研究過程中主要創(chuàng)新之處為:1)發(fā)展了一階速度應(yīng)力方程能量范數(shù)成像條件彈性波逆時偏移方法。該方法能夠有效壓制彈性波逆時偏移過程中低頻噪音,取得較好成像效果。2)構(gòu)造了相位與振幅變化的黏彈波動方程波場延拓算子。構(gòu)建的算子能夠獨立用于補償黏彈逆時偏移正反傳波場延拓中相位與振幅的解耦校正,解決了普通波場延拓算子難以進行黏彈波場中相位振幅校正補償?shù)膯栴}。3)發(fā)展了基于Q補償?shù)酿椊橘|(zhì)波動方程逆時偏移方法。該方法同時考慮黏彈波場多分量復(fù)雜性和黏滯性衰減效應(yīng),在波場延拓過程中對地震波相位和振幅準確校正補償,減小因黏彈介質(zhì)衰減效應(yīng)導(dǎo)致偏移結(jié)果不準確的影響,得到更高分辨率偏移結(jié)果。
[Abstract]:Migration imaging is a key link in seismic exploration. Reverse-time migration imaging technology can adapt to complex velocity model and its imaging results are more accurate. It has become one of the hot topics in seismic exploration. Viscoelastic model can more accurately characterize the seismic wave propagation process in actual underground media, and reverse-time migration method considering the absorption and attenuation law of viscoelastic media and wave field propagation characteristics can guide. Seismic exploration data processing and interpretation help to solve the complex problems faced by oil and gas exploration and improve exploration accuracy. Based on this paper, the inverse time migration imaging of viscoelastic media based on Q compensation is studied. The attenuation effect of viscoelastic media is processed in the process of inverse time migration wave field continuation, and the attenuation and phase distortion are corrected accurately. The main research contents and achievements are as follows: 1) Based on the more accurate constant Q viscoelastic model in seismic exploration frequency band, the wave equation of fractional order time-derivative viscoelastic medium is deduced, and the fractional order time-derivative is converted to spatial derivative by Laplace operator, and the pseudospectral constant Q is obtained. The forward numerical simulation of wave equation in viscoelastic media using pseudo-spectral method is realized by using Laplace fractional wave field continuation operator in the model.The viscoelastic wave equation based on constant Q model can simulate the attenuation of phase delay amplitude in seismic viscoelastic media and the migration of seismic records affected by strong viscoelastic attenuation. The problem of multi-component wavefield mixing and weak attenuated signal of effective reflection wave is existed. 2) The elastic wave inverse migration imaging method is studied to deal with the problem of multi-component wavefield mixing of viscoelastic seismic data. In low-frequency noise, an elastic wave inverse migration method based on first-order velocity-stress energy norm imaging condition is developed, and the intermediate strain component is introduced. A new method for updating the wave field components of first-order velocity-stress equation inverse migration based on energy norm imaging condition is developed. High-resolution elastic wave inverse migration results. 3) The idea of constructing wave field continuation operator for compensating viscoelastic inverse migration is developed in this paper. The wave equation of viscoelastic medium decoupled by phase and amplitude attenuation is modified, and the phase-amplitude variation wave-like field continuation operator is reconstructed. The reconstructed phase-corrected wave field continuation process of viscoelastic operator can be independent. The phase delay correction and amplitude attenuation compensation operators can be combined to form a flexible viscoelastic wave field continuation technique sequence. 4) A Q-compensated viscoelastic inverse time migration framework is proposed to solve the problem. Viscoelastic seismic migration requires correction of phase amplitude.A compensated viscoelastic inverse-time migration method is developed.A strategy of compensating viscoelastic inverse-time migration based on mutual imaging conditions and source normalization imaging conditions is proposed.The phase attenuation amplitude wave field viscoelastic operator is used for mutual imaging conditions and source normalization imaging conditions are used only. The phase delay and phase delay amplitude attenuation viscoelastic operator are used to compensate viscoelastic inverse time migration imaging algorithm, and the compensation viscoelastic inverse time migration technique for P-S wave field separation is developed. The main innovations in this paper are as follows: 1) An energy norm imaging conditional elastic wave inverse migration method for the first-order velocity-stress equation is developed, which can suppress the low-frequency noise effectively during the inverse migration of elastic wave and obtain better imaging effect. 2) The wave field continuation operator for the viscoelastic wave equation with phase and amplitude variations is constructed. The constructed operator can be used independently to compensate for the decoupling correction of phase and amplitude in forward and backward propagation wave field continuation of viscoelastic inverse time migration. The problem that ordinary wave field continuation operators are difficult to compensate for phase and amplitude in viscoelastic wave field is solved. The multi-component complexity and viscous attenuation effect of viscoelastic wave field can correct and compensate the phase and amplitude of seismic wave accurately in the process of wave field continuation, reduce the influence of inaccurate migration result caused by the attenuation effect of viscoelastic medium, and obtain higher resolution migration result.
【學位授予單位】：中國地質(zhì)大學
【學位級別】：博士
【學位授予年份】：2017
【分類號】：P631.4

【相似文獻】

相關(guān)期刊論文 前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];地球物理學進展;2012年04期

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

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

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

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

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

相關(guān)會議論文 前10條

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

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

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

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

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

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

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

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

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

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

相關(guān)重要報紙文章 前2條

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

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

相關(guān)博士學位論文 前10條

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

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

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

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

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

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

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

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

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

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

相關(guān)碩士學位論文 前10條

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

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

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

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

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

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

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

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

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

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

，

本文編號：2195432