====== 水面形の計算 ======

<fc #ff00ff>**TYS-RBA**</fc>モデルでは、任意断面形、台形、複断面形の３種類の断面形状を対象としている。

各断面の流れの計算は、任意断面形と台形断面についてはレベル１の平均流公式を用い、複断面形についてはレベル２の平均流公式を用いている。

ここでは複断面形の水面形の計算について詳述する。

===== 断面形状 =====

本モデルは、次の任意・台形・複断面の３種類の断面形状を対象としている。

==== 任意断面形 ====

具体的な断面形状を特定しないが、水深の変化による流積・径深の関係を次の近似式で表現する。

流積：$A=C_1*h^{C_2}$(h：水深)

径深：$R=C_3*h+C_4$>

係数$C_1$～$C_4$の４つのパラメータを事前に最小自乗法で求めておき、それをプログラム実行時に読み込む。

=== 任意断面形の流積と径深の求め方 ===
==== 台形断面 ====

下図のように河床幅(B)と左右岸勾配(B1,B2)の３つのパラメータを読み込む。

;#;
{{:tysrba:image05.jpg|}}

図-3.1　台形断面
;#;


==== 複断面形 ====

河床断面をできるだけ簡略した形にモデル化し、下のような複数の矩形の組み合わせで複断面形を表現する。

;#;

{{ :tysrba:image06.jpg?nolink& |}}

図-3.2　複断面
;#;

 
===== 水理諸元 =====

水深の変化に応じて、河道断面の水理諸元は変化する。河道の断面形状に応じた流積・径深を以下にまとめる。

==== 任意断面形 ====

任意断面形の流積・径深は、水深：$h$とすると、以下のように算出される。

流積：$A=C_1*h^{C_2}$　・・・(3.1)

径深：$R*C_3*h+C_4$　・・・(3.2)

==== 台形断面 ====

台形断面の流積・径深を河床幅：$B$、左右岸勾配：$B1$・$B2$を用いて表すと以下の通りである。

流積：$\large A=B h + \frac{(B1+B2)}{2}h^2$ ・・・(3.3)

径深：$\large R=\frac {B h + \frac{(B1+B2)}{2}h^2} {B+h( \sqrt{1+B1^2} + \sqrt{1+b2^2})}$　・・・(3.4)

==== 複断面 ====

任意・台形断面は水深の変化に対する流積・径深の変化は線形であるが、複断面の場合は、低水敷～高水敷の境界で断面の変化点があるため、流積・径深の変化は線形ではない。

図-3.2のような複断面において、河床変動による河床高に変化に伴って各分割断面の流積$(A)$・潤辺$(S)$・径深$(R)$の算出式は異なる。河床位置に応じた各分割断面の水理諸元を整理する。　ここで、現況の河床高：$Z$、初期の河床高：$Z0$、水深：$h$、高水敷の高さは$Z1>Z2$とする。

=== $Z<Z0+Z2$の時(河床高が低水敷の範囲) ===

== ①$Z+h<Z0+Z2$ ==

$A0=h \cdot B$

$S0＝B+2h$

== ②$Z0+Z2≦Z+h<Z0+Z1$ ==

$A0=h \cdot B$

$S0=B+h＋(Z2+Z0-Z)$

$A2=(h-Z2-Z0+Z) \cdot B2$

$S2=B2+(h-Z2-Z0-Z)$

== ③$Z0+z1≦Z+h$ ==

$A0=h \cdot B$

$S0=B+Z1+Z2+2Z0-2Z$

$A1=(h-Z1-Z0+Z) \cdot B1$

$S1=B1+(h-Z1-Z0+Z)$

$A2=(h-Z2-Z0+Z) \cdot B2$

$S2=B2+(h-Z2-Z0+Z)$

=== $Z0+Z2≦Z<Z0+Z1$の時(河床高が低い方の高水敷の範囲) ===

