• <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
    夜晚的雷聲
    夜晚的雷聲
    国产女主播在线喷水免费视频网站 | 亚洲电影在线观看av| 亚洲成人精品中文字幕电影| 九九热线精品视视频播放| 久久天躁狠狠躁夜夜2o2o| 欧美色欧美亚洲另类二区| 天堂网av新在线| 欧美日韩综合久久久久久| 国产男靠女视频免费网站| 国产一区亚洲一区在线观看| 国产片特级美女逼逼视频| 精品人妻一区二区三区麻豆 | 国产精品一区二区免费欧美| 国产乱人视频| 一区二区三区四区激情视频 | 岛国在线免费视频观看| 久久精品夜色国产| 日本一二三区视频观看| 欧美国产日韩亚洲一区| 久久鲁丝午夜福利片| 亚洲成人久久性| 色综合亚洲欧美另类图片| 综合色丁香网| 亚洲人与动物交配视频| 亚洲av中文字字幕乱码综合| 插逼视频在线观看| 久久精品国产亚洲av香蕉五月| 免费大片18禁| 国产69精品久久久久777片| 日本一本二区三区精品| 熟女电影av网| 99热只有精品国产| 午夜免费激情av| 成年av动漫网址| 亚洲成av人片在线播放无| 免费在线观看影片大全网站| 国产高清视频在线观看网站| 国产精品嫩草影院av在线观看| 别揉我奶头 嗯啊视频| 国产色婷婷99| 两性午夜刺激爽爽歪歪视频在线观看| 免费在线观看影片大全网站| 无遮挡黄片免费观看| 久久精品影院6| 高清毛片免费观看视频网站| 国产在线精品亚洲第一网站| 亚洲电影在线观看av| 搡女人真爽免费视频火全软件 | 精品国内亚洲2022精品成人| 永久网站在线| 天堂√8在线中文| 亚洲精品国产av成人精品 | 亚洲内射少妇av| 色综合站精品国产| 亚洲成人av在线免费| 国产精品一区www在线观看| 亚洲激情五月婷婷啪啪| 老司机午夜福利在线观看视频| www.色视频.com| 热99re8久久精品国产| 九九热线精品视视频播放| 91久久精品电影网| 老熟妇乱子伦视频在线观看| 国产aⅴ精品一区二区三区波| 男插女下体视频免费在线播放| 免费高清视频大片| 欧美日韩精品成人综合77777| av.在线天堂| 久久久久久久午夜电影| 99在线人妻在线中文字幕| 男女那种视频在线观看| 韩国av在线不卡| 在线免费十八禁| а√天堂www在线а√下载| 国产精品人妻久久久影院| 18禁黄网站禁片免费观看直播| 欧美一区二区精品小视频在线| 在线天堂最新版资源| 婷婷亚洲欧美| 国产白丝娇喘喷水9色精品| 国产午夜精品论理片| 亚洲av中文字字幕乱码综合| 99久久精品一区二区三区| 又黄又爽又免费观看的视频| 国产精品亚洲一级av第二区| 一区二区三区高清视频在线| 日韩人妻高清精品专区| 美女大奶头视频| 天堂√8在线中文| 欧美日韩乱码在线| 男人的好看免费观看在线视频| 狂野欧美白嫩少妇大欣赏| 乱码一卡2卡4卡精品| 免费av观看视频| 亚洲专区国产一区二区| 亚洲av中文av极速乱| 久久人人精品亚洲av| 国产精品一及| 精品久久国产蜜桃| 少妇裸体淫交视频免费看高清| 午夜福利视频1000在线观看| 国产中年淑女户外野战色| 1024手机看黄色片| 亚洲欧美清纯卡通| 一区二区三区免费毛片| 亚洲国产精品成人综合色| 国产成人影院久久av| 51国产日韩欧美| 久久久色成人| 欧美中文日本在线观看视频| 男女那种视频在线观看| 久久久久国产网址| 悠悠久久av| 免费av毛片视频| 俄罗斯特黄特色一大片| 久久韩国三级中文字幕| a级毛片免费高清观看在线播放| 丰满乱子伦码专区| 俄罗斯特黄特色一大片| 99久久无色码亚洲精品果冻| 美女cb高潮喷水在线观看| 中文字幕精品亚洲无线码一区| 免费看av在线观看网站| 国产精品一区二区免费欧美| 18+在线观看网站| 成年版毛片免费区| 亚洲人成网站在线播放欧美日韩| 欧美不卡视频在线免费观看| 久久久久国产网址| 国产成人freesex在线 | 国产三级中文精品| 成人三级黄色视频| 菩萨蛮人人尽说江南好唐韦庄 | 国产美女午夜福利| 中文字幕免费在线视频6| 99久久成人亚洲精品观看| 成人欧美大片| 少妇高潮的动态图| 日本免费一区二区三区高清不卡| 搡女人真爽免费视频火全软件 | 尤物成人国产欧美一区二区三区| 久久精品人妻少妇| 国产高清有码在线观看视频| 午夜精品一区二区三区免费看| 男女啪啪激烈高潮av片| av视频在线观看入口| 三级毛片av免费| 欧美另类亚洲清纯唯美| 国产精品一区二区免费欧美| 成人高潮视频无遮挡免费网站| 永久网站在线| 久久精品夜夜夜夜夜久久蜜豆| 亚洲精品国产成人久久av| 国产精品女同一区二区软件| 小说图片视频综合网站| 三级男女做爰猛烈吃奶摸视频| 国产精品福利在线免费观看| 国产男人的电影天堂91| 亚洲av第一区精品v没综合| 1024手机看黄色片| 蜜桃久久精品国产亚洲av| 国产色爽女视频免费观看| 在线国产一区二区在线| 国产亚洲精品综合一区在线观看| 久久久久国产精品人妻aⅴ院| ponron亚洲| 菩萨蛮人人尽说江南好唐韦庄 | 精品久久久噜噜| 午夜福利在线在线| 18禁裸乳无遮挡免费网站照片| 可以在线观看毛片的网站| 午夜影院日韩av| 欧美zozozo另类| 中文字幕人妻熟人妻熟丝袜美| 国产午夜福利久久久久久| 色综合站精品国产| 在线播放国产精品三级| 日韩精品青青久久久久久| 两个人视频免费观看高清| 中文字幕人妻熟人妻熟丝袜美| 天天躁夜夜躁狠狠久久av| 69人妻影院| 亚洲久久久久久中文字幕| 日本-黄色视频高清免费观看| 狂野欧美激情性xxxx在线观看| 校园春色视频在线观看| 日产精品乱码卡一卡2卡三| 精品免费久久久久久久清纯| 真实男女啪啪啪动态图| 国产一区二区在线观看日韩| 午夜激情福利司机影院| 精品福利观看| 欧美性猛交╳xxx乱大交人| 精品99又大又爽又粗少妇毛片| 国产午夜福利久久久久久| 欧美性猛交╳xxx乱大交人| 亚洲图色成人| 国产男靠女视频免费网站| 午夜老司机福利剧场| 国产伦精品一区二区三区四那| 日韩中字成人| 久久久色成人| 日韩,欧美,国产一区二区三区 | 男女下面进入的视频免费午夜| 久久午夜亚洲精品久久| 三级经典国产精品| 久久综合国产亚洲精品| a级一级毛片免费在线观看| 久久精品夜夜夜夜夜久久蜜豆| 午夜爱爱视频在线播放| 全区人妻精品视频| 丝袜美腿在线中文| 亚洲性久久影院| 99久久中文字幕三级久久日本| 在线看三级毛片| 亚洲中文字幕日韩| 亚洲国产精品成人综合色| 天堂av国产一区二区熟女人妻| 国产高清激情床上av| 欧美xxxx性猛交bbbb| 三级经典国产精品| 18禁在线播放成人免费| 啦啦啦韩国在线观看视频| 国产一区二区在线观看日韩| 国产亚洲欧美98| 亚洲欧美精品自产自拍| videossex国产| 国产亚洲精品av在线| 丰满人妻一区二区三区视频av| 97超级碰碰碰精品色视频在线观看| 欧美区成人在线视频| 成人漫画全彩无遮挡| 精品久久久久久久久久久久久| 久久久午夜欧美精品| 99riav亚洲国产免费| 国产成人一区二区在线| 久久久久久国产a免费观看| 久久久久久久久久久丰满| 亚洲精品一区av在线观看| 国产精华一区二区三区| h日本视频在线播放| 久久久精品欧美日韩精品| 十八禁国产超污无遮挡网站| 亚洲精品日韩av片在线观看| 免费在线观看成人毛片| 成人特级黄色片久久久久久久| 激情 狠狠 欧美| 老司机福利观看| 人人妻人人看人人澡| 国产69精品久久久久777片| 啦啦啦韩国在线观看视频| 男人舔奶头视频| 九九热线精品视视频播放| 国产黄色视频一区二区在线观看 | 色综合站精品国产| 干丝袜人妻中文字幕| 大香蕉久久网| АⅤ资源中文在线天堂| 久久久久久久久久久丰满| 午夜影院日韩av| 麻豆成人午夜福利视频| 亚洲国产高清在线一区二区三| av在线播放精品| 看片在线看免费视频| 国产成人91sexporn| 成年免费大片在线观看| 神马国产精品三级电影在线观看| 中文字幕av成人在线电影| 国产精品无大码| 日韩欧美免费精品| 午夜福利成人在线免费观看| 亚洲精品一卡2卡三卡4卡5卡| 国产不卡一卡二| 午夜老司机福利剧场| 久久久精品大字幕| 大型黄色视频在线免费观看| 噜噜噜噜噜久久久久久91| 国产精品久久久久久久久免| 国产精品1区2区在线观看.| 五月玫瑰六月丁香| 亚洲一区高清亚洲精品| 日韩精品中文字幕看吧| avwww免费| 性插视频无遮挡在线免费观看| 男人狂女人下面高潮的视频| 国产精品99久久久久久久久| 欧美日本亚洲视频在线播放| 91久久精品国产一区二区成人| 成人亚洲精品av一区二区| 亚洲经典国产精华液单| 久久精品国产亚洲av香蕉五月| 欧美成人一区二区免费高清观看| 国产色婷婷99| 精品人妻偷拍中文字幕| 成人特级av手机在线观看| 夜夜夜夜夜久久久久| 国内精品一区二区在线观看| 日韩制服骚丝袜av| 午夜爱爱视频在线播放| 一个人看视频在线观看www免费| 最新中文字幕久久久久| 毛片一级片免费看久久久久| 成熟少妇高潮喷水视频| 亚洲,欧美,日韩| 日本与韩国留学比较| 一夜夜www| 亚洲国产精品国产精品| 国产高清视频在线播放一区| 在线免费观看的www视频| 日本一本二区三区精品| 97碰自拍视频| 欧美在线一区亚洲| 91精品国产九色| 黄色一级大片看看| 搡女人真爽免费视频火全软件 | 欧美色视频一区免费| 啦啦啦观看免费观看视频高清| 成人毛片a级毛片在线播放| 亚洲精品粉嫩美女一区| 最近2019中文字幕mv第一页| 赤兔流量卡办理| 亚洲真实伦在线观看| 欧美成人a在线观看| 大又大粗又爽又黄少妇毛片口| 亚洲激情五月婷婷啪啪| 嫩草影院精品99| 狂野欧美白嫩少妇大欣赏| 在线国产一区二区在线| av在线观看视频网站免费| 欧美激情在线99| 成人欧美大片| ponron亚洲| 身体一侧抽搐| 床上黄色一级片| 精品无人区乱码1区二区| 亚洲经典国产精华液单| 久久久久久国产a免费观看| 99久久成人亚洲精品观看| 自拍偷自拍亚洲精品老妇| 午夜福利视频1000在线观看| 天天一区二区日本电影三级| 十八禁网站免费在线| 中文字幕免费在线视频6| 国产高潮美女av| 欧美性猛交黑人性爽| 成年女人毛片免费观看观看9| 亚洲中文字幕日韩| 最近视频中文字幕2019在线8| 国产在线精品亚洲第一网站| 国产 一区精品| 精品欧美国产一区二区三| 你懂的网址亚洲精品在线观看 | 免费不卡的大黄色大毛片视频在线观看 | 久久久久性生活片| av中文乱码字幕在线| 午夜激情欧美在线| 亚洲乱码一区二区免费版| 国产精品一二三区在线看| 中文字幕av成人在线电影| 熟女人妻精品中文字幕| 波多野结衣高清无吗| 午夜激情福利司机影院| 久久人人精品亚洲av| 日本黄色片子视频| 一级a爱片免费观看的视频| 成年女人毛片免费观看观看9| 国产 一区精品| 91精品国产九色| 五月伊人婷婷丁香| 久久精品国产亚洲av天美| 草草在线视频免费看| 校园春色视频在线观看| 国产精品久久久久久精品电影| 尤物成人国产欧美一区二区三区| 97超视频在线观看视频| 大又大粗又爽又黄少妇毛片口| АⅤ资源中文在线天堂| 久久久久久久久中文| 国产一区二区三区在线臀色熟女| 男女那种视频在线观看| 五月伊人婷婷丁香| 少妇被粗大猛烈的视频| 97超视频在线观看视频| 秋霞在线观看毛片| 又黄又爽又免费观看的视频| 日韩欧美三级三区| 性欧美人与动物交配| 日本黄色视频三级网站网址| 搡老妇女老女人老熟妇| 欧美最黄视频在线播放免费| 男女下面进入的视频免费午夜| 免费黄网站久久成人精品| 国产精品av视频在线免费观看| 久久久久久久久久久丰满| 亚洲av.av天堂| 十八禁国产超污无遮挡网站| 女生性感内裤真人,穿戴方法视频| 免费av观看视频| 日韩强制内射视频| 两性午夜刺激爽爽歪歪视频在线观看| 免费搜索国产男女视频| 熟女电影av网| 人妻少妇偷人精品九色| 免费看a级黄色片| 国产单亲对白刺激| 日本 av在线| 中文字幕人妻熟人妻熟丝袜美| 亚洲av美国av| 亚洲av中文字字幕乱码综合| 少妇裸体淫交视频免费看高清| 日韩,欧美,国产一区二区三区 | 国产精品免费一区二区三区在线| 大又大粗又爽又黄少妇毛片口| 国产精品三级大全| 国产精品不卡视频一区二区| 日本a在线网址| 婷婷亚洲欧美| 久久国内精品自在自线图片| 男女啪啪激烈高潮av片| 成年女人毛片免费观看观看9| 少妇丰满av| 在线播放国产精品三级| 国产av不卡久久| 真人做人爱边吃奶动态| 搡老妇女老女人老熟妇| 小说图片视频综合网站| 别揉我奶头~嗯~啊~动态视频| 亚洲精品日韩在线中文字幕 | 天天一区二区日本电影三级| 97热精品久久久久久| 不卡一级毛片| a级毛色黄片| 五月伊人婷婷丁香| 最近最新中文字幕大全电影3| 日本精品一区二区三区蜜桃| 亚洲中文日韩欧美视频| 亚洲欧美日韩卡通动漫| 亚洲电影在线观看av| 能在线免费观看的黄片| 少妇人妻一区二区三区视频| 亚洲国产精品国产精品| 国产黄a三级三级三级人| 亚洲国产日韩欧美精品在线观看| 色播亚洲综合网| 欧美日韩一区二区视频在线观看视频在线 | 99国产精品一区二区蜜桃av| 国产视频一区二区在线看| 国产欧美日韩一区二区精品| 日韩人妻高清精品专区| 国产中年淑女户外野战色| 深夜a级毛片| 国内久久婷婷六月综合欲色啪| 日韩一区二区视频免费看| av在线天堂中文字幕| 欧美激情久久久久久爽电影| 99riav亚洲国产免费| 成人毛片a级毛片在线播放| 亚洲成a人片在线一区二区| 日本色播在线视频| 久久人妻av系列| 国产黄a三级三级三级人| 亚洲人成网站在线观看播放| 日日摸夜夜添夜夜爱| 成人特级黄色片久久久久久久| 国产精品美女特级片免费视频播放器| 深夜a级毛片| 久久久久精品国产欧美久久久| 国产精品一区二区三区四区久久| 日韩成人伦理影院| 欧美成人a在线观看| 国产精品一区二区性色av| 97超级碰碰碰精品色视频在线观看| 午夜a级毛片| 又黄又爽又刺激的免费视频.| 久久久久久久午夜电影| 国产精品国产高清国产av| 禁无遮挡网站| 2021天堂中文幕一二区在线观| 亚洲av电影不卡..在线观看| 波多野结衣高清作品| 国产亚洲精品综合一区在线观看| 性插视频无遮挡在线免费观看| 国内精品美女久久久久久| 99热全是精品| 国产伦精品一区二区三区视频9| 青春草视频在线免费观看| 欧美性感艳星| 小说图片视频综合网站| 非洲黑人性xxxx精品又粗又长| 18禁黄网站禁片免费观看直播| 久久人妻av系列| or卡值多少钱| 亚洲精品一卡2卡三卡4卡5卡| 亚洲国产精品成人久久小说 | 91在线精品国自产拍蜜月| 插逼视频在线观看| 久久精品国产亚洲网站| 亚洲欧美成人综合另类久久久 | 久久久a久久爽久久v久久| 综合色av麻豆| av天堂中文字幕网| 精品人妻熟女av久视频| 午夜a级毛片| 国产精品电影一区二区三区| 亚洲欧美日韩卡通动漫| 12—13女人毛片做爰片一| 成年女人毛片免费观看观看9| 在线免费观看的www视频| 欧美日本视频| 一边摸一边抽搐一进一小说| 日日啪夜夜撸| 中文字幕人妻熟人妻熟丝袜美| a级毛片a级免费在线| 亚洲经典国产精华液单| 热99re8久久精品国产| 日韩精品有码人妻一区| 日日摸夜夜添夜夜添小说| 99国产精品一区二区蜜桃av| 一进一出抽搐gif免费好疼| 国产精品伦人一区二区| 国产成年人精品一区二区| 中国美白少妇内射xxxbb| 97热精品久久久久久| 男人舔奶头视频| 国内久久婷婷六月综合欲色啪| 国产精品日韩av在线免费观看| 老熟妇乱子伦视频在线观看| 最近2019中文字幕mv第一页| 久久婷婷人人爽人人干人人爱| 午夜视频国产福利| 国产精品久久视频播放| 久久久午夜欧美精品| 久久久久精品国产欧美久久久| 噜噜噜噜噜久久久久久91| 日本成人三级电影网站| 秋霞在线观看毛片| 我的女老师完整版在线观看| 欧美日本视频| 国产亚洲av嫩草精品影院| 深夜a级毛片| 黄色欧美视频在线观看| av在线播放精品| 成人av一区二区三区在线看| 亚洲久久久久久中文字幕| 欧美一区二区国产精品久久精品| 国语自产精品视频在线第100页| 精品人妻熟女av久视频| 国产伦精品一区二区三区视频9| 看黄色毛片网站| 亚洲精品久久国产高清桃花| 免费看a级黄色片| 男女那种视频在线观看| 最近中文字幕高清免费大全6| 国产免费一级a男人的天堂| 晚上一个人看的免费电影| 九色成人免费人妻av| 国产精品一区二区三区四区免费观看 | 欧美不卡视频在线免费观看| 国产成人aa在线观看| 婷婷色综合大香蕉| 国产av麻豆久久久久久久| 亚洲五月天丁香| av在线观看视频网站免费| 国产在线精品亚洲第一网站| 男女之事视频高清在线观看| 色哟哟哟哟哟哟| 欧美xxxx性猛交bbbb| 亚洲国产精品sss在线观看| 女的被弄到高潮叫床怎么办| 日韩精品有码人妻一区| 日韩欧美国产在线观看| 亚洲av免费在线观看| 国产激情偷乱视频一区二区| 久久精品综合一区二区三区| 亚洲丝袜综合中文字幕| 国产亚洲91精品色在线| 欧美最黄视频在线播放免费| 激情 狠狠 欧美| 看免费成人av毛片| 性色avwww在线观看| 两个人的视频大全免费| 日本撒尿小便嘘嘘汇集6| av国产免费在线观看| 日本在线视频免费播放| 国产亚洲91精品色在线| 日韩成人伦理影院| 18禁在线播放成人免费| 伦理电影大哥的女人| 俺也久久电影网| 一区二区三区高清视频在线| 久久99热6这里只有精品| 国内揄拍国产精品人妻在线| 亚洲在线自拍视频| 久久久久久国产a免费观看| 熟妇人妻久久中文字幕3abv| 天天躁日日操中文字幕| 欧美三级亚洲精品| 久久精品国产亚洲网站| 久久亚洲精品不卡| 又黄又爽又刺激的免费视频.| 久久久久久国产a免费观看| 婷婷亚洲欧美| 麻豆精品久久久久久蜜桃| 亚州av有码| 国产白丝娇喘喷水9色精品| 亚洲成a人片在线一区二区| 天堂√8在线中文| 国产黄片美女视频| 啦啦啦啦在线视频资源| 精品久久国产蜜桃| 国产老妇女一区| 亚洲欧美精品综合久久99| 欧美性猛交黑人性爽| 亚洲中文字幕一区二区三区有码在线看| 能在线免费观看的黄片|