劉弢，張亦飛，祁琪，方欣，郝春玲，李曉燕，楊輝

(國家海洋局第二海洋研究所工程海洋學重點實驗室 杭州 310012)

?

基于相互作用矩陣的象山港圍填海適宜性評價*

劉弢，張亦飛，祁琪，方欣，郝春玲，李曉燕，楊輝

(國家海洋局第二海洋研究所工程海洋學重點實驗室 杭州 310012)

圍填?；顒邮艿缴鐣洕l(fā)展的驅動，其實施會對海洋生態(tài)環(huán)境產生較大的影響，多種因素之間的相互作用，構成了一類復雜系統(tǒng)。針對象山港海域4個工業(yè)與城鎮(zhèn)建設區(qū)的適宜性評價問題，文章引入巖石工程系統(tǒng)(rock engineering systems，RES)理論的相互作用矩陣計算圍填海驅動因素及生態(tài)環(huán)境影響因子的權重，采用加權模型計算評價對象的適宜性。結果表明，基于RES相互作用矩陣的權重確定方法具有合理性和可行性，適宜性評價結果與實際情況相符，可以為其他區(qū)域圍填海適宜性評價提供參考。

象山港；圍填海；相互作用矩陣；指標權重；適宜性評價

隨著我國沿海經濟的高速發(fā)展和人口壓力的不斷增大，土地資源的短缺成為制約沿海經濟發(fā)展的重要“瓶頸”，而填海造地是緩解沿海地區(qū)土地供求矛盾，拓展生存空間和生產空間的有效手段。但圍填?；顒釉趲砩鐣洕б娴耐瑫r，也引起了一系列負面影響，例如：改變海域水動力條件，造成了港口海灣淤積；加速沿海灘涂濕地生態(tài)系統(tǒng)功能退化；導致生物棲息地喪失，生物多樣性降低；造成海水水質惡化等。另外，圍填海工程不同程度地影響了漁業(yè)、港航及旅游等資源效用的發(fā)揮， 使原本就存在的海灣資源利用矛盾和沖突更加激烈，加劇了資源的利用沖突和產業(yè)之間的矛盾[1-3]。因此，綜合考慮圍填海對海洋生態(tài)環(huán)境和社會經濟的影響，科學合理地進行圍海造地活動，對沿海區(qū)域社會、經濟和環(huán)境協(xié)調發(fā)展具有重要意義。

目前，國內相關專家對圍填海適宜性的評價已開展了一些研究工作。于永海等在海岸自然條件、海洋生態(tài)條件綜合分析的基礎上，綜合考慮開發(fā)利用現(xiàn)狀、災害地質、社會經濟等因素，建立了海岸圍填海適宜性定量評估的指標體系和相應的評價方法[4]。王初升等依據(jù)珊瑚礁和紅樹林生物海岸的資源和生態(tài)特點，通過分析識別近岸圍填海對珊瑚礁和紅樹林生態(tài)影響的主要因子和生態(tài)服務功能的影響，分別構建了珊瑚礁生物海岸和紅樹林生物海岸圍填海適宜性的評價指標體系[5-6]。劉大海等基于環(huán)境經濟學和統(tǒng)計學的基礎理論和評估方法，從海岸、動力、沖淤、生物和水質等 5個方面構建了適用于我國海灣的圍填海適宜性評價體系[7]。王靜等運用多目標決策理論與方法， 綜合考慮圍填海對動力泥沙環(huán)境、海洋生態(tài)環(huán)境、資源綜合開發(fā)和社會經濟影響， 建立了適宜的圍填海規(guī)模評價指標體系[8]。圍填?；顒邮艿缴鐣洕l(fā)展的驅動，其施工和運營對海洋環(huán)境產生較大的影響，多種因素相互聯(lián)系、相互作用，構成一類復雜系統(tǒng)。而現(xiàn)有的研究成果在這方面考慮顯得相對不足，存在進一步探討和改進的空間。

