国立防災科学技術セソター研究報告

1975年8月

第12号

551・579・1：551・482．4（282．2）

タンクモデルによる非湿潤地帯河川の流出解析（付：

ビキン河・木津川の流出解析）

菅原正巳

尾崎容子・渡辺一郎・勝山ヨシ子
国立防災科学技術セソター

Rmo肚A皿alysis

of

the

Rivers

in

by

Mode1

l■≡己egioms

Ta皿k

No皿一humid

By

M・Su喫w趾a，E．Ozaki，I．Wat㎜abe㎜dY．Katsuyama

肋肋伽Z地8θακんC刎θけ07跳α8妙肋閉〃o〃，τ物o
Abstract
As

it

rains

structure

theof

soi1isfactorsmakebasins，anddone

a11theriversyear

wet．Innon−humid

analysis

is

the

the

tank

modiicationof

additionsomeofround，it

without

runoffany

ispossibleconsideration

toofregions，however，one

prob1emofmode1forstructureana1yse

of

the
it

is

necessary

ana1ysisofnon−humid

the

rmoff

corresponding

soi1moisture

the

tank

mode1mentioned

above．But

it

is

not

enough

of

the

soil

is

wet

on1y

in

some

a1ong

the

river

whi1e

part

the

soi1is

occur

only

from

can

model

is

〃x刎types

are
be

can

horizontal
dryresults

about

the

wet

area．Toofsimulatethevertica11y

obtained，asshowntransfered

than
some

of

both

in

groundwater

the1owerbasinszone．By
of

in

this

zones

of〃tanks1aid

movement

ear1ierremaining

performed，eachcomposed

water

the

in

to

part

dry．Thenstate，division

being

most

to

state

in

represented

series．Thus，the

the
basin

the

surface

runo丘

of

the

basin

by

the

tank

modelsthis

directions，andthehigher

zone

the

this

model

it

is

possiblebeexpress

horizontal

on

to

of

tank

Fig．8，where〃＝刎＝4．In

importantshou1dvertica1and

non−humidruno丘

soi1moisture．Soforbasin，wherethesoi1moisture，because

fo11owing

into刎zones

1．

Japan

always

for

some

in

Japaneseofmode1，
bythemountain

side

tums

to

fair1ygood

getregions．

はしがき
1・1菅原は長年，図1の型のいわゆるタソクモデルを用いて，わが国の諾河川の流出解析を行ない，多くの河川につき，洪水，渇水を通じて，かなりよい結果を得ている．1）・2）・3）

このモデルを，海外の条件の異なる諸河川に適用できるか，その際どのような修飾が必要であるか，資料を得て検討する機会を得たいと，かねてから考えていた．
一

1一

国立防災科学技術セソター研究報告

1．2

第12号

1975年8月

昭和35年の秋，菅原が英文論文をStanfo・d大学のLinsley教授に送ったとき，折り返し，彼の方でも同じような考え方を始めていると，
C・awfo・dとの共著になる小論文が送られて来た．4）それはStanfo・d大学の近くを流れるLosTrancosCreekという小さい河の流出解析で，毎目の雨
から毎日の流量を算出することを目的としている．菅原・勝山はそれに出ている小さな図から1年分（図は1年分だけであった）の雨と流量とをよみとり，タソクモデルに修飾を加えてその流出解析を行ない，その結果をLins1ey

教授に送るとともに，さらに解析を進めたいから，Los

T・ancosC・eekの資料を送ってほしいと頼んだ．やがてLinsley教授から10年分のPa1o

A1toの雨量，蒸発量の資料が届いた．ただし，流量はU．S．G．S、の印刷資料にあるから，日本で探してほしいということであった．何分にも戦前1931〜

1941年の資料で，ようやく，地質調査所にあることを知り，資料をそろえ（1年欠），解析にかかったのは昭和36年秋であった．
1．3LosTrancos

Creek流域の雨期は冬季で，夏季半年はまったく雨が降らない．そして雨期になると，！00mmぐらい（目雨量100mmではな
い，何日かの総計である）の雨が降っても，すべて欠損雨量となり，河川に

図1

タソク
モデル

水がまったく現われないことがある．これはわが国の河川では見られないことである．

タソクモデルでは，多くの場合図2a）のように，側面の流出孔が底面よりいくらか上に（10mm〜15mm程度のことが多い）ついているから，雨量が小さいと側面流出孔からの流出はない．しかし，この欠損はたかだか！5mm程度で，目雨量15mm以上の雨があれば，河川に増水が見られる（火山灰地
帯で浸透が大きい所はこの限りでない）．それが図2a）のモデ

ルで，側面の流出孔を底面から10mm〜15mm程度の所につo）ける理由である．いま，日雨量数mm程度の雨が断続して降り続き，図2a）のタソク内の水面がいつも側面の流出孔の高さに到達しないとする．その場合，相当量の雨が降っても側面b）

図2 a）貯留型．b）完全欠損つき貯留型．

の流出孔からの流出はないが，底面の浸透孔からの水によって，2番目，3番日のタソク内の水位が上がり，そこからの流出が始まる．つまり，河川の流量は徐々に増大する．このよ

うにして，わが国では100mmの雨がすべて欠損雨量となることは起こらないのである．
1・4浸透も流出もしないような欠損雨量に対して，図2b）のような構造を作ればよいことを，菅原は以前から考えていた．この構造の点線以下に貯えられた水は，浸透も流出もしない．それは完全な欠損雨量で，蒸発散によって失われるだけである．

一2一

タソクモデルによる非湿潤地帯河川の流出解析一菅原・尾崎・渡辺・勝山

1．5菅原・勝山は昭和32年に宝川の洪水解析を行なった際に，図2b）の構造を使ってみた．そして，完全欠損部分を3mm〜5mm程度にすると，幾分よい結果が出るように思った．しかし，この程度の完全欠損であれば，あろうとなかろうと，結果にあまり影響しな
い．そこで，その後，図2b）の構造は利用しないことにした．

1・6LosTrancosCreekの大きな完全欠損を見て，図2b）の構造を利用すればよいと考え，まず図3のモデルを作った．各タソクに大きな完全欠損部分を設ければ，よい結果が出るであろうと期侍したが，うまく行かない．雨期になって，しばらくたつと，小さな表面流出が現われる．それをLins1eyは非浸透地域の考えで説明する．流域の一部が非浸透となり，そこに降った雨は表面流出する．雨期が進み，雨が降り続くに従って，非浸透地域の面積が拡大するとして，表面流山の仕方の変化を説明する．

1．7菅原・勝山は，この非浸透地域の考えをモデルに取り入れ，次のように考えた．非浸透地域は河沿いに発生し，次第に拡大すると考えてよかろう（図4のλ）．たとえば図4のBのように，河からとび離れた地域が非浸透になったとしても，3からの表面流出は周辺地域で浸透し，表面流出としては河に届かない．つまり非浸透地域としてBは無効である．したがって，非浸透地域は河沿いから拡大して行くと考えてよい．

非浸透地域の拡大を連続的に考える代わりに，図5に示すように流域をいくつかの地帯に分割し，乾期になると流域は山側の地帯から順に乾き，雨期になると川沿いの地帯から順に湿ると考えた．それは連続物体の力学を，質点系の力学に置き換えて考えるようなものであり，また近時流行の有限要素

1．8

、

、

ソクモデ

このようにして図6の構造ができた．並列に置かれた1段目のタソ

■

図3完全欠
損つきタ

法的な考えである．

、

1

一一・、＼

、㌧

1

ノレ

一、一一、

．．．．

．
一

一

、
、

一

、

1

、

f

、

1
1

、
、

