考虑采砂影响的鄱阳湖丰水期悬浮泥沙浓度模拟*

李海军，陈晓玲，2，3，陆建忠，张 鹏，齐亨达，陈莉琼

（1：武汉大学测绘遥感信息工程国家重点实验室，武汉430079）

（2：地球空间信息技术协同创新中心，武汉430079）

（3：江西师范大学鄱阳湖湿地与流域研究教育部重点实验室，南昌330022）

（4：水利部中国科学院水工程生态研究所，水利部水工程生态效应与生态修复重点实验室，武汉430079）



考虑采砂影响的鄱阳湖丰水期悬浮泥沙浓度模拟*

李海军1，陈晓玲1，2，3，陆建忠1**，张 鹏4，齐亨达1，陈莉琼1

（1：武汉大学测绘遥感信息工程国家重点实验室，武汉430079）

（2：地球空间信息技术协同创新中心，武汉430079）

（3：江西师范大学鄱阳湖湿地与流域研究教育部重点实验室，南昌330022）

（4：水利部中国科学院水工程生态研究所，水利部水工程生态效应与生态修复重点实验室，武汉430079）

针对受采砂活动影响显著的鄱阳湖高浑浊水体，结合数值模拟和遥感技术，利用已有的鄱阳湖采砂区遥感监测结果，在构建的鄱阳湖水动力-悬浮泥沙输移模型中添加泥沙点源，对2011年7月1-31日采砂影响下的鄱阳湖丰水期悬浮泥沙浓度进行数值模拟.利用悬浮泥沙浓度实测数据和MODIS影像反演结果对模拟结果的有效验证表明，考虑采砂影响后，悬浮泥沙浓度模拟值与实测值具有强相关关系，确定性系数为0.831，均方根误差为15.5 mg/L，悬浮泥沙浓度空间分布趋势与遥感反演结果基本一致.模拟结果显示，采砂活动对鄱阳湖南部主湖区、河流入湖口影响较小，其主要影响由南向北，经棠荫以西和松门山岛以北航道、入江水道延伸到湖口区域，是鄱阳湖北湖区高浑浊水体形成的重要原因.

数值模拟；悬浮泥沙；遥感；Delft3D；鄱阳湖

©2016 by Journal of Lake Sciences

悬浮泥沙携带大量的营养盐和重金属等污染物，其能改变水体透光性，是水环境变化的重要影响因子，对水环境生态系统有着重大影响［1］.采砂活动极大地影响着悬浮泥沙时空变化，从而影响水环境.现有的水体悬浮泥沙研究大多依靠数值模拟［2-5］或遥感技术［6-7］，两者结合的研究则主要以海岸带、河口和湖泊的水动力、水质营养盐和泥沙模拟为主［8-11］.对于内陆湖泊，尤其是针对受采砂活动影响显著的鄱阳湖高浑浊水体来说，结合遥感技术，对采砂活动影响下的悬浮泥沙数值模拟鲜有研究.

自2001年长江河道采砂被禁止后，鄱阳湖采砂业迅速兴起.此后10余年间，除2008年曾禁采外，鄱阳湖采砂活动不断，湖泊悬浮泥沙浓度显著上升，引发一系列环境问题［12］.由于鄱阳湖复杂的水环境特性，水沙时空格局变化剧烈，水环境的精确数值模拟研究一直是难题.部分学者建立鄱阳湖流域的水文水动力模型，对湖泊水动力、水质、氮磷营养盐以及鄱阳湖水利枢纽工程的影响进行模拟研究［13-17］.针对采砂问题，现有的数值模拟主要集中在河流航道疏浚、港口工程以及海岸带人工造陆与倾沙等方面［18-21］.对于鄱阳湖采砂问题，一些学者利用遥感技术，建立鄱阳湖水体透明度或悬浮泥沙浓度与采砂船只数目之间的关系，分析采砂影响范围与泥沙收支平衡［22-24］.也有学者利用MODIS影像反演悬浮泥沙，对数值模型进行初始化和参数率定，间接考虑采砂影响［25］.通过建立相关关系、间接模拟的方法只能粗略地估算采砂影响，难以针对鄱阳湖采砂活动对悬浮泥沙的影响进行连续、直观和准确的分析.

