本文選題：全波形反演　切入點：互相關(guān)　出處：《吉林大學(xué)》2017年碩士論文

【摘要】：在地震勘探數(shù)據(jù)的處理與解釋中,速度是不可或缺的重要參數(shù),對于地下介質(zhì)構(gòu)造的描述和儲層預(yù)測都有重要作用。近年來油氣田的勘探開發(fā)程度逐日提升,勘探目標從簡單含油構(gòu)造轉(zhuǎn)向復(fù)雜構(gòu)造及隱蔽油氣藏,使得地震勘探對速度建模提出了更高的要求。傳統(tǒng)的速度建模方法,例如走時層析成像和偏移速度分析等,往往不能滿足復(fù)雜地質(zhì)條件下地震資料處理解釋所需的精度和分辨率。全波形反演是一種能夠充分利用疊前地震波場全部信息的速度建模方法,它使用最優(yōu)化方法對模擬記錄和觀測記錄進行匹配,使兩者之間的誤差泛函最小化,從而得到相應(yīng)的地下介質(zhì)速度模型。全波形反演在復(fù)雜地質(zhì)條件下也能較好地恢復(fù)出地下構(gòu)造的細節(jié),與傳統(tǒng)速度建模方法相比具有更高的分辨率和精度。本文深入研究了二維時間域全波形反演的基本理論,推導(dǎo)了聲波方程的規(guī)則網(wǎng)格高階有限差分格式和目標函數(shù)的梯度公式,同時對全波形反演的穩(wěn)定性進行了分析。由于常規(guī)全波形反演采用的2l范數(shù)形式的目標函數(shù)具有很強的非線性,加之模型參數(shù)變化復(fù)雜,使得反演過程不穩(wěn)定,在迭代時容易陷入局部極小值。為了使反演更好地收斂到全局極小值,通常需要提供一個足夠精確的初始速度模型或者需要地震資料包含較高信噪比的低頻信息。針對這一問題,本文提出了一種基于互相關(guān)的初始模型構(gòu)建方法,旨在為常規(guī)全波形反演提供一個較好的初始模型,以改善反演的穩(wěn)定性。該方法使用互相關(guān)目標函數(shù)測量模擬記錄與觀測記錄的相似性,側(cè)重地震記錄之間的相位匹配,具有更強的線性性質(zhì),因此更不容易陷入局部極小值�；诨ハ嚓P(guān)目標函數(shù)的全波形反演能夠恢復(fù)出模型的宏觀構(gòu)造,以此作為常規(guī)全波形反演的初始模型,可以有效提高最終反演結(jié)果的精度和穩(wěn)定性�？紤]到全波形反演的計算成本巨大,本文提出了一種時間域高效全波形反演方法,引入混合震源技術(shù)來提高計算效率。組成混合震源的各個獨立震源之間相互影響會產(chǎn)生嚴重的串?dāng)_噪聲,本文提出了一種動態(tài)隨機編碼方法來壓制串?dāng)_噪聲。該方法在每十次迭代后均重新隨機抽取多個震源參與反演,并對各震源的振幅和相位進行隨機編碼。模型試驗表明,基于動態(tài)隨機編碼的混合震源技術(shù)既能有效提高全波形反演的計算效率又能保證反演結(jié)果的質(zhì)量。為了驗證本文所提出的方法在地震數(shù)據(jù)缺少低頻成分時的反演效果,使用高通濾波器濾除地震子波的低頻成分,結(jié)果表明與常規(guī)方法相比,本文方法對低頻信息的依賴性更低。另外還進行了噪聲測試,進一步驗證了本文所提出的方法在地震數(shù)據(jù)含噪聲時也具有較好的反演效果和穩(wěn)定性。
[Abstract]:In the process and interpretation of seismic exploration data, velocity is an indispensable parameter, which plays an important role in the description of underground media structure and reservoir prediction. In recent years, the degree of exploration and development of oil and gas fields has been increasing day by day. From simple oil-bearing structures to complex structures and subtle reservoirs, seismic exploration has put forward higher requirements for velocity modeling. Traditional velocity modeling methods, such as travel time tomography and migration velocity analysis, etc. The precision and resolution of seismic data processing and interpretation under complex geological conditions are often not satisfied. Full waveform inversion is a velocity modeling method which can make full use of all the information of prestack seismic wave field. It uses the optimization method to match the simulated and observed records to minimize the error functional between them. Thus the corresponding velocity model of underground media can be obtained. Full wave inversion can also recover the details of underground structure in complex geological conditions. Compared with the traditional velocity modeling method, this paper has higher resolution and accuracy. In this paper, the basic theory of 2-D time-domain full waveform inversion is studied, and the regular grid high-order finite difference scheme of acoustic equation and the gradient formula of objective function are derived. At the same time, the stability of full waveform inversion is analyzed. The inversion process is unstable because of the strong nonlinearity of the objective function in the form of 2l norm and the complex variation of the model parameters. It is easy to fall into local minima during iteration. In order to make the inversion converge better to the global minimum, It is usually necessary to provide a sufficiently accurate initial velocity model or to contain low frequency information with high signal-to-noise ratio (SNR) in seismic data. To solve this problem, a cross-correlation based initial model construction method is proposed in this paper. In order to provide a better initial model for conventional full-waveform inversion and improve the stability of inversion, this method uses cross-correlation objective function to measure the similarity between simulated and observed records, with emphasis on phase matching between seismic records. The whole waveform inversion based on the cross-correlation objective function can restore the macroscopic structure of the model, which can be used as the initial model of the conventional full-waveform inversion. It can effectively improve the accuracy and stability of the final inversion results. Considering the huge computation cost of the full waveform inversion, a time domain efficient full waveform inversion method is proposed in this paper. The hybrid source technology is introduced to improve the computational efficiency. The interaction between the individual sources that make up the mixed source will produce serious crosstalk noise. In this paper, a dynamic random coding method is proposed to suppress crosstalk noise. After every ten iterations, multiple sources are selected randomly to participate in the inversion, and the amplitude and phase of each source are randomly coded. The hybrid source technique based on dynamic random coding can effectively improve the computation efficiency of full waveform inversion and guarantee the quality of inversion results. In order to verify the inversion effect of the method proposed in this paper, when the seismic data is short of low frequency component, High-pass filter is used to filter the low-frequency components of seismic wavelet. The results show that the method is less dependent on low-frequency information than the conventional method. It is further verified that the proposed method has good inversion effect and stability when seismic data contain noise.
【學(xué)位授予單位】：吉林大學(xué)
【學(xué)位級別】：碩士
【學(xué)位授予年份】：2017
【分類號】：P631.4

