焦若男，劉 琨，孔繁藝，王 婷，韓 雪，李泳江，孫長森

大連理工大學(xué)光電工程與儀器科學(xué)學(xué)院，遼寧 大連 116000

引 言

隨著塑料制品的廣泛應(yīng)用，微塑料逐漸成為地球上一種新興污染物。 微塑料是指尺寸小于5mm的塑料碎片。 由于其形成原因不同，微塑料被分為初級微塑料和次級微塑料，初級微塑料是指直接被制成微觀尺寸的塑料，如化妝品中含有的塑料顆粒、藥物載體等；次級微塑料是指大的塑料碎片經(jīng)過物理、化學(xué)、生物等過程分解所得到的產(chǎn)物。 微塑料的體積小，表面積大，吸附能力強，擴散范圍廣。

起初人們對微塑料的研究主要針對海洋生態(tài)系統(tǒng)。 但近年來研究表明中國黃渤海海域沉積物中存在約100個·kg-1(干重)左右的微塑料[1]。 不僅在我國海域，在世界各大海域[2]、南極羅斯海[3]、甚至北極海冰[4]中都有微塑料的蹤跡。 海洋遭到微塑料的襲擊，海洋生物和海產(chǎn)品也不例外。 微塑料不僅被發(fā)現(xiàn)存在于蝦的體內(nèi)，在溫哥華島東海岸的幼年奇努克鮭魚體內(nèi)[5]、伊比利亞群島常見的海豚胃中[6]均有微塑料的蹤跡。 歐洲鱈魚、螃蟹、海鳥、南極企鵝、雙殼類生物以及市場上購買的冷凍魚中也發(fā)現(xiàn)了或多或少的微塑料，甚至日常生活中不可或缺的海鹽中也有微塑料的蹤跡。 雖然微塑料廣泛存在于海洋生態(tài)系統(tǒng)中，但大部分的微塑料都來源于陸地[7]。 耕地是生產(chǎn)糧食必不可少的，但有的高原耕地卻也存在嚴重的微塑料污染[8]，這就可能導(dǎo)致種植出的糧食中含有微塑料。 2019年，Panebianco等首次在陸生蝸牛體內(nèi)發(fā)現(xiàn)了微塑料[9]。 諸如瓶裝水、罐裝魚等產(chǎn)品中檢測出微塑料的報道表明微塑料這一新興污染物已經(jīng)普遍存在于我們身邊。

雖然小于20 μm的微塑料不能在組織中轉(zhuǎn)移，但更小尺寸的微塑料被吸收轉(zhuǎn)移的概率會大大增加，如2 μm的微塑料就可被牡蠣吸收[10]。 人類位于生物鏈的最頂端，微塑料最終會在人體內(nèi)聚集，引發(fā)呼吸道過敏、呼吸困難、咳嗽等癥狀，進而引發(fā)哮喘、心血管類疾病，甚至可能引發(fā)癌癥，對人體健康產(chǎn)生危害。 微塑料污染日益嚴重，對生物及人類造成的影響也是不可估量的，因此研究微塑料的來源、傳播以及治理問題是刻不容緩的。

目前，人們對于微塑料的提取、檢測及鑒定有幾種不同的方法。 對于微塑料的提取大多數(shù)采用的是浮選法、密度分離法及過濾法，提取出微塑料后再利用顯微鏡進行目視觀察，或運用拉曼光譜技術(shù)[11]、傅里葉紅外光譜技術(shù)[12]、熱解氣相色譜與質(zhì)譜法[13]、高光譜成像等進行檢測與鑒定。 土壤樣品中的微塑料可采用浮選法進行提取，然后利用顯微鏡觀察結(jié)合加熱法進行鑒定，由于土壤樣品中的微塑料受熱熔化后會轉(zhuǎn)化為圓形透明顆粒，而有機物、硅酸鹽等其他雜質(zhì)在受熱后不發(fā)生變化，因此利用加熱法進行處理，再利用顯微鏡進行目視觀察，對土壤中的微塑料的提取率可高達90%[14]。 對環(huán)境中微塑料的提取一般還較為簡單，而提取生物體內(nèi)的微塑料則需要更復(fù)雜的樣品處理，首先要消解生物組織中的有機物，可向樣品中加入H2O2，HNO3[11]和KOH溶液[12]等試劑達到目的，然后再采用過濾的方式進行提取，最后將過濾所得的濾渣放在高倍顯微鏡下觀察[5]，或結(jié)合拉曼光譜、傅里葉紅外光譜、高光譜成像等方法進行識別與鑒定。

