3

水物質循環モデルの構築

3.1
3.1.1

水物質循環モデルの概要
水物質循環とは

ここで、水物質循環とは、地表における巨視的な水移動のシステムを意味する「水循環」と、この水循環に伴って運ばれる物質の流れのシステムを意味する「物質循環」を合わせたものである。
水循環は、降雨、貯留（窪地、湿地、湖沼等）、表面流出（河川等）、地下浸透（地下水の流れ等）、蒸発散といった自然系の水の循環と人工的に構築された上水道（河川からの導水、地下水の汲み上げ等）や下水道による雨水排水、汚水排水といった人工系の水の循環から構成されている。例えば、自然系と人工系の水循環系をそれぞれ模式化すると図3.1.1.1、図- 3.1.1.2のようになり、水が存在形態と存在場所を変えながら循環していることが分かる。また、この水循環に伴う汚濁負荷等の物質移動は、取水・排水等の人為的な現象と希釈、拡散、移流、沈殿、吸着・脱着、溶解・分解、浸透、掃流、生物による取り込み・排泄等の自然現象から構成されており、物質の移動過程で生物化学的な変化も生じている。
この水物質循環は、河川、湖沼等における水質、水量の問題を考えるに当たっての最も基本的かつ普遍的な工学的システムと言える。

図- 3.1.1.1

自然系の水循環系のイメージ 1)

33

図- 3.1.1.2
3.1.2

人工系の水循環系のイメージ 1)

水物質循環モデルとは流域における水物質循環は、流域の状態を捉え、流域の健全性を診断し、流域再生のあり方を探るための最も基本的な現象である。水物質循環モデルには「水循環モデル」と「物質循環モデル」に分類される。水量を推定する水循環モデルには、表- 3.1.2.1 に示すように年間総量ベースの大まかな四則演算により静的な水収支を定量化する簡便な方法（年間水収支法）から、経路毎の水量の時系列変動を推定できる精度の高いものまで種々の方法がある。
水物質循環におけるシミュレーション・モデルとして、物理的な基礎式に基づいているか否かで、物理モデルと概念モデルに分けられる。概念モデルとしては単位図、タンクモ
デル、貯留関数法等が挙げられる。概念モデルは複雑な事象を単純化していることから、
モデルパラメータの同定は、経験的であったり、試行錯誤的となる。一方、物理モデルは、各種提案されてきており、モデルパラメータは、水循環については基本的には対象流域の物性値で与えることができる。このことは、水循環系の将来予測や施策の効果の評価において、物理モデルのシミュレーション結果に一定の信頼性をおくことができる事を意味する。しかしながら、多くの要素をモデルに組み込み、複雑になるほどパラメータの数が膨大になることから、表- 3.1.2.1 で示されているように、物理モデル構築に必要なデータや費用は莫大なものとなり、実務上の課題となっている。

34

表- 3.1.2.1

水循環の定量化手法の特徴 1)

分 類

年間水収支法
概念モデル
物理モデル
降雨量、流量、人工系給 流出の各プロセスをタンク 流出の各プロセスを数理モ
排水量などの観測結果と モデル、貯留関数などの概 デルで表現したもので、計
地目別流出率などの経済 念的なモデルで表現したも 算上の流向の仮定により１的な係数により経路毎の循 の。モデル中のパラメータ 次元モデルから３次元モデ
環量を概略推定するもの。 は観測流量と一致するよう ルに分かれる。モデル中の計算は貯留量の時間変化 に設定する。
定数は計測可能な物性値を無視するため四則演算で 年間水収支法、物理モデ で設定することを前提とす実行される。
ルの中間的な特徴を持つ。 る。
計算は極めて容易である
精度の高い予測が期待さが、計算結果は年間値のれるが、演算に要するトー
目安程度の利用に留まる。
タルコストが高い。
入力情報の多少少中多
出力情報の多少
少中
多演算の難易度低中高結果の厳密性
低中
高計算時間単位
年
日
秒～時間
流域一括
計算格子点毎
計算空間単位流域一括あるいは小分割流域毎
入 力 条 件 （計算空間単位毎）
降雨量
年間値を入力計算単位で入力計算単位で入力
蒸発散量
実蒸発散量を入力
実蒸発散量を入力
可能実蒸発散量を入力
浸透・不浸透の別、
土地利用情報浸透・不浸透の別を入力 浸透・不浸透の別を入力
等価精度係数を入力
人工系給排水量
年間値を入力計算単位で入力計算単位で入力
タンク孔定数、貯留量関数 透水係数、不飽和帯特性
自然系保水特性
考慮せず
等を入力
を入力
可動特性
考慮せず
考慮せず
形状、精度係数を入力
算定項目
経路毎の年循環量
算定可
算定可
算定可
日流量・流況
算定不可
算定可
算定可
高水流量
算定不可
一部のモデルで算定可
算定可
地下水流去量算定不可
一部のモデルで算定可
算定可
地下水位・かん養量
算定不可
一部のモデルで算定可
算定可
概 要注）入力条件、算定項目は一般的な構成、内容を掲げたものであり、個別のモデルにより多少変化する。

水質を定量化するためのモデルには、流域全体を対象とした発生汚濁負荷量を算定し、水量との関係からこれを水質に変換するモデル、また、湖沼や河川における一定の領域を対象に有機物の内部生産を考慮した水質シミュレーションモデル等がある。前者は河川や湖沼に流出する汚濁物質の発生源を流域全体で把握し、水域の将来水質を予測する際に有効である。後者は河川や湖沼に排出された汚濁の集中や拡散等の状態を一定の領域内で把握する場合に有効である。
なお、水物質循環モデルの概論については、「都市小流域における雨水浸透、流出機構の定量的解明」研究会から出版されている「都市域における水循環系の定量化手法」に詳しく解説されているので参考にされたい 2)。

35

3.1.3

水物質循環モデルの開発の意義

流域水物質循環モデルの検討対象としている東京湾とのその流域では、流域の高度な土地利用と閉鎖性水域における沿岸域の開発などにより水環境など悪化を招いている。これに対応するため様々な政策や事業が実施されてきているが、閉鎖性水域を抱える流域については未だ根本的な解決に至っていないケースが少なくない。このため、土地利用等流域
条件や取排水条件の異なる様々な流域の水物質循環現象の再現が可能な分布型水物質循環モデルである「陸域モデル」と、流域からの外部負荷、底泥溶出等の内部負荷及び植物プ
ランクトンの内部生産による水域の水質形成機構を表現できる「水域モデル」の統合モデルを構築することにより干潟などの自然再生や下水道整備などのインフラ整備をはじめとする様々な環境改善施策に対する検討を行うことが可能となる。また既存の知見の範囲内
で過去の人口条件、土地利用状況、気象条件などを設定することにより、過去から現在に至る水環境の変遷を理解することができる。
3.1.4

水物質循環モデルの開発目標

以下の点に留意して水物質循環モデルを構築する。
①

水環境の改善に係わる政策検討を目的とした水物質循環モデル
水問題の中でも特に平常時の水質・水量に係わる水環境問題の解決に資する「水物質循環モデル」を構築する。このため、流域の流出については、日平均流量のレベルで高水及び低水の再現を図るとともに、物質に関しては COD、N、P の汚濁負荷を対象の基本とする。
②

汎用型の水物質循環モデル
検討対象とする東京湾流域のみならず、全国各地域・水域の水物質循環現象の再現や水環境問題への適用が可能な汎用型モデルとする。このため、土地利用等流域条件や取排水条件の異なる様々な流域の水物質循環現象の再現が可能な分布型水物質循環モデルである「陸域モデル」と、流域からの外部負荷、底泥溶出等の内部負荷及び植物プランクトンの内部生産による水域の水質形成機構を表現できる「水域モデル」を構築する。
③

操作性の高い水物質循環モデル
開発する水物質循環モデルについては、今後普及し、全国の水環境問題の解決に資することを開発目的としており、このためには、シミュレーション・モデルについて特別の知識を有しない現場管理者のレベルで、容易に使いこなせることが重要と考えている。このため、図- 3.1.4.1に示すように GIS をベースとしたユーザーインターフェイスを整備し、土地利用の改変等流域条件の変化や水環境に係わる様々な施策の入力や計算結果の出力を容易にする。
④

理論構造や計算内容が分かりやすい水物質循環モデル
開発対象とする水物質循環モデルについては、当該流域の水物質循環現象を限りなく精緻に再現するというよりも、関係者による現象や課題の共通理解、合意形成等を支援することを目的としており、政策議論に必要な精度・信頼性は確保しつつも、モデルの構造や計算の内容が分かりやすいモデルとすることとし、水物質循環の基本的現象に対応するモ36

デルとする。なお、適用する流域や現象の特殊性により、本モデルの原理だけで対応できない場合に対して、その特殊性に対応する個別的・部分的モデルを作成し、本モデルに連結することによる適用性を確保する。

従来の手順

細かな設定・変換作業が必要、扱える人間が限られるモデルに合わせた
データの整形

各種のデータ

施策シナリオに応じたパラメータの設定

a

表

b

c

表示ソフトへの計算結果の変換

シミュレーション

結果の表示

・・・

表・グラフ

手作業で設定Excel
テキスト
GISデータ

Excel

Excel

モデル

入力データ

テキスト
計算結果
Excel

地図情報

GISデータ

データ変換

データ変換

ﾕｰｻﾞｰｲﾝﾀｰﾌｪｲｽ
各種の
データ

地図表示

作業の自動化・簡略化、幅広い普及が可能
施策シナリオの選択
出力形式の指定シミュレーション

結果の表示

シナリオ

表

表・グラフ

メニューから選択

Excel

Excel
GISデータモデル

データの自動変換

指定した形式で表示
GISデータ

地図情報

地図表示

図- 3.1.4.1
3.1.5

ユーザーインターフェイスの意義

本研究における水物質循環モデルの特徴

本モデルの特徴としては、以下のことが挙げられる。
①

分布型モデル：流域をメッシュに分割し、メッシュ単位で解析を行うことから、任意地点（メッシュ）での水・物質の挙動の解析・出力が可能②

非定常モデル：平常時、降雨時を含めた任意時間での解析が可能③

水量・物質連成解析モデル：水の挙動および物質の挙動の同時解析が可能④

統合型モデル：流域内の諸現象（蒸発散・地下水・地表流・河道流・水域）に応じて、各モデルにより解析し、それらを連成することによって流域全体の現象を解析⑤

GIS 連成型モデル：GIS を介しての流域データのモデルへのインプットデータ作成、解析結果の格納、可視化⑥

施策指向型モデル：雨水浸透桝設置や生活排水対策等、様々な水循環健全化施策による改善効果を解析することが可能

37

参考文献
1) 健全な水循環系構築に関する関係省庁連絡会議：健全な水循環系構築のための計画づくりに向けて，pp37-45，2003.
http://www.mlit.go.jp/tochimizushigen/mizsei/junkan/keikakudukuri.html
2) 「都市小流域における雨水浸透、流出機構の定量的解明」研究会：都市域における水循環系の定量化手法-水循環系の再生に向けて-，pp.44-82，，2000.

38

3.2

水物質循環モデルの開発

本研究で構築した水物質循環モデルは流域モデル及び水域モデル（湾モデル）から構成される。以下、3.2.1 では流域モデルの概要を 3.2.2 では水域モデルの概要を説明する。
3.2.1

流域モデルの概要

流域モデルは、流域全体の水物質循環を表現する“流域水循環モデル”と、流域内における汚濁物質流動を表現する“流域物質循環モデル” により構成される。
“流域水循環モデル”は、図-3.2.1.1 に示すように、降雨、蒸発散等の自然系水循環を領域とするモデルと、農業用水・上水等の人工系水循環を領域とするモデルとからなり、両者をあわせることにより流域内の水循環を表現するものである。
また、“流域物質循環モデル”は、降雨や大気由来、及び人為的に地表面に散布される等の汚濁源を特定されない面源負荷（ノンポイントソース）を対象とするモデルと、下水処理場などの生活系や工業系由来などの汚濁源の特定される点源負荷（ポイントソース）を対象とするモデルとからなり、両者にあわせることにより流域内における物質流動を表現するものである。ここで、これら汚濁源の流出は、上述の“流域水循環モデル”により導出される“水の流れ”に規定され、地表面に堆積する物質は掃流され、また水域内に建濁している物質は流達することになる。
以上より、両モデルの関係を示せば「流域水物質循環モデル＝流域水循環モデル＋流域
物質循環モデル」と表すことができる。実流域に流域水物質循環モデルを適用する場合は、①水循環ならびに汚濁物質流動を表現するためのモデルを構築し、②流域水循環モデルによって流域内水理量を解き、③流域内水理量をもとに汚濁物質流動を解く、という手続を踏むことになる。

図- 3.2.1.1

流域水循環モデルの定量化項目のイメージ 1)39本検討での流域モデルは、タンクモデルを対象としその概念図を図- 3.2.1.2、図- 3.2.1.3に示す。

自然系水循環 降雨

Sf0Sf1Sf2

地下水層
モデル

農業用水

表面流
5/3 1/2= Kf2 /N･(h-Sf2) ･i
早い中間流=α0･αri･Kf0･(h-Sf1)
地下浸透 =Kf0･α0･(h-Sf0)

表層
モデル
不飽和層
モデル

人工系水循環灌漑取水の供給（表層）

蒸発散

D

灌漑用水の還元 灌漑

生活用水
工業用水
排水

上水道からの漏水（不飽和層）

地中浸透 =Ks0･Kz
不圧地下水流=kg12･Au2･(h-Sg)2
被圧地下水流=kg0･Ag･h

Sg

取水

河道
モデル

遅い中間流
=Ks1･kx･D･i

下水管路への浸出
地下水取水

N: 等価粗度係数
i : 斜面勾配
α0 : 表層モデルの浸透能
αri :早い中間流の定数i: 斜面勾配
Kf2,Kf0,Ks1,Ks0 :補正係数
Sf2,Sf1,Sf0,Sg: 表面流，早い中間流，地下浸透，不圧地下水流の発生高さ
kx,kz :不飽和層透水係数
D : 層厚Au: 不圧地下水減水定数
Ag :被圧地下水減水係数

図- 3.2.1.2

自然系物質循環
表層
モデル

脱着 Cf
吸着
Snp 掃流Hf

人工系物質循環農業系負荷

表層流出負荷mn

Lf２= knp・Snp ・A・Qf

灌漑用水の還元
取水負荷

畜産系負荷

浸透負荷 Lf1 = Lf2/Qf0・Qf1

不飽和層
モデル
Hsmx 脱着 Cs
吸着
HS
Ss

不飽和層流出負荷
Ls2 = ks・Ssm・A・Qs2n

浸透負荷 Ls1 = Lss・Qs2/Qs1
地下水層
脱着 Cg
モデル
Hgmx
Hg
Sg 吸着

湾モデル

流域水循環モデルの概念図

灌漑取水の供給(表層)降雨

取水

地下水層流出負荷
Lg = kg・Sg・A・Qgn

河道
モデル

取水

生活系負荷
工業系負荷

排水
上水道からの漏水(不飽和層)
下水管路への浸出
地下水取水

Ci : i 層の貯留水濃度
Hi : i 層の貯水位（i = f,s,g）
Qi2 : i 層からの流出量
Qi1 : i 層からの浸透量
Li1 : i 層の浸透負荷量
Li2 : i 層の流出負荷量
湾モデル
Snp : 面源負荷由来の堆積物量
Hsmx,Hgmx : 最大貯水位
Ss,Sg : 土壌内の汚濁負荷蓄積量
knp : 面源負荷由来の堆積物の掃流係数

図- 3.2.1.3

流域物質循環モデルの概念図
40

(1) 流域水循環モデルの概要
分布型の自然系水循環モデルは流域をメッシュ分割して、メッシュ毎で降雨－流出－地下浸透などの計算を行い、表面流・地中流をメッシュ間の移動として計算するものである。
図-3.2.1.4 に、自然系水循環モデルのメッシュ単位での計算概念図を示し、各タンクでの解析方法を次ページ以降に説明する。

降雨
降雨

蒸発

表層

落水線
不飽和帯

粘土層

被圧層

河道

図- 3.2.1.4

自然系水循環モデルのメッシュ単位での計算概念図

41

1) 表層モデル
表面流

表層モデルの構造は右図のとK1f2・(1/n)・(h－Sf1)5/3・i1/2

L (h − S f2 )5/3 i
Nおりである。雨量を表面流、早い中間流、地下浸透流に分離すSf2るモデルである。上から１番目

h

早い中間流
Sf0

Sf1

の穴が表面流の穴で、２番目の地下浸透

穴が早い中間流の穴で、下の穴

Kf0−・α0・(h－Sf0)
A･f0(hS fo )/(Sf2 − S f0 )が地下浸透の穴である。

図- 3.2.1.5

Sf2：表面流の発生する高さ(m)α− S f1)/(S f2 −f1)S f1 )
Knf0･A･f・α0(hri・α0・(h－S

表層モデルの概念図

Sf1：早い中間流の発生する高さ(m)Sf0：地下浸透の発生する高さ(m)h：タンク内水位(m)L：メッシュ長(m)N：等価粗度係数(m-1/3・s)
i：斜面勾配（落水線方向の勾配：メッシュ内は一様に傾いていると仮定）
α ri：早い中間流の定数（地下浸透量に対する比率）
α 0：表面タンクの浸透能(cm/s/m)
Kf2, Kf0：補正係数① 表面流
表面流はManning則が成り立つものとして基礎式を決定している。
・ Manning平均流速公式
V＝(1/n)R2/3i1/2

V：平均流速(m/s)、n：粗度係数(m-1/3・s)、R：径深（m）、i：

勾配（流れは等流と仮定し勾配を水路の河床勾配と同一とする）
上式を表面タンクの構造に当てはめ、他のタンクからの流出量とのバランスをとるための補正係数を考慮して、
qf2＝Kf2×(1/N)×(h－Sf2) 5/3×i1/2
Kf2：補正係数（＝1/L（L：メッシュ長(m)）で取り扱っている）となる。
② 早い中間流出量、地下浸透量
早い中間流出量、地下浸透量は浸透能に比例するものとする。
地下浸透量は浸透能から決定されるものとし、算定式を以下としている。
qf0＝Kf0・f
qf0：地下浸透量高(m/s)、Kf0：補正係数、f：浸透能(m/s)
ここで、表面タンクの浸透能はタンク水位が表面流を発生するときの高さh＝Sf2の時の値が最終浸透能と等しくなり、h＝Sf0の時ゼロとなることとした。
（浸透能は降雨強度により異なり、降雨開始直後にきわめて大きな値を示し、急激に減
42少し、最終的に一定値（最終浸透能）に近づくように変化する。最終浸透能は土壌の性質
により異なり、降雨強度が浸透能を超えると地表流が発生することとなる。）
従って、浸透能fはタンクの水深（h－Sf0)に比例して最終浸透能f0を配分して与えるものとした。
f＝α0×(h－Sf0)：α0＝f0／(Sf2－Sf0)

f0：最終浸透能(m/s)

従って、地下浸透量は以下の式で表現される。
qf0＝Kf0×α 0×(h－Sf0)

一方、早い中間流出量は浸透量の一部が地下に浸透する能力を上回る強度で供給された
場合に表層タンクの底部から横方向に移動する量として与えられるものとし、h－Sf1 に比
例する値により得られるものとして算定する。その比例定数をαriとして、以下のように記述される。
qf1＝α ri×Kf0×α 0×(h－Sf1)
表面タンクの連続式は以下のとおりである。
①

ｈ≧Ｓf2の場合∂h= R − E ps − q f0t − q f1t − q f2t∂t②（式- 3.2.1.1）

Ｓ f1≦ｈ＜Ｓ f2の場合∂h= R − E ps /S f2･h − q f0t − q f1t∂t③（式- 3.2.1.2）

Ｓ f0≦ｈ＜Ｓ f1の場合∂h= R − E ps /S f2･h − q f0t∂t④（式- 3.2.1.3）

ｈ＜Ｓf0の場合∂h= R − E ps /S f2･h∂t

（式- 3.2.1.4）
qf0t：地下浸透量高(m)

ここに、Ｒ ：雨量(m)Ｅps：蒸発量(m) : 5)で詳述