本文针对鄱阳湖高浑浊水体受采砂活动影响显著的特点，借助数值模型对水动力过程、悬浮泥沙输移进行连续模拟的优势，利用采砂区遥感影像监测结果，通过在模型中添加泥沙点源的方法考虑采砂活动影响，对鄱阳湖丰水期悬浮泥沙进行数值模拟，以期更精确地模拟悬浮泥沙，为采砂活动影响下水沙时空变化及其生态效应等进一步研究奠定基础.

1　研究区域与数据

鄱阳湖（28°22′～29°45′N，115°47′～116°45′E）地处江西省北部，长江中下游南岸，是我国最大的淡水湖泊.以松门山为界，鄱阳湖被分为南北2部分，北面为入江水道，长约40 km，最窄处约2.8 km；南面为主湖区，长约133 km，最宽处达74 km.鄱阳湖主要承接修水、赣江、抚河、饶河（分支为乐安河和昌江）、信江“五河”以及博阳河、漳河来水，经调蓄后由湖口汇入长江，是一个过水性、吞吐型、季节性的大型通江湖泊（图1a）.受入湖来水和长江的双重影响，鄱阳湖水情复杂，水沙格局时空变化剧烈.

本文选取2011年7月1-31日时间段对鄱阳湖丰水期水动力及悬浮泥沙进行数值模拟.研究涉及的数据主要包括：（1）鄱阳湖及其流域水文数据：修水的虬津站和万家埠站、赣江的外洲站、抚河的李家渡站、信江的梅港站、乐安河的虎山站和昌江的渡峰坑站以及湖口站的日均流量和泥沙含量数据；星子、都昌、棠荫、湖口4个水文站的日均水位数据；鄱阳湖湖底高程以及鄱阳湖附近波阳气象站的风速、风向等气象数据. （2）现场观测数据：2011年7月15-23日鄱阳湖野外观测的50个站点（图1b）悬浮泥沙浓度数据.（3）遥感影像数据：从美国国家航空航天局NASA网站（httP：∥ladsweb.nascom.nasa.gov/data/search.html）下载的2011年7月1-31日鄱阳湖区无云或少云MODIS影像250 m分辨率地表反射率产品MOD09GQ；从美国地质调查局USGS网站（httP：∥earthexPlorer.usgs.gov/）下载的2011年7月4日鄱阳湖区Landsat ETM+清晰影像，利用ENVI去除影像条带.

2　研究方法

2.1MODLS影像水体提取与悬浮泥沙反演

对MODIS影像进行云掩膜后［26］，采用归一化植被指数NDVI（normalized difference vegetation index）阈值法进行水体提取［27］，根据水体光谱特性，选取NDVI＜0作为水体的判断条件.根据水体像元个数和分辨率计算水体面积，将水体范围矢量化后得到鄱阳湖水体边界.考虑到反演模型的区域性和精度，采用文献［28］的指数模型反演鄱阳湖悬浮泥沙浓度.

图1　研究区域：（a）鄱阳湖地理位置和相关水文站、气象站；（b）2011年7月15-23日鄱阳湖采样点（HJ1B-CCD1影像，RBG：432，2011年7月21日）Fig.1 Study area：（a）location of Lake Poyang and main hydrologic stations，meteorological station；（b）In-situ stations in Lake Poyang from July 15，2011 to July 23，2011（HJ1B-CCD1 image，RGB：432，July 21，2011）

2.2鄱阳湖悬浮泥沙输移数值模型构建

