土木学会論文集B1(水工学) Vol.74, No.5, I_919-I_924, 2018.

石礫粒子形状の違いが流れ場と河床形状に及ぼす影響の研究

熱海 孝寿 1・福岡
1 学生会員中央大学大学院

捷二 2理工学研究科

都市人間環境学専攻（〒112-8551 東京都文京区春日 1-13-27）
E-mail: a14.c3y6@g.chuo-u.ac.jp
2 フェロー

中央大学 研究開発機構教授 （同上）
E-mail: sfuku@tamacc.chuo-u.ac.jp

石礫河川の流砂現象には粒子形状の影響が大きく，流れと流砂運動に対する影響を明らかにすることは重要である．本研究では，石礫河川の粒子形状に着目し，異なる形状を持つ粒径集団からなる石礫群を用いた数値移動床実験によって主流部と堆積層内の流れ構造および相互作用，また，河床形状，粒子移動量に及ぼす粒子形状の影響について調べた．その結果，主流部の流れは河床表層から水深程度の深さの範囲で，堆積層内の流れに影響を与えること，そして，堆積層内の流れが主流の乱れ構造に影響を与えることを明らかにした．さらに，異なる粒子形状が河床変動の機構に違いをもたらす．過程を詳細に分析し，流砂量，河床の形状，河床変動に及ぼす影響を明らかにしている．

Key Words: gravel bed river, particle shapes, permeable gravel bed , numerical movable-bed channel,
turbulence characteristics, bed geometry

1. 序論を用いた田所ら 6)の検討では，移動する石礫群の動摩擦角の違いが粒子群の運動形態に影響を与え，流砂量に大石礫河川を構成する河床材料は様々な粒子形状の幅のきな違いを生じさせることを示した．しかし，これらの広い粒度分布を有しており，洪水時には，土粒子の間欠研究は，形や大きさの異なる粒径集団による数値移動床的な運動や分級，覆瓦構造等の特徴的な現象が見られる． 実験に基づくため，沢山の粒子が互いに影響し合い，粒石礫河川の土砂移動現象を把握するため，福岡 1)は常願

子移動の結果に及ぼす形の影響が分かりにくい面があっ

寺川の高水敷を掘削して作られた大型礫床水路で実験をた．また，粒子運動の流れ場に与える影響が大きく，平行い，石礫河道の土砂移動機構に関する多くの新しい結均流れと乱れ構造の把握が必要であること，さらに，礫果を得た．主要な結果は，洪水時に露出した大粒径粒子床のように透水性が高い河床上の流れは，主流部と堆積が抵抗となり，大粒径粒子が核となって周囲の中小粒子層内で水の出入りが生じ，これらが流れの抵抗特性 7)やの流出を抑制する機構を明らかにし，これを石礫河道の流砂運動に影響を及ぼす．

流砂機構に取り込み定式化したことである．関根ら は本研究では，異なる形状を持つ一様な粒径集団からな鉛直分級および河床の平衡状態の形成過程について一連る石礫群を用いた数値移動床実験を行い，流れ構造と河の水理実験により興味ある結果を得てきた．しかし，水床変動に及ぼす粒子形状の影響を明確化する．

2)

理実験で扱える粒径や形状に制限があり，得られる成果に限界を伴うことは否めない．福岡ら 3)は粒子の形状と2. 数値移動床実験 3)大きさを任意に取り込み，粒子群の運動と流れの相互作用を考慮した数値移動床水路を構築した．数値移動床水路を用いた大型水路実験は，石礫粒子の移動軌跡や運動

数値移動床実験では，福岡ら 3)の方法を用い，流れを等を適切に説明し，また，固定床土石流水路実験 4)の粒

Euler 的に，石礫粒子群は Lagrange 的に解析する．流れ子速度や水と土砂の輸送量を再現 していることから，の解析では固液混相流の一流体モデルを用いて石礫粒子その解析法の有効性が確認されている．数値移動床水路の固相部分を密度の異なる液体として与え，非圧縮性流5)1I_919