qf2t：表面流出量高(m)Ｓ f2：表面流の発生する高さ(m)qf1t：早い中間流出量高(m)Ｓ f1：早い中間流の発生する高さ(m) ｈ：タンク内水位(m)Ｓ f0：地下浸透の発生する高さ(m)432) 不飽和帯モデル
不飽和帯モデルは不飽和帯の水分移動を考え たモデルで ある(タン ク
内の水分を層厚に対する水分量(比)

遅い中間流
SS2

D

で取り扱う)。モデルの構造は右図のK出
s1 ・Kx ・D・ih

とおりである。遅い中間流は不飽和SS1

透水係数に比例し動水勾配が近似的に地形勾配に等しいとして流出量を地中浸透計算する。地中浸透量は不飽和透水係数に比例し、動水勾配を1として計

Ks0・KZ図- 3.2.1.6

不飽和帯モデルの概念図算する。
（透水係数について：透水係数は多孔質体中の水の流速の大きさを示す指標で、飽和時の透水係数を飽和透水係数、不飽和時は不飽和透水係数といい、不飽和透水係数は飽和透水係数に較べかなり小さい値となる。また、地下水の流れはダルシーの法則により得られるとし、これはフランスの水道技師Henry Darcy(1856)が、砂ろ過の実験から明らかにした地下水流動に関する法則である。地下水の流速は透水係数と動水勾配の積で表わされるとするものである。）
SS2：飽和水分量相当高
SS1：最小水分量相当高
不飽和モデルは水分量で水の移動を考えるため、
h：タンク内水分量相当高
これらは概念的な諸量である。
D：層厚(m)Kx：水平方向の不飽和透水係数(m/s)
Kz：鉛直方向の不飽和透水係数(m/s)
i：斜面勾配（落水線方向の勾配：メッシュ内は一様に傾いていると仮定）
Ks1, Ks0：補正係数① 遅い中間流出
遅い中間流はDarcy則に従うものとし、タンク内の水分量により透水係数が変化するものとする。つまり、不飽和透水係数を水分量に応じて推定して算定することとする。
・ Darcy則
q＝k×i

q：流出高（m/s）、k：透水係数 (m/s)、i：動水勾配

動水勾配が近似的に斜面勾配に等しいと仮定し、上式を不飽和タンクの構造に当てはめ、補正係数を考慮し、以下のとおりとなる。
qs1＝Ks1×Kx×D×i

qs1：流出高(m/s)、Ks1：補正係数（＝1/Lで取り扱っている）、Kx：水平方向の不飽和透水係数(m/s)、D：層厚(m)、i：動水勾配

ここで、不飽和透水係数は水分量に比例するものとし、ｈ＝Ｓ s2で飽和透水係数と等しくなり、ｈ＝Ｓ s1で０となるよう、次式により推定する。
44

Kx = Ksx / 100 ⋅exp(b ⋅ θ ) − exp(b ⋅ θ w )exp(b ⋅ θ s ) − exp(b ⋅ θ w )Kx：水平方向の不飽和透水係数(m/s)，
Ksx：水平方向の飽和透水係数（cm/s,定数として与える），
b：定数，θ：水分量（h /D），θs：飽和水分量（SS2 /D）（定数として与える），
θw：最小水分量（SS1 /D）（定数として与える），
SS2：飽和水分量に相当する貯留高，SS1：最小水分量に相当する貯留高② 地中浸透地中浸透もDarcy則に従うものとし、タンク内の水分量により透水係数が変化するものとする。
・ Darcy則
q＝k×i

q：流出高(m/s)、k：鉛直方向不飽和透水係数(m/s)、i：動水勾配

上式を不飽和タンクの構造に当てはめ、補正係数を考慮すると、以下のとおりとなる。
qs0＝Ks0×kz×i
ここで、不飽和透水係数は水分量に比例するとして、次式のとおりとする。

kz = ksz ⋅exp(b ⋅ θ ) − exp(b ⋅ θ w )exp(b ⋅ θ s ) − exp(b ⋅ θ w )kz：垂直方向の不飽和透水係数(m/s)，
ksz：垂直方向の飽和透水係数（m/s,定数として与える），
b：定数，θ：水分量（h /D），θs：飽和水分量（SS2 /D）（定数として与える），
θ w：最小水分量（SS1 /D）（定数として与える），
SS2：飽和水分量に相当する貯留高，SS1：最小水分量に相当する貯留高※

不飽和透水係数は土壌内の水分量との関係は基本的に実験や回帰によって得られるものであるが、計算を簡便化するために用いた関係式である。水理公式集ではBrooks
and Coreyによって提案された下記の不飽和透水係数と飽和透水係数との関係式を紹介している。

⎛ θ − θw ⎞⎟⎟kz = ksz ⋅ ⎜⎜
−
θ
θ
w ⎠
⎝ s

3+ 2 / λ

λ：土壌の特性による経験的パラメータ・0.1～10程度の範囲で変化する。
※

次ページに定数bを変化させたときの飽和透水係数にかかる比率とBrooks and Corey
式のパラメータλを変化させたときの飽和透水係数にかかる比率との比較図を示した。これより本モデルの定数bを15とした場合の比率がBrooks and Corey式のパラメー
タλの一般値上限である10以下の場合と整合がよいため、定数b=15とした。

45

図- 3.2.1.7

B=1
B=15
λ=0.1
λ=5

定数による透水係数比率の違い
1.0

B=10
B=100
λ=2
λ=10

0.8
0.6
0.4
0.2

0.58

0.54

0.46

0.42

0.38

0.34

0.0
0.30

1
10
15
100
0.30
0.000
0.000
0.000
0.000
0.32
0.058
0.012
0.004
0.000
0.34
0.117
0.026
0.009
0.000
0.36
0.177
0.043
0.016
0.000
0.38
0.238
0.064
0.026
0.000
0.40
0.301
0.090
0.039
0.000
0.42
0.364
0.122
0.057
0.000
0.44
0.430
0.160
0.081
0.000
0.46
0.496
0.207
0.113
0.000
0.48
0.564
0.265
0.156
0.000
0.50
0.633
0.335
0.214
0.000
0.52
0.703
0.420
0.293
0.000
0.54
0.775
0.525
0.400
0.002
0.56
0.849
0.653
0.544
0.018
0.58
0.924
0.809
0.738
0.135
0.60
1.000
1.000
1.000
1.000
不飽和透水係数算出のさいの飽和透水係数にかかる比率

0.50

0.6
0.3

飽和透水係数にかかる比率

飽和水分量θs
最小水分量θw
定数b
タンク内水分量θ

タンク内水分量

不飽和帯モデルの係数と透水係数の関係

不飽和タンク層から地中への浸透は不飽和層での地下水流動の鉛直方向移動量として把握されるべき量として見積もられ、ここでは近似的に動水勾配を1.0とし、不飽和透水係数に比例して地下浸透量が発生するものとして取り扱う。
qs0＝Ks0×Kz
不飽和タンクの連続式は以下のとおりである。
①

ｈ≧Ｓ s2の場合飽和状態であるため上段の表面タンクから不飽和タンクへは浸透しないと考える。両タンク間の動的な結合は計算を複雑にすることから、不飽和タンク水分量相当高がＳ s2を越える場合は表面タンクからの浸透量を落水線下流の表面タンクに流入させる。
②

Ｓ s2＞ｈ≧Ｓ s1の場合

θs
③

∂h
= q int − E ps − q s0t − q s1t∂t（式- 3.2.1.5）

ｈ＜Ｓ s1の場合

θsδh= qin − E ps /Ss1･hδt（式- 3.2.1.6）※タンク内の水分量が最小水分量より小さいときは遅い中間流も地中浸透も発生しない

qs0t

ここに、Ｓs2：飽和状態の貯留高(m)：地中浸透量高(m)Ｓs1：最小水分量に相当する貯留高(m)qs1t ：遅い中間流出量高(m)qint：表面モデルからの流入量高(m)θs

Ｅps：蒸発量(m)ｈ ：不飽和モデルの貯留高(m)46：飽和水分量3) 地下水モデル
地下水モデルの構造は右図のとおり

不圧地下水流不圧地下水流A 2u･(h − S2g )2・A
Kg1・Au ・(h-Sg)2である。横の穴の上から１番目は不圧
地下水の穴で、２番目が被圧地下水のh Sg

彼圧地下水流A被圧地下水流g･h･A
Kg0・Ag・h穴である。地下水の流出は被圧地下水がhに比例し、不圧地下水が（h－Sg）2
図- 3.2.1.8に比例すると考える。

地下水モデルの概念図

h：タンク水位(m)Sg：不圧地下水発生高(m)Au：不圧地下水減水定数(mm-1/2・day-1/2)
Ag：被圧地下水減水定数(1/day)
Kg1, Kg0：補正係数

不圧地下水の貯留量と流出量の対応関係は貯留関数としてCoutagneとDingにより次式のように示されている2)。
Qg = K2 Sg2(ここで、K：定数、Sg：地下水貯留量)

また被圧地下水の流出量は同じく地下水貯留量と線形関係で得られることが同様に示されていることから、本検討での不圧地下水と被圧地下水の流出量は次式に示す通りとした。

2

qg1 = K g1 × A u (h − Sg )2（式- 3.2.1.7）

q g2 = K g 0 × A g h（式- 3.2.1.8）

q g1 ：不圧地下水流出高（m/s）、 q g2 ：被圧地下水流出高（m/s）
地下水タンクの連続式は次のとおりである。
①

ｈ≧Ｓgの場合∂h= qint − q g1t − q g2t∂t②（式- 3.2.1.9）

ｈ＜Ｓgの場合∂h= qint − q g2t∂t（式- 3.2.1.10）

ここに、qint：地中浸透量高(m)ｈ ：地下水モデルの貯留高(m)qg1t：不圧地下水流出量高(m)Ｓg：不圧地下水流出が発生する高さ(m)qg2t：被圧地下水流出量高(m)474) 河道モデル
河道モデルの構造は図- 3.2.1.9のとおりである。流出はManningの式に従うものとする。
基礎式は式-3.2. 1.11のとおりである。

∂h= qin - q r∂t（式- 3.2.1.11）

ここに、qin：表層・不飽和帯・地下水モデル
河道流及び上流河道モデルからの1
B h 5/3 in

流入量高(m)qr ：河道流出量(m3/s)n：粗度係数 (m-1/3・s)h：河道水位(m)図- 3.2.1.9

河道モデルの概念図

河道の流れはマニングの等流公式が成り立つとすると次のとおりである。

q r = (1/n)h5/3i1/2（式- 3.2.1.12）

ここに、i ：河道（河床）勾配※河道勾配iが1/10万以下となる場合は一律勾配を1/10万としている。

河道流出量高から河道流出量を算出する際には、河道長L(m)および河道幅B(m)を考慮する。河道幅Bは対象とする河道の測量結果が得られれば、そのまま用いるが、測量結果が得られない場合は、レジューム則より河道幅を設定するものとした。

B = cA s（式- 3.2.1.13）

ここに、ｃ，ｓ：定数（ｓは一般に１以下である。）
Ａ：流域面積(m2)
河道長L(m)は落水線の方向により下記のように設定するものとした
落水線の方向縦および横（下図２，４，６，８）
斜め（下図１，３，７，９）

河道長設定
メッシュ長
メッシュ長×√２

１

２

３
６

４
７

８

図- 3.2.1.10

モデル上の河道の設定方法

48９5) 蒸発モデル
蒸発量は、可能蒸発散量に係数を乗じる方法とする。

Eac = fE p

ここに（式- 3.2.1.14）

E ac：実蒸発散量Ep：可能蒸発散量f：係数

係数 f は日毎の可能最大日照時間と実際の日照時間との比率により設定する。
また、本モデルでは可能蒸発散量を推定するモデルとして、Penman の式を採用している。
【Penman の式】
a. 可能蒸発散量の推定

Ep =

Δ
γ
(Qa / L) +f (u )(esa − ea )
Δ +γ
Δ +γ

f (u ) = 0.026(0.5 + 0.54u 2 )ここに（式- 3.2.1.15）（式- 3.2.1.16）

Ep： 可能蒸発散量（cm/day）Qa： 正味放射量（cal/day）

Δ： 飽和水蒸気曲線の温度に対する勾配（hpa/℃）

L： 水の気化熱（＝597.31－0.563T）（T：気温；℃）γ： 乾湿計定数esa： 飽和水蒸気圧（hpa）

ea： 実水蒸気圧（hpa）

f (u ) ： 風速関数
： 地表面から 2ｍの高さの風速（m/s）

u2b. 正味放射量Qa = (1 − α ) Rsc (a + b ⋅ n / N ) − σTa (a′ − b′ ea )(a′′ − b′′ ⋅ n / N )ここにα： アルベド（土地利用・季別で設定、表- 3.2.1.1参照）

Rsc ： 1 日の大気外太陽放射量（cal/day）

n/ N

σ

（式- 3.2.1.17）： 日照率（ n ：日照時間， N ：可照時間）： Stefan‐Boltzman 定数（＝8.132×10－11 cal/min/k－4）49

： 絶対温度（ k ）（＝273.15+T）（T：温度；℃）

Taa ， b ， a ′ ， b′ ， a ′′ ， b′′ ：経験的に決まる定数（下表）
a

b

a′

b′

a ′′

b′′

0.18

0.55

0.56

0.080

0.1

0.9

Rsc = Ro (rm / r ) 2 ⋅ I（式- 3.2.1.18）

I = 1 / π (ϖ o sin φ sin δ + cos φ cos δ sin ϖ o )（式- 3.2.1.19）

ここに
r

Ro ：

太陽定数（1.98 cal/min＝2851.2 cal/day）：

太陽と地球との距離rm ：太陽と地球との平均距離

I：

緯度，太陽の傾き，地形勾配を考慮した係数

φ

：

緯度

δ ： 太陽の傾き（赤道上が 0，北半球が正）ϖo ：cos −1 (tan φ tan δ )土地利用別の月毎のアルベドを表-3.2.1.1 に示す。
表- 3.2.1.1

アルベド設定値

季別(月) 山地
水田
畑
1
0.095
0.160
2
0.100
0.150
3
0.100
0.100
4
0.100
0.080
5
0.120
0.080
6
0.135
0.100
7
0.135
0.130
8
0.135
0.220
9
0.135
0.180
10
0.125
0.150
11
0.120
0.150
12
0.100
0.150

0.190
0.190
0.188
0.188
0.185
0.188
0.193
0.193
0.190
0.185
0.190
0.190

市街地
水域
0.290
0.100
0.290
0.090
0.275
0.080
0.275
0.070
0.260
0.060
0.250
0.060
0.265
0.060
0.275
0.060
0.275
0.070
0.275
0.080
0.275
0.100
0.285
0.110

c. 乾湿計定数（ γ ）γ = C p ⋅ P / 0.622 / L（式- 3.2.1.20）

C p = C pd (1 + 0.846q)（式- 3.2.1.21）

q = 0.622ea / P（式- 3.2.1.22）

ここにC p ：水蒸気圧を考慮した定圧比熱50C pd ································：完全乾燥時の定圧比熱（0.240cal/g/℃）
q

：比湿

p

：大気圧（hpa）

d. 飽和水蒸気圧（ e sa ）飽和水蒸気圧と温度の関係は次の Tetens(1930)の式により算定した。

e sa = 6.11 × 10 aT /( b +T )（式- 3.2.1.23）

水面上（気温 T＞0）

a＝7.5

b=237.3

氷面上（気温 T≦0）

a＝9.5

b=265.3

e. 飽和水蒸気圧の勾配（ Δ ）Tetens(1930)の式では、

Δ = m ⋅ ab /(b + T ) 2 e sa （hpa / ℃）

ここにf.

（式- 3.2.1.24）

m = 1n10

実水蒸気圧（ ea ）ea = C rh ⋅ esa
ここに（式- 3.2.1.25）

C rh ： 相対湿度

516) 融雪モデル
入力条件として与えた日雨量が雨か雪かを判断する基準に日最低気温を用いている。日最低気温が 0℃を下回る場合は降雪と見なし、融雪モデルによる計算を実施する。
融雪のモデルは、熱収支法を採用するものとし、顕熱交換量・潜熱交換量の式は、吉田（1962）の式を用いる。
a.

基礎式

QM = QR + Q A + QE + QP
ここにQM

（式- 3.2.1.26）： 融雪熱量QR ： 放射収支量
Q A ： 大気と雪面の温度差に基づく顕熱交換量
QE ： 雪の蒸発または大気中の水蒸気の雪面への凝結による潜熱交換量
QP ： 雨または雪が与える熱量

b.

融雪熱量

QM = Lm ⋅ q
ここにLm（式- 3.2.1.27）： 氷の融解潜熱（79.7cal/g）q ： 融雪量（g/cm2/hr=10mm/hr：時間単位の場合）

式-3.2.1.26 から融雪熱量を算出し、式-3.2.1.27 より融雪量を求め、この融雪量が表層モデルへ流入するものとする。
c.

放射収支量
放射収支は短波放射収支（第 1 項）と長波放射収支（第 2 項）からなる。長波放射収支

は Brunt の実験による。

{}

Q R = k ( z ) I (1 − α ) + σTs4 − σT 4 ( 0 .51 − 0 .0066 e ) (1 − n ⋅ k )52（式- 3.2.1.28）

ここにk (z ) ： 雪に吸収される短波放射量の割合で積雪深 (z ) の関係〔雪の中の日射の衰弱係数を a とすると k ( z ) = 1 − esp(az ) 〕
σ： Stefan-Boltzman 係数（4.986×10-9cm/cm2/hr/ｋ4）： 緯度，太陽の傾き，地形勾配を考慮した係数

Iα： アルベドTs： 雪の表面の温度（ｋ）

T： 表層付近の温度（ｋ）

e： 地表付近の水蒸気圧（hpa）k： 雲の種類により決まる定数（下層雲：0.86n中層雲：0.77

上層雲：0.21、本モデルでは 0.86 を採用）

： 雲量

d. 顕熱交換量
経験的に次の式で表される。
（式- 3.2.1.29）

Q A = f A (u ) ⋅ T

ここにf A (u )：風速関数で風速の一次式

経験式としては吉田の式等がある。
e. 潜熱交換量
経験的に次の式で表される。
QE = L ⋅ E（式- 3.2.1.30）

E = f E (u ) ⋅ Δ e（式- 3.2.1.31）

ここにE：蒸発量（ｇ/ cm2/hr＝10mm/hr）

f E (u ) ：風速関数で風速の一次式
L： 氷の融解潜熱と水の気化熱（667cal/g）

Δe： 雪の表面と表層付近の水蒸気圧の差経験式としては吉田の式等がある。
f.

吉田の式

{

}

QA = 0.306 1 + 2.03(2 − uτ )3 (T1.2 − T0 )0.5 < uτ < 2.0 （m/s）
Q A = 0.306uτ (T1.2 − T0 )uτ ≥ 2.0 （m/s）
ここに Q A ： 顕熱交換量（cal･cm-3･hr-1）

53（式- 3.2.1.32）（式- 3.2.1.33）

Ta ： 雪面上 a m の高さの気温（℃）（T1.2 はﾒｯｼｭ平均気温、T0 は 0℃とした）uτ ： 風速{}E = 0.844 × 10 −3 1 − 2.03(2 − uτ ) 3 (e1.2 − e0 )（式- 3.2.1.34）

