珠江三角洲近30年来洪潮变化的水文学分析.docxVIP

珠江三角洲近30年来洪潮变化的水文学分析 在未来几十年的海拔高度上升的过程中，珠江三角洲有一大片低地（占平原的23.7%）受到洪水的淹没（李曾等人，1993）。若海平面上升40—60cm,按20年一遇洪水水面线分析,受影响的范围一般距河口60—70km,回水终点以外的上游地区不受影响,引起水位抬高的范围主要是佛山、中山、顺德、江门、新会、东莞六市(范锦春,1994)。实际上,珠江三角洲86.7%的土地有堤围防护,故海平面上升的影响范围不宜与无设防的状况相提并论。而且,不同条件下海平面上升的影响程度是不同的。本文按照7种典型年不同的条件(流量、风暴潮、口门延伸),用水文学及水力学方法,计算各站洪潮水位的升幅,从而划分出海平面上升的不同程度的影响范围。珠江三角洲2030年相对海平面的可能上升幅度为22—33cm(黄镇国等,1999),本文取海平面上升幅度30cm来讨论海平面上升的影响范围。 1 海平面上升和变化的影响区域和影响最大区 计算31个站不同重现期(100年、50年、20年一遇)最高洪潮水位的差值,共得62个数据。若洪潮水位抬高5cm,只有6个数据所代表的重现期缩短,占9.7%;而洪潮水位抬高25cm,则有57个数据所代表的重现期缩短,占92%。海平面上升的影响集中体现在各种潮位的全面抬升及其连带效应,因此,可以考虑以海平面上升后最高洪潮水位升幅5cm等值线作为海平面上升影响范围的边界,此界上游区域称为影响很小区。再以25cm等值线划分影响较大区和影响最大区。 对珠江三角洲共54个水文(水位)站(图1),计算了海平面上升30cm后不同典型年最高洪潮水位的升幅(表1)。 1974年7月22日(74·7典型年,图2a),三角洲上游马口站为中等流量(33108m3/s),沿海为一般风暴潮,代表中水和一般风暴潮的洪潮组合,即代表常年的状况。对47个站,由一维、二维不恒定流水力学模型进行计算。计算时,在网河区布设255个河段, 共截取948个计算面,地形资料大部分采用20世纪80年代以后测绘的 1/5000河道水下地形图。计算结果,等值线的展布如图2a,5cm等值线表明,西北江三角洲西北部属于影响很小区(A区)。25cm等值线表明,东江三角洲大部、广州市区南部、东四口门和西四口门的两侧为影响最大区(C区)。上述两区之间为影响较大区(B区)。 2 潮区锋线与潮区的关系 海平面上升的影响范围随海平面上升幅度、径流、风暴潮、口门延伸等条件的改变而有动态变化。最高洪潮水位升幅5cm等值线可称为“洪区锋线”,其上游区域为洪区,汛期水位完全受洪水流量控制。25cm等值线可称为“潮区锋线”,其下游区域为潮区,台风季节的潮水位主要受风暴潮影响。这两条锋线的位置因径流和风暴潮的强弱对比而变动,多数情况下是一进俱进(向上游),一退俱退(向下游)。 如上述,若以74·7典型年代表常年的状况,则可将下列4种典型年与其对比,反映不同条件下海平面上升影响范围的差异。 2.1 根据风速条件,计算各站标准的海平面上升幅度表2 取马口站20年一遇流量40203m3/s代表一般洪水,在较强风暴潮的条件下,对37个站,由水文学方程计算各站最高洪潮水位(黎子浩,1994): Z2=C1Z1+βQ+αZ1Q+C2Ζ2=C1Ζ1+βQ+αΖ1Q+C2 式中,Z2为各站高潮水位;Z1为海平面上升幅度(以口门外的三灶站为代表);Q为马口站流量;α和β为斜率;C1和C2为常数。α、β、C1、C2的值可以通过相关分析计算而求得。 如图2b,由于洪水条件,径流增强,故5cm等值线向下游推移,即海平面上升的影响范围缩小。因风暴潮较强,故在强潮河口(伶仃洋、崖门),25cm等值线尤其是30cm等值线向上游推移,即影响最大区的范围扩大。 2.29 最大洪授信位设计 1994年6月20日(94·6典型年)发生建国以来西、北江的最大洪水,马口站流量为46300m3/s,接近百年一遇(46736m3/s),代表特大洪水和一般风暴潮的洪潮组合。对48个站的最高洪潮水位进行了水力学计算。升幅等值线(图3)表明,由于特大洪水而风暴潮属一般水平,故25cm等值线明显向口门推移,30cm等值线难以表达,即影响最大区的范围缩小。影响很小区比74·7典型年有所扩大。 2.39 口门不延伸条件下的25cm等值线 1993年9月17日(93·9典型年)马口站流量仅为11800m3/s,而沿海出现特大风暴潮,代表特大风暴潮和小水的洪潮组合。对24个站进行了水力学计算。最高洪潮水位升幅等值线表明(图4a),由于小水年径流弱,故影响很小区向上游退缩,但是,特大风暴使25cm等值线明显向上游推移,即影响最大区大为扩展。 同是93·9典型年,若考虑口门向海延伸,即按珠江三角洲整治规划使口门沿治导线定向延长,延伸长度为:蕉门