可見在提取微塑料的過程中不論是哪種方法，都需要較長時間的等待，甚至需要復(fù)雜的樣品處理。 諸多檢測方法中，目視檢測對于檢測大塊的微塑料比較直觀，但受主觀因素的影響較多，而對于小尺寸微塑料的鑒別，拉曼光譜檢測與傅里葉紅外光譜檢測有一定的優(yōu)勢，但也需要對可疑顆粒進行逐個分析。 鑒于此，開發(fā)一種快速、準(zhǔn)確鑒別環(huán)境中是否含有微塑料的技術(shù)是有必要的。 因此提出使用多通道圖像采集的方法檢測環(huán)境中的微塑料，并結(jié)合形態(tài)學(xué)分析與中值濾波的算法直觀顯示出微塑料在環(huán)境中的分布。

選取海水及沙子為研究對象，由于聚苯乙烯和聚乙烯是世界上最常用的塑料，也是海洋環(huán)境中最常見的塑料，所以我們選擇將聚苯乙烯微球摻入到海水及沙子中模擬被微塑料污染的海水和沙子，在不進行任何處理的情況下，利用多通道圖像采集系統(tǒng)同時獲得樣品白光通道下的圖像及相干反斯托克斯拉曼散射(coherent anti-stokes Raman scattering, CARS)光譜成像信號，通過CARS光譜成像結(jié)合形態(tài)學(xué)分析與中值濾波算法，直觀地顯示出海水及沙子中的聚苯乙烯微球的分布，同時與拉曼光譜檢測進行了對照。

1 原 理

1.1 CARS光譜成像原理

1965 年，Terhune和Maker等發(fā)現(xiàn)了CARS現(xiàn)象。 CARS光譜成像是一種基于非線性光學(xué)與分子振動的四波混頻過程，在CARS光譜成像過程中，通常使用泵浦光(頻率為ωp)和斯托克斯光(頻率為ωs)兩束光同時激發(fā)樣品，本研究中采用正向CARS，當(dāng)兩者頻率的差值正好等于所要探測的化學(xué)鍵的振動頻率ωR(即共振條件ωR=ωp-ωs)，同時又滿足匹配條件Ka=2Kp±Ks時，會激發(fā)出頻率為ωas=2ωp-ωs的反斯托克斯光(圖1)。 通過調(diào)節(jié)延時線使兩束光達到匹配從而得到反斯托克斯信號。 CARS光譜成像作為一種成像手段，具有非侵入性、非破壞性、實時性且可三維成像的特點，常被應(yīng)用于化學(xué)、材料、生物和醫(yī)學(xué)等領(lǐng)域，如利用CARS對活體小鼠中C—H化學(xué)鍵進行成像[15]。 CARS光譜成像不需要特殊制備或處理樣品，因此把它用來檢測環(huán)境中的微塑料。

圖1 CARS示意圖Fig.1 The schematic of CARS process

1.2 算法分析

圖2 算法流程圖Fig.2 Algorithmic framework

2 實驗部分

2.1 儀器及材料

使用RENISHAW(INVIA)拉曼光譜儀(激發(fā)波長532 nm，輸出功率5.15 mW)、Olympus FV1200MPE系統(tǒng)的多光子激光掃描顯微鏡(激發(fā)波長800 nm，泵浦光功率149 mW，斯托克斯光功率81 mW)、電子天平(上海梅特勒JB/T)、高速離心機(TG16-WS)。

