近年来, 习近平总书记多次在不同场合提及“海洋是高质量发展战略要地”的重要论述, 充分表明海洋经济高质量发展已经成为我国海洋工作的重要任务之一。然而, 在实践中, 我国海洋经济高质量发展仍面临许多问题, 包括资源协调不足、配套政策不完善、产能过剩、效益低下和环境污染等。因此, 通过研究海洋资源环境经济复合系统演化过程, 分析其演化机理, 进而探索我国海洋经济高质量发展、海洋生态文明建设未来方向。

1 研究概况

对于海洋复合系统相关研究, 国内研究成果比较丰富。在系统构建与评价方面, 高乐华等^[1]根据生态经济系统基本内涵, 展开系统研究; 李帅帅等^[2]研究海洋经济系统构成及其互动规律, 提出了蓝色经济空间拓展路径; 丁黎黎^[3]从“对象–理念–层次”三个维度阐述了海洋经济高质量发展的内涵, 构建海洋经济高质量发展评判体系。在系统协调与优化方面, 孙伯良等^[4]基于浙江省海洋经济、资源、环境现状, 测算其对应系统的协调水平; 盖美等^[5]借助三元协调发展模型对海洋资源环境经济复合系统承载力进行协调发展测度; 高强等^[6]基于协同学理论, 对海南省海洋生态经济系统协调度进行测算, 明确海南省海洋生态经济系统协调发展阶段。在系统动力与稳定方面, 狄乾斌等^[7]将系统动力学模型应用到海洋研究领域; 姜旭朝等^[8]从系统论的角度深化研究海洋经济, 从微观、中观和宏观层面分析不同要素的海洋经济系统动力演进机制; 王泽宇等^[9]采用综合评估模型测算了中国沿海11个省市的海洋经济系统稳定性指数, 引入障碍度诊断模型, 有针对性地提出中国海洋经济发展方向。

综上所述, 国内学者以协同学、系统动力学为研究理论基础, 对我国海洋资源、环境和经济发展进行了深入的探讨和研究。然而, 现阶段大部分研究成果较少涉及海洋复合系统的演化规律、演化路径等领域。本文从资源、环境和经济3个层面构建我国海洋复合系统, 运用熵权TOPSIS方法进行测算, 通过建立Logistic模型, 对我国海洋资源环境经济复合系统演化过程进行拟合分析和趋势预测。

2 海洋资源环境经济复合系统构建 2.1 海洋资源环境经济复合系统

海洋经济发展会对海洋资源和海洋环境带来影响, 海洋资源和海洋环境也会在一定程度上对海洋经济发展产生作用。图 1表示海洋资源环境经济复合系统概念模型, 从相互作用、相互制约关系来说, 海洋经济活动对海洋资源、环境的负面影响是复合系统的压力; 海洋资源、环境为海洋经济发展提供的支持是复合系统的承压; 而海洋资源和海洋环境的自身恢复、海洋经济的正向促进作用是复合系统的弹力。在压力–承压–弹力综合作用下, 复合系统不断进行自我调整, 保持系统稳定性。

深入研究海洋资源环境经济复合系统, 一方面可以探索海洋资源、环境和经济的相互作用、相互制约关系, 进而调控其协调发展模式, 对于影响系统稳定性的干扰因素, 采取合理科学措施, 提高整体协调水平, 实现海洋资源、环境和经济可持续发展; 而另一方面, 借助海洋资源环境经济复合系统能够反映海洋资源、环境和经济演化规律, 从而进行拟合分析和趋势预测。在此基础上, 通过对复合系统的结构进行调整和优化, 转变传统单一的发展方式, 实现海洋资源、环境和经济高质量发展。

2.2 评价指标的选取

在总结苟露峰等^[10]、鲁亚运等^[11]、赵玉杰^[12]等学者现有研究成果的基础上, 结合经济的实际情况和发展阶段, 遵循客观性、可行性和可比性的原则^[13], 构建评价指标体系, 如表 1所示。

在海洋资源方面, 不仅包括自然资源, 而且包括人为资源。选取人均海域面积、人均海洋捕捞量、人均海水养殖面积、人均沿海湿地面积、人均海洋渔业资源、人均海洋盐业资源和人均海洋矿业资源, 反映海洋自然资源情况; 选取港口码头长度、港口泊位个数和涉海就业人员指标, 反映海洋人为资源情况。

