8 v1 x+ \1 F8 P- Q一、引言 裂流是一种起源于破波带的快速离岸水流，是近岸环流重要的组成部分。它可以将泥沙等沉积物带离近岸，是泥沙向海输移的重要途径。裂流不仅会加剧局部海岸线和沙丘的侵蚀，对沿海基础设施造成破坏，而且会严重威胁海滩上游客的生命安全。 陆旭等建立了我国海滩事故数据库，发现除海浪和潮汐外，裂流很可能是造成溺水事故的重要因素。由于裂流没有引起海滩游客的广泛关注，甚至有游客将其当成是快速进入深水区的通道，以致海滩安全事故频繁发生，因此裂流逐渐引起许多学者和专家的重视。国际上已经有许多国家开展了对裂流的观测与研究，如美国国家气象局（NWS）已经开始预报美国所有海滩的裂流情况，韩国气象局建立了一套预报系统以实现对裂流的监测，我国原国家海洋局海洋减灾中心也在2018年联合各部门对全国重点滨海旅游海滩的裂流灾害风险进行了排查和警示工作。 + O3 U/ F" v+ \6 K! V
破波带受到波浪、潮汐、地形等影响随时都在变化，产生于破波带的裂流往往也具有一定的随机性，另外在破波带进行现场观测是困难和昂贵的，而且也具有相当的危险性。因此，使用高空视频或者遥感影像对裂流进行研究有着巨大的优势并受到一些学者的青睐。Lin等使用航空影像和架设在高处的摄像机获得了台湾福隆海滩的图像，目视解译了影像中海滩的沙坝形态和裂流位置，并研究了福隆海滩类型和裂流分布的短期和长期变动。Maryan等比较了Viola-Jones算法、卷积神经网络算法和元学习器对裂流槽图像数据集的识别效果。Silva等在GoogleEarth上分别提取了1740张含裂流和700张不含裂流的图像，使用FasterR-CNN和自定义时间聚合算法对裂流图像进行了识别检测，识别精度高于人类和其他裂流检测算法。本研究利用前人提出的裂流分类理论，基于GoogleEarth的历史影像对广东省沿岸14个主要旅游海滩的裂流类型进行了分类，并开展了海滩形态及海湾形状与裂流关系的初步分析。
; D$ O% t# e* Y3 K: m
二、裂流的分类 Shepard在1936年提出“裂流”这一概念，并将其描述为快速流向大海的水流。Long和Özkan-Haller通过总结前人的研究指出裂流是由于岸滩入射波高的变化和波浪增减水在沿岸方向分布不均匀，导致水位分布不均匀，从而产生辐射应力驱动生成的。 根据裂流发生的原因，Short和Brander提出4种裂流类型：淤积型海滩裂流、侵蚀型海滩裂流、地形裂流和巨型裂流。淤积型海滩裂流通常在某一位置固定出现数日至数周，侵蚀型海滩裂流则倾向于沿海岸迁移，两者统称为开阔海岸裂流。地形裂流是指由地形结构控制的裂流，经常发生在丁坝、海岬等地形附近。巨型裂流是一种大型地形裂流，与开阔海岸裂流相比可以延伸3～4个破波带的宽度。 Dalrymple等把裂流分为两类:水深控制的裂流和水动力控制的裂流，并将常见的海滩形态具体描述为具有下切沟槽的沿岸沙坝−沟槽型海滩和具有下切沟槽的阶地型海滩。另外，在陡峭的海滩滩角处发现了一种常见的裂流类型，被称为冲流带裂流。它是波浪上升过程中由冲流过程驱动的裂流，对海滩游客和海岸的危害有限。水动力控制的裂流仅仅是由水动力机制驱动的，有着更小的离岸强度和范围。 Castelle等认为以是否出现在开阔海岸或者海湾海滩来区分裂流是具有误导性的。实际上许多海湾海滩上发生的裂流并不都靠近岬角或海岸结构，它们表现出与开阔海岸裂流相同的特征。所有的裂流都是由水动力驱动的，但是控制机制的主导程度不同。Castelle等根据控制机制所占的主导地位，将裂流分为三大类。由于物理驱动机制不同，又进一步将每个类别分为两种类型，总共有6种基本类型。在现实场景中，各种驱动机制往往结合在一起，因此Castelle等又定义了4种混合型裂流类型（表1），冲流带裂流因为危害程度有限没有被列入分类中。 Castelle等结合了其他裂流分类的优点，侧重于从控制机制的主导程度对裂流进行区分，考虑到海湾海滩也会产生类似开阔海岸裂流的特征并针对现实情况提出了4种混合型裂流类型，为理解和预测裂流提供了进一步的框架。 表1 Castelle裂流分类类型表 本文选取了广东省14个主要旅游海滩，在GoogleEarth上对每个海滩的历史影像中的裂流进行目视解译，采用Castelle的裂流分类模型对裂流进行了分类。冲流带裂流是这14个海滩明显的特征，因此也被考虑在本文的分类之中。
. f9 s& [ Z$ D6 m1 {+ d
三、方法和数据 ⒈裂流解译标准 不同的裂流种类有着不同的驱动机制，研究海滩上裂流种类的组成对研究海滩的状态有重要的参考价值。目前，对于裂流分类的理论研究非常成熟，裂流的分类模型也在不断完善。本文以Castelle提出的裂流分类框架为基础（图1），根据其对每种裂流形态的描述作为解译的依据，在GoogleEarth上查看各个海滩的历史影像，根据分类理论选取每种裂流的特征影像（图2），提出针对影像的裂流解译标准（表2）。
9 ^# u6 M0 c! c% \0 L z1 s" p4 {3 @图1 裂流分类框架图
. K% r) q3 W: @ Pa.剪切不稳定裂流；b.瞬变裂流；c.沟槽裂流；d.聚焦型裂流；e.偏斜裂流；f.阴影区裂流；g.沟槽–瞬变裂流；h.聚焦型–沟槽裂流；i.边界–沟槽裂流；j.岬湾单元裂流；k.冲流带裂流 图2 各类型裂流的特征影像 表2 影像各类型裂流解译判释表 | | | 位于沿岸均匀的开阔海滩，波浪高度倾斜入射，形成离岸的裂流头和狭窄的裂流颈，或物质离岸输移使颜色变深 | | 位于沿岸均匀的开阔海滩，波浪近垂直入射，形成离岸的裂流头和狭窄的裂流颈，或物质离岸输移使颜色变深 | | 沟槽切割沙坝或阶地并在开阔海滩上规律性地呈现；沟槽两侧波浪破碎形成白沫，影像上表现为白色 | | 入射波峰线由海向岸逐渐弯曲成韵律状，沿岸有白沫规律性凸起，在凸起附近有物质离岸输移使颜色变深 | | | | | | 波浪近垂直入射，沟槽宽度很小，在沟槽附近有细小的离岸白沫 | | 有两层白沫(双沙坝)，外侧白沫相对不明显，内沙坝沟槽切割程度较深，沟槽宽度很大且呈现规律性分布 | | | | 海湾凹入度小于冲浪区宽度，海湾内白沫非常明显，在海岬两侧或中心有离岸白沫 | | |
⒉基于地形动力指标的裂流风险评价模型 Wright等对澳大利亚一系列的海滩进行现场调查后提出了无量纲沉降速率参数（Ω）模型。后来，Masselink和Short在Ω模型的基础上引入了相对潮差（RTR），建立了更系统和全面的分类模型。李志强等结合Ω模型和RTR模型提出了一种操作简单，能定量区分裂流出现的可能性的裂流风险性评价方法。该方法受到人们广泛的认可，现在已经应用在浴场选址、裂流风险排查等项目中。这两个参数的定义为： Q＝Hb/TWs， RTR＝TR/Hb 式中，Hb是局地平均破波波高；T是平均波周期；Ws表示泥沙沉降速率，TR是大潮平均潮差。根据参数Ω和RTR的阈值范围可将海滩划分为4组8类，并且各自对应不同的裂流风险等级（表3）。基于地形动力的海滩状态分类模型基本能够反映海滩对波浪、潮汐、泥沙输移等因素的响应，已经被地质地貌领域专家广泛接受。 表3 基于地形动力指标的滨海旅游海滩裂流风险等级表
⒊数据获取
- Q- i" z, R8 I1 C$ G 本文选取广东省沿岸14个主要旅游海滩，如图3所示，其中有9个是原国家海洋局2018年进行裂流灾害风险排查的广东省重点滨海旅游海滩。GoogleEarth的影像时间跨度为2002–2021年。使用GoogleEarth测量工具，测量每个海滩的长度，测量时以海滩向陆边界为界限来减少因潮汐导致的海滩长度的变化。采用表2中的解译标准通过目视解译法将各海滩出现过的裂流从历史影像中判别出来，再进行分类并统计，共获得7334个不同类型的裂流。需要注意的是在裂流发生时，一些狭长的裂流会参差交错，因此本文认为当多个裂流交叉最后一起向海流出时，视为一个裂流。表4是对14个海滩的GoogleEarth影像进行统计的结果。
! [, w5 G( _1 |4 k n* X* M* U注：平海东、平海西分别是指广东省惠州市平海镇东侧和西侧海湾 图3 14个海滩的位置示意图 表4 14个海滩的影像统计
" b5 J7 k7 b/ z7 w4 i注：由于云量或光线等原因会造成某些影像无法识别
9 A1 Q% }6 Y. Z) _. |$ c 某些海滩的GoogleEarth影像不能完全覆盖，往往是由多年的历史影像拼接而成，所以影像的数量不能表示完整海滩的数量。针对不完整的海滩影像，按照影像中海滩的长度占原始海滩长度的比例计算每景影像的有效影像数量（表5），平均裂流数等于当月的海滩裂流总数与有效影像数量的比值。 # m" E. }7 }8 U! |3 `/ d! Y
5 p3 \% t1 r' H" {5 L0 q* a
( o0 q% ^# n n3 Z! Y2 f
四、裂流的统计特征 ⒈逐月分类结果与裂流风险评价模型对比 本研究收集到青澳湾、广澳湾、红海湾东、西涌海滩、十里银滩、大角湾、月亮湾以及东海岛8个海滩的裂流风险等级，数据来自广东省裂流风险排查报告以及李志强和陈杏文的结果。选取广澳湾、月亮湾和东海岛为例进行讨论，它们分别是这8个海滩中有效影像数量最多、最少和位于中间值的海滩。 广澳湾的有效影像数量最多，影像多分布在1月、9月、10月、11月和12月份，涵盖了风险评价模型对该海滩分类的全部范围，4月和7月缺少影像数据。如图4所示，广澳湾在1–9月裂流风险等级为高，10–12月为中，整体上看属于裂流风险危害性较大的海滩。 6 ^0 T" Q F* N6 D' h
图4 广澳湾裂流逐月分类统计* `0 G7 } _; k6 |2 N2 o
从图中可以看出，裂流风险等级为高的月份，其裂流的组成较为单一，混合型裂流较少，水深控制的裂流组成比例最大，而且其平均数量明显大于裂流风险等级为中的月份，符合中间状态海滩的特征。查看排查报告中的Ω值和RTR值发现10–12月的海滩为沙坝消散型海滩，虽然该海滩风险等级为中，但在RTR值与其他各月在相同范围内的情况下，其Ω值的计算结果接近分界值，裂流风险也相对较高，10–12月的平均裂流数量也反映出这一点，这可能是平均裂流数没有与1月表现出明显差距的一个原因。本研究针对广澳湾的裂流分类在一定程度上符合了其裂流风险型性评价方法。 东海岛的有效影像数量相对广澳湾少了近一半，而且主要分布在风险评价等级为中的月份，在低等级月份无影像，在高等级月份的有效影像数量仅为0.5个。从图5可以发现，3月水深控制的裂流比例最大，其他月份水动力控制的裂流或混合型裂流所占比例较大。除去1月和9月，剩下中等级月份的平均裂流数量接近持平。4月、10–12月的平均裂流数量差距不大，并且水深控制的裂流比例较低，符合中风险等级的海滩评价,但由于缺乏其他风险等级的海滩影像，所以导致对东海岛的裂流分类与评价模型无法完全对应。 0 |5 n% P1 A, C! R4 j
图5 东海岛裂流逐月分类统计) p5 y7 e+ J: N! x! f" l$ e
月亮湾仅有4景影像（图6）。在10月，月亮湾冲流带裂流的平均裂流数量很多，与其他海滩的冲流带裂流不同的是，其波浪的冲流过程更强，对海滩的侵蚀程度更大。不过由于该海滩的有效影像数量少，能解译出的裂流数量非常有限，因此对比效果最差。
1 O1 R2 h, T% Z5 D图6 月亮湾裂流逐月分类统计: y* |( I; }# S. c+ O. v( S' |: F
通过上述研究发现，裂流分类结果与裂流评价模型的对比效果与有效影像数量有直接的关系，对比效果越好，海滩的有效影像数量越大。不过在图4中5月、8月和图5中1月、9月的平均裂流数量与同裂流风险等级的月份差距很大，另外也出现图4中8月份裂流全部由混合裂流组成的情况，这与该月有效影像数量较少有关。 ⒉海滩裂流的类型结构 剪切不稳定裂流和瞬变裂流在影像上主要是通过波浪的入射角度来区分。GoogleEarth的一些影像由于图像质量等原因无法分辨入射波方向。另外水动力控制的裂流在空间和时间上都是短暂的，因此不可预测。受条件的限制，解译上述裂流时不作区分，统一解译为大类裂流。根据各种裂流占总裂流的比例计算百分比，图7为绘制的14个海滩的裂流类型结构图，下面对结构图进行分析。
5 m% A/ _, o% {" `3 I# }' v u. b 图7 14个海滩的裂流类型结构 冲流带裂流被认为是危害程度最低的裂流类型。从结构图可以看出，冲流带裂流频率最高的前3个海滩分别是：小梅沙、月亮湾和大梅沙，均超过了80%，其中小梅沙海滩完全由冲流带裂流组成，大梅沙海滩有少部分由水动力控制的裂流，月亮湾有比例很低的沟槽裂流。大梅沙、小梅沙海滩的裂流主要由冲流带裂流组成，从裂流结构上很容易判断其裂流危险性较低。东海岛无冲流带裂流。 水动力控制的裂流具有不可预测性，而且持续时间相对较短，一般为3～5min，裂流的离岸强度相比于沟槽裂流较弱，但危险性依然很高。青澳湾、平海西和东海岛的水动力控制的裂流的组成比例最高，其他海滩除了小梅沙、大角湾和月亮湾均有这种裂流类型，而且比例低于沟槽裂流。 沟槽裂流与中等海滩状态相关联，是记录最多、理解最透彻的裂流类型，结构图也反映出这一点。除大梅沙、小梅沙外，其他海滩均出现了沟槽裂流，并且是主要的结构类型。沟槽−瞬变裂流、聚焦型−沟槽裂流和边界−沟槽裂流形成时也受到沟槽的影响，加上这些裂流后由沟槽影响的裂流比重将进一步扩大。 聚焦型−沟槽裂流一般出现在双沙坝系统的海滩上，波浪折射并且穿过外沙坝导致波浪破碎后，内沙坝深度的沿岸变化进一步引起破波高沿岸变化，驱动了这种裂流的产生。青澳湾和东海岛的聚焦型−沟槽裂流的组成比例较高，说明相对来说这两个海滩容易出现双沙坝系统。 边界控制的裂流如偏斜裂流和阴影区裂流的比例很小。一方面是由于边界裂流只在刚性边界产生，产生所需条件相对苛刻；另一方面是因为人工建筑少和大多数海滩的岬角弧形过大，不易产生。平海东海滩长14.4km，包括3个岬角距离很小的海湾海滩，因此在平海东海滩会产生边界控制的裂流和岬湾单元裂流，特别是岬湾单元裂流，只出现在海湾海滩中。另外，广澳湾、大角湾也有一定比例的偏斜裂流和阴影区裂流。 其他裂流类型由于比例过小或特征不够典型而不作具体分析。通过上述分析可知，裂流的类型结构与海滩的状态密切相关，受沟槽影响的裂流与水动力控制的裂流是海滩上最容易出现的裂流类型。 ⒊海湾形状与裂流 研究的14个海滩中有11个位于受岬角控制的海湾中，其中十里银滩、水东湾和东海岛是无岬角控制的弧形海滩。两端岬角之间的距离为海湾长度；海滩最凹处至海湾长度线上的垂直距离为海湾凹入度，使用GoogleEarth测量工具测量海滩长度以及海湾长度和海湾凹入度。 从图8的变化趋势看，年平均裂流数随着海滩长度而变化，但在水东湾和东海岛处有较大偏差。广澳湾于2013年开始修建防坡堤，改变了该海湾自然岬角对海滩的控制，因此对广澳湾的海湾长度修正为防坡堤两端的距离，海湾长度也相应改变。排除非岬角海湾后发现，年平均裂流数与海湾长度变化趋势的符合程度更高（图9）。
( F6 g! Y1 s$ M5 y' B+ |图8 14个海滩的海滩长度和年平均裂流数 ) J6 e( A: E0 q8 l4 F8 e9 G
图9 岬角海湾的海湾长度和年平均裂流数( {; z( I2 ]9 v9 q0 b7 \/ \
图10显示，海滩出现冲流带裂流的比例与海湾凹入度存在一定的负相关性，而沟槽裂流则刚好相反。不过也有例外，如海湾凹入度小的青澳湾的冲流带裂流比例也很低。水动力控制的裂流和凹入度在青澳湾、广澳湾、海门湾、红海湾东、平海东、平海西、西涌海滩的变化趋势相同，在大小梅沙、大角湾和月亮湾无明显关系。 P8 j7 l, |# @7 U" t" t
图10 岬角海湾的裂流类型比例和海湾凹入度
" R' G' H# g3 i3 o s+ U: X 上述统计规律显示，海湾长度越长的海滩出现的裂流数量越多，海湾凹入度可以影响裂流的组成比例，海湾凹入度越长的海湾，沟槽裂流所占的比例越高。不过考虑研究条件的限制，这种裂流与海湾长度和凹入度的关系还有待进一步的研究。 U' k% i) y B+ R
五、结论和展望 随着海滩逐渐成为许多人出游热门目的地，海滩事故的发生率也逐渐增加。作为“海滩第一杀手”的裂流成为众多学者关注和研究的对象。裂流的产生原因和分类也一直是众多学者研究的重要课题。各种裂流的离岸强度有所区别，对裂流进行分类可以帮助理解海滩裂流的类型结构，与裂流风险性评价模型互为补充。不同的裂流类型会产生不同的视觉特征，例如沟槽裂流，与两侧的浅沙坝相比，中间波浪破碎相对较少，在视觉上呈现为一个更暗、更平静的水域，缺乏经验的游客通常会误以为这片较平静的水域是最安全的游泳场所。提高公众对各类裂流的认识，有助于降低海滩事故的伤亡概率。 本文选取了广东省14个海滩作为研究对象，根据Castelle提出的裂流分类理论在GoogleEarth上选取特征影像，提出了针对影像的裂流解译标准，随后对各个海滩的历史影像进行解译分类，获得了裂流分类数据。裂流逐月分类结果符合裂流风险性评价模型，从裂流的逐月分布来看，在裂流风险等级为高的月份中，海滩裂流的组成结构较为单一，混合型裂流少，水深控制的裂流比例大于其他裂流类型；当裂流风险等级降为中时，水动力控制的裂流和混合型裂流的组成比例增加，水深控制的裂流的比例相应减少；在裂流风险等级为低的月份出现的裂流很少，甚至无裂流发生。裂流的类型结构与海滩的状态密切相关，受沟槽影响的裂流和水动力控制的裂流是海滩上主要的裂流类型。通过统计还发现海湾形状与裂流产生具有一定的规律：海湾长度和海湾凹入度分别会影响裂流出现的数量和不同类型裂流的组成比例。 在使用GoogleEarth进行研究时遇到了一些限制因素，如：⑴海滩的影像不完整，有些海滩往往是由多景影像拼接而成的；⑵影像的时间不连续，甚至某些月份的影像缺失，以致于每个海滩的历史影像数量不一致；⑶影像质量优劣不一导致识别困难。针对第一点和第三点，本文在解译和统计时提出了一些相应的解决办法：对于海滩影像的不完整，使用有效影像数计算平均裂流数量；解译剪切不稳定裂流和瞬变裂流时不作区分，归类为大类。另外，由于解译裂流受到解译者主观判断的影响较大，以后可以通过现场观测的方式训练解译者的识别能力，并提出一套基于现场观测的裂流解译标准来提高裂流解译的精确性。 对数据进行分析后发现，裂流逐月分类结果与有效影像数量有直接的关系，有效影像数最多的广澳湾的裂流逐月分类结果与裂流风险评价模型吻合相对较好。裂流在破波带产生，受破波带状态和宽度的限制，裂流的识别对影像的时间和空间分辨率要求很高。目前，常用的高分辨率卫星有美国的World-view、法国的SPOT等。我国也在2016年发射了高分三号卫星，其空间分辨率能达到1m。如果增加统计影像的数量，可以提高分类结果与裂流风险性评价模型的吻合程度。裂流的分类对海滩研究具有重要的意义，在后期的工作中进一步提高影像数量和解译精度，将获得更为精确的分类效果。 近年来，图像分类与目标检测技术成为计算机领域研究的热点。已经有学者将卷积神经网络技术应用在沟槽裂流的识别上，识别结果比人类识别结果精度更高。不同类型的裂流具有不同特征，使用图像识别技术对各类裂流进行分类也是未来的研究方向之一。
9 i# p- B1 M4 B+ U1 z6 R: y! C | 
