Microsoft Word - 水資源工学_講義資料1,2,4_.docx

Size: px

Start display at page:

Download "Microsoft Word - 水資源工学_講義資料1,2,4_.docx"

えいしろういまいだ
5 years ago
Views:

水資源工学 ( 立川担当分 1 回目 :12 月 21 日資料 ) 流況評価の方法 (1) 講義の目標 :1 年を通じた河川流量の特徴を河川流況という流況を評価することは河川による水資源を把握することである今日の講義では流況を評価する基本的な手法と流況の量的な特性を表す物理量 ( 水文量 ) を理解する次に水資源確保の基本量である平均して 10 年に 1

1 水資源工学 ( 立川担当分 1 回目 :12 月 21 日資料 ) 流況評価の方法 (1) 講義の目標 :1 年を通じた河川流量の特徴を河川流況という流況を評価することは河川による水資源を把握することである今日の講義では流況を評価する基本的な手法と流況の量的な特性を表す物理量 ( 水文量 ) を理解する次に水資源確保の基本量である平均して 10 年に 1 回程度の割合で発生する渇水流量 ( 基準渇水流量 ) を推定する手法を理解する 1) 河川流量の時間変化 ( 流量ハイドログラフ ) 河川流量の時間変化を図示したものを流量ハイドログラフという流量データは水位流量曲線を用いて水位データを変換して得られることが大半である連続的な流量データから時間平均したデータを時間流量日平均したデータを日流量月平均した流量を月流量という水位流量曲線の例 2) 様々な河川のハイドログラフ最上川の年間の日流量変化吉野川年間の日流量変化 3) ハイドログラフから得られる河川流況の情報 4) 流況曲線一年間の日流量ハイドログラフを流量の多い順に並べ替えて大きい順に左から図示したものを流況 1

曲線という流況曲線から得られる流況の特徴量として豊水流量 ( 一年間を通じて 95 日はこれを下らない流量 ) 平水流量 ( 一年を通じて 185 日はこれを下らない流量 ) 低水流量 ( 一年通じて 275 日はこれを下らない流量 ) 渇水流量 ( 一年を通じて 355 日はこれを下らない流量 ) が得られる流況曲線の例 5) 基準渇水流量平均的に 10 年に 1

2 曲線という流況曲線から得られる流況の特徴量として豊水流量 ( 一年間を通じて 95 日はこれを下らない流量 ) 平水流量 ( 一年を通じて 185 日はこれを下らない流量 ) 低水流量 ( 一年通じて 275 日はこれを下らない流量 ) 渇水流量 ( 一年を通じて 355 日はこれを下らない流量 ) が得られる流況曲線の例 5) 基準渇水流量平均的に 10 年に 1 回の頻度で発生する渇水流量を基準渇水流量という新たに取水を予定する場合は基準渇水流量から河川維持流量と他の水利使用者の双方を満足する水量 ( 正常流量 ) を差し引いた水量の範囲内として決定される河川維持流量とは塩害防止各種排水の希釈浄化河道の維持河口埋塞防止水生動植物の生存繁殖等河川に関する公利の確保公害の除去若しくは軽減のため流水の果す機能を確保するための流量をいう 6) 水文頻度解析手法を用いた基準渇水流量の設定 10 年確率渇水流量 ( 基準渇水流量 ) の求め方は以下の通りである 1 数十年以上の日流量データを収集する 2 毎年の渇水流量を抽出する 3 水文頻度解析手法 * により渇水流量が適合する確率分布モデルとその母数を推定する確率分布モデルは渇水流量についてはワイブル分布を用いることが多い 4 非超過確率が 1/10 となる渇水流量を基準渇水流量とする * 参考図書 : 工学のための確率統計 : 北村堀 ( 編著 ) 朝倉書店例題で学ぶ水文学 : 椎葉立川市川森北出版