0.5 < uτ < 2.0 （m/s）

E = 0.844 × 10 −3 uτ (e1.2 − e0 )ここにE

（式- 3.2.1.35）-2uτ ≥ 2.0 （m/s）： 蒸発量（g･cm ･hr-1）

ea ： 雪面上 a m の高さの水蒸気圧（hpa）（Tetens の式-3.2.1.25 より算出）
g. 雨または雪が与える熱量（式- 3.2.1.36）

Q p = C ⋅ R ⋅ T / 10

ここにC

： 雨または雪の比熱（cal/g）

R： 雨または雪の強度（mm/hr）

T： 雨または雪の温度（℃）7) ダムモデル
ダム地点のメッシュを特定し、下記３種の設定を行うことが可能である。
設定１：ダム無し（流入量＝放流量）
設定２：ダム放流量入力（ファイルで放流量データを与え、当該メッシュ流出量とする）
設定３：ダム調節量入力（ファイルで調節量データを与え、当該メッシュ流入量＋調節
量を当該メッシュ流出量とする）
本モデルでは上記の内、設定２を採用した。考慮したダムは利根川流域で７ダム、荒川流域で２ダム、多摩川流域で１ダムである。
〔モデルに組み込んだダム〕○利根川流域・奈良俣ダム・矢木沢ダム

・藤原ダム

・薗原ダム

・草木ダム

・下久保ダム

○荒川流域
・二瀬ダム

・滝沢ダム

○多摩川流域・小河内ダム

54

・相俣ダム

(2) 人工系水循環モデルの概要
水循環モデルに取り込む人工系水循環は、以下に示す 3 種類の系統とした。
a)

農業用水

b)

上水道・工業用水道

c)

下水道排水

これらのモデル化方法は、以降に示すとおりである。
1) 人工系水循環のモデル化
a) 農業用水
農業用水として対象となるのは、水田かんがい、畑地かんがい、その他用水等があげられる。この中でも年間を通じた大部分が水田かんがい用水である事から、水田かんがい用水を対象にモデル化する。
水田かんがい用水は、表流水（河川水）を取水し、農業用排水路を経由して各水田で利用する方法と、ポンプにより揚水し注水する方法とがある。ここでは、この両方法をモデル化し取り扱う。水田は、季節により表面の状態が変化するとともに、水が貯まるという特徴がある。
地点別の取水量は、水利権量による取水パターンを設定して、その量を対象となるかんがい区域内の各タンクに面積割合に応じて配分するものとする。
なお、取水後の還元水が河川に戻る地点が特定できないため、取水地点で還元量分を割り引いて取水するものとした。

A1
取水路

Qwr1
A2
排
水
路

Qwr2
.
.
.

河

An

川

Qwrn
Qpw
An+1、 Qwr、

図- 3.2.1.11

Qdr

還元

農業用水モデル化の概念図
55Q

○かんがい取水量・揚水量・還元水量
Q wri =Ai
Q wr
A 1 + A 2 + ... + A n

Qwr = ( 1 －α n )・Qtn
Qdr = βn ⋅ Qwr
ここに
Ai：第 i-th メッシュ面積(m2)Qwri：第 i-th メッシュかんがい流量(m3/day)
Qtn ：河川からのかんがい取水量(m3/day)
Qwr：かんがい用水量(m3/day)
Qpw ：ポンプによる揚水量(m3/day)
Qdr：かんがい還元水量(m3/day)
αn：かんがい取水路漏水率
βn：かんがい用水還元率

○還元率

βn =

Qdr1 Qdr 2
QQ==...= dri = dr
Qwr1 Qwr 2
Qwri Qwr

ここに、h:モデルの貯留高(mm)R:雨量(mm)○かんがい水深⎧ (1 − β n )Qwri α nQtn ⎫hwri = 1000 ⋅ ⎨
+
⎬
Ai
Ai ⎭
⎩

Eps:蒸発量(mm)(mm)

Sf2:表面流の発生する高さ(mm)q0:浸透モデルへの浸透量高(mm)qsf:表面流出高(mm)
qri:早い中間流出高(mm)○表層タンク水位変動量（h>=Sf2）

qin:浸透モデルからの流入高(mm)∂h= R + hwri − E ps − q 0 − q sf − q ri∂t

qg1:不圧地下水流出高(mm)qg2:被圧地下水流出高(mm)Sg : 不 圧 地 下 水 流 出 が 発 生 す る 高 さ○地下水タンク水位変動量（h>=Sg）(m3/day)∂h= qin − q g1 − q g 2 − Q pw / Ai∂tqrin:上流河道からの流入高(mm)○河道タンク河道貯留変動量∂h= q rin + ( Q dr − Q tn ) / A i − q r∂t56b) 上水道・工業用水道
都市生活用水量は、地点別の取水実績データもしくは、これらのデータから予測される取水量を用いる。同じく工業用水量は、取水点の実績データもしくはこのデータから予測される取水量を用い、給水区域で消費されるものとしている。
（利根大堰→秋ヶ瀬といった河川間の導水等についても取水地点と排水地点をつなぐことによりモデル化している）○不飽和・地下水の減水深Δhss = α k qtk + α c qtc ：不飽和タンク減水深Δh sg = − q pk − q pc

：地下水タンク減水深○不飽和タンク水位変動量（Ss1 =< h < Ss2）θs∂h= qin − E ps − q s 0 − q s1 + Δh ss∂t○地下水タンク水位変動量（h>=Sg）∂h= qin − q g1 − q g 2 + Δh sg∂t○河道タンク河道貯留変動量∂h= qrin − qr − qtk − qtc∂t

ここに、
αk:工業用水取水管の漏水率
αc:都市生活用水取水管の漏水率
qtk:工業用水取水量高(mm)qtc:都市生活用水取水量高(mm)
qpk:地下水からポンプによる工業用水揚水量高(mm)
qpc: 地下水からポンプによる生活用水揚水量高
(mm)

A1
取水路

qtk1, qtc1
A2
河川

qtk2, qtc2
.
.

B

.

河

An

川

qtkn, qtcn

A

図- 3.2.1.12

Q

上水道・工業用水道のモデル化の概念図
57c) 下水道排水
下水管に地下水が浸出していることは、晴天日における汚水発生量と下水処理量に差があるため、この点を地下水タンク・河道タンクへ反映する。
○地下水の減水深Δhc = β c q dc
○地下水タンク水位変動量（h>=Sg）∂h= qin − q g1 − q g 2 + Δhc∂t○河道タンク河道貯留変動量
ここに、∂h= q rin − q r + Qdc / A∂tβ c: 下 水 管 へ の 地 下 水 の 浸 出 率(0.13)
Qdc:処理場生活汚水排水量(m3/s)
A : 処理場集水面積(m2)

集水エリア※Qdc の流入先：
処理場位置により、河道タンクに流入するメッシュ、落水線に流入する
メッシュ、東京湾に直接流入するメッシュに区分している。

処理場

Qdc

河川

図- 3.2.1.13

下水道排水のモデル化の概念図

58d) 人工系水循環モデルの組み込み
人工系の水循環を自然系水循環モデルに考慮する方法として、以下①～③のゾーン を想定する。
① 表層タンクに位置し、供給された水量は早い中間流と不飽和帯への供給源となる。
② 不飽和帯タンクに位置し、供給された水量は遅い中間流及び不圧地下水の供給源となる。
③ 地下水タンクに位置し、供給された水量は不圧地下水の供給源となる。
これらの仮定をもとに、各人工系の諸量を、表- 3.2.1.2の通り取り扱うものとする。
表- 3.2.1.2

人工系諸量の取り扱い

カテゴリ

取り扱い

農 業 用 水 ： ①に適用する。表面流としての流出はないが、中間流の一部及び地下浸透の源と仮定する。
上水道用水： 漏水量が対象となり、②に適用する。
工 業 用 水 ： 同上。
下水道用水： ③に適用する。下水道は基本的に地下水面より低い位置に埋設される事を考慮する。
ポンプ揚水： ③に適用する。

59

表層タンクの概念図
表面流

K f2・(1/n)・(h - S f2) 5/3・ i1/2

農業用水

中間流と地下浸透へ

早い中間流

Sf 2

K f0・ α ri・ α ・
0 (h - S f1 )

①
Sf 1

Sf 0

地下浸透

K f0・ α ・0 (h - S f0 )上水道・工業用水

漏水が不飽和帯タンクへ

不飽和帯タンクの概念図
h

D

Sf 2

②

遅い中間流

K s1・ K ・
x D・ i

Sf 1
下水道

地中浸透

地下水タンクへK s0・ K Z

地下水タンクの概念図h
不圧地下水流ポンプ揚水③

地下水タンクへ2

K g1・ A u ・(h - S g ) 2
被圧地下水流Sg

K g0・ A ・
g h

落水線へ

図- 3.2.1.14

人工系水循環モデルの組み込み概念図

60

Q 取水
（消費）

河道タンク

表層タンク

Q'
不飽和タンク
Q 取水
（消費）

漏水q

地下タンク

図- 3.2.1.15

上水道用水のモデル化概念図(1) 取水量の推定（Q’m3/日）
各計算ブロックについて、漏水量（q m3/日）を推定するために日取水量を推定し、日取水量に漏水率を乗じて漏水量を算出し、シミュレーションの入力データとした。
漏水量（q m3/日）＝日取水量（m3/日）×漏水率
日取水量は、月別取水量（A m3/日）をもとにし、市町村人口（B 人）割合で市町村別日取水量（C m3/日）を求め、この値を市町村面積（D km2）を用い、1 格子面積（1 km2）に分け、各格子の平均の取水量（E m3/月）とした。さらに、この取水量を月間日数で除して日取水量（Q’m3/日）を求めた。
(2) 漏水率
漏水率は、年間給水量の合計と無効水量の合計の比較より 0.7%とした。
ここで各種資料の出典は、以下の通りである。
月別取水量 A(m3/日)
水道統計
市町村人口 B(人)
都県勢要覧
市町村面積 D(km2)
都県統計年鑑
年間取水量合計
水道統計
H13
茨城
栃木
群馬
埼玉
千葉
東京
神奈川
合計

年間給水量(千 m3)
306,910,000
243,639,000
315,563,000
898,498,000
658,541,000
1,704,373,000
1,164,913,000
5,292,437,000

年間無効水量(千 m3) 漏水率
26,0020.8%
33,258
1.4%
37,1491.2%
62,764
0.7%
39,882
0.6%
116,392
0.7%
77,3180.7%
392,7650.7%61

表流水取水
QR

河道タンク

表層タンク

水道水取水
Qw

不飽和タンク
漏水q

地下水取水
QG

地下タンク

図- 3.2.1.16

工業用水のモデル化概念図

工業用水道として利用される水量の取水内訳は、河川水（地表水）、水道水及び地下水に区分される。
工業用水(QM)＝地表水(QR)＋水道水(QW)＋地下水(QG)
工業統計表では、県別の年間の水源別用水量（QR、QW 、QG）が取りまとめられていることより、
これらの諸量を市町村別に、工業出荷額の規模に比例させて配分する。
次に各市町村毎の水源別工業用水を市町村面積で除して 1km2 当りの用水量を求め、その値を各ブロックの用水量とした。さらに各ブロックの水源別用水量を年間日数で除して日量に換算し、シミュレーションの入力データとした。
水源別の取水方法、漏水量の取り扱いは①上水道用水と同一と仮定する。ここで漏水率は、0.7%
とした。

62

河道タンク

表層タンク

不飽和タンク

排水
地下タンク
浸出q

図- 3.2.1.17

下水道のモデル化概念図(1) 汚水処理量
下水道統計より、県別の年間汚水処理量（QMm3/年）が集計されている。この値をもとに、各市町村の人口比で市町村別汚水処理量を算出した。さらに、市町村別汚水処理量を市町村面積で除し、市町村別にブロック１km2 あたりの年間汚水処理量を求め、この値を年間日数で除して日平均の処理水量とした。
(2) 浸出量
地下水の下水管路への浸出量は、下水道統計を用い、年間の汚水処理量と総有収水量の差から求めた。
地下水浸出量(m3/年)＝年間総汚水処理量(m3/年)―総有収水量(m3/年)
得られた地下水浸出量を市町村人口比で各市町村に配分し、さらに市町村面積及び年間日数で除して各市町村の 1 ブロック当りの地下水浸出量とした。
H13
茨城
栃木
群馬
埼玉
千葉
東京
神奈川
合計

年間汚水処理量(百万 m3)
231,709146,304
125,401
602,699452,956
1,789,244
1,155,183

年間総有収水量(百万 m3)
200,059
104,522
102,704
570,116
377,024
1,646,602
917,010

漏水量(百万 m3)
31,65041,782
22,69732,58375,932
142,642
238,173

漏水率
14%
29%
18%
5%
17%
8%
21%

4,503,496

3,918,037

585,45913%

合計の漏水率 13％を浸出率とした。

63

BQ
水田

表流水取水
Q

河道タンク

(1-B)Q

表層タンク

不飽和タンク

地下水取水
(Q)

地下タンク

図- 3.2.1.18

農業用水のモデル化概念図(1) 農業用水取水量
農業用水取水量(m3/日)＝かんがい面積(km2)×期別減水深(mm/日)×1000/(1-水路損失率)(2) 期別減水深農業土木ハンドブックより、標準的な値として次頁(表-3.2.1.3)の数値を用いる。
(3) かんがい面積
都県農林水産統計年報より、各年の市区町村別水田面積および畑地面積を用いる。
(4) 水路損失15%を見込んだ。
(5) 農業用水取水量
これらの値をもとに、市町村の農業用水取水量を算定し、市町村面積をもとに、シミュレー
ションブロック 1km2 あたりの水量を算出し、農業用取水量として、設定した。
取水水源は地表水（河川水）と地下水（ポンプ取水）に区分される。その振り分けは、平成15 年度水資源白書の地域別用途別地下水依存率から「地表：地下＝81：19」とした。さらに、取水量に対して還元率 B を考慮し、田面からの蒸発及び地下浸透で消費される分は、(1-B)×Qとなり、B×Q は直接河川に還元されるものとした。

64

表- 3.2.1.3

期別減水深一覧［水田］期

別

早期栽培（mm/日）

普通期栽培（mm/日）

4/21～30

代かき（150）

5/01～1025

5/11～2026

5/21～3027

5/31～6/09

28

6/10～6/19

28

代かき（150）

6/20～6/29

27

26

6/30～7/09

28

28

7/10～7/19

27

28

7/20～7/29

29

29

7/30～8/08

27

28

8/09～8/18

26

28

8/19～8/28

27

8/29～9/07

27

9/08～9/17

28

9/18～9/27

26

9/28～10/07

25
9.6

上記以外

9.6

早期・普通期栽培の面積割合は、それぞれ 50%ずつとした。

［畑地］
期

別

減水深（mm/日）

4/1～9/30

4

上記以外

0

65

2) 合流式下水道の取扱い
a) 合流改善に関する評価
合流改善事業は、合流式下水道において降雨の初期段階における汚濁濃度の高い汚水を公共水域へ放流することを極力抑えることを目的としている(貯留施設による対応やポンプ所などの施設運用による対応など)。東京都においては、晴天時汚水の流量(1Q)に対して
ある程度までの初期雨水を貯留するものとしている。合流改善事業の概要について以下に整理した。
¾

東京都の基本対策：下水道は 3Q 以上（総雨量 5mm 程度）で雨水流出となる。（雨水ポンプにて排水を
開始する）

¾

現状では 5mm 程度までは放流が行われない（3Q 対応の場合、処理場で処理されて廃水される）が、そ¾

降雨開始から累加雨量 5mm までは現状の処理場での処理を経由して排水され、ポンプ運転による初期雨れを超えると直接雨水が河川に排水され、生活排水も含まれることとなる。
水の排水は行わない。
¾

累加雨量が 5mm を超過した場合には初期汚濁を含めてポンプにて直接河川に排水される。

¾

合流改善事業では初期雨水の 8mm までは流域貯留を対策案として用いる。また、リアルタイムコントロール(RTC)による対策を推進中。

¾

東京都の下水道区域は 82％が合流式、79 ポンプ所で 100,000m3 / 分（1,670 m3/s）を排水。

¾

合流改善事業で、低平地での事例では、年間約 75 回の降雨に対し、ポンプ運転が 55 回行われており、雨量が 5mm を超過した段階で排水が始まる。つまり、5mm を超過する降雨回数が 55 回見込まれる (江東区での事例流域での実績値)。
→（例えば H13 の東京（大手町）では 1mm 以上が 97 日、5mm 以上が 63 日、20mm 以上が 25 日、最大で 186mm が観測されている。）

b) モデルでの表現方法
合流式下水道が流入する下水処理場について、ある一定以上の降雨があった場合に、汚水が処理場を経由されず直接公共水域へ排水されるという現象をモデルへ反映させる。
したがって、日雨量に対する判断基準値を設定し、基準値を超える際には処理場からの放流負荷濃度を変化させるという方針で合流改善の施策をモデルに反映させることとした。
反映方法を以下に整理した。
¾

下水道区域について、日降雨量 5mm 以下については現状のモデルでの処理場からの処理後の負荷量排出のみとする。排水量も処理場からの排水量で表現する。

¾

日雨量で 5mm を超過する日の当該日負荷量は 5mm 超過分について未処理相当水質分が各メッシ
ュからポンプ場に集計され、ポンプ場から排水される。
（処理水と未処理水の比率に応じ、排水水質を決定する）¾

雨水排水と汚水排水は異なる場所で行われるケースも多いが、ここでは同一地点で排水されるものとする。

¾

合流式下水道雨天時越流水(CSO)対策は東京都を例として、5mm 対策（3Q）であったものをクイ
ックプラン対応の 8mm に引き延ばす対策を第一とし、次いで、RTC 制御などの導入で 25mm 対応
までの対策を適用する。

66

降雨規模（日雨量）毎の晴天時流量に対する雨水流量の比率、処理水と未処理水の割合を図- 3.2.1.19に示した。この処理水と未処理水の割合に応じて下水処理場からの排水水質を変化させるものとした。
降雨時負荷量 ＝ 通常時負荷量 ＋ 通常時負荷量 × 比率A × 比率B
比率A：処理水質と未処理水質の比率（未処理水質／処理水質）
比率B：下記で算定される、降雨規模別の処理水量と未処理水量の比率（未処理水質／処理水質）
ケース①現状：5mm・３Qまでは全て処理し、３Qを超える流量が未処理で排出されると想定。
ケース②クイックプラン対応：8mm・4.2Qまでは全て処理し、4.2Qを超える流量が未処理で排出されると想定。
ケース③RTC制御対応：25mm・11Qまでは全て処理し、11Qを超える流量が未処理で排出されると想定。

1.000
未処理水量の占める割合（比率B）

0.900
0.800
0.700
0.600

現状
クイックプラン
RTC制御

0.500
0.400
0.300
0.200
0.100
0.000
0

図- 3.2.1.1950100
降雨規模(mm)

150

200

合流式下水道の処理場からの排水負荷量算定方法

67

(3) 流域物質循環モデルの概要
1) 物質循環モデルの基礎式
物質流動モデルの基礎式を以下に示す。
a) 表層物質流動モデル
降雨負荷、面源負荷の地表面堆積・掃流、堆積物への吸脱着等を下記のモデルにより表現する。表層に堆積・吸着する汚濁負荷は表層流（表面流＋早い中間流）の掃流力によって流出するものと考える。堆積物の流出負荷量は流量のべき乗に比例するものとした。また、流出負荷量と同じ水質で地中浸透するものとした。
①運動式

L f 2 = k np ⋅ S np m ⋅ A ⋅Q f 2 n(式- 3.2.1.37)

L f 1 = L f 2 / Q f 2 ⋅Q f 1(式- 3.2.1.38)

②連続式

A⋅

dS np
dt

