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

    溴化鑭探測(cè)器探測(cè)效率刻度及影響因素分析

    2021-10-26 01:52:48楊體波代光明廖光輝范新洋
    核技術(shù) 2021年10期
    關(guān)鍵詞:反射層溴化全能

    楊體波 王 敏 代光明 曹 彧 廖光輝 范新洋

    (成都理工大學(xué)核技術(shù)與自動(dòng)化工程學(xué)院 成都 610059)

    計(jì)算地層中各元素的含量、確定礦物成分是地質(zhì)礦產(chǎn)勘探中的一項(xiàng)重要工作。隨著測(cè)井技術(shù)的發(fā)展,基于核探測(cè)與分析的地層元素測(cè)井[1?4]技術(shù)在地層巖性識(shí)別中得到了較為廣泛的應(yīng)用。由于具有能量分辨率高、溫度性能好等諸多優(yōu)點(diǎn)[5?6],溴化鑭探測(cè)器(LaBr3:Ce晶體)在地層元素測(cè)井中取代了傳統(tǒng)的碘化鈉探測(cè)器。在地層元素測(cè)井中,所需要測(cè)量的俘獲γ 射線與非彈性γ 射線的能量范圍較大(100 keV~10 MeV),一般實(shí)驗(yàn)室所用的放射源發(fā)出的γ 射線的能量在2 MeV 以下,如果采用有源效率刻度,難以實(shí)現(xiàn)高能區(qū)的刻度[7]。無(wú)源效率刻度則是一種快捷、方便、準(zhǔn)確的方法,利用蒙特卡羅軟件,可以實(shí)現(xiàn)對(duì)溴化鑭探測(cè)器的無(wú)源效率刻度,確定出γ射線全能峰效率與能量之間的關(guān)系曲線。全能峰效率是評(píng)價(jià)探測(cè)器性能好壞的一項(xiàng)重要指標(biāo),影響全能峰效率的因素有很多,比如γ射線的能量、晶體尺寸、反射層厚度、外殼厚度等。鐘丁生等[8]使用MCNP(Monte Carlo N Particle Transport Code,基于蒙特卡羅方法的用于計(jì)算三維復(fù)雜幾何結(jié)構(gòu)中的中子、光子、電子或者耦合中子/光子/電子輸運(yùn)問(wèn)題的通用軟件包)軟件研究了溴化鑭探測(cè)器幾何因素對(duì)探測(cè)效率的影響,結(jié)果表明:增大晶體的幾何尺寸,探測(cè)效率增大,當(dāng)晶體尺寸增大到一定程度時(shí),探測(cè)效率值趨于穩(wěn)定;謝希成等[9]利用MCNP 軟件模擬了不同尺寸溴化鑭晶體的全能峰探測(cè)效率,與實(shí)驗(yàn)所得值進(jìn)行對(duì)比,其誤差在7%以內(nèi),并且通過(guò)擬合得到了點(diǎn)源效率函數(shù);鄧中華等[10]利用MCNP 軟件通過(guò)模擬獲得了溴化鑭探測(cè)器能量線性曲線、探測(cè)效率曲線,與實(shí)驗(yàn)值的相對(duì)誤差在5%以內(nèi);姚劍峰等[11]通過(guò)MCNP 模擬獲得了不同探測(cè)器厚度、晶體半徑以及源距對(duì)應(yīng)的全能峰效率,實(shí)現(xiàn)了對(duì)探測(cè)器的無(wú)源效率刻度;盧毅等[7]采用蒙特卡羅方法計(jì)算了不同尺寸溴化鑭閃爍體的吸收效率以及峰總比曲線,并 對(duì)?50 mm×10 mm 的LaBr3: Ce 探測(cè)器在0~10 MeV 的能段進(jìn)行了效率刻度;Favalli 等[5]研究了LaBr3: Ce 閃爍探測(cè)器在寬能量范圍(200~5 000 keV)下的全能峰效率;Yang 等[12]基于高純鍺探測(cè)器比較了三種效率刻度方法,即實(shí)驗(yàn)室有源效率刻度、實(shí)驗(yàn)室無(wú)源效率刻度軟件、半經(jīng)驗(yàn)公式與數(shù)值模擬相結(jié)合,證明了有源效率刻度是最準(zhǔn)確的方法,同時(shí),隨著高效計(jì)算技術(shù)的進(jìn)步,無(wú)源效率刻度、半經(jīng)驗(yàn)公式和數(shù)值模擬相結(jié)合的方法仍有望得到更廣泛的運(yùn)用。

    以上研究中,對(duì)探測(cè)器進(jìn)行無(wú)源效率刻度時(shí),只是考慮了晶體尺寸的大小、源距等參數(shù)對(duì)全能峰效率的影響,并沒有考慮反射層厚度以及探測(cè)器外殼厚度對(duì)全能峰效率的影響;在研究晶體尺寸大小對(duì)全能峰效率的影響時(shí),也沒有分析晶體長(zhǎng)度與晶體半徑對(duì)全能峰效率的具體影響。在此基礎(chǔ)上,本文利用MCNP 軟件研究了反射層、探測(cè)器外殼、晶體半徑、晶體長(zhǎng)度對(duì)全能峰效率的影響,并根據(jù)模擬結(jié)果選擇最優(yōu)參數(shù),針對(duì)地層元素測(cè)井中所使用的溴化鑭探測(cè)器(LaBr3: Ce 晶體尺寸為?3.18 cm×10.16 cm)進(jìn)行了無(wú)源效率刻度。

    1 理論基礎(chǔ)與模型建立

    1.1 理論基礎(chǔ)

    全能峰效率εp定義為γ譜全能峰下包含的計(jì)數(shù)與放射源在相同時(shí)間內(nèi)發(fā)射的γ 射線計(jì)數(shù)的比值[13],即:

    式中:Npeak為全能峰計(jì)數(shù),即能譜上全能峰的面積;Nγ為放射源發(fā)出的γ 射線總計(jì)數(shù)。實(shí)測(cè)時(shí),根據(jù)能譜找到所需要的全能峰和對(duì)應(yīng)的道址范圍,計(jì)算其面積,再減去本底計(jì)數(shù)即可得到全能峰計(jì)數(shù),也可以直接通過(guò)測(cè)量軟件讀出全能峰計(jì)數(shù)。利用MCNP 軟件進(jìn)行模擬時(shí),通過(guò)輸出文件讀出歸一化的概率即可得到全能峰效率。全能峰效率是描述探測(cè)器性能的一個(gè)重要指標(biāo),它主要與探測(cè)器的種類、γ射線的能量以及測(cè)量時(shí)的源距等因素有關(guān),利用式(2)[14]可以計(jì)算全能峰效率值:

    式中:G為幾何效率;I表示放射源發(fā)出的γ射線經(jīng)探測(cè)器Al殼、MgO反射層衰減后到達(dá)探測(cè)器晶體表面時(shí)的強(qiáng)度因子;εab為探測(cè)器晶體對(duì)γ 射線的吸收效率,對(duì)于不同能量的γ射線,其吸收效率不同;R為峰總比,即探測(cè)器的輸出能譜中,全能峰計(jì)數(shù)與總計(jì)數(shù)的比值,也可以通過(guò)MCNP軟件模擬得到。I可以通過(guò)式(3)[15]計(jì)算得出;

    式中:μ為線性衰減系數(shù),cm?1;x為衰減物質(zhì)的厚度,cm,在計(jì)算時(shí)μ、x分別對(duì)應(yīng)探測(cè)器Al殼和MgO反射層的線性衰減系數(shù)和厚度。吸收效率εab由式(4)[14]計(jì)算:

    此時(shí),μ、x使用LaBr3:Ce 晶體的線性衰減系數(shù)、厚度;幾何效率G由式(5)[14]計(jì)算得到:

    式中:r為探測(cè)晶體的半徑,cm;d為源距,cm,即源到探測(cè)器之間的距離。線性衰減系數(shù)μ與質(zhì)量衰減系數(shù)μm的關(guān)系可以由式(6)[16]得出:

    式中:ρ為衰減物質(zhì)的密度,g?cm?3;μm為質(zhì)量衰減系數(shù),cm2?g?1,可以通過(guò)WinXcom軟件獲得。

    1.2 模型建立

    參考實(shí)際的溴化鑭探測(cè)器尺寸,在MCNP 中構(gòu)建的探測(cè)器模型如圖1所示。以探測(cè)器Al殼左端面圓心為原點(diǎn)(0,0,0),點(diǎn)源位置設(shè)置為(?5,0,0),源距為5 cm,模擬能量為0.662 MeV的137Cs源,溴化鑭晶體直徑為3.18 cm,晶體長(zhǎng)度為10.16 cm,晶體外面包裹一層MgO反射層,其厚度為0.05 cm,最外層為Al殼,其端面厚度為0.1 cm,側(cè)面厚度為0.25 cm,右側(cè)端面為SiO2光導(dǎo),厚度為0.2 cm。計(jì)數(shù)卡采用F8脈沖計(jì)數(shù)卡,將能量劃分為1 024道來(lái)記錄γ光子的能量分布,模擬運(yùn)輸?shù)牧W訑?shù)為5×107。

    圖1 溴化鑭探測(cè)器模型Fig.1 Model for the lanthanum bromide detector in the MCNP calculation

    2 探測(cè)效率影響因素

    本文主要研究溴化鑭探測(cè)器的端面Al殼厚度、側(cè)面Al 殼厚度、MgO 反射層厚度、晶體長(zhǎng)度以及晶體半徑等參數(shù)對(duì)全能峰效率的影響。

    2.1 端面Al殼厚度對(duì)全能峰效率的影響

    研究端面Al 殼厚度對(duì)全能峰效率的影響時(shí),γ射線能量設(shè)置為0.662 MeV,源距為5 cm,MgO反射層厚度為0.05 cm,側(cè)面Al殼厚度為0.25 cm,只改變端 面Al 殼的 厚度(0.05 cm、0.10 cm、0.15 cm、0.20 cm、0.25 cm、0.30 cm、0.35 cm、0.40 cm),結(jié)果如圖2所示。

    圖2 端面Al殼厚度對(duì)全能峰效率的影響Fig.2 Effect of the thickness of Al shell in end-face of the detector on the full-energy peak efficiency

    由圖2可知,Al殼的厚度會(huì)影響全能峰效率;隨著Al殼厚度的增加全能峰效率降低,設(shè)計(jì)探測(cè)器尺寸時(shí),在考慮保護(hù)作用與屏蔽β射線的前提下,可以盡量降低Al殼的厚度,以提高全能峰探測(cè)效率。

    2.2 側(cè)面Al殼厚度對(duì)全能峰效率的影響

    研究側(cè)面Al 殼厚度對(duì)全能峰效率的影響時(shí),γ射線能量設(shè)置為0.662 MeV,源距為5 cm,MgO反射層厚度為0.05 cm,端面Al 殼厚度為0.1 cm,只改變側(cè) 面Al 殼的 厚度(0.05 cm、0.10 cm、0.15 cm、0.20 cm、0.25 cm、0.30 cm、0.35 cm、0.40 cm),結(jié)果如圖3所示。

    圖3 側(cè)面Al殼厚度對(duì)全能峰效率的影響Fig.3 Effect of the thickness of cylindrical Al shell of the detector on the full-energy peak efficiency

    由圖3可以看出,當(dāng)端面Al殼厚度一定時(shí),改變側(cè)面Al殼的厚度不會(huì)影響全能峰探測(cè)效率,在地層元素測(cè)井中使用溴化鑭探測(cè)器時(shí),可以適當(dāng)增加側(cè)面Al殼的厚度以提高對(duì)探測(cè)器的保護(hù)效率。

    2.3 MgO反射層對(duì)全能峰效率的影響

    模擬MgO 反射層對(duì)全能峰效率的影響時(shí),γ 射線能量設(shè)置為0.662 MeV,源距為5 cm,端面Al殼厚度為0.1 cm,側(cè)面Al 殼厚度為0.25 cm,分別模擬了MgO 反射層的端面和側(cè)面厚度同時(shí)為0.01 cm、0.03 cm、0.05 cm、0.07 cm、0.10 cm、0.15 cm、0.20 cm的情況,模擬結(jié)果如圖4所示。

    圖4 MgO反射層對(duì)全能峰效率的影響Fig.4 Effect of the thickness of the MgO reflective layer on the full-energy peak efficiency

    從圖4 中可以看到,增加MgO 反射層的厚度,全能峰效率隨之下降,其趨勢(shì)呈一條直線,為了提高探測(cè)器的探測(cè)效率,在保證起到反射目的的情況下,可以降低MgO反射層的厚度,一般取0.05 cm。

    2.4 晶體尺寸對(duì)全能峰效率的影響

    LaBr3:Ce晶體尺寸大小對(duì)全能峰效率的影響主要包括兩部分:晶體長(zhǎng)度的影響和晶體半徑的影響,模擬時(shí)γ 射線能量設(shè)置為0.662 MeV,源距為5 cm,端面Al殼厚度為0.1 cm,側(cè)面Al殼厚度為0.25 cm,MgO 反射層厚度為0.05 cm,分別模擬了晶體長(zhǎng)度為2.54 cm、5.08 cm、7.62 cm、10.16 cm,晶體半徑為1.59 cm、3.81 cm、4.77 cm、6.36 cm、7.95 cm、9.54 cm、11.13 cm的情況,模擬結(jié)果如圖5所示。

    從圖5 可以看出,在γ 射線的能量為0.662 MeV、源距為5 cm的情況下,隨著晶體半徑的增加,全能峰效率逐漸增加,但是增加的幅度越來(lái)越小。這是由于當(dāng)晶體半徑較小時(shí),點(diǎn)源相對(duì)探頭的立體角較小,幾何效率較底;隨著晶體半徑的增加,點(diǎn)源相對(duì)探頭的立體角增加,幾何效率逐漸增加,從而全能峰效率也隨之增加;當(dāng)半徑增加到某個(gè)數(shù)值后,再增加半徑,其立體角的變化很小,此時(shí)全能峰效率的變化也較小。同時(shí)也可以看到,當(dāng)晶體半徑一定時(shí),隨著晶體長(zhǎng)度的增加,全能峰效率也會(huì)增加。晶體長(zhǎng)度的增加會(huì)使得γ射線與晶體發(fā)生反應(yīng)的概率增加,更容易使其能量全部沉積于晶體中,進(jìn)而提高全能峰效率。由于增加晶體半徑比增加晶體長(zhǎng)度對(duì)探測(cè)效率的提升更加顯著,在地層元素測(cè)井中,如果地層孔徑的尺寸允許,應(yīng)盡量增加晶體半徑大小以提高探測(cè)效率。

    圖5 晶體尺寸對(duì)全能峰效率的影響Fig.5 Effect of the LaBr3:Ce crystal size on the full- energy peak efficiency

    為了進(jìn)一步探究晶體半徑以及晶體長(zhǎng)度對(duì)全能峰效率的影響,分別改變晶體的長(zhǎng)度與半徑(1.27 cm、2.54 cm、3.81 cm、5.08 cm、6.35 cm、7.62 cm)進(jìn)行模擬,晶體長(zhǎng)度與半徑的變化步長(zhǎng)設(shè)置為1.27 cm,對(duì)比晶體長(zhǎng)度和半徑變化對(duì)全能峰效率的影響,所得結(jié)果如圖6所示。

    圖6 對(duì)比晶體長(zhǎng)度變化和半徑變化對(duì)全能峰效率的影響Fig.6 Comparison of the effects of length and radius of the crystal on the full-energy peak efficiency

    由圖6 可知,晶體長(zhǎng)度變化和半徑變化對(duì)全能峰效率的影響并不相同,當(dāng)半徑較小時(shí)(<3.81 cm),晶體長(zhǎng)度的變化對(duì)全能峰效率的影響較??;當(dāng)晶體半徑較大時(shí)(3.81~7.62 cm),隨著晶體長(zhǎng)度的增加,全能峰效率也隨之增加,但增加的趨勢(shì)逐漸減緩。晶體半徑變化對(duì)全能峰效率的影響較大,隨著晶體半徑的增加全能峰效率也隨之增加;相比于晶體長(zhǎng)度較?。?.27 cm)的情況來(lái)說(shuō),當(dāng)晶體長(zhǎng)度較大時(shí)(>2.54 cm),晶體半徑變化對(duì)全能峰效率的影響更顯著。

    3 效率刻度

    對(duì)影響全能峰效率的因素進(jìn)行研究和分析后,對(duì)溴化鑭探測(cè)器進(jìn)行了無(wú)源效率刻度,獲得全能峰效率和能量之間的關(guān)系曲線。模擬時(shí),源距設(shè)置為5 cm,端面Al 殼厚度為0.1 cm,側(cè)面Al 殼為0.25 cm,MgO 反射層厚度為0.05 cm,LaBr3:Ce 晶體尺寸為?3.18 cm×10.16 cm。對(duì)地層元素測(cè)井中的非彈性散射γ射線、俘獲γ射線進(jìn)行分析后,發(fā)現(xiàn)γ射線的能量主要分布在100 keV~10 MeV之間,因此將模擬的γ 射線能量設(shè)置為:0.1 MeV、0.3 MeV、0.5 MeV、0.8 MeV、1.0 MeV、1.5 MeV、2.0 MeV、2.5 MeV、3.0 MeV、3.5 MeV、4.0 MeV、4.5 MeV、5.0 MeV、6.0 MeV、7.0 MeV、8.0 MeV、9.0 MeV、10.0 MeV,結(jié)果如圖7所示,通過(guò)擬合獲得了刻度曲線及其參數(shù),如式(7)所示:

    式中:y為全能峰效率;x為γ 射線的能量,MeV。其中,參數(shù)A1= 0.020 05,A2= 0.00394,t1= 0.382 65,t2= 2.78106,y0= 3.384 23 × 10?4。從圖7 中可以看出,隨著γ 射線能量的增加,全能峰效率呈下降趨勢(shì),這是由于能量越高的γ射線越容易穿透晶體,并沒有在晶體中沉積能量而被記錄下來(lái),全能峰效率隨之降低。

    圖7 溴化鑭探測(cè)器全能峰效率刻度曲線Fig.7 Curve of full-energy peak efficiency calibration of the LaBr3 detector

    為了確保模擬結(jié)果的準(zhǔn)確性,首先通過(guò)理論公式計(jì)算進(jìn)行驗(yàn)證。圖8 為WinXcom 軟件得到的LaBr3:Ce 晶體和Al 對(duì)不同γ 射線的質(zhì)量衰減系數(shù),由于MgO 的質(zhì)量衰減系數(shù)圖與Al的質(zhì)量衰減系數(shù)圖8(b)比較相似,避免重復(fù)而沒有進(jìn)行展示。通過(guò)式(6)可以求出LaBr3:Ce 晶體和Al 對(duì)不同γ 射線能量的線性衰減系數(shù)。

    圖8 LaBr3:Ce晶體(a)和Al (b)對(duì)不同能量γ射線的質(zhì)量吸收系數(shù)Fig.8 Mass absorption coefficient of the LaBr3:Ce crystal (a) and Al (b) vs. gamma ray energy

    通過(guò)MCNP模擬得到了不同能量的γ射線對(duì)應(yīng)的峰總比,并與文獻(xiàn)[14]中給出的標(biāo)準(zhǔn)數(shù)據(jù)進(jìn)行了對(duì)比,在參考標(biāo)準(zhǔn)數(shù)據(jù)時(shí),某些能量對(duì)應(yīng)的峰總比在文獻(xiàn)中并沒有給出具體數(shù)值,根據(jù)擬合的曲線通過(guò)插值法得到。對(duì)于某些晶體來(lái)說(shuō),峰總比的值可以通過(guò)模擬的方式得到,在保證模型準(zhǔn)確性的情況下,模擬數(shù)據(jù)可以直接用來(lái)計(jì)算全能峰效率。以上工作得到了LaBr3:Ce 晶體、Al、MgO 的質(zhì)量衰減系數(shù)和峰總比,結(jié)合式(2~6)可以計(jì)算全能峰效率值(表1)。

    通過(guò)表1 可以看出,峰總比的模擬值和理論值的誤差在6% 以內(nèi),全能峰效率的模擬值與理論值的誤差在7%以內(nèi),表明了模擬數(shù)據(jù)的可靠性。

    表1 峰總比和全能峰效率的模擬值與理論值Table 1 Simulated and theoretical values of peak to total ratio and full-energy peak efficiency of the LaBr3 detector

    為了進(jìn)一步驗(yàn)證數(shù)據(jù)的準(zhǔn)確性,利用LaBr3:Ce探測(cè)器分別對(duì)152Eu、60Co、137Cs 源進(jìn)行了測(cè)量,源距為5 cm,與模擬時(shí)的參數(shù)保持一致,根據(jù)測(cè)量的能譜計(jì)算了全能峰效率值,進(jìn)行效率刻度。三種放射源發(fā)出的γ 射線能量主要在121.78~1 408.93 keV 之間,將其能量代入模擬得到的刻度曲線中,計(jì)算出全能峰效率值,與利用實(shí)測(cè)譜求得的全能峰效率進(jìn)行對(duì)比,結(jié)果如表2所示。

    從表2可以看出,實(shí)驗(yàn)值與模擬值的誤差在6%以內(nèi),進(jìn)一步證明了模擬的準(zhǔn)確性。對(duì)于不具備有源效率刻度條件的實(shí)驗(yàn)室來(lái)說(shuō),無(wú)源效率刻度是一種可靠的方法。

    表2 全能峰效率實(shí)驗(yàn)值與模擬值Table 2 Simulated and measured values of full-energy peak efficiency of the LaBr3 detector

    4 結(jié)語(yǔ)

    運(yùn)用MCNP5軟件研究了溴化鑭探測(cè)器Al殼厚度、MgO 反射層厚度、晶體長(zhǎng)度和晶體半徑對(duì)全能峰效率的影響,針對(duì)地層元素測(cè)井中使用的溴化鑭探測(cè)器(LaBr3:Ce晶體,?3.18 cm×10.16 cm),獲得了100 keV~10 MeV能區(qū)的全能峰效率刻度曲線,并與理論計(jì)算結(jié)果和實(shí)測(cè)結(jié)果進(jìn)行了對(duì)比,得出以下結(jié)論:

    1)增加Al 殼厚度和MgO 反射層厚度會(huì)降低全能峰效率。

    2)增加晶體的長(zhǎng)度和半徑都會(huì)提高全能峰效率,且晶體半徑的變化對(duì)全能峰效率的影響較大;晶體長(zhǎng)度變化對(duì)全能峰效率的影響效果在不同晶體半徑時(shí)相差較大,當(dāng)晶體半徑較小時(shí),長(zhǎng)度變化對(duì)全能峰效率的影響較小,隨著半徑的增加,長(zhǎng)度變化對(duì)全能峰效率的影響逐漸增加,但是增加的幅度逐漸降低。

    3)為了提高探測(cè)效率,對(duì)于地層元素測(cè)井中使用的LaBr3:Ce探測(cè)器,在孔徑允許的情況下,應(yīng)盡量增加晶體的半徑、降低端面Al 殼和MgO 反射層的厚度;為了更好地保護(hù)探測(cè)器,可以適當(dāng)增加側(cè)面Al殼的厚度。

    4)全能峰效率刻度的模擬結(jié)果與理論計(jì)算結(jié)果和實(shí)測(cè)結(jié)果之間的誤差分別在7%和6%以內(nèi),證明了模擬結(jié)果的可靠性,同時(shí)也為實(shí)驗(yàn)室進(jìn)行無(wú)源效率刻度提供了一定的參考。

    猜你喜歡
    反射層溴化全能
    一種納米材料復(fù)合溴化丁基橡膠及其制備方法
    橡膠科技(2022年9期)2022-12-12 05:26:53
    全能小達(dá)人
    幼兒畫刊(2022年8期)2022-10-18 01:43:58
    全能小達(dá)人
    幼兒畫刊(2022年10期)2022-10-14 07:03:08
    全能小達(dá)人
    幼兒畫刊(2022年9期)2022-09-14 02:52:36
    全能小達(dá)人
    幼兒畫刊(2022年5期)2022-06-06 06:25:54
    反射層對(duì)NaI(Tl)閃爍體探測(cè)器探測(cè)效率的影響
    科技視界(2021年10期)2021-05-20 01:46:42
    基于有限元某壓水堆反射層組件模態(tài)分析的建模方法
    四川建筑(2020年3期)2020-07-18 01:56:24
    芳胺類化合物的原位氧化溴化工藝研究
    隨鉆雙感應(yīng)測(cè)井儀背景影響研究
    溴化丁基橡膠成套技術(shù)升級(jí)
    国产永久视频网站| 九九在线视频观看精品| 国产精品伦人一区二区| 少妇猛男粗大的猛烈进出视频 | 大陆偷拍与自拍| 高清在线视频一区二区三区| 大香蕉97超碰在线| a级毛片免费高清观看在线播放| 亚洲av国产av综合av卡| 国产成人午夜福利电影在线观看| 亚洲精品成人av观看孕妇| 校园人妻丝袜中文字幕| 美女国产视频在线观看| 午夜福利网站1000一区二区三区| 欧美潮喷喷水| 一级黄片播放器| 在线观看免费高清a一片| 亚洲自偷自拍三级| 久久99精品国语久久久| 亚洲色图av天堂| 亚洲在线观看片| 欧美xxxx性猛交bbbb| 搡老妇女老女人老熟妇| 久久99热这里只有精品18| 亚洲欧美一区二区三区黑人 | 能在线免费看毛片的网站| 国产欧美另类精品又又久久亚洲欧美| 日产精品乱码卡一卡2卡三| 国产精品.久久久| 国产精品麻豆人妻色哟哟久久 | 亚洲av成人av| 欧美日韩在线观看h| 国产伦在线观看视频一区| 亚洲av在线观看美女高潮| 男人爽女人下面视频在线观看| 天美传媒精品一区二区| 久久久久久久久久黄片| 国产爱豆传媒在线观看| 国产毛片a区久久久久| 大又大粗又爽又黄少妇毛片口| 亚洲成人一二三区av| 日本一二三区视频观看| 我要看日韩黄色一级片| 久久国产乱子免费精品| 在线天堂最新版资源| a级一级毛片免费在线观看| 极品教师在线视频| 国产精品熟女久久久久浪| 嫩草影院新地址| 欧美精品国产亚洲| freevideosex欧美| 欧美最新免费一区二区三区| 波多野结衣巨乳人妻| 看十八女毛片水多多多| 国产精品一区二区在线观看99 | 街头女战士在线观看网站| 美女大奶头视频| 国产精品人妻久久久影院| 一级毛片 在线播放| 黑人高潮一二区| www.色视频.com| 国产精品人妻久久久影院| 精品久久久久久久久av| 91aial.com中文字幕在线观看| 亚洲成人av在线免费| 97超碰精品成人国产| 别揉我奶头 嗯啊视频| 春色校园在线视频观看| 国产成人精品婷婷| 91av网一区二区| 永久网站在线| 亚洲,欧美,日韩| av在线老鸭窝| 麻豆精品久久久久久蜜桃| 大片免费播放器 马上看| 夫妻午夜视频| 国产精品av视频在线免费观看| 天堂中文最新版在线下载 | 女人久久www免费人成看片| 亚洲成人av在线免费| 日日摸夜夜添夜夜添av毛片| 晚上一个人看的免费电影| 97热精品久久久久久| 七月丁香在线播放| 身体一侧抽搐| 亚洲人成网站在线观看播放| 三级毛片av免费| 国产视频首页在线观看| 午夜激情久久久久久久| 国产黄色小视频在线观看| 好男人视频免费观看在线| 禁无遮挡网站| 国产黄a三级三级三级人| 七月丁香在线播放| 国产精品国产三级国产av玫瑰| or卡值多少钱| 少妇猛男粗大的猛烈进出视频 | 亚洲av电影在线观看一区二区三区 | 毛片一级片免费看久久久久| 亚洲婷婷狠狠爱综合网| 美女内射精品一级片tv| 成人国产麻豆网| 午夜福利高清视频| 亚洲最大成人中文| 国产成人精品一,二区| 国内精品宾馆在线| 舔av片在线| 午夜精品一区二区三区免费看| 欧美成人午夜免费资源| 日韩 亚洲 欧美在线| 超碰av人人做人人爽久久| 精品一区二区免费观看| 午夜福利高清视频| 成人高潮视频无遮挡免费网站| 日本三级黄在线观看| 在线 av 中文字幕| 视频中文字幕在线观看| 亚州av有码| 国产真实伦视频高清在线观看| 人体艺术视频欧美日本| 日韩,欧美,国产一区二区三区| 中国美白少妇内射xxxbb| 久久这里有精品视频免费| 十八禁国产超污无遮挡网站| 亚洲精品456在线播放app| 狂野欧美白嫩少妇大欣赏| 男插女下体视频免费在线播放| 男女视频在线观看网站免费| 又粗又硬又长又爽又黄的视频| 日韩大片免费观看网站| 日韩伦理黄色片| 亚洲精品日韩av片在线观看| 国产午夜精品久久久久久一区二区三区| 亚洲精品视频女| 国产高潮美女av| 午夜免费男女啪啪视频观看| 一级毛片aaaaaa免费看小| av在线老鸭窝| 精品久久久久久久久久久久久| 中文字幕制服av| 啦啦啦啦在线视频资源| 久久人人爽人人片av| 97精品久久久久久久久久精品| 成年人午夜在线观看视频 | 亚洲av不卡在线观看| 超碰97精品在线观看| 最近2019中文字幕mv第一页| 国产老妇伦熟女老妇高清| 一本久久精品| 亚洲婷婷狠狠爱综合网| 久久韩国三级中文字幕| 欧美zozozo另类| 免费观看性生交大片5| 偷拍熟女少妇极品色| 亚洲av免费在线观看| 国内揄拍国产精品人妻在线| 久久99精品国语久久久| 天堂√8在线中文| 国产视频首页在线观看| 美女内射精品一级片tv| 国产亚洲一区二区精品| 秋霞伦理黄片| 欧美激情国产日韩精品一区| 亚洲精品aⅴ在线观看| 国产精品综合久久久久久久免费| 日韩不卡一区二区三区视频在线| 国产男人的电影天堂91| 国产成人a∨麻豆精品| 午夜福利视频1000在线观看| 国产片特级美女逼逼视频| 成人鲁丝片一二三区免费| 欧美高清性xxxxhd video| 国产淫语在线视频| 久久久久久久久久成人| 国产成人freesex在线| 国产精品综合久久久久久久免费| 国产亚洲5aaaaa淫片| 日韩不卡一区二区三区视频在线| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 91精品伊人久久大香线蕉| 一区二区三区免费毛片| 亚洲欧美清纯卡通| 国产美女午夜福利| 秋霞在线观看毛片| 麻豆精品久久久久久蜜桃| 干丝袜人妻中文字幕| 午夜福利成人在线免费观看| 男人狂女人下面高潮的视频| 女人被狂操c到高潮| 精品99又大又爽又粗少妇毛片| 狠狠精品人妻久久久久久综合| 久久久久久久大尺度免费视频| 亚洲精品国产av蜜桃| 亚洲精品乱码久久久v下载方式| 亚洲人成网站在线播| 极品少妇高潮喷水抽搐| 成人综合一区亚洲| 国产在线男女| 伦精品一区二区三区| 成人亚洲精品一区在线观看 | 久久久a久久爽久久v久久| 最近手机中文字幕大全| 高清欧美精品videossex| 熟妇人妻不卡中文字幕| 国产精品麻豆人妻色哟哟久久 | 人妻夜夜爽99麻豆av| 久久热精品热| 蜜臀久久99精品久久宅男| 国产女主播在线喷水免费视频网站 | av福利片在线观看| 成人亚洲精品一区在线观看 | 亚洲精华国产精华液的使用体验| 少妇的逼水好多| 国产毛片a区久久久久| 一区二区三区高清视频在线| 看非洲黑人一级黄片| 国产高清国产精品国产三级 | 日韩成人伦理影院| 久久久久久久久久成人| 99视频精品全部免费 在线| 国产午夜福利久久久久久| 欧美日韩在线观看h| 我的女老师完整版在线观看| 91久久精品国产一区二区三区| 亚洲aⅴ乱码一区二区在线播放| 毛片女人毛片| 久久99热6这里只有精品| 精品久久久久久久久久久久久| 联通29元200g的流量卡| 久久久久久久久久成人| 亚洲电影在线观看av| 欧美bdsm另类| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 久久久久久久久中文| 少妇熟女aⅴ在线视频| 熟女电影av网| 亚洲精品,欧美精品| 爱豆传媒免费全集在线观看| 青青草视频在线视频观看| 简卡轻食公司| av在线老鸭窝| 亚洲,欧美,日韩| 永久网站在线| 国产 一区精品| 亚洲av免费在线观看| 午夜亚洲福利在线播放| 高清欧美精品videossex| 少妇熟女欧美另类| 欧美人与善性xxx| 欧美日韩精品成人综合77777| 日本wwww免费看| 91久久精品国产一区二区三区| av卡一久久| 狠狠精品人妻久久久久久综合| 午夜免费观看性视频| 国产精品.久久久| 午夜日本视频在线| 丰满人妻一区二区三区视频av| 校园人妻丝袜中文字幕| 麻豆成人午夜福利视频| 如何舔出高潮| 欧美激情国产日韩精品一区| 欧美三级亚洲精品| 亚洲熟女精品中文字幕| 国产黄a三级三级三级人| 少妇猛男粗大的猛烈进出视频 | 嫩草影院入口| 免费看日本二区| 久久精品国产亚洲网站| 一本久久精品| 超碰97精品在线观看| 免费人成在线观看视频色| 美女脱内裤让男人舔精品视频| 国产精品福利在线免费观看| 男女边吃奶边做爰视频| 亚洲高清免费不卡视频| 成人鲁丝片一二三区免费| 日韩欧美一区视频在线观看 | 国内精品一区二区在线观看| 日日啪夜夜撸| 国产一区亚洲一区在线观看| 欧美zozozo另类| 一级av片app| 国产午夜精品一二区理论片| 边亲边吃奶的免费视频| 亚洲综合色惰| 免费观看性生交大片5| 99热全是精品| 亚洲人成网站高清观看| 我的老师免费观看完整版| 久久精品夜夜夜夜夜久久蜜豆| 国产欧美另类精品又又久久亚洲欧美| 亚洲最大成人手机在线| 免费看日本二区| 国产伦理片在线播放av一区| 婷婷六月久久综合丁香| 精品欧美国产一区二区三| 欧美成人午夜免费资源| 亚洲在久久综合| 蜜臀久久99精品久久宅男| 国产成人91sexporn| 国产黄色小视频在线观看| 成年免费大片在线观看| av卡一久久| 成人性生交大片免费视频hd| 中文字幕av在线有码专区| 网址你懂的国产日韩在线| av播播在线观看一区| 久久久久久久久大av| 一夜夜www| 日韩一区二区三区影片| 日本一本二区三区精品| 岛国毛片在线播放| 少妇人妻一区二区三区视频| 亚洲熟妇中文字幕五十中出| 精品久久久精品久久久| 欧美 日韩 精品 国产| 久久久久久久久久成人| 亚洲av国产av综合av卡| 有码 亚洲区| 亚洲精品国产av成人精品| 亚洲人成网站高清观看| 亚洲成色77777| 26uuu在线亚洲综合色| 国产综合精华液| 麻豆乱淫一区二区| 亚洲欧美成人精品一区二区| 成人毛片a级毛片在线播放| 国产亚洲最大av| 99久国产av精品| 美女脱内裤让男人舔精品视频| 美女大奶头视频| 一区二区三区乱码不卡18| 日本三级黄在线观看| h日本视频在线播放| 国产不卡一卡二| 天美传媒精品一区二区| 成人二区视频| 日韩在线高清观看一区二区三区| 国产黄色免费在线视频| 激情 狠狠 欧美| 最近的中文字幕免费完整| 亚洲av中文字字幕乱码综合| 纵有疾风起免费观看全集完整版 | 国产一区二区亚洲精品在线观看| 99久国产av精品国产电影| 亚洲精品一二三| 国产黄色免费在线视频| 听说在线观看完整版免费高清| 亚洲精品第二区| 国产伦在线观看视频一区| 国产精品无大码| 美女主播在线视频| 欧美性猛交╳xxx乱大交人| 一级黄片播放器| 麻豆成人av视频| 国产精品无大码| 日本av手机在线免费观看| 久久久久性生活片| 亚洲av中文字字幕乱码综合| 国产永久视频网站| 国产亚洲午夜精品一区二区久久 | 一区二区三区四区激情视频| 久久草成人影院| h日本视频在线播放| 少妇丰满av| 五月天丁香电影| 午夜免费观看性视频| 永久免费av网站大全| 嫩草影院精品99| 韩国高清视频一区二区三区| 国产精品不卡视频一区二区| 亚洲美女视频黄频| 一二三四中文在线观看免费高清| 精品99又大又爽又粗少妇毛片| av天堂中文字幕网| 亚洲综合色惰| 美女cb高潮喷水在线观看| 一级爰片在线观看| 久久久久久九九精品二区国产| 欧美另类一区| 亚洲精品成人av观看孕妇| av在线亚洲专区| 国产不卡一卡二| 成人无遮挡网站| 成人综合一区亚洲| 精品人妻偷拍中文字幕| 综合色av麻豆| 亚洲av在线观看美女高潮| 国产熟女欧美一区二区| 五月玫瑰六月丁香| 91精品国产九色| 久久久久久久久久久免费av| 国产精品一及| 午夜激情欧美在线| 欧美+日韩+精品| 亚洲欧洲国产日韩| 欧美成人一区二区免费高清观看| 日本一本二区三区精品| 国内少妇人妻偷人精品xxx网站| 2018国产大陆天天弄谢| 午夜激情福利司机影院| 日韩一区二区视频免费看| 国产av在哪里看| 国产av国产精品国产| av专区在线播放| 国产精品国产三级国产av玫瑰| 五月天丁香电影| 在线播放无遮挡| 精品一区二区免费观看| 午夜爱爱视频在线播放| 国产伦理片在线播放av一区| 欧美日韩精品成人综合77777| 日本一二三区视频观看| 一级毛片aaaaaa免费看小| 免费高清在线观看视频在线观看| 免费观看无遮挡的男女| 午夜久久久久精精品| 日韩欧美 国产精品| 欧美另类一区| 天堂俺去俺来也www色官网 | 身体一侧抽搐| 在线观看美女被高潮喷水网站| 又爽又黄a免费视频| 国产成人精品久久久久久| 国产亚洲最大av| 九九爱精品视频在线观看| 国内精品美女久久久久久| 舔av片在线| 色尼玛亚洲综合影院| 我的老师免费观看完整版| 中文欧美无线码| av一本久久久久| 天堂影院成人在线观看| 欧美三级亚洲精品| 国产精品女同一区二区软件| 尤物成人国产欧美一区二区三区| 欧美日韩综合久久久久久| 亚洲欧洲日产国产| 建设人人有责人人尽责人人享有的 | 午夜老司机福利剧场| 成人美女网站在线观看视频| 亚洲欧美日韩无卡精品| 久久综合国产亚洲精品| 亚洲天堂国产精品一区在线| 男女边摸边吃奶| 97热精品久久久久久| 久久久精品免费免费高清| 男女下面进入的视频免费午夜| 久久久久精品性色| 欧美激情久久久久久爽电影| 国产欧美另类精品又又久久亚洲欧美| 亚洲av不卡在线观看| 国产免费福利视频在线观看| 美女xxoo啪啪120秒动态图| 国产高清三级在线| 免费高清在线观看视频在线观看| 国产精品久久久久久av不卡| 久久久久久久久久成人| 777米奇影视久久| 麻豆av噜噜一区二区三区| 亚洲第一区二区三区不卡| 国产精品蜜桃在线观看| 网址你懂的国产日韩在线| 国产精品日韩av在线免费观看| 午夜福利成人在线免费观看| 日日啪夜夜撸| 美女cb高潮喷水在线观看| 久久亚洲国产成人精品v| 高清日韩中文字幕在线| 亚洲成人一二三区av| 亚洲欧美精品自产自拍| 校园人妻丝袜中文字幕| 成人一区二区视频在线观看| eeuss影院久久| 亚洲国产精品成人久久小说| 一个人看的www免费观看视频| 亚洲精品色激情综合| 别揉我奶头 嗯啊视频| 日韩欧美 国产精品| 国产精品久久视频播放| 日韩,欧美,国产一区二区三区| xxx大片免费视频| 内地一区二区视频在线| 男的添女的下面高潮视频| 欧美zozozo另类| 日韩成人伦理影院| 我的女老师完整版在线观看| 日韩一区二区三区影片| 人妻夜夜爽99麻豆av| 91久久精品国产一区二区成人| 亚洲av在线观看美女高潮| 黄片无遮挡物在线观看| 麻豆成人午夜福利视频| 免费无遮挡裸体视频| 狂野欧美激情性xxxx在线观看| 99热这里只有是精品50| 国产黄色免费在线视频| 亚洲va在线va天堂va国产| 熟妇人妻久久中文字幕3abv| 少妇的逼好多水| 久久久久久久久久久丰满| 国产三级在线视频| 久久精品国产亚洲网站| 国产免费福利视频在线观看| 少妇猛男粗大的猛烈进出视频 | 亚洲国产成人一精品久久久| 中文字幕av成人在线电影| 99久国产av精品国产电影| 欧美激情久久久久久爽电影| 国国产精品蜜臀av免费| 日韩视频在线欧美| 在线免费观看的www视频| 女的被弄到高潮叫床怎么办| 91午夜精品亚洲一区二区三区| av女优亚洲男人天堂| 22中文网久久字幕| 亚洲精品第二区| 中文天堂在线官网| 99久国产av精品国产电影| 看十八女毛片水多多多| 超碰av人人做人人爽久久| 精品一区二区三卡| 精品一区二区免费观看| 国产精品伦人一区二区| 午夜福利在线在线| 精品一区二区三区视频在线| 成人亚洲欧美一区二区av| 一个人观看的视频www高清免费观看| 91精品一卡2卡3卡4卡| 国产亚洲91精品色在线| 精品人妻一区二区三区麻豆| 精品国产三级普通话版| 又爽又黄无遮挡网站| 国产午夜精品久久久久久一区二区三区| 麻豆成人午夜福利视频| 最近中文字幕2019免费版| 69人妻影院| 人妻一区二区av| 国产精品爽爽va在线观看网站| 夜夜爽夜夜爽视频| 高清日韩中文字幕在线| 高清在线视频一区二区三区| 国产伦精品一区二区三区视频9| 蜜桃亚洲精品一区二区三区| 午夜福利在线在线| 少妇人妻精品综合一区二区| av在线蜜桃| 亚洲va在线va天堂va国产| 在线a可以看的网站| 日本免费a在线| 欧美高清成人免费视频www| 亚洲内射少妇av| 国产综合懂色| av免费在线看不卡| 国产久久久一区二区三区| 26uuu在线亚洲综合色| 国产精品伦人一区二区| 人人妻人人澡欧美一区二区| 亚洲成人久久爱视频| 两个人的视频大全免费| 亚洲av免费高清在线观看| 欧美区成人在线视频| 午夜福利视频精品| 免费av不卡在线播放| 一二三四中文在线观看免费高清| 黄色日韩在线| 性插视频无遮挡在线免费观看| 欧美不卡视频在线免费观看| 最近中文字幕高清免费大全6| 欧美97在线视频| 九色成人免费人妻av| 国产高清三级在线| 国产综合精华液| 久久亚洲国产成人精品v| 午夜亚洲福利在线播放| 91精品一卡2卡3卡4卡| 两个人视频免费观看高清| .国产精品久久| 久久久精品欧美日韩精品| 青春草视频在线免费观看| 午夜激情久久久久久久| 亚洲丝袜综合中文字幕| 乱人视频在线观看| 亚洲经典国产精华液单| 中文乱码字字幕精品一区二区三区 | 欧美日韩在线观看h| 久久久久久久久久久丰满| 国产男人的电影天堂91| 欧美日本视频| 色播亚洲综合网| 亚洲电影在线观看av| 能在线免费观看的黄片| 久久国产乱子免费精品| 欧美 日韩 精品 国产| 亚洲欧洲国产日韩| 成人美女网站在线观看视频| 国产精品三级大全| 高清视频免费观看一区二区 | 色综合色国产| 成人无遮挡网站| 日韩欧美 国产精品| a级毛色黄片| 免费播放大片免费观看视频在线观看| 午夜久久久久精精品| 91久久精品国产一区二区三区| 不卡视频在线观看欧美| 欧美xxxx黑人xx丫x性爽| 午夜福利在线观看免费完整高清在| 亚洲综合精品二区| 中文天堂在线官网| 99热这里只有是精品在线观看| 看非洲黑人一级黄片| 久久精品国产鲁丝片午夜精品| 亚洲熟妇中文字幕五十中出|