在海洋环境方面, 从环境污染、治理、保护和监测角度, 综合反映海洋环境情况。具体而言, 在环境污染上, 根据污染种类分为固体污染和水污染, 所以选取万元海洋生产总值工业废水排放量、万元海洋生产总值工业固体废物排放量等指标; 在环境治理上, 针对不同的污染以项目的形式推进污染的治理, 因此选取治理固体废物当年竣工项目、治理废水当年竣工项目等指标; 在环境保护上, 自然保护区作为环境保护的直接表现形式, 能够反映出环境保护的效果, 因此选取自然保护区总面积、自然保护区个数等指标; 在环境监测上, 水质达标率是海洋环境重要监测指标之一, 因此选取近岸海域功能区水质达标率、近岸海域一、二类海水比例、海洋工业废水排放达标率和海洋工业固体废物综合利用率等指标。

在海洋经济方面, 选取能够反映海洋经济规模情况相关指标, 包括人均海洋生产总值、人均海洋产业产值; 选取能够反映海洋经济变化情况相关指标, 包括海洋生产总值增长率、海洋产业产值增长率; 选取能够反映海洋经济结构情况相关指标, 包括海洋生产总值占地区生产总值比重、第三产业产值占海洋生产总值比重; 选取能够反映海洋经济开放情况相关指标, 包括地区沿海港口货物吞吐量、地区沿海港口旅客吞吐量、地区沿海港口集装箱吞吐量和地区沿海城市国际旅游收入。

2.3 数据的来源

本文以2006—2018年为研究时段, 以我国11个沿海省市为研究对象, 建立海洋资源环境经济复合系统。相关数据主要来源于《中国统计年鉴》《中国海洋统计年鉴》和11个沿海省市的地方统计年鉴、海洋环境状况公报等。个别数据经过相应的公式运算后得出。各指标的统计性描述结果, 如表 2所示。

3 基于熵权TOPSIS的海洋资源环境经济复合系统发展指数测算 3.1 基于熵权法的指标权重确定

熵值法按照各个指标之间数值离散程度, 计算对应权重, 从而进行客观赋权^[14]。本文所选指标均为各沿海省市海洋领域重要指标, 通过运用熵值法可以在实现客观评价的同时, 避免重要信息的遗漏。因此, 本文基于熵值法来确定海洋资源环境经济复合系统评价指标体系中不同评价指标的权重。

1) 标准化矩阵构建

X_ij表示第i项指标的第j年的初始值, i=1, 2, 3, …, n, n表示评价指标数; j=1, 2, 3, …, m, m表示评价年份数。

$ \mathit{\boldsymbol{X}} = \left[ {\begin{array}{*{20}{c}} {{X_{11}}}&{{X_{12}}}&.&.&{{X_{1m}}}\\ {{X_{21}}}&{{X_{22}}}&.&.&{{X_{2m}}}\\ .&.&.&.&.\\ .&.&.&.&.\\ {{X_{n1}}}&{{X_{n2}}}&.&.&{{X_{nm}}} \end{array}} \right]. $

(1)

各个被选取指标有不同量纲, 难以简单比较, 需要分为正向、逆向指标^[15]; 然后, 再运用极值变化法进行标准化处理。具体公式如下:

$ 正向指标: {X'_{ij}} = \left( {{X_{ij}} - {X_{i{\rm{min}}}}} \right)/\left( {{X_{i{\rm{max}}}} - {X_i}_{{\rm{min}}}} \right), $

(2)

$ 负向指标: {X'_{ij}} = \left( {{X_{i{\rm{max}}}} - {X_{ij}}} \right)/\left( {{X_{i{\rm{max}}}} - {X_{i{\rm{min}}}}} \right). $

(3)

其中, X_imax表示第i个指标中所有年份的最大值; X_imin表示在第i个指标中所有年份的最小值, ${X'_{ij}}$表示X_ij经过标准化处理得到的数据。同时, 由于标准化处理后的${X'_{ij}}$可能为零, 而对数的底数应为非负数, 标准化结果需要进一步处理:

$ {R_{ij}} = {X'_{ij}} + {\rm{0}}{\rm{.01}}{\rm{.}} $

(4)

根据标准化的结果, 构建标准化矩阵R。

$ \mathit{\boldsymbol{R}} = \left[ {\begin{array}{*{20}{c}} {{R_{{\rm{11}}}}}&{{R_{{\rm{12}}}}}&.&.&{{R_{{\rm{1}}m}}}\\ {{R_{{\rm{21}}}}}&{{R_{{\rm{22}}}}}&.&.&{{R_{{\rm{2}}m}}}\\ .&.&.&.&.\\ .&.&.&.&.\\ {{R_{n{\rm{1}}}}}&{{R_{n{\rm{2}}}}}&.&.&{{R_{nm}}} \end{array}} \right]. $

(5)

2) 熵值法赋权

计算第i项指标在第j年的样本值比重:

$ {p_{ij}} = \frac{{{R_{ij}}}}{{\sum\limits_{j = 1}^m {{R_{ij}}} }}. $

(6)

计算第i项指标的熵值e_i: m是年数, n是指标总数

$ {e_i} = - k\sum\limits_{j = 1}^m {{p_{ij}}\ln {p_{ij}}} . $

(7)

其中, e_i > 0, k > 0, 假设R_ij全部相等, 令$k = \frac{1}{{\ln m}}$, 则此时的e_i为极大值:

$ {e_i} = - \frac{1}{{\ln m}}\sum\limits_{j = 1}^m {{p_{ij}}\ln {p_{ij}}} . $

(8)

计算信息权重w_i:

$ {w_i} = \frac{{1 - {e_i}}}{{\sum\limits_{i = 1}^n {1 - {e_i}} }}. $

(9)

3) 多层次评价系统的评价

由于熵具有可加性特征, 因此通过计算下层结构指标信息效用值, 得出上层结构权重。用H_k(k=1, 2, 3, …, k)表示指标效用值, 用H表示全部指标效用值总和。

$ H = \sum\limits_{k = 1}^n {{H_k}} . $

(10)

则相应类指数的权重为:

$ w = {H_k}/H. $

(11)

4) 基于熵值的评价矩阵构建

$ \begin{array}{l} {\boldsymbol{Z}} = \left[ {\begin{array}{*{20}{c}} {{z_{{\rm{11}}}}}&{{z_{{\rm{12}}}}}&.&.&{{z_{{\rm{1}}m}}}\\ {{z_{{\rm{21}}}}}&{{z_{{\rm{22}}}}}&.&.&{{z_{{\rm{2}}m}}}\\ .&.&.&.&.\\ .&.&.&.&.\\ {{z_{n{\rm{1}}}}}&{{z_{n{\rm{2}}}}}&.&.&{{z_{nm}}} \end{array}} \right]\\ = \left[ {\begin{array}{*{20}{c}} {{R_{{\rm{11}}}}.{{\rm{w}}_{\rm{1}}}}&{{R_{{\rm{12}}}}.{{\rm{w}}_{\rm{1}}}}&.&.&{{R_{{\rm{1}}m}}.{{\rm{w}}_{\rm{1}}}}\\ {{R_{{\rm{21}}}}.{{\rm{w}}_{\rm{2}}}}&{{R_{{\rm{22}}}}.{{\rm{w}}_{\rm{2}}}}&.&.&{{R_{{\rm{2}}m}}.{{\rm{w}}_{\rm{2}}}}\\ .&.&.&.&.\\ .&.&.&.&.\\ {{R_{n{\rm{1}}}}.{{\rm{w}}_n}}&{{R_{n{\rm{2}}}}.{{\rm{w}}_n}}&.&.&{{R_{nm}}.{{\rm{w}}_n}} \end{array}} \right]. \end{array} $

(12)

经过熵值法客观赋权后, 反映出各个评价指标之间不同的重要程度^[16], 得到不同指标权重下的海洋资源环境经济复合系统评价矩阵。

3.2 基于TOPSIS法的系统评价

TOPSIS模型适用于对多项指标、多个方案进行选择^[17]。通过运用TOPSIS模型, 能够得到较为客观的评价。

1) 正负理想解确定