1

一

一

一

、

、

、

、

、

、

図4非浸透地域

図5地帯分割国立防災科学技術セソター研究報告第12号

1975年8月

図6

クは各地帯を表わす．ここで山の側のタソクか

図7ら，河の側のタソクヘの水分の水平方向移動を考えるならば，1段目のタソクは，乾期になると山の側から乾き，雨期になると河の側から湿る．これについての報告原稿は印刷の途中で紛失し，しかも年度末の多忙のため控えの原稿もなく，現在その細部は不明となった．

1．9
Nam

その後1年ばかりして，昭和39年に菅原・勝山はメコン河の支流でタイ国を流れる

Muneの小さい支川にあるSae機会を得た．このときLos

Trancos

Fa11s（流域面積2，060km2）地点の流出解析を試みる
Creekのことを，あわせて考え直してみることにした．

前図，図6のモデルで，土壌水分（完全欠損部分の水）の水平移動を仮定したのが，不白然であって，気にかかっていたからである．

そこで思いついたのが図7のモデルであった．それは完全欠損を表わす構造を持つタソク（図2b）を縦と横とに置いたものである．
1・10

図7のモデルのそれぞれの縦の列は，図5の各地帯S1，∫・，S・，、8・に対応する．各地

帯は地表および地下の何段かの帯水層を持っている．各地帯の各帯水層は，その地帯の上の帯水層からの浸透水を受け，下の帯水層へ水を送る（水の鉛道方向移動）とともに，山の側の地帯の対応する帯水層から水を受け，河の側の地帯の対応する帯水層へ水を送り出す（水の水平方向移動）．この場合，完全欠損に当たる部分の水は移動しない．以後，完全欠損に当たる部分の水を拘束水，動き得る水を白由水とよぶことにする．
一

4

一

タソクモデルによる非湿潤地帯河川の流出解析一菅原・尾崎・渡辺・勝山

1．11蒸発は，まず白由水が優先する．1段目のタソクの拘束水をそのままにして，2段目以下のタソクの自由水から先に蒸発させるのは一見不思議であるが，土壌水分が蒸発すれば，下層の自由地下水によりその穴が埋められると考えれぱよい．

ある地帯の白由地下水が消失すると，士壌水分からの蒸発が始まり，その地帯は乾き始める．この場合，河の側の隣接地帯の白由地下水からの補給は考えない．白由水の重力による

水平方向移動を考えるのに対し，蒸発について水平方向移動を考えないのは，一見つり合いを欠くようであるが，蒸発による土壌水分の穴を下層自由地下水が埋めるのは，毛管現象による吸い上げである．鉛直方向はm単位，水平方向はkm単位と，長さの程度がけた違いであることを考えれば，毛管現象による地帯問の水の移動を無視し得ることが納得できるで
あろう．

1．12

山の側の地帯の白由地下水の一部は，河の側に水平移動して行くから，蒸発により白由地下水が消失して行けば，山の側から乾き始める．同じ理由により，雨期になれば河の側から湿り始める．したがって，雨期になると非浸透地域（土壌水分が飽和した地域）が河の側から増大し，乾期になると山の側から減少して行く様子を，図7の構造で白動的に表わすことができる．

土壌水分の移動を考えず，自由水の移動だけから，非浸透地域の拡大・縮小を表わせる所が，この方式の長所である．この方式により，Los

Trancos

Creek，Sae

Fa11sに対してかな

りよい結果が得られた．5）

1．13

われわれはこのモデルによる非湿潤地帯河川の解析例をふやしたいと思っていたが，なかなか資料が得られなかった．今回WMOが行なった「水文学における概念的モデ
ルの相互比較」（Intercomparisonofconceptua1mode1s

inhydro1ogy）に関連し，つぎの非湿潤地帯4河川の資料を得ることができた．

Bird

Creek：U1S，A．オクラホマ州北部流域面積2，344km2．

Wo11ombi
HunterSanagaBrook：オーストラリア，シドニーの北方，ニューカッスルの西方にあり，
Va11eyの支川，流域面積1，593km2．

River：アフリカのカメノレーソ，流量観測地点は河口に近く，流域面積131，500km2，⇒ヒ緯5。のあたり．

N・m

Mune：タイ，NamやMeは河の意，メコソの右側支川，流量観測地点ウボソ，流域面積104，000km2，北緯15。のあたり．

これら河川の流出解析にあたり，モデルにはいくらかの修飾，変更が行なわれた．

2．非湿潤地帯の河川に対するタンクモデルの構造

2．1地帯分割
流域を図5のように舳個（今回は舳＝4）の地帯に分割する．その面積比を∫・：S・：一

5一

国立防災科学技術セソター研究報告

S1

xSl

第12号

1975年8月

S肌とする．S1が山の側，S㎜が河沿いである．流

域の種々の状況を知っていれば，それによって面

、

S2

xS2

積比を定めるべきであろうが，何も知らないとき＝＝＝⇒xS。

は，試行錯誤で定めるより致し方ない．この面積

↓
S，

、

一→xS
→
→

比をうまく定めることは，よい結果を得るための大切な条件であるらしい．およその方針として，
S1，S・，…，S㎜が等比数列となるように定める．

たとえば
33：32：3：1＝67．5：22．5：7．5：2．5，
2．53＝2．52：2．5：1＝125：50：20：8

→xS4

≒61

：25：10：4，

23：22：2：1＝53．3：26．7：13．3：6．7．

一→xS4

2．2π×㎜型のタンクモデル

1）各地帯を〃段（今回は〃＝4）の直列貯留型モデルで表わす．このようにして，図8に示す↓〃×舳型のタンクモデルができる．この図で左が

図84x4型タソクモデル

山の側，右が河の側である．水は上から下への鉛直移動と，左から右への水平方向移動をする．

2）第1段目のタンクだけが土壌水分に関する2段の構造を持っている．この所が10年前の図7のモデルと異なっている．以前のモデルでは，各タソクの完全欠損部分の構造は1段で，その代わり2段目，3段目のタソクにも完全欠損の構造を考えていた．土壌水分を2段にするのはLins1eyの考えによるもので，彼は上層土壌水分（upPer

soi1moisture）と下層土壌水分（1OWe・Soi1mOiSture）とを考える．近時のサクラメント・モデル6）毛同様の考え方をしている．

3）1段目のタ：／クの側面の流出孔は2種類あって，a）下の流出孔からの水は右側の1段
目のタソクにはいる．b）上の2個の流出孔からの水は直接河川への流出となる．今回の解析結果では，Bi・dCreek，Wo1lombi

孔だけである．一方，Sanaga

Brookではa）の型の流出孔は無く，b）の型の流出

River，Nam

Muneでは，a）の型の流出孔だけで，b）の型は無い．a）の型の場合，一番右の河沿い地帯のタ：／クでは，1段目タンクの流出孔からの流出は，河にはいるのである．

4）〃x刎型のタソクモデルで，各タンクの構造をすべて独立に変化させると，動き得る

パラメータの数は50個をこえる．このような多数のバラメータを制御することは，人間の能力をこえるであろう．しかしながら，〃x舳型タンクモデルの1段目，2段目，…

のタ

ンクは，それぞれの帯水層に対応すると考えているのであるから，同一の構造を与える．つ一

6一

タソクモデルによる非湿潤地帯河川の流出解析一菅原・尾崎・渡辺・勝山

