鄧桂榮，梅承芳，柳燕貞，張宏濤,4，陳澤芳，曾國驅(qū),*

1. 廣東省微生物分析檢測中心 生態(tài)毒理與環(huán)境安全實驗室，廣州 510070 2. 廣東省微生物研究所 廣東省菌種保藏與應用重點實驗室，廣州 510070 3. 省部共建華南應用微生物國家重點實驗室，廣州 510070 4. 佛山市環(huán)境健康與安全評價研究中心，佛山 528000

穩(wěn)健統(tǒng)計在快速生物降解性試驗質(zhì)量分析中的應用

鄧桂榮1,2,3，梅承芳1,2,3，柳燕貞1,2,3，張宏濤1,2,3,4，陳澤芳1,2,3，曾國驅(qū)1,2,3,*

1. 廣東省微生物分析檢測中心 生態(tài)毒理與環(huán)境安全實驗室，廣州 510070 2. 廣東省微生物研究所 廣東省菌種保藏與應用重點實驗室，廣州 510070 3. 省部共建華南應用微生物國家重點實驗室，廣州 510070 4. 佛山市環(huán)境健康與安全評價研究中心，佛山 528000

快速生物降解性試驗在新化學物質(zhì)的生態(tài)毒理測試中占據(jù)重要地位，其測試數(shù)據(jù)的可靠性直接影響后續(xù)的測試策略和物質(zhì)的申報登記、以及新化學物質(zhì)進入環(huán)境后的風險管理。為研究OECD 301系列試驗(301B、301D和301F)的質(zhì)量，選用本機構近5年來積累的標準化物質(zhì)苯甲酸鈉14 d快速生物降解率數(shù)據(jù)，首次采用穩(wěn)健統(tǒng)計技術，以穩(wěn)健總計統(tǒng)計量的Z比分數(shù)作為評價指標，對各試驗數(shù)據(jù)和試驗系統(tǒng)的質(zhì)量開展評估。結果顯示，當301B、301D和301F試驗的苯甲酸鈉降解率分別在63.8%～88.0%、63.4%～89.2%和66.0%～99.8%范圍內(nèi)時，表明測試結果滿意，該批次試驗數(shù)據(jù)質(zhì)量良好；若苯甲酸鈉降解率不在此范圍，則表明測試結果可能存在問題，甚至得到的可能是錯誤數(shù)據(jù)，需要進一步驗證。據(jù)此提出了對實際測試工作具有指導意義的質(zhì)量分析對照表，以期為測試機構開展新化學物質(zhì)快速生物降解性試驗時的質(zhì)量分析提供參考。

苯甲酸鈉；快速生物降解性試驗；穩(wěn)健統(tǒng)計；Z比分數(shù)；質(zhì)量分析

化合物快速生物降解性是有機物作為唯一碳源在限定時間內(nèi)與微生物接觸表現(xiàn)出的生物可降解能力[1]，是鑒別化學物質(zhì)的環(huán)境危害性和持久性、進行分類和標簽[2]、評估和控制環(huán)境風險[3]的基本指標，也是我國實施新化學物質(zhì)申報登記制度，對其進行環(huán)境風險管理的重要基礎[4-5]。因此快速生物降解性試驗數(shù)據(jù)的科學性、可靠性意義重大。而能基本確保試驗數(shù)據(jù)經(jīng)得起推敲和質(zhì)疑的，是與各快速生物降解性試驗配套的一系列質(zhì)量控制指標，比如14 d程序?qū)φ战M和毒性對照組的生物降解率等。