れとして解析する．LES により流れの計算を行い乱流モデルとして Smagorinsky モデルを用いている．以下に流れの連続式と運動方程式を示す．
ui=0xi

Dui
1 p= gi −+2 ( + t ) Sij Dt xi x j t = ( Cs  )2

2 S kl S kl

(1)
（2）
（3）

ここに，ui は流体計算格子内で質量平均した流速で，

図-1 粒子形状特性

添え字は各方向成分を示す．ρ は格子内で体積平均した密度，gi は重力加速度，p は圧力と SGS 応力の等方成分の和，νは動粘性係数，Sij はひずみ速度テンソル，Cs は(a) 扁平な形状

Smagorinsky定数である．自由水面の変化は VOF 法に基づき，次の流体計算格子内の液相の体積割合 f の連続式を用いて評価した．
f fui+=0t xi(b) 棒状図-2 異なる視点から見た粒子形状（4）

また，石礫粒子群の運動解析では，小球を隙間がない

t=14.8st=14.7s
離脱ように連結させた粒子のモデルを用いており，並進と回t=14.9st=15.0st=23.9st=23.8s
離脱t=24.0s

転の運動方程式を用いて剛体として解析する．個別要素法を用いて粒子の接触を考慮し，流れ場の解析から粒子

t=23.7s
停止

t=14.6s停止に作用する流体力を評価している．流体計算格子サイズ(a) 扁平な形状(Case1_x=10m)

Δ=0.01m である．解析法の詳細は文献 3)を参照されたい．

(b) 棒状(Case2_x=7.7m)

図-3 粒子の離脱表-1 実験条件

3. 数値移動床実験の結果

水路床からの鉛直高さz[m](1) 実験条件
流下方向，水路に平行な横方向および水路に垂直な上向きをそれぞれ x 軸，y 軸，z 軸にとり u，v，w を各方向の速度としている．数値実験水路は，水路長 15m，幅 1m，勾配 1/20 の一様な直線水路で，通水前に粒子をランダムに投下し初期河床を作る．実験に用いる粒子は多摩川支川の浅川で採取された石礫を参考に決めている．浅川
2.3km 右岸の砂州上で線格子法により石礫 50 個を採取し，それぞれの長径 a，中径 b，短径 c を計測した．結果を図-1 の緑色の点で示している．実線は Shape Factor(cCase1

Case2

Case3

Case4

0.70.60.5
0.4
0.3
0.20.1000.2

0.40.6
空隙率0.81

図-4 初期河床の空隙率

方向に向け，めくれ上がるように離脱をし，一方，棒状ab )を示す．浅川の石礫は Shape Factor の値が 0.4～の粒子は長軸を横断方向に向け，転がるように離脱する．

0.7 の範囲に多く見られる．図-2 は Shape Factor が 0.5～

この 2 種類の粒子は，同一の体積を有するようにつく

0.6 の範囲にある粒子の中で代表的であった扁平，棒状られており，換算粒径は同一体積の球の直径で d=70mmの粒子を小球の連続で形をモデル化したものを視点を変と定義した．実験条件を表-1 に示す．Case1～Case3 は粒えて示す．図-2(a)の扁平な粒子は不規則な凹凸を有し子形状の影響を明らかにすることを目的としている．投ており図-2(b)の棒状の粒子は滑らかな形状である．浅入した粒子群の総量はいずれも等しく，Case3 の粒子群は 2 種の粒子を同じ量混合したものである． Case4 は堆川の粒子群との関係を図-1 から見ることが出来る．図3(a)，3(b)に 2 種類の粒子の離脱の瞬間の連続写真を示

積層内の流れの影響範囲を明らかにする目的で河床の厚さを Case1～3 に比して 3 倍にしている．境界条件は，す．扁平な形状の粒子は，離脱の瞬間に平らの面を流下

2I_920水路床からの鉛直距離z[m]流下方向速度u[m/s]0.60.50.40.30.2後流域は見られない．
AB0.1
012.81313.413.213.6

水路 床

