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

        ?

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

        2017-10-13 09:21:58宋英明梁燁葉凱萱張震宇張秋楠王寧顏佳偉朱志超鄒樹梁
        核技術(shù) 2017年5期
        關(guān)鍵詞:劑量優(yōu)化

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

        ?

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

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

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

        輻射場(chǎng)分布和路徑選擇是核設(shè)施退役過程中影響工作人員吸收劑量的主要因素。針對(duì)核設(shè)施退役過程中多源項(xiàng)的三維輻射場(chǎ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è)施退役過程的輻射場(chǎng)重構(gòu)與拆除路徑優(yōu)化功能軟件。

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

        核設(shè)施退役是一項(xiàng)周期長(zhǎ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)境的 危害。

        國(guó)際上對(duì)核設(shè)施退役仿真系統(tǒng)進(jìn)行了廣泛和深入的研究。法國(guó)原子能委員會(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ì)劃等功能。韓國(guó)原子能研究所和忠南大學(xué)[3]聯(lián)合對(duì)韓國(guó)KRR-1&2反應(yīng)堆退役開展工人輻射劑量計(jì)算機(jī)模擬,能實(shí)現(xiàn)基于MCNP-4C[4]輻射輸運(yùn)程序下的人的輻射劑量計(jì)算。隨著我國(guó)核能事業(yè)的高速發(fā)展,對(duì)核設(shè)施退役技術(shù)研究與工程準(zhǔn)備工作的要求更加體系化、專業(yè)化、標(biāo)準(zhǔn)化。目前,國(guó)內(nèi)已有部分單位在輻射場(chǎng)可視化軟件開發(fā)和核設(shè)施虛擬退役等方面開展了相關(guān)研究,取得了較好的成果[5?9],但對(duì)于核設(shè)施退役中拆除路徑最優(yōu)化方法的研究相對(duì)較少且缺乏針對(duì)性。本文驗(yàn)證并采用點(diǎn)核積分法重構(gòu)核設(shè)施退役過程中多源項(xiàng)的三維輻射場(chǎng),然后利用類旅行商問題(Traveling Salesman Problem, TSP)數(shù)學(xué)模型求解核退役拆除路徑問題,構(gòu)造不同拆除路徑下所受外照射劑量矩陣,通過遺傳算法尋找最優(yōu)路徑。設(shè)計(jì)開發(fā)了軟件實(shí)現(xiàn)相應(yīng)功能并進(jìn)行了實(shí)例分析。

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

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

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

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

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

        利用MCNP5驗(yàn)證簡(jiǎn)單源項(xiàng)情況下點(diǎn)核積分的結(jié)果,為簡(jiǎn)化計(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 核退役拆除路徑問題的簡(jiǎn)述

        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ìng)爭(zhēng)的方法進(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ì)算輻射場(chǎng),利用功能軟件進(jìn)行二次尋優(yōu)。

        4 功能軟件開發(fā)

        利用MATLAB實(shí)現(xiàn)了基于遺傳算法求解旅行商問題下的最優(yōu)退役路線,并通過GUI設(shè)計(jì)了面向核設(shè)施退役過程的輻射場(chǎng)重構(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)的三維輻射場(chǎng),發(fā)展了基于點(diǎn)核積分的輻射場(chǎng)快速重構(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è)施退役過程的輻射場(chǎng)重構(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ī)模輻射場(chǎng)景三維實(shí)時(shí)仿真關(guān)鍵技術(shù)研 究[D]. 合肥: 中國(guó)科學(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è)施退役場(chǎng)景的仿真實(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 謝小龍, 胡國(guó)輝. 虛擬現(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

        猜你喜歡
        劑量優(yōu)化
        課堂內(nèi)外·初中版(科學(xué)少年)(2023年10期)2023-12-10 00:43:06
        ·更正·
        中藥的劑量越大、療效就一定越好嗎?
        超限高層建筑結(jié)構(gòu)設(shè)計(jì)與優(yōu)化思考
        民用建筑防煙排煙設(shè)計(jì)優(yōu)化探討
        關(guān)于優(yōu)化消防安全告知承諾的一些思考
        一道優(yōu)化題的幾何解法
        不同濃度營(yíng)養(yǎng)液對(duì)生菜管道水培的影響
        由“形”啟“數(shù)”優(yōu)化運(yùn)算——以2021年解析幾何高考題為例
        90Sr-90Y敷貼治療的EBT3膠片劑量驗(yàn)證方法
        岛国av一区二区三区| 精品亚洲成a人在线观看青青| 国产98色在线 | 国产| 国产精品熟女一区二区| 99精品视频69V精品视频| 乱子伦av无码中文字幕| 青青自拍视频成人免费观看| 国产三级视频不卡在线观看| 日韩人妻少妇一区二区三区| 不卡高清av手机在线观看| 国产亚洲美女精品久久| 亚洲国产成人久久精品美女av| 蜜桃91精品一区二区三区| 一本久久a久久精品vr综合| 日韩插啊免费视频在线观看| 91亚洲人成手机在线观看| 午夜视频在线观看国产| 国产欧美日韩中文久久| 日韩欧美亚洲综合久久影院ds| 国产70老熟女重口小伙子| 人妻AV无码一区二区三区奥田咲| 国产精品亚洲综合色区丝瓜 | 激情亚洲一区国产精品| 99久久国语露脸国产精品| 国产av午夜精品一区二区入口| 日本少妇又色又爽又高潮| 九九99久久精品国产| 人妻av一区二区三区av免费| 五月婷婷开心五月播五月| 色综合久久久久综合99| 欧美丰满熟妇aaaaa片| 亚洲无码视频一区:| 丰满又紧又爽又丰满视频| 99久久超碰中文字幕伊人| 亚洲美女影院| 久久国产精品免费久久久| 中国孕妇变态孕交xxxx| 无码人妻品一区二区三区精99 | 日韩av在线免费观看不卡| 大又大又粗又硬又爽少妇毛片| 国产成人无码一区二区在线播放|