3 水資源工学 ( 立川担当分 2 回目 :1 月 4 日資料 ) 流況評価の方法 (2) 講義の目標 : 流況の特性として量的な特性と時間的な特性がある量的な特性として流況曲線や基準渇水流量がある時間的な特性としては降水量や水位流量の時系列的な特性がある今日の講義では流量データの時系列特性と時系列モデルの基礎を理解することを目的とする 1) 水文時系列時間の経過とともに変動する現象の記録を時系列 (time series) といい時間変化する降雨強度や河川の水位流量などの水文量を総称して水文時系列という 2) 水文時系列解析の目的と手順 1 水文時系列データの確率統計的な特性を分析する 2 時系列データの確率統計的な特性を再現する時系列モデルを構成する 3 構成した時系列モデルを用いて観測時系列と同じ確率統計的な特性を有する長期間の時系列を模擬発生させ水工施設の設計に利用するあるいは時系列モデルを用い現在までに得られた観測時系列から将来の水文量の変動を予測する 3) 定常時系列と非定常時系列時系列はその時間的は変動によって定常時系列と非定常時系列に分類される時間とともに時系列の確率的特性が変化しない時系列を定常時系列という一方ある期間で見た時系列の平均値が時間とともに変動したり平均値の周りのばらつきが時間とともに変化したりする場合など時間とともに時系列の確率統計的特性が変化する時系列を非定常時系列という 4) 定常時系列とその特性値離散的な時系列を,,, としこれらの変量の同時確率分布関数を考える時系列を表す確率分布関数が時間をずらしても変化しない場合すなわち,,, を時系列の同時確率分布関数としを任意の時間間隔を表すインデックスとして,,,,,, が成り立つ場合その時系列は強定常であるという特に二次までのモーメントについて μe E Var Var Cov, E Cov, E が成り立つ場合弱定常あるいは二次定常という E は期待値 Var は分散 Cov は共分散を求める演算記号であるこれらをもとに自己共分散関数自己相関関数が定義される Cov, E Cov, Var Var 3

4 5) 時系列モデル (AR モデル ) 定常時系列を過去の時系列の値と白色ノイズの線形和で表現したモデルを次数の自己回帰モデルあるいは AR モデル (AutoRegressive model) といい AR() と表すまた AR モデルに従う時系列を AR 過程という 1,2,, は自己回帰係数は平均値 0 で分散 σ の正規分布に従う白色雑音でありはそれより前の時系列とは無相関とするもっとも簡単な 1 次の自己回帰モデル AR(1) はである AR(1) のモデルパラメータは,, となるの観測時系列が得られておりそれが AR(1) 過程で表現できるとしてモデルパラメータを求めるまず μe であるまたの分散をとすれば 2 となるので 11 という条件のもとに 1 が得られる次に時間差が 1 の自己共分散関数は時間差が 2 の自己共分散関数はとなる一般に時間差 k > 0 についての関係が得られるので時間差 k > 0 の自己相関関数は / となる観測時系列データから標本自己相関関数を求めればの推定値を定めることができるが定まれば観測時系列データから得られる標本分散用いて 1 となり 1 次の自己回帰モデル AR(1) のパラメータが定まる 4

5 3 : 1 11 AR III AR(1) AR(1) AR(1) y n 3 { 1 f Y (y n ) = (y n c)ζ 2π exp 1 2 cλζ 3 y n ( ) } 2 ln(yn c) λ ζ (1) x n = ln(y n c) λ ζ (2) x n N(0, 1) y n 1 ρ y x n 1 ρ x ρ y = exp(ζ2 ρ x ) 1 exp(ζ 2 ) 1 (3) 1) 1 ρ y (3) ρ x 2) (2) y n x n σ 2 v = (1 a 2 )σ 2 x v n a = ρ x 3) AR(1) x n = ax n 1 + v n x n 4) (2) x n y n (3) X 1 N(m X1, σ 2 X 1 ), X 2 N(m X2, σ 2 X 2 ) Y 1 = e X1 Y 2 = e X2 C Y1,Y 2 Y 1,2 m Y1,2 σ 2 Y 1,2 m Yi ( ) = exp m Xi + σ2 X i, i = 1, 2 (4) 2 σ 2 Y i = m 2 Y i {exp(σ 2 X i ) 1}, i = 1, 2 (5) 5

