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

    小型自然循環(huán)鉛鉍快堆無保護(hù)超功率事故不確定性分析研究

    2022-08-17 02:19:18余清遠(yuǎn)漆靜雯趙鵬程趙亞楠
    核技術(shù) 2022年8期
    關(guān)鍵詞:分析

    余清遠(yuǎn) 漆靜雯 趙鵬程 趙亞楠 于 濤

    (南華大學(xué)核科學(xué)技術(shù)學(xué)院 衡陽421001)

    采用液態(tài)鉛鉍冷卻的鉛鉍堆在嬗變核廢料和增殖核燃料方面具有獨(dú)特優(yōu)勢,擁有良好的經(jīng)濟(jì)性和固有安全性,被第四代國際核能論壇預(yù)計(jì)有望成為首個(gè)實(shí)現(xiàn)工業(yè)示范化的第四代核能系統(tǒng)[1]。事故瞬態(tài)是鉛鉍快堆熱工安全分析的重要內(nèi)容,有必要探究鉛鉍堆瞬態(tài)過程中關(guān)鍵參數(shù)的變化規(guī)律,評估鉛鉍堆的固有安全性。由于設(shè)計(jì)參數(shù)、初始邊界條件等因素的影響,鉛鉍堆瞬態(tài)模擬計(jì)算結(jié)果不可避免地存在不確定性,將不確定性分析方法應(yīng)用鉛鉍堆瞬態(tài)事故分析,可量化不確定性分析結(jié)果,評估瞬態(tài)工況中輸入?yún)?shù)對鉛鉍堆安全的影響,對鉛鉍堆的設(shè)計(jì)分析、安全審評具有現(xiàn)實(shí)意義。

    國際上,國際性合作研究項(xiàng)目OECD(Organization for Economic Co-operation and Development)/NEA LWR UAM(Nuclear Energy Agency,Light Water Reactor,Uncertainty Analysis in Modeling)正在全力開展有關(guān)多物理場、多尺度耦合的不確定性分析[2-5];日本基于京都大學(xué)臨界裝置開展了一系列加速器驅(qū)動系統(tǒng)實(shí)驗(yàn),并采用隨機(jī)采樣技術(shù)對實(shí)驗(yàn)開展了不確定性量化分析[6-9];瑞典皇家理工學(xué)院基于ROAAM+、MELCOR 等程序?qū)Ρ睔W沸水堆開展了不確定性量化和敏感性分析[10-13];由美國阿貢實(shí)驗(yàn)室牽頭的先進(jìn)反應(yīng)堆E級高性能模擬中心開展了堆芯物理計(jì)算的不確定量化分析和敏感性分析[14-15]。近年來,國內(nèi)有關(guān)快堆的不確定性研究主要集中于不確定性分析方法:西安交通大學(xué)先后采用不確定性分析方法開展了鈉冷快堆事故SHRT-17 試驗(yàn)數(shù)據(jù)和BN-600 快堆有效增值因數(shù)的不確定性和敏感性分析[16-17];哈爾濱工程大學(xué)基于抽樣統(tǒng)計(jì)法開展了中國實(shí)驗(yàn)快堆(China Experimental Fast Reactor,CEFR)的不確定性分析和敏感性分析[18];中國原子能研究院采用多種不確定性分析方法開展了CEFR堆芯設(shè)計(jì)計(jì)算的不確定度定量評價(jià)[19]。中國科學(xué)技術(shù)大學(xué)針對鉛鉍快堆開發(fā)了多物理多尺度程序Fluent/KMC-sub/NDK[20],并對SNCLFR-100開展了啟堆條件下的流動不穩(wěn)定性分析[21]。而鉛鉍快堆作為具有發(fā)展?jié)摿Φ暮蜻x堆型之一,其瞬態(tài)安全特性對反應(yīng)堆運(yùn)行至關(guān)重要,開展鉛鉍快堆的瞬態(tài)事故的不確定性量化分析,研究瞬態(tài)事故的不確定性來源,對鉛鉍快堆的設(shè)計(jì)、安全特性評估有重要意義。

    本工作是針對中國科學(xué)技術(shù)大學(xué)和南華大學(xué)聯(lián)合設(shè)計(jì)的小型自然循環(huán)鉛基鉛鉍快堆(Small Natural Circulation Lead-bismuth Cooled Fast Reactor 10 MWth,SNCLFR-10),首先,通過文獻(xiàn)資料以及工程經(jīng)驗(yàn)判斷,考慮三類不確定性輸入?yún)?shù):熱工水力模型參數(shù)、中子物理模型參數(shù)以及燃料制造公差;然后,基于課題組自主開發(fā)的物理/熱工耦合不確定性量化分析程序CFD/PFS 程序開展無保護(hù)超功率事故的不確定性分析,采用拉丁超立方(Latin Hypercube Sampling,LHS)對輸入?yún)?shù)進(jìn)行抽樣,對瞬態(tài)安全參數(shù)進(jìn)行不確定性量化分析,并通過相關(guān)系數(shù)分析法評估輸入?yún)?shù)對瞬態(tài)安全參數(shù)的敏感度。

    1 計(jì)算模型

    1.1 分析模型

    SNCLFR-10 堆是由中國科學(xué)技術(shù)大學(xué)和南華大學(xué)聯(lián)合設(shè)計(jì)的小型自然循環(huán)鉛鉍快堆,具有良好的現(xiàn)實(shí)可行性、安全可靠性、實(shí)驗(yàn)靈活性和技術(shù)延續(xù)性[22]。圖1 為SNCLFR-10 堆的結(jié)構(gòu)設(shè)計(jì)圖,表1、2給出了反應(yīng)堆設(shè)計(jì)參數(shù)及堆芯主要設(shè)計(jì)參數(shù)。SNCLFR-10堆的熱功率為10 MW,燃料采用富集度為19.75%的UO2,燃料包殼材料采用316Ti不銹鋼。堆芯活性區(qū)由94盒組件組成,其中86盒燃料組件,8盒控制棒組件。一回路冷卻劑系統(tǒng)采用自然循環(huán)驅(qū)動液態(tài)鉛鉍冷卻,冷卻劑入口溫度為533.15 K,出口溫度為663.15 K;二回路冷卻劑采用4 MPa 的加壓液態(tài)水。

    表1 SNCLFR-10主要設(shè)計(jì)參數(shù)Table 1 Main design parameters of SNCLFR-10

    圖1 SNCLFR-10堆結(jié)構(gòu)示意圖Fig.1 Structure diagram of SNCLFR-10

    無保護(hù)超功率事故(Unprotected Transient Overpower,UTOP)是指向反應(yīng)堆內(nèi)突然引入一個(gè)正反應(yīng)性,導(dǎo)致功率急劇升高的瞬態(tài)過程,此過程中反應(yīng)堆停堆保護(hù)系統(tǒng)失效。UTOP 事故的瞬態(tài)條件:反應(yīng)堆初始工況為穩(wěn)態(tài)運(yùn)行,假設(shè)最大價(jià)值控制棒失去控制而不可抽出,故在2 s 內(nèi)迅速引入1$正反應(yīng)性,整個(gè)過程無停堆保護(hù)動作。

    1.2 CFD/PFS程序

    CFD/PFS 是基于CFD 程序FLUENT 的用戶自定義函數(shù)(User Defined Function,UDF)耦合點(diǎn)堆中子動力學(xué)模型(Point Kinetics Model,PKM)、燃料棒熱傳導(dǎo)模型(Fuel Pin Heat Transfer Model,F(xiàn)PH)以及不確定性分析程序SIMLAB開發(fā)的一套物理熱工耦合不確定性分析程序,可用于分析小型鉛鉍快堆的瞬態(tài)事故。SIMLAB 是研究不確定性的開源程序,由歐洲聯(lián)合研究中心主持研究開發(fā),配備可視化操作界面,其中含有統(tǒng)計(jì)預(yù)處理子系統(tǒng)、模型執(zhí)行子系統(tǒng)、統(tǒng)計(jì)后處理子系統(tǒng),廣泛用于不確定性和整體敏感性研究[23]。采用與SIMMER 程序code-to-code的方法進(jìn)行CFD/PFS 的對比驗(yàn)證,CFD/PFS 程序開發(fā)與詳細(xì)驗(yàn)證見參考文獻(xiàn)[24]。

    CFD/PFS程序分為物理熱工耦合計(jì)算模型與不確定性分析計(jì)算,其計(jì)算流程如圖2 所示。物理熱工耦合計(jì)算模型包括CFD模型、點(diǎn)堆中子動力學(xué)模型、燃料棒熱傳導(dǎo)模型。PKM是三維時(shí)空動力學(xué)的簡化模型,采用6 組緩發(fā)中子的點(diǎn)堆動力學(xué)方程主要用于計(jì)算反應(yīng)性變化:

    圖2 CFD/PFS數(shù)據(jù)傳遞流程Fig.2 The data transfer route of CFD/PFS

    式中:n(t)是中子密度;ρ(t)是反應(yīng)性;β是緩發(fā)中子總份額;βi是i組緩發(fā)中子的份額;Λ是中子代時(shí)間;λi是第i組緩發(fā)中子的衰變常數(shù);Ci是i組緩發(fā)中子先驅(qū)核的濃度。

    再根據(jù)反應(yīng)性通過點(diǎn)堆動力學(xué)方程求解相對中子密度n(t),則總功率為:

    反應(yīng)性反饋源于反應(yīng)堆溫度、壓力或者流量的變化,而溫度對反應(yīng)性的影響是主要反饋,故所有反應(yīng)性反饋、堆內(nèi)溫度總反應(yīng)性公式可表示為:

    式中:KD為多普勒系數(shù);αc,a為堆芯軸向膨脹反饋系數(shù);αc,γ為堆芯徑向膨脹反饋系數(shù);αc為冷卻劑溫度反饋系數(shù);Tf(t)和Tf(0)分別為t時(shí)刻的燃料、冷卻劑平均溫度。

    FTM包括燃料棒穩(wěn)態(tài)和瞬態(tài)熱傳導(dǎo)模型,主要用于計(jì)算燃料棒的溫度分布,包括燃料芯塊、氣隙以及包殼的溫度,提供熱流密度給CFD 模型,反饋芯塊平均溫度給PKM;CFD模型主要用于計(jì)算堆芯熱工水力參數(shù),包括冷卻劑的溫度、速度等,反饋冷卻劑平均溫度給PKM,提供冷卻劑溫度和速度給FTM。不確定性分析計(jì)算流程如下:1)確定重要輸入?yún)?shù)及其分布范圍;2)利用SimLab程序?qū)Σ淮_定性輸入?yún)?shù)抽樣,生成輸入工況;3)將輸入工況傳遞給物理熱工耦合計(jì)算模型進(jìn)行計(jì)算,生成輸出工況;4)將輸出工況傳遞給SimLab程序,對結(jié)果進(jìn)行不確 定性分析和敏感性分析。

    表2 SNCLFR-10堆芯主要設(shè)計(jì)參數(shù)Table 2 Core design parameters of SNCLFR-10

    2 不確定性方法及關(guān)鍵參數(shù)識別

    2.1 不確定性分析方法

    本工作基于拉丁超立方抽樣的非參數(shù)統(tǒng)計(jì)法進(jìn)行不確定性分析,拉丁超立方采取分層抽樣,LHS抽樣過程如下:1)根據(jù)輸入?yún)?shù)的數(shù)量及分布范圍等概率的分成n個(gè)區(qū)域;2)在每個(gè)區(qū)域內(nèi)進(jìn)行隨機(jī)抽樣;3)將抽樣產(chǎn)生的樣本進(jìn)行隨機(jī)組合;4)挑選滿足約束條件的組合樣本作為輸入?yún)?shù)樣本。與隨機(jī)抽樣的蒙特卡羅方法比較,拉丁超立方抽樣可用較少樣本完成不確定性量化,且在均值和方差的估計(jì)上有顯著改善。

    非參數(shù)統(tǒng)計(jì)法(wilks 公式)根據(jù)分布區(qū)間的置信度β以及輸出參數(shù)的概率γ確定所需最小抽樣數(shù)目n。wilks公式如下:

    其中,式(4)適用于單側(cè)限值分布,式(5)適用于雙側(cè)限值分布。

    由wilks 公式計(jì)算出的滿足不同置信度和容許區(qū)間的最小樣本量n值如表3 所示。根據(jù)非參數(shù)統(tǒng)計(jì)理論,最小樣本量選取59 時(shí),目標(biāo)參數(shù)的容忍上限值取最大值;最小樣本量選取93 時(shí),目標(biāo)參數(shù)的容忍上限值取次大值[25]。

    表3 不同置信度和容許區(qū)間下最小樣本信息分布Table 3 The minimum number of samples under different confidence and tolerance intervals

    2.2 關(guān)鍵參數(shù)的選取

    反應(yīng)堆瞬態(tài)計(jì)算中的不確定性輸入?yún)?shù)主要考慮三部分:1)物理模型輸入?yún)?shù);2)熱工水力模型輸入?yún)?shù);3)燃料制造公差輸入?yún)?shù)。針對SNCLFR-10堆的UTOP事故,根據(jù)始發(fā)事件選取清單、文獻(xiàn)資料以及工程經(jīng)驗(yàn)判斷等方式選取不確定性輸入?yún)?shù),并確定其范圍和概率密度分布。選取的不確定輸入?yún)?shù)分為三個(gè)部分:首先是幾何輸入?yún)?shù),一些參數(shù)測量有一定的不安全性故缺乏,只能由專家假設(shè)估算;其次是中子學(xué)參數(shù),因?yàn)閷?shí)驗(yàn)測量系統(tǒng)及統(tǒng)計(jì)誤差,故其引入的不確定性需考慮;最后是熱工水力參數(shù),程序、模型及經(jīng)驗(yàn)關(guān)系式等都會引入不確定性,故選取其中最典型5 個(gè)參數(shù)。表4 列出了25 個(gè)不確定輸入?yún)?shù)信息,輸入?yún)?shù)服從均勻分布或正態(tài)分布,其中不確定性度表示參數(shù)名義值(設(shè)計(jì)值)和參數(shù)上下限值的相對偏差。

    表4 不確定性輸入?yún)?shù)信息Table 4 Uncertainty distribution of input parameter

    由于測量儀器精準(zhǔn)、技術(shù)先進(jìn),反應(yīng)堆部件制造精度高,因此燃料制造公差輸入?yún)?shù)的不確定度設(shè)為±0.5%。中子代時(shí)間、緩發(fā)中子份額、衰變常數(shù)以及初始功率的不確定度參考其他類似反應(yīng)堆的不確定度;對于反應(yīng)性反饋系數(shù),考慮它包含了一系列不確定性(材料屬性、幾何形狀、截面數(shù)據(jù)以及蒙特卡羅計(jì)算),為保證其保守性,將反應(yīng)性反饋系數(shù)的不確定度設(shè)為10%[26]。熱工水力現(xiàn)象極其復(fù)雜,熱工物性參數(shù)大多是根據(jù)實(shí)驗(yàn)數(shù)據(jù)擬合以及經(jīng)驗(yàn)關(guān)系式所得,因此熱工水力模型輸入?yún)?shù)通常存在較大的不確定度[27]。

    2.3 關(guān)鍵參數(shù)的抽樣

    基于SIMLAB程序采用拉丁超立方方法實(shí)現(xiàn)不確定性輸入?yún)?shù)抽樣,種子數(shù)設(shè)定為98 765。根據(jù)wilks 非參數(shù)統(tǒng)計(jì)法,在95%置信水平、95%概率水平下,所需最小計(jì)算次數(shù)為93次??紤]到抽樣可能存在壞點(diǎn)、程序可能運(yùn)行失效等原因,抽樣數(shù)目設(shè)定為100 組,其中每一組工況包括25 個(gè)輸入?yún)?shù)。為直觀地觀測輸入?yún)?shù)的概率分布特征、檢測其抽樣的合理性,對不確定性輸入?yún)?shù)進(jìn)行歸一化處理。圖3為不確定性輸入?yún)?shù)的拉丁超立方抽樣樣本歸一化結(jié)果。

    圖3 不確定性輸入?yún)?shù)抽樣歸一化結(jié)果Fig.3 Normalized result of input uncertainty parameter sampling

    當(dāng)輸入?yún)?shù)分布類型為正態(tài)分布,抽樣樣本在分布范圍內(nèi)呈現(xiàn)為中間密集,兩邊稀疏;當(dāng)輸入?yún)?shù)分布類型為均勻分布,抽樣樣本在分布范圍內(nèi)呈均勻態(tài)勢。由圖3對比表4發(fā)現(xiàn),選定各輸入?yún)?shù)的概率分布類型準(zhǔn)確無誤。

    3 不確定性與敏感性分析

    3.1 不確定性量化分析

    針對SNCLFR-10 堆的UTOP 事故不確定性分析,通過不確定性輸入?yún)?shù)的組合量化反應(yīng)性、堆芯功率、燃料中心溫度以及包殼溫度的不確定性。圖6~9 為不同目標(biāo)參數(shù)的變化情況,目標(biāo)參數(shù)的名義值均處于95/95不確定包絡(luò)帶內(nèi),不確定性分析結(jié)果具有合理性,且不確定度上限沒有出現(xiàn)過于保守的情況;對比不同目標(biāo)參數(shù)的不確定性包絡(luò)帶,可知不同時(shí)刻輸入?yún)?shù)的不確定性傳播到目標(biāo)參數(shù)的不確定性不同。

    圖4、5 為100 組工況下反應(yīng)性與功率的變化情況,發(fā)生UTOP事故的2 s內(nèi),總反應(yīng)性、總功率急劇增強(qiáng)。此階段不確定帶較小,隨后由于燃料和冷卻劑溫度上升,負(fù)反應(yīng)性反饋增強(qiáng),且反應(yīng)性反饋系數(shù)不確定度較大,促使反應(yīng)性、功率的不確定帶相對變寬,隨著負(fù)反應(yīng)性反饋?zhàn)饔?,反?yīng)性與功率逐漸進(jìn)入一個(gè)新的穩(wěn)態(tài)。圖6、7為100組工況下包殼溫度、燃料溫度的變化情況,UTOP事故發(fā)生后,隨著功率的增大,燃料與包殼的溫度迅速上升,之后由于負(fù)反饋效應(yīng)增強(qiáng),且反應(yīng)性反饋系數(shù)不確定度較大,燃料與包殼溫度的不確定帶愈加明顯,隨著功率下降,燃料 與包殼溫度逐漸進(jìn)入穩(wěn)態(tài)。

    圖4 反應(yīng)性變化情況(100組工況)(a)與不確定性包絡(luò)帶(b)Fig.4 Changes of reactivity(100 sets of outputs)(a)and uncertainty bands(b)

    圖5 功率變化情況(100組工況)(a)與不確定性包絡(luò)帶(b)Fig.5 Changes of power(100 sets of outputs)(a)and uncertainty bands(b)

    圖6 包殼溫度變化情況(100組工況)(a)與不確定性包絡(luò)帶(b)Fig.6 Changes of cladding temperature(100 sets of outputs)(a)and uncertainty bands(b)

    圖7 燃料溫度變化情況(100組工況)(a)與不確定性包絡(luò)帶(b)Fig.7 Changes of fuel temperature(100 sets of outputs)(a)and uncertainty bands(b)

    基于非參數(shù)統(tǒng)計(jì)理論,在SNCLFR-10 堆的UTOP 事故中取100 組工況的次大值作為目標(biāo)參數(shù)容忍上限值。計(jì)算可得,包殼峰值溫度的95/95容忍上限值為1 208.66 K,燃料峰值溫度的95/95 容忍上限值為2 757.25 K。

    3.2 敏感性分析

    針對SNCLFR-10 堆UTOP 事故的敏感性分析,本工作基于SIMLAB 程序采用Spearman 秩相關(guān)分析方法研究輸入?yún)?shù)對燃料最高溫度以及包殼最高溫度的敏感度。Spearman 秩相關(guān)分析法用于全局敏感性分析,可分析多個(gè)輸入?yún)?shù)對目標(biāo)參數(shù)的影響。相關(guān)系數(shù)取值范圍為-1~+1,正號表示正相關(guān),負(fù)號表示負(fù)相關(guān),零表示不相關(guān),相關(guān)系數(shù)絕對值越靠近1,敏感性越強(qiáng)。Spearman 秩相關(guān)系數(shù)表達(dá)式如下:

    式中:ρs為Spearman秩相關(guān)系數(shù);RXi為Xi在X中的大小排序;RYi為Yi在Y中的大小排序;n為樣本數(shù)。

    圖8、9 分別為輸入?yún)?shù)與包殼最高溫度、燃料最高溫度的Spearman秩相關(guān)系數(shù),深灰色柱狀條表示輸入?yún)?shù)與目標(biāo)參數(shù)呈正相關(guān),淺灰色柱狀條表示輸入?yún)?shù)與目標(biāo)參數(shù)呈負(fù)相關(guān)。由圖8、9 可得,燃料制造公差輸入?yún)?shù)中,包殼內(nèi)直徑與燃料最高溫度、包殼最高溫度的相關(guān)系數(shù)分別為0.345 1、-0.352 2,呈中等相關(guān);熱工水力輸入?yún)?shù)中,燃料平均溫度與燃料最高溫度、包殼最高溫度的相關(guān)系數(shù)分別為0.912、0.908 6,呈極強(qiáng)相關(guān),這是由于燃料平均溫度自身存在較大的不確定度,且對反應(yīng)性反饋有直接影響;中子物理輸入?yún)?shù)中,第1組衰變常數(shù)與目標(biāo)參數(shù)的相關(guān)系數(shù)分別為-0.619 4、0.565 2,功率與目標(biāo)參數(shù)的相關(guān)系數(shù)分別為0.386、-0.241 1,均呈中度相關(guān);燃料多普勒常數(shù)與目標(biāo)參數(shù)的的相關(guān)系數(shù)分別為0.934 5、0.945 6,徑向膨脹系數(shù)與目標(biāo)參數(shù)的的相關(guān)系數(shù)分別為0.703、0.742 2,軸向膨脹系數(shù)與目標(biāo)參數(shù)的的相關(guān)系數(shù)分別為0.477 99、0.554,說明反應(yīng)性反饋系數(shù)對目標(biāo)參數(shù)具有較高的靈敏度,這是由于其自身存在較大的不確定度,且反應(yīng)性反饋系數(shù)是引發(fā)UTOP工況的主要參數(shù)。

    圖8 不確定性輸入?yún)?shù)與包殼最高溫度的Spearman秩相關(guān)系數(shù)Fig.8 Correlation coefficient of Spearman rank between the input parameters and peak cladding temperature

    圖9 不確定性輸入?yún)?shù)與燃料最高溫度的Spearman秩相關(guān)系數(shù)Fig.9 Correlation coefficient of Spearman rank between the input parameters and peak fuel temperature

    綜上所述,在SNCLFR-10 堆的UTOP 事故中,顯著影響燃料、包殼峰值溫度的輸入?yún)?shù)包括包殼直徑、穩(wěn)態(tài)燃料平均溫度、初始功率、第1 組衰變常數(shù)以及反應(yīng)性反饋系數(shù);其中燃料平均溫度以及反應(yīng)性反饋系數(shù)對瞬態(tài)安全參數(shù)的影響最大。

    4 結(jié)語

    本工作基于CFD/PFS 程序開展了小型自然循環(huán)鉛鉍快堆SNCLFR-10 的無保護(hù)超功率瞬態(tài)事故不確定性分析,并研究了輸入?yún)?shù)對瞬態(tài)安全參數(shù)的影響程度。結(jié)論如下:

    1)基于非參數(shù)統(tǒng)計(jì)法的不確定性計(jì)算,得到的目標(biāo)參數(shù)95/95不確定帶能夠較好地包絡(luò)名義值;對比目標(biāo)參數(shù)不確定性包絡(luò)帶,可得不確定性輸入?yún)?shù)傳播到輸出參數(shù)的不確定性隨時(shí)間而變化;包殼峰值溫度95/95 容忍上限值為1 208.66 K,燃料峰值溫度的95/95容忍上限值為2 757.25 K。

    2)通過Spearman 相關(guān)系數(shù)法對SNCLFR-10 的UTOP 事故進(jìn)行敏感性分析。結(jié)果表明:包殼內(nèi)直徑、穩(wěn)態(tài)燃料平均溫度、初始功率、第1 組衰變常數(shù)以及反應(yīng)性反饋系數(shù)顯著影響燃料和包殼峰值溫度,其中燃料平均溫度和反應(yīng)性反饋系數(shù)對瞬態(tài)安全參數(shù)的敏感度最強(qiáng)。

    作者貢獻(xiàn)聲明余清遠(yuǎn):實(shí)施研究,文章撰寫;漆靜雯:采集數(shù)據(jù),分析/解釋數(shù)據(jù);趙鵬程:統(tǒng)計(jì)分析,獲取研究經(jīng)費(fèi);趙亞楠:對文章的知識性內(nèi)容作批評性審閱;于濤:行政、技術(shù)支持,指導(dǎo),支持性貢獻(xiàn)。

    猜你喜歡
    分析
    禽大腸桿菌病的分析、診斷和防治
    隱蔽失效適航要求符合性驗(yàn)證分析
    電力系統(tǒng)不平衡分析
    電子制作(2018年18期)2018-11-14 01:48:24
    電力系統(tǒng)及其自動化發(fā)展趨勢分析
    經(jīng)濟(jì)危機(jī)下的均衡與非均衡分析
    對計(jì)劃生育必要性以及其貫徹實(shí)施的分析
    GB/T 7714-2015 與GB/T 7714-2005對比分析
    出版與印刷(2016年3期)2016-02-02 01:20:11
    網(wǎng)購中不良現(xiàn)象分析與應(yīng)對
    中西醫(yī)結(jié)合治療抑郁癥100例分析
    偽造有價(jià)證券罪立法比較分析
    一级黄色大片毛片| a级毛片免费高清观看在线播放| 欧美一区二区国产精品久久精品| 精品一区二区三区视频在线| 中文字幕人妻熟人妻熟丝袜美| netflix在线观看网站| 国产一区二区三区视频了| 精品一区二区免费观看| 日韩av在线大香蕉| 校园春色视频在线观看| 婷婷亚洲欧美| 日本 av在线| 丰满乱子伦码专区| 国内精品久久久久久久电影| 少妇的逼水好多| 亚洲国产高清在线一区二区三| 观看美女的网站| 日韩欧美精品免费久久| 乱码一卡2卡4卡精品| 免费无遮挡裸体视频| 亚洲精品成人久久久久久| 一级毛片久久久久久久久女| 一进一出抽搐动态| 国产伦精品一区二区三区四那| 精品一区二区三区视频在线观看免费| 亚洲18禁久久av| 亚洲精品久久国产高清桃花| 午夜久久久久精精品| 精品福利观看| 国产精品av视频在线免费观看| 一夜夜www| 日本黄色片子视频| 亚洲最大成人中文| 男女下面进入的视频免费午夜| 亚洲av一区综合| xxxwww97欧美| 两性午夜刺激爽爽歪歪视频在线观看| 国产精品无大码| 最近视频中文字幕2019在线8| 99精品久久久久人妻精品| 亚洲在线观看片| 精品不卡国产一区二区三区| 精品一区二区免费观看| 久久国内精品自在自线图片| 免费在线观看日本一区| 在线观看av片永久免费下载| 亚洲精品国产成人久久av| 国产精品一区www在线观看 | 欧美bdsm另类| 久久午夜福利片| 欧美高清性xxxxhd video| 欧美日韩瑟瑟在线播放| 亚洲 国产 在线| 神马国产精品三级电影在线观看| 老女人水多毛片| 免费黄网站久久成人精品| 一卡2卡三卡四卡精品乱码亚洲| 日本三级黄在线观看| 欧美一级a爱片免费观看看| 欧美最黄视频在线播放免费| 国产亚洲精品av在线| 日韩,欧美,国产一区二区三区 | 中文字幕av成人在线电影| 久久精品国产自在天天线| 亚洲av第一区精品v没综合| 成年女人毛片免费观看观看9| 18+在线观看网站| 99国产极品粉嫩在线观看| 舔av片在线| 免费在线观看成人毛片| 一个人看视频在线观看www免费| 深夜a级毛片| 精品免费久久久久久久清纯| 久久久久久久久久久丰满 | 制服丝袜大香蕉在线| 男人狂女人下面高潮的视频| 精品福利观看| 亚洲av成人av| xxxwww97欧美| 波野结衣二区三区在线| 欧美黑人欧美精品刺激| 99热6这里只有精品| 色综合婷婷激情| 亚洲真实伦在线观看| 99九九线精品视频在线观看视频| 老司机午夜福利在线观看视频| 不卡视频在线观看欧美| 我的女老师完整版在线观看| 国产精华一区二区三区| 欧美最黄视频在线播放免费| 亚洲自拍偷在线| 精品免费久久久久久久清纯| 国产综合懂色| 成人国产一区最新在线观看| 国产高清激情床上av| 联通29元200g的流量卡| 少妇猛男粗大的猛烈进出视频 | 亚洲av一区综合| 午夜福利在线观看吧| 国产淫片久久久久久久久| 美女高潮的动态| 中文字幕人妻熟人妻熟丝袜美| 中文字幕av在线有码专区| 欧美精品啪啪一区二区三区| 亚洲中文字幕日韩| 91久久精品电影网| 亚洲成人中文字幕在线播放| 国产精品,欧美在线| 一级毛片久久久久久久久女| 亚洲久久久久久中文字幕| 女同久久另类99精品国产91| 精品午夜福利在线看| av女优亚洲男人天堂| 国产激情偷乱视频一区二区| 亚洲成人中文字幕在线播放| 在线观看午夜福利视频| a级一级毛片免费在线观看| 此物有八面人人有两片| 最近最新免费中文字幕在线| 啦啦啦韩国在线观看视频| 成人毛片a级毛片在线播放| 3wmmmm亚洲av在线观看| 久久午夜福利片| 免费看日本二区| 99久久精品一区二区三区| 亚洲精品在线观看二区| 久久久久国内视频| 亚洲中文日韩欧美视频| 精华霜和精华液先用哪个| videossex国产| 黄色欧美视频在线观看| 亚洲成人精品中文字幕电影| 亚洲va在线va天堂va国产| 亚洲欧美日韩高清在线视频| 成人国产综合亚洲| 男女做爰动态图高潮gif福利片| 亚洲欧美清纯卡通| 国产亚洲精品综合一区在线观看| 少妇丰满av| 亚洲va在线va天堂va国产| 成人欧美大片| 日韩中字成人| 简卡轻食公司| 两性午夜刺激爽爽歪歪视频在线观看| 国产精品一区二区三区四区久久| 无人区码免费观看不卡| 精品福利观看| a级毛片免费高清观看在线播放| 91精品国产九色| 波野结衣二区三区在线| 久久精品综合一区二区三区| 九九爱精品视频在线观看| 亚洲中文字幕一区二区三区有码在线看| 久久久久久久午夜电影| 日韩欧美在线二视频| av在线蜜桃| 国产精品自产拍在线观看55亚洲| 99久久成人亚洲精品观看| 国产免费一级a男人的天堂| 亚洲无线在线观看| 高清在线国产一区| 国产男靠女视频免费网站| 欧美激情在线99| 观看免费一级毛片| 亚洲成人精品中文字幕电影| 成年人黄色毛片网站| 人人妻人人看人人澡| 午夜福利欧美成人| 国产淫片久久久久久久久| 国产精品99久久久久久久久| 日韩欧美一区二区三区在线观看| 非洲黑人性xxxx精品又粗又长| 婷婷亚洲欧美| 少妇的逼好多水| 成年免费大片在线观看| 亚洲电影在线观看av| 狂野欧美激情性xxxx在线观看| 老司机午夜福利在线观看视频| 久久久久久久久久成人| 日本与韩国留学比较| 中文字幕人妻熟人妻熟丝袜美| 中国美白少妇内射xxxbb| 亚洲中文字幕日韩| 国产一区二区三区av在线 | 无人区码免费观看不卡| 国产色婷婷99| 熟女电影av网| 少妇被粗大猛烈的视频| 18禁在线播放成人免费| av女优亚洲男人天堂| 亚洲欧美精品综合久久99| 少妇的逼水好多| 国产视频一区二区在线看| 久久国产乱子免费精品| 综合色av麻豆| 中出人妻视频一区二区| 国产熟女欧美一区二区| a级毛片免费高清观看在线播放| 99热这里只有精品一区| 国产精品久久久久久亚洲av鲁大| a在线观看视频网站| 久久久久久久久中文| 美女大奶头视频| 桃红色精品国产亚洲av| 少妇熟女aⅴ在线视频| 88av欧美| 欧美另类亚洲清纯唯美| 欧美一区二区精品小视频在线| 欧美成人性av电影在线观看| 日韩一区二区视频免费看| 18禁黄网站禁片免费观看直播| 国产精华一区二区三区| 久久精品国产亚洲av涩爱 | 美女免费视频网站| 亚洲成人久久爱视频| 白带黄色成豆腐渣| 欧美性猛交╳xxx乱大交人| 午夜久久久久精精品| av在线观看视频网站免费| 久久午夜福利片| 亚洲中文日韩欧美视频| 88av欧美| 亚洲不卡免费看| 日韩精品青青久久久久久| 国产精品久久久久久av不卡| 内地一区二区视频在线| 99视频精品全部免费 在线| 国产欧美日韩精品亚洲av| 精品无人区乱码1区二区| 亚洲最大成人中文| 亚洲国产欧洲综合997久久,| 简卡轻食公司| 看黄色毛片网站| 精品人妻视频免费看| 一个人观看的视频www高清免费观看| 国产综合懂色| 久久精品人妻少妇| 欧美性猛交╳xxx乱大交人| 亚洲欧美精品综合久久99| 精品不卡国产一区二区三区| 黄色女人牲交| 最近最新免费中文字幕在线| 亚洲第一电影网av| 国产黄片美女视频| 在线观看一区二区三区| 嫩草影院新地址| 如何舔出高潮| 两性午夜刺激爽爽歪歪视频在线观看| 国产精品久久久久久亚洲av鲁大| 亚洲狠狠婷婷综合久久图片| 精品一区二区三区视频在线| 国产伦精品一区二区三区四那| 国产精品一区二区性色av| 亚洲av成人av| 别揉我奶头~嗯~啊~动态视频| 中文在线观看免费www的网站| 日韩中文字幕欧美一区二区| 夜夜爽天天搞| 亚洲va日本ⅴa欧美va伊人久久| 国产成人福利小说| 久久久久久久亚洲中文字幕| 美女免费视频网站| 成人鲁丝片一二三区免费| 亚洲午夜理论影院| 久久久成人免费电影| 美女高潮的动态| 少妇熟女aⅴ在线视频| 久久久久久伊人网av| 国产单亲对白刺激| av中文乱码字幕在线| 97热精品久久久久久| 欧美激情国产日韩精品一区| 日韩强制内射视频| 国产又黄又爽又无遮挡在线| 草草在线视频免费看| 在线看三级毛片| 1000部很黄的大片| 国内精品宾馆在线| 夜夜夜夜夜久久久久| 成人特级黄色片久久久久久久| 亚洲av熟女| 小蜜桃在线观看免费完整版高清| 国产av不卡久久| 久久久久久久久久久丰满 | 99久久精品国产国产毛片| 两性午夜刺激爽爽歪歪视频在线观看| 最近视频中文字幕2019在线8| 啦啦啦观看免费观看视频高清| 久久婷婷人人爽人人干人人爱| 欧美丝袜亚洲另类 | 美女高潮的动态| 一区二区三区高清视频在线| 国内精品宾馆在线| 99视频精品全部免费 在线| 两个人视频免费观看高清| 看黄色毛片网站| 十八禁网站免费在线| 高清在线国产一区| www日本黄色视频网| 久久久久性生活片| 91午夜精品亚洲一区二区三区 | 欧美又色又爽又黄视频| 亚洲人与动物交配视频| 91久久精品电影网| 能在线免费观看的黄片| 国产午夜福利久久久久久| 久久午夜福利片| 天天一区二区日本电影三级| 日本a在线网址| 日韩欧美在线乱码| 国产欧美日韩精品亚洲av| 中文字幕人妻熟人妻熟丝袜美| eeuss影院久久| 麻豆精品久久久久久蜜桃| 国产免费一级a男人的天堂| 国产成人影院久久av| 精品福利观看| 一个人观看的视频www高清免费观看| 熟妇人妻久久中文字幕3abv| 我要看日韩黄色一级片| 91在线观看av| 99九九线精品视频在线观看视频| 国产午夜精品论理片| 亚洲精品国产成人久久av| 日韩大尺度精品在线看网址| 欧美精品国产亚洲| 成人美女网站在线观看视频| 91麻豆精品激情在线观看国产| 日本撒尿小便嘘嘘汇集6| 亚洲av二区三区四区| 国产高清视频在线观看网站| 久久久色成人| 在线播放无遮挡| 18禁黄网站禁片免费观看直播| 午夜福利在线观看吧| 免费在线观看日本一区| 99国产极品粉嫩在线观看| 一区二区三区四区激情视频 | 乱人视频在线观看| 91麻豆av在线| 亚洲中文字幕一区二区三区有码在线看| a级毛片a级免费在线| 人妻久久中文字幕网| 亚洲电影在线观看av| 综合色av麻豆| 特大巨黑吊av在线直播| 国内精品美女久久久久久| 午夜亚洲福利在线播放| 欧美人与善性xxx| a级毛片a级免费在线| 国产激情偷乱视频一区二区| 日韩强制内射视频| 日日摸夜夜添夜夜添av毛片 | 欧美激情久久久久久爽电影| 69人妻影院| 国产精品99久久久久久久久| 91在线精品国自产拍蜜月| 熟妇人妻久久中文字幕3abv| 日日干狠狠操夜夜爽| 亚洲欧美日韩东京热| 天美传媒精品一区二区| 欧美+日韩+精品| 深爱激情五月婷婷| 欧美人与善性xxx| 偷拍熟女少妇极品色| 亚洲无线观看免费| 乱人视频在线观看| 欧美激情国产日韩精品一区| а√天堂www在线а√下载| 免费看光身美女| 日本a在线网址| 欧美精品啪啪一区二区三区| 色播亚洲综合网| 亚洲国产精品合色在线| 99久久无色码亚洲精品果冻| 又紧又爽又黄一区二区| 在线观看午夜福利视频| 日本-黄色视频高清免费观看| 高清毛片免费观看视频网站| 午夜福利在线观看免费完整高清在 | 国产精品国产高清国产av| 一夜夜www| 国产精品乱码一区二三区的特点| av视频在线观看入口| 国产毛片a区久久久久| 99国产精品一区二区蜜桃av| 在线a可以看的网站| 日韩高清综合在线| 精品一区二区三区视频在线观看免费| 午夜爱爱视频在线播放| 亚洲精品在线观看二区| 免费看av在线观看网站| 日韩精品有码人妻一区| 亚洲美女搞黄在线观看 | 久久午夜福利片| 在线a可以看的网站| 性色avwww在线观看| 久久精品91蜜桃| 真人一进一出gif抽搐免费| 不卡一级毛片| 99精品久久久久人妻精品| 久久热精品热| 亚洲人成伊人成综合网2020| 国产 一区精品| 欧美日韩亚洲国产一区二区在线观看| 精品国产三级普通话版| 精品久久久久久久久久久久久| 亚洲一级一片aⅴ在线观看| 男女那种视频在线观看| 男女边吃奶边做爰视频| 国产伦一二天堂av在线观看| 两个人视频免费观看高清| 日韩精品中文字幕看吧| 日韩高清综合在线| 99久久久亚洲精品蜜臀av| 狂野欧美白嫩少妇大欣赏| 午夜免费成人在线视频| 精品午夜福利视频在线观看一区| 国内揄拍国产精品人妻在线| 最近最新中文字幕大全电影3| 久久中文看片网| 日韩大尺度精品在线看网址| 亚洲无线在线观看| 亚洲va日本ⅴa欧美va伊人久久| 夜夜爽天天搞| 99视频精品全部免费 在线| 亚洲四区av| ponron亚洲| 欧美色欧美亚洲另类二区| 精品人妻视频免费看| 校园春色视频在线观看| 免费不卡的大黄色大毛片视频在线观看 | 欧美潮喷喷水| 有码 亚洲区| 淫妇啪啪啪对白视频| 日韩欧美一区二区三区在线观看| 欧美在线一区亚洲| 亚洲自偷自拍三级| 一本精品99久久精品77| 国产伦精品一区二区三区视频9| 日日撸夜夜添| 亚洲av第一区精品v没综合| 欧美一级a爱片免费观看看| 久久亚洲精品不卡| 亚洲人成伊人成综合网2020| 99视频精品全部免费 在线| 69人妻影院| 国产在线男女| 国产精品人妻久久久影院| 色尼玛亚洲综合影院| 精品人妻1区二区| 亚洲欧美日韩高清专用| 日本成人三级电影网站| 91久久精品国产一区二区三区| av在线老鸭窝| 亚洲美女视频黄频| 久久精品国产99精品国产亚洲性色| 悠悠久久av| 一本久久中文字幕| 国产精品国产高清国产av| 狂野欧美激情性xxxx在线观看| 国产综合懂色| 久久精品国产清高在天天线| 一卡2卡三卡四卡精品乱码亚洲| 2021天堂中文幕一二区在线观| 一边摸一边抽搐一进一小说| 亚洲,欧美,日韩| 十八禁网站免费在线| 亚洲三级黄色毛片| 国国产精品蜜臀av免费| 99久久精品一区二区三区| 他把我摸到了高潮在线观看| 九九爱精品视频在线观看| 国产三级在线视频| 久久久色成人| 免费无遮挡裸体视频| 琪琪午夜伦伦电影理论片6080| 一本精品99久久精品77| 一区福利在线观看| 成年免费大片在线观看| 又粗又爽又猛毛片免费看| 日韩欧美精品免费久久| 欧美激情久久久久久爽电影| av在线天堂中文字幕| 波多野结衣巨乳人妻| 国产亚洲91精品色在线| 九色国产91popny在线| 久久久久久伊人网av| 免费黄网站久久成人精品| 五月伊人婷婷丁香| 女人被狂操c到高潮| 日韩一本色道免费dvd| 真人一进一出gif抽搐免费| 国产精品久久久久久久电影| 精品久久久久久久末码| 婷婷亚洲欧美| 岛国在线免费视频观看| 老司机深夜福利视频在线观看| 国产高清有码在线观看视频| 日韩精品中文字幕看吧| 午夜福利在线观看吧| 国产精品亚洲一级av第二区| 亚洲精华国产精华精| 看片在线看免费视频| 午夜日韩欧美国产| 久久精品影院6| 十八禁国产超污无遮挡网站| 亚洲美女视频黄频| 人妻久久中文字幕网| 日韩中字成人| 精品日产1卡2卡| 亚洲成人中文字幕在线播放| 91精品国产九色| 色噜噜av男人的天堂激情| 欧美另类亚洲清纯唯美| 久久精品国产自在天天线| 色av中文字幕| 国产主播在线观看一区二区| 日韩欧美在线二视频| 久久精品国产亚洲av香蕉五月| 国产久久久一区二区三区| 久久亚洲精品不卡| 亚洲精品乱码久久久v下载方式| 久久婷婷人人爽人人干人人爱| 欧美国产日韩亚洲一区| 国产蜜桃级精品一区二区三区| 成人亚洲精品av一区二区| 男女下面进入的视频免费午夜| 九九爱精品视频在线观看| 久久精品国产亚洲av香蕉五月| 亚洲天堂国产精品一区在线| 欧美性感艳星| 久久精品91蜜桃| 成人鲁丝片一二三区免费| 国产精品国产三级国产av玫瑰| 亚洲无线观看免费| 国产亚洲精品久久久com| 无遮挡黄片免费观看| 国产色爽女视频免费观看| 在线观看美女被高潮喷水网站| 三级男女做爰猛烈吃奶摸视频| 99热6这里只有精品| 九九久久精品国产亚洲av麻豆| 看黄色毛片网站| 色精品久久人妻99蜜桃| 亚洲av中文av极速乱 | 国产 一区精品| 成人美女网站在线观看视频| 中文字幕熟女人妻在线| 中文亚洲av片在线观看爽| 最新在线观看一区二区三区| 亚洲av免费在线观看| 亚洲va日本ⅴa欧美va伊人久久| 久久久色成人| 少妇丰满av| 国产伦精品一区二区三区四那| 久久久久久久精品吃奶| 女同久久另类99精品国产91| 精品欧美国产一区二区三| 午夜福利在线在线| 真实男女啪啪啪动态图| 狂野欧美激情性xxxx在线观看| 久久精品国产亚洲av涩爱 | 黄色视频,在线免费观看| 久久香蕉精品热| 又爽又黄无遮挡网站| 欧美精品国产亚洲| 看十八女毛片水多多多| 精品无人区乱码1区二区| 精品一区二区免费观看| 在线观看舔阴道视频| 日韩高清综合在线| 精品人妻视频免费看| 日韩大尺度精品在线看网址| 一个人观看的视频www高清免费观看| 成人国产一区最新在线观看| 亚洲天堂国产精品一区在线| 制服丝袜大香蕉在线| 国产女主播在线喷水免费视频网站 | 精品日产1卡2卡| 琪琪午夜伦伦电影理论片6080| 变态另类成人亚洲欧美熟女| 最新在线观看一区二区三区| 亚洲av成人av| 亚洲色图av天堂| 日本撒尿小便嘘嘘汇集6| 亚州av有码| 又黄又爽又免费观看的视频| 熟女人妻精品中文字幕| 国产精品一及| 国产伦在线观看视频一区| 精品一区二区三区视频在线| 男女做爰动态图高潮gif福利片| 韩国av一区二区三区四区| 国产精华一区二区三区| 欧美bdsm另类| 久久久久久久久中文| 国产亚洲av嫩草精品影院| 日本色播在线视频| 波多野结衣巨乳人妻| 午夜免费激情av| 欧美bdsm另类| 亚洲午夜理论影院| 久久久色成人| 精品国产三级普通话版| 久久久久免费精品人妻一区二区| 麻豆国产av国片精品| 九九久久精品国产亚洲av麻豆| 亚洲,欧美,日韩| 熟女人妻精品中文字幕| 丰满人妻一区二区三区视频av| 国内精品一区二区在线观看| x7x7x7水蜜桃| 十八禁网站免费在线| 又紧又爽又黄一区二区|