999精品在线视频,手机成人午夜在线视频,久久不卡国产精品无码,中日无码在线观看,成人av手机在线观看,日韩精品亚洲一区中文字幕,亚洲av无码人妻,四虎国产在线观看 ?

反潛直升機吊放聲納水面艦艇編隊尾后鋸齒形法搜潛

2017-11-23 05:57:57羅木生曾家有侯學隆
航空學報 2017年1期

羅木生,曾家有,侯學隆

反潛直升機吊放聲納水面艦艇編隊尾后鋸齒形法搜潛

羅木生*,曾家有,侯學隆

海軍航空工程學院 指揮系,煙臺 264001

針對水面艦艇編隊尾后方向的潛艇威脅,研究反潛直升機使用吊放聲納搜潛中的警戒扇面、兵力需求、參數決策等問題。通過分析鋸齒形法特點、吊放聲納搜潛過程,建立了多機艦艇編隊尾后搜潛的警戒扇面定量計算模型。研究了齒形法相鄰探測點間距系數、搜索速度等制約因素,探索了吊放聲納作用距離與潛艇航速的定量關系,給出了已探測區等效圓半徑的計算方法。仿真計算了艦艇編隊尾后搜潛中實際對潛警戒扇面、最優間距系數、不同任務要求下的兵力需求。結果表明:配置2架直升機于艦艇編隊尾后15海里處,可建立80.9°的對潛警戒扇面;齒形法最佳間距系數介于1.1~1.3之間。

艦艇編隊尾后搜潛;鋸齒形法;對潛警戒扇面;吊放聲納;反潛直升機;尾隨跟蹤潛艇

水面艦艇以編隊形式遂行海上作戰任務已漸成常態。現代潛艇以其優良的靜音性、隱蔽性,對水面艦艇編隊(如航母編隊)構成的威脅日趨嚴峻。尤其是推進技術日新月異,潛艇水下航速較低、機動能力弱的不足得到較大改善,使潛艇采用追擊戰術攻擊水面艦艇編隊成為可能。因而,水面艦艇編隊不僅要抗擊來自前方、側翼的水下威脅,還需防御來自尾后的水下威脅。實施尾后搜潛已成為編隊反潛作戰的重要方面[1]。

搜潛相關研究較豐富[2-5],然針對編隊尾后方向的搜潛問題,目前的研究集中于水面艦艇使用拖曳線列陣聲納[6]、反潛直升機使用浮標聲納搜潛等方面[7]。拖曳線列陣聲納使用時對航速、航向有要求,制約編隊機動;浮標聲納[8-10]是消耗品,難以滿足編隊航渡中多次布設需求。吊放聲納[11-13]可重復使用,是反潛直升機特有的搜潛設備。因此,本文對反潛直升機配置在水面艦艇編隊尾后,使用吊放聲納搜潛展開研究,以期較好解決作戰決策與實施過程中存在的關鍵問題。

1 水面艦艇編隊尾后搜潛特點

水面艦艇編隊尾后搜潛,是指將反潛兵力配置于艦艇編隊尾后適當位置,建立攔截線,實施對潛搜索的活動方法,以發現尾隨跟蹤或采用追擊戰術的敵方潛艇。如圖1所示。與水面艦艇相比,反潛直升機遂行尾后搜潛任務,可充分發揮其飛行速度快、機動性好、隱蔽性優、不易遭敵潛艇攻擊等優點,從而實現與編隊其他兵力一道,快速構建起360°的對潛搜索防線。實施編隊尾后搜潛,尤其是編隊航渡過程中,必須能夠較長時間、較長距離上對編隊尾后方向的水下目標實施搜索。反潛直升機雖可以使用反潛雷達、浮標聲納、吊放聲納等諸多器材實施搜索。但反潛雷達難以發現水下的潛艇;而浮標聲納屬于一次性消耗品,且投布之后位置固定,編隊航渡中尾后搜索需投布的數量大,難以滿足尾后持續搜索需求。吊放聲納可重復使用、通過直升機機動可到最佳位置進行探測,且作用距離大、搜潛效率高,能較好地滿足編隊尾后搜潛的需要[14-15]。

2 吊放聲納艦艇編隊尾后搜潛方法

由于艦艇編隊處于航行中,動態變化的特點對反潛直升機吊放聲納的搜索方法提出了嚴格要求。

