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

    當(dāng)量和爆高(埋深)反演的聲震分析方法

    2023-03-11 10:26:54張亮永張德志肖衛(wèi)國梁旭斌王同東郭權(quán)勢(shì)
    現(xiàn)代應(yīng)用物理 2023年4期
    關(guān)鍵詞:爆源當(dāng)量聲學(xué)

    張亮永,張德志,肖衛(wèi)國,梁旭斌,王同東,郭權(quán)勢(shì),李 翱

    (強(qiáng)脈沖輻射環(huán)境模擬與效應(yīng)全國重點(diǎn)實(shí)驗(yàn)室;西北核技術(shù)研究所:西安710024)

    在國防和民用領(lǐng)域中,近地面爆源參數(shù)反演是非常重要的研究內(nèi)容,快速預(yù)測(cè)爆炸當(dāng)量和爆高(埋深)(height of burnst, HOB)等源參數(shù),對(duì)爆炸事故和恐怖襲擊預(yù)警[1-4]、地震和火山噴發(fā)事件預(yù)報(bào)[5-7]、核查監(jiān)測(cè)[8-10]及武器性能評(píng)估[11-12]等方面具有重要意義。采用聲學(xué)數(shù)據(jù)(近區(qū)為沖擊波)[13-21]或地震波數(shù)據(jù)(近區(qū)為地運(yùn)動(dòng))[2-3, 22-24]預(yù)測(cè)爆源參數(shù)是常見的兩種近地面爆源參數(shù)反演方法,但采用這兩種方法反演的爆炸當(dāng)量和爆高(埋深)分布在特定范圍內(nèi)存在多個(gè)極值點(diǎn)(稱為折中關(guān)系[22-27]),給源參數(shù)估計(jì)造成了極大的不確定性。為提高源參數(shù)預(yù)測(cè)的精度,研究人員提出了聲震分析方法[22-27],該方法基于數(shù)據(jù)融合算法綜合地震波(地運(yùn)動(dòng))數(shù)據(jù)和聲學(xué)(空氣超壓)數(shù)據(jù)反演爆源參數(shù),通過同時(shí)利用聲震波動(dòng)信息對(duì)源參數(shù)進(jìn)行多重約束,減小折中關(guān)系的影響。聲震分析方法分為聲學(xué)模型、地震波模型和數(shù)據(jù)融合算法3個(gè)部分,本文將對(duì)這3部分目前的進(jìn)展展開論述。

    1 聲學(xué)模型

    由于近地面爆炸能量會(huì)耦合到空氣中產(chǎn)生明顯的聲擾動(dòng),一般通過遠(yuǎn)場(chǎng)聲學(xué)數(shù)據(jù),利用超壓峰值、正向聲沖量及正向脈寬等波形特征量來預(yù)測(cè)地面爆炸當(dāng)量[22-27],研究人員相繼提出適合快速計(jì)算,以超壓、聲沖量和持續(xù)時(shí)間為特征量的半經(jīng)驗(yàn)聲學(xué)模型及描述波形細(xì)節(jié)的全波形聲學(xué)反演方法。半經(jīng)驗(yàn)聲學(xué)模型包括ANSI模型[28, 29]、KG85標(biāo)準(zhǔn)模型[30]、BOOM模型[15, 31]和IPM聲學(xué)模型[19]等自由場(chǎng)聲學(xué)模型及地面爆炸聲學(xué)模型[2]和近地面爆炸聲學(xué)模型[22-26]。半經(jīng)驗(yàn)自由場(chǎng)聲學(xué)模型未考慮爆炸能量耦合到地下的影響,對(duì)地面爆炸事件通常假設(shè)兩倍等效爆炸當(dāng)量進(jìn)行爆源參數(shù)估計(jì)[20, 28-30, 32];而對(duì)于近地面爆炸事件,Zhang等[33-34]通過考慮地面反射、耦合作用及大氣環(huán)境因素,建立了基于自由場(chǎng)聲學(xué)模型的近地面聲學(xué)模型(簡稱為自由場(chǎng)近地面聲學(xué)模型),進(jìn)行了爆炸當(dāng)量和爆高估計(jì)。表1為不同的自由場(chǎng)聲學(xué)模型。地面爆炸聲學(xué)模型[2]是基于地面爆炸聲學(xué)數(shù)據(jù)建立的模型,考慮了地面爆炸能量耦合情況,用于地面爆炸事件的當(dāng)量估計(jì),但不適用于偏離地面位置的爆炸事件當(dāng)量估計(jì)。而近地面爆炸聲學(xué)模型考慮了近地面聲震能量耦合特性,可用于近地面爆炸事件的當(dāng)量和爆高(埋深)估計(jì),但不適用比爆高為-2.4 ~-0.6 m·kg-1/3的淺埋情況。上述兩種模型需根據(jù)不同場(chǎng)地介質(zhì)提前確定多個(gè)待定系數(shù),適用的場(chǎng)地介質(zhì)有限。Zhang等[25]在近地面爆炸聲來源機(jī)制[35-36]的基礎(chǔ)上,通過考慮氣體噴出地表和地沖擊耦合影響對(duì)近地面聲學(xué)模型進(jìn)行改進(jìn),拓寬了比例爆高適用范圍,同時(shí)確定了花崗巖夾雜變質(zhì)巖場(chǎng)地的模型待定系數(shù)。圖1為近地面爆炸聲來源機(jī)制[35-36]和改進(jìn)前后近地面聲學(xué)模型[25]。

    (a)Mechanism of near-surface explosion producing sound [35-36]

    (b)Unimproved and improved near-surface acoustic models[25]

    表1 自由場(chǎng)聲學(xué)模型

    全波形聲學(xué)反演方法包括超壓波形模板匹配法[16]、經(jīng)驗(yàn)聲源模型法[20]和折合聲沖量法[14, 17, 24]等。超壓波形模板匹配法[16]以不同比距離下隨比時(shí)間變化的大量壓力波形數(shù)據(jù)為模板,將未知當(dāng)量的全壓力波形和模板壓力波形進(jìn)行匹配,對(duì)應(yīng)匹配偏差最小的當(dāng)量為估計(jì)當(dāng)量,但該方法需大量不同當(dāng)量不同比距離的壓力波形組成的樣本庫,但目前樣本數(shù)據(jù)庫還不夠充分;經(jīng)驗(yàn)聲源模型法[20]以G17HE為標(biāo)準(zhǔn)源模型,通過大量實(shí)測(cè)視壓力數(shù)據(jù)并假設(shè)標(biāo)準(zhǔn)高斯分布,獲取源模型的峰值壓力和持續(xù)時(shí)間平均值和方差(偏差分布),由此計(jì)算爆炸當(dāng)量,但該模型聲測(cè)點(diǎn)數(shù)據(jù)需選取20~100 m·kg-1/3的比距離范圍,在未知當(dāng)量時(shí)不好提前預(yù)設(shè)測(cè)點(diǎn)距離。圖2為超壓波形樣本庫[16]以及G17HE源模型和實(shí)測(cè)波形[20]。Kim等[14, 17]基于線性聲學(xué)理論反演爆源等效源時(shí)間函數(shù),以KG85模型折合聲沖量作為參考函數(shù),提出一種地面爆炸當(dāng)量反演思路(本文稱為折合聲沖量法),并將其應(yīng)用于剛性地面爆炸當(dāng)量反演,具有較高的當(dāng)量預(yù)測(cè)精度,但該方法只適用剛性地面爆炸當(dāng)量估計(jì),沒有考慮波傳播過程中地表結(jié)構(gòu)的影響。圖3為等效源時(shí)間函數(shù)及預(yù)測(cè)和實(shí)測(cè)波形。張亮永等[37]對(duì)該方法進(jìn)行了發(fā)展,通過考慮波傳播過程中地表結(jié)構(gòu)的吸收衰減影響及近地面爆炸聲震能量分配過程,建立了多孔彈性場(chǎng)地近地面折合聲沖量模型,解決了多孔彈性場(chǎng)地的當(dāng)量估計(jì)問題,給出當(dāng)量相對(duì)偏差為7%,爆高相對(duì)偏差不超過19%。圖4為多孔彈性場(chǎng)地聲傳播示意圖和爆源參數(shù)反演結(jié)果。

    (a)Pressure wave template[16] producing sound [35-36]

    (b)Source model of G17HE with measured waveforms [20] acoustic models[25]

    圖3 等效源時(shí)間函數(shù)(頂端黑色曲線)及預(yù)測(cè)(紅色曲線)和實(shí)測(cè)(黑色曲線)波形[14]

    (a)Sound propagation diagram in three-layer medium model

    (b)Yield and height of burst predicted by seismoacoustic analysis

    2 地震波模型

    近地面爆炸耦合的地震波(近距離稱為地沖擊或地運(yùn)動(dòng))來源于爆炸直接耦合方式和空氣超壓(聲波)耦合方式,對(duì)近距離前者稱為直接地沖擊,后者稱為感生地沖擊[12, 38-39],而對(duì)遠(yuǎn)距離,Albert等[40]將前者耦合的地震波稱為前驅(qū)地震波,后者耦合的地震波稱為聲耦合地震波。圖5為近地面爆炸耦合的地震波來源方式示意圖。

    圖5 近地面爆炸耦合的地震波來源方式示意圖[40]

    目前,近地面地震波模型以前驅(qū)地震波模型為主,隨著DM/DB[2],HRI-III[16, 23, 41],SAY[22, 42],FSE和SPE[24]等系列化爆試驗(yàn)的開展,相繼建立了沖積土、沉積巖和花崗巖等場(chǎng)地介質(zhì)的近地面前驅(qū)地震波爆炸當(dāng)量估計(jì)模型[23-24, 26, 42]。其中,Templeton 等[42]基于HRI-III 和SAY 化爆試驗(yàn)地震波數(shù)據(jù)對(duì)地震波模型進(jìn)行了深入討論,對(duì)比和分析了沖積土和沉積巖兩種地質(zhì)成分下的近地面爆炸地震波模型,建立了地震波數(shù)據(jù)分析和模型建立的基本方法,對(duì)近地面化爆地震波建模具有指導(dǎo)意義。圖6為不同場(chǎng)地介質(zhì)的近地面地震波模型和實(shí)測(cè)數(shù)據(jù)。但上述化爆試驗(yàn)場(chǎng)地介質(zhì)類型主要是沖積土、沉積巖和花崗巖等,受試驗(yàn)地質(zhì)成分的限制,地震波模型局限在以上幾種地質(zhì)成分,數(shù)據(jù)還不夠充分,模型還不夠完善[22, 24, 42-43],且未考慮爆源附近地質(zhì)和其余位置地質(zhì)不一樣的情況。

    (a)Alluvium

    (b)Sedimentary rock

    Zhang等[25, 34, 44]基于試驗(yàn)數(shù)據(jù)探討了花崗巖夾雜變質(zhì)巖場(chǎng)地和土石混合場(chǎng)地的近地面地震波實(shí)測(cè)數(shù)據(jù)和模型,分別如圖7和圖8所示。

    (a)Seismic data

    (b)Seismic model

    (a)Model one

    (b)Model two

    (c)Model three

    (d)Model four

    (e)Model five

    為解決爆源附近地質(zhì)和其余位置地質(zhì)不一樣的混合場(chǎng)地地震波模型建模問題,給出了模型的待定系數(shù),并通過引入放大系數(shù)。圖9為放大系數(shù)K不同取值對(duì)反演結(jié)果的影響。圖9中,綠色、紅色虛線分別為聲耦合地震波模型和地震波位移首峰值模型的折中曲線(最小偏差對(duì)應(yīng)的當(dāng)量和埋深),白色五角星為交叉和最小偏差位置。表2為土石混合場(chǎng)地不同模型的定義。

    (a)K=2

    (b)K=4

    (c)K=6

    (d)K=10

    表2 土石混合場(chǎng)地模型的定義[44]

    上述近地面前驅(qū)地震波模型,基于初至波首峰值建立當(dāng)量和爆高(埋深)關(guān)系,適用中近區(qū)爆炸當(dāng)量估計(jì)。對(duì)于遠(yuǎn)區(qū)情況,Pasyanos等[10, 24]基于近地面爆炸地震波能量耦合關(guān)系,采用地震波尾波包絡(luò)法,建立近地面遠(yuǎn)區(qū)地震波模型,進(jìn)一步拓展了近地面地震波模型的適用范圍。圖10為尾波包絡(luò)法示意圖[46]和典型反演結(jié)果[24]。

    (a)Schematic diagram of coda-wave method

    (b)Typical inversion results

    聲耦合地震波模型可解決單一地震計(jì)的近地面爆源當(dāng)量估計(jì)問題,為爆源參數(shù)反演提供除前驅(qū)地震波外的新途徑,但相關(guān)研究較少。Zhang等[46-47]基于聲耦合系數(shù)建立聲耦合地震幅值和聲學(xué)幅值之間的關(guān)系,確定沙土場(chǎng)地的聲震耦合系數(shù),并聯(lián)合近地面聲學(xué)模型建立聲耦合地震波模型。沙場(chǎng)地聲耦合地震波典型波形及耦合系數(shù)衰減特性分別如圖11和圖12所示[46]。表3為聲耦合地震波特征值的定義[46]。

    (a)East

    (b)North

    (c)Verticle

    (a)Z1

    (b)ENZ1 1&2

    (c)ENZ1 0&1

    (d)Z2

    (e)Zhalf

    (f)ENZ2 1&2

    (g)ENZhalf1 1&2

    (h)ENZ2 0&1

    (i)ENZhalf1 0&1

    表3 聲耦合地震波特征值的定義[46]

    3 數(shù)據(jù)融合算法

    近地面爆炸能量會(huì)耦合到大氣和地介質(zhì)中,產(chǎn)生大氣擾動(dòng)(遠(yuǎn)場(chǎng)為聲波)和地運(yùn)動(dòng)(遠(yuǎn)場(chǎng)為地震波)。圖13為近地面爆炸產(chǎn)生的典型聲震信號(hào)[23]。

    (a)Seismo-acoustic energy distribution

    (b)Seismic data of HRI and HRII

    (c)Acoustic data of HRI and HRII

    采用聲學(xué)數(shù)據(jù)或地震波數(shù)據(jù)對(duì)源參數(shù)進(jìn)行反演,當(dāng)量和爆高(埋深)之間通常會(huì)存在明顯的折中關(guān)系使求解不適定。而聲震分析方法通過融合地震波(地運(yùn)動(dòng))數(shù)據(jù)和聲學(xué)(空氣超壓)數(shù)據(jù)使得反演問題變得適定,從而得到當(dāng)量和爆高(埋深)估計(jì)值。圖14為典型試驗(yàn) (FSE-4) 的折中關(guān)系和聲震分析結(jié)果[24]。圖14中,五角星為真實(shí)值,空心三角形為地表處的最優(yōu)點(diǎn),實(shí)心三角形為整個(gè)區(qū)域的最優(yōu)點(diǎn)。

    (a)Results by seismic data

    (b)Results by acoustic data

    (c)Results by seismic and acoustic data

    聲震融合算法是聲震分析的關(guān)鍵組成,隨著聲震分析方法的不斷改進(jìn),聲震融合算法從相對(duì)偏差法[22]、對(duì)數(shù)偏差法[23]及折中曲線交叉法[23]到似然函數(shù)法[24]和貝葉斯MCMC方法[25-26, 48]再到機(jī)器學(xué)習(xí)方法[27],得到了不斷發(fā)展。表4為數(shù)據(jù)融合算法特點(diǎn)和計(jì)算公式。

    表4 數(shù)據(jù)融合算法特點(diǎn)和計(jì)算公式

    圖15為不同數(shù)據(jù)融合方法的典型反演結(jié)果。表4列出的前4種方法沒有考慮模型參數(shù)的先驗(yàn)信息,源參數(shù)缺少進(jìn)一步約束,且計(jì)算采用格點(diǎn)搜索法,效率較低。貝葉斯馬爾可夫鏈蒙特卡羅 (Markov chain Monte Carlo,MCMC)方法[25-26, 48]同時(shí)考慮了偏差信息和先驗(yàn)信息,采用統(tǒng)計(jì)的方法對(duì)爆源參數(shù)進(jìn)行估計(jì)可給出更合理的統(tǒng)計(jì)解釋,同時(shí)采用MCMC方法獲得較高的求解效率。Johannesson等[48]給出了聲震融合的貝葉斯MCMC反演框架,Zhang等[25]和Ford等[26]進(jìn)一步發(fā)展了該方法,基于MCMC方法求解源參數(shù),并應(yīng)用于近地面爆源參數(shù)估計(jì)。圖16為貝葉斯MCMC方法求解流程。以上數(shù)據(jù)融合方法基于已有的聲震模型反演爆源參數(shù)和不同場(chǎng)地巖土類型具有不同的聲震模型。數(shù)據(jù)融合方法的反演精度和聲震模型密切相關(guān)。研究表明[42, 45],聲震模型和場(chǎng)地條件的不匹配將導(dǎo)致較大的反演偏差。為減小因聲震模型匹配程度帶來的反演偏差,提高反演模型在不同場(chǎng)地的適用性,Stroujkova等[27]采用人工神經(jīng)網(wǎng)絡(luò)(artificial neural network,ANN)和貝葉斯正則化方法融合小樣本聲震數(shù)據(jù)來估計(jì)當(dāng)量和爆高(埋深),嘗試將機(jī)器學(xué)習(xí)算法應(yīng)用在爆源參數(shù)反演。圖17為基于深度學(xué)習(xí)的爆源參數(shù)反演流程。

    (a)Method of relative deviation[22]

    (b)Method of logarithmic deviation[23]

    (c)Tradeoff curve crossing method[23]

    (d)Likelihood function method[24]

    (e)Bayesian MCMC method[25, 26, 48]

    (f)Machine learning method[27]

    圖16 貝葉斯MCMC方法求解流程[25]

    圖17 基于深度學(xué)習(xí)的爆源參數(shù)反演流程[27]

    4 總結(jié)與展望

    經(jīng)過十幾年的研究,發(fā)展了多種近地面聲學(xué)模型和地震波模型,適用沖積土、沙土、沉積巖、花崗巖和土石混合體等多種地質(zhì)條件,聯(lián)合相對(duì)偏差法、折中曲線交叉法和貝葉斯MCMC等數(shù)據(jù)融合算法,或基于聲震數(shù)據(jù)采用機(jī)器學(xué)習(xí)方法,可有效解決多種場(chǎng)地近地面爆源的爆炸當(dāng)量和爆高(埋深)反演問題。聲震分析方法雖得到了快速發(fā)展,但還需重點(diǎn)解決以下問題:

    (1)近區(qū)聲/震模型修正完善。受到場(chǎng)地條件限制,不同地質(zhì)、不同爆高和不同距離的聲震數(shù)據(jù)有限,還需更多的試驗(yàn)數(shù)據(jù)對(duì)近地面聲/震模型進(jìn)行補(bǔ)充、修正和完善,從而獲取適用不同地質(zhì)成分的通用模型。目前,聲/震數(shù)據(jù)主要通過外場(chǎng)試驗(yàn)獲取,需進(jìn)一步發(fā)展基于微藥量化爆模擬及超重力離心機(jī)等平臺(tái)的實(shí)驗(yàn)室模擬技術(shù)。

    (2)遠(yuǎn)區(qū)聲/震模型快速建模。區(qū)域較大時(shí),受地質(zhì)地形和大氣環(huán)境等因素的影響,近地面爆炸耦合的地震波和聲波表現(xiàn)出復(fù)雜的傳播特性,給遠(yuǎn)區(qū)爆源參數(shù)快速準(zhǔn)確反演帶來了巨大挑戰(zhàn)。目前,聲震分析技術(shù)以近區(qū)為主,需進(jìn)一步發(fā)展大區(qū)域快速建模技術(shù),實(shí)現(xiàn)遠(yuǎn)區(qū)爆源參數(shù)的快速準(zhǔn)確反演。

    (3)基于機(jī)器學(xué)習(xí)的快速反演方法。機(jī)器學(xué)習(xí)可有效解決數(shù)據(jù)融合和快速計(jì)算問題,但該方法在聲震分析應(yīng)用較少,需基于機(jī)器學(xué)習(xí)算法,進(jìn)一步發(fā)展小樣本聲震數(shù)據(jù)反演技術(shù)及基于物理機(jī)制的大區(qū)域爆源參數(shù)反演技術(shù)。

    (4)未知爆源位置的聲震分析方法。目前聲震分析研究是在爆源位置已知條件下開展的,未考慮爆源位置估計(jì)偏差的影響,需進(jìn)一步開展基于聲震數(shù)據(jù)的爆源位置、當(dāng)量和爆高(埋深)的綜合反演方法研究??赏ㄟ^貝葉斯MCMC或機(jī)器學(xué)習(xí)方法建立爆源多參數(shù)的綜合反演框架。

    猜你喜歡
    爆源當(dāng)量聲學(xué)
    爆源和測(cè)點(diǎn)深度對(duì)水下爆炸沖擊波載荷的影響
    愛的就是這股Hi-Fi味 Davis Acoustics(戴維斯聲學(xué))Balthus 70
    Acoustical Treatment Primer:Diffusion談?wù)劼晫W(xué)處理中的“擴(kuò)散”
    Acoustical Treatment Primer:Absorption談?wù)劼晫W(xué)處理中的“吸聲”(二)
    中型浮動(dòng)沖擊平臺(tái)試驗(yàn)爆源定位方法研究
    Acoustical Treatment Primer:Absorption 談?wù)劼晫W(xué)處理中的“吸聲”
    黃河之聲(2016年24期)2016-02-03 09:01:52
    不同攻角對(duì)實(shí)船爆炸試驗(yàn)沖擊響應(yīng)影響研究
    某型艦抗沖擊試驗(yàn)爆源定位方法研究*
    超壓測(cè)試方法對(duì)炸藥TNT當(dāng)量計(jì)算結(jié)果的影響
    国产免费av片在线观看野外av| 欧美激情高清一区二区三区| 久久精品国产亚洲av高清一级| 69av精品久久久久久 | 午夜免费成人在线视频| 高清视频免费观看一区二区| a级毛片在线看网站| 国产97色在线日韩免费| 欧美乱码精品一区二区三区| 午夜精品国产一区二区电影| 亚洲av美国av| 国产色视频综合| 精品久久蜜臀av无| 免费在线观看视频国产中文字幕亚洲| 国产亚洲精品第一综合不卡| 叶爱在线成人免费视频播放| 国产97色在线日韩免费| 亚洲,欧美精品.| 免费女性裸体啪啪无遮挡网站| 99热国产这里只有精品6| 久久国产精品人妻蜜桃| 精品一区二区三区视频在线观看免费 | 成年女人毛片免费观看观看9 | 97人妻天天添夜夜摸| 啦啦啦免费观看视频1| 国产有黄有色有爽视频| 一区二区av电影网| 丰满少妇做爰视频| 丁香欧美五月| 性色av乱码一区二区三区2| 女警被强在线播放| 亚洲欧美一区二区三区久久| 免费看a级黄色片| 日本黄色日本黄色录像| 99香蕉大伊视频| 亚洲精品久久午夜乱码| 超碰97精品在线观看| 国产欧美日韩一区二区三| 人成视频在线观看免费观看| 日韩中文字幕视频在线看片| 一边摸一边抽搐一进一出视频| 国产成人系列免费观看| 国产精品久久久久久人妻精品电影 | 丰满饥渴人妻一区二区三| h视频一区二区三区| 满18在线观看网站| av又黄又爽大尺度在线免费看| 女人爽到高潮嗷嗷叫在线视频| 久久精品国产a三级三级三级| 久久久国产精品麻豆| 一进一出抽搐动态| 一区二区三区乱码不卡18| 精品人妻1区二区| 成人亚洲精品一区在线观看| 十分钟在线观看高清视频www| 香蕉久久夜色| 国产亚洲精品第一综合不卡| 热99久久久久精品小说推荐| 国产av国产精品国产| 日韩中文字幕视频在线看片| 91大片在线观看| 大码成人一级视频| 亚洲精品乱久久久久久| 中文字幕高清在线视频| videos熟女内射| 国产有黄有色有爽视频| 成人18禁在线播放| 蜜桃国产av成人99| 午夜福利免费观看在线| 精品久久蜜臀av无| 女同久久另类99精品国产91| 久久久久精品人妻al黑| 老司机午夜十八禁免费视频| 免费日韩欧美在线观看| 亚洲视频免费观看视频| 一进一出抽搐动态| 午夜日韩欧美国产| 国产xxxxx性猛交| 交换朋友夫妻互换小说| 午夜激情久久久久久久| 黄色毛片三级朝国网站| 不卡一级毛片| 国产精品一区二区在线观看99| 国产精品秋霞免费鲁丝片| 一区二区三区国产精品乱码| 国产日韩欧美在线精品| 欧美亚洲 丝袜 人妻 在线| 最近最新中文字幕大全电影3 | 侵犯人妻中文字幕一二三四区| av视频免费观看在线观看| 久久天躁狠狠躁夜夜2o2o| 久久久国产成人免费| 欧美久久黑人一区二区| 国产精品1区2区在线观看. | 伊人久久大香线蕉亚洲五| 91成年电影在线观看| 国产精品美女特级片免费视频播放器 | 日韩欧美一区二区三区在线观看 | 久久久精品国产亚洲av高清涩受| 久久精品亚洲精品国产色婷小说| 亚洲精品久久成人aⅴ小说| 最近最新免费中文字幕在线| 久久午夜综合久久蜜桃| 伊人久久大香线蕉亚洲五| 99国产综合亚洲精品| 成人国语在线视频| 一边摸一边做爽爽视频免费| 亚洲av日韩精品久久久久久密| 亚洲伊人久久精品综合| 亚洲色图av天堂| 欧美午夜高清在线| 国产国语露脸激情在线看| 王馨瑶露胸无遮挡在线观看| 日韩制服丝袜自拍偷拍| 精品福利永久在线观看| 亚洲伊人久久精品综合| 午夜激情久久久久久久| 欧美人与性动交α欧美软件| 亚洲国产看品久久| 汤姆久久久久久久影院中文字幕| 大片免费播放器 马上看| 国产亚洲精品第一综合不卡| 欧美亚洲 丝袜 人妻 在线| av国产精品久久久久影院| 亚洲av成人不卡在线观看播放网| 亚洲黑人精品在线| 亚洲 欧美一区二区三区| 日韩欧美国产一区二区入口| 日韩人妻精品一区2区三区| 91老司机精品| 极品少妇高潮喷水抽搐| 一边摸一边抽搐一进一小说 | 少妇粗大呻吟视频| 国产男女超爽视频在线观看| 中文字幕最新亚洲高清| 人人妻,人人澡人人爽秒播| 18禁美女被吸乳视频| 国产成人啪精品午夜网站| 黄色片一级片一级黄色片| 啦啦啦视频在线资源免费观看| 男女午夜视频在线观看| 亚洲一码二码三码区别大吗| 黄色a级毛片大全视频| 丝瓜视频免费看黄片| 天堂中文最新版在线下载| 高清av免费在线| 在线观看www视频免费| 女人久久www免费人成看片| 成人国语在线视频| 老司机靠b影院| 黄色视频,在线免费观看| 午夜久久久在线观看| 精品一区二区三区视频在线观看免费 | 久久国产亚洲av麻豆专区| 国产成人系列免费观看| 国产熟女午夜一区二区三区| 老熟妇乱子伦视频在线观看| 日韩熟女老妇一区二区性免费视频| 日韩中文字幕视频在线看片| 亚洲精品中文字幕在线视频| 久久精品国产亚洲av高清一级| 国产人伦9x9x在线观看| 免费黄频网站在线观看国产| 免费观看a级毛片全部| 日韩免费高清中文字幕av| 99国产极品粉嫩在线观看| 99国产精品一区二区蜜桃av | 久久婷婷成人综合色麻豆| 国产av精品麻豆| 欧美中文综合在线视频| 亚洲精品在线观看二区| 青草久久国产| 不卡av一区二区三区| 亚洲视频免费观看视频| 亚洲人成伊人成综合网2020| 国产日韩欧美视频二区| 免费观看av网站的网址| 国产精品一区二区在线不卡| 亚洲第一av免费看| 国产国语露脸激情在线看| a级毛片黄视频| 亚洲,欧美精品.| 麻豆国产av国片精品| 久久精品亚洲精品国产色婷小说| 99在线人妻在线中文字幕 | 夜夜骑夜夜射夜夜干| 国产一卡二卡三卡精品| 少妇的丰满在线观看| 建设人人有责人人尽责人人享有的| 90打野战视频偷拍视频| 久久国产精品人妻蜜桃| 精品视频人人做人人爽| 99久久精品国产亚洲精品| 国产精品一区二区精品视频观看| 岛国毛片在线播放| 国产亚洲午夜精品一区二区久久| 飞空精品影院首页| 日韩人妻精品一区2区三区| 少妇猛男粗大的猛烈进出视频| 欧美av亚洲av综合av国产av| 最黄视频免费看| 女性被躁到高潮视频| 中文字幕人妻丝袜一区二区| 国产激情久久老熟女| 12—13女人毛片做爰片一| 性色av乱码一区二区三区2| 9色porny在线观看| 国产不卡一卡二| 精品国产乱码久久久久久小说| 老熟妇乱子伦视频在线观看| 国产精品av久久久久免费| 国产精品一区二区在线不卡| 日韩 欧美 亚洲 中文字幕| 午夜精品久久久久久毛片777| 热99国产精品久久久久久7| 亚洲av美国av| 亚洲第一欧美日韩一区二区三区 | a级毛片在线看网站| 热re99久久国产66热| 亚洲精品久久成人aⅴ小说| 人人妻人人澡人人看| 成年版毛片免费区| 老司机影院毛片| 国产片内射在线| 国产男靠女视频免费网站| 亚洲av国产av综合av卡| 欧美激情久久久久久爽电影 | 曰老女人黄片| 人人妻人人澡人人爽人人夜夜| 国产高清videossex| 人人妻,人人澡人人爽秒播| 汤姆久久久久久久影院中文字幕| 国产精品国产高清国产av | 精品国产一区二区久久| 热re99久久精品国产66热6| 一夜夜www| 啦啦啦中文免费视频观看日本| 巨乳人妻的诱惑在线观看| 午夜精品久久久久久毛片777| 亚洲色图av天堂| 精品久久久久久久毛片微露脸| 伦理电影免费视频| 丝袜喷水一区| 国产成人一区二区三区免费视频网站| 亚洲三区欧美一区| 久久久久久久久久久久大奶| 欧美人与性动交α欧美软件| 黑丝袜美女国产一区| 50天的宝宝边吃奶边哭怎么回事| 他把我摸到了高潮在线观看 | 一进一出好大好爽视频| 欧美另类亚洲清纯唯美| 亚洲成人免费av在线播放| 成人手机av| 老司机亚洲免费影院| 19禁男女啪啪无遮挡网站| 一级毛片电影观看| 精品一区二区三卡| 中文字幕色久视频| 99国产极品粉嫩在线观看| 欧美成狂野欧美在线观看| 女性被躁到高潮视频| 色94色欧美一区二区| 欧美国产精品一级二级三级| 国产又色又爽无遮挡免费看| 18禁裸乳无遮挡动漫免费视频| 少妇粗大呻吟视频| 亚洲色图 男人天堂 中文字幕| 亚洲一区中文字幕在线| 亚洲欧美一区二区三区黑人| 欧美乱码精品一区二区三区| 露出奶头的视频| 最新美女视频免费是黄的| 99re在线观看精品视频| 日韩中文字幕视频在线看片| 亚洲成人免费av在线播放| 99香蕉大伊视频| 国产精品亚洲一级av第二区| 夜夜爽天天搞| 岛国在线观看网站| 久久人妻熟女aⅴ| 少妇猛男粗大的猛烈进出视频| 国产精品免费一区二区三区在线 | 精品少妇久久久久久888优播| 欧美变态另类bdsm刘玥| 国产男女内射视频| 巨乳人妻的诱惑在线观看| aaaaa片日本免费| 妹子高潮喷水视频| 欧美日韩成人在线一区二区| 麻豆乱淫一区二区| 午夜福利在线观看吧| 中文字幕精品免费在线观看视频| 久久国产精品大桥未久av| 国产片内射在线| 国产亚洲精品第一综合不卡| 日本vs欧美在线观看视频| 老汉色∧v一级毛片| 这个男人来自地球电影免费观看| 变态另类成人亚洲欧美熟女 | 热re99久久精品国产66热6| 日韩人妻精品一区2区三区| av不卡在线播放| 免费高清在线观看日韩| 亚洲中文日韩欧美视频| 91成年电影在线观看| 免费在线观看日本一区| 日本黄色日本黄色录像| av国产精品久久久久影院| a级片在线免费高清观看视频| 熟女少妇亚洲综合色aaa.| 国产男靠女视频免费网站| 欧美日韩视频精品一区| 动漫黄色视频在线观看| 大香蕉久久成人网| 建设人人有责人人尽责人人享有的| 电影成人av| 岛国在线观看网站| 欧美精品一区二区大全| 国产精品国产高清国产av | 亚洲综合色网址| 99国产极品粉嫩在线观看| 国产极品粉嫩免费观看在线| 纯流量卡能插随身wifi吗| 如日韩欧美国产精品一区二区三区| 一二三四在线观看免费中文在| 国产有黄有色有爽视频| 欧美人与性动交α欧美软件| 51午夜福利影视在线观看| 99九九在线精品视频| 日韩熟女老妇一区二区性免费视频| 久久99热这里只频精品6学生| av视频免费观看在线观看| 午夜久久久在线观看| 国产淫语在线视频| 欧美国产精品一级二级三级| 久久狼人影院| √禁漫天堂资源中文www| 99re在线观看精品视频| 国产精品久久电影中文字幕 | 精品少妇黑人巨大在线播放| 黄色a级毛片大全视频| 18禁国产床啪视频网站| 亚洲情色 制服丝袜| 久久天堂一区二区三区四区| www.熟女人妻精品国产| 欧美老熟妇乱子伦牲交| 91大片在线观看| 国产欧美日韩一区二区三| 成人18禁在线播放| 久久中文看片网| 久久ye,这里只有精品| 国产高清激情床上av| 国产激情久久老熟女| 国产又色又爽无遮挡免费看| 午夜成年电影在线免费观看| kizo精华| 50天的宝宝边吃奶边哭怎么回事| 黄色视频不卡| 亚洲精品av麻豆狂野| 亚洲成人国产一区在线观看| 天堂8中文在线网| 国产欧美日韩一区二区三| 国产高清videossex| 一进一出好大好爽视频| 天堂8中文在线网| 人人妻人人爽人人添夜夜欢视频| 精品国内亚洲2022精品成人 | 国产免费福利视频在线观看| 老熟妇乱子伦视频在线观看| 丝袜人妻中文字幕| 久久久久久免费高清国产稀缺| 欧美成狂野欧美在线观看| 欧美日韩黄片免| 亚洲天堂av无毛| 午夜两性在线视频| 国产欧美日韩一区二区精品| 好男人电影高清在线观看| 日本wwww免费看| 伊人久久大香线蕉亚洲五| 国产精品秋霞免费鲁丝片| 久久精品国产a三级三级三级| 亚洲伊人久久精品综合| 无遮挡黄片免费观看| 日韩大片免费观看网站| 久久精品熟女亚洲av麻豆精品| 人妻久久中文字幕网| 精品国产一区二区三区四区第35| 在线播放国产精品三级| 亚洲av美国av| 国产免费视频播放在线视频| 精品国产乱码久久久久久男人| 老司机深夜福利视频在线观看| 男女免费视频国产| 亚洲五月色婷婷综合| 国产黄色免费在线视频| 婷婷成人精品国产| 女警被强在线播放| 丝袜喷水一区| 国产精品 欧美亚洲| 热99久久久久精品小说推荐| 高潮久久久久久久久久久不卡| 91精品三级在线观看| 看免费av毛片| 老司机影院毛片| 免费在线观看视频国产中文字幕亚洲| 国产成人精品久久二区二区91| 欧美日韩精品网址| 国产精品偷伦视频观看了| 999久久久精品免费观看国产| 久久香蕉激情| 国产野战对白在线观看| 大片电影免费在线观看免费| 亚洲成人国产一区在线观看| 一级毛片精品| 大陆偷拍与自拍| 一本综合久久免费| 69av精品久久久久久 | 三级毛片av免费| 精品亚洲成a人片在线观看| 一本—道久久a久久精品蜜桃钙片| 99国产极品粉嫩在线观看| 18禁裸乳无遮挡动漫免费视频| 欧美成狂野欧美在线观看| 久久久久久免费高清国产稀缺| 午夜精品久久久久久毛片777| www.自偷自拍.com| 中文字幕色久视频| 欧美午夜高清在线| 啦啦啦视频在线资源免费观看| 亚洲性夜色夜夜综合| 欧美性长视频在线观看| 亚洲精华国产精华精| 午夜两性在线视频| 人人妻人人澡人人看| 波多野结衣av一区二区av| 色综合欧美亚洲国产小说| 男女无遮挡免费网站观看| 久久ye,这里只有精品| 国产欧美日韩精品亚洲av| 国产1区2区3区精品| 国产有黄有色有爽视频| 亚洲欧美精品综合一区二区三区| 久久九九热精品免费| 女人被躁到高潮嗷嗷叫费观| 久久久久久亚洲精品国产蜜桃av| 久久久精品免费免费高清| 青青草视频在线视频观看| 亚洲精品国产精品久久久不卡| 午夜91福利影院| 在线看a的网站| 国产精品免费视频内射| 美女主播在线视频| 国产精品麻豆人妻色哟哟久久| 午夜久久久在线观看| 亚洲av美国av| 9191精品国产免费久久| 电影成人av| 亚洲精品中文字幕在线视频| 两人在一起打扑克的视频| 国产不卡一卡二| 成人影院久久| av有码第一页| 五月天丁香电影| 视频在线观看一区二区三区| 老熟妇仑乱视频hdxx| 成人国产av品久久久| 一级片'在线观看视频| 久久毛片免费看一区二区三区| 色视频在线一区二区三区| 一个人免费在线观看的高清视频| 9热在线视频观看99| 大香蕉久久成人网| 97在线人人人人妻| 在线播放国产精品三级| 国产精品.久久久| 超碰成人久久| 十八禁人妻一区二区| 久久久久久久大尺度免费视频| 亚洲成a人片在线一区二区| 久久精品aⅴ一区二区三区四区| 我的亚洲天堂| 国产有黄有色有爽视频| 国产日韩欧美视频二区| 极品人妻少妇av视频| 男女无遮挡免费网站观看| 天堂8中文在线网| 亚洲性夜色夜夜综合| 777久久人妻少妇嫩草av网站| videosex国产| 一边摸一边抽搐一进一出视频| 成人手机av| 国产日韩欧美亚洲二区| www.熟女人妻精品国产| 午夜福利免费观看在线| 国产精品一区二区精品视频观看| 日韩有码中文字幕| 久久狼人影院| 国产精品电影一区二区三区 | 精品欧美一区二区三区在线| 一本一本久久a久久精品综合妖精| 99久久99久久久精品蜜桃| 久久亚洲真实| 午夜视频精品福利| 高清av免费在线| 久久久久精品人妻al黑| 免费日韩欧美在线观看| 亚洲熟女精品中文字幕| 天天影视国产精品| 我的亚洲天堂| 久久精品国产亚洲av高清一级| 久久久久久免费高清国产稀缺| 女人精品久久久久毛片| 国产精品免费一区二区三区在线 | 亚洲专区字幕在线| 久久精品亚洲av国产电影网| 午夜福利视频在线观看免费| 91大片在线观看| 大陆偷拍与自拍| 女同久久另类99精品国产91| 亚洲伊人久久精品综合| 日韩欧美一区视频在线观看| 日本撒尿小便嘘嘘汇集6| 精品一区二区三区视频在线观看免费 | 在线观看人妻少妇| 久久久久精品国产欧美久久久| 色94色欧美一区二区| 国产精品一区二区精品视频观看| 黄色片一级片一级黄色片| 精品一品国产午夜福利视频| 国精品久久久久久国模美| 大型av网站在线播放| 999精品在线视频| 精品欧美一区二区三区在线| 亚洲精品美女久久av网站| 不卡av一区二区三区| 亚洲色图 男人天堂 中文字幕| 亚洲欧美色中文字幕在线| 久久人妻av系列| 中文字幕另类日韩欧美亚洲嫩草| 丰满饥渴人妻一区二区三| 国产精品亚洲av一区麻豆| 一本一本久久a久久精品综合妖精| 久久精品亚洲精品国产色婷小说| 国产无遮挡羞羞视频在线观看| 国产精品亚洲av一区麻豆| 人人妻人人爽人人添夜夜欢视频| 两性午夜刺激爽爽歪歪视频在线观看 | 新久久久久国产一级毛片| 91大片在线观看| 日韩欧美免费精品| 一边摸一边做爽爽视频免费| 成人国产av品久久久| 50天的宝宝边吃奶边哭怎么回事| 亚洲欧美一区二区三区黑人| 国产高清视频在线播放一区| 国产成人啪精品午夜网站| 91成年电影在线观看| 国产视频一区二区在线看| 色婷婷av一区二区三区视频| 欧美黄色片欧美黄色片| 不卡av一区二区三区| 99精品久久久久人妻精品| 欧美黑人精品巨大| 热re99久久国产66热| 汤姆久久久久久久影院中文字幕| 免费黄频网站在线观看国产| 免费在线观看日本一区| 母亲3免费完整高清在线观看| 黑人巨大精品欧美一区二区蜜桃| 正在播放国产对白刺激| 日韩欧美免费精品| 成人av一区二区三区在线看| 美女扒开内裤让男人捅视频| 国产精品一区二区在线观看99| 中文亚洲av片在线观看爽 | 老熟妇乱子伦视频在线观看| 国产极品粉嫩免费观看在线| 欧美午夜高清在线| 欧美 亚洲 国产 日韩一| 后天国语完整版免费观看| 日韩视频在线欧美| 真人做人爱边吃奶动态| 久久久久久久久久久久大奶| 亚洲少妇的诱惑av| 怎么达到女性高潮| 老鸭窝网址在线观看| 国产男靠女视频免费网站| 我的亚洲天堂| 亚洲国产欧美一区二区综合| 亚洲伊人色综图| 国产精品一区二区免费欧美| 纵有疾风起免费观看全集完整版| av线在线观看网站| 黄色a级毛片大全视频| 啪啪无遮挡十八禁网站| 国产男靠女视频免费网站| 亚洲精品久久午夜乱码| 最近最新免费中文字幕在线| 免费黄频网站在线观看国产| 日日爽夜夜爽网站| 十八禁人妻一区二区| 日韩视频在线欧美| 亚洲午夜理论影院| 日本av免费视频播放| 国产主播在线观看一区二区| netflix在线观看网站| 日本av免费视频播放| 久久精品aⅴ一区二区三区四区| 国产亚洲精品一区二区www | 91精品国产国语对白视频| 超色免费av| 日韩精品免费视频一区二区三区| 亚洲av片天天在线观看| 成人手机av|