== ①$Z0+Z2≦Z+h<Z0+Z1$ ==

$A0=h \cdot B$

$S0=B+h$

$A2=h \cdot B2$

$S2=B2+h$

== ②$Z0+Z1≦Z+h$ ==

$A0=h \cdot B$

$S0=B＋Z1+Z0-Z$

$A1=(h-Z1-Z0+Z) \cdot B1$

$S1=B1+(h-Z1-Z0+Z)$

$A2=(h-Z2-Z0+Z) \cdot B2$

$S2=B2+(h-Z2-Z0+Z)$

=== $Z0+Z1≦Z$の時(河床高が高水敷よりも上) ===

== ①$Z0+Z1≦Z+h$ ==

$A0=h \cdot B$

$S0=B$

$A1=h \cdot B1$

$S1=B1+h$

$A2=h \cdot B2$

$S2=B2+h$
 
===== 等流水深 =====

==== 任意断面・台形断面 ====

任意断面・台形断面の等流水深は、次式のマニングの平均流公式が成立する水深を示す。

$\large u=\frac {Q}{A} = \frac {1}{n}R^{\frac{2}{3}}I_b^{\frac{1}{2}}$　・・・(3.5)

ここで、$Q$流量、$A$：流積、$R$：径深、$n$：粗度係数、$I_b$：河床勾配である。

==== 複断面 ====

複断面を高水敷・低水敷で複数の断面に分割して考えた場合、任意の分割断面iの平均流速公式は、

$\large u_i=\frac {Q_i}{A_i} = \frac {1}{n_i}R_i^{\frac{2}{3}}I_b^{\frac{1}{2}}$　・・・(3.6)

ここで等流水深は(3.6)式が成立する水深を示す。

$Q_i$：分割断面の流量、$A_i$：分割断面の流積、$R_i$：分割断面の径深、$n_i$：分割断面の粗度係数である。

河道断面ごとに断面の合計流量$Q$と河床の縦断勾配$I_b$が既知であるので、水深：$h$を仮定して各分割断面の流積$(A_i)$・径深$(R_i)$を求め、(3.6)式より各分割断面の流量$Q_i$を計算する。そして、合計の流量が$Q$に等しくなるまで（実際の計算では所定の許容誤差内収まるまで）水深を変えて収束計算を行う。

また、(3.6)式において河床勾配$<0$になる状況では当該断面の等流計算が行えないため、本プログラムでは、河床勾配の算出方法として

  -入力データで区間を定め、区間内の平均河床勾配＝当該区間の断面の河床勾配とする。
  -各断面の河床勾配をそのまま求め、河床勾配が負になるときは、等流水深$=$とする。

の２つの処理を行えるようにした。

1.の方法は、計算水深$<$限界水深もしくは、水深が求められない($=$収束計算が収束しない)場合、水深$=$限界水深に置き換える方法(一般的な不等流計算の常流近似)に適応される。

2.の方法は水深計算の等流水深への置き換え処理(詳細は3.6参照)に用いるものとする。これらの切り替えは入力データではなく、プログラム内に記述したフラグの選択で行うようになっているが、通常は①の処理を選択する。

===== 限界水深 =====

==== 任意断面・台形断面 ====

限界水深は、次式が成立する水深である。

$\large \frac{Q^2}{gA^3} \frac{\partial A}{\partial h} =1$　・・・(3.7)

任意断面・台形断面の限界水深は、3-2の水理諸元を用いて求める。

==== 複断面 ====

建設省河川砂防技術基準（案）同解説「調査編」P119の方法より、複断面の限界水深は以下の式が成立する水深を示す。

$\large \frac {U}{\sqrt {gR_c}} =1$　・・・(3.8)

ここで、$U$：全断面の平均流速、$R_c$：井田による合成径深で次式により求められる。

$\large R_i= \{ {\frac{\sum R_i^{\frac{2}{3}} A_i} {A}} \} ^{\frac {3}{2}}$　・・・(3.9)

