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

    輻射環(huán)境監(jiān)測數(shù)據(jù)合理性評價中的問題探討

    2017-06-22 14:42:22沙連茂
    關(guān)鍵詞:比活度活度放射性

    沙連茂

    (中國輻射防護(hù)研究院, 山西 太原 030006)

    ·特約來稿·

    輻射環(huán)境監(jiān)測數(shù)據(jù)合理性評價中的問題探討

    沙連茂

    (中國輻射防護(hù)研究院, 山西 太原 030006)

    對輻射環(huán)境監(jiān)測數(shù)據(jù)合理性評價中探測限附近的數(shù)據(jù)、監(jiān)測結(jié)果趨勢的分析和實驗室分析過程控制中的相關(guān)問題進(jìn)行分析和探討,提出了意見和建議。

    環(huán)境監(jiān)測;輻射;合理性評價;測量數(shù)據(jù)

    在輻射環(huán)境監(jiān)測中按照方案設(shè)計可以獲得一定數(shù)量的樣本數(shù)據(jù),并據(jù)此對總體做出推斷。判斷環(huán)境監(jiān)測數(shù)據(jù)及其推斷結(jié)論合理性的依據(jù)很多,主要是看其在統(tǒng)計學(xué)上和專業(yè)知識上是否合理。不同的領(lǐng)域,衡量的側(cè)重點會有所不同。對計量工作來說,準(zhǔn)確性、一致性、溯源性和法制性是其基本特性?,F(xiàn)就有關(guān)輻射環(huán)境監(jiān)測數(shù)據(jù)合理性評價中的有關(guān)問題進(jìn)行探討。

    1 探測限附近的數(shù)據(jù)

    在環(huán)境放射性監(jiān)測過程中,會獲得大量探測限附近、甚至小于探測限的數(shù)據(jù)。美國能源部設(shè)施環(huán)境放射性監(jiān)督導(dǎo)則[1]將這類環(huán)境監(jiān)測數(shù)據(jù)描述為:“不對稱和混合分布,連同極端值和低于探測限值的存在是環(huán)境放射性監(jiān)督數(shù)據(jù)的通常特征。具有這種特征的監(jiān)督數(shù)據(jù)的分析和隨后的報告是不容易的,或者說是不能按常規(guī)處理的。”

    1.1 探測限附近測量結(jié)果的不確定度

    圖1示出探測限附近測量結(jié)果在置信度95%和99%時的不確定度曲線, 圖中LQ為定量測定限,Ld為探測限。

    圖1 置信度95%和99%時探測限附近測量結(jié)果的不確定度

    由圖1可見,Ld附近測量結(jié)果的不確定度可以達(dá)到±100%,所以只要數(shù)據(jù)的偏差是不帶方向性的,測量結(jié)果在較大的離散度下漲落也不足為奇。有的質(zhì)保計劃規(guī)定輻射環(huán)境監(jiān)測數(shù)據(jù)的誤差要控制在5%~10%以內(nèi),在一般情況下難以實現(xiàn)。在同樣條件下獲得的數(shù)據(jù)組中未發(fā)現(xiàn)異常值,而發(fā)現(xiàn)有的數(shù)據(jù)組精度太好的情況,所得χ2(卡方)值明顯偏低,反映出有的監(jiān)測人員可能對低水平測量結(jié)果在較大的離散度下漲落不太理解,誤以為這是不正?,F(xiàn)象,有的甚至對測量數(shù)據(jù)進(jìn)行“修正”。 例如在輻射環(huán)境監(jiān)測比對活動中,曾不時出現(xiàn)“有選擇性的記錄變化不大的數(shù)據(jù)值”的問題,這是缺乏客觀、隨機(jī)地記錄測量數(shù)據(jù)的基本素質(zhì)和實事求是科學(xué)態(tài)度的表現(xiàn)。

    1.2 采用判斷限作為測量結(jié)果的評價標(biāo)準(zhǔn)

    探測限Ld是反映一種特定測量(包括方法、儀器、樣品特征等)的技術(shù)指標(biāo),用于評述一種測量的檢測能力,是評價輻射環(huán)境監(jiān)測方法合理性的一個重要參數(shù),而不是用于表述測量結(jié)果。判斷限Lc通常用于檢驗測量結(jié)果是否在統(tǒng)計上與本底有顯著性差異,又稱判斷閾(decision threshold)。探測限Ld可以判“有”,但不能判”無”。如果有一個活度為Ld的樣品,由于測量的統(tǒng)計誤差,測量結(jié)果>Ld的概率只有約50%,

    根據(jù)Ld的定義,若取β=0.05,可理解為如果樣品的放射性期望值達(dá)到了Ld,那么可以(事先)斷言若采用該測量裝置,樣品的計數(shù)結(jié)果會有95%的概率≥Lc(注意:不是有95%的概率≥Ld),而

    實際操作中可保守地認(rèn)為判斷限等于探測限的一半。鑒于Lc在數(shù)值上比Ld小,因此,采用Lc作為評判測量結(jié)果的標(biāo)準(zhǔn)比單純用

    還要明確,測量一個樣品的探測限、判斷限和測量若干個樣品平均值的探測限、判斷限是有差別的。例如,對同一樣品重復(fù)測量100次,其探測限、判斷限便降低為原來的1/10。

    1.3 小于判斷限的數(shù)據(jù)的報告

    如何報告

    美國NBS提出報小于水平Lt[4]:

    Lt=R+Kσs

    (1)

    式中:R——樣品的凈計數(shù)率;K——單側(cè)置信因子,在95%的置信度下K=1.645;σs——樣品測量的標(biāo)準(zhǔn)差。

    此法具有嚴(yán)格的統(tǒng)計學(xué)依據(jù),報出的水平與測量結(jié)果“掛鉤”,在20世紀(jì)80年代我國進(jìn)行的食品放射性水平調(diào)查處理數(shù)據(jù)時,曾采用了這種方法。

    近年來歐盟對小于判斷限的核電廠排放流出物的監(jiān)測結(jié)果提出新的建議書[2004/2/Euratom],提出測量結(jié)果低于判斷閾的,結(jié)果可以保守地以判斷閾的一半來替代。如不斷重復(fù)測量的結(jié)果都低于判斷閾時,可合理地認(rèn)為真值就是零,即排放水中不存在此種放射性核素。我國在核電廠流出物監(jiān)測中正在參考此建議書[5]。

    1.4 提供樣本容量和標(biāo)準(zhǔn)差(或不確定度)信息

    輻射環(huán)境監(jiān)測數(shù)據(jù)相當(dāng)一部分都是處于探測限附近,為了反映測量的精確度,報告測量結(jié)果時,不僅要提供結(jié)果的平均值,還要提供樣本容量和標(biāo)準(zhǔn)差(或不確定度)。對于多個樣品重復(fù)測量的結(jié)果,還要提供超出判斷限的樣品數(shù)占總樣品數(shù)比率的信息。通過樣本估計總體時,樣本容量n的大小反映了估計的可靠程度。

    (1)樣本容量大小不僅關(guān)系到結(jié)果的偏差能否接受,且影響統(tǒng)計檢驗的實際意義。在統(tǒng)計檢驗中,n大小也直接影響檢驗結(jié)果的可靠性;

    (2)在等精度的重復(fù)測量中樣本容量的大小反映樣本均值的精度。

    標(biāo)準(zhǔn)差,也是測量結(jié)果的不確定度,又稱標(biāo)準(zhǔn)不確定度。測量的不確定度是表征合理地賦予被測量值的分散性、與測量結(jié)果相聯(lián)系的參數(shù)。從詞義上理解,意味著對測量結(jié)果可信性、有效性的懷疑程度或不肯定程度,是定量說明測量結(jié)果質(zhì)量的一個參數(shù)。放射性測量的泊松均方差和高斯標(biāo)準(zhǔn)差的計算較為簡單。關(guān)于化學(xué)分析結(jié)果不確定度的評定,中國實驗室認(rèn)可委員會已于2001年頒布了不確定度的評估指南[6],應(yīng)按照指南的要求去進(jìn)行不確定度的評定。對照“指南”的要求,一些文章中存在的問題有:數(shù)學(xué)模式與不確定度的來源相互之間缺乏聯(lián)系,存在“兩張皮”的現(xiàn)象,即有的分量脫離數(shù)學(xué)模式中的參數(shù)隨意評定結(jié)果的不確定度。應(yīng)當(dāng)緊緊圍繞數(shù)學(xué)模式的參數(shù),逐個分析、評定,最后按模式中各參數(shù)的函數(shù)關(guān)系進(jìn)行不確定度的合并。

    圖2 平均值的標(biāo)準(zhǔn)差隨樣本容量的變化

    1.5 探測限附近測量數(shù)據(jù)的有效性

    測量數(shù)據(jù)的有效性可以用某測量系統(tǒng)在當(dāng)時的測量條件下預(yù)期能測到最小凈計數(shù)期望值μN(yùn)=Ld的把握度(1-β)來評價,(1-β)被稱之為檢驗規(guī)則的功效。β在統(tǒng)計學(xué)上稱為犯第二類錯誤的概率,即樣品實際上含有放射性(μN(yùn)>0),卻被判為不含放射性(μN(yùn)=0)的概率。實際計算時可按下式獲得Uβ值。

    Uβ=(ns-Lc)/σs

    (2)

    式中:Uβ——標(biāo)準(zhǔn)正態(tài)變量;ns——樣品的凈計數(shù)率;Lc——判斷限;σs——樣品計數(shù)測量的標(biāo)準(zhǔn)差。

    按式(2)計算出Uβ值后,便可以在數(shù)理統(tǒng)計表中查得β值,從而做出把握度1-β的估計。表1列出采用3種不同測量方案對一個樣品進(jìn)行計數(shù)時獲得的Lc、Uβ和(1-β)值。從中可看出,即使其獲得的凈計數(shù)率相同,均為0.15 cpm,但由于本底或測量時間不同,所得結(jié)果的把握度是有差別的,有的甚至差別很大。

    表1 凈計數(shù)率相同、測量方案不同時的Lc、Uβ和(1-β)值

    由于必然性因素(核設(shè)施的排放)的影響總是和偶然性(包括低水平放射性環(huán)境樣品的采樣、制樣、分析和測量過程中存在的誤差)因素的影響糾纏在一起,只有排除偶然性因素的影響,才能做出是否受核設(shè)施的排放這一必然性因素影響的判斷。

    因此需要強(qiáng)調(diào),從事核設(shè)施本底調(diào)查工作所采用的方法和技術(shù)應(yīng)當(dāng)具有比常規(guī)工作更好的靈敏度和可靠性,獲得更為有效的測量數(shù)據(jù),而不要僅滿足于 “測量結(jié)果

    2 監(jiān)測結(jié)果趨勢的分析

    首先考慮,一組數(shù)據(jù)的范圍值是否在合理范圍內(nèi)。在沒有外來源或其他特殊因素存在時,將范圍值與文獻(xiàn)報道的環(huán)境和地理條件類似的同類數(shù)據(jù)范圍比較,如發(fā)現(xiàn)超出合理范圍,便有可能在監(jiān)測中存在缺陷。這類檢查,還應(yīng)包括其他不合理數(shù)據(jù)的檢查,如數(shù)據(jù)的有效數(shù)字的位數(shù)偏差太大,平均值超出范圍值;單個核素的活度濃度超過總放射性的活度濃度;在核電站運行期間,人工核素在遠(yuǎn)區(qū)明顯大于近區(qū);鈾/鐳比例不合理;全球沉降的人工核素和某些大氣沉降的天然核素在樣品中的活度濃度比不合理等。總之,還要根據(jù)專業(yè)知識,更有效地利用數(shù)理統(tǒng)計方法來發(fā)現(xiàn)問題。

    2.1 判斷監(jiān)測結(jié)果之間是否存在顯著差異的問題

    在測量條件相對穩(wěn)定的情況下,同一點不同時間的樣品數(shù)據(jù)應(yīng)當(dāng)滿足正態(tài)或?qū)?shù)正態(tài)分布,對歷年周期為月或季的測量數(shù)據(jù),可采用質(zhì)控圖方法檢驗。此外,也可檢查各年間的成對比較是否有顯著性差異,若發(fā)現(xiàn)某些年度間或與本底調(diào)查結(jié)果比較存在顯著性差異,則說明這些年度間或存在不同污染源的影響,或監(jiān)測中存在缺陷。應(yīng)當(dāng)注意比較不同樣本結(jié)果之間的差異是否顯著,不能僅看其平均值,還要看其不確定度。有的文章在比較不同江河、或不同采樣點、或同一采樣點不同時間的環(huán)境介質(zhì)放射性核素的活度濃度時,僅根據(jù)其平均值就排列出其活度濃度的高低順序,事實上如果考慮不確定度,其置信區(qū)間是深度重復(fù)交叉的,看不出它們之間存在顯著性差異。檢查2組數(shù)據(jù)平均值的置信區(qū)間是否相重疊,在樣品數(shù)較少時很有用。2組數(shù)據(jù)若能看作來自同一總體,則其均值的置信區(qū)間不能是分離的。

    某篇文章對2007—2015年某主要江河的枯水期和豐水期放射性核素水平進(jìn)行比較,所得數(shù)據(jù)見表2。作者僅僅依據(jù)平均值的差別就認(rèn)為226Ra、40K、90Sr、總α枯水期和豐水期之間差別非常顯著。但是,根據(jù)表2提供的枯水期和豐水期的平均值和標(biāo)準(zhǔn)差,其置信區(qū)間都是高度重合的,看不出“差別非常顯著”。

    表2 枯水期和豐水期放射性核素水平比較

    另篇關(guān)于土壤中232Th含量的文章,說某地區(qū)的測量值大于全國調(diào)研值,但該地區(qū)3個點平均測量值為49.7 Bq/kg, 單次測量標(biāo)準(zhǔn)差11.6 Bq/kg。而全國土壤中232Th的含量為(49.1±27.6) Bq/kg。顯然,該地區(qū)的測量值與全國調(diào)研值并無顯著性差異。

    當(dāng)各年間測量條件及環(huán)境條件相同時,各年間的成對比較應(yīng)該不會有顯著性差異,一旦發(fā)現(xiàn)存在顯著性差異,說明這些年度間或有新的污染源影響,或監(jiān)測中存在缺陷。在進(jìn)行統(tǒng)計檢驗時,成對比較法比成組比較法要靈敏得多;總體均值未見差異時,還要觀測各點各年間是否有顯著性差異。當(dāng)發(fā)現(xiàn)各年間存在顯著性差異時,應(yīng)當(dāng)查明原因,不要輕易放過。曾有一篇水體調(diào)查的文章,發(fā)現(xiàn)多年來某河水中238U、232Th、226Ra、40K、137Cs活度濃度均有上升趨勢, 尤其是鈾,與全省江河相比,該河段不僅高于全省江河,比1987年此點位參考范圍水平高出近9倍。但作者只滿足于其遠(yuǎn)低于我國規(guī)定的露天水源的限制濃度,并未查明原因,給公眾一個交代。

    還有一篇關(guān)于地表水中238U測量的文章,檢測數(shù)值很不穩(wěn)定,同一斷面不同年份的數(shù)值最大相差41倍。作者認(rèn)為其原因可能是采集到較多的懸浮物,這有可能,但作者又以“檢測數(shù)據(jù)最大值沒有超過本底水平的3倍以上”為由,試圖大事化小。從238U的監(jiān)測結(jié)果看,有幾個特別低的結(jié)果基本上都出現(xiàn)在2011年。例如,某壩2010年為46.35 mBq/L,而2011年上半年突變?yōu)?.7mBq/L; 某大橋2010年為79 mBq/L, 2011年突變?yōu)?.9 mBq/L; 還有一橋2010年為67 mBq/L, 2011年突變?yōu)?.6 mBq/L。這里是否存在監(jiān)測技術(shù)方面的原因,也需要作者認(rèn)真自查。

    2.2 全球沉降的人工核素和某些大氣沉降的天然核素在樣品中的活度濃度比是否合理

    一些大氣沉降的天然放射性核素,如210Pb和7Be等,在一定條件下和全球沉降的人工核素在雨水、沉降灰、土壤和茶葉等樣品中的活度濃度往往呈一定的比例關(guān)系。

    210Pb是天然鈾-鐳放射系中的氡子體,氡從地殼逸出到大氣中衰變生成一系列子體,其中210Pb是壽命較長的核素,與全球沉降的人工核素90Sr、137Cs和239Pu等一起沉降到地表。1990年曾分析過北京、太原、石家莊和濟(jì)南表層土壤中137Cs、239Pu、241Am 與210Pb的沉積量[7],圖3中僅顯示出137Cs 與210Pb 在土壤表面的沉積量關(guān)系圖。

    圖3 137Cs 與210Pb 在土壤表面的沉積量

    分析結(jié)果表明,地表土壤中210Pb、137Cs、239Pu、241Am沉積量的比例大約為4∶1∶0.02∶0.007。了解樣品中核素間的大致比值,除了用來評價測量結(jié)果的可靠性外,也有助于采樣計劃的制定。例如,可以通過γ譜法對土壤樣品用非破壞性測量方法獲得210Pb和137Cs的含量,按照上述比例關(guān)系估計土壤中241Am和239Pu的大致含量,以便采取足夠量的樣品用于超鈾核素的分析,才能獲得有效的數(shù)據(jù)。

    采樣的代表性往往很難被定量,因為不知道“真值”,但在實際測量數(shù)據(jù)的支持下,應(yīng)用科學(xué)的方法可以對數(shù)據(jù)的合理性進(jìn)行估計,做出定性說明。某些情況下,像7Be這樣的天然放射性核素也能助一臂之力。

    7Be 為大氣層中宇宙射線與氮氧相互作用的產(chǎn)物,半衰期53.3 d, 產(chǎn)率: 0.081 原子/cm2· s。7Be在對流層空氣中濃度為1.04×10-2Bq/m3,雨水中平均濃度約0.7×10-2Bq/L。它能附著在落下灰和雨水上而成為樣品時效性的有用指示劑。一般情況下,在同一時期現(xiàn)場采集的茶葉中7Be 的含量應(yīng)當(dāng)如人工核素137Cs那樣有相關(guān)關(guān)系。在田灣核電站申請裝料許可證階段環(huán)境輻射本底調(diào)查期間,測得2001年6月8日采集的云臺山和花果山的云霧茶中137Cs的含量均為(0.23±0.01) Bq/kg,但是,7Be的含量差別很大,云臺山云霧茶中7Be的含量為(44±10)Bq/kg,而花果山的云霧茶中7Be卻異常得低,含量<0.3 Bq/kg。質(zhì)量保證小組對其時效性提出質(zhì)疑。后經(jīng)調(diào)查確認(rèn),該批花果山的云霧茶不是采樣人員親臨茶場摘采,而是為了“省事”,從茶葉店中購買的一年前的茶葉包裝商品。由于137Cs的半衰期較長,活度濃度幾乎不變;而7Be的半衰期比137Cs短得多,只有53.3 d,經(jīng)過1年多的時間,已基本衰變完了。這樣,用γ譜儀未能測出7Be,而137Cs的含量基本正常。這表明該茶葉樣品不能代表此時期的茶葉產(chǎn)品。

    2.3 樣品中個別核素的活度與總放射性水平關(guān)系

    一般而言,個別核素的活度濃度是不會超出總活度濃度水平,否則,其中2者之一必定存在問題。有一篇關(guān)于稀土廠排放廢水中總α和天然放射性核素的報告,2008年處理后水樣的總α放射性濃度為0.026 Bq/L,還沒有單個核素226Ra的濃度水平(0.056 Bq/L)高;同樣,某廠2009年處理后水樣的總α放射性濃度為0.027 Bq/L, 還沒有單個核素238U的水平(0.067 Bq/L)高;2009年削減斷面水樣的總α放射性濃度為0.042 Bq/L, 比單個核素238U水平(0.063 Bq/L)還低。

    但是,在考慮含鉀量高的樣品中的總β水平與40K的活度濃度時會出現(xiàn)例外情況。在不少生物樣品中,能放射出β輻射的核素有鈾釷的衰變子體、40K及微量的90Sr和137Cs等,其中40K的含量最高,其活度濃度比其他核素大約高出2個數(shù)量級以上。從40K衰變綱圖可知,每衰變一次,有約89%的概率放射出β輻射。通常40K的活度濃度是通過γ譜儀、原子吸收分光光度計或化學(xué)方法測定,而其總β放射性的活度是通過β計數(shù)器測定的,因此,40K的活度濃度與總β輻射水平的比值為1.12,即總β放射性的活度僅為40K的活度濃度的89%。例如,一些南方水果、蔬菜的總β與40K比值為0.75~1.05,平均值為0.88[8]。

    2.4 變量之間相關(guān)性的判斷

    在環(huán)境放射性監(jiān)測工作中,經(jīng)常會遇到相互間存在一定聯(lián)系的變量?;貧w分析和相關(guān)分析是處理變量間相關(guān)關(guān)系的有力工具。其中的相關(guān)分析是研究變量之間的相關(guān)性及其密切程度。但在通過測量獲得各種變量之后,要判斷其是否相關(guān),相關(guān)密切程度如何,并不是僅僅用肉眼的觀察所能解決的,而是需要運用數(shù)理統(tǒng)計方法。

    表3是某篇文章作者報告的關(guān)于土壤氡濃度、表面氡析出率和土壤含水率檢測數(shù)據(jù)。作者通過表3數(shù)據(jù)發(fā)現(xiàn)土壤表面氡析出率與土壤含水率存在著負(fù)相關(guān)性,即土壤含水率越大,氡析出率越小,這是正確的。但又認(rèn)為,土壤氡濃度與氡析出率沒有明顯的相關(guān)性,這是欠妥的。如果按表3的數(shù)據(jù)作土壤氡濃度和氡析出率之間的相關(guān)圖(見圖4),2者之間的相關(guān)性非常明顯。除2010年12月8日的氡濃度(3 740 Bq/m3)太高外,其他各點的負(fù)相關(guān)性都很明顯。反映出氡析出率越高,留在土壤中的氡濃度就越低。其相關(guān)系數(shù)已大于顯著性相關(guān)臨界值r0.05/2,6-2。

    表3 土壤氡濃度、表面氡析出率和土壤含水率的檢測數(shù)據(jù)

    圖4 土壤氡濃度和氡析出率之間的相關(guān)性

    其次,作者還認(rèn)為,土壤氡濃度與含水率沒有明顯的相關(guān)性,這也欠妥。如果按表3的數(shù)據(jù)作圖(見圖5),土壤氡濃度和土壤含水率之間的正相關(guān)性也很明顯, 隨著土壤含水率的增加,土壤中氡濃度也相應(yīng)增加。其實這一趨勢與作者曾經(jīng)發(fā)現(xiàn)的“土壤表面氡析出率與土壤含水率存在著負(fù)相關(guān)性”有密切聯(lián)系,即土壤含水率越大,氡析出率越小,使得土壤中氡濃度也相應(yīng)增加。同樣,其相關(guān)系數(shù)也已大于顯著性相關(guān)臨界值r0.05/2,6-2。

    圖5 土壤氡濃度和含水率之間的相關(guān)性

    2.5 測量項目之間相關(guān)性的分析

    環(huán)境樣品中不同測得量之間或不同同位素或核素活度之間的比值,在一定條件下具有相關(guān)性。當(dāng)時間或空間等環(huán)境條件發(fā)生變化后,這一比值也可能會隨之發(fā)生變化。通過測量它們之間比值及其變化,有助于輻射環(huán)境監(jiān)測結(jié)果合理性的評價。

    2.5.1 原野γ輻射劑量率與土壤中天然放射性核素40K、238U和232Th含量的關(guān)系

    扣除宇宙射線貢獻(xiàn)后的原野γ輻射水平,是由陸地土壤中的放射性核素產(chǎn)生。為了通過土壤中的放射性核素來估算γ輻射的空氣吸收劑量率,文獻(xiàn)[9]推薦過不同形式的Beck公式,當(dāng)土壤采樣布點具有充分的代表性時,估算陸地γ輻射水平的1993年Beck公式是:

    Dγ=CKfK+CUfU+CThfTh

    (3)

    式中:Dγ——γ輻射空氣吸收劑量率,nGy/h;CK、CU、CTh——40K、238U和232Th在土壤中的比活度,Bq/kg;fK、fU、fTh——40K、238U和232Th的劑量系數(shù),(nGy·h-1)/(Bq·kg-1),其值分別為0.041 7,0.462,0.060 4。

    由于在自然環(huán)境中實際情況與假設(shè)的條件并不完全相同,估算結(jié)果與實際結(jié)果會存在一定的差異。在基本滿足公式的條件下,若實測值與計算值相差超過30%,則至少有一種方法不可靠。表4為2000年4月—2002年3月在某核電站周圍9處測得的原野γ輻射的空氣吸收劑量率[10]。

    表4 原野γ輻射的空氣吸收劑量率實測結(jié)果(D實測) nGy/h

    為了核查上述實測結(jié)果D測,通過相應(yīng)9個測點土壤中40K、238U和232Th的比活度(Bq/kg),分別計算各點γ輻射的空氣吸收劑量率,并與表4的實測值進(jìn)行比較,結(jié)果示于表5。結(jié)果表明,2者相差未超過25%。

    表5 原野γ輻射的空氣吸收劑量率實測結(jié)果與計算結(jié)果的比較(D計算)

    2.5.2 土壤樣品中238U和226Ra比活度的關(guān)系

    238U和226Ra同屬238U衰變系列,在無外部因素作用的情況下,這2種放射性核素處于放射性平衡,如原礦石,其中238U和226Ra的比活度相等。

    自然界中土壤由成土母巖風(fēng)化形成,或由于河水沖積而成。盡管土壤形成初期階段238U和226Ra的平衡關(guān)系可能受到破壞,但如果土壤形成后長期未受到雨水、地表水或地下水侵蝕等外部因素的影響,這種平衡關(guān)系會被重新建立。

    對于地表土壤、水庫底泥或海底泥,由于雨水的沖刷和淋洗、地表水或海水的浸泡,使這種平衡關(guān)系受到不同程度的影響,因而這些樣品中238U和226Ra的比活度不一定相等。由于這些外在因素的作用差異很大,在不同地區(qū)的土壤樣品中238U和226Ra的比活度的比值也很不相同。但總的來說,238U的比活度高于226Ra的比活度。就地表土而言,雨水沖刷和淋洗越多,238U和226Ra的比活度可能差別越大。如果有酸雨的淋洗,也可能出現(xiàn)226Ra的比活度高于238U的情況。由于海底泥比水庫底泥所經(jīng)歷的化學(xué)過程更復(fù)雜,因此前者的比值一般比后者大。表6列出部分土壤中238U和226Ra的比活度及其比值。

    2.5.3 陸地介質(zhì)樣品中137Cs和90Sr比活度的差異

    137Cs和90Sr均是裂變產(chǎn)物,半衰期分別為30年和28年,環(huán)境中的殘留量來自全球進(jìn)行的核試驗的散落物。在核爆中137Cs的產(chǎn)額比90Sr高,產(chǎn)額比為1.6∶1.0。137Cs與土壤中的礦物質(zhì)結(jié)合得非常牢固,植物根部對其吸收很少,而90Sr則易于被吸收。因此地表土壤中137Cs的含量會明顯高于90Sr,

    表6 某核電站周圍部分土壤中238U和226Ra的比活度和比值(2000年)

    隨著時間的推移,這種差距會不斷擴(kuò)大。表7列出某核電站周圍部分土壤中137Cs和90Sr的比活度及其比值,為了比較,同時列出1個水庫底泥和1個海洋沉積物中137Cs和90Sr的比活度及其比值。由表7可見,表層土壤中137Cs的含量明顯高于90Sr,其比值為3.1~7.7。水庫底泥和海洋沉積物中137Cs和90Sr的比活度差別較小,可能是在水和海水的作用下,部分137Cs溶于水中,使137Cs的比活度減少,而90Sr在沉積物中的比活度卻變化不大。

    表7 某核電站周圍部分土壤和沉積物中137Cs和90Sr的比活度和比值 (2000—2001年)

    由于植物易于吸收90Sr,因此一般情況下,植物樣品中137Cs的含量會低于90Sr;但是另一方面,由于土壤中90Sr比137Cs衰減得快,因此當(dāng)土壤中的90Sr衰減到很小時,盡管植物易于吸收90Sr,但在長期未進(jìn)行核試驗和未發(fā)生核事故的情況下,隨著時間的推移,同種植物中90Sr的比活度和137Cs的比活度會逐漸接近,在某一時刻會相等,甚至?xí)l(fā)生反向變化,即90Sr的比活度低于137Cs的比活度。尤其是從水田中生長的水稻,由于稻田中的90Sr被水反復(fù)浸取和洗脫,其含量比137Cs減少得更快,所以常常發(fā)現(xiàn)不少地區(qū)的大米90Sr的比活度低于137Cs。 表8 列出部分植物中90Sr和137Cs的比活度及其比值[11]。

    表8 我國部分植物中90Sr和137Cs的比活度及其比值(1982—1985年)

    3 實驗室分析過程控制中的一些問題

    在分析實驗室管理中,要求首先應(yīng)識別并確定過程,然后及時發(fā)現(xiàn)和排除檢測過程中的異常變異,使上一過程的變異不要帶到下一過程中去,以保證過程的穩(wěn)定性和檢測質(zhì)量的一致性。在控制技術(shù)中,統(tǒng)計技術(shù)是識別、分析和控制異常變異的重要手段。用控制技術(shù)進(jìn)行過程控制稱為統(tǒng)計過程控制,其重要手段是控制圖,它一直是各實驗室普遍存在的薄弱環(huán)節(jié)。

    3.1 質(zhì)量控制圖的繪制

    質(zhì)量控制圖常用來反映分析測量過程中所用儀器和分析測量操作是否處于質(zhì)量控制狀態(tài),例如,對測量裝置長期穩(wěn)定性的檢驗,為繪制好控制圖,必須有一定數(shù)量的預(yù)備數(shù)據(jù)。按輻射環(huán)境監(jiān)測技術(shù)規(guī)范[12]規(guī)定,建立一張本底或效率的控制圖,需要20個以上的常規(guī)測量數(shù)據(jù),由這些數(shù)據(jù)計算平均值和標(biāo)準(zhǔn)差。然后,以計數(shù)率為縱坐標(biāo),以日期為橫坐標(biāo)作圖,在平均值處引一條水平線,再分別在平均值加、減兩倍標(biāo)準(zhǔn)差的地方引出兩條水平的控制線。以此控制圖檢查以后在該測量條件下獲得的本底或效率是否在控制范圍內(nèi)。當(dāng)發(fā)現(xiàn)異常時,對儀器進(jìn)行調(diào)整,重新繪制控制圖。按照類似的方法,同樣可以繪制出空白控制圖、準(zhǔn)確度控制圖、平行樣控制圖和回收率控制圖,等等。但是在實際工作中有時會發(fā)現(xiàn)如下問題:

    (1)一些單位為了應(yīng)付檢查,沒有按規(guī)范要求,在常規(guī)分析測量之前積累20個以上的測量數(shù)據(jù),繪制出控制圖。而是將儀器或分析測量運行中獲得的實際測量數(shù)據(jù)計算平均值和標(biāo)準(zhǔn)差,計算控制限。這種把運行中獲得的數(shù)據(jù)既當(dāng)“運動員”又當(dāng)“裁判員”的方法,除了制造“符合質(zhì)控要求”的假象之外,對質(zhì)量控制起不到任何作用;

    (2)有些“質(zhì)量控制圖”沒有控制限、警告限等界線。例如有一個平行樣偏差控制圖,實際上只是一個“平行樣偏差變化圖”。圖中的縱坐標(biāo)是“相對標(biāo)準(zhǔn)偏差”,含義不明確,應(yīng)是相對偏差;橫坐標(biāo)叫“平行樣個數(shù)”也不合適,應(yīng)是平行樣序號。例如偏差達(dá)18.3%的是第8號平行樣,而不是平行樣個數(shù)為8,因為每一次的平行樣個數(shù)只為2。

    3.2 化學(xué)回收率與加標(biāo)回收率

    在大多數(shù)的放化分析過程中,待測核素不可避免地會受到一定的損失。為此必須建立一種測定化學(xué)回收率的方法來估計這種損失,并對測量結(jié)果進(jìn)行校正。

    化學(xué)回收率又稱化學(xué)產(chǎn)額,其主要用于分析過程中對待測核素的損失量進(jìn)行校正,通常要求其不低于50%。但是,最重要的是要求回收率必須穩(wěn)定,同時產(chǎn)品必須達(dá)到所需要的化學(xué)純度與放射性純度。對于一個訓(xùn)練有素的分析人員,經(jīng)多次操作后回收率會穩(wěn)定在一定的范圍內(nèi)。如果某次分析中獲得的回收率太低,應(yīng)當(dāng)舍棄這次分析結(jié)果。反之,超過正常范圍的回收率也是不可取的,這是操作失誤或者在原始樣品中有載體物質(zhì)存在的先兆。雖然人們總是希望獲得高的化學(xué)回收率和高的產(chǎn)品純度,但是往往很難兩全其美。當(dāng)兩者發(fā)生矛盾時,常常是稍稍犧牲高的回收率來達(dá)到高的純度。

    測定化學(xué)回收率時要根據(jù)待測核素的特性和實驗室的條件選擇適當(dāng)?shù)姆椒āD壳俺S玫氖欠€(wěn)定性同位素載體法和放射性同位素作為示蹤劑的放射測量法。

    采用事前預(yù)定的固定回收率去校正測量結(jié)果會帶來很大的系統(tǒng)誤差。在分析土壤中238Pu時,初期由于缺乏238Pu的同位素,采用事前實驗確定的0.70作為固定回收率對分析結(jié)果進(jìn)行校正。后來獲得236Pu后,采用其做示蹤劑測定了628個試樣的全程回收率,并將其繪制成直方圖,結(jié)果示于圖6。由圖6可見,全程回收率近似正態(tài)分布,最小值為0.10,最大值為0.91,平均值為0.48,標(biāo)準(zhǔn)偏差0.20。如果采用原先預(yù)定的回收率0.70去校正分析結(jié)果,將會帶來很大的誤差。

    圖6 用236Pu做產(chǎn)額示蹤劑測得的238Pu回收率

    國內(nèi)環(huán)境放化分析領(lǐng)域較為流行一種叫“加標(biāo)回收率”,在一些質(zhì)量保證計劃中要求“加標(biāo)回收率控制在95%~105%以內(nèi)”。這種提法和要求不嚴(yán)謹(jǐn)和科學(xué)。在放射化學(xué)中,所謂“加標(biāo)回收率”,意即回收的放射性活度與加入的放射性活度之比值,將其控制在95%~105%以內(nèi)。在90Sr和137Cs分析中,一般用載體回收率或放射性示蹤劑的回收率表示,并不要求一定要達(dá)到95%~105%,只要不太低但穩(wěn)定即可。

    要求“加標(biāo)回收率控制在95%~105%以內(nèi)”,其“回收率”并不是真正的回收率,而是經(jīng)回收率校正后測得的放射性活度與加入的放射性活度之比值。這種要求,在國際放射化學(xué)領(lǐng)域比較罕見。在環(huán)境放化分析中,建議用準(zhǔn)確度(或與標(biāo)準(zhǔn)值的相對偏差)可能更好,“加標(biāo)回收率”95%~105%相當(dāng)于準(zhǔn)確度為±5%。

    3.3 分析測量結(jié)果的比對

    分析測量結(jié)果的比對包括實驗室內(nèi)不同方法、不同儀器和不同時期所得結(jié)果的內(nèi)部比對,以及不同實驗室間測量結(jié)果的外部比對。

    對測量結(jié)果有重要影響的計量器具和測量儀表,應(yīng)按計量認(rèn)證要求的周期進(jìn)行檢定、校準(zhǔn),在運行過程中繪制本底和探測效率質(zhì)控圖。還要特別注意,多數(shù)單位的輻射環(huán)境監(jiān)測的結(jié)果常常是用不同探測器、不同方法連續(xù)進(jìn)行多年獲得的。隨著儀器和方法的進(jìn)步,或儀器的更新?lián)Q代,改用新儀器或新方法是難免的,但其常給各年間的數(shù)據(jù)分析增加一個新的不一致因素。一個有效地防止結(jié)果不一致方法就是選擇一定數(shù)量的有代表性的點或保留樣品,作新、舊的比對測量,以保持?jǐn)?shù)據(jù)的連續(xù)可比性。如果在發(fā)現(xiàn)新、舊數(shù)據(jù)存在顯著差異后,才想到可能是儀器或元件、方法變化所致,再作這類比對測量已幾乎不可能了。因此,保存一部分樣品或建立相對穩(wěn)定的環(huán)境輻射場,對于保持監(jiān)測結(jié)果的穩(wěn)定性和溯源性是十分必要的。

    實驗室能力驗證是利用實驗室間比對確定實驗室的檢測能力的活動,實際上它是為確保實驗室維持較高的檢測水平而對其能力進(jìn)行考核、監(jiān)督和確認(rèn)的一種驗證活動。要想發(fā)現(xiàn)和確定實驗室分析測量所產(chǎn)生的系統(tǒng)誤差,僅靠實驗室自身的質(zhì)量控制,往往是十分困難的,參加能力驗證計劃,可為實驗室提供評價其出具數(shù)據(jù)可靠性和有效性的客觀證據(jù),從而增加公眾對其能力的信任。作為實驗室來說,還可通過能力驗證的外部措施,補(bǔ)充實驗室內(nèi)部的質(zhì)量控制程序。

    我國開展不同輻射環(huán)境測量實驗室間測量結(jié)果的比對工作大約可追溯到20世紀(jì)80年代,對提高實驗室檢測能力起到了重要作用。這類比對工作最好由權(quán)威計量機(jī)構(gòu)組織,要求比較嚴(yán)格,以達(dá)到實驗室的測量結(jié)果與國際權(quán)威實驗室間具有可比性。

    1985年以后,中國輻射防護(hù)研究院在李德平院士的指導(dǎo)下參加了由世界衛(wèi)生組織國際放射性參考中心(WHO IRC)、國際原子能機(jī)構(gòu)分析質(zhì)量控制服務(wù)處(IAEA AQCS)等部門組織的國際比對,受到國際組織的好評,被IRC 主席譽(yù)為一流實驗室。此后,中輻院低本底放射性實驗室的歷屆領(lǐng)導(dǎo),都堅持了這一傳統(tǒng),2000年以來每年繼續(xù)參加美國能源部環(huán)境測量實驗室質(zhì)保項目(EML-QAP)和愛達(dá)華放射性環(huán)境科學(xué)實驗室混合分析能力評估項目(RESL-MAPEP)質(zhì)控樣品測量與分析活動,比對樣品包括土壤、蔬菜和空氣濾膜等,對保持實驗室長期處于質(zhì)量控制狀態(tài),起到了很好的作用。

    參加國內(nèi)外的比對活動,要求參比單位和人員具有長遠(yuǎn)的眼光,不要為短期的功利主義或經(jīng)濟(jì)效益所左右。此外,比對的結(jié)果不太可能總是百分之百滿意,暴露出問題是好事,不要輕易地追查個人責(zé)任或與某種利益掛鉤,其不利于調(diào)動實驗室人員參加比對的積極性,相反還會滋長掩蓋問題或弄虛作假的壞作風(fēng)。

    曾經(jīng)發(fā)現(xiàn)國內(nèi)某行業(yè)組織的一次有幾十個實驗室參加的測量比對中,對開展測量比對的目的缺乏正確的認(rèn)識,不是為了發(fā)現(xiàn)問題,找到解決問題的方法,改進(jìn)今后的工作,而是為了顯示各單位的結(jié)果都很“滿意”,彼此“高度一致”。 從比對報告可以看出,比對測量并不是由各家獨立完成,而是“各單位監(jiān)測數(shù)據(jù)一律經(jīng)比對技術(shù)組一一校核匯總完成”。這有點像考生的考卷都要一一經(jīng)過監(jiān)考老師校核無誤后才能交卷一樣。難怪列出的結(jié)果都能做到100%的“滿意”!

    實驗室間測量結(jié)果的比對是分析過程控制的一種形式,其出發(fā)點是預(yù)防不合格,而不是為了證明過程的正確。2000年中國輻射防護(hù)研究院與某分包方共同承擔(dān)了核電站運行前本底調(diào)查任務(wù)。調(diào)查工作開始前,雙方進(jìn)行過土壤中90Sr的分析比對,但是采用的比對樣品是摻標(biāo)樣品,90Sr的比活度較高(0.75 Bq/g),分析的取樣量必然較少,因而實際土壤中大量存在的210Bi對分析的干擾未得到充分暴露。這說明當(dāng)時在指導(dǎo)思想上并沒有把發(fā)現(xiàn)問題和解決問題放在首位,而是把共同取得業(yè)主的認(rèn)可、通過比對證明雙方的方法和結(jié)果都沒有問題看得很重要。這種指導(dǎo)思想的不恰當(dāng)便造成實施方法不正確,因而掩蓋了矛盾,錯過了時機(jī)。后來雙方與中國原子能院合作,共同采集瞬時平行的土壤作為比對樣品,問題才被發(fā)現(xiàn)并得到解決。

    [1] 國務(wù)院環(huán)境保護(hù)委員會.國外輻射防護(hù)規(guī)程匯編5:環(huán)境放射性監(jiān)測規(guī)定[M].北京:中國環(huán)境科學(xué)出版社,1981:122-332.

    [2] 國防科學(xué)技術(shù)工業(yè)委員會. 電離輻射測量探測限和判斷閾的確定:EJ/T 1204—2006[S]. 北京:中國標(biāo)準(zhǔn)出版社,2006.

    [3] HERBERT L, VOLCHOK E D. EML Procedures Manual[M].26thEdition HASL-300,1983.

    [4] LOCHAMV J C.The minimum detectable concept[M]. NBS 75thanniversary symposium,NBS,1981.

    [5] 袁之倫,潘自強(qiáng), 張艷霞,等. 核電廠流出物低水平放射性核素監(jiān)測問題的探討[J]. 輻射防護(hù),2015,35(1):1-7.

    [6] 中國實驗室認(rèn)可委員會,化學(xué)分析中不確定度的評估指南[M].北京:中國計量出版社,2001.

    [7] 沙連茂,山本政儀, 小村和久,等.中國幾個地區(qū)土壤中239,240Pu、241Am和137Cs 等放射性核素的測定[J].環(huán)境科學(xué),1991:12(6):58.

    [8] 莫光華,黃乃明.環(huán)境輻射監(jiān)測數(shù)據(jù)的相關(guān)性[C]//核安全與輻射環(huán)境安全專業(yè)委員會,中國核學(xué)會輻射防護(hù)學(xué)會.全國放射性流出物和環(huán)境監(jiān)測與評價研討會論文匯編.杭州,2003:459-464.

    [9] 聯(lián)合國原子輻射效應(yīng)科學(xué)委員會.電離輻射源與效應(yīng)[M].北京:原子能出版社,1993.

    [10] 中國輻射防護(hù)研究院.田灣核電站申請資料許可證階段環(huán)境輻射本底調(diào)查報告[R].2002.

    [11] 張景源,諸洪達(dá).中國食品放射性及其所致內(nèi)劑量[M].北京:中國環(huán)境科學(xué)出版社,1989.

    [12] 國家環(huán)境保護(hù)總局.輻射環(huán)境監(jiān)測技術(shù)規(guī)范:HJ/T 61—2001[S].北京:中國標(biāo)準(zhǔn)出版社,2001.

    欄目編輯 李文峻

    Some Problems of Reasonableness Evaluation on Radiation Environmental Monitoring Data

    SHA Lian-mao

    (ChinaInstituteforRadiationProtection,Taiyuan,Shanxi030006,China)

    Some problems of reasonableness evaluation on radiation environmental monitoring data were discussed in this paper. The discussion would be centre on date near by limit of detection, trend analyses of monitoring results and analytical process control in laboratories.

    Environmental monitoring; Radiation; Reasonableness evaluation;Measured date

    10.3969/j.issn.1674-6732.2017.03.001

    2017-01-18

    沙連茂(1937—),男,研究員,大學(xué)本科,從事輻射環(huán)境監(jiān)測工作。

    X830.3;X837

    B

    1674-6732(2017)03-0001-10

    猜你喜歡
    比活度活度放射性
    2021年海陽核電站周邊飲用水中 總α、β放射性水平分析
    居里夫人發(fā)現(xiàn)放射性
    2020—2021年福建運行核電廠周圍生物中90Sr放射性水平調(diào)查與評價*
    淺析大理石礦的放射性檢測
    建材與裝飾(2018年5期)2018-02-13 23:12:02
    CaO-SiO2-FeO-P2O5-Al2O3脫磷渣系中組元活度的計算
    放射性家族個性有不同 十面埋“輻”你知多少
    核電廠惰性氣體排放活度濃度的估算
    來自放射性的電力
    太空探索(2015年10期)2015-07-18 10:59:20
    某地矽卡巖的γ能譜測量
    鉆井液處理劑溶液活度測量方法對比
    丝袜喷水一区| 人人妻人人澡人人看| 丰满人妻熟妇乱又伦精品不卡| 欧美日韩福利视频一区二区| 在线观看免费日韩欧美大片| 大香蕉久久网| 精品人妻熟女毛片av久久网站| 国产免费又黄又爽又色| 国产淫语在线视频| 亚洲国产精品999| videosex国产| 中国美女看黄片| 香蕉国产在线看| 国产麻豆69| av网站在线播放免费| 欧美另类一区| 亚洲国产欧美在线一区| 久久久精品94久久精品| 天堂中文最新版在线下载| 国产成人精品久久久久久| 国产女主播在线喷水免费视频网站| 老鸭窝网址在线观看| 久久天躁狠狠躁夜夜2o2o | 99久久精品国产亚洲精品| 国产成人av教育| 亚洲成国产人片在线观看| 国产精品免费大片| 国产精品一区二区免费欧美 | 2018国产大陆天天弄谢| 亚洲欧美一区二区三区久久| 99久久人妻综合| 99re6热这里在线精品视频| 亚洲av日韩在线播放| 少妇人妻 视频| 黑人欧美特级aaaaaa片| 亚洲精品第二区| 久久国产精品男人的天堂亚洲| 成人亚洲精品一区在线观看| 欧美日韩福利视频一区二区| 国产成人免费无遮挡视频| 欧美成人午夜精品| 国产午夜精品一二区理论片| 亚洲精品自拍成人| 国产精品免费大片| 一级黄色大片毛片| 精品一区二区三卡| 国产有黄有色有爽视频| 狂野欧美激情性xxxx| 91精品国产国语对白视频| 日日爽夜夜爽网站| 晚上一个人看的免费电影| 三上悠亚av全集在线观看| 国产精品国产av在线观看| 精品免费久久久久久久清纯 | 考比视频在线观看| 日韩av不卡免费在线播放| 一级片'在线观看视频| 国产精品久久久久成人av| 丝瓜视频免费看黄片| videosex国产| 操美女的视频在线观看| av电影中文网址| 丁香六月天网| 美女脱内裤让男人舔精品视频| 亚洲av成人精品一二三区| 高清不卡的av网站| 黄网站色视频无遮挡免费观看| 久久天堂一区二区三区四区| 久久久久久久久久久久大奶| 欧美久久黑人一区二区| 亚洲国产成人一精品久久久| 欧美日本中文国产一区发布| 免费人妻精品一区二区三区视频| 搡老乐熟女国产| 美女国产高潮福利片在线看| 久久精品aⅴ一区二区三区四区| 最近最新中文字幕大全免费视频 | 久久精品aⅴ一区二区三区四区| 男女下面插进去视频免费观看| 最新在线观看一区二区三区 | 蜜桃在线观看..| 伊人亚洲综合成人网| 成人三级做爰电影| 欧美精品亚洲一区二区| 亚洲男人天堂网一区| 亚洲精品久久久久久婷婷小说| 夫妻午夜视频| 一级黄色大片毛片| 亚洲第一青青草原| 日本色播在线视频| 高清不卡的av网站| 国产深夜福利视频在线观看| 亚洲人成77777在线视频| 99香蕉大伊视频| 天堂8中文在线网| 一边摸一边做爽爽视频免费| 免费看不卡的av| 十八禁人妻一区二区| 多毛熟女@视频| 亚洲精品日本国产第一区| 精品亚洲成国产av| 性色av乱码一区二区三区2| 日本vs欧美在线观看视频| 精品高清国产在线一区| 久久精品国产a三级三级三级| 视频区欧美日本亚洲| 99热国产这里只有精品6| 国产成人精品久久二区二区免费| 日本91视频免费播放| 午夜福利免费观看在线| 一区二区av电影网| netflix在线观看网站| 人人妻,人人澡人人爽秒播 | 大香蕉久久网| 日本wwww免费看| 亚洲国产欧美网| 精品人妻一区二区三区麻豆| cao死你这个sao货| 一级片'在线观看视频| 中文乱码字字幕精品一区二区三区| 少妇精品久久久久久久| 一级,二级,三级黄色视频| 激情五月婷婷亚洲| 韩国高清视频一区二区三区| 男女床上黄色一级片免费看| 后天国语完整版免费观看| 精品一区二区三区av网在线观看 | 久久99热这里只频精品6学生| 中文字幕制服av| 人妻一区二区av| 日本av免费视频播放| 亚洲,一卡二卡三卡| 91字幕亚洲| 亚洲av电影在线进入| 女人精品久久久久毛片| 人妻人人澡人人爽人人| 日韩电影二区| 久久天躁狠狠躁夜夜2o2o | 99精国产麻豆久久婷婷| 大码成人一级视频| 精品少妇一区二区三区视频日本电影| 成人免费观看视频高清| 一本—道久久a久久精品蜜桃钙片| 在现免费观看毛片| 黄片小视频在线播放| 久9热在线精品视频| 真人做人爱边吃奶动态| 亚洲图色成人| cao死你这个sao货| 黑人欧美特级aaaaaa片| 国产亚洲av片在线观看秒播厂| 成在线人永久免费视频| 夫妻午夜视频| 少妇粗大呻吟视频| 夫妻性生交免费视频一级片| 美女扒开内裤让男人捅视频| 美女视频免费永久观看网站| 岛国毛片在线播放| 女警被强在线播放| 男人爽女人下面视频在线观看| 男女午夜视频在线观看| 天天影视国产精品| 九色亚洲精品在线播放| 免费不卡黄色视频| 欧美精品av麻豆av| 日日夜夜操网爽| 一本综合久久免费| 欧美 亚洲 国产 日韩一| 国产精品久久久久久精品电影小说| 久久久久久久国产电影| 一区福利在线观看| 久久久国产精品麻豆| 日韩大码丰满熟妇| 在线av久久热| 丁香六月欧美| 亚洲欧美日韩另类电影网站| 欧美+亚洲+日韩+国产| av欧美777| 亚洲国产日韩一区二区| 亚洲av欧美aⅴ国产| 亚洲伊人久久精品综合| 国产老妇伦熟女老妇高清| 一区二区日韩欧美中文字幕| 汤姆久久久久久久影院中文字幕| netflix在线观看网站| 精品人妻一区二区三区麻豆| a级毛片在线看网站| 黄色a级毛片大全视频| 午夜视频精品福利| 免费看不卡的av| 视频区图区小说| 日韩免费高清中文字幕av| 久久精品熟女亚洲av麻豆精品| 免费高清在线观看视频在线观看| 无限看片的www在线观看| 国产精品香港三级国产av潘金莲 | 不卡av一区二区三区| 亚洲视频免费观看视频| 精品一区二区三区四区五区乱码 | 视频区欧美日本亚洲| 女人精品久久久久毛片| 欧美日韩黄片免| 国产高清不卡午夜福利| 熟女少妇亚洲综合色aaa.| 亚洲熟女毛片儿| 女性生殖器流出的白浆| 欧美另类一区| 亚洲精品国产av蜜桃| 精品免费久久久久久久清纯 | 国产一区二区三区综合在线观看| 久久久精品国产亚洲av高清涩受| 国产精品一区二区精品视频观看| 色94色欧美一区二区| 亚洲精品国产一区二区精华液| 建设人人有责人人尽责人人享有的| 国产一区二区激情短视频 | netflix在线观看网站| 免费看av在线观看网站| 精品亚洲成a人片在线观看| 在线av久久热| 亚洲午夜精品一区,二区,三区| 国产主播在线观看一区二区 | 国产片内射在线| 国产精品久久久久成人av| av天堂在线播放| 成人亚洲精品一区在线观看| 免费在线观看黄色视频的| 黄色a级毛片大全视频| 国产高清视频在线播放一区 | 久久鲁丝午夜福利片| 1024视频免费在线观看| 久久久久视频综合| 啦啦啦中文免费视频观看日本| 午夜两性在线视频| 丝袜喷水一区| 女性生殖器流出的白浆| 这个男人来自地球电影免费观看| 亚洲国产av新网站| 午夜福利,免费看| 国产女主播在线喷水免费视频网站| 777米奇影视久久| 最近手机中文字幕大全| 精品亚洲成a人片在线观看| av有码第一页| 人成视频在线观看免费观看| 又大又黄又爽视频免费| 在线观看国产h片| 大片免费播放器 马上看| 少妇粗大呻吟视频| av视频免费观看在线观看| 亚洲黑人精品在线| 各种免费的搞黄视频| 欧美精品啪啪一区二区三区 | 精品少妇久久久久久888优播| 极品少妇高潮喷水抽搐| 制服诱惑二区| 日本欧美国产在线视频| 一本一本久久a久久精品综合妖精| 国产麻豆69| 亚洲精品第二区| 国产熟女欧美一区二区| 国产精品一区二区在线不卡| 午夜影院在线不卡| 国产成人91sexporn| 国产精品熟女久久久久浪| 久久久久精品人妻al黑| 操出白浆在线播放| 91精品伊人久久大香线蕉| 日韩 亚洲 欧美在线| 中文字幕精品免费在线观看视频| 亚洲av在线观看美女高潮| 叶爱在线成人免费视频播放| 成人国语在线视频| 亚洲自偷自拍图片 自拍| 男人添女人高潮全过程视频| 汤姆久久久久久久影院中文字幕| 久久久国产一区二区| 男人爽女人下面视频在线观看| 日韩 欧美 亚洲 中文字幕| 狠狠婷婷综合久久久久久88av| 桃花免费在线播放| 97人妻天天添夜夜摸| 免费看不卡的av| 国产麻豆69| 丝袜脚勾引网站| 亚洲精品美女久久av网站| 色婷婷av一区二区三区视频| 精品福利观看| 久久精品人人爽人人爽视色| 2018国产大陆天天弄谢| 韩国高清视频一区二区三区| 国产精品一区二区在线观看99| 免费久久久久久久精品成人欧美视频| 国产精品久久久久久精品古装| 久久女婷五月综合色啪小说| 国产精品免费大片| 丝袜在线中文字幕| 你懂的网址亚洲精品在线观看| 激情五月婷婷亚洲| 涩涩av久久男人的天堂| 国产精品 国内视频| 亚洲欧洲精品一区二区精品久久久| 欧美亚洲日本最大视频资源| 亚洲,欧美,日韩| 超色免费av| 久久99热这里只频精品6学生| 下体分泌物呈黄色| 亚洲久久久国产精品| 人人妻人人澡人人看| 女性被躁到高潮视频| kizo精华| 又大又黄又爽视频免费| 午夜福利影视在线免费观看| 国产又爽黄色视频| 国产成人精品久久二区二区91| 欧美日韩亚洲国产一区二区在线观看 | 国产av精品麻豆| 亚洲一卡2卡3卡4卡5卡精品中文| 亚洲成国产人片在线观看| 美女福利国产在线| 久久人人97超碰香蕉20202| 亚洲精品美女久久久久99蜜臀 | 99热全是精品| 国产免费现黄频在线看| 亚洲国产av影院在线观看| 午夜日韩欧美国产| 精品国产一区二区三区四区第35| 电影成人av| 国产精品一区二区在线不卡| 欧美 日韩 精品 国产| 久久人妻福利社区极品人妻图片 | tube8黄色片| 国产熟女欧美一区二区| 国产人伦9x9x在线观看| 久久国产精品人妻蜜桃| 男人添女人高潮全过程视频| 国产男人的电影天堂91| 国产精品麻豆人妻色哟哟久久| 久久精品亚洲av国产电影网| 一边摸一边做爽爽视频免费| 麻豆av在线久日| 老司机亚洲免费影院| 视频区欧美日本亚洲| 精品一区二区三区av网在线观看 | 好男人视频免费观看在线| 一二三四社区在线视频社区8| 一级毛片电影观看| 人妻 亚洲 视频| 亚洲精品乱久久久久久| 国产在线观看jvid| 久久这里只有精品19| 丰满饥渴人妻一区二区三| 精品福利永久在线观看| 欧美精品亚洲一区二区| 纵有疾风起免费观看全集完整版| 亚洲欧美中文字幕日韩二区| 肉色欧美久久久久久久蜜桃| 色播在线永久视频| 色94色欧美一区二区| 精品少妇久久久久久888优播| 国产av一区二区精品久久| xxx大片免费视频| av在线app专区| 视频区图区小说| 亚洲国产欧美日韩在线播放| 欧美日韩成人在线一区二区| 亚洲成人免费电影在线观看 | 久久久精品区二区三区| 亚洲av国产av综合av卡| 一二三四在线观看免费中文在| 国产高清不卡午夜福利| 大片免费播放器 马上看| 超碰97精品在线观看| 久久精品亚洲av国产电影网| 国产亚洲一区二区精品| 亚洲一卡2卡3卡4卡5卡精品中文| 亚洲精品国产av蜜桃| 成年人黄色毛片网站| 国产亚洲一区二区精品| 男人操女人黄网站| 好男人电影高清在线观看| 精品国产乱码久久久久久小说| 男人舔女人的私密视频| 操美女的视频在线观看| 美女午夜性视频免费| 久久久国产一区二区| h视频一区二区三区| 国产激情久久老熟女| 80岁老熟妇乱子伦牲交| 久久人人爽av亚洲精品天堂| 国产欧美亚洲国产| 天天躁日日躁夜夜躁夜夜| 丝袜美腿诱惑在线| 老司机亚洲免费影院| 永久免费av网站大全| 国产精品一区二区精品视频观看| 久久精品亚洲av国产电影网| 999久久久国产精品视频| 午夜老司机福利片| 日本91视频免费播放| xxx大片免费视频| av欧美777| 国产片内射在线| 免费在线观看完整版高清| 高潮久久久久久久久久久不卡| 亚洲国产日韩一区二区| netflix在线观看网站| 一级毛片黄色毛片免费观看视频| 亚洲一区中文字幕在线| 免费看十八禁软件| 欧美 日韩 精品 国产| 在线观看免费午夜福利视频| 极品少妇高潮喷水抽搐| 熟女少妇亚洲综合色aaa.| 久久毛片免费看一区二区三区| 精品一区在线观看国产| 亚洲精品美女久久av网站| 美女福利国产在线| 成年美女黄网站色视频大全免费| 97人妻天天添夜夜摸| 超色免费av| 欧美日韩亚洲高清精品| 成人国产av品久久久| 国产精品免费视频内射| 91国产中文字幕| 美女视频免费永久观看网站| 色婷婷av一区二区三区视频| 日韩精品免费视频一区二区三区| 国产日韩一区二区三区精品不卡| www.av在线官网国产| 91精品伊人久久大香线蕉| 香蕉丝袜av| a级毛片黄视频| 日韩av不卡免费在线播放| 亚洲三区欧美一区| 久久精品久久久久久噜噜老黄| 老汉色∧v一级毛片| 欧美日韩一级在线毛片| tube8黄色片| 国产国语露脸激情在线看| 蜜桃国产av成人99| 亚洲成av片中文字幕在线观看| 国产黄频视频在线观看| 侵犯人妻中文字幕一二三四区| 欧美在线一区亚洲| 一级黄片播放器| 亚洲av日韩在线播放| 国产在视频线精品| 久久99一区二区三区| av国产精品久久久久影院| 亚洲人成77777在线视频| 午夜免费鲁丝| 亚洲人成77777在线视频| 99国产精品99久久久久| 久久99一区二区三区| 大片免费播放器 马上看| 搡老岳熟女国产| 激情五月婷婷亚洲| 亚洲av在线观看美女高潮| 这个男人来自地球电影免费观看| 免费在线观看影片大全网站 | 91成人精品电影| 又粗又硬又长又爽又黄的视频| 亚洲精品一区蜜桃| 国产精品欧美亚洲77777| 欧美在线黄色| 丰满人妻熟妇乱又伦精品不卡| 欧美黑人精品巨大| av在线播放精品| 男女床上黄色一级片免费看| 精品久久久精品久久久| 黑丝袜美女国产一区| 欧美日韩av久久| 99国产综合亚洲精品| av电影中文网址| 欧美日韩综合久久久久久| 丰满少妇做爰视频| 亚洲久久久国产精品| 男女下面插进去视频免费观看| 日韩中文字幕欧美一区二区 | 久久久久国产一级毛片高清牌| 男人舔女人的私密视频| 一边亲一边摸免费视频| 亚洲,欧美,日韩| svipshipincom国产片| 日韩,欧美,国产一区二区三区| 精品卡一卡二卡四卡免费| 国产又爽黄色视频| 国产成人一区二区三区免费视频网站 | 日本vs欧美在线观看视频| 亚洲精品久久成人aⅴ小说| 黄色怎么调成土黄色| 午夜免费男女啪啪视频观看| 丁香六月欧美| 欧美在线一区亚洲| 在线 av 中文字幕| 欧美成狂野欧美在线观看| 国产亚洲精品第一综合不卡| 精品一区二区三区四区五区乱码 | 亚洲av电影在线进入| 中文乱码字字幕精品一区二区三区| 久久九九热精品免费| 侵犯人妻中文字幕一二三四区| 国产极品粉嫩免费观看在线| 无遮挡黄片免费观看| 日本vs欧美在线观看视频| 成人午夜精彩视频在线观看| av一本久久久久| 国产三级黄色录像| xxx大片免费视频| 99re6热这里在线精品视频| 黑人猛操日本美女一级片| 女人久久www免费人成看片| av一本久久久久| 久久亚洲精品不卡| 一区二区三区四区激情视频| 亚洲久久久国产精品| 夜夜骑夜夜射夜夜干| 亚洲一码二码三码区别大吗| 午夜视频精品福利| 视频在线观看一区二区三区| 91字幕亚洲| 国产精品偷伦视频观看了| 好男人电影高清在线观看| 国产亚洲欧美精品永久| 久久狼人影院| 久久精品久久精品一区二区三区| 美女高潮到喷水免费观看| 在线 av 中文字幕| 国产黄色免费在线视频| 国产爽快片一区二区三区| 在线观看免费午夜福利视频| 亚洲 国产 在线| 黄色片一级片一级黄色片| 嫩草影视91久久| 国产精品三级大全| 国产av一区二区精品久久| 大片免费播放器 马上看| 欧美黄色片欧美黄色片| 9191精品国产免费久久| 国产主播在线观看一区二区 | 国产一级毛片在线| 久久久久国产一级毛片高清牌| 自拍欧美九色日韩亚洲蝌蚪91| 欧美黑人欧美精品刺激| 免费日韩欧美在线观看| 中文字幕人妻熟女乱码| 人妻 亚洲 视频| 女人精品久久久久毛片| 日韩精品免费视频一区二区三区| 少妇猛男粗大的猛烈进出视频| 久久狼人影院| 精品少妇内射三级| 成年美女黄网站色视频大全免费| 亚洲午夜精品一区,二区,三区| 国产极品粉嫩免费观看在线| 精品人妻1区二区| 国产一区二区 视频在线| 免费看不卡的av| 最近手机中文字幕大全| 电影成人av| 午夜免费成人在线视频| 久久影院123| 视频区欧美日本亚洲| 亚洲精品一卡2卡三卡4卡5卡 | 国产一区二区三区综合在线观看| 五月天丁香电影| 国产伦人伦偷精品视频| 啦啦啦啦在线视频资源| 亚洲美女黄色视频免费看| 日韩 欧美 亚洲 中文字幕| www.熟女人妻精品国产| 久久久久精品人妻al黑| 性少妇av在线| 久久av网站| 91老司机精品| 男女免费视频国产| 欧美日韩av久久| 性少妇av在线| 午夜福利,免费看| 国产成人免费无遮挡视频| 视频在线观看一区二区三区| 亚洲av综合色区一区| 国产亚洲欧美精品永久| 免费在线观看完整版高清| 欧美精品啪啪一区二区三区 | 91成人精品电影| 国产精品免费大片| 免费少妇av软件| av片东京热男人的天堂| 亚洲精品国产av蜜桃| 久久这里只有精品19| 女人被躁到高潮嗷嗷叫费观| 操美女的视频在线观看| 久久久久视频综合| 男人爽女人下面视频在线观看| a级毛片在线看网站| 精品久久久精品久久久| 一级,二级,三级黄色视频| 免费在线观看影片大全网站 | 一级片免费观看大全| 老司机深夜福利视频在线观看 | 啦啦啦啦在线视频资源| 成人黄色视频免费在线看| 日韩 亚洲 欧美在线| 只有这里有精品99| 欧美精品啪啪一区二区三区 | 亚洲一区中文字幕在线| 免费人妻精品一区二区三区视频| 嫩草影视91久久| 亚洲国产欧美在线一区| 国产又色又爽无遮挡免| 丝袜美足系列| 亚洲欧美精品自产自拍| 九色亚洲精品在线播放|