莫俊超,姚洪偉,吳孝槐,郭丹丹,趙方舟,舒耀皋*

頂空被動加標(biāo)系統(tǒng)的多介質(zhì)動力學(xué)模型

莫俊超1,2,姚洪偉1,2,吳孝槐3,郭丹丹1,2,趙方舟1,2,舒耀皋1,2*

(1.上?；ぱ芯吭河邢薰?上海 200062；2.上?；瘜W(xué)品公共安全工程技術(shù)研究中心,上海 200062；3.中國合格評定國家認(rèn)可委員會,北京 100062)

使用多介質(zhì)模型建模方法,構(gòu)建了針對疏水性液體物質(zhì)的頂空被動加標(biāo)系統(tǒng)的動力學(xué)模型,以浮萍生長抑制試驗系統(tǒng)和魚類胚胎毒性試驗系統(tǒng)為例,得出了3種不同性質(zhì)的疏水性液體物質(zhì)正壬烷、十甲基環(huán)五硅氧烷(D5)和鄰苯二甲酸(2-乙基己基)酯(DEHP)在頂空被動加標(biāo)系統(tǒng)中的濃度變化,并使用實測濃度對模型進行驗證.結(jié)果表明,正壬烷和D5在兩類系統(tǒng)中平衡24h后,水中濃度便已達(dá)到完全平衡時濃度的80%以上,而DEHP的平衡時間很長,不適合使用頂空被動加標(biāo)系統(tǒng)開展試驗,原因是DEHP的揮發(fā)性極低.實測濃度和模型預(yù)測結(jié)果偏差較小.敏感性分析表明,試驗系統(tǒng)參數(shù)中容器直徑對待測物質(zhì)在水中濃度影響較大,因此可通過調(diào)整容器直徑來縮短平衡時間,或提高平衡后化學(xué)物質(zhì)在水中的濃度.使用所構(gòu)建的模型可確定試驗系統(tǒng)平衡時間、水中待測物質(zhì)濃度和待測物質(zhì)添加體積是否滿足試驗需求,并可通過敏感性分析來優(yōu)化試驗系統(tǒng)參數(shù),提升試驗效率和質(zhì)量.以上結(jié)果完善了頂空被動加標(biāo)系統(tǒng)應(yīng)用的理論基礎(chǔ).

頂空被動加標(biāo)；疏水性物質(zhì)；動力學(xué)；多介質(zhì)模型；敏感性

在新化學(xué)物質(zhì)和農(nóng)藥等化學(xué)物質(zhì)的環(huán)境風(fēng)險評估工作中,評估數(shù)據(jù)需滿足相關(guān)性、可靠性、科學(xué)性和充分性等要求,應(yīng)優(yōu)先通過標(biāo)準(zhǔn)的試驗方法來獲取[1-4].經(jīng)濟合作發(fā)展組織(OECD)和中國生態(tài)環(huán)境部已發(fā)布了一系列試驗導(dǎo)則,如藻類生長抑制試驗[5-6]、溞類急性活動抑制試驗[7-8]和浮萍生長抑制試驗[9-10]等,來支持化學(xué)物質(zhì)的環(huán)境風(fēng)險評估.但對于疏水性高、易揮發(fā)和易降解等難處理物質(zhì),目前的水環(huán)境生態(tài)毒理學(xué)試驗方法面臨諸多挑戰(zhàn).此類物質(zhì)的水溶液制備困難,溶液濃度難以維持恒定,造成試驗結(jié)果不準(zhǔn)確,進而低估其生態(tài)毒性和環(huán)境風(fēng)險[11-12].

OECD在2019年更新了難處理物質(zhì)的水環(huán)境生態(tài)毒性試驗指南,新版指南中補充了近年來的新測試技術(shù),其中附錄6中詳細(xì)描述了適用于疏水性物質(zhì)的被動加標(biāo)法,即使用一個載有待測物質(zhì)的載體,通過待測物質(zhì)在載體和水環(huán)境中再分配,來補充試驗體系損失的量,對疏水性物質(zhì)而言比傳統(tǒng)的溶劑加標(biāo)法更準(zhǔn)確,也可降低試驗成本,減少試驗動物引入[13].但由于在試驗系統(tǒng)中引入了載體,也帶來了一些問題,如載體的清理和待測化學(xué)物質(zhì)的加載會帶來額外的工作量,載體的直接加入也可能會向溶液中引入一些小液滴,影響試驗結(jié)果.因此,相關(guān)研究開發(fā)了適用于疏水性液體物質(zhì)的頂空被動加標(biāo)系統(tǒng),在此系統(tǒng)中待測物質(zhì)先揮發(fā)充滿容器內(nèi)頂空部分,從頂空向水中再分配達(dá)到平衡[14-16].試驗操作簡便,避免了載體和水溶液直接接觸,在水環(huán)境生態(tài)毒理試驗中的應(yīng)用前景較廣,但相關(guān)研究仍處于起步階段.