在經(jīng)濟合作與發(fā)展組織(OECD)的6種經(jīng)典快速生物降解性測試方法[6]中，判定試驗系統(tǒng)是否有效的重要觀測點是陽性參比物14 d的生物降解率是否達到了“通過水平”，即生物降解率大于60%或70%。然而，即使參比物的降解率達到方法的最低要求，試驗數(shù)據(jù)依然具有較大的可變性，嚴格來說并不能因此就斷定所得數(shù)據(jù)是科學、可靠的，有必要開展進一步的質(zhì)量分析，評估該數(shù)據(jù)的可接受程度，進而對試驗系統(tǒng)的運行狀況做出有效判斷。近年來我國環(huán)境壓力不斷加大，新化學物質(zhì)的管理也日趨嚴格，但主管部門的關注點僅限于所提供數(shù)據(jù)之于測試標準的遵循性和符合性，對數(shù)據(jù)質(zhì)量未做額外要求。作為提供新化學物質(zhì)測試數(shù)據(jù)且符合良好實驗室管理規(guī)范(GLP)的測試機構，采用科學的統(tǒng)計方法對即將用于申報的快速生物降解性試驗數(shù)據(jù)開展進一步的質(zhì)量評估具有重要的科學意義和應用價值。目前，國內(nèi)尚無關于快速生物降解性數(shù)據(jù)質(zhì)量分析的相關研究和報道。

本文利用近5年來采集的OECD快速生物降解性301系列試驗[7](301B：CO2產(chǎn)生試驗[8]，301D：密閉瓶法試驗[9]和301F：測壓呼吸計量法試驗[10])中苯甲酸鈉的128個標準化的快速生物降解率數(shù)據(jù)，采用穩(wěn)健統(tǒng)計分析技術，以基于穩(wěn)健總計統(tǒng)計量的Z比分數(shù)作為表征和評價指標，評估試驗數(shù)據(jù)和試驗系統(tǒng)的質(zhì)量，并給出對實際測試工作具有指導意義的質(zhì)量分析對照表，以期為測試機構開展新化學物質(zhì)快速生物降解性試驗時進行質(zhì)量分析提供參考。

1 材料與方法(Materials and methods)

1.1 試劑

苯甲酸鈉(≥99.5%)、磷酸二氫鉀(≥99.0%)、三水合磷酸氫二鉀(≥99.0%)、十二水合磷酸氫二鈉(≥99.0%)、氯化銨(≥99.5%)、無水氯化鈣(≥96.0%)、七水合硫酸鎂(≥99.0%)、六水合氯化鐵(≥99.0%)。

以上試劑中，苯甲酸鈉作為標準推薦使用的參比物之一，加入到試驗系統(tǒng)，其余用于配制試驗培養(yǎng)基。所有試劑均為分析純，購自廣州化學試劑廠。

1.2 儀器設備

301B：3 L燒瓶(2 L反應液)，數(shù)顯滴定儀(Titrette，德國，BRAND GMBH + CO KG)，脫CO2裝置，空氣流量計；

301D：300 mL BOD瓶(充滿反應液)，溶解氧測定儀(Multi 3430，德國，WTW)；

301F：500 mL試驗瓶(一定體積反應液，留有上層頂空)，測試系統(tǒng)(OxiTop?Control，德國，WTW)，生化培養(yǎng)箱(SHP-350，中國，上海精宏實驗設備有限公司)。

1.3 接種物的采集與準備

試驗所用的接種物取自廣州地區(qū)2家主要處理生活污水且運行狀況良好的污水處理廠。

根據(jù)測試標準推薦，301B和301F試驗取用曝氣池末端的活性污泥，采集后去除粗糙顆粒和浮在表面上的雜質(zhì)，再用試驗培養(yǎng)基離心清洗，棄除上清液，并將此接種物重新懸浮于試驗培養(yǎng)基中，在均勻狀態(tài)下取少量懸浮液測定混合液懸浮固體濃度(MLSS)[11]。試驗前，使用試驗培養(yǎng)基將其濃縮或稀釋，得到活性污泥濃度為4.0 g·L-1(以干重計)的接種物。保持有氧狀態(tài)直至使用。

301D試驗取用二級出水，不含有菌膠團。使用前，采用營養(yǎng)瓊脂培養(yǎng)基測定其可培養(yǎng)細胞總數(shù)[12]，達到104～106cell·mL-1培養(yǎng)基，方能滿足測試標準的要求。

1.4 試驗條件

所開展的盡管都是快速生物降解性試驗，但各個類型試驗的條件存在一定差異[13]，具體見表1。