縦断距離x[m](a) 跳躍時(b) 離脱時
図-6 粒子運動と流れ場(Case1)

図-5 河床近傍の流れ場(Case1)

上流端で一定流量 0.5 m3 / s を与え，下流端の自由落下水脈内はほぼ大気圧となることから，下流端は圧力を 0 としている．給砂は下流端を通過した粒子を同時刻に上流x = 1～3 mの範囲でランダムに投下させている．図-4 は初期の堆積層内の空隙率の鉛直分布を示している．いずれの実験ケースも平らな水路底面に接して存在する礫層の空隙率は大きくなる．河床堆積層の大きい Case4 では，(a) 流速の鉛直分布(b) 対数分布則
図-7 流速分布(Case1_x=12m)

換算粒径より大きい水路床から 0.1 m 以上の礫の空隙率は約 0.36 と一定に近い．
(2) 主流と堆積層内の流れ
a) 礫床河川の流れ場
図-5 は Case1 の実験水路中心線上縦断面の流れの瞬間値を示している．式(1)～(3)で示す固液混相一流体モデ
ルの密度コンター図に，流下方向速度 u と鉛直方向速度w から成る 2 次元速度ベクトルを示す．図より主流部から堆積層内への吸い込みと湧き出しが確認できる．用いた石礫層に連続した空隙が存在することで，堆積層内で

図-8 乱れ強度分布(Case1_x=12m)

0. 2 m/s に近い流速が見られる．図-6 は実験水路中心線上縦断面における流下方向流速コンター図に，移動粒子と流れ場の関係を示している．図の白色の領域は固相を示している．本解析では各計算メッシュにおける密度を液相と固相の体積割合を用いて評価しているため，液相と固相の境界部では，密度の空間的な勾配が生じ密度コンター線の間隔が狭く多数描かれる．図-6(a)は跳躍中の粒子周りの流下方向速度 u を示しており，粒子跳躍中は粒子の移動速度と流れの速度差が小さいため粒子背後の後流域はあまり見られない．一方，図-6(b)に見られる停止していた粒子が離脱する瞬間には，粒子がめくれ

図-9 レイノルズ応力分布(Case1_ x=12m)

上がると粒子の背後に粒径ほどの高さの後流域が見られ，b) 主流部の流れ場その後遅い流体が輸送されていることが確認できる．河床表層上の粒子の運動は，転動や跳躍，回転運動等を伴Case1 で計算時間 t = 90～100s，x = 12mを対象に主流部った三次元的な運動のため，粒子運動の流れ場へ与えるの流れ場を分析した．各格子点の流速を 0.1 秒ごとに液影響は大きい．

相のみ抽出し気相および固相は母数に含めない．また，各格子点でサンプル数が 50 個以上得られた場合のみ平礫床河川の流れ場をさらに詳細に明らかにするため主流部および堆積層内について流速分布，レイノルズ応力， 均流速を求め，これよりレイノルズ応力および乱れ強度を算出した．また，各値の鉛直分布を調べる際は，水路乱れ強度について調べる．
3I_921

中心線上から横断方向に粒径 0.07m 分空間平均している．
図-7(a)は速度の 3 成分の鉛直分布を示している．t = 100s
時の x = 12 m水路中央地点の水位と河床高は，それぞれz = 0.45 m，0.19 mである．図-7(b)では，粗面の対数則と解析で求めた流下方向速度 u (z)を比較している．摩擦速度 U*は重力加速度 g，径深 R，エネルギー勾配 I を用いて(gRI)1/2 で算出し，値は U* =0.29 m/s であった．エネルギ
ー勾配 I は，x =5 m,11 m地点の平均水位に断面平均流速から速度水頭を求め全エネルギーの差をとり評価した．

図-10 堆積層内の流速分布(Case4_x=10m)

