• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    基于VSP地震資料的品質(zhì)因子提取方法

    2016-10-18 06:23:27紀(jì)春玲周安聘張環(huán)曦趙雨晨牛淑瑜
    華北地震科學(xué) 2016年3期
    關(guān)鍵詞:時(shí)窗比法計(jì)算結(jié)果

    董 博,紀(jì)春玲,周安聘,李 鳳,張環(huán)曦,趙雨晨,劉 靜,牛淑瑜

    (1.河北地質(zhì)大學(xué),石家莊 050031;2.河北省地震局石家莊中心臺(tái),石家莊 050021)

    ?

    基于VSP地震資料的品質(zhì)因子提取方法

    董博1,2,紀(jì)春玲1,周安聘1,李鳳1,張環(huán)曦1,趙雨晨1,劉靜1,牛淑瑜1

    (1.河北地質(zhì)大學(xué),石家莊050031;2.河北省地震局石家莊中心臺(tái),石家莊050021)

    介紹了改進(jìn)的頻譜比方法原理、影響因素及其優(yōu)點(diǎn),并利用改進(jìn)的譜比法對(duì)勝利油田某工區(qū)采集到的VSP地震資料進(jìn)行測(cè)試,對(duì)所有下行反射波進(jìn)行振幅譜計(jì)算。結(jié)果顯示,該方法使得提取出來的品質(zhì)因子Q值計(jì)算精度得到了提高,計(jì)算結(jié)果更加真實(shí)可靠。

    品質(zhì)因子;改進(jìn)的譜比法;吸收衰減;傅里葉變換;影響因素

    0 引言

    品質(zhì)因子Q能夠量化為介質(zhì)的粘彈性所導(dǎo)致的地震波傳播過程中能量的衰減和頻散,它是用來表征巖層吸收衰減特征的參數(shù)[1]。彈性介質(zhì)中傳播的地震波沒有能量的衰減,因?yàn)榻橘|(zhì)的Q值趨于無窮大。介質(zhì)的粘彈性和Q值呈反比例關(guān)系,隨著Q值的增大,介質(zhì)的粘彈性減小,介質(zhì)的吸收衰減作用也就越小,反之亦成立[2]。地下介質(zhì)對(duì)地震波的吸收衰減作用包括巖層的非固有衰減和固有衰減[2]。在計(jì)算地層的固有衰減時(shí),應(yīng)該首先消除幾何擴(kuò)散、非彈性以外的散射、與頻率有關(guān)的透射和反射等非固有衰減因素的影響。因?yàn)闊o法完全地分離并消除這些因素,在計(jì)算等效Q值時(shí)還是會(huì)受到這些因素的影響[3]。

    1 常用品質(zhì)因子的計(jì)算方法

    根據(jù)計(jì)算方法原理的不同,可以把傳統(tǒng)品質(zhì)因子Q的計(jì)算方法分為時(shí)間域和頻率域方法。因?yàn)槊糠N計(jì)算方法的適用條件和對(duì)地震資料信噪比的要求不一樣[2],所有方法均不具有普遍適用性。

    時(shí)間域中計(jì)算品質(zhì)因子Q的方法主要包括:相位模擬法、上升時(shí)間法、脈沖振幅法、解析信號(hào)法、子波模擬法、瞬時(shí)頻率模擬法和振幅衰減法等。時(shí)間域求取Q值的方法所面臨的困難是地震資料的信噪比較低,振幅信息不保真,還不能通過有效的技術(shù)方法進(jìn)行準(zhǔn)確校正,因而時(shí)間域方法計(jì)算精度不高,在實(shí)際資料的處理中應(yīng)用很少。

    頻率域計(jì)算品質(zhì)因子Q的方法主要有:頻譜比法、譜模擬法、質(zhì)心頻率偏移法、匹配擬合技術(shù)等。但是每種方法的適用條件和計(jì)算精度都不太一樣。例如譜模擬方法[4-5],是Rosa等人提出在地震子波振幅譜的基礎(chǔ)上進(jìn)行計(jì)算的[6],主要是考慮了振幅的變化,而沒有考慮相位的變化。頻譜比法同時(shí)考慮了相位變化和振幅變化。

    頻譜比法在計(jì)算中涉及到的振幅譜必須是在目的層的上層位和下層位中分別選取得到的,然后利用選取的地震子波記錄段進(jìn)行頻譜分析和Q值的計(jì)算。假設(shè)地震波在t1和t2處的振幅譜分別為[7]:

    (1)

    (2)

    式中:A(t)是與頻率無關(guān)的函數(shù),代表其它因素對(duì)地震波衰減的影響。B(f)是初始時(shí)刻地震子波的振幅譜。觀察可以發(fā)現(xiàn),譜比法中頻譜比的斜率是Q的函數(shù),對(duì)頻譜比等式兩邊取對(duì)數(shù)得:

    (3)

    (4)

    由式(4)可以計(jì)算出地層等效Q值[8]。目的層厚度隨時(shí)間間隔的增大而變大,譜比法的計(jì)算結(jié)果也逐漸趨近于一個(gè)常數(shù)。如果排除了非固有衰減對(duì)Q值估計(jì)的影響,那么等效Q值只與地層介質(zhì)的吸收性質(zhì)有關(guān)。

    2 改進(jìn)的頻譜比法基本原理

    零偏移距VSP資料可以為估計(jì)品質(zhì)因子值提供很好的數(shù)據(jù),我們利用近偏移距VSP資料估計(jì)品質(zhì)因子值。改進(jìn)的譜比法是利用所有的下行反射波而不是只利用直達(dá)波進(jìn)行計(jì)算[9]。首先,我們利用波場(chǎng)分離將下行波場(chǎng)從近偏移距VSP資料中分離出來;然后,在頻率域中利用下行波場(chǎng)進(jìn)行計(jì)算。將地震數(shù)據(jù)從時(shí)間域變換到頻率域常用的方法是S變換和傅立葉變換。

    為了提高Q值的提取精度,需要在頻率域進(jìn)行計(jì)算,將地震數(shù)據(jù)變換到頻率域常用的有傅立葉變換和S變換方法,尤其是傅立葉變換。利用傅立葉變換將地震記錄變換到頻率域,在頻率域提取Q值[9],近偏移距VSP資料Q值估計(jì)通常是假設(shè)在頻率域中吸收是線性相關(guān)的(常數(shù)Q),譜比法的表達(dá)形式是:

    S2(f)=C×S1(f)e-α1,2t1,2f。

    (5)

    其中,吸收系數(shù)與Q的關(guān)系為:

    Q=π/α1,2,

    (6)

    對(duì)方程(5)取對(duì)數(shù):

    log(S2(f)/S1(f))=b-α1,2t1,2f 。

    (7)

    其中:b=log(C),為了估計(jì)吸收系數(shù)α,并且避免產(chǎn)生較大的誤差,最小化函數(shù)G:

    (8)

    式中:

    (9)

    將(8)帶入(9)可以得到目標(biāo)函數(shù):

    (10)

    在給定的α范圍內(nèi),找出G1(α)的最小值,從而可以得到α值。由于VSP地震數(shù)據(jù)信噪比較高,利用VSP地震數(shù)據(jù)進(jìn)行計(jì)算時(shí)沒有必要再進(jìn)行道集歸一化或者是幾何擴(kuò)散補(bǔ)償。從公式(5)可以推出zm,zk層振幅之間的關(guān)系表達(dá)式:

    Sm(f)=Const×Sk(f)e-Ak,m(Tm-Tk)f。

    (11)

    (12)

    從方程(11)和(12)可以得到平均速度和層間速度關(guān)系以及有效吸收系數(shù)與層間吸收系數(shù)之間的關(guān)系。不同的是,我們知道從地表到最淺層檢波器的平均速度,而不知道從地表到最淺層檢波器的吸收系數(shù),有效吸收系數(shù)和平均速度的表達(dá)式如下:

    (13)

    在方程式(13)中,只保存那些在給定范圍之內(nèi)的Ak,m。如果所有的Qk,m都在合理的范圍內(nèi),則應(yīng)該有N(N-1)/2個(gè)值。通常情況下,大概有總數(shù)的1/4,這為穩(wěn)定的Q值估計(jì)提供了良好的統(tǒng)計(jì)條件。

    3 實(shí)際資料處理

    為保證結(jié)果的可靠性,本文3.1和3.2節(jié)所選用數(shù)據(jù)均來自于勝利油田某工區(qū)周邊的VSP炮記錄,與3.3節(jié)所用VSP炮記錄場(chǎng)地相鄰,數(shù)據(jù)共100道,采樣間隔為1 ms,資料處理模型選用層狀衰減—頻率曲線的粘彈性VSP地震響應(yīng)模型。

    3.1改進(jìn)的譜比方法的影響因素

    地下介質(zhì)中地震波能量的衰減受到很多因素的影響,其中有些因素與地震波的頻率無關(guān),例如反射透射損失、波前發(fā)散等。除了地層吸收作用的影響以外,其它因素在地震波的頻率發(fā)生變化時(shí)對(duì)地震波的能量產(chǎn)生的影響很小。從地震波衰減機(jī)理上講,采用譜分析方法估算地層的吸收衰減性質(zhì)是可行的[7]。

    譜比法是計(jì)算品質(zhì)因子Q值的最常用方法,而改進(jìn)的譜比法雖然VSP地震記錄的信噪比較高,可以計(jì)算出較精確的品質(zhì)因子Q值,但是仍然受到如時(shí)窗的長(zhǎng)度、形狀等諸多因素的影響。

    層狀衰減—頻率曲線的粘彈性VSP地震響應(yīng)模型分3層(圖1),第1層Q值為200,第2層的Q值為30,第3層的Q值為40。上行反射波能量相對(duì)于下行反射波能量小的多,因?yàn)橄滦胁芰刻珡?qiáng),導(dǎo)致上行反射波不清楚。在這里估算品質(zhì)因子Q時(shí),采用下行波(圖2)。在求取Q的過程中,以第1層的Q值計(jì)算為例,分析時(shí)窗長(zhǎng)度、類型、噪聲等各種因素對(duì)Q值計(jì)算結(jié)果的影響。因?yàn)槊?道數(shù)據(jù)就可以計(jì)算出一個(gè)Q值,為了減少不必要的誤差,在分析結(jié)果時(shí)利用等效品質(zhì)因子值進(jìn)行比較。從計(jì)算結(jié)果(圖3~11)可以看出,在近道和遠(yuǎn)道存在誤差。近道誤差主要是因?yàn)槭艿秸鹪吹挠绊?,遠(yuǎn)道誤差是因?yàn)榻咏鼜?qiáng)吸收層,在2層分界面附近計(jì)算結(jié)果會(huì)受到一些影響。

    圖1 粘彈性VSP地震響應(yīng)

    圖2 下行波地震響應(yīng)

    3.1.1時(shí)窗長(zhǎng)度的影響

    等效品質(zhì)因子值為計(jì)算品質(zhì)因子值的平均值。對(duì)比圖3和圖4,時(shí)窗長(zhǎng)度為128 ms的等效品質(zhì)因子值為214.557 8,時(shí)窗長(zhǎng)度為256 ms的等效品質(zhì)因子值為206.390 6。對(duì)比結(jié)果說明,時(shí)窗長(zhǎng)度的適當(dāng)增加,對(duì)改進(jìn)的譜比法結(jié)果有一定的改善作用。

    注:加矩形時(shí)窗,時(shí)窗長(zhǎng)度128 ms,起止頻率為1~90 Hz圖3 時(shí)窗長(zhǎng)度為128 ms的等效Q值

    注:加矩形時(shí)窗,時(shí)窗長(zhǎng)度256 ms,起止頻率為1~90 Hz圖4 時(shí)窗長(zhǎng)度為256 ms的等效Q值

    3.1.2不同時(shí)窗函數(shù)的影響

    通過計(jì)算,加矩形窗的等效品質(zhì)因子值為206.390 6,加漢寧窗的值為226.013 3。對(duì)比圖4和圖5,說明時(shí)窗類型的改變對(duì)改進(jìn)的譜比法結(jié)果具有一定的影響,在計(jì)算時(shí)應(yīng)選擇合適的時(shí)窗函數(shù)。通過對(duì)各種窗函數(shù)的驗(yàn)證比較,因矩形窗主瓣比較集中,能量相對(duì)集中在主瓣,較接近真實(shí)的頻譜。在本章的計(jì)算中我們選擇矩形時(shí)窗。

    注:添加漢寧窗,時(shí)窗長(zhǎng)度為256 ms,起止頻率為1~90 Hz圖5 添加漢寧窗后的等效Q值

    3.1.3起止頻率的影響

    計(jì)算時(shí)起止頻率的選擇是多次試驗(yàn)總結(jié)得出的最佳頻率選擇范圍,圖4和圖6進(jìn)行比較說明,計(jì)算結(jié)果穩(wěn)定的前提是頻率范圍的選擇要合理。

    注:加矩形時(shí)窗,時(shí)窗長(zhǎng)度256 ms,起止頻率為1~70 Hz圖6 起止頻率為1~70 Hz的等效Q值

    3.1.4初始Q值的影響

    比較圖7和圖8可以看出,只有選擇合理的初始Q值,最終的計(jì)算結(jié)果才會(huì)比較穩(wěn)定。模型給出的理論Q值為200,該Q值說明地震波的衰減相對(duì)較小。當(dāng)理論模型為一般衰減層而給出的初始Q值為強(qiáng)衰減層時(shí),計(jì)算結(jié)果不穩(wěn)定;只有當(dāng)初始值與理論模型衰減程度相當(dāng)時(shí),計(jì)算結(jié)果才比較穩(wěn)定。

    注:時(shí)窗長(zhǎng)度256 ms,起止頻率為1~90 Hz,初始Q=314圖7 初始Q=314的等效Q值

    注:時(shí)窗長(zhǎng)度256 ms,起止頻率為1~90 Hz,初始Q=31.4圖8 初始Q=31.4的等效Q值

    3.1.5噪聲的影響

    圖9~10說明改進(jìn)的譜比法的計(jì)算結(jié)果隨地震資料中噪聲成分的增加越來越不穩(wěn)定,因此,改進(jìn)的譜比法對(duì)地震資料的品質(zhì)要求比較高。在沒有噪聲和干擾波干涉的情況下,該方法計(jì)算得到的等效品質(zhì)因子與品質(zhì)因子的理論值極其接近。當(dāng)?shù)卣鹩涗浿械碾S機(jī)噪聲較少時(shí),計(jì)算出來的Q值還是可靠的;當(dāng)?shù)卣鹌拭娴男旁氡容^低時(shí),計(jì)算結(jié)果就會(huì)變得很不穩(wěn)定,甚至還會(huì)出現(xiàn)負(fù)Q值的情況。這說明改進(jìn)的譜比法計(jì)算結(jié)果的準(zhǔn)確性在很大程度上依賴于地震資料信噪比的高低,因此,在進(jìn)行計(jì)算品質(zhì)因子之前,應(yīng)該對(duì)地震資料進(jìn)行去噪處理。

    注:加噪聲0.5%,時(shí)窗長(zhǎng)度為256 ms,起止頻率為1~90 Hz圖9 加噪聲0.5%的等效Q值

    注:加噪聲5%,時(shí)窗長(zhǎng)度256 ms,起止頻率為1~90 Hz圖10 加噪聲5%的等效Q值

    3.2改進(jìn)的譜比法的優(yōu)點(diǎn)

    短時(shí)傅立葉變換方法因?yàn)椴荒芨鶕?jù)頻率的變化進(jìn)行時(shí)窗大小的調(diào)節(jié),在使用時(shí)受到限制。S變換[10]采用的是可變高斯窗函數(shù),有多尺度聚焦性,可以保持頻率的絕對(duì)相位,其基本變換函數(shù)可以不滿足容許性條件,并且可以實(shí)現(xiàn)地震信號(hào)在時(shí)間域和頻率域的無損轉(zhuǎn)換。但是S變換基本函數(shù)的形態(tài)是固定不變的,導(dǎo)致其在應(yīng)用中受到限制,為此學(xué)者們將S變換推廣到廣義S變換??紤]到傅立葉變換的不足和廣義S變換的優(yōu)點(diǎn),在這里我們將改進(jìn)的譜比法中的傅立葉變換改為廣義S變換。

    (14)

    其中:k是大于零的常數(shù)。

    在廣義S變換中,可以通過控制k值來提高時(shí)間-頻率譜(以下簡(jiǎn)稱“時(shí)-頻譜”)的時(shí)間(或頻率)分辨率。但是由于k是常數(shù),所以它若提高了頻率分辨率,那么一定是以犧牲時(shí)間分辨率為代價(jià)的,反之亦然。從總體上來說,S變換得到的時(shí)-頻譜是基于時(shí)頻分辨率不變的前提;雖然在廣義S變換中,調(diào)節(jié)k值可以得到相應(yīng)的多個(gè)時(shí)間頻率譜,但當(dāng)k值固定不變的時(shí)候,由廣義S變換得到的時(shí)-頻譜分辨率具有唯一的變化趨勢(shì)。為了使時(shí)頻譜能夠更好地反映全時(shí)段內(nèi)所有信號(hào)的局部頻段的細(xì)節(jié),我們引入一個(gè)和頻率有關(guān)的可變因子σf,對(duì)Gauss窗函數(shù)進(jìn)行改進(jìn)[12],那么Gauss窗函數(shù)可以修正為:

    (15)

    (16)

    在這里,我們將σf定義為σf(f)=(kf+b)a。

    本文使用的改進(jìn)的廣義S變換方法主要是針對(duì)非平穩(wěn)信號(hào)時(shí)間頻率的特征,該方法含有一個(gè)可變的因子函數(shù),函數(shù)的自變量是頻率。使用時(shí)窗函數(shù)的類型以及其他參數(shù)的調(diào)節(jié)都是靈活多變的,從而可以避免S變換和廣義S變換中時(shí)頻分辨率的變化趨勢(shì)問題,增強(qiáng)其適應(yīng)性和針對(duì)性[12]。

    注:廣義S變換,起止頻率為1~90 Hz 圖11 相鄰檢波器之間的Q值

    對(duì)比圖4中的層Q值206.390 6和圖11中的層Q值197.074 9,可以看出,改進(jìn)的譜比法中利用廣義S變換進(jìn)行時(shí)頻變換計(jì)算得到的地層Q值與理論值誤差更小。

    利用改進(jìn)的譜比法對(duì)3.1節(jié)用的粘彈性VSP地震記錄提取品質(zhì)因子,對(duì)波場(chǎng)分離得到的下行波場(chǎng)進(jìn)行等效品質(zhì)因子的提取(圖12)。提取結(jié)果與理論值以及誤差見表1。

    圖12 估計(jì)Q值與理論Q值的比較

    層數(shù)123理論Q值2003040基于廣義S變換改進(jìn)的譜比法201.826230.902640.9690基于傅立葉變換改進(jìn)的譜比法207.616025.307948.1033基于廣義S變換的計(jì)算誤差0.9%3.0%2.4%基于傅里葉變換的計(jì)算誤差3.8%15.6%20.3%

    通過表1和圖12可以看出,基于廣義S變換的改進(jìn)的譜比法提取的Q值與理論Q值匹配效果較好,而基于傅立葉變換的改進(jìn)的譜比法提取的Q值相對(duì)較差。不過總體而言,對(duì)于低品質(zhì)因子的層估計(jì)的Q值誤差較大,品質(zhì)因子較大的層估算的Q值誤差相對(duì)較小,主要原因?yàn)镼值小的地層吸收大,Q值較大的地層吸收小。

    基于上述結(jié)論,在下面的實(shí)例中,采用基于廣義S變換的改進(jìn)的譜比法進(jìn)行品質(zhì)因子的估算。

    3.3改進(jìn)譜比法的應(yīng)用實(shí)例

    選用勝利油田某工區(qū)采集到的VSP炮記錄進(jìn)行測(cè)試。具體采集參數(shù)為:采用零井源距接收觀測(cè)系統(tǒng),震源的沉放深度為15 m,接收道數(shù)為315道,道間距為5 m,第一個(gè)檢波點(diǎn)位于230 m,采樣點(diǎn)為2 000,時(shí)間采樣間隔為0.5 mm(圖13)。由于地震資料在采集中受到很多因素的影響,例如多次波和自由噪聲等,對(duì)該炮記錄進(jìn)行波場(chǎng)分離,得到上行波場(chǎng)和下行場(chǎng)(圖14)。為了提取更加準(zhǔn)確的等效品質(zhì)因子,首先對(duì)下行波場(chǎng)進(jìn)行去除噪聲處理,采用改進(jìn)的譜比法進(jìn)行提取(圖15)。由圖15看出,該資料存在強(qiáng)吸收層和弱吸收層。在接近地表部位為強(qiáng)吸收層,這與我們常識(shí)認(rèn)知是相符的,因此該計(jì)算結(jié)果具有一定的可靠性。

    圖13 VSP炮記錄

    a 下行波場(chǎng)b 上行波場(chǎng)圖14 分離得到的波場(chǎng)

    圖15 提取的品質(zhì)因子

    上述實(shí)例表明,利用改進(jìn)的譜比法對(duì)VSP地震資料進(jìn)行研究,提取出的品質(zhì)因子計(jì)算精度得到了提高。譜比法只是利用VSP資料的直達(dá)波進(jìn)行計(jì)算,改進(jìn)的譜比法利用所有下行反射波進(jìn)行計(jì)算振幅譜,利用改進(jìn)的譜比法進(jìn)行計(jì)算的結(jié)果更加真實(shí)可靠。

    4 討論與結(jié)論

    針對(duì)振幅衰減法和解析信號(hào)法在計(jì)算時(shí)需要真振幅數(shù)據(jù),而在實(shí)際采集中,真振幅數(shù)據(jù)很難獲得,以及譜比法計(jì)算時(shí)目的層的厚度隨時(shí)間間隔的增大而變大,計(jì)算結(jié)果也逐漸趨近于一個(gè)常數(shù)的問題,利用改進(jìn)的譜比法求取介質(zhì)Q值,更符合實(shí)際需要。

    根據(jù)理論分析和實(shí)際VSP地震記錄數(shù)據(jù)計(jì)算表明:

    1)改進(jìn)的譜比法充分考慮了地層對(duì)地震波的吸收作用,使得求得的Q值更符合實(shí)際情況。

    2)改進(jìn)的譜比法中利用廣義的S變換求取Q值,使得得到的Q值的物理意義更加清晰。

    3)采用改進(jìn)的譜比法求取品質(zhì)因子時(shí),綜合考慮了算法的一些影響因素,對(duì)算法進(jìn)行了相應(yīng)的改進(jìn),使得高頻段的分辨率得到了提高,Q值的計(jì)算精度也相應(yīng)提高。

    當(dāng)然Q值計(jì)算精度的提高不只是依賴于算法本身的計(jì)算精度和地震資料的高信噪比,在實(shí)際資料處理時(shí),應(yīng)該對(duì)具體問題進(jìn)行具體分析。

    [1]劉祖沅. 巖石中地震波衰減的實(shí)驗(yàn)研究[J]. 地震科學(xué)研究. 1983, 6:34-40.

    [2]Rainer Tonn. The deternation of the seismic quality factor Q from VSP data:A comparison of different computational methods[J].Geophysical Prospecting, 1991,39(1):1-27.

    [3]Aki k,Richards P G. Quantitative Seismology: Theory and Methods[M]. San Francisco: W H Freeman and Company,1980.

    [4]王慧茹. 譜模擬頻譜比法反演介質(zhì)品質(zhì)因子穩(wěn)定性分析[J]. 物探與化探,2009, 33(1):85-87.

    [5]孫成禹. 譜模擬方法及其在提高地震資料分辨率中的應(yīng)用[J]. 石油地球物理勘探, 2000, 35(1):27-35.

    [6]Ricker N H. Transient waves in visco-elastic media [M]. New York: Elsevier Science Publishing Company,1977.

    [7]宮同舉, 孫成禹, 彭洪超,等. 幾種提取品質(zhì)因子方法的對(duì)比分析[J]. 勘探地球物理進(jìn)展, 2009, 32(4):252-257.

    [8]吳琳. 品質(zhì)因子求取與補(bǔ)償方法研究[D]. 山東:中國石油大學(xué)(華東),2011.

    [9]Blias F. Q-factor estimation through optimization approach to near-offset VSP data[J]. SEG Technical Program Expanded Abstracts , 2011:4278-4282.

    [10]Stockwell R G,Mansinha L, Lowe R P.Localization of the complex spectrum:The S transform[J].IEEE transactions on Signal Process,1996,44(4):998-1001.

    [11]Mansinha L, Stockwell R G. Local S-spectrum analysis of 1-D and 2-D data[J]. Physics of the Earth and Plantetary Interiors, 1997, 103(3):329-336.

    [12]周竹生, 陳友良. 含可變因子的廣義S變換及其時(shí)頻濾波[J]. 煤田地質(zhì)與勘探, 2011, 39(6):63-67.

    Extraction Method of Quality Factor based on VSP Seismic Data

    DONG Bo1,2, JI Chun-ling1, ZHOU An-pin1, LI Feng1,ZHANG Huan-xi1, ZHAO Yu-chen1, LIU Jing1, NIU Shu-yu1

    (1. Hebei GEO University, Shijiazhuang 050031,China;2. Shijiazhuang Central Seismic Station, Earthquake Administration of Hebei province, Shijiazhuang 050021,China)

    This paper studies VSP data in shengli oilfield with advanced spectral ratio method to improve the calculation accuracy of quality factors. Different from traditional spectral ratio method in which only direct waves from VSP data are used to calculate, the advanced one takes all downgoing reflected waves to calculate amplitude spectrum, whose results will be more authentic and reliable.

    quality factor; improved spectral ratio method; attenuation; Fourier transform; influencing factors

    2016-02-16

    2015年度震情跟蹤定向工作任務(wù)“利用背景噪聲監(jiān)測(cè)波速變化系統(tǒng)初步研究”(2015010117)

    董博(1986—),男,河北鹿泉人,助理工程師,現(xiàn)主要從事地震監(jiān)測(cè)與分析預(yù)報(bào)工作.E-mail:dongbo0002@126.com

    P315.6

    A

    1003-1375(2016)03-0061-07

    10.3969/j.issn.1003-1375.2016.03.011

    董博,紀(jì)春玲,周安聘,等.基于VSP地震資料的品質(zhì)因子提取方法[J].華北地震科學(xué),2016,34(3):61-67.

    猜你喜歡
    時(shí)窗比法計(jì)算結(jié)果
    化虛為實(shí) 觸摸物理——物理方法之類比法
    加權(quán)譜比法Q值估計(jì)
    物理方法之類比法
    GRAPES-GFS模式2 m溫度預(yù)報(bào)的最優(yōu)時(shí)窗滑動(dòng)訂正方法
    氣象科技(2021年5期)2021-11-02 01:36:28
    一種基于改進(jìn)時(shí)窗法的爆炸沖擊波檢測(cè)方法
    最好的比較
    不等高軟橫跨橫向承力索計(jì)算及計(jì)算結(jié)果判斷研究
    甘肅科技(2020年20期)2020-04-13 00:30:40
    不同介入時(shí)窗和療程對(duì)高壓氧治療急性脊髓損傷的影響研究
    UWB搜救生物雷達(dá)分段時(shí)窗探測(cè)系統(tǒng)的實(shí)現(xiàn)
    超壓測(cè)試方法對(duì)炸藥TNT當(dāng)量計(jì)算結(jié)果的影響
    日本免费a在线| 色哟哟·www| 白带黄色成豆腐渣| 精品无人区乱码1区二区| 观看美女的网站| 琪琪午夜伦伦电影理论片6080| 国产91精品成人一区二区三区| 男女做爰动态图高潮gif福利片| 日韩成人在线观看一区二区三区| 人人妻人人澡欧美一区二区| 国产亚洲精品久久久久久毛片| 成人国产综合亚洲| 婷婷亚洲欧美| 两性午夜刺激爽爽歪歪视频在线观看| a级一级毛片免费在线观看| 在线播放无遮挡| 久久精品国产亚洲av天美| 老司机午夜十八禁免费视频| 麻豆一二三区av精品| 中文字幕人妻熟人妻熟丝袜美| 亚洲国产精品久久男人天堂| 成人美女网站在线观看视频| 国产成人a区在线观看| 中文字幕久久专区| 91九色精品人成在线观看| 欧美日韩国产亚洲二区| 少妇的逼水好多| 中文字幕高清在线视频| 在线播放无遮挡| 国产精品国产高清国产av| 可以在线观看毛片的网站| 久久国产精品人妻蜜桃| 一级毛片久久久久久久久女| 久久久久性生活片| 在线免费观看的www视频| 欧美精品国产亚洲| 美女xxoo啪啪120秒动态图 | 日本免费a在线| 久久精品国产清高在天天线| 久久国产乱子免费精品| 欧美一区二区精品小视频在线| av福利片在线观看| 伊人久久精品亚洲午夜| 日日夜夜操网爽| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 日韩欧美在线乱码| 午夜日韩欧美国产| 波多野结衣巨乳人妻| 成人一区二区视频在线观看| 真人做人爱边吃奶动态| 91九色精品人成在线观看| 国产精品98久久久久久宅男小说| 麻豆久久精品国产亚洲av| 亚洲不卡免费看| 天天一区二区日本电影三级| 精品人妻一区二区三区麻豆 | 亚洲电影在线观看av| 三级男女做爰猛烈吃奶摸视频| 久久国产乱子伦精品免费另类| 九九热线精品视视频播放| 嫁个100分男人电影在线观看| 免费黄网站久久成人精品 | 亚洲精华国产精华精| 亚洲欧美精品综合久久99| 在线国产一区二区在线| 久久精品人妻少妇| 可以在线观看的亚洲视频| 一卡2卡三卡四卡精品乱码亚洲| 9191精品国产免费久久| 精品人妻1区二区| 亚洲内射少妇av| 男插女下体视频免费在线播放| 亚洲av不卡在线观看| 国产精品不卡视频一区二区 | 中出人妻视频一区二区| 久久精品国产亚洲av涩爱 | 日本三级黄在线观看| 欧美激情国产日韩精品一区| 成人三级黄色视频| 91在线精品国自产拍蜜月| xxxwww97欧美| 亚洲,欧美,日韩| 男人舔奶头视频| 天堂√8在线中文| 午夜福利视频1000在线观看| 天堂影院成人在线观看| 亚洲电影在线观看av| 久久精品国产清高在天天线| 日韩欧美精品免费久久 | 精品午夜福利在线看| 亚洲五月婷婷丁香| 男女那种视频在线观看| 国产精品一区二区性色av| 深夜精品福利| 国产淫片久久久久久久久 | 亚洲美女视频黄频| 99riav亚洲国产免费| 日本在线视频免费播放| 亚洲电影在线观看av| 黄片小视频在线播放| 欧美xxxx黑人xx丫x性爽| 亚洲不卡免费看| 亚洲色图av天堂| 欧美在线黄色| 老鸭窝网址在线观看| 日韩欧美在线二视频| 网址你懂的国产日韩在线| 亚洲自拍偷在线| 欧美性感艳星| 国产精品久久电影中文字幕| 麻豆成人av在线观看| 国产午夜精品久久久久久一区二区三区 | 亚洲最大成人av| 深夜a级毛片| 成人三级黄色视频| 精品国产三级普通话版| 在线免费观看的www视频| 90打野战视频偷拍视频| 国产黄色小视频在线观看| 夜夜爽天天搞| 一级作爱视频免费观看| 国产精品亚洲av一区麻豆| 别揉我奶头 嗯啊视频| 中文字幕久久专区| 日日摸夜夜添夜夜添av毛片 | 午夜影院日韩av| 国产精品精品国产色婷婷| 国产精品美女特级片免费视频播放器| 韩国av一区二区三区四区| 99国产极品粉嫩在线观看| 久久久久久大精品| 国产69精品久久久久777片| 欧美日本视频| 亚洲欧美日韩无卡精品| 一区二区三区免费毛片| 美女被艹到高潮喷水动态| 99久久久亚洲精品蜜臀av| 国产一区二区激情短视频| 美女 人体艺术 gogo| 直男gayav资源| 人人妻人人看人人澡| 在线天堂最新版资源| 最好的美女福利视频网| 国产精品女同一区二区软件 | 精品久久久久久成人av| 欧美国产日韩亚洲一区| 成人美女网站在线观看视频| 特级一级黄色大片| 国产亚洲精品久久久久久毛片| 午夜福利在线在线| 中文字幕人成人乱码亚洲影| 国产精品,欧美在线| 亚洲最大成人av| 久久久久九九精品影院| 久99久视频精品免费| 精品99又大又爽又粗少妇毛片 | 简卡轻食公司| 国语自产精品视频在线第100页| 国产精品国产高清国产av| 三级国产精品欧美在线观看| 色噜噜av男人的天堂激情| 香蕉av资源在线| 久久久久久国产a免费观看| 一个人免费在线观看电影| 青草久久国产| 国产亚洲欧美98| 国产成人a区在线观看| 国产精品99久久久久久久久| 国产精品亚洲一级av第二区| 午夜福利高清视频| 久久这里只有精品中国| 国产精品爽爽va在线观看网站| 亚洲欧美激情综合另类| 国产精品一区二区免费欧美| 身体一侧抽搐| 99国产综合亚洲精品| 如何舔出高潮| 动漫黄色视频在线观看| 偷拍熟女少妇极品色| 91麻豆av在线| 精品不卡国产一区二区三区| 国产高清视频在线观看网站| 国产单亲对白刺激| 日本黄色视频三级网站网址| 中文字幕熟女人妻在线| 日本三级黄在线观看| 动漫黄色视频在线观看| 免费观看人在逋| 亚洲avbb在线观看| 九九在线视频观看精品| 中文字幕熟女人妻在线| 天堂影院成人在线观看| 婷婷精品国产亚洲av在线| 观看美女的网站| 欧美国产日韩亚洲一区| 麻豆久久精品国产亚洲av| or卡值多少钱| 精品午夜福利视频在线观看一区| 亚洲第一区二区三区不卡| 欧美一区二区国产精品久久精品| 国产成人啪精品午夜网站| 成年女人看的毛片在线观看| 欧美日韩乱码在线| 亚洲欧美日韩高清在线视频| 久久久久亚洲av毛片大全| 欧美激情在线99| 国产av一区在线观看免费| 1000部很黄的大片| 69av精品久久久久久| 狂野欧美白嫩少妇大欣赏| 亚洲av免费高清在线观看| 精品人妻一区二区三区麻豆 | 夜夜夜夜夜久久久久| 亚洲色图av天堂| www.999成人在线观看| 国产麻豆成人av免费视频| 91狼人影院| 国产高清激情床上av| 免费观看的影片在线观看| 可以在线观看毛片的网站| 亚洲欧美日韩无卡精品| 久9热在线精品视频| 国产精品人妻久久久久久| 18禁黄网站禁片免费观看直播| 日本 欧美在线| 十八禁国产超污无遮挡网站| 日韩精品青青久久久久久| 日本a在线网址| 亚洲国产精品久久男人天堂| 我的老师免费观看完整版| 九色成人免费人妻av| 精品无人区乱码1区二区| 亚洲欧美日韩东京热| 亚洲激情在线av| 国产真实乱freesex| 欧美一级a爱片免费观看看| 每晚都被弄得嗷嗷叫到高潮| 久久久精品欧美日韩精品| 国产精品伦人一区二区| 动漫黄色视频在线观看| 岛国在线免费视频观看| 成人精品一区二区免费| 香蕉av资源在线| 精品无人区乱码1区二区| 看十八女毛片水多多多| 久久久久久久午夜电影| 丰满人妻熟妇乱又伦精品不卡| 999久久久精品免费观看国产| 变态另类丝袜制服| 日本在线视频免费播放| 亚洲美女黄片视频| 久久九九热精品免费| 一本综合久久免费| 国产高清视频在线观看网站| 欧美另类亚洲清纯唯美| 亚洲精品亚洲一区二区| 国产精华一区二区三区| 欧美+日韩+精品| 亚洲av.av天堂| 十八禁人妻一区二区| 国产精品久久电影中文字幕| 日韩欧美国产一区二区入口| 18禁裸乳无遮挡免费网站照片| 赤兔流量卡办理| 热99在线观看视频| 国产人妻一区二区三区在| 亚洲最大成人中文| 亚洲av日韩精品久久久久久密| 日韩欧美在线二视频| 亚洲精品亚洲一区二区| 嫁个100分男人电影在线观看| 国产大屁股一区二区在线视频| 波多野结衣巨乳人妻| 十八禁人妻一区二区| 欧美成人免费av一区二区三区| 99精品在免费线老司机午夜| 久久婷婷人人爽人人干人人爱| 国产精品美女特级片免费视频播放器| 国产精品亚洲一级av第二区| av视频在线观看入口| 天堂网av新在线| 日韩中字成人| 日韩亚洲欧美综合| 国产乱人伦免费视频| 欧美日韩中文字幕国产精品一区二区三区| 尤物成人国产欧美一区二区三区| 亚洲精品亚洲一区二区| 夜夜爽天天搞| 国内久久婷婷六月综合欲色啪| 少妇人妻精品综合一区二区 | 搡女人真爽免费视频火全软件 | 国产成人福利小说| 国产成人啪精品午夜网站| 成人国产一区最新在线观看| 性色avwww在线观看| 欧美日本视频| 精品久久久久久久久av| 久久久久精品国产欧美久久久| 大型黄色视频在线免费观看| 亚洲精品久久国产高清桃花| 91九色精品人成在线观看| 亚洲美女黄片视频| 久99久视频精品免费| 在线观看66精品国产| 床上黄色一级片| 蜜桃久久精品国产亚洲av| 日本免费一区二区三区高清不卡| 99riav亚洲国产免费| 99国产精品一区二区三区| 麻豆国产97在线/欧美| 亚洲av不卡在线观看| 国产三级中文精品| 久久久国产成人免费| 国产精品乱码一区二三区的特点| 亚洲狠狠婷婷综合久久图片| 婷婷六月久久综合丁香| 99热这里只有是精品在线观看 | 国产麻豆成人av免费视频| 国产视频一区二区在线看| 国内毛片毛片毛片毛片毛片| 国产伦精品一区二区三区四那| 亚洲综合色惰| 夜夜夜夜夜久久久久| 国产免费一级a男人的天堂| 久久6这里有精品| 国产精品人妻久久久久久| 在线免费观看的www视频| 亚洲av电影不卡..在线观看| 亚洲自拍偷在线| 久99久视频精品免费| 九九在线视频观看精品| 最近最新免费中文字幕在线| 久久精品国产亚洲av涩爱 | 黄色日韩在线| 精品免费久久久久久久清纯| 亚洲国产日韩欧美精品在线观看| 床上黄色一级片| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 一级作爱视频免费观看| 精品国产亚洲在线| 国产成人影院久久av| 又黄又爽又免费观看的视频| 欧美一区二区亚洲| 中文字幕av在线有码专区| 国产高清视频在线观看网站| 久久国产精品人妻蜜桃| 香蕉av资源在线| 成熟少妇高潮喷水视频| 精品人妻视频免费看| 可以在线观看的亚洲视频| 亚洲国产高清在线一区二区三| 一区二区三区四区激情视频 | 97热精品久久久久久| 熟妇人妻久久中文字幕3abv| av在线天堂中文字幕| 窝窝影院91人妻| 99热只有精品国产| 尤物成人国产欧美一区二区三区| 国产伦精品一区二区三区视频9| 桃色一区二区三区在线观看| 99国产极品粉嫩在线观看| 国产人妻一区二区三区在| 欧美三级亚洲精品| 男人的好看免费观看在线视频| 人妻制服诱惑在线中文字幕| 天美传媒精品一区二区| 欧美乱色亚洲激情| 日韩欧美免费精品| 国产精品野战在线观看| 精品久久国产蜜桃| 久久久久久久亚洲中文字幕 | 国产精品久久久久久久久免 | 久9热在线精品视频| 国产精品久久久久久久电影| 18禁裸乳无遮挡免费网站照片| 欧美乱色亚洲激情| 一区福利在线观看| 国内毛片毛片毛片毛片毛片| 99视频精品全部免费 在线| 黄色女人牲交| 国产美女午夜福利| 真实男女啪啪啪动态图| 国产伦在线观看视频一区| 久久国产乱子免费精品| 黄色视频,在线免费观看| 欧美成人免费av一区二区三区| 在线观看免费视频日本深夜| 欧美性猛交黑人性爽| 男人狂女人下面高潮的视频| 久久香蕉精品热| 美女高潮喷水抽搐中文字幕| 少妇人妻精品综合一区二区 | 人妻久久中文字幕网| 精品国产亚洲在线| 深爱激情五月婷婷| 尤物成人国产欧美一区二区三区| 久久久久久九九精品二区国产| 国产伦人伦偷精品视频| 黄色丝袜av网址大全| 又粗又爽又猛毛片免费看| 男人的好看免费观看在线视频| 成年版毛片免费区| 欧美黑人巨大hd| 在线十欧美十亚洲十日本专区| 黄色视频,在线免费观看| 日韩欧美在线乱码| 国产精品精品国产色婷婷| 9191精品国产免费久久| 亚洲av电影在线进入| 国产精品野战在线观看| 欧美性猛交黑人性爽| 精品久久久久久久人妻蜜臀av| 国产黄a三级三级三级人| 亚洲av成人精品一区久久| 久久中文看片网| 国产大屁股一区二区在线视频| 中文字幕人成人乱码亚洲影| 国产真实乱freesex| 国产探花在线观看一区二区| 成人性生交大片免费视频hd| 人人妻人人看人人澡| 亚洲成人精品中文字幕电影| 日日摸夜夜添夜夜添小说| 国内精品美女久久久久久| 淫秽高清视频在线观看| 直男gayav资源| 国产主播在线观看一区二区| 久久久久精品国产欧美久久久| 99久久精品热视频| 午夜久久久久精精品| 久久婷婷人人爽人人干人人爱| 别揉我奶头 嗯啊视频| 九九久久精品国产亚洲av麻豆| 日本在线视频免费播放| 国产精品一及| 日本一本二区三区精品| 99在线人妻在线中文字幕| 成人欧美大片| 国产精品日韩av在线免费观看| 欧美潮喷喷水| 女生性感内裤真人,穿戴方法视频| 高潮久久久久久久久久久不卡| 欧美精品国产亚洲| 真人做人爱边吃奶动态| 欧美性猛交╳xxx乱大交人| 一本综合久久免费| 免费看美女性在线毛片视频| 757午夜福利合集在线观看| 丝袜美腿在线中文| 小说图片视频综合网站| 亚洲欧美清纯卡通| 精品国产亚洲在线| 1024手机看黄色片| 97碰自拍视频| 能在线免费观看的黄片| 成人高潮视频无遮挡免费网站| 国产一区二区亚洲精品在线观看| 免费一级毛片在线播放高清视频| 亚洲狠狠婷婷综合久久图片| 久久精品国产亚洲av香蕉五月| 国产久久久一区二区三区| 国产欧美日韩一区二区三| 九色国产91popny在线| 精品午夜福利在线看| 脱女人内裤的视频| 亚洲av日韩精品久久久久久密| 9191精品国产免费久久| 91在线精品国自产拍蜜月| 亚洲av五月六月丁香网| 在线观看美女被高潮喷水网站 | 亚洲无线观看免费| 精品一区二区免费观看| 精品人妻熟女av久视频| 国产精品一区二区三区四区免费观看 | 亚洲成a人片在线一区二区| 国内精品美女久久久久久| 亚洲五月天丁香| 99久久精品一区二区三区| 久久99热这里只有精品18| 性欧美人与动物交配| 我的老师免费观看完整版| 欧美黑人巨大hd| 欧美+日韩+精品| 欧美性猛交黑人性爽| 亚洲av成人精品一区久久| 国产精品久久久久久精品电影| 在线看三级毛片| 国产精品98久久久久久宅男小说| 国内毛片毛片毛片毛片毛片| 99久久精品一区二区三区| 日韩欧美一区二区三区在线观看| 国产视频内射| 三级国产精品欧美在线观看| 欧美成人a在线观看| 欧美日本视频| 欧美黄色片欧美黄色片| 欧美极品一区二区三区四区| 中文字幕av成人在线电影| 午夜福利视频1000在线观看| 久久亚洲真实| 久久久久久久久中文| 国产国拍精品亚洲av在线观看| 99热这里只有精品一区| 欧美日本亚洲视频在线播放| 亚洲成人久久爱视频| 女同久久另类99精品国产91| 国产黄色小视频在线观看| a级一级毛片免费在线观看| 久久精品国产清高在天天线| 夜夜夜夜夜久久久久| 成人毛片a级毛片在线播放| 久久久久国内视频| 757午夜福利合集在线观看| 一个人免费在线观看的高清视频| 九色国产91popny在线| 在线观看一区二区三区| 麻豆成人午夜福利视频| 在线观看免费视频日本深夜| 亚洲av熟女| 国产精品不卡视频一区二区 | 极品教师在线视频| 久久久午夜欧美精品| 日本欧美国产在线视频| av福利片在线观看| 成人漫画全彩无遮挡| 赤兔流量卡办理| 九色成人免费人妻av| 国产成人a区在线观看| 丰满人妻一区二区三区视频av| 亚洲不卡免费看| 高清日韩中文字幕在线| av在线老鸭窝| 国产成人91sexporn| 国产淫片久久久久久久久| 国产熟女欧美一区二区| 男插女下体视频免费在线播放| 免费观看av网站的网址| 亚洲欧美日韩卡通动漫| 欧美+日韩+精品| 国产精品精品国产色婷婷| 亚洲无线观看免费| 黄色怎么调成土黄色| 久久久久国产网址| 免费播放大片免费观看视频在线观看| 人妻系列 视频| 亚洲av国产av综合av卡| 午夜激情福利司机影院| 99视频精品全部免费 在线| 欧美xxxx黑人xx丫x性爽| 国产毛片a区久久久久| 日韩强制内射视频| 国产在线男女| 深爱激情五月婷婷| 国产成人freesex在线| 日本爱情动作片www.在线观看| 黄色一级大片看看| 精品久久久久久久末码| 夫妻性生交免费视频一级片| 美女内射精品一级片tv| 三级经典国产精品| 在线观看免费高清a一片| 波野结衣二区三区在线| 一级毛片我不卡| 大话2 男鬼变身卡| av卡一久久| 国产有黄有色有爽视频| 国内精品宾馆在线| 伊人久久精品亚洲午夜| 成人特级av手机在线观看| 极品教师在线视频| 日韩精品有码人妻一区| 简卡轻食公司| 国产免费一级a男人的天堂| 日韩亚洲欧美综合| 国产免费福利视频在线观看| 国产精品一区www在线观看| 午夜福利网站1000一区二区三区| 男插女下体视频免费在线播放| 亚洲欧美精品自产自拍| 久久久久精品性色| 交换朋友夫妻互换小说| 麻豆成人午夜福利视频| 国产有黄有色有爽视频| 最后的刺客免费高清国语| 夫妻性生交免费视频一级片| 91精品伊人久久大香线蕉| 亚洲国产高清在线一区二区三| 美女cb高潮喷水在线观看| 欧美精品一区二区大全| 亚洲精品第二区| 久久热精品热| 久久久久网色| 免费高清在线观看视频在线观看| 高清欧美精品videossex| 人人妻人人澡人人爽人人夜夜| 一区二区三区免费毛片| 伦精品一区二区三区| 嫩草影院新地址| 91精品一卡2卡3卡4卡| 麻豆国产97在线/欧美| 欧美97在线视频| 亚洲欧洲日产国产| 久久久a久久爽久久v久久| 自拍欧美九色日韩亚洲蝌蚪91 | 亚洲精品日韩av片在线观看| 久热久热在线精品观看| 亚洲精品成人久久久久久| 香蕉精品网在线| 亚洲av成人精品一二三区| 中文字幕制服av| 精品国产三级普通话版| 久久久久久久亚洲中文字幕| 久久人人爽人人爽人人片va| 联通29元200g的流量卡| 嫩草影院精品99|