第 48 巻第 10 号 2006 17 丸亀平野における地下水流動の調査と数値解析 * 濱崎修 ** 河原能久 1. はじめに丸亀平野は, 瀬戸内気候区に位置し, 自己流域の降水量が少ないため, しばしば渇水被害を受けてきた 近年でも1994 年の渇水 以降平 6 渇水と略記 は厳しく, 甚大な被害を与えた このため, 丸亀平野では, 満濃池などの溜池の造成に代表されるように水の確保に苦心するとともに, 水利用について種々の工夫を凝らしてきた 1974 年に香川用水 吉野川流域からの導水 の運用が開始され, 農業用水及び都市用水の取水形態は大きく変貌した また, 流域の社会 経済情勢の変遷に伴う灌漑面積の減少や都市用水 排水量の増大等により, 流域の水循環系も大きく変化してきた 国土地理院の水準測量結果 1991 年 によれば, 近年, 丸亀平野では海岸部での地盤沈下 地下水塩水化が認められている また, 平野部においても, 水源の変更 不圧地下水から被圧地下水に変更 による被圧地下水位の低下がみられる このような状況に対応するため, 現在, 地下水利用者等が香川中央地域地下水利用対策協議会を組織し, 自主規制を行っている 1 また, 香川県が地下水利用を推進するための調査を行っているが 2, 地下水流動の実態は十分には明らかにされていない 丸亀平野においても, 限られた水資源を有効に利用するとともに, 流域の健全な水循環系を構築していくことが急務である 本研究では, 丸亀平野の適正な水管理システムの構築を目標として, 水循環において重要な役割を果たしている地下水や農業用水の実態を調査するとともにMIKE-SHEモデルによる水循環系の解析を行った ここでは, 地下水流動に関して得られた成果を中心に報告する 2. 丸亀平野の概要 2.1 対象領域本研究の対象流域である丸亀平野は, 瀬戸内海に面し, 讃岐平野の中央に位置する平野である 地形的には平野上流部の山地とその下方に発達する段丘 台地および中 下流域に分布する沖積扇状地から成る 丸亀平野流域 丸亀平野とその涵養域を包括する流域 の面積は約 320km 2 である 平野中央部には, 一級河川の土器川と二級河川の大束川, 金倉川, 弘田川の合計 4 河川が並行して流れている 図 -1 参照 大束川弘田川金倉川土器川 図 1 丸亀平野と観測地点丸亀平野の土地利用については,1976 年以降, 減反や農作物転換により, 水田面積が減少し, 宅地化が急速に進展している 現在, 流域内の土地利用は, 山林が約 43%, 水田が約 28%, 畑が約 8%, 宅地 市街地が約 10%, その他が約 11% となっている 本研究に関連する観測項目および観測地点を表 -1に示す * はまざきおさむ国土交通省香川河川国道事務所 ** かわはらよしひさ広島大学大学院工学研究科
18 地下水技術 観測項目河川流量 不圧地下水位 被圧地下水位 2.2 水理地質 1 帯水層 表 -1 観測所観測所丸亀橋, 祓川橋, 常包橋丸亀 No.2, 丸亀 No.4, 丸亀 No.5, 丸亀 浅井戸, 左岸 6.8k, 右岸 9.6k, 坂出 川津, 西坂元丸亀 深井戸, 豊原小学校, 高憧神社 既往の地質調査資料を総合的に分析して丸亀平野の水文地質構造図 図 -2 を作成した 丸亀平野の水文地質構造は4 層に大別できる 最上位の沖積層 A 層 は, 砂礫質主体で約 10~20m 層厚で不圧地下水が分布し, その下位には粘土 シルト主体の難透水層の段丘堆積層 B 層 が約 20~40m 分布し, 更にその下位には被圧地下水が分布する砂泥互層の三豊層群 C 層 が約 40 ~80m 現れ, 平野部の基盤は領家花崗岩類で構成される 2 透水係数及び比貯留係数既往の現場透水試験資料に基づいて各帯水層の透水係数を一次設定した A 層の透水係数は 1 10-4 ~10-5 m/s であり, 砂礫主体の層に帯水 すると考えられる B 層の透水係数は1 10-8 ~ 10-9 m/s で粘性土が主体の不透水層である C 層の透水係数は, 透水試験資料がないため, 深井戸の揚水実績資料から計算により求めた C 層の透水係数は5 10-5 m/s で砂泥互層中の砂層中に帯水すると考えられる 比貯留係数については, 既往の揚水試験結果 不圧帯水層 :2.5 10-4 に基づき設定した このような透水係数及び比貯留係数の値は, 後述するシミュレーション結果と観測結果との比較を通して, 調整された 最終的に設定した係数の値は後述の表 -3, 表 -4に示されている 各層の透水係数の分布を図 -3に示す A 層の透水係数は土器川下流の沿川付近で最も高い B 層の透水係数は各河川沿いで高い値を示すのに対して,C 層の透水係数は丸亀平野下流域一帯で大きい 3. 丸亀平野における水循環系の実態調査 3.1 地下水位の経年変化 1 不圧地下水位図 -4は, 図 -1に示す観測地点, 丸亀 No.4, 丸亀 No.5, 丸亀 浅井戸, 左岸 6.8k, 西坂元, 右岸 9.6k 地点における地下水位の経年変化を表している 図より, 各地点での地下水位の変動特性として, 以下のことが挙げられる 図 -2 丸亀平野の水文地質構造
第 48 巻第 10 号 2006 19 a A 層 a 丸亀 No.4, 丸亀 No.5 b B 層 b 丸亀 浅井戸, 左岸 6.8k C C 層 図 -3 各層の透水係数 c 右岸 6.8k, 右岸 9.6k 図 -4 不圧地下水位の経年変化
20 地下水技術 1 丸亀 No.4は土器川河口部の堤内側の砂州上にあり, 地下水位は夏期に約 0.5m~0.7m 上昇する 平 6 渇水の影響を受けて水位はやや低下したが, その後ゆるやかに回復した 2 丸亀 No.5は丸亀平野下流部の河川沿いの砂州上に位置する 地下水位はやや大きく, 夏期に約 1~2m 上昇する 平 6 渇水の影響は最低水位の顕著な低下という形で表れた 3 左岸 3.8k 地点に位置する丸亀 浅井戸 は丸亀平野下流部の砂州上に位置し, 周辺には旧河道の痕跡がみられる 地下水位は, 夏期に約 1.0m~1.5m 水位上昇する 経年的にみると, 観測開始時点より水位は低下傾向にある 最低水位は, 平 6 渇水時とその後に低下し, 回復が鈍い 4 左岸 6.8kは土器川河川沿いの氾濫源にあり, 地下水位は夏期に約 2m 水位上昇する 経年的な傾向は丸亀 浅井戸 と類似している 5 右岸 6.8k 地点に位置する西坂元は扇状地に位置し, 地下水位は夏期に1~2m 水位上昇する 平 6 渇水の影響は限定的であり, 長期的な変動傾向は認められない 6 右岸 9.6kは土器川河川沿いの氾濫源に位置し, 地下水位は夏期に1~2m 上昇する 平 6 渇水の影響は少なく, 長期的な変動はみられない 以上の観測結果と地形分類図から推察すると, 土器川の左岸側にある丸亀 浅井戸 及び左岸 6.8kは土器川沿いの丸亀平野の地下水盆のほぼ中央部に位置する 図 -5の不圧地下水面の深度や旧河道等の情報 図 -5の作成については後述 から知られるように, 当該地点は地下水脈 2の影響下にあり, 地下水脈 1 及び3に比べて, 地下水涵養量は少ないものと推察される 一方, 右岸側 6.8k にある西坂元は, 基盤岩が浅く地下水盆の縁に位置しており, 周辺には旧河道の痕跡がみられる 当該地点は土器川左岸側から右岸側を流れる地下水脈 1の影響を受けて, 地下水涵養がなされているものと推察される 右岸 9.6kの地点は, 土器川沿いの氾濫原内に位置するが, 沖積層が薄い岡田台地に隣接しているという, 他の観測地点とは異なる特徴を有している 当該地点は, 土器川の伏流水が左岸側か ら右岸側に流れる地域に位置し, 地下水涵養がなされているものと推察される 図 -5 不圧地下水面の地表面からの深度 2 被圧地下水位丸亀 深井戸, 高憧神社, 豊原小学校における被圧地下水の経年変化を図 -6に示す 図より次の特徴を指摘することができる 地下水位 T.P. m 地下水位 T.P.m 10.0 5.0 0.0-5.0-10.0-15.0-20.0 丸亀 深井戸 丸亀 浅井戸 -25.0 S S S S S S S S S S S S S H H H H H H H H H H H H H H H 52 53 54 55 55 56 57 58 59 60 61 62 63 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15-8.0-10.0-12.0-14.0-16.0-18.0-20.0 豊原小学校 50 年度 55 年度 60 年度 2 年度 6 年度 高憧神社 7 年度 8 年度 9 年度 10 年度 11 年度 12 年度 図 -6 被圧地下水の経年変化 13 年度 14 年度 15 年度 16 年度 17 年度 1 丸亀 深井戸 では大きな水位変化が記録さ れている 詳細に検討すると,1977 年 ~1979 年の水位低下は, 灌漑期に発生しており, 農業用水の揚水に起因すると推測される 一方, 1985 年 ~1987 年の水位低下は, 非灌漑期であ
第 48 巻第 10 号 2006 21 り, 上水道の揚水によるものと考えられる 1988 年 ~1993 年において, 水位がほぼ安定しており, この期間中の揚水量が, この地域における安定揚水量であると考えられる 1994 年の渇水 平 6 渇水 では, 深層地下水位は過剰な汲み上げにより著しく低下した 1996 年から2001 年にかけて3m 程度の回復が見られるが, これは揚水量の減少によると推測されている この地域の被圧地下水は, 周辺での揚水の影響を強く受けるが, 汲み上げが停止すれば水位は比較的早く回復する 2 高憧神社は丸亀平野下流のほぼ中央に位置し, 旧河道の痕跡に近接している 地下水位は 1980 年をピークに, その後低下し, 平 6 渇水の影響を受けて最低水位に達した 水位は2001 年 平成 13 年 まで徐々に回復したが,2002 年の冬期渇水以降, 低下し続けている 3 豊原小学校は丸亀平野下流の金倉川左岸に位置し, 旧河道の痕跡がみられるところにある 水位は,1980 年のピーク後に低下したが, 高憧神社に比べて低下量は小さい 平 6 渇水の影響も顕著ではなく,1999 年 平成 11 年 まで水位はほぼ横這いで推移した その後若干の上昇傾向を示していたが,2002 年の冬期渇水を経て, 水位は低下を続けている 水位の変化は, 平 6 渇水を除けば, 高憧神社のそれと類似している 3.2 地下水位の平面分布と地下水の流動丸亀平野における地下水位の平面分布と地下水流動状況を把握するため, 平成 15 年から17 年の灌漑期と非灌漑期に地下水調査を実施した 浅井戸 105 箇所, 深井戸 22 箇所 浅層地下水面の地表面からの深度を前述の図 -5に, 浅層および深層地下水の水位の等高線を図 -7に示す 丸亀平野では, これまでに旧河道とその周囲の透水性の高い氾濫原が確認されている 今回, 丸亀平野で実施されている地質調査ボーリング資料 1,160 本 を分析することにより旧河道と氾濫源の痕跡を整理した 図 -5にはその結果を不圧地下水面の地表面からの深さとともに示している その結果,3 系統の水脈, すなわち, 土 器川の伏流水が左岸側から右岸側へ流動する水脈 図中の1, 平野中央部を流れ再び土器川下流に環流する水脈 図中の2, 金倉川方向へ多く流れ旧河道沿いを流下する水脈 図中の 3 の存在が推測された 図 -7 地下水位の等高線 標高表記 浅層地下水位は, 概ね地表の地形に沿った分布であるが, 平野の東西に位置する氾濫原では, 地表面から地下水位までの深度が平野中央部に比べ2~3m 深い また, 土器川の瀬切れ発生区間と氾濫原や旧河道の位置が一致していることから, 河川水と浅層地下水は, 密接に関連していると考えられる 図 -7より, 灌漑期と非灌漑期の浅層地下水を比較すると, 地下水位の標高 15~20m,20~25m の範囲において約 1.0~1.5mの変動幅がみられる 一方, 深層地下水は, 浅層地下水に比較して地下水位の勾配が緩く, 地表の地形に沿った分布となっている 3.3 出水 の流出量調査丸亀平野の農業用水の特徴的な取水形態として, 出水 ですい と呼ばれる地下水取水がある 土器川下流の沿川で約 60 箇所 ポンプまたは自然流下 の出水があり, その中から自然流下で送水している出水 28 箇所 について, 流出量の観測 2003 年 9 月 を行った そのうち, 土器川左岸 6k/380に位置する竜王出水について, 流出量と前 10 日間降雨量の間の関係を図 -8に示す バラツキはあるものの, 両者に相関が認められる
22 地下水技術 前 10 日間降雨量 mm 100 90 80 70 60 50 40 30 20 10 0 0.00 0.02 0.04 0.06 0.08 0.10 出水流出量 m3/s 図 -8 出水流出量と前 10 日間降雨量の関係 竜王出水 3.4 農業用排水系統丸亀平野においては, 灌漑期における水田からの地下水涵養が重要となることから, 本研究では, 対象流域における灌漑実態を明確にするため, 現地において農業用排水系統を調査するとともに, 関係市町並びに土地改良区にヒアリング アンケート調査を実施した 丸亀平野における用排水系統を図 -9に示す 対象流域では, 香川用水, 河川水, 溜池 満濃池含む, 出水 湧水, 井戸 浅井戸 深井戸 が 農業用水の水源として利用されている 土器川等の4 河川に対し, 香川用水が横断しており, 香川用水の導水路と約 80 箇所の溜池がネットワーク化された複雑な水利用形態となっている 4. 水循環解析モデルによる解析 4.1 水循環解析モデルの概要本研究では, 対象流域における水循環系を定量的に把握するため, 水循環素過程を物理式で表現し, 表流水と地下水を連動して解析できる MIKE-SHEモデルを適用した MIKE-SHEモデルは, デンマーク水理研究所 DHI 等で開発されたグリッド型の水循環解析モデルであり,1 遮断 蒸発散,2 表面流 河道流,3 不飽和流,4 飽和流,5 帯水層と河川との水交換,6 融雪の6つの部分から構成されている 図 -10にMIKE-SHEモデルの基礎式を示す 表面流と河道流は,St. Venant 方程式を簡略化したdiffusive waveモデル 河道流には dynamic waveモデルも適用可 と連続の式により解析され, 表面流は2 次元モデル, 河道流は1 次元モデルで表現されている 土壌層における不飽和流には1 次元 Richards 方程式を用い, 飽和 海 海 海 多度津 幸町 多度津 桜川 多度津 青木 多度津 青木 海 聖池 天満池 柳池 西汐入川 長太夫池 北鴨幹線 道福寺幹線 庄幹線 青木幹線 弘田川 金倉川 大束川 飯野幹線水路 源田池 土器川 道池金丸池上池八丈池 中原池宮池菱池 丸亀市中府町 新池 前池 藤高池 蓮池 山崎池 下池 大林池 馬池庄之池宮池大池籠池 北原池 先代池田村池大井池新池中池伊予勢池 水附池 蓮池 瓢池平池道池前池 白方池 海 古子川 宮池 滝川幹線 村上池 金倉寺幹線 蓮池幹線 丸亀幹線 金倉幹線 宝憧寺幹線 上池 図 -9 丸亀平野における農業用排水系統 皿池 上池 宮池 菰池新池上池千代池新池 要池 買地池 清水川 丸亀市 丸亀市配水池 善通寺市 多度津町 綾歌町 西山受水池 丸亀市 宇多津町 坂出市 新池 買田池 牛池 綾川 宮後池新池藤波池常磐池 水橋池 仁池 岡池導水路 為久導水路 岡 1~4 号為久池 上水 工水 大窪池 宝憧寺幹線川西支線 矢野池 新池 宝憧寺池 桝池 買田池 土路池 椎尾池 仁池導水路 天神池 大窪導水路 熊ヶ池 瓢箪池 竜社池 中部浄水場 豊原幹線 満賀池 西池 金倉川左岸幹線 水橋池 仁池 天神池 大窪池 西池 飯野幹線 宝憧寺幹線 丸亀幹線 中部上工水 金倉川畑地かんがい施設 善通寺 善通寺 弘田町 金倉支線 真野上井水路 買田幹線 管池 香川用水 小津森 打越池 左岸地区 小津森池 今滝池 小津森導水路 吉野方面 羽間導水路 羽間池 吉野幹線 打越導水路 大川頭首工 満濃池 亀越池 備中地池 天川頭首工 満濃池導水路 財田川流域 野口ダム 凡例出水香川用水分水地点
第 48 巻第 10 号 2006 23 降雨 降雪 蒸発散 樹冠遮断モデル 遮断降雨から 土壌 水面から 根系層から 有効降雨 融雪モデル 人工取水と涵養 表面流出 浸透 2 次元表面流 河道流モデル 水面昇降 湖沼 貯水池 根系層 1 次元不飽和流モデル浸潤面を通しての水移動 3 次元飽和流 地下水 モデル 境界条件を通しての水移動 地下水位 i 層 i+1 層 河川潤辺 w 河川水位 地下水位 図 -10 MIKE-SHE モデル中の基礎式 流はダルシー式と連続式に基づく3 次元モデルにより解析される また, 不飽和層底部の土壌水分と地下水位の変動が連成解析される 帯水層と河川の水交換は, ダルシー則を用いて算定され, 河川水位に対して地下水位が高い場合には河川へ流出し, 地下水位が低い場合には河川水は地下へ涵養される 解析手法の詳細については文献 3-5 を参照されたい 4.2 水循環解析モデルの構築地形条件の再現性及び計算時間を考慮して, 計算格子の大きさを 250 mの正方形とし, 対象流域を5,518グリッドに分割した 表層モデルの作成にあたっては, 数値地図 50mメッシュ標高を抽出し, グリッド毎に平均化してグリッド標高を設定した 地下水層のモデル化にあたっては, 図 -2に示すA 層を不圧帯水層,B 層を難透水層,C 層を被圧帯水層とし, 既存のボーリング柱状図 1,160 箇所 をもとに内挿法によりグリッド毎の各層の基底標高を設定した また, 河道縦横断図から横断形状を設定し, 河道モデルを作成した 測量成果が存在しない区間については国土数値情報の河道位置, 河道長, 河床高を取得し, 矩形断面 幅 5m 深さ2 m を設定した 河川流量の再現性向上を図るために, 再現期間中, 降雨によらず地下水が常時飽和状態にな ることを避けるため, 河川基底流量の計算にドレーン モジュールを適用した 4.3 水田のモデル化の方法水田は, 人為的な灌漑により水位 水量がコントロールされ, また, その構造上, 大きな貯留効果を持つなど, 自然地とは異なった特性を有する さらに一般的に, 農業用排水の系統が複雑であり, 灌漑の実態が不明であるため, 農業用排水の流動に対する物理的なモデルの構築は容易ではない 本研究では, 灌漑による水循環系を可能な限り定量化するため, 下記に示すように, 降雨, 蒸発散, 灌漑, 浸透等を表現する水田モデルの開発を試みた 水田の水管理状況をモデル化するため, 水田の減水深 蒸発散 浸透 に一致するように, 不飽和浸透の特性を設定した 灌漑区域毎の水源位置 香川用水, 溜池, 出水, 井戸, 用排水系統の調査結果をもとに取水地点, 供給区域を設定した 対象流域では灌漑区域毎に複数の水源を有しているため, 各灌漑区域の取水地点を設定するにあたって, 使用水源タイプ別に区分し, 水源の優先順位を設定した 減水深法により減水深, 田面降水量等から単位面積あたりの必要水量を算定し, 各グリッドの
24 地下水技術 灌漑面積を乗じて, 農業用水使用量を推定した グリッド毎の必要水量に対して, 実績取水量が既知の水源からの供給量を差し引き, 残分については地下水取水量をパラメータとして地下水位を再現するように, 表流水と地下水の供給量を配分した 4.4 境界条件 初期条件の設定境界条件として, 表 -2に示す条件を各モデルに設定した 各グリッドの生活用水, 工業用水の人工系の水量を定量化するにあたっては, 各種統計書および関連市町へのヒアリング等から水源別の取水量を整理し, 境界条件として設定した 6 また, 本研究ではPenman 式を用いて可能蒸発散量を算定し, 古藤田による補正係数 7 を乗じて実蒸発散量を推定した 初期条件として, 非灌漑期 1 月 の地下水一斉観測より得られた浅層 深層地下水位の面的な分布を各グリッドに与えた 表 -2 境界条件の設定 を実施し, 河川流量及び地下水位の実測値と計算値を比較して, 各種パラメータの同定を実施した 表層モデルでは, 土地利用分類 国土数値情報 1/10 細分区画土地利用分類データから抽出 等により各グリッドを分類し, 等価粗度係数および窪地貯留量を設定した 不飽和層モデルでは, 土壌分類 国土数値情報 土壌分類データから抽出 や土地利用分類等により各グリッドを浸透度別 8 に分類し, 不飽和浸透特性を設定した 不飽和浸透特性としては, 不飽和透水係数と水分特性曲線が必要であり, これらについては, 西垣ら 9 により整理されたvan Genuchten 式への適用結果を参考にパラメータを設定した 地下水層モデルでは, 前述のように, 各層の透水係数については, 現場で行われた透水試験結果を基に一次設定し, 再現計算により最適値を求めた 表 -3, 表 -4に最終結果を示す 比貯留係数に表 -3 A 層のパラメータ 表 -4 B 層,C 層のパラメータ 5. 数値解析結果 5.1 モデルパラメ-タの設定と同定近年の渇水年 平水年 豊水年を含むよう, シミュレーション期間を1994 年 1 月 1 日 ~1998 年 12 月 31 日の5 年間とした 日単位のシミュレーション 表 -5 河道モデルのパラメータ
第 48 巻第 10 号 2006 25 ついては, 文献 10 による一般的な値を適用した また, 地下水取水量 農業用水 をパラメータとして, 地下水位を再現するように設定した 河道モデルに必要な粗度係数及び漏水係数の値は, 河川流量の再現計算により同定した 最終的に設定した値を表 -5に示す 5.2 水循環解析の検証設定したモデルパラメータを用いで, 河川流量, 不圧地下水, 被圧地下水の水位を計算し, 解析モデルやパラメータの妥当性を検討した 河川流量 日流量 の計算結果を図 -11に示す 丸亀橋と祓川橋の両地点において, 平常時の流量については概ね実測値を再現できているが, やや規模の大きな出水に対しては改善の余地がある 図 -12に不圧地下水位の計算結果を示す 浅層地下水位については, 長期的な再現性は概ね良好であるが, 日変動については実績水位との間に若干の乖離がみられる 図 -13では被圧地下水位について, 計算結果と観測結果の比較を行っている 丸亀深井戸及び高憧神社では, 取水の影響による水位の低下現象がみられるが, 低下の度合いも精度良く再現できている 豊原小学校については, 他観測所でみられるような取水の影響を受けた水位低下部の実績値がなく, 水位低下の評価ができないが, 欠測を除いた実績値に対する再現性は概ね良好である いずれの地下水も長期的な傾向を概ね再現できており, 設定した透水係数及び比貯留係数の値が妥当であることを示している ハイエトグラフ 丸亀 降雨量 mm/ 日 0.0 50.0 100.0 150.0 200.0 降雨量 ハイドログラフ 丸亀橋 河川流量 m3/s 100.0 90.0 80.0 70.0 60.0 50.0 40.0 30.0 20.0 10.0 0.0 1994 1995 1996 1997 1998 実測値計算値 a 丸亀橋地点 ハイエトグラフ 祓川橋 mm/ 日 降雨量 0.0 50.0 100.0 150.0 200.0 降雨量 ハイドログラフ 祓川橋 河川流量 m3/s 100.0 90.0 80.0 70.0 60.0 50.0 40.0 30.0 20.0 10.0 0.0 1994 1995 1996 1997 1998 実測値計算値 b 祓川橋地点図 -11 日流量に関する計算結果と観測結果の比較
26 地下水技術 a 丸亀 No.4 b 丸亀 No.5 c 左岸 6.8k d 坂出 川津 e 右岸 9.6k f 西坂元 図 -12 不圧地下水位の検証
第 48 巻第 10 号 2006 27 a 丸亀 深井戸 b 豊原小学校 c 高憧神社 図 -13 被圧地下水位の検証 図 -14に非灌漑期における不圧地下水位の等高線を示す 1994 年 12 月と1998 年 12 月における計算水位 実線 と一斉観測結果 破線,2004 年 1 月 を併記している 詳細に見ると,1994 年 渇水年 では, 実績水位に比べ, 平野部全体の地下水が低下している その後, 地下水位は回復傾向にあり,1998 年 豊水年 においてようやく2004 年 平水年 の実績値と同程度まで地下水位が回復している なお, 計算水位は沿岸部の旧河道に沿って実測水位とやや異なっているが, 面的な再現性は比較的良好である その旧河道沿いの地下水位の差違は, その地域の透水係数を過大評価しているためかと推測している a 1994 年 12 月 渇水年 b 1998 年 12 月 豊水年 図 -14 不圧地下水の水位分布
28 地下水技術 以上のように,MIKE-SHEモデルを用いて丸亀平野の水循環系の大略は把握できるようになった 結果を図示しないが, このモデルを用いて, 現在, 丸亀平野における地下水涵養域の特定や水収支等を検討している 例えば,1996 年において, 年間降水量 969mmのうち, 約 56% が蒸発 蒸発散で大気に戻り, 残りの地表面に到達した水量のうち約 62% が流出し, 約 38% が地下に降下浸透したことが明らかとなった また, 人工系の使用水量の内訳は, 河川水 : 香川用水 : 地下水 = 約 1.0:1.2:1.5であり, 丸亀平野では香川用水と地下水に大きく依存していることが確認されている 6. おわりに本研究では, 丸亀平野における水循環系の実態を明らかにするために, 旧河道の特定, 地下水や農業用水の現地調査を行うとともに, 水循環解析モデル MIKE-SHE を用いた数値解析を行った その結果, 透水係数等のモデル定数値を適切に評価することによって, モデルが不圧 被圧地下水や河川の日流量を良好に算出することを確認した 今後, モデルの精度の向上をはかるとともに, 流域の水循環に関する課題を取り上げ, そのメカニズムの解明と対応策の検討を進める予定である 謝辞 : 用排水系統調査に御協力下さり, 貴重な資料を提供して頂いた関係県 市 町並びに満濃池土地改良区, 土器川右岸土地改良区連合に感謝の意を表します 参考文献 1 吉野文雄 : 香川中央地域の地下水, 香川中央地域地下水利用対策協議会, 2003. 2 香川県企画部水資源対策課 : 平成 10 年度香川県地下水利用推進調査総括報告書, 1999. 3 Singh, V. P.: Computer Models of Watershed Hydrology, Water Resources Publications, 1995. 4 Abbott, M. B. et al: An Introduction to the European Hydrological System - Systeme Hydrologique Europeen SHE 1: History and philosophy of a physically based distributed modelling system, Journal of Hydrology, Vol.87, pp.45-59, 1986. 5 Abbott, M. B. et al: An Introduction to the European Hydrological System - Systeme Hydrologique Europeen SHE 2: Structure of a physically based distributed modelling system, Journal of Hydrology, Vol.87, pp.61-77, 1986. 6 安陪和雄, 大八木豊 : 分布型流出モデルの広域的適用, 水工学論文集,Vol.46,pp.247-252, 2002. 7 古藤田一雄 :Estimation of River Basin Evapotranspiration, Environmental Research Center Papers, No.8, pp.1-66, 1986. 8 中野秀章 : 森林水文学, 共立出版,1980. 9 西垣誠, 竹下祐二 : 室内及び原位置における不飽和透水特性の試験及び調査法に関する研究, 岡山大学土木工学科,pp.30-95, 1993. 10 Anderson, M.P. and Woessner, W.W.: Applied Groundwater Model Simulation of Flow and Advective Transport, Academic Press, p.41, 1992.