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

    瞬時(shí)與累積γ劑量率測(cè)量數(shù)據(jù)統(tǒng)計(jì)學(xué)檢驗(yàn)方法

    2022-07-14 11:39:14王宇健陳志祥
    關(guān)鍵詞:劑量率正態(tài)分布點(diǎn)位

    王宇健 陳志祥

    (1.山東省煙臺(tái)生態(tài)環(huán)境監(jiān)測(cè)中心,山東 煙臺(tái) 264001;2.山東省核與輻射安全監(jiān)測(cè)中心,山東 濟(jì)南 250000)

    0 引言

    環(huán)境γ 劑量率監(jiān)測(cè)在輻射環(huán)境監(jiān)測(cè)領(lǐng)域有重要意義,可用于積累天然放射性本底水平數(shù)據(jù),為異常情況和輻射應(yīng)急決策提供參考;監(jiān)測(cè)輻射源放射性水平、了解γ輻射變化趨勢(shì),評(píng)估輻射源的影響。瞬時(shí)γ 劑量率監(jiān)測(cè)、累積γ 劑量率監(jiān)測(cè)是環(huán)境γ 劑量率監(jiān)測(cè)常用的監(jiān)測(cè)方式,兩種方法各有特點(diǎn),互為補(bǔ)充。

    目前,對(duì)于輻射環(huán)境γ 劑量率測(cè)量結(jié)果準(zhǔn)確性的保證多側(cè)重于測(cè)量過(guò)程的質(zhì)量控制方面,包括規(guī)范儀器操作使用方法、按時(shí)檢定、繪制質(zhì)控圖、儀器比對(duì)等方法以保證儀器測(cè)量穩(wěn)定性、準(zhǔn)確性,以上方法可保證測(cè)量結(jié)果的可信度。對(duì)于測(cè)得的數(shù)據(jù)本身是否在合理水平缺少必要的甄別手段,常通過(guò)簡(jiǎn)單的數(shù)據(jù)對(duì)比,利用技術(shù)人員的經(jīng)驗(yàn)判定,這種基于簡(jiǎn)單對(duì)比和經(jīng)驗(yàn)判斷的評(píng)價(jià)方式,說(shuō)服力相對(duì)有限。也有人采用3σ 判定法則進(jìn)行異常值判定,該方法使用較為簡(jiǎn)單,但檢驗(yàn)效果較差。所以,找到一種基于統(tǒng)計(jì)學(xué)的,通過(guò)定量計(jì)算對(duì)γ 劑量率(瞬時(shí)、TLD)測(cè)量結(jié)果進(jìn)行判定的合理檢驗(yàn)方法,有助于提升數(shù)據(jù)質(zhì)量和輻射環(huán)境水平描述的合理性、說(shuō)服性。

    1 統(tǒng)計(jì)分布假設(shè)與檢驗(yàn)方法的選擇

    1.1 統(tǒng)計(jì)分布假設(shè)的選擇

    環(huán)境統(tǒng)計(jì)學(xué)作為統(tǒng)計(jì)學(xué)在環(huán)境領(lǐng)域的應(yīng)用,要求參與統(tǒng)計(jì)的測(cè)量數(shù)據(jù)能夠代表總體環(huán)境水平,可以通過(guò)合理的假設(shè)找到一種適合的統(tǒng)計(jì)學(xué)分布并加以檢驗(yàn)。按照以上原則,該文對(duì)試驗(yàn)條件設(shè)置如下:1)參考《環(huán)境空氣質(zhì)量監(jiān)測(cè)點(diǎn)位布設(shè)技術(shù)規(guī)范(試行)》(HJ 664-2013)和TOPSIS方法,在半徑約10km 范圍內(nèi)各布設(shè)6 個(gè)γ 劑量率瞬時(shí)和TLD 測(cè)量點(diǎn)位,且實(shí)驗(yàn)時(shí)間段內(nèi)點(diǎn)位固定不變,以保證區(qū)域覆蓋性和代表性。2)測(cè)量全過(guò)程、各要素嚴(yán)格執(zhí)行相關(guān)標(biāo)準(zhǔn)及質(zhì)控要求,確保數(shù)據(jù)真實(shí)、可靠,在理論上保證測(cè)量數(shù)據(jù)—檢驗(yàn)結(jié)果—輻射環(huán)境水平三者的一致性。3)有一定的測(cè)量頻次,且測(cè)量時(shí)間相對(duì)固定,在一定時(shí)間段內(nèi)可積累足夠數(shù)量并相對(duì)連續(xù)的數(shù)據(jù)以供統(tǒng)計(jì)分析,保證統(tǒng)計(jì)分布假設(shè)與檢驗(yàn)的準(zhǔn)確性;四是γ 劑量率瞬時(shí)和TLD 測(cè)量同點(diǎn)位、同頻次進(jìn)行,且數(shù)據(jù)經(jīng)過(guò)計(jì)算轉(zhuǎn)換均以nGy/h 為單位,作為相對(duì)獨(dú)立,但理論上又能反映同一對(duì)象及其變化趨勢(shì)的兩種樣本數(shù)據(jù),其樣本間比較結(jié)果具有一定的說(shuō)明性和驗(yàn)證性;五是本底情況下,環(huán)境γ 劑量率瞬時(shí)與TLD 測(cè)量結(jié)果一般相對(duì)穩(wěn)定,出現(xiàn)離群值或整體水平較大波動(dòng)的原因較易辨識(shí)。

    同時(shí),作為環(huán)境級(jí)測(cè)量數(shù)據(jù),γ 劑量率測(cè)量受多種因素影響,在測(cè)量過(guò)程中也應(yīng)注意和盡量避免以下問(wèn)題對(duì)測(cè)量結(jié)果及統(tǒng)計(jì)分布的影響:1)TLD 測(cè)量過(guò)程中重復(fù)測(cè)量引入的A 類不確定度以及刻度因子引入的B 類不確定度,以及瞬時(shí)測(cè)量時(shí),遇到施工、修路等本底環(huán)境臨時(shí)或永久改變等;2)有研究表明,TLD 測(cè)量結(jié)果在一定程度上受環(huán)境溫、濕度影響,具體影響情況為,由于低溫峰和主計(jì)量峰部分重疊,熱釋光測(cè)量?jī)x在計(jì)數(shù)前的“退火”階段,難以完全排除低溫峰對(duì)測(cè)量計(jì)數(shù)的影響,從而難以排除布片過(guò)程中環(huán)境溫度對(duì)最終測(cè)量結(jié)果的影響,另外長(zhǎng)時(shí)間連續(xù)降水導(dǎo)致的環(huán)境濕度過(guò)高對(duì)劑量片吸收情況也有影響;3)一般情況下,瞬時(shí)γ 劑量率的測(cè)量結(jié)果減去了宇宙射線響應(yīng)值,而TLD 難以準(zhǔn)確測(cè)量宇宙射線響應(yīng)值而未予去除,因此兩種測(cè)量結(jié)果存在一定差值(正常情況下瞬時(shí)γ 劑量率小于TLD 測(cè)量數(shù)據(jù)),且差值難以準(zhǔn)確計(jì)量,一定程度上影響二者比對(duì);4)TLD 測(cè)量數(shù)據(jù)反映一個(gè)時(shí)段γ 劑量率的平均值,瞬時(shí)測(cè)量數(shù)據(jù)反映測(cè)量當(dāng)時(shí)的γ 劑量率水平,二者趨勢(shì)一致,但單點(diǎn)位、單次測(cè)量值對(duì)比結(jié)果僅供參考。

    綜上所述,在做好測(cè)量過(guò)程質(zhì)量控制的基礎(chǔ)上,可暫時(shí)忽略上述不利因素對(duì)統(tǒng)計(jì)分布假設(shè)與檢驗(yàn)的影響,而認(rèn)為輻射環(huán)境γ 劑量率瞬時(shí)和TLD 測(cè)量數(shù)據(jù)具備準(zhǔn)確、穩(wěn)定、連續(xù)、受環(huán)境本身非輻射水平因素影響較少,且可以積累足夠數(shù)量的數(shù)據(jù)。所以,其數(shù)據(jù)作為樣本,具有一定的代表性,可視為連續(xù)隨機(jī)變量,從而對(duì)其進(jìn)行正態(tài)分布假設(shè),并對(duì)該假設(shè)進(jìn)行檢驗(yàn)。

    1.2 檢驗(yàn)方法的選擇

    正態(tài)分布是統(tǒng)計(jì)分析中最重要的分布,檢驗(yàn)一組數(shù)據(jù)是否服從正態(tài)分布通常有3 種方法:計(jì)算偏度系數(shù)和峰度系數(shù)檢驗(yàn);夏皮洛-威爾克(S-W)正態(tài)分布檢驗(yàn)(樣本容量較小時(shí));P-P 圖及Q-Q 圖圖形檢驗(yàn)。該文選擇前兩種定量分析方法進(jìn)行數(shù)據(jù)分析,兩種方法的檢驗(yàn)結(jié)果在表達(dá)正態(tài)分布不同特性方面各有側(cè)重,為更全面分析試驗(yàn)數(shù)據(jù)的特性,該文同時(shí)選用2 種方法進(jìn)行檢驗(yàn)。

    2 統(tǒng)計(jì)分布假設(shè)檢驗(yàn)過(guò)程與結(jié)果分析

    該文結(jié)合試驗(yàn)數(shù)據(jù)(瞬時(shí)、TLD 各6 組共址測(cè)量數(shù)據(jù),實(shí)驗(yàn)時(shí)間2017 年三季度~2021 年二季度,測(cè)量頻次一季度/次,不考慮缺失值情況下單組樣本數(shù)量16 個(gè)),重點(diǎn)研究統(tǒng)計(jì)檢驗(yàn)方法在γ 劑量率瞬時(shí)和TLD 測(cè)量結(jié)果分析中的應(yīng)用。監(jiān)測(cè)數(shù)據(jù)中出現(xiàn)個(gè)別缺失值,是環(huán)境監(jiān)測(cè)及環(huán)境統(tǒng)計(jì)學(xué)中經(jīng)常遇到的情況。尤其在TLD 測(cè)量中,劑量片經(jīng)常會(huì)因?yàn)槠洳荚O(shè)樹(shù)木被修剪砍伐而導(dǎo)致數(shù)據(jù)缺失。

    在某些環(huán)境水平分析應(yīng)用中,為更好體現(xiàn)整體環(huán)境的連續(xù)水平,且考慮測(cè)量數(shù)據(jù)具有序列特征,常對(duì)原始測(cè)量數(shù)據(jù)中一組按照時(shí)間排列樣本中的缺失值進(jìn)行統(tǒng)計(jì)學(xué)補(bǔ)全:對(duì)于單個(gè)缺失值,可使用臨近點(diǎn)均值插補(bǔ),臨近值個(gè)數(shù)一般選1~3;若測(cè)量數(shù)據(jù)中出現(xiàn)連續(xù)2 個(gè)缺失值,可進(jìn)行線性插值插補(bǔ);若單組樣本首、末位序數(shù)樣本出現(xiàn)缺失值,可進(jìn)行臨近點(diǎn)線性趨勢(shì)插補(bǔ)。

    需要注意的是,統(tǒng)計(jì)學(xué)補(bǔ)全的數(shù)據(jù)相當(dāng)于人為對(duì)該組數(shù)據(jù)的分布特性加入了一定程度的干擾,會(huì)讓該組數(shù)據(jù)更加符合某種統(tǒng)計(jì)學(xué)分布,其檢驗(yàn)結(jié)果更傾向于統(tǒng)計(jì)學(xué)預(yù)期。為更直接地說(shuō)明正態(tài)分布假設(shè)及其檢驗(yàn)方法在該試驗(yàn)中的應(yīng)用,所以不對(duì)缺失值進(jìn)行補(bǔ)全處理。原始數(shù)據(jù)見(jiàn)表1。

    2.1 正態(tài)分布檢驗(yàn)方法與結(jié)果

    使用以下2 種方法對(duì)表1 中γ 劑量率瞬時(shí)與TLD 測(cè)量數(shù)據(jù)的正態(tài)分布假設(shè)進(jìn)行檢驗(yàn)。

    計(jì)算偏度和峰度系數(shù)定量檢驗(yàn)。偏度是統(tǒng)計(jì)數(shù)據(jù)分布偏斜方向和程度的度量,定義為樣本的三階標(biāo)準(zhǔn)化距,樣本偏度系數(shù)如公式(1)所示。

    峰度是研究數(shù)據(jù)分布陡峭或平滑的統(tǒng)計(jì)量,定義為四階中心距除以方差的平方減3,樣本峰度系數(shù)如公式(2)所示。

    以Z 評(píng)分判定其是否服從正態(tài)分布,Z=偏度系數(shù)or 峰度系數(shù)/偏度系數(shù)標(biāo)準(zhǔn)誤,在α=0.05 的檢驗(yàn)水平下,如果Z在±1.96 之間,則可以說(shuō)明數(shù)據(jù)服從正態(tài)分布。

    因公式中各參量可直接讀取或通過(guò)簡(jiǎn)單計(jì)算獲取,不再贅述計(jì)算過(guò)程,計(jì)算結(jié)果見(jiàn)表2。

    夏皮洛-威爾克(S-W)正態(tài)分布檢驗(yàn)。S-W 屬于非參數(shù)檢驗(yàn),用來(lái)檢驗(yàn)數(shù)據(jù)是否符合正態(tài)分布,樣本容量較?。颖救萘?-50)時(shí)適用。具體檢驗(yàn)步驟如公式(3)所示。

    該文中樣本數(shù)據(jù)假設(shè)檢驗(yàn)顯著性水平α 取0.05,選取點(diǎn)位1 瞬時(shí)γ 劑量率為例進(jìn)行假設(shè)檢驗(yàn)。

    作出假設(shè)H:點(diǎn)位1 瞬時(shí)γ 劑量率數(shù)據(jù)服從正態(tài)分布,備擇假設(shè)為H。

    a:將歷次測(cè)量瞬時(shí)γ 劑量率數(shù)據(jù)從小到大排列成:72.4,73.1,74.4,74.5,77.5,79.5,79.9,80.7,82.7,83.1,85.6,86.0,86.9,93.1,93.6,97.5。

    b:查閱GB4882 中a值表a~a分別為0.5056,0.3290,0.2521,0.1939,0.1447,0.1005,0.0593,0.0196,同時(shí)計(jì)算樣本均值,帶入公式,求得W=0.946。

    c:對(duì)α=0.05,查表得n=16 時(shí),Z=0.887。

    d:由于0.946 >0.887,所以不拒絕正態(tài)性原假設(shè)H0。

    所有點(diǎn)位瞬時(shí)及累積γ 劑量率S-W 檢驗(yàn)結(jié)果見(jiàn)表2。

    正態(tài)分布假設(shè)的檢驗(yàn)結(jié)果。從上述兩種正態(tài)分布定量檢驗(yàn)方法對(duì)γ 劑量率瞬時(shí)和TLD 測(cè)量各6 組樣本數(shù)據(jù)的檢驗(yàn)情況來(lái)看,兩種方法的檢驗(yàn)結(jié)果一致度較高,即:瞬時(shí)測(cè)量全部6 組樣本數(shù)據(jù)及TLD 測(cè)量點(diǎn)位1、2、3、6 四組樣本數(shù)據(jù)服從正態(tài)分布假設(shè),樣本數(shù)據(jù)其偏(峰)度系數(shù)/偏(峰)度系數(shù)標(biāo)準(zhǔn)誤的絕對(duì)值均小于1.96 的同時(shí),其夏皮洛-威爾克(S-W)檢驗(yàn)結(jié)果也服從正態(tài)分布;TLD 測(cè)量點(diǎn)位4、5兩組樣本數(shù)據(jù)不服從正態(tài)分布假設(shè),其中點(diǎn)位4、5 偏度系數(shù)/偏度系數(shù)標(biāo)準(zhǔn)誤的絕對(duì)值均大于1.96,同時(shí)點(diǎn)位5 夏皮洛-威爾克(S-W)檢驗(yàn)W 值小于臨界值。

    另外不難看出,雖然TLD 測(cè)量點(diǎn)位1、2、3、6 四組數(shù)據(jù)樣本總體服從正態(tài)分布,但也存在一定的“離散”傾向,主要體現(xiàn)在其偏度系數(shù)/偏度系數(shù)標(biāo)準(zhǔn)誤的絕對(duì)值比較接近1.96、夏皮洛-威爾克(S-W)檢驗(yàn)結(jié)果比較接近臨界值。對(duì)于不符合正態(tài)分布和有一定程度“離散”傾向的樣本數(shù)據(jù)組,均應(yīng)懷疑其存在離群值。

    2.2 數(shù)據(jù)分析

    對(duì)TLD 數(shù)據(jù)出現(xiàn)的不符合正態(tài)分布檢驗(yàn)的結(jié)果進(jìn)行追溯和分析發(fā)現(xiàn),TLD 測(cè)量過(guò)程中TLD 測(cè)量設(shè)備進(jìn)行過(guò)多次校準(zhǔn),以2018 年1 季度為界,在此之前,各點(diǎn)位TLD 數(shù)據(jù)整體高于后續(xù)組內(nèi)其他數(shù)據(jù)。經(jīng)分析TLD 測(cè)量設(shè)備經(jīng)校準(zhǔn)后刻度因子較高,導(dǎo)致了測(cè)量數(shù)據(jù)整體正態(tài)性變差,這也是TLD 各點(diǎn)位出現(xiàn)檢驗(yàn)結(jié)果不符合正態(tài)分布或接近臨界值的根本原因,2018 年1 季度重新檢定并使用新的刻度因子后,測(cè)量數(shù)據(jù)整體穩(wěn)定。

    表1 γ 劑量率瞬時(shí)與TLD 測(cè)量數(shù)據(jù)(單位:nGy/h)

    表2 正態(tài)分布假設(shè)檢驗(yàn)結(jié)果

    雖然γ 劑量率瞬時(shí)與TLD 測(cè)量有諸多不同之處,外界環(huán)境影響以及儀器設(shè)備本身穩(wěn)定性差異必然導(dǎo)致最終測(cè)量結(jié)果的不同,但二者在同點(diǎn)位較長(zhǎng)時(shí)間段內(nèi)測(cè)量數(shù)據(jù)的變化趨勢(shì)總體保持一致。

    3 結(jié)論

    該文對(duì)試驗(yàn)區(qū)域內(nèi)(同點(diǎn)位、同頻次)環(huán)境γ 輻射空氣吸收劑量率瞬時(shí)和TLD 測(cè)量結(jié)果進(jìn)行了正態(tài)分布假設(shè)與兩種不同方法檢驗(yàn),檢驗(yàn)結(jié)果表明,該試驗(yàn)數(shù)據(jù)總體質(zhì)量良好,其中瞬時(shí)測(cè)量數(shù)據(jù)的穩(wěn)定性優(yōu)于TLD 測(cè)量數(shù)據(jù),試驗(yàn)區(qū)域內(nèi)γ 劑量率輻射水平發(fā)展趨勢(shì)在合理區(qū)間內(nèi)正常波動(dòng)。另外,可通過(guò)繪制瞬時(shí)與累積測(cè)量數(shù)據(jù)趨勢(shì)圖的方式,進(jìn)一步驗(yàn)證二者測(cè)量結(jié)果的合理性。

    該文選用的相關(guān)假設(shè)、檢驗(yàn)及相應(yīng)的數(shù)據(jù)追溯方法,可以作為對(duì)環(huán)境γ 劑量率瞬時(shí)、TLD 測(cè)量數(shù)據(jù)質(zhì)量定量評(píng)價(jià),以及環(huán)境γ 劑量率輻射水平描述的有效手段。通過(guò)對(duì)數(shù)據(jù)本身的統(tǒng)計(jì)假設(shè)檢驗(yàn),也對(duì)儀器設(shè)備的穩(wěn)定性,運(yùn)行狀況提供了參考資料。

    兩種不同的假設(shè)檢驗(yàn)方法不僅適用于環(huán)境γ 劑量率數(shù)據(jù)分析判定,也可在輻射環(huán)境監(jiān)測(cè)的其他項(xiàng)目得以應(yīng)用,特別是正常情況下監(jiān)測(cè)結(jié)果穩(wěn)定在探測(cè)下限(MDC)以上的環(huán)境樣品,例如核電基地周邊環(huán)境介質(zhì)中總放射性分析、環(huán)境介質(zhì)中氚、碳-14 的分析等。隨著核電廠放射性流出物的排放,采用兩種不同的假設(shè)檢驗(yàn)方法,可有效甄別一組數(shù)據(jù)是否存在異常,判別核電廠對(duì)周邊輻射環(huán)境帶來(lái)的影響,為輻射環(huán)境相關(guān)決策提供資料。

    猜你喜歡
    劑量率正態(tài)分布點(diǎn)位
    1年期國(guó)債收益率跌至1%
    銀行家(2025年1期)2025-02-08 00:00:00
    機(jī)器人快速示教方法及示教點(diǎn)位姿變換的研究
    甲狀腺乳頭狀癌患者術(shù)后首次131Ⅰ治療后輻射劑量率的影響因素及出院時(shí)間的探討
    基于對(duì)數(shù)正態(tài)分布的出行時(shí)長(zhǎng)可靠性計(jì)算
    正態(tài)分布及其應(yīng)用
    垂直面內(nèi)建立基線的特殊點(diǎn)位高程測(cè)量法
    X線照射劑量率對(duì)A549肺癌細(xì)胞周期的影響
    正態(tài)分布題型剖析
    χ2分布、t 分布、F 分布與正態(tài)分布間的關(guān)系
    ITER極向場(chǎng)線圈PF4維修區(qū)域停堆劑量率分析
    国产午夜精品论理片| 精品久久国产蜜桃| 亚洲aⅴ乱码一区二区在线播放| 亚洲人与动物交配视频| 国产精品久久电影中文字幕| 欧美日韩国产亚洲二区| 精品人妻视频免费看| 88av欧美| 嫁个100分男人电影在线观看| 亚洲美女搞黄在线观看 | 91麻豆av在线| 国产免费男女视频| 亚洲美女黄片视频| 精品久久久久久久久av| 亚洲人与动物交配视频| 国产熟女xx| 老熟妇仑乱视频hdxx| 久久99热这里只有精品18| 亚洲av第一区精品v没综合| 国产真实乱freesex| 内地一区二区视频在线| 小说图片视频综合网站| 一进一出好大好爽视频| 毛片女人毛片| 精品久久久久久久久久免费视频| 亚洲成人久久性| 欧美黑人欧美精品刺激| 成人性生交大片免费视频hd| 脱女人内裤的视频| 又粗又爽又猛毛片免费看| 少妇的逼好多水| 青草久久国产| 免费在线观看日本一区| 午夜福利高清视频| 最近视频中文字幕2019在线8| 久99久视频精品免费| 国产色婷婷99| 国产高清视频在线播放一区| 成人特级av手机在线观看| 麻豆成人午夜福利视频| 国产激情偷乱视频一区二区| 国产野战对白在线观看| 久久6这里有精品| 国产v大片淫在线免费观看| 亚洲,欧美,日韩| 国产午夜精品久久久久久一区二区三区 | 中文字幕久久专区| 在线观看av片永久免费下载| 国产一区二区亚洲精品在线观看| 在线免费观看不下载黄p国产 | 日韩免费av在线播放| 欧美黄色片欧美黄色片| 欧美日韩福利视频一区二区| 小蜜桃在线观看免费完整版高清| 国产三级在线视频| 久久久久久久精品吃奶| 级片在线观看| 欧美三级亚洲精品| 色哟哟·www| 国产精品美女特级片免费视频播放器| 99久久精品一区二区三区| 免费在线观看影片大全网站| 欧美一区二区国产精品久久精品| 很黄的视频免费| 成人精品一区二区免费| 2021天堂中文幕一二区在线观| a级毛片免费高清观看在线播放| 一本一本综合久久| 国产aⅴ精品一区二区三区波| 国产精品98久久久久久宅男小说| 精品午夜福利在线看| 欧美乱妇无乱码| 亚洲人成网站高清观看| 两个人的视频大全免费| 欧美性感艳星| 精品99又大又爽又粗少妇毛片 | 国产伦精品一区二区三区视频9| 69av精品久久久久久| 亚洲久久久久久中文字幕| 国产精品1区2区在线观看.| 丰满的人妻完整版| 别揉我奶头~嗯~啊~动态视频| 亚洲三级黄色毛片| 色在线成人网| 日本 欧美在线| 精品久久久久久久久久久久久| 精品不卡国产一区二区三区| 国产成人福利小说| 老司机午夜福利在线观看视频| 欧美日韩瑟瑟在线播放| 99热只有精品国产| 精品无人区乱码1区二区| 亚洲一区二区三区不卡视频| eeuss影院久久| 男女视频在线观看网站免费| 午夜精品一区二区三区免费看| 国产精品一区二区性色av| 变态另类丝袜制服| 欧美中文日本在线观看视频| 黄色一级大片看看| 欧美精品国产亚洲| 成年免费大片在线观看| 日本成人三级电影网站| 美女黄网站色视频| 在线看三级毛片| 亚洲最大成人中文| 精品国产三级普通话版| 嫁个100分男人电影在线观看| 怎么达到女性高潮| 神马国产精品三级电影在线观看| 丝袜美腿在线中文| 欧美最黄视频在线播放免费| 久久国产乱子免费精品| 禁无遮挡网站| 国产精品99久久久久久久久| 他把我摸到了高潮在线观看| 无人区码免费观看不卡| 在线十欧美十亚洲十日本专区| av欧美777| 午夜老司机福利剧场| 久久久精品大字幕| 精品人妻视频免费看| 亚洲精品影视一区二区三区av| 真实男女啪啪啪动态图| av在线观看视频网站免费| 亚洲av第一区精品v没综合| 日本在线视频免费播放| 我的老师免费观看完整版| 淫秽高清视频在线观看| 国产精品永久免费网站| 91午夜精品亚洲一区二区三区 | 色av中文字幕| av在线蜜桃| 国产精品爽爽va在线观看网站| 99热这里只有是精品50| 亚洲成人精品中文字幕电影| 色播亚洲综合网| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 色精品久久人妻99蜜桃| 欧美最黄视频在线播放免费| 久久国产乱子免费精品| 好男人在线观看高清免费视频| 亚洲中文日韩欧美视频| 男女做爰动态图高潮gif福利片| a在线观看视频网站| 99热精品在线国产| 人人妻人人澡欧美一区二区| 中国美女看黄片| 淫妇啪啪啪对白视频| 嫩草影视91久久| 黄色视频,在线免费观看| 91麻豆av在线| 日韩免费av在线播放| 精品免费久久久久久久清纯| 高潮久久久久久久久久久不卡| 91字幕亚洲| 亚洲七黄色美女视频| 尤物成人国产欧美一区二区三区| 日韩欧美国产一区二区入口| 亚洲精品久久国产高清桃花| 在线播放无遮挡| 亚洲美女视频黄频| 亚洲成人久久爱视频| 一个人看视频在线观看www免费| 亚洲综合色惰| 热99re8久久精品国产| 日本熟妇午夜| 变态另类丝袜制服| 极品教师在线视频| 三级国产精品欧美在线观看| 亚洲欧美日韩卡通动漫| 精品一区二区三区av网在线观看| 亚洲精品成人久久久久久| 最后的刺客免费高清国语| 小蜜桃在线观看免费完整版高清| 亚洲美女视频黄频| xxxwww97欧美| 亚洲精品在线美女| 岛国在线免费视频观看| 国产av麻豆久久久久久久| 美女免费视频网站| 久久午夜亚洲精品久久| 亚洲最大成人中文| 丰满乱子伦码专区| 午夜久久久久精精品| 欧美+亚洲+日韩+国产| 国产精品久久久久久人妻精品电影| 亚洲第一电影网av| 国产av在哪里看| 欧美xxxx黑人xx丫x性爽| 99精品在免费线老司机午夜| 国产精品av视频在线免费观看| 亚洲人成电影免费在线| 黄色日韩在线| 免费大片18禁| 窝窝影院91人妻| 3wmmmm亚洲av在线观看| 非洲黑人性xxxx精品又粗又长| 国产精品98久久久久久宅男小说| 日韩欧美免费精品| 免费一级毛片在线播放高清视频| 午夜a级毛片| 可以在线观看的亚洲视频| 欧美激情久久久久久爽电影| 淫秽高清视频在线观看| 亚洲欧美日韩高清在线视频| 最新中文字幕久久久久| 免费在线观看亚洲国产| 永久网站在线| 日本 欧美在线| 中文字幕av在线有码专区| 午夜a级毛片| 一进一出抽搐动态| 国产一级毛片七仙女欲春2| 一级毛片久久久久久久久女| 看片在线看免费视频| 国产主播在线观看一区二区| 女人被狂操c到高潮| 亚洲中文字幕一区二区三区有码在线看| av在线天堂中文字幕| 精品久久久久久久久久久久久| 欧美又色又爽又黄视频| 禁无遮挡网站| 国产精品98久久久久久宅男小说| 观看美女的网站| 亚洲在线观看片| 丁香欧美五月| aaaaa片日本免费| 国产视频一区二区在线看| 偷拍熟女少妇极品色| 久久人妻av系列| 亚洲乱码一区二区免费版| 久久久久九九精品影院| 免费无遮挡裸体视频| 欧美一区二区国产精品久久精品| 天堂网av新在线| 3wmmmm亚洲av在线观看| 色视频www国产| 亚洲国产高清在线一区二区三| 国产黄色小视频在线观看| 亚洲在线观看片| avwww免费| 美女 人体艺术 gogo| 一级黄色大片毛片| 欧美成人性av电影在线观看| 国产亚洲精品av在线| 高清日韩中文字幕在线| 欧美高清性xxxxhd video| 亚洲精品在线观看二区| 午夜精品久久久久久毛片777| 国产毛片a区久久久久| 国产精品日韩av在线免费观看| 直男gayav资源| 亚洲自拍偷在线| 好男人在线观看高清免费视频| 白带黄色成豆腐渣| 18禁在线播放成人免费| 精品久久久久久久人妻蜜臀av| 国产蜜桃级精品一区二区三区| 夜夜夜夜夜久久久久| 国内精品美女久久久久久| 女人十人毛片免费观看3o分钟| 精品久久久久久久人妻蜜臀av| 国产精品,欧美在线| 国产成人欧美在线观看| 精品人妻偷拍中文字幕| 欧美乱妇无乱码| 看十八女毛片水多多多| av视频在线观看入口| 久久久久久久久中文| 国产精品一及| 国产精品永久免费网站| 国内揄拍国产精品人妻在线| 如何舔出高潮| 久久国产乱子伦精品免费另类| 成人性生交大片免费视频hd| 日本黄大片高清| 69av精品久久久久久| 亚洲国产日韩欧美精品在线观看| 久久精品综合一区二区三区| 久久国产精品人妻蜜桃| 精品一区二区免费观看| 一卡2卡三卡四卡精品乱码亚洲| 久久精品人妻少妇| 欧美日韩福利视频一区二区| 亚洲成人久久性| 日韩高清综合在线| 网址你懂的国产日韩在线| 黄色视频,在线免费观看| 精品国内亚洲2022精品成人| 深夜精品福利| 久久人妻av系列| 日本 av在线| 久久久国产成人精品二区| 老熟妇仑乱视频hdxx| 久久久久九九精品影院| 九九在线视频观看精品| 色精品久久人妻99蜜桃| 黄色一级大片看看| 一级毛片久久久久久久久女| 中文字幕高清在线视频| 国产一区二区在线av高清观看| 九九在线视频观看精品| 中国美女看黄片| 国产精品国产高清国产av| 99国产综合亚洲精品| a级一级毛片免费在线观看| 在线免费观看的www视频| 校园春色视频在线观看| av国产免费在线观看| 国产视频内射| 免费看a级黄色片| 日韩有码中文字幕| 中文字幕av成人在线电影| 亚洲最大成人手机在线| 夜夜看夜夜爽夜夜摸| 男女视频在线观看网站免费| 人妻制服诱惑在线中文字幕| 90打野战视频偷拍视频| 成人精品一区二区免费| 亚洲无线在线观看| 小蜜桃在线观看免费完整版高清| 国产欧美日韩精品亚洲av| 国产在视频线在精品| 男女下面进入的视频免费午夜| 一级av片app| 国产男靠女视频免费网站| 国产麻豆成人av免费视频| 国内精品美女久久久久久| 午夜福利视频1000在线观看| 永久网站在线| 一本久久中文字幕| 日本黄色视频三级网站网址| 成年人黄色毛片网站| 又紧又爽又黄一区二区| 高清毛片免费观看视频网站| 亚洲中文字幕日韩| 亚洲av免费高清在线观看| 中文亚洲av片在线观看爽| 亚洲精品影视一区二区三区av| 日韩精品中文字幕看吧| 欧美黑人巨大hd| 国产毛片a区久久久久| 一夜夜www| 国产午夜精品论理片| 乱人视频在线观看| 国产熟女xx| 乱人视频在线观看| 国产免费男女视频| 欧美激情在线99| 欧美日韩瑟瑟在线播放| 一个人观看的视频www高清免费观看| 人妻制服诱惑在线中文字幕| 免费搜索国产男女视频| 久久久久国产精品人妻aⅴ院| 国产伦人伦偷精品视频| 一个人看视频在线观看www免费| 在线观看美女被高潮喷水网站 | 色哟哟·www| 啦啦啦韩国在线观看视频| 天堂动漫精品| 成人欧美大片| or卡值多少钱| 亚洲,欧美,日韩| 午夜亚洲福利在线播放| 99国产精品一区二区三区| 免费看美女性在线毛片视频| 99国产精品一区二区三区| 性色av乱码一区二区三区2| 好看av亚洲va欧美ⅴa在| 深爱激情五月婷婷| 两性午夜刺激爽爽歪歪视频在线观看| 91在线精品国自产拍蜜月| 亚洲最大成人中文| 日本 欧美在线| 欧美最新免费一区二区三区 | 欧美绝顶高潮抽搐喷水| 女人十人毛片免费观看3o分钟| 日韩精品中文字幕看吧| 一区二区三区免费毛片| 国产极品精品免费视频能看的| 免费观看的影片在线观看| 成人特级黄色片久久久久久久| 一本久久中文字幕| 成人av一区二区三区在线看| 三级毛片av免费| 精品一区二区三区人妻视频| 麻豆国产av国片精品| 亚洲精品在线美女| 成年人黄色毛片网站| 亚洲成人免费电影在线观看| 给我免费播放毛片高清在线观看| 精品一区二区三区人妻视频| 午夜激情欧美在线| 久久国产乱子免费精品| 老司机福利观看| 无遮挡黄片免费观看| 日本撒尿小便嘘嘘汇集6| 男人舔奶头视频| 精品熟女少妇八av免费久了| 欧美一区二区精品小视频在线| 国产亚洲av嫩草精品影院| 熟女电影av网| 亚洲五月婷婷丁香| 亚洲 欧美 日韩 在线 免费| 亚洲国产色片| 成人av一区二区三区在线看| bbb黄色大片| 欧美一区二区精品小视频在线| 好男人在线观看高清免费视频| 不卡一级毛片| 欧美色视频一区免费| 国产精品一区二区性色av| 男人和女人高潮做爰伦理| 国产亚洲精品久久久久久毛片| 最近最新免费中文字幕在线| 又紧又爽又黄一区二区| 亚洲av免费高清在线观看| 永久网站在线| 精品熟女少妇八av免费久了| 久久久久久九九精品二区国产| 国产人妻一区二区三区在| 国产精品综合久久久久久久免费| 精品熟女少妇八av免费久了| av中文乱码字幕在线| 99久久精品国产亚洲精品| av中文乱码字幕在线| 国产视频内射| 久久国产乱子免费精品| 亚洲avbb在线观看| 深夜a级毛片| 久久午夜亚洲精品久久| 老熟妇乱子伦视频在线观看| 国内毛片毛片毛片毛片毛片| 成人精品一区二区免费| 国产三级黄色录像| 又爽又黄a免费视频| 精品久久久久久久久亚洲 | 国产综合懂色| 久久久成人免费电影| 老司机午夜十八禁免费视频| 成人特级黄色片久久久久久久| 亚洲av第一区精品v没综合| 色哟哟哟哟哟哟| 一个人观看的视频www高清免费观看| 国产精品亚洲美女久久久| 欧美一区二区国产精品久久精品| 欧美一级a爱片免费观看看| 香蕉av资源在线| 成熟少妇高潮喷水视频| 色在线成人网| 免费在线观看成人毛片| 久久久久久九九精品二区国产| 青草久久国产| 中文亚洲av片在线观看爽| 最近最新免费中文字幕在线| 欧美不卡视频在线免费观看| 制服丝袜大香蕉在线| 色5月婷婷丁香| 欧美日韩中文字幕国产精品一区二区三区| 夜夜躁狠狠躁天天躁| 久久婷婷人人爽人人干人人爱| 久久久国产成人精品二区| 亚洲真实伦在线观看| 老司机午夜十八禁免费视频| 91字幕亚洲| 一进一出抽搐动态| 老熟妇仑乱视频hdxx| 精品一区二区三区视频在线| 三级国产精品欧美在线观看| 日韩高清综合在线| 国内精品一区二区在线观看| 超碰av人人做人人爽久久| 韩国av一区二区三区四区| 精品久久久久久久久av| 精品久久久久久久久久免费视频| 高清日韩中文字幕在线| 亚洲美女搞黄在线观看 | 午夜福利视频1000在线观看| 婷婷精品国产亚洲av| 亚洲精品影视一区二区三区av| 国产毛片a区久久久久| 不卡一级毛片| 午夜福利视频1000在线观看| 精品不卡国产一区二区三区| 亚洲精华国产精华精| 白带黄色成豆腐渣| 亚洲黑人精品在线| 国产精品乱码一区二三区的特点| 俄罗斯特黄特色一大片| 18+在线观看网站| 亚洲最大成人av| 99久久99久久久精品蜜桃| 国产三级黄色录像| av专区在线播放| 有码 亚洲区| 亚洲无线观看免费| 午夜福利在线观看吧| 老熟妇仑乱视频hdxx| 亚洲第一欧美日韩一区二区三区| 在线免费观看的www视频| 午夜日韩欧美国产| 久久天躁狠狠躁夜夜2o2o| 久久热精品热| 国产成人啪精品午夜网站| 亚洲熟妇熟女久久| 日日摸夜夜添夜夜添av毛片 | 乱码一卡2卡4卡精品| 亚洲无线在线观看| 黄色日韩在线| 非洲黑人性xxxx精品又粗又长| 国内揄拍国产精品人妻在线| 丁香六月欧美| 国产精品久久久久久人妻精品电影| 人妻夜夜爽99麻豆av| 男人狂女人下面高潮的视频| 亚洲av免费在线观看| 国产伦一二天堂av在线观看| 嫁个100分男人电影在线观看| 婷婷精品国产亚洲av在线| 亚洲男人的天堂狠狠| 亚洲精品粉嫩美女一区| 伦理电影大哥的女人| 好男人在线观看高清免费视频| 亚洲精品日韩av片在线观看| 国内精品久久久久精免费| 男女之事视频高清在线观看| 色哟哟哟哟哟哟| 中出人妻视频一区二区| 长腿黑丝高跟| 国产精品一区二区三区四区免费观看 | 婷婷六月久久综合丁香| 亚洲在线观看片| 搡老妇女老女人老熟妇| 国产极品精品免费视频能看的| 久久久国产成人精品二区| 亚洲精品色激情综合| 国产成+人综合+亚洲专区| 毛片女人毛片| 婷婷色综合大香蕉| 国产精品野战在线观看| 色精品久久人妻99蜜桃| 99国产极品粉嫩在线观看| 国产av一区在线观看免费| 午夜福利免费观看在线| 精品不卡国产一区二区三区| 两人在一起打扑克的视频| 久久中文看片网| 亚洲一区二区三区不卡视频| 亚洲最大成人手机在线| 亚洲精品成人久久久久久| 成人一区二区视频在线观看| 久久亚洲真实| 又粗又爽又猛毛片免费看| 高清在线国产一区| 国产一区二区三区在线臀色熟女| 全区人妻精品视频| 蜜桃亚洲精品一区二区三区| 99在线人妻在线中文字幕| 看十八女毛片水多多多| x7x7x7水蜜桃| 91麻豆精品激情在线观看国产| 亚洲一区二区三区不卡视频| 欧美成狂野欧美在线观看| 婷婷色综合大香蕉| 天堂av国产一区二区熟女人妻| 夜夜爽天天搞| 日本五十路高清| 久久久久久久亚洲中文字幕 | 亚洲中文字幕一区二区三区有码在线看| 免费看日本二区| 国产黄a三级三级三级人| 国产乱人伦免费视频| 亚洲欧美激情综合另类| 国产精品免费一区二区三区在线| 欧美激情国产日韩精品一区| 午夜精品久久久久久毛片777| 国产主播在线观看一区二区| 91在线观看av| 男插女下体视频免费在线播放| 人人妻,人人澡人人爽秒播| 亚洲欧美精品综合久久99| 一进一出好大好爽视频| 国产av麻豆久久久久久久| 日本熟妇午夜| 中文字幕熟女人妻在线| 伊人久久精品亚洲午夜| 精品久久久久久久末码| 1000部很黄的大片| 国产单亲对白刺激| 色综合亚洲欧美另类图片| 国产高清三级在线| 欧美性感艳星| 又紧又爽又黄一区二区| 国产真实乱freesex| 无遮挡黄片免费观看| 国产精品永久免费网站| 精品久久久久久久久av| 国产伦在线观看视频一区| 最后的刺客免费高清国语| 欧美一区二区精品小视频在线| 桃色一区二区三区在线观看| 国产伦精品一区二区三区视频9| 精品国内亚洲2022精品成人| 偷拍熟女少妇极品色| 一个人看的www免费观看视频| 1024手机看黄色片| 自拍偷自拍亚洲精品老妇| 欧洲精品卡2卡3卡4卡5卡区| 十八禁网站免费在线| 午夜视频国产福利| 国产伦人伦偷精品视频| 中文亚洲av片在线观看爽| 亚洲人与动物交配视频| 亚洲精品日韩av片在线观看|