本文引入巖石工程系統(tǒng)(Rock Engineering Systems，RES)理論，并利用RES理論的相互作用矩陣，分析各評價因子的相互作用關系及作用程度，計算各評價指標的權重，同時，以《浙江省海洋功能區(qū)劃(2011—2020)》確定的象山港海域4個工業(yè)與城鎮(zhèn)建設區(qū)作為評價對象，進行圍填海適宜性評估。

1 相互作用矩陣的基本原理

巖石工程系統(tǒng)(RES)理論由J.A. Hudson 教授于1992年提出，其基本思想是把巖石工程看做是完整系統(tǒng)，通過列舉與系統(tǒng)相關的變量構造相互作用矩陣，并通過對矩陣編碼定量評價各變量間的相互作用[9]。該方法不僅考慮了系統(tǒng)中每個參數(shù)對系統(tǒng)的影響，同時考慮了參數(shù)之間的相互作用對整個系統(tǒng)的貢獻，特別適合于解決多種因素之間的相互作用和相互耦合復雜問題的研究[10]。

以圖1所示的3×3矩陣為例，將所有與問題相關的主要因素列于方陣的主對角線上(A→C)，這些因素的狀態(tài)變量具有概念化屬性。矩陣非對角線上的空格被賦予一個數(shù)值以描述兩個參數(shù)間的相互影響程度，這些值可以通過數(shù)值分析、經驗判斷或試驗研究等方法獲得。相互作用矩陣具有非對稱性，如A對B的影響可以不同于B對A的影響。矩陣建立后嚴格按照順時針方向進行影響的分析與賦值。

圖1 相互作用矩陣的基本原理

相互作用矩陣建立和編碼后，將矩陣行和列的編碼值相加就可以確定每個變量在系統(tǒng)中的相互作用程度。每行的總和表示位于該行主對角線上的參數(shù)對系統(tǒng)產生的影響，每列的總和表示系統(tǒng)對位于該列主對角線上的參數(shù)的影響。相互作用程度越高，說明該參數(shù)在系統(tǒng)中的重要性越大，因此可以利用相互作用矩陣確定所研究系統(tǒng)中各相關參數(shù)的權重。參數(shù)的權重Wi按照式(1)至式(3)計算：

(1)

式(1)中：

(2)

(3)

式(2)和式(3)中：Vi，j表示參數(shù)i對參數(shù)j的作用程度。

目前，圍填海適宜性評價過程中指標權重的確定多使用層次分析法(Analytic Hierarchy Process，AHP)。AHP方法是建立在判斷矩陣為一致性矩陣基礎上的，實際操作中當判斷矩陣階數(shù)n>3時，往往難以滿足一致性的要求，此時應用AHP方法就顯得較困難。相互作用矩陣通過兩個之間的相互作用分析，既考慮了單一因素對所研究系統(tǒng)的影響，也考慮了參數(shù)之間的相互作用對整個系統(tǒng)的貢獻，同時能夠大大降低計算工作量，為復雜系統(tǒng)分析提供了一種更為快捷、便利的權重確定方法，與AHP方法相比具有一定的優(yōu)越性。

2 研究區(qū)域和評價對象概況

象山港位于浙江北部沿海， 北鄰杭州灣， 南靠三門灣， 東側為舟山群島， 是一個NE-SW走向的狹長形半封閉海灣。沿灣分布有寧波市的象山、寧海、奉化、鄞州和北侖5個縣(區(qū))，海域面積391.76 km2，海域大部分水深在10～20 m 之間，最深達到47 m，灘涂面積約171 km2。該海灣是浙江省重要的海灣之一，具有良好的區(qū)位優(yōu)勢和資源優(yōu)勢。隨著區(qū)域經濟社會的快速發(fā)展，多個涉海行業(yè)參與海灣資源的開發(fā)，開發(fā)強度不斷增大，用?；顒又g的矛盾突出。同時由于象山港獨特的地形條件和水動力條件，海水水體交換能力較弱，環(huán)境容量有限，海域生態(tài)環(huán)境日趨惡化。選擇《浙江省海洋功能區(qū)劃(2011—2020)》確定的象山港海域4個工業(yè)與城鎮(zhèn)建設區(qū)作為對象，進行適宜性評價。評價對象基本情況見表1。