また，κ = 0.4，Ar = 8.5 として，ks は換算粒径 d と等しい
70 mm を用いた．流速分布の原点位置は河床高さから0.4d = 0.028 m上方に補正すると(b)の結果が得られ，対数則と概ね一致している．
図-8 は x =12m 水路中央地点の t =100s 時の河床高を基準に測られた乱れ強度の鉛直分布を示す．祢津 8)は，固定床の開水路乱流場の乱れ強度の鉛直分布を図-8 の実線で示している．図から z/h = 0.0～0.3 (z = 0.0～0.075 m)の河床面近傍では乱れ強度の三成分 u2 , v2 , w2 は，ほぼ

図-11 堆積層内の乱れ強度分布(Case4_x=10m)
水路床からの鉛直高さZ [m]

水面方向に直線的に発達し，固定床の乱れ構造と大きく異なる．また，z/h = 0.3 以上で u2 は減衰していくが，v2 , w2 はほぼ一定の値を示している．これらの原因として河床堆積層内の流れや河床の凹凸が考えられるが，乱れ構造に関して更なる検討が必要である．
図-9 は流下方向に作用するレイノルズ応力の鉛直分布を示す．z/h = 0.4 以下で等流のレイノルズ応力分布か平均水位（Case1_t=100s)
平均水位（Case2_t=100s）
平均水位（Case3_t=100s）

0.7
0.6
0.5
0.4
0.3
0.20.1

礫止めの堰

04らずれている．この高さは乱れ強度の鉛直分布が直線的

平均河床（Case1_t=100s)
平均河床（Case2_t=100s）
平均河床（Case3_t=100s）

5678

910 11
縦断距離x[m]12131415

図-12 水位と河床高の縦断形に発達していた領域とほぼ一致する．

t=50.0sc) 堆積層内の流れ場
次に堆積層を厚くした Case4 の堆積層内の流れ場につ核になる粒子いて調べる．図-10 は x = 10m の t = 25～35s の 10 秒間でx=12mの堆積層内での平均流速の鉛直分布を示す．v，w に着目すると正負，様々に分布しており，堆積層内の流れ場が複雑であることが分かる．図-11 に堆積層内の乱れ強1m度分布を示す．図中の①～④は図-10 の番号と対応して(a) 亀甲構造をなす粒子群 (b) 河床の平面状況(Case1_x=12m)
図-13 粒子の配置(Case1_t=50s)

いる．領域①では他の②～④に比べ乱れ強度が大きく主流の影響によるものと考えられる．主流部の流れが直接堆積層の流れに与える影響は z/h = -1 までの水深程度の範囲である．今回得られた結果は主流部の大きさ，堆積層の空隙率や粒径など河床材料特性によっても異なり，実験条件を変化させて堆積層内の流れについて，さらに検討が必要である．
(3) 河床変動a) 河床高と水位の時間変化
図-12 は Case1～3 の t = 100s 時の水路床からの縦断水位と河床高を示す．水位および河床高は，水路中央から横断方向に 1 mを空間平均して求めている．初期河床は平坦に作られている．通水により振幅が 0.1 m 以下の緩図-14 粒子衝突による礫の離脱(Case1_x=11.6m)
4

I_922

やかな波が河床面に形成され，水面と同位相となっている．x = 5～11m ，t = 100s 時での平均フルード数 Fr は，Case1～3 それぞれ Fr = 1.27 ,1.35 ,1.31 である．水路下流端
x = 15ｍには水路床から高さ 0.2 m の礫止め堰が設置されている．このため，3 ケースとも下流端付近では堆積傾向にある．
b) 粒子形状の効果と河床変動

図-15 河床高の時間変化(Case1，粒子形状:扁平)

河床表層粒子は流れの作用により姿勢や配置を変化させる．図-13(a)は実験で見られた河床表層面での粒子群が形作る特徴的な構造を示している．オレンジ色の粒子は，長軸を鉛直方向に向け礫床に長軸の半分程度埋没し(a) Case1(粒子形状:扁平)ており，容易には動かない．この粒子を核にして，周囲に粒子が堆積し 6 個程度の粒子集団を形成する．周囲に堆積する粒子は長軸を流下方向に向け，高さを増大させ，かつ平らな面を上に向けた覆瓦構造を取る．集団を構成(b) Case2(粒子形状:棒状)