本文选取荷兰Delft3D水动力模型中的Delft3D-Flow对2011年7月1-31日鄱阳湖丰水期水动力和悬浮泥沙进行数值模拟.Delft3D-Flow采用正交曲线网格，其数值格式基于有限差分法采用ADI时间积分［29］.其耦合了泥沙模块，可以有效模拟风场驱动下的非恒定流，在模拟水位、面积和流速等水动力条件的同时，对泥沙输移进行同步模拟.

2.2.1计算网格 由于2011年7月1日没有鄱阳湖区无云MODIS影像，分析水文数据得知，都昌站2011年6月20-30日实测水位均高于模拟时间段内水位，因此采用2011年6月27日MODIS无云影像，提取鄱阳湖区水体范围，剔除湖区周边未与主湖连通的水域，以鄱阳湖水体边界为模拟范围.采用矩形网格对模拟范围进行网格划分.考虑到计算机性能和数值模型计算收敛性的要求，网格尺寸大小约为300～500 m，将湖底高程插值到网格，生成模型所需地形数据（图2）.

2.2.2边界条件和初始条件 模型边界条件为修水、赣江、抚河、信江、饶河入湖流量和悬浮泥沙浓度.修水、抚河、信江各设置1个入湖口，饶河分为昌江和乐安河2支各设1个入湖口，入湖流量和悬浮泥沙浓度采用各站点日均数据；赣江设置4个入湖口，其中西支在吴城附近入湖，另外3支在鄱阳湖西南部入湖，赣江总支入湖流量和悬浮泥沙浓度采用外洲站日均数据，根据外洲站流量设定动态分流比［30］，计算4个分支入湖口的流量和悬浮泥沙浓度.以鄱阳湖流入长江湖口为开边界，开边界水位为湖口站日均水位，开边界悬浮泥沙浓度采用Neumann边界条件［31］.模型考虑风场影响，采用波阳气象站实测风速及风向数据，作为整个湖区统一的风应力条件.模型采用1个月的热启动，不考虑水文站至河流入湖口的水量变化，以各水文站实测数据作为边界条件；以主湖区内都昌站2011年6月1日实测水位作为初始水位，初始悬浮泥沙浓度和初始流速均设为0，对2011年7月1-31日的模拟结果进行分析.

2.2.3其他参数 模拟计算的时间步长为5 min，其他参数如底部糙率、临界水深、粘滞系数、泥沙沉降速度、泥沙沉降临界剪切力、泥沙侵蚀沉降剪切力、侵蚀常量等参照文献［25］给定值进行设置.

图2　鄱阳湖计算网格（a）和湖底高程示意（b）Fig.2 ComPutational grids（a）and bottom elevation（b）of Lake Poyang

2.3考虑采砂活动的鄱阳湖悬浮泥沙浓度模拟

（1）不考虑采砂的悬浮泥沙模拟：模型配置和上文描述的一致.（2）考虑采砂的悬浮泥沙模拟：2011年7 月15-23日鄱阳湖现场观测时发现多个采砂作业区.根据Li等［32］的算法，采用2011年7月4日鄱阳湖Landsat ETM+影像对采砂船进行监测，确定3个主要采砂区，分别位于松门山岛以北（图3A）、松门山岛以东（图3B）和棠荫以西（图3C）.根据PennekamP等［33］的研究，采砂效率5500 m3/h的采砂船，采砂中心悬浮泥沙浓度约为150 mg/L.根据现场观测资料，鄱阳湖采砂船的采砂效率一般为10000 m3/h，其外围表层水体的悬浮泥沙浓度达到277.4 mg/L，因此取采砂中心悬浮泥沙浓度约为300 mg/L.遥感影像解译的结果中，剔除航道内运砂船后，A、B、C采砂区内采砂船数目分别为5、4、8.考虑到计算网格尺寸较大、同一个网格内存在多艘采砂船的情况，对模型适当简化.分别以1个泥沙点源代替各采砂区内多艘采砂船的共同影响，则A、B、C采砂区内的泥沙点源浓度分别为1500、1200和2400 mg/L，符合2011年丰水期A、B、C主要采砂区的采砂强度差异规律［32］.由于采砂船一般长时间定点作业，本文假设在整个模拟时间段内，采砂船的位置不发生变化，按照每天8：00至18：00工作10 h，模拟采砂活动对悬浮泥沙浓度的影响.

