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

    吉林松原流星進入大氣層過程及對地面的影響分析

    2021-03-25 12:30:02黨雷寧梁世昌黃潔柳森
    空間碎片研究 2021年4期
    關鍵詞:松原解體隕石

    黨雷寧,梁世昌,黃潔,柳森

    (1.中國空氣動力研究與發(fā)展中心超高速所,綿陽 621000;2.中國空氣動力研究與發(fā)展中心超高速碰撞研究中心,綿陽 621000)

    1 引言

    地球經(jīng)常遭受小行星和彗星等近地天體的撞擊。這些撞擊事件對地球環(huán)境和生命演化產(chǎn)生重要影響,也是人類生存和發(fā)展的潛在威脅之一[1]。著名的小行星撞擊事件有6500萬年前導致恐龍滅絕的墨西哥??颂K魯伯地區(qū)小行星撞擊事件[2]、1908年俄羅斯的通古斯 (Tunguska)大爆炸事件[3]、2013年俄羅斯的車里雅賓斯克 (Chelyabinsk)小行星撞擊事件[4]。這三次撞擊事件由于撞擊體直徑較大而發(fā)生頻率較低,分別約為1億年、1000年和100年[5]。人們頻繁遇到的是小直徑撞擊體進入地球大氣的流星事件。研究流星現(xiàn)象不僅可對近地天體的起源和演化等行星科學問題提供幫助[6],還能對小行星進入大氣效應等行星防御科學技術問題研究提供支撐[7]。

    國外自20世紀50年代開始研究流星現(xiàn)象[6]。經(jīng)過幾十年的發(fā)展,美國、加拿大、捷克、澳大利亞等國已經(jīng)建立了完善的流星觀測網(wǎng)[6-8],如美國的PN觀測網(wǎng),歐洲的EN觀測網(wǎng),澳大利亞的DFN觀測網(wǎng)。美國、歐洲還利用監(jiān)測核試驗的次聲網(wǎng)絡以及民用、軍用衛(wèi)星[9,10],獲得了大量流星進入事件的觀測數(shù)據(jù)。美國NASA的CNEOS網(wǎng)站統(tǒng)計了自1988年以來所有流星的觀測記錄[9]。美國、歐洲和俄羅斯等還發(fā)展了從觀測數(shù)據(jù)分析重構進入體質量、進入彈道、解體時序、空中爆炸情況、燒蝕系數(shù)、流明系數(shù)、能量沉積等重要參數(shù)的模型和方法[4,11-13]。捷克科學院的Ceplecha和Borovi?ka等,從流星觀測網(wǎng)、安防視頻和次聲網(wǎng)絡觀測數(shù)據(jù)出發(fā),獲得了Ko?ice流星[11](2010年2月28日,斯洛文尼亞)、Mor'avka流星[12](2000年5月6日,捷克)的進入大氣彈道、解體情況、能量沉積、隕石落點和質量。加拿大西安大略大學的Brown[14]和俄羅斯科學院的Popova等[4],基于觀測數(shù)據(jù)對車里雅賓斯克 (Chelyabinsk)流星開展了系統(tǒng)性的研究,獲得了流星進入大氣前的軌道、母體質量、進入大氣彈道、解體、能量沉積、隕石落區(qū)與質量、沖擊波效應等結果。

    我國正著手發(fā)展流星觀測網(wǎng),但目前覆蓋面低,對流星現(xiàn)象的研究較少且研究水平落后國外。1976-1986年,我國科學家通過現(xiàn)場考察獲得了吉林隕石雨 (1976年3月8日)進入大氣方向與空爆高度。1998年,中國-荷蘭獅子座流星雨聯(lián)合觀測隊在紫金山天文臺德令哈觀測站對獅子座流星雨進行了目視觀測[15],獲得了流星雨的分布特征。2020-2021年,針對青海玉樹流星 (2021年12月20日),紫金山天文臺基于安防視頻獲得了進入彈道,并搜索到隕石[16];中科院空間中心開展了進入彈道與落點的計算分析[17]。中國空氣動力研究與發(fā)展中心超高速所自2018年以來發(fā)展了小行星進入與撞擊效應分析評估軟件AICA(Asteroid Impact Consequence Analysis)[18,19],能夠開展小行星進入大氣過程及對地面危害的分析,并應用該軟件對車里雅賓斯克 (Chelyabinsk)、通古斯 (Tunguska)、青海玉樹等流星事件開展了計算分析。

    北京時間2019年10月11日0時16分,我國吉林松原地區(qū)發(fā)生流星進入大氣層事件。針對該事件,本文首先分析了流星進入過程的觀測數(shù)據(jù),獲得了流星母體質量、材料、飛行方向、亮度變化曲線、隕石搜索情況等,然后利用項目組建立的小行星進入與撞擊效應分析評估軟件AICA,分析了松原流星進入大氣過程及對地面的影響,給出流星母體結構的合理推測。本文的計算分析與觀測數(shù)據(jù)基本吻合,是對松原流星事件的合理解釋,為分析同類現(xiàn)象提供了借鑒。

    2 觀測數(shù)據(jù)及分析

    觀測數(shù)據(jù)包括美國NASA CNEOS網(wǎng)站的數(shù)據(jù)、現(xiàn)場視頻以及隕石搜尋結果。

    2.1 CNEOS網(wǎng)站數(shù)據(jù)

    美國利用軍用衛(wèi)星,通過可見光波段及紅外波段的觀測,在CNEOS網(wǎng)站[9]給出了松原流星最大亮度 (空爆點)位置、速度以及總輻射能,并根據(jù)Brown的經(jīng)驗公式[20]計算出進入動能。CNEOS網(wǎng)站的數(shù)據(jù)如表1所示,表中的速度分量相對于WGS84地固坐標系。表1中進入動能單位是kt,即千噸TNT當量,1kt=4.184×1012J。

    表1 吉林松原流星CNEOS網(wǎng)站數(shù)據(jù)[9]Table 1 Observation data of Songyuan meteor at CNEOS website[9]

    根據(jù)表1中的空爆點高度和速度,可以大致推斷松原流星的材料。松原流星的母體可能包括小行星和彗星兩類。由表1可見,松原流星在空中爆炸的高度為47.3km,此時流星的駐點壓力ρairV2=0.27MPa,而解體早于空爆,解體時的駐點壓力更小。根據(jù)目前天文學廣泛采用的 “駐點壓力大于強度即解體”的解體判據(jù)[1],流星的解體強度應小于0.27MPa。鐵質小行星材料的強度為40~200MPa[21],遠大于松原流星的解體強度。彗星材料的強度為0.001~0.01MPa[21],小于松原流星解體強度一個量級。長周期彗星撞擊地球的平均速度為55km/s[21],遠高于松原流星空爆點的速度14.1km/s。短周期彗星平均撞擊地球速度與松原流星空爆點速度偏差相對小,但較容易被發(fā)現(xiàn)[5],而且沒有任何公開報道和文獻給出松原流星母體為彗星的線索。因此,松原流星母體材料是鐵質小行星和彗星的可能性很小。Popova統(tǒng)計了大量流星進入大氣層數(shù)據(jù),認為石質小行星的解體強度為0.1~10MPa[22],松原流星的解體強度恰好在此范圍。根據(jù)以上分析,本文認為松原流星的母體是一顆結構松散 (解體強度小)的石質小行星。

    由表1的總撞擊能量和空爆點速度,假定進入點 (高度100km)速度與空爆點速度偏差小(4.2節(jié)證明了這一假定的合理性),可得到流星母體質量為24174.4kg。研究指出,石質小行星的密度與具體組成 (含水、不含水)、隕落概率、孔隙率和強度相關[5,23]。我們?nèi)∶芏确秶母怕首畲笾?500kg/m3為松原流星材料密度,則松原流星母體直徑為2.642m。

    根據(jù)表1中空爆點的速度分量,可計算得到流星空爆點航跡角為 -59.4°,航向角為36.6°。飛行航跡角定義[24]為速度矢量與當?shù)厮矫娴膴A角,速度矢量在當?shù)厮矫嫔戏綖檎?。飛行航向角定義[24]為速度矢量在當?shù)厮矫娴耐队芭c當?shù)卣龞|方向的夾角,偏北為正。因此松原流星的彈道陡峭,自西南往東北方向進入大氣層。

    2.2 視頻數(shù)據(jù)

    松原當?shù)匕卜涝O備與個人手機拍攝到流星進入大氣的視頻,并上傳到互聯(lián)網(wǎng)。圖1為某視頻截圖[25],該視頻是能夠獲得的從固定設備拍攝的唯一視頻,其他視頻都存在設備搖晃而不便分析。圖中各子圖的時間為相對視頻初始時刻。由圖可見,流星進入過程中亮度隨時間先增大后減小,彈道整體上比較陡,航跡角絕對值粗略看來應大于45°,與上文從CNEOS數(shù)據(jù)分析得到的航跡角在定性上符合。

    圖1 吉林松原流星視頻截圖Fig.1 Images from video of Songyuan meteor

    本文對視頻每幀圖像的灰度值求和,得到總灰度值的序列。假定從進入點到空爆點流星飛行速度保持不變且彈道為直線,得到了總灰度值隨時間的變化曲線 (圖2),其中時間起點為進入點。需要說明的是,從進入點到空爆點彈道為直線是從視頻數(shù)據(jù)分析進入彈道的合理假設,已經(jīng)應用到國外很多流星事件的分析中[11,12,14]。Borovi?ka等[26]在 Jesenice(2009年4月9日, 斯洛文尼亞)流星事件的視頻數(shù)據(jù)重構亮度曲線中,通過與月亮亮度比較,得到流星的視星等,計算公式為

    圖2 松原流星進入視頻總灰度值隨時間變化(時間起點為進入點)Fig.2 Total grey value of Songyuan meteor video(time begins at entry point)

    式中,M是火球的視星等,Mref是月亮的視星等,B是每幀視頻中的像素灰度值之和,Bm是當月亮單獨存在時像素灰度值的和,Bd是無月亮、無火球時像素灰度值的和,B、Bm、Bd都是經(jīng)過對拍攝設備校核而得到的數(shù)據(jù)。松原流星視頻無月亮或其它參考天體,且無法對視頻拍攝設備進行校核,因此無法利用式 (1)得到視星等。天文學中著名的Pogson公式[27]為

    式中,m為視星等,E為光輻射能量。小行星進入動能Et與光輻射能量Er、流明效率η存在如下關系

    對比 (2.1)和 (2.2),流星視頻圖像總灰度值與光輻射能量功率dE/dt成正比。再根據(jù)(2.3),本文認為流星視頻圖像總灰度值與動能功率或能量沉積功率dEt/dt成正比。由圖2可見,流星視頻總灰度值先緩慢變化,在3.3 s急劇增加且在4.35 s達到峰值,隨后急劇下降至5.7 s左右不再變化?;叶惹€的下降段還出現(xiàn)多個局部峰值,說明流星在主空爆后還經(jīng)歷了數(shù)次小空爆。

    2.3 隕石搜尋結果

    根據(jù)新聞報道,吉林博物館和其他隕石搜尋者曾組織隕石搜尋,未發(fā)現(xiàn)隕石。作者和吉林博物館聯(lián)系并通過其他渠道,證實了這一消息。

    3 進入過程及對地面影響計算方法

    本文采用作者團隊發(fā)展的小行星進入與撞擊效應分析評估軟件AICA,對松原流星進入大氣過程及對地面影響開展計算分析。本節(jié)簡要介紹AICA軟件中的模型,對模型的詳細論述可參見文獻[18,19]。

    3.1 彈道方程

    圓球自轉地球模型下,僅考慮飛行中的阻力,單個飛行體的彈道方程為:

    式中,h為飛行高度,θ和φ分別是經(jīng)度、緯度,V是飛行速度大小,γ是飛行航跡角或進入角,ψ是航向角,其定義已經(jīng)在2.1節(jié)中闡述;Cd是阻力系數(shù),計算中取1.0;A是橫截面積,在計算中假定小行星為球形,則A=πr2;rt是小行星到地心的距離;采用美國1976年標準大氣計算大氣密度ρair。

    3.2 質量損失方程

    小行星由于燒蝕而不斷損失質量,質量損失方程為:

    式中,m是小行星瞬時質量,σ是燒蝕系數(shù),CH是熱流系數(shù),Q是燒蝕熱。

    3.3 解體判據(jù)

    使用如下解體判據(jù):

    式中,P是小行星駐點壓力,S是解體強度。即當駐點壓力超過強度,小行星解體。

    小行星解體后,碎片的強度和質量服從Weibull定律:

    式中,α為強度指數(shù),下標 “c”和 “p”分別代表子碎片和父碎片。

    3.4 解體模型

    采用NASA的FCM解體模型。FCM模型有2016和2018兩個版本,本文分別稱為FCM(2016)和FCM(2018)。FCM(2016)模型[28]認為,小行星每次解體為一系列大尺寸的碎片和一個由小碎片、液滴組成的碎片云。所有的碎片都是獨立的飛行體,其運動服從方程組 (4)和 (5)。碎片云的運動按照Pancake模型[21]計算。FCM(2016)模型中強度指數(shù)、碎片云質量分數(shù)以及每次解體后的子碎片數(shù)都保持不變,限制了其應用。在FCM(2018)模型[29](圖3)中,小行星的初始結構由若干結構組和一個碎片云組成,每個結構組都包含數(shù)量不等的、尺寸和力學性質相同的子結構,每個子結構在達到解體條件時都按照FCM(2016)模型解體。FCM(2018)模型對小行星母體結構靈活的處理方式,使其可以應用于碎石堆復雜結構。Wheeler等[29]應用FCM(2018)模型,研究了 Ko?ice、 Bene?ov(1991年5月7日, 捷克)、Tagish Lake(2000年1月18日,加拿大)等小行星進入事件,獲得與觀測數(shù)據(jù)一致的能量沉積和地面隕石質量計算結果。

    圖3 FCM(2018)模型示意圖[29]Fig.3 Sketch of FCM(2018)

    3.5 地面損傷模型

    小行星進入大氣過程對地面的損傷包括兩類:一是沖擊波引起的地面超壓損傷,二是火球產(chǎn)生的熱輻射損傷。

    應用靜態(tài)點爆炸模型和基于美國核試驗數(shù)據(jù)的經(jīng)驗公式,評估小行星進入大氣對地面的損傷。靜態(tài)點爆炸模型把小行星的空中爆炸等效于發(fā)生在最大能量沉積高度處的點爆炸[1]。對于這個點爆炸,采用Collins的經(jīng)驗公式[1]計算沖擊波引起的地面超壓。Collins的經(jīng)驗公式通過擬合美國核試驗數(shù)據(jù),把自空爆點正下方 (文獻中稱為ground zero)沖擊波的影響范圍分為正則反射與馬赫反射兩個區(qū)域分別計算,公式的具體形式詳見文獻[1]、[19]。

    對于熱輻射損傷,Collins基于美國核試驗數(shù)據(jù)和理論分析,推導出在地面空爆時火球熱輻射損傷范圍的經(jīng)驗公式[1],Wheeler等將其修改到適用于在某高度的空爆情形。具體公式詳見文獻[1]、[19]。

    4 進入過程及對地面影響計算結果和分析

    4.1 流星進入點與結構參數(shù)

    AICA軟件的輸入?yún)?shù)包括流星進入點參數(shù)、母體結構參數(shù)與模型參數(shù)。與2.1和2.2相同,仍舊假設流星從進入點下降到空爆點速度變化小且彈道為直線,根據(jù)CNEOS的空爆點位置和速度,得到了進入點的位置、速度、航跡角和航向角,如表2所示。

    表2 松原流星進入點參數(shù)Table 2 Kinetic parameters of Songyuan meteor at entry point

    根據(jù)2.1節(jié)的結果,流星母體質量為24174.4 kg,密度取為2500 kg/m3,母體直徑為2.642 m。燒蝕系數(shù)取為石質小行星材料的典型量值[1,23],其中碎片燒蝕系數(shù)為1.0×10-8kg/J,碎片云燒蝕系數(shù)為7.0×10-9kg/J。流明效率根據(jù)俄羅斯科學院的數(shù)值仿真結果取為17%[30]。計算過程中,首先假設松原流星母體結構為完整巖石,計算得到能量沉積曲線 (圖4),能夠復現(xiàn)總灰度曲線 (圖2)大致變化趨勢,但無論怎樣調(diào)整輸入?yún)?shù)都無法獲得主空爆點后的幾次空爆點特征。

    圖4 假設松原流星為完整巖石的能量沉積計算結果與視頻總灰度曲線對比Fig.4 Comparison between computational results of energy deposition and total grey value of video assuming that the parent body is a monolith

    本文假設松原流星為碎石堆結構,包含風化層以及4個由完整巖石組成的結構組,風化層質量分數(shù)為70%,結構組質量分數(shù)為30%,母體初次解體強度為0.018 MPa,各結構組的子結構數(shù)量以及子結構解體強度、強度指數(shù)、碎片云質量分數(shù)、子碎片數(shù)與質量分布詳見表3。子結構碎片云質量分數(shù)對能量沉積與墜落地面隕石質量有著重要影響,為了得到與觀測數(shù)據(jù)一致的計算結果,本文對前3塊完整巖石取80%,第4塊完整巖石取100%。通古斯 (Tunguska)大爆炸地面無隕石以及Chyba[21]采用Pancake模型對此事件的計算,實際上說明了碎片云質量分數(shù)取較大值甚至100%的合理性。Wheeler等[29]采用FCM(2018)模型對 Chelyabinsk、 Bene?ov、 Tagish Lake 等流星事件能量沉積與地面隕石質量的計算中,碎片云質量分數(shù)都取到了50%以上甚至100%的量值。

    表3 松原流星母體結構參數(shù) (初次解體強度為0.018 MPa)Table 3 Parent body's structure of Songyuan meteor(strength of 1stbreakup is 0.018 MPa)

    在表3所示的流星母體結構參數(shù)下,本文完全復現(xiàn)了總灰度曲線的變化趨勢 (圖5),包括主空爆點之后的二次、三次、四次空爆點特征。因此,本文認為表3所示的流星母體結構參數(shù)是松原流星結構的合理推斷。對圖4~5說明如下:紅色曲線是計算結果,縱坐標取右側縱軸;黑色曲線是總灰度曲線,縱坐標取左側縱軸;因無法把總灰度轉換為能量沉積 (見2.2節(jié)),僅對比計算和觀測曲線的形狀,而不是量值。

    圖5 假設松原流星為表3所示碎石堆結構的能量沉積計算結果與視頻總灰度曲線對比Fig.5 Comparison between computational results of energy deposition and total grey value of video assuming that the parent body is rubble pile with structure of Table 3

    4.2 進入過程分析

    圖6給出在表3所示碎石堆結構參數(shù)下的能量沉積組成,圖中 “regolith”表示表3中的風化層,“rubble 1”~“rubble 4” 分別表示表3中的4個完整巖石結構組。圖7是松原流星的進入彈道,圖中“asteroid”表示未解體的流星母體,“rubble 1”~“rubble 4”與圖6的意義相同,各實心符號表示結構組中子結構的初次解體點,圖7(a)是速度-高速曲線,圖7(b)是縱向航程-高度曲線。綜合圖6和圖7,流星母體進入大氣后首先在高度69.8km解體,釋放出的風化層和4個結構組;占母體質量70%風化層的減速和燒蝕使得流星能量沉積在主空爆點達到最大;結構組 “rubble 1”~“rubble 4”基于自身的強度 (見表3)分別在高度47.9km、44.9km、39.1km、40.6km繼續(xù)解體,并分別釋放出占其子結構質量分數(shù)80%、80%、80%、100%的碎片云,從而對主空爆點后的能量沉積做出貢獻,使得主空爆點后出現(xiàn)二次、三次、四次空爆。進入點 (高度100km)與空爆點 (47.3km)的速度偏差小 (圖7a),進入彈道近似為直線 (圖7b),說明2.1、2.2、4.1節(jié) “從進入點到空爆點速度變化小以及彈道為直線”的假設正確。

    圖6 松原流星進入大氣能量沉積組成Fig.6 Constitution of energy deposition

    圖7 松原流星進入彈道計算結果Fig.7 Computational results of entry trajectory

    圖8給出松原流星總質量以及大尺寸碎片總質量隨高度的變化情況??梢钥吹?當飛行至高度30 km時,所有大尺寸碎片都通過解體或燒蝕的方式消耗完質量,大氣中留下由小碎片、液滴、灰塵組成的碎片云。在表3的計算條件下,對質量大于0.01 kg(直徑約2 cm)的隕石進行統(tǒng)計,地面隕石總質量為1.04 kg,其中最大隕石質量0.065 kg(直徑約3.4 cm)。隕石總質量較小且考慮到高空風對隕石散布影響較大,搜索到隕石的可能性較小。這是地面搜索不到隕石的重要原因。

    4.3 對地面影響分析

    圖9給出松原流星進入大氣沖擊波引起的地面超壓曲線,圖中橫軸是距空爆點星下點的距離??梢钥吹?地面超壓隨距離快速遞減,空爆點星下點的地面超壓最大 (量值為76 Pa),遠小于導致玻璃破碎的最小超壓500 Pa[4]。因此,松原流星的沖擊波不會對地面造成損傷。

    圖9 松原流星空中爆炸沖擊波地面超壓計算結果Fig.9 Computational results of overpressure on ground caused by airburst

    圖10給出進入過程熱輻射1級燒傷損傷半徑與空爆高度的關系。由圖可見,在松原流星的進入動能下,只有當空爆高度低于1.3 km時才能對地面造成1級燒傷。松原流星的實際空爆高度為47.3 km,因此不會對地面造成熱輻射損傷。

    圖10 松原流星進入大氣引起地面1級燒傷的損傷半徑與空爆高度的關系Fig.10 Radius of 1st degree burn versus height of airburst

    5 結論

    本文針對2019年10月11日發(fā)生在我國吉林松原地區(qū)的流星進入事件,結合NASA CNEOS網(wǎng)站的數(shù)據(jù)、現(xiàn)場視頻以及隕石搜尋結果等觀測數(shù)據(jù),利用項目組建立的小行星進入與撞擊效應分析評估軟件AICA,計算分析了松原流星進入大氣過程及對地面的影響,給出流星母體結構的合理推測。本文的計算分析指出:

    (1)松原流星的母體質量約為2.4 t,是具有碎石堆結構的石質小行星。流星母體結構的合理推測為:風化層質量占70%,4個結構組共5塊碎石占30%。

    (2)松原流星進入大氣過程中,在69.8 km高度初次解體,在47.3 km高度達到最大能量沉積即發(fā)生主空爆,隨后又發(fā)生了若干次小空爆,在30 km高度以上所有大尺寸碎片都通過解體或燒蝕的方式消耗完質量,沒有大尺寸隕石降落到地面。

    (3)流星進入大氣沖擊波引起的地面超壓最大為76Pa,且實際空爆高度大于引起1級燒傷的最大空爆高度,不會對地面人員和實施造成損傷。

    猜你喜歡
    松原解體隕石
    南極為何多隕石
    隕石獵人
    河南松原生物科技有限公司
    中國釀造(2019年9期)2019-10-08 05:44:02
    愛的花旗松素 走進吉林松原市場
    隕石
    幼兒畫刊(2018年7期)2018-07-24 08:26:14
    幸運隕石
    蘇聯(lián)1991年解體前的最后時光
    中外文摘(2017年14期)2017-07-31 16:16:48
    啊,松原
    黃河之聲(2016年13期)2016-11-09 01:42:16
    “娃娃親”因兩家發(fā)展不同而解體
    美空軍又一退役氣象衛(wèi)星在軌解體
    太空探索(2016年12期)2016-07-18 11:13:43
    久久精品夜夜夜夜夜久久蜜豆| 床上黄色一级片| 91aial.com中文字幕在线观看| 久久99热这里只有精品18| 麻豆国产97在线/欧美| av在线观看视频网站免费| 麻豆成人午夜福利视频| 人妻久久中文字幕网| 丰满的人妻完整版| 久久久午夜欧美精品| 亚洲国产日韩欧美精品在线观看| 日韩一区二区三区影片| 男人的好看免费观看在线视频| 国产精品女同一区二区软件| 亚洲国产欧美在线一区| 久久国产乱子免费精品| 在线观看av片永久免费下载| 夜夜看夜夜爽夜夜摸| 国产伦精品一区二区三区四那| 日韩视频在线欧美| 国产精品一区二区性色av| 成人一区二区视频在线观看| 国产精品久久久久久精品电影| 日韩欧美精品免费久久| 日韩欧美精品v在线| 亚洲经典国产精华液单| 午夜免费男女啪啪视频观看| 99久国产av精品| 在线免费十八禁| 久久韩国三级中文字幕| 97超视频在线观看视频| 少妇丰满av| www.av在线官网国产| 麻豆一二三区av精品| 老司机福利观看| 日本一二三区视频观看| 嫩草影院精品99| 国产精品人妻久久久影院| 国产精品久久久久久精品电影| 日韩中字成人| 亚洲无线观看免费| 男女那种视频在线观看| 亚洲第一电影网av| 春色校园在线视频观看| 婷婷色av中文字幕| 成人高潮视频无遮挡免费网站| 国产极品天堂在线| av在线天堂中文字幕| 国产亚洲精品久久久com| 九九爱精品视频在线观看| 国产午夜精品论理片| 久久热精品热| 久久久精品欧美日韩精品| 麻豆国产av国片精品| 国产 一区 欧美 日韩| 国产精品.久久久| 最好的美女福利视频网| 欧美日韩乱码在线| 婷婷色av中文字幕| .国产精品久久| 久久精品夜夜夜夜夜久久蜜豆| 成人午夜高清在线视频| 亚洲国产精品国产精品| 插逼视频在线观看| 精品午夜福利在线看| av卡一久久| 亚洲久久久久久中文字幕| 国产毛片a区久久久久| 免费观看a级毛片全部| 超碰av人人做人人爽久久| 麻豆精品久久久久久蜜桃| 免费av观看视频| 亚洲精品自拍成人| 在线观看美女被高潮喷水网站| 久久99热6这里只有精品| a级一级毛片免费在线观看| 欧美一区二区精品小视频在线| 亚洲av成人av| 久久久久性生活片| 国内精品一区二区在线观看| avwww免费| 丝袜美腿在线中文| 欧美成人免费av一区二区三区| 一个人看视频在线观看www免费| 日韩欧美精品v在线| 日本免费a在线| 美女高潮的动态| 中国美白少妇内射xxxbb| 亚洲欧美精品专区久久| 性插视频无遮挡在线免费观看| 日本三级黄在线观看| 少妇高潮的动态图| 欧美精品一区二区大全| 成人特级黄色片久久久久久久| 国产亚洲5aaaaa淫片| 午夜激情福利司机影院| 女人被狂操c到高潮| 国产日韩欧美在线精品| 国内精品一区二区在线观看| 22中文网久久字幕| 高清日韩中文字幕在线| 亚洲七黄色美女视频| 久久99热6这里只有精品| 免费人成在线观看视频色| 亚洲成av人片在线播放无| 日韩成人av中文字幕在线观看| 日韩制服骚丝袜av| 精品国产三级普通话版| av国产免费在线观看| 亚洲精品日韩av片在线观看| 简卡轻食公司| 最近视频中文字幕2019在线8| 日本免费a在线| 人妻夜夜爽99麻豆av| 18禁在线无遮挡免费观看视频| 天堂av国产一区二区熟女人妻| 九九爱精品视频在线观看| 可以在线观看毛片的网站| 变态另类成人亚洲欧美熟女| 亚洲精品日韩在线中文字幕 | 亚洲精品色激情综合| 只有这里有精品99| 桃色一区二区三区在线观看| 欧美丝袜亚洲另类| 日本欧美国产在线视频| 桃色一区二区三区在线观看| 直男gayav资源| 日韩成人伦理影院| 不卡一级毛片| 久久这里只有精品中国| 国产精品精品国产色婷婷| 午夜精品在线福利| 高清日韩中文字幕在线| 观看美女的网站| 国产麻豆成人av免费视频| 啦啦啦观看免费观看视频高清| 日韩国内少妇激情av| 亚洲美女视频黄频| 久久午夜亚洲精品久久| 欧美激情国产日韩精品一区| 久久久久久久久久黄片| 免费看av在线观看网站| 国产精品久久久久久久电影| 国产精品女同一区二区软件| 亚洲人成网站高清观看| 国产精品无大码| 好男人视频免费观看在线| 成人美女网站在线观看视频| 97人妻精品一区二区三区麻豆| 嫩草影院入口| 美女被艹到高潮喷水动态| 自拍偷自拍亚洲精品老妇| 免费搜索国产男女视频| 99热这里只有是精品50| 欧美日韩在线观看h| 日韩制服骚丝袜av| 99热这里只有是精品在线观看| 99久久九九国产精品国产免费| 免费不卡的大黄色大毛片视频在线观看 | 非洲黑人性xxxx精品又粗又长| 国产免费男女视频| 久久精品国产清高在天天线| 日韩 亚洲 欧美在线| 成人永久免费在线观看视频| 丰满人妻一区二区三区视频av| 中国美女看黄片| 亚洲第一电影网av| 看黄色毛片网站| 男人舔奶头视频| 久久精品人妻少妇| 校园人妻丝袜中文字幕| 午夜精品一区二区三区免费看| 免费看a级黄色片| 舔av片在线| 一卡2卡三卡四卡精品乱码亚洲| 亚洲国产精品合色在线| 免费不卡的大黄色大毛片视频在线观看 | 免费看光身美女| 亚洲真实伦在线观看| 一个人观看的视频www高清免费观看| 极品教师在线视频| 国产伦在线观看视频一区| 久久久成人免费电影| 日韩av在线大香蕉| 国产精品野战在线观看| 国产精品福利在线免费观看| 99久久精品国产国产毛片| 中文字幕精品亚洲无线码一区| 国产精品一及| 日韩一区二区视频免费看| 99热这里只有是精品50| 久久久久久久久久成人| 国产熟女欧美一区二区| 男人舔女人下体高潮全视频| 国产国拍精品亚洲av在线观看| 国产精品女同一区二区软件| 中国美女看黄片| 午夜福利在线在线| h日本视频在线播放| 国产爱豆传媒在线观看| 色综合亚洲欧美另类图片| 国产精品av视频在线免费观看| 美女大奶头视频| 国产熟女欧美一区二区| 亚洲久久久久久中文字幕| 亚洲欧美日韩卡通动漫| 欧美区成人在线视频| 亚洲一级一片aⅴ在线观看| 1000部很黄的大片| 18禁裸乳无遮挡免费网站照片| 国产黄色视频一区二区在线观看 | 国产不卡一卡二| 国产女主播在线喷水免费视频网站 | 久久久a久久爽久久v久久| 美女内射精品一级片tv| 久久久成人免费电影| 直男gayav资源| 日日撸夜夜添| 性插视频无遮挡在线免费观看| 99国产精品一区二区蜜桃av| 色哟哟哟哟哟哟| 网址你懂的国产日韩在线| 久久人人精品亚洲av| 国产精品99久久久久久久久| 欧美激情久久久久久爽电影| 晚上一个人看的免费电影| 狂野欧美白嫩少妇大欣赏| 亚洲精品国产成人久久av| 精品不卡国产一区二区三区| 欧美精品国产亚洲| 久久精品综合一区二区三区| 青青草视频在线视频观看| 插阴视频在线观看视频| 非洲黑人性xxxx精品又粗又长| 国产久久久一区二区三区| 日韩大尺度精品在线看网址| 老师上课跳d突然被开到最大视频| 小蜜桃在线观看免费完整版高清| 伦理电影大哥的女人| 夜夜看夜夜爽夜夜摸| 国产精品综合久久久久久久免费| 日韩欧美精品v在线| 婷婷六月久久综合丁香| 99久久无色码亚洲精品果冻| 久久久精品大字幕| 黄片无遮挡物在线观看| 搡老妇女老女人老熟妇| 亚洲第一电影网av| 不卡一级毛片| 狂野欧美激情性xxxx在线观看| 亚洲第一电影网av| 国产午夜福利久久久久久| 亚洲内射少妇av| 国产中年淑女户外野战色| 国产午夜精品一二区理论片| 亚洲欧洲国产日韩| 色播亚洲综合网| 免费一级毛片在线播放高清视频| 亚洲精品日韩在线中文字幕 | 日本撒尿小便嘘嘘汇集6| 人人妻人人看人人澡| 又粗又爽又猛毛片免费看| 噜噜噜噜噜久久久久久91| 亚洲最大成人手机在线| 国产精品国产高清国产av| 亚洲欧美成人综合另类久久久 | 欧美区成人在线视频| 欧美+亚洲+日韩+国产| 久久国产乱子免费精品| 国产一区二区激情短视频| av卡一久久| 桃色一区二区三区在线观看| 啦啦啦观看免费观看视频高清| 99久久九九国产精品国产免费| 国产精品久久久久久亚洲av鲁大| 久久久久久久久大av| 一个人免费在线观看电影| 午夜爱爱视频在线播放| 又爽又黄无遮挡网站| av黄色大香蕉| 深爱激情五月婷婷| 小说图片视频综合网站| 国模一区二区三区四区视频| 亚洲欧美日韩无卡精品| 高清日韩中文字幕在线| 三级经典国产精品| 91在线精品国自产拍蜜月| 国产av在哪里看| 欧美日韩一区二区视频在线观看视频在线 | 久久99蜜桃精品久久| 97超碰精品成人国产| 人人妻人人看人人澡| 婷婷色av中文字幕| 久久人人爽人人片av| 男人的好看免费观看在线视频| 色吧在线观看| 国产成人91sexporn| 看非洲黑人一级黄片| 2022亚洲国产成人精品| 欧美一区二区精品小视频在线| 伦理电影大哥的女人| 亚洲第一电影网av| a级毛色黄片| 国产精品久久久久久久久免| 国内精品久久久久精免费| 国产免费男女视频| 精品熟女少妇av免费看| 亚洲成人精品中文字幕电影| 深夜精品福利| 国产精品一区二区三区四区久久| 久久人人爽人人爽人人片va| 国产高清有码在线观看视频| 国产成人aa在线观看| 国产精品久久久久久久久免| 日韩视频在线欧美| 亚洲在久久综合| 日韩av不卡免费在线播放| 搡女人真爽免费视频火全软件| 国产亚洲91精品色在线| 亚洲18禁久久av| 熟女人妻精品中文字幕| 在线观看午夜福利视频| 有码 亚洲区| 成人午夜高清在线视频| 成人无遮挡网站| 99久久无色码亚洲精品果冻| 18禁黄网站禁片免费观看直播| 久久婷婷人人爽人人干人人爱| 麻豆国产97在线/欧美| 舔av片在线| 三级男女做爰猛烈吃奶摸视频| 亚洲人成网站在线播| 直男gayav资源| 人人妻人人看人人澡| 亚洲成人av在线免费| 亚洲成人精品中文字幕电影| 最新中文字幕久久久久| 国国产精品蜜臀av免费| 又粗又硬又长又爽又黄的视频 | 欧美变态另类bdsm刘玥| 国产精品久久久久久久电影| 国产一级毛片七仙女欲春2| 3wmmmm亚洲av在线观看| 综合色av麻豆| 欧美日韩乱码在线| 麻豆国产av国片精品| 国产视频内射| 卡戴珊不雅视频在线播放| 成人无遮挡网站| 少妇熟女欧美另类| 婷婷色av中文字幕| 三级经典国产精品| 黄片wwwwww| 国产黄片视频在线免费观看| 高清毛片免费观看视频网站| 在线免费十八禁| 国产午夜精品一二区理论片| 国产 一区精品| 校园人妻丝袜中文字幕| 中文字幕制服av| 女人被狂操c到高潮| 99视频精品全部免费 在线| av又黄又爽大尺度在线免费看 | 国产午夜精品久久久久久一区二区三区| or卡值多少钱| 床上黄色一级片| 国产精品美女特级片免费视频播放器| 欧美一级a爱片免费观看看| 国产精品麻豆人妻色哟哟久久 | 免费观看在线日韩| 五月伊人婷婷丁香| 国产精品国产三级国产av玫瑰| 99riav亚洲国产免费| 国产久久久一区二区三区| 国产精品99久久久久久久久| 男女啪啪激烈高潮av片| 在线观看一区二区三区| 偷拍熟女少妇极品色| 秋霞在线观看毛片| 在线观看免费视频日本深夜| 91aial.com中文字幕在线观看| 天美传媒精品一区二区| 精品久久久久久久久久免费视频| 日本-黄色视频高清免费观看| 亚洲一级一片aⅴ在线观看| 性插视频无遮挡在线免费观看| 成人性生交大片免费视频hd| 精品久久久久久成人av| 日韩在线高清观看一区二区三区| 波多野结衣高清无吗| 少妇被粗大猛烈的视频| 亚洲精品色激情综合| 久久久久久久久久黄片| 国产高清有码在线观看视频| 国产亚洲av嫩草精品影院| 超碰av人人做人人爽久久| 麻豆国产av国片精品| 亚洲人成网站高清观看| 黄色日韩在线| 午夜福利在线在线| 久久久久久久久久黄片| 亚洲va在线va天堂va国产| 国产乱人偷精品视频| 大香蕉久久网| 国产爱豆传媒在线观看| 精品人妻视频免费看| 国产精品乱码一区二三区的特点| 亚洲人成网站在线观看播放| 国产毛片a区久久久久| 99热6这里只有精品| 性色avwww在线观看| 国产乱人偷精品视频| 国产精品电影一区二区三区| 精品一区二区三区人妻视频| 国产男人的电影天堂91| 国产人妻一区二区三区在| 熟妇人妻久久中文字幕3abv| 亚洲国产日韩欧美精品在线观看| 国内少妇人妻偷人精品xxx网站| 久久精品国产亚洲av香蕉五月| 亚洲av中文av极速乱| 日韩欧美一区二区三区在线观看| 深爱激情五月婷婷| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 美女cb高潮喷水在线观看| 国产高清三级在线| 亚洲电影在线观看av| 日韩中字成人| 色播亚洲综合网| 久久精品国产亚洲网站| 国产白丝娇喘喷水9色精品| 九色成人免费人妻av| 变态另类丝袜制服| 中文字幕人妻熟人妻熟丝袜美| 麻豆久久精品国产亚洲av| 亚洲最大成人中文| 久久久成人免费电影| 亚洲av免费高清在线观看| 国产乱人偷精品视频| 久久精品国产亚洲网站| 可以在线观看的亚洲视频| 一级av片app| 亚洲在线自拍视频| 国产免费一级a男人的天堂| 久久久久久国产a免费观看| 97在线视频观看| 成人av在线播放网站| 色噜噜av男人的天堂激情| 成人三级黄色视频| 99riav亚洲国产免费| 日本成人三级电影网站| 免费搜索国产男女视频| 级片在线观看| 日韩三级伦理在线观看| 99久国产av精品| 少妇被粗大猛烈的视频| 久久人妻av系列| 久久精品影院6| 大又大粗又爽又黄少妇毛片口| 亚洲欧洲日产国产| 欧美日韩国产亚洲二区| 日本爱情动作片www.在线观看| 国产午夜精品一二区理论片| 九九热线精品视视频播放| 一卡2卡三卡四卡精品乱码亚洲| 22中文网久久字幕| 人妻少妇偷人精品九色| 亚洲色图av天堂| 日韩欧美一区二区三区在线观看| 一边亲一边摸免费视频| 色哟哟·www| 极品教师在线视频| 在线a可以看的网站| 免费av毛片视频| 亚洲成人久久爱视频| 成年版毛片免费区| 国产精品三级大全| 别揉我奶头 嗯啊视频| 国产乱人视频| 六月丁香七月| 波多野结衣巨乳人妻| 国产 一区精品| 麻豆国产av国片精品| 少妇人妻精品综合一区二区 | 老司机福利观看| 美女脱内裤让男人舔精品视频 | 大型黄色视频在线免费观看| 美女xxoo啪啪120秒动态图| 亚洲国产精品国产精品| 亚洲av中文av极速乱| 一区二区三区免费毛片| 亚洲人成网站在线观看播放| 久久婷婷人人爽人人干人人爱| 男人的好看免费观看在线视频| 在线播放无遮挡| 亚洲精品成人久久久久久| 卡戴珊不雅视频在线播放| 可以在线观看毛片的网站| 国产午夜精品久久久久久一区二区三区| 日韩 亚洲 欧美在线| 白带黄色成豆腐渣| 亚洲av成人av| 成人亚洲精品av一区二区| 12—13女人毛片做爰片一| 18禁在线播放成人免费| 欧美一级a爱片免费观看看| 精品一区二区三区人妻视频| 亚洲国产欧洲综合997久久,| 成人鲁丝片一二三区免费| 亚洲精品日韩在线中文字幕 | 久久久精品94久久精品| 日韩一区二区三区影片| 国产白丝娇喘喷水9色精品| 国产成人精品婷婷| 国产中年淑女户外野战色| 成年版毛片免费区| 五月玫瑰六月丁香| a级毛色黄片| 久久精品久久久久久久性| 成年版毛片免费区| 舔av片在线| 亚洲久久久久久中文字幕| 午夜福利在线在线| 18禁裸乳无遮挡免费网站照片| 国产伦理片在线播放av一区 | 国产伦精品一区二区三区视频9| 久久久久网色| 天天一区二区日本电影三级| 国产精品精品国产色婷婷| 免费黄网站久久成人精品| 国产一级毛片七仙女欲春2| 少妇裸体淫交视频免费看高清| 国产 一区 欧美 日韩| 欧美日韩精品成人综合77777| 国产美女午夜福利| 亚洲中文字幕一区二区三区有码在线看| 男人和女人高潮做爰伦理| 深夜a级毛片| 国产中年淑女户外野战色| 精品久久久久久久人妻蜜臀av| 国产av不卡久久| 亚洲欧美中文字幕日韩二区| 亚洲精品自拍成人| 99热全是精品| 日韩中字成人| 欧美3d第一页| 99久久久亚洲精品蜜臀av| 99在线视频只有这里精品首页| 五月伊人婷婷丁香| 老司机福利观看| 日本五十路高清| 一本精品99久久精品77| 我的女老师完整版在线观看| 国产精品蜜桃在线观看 | 99九九线精品视频在线观看视频| 黄色日韩在线| 日韩一区二区三区影片| 国产精品免费一区二区三区在线| 国产伦理片在线播放av一区 | 成人高潮视频无遮挡免费网站| 老司机影院成人| 亚洲欧美精品自产自拍| 国产精品久久久久久久久免| 99在线人妻在线中文字幕| 极品教师在线视频| 日本免费一区二区三区高清不卡| 国产淫片久久久久久久久| 一区二区三区高清视频在线| 麻豆精品久久久久久蜜桃| 美女黄网站色视频| 国产精品乱码一区二三区的特点| 亚洲内射少妇av| 久久精品国产亚洲av天美| 国产成年人精品一区二区| 免费看日本二区| 亚洲在久久综合| 亚洲av熟女| 中文字幕久久专区| 人人妻人人澡欧美一区二区| 禁无遮挡网站| 老熟妇乱子伦视频在线观看| 国产精品女同一区二区软件| 亚洲丝袜综合中文字幕| 99国产极品粉嫩在线观看| 春色校园在线视频观看| 熟妇人妻久久中文字幕3abv| 国产精品无大码| 国产淫片久久久久久久久| 亚洲中文字幕日韩| 日韩成人av中文字幕在线观看| 国产探花在线观看一区二区| 国产私拍福利视频在线观看| 校园春色视频在线观看| 亚洲乱码一区二区免费版| 久久久国产成人免费| 一级黄色大片毛片| 91在线精品国自产拍蜜月| 欧美最黄视频在线播放免费| 久久精品国产清高在天天线| 中文字幕免费在线视频6| 国产真实伦视频高清在线观看| 久久人人爽人人片av| 寂寞人妻少妇视频99o| 国产精华一区二区三区| 国产日本99.免费观看| 你懂的网址亚洲精品在线观看 | 国产精品美女特级片免费视频播放器| 免费看a级黄色片| 蜜桃久久精品国产亚洲av| 美女黄网站色视频| 久久这里有精品视频免费| 国产精品精品国产色婷婷| 波野结衣二区三区在线| 麻豆成人午夜福利视频| 国产 一区精品|