【相似文獻】

相關(guān)期刊論文 前10條

1 陳勇;李揚;陳小宏;;地震波形反演綜述[J];黑龍江大學(xué)自然科學(xué)學(xué)報;2012年03期

2 王薇;韓波;唐錦萍;;地震波形反演的稀疏約束正則化方法[J];地球物理學(xué)報;2013年01期

3 胡光輝;賈春梅;夏洪瑞;賀劍波;宋林;沈忠秋;;三維聲波全波形反演的實現(xiàn)與驗證[J];石油物探;2013年04期

4 C. Bunks;多尺度地震波形反演[J];石油物探譯叢;1996年05期

5 楊勤勇;胡光輝;王立歆;;全波形反演研究現(xiàn)狀及發(fā)展趨勢[J];石油物探;2014年01期

6 丁繼才;常旭;劉伊克;汪長永;;反射地震數(shù)據(jù)的逐層波形反演[J];地球物理學(xué)報;2007年02期

7 肖盈;薛明星;賀振華;黃德濟;;高斯—牛頓法與梯度法波形反演比較研究[J];內(nèi)蒙古石油化工;2008年18期

8 孫曉琳;高建;薛冰;;拉普拉斯域頻變衰減常數(shù)波形反演方法[J];山東科技大學(xué)學(xué)報(自然科學(xué)版);2011年01期

9 汪超;趙偉;高靜懷;;一種用于波形反演的改進差分進化算法[J];石油地球物理勘探;2012年02期

10 楊午陽;王西文;雍學(xué)善;陳啟燕;;地震全波形反演方法研究綜述[J];地球物理學(xué)進展;2013年02期

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

1 陳永芮;李振春;張凱;;一種改進的正則化——測井約束全波形反演方法[A];中國地球物理2013——第十九專題論文集[C];2013年

2 王一博;常旭;劉伊克;;多尺度波形反演方法[A];中國地球物理2010——中國地球物理學(xué)會第二十六屆年會、中國地震學(xué)會第十三次學(xué)術(shù)大會論文集[C];2010年

