河川復元・バンクフル流量設計シミュレーター 戻る
水文・自然河川工学

河川復元・バンクフル流量設計シミュレーター

人工改修された河川を自然形態へ戻す復元設計のツールです。流域面積・川幅・水深・河床勾配・粒径・マニング係数を変えると、バンクフル流量、流速、幅深さ比、Rosgen河型、河床粒子の移動限界がリアルタイムで分かります。

パラメータ設定
流域面積 A
km²
川幅 W
m
川深さ D
m
河床勾配 S
m/m
中央粒径 d50
mm
河床材料の代表粒径
マニング n
自然河川で 0.030〜0.060 が標準
Rosgen 河型
復元の目標形態
計算結果
流積 A_xs (m²)
平均流速 V (m/s)
バンクフル流量 Q_bf (m³/s)
地域曲線推定 Q_reg (m³/s)
幅深さ比 W/D
推奨 Rosgen 型
河川断面・氾濫原 — Rosgen 型可視化

中央が河道、両側が岸と氾濫原。青の濃さが流速、黄色の点が粒径 d50、上部のラベルが推奨 Rosgen 型。

地域回帰曲線 — Q_bf vs 流域面積
Rosgen 型ごとの W/D 比の目安
理論・主要公式

$$Q_{bf} = \frac{1}{n}\,A_{xs}\,R^{2/3}\,S^{1/2},\qquad \tau_b = \rho\,g\,R\,S$$

マニング式によるバンクフル流量 Q_bf と境界せん断応力 τ_b。n:マニング粗度係数、A_xs:流積、R:水力半径(=A_xs/P、P=W+2D:潤辺)、S:河床勾配、ρ=1000 kg/m³、g=9.81 m/s²。

$$Q_{reg} = 2.5\,A^{0.85},\qquad \tau_c = 0.045\,(\rho_s-\rho)\,g\,d_{50}$$

地域回帰曲線(Q_reg:m³/s、A:km²)と Shields の臨界せん断応力 τ_c。粒径 d50 が動くのは τ_b > τ_c のとき。本ツールでは石英質粒子の水中比重から ρs−ρ = 1650 kg/m³ を採用。

河川復元・バンクフル流量設計 — Rosgen 分類

🙋
「河川復元」って、何を「復元」するんですか?コンクリート護岸を剥がしてもとに戻すこと?
🎓
それも含まれるけど、もう少し広い概念だね。20世紀の治水で多くの河川は直線化・三面張りにされ、蛇行も氾濫原も失った。河川復元(Stream Restoration)は「自然の河川がもっていた形(蛇行・瀬と淵・氾濫原)と機能(生態系・地下水涵養・洪水緩衝)」を取り戻す環境工学のことだ。米国ではNRCS や EPA、欧州ではWFD(水枠組み指令)の枠組みで進められているよ。
🙋
なるほど。デフォルトの設定で見ると「Q_bf=35.9 m³/s」と「Q_reg=69.5 m³/s」が出ています。この2つの違いは何ですか?
🎓
いいところに気がついたね。Q_bf はいま入力した「川幅 12m・水深 1.5m・勾配 0.5%・マニング 0.04」の断面でマニング式から計算した流量。一方 Q_reg は「流域面積 50km² の似たような地域なら、過去の観測ではバンクフル流量はだいたい 69.5 m³/s だよ」という統計値だ。両者が大きく食い違うと、いま設計しようとしている断面が「自然な大きさ」から外れていることを示すサイン。デフォルトでは Q_bf が Q_reg のおおむね半分、つまり川を細く絞り込みすぎている可能性がある。
🙋
推奨 Rosgen 型に「E」と出ていて、選んだ「C」と違うのも気になります。これは何を意味しますか?
🎓
Rosgen 分類は「幅深さ比 W/D」と「勾配」から河川の自然な形を決める枠組みだ。デフォルトは W/D=8(12÷1.5)で勾配 0.5% だから、判定式では「W/D<12 かつ 勾配<2%」→ E 型(深く狭い小川)が自然形になる。一方ユーザーは C 型(蛇行幅広、W/D≥12)を目標にしている。これは「現在の断面は E 型寄りなのに C 型に作りたい」という不整合で、復元設計でよくある落とし穴なんだ。C 型にしたいなら、川幅をもっと広げて W/D≥12 にする必要がある。
🙋
河床粒子の移動判定(Shields数)は復元設計でどう使うんですか?
🎓
復元では「自然の力で河床が自分自身を整える」のが理想なので、バンクフル流量で粒子がギリギリ動くくらいに粒径を選ぶのがセオリーだ。動かなさすぎると河床が固まって生態系が単調になり、動きすぎると護床材が流出してしまう。Shields数 τ* が 0.045 を超えるか否かが目安で、本ツールの「is mobile」フラグはまさにそこを見ている。例えば渓流の復元で d50=200mm の玉石を敷くと τ_c が大きくなりすぎて全く動かない、というのは典型的な失敗パターンだよ。
🙋
最後にひとつ。なぜ「バンクフル」が再現期間 1.5〜2 年なんですか?もっと大きな洪水のほうが河川を変えそうなのに。
🎓
直感に反するけど、河床形態を最も支配するのは「中程度の流量×それなりの頻度」の組み合わせなんだ。100年洪水は強烈だけど100年に1回しか来ない。雨が降るたびの小流量は頻繁でも仕事率が小さい。両者の積(仕事率×頻度)が最大になるのが、ちょうど1.5〜2年再現期間の流量で、これが Wolman & Miller (1960) 以降「形成流量=バンクフル流量」と呼ばれる根拠だ。だから護岸・粒径・氾濫原幅は全部この流量で決める。

