摘 要：泥沙运动规律的研究是钱塘江河口河床演变研究的核心问题之一。借助卫星遥感新技术，建立了适合于钱塘江河口水域的含沙量遥感定量模式，经计算、处理、转化卫星数据成悬沙分布状况图，为定量分析钱塘江河口的悬沙特征提供帮助。
关键词：含沙量；etm图像；钱塘江河口
1 问题的由来 钱塘江河口悬沙的时空分布是研究钱塘江河口冲、淤变化，泥沙运动的动力学机制和各种主控因子的主要依据，钱塘江河口悬沙的运移和沉降过程，和潮流、地形、地貌形态的变化有直接关系，研究钱塘江河口的悬沙分布变化规律，对于岸滩演变，滩涂的开发利用、河口规划、航道、港口建设等都具有重要意义。若干年来，钱塘江河口悬沙浓度的调查研究、分析模拟所采用的主要手段是现场取样分析，其调查速度慢、周期长，且只能获得在时间、空间分布上都很离散的少量数据，而河口海岸地区水流情况复杂多变，含沙量的时间变化率很大，这种在时、空分布上很离散的采样数据对比精度很差。而用水力学、河流动力学等理论计算，也难以得到满意结果。卫星遥感（rs）技术的发展使这一状况得到了改观，采用遥感定量技术瞬间同步性极好，重复获取数据的周期短，能迅速地获得大面积水域含沙量的资料，能有效地监测含沙量的分布和动态变化，克服了现有常规方法的不足，结合现有的资料分析、数学模型、物理模型等研究方法，相互补充、相互印证，能更好地满足理论研究和海域规划、开发的需要，无疑地，rs加快了悬沙分布资料的积累，势必大大加快悬沙运移规律的研究和认识。根据etm图像和准同步的实测资料探讨了钱塘江河口遥感悬沙定量模式，将卫星数据转化成为了大面积、同步性较好的表层悬沙浓度分布图，进而分析了钱塘江河口泥沙的运移扩散规律。
2 资料获取与处理
2．1 资料选取
根据钱塘江河口的特点，其变化主要为冬夏之差，年际变化相对较小，并且洪潮比小，因此在选择卫片时，主要考虑不同潮情（大潮、中潮、小潮、涨潮、落潮）、风情（强风、常风）和季节的组合状况，并且还要有对应的实测资料，钱塘江河口局部实测水沙资料较多，大范围资料较少，能与遥感影像相对应的大范围实测资料更少。通过调查分析，本文选取landsat—7卫星，2000．09．02的etm图像进行试验，对应的潮汐情况见表1，测点分布见图1。
2．2 工作原理和图象处理
悬沙水体的光谱反射特征是建立遥感信息—悬沙浓度定量模型的重要基础，已有的实验表明含沙水体的反射率光谱有两个峰值：位于黄光波段（560～590nm）的主峰和位于近红外波段（760～1100nm）的次级峰。许多研究者研究了两个峰值所在波段海面光谱反射率（rw）和表层含沙量（s）的关系，发现rw－s关系具有下列特点：
（1）rw随着s的增加而增加，即drw／ds＞0；
（2）变化率drw／ds不是常量，而是随着s的增加而减少，即d2rw／ds2＜0；
（3）s＝0时，rw为一大于0的常量；s较大时，rw随s的增加迅速趋于一小于1的极值，即lims→0，rw（s）＞0和1ims→＋∞，rw（s）＜1；
（4）s＝0时，黄光波段光谱反射率rw1和近红外波段光谱反射串rw2与清水相同，趋于一个接近0的常量，此时rw1≈rw2；
（6）s＝st时，rw1＝rw2，rw1－s关系曲线和rw2－s关系曲线有1个对应于st的交点。
2．3 遥感模式的分类
上述特点为建立含沙量与反射比的关系提供了指导，为悬沙扩散的定量分析奠定了理论基础。根据这些特点，众多的学者提出了不同的模式来模拟悬浮泥沙的含量与遥感数据的关系，这些模式的特点一般都是半理论半经验，它们在形式上可分成4类：
（1）线性关系式：r＝a＋bs
该式是最早被学者提出来作为模拟悬浮泥沙含量与遥感数据关系的最简单的关系式。r为水面光谱反射率，a和b为与水体光学性质和太阳与传感器几何关系有关的经验参数，线性关系式仅满足实测关系中的第1和第4个特征，属有限的线性区间内近似表达式。
（2）对数关系式：r＝a＋blogs
式中r为水面光谱反射率，a和b为与水体光学性质和太阳与传感器几何关系有关的经验参数，对数关系式仅满足实测关系中的第1，2和第4个特征，但反映高含沙量特性的第3和第5，6个特征差别较大，因此在悬浮泥沙浓度不高的情况下，该式能较真实地反映悬浮泥沙与卫星数据的相关关系，但当悬浮泥沙浓度很高时，该式有很大的误差。
（3）gordon关系式：r＝c＋s／（a＋bs）
式中a、b、c为系数，它是根据准单散射近似公式得到，在实际应用中用得不多，但较有名，cordon关系式满足实测关系中的第1～6个特征，适用区间包括低含沙量和高含沙量区。
（4）负指数关系式：r＝a＋b（1－e－ds）或in（d－l）＝a＋bs