ここで、合計流量$Q$が既知であるので、水深：$h$を仮定して、各分割断面の流積を計算し、(3.8)式が成立するまで（実際の計算では所定の許容誤差内収まるまで）水深を変えて収束計算を行う。
 
===== 不等流計算 =====

本プログラムは常流近似の不等流計算によって水面形を決定する。すなわち下流側の水深を既知として上流側の水深を収束計算によって求める。

不等流の運動方程式は以下の式で表される。（建設省河川砂防技術基準（案）同解説「調査編」P115より）

$\large \frac {1} {gA} \frac{\partial}{\partial x}(\int u^2 dA) + \frac {\partial H}{\partial x} + \frac{T_r}{\rho g A} = 0$　・・・(3.10)


==== 単断面の場合（台形断面） ====

単断面（レベル1の平均流速公式）の場合、運動量補正係数：$\alpha = \frac{1}{A} \int \frac{u^2}{U^2}dA$ （一定）とすると、

$\large \frac {\alpha}{2g} \frac{\partial}{\partial x}(\frac{Q}{A})^2 + \frac{\partial H}{\partial x} + \frac{T_r}{\rho g A}$ =0　・・・(3.11)

ここで、摩擦勾配：$\frac{T_r}{\rho g A}= \frac{n^2Q^2}{A^2R^{4/3}}$とおき、(3.11)式を標準逐次計算法による表記で表すと以下のようになる。

$\large \{H_U + \frac {\alpha}{2g} (\frac {Q_U}{A_U})^2 \} - \{H_D + \frac {\alpha}{2g} (\frac {Q_D}{A_D})^2 \} = \frac{1}{2}\{ \frac{n_D^2 Q_D^2}{A_D^2 R_D^{4/3}} + \frac{n_U^2 Q_U^2}{A_U^2 R_U^{4/3}} \} $　・・・(3.12)

ここで、$h$：水位、$\Delta x$：区間距離を示し、変数の添え字$U$・$D$はそれぞれ上流側・下流側を示す。(3.12)式が通常の矩形・台形断面で用いられている運動方程式の形式である。

下流側の断面の水理諸元と上流側の断面流量が既知として、上流側の水深を仮定して、(3.12)式が成立するまで（実際の計算では所定の許容誤差内収まるまで）水深を変えて収束計算を行う。

 
==== 複断面の場合 ====

複断面の場合、レベル２の平均流速の公式より、

$\large \frac{T_r}{\rho g A} = \frac{A^2(Q/A)^2}{(\sum \frac{1}{n_i}R_i^{2/3}A_i)^2}$　・・・(3.13)

$\large \int u^2 dA = Q^2 \frac{\alpha_i \sum \frac{1}{n_i}R_i^{4/3}A_i }{(\sum \frac{1}{n_i}R_i^{2/3}A_i)^2}$　・・・(3.14)


(3.13)・(3.14)式を(3.10)式に代入する。

これを標準逐次計算法で表すと以下のようになる。

$\large \{H_U + \frac {1}{gA_U} (\int u^2 dA)_U \} - \{H_D + \frac {1}{gA_D} (\int u^2 dA)_D \} = \frac{\Delta x}{2}\{ \frac{T_{rD}}{\rho g A_D} + \frac{T_{rU}}{\rho g A_U} \} $　・・・(3.15)

以上より、(3.12)・(3.15)式を組み合わせれば、単断面・複断面が混在する河道についても不等流計算が実施できる。

尚、実際のプログラムは(3.12) ・(3.15)の右辺における1/2の部分が岡部らによる区間内平均エネルギー勾配の近似法によって、上流側と下流側でそれぞれ異なる係数で計算を行っている。

区間内平均エネルギー勾配の近似法：\\
「急流河川の１次元河床変動(その１)」：岡部　健士、砂防学会誌Vol.50,No.3,pp58～65,1997
 