材料為直徑10 μm，濃度為0.05 g·mL-1的聚苯乙烯微球(上海麥克林生化科技有限公司)、海水(采集于大連市東港)、沙子(采集于中國大連星海浴場)。

2.2 方法

將10 μL的聚苯乙烯微球與3 mL海水均勻混合作為海水的實驗樣品。 將聚苯乙烯微球加去離子水稀釋至0.001 g·mL-1，取0.5 mL加入到2 g干沙中，混合均勻(用離心機以8 000 r·min-1離心3 min)作為沙子的實驗樣品。

拉曼光譜在檢測環(huán)境中的微塑料方面已經(jīng)是一種成熟的手段，因此將拉曼光譜檢測作為一種對照方法。 對樣品進行拉曼光譜檢測與多通道的圖像采集，并進行分析。

采用Origin9.0軟件繪制曲線圖，采用FV10-ASW 4.2 Viewer進行多通道采集圖像的簡單處理，采用Matlab R2017a進行算法處理。

3 結(jié)果與討論

3.1 聚苯乙烯微球的拉曼光譜與多通道圖像采集

為了能夠準(zhǔn)確辨別海水及沙子中的聚苯乙烯微球，首先對聚苯乙烯進行了拉曼光譜檢測及多通道圖像采集。 聚苯乙烯在998，1 597和3 047 cm-1等位置處有明顯的拉曼頻移，且拉曼頻移信號的強度隨著微球濃度的降低而減弱(圖3)。 選取聚苯乙烯中由苯環(huán)內(nèi)碳原子之間非對稱收縮引起的1 597 cm-1處的拉曼頻移進行多通道圖像采集。 在多通道圖像采集中使用800 nm的飛秒激光作為激發(fā)光源，通過分光鏡將激光分為兩束，其中一束光作為泵浦光，另一束光經(jīng)過光子晶體拓寬為一束寬譜帶的光，兩束光同時照射在樣品上，激發(fā)出反斯托克斯信號，通過放置合適的濾光片，得到特定波數(shù)的CARS光譜成像信號，此處采用的為正向的。 在CARS光譜成像通道中放置(710±6.5) nm的濾光片可獲取1 597 cm-1拉曼頻移處的信號，調(diào)節(jié)延時線使其相位匹配，得到了聚苯乙烯的以1 597 cm-1為中心的多通道圖像采集信號(圖4)，其CARS光譜成像信號與白光通道下聚苯乙烯的輪廓可以一一對應(yīng)，且不會隨著小球的濃度變化而發(fā)生變化。

圖3 不同濃度聚苯乙烯微球的拉曼光譜圖Fig.3 Raman spectrum of polystyreneat different concentrations

3.2 海水中微塑料的檢測

在對海水中的聚苯乙烯進行多通道圖像采集時，信號明顯且與白光通道下觀察到的聚苯乙烯輪廓一一對應(yīng)。 由于激光照射加快了海水蒸發(fā)的過程并析出鹽顆粒，海水發(fā)生了結(jié)晶，雖然結(jié)晶部位有微弱的CARS光譜成像信號，但經(jīng)過一次開運算算法，可以直觀地看到聚苯乙烯在海水中的分布(圖5)。

圖4 以1 597 cm-1為中心的聚苯乙烯微球的多通道成像(a)： 白光通道成像；(b): CARS光譜成像Fig.4 Multi-channel imaging of polystyrenecentered at 1 597 cm-1(a)： The white light imaging；(b): The CARS imaging

圖5 以1 597 cm-1為中心的聚苯乙烯在海水中的多通道成像(a)： 白光通道成像；(b): CARS光譜成像Fig.5 Multi-channel imaging of polystyrenecentered at 1 597 cm-1 in seawater(a)： The white light imaging；(b): The CARS imaging