= S u ⋅ (1 − e − Ks⋅T ) + C r ⋅ R ⋅ A − k np ⋅ S np ⋅ A ⋅ Q f 2 − L f 2 / Q f 2 ⋅ Q f 1m

ここに、 L f 2n： 表層モデルの流出負荷量（g/d）

Lf 1： 表層モデルから不飽和タンクへの浸透負荷量（g/d）

Qf 1： 表層モデルの流出量（m3/d）

Q f 2 ： 表層モデルから不飽和タンクへの浸透量（m3/d）
k np： 面源由来堆積物の掃流係数（d/m6）

S np： 面源由来堆積物量（g/m2）

Su： 極限堆積物量（g）

Ks： 堆積速度係数(1/d)

T： 堆積 0 からの経過日数(d)Cr： 降雨濃度（g/m3）

R： 降雨量（m/d）

A： 面積（m2）

m, n ： 定数68(式- 3.2.1.39)

b) 不飽和層物質流動モデル
表層モデルより地中浸透した汚濁負荷は土壌内に蓄積し、蓄積物に吸脱着する過程を下式により表現する。蓄積物から溶脱した汚濁負荷は中間流によって流出する。また、流出負荷量と同じ水質で地中浸透するものとした。
①運動式

L s 2 = k s ⋅ S s m ⋅ A ⋅Q s 2 n(式- 3.2.1.40)

L s 1 = L s 2 / Q s 2 ⋅Q s 1(式- 3.2.1.41)

②連続式

A⋅

dS s
m
n
= S su ⋅ (1 − e − Kss⋅T ) − k s ⋅ S s ⋅ A ⋅ Qs 2 − Ls 2 / Qs 2 ⋅ Qs1 + L f 2 / Q f 2 ⋅ Q f 1 (式- 3.2.1.42)dt

ここに、 L s 2： 不飽和層モデルの流出負荷量（g/d）

Ls1： 不飽和層モデルから地下水モデルへの浸透負荷量（g/d）

Qs 2： 不飽和層モデルの流出量（m3/d）

Q s1： 不飽和層モデルから地下水モデルへの浸透量（m3/d）

ks： 不飽和層モデルの吸脱着速度係数（1/d）

A： 面積（m2）Ss： 不飽和層モデルの汚濁負荷の土壌内蓄積量（g/m2）

S su： 極限溶脱物量（g）

Kss： 溶脱速度係数(1/d)

T： 溶脱 0 からの経過日数(d)m, n ： 定数69c) 地下水層物質流動モデル
不飽和層モデルより地下浸透した汚濁負荷は土壌内に蓄積し、蓄積物に吸脱着する過程
を下式により表現する。蓄積物から溶脱した汚濁負荷は基底流によって流出する。
①運動式

Lg = k g ⋅ S g ⋅ A ⋅ Qgmn

(式- 3.2.1.43)

②連続式

A⋅

dS g
dt

ここに、 L g= S gu ⋅ (1 − e − Kgs⋅T ) − k g ⋅ S g ⋅ A ⋅ Q g + Ls 2 / Qs 2 ⋅ Qs1mn： 地下水層モデルの流出負荷量（g/d）

Qg： 地下水層モデルの流出量（m3/d）kg： 地下水層モデルの吸脱着速度係数（1/d）

A： 面積（m2）

Sg： 地下水層モデルの汚濁負荷の土壌内蓄積量（g/m2）

S gu： 極限溶脱物量（g）

Kgs ： 溶脱速度係数(1/d)
T： 溶脱 0 からの経過日数(d)m, n ： 定数70(式- 3.2.1.44)

d) 河道物質流動モデル
各層のタンクからの流出負荷の自浄作用による河床への吸着等を含む河床堆積物が流量に比例して次メッシュへ流出するものとする。また、点源負荷は河道に直接排出され（河道タンクがないメッシュについては落水線に排出され）、各層からの流出負荷と共に流量に応じて次メッシュへ流出するものとする。
①運動式

LRO1 = C R ⋅ QRO(式- 3.2.1.45)②連続式

CR = ( LRI 1 + LI 1 + L p ) / QRI ⋅ e − ( k1+ k 2 )・xt / uu + Cka ⋅ (1 − e − ( k1+ k 2 )・xt / uu ) /( k1 + k 2)

(式- 3.2.1.46)

ここに、 LRO1 ： 次メッシュへ流出負荷量（g/d）

CR： 河道モデルの流出水濃度（g/m3）

Q RO ： 河道モデルの流出量（m3/d）
k1： 脱酸素係数（1/d）

k2： 沈降・吸着による除去速度係数（1/d）

QRI： 河道モデルへの流入量（m3/d）

xt： 再懸濁による負荷速度定数（1/d）

Cka： 河道内汚濁負荷濃度（g/m3）uu： 河道水位（m）

QRO / B / L （B:河道幅, L:河道長）

LRI 1 ： 前メッシュからの流入負荷量（g/d）
LI 1： 流域（各層モデル）からの流入負荷量（g/d）

Lp： 点源（生活系・工業系・畜産系）排出負荷量（g/d）

71

2) 人工系物質循環モデルの概要
家庭排水による負荷量は、生活者の生活習慣・食習慣などを反映するので、地域（たとえば都市と農村）によって必ずしも同じではない。家庭排水は、し尿とその他の雑排水に分けられるが、下水道整備地域では両方とも下水処理場から排出される。しかし、未整備地域では雑排水は発生源で排出されるが、し尿は処理方式によって排出量・排出点が異なる。汲取り処分地
域では、し尿は処理場から排出されることになり、処理場の位置によっては他流域へ持ち出されることもある。農地還元されている場合は、農地からの流出としてカウントされなければならない。し尿単独処理浄化槽で水洗化されている家庭では、一部浄化されて発生源で排出される。合併処理浄化槽によって雑排水とし尿が処理されている場合（集中浄化槽・コミュニティ
プラント・農村下水道など）は、両方とも処理されて発生源の近くから排出される。さらに、都市域では昼間人口の評価、観光地では入込み人口の評価などが必要になる（ 図- 3.2.1.20参
照）。

下水処理場
発生源排出
発生源

小型合併処理浄化槽

集中排出
非特定発生源排出

雑排水

し尿

し尿処理場
無処理

し尿単独処理
浄化槽

汲取り

海洋投棄
山林投棄

農地還元

図- 3.2.1.20

家庭排水の発生形態

ここでは、生活系排水、事業所系排水、工場系排水及び畜産系排水の汚染源別に、それぞれ
点源負荷系汚濁物質流動モデルを構築する。人工系物質循環モデルの設定方法を以下に示す。

72a) 生活系排水モデル
生活系排水の汚濁負荷でモデルに取り込む必要のある汚濁負荷源は、下記に示す通りである。
① 下水処理場② 農業集落排水処理場③ 合併浄化槽④ し尿処理場⑤ 単独浄化槽⑥ 自家処理⑦ 雑排水処理
これらの負荷は、年間を通じて大きな変動のないものと考えられるため、時間的に一定の負荷量が近傍の落水線に流出するものとして、物質循環モデルを構築する。
① 下水処理場
下水処理場のモデル化は、処理場の位置、処理人口実績及び 1 人当たりの負荷量原単位をもとにしてモデルに組み込むものとする。モデル上では、下水処理場の排水口が位置しているメ
ッシュの河道（河道タンクがないメッシュでは落水線）に汚濁負荷が流出するものと考える。
② 農業集落排水処理場
処理人口実績及び 1 人当たりの負荷量原単位をもとにして、モデルに組み込むものとする。
モデル上では、それぞれの市町村を構成する全メッシュに平均的に汚濁負荷を与え、近傍の落水線に汚濁負荷が直接流出するものと考える。
③ 合併浄化槽
各合併浄化槽について、市町村毎の処理人口実績および 1 人当たりの負荷量原単位をもとに、モデル化を行う。モデル上の設定方法は②と同様とする。
④ し尿処理場
処理場の位置、稼働実績および市町村毎の処理人口実績をもとにモデルに組み込むものとする。モデル上の設定方法は②と同様とする。
⑤ 単独浄化槽
各単独浄化槽の市町村毎の処理人口実績、および 1 人当たりの負荷量原単位をもとにモデルに組み込むものとする。モデル上の設定方法は②と同様とする。
⑥ 自家処理
市町村毎の自家処理人口実績、および 1 人当たりの負荷量原単位をもとにモデルに組み込むものとする。モデル上の設定方法は②と同様とする。
⑦ 雑排水処理流域で実施されている市町村毎の雑排水処理人口実績および 1 人当たりの負荷量原単位をもとに、モデルに組み込むものとする。モデル上の設定方法は②と同様とする。
なお、し尿処理場、単独浄化槽及び自家処理は、し尿のみを対象としているため、別途家庭排水を対象に雑排水処理を行うこととなる。したがって、市町村毎のし尿処理場対象人口、単独浄化槽処理人口及び自家処理人口の和は、雑排水処理人口に等しくなる。

73b) 事業所系排水モデル
事業所系排水の汚濁物負荷についても生活系排水と同様の方法で物質循環モデルを構築する。データとしては、市町村別に排水量や汚濁負荷濃度が存在するので、それぞれの市町村を構成する全メッシュに平均的にこれらの負荷量を入力し、近傍の落水線に汚濁負荷が直接流出するものと考える。

c) 工業系排水モデル
工場系排水の汚濁負荷についても生活系排水と同様の方法で物質循環モデルを構築する。データとしては、市町村別、業種別に排水量や汚濁負荷濃度が存在するので、それぞれの市町村を構成する全メッシュに平均的にこれらの負荷量を入力し、近傍の落水線に汚濁負荷が直接流出するものと考える。

d) 畜産系排水モデル
畜産系排水の汚濁物負荷については、生活系排水と同様の方法で物質循環モデルを構築する。
データとしては、市町村別、家畜別に汚濁負荷量が存在するので、それぞれの市町村を構成する全メッシュに平均的にこれらの負荷量を入力し、近傍の落水線に汚濁負荷が直接流出するものと考える。

74

3.2.2水域モデル（流動・水質モデル）の概要(1) 流域モデルと水域モデル接続の目的
1) 対象とする時空間スケールについて
東京湾のような閉鎖性水域の水環境の再生について検討する場合、また、そのための流動・水質モデルを選定する場合、取り扱う現象やその時空間スケールを明らかにする必要がある。
例えば、河口沿岸域や汽水域の干潟、藻場等の水質（濁り、貧酸素化）や底質（粒径、Eh、IL 等）の改善を議論する場合には、塩淡成層、河口流出流、波浪変形等による物質移動を把握する事が望ましい。
一方、東京湾の広域的な水質レベルの改善を議論する場合には、湾内の富栄養化現象を扱うことが重要である。
前者では、数 m～数十 m の空間スケールで数時間～数日間の時間スケールを扱うモデルが必要であり、後者は、数百 m～数 km の空間スケールで数ヶ月間の時間を扱うモデルが必要である。 今 回 の 検 討 で は 東 京 湾 全 体 の 流 動 や 水 質 の 評 価 と 通 年 を 通 じ て検討を行うことを念頭にモデルの構築を行った。
表- 3.2.2.1
時
間
スケール
0

10

1

10
1分

エスチャリーにおける水質変動に関わる現象の時間スケール 2)
流体運動

自然現象
気象・地象
生物的変化

渦、乱流

人為現象
化学的変化
懸濁物質による
吸脱着
酸化・還元

底泥巻上げ
風波
うねり

2

10

BOD菌の増殖
脱窒菌の増殖

懸濁物質の沈降
静振

酸素の溶解

3

10
1時間

4

10
1日

潮位変化

気温変化
日射量変化
水温変動
風
降雨
河川水量変化

5

10
1周

6

10

硝化菌の増殖植物プランクトンの上下動
植物プランクトンの増殖
魚群の衆参

温度躍層の季節
変化
7

貧酸素化

大潮・小潮

1月
10

生活排水による
負荷

肥料流亡
農業利水

大大潮・小小潮

1年
8

10

突発的

気象長期変動地盤変動

生物資源の長期変化

台風

土地利用変化
汚濁物質の排出

75

広域的
水質レベル

2) 水域モデルにおいて求められる要件
閉鎖性水域の水物質循環に関わる施策（流域負荷削減、干潟再生など）による湾内の水質改善効果を評価するには、次の要件を満たす必要があると考えられる。
a) 流域水物質循環に係わる施策の評価
水質季節変化を表現できるモデルである必要がある。
b) 干潟の再生等の施策の評価
底泥の溶出や浅場の浄化機能を表現する項目が必要である。
c) 今後、下水道の合流改善等に関わる施策にも使える
大腸菌群数の移流拡散計算が可能なモデルである必要がある。
d) 水質改善効果をわかりやすい指標で表現
水質改善効果を透明度で表現する。
(2) 流動モデルの構造
1) 流れの基礎方程式a) 連続方程式
空間的に密度が変化する水域における連続方程式は、均質流体の質量保存式
Dρ+ ρdivU = 0
Dt（式- 3.2.2.1）

に対して、溶質の拡散を考慮した以下の式が成立する。
Dρ+ ρdivU = K∇ 2C
Dt（式- 3.2.2.2）

また、濃度 C の保存式は以下のようになる。
DC+ CdivU = K∇ 2C
Dt（式- 3.2.2.3）

ただし、記号はそれぞれ以下に示すとおりである。
∂
∂
D
∂
∂
= +u +v +w
∂z
∂x
∂y
Dt ∂t
∂2
∂2
∂2
+
+
∂x 2 ∂y 2 ∂z 2
∂
∂
∂
div =
+
+
∂x ∂y ∂z

∇2 =

t

：時間

x，y，z

：空間座標（直交座標系）

ρ

：密度

U

：速度ベクトル

u，v，w：速度ベクトルの x，y，z 方向成分

K

：拡散係数

C

：溶質濃度
76

ここで、密度変化の原因となる温度及び濃度の変化がそれほど大きくなく、かつ水は非圧縮性を仮定し、圧力による変化を無視すると、式-3.2.2.2 と式-3.2.2.3 から式-3.2.2.4 が導かれる。
（式- 3.2.2.4）

divU = 0

したがって、密度差が存在する流れの場においても密度が一様な非圧縮性流体と同様の連続式が適応可能となる。
b) 運動方程式
x，y 軸を水平方向、z 軸を鉛直方向に設定し、鉛直方向には速度変化及び加速度が小さ
く、静力学平衡が成り立つと仮定すると、x，y 方向の運動方程式（flux form）は以下のようになる。
∂u∂ 2∂∂∂ζ g 0 ∂ρ
1 ∂p
=−
− ∫
u − (uv ) − (uw) + f 0 v − gdz − ⋅ 0
ρ ∂x
∂t
∂x
∂y
∂z
∂x ρ z ∂x( )∂u ⎞ ∂ ⎛ ∂u ⎞∂ ⎛∂u ⎞ ∂ ⎛⎟ + ⎜ Nz
+ ⎜ N x ⎟ + ⎜⎜ N y⎟∂z ⎠∂y ⎟⎠ ∂z ⎝∂x ⎝∂x ⎠ ∂y ⎝
∂u
∂
∂ 2 ∂
∂ζ g 0 ∂ρ
1 ∂p= − (uv ) −− ∫v − (vw) − f 0u − g
dz − ⋅ 0
z
ρ ∂y
∂t
∂x
∂y
∂z
∂y ρ ∂y( )∂v ⎞ ∂ ⎛ ∂v ⎞∂ ⎛ ∂v ⎞ ∂ ⎛+ ⎜ N x ⎟ + ⎜⎜ N y ⎟⎟ + ⎜ N z ⎟∂y ⎠ ∂z ⎝ ∂z ⎠∂x ⎝ ∂x ⎠ ∂y ⎝（式- 3.2.2.5）（式- 3.2.2.6）

ここで、
f0

：コリオリパラメータζ：平均水面から自由水面までの高さg：重力加速度p0：大気圧

Nx，Ny，Nz：x，y，z 方向の渦動粘性係数c) 熱収支の式
熱についての移流拡散方程式を用いる。
∂T∂∂∂= − (uT ) − (vT ) − (wT )
∂t
∂x
∂y
∂z
∂ ⎛ ∂T ⎞ ∂ ⎛ ∂T ⎞ ∂ ⎛ ∂T ⎞⎟ + ⎜ kz
+ ⎜ kx⎟ + ⎜ky⎟∂x ⎝ ∂x ⎠ ∂y ⎜⎝ ∂y ⎟⎠ ∂z ⎝ ∂z ⎠

ここで、
T：水温
kx，ky，kz：x，y，z 方向の熱に関する乱流拡散係数77（式- 3.2.2.7）

d) 塩素量収支の式
塩素量についての移流拡散方程式を用いる。
∂Cl
∂
∂
∂
= − (uCl ) − (vCl ) − (wCl )
∂t
∂x
∂y
∂z
∂Cl ⎞ ∂ ⎛ ∂Cl ⎞∂ ⎛ ∂Cl ⎞ ∂ ⎛⎟ + ⎜ Kz
+ ⎜ Kx⎟ + ⎜⎜ K y⎟∂x ⎝∂x ⎠ ∂y ⎝∂y ⎟⎠ ∂z ⎝∂z ⎠（式- 3.2.2.8）

ここで、
Cl：塩素量Kx，Ky，Kz：x，y，z 方向の塩素量に関する乱流拡散係数e) 状態方程式
塩素量と水温から密度を計算するための状態方程式として Knudsen の式を用いる。

ρ=

σ

t
+1
1000
σ t = Σ t + (σ 0 + 0.1324){1 − At + Bt (σ 0 − 0.1324)}

σ 0 = −0.069 + 1.4708Cl − 0.001570Cl 2 + 0.0000398Cl 3
Σt = −(T − 3.98)2 ⋅ T + 28.30

503.570 T + 67.26
At = T (4.7867 − 0.098185T + 0.0010843T 2 )×10 −3

Bt = T (18.030 − 0.8164T + 0.01667T 2 )×10 −6782) マルチレベルモデル密度流モデル
基礎方程式、式-3.2.2.5～式-3.2.2.8 をセル内で Z 方向に積分すると、以下のようなマル
チレベルモデルの基礎式が得られる。
a) 運動方程式∂∂∂M k≅ − (M k vk ) − (M k vk ) − (uw) z = − H k −1 + (uw) z = − H k + f 0 N k
∂y
∂x
∂t

[ ]

1
∂M k ⎞
hk ⎧ ~∂ρ k ⎫ ∂ ⎛ ∂M k ⎞ ∂ ⎛⎟⎟ + ⎜⎜ N y⎨ Px k − ghk⎬ + ⎜ Nx
ρk ⎩
2
∂y ⎟⎠
∂x ⎠ ∂y ⎝∂x ⎭ ∂x ⎝1
1+ τ xk −1,k − τ xk ,k +1+

ρ（式- 3.2.2.9）

ρ

∂N k
∂
∂
≅ − ( N k vk ) − (N k vk ) − (uw) z = − H k −1 + (uw) z = − H k − f 0 M k
∂t
∂x
∂y

[ ]

hk ⎧ ~∂N k ⎞∂ρ k ⎫ ∂ ⎛∂N k ⎞ ∂ ⎛1⎟⎟ + ⎜⎜ N y⎨ Py k − ghk⎬ + ⎜ Nx
ρk ⎩
∂y ⎟⎠
2
∂y ⎭ ∂x ⎝∂x ⎠ ∂y ⎝11+ τ yk −1,k − τ yk +1,k+ρ

（式- 3.2.2.10）

ρb) 連続方程式（式- 3.2.2.11）

wk = 0 (k = K )wk +1 = wk −∂M k ∂N k−∂x∂y(k = 2,3LL, K )∂ζ
∂M 1 ∂N1= w1 −−∂t∂x∂x（式- 3.2.2.12）（式- 3.2.2.13）