まり従来のタソクモデル（〃×！型）を，パラメータを固定したまま舳列に並べたものを，〃×〃型のタソクモデルとする．ただし1段目の土壌水分の構造だけは，地帯別に少し変えることがある．したがって，〃×刎型であるからといって，従来の〃×1型のタンクモデルに比べ，それほどパラメータの数が増すわけではない．地帯面積比，±壌水分についてのパラメータが増すだけである．目本の河川で，積雪・融雪の機構を導入するには，高度（寒さ暖かさという方がより適しているであろう）による地帯分割が必要で，それにはかなりのパラメータが必要である．雪の場合に比べて，〃×舳型のモデルのパラメータはそれほど多数
ではない．

5）〃×舳型のタソクモデルのパラメータがそれほど多くならないのは上述のとおりであるが，〃×舳型のタンクモデルのパラメータの数がかなり多いから，それを適当に加減するならば，実測流量に合わせることができるのは当然であると考える人もあろう．実際には，

そうならない．〃×刎型のタソクモデルという，一つの型を想定すれば，それによる強い制約が生ずる．パラメータを変えるだけで，任意のハイドログラフを出せるわけではない．それは従来の〃×1型のタソクモデルについてもいえることである．一つのモデルの型を定めることにより，自由さはきわめて限定される1

たとえば〃×1型のタソクモデルで日流量解析をする場合についていえば，1段目のタソクの半減期は1目〜3目程度，2段目のは5日〜10目程度，3段目は1月〜3月程度，4段目は半年〜2年程度にとることが多く，それはそれぞれの流出成分に対応している．もちろんタ：／クモデルは有限要素法的な近似の意味を持つもので，現実の流出成分や地下の帯水層との対1一ごがはっきりつくものとは思わないが，ともかく各パラメータに対する白由度はかな

り限られている．逆にいえば，白由度が限られているから，パラメータを深し求めることができるのである．

タソクモデルにかなり強い制約があり，どんなハイドログラフでも出せるわけではない所に，タソクモデルの一つのよい点があると考えられる．与えられた雨量から与えられた実測

流量がどうしてもうまく出て来ないときは，流量に誤りがあるのではなかろうか（水位流量曲線に欠陥があることもあるし，水位計が故障していることもある），雨量に問題があるのではなかろうか，流量や水位に何かの人為的操作が加えられているのではなかろうか等と考え，その原因が明らかになったこともしばしばあった．与えられた雨量から与えられた流量を出すのだけが，流出モデルの任務ではあるまい．流出モデルの側から，雨量・流量のチェックをするのも，流出モデルの任務の一つである．そのとき，モデルの持つ不白由さ，制約が重大な役害■1を演ずる．それがモデルの重要な意義の一つではあるまいか．

6）今回得られたモデルのパラメータは，従来口本の河川で得られたものと，似たりよったりであった．もちろんSanaga

River，Nam

Muneのような川は13万km2，10万km2という，わが国では経験できない大流域であり，かかる大流域では第1段のタソクの半減期

一7一

国立防災科学技術セソター研究報告第12号

1975年8月

は長くしなければならないのであるが，それについても，半滅期と流域面積の平方根はほぼ比例関係にあるという従来の経験則から，あまり離れていなかった．流域を地帯に分割し，ある地帯では地下水が消滅し，土壌水分が乾き切るという構造さえ導入すれば，あとは日本の河と似た構造でよいのである．なお，わが国では，四季を通じて雨に恵まれているから，流域の全体がつねに湿っており，したがって〃x1型のタソクモデルで近似できるのである．

2．3土壌水分
1）地面の表層部分に，土の粒子間に空気と水とが共存する層，いわゆる不飽和帯がある．

これの状態が，水の浸透，表面流出に深く関係するのであるが，それの簡単なモデルを作るのが，ここの目的である．

2）雨が降って不飽和帯に水が供給されたとき，水の一部は土の粒子間の狭い空げきに入
りこみ，一部は粒子問の広い空げきの部分をみたす．この両者が完全に区別できるかどうかには疑問があろうが，前者が拘束水，後者が自由水である．白由水は重力により下方に浸透する．

3）拘束水を1次と2次とに分ける．不飽和帯に供給された水は，まず入りやすい狭い空
げきを満たす．これが1次拘束水である．これが飽和すると，余った水は白由水となりその一部は浸透または表面流出する．

4）

1次拘束水は徐々に，より入りにくい狭い空げきに移動して，2次拘束水となる．今回のモデルでは，！次拘束水は50mm程度，2次拘束水は250mm程度で飽和すると置いた場合が多い．

1次拘束水から2次拘束水への移動速度ムは2次拘束水が少ないとき（乾いているとき）に大きく，それと2次拘束水の貯留高X。との関係は直線的であると仮定し（図9a）），次式
で与えられるとする．

ム＝oo＋c（1−X2伽）．

すなわち，2次拘束水oのとき（完全に乾いているとき）移動速度はム＝c。十6（mm／目），
2次拘束水飽和のとき（X2＝ん。），ム＝ooである．今回はo。＝0．5，c：1とした場合が多い．

もちろん2次拘束水が飽和すれば移動は止まる．

5）

1段目のタソクに自由水があれば，1次拘束水の穴は自由水で埋められる．1段目タソクに白由水が無く，2段目のタソク以下に白由水があるときは，1次拘束水に生じた穴は，下層の白由水の吸い上げによって埋められる．その移動速度ムは，1次拘束水の貯留高X・の1次関数として，次式で与えられると仮定する（図g，b））．
ム＝6o＋6（1−X1／乃1）．

すなわち，1次拘束水がO（完全に乾いている）のとき，下層からの吸い上げは伽十6（mm／目），1次拘束水が飽和のとき（x1＝加），あ。（mm／日）である．もちろん飽和すれば吸い上げ

一8

一タ：／クモデルによる非湿潤地帯河川の流出解析一菅原・尾崎・渡辺・勝山

I2、．lll…ザニ
X2h2｝b．L・．．．．．、㌧

…
●．

X1
O

図9

h、

土壌水分についてのモデルは止まる．今回は仇＝5：3（mm／日）を用いた．

6）1次拘束水，2次拘束水はLins1eyその他の上層土壌水分，下層土壌水分に相当するものである．そもそも土壌水分を2層に分けることも，有限要素法的な近似であろうが，雨水が土壌水分（その主要部分は2次拘束水である）に移行する際に，中間段階が必要なようで，それが上層土壌水分，ここにいう1次拘束水である．ここで上層・下層の名を避け，1次・2次を用いたのは，それが必ずしも上下の関係にはないと考えるからである．雨が降って上層の土壌が十分に湿れば，下層が乾いていても表面流出が生ずる．その意味で上層土壌

水分の名がつけられたのであろう．しかし表面流出以外に浸透の問題がある．浸透は土粒子の間げきを通る水の下方移動だけではあるまい．それだけであれば，浸透の速度はかなり遅いものであろう．しかし，地面にはひび割れや，孔や，断層など，上下をつなぐ大きな抜け道がある．水のかなりの部分が，これらの抜け道を通って下方へ移動するのではないだろう

か．乾いた地面に雨が降ったとき，これらの抜け道を通って雨水が浸透できるのは，抜け道の周りが十分に湿ってからであろう．それは土地の表面が十分に湿ってから表面流出が生ずるのと同様のことである．これが上層・下層の名を避けた理由である．

1次拘束水は自由水との接触部分であるという考え方から，下層の白由水の吸い上げも1
次拘束水を経由することにした．

わが国は四季雨が降るから，土壌水分はつねに飽和に近い状態にあるとみてよい．日本で一

9一

国立防災科学技術セ：■ター研究報告第12号

1975年8月

畑が乾いたというのと，乾燥地帯で畑や牧野が乾いたというのとは，かなり程度が異なるものらしい．その意味で，われわれは土壌水分の欠損による，大きな初期欠損雨量の実例を見る機会がない．不幸にして，今回資料として用いたBi・dCreek，Wo11ombi