选取2011年7月4日和7月20日2景MODIS无云或少云影像反演悬浮泥沙浓度，与最接近MODIS过境时间的悬浮泥沙模拟结果进行对比.根据鄱阳湖由南向北的流场规律和地理位置特征，选取5个兴趣点（图3），分析其在不同情景下悬浮泥沙浓度变化和采砂活动影响大小.其中，a点位于棠荫以西采砂区下游，是鄱阳湖采砂活动的最上游；b点位于松门山岛以北的狭窄航道内，是鄱阳湖南北湖分界处；c点位于鄱阳湖主湖区与入江水道的连通处，鄱阳湖经此处汇入狭长的入江水道；d点位于星子附近的入江水道内；e点位于湖口，鄱阳湖经该处流入长江.

图3　鄱阳湖3个主要采砂区和5个兴趣点（Landsat ETM+第7波段，2011年7月4日）Fig.3 Three main dredging areas and five interest Points in Lake Poyang（Landsat ETM+，band 7，July 4，2011）

3　结果与讨论

3.1鄱阳湖水动力模拟验证

将星子、都昌、棠荫3个水文站的水位模拟值与实测值进行对比分析.模拟时间段内3个站点水位模拟值与实测值的均方根误差（RMSE）较小，RMSE均小于0.19 m，确定性系数R2在0.98以上（图4）.不同位置的3个站点水位模拟值与实测值比较吻合，均与实测变化趋势一致.同时间段3个站点的水位差异较小，表明鄱阳湖丰水期南北湖区水位梯度小，模型水位模拟精度较高.

由于水体面积精确值难以获取，本文利用模拟时间段内鄱阳湖区5景MODIS无云或少云影像，提取模拟范围内水体面积，对模拟结果进行验证.鄱阳湖水体面积模拟值和MODIS影像提取值比较吻合，多数日期的模拟面积与影像提取面积相对误差在4%以内，水体面积的变化趋势与水位验证趋势一致.验证结果表明本文选取Delft3D模型能够对鄱阳湖丰水期水环境进行较为准确的模拟（表1）.

图4　星子、都昌、棠荫水文站模拟水位与实测水位对比Fig.4 ComParison of simulated and measured water level at Xingzi，Duchang and Tangyin Stations

表1　水体面积模拟与MODIS影像提取值对比Tab.1 Simulated water area and MODIS images extracted area

3.2采砂影响下的鄱阳湖悬浮泥沙浓度模拟结果分析

分析不考虑采砂的鄱阳湖悬浮泥沙模拟结果发现，模拟值与实测值之间相关性较差.在18#、20#、24#、25#采样点处，悬浮泥沙浓度实测值均超过100 mg/L，模拟值远小于实测值（图5）.由此可见，若不考虑采砂的影响，鄱阳湖悬浮泥沙的浓度模拟结果很差，精度不高.考虑采砂影响后，总的来说，模拟值与实测值的变化趋势比较一致，二者数值上比较接近，R2达到0.831，均方根误差为15.5 mg/L，达到较高的模拟精度（图6）.