79c) 熱収支の式（top-レベル）∂(h1T1 ) = − ∂ (M 1T1 ) − ∂ (N1T1 ) + (wT ) z =− H1
∂t
∂x
∂y
∂ ⎛
∂T ⎞ ∂ ⎛∂T ⎞ Q ⎛ ∂T ⎞+ ⎜ h1k x 1 ⎟ + ⎜⎜ h1k y 1 ⎟⎟ − s − ⎜ k z⎟ z = − H1∂x ⎝∂x ⎠ ∂y ⎝∂y ⎠ cv ρ ⎝ ∂z ⎠（式- 3.2.2.14）（他のレベル）∂(hkTk ) = − ∂ (M kTk ) − ∂ (N k Tk ) − (wT ) z =− H k−1 + (wT ) z =− H k
∂t
∂x
∂y
+

∂T ⎞
∂ ⎛
∂T ⎞ ∂ ⎛⎜ hk k x 1 ⎟ + ⎜⎜ hk k y k ⎟⎟∂x ⎝∂x ⎠ ∂y ⎝∂y ⎠（式- 3.2.2.15）⎛ ∂T ⎞⎛ ∂T ⎞+ ⎜ kz⎟ z = − H k −1 − ⎜ k z⎟ z =− H k⎝ ∂z ⎠⎝ ∂z ⎠d) 塩素量収支の式∂(hk Clk ) = − ∂ (M k Clk ) − ∂ (N k CLk )∂t∂y∂x− (wCl ) z = − H k −1 + (wCl ) z = − H k∂Clk ⎞∂ ⎛∂Cl1 ⎞ ∂ ⎛⎟+ ⎜ hk K x⎟ + ⎜⎜ hk K y∂x ⎝∂x ⎠ ∂y ⎝∂y ⎟⎠（式- 3.2.2.16）⎛ ∂Cl ⎞⎛ ∂Cl ⎞
+ ⎜ Kz⎟ z =− H k⎟ z = − H k −1 − ⎜ K z
∂z ⎠
∂z ⎠
⎝
⎝

ここで変数の添え字 k は、第 k 番目のセルの値を意味する。またセル内の線流量 M，Nを以下のように定義する。
Mk ≡ ∫− H k −1

−Hkudz , N k ≡ ∫− H k −1

−Hk

[P~ ] :
[P~ ] = − ∂∂px − ρ g ∂∂ζx

vdz

x k

0

x 1

1

[P~ ] = [P~ ] − gh ∂ρ∂xx k

x k −1k −1

k −1( k > 2)y 方向も同様
hk：第 k 番目のセルの高さ（ = H k − H k −1 ）τ xk −1, k ,τ yk −1, k ：第 1 番目セルと第 k 番目セルの間に働く
x，y 方向の剪断力ｱ) 自由水面上に働く剪断力（風による粘性応力）

τ x0.1 = ρ aγ a 2Wx Wx 2 + Wy 2
τ y0.1 = ργ a 2Wy Wx 2 + Wy 280ρ a：空気密度
γ a2：海面摩擦係数
Wx，Wy：風速ベクトルの x，y 方向成分
ｲ) 海底面に働く剪断力

τ xK ,K +1 = ργ b 2uK uK 2 + vK 2
τ yK ,K +1 = ργ b 2vK uK 2 + vK 2

γb2：底面摩擦係数ｳ) レベル間に働く剪断力τ xk −1,k = ργ i 2 (u k −1 − uk ) ⋅(uk −1 − uk )2 + (vk −1 − vk )2

τ yk −1,k = ργ i 2 (vk −1 − vk ) ⋅(uk −1 − uk )2 + (vk −1 − vk )2(k = 2,3LL , K )γi2：内部摩擦係数
Qs：大気への放熱係数
Cv：比熱

81

y

z

ζ
k=

1
2

v
u

k =1

3
k=
2

x
Δz = h1wvuz = − H1
u

v

k=2
k=

u
v

Δz = h2w

52

z = −H 2
Δy
Δx

k=K−

1
2

k=K

k=K+

w

v

u

u
v

Δz = hk

12

図- 3.2.2.1

z = −H k

マルチレベルモデルの座標系と変数の配置

823) 拡散の基礎方程式
3 次元の物質拡散の基礎方程式は、
∂c
∂c
∂c
∂c
+u
+v
+w
∂t
∂x
∂y
∂z
∂c ⎞
∂c ⎞ ∂ ⎛∂⎛∂c ⎞ ∂ ⎛⎟⎟ + ⎜ Dv= ⎜ D hx
⎟ + Sc⎟ + ⎜⎜ Dhyx⎝∂z ⎠∂y ⎠ ∂z ⎝∂x ⎠ ∂y ⎝（式- 3.2.2.17）

である。ここで、Sc：内部変化項※大腸菌群数の計算ではゼロとする

c：物質濃度（g/m3）h：水深（m）

u，v，w ：x，y，z 方向の平均流速成分（m/s）
Dh，Dvt：水平、鉛直渦動拡散係数（m2/s）：時間（s）

・境界条件は水底（ z = − h ）で、
− Dv

∂c
= Rb
∂z

（式- 3.2.2.18）

ここで、Rb：底面からの溶出量（g/m2/s）境界層（ z = − hi ）で、
Dv

ci − ci +1∂cz = − hi = Dv
hi∂z

（式- 3.2.2.19）

水面（ z = ζ ）で、
Dv

∂c
=0
∂z

（式- 3.2.2.20）式-3.2.2.17 をセル内で Z 方向に積分すると次式になる。
∂(hk Ck ) = − ∂ (M k Ck ) − ∂ (N k Ck )∂t∂x∂y
− (wC ) z = − H k −1 + (wC ) z = − H k
∂Ck ⎞
∂ ⎛
∂C ⎞ ∂ ⎛⎟+ ⎜ hk K x 1 ⎟ + ⎜⎜ hk K y
∂y ⎟⎠
∂x ⎝
∂x ⎠ ∂y ⎝⎛ ∂C ⎞⎛ ∂C ⎞+ ⎜ Kz⎟ z = − H k + S c ⋅ hk
⎟ z = − H k −1 − ⎜ K z
∂z ⎠
∂z ⎠
⎝
⎝

83

（式- 3.2.2.21）

(3) 生態系（富栄養化）モデルの考え方
本モデルで取り扱う水質項目は表- 3.2.2.2 に示す 8 項目である。生態系モデルの概念図
は、図- 3.2.2.2のとおりである。
表- 3.2.2.2
項目

水質項目及びその主成分、単位
変 数
単
主 成 分
名

位

植物プランクトン（Phytoplankton）

クロロフィル a（Chlorophyll-a）Pμg − chla / l

動物プランクトン（Zooplankton）

動物炭素量―
Zooplankton Carbon

Z

mgC / lAmmonia Nitrogen、
無機態窒素
Nitrite Nitrogen、（ IN ―Inorganic Nitrogen）
Nitrate Nitrogen

CIN

mgN / l非生物（命）体
有機態窒素（ ON ―Organic Nitrogen） Organic NitrogenCON

mgN / l

無機態リン
IP
―
（
Phosphorus）
有機態リン
OP
―
（
Phosphorus）

Inorganic

リン酸態リン―
Orthophosphate

CIP

mgP / lOrganic

非生物（命）体
Organic Phosphorus

COP

mgP / l

化学的酸素要求量（ COD ）

非生物（命）体有機物
Chemical
OxygenDemandCCOD

mgO2 / l

溶存酸素（ DO ）Dissolved Oxygen

CDO

mgO2 / l

有機態の窒素及びリンの定式化には、生物（命）体を含める場合と含めない場合がある
が、本モデルでは、プランクトンの死骸や排泄物などを非生物体の有機物とし、生物体の植物プランクトン・動物プランクトンと分離して計算する方法をとる。すなわち本モデル
における有機態窒素、有機態リン及び COD の算出は、まず河川などから流入するのが非生物体であるとして、デトリタスなどによって生産される非生物体物質を微分方程式によって計算する。次に全有機態窒素、全有機態リン及び全 COD は、この非生物体と別途計算される植物プランクトン、動物プランクトンによって生産される生物体との合計で構成されるものとする。このようなモデル化により、有機態窒素、有機態リン及び有機物 COD が生
物体と非生物体の両方を含むことによる物質生産を、合理的に計算することが可能である。
何が有機物の生産や分解に関与するかが明確になる。
したがって、全有機態窒素 TON（Total Organic Nitrogen）は植物プランクトン及び動物プ
ランクトンの生物体有機態窒素と非生物体有機態窒素との加算となり、全有機態窒素 TONと無機態窒素との累加は全窒素 TN（Total Nitrogen）となる。同じく、全有機態リン TOP（Total
Organic Phosphorus）は植物プランクトン及び動物プランクトンの生物体有機態リンと非84

生物体有機態リンとを加算したものとなり、全有機態リン TOP と無機態リンとの累加は全リン TP（Total Phosphorus）となる。また、全 COD（ TCOD ―Total Chemical Oxygen Demand）
は植物プランクトン及び動物プランクトンの生物体 COD と非生物体の有機物 COD との累
加となる。

鉛直方向

日射量

：シミュレーション対象項目

移流
拡散（表層の場合なし）

再曝気（表層のみ）
水温

分
解
内部生産

COD

移流

無機態リン

摂取

光合成

植物プランクトン

流下方向

溶 沈
出 降

捕
食

呼吸
死滅

非生命体
有機態リン

拡散

呼吸
死滅

呼吸
死滅

動物プランクトン

有機態リン
溶 沈
出 降

無機態窒素

消費（底泥による）（底層のみ）

移流

排泄

排泄
分
解

摂取

消費（有機物分解に伴う）

D O

呼吸
死滅

分
解

濁質

流下方向

非生命体
有機態窒素

拡散

有機態窒素

沈
降

沈
降

沈
降

溶
出

沈
降

消
費

移流
拡散
（底層の場合なし）
底泥

鉛直方向

図- 3.2.2.2

底泥生態系モデルの概念図

生態系モデルによる内部領域における各水質の発生及び消滅濃度率は、濃度収支式（式-3.2.2.17）の右辺 S c 項となる。以下、図- 3.2.2.2で示したモデルの概念図に基づき、各水質項目ごとの発生及び消滅濃度率の算定法を示す。

85

1) 植物プランクトン（Phytoplankton）
植物プランクトンはクロロフィル a で表示される。これは、植物プランクトン個体数を計算する方法に比べ測定者による誤差が少なく、また、有機物生産などを考える場合には光合成の直接の担い手であるクロロフィル a がよい指標となりうるからである。
細胞が無限に増殖しない理由は多く存在するが、一般の環境では、基質や栄養塩の枯渇によるものである。ダム湖、湖沼や湾のように閉鎖性が強く、水が長期に滞留する水域では、植物プランクトンが、律速栄養塩として窒素やリンなどを摂取して、光合成により水中の無機炭素より有機物を合成し増殖する。水域で窒素及びリンの濃度が高くなり植物プ
ランクトンが異常増殖し、これに伴って種々の水質障害が生じる現象は富栄養化とよばれ、水質保全上の重要な問題となっている。
植物プランクトンの細胞（又は濃度）の増加速度 dP dt は、単位濃度当たりの増殖速度である比増殖速度 μ P を用いて表示される。
dP= μP Pdt

ここで、 μ P は植物プランクトンの比増殖速度、 P は植物プランクトン濃度 ( μg − chla / l)である。植物プランクトンの比増殖速度に影響する因子としては、日照強度、水温、無機
態窒素濃度、無機態リン濃度及びスペース効果を考慮する。
μ P = μ P max ⋅ FT ⋅ FI ⋅ FN
ここで、 μ P max はその最大比増殖速度、 FI 、 FN 、 FT はそれぞれ日照強度、栄養塩濃度
及び温度に関する影響関数である。
光の影響関数 FI としては照度（あるいは日射量） I に関して直線、飽和型又は指数型の近似法のいずれかがよく用いられる。ここで、日照強度効果は植物プランクトン増殖に対する強光阻害を考え、Ditoro らの考えをもとに以下の式（強光阻害も連続表示）で与える。
FI =

Iz
I
exp(1 − z )ISIS

ここで、 F1 は日照強度とスペース効果による影響の修正係数、 I S は最大増殖速度を与える最適照度（ lux 、あるいは最適日射量）、 I z は水深 z 位置における照度（ lux 、あるいは日射量）である。最適照度あるいは最適日射量 I S は種によって異なるが、2000～20000 lux あ
るいは 300～350 cal /(cm 2 ⋅ day ) 程度の値が日本ではよく用いられる。
光は水中で Lambert-Beer（ランベルト‐ベール）の法則に従って減衰する。したがって、水深 z 位置における照度 I z は
I z = I 0 exp(−λ (ζ − z ))λ = λ g P + λh
ここで、 I 0 は水面における照度（ lux 、あるいは日射量）、 ζ は水位、 λ は光の減衰係数であり、植物プランクトンによるもの（係数 λg）と、その他のもの（ λ h ）とに分ける必要がある。
86

栄養塩類濃度による影響 FN は律速栄養塩である無機態窒素及び無機態リンの濃度 C IN
及び C IP に関する Michaelis－Menten 型の式で近似する。
生物の生育収量は、必須栄養分の中で不足する養分量に支配されるため、最小値が選ばれてFN = MIN (

C IN
CIP
,
)、
K IN + CIN K IP + CIP

又は複数乗じられてFN =（式- 3.2.2.22）

C IN
C IP
K IN + C IN K IP + C IP

（式- 3.2.2.23）の二つの計算法がある。ここで、 K IN 及び K IP は半飽和定数である。
本モデルでは式-3.2.2.22 を採用する。
したがって、

μ P = μ P max ⋅ FT ⋅ FI ⋅ MIN (
又は

C IN
C IP
,
)、
K IN + C IN K IP + C IP

（式- 3.2.2.24）

C IN
C IP
K IN + C IN K IP + C IP

（式- 3.2.2.25）

μ P = μ P max ⋅ FT ⋅ FI ⋅がある。本モデルでは式-3.2.2.24 を採用する。
温度条件は微生物増殖・分解や物質変換経路にとって重要である。微生物による変換反応に関与する酵素は蛋白質であり、温度が高いほどその活性は高い一方で蛋白質の変質が生じるため温度の上限があり、最適温度範囲が存在する。水温による影響 FT は以下の式を
適用する。
FT =

T
T
exp(1 − )Ts
Ts

（式- 3.2.2.26）

ここで、θ （＝1.02～1.06）：定数、 T ：水温、 Ts ：最適温度、 Ts ： μ P max を決定したときの基準温度、一般的には Ts = 20 o C 。
ここに、

T ：

μP ：

μ P max ：水温(℃)
植物プランクトン増殖速度 (1 / day ) ；20 ℃ に お け る 植 物 プ ラ ン ク ト ン 最 大 増 殖 速 度(1 / day ) ；C IN ：

無機態窒素濃度 (mgN / l) ；C IP ：

無機態リン濃度 (mgP / l ) ；K IN ：

無機態窒素のMichaelis定数 (mgN / l) ；K IP ：

無機態リンのMichaelis定数 (mgP / l ) ；FT ：水温による影響の修正項；

Ｔs ：最適水温（℃）；

である。

87

植物プランクトンの変化速度は増殖速度と死滅速度とを加算した形で示される。したがって、植物プランクトンの増殖活動は一次反応型として次式で表される。
S P = μ P P − DP

ここで、 D P は植物プランクトンの減少である。植物プランクトンの減少項としては、自己分解項、動物プランクトンによる捕食項及び沈降項を考慮する。動物プランクトンはフィルターフィーダーであり、捕食速度については、動物プランクトン濃度と捕食率（濾水速度 Cg）に比例するものとし、飽食効果を併せて考慮した。沈降項については現存量に比例して沈殿するものと考えた。
D P = ( R Pθ PT − 20 + C g Z ) P − ω P∂P∂z

S P = μ P P − ( R Pθ PT − 20 + C g Z ) P + ω P（式- 3.2.2.27）∂P∂z（式- 3.2.2.28）

動物プランクトンの濾水速度 Cg は温度が上昇するとともに大きくなり、また植物プラン
クトン濃度が高くなるにつれて減少する。したがって、
T −20
C g = C g maxθ Pg

K Pg（式- 3.2.2.29）

K Pg + P

ここで、

D P ： 植物プランクトン減少速度 (1 / day ) ；R P ： 20℃における植物プランクトン死滅率 (1 / day ) ；θ P ： 植物プランクトン死滅率の温度補正係数；
C g ： 動物プランクトンの濾過（捕食）率 (l / mgC.day ) ；C g max ： 20℃における最大濾過率；

θ Pg ： 動物プランクトンの濾過率の温度補正係数；
Z ： 動物プランクトン濃度（ mgC / l ）；K Pg ： 飽食効果に対するMichaelis定数 ( μg − chla / l) 、動物プラ
ンクトンによる植物プランクトンの捕食の半飽和定数；

ω P ： 動物プランクトンの沈降速度 (m / day )
である。

88

2) 動物プランクトン（Zooplankton）
動物プランクトン（濃度 Z ）は動物炭素量で表示する。増殖項 μ Z については、植物プ
ランクトン捕食によるもののみを考える。捕食では前述のとおり、捕食速度が動物プラン
クトン濃度と濾水速度に比例するものとして、飽食効果も考慮する。減少項としては自己分解項及び自然死滅を考慮し、大型動物による摂取は考慮していない。また、自泳性を有することから沈降項も考慮しない。
したがって、動物プランクトンの変化速度は次式で表される。

S Z = ( μ z − R zθ zT − 20 − D z ) Z（式- 3.2.2.30）

μ z = γ PZ a z C g P（式- 3.2.2.31）

ここで、

μ z ： 動物プランクトン生成速度 (1 / day ) ；γ PZ ： C / chla 比 (mgC / μg − chla) 、植物プランクトン量からの動物プランクトン炭素量への換算係数；

a z ： 植物プランクトン同化率（摂食効率、0.40～0.85）；
R z ： 動物プランクトンの呼吸による分解速度定数 (1 / day ) ；θ z ： 動物プランクトンの呼吸による分解速度の温度補正係数
D z ： 動物プランクトンの自然死滅速度 (1 / day ) ；
である。

89

3) 無機態窒素（ IN ：Inorganic Nitrogen）
無機態窒素濃度は河川などからの流入量の他、増加項として有機窒素の分解による増加、底泥からの溶出（最下層コントロールボリュームのみ）による増加、さらに植物プランク
トンの自己分解による回帰、動物プランクトンの自己分解による回帰があると考えられる。
減少項としては、植物プランクトン増殖に伴う摂取、吸収による減少を考慮する。したがって、無機態窒素の変化速度は次式で表される。

S IN = －（植物プランクトンによる摂取）＋（植物プランクトンの自己分解による回帰）＋（動物プランクトンの自己分解による回帰）＋（非生物体有機態
窒素のバクテリアによる分解）＋（底泥からの溶出）= γ PN (− μ P + aPI RPθ PT − 20 ) P + γ ZN aZI RZθ ZT − 20 Z
T − 20T − 20CON + WINθ IN+ k Nθ ON

A
V（式- 3.2.2.32）

ここで、

C IN ： 無機態窒素濃度 (mgN / l) ；C ON ： 非生物体有機態窒素濃度 ( mgN / l) ；γ PN ： 植物プランクトン中の窒素クロロフィルa ( N / chla) 比(mgN / μg − chla) 、植物プランクトン量からの窒素量への換算係数；