多介質(zhì)模型是研究化學(xué)物質(zhì)在各介質(zhì)中遷移和轉(zhuǎn)化的常用方法,自提出以來[17],已廣泛用于實驗室和自然環(huán)境系統(tǒng)的研究中[18-23].其中IV級模型可模擬系統(tǒng)的動力學(xué)過程,也稱動力學(xué)模型[24].

為了完善頂空被動加標(biāo)系統(tǒng)應(yīng)用的理論基礎(chǔ),本研究擬使用多介質(zhì)動力學(xué)模型來構(gòu)建試驗系統(tǒng)模型,對三種性質(zhì)差異較大的物質(zhì)在試驗系統(tǒng)中濃度變化進行模擬,并使用實測濃度對模型進行驗證.同時對模型參數(shù)的敏感性進行分析,討論影響頂空被動加標(biāo)系統(tǒng)試驗結(jié)果的因素,得出試驗優(yōu)化方法,以促進頂空被動加標(biāo)方法的推廣應(yīng)用和標(biāo)準(zhǔn)化.

1 材料與方法

1.1 模型構(gòu)建

頂空被動加標(biāo)系統(tǒng)中包括三種介質(zhì),即待測物質(zhì)、空氣和水[14-16],系統(tǒng)圖示見圖1.每個介質(zhì)視為均勻混合,對每個介質(zhì)編寫微分質(zhì)量平衡方程,如下.

待測物質(zhì)相(下標(biāo)S):

空氣相(下標(biāo)A):

(3)

圖1 頂空被動加標(biāo)系統(tǒng)示意

Fig.1 Schematic diagram of headspace passive dosing system

水相:

式中各參數(shù)定義和計算方法見表1,計算方法均引自Mackay的著作[24].式(2)、(3)和(4)構(gòu)成了一個水環(huán)境頂空被動加標(biāo)系統(tǒng)的動力學(xué)模型.

1.2 試驗系統(tǒng)參數(shù)

頂空被動加標(biāo)系統(tǒng)適用的試驗導(dǎo)則及常用的試驗系統(tǒng)參數(shù)見表2,這些參數(shù)均根據(jù)試驗導(dǎo)則要求和經(jīng)驗設(shè)定.體積更大的系統(tǒng),如魚類急性毒性試驗中需要用到超過2L的試驗溶液,不便于使用頂空被動加標(biāo)方法進行配制.本文選擇表2中體積最大的浮萍生長抑制試驗系統(tǒng)和最小的魚類胚胎毒性試驗系統(tǒng)進行討論.待測物質(zhì)與空氣接觸的面積可根據(jù)需要設(shè)定,這里先按照直徑3mm的圓形計算.

1.3 物質(zhì)選擇

頂空被動加標(biāo)系統(tǒng)適用于試驗溫度下為液體且具有一定揮發(fā)性的物質(zhì).選擇常溫下為液體且水溶解度低的正壬烷、十甲基環(huán)五硅氧烷(D5)和鄰苯二甲酸二異辛酯(DEHP)進行研究,可分別代表揮發(fā)性高且降解快、揮發(fā)性低且降解慢和揮發(fā)性極低且降解快三類物質(zhì).另外,由于暫未見水溶解度低、揮發(fā)性高且降解慢的液體物質(zhì),將正壬烷的半衰期設(shè)為原值的10倍來模擬此類物質(zhì)的分配.所選擇物質(zhì)的理化性質(zhì)和半衰期見表3.

1.4 模型計算和敏感性分析

使用自由軟件R進行模型計算和敏感性分析[35],版本4.0.3.使用deSolve包(版本1.28)中ode函數(shù)求解,方法選擇lsoda.常用的敏感性分析方法有morris法和sobol法等,其中morris法僅能簡單篩選出對模型輸出影響高的參數(shù),而sobol法可對參數(shù)敏感性定量分析[36-37],本文選擇sobol法對模型參數(shù)敏感性進行分析,使用ODEsensitivity包(版本1.1.2)中ODEsobol函數(shù),模擬次數(shù)選擇1000次.

表1 參數(shù)定義及計算方法

表2 頂空被動加標(biāo)系統(tǒng)可適用的OECD生物系統(tǒng)效應(yīng)試驗導(dǎo)則及常用的試驗系統(tǒng)參數(shù)

表3 所選物質(zhì)的理化性質(zhì)和半衰期

