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

    基于L-M算法和聲波融合的TDOA雷電定位改進方法

    2022-03-29 07:30:04駱瑤瑩卞宏志劉全楨劉保全傅正財張建勛劉亞坤
    上海交通大學學報 2022年3期
    關鍵詞:子站雷聲聲波

    駱瑤瑩, 卞宏志, 劉全楨, 劉保全, 傅正財, 張建勛, 劉亞坤

    (1. 上海交通大學 電力傳輸與功率變換控制教育部重點實驗室, 上海 200030;2. 中國石化青島安全工程研究院 化學品安全控制國家重點實驗室, 山東 青島 266071;3. 國網福建省電力有限公司建設分公司, 福州 350001)

    雷電作為一種長間隙強放電現(xiàn)象,產生時伴有強烈的電磁和聲波信號[1],嚴重影響油罐儲運、電力設備、航空航天活動和森林安全等[2-4].雷電定位技術基于雷電產生時的電、磁、聲和光等信號特征,通過分析不同信號的傳播特性,反演計算得到雷電的位置信息和產生時間[5].雷電定位技術是目前用于研究雷電活動最主要的手段之一[6],為關鍵系統(tǒng)的雷電防護和雷電災害分析提供數據保障[7].

    現(xiàn)有雷電定位系統(tǒng)多測量雷電放電時發(fā)出的電磁脈沖信號,其測量探頭可工作在甚高頻、甚低頻/低頻或極低頻等頻段,依據多個探頭測量得到的雷電電磁脈沖的波峰到達時間和各探頭間的時間差,計算確定雷電電磁輻射源的位置和產生時間.由該方法測得的雷電信息主要受探測子站探頭的時鐘精度、電磁環(huán)境和反演定位算法等影響.在雷電定位算法研究中,趙文光等[8]綜合利用到達時差定位法(Time Difference of Arrival, TDOA)和方向定位法,通過考慮平差計算剔除粗差的算法,提高了雷電定位的精度.胡志祥等[9]在雷電定位計算中引入粒子群優(yōu)化算法,得到了更精確的雷電發(fā)生位置.郭小紅等[10]提出優(yōu)化果蠅算法,實現(xiàn)了雷電定位結果的快速計算.李濤等[11]利用改進密度聚類算法對雷電位置進行聚類解算,提高了雷電定位的抗誤差干擾能力.梁麗等[12]綜合密度聚類算法和網格搜索法,降低了電磁探測數據誤差等對雷電定位的影響.在實際的TDOA雷電定位系統(tǒng)中,探測子站存在測量誤差,因此采用3個及3個以上探測子站的雷電定位系統(tǒng)會出現(xiàn)雷電位置信息難以求解或求解不唯一的問題.如何解決雷電定位計算中的求解問題,且減小探測站與雷電相對位置以及定位網絡幾何尺寸對計算結果的影響,具有實際應用意義.

    安裝于我國油罐等系統(tǒng)中的雷電定位系統(tǒng)具備同時采集聲波和電磁信息的能力.為充分利用聲波和電磁信號, 李皖等[13]利用人工識別雷擊電磁脈沖信號法得到了電磁信號與雷聲信號的時間差;湯偉等[14]將雷聲按點聲源處理,利用坐標搜索法進行雷聲方向定位.綜合已有研究發(fā)現(xiàn),利用雷電放電時的聲波和電磁信號同時進行雷電定位可以提高定位精度,但在雷電電磁和聲波信號的融合交互方面還需進一步研究.因此,本文提出一種基于Levenberg-Marquardt算法(簡稱L-M算法)和聲波信息融合的TDOA雷電定位改進方法.根據已有探測站點的布站信息,將目標雷電監(jiān)測區(qū)域劃分為16個子區(qū)域;利用電磁信號的到達時差,在整個監(jiān)測區(qū)域內計算雷電初定位結果,并將其作為依據,動態(tài)納入雷聲站點測得的聲波信息,避免聲波信號在長距離傳播時的大幅衰減和畸變.在計算求解雷電發(fā)生位置時,利用可融合聲波信息分析的L-M迭代算法,得到更精確的最優(yōu)解.所提方法在TDOA電磁信息中納入聲波信息,可以幫助解決雷電定位計算中的求解問題,提升雷電定位系統(tǒng)的定位精度和抗誤差擾動能力.

    1 TDOA雷電定位方法

    基于TDOA雷電定位系統(tǒng),利用3個及3個以上探測子站測得的電磁信息,計算雷電電磁脈沖到達各探測子站的時間差,以確定電磁輻射源的位置.TDOA雷電定位系統(tǒng)一般由連接在同一個通信網絡上的若干探測站、中心站和用戶終端組成[15].其中,探測站主要包括電磁脈沖接收天線、信號處理單元、高精度時鐘單元、電源和通信等模塊.伴隨雷電產生的電磁信號到達接收天線,并由信號處理單元處理獲得同步的時間標簽后,通信口將電磁脈沖波峰到達探測站的精確時間實時發(fā)送到中心站.中心站的雷電位置分析儀根據電磁脈沖到達各探頭的時間差來確定雷電電磁輻射源的位置,并由中心站將定位結果輸出到各用戶工作站.在保證系統(tǒng)運行可靠性的基礎上,多個定位系統(tǒng)聯(lián)網可以實現(xiàn)數據共享,滿足雷電定位系統(tǒng)的擴展需求[16].

    在雷電電磁信號測量中,雷電輻射源產生的電磁波到達探測站#1和探測站#2的時間差為Δt12,結合兩個探測站的地理位置信息,可以確定一條雙曲線,其計算式為

    (1)

    式中:(x1,y1)和(x2,y2)分別為TDOA探測子站T1和T2的位置坐標;c為電磁波傳播速度,且c=3×108m/s.根據雙曲線交匯原理,電磁波輻射源可以由交匯點位置坐標確定.在3個及3個以上探測站的情況中,不同雙曲線相互交匯并產生公共交點,據此可以確定雷電發(fā)生位置.

    在TDOA雷電定位的實際計算中,雙曲線發(fā)散的幾何形狀特性導致整個雷電監(jiān)測區(qū)域的定位精度不均勻[17].采用3站TDOA定位可能會出現(xiàn)兩條雙曲線交于兩點的情況,導致雷電定位系統(tǒng)作出偽雷擊位置判斷[18].對于3站以上的TDOA定位系統(tǒng),電磁波探測站在測量時存在誤差,因此利用探測站之間的時間差信息得到的多條雙曲線不易交于一點,從而影響基于到達時差信息的雷電定位方程組求解.對此,通常利用優(yōu)化算法來減小測量誤差的影響并解決方程組的求解問題[8-12].

    2 基于L-M算法和聲波信息融合的TDOA雷電定位改進方法

    2.1 L-M迭代算法

    L-M算法是牛頓法的一種改進.在計算時,其能夠避免牛頓法由于Hessian矩陣奇異而導致算法無法繼續(xù)迭代的情況[19],常用于非線性最小二乘問題的最優(yōu)化實現(xiàn),能夠有效解決時差信息冗余時雷電位置最優(yōu)解的求解問題.

    L-M算法用于雷電定位計算時的流程如圖1所示.其中,L0為雷電初定位結果,即L-M算法的迭代初值.在求解由多個探測子站的聲波信息得到的非線性方程組f(x)=0時,可轉換為最小值尋找問題:

    圖1 利用L-M迭代算法進行雷電定位計算的流程圖Fig.1 Flow chart of lightning location calculation using L-M iterative algorithm

    (2)

    式中:x為需求解的雷電定位結果.高斯-牛頓法將迭代步長設定為

    (3)

    式中:I為單位矩陣;μ>0保證系數矩陣正定,從而保證迭代的下降方向,克服了因Hk不可逆而導致的迭代無法繼續(xù)問題.為保證算法在接近最優(yōu)解時能夠快速收斂[20],當μ較大時,

    當μ較小時,

    hlm≈hgn

    其中:τ為限定μ0在適當數值的系數.μ的數值更新由系數?決定:

    (4)

    式中:F(xk)-F(xk+hlm)表示目標函數在步長hlm下的實際變化;L(0)-L(hlm)表示目標函數的預估變化值,其計算式為

    (5)

    (6)

    ‖xnew-x‖<ε2(‖x‖+ε2)

    (7)

    k>kmax

    (8)

    式中:ε1和ε2為迭代終止的條件參數;kmax為最大迭代步數,合適的kmax值可以防止迭代無限循環(huán)引起的算法崩潰.ε1、ε2和kmax主要根據L-M算法的具體應用情況進行選取,對算法最終收斂結果影響不大.在本算法中,設ε1=ε2=10-12,kmax=200.在計算過程中,3個參數的取值滿足上述條件之一即終止算法,輸出雷電定位計算結果.當探測站測得的信號存在誤差或信息冗余時,會導致雷電定位計算無法得到非線性方程的唯一解,而L-M迭代算法可以為尋找雷電位置最優(yōu)解提供可靠的算法支持.

    2.2 融合聲波信息的TDOA定位

    TDOA雷電定位系統(tǒng)對測試精度要求較高,其中,探測站布局方式會引起時間誤差[21],雷電電磁波峰值點因傳播路徑和距離的不同而發(fā)生漂移或畸變,導致時間測量誤差[22],使得TDOA的定位誤差幅值達到幾百米或幾千米.此外,整個雷電監(jiān)測目標區(qū)域的TDOA雷電定位精度不均勻,存在局部定位誤差較大的區(qū)域塊.雷電發(fā)生時會發(fā)出強烈的電磁信號和聲波信號.聲波信號的傳播特性與信號頻率有關,在相同距離和環(huán)境背景下,信號頻率越低則聲波傳播性能越優(yōu)異[23].隨著傳播距離的增大,低于 100 Hz 的聲波信號幅度衰減速度較慢,可以作為有效的雷聲信息進行雷電定位研究.因此,可以考慮在進行雷電定位時納入雷聲聲波信息來減小局部區(qū)域塊的定位誤差,從而提升整個雷電監(jiān)測目標區(qū)域的雷電定位效果.

    在TDOA雷電定位系統(tǒng)融合聲波信息時,利用波形延時壓縮處理技術可以有效發(fā)送聲波信號,減小聲波信號傳遞帶來的影響[14].為提高所納入雷聲聲波信息的準確性,綜合根據已有TDOA探測子站和雷聲探測站的位置信息,將整個雷電監(jiān)測區(qū)域進行子區(qū)域劃分處理,如圖2所示.其中,T3和T4為TDOA探測子站,N1~N4為4個雷聲探測站.共劃分為16個雷電監(jiān)測子區(qū)域.其中,子區(qū)域#1、#4、#13和#16為近時差站子區(qū)域,子區(qū)域#6、#7、#10 和#11為近雷聲站子區(qū)域,子區(qū)域#2、#3、#5、#8、#9、#12、#14和#15為時差-雷聲站點近似等距子區(qū)域.以TDOA探測子站與子站之間的連線(藍色線)為邊界,將線內區(qū)域定義為探測網內區(qū)域,線外區(qū)域即監(jiān)測邊緣區(qū)域定義為探測網外區(qū)域.在處理雷聲信號時,將極短時間內連續(xù)的多點聲源近似處理為單個定位計算所需要的雷聲點聲源信號[13].在接收雷電放電時,聲波信號雷聲站中的聲波麥克風(M0~M4)陣列布置方式為十字陣列排布,如圖3所示.

    圖2 TDOA探測子站和雷聲探測站布局Fig.2 Site location of TDOA lightning location sensors and acoustic sensors

    圖3 十字麥克風陣列示意圖Fig.3 A diagram of crossed microphone array

    當雷電發(fā)生在目標監(jiān)測區(qū)域內時,4個TDOA探測子站會采集到一套相應的雷電電磁信息,依據測得的雷電電磁時間信息,可列出由3個互相獨立方程組成的方程組:

    (9)

    式中:(x3,y3)和(x4,y4)分別為T3和T4的位置坐標.

    利用式(9)計算得到雷電初定位結果(單一點),記為L0=(x0,y0).根據初定位坐標判斷雷電發(fā)生位置所在的子區(qū)域.同時,利用子區(qū)域最小距離判定并選取合適的麥克風陣列,進行麥克風陣列所探雷聲信息的融合和求解.根據測得的雷聲聲波信息、麥克風陣列相對位置和聲波傳播特性,得到

    (10)

    式中:(xi,yi)和(xj,yj)分別為麥克風Mi和Mj的坐標;Δtij為雷聲到達時差;v為聲波傳播速度.

    由此可知,對于包含多個麥克風的測量陣列系統(tǒng),如由5個麥克風組成的測量陣列系統(tǒng)可以得到10個約束方程,即使在少數麥克風傳感器工作失效的特殊情況下,其仍能夠獲得足夠的雷聲信息.納入雷聲聲波信息后,方程組中的方程數目大于未知數數目,可以采用L-M算法進行迭代計算,求出雷電位置的最優(yōu)解.

    基于L-M算法和聲波信息融合的TDOA雷電定位算法的整體流程如圖4所示.利用TDOA探測子站測得的電磁信號及其時間信息,判定雷電初定位結果.根據L0所屬子區(qū)域k′的判斷結果,調取距離子區(qū)域k′最近的雷聲站n′測得的聲波信息,并將其納入非線性方程組,利用L-M迭代算法計算更精確的雷電位置.

    圖4 基于L-M算法和聲波信息融合的TDOA雷電定位算法流程Fig.4 Flow chart of TDOA lightning location approach considering acoustic information and L-M algorithm

    2.3 算法驗證

    在MATLAB中完成所提方法的算法編寫,并利用蒙特卡羅法[24]對雷電定位法進行精度評估.以邊長為40 km的正方形雷電監(jiān)測目標區(qū)域為例,據圖2中x軸和y軸建立坐標系4個時差法探測子站 T1~T4探頭位置的坐標分別為(15, 15)(-15, 15)(-15,-15)和(15,-15);4個雷聲站N1~N4的中心麥克風位置的坐標分別為(8, 8)(-8, 8)(-8,-8)和(8,-8);圖3中M1~M4與M0的距離均為20 m.

    依據蒙特卡羅法,以2 km為空間步長將 40 km×40 km的雷電監(jiān)測區(qū)域進行網格劃分,將區(qū)域內的441(21×21)個格點分別作為雷電發(fā)生位置進行雷電定位精度評估,具體步驟如下:

    步驟1計算在某格點處產生的雷電電磁波和聲波信號到達每個探頭的準確時間.

    步驟3分別利用TDOA定位的傳統(tǒng)法和所提改進法進行雷電位置計算,得到相應的雷電定位結果.

    步驟4評估雷電定位結果的誤差,即雷電發(fā)生的實際位置與定位計算結果之間的絕對差值.

    步驟5重復步驟2至步驟4共100次,得到以該格點作為雷電發(fā)生位置的100次定位結果和定位誤差,并將100次定位誤差的均值作為該格點的最終定位誤差值.

    圖5 雷電定位誤差分布Fig.5 Error distribution of lightning location

    表1 不同子區(qū)域的雷電定位誤差均值

    TDOA雷電定位傳統(tǒng)法和改進法在不同子區(qū)域內的誤差均值分布如圖6所示.結合圖5可知,改進法利用雷電監(jiān)測區(qū)域測得的電磁時差信息開展雷電活動的初步定位,根據初步定位結果動態(tài)選取距離最近的雷聲站信息,同時融合電磁和聲波信息進行雷電活動精準定位,以改善整個監(jiān)測區(qū)域的雷電定位效果.傳統(tǒng)法在探測網內的定位誤差為120~200 m,在探測網外的定位誤差接近300 m,在整個目標監(jiān)測區(qū)域內的定位誤差均值為203.2 m,均方差為82.8;改進改法在探測網內的定位誤差為90~110 m,在整個目標監(jiān)測區(qū)域內的定位誤差均值為108.4 m,定位精度提高了46.6%,均方差為27.3.此外,利用傳統(tǒng)法定位時,位于監(jiān)測區(qū)域邊緣的雷電定位誤差明顯增大;而利用改進法進行定位時,位于監(jiān)測區(qū)域邊緣的雷電定位誤差得到改善,雷電定位精度比傳統(tǒng)法提高了51.2%.

    圖6 不同子區(qū)域的雷電定位誤差均值分布Fig.6 Average lightning location error distribution in different sub-regions

    3 結語

    本文提出一種基于L-M算法和聲波信息融合的TDOA雷電定位改進方法,根據TDOA探測子站和雷聲探測站的位置信息,將目標監(jiān)測區(qū)域劃分為16個子區(qū)域.利用電磁時差信息完成雷電的初步定位計算,并根據初定位結果動態(tài)納入雷聲站探測得到的聲波信息.同時,利用L-M算法對融合的聲波信息進行計算,求解雷電位置最優(yōu)值.采用蒙特卡羅法分別對TDOA雷電定位傳統(tǒng)法和所提改進法進行精度評估,發(fā)現(xiàn)傳統(tǒng)TDOA雷電定位系統(tǒng)的定位誤差在整個監(jiān)測區(qū)域分布不均,在目標區(qū)域邊緣的雷電定位精度下降明顯,誤差均值為203.2 m,所算結果均方差為82.8;所提方法經過聲波信息融合和L-M算法改進后,能夠有效降低區(qū)域誤差不均勻的程度并提高雷電監(jiān)測區(qū)域邊緣的定位精度,誤差均值為108.4 m,目標監(jiān)測區(qū)域整體定位精度提高了46.6%,所算結果均方差為27.3.所提方法能夠對已有雷電定位系統(tǒng)的改進和聲波信息的融合利用提供參考.

    猜你喜歡
    子站雷聲聲波
    液壓平推CNG子站改造為標準CNG子站的建議
    愛的聲波 將愛留在她身邊
    中國寶玉石(2018年3期)2018-07-09 03:13:58
    聲波殺手
    淺談10kV配電自動化系統(tǒng)設計
    自適應BPSK在井下鉆柱聲波傳輸中的應用
    “聲波驅蚊”靠譜嗎
    配電自動化建設方案研究
    科學家(2016年4期)2016-07-25 03:07:01
    悲愴者
    文學港(2016年7期)2016-07-06 15:06:09
    夜晚的雷聲
    夜晚的雷聲
    一级,二级,三级黄色视频| 五月天丁香电影| 亚洲一码二码三码区别大吗| 日韩欧美三级三区| www.自偷自拍.com| 欧美日韩黄片免| 一级a爱视频在线免费观看| 99九九在线精品视频| 美女福利国产在线| 亚洲国产中文字幕在线视频| 欧美人与性动交α欧美精品济南到| 视频在线观看一区二区三区| 51午夜福利影视在线观看| 一区二区日韩欧美中文字幕| 波多野结衣av一区二区av| 天天躁夜夜躁狠狠躁躁| 精品少妇一区二区三区视频日本电影| 欧美日韩国产mv在线观看视频| 丁香六月欧美| 我的亚洲天堂| 国产免费视频播放在线视频| 91老司机精品| 最近最新免费中文字幕在线| 精品国产一区二区三区四区第35| 老熟女久久久| 日韩制服丝袜自拍偷拍| 在线观看人妻少妇| 少妇精品久久久久久久| 黑人巨大精品欧美一区二区蜜桃| 18禁观看日本| videos熟女内射| 国产男靠女视频免费网站| 亚洲 国产 在线| 欧美+亚洲+日韩+国产| 巨乳人妻的诱惑在线观看| 巨乳人妻的诱惑在线观看| 天天躁夜夜躁狠狠躁躁| 欧美日韩亚洲高清精品| 男人舔女人的私密视频| 亚洲国产成人一精品久久久| 国产精品熟女久久久久浪| 黄色视频在线播放观看不卡| 国产一区二区三区视频了| 日韩一区二区三区影片| 免费在线观看影片大全网站| 中文字幕制服av| 在线观看免费视频日本深夜| 精品国产国语对白av| 欧美性长视频在线观看| 老司机亚洲免费影院| 亚洲色图av天堂| 999精品在线视频| 建设人人有责人人尽责人人享有的| 老熟妇仑乱视频hdxx| 亚洲少妇的诱惑av| 99re在线观看精品视频| 久久狼人影院| 亚洲精品粉嫩美女一区| 性高湖久久久久久久久免费观看| 成在线人永久免费视频| 性色av乱码一区二区三区2| 国产精品 国内视频| 夜夜爽天天搞| 制服诱惑二区| 99久久国产精品久久久| 成人国语在线视频| 黄色丝袜av网址大全| 日日夜夜操网爽| 亚洲色图av天堂| 欧美 日韩 精品 国产| 少妇精品久久久久久久| 久久精品aⅴ一区二区三区四区| 亚洲av美国av| 99在线人妻在线中文字幕 | 亚洲欧美色中文字幕在线| 久久人人97超碰香蕉20202| 亚洲少妇的诱惑av| 久久中文看片网| 深夜精品福利| 国产精品.久久久| 69av精品久久久久久 | av不卡在线播放| 国产日韩欧美视频二区| 成人国产一区最新在线观看| 国产xxxxx性猛交| 一区二区日韩欧美中文字幕| 99re6热这里在线精品视频| 久久天堂一区二区三区四区| 国产在线一区二区三区精| 国产老妇伦熟女老妇高清| 免费在线观看视频国产中文字幕亚洲| 午夜两性在线视频| 悠悠久久av| 久久久国产欧美日韩av| 久久久精品免费免费高清| 亚洲精品美女久久久久99蜜臀| 欧美日本中文国产一区发布| 男女之事视频高清在线观看| 亚洲一区二区三区欧美精品| 别揉我奶头~嗯~啊~动态视频| 久久香蕉激情| 国产精品美女特级片免费视频播放器 | 亚洲精华国产精华精| av福利片在线| 午夜激情久久久久久久| 自拍欧美九色日韩亚洲蝌蚪91| av天堂久久9| 啦啦啦视频在线资源免费观看| 国产精品久久久久成人av| 国产精品久久久久久精品电影小说| 精品一区二区三区视频在线观看免费 | 黄网站色视频无遮挡免费观看| 天天躁夜夜躁狠狠躁躁| 国产一区有黄有色的免费视频| 国产不卡av网站在线观看| 国产伦理片在线播放av一区| 成年人免费黄色播放视频| 欧美激情极品国产一区二区三区| 亚洲av电影在线进入| 亚洲九九香蕉| 欧美日韩成人在线一区二区| 日本wwww免费看| 久久久久精品国产欧美久久久| 正在播放国产对白刺激| av国产精品久久久久影院| 国产精品久久久久成人av| 亚洲熟女毛片儿| 美女视频免费永久观看网站| 亚洲美女黄片视频| 一进一出好大好爽视频| 一个人免费在线观看的高清视频| 大陆偷拍与自拍| 中亚洲国语对白在线视频| 亚洲国产av新网站| 国产精品免费大片| 亚洲精品国产色婷婷电影| 国产精品一区二区免费欧美| 又黄又粗又硬又大视频| 久久久久久久大尺度免费视频| 精品一区二区三卡| 欧美日本中文国产一区发布| 国产亚洲精品久久久久5区| av超薄肉色丝袜交足视频| 99热网站在线观看| 高潮久久久久久久久久久不卡| 不卡av一区二区三区| 欧美性长视频在线观看| 黑人巨大精品欧美一区二区mp4| 18禁观看日本| 三上悠亚av全集在线观看| av超薄肉色丝袜交足视频| 国产av国产精品国产| a在线观看视频网站| 一本久久精品| 无限看片的www在线观看| 亚洲精品在线美女| 精品一品国产午夜福利视频| 一个人免费在线观看的高清视频| 高清毛片免费观看视频网站 | 午夜免费鲁丝| 成人永久免费在线观看视频 | 亚洲av国产av综合av卡| 成年女人毛片免费观看观看9 | 少妇的丰满在线观看| 国产av国产精品国产| 手机成人av网站| 女人高潮潮喷娇喘18禁视频| 国产精品一区二区精品视频观看| 狠狠精品人妻久久久久久综合| 亚洲,欧美精品.| 天堂中文最新版在线下载| 午夜老司机福利片| 亚洲精品国产色婷婷电影| 精品久久久精品久久久| 999久久久国产精品视频| 国产av一区二区精品久久| 少妇裸体淫交视频免费看高清 | aaaaa片日本免费| 男女高潮啪啪啪动态图| 中文字幕av电影在线播放| 久久国产精品影院| 亚洲精品在线观看二区| 18禁国产床啪视频网站| 纵有疾风起免费观看全集完整版| 性少妇av在线| 最近最新中文字幕大全电影3 | 亚洲熟妇熟女久久| 大型av网站在线播放| 久久久精品国产亚洲av高清涩受| 亚洲情色 制服丝袜| 超碰成人久久| 韩国精品一区二区三区| 亚洲精品一二三| 人成视频在线观看免费观看| 一本—道久久a久久精品蜜桃钙片| 国产99久久九九免费精品| 久久国产精品影院| 日本欧美视频一区| 女性生殖器流出的白浆| 18在线观看网站| 免费av中文字幕在线| 美国免费a级毛片| 妹子高潮喷水视频| 人妻 亚洲 视频| 精品免费久久久久久久清纯 | 男人操女人黄网站| 国产亚洲午夜精品一区二区久久| 99九九在线精品视频| 激情视频va一区二区三区| 美女午夜性视频免费| 国产有黄有色有爽视频| 国产欧美日韩综合在线一区二区| 丰满人妻熟妇乱又伦精品不卡| 精品国产国语对白av| 亚洲 国产 在线| 国产在线视频一区二区| 国产黄频视频在线观看| 国产精品亚洲av一区麻豆| 国产精品一区二区免费欧美| 欧美国产精品va在线观看不卡| 国产深夜福利视频在线观看| 精品少妇黑人巨大在线播放| 午夜福利在线观看吧| av网站在线播放免费| 一本综合久久免费| 亚洲伊人色综图| 亚洲视频免费观看视频| 久久免费观看电影| 自拍欧美九色日韩亚洲蝌蚪91| 黄色毛片三级朝国网站| 国产精品1区2区在线观看. | 两个人看的免费小视频| 欧美 亚洲 国产 日韩一| 人人妻,人人澡人人爽秒播| av有码第一页| 91国产中文字幕| 一级片'在线观看视频| 一级毛片精品| 国产视频一区二区在线看| 久久精品国产亚洲av香蕉五月 | 精品免费久久久久久久清纯 | 精品熟女少妇八av免费久了| 精品第一国产精品| 手机成人av网站| 十八禁高潮呻吟视频| 国产av一区二区精品久久| av网站免费在线观看视频| 熟女少妇亚洲综合色aaa.| 成人永久免费在线观看视频 | 成人黄色视频免费在线看| 人人澡人人妻人| 一区二区三区激情视频| 日韩欧美三级三区| 国产成人影院久久av| 欧美中文综合在线视频| 两性午夜刺激爽爽歪歪视频在线观看 | 麻豆国产av国片精品| 久久久久精品国产欧美久久久| 手机成人av网站| 国产在线观看jvid| 亚洲性夜色夜夜综合| 欧美人与性动交α欧美精品济南到| 国产一区二区 视频在线| 国产精品98久久久久久宅男小说| 精品一品国产午夜福利视频| 中文字幕另类日韩欧美亚洲嫩草| 久久这里只有精品19| 一本久久精品| 日韩制服丝袜自拍偷拍| 国产精品自产拍在线观看55亚洲 | 老司机午夜十八禁免费视频| 99国产精品99久久久久| 国产精品一区二区在线观看99| 别揉我奶头~嗯~啊~动态视频| 精品国产乱子伦一区二区三区| 999久久久精品免费观看国产| 国产欧美日韩一区二区三| tocl精华| 亚洲精品av麻豆狂野| 女人被躁到高潮嗷嗷叫费观| 国产在线免费精品| 美女主播在线视频| 黄色怎么调成土黄色| 久久婷婷成人综合色麻豆| 美女高潮喷水抽搐中文字幕| 99riav亚洲国产免费| 一级毛片女人18水好多| 纯流量卡能插随身wifi吗| 国产精品.久久久| 视频区图区小说| 天天躁日日躁夜夜躁夜夜| 丝袜喷水一区| 欧美精品一区二区大全| 国产精品欧美亚洲77777| 精品少妇一区二区三区视频日本电影| 人人妻人人爽人人添夜夜欢视频| 国产黄频视频在线观看| xxxhd国产人妻xxx| 老熟女久久久| 丝袜美腿诱惑在线| 高清av免费在线| 亚洲精品一卡2卡三卡4卡5卡| a级片在线免费高清观看视频| 欧美日韩视频精品一区| 中文字幕av电影在线播放| 91成人精品电影| 国产精品久久久久久人妻精品电影 | 一区二区三区国产精品乱码| 精品人妻1区二区| 国产一区二区三区在线臀色熟女 | 久久毛片免费看一区二区三区| 日本av免费视频播放| 午夜日韩欧美国产| 波多野结衣av一区二区av| av网站免费在线观看视频| 亚洲情色 制服丝袜| 国产深夜福利视频在线观看| 免费观看人在逋| 国产成人欧美在线观看 | 国产精品免费视频内射| 久久久欧美国产精品| 一区在线观看完整版| 国产成人精品久久二区二区免费| 精品人妻1区二区| 99国产精品一区二区蜜桃av | 精品国产乱子伦一区二区三区| 嫩草影视91久久| 国内毛片毛片毛片毛片毛片| 久久久国产一区二区| 亚洲中文字幕日韩| 欧美日本中文国产一区发布| 国产一区二区在线观看av| 国产精品99久久99久久久不卡| 日韩 欧美 亚洲 中文字幕| 久久亚洲真实| 国产在线精品亚洲第一网站| 老熟女久久久| 色94色欧美一区二区| 老司机午夜十八禁免费视频| 免费在线观看视频国产中文字幕亚洲| 日韩精品免费视频一区二区三区| a级毛片在线看网站| 黄色怎么调成土黄色| 免费黄频网站在线观看国产| 国产精品自产拍在线观看55亚洲 | 国产精品久久电影中文字幕 | 久久99热这里只频精品6学生| 亚洲专区中文字幕在线| 亚洲九九香蕉| 成人18禁在线播放| 久久久久久久大尺度免费视频| 国产一区二区三区在线臀色熟女 | 巨乳人妻的诱惑在线观看| a级毛片在线看网站| 日日爽夜夜爽网站| videos熟女内射| 亚洲七黄色美女视频| 两个人看的免费小视频| 极品少妇高潮喷水抽搐| 高清黄色对白视频在线免费看| 久久热在线av| 精品亚洲乱码少妇综合久久| 国产成人精品久久二区二区91| 黑人欧美特级aaaaaa片| 人人妻,人人澡人人爽秒播| 黄色视频,在线免费观看| 2018国产大陆天天弄谢| 婷婷成人精品国产| 精品国产乱码久久久久久男人| 国产片内射在线| 好男人电影高清在线观看| 51午夜福利影视在线观看| 久久ye,这里只有精品| 999久久久国产精品视频| 日韩人妻精品一区2区三区| 国产成人免费无遮挡视频| 大型黄色视频在线免费观看| 精品福利永久在线观看| 我的亚洲天堂| 丁香欧美五月| 欧美日韩视频精品一区| 精品高清国产在线一区| 欧美在线黄色| 亚洲天堂av无毛| 欧美国产精品va在线观看不卡| 精品少妇黑人巨大在线播放| 99热网站在线观看| 在线观看免费视频日本深夜| 天天添夜夜摸| 中文字幕最新亚洲高清| 久久精品国产a三级三级三级| 中国美女看黄片| 国产1区2区3区精品| avwww免费| 99热国产这里只有精品6| 国产精品熟女久久久久浪| 久久久久精品人妻al黑| 亚洲国产欧美在线一区| 动漫黄色视频在线观看| 久久中文字幕人妻熟女| 欧美精品av麻豆av| 十八禁高潮呻吟视频| 欧美精品一区二区免费开放| 欧美日韩国产mv在线观看视频| 两个人免费观看高清视频| 国产一区二区在线观看av| 99re6热这里在线精品视频| e午夜精品久久久久久久| 国产精品亚洲一级av第二区| 国产精品自产拍在线观看55亚洲 | 欧美精品av麻豆av| 人人妻,人人澡人人爽秒播| 中国美女看黄片| 国产亚洲精品一区二区www | a级片在线免费高清观看视频| 日本wwww免费看| 亚洲成人免费电影在线观看| 国产精品自产拍在线观看55亚洲 | 成人18禁高潮啪啪吃奶动态图| 精品人妻1区二区| 男人操女人黄网站| 曰老女人黄片| 90打野战视频偷拍视频| 国产精品 欧美亚洲| 日本五十路高清| 久久精品亚洲av国产电影网| 精品久久久精品久久久| 精品乱码久久久久久99久播| 欧美黄色淫秽网站| 丝袜在线中文字幕| 夜夜爽天天搞| 丰满少妇做爰视频| 精品人妻熟女毛片av久久网站| 精品福利永久在线观看| 日日夜夜操网爽| 精品一品国产午夜福利视频| 99久久精品国产亚洲精品| 一区二区三区激情视频| 99精国产麻豆久久婷婷| 天堂8中文在线网| 亚洲精品久久成人aⅴ小说| 日本撒尿小便嘘嘘汇集6| 精品国产超薄肉色丝袜足j| 欧美在线一区亚洲| 啦啦啦在线免费观看视频4| 国产欧美日韩一区二区三| 欧美一级毛片孕妇| 午夜精品国产一区二区电影| 99riav亚洲国产免费| 在线十欧美十亚洲十日本专区| 国产精品.久久久| 最新美女视频免费是黄的| 99re在线观看精品视频| 热99re8久久精品国产| 侵犯人妻中文字幕一二三四区| 香蕉丝袜av| 亚洲色图综合在线观看| 久久久精品免费免费高清| 777米奇影视久久| 黑人巨大精品欧美一区二区mp4| 中文字幕另类日韩欧美亚洲嫩草| 久久人人爽av亚洲精品天堂| 国产真人三级小视频在线观看| 极品人妻少妇av视频| 天堂俺去俺来也www色官网| 国产成人精品无人区| 精品久久久精品久久久| 国产在线视频一区二区| 午夜久久久在线观看| 在线观看人妻少妇| 欧美日韩一级在线毛片| 丰满饥渴人妻一区二区三| 国产精品 国内视频| 免费在线观看视频国产中文字幕亚洲| 日韩三级视频一区二区三区| 久久久久精品国产欧美久久久| 国产又色又爽无遮挡免费看| 日韩大片免费观看网站| 午夜视频精品福利| 亚洲一区二区三区欧美精品| 曰老女人黄片| 亚洲av第一区精品v没综合| 女人爽到高潮嗷嗷叫在线视频| 精品人妻熟女毛片av久久网站| 欧美中文综合在线视频| 国产野战对白在线观看| 国产一区二区在线观看av| 12—13女人毛片做爰片一| 嫁个100分男人电影在线观看| av免费在线观看网站| 亚洲七黄色美女视频| 女人爽到高潮嗷嗷叫在线视频| 成人亚洲精品一区在线观看| 黄色a级毛片大全视频| 久久毛片免费看一区二区三区| 日韩中文字幕视频在线看片| 日本五十路高清| 国产熟女午夜一区二区三区| 999久久久精品免费观看国产| 中亚洲国语对白在线视频| 757午夜福利合集在线观看| 国产免费现黄频在线看| 高潮久久久久久久久久久不卡| 欧美变态另类bdsm刘玥| 亚洲精品国产一区二区精华液| av在线播放免费不卡| 91麻豆av在线| 日韩视频在线欧美| 18禁裸乳无遮挡动漫免费视频| 午夜免费成人在线视频| 99香蕉大伊视频| 老司机靠b影院| 成人黄色视频免费在线看| 日韩一区二区三区影片| 纵有疾风起免费观看全集完整版| 免费在线观看视频国产中文字幕亚洲| 成在线人永久免费视频| 色播在线永久视频| 色94色欧美一区二区| 男男h啪啪无遮挡| 久久久国产精品麻豆| 免费在线观看日本一区| 十八禁网站网址无遮挡| 人人澡人人妻人| 免费在线观看完整版高清| 国产日韩一区二区三区精品不卡| 国产精品久久久人人做人人爽| 亚洲欧洲日产国产| 亚洲国产毛片av蜜桃av| 一本色道久久久久久精品综合| 色在线成人网| 少妇的丰满在线观看| 午夜福利视频在线观看免费| 女警被强在线播放| 欧美日韩亚洲高清精品| 视频区图区小说| 午夜免费成人在线视频| 菩萨蛮人人尽说江南好唐韦庄| 丰满饥渴人妻一区二区三| 超色免费av| 动漫黄色视频在线观看| tocl精华| 欧美亚洲日本最大视频资源| 美女主播在线视频| 亚洲 国产 在线| 99热国产这里只有精品6| 91字幕亚洲| 日韩一卡2卡3卡4卡2021年| 国产精品影院久久| 极品人妻少妇av视频| 午夜免费成人在线视频| 极品人妻少妇av视频| 日韩成人在线观看一区二区三区| 亚洲国产av影院在线观看| 纯流量卡能插随身wifi吗| 中文字幕人妻熟女乱码| 免费av中文字幕在线| 午夜福利欧美成人| 久久精品国产99精品国产亚洲性色 | 热99re8久久精品国产| videos熟女内射| 精品欧美一区二区三区在线| 久久精品熟女亚洲av麻豆精品| 国产精品免费大片| 久久久精品免费免费高清| 国产高清videossex| 国产不卡一卡二| 亚洲欧洲精品一区二区精品久久久| 精品乱码久久久久久99久播| 窝窝影院91人妻| 免费高清在线观看日韩| 丰满少妇做爰视频| 天堂中文最新版在线下载| 侵犯人妻中文字幕一二三四区| 成人手机av| 纯流量卡能插随身wifi吗| 99热网站在线观看| 色在线成人网| 高清av免费在线| 欧美精品高潮呻吟av久久| 操出白浆在线播放| 欧美成人午夜精品| 高清毛片免费观看视频网站 | 国产野战对白在线观看| 欧美国产精品一级二级三级| 亚洲成av片中文字幕在线观看| 王馨瑶露胸无遮挡在线观看| 亚洲av日韩在线播放| 美女国产高潮福利片在线看| 午夜福利在线观看吧| 黑人巨大精品欧美一区二区mp4| 久久午夜亚洲精品久久| 久久久久久久国产电影| 色婷婷av一区二区三区视频| 日韩视频在线欧美| 在线播放国产精品三级| 人人妻人人爽人人添夜夜欢视频| 91精品国产国语对白视频| 狠狠婷婷综合久久久久久88av| 国产精品久久久久久精品电影小说| 国产精品偷伦视频观看了| 欧美日韩亚洲高清精品| 操美女的视频在线观看| 国产又爽黄色视频| 欧美精品av麻豆av| 岛国毛片在线播放| av电影中文网址| 久久久久久久久久久久大奶| 欧美日韩黄片免| 性少妇av在线| 久久精品aⅴ一区二区三区四区| 黄色a级毛片大全视频| 久久久国产精品麻豆| 国产精品熟女久久久久浪| 亚洲av国产av综合av卡|