膠州灣石油類污染狀況及其環境自凈能力評估?

李　巖， 李克強， 馬云鵬， 梁生康， 王修林

(1. 中國海洋大學環境科學與工程學院， 山東 青島 266100；

2. 中國海洋大學化學化工學院，海洋化學理論與工程技術教育部重點實驗室， 山東 青島266100)

?

膠州灣石油類污染狀況及其環境自凈能力評估?

李巖1， 李克強2??， 馬云鵬2， 梁生康2， 王修林2

(1. 中國海洋大學環境科學與工程學院， 山東 青島 266100；

2. 中國海洋大學化學化工學院，海洋化學理論與工程技術教育部重點實驗室， 山東 青島266100)

摘要：本文根據調查分析了膠州灣2007和2015年石油類污染狀況，并應用膠州灣石油類污染物三維水動力-水質模型，評估了膠州灣石油類環境自凈能力。結果表明，2007年8和11月膠州灣石油類平均濃度分別約為(72±32)和(70±41) mg·m-3，超標率分別為68%和57%，而2015年5和8月則分別約為(38±17)和(150±46) mg·m-3，超標率分別為60%和100%，分布上均在東北部近岸海域污染較嚴重。膠州灣石油類環境自凈能力約為2 100 t·a-1，以水動力自凈為主，約占54%，其次為大氣揮發自凈和微生物降解自凈。膠州灣石油類自凈能力有明顯的季節變化，秋季自凈能力較大，夏季較小，主要受潮汐和季風等因素影響。

關鍵詞：石油類；自凈能力；污染狀況；模型；膠州灣

Li Yan,Li Ke-Qiang,Ma Yun-Peng,et al. Pollution condition and estimation of self-purification capacity of petroleum in Jiaozhou Bay[J].Periodical of Ocean University of China, 2016, 46(2): 83-89.

膠州灣位于黃海之濱，山東半島南岸，是中國北方海區典型的半封閉海灣，水域面積約320km2，灣口狹小，最窄處僅有3.1km，整個海灣水深較淺，平均8.8m[1]。隨著青島市經濟高速發展、海上活動日益增強等，海上船舶活動和陸源排污致使石油類成為膠州灣海域重要的污染物之一，甚至在1980年代中后期致使膠州灣成為中國近海海域石油類污染最為嚴重的海域之一[2-4]，近年來盡管有所降低，但某些年份仍然是主要污染指標[5-7]，對于海域生態系統健康影響不容忽視，特別是溢油事故所引發的生態災害問題[3]。至今很少見到有關膠州灣海域石油類污染分布的報道，這在很大程度上制約了人們對膠州灣海洋生態環境的整體了解。同時，為了對膠州灣石油類污染進行有效管理，迫切需要了解其海洋自凈能力。

水動力輸運擴散過程、大氣揮發自凈過程和微生物降解自凈過程是石油類自凈的主要途徑，通過水動力-水質模型可以評估石油類環境自凈能力[3]。其中，水動力模型主要包括箱式潮交換模型[8-9]、標識質點跟蹤模型[10-11]和對流-擴散輸運模型[12-13]等，石油類水質模型主要包括箱式模型[3，14]和三維水動力-水質模型[15]。其中，箱式和標識質點跟蹤模型因不能刻畫水交換時空結構或擴散過程，往往高估或低估水動力自凈能力，而水動力對流-擴散輸運模型只能評估水動力自凈能力，只有三維水動力-水質模型能夠比較精確、完備地描述由于對流、擴散等所產生的污染物水動力輸運過程以及其他生物、化學自凈過程，廣泛應用于海灣、河口等近海海域自凈能力研究[15-18]。應用上述模型，膠州灣海域開展了水半交換時間、存留時間等水動力自凈能力等方面的研究[11,18-20]，也有分別基于箱式模型的石油類環境容量研究[14]和基于三維水動力-水質模型的石油類分配容量研究[15]，前者因水動力過程過于簡單而不能準確評估水動力自凈能力，而后者只有排污口分配容量結果，沒有對自凈能力進行研究。為此，分別于2007年8、11月和2015年5、8月在膠州灣進行了石油類分布調查，并在此基礎上，根據膠州灣石油類三維水動力-水質模型，評估了膠州灣石油類環境自凈能力。可為進一步加強膠州灣海域環境管理，實施“污染物排放總量控制”提供科學依據。

1材料與方法

1.1 調查站位和方法