$z_i^ + $是第i个指标在j年期间的最大值, 将$z_i^ + $设定为正最理想解; $z_i^ - $是第i个指标在j年期间的最小值, 将$z_i^ - $设定为负最理想解, 对应具体公式为:

$ {z^ + } = \left\{ {\mathop {\max }\limits_{1 \le i \le n} {\rm{ }}{z_{ij}}|i = 1, 2, \cdots , n} \right\} = \left\{ {z_1^ + , z_2^ + , \cdots , z_n^ + } \right\}. $

(13)

$ {z^ - } = \left\{ {\mathop {\min }\limits_{1 \le i \le n} {\rm{ }}{z_{ij}}|i = 1, 2, \cdots , n} \right\} = \left\{ {z_1^ - , z_2^ - , \cdots , z_n^ - } \right\}. $

(14)

2) TOPSIS确定指标到正、负理想值之间的距离

采用欧氏计算法, 确定指标到正、负理想值的距离, 令$D_j^ + $为第i个指标与$z_i^ + $的距离, 令$D_j^ - $为第i个指标与$z_i^ - $的距离, 计算公式为:

到正理想解之间的距离:

$ D_j^ + = \sqrt {\sum\limits_{i = 1}^n {{{(z_i^ + - {z_{ij}})}^2}} } . $

(15)

到负理想解之间的距离:

$ D_j^ - = \sqrt {\sum\limits_{i = 1}^n {{{\left( {z_i^ - - {z_{ij}}} \right)}^2}} } . $

(16)

3) 计算综合评价指数

M_j代表第j年综合评价指数, 取值区间: (0, 1]。指数数值越接近1表示综合评价的得分越高, 而越接近0则表示综合评价得分越低。评价指数的计算公式如下所示:

$ {M_j} = \frac{{D_j^ - }}{{D_j^ + + D_j^ - }}. $

(17)

3.3 海洋资源环境经济复合系统综合指数测算

由于各个系统的综合指数对等, 可以进一步计算出海洋资源环境经济复合系统综合指数, 具体情况见表 3。

从上述海洋资源环境经济复合系统及其子系统综合指数的变化情况可知, 2006—2018年期间, 我国海洋资源、经济子系统综合指数总体水平较高, 呈现明显上升趋势, 这表明由于我国各项海洋活动、事业快速推进, 海洋资源、经济发展取得显著成效。相比之下, 海洋环境子系统综合指数总体水平较低, 而且存在波动情况, 这表明我国海洋环境污染问题较为严重, 海洋生态调节功能弱化。因此, 在3个子系统相互作用下, 我国海洋资源环境经济复合系统综合指数不高, 状态较为稳定, 仍有一定的提升空间。

4 海洋资源环境经济复合系统演化实证研究 4.1 海洋资源环境经济复合系统的Logistic模型

Logistic模型作为一种重要分析工具, 能够描述一般系统发展演化过程^[18]。依据Verhulst^[19]所构建的Logistic增长模型, 参照陈海波等^[20]、周凌云等^[21]、周韬^[22]等学者的研究成果, 海洋资源环境经济复合系统演化路径数学模型为:

$ \left\{ \begin{array}{l} \frac{{{\rm{d}}X}}{{{\rm{d}}t}} = kX \times \left( {{\rm{1}} - \frac{X}{N}} \right)\\ {X_{{\rm{(0)}}}} = {X_{\rm{0}}} \end{array} \right., $

(18)

其中, dX/dt表示海洋资源环境经济复合系统在t时刻下发展速度; k表示为复合系统内在增长率, 与内部耦合度具有关联性; N为海洋资源经济复合系统最大承载值, 与地区海洋经济、资源水平存在一定关系; 每个个体平均所占有资源量为1/N, X/N为消耗总资源; [1–X/N]为剩余资源, 代表Logistic系数; X₀为初始时刻系统总量。

对其系数分析可知: 若海洋资源环境经济复合系统总量趋近于0, 则[1–X/N]就接近1, 表明复合系统中资源尚未被利用, 系统演化趋势呈指数增长; 若复合系统总量X趋近于N, 则[1–X/N]就接近0, 表明复合系统中资源被充分利用, 此时复合系统演化趋势呈饱和状态; 当复合系统总量由0逐渐上升到N时, [1–X/N]由1逐渐下降到0, 表明剩余资源逐渐变小。

对公式(18)进行求解, 可得:

$ X = \frac{N}{{{\rm{1}} + C \times {\rm{exp}}\left( { - kNt} \right)}}. $

(19)

其中, C为常数, 数值随系统演化阶段的变化而变化。

设X(0)=α为初始状态, 其中0 < α < M, 则:

$ X = \frac{N}{{{\rm{1}} + \frac{N}{{{\rm{ { α} }} - {\rm{1}}}}{\rm{exp}}\left( { - kNt} \right)}}. $

(20)

公式(18)表示海洋资源环境经济复合系统在任一时刻的增长速度, 可称之为复合系统成长速度方程。公式(19)表示海洋资源环境经济复合系统演化动态变化轨迹, 是复合系统状态演化方程。

对公式(18)求导, 可得:

$ \frac{{{{\rm{d}}^2}X}}{{{\rm{d}}{t^2}}} = {k^2}X\left( {1 - \frac{X}{N}} \right)\left( {1 - \frac{{2X}}{N}} \right). $

(21)

公式(21)用来表示海洋资源环境经济复合系统任何一时刻加速度。令d²X/dt²=0, 可以得到复合系统状态演化曲线拐点: X₁=0, X₂=N/2和X₃=N。由于0 < X < N, 所以拐点为X=N/2。此时, 复合系统成长速度曲线达到最大值kN/4。

对公式(21)求导, 可得:

$ \frac{{{{\rm{d}}^3}X}}{{{\rm{d}}{t^3}}} = {k^3}X\left( {1 - \frac{X}{N}} \right)\left( {\frac{{6{X^2}}}{{{N^2}}} - \frac{{6X}}{N} + 1} \right). $

(22)

令d³X/dt³=0, 得到拐点$\left({3 - \sqrt 3 } \right)N/6$和$\left({3 + \sqrt 3 } \right)N/6$。此时, 复合系统成长速度曲线数值为kN/6。经过推导, 复合系统成长速度曲线和状态演化曲线的特征, 如图 2所示。

由图 2可得, 海洋资源环境经济复合系统随着时间以S型曲线增长, 上界渐进线X=N(t, ∞), 根据增长情况, 可以将演化过程可以分为四个阶段, 其具体演化特征如表 4所示。

1) 第一阶段(0 < t≤t₁): 海洋资源环境经济复合系统限制较小, 虽然总体发展程度较低, 但发展速度逐步提升。在这一阶段, 复合系统成长量呈现指数增长。在t₁处, 速度为kN/6, 此时加速度达到最大值, 影响海洋资源环境经济复合系统的各种因素累计最大, 系统成长量达到$\left({3 - \sqrt 3 } \right)N/6$, 系统处于发展的初步阶段。

2) 第二阶段(t₁ < t≤t₂): 海洋资源、环境、经济规模增加迅速, 复合系统活力增强、发展空间大。在这一阶段, 复合系统发展速度递增, 加速度减少, 属于准线性增长。在t₂处, 速度为kN/4, 成长量为N/2, 成长速度达到最大, 系统处于发展的成长阶段。

3) 第三阶段(t₂ < t≤t₃): 海洋资源环境经济复合系统增长动力减弱, 系统受到资源、空间制约。在这一阶段, 复合系统成长速度和加速度均减小, 但仍属于准线性增长。在t₃处, 速度为kN/6, 成长量为$\left({3 + \sqrt 3 } \right)N{\rm{/6}}$, 系统处于发展的成熟阶段。

4) 第四阶段(t₃ < t < +∞): 海洋资源环境经济复合系统趋于平稳, 接近市场需求最大值。在这一阶段, 复合系统成长速度递减, 加速度递增。在t₃处, 速度小于kN/6, 成长量大于$(3 + \sqrt 3)N{\rm{/6}}$, 逐渐达到极限值N, 系统处于发展的衰退阶段。

4.2 海洋资源环境经济复合系统演化方程分析

为了进行Logistic法则参数估计, 将公式(20)进一步转换成:

$ E = a/\left[ {1 + b \exp\left( { - ct} \right)} \right]. $

(23)