する粒子群の粒子間の距離は近く，粒子同士，力を伝達し合っている．これは亀甲積みと呼ばれる安定な石の積se2

み方と類似の構造である．図-13(b)は Case1 の通水開始から 50 秒後の河床表面の状況を示す．白線で囲まれた(c) Case3(粒子形状:扁平,棒状の 2 種類)
図-16 河床コンター(Case1～Case3)

粒子は核となる粒子で，黒線で囲う粒子が粒子集団を示す．核になる石礫が多数存在するところでは，河床は安定している．粒子の形の効果は，粒子自身の姿勢だけでなく周囲の粒子と色々な配置を取りながら粒子群の移動に抵抗することが分かる．
次に河床の洗掘過程を説明する．図-14 は粒子衝突により河床粒子群が移動し，澪筋が形成されていく過程を示す．t = 76.2s の白色の破線は周囲より低い場所を示す．
着目している粒子に①～③の番号を付している．t =76.3s

図-17 通過礫量

では移動粒子①が流下し t = 76.5s では移動粒子①と停止t=48.5s

粒子③が衝突する．また，停止粒子②は移動粒子の通過

t=48.8sにより姿勢を崩し離脱している．このように河床の洗掘は，河床高が周囲よりも低くなるとその場所に移動して来る粒子との衝突により河床粒子の移動が起こりやすく，流下方向に澪筋が形成されるようになる．平坦に近い河床を粒子が流下する場合，堆積している粒子の上を流下(a) Case1(x=12m)(b) Case2(x=12m)
図-18 河床表層の礫の状況するため，河床は洗掘されにくい．図-15 は Case1 で洗を示す．流下方向速度が 0.05 m/s 以下の粒子を河床に堆c) 通過礫量
図-17 は Case1～3 について通水開始から 100 秒間に x =積している粒子と定義している．t = 75s および t = 100s か11m の断面を通過した礫量を示す．Case1 と Case2 の通ら流下方向に澪筋が発達していく様子が見られる．図-過礫量を比較すると，Case2 は Case1 の約 1.85 倍となり16 は Case1～Case3 の t = 100s の河床コンター図を示す．
Case1 では，x = 6～8m 地点の洗掘箇所でも黒の鎖線に示

棒状の粒子は移動しやすい形状といえる．一方，Case3（2 種類）の通過礫量は，Case1 と Case2 の通過礫量のほすように澪筋の発達が見られる．Case2 の棒状の粒子群

ぼ中間の量であるが，扁平な粒子と棒状の粒子の通過礫では，澪筋は横断方向にも広がったため両岸の水路壁まで一様に洗掘しており，より動き易い河床であることが量の内訳見ると，両者にあまり差が現れていない．これは，それぞれ単独な形状の粒子に対して Case3 では，扁分かる．Case3 では図-16(c)の黒の鎖線に示すように水路平な粒子と棒状の粒子が互いの運動に影響を及ぼし合っ中央が一様に洗掘している．このように粒子形状によっているためである．図-18 で河床表層の礫の状況を比較て粒子群の動き易さや河床の変動に違いが見られる．

する．図-18(a)，(b)の白色の鎖線は図-13(a)で示した粒掘が見られた x = 11～13ｍの河床高コンターの時間変化

5I_923

子集団の核になる粒子を示している．Case1 の方が核に本検討では，グリッドスケール以下の渦を，スマゴリなる粒子が多く，より安定した河床となっている．扁平

ンスキーモデルでモデル化することで河床の凹凸や移動な粒子は平らな面が大きいため起立し，核になる粒子に粒子による乱れ構造を概ね解析できているとの立場で分なりやすい．一方，棒状の粒子が核になるためには，周析を行った．今後石礫河道の粒子運動と乱れの構造につ囲の複数の方向から支持される必要があり，核粒子にはいてさらに調査・検討する予定である．

