GIS支持下的黄河口近期淤、蚀动态研究

　　中图分类号：P208/TV148 文献标识码：A

　　黄河以泥沙含量高而著称于世，含沙量高达36kg/m[3][1]。自从1855年铜瓦厢决口夺大清河道注入渤海，行水百余年，尾闾河道改道十余次，形成了现今由10个三角州舌体组成的黄河三角洲[2～4]。现河道是1976年5月由清水沟入海，至今已行水22a。一般认为，黄河三角州每年造陆速率达20余km[2][5～7]，但是，在黄河三角州的形成和发育过程中也遭受着海洋潮流的侵蚀，造成陆地面积减少。由此可见，黄河三角州是黄河泥沙淤积与海洋潮流侵蚀相互作用的结果。本文采用遥感和地理信息系统(GIS)技术，研究了黄河口1986年5月至1996年10月十余年间的淤积与侵蚀面积变化及其空间动态，旨在准确掌握黄河口淤、蚀动态变化及趋势，为长期监测打下基础。

　　1 供试资料与监测范围

　　为了保证研究结果准确可靠，全部采用不同时相的黄河口卫星图像资料（表1），以目视解译方法获得各时相黄河口岸线图。通过岸线变化的对比分析，确定的研究监测范围为：以黄河故道与现河道交汇处为项点，北部以黄河大坝为界，南部以垦东水库西至垦利县养虾池、小岛河，形成一个小的扇形区域，总面积1106.86km[2](1996.10)，从入海口至最北端约48km，至最南端约33km。

　　

　　2 监测方法

　　水陆界线在卫星像片上有较清晰的反映，首先以目视解译方法获得了黄河口各时相的岸线解译图，采用PC ARC/INFO GIS软件，将各时相黄河口目视解译图手扶数字化输入，不同时相图赋不同的用户标识码(USER-ID)，以便显示识别，将图面数字化仪座标转换为统一的实地方里网座标，经编辑处理形成黄河口图形库。在ARC/INFO支持下将不同时相的解译图按确定的图幅控制点（TIC点）进行准确叠加，通过分析显示其空间变化图并从其属性库中自动统计出各时期的淤积和侵蚀面积[8]。

　　3 结果分析

　　3.1 黄河口淤、蚀面积变化分析

　　3.1.1 研究区域总面积变化分析

　　表2为1986年5月至1996年10月5个时相的黄河口总面积。可以看出，从1986年5月至1996年10月十余年间，总面积总体上是增加的，但并非每年递增。以1990年1月分为前后两个阶段，前段总面积为显著增加，增长了91.51km[2]，而后段总面积却减少了21.16km[2]，从1995年3月至1996年10月，总面积有所增加，但仍未达到1990年1月的总面积水平。

　　

　　3.1.2 黄河口不同时段淤、蚀面积变化分析

　　从各不同时段黄河口面积的增减变化情况分析（表3，图1，此面积为淤蚀陆地面积，不含其中黄河本身水域面积），可以看出：①从1986年5月至1996年10月的10a5个月中，黄河淤积造陆面积是95.14km[2]，而被侵蚀面积是26.44km[2]，净增加面积为68.70km[2]，平均年净增6.59km[2]；②从各时段淤积面积变化情况看，以1988年12月至1990年1月时段淤积面积最大，平均44.58km[2]/a，其次是1986年5月至1988年12月时段和1995年3月至1996年10月时段，分别为26.46km[2]/a和21.56km[2]/a，1990年1月至1995年3月时段最小，为7.56km[2]；③从各时段侵蚀面积变化分析，以1995年3月至1996年10月时段最多，达17.72km[2]/a，其它时段相差不大；④不同时段净增减面积差异较大，1990年1月以前面积为净增加，平均24.93km[2]/a，而1990年1月以后面积转为净减少，平均年递减2.40km[2]。也就是说，黄河口1990年1月以前以淤积为主，净面积增加，而此后则以侵蚀为主，净面积减少。

　　3.2 黄河口淤、蚀变化的空间分析

　　根据黄河口不同部位的淤、蚀状况，可分为四个不同的淤蚀部位：河口北部、河口嘴部、河口南部以及自垦东水库部位沿防潮坝向南急转至垦利县养虾池的河口西南部（图2），1986年5月至1988年12月，河口的主要变化是由短粗变为细长。淤积主要在河口嘴部和与河道近直交的河口西南部，而侵蚀则集中在河口北部，使河嘴形状由钝园变为尖细；1988年12月至1990年1月，河口北部和西南部受到轻微的侵蚀，变化不大，泥沙则集中沉积于河口南部，使河道下部（南侧）加厚；1990年1月至1995年3月，除嘴部延伸面积增加外，河口北部和西南部均受到较为强烈的侵蚀，使河口嘴形更显细长；至1996年10月黄河泥沙改由靠近嘴部斜向东北的一条支流入海，泥沙淤积在1990年1月至1995年3月时段受到较大侵蚀的部位，从而形成十分明显的凸起。同时，由于河水的转向，河口南部受到较大的侵蚀，此时段侵蚀面积较大的原因与黄河入海流路的改变有直接关系。