B・ookはともに，流出解析用の資料としては最適でないように感じられる．したがって，これらの資料か
ら作った今回の土壌水分に関するモデルは，そこで仮定したバラメータの値とともに，信頼性にはやや欠ける所があろう．
2．4

蒸

発

散

1）蒸発散は地表から行なわれると考え，第1段のタソクから引き去る．蒸発計による目
蒸発量の観測値Eが与えられているとき，第1段タソクの白由水からはO．8Eを差し引き，
白由水が無く拘束水から差し引くときはO．6Eを引く．白由水の残りが少なく，それから
o．8Eを引き切れないときは，不足分の3／4倍を拘束水から差し引く．差し引きは，白由水，

1次拘束水，2次拘束水の順である．

2）

2段目以下のタソクに自由水があるときは，蒸発散の差し引きにより1次拘束水に生じた穴は，下層の白由水の吸い上げにより埋められる．したがって，拘束水より先に，下層の白由水が蒸発するのと似た結果が現われる．

3）ある地帯がすっかり乾けば，その地帯からの蒸発は無くなる．上記のように，蒸発計の値に対し，白由水からはO．8E，拘束水からはO．6Eと地面からの蒸発を小さくしてあるのであるが，流域の一部が乾いてしまえぱそこからの蒸発は0となるから，流域全体からの蒸発はEに比べはるかに小さくなる場合が生ずる．
2−5

河道による変形1）〃×舳型のタソクモデルで雨量から変換された流出高は，河道に入ってから，変形を受ける．

2）河道による変形は，図10の貯留型構造により与える．Sanaga

River，Nam

Muneに

ツ＝oX▲

2．

y：2oX．X一αX。
◎く2

・↑［改．

o）

え

y：αX

八h「

x
O

b）

X．
C）

図10 河道による変形に対するモデル

ー10一

タソクモデルによる非湿潤地帯河川の流出解析一菅原・尾崎・渡辺・勝山

はa）を用い，BirdCreek，Wo11ombi

Brookにはb）を用いた．b）は底部に初期欠損を持ち，その部分を除けば，流出高は貯留高の2乗に比例するものである．

3）流出を貯留高の2乗に比例させると，計算技術上の難点が出て来る．たとえばμ＝0．1×2のとき，X＝10でひ＝10となり，X＞10でμ＞Xとなる．この式をそのまま用いると，貯留高以上に流出が生ずるという不合理になる．流出が貯留高の2乗に比例するアナログ構造を作り，実際に運転すると，流出して貯留高が減れば白然に流出が減るから，貯留高以上に流出することはない．ディジタルで計算しても，適当に時間区間を細分すれば，アナログと同様の答が出て来るのであるが，河道による変形の計算に，それほどの手間をかけることもない．貯留高が大きいときは，貯留高の大部分をその日のうちに流出させればよ
いので，そのとき河道による変形効果が小さくなるようにして置けばよいのである．
そこで，便宜上，次の方式を用いる．放物線ひ＝αX2上の1点（X・，αX・2）で接線を引き，この点で放物線と接線とをつなぎ合わせた図形（図10，c））で，ひとXの関係を与えることにするのである．

1−／；㍍

：㌶

たとえば2αX。＝1．0によりXoを定めれば，X＞Xlでひ＝X一αX12であるから，入力が大きいとき，入力はそのまま出力となり，貯留効果はまったくない．いくらか貯留効果を持
たせ，ハイドログラフの平滑化作用を持たせるためには，2αX。を0．9とカ㍉0．8とかに置けばよい．今回は2αX・＝O．8と置いて，図10c）の曲線を作った．

4）Bi・dC・eek，Wo11ombi

Brookでは河遣（河に沿った帯状の狭い地域で，河川表流水の供給により地面が水分を持っている地域を考える）への降雨，河道からの蒸発の影響を計算に入れる．今回の計算では，両河川とも，河道の面積を流域のO．4％とした．これを0．2％〜O．5％程度にすると，河川が極度に渇水になったとき，河道への降雨，河道からの蒸発が，河川表流量にかなりの影響を及ぼす．とくに蒸発の影響が大きい．

たとえば，地下水から河へo．05mm／日の水が出て来たとする．これが河道に水分を供給し，そこから5mm／日が蒸発するとする．河道の面積が流域のo．4％とすると，5mm／日の蒸発は流域全体で割ればO．02mm／目で，これを0．05mm／目の表流から引けば，残りはo．03mm／日となる．もしも蒸発が8mm／日であれば，o．018mm／日しか残らない．このように，渇水になると，表流は蒸発の影響を受けて大きく変動する．影響の受け方は，河道の流域に対する面積比を0．5％にするか，O．1％にするかで大きく変わるが，河道への流出高がo．05〜o．01mm／日程度になると，河道からの蒸発の影響は大きい．Bird
Wo11ombi

Sanagacreek，

Brookでは，この影響と思われる渇水流量の変動が見られる．

River，Nam

Muneの渇水はそれほど小さくないから，河道に対する降水，蒸発は無視する．河道による変形を図10a）の簡単な構造で与えることにしたのも，渇水があま

一11一

国立防災科学技術セソター研究報告第12号

1975年8月

りノ」・さくならないからである．
マ

O．35

3．計算の仕方

等4〔。。

3．1計算の仕方や，タソクモデルの図に記

或入されたパラメータの意味は，ふつうのタソク

15

モデル（〃×1型）の場合と同じである．計算はすべて雨量に換算した水深，mm単位で行な

0

60

mm

う

たとえは図11の構造で・拘束水は飽和・

α）

その上に20mmの白由水がある所に，30mm

0．25（≦15）の雨が降り，蒸発量が5mmであるとする．蒸図11
え，

発散5xO．8＝4（mm）を引き，雨30mmを加白由水の貯留高は46mmになる．これから流出高，浸透高を出すには，流出孔，浸透孔にかかる水圧（水深換算mm）に係数を掛けて，次のように得られる．
流出高：（46一（25＋10））×0．35＋（46−10）×O．25＝12．85（mm／日）．

浸透高：46×o．25＝11．5（mm／日）．
貯留残高：46−12．85−11．5＝21．65（mm）．

図11の浸透孔の所に≦15とあるのは，浸透高が15mm／目で限界に達することを表わ
す．すなわち，貯留高60mmまでは，浸透高が貯留高に比例するが，以後は15mm／日の一定値となる（図11，a））．

3．2今回の計算方式で，計算は次の手順で行なう．1）蒸発散の差し引き，2）下層からの白由水の吸い上げ，2次拘束水への移行，3）雨量の加算，4）流出，浸透の算出，およびその差し引き．

以上は1段目のタソクの場合である．2段目以下のタソクでは1げの差し引き，2）上段，山の側のタソクからの水の受け入れ）1次拘束水への吸い上31）下段，河の側のタソク

ヘの水の移行，の順になる．

これは計算上の約束にすぎない．手順を変更すれば，答はいくらか異なる．どの手順がよい■とも言えないが，計算を実行する以上，規則を定めて置かなければならない．

．
1

1』・・．■1・

1

■

、 ：
1L ・

一
・I

・

3．3〃x舳型のモデルの計算で注意を要す、

るのは，隣接地帯との水のやりとりである・各
地帯のタソクの水はそれぞれ水深（mm）で計算している．それは単位面積あたり水深である．

同一地帯内の上下の水の移動はそのままでよい

一12一

。且1
↓

：
L

S1

←二圭

・…
l1

1

図12地帯問の水の移行タソクモデルによる非湿潤地帯河川の流出解析一菅原・尾崎・渡辺・勝山

