21世纪是海洋的世纪，建设海洋工程、发展海洋经济、开展海洋研究是推动强国战略的重要内容。目前，我国海洋事业正逐步向深远海迈进，在资源开发、工程建设和捍卫海洋主权等方面已取得一定成果^[1-3]。海底稳定性是海洋开发建设的重要制约因素，海底地质稳定性差，制约程度高，开发风险高；海底地质稳定性好，制约程度低，开发风险低^[4]。如何准确对深远海海底稳定性进行评价是开展深远海开发利用的重点工作之一，由于影响海底稳定性的因素具有复杂性和不确定性，因此，评价结果并不具有绝对含义^[5]。

目前，涉及海底稳定性的研究成果多集中于近海且常以单一地貌为研究对象，缺少对深远海海底稳定性的宏观性评价。国外学者侧重于海底边坡、斜坡等单一地貌，常采用概率、数值模拟、半解析等数学方法进行分析^[6-7]，对黑海高加索大陆架^[8]、高纬度冰期斯瓦尔巴-巴伦支海^[9]等区域的海底斜坡稳定性进行了研究，探索了地震和沉积速率对海底斜坡稳定性的影响^[10]，海底沉积斜坡和增生楔的稳定性^[11]等问题，成果大多只聚焦于目标区的稳定性分析，缺乏定量化的分区评价。我国学者除对海底斜坡进行了大量研究外^[12]，近年来以地貌类型、地形坡度、地震数据、工程地质参数等因素作为主要评价因子，利用GIS、层次分析和模糊综合评价等方法进行海底稳定性评价和分区，在南海北部陆坡^[13]、南堡-曹妃甸海域^[14]和莱州湾^[15]等区域进行了小范围的宏观评价。然而，这些区域大多位于近海，深远海海底稳定性评价研究还存在较多空白。相较于近海，深远海地质环境复杂，具有勘查难度大、数据精度低和研究成果少等问题，如何基于现有数据资料正确、全面地评价深远海海底稳定性，对于深远海资源开发和工程建设具有深远意义。

海底微地貌是影响深远海海底稳定性的重要因素，不仅其本身的起伏变化会对区域稳定性产生影响，而且常为海底滑坡、陡坎等致灾因子的发育提供条件。基于研究团队前期工作基础，本文选取西太平洋菲律宾海中南部某区域为研究区，利用ArcGIS平台建立研究区DEM（数字高程模型，Digital Elevation Model），提取宏、微观地貌因子（坡度、地形起伏度），结合海底地震、海底底质类型、海底地质灾害等要素，通过模糊数学评价方法将以上要素数字化，综合分析研究区海底稳定性并进行分区，提出了基于微地貌特征的深远海海底稳定性评价方法。研究成果能够为深远海海底稳定性评价提供技术方法参考，丰富深远海海洋地质研究。

1.   研究区概况

研究区地处欧亚板块、印度-澳大利亚和太平洋板块交汇处，位于西太平洋菲律宾海中南部，横跨西菲律宾海盆和帕里西维拉海盆，九州-帕劳海岭从中间穿过，范围为14.5°～17.0°N、132.5°～136.0°E（图1），面积约为1.0×10⁵ km²，海底标高为−6712.8～−1692.2 m。

图  1  菲律宾海板块海底地形地貌与研究区位置图^[22]

底图数据来源于全球水深数据。

Figure  1.  Seabed topography and geomorphology of the Philippine Sea plate and location of the study area ^[22]

DEM are derived from global bathymetric data.

下载: 全尺寸图片 幻灯片

菲律宾海四周被呈链状的岛弧、海沟包围，中部宽阔，南北窄小，近菱形展布，具有典型的沟-弧-盆体系特征^[16]。现有的研究表明，菲律宾海板块自61 Ma开始扩张，期间经历多次旋转、扩张，其形成与演化与太平洋板块运动息息相关^[17]。在43 Ma左右，受太平洋板块俯冲影响形成古伊豆-小笠原岛弧^[11]。其后，因四国-帕里西维拉海盆扩张影响，古伊豆-小笠原岛弧裂解为九州-帕劳海岭和伊豆-小笠原-马里亚纳弧^[18-21]。

研究区受复杂而剧烈的地质作用影响，微地貌单元广泛分布，发育有海山、海丘、山间谷地等多种微地貌单元，其中海山主要发育在九州-帕劳海岭及其西部海盆，海丘、山间谷地则集中分布在海岭东侧的帕里西维拉海盆。研究区能够较为明显地反映出菲律宾海海底地貌特征，利于对其海底稳定性进行研究，具有一定的代表性。

2.   数据来源和分析方法

2.1   数据来源