a PI ： 植物プランクトンの回帰率（ ≈ 0.5 ）；k N ： 有機態窒素の無機化分解速度（ 1 day ）；θ ON ： 窒素分解速度の温度補正係数；
γ ZN ： 動 物 プ ラ ン ク ト ン 中 の 窒 素 炭 素 量 ( N / C ) 比(mgN / mgC ) ；a ZI ： 動物プランクトンの回帰率（ ≈ 0.5 ）；W IN ： 底泥からの無機態窒素溶出速度 ( gN / m 2 ⋅ day ) ；θ IN ： 底泥からの無機態窒素溶出速度の温度補正係数；2
A ： 最下層コントロールボリュームに対する底面積 (m ) ；3
V ： 最下層コントロールボリューム (m )

である。

90

4) 有機態窒素（ ON ：Organic Nitrogen）
本モデルでは、非生物体有機態窒素を取り扱う。非生物体有機態窒素濃度の変化は流入量及び放流量などからの影響のほかに、増加項として植物プランクトン及び動物プランク
トンの自己分解による回帰、動物プランクトンからの排泄、動物プランクトンの自然死滅であり、減少量としては、バクテリアによる無機態窒素への交換、非溶存態の非生物体有機態窒素の沈殿であると考えられる。なお、底泥からの溶出については有機態窒素としての溶出はないものとする。したがって、有機態窒素の変化速度は次式で表される。

S ON = （植物プランクトンの自己分解による回帰）＋（動物プランクトンからの排泄）＋（動物プランクトンの自己分解による回帰）＋（動物プランクトンの自然死滅）－（非生物体有機態窒素のバクテリアによる分解）＋（沈降）= γ PN ((1 − a PI ) RPθ PT − 20 + (1 − a z )C g Z ) P
T − 20+ γ ZN ((1 − a ZI ) RZ θ ZT − 20 + DZ ) Z − k N θ ON
C ON

+ (1 − f ON )ω ON（式- 3.2.2.33）

∂C ON
∂z

ここで、 ω ON は有機態窒素沈降速度 (m / day ) 、 f ON は溶存態非生物体の有機態窒素成分
である。
以上から、全有機態窒素および全窒素の濃度は下式により与えられる。
全有機態窒素：TON =（植物プランクトン態窒素）＋（動物プランクトン態窒素）+（非生物体有機態窒素）= γ PN P + γ ZN Z + C ON
全窒素：

TN = （植物プランクトン態窒素）＋（動物プランクトン態窒素）+（非生物体有機態窒素）＋（無機態窒素）= γ PN P + γ ZN Z + C ON + C IN915) 無機態リン（IP：Inorganic Phosphorus）
無機態リン濃度変化については河川からの流入量及び堰から放流量などの影響のほか、有機態リンのバクテリアによる無機態リンへの変換量、底泥からの溶出量、さらに動物プ
ランクトンからの排泄などにより増加する。また減少量については、植物プランクトンの摂取及び、懸濁物に吸着され、沈殿による減少を考える。したがって､無機態リンの変化速度は次式で表される。

S IP = －（植物プランクトンによる摂取）＋（植物プランクトンの自己分解による回帰）＋（動物プランクトンの自己分解による回帰）＋（非生物体有機態リンのバクテリアによる分解）＋（底泥からの溶出）＋（沈降）= γ PP (− μ P + a PI RPθ PT − 20 ) P + γ ZP a ZI RZ θ ZT − 20 Z
T − 20T − 20C OP + WIPθ IP+ k Pθ OP∂C IP
A+ (1 − f IP )ω IP
V∂z（式- 3.2.2.34）

ここで、

C IP ： 無機態リン濃度 (mgP / l) ；C OP ： 有機態リン濃度 (mgP / l) ；k P ： 有機態リンの無機化分解速度（20℃） (1 / day) ；θ OP ： 有機態リンの無機化速度の温度補正係数；

W IP ： 底泥からの無機態リン溶出速度 ( gP / m 2 ⋅ day ) ；θ IP ： 底泥からの無機態リン溶出速度の温度補正係数；
ω IP ： 無機態リンの沈降速度 (m / day )f IP ： 溶存態又は吸着されない無機態リン成分
である。
無機態リンが懸濁物に吸着されやすいため、本モデルでは吸着沈殿過程による無機態リンの減少も考慮する。

92

6) 有機態リン（ OP ：Organic Phosphorus）
本モデルでは、非生物体有機態リンを取り扱う。非生物体の有機態リン濃度の変化については河川からの流入・堰及び水門からの放流などの影響のほか、増加項として植物プラ
ンクトン及び動物プランクトンの自己分解による回帰、動物プランクトンからの排泄、動物プランクトンの自然死滅であると考えられる。減少項については、有機態窒素と同様にバクテリアによる無機態リンへの変換量、非溶存態の非生物体有機態リンの沈殿を考慮する。なお、底泥からの溶出については有機態リンとしての溶出はないものとする。したがって、有機態リンの変化速度は次式で表される。

S OP = （植物プランクトンの自己分解による回帰）＋（動物プランクトンからの排泄）＋（動物プランクトンの自己分解による回帰）＋（動物プランクトンの自然死滅）－（非生物体有機態リンのバクテリアによる分解）＋（沈降）= γ PP ((1 − a PI ) RPθ PT − 20 + (1 − a z )C g Z ) P
T − 20COP+ γ ZP ((1 − a ZI ) RZ θ ZT − 20 + DZ ) Z − k Pθ OP+ (1 − f OP )ωOP

（式- 3.2.2.35）

∂COP
∂z

ここで、 ω OP は有機態リン沈降速度 (m / day ) 、 f OP は溶存態非生物体有機態リン成分である。
以上から、全有機態リンおよび全リンの濃度は下式により与えられる。
全有機態リン：TOP =（植物プランクトン態リン）＋（動物プランクトン態リン）＋（非生物体有機態リン）= γ PP P + γ ZP Z + C OP
全リン：

TP = （植物プランクトン態リン）＋（動物プランクトン態リン）＋（非生物体有機態リン）＋（無機態リン）= γ PP P + γ ZP Z + C OP + C IP

937) 化学的酸素要求量（ COD ：Chemical Oxygen Demand）
有機物の化学的酸素要求量 COD は水中の易分解性有機物の存在量を示す指標であり、生物では植物、動物のプランクトン、非生物では懸濁態有機物や溶存態有機物のコンパー
トメントの挙動を知れば COD の時間変化を知ることができる。ここでは非生物体有機物の COD の変化速度を次式で表す。

S COD = （植物プランクトンの自己分解による回帰）＋（動物プランクトンからの排泄）＋（動物プランクトンの自己分解による回帰）＋（動物プランクトンの自然死滅）－（非生物体有機物のバクテリアによる分解）＋（底泥からの溶出）＋（沈降）= γ PC ((1 − a PI ) RPθ PT − 20 + (1 − a z )C g Z ) P
T − 20+ γ ZC ((1 − a ZI ) RZ θ ZT − 20 + DZ ) Z − k Cθ CT − 20 C COD + WCODθ COD

+ ωC

A
V

（式- 3.2.2.36）∂C COD∂z

ここで、

C COD ： 化学的酸素要求量 (mgO2 / l) ；γ PC ： 植 物 プ ラ ン ク ト ン の COD へ の 換 算 係 数(mgO2 / μg − chla ) ；γ ZC ： 動物プランクトンの COD への換算係数 (mgO2 / mgC ) ；k C ： 有機物分解速度定数 (1 / day ) ；θ C ： 有機物分解速度定数の温度補正係数；

WCOD ： 底泥からの COD 溶出速度 ( gO2 / m 2 .day ) ；θ COD ： 底泥からの COD 溶出速度の温度補正係数；
ω COD ： COD の沈降速度 (m / day )である。
以上から、全 COD の濃度は下式により与えられる。
全 COD ：
TCOD = （植物プランクトン態 COD ）＋（動物プランクトン態 COD ）+（非生物体有機物 COD ）= γ PC P + γ ZC Z + C COD

948) 溶存酸素（ DO ：Dissolved Oxygen）
溶存酸素 DO の変化は、水塊にすでに含まれている溶存酸素と飽和酸素との差を推進力とする水面からの再曝気、植物プランクトンの光合成作用に伴う溶存酸素の生成、また、有機物分解に伴う消費、底泥における消費などより表すことができる。
植物プランクトンによる光合成作用が活発な場合には、過飽和となる状況がみられ、また、水塊中に有機物が多量に含まれている場合には、低濃度となる富栄養化の影響側面を考える上で重要な水質項目である。DO の変化速度は下式により与えられる。

S DO = （植物プランクトンによる生産）－（植物プランクトン呼吸）－（動物プランクトン呼吸）－（非生物体有機物のバクテリアによる分解）＋（再曝気）－（底泥による消費）= γ PO ( μ P − RPθ PT − 20 ) P − γ ZO RZ θ ZT − 20 Z − k Cθ CT − 20 C COD
T − 20+ k Osθ Os(C DOs − C DO )As
T − 20 Ab− SODθ SOD
Vs
Vb（式- 3.2.2.37）

ここで、

C DO ： 溶存酸素濃度 ( mgO 2 / l ) ；γ PO ： 光合成による溶存酸素生産速度又は植物プランクトンのDO への換算係数 (mgO2 / μg − chla ) ；γ ZO ： 動物プランクトンの DO への換算係数 (mgO2 / mgC ) ；SOD ： Sediment Oxygen Demand―底泥による溶存酸素消費速度 ( gO 2 / m 2 ⋅ day ) ；底面のみ

θ SOD ： 底泥による溶存酸素消費速度温度補正係数；k Os ： 再曝気係数（m/day）；水面のみ

θ Os ： 再曝気係数温度補正係数；
C DOs ： 飽和溶存酸素濃度 (mgO2 / l) ；C DOs = (10.291 − 0.2809T + 0.006009T 2 − 0.000063T 3 ) × 32.0 22.4

As , V s ： 表層elementの表面積・体積；
Ab , Vb ： 底層elementの表面積・体積；
である。

95

(4) 干潟の浄化機能について
1) 干潟・浅場の機能
干潟の機能として、生息生育基盤としての機能、水質浄化機能および親水機能等があげられる。この内、水質浄化機能には次のものがあり、これらの機能による沿岸域の水質浄化効果を評価することが一般的な技術的課題となっている。
¾

海水中の懸濁物質の濾過（物理的作用）¾

曝気作用（物理的作用）¾

底生動物等による栄養塩、有機物の摂取（生物的作用）¾

底生動物→魚類・鳥類の捕食→系外への移動（生物的作用）¾

微生物による分解、無機化、脱窒（生物的作用）

2) モデルの要件
海水中の懸濁物質の濾過（物理的作用）
干潟（潮間帯）の微地形や潮位変化に伴う干出域の変化を表現するような空間分解能が必要である。干出⇔冠水の変化の計算技術も必要である。
曝気作用（物理的作用）
干潟、浅場の流れ場と波浪場を表現するモデルが必要である。また、貧酸素水の湧昇に対する曝気効果を評価するには密度成層を表現するような鉛直方向の空間分解能や計算技
術（k-ε、CIP 等）を導入する必要がある。
底生動物等による栄養塩等の摂取（生物的作用）
モデル構造としては、海底に単位面積当たりの摂取量（mg/day/m2）を与えれば、水中の栄養塩・有機物の収支は計算できる。（事例はあるものの摂取量に関する知見は少ない。）
底生動物→魚類、鳥類の捕食→系外への移動（生物的作用）
基本的に食物連鎖による物質収支を表現するモデルが必要であると考えられるが、その計算手法は確立されていない。
微生物による分解、無機化、脱窒（生物的作用）
分解、無機化及び脱窒は水温関数で与えるのが通常の方法であるが、干潟の評価のためには、干潟・浅場と沖合での異なる分解・無機化・脱窒速度を与える必要がある。干潟と沖合で異なる速度を与えた事例は少ない。
以上のような作用の定量的評価法については、現時点では未だ研究途上にあるが、三番瀬について一定の評価事例 3)があるため、これを参考に、(3)の生態系モデルにより浄化効
果を組み込んだ。

96

図- 3.2.2.1

ボックス
項 目
面 積 (ha)①浅海域→系外
二枚貝の漁獲脱窒
鳥の採餌
堆積・不活性化
ノリの漁獲(合計)(tonN/年)②単位面積当たりのT-N浄化量(mgN/㎡/日)③補正したT-N
浄化量(mgN/㎡/日)④系外→浮遊系
O-N (tonN/年)
COD (ton/年)⑤単位面積当たりのCOD浄化量(mg/㎡/日) 〔注1〕

ボックス分割及び窒素収支の算出範囲

1

2

3

4

5

6

7

1～7

260

201

277

341

141

185

162

1567

0.89

2.27

9.90

13.23

1.22

6.72

2.41

36.6

1.83
1.59
14.39
0.00
18.70

20.45 110.30 138.80
3.47 14.05 18.47
3.71
7.61 32.20
0.55
3.70
1.57
30.45 145.56 204.27

44.37
1.80
3.60
0.00
50.99

48.12
9.48
7.91
1.81
74.04

35.52
3.34
8.67
0.24
50.18

399.4
52.2
78.4
7.9
574.5

19.7

41.5 144.0

164.1

99.1

109.6

84.9

100.4

47.7

69.5 130.0

145.6

97.4

100.1

81.6

100.4

69.0

27.0

120.2

81.4

47.5

107.7

7.5

462.0

335.3
353.3

131.2
178.8

584.2
577.8

395.6
317.8

230.9
448.7

523.4
775.1

36.5 2245.3
61.7 392.6(注) 1．④系外→浮遊系のO-N浄化量は浅海域がO-Nの吸収源(Sink)となっていることを示し、またCOD浄化量は以下の換算式によって算出したものである。
[COD浄化量]＝［O-N浄化量］×［COD/PON比＝4.86］

図- 3.2.2.3

干潟の浄化機能に関する文献 3)

参考文献
1)都市水文研究会：多摩ニュータウン流出試験地調査報告，p333，1986
2)楠田哲也：感潮河川における物質輸送と変換，（財）河川環境管理財団

河川水質勉強会

資料， p36，2000
3) 環境庁水質保全局：平成10年度

藻場・干潟等の環境保全機能定量評価基礎調査書，1999
97

報告

3.3
3.3.1

モデルの構造およびインターフェイスの構築
システム構築の基本方針

既存の水物質循環モデルのプログラムは、データを入力しシミュレートする一連作業に手間と経験を要する等の問題が存在する。この状況を鑑み、一般のユーザーが容易に操作できるように、ユーザーインターフェース（システム）を作成する。
ここでシステム作成にあたっての基本方針は、以下によるものとする。
シミュレーションシステム構築方針① 開発は、Esri 社が販売する ArcGIS9.0 を利用する。
② 動作環境は、PC を考える。また OS は、WindowsXP を考える。
③ システムには、以下に示す機能を持たせる。
・ 入力データ（タンク情報、落水線網、気象関係データなど）は、既存ファイルの選択と外部アプリケーション（テキストエディタ）による修正を可能とする。
・ 計算結果を左右する基礎パラメータについては、システム上からの入力を可能とする。
・ 計算の実行は、システム上から操作可能とする。
・ 計算に入力するデータセットを識別し、既往計算の再現が可能とする。
・ 計算結果の表示は、グラフ表示及び数値表表示を可能とする。
モデルを実行するにあたり必要なソフト・作業環境は下記の通りである。
ArcView9.0、VB ランタイム(ArcView9.0 をインストールするために、以下動作環境が必要)
基本ソフトウェア
Windows 2000 Professional、Windows NT 4.0、 Windows XP Home Edition、Windows XP Professional Edition、Windows 2000 Advanced Server、Windows Server
2003
サービスパック／パッチ
Windows 2000：SP1、SP2、SP3、SP4（オプション）
Windows NT 4.0：SP 6a
Windows XP：SP1、SP2（制限付きサポート*）プロセッサ Intel Pentium または Intel Xeon プロセッサ（32 bit）
CPU の速度 800 MHz（必須）

1.0 GHz（推奨）

メモリ/RAM 256 MB （必須）

512 MB 以上（推奨）ディスク容量 605MB NTFS
3.3.2

695MB FAT32

システム構成本システムは、以下の図- 3.3.2.1に示す 構成を考える。ユーザー側の操作は、計算に使用する基礎データの選択と、モデル定数の設定、及び計算結果の表示とし、残りの作業はすべてシステム側に持たせるものとする。

98

ユーザー

システム

データ

Inportシステムが使用するファイル

Convert

図- 3.3.2.1
3.3.3

システム

計算条
件入力

計算結果表示Output水循環
プログラム

計算結果表示Output物質流動プログラム

システムを使用した操作状況図

機能概要

今回構築するシミュレーションシステムは、以下の機能を有するものである。

・施策設定
人口や土地利用などのフレーム設定および流域での対策による境界条件データの設定を行う。以前の計算において作成した施策の境界条件については、政策を選択することにより条件設定できるものとする。フレーム設定の変更があった場合には、境界条件の平面分布を ArcView 上で表示する。
・

計算パラメータの設定計算に用いる水文データ、計算モデルの自然流出系パラメータおよび湾内モデルのパラメータを画面に表示し、変更可能な機能を持たせる。
・陸域モデル、湾内モデルの実行
入力データおよびパラメータの設定後、計算実行ボタン等により、自動的に陸域・湾内の水・物質流動計算を行う。計算終了後、結果の保存ができる機能も持たせる。
・

計算結果の表示
下記の形態で計算結果を表示する。

A) 陸域特定地点時系列グラフ（流量・濃度）表示
B)

陸域の河道メッシュ年平均濃度の面表示

C)

湾内特定地点時系列グラフ(濃度・透明度)表示

D) 湾内指定日の面表示（濃度）993.3.4システムのフォルダ構成

本シミュレーションシステムは，データセット毎のファイル管理を容易にするため，下記に示すディレクトリ構成で構築するものとする。

(1) フォルダ全体構成
フォルダ構成ファイル・フォルダ内容の概要%Home%¥東京湾.mxd

シミュレーションシステムの実行ファイル¥CaseDir..txt

計算時の参照・格納フォルダを規定するファイル¥BASEDATA 変更

メッシュ基本情報などを格納するフォルダ¥dataedit

施策変更結果変換プログラムを格納するフォルダ¥System

ArcView 関連のファイルを格納するフォルダ¥ORG

フレームや施策変更前のオリジナルセット¥Folder01

ユーザが指定する計算ケースフォルダ¥Folder02

ユーザが指定する計算ケースフォルダ¥Folder03

ユーザが指定する計算ケースフォルダ

・・・・・(2) ¥ORG や¥Folder**以下のフォルダ構成
フォルダ構成ファイル・フォルダ内容の概要¥Folder**

陸域モデル水循環計算¥Hsimulation¥bin

実行形式(***.exe, ***.bat)¥cntl

計算条件ファイル

¥data

気象データ・人工系データ¥result

計算結果データ¥src

ソースプログラム¥Psimulation

陸域モデル物質循環計算¥bin

実行形式(***.exe, ***.bat)¥cntl

計算条件ファイル

¥data

人工系データ¥result

計算結果データ¥src

ソースプログラム¥Tokyo_bay

湾内モデル計算¥bin

実行形式(***.exe, ***.bat)¥cntl

計算条件ファイル

¥data

陸域モデル計算結果・境界データ¥result

計算結果データ¥src

ソースプログラム¥BASEDATA

計算ｹｰｽ毎のﾒｯｼｭ基本情報などを格納するﾌｫﾙﾀﾞ¥Flame

計算ｹｰｽ毎のﾌﾚｰﾑﾃﾞｰﾀを格納するﾌｫﾙﾀﾞ¥気象ﾃﾞｰﾀ

計算ｹｰｽ毎の気象ﾃﾞｰﾀを格納するﾌｫﾙﾀﾞ100

(3) 陸域水循環モデルのフォルダ構成およびファイル構造など① フォルダ構成

Bin：プログラム実行形式、バッチファイルを格納

Cntl：計算条件ファイルを格納

Data：計算元データファイルを格納

Result：計算結果ファイルを格納