不同位置采样点（图1）的悬浮泥沙浓度呈现不同的变化规律.1#～17#采样点位于Ⅰ采样区入江水道，悬浮泥沙随湖流汇集，受采砂活动的综合影响，泥沙浓度上升；模拟值与实测值比较吻合.18#～25#采样点位于Ⅱ采样区，靠近3个主要采砂区，受采砂影响较大；18#、20#采样点位于入江水道与主湖区连通处，受较近的A和B采砂区的双重影响，悬浮泥沙浓度大幅增加；24#、25#采样点位于最大采砂强度的C采砂区下游，受到的影响最大.同样处于Ⅱ采样区的21#、22#、23#采样点均靠近B采砂区，却呈现出不同的变化规律：整个模拟时间段内，模型中都有采砂作业，而样点采集时间不同，实际情况中采砂作业可能中途停止，导致21#采样点模拟值高于实测值；另外，悬浮泥沙输移扩散主要受湖流影响，从B采砂区附近流场来说（图7），23#采样点位于B采砂区上游，且B采砂区采砂强度相对较小，该点受采砂影响较小；22#采样点受B、C采砂区共同影响，悬浮泥沙浓度上升.对其他采样点而言，Ⅲ采样区位于松门山岛以西的航道，由于松门山岛的阻隔，采砂活动影响较小，悬浮泥沙浓度模拟值几乎没有变化；Ⅳ、V采样区分别位于松门山岛以南主湖区和鄱阳湖东部主湖区，受下游采砂活动的影响小，悬浮泥沙浓度在2种情景下较为稳定一致.

图6　考虑采砂的悬浮泥沙浓度模拟值与实测值的线性关系Fig.6 Regression between in-situ susPended sediment concentration and simulated susPended sediment concentration with dredging

图7　B采砂区附近流场示意（2011年7月18日）Fig.7 Velocity field in dredging area B（July 18，2011）

3.3鄱阳湖采砂活动对悬浮泥沙分布的影响范围分析

分别选取2011年7月4日和2011年7月20日2组不考虑采砂和考虑采砂的模拟结果与遥感反演结果进行对比分析.2组鄱阳湖悬浮泥沙浓度空间分布可以看出（图8），不考虑采砂和考虑采砂的模拟结果存在较大差异.不考虑采砂的结果中，悬浮泥沙浓度总体较低，最高值约为80 mg/L，浓度相对较高的区域主要集中在松门山岛北部以及北部入江水道.考虑采砂影响后，悬浮泥沙浓度的空间分布更加合理，其空间相对变化趋势与MODIS遥感反演结果趋于一致.在松门山岛以北、以东和棠荫以西的3个采砂区范围内，悬浮泥沙浓度明显上升，入江水道的悬浮泥沙浓度整体上升，总体呈现由南向北浓度由高到低的空间分布趋势，与遥感反演的结果比较一致.2组模拟结果相隔时间较长，存在比较一致的空间分布规律，考虑采砂模拟结果均优于不考虑采砂模拟结果.2组模拟结果出现少部分区域异常高值的情况，可能的原因是模型在3个主要采砂区添加泥沙点源，按照采砂效率以及采砂强度设置泥沙点源浓度，导致采砂中心的网格点悬浮泥沙浓度值较高.泥沙点源浓度的估算和设置误差可能导致异常高值的出现.

图8　2011年7月4日和7月20日鄱阳湖悬浮泥沙浓度空间分布Fig.8 SPatial distribution of susPended sediment concentration on July 4 and July 20，2011

5个兴趣点悬浮泥沙浓度变化对比可以看出（图9），整个模拟时间段内，由于“五河”水沙的汇入和采砂区影响，同一个兴趣点悬浮泥沙浓度在不同情景下虽然有轻微波动，但是总体上呈现基本一致的增加趋势.不考虑采砂影响时，5个兴趣点悬浮泥沙浓度总体较低，其中a兴趣点前20 d悬浮泥沙浓度极低，至模拟后期也仅为40 mg/L左右，这与现场实测以及遥感反演的结果很不相符.由此可见，若不考虑采砂影响，鄱阳湖悬浮泥沙浓度模拟精度将受到极大影响.考虑采砂影响后，5个兴趣点悬浮泥沙浓度增加，反映出由南向北采砂活动影响由强到弱的空间变化趋势：a兴趣点受采砂强度最大的C采砂区影响，悬浮泥沙浓度变化最大；b兴趣点主要受采砂强度相对较小的B采砂区影响，悬浮泥沙浓度小幅上升；c兴趣点和d兴趣点受A、B采砂区的双重影响，相较于其上游的b兴趣点，悬浮泥沙浓度变化更大；位于湖口的e兴趣点距离3个主要采砂区虽比较远，悬浮泥沙浓度亦有所上升.

