牛羊無(wú)布魯氏菌病區(qū)建設(shè) 第1部分：傳入風(fēng)險(xiǎn)定量評(píng)估

III牛羊無(wú)布魯氏菌病區(qū)建設(shè)第1部分：傳入風(fēng)險(xiǎn)定量評(píng)估1范圍本文件規(guī)定了牛羊布魯氏菌病（以下簡(jiǎn)稱“布病”）傳入風(fēng)險(xiǎn)定量評(píng)估的評(píng)估步驟。本文件適用于省內(nèi)引進(jìn)牛羊傳入布病的風(fēng)險(xiǎn)評(píng)估工作。2規(guī)范性文件引用 下列文件對(duì)于本文件的應(yīng)用是必不可少的。凡是注日期的引用文件，僅所注日期的版本適用于本文件。凡是不注日期的引用文件，其最新版本（包括所有的修改單）適用于本文件。GB/T18646動(dòng)物布魯氏菌病診斷技術(shù)GB/T27921風(fēng)險(xiǎn)管理風(fēng)險(xiǎn)評(píng)估技術(shù)世界動(dòng)物衛(wèi)生組織陸生動(dòng)物衛(wèi)生法典國(guó)際動(dòng)物衛(wèi)生法典3術(shù)語(yǔ)定義下列術(shù)語(yǔ)和定義適用于本文件。3.1風(fēng)險(xiǎn)因子risk可能影響疫病發(fā)生的不確定性因素。3.2傳入評(píng)估entryassessment指動(dòng)物疫病在特定條件下傳入特定區(qū)域的途徑和可能性。3.3定量風(fēng)險(xiǎn)評(píng)估quantitativeriskassessment定量研究是在定性研究基礎(chǔ)上，確定某動(dòng)物疫病發(fā)生的關(guān)鍵風(fēng)險(xiǎn)路徑，運(yùn)用情景樹(shù)法結(jié)合數(shù)理統(tǒng)計(jì)方法，逐步確定各關(guān)鍵風(fēng)險(xiǎn)點(diǎn)發(fā)生風(fēng)險(xiǎn)概率，提出特定情景下該疫病通過(guò)該風(fēng)險(xiǎn)路徑發(fā)生的風(fēng)險(xiǎn)概率，最終根據(jù)風(fēng)險(xiǎn)概率提出風(fēng)險(xiǎn)管理措施及建議。3.4風(fēng)險(xiǎn)管理riskmanagement在風(fēng)險(xiǎn)評(píng)估的基礎(chǔ)上，選擇、確定和實(shí)施能夠降低風(fēng)險(xiǎn)水平的措施及對(duì)措施開(kāi)展評(píng)價(jià)的過(guò)程。3.5養(yǎng)殖單元breedingunit把一些特征相似的動(dòng)物劃為同一個(gè)群，該群所處的地域定義為一個(gè)養(yǎng)殖單元。4評(píng)估步驟4.1數(shù)據(jù)來(lái)源與農(nóng)業(yè)農(nóng)村部門(mén)、動(dòng)物疫病預(yù)防控制機(jī)構(gòu)、動(dòng)物衛(wèi)生監(jiān)督機(jī)構(gòu)、養(yǎng)殖場(chǎng)（戶）、活牛羊育肥場(chǎng)場(chǎng)主、活牛羊販運(yùn)人和地方基層獸醫(yī)等有關(guān)人員開(kāi)展座談交流，并查閱相關(guān)實(shí)驗(yàn)室檢測(cè)數(shù)據(jù)，獲取周邊地區(qū)活牛羊調(diào)入途徑、路線、數(shù)量等數(shù)據(jù)。4.2假設(shè)條件成立假設(shè)省內(nèi)引進(jìn)牛羊（群）均為布病易感動(dòng)物，在裝車(chē)、運(yùn)輸?shù)拳h(huán)節(jié)中假定不存在被感染的情況。研究中的檢測(cè)方法為生產(chǎn)實(shí)踐中最常用的兩種垂直血清學(xué)檢測(cè)方法，即虎紅平板凝集試驗(yàn)與試管凝集試驗(yàn)垂直檢測(cè)（敏感性81%、特異性98%）、虎紅平板凝集試驗(yàn)與cELISA垂直檢測(cè)方法（敏感性89%、特異性97%）。試驗(yàn)步驟依照標(biāo)準(zhǔn)GB/T18646-2018執(zhí)行，假定兩種檢測(cè)方法隨機(jī)使用。4.3傳入風(fēng)險(xiǎn)路徑模型本文件針對(duì)省內(nèi)引進(jìn)牛羊傳入布病的可能性開(kāi)展相關(guān)風(fēng)險(xiǎn)評(píng)估。以引進(jìn)羊?yàn)槔?，根?jù)實(shí)際羊只引進(jìn)流程，繪制了布病傳播路徑“情景樹(shù)”模型（圖1）。圖1布病傳播路徑“情景樹(shù)”模型4.4模型參數(shù)取值（1）檢測(cè)試驗(yàn)敏感性（Se）與特異性（Sp）已知檢測(cè)方法1中，虎紅平板凝集試驗(yàn)與試管凝集試驗(yàn)為垂直檢測(cè)，其敏感性為81%，特異性為98%；檢測(cè)方法2中，虎紅平板凝集試驗(yàn)與cELISA為垂直檢測(cè)，其敏感性為89%，特異性為97%。兩種檢測(cè)方法在實(shí)際工作中均普遍應(yīng)用，所以在模型中應(yīng)用均勻分布模擬試驗(yàn)的敏感性（Se）和特異性（Sp）（2）場(chǎng)群流行率（HP）、個(gè)體流行率（P1）、產(chǎn)地檢疫發(fā)現(xiàn)陽(yáng)性動(dòng)物概率（P2）本文件均采用Beta分布模擬，即HP&P1&P2=Beta（M+1，N-M+1），其中M取值為檢測(cè)或排查出陽(yáng)性數(shù)量，N為總樣品數(shù)量或排查總數(shù)量。（3）入場(chǎng)前隔離飼養(yǎng)發(fā)現(xiàn)陽(yáng)性動(dòng)物概率（P3）通過(guò)走訪問(wèn)詢、查閱實(shí)驗(yàn)室相關(guān)檢測(cè)記錄、相關(guān)文獻(xiàn)、向羊（牛）養(yǎng)殖場(chǎng)（戶）發(fā)放調(diào)查問(wèn)卷，綜合分析得出假設(shè)引進(jìn)100只羊（牛），通過(guò)產(chǎn)地檢疫、隔離飼養(yǎng)方式能發(fā)現(xiàn)幾只患病羊（牛），最多能發(fā)現(xiàn)幾只以及最少能發(fā)現(xiàn)幾只患病羊（牛）等問(wèn)題。分析調(diào)查結(jié)果，運(yùn)用Pert分布模擬隔離飼養(yǎng)排查出陽(yáng)性動(dòng)物概率（P3），即P3=Pert（Min，MostL，Max），其中Min為最小取值，MostL為最可能取值，Max為最大取值。4.5模型參數(shù)描述參見(jiàn)附錄A。4.6模型概率計(jì)算隨機(jī)引進(jìn)1只羊（牛）經(jīng)產(chǎn)地檢疫、隔離飼養(yǎng)、混群前檢測(cè)篩選陽(yáng)性動(dòng)物后，傳入布病的概率為P（D+|T-）。按照引進(jìn)羊（牛）只傳入風(fēng)險(xiǎn)路徑模型，羊（牛）場(chǎng)群流行率（HP）、場(chǎng)內(nèi)個(gè)體流行率（P1）、產(chǎn)地檢疫發(fā)現(xiàn)陽(yáng)性動(dòng)物概率（P2）、隔離飼養(yǎng)發(fā)現(xiàn)陽(yáng)性動(dòng)物概率（P3）以及兩次診斷試驗(yàn)的敏感性（Se）與特異性（Sp）都是影響隨機(jī)引進(jìn)1只羊（牛）經(jīng)產(chǎn)地檢疫、混群前檢測(cè)后傳入布病概率的主要影響因素。P（D+|T-）=4.7風(fēng)險(xiǎn)計(jì)算通過(guò)上述風(fēng)險(xiǎn)路徑“情景樹(shù)”模型，將模型參數(shù)的取值及估計(jì)值帶入公式（1），通過(guò)@Risk8.0風(fēng)險(xiǎn)分析軟件，采用蒙特卡洛仿真模擬方法，對(duì)模型迭代10000次仿真分析，計(jì)算假定陜西省內(nèi)某羊（牛）場(chǎng)引進(jìn)100只羊（牛），至少有1只為感染羊（牛）而導(dǎo)致羊（牛）群發(fā)病的可能性為α。α（x≥1）=1-[1-P（D+|T-）]100（1）4.8敏感性分析采用@Risk8.0軟件進(jìn)行模型參數(shù)敏感性分析，根據(jù)敏感性分析結(jié)果確定主要的保護(hù)性因素和風(fēng)險(xiǎn)因素有哪些。4.9評(píng)估報(bào)告（1）評(píng)估目的評(píng)估省內(nèi)引進(jìn)牛羊傳入布病的風(fēng)險(xiǎn)大小，從而為陜西省布病無(wú)疫區(qū)建設(shè)提供支撐，為建立動(dòng)物及產(chǎn)品風(fēng)險(xiǎn)預(yù)警系統(tǒng)提供基礎(chǔ)。（2）評(píng)估組織評(píng)估單位應(yīng)成立3-5人風(fēng)險(xiǎn)評(píng)估小組，成員應(yīng)包括獸醫(yī)、公共衛(wèi)生、流行病學(xué)調(diào)查等領(lǐng)域?qū)I(yè)人才。風(fēng)險(xiǎn)評(píng)估小組人員名單如下：分類分類表2風(fēng)險(xiǎn)評(píng)估小組組別組別（3）評(píng)估方法省內(nèi)引進(jìn)牛羊傳入布病風(fēng)險(xiǎn)采取定量分析方法，通過(guò)繪制傳入風(fēng)險(xiǎn)路徑“情景樹(shù)”模型，將模型參數(shù)的取值及估計(jì)值帶入公式，通過(guò)@Risk8.0風(fēng)險(xiǎn)分析軟件，采用蒙特卡洛仿真模擬方法，對(duì)模型迭代10000次仿真分析，計(jì)算假定在省內(nèi)某羊（牛）場(chǎng)引進(jìn)100只羊（牛），至少有1只為感染羊（牛）而導(dǎo)致羊（牛）群發(fā)病的可能性為α。（4）評(píng)估結(jié)論和建議經(jīng)動(dòng)物疫病風(fēng)險(xiǎn)評(píng)估，并根據(jù)評(píng)估支撐信息的完整性、方法的科學(xué)性等，對(duì)存在的不確定性進(jìn)行分析，形成動(dòng)物疫病風(fēng)險(xiǎn)評(píng)估報(bào)告，提出風(fēng)險(xiǎn)管理措施建議。風(fēng)險(xiǎn)管理措施建議按下列歸類：1）限制特定區(qū)域（場(chǎng)所）動(dòng)物及動(dòng)物產(chǎn)品調(diào)運(yùn)；2）加強(qiáng)特定動(dòng)物疫病免疫、監(jiān)測(cè)、流行病學(xué)調(diào)查、檢疫監(jiān)管應(yīng)急處置、疫病凈化等綜合防控措施；3）改善基礎(chǔ)動(dòng)物防疫條件，提高生物安全防控水平；4）加強(qiáng)人畜共患傳染病人畜同步協(xié)同防控；5）完善動(dòng)物防疫體系，確保與承擔(dān)職能任務(wù)相匹配；6）需要落實(shí)風(fēng)險(xiǎn)閉環(huán)管控的其他建議。附錄A（資料性附錄）表A模型參數(shù)取值描述參數(shù)描述分布參數(shù)值群流行率（HP）%抽樣場(chǎng)戶中至少有1只布病陽(yáng)性動(dòng)物的規(guī)模場(chǎng)所占比例Beta個(gè)體流行率（P1）%抽樣個(gè)體中檢測(cè)結(jié)果為感染抗體陽(yáng)性動(dòng)物數(shù)所占比例Beta產(chǎn)地檢疫發(fā)現(xiàn)陽(yáng)性