Src：計算プログラムを格納

101

② 使用ファイルの一覧
Unit No.
99
1
2
3
5
41
42
43
44
45
46
47
48
51
52
53
60
61
62
65
66
67
68
71
72
7
11
12
13
14
15
16
17
18
21
22
23
24
25
26
27
28
31
32
33
34
35
36
37
38
39
75
76
80
85
90
86
87
88

File I/O
ファイル格納ファイル
説明Read Write フ ォ ル ダ○
入出力ファイルリスト
Hfilelist.dat
CNTL○
タンク構成情報を記述したファイル
tran.dat
CNTL○
湖沼周辺メッシュと計算メッシュとの対応をとるファイル
tnlake.dat
CNTL○
シミュレーションを制御するコントロールデータ
contrl.dat
CNTL○
シミュレーションの基礎定数データファイル
tntank.dat
CNTL○
降水量データファイル
precip.dat
DATA○
平均気温データファイル
tmmean.dat
DATA○
最低気温データファイル
tmmin.dat
DATA○
湿度データファイルhumity.dat
DATA○
日射量データファイル
suntim.dat
DATA○
風速データファイル
DATA
windsp.dat○
ダム流入量データファイル
daminf.dat
DATA○
hamon.dat
DATA
HAMONモデル定数を指定するファイル○
人工系水循環シミュレーション基本定数
manwat.dat
DATA○
人工系取水量DATA
supwat.dat○
人工系水循環用水番号
DATA
supcod.dat○
農業用水取水量データ
INTAKE.dat
DATA○
都市用水取水量データ
WATER.dat
DATA○
下水排水量データ
SUWAGE.dat
DATA○
メッシュ番号毎属性データ
FORT_MESHDATA.dat
DATA○
生活用水データ
FORT_CITYDATA.dat
DATA
○
県別データ
FORT_PREFDATA.dat
DATA○
農業用水データ
FORT_IRRIDATA.dat
DATA
流域メッシュの落水線番号とモデル領域X,Y座標を記述したファイル○
tnline.dat
CNTL○
河道タンクがあるメッシュ番号とモデル領域X,Y座標を記述したファイル
tnflow.dat
CNTL○
タンクからの流出量データ（物質流動モデルへとインプットされる）
tnkflow.dat
RESULT○
シミュレーションに使用された降水量のリスト
RESULThypre.dat○
シミュレーションに使用された最低気温のリスト
RESULThytma.dat○
シミュレーションに使用された平均気温のリスト
hytmm.dat
RESULT○
シミュレーションに使用された湿度のリスト
hyhum.dat
RESULT○
シミュレーションに使用された日射量のリスト
hysun.dat
RESULT○
シミュレーションに使用された風速のリスト
RESULThywin.dat○
シミュレーションに使用されたダム流入量のリスト
hyflw.dat
RESULT○
hyevp.dat
RESULT
シミュレーションに使用された蒸発散量データ（時間計算のみ）○
level1.dat
RESULT
指定メッシュ1・日付の水位（表層・不飽和・地下・川・融雪・融雪＋雨）○
level2.dat
RESULT
指定メッシュ2・日付の水位（表層・不飽和・地下・川・融雪・融雪＋雨）○
指定メッシュ河道流量
rivflw.dat
RESULT○
指定メッシュ成分流量
flwtot.dat
RESULT○
指定メッシュ（804）からの流出量（使用していない）
flw804.dat
RESULT○
指定メッシュ（614）からの流出量（使用していない）
flw614.dat
RESULT○
指定メッシュ（698）からの流出量（使用していない）
flw698.dat
RESULT○
指定メッシュ河道貯留量
rivstr.dat
RESULT○
計算対象モデルの基礎情報リスト
parain.dat
RESULT○
タンクのパラメータ出力リスト
parout.dat
RESULT○
流出計算結果の出力リスト
liscon.dat
RESULT○
流出計算結果の出力リスト
listres.dat
RESULT○
タンクの流量リスト（タンク毎に時系列出力）
RESULT
nenpyq.dat○
河道タンクの流出量
riverq.dat
RESULT○（使用していない）
netflw.dat
RESULT○
表層・不飽和帯・地下水タンクの貯留高
tnstor.dat
RESULT○
湖沼メッシュへの流出量
lakflw.dat
RESULT○
指定タンクからの流出量
tnline.dat
RESULT○
指定タンクからの流出量（詳細出力）
tnflow.dat
RESULT○
テンポラリ
RESULT
tmp.dat○（使用していない）
sum.dat
RESULT○
テンポラリ
tmp1.dat
RESULT○
DATA
メッシュのXY座標とメッシュ番号対応
X_Y_No.prn○
DATA
大口農業取水の受益エリア
Z_Agr_Area.prn○
DATA
大口生活取水の受益エリア
Z_Cit_Area.prn※必要に応じファイル名に流域番号を付けているFile Name※102

③ 入力ファイルフォーマット
FILE－1：TRAN．DAT（機番 1）【1 行目】NTANK －タンクの数（流域内メッシュの格子点の数）（ I5）
NSTORE－流出計算結果を記憶させておく領域の大きさ (表層・不飽和・地下・河道 )（ 4I5）
NLTANK－流域モデルの数＝タンクの数－ダム湖のタンク数（ I5）【2 行目以降】(1 行が 1 メッシュに関する情報)(読み取りなし )－計算順序（ I5）

ITANK－タンクの番号（ I5）

IRSW－フラグ： ＝ 0，河道のないモデル＝ 1，表層タンクと不飽和帯タンクの流出量が河道に流入する。
＝ 2，表層タンクと不飽和帯タンクと地下水タンクの流出量が河道に流入する。
＝ 3， 2 と同じ流出量が河道に流入し特性曲線法により洪水追跡をする。
＝ 4，ダム湖の計算

IN1－表層タンクへの流入量の番地（ I3）IN2－不飽和帯タンクへの流入量の番地（ I3）

IN3－地下水タンクへの流入量の番地（ I3）IN4－河道タンクへの流入量の番地（ I3）

I01－表層タンクの流出量の番地（ I3）I02－不飽和帯タンクの流出量の番地（ I3）I03－地下水タンクの流出量の番地（ I3）

I04－河道タンクの流出量の番地（ I3）IA－流域面積＝数値×メッシュ面積／ 8（ I7）

IPRC－雨量観測所番号（ I3）

ITMP－気温の観測所番号（ I3）

IHUM－湿度の観測所番号（ I3）

ISSD－日照時間の観測所番号（ I3）

IWS－風速の観測所番号（ I3）

MDF－表層タンクの識別番号（ I2）

MDS－不飽和帯タンクの識別番号（ I2）

MDG－地下水タンクの識別番号（ I2）

MDR－河道タンクの識別番号（ I2）MSE－（使用していない）（ I2）

IDAMC－ダムの番号（ I2）
RSLOP －河道勾配（ F10.6）
SSLOP －斜面勾配（ F10.6）
ALPHA－係数，係数×メッシュの面積＝タンクの面積（ F10.6）
ALTIT －標高（ F10.6）
WDTH －河道幅（ F10.6）
IWET－降雨割増係数の識別番号（ I2）103

FILE－2：TNLAKE．DAT（機番 2）
内容：湖の周辺のメッシュに対応する計算区間の番号（東京湾では使用していない）

FILE－3：CONTRL．DAT（機番 3）【1 行目】IDE－【2 行目】IMW－【3 行目】JMW【4 行目】IGWAT【5 行目】KINEM

蒸発散モデル識別番号 ( I5)=1,

Penman モデル=2,

Hamon モデル

=3,

Morton モデル=4,

Brutsaert-Stricker モデル

人工系水循環シミュレーションフラグ ( I5)=0,

人工系水循環シミュレーションしない=0,

人工系水循環シミュレーション－

人工系水循環取水地点の最大数 ( I5)－

地下水タンクコントロール変量 ( I5)－特性曲線法コントロール変量 ( I5)

DTK－

標準計算時間、単位は分 ( I5)

NTK1

－

時間分割数 ( I5)

NTK2

－

解の安定のための時間分割数 ( I5)【6 行目】ITDEL－

日計算及び時間計算を指定するフラグ（ 0:日計算， 1:時間計算）

FILE－4：TNTANK．DAT（機番 5）【1 行目】
NENS－計算開始年（西暦）（ I5）MONS－計算開始月（ I5）

NDYS－計算開始日（ I5）

NENE－計算終了年（西暦）（ I5）MONE－計算終了月（ I5）

NDYE－計算終了日（ I5）

NENDS －雨量・気象データの開始の年（ I5）MONDS－雨量・気象データの開始月（ I5）

NDYDS－雨量・気象データの開始日（ I5）

NENDE －雨量・気象データの最後の年（ I5）MONDE－雨量・気象データの最後の月（ I5）

NDYDE－雨量・気象データの最後の日（ I5）【2 行目】
NDAM－ダムの数（ I5）MAXS－不飽和帯のタンクの鉛直方向の最大個数（ I5）

LDPRINT－フラグ (I5)：＝ 0，雨量・気象データを出力しない。
＝ 1，雨量・気象データを出力する。

ISWINP －フラグ (I5)： ＝ 0，初期水位を仮定値で与える。
＝ 1，初期水位を入力する。

ISWOUT－フラグ (I5)： ＝ 0，ある月日の水位を DISK に出力しない。
＝ 1，ある月日の水位を DISK に出力するが、 PRINTER には出力しな

104

い。
＝ 2，ある月日の水位を DISK に出力し、さらに PRINTER に出力する。

ISWICK －フラグ (I5)： ＝ 0，計算の進行状況を出力しない。
＝ 1，計算の進行状況を出力する。

ISPDAM －フラグ (I5)： ＝ 0，ダム調節計算の結果を出力しない。
＝ 1，ダム調節計算の結果を出力する。

【3 行目】NPRC－雨量のデータの数（ I5）NTMP－気温のデータの数（ I5）

NHUM－湿度のデータ数（ I5）

NSSD－日照時間のデータの数（ I5）NWS－風速のデータの数（ I5）NIF－流入量のデータの数（ I5）【4 行目】
IBAS－フラグ（ I5）

【5 行目】NMDF －表層タンクの定数の種類（ I5）【6 行目】
SKF－表層タンクの定数：最終浸透能 (cm/s)（ F10.0）

HFMXD－表層タンクの定数：表面の最大貯留高 (m)（ F10.0）

HFMND－表層タンクの定数：早い流出が発生する高さ ( m)（ F10.0）

HFOD－表層タンクの定数：地下浸透が発生する高さ ( m)（ F10.0）

SNF－表層タンクの定数：地表面の粗度係数（ m -1/3s-1）（ F10.0）

FALFX－表層タンクの定数：早い中間流の流出量を規定する係数（ F10.0）

HIFD－表層タンクの定数：計算のための初期値 ( mm)（ F10.0）【6+NMDF 行目】
NMDS－不飽和帯タンクの定数の種類（ I5）【7+NMDF 行目】
NSTAND－不飽和帯タンクの段数（ I5）【8+NMDF 行目】
SKD－不飽和帯タンクの定数：鉛直方向の透水係数 (cm/s)（ F10.0）【9+NMDF 行目】
SKX－不飽和帯タンクの定数：水平方向の透水係数 (cm/s)（ F10.0）【10+NMDF 行目】
HMXSD－不飽和帯タンクの定数：層厚 ( m)（ F10.0）【11+NMDF 行目】
STS－不飽和帯タンクの定数：飽和水分量を想定する（ F10.0）【12+NMDF 行目】
STW－不飽和帯タンクの定数：最小水分量を想定する（ F10.0）【13+NMDF 行目】
SBD－不飽和帯タンクの定数：ψ～θ関係を設定する定数（ F10.0）【14+NMDF 行目】
HISD－不飽和帯タンクの定数：計算のための初期値（水分量で与える）（ F10.0）【7+NMDF+NMDS×8 行目】
NMDG－地下水タンクの定数の種類（ I5）105【8+NMDF+NMDS×8 行目】
AUD－不圧地下水の流出係数 (1/mm/day)（ F10.0）

AGO－被圧地下水の流出係数 (1/day)（ F10.0）

HCGD－不圧地下水の流出する貯留高 ( mm)（ F10.0）

HIGD－計算のための初期値 ( mm)（ F10.0）【8+NMDF+NMDS×8+NMDG 行目】
NMDR－河道タンクの定数の種類（ I5）【9+NMDF+NMDS×8+NMDG 行目】
RBW－係数、実際の川幅から設定する（ F10.0）

RBS－係数、 0.3～ 0.5 程度である。（ F10.0）

RNS－ Manning の粗度係数（ m -1/3s-1）（ F10.0）

RRID－計算のための初期値 (m)（ F10.0）【10+NMDF+NMDS×8+NMDG 行目】
使用しない【9+NMDF+NMDS×8+NMDG+NMDR×2 行目】
ISDAM1－ダムの機能の番号（ I5）

ISDAM2－ダムの機能の番号（ I5）

DAMNAM－ダムの名前（ A20）

【9+NMDF+NMDS×8+NMDG+NMDR×2+NDAM 行目】
DAA－メッシュの面積 (m2)（ F10.0）【10+NMDF+NMDS×8+NMDG+NMDR×2+NDAM 行目】
RRAT－雨の割増係数（ F10.0）

TEMPGR－気温の増加量 (℃ )（ F10.0）
TDCOF－気温の減率（ ℃ ／ｍ）（ F10.0）

EVC1

2π

EVC2

EVC=EVC1+EVC2

Sin［ 365.25(T－ EVC3)］EVC3で計算する蒸発量の補正係数（ F10.0）【11+NMDF+NMDS×8+NMDG+NMDR×2+NDAM 行目】
NWET－降雪の割増係数の個数（ I5）【12+NMDF+NMDS×8+NMDG+NMDR×2+NDAM 行目】
TICOF－降雪の割増係数（ F10.0）

COMELT－降雪のおくれ係数（ F10.0）【12+NMDF+NMDS×8+NMDG+NMDR×2+NDAM+NWET 行目】
NPOINT－計算結果の詳細リストを出力するタンクの個数【13+NMDF+NMDS×8+NMDG+NMDR×2+NDAM+NWET 行目】
LP－出力するタンクの番号（ NPOINT 分）（ 6I5）106

【14+NMDF+NMDS×8+NMDG+NMDR×2+NDAM+NWET 行目】
NPOINT－河道タンクの流出量を年表で出力するタンクの個数【15+NMDF+NMDS×8+NMDG+NMDR×2+NDAM+NWET 行目】
LP－出力する河道タンクの番号（ NPOINT 分）（ 6I5）

【16+NMDF+NMDS×8+NMDG+NMDR×2+NDAM+NWET 行目】
HOSEI2－年表形式の出力データにかかる係数。 HOSEI2 倍のデータが出力される。（ 6F10.0）【17+NMDF+NMDS×8+NMDG+NMDR×2+NDAM+NWET 行目】
NPOINT－河道タンクの流出量を DISK に出力するタンクの個数【18+NMDF+NMDS×8+NMDG+NMDR×2+NDAM+NWET 行目】
LP－出力する河道タンクの番号を入力する。（ NPOINT 分）（ 6I5）

【19+NMDF+NMDS×8+NMDG+NMDR×2+NDAM+NWET 行目】
HOSEI3－ DISK 出力データにかかる係数。 HOSEI3 倍にデータが出力される。（ 6F10.0）【20+NMDF+NMDS×8+NMDG+NMDR×2+NDAM+NWET 行目】
NPOINT－タンクの定数を出力するタンクの個数【21+NMDF+NMDS×8+NMDG+NMDR×2+NDAM+NWET 行目】
LP－タンクの定数を出力するタンクの番号を入力する。（ NPOINT 分）（ 6I5）

107

(4) 陸域物質循環モデルのフォルダ構成およびファイル構造など① フォルダ構成

Bin：プログラム実行形式、バッチファイルを格納

Cntl：計算条件ファイルを格納

Data：計算元データファイルを格納

Result：計算結果ファイルを格納

Src：計算プログラムを格納② 使用ファイルの一覧
Unit No.
99
1
2
3
4
5
23
31
32
33
34
51
52
53
54
55
56
60
80
90
61
43
91
62

File Name
Pfilelist.dat
TRAN.DATLOADCOD.DAT
TNKFLOW.DATPARAMET.DATTNTANK.DAT
SUWAGE.DAT
FORT_MESHDATA.DAT
FORT_CITYDATA.DAT
FORT_PREFDATA.DAT
FORT_LOADDATA.DATPARAIN.DATPAROUT.DATFLOWLS.DATLOADLS.DATTABLEP.DAT
OutDiskp00.dat
PLTOT.DAT
tmp.dat
tmp1.dat
A_Seibun.dat
湾直接.csv
Syori2bay.dat
paverage.dat

File I/O
入出力箇所
フ
ァ
イ
ル
説
明
Read Write○
入出力ファイルリスト
CNTL○
タンク構成情報を記述したファイル
CNTL(水)○
人工系物質循環データ番号
DATA○
RESULT(水) シミュレーションの基礎定数データファイル○
ノンポイントソース系物質循環モデル定数
CNTL○
構築プログラムの出力ファイル
CNTL(水)○
下水排水量データ
DATA○
メッシュ番号毎属性データ
DATA○
生活用水データ
DATA
○
県別データ
DATA
○
人工系汚濁濃度データ
DATA○
入力した水循環モデルの定数
RESULT○
入力した物質循環モデルの定数
RESULT○
流出計算結果RESULT○
負荷量計算結果
RESULT○RESULT
年表形式で負荷量計算結果○
RESULTDISKに出力した負荷量○
流域出口においてタンクの成分出力
RESULT○
テンポラリ
RESULT○
テンポラリ
RESULT○
土地利用別流出負荷量出力ファイル
RESULT○
湾直接放流下水処理場判別ファイル
DATA○
湾モデル計算のための中間ファイル
RESULT○
メッシュ毎年平均水質データ
RESULT

モデル計算後湾モデル計算の為のプログラムで使用されるファイル
Temp_tokyo.txt
To－Baymodel.cnt
E_kasen_hokan.dat
CHECK.csv
v_kasen2001_day-c1.prn
e_kasen2001_day_**.prn

○
○
○
○
○
○

CNTL
CNTL
CNTL
RESULTRESULTRESULT108

東京(大手町)日平均気温
湾計算境界データ作成時の設定ファイル
湾計算境界データ時刻内挿時の中間ファイル

テンポラリ
湾モデル計算境界条件(河川流量)中間ファイル
湾モデル計算境界条件(河川水質)中間ファイル

(5) 東京湾モデルのフォルダ構成およびファイル構造など① フォルダ構成

Bin：プログラム実行形式、バッチファイルを格納

Cntl：計算条件ファイルを格納

Data：計算元データファイルを格納

Result：計算結果ファイルを格納

Src：計算プログラムを格納② 使用ファイルの一覧
入力データ１(¥cntl)
ファイル名
1-1 coefficient.txt
1-2 coefficient_area.dep

1-3
1-4
1-5
1-6

coefficient_area.txt
ml3eco.dep
ml3eco.drw
ml3eco.jcl

内 容
水質モデルのパラメータ
・水深データ形式の海底面の境界条件設定ファイル
・各計算格子毎に coefficient_area.txt で定義されるどの番号の値を設定するか定義されている。
溶出速度、DO 消費速度の設定値(単位：g/m2/day)
河川の流入位置、水深値等定義
モニターする点の指定
入力データや出力ファイルの指定

Unit 番号
(16)
(51)

(16)
(1)
(5)
(4)

入力データ 2(¥data)

2-1
2-2
2-3
2-4
2-5

ファイル名

内

容

v_kasen2001_mit.txt
e_kasen2001_mit.txt
e_kyokai2001_mit.txt
TIDE2001_mit.txt
k_kyokai2001_mit.txt