1)搜索速度要同步、可調。為確保對艦艇編隊尾后一定扇面實施警戒,應使得搜索陣位相對編隊保持穩定,這就要求搜索速度與編隊航行速度同步。另一方面,編隊航行速度可能因作戰需要而調整,這就要求吊放聲納的搜索速度能夠適應性調整。

2)搜索寬度要可控。搜索寬度即攔截線長度,決定了搜索扇面的大小。戰場態勢的變化將導致尾后搜索扇面需求發生變化,因而要求搜索寬度能夠靈活控制。

在吊放聲納諸多搜索方法中,鋸齒形法最能滿足上述要求。鋸齒形法[16],又稱曲折機動搜索法,指反潛直升機使用吊放聲納沿一個主航向按曲折航線同步變換探測點探測潛艇的一種方法,如圖2所示。圖中,dhel,jg為相鄰反潛直升機之間的間距,djg為相鄰探測點的間距。反潛直升機通常采取S形航線從當前懸停點逆風爬升飛行到下一懸停點,逆風懸停后使用吊放聲納搜索。但在一個搜索周期 Tdipson,cyc內[17],直升機運動的位移為相鄰探測點的間距djg,將其投影到主航向上則為djgcosα。因此,鋸齒形法主航向平均搜索速度Vhel,search為

式中:Rdipson為吊放聲納有效作用距離;ξdipson為間距系數;α為反潛直升機飛行方向與主航向(即水面艦艇編隊航向)之間的夾角。

由式(1)可知,α、ξdipson等參數可以調整搜索速度從而實現與編隊航速同步;也可結合dhel,jg,對搜索寬度進行控制。此外,鋸齒形法還提高了搜索的隱蔽性,即使吊放聲納以主動方式工作,潛艇也難以判斷直升機搜索的主航向。

3 吊放聲納艦艇編隊尾后鋸齒形法搜潛建模

鋸齒形法搜索區域覆蓋尾后扇面大小是編隊尾后搜潛決策的關鍵問題之一。由于艦艇編隊以一定速度航行,故尾后搜潛區域相對編隊的距離、方位處于動態變化中。而且,直升機使用吊放聲納搜潛[18-19]屬于離散式的搜索方法,使得尾后搜索區域覆蓋艦艇編隊尾后扇面的大小隨時間不斷變化,故需建立動態求解方法。

首先建立直角坐標系:以編隊航向為X軸,以過吊放聲納首個探測點、垂直X軸方向為Y軸,兩軸交點為坐標原點,如圖3(a)所示。設出動Nhel架反潛直升機采用鋸齒形法實施尾后搜潛,首次使用吊放聲納開始聽測時刻為0時刻,D0為此時艦艇編隊中心位置的橫坐標,則第k架直升機第1個探測點的坐標

為使得搜索區域均勻分布在編隊航線兩側,則第1架直升機第1個探測點的坐標O1,1(x1,1,y1,1)需滿足:

根據鋸齒形法的特點,結合圖3(b),可得第k架直升機第i (i>1)個探測點的坐標Ok,i(xk,i,yk,i)為

令Ok,i與編隊中心的距離為dk,i、相對編隊中心的舷角為ψk,i,cent;rk,i()t為第k架直升機第i個探測點不存在潛艇目標區域的等效圓半徑,βk,i為此區域覆蓋艦艇編隊尾后扇面的圓心角,則有

式中:dk,i(t)=槡(xk,i(t)-xship(t ))2+( Δyk,i(t))2,Δyk,i(t)=yk,i-yship(t);(xship,yship)為艦艇編隊中心位置的坐標,令Vship為水面艦艇編隊航速,則有

為便于計算,可取

若ψk,i,cent>π,表示Ok,i位于右舷,取ψk,i,cent(t)-2π即為實際舷角。

令ψk,i,start()t、ψk,i,end()t為rk,i()t覆蓋編隊尾后扇面的起始、終止舷角,則有

根據式(7),可計算出任意時刻t每架反潛直升機、每個探測點未發現目標區域覆蓋編隊尾后扇面的起始、終止舷角;若扇面之間存在交疊,則合并為一個扇面;從中選取圓心角最大的扇面,作為該時刻反潛直升機搜索區域覆蓋尾后的扇面。

1)rk,i(t)的解算