1.5 方法

按照我國《化學品測試方法》[1,14]和OECD化學品測試導則[7]中的標準條件進行試驗，所有試驗均遵循GLP規(guī)范的要求[15-16]，操作步驟按照標準操作程序(SOP)的規(guī)定實施。

1.5.1 301B試驗

在已接種的無機培養(yǎng)基中加入適量的苯甲酸鈉作為唯一有機碳源。在黑暗條件下，用脫CO2的空氣以受控速率對試驗培養(yǎng)基進行曝氣。通過測定CO2的產(chǎn)生量來確定苯甲酸鈉的降解率。CO2用氫氧化鋇溶液吸收，使用Titrette數(shù)顯滴定儀滴定剩余的氫氧化鋇來測定CO2產(chǎn)生量。苯甲酸鈉產(chǎn)生的CO2數(shù)量用理論二氧化碳(ThCO2)的百分率來表示。

1.5.2 301D試驗

將苯甲酸鈉溶于試驗培養(yǎng)基，接種適量二級出水后，分裝至BOD瓶中，使BOD瓶完全裝滿，輕輕拍打BOD瓶以去除任何氣泡。密閉BOD瓶后置于試驗條件下培養(yǎng)，并在試驗期間使用Multi 3430溶解氧測定儀追蹤測定溶解氧消耗[17]來計算生物降解率。

1.5.3 301F試驗

向一定體積已接種的無機培養(yǎng)基中加入已知濃度的苯甲酸鈉貯備液作為唯一有機碳源，密閉試驗瓶后在試驗條件下連續(xù)攪拌培養(yǎng)。通過OxiTop?Control BOD測試系統(tǒng)測定瓶內(nèi)壓力的變化，可得到耗氧量。釋放的CO2用NaOH顆粒吸收。苯甲酸鈉生物降解過程中的耗氧量用空白對照組校正，微生物種群吸收氧的量以理論需氧量(ThOD)的百分率來表示。

1.6 數(shù)據(jù)統(tǒng)計

1.6.1 301B試驗降解率的計算

本試驗以氫氧化鋇溶液作吸收劑，根據(jù)HCl滴定體積，計算空白對照組和苯甲酸鈉試驗組所產(chǎn)生的CO2量，二者之差即為苯甲酸鈉生物降解產(chǎn)生的CO2量。將CO2量除以加入試驗體系的苯甲酸鈉的TOC含量與以碳轉(zhuǎn)換CO2系數(shù)3.67的乘積，即為苯甲酸鈉降解百分率[8]。

1.6.2 301D和301F試驗降解率的計算

用第14天的苯甲酸鈉試驗組耗氧量減去相應的空白對照組耗氧量，再除以苯甲酸鈉濃度，得到用毫克氧比毫克苯甲酸鈉表示的BOD，結合ThOD可計算301D試驗中苯甲酸鈉降解百分率[9]。

301F試驗中苯甲酸鈉降解百分率計算方法類似于301D試驗，只需用測得的BOD替代耗氧量計算即可[10]。

1.6.3 降解百分率的穩(wěn)健統(tǒng)計

穩(wěn)健統(tǒng)計前，應用SPSS Statistics 16.0軟件的單個樣本的K-S檢驗(One-Sample Kolmogorov-Smirnov Test)，對各試驗的苯甲酸鈉降解率數(shù)據(jù)進行正態(tài)分布檢驗。

針對2010年至2014年在本機構開展的128個標準化的苯甲酸鈉14 d快速生物降解性結果，使用通用的穩(wěn)健統(tǒng)計方法[18-22]處理，計算穩(wěn)健Z比分數(shù)，分析各數(shù)據(jù)在總體數(shù)據(jù)中所處的位置。

計算公式：

Z=(X-Xm)/標準化IQR(X)

式中：

Z—穩(wěn)健比分數(shù)；X—降解率結果；Xm—各類型試驗所有降解率數(shù)據(jù)的中位值；標準化IQR—即IQR與因子0.7413的乘積。