3 盧回憶;劉伊克;常旭;;波形反演多尺度方法研究[A];中國地球物理學(xué)會第二十七屆年會論文集[C];2011年

4 孟鴻鷹;劉貴忠;;地震波形反演的小波多尺度分解方法[A];1996年中國地球物理學(xué)會第十二屆學(xué)術(shù)年會論文集[C];1996年

5 盧回憶;劉伊克;常旭;白蘭淑;;彈性介質(zhì)多尺度全波形反演方法研究[A];中國地球物理2013——第二十二專題論文集[C];2013年

6 丁繼才;常旭;劉伊克;汪長永;;反射地震數(shù)據(jù)的逐層波形反演[A];中國科學(xué)院地質(zhì)與地球物理研究所2007學(xué)術(shù)論文匯編(第四卷)[C];2008年

7 劉璐;劉洪;張衡;崔永福;李飛;段文勝;彭更新;;基于修正擬牛頓公式的全波形反演[A];中國科學(xué)院地質(zhì)與地球物理研究所2013年度（第13屆）學(xué)術(shù)論文匯編——油氣資源研究室[C];2014年

8 吳建平;明躍紅;;區(qū)域地震波形反演及其在川滇地區(qū)的應(yīng)用[A];中國地球物理學(xué)會年刊2002——中國地球物理學(xué)會第十八屆年會論文集[C];2002年

9 何樵登;周輝;;人工神經(jīng)網(wǎng)絡(luò)在地震波形反演中的應(yīng)用[A];1995年中國地球物理學(xué)會第十一屆學(xué)術(shù)年會論文集[C];1995年

10 李志曄;李振春;張凱;;基于共散射點道集的全波形反演[A];中國地球物理2013——第二十二專題論文集[C];2013年

相關(guān)博士學(xué)位論文 前9條

1 陳宇澍;地震波場正演與集合震源編碼全波形反演算法[D];清華大學(xué);2015年

2 何清龍;基于有限差分—對比源方法的波動方程全波形反演研究[D];哈爾濱工業(yè)大學(xué);2016年

3 岳蕾;電磁波在煤層中的傳播規(guī)律與全波形概率反演方法研究[D];中國礦業(yè)大學(xué);2016年

4 張海如;探地雷達信號后處理關(guān)鍵技術(shù)研究[D];西安電子科技大學(xué);2016年

5 楊賀龍;基于波場優(yōu)化匹配的高精度全波形反演方法研究[D];吉林大學(xué);2016年

6 韓淼;深水區(qū)地震全波形反演策略與應(yīng)用[D];吉林大學(xué);2014年

7 秦寧;地震走時層析與波形反演方法研究[D];中國石油大學(xué)（華東）;2013年

8 郭振波;彈性介質(zhì)波形反演方法研究[D];中國石油大學(xué)(華東);2014年

9 吳俊軍;跨孔雷達全波形層析成像反演方法的研究[D];吉林大學(xué);2012年

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

1 張偉;GPU/CPU協(xié)同并行計算全波形反演理論及應(yīng)用[D];電子科技大學(xué);2014年

2 陳永芮;時間域全波形反演方法研究[D];中國石油大學(xué)(華東);2014年

3 王慶;二維地震時間域多尺度全波形反演[D];中國海洋大學(xué);2015年

4 孫瑞雪;基于雙參數(shù)正則化的頻率域解纏相位全波形反演研究[D];吉林大學(xué);2016年

5 魯曉滿;基于逐減隨機震源采樣法及改進MLQN算法的頻率域全波形反演方法研究[D];吉林大學(xué);2016年

6 孫慧秋;全波形反演地層品質(zhì)因子方法研究[D];吉林大學(xué);2016年

7 付振武;地震勘探全波形反演的近似Newton型方法[D];哈爾濱工業(yè)大學(xué);2016年

8 周鐵軍;近地表速度模型地震全波形反演技術(shù)研究[D];東北石油大學(xué);2016年

9 李志曄;基于地表參數(shù)擴展模型的波形反演方法研究[D];中國石油大學(xué)(華東);2015年

10 尹力;低頻缺失下全波形反演應(yīng)用研究[D];中國石油大學(xué)(華東);2015年

，

本文編號：1662768