2007年8月21日、11月24日和2015年5月16日、8月19日在膠州灣進行了4個航次污染物陸海同步調查，對石油類的含量進行了測定。海上共調查了7個站位，陸上共調查了7個入海河口和2個污水處理廠排污口(見圖1)。調查中采用拋浮式無油玻璃采水器采集表層海水樣品用于石油類分析，樣品采集后立即用1∶3 H2SO4調節樣品至弱酸性(pH ≈ 4)，然后用10 mL正己烷分2次萃取，萃取液密封后在(5±2)℃條件下避光保存，在實驗室中應用紫外分光光度法[17]分析石油類濃度，平均相對標準偏差為16%(n=90)。

圖1　膠州灣地理位置和海上調查站位(+)及

1.2 石油類自凈能力評估方法

基于普林斯頓大學河口、陸架和海洋模式(ECOM)[18-19]所構建的膠州灣三維對流擴散數值模型[20-22]，結合膠州灣石油類多介質動力學模型[10]，可構建膠州灣石油類三維水動力-水質模型[11]。在水動力及大氣揮發、生物降解等生物地球化學過程作用下，石油類從水體中去除而自凈[23]，自凈能力可通過積分自凈過程計算[24]。

SPC=

(1)

式中：右邊第一項為石油類大氣揮發自凈量；第二項為生物降解自凈量；第三項為水動力自凈量。

這樣，通過設定石油類陸源岸線水質控制點位置、控制濃度標準，通過方程1，采用積分方法，可計算評估石油類自凈能力。

1.3 邊界條件和初始條件

模式中以大沽河、洋河、墨水河、樓山河、板橋坊河、李村河和海泊河以及廉灣和團島污水處理廠作為陸源源強(見圖1)，其年徑流量和入海石油類量按照實際監測值加入，2007年陸源石油類入海共約687t·a-1(見表1)。海源主要是船舶廢水排放，采用估算值，約600t·a-1[3]，合并到大港碼頭作為海源源強。

模型模擬校驗陸地邊界：河流入海口和排污口采用實際監測通量，陸地邊界污染物通量為零，法向梯度為零。

自凈能力評估陸地邊界：河流和陸地邊界污染物通量均直接給定石油類濃度，具體采用功能區國家海水水質標準。

開邊界：C=0，入流段；

初始條件：t=0時，C=0。

表1　環膠州灣2007年主要河流及排污口徑流量和入海石油類量

2結果與討論

2.1 石油類污染狀況及其分布特征

根據2007年海水監測結果(見圖2)表明，2007年8和11月膠州灣表層海水中石油類平均濃度較為接近，分別約為(72±32)和(70±41)mg·m-3，變化幅度分別為32～135和32～158mg·m-3，以國家一、二類海水水質標準(50mg·m-3)為評價標準，污染指數(Ai)[25]分別平均為1.5±0.6和1.4±0.8，超標率分別為68%和57%，污染較為嚴重，檢出率均為100%。濃度分布總體上均呈現出由膠州灣東北部向西南部、灣口遞減的趨勢，灣內大部分區域石油類濃度大于40mg·m-3，大于100mg·m-3的高值區主要集中在膠州灣東北部附近，其中樓山河口附近海域濃度最高(見圖2)，與公報結果基本一致[5]。

圖2　2007年8月和11月膠州灣海水中石油

2015年海水監測結果(見圖3)表明，2015年5和8月膠州灣表層海水中石油類平均濃度差別較大，分別約為(38±17)和(150±46) mg·m-3，變化幅度分別為26～90和56～230mg·m-3，以國家一、二類海水水質標準為評價標準，污染指數(Ai)分別平均為0.8±0.3和3.0±0.9，超標率分別約為60%和100%，8月份污染嚴重，檢出率均為100%。2015年8月份石油類濃度較2007年有所升高，但濃度分布整體趨勢上與2007年基本一致，只是局部有所區別，如5月份西南部較中部濃度高，而8月中部和西北部出現濃度高值區，灣內大部分區域石油類濃度大于40mg·m-3，特別是8月份均超過國家一、二類海水水質標準，高值區同2007年，也集中在膠州灣東北部附近，其中李村河口附近海域濃度最高(見圖3)。海上石油類的時空分布與其來源密切相關，文獻結果表明，膠州灣石油類主要來自與陸源排污和海上船舶污染，東北部的李村河、墨水河以及西北部的大沽河是石油類的主要陸源，而灣東部的碼頭及航道船舶航運是主要海源[3]。

