亚洲免费av电影一区二区三区,日韩爱爱视频,51精品视频一区二区三区,91视频爱爱,日韩欧美在线播放视频,中文字幕少妇AV,亚洲电影中文字幕,久久久久亚洲av成人网址,久久综合视频网站,国产在线不卡免费播放

        ?

        基于小波變換提取吸收系數(shù)方法及應(yīng)用

        2013-06-17 05:52:38周喆桂志先王寧
        斷塊油氣田 2013年3期
        關(guān)鍵詞:比法斷塊時(shí)頻

        周喆,桂志先,王寧

        (1.長(zhǎng)江大學(xué)地球科學(xué)學(xué)院,湖北 荊州434023;2.油氣資源與勘探技術(shù)教育部重點(diǎn)實(shí)驗(yàn)室(長(zhǎng)江大學(xué)),湖北 荊州434023;3.長(zhǎng)江大學(xué)學(xué)校辦公室,湖北 荊州434023)

        地層是非完全彈性介質(zhì),當(dāng)?shù)卣鸩ㄔ诘貙咏橘|(zhì)中傳播時(shí),地震信號(hào)的能量被吸收,而衰減量包含了該地層的巖性及含流體信息,與其他參數(shù)相比,其具有更高的靈敏度[1-14];所以,研究地層的吸收衰減規(guī)律意義重大。根據(jù)頻譜比斜率確定的吸收系數(shù)比直接利用相鄰反射波振幅衰減比值計(jì)算的結(jié)果更加準(zhǔn)確,這是因?yàn)轭l率和吸收系數(shù)關(guān)系密切[15-19],頻譜能較好地反映出不同頻率的變化強(qiáng)度。

        1 Morlet 小波時(shí)頻分析

        小波變換時(shí)頻分析過(guò)程中,小波基的選擇很重要,在地球物理信號(hào)處理中常用的是復(fù)值Morlet 小波,該次研究選擇了Morlet 小波對(duì)地震資料進(jìn)行時(shí)頻分析。通過(guò)小波變換對(duì)地震數(shù)據(jù)進(jìn)行時(shí)頻分析,能更靈敏地反映地下含油氣特征,為直接進(jìn)行油氣檢測(cè)和油氣分析提供有利依據(jù)。

        Morlet 小波變換原理:

        其傅里葉變換為

        式中:h(t)為小波函數(shù);ω0為小波中心頻率,Hz;t 為時(shí)間,s;h(ω)為小波傅里葉變換頻率域函數(shù);ω 為小波頻率,Hz。

        由式(2)可看出,當(dāng)ω=0 Hz 時(shí),h(0) =0,滿(mǎn)足容許條件,因此Morlet 小波是容許小波。在信號(hào)處理應(yīng)用中取ω0≥5 Hz,可忽略式(1)中的第2 項(xiàng),而此時(shí)依然能滿(mǎn)足容許條件。Morlet 小波是復(fù)值小波,能夠提取信號(hào)中的幅值和相位信息,在地球物理信號(hào)處理中得到廣泛應(yīng)用,Trace 5540 地震剖面數(shù)據(jù)及利用Morlet 小波變換得到的頻譜見(jiàn)圖1。

        圖1 Trace 5540 線(xiàn)剖面

        由圖1可看出,地震信號(hào)主頻約為35 Hz,頻帶寬度20~50 Hz,頻帶較窄,與該區(qū)域提供的分析結(jié)果相符合。小波變換可以提取單一頻率的地震頻譜能量,圖2為T(mén)race 5540 線(xiàn)地震剖面65,35 Hz 的能量譜。

        圖2 Trace 5540 線(xiàn)能量剖面

        由圖2可看出,淺、中層的高頻能量要遠(yuǎn)遠(yuǎn)大于深層的高頻能量,這是地層對(duì)地震信號(hào)高頻成分的吸收衰減引起的。

        2 吸收系數(shù)提取

        利用Morlet 小波變換頻譜比法提取吸收系數(shù)。首先,對(duì)時(shí)間域振幅值進(jìn)行Morlet 小波變換時(shí)頻分析,得到頻率域的振幅值; 然后,利用時(shí)窗滑動(dòng)選擇振幅譜值; 最后,通過(guò)頻譜比法編寫(xiě)程序來(lái)確定頻譜比的斜率,即求出不同層間的吸收衰減系數(shù)。

        2.1 頻譜比法原理

        假設(shè)介質(zhì)吸收系數(shù)為α,則一維波動(dòng)方程的解為

        其中 α=α0f

        式中:c 為相速度,在地震勘探頻率帶內(nèi)可近似認(rèn)為是常數(shù);α 為吸收系數(shù),由介質(zhì)的品質(zhì)因素Q 及相速度c決定;f 為頻率,Hz;x 為波傳播距離,m。

        2.1.1 利用直達(dá)波計(jì)算吸收系數(shù)

        設(shè)F1(ω)為震源傳播距離x1處的脈沖復(fù)變譜,F(xiàn)2(ω)為震源傳播距離x2處的脈沖復(fù)變譜,由式(3)和有關(guān)頻譜定理可以知道,在均勻介質(zhì)中,F(xiàn)1(ω),F(xiàn)2(ω)的關(guān)系為

        如果介質(zhì)是不均勻的,則α 和c 可看成一段距離的平均值。

        2.1.2 利用反射波提取吸收系數(shù)

        設(shè)吸收介質(zhì)模型為水平層狀,各層品質(zhì)因素Q 與頻率無(wú)關(guān),即屬恒Q 介質(zhì);地震波自激自收,第n 層頂、底面反射波分別為Fn-1(ω),F(xiàn)n(ω)。忽略該層內(nèi)吸收系數(shù)對(duì)脈沖波形的改造及吸收作用對(duì)反射系數(shù)的影響,則頂、底反射波振幅譜之比可表示為

        式中:Dn為波前發(fā)射因子;Rn為第n 層反射系數(shù);Zn為第n 層界面深度,m;αn為第n 層吸收系數(shù);n 為層數(shù)。

        式(5)兩邊取對(duì)數(shù)得

        根據(jù)式(6),并考慮到前述假設(shè),則式(6)是頻率f的線(xiàn)性方程,其斜率就是有效吸收系數(shù)α0,式(6)可變換為

        由式(7)可看出,波前發(fā)射因子只影響直線(xiàn)方程截距,并不影響斜率,所以式(7)可用于求取有效吸收系數(shù)α0。

        2.2 計(jì)算方法

        頻譜比法求取吸收系數(shù)按式(7)進(jìn)行,具體步驟:1)拾取反射波;2)利用Morlet 小波變換求取各層的反射波頻譜,進(jìn)而求得相位譜與振幅譜;3)計(jì)算相鄰界面之間振幅譜的比值,同時(shí)求取其自然對(duì)數(shù);4)按式(7)擬合直線(xiàn),所求得的斜率值就是有效吸收系數(shù)。

        另外,在對(duì)直線(xiàn)擬合之前,可以將振幅譜比的自然對(duì)數(shù)做橫向平滑。這是因?yàn)榇嬖诟鞣N干擾,所得振幅譜比值自然對(duì)數(shù)與頻率f 的關(guān)系不可能為真的直線(xiàn)[16]。因此,根據(jù)單個(gè)記錄道求α0誤差較大,故采用橫向平滑等手段以消除干擾影響。

        3 應(yīng)用實(shí)例

        研究區(qū)含氣層主要分布在Es33段,為深—半深湖相的油頁(yè)巖、砂泥與鹽膏巖交替沉積不等厚互層,經(jīng)勘探,Es33段含氣較豐富。

        WL15 井是研究區(qū)的一口評(píng)價(jià)井,在Es33層鉆遇并發(fā)現(xiàn)工業(yè)氣流。WL7 和WL9 井也發(fā)現(xiàn)工業(yè)氣流,從而推斷該區(qū)域富含油氣。為更好了解該區(qū)域油氣分布與儲(chǔ)集層變化的特征,利用吸收系數(shù)算法,對(duì)該區(qū)地震數(shù)據(jù)進(jìn)行處理(見(jiàn)圖3、圖4),分析該區(qū)域的吸收特征。

        圖3 Trace 4890 地震吸收系數(shù)剖面

        圖4 Line 4550 地震吸收系數(shù)剖面

        圖3與圖4為過(guò)WL15 井的2 條交叉測(cè)線(xiàn)。由圖3可看出,WL15 井在目的層Es33層,其吸收特性比附近地層要強(qiáng)得多。鉆探后發(fā)現(xiàn)該井段所處位置是含氣層,獲得工業(yè)氣流,與實(shí)際鉆探結(jié)果吻合較好。由圖4也可看出同樣的吸收特性。實(shí)際結(jié)果表明,該方法能夠較好識(shí)別含氣層。

        4 結(jié)束語(yǔ)

        利用基于Morlet 小波變換的頻譜比法,對(duì)研究區(qū)地震資料進(jìn)行處理。結(jié)果表明,該算法對(duì)氣層識(shí)別的效果較好,算法比較簡(jiǎn)單,比單獨(dú)利用時(shí)間域振幅的識(shí)別效果好。

        [1]劉企英.利用地震信息進(jìn)行油氣預(yù)測(cè)[M].北京:石油工業(yè)出版社,1995:136-145.

        [2]邊立恩,于茜,韓自軍,等.基于譜分解技術(shù)的儲(chǔ)層定量地震解釋?zhuān)跩].石油與天然氣地質(zhì),2011,32(5):718-723.

        [3]朱廣生.地震資料儲(chǔ)層預(yù)測(cè)方法[M].北京:石油工業(yè)出版社,1995:56-73.

        [4]李慶忠.走向精確勘探的道路[M].北京:石油工業(yè)出版社,1993:78-102.

        [5]俞壽朋.高分辨率地震勘探[M].北京:石油工業(yè)出版社,1993:32-45.

        [6]董潔,陳世悅,袁波.時(shí)頻分析精細(xì)劃分歡喜嶺地區(qū)地層層序[J].斷塊油氣田,2010,17(2):146-149.

        [7]袁井菊,賈慶素,施少武.大情字井地區(qū)特征參數(shù)反演儲(chǔ)層預(yù)測(cè)[J].斷塊油氣田,2010,17(4):426-430.

        [8]郭偉,代宗仰,尹宏,等.疊后波阻抗約束反演的若干問(wèn)題[J].斷塊油氣田,2006,13(6):16-17.

        [9]Hough S E,Lees J M,Monastero F.Attenuation and source properties at the Coso Geothermal Area ,California [J ].Bulletin of the Seismological Society of America,1999,89(6):1606-1619.

        [10]Vladislav G Martynov,F(xiàn)rank L Vernon,Robert J Mellors,et al.Highfrequency attenuation in the crust and Upper Mantle of the Northern Tien Shan [J].Bulletin of the Seismological Society of America,1999,89(1):215-238.

        [11]Rahul Dasgupta,Roger A Clark.Estimation of Q from surface seismic reflection data[J].Geophysics,1995,63(6):2120-2128.

        [12]王西文,劉全新,趙應(yīng)成,等.地震信號(hào)瞬時(shí)特征在小波域分頻提取的方法和應(yīng)用[J].石油地球物理勘探,2000,35(4):452-460,478.

        [13]王西文.小波域地震資料處理方法(WSSP)研究[D].北京:中國(guó)科學(xué)院地質(zhì)與地球物理研究所,2000.

        [14]高靜懷,汪文秉,朱光明.小波變換與信號(hào)瞬時(shí)特征分析[J].地球物理學(xué)報(bào),1997,40(6):821-832.

        [15]陳培善.地震波研究在我國(guó)的進(jìn)展[J].地球物理學(xué)報(bào),1994,37(1):231-241.

        [16]辛可峰,李振春,王永剛,等.地層等效吸收系數(shù)反演[J].石油物探,2001,40(4):14-20.

        [17]夏雨,陳世悅,楊俊生.準(zhǔn)噶爾盆地烏夏地區(qū)二疊系層序地層特征研究[J].?dāng)鄩K油氣田,2008,15(1):18-21.

        [18]董臣強(qiáng),王軍,張金偉.時(shí)頻分析技術(shù)在三角洲層序分析中的應(yīng)用[J].斷塊油氣田,2002,9(2):18-20.

        [19]曹健,趙憲生,施澤進(jìn),等.振幅、頻率動(dòng)態(tài)吸收系數(shù)在嘉一段儲(chǔ)層預(yù)測(cè)中的應(yīng)用[J].特種油氣藏,2010,17(6):35-37.

        猜你喜歡
        比法斷塊時(shí)頻
        化虛為實(shí) 觸摸物理——物理方法之類(lèi)比法
        加權(quán)譜比法Q值估計(jì)
        物理方法之類(lèi)比法
        復(fù)雜斷塊油藏三維地質(zhì)模型的多級(jí)定量評(píng)價(jià)
        最好的比較
        斷塊油藏注采耦合物理模擬實(shí)驗(yàn)
        港中油田南一斷塊高含水后期提高開(kāi)發(fā)效果研究
        基于時(shí)頻分析的逆合成孔徑雷達(dá)成像技術(shù)
        對(duì)采樣數(shù)據(jù)序列進(jìn)行時(shí)頻分解法的改進(jìn)
        雙線(xiàn)性時(shí)頻分布交叉項(xiàng)提取及損傷識(shí)別應(yīng)用
        国产精品51麻豆cm传媒| 日韩一区二区中文字幕视频| 在线亚洲妇色中文色综合| 国产影片一区二区三区| 亚洲av成人片色在线观看| 久久亚洲av无码西西人体 | 亚洲色欲久久久综合网| 伊伊人成亚洲综合人网7777| 日日躁欧美老妇| 国产精品人成在线观看不卡| 亚洲中文字幕在线一区| 精品久久久久久久无码人妻热| 人妻在线日韩免费视频| 成人无码区免费AⅤ片WWW| 青春草在线观看免费视频| 亚洲乱码中文字幕一线区| 天天做天天爱夜夜爽女人爽 | 亚洲av无码精品色午夜| 高清国产一级毛片国语| 天堂av中文在线官网| 亚洲精品视频1区2区| 丰满少妇a级毛片| www国产亚洲精品久久网站| 91精品国产91热久久p| 亚洲国产精品成人一区二区三区| 亚洲综合偷自成人网第页色| 内射人妻少妇无码一本一道| 少妇人妻偷人精品视蜜桃| 久久久久久免费播放一级毛片| 成人亚洲av网站在线看| 日韩人妻无码精品一专区二区三区 | 91情侣视频| 激情乱码一区二区三区| 精品国内日本一区二区| 疯狂的欧美乱大交| 亚洲暴爽av天天爽日日碰| 中文字幕一区二区人妻痴汉电车| 中文字幕人妻精品一区| 无码中文字幕人妻在线一区| 欧美两根一起进3p做受视频| 国产国拍亚洲精品福利|