表1 快速生物降解性試驗條件Table 1 Test conditions of ready biodegradability tests

所涉及的統(tǒng)計量還有四分位數(shù)值、變異系數(shù)、最小值、最大值和極差等。

應用Microsoft Office Excel 2003軟件[23]作穩(wěn)健Z比分數(shù)序列圖，圖中按照大小的順序顯示出各試驗數(shù)據(jù)的穩(wěn)健Z比分數(shù)，并標有相應的試驗編號。

將各降解率結果換算為概率單位[24]后，使用origin 8軟件以概率單位對試驗編號作圖。

2 結果與分析(Results and analysis)

表2給出了各類型試驗數(shù)據(jù)統(tǒng)計結果。

K-S檢驗結果表明，301B、301D和301F試驗的苯甲酸鈉降解率數(shù)據(jù)的雙尾漸進顯著性概率分別為0.723、0.705和0.830，均大于0.05，因此所得數(shù)據(jù)服從正態(tài)分布，適用基于穩(wěn)健總計統(tǒng)計量的Z比分數(shù)統(tǒng)計。

評價數(shù)據(jù)時，采用穩(wěn)健Z比分數(shù)作為評價指標，該值受數(shù)據(jù)中離群值影響較小。Z比分數(shù)的“+”與“-”分別代表了某個測試值偏離中位值的方向，也反映該測試值的系統(tǒng)誤差，同時還反映測試值與穩(wěn)健值的偏離程度，Z比分數(shù)越接近X軸，表示當次苯甲酸鈉降解率結果偏離程度越小，數(shù)據(jù)的可接受程度越高，相反，偏離程度越大，相應批次試驗所得數(shù)據(jù)越趨于異常。依據(jù)評價原則[18,25]，當|Z|≤2時表示測試結果滿意，該批次試驗數(shù)據(jù)質(zhì)量良好；當2<|Z|<3時表示結果可能存在問題，項目負責人(SD)應引起重視，并進行復查，尋找出現(xiàn)較大偏差的原因，以進一步提高數(shù)據(jù)質(zhì)量；當|Z|≥3時表示結果不滿意，屬于統(tǒng)計分析中的離群值，SD應慎密檢查，找出問題所在，并制定有效的糾正措施，糾正偏離。如有必要，應重新啟動試驗加以驗證。

表2 苯甲酸鈉降解率數(shù)據(jù)統(tǒng)計結果Table 2 Statistical results of sodium benzoate biodegradation data

注：*分別取最小降解率和最大降解率，以及與其相鄰的一個值。

Note: The min, max value and its adjacent were obtained respectively.

由數(shù)據(jù)統(tǒng)計結果可知，301B試驗中，降解率的最值處均出現(xiàn)了可疑值，其穩(wěn)健Z比分數(shù)分別為-2.400、-2.065和2.007、2.836，超出了|Z|≤2范圍，對應的降解率分別為60.0%、62.3%和90.3%、96.0%，表明當降解率≤62.3%或≥90.3%時，結果可能存在問題。其余結果均為滿意，數(shù)據(jù)質(zhì)量良好。301D試驗中，可疑值僅出現(xiàn)在最大降解率和次最大降解率處，其穩(wěn)健Z比分數(shù)分別為2.463和2.044，對應的降解率為95.0%和98.3%，即當降解率≥95.0%時，試驗結果可能存在問題甚至離群，應復查或驗證。對于301F試驗，降解率62.8%～66.0%是存在問題的結果，其對應的穩(wěn)健Z比分數(shù)分別為-2.345～-1.977，因此應關注并分析降解率≤66.0%的結果。另外，盡管301D試驗最小降解率對應的Z比分數(shù)為-1.968，301F試驗最大降解率對應的Z比分數(shù)為1.920，其絕對值均小于2.0，未指示結果異常，但在具體的實踐過程中，若遇到該接近臨界的數(shù)據(jù)，即當301D試驗降解率小于63.4%和301F試驗降解率大于99.8%時，SD也應注意，慎重分析，因為就目前的數(shù)據(jù)仍無法判斷該數(shù)據(jù)的質(zhì)量水平。3類試驗的降解率對應的穩(wěn)健Z比分數(shù)均未出現(xiàn)|Z|≥3，即苯甲酸鈉降解率合格程度達100%，無不滿意或離群的情況。