注:1)密度數(shù)據(jù)取自Sigma-Aldrich官網(wǎng)(https://www.sigmaaldrich.com/);2)正壬烷和DEHP水溶解度數(shù)據(jù)取自文獻[32];3)除密度和水溶解度外,正壬烷的其他理化性質(zhì)取自EPI Suite 4.11[33]中試驗數(shù)據(jù);正壬烷空氣半衰期取EPI Suite 4.11中AOPWIN模塊計算值;正壬烷水中半衰期取EPI Suite 4.11中BioHCwin模塊計算值;4)除密度和水溶解度外,DEHP的理化性質(zhì)和半衰期取自文獻[24];5)除密度外,D5的理化性質(zhì)和半衰期取自文獻[34];6)*表示半衰期經(jīng)過調(diào)整.

1.5 模型驗證

Trac等[14-15]已研究了較小體積中正壬烷和D5的頂空被動加標(biāo)系統(tǒng)中質(zhì)量平衡,證明了24h已達(dá)到平衡狀態(tài).為了對模型進行進一步驗證,設(shè)計正壬烷和DEHP的浮萍生長抑制試驗系統(tǒng)及DEHP的魚類胚胎毒性試驗系統(tǒng).

所用試劑:正壬烷(純度98.5%),德國Dr. Ehrenstorfer公司;DEHP(純度99%),上海泰坦科技股份有限公司;甲醇(HPLC級),默克股份兩合公司;甲酸(HPLC級,純度398%),上海安譜實驗科技股份有限公司;正己烷(HPLC級),美國Thermo Fisher Scientific公司.

分別將100μL正壬烷和DEHP加入樣品管中,再將樣品管放入表2中的浮萍生長抑制試驗系統(tǒng)中;將20 μL DEHP加入樣品管中,再將樣品管放入表2中的魚類胚胎毒性試驗系統(tǒng)中.密封并設(shè)置多個平行樣.在24和48h后取水樣測定,由于取樣時破壞了試驗系統(tǒng)的密封性,取樣后的試驗系統(tǒng)均不再使用.

正壬烷在水中的濃度使用賽默飛Trace 1310/ISQ QD型氣質(zhì)聯(lián)用儀測定:高純氮氣作為載氣,載氣流速1mL/min;進樣口溫度260℃,分流比5:1;40℃保持2min,以15℃/min升至250℃,保持1min;EI離子源溫度為250℃,傳輸線溫度為300℃;賽默飛TG-5MS色譜柱(30m×0.25mm,0.25μm),進樣量為1μL.定量離子為43,保留時間約為5.8min.取水相樣品使用1/5體積正己烷萃取后進樣.

DEHP在水中的濃度使用安捷倫1290液相色譜儀串聯(lián)AB QTRAP 4500質(zhì)譜儀測定:流動相為99%甲醇和1%水,水中加入0.1%甲酸,流速為0.3mL/min;Waters ACQUITY UPLC BEH C18色譜柱(2.1×50mm,1.7μm),進樣量為5μL;MRM正離子模式,定量離子對為391.2/149.0,去簇電壓為50,碰撞電壓為25,離子源溫度500℃,保留時間約為3.2min.取水相樣品直接進樣.

另外,由于D5的水溶解度低于上述儀器檢出限(LOD),暫未對D5在頂空被動加標(biāo)系統(tǒng)水相中濃度進行驗證.

2 結(jié)果和討論

2.1 平衡時間和水中濃度

表4 典型試驗系統(tǒng)中平衡24h后水中濃度和完全平衡狀態(tài)水中濃度

一般來說,頂空被動加標(biāo)系統(tǒng)的平衡時間設(shè)為24h[14-15].由模型計算結(jié)果得出在浮萍生長抑制試驗系統(tǒng)和魚類胚胎毒性試驗系統(tǒng)中,平衡24h后水中濃度和完全平衡后水中濃度,結(jié)果見表4.可見,對于正壬烷和D5,平衡24h后水中濃度已達(dá)到完全平衡后的80%以上,甚至在體積較小的魚類胚胎毒性試驗系統(tǒng)中,1h的平衡時間便可令系統(tǒng)達(dá)到平衡.但對于DEHP,24h的平衡時間僅可達(dá)到完全平衡狀態(tài)下水中濃度的9%,達(dá)到平衡需要更長的時間.

對于待測物質(zhì)在水中的濃度,對比表3和表4可見,完全平衡后,并非所有物質(zhì)在水中的濃度都能接近其水溶解度.正壬烷(物質(zhì)編號1)在浮萍生長抑制試驗系統(tǒng)中,以及DEHP在兩類試驗系統(tǒng)中平衡后水中濃度與其水溶解度相差較大.體積較小的魚類胚胎毒性試驗系統(tǒng)中待測物質(zhì)在水中的濃度更高.另外,半衰期更長的物質(zhì)在水中的濃度也更接近水溶解度.

2.2 模型驗證

以10倍信噪比計算,正壬烷的定量限(LOQ)為0.087mg/L,在0.1～5mg/L濃度范圍內(nèi)有線性響應(yīng),方法回收率為73%～92%.在浮萍生長抑制試驗系統(tǒng)平衡24和48h后,分別測得水中濃度為0.040和0.038mg/L,表明24h已可達(dá)到平衡狀態(tài).表4的計算結(jié)果中平衡24和48h的濃度分別為0.068和0.072mg/L,較實測值稍高.

以3倍信噪比計算,DEHP的LOD為0.0012mg/L.在使用頂空被動加標(biāo)方法暴露24和48h后將水相直接進樣,測得浮萍生長抑制試驗系統(tǒng)和魚類胚胎毒性試驗系統(tǒng)中DEHP的濃度均

從以上結(jié)果,結(jié)合Trac等[14-15]開展的小體積頂空被動加標(biāo)系統(tǒng)測試結(jié)果,可見模型可較好地模擬頂空被動加標(biāo)系統(tǒng)中水相濃度的變化.引起模擬結(jié)果和實測結(jié)果存在偏差的原因有:1)生態(tài)毒理試驗系統(tǒng)的水相含有營養(yǎng)鹽等組分,比水溶解度測定中使用的蒸餾水組分復(fù)雜,因此疏水性有機物濃度總是低于其水溶解度[13],而模型中使用的參數(shù)為蒸餾水中有機物水溶解度數(shù)據(jù).2)疏水性液體有機物的水溶解度數(shù)據(jù)有一定不確定性.例如,不同研究者得出的正壬烷水溶解度數(shù)據(jù)相差超過1倍[32,38].經(jīng)典的水溶解度測試方法中,有機物和水直接接觸,其在水中乳化可能會引起水中溶解度實測值高于實際溶解的濃度值,測試方法有待進一步研究[38].水溶解度對模型計算結(jié)果影響很大(見2.3節(jié)),獲取準(zhǔn)確的水溶解度數(shù)據(jù)可提高模型計算結(jié)果準(zhǔn)確性.3)檢測方法回收率很難達(dá)到100%,引起實測濃度低于真實值.4)模型中半衰期的不確定性較大[39],也會影響模型計算結(jié)果的準(zhǔn)確性.