==== 複断面と台形断面が混在する部分の不等流計算 ====

不等流計算は、下流側断面の水深を既知として、上流側の水深を計算していくが、同じ河道上で複断面と台形断面が混在する場合、平均流速公式のレベル１と２が混在することになる。

平均流公式によって、運動方程式の各項の取り扱いは異なるため、同じ項でも計算値のオーダーが全く異なる場合もあり、レベルの違う両者をそのまま１つの運動方程式で混在して取り扱うことは水理学的にも無理がある。

よって、以下のように平均流公式のレベルをそろえた形の取り扱いを行う。

=== 上下共に台形断面 ===

レベル１の平均流公式を用いて計算を行う。

=== 上下共に複断面 ===

レベル２の平均流公式を用いて計算を行う。

=== 下流が台形断面・上流が複断面 ===

レベル２の平均流公式を用いて計算を行う。下流の台形断面の箇所についても(6)・(7)式を適用する。

=== 下流が複断面・上流が台形断面 ===

レベル１の平均流公式を用いて計算を行う。下流側の複断面の流積：$A= \sum A_i$、径深：$R=R_c$（井田の合成径深）として取り扱うものとする。

===== 射流区間の取り扱い =====

本プログラムは、常流近似の不等流計算によって水面形の計算を行っている。当該断面の水深を計算するためには、下流側の水深(既知)を用いて運動方程式が許容の誤差範囲内で成立するように収束計算を行う。この時、求められた水深の値が限界水深より小さいもしくは収束しない場合、以下の２種類の処理によって当該断面の水深を決定する。

==== 一般的な常流近似の処理 ====

　当該水深を限界水深に置き換える。それによって、各断面の水深は常に限界水深以上の水深となる。

==== 等流水深への置き換え ====

等流水深は、河床勾配によって大きく変化する可能性がある。特に河床勾配が小さくなった断面では、前後の断面と比べて等流水深が不自然に大きな値をとることとなり、射流区間でそのまま等流水深に置き換えると、下流断面水位と比較して不連続な水面形になってしまう。よって、以下のような処理を行う。

=== 当該断面で水深計算の結果 ===

  *解が収束しなかった時、もしくは解が収束するが収束解＜限界水深の場合
    *a.等流水深が値を持つ時：限界水深と等流水深の小さい方の値を用いる(射流域フラグを立てる)
    *b.等流水深が値を持たない時：「下流断面の水位と同じになるような水深」と「限界水深」を比べ、大きな用の値を用いる(射流域フラグ立てない)

  *解が収束し、かつ収束解≧限界水深の場合(射流域フラグ立てない)
    *c.そのままの解を用いる

=== 上流断面を計算する場合 ===

  *（下流側に射流フラグが立っているとき）→下流断面の水深を限界水深に置換。
  *（下流側に射流フラグが立っていないとき）→下流側の水深はそのまま。

=== 上流側の水深計算の結果====
  *上流の水深の扱いはa～c
  *下流水深の扱い（下流側に射流フラグが立っているとき）
    *上流がaの処理：水深を限界水深に置き換える前の値に戻す
    *上流がbの処理：水深＝限界水深のまま
    *上流がcの処理：水深＝限界水深のまま

　これら２種類の方法は、入力データではなくプログラム内のフラグ（変数名：iDep）の定義で操作の切り替えを行う。デフォルトでは①の方法を選択するようになっており、通常は①の方法を用いる。

===== ダム地点の取り扱い =====

通常のダム地点の取り扱いは、河道断面データ：KGS2とダム諸元のパラメータの指定によって行う。ダム・堰堤のみではなく、本川最下流端の水位を一定値で与えたい場合や支川の最下流端水深を限界水深で与えたい場合等、不等流計算以外の方法で水深を設定する場合に用いる。


==== ダムのタイプ ====

ダムのタイプには、以下の４種類がある。
    *0:水位の指定のみ
    *1:不透過ダム
    *2:スリットダム
    *3:限界水深
    *4:限界水深と等流水深の比較