rk,i(t)與吊放聲納的使用密切相關。在每一個周期內,吊放聲納真正探測潛艇目標的時間僅為換能器開機聽測 時 間Tdipson,listen[19]。設第k 架直升機在第i個探測點使用吊放聲納開始聽測與結束聽測的時刻分別為tk,i,start、tk,i,end,則有

在聽測期間[tk,i,start,tk,i,end],吊放聲納可以發現其有效作用距離Rdipson內的水下目標。由于潛艇目標處于不斷運動中,也就意味著第i個探測點不存在目標的區域,在聽測結束后會隨著時間的推移不斷縮小直至為0[19-20],如圖3(b)所示。rk,i()t縮小的速率取決于目標可能的運動速度。令珔vsub為目標運動速度的均值,有

2)ξdipson的解算

鋸齒形法搜索時,最佳的間距可獲得較高的搜索效率和發現概率。式(1)變換后有

將djg=ξdipsonRdipson代入式(10)變換得

由于cosα≤1,則

若α=0,即鋸齒形法變為平行搜索法,有

3)搜索時間的解算

搜索時間為探測點數量與周期Tdipson,cyc的乘積。設單架反潛直升機完成第nsearch個探測點的搜索后準備返航時,與編隊的距離為Dback,則

由于反潛直升機向尾后搜潛陣位航渡、實施搜潛、返回載艦的時間之和應不超過反潛直升機留空時間。故有

式中:Thel,tot為反潛直升機留空時間;Vhel,fly為反潛直升機飛行速度;Tdipson,cyc為相鄰探測點之間過渡飛行、收放吊放聲納與聽測三者時間之和。

將式(14)代入式(15),變化后可得單架反潛直升機一次出動最多可搜索的探測點數量nsearch,即

4)潛艇航速及其影響分析

潛艇為了實現對水面艦艇編隊的尾隨跟蹤,其平均航速v珔sub必須與艦艇編隊航速Vship相近。但當潛艇航速超過最大靜音航速時,其產生的噪聲急劇增加,容易在更遠的距離上被吊放聲納發現。

由于潛艇航速vsub與吊放聲納作用距離Rdipson之間的關系極為復雜,應根據實際試驗測量得出。為簡化分析,假定當vsub∈ [0 ,vsub,two],則

式中:vsub,one、vsub,two分別為潛艇的靜音航速、最大靜音航速;Rdipson,0為吊放聲納對以vsub,one航速航行的潛艇的有效探測距離。

4 仿真與分析

設水面艦艇編隊以20節(1節=1.852km/h)速度航行;3架反潛直升機可用,其留空時間Thel,tot為4h、飛行速度200km/h;吊放聲納對以6節航速航行的潛艇的有效探測距離Rdipson,0為6km,搜索周期Tdipson,cyc為11min,在每個探測點的聽測時間Tdipson,listen為6min;潛艇最大航速28節,靜音航速vsub,one為6節、最大靜音航速vsub,two為12節;潛艇為實現對水面艦艇編隊的持續跟蹤,其平均航速與編隊航速相同。

4.1 尾后搜潛陣位相對艦艇編隊固定

為確保尾后搜潛陣位相對編隊的固定,則要求Vhel,search=Vship。為取得最佳效能,需對兵力數量、間距系數ξdipson、D0等參數進行決策。

1)兵力數量與間距系數ξdipson的決策

若任務要求在艦艇編隊尾后20海里(1海里=1.852km)處建立攔截線,且保持陣位與編隊的距離,即D0取20海里。對于不同的ξdipson取值進行仿真,結果如圖4和表1所示。

由圖4(a)可知,兵力數量的增加可擴大尾后搜索扇面,由表1知,當出動3架反潛直升機時,至少可對尾后89.2°的扇面進行搜索。存在最佳間距系數,使得尾后最大扇面的中心角最大,最佳值如表1所示;但超過最佳值后急劇變小,這是因為已探測區不存在潛艇的區域衰減,使得原本重疊的已探測區分離,導致尾后最大扇面急劇變小。

表1 D0不變時尾后搜潛扇面最佳參數Table 1 Optimal parameters searching stern of ship formation with same D0

圖4(b)為最佳間距時,不同兵力數量下尾后最大扇面的中心角隨時間的變化情況。在第1~3個探測點搜索期間,由于已探測區數量少,未形成較為穩定的尾后扇面,但隨后尾后扇面中心角在一定區間內變化不大。