6 X 1 X 2 C X1,X 2 X 1 + X 2 N(m X1 + m X2, σ 2 X 1 + σ 2 X 2 + 2C X1,X 2 ) (6) C Y1,Y 2 = E[Y 1 Y 2 ] E[Y 1 ]E[Y 2 ] = exp{m X1 + m X (σ2 X1 + σx C X1,X 2 )} m Y1 m Y2 = m Y1 m Y2 {exp(c X1,X 2 ) 1} (7) ρ(x 1, X 2 ) = C X1,X 2 /(σ X1 σ X2 )ρ(y 1, Y 2 ) = C Y1,Y 2 /(σ Y1 σ Y2 ) ρ(y 1, Y 2 ) = C Y 1,Y 2 = m Y 1 m Y2 (exp(c X1,X 2 ) 1) σ Y1 σ Y2 σ Y1 σ Y2 exp(c X1,X = 2 ) 1 = exp(σ X 1 σ X2 ρ(x 1, X 2 )) 1 exp(σx 2 1 ) 1 exp(σx 2 2 ) 1 exp(σx 2 1 ) 1 exp(σx 2 2 ) 1 (8) (3) AR(1) n y n x n = (y n µ n )/σ n (9) x n µ n σ n n x n N(0,1) AR(1) x n = a n,n 1 x n 1 + v n (10) a n,n 1 n n 1 x n σ 2 x 1 σ 2 v = (1 a n,n 1 2 )σ 2 x v n N(0, 1 a 2 n,n 1) (10) x n (9) y n a n,n 1 v n Thomas Fiering Thomas-Fiering n y n 3 { 1 f Y (y n ) = exp 1 ( ) } 2 ln(yn c n ) λ n (y n c n )ζ n 2π 2 ζ n y n (11) x n = ln(y n c n ) λ n ζ n (12) (11) (8) y n ρ y,n,n 1 x n ρ x,n,n 1 exp[ζ n 1 ζ n ρ x,n,n 1 ] 1 ρ y,n,n 1 = exp(ζn 1 2 ) 1 exp(ζn) 2 1 (13) 6

水資源工学 ( 立川担当分 4 回目 :1 月 18 日資料 ) 流況予測のためのモデル (2) 講義の目標 : 河川流量をシミュレーション発生させる手法は前回の講義で示した確率統計的な時系列モデルによる方法だけでなく物理的なモデルを用いて降雨を発生させそれを流出モデルを介して河川流量に変換することが考えられる地球温暖化によって将来の降水量が予測され

結果としての観測データを分析して現象の確率的な変動を再現しようとする一方で自然現象を物理的な方程式で表現するモデルも時系列データを発生させるために用いられるたとえば気候変動にともなう将来の気候推計データは物理法則にしたがう支配方程式を解いて将来の気温や降水量などの時系列データを生成し

その出力データはわが国の河川流域を対象とする水工シミュレーションモデルの入力データとしてそのまま利用できる解像度を有するようになってきている気象庁気象研究所の全球大気モデルの空間分解能は約 20km でありアメダス観測網に匹敵する空間分解能となっている全球を対象とする気候モデルとしては現時点では世界最高の空間分解能を有する気候モデルである

7 水資源工学 ( 立川担当分 4 回目 :1 月 18 日資料 ) 流況予測のためのモデル (2) 講義の目標 : 河川流量をシミュレーション発生させる手法は前回の講義で示した確率統計的な時系列モデルによる方法だけでなく物理的なモデルを用いて降雨を発生させそれを流出モデルを介して河川流量に変換することが考えられる地球温暖化によって将来の降水量が予測され流出モデルや河川流モデルを介して河川流量の変化が分析されている今日の講義では温暖化時の流量変化予測を例として河川流況のシミュレーション手法を理解することを目標とする時系列モデルは観測された時系列データに内在する確率統計的な特性を分析してその特性を再現する時系列を模擬発生させることに主眼が置かれる時系列データが生み出される物理的背景よりも結果としての観測データを分析して現象の確率的な変動を再現しようとする一方で自然現象を物理的な方程式で表現するモデルも時系列データを発生させるために用いられるたとえば気候変動にともなう将来の気候推計データは物理法則にしたがう支配方程式を解いて将来の気温や降水量などの時系列データを生成しそれを流出モデルを介して河川流量に変換して将来の流況が評価されている 1) 全球大気大循環モデル将来の気候を推計する大気大循環モデル (GCMs, General Circulation Models) の空間分解能が著しく高解像度化しその出力データはわが国の河川流域を対象とする水工シミュレーションモデルの入力データとしてそのまま利用できる解像度を有するようになってきている気象庁気象研究所の全球大気モデルの空間分解能は約 20km でありアメダス観測網に匹敵する空間分解能となっている全球を対象とする気候モデルとしては現時点では世界最高の空間分解能を有する気候モデルである将来気候の推計計算は温室効果ガス排出量の変化シナリオに従って温室効果ガスの濃度の時間変化が設定されその条件のもとで将来の気候が予測される気象庁気象研究所の全球 20km 格子大気モデルは A1B シナリオに従って温室効果気体の濃度変化が設定されており以下の期間のデータが提供されている 1979 年 1 月 ~2003 年 12 月 : 現在気候実験 2015 年 1 月 ~2039 年 12 月 : 近未来気候実験 2075 年 1 月 ~2099 年 12 月 : 21 世紀末気候実験図 1: 全球大気大循環気候モデルの模式図 ( 左 ) とそれが出力するデータ ( 右 ) 7