本文数据主要来源于公共平台公开的海底高程数据，主要有SRTM、ASTER等。利用ArcGIS平台建立研究区海底DEM（DEM即数字高程模型，是通过有限的地形高程数据实现对地面地形的数字化模拟，是用一组有序数值阵列形式表示地面高程的一种实体地面模型），从微观和宏观两个角度提取关键地貌因子；基于全球地震数据（https://earthquake.usgs.gov/earthquakes/search/），运用Kriging插值法对研究区地震区划进行分析；海洋地质灾害分布和底质类型数据来源于对公开数据的处理和识别。

2.2   典型地貌因子计算

坡度（微观地貌因子）和地形起伏度（宏观地貌因子）作为典型的地貌因子，能够代表性地反映海底微地貌特征，参考已有研究成果，本文选择坡度和地形起伏度两个地貌因子进行海底稳定性评价^[23-24]。将研究区DEM在ArcGIS平台中以“标准差”的形式进行分类，根据表1中的计算方式，获取研究区坡度和地形起伏度分布特征。

表  1  典型地貌因子概念及算法^[25-26]

Table  1.  Concept and algorithm of typical geomorphic factor ^[25-26]

地形因子	概念	公式
坡度	描述地表单元陡缓程度，表示地面在某点的倾斜度	tanα=△h/△dα—坡度，△h—高程差（m），△d—水平距离（m）
地形起伏度	区域最高点与最低点海拔高度的差值，反映地形起伏特征[25]	RAi=Zimax—Zimin，RAi—地形起伏度（m），Zimax—区域内最大高程值（m），Zimin—区域内最小高程值（m）

下载: 导出CSV 

| 显示表格

2.3   模糊数学评价

由于影响海底稳定性的因素复杂，稳定性分级界限模糊，加之人们主观上的差异，不易对海底稳定性得出一致评价。因此，综合前人的经验和相关的评价结果，采用模糊数学方法来进行海底稳定性评价工作是可行的。

模糊数学评价是地质评价中常用的一种数学评价方法，可以实现将相对复杂且不够确定的影响因素数学化。研究区位于西太平洋菲律宾海中南部，受制于地理位置和技术条件，很难对其海底稳定性进行精确评价，模糊数学评价法可以利用模糊数学的隶属度和隶属函数理论把定性评价转化为定量评价，从而对研究对象做出一个综合的评价^[27]。通过数字化评价指标，综合评定评价指标等级，可以有效减少评价主体间的主观评定误差（图2）。

图  2  评价方法的技术路线

Figure  2.  The technical flowchart for evaluation methods

下载: 全尺寸图片 幻灯片

3.   微地貌特征

3.1   微观地貌因子——坡度

研究区坡度变化范围为0°～58°，平均坡度6.9°，标准差为7.2°，坡度主要集中于0°～15°，表明研究区坡度变化较大，但大部分区域坡度较为平缓。坡度的低值主要集中在研究区的中北部，该区域地形较为平坦，坡度变化较小；研究区东部和西南部地形相对复杂，海山、海丘等地貌单元发育，坡度较高且变化频繁；坡度最大值分布在九州-帕劳海岭，该区域具有较高的高程差，地形起伏大，坡度较高。

3.2   宏观地貌因子——地形起伏度

研究区地形起伏度相对较小，主要集中在44.6 m以下，平均值为31.9 m，最大值仅为344.6 m。低值主要集中在研究区中北部，该区域地形较为平坦，地貌单元类型相对单一；高值集中在九州-帕劳海岭和海山等区域。

3.3   海底底质类型

研究区内底质类型丰富，除常见的远洋黏土外，区域内广泛发育了基岩、铁锰结核和硅藻软泥。基岩硬度相对较高，分布的区域通常具有较大的土体承载力；硅藻软泥则相反，其相对较软，对土体承载力会产生完全不同的影响。

4.   潜在地质灾害分布特征

4.1   海底滑坡

海底滑坡在主、被动大陆边缘广泛发育，可以将失稳的沉积物带到深海盆地^[28-33]。九州-帕劳海岭边缘地带不仅有较高的坡度，而且具有很大的高程差，极易发生海底滑坡。研究区中北部地形较为平坦，西部海盆区虽然坡度也比较大，但是高程差相对较低，地貌起伏变化不大，为海底滑坡的低易发区。

4.2   陡坎

陡坎是研究区内一种常见且分布广泛的致灾因子，是指与海底斜坡地形走向平行或近于平行且坡度较大的海底陡坡^[34]。受九州-帕劳海岭和大型地貌单元影响，陡坎主要发育在研究区中部和大规模海山、海丘群处，这些区域都具有比较高的高程差和坡度。

