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

    偽隨機(jī)編碼源激發(fā)下的時域電磁信號合成

    2016-12-07 08:12:51王若王妙月底青云薛國強(qiáng)殷長春雷達(dá)
    地球物理學(xué)報 2016年12期
    關(guān)鍵詞:信號

    王若, 王妙月, 底青云, 薛國強(qiáng), 殷長春, 雷達(dá)

    1 中國科學(xué)院頁巖氣與地質(zhì)工程重點實驗室,中國科學(xué)院地質(zhì)與地球物理研究所, 北京 100029 2 吉林大學(xué)地球探測科學(xué)與技術(shù)學(xué)院, 長春 130026

    ?

    偽隨機(jī)編碼源激發(fā)下的時域電磁信號合成

    王若1, 王妙月1, 底青云1, 薛國強(qiáng)1, 殷長春2, 雷達(dá)1

    1 中國科學(xué)院頁巖氣與地質(zhì)工程重點實驗室,中國科學(xué)院地質(zhì)與地球物理研究所, 北京 100029 2 吉林大學(xué)地球探測科學(xué)與技術(shù)學(xué)院, 長春 130026

    將偽隨機(jī)編碼技術(shù)引入到人工源電磁法后,可以通過加大發(fā)射功率以及應(yīng)用后續(xù)的相關(guān)處理技術(shù)來達(dá)到壓制噪聲、加大探測深度及提高分辨率的目的,因此引起了越來越多學(xué)者的關(guān)注及研究,但大多數(shù)研究集中在資料處理的相關(guān)技術(shù)上,對模擬電磁信號關(guān)注較少.然而,資料處理工作大多是從電磁信號出發(fā)的,模擬偽隨機(jī)編碼源激發(fā)下的電磁信號不但可以為資料處理環(huán)節(jié)提供理論數(shù)據(jù),而且可以為檢測資料處理的效果提供中間結(jié)果,因此,模擬電磁信號工作必不可少.本文根據(jù)獲得接收信號的物理過程來實現(xiàn)偽隨機(jī)編碼源激發(fā)下的電磁信號合成.首先用解析公式獲得特定地電結(jié)構(gòu)的大地頻率域響應(yīng),然后通過余弦變換得到時間域階躍響應(yīng),接下來用階躍響應(yīng)的時間導(dǎo)數(shù)得到大地脈沖響應(yīng),通過將大地脈沖響應(yīng)與偽隨機(jī)編碼源的褶積得到理想接收信號,最后,用低通濾波器來模擬發(fā)射設(shè)備和接收設(shè)備的頻帶限制,將之和噪聲一起加到理想接收信號上,最終模擬出仿真的合成信號.通過和野外實際接收信號對比發(fā)現(xiàn)本文合成信號仿真度較高,可以服務(wù)于后續(xù)的數(shù)據(jù)處理環(huán)節(jié).

    偽隨機(jī)編碼技術(shù); 時間域; 大地脈沖響應(yīng); 接收信號; 合成

    1 引言

    近年來,傳統(tǒng)的電磁方法在能源、礦產(chǎn)資源、水資源、環(huán)境地質(zhì)及工程地質(zhì)勘察中都得到了廣泛應(yīng)用且發(fā)揮了重要作用(底青云等,2006;Xue et al., 2014).隨著國民經(jīng)濟(jì)的發(fā)展,需要探測埋深更深、精度更高的目標(biāo)體時,這些方法的應(yīng)用受到了一定的限制,亟待開發(fā)新的技術(shù)來提高探測精度與探測深度.

    近年來,一些新技術(shù)新方法被陸續(xù)提出來,英國愛丁堡大學(xué)的Wright(2002)提出了多道瞬變電磁法的概念和探測油氣目標(biāo)體的處理技術(shù),何繼善(2010)提出了可提高分辨率及加大勘探深度的廣域電磁法及相關(guān)技術(shù).

    在這些新方法中,都引入了偽隨機(jī)編碼發(fā)射技術(shù),其全新的處理方式及應(yīng)用效果引起各個研究機(jī)構(gòu)的關(guān)注.近幾年,中國科學(xué)院、中南大學(xué)、中國地質(zhì)大學(xué)、國土資源部以及相關(guān)單位對這種發(fā)射技術(shù)(趙璧如等,2006;羅先中等,2014;羅延鐘等,2015)及后續(xù)的處理工作都進(jìn)行了不同程度的研究.湯井田和羅維斌(2008)分析了逆重復(fù)m序列的波形特征、相關(guān)函數(shù)及頻譜特性,討論了對各種干擾的壓制方法及測量精度,薛國強(qiáng)等(2015)分析了多通道瞬變電磁法(MTEM)偽隨機(jī)編碼發(fā)射的響應(yīng)特性和方法的關(guān)鍵技術(shù),齊彥福等(2015)利用方波響應(yīng)移位疊加和電流導(dǎo)數(shù)與階躍響應(yīng)褶積兩種方法實現(xiàn)了理論m序列和實際發(fā)射波形的全時正演模擬,武欣等(2015)研究了m序列偽隨機(jī)編碼源電磁響應(yīng)的精細(xì)辨識問題.

    學(xué)者多把研究重點放在相關(guān)辨識及對噪聲的壓制上,對相關(guān)辨識的方法理論做了充分的探討,但在模擬電磁信號方面卻著墨較少.然而,采用偽隨機(jī)編碼源的資料處理工作大多是從電磁信號出發(fā)的,模擬偽隨機(jī)編碼源激發(fā)下的電磁信號不但可以為資料處理環(huán)節(jié)提供理論數(shù)據(jù),而且可以為檢測資料處理的質(zhì)量提供中間結(jié)果,因此,模擬電磁信號工作是必不可少的.

    本文根據(jù)獲得接收信號的物理過程來實現(xiàn)偽隨機(jī)編碼源激發(fā)下的電磁信號合成.首先簡述偽隨機(jī)源激發(fā)下的電磁信號合成基本理論,然后根據(jù)合成公式中的大地脈沖響應(yīng)、偽隨機(jī)編碼源、儀器帶限影響及噪聲逐一進(jìn)行研究,最后給出信號合成結(jié)果,并與一個野外實測信號作對比,來驗證信號合成的正確性.

    2 基本理論

    偽隨機(jī)源激發(fā)下的電磁信號合成工作是指給定一個m序列偽隨機(jī)編碼源,在編碼源的激發(fā)下,模擬大地及發(fā)射、接收系統(tǒng)和噪聲的共同響應(yīng),繼而合成一個仿真的野外觀測信號的過程.

    在野外工作中,發(fā)射機(jī)通過兩個接地電極向地下供入偽隨機(jī)交變電流信號,在離發(fā)射電極一定距離處用不極化電極排列接收相鄰電極間的電位差,根據(jù)這種發(fā)射接收裝置,發(fā)射機(jī)發(fā)射的電磁信號在大地中傳播,然后通過電極與導(dǎo)線、接收系統(tǒng)接收.在勘探工作中,大地是探測目標(biāo),經(jīng)過大地的信號為有用信號.對有用信號產(chǎn)生影響的因素主要為噪聲,除此之外,電極附近的極化效應(yīng)、接收器的頻率響應(yīng)限制(即帶限)、電線的感抗等都對有用信號或多或少產(chǎn)生影響.在這些影響因素中,除噪聲外,每一部分的影響都相當(dāng)于在有用信號上加載了一個濾波器,最后接收到的信號是對有用信號經(jīng)過多次濾波后的綜合信號.因此,模擬電磁信號工作應(yīng)當(dāng)考慮各種因素的影響,然而,現(xiàn)實中有些因素的影響不易模擬或相對較小,本文只考慮噪聲及儀器帶限影響,并且將發(fā)射和接收系統(tǒng)的帶限影響綜合考慮.

    接收信號的合成公式可以用(1)式表示(Ziolkowski et al., 2007):

    v(t)=i(t)*g(t)*r(t)+n(t),

    (1)

    其中,t表示接收時間,v(t)表示接收機(jī)在接收點接收的信號,i(t)表示發(fā)射系統(tǒng)向大地發(fā)射的電流,g(t)表示經(jīng)由大地濾波后在接收點處的大地脈沖響應(yīng),r(t)表示儀器頻率響應(yīng).n(t)為接收點處的噪聲.

    從公式(1)可以看出,只要獲得了模型地電結(jié)構(gòu)的大地脈沖響應(yīng)g(t)、偽隨機(jī)發(fā)射源信號i(t)、儀器帶限響應(yīng)r(t)及噪聲n(t),就可合成出仿真的電磁信號v(t).下面對公式(1)中各項的實現(xiàn)過程逐一闡述.

    3 大地脈沖響應(yīng)

    理想情況下,向地下輸入脈沖電源,在不計及各種影響因素時,接收機(jī)測量的輸出即大地脈沖響應(yīng).但在野外工作中,一般采用方波作為輸入.在數(shù)值模擬時,先在頻率域獲得不同頻率的電磁響應(yīng),然后通過余弦變換轉(zhuǎn)換到時間域,得到大地階躍響應(yīng),再通過計算階躍響應(yīng)的時間導(dǎo)數(shù)來得到大地脈沖響應(yīng).

    3.1 頻率域響應(yīng)

    本文用均勻半空間地電模型來模擬地表電場信號的合成過程.根據(jù)文獻(xiàn)(Wright,2004),只有與源同方向的電場可提供有用信息,其他方向的場不能提供更多的信息,因此本文只用與源同方向的電場來合成所需信號.假設(shè)源的方向為x方向,則頻率域中同方向的電場可寫為(樸化榮,1990):

    [3cos2φ-2+(1+k1r)e-k1r],

    (2)

    用于頻率域到時間域變換的方法較多,如正弦變換或余弦變換(Anderson,1982;Newman et al., 1986; Li and Constable,2010)、折線法(Li et al., 2011),這三種方法同屬于傅里葉逆變換.在這三種方法中,正弦變換與余弦變換需要計算較寬頻帶內(nèi)的頻率響應(yīng),其頻帶寬度與時間點有關(guān);折線法相對來說計算速度快,但相對粗糙.除此之外,還有基于拉普拉斯逆變換的系列算法(Li and Farquharson,2015)可以實現(xiàn)信號的頻時域轉(zhuǎn)換,這類算法的優(yōu)勢在于可以用較少的濾波系數(shù)得到較精確的結(jié)果,但其對計算機(jī)精度的要求太高(4倍精度).本文最終選擇了余弦變換作為把信號從頻率域轉(zhuǎn)換到時間域的方法.

    余弦變換的表達(dá)式為(樸化榮,1990):

    (3)

    其中,f(t)為時間域的場,F(xiàn)(ω)為頻率域的場.

    參考Anderson (1982)的方法,選用函數(shù)對

    (4)

    來求余弦變換的濾波權(quán)系數(shù),其中,a為大于0的實數(shù).利用濾波系數(shù)可將(3)式寫為:

    }/t,

    (5)

    其中Ci為余弦變換的濾波權(quán)系數(shù),Ai-x為移動的橫坐標(biāo),N1和N2的大小由積分時間自動調(diào)節(jié),本文所用的余弦變換子程序為公開的代碼(方文藻等,1993).

    3.3 大地脈沖響應(yīng)

    殷長春等(2013)利用階躍響應(yīng)與電流的時間導(dǎo)數(shù)褶積代替脈沖響應(yīng)與電流褶積的方法成功實現(xiàn)了任意波形全時響應(yīng)的計算,避免了記錄時間t→0時大地脈沖奇異值的出現(xiàn),但只考慮了電流變化感應(yīng)出的電磁場,并沒有考慮電流自身產(chǎn)生的直流場,所以還需再將直流場加進(jìn)感應(yīng)的場中.本文采用Strack(1992)的提法,直接計算階躍響應(yīng)的時間導(dǎo)數(shù)來得到大地脈沖響應(yīng).t→0處的值用第一個階躍值與該處時間的比值來代替,通過將起始時間設(shè)計得較小(如0.01 ms)來減小t→0時的誤差.

    設(shè)計一個電阻率是100 Ωm的均勻半空間模型,假設(shè)源的方向為x方向,源的長度為200 m,計算點位于源的延長線上且離源中心點500 m處,電場的頻率域響應(yīng)曲線如圖1a所示.圖1a顯示了從10-7Hz到107Hz的計算結(jié)果,這個頻帶寬度可滿足1 s內(nèi)所有時間點上進(jìn)行余弦變換時所需的頻率范圍.從圖1a可以看出,在高頻和在低頻一定范圍內(nèi),電場曲線保持水平,不受頻率的影響.通過余弦變換將頻率域結(jié)果轉(zhuǎn)換到時間域,得到如圖1b所示的階躍響應(yīng),將其與解析公式計算結(jié)果對比,發(fā)現(xiàn)二者的吻合程度非常高,驗證了余弦變換過程的正確性.采用3.3節(jié)中提到的階躍響應(yīng)的時間導(dǎo)數(shù)來計算脈沖響應(yīng),計算結(jié)果如圖1c所示,對圖中所示的x方向電場來說,曲線中出現(xiàn)了一個峰值,這個峰值對后續(xù)的處理非常重要,因不是本文的研究內(nèi)容,所以不進(jìn)行闡述.這里為了突出峰值的形態(tài),對曲線進(jìn)行了截斷.所用的時間域解析公式來自文獻(xiàn)(樸化榮,1990).

    圖1 均勻半空間模型的頻率域與時間域響應(yīng)(a) 頻率域大地響應(yīng)曲線; (b) 時間域大地階躍響應(yīng); (c) 時間域大地脈沖響應(yīng).Fig.1 Responses of homogenous half space in frequency domain and time domain(a) Response curve in frequency domain; (b) Step response of earth in time domain; (c) Impulse response of earth in time domain.

    4 m序列偽隨機(jī)源

    4.1 m序列偽隨機(jī)編碼

    偽隨機(jī)序列是由移位寄存器產(chǎn)生的, 根據(jù)文獻(xiàn)(林可祥和汪一飛,1978),可通過查表的方式得到不同級移位寄存器產(chǎn)生的偽隨機(jī)序列的本原多項式F(x).通過F(x)與序列多項式G(x)的關(guān)系式G(x)=1/F(x),用多項式的長除算法可計算出相應(yīng)的序列多項式G(x),從G(x)中可直接獲得m序列偽隨機(jī)編碼.

    序列多項式G(x)最終可寫為下面的表達(dá)形式:

    (6)

    其中⊕表示模二加法運算符,x為移位運算符,其上標(biāo)的數(shù)字表示移位的位數(shù),例如:akx2表示將ak位移兩位,an所組成的序列{an}=a0,a1,a2,a3,…即為m序列偽隨機(jī)編碼.(5)式中的N為偽隨機(jī)碼一個周期的碼元個數(shù),N=2r-1,r為移位寄存器級數(shù).

    圖2是4級和6級的m序列偽隨機(jī)編碼一個周期的波形,它是通過文獻(xiàn)(林可祥和汪一飛,1978)本原多項式列表中4級和6級第一個本原多項式得到的m序列的編碼波形,從圖中看出其隨機(jī)性隨著移位寄存器級數(shù)的增多而愈發(fā)明顯.

    圖2 m序列偽隨機(jī)編碼一個周期的波形(a) 4級移位寄存器產(chǎn)生的波形; (b) 6級移位寄存器產(chǎn)生的波形.Fig.2 Waveforms of pseudo random binary sequence(a) Generated by 4-rank shift register; (b) Generated by 6-rank shift register.

    4.2 m序列偽隨機(jī)編碼源的頻譜

    對m序列偽隨機(jī)編碼應(yīng)用快速傅氏變換,獲得了其作為源時相應(yīng)的頻譜.通過對比不同級數(shù)的頻譜幅值,可了解不同級數(shù)的偽隨機(jī)信號所包含的頻率成分及幅頻特性.本文給出了4、6、8、10、12級偽隨機(jī)編碼源的頻譜,其中與偽隨機(jī)編碼相關(guān)的參數(shù)見表1.為了對比不同源譜的寬度,將所有偽隨機(jī)編碼周期與采樣率都定為0.25 s及16384 Hz.偽隨機(jī)編碼源的頻譜如圖3所示.

    表1 偽隨機(jī)編碼相關(guān)參數(shù)

    從圖3可以看出,不同級移位寄存器得到的頻譜形態(tài)不同,級數(shù)越高,主瓣的寬度越寬,當(dāng)移位寄存器的級數(shù)r=12時,從0 Hz到8192 Hz范圍內(nèi)的頻譜豐富,包絡(luò)線的形態(tài)接近水平,接近于脈沖源的頻譜,由此源激發(fā)得到的大地響應(yīng)會接近于脈沖源

    圖3 不同階m序列偽隨機(jī)編碼源發(fā)射波形的頻譜特性(a) 4級移位寄存器; (b) 6級移位寄存器; (c) 8級移位寄存器; (d) 10級移位寄存器; (e) 12級移位寄存器.Fig.3 Spectral properties of pseudo-random binary sequence generated by shift registers of different ranks(a) 4-rank; (b) 6-rank; (c) 8-rank; (d) 10-rank; (e) 12-rank.

    激發(fā)得到的響應(yīng).圖3中從10級到12級偽隨機(jī)編碼源的頻譜外包絡(luò)線不光滑,這是采樣率不夠?qū)е碌?若提高采樣率,則能得到外包絡(luò)線較光滑、幅值較均勻的譜圖.從圖3還可以看出,隨著寄存器級數(shù)的增加,主瓣的幅值也隨之下降,說明信號的能量降低.對地球物理勘探來說,信號能量降低會導(dǎo)致電磁信號的穿透深度變淺.換句話說,雖然主瓣頻率成分提高,分辨率也提高,但信號強(qiáng)度因能量分散而降低,導(dǎo)致探測深度淺.若想解決這個問題,需在加大發(fā)射功率的同時增加垂直疊加的時間,而且在后續(xù)處理采用類地震水平疊加的技術(shù).

    5 其他影響因素

    5.1 儀器頻帶限制的影響

    發(fā)射機(jī)和接收器自身的頻帶寬度是有限的,作為理論研究,可以用有一定帶寬的低通濾波器來模擬儀器的頻帶寬度(即帶限問題).

    低通濾波器的設(shè)計方法有多種,不同方法得到的濾波器的旁瓣幅值大小以及從主瓣到旁瓣的過渡帶寬度不同.理想濾波器要求過渡帶盡可能窄,旁瓣幅值盡可能小,以減小能量的泄露,使能量盡可能多地保留在主瓣內(nèi).本文選用海明窗函數(shù)法作為濾波器的設(shè)計方法,因為海明窗可將99.963%的能量集中在窗譜的主瓣內(nèi),旁瓣的峰值小于主瓣峰值的1%(程佩青,2007),符合理想濾波器的要求.

    海明窗函數(shù)為:

    (7)

    借鑒頻率域電磁儀器的頻帶寬度,設(shè)計一個低通濾波器,其通帶為8192Hz,阻帶為9600Hz.理想的低通濾波器的頻帶如圖4a所示,在通帶內(nèi),信號可以通過,通帶外是阻帶,阻帶內(nèi)的信號全部被濾除.理想的低通濾波器在時間域內(nèi)如圖4c中的黑實線所示(實際為無限信號,為了和加窗后的信號相比,只顯示了其中的一部分).圖4b是海明窗在頻率域的表現(xiàn)形式,用海明窗對圖4c中的信號進(jìn)行截斷,截斷后的信號如圖4c中的離散點所示,將加窗后的時間域濾波器變換到頻率域,如圖4d所示,相對于理想低通濾波器的頻帶(圖4a),圖4d出現(xiàn)了通帶與阻帶之間的過渡帶,但較窄,阻帶出現(xiàn)了小幅

    圖4 低通濾波器的特性(a) 理想的低通濾波器的譜; (b) 海明窗; (c) 理想的低通濾波器與經(jīng)過海明窗截斷后的濾波器對比; (d) 加窗后的信號對應(yīng)的低通濾波器的譜.Fig.4 Properties of low-pass filter(a) Ideal spectrum; (b) Hamming window; (c) Comparison of ideal low-pass filter and filter cut-off by Hamming window; (d) Spectrum of low-pass filter after cutting.

    振蕩(因其相對于主瓣來說,能量較小,振蕩不明顯,若放大該部分,可以看到振蕩現(xiàn)象),說明所設(shè)計的濾波器過渡帶較窄,阻帶能量泄漏很小,是一個性能優(yōu)良的濾波器.

    5.2 噪聲的加入

    當(dāng)用電磁法在野外工作時,常遇到的噪聲類型有兩種,一是白噪聲,二是50 Hz的工頻噪聲.

    (1) 白噪聲

    本文用正態(tài)分布隨機(jī)數(shù)來模擬噪聲:

    (8)

    其中,n1為所求的正態(tài)分布隨機(jī)數(shù),μ為正態(tài)分布的均值,本文為觀測資料,σ2為正態(tài)分布的方差,在本文中σ為觀測值與噪聲水平之積,Ni為0到1之間均勻分布的隨機(jī)數(shù),m為隨機(jī)數(shù)的個數(shù).

    (2) 50 Hz噪聲

    用一個頻率是50 Hz的正弦波來模擬工頻噪聲:

    n2=Asin(100πt),

    (9)

    其中n2為噪聲值,A為50 Hz工頻噪聲的振幅.

    6 偽隨機(jī)源激發(fā)下的電磁信號合成

    6.1 信號合成

    通過(1)式將以上內(nèi)容進(jìn)行組合,便可分別得到理想接收信號、加帶限的信號和帶噪聲的信號.

    以4階m序列偽隨機(jī)編碼源作為激發(fā)源(見圖5a),將源和均勻半空間大地脈沖響應(yīng)(見圖1c)做褶積,得到電場的理想信號,如圖5b所示.將理想信號與低通濾波器(儀器帶限)做褶積,便可得到考慮了帶限影響的電場信號,如圖5c所示.從這些圖可以看出,不加帶限前,信號比較平滑,加了帶限后,信號出現(xiàn)了一些毛刺狀的干擾.

    對加了帶限的信號分別加上了5%及10%的白噪聲,如圖6a和6b所示,噪聲疊加在原來的信號之上,引起信號的畸變,當(dāng)噪聲達(dá)到10%時,形成的干擾已非常明顯.圖6c為幅值是4倍于有用信號的50 Hz干擾信號,干擾信號和有用信號疊加后的信號如圖6d所示,可以看出,有用信號湮沒在了50 Hz干擾信號中.

    6.2 合成信號與實測信號的對比

    中國科學(xué)院項目組2015年在河北張北地區(qū)進(jìn)行了野外試驗.用一個長度是300 m的接地電偶極發(fā)射一個12階的偽隨機(jī)編碼源,采用軸向裝置,用一系列極距是60 m的電偶極子一字排開,在離發(fā)射源不同距離處的接收端接收耦合了大地響應(yīng)及各種影響的信號.為了和野外信號對比,在室內(nèi)也進(jìn)行了模擬.根據(jù)測區(qū)已知資料,在室內(nèi)用電阻率是100 Ωm的均勻半空間模型來近似地下介質(zhì),發(fā)射源也用12級移位寄存器生成的m序列偽隨機(jī)編碼序列,編碼的部分形態(tài)如圖7a所示.采用了和野外一致的碼元長度及采樣率(碼元長度為1/1024 s,采樣率為16 kHz,發(fā)射一個周期的信號大約需要4 s),因信號太長,為了突出細(xì)節(jié),我們僅顯示了其中的一小部分.值得注意的是,室內(nèi)的編碼序列和室外的編碼序列是不一致的,因為對于任何級的移位寄存器來說,可以產(chǎn)生很多種形態(tài)編碼,級數(shù)越多,能產(chǎn)生的編碼形態(tài)也越多,所以并未刻意和野外用相同的編碼形式.模擬時,分別在偏移距由近及遠(yuǎn)的位置接收信號,300 m和2700 m處的脈沖響應(yīng)分別如圖7b和圖7c所示,離源較近的接收點脈沖響應(yīng)的信號很強(qiáng),并且和離源較遠(yuǎn)時接收到的脈沖響應(yīng)形態(tài)有差異,這是因為離源近,信號接近飽和,離源越遠(yuǎn),得到的脈沖響應(yīng)越平緩.

    圖5 在4階偽隨機(jī)源激發(fā)下的接收信號(a) 4階偽隨機(jī)源; (b) 不考慮帶限的電場理想信號; (c) 考慮帶限的電場信號.Fig.5 Signal excited by 4-rank pseudo-random binary sequence(a) Source; (b) Ideal signal; (c) Signal considering band limit.

    圖6 在4階偽隨機(jī)源激發(fā)下考慮噪聲的信號(a) 加5%白噪聲; (b) 加10%白噪聲; (c) 50 Hz干擾噪聲; (d) 加50 Hz干擾噪聲.Fig.6 Signal with different noise type(a) Adding 5% white noise; (b) Adding 10% white noise; (c) 50 Hz industry noise; (d) Adding 50 Hz industrial noise.

    圖7 模擬信號和野外實測信號的對比(a) 12級m序列偽隨機(jī)信號部分波形; (b) 偏移距是300 m時的電場脈沖響應(yīng); (c) 偏移距是2700 m時的電場脈沖響應(yīng); (d) 偏移距是300 m時的電場計算信號; (e) 偏移距是300 m時的電場實測信號.Fig.7 Comparison of simulated signal to survey signal(a) Source; (b) Impulse response at offset 300 m; (c) Impulse response at offset 2700 m; (d) Simulated signal at offset 300 m; (e) Measured signal at offset 300 m.

    由于本文所依據(jù)的理論是電偶極子理論,當(dāng)發(fā)射偶極與接收偶極子距離較近時,對發(fā)射偶極源進(jìn)行了離散疊加.對300 m的收發(fā)距,將偶極源劃分成一系列5 m長的小偶極子,離散后的小偶極子是300 m收發(fā)距的1/60,滿足偶極子理論,同時大大減小了線間的感應(yīng)耦合.

    我們計算了部分點的接收信號,并和野外接收信號做了對比,發(fā)現(xiàn)二者的相似程度較高,這里只展示偏移距是300 m時的合成信號,如圖7c所示,通過和野外300 m處接收信號的屏幕截圖7e對比,可以看出信號的相似度較高.二者不盡相同的原因有三個:一是所用的編碼形式不同,二是地下介質(zhì)的電阻率不同,這兩個在前面也有所描述,還有一個原因,就是在野外實測中所受的干擾更多,如電極附近的極化效應(yīng),接收電線的感抗,而在數(shù)值模擬中只是模擬了較為典型的干擾.即便如此,二者較高的相似度說明了本文信號合成過程是正確的.

    7 結(jié)論

    本文根據(jù)偽隨機(jī)編碼源信號的傳播物理過程對接收信號進(jìn)行了合成,通過研究,發(fā)現(xiàn)大地脈沖響應(yīng)離源較近時會出現(xiàn)飽和現(xiàn)象,是一條單調(diào)衰減曲線,隨著測點遠(yuǎn)離源,峰值逐漸出現(xiàn).

    本文通過計算階躍響應(yīng)的時間導(dǎo)數(shù)來得到大地脈沖響應(yīng),t→0處的值用第一個階躍值與該處時間的比值來代替,通過將起始時間設(shè)計得較小(如0.01 ms)來減小t→0時的誤差,通過數(shù)值模擬發(fā)現(xiàn),這種近似方法可行.

    當(dāng)用移位寄存器得到m序列偽隨機(jī)編碼時,移位寄存器級數(shù)越多,其頻譜主瓣的寬度越寬,當(dāng)移位寄存器的級數(shù)r=12時,從0~8192 Hz范圍內(nèi)的各個頻率點上的幅值外包絡(luò)線形態(tài)接近水平,接近于脈沖響應(yīng)的頻譜,由此源激發(fā)而得到大地響應(yīng)會接近于脈沖源激發(fā)得到的響應(yīng),從而有較高的電性結(jié)構(gòu)分辨能力.然而隨著寄存器級數(shù)的增加,主瓣的各個頻率信號的幅值也隨之下降,影響了勘探深度的提高,為此需通過增加源的發(fā)射時間(便于加大垂直疊加次數(shù))和增加擬地震水平疊加次數(shù)來彌補(bǔ)因增加頻率成分而導(dǎo)致的信號強(qiáng)度損失,以達(dá)到既提高分辨率,又提高探測深度的雙重目標(biāo).

    由本文的方法得到的大地脈沖響應(yīng)與理論值吻合,合成信號與野外信號非常接近,因此,本文方法所產(chǎn)生的中間結(jié)果及理論信號合理可靠,可用在后續(xù)的數(shù)據(jù)處理中,并可用于檢驗后續(xù)數(shù)據(jù)處理方法的可靠性和有效性.

    Anderson W L. 1982. Calculation of Transient Soundings for a Coincident Loop System. Denver, Colo.: U. S. Geological Survey.Cheng P Q. 2007. Digital Signal Processor Tutorial. 3rd ed. (in Chinese). Beijing: Tsinghua University Press.

    Di Q Y, Wang G J, An Z G, et al. 2006. Geophysical exploration of a long deep tunnel on the west route of South to North Water Diversion Project.ChineseJ.Geophys. (in Chinese), 49(6): 1836-1842.Fang W Z, Li Y G, Li X. 1993. Theory of Transient Electromagnetic Sounding (in Chinese). Xi′an: Northwestern Polytechnical University Press. He J S. 2010. Full Domain Electromagnetic Method and Pseudo-Random Signal Electric Method (in Chinese). Beijing: China Higher Education Press.

    Li J H, Farquharson C G. 2015. Two effective inverse Laplace transform algorithms for computing time-domain electromagnetic responses.∥ SEG Technical Program Expanded Abstracts. SEG: 957-962.Li J H, Zhu Z Q, Liu S C, et al. 2011. 3D numerical simulation for the transient electromagnetic field excited by the central loop based on the vector finite-element method.JournalofGeophysicsandEngineering, 8(4): 560-567.Li Y G, Constable S. 2010. Transient electromagnetic in shallow water: insights from 1D modeling.ChineseJ.Geophys., 53(3): 737-742, doi: 10.3969/j.issn.0001-5733.2010.03.029.

    Lin K X, Wang Y F. 1978. Theory and Application of Pseudo-Random Code (in Chinese). Beijing: Posts and Telecom Press.

    Luo X Z, Li D W, Peng F P, et al. 2014. Implementation and applications of an coded electrical instrument with anti-interference ability.ProgressinGeophysics(in Chinese), 29(2): 944-951, doi: 10.6038/pg20140263.

    Luo Y Z, Lu Z G, Sun G L, et al. 2015. New generation of instruments for electrical and electromagnetic prospecting-instruments using pseudo random signal.ProgressinGeophysics(in Chinese), 30(1): 411-415, doi: 10.6038/pg20150159.Newman G A, Hohmann G W, Anderson W L. 1986. Transient electromagnetic response of a three-dimensional body in a layered earth.Geophysics, 51(8): 1608-1627.

    Piao H R. 1990. Theory of Electromagnetic Sounding (in Chinese). Beijing: Geological Publishing House.

    Qi Y F, Yin C C, Wang R, et al. 2015. Multi-transient EM full-time forward modeling and inversion of m-sequences.ChineseJ.Geophys. (in Chinese), 58(7): 2566-2577, doi: 10.6038/cjg20150731.Strack K M. 1992. Exploration with Deep Transient Electromagnetics. Netherlands: Elsevier Science Publishers.Tang J T, Luo W B. 2008. Pseudo-random electromagnetic exploration based on invert-repeated m-Sequence correlation identification.ChineseJ.Geophys. (in Chinese), 51(4): 1226-1233.

    Wright D, Ziolkowski A, Hobbs B A. 2002. Hydrocarbon detection and monitoring with a multicomponent transient electromagnetic (MTEM) survey.TheLeadingEdge, 21(9): 852-864. Wright D A. 2004. Detection of hydrocarbons and their movement in a reservoir using time-lapse multichannel transient electromagnetic (MTEM) data [Ph. D. thesis]. Edinburgh: University of Edinburgh.Wu X, Xue G Q, Di Q Y, et al. 2015. Accurate identification for the electromagnetic impulse response of the earth with pseudo random coded waveforms.ChineseJ.Geophys. (in Chinese), 58(8): 2792-2802, doi: 10.6038/cjg20150814.

    Xue G Q, Gelius L J, Sakyi P A, et al. 2014. Discovery of a hidden BIF deposit in Anhui province, China by integrated geological and geophysical investigations.OreGeologyReview, 63: 470-477.

    Xue G Q, Yan S, Di Q Y, et al. 2015. Technical analysis of multi-transient electromagnetic method.JournalofEarthSciencesandEnvironment(in Chinese), 37(1): 94-100.

    Yin C C, Huang W, Ben F. 2013. The full-time electromagnetic modeling for time-domain airborne electromagnetic systems.ChineseJ.Geophys. (in Chinese), 56(9): 3153-3162, doi: 10.6038/cjg20130928.

    Zhao B R, Zhao J, Zhang H K, et al. 2006. The PS100 high precision earth electricity instrument system (IP to IP) with controllable source-Application of CDMA technology to the measurement of earth resistivity for the first time.ProgressinGeophysics(in Chinese), 21(2): 675-682.

    Ziolkowski A, Hobbs B A, Wright D. 2007. Multitransient electromagnetic demonstration survey in France.Geophysics, 72(4): F197-F209.

    附中文參考文獻(xiàn)

    程佩青. 2007. 數(shù)字信號處理教程(第三版). 北京: 清華大學(xué)出版社.

    底青云, 王光杰, 安志國等. 2006. 南水北調(diào)西線千米深長隧洞圍巖構(gòu)造地球物理勘探. 地球物理學(xué)報, 49(6): 1836-1842.

    方文藻, 李予國, 李貅. 1993. 瞬變電磁測深法原理. 西安: 西北工業(yè)大學(xué)出版社.

    何繼善. 2010. 廣域電磁法和偽隨機(jī)信號電法. 北京: 高等教育出版社.

    林可祥, 汪一飛. 1978. 偽隨機(jī)碼的原理與應(yīng)用. 北京: 人民郵電出版社.

    羅先中, 李達(dá)為, 彭芳蘋等. 2014. 抗干擾編碼電法儀的實現(xiàn)及應(yīng)用. 地球物理學(xué)進(jìn)展, 29(2): 944-951, doi: 10.6038/pg20140263.

    羅延鐘, 陸占國, 孫國良等. 2015. 新一代電法勘查儀器——偽隨機(jī)信號電法儀. 地球物理學(xué)進(jìn)展, 30(1): 411-415, doi: 10.6038/pg20150159.

    樸化榮. 1990. 電磁測深法原理. 北京: 地質(zhì)出版社.

    齊彥福, 殷長春, 王若等. 2015. 多通道瞬變電磁m序列全時正演模擬與反演. 地球物理學(xué)報, 58(7): 2566-2577, doi: 10.6038/cjg20150731.

    湯井田, 羅維斌. 2008. 基于相關(guān)辨識的逆重復(fù)m序列偽隨機(jī)電磁法. 地球物理學(xué)報, 51(4): 1226-1233.

    武欣, 薛國強(qiáng), 底青云等. 2015. 偽隨機(jī)編碼源電磁響應(yīng)的精細(xì)辨識. 地球物理學(xué)報, 58(8): 2792-2802, doi: 10.6038/cjg20150814.

    薛國強(qiáng), 閆述, 底青云等. 2015. 多道瞬變電磁法(MTEM)技術(shù)分析. 地球科學(xué)與環(huán)境學(xué)報, 37(1): 94-100.

    殷長春, 黃威, 賁放. 2013. 時間域航空電磁系統(tǒng)瞬變?nèi)珪r響應(yīng)正演模擬. 地球物理學(xué)報, 56(9): 3153-3162, doi: 10.6038/cjg20130928.

    趙璧如, 趙健, 張洪魁等. 2006. PS100型IP到端可控源高精度大地電測儀系統(tǒng)——CDMA技術(shù)首次在地電阻率測量中的應(yīng)用. 地球物理學(xué)進(jìn)展, 21(2): 675-682.

    (本文編輯 何燕)

    Signal synthesis of time-domain electromagnetic response excited by a pseudo random binary sequence source

    WANG Ruo1, WANG Miao-Yue1, DI Qing-Yun1, XUE Guo-Qiang1, YIN Chang-Chun2, LEI Da1

    1KeyLaboratoryofShaleGasandEngineeringGeology,InstituteofGeologyandGeophysics,ChineseAcademyofSciences,Beijing100029,China2CollegeofGeo-explorationScienceandTechnology,JilinUniversity,Changchun130026,China

    After the pseudo-random binary sequence source is introduced into the artificial electromagnetic method, the exploration depth and the resolution are improved, because in this approach correlation technique can be used in the data processing to suppress noise. Most researchers didn’t focus on the synthetic signal rather than data processing. However, some data processing techniques have to use the electromagnetic signal, therefore it is indispensable to rely on synthesis of electromagnetic signal. In this paper, the electromagnetic signal is synthesized according to the propagating routine of the electromagnetic field. Firstly, electromagnetic response of a simulated geo-electric model in the frequency domain is calculated according to analytic formula. Secondly, the step electromagnetic response in the time domain is obtained by applying cosine transform to the response in the frequency domain, and then impulse response of the earth is obtained from time derivative of the step response. Through convoluting the impulse response of the earth and with PRBS source, the ideal electromagnetic signal can be composed. Finally, a low-pass filter is used to simulate the synthetic system response of the transmitter and receiver because all of the instruments are band-limited. Then the simulated signal is obtained by convoluting a low pass filter to the ideal signal added with noise. Comparing with actual field received signal proves that the composed signal is correct, and the signal of high simulation can be used to serve the data processing technique.

    Pseudo randomly coded technology; Time domain; Impulse response; Received signal; Compose

    10.6038/cjg20161203.

    國家重大科研裝備研制項目“深部資源探測核心裝備研發(fā)”(ZDYZ2012-1-05-04),國家重點研發(fā)計劃(2016YFC0600507),國家自然科學(xué)基金(41174111)資助.

    王若,女,博士,主要從事電磁法正反演研究.E-mail:wangruo@mail.iggcas.ac.cn

    10.6038/cjg20161203

    P631

    2016-02-01,2016-11-14收修定稿

    王若, 王妙月, 底青云等. 2016. 偽隨機(jī)編碼源激發(fā)下的時域電磁信號合成. 地球物理學(xué)報,59(12):4414-4423,

    Wang R, Wang M Y, Di Q Y, et al. 2016. Signal synthesis of time-domain electromagnetic response excited by a pseudo random binary sequence source.ChineseJ.Geophys. (in Chinese),59(12):4414-4423,doi:10.6038/cjg20161203.

    猜你喜歡
    信號
    信號
    鴨綠江(2021年35期)2021-04-19 12:24:18
    完形填空二則
    7個信號,警惕寶寶要感冒
    媽媽寶寶(2019年10期)2019-10-26 02:45:34
    孩子停止長個的信號
    《鐵道通信信號》訂閱單
    基于FPGA的多功能信號發(fā)生器的設(shè)計
    電子制作(2018年11期)2018-08-04 03:25:42
    基于Arduino的聯(lián)鎖信號控制接口研究
    《鐵道通信信號》訂閱單
    基于LabVIEW的力加載信號采集與PID控制
    Kisspeptin/GPR54信號通路促使性早熟形成的作用觀察
    精品人妻视频免费看| 日韩欧美免费精品| 日本色播在线视频| 国内揄拍国产精品人妻在线| 国产蜜桃级精品一区二区三区| 99热只有精品国产| 欧美xxxx黑人xx丫x性爽| 国产亚洲欧美98| 欧美一级a爱片免费观看看| 国产精品98久久久久久宅男小说| 97人妻精品一区二区三区麻豆| 亚洲国产欧美人成| 午夜福利视频1000在线观看| 日韩欧美在线二视频| 人妻制服诱惑在线中文字幕| 欧美日韩精品成人综合77777| 最近视频中文字幕2019在线8| 在现免费观看毛片| 免费在线观看日本一区| 99久久久亚洲精品蜜臀av| 国产精品,欧美在线| 免费搜索国产男女视频| 麻豆av噜噜一区二区三区| 国产淫片久久久久久久久| 国产伦在线观看视频一区| 精品一区二区三区av网在线观看| 精品人妻熟女av久视频| 无遮挡黄片免费观看| 中文亚洲av片在线观看爽| 一卡2卡三卡四卡精品乱码亚洲| 嫁个100分男人电影在线观看| 午夜老司机福利剧场| 亚洲av免费在线观看| 欧美日韩黄片免| 最近最新中文字幕大全电影3| 干丝袜人妻中文字幕| 欧美一区二区亚洲| 欧美日本视频| 熟女人妻精品中文字幕| 国产一区二区三区在线臀色熟女| 男女做爰动态图高潮gif福利片| 99视频精品全部免费 在线| 伦精品一区二区三区| 亚洲avbb在线观看| 亚洲成人久久爱视频| 亚洲最大成人手机在线| 免费搜索国产男女视频| 他把我摸到了高潮在线观看| 国产精品无大码| 大又大粗又爽又黄少妇毛片口| 午夜老司机福利剧场| 成人国产麻豆网| 波野结衣二区三区在线| 国产视频一区二区在线看| 91狼人影院| 亚洲,欧美,日韩| 熟妇人妻久久中文字幕3abv| 日日摸夜夜添夜夜添av毛片 | 美女黄网站色视频| 99九九线精品视频在线观看视频| 欧美xxxx黑人xx丫x性爽| 亚洲美女搞黄在线观看 | 亚洲精品成人久久久久久| 成年女人永久免费观看视频| 欧美色欧美亚洲另类二区| 18+在线观看网站| 亚洲精品一区av在线观看| a级毛片a级免费在线| 国产成人av教育| 麻豆一二三区av精品| 22中文网久久字幕| 超碰av人人做人人爽久久| 两人在一起打扑克的视频| 中文资源天堂在线| 一级a爱片免费观看的视频| 嫩草影院新地址| 国产精品伦人一区二区| 在线观看免费视频日本深夜| 亚洲熟妇熟女久久| 欧美3d第一页| 国产黄片美女视频| 精品午夜福利在线看| 国内精品宾馆在线| 亚洲中文日韩欧美视频| 精品久久久久久久久亚洲 | 成年女人看的毛片在线观看| 五月玫瑰六月丁香| 一级黄片播放器| 亚洲精华国产精华液的使用体验 | 成年女人看的毛片在线观看| 成人特级av手机在线观看| 狂野欧美白嫩少妇大欣赏| 中文在线观看免费www的网站| 日本成人三级电影网站| 亚洲av美国av| 亚洲四区av| 免费观看精品视频网站| 午夜影院日韩av| 午夜亚洲福利在线播放| 内地一区二区视频在线| 色哟哟哟哟哟哟| 变态另类丝袜制服| 成人二区视频| 女同久久另类99精品国产91| 国产一区二区激情短视频| 很黄的视频免费| 久久久久久久久久成人| 免费不卡的大黄色大毛片视频在线观看 | 国产精品久久久久久精品电影| 欧美xxxx黑人xx丫x性爽| 亚洲人成网站高清观看| 97碰自拍视频| 日本免费a在线| 午夜福利成人在线免费观看| 成人av一区二区三区在线看| 欧美性猛交╳xxx乱大交人| 欧美中文日本在线观看视频| 麻豆成人av在线观看| 88av欧美| 97超级碰碰碰精品色视频在线观看| 亚洲性久久影院| av国产免费在线观看| 亚洲黑人精品在线| 99热精品在线国产| 五月伊人婷婷丁香| 麻豆成人午夜福利视频| 欧美日韩瑟瑟在线播放| 波多野结衣高清作品| 能在线免费观看的黄片| 99热精品在线国产| 最新中文字幕久久久久| 国产成人aa在线观看| 在线观看舔阴道视频| 97人妻精品一区二区三区麻豆| 最近最新免费中文字幕在线| 国产精品野战在线观看| 3wmmmm亚洲av在线观看| 1024手机看黄色片| 亚洲最大成人中文| 亚洲av成人av| 91麻豆av在线| 色哟哟哟哟哟哟| 亚洲性夜色夜夜综合| 国产精品久久久久久久久免| 国产极品精品免费视频能看的| 国产一区二区在线av高清观看| 亚洲人成伊人成综合网2020| 九九在线视频观看精品| 哪里可以看免费的av片| 欧洲精品卡2卡3卡4卡5卡区| 亚洲专区中文字幕在线| 国产精品久久久久久av不卡| 在线国产一区二区在线| 久久99热这里只有精品18| 亚洲国产色片| 亚洲精品一卡2卡三卡4卡5卡| 97超视频在线观看视频| 欧美丝袜亚洲另类 | 噜噜噜噜噜久久久久久91| 亚洲精华国产精华液的使用体验 | 国产精品一区二区性色av| a级毛片免费高清观看在线播放| 亚洲精品456在线播放app | 美女 人体艺术 gogo| 国产免费男女视频| 深夜精品福利| 国产精品久久久久久久久免| 久久久久久九九精品二区国产| 99九九线精品视频在线观看视频| 成人美女网站在线观看视频| 婷婷色综合大香蕉| 在线免费十八禁| 在线播放无遮挡| 草草在线视频免费看| 男人和女人高潮做爰伦理| 中国美白少妇内射xxxbb| 亚洲精品一区av在线观看| 亚洲精华国产精华精| 日韩大尺度精品在线看网址| 久久精品国产自在天天线| 一本一本综合久久| 国产精品国产高清国产av| 亚洲国产欧洲综合997久久,| 蜜桃久久精品国产亚洲av| 成人特级黄色片久久久久久久| 男女边吃奶边做爰视频| 中文字幕人妻熟人妻熟丝袜美| 国产精品综合久久久久久久免费| 久久久久久久精品吃奶| 日韩欧美在线乱码| 久久久成人免费电影| 国产人妻一区二区三区在| 熟妇人妻久久中文字幕3abv| 亚洲一级一片aⅴ在线观看| 久久精品国产亚洲网站| 成人特级黄色片久久久久久久| 国产精品精品国产色婷婷| 亚洲不卡免费看| 欧美一区二区精品小视频在线| 天美传媒精品一区二区| 免费观看精品视频网站| 国产一区二区三区av在线 | 成人午夜高清在线视频| 欧美潮喷喷水| 久久久成人免费电影| 国产中年淑女户外野战色| 免费无遮挡裸体视频| 国产国拍精品亚洲av在线观看| 干丝袜人妻中文字幕| 欧美三级亚洲精品| 天天躁日日操中文字幕| 久久久久久久久久成人| 床上黄色一级片| 干丝袜人妻中文字幕| 国产午夜福利久久久久久| 免费不卡的大黄色大毛片视频在线观看 | 99久国产av精品| 伦理电影大哥的女人| 亚洲精品色激情综合| 美女 人体艺术 gogo| 久久久久国产精品人妻aⅴ院| 我的女老师完整版在线观看| 久久草成人影院| 一区二区三区激情视频| 午夜亚洲福利在线播放| 黄色欧美视频在线观看| 婷婷色综合大香蕉| 精品免费久久久久久久清纯| 欧美中文日本在线观看视频| 天堂av国产一区二区熟女人妻| 欧美激情久久久久久爽电影| 欧美日本亚洲视频在线播放| 久久精品综合一区二区三区| 真人一进一出gif抽搐免费| 真实男女啪啪啪动态图| av福利片在线观看| 精品一区二区免费观看| 午夜福利视频1000在线观看| 中文字幕人妻熟人妻熟丝袜美| av福利片在线观看| 日韩中字成人| 我要看日韩黄色一级片| 欧美色视频一区免费| 午夜久久久久精精品| 熟妇人妻久久中文字幕3abv| 又粗又爽又猛毛片免费看| 哪里可以看免费的av片| 一本一本综合久久| 精品一区二区三区av网在线观看| 亚洲国产高清在线一区二区三| 久久精品夜夜夜夜夜久久蜜豆| 天美传媒精品一区二区| 99久国产av精品| 久久精品国产亚洲av涩爱 | 国产欧美日韩一区二区精品| 窝窝影院91人妻| 三级国产精品欧美在线观看| 天堂动漫精品| 啪啪无遮挡十八禁网站| 热99re8久久精品国产| 国产色爽女视频免费观看| 亚洲经典国产精华液单| 国产真实伦视频高清在线观看 | 色av中文字幕| 日韩,欧美,国产一区二区三区 | 国产精品99久久久久久久久| 窝窝影院91人妻| 欧美绝顶高潮抽搐喷水| 日韩精品青青久久久久久| 国产三级中文精品| 免费看a级黄色片| 国内少妇人妻偷人精品xxx网站| 中文字幕av成人在线电影| av在线观看视频网站免费| 久久精品国产清高在天天线| 欧美bdsm另类| 午夜福利欧美成人| 黄片wwwwww| 精品人妻视频免费看| 春色校园在线视频观看| 69人妻影院| 午夜免费激情av| 两人在一起打扑克的视频| 精品国内亚洲2022精品成人| 亚洲av五月六月丁香网| 自拍偷自拍亚洲精品老妇| 在线观看美女被高潮喷水网站| 国产精品久久久久久久久免| 天天一区二区日本电影三级| 少妇被粗大猛烈的视频| 性色avwww在线观看| 日日摸夜夜添夜夜添av毛片 | 色综合婷婷激情| 免费电影在线观看免费观看| 99久国产av精品| 97碰自拍视频| 国产成人aa在线观看| 精品一区二区三区视频在线观看免费| 欧美成人一区二区免费高清观看| 老熟妇乱子伦视频在线观看| 午夜亚洲福利在线播放| 午夜免费男女啪啪视频观看 | 久久久久国内视频| 亚洲欧美日韩东京热| 最近视频中文字幕2019在线8| 欧美潮喷喷水| 精品久久国产蜜桃| 大又大粗又爽又黄少妇毛片口| 亚洲美女视频黄频| 一卡2卡三卡四卡精品乱码亚洲| 久久人人精品亚洲av| 在线免费十八禁| 女同久久另类99精品国产91| 日本 欧美在线| 搞女人的毛片| 日本a在线网址| 尾随美女入室| 免费看美女性在线毛片视频| 亚洲精华国产精华液的使用体验 | 久久久久久久久大av| 婷婷丁香在线五月| 美女大奶头视频| 可以在线观看的亚洲视频| 国产亚洲精品久久久久久毛片| 超碰av人人做人人爽久久| 在线国产一区二区在线| 日韩强制内射视频| 91在线观看av| 最近视频中文字幕2019在线8| a级一级毛片免费在线观看| 亚洲熟妇中文字幕五十中出| 成人美女网站在线观看视频| 久久精品影院6| av黄色大香蕉| 久久人人精品亚洲av| 99久久无色码亚洲精品果冻| 老司机福利观看| 国产精品一区二区免费欧美| 成人美女网站在线观看视频| 日日夜夜操网爽| 亚洲欧美日韩无卡精品| 一本一本综合久久| 日韩国内少妇激情av| 久久人人爽人人爽人人片va| 国产精品伦人一区二区| 欧美三级亚洲精品| 嫩草影院新地址| 欧美最黄视频在线播放免费| 婷婷色综合大香蕉| 成人性生交大片免费视频hd| 成人高潮视频无遮挡免费网站| 亚洲av免费在线观看| 欧美成人性av电影在线观看| av专区在线播放| 国产一区二区三区av在线 | 在线天堂最新版资源| 我要看日韩黄色一级片| 亚洲第一电影网av| 男女啪啪激烈高潮av片| videossex国产| 国产黄片美女视频| 日韩高清综合在线| 在线天堂最新版资源| 午夜免费激情av| 全区人妻精品视频| a在线观看视频网站| 熟妇人妻久久中文字幕3abv| 亚洲av第一区精品v没综合| 啦啦啦观看免费观看视频高清| 亚洲成人中文字幕在线播放| 一区二区三区免费毛片| 国产v大片淫在线免费观看| 亚洲三级黄色毛片| 日本免费一区二区三区高清不卡| 99九九线精品视频在线观看视频| 欧美日韩亚洲国产一区二区在线观看| 两个人的视频大全免费| 校园人妻丝袜中文字幕| 精品久久久久久久久久免费视频| 搡老岳熟女国产| 久99久视频精品免费| 欧美激情久久久久久爽电影| 欧美性感艳星| 两性午夜刺激爽爽歪歪视频在线观看| 女人十人毛片免费观看3o分钟| 免费看光身美女| 久久久精品欧美日韩精品| 91在线观看av| 一区二区三区四区激情视频 | 特大巨黑吊av在线直播| 国内揄拍国产精品人妻在线| 老师上课跳d突然被开到最大视频| 国产精品久久久久久精品电影| 精品久久久久久久久久久久久| 国产私拍福利视频在线观看| 日韩大尺度精品在线看网址| 男女边吃奶边做爰视频| 亚洲无线观看免费| 日本爱情动作片www.在线观看 | 成人三级黄色视频| 中文字幕精品亚洲无线码一区| 亚洲专区国产一区二区| 成人特级黄色片久久久久久久| 九九久久精品国产亚洲av麻豆| 人妻少妇偷人精品九色| 久久精品国产亚洲av涩爱 | 啦啦啦观看免费观看视频高清| 婷婷六月久久综合丁香| 麻豆精品久久久久久蜜桃| 国产精品不卡视频一区二区| 国产精品久久久久久亚洲av鲁大| 九九热线精品视视频播放| 性欧美人与动物交配| 成人国产综合亚洲| 亚洲最大成人av| 国内精品久久久久精免费| 97碰自拍视频| 日本a在线网址| 亚洲国产精品久久男人天堂| 91久久精品电影网| 色哟哟·www| aaaaa片日本免费| 国产精品乱码一区二三区的特点| 亚洲18禁久久av| 欧美国产日韩亚洲一区| 中文字幕av成人在线电影| 中国美女看黄片| 国产熟女欧美一区二区| 听说在线观看完整版免费高清| 精品日产1卡2卡| 亚洲人成网站在线播放欧美日韩| 亚洲av.av天堂| 97超级碰碰碰精品色视频在线观看| 久久久久久久久久成人| 美女大奶头视频| 三级毛片av免费| 国产伦一二天堂av在线观看| 国产午夜精品久久久久久一区二区三区 | av天堂中文字幕网| 成年女人永久免费观看视频| 久久久久久久亚洲中文字幕| 一a级毛片在线观看| 国产av一区在线观看免费| 99热只有精品国产| 色综合婷婷激情| 搡老妇女老女人老熟妇| 在线观看免费视频日本深夜| av女优亚洲男人天堂| 久久久久久九九精品二区国产| 亚洲欧美清纯卡通| 免费av不卡在线播放| 91在线精品国自产拍蜜月| 亚洲中文日韩欧美视频| 国产精品免费一区二区三区在线| 麻豆av噜噜一区二区三区| 欧美一区二区精品小视频在线| 久久国内精品自在自线图片| 欧美日韩综合久久久久久 | 国产白丝娇喘喷水9色精品| 中亚洲国语对白在线视频| 久久香蕉精品热| 午夜爱爱视频在线播放| 亚洲国产精品sss在线观看| 亚洲图色成人| 国产精品人妻久久久影院| av专区在线播放| 免费观看人在逋| 黄色日韩在线| 欧美色视频一区免费| 久久久午夜欧美精品| 特大巨黑吊av在线直播| 88av欧美| 中国美白少妇内射xxxbb| 一个人观看的视频www高清免费观看| 夜夜看夜夜爽夜夜摸| 午夜激情欧美在线| www日本黄色视频网| 精品人妻一区二区三区麻豆 | 亚洲无线在线观看| 欧美绝顶高潮抽搐喷水| 亚洲va在线va天堂va国产| 精品人妻偷拍中文字幕| 精品午夜福利视频在线观看一区| 亚洲av免费在线观看| 人人妻人人看人人澡| 天美传媒精品一区二区| 全区人妻精品视频| 日本撒尿小便嘘嘘汇集6| 一区福利在线观看| 99热精品在线国产| 国产高清有码在线观看视频| 熟女电影av网| aaaaa片日本免费| 成人国产麻豆网| 国产探花在线观看一区二区| 制服丝袜大香蕉在线| 变态另类成人亚洲欧美熟女| 国产成人aa在线观看| 美女免费视频网站| 日本一本二区三区精品| av在线天堂中文字幕| av天堂中文字幕网| 国产成人av教育| 男人舔女人下体高潮全视频| 啦啦啦观看免费观看视频高清| 免费观看人在逋| 女同久久另类99精品国产91| 少妇人妻精品综合一区二区 | 国产探花极品一区二区| 春色校园在线视频观看| 2021天堂中文幕一二区在线观| 欧美国产日韩亚洲一区| a级毛片免费高清观看在线播放| 成人高潮视频无遮挡免费网站| 亚洲欧美清纯卡通| 2021天堂中文幕一二区在线观| 色视频www国产| 国产一区二区三区av在线 | 色综合亚洲欧美另类图片| 欧洲精品卡2卡3卡4卡5卡区| 午夜a级毛片| 午夜福利成人在线免费观看| 最新中文字幕久久久久| 欧美另类亚洲清纯唯美| 亚洲av免费在线观看| 成人永久免费在线观看视频| 高清毛片免费观看视频网站| 97碰自拍视频| 直男gayav资源| 成年女人永久免费观看视频| 亚洲av五月六月丁香网| 久99久视频精品免费| 男女下面进入的视频免费午夜| 亚洲三级黄色毛片| 亚洲欧美日韩卡通动漫| 日本黄大片高清| 中出人妻视频一区二区| 亚洲av美国av| 国产免费男女视频| 好男人在线观看高清免费视频| 日韩欧美精品v在线| 日本一本二区三区精品| 免费观看人在逋| 999久久久精品免费观看国产| 欧美3d第一页| 十八禁网站免费在线| 天天一区二区日本电影三级| 亚洲四区av| 精华霜和精华液先用哪个| 国内精品久久久久精免费| 嫩草影院新地址| 国产精品98久久久久久宅男小说| 春色校园在线视频观看| 国产精品国产高清国产av| 欧美人与善性xxx| 男女之事视频高清在线观看| 欧美激情国产日韩精品一区| 国产黄色小视频在线观看| 老熟妇乱子伦视频在线观看| 18禁裸乳无遮挡免费网站照片| 91久久精品电影网| 干丝袜人妻中文字幕| 午夜a级毛片| 人妻丰满熟妇av一区二区三区| 亚洲性夜色夜夜综合| h日本视频在线播放| 真人一进一出gif抽搐免费| 日本五十路高清| 亚洲人成网站在线播| 久久久久久久久大av| 国产色婷婷99| 亚洲熟妇熟女久久| 国产精品人妻久久久久久| 国内毛片毛片毛片毛片毛片| 天天一区二区日本电影三级| 色播亚洲综合网| 999久久久精品免费观看国产| 可以在线观看的亚洲视频| ponron亚洲| 在线免费观看不下载黄p国产 | 最后的刺客免费高清国语| 我的老师免费观看完整版| 一进一出抽搐动态| 午夜福利18| 内地一区二区视频在线| 我的女老师完整版在线观看| 成人欧美大片| 国产伦人伦偷精品视频| 欧美一区二区精品小视频在线| 亚洲美女搞黄在线观看 | 99热这里只有是精品在线观看| 久久精品久久久久久噜噜老黄 | 亚洲va日本ⅴa欧美va伊人久久| 亚洲黑人精品在线| 永久网站在线| 午夜激情福利司机影院| 国产真实伦视频高清在线观看 | 久久久成人免费电影| 成人综合一区亚洲| 我要搜黄色片| a级毛片免费高清观看在线播放| 赤兔流量卡办理| 男女之事视频高清在线观看| 亚洲专区国产一区二区| 在现免费观看毛片| 男人舔奶头视频| 人妻久久中文字幕网| 精品乱码久久久久久99久播| 亚洲av中文字字幕乱码综合| 男女之事视频高清在线观看| 三级男女做爰猛烈吃奶摸视频| 久久九九热精品免费| 国产av在哪里看| 99热这里只有精品一区| 哪里可以看免费的av片|