表3為快速生物降解性試驗質(zhì)量分析對照表，是對上述結果的總結。

圖1為概率單位趨勢圖，從中可以看出各類型試驗生物降解率的趨勢。圖2至圖4為穩(wěn)健Z比分數(shù)序列圖，直觀地顯示3類試驗中苯甲酸鈉生物降解率的穩(wěn)健Z比分數(shù)分布情況，可清楚地看出每一個試驗苯甲酸鈉降解率與其他試驗的比較結果，且可清晰地辨別出可疑值和離群值。

圖2 OECD 301B試驗穩(wěn)健Z比分數(shù)序列圖Fig. 2 Z-score sequence diagram of OECD 301B test

圖3 OECD 301D試驗穩(wěn)健Z比分數(shù)序列圖Fig. 3 Z-score sequence diagram of OECD 301D test

圖4 OECD 301F試驗穩(wěn)健Z比分數(shù)序列圖Fig. 4 Z-score sequence diagram of OECD 301F test

表3 快速生物降解性試驗質(zhì)量分析對照表Table 3 Comparison table of quality analysis of ready biodegradability test

3 討論與結論(Discussion and conclusion)

穩(wěn)健統(tǒng)計技術通常用于實驗室間或?qū)嶒炇覂?nèi)的能力驗證或結果比對活動，較多是分析樣品對的結果，但也適用于單個樣品的驗證計劃[18]。鑒于本文所涉及的試驗多數(shù)由不同的技術人員操作，因此，在一定程度上可將這些試驗結果看作實驗室內(nèi)部人員比對。

對各類型試驗中苯甲酸鈉的14 d降解率進行穩(wěn)健統(tǒng)計分析的意義在于，通過對過往結果的質(zhì)量分析，測試機構可知在現(xiàn)有的試驗條件下，哪些數(shù)據(jù)是可疑值或離群值。例如，當降解率≤62.3%和≥90.3%(301B)，或≥95.0%(301D)，或≤66.0%(301F)時，SD應提高警覺。若試驗正在進行，應密切關注試驗進程，發(fā)現(xiàn)問題后及時采取糾正措施，避免試驗失敗。若試驗已經(jīng)完成，應仔細檢查分析測試過程的各個環(huán)節(jié)，如溶液配制方法、體系準備、儀器設備的穩(wěn)定性和可靠性及誤差來源等，從技術層面上進行全方位評價，找出產(chǎn)生“可疑”或“離群”的原因，評估該批次試驗數(shù)據(jù)的可接受程度。基于測試機構大量數(shù)據(jù)而開展的穩(wěn)健統(tǒng)計分析和評估，可以反映測試機構的整體數(shù)據(jù)質(zhì)量，反映該測試機構的技術水平，進而獲得提高數(shù)據(jù)質(zhì)量的方法。由此可見，穩(wěn)健統(tǒng)計技術Z比分數(shù)的應用，對測試機構及時發(fā)現(xiàn)質(zhì)量問題、提升數(shù)據(jù)質(zhì)量、提高工作效率，具有重要的應用價值。