が，4番目の地帯から，河の側の隣接地帯（ゴ十1番目）のタソクヘの移行では，S包の面積
から出る水であるから，量としてはS。倍しなければならない．それが面積S｛。1の地帯に入るから，そこの水深に直すにはS。。1で割らなければならない．結局，地帯クから地帯（ゴ十1）への移行にはS1／S・十1を掛ける必要がある．

同様の理由から，ク番目の地帯のタソクから河道への流出には8。を掛ける．

3・4

河道による変形を与える構造では，流域あたりの水深mmで計算している．それに対し，降雨，蒸発の影響を計算に入れるときは，河道の面積∫。（流域面積を1とした比）を
降水量，蒸発量に掛けたものを，加え，または差し引く．

4．得られた結果，流出機構探求の経過，および残された問題

4．1

表1は各河川について得られた流出機構のバラメータを示す．パラメータの意味は

図13に示されている．

図14は得られた結果の一部を示している．結果のすべては英文報告7）に図示してあるか
ら，関心のある方は，それを見ていただきたい．

図14の縦軸は対数尺になっている．Bird

Creek，Wo11ombi

Brookでは流量が0になることがあり，そのまま対数尺に刻むと，負の無限大になり，流量の小さい所でハイドログ

ラフが拡大されすぎて，具合が悪い．この欠点を除くために，ある小さな2・を流出高σに加え，（σ十ωを対数尺に刻むことにした．卯はBird

creekで103mm／目，Wo11ombi

Brookで10−2mm／日である．その上，図14では1og（σ十σ1）の関数尺の上に，4の目盛が刻んである．したがって，σが大きい所ではほぼ対数尺，αがOに近い所で普通尺に近づくという目盛になっている．

4・2解析はBird

Creekから始められた．最初は図7のモデルを用いたが，1段目のタンクに2層の土壌水分を置く構造の方がよいらしいことが分かり，図8の形となった．

4．3

かなり長い間，図8の構造で，1次拘束水から2次拘束水への移動は考えるが，下層の白由水から1次拘束水への吸い上げを考えない方式を用いていた．そして，蒸発を引くとき，白由水がある限り，拘束水には手を付けず，白由水から蒸発させる方式をとった．この方式は物理的にやや合理性を欠くと考え，下層の白由地下水から1次拘束水への吸い上げを考えることにし，蒸発は1段目からだけ差し引くことにした．この変更による実質的影響，
つまり結果に及ぼす数値的変化は，わずかなものである．

4・4土壌水分の飽和限界，水の移行速度等に関する定数の大部分はBird

Creekの解析から定めた．これらの定数のあるものは，その値をかなり変化させても，計算結果にあまり大きな影響を与えないようである．つまり鈍感な定数である．たとえば，拘束水間の移動，
白由水から拘束水への吸い上げの速度を与える図9a），b）の直線の代わりに，これを定数と置いても，結果はそう変わらない．

一13一

◎H．o

山⑩o．o

○つo．o

寸o．o

雪o．◎

⑩o◎．o

寸o．o

○つoo．o○つooo．o

山oo．o

寸oo．o

旧o．o

Hoo．o

寸o◎．◎

害︑o

Hoo．o

○寸

寸

◎oH

oo

謎

o0N

○寸

○寸1

oo
oLΩ

○旧

囚＝

轟漣斗く
謎

彗目凄昌Φz

一︹一〇目◎＝o奉−ω＞一串伺b0伺⊆碍αつよoo−串よ①ω−O勺占串

ω；雲昌由zo︸︒つ卜㊤○つ○つH卜㊤N○つ○つω

一①．≧串伺b0由目祠ω

HO
蛍

○寸

寸oHLONH㊤

○寸

○山
○旧

よoo﹄囚石昌2iき○卜

寸oH0NH㊤o⑩oo﹂﹂寸⊥十の辛風峠↑．1厘6コ剖コー．本6幸・・尋6

わ︷hr

■

匡

一

ω目コー﹀﹇F口碍之 ︑H唄μ

岬一♪U﹄兵ぶぺφ州枳蝸匿掻倶〜．o．oHo．へーへ心γ︑︵︶▲︑ト山Tfコ糧＝﹁匡ヰ巾﹄o目◎＝◎峯ム①巴O勺上︸よ◎o﹂mH−①＞⁝︸ob0宙⊆伺oつ

当自妻昌田7︷ω嘗○巾．自因○旧Nω〜一一の︷︵N︵く怯黒︶匝鑓H︑鳩貞︵H

三︷・一さ童芸墓︑妄筆一一嚢︑景婁蓬一一γ︑ψLT匿︺︷粛漣く皆黒無川収︵︶U﹄・♪︵ぺ︺︷以担綿機・串置

H︷︑−へ

旧

N①ooH寸ト

儘漣汁く

雪o．o

ooo．o

害︑o

ま寸．o

OH．O

N．o

oH
◎H

○山N

○血

o﹁○つ〜

糟

氏蜜Qく怯黒起H

寸oo︷

◎0N

皆氏黛豊oo起QoつNく匝︒つ上︷×︵︵8蟹8竃︶十︵慢蟹一一十︵︷易︶十一圭蟹︶︸以榊⁝択Q▲∴ト山Hへ＼川祇

H饅

oH．o

O＝

N一．o

○つo

oO

一〇．o

ρ＝

旧Ho．o

〇一

H匡O

く

Ho

Hρ

ωooo．◎

oρ

皿

旧N．o

N葭O

o○つ

◎H

癒漣片く

ざ寸．o

副

oo◎

一︑o

＝Oω

副 因

血oo

No
血No

H．H

トo．H

閉

漫厘抱匡

○ト

旧o

o○つ

0No
○つく■く葭I

糧

HH〜πへ①︶︵囚ω〜囚H︶
薫蜂山H任牝噌僅掻u応◎．HNく石く葭 Hく；く葭
＜損匡
＜幸碧鑓整

旧○つo

0NH0NooH
旧H糟幸 妻些■一†滅o︷撃く！

1︵口 NH ︑﹇へ oo〜rへ㊤︶︵口﹁コ︵く佃口︶皿 鐘 H

因㊤．o ぶるく嵜黒
ぶるく︷皿

串◎o．o

工＼目︺⊆−㊤．oつ心＾倶︶マ缶黛

因①．oぶ月恨yマ﹇t而﹂

輿卜．◎ぶ■債︶べ音黒

串⑩．◎

ぶ会く嵜黒
ぶるく丑皿

絆

串oo︐o

，

︑o．HHo ︐o．◎Hoo

14

榊

HHO
⊂oooo

︹︑昌昌oo︑寸ぶ月債︶﹁﹇血正墨．心嶋p oつH︷Ho亀

榊

1 1

匡

1 1

脳

匿

■

蝦
司
串

■

糟

1 ■

■

綱

■
織

1975年8月
第12号

国立防災科学技術セソター研究報告タソクモデルによる非湿潤地帯河川の流出解析一菅原・尾崎・渡辺・勝山

1通

↓

↓

冠

A

type

↓

A3→
A2→
A1→

8type

B1→

∴L

（2AXξ0．8）

D1→
DO
図13

4．5

〃×刎型タ：■クモデルのパラメータの意味

われわれのように湿潤地帯の経験しか持たない人問には，土壌水分についてのよい

モデルは作りにくい．1次拘束水の限度を50mm程度，2次拘束水の限度を250mm程度にしたのはLinsleyに従ったのである．場合により，河川沿いの地帯の1次拘束水の限界
を，山の側より小さくしたのは，実測流量に一合わせる必要があるからである．しかし，河川沿い地帯の方が，上層土壌水分の容量が小さいというのは，少し気になる．河沿いの方が，土壌層は厚いのが常識的と思うからである．あるいは，1次拘束水の限界を河川沿いで小さ
一ユ5一