　ここで、「限界水深と等流水深の比較」とは、水深を限界水深とし、上流断面水深＝等流水深の時、限界水深＞等流水深であれば、等流水深を用いる。

==== 起算水位or水通し標高（水位の指定のみ・不透過ダム・スリットダムの場合） ====

 「水位指定のみ」選択時には、限界水深を下回らない範囲で、「断面の水深＝起算水位－河床高」で当該断面の水深を設定する。

　「不透過ダム・スリットダム」選択時には、「設定値＝ダムの水通し標高」とし、他のパラメータ（水通し幅・流量係数等）と共に堰堤の式より、流量から水深を逆算し、「その値＋水通し標高－断面標高」を当該断面の水深とする。

==== 水通し幅（不透過ダム・スリットダムの場合のみ）====
　「不透過ダム・スリットダム」水通し幅を設定する。

○	水通し部側岸勾配（不透過ダム・スリットダムの場合のみ）
水通し部分の形状が台形の場合の側岸勾配を設定する。勾配の設定は河道断面が台形形状の場合と同じである。

○	スリット部高さ（スリットダムの場合のみ）
「スリットダム」の場合、水通し部分とスリットの下端の高さを設定する。

○	スリット部幅（スリットダムの場合のみ）
「スリットダム」の場合、スリット部の幅を設定する。

○	流量係数（不透過ダム・スリットダムの場合のみ）
　「不透過ダム・スリットダム」の堰堤の式に用いる流量係数を設定する。

==== 流砂量の取り扱い（不透過ダム）====

「不透過ダム」選択時には、以下のように掃流砂を調節する。
a.河床高＜水通し標高　：　掃流砂量をカット
b.河床高≧水通し標高　：　掃流砂量を流す
すなわち、ダムが満砂状態になるまでは、浮遊砂・ウォッシュロードのみを流下させる。

==== スリットダムの取り扱い ====


「スリットダム」選択時には、上流断面との区間距離の長短で水深・流砂量の取り扱いの方法が異なる。

/*「スリットダム」選択時には、上流断面との区間距離の長短で水深・流砂量の取り扱いの方法が異なる。*/

