4. 浮遊砂と河床変勤

荒井信行・清水康行

 第 3 章で述べたように, 浮遊砂は乱流拡散現象によって輸送される形式の流砂である。 浮遊砂を含む流れは, 例えばダム貯水池の堆砂や河口部浚渫個所などの洪水による多量の土砂堆積等に見られるように, 河床変動に大きな影響を与え,実際上も重要な問題であるため, これまで浮遊砂に関して数多くの研究がなされてきた。 それらの研究は,大別すると浮遊砂流の特性に関するものと土砂の輸送量に関するものがある。
 今回は,主に均一粒径砂を対象とした浮遊砂量式と浮遊砂による 1 次元単断面の河床変動の問題に重点をおいて述べることにする。
 なお,浮遊砂に関する諸問題は,複雑, かつ広範多岐にわたるため, 上記の研究成果は多くの仮定や近似的取扱いに基づいており, 未解明の問題も多く残されているのが現状である。 このような状況から,本章では実際間題として一般的な適用性を失なわない程度に問題を簡略化して解説するつもりである。

4.1 浮遊砂の基準点濃度と砂粒の浮上量式

 浮遊砂の濃度 c は,容積表示で次のように定義される。 いま,流量を Q (m3/s), その流量中の浮遊砂量を Qs (m3/s) とすると, c=Qs/Q である。 図-4.1 のように,定常 2 次元等流を考えた場合の単位幅当たりの浮遊砂量 qs (m3/s/m) は, 濃度分布 c(z) および流速分 u(z) が決定されれば, 次式で与えられる。

(1)qs=bhc(z)u(z)dz

 ここに,h:水深, c(z):河床面上 z の点における濃度, u(z):河床面上 z の点における流速, b:濃度の基準点とされる河床面からの高さであり, b のとり方により濃度分布が若干異なってくる。

_images/fig_6.1.png


図-4.1 流速分布 u(z) および濃度分布 c(z) の定義図

 浮遊砂量の算定は, 濃度の基準点 b の値とその点の濃度 cb の算定の問題に帰着する。 一般的には,b の値は水深の 5 %,すなわち b=0.05h が採用されている。 また,浮遊砂の基準点濃度および砂粒の浮上量式としては, 以下の板倉の式1)を紹介する。 板倉は,河床が平担な場合の浮遊砂を含む開水路流れについて, 流速分布,浮遊砂の濃度分布,浮遊砂の基準点濃度,浮遊砂流の抵抗特性, 浮遊砂量式,堆積限界を一貫した理論によって系統的に導いている。 詳細は原論文を参照されたい。
 いま,河床が均一砂粒で構成されている場合を考える。 河床が平衡状態にあるとき, 河床の単位面積から単位時間に巻き上がる粒子群の体積 qsu (m/s) は, 単位面積当たり単位時間に河床へ沈降する粒子群の体積 wfcb (m/s) に等しいので,次式がなりたつ。

(2)qsu=wfcb

 ここで,浮遊砂の底面付近 (z=b) の濃度, すなわち基準点濃度 cb は,

(3)cb=K(αρsρρsgduwfΩ1)

と書き表わされ,qsu は次式により与えられる。

(4)qsusgd=K(αρρsΩτwfsgd)
(5)Ω=τBaξ1πexp(ξ2)dξa1πexp(ξ2)dξ+τBη01
a=Bτ1η0,η0=0.5,α=0.14,K=0.008,B=0.143

 ただし, ρs:砂粒の密度,ρ:水の密度,s:砂の水中比重, g:重力加速度,d:砂粒の粒径,wf:砂粒の沈降速度
 また,(4)式の近似式として次式がある。

(6)qsusgd=0.008[0.14ρρs(14τ0.9τ)wfsgd]

4.2 浮遊砂濃度の連続式と河床変動の連統式

 問題を簡単にするために, 以後浮遊砂の濃度 c として次式で定義される水深方向の平均濃度 c¯ を用いて議論を進める。 すなわち,c=c¯ とする。