開展快速生物降解性試驗時，可將第14天苯甲酸鈉的降解率數(shù)據(jù)與表3進行對比，對該試驗數(shù)據(jù)的質(zhì)量情況做出一個初步判斷。然而，盡管文中涉及的試驗均始終按照標準方法操作，項目實施和原始記錄等都經(jīng)質(zhì)量保證部門(QAU)的檢查并確認符合GLP規(guī)范和SOP規(guī)定，但所得數(shù)據(jù)仍然具有局限性，比如在苯甲酸鈉生物降解過程中起關鍵性作用的接種物僅取自廣州2家生活污水處理廠，缺乏以廣州地區(qū)其他生活污水處理廠和廣州以外地區(qū)生活污水處理廠為接種物來源的相應苯甲酸鈉降解率數(shù)據(jù)。而楊倩等[26]和楊婧等[27]研究指出，接種物易受到污水處理廠本身運行工藝、不同地區(qū)的自然經(jīng)濟環(huán)境以及人們生活習慣等因素影響，微生物群落結構呈現(xiàn)差異性，進而影響到試驗的重復性及結果間的可比性。Ginkel等[28]、Vázquez-Rodríguez等[29]和Dechrugsa等[30]也指出，接種物來源不同，同一個化合物的快速生物降解性會呈現(xiàn)不一致性，因為接種物不同，其組分尤其是微生物種群存在較大差異。另外，上述苯甲酸鈉快速生物降解性試驗僅在本機構完成，未與其他機構進行實驗室間結果比對，不具備普遍適用性。因此如需進行類似數(shù)據(jù)評估，建議根據(jù)實際情況設計一個完整的穩(wěn)健統(tǒng)計分析，進而得到更有針對性和指導性的質(zhì)量分析對照表，從而有效避免出現(xiàn)誤判的情況。

特別指出，由穩(wěn)健Z比分數(shù)分析得出試驗結果可疑甚至離群時，并不等同于該試驗無效，只是指示該試驗系統(tǒng)可能存在嚴重的問題或發(fā)生偏離，此時應針對試驗的各個環(huán)節(jié)開展全面的檢查和深入分析。只有當該試驗結果滿足所有相應的質(zhì)量控制要求時，方能判定試驗為有效。

[1] 環(huán)境保護部化學品登記中心《化學品測試方法》編委會. 化學品測試方法 降解與蓄積卷[M]. 第二版. 北京: 中國環(huán)境出版社, 2013: 25-61

[2] United Nations. Globally Harmonized System of Classification and Labeling of Chemicals (GHS) [M]. 3rd ed. New York: United Nations Publication, 2009: 215-240

[3] Leeuwen C J V, Vermeire T G. Risk Assessment of Chemicals an Introduction [M]. 2nd ed. Dordrecht: Springer, 2007: 73-409

[4] 聶晶磊, 周紅, 刁樹林, 等. 新化學物質(zhì)環(huán)境管理與生態(tài)毒理學[J]. 毒理學雜志, 2005, 19(1): 12-14

[5] 周紅, 沈英娃, 羅桂玲, 等. 加強我國對新化學物質(zhì)的管理[J]. 環(huán)境科學研究, 1999, 12(6): 57-59

Zhou H, Chen Y W, Luo G L, et al. To strengthen the management on new chemicals in China [J]. Research of Environmental Sciences, 1999, 12(6): 57-59 (in Chinese)

[6] Organization for Economic Co-operation and Development. OECD Guideline for Testing of Chemicals [M]. Paris: OECD, 1992: 9-55

[7] OECD. OECD Guidelines for Testing of Chemicals, 301 Ready Biodegradability [S]. Paris: OECD, 1992

[8] 中華人民共和國國家質(zhì)量監(jiān)督檢驗檢疫總局. 中國國家標準化管理委員會. GB/T 21856-2008化學品 快速生物降解性 CO2產(chǎn)生試驗[S]. 北京: 中國標準出版社, 2008

[9] 中華人民共和國國家質(zhì)量監(jiān)督檢驗檢疫總局. 中國國家標準化管理委員會. GB/T 21831-2008 化學品 快速生物降解 密閉瓶法試驗[S]. 北京: 中國標準出版社, 2008

[10] 中華人民共和國國家質(zhì)量監(jiān)督檢驗檢疫總局. 中國國家標準化管理委員會. GB/T 21801-2008化學品 快速生物降解性: 呼吸計量法試驗[S]. 北京: 中國標準出版社, 2008

[11] 奚旦立, 孫裕生, 劉秀英. 環(huán)境監(jiān)測[M]. 第三版. 北京: 高等教育出版社, 2004: 142

