亚洲免费av电影一区二区三区,日韩爱爱视频,51精品视频一区二区三区,91视频爱爱,日韩欧美在线播放视频,中文字幕少妇AV,亚洲电影中文字幕,久久久久亚洲av成人网址,久久综合视频网站,国产在线不卡免费播放

        ?

        最佳估算加不確定性分析方法在我國核安全審評中的應(yīng)用

        2016-03-13 03:26:47,,,,,
        核安全 2016年4期
        關(guān)鍵詞:包殼破口堆芯

        ,,,,,

        (環(huán)境保護(hù)部核與輻射安全中心,北京 100082)

        最佳估算加不確定性分析方法在我國核安全審評中的應(yīng)用

        靖劍平,賈斌,高新力,畢金生,孫微?,張春明

        (環(huán)境保護(hù)部核與輻射安全中心,北京 100082)

        目前核電廠安全分析用計算機程序多是基于保守方法開發(fā)的,給核電廠的設(shè)計和分析帶來了過量裕度,增加了核電廠優(yōu)化和改進(jìn)的難度,使用最佳估算加不確定性分析方法可以減少或消除這些不必要的限制。在AP1000和CAP1400的審評過程中,國家核安全局采用最佳估算加不確定性方法對大破口失水事故進(jìn)行了審查。本文介紹了四種最佳估算加不確定性分析方法,對不確定性的來源和不確定性統(tǒng)計方法進(jìn)行了論述?;贏STRUM方法,利用RELAP5程序?qū)P1000核電廠大破口失水事故進(jìn)行了獨立審核計算,經(jīng)59組抽樣計算后,最大的燃料包殼溫度值為1070℃,滿足驗收準(zhǔn)則要求。

        最佳估算;不確定性分析;核安全審評;AP1000;LOCA

        在過去的三十多年中,與核安全相關(guān)的現(xiàn)實環(huán)境一直在演變,人們對于事故進(jìn)程和物理現(xiàn)象方面的知識以及相應(yīng)的模擬能力也一直在進(jìn)步。與此同時,計算機的運算能力也取得了革命性的發(fā)展,這給計算分析者提供了更高的能力和更多的功能。美國核管會在1989年5月發(fā)布了RG1.157。依據(jù)導(dǎo)則規(guī)定,在確定論分析中,可以使用最佳估算方法,但同時要考慮不確定性。該導(dǎo)則的發(fā)布具有重要意義,使得安全分析方法向 “最佳估算方法加不確定性分析” (Best-Estimate Plus Uncertianty analyses,BEPU)的方向發(fā)展[1]。目前最佳估算方法加不確定性分析方法主要用于壓水堆核電廠大破口失水事故 (LB LOCA)分析。

        相比于國際上的研究,我國針對最佳估算方法加不確定性分析方法的研究起步晚,經(jīng)驗少。隨著AP1000技術(shù)的引進(jìn)以及我國自主知識產(chǎn)權(quán)的三代核電機組開發(fā),最佳估算方法加不確定性分析方法已用于AP1000和CAP1400大破口失水事故的分析。為了完成對AP1000和CAP1400的安全審評,國家核安全局組織開展了最佳估算方法加不確定性分析方法研究,并依據(jù)最佳估算加不確定性分析方法,利用RELAP5程序完成了AP1000大破口失水事故的分析,支撐了AP1000和CAP1400的安全審評工作,也為后續(xù)建立適用于我國的最佳估算方法加不確定性分析方法奠定了基礎(chǔ)。

        l 最佳估算方法加不確定性分析方法

        最佳估算要求估算盡可能地反應(yīng)真實的物理現(xiàn)象,既不偏悲觀,也不偏樂觀。國際原子能機構(gòu)認(rèn)為最佳估算事故分析應(yīng)滿足下述三個條件:(1)應(yīng)根據(jù)所選擇驗收準(zhǔn)則的要求,在事故分析中不引入過分的保守性; (2)使用最佳估算程序;(3)進(jìn)行不確定性分析[2-4]。

        在過去的近四十年中,國外相關(guān)研究機構(gòu)開展了大量的研究工作,并由此開發(fā)出多種最佳估算分析方法,包括 CSAU、ASTRUM、GRS、CIAU等方法[5]。

        l.l CSAU方法

        美國核管會 (NRC)在1989年頒布的RG l.157中,明確規(guī)定了最佳估算分析可采用經(jīng)驗關(guān)系式、模型、數(shù)據(jù)以及模型的評估方法和程序。并在附錄中提出了程序的比例模擬、適用性和不確定性分析 (Code Scaling,Applicability and Uncertainty,CSAU)方法,形成了一套用于開發(fā)最佳估算方法加不確定性分析方法的CSAU思想方法,并在Nuclear Engineering and Design雜志發(fā)表了六篇論文來對該方法及其在大破口失水事故分析中的應(yīng)用進(jìn)行描述[6-11]。

        CSAU的方法流程如圖1所示,分三大步14個小步。

        圖l CSAU方法流程圖Fig.l Flowchart of CSAU method

        l.2 ASTRUM方法

        西屋公司在CSAU方法的基礎(chǔ)上,開發(fā)了ASTRUM (Advanced Statistical Treatment of Uncertainty Method)方法,表示一種 “不確定自動統(tǒng)計處理”方法。在CQD方法基礎(chǔ)上,西屋公司更進(jìn)一步地完善了ASTRUM方法,并于2004年得到NRC的認(rèn)可[12]。該種方法遵循CSAU方法的步驟,并將CSAU方法中響應(yīng)面法替換為基于Wilk's公式的非參數(shù)統(tǒng)計方法。

        l.3 GRS方法

        德國核安全技術(shù)咨詢和安全分析中心(GRS)開發(fā)的最佳估算方法加不確定性分析方法稱為GRS方法。GRS方法主要著重于不確定性的統(tǒng)計,通過采用特定的不確定性研究工具,對參數(shù)值分布的重要性進(jìn)行分析。該方法的一個很明顯的優(yōu)勢是不人為地減少考慮不確定性的參數(shù)個數(shù),這一點與CSAU有顯著的區(qū)別。所有具有重要潛在影響的參數(shù)都在不確定性分析中得到評價,這是通過專家判斷無法達(dá)到的。采用非參數(shù)統(tǒng)計的方法,使得考慮不確定性影響的參數(shù)個數(shù)與計算次數(shù)相獨立。

        GRS方法的另外一個重要輸出是參數(shù)不確定性影響重要性排序,這些信息有助于通過進(jìn)一步把握狀態(tài)參數(shù)特性和提高程序模型減少輸出的不確定性。與其他的不確定性方法相區(qū)別的是:這種排序是程序輸出的結(jié)果,而不是預(yù)先的專家判斷[13]。

        l.4 CIAU方法

        CIAU (The Code with the capability of Internal Assessment of Uncertainty)方法是基于不確定性分析方法UMAE的一種最佳估算方法加不確定性分析方法,他們都是由意大利比薩大學(xué)F.D'Auria、A.Petruzzi等人開發(fā)的[14,15]。該方法是一種基于輸出不確定性方法,即注重程序計算結(jié)果與試驗結(jié)果的對比分析情況,而不十分關(guān)注單個輸入?yún)?shù)的不確定性傳遞。完整的CIAU主要有四個步驟,詳細(xì)流程如圖2所示。

        圖2 CIAU流程圖Fig.2 Flowchart of CIAU method

        2 不確定性分析方法

        采用最佳估算方法進(jìn)行事故模擬分析時,計算結(jié)果也帶有不確定性。實現(xiàn)結(jié)果的不確定性量化,對最佳估算分析有著重要意義。

        2.l 不確定性的來源

        不確定性分析應(yīng)全面的考慮各個環(huán)節(jié)的不確定性,包括程序模型不確定性、輸入?yún)?shù)不確定性以及基于整體性能信息的輸出結(jié)果不確定性。由于認(rèn)知和計算資源的限制,我們無法獲得所有熱工水力現(xiàn)象對輸出結(jié)果不確定性的影響程度,執(zhí)行過程中可能會出現(xiàn)遺漏問題。通常來講,不確定性的來源可以分為三部分[16]。

        第一部分來源于程序中的計算模型以及求解方法采用了近似,這是不確定性的根本來源。

        第二部分來源于建模時的簡化,這是不確定性的重要來源。

        第三部分來源于使用者的經(jīng)驗差異和業(yè)務(wù)水平等。

        表2總結(jié)了這三部分不確定性的來源并指明其量化的可能性。

        表2 主要不確定性來源及其量化可能性Table 2 The Main uncertainty source and quantization possibility

        2.2 不確定性統(tǒng)計方法

        對于不確定性計算,目前在最佳估算方法加不確定性分析方法中使用較多的統(tǒng)計方法有以下兩種。

        2.2.1 響應(yīng)面法

        響應(yīng)面法是將統(tǒng)計方法和數(shù)學(xué)方法相結(jié)合,通過對所關(guān)注的響應(yīng)進(jìn)行建模和分析,最終得到響應(yīng)值。由于響應(yīng)面法能夠較簡單地把確定性仿真和隨機仿真結(jié)合起來,因此在不確定性計算中得到了廣泛的應(yīng)用。

        2.2.2 非參數(shù)統(tǒng)計法

        在采用數(shù)理統(tǒng)計進(jìn)行不確定性的分析方面,主要有 “非參數(shù)統(tǒng)計”和 “參數(shù)統(tǒng)計”兩種方法。參數(shù)統(tǒng)計方法一般要求所分析的對象是總體分布函數(shù)已知 (如之前提及的響應(yīng)面的不確定性分析方法),而非參數(shù)統(tǒng)計方法對總體分布函數(shù)沒有明確要求,它與參數(shù)統(tǒng)計是相對應(yīng)的。與參數(shù)統(tǒng)計法相比,非參數(shù)統(tǒng)計法具有適用面廣、假設(shè)條件少、計算穩(wěn)健等特點。

        在不確定性分析中,程序的輸出結(jié)果一般為未知分布的隨機變量,因此更適合采用非參數(shù)統(tǒng)計方法進(jìn)行量化計算。非參數(shù)統(tǒng)計法通過抽樣計算得到總體在某概率 (γ)下一定置信水平 (β)的容許限,抽樣數(shù)量取決于概率和置信水平。

        Wilk's是最早開展容許限方面研究工作的學(xué)者之一,他提出了Wilk's公式 (見公式1和2)。該公式大大減少了樣本數(shù)量,并實現(xiàn)了抽樣樣本與輸入?yún)?shù)數(shù)量無關(guān)。公式1給出了單個輸出變量的單側(cè)容許限的Wilk's公式展開,公式2給出了單個輸出變量的雙側(cè)容許限的Wilk's公式展開。

        3 APl000核電廠大破口失水事故分析

        美國西屋公司在AP600的基礎(chǔ)上,研發(fā)了采用非能動技術(shù)的AP1000核電機組。與傳統(tǒng)壓水堆相比,該核電機組簡化了核電廠系統(tǒng),提高了核電廠的安全性和可靠性?;诖罅康膶嶒灁?shù)據(jù)和運行經(jīng)驗反饋,美國核管會于2001年發(fā)布了RELAP5/MOD 3.3程序,該程序在國內(nèi)外均得到廣泛應(yīng)用[17]。

        為了支持AP1000的安全審評,國家核安全局組織開展了獨立審核計算工作?;贏STRUM方法,利用RELAP5程序?qū)P1000核電廠大破口失水事故進(jìn)行了分析,抽樣程序采用的是DAKOTA程序。

        圖3 APl000系統(tǒng)節(jié)點模型Fig.3 Node model for APl000

        圖3為AP1000系統(tǒng)結(jié)構(gòu)節(jié)點圖。模型包含了反應(yīng)堆堆芯、一回路管道、蒸汽發(fā)生器、穩(wěn)壓器、主泵、堆芯補水箱、安注箱、非能動余熱排出系統(tǒng)以及邏輯控制系統(tǒng)等[18]。在大破口失水事故分析中,冷管段雙端剪切斷裂事故 (Double Ended Cold Leg,簡稱DECLG)是最惡劣的破口事故,并且當(dāng)破口發(fā)生在接有堆芯補水箱的環(huán)路上時,計算發(fā)生的燃料包殼溫度最高,本文即對該工況進(jìn)行模擬分析。

        圖4給出了冷卻劑系統(tǒng)壓力、堆芯補水箱注入流量、安注箱注入流量以及燃料包殼溫度(PCT)隨時間變化曲線。0秒時,冷管段雙端剪切斷裂事故發(fā)生,一回路冷卻劑通過破口噴出,冷卻劑系統(tǒng)的壓力迅速降低,并觸發(fā)停堆。隨著堆芯水位的降低,堆芯內(nèi)的冷卻劑被加熱汽化,此時流過堆芯的冷卻劑也由單相變成兩相。同時,由于堆芯流量和壓力的同時降低,燃料組件的冷卻狀況急劇惡化。在16s時,燃料包殼溫度達(dá)到第一個峰值,隨著功率的下降和堆芯補水箱的注入,燃料包殼溫度出現(xiàn)下降。由于在安注箱的水灌入到下腔室前,堆芯已處于裸露狀體,在衰變熱的加熱下,燃料組件包殼的溫度迅速上升,在97s時達(dá)到第二個峰值。隨著安注箱和堆芯補水箱的注入以及破口處流量的不斷減小,反應(yīng)堆冷卻劑系統(tǒng)裝量增加。隨著事故的發(fā)展,堆芯內(nèi)水位不斷增加,堆芯得到持續(xù)冷卻,堆芯下部空泡份額降低。當(dāng)堆芯水位達(dá)到足夠高時,堆芯得到充分冷卻,燃料包殼溫度降低。當(dāng)自動泄壓系統(tǒng)第四級閥門和內(nèi)置換料水箱隔離閥打開后,內(nèi)置換料水箱內(nèi)的水注入到一回路冷卻劑系統(tǒng)中,堆芯獲得了長期冷卻,燃料包殼溫度逐漸趨于穩(wěn)定,大破口失水事故結(jié)束。整個過程,燃料包殼溫度的最大值為1019.4℃。

        圖4 APl000大破口失水事故主要參數(shù)變化Fig.4 Primary parameter variations of the APl000 large-break LOCA accident

        在AP1000大破口失水事故分析的基礎(chǔ)上,本文重點對核電廠運行參數(shù)引入的不確定性進(jìn)行分析?;贏P1000的設(shè)計文件,選擇表3中的運行參數(shù)進(jìn)行不確定性分析,表中給出了各參數(shù)的平均值、最大值、最小值以及分布類型。由于輸出結(jié)果只關(guān)注燃料包殼溫度,基于Wilk's公式,共進(jìn)行了59組抽樣,取最大值作為不確定性分析后的燃料包殼溫度。

        表3 核電廠設(shè)計參數(shù)取值范圍Table 3 Range of APl000 design parameter valuesy

        圖5給出了前5個燃料包殼溫度最高的抽樣計算結(jié)果,由于輸入?yún)?shù)變化范圍不大,對整個事故進(jìn)程影響較小,所以不同抽樣工況的燃料包殼溫度變化規(guī)律基本接近,依據(jù)Wilk's理論,經(jīng)小于驗收準(zhǔn)則為1024℃。該值作為59組抽樣中的許用限值,如果該許用限值滿足驗收準(zhǔn)則,則說明該事故通過最佳估算及不確定性分析方法審評后,滿足驗收準(zhǔn)則。

        圖5 燃料包殼溫度隨時間變化Fig.5 PCT with time

        小于驗收準(zhǔn)則為1024℃。該值作為59組抽樣中的許用限值,如果該許用限值滿足驗收準(zhǔn)則,則說明該事故通過最佳估算及不確定性分析方法審評后,滿足驗收準(zhǔn)則。

        4 結(jié)論

        (1)隨著核安全要求的提高和發(fā)展,安全分析方法從保守的方法逐步轉(zhuǎn)為最佳估算方法加不確定性分析方法,典型的最佳估算方法加不確定性分析方法包括CSAU、GRS、ASTRUM、CIAE等方法,不確定性來源主要來自最佳估算程序模型、輸入?yún)?shù)和用戶因素。

        (2)基于ASTRUM方法,利用RELAP5程序?qū)P1000核電廠大破口失水事故進(jìn)行了獨立審核計算,重點考慮對核電廠運行參數(shù)引入的不確定性,經(jīng)59組抽樣計算后,最大的燃料包殼溫度值為1070℃,滿足驗收準(zhǔn)則要求。

        (3)將最佳估算方法應(yīng)用于我國的核電廠執(zhí)照申請級安全分析,不僅可以提升我國核電廠安全分析水平,而且可以有效提升核電廠的經(jīng)濟性,從而促進(jìn)我國核電行業(yè)和核安全分析技術(shù)發(fā)展。因此我國核電領(lǐng)域應(yīng)對最佳估算方法加不確定性分析方法給予高度重視,并組織力量開展深入研究,特別是核安全監(jiān)管部門更應(yīng)全面理解并掌握其中的關(guān)鍵技術(shù)內(nèi)容,參考目前廣泛應(yīng)用的最佳估算方法加不確定性分析方法,逐步建立適用于我國最佳估算方法加不確定性分析的技術(shù)文件體系,為大破口失水事故分析的審評和獨立審核計算提供指導(dǎo),也為國內(nèi)自主化軟件的開發(fā)、驗證和審查提供支持。

        [1]路璐.認(rèn)證級LOCA分析工具開發(fā)中相關(guān)模型修改及驗證[D].上海:上海交通大學(xué),2010.

        [2]IAEA.Implications of power prates on safe margins of nuclear power plant[R],Wienna:IAEA,2004.

        [3]靖劍平,孫微,賈斌,等.最佳估算程序不確定性分析方法研究[C].中國核學(xué)會 2013年學(xué)術(shù)年會論文集第 3冊,2013.

        [4]林誠格,劉志弢,趙瑞昌.壓水堆失水事故最佳估算方法研究[J].核安全,2010,30(1):1-12.

        [5]倪超.AP1000核電廠大破口失水事故最佳估算分析建模與不確定性研究[D].上海:上海交通大學(xué),2011.

        [6]B.E.Boyack,et al.Quantifying reactor safety margins-Part 1:An overview of the CSAU evaluation methodology[J].Nuclear Engineering and Design,1990,119(1):1-15.

        [7]G.E.Wilson,et al.Quantifying reactor safety margins-Part 2:Characterization of important contributors to uncertainty[J].Nuclear Engineering and Design,1990,119(1):17-31.

        [8]W.Wulff,et al.Quantifying reactor safety margins-Part 3:Assessment and ranging of parameters[J].Nuclear Engineering and Design,1990,119(1):33-65.

        [9]G.S.Lellouche,S.Levy,et al.Quantifying reactor safety margins-Part 4:Uncertainty evaluation of LB LOCA analysis based on TRAC-PF1/MOD1[J].Nuclear Engineering and Design,1990,119(1):67-95.

        [10]N.Zuber,G.E.Wilson,et al.Quantifying reactor safety margins-Part 5:Evaluation of scale-upcapabilities of best estimate codes[J].Nuclear Engineering and Design,1990,119(1):97-107.

        [11]I.Catton,R.B.Duffey,R.A.Shaw,et al.Quantifying reactor safety margins-Part 6:A physically-based method of estimating LB LOCA PCT[J].Nuclear Engineering and Design,1990,119(1):109-117.

        [12]Westighouse Non-proprietary class 3.Realistic Large-break LOCA evaluation methodology using the automated statistical treatment of uncertainty method(ASTRUM)[R],Washington:Westighouse,2005.

        [13]Horst Glaeser.GRS Method for Uncertainty and Sensitivity E-valuation of Code Results and Applications[R].Science and Technology of Nuclear Installations,2008,vol.2008,7 pages,Article ID 798901,doi:10.1155/2008/798901.

        [14]A.Petruzzi,F(xiàn).D'Auria.CIAU method for uncertainty evaluation for system thermal-hydraulic code calculations[C]//Proceedings of ICAPP.Nice,F(xiàn)rance:ReportSociete Francaise d'Energie Nucleaire,2011:2045-2054.

        [15]F.D'Auria.Approach and methods to evaluate the uncertainty in system thermalhydrulic calculations[C]//International conference on simulation methods in nuclear engineering,Ontario,Canada:ReportCanadian Nuclear Society,2004:12-17.

        [16]T.Wickett,et al.Report of the uncertainty methods study for advanced best estimate thermal hydraulic code applications[R].Washington:NEA/CSNI/OECD,1998.

        [17]靖劍平,張春明,陳妍,等.淺談核電領(lǐng)域中的熱工水力分析程序[J].核安全,2012,33(3):70-74.

        [18]靖劍平,喬雪冬,賈斌,等.基于RELAP5程序的AP1000典型事故瞬態(tài)特性研究[J].原子能科學(xué)與技術(shù),2015,49(4):646-653.

        Study on the Nuclear Safety Review through Best-Estimate Plus Uncertainty Analyses Methods

        JING Jianping,JIA Bin,GAO Xinli,BI Jinsheng,SUN Wei?,ZHANG Chunming
        (Nuclear and Radiation Safety Center,MEP,Beijing 100082,China)

        The methods that used in most conventional computer codes are conservative,which will lead to unnecessary excessive margins and more difficulty for nuclear power plants design and improvement.U-sing best-estimate plus uncertainty analysis method may reduce or eliminate these unnecessary limits.In the review process of AP1000 and CAP1400,National Nuclear Safety Agency has used the best-estimate method for large-break LOCA accident.In this paper the four best-estimate plus uncertainty analysis methods are introduced,the source of uncertainty and uncertainty statistic methods are described.Then the independent audit calculation of AP1000 NPP large-break LOCA accident is carried out with RELAP5 code based on ASTRUM method.Through calculations of 59 groups of samples,the maximum PCT is 1 070℃,and it can meet the acceptance criteria.

        Best-estimate;Uncertainty analyses;Nuclear safety review;AP1000;LOCA

        TL99

        :A

        :1672-5360(2016)04-0011-07

        2016-05-13

        2016-08-16

        大型先進(jìn)壓水堆及高溫氣冷堆核電站國家科技重大專項 “CAP1400安全審評技術(shù)及獨立驗證試驗” (2011ZX06002-010)和 “CAP1400安全審評關(guān)鍵技術(shù)研究”(2013ZX06002001)

        靖劍平 (1983),男,漢,黑龍江綏棱人,高級工程師,博士,現(xiàn)主要從事反應(yīng)堆熱工水力及安全分析方面的研究工作

        ?通訊作者:孫 微,E-mail:sunwei@chinansc.cn

        猜你喜歡
        包殼破口堆芯
        LOCA事故下碳化硅復(fù)合包殼失效概率計算
        核技術(shù)(2023年9期)2023-09-21 09:21:32
        華龍一號蒸汽發(fā)生器傳熱管6mm破口事故放射性后果分析
        核安全(2022年2期)2022-05-05 06:55:38
        基于“華龍一號”大破口事故先進(jìn)安注箱研究
        碳化硅復(fù)合包殼穩(wěn)態(tài)應(yīng)力與失效概率分析
        耐事故包殼中子經(jīng)濟性分析*
        中國核電(2021年2期)2021-06-04 06:23:16
        破口
        應(yīng)用CDAG方法進(jìn)行EPR機組的嚴(yán)重事故堆芯損傷研究
        基于Hoogenboom基準(zhǔn)模型的SuperMC全堆芯計算能力校驗
        核技術(shù)(2016年4期)2016-08-22 09:05:32
        改善研究堆用鋁合金包殼抗腐蝕性能的研究
        壓水堆堆芯中應(yīng)用可燃毒物的兩個重要實驗
        色橹橹欧美在线观看视频高清| 百合av一区二区三区| 人人爽亚洲aⅴ人人爽av人人片| 亚洲精品美女久久久久网站| 久久精品国产亚洲av四区| 国产中文三级全黄| 无码av免费一区二区三区| 亚洲成在人线电影天堂色| 中文字幕亚洲中文第一| 日韩国产人妻一区二区三区| 蜜臀av无码精品人妻色欲| 日韩av一区二区三区四区av| 亚洲国产av高清一区二区三区| 黄片视频免费观看蜜桃| 又黄又硬又湿又刺激视频免费| 色欲aⅴ亚洲情无码av蜜桃| 日韩亚洲国产av自拍| 久久精品国产自产对白一区| 丰满人妻一区二区三区免费视频| 高清偷自拍第1页| 免费国产一级特黄aa大片在线| 久久爱91精品国产一区| 日韩亚洲一区二区三区四区| 国产免费人成视频在线观看| 日韩一区二区肥| 蜜桃激情视频一区二区| 欧美大胆性生话| 国产精品区一区第一页| 亚洲精品白浆高清久久| 在线观看 国产一区二区三区| 乱子伦一区二区三区| 熟妇人妻AV中文字幕老熟妇 | 亚洲欧美日韩国产一区二区精品| 国产精品亚洲av国产| 亚洲黄片av在线播放| 国产一区二区三精品久久久无广告 | 99精品热这里只有精品| 大屁股少妇一区二区无码| 精品一区二区av在线| 亚洲裸男gv网站| 久久亚洲道色宗和久久|