总的来说，3个主要采砂区对南部主湖区、河流入湖口的悬浮泥沙浓度影响很小，其影响主要集中在采砂区下游，在湖流作用下，由南向北经棠荫以西和松门山岛以北航道、入江水道一直延伸到湖口区域，是鄱阳湖北湖区水体相较于南湖区一直处于高浑浊状态的重要原因.

图9　5个兴趣点考虑和不考虑采砂模拟的悬浮泥沙浓度变化对比Fig.9 ComParison of simulated susPended sediment concentration for five Points of interests with and without dredging

4　结语

本文借助遥感技术和数值模拟，针对采砂活动影响下的鄱阳湖高浑浊水体，提出了一种利用已有的鄱阳湖采砂区遥感监测结果、综合考虑采砂功率和采砂强度在模型中添加泥沙点源，对采砂影响下的悬浮泥沙进行浓度模拟的方法，并对鄱阳湖3个主要采砂区的影响范围作了初步分析.利用悬浮泥沙现场实测数据和遥感反演结果对模拟结果进行的有效验证表明，考虑采砂影响后，悬浮泥沙浓度模拟值与实测值之间R2为0.831，RMSE为15.5 mg/L，达到较高的精度，其空间分布趋势与遥感反演结果趋于一致.3个主要采砂区对鄱阳湖南部主湖区、河流入湖口的影响较小，其主要影响由南向北，经棠荫以西和松门山岛以北航道、入江水道一直延伸到湖口区域，是鄱阳湖北湖区高浑浊水体形成的重要原因.研究表明，利用此方法能够有效考虑采砂活动影响，提高鄱阳湖悬浮泥沙模拟精度并对采砂影响范围进行有效分析.本文综合考虑模型要求和计算性能，采取较大尺寸的计算网格，不利于鄱阳湖复杂岸线边界内悬浮泥沙的高精度模拟，网格的改进将进一步提高悬浮泥沙模拟精度.基于现场观测资料和遥感影像确定采砂区设置了情景模拟，对数值模型模拟人为活动影响作了初步探索.利用更精确的现场观测资料和遥感影像确定采砂区域、时间、强度等，对鄱阳湖采砂影响下的悬浮泥沙进行更高精度数值模拟，分析悬浮泥沙时空变化机制，是未来进一步的研究方向.

［1］ 陈晓玲，袁中智，李毓湘等.基于遥感反演结果的悬浮泥沙时空动动态规律研究——以珠江河口及邻近海域为例.武汉大学学报：信息科学版，2005，30（8）：677-681.

［2］ 陈界仁，沙捞·巴里，陈国祥.二维水库水流泥沙数值模拟.河海大学学报，2000，28（5）：11-15.

［3］ 李国胜，王海龙，董 超.黄河入海泥沙输运及沉积过程的数值模拟.地理学报，2005，60（5）：707-716.

［4］ Elias EPL，Cleveringa J，Buijsman MC et al.Field and model data analysis of sand transPort Patterns in Texel Tidal inlet （the Netherlands）.Coastal Engineering，2006，53（5）：505-529.

［5］ Wang ZB，Tse ML，Lau SC.A study on sedimentation of tidal rivers and channels flowing into deeP bay with a Delft3d model.Proceedings of the Sixth International Conference on Asian and Pacific Coasts，2011.

［6］ 陈晓玲，吴忠宜，田礼乔等.水体悬浮泥沙动态监测的遥感反演模型对比分析——以鄱阳湖为例.科技导报，2007，25（6）：19-22.

［7］ 陶 菲，张 鹰，王晶晶等.悬浮泥沙浓度遥感反演模式研究.海洋工程，2007，25（4）：96-101.