4.3   崩塌

崩塌是海底较陡斜坡上的岩土体突然受力而脱离崩落的地质现象，具有较大的破坏性。研究区地形起伏较大，陡坡广泛发育，为崩塌的产生提供了有利条件。区域内大规模海山、山间盆地以及九州-帕劳海岭附近往往都具有较大的地形起伏度和坡度，是崩塌的易发区域。

5.   海底稳定性评价

本文采用模糊数学原理，建立由多种影响因素组成的评价因素集A，构建与影响因素集相匹配的模糊评语集X。通过模糊信息化确定各单一因素对评审等级的归属程度，建立评价因素集A到评语集X的模糊评价矩阵R。利用模糊权值矩阵计算各因素在评价目标中的权重值W，最后运用模糊运算矩阵计算出定量解。

5.1   评价指标

评价指标主要由孕灾背景和潜在地质灾害两个指标组成，其中孕灾背景指标又包括地震区划、微观地貌因子、宏观地貌因子和底质类型四种指标^[5]。研究区位于菲律宾海盆中南部，水深为–6712.8～–1692.2 m，区域内无大型地震，受热带气旋影响不大，基于以上特点，综合考虑各个指标特性，将5类指标划分为5个等级，建立海底稳定性评价指标体系（表2）。数据来源于海底地震烈度区划图、海底地质灾害分布图、微观与宏观地貌因子分布图（通过对坡度、坡度变率、坡向、坡向变率、曲率5种微观地貌因子和地表粗糙度、地形起伏度、高程变异系数、地表切割深度4种宏观地貌因子进行综合分析，参考已有研究成果，分别选择坡度和地形起伏度作为微观与宏观地貌因子的指标）以及底质类型分布图。

表  2  海底稳定性评价指标体系

Table  2.  Evaluation index system for seabed stability

指标体系	地震区划	灾害地质	微观地貌因子（坡度）	宏观地貌因子（地形起伏度）	底质类型
1级	地震动峰加速度值=0.05 g，相当于5级地震区	海底火山	0°～3°	0～19 m	基岩
2级	地震动峰加速度值=0.1 g，相当于6级地震区	裂谷	3°～7°	19～45 m	铁锰结壳
3级	地震动峰加速度值=0.15 g，相当于7级地震区	海底滑坡	7°～15°	45～77 m	含铁锰结核的远洋黏土
4级	地震动峰加速度值=0.2 g，相当于8级地震区	陡坎	15°～25°	77～120 m	远洋黏土
5级	地震动峰加速度值≥0.3 g，相当于9级地震区	崩塌	＞25°	＞120 m	硅藻软泥

下载: 导出CSV 

| 显示表格

根据选取的5类评价指标，可以得到稳定性评价因素集A（式1）。

a₁—地震区划指标，a₂—灾害地质指标，a₃—微观地貌因子指标，a₄—宏观地貌因子指标，a₅—底质类型指标。

5.2   海底稳定程度分级

逻辑信息分类法和特征分类法是常用的分级方法，通常将级别划分为3级和5级^[35]。考虑到研究区具体情况，本文选择五级分类体系，分为不稳定、较不稳定、中等稳定、较稳定和基本稳定5个等级区，得到模糊评语集X。

v₁—基本稳定，v₂—较稳定，v₃—中等稳定，v₄—较不稳定性，v₅—不稳定。

5.3   确定评价因子归属度

评价因子归属度反映了评价因子隶属于评语集的程度，通常由隶属函数确定。由于研究区位于深海，难以用具有连续变化特点的函数对其进行模糊信息化，因此选取具有正态分布特点的“0，0.2，0.5，0.8，1.0，0.8，0.5，0.2，0”作为隶属度函数，对海底稳定性评价因子的影响程度进行表示（表3）。通过例证法、正太型函数分布规律以及专家经验法得到隶属度值，各评价因子的隶属程度应基本相等^{[23, 36]}。

表  3  海底稳定性评价指标隶属度值

Table  3.  Membership value of evaluation index of seabed stability

稳定性 等级	评价指标
	地震区划	灾害地质	微观地貌因子	宏观地貌因子	底质类型
基本稳定	1级	1级	1级	1级	1级
v1	0.15	0.15	0.15	0.15	0.15
较稳定	2级	2级	2级	2级	2级
v2	0.30	0.30	0.30	0.30	0.30
中等稳定	3级	3级	3级	3级	3级
v3	0.55	0.55	0.55	0.55	0.55
较不稳定	4级	4级	4级	4级	4级
v4	0.80	0.80	0.80	0.80	0.80
不稳定	5级	5级	5级	5级	5级
v5	0.85	0.85	0.85	0.85	0.85