根據表1,可根據任務要求,確定合適的兵力出動數量。例如:要求對尾后65°扇面進行搜索,則出動2架反潛直升機即可完成任務。

2)D0的決策

不同的D0取值,所能覆蓋尾后扇面的大小也不同;實際作戰中指揮員應根據編隊反潛體系構建需要,確定合適的尾后搜潛陣位與艦艇編隊中心位置的距離D0。以出動2架反潛直升機為例,當D0取不同值時,仿真結果如圖5所示。由圖5可知,隨著D0的增加,尾后最大扇面的中心角迅速減少;間距系數ξdipson雖有所增加,但變化不大,均在1.14~1.20之間。例如:若任務要求出動2架直升機對尾后80°扇面進行搜索,則取D0為15海里、取ξdipson為1.15,可建立尾后80.9°的對潛警戒扇面。實際作戰中,應根據任務要求、敵潛艇可能出現的尾后陣位等因素,結合圖5所示的仿真結果以確定D0的取值。

4.2 尾后搜潛陣位相對艦艇編隊變化

為查明尾后是否有潛艇跟蹤,水面艦艇編隊以30節航速航行;任務要求直升機主航向搜索速度為15節、在艦艇編隊尾后10海里處建立攔截線,然后逐步增大與編隊的距離。

尾后對潛警戒扇面將隨著直升機與編隊距離的增大而減少。為研究這種變化趨勢,仿真計算了各探測點對應的尾后最大扇面的中心角,如圖6所示。根據圖6可知,由于各探測點距離編隊的距離逐漸增加,相同兵力數量下,所能形成的尾后最大扇面的中心角逐漸減小;完成第17個點的探測后,與編隊相距100km,即使出動4架直升機,此時最大尾后扇面也只有50.6°。

根據海域水聲環境,預測敵潛艇傳感器最大作用距離不超過60km,則直升機只需搜索9個探測點;完成搜索后可飛行至編隊尾后10海里處重新開始搜索;若出動3架直升機,則至少可覆蓋尾后51°的扇面。

5 結 論

反潛直升機配置于水面艦艇編隊尾后,使用吊放聲納采用鋸齒形法搜潛可達成抗擊尾后水下威脅的戰術意圖。通過建模與仿真分析得出:

1)若在尾后15海里處實施搜潛,則2架反潛直升機可建立80.9°的對潛警戒扇面。

2)吊放聲納相鄰探測點間距系數對尾后扇面大小影響較大,鋸齒形法搜索時的最佳間距系數介于1.1~1.3之間。

3)出動3架反潛直升機,可查明尾后51°扇面、60km范圍內是否有跟蹤潛艇。

由此可見,本文建立的數學模型較好地解決了反潛直升機使用吊放聲納尾后搜潛中實際對潛警戒扇面、出動兵力數量的計算問題,可為反潛直升機執行艦艇編隊尾后搜潛任務,提供一定的理論依據和實踐參考。

[1] 羅木生,姜青山,侯學隆.反潛直升機艦艇編隊尾后截擊反潛建模[J].電光與控制,2012,19(12):6-9.LUO M S,JIANG Q S,HOU X L.Modeling of helicopter stern interception antisubmarine for surface ship formation[J].Electronics Optics & Control,2012,19(12):6-9(in Chinese).

[2] SCHERER S D.Game-theoretic anti-submarine warfare mission planner[D].Monterey,CA:Naval Postgraduate School,2009.

[3] MANVI S S,MANJULA R B.Issues in underwater acoustic sensor networks[J].International Journal on Computer and Electrical Engineering,2011,3(1):101-111.

[4] SIMAKOV S.Localisation in airborne multistatic sonars[J].IEEE Journal of Oceanic Engineering,2008,3(3):278-288.

[5] KROUT D W,FOX W J L,EL-SHARKAWI M A.Probability of target presence for multi-static sonar ping sequencing[J].IEEE Journal of Oceanic Engineering,2009,34(4):603-609.

[6] 郭傳福,王洪勝,曲延明.驅護艦編隊使用拖曳聲納反潛[J].艦船電子工程,2015,35(2):46-48.GUO C F,WANG H S,QU Y M.Research of configuaration of ship formation using towed sonar to search submarine[J].Ship Electronic Engineering,2015,35(2):46-48(in Chinese).

