泥炭草地圃場における排水改良効果の定量評価
Quantitative Evaluation of Improvement of Drainage for Grass Field on peat
植屋 賢祐

松中 慧太

（UEYA Kensuke）（MATSUNAKA Keita）

25

I. はじめに

図-1は，暗渠排水と周辺排水路の整備が行われた
圃場における地下水位の観測事例である。図-1(a)は

泥炭は，未分解の植物遺体を多く含む特殊土であり，改良前，図-1(b)は改良後であり，改良前に比べて改わが国では，その大部分が北海道に分布する。

良後は雨量が多く，地下水位は全般に高く推移した。

泥炭農地では，泥炭地の形成過程とその理工学的このケースのように，時系列データを単純に比較し5 特性から，経年により不陸や地盤の低下を生じ，過湿30 ただけでは排水改良の効果が判然としない。事業期間障害など農用地としての機能低下が生じている事例が中に事業効果の発現を確認したい場合など，限られた見られる。北海道内各地の泥炭農地では，そのような期間の調査から改善効果を定量評価する必要がある

農用地の機能回復を目的として，国営総合農地防災

場合には，期間によってばらつきのある降雨の影響を事業が展開されている。

どう排除するかが課題となる。

10

本報告は，国営総合農地防災事業が実施された泥炭地の草地圃場を対象に，排水改良工事の効果を理

35

III. 排水改良効果の定量評価手法

工学的な側面から定量評価することを試みた。すなわ

1. 調査圃場の概要

ち，泥炭地の代表的な水文環境指標とされる地下水調査対象の圃場は，河川下流部に位置した標高3m，位に着目し，地下水位変動タンクモデルによる定量化

面積6ha程度の低平な泥炭地に位置する牧草地であり，15 手法を適用して，事業の有効性を考察した事例を紹介する。

暗渠排水の整備と付帯排水路の改修が行われた。
40

2. 排水改良の定量評価(1) 定量化の方法

II. 定量評価の必要性と課題

排水改良を定量評価する手法として，地下水位変対象地区では，泥炭に起因した過湿被害の防止な動を係数化して表現した既往の研究1)に倣い，地下水どを目的として，排水路整備を行う農地防災工，およ位変動タンクモデル（以下，タンクモデル）を用いた。

20 び暗渠排水・整地等の農地保全工による排水改良が45(2) タンクモデルの設定進められている。事業の実施によって，事業目的に見対象の牧草地は，前歴の国営農地開発事業により合った効果が発現されることが期待されるが，そのこと排水系が整備された。このため周辺流入はないと考え，を明確にするには，排水改良の効果を定量的に評価

地下水の供給源は降雨のみを想定した2)。流出孔は，することが必要となる。

比較的早い地下水位変化に係わる上段孔と，それ以50 外の比較的緩やかで定常的な逓減に係わる下段孔の2孔を想定し，図-2に示す1段タンクを設定した。
hn = h( n −1) − A1 × (h( n −1) − H A1 ) − A2 × (h( h −1) − H A 2 ) + C × R − E p
降雨 R55

蒸発散 E p

地下水位 hn
流出孔(上段)(a) 排水改良前

流出孔(下段)(b) 排水改良後

図-1 地下水位経時変化図

図-2
1※C：水位上昇係数
暗渠(施工後)HA2

A2 (流出係数)
早い逓減(施工後の暗渠など)

HA1

A1 (流出係数)
遅い逓減，地下浸透など

タンクモデルの概念図

(3) 排水改良の検証

排水改良前の地下水位観測値と，当該観測値より最適化したタンクモデルを用いて降雨をもとに計算した地下水位の計算値との比較を図-3(a)に示す。モデ
5 ル精度は相対誤差で表現すると10.0%と良好である。