2.3 敏感性分析

表5 平衡24h時水中待測物質(zhì)濃度的模型參數(shù)總階敏感性指數(shù)

如前文所述,不同的試驗系統(tǒng)參數(shù)和待測物質(zhì)性質(zhì)會影響水中待測物質(zhì)濃度和平衡時間,為了進一步研究其影響因素,對試驗系統(tǒng)和待測物質(zhì)性質(zhì)共14個參數(shù)進行敏感性分析,每個參數(shù)浮動±10%,且在此區(qū)間設(shè)為均勻分布,得出其敏感性指數(shù), 見表5.由于一階敏感性和總階敏感性指數(shù)差別很小,說明每個參數(shù)與其他參數(shù)的交互作用很弱[40],因此表5中僅列出了總階敏感性指數(shù).

從表4和表5可見,在平衡24h后,對于已達(dá)到完全平衡狀態(tài)的情況,即編號1、2和3的物質(zhì)在魚類胚胎毒性試驗系統(tǒng),敏感性指數(shù)最高的參數(shù)為25 ℃下水溶解度,其次為試驗溫度,其他參數(shù)的敏感性指數(shù)很低;對于接近平衡狀態(tài)但未達(dá)到完全平衡狀態(tài)的情況,即編號1、2和3的物質(zhì)在浮萍生長抑制試驗系統(tǒng),敏感性指數(shù)最高的參數(shù)為25 ℃下水溶解度,其次為容器直徑,其他參數(shù)的敏感性指數(shù)均較低;對于遠(yuǎn)未達(dá)到平衡狀態(tài)的情況,即物質(zhì)4在兩類系統(tǒng)中,敏感性指數(shù)最高的參數(shù)為試驗溫度和容器直徑,其他參數(shù)的敏感性指數(shù)較低.

2.4 模型應(yīng)用

2.4.1平衡時間 通過敏感性分析的結(jié)果,可優(yōu)化驗系統(tǒng)參數(shù),縮短平衡時間.較敏感的參數(shù)中,25 ℃下水溶解度是待測物質(zhì)的理化參數(shù),不能進行調(diào)整;而試驗溫度的調(diào)整范圍很小.因此對于平衡時間較長的物質(zhì),可通過改變試驗系統(tǒng)參數(shù)中容器直徑、水相高度、空氣相高度和待測物質(zhì)與空氣接觸的面積等方式進行嘗試.