＼1

Um

U
0
Q＿

Σ

L］

L］

oo

〉

0
Z

（j

0二

⊃

ト

0u

︵皇鉢走
トu0

■

■

1

Z⊃﹁

、」

⊃
「

o⊃⊂

卜圧Σ

■

﹄o8因
Z︺﹁

口止H口

口山L

プ⊂Σ

一﹄冒o＝o奉．ノ∵︑ヒσΣHmΣ0﹂﹂0≧

図掴糧皿
﹂⊃﹁

ヒL圧

y00ヒm

包；灘蝋一繁w岬︶

0一繁
出竃

︑ξ

y山山叱U

．︷①賞O巾上因

止圧Σ

一．ミ図

OU口

旺し圧

Z⊂﹁

へη㎝一

口山L

＞0Z

Z⊂﹁

勺ωm一

トU0

①川㎝一

σ

〇一

−N−

ヨ

三

丁

〇

〇一

〇一

N

卜圧口＼ΣΣ

し
○叶

or

m−−

NI

〇一

F−

〇一
〇

〇﹁

一1

〇一
〜

＞⊂0＼ΣΣ

16

1975年8月
第12号
国立防災科学技術セソター研究報告タソクモデルによる非湿潤地帯河川の流出解析一菅原・尾崎・渡辺・勝山■㌔・「「π■□1「「「「「□□「「

←
一
」〔

○山0

U
0
Σ

L⊥」

﹂⊃⊃

山Z︺Σ

﹂⊃﹁

b0︸目伺ω

﹃し江

一㌃U

ヒL工

N．寸一﹇幽

此ぴΣ

ヒ隼一

○山L

︹山L

Z⊆j

フ匡﹂

σ﹁

；−

卜皿0＼ΣΣ

0F

0

0r

一−

O

プ〇ニロ︑︑τΣ

三

じ一

丁

已

ナ

．5主串

戸圧Σ

＼︑︑

卜圧Σ

＼

Φ目目﹈≧ 昌伺宕

Z⊃つ

Σ圧Z

Z⊃﹁

㌔・、

一担；黙蝋一誰一蝋︶ 図蝸填匝

o⊃⊂

o⊃圧

圧O圧Z圧ω

＼

︵運鉢柘
トu0
L］

ω竃

L山oつ

ヒU﹀Hヒ

17

＞0Z

トU0

／︸1 ・＾︹リス︹いい＼﹀0Z／国立防災科学技術セ1／ター研究報告第ユ2号

1975年8月

くする代わりに，2次拘束水の限度を大きくして，全体として，河川沿いで土壌水分を大きくすればよかったのかもしれない．そうすれば，河沿いの方が土壌層が厚いことになるし，河沿いで！次拘束水を小さくしたのは，土の粒子が細かいことで説明できるかもしれない．

元来，2次拘束水の限界を250mmに固定したのは（NamNuneでは150mmにした），

これが鈍感なパラメータだからである．かかる鈍感なパラメータは，流出解析以外の，不飽和帯に関する知見を考慮に入れて定める必要がある．不幸にして，われわれにはそれが欠けているのだから致し方ない．

4．6土壌水分に関する経験がないためにモデルが作りにくいのは事実であるが，土壌水分のことを知っていたとしても，それがそのまま流出モデルに役立つとは限らない．流域の流出モデルはマクロのモデルで，ミクロの性質がそのまま反映するとは限らない．乾燥地帯の多くの河川の流出解析例を重ね，流出モデルを作る側から，土壌水分のマクロ的性質に近づくことも必要であろう．今回の流出解析を一つのきっかけとして，今後世界各地の資料が入手できることを望んでいる．

4．7BirdCreekの解析がある程度進行した時点で，Wo11ombi

Brookに手をつけた．両

河川とも，雨量は地点雨量にある種の加工をし，あるときは雪どけまでを考慮して，流域雨量が与えられていたが，地点雨量がそのまま与えられていることが望ましかった．

雨量の地点による変動は，かなり大きいものであろう．1地点で大雨が降り，他の地点で降らないときは，局地的に表面流出が生ずる．それは平均雨量からでは算出できない．1段目モデルだけは地点雨量で別々に計算するという方式をとれば，地点雨量を用いることによって，もう少しよい結果が出るであろうし，土壌水分の構造によると考えるべきか，雨量の地域変動によると考えるべきか，判断に迷った所も，地点雨量がそのまま与えられていれば，それなりの判断材料として有効であったろう．

4．8地点雨量がほしいと思った，もう一つの理由は，雨量の害11増率を考えるためである．

雨量地点の数がかなり多くても，その大部分が平地にあることが多いから，流域雨量を出すためには，f可かの補正を必要とすることが多い．わが国の多雪地域を流域とする河川では，冬期に大きな割増率を掛け，春から夏にかけては割増率なしにしないと，実測流量に合わせられないのがふつうである．それに加えて，秋には，山地の方が平地より降雨が少ないとみ
なす以外に，実測流量との対比がつかない例もいくつか出て来た．8）山地の雨量（とくに冬期の降雪量）が実測されていない以上，流出解析の立場から得られたこの推測が正しいかど
うか，実証できないのであるが，平地内の観測値の相互比較から，この傾向を証明できる例もある．Bi・dCreekでは，季節別に割増率を変化させる方式をとった．割増，割引きの程度，その季節変化の仕方は，わが国でみられるのと同程度のもので，その点ほぼ妥当であろうと思うが，地点雨量の変動の仕方，その法則性を見ずに，これらを仮定したことが気がか
りである．

一18一

タソクモデルによる非湿潤地帯河川の流出解析一菅原・尾崎・渡辺・勝山

4・9Sanaga

Riverでは13地点，Nam

Muneでは4地点の地点雨量が与えられた．流域面積がそれぞれ13万km2，10万km2であるから，地点により，年雨量の大きさにも，
降雨の季節的変化のパターンにも，かなりの相異がある．

ある地点は雨期の始めに多く降り，他の地点は雨期の終りごろに多く降る．これらの地点の単純平均を流域雨量として用い

ると，あまりよい結果にならない．そこで，流量の季節的変化と，各地点雨量の季節的変化とをにらみ合わせて，いくつかの地点には1／2の荷重を与え，荷重平均をとることにした．

Nam

Muneにおいては，流域外のKhon

Kaen流域の奥のKo・atに1／2，流域の中ほどにあるRoiet，Surinに1の荷重を与えた荷重平均を流域雨量とした．

SanagaRiverでは，13地点のうち，流域外と流域境界上にある4地点に1／2，他の9地

点に1の荷重を与えた荷重平均を用いた．
荷重の与え方にもう少し工夫すれば，結果はいくらかよくなったであろう．Sanaga

Riverでは，流域内部の地点で，そこの荷重を小さくすれば，結果はよくなると思われる地点もあったが，十分の根拠もなしに，荷重を小さくしたり，除いたりする訳にも行かないから，地図の上でながめて，一応無理のない荷重を用いることにした．

4・10NamMune，Sanaga

Riverの解析は，BirdCreek，Wo11ombi

Brookである程度の結果が出た時点で始められた．後から着手した2河川は，解析が楽で，試算の回数が少ない割によい結果が出て来た．とくにSanaga

Riverではよい結果が出て来た．雨量地点数が13と多かったこと，河川が白然状態にあることによるものではあるまいか．

Nam