［8］ 钱一婧.基于遥感和数值模拟技术的河口水流及水质过程模拟研究［学位论文］.南京：南京师范大学，2008.

［9］ 李 宁.近岸水质的遥感监测和数值模拟研究［学位论文］.天津：天津大学，2009.

［10］ 张 鹏，陆建忠，陈晓玲等.MODIS遥感数据辅助的鄱阳湖水体范围变化数值模拟.武汉大学学报：信息科学版，2012，37（9）：1087-1091.

［11］ Oey LY，Ezer T，Hu CM et al.Baroclinic tidal flows and inundation Processes in Cook Inlet，Alaska：numerical modeling and satellite observations.Ocean Dynamics，2007，57（3）：205-221.

［12］ Feng L，Hu CM，Chen XL et al.Human induced turbidity changes in Lake Poyang between 2000 and 2010：Observations from MODIS.Journal of Geophysical Research，2012，117（C7）.

［13］ 胡春华，施 伟，胡龙飞等.鄱阳湖水利枢纽工程对湖区氮磷营养盐影响的模拟研究.长江流域资源与环境，2012，21（6）：749-755.

［14］ 赖锡军，姜加虎，黄 群等.鄱阳湖二维水动力和水质耦合数值模拟.湖泊科学，2011，23（6）：893-902.DOI 10. 18307/2011.0611.

［15］ 李云良，张 奇，姚 静等.鄱阳湖湖泊流域系统水文水动力联合模拟.湖泊科学，2013，25（2）：227-235.DOI 10. 18307/2013.0208.

［16］ Lin Y，Cao Q.Lake Poyang water quality model for dynamic Prediction.2012 Fourth International Conference on ComPutational and Information Sciences（ICCIS）IEEE，2012：1214-1216.

［17］ 王 鹏，赖格英，黄小兰.鄱阳湖水利枢纽工程对湖泊水位变化影响的模拟.湖泊科学，2014，26（1）：29-36.DOI 10.18307/2014.0104.

［18］ 张 强，张小峰.采砂深度对分汊河道影响分析.华东师范大学学报：自然科学版，2009，（3）：127-135.

［19］ Hibma A，Wang ZB，Stive MJF et al.Modelling imPact of dredging and dumPing in ebb-flood channel systems.Transaction of Tianjin University，2008，14：271-281.

［20］ Jeroen S.The influence of dredging activities on the morPhological develoPment of the Columbia River mouth［Thesis］. Delft University of Technology，DePartment of Hydraulic Engineering Section of Coastal Engineering，2012.

［21］ 李 佳，姚炎明，孙志林等.大型海洋倾倒区悬浮物迁移扩散的数值模拟.浙江大学学报：工学版，2011，45（7）：1319-1328.

［22］ 邬国锋，崔丽娟.基于遥感技术的鄱阳湖采砂对水体透明度的影响.生态学报，2008，28（12）：6113-6120.

［23］ 崔丽娟，翟彦放，邬国锋.鄱阳湖采砂南移扩大影响范围——多源遥感的证据.生态学报，2013，33（11）：3520-3525.

［24］ Leeuw J，Shankman D，Wu GF et al.Strategic assessment of the magnitude and imPacts of sand mining in Lake Poyang，China.Regional Environmental Change，2009，10（2）：95-102.

［25］ 张 鹏.基于遥感的鄱阳湖悬浮泥沙数值模拟与同化研究［学位论文］.武汉：武汉大学，2014.

［26］ Ackerman S，Kathleen S，Menzel P et al.Discriminating clear-sky from cloud with MODIS algorithm theoretical basis document（MOD35）.MODIS Cloud Mask Team，CooPerative Institute for Meteorological Satellite Studies，University of Wisconsin，2010.

［27］ Huet A，Didan K，Miura T et al.Overview of the radiometric and bioPhysical Performance of the MODIS vegetation indices.Remote Sensing of Environment，2002，83（1）：195-213.