從2.2節(jié)討論可知,系統(tǒng)未達(dá)到平衡時,試驗系統(tǒng)參數(shù)中容器直徑是最敏感的參數(shù),可優(yōu)先對其進行優(yōu)化.如對于本研究中使用的正壬烷(物質(zhì)編號2),在適當(dāng)減少浮萍生長抑制試驗系統(tǒng)容器直徑后,可縮短平衡時間:當(dāng)容器直徑從9cm減小到6cm,達(dá)到平衡濃度的80%以上的時間便可從18h縮短至11h.但對于DEHP,將魚類胚胎毒性試驗系統(tǒng)的直徑和水相高度減至1cm,將空氣相高度減至3mm,且將待測物質(zhì)與空氣接觸的直徑增加至5mm(這已是開展此試驗的極限條件),此時達(dá)到平衡濃度80%的時間僅從387h減至332h,這個時間仍然較長.過長的平衡時間不僅降低了試驗效率,也會引起系統(tǒng)內(nèi)微生物增長,同時增加了待測物質(zhì)在水中由水解等途徑轉(zhuǎn)化的可能性,從而影響試驗結(jié)果.因此對于DEHP這類揮發(fā)性極低的物質(zhì),不建議使用頂空被動加標(biāo)法開展試驗.

雖然待測物質(zhì)的揮發(fā)性決定了是否可通過優(yōu)化頂空被動加標(biāo)試驗系統(tǒng)參數(shù)來縮短平衡時間,但由于平衡時間還由半衰期等參數(shù)綜合決定,所以并不能計算得出頂空被動加標(biāo)系統(tǒng)適用的物質(zhì)蒸汽壓或亨利常數(shù)范圍.試驗前可通過1.1節(jié)構(gòu)建的模型來計算平衡時間,確定頂空被動加標(biāo)系統(tǒng)對待測物質(zhì)是否適用.

2.4.2 水中濃度 從表4可見,完全平衡后,正壬烷(物質(zhì)編號1)在體積較大的浮萍生長抑制試驗系統(tǒng)水中的濃度(0.072mg/L)并未達(dá)到其水溶解度(0.122mg/L)的80%以上,按此條件開展試驗不能準(zhǔn)確評估其生態(tài)毒性,因此需探討如何提升這種情況下水中待測物質(zhì)的濃度.

從表5的敏感性分析結(jié)果可見,在浮萍生長抑制試驗系統(tǒng)中,試驗系統(tǒng)參數(shù)中容器直徑的敏感性最高,對水中濃度影響最大.因此可調(diào)整容器直徑,提升平衡后待測物質(zhì)在水中的濃度.如將容器直徑從9cm減至5.2cm以下,便可將正壬烷(物質(zhì)編號1)在平衡24h后的水中濃度提升至其水溶解度的80%以上.

雖然Trac等[14-15]的初步研究認(rèn)為頂空被動加標(biāo)系統(tǒng)可獲取濃度接近待測物質(zhì)水溶解度的水溶液,但其研究結(jié)論是基于較小的試驗系統(tǒng)下得出的,適用范圍有限.本研究得出了試驗系統(tǒng)參數(shù)優(yōu)化方法,拓展了頂空被動加標(biāo)法的應(yīng)用范圍.

2.4.3 待測物質(zhì)添加體積 在試驗前可通過1.1節(jié)構(gòu)建的模型計算結(jié)果,來確定頂空被動加標(biāo)系統(tǒng)中待測物質(zhì)需要加入的體積.計算得到的S應(yīng)始終為正值,否則應(yīng)增大S的初始值.雖然試驗期間也可觀察補充待測物質(zhì),但此操作風(fēng)險較大,有可能破壞試驗系統(tǒng)密封性和平衡狀態(tài).應(yīng)盡量在試驗開始時加入足量的待測物質(zhì).

3 結(jié)論

3.1 構(gòu)建了頂空被動加標(biāo)系統(tǒng)的多介質(zhì)動力學(xué)模型,根據(jù)模型計算結(jié)果,正壬烷和D5在浮萍生長抑制試驗系統(tǒng)和魚類胚胎毒性試驗系統(tǒng)中很快達(dá)到平衡,平衡24h已達(dá)到完全平衡時濃度的80%以上;但DEHP的平衡時間很長,不適合使用頂空被動加標(biāo)系統(tǒng)開展試驗.另外,并非所有物質(zhì)在頂空被動加標(biāo)系統(tǒng)中水中的濃度都能接近其水溶解度.使用實測濃度對模型進行驗證,表明模型可較好地模擬實際系統(tǒng).