Mmeで実測，推定がよく合わない所があるのは，雨量地点が4地点であることによると思われる．この地域の雨は，局地的スコールだから，よい結果を出すには，4地点では不足であろう．

5．

ビキン河の流出解析

5．1のBikin流出モデルの相互比較のための対象河川として，上記の4河川のほかに，シベリアRiverと，日本の木津川が選ばれている．どちらも湿潤地帯河川で，ふつうのタソクモデルで解析可能である．その結果もここに付記する．

5．2BikinRiverはハバロフスクとナホトカの中間あたりにある，アムール河の支川であるらしい．流域面積13，100km2，流二㌫ぷ㌫㌶ぷ㌫鴛㍗㎏

資料として，夏期6月〜9月の4か月だけ与えられている．
5・3

われわれは日本でかなりよい結果を得ている積雪・融雪の計算方式をシベリアの河川で試みたいと期待していたので，夏期

ー19一

図15BikinRiverに対するタソクモデル

u山口

〇一

OUω＼Σ
①

﹀0Z

L山ω

Fu0

o⊃圧

︑ノ

﹂ユ﹁

LUω

Z弓

畠江

〇一

〇一

ヒ山﹀Hヒ

﹂⊃﹁

〇一

彗
ヒし圧

石﹄蕾

︑5＞蟹

卜圧Σ

o︺圧

﹁﹁II1−﹁

出ω

ZHyHm

○山ω＼Σ
m

︹

︑㌻

︺ωトHy

ツ

Z⊃﹁

︵埋鉢左一誰一鴬∴埋目灘蝋一繁蝋︶ 図蝸握口＝﹃朴〜共

︑

Nω3

F
卜⊂口＼ΣΣ

ヒ山﹀H此

ヒσΣ

﹂一﹂﹁

㌧r

口山L

75﹁

δ㎝F

z圧﹁

η山2

口﹁

〇一

−

ト⊂口＼ΣΣ

σ一

一﹁−

〇一

〇

○ト

ー

〇一

〜

戸エロ＼τΣ

20〇一m○山ゼ＼．二

ぺ
〇一

N
［

へ

m

1975年8月
第12号

国立防災科学技術セソター研究報告タソクモデルによる非湿潤地帯河川の流出解析一菅原・尾崎・渡辺・勝山のみの資料に失望した．夏期4か月の資料というこ

表2ビキン河初期値
1段

2段

3段

河道

とは，長期流出成分の解析ができないということ

1961

5

ユ5

300

10

で，それは日雨量から日流量を出す計算の難しさ面

1962

10

30

500

15

白さの大半が消えることである．基底流量の半減期

1963

0
O
O
O

30

400

13

20

400

11

は年単位程度のものであるらしいから，何年も連続

70

300

16

して解析してみない限り，渇水の解析はできないの

50

300

14

1964
ユ965
1966

である．

5．4

困ったのは蒸発散である．この流域は，あまり高温でなく，雨が多いから，蒸発散の影響は小さいと考え，月ごとに3mm／日〜4mm／日程度に固定してもよいかと考えたが，せっかく与えられた飽和水蒸気圧欠損や風速を利用しないのも失礼であろう．ところが，飽和水蒸気圧欠損（θrθ），風速〃と，蒸発量Eとの関係を与える実験式
E＝（α十肋）（θ。一θ）の係数が人により，国によりかなり変わるらしい．しかも，この場合，流域の出口で測った〃，（θrθ）から，流域全域についての蒸発散を出す式がほしいのであるから，だれかの式をそのまま用いるわけにもいかない．
（α，6）の値として（0．15，0．8），（O．1，0．4），（0．05，0．2），（0．07，0．3）等を試み，結局

E＝（O．1＋0．4〃）（θ3一θ）を用いることにした．別に根拠はない．

5．5図15は得られた流出機構を，図16は結果の一部を示す．
毎年6月1目からの計算開始であるから，その度に初期条件を与えなければならない．それは実測流量をながめて，およその見当で定めた．表2はそれを示す．

6．

木津川の流出解析

6．1

かつて行なった，加茂地点の流出解析

1）木津川加茂は，菅原・勝山が昭和31年に流出解析を行なっている．それが直列貯留I型のタソクモデル（3×1型）を用いた最初であった．

木津川の加茂は河床低下が大きい所で，しかも洪水が出ると一時的に河床は上がる．その結果，流量の信頼性はあまり高いといえない．そこで，雨量から推定流量を出し，その平方根にある定数を掛けて推定水深を出し，それと実測水位とを比較することにより，流出機構と，加茂の河床高の長期的変化との両方を求めようとした．

それには，水位流量曲線
Q＝α（ゐ一乃・）2

において，河床変動により河床高伽は変動するが，αは変わらないと仮定する．この仮定
一21一

国立防災科学技術セソター研究報告

表3

木津川農業用水量
農業用水（m3／・）

常

時

最

大

αと，昭和24年〜29年の資料から菅原が定めた

（ha）

木津川

88．6

104．8

12，718

長剛11

33．4

36，8

3，609

名張川

38．3

47，5

ユ975年8月は大正6年〜9年の資料から坂本助太郎氏の定めた

水田面積

■

第12号

αとが，あまり相違しないことから，ほぼ妥当であると考えた．

2）河床高の長期的変動と，流出機構とを同時に求めようとするのが，無理な試みであることが，やがて明らかになった．

最初の木津川の流出解析では，農業用水取水の影響は無視した．当時，木津川の農業水利権を調べた結果が表3である．加茂の低渇水流量が15m3／s程度であることを考えると，これだけの水が取水できるはずがない．もし取水できるとすれば，取水した水がすぐ遠元し，

くり返し利用されているからであろう．すぐ還元するとすれば，農業用水の影響を無視しても大過あるまいと考えた．

農業用水を無視して計算すると，6月ごろには明らかに推定値が過大に出るが，これは農業用水を無視したからで，それは致し方ない．しかし，その後は取水と還元とが消し合って，大きな相違は出ないだろうと楽観していた．その後3年ばかりして，別の欠点が分かって来た．

養老川の日流量解析では，農業用水を考慮に入れなければ，実測と合わなかった．農業用

水を計算に入れると，農業用水取水期に河の水が滅るだけでなく，秋から冬にかけて，還元
水により流量がふえる．それを見て，木津川加茂の解析結果の欠陥に思い当たった．

加茂の水位解析の際，秋になると河床が上昇する傾向が見られた．それは洪水による河床の上昇であると解釈したが，実は農業用水還元による基底流量増大によるものであろう．流
出解析と河床高変動とを同時に行なう立場からでは，秋の基底流量増加と，河床上昇とを区別することは不可能である．

3）農業用水の影響を計算に入れて，基底流量の変化を調べるためには，河床変動が小さ
く，渇水時流量の信頼性の高い地点の資料を用いる必要がある．そこで木津／11加茂の計算をやり直す代わりに，名張川月ケ瀬の流量で流出解析を行なうことにした．昭和35年のこと
である．

6．2

名張川月ケ瀬の流出解析

1）養老川の流出解析で行なった農業用水の計算は簡単である．雨をタンクモデルに投入
して流出高に換算する．それから農業用水を差し引いた残りが推定流出高になる．農業用水は浸透すると考え，4段のタンクモデルの3段目に戻す．
この方式は，多くの河川で，かなりよい結果を与えている．問題は農業用水量の決定で，

それは農業水利権として認められている取水量よりかなり小さいものであるのがふつうである．農業用水として坂水された水の一部には，用水路からそのまま排水路を通って河に戻る
＿

フリ

＿

タ：／クモデルによる非湿潤地帯河川の流出解析一菅原・尾崎一渡辺・勝山
表4減水深（mm／日）