圖3　2015年5月和8月膠州灣海水中石油

2.2 石油類環境自凈能力評估

2.2.1 膠州灣石油類濃度分布模擬應用膠州灣石油類三維水動力-水質模型，經過180天左右的計算，膠州灣內的石油類濃度基本達到動態平衡。達到平衡后即為1月1日，然后模擬計算1a，石油類按照2007年陸源監測值輸入，取同一潮周期的石油類濃度時間序列進行潮周期平均，從而得到膠州灣石油類2007年平均濃度的分布(見圖4)，與監測值相比分布趨勢基本一致。這說明，該模式能近似反映石油類的輸運、擴散過程，可以用來評估石油類環境自凈能力。

圖4　膠州灣表層石油類濃度模擬結果分布(單位：mg·m-3)

2.2.2 膠州灣石油類環境自凈能力通過模式運算，以國家一、二類海水水質為控制標準，膠州灣石油類環境自凈能力約為2100t·a-1，與文獻結果基本一致[10-11]，較箱式模型計算結果高約600t·a-1[10]，主要來自水動力自凈貢獻的增量。其中，水動力自凈能力最大，約為1100t·a-1，約占54%，其次為大氣揮發，約700t·a-1，約占33%，而生物降解最小，僅占13%(見圖5)。受水動力輸運過程影響，膠州灣石油類自凈能力呈現明顯的季節變化(見圖6)，其中，秋季自凈能力較大，平均可達6.7t·d-1，夏季自凈能力較小，平均約為5.6t·d-1；每月隨著膠州灣大小潮的不同，自凈容量也有大小之別，大潮自凈容量大，小潮自凈容量小。這可能與膠州灣的潮汐和季風等因素季節變化有關。

2.2.3 影響膠州灣石油類環境自凈能力的因素受潮汐的月不等影響，膠州灣海域潮汐每月有大小潮之分，通過模式可以計算一個月的石油類水動力自凈量，將一個月中每天的水動力自凈量除以日平均值，可定義潮汐影響因子(KT)，即潮汐對膠州灣石油類水動力自凈能力的影響。結果表明，水動力自凈能力隨著潮汐月變化而變化，存在月極大值和極小值(見圖7)。其中，大潮時期膠州灣納潮量大，海水交換量大，自凈能力大；小潮時期海水交換量小，自凈能力小。結果表明，KT對石油類水動力自凈能力的影響變化范圍約為0.4～1.5，可用三角函數描述：

KT=0.99+0.51·sin(π·(t+9.9)/7.5)。

(2)

圖5　膠州灣石油類生物降解、大氣揮發和

圖6　膠州灣石油類環境自凈能力(SPC)季節變化

季風取團島歷年統計資料[26]，通過模式中引入季風和無季風影響下對石油類水動力自凈量的計算，將每日季風影響下的石油類水動力自凈量與無季風影響下的石油類水動力自凈量相比，定義為季風影響因子(KW)。當KW>1時，表示季風對水動力自凈起促進作用，當KW<1時，表示季風對水動力自凈起阻礙作用。圖8為季風影響因子的季節變化，結果表明，秋冬季節KW>1，說明石油類水動力自凈受到了促進，夏季KW<1，說明石油類水動力自凈受到阻礙。進一步分析表明，膠州灣口朝向東南，膠州灣常年秋冬季節以西北風為主，在風力的驅動下加強了灣內水與黃海水的交換，而夏季以南風為主，不利于海水交換。結果表明，季風影響因子(KW)對石油類水動力自凈能力的影響變化范圍約為0.98～1.12，可用三角函數描述：

KW=0.988+0.13·(sin(π·(t-165)/360))2-0.02·cos(2·π·(t-200)/160)。

(3)

圖7　潮汐影響因子的變化

圖8　季風影響因子的季節變化

4結語

本文利用2007和2015年監測數據對膠州灣石油類污染狀況進行了評價，結果表明，2015年石油類污染較2007年有所增加，超過國家海水一、二類海水水質標準的范圍由約60%增加到100%。應用膠州灣石油類三維水動力-水質模型，對膠州灣石油類自凈能力進行了評估。結果表明，膠州灣石油類自凈能力約為2100t·a-1，以水動力自凈為主，并有明顯的季節變化，受潮汐變化和季風的影響，表現出秋季石油類自凈能力大，夏季自凈能力小的特征，每月石油類水動力自凈容量有大小潮之分，大潮水動力自凈能力大，小潮水動力自凈能力小。其結果可為膠州灣石油類環境管理提供科學依據。