よくある質問

バンクフル流量(Bankfull Discharge)は、河川の水面が両岸の上端、つまり氾濫原(floodplain)に溢れ出す直前まで達したときの流量です。再現期間でいうとおおむね1.5〜2年に一度発生する流量に相当し、河床形態を最も支配する「形成流量(dominant discharge)」として扱われます。本ツールではマニング式 Q = (1/n)·A·R^(2/3)·S^(1/2) で計算します。
Rosgen分類(Rosgen 1994/1996)は河川の形態を勾配・幅深さ比・蛇行度・河床材料からA〜Gの7タイプ+細分類に分けます。代表的にはB型(半閉鎖・中勾配)、C型(蛇行・氾濫原が広く幅深さ比≥12)、D型(編流・砂礫床)、E型(深く幅深さ比<12の小川)があります。復元設計では「いま強制改修されている河川」と「自然に近い形態」を比較し、目標形を決めるために使います。
地域回帰曲線(Regional Curve)は、観測流域でバンクフル流量と流域面積の関係を統計的に回帰した式で、一般に Q_bf = a·A^b(aは地域定数、b≈0.7〜0.9)の形をとります。本ツールでは便宜上 Q_bf = 2.5·A^0.85(A:km²、Q:m³/s)の標準的な値を使用します。実務では地域ごとに係数が大きく異なるため、対象流域の観測曲線を必ず確認してください。
Shields数(無次元せん断応力)τ* = τ/((ρs−ρ)·g·d) が臨界値(砂礫で約0.045)を超えると、河床粒子は流れにより動き出します。本ツールでは τ = ρgRS(境界せん断応力)と τc = 0.045·(ρs−ρ)·g·d50 を比較し、バンクフル流量で中央粒径が移動するか否かを判定します。「動かない=硬い河床」「動く=活発な河床」となり、復元工事での粒径選定の根拠になります。

実世界での応用

都市河川の自然再生:東京の野川、横浜のいたち川、大阪の長居公園調節池など、コンクリート三面張りから自然形に戻す事業で本ツールのような評価が前提になります。流域面積から地域回帰曲線でバンクフル流量を推定し、Rosgen型に合わせて目標断面を決定。マニング n は植生・河床粒径から 0.035〜0.055 を使い、計算流量と推定流量の整合を確認します。

米国 NRCS / EPA の復元プロジェクト:米国農務省自然資源保全局(NRCS)とEPAは1990年代から Rosgen 分類をベースにした河川復元(Natural Channel Design, NCD)を標準化し、ノースカロライナ州など東部の何百もの河川で適用されています。設計の出発点は「参照流域(reference reach)」の幅深さ比・蛇行波長・粒径分布を測り、それを目標流域に転写する考え方です。

