• <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ū)域停堆劑量率分析
    99国产精品免费福利视频| 自拍欧美九色日韩亚洲蝌蚪91| 成人漫画全彩无遮挡| 99久久精品国产国产毛片| 欧美 日韩 精品 国产| av在线老鸭窝| 少妇的逼水好多| 高清在线视频一区二区三区| 在线观看一区二区三区激情| 日本欧美视频一区| 久久鲁丝午夜福利片| 少妇高潮的动态图| 国产精品99久久99久久久不卡 | 视频中文字幕在线观看| 欧美精品国产亚洲| 91久久精品国产一区二区三区| 日本av手机在线免费观看| 97在线视频观看| 91久久精品国产一区二区成人| 亚洲国产精品一区三区| 最后的刺客免费高清国语| 男人爽女人下面视频在线观看| 最近2019中文字幕mv第一页| 成人无遮挡网站| 最近2019中文字幕mv第一页| 精品亚洲成a人片在线观看| 欧美+日韩+精品| 91久久精品电影网| av免费观看日本| 丰满乱子伦码专区| 夜夜骑夜夜射夜夜干| 亚洲av不卡在线观看| 考比视频在线观看| 最近中文字幕2019免费版| 日韩av不卡免费在线播放| 亚洲精品日韩在线中文字幕| 日韩中字成人| 国产一区二区在线观看av| 日本猛色少妇xxxxx猛交久久| 有码 亚洲区| 天美传媒精品一区二区| 亚洲欧美一区二区三区国产| 新久久久久国产一级毛片| 亚洲久久久国产精品| 国产黄色免费在线视频| 久久久久久伊人网av| 91久久精品国产一区二区成人| 久久人人爽人人片av| 亚洲精品一区蜜桃| 国产成人一区二区在线| a级毛片在线看网站| 国产精品一国产av| 国产精品欧美亚洲77777| 一本大道久久a久久精品| 欧美成人午夜免费资源| 少妇被粗大的猛进出69影院 | 久久久久久久大尺度免费视频| 国产免费又黄又爽又色| 日韩精品免费视频一区二区三区 | 91成人精品电影| h视频一区二区三区| 亚洲人成77777在线视频| 这个男人来自地球电影免费观看 | 国产免费一区二区三区四区乱码| 欧美丝袜亚洲另类| 精品视频人人做人人爽| 99国产综合亚洲精品| 国产黄色免费在线视频| 亚洲精品日韩av片在线观看| 极品人妻少妇av视频| 久久99热这里只频精品6学生| 欧美日韩av久久| 久久精品国产鲁丝片午夜精品| 亚洲国产精品999| 91久久精品电影网| 日日撸夜夜添| 啦啦啦视频在线资源免费观看| av女优亚洲男人天堂| 婷婷成人精品国产| 欧美一级a爱片免费观看看| 麻豆成人av视频| 国产女主播在线喷水免费视频网站| 99久久中文字幕三级久久日本| 校园人妻丝袜中文字幕| 狂野欧美白嫩少妇大欣赏| 亚洲国产精品国产精品| 精品酒店卫生间| 国产av精品麻豆| 亚洲av成人精品一二三区| 黄色怎么调成土黄色| 国内精品宾馆在线| 丰满饥渴人妻一区二区三| 欧美精品高潮呻吟av久久| 免费看不卡的av| 中文字幕人妻丝袜制服| 99热国产这里只有精品6| 人人澡人人妻人| a级毛片黄视频| 特大巨黑吊av在线直播| 国产午夜精品一二区理论片| 国产成人freesex在线| 亚洲久久久国产精品| 日本-黄色视频高清免费观看| 少妇人妻 视频| 欧美性感艳星| 人妻人人澡人人爽人人| 欧美亚洲日本最大视频资源| 91精品一卡2卡3卡4卡| 亚洲av欧美aⅴ国产| 黑人高潮一二区| 国内精品宾馆在线| av在线app专区| 久久久国产精品麻豆| 成人黄色视频免费在线看| 丰满迷人的少妇在线观看| 亚洲精品第二区| 青春草国产在线视频| 男女无遮挡免费网站观看| 夜夜爽夜夜爽视频| 91精品三级在线观看| 欧美老熟妇乱子伦牲交| www.av在线官网国产| 午夜福利在线观看免费完整高清在| 成年美女黄网站色视频大全免费 | 亚洲国产精品一区三区| 青春草亚洲视频在线观看| 精品久久久精品久久久| 超碰97精品在线观看| 国产免费又黄又爽又色| 国产在线免费精品| 国产精品一区二区在线观看99| 夜夜看夜夜爽夜夜摸| 黑丝袜美女国产一区| 久久av网站| 亚洲av日韩在线播放| 桃花免费在线播放| 日日爽夜夜爽网站| 少妇 在线观看| 一级爰片在线观看| 亚洲成人av在线免费| 一级,二级,三级黄色视频| 最黄视频免费看| 免费观看a级毛片全部| 国产极品粉嫩免费观看在线 | 精品久久久精品久久久| 最近2019中文字幕mv第一页| 99热6这里只有精品| 黑人巨大精品欧美一区二区蜜桃 | 中国美白少妇内射xxxbb| 大话2 男鬼变身卡| 婷婷色综合大香蕉| 2018国产大陆天天弄谢| 99久国产av精品国产电影| 如何舔出高潮| 亚洲久久久国产精品| 91久久精品电影网| 国产精品秋霞免费鲁丝片| 免费观看性生交大片5| 少妇人妻精品综合一区二区| 亚洲不卡免费看| 国产毛片在线视频| 纵有疾风起免费观看全集完整版| 我的女老师完整版在线观看| 久久久久久久精品精品| 亚洲五月色婷婷综合| 乱码一卡2卡4卡精品| 黄色配什么色好看| 亚洲精品乱码久久久久久按摩| 中文精品一卡2卡3卡4更新| 91精品伊人久久大香线蕉| 在线观看www视频免费| 欧美变态另类bdsm刘玥| 国产av精品麻豆| 国产极品粉嫩免费观看在线 | 三上悠亚av全集在线观看| 亚洲精品一区蜜桃| 久久鲁丝午夜福利片| 欧美日韩在线观看h| 精品一区在线观看国产| 女性被躁到高潮视频| 涩涩av久久男人的天堂| av有码第一页| 色婷婷av一区二区三区视频| 99久国产av精品国产电影| 精品一区二区三卡| 婷婷色麻豆天堂久久| 久久午夜综合久久蜜桃| 亚洲欧美一区二区三区国产| 在线观看一区二区三区激情| 国产男女超爽视频在线观看| 国产一区二区在线观看日韩| av国产精品久久久久影院| 麻豆乱淫一区二区| 免费少妇av软件| 欧美日韩一区二区视频在线观看视频在线| 日本wwww免费看| 亚洲四区av| 日韩视频在线欧美| 婷婷色综合大香蕉| 亚洲人与动物交配视频| 五月开心婷婷网| 午夜久久久在线观看| 美女视频免费永久观看网站| 天天躁夜夜躁狠狠久久av| 亚洲欧美日韩卡通动漫| 国产精品国产av在线观看| 国产欧美日韩综合在线一区二区| 少妇熟女欧美另类| 国产免费一级a男人的天堂| 午夜免费鲁丝| 国产精品国产av在线观看| 国产精品秋霞免费鲁丝片| 狂野欧美激情性bbbbbb| 精品熟女少妇av免费看| 最近中文字幕高清免费大全6| 免费少妇av软件| 观看av在线不卡| 老司机影院成人| 亚洲国产日韩一区二区| a级毛片免费高清观看在线播放| 久久久久久人妻| 日韩人妻高清精品专区| 一区二区三区乱码不卡18| 国产精品不卡视频一区二区| 久久久久久伊人网av| 国产乱来视频区| 99视频精品全部免费 在线| 久久久久人妻精品一区果冻| 免费人成在线观看视频色| 欧美bdsm另类| 天堂俺去俺来也www色官网| 99视频精品全部免费 在线| 黄片播放在线免费| 另类亚洲欧美激情| 各种免费的搞黄视频| 超碰97精品在线观看| 男女免费视频国产| 欧美日韩国产mv在线观看视频| 日日摸夜夜添夜夜添av毛片| 七月丁香在线播放| 久久久久久久国产电影| 黄片无遮挡物在线观看| 国产精品99久久99久久久不卡 | 人妻少妇偷人精品九色| 国产精品久久久久久久电影| 国产亚洲午夜精品一区二区久久| 水蜜桃什么品种好| 成年女人在线观看亚洲视频| 亚洲国产精品一区二区三区在线| 日韩大片免费观看网站| 日韩亚洲欧美综合| 伊人久久国产一区二区| 黑丝袜美女国产一区| 久久久久久伊人网av| 亚洲av免费高清在线观看| 久久国产精品男人的天堂亚洲 | 在线观看免费视频网站a站| 日韩电影二区| 中国美白少妇内射xxxbb| 国产在线免费精品| 国产伦精品一区二区三区视频9| 最后的刺客免费高清国语| 91成人精品电影| 亚洲综合精品二区| 国产日韩欧美亚洲二区| 岛国毛片在线播放| 亚洲激情五月婷婷啪啪| 十分钟在线观看高清视频www| 少妇猛男粗大的猛烈进出视频| 丝袜在线中文字幕| 国产av一区二区精品久久| 久久久久久久久久久久大奶| 欧美日韩在线观看h| 中国国产av一级| 国产精品一区www在线观看| 国产精品国产三级国产专区5o| 国产不卡av网站在线观看| 亚洲欧美成人综合另类久久久| 国产精品久久久久成人av| 丝袜喷水一区| 久久免费观看电影| 国产男女内射视频| 国产精品国产三级国产专区5o| 午夜视频国产福利| 母亲3免费完整高清在线观看 | 人妻少妇偷人精品九色| 男女无遮挡免费网站观看| 寂寞人妻少妇视频99o| 日韩制服骚丝袜av| 一本久久精品| 丝袜喷水一区| 肉色欧美久久久久久久蜜桃| 一本色道久久久久久精品综合| 丰满少妇做爰视频| 97在线人人人人妻| 纵有疾风起免费观看全集完整版| 最后的刺客免费高清国语| 免费大片黄手机在线观看| 日韩强制内射视频| 国产成人午夜福利电影在线观看| 欧美另类一区| 亚洲人成网站在线观看播放| 国产男女内射视频| 美女中出高潮动态图| 女人精品久久久久毛片| 18禁观看日本| 哪个播放器可以免费观看大片| 在线免费观看不下载黄p国产| 狠狠婷婷综合久久久久久88av| 91久久精品国产一区二区成人| 高清视频免费观看一区二区| 美女xxoo啪啪120秒动态图| 亚洲av成人精品一区久久| 九九久久精品国产亚洲av麻豆| 欧美 日韩 精品 国产| 只有这里有精品99| 久久久精品区二区三区| 亚洲精品,欧美精品| 欧美xxxx性猛交bbbb| 男女免费视频国产| 男人爽女人下面视频在线观看| 欧美日韩一区二区视频在线观看视频在线| 一本久久精品| 欧美三级亚洲精品| 一本一本综合久久| 成年美女黄网站色视频大全免费 | a级毛片在线看网站| 热99国产精品久久久久久7| 国产高清不卡午夜福利| 最后的刺客免费高清国语| 亚洲经典国产精华液单| 亚洲成人手机| 99精国产麻豆久久婷婷| 精品久久久精品久久久| 肉色欧美久久久久久久蜜桃| 2021少妇久久久久久久久久久| 免费高清在线观看日韩| 能在线免费看毛片的网站| 国产在线免费精品| 热re99久久国产66热| 亚洲中文av在线| 免费观看性生交大片5| 国产精品99久久久久久久久| 人人妻人人添人人爽欧美一区卜| 伊人久久精品亚洲午夜| 日日摸夜夜添夜夜爱| 日韩,欧美,国产一区二区三区| 综合色丁香网| 欧美激情极品国产一区二区三区 | 欧美 日韩 精品 国产| 欧美日韩视频精品一区| 亚洲成人一二三区av| 国产精品成人在线| 飞空精品影院首页| 少妇人妻 视频| 熟妇人妻不卡中文字幕| 国产一区二区在线观看av| 九九爱精品视频在线观看| 最新中文字幕久久久久| 国产深夜福利视频在线观看| 色婷婷av一区二区三区视频| 在线观看人妻少妇| 日日撸夜夜添| 在线观看一区二区三区激情| 2018国产大陆天天弄谢| 国产精品国产av在线观看| 黑人猛操日本美女一级片| 日韩精品免费视频一区二区三区 | 国产黄片视频在线免费观看| 人妻制服诱惑在线中文字幕| 午夜福利视频精品| 精品国产露脸久久av麻豆| av电影中文网址| 插逼视频在线观看| 丰满少妇做爰视频| 亚洲国产精品一区二区三区在线| 黄色视频在线播放观看不卡| 汤姆久久久久久久影院中文字幕| 婷婷色av中文字幕| 国产国语露脸激情在线看| 亚洲不卡免费看| 亚洲国产av新网站| 日本91视频免费播放| 国产成人免费观看mmmm| 国产日韩欧美在线精品| 国产精品欧美亚洲77777| 考比视频在线观看| 在线观看美女被高潮喷水网站| 一个人看视频在线观看www免费| 日本爱情动作片www.在线观看| 亚洲精品美女久久av网站| 国精品久久久久久国模美| 性高湖久久久久久久久免费观看| 精品久久久精品久久久| 成人毛片60女人毛片免费| 亚洲国产色片| 欧美日韩成人在线一区二区| 热re99久久精品国产66热6| 亚洲精品日韩av片在线观看| 国产深夜福利视频在线观看| 久久久久视频综合| 一区二区日韩欧美中文字幕 | 午夜福利视频精品| 国产有黄有色有爽视频| a级毛片黄视频| 日韩av不卡免费在线播放| 一级黄片播放器| 中文欧美无线码| 如何舔出高潮| 色婷婷av一区二区三区视频| 国产色婷婷99| 亚洲欧洲精品一区二区精品久久久 | 又黄又爽又刺激的免费视频.| 一本—道久久a久久精品蜜桃钙片| 亚洲第一av免费看| 久久久精品94久久精品| 麻豆成人av视频| 国产av国产精品国产| 亚洲av二区三区四区| 亚洲精品国产av成人精品| 高清黄色对白视频在线免费看| 水蜜桃什么品种好| 欧美日韩在线观看h| 在线亚洲精品国产二区图片欧美 | 亚洲丝袜综合中文字幕| 9色porny在线观看| 大香蕉久久网| 久久国内精品自在自线图片| 日韩av在线免费看完整版不卡| 国产精品一二三区在线看| 国产 一区精品| 自线自在国产av| 桃花免费在线播放| 人妻 亚洲 视频| 高清欧美精品videossex| 熟女电影av网| 曰老女人黄片| 天天躁夜夜躁狠狠久久av| 国产亚洲一区二区精品| 免费播放大片免费观看视频在线观看| tube8黄色片| 2018国产大陆天天弄谢| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 亚洲av国产av综合av卡| 中文字幕亚洲精品专区| 国产精品国产三级国产av玫瑰| 精品亚洲乱码少妇综合久久| 午夜激情久久久久久久| 亚洲久久久国产精品| 91在线精品国自产拍蜜月| .国产精品久久| 国产精品嫩草影院av在线观看| 成人毛片a级毛片在线播放| 久久精品夜色国产| 精品国产乱码久久久久久小说| 涩涩av久久男人的天堂| 黑人猛操日本美女一级片| 亚洲,欧美,日韩| 美女内射精品一级片tv| 午夜福利网站1000一区二区三区| 久久精品久久久久久久性| 视频区图区小说| 在线观看一区二区三区激情| 久久午夜综合久久蜜桃| 久久精品人人爽人人爽视色| 老司机影院成人| 汤姆久久久久久久影院中文字幕| 蜜桃久久精品国产亚洲av| av又黄又爽大尺度在线免费看| 少妇的逼水好多| 一区二区三区四区激情视频| 狠狠婷婷综合久久久久久88av| 狂野欧美激情性bbbbbb| 欧美bdsm另类| 精品一区二区三卡| 人人妻人人添人人爽欧美一区卜| 99视频精品全部免费 在线| 精品人妻熟女av久视频| 日韩中字成人| 久久久久久久亚洲中文字幕| av天堂久久9| 伦理电影大哥的女人| 久久国内精品自在自线图片| 只有这里有精品99| 十八禁高潮呻吟视频| 一边亲一边摸免费视频| 伦精品一区二区三区| 美女内射精品一级片tv| 久久免费观看电影| 中文字幕免费在线视频6| 97超视频在线观看视频| 91精品一卡2卡3卡4卡| av卡一久久| 满18在线观看网站| 亚洲av综合色区一区| 久久久久久久久大av| 日本欧美视频一区| 日本猛色少妇xxxxx猛交久久| 午夜免费男女啪啪视频观看| av天堂久久9| 人妻一区二区av| 日韩成人av中文字幕在线观看| 亚洲少妇的诱惑av| 免费看不卡的av| 少妇的逼好多水| 精品人妻在线不人妻| 99九九线精品视频在线观看视频| 国产在线视频一区二区| 王馨瑶露胸无遮挡在线观看| 天堂中文最新版在线下载| 亚洲欧美一区二区三区国产| 超色免费av| 特大巨黑吊av在线直播| 亚洲综合色惰| 国产日韩欧美在线精品| 亚洲精品国产色婷婷电影| 蜜桃在线观看..| 天美传媒精品一区二区| 久久久久久久久久久免费av| 美女内射精品一级片tv| 久久久久久久久久人人人人人人| 精品少妇内射三级| 尾随美女入室| 高清不卡的av网站| 大陆偷拍与自拍| 亚洲精品日本国产第一区| 亚洲精品久久久久久婷婷小说| 晚上一个人看的免费电影| 菩萨蛮人人尽说江南好唐韦庄| 日本欧美视频一区| a级毛片免费高清观看在线播放| 亚洲第一av免费看| 激情五月婷婷亚洲| 在线精品无人区一区二区三| 午夜福利视频精品| 肉色欧美久久久久久久蜜桃| a级毛片黄视频| 最新的欧美精品一区二区| 丰满饥渴人妻一区二区三| 街头女战士在线观看网站| 超碰97精品在线观看| 91国产中文字幕| 女的被弄到高潮叫床怎么办| 欧美日韩视频精品一区| 免费黄网站久久成人精品| 国产成人免费观看mmmm| 国产免费又黄又爽又色| 男女高潮啪啪啪动态图| 日韩成人伦理影院| 久久97久久精品| 欧美日韩一区二区视频在线观看视频在线| 欧美日韩av久久| 只有这里有精品99| 18禁观看日本| 满18在线观看网站| 麻豆精品久久久久久蜜桃| 日本色播在线视频| 亚洲不卡免费看| 久久人人爽人人爽人人片va| 国产不卡av网站在线观看| 亚洲av电影在线观看一区二区三区| 人妻一区二区av| 日韩强制内射视频| 精品人妻熟女毛片av久久网站| 黄色毛片三级朝国网站| 在线天堂最新版资源| 久久久久国产精品人妻一区二区| 欧美日韩在线观看h| 丝瓜视频免费看黄片| 日韩中文字幕视频在线看片| 80岁老熟妇乱子伦牲交| 午夜免费鲁丝| 精品酒店卫生间| 日韩免费高清中文字幕av| 天天操日日干夜夜撸| 下体分泌物呈黄色| 午夜日本视频在线| 亚洲成人一二三区av| 搡女人真爽免费视频火全软件| 国产亚洲精品第一综合不卡 | 少妇被粗大猛烈的视频| 高清毛片免费看| 国产一区亚洲一区在线观看| 久热这里只有精品99| 在线观看国产h片| 十八禁网站网址无遮挡| 一级毛片我不卡| 男女免费视频国产| 99热全是精品| 久久精品国产自在天天线| av在线老鸭窝| 日韩强制内射视频| 国产成人精品福利久久| 99久久精品国产国产毛片| 视频在线观看一区二区三区| 如何舔出高潮| 成人国语在线视频| 两个人免费观看高清视频| 久久久久久久精品精品| 国产有黄有色有爽视频| 97精品久久久久久久久久精品| 日本午夜av视频| 老女人水多毛片| 国产高清国产精品国产三级| 免费观看在线日韩| 亚洲精品国产av蜜桃| 夫妻午夜视频| 免费观看av网站的网址| 蜜桃久久精品国产亚洲av| 18禁在线播放成人免费| 91久久精品电影网| 涩涩av久久男人的天堂| 国精品久久久久久国模美| 国产极品天堂在线| 日韩视频在线欧美|