按键盘上方向键 ← 或 → 可快速上下翻页，按键盘上的 Enter 键可回到本书目录页，按键盘上方向键 ↑ 可回到本页顶部！
————未阅读完？加入书签已便下次继续阅读！
关于长江三峡砾卵石输移量的讨论 …………………………………… 240
关于长江三峡砾卵石输移量的讨论（续） …………………………… 250
吁请长江三峡大坝即日停工！此坝决不可修！ ……………………… 265
关于长江三峡修建高坝可行性问题 …………………………………… 267
诗文拾零
花丛小语 ………………………………………………………………… 273
《治水吟草》自序 ……………………………………………………… 277
右冠残草 30 首 ………………………………………………………… 279
治河咏怀 15 首 ………………………………………………………… 289
忆旧感怀 26 首 ………………………………………………………… 299
漫游闲咏 29 首 ………………………………………………………… 308
许伯伦先生墓志铭 ……………………………………………………… 318
敬和克木先生惠诗 ……………………………………………………… 319
魏明初先生墓志铭 ……………………………………………………… 320
先师罗公建侯讳忠忱廿年祭 …………………………………………… 321
哭长江三峡大坝开工 …………………………………………………… 323
黄万里自述 ……………………………………………………………… 324
黄万里简历 ……………………………………………………………… 326
记者访谈录
不重水文，何来水利？（宜宏） ………………………………………… 327
黄万里：洪灾过后谈治水 （伍诗一） ………………………………… 330
“右派教授和他的江河缘”…访清华大学水利系教授黄万里先生… 333
但教莫绝广陵散—记水利专家黄万里 （赵诚）………………………… 340
补遗
水利工程学理论
沙流连续方程意义的简释 *
目的：本文通过沙流连续方程这一课题，企图说明近代水沙合体流（固体和
液体混合流 Two Phase Flow）力学和长流程河工模型试验理论还很幼稚，甚至还
没有能成立。内容采用了浅易的方法解释，为了帮助年青同志普遍地理解；也请 大家指正，提出不同意见，公开讨论。
沙流连续方程的释义
在水沙流动力学中，在水沙流模型理论中，沙流的连续方程普遍地采用了拉 波拉斯 Laplace 方程：
？qs
？x
＋ ？ ？
？H b = 0
？t
（1）
见 Великанов，ДинамикаРусловыхПомоков Ⅱ304 页 Ливи；ДинамикаРусловыхПомоков 243 页
Yalin； Theory of Hydraulic Modeds 183 页
《河流动力学》武汉水利学院讲义，及窦国仁讲义等。
式中 qs 为单宽输沙率，Hb 河床底面平均高程，？ ？ 床沙带有空隙的干容重，x 流程，
t 时间。
这一方程假设 Q 和水流中含沙浓度在各断面上并不随时变化，是错误的；它
在中外至今普遍地应用，产生了并正不断地产生着不良的结果。
这是水动力学中最基本最简单的第一个方程。下面作简明的解释。
这个方程是仿照了水流或水沙混合流的连续方程而成立的。公式
* 讲义，1975 年初稿，1977 年重写。
1
？q ？h
＋ = 0
（2）
？x ？t
（2）对水流是正确的，但准确的沙流方程应是泊桑 Poisson 连续方程：
？qs ＋ ？ ？ ？H b = ？ =？（sh）
（3）
？x ？t ？t
（式中 q — 单宽水流率，h 水深）其意义无非是，根据物质不生不灭的公理：当
？qs
一定时刻 t，在 dx 流程内，输沙率的沿程增值 dx 和河床在同一地点堆高起来
？x
的沙量增率 ？ ？ ？H b dx ，一定等于这一时刻、这一地点来自水流中挟沙率的减值
？t
？（sh）
？ dx 。（式中 s 为含沙浓度）。在现行沙流连续方程（1）中，却把末项假设
？t
为 0：
？（sh）
= 0
？t
？qs
也就 是假 设输 沙率 qs 的沿程增 值 dx 是全部 从河 床上 冲起 来的 沙量增 率
？x
dx
？ ？ ？ ？H b dx ，而水中的挟沙率 =？（sh） 则为 0，或 sh 始终随时不变。这不符合任
？t ？t
何河中水沙流的实际情况，是没有根据的。
？H b
试想，当河床淤高时（
？t
》 0 ），沙是从哪里来的？不是先得从河床上的水
流含沙 sh 里供给的吗？那么怎么可能 sh 本身却不也随时而变呢？怎可假设完全
？H b
来自上游多送了的沙所淤的呢？又当河床冲深时
？t
《 0 ，当然冲起来的沙首先
去加大其顶头上的 sh，然后再靠水流把它部分输送下去，怎可假设全部冲起来的
沙恰恰都输往下游呢？作者曾绘过黄河、渭河、丹江上许多 Q（流率），V（流速），
s，Qs ，Hb 等的时程线，把它们套在同一 t 轴（见附图 1），没有一张图不显示 s
？s
是随着 Q，Qs ，V，Hb 等一起改变着的，而且其变率 特殊地大，决不能把它假
？t
设为 0。
2
由图可见：
？？ ？？ ？Q ？V
（1） ？ 0 ， ~ ， ，……（后段 ？ ~ t 线可能因 ？ 未测准，有问题。）
？t ？t ？t ？t
（2）当 H，Q，V 增时， ？ = s ， d 50 跟着加大，而 H b 减小，即河底冲深；
当 H，Q，V 增时， ？ = s ， d 50 跟着减小，而 H b 加大，即河底淤高；
（3）洪水发时，粗沙冲起；洪水过后，细沙铺面。
图 1 渭河临潼站各水力因子的变化
（1966 年 8 月 31 日~9月 10 日）
3
从实测资料的图线里，可以看到，只有在枯水期或汛期的尾水段，在近乎恒
定流的情形下，当 Q，Qs ，Hb 变化极小时，s 才可能也变化较小，但其时河床也
根本并没有显著的冲淤。
约一个世纪以来，数学家们致力于解出拉波拉斯的偏微分方程，作了一些简
化的假设就可能解出。但人们没有去考虑这种假设所引起的脱离实际情形的错误。 下面是作者以前分析沙流连续方程的结果：
s
按 q = qs ， q = Vh ， V = dx ，代入式（3）得
dt
？s ？q ？H ？h ？s
q ＋ s ＋ ？ ？ b = ？s ？ h
？x ？x ？t
？t ？t
？s
移项， Vh
？s ？H
？
＋ h ＋ ？ ？ b
= ？s？ ？q
？h ？
＋ ？
？x ？t ？t
？ ？x
？t ？
？q ？h ？x
按式（2）， ＋ = 0 ，双代入 V =
？x ？t ？t
h
？ ？s ？s dx
乃得 ＋
= ？？ ？ ？H b
？？ ？t
？x dt ？？ ？t
h b
ds ？H
= ？？ ？
（4）
dt ？t
这一简单方程的物理概念可以这样解释：令 b = 流宽
？ ？ ？ ？
b × h ds × dx = ？b？？
H b ？dx
dt ？
dt ？
或 （b × h × ？x）？
x；t
t
b
s = ？（b × ？ H
× ？x）？ ？
（5）
此式的意义是，当一定时刻含沙浓度 s 在水体 （bh？x）里沿程 ？x 过了 ？t 时间所增
加的输沙重量微量 （bh？x） ？
x ；t
s 一定就是从这一地点 x 处河床里所冲起来的床沙
t
b
微量 ？（b ？ ？ H
？ ？x） 乘以其具有空隙的床沙容重 ？ ？ 所得的沙重增微量
t
b
？（b； ？ H
？ ？x）？ ？ 。换句话说，就是，当一定时刻 t 从一定地段 ？x 在相应时段 ？t 冲
起来的沙重增微量，就是当时当地水流中本身所加浓的沙重，以及沙重沿程随时
增加的输沙微量；否则沙到哪里去了呢？或沙从哪里来的呢？下面再作一图例说
明，使读者心中完全明确，目前流行的沙流连续方程是极其错误的，它错误地引
4
导了连续方程在水动力学中的应用途径，连带模型比尺的理论一起导入了错误的
途径。
设在某段水沙流中水宽 1 米，水深 h = 1 米，每粒沙的直径和容重都是一样
的，如图 2 所示。设开始水沙流是恒定的、均匀的，在时段 ？t 内有 4 粒单位沙
（q ？t = 4） 通过断面。接着河底平均冲深了 ？ H
= ？0。03 米。在河底上每米河宽、
s t b
沿程 ？x 的平原上原来铺着 100 粒单位沙。（图中画了 10 粒，每米河宽有 10 行）
根据式（5）
（1×1× ？x）？
x；t
s = ？（1× （0。03）× ？x）×100
因此有 ？ x ；t ，s = 3 粒单位沙冲了起来。图中画的是 2 粒冲到顶上水里加浓了水中
？（hs ）
？ ？q ？
含沙浓度，= ？t？x = 2 粒单位沙；另外 1 粒顺流而下，就是 ？q ？t； ？ s dxdt ？ = 1
？t
粒单位沙。
s ？ ？x ？
其次，我们来讨论连续方程在水动力学中的作用。这方程仅仅是水沙质量在
运 动中的 守恒 性质， 还没 有触及 力和 能的守 恒分 析，后 者将 依靠运 动方 程 。
？ ？ ？ ？H b dtdx = ？？ ？（1 ？ ？ H？x） 代表在床面 b × ？x（b = 1） 上在一定时间 ？t 内冲起的
？t t
5
？ ？ ？ ？
沙重。 h ds dtdx = h？ s ＋ s dx ？dtdx ？ （？ s ）＇1× h？x＇ ＋ （？
s ）＇1× h？x＇，代表这床面
dt ？ ？t
？x dt ？ t x
上 h 水深中同时间内加浓了沙重及当水沙流穿过这一体积时增多了的沙重。当此
式和运动方程联解时，若忽略掉最后这一项，结果就不合理了。 沙流连续方程在模型理论中的应用
拿式（4）、（5）和错误的式（1）比较，模比的关系就大大分歧了。按式（1） 得出错误的模比关系是
s
？ = ？t ？q
？H b
？？ ？s ？ x
（见 Yalin 183 页） （6）
在长江模型中则是
？t =
？ x ？ ？？ ？
？1 2 ？ ？