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

    基于信息擴(kuò)散理論的塵肺病風(fēng)險(xiǎn)評(píng)估研究

    2016-07-19 06:23:10束堯宸
    安全與環(huán)境工程 2016年3期
    關(guān)鍵詞:職業(yè)危害塵肺病風(fēng)險(xiǎn)評(píng)估

    束堯宸,邵 輝,李 展,趙 旭

    (常州大學(xué)環(huán)境與安全工程學(xué)院,江蘇 常州 213164)

    ?

    基于信息擴(kuò)散理論的塵肺病風(fēng)險(xiǎn)評(píng)估研究

    束堯宸,邵輝,李展,趙旭

    (常州大學(xué)環(huán)境與安全工程學(xué)院,江蘇 常州 213164)

    摘要:在進(jìn)行塵肺病職業(yè)危害風(fēng)險(xiǎn)評(píng)估時(shí),由于相關(guān)信息獲得難度大、樣本數(shù)量少,傳統(tǒng)的統(tǒng)計(jì)模型基本上無(wú)法滿足風(fēng)險(xiǎn)分析的需求?;谛畔U(kuò)散原理,通過(guò)擴(kuò)散函數(shù)對(duì)不完備信息進(jìn)行適當(dāng)?shù)臄U(kuò)充,建立起觀測(cè)樣本和風(fēng)險(xiǎn)概率分布的映射關(guān)系,構(gòu)建了評(píng)估區(qū)域塵肺病風(fēng)險(xiǎn)分析模型,并利用某省2009—2013年新發(fā)塵肺病病例數(shù)以及新增的Ⅰ、Ⅱ、Ⅲ期塵肺病病例數(shù)作為塵肺病風(fēng)險(xiǎn)分析樣本,應(yīng)用構(gòu)建的塵肺病風(fēng)險(xiǎn)分析模型,對(duì)上述各類塵肺病的潛在風(fēng)險(xiǎn)進(jìn)行評(píng)估,以為塵肺病預(yù)防和治理提供依據(jù)。

    關(guān)鍵詞:塵肺?。恍畔U(kuò)散理論;職業(yè)危害;風(fēng)險(xiǎn)評(píng)估

    塵肺病已經(jīng)成為職業(yè)病中影響最嚴(yán)重的一種職業(yè)危害,每年罹患?jí)m肺病的人數(shù)居高不下?,F(xiàn)階段對(duì)于塵肺病的研究大多處于理論探究、監(jiān)測(cè)與數(shù)值模擬相結(jié)合等方式[1-9],然而區(qū)域內(nèi)罹患?jí)m肺病是小樣本事件,想要預(yù)測(cè)粉塵的潛在職業(yè)危害性,事先給出一個(gè)合理的統(tǒng)計(jì)模型進(jìn)行風(fēng)險(xiǎn)評(píng)估幾乎是不可能的。以信息分配和信息擴(kuò)散為核心的模糊信息優(yōu)化處理技術(shù)是一門新興的數(shù)據(jù)處理技術(shù),可以彌補(bǔ)樣本信息不足的缺陷,但目前信息擴(kuò)散技術(shù)大多數(shù)應(yīng)用于自然災(zāi)害風(fēng)險(xiǎn)評(píng)估等方面,與職業(yè)危害方面相結(jié)合的研究并不多。本文基于信息擴(kuò)散理論,根據(jù)某省2009—2013年塵肺病的統(tǒng)計(jì)數(shù)據(jù),通過(guò)構(gòu)建塵肺病風(fēng)險(xiǎn)分析模型,綜合分析出該區(qū)域內(nèi)塵肺病的潛在危害性,并對(duì)其風(fēng)險(xiǎn)進(jìn)行了評(píng)估,以為塵肺病預(yù)防和治理提供依據(jù)。

    1信息擴(kuò)散的定義

    定義1:令X={x1,x2,…,xn}是一個(gè)隨機(jī)樣本,令U代表隨機(jī)變量X所有可能值的集,稱為X的論域,即U是隨機(jī)變量X的定義域。X和U分別簡(jiǎn)稱為樣本和論域[10]。

    定義2:令X是一個(gè)樣本,V是論域U的一個(gè)子集,從X×V到[0,1]的一個(gè)映射為

    μ:X×V→[0,1]

    (1)

    (2)

    稱為X在V上的一個(gè)信息擴(kuò)散,如果它是遞減的,即

    ?x∈X,?v′,v″∈V

    (3)

    如果‖v′-x‖≤‖v″-x‖,則μ(x,v′)≥μ(x,v″).μ稱為一個(gè)擴(kuò)散函數(shù),簡(jiǎn)稱為擴(kuò)散,V稱為一個(gè)監(jiān)控空間。由此我們得到:一個(gè)分配函數(shù)是一個(gè)擴(kuò)散函數(shù)。

    定義3:μ(x,u)是守恒的,當(dāng)且僅當(dāng)?x∈X,在其論域U上的積分值是1,即

    ∫Uμ(x,u)du=1

    (4)

    如果隨機(jī)變量的定義域U是離散的,U={u1,u2,…,um},則守恒條件為

    (5)

    2基于信息擴(kuò)散理論的塵肺病風(fēng)險(xiǎn)分析模型的建立

    信息擴(kuò)散是一種為了彌補(bǔ)信息不足,通過(guò)優(yōu)化利用小樣本信息,并對(duì)樣本進(jìn)行集值化模糊數(shù)學(xué)處理的方法,該方法可以將少量的觀測(cè)值樣本,擴(kuò)散為一個(gè)模糊集,即把單值樣本變成集值樣本。根據(jù)信息擴(kuò)散原理,在一個(gè)封閉系統(tǒng)中,信息擴(kuò)散函數(shù)是遞減的,也就是說(shuō)信息是沿著濃度減小的方向來(lái)移動(dòng)。分子的自然擴(kuò)散和信息擴(kuò)散都具有填補(bǔ)空白的功能,并且都是沿著濃度減小的方向進(jìn)行,因此當(dāng)信息擴(kuò)散不借助任何中介時(shí),可以認(rèn)為信息擴(kuò)散過(guò)程遵循分子擴(kuò)散過(guò)程的菲克第一定律和菲克第二定律[11-14]。在這個(gè)前提下,可得信息擴(kuò)散過(guò)程的正態(tài)擴(kuò)散函數(shù)f(ui)為

    (6)

    式中:h為正態(tài)擴(kuò)散的擴(kuò)散系數(shù),該擴(kuò)散系數(shù)是在研究正態(tài)擴(kuò)散的擇近原則時(shí)人為提出的參數(shù),其取值方式是根據(jù)黃金分割法推導(dǎo)得出的,在《自然災(zāi)害風(fēng)險(xiǎn)評(píng)價(jià)理論與實(shí)踐》[10]一書中對(duì)其有詳細(xì)的推導(dǎo)過(guò)程,由于篇幅所限在此不予證明,直接運(yùn)用即可。

    設(shè)Y={y1,y2,…,ym}為觀測(cè)樣本集合,是風(fēng)險(xiǎn)因素在n年里的實(shí)際觀測(cè)值的集合,指標(biāo)論域U={u1,u2,…,un}。

    根據(jù)式(6),第j個(gè)樣本的擴(kuò)散函數(shù)可寫為fj(ui),即

    (7)

    式中:fj(ui)為第j個(gè)樣本的擴(kuò)散函數(shù)(j=1,2,…,m;i=1,2,…,n)。

    (8)

    相應(yīng)的模糊子集的隸屬函數(shù)為

    (9)

    式(9)稱為樣本的信息歸一化分布,通過(guò)對(duì)其進(jìn)行處理,可以得到一種較好的風(fēng)險(xiǎn)評(píng)估結(jié)果,即令

    (10)

    式(10)的物理意義是:由觀測(cè)樣本集合Y={y1,y2,…,ym},經(jīng)信息擴(kuò)散推斷出,如果職業(yè)病觀測(cè)值只能取u1,u2,…,un中的一個(gè),那么在將yj均看作是樣本代表時(shí),觀測(cè)值為ui的樣本個(gè)數(shù)為q(ui)個(gè)。顯然,該值通常是一個(gè)正整數(shù),但一定是個(gè)不小于0的數(shù)字,令

    (11)

    式(11)就是樣本落在ui處的頻率值,可以作為概率估計(jì)值。超過(guò)ui的概率估計(jì)值,即為超越概率p(u≥ui),有

    (12)

    式(12)就是職業(yè)病風(fēng)險(xiǎn)估計(jì)值。

    3實(shí)例應(yīng)用與分析

    3.1塵肺病風(fēng)險(xiǎn)分析數(shù)據(jù)的調(diào)查統(tǒng)計(jì)

    某省以冶金、電子、機(jī)械等行業(yè)作為支柱產(chǎn)業(yè),近幾年塵肺病在該省份表現(xiàn)較為顯著。該區(qū)域塵肺病按其嚴(yán)重程度可以劃分為Ⅰ、Ⅱ、Ⅲ期,本文遂以該區(qū)域內(nèi)的塵肺職業(yè)病風(fēng)險(xiǎn)作為研究對(duì)象,在調(diào)研了大量的數(shù)據(jù)記錄以及相關(guān)的資料的情況下,以該省2009—2013年新發(fā)塵肺病病例數(shù)為樣本,統(tǒng)計(jì)分析出這五年該區(qū)域新發(fā)塵肺病病例數(shù)以及Ⅰ、Ⅱ、Ⅲ期塵肺病病例數(shù),詳見表1。

    3.2塵肺病風(fēng)險(xiǎn)分析觀測(cè)樣本的建立

    區(qū)域內(nèi)的塵肺職業(yè)危害風(fēng)險(xiǎn)評(píng)估是一個(gè)極其復(fù)雜的系統(tǒng)工程,涉及的因素多、變化大,可以利用的信息較少,因此本文采用信息擴(kuò)散理論構(gòu)建一個(gè)基于小樣本的區(qū)域塵肺病風(fēng)險(xiǎn)評(píng)估模型。根據(jù)表1統(tǒng)計(jì)得到的數(shù)據(jù),建立了4個(gè)觀測(cè)樣本,即新發(fā)塵肺病病例數(shù)和Ⅰ、Ⅱ、Ⅲ期新增塵肺病病例數(shù),這4個(gè)樣本分別為Y1、Y2、Y3、Y4,記為: Y1={760,1 120,753,803,615};Y2={623,982,632,663,498};Y3={102,89,70,95,76};Y4={35,49,51,45,41}。

    表1 某省2009—2013年新發(fā)塵肺病病例數(shù)

    指標(biāo)論域是根據(jù)觀測(cè)樣本數(shù)值的大小建立的,以包含所有觀測(cè)樣本。危險(xiǎn)等級(jí)(步長(zhǎng))的設(shè)定并沒(méi)有特定的標(biāo)準(zhǔn),在此要使觀測(cè)樣本點(diǎn)盡量多地落在控制點(diǎn)附近,以適當(dāng)簡(jiǎn)化計(jì)算工作量。

    根據(jù)樣本Y1,確定指標(biāo)論域U1[0,1 200],取步長(zhǎng)150,控制點(diǎn){0,150,300,450,600,750,900,1 050,1 200}共9個(gè),劃分9個(gè)風(fēng)險(xiǎn)水平;同理,根據(jù)樣本Y2,確定指標(biāo)論域U2[0,1 000],取步長(zhǎng)100,控制點(diǎn){0,100,200,300,400,500,600,700,800,900,1 000}共11個(gè),劃分11個(gè)風(fēng)險(xiǎn)水平;根據(jù)樣本Y3,確定指標(biāo)論域U3[0,105],取步長(zhǎng)15,控制點(diǎn){0,15,30,45,60,75,90,105}共8個(gè),劃分8個(gè)風(fēng)險(xiǎn)水平;根據(jù)樣本Y4,確定指標(biāo)論域U4[0,60],取步長(zhǎng)10,控制點(diǎn){0,10,20,30,40,50,60}共7個(gè),劃分7個(gè)風(fēng)險(xiǎn)水平。

    3.3擴(kuò)散系數(shù)h的取值

    根據(jù)樣本數(shù)據(jù),可以進(jìn)一步得到擴(kuò)散系數(shù)h的取值,見表2。

    表2 擴(kuò)散系數(shù)的取值

    注:表中n表示樣本容量;a表示樣本觀測(cè)值中的最小值;b表示樣本觀測(cè)值中的最大值;h為擴(kuò)散系數(shù)。

    樣本Y1的擴(kuò)散系數(shù)h1:樣本容量為5,樣本觀測(cè)值中最大值為1 120 ,最小值為615,根據(jù)表2計(jì)算,h1=0.814 6×(1 120-615)=411.373;同理,樣本Y2的擴(kuò)散系數(shù)h2:樣本容量為5,樣本觀測(cè)值中最大值為982,最小值為498,h2=0.814 6×(982-498)=394.27;樣本Y3的擴(kuò)散系數(shù)h3:樣本容量為5,樣本觀測(cè)值中最大值為102,最小值為70,h3=0.814 6×(102-70)=26.07;樣本Y4的擴(kuò)散系數(shù)h4:樣本容量為5,樣本觀測(cè)值中最大值為51,最小值為35,h4=0.814 6×(51-35)=13.033 6。

    3.4模糊隸屬函數(shù)計(jì)算

    根據(jù)公式(7)至(9),可得到4個(gè)樣本對(duì)應(yīng)的模糊隸屬函數(shù)如下:

    3.5塵肺病風(fēng)險(xiǎn)概率和超越概率的確定

    根據(jù)樣本(這五年來(lái)新患?jí)m肺病人的數(shù)量),可計(jì)算出2009—2013年新發(fā)塵肺病不同的風(fēng)險(xiǎn)等級(jí)及其所對(duì)應(yīng)的概率和超越概率,其計(jì)算結(jié)果見表3。

    表3 2009—2013年新發(fā)塵肺病不同的風(fēng)險(xiǎn)等級(jí)及其

    根據(jù)樣本Y2,可計(jì)算出2009—2013年新增Ⅰ期塵肺病風(fēng)險(xiǎn)等級(jí)及其所對(duì)應(yīng)的概率和超越概率,其計(jì)算結(jié)果見表4。

    表4 2009—2013年新增Ⅰ期塵肺病風(fēng)險(xiǎn)等級(jí)及其

    根據(jù)樣本Y3,可計(jì)算出2009—2013年新增Ⅱ期塵肺病風(fēng)險(xiǎn)等級(jí)及其所對(duì)應(yīng)的概率和超越概率,其計(jì)算結(jié)果見表5。

    表5 2009—2013年新增Ⅱ期塵肺病風(fēng)險(xiǎn)等級(jí)及其

    根據(jù)樣本Y4,可計(jì)算出2009—2013年新增Ⅲ期塵肺病風(fēng)險(xiǎn)等級(jí)及其對(duì)應(yīng)的概率和超越概率,其計(jì)算結(jié)果見表6。

    表6 2009—2013年新增Ⅲ期塵肺病風(fēng)險(xiǎn)等級(jí)及其

    3.6結(jié)果分析

    根據(jù)表3至表6,可得到該省近五年塵肺病的發(fā)生概率隨其風(fēng)險(xiǎn)(危害)等級(jí)的變化規(guī)律,見圖1至圖4。

    圖1 新發(fā)塵肺病的風(fēng)險(xiǎn)等級(jí)和概率分布圖Fig.1 Risk levels and probability distribution of the   newly diagnosed pneumoconiosis

    圖2 新增Ⅰ期塵肺病的風(fēng)險(xiǎn)等級(jí)和概率分布圖Fig.2 Risk levels and probability distribution of the   newly diagnosed phase I pneumoconiosis

    圖3 新增Ⅱ期塵肺病風(fēng)險(xiǎn)等級(jí)和概率分布圖Fig.3 Risk levels and probability distribution of the   newly diagnosed phase II pneumoconiosis

    圖4 新增Ⅲ期塵肺病風(fēng)險(xiǎn)等級(jí)和概率分布圖Fig.4 Risk levels and probability distribution of the   newly diagnosed phase III pneumoconiosis

    由圖1至圖4可見,該省新發(fā)塵肺病總?cè)藬?shù)以及新增Ⅰ、Ⅱ、Ⅲ期塵肺病的發(fā)生概率隨著其風(fēng)險(xiǎn)等級(jí)的上升均呈現(xiàn)出先上升后下降的趨勢(shì)。具體來(lái)看,圖1顯示出該省每年新發(fā)塵肺病Y1總?cè)藬?shù)在風(fēng)險(xiǎn)等級(jí)為450~1 200人之間的概率都在11%以上,在750~900人之間的概率達(dá)到了峰值,在1 200人左右其概率依舊維持在一個(gè)較高的水平上,因此應(yīng)該注重控制新發(fā)塵肺病總?cè)藬?shù),并通過(guò)相關(guān)措施將新發(fā)塵肺病人數(shù)控制在一個(gè)可接受的范圍內(nèi);圖2顯示出新增Ⅰ期塵肺病Y2病例人數(shù)在300~1 000人之間的概率達(dá)到了10%以上,在700人時(shí)的概率達(dá)到了峰值12.3%,由此看出新增Ⅰ期塵肺病病例人數(shù)在將來(lái)很大程度上將保持近幾年來(lái)的趨勢(shì);圖3則說(shuō)明了風(fēng)險(xiǎn)等級(jí)較高時(shí),新增Ⅱ期塵肺病Y3病例人數(shù)的概率也偏大,新增Ⅱ期塵肺病病例人數(shù)在風(fēng)險(xiǎn)等級(jí)為60~110人之間的概率在15%以上,在90人左右時(shí)的概率達(dá)到了峰值25.7%;圖4則說(shuō)明了新增Ⅲ期塵肺病Y4病例人數(shù)在風(fēng)險(xiǎn)等級(jí)為40~50人之間的概率最高,達(dá)到了將近30%,由于Ⅲ期塵肺病的危害程度最高,因此應(yīng)該將新發(fā)Ⅲ期塵肺病病例人數(shù)控制在一個(gè)合理的可接受范圍內(nèi)。

    本文將計(jì)算得到4組超越概率繪制成曲線,該省近五年的塵肺病職業(yè)危害以及新增Ⅰ、Ⅱ、Ⅲ期塵肺病風(fēng)險(xiǎn)等級(jí)與超越概率的關(guān)系見圖5和圖6。

    由圖5可見,該省新發(fā)塵肺病Y1的風(fēng)險(xiǎn)估計(jì)值隨著風(fēng)險(xiǎn)水平的增大而減小,每年新增塵肺病總?cè)藬?shù)在750人左右的概率為57.72%,在1 200人時(shí)的概率為11.42%;每年新增Ⅰ期塵肺病Y2風(fēng)險(xiǎn)估計(jì)值隨著風(fēng)險(xiǎn)水平的增大而減小,新增Ⅰ期塵肺病病例人數(shù)在0~600人之間的概率達(dá)到了55.66%,在900人以上的概率較小,約為20.45%。

    由圖6可見,該省新增Ⅱ期塵肺病Y3的風(fēng)險(xiǎn)估計(jì)值隨著風(fēng)險(xiǎn)水平的增大而減小,每年新增Ⅱ期塵肺病病例人數(shù)在0~75人之間的概率達(dá)到了70.87%,而在105人左右的概率為21.89%;新增Ⅲ期塵肺病Y4的風(fēng)險(xiǎn)估計(jì)值也是隨著風(fēng)險(xiǎn)水平的增大而降低,每年新增Ⅲ期塵肺病病例人數(shù)在0~40人之間的概率達(dá)到了73.42%,在60人左右的概率為16.9%。

    圖5 新發(fā)塵肺病和新增Ⅰ期塵肺病風(fēng)險(xiǎn)等級(jí)和超越   概率分布圖Fig.5 Risk levels and exceedance probability   distributions of the newly diagnosed and   phase I pneumoconiosis

    圖6 新增Ⅱ、Ⅲ期塵肺病風(fēng)險(xiǎn)等級(jí)和超越概率分布圖Fig.6 Risk levels and exceedance probability   distributions of the newly diagnosed and   phase Ⅱand Ⅲ pneumoconiosis

    由此可見,該省近五年塵肺病的發(fā)病現(xiàn)狀與通過(guò)信息擴(kuò)散理論構(gòu)建的數(shù)學(xué)模型所計(jì)算得到的結(jié)果保持了高度的一致性,基本符合實(shí)際情況。根據(jù)計(jì)算得到的結(jié)果,該省每年有50%以上的可能性產(chǎn)生750例左右的新發(fā)塵肺病患者,Ⅰ、Ⅱ、Ⅲ期塵肺病患的新增病例人數(shù)在600人、75人、40人這三個(gè)風(fēng)險(xiǎn)等級(jí)上的概率都在50%以上,說(shuō)明這三類新增塵肺病均處于一個(gè)較高的水平,特別是Ⅱ、Ⅲ期塵肺病為更為嚴(yán)重的塵肺病,因此這就提示我們要注重作業(yè)現(xiàn)場(chǎng)的粉塵控制,加強(qiáng)相關(guān)部門的監(jiān)管力度,做好除塵防塵措施以及為相應(yīng)的人員配備防塵用品等工作,把塵肺病風(fēng)險(xiǎn)降至最低水平,以確保作業(yè)人員的生命健康和安全。

    4結(jié)論

    (1) 基于區(qū)域內(nèi)的塵肺病病例是小樣本事件,要對(duì)區(qū)域內(nèi)的塵肺病風(fēng)險(xiǎn)進(jìn)行評(píng)估,現(xiàn)有的統(tǒng)計(jì)模型并不能滿足風(fēng)險(xiǎn)分析的需求,因此本文嘗試將職業(yè)衛(wèi)生問(wèn)題和信息擴(kuò)散方法相結(jié)合,對(duì)樣本點(diǎn)進(jìn)行集值化處理,將信息擴(kuò)散方法延伸到了職業(yè)衛(wèi)生問(wèn)題,構(gòu)建了評(píng)估區(qū)域塵肺病風(fēng)險(xiǎn)的分析模型。

    (2) 利用某省2009—2013年間的塵肺病統(tǒng)計(jì)資料,應(yīng)用構(gòu)建的塵肺病風(fēng)險(xiǎn)分析模型對(duì)該區(qū)域內(nèi)的塵肺病風(fēng)險(xiǎn)進(jìn)行了評(píng)估,結(jié)果表明:評(píng)估結(jié)果與實(shí)際情況基本保持一致,可以為該省的粉塵防治與控制提供參考。

    (3) 應(yīng)用結(jié)果證明將信息擴(kuò)散方法與職業(yè)衛(wèi)生問(wèn)題相結(jié)合是切實(shí)可行的,但是論域、步長(zhǎng)的選擇和風(fēng)險(xiǎn)等級(jí)的劃分對(duì)于結(jié)果產(chǎn)生的影響以及步長(zhǎng)的設(shè)置是否具有特定的標(biāo)準(zhǔn)還有待進(jìn)一步的研究。

    參考文獻(xiàn):

    [1] 張忠彬.作業(yè)場(chǎng)所職業(yè)危害風(fēng)險(xiǎn)綜合評(píng)價(jià)研究[J].中國(guó)安全生產(chǎn)科學(xué)技術(shù),2010,6(4):124-127.

    [2] 馬超.作業(yè)環(huán)境中粉塵的危害、傳播機(jī)理和防治措施[J].安防科技,2011(2):57-59.

    [3] 李明,吳超.粉塵點(diǎn)污染擴(kuò)散模型的可視化研究[J].環(huán)境科學(xué)與技術(shù),2006,11(29):12-14.

    [4] 魏宗蘇,魏傳光,金龍哲.全工班呼吸性粉塵測(cè)定及防塵措施研究[J].中國(guó)安全科學(xué)學(xué)報(bào),2008,18(9):106-110.

    [5]HuXM,WangB,FanBB.ResearchonKMC-basedevolutionsimulationofdustparticlesinvirtualcampusenvironment[J].Simulation Modelling Practice and Theory,2014,40:28-38.

    [6]XuZ,TravisJR,BreitungW,etal.Verificationofadusttransportmodelagainsttheoreticalsolutionsinmultidimensionaladvectiondiffusionproblems[J].Fusion Engineering and Design,2010,85:1935-1940.

    [7]QinY,WangT,ChenY,etal.Investigationandponderingofdustandpoisonhazardsin44smallwoodenfurnitureenterprises[J].Procedia Engineering,2012,43:484-488.

    [8]DatsonH,HallD,BirchB.Validationofanewmethodfordirectionaldustmonitoring[J].Atmospheric Environment,2012,50:1-8.

    [9]KhalajiM,RoshanzadehB,MansooriA,etal.Continuousdustmonitoringandanalysisbysparkinducedbreakdownspectroscopy[J].Optics and Lasers in Engineering,2012,50:110-113.

    [10]黃崇福.自然災(zāi)害風(fēng)險(xiǎn)評(píng)價(jià)[M].北京:科學(xué)出版社,2005.

    [11]邵輝,朱岳清,邵峰.基于信息擴(kuò)散的城市區(qū)域風(fēng)險(xiǎn)研究[J].中國(guó)安全科學(xué)學(xué)報(bào),2011,21(5):166-170.

    [12]楊旭,李春晨.基于信息擴(kuò)散理論的火災(zāi)風(fēng)險(xiǎn)評(píng)估模型研究及其應(yīng)用[J].工業(yè)安全與環(huán)保, 2010,36(1):42-43.

    [13]劉桂友,徐琳瑜.一種區(qū)域環(huán)境風(fēng)險(xiǎn)評(píng)價(jià)方法—信息擴(kuò)散法[J].環(huán)境科學(xué)學(xué)報(bào),2007,27(9):1549-2556.

    [14]王洪德,鄭玉錢.基于網(wǎng)格劃分及信息擴(kuò)散的化工園區(qū)安全風(fēng)險(xiǎn)評(píng)價(jià)技術(shù)[J].系統(tǒng)工程理論與實(shí)踐,2010,30(7):1286-1292.

    [15]鄭雙忠,鄧云峰,蔣清華.基于火災(zāi)統(tǒng)計(jì)災(zāi)情數(shù)據(jù)的城市火災(zāi)風(fēng)險(xiǎn)分析[J].中國(guó)安全生產(chǎn)科學(xué)技術(shù),2005,1(3):15-18.

    Pneumoconiosis Risk Assessment Based on Information Diffusion Theory

    SHU Yaochen,SHAO Hui,LI Zhan,ZHAO Xu

    (SchoolofEnvironment&SafetyEngineering,ChangzhouUniversity,Changzhou213164,China)

    Abstract:Due to the difficulty of acquiring relevant information and lack of samples in the occupational hazard risk evaluation of pneumoconiosis,traditional statistic models fail to meet the need for risk evaluation.Based on the principle of information diffusion,this paper uses spread function to extend the incomplete gathered information,establishes the mapping relationship between observation samples and risk probability distributions,and constructs the calculating model of risk evaluation of regional pneumoconiosis.By using the constructed model,the paper takes newly diagnosed pneumoconiosis patients and newly diagnosed phase I,II and III pneumoconiosis patients in a certain province from 2009—2013 as samples to evaluate the respective potential risk of above-mentioned pneumoconiosis,which provides reference for preventing and curing pneumoconiosis.

    Key words:pneumoconiosis;information diffusion theory;occupational hazards;risk assessment

    文章編號(hào):1671-1556(2016)03-0030-06

    收稿日期:2015-10-13修回日期:2016-02-01

    基金項(xiàng)目:國(guó)家安全監(jiān)管總局2013年安全生產(chǎn)重大事故防治關(guān)鍵技術(shù)科技項(xiàng)目[安監(jiān)總廳科技(2013)140號(hào)]

    作者簡(jiǎn)介:束堯宸(1991—),男,碩士研究生,主要研究方向?yàn)槌鞘泄舶踩c職業(yè)衛(wèi)生。E-mail:yaochenshu@126.com

    中圖分類號(hào):X18;R135

    文獻(xiàn)標(biāo)識(shí)碼:A

    DOI:10.13578/j.cnki.issn.1671-1556.2016.03.005

    通訊作者:邵輝(1955—),男,教授,主要從事安全工程等方面的教學(xué)與研究工作。E-mail:hshao@cczu.edu.cn

    猜你喜歡
    職業(yè)危害塵肺病風(fēng)險(xiǎn)評(píng)估
    國(guó)家塵肺病診療中心組建
    一例塵肺病患者的呼吸康復(fù)治療體會(huì)
    痛瀉要方加味湯治療塵肺病驗(yàn)案2則
    介入手術(shù)室護(hù)士職業(yè)危害因素分析及防護(hù)措施
    關(guān)于焦化廠職業(yè)危害分析及技術(shù)控制的幾點(diǎn)思考
    我國(guó)養(yǎng)老保險(xiǎn)基金投資運(yùn)營(yíng)的風(fēng)險(xiǎn)評(píng)估
    天然氣開采危險(xiǎn)有害因素辨識(shí)及風(fēng)險(xiǎn)評(píng)估管理
    現(xiàn)代風(fēng)險(xiǎn)導(dǎo)向?qū)徲?jì)局限性及其對(duì)策研究
    中小企業(yè)財(cái)務(wù)管理問(wèn)題研究
    寶石加工廠防塵設(shè)施綜合評(píng)價(jià)的層次—模糊綜合評(píng)價(jià)法及其應(yīng)用
    18禁黄网站禁片午夜丰满| 波多野结衣巨乳人妻| 韩国av一区二区三区四区| 757午夜福利合集在线观看| 老汉色∧v一级毛片| avwww免费| 久久中文字幕人妻熟女| 在线观看日韩欧美| 一级毛片精品| 999精品在线视频| 亚洲激情在线av| 90打野战视频偷拍视频| av福利片在线| 国产av麻豆久久久久久久| 国产黄片美女视频| 久久这里只有精品中国| 国产成人精品无人区| 欧美zozozo另类| 女警被强在线播放| 国产激情久久老熟女| 中文字幕久久专区| 1024视频免费在线观看| 99热这里只有精品一区 | 高潮久久久久久久久久久不卡| 婷婷精品国产亚洲av| 久久人人精品亚洲av| 欧美一区二区国产精品久久精品 | 不卡av一区二区三区| 成人国产一区最新在线观看| 老汉色av国产亚洲站长工具| 后天国语完整版免费观看| 成年女人毛片免费观看观看9| 久久热在线av| 麻豆久久精品国产亚洲av| 国产精品,欧美在线| 啪啪无遮挡十八禁网站| 亚洲第一欧美日韩一区二区三区| 香蕉久久夜色| 欧美在线黄色| 99国产精品一区二区蜜桃av| 亚洲avbb在线观看| 亚洲一区中文字幕在线| 两个人免费观看高清视频| 一夜夜www| 夜夜夜夜夜久久久久| 亚洲中文日韩欧美视频| 日本免费一区二区三区高清不卡| 国产野战对白在线观看| 美女大奶头视频| 国产蜜桃级精品一区二区三区| 日韩欧美一区二区三区在线观看| 丁香欧美五月| 97碰自拍视频| 妹子高潮喷水视频| 国产午夜福利久久久久久| 日韩欧美 国产精品| 麻豆成人av在线观看| 精品少妇一区二区三区视频日本电影| 搞女人的毛片| 波多野结衣巨乳人妻| 9191精品国产免费久久| 色综合婷婷激情| 国产精品一区二区免费欧美| 久久国产乱子伦精品免费另类| 亚洲成人精品中文字幕电影| 男人舔女人的私密视频| 不卡一级毛片| 日韩精品青青久久久久久| 免费在线观看影片大全网站| 老熟妇仑乱视频hdxx| 男女下面进入的视频免费午夜| 国产高清有码在线观看视频 | 91在线观看av| 99re在线观看精品视频| 亚洲五月婷婷丁香| 天堂av国产一区二区熟女人妻 | 99久久无色码亚洲精品果冻| 精品久久蜜臀av无| 一二三四在线观看免费中文在| 亚洲精品久久成人aⅴ小说| 日日干狠狠操夜夜爽| www国产在线视频色| 熟妇人妻久久中文字幕3abv| avwww免费| 午夜福利在线观看吧| 校园春色视频在线观看| 小说图片视频综合网站| 99精品欧美一区二区三区四区| 午夜福利在线观看吧| 久久人妻av系列| 三级国产精品欧美在线观看 | 欧美色视频一区免费| 看免费av毛片| 午夜福利18| 性色av乱码一区二区三区2| www国产在线视频色| 欧美激情久久久久久爽电影| 国产蜜桃级精品一区二区三区| 18禁裸乳无遮挡免费网站照片| 国产精品一区二区免费欧美| 午夜免费成人在线视频| 最近最新免费中文字幕在线| 精品国内亚洲2022精品成人| 日韩精品青青久久久久久| 久久中文字幕人妻熟女| 级片在线观看| 一区二区三区高清视频在线| videosex国产| 久久久国产成人精品二区| 久久午夜综合久久蜜桃| 欧美成人午夜精品| 亚洲国产欧美人成| 国产成人欧美在线观看| 在线观看免费视频日本深夜| 成人av一区二区三区在线看| 又粗又爽又猛毛片免费看| 国产一区二区在线观看日韩 | 婷婷精品国产亚洲av在线| 亚洲精品粉嫩美女一区| 黄色视频,在线免费观看| 可以在线观看的亚洲视频| 久久久久久久久免费视频了| 熟女电影av网| 国产又色又爽无遮挡免费看| 国产一区二区在线av高清观看| 午夜福利免费观看在线| 久久久久久久午夜电影| 狠狠狠狠99中文字幕| 日本一二三区视频观看| 18禁黄网站禁片午夜丰满| 国产激情久久老熟女| 波多野结衣巨乳人妻| 成年版毛片免费区| 男人舔女人的私密视频| 国产不卡一卡二| 1024香蕉在线观看| 51午夜福利影视在线观看| 国产欧美日韩一区二区精品| 欧美色视频一区免费| 国产爱豆传媒在线观看 | 两个人免费观看高清视频| 啦啦啦免费观看视频1| 久久精品国产清高在天天线| 亚洲熟妇中文字幕五十中出| 日本 av在线| 国产一区二区激情短视频| 男女之事视频高清在线观看| 亚洲五月天丁香| 午夜精品久久久久久毛片777| 99久久综合精品五月天人人| 91大片在线观看| 国产午夜精品论理片| 久久久久久人人人人人| 午夜福利欧美成人| 美女高潮喷水抽搐中文字幕| 很黄的视频免费| 老司机深夜福利视频在线观看| 在线永久观看黄色视频| 非洲黑人性xxxx精品又粗又长| 亚洲专区字幕在线| 欧美黑人欧美精品刺激| 欧美丝袜亚洲另类 | 人人妻人人澡欧美一区二区| 69av精品久久久久久| 午夜久久久久精精品| 两个人的视频大全免费| 日本a在线网址| 天堂影院成人在线观看| 女人被狂操c到高潮| 曰老女人黄片| 国产精品日韩av在线免费观看| 欧美丝袜亚洲另类 | 一本久久中文字幕| 免费在线观看黄色视频的| 午夜福利成人在线免费观看| 国产aⅴ精品一区二区三区波| 亚洲精品av麻豆狂野| 亚洲va日本ⅴa欧美va伊人久久| 麻豆久久精品国产亚洲av| 一区福利在线观看| 亚洲av美国av| 1024视频免费在线观看| 国产av麻豆久久久久久久| 国产成人系列免费观看| 给我免费播放毛片高清在线观看| 久久久久免费精品人妻一区二区| 亚洲,欧美精品.| 老司机午夜十八禁免费视频| av免费在线观看网站| 久久精品影院6| 久久草成人影院| 国产精品亚洲一级av第二区| 99国产综合亚洲精品| 国产成人啪精品午夜网站| 亚洲成人久久性| 淫秽高清视频在线观看| 制服人妻中文乱码| 搡老岳熟女国产| 在线观看免费日韩欧美大片| 欧美乱色亚洲激情| 丰满的人妻完整版| 国产日本99.免费观看| 亚洲午夜精品一区,二区,三区| 亚洲av美国av| 午夜成年电影在线免费观看| 99久久久亚洲精品蜜臀av| 国产亚洲精品久久久久5区| 天堂影院成人在线观看| 精华霜和精华液先用哪个| 在线观看午夜福利视频| 韩国av一区二区三区四区| 色综合亚洲欧美另类图片| 99精品久久久久人妻精品| 欧美乱码精品一区二区三区| 亚洲精品美女久久av网站| 国产97色在线日韩免费| 午夜视频精品福利| 一个人观看的视频www高清免费观看 | 青草久久国产| 老熟妇乱子伦视频在线观看| 亚洲九九香蕉| 亚洲欧美日韩高清在线视频| 午夜福利高清视频| 国产私拍福利视频在线观看| 久久久久久亚洲精品国产蜜桃av| 不卡一级毛片| 日本 欧美在线| 波多野结衣巨乳人妻| 男女做爰动态图高潮gif福利片| 日本一区二区免费在线视频| 亚洲精品在线美女| 1024视频免费在线观看| 成人高潮视频无遮挡免费网站| 国产精品一区二区三区四区久久| www.精华液| 国产精品综合久久久久久久免费| 欧美成人一区二区免费高清观看 | 欧美又色又爽又黄视频| 香蕉国产在线看| 老司机深夜福利视频在线观看| 亚洲18禁久久av| 成人三级做爰电影| 大型黄色视频在线免费观看| 制服诱惑二区| 久热爱精品视频在线9| 欧美日韩中文字幕国产精品一区二区三区| 男女那种视频在线观看| 美女大奶头视频| 欧美人与性动交α欧美精品济南到| 在线十欧美十亚洲十日本专区| 亚洲av成人av| 99热这里只有精品一区 | 在线播放国产精品三级| 国产一区二区在线av高清观看| 国产av一区二区精品久久| 男女床上黄色一级片免费看| www.999成人在线观看| 色精品久久人妻99蜜桃| 国产亚洲欧美98| 午夜福利在线观看吧| 国内精品久久久久久久电影| 真人做人爱边吃奶动态| 悠悠久久av| 国产av不卡久久| 日本一二三区视频观看| 精品一区二区三区av网在线观看| 99国产精品一区二区三区| 亚洲自拍偷在线| 男女之事视频高清在线观看| 亚洲国产欧美网| 两个人免费观看高清视频| 国产av又大| 欧美高清成人免费视频www| 久久久久久久久久黄片| 人妻久久中文字幕网| 此物有八面人人有两片| 亚洲欧美日韩高清专用| 熟女电影av网| 亚洲精品久久国产高清桃花| 老司机深夜福利视频在线观看| 国产精品九九99| 制服人妻中文乱码| 一a级毛片在线观看| 精品无人区乱码1区二区| 国产99白浆流出| 久久亚洲真实| 婷婷亚洲欧美| av视频在线观看入口| 亚洲中文日韩欧美视频| 亚洲中文av在线| 午夜久久久久精精品| 国产真人三级小视频在线观看| 日本一区二区免费在线视频| 12—13女人毛片做爰片一| 中亚洲国语对白在线视频| 无限看片的www在线观看| 色综合站精品国产| 亚洲 欧美一区二区三区| 午夜免费观看网址| 亚洲avbb在线观看| 男女下面进入的视频免费午夜| 免费在线观看成人毛片| 一进一出抽搐gif免费好疼| 日本 欧美在线| 亚洲欧美日韩高清在线视频| 午夜a级毛片| 制服丝袜大香蕉在线| 99精品在免费线老司机午夜| 五月玫瑰六月丁香| 观看免费一级毛片| 99热只有精品国产| 亚洲第一电影网av| 两个人免费观看高清视频| 久久久精品国产亚洲av高清涩受| 中文在线观看免费www的网站 | 中文字幕久久专区| 老司机深夜福利视频在线观看| avwww免费| 精品久久久久久久末码| 观看免费一级毛片| 日韩精品青青久久久久久| 天堂√8在线中文| 久久精品人妻少妇| 国产av不卡久久| 亚洲人成伊人成综合网2020| 在线观看舔阴道视频| 亚洲乱码一区二区免费版| 午夜福利免费观看在线| 日本撒尿小便嘘嘘汇集6| 久久久久国产精品人妻aⅴ院| 国产成人系列免费观看| 麻豆av在线久日| 50天的宝宝边吃奶边哭怎么回事| 日韩免费av在线播放| 欧美日韩亚洲综合一区二区三区_| 中文字幕精品亚洲无线码一区| 久久久久国产精品人妻aⅴ院| 在线永久观看黄色视频| 麻豆国产av国片精品| 成人永久免费在线观看视频| 亚洲中文av在线| 人妻丰满熟妇av一区二区三区| 亚洲激情在线av| 丰满人妻熟妇乱又伦精品不卡| 午夜a级毛片| 色av中文字幕| 国产真人三级小视频在线观看| 久久久国产成人精品二区| 亚洲精品一区av在线观看| 国产成人av教育| 免费在线观看影片大全网站| 在线观看免费午夜福利视频| 嫩草影视91久久| 99精品久久久久人妻精品| 老汉色∧v一级毛片| 午夜精品在线福利| 亚洲 欧美 日韩 在线 免费| 亚洲欧美日韩高清专用| 国内少妇人妻偷人精品xxx网站 | 亚洲欧美精品综合久久99| 国产三级在线视频| 最新在线观看一区二区三区| 村上凉子中文字幕在线| 欧美黄色淫秽网站| 国产爱豆传媒在线观看 | 精华霜和精华液先用哪个| 国产单亲对白刺激| 在线观看一区二区三区| 国产成人啪精品午夜网站| 国产主播在线观看一区二区| 国产精品 欧美亚洲| 日韩精品免费视频一区二区三区| 熟女少妇亚洲综合色aaa.| 一a级毛片在线观看| 精品欧美一区二区三区在线| 亚洲欧美日韩无卡精品| 久久久久性生活片| 午夜视频精品福利| 欧美色视频一区免费| 在线观看www视频免费| 又大又爽又粗| 999久久久精品免费观看国产| a在线观看视频网站| 日韩中文字幕欧美一区二区| 国产av又大| 三级国产精品欧美在线观看 | 亚洲国产中文字幕在线视频| 男女做爰动态图高潮gif福利片| 国产久久久一区二区三区| 国产又色又爽无遮挡免费看| 国产av一区在线观看免费| 91麻豆精品激情在线观看国产| 级片在线观看| 岛国视频午夜一区免费看| 亚洲人成网站高清观看| 亚洲国产欧美人成| 久久久久免费精品人妻一区二区| 51午夜福利影视在线观看| 欧美日韩瑟瑟在线播放| 夜夜爽天天搞| 一级黄色大片毛片| av超薄肉色丝袜交足视频| 亚洲av成人一区二区三| 香蕉av资源在线| 久久久久久亚洲精品国产蜜桃av| 精品高清国产在线一区| av中文乱码字幕在线| 亚洲色图 男人天堂 中文字幕| 级片在线观看| 人妻夜夜爽99麻豆av| 久久人妻福利社区极品人妻图片| 精品第一国产精品| 久久欧美精品欧美久久欧美| 久久久水蜜桃国产精品网| 久久久久国产一级毛片高清牌| 熟女电影av网| 丝袜美腿诱惑在线| 欧美极品一区二区三区四区| 一进一出抽搐gif免费好疼| 听说在线观看完整版免费高清| 99久久精品国产亚洲精品| 一二三四在线观看免费中文在| 女生性感内裤真人,穿戴方法视频| 久久久久国产一级毛片高清牌| 熟女电影av网| 天堂√8在线中文| 婷婷丁香在线五月| 美女午夜性视频免费| 亚洲男人的天堂狠狠| 男女做爰动态图高潮gif福利片| 黄色丝袜av网址大全| 欧美黑人欧美精品刺激| 波多野结衣巨乳人妻| 久久久精品欧美日韩精品| 搞女人的毛片| 最近最新免费中文字幕在线| 精品一区二区三区视频在线观看免费| 精品熟女少妇八av免费久了| 神马国产精品三级电影在线观看 | www日本黄色视频网| 亚洲avbb在线观看| 深夜精品福利| 国产1区2区3区精品| av在线天堂中文字幕| 成人高潮视频无遮挡免费网站| 国产精品久久久久久人妻精品电影| 白带黄色成豆腐渣| 一本一本综合久久| 婷婷精品国产亚洲av在线| 91麻豆av在线| 免费电影在线观看免费观看| 欧美一级a爱片免费观看看 | 激情在线观看视频在线高清| 亚洲五月婷婷丁香| 嫩草影视91久久| 性欧美人与动物交配| 日韩精品免费视频一区二区三区| 露出奶头的视频| 一级毛片精品| а√天堂www在线а√下载| 国产黄片美女视频| 精品欧美国产一区二区三| 男人舔奶头视频| 久久草成人影院| 日韩欧美一区二区三区在线观看| 淫妇啪啪啪对白视频| 国产精品日韩av在线免费观看| 欧美乱色亚洲激情| 在线观看午夜福利视频| 天天躁夜夜躁狠狠躁躁| 老司机午夜福利在线观看视频| 少妇被粗大的猛进出69影院| 久久久精品大字幕| 免费看美女性在线毛片视频| 黄色a级毛片大全视频| 国产一区二区在线av高清观看| 国产午夜精品论理片| 大型av网站在线播放| 亚洲 欧美 日韩 在线 免费| 别揉我奶头~嗯~啊~动态视频| av免费在线观看网站| 亚洲色图av天堂| 听说在线观看完整版免费高清| 一级毛片女人18水好多| 香蕉av资源在线| 熟女少妇亚洲综合色aaa.| 丁香欧美五月| e午夜精品久久久久久久| 久久久久久人人人人人| 欧美乱妇无乱码| 国产亚洲av高清不卡| 人成视频在线观看免费观看| 757午夜福利合集在线观看| 午夜亚洲福利在线播放| 国产伦在线观看视频一区| 淫秽高清视频在线观看| 国产1区2区3区精品| 亚洲中文av在线| 国产精品电影一区二区三区| 国产av不卡久久| 50天的宝宝边吃奶边哭怎么回事| 亚洲,欧美精品.| 天堂影院成人在线观看| 国产主播在线观看一区二区| 国产精品电影一区二区三区| 国产成人精品久久二区二区91| 精品国产亚洲在线| 欧美3d第一页| 久久天堂一区二区三区四区| 午夜日韩欧美国产| 最近最新免费中文字幕在线| 日本免费一区二区三区高清不卡| 校园春色视频在线观看| 国产av不卡久久| 天天一区二区日本电影三级| 两性夫妻黄色片| 女人高潮潮喷娇喘18禁视频| 国产av一区在线观看免费| 女生性感内裤真人,穿戴方法视频| 国产精品乱码一区二三区的特点| 99riav亚洲国产免费| 此物有八面人人有两片| 国产久久久一区二区三区| 国产黄a三级三级三级人| 欧美性长视频在线观看| 两性夫妻黄色片| 亚洲欧美日韩无卡精品| 久久国产精品人妻蜜桃| 听说在线观看完整版免费高清| av免费在线观看网站| 狂野欧美白嫩少妇大欣赏| av在线播放免费不卡| 午夜精品一区二区三区免费看| 国产精品精品国产色婷婷| 三级男女做爰猛烈吃奶摸视频| 在线观看日韩欧美| 色噜噜av男人的天堂激情| 精品国内亚洲2022精品成人| 国产高清有码在线观看视频 | 黄色视频不卡| 岛国在线观看网站| 国产真实乱freesex| 一级a爱片免费观看的视频| 亚洲在线自拍视频| 精品免费久久久久久久清纯| 亚洲电影在线观看av| 亚洲免费av在线视频| 91麻豆精品激情在线观看国产| 亚洲午夜理论影院| 桃红色精品国产亚洲av| 免费看a级黄色片| 91大片在线观看| 亚洲九九香蕉| 身体一侧抽搐| 人妻久久中文字幕网| 九色成人免费人妻av| 亚洲午夜精品一区,二区,三区| 国产真实乱freesex| 久久亚洲精品不卡| 在线国产一区二区在线| 日韩高清综合在线| 午夜影院日韩av| 夜夜躁狠狠躁天天躁| 狠狠狠狠99中文字幕| 国产爱豆传媒在线观看 | 变态另类丝袜制服| 手机成人av网站| 成人av在线播放网站| 亚洲国产看品久久| 国产av麻豆久久久久久久| 亚洲五月天丁香| 亚洲精品在线观看二区| 欧美久久黑人一区二区| 男女之事视频高清在线观看| 久久久久久人人人人人| svipshipincom国产片| 视频区欧美日本亚洲| 一二三四在线观看免费中文在| 国内精品久久久久久久电影| 亚洲,欧美精品.| 免费在线观看视频国产中文字幕亚洲| 久久久国产成人免费| 久久精品夜夜夜夜夜久久蜜豆 | 在线免费观看的www视频| 无限看片的www在线观看| 欧美色视频一区免费| av欧美777| 久久精品国产综合久久久| 在线观看午夜福利视频| 在线视频色国产色| 国产精品电影一区二区三区| 精品久久久久久久末码| 国产野战对白在线观看| 夜夜躁狠狠躁天天躁| 啪啪无遮挡十八禁网站| 国产精品影院久久| 日本一本二区三区精品| 岛国在线免费视频观看| 精品电影一区二区在线| 日本三级黄在线观看| 免费看美女性在线毛片视频| 麻豆久久精品国产亚洲av| 国产久久久一区二区三区| 俄罗斯特黄特色一大片| 老汉色∧v一级毛片| 日韩欧美在线二视频| 最近视频中文字幕2019在线8| 国产一区二区三区视频了| 99热这里只有精品一区 | 99久久无色码亚洲精品果冻| 午夜福利在线在线| 精品久久久久久久久久免费视频| 欧美zozozo另类| 国产男靠女视频免费网站| 国产精品一区二区免费欧美|