砂防・治山事業との接続:渓流復元では本ツールのような計算で「自然の輸送能力を超えない護床材径」を選びます。τ_c が τ_b の2倍以上だと護床は動かず固結し、生態系も貧弱になります。逆に τ_c < τ_b なら一回の出水で粒径が抜けてしまう。両者がほぼ等しくなる粒径を選ぶのがプロのコツです。

魚道・河川改修の事前検討:サケ・マス類の遡上を考慮する河川では、バンクフル流量時に水深が一定以下(0.3m以下)となる「瀬」が必要で、Rosgen C/E 型の特徴的な瀬-淵構造を再現することが目標になります。本ツールで W/D 比と勾配を変えながら、瀬の流速・水深を試算するのが典型的な使い方です。

よくある誤解と注意点

まず最大の落とし穴が、「バンクフル流量=計画洪水流量」と混同することです。バンクフル流量は再現期間 1.5〜2 年の「形成流量」であり、治水計画でいう 50〜100 年洪水とはまったく別物です。Rosgen 型の判定や粒径選定はバンクフル基準で行い、護岸高さ・堤防天端高は計画洪水で別途決めます。両者を混同して「バンクフルで設計した断面で 100 年洪水を流そうとする」と、堤防越流の重大事故につながります。

次に、「Rosgen 型を機械的に当てはめると不自然な河川になる」こと。Rosgen 分類は北米の温帯渓流の観測から構築された経験的な体系で、日本のように急峻で粒径分布が幅広い河川にそのまま適用すると、特に D 型(編流)と F 型(深溝)の境界などで判定がぶれます。本ツールの判定はあくまで目安で、実務では参照流域の現地調査(reference reach survey)を必ず行い、ローカルなキャリブレーションを掛けてください。

最後に、「地域回帰曲線をどの流域にも使える」と誤解すること。本ツールが使う Q_bf = 2.5·A^0.85 は標準的な係数ですが、地域定数 a は降水量・流出率・地形で大きく変わります。例えば乾燥地(米国南西部)では a≈0.5、湿潤多雨地域(米国太平洋岸北西部)では a≈8 と1桁以上違うことも珍しくありません。日本では国土交通省の流出特性解析や、近隣の観測点での実測値で必ず補正してから使用してください。

使い方ガイド

  1. 流域面積(aNum)を設定します。例えば関東平野の中規模河川で50~100km²の範囲を入力し、aRangeスライダーで変動幅を±20km²に調整
  2. 川幅(wNum)と平均水深(dNum)を入力します。利根川支流の復元設計例では川幅25m、水深1.8mが標準値
  3. 河床勾配(sNum)をパーセント単位で設定し、sRangeで勾配変化シナリオを検討。山地河川0.5~1.5%、平野部0.01~0.05%が一般的
  4. 計算実行後、バンクフル流量Q_bf(m³/s)と幅深さ比W/Dからデザイン河型(Rosgen分類)が自動決定される

具体的な計算例

多摩川支流の河川復元プロジェクト:流域面積80km²、川幅22m、平均水深1.6m、河床勾配0.08%の条件下。マニング粗度係数n=0.035(玉石混合河床)を適用すると、平均流速V=0.92m/s、流積A_xs=35.2m²が算出され、バンクフル流量Q_bf≒32.4m³/sとなります。地域曲線推定値Q_reg≒35.8m³/s、幅深さ比W/D≒13.75から推奨Rosgen型はC型(蛇行傾向、中勾配)に分類されます

実務での注意点

  1. 出水期(6~9月)の観測流量データと照合し、シミュレーション値の±15%範囲内に納まることを確認。河床材料変化(砂礫→玉石)でnが0.03~0.05に変動するため複数シナリオ計算
  2. 護岸復元設計では計算流速に安全係数1.3を乗じた値で侵食検討を実施。例えば護岸ブロックの許容流速3.0m/sに対し、現地V=0.92m/sであれば余裕が確保
  3. 気候変動対応で流域面積や降雨強度が増加する場合、aNum+20%、勾配シナリオを併行検討し、河積拡大が必要か判断