(7)c¯=qsuq=bhc(z)u(z)dzbhu(z)dz

 まず,浮遊砂濃度の連続式を誘導してみよう。 ただし,本節では濃度 c はある時間 Δt の間では一定として扱い, 拡散作用を無視する。

_images/fig_6.2.png


図-4.2 浮遊砂濃度の連続式

 いま,図-4.2 に示すように主流方向に x 軸をとると, 微小時間 Δt の間に断面 ① から入ってくる浮遊砂量は, Qs1Δt=Qc1Δt であり, 同様に Δt の間に断面 ② からでていく浮遊砂量は, Qs2Δt=Qc2Δt である。 次に,断面 ①,② 間で Δt の間に浮上する砂量および沈降する砂量は, それぞれ qsuBΔxΔt および wfcBΔxΔt である。
 したがって,断面 ①,② 間でこれらの砂量の収支が 0 となるためには,

(8)Qc1ΔtQc2Δt+qsuBΔxΔtwfcBΔxΔt=0

 ゆえに,

Q(c2c1)BΔx=qsuwfc

となる。上式を微分方程式で表わすと,

(9)QBdcdx=qsuwfc

となる。(9)式が浮遊砂濃度の連続式である。
 ここで,特別な場合として流下方向に濃度が一定 (平衡状態), すなわち dc/dx=0 とすると,(9)式の右辺 = 0 であるから, c=qsu/wf となる。
 実際に(9)式により流下方向の濃度分布を計算する場合には, 左辺の dc/dx の差分形が問題となる。 いま,(9)式を後進差分で計算することを考えてみる。
(9)式の差分式は,

QB(x)c(x)c(xΔx)Δx=qsu(x)wfc(x)

であるから,c(x) について整理すると,次式となる。

(10)c(x)=qsu(x)+c(xΔx)QB(x)ΔxQB(x)Δx+wf

 したがって,各断面における砂粒の浮上量 qsu(x), 砂粒の沈降速度 wf が求まり, 境界条件として上流端の濃度 c(0) が与えられれば, (10)式により各断面の濃度 c(x) を計算することができる。
 次に,浮遊砂による河床変動の連統式を誘導する。

_images/fig_6.3.png


図-4.3 浮遊砂による河床変動の連続式

図-4.3 に示すように, ある時刻 t から微少時間 Δt の間に河床から浮上する砂量および河床へ沈降する砂量は, それぞれ qsuBΔxΔt および wfcBΔxΔt である。 いま河床が低下する場合を考えると,流水中の砂粒の増加量は,

(qsuwfc)BΔxΔt

であり,これは河床から流水中にでていく砂量,すなわち河床低下量に等しい。
 河床低下量は,河床変動高を Δz, 空隙率を λ とすると,

ΔzBΔx(1λ)

で表わされるから,上記の 2 式を等しいとおくと,次式が得られる。

ΔzΔt=11λ(qsuwfc)

 上式を微分方程式で表わすと,

(11)dzdt=11λ(qsuwfc)

となる。 (11)式が浮遊砂による河床変動の連続式である。

4.3 非定常性を考慮した浮遊砂濃度の連続式

 前出の浮遊砂濃度の連統式(9)式は, 微小時間 Δt の間で濃度 c を一定として取扱った場合の式である。 しかしながら,実際問題では多くの場合濃度 c が時間的に一定とは考えにくく, 例えば演習問題 18で示すように, ダム貯水池へ浮遊砂が流入するような場合には, 前記の計算法では初期条件, すなわち時刻 t=0 の段階で全断面で濃度 c があることとなるのに対して, 実際には時刻 t=0 ではすべての断面で c=0, すなわち清水という条件があり得る。
 そこで,本節では濃度 c の時間変化を考慮した浮遊砂濃度の運続式を誘導する。 前出の 図-4.2 より,ある時刻 t における濃度を c(t), Δt 時間後の濃度を c(t+Δt) とすると, Δt 時間の間の断面 ①,② 間における流水中の浮遊砂の増加量は, 次のように表わされる。

