张 杰 ,张坤军 ,陈佳兵,李京兵
(1.浙江省水利河口研究院,浙江 杭州 310020;2.浙江省水利防灾减灾重点实验室,浙江 杭州 310020)
底泥是河湖库塘等水域底部淤积土的总称,主要为淤泥土。淤泥土类是指在静水或缓慢的流水环境中沉积,或伴有生物化学作用形成的粘性土,其天然含水量大于液限,天然孔隙比大于或等于1.0[1]。通常,当河湖库塘淤积存量较大时,它将直接影响水环境质量、行洪能力和通航安全,因此对底泥进行探测调查,对于海岸工程、水利工程等基础建设具有重要意义。目前,底泥探测方法主要有静力触探法、钻孔取样法、声学探测法和放射线测量法[2]。静力触探法、钻孔取样法受限于单点测定,无法进行连续测量,在效率和作业成本上具有很大劣势。声学探测法具有低成本、高效率、高精度的优点,是目前水下勘探和调查中应用最为广泛的探测方法[3]。声波相较于电磁波在水中的衰减程度最小[4],将其应用于近海、江河、湖泊的底泥探测有助于提高作业效率与精度,为海洋工程建设、港口航道清淤、水利工程建设等提供基础技术支撑[5]。
底泥探测一直是水下底质探测领域的一个难题,国外相关文献研究资料较少,国内研究人员主要利用声学设备进行淤泥厚度的探测。谢津平[2]利用Silas系统实现对太湖湖底淤泥密度的划分,进而测量了太湖湖底淤泥量;张惟河等[6]将侧扫声呐和浅地层剖面仪相结合,探索在淤泥厚度和密度较大的科伦坡进行淤泥测量的方法,结果表明,该方法对于表层厚度较小的淤泥分辨能力较弱;张俊等[7]利用双频测深仪与泥沙溶重测试仪配合使用,对长江口深水航道地区的底部淤泥质土分布情况进行了实验室测试和实地验证测试研究,结果表明低频(33 kHz)声波相较高频(210 kHz)声波具有更强的穿透力,但在遇到14 kN/m3的重土时也无法穿透;郭发滨等[8]通过对双频测深仪的声波透射信号的分析和地质取样资料的对比,验证了利用双频测深仪进行海底浅层物质类型判别的可行性;王宝成等[9],柯勰等[10]在水库利用双频测深仪实现了水底淤泥厚度的测量,并验证了测量结果的准确性;杨姗[11]采用SES-96参量阵测深-浅地层剖面仪实现了港池的浮泥清淤测量。目前,国内研究人员对声学设备的不同声波频率与其底泥穿透能力间的关系研究较少,本文将通过实例数据进行对比分析,研究两者间的相互关系,为典型水域(滨海、河道、水库)的底泥探测提供技术参考。
声波由发射换能器发射穿过水体触及河床时,因为声阻率在水体中和水底界面的不同,造成了声波的反射,由接收换能器接受反射波,通过测定声波在水中的传播时间[12],按式(1)计算反射界面相对于换能器的距离(H):
式中:t1,2,…,n表示各反射界面对应的声波传播时间;ν表示声波传播速度。
高频声波受散射衰减影响大,穿透力弱,遇到河床容易产生反射,低频声波在软质水底中传递则相对复杂,其受散射衰减影响较小,且其能量远大于高频声波,因而具有更强的穿透性,能穿透松软底泥并产生反射,所以低频声波在垂直水底方向可出现多个反射波。测定不同反射波在地层内的传播时间,可确定每一层的厚度[13],如图1所示。到达分界面的声波能量与到达该分界面之前的声波的衰减有关,反射能量大小由反射系数决定,反射系数R为:
式(2)中:ρ1ν1,ρ2ν2分别表示上、下两层介质密度和声速的乘积(即声阻率),当相邻两层介质存在一定声阻率量差,就容易形成反射界面[13]。在一定范围内,底泥的天然含水率ω与密度ρ和声速ν存在关联性,通常当含水率越大时,密度和声速越小,当含水率越小时,密度和声速越大[14]。本文采用含水率作为综合反映各底泥分层密度和声速变化的指标。
图1 声学探测设备工作原理示意图
本文分别在滨海(浙江省舟山市金塘岛灰鳖洋附近)、河道(浙江省杭州市区中河附近)、水库(浙江省新昌县长诏水库),利用不同频率声学设备进行定点底泥穿透试验数据的采集。
底泥主要由淤泥土类组成,淤泥土类根据孔隙比或含水量分为淤泥质土、淤泥、流泥、浮泥。淤泥土类的分类见表1[1]。本文选择12 kHz,24 kHz,33 kHz和200 kHz共4组声波频率进行底泥穿透试验,各声波频率对的应声学设备见表2。以200 kHz频率作为测深参考基准,计算各频率测深数据与200 kHz频率测深数据的差值,即为各频率声波穿透深度,并采用测深杆校核仪器200 kHz测深数据,校核结果显示,测深杆读数与仪器200 kHz高频测深读数误差不超过3 cm。每次试验均选择河床平坦、底泥状态稳定的水域,根据测点水深情况统一确定测深仪换能器安装方式,将双频测深仪(频率12/200 kHz,24/200 kHz,33/200 kHz)换能器垂直安装,对同一点进行深度数据采集,之后对同一地点利用定深采样钻采集底泥柱状样,并在实验室按5~10 cm间隔分层对采集的底泥柱状样进行含水率检测。
表1 淤泥土类的分类
表2 试验参数及设备
表3为舟山市金塘岛灰鳖洋附近底泥柱状样室内试验结果。根据所取的底泥样品目力鉴别,底泥整体呈灰黄色,表部10 cm左右含水率极高,明显沁水呈油脂光泽,12 cm以下局部含条带状黑色有机物,厚度0.2~0.5 cm左右,其中50~55 cm处有机质含量高,底泥颗粒细腻、手感粘糊,随着深度增加,硬度明显有所增加。根据室内试验结果,0.5 m以上主要以流泥为主,0.5 m以下主要以淤泥为主,且均为超软土(ω>70%以上)。
表3 滨海淤泥室内试验结果
表4为杭州市区中河附近底泥柱状样室内试验结果。根据目力鉴别,底泥整体为灰色夹黑色有机物,沉积粒径较细,有气孔,比较松软。土层0~0.25 m为污染底泥层,以灰黑色为主,土质非常松软,土性以浮泥为主,含较多腐殖质,有腐臭味;0.25~0.60 m为污染过渡层,灰色夹条带状黑色有机物,具明显层理,有气孔,比较松软,土性以流泥为主;0.60~0.95 m为无污染层,主要以灰黄色为主,具层理,局部有孔洞,土性以流淤泥为主。总体上土层为超软土。
表4 河道泥室内试验结果
表5为新昌县长诏水库底泥柱状样室内试验结果。根据目力鉴别,底泥整体为灰色夹黑色有机物,沉积粒径较细,有气孔,比较松软。土层0~0.30 m为污染底泥层,以灰黑色为主,土质非常松软,土性以浮泥为主,含较多腐殖质,有腐臭味,表层含絮状悬浮物质;0.30~1.00 m为污染过渡层,灰色夹条带状黑色有机物,具明显层理,有气孔,比较松软,土性以流泥为主。
表5 水库底泥室内试验结果
选取 12 kHz,24 kHz,33 kHz 3 组频率声波的部分测深数据样本,以200 kHz测深数据为基准,分别计算各低频测深数据与它的差值,进行样本统计,即得到12 kHz,24 kHz,33 kHz声波相对于200 kHz声波的穿透深度(见表6)。
表6 各试验区域底泥穿透深度
在滨海水域,33 kHz声波穿透底泥深度约0.09~0.17 m,平均值为0.14 m;24 kHz声波穿透底泥深度约0.14~0.29 m,平均值为0.21 m;12 kHz声波穿透底泥深度约0.59~0.69 m,平均值为0.64 m。
在河道区域,33 kHz声波穿透底泥深度约0.04~0.06 m,平均值为0.05 m;24 kHz声波穿透底泥深度约0.11~0.17 m,平均值为0.14 m;12 kHz声波穿透底泥深度约0.56~0.57 m,平均值为0.57 m。
在水库区域,33 kHz声波穿透底泥深度约0.13~0.24 m,平均值为0.18 m;24 kHz声波穿透底泥深度约0.32~0.48 m,平均值为0.42 m;12 kHz声波穿透底泥深度约0.55~0.96 m,平均值为0.70 m。
各试验区底泥穿透深度与柱状样检测结果对比见图2~图4。图中3处垂向黑色虚线分别对应含水率数值55、85和150,分别为文献[1]中淤泥质土、淤泥、流泥、浮泥之间的分类界限,红色曲线为底泥柱状样分层检测得出的底泥含水率与底泥穿透深度关系曲线。
图2 滨海底泥穿透深度与柱状样检测结果对比
图3 河道底泥穿透深度与柱状样检测结果比对
图4 水库底泥穿透深度与柱状样检测结果比对
分析图2~图4可知,滨海区域底泥含水率位于70%~120%之间,总体变化幅度较小,但分层之间含水率波动较大。河道区域底泥含水率位于70%~270%之间,总体变化幅度较大,但分层之间含水率波动相对较小。水库区域底泥含水率位于90%~360%之间,总体变化幅度较大,但曲线较平滑,表明分层之间含水率几乎无波动。而含水率又与界面声阻率相关,进而直接影响不同分层间的反射系数,对声学设备在不同水域底泥探测的结果产生影响。
通过实验数据可以发现,在不同区域,同一频率声波的反射区间所对应的底泥含水率差异较大,在滨海区域对应的底泥含水率最低,水库区域次之,河道区域最高。在同一区域,随着声波频率的降低,声波反射区间对应的底泥含水率也相应减小,说明声波频率越低,其穿透能力越强。
本文利用不同频率的声波在3种典型水域进行底泥穿透试验,以实例数据对比分析了在不同区域不同频率声波的穿透能力,可以看到,随着声波频率的降低,其底泥穿透能力随之增强,符合声波物理原理。12~33 kHz频段声波可穿透河床到达浮泥、流泥层,在滨海水域中可穿透部分淤泥层。
滨海水域垂直方向上底泥含水率在分层之间波动最大,水库中底泥含水率在分层之间波动最小,因此声学设备更容易在滨海水域中获得清晰的反射层理,河道次之,在水库中则不明显。对于水下底泥探测作业,12 kHz频率的双频换能器相比24 kHz和33 kHz具有更强的穿透能力,可满足浮泥、流泥及部分淤泥层的探测。