下载: 导出CSV 

| 显示表格

5.4   构建模糊评价矩阵

利用获取的评价指标隶属度值对各评价指标进行量化，得到稳定性评价因素集A到模糊评语集X的模糊评价矩阵R。

v₁，v₂，v₃，v₄，v₅—各评价因子的隶属度值；r₁—地震区划指标；r₂—灾害地质指标；r₃—微观地貌因子指标；r₄—宏观地貌因子指标；r₅—底质类型指标。

5.5   确定指标权重

评价指标的权重值反映了评价因子间的相对重要性，对海底稳定性评价具有重要意义。本文采用专家评判与层次分析决策相结合的方法确定海底稳定性评价指标的权重值^[37]。首先通过专家对各指标的相对重要性进行判断，然后利用层次分析过程建立权重值判别表，得到权重值判别矩阵，最终计算出各评价指标的权重值分配集W（表4）。

表  4  评价指标权重值

Table  4.  The weight value of the evaluation index

评价指标	地震区划	灾害地质	微观地貌因子	宏观地貌因子	底质类型
权重值	0.0882	0.4412	0.2206	0.1471	0.1029

下载: 导出CSV 

| 显示表格

5.6   计算结果

通过建立的模糊关系矩阵R和权值分配集W，以3'×3'大小的网格单元对研究区海底稳定性进行评价，利用模糊运算矩阵得到综合评价结果B。

W—权重值分配集，R—模糊评价矩阵，*—模糊变换算子。

评价结果B是基于评价因素集A的综合评价结果，根据模糊数学的最大隶属度原则，取隶属度最大者所对应的等级作为评价单元海底稳定性等级。

5.7   海底稳定性区划

以3'×3'大小的网格为评价单元，依据模糊数学原理构建的模糊综合评价系统进行计算，采用优势合并方法，将优势的等级值作为该区的稳定性等级，得到研究区海底稳定性区划图（图3）。

图  3  海底稳定性区划图

Figure  3.  Division map of submarine stability

下载: 全尺寸图片 幻灯片

（1）基本稳定区

研究区基本稳定区主要分布在中北部。该区域地形地貌相对平坦，缺乏大规模地貌单元发育，受宏观、微观地貌因子影响较小。此外，由于构造、岩浆热液等地质运动相对匮乏，灾害地质因子并不发育，仅有局部地区分布少量海底火山，因此具有较高的稳定性。

（2）较稳定区

较稳定区广泛分布在各地貌单元过渡区。该区地形地貌的坡度和地形起伏度通常不大，底质类型多为基岩和铁锰结核，一般受地质灾害影响较小，较为稳定。

（3）中等稳定区

中等稳定区在研究区各个范围内都有发育。该区通常具有一定的坡度和地形起伏度，多位于大规模地貌的过渡区或较为陡峭的微地貌单元处。受地形起伏变化和底质类型影响，局部可能发育海底滑坡等较严重的地质灾害，但总体还是比较稳定。

（4）较不稳定区

较不稳定区零散分布在研究区内。该分区主要位于地形地貌变化较大的区域，可能受海底滑坡、陡坎等地质灾害影响。此外，局部地区的不稳定也可能与硅藻软泥的分布有关。

（5）不稳定区

不稳定区大多位于海山、山间盆地等地形地貌起伏变化较大、坡度较陡的区域。该分区通常受灾害地质因子影响较大，是海底滑坡、陡坎、崩塌等地质灾害的易发区域。这些区域的底质类型多为远洋黏土和硅藻软泥，稳定性相对较低。

6.   结论和建议

（1）西太平洋菲律宾海微地貌单元发育，其海底稳定性与地形地貌变化密切相关。稳定区通常较为平坦，缺乏地貌单元和灾害地质干扰；不稳定区主要集中在九州-帕劳海岭及其西部的海山、山间盆地等大规模地貌单元发育区；较稳定、中等稳定、较不稳定区广泛分布在各地貌单元及地貌单元间的过渡地带。

（2）基于研究区自身特点，选取地震区划、灾害地质、微观地貌因子、宏观地貌因子和底质类型作为评价因子，运用模糊数学原理构建稳定性评价模型，建立了5类5级评价指标体系。评价结果能较为合理地显示研究区海底稳定性，与实际情况基本吻合。

（3）深远海海底稳定性评价对于我国未来的海洋资源开发、海洋经济发展具有重要指导意义。然而，该领域目前还存在很大空白，技术规程有待完善，研究成果也大多集中在近海。建议规范深远海海底稳定性评价体系，为进一步建设深远海海洋工程奠定基础。