E代表我国t年的海洋资源环境经济系统演化程度, a、b、c为待估参数: a为复合系统演化度能够达到的最大值; b为积分常数; c为复合系统的增长率。通过三段和值法确定参数a、b和c的初始值, 使用计量经济软件Origin进行拟合和参数估计, 计算结果如表 5。

根据表 4中的Logistic法则参数估计结果, 得到海洋资源环境经济复合系统及其子系统的估计方程。

海洋资源子系统:

$ X = \frac{{0.709\;7}}{{1 + 0.936\;4 \times \exp\left[ { - 0.289\;3 \times \left( {t - 2\;005} \right)} \right]}}. $

(24)

海洋环境子系统:

$ X = \frac{{0.632\;9}}{{1 + 0.785\;3 \times \exp\left[ { - 1.163\;0 \times \left( {t - 2\;005} \right)} \right]}}. $

(25)

海洋经济子系统:

$ X = \frac{{0.632\;9}}{{1 + 0.785\;3 \times \exp\left[ { - 1.163\;0 \times \left( {t - 2\;005} \right)} \right]}}. $

(26)

海洋资源环境经济复合系统:

$ X = \frac{{0.647\;8}}{{1 + 0.512\;0 \times \exp\left[ { - 0.391\;2 \times \left( {t - 2\;005} \right)} \right]}}. $

(27)

由参数估计结果可知, 系统拟合度R分别为: 0.874 4、0.848 6、0.844 9、0.892 3, 参数检验比较显著, 拟合度比较理想, 分别绘制各系统的Logistic演化发展曲线。

在图 3中, 实线表示2006—2018年海洋资源环境经济复合系统及其子系统演化实际曲线, 反映出我国海洋资源、环境、经济实际演化过程; 而虚线则表示拟合曲线, 反映出海洋资源环境经济复合系统及其子系统的Logistic方程拟合结果。演化实际曲线、拟合曲线基本吻合, 这表明海洋资源环境经济复合系统及其子系统的演化规律符合Logistic法则。

通过比较可以发现, 2011—2013年期间, 海洋资源子系统中的演化实际曲线明显高于演化拟合曲线, 这可能是我国2011年制定的《全国海洋功能区划》有效地指导沿海地区海洋资源开发活动。我国11个沿海省市陆续制定和公布地方海洋功能区划, 明确海洋资源开发利用的相关法律法规、管理制度、技术手段和监督评价机制, 海洋资源利用效率有所提高。因而, 在此期间我国海洋资源子系统的演化实际曲线高于拟合曲线。然而, 受制于传统治理观念、治理手段单一等因素的影响, 海洋功能区划的措施实施成效有待提高, 同时我国沿海地区海洋渔业过度捕捞问题越来越严峻, 大型海洋钻井平台数量不断增加。所以, 2014—2017年海洋资源子系统的演化实际曲线低于拟合曲线。2018年, 我国沿海省市加强对海洋生态红线监管力度, 有效引导海洋资源开发活动, 海洋资源子系统演化实际曲线呈现上升趋势。总体上看, 海洋资源开发利用仍有一定改善空间。

从2008年开始, 海洋环境子系统中的演化实际曲线明显波动, 这可能是由于赤潮等海洋环境灾害增多, 影响我国沿海海域生态健康。近年来, 我国近海海洋环境灾害呈现出增加趋势, 海上重大污染事故频发, 近海灾害性风暴潮偏多。2008—2013年我国近岸海域四类和劣四类海水比例从18.3%上升到25.6%, 海洋环境灾害造成部分沿海海域海水富营养化程度较为严重, 降低我国沿海海域水体质量。除此之外, 海洋污染过度排放等问题日益突显, 也造成了系统演化实际曲线的波动。2008年直排海污染源污水排放总量为45.65亿吨; 2013年为63.84亿吨, 上升幅度较大。海洋污染物过度排放超过海洋环境承载力, 造成海洋生物大量死亡, 降低海洋环境生态修复能力。从整体情况上看, 海洋环境问题不断凸显。