c(t+Δt)BhΔxc(t)BhΔx

 ここで,濃度 c の変化量を Δc=c(t+Δt)c(t) とおくと,上式は,

ΔcBhΔx

となる。
 断面 ①,② の間の砂量の収支を考えると,前出の(8)式の左辺が ΔcBhΔx に等しいことから, 次式が得られる。

(12)Qc1ΔtQc2Δt+qsuBΔxΔtwfcBΔxΔt=ΔcBhΔx

(12)式の両辺を BΔxΔt で割ることにより,次式が導かれる。

hΔcΔt+QBc2c1Δx=qsuwfc

 上式を微分方程式で表わすと,

(13)hct+QBcx=qsuwfc

となり,上式が非定常性を考慮した浮遊砂濃度の連続式である。
 ここで,(13)式を時間 t,距離 x に関して差分表示することを考える。
 まず,(13)式より Δt 時間後の濃度 c(t+Δt) を求める式として,次式が得られる。

(14)c(t+Δt)=c(t)+Δth(qsuwfcQBcx)

(14)式により,時刻 t における流下方向の濃度分布から Δt 時間後の濃度分布を求める際に問題となるのは, c/x の差分形のとり方である。
 いま,後進差分の場合を考えると, (14)式における濃度 cxt の関数なのでより丁寧に書くと,(14)式は次式となる。

(15)c(x,t+Δt)=c(x,t)+Δth{qsuwfcQBc(x,t)c(xΔx,t)Δx}

(15)式が求める(13)式の差分式である。
 また,定常等流状態の流れを考えると, (13)式は次式のように変形される。 (13)式の両辺を水深 h で割り, x 軸方向の平均流速 u=Q/Bh (一定) とおくと,

(16)ct+ucx+wfch=qsuh

となる。(16)式を移流方程式と呼んでいる。
 一般に,(16)式のような形式の偏微分方程式においては, 特性曲線は dx/dt=u となり, c/x の差分形式は流速 u と同方向, すなわち流下方向に向いた差分形である後進差分とする必要がある。

4.4 拡散項を含む浮遊砂濃度の連統式

図-4.4 において,断面 ① における濃度 c1, 断面 ② における濃度を c2 とすると, (c2c1)/Δx=dc/dx を濃度勾配という。
 一般に,2 地点間に濃度勾配があるとき, すなわち 2 地点間で濃度に差があるときには, 濃度勾配に比例した物質移動が,その差をなくす方向に生ずる。 その移動量は,比例定数を D とすると Ddc/dx と表わされる。 この D を拡散係数という。

_images/fig_6.4.png


図-4.4 濃度勾配 (矢印は浮遊粒子の移動方向)

_images/fig_6.5.png


図-4.5 拡散作用による浮遊砂の濃度変化

 いま,図-4.5に示すように流速ゼロの状態, すなわち静水中を考える。
 断面 ① から微小時間 Δt の間に拡散作用により流入する浮遊砂量は,

D(dcdx)1hBΔt

であり,断面 ② から Δt 時間に流出する浮遊砂量は,

D(dcdx)2hBΔt

である。一方,断面 ① ,② の間で Δt 時間に濃度が c(t) から c(t+Δt) に変化するとすれば, その変化量は,

{c(t+Δt)c(t)}hBΔx

で与えられる。ここで,断面 ①,② 間の砂量の収支を考えると,

{c(t+Δt)c(t)}hBΔx=D(dcdx)2hBΔtD(dcdx)1hBΔt

 ゆえに,

c(t+Δt)c(t)Δt=D(dcdx)2(dcdx)1Δx

 上式を微分方程式で表わすと,

(17)ct=D2cx2

となる。(17)式を拡散方程式という。
 したがって,実際河川における浮遊砂濃度の連統式は,(16)式と(17)式をあわせて考えると, 次式により表わされる。

(18)ct+ucx+wfch=D2cx2+qsuh