[12] 中華人民共和國建設部. CJ/T 221-2005城市污水處理廠污泥檢驗方法[S]. 北京: 中國標準出版社, 2006

[13] Mei C F, Liu Y Z, Long W N, et al. A comparative study of biodegradability of acarcinogenic aromatic amine (4,4′-diamino dipheny lmethane) with OECD 301 test methods [J]. Ecotoxicology and Environmental Safety, 2015, 111: 123-130

[14] 國家環(huán)境保護總局《化學品測試方法》編委會. 化學品測試方法[M]. 北京: 中國環(huán)境科學出版社, 2004: 329-369

[15] OECD. OECD Principles on Good Laboratory Practice [S]. Paris: OECD, 1998

[16] 國家環(huán)境保護總局科技標準司. HJ/T 155-2004化學品測試合格實驗室導則[S]. 北京: 中國環(huán)境科學出版社, 2004

[17] 環(huán)境保護部. HJ 506-2009 水質(zhì) 溶解氧的測定 電化學探頭法[S]. 北京: 中國環(huán)境科學出版社, 2009

[18] 中國合格評定國家認可委員會. 能力驗證結果的統(tǒng)計處理和能力評價指南[S]. 北京:中國合格評定國家認可委員會, 2006

[19] 吳志清. 淺議實驗室對比數(shù)據(jù)分析方法——穩(wěn)健統(tǒng)計Z比分數(shù)[J]. 云南冶金, 2011, 40(增刊): 163-167

Wu Z Q. Discussion on analytical method of lab data comparison -Robust statistic Z-score [J]. Yunnan Metallurgy, 2011, 40(Suppl.): 163-167 (in Chinese)

[20] 張媛, 張磊, 李建, 等. 穩(wěn)健統(tǒng)計方法及t檢驗法在活菌計數(shù)能力比對中的應用[J]. 中國獸藥雜志, 2013, 47(4): 19-22

Zhang Y, Zhang L, Li J, et al. Application of robust statistic method and t test of proficiency testing for live bacteria counting [J]. Chinese Journal of Veterinary Drug, 2013, 47(4): 19-22 (in Chinese)

[21] 邢小茹, 馬小爽, 田文, 等. 實驗室間比對能力驗證中的兩種穩(wěn)健統(tǒng)計技術探討[J]. 中國環(huán)境監(jiān)測, 2011, 27(4): 4-8

Xing X R, Ma X S,Tian W, et al. Two robust statistic techniques in proficiency testing by interlaboratory comparisons [J]. Environmental Monitoring of China, 2011, 27(4): 4-8 (in Chinese)

[22] Todorov V, Filzmoser P. Robust statistic for the one-way MANOVA [J]. Computational Statistics and Data Analysis, 2010, 54: 37-48

[23] 徐楓, 周艷茁. Excel在穩(wěn)健統(tǒng)計Z比分數(shù)計算中的應用[J]. 中國衛(wèi)生檢驗雜志, 2008, 18(8): 1668-1669

[24] 康國強. EXCEL在概率計算中的應用[J]. 安陽大學學報, 2003(1): 55-57

Kang G Q. Excel’application to probabilistic calculation [J]. Journal of Anyang University, 2003(1): 55-57 (in Chinese)

[25] 封躍鵬. 穩(wěn)健統(tǒng)計技術及其在實驗室能力驗證數(shù)據(jù)處理中的應用[J]. 化學分析計量, 2007, 16(2): 55-57

Feng Y P. Robust statistical techniques and its application on proficiency testing inter-laboratory comparisons [J]. Chemical Analysis and Meterage, 2007, 16(2): 55-57 (in Chinese)

[26] 楊倩, 蔣陽月, 王小軍, 等. 不同區(qū)域污水處理廠活性污泥中微生物菌落結構分析[J]. 化工進展, 2014, 33(12): 3329-3336

Yang Q, Jiang Y Y, Wang X J, et al. Diversity of bacterial groups in cultivated sludge samples from different areas of China [J]. Chemical Industry and Engineering Progress, 2014, 33(12): 3329-3336 (in Chinese)