表1 象山港工業(yè)與城鎮(zhèn)建設區(qū)基本情況

3 評價指標選擇和模型構建

3.1 指標選取

遵循指標選取的科學性、系統(tǒng)性、綜合性、層次性、區(qū)域性、動態(tài)性原則，從社會經濟、生態(tài)、環(huán)境和對周邊海洋功能區(qū)的影響4個層面，結合象山港海域資源的環(huán)境特點，選擇影響圍填海適宜性的主要指標。

(1)區(qū)域經濟社會指標。按照資源經濟學的觀點，自然資源消耗的根本驅動因素是人口增加與經濟發(fā)展[11]，為此，選擇功能區(qū)所在縣(市、區(qū))的經濟規(guī)模(GDP)、人口數(shù)量作為評價指標，采用2012年寧波市經濟社會發(fā)展統(tǒng)計數(shù)據(jù)作為資料來源。

(2)生態(tài)指標。眾所周知，圍填海是對海域生態(tài)環(huán)境影響最大的一種用?；顒樱瑢Ｑ笊锂a生嚴重影響，特別是潮間帶底棲生物。圍填海占用的海域，其底質環(huán)境被覆蓋，絕大部分底棲生物被掩埋、覆蓋而死亡，對潮間帶和底棲生物群落的破壞具有不可逆轉性[12]，為此選取城鎮(zhèn)與工業(yè)建設區(qū)所在海域的底棲生物生物量作為評價指標。

(3)環(huán)境指標。海灣圍填?；顒又苯佑绊懼苓吅Ｓ虻沫h(huán)境容量和泥沙沖淤，并且研究區(qū)域地處狹長形半封閉港灣，水交換能力弱，環(huán)境容量有限，多年來氮、磷等指標嚴重超標，所以環(huán)境指標方面選取納潮量、無機氮含量、活性磷酸鹽含量以及沖淤速率的變化量。評價時，采用國家海洋公益性行業(yè)科研專項項目“東海狹長形海灣綜合整治集成技術及示范應用研究”現(xiàn)場補充調查得到的相關數(shù)據(jù)。

(4)對周邊功能區(qū)的影響。除了對直接覆蓋海域生態(tài)環(huán)境產生影響外，圍填?；顒右矔χ苓吅Ｓ蛎舾袇^(qū)域及利益相關者產生影響，比如碼頭群、航道、海灣口門，濱海旅游資源、海洋保護區(qū)等。如果周邊海域越敏感、利益相關者越重要，越需要采用多種措施控制圍填?；顒訉λ鼈儺a生的負面影響。評價時，根據(jù)象山港海域資源生態(tài)環(huán)境現(xiàn)狀、開發(fā)利用情況及相關區(qū)劃、規(guī)劃，分析圍填海活動對周邊功能區(qū)、利益相關者的影響，主要包括漁業(yè)、航運、旅游業(yè)和保護區(qū)等。

(5)圍填海工程規(guī)模。圍填?；顒幼鳛樗芯繂栴}系統(tǒng)的原始驅動因素，其規(guī)模(面積)大小、平面設計方案、所在海域水深條件等因素與可能造成的生態(tài)環(huán)境影響有著密切關系[13]。本研究中選取圍填海面積作為代表進入評價指標體系。

3.2 確定指標權重

評價指標體系選定后，構建圍填海適宜性的相互作用矩陣。將評價指標之間的影響程度分為5級，分別用無量綱數(shù)值0至4表示：0表示沒有影響；1表示微弱影響；2表示有較大影響；3表示有強烈的影響；4表示有極其強烈的影響。圍填海適宜性相互作用矩陣如表2所示。

表2中，各評價指標的意義如下：P1為功能區(qū)所在縣(市、區(qū))的經濟規(guī)模(GDP)；P2為功能區(qū)所在縣(市、區(qū))的人口數(shù)量；P3為底棲生物量；P4為海域海水中無機氮含量；P5為海域海水中活性磷酸鹽含量；P6為納潮量；P7為沖淤速率；P8為對漁業(yè)區(qū)的影響；P9為對航運區(qū)的影響；P10為對旅游區(qū)的影響；P11為對保護區(qū)的影響；P12為工業(yè)與城鎮(zhèn)建設區(qū)圍填海面積。