致謝本文水動力模型由中國海洋大學海洋環境學院鮑獻文老師提供，在此對鮑獻文老師及萬修全博士表示感謝！相關計算工作在中國海洋大學計算服務中心完成。

參考文獻：

[1]《中國海灣志》編纂委員會．《中國海灣志》(第四分冊)[M]．北京：海洋出版社, 1993： 157-260.

Editorial board of China bay survey. Survey of China bays (Vol. 4) [M]. Beijing: China Ocean Press, 1993： 157-260.

[2]國家海洋局. 二十世紀末中國海洋環境質量公報[C].北京：國家海洋局. 2000： 9-10

State Oceanic Administration. Marine environmental quality communique of China at the end of 20th century [C]. 2000： 9-10.

[3]王修林, 李克強, 石曉勇. 膠州灣主要化學污染物海洋環境容量[M]. 北京：科學出版社, 2006.

WANG Xiulin, LI Keqiang, SHI Xiaoyong. Marine environmental capacity of pollutants in Jiaozhou Bay [M]. Beijing: Science Press, 2006.

[4]孫耀, 崔毅, 于宏, 等. 膠州灣表層水中石油烴、化學需氧量和金屬鉻的分布及污染現狀分析[J]. 海洋水產研究, 1992, 13：123-130.

SUN Yao, CUI Yi, YU Hong, et al. Distribution and pollution status for petroleum, chemical oxygen demand and chromium in the surface water of Jiaozhou Bay [J]. Marine Fisheries Research, 1992, 13：123-130.

[5]國家海洋局. 2001-2014年中國海洋環境狀況公報 [R]. 2002-2015.

State Oceanic Administration. Marine environmental quality communique of China in 2001-2014 [R]. 2002-2015, 3.

[6]青島市環境保護局. 2001-2014年青島市環境狀況公報 [R].青島：青島市環境保護局， 2002-2015.

Environmental Protection Bureau of Qingdao. Report on the State of Qingdao's Environment in 2001-2014 [R].2002-2015, 5.

[7]張志強. 青島膠州灣石油烴和其它環境因子的時空分布及其相關性研究[D]. 青島： 中國海洋大學， 2005.

ZHANG Zhiqiang. Temporal and Spatial Distribution of Petroleum Hydrocarbons and Other Environmental Factors and Their Relativity in Jiaozhou Bay [D]. Qingdao： Ocean University of China, 2005.

[8]Zimmerman J T F. Mixing and flushing of tidal embayments in the Western Duch Wadden Sea [J]. Neth J Sea Res, 1976, 10(4): 397-439.

[9]匡國瑞, 楊殿榮, 喻祖祥, 等. 海灣水交換的研究-乳山東灣環境容量初步探討[J]. 海洋環境科學, 1987, 6(1)：13-23.

KUANG Guorui, YANG Dianrong, Yu Zuxiang, et al. Study on the water exchange of bay-the initial exploration on the environmental capacity in East Rushan Bay [J]. Marine Environmental Science, 1987, 6(1): 13-23.

[10]Signell R P, Butman B. Modeling tidal exchange and dispersion in Bostan Harbor [J]. Geophys Res, 1992, 97(C10): 15, 591-606.

[11]趙亮, 魏皓, 趙建中. 膠州灣水交換的數值研究[J]．海洋與湖沼, 2002, 33(1)：23-29.

ZHAO Lian, WEI Hao, ZHAO Jianzhong. Numerical study on water exchange in Jiaozhou Bay [J]．Oceanologia Et Limnologia Sinic, 2002, 33(1): 23-29.

[12]董禮先, 蘇紀蘭. 象山港水交換數值研究Ⅱ-模型應用和水交換研究[J]．海洋與湖沼, 1999, 30(5)：465-470.

DONG Lixian, SU Jilan. Numerical study of water exchange in Xiangshangang Bay II. Model application and water exchange study [J]．Oceanologia Et Limnologia Sinic, 1999, 30(5): 465-470.

[13]孫英蘭, 張越美. 丁字灣物質輸運及水交換能力研究[J]．青島海洋大學學報, 2003, 33(1)：1-6.

