稀疏脈沖反演方法及Jason軟件實(shí)現(xiàn)

1、第四篇約束稀疏脈沖波阻抗反演1、概述：約束稀疏脈沖波阻抗反演方法(簡稱CSSI)是目前應(yīng)用比擬廣泛的一種方法,它是以趨勢約束稀疏脈沖反褶積算法為根底.對地下的反射吸收序列作如下假設(shè)：它的分布是稀疏分布的,而不是像其它的反演方法那樣假設(shè)其是連續(xù)分布的.在反演過程中,它是通過尋找一個(gè)使目標(biāo)函數(shù)最小的反射系數(shù)脈沖數(shù)目來得到波阻抗數(shù)據(jù)的,此時(shí)合成的地震記錄可以和實(shí)際地震記錄匹配的最正確.2、根本假設(shè)和方法原理約束稀疏脈沖反演(CSSI)是一種遞推地震波阻抗反演方法,基于范數(shù)反褶積和最大似然反褶積.它假設(shè)地層的反射系數(shù)序列是由一系列稀疏且不連續(xù)的大反射系數(shù)和高斯分布的小反射系數(shù)疊加而成.從地質(zhì)角都解釋,

2、大反射系數(shù)代表地下不連續(xù)界面、巖性分界面.首先由約束稀疏脈沖的根本假設(shè),得到目標(biāo)函數(shù)：1Lo1LJ=至£(K)+m£n2(K)-2Mln(?.)-2L-M)ln(1-Z)(1)RKwNKR2為反射系數(shù)均方值,N2為噪音均方值；第K個(gè)采樣點(diǎn)的反射系數(shù)為r(K),n(K)為噪音；M為反射層數(shù)；L為總采樣數(shù)；K為既定反射系數(shù)似然值；可經(jīng)迭代計(jì)算出反射系數(shù).其次用最大似然反褶積反演得到反射系數(shù).采用遞推算法結(jié)合初始波阻抗趨勢,得到波阻抗模型.Z(i)=Z(iZ(i)為第i層的波阻抗值；R(i)為第i層的反射系數(shù).最后約束稀疏脈沖反演依據(jù)目標(biāo)函數(shù)逐一對地震道計(jì)算的初始波阻抗進(jìn)行調(diào)整,

3、并調(diào)整反射系數(shù),得到目標(biāo)函數(shù)：F=Lp(r)+人(s-d)+aZ(3)r是反射系數(shù),AZ是阻抗趨勢差；d是地震道；s是合成地震道；兒是殘差權(quán)重因子；a是趨勢權(quán)重因子；p,q是L模因子.式(3)右邊分別是反射系數(shù)的絕對值和,合成聲波記錄與原始地震數(shù)據(jù)的差值,趨勢約束項(xiàng).約束稀疏脈沖反演基于地震數(shù)據(jù),用波阻抗趨勢進(jìn)行約束,以期用最少的反射系數(shù)脈沖,對合成記錄與地震道進(jìn)行最正確匹配.3、CSSI反演流程先對地震資料細(xì)致處理并提取子波；對地震層位精細(xì)解釋,建立合理地質(zhì)地層模型；關(guān)鍵是對測井資料做標(biāo)準(zhǔn)化分析,結(jié)合地震子波做合成記錄,以保證時(shí)間深度域準(zhǔn)確對應(yīng),結(jié)合地層框架模型和合成記錄,以三維疊后地震數(shù)據(jù)

4、為根底,利用地震子波以及初始模型趨勢模型,進(jìn)行波阻抗反演,最終得到精細(xì)反演結(jié)果.反演流程如下：塞波浦橙科¥液拄取限制學(xué)數(shù)mi道合并T4、二維模型測試選用二維mamous模型,利用Jason軟件進(jìn)行稀疏脈沖反演,目標(biāo)油氣層大約在0.35s處.其根本流程如下:圖4-1:原始地震道圖4-2:解釋層位rAc-<_圖4-3:井震標(biāo)定：相關(guān)系數(shù)為0.82圖4-5:反演結(jié)果5、實(shí)際資料反演Jason軟件來選用某工區(qū)的實(shí)際疊后資料來進(jìn)行疊后的稀疏約束脈沖反演,本文采用的是進(jìn)行反演.其根本流程如下：(1)確定反演的范圍,解釋層位：5-1:原始地震剖面及解釋的層位5-2:原始地震的連井剖面5-2:原始地震剖面在頻率域的顯示(2)對反演范圍內(nèi)的井進(jìn)行井震標(biāo)定:5-2:對工區(qū)CB32井標(biāo)定,相關(guān)系數(shù)0.41,=,£5-3:對工區(qū)CB804井標(biāo)定,相關(guān)系數(shù)0.76(3)構(gòu)建低頻模型5-4:固體模型5-5:低頻模型(4)