3.2 模型敏感性分析結(jié)果表明,水溶解度、容器直徑和試驗溫度對未達(dá)到平衡時水中濃度影響較大,而達(dá)到平衡后水中濃度主要受水溶解度影響.

3.3 模型可在三個方面得到應(yīng)用,即確定待測物質(zhì)添加體積、系統(tǒng)平衡時間和水中待測物質(zhì)濃度是否滿足試驗需求.通過敏感性分析結(jié)果,可對試驗系統(tǒng)進行優(yōu)化,提升試驗效率和質(zhì)量.

[1] 中華人民共和國生態(tài)環(huán)境部.新化學(xué)物質(zhì)環(huán)境管理登記指南[EB/OL]. http://www.mee.gov.cn/xxgk2018/xxgk/xxgk01/202011/ W020201119648157169625.pdf.

Ministry of Ecology and Environment of the People’s Republic of China. Guidelines on environmental management registration for new chemical substances [EB/OL]. http://www.mee.gov.cn/xxgk2018/ xxgk/xxgk01/202011/W020201119648157169625.pdf.

[2] 中華人民共和國生態(tài)環(huán)境部.化學(xué)物質(zhì)環(huán)境與健康危害評估技術(shù)導(dǎo)則(試行) [EB/OL]. http://www.mee.gov.cn/xxgk2018/xxgk/ xxgk01/ 202012/W020201225515743628373.pdf.

Ministry of Ecology and Environment of the People’s Republic of China. Technical guidelines on environmental and health hazard assessment of chemical substances (Trial) [EB/OL]. http://www.mee. gov.cn/xxgk2018/xxgk/xxgk01/202012/W020201225515743628373.pdf.

[3] NY/T 2882.2-2016 農(nóng)藥登記環(huán)境風(fēng)險評估指南第2部分:水生生態(tài)系統(tǒng)[S].

NY/T 2882.2-2016 Guidance on environmental risk assessment for pesticide registratioin-part 2: aquatic ecosystem [S].

[4] GB/T 34033.1-2017/ISO 13073-1:2012船舶與海上技術(shù)船舶防污底系統(tǒng)風(fēng)險評估第1部分:船舶防污底系統(tǒng)用防污活性物質(zhì)的海洋環(huán)境風(fēng)險評估方法[S].

GB/T 34033.1-2017/ISO 13073-1:2012 Ships and marine technology-Risk assessment on anti-fouling systems on ships-part1: marine environmental risk assessment method of biocidally active substances used for anti-fouling systems on ships [S].

[5] Test No. 201: Freshwater alga and cyanobacteria, growth inhibition Test [S].

[6] 環(huán)境保護部化學(xué)品登記中心.化學(xué)品測試方法生物系統(tǒng)效應(yīng)卷[M]. 北京:中國環(huán)境出版社, 2013:6-21.

Chemical Registration Center of Ministry of Environmental Protection. The guidelines for the testing of chemicals: effects on biotic systems [M]. Beijing: China Environmental Science Press, 2013:6-21.

[7] Test NO.202:sp. Acute Immobilization Test [S].

[8] 環(huán)境保護部化學(xué)品登記中心.化學(xué)品測試方法生物系統(tǒng)效應(yīng)卷[M]. 北京:中國環(huán)境出版社, 2013:22-29.

Chemical Registration Center of Ministry of Environmental Protection. The guidelines for the testing of chemicals: effects on biotic systems [M]. Beijing: China Environmental Science Press, 2013:22-29.

[9] Test No. 221:sp.Growth Inhibition Test [S].

[10] 環(huán)境保護部化學(xué)品登記中心.化學(xué)品測試方法生物系統(tǒng)效應(yīng)卷[M]. 北京:中國環(huán)境出版社, 2013:206-217.

Chemical Registration Center of Ministry of Environmental Protection. The guidelines for the testing of chemicals: effects on biotic systems [M]. Beijing: China Environmental Science Press, 2013:206-217.

[11] Hodson P V, Adams J, Brown R S. Oil toxicity test methods must be improved [J]. Environmental Toxicology and Chemistry, 2019, 38(2):302-311.

[12] Birch H, Redman A D, Letinski D J, et al. Determining the water solubility of difficult-to-test substances: a tutorial review [J]. Analytica Chimica Acta, 2019,1086:16-28.

[13] Series on testing and assessment No. 23 (second edition) Guidance document on aquatic toxicity testing of difficult substances and mixtures [S].

[14] Trac L N, Schmidt S N, Mayer P. Headspace passive dosing of volatile hydrophobic chemicals – Aquatic toxicity testing exactly at the saturation level [J]. Chemosphere, 2018,211:694-700.