図-3(b)は，図-3(a)のタンクモデルを用いて多雨年の降雨条件で計算した地下水位，すなわち改良前の地下水位をタンクモデルで推定した結果を，排水改良※GL-0.5m 以上が 3 日以上：一時過湿，7 日以上：常時過湿後の観測地下水位と比較した図である。この場合，相図-510 対誤差は33.3%となり，モデルの適合性は当然ながら低下した。しかし，同図に斜線で示した観測値と計算えた地下水位変動のシミュレーションである。地下値との乖離は，農作物の生産性に影響を及ぼす過湿水位の上昇後，GL-0.5m までの低下に要する時間は，状態の改善を定量的に示している。

改良前の 7 日に対して，改良後は 3 日に短縮した。
∑R=185.5mmm

3.639281【相対誤差393863.410.0%】
雨量(アメダス)

地表面 GL=3.26m

観測値

3.2

地下水位WL(m)

この結果は，地下水の早期低下による農作業効率の8070

計算値

40 回復の可能性を示唆している。

60

3.050
GL=-0.5m

2.8

402.6302.4202.2

雨量(mm/日)

15IV. おわりに

102.0

地下水位変動は降雨の多寡に影響を受けるため，06.16.116.21

7.1

7.11

7.217.31
日

8.108.20

8.30

9.99.199.29

排水改良前後の単純な比較による排水改良効果の定付(a) 排水改良前の最適モデルを排水改良前に適用∑R=471.0mmm

3.6【相対誤差 33.3%】

3.4

雨量(アメダス)

観測値

計算値

3.2
地下水位WL(m)

量把握は困難である。そこで，泥炭地における地下水

80

45 位変動を地下水位タンクモデルで表現することにより，70

地表面 GL=3.26m

60

3.050
GL=-0.5m

2.840

2.630

2.420

2.2

102.0

異なる係数設定や任意の降雨を入力とした排水改良

雨量(mm/日)2025

一雨を対象とした排水改良の定量評価

前後の最適モデルによる地下水位変動のシミュレーシ
ョンを試みた。泥炭農地において排水改良による圃場

06.1

6.11

6.21

7.1

7.117.217.31日

8.10

8.20

8.30

9.9

9.19付

(b) 排水改良前の最適モデルを排水改良後に適用

図-3の機能回復を定量的に評価する一手法として，有効な9.29

排水改良前の最適モデルの検証

50 手法であることが確認できた。

またこの手法により，継続的にモニタリングされた地次に，排水改良後の地下水位観測値と，当該観測下水位変動を係数化することで，経年による排水機能値より最適化したタンクモデルを用いて降雨をもとに算の低下を定量評価することも可能である。排水性能の30 出した地下水位計算値との比較を図-4に示す。相対照査および排水設備の更新適期の判断指標など，スト誤差は8.0%であり，図-3(a)と同等の精度が得られた。

55 ックマネジメントの視点からの適用も視野に入る。

以上から，排水改良前後の最適モデルはそれぞれ

参 考 文 献異なっている。各モデル係数を表-1に示す。

1) 梅田安治：泥炭地の水文特性に関する研究 昭和 58・59 年

図-5 は，排水改良前・後の両最適モデルに暗渠度科学研究費助成金(一般研究 B)研究成果報告書，pp.12～

35 の計画排水量に相当する 30mm/日の一雨降雨を与24，(1985)

60 2) 井上 京：泥炭地の地下水変動による水文環境評価に関する研究，pp.42～44，(1996)
植屋 賢祐
1969 年 北海道に生まれる

65

1992 年 弘前大学農学 部農業工学科 卒業
1992 年 (株)アルファ技研入社

図-4 排水改良後の最適モデルの検証

表-1係数

HA1 EL(m)HA2 EL(m)
A1
A2C

現在に至る

最適モデル係数の比較
排水改良前
0.00
2.50
0.00010.0037
14.5

排水改良後
0.00
2.32
0.00010.00437.070

松中 慧太
1984 年 広島県に生まれる
2007 年 帯広畜産大学 畜産学部畜産 科学 科卒業
2007 年 (株)アルファ技研入社
現在に至る

2