ものもあろうし，田面を通過するのみで排水路に戻るものもあ

水田面積比率ろうし，浸透しても直ちに河に戻るものもあろう．かかる短期（％）

還元水は無視しても，流出解析の結果にはあまり影響しない．

0〜10

20

10〜50

40

50〜10030

100〜200

10

閉題になるのは，長期還元成分であるが，どれだけが長期であるかを知るには，試行錯誤法で探し求める以外に方法がない．

2）名張川では，農業用水の影響を上の簡単なモデルで計算できないようである．月ケ瀬の流域面積は615km2であるから，表3の38．3m3／sの水が常時取水できないのは当然であろうが，渇水時にもかなりの流量が残っている所をみると，渇水時はわずかな坂水で我慢しているらしい．かなりの流量が河に残るのは，取水地点が上流にあり，そこから下流にある地域からの流出は，農業用水坂

水の対象にならないからであろう．渇水時に坂水を我慢している代わりに，雨が降ったとき，河にあまり水が出て来ないことがある．増水を水田に坂り入れるからであろう．したがって，水田に水がなく，川に水があるとき，すなわち取水の必要があり，それが可能なときは，かなり大量の水を坂水しているらしい．その意味で表3の取水量は現実的なのであろう．

農業用水の影響の出方があまりに不思議なので，現地に行き，名張市役所で聞いた，水田の減水深調査の結果が表4である．

3）以上の事情を考え，次の方式で農業用水の計算をすることにした．
a）水田面積は流域の2，5％とする．

b）o．4mm／日の農業用水は常時取水する．これは水田に対し16mm／日の供給である．

表4の滅水深の平均は約50mm／日であるから，仮に水田の滅水深はすべて50mm／日とする．上記16mm／日の供給の目が続けば，水田に100mmの水がはってあっても，やがて水は無くなる．しかし，16mm／目の水が供給される限り稲は枯れない．

O．23◎c）河に水があれば，水田に水深100mmの水をはるまでは取水する．ただし坂水地点の集水面積は流域の50％と考え，河川流量の1／2までしか取

N 012
n ．30

水できない．

以上の規則には不白然な所もあるが，当時リレー計算機で計算していて，1年分の計算に1時間程度かかったと記憶している．したがって，規則が簡単で，手間がかからないことが必要であった．

34）図17はこのとき得られた名張川月ケ瀬の流量に対する流出機構であ

図17名張川

る．

6．3

月ケ瀬に
対するタ
ソクモデ

木津川加茂の流出モデル

1）今回，木津川加茂（流域面積1，456km2）については，農業用水のモ

ー23一

ノレ

国立防災科学技術セソター研究報告第12号

1975年8月

デルを次のように考える．

表5a）表3に従い，水田面積を127km2とする．表4を参考に

減水深（mm／日）

水田面積比率

ユ0

20

30

40

mm／目，120mm／日の水田が，表5の比率で存在すると仮定

75

30

する．

120

10

（％）

し，加茂全流域に対して，減水深10mm／日，30mm／目，75

この仮定に従えば，これら水田の減水深に対して必要な補給
水量は（10×O．2＋30×0．4＋75×O．3＋120×O．1）×10−3x127×106x（1／86，400）＝71．3（m3／s）

である．これは表3の88．6m3／sと大体合っているとみてよい．

b）以下の計算方式を図18を見ながら理解していただきたい．
タソクモデルからの流出高μをんμと（1一ん）ひとに分ける．切が農業用水の取水対象となる流量，すなわち取水地点より上流の地域からの流出である．切から，表5の4種類の水田にそれぞれ21，幼，2・，忽の水が供給され，余りは（1一ん）ひに加えられ，推定流出高σとなる．ここではん＝1／2とした．

ヲ■o・k）y」」一

止＼・。→〜一Σ・己

L．L

1マ＼＼
子 テテ㍗

団団出出
↓

W1

↓

↓

↓

W2

W3

W4

図18農業用水の影響の計算の流れ

一24一

■

タソクモデルによる非湿潤地帯河川の流出解析一菅原・尾崎・渡辺・勝山

c）水田には取水できるならば100mmの水をはる．なお減水深120mm
O．23

No O．12nの水田には，計算技術の便宜から120mmの水をはる．この限度以上には取水しない．

d）21，…，2・は次のようにして定める．まず，一様にんμを配分する．

すなわち4種類の水田に面積に比例して配分する．ここで，減水深10mmの水田の水深が100mmに達すれば，余った水は残りの水田に一様に配分する．この段階で，減水深30mmの水田の水深が100mmに達すれば，余った水を残りの水田に一様に配分する．以下同様で，すべての水田が限度に達すれば，余った水は河に戻す．

唐1
図19

現実の水配分がこのようであるとは思わないが，水をなるべくむだなく配

木津川
加茂に対
するタソ
クモデル

分しようとすれば，このようなものになるであろう．そして，現実の姿も，案外これに近いかもしれない．

e）

4種類の水日ヨからは，それぞれの滅水深に応じて水が浸透する．ただし，水田に水が無くなれば，それ以上は浸透しない．

f）

4種類の水田からの浸透の和が，タソクモデルの3段目に戻される．

2）以上の農業用水の計算方式と，図17の名張川の流出モデルとを用いて加茂の流出高を算出し，実測と比べると，基底流量がいくらか大きすぎるようである．加茂地点では伏流がいくらかあると思われる．そこで図17のモデルの4段目のタソクの底に浸透孔をつけ，基底流量の一部（約1／4）を伏流にまわすことにし，図19のモデルを作った．これでかな
り合うようである．

3）図16は得られた結果の一部を示す．実測との一致は必ずしもよくないが，加茂では河床変動が大きいことを考えると，あまりよい一致は望めない気がする．

7．謝辞

本研究の資料入手につき，WMOのJ1Nさmec教授，アメリカエ兵隊のD．Rockwood
氏，建設省河川局河川計画課，近畿地方建設局淀川ダム統合管理事務所，当セソター大型実験研究部長木下武雄氏に，お世話になった．ここに厚く感謝する．

参1）Sugawara，M。（1961）l

On

the

analysis

考
of

文

献

runo鉦structureabout

severa1Japaneserivers．Jψ、∫G60ψ1リ∫．，2，No．4．

2）水資源の循環機構に関する調査報告，科学技術庁資源調査会報告，第47号，昭和44年（1969）．
3）菅原正巳（1972）：流出解析法1共立出版株式会社．
4）Lins1ey，R．K．and

digitalcomputer．

Crawford，N．H．（！960）：

ComputationofasyntheticstreamHow

recordon

aP肋Z．」Mo．51，1λ8H．

5）菅原正已，勝山ヨシ子（1964）：半年間の乾季を持つ地域の河川の流出について，資源調査会資料．
（上記2）に収録）

一25一

国立防災科学技術セソター研究轍告6）Bumash，R．J．C．，Ferral，R．L．and
system；conceptualmodelingfor

McGure，R．A．（1973）：A

digitaltoBirdCreek，Wo1lombi

1975年8月
genera1ized

streamHow

simulation

computers，7）Sugawara，M．，Ozaki，E．，Watanabe，I．and
cations

第12号

Katsuyama，Y一（1974）：Tank

Brook，Bikin

River，Kitsu

mode1andRiver，Sanaga

River

and

its
Nam

appli−
Mune．

Rω．八b姑ハ肋1．地∫、C6〃．一α∫α∫肋・P1洲．，No．11，1−64−

8）菅原正巳（1973）：秋の雨はおもに平地に降る．水利科学，No・92，17（3）・（1974年5月1日原稿受理）

一26一