[15] Trac L N, Schmidt S N, Holmstrup M, et al. Headspace passive dosing of volatile hydrophobic organic chemicals from a lipid donor—linking their toxicity to well-defined exposure for an improved risk assessment [J]. Environmental Science & Technology, 2019,53(22): 13468-13476.

[16] Trac L N, Sjo?holm K K, Birch H, et al. Passive dosing of petroleum and essential oil UVCBs—whole mixture toxicity testing at controlled exposure [J]. Environmental Science & Technology, 2021,55(9): 6150-6159.

[17] Mackay D. Finding fugacity feasible [J]. Environmental Science & Technology, 1979,13(10):1218-1223.

[18] 莫俊超,舒耀皋,董 晶,等.預(yù)測沉積物-水微宇宙系統(tǒng)中化學(xué)品濃度變化的多介質(zhì)模型[J]. 生態(tài)毒理學(xué)報, 2015,10(6):101-109.

Mo J, Shu Y, Dong J, et al. A multimedia model for prediction of chemicals concentration changes in sediment-water microcosm systems [J]. Asian Journal of Ecotoxicology, 2015,10(6):101-109.

[19] O’Driscoll K, Robinson J, Chiang W S, et al. The environmental fate of polybrominated diphenyl ethers (PBDEs) in western Taiwan and coastal waters: evaluation with a fugacity-based model [J]. Environmental Science and Pollution Research, 2016,23(13):13222- 13234.

[20] Nurmi T M A, Kiljunen T K, Knuutinen J S. A fugacity model assessment of ibuprofen, diclofenac, carbamazepine, and their transformation product concentrations in an aquatic environment [J]. Environmental Science and Pollution Research, 2019,26(1):328-341.

[21] 程浩淼,陳玉茹,趙永嶺,等.巢湖水域四溴雙酚A的多介質(zhì)遷移與歸趨模擬[J]. 中國環(huán)境科學(xué), 2019,39(1):314-320.

Cheng H, Chen Y, Zhao Y, et al. Simulation of multimedia transfer and fate of tetrabromobisphenol A in Lake Chaohu [J]. China Environmental Sciencece, 2019,39(1):314-320.

[22] Chen Y, Liu X, Dong F, et al. Characterization of the fate and distribution of methoxyfenozide in a water-plant-fish-sediment microcosm using a multimedia fugacity model [J]. Science of the Total Environment, 2021,755:142482.

[23] 毛瀟萱,宋世杰,李子璇,等.我國半干旱區(qū)工業(yè)化城鎮(zhèn)PCBs干沉降及歸趨[J]. 中國環(huán)境科學(xué), 2021,41(4):1521-1529.

Mao X, Song S, Li Z, et al. Dry deposition and fate of PCBs in a semi-arid and industrialized city, China [J]. China Environmental Science, 2021,41(4):1521-1529.

[24] Mackay D. Multimedia environmental model: the fugacity approach [M]. BOCA Raton: Lewis Publisher, 2001.

[25] 中華人民共和國生態(tài)環(huán)境部.化學(xué)物質(zhì)環(huán)境與健康暴露評估技術(shù)導(dǎo)則(試行) [EB/OL]. http://www.mee.gov.cn/xxgk2018/xxgk/ xxgk01/ 202012/W020201225515745341837.pdf.

Ministry of Ecology and Environment of the People’s Republic of China. Technical guidelines on environmental and health exposure assessment of chemical substances (Trial) [EB/OL]. http://www. mee.gov.cn/xxgk2018/xxgk/xxgk01/202012/W020201225515745341837.pdf.

[26] Mackay D, Guardo A D, Paterson S, et al. Evaluating the environmental fate of a variety of types of chemicals using the EQC model [J]. Environmental Toxicology and Chemistry, 1996,15(9): 1627-1637.

[27] Test No. 210: Fish, Early-life Stage Toxicity Test [S].

[28] Test No. 211:Reproduction Test [S].

[29] Test No. 212: Fish, Short-term Toxicity Test on Embryo and Sac-Fry Stages [S].

[30] Test No. 235:sp., Acute Immobilisation Test [S].

[31] Test No. 236: Fish Embryo Acute Toxicity (FET) Test [S].

[32] Yalkowsky S H, Yan H, Jain P. Handbook of aqueous solubility data [M]. Boca Raton.: CRC Press, 2010:632,1217.

[33] EPA U S. EPI Suite: United States [Z/OL]. https://www.epa.gov/ tsca-screening-tools/download-epi-suitetm-estimation-program-interface-v411.

[34] Hughes L, Mackay D, Powell D E, et al. An updated state of the science EQC model for evaluating chemical fate in the environment: application to D5 (decamethylcyclopentasiloxane) [J]. Chemosphere, 2012,87(2):118-124.