4.5 水深方向の浮遊砂濃度分布

 これまでは,浮遊砂量 Qs を流量 Q と水深方向の平均濃度 c(=c¯) の積として求めてきたが, 実際には 第 1 節で述べたように, 濃度 c は水深方向の分布を考える必要がある。 例えば,前出の 図-4.2 のように河床からの砂の出入りを考える場合の沈降量 wfc における c は, 平均濃度 c¯ ではなく,河床底面付近の濃度, すなわち基準点濃度 cb を用いなければならない。
 いま,定常 2 次元等流状態の流れを考え, 流下方向には濃度変化はないものとする。

_images/fig_6.6.png


図-4.6 水深方向の浮遊砂の濃度変化

図-4.6 に示したような微小要素 ΔzΔx 内に, 微小時間 Δt 間に拡散作用により流入する浮遊砂量は,

(Dcz)2ΔxΔt

であり,同様に流出する浮遊砂量は

(Dcz)1ΔxΔt

である。
 また,微小要素 ΔzΔx 内に Δt 間に流入する沈降量は,

wfc2ΔxΔt

であり,微小要素 ΔzΔx から流出する沈降量は,

wfc1ΔxΔt

である。
 以上の砂量の収支を考えると,

{(Dcz)2(Dcz)1}+wf(c2c1)=0

となるので,両辺を Δz で割って微分方程式で表示すると, 次式が得られる。

(19)Dddz(dcdz)+wfdcdz=0

 通常,平担河床では(19)式中の拡散係数 D は渦動粘性係数 ε に等しく考える。 この ε を次式で定義する。

(20)ε=κ6uh

 ただし,κ:カルマン定数, h:水深,u:摩擦速度
 したがって,(19)式は次式に書換えられる。

(21)εddz(dcdz)+wfdcdz=0

 上式を z で積分すると,

εdcdz+wfc+A1=0

であり,積分定数 A1 は水面における条件, すなわち水面 z=h では水面から浮遊砂粒子が飛びだすことはないので, dc/dz=0wfc=0 であるから,A1=0 となる。
 ゆえに,上式は次式で表わされる。

(22)εdcdz+wfc=0

 上式を積分すると,

dcc=(wfε)dz

 ゆえに,

lnc=wfεz+A2

 または,

c=A3exp(wfεz)

となる。ここで,基準点濃度 cb を近似的に河床面上 (z0) の濃度と考えると, A3=cb であるから, 濃度 c は次式で与えられることになる。

(23)c=cbexp(wfεz)

 上式が浮遊砂の濃度分布式である。
 これに対して, 前節までで使用していた濃度 c は河床から水面までの平均濃度である。 (23)式により平均濃度 c¯ を求めると,

c¯=1h0hcdz=1h0hcbexp(wfεz)dz=cbh[εwfexp(wfεz)]0h=εcbwfh{exp(wfεh)exp(wfε0)}=εcbwfh{1exp(wfhε)}

 ここで,β=wfh/ε とおくと, 次式が得られる。

(24)c¯=cbβ(1eβ)

 上式が平均濃度 c¯ と基準点濃度 cb の関係式である。 したがって,上式より基準点濃度は次式のように書き表わされる。

(25)cb=β1eβc¯

4.6 実際河川の浮遊砂による河床変動計算法

 前節で述べたように, これまで河床への沈降量は wfc としてきたが, 実際には wfcb であるので, (18)式および(25)式より実際河川における浮遊砂濃度の連続式として最終的に次式が得られる。

(26)ct+ucx+wfhβ1eβc=D2cx2+qsuh

 既述のように,上式で拡散項を無視すると,次式が得られ,

(27)ct+ucx+wfhβ1eβc=qsuh

 さらに,濃度の時間変化を無視すると,次式となる。

(28)ucx+wfhβ1eβc=qsuh

 また,浮遊砂による 1 次元の河床変動を考える場合の河床変動の連続式(11)式も,wfc の代わりに wfcb とする必要がある。 したがって,浮遊砂による河床変動の連統式は次式で与えられることになる。

