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

    LaBr3(Ce)探測器探測效率影響因素研究

    2022-01-20 00:56:56陳志強陳晶晶舒雙寶余子樵張育中郎賢禮
    核技術(shù) 2022年1期
    關(guān)鍵詞:光子射線晶體

    陳志強,陳晶晶,舒雙寶,余子樵,張育中,郎賢禮

    (合肥工業(yè)大學(xué)儀器科學(xué)與光電工程學(xué)院合肥230009)

    LaBr3閃爍體探測器是21世紀(jì)初開發(fā)的高性能輻射探測器,由于其優(yōu)異的性能引起了研究者的廣泛關(guān)注。從LaBr3晶體的最初發(fā)現(xiàn)和科研人員的發(fā)展,其性能得到了顯著的提高,LaBr3閃爍體探測器已被廣泛應(yīng)用于許多領(lǐng)域[1]。已有研究報道,摻雜一定比例的金屬離子可以使LaBr3晶體具有更高的能量分辨率和更好的探測性能。其中摻雜一定比例的Ce3+晶體具有較高的能量分辨率和較好的衰減時間系數(shù)[2-3]。本文中的LaBr3(Ce)晶體是指含5%Ce3+的LaBr3晶體。

    表1為幾種常見的無機閃爍晶體探測器性能參數(shù)比較,對于LaBr3(Ce)晶體,由表1可以看出,與其他無機閃爍體材料相比,LaBr3(Ce)晶體在一些方面具有明顯優(yōu)勢。LaBr3(Ce)閃爍晶體的突出特性包括優(yōu)異的能量分辨率(約3%@662 keV γ射線)、高光產(chǎn)額(約63光子/keV)、衰減時間短(約20 ns)以及良好的線性響應(yīng)(在60~1 300 keV能量范圍小于7%)[4-7],當(dāng)性能要求較高時,LaBr3(Ce)通常廣泛用作NaI(TI)晶體的替代品[8]。但是這種晶體在空氣中容易潮解,難以生長和加工出大尺寸的晶體[7,9]。目前,獲得大尺寸LaBr3(Ce)晶體的主要方法是Bridgeman-Stockbarge法[10]。LaBr3大尺寸晶體的發(fā)展使其商業(yè)應(yīng)用成為現(xiàn)實[11]。LaBr3(Ce)晶體具有良好的性能,已成為許多領(lǐng)域的常用材料,包括核測井[12-13]、環(huán)境監(jiān)測[14-15]、核輻射防護[16],以及核醫(yī)學(xué)等[17-18]。

    由于LaBr3(Ce)晶體是一種極具發(fā)展?jié)摿Φ臒o機閃爍體,鐘丁生等針對幾何因素對LaBr3晶體的探測效率進行了研究[1],發(fā)現(xiàn)不斷增大溴化鑭晶體的幾何尺寸,探測效率會明顯非線性增大,但是并沒有對晶體形狀、能量等因素進行分析。本文增加了晶體形狀、能量因素等對探測效率影響的研究。本文以LaBr3(Ce)閃爍體探測器為研究對象分析了影響探測器探測效率的因素。根據(jù)定義,射線在閃爍體中產(chǎn)生的閃爍脈沖數(shù)NC與入射粒子N的比值稱為閃爍體對射線的探測效率,即C=NC/N。

    本文通過使用Geant4工具箱建立LaBr3(Ce)閃爍體探測器的探測模型,研究影響LaBr3閃爍體探測效率的因素。介紹了Geant4模型的構(gòu)建過程以及模擬過程中所需要的物理過程。對模擬數(shù)據(jù)的處理分析,并給出了晶體的形狀、尺寸以及射線的能量和放射源位置與探測效率的關(guān)系。

    1 Geant4模擬實驗步驟

    計算機程序模擬實驗可以大大降低實驗成本,減少不必要的時間浪費,使研究人員能夠立即從事其他重要方面的研究工作。例如,文獻[19]采用MCNP(Monte Carlo N Particle Transport Code)代碼進行城市范圍內(nèi)輻射傳輸?shù)母弑U婺M。Geant4工具包是歐洲核研究組織基于C++面向?qū)ο笳Z言開發(fā)的蒙特卡羅應(yīng)用軟件工具包。經(jīng)過多年的仿真實踐,Geant4的仿真可靠性得到了很好的驗證,使得Geant4工具箱在許多領(lǐng)域得到了廣泛的應(yīng)用[20]。比如,利用Geant4工具箱進行缺陷類型的后向散射分析具有良好的應(yīng)用前景[21];Geant4蒙特卡羅模擬解決了低能γ能譜模擬實驗[22]。

    本研究是基于Geant4代碼構(gòu)建的模型。在圖1中的框圖展示了Geant4工具箱的具體組成部分,包括粒子源、物理幾何、信息獲取與統(tǒng)計及物理過程等。Geant4模擬的具體過程如下:粒子發(fā)生器產(chǎn)生的隨機射線粒子被輻射到預(yù)先構(gòu)造的物理幾何中,射線粒子與物理幾何材料發(fā)生相應(yīng)的物理過程。最后,得到用戶需要的粒子的相關(guān)信息。在這里,我們主要介紹物理幾何模型的構(gòu)建和物理過程的選擇兩部分。

    圖1 Geant4的模擬過程Fig.1 Diagram of simulation process of Geant4

    1.1 模型構(gòu)建

    在Geant4工具包的基礎(chǔ)上,通過C++編程構(gòu)建了LaBr3(Ce)探測器的幾何模型。具體模型結(jié)構(gòu)如圖2所示。整個模型置于真空環(huán)境中,入射射線(輻射源)的起始位置與探測器軸線在同一條直線上。假設(shè)入射粒子為662 keV單能γ射線束,事件過程大致過程為:首先,入射射線沿軸線發(fā)出,入射γ粒子打到探測器前端面的中心處。然后,射線粒子和探測器材料相互作用,經(jīng)歷相應(yīng)的物理過程包括康普頓散射、光電效應(yīng)、閃爍等。在這些過程中射線粒子沉積能量,并在閃爍體中會產(chǎn)生相應(yīng)數(shù)目的閃爍光子。閃爍光子在晶體中傳播可能發(fā)生多次折射和反射,最后從閃爍體出射的光子被光電倍增管的光電陰極收集。最后,程序收集光子數(shù)量和其他相關(guān)數(shù)據(jù)。在Geant4程序中是利用粒子與光電陰極之間的反應(yīng)來確定光子是否被記錄,并獲得相應(yīng)的光子數(shù)據(jù)。如果放射源是各向同性γ源,所不同的是該源的輻射面是球狀,放射的γ射線束與探測器的前端面形成一個錐角[23],也就是理論上放射的射線粒子可以到達探測器前端的任意位置。通過上述公式處理得到的數(shù)據(jù)就可以得到探測器的絕對探測效率。

    圖2 LaBr3(Ce)探測器的Geant4模型Fig.2 Geant4 model of the LaBr3(Ce)detector

    LaBr3(Ce)探測器的結(jié)構(gòu)如圖2模型的右側(cè)所示。該模型結(jié)構(gòu)外部有一層鋁涂層,中間是MgO反射層,然后是LaBr3(Ce)晶體,其后端通過光學(xué)玻璃與PMT光陰極進行光學(xué)耦合(從左到右排列)。在設(shè)置模型數(shù)據(jù)時,我們參考之前文獻中LaBr3(Ce)探測器的相關(guān)信息設(shè)置模型的結(jié)構(gòu)尺寸和相關(guān)物理參數(shù)[24-25]。由于需要模擬不同截面的LaBr3(Ce)晶體,因此模擬晶體的橫截面考慮了以下三種形狀:圓形、正方形和長方形(長邊是短邊的兩倍)。在Geant4工具箱中,這三種結(jié)構(gòu)的探測器模型如圖3(a)、(b)、(c)所示。在對晶體尺寸進行仿真研究時,需要考慮半徑和長度變化對仿真的影響。因此,Geant4模型中晶體的長度也需要設(shè)定以滿足研究的需要。

    圖3 不同截面形狀的探測器幾何模型Fig.3 Geometric models of different cross-sectional shapes

    1.2 反應(yīng)過程及參數(shù)設(shè)定

    粒子的反應(yīng)過程主要考慮電磁相互作用、閃爍光透射等物理過程。這些物理過程可以在Geant4程序中提供的類中獲得。因此,在選擇物理過程時,需要調(diào)用Geant4內(nèi)置函數(shù)庫中標(biāo)準(zhǔn)物理過程的G4EmStandardPhysics類和光學(xué)反應(yīng)過程的G4OpticalPhysics類,并進行相應(yīng)的設(shè)置。具體設(shè)置包括確定粒子截斷值為0.1 mm,對物理過程進行注冊以及構(gòu)建物理過程需要構(gòu)建這些過程中所需的粒子類型,如電子、光子等。

    用于光學(xué)反應(yīng)過程的G4OpticalPhysics類需要設(shè)置與閃爍體相關(guān)的光學(xué)參數(shù)。設(shè)置LaBr3(Ce)閃爍體的折射率為1.9;光產(chǎn)額率為63光子·keV-1;快時間常數(shù)是16。石英玻璃窗的折射率為1.47。MgO表面反射率將波瓣反射設(shè)為0.94,鏡面反射設(shè)為0.03,漫反射設(shè)為0.03,物質(zhì)之間的邊界類型設(shè)為 dielectric_metal。

    2 結(jié)果和討論

    2.1 截面形狀對探測效率的影響

    目前,閃爍探測器常用的LaBr3(Ce)晶體截面為圓形。本研究通過模擬不同形狀的晶體截面,分析截面形狀對探測器探測效率的影響。為了實現(xiàn)這一目標(biāo),本文設(shè)置了三個具有代表性的截面,分別是圓形截面、方形截面和矩形截面。具體模型結(jié)構(gòu)如圖3所示。設(shè)置矩形截面來解釋x軸和y軸尺寸差異對探測效率的影響。同時,為簡化仿真程序和設(shè)計的復(fù)雜性,在設(shè)定矩形邊長時,將矩形的長邊長定為短邊長的兩倍。為了消除實驗中其他因素的干擾,采用變量控制方法使三種情況的晶體截面積相等。

    模擬過程的射線粒子由Geant4程序中的粒子發(fā)生器進行生成,共模擬105個事件。由此產(chǎn)生的662 keV單能γ射線粒子從起始點穿過探測器前端中心沿探測器軸線傳播。入射到LaBr3(Ce)晶體中的射線粒子與探測介質(zhì)相互作用并沉積能量,從而產(chǎn)生相應(yīng)的粒子,如閃爍光子。模擬程序根據(jù)用戶選擇的數(shù)據(jù)采集條件采集光電陰極的粒子數(shù)據(jù),并對采集的數(shù)據(jù)信息進行處理。LaBr3(Ce)晶體三種不同截面的探測效率仿真結(jié)果如圖4所示,其中晶體的面積范圍為52π~802π mm2,面積對應(yīng)的尺寸間隔為5 mm。圖中橫坐標(biāo)圓形截面所對應(yīng)的半徑來進行標(biāo)注,縱坐標(biāo)表示探測效率。

    圖4 LaBr3(Ce)閃爍體探測器在不同截面形狀下的探測效率Fig.4 Detection efficiency of LaBr3(Ce)crystals detector with different cross-section shapes

    從圖4中可以看出,三種形狀截面晶體的探測效率曲線都有一個共同的特點,即探測效率隨晶體截面面積的增加而增加并且增速在不斷減小。但是,從圖4可以看出,圓形晶體截面與正方形晶體截面代表的探測效率直線幾乎重合,說明探測效率差異不顯著,且這兩者都比矩形截面閃爍晶體的探測效率高。產(chǎn)生這個結(jié)果的原因是三個截面面積相同,但由于矩形在x軸和y軸方向上晶體的尺寸差異較大,放射源進行放射時,由于達到截面的粒子密度越靠近中心越高,所以矩形截面晶體接收的粒子數(shù)較少并且射線在閃爍體中的運動在較短的軸方向上受到限制,這些粒子不能順利達到光陰極產(chǎn)生有效光子計數(shù),導(dǎo)致矩形截面LaBr3(Ce)閃爍晶體的探測效率較低。此外,從圖4中可以看出,隨著閃爍體半徑的增大,矩形截面LaBr3(Ce)閃爍體晶體與其他兩種閃爍體的探測效率差異逐漸增大。因此,若想要得到更好的探測效率,可以考慮使用圓形或正方形截面的LaBr3(Ce)閃爍體晶體。

    2.2 LaBr3(Ce)晶體幾何尺寸對探測效率的影響

    從以上結(jié)果可以看出,圓形截面是閃爍體的最佳截面形狀,具有較高的探測效率。因此,本節(jié)使用一個具有圓形截面的閃爍體模型??芍?,圓柱形晶體的尺寸大小由兩個方面長度組成:閃爍體的橫截面半徑和閃爍體的長度。在模擬過程中,采用控制變量方法模擬了105個事件。根據(jù)尺寸要求,我們的半徑范圍為5~40 mm,尺寸間隔為5 mm;同樣,長度范圍為10~80 mm,間隔10 mm。這就得到了8×8組數(shù)據(jù)。然后根據(jù)探測效率計算公式計算這些數(shù)據(jù),得到相應(yīng)的探測效率,數(shù)據(jù)處理結(jié)果如圖5所示。

    圖5 探測效率與閃爍體長度(a)和半徑(b)的關(guān)系Fig.5 The relationship between detection efficiency and length(a)and radius(b)of the scintillator

    從圖5(a)中可以看出,整體上隨著LaBr3(Ce)閃爍體長度的增加,探測效率逐漸增加。閃爍體半徑與探測效率的關(guān)系如圖5(b)所示。隨著閃爍體半徑的增大,閃爍體的探測效率也增大,半徑的變化對探測效率的影響比長度更顯著。本文得到的數(shù)據(jù)與文獻[1]結(jié)果吻合。圖6為同時考慮長度和半徑尺寸對探測效率的影響的二維曲面圖??梢酝茢?,LaBr3(Ce)閃爍體的探測效率取決于閃爍體的尺寸,且閃爍體的尺寸越大,其探測效率越高。這種可能是因為晶體的尺寸越大,與閃爍體材料在閃爍體中的輻射反應(yīng)時間越長,將增加產(chǎn)生的閃爍光子的數(shù)量。但增加閃爍體的長度并不能顯著提高探測效率,一方面是因為輻射與閃爍體材料的反應(yīng)時間增加;另一方面,由于傳輸距離的增加,部分反應(yīng)產(chǎn)生的粒子在閃爍體中湮沒,造成了在閃爍體中無法達到的光電陰極的損失,從而減緩了增長。

    圖6 探測效率與晶體的半徑和長度的關(guān)系Fig.6 The relationship between detection efficiency and the combination of radius and length of LaBr3(Ce)crystals

    2.3 射線能量與閃爍效率的關(guān)系

    先前的研究表明,輻射的能量也會影響探測器的輸出響應(yīng)。那么,輻射能量的大小是否也會影響晶體的探測效率。為了驗證這一假設(shè),在模擬中使用了能量范圍為0.1~1 MeV(能量區(qū)間為0.1 MeV)的單能γ射線LaBr3(Ce)閃爍體探測器,共105個事件。探測器模型仍然采用30×15 mm2的圓形截面結(jié)構(gòu)。收集并處理數(shù)據(jù)結(jié)果,數(shù)據(jù)結(jié)果如圖7所示。

    圖7 探測效率與γ射線能量的關(guān)系Fig.7 The relationship between the detection efficiency and the energy of γ-ray

    在圖7中,LaBr3(Ce)閃爍體的探測效率隨著輻射能量的增加而降低,然后隨能量的增加而趨于穩(wěn)定。因為當(dāng)使用單能量γ射線束進行探測時,射線的能量越低,與閃爍體材料充分反應(yīng)的粒子越多,產(chǎn)生的閃爍光子也就越多。由探測效率公式計算的探測效率可以看出,探測效率隨入射能量的增加而降低。因此,可以看出LaBr3(Ce)探測器用來探測能量較低的γ射線輻射源的光譜,可以取得較好的結(jié)果。

    2.4 輻射位置對探測效率的影響

    在上述三個部分,我們討論都是對單一能量γ射線束的探測,但是考慮到在現(xiàn)實世界中探測到輻射的來源時的情況,發(fā)散的放射源發(fā)出的位置對探測效果有一定的影響。在文獻[26]中就討論了探測效率與污染表面之間的距離和表面形狀之間的關(guān)系。因此,我們用Geant4工具包模擬了LaBr3(Ce)探測器對各向同性γ輻射源的探測。模型結(jié)構(gòu)參照上節(jié)模型,輻射能量為662 keV。然后,在每次模擬開始之前,對放射源和探測器之間的距離進行變化。距離范圍為5~100 mm,間隔5 mm。得到20組數(shù)據(jù),結(jié)果如圖8所示。

    圖8 探測效率與放射源位置的關(guān)系Fig.8 The relationship between the detection efficiency and the position of the radioactive source

    從圖8的數(shù)據(jù)結(jié)果可以看出,LaBr3(Ce)閃爍體的探測效率確實與反射鏡與閃爍體之間的距離有關(guān)。閃爍體的探測效率隨著距離的增加而降低。從圖中還可以看出,當(dāng)距離大于70 mm時,探測效率逐漸趨于零。原因是輻射源發(fā)散得越遠,射線越發(fā)散,擊中探測器前端的射線就越稀疏。因此,收集得到的脈沖數(shù)量較少,導(dǎo)致探測效率降低。當(dāng)放射源與閃爍體之間的距離過大時,探測效率過低,會影響探測精度。

    3 結(jié)語

    本文基于Geant4工具包利用蒙特卡羅方法研究了LaBr3(Ce)閃爍體的探測效率。研究發(fā)現(xiàn):LaBr3(Ce)閃爍體的探測效率與晶體截面形狀、晶體尺寸、輻射源能量大小和位置有關(guān)。當(dāng)閃爍體面積相同時,正方形和圓形閃爍體的探測效率高于矩形截面的探測效率。同時,閃爍體尺寸越大,其探測效率越高。從能量角度可知,入射輻射能量越高,LaBr3(Ce)閃爍體的探測效率反而越低。此外,用各向同性γ粒子源模擬放射源時,發(fā)現(xiàn)閃爍體的探測效率隨著距離的增加而降低。因此,在設(shè)計LaBr3(Ce)閃爍體探測器時,如果想獲得更好的探測效率,可以考慮使用更大的晶體。同時,對輻射源的探測不能將探測器放置得離輻射源太遠,這樣會導(dǎo)致探測效率低,造成誤差。通過本次模擬實驗,我們可以從以上幾個方面了解LaBr3(Ce)閃爍體探測效率的影響因素及相關(guān)規(guī)律,從而更好地了解LaBr3(Ce)晶體特性,這對探測器的實際應(yīng)用和閃爍體的外觀設(shè)計具有指導(dǎo)意義,有利于在進行傳感器設(shè)計時選擇合適的晶體尺寸避免因尺寸不匹配而造成的設(shè)計誤差從而影響儀器精度,達到節(jié)約生產(chǎn)成本和時間成本的目的。本文模擬過程是在理想真空環(huán)境中進行的,沒有考慮真實環(huán)境中的宇宙射線等因素影響。在未來的工作中,這些因素將被考慮以達到全面地了解LaBr3(Ce)探測器探測效率的影響因素。

    作者貢獻聲明陳志強負(fù)責(zé)文章的起草及最終版本修訂;陳晶晶負(fù)責(zé)資料的收集及整理;舒雙寶負(fù)責(zé)研究的提出、設(shè)計及文章的修改;余子樵負(fù)責(zé)資料的收集及整理;張育中負(fù)責(zé)資料的收集,郎賢禮負(fù)責(zé)論文的修改。

    猜你喜歡
    光子射線晶體
    《光子學(xué)報》征稿簡則
    “輻射探測晶體”專題
    “直線、射線、線段”檢測題
    『直線、射線、線段』檢測題
    赤石脂X-射線衍射指紋圖譜
    中成藥(2017年3期)2017-05-17 06:09:16
    在光子帶隙中原子的自發(fā)衰減
    光子晶體在兼容隱身中的應(yīng)用概述
    多光子Jaynes-Cummings模型中與Glauber-Lachs態(tài)相互作用原子的熵壓縮
    放大鏡下的晶體
    把晶體找出來
    国产男人的电影天堂91| 校园人妻丝袜中文字幕| 午夜日本视频在线| 人妻 亚洲 视频| 999久久久国产精品视频| 啦啦啦视频在线资源免费观看| 搡女人真爽免费视频火全软件| 天堂中文最新版在线下载| 男女午夜视频在线观看| 黄色视频在线播放观看不卡| 国产视频首页在线观看| 国产成人精品无人区| 国产精品国产av在线观看| 国产深夜福利视频在线观看| 99久国产av精品国产电影| 亚洲综合色网址| 国产精品偷伦视频观看了| 母亲3免费完整高清在线观看 | 国产老妇伦熟女老妇高清| 少妇的丰满在线观看| 亚洲欧美一区二区三区黑人 | av免费在线看不卡| 日本-黄色视频高清免费观看| 精品一区二区三卡| 麻豆精品久久久久久蜜桃| 日韩免费高清中文字幕av| 国产精品 欧美亚洲| 十分钟在线观看高清视频www| 一级,二级,三级黄色视频| 精品人妻在线不人妻| 蜜桃国产av成人99| 一本久久精品| 国产视频首页在线观看| 亚洲国产av影院在线观看| 伊人亚洲综合成人网| 九草在线视频观看| 男女无遮挡免费网站观看| 亚洲精品中文字幕在线视频| 好男人视频免费观看在线| 青春草亚洲视频在线观看| 久久精品久久久久久噜噜老黄| 久久精品国产亚洲av高清一级| 黄片小视频在线播放| 国产精品香港三级国产av潘金莲 | 国产在线视频一区二区| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 亚洲激情五月婷婷啪啪| 国产黄色免费在线视频| 久久99一区二区三区| 日本免费在线观看一区| 熟女电影av网| 在线看a的网站| 天天躁夜夜躁狠狠久久av| 伊人亚洲综合成人网| 成人毛片a级毛片在线播放| 免费观看在线日韩| 久久国产精品大桥未久av| 国产一区二区激情短视频 | 午夜福利,免费看| 人妻一区二区av| 亚洲,欧美精品.| 精品一区二区三卡| 亚洲第一av免费看| 韩国精品一区二区三区| 欧美亚洲日本最大视频资源| 亚洲国产色片| 嫩草影院入口| 国产免费一区二区三区四区乱码| 国产精品熟女久久久久浪| 好男人视频免费观看在线| 美女xxoo啪啪120秒动态图| 国产男人的电影天堂91| 日本vs欧美在线观看视频| 欧美人与性动交α欧美精品济南到 | 韩国av在线不卡| 国产片特级美女逼逼视频| 国产免费视频播放在线视频| 久久人人爽av亚洲精品天堂| 国产精品三级大全| 精品午夜福利在线看| 中文字幕人妻熟女乱码| 一区在线观看完整版| 视频在线观看一区二区三区| 黑丝袜美女国产一区| 在线精品无人区一区二区三| 少妇的逼水好多| 精品视频人人做人人爽| 久久久a久久爽久久v久久| videos熟女内射| 久久精品国产亚洲av天美| 女的被弄到高潮叫床怎么办| 黄色视频在线播放观看不卡| 一级毛片我不卡| 国产成人免费无遮挡视频| 亚洲精品国产色婷婷电影| 少妇的丰满在线观看| 亚洲国产av影院在线观看| 国产黄频视频在线观看| 老鸭窝网址在线观看| 香蕉丝袜av| 免费观看在线日韩| 国产在线视频一区二区| 国产不卡av网站在线观看| 亚洲人成77777在线视频| 日韩,欧美,国产一区二区三区| 日韩在线高清观看一区二区三区| 国产成人精品一,二区| 久久精品国产鲁丝片午夜精品| 国产成人欧美| 久久人人爽人人片av| 国产 精品1| 色婷婷久久久亚洲欧美| 高清欧美精品videossex| 丝袜美足系列| 日韩成人av中文字幕在线观看| 中文字幕色久视频| 亚洲成国产人片在线观看| 一区二区av电影网| 国产精品欧美亚洲77777| 大片电影免费在线观看免费| 久久韩国三级中文字幕| a 毛片基地| 狠狠精品人妻久久久久久综合| 久久国产精品男人的天堂亚洲| 久久久久久久亚洲中文字幕| 精品亚洲成a人片在线观看| 波多野结衣av一区二区av| av线在线观看网站| 自线自在国产av| 97在线人人人人妻| 三上悠亚av全集在线观看| 在线免费观看不下载黄p国产| 亚洲精品日本国产第一区| 欧美国产精品va在线观看不卡| 久久韩国三级中文字幕| 免费在线观看视频国产中文字幕亚洲 | 男女啪啪激烈高潮av片| 欧美少妇被猛烈插入视频| 免费播放大片免费观看视频在线观看| 久久久久视频综合| av一本久久久久| 国产1区2区3区精品| 午夜老司机福利剧场| 日韩 亚洲 欧美在线| 黄色 视频免费看| 久久久国产欧美日韩av| 国产精品香港三级国产av潘金莲 | 七月丁香在线播放| 日本wwww免费看| 在线观看三级黄色| 亚洲精品国产色婷婷电影| 午夜老司机福利剧场| 亚洲欧美中文字幕日韩二区| 亚洲五月色婷婷综合| 日韩一区二区三区影片| 狠狠精品人妻久久久久久综合| 一个人免费看片子| 亚洲三级黄色毛片| 十八禁网站网址无遮挡| 伦理电影大哥的女人| 香蕉丝袜av| 只有这里有精品99| a级毛片在线看网站| 国产黄频视频在线观看| 欧美bdsm另类| 精品午夜福利在线看| av网站在线播放免费| 丰满迷人的少妇在线观看| 热99国产精品久久久久久7| 侵犯人妻中文字幕一二三四区| 韩国av在线不卡| 欧美成人精品欧美一级黄| 九九爱精品视频在线观看| 91成人精品电影| 在线观看一区二区三区激情| 成年女人在线观看亚洲视频| 亚洲第一青青草原| 国产在线视频一区二区| 久久久久久久久久人人人人人人| 久久精品国产亚洲av高清一级| 天天躁夜夜躁狠狠久久av| 久久99精品国语久久久| 欧美bdsm另类| av又黄又爽大尺度在线免费看| 国产亚洲午夜精品一区二区久久| 成人18禁高潮啪啪吃奶动态图| 国产黄色视频一区二区在线观看| 欧美亚洲日本最大视频资源| 国产又爽黄色视频| 国产男女超爽视频在线观看| 国产麻豆69| 黑人欧美特级aaaaaa片| 麻豆乱淫一区二区| 亚洲精品一区蜜桃| 97在线人人人人妻| 熟女av电影| 一级,二级,三级黄色视频| 亚洲av.av天堂| 少妇猛男粗大的猛烈进出视频| 人妻少妇偷人精品九色| 国产精品一二三区在线看| 久久精品aⅴ一区二区三区四区 | 制服丝袜香蕉在线| 国产 一区精品| 一本—道久久a久久精品蜜桃钙片| 久久久精品94久久精品| 一区二区三区四区激情视频| 国精品久久久久久国模美| 国语对白做爰xxxⅹ性视频网站| 亚洲四区av| 三级国产精品片| 97在线人人人人妻| 午夜福利乱码中文字幕| 亚洲国产精品999| 久久国产精品大桥未久av| 欧美xxⅹ黑人| 91午夜精品亚洲一区二区三区| 中文欧美无线码| 国产黄色免费在线视频| 中文欧美无线码| 熟女少妇亚洲综合色aaa.| 亚洲av免费高清在线观看| 伊人久久大香线蕉亚洲五| 久久久国产精品麻豆| 少妇被粗大猛烈的视频| 日本av手机在线免费观看| 欧美日韩亚洲高清精品| 久久毛片免费看一区二区三区| 99国产综合亚洲精品| 国产成人91sexporn| 国产精品久久久久久av不卡| 亚洲精品中文字幕在线视频| 午夜福利一区二区在线看| 99精国产麻豆久久婷婷| 26uuu在线亚洲综合色| 制服人妻中文乱码| 建设人人有责人人尽责人人享有的| 国产精品麻豆人妻色哟哟久久| 精品酒店卫生间| 日本色播在线视频| 日韩大片免费观看网站| 欧美精品av麻豆av| 美国免费a级毛片| 国产乱人偷精品视频| 美女高潮到喷水免费观看| 青草久久国产| 精品午夜福利在线看| 亚洲国产欧美网| 如何舔出高潮| 国产又色又爽无遮挡免| 国产有黄有色有爽视频| 国产欧美日韩综合在线一区二区| av免费在线看不卡| www.精华液| 日本欧美国产在线视频| 日韩制服丝袜自拍偷拍| 亚洲成av片中文字幕在线观看 | 亚洲精品久久久久久婷婷小说| 777米奇影视久久| 亚洲国产精品成人久久小说| videos熟女内射| 看免费成人av毛片| 18在线观看网站| 高清av免费在线| 美女国产高潮福利片在线看| 久久久精品免费免费高清| 国产极品天堂在线| 人人妻人人澡人人看| 黄色毛片三级朝国网站| 欧美精品亚洲一区二区| 纯流量卡能插随身wifi吗| 热99久久久久精品小说推荐| 美女高潮到喷水免费观看| 免费观看av网站的网址| 最黄视频免费看| 午夜久久久在线观看| videossex国产| 超碰97精品在线观看| 少妇的逼水好多| 免费观看在线日韩| 91aial.com中文字幕在线观看| 日韩三级伦理在线观看| 精品少妇久久久久久888优播| 免费在线观看视频国产中文字幕亚洲 | 日韩av不卡免费在线播放| 看非洲黑人一级黄片| √禁漫天堂资源中文www| 午夜日韩欧美国产| 在线免费观看不下载黄p国产| 国产亚洲av片在线观看秒播厂| 欧美黄色片欧美黄色片| 免费女性裸体啪啪无遮挡网站| 欧美日韩成人在线一区二区| 爱豆传媒免费全集在线观看| 看免费av毛片| 水蜜桃什么品种好| 999久久久国产精品视频| 久久久精品免费免费高清| 国产有黄有色有爽视频| 91国产中文字幕| 国产一区亚洲一区在线观看| 黄色毛片三级朝国网站| 韩国av在线不卡| 国产精品麻豆人妻色哟哟久久| 在线精品无人区一区二区三| 国产成人精品久久二区二区91 | 精品人妻在线不人妻| 国产女主播在线喷水免费视频网站| 色哟哟·www| 亚洲欧美一区二区三区久久| 女人被躁到高潮嗷嗷叫费观| 考比视频在线观看| 波野结衣二区三区在线| 亚洲色图 男人天堂 中文字幕| 日韩av免费高清视频| 精品少妇内射三级| 亚洲av综合色区一区| av视频免费观看在线观看| 狠狠婷婷综合久久久久久88av| 看非洲黑人一级黄片| 人人澡人人妻人| 亚洲欧美一区二区三区国产| 新久久久久国产一级毛片| 99re6热这里在线精品视频| 老鸭窝网址在线观看| av免费在线看不卡| 国产1区2区3区精品| 在线观看www视频免费| 午夜福利乱码中文字幕| 免费观看av网站的网址| 99香蕉大伊视频| 制服人妻中文乱码| 中文乱码字字幕精品一区二区三区| 精品亚洲乱码少妇综合久久| 不卡视频在线观看欧美| 2022亚洲国产成人精品| 有码 亚洲区| 老司机亚洲免费影院| 极品少妇高潮喷水抽搐| 精品酒店卫生间| 亚洲熟女精品中文字幕| 国产片内射在线| 亚洲欧美清纯卡通| 99久久人妻综合| 看十八女毛片水多多多| 亚洲精品,欧美精品| 亚洲内射少妇av| 免费在线观看视频国产中文字幕亚洲 | 一区二区三区四区激情视频| 在线免费观看不下载黄p国产| 国产xxxxx性猛交| 国产熟女欧美一区二区| 精品国产一区二区久久| av在线播放精品| 成人二区视频| 免费在线观看完整版高清| av福利片在线| 欧美日韩视频高清一区二区三区二| 国产精品免费视频内射| 男人添女人高潮全过程视频| 免费观看av网站的网址| 七月丁香在线播放| 欧美日韩国产mv在线观看视频| 天天躁夜夜躁狠狠躁躁| 成年人免费黄色播放视频| 亚洲熟女精品中文字幕| 美女xxoo啪啪120秒动态图| 欧美日韩精品网址| 国产精品久久久av美女十八| 国产高清不卡午夜福利| 免费播放大片免费观看视频在线观看| 黄色怎么调成土黄色| 女人被躁到高潮嗷嗷叫费观| 亚洲国产精品999| 亚洲人成电影观看| 欧美精品av麻豆av| 狂野欧美激情性bbbbbb| 亚洲三级黄色毛片| 亚洲一码二码三码区别大吗| 啦啦啦在线观看免费高清www| 一级毛片我不卡| 国产亚洲av片在线观看秒播厂| 欧美日韩综合久久久久久| 色吧在线观看| 高清欧美精品videossex| 国产亚洲av片在线观看秒播厂| 美女脱内裤让男人舔精品视频| 日本vs欧美在线观看视频| 女性被躁到高潮视频| 国产精品香港三级国产av潘金莲 | 免费在线观看视频国产中文字幕亚洲 | 久久人人爽人人片av| 在线精品无人区一区二区三| 日韩不卡一区二区三区视频在线| av国产久精品久网站免费入址| 久久精品国产亚洲av天美| 日本-黄色视频高清免费观看| 久久久久国产一级毛片高清牌| 下体分泌物呈黄色| 久久久精品免费免费高清| 2022亚洲国产成人精品| 久久久久久久久免费视频了| 国产成人欧美| 90打野战视频偷拍视频| 成人国产av品久久久| 久久精品久久久久久久性| 免费看不卡的av| 亚洲人成电影观看| 色视频在线一区二区三区| 中文字幕亚洲精品专区| 成人亚洲精品一区在线观看| 一级毛片 在线播放| 欧美日韩av久久| 赤兔流量卡办理| 只有这里有精品99| 久久国产亚洲av麻豆专区| 伊人久久大香线蕉亚洲五| 99国产综合亚洲精品| 久久久久国产一级毛片高清牌| 亚洲国产精品成人久久小说| 成人午夜精彩视频在线观看| 大香蕉久久成人网| 最近中文字幕高清免费大全6| 成人二区视频| 大香蕉久久网| 亚洲伊人久久精品综合| 一边摸一边做爽爽视频免费| 日韩熟女老妇一区二区性免费视频| 高清欧美精品videossex| 我的亚洲天堂| 欧美av亚洲av综合av国产av | 精品人妻一区二区三区麻豆| 日韩人妻精品一区2区三区| 久久午夜福利片| 国产淫语在线视频| 国产极品天堂在线| 人人妻人人添人人爽欧美一区卜| 亚洲伊人色综图| 亚洲国产精品一区二区三区在线| 搡老乐熟女国产| 亚洲国产欧美日韩在线播放| 欧美激情极品国产一区二区三区| 性少妇av在线| 老司机亚洲免费影院| 亚洲视频免费观看视频| 久久精品亚洲av国产电影网| 国精品久久久久久国模美| 涩涩av久久男人的天堂| 边亲边吃奶的免费视频| 天美传媒精品一区二区| av在线app专区| 国产极品天堂在线| 一二三四中文在线观看免费高清| av福利片在线| 色播在线永久视频| 精品国产一区二区三区四区第35| 这个男人来自地球电影免费观看 | 桃花免费在线播放| 国产精品久久久久久精品古装| 欧美精品一区二区大全| 99九九在线精品视频| 国产又色又爽无遮挡免| av片东京热男人的天堂| 麻豆精品久久久久久蜜桃| 国产有黄有色有爽视频| 青青草视频在线视频观看| 国产日韩欧美在线精品| 美女主播在线视频| 最近2019中文字幕mv第一页| 99精国产麻豆久久婷婷| 久久热在线av| 美女午夜性视频免费| 波多野结衣一区麻豆| av女优亚洲男人天堂| 国产成人精品在线电影| 青春草亚洲视频在线观看| 精品国产一区二区三区四区第35| 哪个播放器可以免费观看大片| 老司机影院成人| 久久精品国产亚洲av涩爱| 中文精品一卡2卡3卡4更新| 超碰97精品在线观看| 亚洲欧美一区二区三区黑人 | 80岁老熟妇乱子伦牲交| 99热国产这里只有精品6| 午夜激情av网站| 成人亚洲欧美一区二区av| 青春草亚洲视频在线观看| 精品亚洲成国产av| 街头女战士在线观看网站| 国产 精品1| 在线天堂最新版资源| 毛片一级片免费看久久久久| 欧美激情极品国产一区二区三区| 搡老乐熟女国产| 亚洲国产看品久久| 桃花免费在线播放| 狂野欧美激情性bbbbbb| 亚洲成色77777| av天堂久久9| tube8黄色片| 婷婷色综合大香蕉| 欧美日韩精品网址| 久久女婷五月综合色啪小说| 久久ye,这里只有精品| 精品视频人人做人人爽| 美女国产视频在线观看| 欧美成人精品欧美一级黄| 电影成人av| 美国免费a级毛片| 免费在线观看完整版高清| 亚洲成国产人片在线观看| 亚洲美女搞黄在线观看| 久久久亚洲精品成人影院| 性色avwww在线观看| 黄色配什么色好看| 国产亚洲午夜精品一区二区久久| 国产在线一区二区三区精| 久久99精品国语久久久| 色婷婷久久久亚洲欧美| 国产白丝娇喘喷水9色精品| 校园人妻丝袜中文字幕| 一级爰片在线观看| 婷婷色综合大香蕉| 黄色视频在线播放观看不卡| 亚洲av电影在线进入| 国产高清国产精品国产三级| 老女人水多毛片| 美女主播在线视频| 国产黄色视频一区二区在线观看| 人人妻人人澡人人爽人人夜夜| av一本久久久久| 国产亚洲av片在线观看秒播厂| 波多野结衣一区麻豆| 最近2019中文字幕mv第一页| 熟女电影av网| 黄色毛片三级朝国网站| 1024香蕉在线观看| 在线观看三级黄色| 久久精品国产a三级三级三级| 街头女战士在线观看网站| 成年人午夜在线观看视频| 人人妻人人澡人人爽人人夜夜| 国产一区二区 视频在线| 天天影视国产精品| 精品国产乱码久久久久久小说| 在线免费观看不下载黄p国产| 狠狠婷婷综合久久久久久88av| 亚洲,一卡二卡三卡| 1024香蕉在线观看| av网站在线播放免费| 精品人妻一区二区三区麻豆| 69精品国产乱码久久久| 亚洲精品成人av观看孕妇| 永久网站在线| 毛片一级片免费看久久久久| 国产精品二区激情视频| www.熟女人妻精品国产| 一级爰片在线观看| 极品人妻少妇av视频| 国产精品久久久久久久久免| 麻豆av在线久日| 三级国产精品片| 亚洲av国产av综合av卡| 久久国内精品自在自线图片| 边亲边吃奶的免费视频| 九色亚洲精品在线播放| 久久久久久久久久久久大奶| www日本在线高清视频| 日韩欧美一区视频在线观看| 成年动漫av网址| 欧美少妇被猛烈插入视频| 国产av一区二区精品久久| 青青草视频在线视频观看| 免费大片黄手机在线观看| 好男人视频免费观看在线| 欧美激情高清一区二区三区 | 少妇人妻久久综合中文| 欧美成人午夜免费资源| 99久久中文字幕三级久久日本| 国产免费又黄又爽又色| 久久久精品94久久精品| 久久亚洲国产成人精品v| 黑人欧美特级aaaaaa片| 国产一区二区激情短视频 | 日韩中文字幕欧美一区二区 | 女的被弄到高潮叫床怎么办| 国产熟女午夜一区二区三区| 熟妇人妻不卡中文字幕| 性少妇av在线| 国产色婷婷99| 欧美国产精品一级二级三级| 1024香蕉在线观看| 色视频在线一区二区三区| www.av在线官网国产| av在线播放精品| 人人妻人人澡人人爽人人夜夜| 国产在视频线精品| 成人毛片a级毛片在线播放| a 毛片基地| 欧美日韩成人在线一区二区| 亚洲国产欧美日韩在线播放| 日韩不卡一区二区三区视频在线| 国产乱人偷精品视频| 各种免费的搞黄视频| 成年女人在线观看亚洲视频| 成人影院久久| 黄色配什么色好看| 9191精品国产免费久久| 国产成人精品久久久久久| 欧美日韩av久久| 亚洲精品,欧美精品| 久久久久久久国产电影| 美女国产视频在线观看| 精品国产一区二区三区四区第35|