/*=== 区間距離（断面ピッチ）の長短によって場合分け ===

【図-3.3参照】

区間距離l1（砂防えん堤断面＝断面0～砂防えん堤直上流断面＝断面1までの距離）が、断面0から堆砂肩までの距離l0のn倍（当面1.5倍）より長いか短いかで場合分け
l1＞n・l0・・・区間距離が長く堆砂肩を表現できない場合　⇒（2）簡便法へ
l1≦n・l0・・・区間距離が短く堆砂肩を表現できる場合　　⇒（3）まじめ法へ

※断面0から堆砂肩までの距離l0
堆砂肩がダム高と同じ高さまで堆積したと仮定した場合（＝ダム高以上には堆積しないと仮定）における砂防えん堤地点から堆砂肩までの距離（想定される最大の長さ）。この値をもとに、"簡便法"あるいは"まじめ法"の適用を判断する。
具体的には、ダム高（初期河床高から水通しまでの高さ（スリット下端を初期河床にすりつけた場合はスリット高））hdmをtanφ（φ：水中安息角30～35°　当面35°？）で除した値　　⇒　l0＝hdm／tanφ

=== 簡便法 ===

【図-3.4参照】

●堆砂肩は架空の断面上に形成される　⇒　計算上は存在しない
●砂防えん堤地点：断面0において、２つの河床高（堆砂肩高hz0、断面0の河床高hd0）が算出されることになる。
●堆砂肩高hz0は、これまでの計算で断面0の河床高として算出されていたもの
断面0の流砂量V0と断面1の流砂量V1から算出される河床高なので、これが堆砂肩高をあらわすと考えた。
●断面0の河床高hd0は、直下流の断面（断面－1）の河床高の上昇分⊿hd-1と断面0の初期河床高hd0初期値より算出する
断面－1の河床高hd－1は、断面－1の流砂量V－1と断面0の流砂量V0から算出される断面－1～断面0間の河床高であるため、出水後半に透過型砂防えん堤から流出する土砂に伴う河床上昇を表現できると考えた。
なお、断面－1の河床高hd－1は、これまでどおり、断面－1の流砂量V－1と断面0の流砂量V0から算出する。

①砂防えん堤直上流断面（断面1）の流砂量V1、河床高hd1を算出
②断面0の水深hw0と前ステップの堆砂肩高hz0' を比較
　　hw0＞hz0'　の場合　③へ
　　hw0≦hz0'　の場合　⑤へ
　　※水深hw0はスリット形状を用いて算出
－ － － －－ － －－ － －－ － －－ － －－ － －－ － －－ － －－ － －－ －
③断面0のエネルギー勾配i0を算出
　　※（流量／断面0の砂防えん堤がない場合の断面積）の流速を用いる
④断面0の流砂量V0、堆砂肩高hz0を算出
－ － － －－ － －－ － －－ － －－ － －－ － －－ － －－ － －－ － －－ －
⑤堆砂肩までの距離l0、勾配i0aを算出
l0＝hz0／tanφ　　　　φ：水中安息角（30～35°）
i0a＝（hz0－hw0）／l0
⑥断面1の限界水深でのエネルギー勾配i0bを算出
⑦勾配i0aと勾配i0bを比較し、大きいほうを断面0のエネルギー勾配i0とする
⑧断面0の流砂量V0を算出
⑨堆砂肩高hz0を算出
－ － － －－ － －－ － －－ － －－ － －－ － －－ － －－ － －－ － －－ －
⑩堆砂肩高hz0と砂防えん堤のダム高hdmを比較
hz0＞hdm　の場合　⑪へ
hz0≦hdm　の場合　⑬へ
－ － － －－ － －－ － －－ － －－ － －－ － －－ － －－ － －－ － －－ －
⑪堆砂肩高hz0をダム高hdmと置き換える
⑫流砂量V0の補正
－ － － －－ － －－ － －－ － －－ － －－ － －－ － －－ － －－ － －－ －
⑬堆砂肩高hz0はそのまま
－ － － －－ － －－ － －－ － －－ － －－ － －－ － －－ － －－ － －－ －
⑭砂防えん堤直下流断面：断面-1の流砂量V-1、河床高hd-1を算出
⑮河床高hd-1の差分⊿hd-1の正負を確認
　⊿hd-1＞0　の場合　⑯へ
　⊿hd-1≦0　の場合　⑰へ
－ － － －－ － －－ － －－ － －－ － －－ － －－ － －－ － －－ － －－ －
⑯断面0の河床高hd0の差分⊿hd0を⊿hd-1と置き換える
－ － － －－ － －－ － －－ － －－ － －－ － －－ － －－ － －－ － －－ －
⑰断面0の河床高の差分は0　
－ － － －－ － －－ － －－ － －－ － －－ － －－ － －－ － －－ － －－ －
⑱断面0の河床高hd0はhd0の初期値＋差分⊿hd0




hw0：スリット形状を用いて算出（スリットの底面高さは前ステップのhd0を利用）i0（エネルギー勾配）、Ｖ0、hz0：河道形状を用いて算出
図-3.4　簡便法の概念図
*/

===== 支川の最下流端水位 =====

河道の最下流端水位の設定は、ダム地点の設定方法を利用できる。しかし、別の河道に合流する支川の最下流端水位に合流河道の水位を反映したい場合がある。

この時には、断面データ作成用のマクロの「合流状況データ」シートで、合流する河道の水位考慮の有無を指定することができる。この時、合流河道の水位を考慮するように指定した場合、支川の最下流端水位は、ダム地点設定よりも優先される。