基于疊前地震數(shù)據(jù)振幅隨入射角變化特征的疊前反演方法是獲得儲層彈性參數(shù)估計(jì)的重要手段。常規(guī)貝葉斯疊前反演假設(shè)彈性參數(shù)反射系數(shù)滿足某一特定分布,利用該分布作為先驗(yàn)信息約束反演,且分布一經(jīng)確定不再隨不同地震道變化,但在彈性參數(shù)反射系數(shù)的分布橫向變化較大或不滿足既定分布時會產(chǎn)生誤差。為此,通過引入自動相關(guān)判別先驗(yàn),提出了一種貝葉斯疊前同時反演方法,將彈性參數(shù)反射系數(shù)的先驗(yàn)信息融合到反演中,該先驗(yàn)信息不再對彈性參數(shù)反射系數(shù)的總體分布作出假設(shè),且隨著地震數(shù)據(jù)的變化自適應(yīng)改變,使彈性參數(shù)反射系數(shù)的反演結(jié)果具有更好的準(zhǔn)確性、更符合橫向變化的地質(zhì)特征,并且更好地反映巖性邊界;在構(gòu)建反演目標(biāo)函數(shù)時,加入了趨勢約束算子,提高了反演的穩(wěn)定性;采用同時反演策略,避免了基于彈性阻抗的疊前反演過程中存在轉(zhuǎn)換誤差。單道彈性參數(shù)模型和經(jīng)典推覆體2D模型測試結(jié)果證實(shí)了方法的可行性和正確性,實(shí)際地震數(shù)據(jù)測試結(jié)果證實(shí)了方法在提高分辨率和反映地質(zhì)體特征及邊界上的優(yōu)勢。 

【文章來源】：石油物探. 2020,59(04)北大核心CSCD

【文章頁數(shù)】：11 頁

【部分圖文】：

無噪和含噪模擬數(shù)據(jù)對本文方法的測試結(jié)果

采用經(jīng)典的推覆體2D模型進(jìn)一步測試本文方法,利用30 Hz雷克子波并加入10%的隨機(jī)噪聲獲得了含噪聲的、采樣間隔為1 ms的不同角度模擬數(shù)據(jù)(入射角度分別為10°、20°和30°),模擬數(shù)據(jù)的縱橫波阻抗參數(shù)模型如圖2a和圖2b所示;利用ARD先驗(yàn),但未加入趨勢約束(參數(shù)模型的0～6 Hz成分)的反演結(jié)果如圖2c和圖2d所示;加入趨勢約束的反演結(jié)果如圖2e和圖2f所示�？傮w上看,反演獲得的縱波阻抗結(jié)果優(yōu)于橫波阻抗結(jié)果;趨勢約束的加入提高了反演結(jié)果的橫向連續(xù)性,并且減輕了隨機(jī)噪聲引起的“掛面條”現(xiàn)象,這一提高對于穩(wěn)定性相對較差的橫波阻抗反演結(jié)果更為明顯。2D模型的測試結(jié)果進(jìn)一步證實(shí)了本文方法在獲得高分辨率反演結(jié)果和保持反演結(jié)果橫向穩(wěn)定性方面的有效性和優(yōu)勢。3 實(shí)際數(shù)據(jù)應(yīng)用

選取實(shí)際數(shù)據(jù)中某二維剖面作為實(shí)驗(yàn)數(shù)據(jù)。不同入射角的疊前部分疊加地震數(shù)據(jù)如圖4所示,采用的入射角分別為7°、14°和21°。采用的低頻趨勢約束由提取測井彈性參數(shù)曲線插值模型的低頻成分(0～6 Hz)獲得。通過井震標(biāo)定獲得不同角度子波進(jìn)行反演。反演時選取10～70 Hz作為輸入數(shù)據(jù)的頻帶范圍,反演后將結(jié)果中高于150 Hz的成分濾除。圖5a和圖5b分別給出了采用本文方法反演得到的縱波阻抗和橫波阻抗。由圖5a和圖5b可見,本文方法得到的反演結(jié)果分辨率較高,地層界面清晰,與地質(zhì)認(rèn)識較為吻合,為后續(xù)儲層和地質(zhì)解釋提供了較好的基礎(chǔ)資料。圖5c和圖5d分別給出了采用某商業(yè)軟件得到的縱波阻抗和橫波阻抗反演結(jié)果。由圖5可見,本文方法具有更高的分辨率,尤其是橫波阻抗的反演結(jié)果優(yōu)勢明顯,具有更好的可解釋性,有利于后續(xù)儲層與地質(zhì)解釋。不足的是,盡管橫波阻抗反演結(jié)果的改善很大,但剖面局部仍存在一些由于噪聲和振幅畸變造成的不穩(wěn)定現(xiàn)象,后續(xù)研究可以重點(diǎn)嘗試?yán)枚嗟啦呗约訌?qiáng)反演結(jié)果的橫向穩(wěn)定性。本文方法的計(jì)算效率在可接受的范圍內(nèi),但略遜于商業(yè)軟件,后續(xù)研究可以進(jìn)行針對性的算法優(yōu)化以提高效率。圖4 具有不同入射角的疊前部分疊加地震數(shù)據(jù)


本文編號：3110341