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

    基于DBN的掘進(jìn)工作面煤塵爆炸事故風(fēng)險預(yù)測

    2023-02-21 09:14:04李泓鯤王孝東陳鐘輝
    中國煤炭 2023年1期
    關(guān)鍵詞:煤塵后驗通風(fēng)

    李泓鯤,王孝東,李 陽,肖 露,陳鐘輝,劉 杰

    (1.昆明理工大學(xué)國土資源工程學(xué)院,云南省昆明市,650093;2.大唐華銀電力股份有限公司耒陽分公司,湖南省衡陽市,421800;3.云南農(nóng)業(yè)大學(xué)經(jīng)濟管理學(xué)院,云南省昆明市,650201;4.昆明理工大學(xué)公共安全與應(yīng)急管理學(xué)院,云南省昆明市,650093)

    0 引言

    在煤礦掘進(jìn)工作面生產(chǎn)施工過程中,往往會產(chǎn)生大量的煤塵。這些煤塵是煤礦掘進(jìn)工作面作業(yè)的嚴(yán)重災(zāi)害源之一,具有不可忽視的潛在爆炸風(fēng)險,嚴(yán)重威脅煤礦掘進(jìn)工作面的安全生產(chǎn)[1]。因此,研究煤礦掘進(jìn)工作面煤塵爆炸事故,尋找造成其煤塵爆炸事故的根本因素,具有重要意義。景國勛等[2]通過研究煤質(zhì)指標(biāo)對煤塵爆炸的影響,利用20 L球形爆炸容器對具有不同工業(yè)成分的煤塵爆炸最大壓力上升速率進(jìn)行測定,建立煤質(zhì)指標(biāo)與爆炸最大壓力上升速率的計算模型;裴蓓等[3]通過對煤塵爆炸初期與瓦斯的耦合爆炸災(zāi)害強化的產(chǎn)生機制探究,認(rèn)為初始爆炸強度對瓦斯/煤塵復(fù)合體系有重要影響;ZHANG Leilin等[4]通過建立不同彎曲角度的管道爆炸系統(tǒng),研究了管道彎曲角度對煤塵爆炸傳播規(guī)律的影響,認(rèn)為爆炸后沖擊波的峰值超壓急劇增大,拐角后又緩慢減小。

    目前隨著人工智能技術(shù)的發(fā)展,動態(tài)貝葉斯網(wǎng)絡(luò)(DBN)被大量運用到事故風(fēng)險分析和預(yù)測方面。王文和等[5]利用DBN對城鎮(zhèn)燃?xì)夤芫W(wǎng)進(jìn)行動態(tài)分析,認(rèn)為該方法克服了傳統(tǒng)風(fēng)險分析方法的不足,可動態(tài)反映燃?xì)夤芫W(wǎng)失效和事故后果發(fā)生概率隨時間變化的特征。同時,DBN還可用于水下連接器故障診斷、煤礦安全態(tài)勢感知[6-8]等諸多方面。

    由于前人對煤塵爆炸的相關(guān)分析多集中于煤塵爆炸機理實驗探究或單一因素對煤塵爆炸的影響[1-4],而對于造成煤塵爆炸事故風(fēng)險分析及預(yù)測的相關(guān)研究較少,因此其分析結(jié)果具有一定局限性。DBN作為目前被廣泛應(yīng)用于事故風(fēng)險分析和預(yù)測的手段之一,通過運用條件概率和時間序列能夠很好反映造成事故的各個風(fēng)險因素概率隨時間變化的情況,因此筆者將通過DBN對掘進(jìn)工作面煤塵爆炸事故進(jìn)行分析,尋找掘進(jìn)工作面煤塵爆炸事故中應(yīng)重點關(guān)注的基本風(fēng)險因素,為預(yù)防掘進(jìn)工作面煤塵爆炸事故提供客觀科學(xué)的依據(jù)。

    1 事故風(fēng)險分析方法

    1.1 動態(tài)貝葉斯網(wǎng)絡(luò)

    DBN是靜態(tài)貝葉斯網(wǎng)絡(luò)(BN)的擴展,可以容納時間序列數(shù)據(jù),適用于對動態(tài)系統(tǒng)進(jìn)行建模[9]。DBN主要由初始網(wǎng)絡(luò)U0和狀態(tài)轉(zhuǎn)移網(wǎng)絡(luò)Ut構(gòu)成[10],初始網(wǎng)絡(luò)為BN,而狀態(tài)轉(zhuǎn)移網(wǎng)絡(luò)則是2個相鄰時間片上網(wǎng)絡(luò)模型。

    由于掘進(jìn)工作面煤塵爆炸事故的發(fā)生是一個動態(tài)過程,且掘進(jìn)工作面煤塵爆炸事故具有不確定性、不完整性和時間相關(guān)性的特點,因此選用DBN作為掘進(jìn)工作面煤塵爆炸事故風(fēng)險分析方法。對于DBN事故分析方法,通常有以下2條假設(shè)[11]:

    (1)穩(wěn)定性。對于所有時間切片t,其轉(zhuǎn)移概率PU(t)(Xt|Xt-1)總是相同的。

    (2)馬爾科夫性。對于t+1時間切片中的節(jié)點,其節(jié)點狀態(tài)僅與t時間切片中節(jié)點狀態(tài)有關(guān),與t-1時間切片中節(jié)點狀態(tài)無關(guān)。

    對于滿足上述2個條件所構(gòu)建的用于掘進(jìn)工作面煤塵爆炸事故分析的DBN中,任意風(fēng)險因素概率分布公式為:

    (1)

    式中:P(Xi(t))——t時刻風(fēng)險因素的概率;

    P(Ut|Xi(t))——任意風(fēng)險因素在轉(zhuǎn)移網(wǎng)絡(luò)中的條件概率。

    DBN推理的依據(jù)與BN推理依據(jù)相同,均依據(jù)貝葉斯公式進(jìn)行推理:

    (2)

    1.2 DBN模型要素

    由于影響掘進(jìn)工作面煤塵爆炸事故的風(fēng)險因素有很多,在結(jié)合國家礦山安全監(jiān)察局(原國家煤礦安全監(jiān)察局)關(guān)于山東能源肥城礦業(yè)集團梁寶寺能源有限責(zé)任公司“8·20”煤塵爆炸事故通報[12]與前人相關(guān)分析的基礎(chǔ)上[13],從煤塵性質(zhì)、環(huán)境因素影響、管理缺陷、通風(fēng)時間與通風(fēng)設(shè)計不合理這4個方面對掘進(jìn)工作面煤塵爆炸事故基本風(fēng)險因素進(jìn)行選取,共選取15個造成掘進(jìn)工作面煤塵爆炸事故的基本風(fēng)險因素,分別為:煤塵粒徑分布X0、煤塵揮發(fā)分X1、煤塵爆炸性指數(shù)X2、煤塵水分及灰分X3、煤塵濃度在爆炸范圍X4、瓦斯?jié)舛萖5、工作面風(fēng)速不合理X6、氧氣濃度X7、環(huán)境火源X8、通風(fēng)時間不合理X9、局部通風(fēng)方案設(shè)計不合理X10、整體通風(fēng)方案設(shè)計不合理X11、除塵措施未落實X12、無相關(guān)安全培訓(xùn)X13、規(guī)章制度不健全X14,并根據(jù)事故因果關(guān)系建立初始貝葉斯網(wǎng)絡(luò)。將掘進(jìn)工作面煤塵爆炸事故作為初始貝葉斯網(wǎng)絡(luò)葉節(jié)點,二級風(fēng)險因素作為初始貝葉斯網(wǎng)絡(luò)中間節(jié)點,15種基本風(fēng)險因素作為初始貝葉斯網(wǎng)絡(luò)根節(jié)點。

    初始貝葉斯網(wǎng)絡(luò)中的中間節(jié)點與根節(jié)點的設(shè)置狀態(tài)如下:

    (3)

    將中間節(jié)點與根節(jié)點設(shè)置為2種狀態(tài),用于簡化網(wǎng)絡(luò)和減少中間節(jié)點數(shù)量以及優(yōu)化網(wǎng)絡(luò)計算。

    將葉節(jié)點T按發(fā)生事故可能造成的后果嚴(yán)重程度設(shè)置狀態(tài)如下:

    (4)

    將葉節(jié)點T按事故后果的5個等級設(shè)置為3種狀態(tài),貼近于實際情況,同時可簡化計算過程,減小樣本所需數(shù)據(jù)量,以實現(xiàn)DBN在有限信息條件下的定量推理。

    1.3 基于D-S證據(jù)理論與G1法確定專家權(quán)重

    在應(yīng)用DBN對掘進(jìn)工作面煤塵爆炸事故進(jìn)行推理時,需要大量事故演變過程的數(shù)據(jù)或?qū)<揖C合評價結(jié)果作為推理的先驗依據(jù)。由于缺乏大量的掘進(jìn)工作面煤塵爆炸事故演變過程數(shù)據(jù),且直接運用專家綜合評價結(jié)果作為推理的先驗依據(jù)主觀性太強,因此將基于專家綜合評價結(jié)果,應(yīng)用D-S證據(jù)理論預(yù)處理專家綜合評價結(jié)果,再通過應(yīng)用G1算法將預(yù)處理結(jié)果轉(zhuǎn)化為專家權(quán)重,根據(jù)相應(yīng)的專家權(quán)重確定掘進(jìn)工作面煤塵爆炸事故DBN先驗依據(jù)。

    對D-S證據(jù)理論簡要描述如下[14-15]:

    設(shè)Ω為一個識別框架,在識別框架上的基本概率賦值2Ω→[0,1]的函數(shù)m,成為mass函數(shù),并且滿足:

    (5)

    對于任意B∈Ω,其識別框架內(nèi)的有限個mass函數(shù)m1…mn的合成規(guī)則為:

    (6)

    式中:K——歸一化常數(shù)。

    利用D-S證據(jù)理論處理后,應(yīng)用G1算法確定專家權(quán)重步驟如下[16]:

    (1)確定評價指標(biāo)重要性排序Cn;

    (2)確定評價指標(biāo)相對重要程度Ri,其中Ri為Cn-1與Cn之比;

    (3)計算專家權(quán)重Vj:

    (7)

    上述過程均由煤礦安全領(lǐng)域相關(guān)專家進(jìn)行確定和評價,式(7)中n取15。3類不同專家的權(quán)重結(jié)果計算見表1。

    表1 專家權(quán)重說明

    2 動態(tài)風(fēng)險分析

    2.1 節(jié)點初始先驗概率

    在掘進(jìn)工作面煤塵爆炸事故DBN中,其各節(jié)點的條件概率用來描述各節(jié)點之間的關(guān)聯(lián)程度。根節(jié)點先驗概率是將專家綜合評價所得分?jǐn)?shù)轉(zhuǎn)化為專家權(quán)重并依據(jù)式(8)根節(jié)點先驗概率計算公式計算得出:

    (8)

    i——風(fēng)險因素的第i個狀態(tài);

    n——風(fēng)險因素的狀態(tài)個數(shù);

    j——風(fēng)險因素個數(shù);

    wijk——風(fēng)險因素j的狀態(tài)i發(fā)生的第k個專家對應(yīng)的權(quán)值系數(shù)。

    中間節(jié)點先驗概率由專家綜合評價結(jié)果轉(zhuǎn)化為專家權(quán)重后,結(jié)合式(9)聯(lián)合概率公式計算得出:

    (9)

    式中:P(Y)——中間節(jié)點先驗概率;

    P(Y|An)——An條件下發(fā)生Y的概率。

    對于掘進(jìn)工作面煤塵爆炸事故DBN,其各節(jié)點在初始網(wǎng)絡(luò)的先驗概率稱為初始先驗概率。掘進(jìn)工作面煤塵爆炸事故DBN初始先驗概率如圖1所示。

    圖1 掘進(jìn)工作面煤塵爆炸事故DBN初始先驗概率

    由圖1可知,掘進(jìn)工作面煤塵爆炸事故造成嚴(yán)重后果的初始先驗概率為0.090,造成較輕微后果的先驗概率為0.280,后果可忽略不計的先驗概率為0.630。對于各根節(jié)點先驗概率,最大的為煤塵濃度在爆炸范圍X4(0.282),最小的為無相關(guān)安全培訓(xùn)X13(0.114)。對圖1中各節(jié)點先驗概率數(shù)據(jù)進(jìn)行分析可知,掘進(jìn)工作面煤塵爆炸事故嚴(yán)重后果,輕微后果與后果忽略不計的概率分布與海因里希安全法則——即大量后果忽略不計和輕微后果的意外事件必將導(dǎo)致重大傷亡事故發(fā)生相符合。同時,依據(jù)先驗概率初步判斷影響掘進(jìn)工作面煤塵爆炸事故發(fā)生的最大二級節(jié)點為環(huán)境因素影響B(tài)2,最小為管理缺陷B4。在4個不同二級風(fēng)險因素下,對影響各個二級風(fēng)險因素發(fā)生的最大的基本風(fēng)險因素分別為煤塵爆炸性指數(shù)、煤塵濃度在爆炸范圍、局部通風(fēng)方案設(shè)計不合理、防塵措施未落實。

    從初始先驗概率角度對這4個基本風(fēng)險因素進(jìn)行分析,煤塵爆炸性指數(shù)越高,煤塵爆炸的風(fēng)險越大。而掘進(jìn)工作面煤塵濃度如果在爆炸范圍,往往是導(dǎo)致掘進(jìn)工作面發(fā)生煤塵爆炸事故的直接原因之一。局部通風(fēng)方案設(shè)計一般是在掘進(jìn)工作面施工之前就已經(jīng)完成,但由于掘進(jìn)工作面施工是一個動態(tài)過程,因此局部通風(fēng)方案設(shè)計與掘進(jìn)工作面實際生產(chǎn)施工之間存在滯后性,進(jìn)而導(dǎo)致局部通風(fēng)方案設(shè)計不合理,從而可能導(dǎo)致掘進(jìn)工作面煤塵爆炸事故的發(fā)生。防塵措施未落實會進(jìn)一步導(dǎo)致掘進(jìn)工作面附近煤塵積聚,在一定條件下可引發(fā)煤塵爆炸事故。

    2.2 節(jié)點先驗概率預(yù)測

    由于掘進(jìn)工作面煤塵爆炸事故DBN中的各節(jié)點概率會隨時間發(fā)生變化,且各節(jié)點的轉(zhuǎn)移過程符合馬爾科夫過程[17],因此對于掘進(jìn)工作面煤塵爆炸事故的先驗概率預(yù)測結(jié)果,可將專家經(jīng)驗轉(zhuǎn)化為時間權(quán)重結(jié)合馬爾科夫過程進(jìn)行判定[18-19]。掘進(jìn)工作面煤塵爆炸事故DBN轉(zhuǎn)移概率公式為:

    (10)

    式中:t——時間片;

    Si——節(jié)點X處于狀態(tài)i的個數(shù);

    Si,k——節(jié)點X有狀態(tài)i到狀態(tài)k的個數(shù)。

    時間權(quán)重計算公式為:

    (12)

    式中:τt——第t個時刻時間權(quán)重;

    δt——時間折現(xiàn)系數(shù)。

    將節(jié)點的轉(zhuǎn)移過程設(shè)置為6個時間片(t=0~5),葉節(jié)點T即掘進(jìn)工作面煤塵爆炸事故的3種狀態(tài)在6個時間片內(nèi)概率變化過程,如圖2所示。

    圖2 葉節(jié)點T的3種狀態(tài)概率變化

    由圖2可以看出,在6個時間片內(nèi),事故造成可忽略不計后果和較輕微后果的狀態(tài)概率均有所下降。其中造成可忽略不計后果的狀態(tài)概率在t=0~1時間片內(nèi)顯著下降,下降了0.128;在t=1~5時間片內(nèi)下降緩慢,下降了0.019,并最終趨于0.485這一定值。造成較輕微后果的狀態(tài)概率在t=0~1時間片內(nèi)下降了0.008;在t=1~5時間片內(nèi)下降了0.014,并最終趨于0.261這一定值。造成嚴(yán)重后果的狀態(tài)概率則在t=0~1時間片內(nèi)顯著上升,上升了0.112;在t=1~5時間片內(nèi)上升緩慢,上升了0.050,并最終趨于0.254這一定值。同時,與不同狀態(tài)的初始先驗概率進(jìn)行對比,可知掘進(jìn)工作面煤塵爆炸事故發(fā)生嚴(yán)重后果的狀態(tài)概率相較于初始先驗概率上升了0.160,發(fā)生較輕微后果的狀態(tài)概率相較于初始先驗概率下降0.022,事故后果忽略不計的狀態(tài)概率相較于初始先驗概率下降了0.147,這與圖2中所示變化趨勢相對應(yīng)。

    利用DBN正向預(yù)測各根節(jié)點事件(X0~X14)發(fā)生先驗概率在6個時間片后的結(jié)果分別為0.221、0.181、0.268、0.259、0.293、0.304、0.352、0.301、0.350、0.323、0.329、0.285、0.274、0.170、0.154。預(yù)測對掘進(jìn)工作面煤塵爆炸事故造成嚴(yán)重后果和較輕微后果的狀態(tài)概率影響最大的2個基本風(fēng)險因素是工作面風(fēng)速不合理(0.352)和環(huán)境火源(0.350),影響最小的基本風(fēng)險因素是規(guī)章制度不健全(0.154)。在4個不同二級風(fēng)險因素下,對各個二級風(fēng)險因素事件發(fā)生概率影響最大的基本風(fēng)險因素分別是煤塵爆炸性指數(shù)(0.268)、工作面風(fēng)速不合理(0.352)、局部通風(fēng)方案設(shè)計不合理(0.329)、防塵措施未落實(0.274)。

    通過對圖2與根節(jié)點先驗概率結(jié)果的分析,綜合考慮各基本風(fēng)險因素預(yù)測結(jié)果與掘進(jìn)工作面煤塵爆炸事故3種狀態(tài)概率變化情況,初步確定對掘進(jìn)工作面煤塵爆炸事故造成嚴(yán)重后果影響最大的6個基本風(fēng)險因素分別為煤塵爆炸性指數(shù)、工作面風(fēng)速不合理、瓦斯?jié)舛取h(huán)境火源、局部通風(fēng)方案設(shè)計不合理、防塵措施未落實。

    3 事故預(yù)防決策分析

    3.1 根節(jié)點后驗概率預(yù)測

    直接由初始先驗概率經(jīng)轉(zhuǎn)移網(wǎng)絡(luò)預(yù)測所得出的各個風(fēng)險因素的事件發(fā)生概率結(jié)果仍具有一定的主觀性,為進(jìn)一步減小主觀性和對預(yù)測結(jié)果進(jìn)行驗證,利用后驗概率公式并結(jié)合DBN進(jìn)行逆推理,得出各根節(jié)點后驗概率并對預(yù)測結(jié)果進(jìn)行修正。DBN根節(jié)點后驗概率計算公式為:

    (13)

    根據(jù)圖2所示,在未采取任何預(yù)防措施前提下,掘進(jìn)工作面煤塵爆炸事故后果忽略不計狀態(tài)概率與較輕微概率均有所下降,而造成嚴(yán)重后果狀態(tài)概率將上升。因此后驗概率主要討論掘進(jìn)工作面煤塵爆炸事故造成嚴(yán)重后果狀態(tài)時各節(jié)點后驗概率。

    將葉節(jié)點即掘進(jìn)工作面煤塵爆炸事故T造成嚴(yán)重后果狀態(tài)概率值設(shè)置為1.000,結(jié)合式(10)~式(13),預(yù)測各根節(jié)點(X0~X14)在6個時間片(t=0~5)后在煤塵爆炸事故T造成嚴(yán)重后果的狀態(tài)下,事件發(fā)生的后驗概率結(jié)果分別為0.231、0.183、0.272、0.260、0.289、0.302、0.348、0.304、0.353、0.318、0.332、0.288、0.271、0.173、0.154。

    由節(jié)點后驗概率結(jié)果可知,根節(jié)點煤塵粒徑分布X0與整體通風(fēng)方案設(shè)計不合理X11后驗概率相較于先驗概率上升了0.010,局部通風(fēng)方案設(shè)計不合理X10與環(huán)境火源X8后驗概率相較于先驗概率提升了0.003。各根節(jié)點后驗概率預(yù)測值與先驗概率預(yù)測結(jié)果基本一致,驗證了先驗概率對于各根節(jié)點預(yù)測結(jié)果的準(zhǔn)確性。

    綜合對各根節(jié)點先驗概率與后驗概率預(yù)測結(jié)果分析,對初步選取的6個基本風(fēng)險因素進(jìn)行修正,最終確定對掘進(jìn)工作面煤塵爆炸事故造成嚴(yán)重后果的5種基本風(fēng)險因素分別為瓦斯?jié)舛?、環(huán)境火源、通風(fēng)時間不合理、局部通風(fēng)方案設(shè)計不合理與防塵措施未落實。

    3.2 單因素作用分析

    由于DBN后驗概率公式具有雙向性,因此通過DBN后驗概率公式的雙向性可知,葉節(jié)點T在單個基本風(fēng)險因素發(fā)生條件下的后驗概率計算公式為:

    (14)

    由DBN推理雙向性,可通過已獲得的最終確定的5種基本風(fēng)險因素先驗概率、后驗概率,結(jié)合式(10)~式(12)和式(14),對在單因素發(fā)生條件下導(dǎo)致掘進(jìn)工作面煤塵爆炸事故造成嚴(yán)重后果的后驗概率(P(T=0|Xi=0))結(jié)果進(jìn)行計算和預(yù)測,6個時間片內(nèi)預(yù)測結(jié)果:X5為0.252、X8為0.251、X9為0.254、X10為0.255、X12為0.253。不同單因素作用條件下在6個時間片內(nèi)造成掘進(jìn)工作面煤塵爆炸事故嚴(yán)重后果的概率變化過程如圖3所示。

    圖3 單因素發(fā)生條件下嚴(yán)重后果后驗概率變化情況

    綜合圖3和6個時間片內(nèi)預(yù)測結(jié)果可知,對所選取的5種基本風(fēng)險因素,在單個風(fēng)險因素發(fā)生的條件下,預(yù)測在6個時間片后,掘進(jìn)工作面煤塵爆炸嚴(yán)重后果后驗概率最大為0.255,最小為0.251。且掘進(jìn)工作面煤塵爆炸事故造成嚴(yán)重后果的后驗概率變化趨勢均是初期先急劇上升(t=0~2),后緩慢上升(t=2~4),最后趨于某一定值(t=4~5)。局部通風(fēng)方案設(shè)計不合理的發(fā)生對于掘進(jìn)工作面煤塵爆炸事故造成嚴(yán)重后果的影響最大,使其后驗概率預(yù)測結(jié)果相對先驗概率預(yù)測結(jié)果提升了0.005。通過DBN雙向綜合分析,最終確定對掘進(jìn)工作面煤塵爆炸事故造成嚴(yán)重后果狀態(tài)影響最大的基本風(fēng)險因素為局部通風(fēng)方案設(shè)計不合理。

    3.3 多因素耦合作用分析

    不同組合條件下多個基本風(fēng)險因素耦合作用結(jié)果計算公式:

    (15)

    由于掘進(jìn)工作面煤塵爆炸事故造成嚴(yán)重后果往往是多個基本風(fēng)險因素耦合作用的結(jié)果,且每個基本風(fēng)險因素是否發(fā)生都是相互獨立的。因此將所選取的5種基本風(fēng)險因素進(jìn)行組合,結(jié)合DBN雙向推理性,通過已獲得的最終確定的5種基本風(fēng)險因素的先驗概率、后驗概率,并結(jié)合式(10)~式(12)與式(15),對多因素耦合作用條件下掘進(jìn)工作面煤塵爆炸事故造成嚴(yán)重后果狀態(tài)的后驗概率進(jìn)行計算與預(yù)測,6個時間片內(nèi)不同組合條件下的多種基本風(fēng)險因素耦合作用導(dǎo)致掘進(jìn)工作面煤塵爆炸事故造成嚴(yán)重后果的預(yù)測結(jié)果的最大值見表2,基本風(fēng)險因素耦合作用結(jié)果最大值的概率變化過程如圖4所示。

    表2 不同組合條件下基本風(fēng)險因素耦合作用結(jié)果最大值

    圖4 基本風(fēng)險因素耦合作用結(jié)果最大值的概率變化

    由表2可知,預(yù)測在6個時間片內(nèi),2種因素組合條件下通過基本風(fēng)險因素耦合作用,導(dǎo)致掘進(jìn)工作面煤塵爆炸事故發(fā)生嚴(yán)重后果后驗概率最小,為0.256;在5種因素組合條件下通過基本風(fēng)險因素耦合作用,導(dǎo)致掘進(jìn)工作面煤塵爆炸事故發(fā)生嚴(yán)重后果后驗概率最大,為0.266。由圖4可知,不同組合條件下基本風(fēng)險因素耦合作用導(dǎo)致掘進(jìn)工作面煤塵爆炸事故發(fā)生嚴(yán)重后果的初始后驗概率也略有不同,其中5種因素組合條件下初始后驗概率最大,而2種因素組合條件下初始后驗概率最小。

    由表2和圖4綜合分析可知,掘進(jìn)工作面煤塵爆炸事故發(fā)生嚴(yán)重后果的后驗概率變化情況均是初期先急劇上升(t=0~2),后緩慢上升(t=2~4),最后趨于某一定值(t=4~5)。因此對掘進(jìn)工作面煤塵爆炸事故,其事故預(yù)防重點應(yīng)是在事故發(fā)展初期對瓦斯?jié)舛?、環(huán)境火源、通風(fēng)時間不合理、局部通風(fēng)方案設(shè)計不合理與防塵措施未落實這5種基本風(fēng)險因素采取相關(guān)措施,防止這5種基本風(fēng)險因素同時發(fā)生。若在事故初期無法采取相關(guān)措施防止這5種基本風(fēng)險因素同時發(fā)生,則應(yīng)在事故發(fā)展初期采取相關(guān)措施防止表2中所列出的其他3種因素組合中的基本風(fēng)險因素同時發(fā)生,以降低掘進(jìn)工作面煤塵爆炸事故最終發(fā)展成嚴(yán)重后果的風(fēng)險。

    4 結(jié)論

    (1)結(jié)合實際生產(chǎn)事故案例,建立掘進(jìn)工作面煤塵爆炸DBN事故分析模型,并利用DBN進(jìn)行正向推理,得出掘進(jìn)工作面煤塵爆炸事故3種不同狀態(tài)的初始先驗概率分別為0.090、0.280、0.630;由初始先驗概率結(jié)果結(jié)合式(10)~式(12)計算可知,在6個時間片后3種狀態(tài)概率分別為0.254、0.261、0.485。

    (2)利用DBN進(jìn)行逆向推理,最終確定對掘進(jìn)工作面煤塵爆炸事故造成嚴(yán)重后果影響最大的5種基本風(fēng)險因素分別為瓦斯?jié)舛?、環(huán)境火源、通風(fēng)時間不合理、局部通風(fēng)方案設(shè)計不合理與防塵措施未落實。由計算可知,其后驗概率預(yù)測結(jié)果與先驗概率預(yù)測結(jié)果基本一致,分別為0.302、0.353、0.318、0.332、0.271。同時,在單因素作用條件下,最終確定對掘進(jìn)工作面煤塵爆炸事故造成嚴(yán)重后果影響最大的基本風(fēng)險因素為局部通風(fēng)方案設(shè)計不合理。

    (3)通過多因素耦合作用分析,預(yù)測掘進(jìn)工作面煤塵爆炸事故造成嚴(yán)重后果狀態(tài)概率變化情況為先急劇上升,后緩慢上升,最后趨于某一定值;在5種因素組合條件下耦合作用結(jié)果導(dǎo)致掘進(jìn)工作面事故造成嚴(yán)重后果的后驗概率最大,為0.266;掘進(jìn)工作面煤塵爆炸事故造成嚴(yán)重后果的預(yù)防重點是在事故發(fā)展初期對所選取的5種因素采取相關(guān)措施防止這5種基本風(fēng)險因素同時發(fā)生。

    猜你喜歡
    煤塵后驗通風(fēng)
    燒爐取暖要注意通風(fēng)
    溫度對煤塵水分蒸發(fā)特性的影響規(guī)律研究
    煤礦安全(2022年10期)2022-11-09 02:29:22
    一例育雛室通風(fēng)不良造成雞苗慢性死亡的調(diào)查報告
    煤塵云著火敏感性影響因素的實驗研究
    基于對偶理論的橢圓變分不等式的后驗誤差分析(英)
    貝葉斯統(tǒng)計中單參數(shù)后驗分布的精確計算方法
    一種基于最大后驗框架的聚類分析多基線干涉SAR高度重建算法
    淺圓倉智能通風(fēng)系統(tǒng)在冬季通風(fēng)中的應(yīng)用分析
    火力發(fā)電廠輸煤系統(tǒng)煤塵綜合治理探討
    基于貝葉斯后驗?zāi)P偷木植可鐖F發(fā)現(xiàn)
    可以在线观看的亚洲视频| 白带黄色成豆腐渣| 一区二区三区高清视频在线| 99热网站在线观看| 免费搜索国产男女视频| 极品教师在线视频| 大香蕉久久网| 日韩欧美精品v在线| 国产一区二区亚洲精品在线观看| 国产黄色视频一区二区在线观看 | 午夜激情福利司机影院| 精品久久久久久久久久免费视频| 亚洲欧美精品自产自拍| 亚洲熟妇中文字幕五十中出| 成人av在线播放网站| 男女那种视频在线观看| 亚洲一区二区三区色噜噜| www.色视频.com| 男女视频在线观看网站免费| 人妻久久中文字幕网| 免费av观看视频| 午夜福利高清视频| 91久久精品国产一区二区三区| 国产淫片久久久久久久久| av黄色大香蕉| 中文亚洲av片在线观看爽| 日本 av在线| 一夜夜www| 最新在线观看一区二区三区| 久久鲁丝午夜福利片| 自拍偷自拍亚洲精品老妇| 久久久久久大精品| 久久久久性生活片| 久久草成人影院| 搡老熟女国产l中国老女人| 97碰自拍视频| 露出奶头的视频| 狂野欧美白嫩少妇大欣赏| 天堂av国产一区二区熟女人妻| 国产女主播在线喷水免费视频网站 | 日韩精品青青久久久久久| 女生性感内裤真人,穿戴方法视频| 欧美色视频一区免费| 十八禁国产超污无遮挡网站| 亚洲中文字幕一区二区三区有码在线看| 久久久久久伊人网av| 久久天躁狠狠躁夜夜2o2o| 亚洲久久久久久中文字幕| 少妇熟女aⅴ在线视频| 久久人人精品亚洲av| 国产一区二区三区在线臀色熟女| 欧美色欧美亚洲另类二区| 少妇的逼水好多| 午夜福利高清视频| 欧美性猛交╳xxx乱大交人| 久久精品国产亚洲网站| 内射极品少妇av片p| 国产 一区精品| 我要搜黄色片| 精品99又大又爽又粗少妇毛片| 久久鲁丝午夜福利片| 3wmmmm亚洲av在线观看| 欧美xxxx黑人xx丫x性爽| 亚洲国产精品久久男人天堂| 高清午夜精品一区二区三区 | 超碰av人人做人人爽久久| 在线观看午夜福利视频| 99热这里只有是精品在线观看| 久久久久精品国产欧美久久久| 国产亚洲精品av在线| av天堂在线播放| 国产精品爽爽va在线观看网站| 中文字幕精品亚洲无线码一区| 午夜免费男女啪啪视频观看 | 亚洲不卡免费看| 网址你懂的国产日韩在线| 菩萨蛮人人尽说江南好唐韦庄 | 国产精品一及| 亚洲成人久久爱视频| 赤兔流量卡办理| 又粗又爽又猛毛片免费看| 亚洲无线在线观看| 极品教师在线视频| 亚洲人成网站在线播放欧美日韩| 精品午夜福利在线看| 欧美xxxx黑人xx丫x性爽| 亚洲国产色片| 色尼玛亚洲综合影院| 成人精品一区二区免费| 亚洲精品一卡2卡三卡4卡5卡| 变态另类成人亚洲欧美熟女| 男人舔奶头视频| 欧美性猛交╳xxx乱大交人| 女人被狂操c到高潮| 不卡视频在线观看欧美| 99久久精品热视频| 精品人妻熟女av久视频| 日韩,欧美,国产一区二区三区 | 国产乱人视频| 综合色丁香网| 国产av一区在线观看免费| 欧美日韩一区二区视频在线观看视频在线 | av天堂在线播放| 一级毛片久久久久久久久女| 校园春色视频在线观看| 亚洲欧美精品自产自拍| 偷拍熟女少妇极品色| 国产v大片淫在线免费观看| 三级经典国产精品| 亚洲av成人av| 成人二区视频| 女同久久另类99精品国产91| 丰满乱子伦码专区| 夜夜夜夜夜久久久久| 人妻制服诱惑在线中文字幕| 国产高清激情床上av| 欧美xxxx黑人xx丫x性爽| 日本一本二区三区精品| 国产精品福利在线免费观看| 97热精品久久久久久| 成人美女网站在线观看视频| 日本免费一区二区三区高清不卡| 久久6这里有精品| 精品久久久久久久久久久久久| 国产精品久久久久久久电影| 国产精品爽爽va在线观看网站| 干丝袜人妻中文字幕| 欧美激情久久久久久爽电影| 国产精品人妻久久久影院| 能在线免费观看的黄片| 国产精品一区二区三区四区免费观看 | 三级毛片av免费| 国产成人一区二区在线| 亚洲av.av天堂| av女优亚洲男人天堂| 亚洲欧美成人精品一区二区| 国产精品久久久久久精品电影| 午夜老司机福利剧场| 精品久久久噜噜| 波多野结衣高清无吗| 69人妻影院| 天堂影院成人在线观看| 91午夜精品亚洲一区二区三区| 综合色丁香网| 久久精品国产鲁丝片午夜精品| 不卡一级毛片| 日本a在线网址| 直男gayav资源| 寂寞人妻少妇视频99o| 亚洲精品国产成人久久av| 欧美3d第一页| 亚洲欧美精品自产自拍| 欧美性感艳星| 免费黄网站久久成人精品| 亚洲国产高清在线一区二区三| 麻豆国产97在线/欧美| 日本黄色片子视频| 免费观看在线日韩| 国产精品精品国产色婷婷| 成熟少妇高潮喷水视频| 欧美一区二区国产精品久久精品| 美女 人体艺术 gogo| 搞女人的毛片| 精品午夜福利在线看| 97超级碰碰碰精品色视频在线观看| 99久久九九国产精品国产免费| 久久中文看片网| 亚洲人与动物交配视频| 欧美日韩综合久久久久久| 亚洲国产欧美人成| 国产av一区在线观看免费| 成人毛片a级毛片在线播放| 免费大片18禁| 级片在线观看| 此物有八面人人有两片| 高清毛片免费观看视频网站| 亚洲性久久影院| 菩萨蛮人人尽说江南好唐韦庄 | 亚洲av美国av| 欧美中文日本在线观看视频| 你懂的网址亚洲精品在线观看 | 日韩欧美三级三区| 老司机午夜福利在线观看视频| 一级黄色大片毛片| 18禁在线播放成人免费| 成人综合一区亚洲| 久久99热这里只有精品18| 精品人妻熟女av久视频| 亚洲成人av在线免费| 久久精品综合一区二区三区| 亚洲av电影不卡..在线观看| 国产精品久久久久久亚洲av鲁大| 午夜老司机福利剧场| 中文字幕av在线有码专区| 日韩精品有码人妻一区| 麻豆精品久久久久久蜜桃| 在线观看66精品国产| 亚洲av成人av| 亚洲图色成人| 99国产极品粉嫩在线观看| 亚洲五月天丁香| 麻豆一二三区av精品| 一区福利在线观看| 亚洲人与动物交配视频| 亚洲电影在线观看av| 亚洲精品国产成人久久av| 亚洲不卡免费看| 国内少妇人妻偷人精品xxx网站| 成人午夜高清在线视频| 青春草视频在线免费观看| 成人三级黄色视频| 女生性感内裤真人,穿戴方法视频| 桃色一区二区三区在线观看| 国产成人精品久久久久久| 免费观看人在逋| 十八禁国产超污无遮挡网站| 欧美日韩在线观看h| 亚洲最大成人手机在线| 婷婷精品国产亚洲av| 长腿黑丝高跟| 精品人妻一区二区三区麻豆 | 给我免费播放毛片高清在线观看| 国产精品一区二区三区四区免费观看 | 成人亚洲精品av一区二区| 寂寞人妻少妇视频99o| 国产黄色小视频在线观看| 久久久久久久亚洲中文字幕| 婷婷色综合大香蕉| a级毛色黄片| 夜夜看夜夜爽夜夜摸| 国国产精品蜜臀av免费| 99热6这里只有精品| 国产亚洲精品av在线| 自拍偷自拍亚洲精品老妇| 欧美一区二区国产精品久久精品| 久久久国产成人免费| 人妻丰满熟妇av一区二区三区| 91av网一区二区| 1024手机看黄色片| 内射极品少妇av片p| 亚洲国产欧洲综合997久久,| 欧美高清性xxxxhd video| 天天躁夜夜躁狠狠久久av| 欧美成人a在线观看| 18+在线观看网站| 色吧在线观看| 中出人妻视频一区二区| 波多野结衣高清无吗| 高清毛片免费观看视频网站| 国产黄片美女视频| 日本爱情动作片www.在线观看 | 国产精品女同一区二区软件| 国内精品久久久久精免费| 色噜噜av男人的天堂激情| 可以在线观看的亚洲视频| 欧美最新免费一区二区三区| 亚洲七黄色美女视频| 国产欧美日韩精品一区二区| 亚洲一级一片aⅴ在线观看| 国产v大片淫在线免费观看| 一进一出抽搐gif免费好疼| 亚洲精品影视一区二区三区av| 人人妻人人看人人澡| 精品乱码久久久久久99久播| 特大巨黑吊av在线直播| 日韩精品青青久久久久久| 亚洲va在线va天堂va国产| 欧美日本视频| 久久久久久大精品| 99精品在免费线老司机午夜| 99热网站在线观看| 中文字幕免费在线视频6| 日韩av不卡免费在线播放| 99riav亚洲国产免费| 国产精品乱码一区二三区的特点| 我要看日韩黄色一级片| 麻豆一二三区av精品| av免费在线看不卡| 国产av一区在线观看免费| 亚洲成人精品中文字幕电影| 一本久久中文字幕| 国产成人a区在线观看| 寂寞人妻少妇视频99o| 国产精品日韩av在线免费观看| 亚洲国产精品成人久久小说 | 在线观看免费视频日本深夜| 国产高清有码在线观看视频| 国产午夜精品论理片| 亚洲国产欧美人成| 国产精品野战在线观看| 亚洲内射少妇av| 观看免费一级毛片| 国产精品免费一区二区三区在线| 国产在视频线在精品| 97人妻精品一区二区三区麻豆| 狂野欧美白嫩少妇大欣赏| 永久网站在线| 国内精品宾馆在线| 亚洲中文字幕一区二区三区有码在线看| 男女视频在线观看网站免费| 婷婷亚洲欧美| 成人午夜高清在线视频| 亚洲欧美日韩东京热| 噜噜噜噜噜久久久久久91| 国产色婷婷99| 国产亚洲91精品色在线| 观看免费一级毛片| 国产精品,欧美在线| 国产精品久久电影中文字幕| 亚洲成人精品中文字幕电影| 麻豆av噜噜一区二区三区| 国产乱人视频| 熟女人妻精品中文字幕| 12—13女人毛片做爰片一| 高清毛片免费观看视频网站| 熟女电影av网| 日韩欧美精品免费久久| 亚洲精品国产av成人精品 | av卡一久久| 两性午夜刺激爽爽歪歪视频在线观看| 听说在线观看完整版免费高清| 人妻丰满熟妇av一区二区三区| 亚洲av.av天堂| 免费黄网站久久成人精品| 午夜福利在线在线| 热99re8久久精品国产| 国产精品一区二区性色av| 高清毛片免费观看视频网站| 国产精品1区2区在线观看.| 岛国在线免费视频观看| 又黄又爽又免费观看的视频| 欧美zozozo另类| 久久精品国产亚洲av涩爱 | 成年免费大片在线观看| 国产精品福利在线免费观看| 在线看三级毛片| 亚洲欧美精品自产自拍| 亚洲精品国产av成人精品 | 国产真实乱freesex| 国产极品精品免费视频能看的| 久久久久国内视频| 插逼视频在线观看| 97热精品久久久久久| 人妻制服诱惑在线中文字幕| 亚洲av熟女| 国产精品亚洲美女久久久| 噜噜噜噜噜久久久久久91| 精品久久久久久久人妻蜜臀av| 99在线人妻在线中文字幕| 精品一区二区三区av网在线观看| 日韩三级伦理在线观看| 最近最新中文字幕大全电影3| 男女做爰动态图高潮gif福利片| 成人亚洲精品av一区二区| 久久国产乱子免费精品| 成人二区视频| 美女免费视频网站| 国产91av在线免费观看| 婷婷六月久久综合丁香| 亚洲图色成人| 久久久久久久久大av| 欧美国产日韩亚洲一区| 亚洲av.av天堂| av黄色大香蕉| 日韩三级伦理在线观看| 观看免费一级毛片| 国产高清视频在线观看网站| 99热这里只有是精品50| 亚洲精品色激情综合| 国产高清三级在线| 天堂√8在线中文| 国产色婷婷99| 少妇被粗大猛烈的视频| 最近手机中文字幕大全| av在线亚洲专区| 久久精品国产清高在天天线| 午夜精品在线福利| 国产激情偷乱视频一区二区| 免费av毛片视频| 精品人妻一区二区三区麻豆 | 日本黄色片子视频| 美女免费视频网站| 精品日产1卡2卡| 久久精品国产自在天天线| 人妻久久中文字幕网| 精品久久久久久久久亚洲| 午夜福利成人在线免费观看| 国产精品av视频在线免费观看| 欧洲精品卡2卡3卡4卡5卡区| 色综合亚洲欧美另类图片| 尤物成人国产欧美一区二区三区| 亚州av有码| 精品久久久噜噜| 久久人人爽人人爽人人片va| 成人精品一区二区免费| 深爱激情五月婷婷| 特级一级黄色大片| 免费搜索国产男女视频| 高清午夜精品一区二区三区 | 联通29元200g的流量卡| 91久久精品国产一区二区三区| 在线国产一区二区在线| 国产精品精品国产色婷婷| 国产伦一二天堂av在线观看| 岛国在线免费视频观看| 波野结衣二区三区在线| 午夜福利视频1000在线观看| 日本黄大片高清| 国产真实乱freesex| 午夜福利高清视频| 99热全是精品| 婷婷色综合大香蕉| 国产精品久久电影中文字幕| 国产私拍福利视频在线观看| 中文字幕久久专区| 露出奶头的视频| 国产精品亚洲美女久久久| 天美传媒精品一区二区| 网址你懂的国产日韩在线| 最后的刺客免费高清国语| 免费观看在线日韩| 国产黄色视频一区二区在线观看 | 久久精品91蜜桃| 美女cb高潮喷水在线观看| 亚洲中文字幕一区二区三区有码在线看| 男人和女人高潮做爰伦理| 老熟妇仑乱视频hdxx| 一级黄片播放器| 欧美日韩在线观看h| 亚洲av电影不卡..在线观看| 免费电影在线观看免费观看| 黄色欧美视频在线观看| 国产精品爽爽va在线观看网站| 色5月婷婷丁香| 嫩草影院新地址| 国产黄片美女视频| 性色avwww在线观看| 国产极品精品免费视频能看的| 国产人妻一区二区三区在| 99国产极品粉嫩在线观看| 国产v大片淫在线免费观看| 日韩制服骚丝袜av| 亚洲欧美日韩高清在线视频| 听说在线观看完整版免费高清| 国内精品一区二区在线观看| 男人的好看免费观看在线视频| 精品无人区乱码1区二区| 欧美又色又爽又黄视频| 色av中文字幕| 日韩三级伦理在线观看| 97超碰精品成人国产| 69人妻影院| 欧美日韩综合久久久久久| 国产欧美日韩一区二区精品| 人人妻,人人澡人人爽秒播| 国产v大片淫在线免费观看| 国产成人91sexporn| 99热精品在线国产| 午夜亚洲福利在线播放| 国产熟女欧美一区二区| 少妇人妻精品综合一区二区 | 国产精品国产三级国产av玫瑰| 成人一区二区视频在线观看| 国产精品女同一区二区软件| 在线观看66精品国产| 精品一区二区三区视频在线| 国产午夜精品久久久久久一区二区三区 | 成人高潮视频无遮挡免费网站| 一边摸一边抽搐一进一小说| 日韩人妻高清精品专区| 人妻久久中文字幕网| 日韩欧美三级三区| 欧美日本视频| 欧美xxxx黑人xx丫x性爽| 成人av一区二区三区在线看| 亚洲一区高清亚洲精品| 我的女老师完整版在线观看| 看黄色毛片网站| 亚洲av熟女| 免费观看在线日韩| a级毛色黄片| 国产精品av视频在线免费观看| 日韩大尺度精品在线看网址| 日韩欧美免费精品| 国产亚洲精品久久久久久毛片| 精品人妻熟女av久视频| 国产一区二区在线观看日韩| 人妻丰满熟妇av一区二区三区| 国产精品精品国产色婷婷| 亚洲av二区三区四区| 亚洲国产色片| eeuss影院久久| 精品一区二区三区av网在线观看| 欧美一区二区亚洲| 熟女电影av网| 欧美色欧美亚洲另类二区| 久久这里只有精品中国| 真人做人爱边吃奶动态| 夜夜爽天天搞| 日韩av不卡免费在线播放| 午夜福利成人在线免费观看| 哪里可以看免费的av片| 人妻制服诱惑在线中文字幕| 亚洲av免费高清在线观看| 亚洲中文字幕一区二区三区有码在线看| 美女cb高潮喷水在线观看| 成人毛片a级毛片在线播放| 国产v大片淫在线免费观看| 国产一区二区在线观看日韩| 国产精品久久电影中文字幕| 免费电影在线观看免费观看| 国产亚洲精品av在线| 又爽又黄无遮挡网站| 我的女老师完整版在线观看| 中国美白少妇内射xxxbb| 97热精品久久久久久| 亚洲av中文av极速乱| 日本在线视频免费播放| 亚洲无线在线观看| 一级黄色大片毛片| 国产日本99.免费观看| www.色视频.com| 中文字幕av在线有码专区| 性插视频无遮挡在线免费观看| 91久久精品电影网| 精品久久久久久成人av| 亚洲国产欧美人成| 亚洲av二区三区四区| 亚洲人与动物交配视频| 久久婷婷人人爽人人干人人爱| 精品久久久久久久人妻蜜臀av| 免费搜索国产男女视频| 91午夜精品亚洲一区二区三区| 蜜桃久久精品国产亚洲av| 国产一区二区激情短视频| 精品福利观看| 一级黄片播放器| 久久久久久伊人网av| 亚洲欧美成人精品一区二区| 一级毛片我不卡| 自拍偷自拍亚洲精品老妇| 老女人水多毛片| 嫩草影院精品99| .国产精品久久| 韩国av在线不卡| 国产精品一及| 欧美+日韩+精品| 在现免费观看毛片| 国产一区二区在线观看日韩| 日韩一区二区视频免费看| 毛片一级片免费看久久久久| 此物有八面人人有两片| 日韩欧美国产在线观看| 久久综合国产亚洲精品| 国产一级毛片七仙女欲春2| 老司机午夜福利在线观看视频| 十八禁网站免费在线| 国产精品一区二区三区四区免费观看 | 亚洲专区国产一区二区| 热99re8久久精品国产| 99热网站在线观看| 亚洲性久久影院| 日韩亚洲欧美综合| 午夜a级毛片| 人人妻人人看人人澡| 18禁裸乳无遮挡免费网站照片| 村上凉子中文字幕在线| 日韩大尺度精品在线看网址| 色播亚洲综合网| 国产精品久久久久久亚洲av鲁大| 老女人水多毛片| 一夜夜www| 国内少妇人妻偷人精品xxx网站| 直男gayav资源| 偷拍熟女少妇极品色| 国产精品一区二区性色av| 欧美另类亚洲清纯唯美| 久久久久久大精品| 欧美日韩一区二区视频在线观看视频在线 | 欧美色视频一区免费| 午夜a级毛片| 国产亚洲精品综合一区在线观看| 春色校园在线视频观看| www日本黄色视频网| 午夜视频国产福利| 亚洲精品国产av成人精品 | 亚洲欧美日韩高清在线视频| 色播亚洲综合网| 亚洲精品在线观看二区| 久久久a久久爽久久v久久| 日日干狠狠操夜夜爽| 亚洲激情五月婷婷啪啪| 真实男女啪啪啪动态图| 99国产精品一区二区蜜桃av| 久久久午夜欧美精品| 精品久久久久久久久久免费视频| 日韩av在线大香蕉| 亚洲激情五月婷婷啪啪| 久久人妻av系列| 日日摸夜夜添夜夜爱| 人人妻,人人澡人人爽秒播| 久久久精品大字幕| 99久久无色码亚洲精品果冻| av在线观看视频网站免费| 国产色婷婷99| 精品久久久久久久末码| 深夜a级毛片| 国产在线男女| 久久国内精品自在自线图片| 少妇丰满av| av在线亚洲专区| 午夜爱爱视频在线播放| 高清日韩中文字幕在线| 亚洲欧美日韩东京热| 一级黄片播放器|