基于自動(dòng)相關(guān)判別先驗(yàn)的疊前同時(shí)反演方法研究
發(fā)布時(shí)間:2021-03-30 23:31
基于疊前地震數(shù)據(jù)振幅隨入射角變化特征的疊前反演方法是獲得儲(chǔ)層彈性參數(shù)估計(jì)的重要手段。常規(guī)貝葉斯疊前反演假設(shè)彈性參數(shù)反射系數(shù)滿足某一特定分布,利用該分布作為先驗(yàn)信息約束反演,且分布一經(jīng)確定不再隨不同地震道變化,但在彈性參數(shù)反射系數(shù)的分布橫向變化較大或不滿足既定分布時(shí)會(huì)產(chǎn)生誤差。為此,通過(guò)引入自動(dòng)相關(guān)判別先驗(yàn),提出了一種貝葉斯疊前同時(shí)反演方法,將彈性參數(shù)反射系數(shù)的先驗(yàn)信息融合到反演中,該先驗(yàn)信息不再對(duì)彈性參數(shù)反射系數(shù)的總體分布作出假設(shè),且隨著地震數(shù)據(jù)的變化自適應(yīng)改變,使彈性參數(shù)反射系數(shù)的反演結(jié)果具有更好的準(zhǔn)確性、更符合橫向變化的地質(zhì)特征,并且更好地反映巖性邊界;在構(gòu)建反演目標(biāo)函數(shù)時(shí),加入了趨勢(shì)約束算子,提高了反演的穩(wěn)定性;采用同時(shí)反演策略,避免了基于彈性阻抗的疊前反演過(guò)程中存在轉(zhuǎn)換誤差。單道彈性參數(shù)模型和經(jīng)典推覆體2D模型測(cè)試結(jié)果證實(shí)了方法的可行性和正確性,實(shí)際地震數(shù)據(jù)測(cè)試結(jié)果證實(shí)了方法在提高分辨率和反映地質(zhì)體特征及邊界上的優(yōu)勢(shì)。
【文章來(lái)源】:石油物探. 2020,59(04)北大核心CSCD
【文章頁(yè)數(shù)】:11 頁(yè)
【部分圖文】:
無(wú)噪和含噪模擬數(shù)據(jù)對(duì)本文方法的測(cè)試結(jié)果
采用經(jīng)典的推覆體2D模型進(jìn)一步測(cè)試本文方法,利用30 Hz雷克子波并加入10%的隨機(jī)噪聲獲得了含噪聲的、采樣間隔為1 ms的不同角度模擬數(shù)據(jù)(入射角度分別為10°、20°和30°),模擬數(shù)據(jù)的縱橫波阻抗參數(shù)模型如圖2a和圖2b所示;利用ARD先驗(yàn),但未加入趨勢(shì)約束(參數(shù)模型的0~6 Hz成分)的反演結(jié)果如圖2c和圖2d所示;加入趨勢(shì)約束的反演結(jié)果如圖2e和圖2f所示。總體上看,反演獲得的縱波阻抗結(jié)果優(yōu)于橫波阻抗結(jié)果;趨勢(shì)約束的加入提高了反演結(jié)果的橫向連續(xù)性,并且減輕了隨機(jī)噪聲引起的“掛面條”現(xiàn)象,這一提高對(duì)于穩(wěn)定性相對(duì)較差的橫波阻抗反演結(jié)果更為明顯。2D模型的測(cè)試結(jié)果進(jìn)一步證實(shí)了本文方法在獲得高分辨率反演結(jié)果和保持反演結(jié)果橫向穩(wěn)定性方面的有效性和優(yōu)勢(shì)。3 實(shí)際數(shù)據(jù)應(yīng)用
選取實(shí)際數(shù)據(jù)中某二維剖面作為實(shí)驗(yàn)數(shù)據(jù)。不同入射角的疊前部分疊加地震數(shù)據(jù)如圖4所示,采用的入射角分別為7°、14°和21°。采用的低頻趨勢(shì)約束由提取測(cè)井彈性參數(shù)曲線插值模型的低頻成分(0~6 Hz)獲得。通過(guò)井震標(biāo)定獲得不同角度子波進(jìn)行反演。反演時(shí)選取10~70 Hz作為輸入數(shù)據(jù)的頻帶范圍,反演后將結(jié)果中高于150 Hz的成分濾除。圖5a和圖5b分別給出了采用本文方法反演得到的縱波阻抗和橫波阻抗。由圖5a和圖5b可見,本文方法得到的反演結(jié)果分辨率較高,地層界面清晰,與地質(zhì)認(rèn)識(shí)較為吻合,為后續(xù)儲(chǔ)層和地質(zhì)解釋提供了較好的基礎(chǔ)資料。圖5c和圖5d分別給出了采用某商業(yè)軟件得到的縱波阻抗和橫波阻抗反演結(jié)果。由圖5可見,本文方法具有更高的分辨率,尤其是橫波阻抗的反演結(jié)果優(yōu)勢(shì)明顯,具有更好的可解釋性,有利于后續(xù)儲(chǔ)層與地質(zhì)解釋。不足的是,盡管橫波阻抗反演結(jié)果的改善很大,但剖面局部仍存在一些由于噪聲和振幅畸變?cè)斐傻牟环(wěn)定現(xiàn)象,后續(xù)研究可以重點(diǎn)嘗試?yán)枚嗟啦呗约訌?qiáng)反演結(jié)果的橫向穩(wěn)定性。本文方法的計(jì)算效率在可接受的范圍內(nèi),但略遜于商業(yè)軟件,后續(xù)研究可以進(jìn)行針對(duì)性的算法優(yōu)化以提高效率。圖4 具有不同入射角的疊前部分疊加地震數(shù)據(jù)
本文編號(hào):3110341
【文章來(lái)源】:石油物探. 2020,59(04)北大核心CSCD
【文章頁(yè)數(shù)】:11 頁(yè)
【部分圖文】:
無(wú)噪和含噪模擬數(shù)據(jù)對(duì)本文方法的測(cè)試結(jié)果
采用經(jīng)典的推覆體2D模型進(jìn)一步測(cè)試本文方法,利用30 Hz雷克子波并加入10%的隨機(jī)噪聲獲得了含噪聲的、采樣間隔為1 ms的不同角度模擬數(shù)據(jù)(入射角度分別為10°、20°和30°),模擬數(shù)據(jù)的縱橫波阻抗參數(shù)模型如圖2a和圖2b所示;利用ARD先驗(yàn),但未加入趨勢(shì)約束(參數(shù)模型的0~6 Hz成分)的反演結(jié)果如圖2c和圖2d所示;加入趨勢(shì)約束的反演結(jié)果如圖2e和圖2f所示。總體上看,反演獲得的縱波阻抗結(jié)果優(yōu)于橫波阻抗結(jié)果;趨勢(shì)約束的加入提高了反演結(jié)果的橫向連續(xù)性,并且減輕了隨機(jī)噪聲引起的“掛面條”現(xiàn)象,這一提高對(duì)于穩(wěn)定性相對(duì)較差的橫波阻抗反演結(jié)果更為明顯。2D模型的測(cè)試結(jié)果進(jìn)一步證實(shí)了本文方法在獲得高分辨率反演結(jié)果和保持反演結(jié)果橫向穩(wěn)定性方面的有效性和優(yōu)勢(shì)。3 實(shí)際數(shù)據(jù)應(yīng)用
選取實(shí)際數(shù)據(jù)中某二維剖面作為實(shí)驗(yàn)數(shù)據(jù)。不同入射角的疊前部分疊加地震數(shù)據(jù)如圖4所示,采用的入射角分別為7°、14°和21°。采用的低頻趨勢(shì)約束由提取測(cè)井彈性參數(shù)曲線插值模型的低頻成分(0~6 Hz)獲得。通過(guò)井震標(biāo)定獲得不同角度子波進(jìn)行反演。反演時(shí)選取10~70 Hz作為輸入數(shù)據(jù)的頻帶范圍,反演后將結(jié)果中高于150 Hz的成分濾除。圖5a和圖5b分別給出了采用本文方法反演得到的縱波阻抗和橫波阻抗。由圖5a和圖5b可見,本文方法得到的反演結(jié)果分辨率較高,地層界面清晰,與地質(zhì)認(rèn)識(shí)較為吻合,為后續(xù)儲(chǔ)層和地質(zhì)解釋提供了較好的基礎(chǔ)資料。圖5c和圖5d分別給出了采用某商業(yè)軟件得到的縱波阻抗和橫波阻抗反演結(jié)果。由圖5可見,本文方法具有更高的分辨率,尤其是橫波阻抗的反演結(jié)果優(yōu)勢(shì)明顯,具有更好的可解釋性,有利于后續(xù)儲(chǔ)層與地質(zhì)解釋。不足的是,盡管橫波阻抗反演結(jié)果的改善很大,但剖面局部仍存在一些由于噪聲和振幅畸變?cè)斐傻牟环(wěn)定現(xiàn)象,后續(xù)研究可以重點(diǎn)嘗試?yán)枚嗟啦呗约訌?qiáng)反演結(jié)果的橫向穩(wěn)定性。本文方法的計(jì)算效率在可接受的范圍內(nèi),但略遜于商業(yè)軟件,后續(xù)研究可以進(jìn)行針對(duì)性的算法優(yōu)化以提高效率。圖4 具有不同入射角的疊前部分疊加地震數(shù)據(jù)
本文編號(hào):3110341
本文鏈接:http://sikaile.net/projectlw/dqwllw/3110341.html
最近更新
教材專著