[7] 羅木生,姜青山,沈培志,等.反潛直升機側翼伴隨巡邏搜潛建模與仿真[J].航空學報,2015,36(3):921-928.LUO M S,JIANG Q S,SHEN P Z,et al.Modeling and simulation of antisubmarine helicopter flank-accompany patrol search[J].Acta Aeronautica et Astronautica Sinica,2015,36(3):921-928(in Chinese).

[8] CESAR J R.A reactive target active ASW sonar search tactical decision aid[D].Monterey,CA:Naval Postgraduate School,1996.

[9] OZOLS S,FEWELL M P,THREDGOLD J M.Track-initiation probability for multistatic sonar fields:DSTO-TN-1021[R].Edinburgh:Defence Science and Technology Organization,2011.

[10] FIALKOWSKI J M,GAUSS R.Methods for identifying and controlling sonar clutter[J].IEEE Journal of Oceanic Engineering,2011,35(2):330-354.

[11] 羅木生,姜青山,侯學隆.正態圓分布下的雙機吊放聲納擴展方形搜潛[J].電光與控制,2012,19(7):1-4.LUO M S,JIANG Q S,HOU X L.Modeling of search submarine distributed in a circular bivariate normal distribution by two helicopters in extended squared search using dipping sonar[J].Electronics Optics &Control,2012,19(7):1-4(in Chinese).

[12] 吳芳,楊日杰.多機吊放聲納應召搜潛建模與仿真[J].航空學報,2009,30(10):1948-1953.WU F,YANG R J.Model building and simulation based on definite second time submarine search by multi-aircraft dipping sonar[J].Acta Aeronautica et Astronautica Sinica,2009,30(10):1948-1953(in Chinese).

[13] 戚學文,嚴建鋼,謝宇,等.直升機吊放聲納應召搜潛模[J].火力與指揮控制,2014,39(3):417-420.QI X W,YAN J G,XIE Y,et al.Study on antisubmarine searching models of dipping sonar carried by helicopters[J].Fire Control & Command Control,2014,39(3):417-420(in Chinese).

[14] 白晶,韓亮.艦載直升機吊聲的搜潛方法探討[J].北京航天航空大學學報,2007,33(3):282-285.BAI J,HAN L.Search method of dipping sonar on ship borne antisubmarine aircraft[J].Journal of Beijing University of Aeronautics and Astronautics,2007,33(3):282-285(in Chinese).

[15] 吳芳,楊日杰,徐俊艷.對潛的吊放聲納應召搜索技術仿真研究[J].系統仿真學報,2009,21(13):3989-3992.WU F,YANG R J,XUE J Y.Simulation research on methods of definite second time search of dipped sonar against submarine[J].Journal of System Simulation,2009,21(13):3989-3992(in Chinese).

[16] 羅光成,張丹.反潛直升機鋸齒形檢查搜潛建模與仿真技術研究[J].電子測量技術,2011,34(10):34-36.LUO G C,ZHANG D.Modeling and simulation of the anti-submarine helicopter based on serrated check searching[J].Electronic Measurement Technology,2011,34(10):34-36(in Chinese).

[17] 羅木生,侯學隆,徐君明.基于Vega Prime的直升機反潛機動三維視景仿真[J].火力與指揮控制,2012,37(8):143-146.LUO M S,HOU X L,XU J M.3Dscene simulation of helicopter antisubmarine search based on vegaprime[J].Fire Control &Command Control,2012,37(8):143-146(in Chinese).

[18] MARTIN J.Multiplying the effectiveness of helicopter ASW sensors[J].Sea Technology,2006,47(11):33-36.

[19] 羅木生,姜青山,侯學隆.直升機使用吊聲應召反潛兵力需求仿真[J].系統仿真學報,2012,24(6):1277-1281.LUO M S,JIANG Q S,HOU X L.Simulation of the optimum helicopter force in definite second time submarine search by dipping sonar[J].Journal of System Simulation,2012,24(6):1277-1281(in Chinese).

[20] DONALD R D,HEMSTETER K P.GRASP multi-sensor search tactics against evading targets[R].Washington,D.C.:Naval Research Laboratory,2002.

Antisubmarine helicopter searching stern of ship formation adopting sawtooth method with dipping sonar

LUO Musheng*,ZENG Jiayou,HOU Xuelong
Department of Command,Naval Aeronautical and Astronautical University,Yantai 264001,China

