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

    基于天然μ子散射特征評(píng)估鈾板內(nèi)水平狹縫寬度的方法

    2017-05-22 09:58:25肖灑何偉波帥茂兵蘭明聰韋孟伏
    核技術(shù) 2017年5期
    關(guān)鍵詞:通量分類(lèi)器寬度

    肖灑 何偉波 帥茂兵 蘭明聰 韋孟伏

    ?

    基于天然μ子散射特征評(píng)估鈾板內(nèi)水平狹縫寬度的方法

    肖灑 何偉波 帥茂兵 蘭明聰 韋孟伏

    (中國(guó)工程物理研究院材料研究所 江油 621907)

    針對(duì)經(jīng)典μ子成像方法用于特殊對(duì)象結(jié)構(gòu)探測(cè)的技術(shù)難點(diǎn),提出了一種通過(guò)直接比較對(duì)象發(fā)生結(jié)構(gòu)變化前后μ子散射特征的差異來(lái)進(jìn)行形變類(lèi)型和尺度判定的新思路,稱(chēng)為正向參比法。相比于經(jīng)典的逆向反演法,該方法在原理上對(duì)于初始結(jié)構(gòu)已知的特殊核對(duì)象的結(jié)構(gòu)變化探測(cè)更有優(yōu)勢(shì)。從初步模擬結(jié)果來(lái)看,在合理的天然μ子通量下,該方法可實(shí)現(xiàn)鈾板內(nèi)亞毫米尺度水平狹縫的存在及其尺度的準(zhǔn)確判定。

    μ子,重核材料,多重庫(kù)侖散射,非參數(shù)檢驗(yàn),接受者操作特征分析

    μ子斷層成像技術(shù)憑借其無(wú)需額外放射源、安全系數(shù)高、穿透能力強(qiáng)、對(duì)重核材料靈敏等眾多優(yōu)勢(shì)成為了輻射成像領(lǐng)域一個(gè)新興的研究熱點(diǎn)[1?3]。在利用輻射成像手段進(jìn)行目標(biāo)對(duì)象結(jié)構(gòu)診斷相關(guān)領(lǐng)域中,常需要研究重核材料(例如鈾)及其元件經(jīng)過(guò)一段時(shí)間放置或特定工藝處理后的內(nèi)部結(jié)構(gòu)變化問(wèn)題,并且目標(biāo)對(duì)象的初始結(jié)構(gòu)組成通常是可以準(zhǔn)確預(yù)知的。對(duì)于此類(lèi)問(wèn)題,傳統(tǒng)的μ子散射成像反演法的基本思路是:分別對(duì)變化前、后兩個(gè)目標(biāo)對(duì)象進(jìn)行μ子散射特征測(cè)量并利用特定算法反演出其內(nèi)部三維結(jié)構(gòu)圖像,通過(guò)對(duì)比二者圖像的差異即可獲知目標(biāo)對(duì)象之間的結(jié)構(gòu)差異。這是一種直接、自然的方法,對(duì)于某些類(lèi)型的目標(biāo)對(duì)象也是行之有效的,但該方法往往涉及不確定性大、復(fù)雜性高的圖像反演過(guò)程。如果目標(biāo)對(duì)象組成結(jié)構(gòu)較為復(fù)雜、尤其是具有大量孔洞、狹縫或低密度內(nèi)填充物等各種類(lèi)空心結(jié)構(gòu)時(shí),傳統(tǒng)的成像反演方法極有可能受到嚴(yán)重的偽影干擾而難以給出清晰、準(zhǔn)確的圖像,因而使其應(yīng)用范圍大為受限?;谏鲜鲈?,亟需發(fā)展一種適用范圍更廣、實(shí)現(xiàn)難度更低的基于μ子散射特征來(lái)識(shí)別目標(biāo)對(duì)象之間結(jié)構(gòu)差異的新方法。

    在這種背景下,我們提出了一種基于目標(biāo)對(duì)象發(fā)生結(jié)構(gòu)變化前后對(duì)天然宇宙射線μ子的散射特征變化來(lái)識(shí)別目標(biāo)對(duì)象所發(fā)生的結(jié)構(gòu)變化類(lèi)型以及評(píng)估該特定結(jié)構(gòu)變化的尺度的新方法。由于該方法不涉及逆向的圖像反演(重建)過(guò)程,而是采用了與預(yù)先建立的參比數(shù)據(jù)庫(kù)進(jìn)行比對(duì)的方法。并且,由于該方法規(guī)避了繁瑣且易被干擾的圖像重建過(guò)程,與目標(biāo)對(duì)象的初始結(jié)構(gòu)復(fù)雜程度無(wú)關(guān),因而對(duì)于重核材料對(duì)象結(jié)構(gòu)變化診斷相關(guān)領(lǐng)域而言,其適用性要比傳統(tǒng)成像反演方法更為廣泛。

    1 方法原理

    為了解決傳統(tǒng)的μ子散射成像反演法用于特殊目標(biāo)對(duì)象結(jié)構(gòu)差異診斷上的難題,我們嘗試了一種新的方法,如圖1所示。其基本思想是:首先,建立目標(biāo)對(duì)象標(biāo)準(zhǔn)狀態(tài)下的理論模型(P1),然后充分考慮目標(biāo)對(duì)象所有可能發(fā)生的一系列結(jié)構(gòu)變化(包括形變類(lèi)型、尺度、位置等多個(gè)參量),建立一組或多組與之分別對(duì)應(yīng)的微擾模型(P2)。在這一步驟中,微擾模型體系必須具有很好的完備性。其次,通過(guò)蒙特卡羅模擬獲得上述各個(gè)微擾模型對(duì)應(yīng)的μ子散射特征量(P3),然后將這些散射特征量與相應(yīng)的結(jié)構(gòu)變化形成一一對(duì)應(yīng)的映射關(guān)系,從而組成一套完備的參比數(shù)據(jù)庫(kù)待用(P4)。P1?P4這幾個(gè)步驟(圖1中虛線框內(nèi))構(gòu)成了該方法的先決條件。而在日常應(yīng)用(步驟1?3)中,對(duì)于實(shí)際待測(cè)目標(biāo)對(duì)象,采用μ子成像平臺(tái)實(shí)驗(yàn)測(cè)量方法獲得其μ子散射特征量,然后將該散射特征與數(shù)據(jù)庫(kù)中的數(shù)據(jù)進(jìn)行比對(duì),尋找出相似度最高的條目,即可在給定的顯著性水平下,判定待測(cè)目標(biāo)對(duì)象發(fā)生了相應(yīng)的結(jié)構(gòu)變化。由于該方法摒棄了較為繁瑣、易引入誤差的圖像反演過(guò)程,試圖通過(guò)直接比較兩種狀態(tài)下的目標(biāo)對(duì)象對(duì)μ子的散射特征的差異來(lái)獲取目標(biāo)對(duì)象的結(jié)構(gòu)變化信息,因而僅涉及正向的散射測(cè)量過(guò)程,故稱(chēng)之為“正向參比法”。

    圖1 正向參比法的實(shí)現(xiàn)及其所涉及的關(guān)鍵技術(shù)

    2 正向參比法的實(shí)現(xiàn)

    以水平放置的鈾板的GEANT4模擬結(jié)果為例,對(duì)正向參比法用于鈾板內(nèi)部水平狹縫寬度識(shí)別的實(shí)現(xiàn)以及其中的關(guān)鍵技術(shù)作詳細(xì)說(shuō)明。

    2.1 μ子散射特征量參比數(shù)據(jù)庫(kù)的建立

    2.1.1 微擾模型的建立

    本文采用的鈾板狹縫模型如圖2所示,左、右兩邊分別為理想鈾板模型和基于該模型擴(kuò)展得到的一組分別帶有寬度為d(1≤≤,為指定的有限正整數(shù))的狹縫的鈾板模型。為簡(jiǎn)單起見(jiàn),假定鈾板只可能產(chǎn)生水平方向的狹縫,并且所有可能的狹縫寬度的集合={d}為有限集,至少是可以近似看作有限集。由于微擾模型必須具有完備性,因此需要針對(duì)集合中的每一個(gè)元素d分別設(shè)計(jì)相應(yīng)的狹縫模型,然后進(jìn)行后續(xù)的μ子散射蒙特卡羅模擬。

    圖2 鈾板狹縫模型

    2.1.2 μ子散射的蒙特卡羅模擬

    為了獲得足夠多的μ子散射徑跡,我們編制了一套基于GEANT4的μ子成像系統(tǒng)來(lái)模擬μ子在目標(biāo)對(duì)象中的多重庫(kù)侖散射(Multiple Coulomb Scattering, MCS)行為。天然宇宙射線μ子源項(xiàng)采用CRY軟件包[4]來(lái)生成。模擬計(jì)算中考慮了μ子與材料的庫(kù)侖散射、電離以及軔致輻射等主要物理作用。

    2.1.3 數(shù)據(jù)庫(kù)的構(gòu)建

    根據(jù)Moliere經(jīng)驗(yàn)公式,μ子散射極角概率分布函數(shù)()對(duì)不同目標(biāo)對(duì)象具有較好的特異性,故本文選取()作為μ子散射特征量。然而,由于天然μ子能譜范圍極廣、入射方向各異,()往往難以采用精確的解析函數(shù)進(jìn)行描述。并且由于μ子與物質(zhì)多重庫(kù)侖散射的隨機(jī)性,即使對(duì)于同一目標(biāo),兩次實(shí)驗(yàn)得到的散射角分布也可能存在差異??紤]這兩方面的因素,實(shí)際應(yīng)用中,需要進(jìn)行同一條件下的多次重復(fù)實(shí)驗(yàn)來(lái)構(gòu)建總體()的多個(gè)獨(dú)立樣本。

    2.2 目標(biāo)對(duì)象結(jié)構(gòu)變化的判定——基于假設(shè)檢驗(yàn)的二元分類(lèi)器

    根據(jù)上述分析,利用正向參比法評(píng)估目標(biāo)物體結(jié)構(gòu)變化的核心問(wèn)題可歸結(jié)為待測(cè)目標(biāo)對(duì)象μ子散射特征量與預(yù)建數(shù)據(jù)庫(kù)中μ子散射特征量之間的比對(duì)問(wèn)題。考慮到實(shí)際采用的是多次重復(fù)實(shí)驗(yàn)獲得的多個(gè)相互獨(dú)立的散射角分布樣本,因此問(wèn)題的實(shí)質(zhì)是:在總體分布形式()未知的條件下,檢驗(yàn)兩個(gè)獨(dú)立樣本是否取自同一總體的非參數(shù)檢驗(yàn)問(wèn)題。這里采用比較累積頻數(shù)分布函數(shù)的Kolmogorov- Smirnov (KS)檢驗(yàn)方法[5],其具體過(guò)程如下:

    令為單個(gè)μ子的散射角度量:

    = (Δθ)2+ (Δθ)2(1)

    式中:Δθ、Δθ分別為和平面上的累積偏轉(zhuǎn)角。假定兩次獨(dú)立實(shí)驗(yàn)的目標(biāo)物體對(duì)應(yīng)的天然μ子散射角總體分別服從1()和2()分布,則兩次實(shí)驗(yàn)測(cè)得的散射角可以看作分別來(lái)自這兩個(gè)總體的樣本??紤]假設(shè)檢驗(yàn)問(wèn)題:

    在給定的顯著性水平下,計(jì)算KS檢驗(yàn)的值(零假設(shè)0可被拒絕的最小顯著性水平)。如果≤,則拒絕0,認(rèn)為兩次實(shí)驗(yàn)?zāi)繕?biāo)對(duì)象存在差異;反之,如果>,則接受0,認(rèn)為兩次實(shí)驗(yàn)?zāi)繕?biāo)對(duì)象完全相同。

    基于上述KS檢驗(yàn)對(duì)0拒絕/接受的結(jié)論,可以構(gòu)建一個(gè)判定目標(biāo)對(duì)象是否發(fā)生變化的二元分類(lèi)器。令顯著性水平在一定范圍內(nèi)變化,即可在接受者操作特征(Receiver Operating Characteristic, ROC)空間內(nèi)繪制出相應(yīng)的ROC曲線。ROC曲線可以對(duì)該分類(lèi)器的性能進(jìn)行可視化,給出分類(lèi)器的真陽(yáng)性率(True Positive Rate, TPR)(收益)與假陽(yáng)性率(False Positive Rate, FPR)(開(kāi)銷(xiāo))之間的相對(duì)平衡關(guān)系[6?7]。下文中,我們會(huì)通過(guò)ROC曲線來(lái)評(píng)估該基于KS檢驗(yàn)的二元分類(lèi)器性能。

    3 結(jié)果與討論

    3.1 不同μ子通量下可識(shí)別的最小狹縫寬度

    由于所用μ子通量不同會(huì)導(dǎo)致散射樣本容量不同,進(jìn)而可能影響假設(shè)檢驗(yàn)的判定結(jié)果,因此有必要對(duì)上述方法在各種不同μ子通量條件下對(duì)目標(biāo)對(duì)象是否存在差異的識(shí)別能力進(jìn)行研究。

    這里以厚度1 cm的鈾質(zhì)均勻平板作為研究對(duì)象,在平板中心水平面上分別設(shè)置寬度為1= 0.1mm、2=0.2 mm、3=0.5 mm和4=1.0 mm的4種水平狹縫。以40組理想鈾板(狹縫寬度0=0)的μ子散射模擬結(jié)果獲取的40個(gè)散射角分布樣本,以及40組帶有寬度d(為[1, 4]范圍內(nèi)某一指定正整數(shù))狹縫的鈾板的μ子散射模擬結(jié)果獲取的40個(gè)散射角分布樣本共同組成參比數(shù)據(jù)庫(kù),以帶有寬度d狹縫的鈾板作為待測(cè)目標(biāo)對(duì)象,依次選取μ子通量為50 cm?2、125 cm?2、250 cm?2和1000 cm?2,分別采用上述正向參比法考察在一定的顯著性水平下,各種μ子通量條件下該方法對(duì)狹縫的識(shí)別能力。對(duì)于給定的,一對(duì)(TPR, FPR)就對(duì)應(yīng)了ROC空間的一個(gè)點(diǎn)。在(0, 1)范圍內(nèi)改變,就可以得到一系列ROC空間點(diǎn)組成的一條ROC曲線,如圖3?5所示。根據(jù)ROC曲線的含義可知,曲線越靠近左上角,即曲線下面積(Area Under Curve, AUC)越接近1,相應(yīng)的分類(lèi)器性能也越好。由圖3可知,50cm?2的μ子通量下該分類(lèi)器已可實(shí)現(xiàn)對(duì)寬度低至1.0 mm的水平狹縫的完美分辨(AUC約為1),相應(yīng)的閾值≈ 0。

    圖3 μ子通量為50 cm?2時(shí),分類(lèi)器對(duì)于帶有寬度di的水平狹縫的鈾板的識(shí)別性能ROC曲線 (a) d1=0.1 mm,AUC 0.60,(b) d2=0.2 mm,AUC 0.75, (c) d3=0.5 mm,AUC 0.96,(d) d4=1.0 mm,AUC 1.00

    隨著進(jìn)一步增大入射μ子通量,每個(gè)散射樣本的容量也會(huì)增大,從而對(duì)判定結(jié)果也可能產(chǎn)生影響。圖4中結(jié)果顯示,μ子通量增大到125cm?2時(shí),選擇適當(dāng)?shù)拈撝?,該分?lèi)器可近似完美識(shí)別最小寬度為0.5mm的狹縫(AUC約為1)。

    圖4 μ子通量為125 cm?2時(shí),分類(lèi)器對(duì)于帶有寬度di的水平狹縫的鈾板的識(shí)別性能ROC曲線 (a) d1=0.1 mm,AUC 0.55,(b) d2=0.2 mm,AUC 0.86, (c) d3=0.5 mm,AUC 1.00

    圖5表明,當(dāng)μ子通量增大到1000 cm?2時(shí),選擇適當(dāng)?shù)拈撝?,該分?lèi)器可較完美識(shí)別最小寬度為0.2 mm的狹縫(AUC約為0.99)。

    圖5 μ子通量為1000 cm?2時(shí),分類(lèi)器對(duì)于帶有寬度d的水平狹縫的鈾板的識(shí)別性能ROC曲線 (a) d1=0.1 mm,AUC 0.83,(b) d2=0.2 mm,AUC 0.99

    3.2 多種狹縫寬度可能性情況下的識(shí)別性能評(píng)估

    根據(jù)§3.1的結(jié)果,在μ子通量1000cm?2的條件下,采用正向參比法可完美分辨出1cm鈾板內(nèi)寬度最小可達(dá)0.2mm的狹縫。然而,在實(shí)際應(yīng)用中,除了需要對(duì)目標(biāo)對(duì)象是否產(chǎn)生狹縫進(jìn)行簡(jiǎn)單定性的Y/N判斷以外,目標(biāo)對(duì)象的結(jié)構(gòu)變化可能性比較多,即使對(duì)于同一類(lèi)型和位置的變化(如中心層的水平狹縫),往往同時(shí)存在多種尺度(如狹縫寬度)的可能性。為了評(píng)估這種情況下正向參比法對(duì)目標(biāo)對(duì)象結(jié)構(gòu)變化的識(shí)別性能,構(gòu)建如下場(chǎng)景:設(shè)計(jì)狹縫寬度分別為234的三個(gè)鈾板作為待測(cè)目標(biāo)對(duì)象,以狹縫寬度分別為234的鈾板各40組(共計(jì)120組)的模擬結(jié)果獲取的120個(gè)散射角分布樣本組成參比數(shù)據(jù)庫(kù),選取μ子通量為1000 cm?2,考察該情形下正向參比法對(duì)于不同狹縫寬度的識(shí)別能力。其結(jié)果如圖6所示。

    圖6 μ子通量1000 cm?2條件下,鈾板目標(biāo)對(duì)象狹縫寬度同時(shí)存在圖中三種可能時(shí),正向參比法對(duì)狹縫寬度di的目標(biāo)對(duì)象的識(shí)別性能ROC曲線 (a) d2=0.2 mm,AUC 0.99,(b) d3= 0.5 mm,AUC 1.00, (c) d4=1.0 mm,AUC 1.00

    由圖6可知,在μ子通量達(dá)到1000 cm?2的條件下,對(duì)于任意一種狹縫寬度為d(=2, 3, 4)的鈾板對(duì)象來(lái)說(shuō),其ROC曲線的AUC值(分別為0.99、1和1)都接近1的理想值。由此可見(jiàn),在該鈾板對(duì)象結(jié)構(gòu)變化的可能性被限定為上述三種寬度的水平狹縫的前提下,正向參比法可以較為理想地實(shí)現(xiàn)其中任意一種寬度的狹縫與另外兩種狹縫的區(qū)分與識(shí)別。這表明,在適當(dāng)?shù)南薅l件下,正向參比法可以實(shí)現(xiàn)目標(biāo)對(duì)象特定結(jié)構(gòu)變化尺度的準(zhǔn)確評(píng)估。

    4 結(jié)語(yǔ)

    相對(duì)于經(jīng)典μ子成像方法用于特殊對(duì)象結(jié)構(gòu)探測(cè)的技術(shù)難點(diǎn),本文提出的正向參比法更適用于初始結(jié)構(gòu)已知的特殊對(duì)象內(nèi)部結(jié)構(gòu)變化探測(cè),有望實(shí)現(xiàn)更好的空間分辨率和時(shí)間效率。從帶有0.1?1.0mm范圍內(nèi)不同寬度水平狹縫的1 cm厚度鈾板模型的GEANT4模擬結(jié)果來(lái)看,在合理的天然μ子通量條件下,該方法對(duì)于目標(biāo)對(duì)象的微小結(jié)構(gòu)變化尺度進(jìn)行準(zhǔn)確判定和評(píng)估,未來(lái)有望在復(fù)雜大型核對(duì)象的無(wú)損探測(cè)等相關(guān)領(lǐng)域取得顯著成效。

    1 Borozdin K, Hogan G, Morris C,. Surveillance: radiographic imaging with cosmic-ray muons[J]. Nature, 2003, 422(6929): 277. DOI: 10.1038/422277a.

    2 Schultz L, Borozdin K, Gomez J,. Image reconstruction and materialdiscriminationcosmic ray muon radiography[J]. Nuclear Instruments and Methods in Physics ResearchA, 2004, 519(3): 687?694. DOI: 10.1016/j.nima.2003.11.035.

    3 Morris C, Alexander C, Bacon J,. Tomographic imaging with cosmic ray muons[J]. Science & Global Security, 2008, 16(1?2): 37?53. DOI:10.1080/ 08929880802335758.

    4 Hagmann C, Lange D, Wright D. Monte Carlo simulation of proton-induced cosmic-ray cascades in the atmosphere (UCRL-TM-229452)[R]. USA: Lawrence Livermore National Laboratory, 2007.

    5 以恒冠. 面向特殊對(duì)象的宇宙線繆子散射成像方法研究[D]. 北京: 清華大學(xué), 2015. YI Hengguan. Research on cosmic ray muon scattering tomography methods for special objects inspection[D]. Beijing: Tsinghua University, 2015.

    6 Fawcett T. An introduction to ROC analysis[J]. Pattern Recognition Letters, 2006, 27: 861?874. DOI:10.1016/j. patrec.2005.10.010.

    7 鄒洪俠, 秦鋒, 程澤凱, 等. 二類(lèi)分類(lèi)器的ROC曲線生成算法[J]. 計(jì)算機(jī)技術(shù)與發(fā)展, 2009, 19(6): 109?112. ZOU Hongxia, QIN Feng, CHENG Zekai,. Algorithm for generating ROC curve of two-classifier[J]. Computer Technology & Development, 2009, 19(6): 109?112.

    Width evaluation of horizontal slit in uranium slabs using a forward comparison method based on cosmic-ray muon scattering

    XIAO Sa HE Weibo SHUAI Maobing LAN Mingcong WEI Mengfu

    (Institute of Materials, China Academy of Engineering Physics, Jiangyou 621907, China)

    Background: Muon tomography based on multiple Coulomb scattering (MCS) has been widely used in structural detection of nuclear objects. However, for some particular objects, the tomographic methodwould be thwarted by failing to imaging them. Purpose:In view of the technical difficulties in structural detection of those objects, a novel method named “Forward Comparison Method (FCM)” was proposed. Methods:Both the deformation and its scale were identified effectively by directly evaluating the evolution of muon scattering properties with respect to the object. This method, instead of the traditional inverse one, showed a methodological superiority in detection of nuclear objects with prior knowledge of the original structure.Results:Preliminary results based on a computational simulation demonstrated that FCM could be used to identify the horizontal slit in a uranium slab and evaluate its width on the magnitude of sub-millimeters, on condition that the muon flux was reasonable.Conclusion:The FCM method can be effectively applied in structural detection of particular nuclear objects of hollow structure on a satisfactory timescale.

    Muon, High-Z material, MCS, Non-parametric test, Receiver operating characteristic analysis

    TP312,TL81

    10.11889/j.0253-3219.2017.hjs.40.050202

    肖灑,男,1983年出生,2011年畢業(yè)于中國(guó)工程物理研究院,現(xiàn)為博士研究生,主要從事射線成像技術(shù)研究

    帥茂兵,E-mail: shuaimaobing@caep.cn

    2016-11-15,

    2017-03-17

    Supported by Science and Technology Development Fund of China Academy of Engineering Physics (No.2015B0103014)

    中國(guó)工程物理研究院科學(xué)技術(shù)發(fā)展基金(No.2015B0103014)資助

    XIAO Sa, male, born in 1983, graduated from China Academy of Engineering Physics in 2011, doctoral student, mainly engaged in the research of radiography

    SHUAI Maobing, E-mail: shuaimaobing@caep.cn

    2016-11-15, accepted date: 2017-03-17

    猜你喜歡
    通量分類(lèi)器寬度
    冬小麥田N2O通量研究
    BP-GA光照分類(lèi)器在車(chē)道線識(shí)別中的應(yīng)用
    馬屁股的寬度
    加權(quán)空-譜與最近鄰分類(lèi)器相結(jié)合的高光譜圖像分類(lèi)
    結(jié)合模糊(C+P)均值聚類(lèi)和SP-V-支持向量機(jī)的TSK分類(lèi)器
    緩釋型固體二氧化氯的制備及其釋放通量的影響因素
    紅細(xì)胞分布寬度與血栓的關(guān)系
    孩子成長(zhǎng)中,對(duì)寬度的追求更重要
    人生十六七(2015年5期)2015-02-28 13:08:24
    基于LLE降維和BP_Adaboost分類(lèi)器的GIS局部放電模式識(shí)別
    春、夏季長(zhǎng)江口及鄰近海域溶解甲烷的分布與釋放通量
    高清av免费在线| 国产精品一区二区精品视频观看| 另类精品久久| 国产欧美日韩一区二区三| 99国产综合亚洲精品| 免费观看a级毛片全部| 高清在线国产一区| 下体分泌物呈黄色| 国产精品国产av在线观看| 免费女性裸体啪啪无遮挡网站| 国产av又大| 制服人妻中文乱码| 国产精品98久久久久久宅男小说| 国产免费av片在线观看野外av| 51午夜福利影视在线观看| 亚洲伊人久久精品综合| 香蕉久久夜色| 亚洲精品国产区一区二| 欧美亚洲日本最大视频资源| 两个人免费观看高清视频| 亚洲精品中文字幕一二三四区 | 大型黄色视频在线免费观看| 国产精品一区二区在线观看99| 1024视频免费在线观看| 午夜成年电影在线免费观看| 国产精品一区二区精品视频观看| 精品一品国产午夜福利视频| 这个男人来自地球电影免费观看| 国产高清国产精品国产三级| 成人特级黄色片久久久久久久 | 精品国产一区二区久久| 丰满少妇做爰视频| 国产真人三级小视频在线观看| 午夜福利在线免费观看网站| 国产极品粉嫩免费观看在线| 色综合婷婷激情| 欧美日韩亚洲高清精品| 国产成人精品在线电影| 少妇的丰满在线观看| 咕卡用的链子| 日本欧美视频一区| 天天操日日干夜夜撸| av又黄又爽大尺度在线免费看| 亚洲欧美一区二区三区久久| 久久精品国产a三级三级三级| 国产成人精品在线电影| 久久久久久亚洲精品国产蜜桃av| 免费在线观看黄色视频的| 国产亚洲欧美精品永久| 国产精品一区二区精品视频观看| 国产欧美日韩一区二区三| 中亚洲国语对白在线视频| 欧美日韩亚洲国产一区二区在线观看 | 汤姆久久久久久久影院中文字幕| 久久久久精品国产欧美久久久| 午夜福利视频精品| www日本在线高清视频| 九色亚洲精品在线播放| 国产av一区二区精品久久| 交换朋友夫妻互换小说| 2018国产大陆天天弄谢| 国产视频一区二区在线看| 久久中文看片网| 欧美日韩成人在线一区二区| 我的亚洲天堂| 99久久人妻综合| 日本黄色日本黄色录像| 亚洲av日韩在线播放| 亚洲熟女精品中文字幕| 老司机午夜十八禁免费视频| 国产99久久九九免费精品| 黄网站色视频无遮挡免费观看| 亚洲欧美色中文字幕在线| 少妇裸体淫交视频免费看高清 | 亚洲熟女毛片儿| 国产欧美日韩一区二区三| 这个男人来自地球电影免费观看| 成人18禁高潮啪啪吃奶动态图| 久久国产精品影院| 亚洲男人天堂网一区| 久久午夜亚洲精品久久| netflix在线观看网站| 精品人妻1区二区| 久久人妻福利社区极品人妻图片| 美女主播在线视频| 国产精品欧美亚洲77777| 老司机在亚洲福利影院| 亚洲av美国av| 国产麻豆69| 亚洲av成人一区二区三| 97在线人人人人妻| av网站免费在线观看视频| 国产片内射在线| 精品国产亚洲在线| 亚洲精品粉嫩美女一区| 最新在线观看一区二区三区| 不卡av一区二区三区| 精品一区二区三卡| 在线观看免费日韩欧美大片| 亚洲成人手机| av欧美777| 国产黄色免费在线视频| 啦啦啦视频在线资源免费观看| 伊人久久大香线蕉亚洲五| 极品教师在线免费播放| 国产成人啪精品午夜网站| 久热爱精品视频在线9| 欧美激情 高清一区二区三区| 国产av又大| 一区二区三区国产精品乱码| 国产不卡一卡二| 亚洲欧美一区二区三区黑人| 日韩有码中文字幕| 大码成人一级视频| e午夜精品久久久久久久| 欧美日韩亚洲国产一区二区在线观看 | 国产一区二区激情短视频| 久久人妻熟女aⅴ| 桃红色精品国产亚洲av| 国产成人一区二区三区免费视频网站| 人妻久久中文字幕网| 久久99热这里只频精品6学生| 午夜福利影视在线免费观看| 国产日韩欧美视频二区| 99久久人妻综合| 十八禁网站免费在线| 别揉我奶头~嗯~啊~动态视频| 9191精品国产免费久久| 日韩 欧美 亚洲 中文字幕| 亚洲精品自拍成人| 真人做人爱边吃奶动态| 在线观看免费视频网站a站| 99riav亚洲国产免费| 黑人欧美特级aaaaaa片| 丝袜喷水一区| 青草久久国产| 精品国产国语对白av| 天天影视国产精品| 国产欧美日韩一区二区精品| 高清欧美精品videossex| av欧美777| 久久九九热精品免费| 久久 成人 亚洲| 亚洲美女黄片视频| 欧美精品一区二区免费开放| 不卡av一区二区三区| 国产一卡二卡三卡精品| av福利片在线| 777久久人妻少妇嫩草av网站| 精品乱码久久久久久99久播| 欧美日韩一级在线毛片| 下体分泌物呈黄色| 90打野战视频偷拍视频| 久久久精品94久久精品| 18禁观看日本| a级毛片在线看网站| 黄色怎么调成土黄色| 国产日韩欧美亚洲二区| 久久精品亚洲av国产电影网| 大片免费播放器 马上看| 搡老岳熟女国产| 亚洲精品自拍成人| 美国免费a级毛片| 国产精品久久久人人做人人爽| 男人操女人黄网站| 不卡av一区二区三区| 国产精品免费一区二区三区在线 | 日韩中文字幕视频在线看片| 亚洲精品国产区一区二| 精品福利永久在线观看| 午夜免费鲁丝| 亚洲国产毛片av蜜桃av| 一区在线观看完整版| 热re99久久国产66热| 国产真人三级小视频在线观看| 高清毛片免费观看视频网站 | 国产精品1区2区在线观看. | 日本a在线网址| 久久ye,这里只有精品| 91字幕亚洲| 亚洲av日韩在线播放| 久久人人爽av亚洲精品天堂| 亚洲第一av免费看| 99香蕉大伊视频| 夜夜骑夜夜射夜夜干| 久久天堂一区二区三区四区| 精品久久蜜臀av无| 999久久久国产精品视频| 51午夜福利影视在线观看| 精品第一国产精品| 亚洲第一av免费看| 国产精品秋霞免费鲁丝片| av不卡在线播放| 欧美激情 高清一区二区三区| aaaaa片日本免费| 欧美日韩精品网址| 黑丝袜美女国产一区| 法律面前人人平等表现在哪些方面| 美女视频免费永久观看网站| 最近最新免费中文字幕在线| 老司机深夜福利视频在线观看| 亚洲,欧美精品.| 十分钟在线观看高清视频www| 麻豆成人av在线观看| 亚洲av电影在线进入| 欧美国产精品一级二级三级| 免费黄频网站在线观看国产| 中文欧美无线码| 香蕉久久夜色| 国产三级黄色录像| 2018国产大陆天天弄谢| 国产男靠女视频免费网站| 这个男人来自地球电影免费观看| 高清欧美精品videossex| 国产精品二区激情视频| 国产精品 欧美亚洲| 午夜视频精品福利| 1024视频免费在线观看| 不卡av一区二区三区| 国产一区二区三区综合在线观看| 黑人猛操日本美女一级片| 嫩草影视91久久| 国产精品 国内视频| 色综合婷婷激情| kizo精华| 欧美中文综合在线视频| av线在线观看网站| 国产熟女午夜一区二区三区| 亚洲一区中文字幕在线| 色94色欧美一区二区| 久久久精品94久久精品| 美女午夜性视频免费| 亚洲人成电影免费在线| 热re99久久精品国产66热6| 中文字幕人妻丝袜制服| 黄色毛片三级朝国网站| 人人妻人人爽人人添夜夜欢视频| 多毛熟女@视频| 精品人妻1区二区| 久久精品亚洲av国产电影网| 亚洲午夜精品一区,二区,三区| 欧美日本中文国产一区发布| 久久ye,这里只有精品| 色在线成人网| 69精品国产乱码久久久| 黄色视频在线播放观看不卡| 肉色欧美久久久久久久蜜桃| 99国产精品99久久久久| 欧美激情 高清一区二区三区| 狂野欧美激情性xxxx| 成人av一区二区三区在线看| 黄片大片在线免费观看| 亚洲欧美精品综合一区二区三区| 中文字幕另类日韩欧美亚洲嫩草| 99国产精品99久久久久| 亚洲精品美女久久av网站| svipshipincom国产片| av天堂在线播放| 久久性视频一级片| 亚洲欧美日韩另类电影网站| 大型av网站在线播放| 成人18禁在线播放| 在线观看66精品国产| 在线观看一区二区三区激情| 国产精品秋霞免费鲁丝片| 午夜日韩欧美国产| 国产在线一区二区三区精| 欧美日韩亚洲综合一区二区三区_| 日本vs欧美在线观看视频| 午夜日韩欧美国产| 91成人精品电影| 丝袜美腿诱惑在线| 精品午夜福利视频在线观看一区 | 国产精品久久久人人做人人爽| 中亚洲国语对白在线视频| 亚洲欧洲精品一区二区精品久久久| 一区二区三区激情视频| 在线天堂中文资源库| av福利片在线| 精品高清国产在线一区| 国产精品久久久av美女十八| 桃花免费在线播放| 12—13女人毛片做爰片一| 亚洲精华国产精华精| 亚洲自偷自拍图片 自拍| 欧美激情极品国产一区二区三区| 久久午夜综合久久蜜桃| 久久毛片免费看一区二区三区| 久久久久久亚洲精品国产蜜桃av| 97在线人人人人妻| 精品一区二区三区四区五区乱码| 自线自在国产av| 如日韩欧美国产精品一区二区三区| a级毛片黄视频| 亚洲性夜色夜夜综合| 怎么达到女性高潮| 国产精品免费一区二区三区在线 | 精品国产乱子伦一区二区三区| 久热这里只有精品99| 成年女人毛片免费观看观看9 | 久久亚洲精品不卡| 久久久久国产一级毛片高清牌| 99国产极品粉嫩在线观看| 欧美日韩视频精品一区| 激情视频va一区二区三区| 精品福利永久在线观看| 激情在线观看视频在线高清 | 亚洲成国产人片在线观看| 狠狠精品人妻久久久久久综合| 黄片大片在线免费观看| 亚洲欧美日韩另类电影网站| 视频在线观看一区二区三区| 亚洲精品成人av观看孕妇| 最近最新中文字幕大全免费视频| 在线十欧美十亚洲十日本专区| 极品教师在线免费播放| 精品久久久精品久久久| 国产精品自产拍在线观看55亚洲 | 精品国产亚洲在线| 国产精品免费一区二区三区在线 | 亚洲人成伊人成综合网2020| 久久婷婷成人综合色麻豆| 在线观看免费日韩欧美大片| 成人特级黄色片久久久久久久 | 免费在线观看视频国产中文字幕亚洲| 狠狠婷婷综合久久久久久88av| 成人精品一区二区免费| 欧美日本中文国产一区发布| 亚洲色图综合在线观看| 亚洲久久久国产精品| av天堂久久9| 亚洲精品久久成人aⅴ小说| 中文字幕色久视频| 久久精品国产99精品国产亚洲性色 | 人人妻人人爽人人添夜夜欢视频| 9色porny在线观看| 亚洲第一青青草原| 国产精品一区二区免费欧美| 国产精品二区激情视频| 精品少妇久久久久久888优播| 亚洲第一青青草原| 我的亚洲天堂| 在线天堂中文资源库| 人人妻人人添人人爽欧美一区卜| 99精品欧美一区二区三区四区| 久久影院123| 久久精品国产亚洲av香蕉五月 | 麻豆成人av在线观看| 老鸭窝网址在线观看| 欧美日韩国产mv在线观看视频| 日韩熟女老妇一区二区性免费视频| 日日夜夜操网爽| 国产欧美日韩精品亚洲av| 国产精品久久电影中文字幕 | 男女床上黄色一级片免费看| 欧美亚洲 丝袜 人妻 在线| 啦啦啦视频在线资源免费观看| 久久久精品免费免费高清| 日本vs欧美在线观看视频| 999久久久精品免费观看国产| 国产成人精品无人区| 国产1区2区3区精品| 狠狠精品人妻久久久久久综合| 一级,二级,三级黄色视频| 9191精品国产免费久久| 国产精品二区激情视频| 国产单亲对白刺激| 999久久久精品免费观看国产| 亚洲国产欧美一区二区综合| 午夜91福利影院| 欧美在线黄色| 日韩大码丰满熟妇| 男男h啪啪无遮挡| 亚洲精品美女久久久久99蜜臀| 欧美久久黑人一区二区| 99国产精品免费福利视频| 国产成人av激情在线播放| 免费在线观看完整版高清| 国产在线精品亚洲第一网站| 国产在线视频一区二区| 窝窝影院91人妻| av有码第一页| 亚洲av成人一区二区三| 中文字幕另类日韩欧美亚洲嫩草| 欧美激情久久久久久爽电影 | 国产日韩欧美视频二区| 亚洲天堂av无毛| 999精品在线视频| 别揉我奶头~嗯~啊~动态视频| 女人被躁到高潮嗷嗷叫费观| 菩萨蛮人人尽说江南好唐韦庄| 亚洲伊人色综图| 国产欧美日韩综合在线一区二区| 精品亚洲乱码少妇综合久久| 一本色道久久久久久精品综合| 色在线成人网| 法律面前人人平等表现在哪些方面| av视频免费观看在线观看| 欧美国产精品一级二级三级| 18禁国产床啪视频网站| 日日夜夜操网爽| 国产精品.久久久| 王馨瑶露胸无遮挡在线观看| 狠狠婷婷综合久久久久久88av| 色婷婷久久久亚洲欧美| 久久天躁狠狠躁夜夜2o2o| 中文亚洲av片在线观看爽 | 国产成人av激情在线播放| 捣出白浆h1v1| 色综合婷婷激情| 香蕉丝袜av| 久久精品亚洲熟妇少妇任你| 久久精品成人免费网站| 久久久精品区二区三区| 美女主播在线视频| 99久久人妻综合| 亚洲精品乱久久久久久| 午夜福利乱码中文字幕| 菩萨蛮人人尽说江南好唐韦庄| 国产精品亚洲一级av第二区| 欧美久久黑人一区二区| 中文字幕制服av| 色精品久久人妻99蜜桃| 精品卡一卡二卡四卡免费| 精品视频人人做人人爽| 国产免费福利视频在线观看| 亚洲人成电影观看| 成年人免费黄色播放视频| 大码成人一级视频| 久久久国产欧美日韩av| 久久久国产成人免费| 亚洲欧美日韩另类电影网站| 免费看十八禁软件| 国产精品久久久久久精品古装| 欧美日韩一级在线毛片| 波多野结衣av一区二区av| 欧美一级毛片孕妇| 80岁老熟妇乱子伦牲交| 黄色a级毛片大全视频| 女人久久www免费人成看片| 美国免费a级毛片| 纵有疾风起免费观看全集完整版| 亚洲中文字幕日韩| 如日韩欧美国产精品一区二区三区| 国产高清视频在线播放一区| 成人国语在线视频| 国产精品国产av在线观看| 怎么达到女性高潮| 99精国产麻豆久久婷婷| 欧美精品亚洲一区二区| 国产熟女午夜一区二区三区| 极品少妇高潮喷水抽搐| 免费久久久久久久精品成人欧美视频| 日韩欧美三级三区| 美女福利国产在线| 国产免费福利视频在线观看| 亚洲中文字幕日韩| 久久这里只有精品19| 国产又爽黄色视频| 久久人人爽av亚洲精品天堂| 成在线人永久免费视频| 亚洲精品一卡2卡三卡4卡5卡| 国产av又大| 十八禁高潮呻吟视频| 成人三级做爰电影| 久久中文看片网| 日日摸夜夜添夜夜添小说| 久久人妻福利社区极品人妻图片| 侵犯人妻中文字幕一二三四区| 别揉我奶头~嗯~啊~动态视频| videos熟女内射| 国产男女内射视频| 国产免费av片在线观看野外av| 人妻久久中文字幕网| 免费少妇av软件| 人妻久久中文字幕网| 精品一品国产午夜福利视频| 人人妻人人添人人爽欧美一区卜| 国产精品av久久久久免费| 国产精品国产高清国产av | 日本撒尿小便嘘嘘汇集6| 精品一品国产午夜福利视频| 精品一区二区三区av网在线观看 | 国产99久久九九免费精品| 精品国产一区二区三区久久久樱花| 精品福利观看| 久9热在线精品视频| 男女午夜视频在线观看| 黄色丝袜av网址大全| 亚洲中文av在线| 老司机深夜福利视频在线观看| 汤姆久久久久久久影院中文字幕| 国产精品影院久久| 欧美中文综合在线视频| 在线天堂中文资源库| 黄色怎么调成土黄色| 国产精品自产拍在线观看55亚洲 | 久久久欧美国产精品| 99精品在免费线老司机午夜| 男女免费视频国产| 老汉色∧v一级毛片| 18禁裸乳无遮挡动漫免费视频| 久久久精品免费免费高清| 啦啦啦在线免费观看视频4| 久久人妻熟女aⅴ| 建设人人有责人人尽责人人享有的| 美女主播在线视频| 国产又色又爽无遮挡免费看| 我的亚洲天堂| 91av网站免费观看| 亚洲成人免费电影在线观看| 亚洲成a人片在线一区二区| 国产精品国产av在线观看| 久久 成人 亚洲| 亚洲一码二码三码区别大吗| 欧美 日韩 精品 国产| 欧美一级毛片孕妇| 捣出白浆h1v1| 国产精品国产高清国产av | 久久国产亚洲av麻豆专区| 纵有疾风起免费观看全集完整版| 亚洲国产成人一精品久久久| 777米奇影视久久| 亚洲男人天堂网一区| 国产在视频线精品| 久久精品亚洲av国产电影网| 亚洲成人免费av在线播放| 国产一区二区 视频在线| 亚洲九九香蕉| 久久久久久免费高清国产稀缺| 桃花免费在线播放| 人妻一区二区av| 一区二区三区乱码不卡18| 这个男人来自地球电影免费观看| 女人高潮潮喷娇喘18禁视频| 91大片在线观看| 国内毛片毛片毛片毛片毛片| 久久久久久免费高清国产稀缺| 国产单亲对白刺激| 一二三四社区在线视频社区8| 欧美人与性动交α欧美软件| 亚洲天堂av无毛| 国产又爽黄色视频| 一级毛片精品| 亚洲精华国产精华精| 男人舔女人的私密视频| 午夜久久久在线观看| 亚洲av日韩精品久久久久久密| 亚洲精品一卡2卡三卡4卡5卡| 精品一区二区三区四区五区乱码| 国产91精品成人一区二区三区 | 亚洲精品自拍成人| 男女午夜视频在线观看| 国产在线精品亚洲第一网站| 日本a在线网址| 久久婷婷成人综合色麻豆| 大片电影免费在线观看免费| 日韩熟女老妇一区二区性免费视频| 九色亚洲精品在线播放| 日本撒尿小便嘘嘘汇集6| 日韩欧美三级三区| 亚洲va日本ⅴa欧美va伊人久久| 69精品国产乱码久久久| 久久精品aⅴ一区二区三区四区| 色播在线永久视频| 丝袜人妻中文字幕| 天天躁狠狠躁夜夜躁狠狠躁| 人人妻人人添人人爽欧美一区卜| 亚洲一码二码三码区别大吗| 日韩制服丝袜自拍偷拍| 最新美女视频免费是黄的| 黄频高清免费视频| 美女视频免费永久观看网站| 午夜久久久在线观看| 天天添夜夜摸| 国产无遮挡羞羞视频在线观看| 精品午夜福利视频在线观看一区 | 黄色a级毛片大全视频| 亚洲一码二码三码区别大吗| 国产深夜福利视频在线观看| 欧美久久黑人一区二区| 最近最新中文字幕大全免费视频| 欧美日韩亚洲高清精品| 黄色丝袜av网址大全| 91av网站免费观看| 中文字幕制服av| 亚洲精品国产色婷婷电影| 肉色欧美久久久久久久蜜桃| 在线亚洲精品国产二区图片欧美| svipshipincom国产片| 国内毛片毛片毛片毛片毛片| 悠悠久久av| 一级毛片电影观看| 国产有黄有色有爽视频| 少妇粗大呻吟视频| 国产人伦9x9x在线观看| 免费女性裸体啪啪无遮挡网站| 三上悠亚av全集在线观看| www.熟女人妻精品国产| 叶爱在线成人免费视频播放| 少妇粗大呻吟视频| 亚洲欧洲精品一区二区精品久久久| 男女之事视频高清在线观看| 999久久久国产精品视频| 最新美女视频免费是黄的| 国产福利在线免费观看视频| 人人妻,人人澡人人爽秒播| 色视频在线一区二区三区| 国产精品 国内视频| 免费黄频网站在线观看国产| 啦啦啦视频在线资源免费观看| 高清毛片免费观看视频网站 | 一二三四在线观看免费中文在| 在线av久久热| av有码第一页|