2014—2016年, 海洋经济子系统的演化实际曲线明显低于演化拟合曲线, 这可能是我国海洋经济进入新常态, 增速有所放缓, 海洋产业转型升级紧迫, 海洋经济发展面临资源、劳动力、资金、环境等方面的压力。我国海洋产业主体仍是传统海洋产业, 海洋战略性产业发展规模还比较小, 产值较低。过于依赖传统海洋产业的“粗放”型海洋经济发展模式, 导致海洋经济子系统演化实际曲线偏低。2016年, 国家战略性新兴产业规划出台, 明确海水淡化、海洋生物医药、海洋工程装备等海洋战略性新兴产业为今后海洋经济重点发展方向; 各沿海省市如广东省、福建省、上海市等也提出涉及海洋战略性新兴产业发展的相关规划和实施方案。2017—2018年, 由于海洋战略性新兴产业的发展为我国海洋经济提供增长空间, 海洋经济子系统演化实际曲线上升。

海洋资源环境经济复合系统的演化实际曲线在2006—2016年期间, 总体低于演化拟合曲线。在此期间, 我国海洋经济的发展还处在初级阶段, 存在着海洋资源消耗较大和海洋环境污染严重等问题, 海洋生态环境承载能力整体较差。从2016年开始, 海洋资源环境经济复合系统中的演化实际曲线呈现出缓慢上升的趋势, 这可能是由于进入“十三五”时期, 各沿海省市合理配置海域资源、大力推进海洋生态文明建设, 加强完善近岸海域使用和管理, 通过建设“海洋生态文明示范区”和打造“美丽海湾”等一系列的措施和方法, 加大对海洋环境的保护和治理, 改善海洋环境质量, 有效促进海洋资源、环境、经济统筹协调发展。

4.3 海洋资源环境经济复合系统演化预测分析

根据海洋资源环境经济复合系统及其子系统的演化方程, 得出各系统到2050年的演化预测曲线, 如图 4所示。

由海洋资源子系统演化预测曲线可知, 我国传统海洋资源开发已经进入成熟期, 总体发展速度缓慢; 在2020—2050年期间, 传统海洋资源开发利用将会达到饱和状态。目前, 我国各科研院所正积极对“可燃冰”等海洋蓝色能源的开采和利用展开技术攻关。2017年, 我国已经掌握海洋“可燃冰”固态流化开采技术。但是, 现阶段的技术水平难以达到大量、高效、安全开采“可燃冰”的目标。随着未来技术发展和进步, 我国海洋资源开发将会从传统海洋资源转向蓝色海洋资源, “可燃冰”等海洋蓝色能源产业发展将会上升到国家战略高度, 最终实现工业化开采, 形成由浅海至深海的多层次立体海洋资源开发格局。

由海洋环境子系统演化预测曲线可知, 从2010年开始我国海洋环境子系统发展速度基本停滞, 这反映出我国海洋环境污染问题十分严峻、海洋生态调节功能弱化的现实情况; 在2020—2050年期间, 我国海洋环境子系统仍会处于成熟阶段。近些年来, 中央和地方政府积极加大海洋环境治理力度, 通过制定海洋功能区划、划定海洋生态红线等一系列措施来改善海洋环境, 但海洋生态环境承载能力整体仍然较差。随着政策工具和手段的丰富和发展, 中央和地方将会继续出台海洋环境治理相关措施制度, 以海洋环境治理修复工程为抓手, 借助海洋环境税费、污染控制指标转让等市场激励型海洋环境规制工具, 显著提升海洋环境治理成效。

由海洋经济子系统演化预测曲线可知, 在新常态背景下, 我国海洋经济发展增速放缓, 海洋经济子系统也处于成熟阶段; 在2020—2050年期间, 我国传统海洋产业的发展已经达到饱和状态。“十二五”期间, 我国已经基本形成海洋新兴产业体系, 海洋战略性新兴产业链条不断延伸、产业规模日益扩大。在未来一段时间内, 我国海洋生物医药、海洋电力、滨海旅游等海洋战略性新兴产业产值将会大幅度提高, 同比增速将会维持在较高水平。借助海洋产业转型升级和现代海洋产业体系构建, 依靠创新驱动、技术引领, 我国海洋经济将会迎来高质量发展时期。