[35] 莫俊超,吳孝槐,舒耀皋.Julia和R軟件在多介質(zhì)動力學(xué)模型建模中的比較 [J]. 化工環(huán)保, 2021,41(4):494-499.

Mo J, Wu X, Shu Y. Comparison of software Julia and R in modeling of multimedia dynamic models [J]. Environmental Protection of Chemical Industry, 2021,41(4):494-499.

[36] 陳衛(wèi)平,涂宏志,彭 馳,等.環(huán)境模型中敏感性分析方法評述[J]. 環(huán)境科學(xué), 2017,38(11):4889-4896.

Chen W, Tu H, Peng C, et al. Comment on sensitivity analysis methods for environmental models [J]. Environmental Science, 2017,38(11):4889-4896.

[37] 歐陽琦,盧文喜,侯澤宇,等.基于替代模型的地下水溶質(zhì)運移不確定性分析[J]. 中國環(huán)境科學(xué), 2016,36(4):1119-1124.

Ouyang Q, Lu W, Hou Z, et al. Uncertainty analysis of groundwater solute transport based on surrogate model [J]. China Environmental Science, 2016,36(4):1119-1124.

[38] Letinski D J, Parkerton T F, Redman A D, et al. Water solubility of selected C9-C18 alkanes using a slow-stir technique: comparison to structure-property models [J]. Chemosphere, 2016,150:416-423.

[39] Citra M J. Incorporating Monte Carlo analysis into multimedia environmental fate models [J]. Environmental Toxicology and Chemistry, 2004,23(7):1629-1633.

[40] 史鑫蕊,梁 浩,周 豐,等.稻田土壤-作物系統(tǒng)模型參數(shù)敏感性分析與模型驗證[J]. 農(nóng)業(yè)機械學(xué)報, 2020,51(5):252-262.

Shi X, Liang H, Zhou F, et al. Sensitivity Analysis and Parameter Estimation for Soil-Rice System Model [J]. Transactions of the Chinese Society for Agricultural Machinery, 2020,51(5):252-262.

Multimedia dynamic model for headspace passive dosing system.

MO Jun-chao1,2, YAO Hong-wei1,2, WU Xiao-huai3, GUO Dan-dan1,2, ZHAO Fang-zhou1,2, SHU Yao-gao1,2*

(1.Shanghai Research Institute of Chemical Industry Co.,Ltd., Shanghai 200062, China；2.Shanghai Engineering Research Center of Chemicals Public Safety, Shanghai 200062, China；3.China National Accreditation Service for Conformity Assessment, Beijing 100062, China)., 2021,41(11)：5353～5360

Adynamic model for headspace passive dosing system (HPDS) on hydrophobic liquid substances was constructed using multimedia modeling method. Taking the lemna growth inhibition testing system and the fish embryo toxicity testing system as examples, the concentration changes of three hydrophobic liquid substances-nonane, decamethylcyclopentasiloxane (D5) and di (2-ethylhexyl) phthalate (DEHP) with different properties in the HPDS were obtained. The model was verified by the measured concentrations. The results revealed that after 24h equilibrium of-nonane and D5 in two systems, the aqueous concentrations reached more than 80% of the concentrations at complete equilibrium. However, considering the long equilibrium time of DEHP, HPDS was not applicable for performing tests because of the extremely low volatility of DEHP. The deviations between the measured concentrations and the results predicted by the model were small. Sensitivity analysis indicated that the diameter of the container in the testing system parameters had a great impact on the aqueous concentration of test substance. Therefore, it is feasible to adjust the diameter of the container for shortening equilibrium time, or increasing the aqueous concentrations of the chemicals after the equilibrium. The model established in this paper can be used to determine whether the equilibrium time of the testing system, the aqueous concentrations of the test substance and the added volume of the test substance can meet the testing requirements. In addition, the parameters of the testing system can be optimized through sensitivity analysis, thereby enhancing the testing efficiency and quality. The above results are expected to strengthen the theoretical basis for the application of HPDS.

headspace passive dosing；hydrophobic substance；dynamic；multimedia model；sensitivity

X502;X592

A

1000-6923(2021)11-5353-08

莫俊超(1986-),男,陜西咸陽人,高級工程師,碩士,主要從事化學(xué)物質(zhì)環(huán)境安全評價和環(huán)境系統(tǒng)模型研究.發(fā)表論文7篇.

2021-03-31

上海市科委上?；瘜W(xué)品公共安全工程技術(shù)研究中心項目(18DZ2280700);上海化工研究院科研計劃項目(N307-2016jczx-01)

* 責(zé)任作者, 高級工程師, syg@ghs.cn