各河川の流量
各河川の水質・水温
境界の水温、塩分、水質
推算潮位データ
海域の気象条件(気温・北方，東方風速，日射量)109

Unit 番
号
(72)
(73)
(2)
(3)
(2)

出力データ(¥result)
3-1
3-2
3-3
3-4
3-5
3-6
3-7
3-8
3-9
3-10
3-11
3-12
3-13
3-14
3-15
3-16
3-17
3-18
3-19
3-20
3-21
3-22
3-23
3-24
3-25
3-26
3-27
3-28
3-29

ファイル名
choryu.cc1
choryu.cc2
choryu.cc3
choryu.cc4
choryu.cc5
choryu.cc6
choryu.cc7
choryu.cc8
choryu.cc9
choryu.cc0
choryu.den
choryu.flw
choryu.prt
choryu.sal
choryu.temp
choryu_cc1.csv
choryu_cc2.csv
choryu_cc3.csv
choryu_cc4.csv
choryu_cc5.csv
choryu_cc6.csv
choryu_cc7.csv
choryu_cc8.csv
choryu_cc9.csv
choryu_cc0.csv
choryu_flw.csv
choryuv_sal.csv
choryuv_temp.csv
monitor.csv

内 容
メッシュ毎植物プランクトン濃度計算結果(バイナリ)
メッシュ毎動物プランクトン濃度計算結果(バイナリ)
メッシュ毎 I-N 濃度計算結果(バイナリ)
メッシュ毎 O-N 計算結果(バイナリ)
メッシュ毎 I-P 濃度計算結果(バイナリ)
メッシュ毎 O-P 計算結果(バイナリ)
メッシュ毎 COD 濃度計算結果(バイナリ)
メッシュ毎 DO 計算結果(バイナリ)
メッシュ毎 T-N 濃度計算結果(バイナリ)
メッシュ毎 T-P 計算結果(バイナリ)
メッシュ毎密度計算結果(バイナリ)
メッシュ毎潮流計算結果(バイナリ)
計算結果リスト(テキスト)
メッシュ毎塩分濃度計算結果(バイナリ)
メッシュ毎水温計算結果(バイナリ)
メッシュ毎植物プランクトン濃度計算結果(CSV)
メッシュ毎動物プランクトン濃度計算結果(CSV)
メッシュ毎 I-N 濃度計算結果(CSV)
メッシュ毎 O-N 計算結果(CSV)
メッシュ毎 I-P 濃度計算結果(CSV)
メッシュ毎 O-P 計算結果(CSV)
メッシュ毎 COD 濃度計算結果(CSV)
メッシュ毎 DO 計算結果(CSV)
メッシュ毎 T-N 濃度計算結果(CSV)
メッシュ毎 T-P 計算結果(CSV)
メッシュ毎潮流計算結果(CSV)
メッシュ毎塩分濃度計算結果(CSV)
メッシュ毎水温計算結果(CSV)
特定地点計算結果リスト(CSV)110

Unit 番号
(21)
(22)
(23)
(24)
(25)
(26)
(27)
(28)
(29)
(30)
(14)
(11)
(7)
(13)
(12)
(41)
(42)
(43)
(44)
(45)
(46)
(47)
(48)
(49)
(50)
(31)
(33)
(32)
(35)

③ 入力ファイルフォーマット
ml3eco.dep】

【1-4
１行目
カラム

1～

11～

21～

31～

変数名
FORMAT
単位
内容

ＤＴ
F10.0
秒
潮流計算のタイムステップ

ＤＴＤ

ＤＸ

ＴＩＮＲ

Å

Å

Å

*1

秒
ｍ
拡 散 計 算 ( 内 格子間隔
部生産)のタイムステップ

秒
移流項のタイムステップ

DT､DTD は以下の C.F.L(Courant-Friedrichs-Lewy）条件を満足しなければならない。

DT <
DTD <

DX
2 gH max( DX 2 )2Kg: 重力加速度
Hmax : 領域の最大水深
K: 水温・塩分の拡散係数(K は次頁の THD･SHD のうち大きい値を用いる。)*2 移流項は他の圧力項・粘性項等に較べてオーダーが低いため毎タイムステップ計算しない。
2 行目
カラム

1～

11～

21～

31～

41～

51～

変数名
FORMAT
単位
内容

START
F10.0
時
計算開始
時刻

ENDD

FMAP

TMAP

FFLO

TFLO

Å

Å

Å

Å

Å

時
時
秒
時
計算終了 リスト出力 リスト出力 図 化 フ ァ
時刻
開始時刻 間隔
イル出力
開始時刻

秒
図化ファ
イル出力
間隔

3 行目
カラム

1～

21～

31～

変数名
FORMAT
単位
内容

STARTD
F10.0
時
水温・塩分の拡散計算開
始時刻

ＦＤＩＦ(未使用）

ＴＤＩＦ(未使用)

ÅÅ時
拡散プログラム(DIF3D) 用 フ ァ イ
ル出力開始時刻秒
拡散プログラム(DIF3D) 用 フ ァ イ
ル出力間隔*1 拡散計算用ファイルの出力間隔は、拡散計算に用いる拡散係数から C.F.L.を用いて決定す
る。

111

４行目
カラム

1～

11～

21～

変数名
FORMAT
単位
内容

GAV
F10.0
ｃｍ/Ｓ2
重力加速度

COR

GAM

Å

Å

度
緯度

無次元
海底摩擦係数

５行目
カラム

1～

11～

変数名
FORMAT
単位
内容

HVS
F10.0
ｃｍ2/Ｓ
水平渦動粘性係数

VVS
F10.0
ｃｍ2/Ｓ
鉛直摩擦係数*1 拡散計算用ファイルの出力間隔は、拡散計算に用いる拡散係数から C.F.L.を用いて決定す
る。
６行目
カラム

1～

11～

変数名
FORMAT
単位
内容

THD
F10.0
ｃｍ2/Ｓ
水温の水平拡散係数

TVD
F10.0
ｃｍ2/Ｓ
水温の鉛直拡散係数

カラム

1～

11～

変数名
FORMAT
単位
内容

SHD
F10.0
ｃｍ2/Ｓ
塩分の水平拡散係数

SVD
F10.0
ｃｍ2/Ｓ
塩分の鉛直拡散係数

７行目

8 行目
カラム

1～

変数名
FORMAT
単位
内容

KV
I5
層数

112

９行目
カラム

1～

11～

21～

31～

101～

変数名
FORMAT単位
内容

HH(1)10F10.0cm
層厚(1 層)

HH(2)

HH(3)

HH(4)

HH(10)

Å

Å

Å

Å

Å

Å

Å

Å

Å

Å

層厚(2 層)層厚(3 層)層厚(4 層)

*1

11～

21～

31～

101～

１０行目
カラム

1～

変数名
FORMAT
単位
内容

TINT(1)
TINT(2)
Å
10F10.0
Å
℃
初期水温(1 層)

TINT(3)

TINT(4)ÅÅÅ

TINT(10)
Å

Å

Å

Å

Å

*1

１１行目
カラム

1～

11～

変数名
FORMAT
単位
内容

SINT(1)
SINT(2)
Å
10F10.0Å
PSU
初期塩分(1 層)

21～

31～

SINT(3)

SINT(4)

Å
Å

Å
Å

101～
SINT(10)
Å
Å

Å
Å

*1

１２～１９行目（８行分）
カラム

1～

11～

21～

31～

71～

変数名
FORMAT単位
内容

C?INT(1)10F10.0
PSU
初期水質濃度(1 層)

C?INT(2)

C?INT(3)

C?INT(4)

Å

Å

Å

Å

Å

Å

Å

Å

Å

Å

C?INT(10)*1*1 C1INT～C8INT は植物プランクトン,動物プランクトン, 無機態窒素，有機態窒素，無機態リン，有機態リン，COD，DO の順*2 ＫＶ層分定義する。（ ＫＶ＝１（単層）の場合は最大水深より大きい値を入力する。）
▽
HH(1)
HH(2)
HH(3)
HH( )
HH(KV)

上層から１、２、３、･･･ＫＶ層と数える。
k=1
k=2
k=3
k=KV

5.00
１層目
2 層目

図- 3.3.4.1
113

3.00

初期濃度の設定

２０行目
カラム

1～

変数名
KAS
FORMAT I10
内容
開境界数
２１行目
カラム

1～

6～

11～

16～

変数名
FORMAT
内容

IAS
I5
始点開境
界のＸ座
標

JAS

IBS

JBS

Å

Å

Å

始点開境 終点界境 終点開境
界のＹ座 界のＸ座 界のＹ座標
標
標

*1 次ページのデータを入力する。（KAS=0 の場合は不用）*2 通常潮位境界と同じ位置となる。(17/22 参照)
層数分の行
カラム

1～

11～

21～

31～
SKB(1)
F10.0

変数名
FORMA
T
単位
内容

TKS(1)
F10.0

SKS(1)F10.0

TKB(1)
F10.0℃
始点境界水温PSU℃
PSU
始点境界塩分 終点境界水温 終点境界塩分

1* 計算層の数必要,計算全域に設定した層数を持たない境界も定義（0.0 を）する．
層数分の行×８
カラム

1～

11～

変数名
FORMA
T
単位
内容

C?KS(1)
F10.0

C?KB(1)F10.0℃
始点境界水質PSU
終点境界水質・・・・・・・・KV 層分×8 項目分 必要

1 層目

水温

計算領域

塩分

20.0

21.0

32.0

32.0

19.0

19.0

33.0

33.0

2 層目
J
1 層目
(2 ,2 )

水温
( 3 0 ,2 )

2 層目
I

塩分

境界条件

図- 3.3.4.2
114

開境界条件の設定

２１＋（KV*９）＋１行目
カラム

1～

変数名
FORMAT
内容

UNH
F10.0
水深の出力パラメータ

２１＋（KV*９）＋２行目
カラム

1～

11～

21～

変数名
FORMAT
内容

UNZ
F10.0
水位の出力
パラメータ

UNUV

UNW

Å

Å

水平流速の出 鉛直流速の出力
力パラメータパラメータ

２１＋（KV*９）＋３行目
カラム

1～

11～

21～

変数名
FORMAT
内容

UNSAL
F10.0
塩分の出力
パラメータ

UNTMP

UNDEN

Å

Å

水温の出力パ 密度の出力パラ
ラメータ
メータ

２１＋（KV*９）＋４行目
カラム

1～

11～

21～

変数名
FORMAT
単位
内容

WDEN
F10.0

TAWX

TAWY

Å

Å

空気の密度

cm/S
cm/S
X 方向の風速成分 Y 方向の風速成分*1 風は吹き去る方向で扱うことに注意*2 北方向とメッシュを分割方向にも注意

北西の風 3 ｍ/秒の場合
2.21 ｍ/秒

3.00 ｍ/秒
2.21 ｍ/秒
図- 3.3.4.3

風速成分の設定方法

115

２１＋（KV*９）＋５行目（未使用）
カラム

1～

変数名
FORMAT
単位
内容

OMG（未使用）
F10.0
時
境界潮汐の周期

半日周潮の場合は 12.5 時間（45000 秒）
２１＋（KV*９）＋６行目（未使用）
カラム

1～

変数名
FORMAT単位
内容

KAK（未使用）I10
開境界のライン数*1 KAK は１にする事、境界のラインは今のところ１つしか設定できない。
２１＋（KV*９）＋７行目（未使用）
カラム

1～

6～

11～

16～

21～

26～

31～

36～

41～

46～

変数名
FORMAT
単位
内容

IA
I5

JA
I5

APA
F5.0
cm
(IA,J
A) の
潮 汐
振幅

PHA

ＳＣＡ

Å

Å

度
(IA,J
A) の
遅角

cm
(IA,J
A) の
平 均
水面

APB
F5.0
cm
(IB,JB
)の潮
汐 振
幅

ＳＣＢ

Å

JB
I5

PHB

Å

IB
I5

度
(IB,JB)の遅角cm(IB,JB
)の平
均 水
面

始 点
開 境
界 の
Ｘ 座
標

始 点
開 境
界 の
Ｙ 座
標

終 点
開 境
界 の
Ｘ 座
標

終 点
開 境
界 の
Ｙ 座
標

*1 IA,JA,IB,JB 以外は使用しない。境界の潮位は別のファイルでステップ毎のデータが入力さ
れる。

陸
陸

③

湾

海

島

②
陸

①

①

ライン数=１

ライン数＝3図- 3.3.4.4

開境界での潮位振幅の定義方法

116

２１＋（KV*９）＋８行目
カラム

1～

6～

*1 河川流入量が大きい場合、急激にそ

変数名
FORMAT
単位
内容

NX
I5

ＦＬＯＷＮ
F10.0
時
Ｘ 方 向 に 流 河川流量の増加時間
入河川数の値を境界に与えると計算が発散することがある。
そのために徐々に流入量を増加させる。（未定義の場合は 1 時間とする）

２１＋（KV*９）＋９行目～
カラム

1～

6～

11～

16～

21～

26～

変数名
FORMAT
単位
内容

ILNF
I5

IINF

IJNF

IKNF

IXFUGO

IXCH

Å

Å

Å

Å

流入地点
の Y 座標流 入 地 流入方向
点の K 層海 域 流入地点
段階の X 座標31～

UTR（未使用）
Å
F10.0
ｍ3/秒
流 入 流入量量 の
番号２１＋（KV*９）＋９行目～のつづき
カラム

41～

51～

変数名
FORMAT
単位
内容

UTB（未使用）
F10.0
℃
水温

USB（未使用）Å・・・・・・・・NX 回分必要

PSU
塩分

*1 流量の符号(IXFUGO)は流入・流出でなく座標系に対しての方向であることに注意*2 IXCH は河川流入のファイルのケース番号を指定する。
*3 UTR，UTB,USB はステップ毎のデータが与えられる。(ここでの設定は必要ない)
２１＋（KV*９）＋９＋NX＋１行目
カラム

1～

変数名
FORMAT
内容NYI10
Y 方向に流入河川数117

２１＋（KV*９）＋９＋NX＋２行目～
カラム

1～

6～

11～

16～

21～

26～

変数名
FORMAT
単位
内容

JLNF
I5

JINF

JJNF

JKNF

JYFUGO

JYCH

Å

Å

Å

Å

海 域 段 流 入 地 流入地点 流入地点 流入方向
階
点 の X の Y 座標 の K 層
座標

31～

VTR（未使用）
Å
F10.0
ｍ3/秒
流 入 流入量量 の
番号２１＋（KV*９）＋９＋NX＋２行目～のつづき
カラム

41～

51～

変数名
FORMAT単位
内容

VTB（未使用）
F10.0
℃
水温

VSB（未使用）Å・・・・・・・・NY 回分必要

PSU
塩分

*1 流量の符号(JYFUGO)は流入・流出でなく座標系に対しての方向であることに注意*2 JYCH は河川流入のファイルのケース番号を指定する。
*3 VTR，VTB,VSB はステップ毎のデータが与えられる。(ここでの設定は必要ない)
２１＋（KV*９）＋１０＋NX＋NY＋１行目
カラム

1～

変数名
FORMAT
内容

MBAR
I5
線境界の X 座標

２１＋（KV*９）＋１０＋NX＋NY＋２行目～
カラム

1～

6～

11～

16～

21～

変数名
FORMAT
内容

MBL
I5
海域段階

MBX

MBY

MBZ

MBT

Å

Å

Å

Å

X 座標

Y 座標

Z 座標

タイプ

118

・・MBAR 回分必要２１＋（KV*９）＋１１＋NX＋NY＋MBAR＋１行目
カラム

1～

変数名
FORMAT
内容

MARK
A3
データ区切り'END'

２１＋（KV*９）＋１１＋NX＋NY＋MBAR＋２行目
カラム

1～

6～

11～

16～

変数名
FORMAT
単位
内容

ＩＡＲＥＡ
I5

ＩＢＵＮ
I5

UNBR
F10.0

海域数

分割数

水深変換係数

AMSL
F10.0
cm
加算水深

*1 入力する水深値を cm 単位に変換するパラメータ，ここで入力した値を計算結果に乗じて出
力する。
例

入力値
UNBR
変換後(cm)
1.0
1
1
1.0
10
10
1.0
100
100
*2 一般的に水深値の読取りは海図を用いて行なうが、海図の基準面は基本水準面となっているため、平均水面に補正する値。
２１＋（KV*９）＋１１＋NX＋NY＋MBAR＋３行目
カラム

1～

6

11～

16～

変数名
FORMAT
内容

MMX
I5
Ｘ軸のメッシュ
数

MMY

IMOX

IMOY

Å

Å

Å

Ｙ軸 のメッシュ 大海域上のＸ原 大海域上のＹ原数
点
点

２１＋（KV*９）＋１１＋NX＋NY＋MBAR＋４行目
カラム

1～

変数名
FORMAT
内容

FMT
A50
データ読込みFORMAT119

２１＋（KV*９）＋１１＋NX＋NY＋MBAR＋５行目
カラム

1～

6～

61～

66～

変数名
FORMA
T
内容

H（１）
H(2)
前レコー
ドに従う
水深

H(13)

H(14)

カラム

1～

変数名
FORMAT
内容

MARK
A3
データ区切り
'END'

【1-5

水深

ml3eco.drw】

上３行はダミーである。（作図用のヘッダー）
行
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15

....+....1....+....2....+....3....+....4....+....5....+....6....+....7..
60
3.0
0.5
Å----1～３行目は未使用
11Å----モニター点の出力地点数
2871 東京都（st-5 隅田川河口)
Å----モニター点の格子番号（X,Y）、地点名---個数分
42
70 東京都（st-8 荒川河口）
68
78 千葉県（東京湾 3 三番瀬沖）
81
68 千葉県（東京湾 5 いなげの浜沖）
53
57 千葉県（東京湾 8 湾中央）73
54 千葉県（東京湾 9 五井沖）
46
23 千葉県（東京湾 15 木更津沿岸）
35
46 神奈川県（浮島沖）
24
36 神奈川県（扇島沖）
1527 神奈川県（本牧沖）26
2 神奈川県（第三海保東）

120

3.3.5システム利用時の流れ

改良システム利用時に必要なユーザーの操作は、入力条件変更(変更が必要なときのみ)、計算ケース（施策別）の指定、パラメータ設定、計算の実施、計算結果表示である。その他の手
続きは、システム側で実行されることになる。
ユーザー処理シ

システム処理設定項目選択

外部アプリケーションの起動

（200）

条件入力

ス

施策設定
（300）

テ

既往設定施策の選択

設定施策の表示

施策の指定

（301）

施策データのセット
（302）

ム

表示指定

フレーム表示*コ

計算パラメータ
設 定（400）

水循環モデル
パラメータ設定*

・準備データ(tran.dat)
・定数データ(tntank.dat)
・ｺﾝﾄﾛｰﾙデータ(cntrol.dat)
・雨量等水文データ

（410）

平面分
布

Import

外部作成
ファイル

ン

パラメータ修正

Import物質流動モデル
パラメータ設定

外部作成
ファイル

・定数データ(paramet.dat)

ト

（420）
パラメータ修正

東京湾モデルパ
ラメータ設定条件(ml3eco.dep)
係数(coeffcient.txt)

外部作成
ファイル

ロ
パラメータ修正計算のコントロール設定（計算
年度、実行モデル、対象河川）

モデル実行
ー

（500）

ル
画

水循環モデルの実行

Disk保存物質流動モデルの実行

Disk保存東京湾モデルの実行

Disk保存保存先
フォルダ指定計算結果保存面*

結果の可視化

*

（600）

表示指定
グラフ

（100）

終 了
平面分
布

凡例

定義済み処理通常処理

図- 3.3.5.1

外部
ファイル

手入力

システムフロー図

なお、操作マニュアルを巻末に参考資料として整理した。
121

表示（XXX）画面No