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

    一種基于改進差分的測井數(shù)據(jù)可逆變長碼壓縮方法

    2018-07-28 07:19林永強唐露新劉海熊陵
    電腦知識與技術 2018年11期
    關鍵詞:基準值

    林永強 唐露新 劉?!⌒芰?/p>

    摘要:隨鉆測井系統(tǒng)需要提高速率和可靠性,常用的差分脈碼調制可去除數(shù)據(jù)時域冗余,但受干擾時會引發(fā)錯誤擴散。本文根據(jù)測井數(shù)據(jù)及信道傳輸特性,提出一種改進差分測井數(shù)據(jù)壓縮編碼方法,先將測井數(shù)據(jù)進行基準值的動態(tài)分段求差預處理,采用不等重保護策略,對處理的差值數(shù)據(jù)進行量化和可逆變長碼信源編碼、對基準值進行糾錯信道編碼,并采用非等重保護的幀結構和脈沖間隔調制進行傳輸。通過仿真實驗表明,數(shù)據(jù)壓縮可以達到48%,數(shù)據(jù)傳輸?shù)目煽啃悦黠@提升。

    關鍵詞:改進差分;基準值;可逆變長碼

    中圖分類號:TP3 文獻標志碼:A 文章編號:1009-3044(2018)11-0254-04

    A Kind of Inverse Variable Length Code Compression Method for Logging Data Based on Improved Difference

    LIU Yong-qiang1, TANG Lu-xin1, LIU Hai2, XIONG Ling2

    (1.Guangdong University of Technology, School of Electromechanical Engineering, Guangzhou 510006, China; 2.Zhongtian Qiming Petroleum Technology Co, Guangzhou 510535, China)

    Abstract: While logging-while-drilling systems require increased speed and reliability, commonly used differential pulse code modulation removes data time-domain redundancy but causes false proliferation when disturbed. In this paper, based on the logging data and channel transmission characteristics, an improved differential logging data compression coding method is proposed. First, the logging data is pre-processed by dynamic subsection difference of reference value. The unequal weight protection strategy is adopted to deal with the difference The value data is quantized and reversibly variable length code source code, the error correction channel code is used for the reference value, and the frame structure and the pulse interval modulation which are not equal weight are used for protection. Simulation results show that the data compression can reach 48%, and the reliability of data transmission is obviously improved.

    Keywords: improved difference; reference value; reversible variable length code

    1 引言

    目前最常用鉆井是隨鉆系統(tǒng),隨鉆測井系統(tǒng)需將數(shù)千米的井下數(shù)據(jù)傳輸?shù)骄?,常用?shù)據(jù)傳輸主要有鉆井液脈沖傳輸與電磁波數(shù)據(jù)傳輸,前者利用鉆井液壓力正/負脈沖傳輸信號比較穩(wěn)定可靠,國內(nèi)研發(fā)的鉆井液正脈沖發(fā)生器傳輸速率在0.5~3.0bit/s,國際主流的鉆井液正脈沖發(fā)生器傳輸速率一般在10bit/s[1]。目前油氣勘探開發(fā)向復雜地層發(fā)展,需要快速傳輸更多、更精確的測井數(shù)據(jù),但目前鉆井液脈沖傳輸難以滿足數(shù)據(jù)傳輸實時性和抗干擾的更高要求,因此提高測井數(shù)據(jù)傳輸性能成為隨鉆測量技術的一個重點。

    為提高隨鉆測量信息傳輸能力,發(fā)達國家已經(jīng)將數(shù)據(jù)壓縮技術、信息傳輸技術引入地質導向鉆井領域,取得較好效果,為了保密較少發(fā)表論文。斯倫貝謝公司利用圖像壓縮技術將井下井眼成像圖像壓縮50倍后傳輸至地面,很好地解決了信息傳輸能力不足的難題,保證了數(shù)據(jù)實時性[2]。國內(nèi)鉆井液脈沖傳輸?shù)臄?shù)據(jù)編碼研究開展較晚,2010年,張煜提出了基于差分脈碼調制的數(shù)據(jù)壓縮編碼算法,可將實時傳輸數(shù)據(jù)量壓縮為原始數(shù)據(jù)量的一半以下,將數(shù)據(jù)傳輸效率提高一倍以上,但信道噪聲導致數(shù)據(jù)錯誤時,會引起誤差擴散導致數(shù)據(jù)錯誤傳播[3]。2016年,李婷提出一種LZW壓縮編碼和RS糾錯碼聯(lián)合編碼的算法LZW-RS,實現(xiàn)壓縮比為2:1~3:1的無損數(shù)據(jù)壓縮,在信噪比為5.5dB時,LZW解碼能達到良好的重構性能,使數(shù)據(jù)傳輸效率最高提升2/3倍,誤碼率降低為LZW編碼的1/4倍~1/8倍,但編譯碼表是動態(tài)生成的,數(shù)據(jù)傳輸錯誤會導致譯碼錯誤[4]。同年,江佩佩提出一種聯(lián)合卷織編碼與脈沖間隔調制(PIM)的網(wǎng)格編碼脈沖間隔調制(TC-PIM)方法,在相同信息速率條件下,在誤碼率為10-5的時候,3比特TC-PIM與2比特PIM相比有大約2.5dB的增益,相同速率下,在信噪比為-2dB時,3比特TC-PIM誤碼率降低為2比特PIM的1/102倍左右,但卷積編碼和Viterbi算法會造成傳輸延遲[5]。

    針對上述差分易引起錯誤傳播、動態(tài)編譯碼表不一致和糾錯碼帶來傳輸延遲問題,本文提出一種改進差分的可逆變長碼壓縮編碼方法。對測井數(shù)據(jù)進行基于基準值的求差預處理,采用不等重保護策略,對處理后的差值數(shù)據(jù)進行可逆變長碼信源編碼,對基準值進行糾錯信道編碼,然后封幀傳輸。解碼端對數(shù)據(jù)進行解碼,對于錯誤的數(shù)據(jù)利用預測方式進行恢復,實現(xiàn)測井數(shù)據(jù)的實時可靠傳輸。

    2 隨鉆測井數(shù)據(jù)編碼及傳輸方式分析

    隨鉆測井無需中斷正常鉆進操作而將在鉆頭附近測量的多種信息傳輸?shù)降孛?,其?shù)據(jù)測量和環(huán)境具有相應特性,數(shù)據(jù)編碼及傳輸方式應根據(jù)數(shù)據(jù)特點和傳輸因素進行設計。

    2.1測井數(shù)據(jù)及傳輸特性分析

    井下地質測量參數(shù)有電阻率和方位伽馬,鉆井軌跡參數(shù)有井斜、方位、工具面,以8比特或12比特的數(shù)據(jù)形式進行編碼。井下發(fā)射器將傳感器測量的地質、軌跡等參數(shù)以0.5~3.0bit/s的傳輸速度通過泥漿信道傳輸?shù)降孛?。井下傳感器測量的每種數(shù)據(jù)都有其獨立性,但根據(jù)文獻[6]對測井數(shù)據(jù)的分析,井下大部分參數(shù)變化趨勢相對保持平穩(wěn),

    泥漿脈沖信號在傳輸過程主要受到信號衰減和噪聲干擾的影響。泥漿信號的衰減與泥漿信號的工作頻率、井深深度以及泥漿的密度等因素有關。泥漿脈沖信號符合指數(shù)衰減規(guī)律,泥漿壓力脈沖波幅的衰減模型為:

    [P=P0*exp(-xL)] (1)

    [L=Vd2ρπfu] (2)

    式(1)中:P0為信源信號的脈沖幅度,單位為Pa;x為信號傳輸距離,單位為m;L為信號源衰減到1/e時的傳輸距離,單位為m;式(2)中:P為傳輸x距離后的信號強度,單位是Pa;V為信號的傳輸速度,單位為m/s;d為管路內(nèi)徑,單位為m;ρ為鉆井液密度,單位是kg/m3;u為鉆井液粘度,Pa*s;f為脈沖信號的頻率,單位是Hz。由上式分析可知信號頻率與鉆井液粘度對信號衰減有顯著的影響,尤其是信號頻率,頻率越高信號衰減越大,隨鉆測量系統(tǒng)需要提供更多的能量來補償信號衰減。

    誤碼主要由噪聲產(chǎn)生,通常分以下幾類:1)泵沖噪聲。其為多頻噪聲,其基波頻率與有用信號的基頻很接近,具有周期性。2)扭矩噪聲。由堅硬的巖層、鉆頭、穩(wěn)定器等產(chǎn)生,可以通過改變鉆進參數(shù)使其噪聲頻率遠離有用信號頻率而消除。3)其他噪聲,其他噪聲可以近似為白噪聲。

    2.2數(shù)據(jù)采集編碼和傳輸方式

    為提高速率和可靠性,可以分別從數(shù)據(jù)采集編碼和傳輸過程采取措施。

    數(shù)據(jù)采集編碼采用基于基準值的動態(tài)分段差分處理方式,經(jīng)進一步量化處理后進行可逆變長碼壓縮編碼。因參數(shù)變化相對平穩(wěn)考慮采用鄰值差分方式去除數(shù)據(jù)時域冗余,但泥漿信道噪聲復雜易造成錯誤傳播,因而需對差分方式進行改進,采用基于基準值的動態(tài)分段差分處理。壓縮編碼中LZW編碼和動態(tài)統(tǒng)計編碼的編解碼表是動態(tài)生成的,數(shù)據(jù)傳輸錯誤會導致譯碼錯誤,因而考慮使用固定的編解碼表,算術編碼實現(xiàn)復雜,變長碼有唯一可譯的特性因而可以做錯誤自檢,不用對每一個數(shù)據(jù)加額外一個比特進行奇偶校驗,傳統(tǒng)的變長碼抗誤碼性能較差,因而考慮采用可逆變長碼方式進行壓縮編碼。

    傳輸采用非等重保護的幀結構和脈沖間隔調制。考慮到基準值和差值的重要性不均等,因而采用不等重保護策略對幀結構進行設計,設計關鍵幀和非關鍵幀兩種幀結構,其中,關鍵幀用前向糾錯碼方式保護基準值,而非關鍵幀利用可逆變長碼自身的抗誤碼能力對差值進行保護。目前泥漿波隨鉆傳輸系統(tǒng)中井下采用的是Drilog脈沖發(fā)生器,該發(fā)生器采用的編碼方式主要有曼切斯特編碼和脈沖間隔編碼[7]??紤]到曼切斯特編碼傳輸速率只能達到0.5~0.8bit/s,傳輸速率過低,因而采用脈沖間隔編碼進行調制。

    3基于改進差分的可逆變長碼壓縮方法

    改進差分采用基于基準值的動態(tài)分段差分方式對測井數(shù)據(jù)預處理,經(jīng)進一步對差值量化處理后進行可逆變長碼壓縮。

    3.1測井數(shù)據(jù)預處理

    為了減少數(shù)據(jù)變化范圍,在編碼前對測井數(shù)據(jù)進行預處理,主要包括基于基準值的差分預處理和差值量化兩部分。

    3.1.1基于基準值的差分預處理

    鄰值差分在數(shù)據(jù)傳輸錯誤時會引起錯誤傳播,除此之外,如果對差值量化,還會存在量化誤差,其積累可能導致解碼值偏離原始數(shù)據(jù)過大。因而對鄰值差分進行改進,提出一種基于基準值的動態(tài)分段差分處理方法,預設基準值和閾值,以后續(xù)數(shù)據(jù)和基準值的差值替代原始數(shù)據(jù)進行傳輸,當差值超過設定的閾值時對基準值進行更新處理。這里的基準值根據(jù)不同參數(shù)變化特點選取。

    對于升/降、平穩(wěn)趨勢的數(shù)據(jù),如圖1(a),存10個相鄰的歷史數(shù)據(jù),以最鄰近3個數(shù)據(jù)為起始數(shù)據(jù),采用滑動窗口法按時間點從近到遠依次查找最值點,當查找到最值點時停止查找,以時間鄰近值點和最值點計算斜率,當大于45度時判為升/降趨勢,以時間鄰近值加/減半?yún)^(qū)間為基準值,否則判為平穩(wěn)趨勢,對窗內(nèi)數(shù)據(jù)取平均作為基準值,當下一數(shù)據(jù)點不在基準值區(qū)間時視為趨勢誤判,以另外兩種方式取基準值,當都不在區(qū)間時,視為斷層做進一步處理。對于斷層,如圖1(b)所示,初判此點為異常點,以時間相鄰值替代輸出。當下一個數(shù)據(jù)在基準值區(qū)間時,此點判定為異常點,否則視為斷層點。下一個數(shù)據(jù)與此點絕對差值少于半?yún)^(qū)間時,以下一數(shù)據(jù)與此數(shù)據(jù)求取均值作為基準值。當不在半?yún)^(qū)間,但與此點值絕對差值大于半?yún)^(qū)間時,以下一個數(shù)據(jù)作為基準值。

    對于基本呈上升趨勢如圖1(c)的數(shù)據(jù),根據(jù)鄰值加/減半?yún)^(qū)間范圍為基準值對此時數(shù)據(jù)是否在基準值區(qū)間范圍進行判斷以確定基準值,如果都不在區(qū)間范圍,以斷層方式進行處理。

    3.1.2 差值量化

    不同參數(shù)精度要求不一,因而可以通過量化方式變更差分范圍。設定量化因子k,對于不同的輸入x有不同的輸出y,如式(3)所示。

    y=kx (3)

    3.2可逆變長碼

    可逆變長碼(RVLC)是指在滿足前后綴條件下可正反向解碼的變長碼,根據(jù)碼字是否滿足對稱條件,其可分為對稱和非對稱形式,對稱可逆變長碼在進行雙向解碼時只需要單一碼表,而非對稱可逆變長碼雙向解碼時需要兩套碼表,需要占用多一倍的存儲容量[8]。在目前的測井數(shù)據(jù)傳輸中,最大的瓶頸是傳輸帶寬而并非編解碼器設備的存儲容量,相對于對稱可逆變長碼,非對稱可逆變長碼的構造更為靈活,效率也更高,因此采用非對稱可逆變長碼對數(shù)據(jù)進行編碼。

    在RVLC碼樹的自頂向下的構造過程中,核心問題有兩個:第一個是各層應選取多少個葉子節(jié)點,第二個是葉子節(jié)點位置又該如何選取[9]。

    考慮到要實現(xiàn)8比特和12比特數(shù)據(jù)的碼表統(tǒng)一以及壓縮性能,只根據(jù)非對稱可逆變長碼的前后綴條件對-12~+12差值進行編碼,如表1所示。對于8位數(shù)據(jù)只用到-7~+7范圍的差值,而對于12位數(shù)據(jù)用到-12~+12范圍的差值。

    4幀結構與信道調制

    4.1幀結構設計

    根據(jù)不等重保護策略設置非關鍵幀和關鍵幀兩種不同的幀結構類型。非關鍵幀的結構如圖2所示,同步字節(jié)(FSB)包含1bit關鍵幀標志(KF)、1bit帶關鍵幀標志(CKF)、4bit序列號(SN)和2bit幀未填零個數(shù)(FTZN)四部分,其中關鍵幀標志位為0時表示非關鍵幀,為1時表示關鍵幀,帶關鍵幀標志為0時表示下一幀不是關鍵幀,否則下一幀為關鍵幀,序列號表示不同的參數(shù)組合,最大組合類型為16種,幀未填零個數(shù)是由于變長碼碼字長度的不確定,由于幀長是4的倍數(shù),因而需要在尾部填零,從而有利于可逆變長碼的反向解碼,接下來的可變長傳輸數(shù)據(jù)(VLTD)是不同參數(shù)的差值的變長碼字的組合,DC1~DCn各表示4bit數(shù)據(jù) ,F(xiàn)B表示的是填零比特數(shù),范圍從1~4。關鍵幀的結構如圖3所示,同步字節(jié)包含1bit的關鍵幀標志(KF)、3bit的基準值個數(shù)(NDV)、2bit的CRC校驗位(CRC)、2bit的幀未填零個數(shù)(FTZN)四部分,3bit的基準值個數(shù)代表幀內(nèi)數(shù)據(jù)的基準值個數(shù),2bit的CRC校驗位用于檢驗幀內(nèi)數(shù)據(jù)正確性,接下來的數(shù)據(jù)是多個基準值在參數(shù)組合中的位置和對應的值的組合(PBC),糾錯碼(FEC)和填零比特數(shù)(FB),其中糾錯碼用于保護基準值的數(shù)值和位置。

    4.2 脈沖間隔調制

    脈沖間隔調制(PIM)是一種利用兩個脈沖之間的時間間隔攜帶信息的調制方式,其符號如圖4所示。

    在圖4中,TPW( pulse width)表示脈沖寬度;TMP(min pulse time)表示最小脈沖間隔時間,用于減少碼間干擾,對應的數(shù)值M為0,TBW(bit width)表示位寬度,是數(shù)據(jù)編碼的時間間隔最小分辨值。PI碼編碼中脈沖間隔的計算公式如下:

    Interval=TMP+M×TBW (4)

    式4表明,當傳輸一個數(shù)值為M的數(shù)據(jù)時,經(jīng)過PI編碼后,數(shù)據(jù)對應的脈沖間隔的時長Interval的值為TMP與M個TBW時長的和。脈沖間隔通常分為兩類,一類是普通的數(shù)據(jù)位,用4比特表示,數(shù)據(jù)的取值范圍為0~15,可看作一個十六進制數(shù),另一類是同步字節(jié),用8比特表示,取值范圍為16~255。

    假設M在0,1,...,2N-1中等概率分布,N比特PIM的傳輸碼率為

    Tr(N,TMP,TBW)=[NTMP+2N-12TBW] (5)

    非等重的幀結構實現(xiàn)了差值和基準值的不等重保護,以脈沖間隔調制方式對數(shù)據(jù)進行傳輸,與曼切斯碼相比提升了數(shù)據(jù)的傳輸速率。

    5仿真與分析

    5.1壓縮比仿真與分析

    在Matlab仿真平臺上,通過基于改進差分的可逆變長碼壓縮算法分別對三口井的井斜角、方位角、電阻率和方位伽馬各2000個左右數(shù)據(jù)進行壓縮,從而對壓縮比和壓縮時間進行性能分析。

    通過表2可以看出,8bit的電阻率和伽馬數(shù)據(jù)壓縮比在60%左右,12bit井斜角和方位角的壓縮比在50%左右,表明壓縮性能隨著數(shù)據(jù)比特位數(shù)的增大而有所提升。文獻[3]中提出的差分脈碼調制壓縮編碼,8比特數(shù)據(jù)以3比特、4比特、5比特編碼時壓縮率分別是37.5%、50%、62.5%,12比特數(shù)據(jù)以5比特、6比特、7比特編碼時壓縮率分別是41.7%,50%,58.3%。與文獻[3]提出的方法對比,本壓縮算法在8比特數(shù)據(jù)時,和5比特編碼壓縮性能相當,12比特數(shù)據(jù)時,和6比特編碼壓縮性能相當。在壓縮時間方面,數(shù)據(jù)壓縮時間和壓縮比是成反比關系的,這是由于壓縮比越高,其基準值個數(shù)越少,從而減少了基準值變更所消耗的時間,由于井斜角是基本呈上升趨勢的,其基準值選取耗時較低,遠小于方位角、電阻率和伽馬值參數(shù)數(shù)據(jù)壓縮時間。

    5.2誤碼率仿真與分析

    為驗證非對稱可逆變長碼的抗誤碼性能,將非對稱可逆變長碼和常規(guī)變長碼編碼方案進行對比。發(fā)送端把可逆變長碼的4比特碼字轉換成十進制,組合一個數(shù)字序列,每次隨機從序列中抽取一個數(shù)據(jù)轉換成4比特數(shù),每抽取32次為一幀,共抽取500幀數(shù)據(jù),信道模型采用高斯白噪聲(AWGN)疊加周期為5Hz的泵沖噪聲信道,接收濾波器采用FIR濾波器。驗證在不同信噪比下的噪聲對常規(guī)變長碼和可逆變長碼編碼方式下的誤符號率,最終的抗誤碼性能用誤符號率(Pe)來衡量。

    從圖5可以看出,非稱可逆變長碼的誤符號性能明顯優(yōu)于常規(guī)的變長碼的。在疊加噪聲信道下,相同誤符號率下,非對稱可逆變長碼的信噪比變長碼低1~2dB,也就是說在相同誤符號率要求的條件下,非對稱可逆變長碼編碼可以增加傳輸距離,或者在相同的傳輸距離條件下,其可以降低誤符號率。

    6結論

    測井數(shù)據(jù)采集編碼采用基于基準值的動態(tài)分段差分的可逆變長碼壓縮方法,克服了鄰值差分的錯誤傳播和動態(tài)編譯碼的譯碼錯誤,經(jīng)仿真實驗得出壓縮率最大可達48%,提升了數(shù)據(jù)的傳輸速率。傳輸采用不等重保護策略對基準值進行糾錯編碼,對差分值利用可逆變長碼的自身抗誤碼能力編碼,減少了數(shù)據(jù)的傳輸延遲,利用脈沖間隔調制提升了數(shù)據(jù)的傳輸速率。雖然變長碼有錯誤自檢的功能,但也存在抗誤碼性能差的特性,經(jīng)過對變長碼和可逆變長碼抗誤碼性能分析可知,在同誤符號率下,非對稱可逆變長碼降低1~2dB的信噪比。因此,本算法在提高數(shù)據(jù)傳輸速率的同時也提高了數(shù)據(jù)傳輸?shù)目煽啃浴?/p>

    參考文獻:

    [1] 張煜.基于新型編碼技術的隨鉆測井實時數(shù)據(jù)傳輸方法研究[D].北京交通大學,2013.

    [2] 秦緒英,肖立志,索佰峰.隨鉆測井技術最新進展及其應用[J].勘探地球物理進展,2003,26(4):313-322.

    [3] 張煜,裘正定,熊軻,等.基于差分脈碼調制的隨鉆測量數(shù)據(jù)壓縮編碼算法[J].石油勘探與開發(fā),2010,37(6):748-754.

    [4] 李婷,史曉鋒.泥漿信道數(shù)據(jù)壓縮和糾錯聯(lián)合編碼研究[J].電子測量技術,2016,39(4);133-135.

    [5]陳偉,江佩佩,顧慶水,等.用于MWD遙傳系統(tǒng)的聯(lián)合編碼調制方法[J].測控技術,2016,35(9):33-35.

    [6]李忠偉,王瑞和,房軍.無線隨鉆測量系統(tǒng)的數(shù)據(jù)報協(xié)議[J].石油學報,2011,32(5);887-892.

    [7]江佩佩.泥漿波隨鉆傳輸系統(tǒng)中信號處理方法的研究[D].電子科技大學,2016.

    [8]S Kumar,L Xu.RVLC decoding scheme for improved data recovery in MPEG-4 video coding standard[J]Real-Time Imaging , 2004 ,10(5) :315-323.

    [9]范晨,崔慧娟,唐昆.一種高效可逆變長碼的構造算法[J].清華大學學報(自然科學版),2002,42(7):860-862.

    猜你喜歡
    基準值
    輪胎實車磨耗試驗數(shù)據(jù)分析方法的研究
    更正通知
    地表水環(huán)境質量標準綜述(三)
    ——標準銜接中常見問題及分析
    上海海岸帶沉積物地球化學基準值特征
    發(fā)電機組勵磁回路電阻基準值選取研究
    廣東韶關地區(qū)土壤地球化學基準值研究
    基于加權模糊C-均值聚類的鍋爐運行參數(shù)基準值建模
    馬拉硫磷水生生物基準推導研究
    復合材料基準值計算程序及其與可靠性關系
    復合材料B基準值計算程序
    国产精品偷伦视频观看了| 午夜福利在线免费观看网站| 飞空精品影院首页| 国产人伦9x9x在线观看| 国产熟女午夜一区二区三区| 99精国产麻豆久久婷婷| 男女午夜视频在线观看| 久久久久久免费高清国产稀缺| 久久综合国产亚洲精品| 亚洲成国产人片在线观看| 男女边摸边吃奶| cao死你这个sao货| 久久久水蜜桃国产精品网| 国产伦理片在线播放av一区| 人妻一区二区av| 在线永久观看黄色视频| 黑人巨大精品欧美一区二区蜜桃| 国产精品国产三级国产专区5o| 在线观看免费高清a一片| 免费在线观看日本一区| 亚洲欧洲日产国产| 日韩大片免费观看网站| 久久久精品国产亚洲av高清涩受| 最新的欧美精品一区二区| 性色av乱码一区二区三区2| 日韩三级视频一区二区三区| 1024香蕉在线观看| 国产精品久久久av美女十八| 99热全是精品| 亚洲 欧美一区二区三区| 99香蕉大伊视频| 少妇裸体淫交视频免费看高清 | 国产在线视频一区二区| 精品高清国产在线一区| 99久久国产精品久久久| 精品福利永久在线观看| 在线天堂中文资源库| 搡老岳熟女国产| 国产无遮挡羞羞视频在线观看| 夜夜骑夜夜射夜夜干| 国产老妇伦熟女老妇高清| 亚洲午夜精品一区,二区,三区| 精品久久久精品久久久| 久久热在线av| 亚洲综合色网址| 欧美日韩视频精品一区| 伊人亚洲综合成人网| 日韩视频一区二区在线观看| 视频区欧美日本亚洲| 1024香蕉在线观看| www.999成人在线观看| 午夜视频精品福利| 国产亚洲av高清不卡| 久久久久国产精品人妻一区二区| 欧美性长视频在线观看| 女人爽到高潮嗷嗷叫在线视频| 亚洲精品久久午夜乱码| 天堂8中文在线网| 国产精品一区二区免费欧美 | 丝袜在线中文字幕| 国产成人a∨麻豆精品| 亚洲欧洲日产国产| 性少妇av在线| 制服诱惑二区| 欧美xxⅹ黑人| 国产在线一区二区三区精| 精品少妇一区二区三区视频日本电影| 亚洲av成人一区二区三| 国产成人免费无遮挡视频| 正在播放国产对白刺激| 精品欧美一区二区三区在线| 国产人伦9x9x在线观看| 女人高潮潮喷娇喘18禁视频| 制服人妻中文乱码| 少妇的丰满在线观看| 热99re8久久精品国产| 两人在一起打扑克的视频| 波多野结衣一区麻豆| 亚洲七黄色美女视频| 国产亚洲精品久久久久5区| 亚洲人成77777在线视频| 欧美黄色淫秽网站| 国产极品粉嫩免费观看在线| 久久精品成人免费网站| 免费观看人在逋| 天天操日日干夜夜撸| 午夜福利在线免费观看网站| 亚洲精品美女久久av网站| 他把我摸到了高潮在线观看 | 久久国产亚洲av麻豆专区| 国产精品一区二区免费欧美 | 久久狼人影院| 少妇精品久久久久久久| 久久人人爽人人片av| 亚洲色图综合在线观看| 午夜久久久在线观看| 国产精品偷伦视频观看了| 人人妻人人澡人人爽人人夜夜| 韩国精品一区二区三区| 亚洲全国av大片| 97在线人人人人妻| 在线观看人妻少妇| 蜜桃国产av成人99| 亚洲专区中文字幕在线| 天天躁夜夜躁狠狠躁躁| 老熟女久久久| 91成人精品电影| 丝袜人妻中文字幕| 亚洲精品中文字幕一二三四区 | 深夜精品福利| 国产一卡二卡三卡精品| 19禁男女啪啪无遮挡网站| 51午夜福利影视在线观看| 两性夫妻黄色片| 一进一出抽搐动态| 国产精品秋霞免费鲁丝片| videosex国产| 精品国产乱码久久久久久男人| 最近中文字幕2019免费版| 国产精品av久久久久免费| 欧美在线一区亚洲| 一级,二级,三级黄色视频| 欧美日韩成人在线一区二区| 麻豆国产av国片精品| 美女主播在线视频| 18禁黄网站禁片午夜丰满| 一区在线观看完整版| 亚洲一卡2卡3卡4卡5卡精品中文| 国产欧美日韩综合在线一区二区| 夜夜夜夜夜久久久久| 亚洲成国产人片在线观看| 多毛熟女@视频| 国精品久久久久久国模美| 岛国在线观看网站| 亚洲精品国产av蜜桃| 国产精品久久久久久精品电影小说| 搡老乐熟女国产| 欧美大码av| av天堂久久9| 黄色a级毛片大全视频| 日日夜夜操网爽| 女人爽到高潮嗷嗷叫在线视频| 久久久久网色| 亚洲精品第二区| 精品亚洲成a人片在线观看| 久久亚洲精品不卡| 久久久国产精品麻豆| 国产三级黄色录像| 90打野战视频偷拍视频| 色94色欧美一区二区| 国产av精品麻豆| 午夜影院在线不卡| 国产真人三级小视频在线观看| 欧美在线黄色| 男女床上黄色一级片免费看| 大陆偷拍与自拍| 欧美精品一区二区大全| 国产成人影院久久av| 精品久久久久久电影网| 国产一区二区 视频在线| 黄色毛片三级朝国网站| 亚洲一区中文字幕在线| 婷婷色av中文字幕| 丝瓜视频免费看黄片| 国产主播在线观看一区二区| 90打野战视频偷拍视频| 日韩视频一区二区在线观看| 欧美日韩成人在线一区二区| 亚洲精品日韩在线中文字幕| 日韩熟女老妇一区二区性免费视频| 美女福利国产在线| 国产免费视频播放在线视频| 日韩免费高清中文字幕av| 国产视频一区二区在线看| 亚洲精品国产av蜜桃| 国产成人欧美| 国产真人三级小视频在线观看| 久久亚洲精品不卡| 国产欧美亚洲国产| 欧美精品一区二区大全| 超色免费av| 精品第一国产精品| 亚洲精品一卡2卡三卡4卡5卡 | 三上悠亚av全集在线观看| 捣出白浆h1v1| 国产色视频综合| 国产成人啪精品午夜网站| 精品国产国语对白av| 亚洲精品成人av观看孕妇| 性色av一级| 考比视频在线观看| 香蕉丝袜av| 最近中文字幕2019免费版| 啪啪无遮挡十八禁网站| 最近最新免费中文字幕在线| 欧美日韩视频精品一区| www.精华液| 国产一区二区 视频在线| 我要看黄色一级片免费的| 亚洲欧美一区二区三区久久| 两性午夜刺激爽爽歪歪视频在线观看 | 久9热在线精品视频| 亚洲国产av影院在线观看| 成人影院久久| 精品一区二区三区av网在线观看 | 日本一区二区免费在线视频| 伦理电影免费视频| 日本av手机在线免费观看| 女人爽到高潮嗷嗷叫在线视频| 欧美精品亚洲一区二区| 麻豆乱淫一区二区| 久久精品熟女亚洲av麻豆精品| 在线精品无人区一区二区三| 香蕉丝袜av| 丝袜人妻中文字幕| 日本一区二区免费在线视频| 亚洲成国产人片在线观看| 啦啦啦中文免费视频观看日本| 亚洲欧美清纯卡通| 激情视频va一区二区三区| 国产精品久久久av美女十八| 在线精品无人区一区二区三| 女人久久www免费人成看片| 欧美xxⅹ黑人| 国产色视频综合| 日本wwww免费看| 国产男女内射视频| 久久国产精品人妻蜜桃| 亚洲欧美一区二区三区久久| 亚洲av片天天在线观看| 美女国产高潮福利片在线看| 人人妻,人人澡人人爽秒播| 国产免费av片在线观看野外av| 99re6热这里在线精品视频| 人妻一区二区av| 99国产精品一区二区蜜桃av | 免费一级毛片在线播放高清视频 | 一本一本久久a久久精品综合妖精| 午夜福利免费观看在线| 后天国语完整版免费观看| 中文字幕最新亚洲高清| 男人爽女人下面视频在线观看| 亚洲 国产 在线| h视频一区二区三区| 岛国毛片在线播放| 亚洲精品久久久久久婷婷小说| 国产精品亚洲av一区麻豆| 精品一区二区三区av网在线观看 | 久久99一区二区三区| 母亲3免费完整高清在线观看| 性少妇av在线| 亚洲 欧美一区二区三区| 日韩 亚洲 欧美在线| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲精品美女久久av网站| 免费女性裸体啪啪无遮挡网站| 午夜福利在线观看吧| 两人在一起打扑克的视频| 黄片播放在线免费| www日本在线高清视频| 纯流量卡能插随身wifi吗| 国产真人三级小视频在线观看| 久久久久精品人妻al黑| 天天操日日干夜夜撸| 国产精品秋霞免费鲁丝片| 国产精品一区二区免费欧美 | 岛国毛片在线播放| 久久久久国内视频| 亚洲成国产人片在线观看| 美女福利国产在线| 美女视频免费永久观看网站| 人妻 亚洲 视频| 欧美xxⅹ黑人| 亚洲欧美成人综合另类久久久| 国产高清videossex| 免费日韩欧美在线观看| 亚洲全国av大片| 亚洲精品一二三| cao死你这个sao货| 欧美精品一区二区免费开放| av福利片在线| 精品国内亚洲2022精品成人 | 免费日韩欧美在线观看| 水蜜桃什么品种好| 久久久久久免费高清国产稀缺| 亚洲,欧美精品.| 美女脱内裤让男人舔精品视频| 欧美日韩一级在线毛片| 高清在线国产一区| 欧美人与性动交α欧美软件| 男女无遮挡免费网站观看| 成人手机av| 久久久久久久精品精品| 成人国产一区最新在线观看| 亚洲欧美精品自产自拍| 免费看十八禁软件| 国产成人啪精品午夜网站| 建设人人有责人人尽责人人享有的| 亚洲,欧美精品.| 国产成人影院久久av| 欧美精品av麻豆av| 老司机影院成人| 欧美在线黄色| 黄频高清免费视频| 成人国产av品久久久| 精品久久久久久电影网| 91九色精品人成在线观看| 黑人巨大精品欧美一区二区蜜桃| 亚洲久久久国产精品| 99国产精品一区二区蜜桃av | 一区福利在线观看| 亚洲欧美精品自产自拍| 精品一区二区三区四区五区乱码| 黄色怎么调成土黄色| 亚洲精品成人av观看孕妇| 亚洲av美国av| 满18在线观看网站| 国产97色在线日韩免费| 超碰成人久久| 老熟妇仑乱视频hdxx| 人妻一区二区av| 久久香蕉激情| 日日摸夜夜添夜夜添小说| 久久人人97超碰香蕉20202| 久久久精品免费免费高清| 国产一区二区激情短视频 | 精品久久久精品久久久| 人人妻,人人澡人人爽秒播| 黄色视频不卡| 欧美日韩成人在线一区二区| 国产精品免费大片| 欧美激情久久久久久爽电影 | 精品久久蜜臀av无| 啦啦啦视频在线资源免费观看| 妹子高潮喷水视频| 久久av网站| 欧美成狂野欧美在线观看| 91大片在线观看| 九色亚洲精品在线播放| 国产精品麻豆人妻色哟哟久久| 国产精品av久久久久免费| 国产一区二区在线观看av| 久久精品亚洲熟妇少妇任你| 亚洲精品久久成人aⅴ小说| 亚洲国产欧美一区二区综合| 亚洲第一欧美日韩一区二区三区 | 亚洲精品自拍成人| 国产高清国产精品国产三级| 少妇粗大呻吟视频| 人妻一区二区av| 最新的欧美精品一区二区| 国产成人av激情在线播放| 秋霞在线观看毛片| 一二三四社区在线视频社区8| 视频区图区小说| 啦啦啦视频在线资源免费观看| 另类精品久久| 亚洲男人天堂网一区| 中文字幕人妻熟女乱码| 亚洲欧美日韩高清在线视频 | 婷婷色av中文字幕| 99国产精品免费福利视频| 欧美人与性动交α欧美精品济南到| 人人妻人人添人人爽欧美一区卜| 丝袜在线中文字幕| 美女高潮到喷水免费观看| 99香蕉大伊视频| 色视频在线一区二区三区| 18禁观看日本| 国产精品av久久久久免费| 18禁裸乳无遮挡动漫免费视频| 精品少妇内射三级| 久久久久久久国产电影| 中文字幕另类日韩欧美亚洲嫩草| 视频区图区小说| 丁香六月欧美| 国产亚洲一区二区精品| 久热这里只有精品99| www.av在线官网国产| av福利片在线| 亚洲欧美一区二区三区久久| 国产精品麻豆人妻色哟哟久久| 亚洲av日韩在线播放| 久久久久精品国产欧美久久久 | 国产日韩欧美亚洲二区| 国产免费现黄频在线看| 夜夜夜夜夜久久久久| 欧美精品啪啪一区二区三区 | 免费高清在线观看视频在线观看| 亚洲国产中文字幕在线视频| 侵犯人妻中文字幕一二三四区| 久久亚洲国产成人精品v| 国产精品成人在线| 国产av国产精品国产| 欧美亚洲日本最大视频资源| 国产一区二区三区综合在线观看| 美女主播在线视频| 亚洲av成人一区二区三| 欧美在线黄色| 中国美女看黄片| 国精品久久久久久国模美| 99国产极品粉嫩在线观看| 欧美在线一区亚洲| 亚洲欧洲精品一区二区精品久久久| 秋霞在线观看毛片| 亚洲成av片中文字幕在线观看| 99香蕉大伊视频| 美女福利国产在线| 精品亚洲成a人片在线观看| 电影成人av| av在线播放精品| 国产亚洲精品第一综合不卡| 久久免费观看电影| 黄色怎么调成土黄色| 免费黄频网站在线观看国产| 1024视频免费在线观看| 亚洲中文日韩欧美视频| 91成人精品电影| 波多野结衣一区麻豆| 精品少妇内射三级| 人妻 亚洲 视频| 日本一区二区免费在线视频| 五月天丁香电影| 国产在线免费精品| 高清在线国产一区| 欧美在线黄色| 国产成人av教育| a 毛片基地| 亚洲男人天堂网一区| 亚洲熟女精品中文字幕| 欧美性长视频在线观看| 欧美在线黄色| 亚洲一卡2卡3卡4卡5卡精品中文| 美女国产高潮福利片在线看| 国产欧美日韩一区二区三 | 菩萨蛮人人尽说江南好唐韦庄| 秋霞在线观看毛片| 亚洲久久久国产精品| 国产麻豆69| 丝袜脚勾引网站| 99re6热这里在线精品视频| 久久久久国产一级毛片高清牌| 女人被躁到高潮嗷嗷叫费观| 亚洲精品中文字幕一二三四区 | 在线观看舔阴道视频| 国产xxxxx性猛交| 欧美成狂野欧美在线观看| a级毛片黄视频| 男女床上黄色一级片免费看| 18禁裸乳无遮挡动漫免费视频| 日韩一区二区三区影片| 亚洲一码二码三码区别大吗| netflix在线观看网站| 欧美激情 高清一区二区三区| 久久久精品国产亚洲av高清涩受| 老司机在亚洲福利影院| 亚洲国产日韩一区二区| 2018国产大陆天天弄谢| 91字幕亚洲| 三级毛片av免费| 亚洲精品一区蜜桃| 久久这里只有精品19| 自线自在国产av| 老汉色∧v一级毛片| 国产免费现黄频在线看| www.熟女人妻精品国产| 在线观看人妻少妇| 成年人黄色毛片网站| 老司机亚洲免费影院| 俄罗斯特黄特色一大片| 人人澡人人妻人| 宅男免费午夜| 欧美少妇被猛烈插入视频| 亚洲伊人色综图| 亚洲人成电影观看| 美国免费a级毛片| 久久精品人人爽人人爽视色| 伦理电影免费视频| 在线观看人妻少妇| 午夜精品久久久久久毛片777| 女性被躁到高潮视频| 一个人免费在线观看的高清视频 | 国产1区2区3区精品| 18禁国产床啪视频网站| 国产精品欧美亚洲77777| 久久精品熟女亚洲av麻豆精品| 少妇人妻久久综合中文| 纯流量卡能插随身wifi吗| 两性夫妻黄色片| 午夜两性在线视频| 人成视频在线观看免费观看| 99国产精品一区二区三区| 成在线人永久免费视频| 9191精品国产免费久久| 亚洲性夜色夜夜综合| 国产成人a∨麻豆精品| 久久久水蜜桃国产精品网| 黄片播放在线免费| 久久久精品免费免费高清| 成人三级做爰电影| 日本欧美视频一区| 中文字幕色久视频| 国产又爽黄色视频| 亚洲国产av新网站| 国产在线一区二区三区精| 午夜老司机福利片| 亚洲 国产 在线| 日韩制服骚丝袜av| 12—13女人毛片做爰片一| 18在线观看网站| 久久人人爽人人片av| 母亲3免费完整高清在线观看| 国产老妇伦熟女老妇高清| 成年美女黄网站色视频大全免费| 欧美国产精品va在线观看不卡| 狠狠婷婷综合久久久久久88av| 99九九在线精品视频| 国产成人精品在线电影| 国产av国产精品国产| www.av在线官网国产| 国精品久久久久久国模美| 在线精品无人区一区二区三| 9色porny在线观看| 国产伦理片在线播放av一区| 免费久久久久久久精品成人欧美视频| 日韩制服丝袜自拍偷拍| 后天国语完整版免费观看| 欧美成人午夜精品| 久久天躁狠狠躁夜夜2o2o| 一二三四在线观看免费中文在| 亚洲一区二区三区欧美精品| 狂野欧美激情性bbbbbb| 国产精品久久久久久人妻精品电影 | 黄片小视频在线播放| 中文欧美无线码| 在线 av 中文字幕| 日韩视频一区二区在线观看| 成年女人毛片免费观看观看9 | 自线自在国产av| 亚洲一卡2卡3卡4卡5卡精品中文| 亚洲精品国产区一区二| a在线观看视频网站| 在线 av 中文字幕| 国产亚洲精品久久久久5区| 午夜日韩欧美国产| 黄色视频不卡| 国产99久久九九免费精品| 一二三四社区在线视频社区8| 色视频在线一区二区三区| 国产欧美日韩精品亚洲av| 精品国产乱子伦一区二区三区 | 老汉色av国产亚洲站长工具| 久久午夜综合久久蜜桃| 他把我摸到了高潮在线观看 | 亚洲精品av麻豆狂野| 高清视频免费观看一区二区| 制服人妻中文乱码| 另类亚洲欧美激情| 欧美国产精品一级二级三级| 欧美变态另类bdsm刘玥| 国产男女内射视频| 欧美在线黄色| 免费在线观看视频国产中文字幕亚洲 | 精品久久久久久久毛片微露脸 | 超碰97精品在线观看| 亚洲精品久久久久久婷婷小说| 亚洲熟女毛片儿| 久久精品国产亚洲av高清一级| 两性午夜刺激爽爽歪歪视频在线观看 | 成人三级做爰电影| 午夜影院在线不卡| 国产精品亚洲av一区麻豆| av在线app专区| 亚洲免费av在线视频| 飞空精品影院首页| 国产熟女午夜一区二区三区| 国产一级毛片在线| 国产av又大| 飞空精品影院首页| 最黄视频免费看| 女人被躁到高潮嗷嗷叫费观| 视频区图区小说| 国产精品影院久久| 日韩制服骚丝袜av| 亚洲成人免费av在线播放| 国产区一区二久久| 涩涩av久久男人的天堂| e午夜精品久久久久久久| 亚洲精品一区蜜桃| 性少妇av在线| 免费黄频网站在线观看国产| 爱豆传媒免费全集在线观看| 亚洲国产毛片av蜜桃av| 国产男人的电影天堂91| 日韩精品免费视频一区二区三区| 2018国产大陆天天弄谢| 国产精品久久久人人做人人爽| 精品一区在线观看国产| 91麻豆精品激情在线观看国产 | 50天的宝宝边吃奶边哭怎么回事| 亚洲久久久国产精品| 一区二区三区四区激情视频| 日韩熟女老妇一区二区性免费视频| 亚洲欧美一区二区三区久久| 精品国产乱码久久久久久男人| 不卡av一区二区三区| 啪啪无遮挡十八禁网站| 亚洲欧美精品综合一区二区三区| 丝袜人妻中文字幕| 亚洲精品在线美女| 操出白浆在线播放| 亚洲精品一区蜜桃| 操美女的视频在线观看| 亚洲欧美色中文字幕在线| 免费av中文字幕在线|