SUN Yinglan, ZHANG Yuemei. A numerical model of pollutant transport and seawater exhange in Dingzi Bay [J]．Journal of Ocean University Of Qingdao, 2003, 33(1): 1-6.

[14]李克強, 王修林, 石曉勇, 等. 石油烴在膠州灣多介質中遷移-轉化模型研究[J]. 海洋環境科學, 2009, 28(1): 12-16.

LI Keqiang, WANG Xiulin, SHI Xiaoyong, et al. Study on model of transport and transformation of petroleum hydrocarbons in multimedia environment in Jiaozhou Bay [J]. Marine Environmental Science, 2009, 28(1): 12-16.

[15]QIAN Guodong, SHI Xiaoyong, LI Keqiang, et al. Estimation of environmental capacity of petroleum hydrocarbons in jiaozhou bay, China [J]. Fresenius Environmental Bulletin， 2012, 21(1): 48-53.

[16]LI Keqiang, SU Ying, YING Jun, et al. Environmental capacity of petroleum hydrocarbon pollutants in Jiaozhou Bay, China: modeling and calculation[J]. J Ocean Univ China, 2013, 12 (1): 70-76.

[17]姜太良, 宋萬先, 房憲英.萊州灣西南部的物理自凈能力[J]. 海洋通報, 1991, 10(2)：53-79.

JIANG Tailiang, SONG Wanxian, FANG Xianying. Physical self-purification capacity at the southwest part of Laizhou Bay [J]. Marine Science Bulletin, 1991, 10(2): 53-79.

[18]LIU Zhe, WEI Hao, LIU Guangshan, et al. Simulation of water exchange in Jiaozhou Bay by average residence time approach [J]. Estuarine, Coastal and Shelf Science, 2004, 61: 25-35.

[19]WANG Cui, ZHANG Xueqing, SUN Yinglan. Numerical simulation of water exchange characteristics of the Jiaozhou Bay based on A three-dimensional lagrangian [J]. Model China Ocean Engineering, 2009, 23(2): 277-290.

[20]吳永成, 王從敏. 膠州灣海水交換和混合擴散[M]. //劉瑞玉主編. 膠州灣生態學和生物資源. 北京: 科學出版社, 1992, 57-72.

WU Yongcheng, WANG Congmin. Water change and turbulent mixing and diffusion of Jiaozhou Bay. In LIU Ruiyu, ed., Ecology and Living Resources of Jiaozhou Bay[M]. Beijing: Science Press, 1992： 57-72.

[21]國家海洋環境監測中心. 海洋監測規范第4部分: 海水分析GB 17378.4-2007[S]． 北京： 中國標準出版社， 2008.

National Marine Environmental Monitoring Center of China. The specification for marine monitoring Part 4: seawater analysis GB 17378. 4-2007[S]． Beijing： Standards Press of China， 2008.

[22]Blumberg A F, Herring H J. Circulation Medal using orthogonal curvilinear coordinates [M]. // Jacqes J, Nihoul B J, Eds.Three-Dimension Models of Marine and Estuarine Dynamics. Amsterdam: Elsevier Scinence Publishers B V, 1987: 55-88.

[23]Blumberg A F, Mellor G L. A description of a three-dimensional coastal ocean circulation model [M]. // Heaps N S Ed. Three-dimensional Coastal Ocean Models. Coastal and Estuarine Sciences, vol. 4. Washington: American Geophys Union, 1987: 1-16.

[24]鮑獻文, 閆菊, 趙亮,等. ECOM模式在膠州灣潮流中的應用[J]. 海洋科學, 1999(5): 57-60.

BAO Xianwen, YAN Ju, ZHAO Lian, et al. Application of ECOM to simulate tidal currents in Jiaozhou Bay [J]. Marine Science, 1999(5): 57-60.

[25]萬修全, 鮑獻文, 吳德星, 等. 膠州灣及鄰近海域潮流和污染物擴散的數值模擬[J]. 海洋科學, 2003, 27(5)：31-36.

WAN Xiuquan, BAO Xianwen, WU Dexing, et al. Numerical simulation of the tidal currents and the pollutant diffusion in Jiaozhou Bay [J]. Marine Science, 2003, 27(5): 31-36.

[26]閆菊, 王海, 鮑獻文. 膠州灣三維潮流及潮致余環流的數值模擬[J]．地球科學進展, 2001, 16(2)：172-177.