(29)dzdt=11λ(qsuwfβ1eβc)

 以上で, 浮遊砂によるI 次元の河床変動計算を行う準備が整ったので, ここでは河床変動の基礎式を差分法により解く手順を示す。 図-4.7 は浮遊砂による河床変動計算のフローである。

  1. 現河床形について,与えられた流量条件のもとに不等流計算を行い,掃流力の縦断分布を求める。

  2. 砂粒の浮上量式を用いて,浮上粒子量の縦断分布を求める。

  3. 浮遊砂濃度の連続式により,浮遊砂濃度の縦断分布を求める。

  4. 河床変動の連統式により,河床変動量を求める。

 以上のプロセスを繰り返し行うことにより, 河床高の時間変化を予測することができる。
 差分法により安定な数値解を得るためには, 差分スキームの与え方,差分幅 ΔxΔt のとり方, 境界条件の設定に注意が必要である。これらの点については,以下の 演習問題をとおして解説することにする。

_images/fig_6.7.svg


図-4.7 浮遊砂による河床変動計算のフ口ー

演習問題 16

 砂粒の粒径 d = 0.05, 0.1, 1 mm の場合について, おのおの無次元掃流力 τ = 0.1, 1, 2, 5, 10 の場合の砂粒の浮上量 qsu(4)式および(6)式で計算し,計算結果を比較せよ。

解答例

演習問題 17

 前出の演習問題 3で示した河川において, 河床砂の粒径 d = 0.1 mm とした場合について,以下の設問に答えよ。

  1. 砂粒の沈降速度 wf を Rubey の式により求めよ。

  2. 流量 Q = 4000 m3/s,下流端水位 H = 5 m のとき,水面形を求めよ。ただし,Manning の粗度係数 n = 0.025 とする。

  3. 各断面の摩擦速度 u および無次元掃流力 τ を求めよ。

  4. 各断面における砂粒の流砂形態を判定せよ。

  5. 各断面の砂粒の浮上量 qsu を板倉の式(4)式,あるいは(6)式により求めよ。

  6. 各断面の浮遊砂濃度 c(10)式 (後進差分) により求めよ。ただし,境界条件として与える c(9)式の右辺 = 0 として得られる c=qsu/wf とする。

    1. で得られた c と各断面で c=qsu/wf として計算した c を比較せよ。

解答例

演習問題 18

 河床勾配 ib = 1/700,河幅 B = 100 m, Manningの粗度係数 n = 0.02,河床砂の粒径 d = 0.1 mm の広矩形断面水路に流量 Q = 500 m3/s が流下する場合について, 以下の設問に答えよ。 ただし,下流端の水位 H = 12 m, 同河床高 Z = 0 m (標高) とし,計算延長 L = 10 km, 計算区間距離 Δx = 200 m とする。

  1. 初期断面における水面形を求めよ。

  2. 各断面の浮遊砂濃度 c(10)式で求めよ。ただし,差分計算は後進差分で行い,上流端の濃度 c0=qsu/wf とせよ。

  3. (11)式を用いて,各断面の 10 秒後の河床変動高 Δz を後進差分で求めよ。ただし,境界条件として上流端の z = 0 を与えよ。

  4. 計算時間間隔 Δt = 10 秒として,24 時間分の河床変動の過程を縦断図に示せ。なお,縦断図には 1, 4, 12, 24 時間後の河床形状を描くこと。

解答例

演習問題 19

演習問題 18の条件で, 浮遊砂濃度 c の計算を(13)式で行った場合の河床変動を計算せよ。

解答例

演習問題 20

 前出の演習問題 18の条件で, 下記の 3 ケースについて河床変動計算を行ったときの計算法,差分法,プログラムおよび計算結果を示せ。 ただし,拡散係数 D は 0.001 とする。

  1. 浮遊砂の濃度 c の計算を(26)式,河床変動高の算定を(29)式で行った場合

  2. 浮遊砂の濃度 c の計算を(27)式,河床変動高の算定を(29)式で行った場合

  3. 浮遊砂の濃度 c の計算を(28)式,河床変動高の算定を(29)式で行った場合

解答例


参考文献

1) 板倉忠興:河川における乱流拡散現象に関する研究,土木試験所報告,第 83 号,昭和 59 年 7 月