The warning sector,demand of force and parameter decisions in helicopter antisubmarine search for threat from stern of ship formation with dipping sonar are investigated.The process of antisubmarine search with dipping sonar and the features of saw-tooth method are analyzed.A quantitative model for warning sector is built for multiple helicopters antisubmarine search at stern of ship formation.The influencing factors of saw-tooth method are analyzed,such as search velocity,and the spacing coefficient of adjacent detecting points.The relationship between submarine velocity and detection range of dipping sonar is investigated.The method for calculating equivalent circle radius of the detected area is provided.The actual warning sector against submarine,the optimum spacing coefficient,and the force demand when antisubmarine helicopters search the stern of ship formation using dipping sonar are simulated.Results show that deploying 2helicopters to 15nmiles at the stern of ship formation could cover warning sector of 80.9°,and the optimum spacing coefficient is between 1.1and 1.3.

antisubmarine search at stern of ship formation;saw-tooth method;warning sector against submarine;dipping sonar;antisubmarine helicopter;followingsubmarine

2016-02-23;Revised:2016-03-27;Accepted:2016-05-17;Published online:2016-05-30 15:02

URL:www.cnki.net/kcms/detail/11.1929.V.20160530.1502.006.html

National Social Science Foundation of China(14GJ003-154)

V275+.1;E843

A

1000-6893(2017)01-320154-07

http:/hkxb.buaa.edu.cn hkxb@buaa.edu.cn

10.7527/S1000-6893.2016.0156

(責任編輯:蘇磊)

*Corresponding author.E-mail:luosen_2008@163.com

主站蜘蛛池模板: 精品无码视频在线观看| 欧美综合区自拍亚洲综合天堂| 91精品国产情侣高潮露脸| 99久久精品视香蕉蕉| 国产成人精品男人的天堂| 亚洲第一极品精品无码| 一级香蕉人体视频| 人妻丰满熟妇αv无码| 小13箩利洗澡无码视频免费网站| 一区二区三区精品视频在线观看| 中文字幕首页系列人妻| 亚洲黄网视频| 青草视频网站在线观看| 欧美综合中文字幕久久| 国产欧美视频在线| 欧美日韩国产高清一区二区三区| 欧美精品v欧洲精品| 日本三级欧美三级| 国产肉感大码AV无码| 天天躁夜夜躁狠狠躁图片| 国产黄在线免费观看| 精品無碼一區在線觀看 | 欧美国产日韩在线| 国产v精品成人免费视频71pao| 日韩乱码免费一区二区三区| 永久在线精品免费视频观看| 欧美日韩精品一区二区在线线| 日韩av电影一区二区三区四区| 久久综合婷婷| 三级毛片在线播放| 国产91熟女高潮一区二区| 高清不卡毛片| 久久婷婷综合色一区二区| 久久伊人久久亚洲综合| 黄色免费在线网址| 正在播放久久| 91精品日韩人妻无码久久| 国产精品免费福利久久播放| www亚洲天堂| 亚洲国产无码有码| 精品人妻一区二区三区蜜桃AⅤ| 成人在线视频一区| A级毛片高清免费视频就| 国产香蕉在线视频| 亚洲激情99| 日韩精品无码免费一区二区三区 | 中文字幕有乳无码| 91麻豆精品国产高清在线| 91成人在线观看| 高清免费毛片| 亚洲国产日韩欧美在线| 狠狠色噜噜狠狠狠狠色综合久| 青草免费在线观看| 狠狠做深爱婷婷综合一区| 在线观看国产精品第一区免费 | 天堂成人在线视频| 国产精品毛片在线直播完整版| 综合色区亚洲熟妇在线| 精品国产一二三区| 色综合久久88色综合天天提莫| 91丝袜乱伦| 青青热久麻豆精品视频在线观看| 国产成人1024精品| 一区二区三区在线不卡免费| 92精品国产自产在线观看| 欧美啪啪精品| 欧美日韩中文国产va另类| 国产精品尤物在线| 久久国产精品嫖妓| 99视频全部免费| 亚洲中文无码h在线观看 | 国产无遮挡裸体免费视频| 精品国产中文一级毛片在线看 | 久久综合色视频| 丁香六月综合网| 日韩高清中文字幕| 亚洲 欧美 偷自乱 图片| 亚洲国产精品一区二区高清无码久久 | 国产成人超碰无码| jizz亚洲高清在线观看| 久久这里只精品热免费99| 亚洲成人一区二区三区|