YAN Ju, WANG Hai, BAO Xianwen. The simulation of3d tidal and residual current in the Jiaozhou Bay [J]．Advance in Earth Sciences, 2001, 16(2): 172-177.

[27]王修林, 石曉勇, 江玉, 等. 影響海洋環境質量的重要海洋過程研究[C]//98青年海洋論壇-海洋可持續發展論文集. 北京： 海洋出版社, 1998, 121-125.

WANG Xiulin, SHI Xiaoyong, JIANG Yu, et al. Study on the important marine processes affecting the marine environment quality [C]. The 98 Youth Marine Forum-Marine Sustainable Development Symposium. Beijing: Marine Press, 1998： 121-125.

[28]LI Keqiang, WANG Xiulin. Calculation methodology of marine environmental capacity for heavy metal: A case study in Jiaozhou Bay, China[J]. Chinese Science Bulletin, 2013, 58(2): 282-287.

[29]盧冰, 熊健, 鄭士龍．象山灣海域中石油烴含量分布特征及其評價[J]．東海海洋, 1993, 11(4)： 67-74.

LU Bing, XIONG Jian, ZHENG Shilong. The distribution of petroleum carbohydrate in the region of the Xiangshan Bay [J]．Donghai Marine Science, 1993, 11(4)： 67-74.

[30]國家海洋局第一海洋研究所港灣室《膠州灣自然環境》編寫組. 膠州灣自然環境[M]. 北京：海洋出版社, 1984: 3-21.

Editorial board of Natural Environment of Jiaozhou Bay, the harbor group of the First Institute of Oceanography, State Oceanic Administration. Natural Environment of Jiaozhou Bay [M]. Beijing: Marine Press, 1984: 3-21.

責任編輯徐環

Pollution Condition and Estimation of Self-Purification Capacity of Petroleum in Jiaozhou Bay

LI Yan1, LI Ke-Qiang2, MA Yun-Peng2, LIANG Sheng-Kang2, WANG Xiu-Lin2

(1. College of Environmental Science and Engineering, Ocean University of China, Qingdao 266100, China; 2. The Key Laboratory of Marine Chemistry Theory and Technology, MOE, College of Chemistry and Chemical Engineering, Ocean University of China, Qingdao 266100, China)

Abstract:Historically, Jiaozhou Bay was badly contaminated by petroleum, and it was a slight improvement recently. The pollution condition of petroleum was studied based on the survey in 2007 and 2015 in Jiaozhou Bay. The results show that the mean concentrations of petroleum were (72±32) and (70±41)mg·m-3in August and November, 2007, respectively, excessing the seawater quality criteria with 68% and 57%. They were (38±17) and (150±46)mg·m-3in May and August, 2015, respectively, excessing the seawater quality criteria with 60% and 100%, where northeastern part of the bay was contaminated badly. The self-purification capacity of petroleum in Jiaozhou Bay was estimated with a 3 dimensional hydrodynamical petroleum transportation model. The results illustrate that the self-purification capacity of petroleum in Jiaozhou Bay was ～2100t·a-1, higher for hydrodynamical self-purification, accounting for ～54%, and lower for volatilization and biodegradation self-purification. There was seasonal variation about hydrodynamical self-purification capacity, and it occured more in autumn, with average 6.7t·d-1, less in summer, with average 5.6t·d-1， which was mainly affected by tide and wind, with influencing factor ranging from 0.4 to 1.5 and from 0.98 to 1.12, respectively. The results can suplly scientific basis for environmental management in Jiaozhou Bay.

Key words:petroleum; self-purification capacity; pollution; dynamical model; Jiaozhou Bay

DOI：10.16441/j.cnki.hdxb.20150312

中圖法分類號：X820.3

文獻標志碼：A

文章編號：1672-5174(2016)02-083-07

作者簡介：李巖(1980-)，男，博士生，海洋環境規劃與管理方向。E-mail：liyan12@ouc.edu.cn??通訊作者: E-mail:likeqiang@ouc.edu.cn

收稿日期：2015-09-07；

修訂日期：2015-10-21

基金項目：?國家自然科學基金委員會-山東省人民政府聯合資助海洋科學研究中心項目(U1406403)資助

引用格式：李巖，李克強，馬云鵬，等.膠州灣石油類污染狀況及其環境自凈能力評估[J].中國海洋大學學報(自然科學版)， 2016， 46(2)： 83-89.

Supported by NSFC-Shandong Joint Fund for Marine Science Research Centers (U1406403 )