表2 圍填海適宜性相互作用矩陣

利用式(1)至式(3)計算得到各評價指標的權重值見表3。

表3 評價指標權重值計算結果

從表3中可以看出，圍填?；顒幼鳛橄到y(tǒng)的動因，與其他因素存在密切的聯(lián)系，其工程規(guī)模(面積)大小關系到影響程度的大小，所以權重值最高。保護區(qū)、漁業(yè)區(qū)和旅游區(qū)是對圍填?；顒臃浅Ｃ舾械墓δ軈^(qū)，應特別注意圍填海活動對它們可能造成的不良影響，所以權重值次之。社會經濟指標和生態(tài)環(huán)境指標處于適中的位置，且權重值相差不大。經濟社會的發(fā)展是圍填?；顒拥牡谝或寗恿?，而在取得經濟社會發(fā)展的同時，盡量減小圍填海活動對生態(tài)環(huán)境的不良影響，保護生態(tài)環(huán)境的健康也至關重要。由此可見，通過相互作用矩陣的確定的指標權重向量與實際經驗判斷具有較好的一致性。

3.3 評價指標量化賦分

依據(jù)對評價目標圍填海適宜性的影響程度，將所有評價指標分為3個等級，分別取值為3、2、1。如果指標對適宜性為正影響，分別取值3(影響較大)、2(影響一般)、1(影響較小)。如果指標對適宜性為負影響，分別取3(影響較小)、2(影響一般)、1(影響較大)。評價指標賦分情況見表4，各評價指標的物理意義如上文所述。

表4 圍填海適宜性評價指標賦分表

3.4 評價結果及分析

適宜性評價綜合指標計算采用加權模型為

(4)

式中：S為圍填海適宜性綜合評價指數(shù)；Bi為第i種評價因素的得分(無量綱)；Wi為第i種評價因素的權重；n為參加評價的因素數(shù)量。

由式(4)計算得到的適宜性綜合指標為：西店工業(yè)與城鎮(zhèn)建設區(qū)得分1.55，排名第3；西滬港底部工業(yè)與城鎮(zhèn)建設區(qū)得分1.52，排名第4；大港口工業(yè)與城鎮(zhèn)建設區(qū)得分2.78，排名第1；鄞州工業(yè)與城鎮(zhèn)建設區(qū)得分2.51，排名第2。4個評價對象圍填海適宜性從大到小的排序結果依次為：大港口工業(yè)與城鎮(zhèn)建設區(qū)、鄞州工業(yè)與城鎮(zhèn)建設區(qū)、西店工業(yè)與城鎮(zhèn)建設區(qū)、西滬港底工業(yè)與城鎮(zhèn)建設區(qū)。

大港口工業(yè)與城鎮(zhèn)建設區(qū)位于象山港口門附近海域，該海域水體交換較快，環(huán)境質量狀況較好，周邊海域無其他重要利益相關者及新設置的海洋功能區(qū)劃，所以該功能區(qū)的適宜性得分較高。但該圍填海項目的實施可能對象山港航道外段的流場、沖淤條件產生一定影響。鄞州工業(yè)與城鎮(zhèn)建設區(qū)同樣位于象山港口門處，但該圍填海的實施除考慮對象山港航道影響外，還需考慮對鄰近梅山旅游娛樂區(qū)的影響，所以該區(qū)域的得分較大港口工業(yè)與城鎮(zhèn)建設區(qū)低。西店工業(yè)與城鎮(zhèn)建設區(qū)位于象山港次級海灣鐵港底部，該區(qū)域生態(tài)環(huán)境脆弱，圍填海的實施對該次級海灣納潮量、水動力、沖淤變化有較大的影響，所以該功能區(qū)適宜性得分靠后。西滬港底工業(yè)與城鎮(zhèn)建設區(qū)位于次級海灣西滬港內，該圍填海的實施同樣對次級海灣納潮量、局部海域水動力條件變化有較明顯的影響。同時，該次級海灣內生物多樣性指數(shù)較高，灣內還設置了象山港海岸濕地海洋保護區(qū)，可能引起的生態(tài)環(huán)境負效應較大，所以該功能區(qū)的適宜性得分最低。