適當(dāng)增加多通道圖像采集的曝光時間，可以使以1 597 cm-1為中心的聚苯乙烯信號更加清晰(圖6)。 雖然結(jié)晶部位的信號也得到了增強，但結(jié)晶的信號普遍位于結(jié)晶部位附近，對檢測結(jié)果影響不大。 進行長時間的連續(xù)掃描時，以1 597cm-1為中心的聚苯乙烯信號強度減弱，原因可能是長時間的照射，激光的熱量過高導(dǎo)致小球熔化，如圖7所示，箭頭所指位置為含有聚苯乙烯的位置。 可以看出，多通道圖像采集在檢測海水中聚苯乙烯時不受其他成分干擾，可以快速直觀地看到聚苯乙烯的分布情況，雖然隨著照射時間增長，局部熱量導(dǎo)致水分蒸發(fā)，海水結(jié)晶，會對成像有一定的影響，CARS系統(tǒng)的單幅掃描成像時間可以在幾秒鐘完成，完全可以避免海水結(jié)晶的干擾。

圖6 不同曝光時間下海水中聚苯乙烯的CARS光譜成像信號(a)： 曝光時間2 μs·像素-1；(b): 曝光時間8 μs·像素-1

圖7 連續(xù)掃描下海水中聚苯乙烯的CARS光譜成像信號Fig.7 CARS imaging of polystyrene in seawaterunder continuous scanning

3.3 沙子中的微塑料檢測

在對沙子中的微塑料進行檢測之前，由于沙子中有較大顆粒，其大小遠遠超出聚苯乙烯微球的大小，為了避免沙子本身信號形成的干擾，首先對沙子進行了檢測。

為了檢測沙子與聚苯乙烯是否有重疊的拉曼頻移，我們對沙子進行了拉曼光譜檢測。 由于沙子普遍有較強的熒光信號，在700～3 200 cm-1波數(shù)范圍內(nèi)只有一個大的熒光峰，沒有明顯的拉曼頻移。 對所得拉曼光譜圖進行了去基線處理，取100～700 cm-1波數(shù)范圍的拉曼光譜，如圖8所示(為了直觀地區(qū)分沙粒的顏色，圖8中同時展示了顯微鏡下沙粒的形態(tài))。 沙子在479 cm-1處有較小的拉曼頻移，隨著沙粒顏色逐漸變淺，該拉曼頻移的信號強度逐漸增強，同時在125，201和461 cm-1位置處出現(xiàn)了新的拉曼頻移。 由于沙子主要是由礦物和微小的巖石碎片組成，且具體成分因地而異，我們認為拉曼頻移及強度的差異是由于不同沙粒組成成分及含量的不同導(dǎo)致的。 在對沙子進行多通道圖像采集時，淺色的接下來對沙子與聚苯乙烯混合后的樣品進行拉曼光譜檢測，沿一直線取十個點依次進行拉曼光譜檢測，沙子有強烈的熒光，對所得數(shù)據(jù)進行去基線處理后，發(fā)現(xiàn)只有在激光聚焦在聚苯乙烯存在的位置時，才能檢測到微弱的信號(圖10)。 對沙子中的單個聚苯乙烯進行檢測，聚苯乙烯的大部分信號會被沙子的強熒光遮住，進行去基線及平滑處理后，可以看到較小的信號(圖11)。 由此可見，運用拉曼光譜檢測的方法來檢測沙子中的聚苯乙烯具有一定的局限性，若要準(zhǔn)確的測得沙子中是否含有聚苯乙烯就要進行拉曼逐點掃描成像，這一過程非常耗時，且沙子的熒光對聚苯乙烯的檢測有較大影響，需要進一步的數(shù)據(jù)處理才能得到想要的結(jié)果。

圖8 不同沙粒的拉曼光譜圖Fig.8 Raman spectrum of different sand

