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

    基于GeTHU測量確定natPb(p,x)207Bi和natPb(p,x)194Hg的反應(yīng)截面

    2023-09-21 09:20:52劉奇澤歐宇福景明坤代文翰馬豪TitarenkoYuPavlovTitarenkoYuZhivunKovalishinKulevoyKhalikov
    核技術(shù) 2023年9期
    關(guān)鍵詞:核素活度質(zhì)子

    劉奇澤 歐宇福 景明坤 代文翰 曾 志 薛 濤 田 陽 曾 鳴 馬豪 Titarenko Yu E Pavlov K V Titarenko A Yu Zhivun V M Kovalishin A A Kulevoy T V Khalikov R S

    1(清華大學(xué) 工程物理系 北京 100084)

    2(NRC "Kurchatov Institute" 俄羅斯 莫斯科 123182)

    加速器驅(qū)動次臨界系統(tǒng)(Accelerator Driven Sub-critical System,ADS)具備提升反應(yīng)堆安全性和實現(xiàn)放射性廢物嬗變的固有能力[1-2],因此成為下一代核能系統(tǒng)最具潛力的選型之一。ADS 有三個主要部件用于在次臨界狀態(tài)下驅(qū)動系統(tǒng)運行,包括高能質(zhì)子加速器、散裂靶和次臨界堆芯,其中,散裂靶通常由鉛(Pb)和鉍(Bi)構(gòu)成。當(dāng)高能質(zhì)子轟擊散裂靶時,在散裂靶中會不斷產(chǎn)生和積累長壽命放射性核素,一些由Pb(p,x)/Bi(p,x)反應(yīng)所形成的放射性核素產(chǎn)物的壽命可達(dá)數(shù)百年。在儲存和管理放射性廢物時,必須考慮這些長壽命放射性核素的產(chǎn)生。在2000~2002 年期間,俄羅斯理論和實驗物理研究所(NRC "Kurchatov Institute")對此開展了一系列相關(guān)的實驗研究,獲得了一系列較短壽命核素(如24Na、59Fe、101mRh 等)的截面結(jié)果[3],并總結(jié)為IAEA(International Atomic Energy Agency)的INDC(International Nuclear Data Committee)報告[4]。但由于受照后的鉛靶內(nèi)存在較多來自短壽命放射性核素的本底,一部分截面較小、活度較低的長壽命產(chǎn)物(如207Bi(T1/2=31.55 a)和194Hg(T1/2=447 a)等)難以準(zhǔn)確計算反應(yīng)截面,而這些長壽命核素的產(chǎn)生截面[5-9]在EXFOR 等核數(shù)據(jù)庫中又較為匱乏,亟待實驗測定。這些受照鉛靶在經(jīng)過約20 a 的冷卻后,來自短壽命核素的本底已大幅降低,為測量這部分長壽命產(chǎn)物的截面提供了可能,但是對這類低活度長壽命產(chǎn)物的測量易受地面譜儀較高的本底水平和較復(fù)雜的本底環(huán)境干擾,為此中俄雙方作者展開合作,分別在地上實驗室(Institute for Theoretical and Experimental Physics,ITEP)和中國錦屏地下實驗室(China Jinping Underground Laboratory,CJPL)對相同的質(zhì)子轟擊鉛靶進(jìn)行測量:由俄方Titarenko 等在ITEP 利用高純鍺γ 譜儀對樣品進(jìn)行測量,分別給出了產(chǎn)物為207Bi[10]和194Hg[11]對應(yīng)的反應(yīng)截面;本工作則依托CJPL 中的超低本底γ 譜儀對Titarenko 等工作中測量的4 個標(biāo)記為natPb400、natPb100、natPb70和natPb40 的鉛靶樣品進(jìn)行平行測量,在更低的本底水平[11]下分析了207Bi和194Hg的產(chǎn)生截面,并同ITEP的發(fā)表結(jié)果進(jìn)行了比較。

    1 測量與分析方法

    1.1 反應(yīng)截面計算

    為了計算natPb(p,x)207Bi 和natPb(p,x)194Hg 的反應(yīng)截面,使用以下公式:

    式中:σ表示特定產(chǎn)物核素的反應(yīng)截面;R為產(chǎn)物核素的產(chǎn)生率;Φ為質(zhì)子的通量密度。

    質(zhì)子的Φ已在實驗中測定,產(chǎn)物核素的R可通過式(2)計算。

    對于活度可由其自身γ 線確定的反應(yīng)產(chǎn)物(如207Bi),根據(jù)衰變規(guī)律使用以下公式計算反應(yīng)速率R:

    式中:A為特定產(chǎn)物核素的活度;t為質(zhì)子輻照實驗與樣品測量實驗之間的時間間隔(冷卻時間);λ為產(chǎn)物核素的衰變常數(shù);tirr為樣品的質(zhì)子輻照時間;nPb為樣品中的鉛原子數(shù)。

    對于活度不直接由其γ 線確定的產(chǎn)物(如194Hg),根據(jù)與測量相關(guān)的衰變過程計算反應(yīng)速率R。對于194Hg,其產(chǎn)生速率可以根據(jù)遞次衰變規(guī)則通過以下公式計算:

    式中:AAu為194Hg 的子核194Au 的活度;λHg為194Hg 的衰變常數(shù);λAu為194Au的衰變常數(shù)。

    考慮到λHg?λAu,式(3)以式(2)的形式近似表達(dá),此時,A為194Au 的活度值,λ為194Hg 的衰變常數(shù)值。因此,在近似之后,對于natPb(p,x)207Bi 和natPb(p,x)194Hg 兩類反應(yīng),均可用式(2)描述核素的產(chǎn)生速率。

    通過γ 譜分析可以獲得計算反應(yīng)速率R所需的核素活度A,其公式如下:

    式中:A為感興趣核素的活度;Nnet為實測譜中某一特定全能峰下的凈計數(shù);ε為對應(yīng)全能峰的模擬探測效率;t為樣品測量過程的活時間。

    1.2 樣品

    測量所使用的樣品由NRC "Kurchatov Institute"所提供,4個質(zhì)子輻照天然鉛樣品分別命名為natPb40(使用40 MeV質(zhì)子輻照)、natPb70(使用70 MeV質(zhì)子輻照)、natPb100(使用100 MeV 質(zhì)子輻照)和natPb400(使用400 MeV質(zhì)子輻照)。所有的待測樣品均單獨密封在塑料樣品袋中,樣品袋的整體厚度較?。s0.1 mm),經(jīng)模擬對兩種感興趣核素能區(qū)的影響小于3%,可以忽略樣品袋造成的衰減。

    用于輻照的鉛樣品使用高富集度的天然鉛靶丸壓制而成,直徑均為10.5 mm,其純度使用質(zhì)譜儀檢定,其雜質(zhì)水平列于表1,未列出的其他雜質(zhì)在樣品中的占比更低。表2列出了樣品在質(zhì)子輻照實驗中的具體輻照參數(shù),這些參數(shù)均由文獻(xiàn)[4]給出。

    表1 樣品雜質(zhì)水平Table 1 Chemical impurity of irradiated samples

    表2 樣品輻照參數(shù)Table 2 Parameters of irradiated samples

    1.3 GeTHU-I超低本底γ譜儀

    CJPL 是世界上最深的地下實驗室[12-14],位于中國四川省,巖石埋深最高可達(dá)2 400 m(等效水厚6 720 m),得益于良好的低本底環(huán)境,CJPL 成為諸多重要物理前沿實驗的研究平臺[15-17]。在2020 年10 月至2021 年3 月期間,所有樣品均在CJPL 使用GeTHU-I探測器進(jìn)行測量。GeTHU-I探測器是由法國CANBERRA 公司制造的n 型同軸HPGe 探測器,相對探測效率為40%。鍺晶體的高度和直徑分別為59.8 mm和59.9 mm。為了降低探測器樣品腔內(nèi)的γ和中子本底,從內(nèi)到外分別使用5 cm 高純無氧銅、15 cm 鉛和10 cm 硼化聚乙烯板對探測器四周進(jìn)行屏蔽。關(guān)于GeTHU-I 系統(tǒng)的更多細(xì)節(jié)見文獻(xiàn)[18]。下文中提及的GeTHU 均指GeTHU-I 超低本底γ譜儀。

    各樣品分別放置在距離GeTHU 探頭低溫腔0 mm 處和40 mm 處進(jìn)行測量,此外GeTHU 樣品腔內(nèi)的本底也在樣品測量前進(jìn)行了預(yù)先測量。natPb400 樣品在探頭低溫腔表面0 mm 高度處測量7 d 和在探頭低溫腔表面40 mm 高度處測量14 d 的全能譜以及使用GeTHU 測量15 d 的本底能譜如圖1 所示。測量7 d 時間的natPb100、natPb70 和natPb40 樣品的全能譜以及本底譜如圖2所示。測量到的感興趣核素及其半衰期與全能峰峰位列于表3[19]。感興趣核素及其他核素在能譜中的特征峰已在圖3 中標(biāo)記。

    圖2 GeTHU測量natPb100、natPb70、natPb40樣品的全能譜與GeTHU本底譜的對比(彩圖見網(wǎng)絡(luò)版)Fig.2 Full gamma-spectra of natPb100, natPb70 and natPb40 samples measured by GeTHU νs. background spectrum of GeTHU (color online)

    表3 樣品中感興趣核素及全能峰能量Table 3 Peak energy of nuclides of interest in samples

    在natPb400 樣品的γ 測量譜中,207Bi 可以由其自身衰變釋放的γ 射線確定,194Hg 也可通過其衰變子代194Au的γ射線識別。與之相比,在受到100 MeV、70 MeV 和40 MeV 質(zhì) 子 照 射 的natPb100、natPb70和natPb40樣品的γ測量譜中,由于沒有對應(yīng)γ射線存在,無法識別出其中是否存在194Hg。但對于這些樣品而言,207Bi衰變釋放的γ射線依然顯著存在于樣品譜中。

    1.4 探測效率模擬

    GeTHU 的 幾 何 模 型 通 過SAGE[20]和Geant4 11.0.1[21-23]建立,SAGE 是用于HPGe 探測器實驗的蒙特卡羅模擬框架,在CJPL內(nèi)使用的多種探測器模型已在SAGE 框架內(nèi)預(yù)置,Geant4 則是歐洲核子研究 中 心(Conseil Européen pour la Recherche Nucléaire,CERN)開發(fā)的通用蒙特卡羅模擬工具包,用于模擬粒子在物質(zhì)中的相互作用和粒子輸運過程。GeTHU 中的樣品位置和GeTHU 的HPGe 探測器結(jié)構(gòu)分別如圖4 和圖5 所示。圖4 展示了在模擬中探測器的結(jié)構(gòu)和與樣品之間的位置關(guān)系,在圖4的左圖中,樣品直接置于探測器的低溫腔表面上,在圖4 的右圖中,樣品置于空心支架上距離探測器低溫腔表面40 mm高度處。所有的待測樣品均在這兩個不同的位置進(jìn)行了模擬,模擬還考慮了探測器的周圍材料及natPb400、natPb100、natPb70、natPb40 樣品的形狀參數(shù)[4]。

    圖4 探測器結(jié)構(gòu)及與樣品的相對位置Fig.4 Detector structure and relative position to the sample

    圖5 GeTHU幾何的模擬截面圖Fig.5 Cross-section view of GeTHU in Geant4

    SAGE模擬框架為模擬中的幾何和源項配置提供了靈活的接口。在SAGEJson文件(SAGE模擬框架的輸入文件)中設(shè)置好探測器和樣品的幾何參數(shù)以及源項的放射性核素類型與數(shù)量后,使用SAGE模擬了各樣品中感興趣核素的特征γ射線的探測效率。為獲得較好的統(tǒng)計性,每次模擬的粒子數(shù)量為5×106,模擬特征峰的統(tǒng)計漲落小于2%。

    使用下式計算蒙特卡羅模擬得到的探測效率:

    式中:ε表示某一特征峰的探測效率;Anet表示模擬譜中對應(yīng)全能峰的凈計數(shù);N表示模擬的粒子數(shù)。

    使用SAGE 中的衰變模擬模式分別對207Bi和194Hg兩種核素進(jìn)行了模擬,模擬假設(shè)核素均勻分布在鉛樣品內(nèi)。衰變模式的模擬從衰變母核開始,因此在蒙特卡羅模擬中已經(jīng)考慮了γ射線級聯(lián)符合的影響。與特定γ 射線的直接模擬相比,在該模式下獲得的結(jié)果已經(jīng)包含了分支比和級聯(lián)符合的校正。此外,樣品位置不確定度所產(chǎn)生的影響也在模擬中進(jìn)行了評估,并將這一影響包含在了效率的不確定度中。

    表4 和表5 分別列出了使用衰變模擬模式在0 mm和40 mm高度下的模擬探測效率。

    表4 0 mm高度處的模擬探測效率Table 4 Detection efficiency at 0 mm height

    表5 40 mm高度處的模擬探測效率Table 5 Detection efficiency at 40 mm height

    1.5 不確定度評估

    實驗測量的根本目的是確定待測物理量的值(此處為感興趣產(chǎn)物核素的反應(yīng)截面),但測量結(jié)果僅是待測物理量的近似值或估計值。因此,只有在聲明了待測量的不確定度后測量結(jié)果的報告才最終完成。

    根據(jù)§1.1 中所提到的計算過程,結(jié)果的不確定度可以通過各已知量的不確定度傳遞來計算:

    式中:ΔA為核素活度的不確定度;ΔN為全能峰凈計數(shù)的不確定度;Δε為模擬探測效率的不確定度;ΔR為反應(yīng)產(chǎn)率的不確定度:ΔT為核素半衰期的不確定度:Δσ為反應(yīng)截面的不確定度;ΔΦ為質(zhì)子通量密度的不確定度。

    2 結(jié)果與討論

    2.1 natPb(p,x)207Bi的反應(yīng)截面

    使用能量為569.7 keV(97.75%)和1 063.7 keV(74.5%)的γ 射 線,對natPb400、natPb100、natPb70和natPb40 樣品測量譜中的207Bi 進(jìn)行了測定。表6 中列出了根據(jù)式(2)計算的207Bi核素活度。

    表6 207Bi核素的活度Table 6 Activity of 207Bi nuclide

    基于207Bi 的活度值和對應(yīng)樣品的模擬結(jié)果,根據(jù)式(1)計算了不同質(zhì)子能量下natPb(p,x)207Bi 的反應(yīng)截面。計算結(jié)果見表7。

    其 他 作 者 的 實 驗 結(jié) 果[10]、Geant4(INCL++/ABLA)的計算結(jié)果以及質(zhì)子能量低于200 MeV 的TENDL-2019 庫和TENDL-2021 庫的截面評價結(jié)果[24]與本工作獲得的結(jié)果進(jìn)行了比較,如圖6所示。

    圖6 本工作natPb(p,x)207Bi截面計算結(jié)果(紅色)同其他作者的實驗結(jié)果(藍(lán)色)、TENDL-2019的評價結(jié)果(品紅色)、TENDL-2021的評價結(jié)果(綠色)以及Geant4的模擬計算結(jié)果(灰色)的比較(彩圖見網(wǎng)絡(luò)版)Fig.6 Calculated cross-sections of natPb(p,x)207Bi reaction in this research (red) νs. the results of other authors (blue), the evaluated cross sections from TENDL-2019 (magenta) and TENDL-2021 (green), and simulated data by Geant4 (gray)(color online)

    2.2 natPb(p,x)194Hg的反應(yīng)截面

    使用能量為328.5 keV(62.8%)和293.5 keV(10.9%)的γ 射線,對natPb400 的測量譜中的194Hg(194Au)進(jìn)行了測定。如§1.3所述,在natPb100、natPb70和natPb40 樣品的測量譜中,未發(fā)現(xiàn)194Hg(194Au)核素的存在。表8 中列出了根據(jù)式(2)近似計算的194Hg核素活度。

    表8 194Hg核素的活度Table 8 194Hg nuclide activity

    基于194Hg(194Au)的活度和對應(yīng)樣品的模擬結(jié)果,根據(jù)式(1)計算了不同質(zhì)子能量下natPb(p,x)194Hg的反應(yīng)截面。計算結(jié)果見表9。

    表9 natPb(p,x)194Hg的反應(yīng)截面Table 9 Cross-sections of natPb(p,x)194Hg

    類似的,其他作者的實驗結(jié)果[11]、Geant4(INCL++/ABLA)的計算結(jié)果以及質(zhì)子能量低于200 MeV 的TENDL-2019 庫 和TENDL-2021 庫 的 截面評價結(jié)果與本工作獲得的結(jié)果進(jìn)行了比較,如圖7所示。

    圖7 本工作natPb(p,x)194Hg截面計算結(jié)果(紅色)同其他作者的實驗結(jié)果(藍(lán)色)、TENDL-2019的評價結(jié)果(品紅色)、TENDL-2021的評價結(jié)果(綠色)以及Geant4的模擬計算結(jié)果(灰色)的比較(彩圖見網(wǎng)絡(luò)版)Fig.7 Calculated cross-sections of natPb(p,x)194Hg reaction in this research (red) νs. the results of other authors (blue), the evaluated cross sections from TENDL-2019 (magenta) and TENDL-2021 (green), and simulated data by Geant4 (gray)(color online)

    2.3 不同來源的不確定度貢獻(xiàn)

    在§1.5 中描述了結(jié)果不確定度的評估方法,在本工作中,截面的不確定度主要包括活度的不確定度ΔσA、半衰期的不確定度ΔσT和質(zhì)子通量密度的不確定度ΔσΦ。為了更容易與其他實驗結(jié)果進(jìn)行比較,將結(jié)果的不確定度劃分為不同的不確定度來源,如表10和表11所示??偛淮_定度Δσsum由各不確定度分量的平方和表示。

    表10 natPb(p,x)207Bi截面不確定度的組成Table 10 Uncertainty composition of results for natPb(p,x)207Bi

    表11 natPb(p,x)194Hg截面不確定度的組成Table 11 Uncertainty composition of results for natPb(p,x)194Hg

    對于natPb(p,x)207Bi 的截面,質(zhì)子通量密度的不確定度是總不確定度的主要來源。此外,活度的不確定度也占據(jù)了總不確定度的一部分。由于207Bi的半衰期數(shù)據(jù)較為準(zhǔn)確,半衰期帶來的不確定度幾乎可以忽略不計。natPb(p,x)194Hg 的截面不確定度與207Bi的情況有所不同。194Hg的半衰期具有較大的相對不確定度,這使其成為總不確定度的主要來源。質(zhì)子通量密度和活度的不確定度也占據(jù)總不確定度的一部分。

    活度的不確定度在本工作中是來自于實驗和分析的關(guān)鍵不確定度來源,有必要根據(jù)式(4)進(jìn)一步分析來自于全能峰凈計數(shù)和探測效率的不確定度來源(分別表示為δN和δε),其中探測效率的不確定度分為統(tǒng)計不確定度(statistical uncertainty,stat.)和系統(tǒng)不確定度(systematic uncertainty,syst.)兩部分,表12和表13 分別列出了在兩種核素的分析中活度不確定度的各組成部分。

    表12 natPb(p,x)194Hg活度不確定度的組成Table 12 Uncertainty composition of activity for natPb(p,x)194Hg

    表13 natPb(p,x)207Bi活度不確定度的組成Table 13 Uncertainty composition of activity for natPb(p,x)207Bi

    3 結(jié)語

    使 用GeTHU 超低本 底γ 譜儀在CJPL 測 量了natPb400、natPb100、natPb70和natPb40 4個質(zhì)子轟擊天然鉛樣品的能譜,所有樣品的能譜測量分別在兩個位置進(jìn)行:在探測器低溫腔0 mm高度處和探測器低溫腔上方40 mm 高度處。使用蒙特卡羅模擬框架SAGE中的衰變模擬模式計算對應(yīng)放射性核素的探測效率,通過γ 譜分析計算natPb400、natPb100、natPb70和natPb40樣品中感興趣核素的活度。在4種不同能量質(zhì)子(40 MeV、70 MeV、100 MeV 和 400 MeV)轟擊下,天然鉛靶中natPb(p,x)207Bi 反應(yīng)的總產(chǎn)生截面分別為(40.70±3.59) mb、(19.31±1.43) mb、(13.15±0.96) mb 和(2.90±0.22) mb。在400 MeV 能量的質(zhì)子轟擊下,天然鉛靶中natPb(p,x)194Hg反應(yīng)的總產(chǎn)生截面為(57.07±7.83) mb,在100 MeV 能量以下的質(zhì)子轟擊天然鉛靶中沒有顯著的194Hg產(chǎn)物信號發(fā)現(xiàn)。

    本文所進(jìn)行的實驗和分析表明,使用CJPL的超低本底γ 譜儀GeTHU,可以探測和識別樣品中活度較低的長壽命核素,并計算相應(yīng)核素的反應(yīng)截面。通過本實驗的實驗結(jié)果,可以推測natPb 與高能質(zhì)子反應(yīng)最終生成194Hg 的各類反應(yīng)在100 MeV 的質(zhì)子入射能以下的反應(yīng)截面小于當(dāng)前探測器的靈敏度,部分反應(yīng)的閾值能量可能超過100 MeV。

    致謝感謝曾煒赫對蒙特卡羅模擬工作的幫助。

    作者貢獻(xiàn)聲明劉奇澤負(fù)責(zé)實施研究,分析/解釋數(shù)據(jù),起草文章;歐宇福負(fù)責(zé)醞釀和設(shè)計實驗,實施研究,起草文章;景明坤負(fù)責(zé)實施研究,采集數(shù)據(jù);代文翰負(fù)責(zé)實施研究,分析/解釋數(shù)據(jù);曾志負(fù)責(zé)對文章的知識性內(nèi)容作批評性審閱,獲取研究經(jīng)費,行政、技術(shù)或材料支持,指導(dǎo);薛濤負(fù)責(zé)指導(dǎo),支持性貢獻(xiàn);田陽負(fù)責(zé)指導(dǎo),支持性貢獻(xiàn);曾鳴負(fù)責(zé)指導(dǎo),支持性貢獻(xiàn);馬豪負(fù)責(zé)指導(dǎo),支持性貢獻(xiàn);Titarenko Yu E負(fù)責(zé)對文章的知識性內(nèi)容作批評性審閱,獲取研究經(jīng)費,行政、技術(shù)或材料支持,指導(dǎo);Pavlov K V、Titarenko A Yu、Zhivun V M、Kovalishin A A、Kulevoy T V和Khalikov R S負(fù)責(zé)指導(dǎo),支持性貢獻(xiàn)。

    猜你喜歡
    核素活度質(zhì)子
    核素分類開始部分的6種7核素小片分布
    質(zhì)子束放療在腫瘤中的研究新進(jìn)展
    核素分類的4量子數(shù)
    CaO-SiO2-FeO-P2O5-Al2O3脫磷渣系中組元活度的計算
    淺談質(zhì)子守恒
    核電廠惰性氣體排放活度濃度的估算
    “質(zhì)子”號一箭發(fā)雙星
    太空探索(2014年6期)2014-07-10 13:06:11
    鉆井液處理劑溶液活度測量方法對比
    植物對核素鍶的吸附與富集作用研究現(xiàn)狀
    核素分布的三循環(huán)
    人妻夜夜爽99麻豆av| 国产中年淑女户外野战色| 国产精品一及| 色综合色国产| 欧洲精品卡2卡3卡4卡5卡区| 国内毛片毛片毛片毛片毛片| 精品人妻偷拍中文字幕| 99久久无色码亚洲精品果冻| 欧美极品一区二区三区四区| 欧美性猛交黑人性爽| 成人美女网站在线观看视频| 国产高清视频在线观看网站| 色哟哟哟哟哟哟| 精品久久久噜噜| 色视频www国产| 免费人成视频x8x8入口观看| 国产伦精品一区二区三区视频9| 在线观看66精品国产| 九色国产91popny在线| 非洲黑人性xxxx精品又粗又长| 欧美性感艳星| 一区福利在线观看| ponron亚洲| 欧美精品啪啪一区二区三区| 欧美另类亚洲清纯唯美| 午夜福利在线观看吧| 丰满人妻一区二区三区视频av| .国产精品久久| 国内精品久久久久精免费| 男人和女人高潮做爰伦理| 亚洲内射少妇av| 免费看a级黄色片| 特级一级黄色大片| 亚洲不卡免费看| 国产淫片久久久久久久久| 亚洲精品在线观看二区| 久久99热6这里只有精品| 亚洲国产精品成人综合色| 搡老岳熟女国产| 婷婷精品国产亚洲av| 中文字幕人妻熟人妻熟丝袜美| 在现免费观看毛片| 乱人视频在线观看| 国产精品人妻久久久影院| 亚洲第一区二区三区不卡| 在线播放国产精品三级| 久9热在线精品视频| 免费看美女性在线毛片视频| 人人妻人人澡欧美一区二区| av.在线天堂| 国产精品人妻久久久久久| 日本黄色片子视频| 性色avwww在线观看| 亚洲av日韩精品久久久久久密| 国产精品久久久久久精品电影| 国产视频一区二区在线看| 成人综合一区亚洲| 一级av片app| 嫁个100分男人电影在线观看| 亚洲内射少妇av| 国产精品无大码| 人妻制服诱惑在线中文字幕| 尤物成人国产欧美一区二区三区| 日韩 亚洲 欧美在线| 亚洲aⅴ乱码一区二区在线播放| 久久中文看片网| 午夜福利在线在线| 久久久久国产精品人妻aⅴ院| 国产aⅴ精品一区二区三区波| 少妇人妻精品综合一区二区 | 久久久久国内视频| 久久婷婷人人爽人人干人人爱| 一个人看视频在线观看www免费| 高清在线国产一区| 乱系列少妇在线播放| 少妇猛男粗大的猛烈进出视频 | 午夜免费成人在线视频| 国产黄片美女视频| 中文字幕高清在线视频| 国产av一区在线观看免费| 桃色一区二区三区在线观看| 国产高清视频在线观看网站| 天天一区二区日本电影三级| 午夜老司机福利剧场| 男人舔奶头视频| 亚洲avbb在线观看| 精品午夜福利在线看| 内射极品少妇av片p| 在现免费观看毛片| 少妇丰满av| 日日摸夜夜添夜夜添小说| 少妇猛男粗大的猛烈进出视频 | 国产av在哪里看| 亚洲在线自拍视频| а√天堂www在线а√下载| 三级国产精品欧美在线观看| 黄色欧美视频在线观看| 欧美另类亚洲清纯唯美| x7x7x7水蜜桃| 禁无遮挡网站| 男女视频在线观看网站免费| 91久久精品国产一区二区三区| 亚洲精品色激情综合| 黄色女人牲交| 亚洲美女视频黄频| 一a级毛片在线观看| 少妇人妻一区二区三区视频| 我的老师免费观看完整版| 熟女电影av网| 国内精品一区二区在线观看| 亚洲无线观看免费| 真人做人爱边吃奶动态| 99国产极品粉嫩在线观看| 国产老妇女一区| 国产精品1区2区在线观看.| 日日摸夜夜添夜夜添av毛片 | 村上凉子中文字幕在线| 国内少妇人妻偷人精品xxx网站| 成人av一区二区三区在线看| 免费观看在线日韩| 国产精品98久久久久久宅男小说| 欧美激情久久久久久爽电影| 国产精品电影一区二区三区| 国产又黄又爽又无遮挡在线| 内地一区二区视频在线| xxxwww97欧美| 久久久久性生活片| 夜夜夜夜夜久久久久| 一本一本综合久久| 成人无遮挡网站| 听说在线观看完整版免费高清| 亚洲av美国av| 久久精品国产自在天天线| 亚洲真实伦在线观看| 欧美一区二区精品小视频在线| 久久香蕉精品热| 免费无遮挡裸体视频| 中文字幕av在线有码专区| 日本三级黄在线观看| 女人被狂操c到高潮| 少妇人妻精品综合一区二区 | 国产精品嫩草影院av在线观看 | 国产淫片久久久久久久久| 99国产精品一区二区蜜桃av| 亚洲最大成人av| 国产精品人妻久久久影院| 香蕉av资源在线| 丰满乱子伦码专区| 我要看日韩黄色一级片| 国产精品久久久久久精品电影| 精品久久久久久,| 国产乱人视频| 波多野结衣高清无吗| 亚洲精品在线观看二区| 国产三级中文精品| 欧美又色又爽又黄视频| 亚洲精品在线观看二区| 最好的美女福利视频网| 国产精品福利在线免费观看| 国产男靠女视频免费网站| 久久热精品热| 91午夜精品亚洲一区二区三区 | 少妇猛男粗大的猛烈进出视频 | 一a级毛片在线观看| 国产三级在线视频| 亚洲18禁久久av| 国产精品日韩av在线免费观看| 少妇高潮的动态图| 1000部很黄的大片| 日韩中字成人| 日本五十路高清| 亚洲色图av天堂| 91午夜精品亚洲一区二区三区 | 欧美一区二区亚洲| 97热精品久久久久久| 99在线视频只有这里精品首页| 51国产日韩欧美| 欧美极品一区二区三区四区| 国产激情偷乱视频一区二区| 亚洲性久久影院| 人妻夜夜爽99麻豆av| 自拍偷自拍亚洲精品老妇| 精品乱码久久久久久99久播| 国国产精品蜜臀av免费| 狠狠狠狠99中文字幕| 国产精品嫩草影院av在线观看 | 看黄色毛片网站| 国产视频一区二区在线看| 18禁黄网站禁片免费观看直播| 少妇人妻一区二区三区视频| 亚洲欧美激情综合另类| 国产精品久久视频播放| 亚洲七黄色美女视频| 亚洲av熟女| 啪啪无遮挡十八禁网站| 中国美女看黄片| 中文字幕熟女人妻在线| 免费大片18禁| 男女之事视频高清在线观看| 国产精品久久久久久av不卡| 午夜福利在线观看吧| 久久久久久久亚洲中文字幕| 哪里可以看免费的av片| 久久精品国产清高在天天线| 国产91精品成人一区二区三区| 国产成年人精品一区二区| 亚洲av日韩精品久久久久久密| 嫩草影院精品99| 嫩草影院新地址| 久久午夜福利片| 亚洲国产精品成人综合色| 色精品久久人妻99蜜桃| 国产精华一区二区三区| 国产人妻一区二区三区在| 亚洲欧美激情综合另类| 美女高潮的动态| 欧美日韩亚洲国产一区二区在线观看| 两人在一起打扑克的视频| 亚洲图色成人| 欧美潮喷喷水| 成人高潮视频无遮挡免费网站| 日韩在线高清观看一区二区三区 | 亚洲avbb在线观看| 国内精品美女久久久久久| 性色avwww在线观看| 亚洲av熟女| av视频在线观看入口| av视频在线观看入口| 美女cb高潮喷水在线观看| 亚洲欧美激情综合另类| 日本三级黄在线观看| 日韩一区二区视频免费看| 三级国产精品欧美在线观看| 婷婷六月久久综合丁香| 国产一区二区在线av高清观看| 免费看美女性在线毛片视频| 夜夜爽天天搞| 国产成人aa在线观看| 欧美高清性xxxxhd video| 日韩欧美国产在线观看| 日韩在线高清观看一区二区三区 | 精品人妻偷拍中文字幕| 精品国产三级普通话版| 女生性感内裤真人,穿戴方法视频| 嫩草影院新地址| 亚洲最大成人手机在线| 婷婷丁香在线五月| 国产精品电影一区二区三区| h日本视频在线播放| 精品人妻偷拍中文字幕| 99久久精品热视频| 最近视频中文字幕2019在线8| 日韩av在线大香蕉| 男人和女人高潮做爰伦理| 老司机午夜福利在线观看视频| 男女那种视频在线观看| 欧美日韩国产亚洲二区| 黄色日韩在线| 悠悠久久av| 成人国产麻豆网| 色哟哟·www| 午夜福利在线观看吧| 国产欧美日韩一区二区精品| 国产精品伦人一区二区| 国产大屁股一区二区在线视频| 国产欧美日韩精品一区二区| 婷婷亚洲欧美| 1000部很黄的大片| 精华霜和精华液先用哪个| 国产av不卡久久| 欧美成人a在线观看| 国内精品久久久久久久电影| 国产久久久一区二区三区| 成人欧美大片| 白带黄色成豆腐渣| 91久久精品电影网| 97人妻精品一区二区三区麻豆| 97超视频在线观看视频| 欧美成人性av电影在线观看| 国产欧美日韩一区二区精品| 97超级碰碰碰精品色视频在线观看| 直男gayav资源| 成人永久免费在线观看视频| 国产亚洲精品久久久久久毛片| av在线观看视频网站免费| 91麻豆精品激情在线观看国产| 中文字幕熟女人妻在线| 欧美xxxx黑人xx丫x性爽| 啦啦啦啦在线视频资源| 国内精品久久久久精免费| 成人高潮视频无遮挡免费网站| 国产精品综合久久久久久久免费| 国产一区二区在线观看日韩| 欧美性猛交黑人性爽| 美女cb高潮喷水在线观看| 欧美日韩亚洲国产一区二区在线观看| 搡老熟女国产l中国老女人| 国产极品精品免费视频能看的| 成人午夜高清在线视频| 亚洲avbb在线观看| 亚洲专区国产一区二区| 亚洲一级一片aⅴ在线观看| 亚洲国产色片| 国产亚洲精品综合一区在线观看| 村上凉子中文字幕在线| 亚洲av熟女| 成人美女网站在线观看视频| 日韩精品青青久久久久久| 一个人免费在线观看电影| 国产精品永久免费网站| 亚洲中文字幕一区二区三区有码在线看| 变态另类丝袜制服| 日本黄色片子视频| 婷婷精品国产亚洲av| 亚洲国产欧洲综合997久久,| 热99re8久久精品国产| 一个人看视频在线观看www免费| 免费av毛片视频| 别揉我奶头~嗯~啊~动态视频| 伦精品一区二区三区| 婷婷六月久久综合丁香| 99精品久久久久人妻精品| 在线a可以看的网站| 日本免费a在线| 99久久精品一区二区三区| 88av欧美| 真实男女啪啪啪动态图| 色吧在线观看| 久久国产精品人妻蜜桃| 国产91精品成人一区二区三区| 校园春色视频在线观看| 精品人妻视频免费看| 亚洲欧美激情综合另类| 国产av在哪里看| 国产成人aa在线观看| 色哟哟哟哟哟哟| 国产成人影院久久av| 99久久九九国产精品国产免费| 免费看美女性在线毛片视频| 国产精品亚洲美女久久久| 综合色av麻豆| 波多野结衣巨乳人妻| 久久精品影院6| 男女那种视频在线观看| 88av欧美| 日韩大尺度精品在线看网址| 欧美精品国产亚洲| 99热精品在线国产| 有码 亚洲区| av在线老鸭窝| 在线免费观看的www视频| 日本三级黄在线观看| 久久精品国产自在天天线| 一级av片app| 人人妻人人澡欧美一区二区| 在线免费观看不下载黄p国产 | 国产中年淑女户外野战色| 国产精品综合久久久久久久免费| 亚洲成人久久爱视频| 日本黄大片高清| 禁无遮挡网站| 色哟哟哟哟哟哟| 免费观看的影片在线观看| 国产老妇女一区| 中国美女看黄片| .国产精品久久| 天天躁日日操中文字幕| av天堂中文字幕网| 中国美女看黄片| .国产精品久久| 亚洲精品影视一区二区三区av| 赤兔流量卡办理| 国产精品日韩av在线免费观看| 大又大粗又爽又黄少妇毛片口| 国产精品人妻久久久影院| 真实男女啪啪啪动态图| 波多野结衣高清作品| 精品久久久噜噜| 亚洲精品色激情综合| 又黄又爽又刺激的免费视频.| 色吧在线观看| 成人综合一区亚洲| 久久精品综合一区二区三区| 在线国产一区二区在线| 一级毛片久久久久久久久女| 国产免费男女视频| 丝袜美腿在线中文| 在线a可以看的网站| 久久久久国内视频| 午夜福利在线观看免费完整高清在 | 精品久久久久久久久久久久久| 国产伦精品一区二区三区视频9| 色av中文字幕| 丰满的人妻完整版| 亚洲在线自拍视频| 女生性感内裤真人,穿戴方法视频| 欧美潮喷喷水| 嫩草影院精品99| 91狼人影院| 偷拍熟女少妇极品色| 国产精品亚洲美女久久久| 免费观看人在逋| 不卡视频在线观看欧美| 久久久久久久精品吃奶| 亚洲专区国产一区二区| 久久九九热精品免费| 我的老师免费观看完整版| 别揉我奶头~嗯~啊~动态视频| 老熟妇仑乱视频hdxx| 一进一出抽搐动态| 亚洲美女黄片视频| 国产精品久久久久久久久免| 国内揄拍国产精品人妻在线| 国产在视频线在精品| 国产精品免费一区二区三区在线| 淫妇啪啪啪对白视频| 九色成人免费人妻av| av专区在线播放| 国产精品日韩av在线免费观看| 色综合婷婷激情| 美女xxoo啪啪120秒动态图| 欧美人与善性xxx| av视频在线观看入口| 黄色日韩在线| 成年人黄色毛片网站| 91久久精品电影网| 亚洲成a人片在线一区二区| 2021天堂中文幕一二区在线观| 国产精品伦人一区二区| 黄色女人牲交| 人妻丰满熟妇av一区二区三区| 国产精品三级大全| 亚洲av中文字字幕乱码综合| 国产国拍精品亚洲av在线观看| 又爽又黄a免费视频| 午夜a级毛片| 极品教师在线视频| 久久久成人免费电影| 精品人妻一区二区三区麻豆 | 欧美日韩中文字幕国产精品一区二区三区| 欧美日本视频| 丝袜美腿在线中文| 欧美三级亚洲精品| 国产中年淑女户外野战色| 国产精品无大码| 亚洲男人的天堂狠狠| 久久婷婷人人爽人人干人人爱| 观看免费一级毛片| 国产一区二区亚洲精品在线观看| 欧美激情国产日韩精品一区| 中亚洲国语对白在线视频| 日韩中文字幕欧美一区二区| 能在线免费观看的黄片| 亚洲乱码一区二区免费版| 欧美极品一区二区三区四区| 国产精品免费一区二区三区在线| 国产大屁股一区二区在线视频| 一个人看的www免费观看视频| 欧美黑人巨大hd| 在线a可以看的网站| 午夜精品在线福利| av中文乱码字幕在线| 国产精品久久久久久久久免| 非洲黑人性xxxx精品又粗又长| 欧美日韩黄片免| 嫁个100分男人电影在线观看| 国产精品1区2区在线观看.| 中文字幕久久专区| 国产乱人视频| 禁无遮挡网站| 在线观看66精品国产| 欧美性猛交黑人性爽| 黄色欧美视频在线观看| 久久久成人免费电影| 毛片女人毛片| 波野结衣二区三区在线| 日本爱情动作片www.在线观看 | 国产美女午夜福利| 在线观看舔阴道视频| 国产精品三级大全| 精品久久久久久久久久免费视频| 亚洲精品成人久久久久久| 亚洲国产色片| 一区二区三区激情视频| 丝袜美腿在线中文| 两性午夜刺激爽爽歪歪视频在线观看| 久久久久久大精品| 精品人妻一区二区三区麻豆 | 久久精品91蜜桃| 国产精品久久久久久av不卡| 国产欧美日韩精品一区二区| 久久中文看片网| 成人亚洲精品av一区二区| ponron亚洲| 亚洲四区av| 不卡视频在线观看欧美| 丰满人妻一区二区三区视频av| 日本在线视频免费播放| 一a级毛片在线观看| 国产v大片淫在线免费观看| av中文乱码字幕在线| 少妇人妻精品综合一区二区 | 制服丝袜大香蕉在线| 联通29元200g的流量卡| 99久久中文字幕三级久久日本| 啦啦啦啦在线视频资源| 97碰自拍视频| 人妻少妇偷人精品九色| av.在线天堂| 中亚洲国语对白在线视频| 黄色配什么色好看| 五月玫瑰六月丁香| 观看免费一级毛片| avwww免费| 精品无人区乱码1区二区| 男人舔奶头视频| 精品无人区乱码1区二区| 亚洲精品在线观看二区| 国产视频内射| av在线观看视频网站免费| 免费av毛片视频| 性欧美人与动物交配| 男人舔奶头视频| 非洲黑人性xxxx精品又粗又长| 午夜福利高清视频| 可以在线观看的亚洲视频| 嫩草影视91久久| 欧美色欧美亚洲另类二区| 日本与韩国留学比较| 麻豆av噜噜一区二区三区| 2021天堂中文幕一二区在线观| 欧美3d第一页| 国产一区二区三区av在线 | 色5月婷婷丁香| 国产一区二区在线观看日韩| 少妇裸体淫交视频免费看高清| 欧美高清成人免费视频www| 国产一区二区在线av高清观看| 亚洲人成网站高清观看| 日韩欧美三级三区| 国产aⅴ精品一区二区三区波| av在线观看视频网站免费| 亚洲人与动物交配视频| 1000部很黄的大片| 好男人在线观看高清免费视频| 精品无人区乱码1区二区| 99视频精品全部免费 在线| 啪啪无遮挡十八禁网站| 国产精品国产高清国产av| 亚洲七黄色美女视频| 国产毛片a区久久久久| 日本a在线网址| 国产精品一区二区三区四区免费观看 | 国语自产精品视频在线第100页| 99久久无色码亚洲精品果冻| 一级黄色大片毛片| 午夜福利在线观看吧| 夜夜看夜夜爽夜夜摸| 久久久国产成人免费| 亚洲精品亚洲一区二区| aaaaa片日本免费| 深夜a级毛片| 精品日产1卡2卡| 久久精品国产亚洲av香蕉五月| 久久人妻av系列| 有码 亚洲区| 国产探花在线观看一区二区| 国产伦精品一区二区三区视频9| 最新在线观看一区二区三区| 日韩欧美 国产精品| 日韩欧美精品免费久久| 中文资源天堂在线| 亚洲内射少妇av| 在线观看午夜福利视频| 国产欧美日韩一区二区精品| 亚洲国产日韩欧美精品在线观看| 亚洲精品色激情综合| 精品午夜福利视频在线观看一区| 免费在线观看影片大全网站| 欧美成人性av电影在线观看| 日日摸夜夜添夜夜添小说| 久久99热这里只有精品18| 国产精品伦人一区二区| 午夜精品久久久久久毛片777| 国产伦一二天堂av在线观看| 哪里可以看免费的av片| 又爽又黄a免费视频| 91av网一区二区| 非洲黑人性xxxx精品又粗又长| 国产乱人视频| 性欧美人与动物交配| 国产高清激情床上av| 十八禁国产超污无遮挡网站| 亚洲内射少妇av| 91在线观看av| 国产高清有码在线观看视频| 黄色日韩在线| 久久久国产成人免费| 性欧美人与动物交配| 91av网一区二区| 别揉我奶头 嗯啊视频| 丰满人妻一区二区三区视频av| 窝窝影院91人妻| 欧美性感艳星| 亚洲最大成人中文| 欧美黑人巨大hd| 不卡一级毛片| 亚洲欧美激情综合另类| 久久精品国产自在天天线| 日韩欧美国产在线观看| 男女边吃奶边做爰视频| 亚洲最大成人中文| 男人的好看免费观看在线视频| 成人国产一区最新在线观看| 老司机福利观看|