随机数学方法用于复杂水文过程模拟研究.pdfVIP

随机数学方法用于复杂水文过程模拟 袁永生1 张绍峰2 朱庆平3 (1．河海大学理学院南京 210024； 2．黄河水利委员会水资源保护局 郑州450004； 3．黄河水利委员会 郑州450004) 一 一一…、 ·…一：《《=口qIx====#。wm 【摘 要】 本文系统分析了影响水位表现的水沙要素之后，结合高等的随机分析理论和方法， 寻找出影响水位的显著性非线性相关项，并给出了高精度的相应水位过程拟合，阐明了只用水 沙资料也可有效研究水位表现，用提高理论模型精度的方法达到减小误差的目的。 【关键词】 黄河下游 水文要素 河道地形 随机分析 显著性相关 由予影响因素的复杂性，黄河下游水位表现具有较强的随机性叭2】，因此需要有效寻找 其中的稳定项，使模型误差的方差减小，以提高精度。因此吻合实际问题背景的随机数学。 由于其特殊的工具作用，对实际问题内部特征的研究将更深入、更本质，研究成果的技术 含量和理论水平将得到深化。本文借助高等的随机分析工具，通过研究水、沙、河道等微 观因素作用下的宏观结果，即水沙要素表现，即能实现有效研究黄河下游相应水位过程。 1影响黄河下游水位的水沙因素系统分析 变动河床的洪水位拟合及预报，有着较大的复杂性，但用合理抽象的眼光来看，对问 题的研究可以抽取其本质，从各组成成分作用结果上反演这些成分的作用规律。尽管从相 对宏观角度得出的结论不一定能精确代替微观因素的实际作用规律，但宏观研究的科学性、 合理性将会使这一代替的误差完全在生产中允许的范围内，因此，有效研究水位过程可以 从水沙要素表现入手。 通过对黄河下游水位表现内部结构的研究，结合历史数据的初步分析，上游站某水位 出现时，演进到下游站将表现为多少水位(Zz)，可能受如下水沙要素的影响： (1)上游水位(流量)：Z。Q。。 (2)下游此时刻的(下游同时)水位(流量)：Z，、Q；。Z，(Q，)反映了河段在该时刻 的底水状况，当该时刻前河道断面发生冲刷等外界因素的变化，抬高或降低了z，(Q，)，对 未来的Z，也必然会带来影响。 (3)区间人流量：q。 (4)上游此时的含沙量：S。。从理论上来说，若稍有含沙量，可以认为河道传输是某种 浆体。若含沙量较高时，河槽断面将不断有较大改变。因此含沙量大小与Z，的关系，根据 具体情况应有不同。 156 中国水利学会首届青年科技论坛论文集 (5)上游自最近特征点到该点的水位涨率：z，。这里最近特征点的取法是处于涨水阶段 时，取最近起涨点，落水阶段时取最近峰顶。 (6)上游自最近特征点到该点的水位涨差：z。。涨率、涨差都能反映洪水波的水面比 降，研究中发现，就对zz的影响而言Z。强于z，，由于Z。代表了zr的信息，因此后续研 究中只用z。。 (7)来水来沙系数：K。一S。／Q。。 (8)下游站断面形态：在纯自然条件下，如果仅是局部的冲淤，由于水和沙都是不可 压缩的，同时沙也是充分浸透水的，不会吸附水分子进入沙体，所以总体积不应有太大的 变化，即表现出的水位不应有较大的变化。如果是大规模的侧冲同时侧淤。使原先在水面 以下的沙淤积到水面以上，或水面以上的沙坍塌进水面以下。或者前峰过后，其高含沙或 高流速造成断面的大规模淤积或冲刷，则水位的变化将会较为明显。考虑到河道断面形态 过程资料难以获得，河道的变化因素已较多地体现在z，(Q；)中，本文使用随机分析方法 处理中可见，无须单独考虑。 2数据采集方式 位，均具代表性。所以本文选取花园口与夹河滩1994、1996年汛期实测水沙资料，使用随 机数学作进一步分析。本区间内无人流量。 根据相应水位过程要求的变量对应，这里按下列方式采集数据：①在1994、1996年相 应上、下游水位过程线上，基本上按峰～峰、谷～谷对应，寻找同一特征点在上游出现的 日期、时刻、水位Z。(流量Q。)、含沙值(S。)和演进到下游时相应的水位(Z，)。有些含 沙值没有记录的位置上，用前后有记录位置上的含沙值经线性插值而得。②按采集到的上 游各时刻摘录下游同时水位z：(流量Q，)。若下游同时刻上没有相应资料，仍采用前后值经 线性插值取得。③计算上游水沙系数K。=S。／Q。。采集结果数据见表1、表2。 表1