［28］ Cui LJ，Qiu Y，Fei T et al.Using remotely sensed susPended sediment concentration variation to imProve management of Lake Poyang，China.Lake and Reservoir Management，2013，29（1）：47-60.

［29］ Roelvink JA，Van Banning G.Design and develoPment of DELFT3D and aPPlication to coastal morPhodynamics.Hydroinformatics，1994，94：451-455.

［30］ 唐立模，肖 洋，周洪都等.赣江东西河分流比影响因素研究.水利水运工程学报，2011，（4）：64-68.

［31］ Roelvink JA，Walstra DJ.KeePing it simPle by using comPlex models.Advances in Hydro-science and Engineering，2004，6：1-11.

［32］ Li J，Tian LQ，Chen XL et al.Remote-sensing monitoring for sPatio-temPoral dynamics of sand dredging activities at Lake Poyang in China.International Journal of Remote Sensing，2014，35（16）：6004-6022.

［33］ PennekamP JGS，EPskamP RJC，Rosenbrand WF et al.Turbidity caused by dredging：viewed in PersPective.Terra et Aqua，1996，64：10-17.

NumericaL simuLation of suspended sediment concentration in Lake Poyang during fLood season considering dredging activities

LI Haijun1，CHEN Xiao1ing1，2，3，LU Jianzhong1**，ZHANG Peng4，QI Hengda1&CHEN Liqiong1
（1：State Key Laboratory of Information Engineering in Surveying，Mapping and Remote Sensing，Wuhan University，Wuhan 430079，P.R.China）
（2：Collaborative Innovation Center of Geospatial Technology，Wuhan 430079，P.R.China）
（3：Key Laboratory of Lake Poyang Wetland and Watershed Research，Ministry of Education，Jiangxi Normal University，Nanchang 330022，P.R.China）
（4：Key Laboratory of Ecological Impacts of Hydraulic-Projects and Restoration of Aquatic Ecosystem of Ministry of Water Resource，Institute of Hydroecology，Ministry of Water Resource and Chinese Academy of Science，Wuhan 430079，P.R.China）

Numerical simulation and remote sensing are used in susPended sediment concentration simulation during flood season in Lake Poyang，in which water turbidity is highly affected by dredging activities.With the aid of Landsat ETM+images，the location of main dredging areas and shiPs are detected.Several sediment Point sources are added in dredging areas in the Delft3D model to simulate the susPended sediment concentration including the effects of dredging of Lake Poyang from July 1，2011 to July 31，2011.Taking in-situ data and susPended sediment concentration retrieved from MODIS images as validation data，there is a good consistence between the simulating and validating results，with R2of 0.831，a root-mean-square error of 15.5 mg/L and a consistent sPatial Pattern of susPended sediment concentration comPared with susPended sediment concentration retrieved from MODIS images.It also shows less dredging effects in the southern main lake and river inlets but more in the northern from channel in the west of Tangyin and the north of Songmen Mountain，water-way in the northern and the lake outlet，which is one of the main reasonscausing high turbidity water in Lake Poyang.This study shows that by making good use of in-situ data and satellite remote sensing images，the Precision of numerical simulation of susPended sediment concentration can be increased by considering the effect of human dredging activities，thus Provided a good method for human affected water environment simulations of lakes.

Numerical simulation；susPended sediment concentration；remote sensing；Delft3D；Lake Poyang

10.18307/2016.0223

*国家自然科学基金项目（41331174，41101415，41301366）、湖北省自然科学基金项目（2015CFB331）、江西省重大生态安全问题监控协同创新中心专项（JXS-EW-08）、测绘地理信息公益性行业科研专项（201512026）和国家科技支撑计划项目（2012BAC06B04）联合资助.2014-12-23收稿；2015-06-08收修改稿.李海军（1992～），男，硕士；E-mail：lihaijun0629@gmail.com.

**

；E-mail：lujzhong@whu.edu.cn.