2) 流域モデルと流出モデル地球上の任意全域を対象とする 1km 空間分解能のグリッド型の分布型流出モデルを構成し,

方向のうちの最急勾配方向を流水方向と定めそれに従って一次元的に流れを追跡する図 2( 左 ) は流水方向データを標高データから定める模式図であり

右 ) 3) 大気大循環モデルの出力を流出モデルを介して河川流量に変換した結果最上川 ( 砂越地点 6500km2)

年間の河川流量を重ねて示している 4) 分析結果の例図 3:100 年確率年最大流量の変化 ( 左 ) と基準渇水流量 (10 年確率渇水流量 )

8 2) 流域モデルと流出モデル地球上の任意全域を対象とする 1km 空間分解能のグリッド型の分布型流出モデルを構成し, キネマティックウェーブモデルを用いて大気大循環モデルによる出力データを河川流量に変換する標高データを用いてグリッドごとに周り 8 方向のうちの最急勾配方向を流水方向と定めそれに従って一次元的に流れを追跡する図 2( 左 ) は流水方向データを標高データから定める模式図であり右図はこれに従って関東地方の流域モデルを示したものである図 2: 標高データを用いた流下方向の決定と流域モデル ( 左 ) と関東地方の流域モデル ( 右 ) 3) 大気大循環モデルの出力を流出モデルを介して河川流量に変換した結果最上川 ( 砂越地点 6500km2) での計算流量の時系列データを以下に示す左から順に現在気候実験近未来気候実験将来気候実験の結果であり一つの図に 25 年間の河川流量を重ねて示している 4) 分析結果の例図 3:100 年確率年最大流量の変化 ( 左 ) と基準渇水流量 (10 年確率渇水流量 ) の変化 ( 右 ) : 現在気候に対する 21 世紀末気候の比率を示している参考文献 : 立川康人, 滝野晶平, 藤岡優子, 萬和明, キムスンミン, 椎葉充晴 : 気候変化が日本の河川流量に及ぼす影響の予測, 土木学会論文集, 67(1), pp. 1-15,

森林水文水資源学 2 2. 水文統計豪雨があった時, 新聞やテレビのニュースで 50 年に一度の大雨だったなどと報告されることがある. 今争点となっている川辺川ダムは,80 年に 1 回の洪水を想定して治水計画が立てられている. 畑地かんがいでは,10 年に 1 回の渇水を対象として計画が立て

森林水文水資源学 2 2. 水文統計豪雨があった時, 新聞やテレビのニュースで 50 年に一度の大雨だったなどと報告されることがある. 今争点となっている川辺川ダムは,80 年に 1 回の洪水を想定して治水計画が立てられている. 畑地かんがいでは,10 年に 1 回の渇水を対象として計画が立て . 水文統計豪雨があった時, 新聞やテレビのニュースで 50 年に一度の大雨だったなどと報告されることがある. 今争点となっている川辺川ダムは,80 年に回の洪水を想定して治水計画が立てられている. 畑地かんがいでは,0 年に回の渇水を対象として計画が立てられる. このように, 水利構造物の設計や, 治水や利水の計画などでは, 年に回起こるような降雨事象 ( 最大降雨強度, 最大連続干天日数など