なりづらい．さらに，扁平な粒子は不規則な凹凸を有しており，かみ合わせ効果が働くため周囲の粒子に力を伝参考文献

え，より離脱しにくい状況にある．これらのことから通1) 福岡捷二: 石礫河川の移動床水理の諸問題と解決への道筋，水工学に関する夏期研修会，Aコース，Vol.44，pp.A-1-A25，
2008.
2) 関根正人，平松裕基，中川裕貴：骨格が大礫により構成される河床上で生じる鉛直分級のメカニズム，土木学会論文
集B1(水工学)，Vol.74， No4， pp.I_1069-I_1074，2018.
3) Fukuoka, S., Fukuda, T.,Uchida, T. : Effects of sizes andshapes of gravel particles on sediment transports and bed variations in a numerical movable-bed channel，Advances in
Water Resources, Vol.72, pp.84-96, 2014.
4) 伊藤隆郭，江頭進治，宮本邦明，竹内宏隆：土石流の固定床から移動床流れへの遷移過程に関する研究，水工学論文集，第43巻，pp.635- 640，1999.
5) Fukuda,T., Fukuoka,S. ： Interface-resolved large eddy
simulations of hyperconcentrated flows using spheres
and gravel particles， Advances in Water Resources ，2017.
https://doi.org/10.1016/j.advwatres.2017.10.037
6) 田所 弾，福田朝生, 福岡捷二：等価な粒度分布を有する球と石礫で構成される二つの数値移動床上の粒子運動の比較検討，土木学会論文集B1(水工学)，Vol.70，No.4，pp.I_961I_966，2014.
7) 中川博次，祢津家久：透水性開水路乱流場における主流と浸透流との相互作用，土木学会論文報告集，No.244，pp.8190，1975.
8) 祢津家久：開水路乱流の乱れ強度に関する研究，土木学会
論文報告集，No.261，pp.67-76，1977.

過礫量に大きな差異が生まれたと考えられる．

4. 結論
本研究では，粒子形状と粒子の混合割合が異なる数値移動床実験を 3 ケース行い，河床形状や石礫移動量に及ぼす形の効果を検討した．また，主流部と堆積層内の流れについて検討した．以下に主要な結論を示す．
(1) 本研究条件では，堆積層内の流れの流速が 0.2m/s であった．この値は，砂などの小粒径の運動に影響を与える量となる．また，乱れ構造については固定床で検討された祢津の結果と異なっている．さらに，主流部と堆積層内の流れの相互作用の及ぶ範囲は，河床表層から鉛直下向きに水深程度である．
(2) 粒子形状の違いによって河床変動の発達速度に差が見られた．複雑な形状（Case1 の実験）の粒子の方が粒子群が塊となって，安定な河床構造を取る．また，河床の洗堀は，移動粒子が河床の静止粒子と衝突・移動をもたらし澪筋を形成していくことが主因と考えられる．
(3) 粒子形状は自身の運動だけでなく周囲の粒子の運(2018.5.31 受付)動にも影響を与える．粒子形状の違いによって河床表層の堆積構造が異なり，輸送される石礫量が大きく異なる．

EFFECTS OF PARTICLE SHAPES OF GRAVELS ON THE FLOW STRUCTURESAND THE BED GEOMETRY IN GRAVEL-BED STREAMS
Takatoshi ATSUMI and Shoji FUKUOKAIt has been recognized that particle shapes of gravels have essential roles for hydraulics in the gravelbed rivers. The numerical movable-bed channel experiments are conducted using two different kinds ofgravels with uniform particle shapes. Flow structures as well as mutual flow interactions between themain flow region and permeable gravel bed region are investigated by the numerical movable-bed channel experiment.
It is demonstrated by the numerical channel experiment that the lateral turbulent mixing near the gravel-bed surface is intense due to chages in the bed geometry and the main flow affects over the range offlow depth on the structure of flow in the permeable gravel-bed region and also the developing process of
bed configulations is different markedly by the particle shapes of gravels and how big blocks of gravels
and subsequent sand waves are formed on the bed.
6I_924