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

    基于泰勒級數(shù)估計(jì)的油井?dāng)?shù)據(jù)無損壓縮

    2016-11-22 11:35:06丁英濤陳欣劉箭言李怡然
    關(guān)鍵詞:冪函數(shù)階數(shù)級數(shù)

    丁英濤, 陳欣, 劉箭言, 李怡然

    (北京理工大學(xué) 信息與電子學(xué)院, 北京 100081)

    ?

    基于泰勒級數(shù)估計(jì)的油井?dāng)?shù)據(jù)無損壓縮

    丁英濤, 陳欣, 劉箭言, 李怡然

    (北京理工大學(xué) 信息與電子學(xué)院, 北京 100081)

    為了實(shí)現(xiàn)油井?dāng)?shù)據(jù)的高效傳輸,提出一種新的無損壓縮算法. 利用泰勒級數(shù)分解擬合出油井?dāng)?shù)據(jù)曲線,進(jìn)行后向估計(jì),通過傳輸擬合值與實(shí)際值的估計(jì)誤差,實(shí)現(xiàn)數(shù)據(jù)的無損壓縮. 實(shí)測油井?dāng)?shù)據(jù)仿真表明該算法壓縮率可達(dá)25%~40%,其整體性能優(yōu)于霍夫曼編碼、LZW編碼等無損壓縮算法至少20%,并具有時(shí)間空間復(fù)雜度低的特點(diǎn). 通過大港油田數(shù)據(jù)遠(yuǎn)程傳輸系統(tǒng)驗(yàn)證,該算法可將傳輸網(wǎng)絡(luò)數(shù)據(jù)負(fù)荷降低至45%.

    無損壓縮;冪函數(shù)擬合;泰勒級數(shù);后向估計(jì);油井?dāng)?shù)據(jù)

    隨著油田信息化程度的加深,油田井場遠(yuǎn)程控制系統(tǒng)采集到的電參數(shù)、載荷、位移等信息的數(shù)據(jù)量不斷增加. 由于受到油田地理環(huán)境的限制,這些信息主要通過無線網(wǎng)絡(luò)傳輸. 一種是內(nèi)部專用的無線網(wǎng)絡(luò),如ZigBee網(wǎng)絡(luò);另一種是第三方運(yùn)營網(wǎng)絡(luò),如GPRS或者3G網(wǎng)絡(luò). 前者帶寬較小,難以負(fù)荷越來越大的傳輸量;后者帶寬大,但屬于收費(fèi)網(wǎng)絡(luò),占用民用資源[1]. 因此,直接傳輸?shù)姆绞?,已?jīng)不適用于數(shù)字化油田大數(shù)據(jù)量的需求. 目前,已經(jīng)有一些文獻(xiàn)致力于研究油井?dāng)?shù)據(jù)壓縮方法,文獻(xiàn)[2-3]提出的Ramer-Douglas-Peucker (RDP)算法,通過只保留曲線上的足以表達(dá)該曲線形狀特征的關(guān)鍵點(diǎn),來實(shí)現(xiàn)數(shù)據(jù)壓縮,但數(shù)據(jù)精度在刪除點(diǎn)處存在非常明顯的誤差,且屬于有損壓縮. 文獻(xiàn)[4]中提出一種對大量數(shù)據(jù)做奇異值分解,用得到的基向量來近似擬合數(shù)據(jù),達(dá)到數(shù)據(jù)壓縮的目的,但算法計(jì)算量大,無法滿足油井?dāng)?shù)據(jù)的實(shí)時(shí)傳輸. 此外,還有一些經(jīng)典通用的數(shù)據(jù)壓縮算法:霍夫曼壓縮編碼、算術(shù)編碼、游程編碼、LZW算法等[5],這些算法都沒有考慮油井?dāng)?shù)據(jù)的特點(diǎn),不能有效實(shí)現(xiàn)油井?dāng)?shù)據(jù)的壓縮.

    油井?dāng)?shù)據(jù)在時(shí)域上表現(xiàn)為平穩(wěn)緩變的信號,數(shù)據(jù)波形具有很強(qiáng)的相關(guān)性. 本文設(shè)計(jì)一種基于泰勒級數(shù)分解和冪函數(shù)擬合理論的數(shù)據(jù)壓縮方法,根據(jù)油井采集到數(shù)據(jù)的特征,快速確定達(dá)到最佳估計(jì)的泰勒級數(shù)階數(shù),再通過后向估計(jì)逐次得到后續(xù)數(shù)據(jù)的估計(jì)值,傳輸過程只需要傳輸估計(jì)誤差,從而實(shí)現(xiàn)數(shù)據(jù)壓縮. 接收端通過相同階數(shù)的泰勒級數(shù)分解和冪函數(shù)擬合,利用接收到的估計(jì)誤差序列,即可無損還原出原始數(shù)據(jù).

    1 算法原理

    1.1 算法原理概述

    在數(shù)學(xué)上,常用一組簡單的基函數(shù)來逼近復(fù)雜函數(shù),實(shí)現(xiàn)級數(shù)分解,泰勒級數(shù)就是以冪函數(shù)為基函數(shù),逼近一個(gè)復(fù)雜函數(shù).

    泰勒級數(shù)[6]:設(shè)實(shí)變(復(fù)變)函數(shù)f(x)在x0的某一鄰域上存在n+1階連續(xù)導(dǎo)數(shù),則對該領(lǐng)域內(nèi)的f(x)可以分解為如下的冪級數(shù):

    (1)

    根據(jù)上述理論,實(shí)際采樣得到的一個(gè)N點(diǎn)序列{yn|n∈[1,N]},認(rèn)為是抽象連續(xù)函數(shù)f(x)的N個(gè)采樣點(diǎn),那么該函數(shù)必然可以實(shí)現(xiàn)一個(gè)M階的泰勒級數(shù)擬合,每個(gè)點(diǎn)yn存在的誤差為Δen,即N點(diǎn)序列滿足下列冪函數(shù)擬合關(guān)系:

    (2)

    式中ai=f(i)(x0)/i!表示冪函數(shù)的系數(shù).

    擬合函數(shù)和實(shí)際函數(shù)之間存在一個(gè)估計(jì)誤差Δe,在擬合的階數(shù)最接近函數(shù)實(shí)際階數(shù)時(shí)Δe最小,否則擬合的階數(shù)增大或減小都會(huì)導(dǎo)致估計(jì)誤差Δe變大.

    在傳輸過程中只傳輸估計(jì)誤差序列{en}及初始數(shù)據(jù)序列{yn|n∈[1:K]}.

    接收端利用K個(gè)初始點(diǎn),根據(jù)相同的函數(shù)系數(shù)進(jìn)行相同階數(shù)的擬合,恢復(fù)估計(jì)值,再用接收到的誤差序列對估計(jì)值進(jìn)行修正,即可無損還原數(shù)據(jù).

    根據(jù)上述原理,所需要傳輸數(shù)據(jù)的幅值位寬將減少,當(dāng)M?N時(shí),壓縮率近似可以達(dá)到

    (3)

    1.2 算法過程詳述

    下面對抽象函數(shù)y=f(x)進(jìn)行M階泰勒級數(shù)分解,推導(dǎo)壓縮端和解壓縮端的詳細(xì)算法原理.

    為了更好地體現(xiàn)yn是N點(diǎn)的時(shí)域采樣序列,將抽象函數(shù)描述為

    式中Ts為采樣周期.

    式(4)為M階冪函數(shù)的表達(dá)式:

    (4)

    其中所需要確定的變量如下:

    ① 函數(shù)的最佳階數(shù)M;

    ② 函數(shù)的系數(shù)aM~a0;為了方便后續(xù)的分析,記做系數(shù)向量:

    (5)

    根據(jù)求解線性方程的Cramer法則,要確定M+1個(gè)系數(shù)需要M+1個(gè)獨(dú)立的線性方程,所以引入假設(shè)條件:

    可以將式(2)修正為

    (6)

    改寫成矩陣形式:

    (7)

    式中y=[y1y2…yM+1]T,

    因?yàn)閤n=nTs,n∈[1,N],所以X為范德蒙矩陣,只與擬合階數(shù)M有關(guān)系,后面將記做XM+1,滿足滿秩條件:r(XM+1)=M+1,存在唯一的逆矩陣,方程組存在唯一的解,即系數(shù)向量

    (8)

    根據(jù)上式求出的系數(shù),對前M+1個(gè)數(shù)據(jù)的估計(jì)誤差為0,對第M+2點(diǎn)的估計(jì)值為

    (9)

    將式(8)帶入式(9),得

    (10)

    (11)

    估計(jì)誤差為

    (12)

    此后估計(jì)第M+i點(diǎn)時(shí),直接由式(11)得到

    (13)

    因此對于一個(gè)有N個(gè)采樣點(diǎn)的數(shù)據(jù)序列,在確定了擬合冪函數(shù)階數(shù)M后,就可以根據(jù)已知的M+1個(gè)初始點(diǎn)對后一點(diǎn)數(shù)據(jù)進(jìn)行估計(jì),然后得到相應(yīng)點(diǎn)的估計(jì)誤差,得估計(jì)誤差序列:

    (14)

    在傳輸過程中只需傳輸初始的M+1個(gè)點(diǎn)的數(shù)據(jù)序列Yinit=[y1y2…yM+1]和估計(jì)誤差序列{E},即待傳輸?shù)臄?shù)據(jù)序列為

    (15)

    (16)

    重復(fù)上述解壓過程,最終可實(shí)現(xiàn)數(shù)據(jù)還原.

    對大量傳輸數(shù)據(jù)來說,最佳擬合階數(shù)M?N,如果擬合冪函數(shù)的階數(shù)選擇合理,則誤差幅度值很小,因此該方法可以有效降低傳輸數(shù)據(jù)的量化位數(shù).

    1.3 最佳階數(shù)的選取

    在實(shí)際工程中,由于外界環(huán)境影響和設(shè)備損耗等問題,采集到的數(shù)據(jù)存在一定的變化. 此時(shí)需要在井下處理端,自適應(yīng)的選擇能達(dá)到最佳估計(jì)的泰勒級數(shù)階數(shù)M.

    從傳感器采集到數(shù)據(jù)開始,截取一個(gè)周期,利用本文方法進(jìn)行壓縮計(jì)算,從M=1開始,得到一個(gè)暫時(shí)的最大估計(jì)誤差幅值:

    (17)

    計(jì)算后續(xù)M+1階的ΔE,依次循環(huán),當(dāng)最大估計(jì)誤差幅值由減小的趨勢變?yōu)樵龃髸r(shí),即可跳出循環(huán),確定最佳泰勒級數(shù)的擬合階數(shù). 并在后續(xù)工作的一段時(shí)間內(nèi)(如:24 h),都采用上述獲得的最佳階數(shù)M進(jìn)行壓縮.

    2 結(jié)果及分析

    2.1 實(shí)測數(shù)據(jù)軟件仿真

    為驗(yàn)證本文方法的有效性,下面對油井傳輸數(shù)據(jù)中數(shù)據(jù)量較大的電參數(shù),利用本文方法做壓縮仿真分析. 油井?dāng)?shù)據(jù)多為低頻信號,采集系統(tǒng)通常是低通過采樣的,含有大量的高頻噪聲以及帶外雜波,因此在利用本文方法進(jìn)行壓縮前應(yīng)該對數(shù)據(jù)進(jìn)行預(yù)濾波處理,如圖1.

    通過對比圖,可以看出電參數(shù)數(shù)據(jù)經(jīng)過低通濾波之后,不存在突變點(diǎn),波形平滑且相關(guān)性強(qiáng),滿足描述的壓縮算法的適用條件.

    下面分別利用不同階數(shù)的冪函數(shù)對電參數(shù)數(shù)據(jù)進(jìn)行壓縮,得到的壓縮結(jié)果如圖2. 可以看出,從M=1到M=4,傳輸估計(jì)誤差的幅值逐漸減小,最小減至-12~12之間,從M≥5開始,傳輸估計(jì)誤差的幅值開始增大,因此M=4為電參數(shù)數(shù)據(jù)的最佳擬合階數(shù),根據(jù)式(3)可得壓縮率R為37.5%. 油井主要傳輸數(shù)據(jù)的壓縮率R和壓縮最佳階數(shù)M總結(jié)為表1,考慮實(shí)際工程應(yīng)用,傳輸數(shù)據(jù)的量化位數(shù)盡量使用2的整數(shù)倍,表2、表3中的壓縮結(jié)果也有相同考慮.

    參數(shù)原始傳輸數(shù)據(jù)/bit壓縮傳輸數(shù)據(jù)/bit壓縮率(R)/%最佳階數(shù)(M)載荷124(±3)33.32功率166(±12)37.54電流124(±4)33.32電壓124(±2)33.31位移124(±2)33.31

    2.2 算法復(fù)雜度分析

    算法效率的度量最常見的是時(shí)間復(fù)雜度和空間復(fù)雜度,時(shí)間復(fù)雜度是指算法執(zhí)行完畢所需要的時(shí)間,空間復(fù)雜度是指算法在運(yùn)行過程中臨時(shí)占用的輔助空間[7].

    本文算法的計(jì)算主要是乘法和加法,那么最壞時(shí)間復(fù)雜度可表示為T(N)=O(N). 時(shí)間復(fù)雜度隨O(N)輸入序列的個(gè)數(shù)N成線性關(guān)系,算法的時(shí)間效率比較好,滿足快速運(yùn)算的要求.

    計(jì)算過程中,算出的估計(jì)誤差,可直接存儲(chǔ)在相應(yīng)輸入點(diǎn)的位置,覆蓋原始數(shù)據(jù)后直接傳輸,所以本文算法屬于原地工作,不需要額外的臨時(shí)存儲(chǔ)空間,那么空間復(fù)雜度可以表示為S(N)=O(1). 空間復(fù)雜度為常數(shù),不隨輸入序列個(gè)數(shù)的變化而變化,執(zhí)行過程中不會(huì)占用過多的資源.

    2.3 算法魯棒性分析及改進(jìn)

    冪函數(shù)擬合是后向估計(jì),每一個(gè)估計(jì)出來的點(diǎn),會(huì)作為已知點(diǎn)估計(jì)下一個(gè)點(diǎn),如果傳輸過程中發(fā)生數(shù)據(jù)錯(cuò)誤,誤差具有累積效應(yīng). 如果擾動(dòng)出現(xiàn)在第一個(gè)計(jì)算點(diǎn),那么整個(gè)輸出估計(jì)誤差序列都會(huì)受影響,并且越往后造成的影響越大,去除已知點(diǎn)數(shù)據(jù),影響范圍可表示為N-K.

    為了改進(jìn)算法的魯棒性,可將輸入序列等分成L份,并且從每一份的中間點(diǎn)開始向兩邊同時(shí)用本文算法進(jìn)行前向和后向估計(jì)壓縮,前向估計(jì)壓縮原理同后向估計(jì),只是數(shù)據(jù)序列的處理方向相反,在接收端倒置即可恢復(fù)原始數(shù)據(jù). 如果擾動(dòng)仍出現(xiàn)在第一個(gè)計(jì)算點(diǎn),對整體的影響范圍可降低到(N/L-K)/2,如圖3所示.

    改進(jìn)之后算法魯棒性有所增強(qiáng),但同時(shí)會(huì)帶來壓縮率下降的問題. 在實(shí)際工程中,需要根據(jù)實(shí)際情況和要求,平衡兩者關(guān)系,選取適合的L,而本文算法均采用原始算法計(jì)算得出.

    2.4 算法壓縮結(jié)果比較

    傳統(tǒng)經(jīng)典的霍夫曼編碼、差分編碼、LZW編碼都可以實(shí)現(xiàn)數(shù)據(jù)的無損壓縮,與這些典型的算法相比,本文算法能實(shí)現(xiàn)更低的壓縮率,對多組數(shù)據(jù)進(jìn)行壓縮,統(tǒng)計(jì)結(jié)果對比如表2.

    表2 壓縮率對比

    2.5 實(shí)際系統(tǒng)應(yīng)用結(jié)果

    2012年在天津大港采油廠,安裝了40套基于Cortex-M3 ARM處理器的油井參數(shù)采集和遠(yuǎn)程傳輸設(shè)備,采集和遠(yuǎn)程傳輸壓力、溫度、載荷、位移和電參數(shù)數(shù)據(jù),在ARM處理器中實(shí)現(xiàn)并且現(xiàn)場應(yīng)用了上述算法,表3列舉了其中10口油井主要數(shù)據(jù)壓縮的結(jié)果.

    表3 10口井實(shí)測數(shù)據(jù)壓縮效果

    3 結(jié) 論

    針對油井?dāng)?shù)據(jù)量越來越大的問題,提出了一種基于泰勒級數(shù)分解和線性估計(jì)的油井?dāng)?shù)據(jù)無損壓縮方法. 該算法對油井?dāng)?shù)據(jù)的每一點(diǎn)利用泰勒級數(shù)分解和冪函數(shù)進(jìn)行擬合估計(jì),在傳輸過程中只傳輸初始序列和估計(jì)誤差序列,達(dá)到數(shù)據(jù)壓縮的目的. 在接收端,利用初始序列和估計(jì)誤差序列對數(shù)據(jù)進(jìn)行還原和修正,實(shí)現(xiàn)數(shù)據(jù)的無損還原. 對實(shí)測數(shù)據(jù)進(jìn)行仿真,壓縮率可達(dá)25%~40%,相比經(jīng)典的霍夫曼編碼、差分編碼以及LZW算法,壓縮率最少20%. 本文算法成功應(yīng)用在天津大港采油廠的實(shí)際系統(tǒng)中,有效降低了采油廠RTU傳輸數(shù)據(jù)量,將總的傳輸網(wǎng)絡(luò)負(fù)荷降低至45%. 此外,該無損壓縮算法也可適用于滿足平滑緩變的其他數(shù)據(jù)類型.

    [1] 馬曉.油井監(jiān)測數(shù)據(jù)傳輸系統(tǒng)的設(shè)計(jì)與實(shí)現(xiàn)[D].西安:西安電子科技大學(xué),2013.

    Ma Xiao. Design and implementation of oil well data transmission system[D]. Xi’an: Xidian University, 2013. (in Chinese)

    [2] Ramer U. An iterative procedure for the polygonal approximation of plane curves[J]. Computer Graphics and Image Processing, 1972(1):244-256.

    [3] Douglas D H, Peucker T K. Algorithms for the reduction of the number of points required to represent a digitized line or its caricature[J]. The Canadian Cartographer, 1973,10(2):1122122.

    [4] Hatami S, Feldmann P, Abbaspour S, et al. Efficient compression and handling of current source model library waveforms[C]∥Proceedings of Design Automantion and Test in Europe Conference and Exhibition. Dresden, Germany: IEEE Press, 2009:1178-1183.

    [5] 傅祖蕓.信息論[M].3版.北京:電子工業(yè)出版社,2011:273-303.

    Fu Zuyun. Information theory[M]. 3rd ed. Beijing: Publishing House of Electronics Industry, 2011:273-303. (in Chinese)

    [6] 毛京中.高等數(shù)學(xué)教程[M].北京:高等教育出版社,2008:198-211.

    Mao Jingzhong. Advanced mathematics[M]. Beijing: Higher Education Press, 2008:198-211. (in Chinese)

    [7] 嚴(yán)蔚敏,吳偉民.數(shù)據(jù)結(jié)構(gòu)[M].C語言版.北京:清華大學(xué)出版社,2009:13-17.

    Yan Weimin, Wu Weimin. Data structure[M]. C edition. Beijing: Tsinghua University Press, 2009:13-17. (in Chinese)

    (責(zé)任編輯:劉芳)

    Oil Well Data Lossless Compression Based on Taylor Series Estimating

    DING Ying-tao, CHEN Xin, LIU Jian-yan, LI Yi-ran

    (School of Information and Electronics, Beijing Institute of Technology, Beijing 100081, China)

    A lossless compression algorithm applied on the oil well data transmission was proposed in this paper. Taylor series expansion was used to fit the curve of the oil well data, then the backward estimation was used for data processing. The lossless compression was realized by transmitting the fitting estimation error and initial sequence. From the simulation results, data compression ratio is up to 25%~40% with a low time-space complexity. Compared with other typical coding, like Huffman coding and LZW coding and Delta coding, this compression ratio has 20% advanced. Based on the verification of the actual oilfield RTU, this design can significantly decrease the burden of the data-transmitting net to 45%.

    lossless compression; power function fitting; Taylor series; backward estimating; oil well data

    2014-01-02

    丁英濤(1972—),女,副教授,E-mail:dyt@bit.edu.cn.

    TN911

    A

    1001-0645(2016)05-0530-05

    10.15918/j.tbit1001-0645.2016.05.017

    猜你喜歡
    冪函數(shù)階數(shù)級數(shù)
    冪函數(shù)、指數(shù)函數(shù)、對數(shù)函數(shù)(2)
    冪函數(shù)、指數(shù)函數(shù)、對數(shù)函數(shù)(1)
    關(guān)于無窮小階數(shù)的幾點(diǎn)注記
    確定有限級數(shù)解的階數(shù)上界的一種n階展開方法
    冪函數(shù)、指數(shù)函數(shù)、對數(shù)函數(shù)(1)
    Dirichlet級數(shù)及其Dirichlet-Hadamard乘積的增長性
    看圖說話,揭開冪函數(shù)的廬山真面目
    幾個(gè)常數(shù)項(xiàng)級數(shù)的和
    p級數(shù)求和的兩種方法
    Dirichlet級數(shù)的Dirichlet-Hadamard乘積
    国产精品一区二区在线不卡| 在线看a的网站| 国产欧美日韩综合在线一区二区| xxxhd国产人妻xxx| 婷婷成人精品国产| 精品人妻熟女毛片av久久网站| 亚洲精华国产精华精| 一区二区av电影网| 丝袜美足系列| 精品人妻1区二区| 岛国毛片在线播放| 欧美日韩福利视频一区二区| 国产一区二区 视频在线| 在线 av 中文字幕| 在线观看免费视频网站a站| 午夜福利在线免费观看网站| 黄色片一级片一级黄色片| 日日夜夜操网爽| 婷婷成人精品国产| 欧美日本中文国产一区发布| 久久精品亚洲av国产电影网| 一区二区三区四区激情视频| 伊人久久大香线蕉亚洲五| 精品少妇黑人巨大在线播放| 最黄视频免费看| 亚洲精品久久午夜乱码| 少妇粗大呻吟视频| 久久精品人人爽人人爽视色| 桃花免费在线播放| 超碰成人久久| av视频免费观看在线观看| 最近最新中文字幕大全免费视频| 亚洲综合色网址| 国产精品久久久久成人av| 免费黄频网站在线观看国产| 精品免费久久久久久久清纯 | 国产精品偷伦视频观看了| 欧美国产精品一级二级三级| 亚洲三区欧美一区| 十八禁网站网址无遮挡| 男人操女人黄网站| 少妇猛男粗大的猛烈进出视频| 不卡一级毛片| 18禁观看日本| 久久免费观看电影| 男女免费视频国产| 男人操女人黄网站| 十分钟在线观看高清视频www| www.熟女人妻精品国产| 又黄又粗又硬又大视频| 黑人巨大精品欧美一区二区蜜桃| 一区福利在线观看| 老汉色av国产亚洲站长工具| 美女脱内裤让男人舔精品视频| 人人澡人人妻人| 亚洲欧美日韩另类电影网站| 中文字幕人妻丝袜一区二区| 国产欧美日韩一区二区精品| 日日夜夜操网爽| 中文字幕制服av| 老汉色av国产亚洲站长工具| 色婷婷av一区二区三区视频| 99久久99久久久精品蜜桃| 成人国产一区最新在线观看| 少妇裸体淫交视频免费看高清 | 啦啦啦在线免费观看视频4| 国产无遮挡羞羞视频在线观看| 亚洲国产精品999| 18在线观看网站| 大片电影免费在线观看免费| 色老头精品视频在线观看| 亚洲第一av免费看| 久久久久国内视频| av在线播放精品| 午夜视频精品福利| 香蕉丝袜av| 日韩中文字幕视频在线看片| 777久久人妻少妇嫩草av网站| av欧美777| 欧美精品一区二区大全| 青春草亚洲视频在线观看| 国产精品一二三区在线看| 亚洲国产精品999| 波多野结衣一区麻豆| 久久青草综合色| 12—13女人毛片做爰片一| 久久久水蜜桃国产精品网| av福利片在线| 欧美国产精品va在线观看不卡| 亚洲第一av免费看| 中文欧美无线码| av在线播放精品| 成人国产一区最新在线观看| 日韩欧美免费精品| 婷婷丁香在线五月| 建设人人有责人人尽责人人享有的| 黄片小视频在线播放| 中文欧美无线码| 中文字幕av电影在线播放| 伊人久久大香线蕉亚洲五| 大陆偷拍与自拍| 老司机午夜福利在线观看视频 | 久久久精品区二区三区| 国产一区二区三区在线臀色熟女 | 国产91精品成人一区二区三区 | 国产av国产精品国产| 嫩草影视91久久| 天天添夜夜摸| 久久ye,这里只有精品| 激情视频va一区二区三区| 亚洲久久久国产精品| 女警被强在线播放| 欧美日韩av久久| 国产又色又爽无遮挡免| 香蕉国产在线看| 一级毛片精品| 欧美97在线视频| 日韩大码丰满熟妇| 国产成人免费无遮挡视频| 亚洲精品一区蜜桃| 最新的欧美精品一区二区| 久久青草综合色| 99国产极品粉嫩在线观看| 日韩欧美免费精品| 久久久久久久久免费视频了| 丝瓜视频免费看黄片| 99精国产麻豆久久婷婷| 亚洲精品国产色婷婷电影| 欧美亚洲日本最大视频资源| 精品免费久久久久久久清纯 | 亚洲成国产人片在线观看| 女性被躁到高潮视频| 国产又色又爽无遮挡免| 99热国产这里只有精品6| 亚洲精品国产av成人精品| h视频一区二区三区| 一区二区三区精品91| av在线播放精品| 精品亚洲成a人片在线观看| 纵有疾风起免费观看全集完整版| 亚洲欧美日韩高清在线视频 | 国产男女超爽视频在线观看| 大码成人一级视频| 亚洲精品中文字幕在线视频| 国产男女超爽视频在线观看| 国产高清国产精品国产三级| 国产亚洲av片在线观看秒播厂| 男女下面插进去视频免费观看| 久久国产精品人妻蜜桃| 日本91视频免费播放| 老鸭窝网址在线观看| 国产一区二区 视频在线| 嫩草影视91久久| 最近最新中文字幕大全免费视频| 免费在线观看日本一区| 日韩电影二区| 精品熟女少妇八av免费久了| 男人爽女人下面视频在线观看| 亚洲视频免费观看视频| 精品一区二区三区四区五区乱码| 亚洲成人免费电影在线观看| 中文字幕另类日韩欧美亚洲嫩草| 欧美久久黑人一区二区| 成年av动漫网址| 亚洲国产精品一区二区三区在线| 精品国产乱子伦一区二区三区 | 男女高潮啪啪啪动态图| 18禁观看日本| 男女边摸边吃奶| 黄色怎么调成土黄色| av线在线观看网站| 母亲3免费完整高清在线观看| 91老司机精品| 新久久久久国产一级毛片| 最近最新免费中文字幕在线| 色播在线永久视频| 国产av一区二区精品久久| 欧美日韩一级在线毛片| 久久国产精品影院| videosex国产| 高清av免费在线| 黄色 视频免费看| 香蕉丝袜av| 成人影院久久| 免费黄频网站在线观看国产| 欧美人与性动交α欧美精品济南到| 咕卡用的链子| 亚洲视频免费观看视频| 中文字幕色久视频| 男人舔女人的私密视频| 亚洲第一av免费看| 欧美人与性动交α欧美精品济南到| av不卡在线播放| 久久中文看片网| 亚洲国产av新网站| 国产亚洲一区二区精品| 精品福利永久在线观看| 欧美日韩国产mv在线观看视频| 久久久欧美国产精品| 国产极品粉嫩免费观看在线| 亚洲五月色婷婷综合| 成人亚洲精品一区在线观看| 青青草视频在线视频观看| 久久久国产欧美日韩av| 波多野结衣一区麻豆| 免费在线观看视频国产中文字幕亚洲 | 亚洲熟女精品中文字幕| 老熟妇乱子伦视频在线观看 | 精品一区在线观看国产| 国产亚洲av片在线观看秒播厂| 超色免费av| 夜夜骑夜夜射夜夜干| 免费在线观看完整版高清| 十分钟在线观看高清视频www| 日韩中文字幕欧美一区二区| 精品国产一区二区三区久久久樱花| 欧美另类一区| 涩涩av久久男人的天堂| 欧美人与性动交α欧美精品济南到| 免费日韩欧美在线观看| 久久国产精品大桥未久av| 久久久精品94久久精品| tube8黄色片| 久久九九热精品免费| 亚洲自偷自拍图片 自拍| 中文字幕人妻丝袜一区二区| 久久久久精品国产欧美久久久 | 国产熟女午夜一区二区三区| 一边摸一边做爽爽视频免费| 男女之事视频高清在线观看| 在线观看免费高清a一片| 国产主播在线观看一区二区| 久久久久久亚洲精品国产蜜桃av| 悠悠久久av| 久久久久久久久久久久大奶| 老司机深夜福利视频在线观看 | 我要看黄色一级片免费的| 91精品三级在线观看| 久久久国产一区二区| 久久久久久人人人人人| 国产色视频综合| 老司机深夜福利视频在线观看 | av不卡在线播放| 亚洲精品美女久久av网站| 人妻久久中文字幕网| 不卡一级毛片| 18在线观看网站| 老司机深夜福利视频在线观看 | 国产精品久久久久成人av| 十八禁网站免费在线| 免费久久久久久久精品成人欧美视频| av在线老鸭窝| 一进一出抽搐动态| 性色av一级| 99国产精品一区二区蜜桃av | 久久影院123| 日本91视频免费播放| 欧美人与性动交α欧美精品济南到| 老司机在亚洲福利影院| 少妇裸体淫交视频免费看高清 | 他把我摸到了高潮在线观看 | 不卡av一区二区三区| 80岁老熟妇乱子伦牲交| 菩萨蛮人人尽说江南好唐韦庄| 日韩一区二区三区影片| 久久久久久久久久久久大奶| 99久久综合免费| 亚洲第一青青草原| 欧美中文综合在线视频| 日本91视频免费播放| 久久久国产成人免费| 黄色视频在线播放观看不卡| 亚洲欧美精品自产自拍| 欧美乱码精品一区二区三区| 国产精品一区二区在线观看99| 搡老乐熟女国产| 欧美97在线视频| 免费不卡黄色视频| 久久精品国产a三级三级三级| tocl精华| 亚洲欧美激情在线| 丝袜人妻中文字幕| a 毛片基地| 国内毛片毛片毛片毛片毛片| 热re99久久国产66热| 久久毛片免费看一区二区三区| 日韩中文字幕视频在线看片| 韩国精品一区二区三区| 69精品国产乱码久久久| 妹子高潮喷水视频| 国产精品成人在线| 动漫黄色视频在线观看| 这个男人来自地球电影免费观看| 亚洲专区国产一区二区| 中亚洲国语对白在线视频| 夜夜骑夜夜射夜夜干| 国产免费现黄频在线看| 热re99久久国产66热| www.999成人在线观看| 高清在线国产一区| 少妇被粗大的猛进出69影院| 欧美国产精品一级二级三级| 精品人妻在线不人妻| 嫁个100分男人电影在线观看| 欧美日韩av久久| 51午夜福利影视在线观看| 中文字幕高清在线视频| 亚洲中文av在线| 黄色怎么调成土黄色| 1024香蕉在线观看| 丰满人妻熟妇乱又伦精品不卡| 久久国产精品男人的天堂亚洲| 亚洲欧美成人综合另类久久久| 丝袜在线中文字幕| 香蕉国产在线看| 99热全是精品| 免费看十八禁软件| 国产精品一区二区精品视频观看| 久久免费观看电影| 又紧又爽又黄一区二区| 午夜福利在线观看吧| 老司机深夜福利视频在线观看 | 色老头精品视频在线观看| 国产精品av久久久久免费| 亚洲成人手机| 国产免费现黄频在线看| 王馨瑶露胸无遮挡在线观看| 99精国产麻豆久久婷婷| 国产亚洲精品一区二区www | 91精品三级在线观看| 久久午夜综合久久蜜桃| 亚洲av美国av| 黄色a级毛片大全视频| a级片在线免费高清观看视频| 美女扒开内裤让男人捅视频| 国产淫语在线视频| 亚洲avbb在线观看| 欧美变态另类bdsm刘玥| 一级毛片电影观看| 欧美变态另类bdsm刘玥| 中国国产av一级| 欧美午夜高清在线| 9色porny在线观看| 国产欧美日韩精品亚洲av| 9色porny在线观看| 一二三四社区在线视频社区8| 午夜久久久在线观看| 99久久人妻综合| 在线观看免费午夜福利视频| 亚洲av成人不卡在线观看播放网 | 在线观看人妻少妇| 可以免费在线观看a视频的电影网站| 国产日韩欧美亚洲二区| 色综合欧美亚洲国产小说| 国产精品麻豆人妻色哟哟久久| 国产97色在线日韩免费| 巨乳人妻的诱惑在线观看| 十分钟在线观看高清视频www| 十八禁高潮呻吟视频| 性少妇av在线| 新久久久久国产一级毛片| 一区二区av电影网| 伊人久久大香线蕉亚洲五| bbb黄色大片| 女人被躁到高潮嗷嗷叫费观| 国产成人av教育| 欧美+亚洲+日韩+国产| 国产成人av教育| 自线自在国产av| 免费在线观看黄色视频的| 久久久久久久国产电影| 国产成人av教育| www.自偷自拍.com| 性少妇av在线| 日本欧美视频一区| 国产精品 欧美亚洲| 国产av一区二区精品久久| 男女高潮啪啪啪动态图| 深夜精品福利| xxxhd国产人妻xxx| 美女高潮到喷水免费观看| 亚洲欧美日韩高清在线视频 | 一区在线观看完整版| 999久久久精品免费观看国产| 欧美成狂野欧美在线观看| 动漫黄色视频在线观看| 精品熟女少妇八av免费久了| 老熟妇仑乱视频hdxx| 亚洲熟女精品中文字幕| www.自偷自拍.com| 一本—道久久a久久精品蜜桃钙片| 12—13女人毛片做爰片一| 精品国产乱码久久久久久小说| 亚洲 国产 在线| 人人妻人人澡人人看| 国产精品影院久久| 精品欧美一区二区三区在线| av欧美777| 精品一区在线观看国产| 99re6热这里在线精品视频| 脱女人内裤的视频| 黑丝袜美女国产一区| 丰满少妇做爰视频| 成年美女黄网站色视频大全免费| 久久精品亚洲熟妇少妇任你| 色精品久久人妻99蜜桃| 一区二区三区激情视频| 免费在线观看黄色视频的| av线在线观看网站| 亚洲国产欧美在线一区| av视频免费观看在线观看| 久久影院123| 老汉色av国产亚洲站长工具| 99久久99久久久精品蜜桃| 中文字幕人妻丝袜制服| 一边摸一边抽搐一进一出视频| 中亚洲国语对白在线视频| 亚洲一区中文字幕在线| 欧美日韩成人在线一区二区| 欧美 日韩 精品 国产| 欧美日韩福利视频一区二区| 午夜免费成人在线视频| 高清av免费在线| 热99久久久久精品小说推荐| 啦啦啦中文免费视频观看日本| 成人亚洲精品一区在线观看| 热99国产精品久久久久久7| cao死你这个sao货| 秋霞在线观看毛片| 男人操女人黄网站| 女人高潮潮喷娇喘18禁视频| 美女扒开内裤让男人捅视频| 深夜精品福利| 国产免费一区二区三区四区乱码| 99国产精品一区二区蜜桃av | 亚洲一区中文字幕在线| 桃红色精品国产亚洲av| 老司机靠b影院| 国产国语露脸激情在线看| 99国产极品粉嫩在线观看| 成年女人毛片免费观看观看9 | 午夜激情av网站| 亚洲精品日韩在线中文字幕| 婷婷成人精品国产| 桃花免费在线播放| 午夜福利一区二区在线看| 青春草亚洲视频在线观看| 啦啦啦 在线观看视频| 国产成人精品久久二区二区免费| 亚洲一码二码三码区别大吗| 中文字幕精品免费在线观看视频| 国产精品.久久久| 天堂8中文在线网| 最近中文字幕2019免费版| 久久久精品94久久精品| h视频一区二区三区| 老司机福利观看| 久久久久久亚洲精品国产蜜桃av| 啦啦啦啦在线视频资源| 国产福利在线免费观看视频| 一本综合久久免费| 免费女性裸体啪啪无遮挡网站| 久久国产亚洲av麻豆专区| 久久久水蜜桃国产精品网| 我的亚洲天堂| 亚洲,欧美精品.| 欧美另类一区| 日韩欧美一区视频在线观看| av一本久久久久| 国产又色又爽无遮挡免| 法律面前人人平等表现在哪些方面 | 日韩精品免费视频一区二区三区| 美女扒开内裤让男人捅视频| avwww免费| 日本猛色少妇xxxxx猛交久久| 国产在线观看jvid| 可以免费在线观看a视频的电影网站| 99国产极品粉嫩在线观看| 99香蕉大伊视频| 久久久久久免费高清国产稀缺| 三级毛片av免费| 久久久国产欧美日韩av| 国产高清国产精品国产三级| 91成年电影在线观看| 一边摸一边抽搐一进一出视频| 青草久久国产| 丁香六月天网| 中国国产av一级| 久久久国产精品麻豆| 久久久久网色| 一区在线观看完整版| 在线观看人妻少妇| 精品一区在线观看国产| 色94色欧美一区二区| 手机成人av网站| 波多野结衣一区麻豆| 免费黄频网站在线观看国产| 国产精品一区二区在线不卡| 久久精品人人爽人人爽视色| 777久久人妻少妇嫩草av网站| 国产欧美日韩一区二区精品| 午夜影院在线不卡| 国产日韩一区二区三区精品不卡| 黑丝袜美女国产一区| 国产精品1区2区在线观看. | 日韩,欧美,国产一区二区三区| 欧美日韩国产mv在线观看视频| 最近最新免费中文字幕在线| 日韩制服丝袜自拍偷拍| 美女扒开内裤让男人捅视频| 岛国毛片在线播放| 国产精品久久久久久精品电影小说| 国产精品偷伦视频观看了| 极品少妇高潮喷水抽搐| 一级毛片电影观看| 国产精品一区二区在线观看99| 成年人午夜在线观看视频| 后天国语完整版免费观看| 亚洲少妇的诱惑av| 青春草亚洲视频在线观看| 久久精品人人爽人人爽视色| 水蜜桃什么品种好| 三级毛片av免费| 一区二区日韩欧美中文字幕| 中文字幕高清在线视频| 男女边摸边吃奶| 久久久久久久国产电影| 亚洲欧洲日产国产| 午夜福利免费观看在线| 日韩 亚洲 欧美在线| bbb黄色大片| 精品人妻在线不人妻| 美女福利国产在线| 大型av网站在线播放| 欧美亚洲日本最大视频资源| 人妻久久中文字幕网| 日韩三级视频一区二区三区| 亚洲国产精品成人久久小说| 中文欧美无线码| 欧美精品高潮呻吟av久久| 我的亚洲天堂| 亚洲一码二码三码区别大吗| 亚洲精品国产av蜜桃| 午夜成年电影在线免费观看| 一个人免费在线观看的高清视频 | 少妇裸体淫交视频免费看高清 | 欧美 亚洲 国产 日韩一| 成人手机av| 国产成人系列免费观看| 黄片播放在线免费| 在线看a的网站| 纯流量卡能插随身wifi吗| 一级片'在线观看视频| 日本黄色日本黄色录像| 9热在线视频观看99| av福利片在线| 考比视频在线观看| 大片免费播放器 马上看| 久久久久国内视频| 丝袜喷水一区| 精品国内亚洲2022精品成人 | 久久精品国产亚洲av香蕉五月 | 亚洲色图综合在线观看| 久久久久久久久久久久大奶| 亚洲伊人色综图| 国产一区二区 视频在线| 精品免费久久久久久久清纯 | 欧美激情久久久久久爽电影 | 真人做人爱边吃奶动态| 视频在线观看一区二区三区| 国产高清videossex| av欧美777| 国精品久久久久久国模美| a级毛片在线看网站| 女人高潮潮喷娇喘18禁视频| 黄色a级毛片大全视频| 国产成人精品无人区| 中文字幕制服av| 伦理电影免费视频| 狠狠婷婷综合久久久久久88av| 妹子高潮喷水视频| 熟女少妇亚洲综合色aaa.| 女人久久www免费人成看片| av免费在线观看网站| 国产成人免费观看mmmm| 久久av网站| 色老头精品视频在线观看| 久久久精品国产亚洲av高清涩受| 国产又色又爽无遮挡免| 王馨瑶露胸无遮挡在线观看| 国产色视频综合| 日韩精品免费视频一区二区三区| 99久久99久久久精品蜜桃| 午夜福利在线观看吧| 韩国精品一区二区三区| 亚洲情色 制服丝袜| 亚洲欧美精品自产自拍| 久久人人97超碰香蕉20202| 久久久久精品国产欧美久久久 | 热re99久久国产66热| 91国产中文字幕| 啦啦啦免费观看视频1| 亚洲成人手机| 国产色视频综合| 在线观看免费视频网站a站| 丰满少妇做爰视频| 国产熟女午夜一区二区三区| 欧美97在线视频| 精品国产国语对白av| 如日韩欧美国产精品一区二区三区| 在线观看免费视频网站a站| 人妻人人澡人人爽人人| 男人操女人黄网站| 国产精品国产av在线观看| 国产成+人综合+亚洲专区|