从整体上看, 在资源、环境、空间等基本条件不改变的前提下, 我国海洋资源环境经济复合系统演化已经达到基本饱和。然而, 随着我国陆海统筹、区域协调发展等国家战略的实施, 陆海联系、区域联动将会得到加强。未来, 为打破各行政区发展经济带的局限, 将从国家层面统筹北部、东部、南部沿海经济带建设, 形成自北向南连续、统一的海洋经济带发展新格局。随着现代化沿海经济带建设的推进, 我国海洋资源、环境和经济的相互作用、协调耦合将会得到加强和改善, 海洋资源环境经济复合系统的演化也会从较低层次向较高层次跃进。

5 结论与建议 5.1 结论

从资源、环境和经济三个层面构建我国海洋复合系统, 运用熵权TOPSIS方法进行测算, 通过建立Logistic模型, 对我国海洋资源环境经济复合系统发展过程进行拟合分析和趋势预测。得到相关研究结论: 1) 海洋资源环境经济复合系统是基于海洋资源、环境、经济相互作用、相互制约而形成的, 其系统演化过程受到经济增长机制和生态平衡机制的影响。2) 海洋资源环境经济复合系统及其子系统的演化规律符合Logistic法则, 海洋资源子系统最大演化度高于海洋环境、经济子系统。3) 海洋资源环境经济复合系统及海洋资源、经济子系统演化度的增长率较大, 正处于成熟阶段, 而海洋环境子系统演化度达到饱和状态, 正处于衰退阶段。

5.2 建议

基于上述研究结论, 围绕均衡海洋资源开发利用、优化海洋空间规划体系、打造海洋特色产业集群和建立海洋互动协作机制这四个方面, 提出以下措施建议:

从近海到深海, 均衡海洋资源开发利用。各沿海省市应在开展海洋活动、推进海洋事业过程中, 坚持合法、合理、科学原则, 加大海岸、海岛等近海资源的开发利用, 将海岸生态修复工程与滨海旅游产业相结合。海岸修复项目应该要结合海岸自然特征, 以海岸空间养护、海岸自然景观恢复等为主, 在修复生态环境的同时, 提升旅游休闲功能。在深海资源方面, 应该加大海洋资源开发技术研发投入力度, 从单项开发转向立体开发、综合开发, 丰富海洋资源利用手段、途径; 从近海走向深海, 通过借助深海勘探装备和工具, 积极探索和开发各类深海资源。

积极探索“多规合一”, 优化海洋空间规划体系。政府部门应该扮演“多规合一”的主导角色, 制定、落实海洋区划配套支持制度, 依靠政策工具和手段, 积极对各类社会力量的涉海活动进行引导。“多规合一”的探索应该以省级规划试点展开, 这样既能够将国家宏观规划具体化, 又能够指导监督基层规划, 将各类海洋规划向基层推进。优化海洋空间规划体系, 应该理清海洋空间规划层级、主体之间关系, 应该坚持海洋生态红线不动摇, 通过上位规划指导下位规划, 构建一个“多规合一”的海洋空间规划体系。

创新驱动海洋经济发展, 打造海洋特色产业集群。各涉海主体应该通过推动海洋科技研发、体制机制、管理制度等方面的有效创新, 为海洋经济注入新的活力。涉海部门、涉海企业应该加大科技研发投入, 针对战略要求和市场需求, 建设和发展一批重点海洋高新技术产业, 打造海洋特色产业集群。一方面, 应该加快海洋经济三大产业协同发展, 打造国际一流海洋服务业基地; 另一方面, 应该根据地方资源禀赋情况, 打造海洋装备、海洋生物制药、滨海旅游等特色产业集群, 形成较为完备的现代海洋产业组织和分工体系。

推动成立地方海洋委员会, 建立海洋互动协作机制。各沿海省市应该成立地方海洋委员会, 从而有效发挥协调和监管职能, 实现区域海洋事务宏观管理, 加强对区域海洋统筹发展的总体协调和战略指导。为了进一步建立海洋互动协作机制, 在基础设施建设、生态环境治理等过程中, 应该建立跨部门、跨省市的行政管理体制; 在人才、资金等资源配置上, 应该建立和完善区域内外交流网络, 以推动区域间的交流与合作。各级地方海洋委员会应该通过相应的政策激励与约束手段, 构建海洋开发深度合作平台, 建立海洋互动协作机制, 实现海洋持续健康发展。