4 結束語

(1)巖石工程系統(tǒng)(RES)理論中的相互作用矩陣方法，能夠同時考慮系統(tǒng)中每個參數(shù)對系統(tǒng)的影響及參數(shù)之間的相互作用，理論上比較嚴密，用以計算圍填?；顒舆m宜性評價指標權重，具有一定的科學性和合理性。與層次分析法相比計算簡便，是復雜系統(tǒng)評價問題權重賦值的一種有效方法，可以在類似問題中推廣應用。

(2)圍填?；顒邮艿缴鐣洕l(fā)展的驅動，其實施會對海洋生態(tài)環(huán)境產生較大的影響，多種因素之間的相互耦合，因此，在評價指標體系構建時，需要同時考慮驅動因素、影響因素以及它們之間的相互作用。案例分析結果表明，基于巖石工程系統(tǒng)(RES)理論的相互作用矩陣計算評價指標權重，采用加權模型計算評價對象的適宜性具有合理性和可行性，評價結果與實際情況相符，可為其他區(qū)域圍填海適宜性評價提供參考。

[1] 劉偉， 劉百橋. 我國圍填?，F(xiàn)狀、問題及調控對策[J]. 廣州環(huán)境科學， 2008， 23(2): 26-30．

[2] 郭偉， 朱大奎. 深圳圍海造地對海洋環(huán)境影響的分析[J]. 南京大學學報: 自然科學版， 2005， 41(3): 286-296．

[3] 戴桂林， 蘭香. 基于海洋產業(yè)角度對圍填海開發(fā)影響的理論分析：以環(huán)渤海地區(qū)為例[J]. 海洋開發(fā)與管理， 2009， 26(7): 24-28．

[4] 于永海， 王延章， 張永華， 等. 圍填海適宜性評估方法研究[J]. 海洋通報， 2011， 30(1): 81-87．

[5] 王初升， 林榮澄， 黃發(fā)明. 珊瑚礁海岸圍填海適宜性的評估方法研究[J]. 海洋通報， 2012， 31(6): 695-699．

[6] 王初升， 黃發(fā)明， 于東升， 等. 紅樹林海岸圍填海適宜性的評估[J]. 亞熱帶資源與環(huán)境學報， 2010， 5(1): 62-67．

[7] 劉大海， 陳小英， 陳勇， 等. 海灣圍填海適宜性評估與示范研究[J]. 海岸工程， 2011， 30(3): 74-81．

[8] 王靜， 徐敏， 陳可鋒. 基于多目標決策模型的如東近岸淺灘適宜圍填規(guī)模研究[J]. 海洋工程， 2010， 28(1): 76-82．

[9] HUDSON J A.Rock mechanics principle in engineering practice. CIRIA Ground Engineering Report: Underground Construction， 1989．

[10] 吳紅剛， 馬惠民， 張紅利. 基于相互作用矩陣的山區(qū)公路展線路基穩(wěn)定性評價[J]. 巖土工程學報， 2011， 33(Suppl 1): 202-206．

[11] 李京梅， 孫晨， 謝恩年. 圍填海造地經濟驅動因素的實證分析[J]. 中國漁業(yè)經濟， 2013， 30(6): 61-68．

[12] 劉述錫， 馬玉艷， 卞正和. 圍填海生態(tài)環(huán)境效應評價方法研究[J]. 海洋通報， 2011， 30(6): 707-711．

[13] 馬軍. 大連圍填海工程對周邊海洋環(huán)境影響研究[D].大連：大連海事大學， 2009．

海洋公益性行業(yè)科研專項“東海沿岸狹長型海灣綜合整治集成技術及示范應用研究”(201105009)；國家海洋局第二海洋研究所基本科研業(yè)務費專項“基巖淤泥質港灣灘槽沖淤演變及其對人類活動響應的研究”(JT1203).

X171.1;U656.2

A

1005-9857(2015)03-0058-05