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

    核設(shè)施退役過程中的輻射場重構(gòu)與拆除路徑優(yōu)化

    2017-10-13 09:21:58宋英明梁燁葉凱萱張震宇張秋楠王寧顏佳偉朱志超鄒樹梁
    核技術(shù) 2017年5期
    關(guān)鍵詞:核設(shè)施器件遺傳算法

    宋英明 梁燁 葉凱萱 張震宇 張秋楠 王寧 顏佳偉 朱志超 鄒樹梁

    ?

    核設(shè)施退役過程中的輻射場重構(gòu)與拆除路徑優(yōu)化

    宋英明1梁燁1葉凱萱1張震宇1張秋楠1王寧2顏佳偉1朱志超1鄒樹梁1

    1(南華大學(xué)核科學(xué)技術(shù)學(xué)院 衡陽 421001) 2(中國核動(dòng)力研究設(shè)計(jì)院信息中心 成都610213)

    輻射場分布和路徑選擇是核設(shè)施退役過程中影響工作人員吸收劑量的主要因素。針對(duì)核設(shè)施退役過程中多源項(xiàng)的三維輻射場,采用點(diǎn)核積分的方法進(jìn)行重構(gòu),重構(gòu)結(jié)果與蒙特卡羅程序計(jì)算結(jié)果吻合很好,驗(yàn)證了點(diǎn)核積分法的可行性。將核退役拆除路徑問題抽象為一種類旅行商問題數(shù)學(xué)模型,構(gòu)造不同拆除路徑下所受外照射劑量對(duì)應(yīng)的劑量矩陣,根據(jù)輻射防護(hù)ALARA (As Low As Reasonably Achievable)原則,利用遺傳算法進(jìn)行尋優(yōu)。對(duì)于多源項(xiàng)退役拆除實(shí)例,計(jì)算給出了最優(yōu)化拆除路徑和三維可視化顯示,并對(duì)優(yōu)化效果進(jìn)行了討論。設(shè)計(jì)開發(fā)了面向核設(shè)施退役過程的輻射場重構(gòu)與拆除路徑優(yōu)化功能軟件。

    核設(shè)施退役,輻射場重構(gòu),拆除路徑優(yōu)化,功能軟件

    核設(shè)施退役是一項(xiàng)周期長、涉及面廣、投資高的系統(tǒng)工程,基于安全性、經(jīng)濟(jì)性、合理性的考慮,在實(shí)施各項(xiàng)核設(shè)施退役工作前,必須借助計(jì)算機(jī)仿真手段對(duì)退役技術(shù)方案和關(guān)鍵步驟進(jìn)行驗(yàn)證、評(píng)估和演練,優(yōu)化拆除工藝,以降低對(duì)人員和環(huán)境的 危害。

    國際上對(duì)核設(shè)施退役仿真系統(tǒng)進(jìn)行了廣泛和深入的研究。法國原子能委員會(huì)(French Atomic Energy Commission, CEA)[1]開發(fā)了基于虛擬現(xiàn)實(shí)技術(shù)框架與點(diǎn)核積分方法的NARVEOS系統(tǒng)。該軟件支持情景設(shè)計(jì)、快速劑量計(jì)算,是適用于核輻射評(píng)估的ALARA (As Low As Reasonably Achievable)分析工具。經(jīng)濟(jì)合作與發(fā)展組織核能局(Nuclear Energy Agency, Organization for Economic Co-operation and Development, OECD/NEA)、日本原子力科學(xué)研究機(jī)構(gòu)(Japan Atomic Energy Agency, JAEA)和哈爾登(Halden)反應(yīng)堆工程中心[2]共同研發(fā)了一套人員劑量評(píng)估系統(tǒng)VRDose,可實(shí)現(xiàn)評(píng)估劑量、優(yōu)化退役計(jì)劃等功能。韓國原子能研究所和忠南大學(xué)[3]聯(lián)合對(duì)韓國KRR-1&2反應(yīng)堆退役開展工人輻射劑量計(jì)算機(jī)模擬,能實(shí)現(xiàn)基于MCNP-4C[4]輻射輸運(yùn)程序下的人的輻射劑量計(jì)算。隨著我國核能事業(yè)的高速發(fā)展,對(duì)核設(shè)施退役技術(shù)研究與工程準(zhǔn)備工作的要求更加體系化、專業(yè)化、標(biāo)準(zhǔn)化。目前,國內(nèi)已有部分單位在輻射場可視化軟件開發(fā)和核設(shè)施虛擬退役等方面開展了相關(guān)研究,取得了較好的成果[5?9],但對(duì)于核設(shè)施退役中拆除路徑最優(yōu)化方法的研究相對(duì)較少且缺乏針對(duì)性。本文驗(yàn)證并采用點(diǎn)核積分法重構(gòu)核設(shè)施退役過程中多源項(xiàng)的三維輻射場,然后利用類旅行商問題(Traveling Salesman Problem, TSP)數(shù)學(xué)模型求解核退役拆除路徑問題,構(gòu)造不同拆除路徑下所受外照射劑量矩陣,通過遺傳算法尋找最優(yōu)路徑。設(shè)計(jì)開發(fā)了軟件實(shí)現(xiàn)相應(yīng)功能并進(jìn)行了實(shí)例分析。

    1 核退役設(shè)施三維輻射場重構(gòu)

    定量和可視化描述核退役過程中工作人員所受輻射照射,需要構(gòu)建核設(shè)施的三維輻射場。由于核退役設(shè)施一般規(guī)模較大、源項(xiàng)較多,利用MCNP軟件求解核設(shè)施內(nèi)部輻射場會(huì)導(dǎo)致運(yùn)算時(shí)間很長,且難以滿足誤差要求,因此發(fā)展基于點(diǎn)核積分法的核退役設(shè)施三維輻射場求解方法。

    1.1 點(diǎn)核積分重構(gòu)三維輻射場

    點(diǎn)核積分計(jì)算核設(shè)施退役場所內(nèi)輻射場分布的方法如下[10]:

    1.2 蒙特卡羅程序校驗(yàn)點(diǎn)核積分結(jié)果

    利用MCNP5驗(yàn)證簡單源項(xiàng)情況下點(diǎn)核積分的結(jié)果,為簡化計(jì)算,將結(jié)果歸一化處理,如圖1所示。其中,圖1(a)是退役設(shè)施內(nèi)為一點(diǎn)源的結(jié)果,圖1(b)是核設(shè)施內(nèi)為一立方體源(10cm×10cm× 10cm)的結(jié)果。從圖1可以發(fā)現(xiàn),點(diǎn)源和體源情況下,點(diǎn)核積分結(jié)果都與MCNP計(jì)算結(jié)果吻合很好。這可在相同計(jì)算精度下,大幅節(jié)省最優(yōu)退役拆除路徑的搜尋時(shí)間。

    2 基于TSP的核退役拆除路徑優(yōu)化模型

    2.1 核退役拆除路徑問題的簡述

    TSP是著名的數(shù)學(xué)圖論問題,應(yīng)用背景廣泛,問題描述如下:旅行商從某城市出發(fā),遍歷各城市后返回到出發(fā)城市,尋找旅行商所走總距離最小的路線。

    在核設(shè)施退役模型中,設(shè)施內(nèi)有個(gè)器件,則拆除所有器件的路徑有!條,路徑數(shù)目隨增加而急劇上升,最優(yōu)路徑的求解難度增加。為得到最優(yōu)退役拆除路徑,將核設(shè)施退役拆除路徑問題假設(shè)成類似TSP問題,特點(diǎn)在于拆除路徑起點(diǎn)可不同、可不為閉合回路,即:

    1) 退役過程中需拆除的器件作為TSP問題中的城市;

    2) 退役過程中的優(yōu)化目標(biāo)為工作人員所受總劑量,作為TSP問題中的總路程。由于在設(shè)施中行走所需的時(shí)間遠(yuǎn)小于拆除過程所需的時(shí)間,忽略工作人員行走中所受劑量。

    2.2 不同路徑下總劑量的求解

    建立時(shí)間矩陣:

    劑量率矩陣和時(shí)間矩陣相乘,得劑量矩陣,矩陣各元素求和即為拆除人員受到的總劑量()。

    2.3 器件單元化的推廣

    實(shí)際操作中,一個(gè)器件可能分多次拆除。該情況下,設(shè)定器件單元數(shù),即一個(gè)器件平均分為個(gè)單元,各單元的活度和產(chǎn)生的劑量率為器件的1/,拆除各單元所需的時(shí)間為拆除器件總時(shí)間的1/設(shè)同一器件中各單元位置相同,則同一器件不同單元輻射效果相同。單元總數(shù)為′個(gè),拆除人員需遍歷所有單元。

    (8)

    3 遺傳算法尋優(yōu)計(jì)算與分析

    旅行商問題有多種解法,目前較為常用的算法有神經(jīng)網(wǎng)絡(luò)法和遺傳算法。遺傳算法[11]是最近幾年發(fā)展起來的全局優(yōu)化算法,借用生物遺傳學(xué)觀點(diǎn),通過選擇、遺傳、變異和免疫等作用機(jī)制,使每個(gè)個(gè)體的適應(yīng)性提高。由于全局搜索的特性,遺傳算法解決TSP問題有特定優(yōu)勢(shì)。故采用遺傳算法求解上述問題,并通過圖形用戶界面(Graphical User Interface, GUI)設(shè)計(jì)了求解本問題的可視化軟件。

    3.1 算法流程

    利用遺傳算法求解核退役拆除路徑優(yōu)化問題的步驟如下:

    1) 種群初始化

    由于核設(shè)施問題假設(shè)成多目標(biāo)旅行商問題,個(gè)體編碼方法采用實(shí)數(shù)編碼。在核設(shè)施退役中,將器件組順序?qū)崝?shù)編碼為1?(設(shè)施內(nèi)的總器件數(shù))或1?(設(shè)施內(nèi)的總單元數(shù))的隨機(jī)排列,初始化的參數(shù)有種群個(gè)數(shù)、染色體基因個(gè)數(shù)、迭代次數(shù)c、交叉概率c以及變異概率m。

    2) 建立適應(yīng)度函數(shù)

    拆除某一器件時(shí)所受劑量率由點(diǎn)核積分已知,器件組隨機(jī)全排列對(duì)應(yīng)的總劑量可求,總劑量越小,種群的適應(yīng)度函數(shù)越好。將總劑量的倒數(shù)作為適應(yīng)度函數(shù),由求min()變?yōu)榍髆ax[1/()]。

    3) 選擇操作

    用輪盤賭法和兩兩競爭的方法進(jìn)行選擇操作,采用基于適應(yīng)度比例的選擇策略,即適應(yīng)度越好的個(gè)體被選擇的概率越大,并保存所選適應(yīng)度最高的個(gè)體。

    4) 交叉操作

    利用部分映射交叉(Partial Mapped Crossover, PMX)算法進(jìn)行交叉操作,先對(duì)父代進(jìn)行常規(guī)的兩點(diǎn)交叉,再根據(jù)交叉區(qū)域內(nèi)各基因值之間的映射關(guān)系來修改交叉區(qū)域之外的各個(gè)基因值,該算法可以較多地產(chǎn)生父代遍歷路徑所沒有的部分新路徑,收斂速度較低,但其全局搜索能力較強(qiáng)。

    5) 變異操作

    對(duì)于變異操作,隨機(jī)選取個(gè)體,同時(shí)隨機(jī)選取個(gè)體的兩個(gè)基因進(jìn)行交換以實(shí)現(xiàn)變異操作。

    算法整體流程如圖2所示。

    圖2 求解核設(shè)施退役拆除路徑優(yōu)化問題的算法整體流程

    3.2 求解實(shí)例

    當(dāng)核設(shè)施內(nèi)源項(xiàng)個(gè)數(shù)為1時(shí),不同的器件單元數(shù)必然對(duì)應(yīng)同一優(yōu)化路徑。經(jīng)程序驗(yàn)證發(fā)現(xiàn),源項(xiàng)較少時(shí),不同的器件單元數(shù)趨于對(duì)應(yīng)同一優(yōu)化路徑,即單元數(shù)為1對(duì)應(yīng)的路徑。故源項(xiàng)較少情況下,器件可分多次拆除無明顯優(yōu)勢(shì);源項(xiàng)較多時(shí),分多次拆除同一器件的路徑較為復(fù)雜,為便于觀看,建立一個(gè)10cm×10cm×10cm的廠房,隨機(jī)放置30個(gè)正方體形狀的放射性器件,器件的尺寸在(1?8)×10?3m3內(nèi)隨機(jī)分布,器件的活度在(1.85?5.5)×1010內(nèi)隨機(jī)分布,假設(shè)拆除時(shí)間與拆除器件的體積成正比:t=/(其中:為體積,m3;為比例系數(shù),以下取=3為例)。

    各器件單元數(shù)=1時(shí),利用MATLAB實(shí)現(xiàn)遺傳算法求解基于旅行商問題下的最優(yōu)退役路線如圖3(a)所示,其中初始種群為100個(gè),迭代次數(shù)為1000次。由圖3(b)可知,該例中拆除人員的總劑量在迭代200次左右達(dá)到收斂,最優(yōu)路線對(duì)應(yīng)的受照總劑量為。

    單元數(shù)設(shè)為5時(shí),拆除人員可分多次對(duì)某一器件進(jìn)行拆除,最優(yōu)退役路徑如圖4(a)所示。由圖4(b)可知,該例中拆除人員的總劑量在迭代900次左右達(dá)到收斂,最優(yōu)路徑對(duì)應(yīng)受照總劑量為。

    由計(jì)算結(jié)果可知,最優(yōu)路徑可能存在于分多次拆除同一器件情況下。另外,針對(duì)核設(shè)施拆除過程中可能造成的二次污染,可通過劑量估算后,再次計(jì)算輻射場,利用功能軟件進(jìn)行二次尋優(yōu)。

    4 功能軟件開發(fā)

    利用MATLAB實(shí)現(xiàn)了基于遺傳算法求解旅行商問題下的最優(yōu)退役路線,并通過GUI設(shè)計(jì)了面向核設(shè)施退役過程的輻射場重構(gòu)與拆除路徑優(yōu)化功能軟件。

    圖5為新建廠房數(shù)據(jù)界面,在該界面可進(jìn)行三維情景設(shè)計(jì),包括放射性器件數(shù)據(jù)如器件位置、形狀、尺寸、活度、拆除所需時(shí)間及廠房相關(guān)數(shù)據(jù)并保存,廠房及任意房間的布局可視。

    圖6為最優(yōu)路徑求解界面,可在該界面更換保存的廠房數(shù)據(jù),調(diào)整器件單元數(shù)、迭代次數(shù)、種群數(shù)量等參數(shù),并顯示遺傳算法求解最優(yōu)路徑的進(jìn)程,所得結(jié)果可保存為文本文件或Excel文件。

    圖5 新建廠房數(shù)據(jù)界面

    圖6 求解最優(yōu)路徑界面

    圖7為最終的結(jié)果,給出了迭代次數(shù)與拆除人員受到總劑量的曲線,可以得到收斂后所受的最小劑量值。圖7還給出了最優(yōu)拆除路徑,最優(yōu)路徑窗口中按最優(yōu)路徑依次顯示了每次要拆除的器件以及所要拆除的百分比。

    圖7 最終結(jié)果

    圖8為最優(yōu)拆除路徑的動(dòng)態(tài)顯示,顯示速度可調(diào)。正被拆除的器件的顏色將不停變化,尺寸將不斷變小。該界面以一條移動(dòng)的短線顯示拆除人員在器件間的運(yùn)動(dòng)過程。

    圖8 最優(yōu)路徑的動(dòng)態(tài)顯示圖

    5 結(jié)語

    針對(duì)核設(shè)施退役過程中多源項(xiàng)的三維輻射場,發(fā)展了基于點(diǎn)核積分的輻射場快速重構(gòu)方法,并利用蒙特卡羅程序驗(yàn)證。核退役拆除路徑問題抽象為類TSP問題數(shù)學(xué)模型,構(gòu)造不同拆除路徑下所受外照射的劑量矩陣,根據(jù)輻射防護(hù)ALARA原則,利用遺傳算法尋優(yōu)計(jì)算。計(jì)算給出了多源項(xiàng)拆除實(shí)例最優(yōu)化拆除路徑和三維顯示,并對(duì)優(yōu)化效果進(jìn)行了討論?;谝陨戏椒ㄔO(shè)計(jì)開發(fā)了面向核設(shè)施退役過程的輻射場重構(gòu)與拆除路徑優(yōu)化功能軟件。據(jù)此可進(jìn)一步研究多旅行商問題下的核退役拆除路徑優(yōu)化問題,并考慮核設(shè)施退役過程中的個(gè)人受照劑量、集體受照劑量、經(jīng)濟(jì)代價(jià)、環(huán)境影響和公眾心理等眾多影響因素,實(shí)現(xiàn)核設(shè)施退役過程中的輻射防護(hù)最優(yōu)化,為核設(shè)施退役提供更合理、有效的決策 方案。

    1 Leservot A, Chodorge L. Chavir: a virtual site simulation environment[C]. European Nacigation Conference, Versailles, French, 2005.

    2 Rodenas J, Zarza I, Burgos M C,. Developing a virtual reality application for training nuclear power plant operators: setting up a database containing dose rates in the refuelling plant[J]. Radiation Protection Dosimetry, 2004, 111(2): 173?180. DOI: 10.1093/rpd/nch043.

    3 Kim Y H, Park W M. Use of simulation technology for prediction of radiation dose in nuclear power plant[M]. Springer Berlin Heidelberg: Computational and Information Science, 2004: 413?418.

    4 Park H S, Kim G H, Lee K W,. Development of animation and simulation module for evaluation of Worker’s dose[C]. Waste Management Symposium, Tucson, United States, 2007.

    5 戴波, 張永領(lǐng), 周斌, 等. 反應(yīng)堆退役仿真技術(shù)研究方案[J]. 核動(dòng)力工程, 2013, 34(3): 168?171. DAI Bo, ZHANG Yongling, ZHOU Bin,. Study scheme of reactor decommissioning simulation technology[J]. Nuclear Power Engineering, 2013, 34(3): 168?171.

    6 楊子輝. 大規(guī)模輻射場景三維實(shí)時(shí)仿真關(guān)鍵技術(shù)研 究[D]. 合肥: 中國科學(xué)技術(shù)大學(xué), 2016. YANG Zihui. Study on key technologies of three-dimensional real-time simulation for large-scale radiation scenario[D]. Hefei: University of Science and Technology of China, 2016.

    7 晁楠, 劉永闊, 李夢(mèng)堃, 等. 核設(shè)施退役場景的仿真實(shí)現(xiàn)技術(shù)研究[J]. 應(yīng)用科技, 2015, 42(6): 62?66. CHAO Nan, LIU Yongkuo, LI Mengkun,. A scene simulation system for decommissioning of nuclear facilities[J]. Applied Science and Technology, 2015, 42(6): 62?66.

    8 謝小龍, 胡國輝. 虛擬現(xiàn)實(shí)技術(shù)在核設(shè)施退役中的應(yīng)用[J]. 科技創(chuàng)新導(dǎo)報(bào), 2013, 260(8): 93?94. XIE Xiaolong, HU Guohui. Application of virtual reality technology in the decommissioning nuclear facilities[J]. Science and Technology Innovation Herald, 2013, 260(8): 93?94.

    9 吳宜燦, 何桃, 胡麗琴, 等. 核與輻射安全仿真系統(tǒng)SuperMC/RVIS2.3研發(fā)與應(yīng)用[J]. 原子能科學(xué)技術(shù), 2015, 49(增1): 7?15. DOI: 10.7538/yzk.2015.49.S0.0007. WU Yican, HE Tao, HU Liqin,. Development of virtual reality-based simulation system for nuclear and radiation safety[J]. Atomic Energy Science and Technology, 2015, 49(Suppl 1): 7?15. DOI: 10.7538/yzk. 2015.49.S0.0007.

    10 席福坤, 方邦城. 計(jì)算g射線點(diǎn)核積分的蒙特卡羅方 法[J]. 核動(dòng)力工程, 1987, 8(2): 66?74. XI Fukun, FANG Bangcheng. A Monte Carlo method of the point kernel integration for gamma rays calculation[J]. Nuclear Power Engineering, 1987, 8(2): 66?74.

    11 Holland J H. Adaptation in natural and artificial systems[M]. Cambridge, Massachusetts: MIT Press, 2015.

    Radiation field reconstruction and dismantling route optimization during nuclear decommissioning

    SONG Yingming1LIANG Ye1YE Kaixuan1ZHANG Zhenyu1ZHANG Qiunan1WANG Ning2YAN Jiawei1ZHU Zhichao1ZOU Shuliang1

    1(School of Nuclear Science and Technology, University of South China, Hengyang 421001, China) 2(Center of Information, Nuclear Power Institute of China, Chengdu 610213, China)

    Background: Decommissioning nuclear facilities with less cost is a key part of nuclear reactor development, which mainly depends on the dismantling route based on the distribution of radiation field. Purpose: Thisstudy attempts to reconstruct the radiation field of decommissioning nuclear facilities and find out the optimal dismantling route. Methods: The method of point-kernel integration is adopted to reconstruct the three-dimensional radiation field of multi-source terms during the decommissioning of nuclear facilities and the results are compared with the Monte Carlo program. A functional software is designed and developed with functions of radiation field reconstruction and removal route optimization during the nuclear facility decommissioning process. The problem of optimizing dismantling route is abstracted to a mathematical model of the Traveling Salesman Problem (TSP), and the Dose Matrix under different routes are constructed. According to ALARA (As Low As Reasonably Achievable) principle of the radiation protection, the methodology is optimized by genetic algorithm. Results: The reconstructed radiation field is in good agreement with the result calculated by the Monte Carlo program. Upon the instance of decommissioning multi-source terms’ removal, the optimal route and three-dimensional visualization are calculated and obtained, and the optimal results are discussed. Conclusion: The method of point-kernel integrationand the TSP mathematical model can be used during nuclear decommissioning.

    Nuclear decommissioning, Radiation field reconstruction, Dismantling route optimization, Functional software

    TL72

    10.11889/j.0253-3219.2017.hjs.40.050502

    宋英明,男,1980年出生,2009年于清華大學(xué)獲博士學(xué)位,研究領(lǐng)域?yàn)楹嗽O(shè)施退役與輻射安全

    2016-09-20,

    2016-11-07

    Supported by Special Research Fund for the Doctoral Program of Higher Education (No.20134324120003), Key Project Fund of Hunan Provincial Department of Education (No.14A120)

    高等學(xué)校博士學(xué)科點(diǎn)專項(xiàng)科研基金(No.20134324120003)、湖南省教育廳重點(diǎn)項(xiàng)目基金(No.14A120)資助

    SONG Yingming, male, born in 1980, graduated from Tsinghua University with a doctoral degree in 2009, focusing on the field of nuclear facilities and radiation safety

    2016-09-20, accepted date: 2016-11-07

    猜你喜歡
    核設(shè)施器件遺傳算法
    核設(shè)施用爆破片ASME與GB/T 567的對(duì)比分析
    核設(shè)施退役產(chǎn)業(yè)全球經(jīng)驗(yàn)探析
    能源(2018年7期)2018-09-21 07:56:36
    基于自適應(yīng)遺傳算法的CSAMT一維反演
    核設(shè)施退役去污技術(shù)的現(xiàn)狀及發(fā)展
    中國核電(2017年2期)2017-08-11 08:01:17
    一種基于遺傳算法的聚類分析方法在DNA序列比較中的應(yīng)用
    基于遺傳算法和LS-SVM的財(cái)務(wù)危機(jī)預(yù)測
    旋涂-蒸鍍工藝制備紅光量子點(diǎn)器件
    基于改進(jìn)的遺傳算法的模糊聚類算法
    面向高速應(yīng)用的GaN基HEMT器件
    一種加載集總器件的可調(diào)三維周期結(jié)構(gòu)
    超碰97精品在线观看| 男女无遮挡免费网站观看| 亚洲精品一二三| 亚洲在久久综合| 1000部很黄的大片| 91aial.com中文字幕在线观看| 91aial.com中文字幕在线观看| 一级片'在线观看视频| 91精品伊人久久大香线蕉| 韩国av在线不卡| 国产免费又黄又爽又色| 夫妻性生交免费视频一级片| 日产精品乱码卡一卡2卡三| 亚洲综合色惰| 亚洲色图综合在线观看| 啦啦啦啦在线视频资源| 精品视频人人做人人爽| 寂寞人妻少妇视频99o| 看非洲黑人一级黄片| 亚洲精品国产成人久久av| 亚洲四区av| 尤物成人国产欧美一区二区三区| 亚洲欧美日韩东京热| 国产永久视频网站| 一边亲一边摸免费视频| 好男人在线观看高清免费视频| 97在线人人人人妻| 国产白丝娇喘喷水9色精品| 亚洲精品亚洲一区二区| 国产成人免费观看mmmm| 91在线精品国自产拍蜜月| 国产黄色视频一区二区在线观看| 又爽又黄a免费视频| 久久久久国产网址| 97人妻精品一区二区三区麻豆| 久久久久久国产a免费观看| 高清欧美精品videossex| 欧美高清性xxxxhd video| 亚洲欧美成人精品一区二区| 国内少妇人妻偷人精品xxx网站| 七月丁香在线播放| av国产久精品久网站免费入址| 免费看av在线观看网站| 中国美白少妇内射xxxbb| 国产乱人视频| 亚洲精品国产成人久久av| 国产精品久久久久久精品电影小说 | 免费大片18禁| 国产精品一二三区在线看| 国产亚洲午夜精品一区二区久久 | 国产免费福利视频在线观看| 菩萨蛮人人尽说江南好唐韦庄| 天美传媒精品一区二区| 成人鲁丝片一二三区免费| 青春草视频在线免费观看| 亚洲成人av在线免费| eeuss影院久久| 深爱激情五月婷婷| 欧美丝袜亚洲另类| 欧美日本视频| 亚洲真实伦在线观看| 日韩av免费高清视频| 免费观看的影片在线观看| 国产欧美日韩一区二区三区在线 | 国模一区二区三区四区视频| 国产精品伦人一区二区| 最近最新中文字幕大全电影3| 极品少妇高潮喷水抽搐| 高清在线视频一区二区三区| 欧美日韩亚洲高清精品| 久久久欧美国产精品| 免费黄网站久久成人精品| 国产91av在线免费观看| 精品人妻一区二区三区麻豆| 欧美区成人在线视频| 精品一区二区免费观看| 午夜福利视频1000在线观看| 欧美一级a爱片免费观看看| 国产熟女欧美一区二区| 国产男人的电影天堂91| 亚洲欧美一区二区三区国产| 久久久精品免费免费高清| 日韩欧美 国产精品| 一边亲一边摸免费视频| 麻豆乱淫一区二区| 亚洲在久久综合| 五月伊人婷婷丁香| 一区二区三区四区激情视频| 国产精品久久久久久久久免| 成人亚洲精品一区在线观看 | 女的被弄到高潮叫床怎么办| 午夜激情久久久久久久| 成人国产麻豆网| 狂野欧美激情性xxxx在线观看| 丰满少妇做爰视频| 特大巨黑吊av在线直播| 婷婷色av中文字幕| 亚洲人成网站高清观看| 波多野结衣巨乳人妻| 亚洲人成网站在线观看播放| 国产色爽女视频免费观看| a级一级毛片免费在线观看| 亚洲精品一二三| 国产伦理片在线播放av一区| 久久精品国产鲁丝片午夜精品| 老师上课跳d突然被开到最大视频| 一级爰片在线观看| 赤兔流量卡办理| 联通29元200g的流量卡| 秋霞伦理黄片| 秋霞伦理黄片| 欧美+日韩+精品| 搞女人的毛片| 亚洲精华国产精华液的使用体验| 美女被艹到高潮喷水动态| 久久精品熟女亚洲av麻豆精品| 精品99又大又爽又粗少妇毛片| 久久久久久久国产电影| av在线app专区| 别揉我奶头 嗯啊视频| 久久久久久久大尺度免费视频| av线在线观看网站| 别揉我奶头 嗯啊视频| 成年av动漫网址| 精品99又大又爽又粗少妇毛片| 国产在线一区二区三区精| 国产欧美日韩一区二区三区在线 | 99热6这里只有精品| 超碰97精品在线观看| 在线观看人妻少妇| 一级毛片aaaaaa免费看小| 国产在线男女| 97超碰精品成人国产| 伊人久久精品亚洲午夜| 国产精品.久久久| 你懂的网址亚洲精品在线观看| 色网站视频免费| 中文天堂在线官网| 国产精品伦人一区二区| 一个人观看的视频www高清免费观看| 久久久成人免费电影| 国产免费一级a男人的天堂| 久热这里只有精品99| 嘟嘟电影网在线观看| 91精品伊人久久大香线蕉| av又黄又爽大尺度在线免费看| eeuss影院久久| 五月伊人婷婷丁香| 国产精品av视频在线免费观看| 国产精品偷伦视频观看了| 欧美另类一区| 亚洲av二区三区四区| 亚洲精品日韩在线中文字幕| 美女国产视频在线观看| 精品久久国产蜜桃| 王馨瑶露胸无遮挡在线观看| 秋霞在线观看毛片| 国产精品伦人一区二区| 别揉我奶头 嗯啊视频| 日韩在线高清观看一区二区三区| 亚洲国产精品999| 日本黄大片高清| 亚洲精品成人久久久久久| 大码成人一级视频| 一区二区三区免费毛片| 成人鲁丝片一二三区免费| 久久久a久久爽久久v久久| 亚洲国产精品999| .国产精品久久| 亚洲真实伦在线观看| 26uuu在线亚洲综合色| 一区二区av电影网| 少妇人妻一区二区三区视频| 亚洲精品乱码久久久久久按摩| 日韩欧美 国产精品| 午夜福利在线在线| 日韩av不卡免费在线播放| 高清午夜精品一区二区三区| 国产欧美日韩一区二区三区在线 | 日日摸夜夜添夜夜爱| 综合色av麻豆| 最近最新中文字幕免费大全7| 最近最新中文字幕免费大全7| 美女主播在线视频| 国产极品天堂在线| 国产白丝娇喘喷水9色精品| a级一级毛片免费在线观看| 国产色爽女视频免费观看| 国产在线一区二区三区精| 香蕉精品网在线| 蜜桃久久精品国产亚洲av| 国产在线一区二区三区精| 91精品国产九色| 亚洲成人一二三区av| 色婷婷久久久亚洲欧美| 性色avwww在线观看| 久久精品久久精品一区二区三区| 亚洲av电影在线观看一区二区三区 | 国产69精品久久久久777片| 亚洲人与动物交配视频| 丰满少妇做爰视频| 六月丁香七月| 男女国产视频网站| 丰满乱子伦码专区| 激情 狠狠 欧美| 色哟哟·www| 色哟哟·www| 又大又黄又爽视频免费| av免费观看日本| 亚洲精品久久久久久婷婷小说| 国产精品国产三级专区第一集| 成人毛片60女人毛片免费| 看十八女毛片水多多多| 精品酒店卫生间| 热re99久久精品国产66热6| 国产淫片久久久久久久久| 午夜福利视频1000在线观看| av国产久精品久网站免费入址| 18禁在线无遮挡免费观看视频| 自拍欧美九色日韩亚洲蝌蚪91 | 菩萨蛮人人尽说江南好唐韦庄| freevideosex欧美| 欧美精品人与动牲交sv欧美| 天堂中文最新版在线下载 | 欧美激情久久久久久爽电影| 日日摸夜夜添夜夜添av毛片| 熟女av电影| 伦理电影大哥的女人| 最近中文字幕2019免费版| 日产精品乱码卡一卡2卡三| 亚洲国产av新网站| 中国三级夫妇交换| 国产成人福利小说| 亚洲精品国产色婷婷电影| 国产黄色视频一区二区在线观看| 国产精品蜜桃在线观看| 午夜精品一区二区三区免费看| 一个人看视频在线观看www免费| 插逼视频在线观看| 肉色欧美久久久久久久蜜桃 | 欧美xxⅹ黑人| 2022亚洲国产成人精品| 自拍偷自拍亚洲精品老妇| 日韩精品有码人妻一区| 国产成人福利小说| av国产精品久久久久影院| 嫩草影院新地址| 毛片一级片免费看久久久久| 亚洲国产最新在线播放| 国产美女午夜福利| 香蕉精品网在线| 搞女人的毛片| 亚洲欧美精品自产自拍| 亚洲精品视频女| 国产精品秋霞免费鲁丝片| 国产国拍精品亚洲av在线观看| 亚洲av福利一区| 亚洲欧洲日产国产| 国产精品熟女久久久久浪| 狂野欧美激情性bbbbbb| 久久精品国产a三级三级三级| 亚洲av在线观看美女高潮| 成人综合一区亚洲| 麻豆成人av视频| 日韩人妻高清精品专区| 少妇人妻久久综合中文| 国产亚洲91精品色在线| 校园人妻丝袜中文字幕| 亚洲精品日韩在线中文字幕| 久久6这里有精品| 欧美激情在线99| 一本久久精品| 国产精品99久久99久久久不卡 | 日日撸夜夜添| 精品久久国产蜜桃| 卡戴珊不雅视频在线播放| 国产淫片久久久久久久久| 超碰97精品在线观看| 九九爱精品视频在线观看| 人妻 亚洲 视频| 黄色一级大片看看| 国产高潮美女av| 男女下面进入的视频免费午夜| 免费不卡的大黄色大毛片视频在线观看| 国产亚洲最大av| 丝袜美腿在线中文| 精品久久久精品久久久| 亚洲av免费高清在线观看| 99久久精品国产国产毛片| 久久久久久久亚洲中文字幕| 欧美人与善性xxx| 国产高潮美女av| 美女主播在线视频| 啦啦啦在线观看免费高清www| 91精品伊人久久大香线蕉| 精品人妻一区二区三区麻豆| 欧美高清成人免费视频www| 欧美区成人在线视频| 一区二区三区精品91| 色吧在线观看| av在线播放精品| 最近2019中文字幕mv第一页| 国产乱人偷精品视频| 久久久久久久久久成人| 国产男女超爽视频在线观看| 我的老师免费观看完整版| 欧美bdsm另类| 色哟哟·www| 日本av手机在线免费观看| 毛片一级片免费看久久久久| 一级毛片我不卡| 亚洲欧美日韩另类电影网站 | 免费黄色在线免费观看| 色婷婷久久久亚洲欧美| 高清毛片免费看| 高清视频免费观看一区二区| 边亲边吃奶的免费视频| 亚洲精品成人久久久久久| 一级毛片电影观看| 校园人妻丝袜中文字幕| 乱系列少妇在线播放| 麻豆精品久久久久久蜜桃| 尤物成人国产欧美一区二区三区| 在线播放无遮挡| 色网站视频免费| 五月伊人婷婷丁香| 男女边吃奶边做爰视频| 亚洲av免费高清在线观看| 亚洲国产精品成人久久小说| 一级黄片播放器| 观看免费一级毛片| 搞女人的毛片| 永久网站在线| 亚洲自拍偷在线| 777米奇影视久久| 久久精品久久精品一区二区三区| 国产精品一区www在线观看| 一区二区三区免费毛片| 乱系列少妇在线播放| 超碰97精品在线观看| 国产精品久久久久久久久免| 大又大粗又爽又黄少妇毛片口| 亚洲性久久影院| 制服丝袜香蕉在线| 3wmmmm亚洲av在线观看| 免费不卡的大黄色大毛片视频在线观看| 国产亚洲最大av| 亚洲伊人久久精品综合| 亚洲欧美成人综合另类久久久| 大码成人一级视频| 精品久久国产蜜桃| 麻豆国产97在线/欧美| 欧美激情在线99| 欧美成人a在线观看| 欧美xxxx黑人xx丫x性爽| 亚洲最大成人手机在线| 大香蕉97超碰在线| 国产 精品1| 日本黄色片子视频| 亚洲av在线观看美女高潮| 午夜激情福利司机影院| 免费看av在线观看网站| 新久久久久国产一级毛片| 国产日韩欧美亚洲二区| 少妇的逼水好多| 日韩,欧美,国产一区二区三区| 成人免费观看视频高清| 夜夜看夜夜爽夜夜摸| 国产精品一区二区性色av| 亚洲精品久久久久久婷婷小说| 一区二区av电影网| 国产成人免费观看mmmm| 男女边吃奶边做爰视频| 亚洲精品国产色婷婷电影| 亚洲精品影视一区二区三区av| 国产黄色视频一区二区在线观看| 人妻系列 视频| 香蕉精品网在线| 亚洲av中文字字幕乱码综合| 夜夜看夜夜爽夜夜摸| 国产 一区精品| 亚洲一级一片aⅴ在线观看| 一级毛片久久久久久久久女| av女优亚洲男人天堂| 免费黄网站久久成人精品| 亚洲,欧美,日韩| 两个人的视频大全免费| 免费av毛片视频| 亚洲欧美日韩卡通动漫| 男插女下体视频免费在线播放| 国产人妻一区二区三区在| 亚洲图色成人| 免费看光身美女| 日本熟妇午夜| 老司机影院毛片| 国产午夜精品久久久久久一区二区三区| 在线观看av片永久免费下载| 纵有疾风起免费观看全集完整版| 日韩av不卡免费在线播放| 高清日韩中文字幕在线| 国产综合精华液| 最近的中文字幕免费完整| 尤物成人国产欧美一区二区三区| 亚洲性久久影院| 边亲边吃奶的免费视频| 久久久精品欧美日韩精品| 热99国产精品久久久久久7| 高清在线视频一区二区三区| 人妻制服诱惑在线中文字幕| 久久久亚洲精品成人影院| 中国三级夫妇交换| 一本久久精品| 真实男女啪啪啪动态图| av天堂中文字幕网| 美女国产视频在线观看| 国产成人免费无遮挡视频| 少妇 在线观看| 久久久久久久久久久免费av| 日日啪夜夜撸| 我要看日韩黄色一级片| 性色av一级| 国产伦理片在线播放av一区| 2022亚洲国产成人精品| 看非洲黑人一级黄片| av.在线天堂| 久久99精品国语久久久| 国产精品国产三级国产av玫瑰| 一区二区三区乱码不卡18| 国产久久久一区二区三区| av天堂中文字幕网| 日韩制服骚丝袜av| 国产黄色视频一区二区在线观看| 男女无遮挡免费网站观看| 日韩av不卡免费在线播放| 亚洲精品国产色婷婷电影| 99热6这里只有精品| 亚洲三级黄色毛片| 大话2 男鬼变身卡| 久久人人爽人人爽人人片va| 麻豆成人午夜福利视频| 亚洲国产成人一精品久久久| 亚洲国产精品999| 自拍偷自拍亚洲精品老妇| 午夜福利视频1000在线观看| 亚洲不卡免费看| 国产黄片美女视频| 国产av国产精品国产| 春色校园在线视频观看| 99热这里只有精品一区| 18+在线观看网站| av免费在线看不卡| 最近2019中文字幕mv第一页| 老司机影院成人| videos熟女内射| 国产有黄有色有爽视频| 国产成人免费观看mmmm| 99久国产av精品国产电影| 天堂网av新在线| 国产乱人视频| 欧美性猛交╳xxx乱大交人| 国产精品久久久久久精品古装| 99re6热这里在线精品视频| av国产免费在线观看| 国产成人精品婷婷| 国产爽快片一区二区三区| 一级毛片黄色毛片免费观看视频| 少妇人妻久久综合中文| 性色avwww在线观看| 99热这里只有是精品50| 边亲边吃奶的免费视频| 最新中文字幕久久久久| 丝袜美腿在线中文| 尾随美女入室| 少妇丰满av| 六月丁香七月| 中文字幕久久专区| 22中文网久久字幕| 九色成人免费人妻av| 亚洲国产精品成人久久小说| 人人妻人人看人人澡| 街头女战士在线观看网站| 欧美日韩视频高清一区二区三区二| 亚洲欧美日韩另类电影网站 | av国产免费在线观看| 嫩草影院新地址| 亚洲av在线观看美女高潮| 三级男女做爰猛烈吃奶摸视频| 99久久中文字幕三级久久日本| 下体分泌物呈黄色| 国产伦精品一区二区三区四那| 国产高潮美女av| 高清欧美精品videossex| 欧美97在线视频| 国产精品国产三级专区第一集| 少妇猛男粗大的猛烈进出视频 | 亚洲av二区三区四区| 免费av毛片视频| 婷婷色综合www| 亚洲丝袜综合中文字幕| 免费不卡的大黄色大毛片视频在线观看| 久久久久网色| 免费av观看视频| 亚洲精品日韩av片在线观看| 卡戴珊不雅视频在线播放| 亚洲av成人精品一区久久| 精品久久久久久久人妻蜜臀av| 只有这里有精品99| 51国产日韩欧美| 久久女婷五月综合色啪小说 | 免费av毛片视频| 午夜福利视频精品| 日韩在线高清观看一区二区三区| 欧美亚洲 丝袜 人妻 在线| 欧美激情在线99| 最近2019中文字幕mv第一页| 黄色一级大片看看| 一级a做视频免费观看| 日本黄色片子视频| 人人妻人人澡人人爽人人夜夜| 我的女老师完整版在线观看| 欧美丝袜亚洲另类| 亚洲在久久综合| 最新中文字幕久久久久| 欧美成人一区二区免费高清观看| 精品久久久噜噜| 蜜臀久久99精品久久宅男| 亚洲av.av天堂| 黄色视频在线播放观看不卡| 免费看不卡的av| 99热国产这里只有精品6| 大香蕉97超碰在线| kizo精华| 成人黄色视频免费在线看| 国产免费一区二区三区四区乱码| 国产亚洲av片在线观看秒播厂| 成人毛片a级毛片在线播放| 精品久久国产蜜桃| 免费看av在线观看网站| 天天一区二区日本电影三级| 亚洲欧美日韩无卡精品| 在线亚洲精品国产二区图片欧美 | 久久久国产一区二区| 成人特级av手机在线观看| 久久久久久久国产电影| 成人高潮视频无遮挡免费网站| 18禁裸乳无遮挡动漫免费视频 | 夜夜爽夜夜爽视频| av在线蜜桃| 女人十人毛片免费观看3o分钟| 男人添女人高潮全过程视频| 久久久午夜欧美精品| 一个人观看的视频www高清免费观看| 人人妻人人爽人人添夜夜欢视频 | 国内少妇人妻偷人精品xxx网站| 青春草亚洲视频在线观看| 久久久久久久久久久丰满| 亚洲精品日韩在线中文字幕| 亚洲婷婷狠狠爱综合网| 亚洲av中文字字幕乱码综合| 亚洲av中文av极速乱| 亚洲av免费高清在线观看| 特级一级黄色大片| 视频区图区小说| 97精品久久久久久久久久精品| av天堂中文字幕网| 综合色av麻豆| 午夜精品国产一区二区电影 | av天堂中文字幕网| 精品久久久精品久久久| av.在线天堂| 麻豆成人午夜福利视频| 日韩av不卡免费在线播放| 青春草国产在线视频| 国产午夜精品久久久久久一区二区三区| 一级av片app| 少妇人妻 视频| 免费黄频网站在线观看国产| 国产男人的电影天堂91| 国产日韩欧美在线精品| 日韩精品有码人妻一区| 男人狂女人下面高潮的视频| 欧美bdsm另类| 婷婷色av中文字幕| 日韩av在线免费看完整版不卡| 熟女人妻精品中文字幕| 男女啪啪激烈高潮av片| 日本黄色片子视频| 日韩伦理黄色片| 久久久久久久久久久丰满| 日韩大片免费观看网站| 午夜精品一区二区三区免费看| 激情五月婷婷亚洲| 久久久久久久久久久丰满| 自拍偷自拍亚洲精品老妇| 亚洲四区av| 国产男女超爽视频在线观看| 日韩av免费高清视频| 国产欧美日韩一区二区三区在线 | 亚洲,欧美,日韩| 亚洲真实伦在线观看| 中文乱码字字幕精品一区二区三区| 国产一区二区三区综合在线观看 | 色婷婷久久久亚洲欧美| 你懂的网址亚洲精品在线观看| 色哟哟·www| 80岁老熟妇乱子伦牲交| www.色视频.com| 色婷婷久久久亚洲欧美| 国产精品99久久久久久久久| 一级毛片 在线播放| 午夜老司机福利剧场| 直男gayav资源| 成人毛片a级毛片在线播放| 日日摸夜夜添夜夜爱| 午夜精品国产一区二区电影 |