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

    基于遺傳算法的可燃毒物設(shè)計(jì)優(yōu)化方法研究

    2021-08-02 03:02:54劉仕倡李滿倉(cāng)周冰燕王連杰陳義學(xué)
    原子能科學(xué)技術(shù) 2021年8期
    關(guān)鍵詞:燃耗毒物堆芯

    肖 鵬,王 健,劉仕倡,*,李滿倉(cāng),周冰燕,王連杰,陳義學(xué)

    (1.中國(guó)核動(dòng)力研究設(shè)計(jì)院 核反應(yīng)堆系統(tǒng)設(shè)計(jì)技術(shù)重點(diǎn)實(shí)驗(yàn)室,四川 成都 610213; 2.華北電力大學(xué) 核科學(xué)與工程學(xué)院,北京 102206)

    核反應(yīng)堆燃料組件設(shè)計(jì)是提高反應(yīng)堆燃料利用率、堆芯功率展平及反應(yīng)性控制的關(guān)鍵。其中,可燃毒物用于降低初始堆芯剩余反應(yīng)性和冷卻劑的硼濃度,保證反應(yīng)堆具有負(fù)的慢化劑溫度系數(shù)。可燃毒物組件與控制棒組件不同,為堆芯不動(dòng)部件??扇级疚锝M件是按物理設(shè)計(jì)分散布置于堆芯適當(dāng)位置的燃料組件中,有利于展平堆芯中子注量率的分布,可提高功率密度[1-2]。通過(guò)可燃毒物材料、含量、毒物棒數(shù)量等參數(shù),在滿足燃料和冷卻劑等反應(yīng)性系數(shù)在壽期內(nèi)為負(fù)、功率峰因子低于限值等約束條件下,實(shí)現(xiàn)降低初始反應(yīng)性、反應(yīng)性變化較平緩、堆芯壽期盡量長(zhǎng)、毒物在壽期內(nèi)消耗速率平緩、壽期末毒物殘留少、反應(yīng)性懲罰較小等優(yōu)化目標(biāo)[3]。對(duì)于核反應(yīng)堆燃料優(yōu)化設(shè)計(jì),由于是燃料棒、燃料組件、堆芯的多尺度問(wèn)題,同時(shí)包含燃料棒、毒物棒、控制棒等多個(gè)對(duì)象,而且涉及到不同溫度、密度等各種工況,甚至還是隨堆芯燃耗深度不斷變化的動(dòng)態(tài)問(wèn)題,因此核反應(yīng)堆燃料優(yōu)化設(shè)計(jì)具有很大的挑戰(zhàn)性。

    目前,常規(guī)的反應(yīng)堆燃料設(shè)計(jì)過(guò)程通常是設(shè)計(jì)者根據(jù)經(jīng)驗(yàn)給出一個(gè)初步方案,再針對(duì)此方案進(jìn)行堆芯物理計(jì)算分析以驗(yàn)證其設(shè)計(jì)的性能,如不滿足設(shè)計(jì)目標(biāo),則調(diào)整設(shè)計(jì)參數(shù)給出新方案并再次進(jìn)行計(jì)算,如此反復(fù)[4-5]。該方法能否得到優(yōu)化方案,很大程度上取決于設(shè)計(jì)者的經(jīng)驗(yàn),對(duì)于設(shè)計(jì)經(jīng)驗(yàn)缺乏人員往往耗時(shí)巨大且很難得到滿意解,即使是經(jīng)驗(yàn)豐富的設(shè)計(jì)者,進(jìn)行優(yōu)化迭代設(shè)計(jì)也非常耗時(shí)耗力。

    核反應(yīng)堆燃料優(yōu)化設(shè)計(jì)是一多輸入多目標(biāo)問(wèn)題。遺傳算法作為一種有效的隨機(jī)搜索算法,具有智能性、并行性、通用性等特點(diǎn)[6]。它在許多領(lǐng)域得到了應(yīng)用,取得了良好的效果。本文將遺傳算法和蒙特卡羅粒子輸運(yùn)方法相結(jié)合,研究遺傳算法在壓水堆燃料組件可燃毒物多目標(biāo)優(yōu)化設(shè)計(jì)中的理論模型和實(shí)現(xiàn)方法。

    1 壓水堆燃料組件中的可燃毒物設(shè)計(jì)

    本文研究的可燃毒物類型為可燃毒物及燃料混合在一起的整體型可燃毒物。優(yōu)化對(duì)象為13×13的壓水堆組件,其中包含3種棒:燃料棒、毒物棒和中心導(dǎo)管。組件結(jié)構(gòu)如圖1所示。

    圖1 壓水堆組件結(jié)構(gòu)模型 Fig.1 PWR assembly model

    1.1 可燃毒物優(yōu)化變量

    優(yōu)化變量包括:1) 毒物材料類型,包括6種氧化物,即Gd2O3、Er2O3、Sm2O3、Eu2O3、Dy2O3、HfO2(表1);2) 毒物含量,采用的范圍是毒物質(zhì)量占含可燃毒物燃料棒整體質(zhì)量(毒物+UO2)的1%~20%共20種質(zhì)量百分比情況;3) 排布類型(圖2);4) 軸向分層,在三維燃耗計(jì)算中將軸向分為5層,排列方式1表示整根棒全為可燃毒物,2表示兩端各有1層Zr層,中間3層為可燃毒物,3表示兩端4層為Zr層,只有中間1層為可燃毒物(圖3)。

    圖2 含可燃毒物燃料組件布置Fig.2 Fuel assembly layouts with burnable poison

    圖3 軸向分區(qū)布置Fig.3 Axial division

    表1 候選可燃毒物的密度和宏觀吸收截面Table 1 Theoretical density and macroscopic neutron absorption cross section of candidate burnable poison

    1.2 可燃毒物設(shè)計(jì)目標(biāo)

    從中子學(xué)角度,可燃毒物設(shè)計(jì)應(yīng)遵循3個(gè)原則:1) 壽期初,可燃毒物引入足夠大的負(fù)反應(yīng)性;2) 壽期內(nèi),可燃毒物緩慢釋放反應(yīng)性,控制功率分布,降低功率峰因子;3) 壽期末,可燃毒物的反應(yīng)性懲罰盡可能小。

    2 多目標(biāo)優(yōu)化遺傳算法

    遺傳算法是一種基于自然選擇的搜索算法[7],它通過(guò)繁殖將最適合的個(gè)體組合起來(lái),從而找到一個(gè)多元問(wèn)題的近似最優(yōu)解[8],并將潛在的新思想和變異性引入突變?nèi)后w中。該算法具有很強(qiáng)的并行性,因此本文將遺傳算法應(yīng)用于毒物選型優(yōu)化這種搜索空間非常復(fù)雜的多目標(biāo)優(yōu)化問(wèn)題[9]。

    針對(duì)多目標(biāo)優(yōu)化問(wèn)題[10],目前采用最多的是帕累托最優(yōu)解方法。帕累托前沿是以設(shè)計(jì)目標(biāo)為軸的優(yōu)化圖,由在圖上創(chuàng)建的線或平面顯示有沖突設(shè)計(jì)目標(biāo)之間的優(yōu)化權(quán)衡。在帕累托前沿,對(duì)于1個(gè)或多個(gè)目標(biāo),每個(gè)單獨(dú)的最優(yōu)值可能較任何其他最優(yōu)值都好,但不是所有的目標(biāo)都是最優(yōu)的。帕累托最優(yōu)的成員也被稱為“非占優(yōu)”[11]。

    本文的多目標(biāo)優(yōu)化采用NSGA-Ⅲ算法。NSGA-Ⅲ隨機(jī)生成含有N個(gè)個(gè)體的初始種群,算法進(jìn)行迭代直至終止條件滿足。在當(dāng)前種群的基礎(chǔ)上,通過(guò)隨機(jī)選擇產(chǎn)生子代種群。為從種群中選擇最好的N個(gè)解進(jìn)入下一代,首先利用基于帕累托前沿的非支配排序?qū)⑵浞譃槿舾刹煌姆侵鋵?。然后算法?gòu)建1個(gè)新的種群,構(gòu)建方法是逐次將各非支配層的解加入到新種群中,直至新種群大小等于N,或首次大于N,丟棄第N層之后的解。對(duì)于最后一層非支配層解的選取,NSGA-Ⅱ用的是基于擁擠距離的方法,而NSGA-Ⅲ用的是基于參考點(diǎn)的方法。擁擠距離度量并不適合求解3個(gè)及更多目標(biāo)的多目標(biāo)優(yōu)化問(wèn)題。因此NSGA-Ⅲ采用了新的選擇機(jī)制,該機(jī)制會(huì)通過(guò)所提供的參考點(diǎn),對(duì)新種群中的個(gè)體進(jìn)行更加系統(tǒng)的分析,以選擇最后的非支配層中部分解進(jìn)入新種群。因此本文選擇NSGA-Ⅲ作為多目標(biāo)優(yōu)化算法[12]。

    同時(shí),為提高計(jì)算效率,本文采用并行遺傳算法,將種群內(nèi)的個(gè)體分配到不同的計(jì)算核心。常規(guī)遺傳算法和并行遺傳算法的基本流程示于圖4。圖4中,COREn表示第n個(gè)計(jì)算核心。不同的計(jì)算核心處理不同的單個(gè)子種群,種群間互相獨(dú)立進(jìn)行進(jìn)化,種群間進(jìn)行個(gè)體遷移和種群競(jìng)爭(zhēng)。因此并行遺傳算法可顯著提高計(jì)算效率。

    圖4 并行遺傳算法基本流程圖Fig.4 Flow chart of parallel genetic algorithm

    3 可燃毒物優(yōu)化數(shù)學(xué)模型

    3.1 決策變量與目標(biāo)函數(shù)

    二維組件燃耗計(jì)算的輸入?yún)?shù)共3個(gè),故定義決策變量維數(shù)為3。第1個(gè)輸入變量為可燃毒物組件的排列方式(共7種),表示含可燃毒物燃料組件布置,具體排布見(jiàn)1.1節(jié);第2個(gè)輸入變量為毒物類型,包括6種不同材料;第3個(gè)輸入變量是可燃毒物的含量,為1%~20%共20種質(zhì)量占比。三維燃耗計(jì)算增加了1個(gè)決策變量,第4個(gè)決策變量是毒物棒的軸向分層,具體裝載情況見(jiàn)1.1節(jié)。4個(gè)變量以隨機(jī)數(shù)的形式拼接在1條染色體基因鏈上。由于各自的范圍不同,4個(gè)變量在基因鏈上的長(zhǎng)度亦不同。因?yàn)榛蜴準(zhǔn)怯啥M(jìn)制編碼的,所以4個(gè)輸入變量在染色體基因鏈上的長(zhǎng)度分別為3、3、7、2字節(jié),基因鏈總長(zhǎng)為15字節(jié)。

    燃耗計(jì)算目標(biāo)函數(shù)共5個(gè),分別表示為F1、F2、F3、F4、F5。F1為初始有效增殖因數(shù)(keff),F(xiàn)2為各燃耗步最大功率峰因子(ppf),F(xiàn)3為可燃毒物在燃耗中期的含量,F(xiàn)4為可燃毒物在燃耗末期的含量,F(xiàn)5為各燃耗步中有效增殖因數(shù)的最小值(表2)。其中,毒物含量為毒物氧化物中吸收截面最大的同位素剩余百分比。總?cè)己臅r(shí)長(zhǎng)為1 000 d,燃耗中期含量設(shè)為第500 d的毒物含量與初始毒物含量的比值;末期含量為第1 000 d毒物含量與初始毒物含量的比值。為降低初始反應(yīng)性,展平徑向功率分布,減少壽期末毒物殘留,分別設(shè)置優(yōu)化目標(biāo)為F1、F2和F4最小化。為使可燃毒物緩慢釋放反應(yīng)性,達(dá)到最大的循環(huán)長(zhǎng)度,分別設(shè)置優(yōu)化目標(biāo)為F3和F5最大化。

    表2 燃耗計(jì)算決策變量和目標(biāo)函數(shù)Table 2 Decision variable and objective function of burnup calculation

    3.2 約束條件

    組件燃耗計(jì)算中選用約束條件的方法對(duì)目標(biāo)結(jié)果加以限制。在燃耗計(jì)算中,設(shè)置了4個(gè)約束條件。首先使初始有效增殖因數(shù)大于1,即F1>1。二維燃耗計(jì)算采用最大功率峰因子F2<1.4的約束條件。而在三維燃耗計(jì)算中,因三維燃耗的加入使得功率不均勻程度增大,因此設(shè)置約束條件使各燃耗步最大功率峰因子F2<3.0。此外,燃耗中期要保證一定量的可燃毒物剩余,燃耗末期剩余的可燃毒物不能過(guò)多,故增加2個(gè)約束使燃耗中期剩余毒物的含量大于1%和燃耗末期剩余毒物的含量小于10%,即F3>0.01和F4<0.1。

    4 計(jì)算結(jié)果

    本文的遺傳算法采用Python的Geatpy遺傳算法庫(kù)[13]。Geatpy是一個(gè)高性能實(shí)用型進(jìn)化算法工具箱,提供許多已實(shí)現(xiàn)的進(jìn)化算法中各項(xiàng)重要操作的庫(kù)函數(shù),并提供高度模塊化、耦合度低的面向?qū)ο蟮倪M(jìn)化算法框架,利用“定義問(wèn)題類+調(diào)用算法模板”的模式進(jìn)行進(jìn)化優(yōu)化,可用于求解單目標(biāo)優(yōu)化、多目標(biāo)優(yōu)化、復(fù)雜約束優(yōu)化、組合優(yōu)化、混合編碼進(jìn)化優(yōu)化等。組件臨界計(jì)算及燃耗計(jì)算采用清華大學(xué)REAL團(tuán)隊(duì)研發(fā)的堆用蒙特卡羅程序RMC[14],RMC具備內(nèi)耦合的燃耗計(jì)算功能[15-17]。

    4.1 二維組件燃耗計(jì)算

    遺傳算法的種群數(shù)為50,迭代次數(shù)為50。RMC計(jì)算的源迭代采用的粒子數(shù)為2 000/代,非活躍代為5代,總代數(shù)為15,keff的統(tǒng)計(jì)標(biāo)準(zhǔn)差約為0.5%。比功率為30 W/gHM,燃料棒中235U的富集度為3%。不同棒數(shù)具體排列方式見(jiàn)1.1節(jié)。燃耗計(jì)算中采用的燃耗步數(shù)為6,步長(zhǎng)分別為100、100、150、150、250、250 d,共1 000 d。最大平均燃耗深度為30 MW·d/kgHM。不考慮帶軸向分層,毒物棒為整根全為燃料和毒物混合的單一材料。優(yōu)化方案及結(jié)果列于表3,最終的優(yōu)化方案中所用的材料均為Er2O3,因此Er2O3較其他材料更優(yōu)。該結(jié)果與通過(guò)人工搜索優(yōu)化得到的結(jié)論一致。二維優(yōu)化方案的keff隨燃耗時(shí)間的變化如圖5所示。

    圖5 二維組件各燃耗步keff隨燃耗時(shí)間的變化Fig.5 Change of keff with burnup time of two-dimensional assembly

    表3 二維組件燃耗計(jì)算優(yōu)化方案及結(jié)果Table 3 Optimization scheme and result of two-dimensional assembly burnup calculation

    4.2 三維組件燃耗計(jì)算

    考慮加入軸向分層后的燃耗計(jì)算,遺傳算法的種群數(shù)和迭代次數(shù)、蒙特卡羅計(jì)算的粒子數(shù)和代數(shù)均與二維算例相同。keff的統(tǒng)計(jì)標(biāo)準(zhǔn)差約為0.4%。

    優(yōu)化后得到40個(gè)非支配解,其中Gd2O3有13個(gè)、Er2O3有21個(gè)、Sm2O3有2個(gè)、Eu2O3有4個(gè)(表4)。與二維燃耗的只有1種材料Er2O3相比,通過(guò)軸向分層使得候選材料種類增加。各材料的keff隨燃耗時(shí)間的變化如圖6所示。

    圖6 三維組件各燃耗步keff隨時(shí)間的變化Fig.6 Change of keff with burnup time of three-dimensional assembly

    表4 三維組件燃耗計(jì)算優(yōu)化方案及結(jié)果Table 4 Optimization scheme and result of three-dimensional assembly burnup calculation

    三維燃耗ppf隨不同軸向分層的分布如圖7所示??煽闯觯謱臃绞?(即兩端各有1層Zr層、中間3層為可燃毒物)的燃耗ppf較分層方式1(即整根棒全為可燃毒物)的小??梢?jiàn)改變分層方式,以部分毒物代替整體全為可燃毒物的布置方式可減小ppf值。

    圖7 三維燃耗ppf隨軸向分層分布的變化Fig.7 Distribution of ppf in three-dimensional burnup with axial division

    5 結(jié)論

    本文將多目標(biāo)并行遺傳算法應(yīng)用于壓水堆組件毒物選型優(yōu)化,以反應(yīng)性控制、功率分布和不同時(shí)期燃耗剩余等為目標(biāo),對(duì)可燃毒物材料類型、含可燃毒物燃料棒排列方式、毒物含量、軸向分層等決策變量進(jìn)行了優(yōu)化。針對(duì)二維和三維燃耗計(jì)算,分別篩選了13和40種優(yōu)化方案。計(jì)算結(jié)果表明:Er2O3用作毒物的綜合效果最好;Gd2O3、Eu2O3和Sm2O3的應(yīng)用需結(jié)合堆芯方案開展進(jìn)一步研究;HfO2和Dy2O3不適合用作可燃毒物。該結(jié)果與通過(guò)人工搜索優(yōu)化得到的結(jié)論基本一致。同時(shí),三維軸向分層可為優(yōu)化提供更多可選的材料種類方案,以部分毒物的分層布置方式可減小ppf值。

    通過(guò)本文計(jì)算,初步驗(yàn)證了智能優(yōu)化算法在進(jìn)行毒物選型優(yōu)化時(shí),具有一定的協(xié)助工程師進(jìn)行優(yōu)化設(shè)計(jì)方案篩選的潛力。

    猜你喜歡
    燃耗毒物堆芯
    快跑!有毒物 絕密毒藥報(bào)告
    毒物——水銀
    應(yīng)用CDAG方法進(jìn)行EPR機(jī)組的嚴(yán)重事故堆芯損傷研究
    基于Hoogenboom基準(zhǔn)模型的SuperMC全堆芯計(jì)算能力校驗(yàn)
    核技術(shù)(2016年4期)2016-08-22 09:05:32
    基于切比雪夫有理逼近方法的蒙特卡羅燃耗計(jì)算研究與驗(yàn)證
    核技術(shù)(2016年4期)2016-08-22 09:05:28
    AMDIS在法醫(yī)毒物分析中的應(yīng)用
    IFBA/WABA 可燃毒物元件的燃耗特性分析
    壓水堆堆芯中應(yīng)用可燃毒物的兩個(gè)重要實(shí)驗(yàn)
    低價(jià)值控制棒中子吸收體材料燃耗相關(guān)數(shù)據(jù)的制作及驗(yàn)證研究
    基于SOP規(guī)程的大亞灣堆芯冷卻監(jiān)測(cè)系統(tǒng)改造
    国产精品免费一区二区三区在线| 搡老熟女国产l中国老女人| 免费在线观看日本一区| 看十八女毛片水多多多| 国产真实伦视频高清在线观看 | 我要搜黄色片| 在线看三级毛片| 美女被艹到高潮喷水动态| 在线十欧美十亚洲十日本专区| 嫩草影视91久久| 免费电影在线观看免费观看| 少妇裸体淫交视频免费看高清| 国产美女午夜福利| 亚洲在线自拍视频| 一进一出抽搐gif免费好疼| avwww免费| 亚洲av一区综合| 99热这里只有是精品在线观看 | 亚洲av第一区精品v没综合| 免费在线观看成人毛片| 欧美日本视频| 神马国产精品三级电影在线观看| 国产精品日韩av在线免费观看| 国产精华一区二区三区| 国产精品久久电影中文字幕| 在线观看66精品国产| 欧美bdsm另类| 久久亚洲真实| 欧美极品一区二区三区四区| 久久精品国产亚洲av涩爱 | 欧洲精品卡2卡3卡4卡5卡区| 三级男女做爰猛烈吃奶摸视频| 中文字幕人成人乱码亚洲影| 日日干狠狠操夜夜爽| 成人鲁丝片一二三区免费| 久久人妻av系列| 国产激情偷乱视频一区二区| 真人做人爱边吃奶动态| 又爽又黄a免费视频| 午夜久久久久精精品| 又爽又黄a免费视频| 成人美女网站在线观看视频| 日日夜夜操网爽| 亚洲国产日韩欧美精品在线观看| 欧美zozozo另类| 两个人视频免费观看高清| 欧美成人性av电影在线观看| 国模一区二区三区四区视频| 亚洲在线自拍视频| 日韩精品青青久久久久久| 亚洲内射少妇av| 国产淫片久久久久久久久 | 国产高清激情床上av| 欧美zozozo另类| 欧美极品一区二区三区四区| 亚洲专区国产一区二区| 国产伦精品一区二区三区视频9| av视频在线观看入口| 精品午夜福利视频在线观看一区| 欧美区成人在线视频| 香蕉av资源在线| 国产精品久久视频播放| 日韩欧美在线二视频| 成人国产综合亚洲| 亚洲av美国av| 久久性视频一级片| 国产欧美日韩精品一区二区| 国产男靠女视频免费网站| 淫秽高清视频在线观看| 国产欧美日韩一区二区精品| 色精品久久人妻99蜜桃| 欧美最黄视频在线播放免费| 日本一本二区三区精品| 亚洲三级黄色毛片| 又紧又爽又黄一区二区| 少妇被粗大猛烈的视频| 久久久久久久久久黄片| 久久精品综合一区二区三区| 88av欧美| 丰满人妻一区二区三区视频av| 大型黄色视频在线免费观看| 欧美日韩国产亚洲二区| 国产精品伦人一区二区| 久久人妻av系列| 成人一区二区视频在线观看| 国产精品电影一区二区三区| 久久久久免费精品人妻一区二区| 亚洲成av人片在线播放无| 天天一区二区日本电影三级| 欧美激情久久久久久爽电影| 在线a可以看的网站| netflix在线观看网站| 国内精品久久久久精免费| 国产69精品久久久久777片| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 亚洲中文字幕日韩| 99久久成人亚洲精品观看| 成年人黄色毛片网站| 少妇裸体淫交视频免费看高清| 久久这里只有精品中国| 午夜福利在线观看免费完整高清在 | 露出奶头的视频| 亚洲av.av天堂| av天堂在线播放| 三级国产精品欧美在线观看| 99久久99久久久精品蜜桃| 熟女人妻精品中文字幕| 男女视频在线观看网站免费| 国产蜜桃级精品一区二区三区| 日韩欧美三级三区| a级毛片a级免费在线| 亚洲无线在线观看| 国产成人影院久久av| 最近在线观看免费完整版| 免费在线观看亚洲国产| 日韩欧美在线二视频| 99热这里只有是精品50| 午夜视频国产福利| 国产精品99久久久久久久久| 九九在线视频观看精品| 国内少妇人妻偷人精品xxx网站| 亚洲成人久久性| 18禁黄网站禁片免费观看直播| 他把我摸到了高潮在线观看| 欧美性猛交黑人性爽| 亚洲aⅴ乱码一区二区在线播放| 好男人在线观看高清免费视频| 51国产日韩欧美| 校园春色视频在线观看| 三级国产精品欧美在线观看| 日本a在线网址| 欧洲精品卡2卡3卡4卡5卡区| 亚洲精品亚洲一区二区| 91久久精品国产一区二区成人| 日本在线视频免费播放| 亚洲av成人不卡在线观看播放网| 欧美黑人巨大hd| 观看免费一级毛片| 午夜福利在线在线| 欧美xxxx黑人xx丫x性爽| 国内精品久久久久精免费| 99国产极品粉嫩在线观看| 亚洲av第一区精品v没综合| 高清在线国产一区| 自拍偷自拍亚洲精品老妇| 亚洲黑人精品在线| 午夜免费男女啪啪视频观看 | 在线观看av片永久免费下载| 在线观看美女被高潮喷水网站 | 久久人人精品亚洲av| 91在线精品国自产拍蜜月| 日韩av在线大香蕉| 级片在线观看| 亚洲乱码一区二区免费版| 性插视频无遮挡在线免费观看| 欧美国产日韩亚洲一区| 国内精品美女久久久久久| 我要看日韩黄色一级片| 九色国产91popny在线| 亚洲专区中文字幕在线| 最近最新免费中文字幕在线| 国产亚洲精品综合一区在线观看| 黄色女人牲交| 亚洲人成伊人成综合网2020| 久久久久久久久久黄片| 99久久精品国产亚洲精品| 国产乱人伦免费视频| 美女高潮喷水抽搐中文字幕| 日本撒尿小便嘘嘘汇集6| 欧美激情久久久久久爽电影| 丰满乱子伦码专区| 成年女人毛片免费观看观看9| 国产精品1区2区在线观看.| 国产国拍精品亚洲av在线观看| aaaaa片日本免费| 久久久久国产精品人妻aⅴ院| 免费av毛片视频| 成人av在线播放网站| 久久久久久久久久成人| 麻豆国产97在线/欧美| 如何舔出高潮| 一区二区三区高清视频在线| 色播亚洲综合网| 国产色婷婷99| 国产精品av视频在线免费观看| av女优亚洲男人天堂| 欧美乱色亚洲激情| 黄色女人牲交| 久久久久九九精品影院| 国产综合懂色| 国产欧美日韩精品一区二区| 999久久久精品免费观看国产| 久久精品国产清高在天天线| 高清毛片免费观看视频网站| 精品99又大又爽又粗少妇毛片 | 激情在线观看视频在线高清| 日本与韩国留学比较| 一二三四社区在线视频社区8| 免费搜索国产男女视频| 麻豆av噜噜一区二区三区| 色吧在线观看| 3wmmmm亚洲av在线观看| 成年女人永久免费观看视频| 成人av在线播放网站| 国产亚洲欧美98| 啦啦啦观看免费观看视频高清| 在线天堂最新版资源| 久久这里只有精品中国| 色播亚洲综合网| 欧美丝袜亚洲另类 | 久久精品国产99精品国产亚洲性色| 欧美高清成人免费视频www| 中文字幕久久专区| h日本视频在线播放| avwww免费| 性插视频无遮挡在线免费观看| 国产成人av教育| 国产探花在线观看一区二区| 精品一区二区三区人妻视频| 床上黄色一级片| 51国产日韩欧美| 国产熟女xx| 久久国产乱子伦精品免费另类| 深爱激情五月婷婷| 热99re8久久精品国产| 窝窝影院91人妻| 国产高潮美女av| 床上黄色一级片| 国产成人av教育| 欧美潮喷喷水| 黄色一级大片看看| 两个人的视频大全免费| 精品一区二区免费观看| 免费搜索国产男女视频| 久久草成人影院| 婷婷精品国产亚洲av| 国产伦精品一区二区三区视频9| 97热精品久久久久久| 久久久精品大字幕| 国产精品一区二区性色av| 国产一区二区在线av高清观看| 成年人黄色毛片网站| 日日摸夜夜添夜夜添小说| 波多野结衣高清无吗| 舔av片在线| 亚洲精品乱码久久久v下载方式| 精品一区二区三区人妻视频| 国产黄色小视频在线观看| 又爽又黄a免费视频| 一进一出抽搐gif免费好疼| 国产精品日韩av在线免费观看| 国产视频内射| 亚洲国产日韩欧美精品在线观看| 午夜a级毛片| 欧美成人一区二区免费高清观看| 夜夜躁狠狠躁天天躁| 久久精品人妻少妇| 久久人人爽人人爽人人片va | 美女高潮喷水抽搐中文字幕| 岛国在线免费视频观看| 欧美黄色淫秽网站| 国产伦一二天堂av在线观看| 精品久久久久久,| 亚洲五月婷婷丁香| 欧美一级a爱片免费观看看| 在线观看66精品国产| www.色视频.com| 欧美成人a在线观看| 亚洲avbb在线观看| 9191精品国产免费久久| 男人舔奶头视频| 18禁在线播放成人免费| 亚洲欧美精品综合久久99| 亚洲经典国产精华液单 | 国产精品久久电影中文字幕| 国产精品亚洲一级av第二区| 精品午夜福利在线看| 99热这里只有是精品在线观看 | 欧美+日韩+精品| 国产精品综合久久久久久久免费| 亚洲av日韩精品久久久久久密| 99久久精品国产亚洲精品| 在线观看66精品国产| 九九热线精品视视频播放| 午夜福利成人在线免费观看| 亚洲一区二区三区色噜噜| 一夜夜www| 国产成人a区在线观看| 久久国产乱子伦精品免费另类| 亚洲五月婷婷丁香| 麻豆一二三区av精品| 国产主播在线观看一区二区| 精品欧美国产一区二区三| 在线播放无遮挡| 国产欧美日韩一区二区三| 一进一出抽搐动态| h日本视频在线播放| 99久久精品热视频| 无人区码免费观看不卡| av天堂中文字幕网| 淫秽高清视频在线观看| 欧洲精品卡2卡3卡4卡5卡区| 久久久久久九九精品二区国产| 久久6这里有精品| av在线老鸭窝| 亚洲国产精品999在线| 一本久久中文字幕| 欧美激情国产日韩精品一区| 亚洲乱码一区二区免费版| 毛片女人毛片| 午夜免费男女啪啪视频观看 | 亚洲,欧美,日韩| 亚洲美女视频黄频| 一级作爱视频免费观看| 人人妻人人看人人澡| 久久亚洲精品不卡| bbb黄色大片| 波野结衣二区三区在线| 一进一出好大好爽视频| 熟女人妻精品中文字幕| 欧美+日韩+精品| 欧美黑人欧美精品刺激| 一级av片app| 永久网站在线| 一区二区三区四区激情视频 | 精品人妻1区二区| 亚洲精品久久国产高清桃花| 看黄色毛片网站| 男女做爰动态图高潮gif福利片| av女优亚洲男人天堂| 免费人成在线观看视频色| 真人做人爱边吃奶动态| 久久欧美精品欧美久久欧美| 十八禁人妻一区二区| 午夜久久久久精精品| 午夜福利在线在线| 91在线精品国自产拍蜜月| 又黄又爽又刺激的免费视频.| 91字幕亚洲| 嫩草影院入口| 亚洲av中文字字幕乱码综合| 少妇被粗大猛烈的视频| 日韩亚洲欧美综合| 日韩国内少妇激情av| 久久人人爽人人爽人人片va | 真人做人爱边吃奶动态| 国内揄拍国产精品人妻在线| av女优亚洲男人天堂| 精品福利观看| 丰满乱子伦码专区| 国产久久久一区二区三区| 久久久久亚洲av毛片大全| 男女视频在线观看网站免费| 中文字幕高清在线视频| 久久这里只有精品中国| 三级毛片av免费| 999久久久精品免费观看国产| 亚洲欧美清纯卡通| 网址你懂的国产日韩在线| 51午夜福利影视在线观看| 特级一级黄色大片| 国产高清有码在线观看视频| 国产黄片美女视频| 十八禁人妻一区二区| 熟女电影av网| 精品一区二区三区视频在线| 久久精品91蜜桃| 毛片女人毛片| 日韩免费av在线播放| 麻豆国产av国片精品| 悠悠久久av| 国产精品嫩草影院av在线观看 | 亚洲真实伦在线观看| 亚洲人与动物交配视频| 2021天堂中文幕一二区在线观| 麻豆av噜噜一区二区三区| 亚洲自拍偷在线| 一级av片app| 亚洲美女搞黄在线观看 | 亚洲成av人片免费观看| 色5月婷婷丁香| 欧美日韩黄片免| 69人妻影院| www.熟女人妻精品国产| 久久人人精品亚洲av| 亚洲av第一区精品v没综合| 欧美最黄视频在线播放免费| 久久婷婷人人爽人人干人人爱| 国产精品伦人一区二区| 每晚都被弄得嗷嗷叫到高潮| 日韩高清综合在线| 国产伦在线观看视频一区| 久久久久久久午夜电影| 亚洲欧美日韩高清在线视频| 亚洲第一区二区三区不卡| av国产免费在线观看| 97超级碰碰碰精品色视频在线观看| 好男人电影高清在线观看| 色播亚洲综合网| 成人欧美大片| 亚洲精品色激情综合| 如何舔出高潮| 国内少妇人妻偷人精品xxx网站| 欧美黑人巨大hd| 欧美激情国产日韩精品一区| 亚洲18禁久久av| 亚洲男人的天堂狠狠| 国产私拍福利视频在线观看| 欧美3d第一页| 一进一出抽搐gif免费好疼| 国模一区二区三区四区视频| 搞女人的毛片| 精品久久久久久久久av| 中文字幕人妻熟人妻熟丝袜美| 亚洲精品日韩av片在线观看| 三级毛片av免费| 身体一侧抽搐| 性插视频无遮挡在线免费观看| 国产伦一二天堂av在线观看| 国产黄a三级三级三级人| 国产精品永久免费网站| 97人妻精品一区二区三区麻豆| 热99re8久久精品国产| 我的老师免费观看完整版| 18美女黄网站色大片免费观看| 久久精品国产亚洲av涩爱 | 麻豆国产av国片精品| 中文字幕av成人在线电影| 99久久精品国产亚洲精品| 成人三级黄色视频| 国产精品永久免费网站| 久久精品影院6| 日韩欧美在线乱码| 无人区码免费观看不卡| 美女大奶头视频| 日本三级黄在线观看| 欧美黄色片欧美黄色片| 久久久久久久久中文| 天美传媒精品一区二区| 国产老妇女一区| 麻豆国产av国片精品| 亚洲精品在线观看二区| 中出人妻视频一区二区| 精品一区二区免费观看| 成人欧美大片| 精品免费久久久久久久清纯| 午夜两性在线视频| 少妇人妻精品综合一区二区 | 老司机午夜十八禁免费视频| 少妇丰满av| 观看免费一级毛片| 午夜福利免费观看在线| 国产精品一区二区三区四区免费观看 | 国产三级中文精品| 老司机福利观看| 在线观看66精品国产| 国产成人影院久久av| 一进一出抽搐gif免费好疼| 亚洲成a人片在线一区二区| 男人舔女人下体高潮全视频| 天堂网av新在线| 亚洲av电影在线进入| 黄色配什么色好看| 午夜精品在线福利| 人人妻,人人澡人人爽秒播| 男人狂女人下面高潮的视频| 国产私拍福利视频在线观看| 午夜免费男女啪啪视频观看 | 18禁在线播放成人免费| av福利片在线观看| 国产美女午夜福利| 午夜福利欧美成人| 国产高清三级在线| 国产精品精品国产色婷婷| 757午夜福利合集在线观看| 俺也久久电影网| 亚洲 国产 在线| 欧美成人免费av一区二区三区| 国产精品免费一区二区三区在线| 真人做人爱边吃奶动态| 久久久久久久精品吃奶| av在线蜜桃| 成人鲁丝片一二三区免费| 丰满的人妻完整版| 中文字幕高清在线视频| 欧美黑人巨大hd| 99久久久亚洲精品蜜臀av| 亚洲成a人片在线一区二区| 婷婷六月久久综合丁香| 免费一级毛片在线播放高清视频| 一夜夜www| 在线a可以看的网站| 熟女人妻精品中文字幕| 少妇丰满av| 中文字幕久久专区| 69人妻影院| 精品一区二区三区视频在线观看免费| 久久久国产成人免费| 国产真实伦视频高清在线观看 | 永久网站在线| 性欧美人与动物交配| 高清日韩中文字幕在线| 赤兔流量卡办理| 免费搜索国产男女视频| 国内少妇人妻偷人精品xxx网站| 国产精品久久久久久久久免 | 国内久久婷婷六月综合欲色啪| www.999成人在线观看| 亚洲国产日韩欧美精品在线观看| 中国美女看黄片| 97碰自拍视频| 色尼玛亚洲综合影院| 哪里可以看免费的av片| 久久精品国产亚洲av涩爱 | 久久国产乱子免费精品| 国产 一区 欧美 日韩| 久久久久久大精品| 看十八女毛片水多多多| 亚洲人与动物交配视频| 国产精品av视频在线免费观看| 99久久精品一区二区三区| 最近中文字幕高清免费大全6 | 亚洲自拍偷在线| 天堂影院成人在线观看| 露出奶头的视频| av福利片在线观看| x7x7x7水蜜桃| 国产精品一及| 在线播放无遮挡| 久久久久国内视频| 国产精品影院久久| 乱码一卡2卡4卡精品| 露出奶头的视频| 一级黄色大片毛片| 亚洲成a人片在线一区二区| 久久精品国产清高在天天线| 国产亚洲精品av在线| 最后的刺客免费高清国语| 精品人妻1区二区| 国产大屁股一区二区在线视频| 90打野战视频偷拍视频| 亚洲专区国产一区二区| 国产视频内射| 男女床上黄色一级片免费看| 亚洲最大成人中文| 成人av在线播放网站| 亚洲片人在线观看| 欧美高清成人免费视频www| 亚洲精品久久国产高清桃花| 国产三级在线视频| av在线老鸭窝| 91在线精品国自产拍蜜月| 一区二区三区高清视频在线| 欧美高清性xxxxhd video| 亚洲av日韩精品久久久久久密| av天堂中文字幕网| x7x7x7水蜜桃| 99热只有精品国产| 熟妇人妻久久中文字幕3abv| 久9热在线精品视频| 一区福利在线观看| 日韩欧美在线乱码| 国产精华一区二区三区| 亚洲成人精品中文字幕电影| 精品久久久久久久久久免费视频| 岛国在线免费视频观看| 久久精品人妻少妇| 国产大屁股一区二区在线视频| 在线观看美女被高潮喷水网站 | 十八禁人妻一区二区| 国产午夜精品久久久久久一区二区三区 | 精品久久久久久久久av| 黄色日韩在线| 亚洲成人久久爱视频| av天堂在线播放| 亚洲欧美精品综合久久99| 日韩欧美一区二区三区在线观看| 国产色爽女视频免费观看| 99在线人妻在线中文字幕| 9191精品国产免费久久| 嫩草影院新地址| 国产精品1区2区在线观看.| 久久性视频一级片| 国产精品影院久久| 国产高清三级在线| 69av精品久久久久久| 此物有八面人人有两片| 国产精品一区二区免费欧美| 久久6这里有精品| 亚洲av成人不卡在线观看播放网| 级片在线观看| 精品免费久久久久久久清纯| 中文字幕av成人在线电影| 中文在线观看免费www的网站| 性色av乱码一区二区三区2| 国产爱豆传媒在线观看| 级片在线观看| а√天堂www在线а√下载| 99久久久亚洲精品蜜臀av| 亚洲第一电影网av| 国产黄色小视频在线观看| 欧美日本视频| 国产精品久久视频播放| 日韩欧美精品v在线| 丁香六月欧美| 桃红色精品国产亚洲av| 美女大奶头视频| 少妇的逼水好多| 久久九九热精品免费| 3wmmmm亚洲av在线观看| 99久久99久久久精品蜜桃| 亚洲欧美日韩高清专用| 一区福利在线观看| 成人永久免费在线观看视频| 日韩有码中文字幕|