[27] 楊婧, 陳曉倩, 殷浩文. 生物降解實驗中的接種物標準化研究[J]. 環(huán)境工程學報, 2014, 8(5): 2146-2152

Yang J, Chen X Q, Yin H W. Study of inoculum standardization for biodegradation test [J]. Chemical Industry and Engineering Progress, 2014, 8(5): 2146-2152 (in Chinese)

[28] Ginkel C G V, Hann A, Luuten M L G C, et al. Influence of the size and source of the inoculum on biodegradation curves in closed-bottle tests [J]. Ecotoxicology and Environmental Safety, 1995, 31: 218-223

[29] Vázquez-Rodríguez G A, Garabétian F, Rols J L. Inocula from activated sludge for ready biodegradability testing: Homogenization by preconditioning [J]. Chemosphere, 2007, 68: 1447-1454

[30] Dechrugsa S, Kantachote D, Chaiprapat S. Effects of inoculum to substrate ratio, substrate mix ratio and inoculum source on batch co-digestion of grass and pig manure [J]. Bioresource Technology, 2013, 146: 101-108

◆

Application of Robust Statistic Method in Quality Analysis of Ready Biodegradability Test

Deng Guirong1,2,3, Mei Chengfang1,2,3, Liu Yanzhen1,2,3, Zhang Hongtao1,2,3,4, Chen Zefang1,2,3, Zeng Guoqu1,2,3,*

1. Laboratory of Ecotoxicity and Environmental Safety, Guangdong Detection Center of Microbiology, Guangzhou 510070, China 2. Guangdong Provincial Key Laboratory of Microbial Culture Collection and Application, Guangdong Institute of Microbiology, Guangzhou 510070, China 3. State Key Laboratory of Applied Microbiology, Southern China, Guangzhou 510070, China 4. Foshan Research Center for Environmental Health and Safety Assessment, Foshan 528000, China

Received 30 March 2016 accepted 5 May 2016

Ready biodegradability test plays an important role in ecotoxicity testing of new chemicals. The reliability of ready biodegradation test results directly determines the subsequent testing strategy and affects the registration and declaration of new chemical substances, as well as the risk management after entering environment. In order to study the data quality of ready biodegradability series tests (OECD 301B, 301D and 301F), Z-score of robust statistics method was firstly used to analyze the 14 d biodegradation data of sodium benzoate which were produced by our laboratory in the past 5 years. The results showed that, when the biodegradation rates of sodium benzoate of 301B, 301D and 301F are within the range of 63.8%-88.0%, 63.4%-89.2% and 66.0%-99.8%, respectively, it indicates that the test results are satisfactory and the data quality are good. Otherwise, it means that there are some problems with the data and further verification is needed. As a result, a comparison table of quality analysis of ready biodegradability test was constructed for providing the references for quality analysis of ready biodegradability test of new chemical substances.

sodium benzoate; ready biodegradability; robust statistic; Z-score; quality analysis

新材料環(huán)境健康安全評價科技創(chuàng)新服務平臺(廣東省科技計劃項目，2013B090800004)

鄧桂榮(1984- )，男，工程師，研究方向為生態(tài)毒理學，E-mail: deng8818@163.com；

*通訊作者(Corresponding author)， E-mail: zenggq@gddcm.com

10.7524/AJE.1673-5897.20160330004

2016-03-30 錄用日期：2016-05-05

1673-5897(2016)6-259-07

X171.5

A

曾國驅(qū)(1966-)，男，博士，研究員，主要研究方向為生態(tài)毒理學，發(fā)表學術論文30余篇。

鄧桂榮, 梅承芳, 柳燕貞, 等. 穩(wěn)健統(tǒng)計在快速生物降解性試驗質(zhì)量分析中的應用[J]. 生態(tài)毒理學報，2016, 11(6): 259-265

Deng G R, Mei C F, Liu Y Z, et al. Application of robust statistic method in quality analysis of ready biodegradability test [J]. Asian Journal of Ecotoxicology, 2016, 11(6): 259-265 (in Chinese)