沙粒幾乎沒有CARS光譜成像信號，而深色沙粒有較強的CARS光譜成像信號(圖9)，這也進一步證實了不同顏色的沙粒組成成分的不同。 對于細小沙粒，CARS光譜成像信號多為小且離散的信號，這與以1 597 cm-1為中心的聚苯乙烯的信號相比有較大的差別。

圖9 沙粒的多通道成像 (a)，(b) 分別為透明沙粒的白光及CARS光譜成像信號； (c)，(d) 分別為深色沙粒的白光及CARS光譜成像信號Fig.9 Multi-channel imaging of sand

對其進行多通道圖像采集，與海水的多通道圖像采集結(jié)果不同，沙子對聚苯乙烯的檢測有一定的干擾。 如圖11所示，CARS光譜成像信號中有大且聚集信號和一些點狀離散信號。 將白光通道成像與CARS光譜成像信號進行對比，發(fā)現(xiàn)較大且聚集的信號與聚苯乙烯的輪廓相符(圖12)，其中黃色圓圈標(biāo)記的位置是可以確定為聚苯乙烯的位置，點狀離散信號可以確定為沙子本身的信號，在該檢測過程中為干擾信號。 由于在CARS光譜成像信號中并不能明顯地從沙子中辨別出聚苯乙烯，對圖像進行算法處理。 因為聚苯乙烯與沙子的信號尺寸差距較大，首先采用濾波算法進行處理。 然而經(jīng)過3×3的中值濾波算法后仍然存在大范圍的沙子信號的干擾[圖13(a)]。 由于聚苯乙烯微球有一定的體積，在進行多通道圖像采集過程中可能存在不同的聚苯乙烯微球在不同的焦平面從而使得聚苯乙烯的信號不夠均勻，邊界不夠平滑，這可能導(dǎo)致在濾波過程中將有效信息濾除，因此提出運用先腐蝕后膨脹的開運算算法，得到的結(jié)果仍然不夠直觀[圖13(b)]。 經(jīng)過多次試驗，發(fā)現(xiàn)腐蝕膨脹與中值濾波相結(jié)合的算法可以得到較好的結(jié)果。 首先通過3×3的中值濾波結(jié)合圖像腐蝕算法去除掉沙子的點狀信號，再經(jīng)過中值濾波結(jié)合膨脹算法填充聚苯乙烯的信號，使其信號均勻，邊界平滑，便于觀察。 經(jīng)此過程可以濾除掉大部分沙子的點狀離散信號而保留聚苯乙烯的大且聚集信號。 雖然所得結(jié)果不夠精確，但能夠直觀地看到聚苯乙烯在沙子中的分布[圖13(c)]。

圖10 含有聚苯乙烯的沙子中不同位置的拉曼光譜圖Fig.10 Raman spectrum of different locationsin sand containing polystyrene

圖12 以1 597 cm-1為中心的聚苯乙烯在沙子中的CARS光譜成像信號Fig.12 CARS imaging of polystyrene insand centered at 1 597 cm-1

圖13 不同算法處理后的以1 597 cm-1為中心的聚苯乙烯在沙子中的CARS光譜成像信號(a)： 僅中值濾波；(b): 開運算；(c): 濾波與開運算結(jié)合Fig.13 CARS imaging of polystyrene in sand after processing by algorithm(a)： Median filter; (b): Erode and dilation; (c): Corrosion and expansion combined with median filtering

4 結(jié) 論

運用多通道圖像采集的方法快速檢測出海水及沙子中的微米級聚苯乙烯。 對于毫米級的微塑料，該檢測手段的準(zhǔn)確率將大幅提高。 對于不同種類的微塑料，如聚乙烯、聚丙烯等，可以以不同的拉曼頻移為中心，選用對應(yīng)的濾光片進行成像。 與目前其他檢測手段相比，在不進行預(yù)處理的情況下，多通道圖像采集對于微米量級的微塑料可以快速識別，也展現(xiàn)出了其對環(huán)境中的微塑料檢測方面應(yīng)用潛力。