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

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

    本工作是針對中國科學(xué)技術(shù)大學(xué)和南華大學(xué)聯(lián)合設(shè)計的小型自然循環(huán)鉛基鉛鉍快堆(Small Natural Circulation Lead-bismuth Cooled Fast Reactor 10 MWth,SNCLFR-10),首先,通過文獻(xiàn)資料以及工程經(jīng)驗判斷,考慮三類不確定性輸入?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 計算模型

    1.1 分析模型

    SNCLFR-10 堆是由中國科學(xué)技術(shù)大學(xué)和南華大學(xué)聯(lián)合設(shè)計的小型自然循環(huán)鉛鉍快堆,具有良好的現(xiàn)實可行性、安全可靠性、實驗靈活性和技術(shù)延續(xù)性[22]。圖1 為SNCLFR-10 堆的結(jié)構(gòu)設(shè)計圖,表1、2給出了反應(yīng)堆設(shè)計參數(shù)及堆芯主要設(shè)計參數(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è)計參數(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)突然引入一個正反應(yīng)性,導(dǎo)致功率急劇升高的瞬態(tài)過程,此過程中反應(yīng)堆停堆保護(hù)系統(tǒng)失效。UTOP 事故的瞬態(tài)條件:反應(yīng)堆初始工況為穩(wěn)態(tài)運行,假設(shè)最大價值控制棒失去控制而不可抽出,故在2 s 內(nèi)迅速引入1$正反應(yīng)性,整個過程無停堆保護(hù)動作。

    1.2 CFD/PFS程序

    CFD/PFS 是基于CFD 程序FLUENT 的用戶自定義函數(shù)(User Defined Function,UDF)耦合點堆中子動力學(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)計預(yù)處理子系統(tǒng)、模型執(zhí)行子系統(tǒng)、統(tǒng)計后處理子系統(tǒng),廣泛用于不確定性和整體敏感性研究[23]。采用與SIMMER 程序code-to-code的方法進(jìn)行CFD/PFS 的對比驗證,CFD/PFS 程序開發(fā)與詳細(xì)驗證見參考文獻(xiàn)[24]。

    CFD/PFS程序分為物理熱工耦合計算模型與不確定性分析計算,其計算流程如圖2 所示。物理熱工耦合計算模型包括CFD模型、點堆中子動力學(xué)模型、燃料棒熱傳導(dǎo)模型。PKM是三維時空動力學(xué)的簡化模型,采用6 組緩發(fā)中子的點堆動力學(xué)方程主要用于計算反應(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ā)中子的份額;Λ是中子代時間;λi是第i組緩發(fā)中子的衰變常數(shù);Ci是i組緩發(fā)中子先驅(qū)核的濃度。

    再根據(jù)反應(yīng)性通過點堆動力學(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時刻的燃料、冷卻劑平均溫度。

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

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

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

    2.1 不確定性分析方法

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

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

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

    由wilks 公式計算出的滿足不同置信度和容許區(qū)間的最小樣本量n值如表3 所示。根據(jù)非參數(shù)統(tǒng)計理論,最小樣本量選取59 時,目標(biāo)參數(shù)的容忍上限值取最大值;最小樣本量選取93 時,目標(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)計算中的不確定性輸入?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)?shù),并確定其范圍和概率密度分布。選取的不確定輸入?yún)?shù)分為三個部分:首先是幾何輸入?yún)?shù),一些參數(shù)測量有一定的不安全性故缺乏,只能由專家假設(shè)估算;其次是中子學(xué)參數(shù),因為實驗測量系統(tǒng)及統(tǒng)計誤差,故其引入的不確定性需考慮;最后是熱工水力參數(shù),程序、模型及經(jīng)驗關(guān)系式等都會引入不確定性,故選取其中最典型5 個參數(shù)。表4 列出了25 個不確定輸入?yún)?shù)信息,輸入?yún)?shù)服從均勻分布或正態(tài)分布,其中不確定性度表示參數(shù)名義值(設(shè)計值)和參數(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%。中子代時間、緩發(fā)中子份額、衰變常數(shù)以及初始功率的不確定度參考其他類似反應(yīng)堆的不確定度;對于反應(yīng)性反饋系數(shù),考慮它包含了一系列不確定性(材料屬性、幾何形狀、截面數(shù)據(jù)以及蒙特卡羅計算),為保證其保守性,將反應(yīng)性反饋系數(shù)的不確定度設(shè)為10%[26]。熱工水力現(xiàn)象極其復(fù)雜,熱工物性參數(shù)大多是根據(jù)實驗數(shù)據(jù)擬合以及經(jīng)驗關(guān)系式所得,因此熱工水力模型輸入?yún)?shù)通常存在較大的不確定度[27]。

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

    基于SIMLAB程序采用拉丁超立方方法實現(xiàn)不確定性輸入?yún)?shù)抽樣,種子數(shù)設(shè)定為98 765。根據(jù)wilks 非參數(shù)統(tǒng)計法,在95%置信水平、95%概率水平下,所需最小計算次數(shù)為93次??紤]到抽樣可能存在壞點、程序可能運行失效等原因,抽樣數(shù)目設(shè)定為100 組,其中每一組工況包括25 個輸入?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ò)帶,可知不同時刻輸入?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)性反饋作用,反應(yīng)性與功率逐漸進(jìn)入一個新的穩(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)計理論,在SNCLFR-10 堆的UTOP 事故中取100 組工況的次大值作為目標(biāo)參數(shù)容忍上限值。計算可得,包殼峰值溫度的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)分析法用于全局敏感性分析,可分析多個輸入?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)計法的不確定性計算,得到的目標(biāo)參數(shù)95/95不確定帶能夠較好地包絡(luò)名義值;對比目標(biāo)參數(shù)不確定性包絡(luò)帶,可得不確定性輸入?yún)?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ù)據(jù),分析/解釋數(shù)據(jù);趙鵬程:統(tǒng)計分析,獲取研究經(jīng)費;趙亞楠:對文章的知識性內(nèi)容作批評性審閱;于濤:行政、技術(shù)支持,指導(dǎo),支持性貢獻(xiàn)。

    猜你喜歡
    分析
    禽大腸桿菌病的分析、診斷和防治
    隱蔽失效適航要求符合性驗證分析
    電力系統(tǒng)不平衡分析
    電子制作(2018年18期)2018-11-14 01:48:24
    電力系統(tǒng)及其自動化發(fā)展趨勢分析
    經(jīng)濟(jì)危機(jī)下的均衡與非均衡分析
    對計劃生育必要性以及其貫徹實施的分析
    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例分析
    偽造有價證券罪立法比較分析
    中文乱码字字幕精品一区二区三区| 久久精品久久精品一区二区三区| 91aial.com中文字幕在线观看| 新久久久久国产一级毛片| 欧美xxxx性猛交bbbb| 国产一区二区在线观看av| 两个人的视频大全免费| 欧美精品高潮呻吟av久久| 久久精品熟女亚洲av麻豆精品| 亚洲精品久久午夜乱码| 日日爽夜夜爽网站| 青青草视频在线视频观看| 一二三四中文在线观看免费高清| 自拍偷自拍亚洲精品老妇| 久久99精品国语久久久| 2022亚洲国产成人精品| .国产精品久久| 熟女av电影| 精品人妻熟女毛片av久久网站| 在线观看免费视频网站a站| 免费观看的影片在线观看| 久久午夜福利片| 自拍欧美九色日韩亚洲蝌蚪91 | 91精品国产九色| 免费观看性生交大片5| 啦啦啦中文免费视频观看日本| 成人漫画全彩无遮挡| 黑人猛操日本美女一级片| 亚洲一级一片aⅴ在线观看| 女人精品久久久久毛片| 亚洲av国产av综合av卡| 一个人免费看片子| av卡一久久| a 毛片基地| 99热6这里只有精品| 五月天丁香电影| 日本爱情动作片www.在线观看| 26uuu在线亚洲综合色| 久久鲁丝午夜福利片| 免费av不卡在线播放| 国产精品国产三级专区第一集| 少妇丰满av| 极品人妻少妇av视频| 中文字幕人妻丝袜制服| 全区人妻精品视频| 男女啪啪激烈高潮av片| 成人亚洲精品一区在线观看| 国产成人一区二区在线| 亚洲第一av免费看| 亚洲国产精品999| 乱码一卡2卡4卡精品| 自线自在国产av| 亚洲国产欧美日韩在线播放 | 这个男人来自地球电影免费观看 | 亚洲电影在线观看av| 亚洲av国产av综合av卡| 欧美性感艳星| 青春草国产在线视频| 久久亚洲国产成人精品v| 伊人亚洲综合成人网| 国产亚洲欧美精品永久| 亚洲精品成人av观看孕妇| 99久久精品一区二区三区| 交换朋友夫妻互换小说| 青春草亚洲视频在线观看| 免费黄色在线免费观看| 成人午夜精彩视频在线观看| 一区二区三区精品91| 中国国产av一级| 少妇 在线观看| 国产精品国产av在线观看| 午夜日本视频在线| av天堂久久9| 免费在线观看成人毛片| 成年人免费黄色播放视频 | 三级国产精品片| 嫩草影院新地址| 亚洲精品中文字幕在线视频 | 男女国产视频网站| 国产69精品久久久久777片| 欧美高清成人免费视频www| 男人舔奶头视频| 91久久精品国产一区二区三区| 久久精品夜色国产| 尾随美女入室| 亚洲欧美精品专区久久| 在线观看免费高清a一片| 麻豆精品久久久久久蜜桃| 色5月婷婷丁香| 国产精品福利在线免费观看| 国产精品久久久久久av不卡| tube8黄色片| 高清黄色对白视频在线免费看 | 伊人久久精品亚洲午夜| 精品国产一区二区久久| 免费观看的影片在线观看| 欧美97在线视频| 中文天堂在线官网| 国产综合精华液| 亚洲国产成人一精品久久久| 国产av码专区亚洲av| 五月玫瑰六月丁香| 免费在线观看成人毛片| www.av在线官网国产| 又大又黄又爽视频免费| 最近手机中文字幕大全| 成人午夜精彩视频在线观看| 看非洲黑人一级黄片| 国产av国产精品国产| 三上悠亚av全集在线观看 | 欧美xxxx性猛交bbbb| 日日摸夜夜添夜夜添av毛片| 国产精品无大码| 嫩草影院新地址| 九九爱精品视频在线观看| 性高湖久久久久久久久免费观看| av又黄又爽大尺度在线免费看| 在线 av 中文字幕| 少妇被粗大的猛进出69影院 | 国产成人freesex在线| 成人国产麻豆网| 中文字幕精品免费在线观看视频 | 精品亚洲成a人片在线观看| 亚洲电影在线观看av| 午夜久久久在线观看| 久久精品国产亚洲网站| 久久久久久久久久久免费av| 国产男人的电影天堂91| 亚洲精品456在线播放app| 久久久久国产网址| 一边亲一边摸免费视频| 婷婷色av中文字幕| 永久网站在线| 人人澡人人妻人| 蜜桃久久精品国产亚洲av| 久久久久久久久久人人人人人人| 午夜激情久久久久久久| 99热全是精品| h日本视频在线播放| 女人精品久久久久毛片| 精品人妻熟女av久视频| 欧美激情极品国产一区二区三区 | 日韩不卡一区二区三区视频在线| 日韩亚洲欧美综合| 精品人妻熟女毛片av久久网站| 极品人妻少妇av视频| 男人爽女人下面视频在线观看| 人人妻人人看人人澡| 日本午夜av视频| 王馨瑶露胸无遮挡在线观看| 在线观看免费视频网站a站| 国产精品成人在线| 性色avwww在线观看| 两个人免费观看高清视频 | 伊人久久国产一区二区| 亚洲精品久久久久久婷婷小说| 色视频在线一区二区三区| 极品少妇高潮喷水抽搐| 日日啪夜夜撸| 国模一区二区三区四区视频| 两个人免费观看高清视频 | 欧美一级a爱片免费观看看| 久久久亚洲精品成人影院| 午夜免费鲁丝| 国产日韩欧美亚洲二区| 男人舔奶头视频| 日日摸夜夜添夜夜爱| 日日啪夜夜爽| 久久久亚洲精品成人影院| 国产黄片视频在线免费观看| 十八禁高潮呻吟视频 | 久热久热在线精品观看| 久久99一区二区三区| 久久人人爽av亚洲精品天堂| 成人综合一区亚洲| 精品国产一区二区久久| 99热这里只有是精品在线观看| av有码第一页| 97在线视频观看| 中文字幕免费在线视频6| 国产免费福利视频在线观看| 在线观看免费高清a一片| 男人添女人高潮全过程视频| 九九久久精品国产亚洲av麻豆| 色婷婷久久久亚洲欧美| 亚洲经典国产精华液单| 久久午夜综合久久蜜桃| 国产亚洲av片在线观看秒播厂| 国产黄色免费在线视频| 精品久久久精品久久久| 欧美人与善性xxx| 中国三级夫妇交换| 亚洲欧美一区二区三区国产| 一级毛片黄色毛片免费观看视频| av在线播放精品| 久久久久久久亚洲中文字幕| 精品国产乱码久久久久久小说| 日韩免费高清中文字幕av| 另类亚洲欧美激情| 免费不卡的大黄色大毛片视频在线观看| 午夜激情福利司机影院| 国产有黄有色有爽视频| 亚洲欧美日韩卡通动漫| a级片在线免费高清观看视频| 搡女人真爽免费视频火全软件| 国产美女午夜福利| 欧美三级亚洲精品| 亚洲成人手机| 国产成人freesex在线| 国产男人的电影天堂91| 久久久久久久久久久丰满| 国产永久视频网站| 亚洲欧美清纯卡通| 高清视频免费观看一区二区| 一级毛片 在线播放| 国产中年淑女户外野战色| 女人久久www免费人成看片| 久久久国产一区二区| 丰满饥渴人妻一区二区三| a级一级毛片免费在线观看| 91aial.com中文字幕在线观看| 免费观看的影片在线观看| 高清视频免费观看一区二区| 亚洲av福利一区| 国产一级毛片在线| 这个男人来自地球电影免费观看 | 亚洲精品一区蜜桃| 亚洲精品,欧美精品| 夜夜骑夜夜射夜夜干| 99re6热这里在线精品视频| 免费不卡的大黄色大毛片视频在线观看| 国产69精品久久久久777片| 肉色欧美久久久久久久蜜桃| 久久 成人 亚洲| 亚洲三级黄色毛片| 视频区图区小说| 水蜜桃什么品种好| 大话2 男鬼变身卡| 日本爱情动作片www.在线观看| 欧美一级a爱片免费观看看| 亚洲美女黄色视频免费看| 中文字幕精品免费在线观看视频 | 免费看不卡的av| 日韩大片免费观看网站| 久久精品国产亚洲av天美| 男人舔奶头视频| 亚洲精品国产av成人精品| 亚洲在久久综合| 九色成人免费人妻av| 免费黄网站久久成人精品| 永久免费av网站大全| 亚洲av男天堂| 国产精品成人在线| 国产黄频视频在线观看| 黑人高潮一二区| 国产亚洲精品久久久com| 日韩精品免费视频一区二区三区 | 如日韩欧美国产精品一区二区三区 | 欧美另类一区| 亚洲av电影在线观看一区二区三区| 黄色视频在线播放观看不卡| 男女无遮挡免费网站观看| 你懂的网址亚洲精品在线观看| xxx大片免费视频| 99久久综合免费| 国产色爽女视频免费观看| 国产成人精品福利久久| 国产极品粉嫩免费观看在线 | 亚洲内射少妇av| 女性被躁到高潮视频| 国产精品国产av在线观看| 国产精品无大码| 久久综合国产亚洲精品| 久久久久久久国产电影| 午夜福利,免费看| 亚洲伊人久久精品综合| 欧美精品亚洲一区二区| 波野结衣二区三区在线| 亚洲欧美日韩卡通动漫| 久久亚洲国产成人精品v| 久久久国产一区二区| 伊人久久国产一区二区| 亚洲国产av新网站| 国产精品伦人一区二区| av在线观看视频网站免费| 麻豆成人午夜福利视频| 亚洲第一区二区三区不卡| 国产av国产精品国产| 欧美老熟妇乱子伦牲交| 多毛熟女@视频| 乱码一卡2卡4卡精品| 国产欧美日韩综合在线一区二区 | 22中文网久久字幕| av女优亚洲男人天堂| 夫妻性生交免费视频一级片| 亚洲欧美日韩卡通动漫| www.av在线官网国产| 丰满迷人的少妇在线观看| av专区在线播放| 一级毛片 在线播放| 99热这里只有精品一区| 又大又黄又爽视频免费| 成年av动漫网址| 如日韩欧美国产精品一区二区三区 | 国产精品久久久久久av不卡| 黄色一级大片看看| 日日爽夜夜爽网站| 大香蕉久久网| 久久韩国三级中文字幕| 免费大片黄手机在线观看| 欧美性感艳星| 人人妻人人爽人人添夜夜欢视频 | 久久国产精品男人的天堂亚洲 | 国内揄拍国产精品人妻在线| 丝袜脚勾引网站| 高清午夜精品一区二区三区| 夜夜骑夜夜射夜夜干| 欧美最新免费一区二区三区| 少妇猛男粗大的猛烈进出视频| 国产伦理片在线播放av一区| 国产高清国产精品国产三级| 少妇被粗大的猛进出69影院 | 这个男人来自地球电影免费观看 | 色5月婷婷丁香| a级毛片免费高清观看在线播放| 亚洲第一av免费看| 各种免费的搞黄视频| 国产深夜福利视频在线观看| 九九久久精品国产亚洲av麻豆| 三上悠亚av全集在线观看 | 欧美精品高潮呻吟av久久| 亚洲国产精品专区欧美| 又爽又黄a免费视频| 国产一区二区在线观看日韩| 国产爽快片一区二区三区| 日本午夜av视频| 超碰97精品在线观看| 亚洲一区二区三区欧美精品| 亚洲人与动物交配视频| 亚洲av国产av综合av卡| 天堂俺去俺来也www色官网| 国产黄色免费在线视频| 久久99热6这里只有精品| 在线观看av片永久免费下载| 99热6这里只有精品| 男女啪啪激烈高潮av片| 人妻系列 视频| 亚洲精品亚洲一区二区| 亚洲精品第二区| 91久久精品电影网| 最近的中文字幕免费完整| 天天躁夜夜躁狠狠久久av| 亚洲国产精品成人久久小说| 日韩人妻高清精品专区| 一级黄片播放器| 在线观看av片永久免费下载| 少妇熟女欧美另类| 欧美精品一区二区免费开放| 成人二区视频| 久久精品久久久久久噜噜老黄| 国产一级毛片在线| 视频区图区小说| 中文字幕免费在线视频6| 欧美bdsm另类| 九九爱精品视频在线观看| 免费播放大片免费观看视频在线观看| 九九爱精品视频在线观看| 亚洲欧美成人综合另类久久久| 性高湖久久久久久久久免费观看| 国产成人免费无遮挡视频| 女人久久www免费人成看片| 少妇猛男粗大的猛烈进出视频| av有码第一页| 能在线免费看毛片的网站| 欧美日韩综合久久久久久| 一级毛片我不卡| 高清午夜精品一区二区三区| 五月天丁香电影| 纵有疾风起免费观看全集完整版| 久久久久精品性色| 亚洲av欧美aⅴ国产| 国产美女午夜福利| 自拍偷自拍亚洲精品老妇| 啦啦啦啦在线视频资源| 少妇丰满av| 亚洲精品乱码久久久v下载方式| 日韩伦理黄色片| 久久久久久久大尺度免费视频| 99热国产这里只有精品6| 久久久久久久大尺度免费视频| 久久精品久久久久久噜噜老黄| 亚洲美女搞黄在线观看| 人人妻人人澡人人看| 午夜精品国产一区二区电影| 色婷婷av一区二区三区视频| 欧美一级a爱片免费观看看| 一边亲一边摸免费视频| 91精品一卡2卡3卡4卡| 亚洲精品国产色婷婷电影| 80岁老熟妇乱子伦牲交| 在线观看美女被高潮喷水网站| 不卡视频在线观看欧美| 一二三四中文在线观看免费高清| 精品国产一区二区三区久久久樱花| 日韩免费高清中文字幕av| 寂寞人妻少妇视频99o| a 毛片基地| 亚洲精品日韩在线中文字幕| 啦啦啦啦在线视频资源| 亚洲欧美成人综合另类久久久| 97在线人人人人妻| 国产男人的电影天堂91| 国产成人午夜福利电影在线观看| 日韩成人av中文字幕在线观看| 国产探花极品一区二区| 午夜福利视频精品| 在线天堂最新版资源| 老司机亚洲免费影院| 男的添女的下面高潮视频| 亚洲欧美日韩另类电影网站| 少妇人妻一区二区三区视频| 国产在线一区二区三区精| 亚洲一级一片aⅴ在线观看| 大又大粗又爽又黄少妇毛片口| 少妇的逼好多水| 欧美高清成人免费视频www| 视频中文字幕在线观看| 国产男女超爽视频在线观看| 亚洲四区av| 久久久亚洲精品成人影院| 在线观看www视频免费| 一区在线观看完整版| 中文字幕免费在线视频6| 国产一区二区在线观看av| 国产一区有黄有色的免费视频| 色哟哟·www| 欧美97在线视频| 久久久久久久久久久免费av| 日韩免费高清中文字幕av| 人妻夜夜爽99麻豆av| 黄色视频在线播放观看不卡| 欧美日韩在线观看h| 亚洲美女搞黄在线观看| 五月玫瑰六月丁香| 欧美老熟妇乱子伦牲交| 少妇猛男粗大的猛烈进出视频| 精品国产露脸久久av麻豆| 色婷婷久久久亚洲欧美| 亚洲精品,欧美精品| 久久久久久久国产电影| 五月玫瑰六月丁香| 中文资源天堂在线| 亚洲高清免费不卡视频| 两个人的视频大全免费| 热re99久久精品国产66热6| 亚洲精品视频女| 午夜激情福利司机影院| 国产欧美日韩综合在线一区二区 | 国产免费福利视频在线观看| 国产高清有码在线观看视频| 午夜免费男女啪啪视频观看| 午夜久久久在线观看| 九九爱精品视频在线观看| 五月伊人婷婷丁香| 婷婷色av中文字幕| 国产黄频视频在线观看| 一级黄片播放器| 国内少妇人妻偷人精品xxx网站| 成年人午夜在线观看视频| .国产精品久久| 老司机影院成人| 蜜臀久久99精品久久宅男| 国产亚洲一区二区精品| 亚洲美女视频黄频| 青春草视频在线免费观看| 国产视频首页在线观看| 久久鲁丝午夜福利片| 春色校园在线视频观看| av在线观看视频网站免费| 亚洲美女搞黄在线观看| 日韩三级伦理在线观看| 欧美精品高潮呻吟av久久| 国产一区亚洲一区在线观看| 久热久热在线精品观看| 成人黄色视频免费在线看| 丝袜在线中文字幕| 欧美日韩视频精品一区| 国产69精品久久久久777片| 亚洲熟女精品中文字幕| 免费看日本二区| 国产 一区精品| 香蕉精品网在线| 亚洲熟女精品中文字幕| 啦啦啦啦在线视频资源| 日产精品乱码卡一卡2卡三| 国产精品久久久久久精品电影小说| 人体艺术视频欧美日本| 少妇猛男粗大的猛烈进出视频| 黄色一级大片看看| 亚洲欧美日韩卡通动漫| 看十八女毛片水多多多| 日韩成人av中文字幕在线观看| 人妻一区二区av| 青春草亚洲视频在线观看| 国产在线男女| 国产亚洲最大av| 亚洲人成网站在线观看播放| 精品一区在线观看国产| 日韩欧美一区视频在线观看 | 男女无遮挡免费网站观看| 2018国产大陆天天弄谢| 男女无遮挡免费网站观看| 午夜久久久在线观看| 日韩一区二区视频免费看| 各种免费的搞黄视频| 国产日韩欧美在线精品| 国产男人的电影天堂91| 久久狼人影院| 一级毛片电影观看| 日韩成人伦理影院| 精品一区在线观看国产| 国产精品不卡视频一区二区| av天堂久久9| 国产毛片在线视频| av在线播放精品| 老司机影院毛片| 欧美日韩国产mv在线观看视频| 99视频精品全部免费 在线| 亚洲欧美清纯卡通| 日韩欧美 国产精品| 女性被躁到高潮视频| 亚洲激情五月婷婷啪啪| 丰满饥渴人妻一区二区三| 精华霜和精华液先用哪个| 免费看av在线观看网站| 观看av在线不卡| 欧美日韩亚洲高清精品| 在线观看免费日韩欧美大片 | 亚洲精品一区蜜桃| 日韩大片免费观看网站| 亚洲三级黄色毛片| 欧美国产精品一级二级三级 | 成人漫画全彩无遮挡| 高清视频免费观看一区二区| 免费观看性生交大片5| 久久久久久久久久久免费av| 精品卡一卡二卡四卡免费| 亚洲电影在线观看av| 国产淫片久久久久久久久| 男女边摸边吃奶| 久久精品久久久久久久性| 国产真实伦视频高清在线观看| 晚上一个人看的免费电影| 极品人妻少妇av视频| 成人美女网站在线观看视频| 欧美国产精品一级二级三级 | 亚洲精品日本国产第一区| 国产熟女午夜一区二区三区 | 亚洲av欧美aⅴ国产| 熟妇人妻不卡中文字幕| a级毛片免费高清观看在线播放| av天堂中文字幕网| 日本黄色片子视频| 日韩大片免费观看网站| 日韩成人伦理影院| 啦啦啦视频在线资源免费观看| 国产男女超爽视频在线观看| 简卡轻食公司| 亚洲av中文av极速乱| 另类精品久久| 亚洲怡红院男人天堂| 一个人免费看片子| 日本猛色少妇xxxxx猛交久久| 97超碰精品成人国产| 国产高清国产精品国产三级| 国产一区二区三区综合在线观看 | 中文欧美无线码| 香蕉精品网在线| 美女国产视频在线观看| 精品一区二区三卡| 成人漫画全彩无遮挡| 性高湖久久久久久久久免费观看| 亚洲精品国产av成人精品| 在线观看av片永久免费下载| 另类精品久久| 极品人妻少妇av视频| 日韩精品免费视频一区二区三区 | 亚洲欧美成人精品一区二区| 亚洲av电影在线观看一区二区三区| 人妻制服诱惑在线中文字幕| 亚洲综合精品二区| 婷婷色麻豆天堂久久| 精品人妻一区二区三区麻豆| 免费av中文字幕在线| 成人亚洲精品一区在线观看| 日本黄色片子视频| 国产精品免费大片| 美女脱内裤让男人舔精品视频| 美女视频免费永久观看网站| 欧美另类一区| av.在线天堂| 成年av动漫网址| 成人漫画全彩无遮挡| 我的老师免费观看完整版| 乱人伦中国视频| 97在线视频观看| 人妻制服诱惑在线中文字幕| 亚洲欧美清纯卡通| 国产精品国产三级专区第一集| 高清毛片免费看| 大码成人一级视频| 18+在线观看网站| 啦啦啦中文免费视频观看日本| 国产av码专区亚洲av| 狂野欧美白嫩少妇大欣赏|