JP4006098B2 - 水田群の水管理システム，水田群の水管理方法，および，水田群の水管理を行うためのプログラムを記録したコンピュータ読み取り可能な媒体 - Google Patents
水田群の水管理システム，水田群の水管理方法，および，水田群の水管理を行うためのプログラムを記録したコンピュータ読み取り可能な媒体 Download PDFInfo
- Publication number
- JP4006098B2 JP4006098B2 JP15416698A JP15416698A JP4006098B2 JP 4006098 B2 JP4006098 B2 JP 4006098B2 JP 15416698 A JP15416698 A JP 15416698A JP 15416698 A JP15416698 A JP 15416698A JP 4006098 B2 JP4006098 B2 JP 4006098B2
- Authority
- JP
- Japan
- Prior art keywords
- paddy field
- water level
- amount
- water
- target
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
Images
Landscapes
- Control Of Non-Electrical Variables (AREA)
Description
【０００１】
【発明の属する技術分野】
本発明は、複数の水田から構成された水田地域において、収量の最大化を図るべく、水田からの地下浸透量を適正浸透量に保つように、給排水バルブ操作、揚水ポンプ操作その他の水管理作業を行うための水管理方法に係り、特に（ほぼ）同一土壌・土層、同一気候を有する水田地域に対して好適な、複数水田の水管理方法に関する。
【０００２】
【従来の技術】
従来、給排水パイプライン、給排水路の整備された水田地域では、稲の生育段階に応じた好ましい水位管理を行うため、各水田に１つの水位計を設置しておき、水位計測値が所定の目標水位に近づくように給排水バルブの操作を行っている。例えば、特開平９−６５７７６によれば、気象観測データ、気象予報データ、稲の生育データとに基づいて、目標水位テーブル及び該目標水位に対応した水位許容範囲テーブルを予め用意しておき、水位センサを介して水田の水位を検知しながら、該水田水位が目標水位に対して上記水位許容範囲内となるように、灌水装置の操作を行っている。
【０００３】
【発明が解決しようとする課題】
上記従来技術によれば、田植えから落水までの間、稲の生育段階に応じた適切な水田水位となるように、水田水位を自動的に変化させながら灌水を行うことができ、農家の人の灌水作業の省力化を図ることができる。また水位管理には、長年の経験が必要であったが、経験の少ない人でも適切な水田水位となるように水田の灌水を行うことができる。
【０００４】
ところで、「新編灌漑排水上巻（養賢堂）」によれば、水田水の地下浸透は、水に含まれる各種養分・酸素の根への供給、及び土壌の酸化促進等の役割を持ち、適正浸透量と呼ばれる水稲生育に好ましい浸透量が存在する。実際、ある地方の浸透量と水稲収量との関係の調査結果によれば、最低収量は減水深（蒸発散量＋浸透量）の増加とともに減少するが、最高収量は、減水深が２０〜３０mm/dayのときピークを持ち、それより多くても少なくても減少することが示されている。
【０００５】
但し、適正浸透量は水稲生育期間中一定ではなく、時期及び各種条件によって変化する性質を持つ。生育初期から最高分けつ中子期までは、土壌中の養分濃度も高く、土壌の還元による障害も少ないため、浸透量の増加は肥料分の流亡につながり、浸透は水稲収量に対してマイナスに働く。中干し後になると、夏期の高温のために有機物の分解が進み、土壌は著しい還元状態となり、各種の有害物質を生じるため、適度の透水と間断灌漑による土壌の酸化が必要となる。
【０００６】
即ち、植付より中干しまでは浸透を抑え、中干し後から落水までは透水と間断灌漑を行う、といったように、稲の生育段階に応じて適切な浸透量となるように目標浸透量管理を行うことが、高収量を得るために必要不可欠となる。
【０００７】
しかるに、従来技術では、浅水、深水、間断潅漑等のように水位を制御量とした目標水位管理を行っているのみで、水田からの浸透量を全く考慮してしない。浸透量は水田水位のみならず地下水位等にも影響されるため、偶然に水田浸透量が適正浸透量になる可能性はありうるが、従来技術では適正浸透量となるよう水田浸透量を能動的に操作できず、積極的な高収量化を図ることができない。
【０００８】
浸透量を制御量とした適正浸透量管理を行うためには、水田の浸透量推定を行い、推定浸透量が適正浸透量になるよう水位管理（給排水バルブ制御）を行うことか必要不可欠となる。しかるに、現在の技術レベルでは、浸透量を精度よく推定することは容易ではない。例えば、浸透量推定のためによく用いられている方法として、水田において成り立つ水収支式：水田貯留変化量＝給水量−排水量＋雨量−蒸発散量−浸透量、に対して、水田貯留変化量、給水量、排水量、雨量を測定し、蒸発散量を何らかの方法により推定することにより、残された浸透量を算出する方法がある。
【０００９】
しかしながら実際は、上記蒸発散量の推定が容易でないため、上記水収支式に基づく方法は有効とは言えない。傾度法、熱収支法等を用いれば蒸発散量推定可能であるが、特殊な測定器を必要とするので一般水田の水管理に用いることは現実的ではない。夜間は蒸発散量が０になることを利用すれば、上記方法でも浸透量の推定は可能となるが、昼間の浸透量推定には利用できず一般的ではない。
【００１０】
また、仮に水田の浸透量を推定できたとしても、推定浸透量が適正浸透量になるように水管理を行うことは容易ではない。浸透量は水田水位のみならず、地下水位、土壌・土層の状態等にも大きく左右されるため、適正浸透量を実現するための水田水位を決定することが困難なためである。
【００１１】
本発明の目的は、水稲収量に大きな影響を及ぼす物理量として浸透量に着目し、水田群の水位トレンド、給排水・雨量情報のみから、水田の土壌・土層、地下水位等の状態に依存することなく水田の浸透量を推定し、各水田の浸透量が稲の生育に応じた適正浸透量となるよう水田水位を適切に管理することができる、水田水管理の技術を提供することにある。
【００１２】
【課題を解決するための手段】
上記目的達成のため、次の手段を設ける。
（手段１）
複数区画からなる水田に対する水管理を行う水田群の水管理システムにおいて、上記水田の存在する水田地域の雨量を把握する手段と、上記各水田の水位を計測する水位計と、各水田に給水および排水を行う給排水装置と、各水田の浸透量の目標値およびその許容範囲を決定する目標浸透量決定手段と、上記把握された雨量および各水田の水位に基づいて、所定アルゴリズムに従って各水田の浸透量を推定する浸透量推定手段と、上記推定された浸透量の目標値の許容範囲に収まるように前記給排水装置を動作させる制御手段とを有することを特徴とする。
【００１３】
（手段２）
複数区画からなる水田に対する水管理を行う水田群の水管理システムにおいて、上記水田の存在する水田地域の雨量を把握する手段と、上記各水田の水位を計測する水位計と、各水田に給水および排水を行う給排水装置と、各水田の浸透量の目標値およびその許容範囲を決定する目標浸透量決定手段と、上記把握された雨量および各水田の水位に基づいて、所定アルゴリズムに従って各水田の浸透量を推定する浸透量推定手段と、上記推定された浸透量の目標値の許容範囲に収まるように各水田の目標水位及びその許容範囲を決定する目標水位決定手段とを有することを特徴とする。
【００１４】
（手段３）
本発明の水田の水管理システムでは、手段１、２における所定のアルゴリズムは、各水田に対して成り立つ水収支式において、蒸発散量、及び地下水位、透水係数等の浸透に影響を及ぼす物理量が上記各式とも共通の値を取ると仮定して、上記水収支式を連立して解くことにより各水田の浸透量を算出すること、を特徴とする。
【００１５】
（手段４）
本発明の水田の水管理システムでは、手段１、２における所定の目標量及びその許容範囲は、稲の生育段階又は日付と、該当希稲種の目標浸透量及びその許容範囲との対応関係を記録したテーブルに基づいて決定すること、を特徴とする。
【００１６】
（手段５）
本発明の水田の水管理システムでは、水田地域の気温、水温、湿度、日射豊、日照時間等の気象情報を把握するための手段を有し、手段１、２における所定の目標量及びその許容範囲は、稲の生育段階又は日付と、該当稲種の目標浸透量及びその許容範囲との対応関係を記録したテーブルに基づいて決定し、上記目標浸透量又はその許容範囲は、上記把握された気象情報に基づいて所定量だけ修正すること、を特徴とする。
【００１７】
（手段６）
本発明の水田の水管理システムでは、手段２における目標水位決定手段は、各水田に対して成り立つ水収支式において、蒸発散量、及び地下水位、透水係数等の浸透に影響を及ぼす物理量が上記各式とも共通の値を取ると仮定し、上記水収支式を連立して解くことにより上記物理量の値を算出し、上記算出された物理量の値及び前記決定された浸透量の目標値及び浸透量の上下限値を上記各水田ごとの水収支式に代入して解くことにより、前記各水田の目標水位及びその許容範囲を決定すること、を特徴とする。
【００１８】
（手段７）
本発明の水田の水管理システムでは、複数区画からなる水田地域において、上記水田地域の雨量を把握するための手段と、上記各水田の水位を計測するための水位計と、上記把握された雨量及び各水田の水位に基づいて、所定のアルゴリズムに従って各水田の水収支に関連する物理量を推定する物理量推定手段と、を有することを特徴とする。
【００１９】
（手段８）
本発明の水田の水管理システムでは、手段７における水収支に関連する物理量は、浸透量、地下水位、透水係数、その他浸透に影響を及ぼす物理量、蒸発散量、給水量、排水量のいずれかであること、を特徴とする。
【００２０】
（手段９）
本発明の水田の水管理システムでは、手段７における所定のアルゴリズムは、各水田に対して成り立つ水収支式において、蒸発散量、及び地下水位、透水係数等の浸透に影響を及ぼす物理量が上記各式とも共通の値を取ると仮定して、上記水収支式を連立立して解くことにより各水田の任意の期間中の水の移動に関連する物理量を算出すること、を特徴とする。
【００２１】
（手段１０）
本発明の水田の水管理システムでは、複数区画からなる水田地域において、上記水田地域の雨量を把握するための手段と、上記各水田の水位を計測するための水位計と、各水田に給水又は排水を行うための給排水装置と、各水田の水位変化量の目標値及びその許容範囲を決定する目標水位変化量決定手段と、上記把握された雨量及び各水田の水位に基づいて、所定のアルゴリズムに従って降雨の影響又は給排水の影響又は蒸発散の影響を取り除いた各水田の水位変化量を推定する水位変化量推定手段と、上記推定された水位変化量が上記水位変化量の目標値の許容範囲内に収まるように各水田の給排水装置を動作させる制御手段と、を有することを特徴とする。
【００２２】
（手段１１）
本発明の水田の水管理システムでは、複数区画からなる水田地域において、上記水田地域の雨量を把握するための手段と、上記客水田の水位を計測するための水位計と、各水田に給水又は排水を行うための給排水装置と、各水田の水位変化量の目標位及びその許容範囲を決定する目標水位変化量決定手段と、上記把握された雨量及び各水田の水位に基づいて、所定のアルゴリズムに従って降雨の影響又は給排水の影響又は蒸発散の影響を取り除いた各水田の水位変化量を推定する水位変化量推定手段と、上記推定された水位変化量が上記水位変化量の目標値の許容範囲内に収まるように各水田の目標水位及びその許容範囲を決定する目標水位決定手段と、を有することを特徴とする。
【００２３】
（手段１２）
本発明の複数区画からなる水田地域における水田の水管理を行うためのコンピュータ読み取り可能な媒体では、上記水画地域の雨量を把握する機能と、上記客水田に設置された水位計からの水位計測値を把握する機能と、各水田の浸透量の目標値及びその許容範囲を決定する目標浸透量決定機能と、上記把握された雨量及び各水田の水位に基づいて、所定のアルゴリズムに従って各水田の浸透量を推定する浸透量推定技能と、上記推定された浸透量が上記浸透量の目標値の許容範囲内に収まるように各水田に設置された給排水装置を動作させる制御機能と、を実行させるためのプログラムを記録したことを特徴とする。
【００２４】
（手段１３）
本発明の複数区画からなる水田地域における水田の水管理を行うためのコンピュータ読み取り可能な媒体では、上記水田地域の雨量を把握する機能と、上記客水田に設置された水位計からの水位計測値を把握する機能と、各水田の浸透量の目標値及びその許容範囲を決定する目標浸透量決定機能と、上記把握された雨量及び各水田の水位に基づいて、所定のアルゴリズムに従って各水田の浸透量を推定する浸透量推定機能と、上記推定された浸透量が上記浸透量の目標値の許容範囲内に収まるように各水田の目標水位及びその許容範囲を決定する目標水位決定機能と、を実行させるためのプログラムを記録したことを特徴とする。
【００２５】
（手段１４）
本発明の複数区画からなる水田地域における水田の水管理を行うためのコンピュータ読み取り可能な媒体では、手段１２、１３における所定のアルゴリズムは、各水田に対して成り立つ水収支式において、蒸発散量、及び地下水位、透水係数等の浸透に影響を及ぼす物理量が上記各式とも共通の値を取ると仮定して、上記水収支式を連立して解くことにより各水田の浸透量を算出すること、を特徴とする。
【００２６】
（手段１５）
本発明の複数区画からなる水田地域における水田の水管理を行うためのコンピュータ読み取り可能な媒体では、手段１２、１３における所定の目標量及びその許容範囲は、稲の生育段階又は日付と、該当稲種の目標浸透量及びその許容範囲との対応関係を記録したテーブルに基づいて決定すること、を特徴とするとする。
【００２７】
（手段１６）
本発明の複数区画からなる水田地域における水田の水管理を行うためのコンピュータ読み取り可能な媒体では、上記水田地域の気温、水温、湿度、日射量、日照時間等の気象情報を把握する機能を実行させるためのプログラムを記録し、手段１２、１３における所定の目標量及びその許容範囲は、稲の生育段階又は日付と、該当稲種の目標浸透量及びその許容範囲との対応関係を記録したテーブルに基づいて決定し、上記目標浸透量又はその許容範囲は、上記把握された気象情報に基づいて所定量だけ修正すること、を特徴とする。
【００２８】
（手段１７）
本発明の複数区画からなる水田地域における水田の水管理を行うためのコンピュータ読み取り可能な媒体では、手段１３における目標水位決定機能は、各水田に対して成り立つ水収支式において、蒸発散量、及び地下水位、透水係数等の浸透に影響饗を及ぼす物理量が上記各式とも共通の値を取ると仮定し、上記水収支式を連立して解くことにより上記物理量の値を算出し、上記算出された物理量の値及び前記決定された浸透量の目標値及び浸透量の上下限値を上記各水田ごとの水収支式に代入して解くことにより、前記各水田の目標水位及びその許容範囲を決定すること、を特徴とする。
【００２９】
（手段１８）
本発明の複数区画からなる水田地域における水田の水管理を行うためのコンピュータ読み取り可能な媒体では、上記水田地域の雨量を把握する機能と、上記各水田に設置された水位計からの水位計測値を把握する機能と、上記把握された雨量及び各水田の水位に基づいて、所定のアルゴリズムに従って各水田の水収支に関連する物理量を推定する物理量推定機能と、を実行させるためのプログラムを記録したことを特徴とする。
【００３０】
（手段１９）
本発明の複数区画からなる水田地域における水田の水管理を行うためのコンピュータ読み取り可能な媒体では、手段１８における水収支に関連する物理量は、浸透量、地下水位、透水係数、その他浸透に影響を及ぼす物理量、蒸発散量、給水量、排水量のいずれかであることを特徴とするコンピュータ読み取り可能な媒体。
【００３１】
（手段２０）
本発明の複数区画からなる水田地域における水田の水管理を行うためのコンピュータ読み取り可能な媒体では、手段１８における所定のアルゴリズムは、各水田に対して成り立つ水収支式において、蒸発散量、及び地下水位、透水係数等の浸透に影響を及ぼす物理量が上記各式とも共通の値を取ると仮定して、上記水収支式を連立して解くことにより各水田の任意の期間中の水の移動に関連する物理量を算出すること、を特徴とする。
【００３２】
（手段２１）
本発明の複数区画からなる水田地域における水田の水管理を行うためのコンピュータ読み取り可能な媒体では、上記水田地域の雨量を把握する機能と、上記各水田に設置された水位計からの水位計測値を把握する機能と、各水田の水位変化量の目標値及びその許容範囲を決定する目標水位変化量決定機能と、上記把握された雨量及び各水田の水位に基づいて、所定のアルゴリズムに従って降雨の影響又は給排水の影響又は蒸発散の影響を取り除いた各水田の水位変化量を推定する水位変化量推定機能と、上記推定された水位変化量が上記水位変化量の目標値の許容範囲内に収まるように各水田に設置された給排水装置を動作させる制御機能と、を実行させるためのプログラムを記録したことを特徴とする。
【００３３】
（手段２２）
本発明の複数区画からなる水田地域における水田の水管理を行うためのコンピュータ読み取り可能な媒体では、上記水田地域の雨量を把握する機能と、上記各水田に設置された水位計からの水位計測値を把握する機能と、各水田の水位変化量の目標値及びその許容範囲を決定する目標水位変化量決定機能と、上記把握された雨量及び各水田の水位に基づいて、所定のアルゴリズムに従って降雨の影響又は給排水の影響又は蒸発散の影響を取り除いた各水田の水位変化量を推定する水位変化量推定機能と、上記推定された水位変化量が上記水位変化量の目標値の許容範囲内に収まるように各水田の目標水位及びその許容範囲を決定する目標水位決定機能と、を実行させるためのプログラムを記録したことを特徴とする。
【００３４】
【発明の実施の形態】
以下、本発明を図面に示す実施形態によって詳細に説明する。
図１は、水田群における各水田の水位を管理する水管理システムの、実施形態を示すシステム構成図である。本実施形態の水管理システムは、水位計１−ｊ（ｊ＝１、…、３０）、給水バルブ２−ｊ、排水バルブ３−ｊ、コンピュータ１０４、データ伝送路１０５、気象情報提供システム１０６、気象ロボット１０７、給水パイプライン１０８、排水パイプライン１０９から構成される。
【００３５】
水管理の対象となる水田地域は、農水省が推進している大区画水田モデル地区と同様に、３０枚の１ｈａ規模の水田４−ｊによって構成される。各水田４−ｊ内には、水位計１−ｊ、給水バルブ２−ｊ、排水バルブ３−ｊが１つづつ設置されており、所定計測周期ごとに各水田４−ｊの水位、給水量、排水量をデータ伝送路１０５を介してコンピュータ１０４に送信している。ここで、データ伝送路１０５は、無線又は有線のどちらでも構わない。上記水田地域には、気象ロボット１０７が設置されており、上記水位計測周期と同じタイミングで、上記水田地域における気温、湿度、日照時間、日射量、風向・風速、雨量を、データ伝送路１０５を介してコンピュータ１０４に送信している。
【００３６】
更に、上記気象データ（気温、湿度、日照時間、日射量、風向・風速、雨量）の所定周期、所定時刻先までの予報値を、日本気象協会、又は民間団体の保有する気象情報提供システム１０６から、所定タイミングごとに、データ伝送路１０５を介してコンピュータ１０４に送信している。
【００３７】
コンピュータ１０４は、図２に示すように、ＣＰＵ２２1、ＲＡＭ２２２、ハードディスク装置２２３、フロッピーディスク装置２２４、ＭＯディスク装置２２５、ＤＶＤ装置２２６、ＣＤ−ＲＯＭ装置２２７、ＲＡＭカード装置２２８、入力装置２２９、表示装置２３０によって構成される。ＣＰＵ２２１は、水管理システム全体の動作を制御して、水稲収量に大きな影響を及ぼす物理量である浸透量を制御量として、給水バルブ２−ｊ、排水バルブ３−ｊの操作を行い、各水田４−ｊの水管理を行う中央処理装置である。ＲＡＭ２２２は、各種処理プログラムや各種データをロードする記憶装置である。
【００３８】
ハードディスク装置２２３は、各種処理プログラム２０１〜２０５、及び各種データベース２１１、２１２を磁気ディスクに格納する記憶装置である。フロッピーディスク装置２２４は、各種処理プログラム２０１〜２０５、及び各種データベース２１１、２１２が記録されたフロッピーディスクの読み書きを行う装置である。ＭＯディスク装置２２５は、各種処理プログラム２０１〜２０５、及び各種データベース２１１、２１２が記録されたＭＯディスク（光磁気ディスク）の読み書きを行う装置である。
【００３９】
ＤＶＤ装置２２６は、各種処理プログラム２０１〜２０５、及び各種データベース２１１、２１２が記録されたＤＶＤ（デジタルビデオディスク）の読み書きを行う装置である。ＣＤ−ＲＯＭ装置２２７は、各種処理プログラム２０１〜２０５、及び各種データベース２１１、２１２が記録されたＣＤ−ＲＯＭ（コンパクトディスク）の読み書きを行う装置である。ＲＡＭカード装置２２７は、各種処理プログラム２０１〜２０５、及び各種データベース２１１、２１２が記録されたＲＡＭカードの読み書きを行う装置である。
【００４０】
入力装置２２９は、水田浸透量推定プログラム２０３実行に必要となる各水田４−ｊの物理形状データや、稲の生育段階推定処理プログラム２０１実行に必要となる各水田４−ｊの稲の品種や田植え日等を入力するための装置である。表示装置２３０は、各種処理プログラム２０１〜２０５によって算出された各水田４−ｊの目標浸透量、水田浸透量、目標水位、バルブ開閉状態／及び稲の生育段階等、並びに気象情報提供システム１０６から送信された上記水田地域の気象データ（気温、湿度、日照時間、日射量、風向・風速、雨量）の予報値、並びに気象ロボット１０７、水位計１−ｊ、給水バルブ２−ｊ、排水バルブ３−ｊから送信された各水田４−ｊの水位、給水量、及び排水量等の表示を行う装置である。
【００４１】
ここで、上記各種処理プログラム２０１〜２０５、及び各種データベース２１１、２１２は上記記憶媒体のいずれに記録していてもかまわないため、以下、記録媒体の読み書きを行う装置がハードディスク装置２２３のみの場合を例にとって説明する。コンピュータ１０４では、各水田４−ｊの浸透量が稲の生育に応じた適正浸透量となるように、水田水位の管理を行う。上記水田水位管理は、ハードディスク装置２２３に記憶された各種処理プログラム２０１〜２０５、及び各種データベース２１１、２１２をＲＡＭ２２２にロード・記憶し、ＣＰＵ２２１によって、各種処理プログラム２０１〜２０５に対応する以下の処理：
（１）各水田４−ｊの現在の稲の生育段階を推定する処理
（２）各水田４−ｊの上記生育段階に応じた適切な目標浸透量を決定する処理
（３）各水田４−ｊの浸透量を推定する処理
（４）各水田４−ｊの上記目標浸透量に対応、する目標水田水位を決定する処理
（５）各水田４−ｊの水位を上記目標水田水位の所定の近傍内に収めるよう、給水バルブ２−ｊ、排水バルブ３−ｊの自動制御を行う処理
を順次実行していくことにより実現される。
【００４２】
以下、ＣＰＵ２２１において実行される上記処理（１）〜（５）の実行方法について説明する。はじめに、生育段階推定処理プログラム２０１によって実現される、各水田４−ｊの現在の稲の生育段階を推定する処理（１）について、図２、図３を用いて説明する。コンピュータ１０４では、気象ロボット１０７より計測された気温日照時間のデータを、所定周期ごとに受信している。処理プログラム２０１では、１日１回日没後、上記所定周期ごとの気温、日照時間データから当日の昼間時の平均気温、日照時間を算出する。例えば、昼間時の平均気温は、昼間時における上記所定周期ごとの気温の平均を取ることにより算出できる。また、昼間時の日照時間は、昼間時における上記所定周期ごとの日照時間を積算することにより算出できる。
【００４３】
生育モデルデータベース２１１では、稲の品種別の生育段階推定モデルを配慮している。上記生育段階推定モデルとしては、国内において広く用いられているＤＶＩモデルを採用する。ＤＶＩモデルによれば、田植え後ｉ日目の生育値ＤＶＩは、出芽日からその日までの生育速度ＤＶＲを積算して、
【００４４】
【数１】
【００４５】
と記述される。ここで、図３に示すように、生育速度ＤＶＲは気温が高いほど、また短日条件ほど大きくなることが知られている。
【００４６】
処理プログラム２０１では、入力装置２２９より入力された各水田４−ｊの稲の品種に基づいて、上記データベース２１１から対応する生育段階推定モデルを選択する。更に、入力装置２２９より入力された各水田４−ｊの稲の田植え日、及び処理プログラム２０１にて上記算出した田植え日から当日までの昼間時における平均気温、日照時間を（数１）に代入することにより、各水田４−ｊの当日までの稲の生育段階（ＤＶＩ値）を算出する。上述のように処理プログラム２０１において、自動的に稲の生育段階を算出してもよいが、もちろん農家において稲の生育状態を観察し、生育段階（ＤＶＩ値）を処理２０１に代わって入力装置２２９からコンピュータ１０４に入力してもよい。上述のようにして、ＣＰＵ２２１において、処理プログラム２０１により、各水田４−ｊの当日までの稲の生育段階の推定が行われる。
【００４７】
次に、目標浸透量決定処理理プログラム２０２によって実現される、現在の稲の生育段階に応じた適切な目標浸透量を決定する処理（２）について、図２、図４を用いて説明する。目標浸透量データベース２１２では、稲の生育段階（ＤＶＩ値）に対応した、稲の品種別・地方別の、田植え日から収穫までの、目標浸透量テーブル、及び当該目標浸透量データに対する許容範囲テーブルを記憶している。
【００４８】
図４に、目標浸透量テーブル、及び許容範囲テーブルの一例を示す。田植え直後では、土壌の養分濃度も高く、土壌の還元障害も少ないため、浸透量を低めに抑え、養分の流亡を防ぐ。有効分けつ期では地温上昇のため浸透を増やし、無効分けつ期では、分けつを抑制するため浸透を減らし、中干しを行う。
【００４９】
幼穂形成期、穂ばらみ期を通して豊熟期の前半までは、夏期の高温のために有機物の分解が盛んで土壌は著しい還元状態となり、各種の有害物質を生じるため、適度の浸透と土壌の酸化を目的とした間断灌漑を行う。収穫期では、出穂２５〜３０日後を目安に落水を行い、浸透を０にし、土壌を固め刈り取りに備える。
【００５０】
上記目標浸透量テーブル、及び許容範囲テーブルは、過去の稲の生育データや水位・浸透量データ、気象データ等に基づいて、稲の生育段階（ＤＶＩ値）から目標浸透量、及び許容範囲が特定できるよう「低温時には目標浸透量を所定量だけ増加させる」といったノウハウを組み込んでおく。また、上記テーブルは固定ではなく、農家のオリジナリティを組み込めるよう、いつでも変更可能なフレキシブルなテーブルとする。
【００５１】
処理プログラム２０２では、処理プログラム２０１の実行周期と同じタイミングで、以下の処理を行う。処理プログラム２０２では、入力装置２２９より入力された各水田４−ｊの稲の品種・地方に基づいて、上記データベース２１２から対応する目標浸透量テーブル、及び許容範囲テーブルを選択する。更に、上記処理プログラム２０１より算出された各水田４−ｊの稲の生育段階と、上記選択された各水田４−ｊの目標浸透量テーブル、及び許容範囲テーブルとを比較して、各水田４−ｊの適切な目標浸透量、及び許容範囲を決定する。決定された目標浸透量、及び許容範囲は、処理プログラム２０２の今回実行時から次回実行時まで有効とする。
【００５２】
また、コンピュータ１０４では、気象情報提供システム１０６より、所定周期、所定時刻先までの気象予報データを、所定タイミングごとに受信している。上記目標浸透量テーブルに、例えば上記の「低温時には目標浸透量を所定量だけ増加させる」といったノウハウが対応している場合は、処理プログラム２０２では、決定された目標浸透量の有効期間（当日の所定時刻〜次の日の該所定時刻）における気温予報データを考慮し、決定された目標浸透量の修正を行う。上述のようにして、ＣＰＵ２２１において、処理プログラム２０２により、各水田４−ｊの目標浸透量、及びその許容範囲の決定が行われる。
【００５３】
次に、水田浸透量推定処理プログラム２０３によって実現される、各水田４−ｊの浸透量を推定する処理（３）について、図２、図５を用いて説明する。コンピュータ１０４では、水位計１−ｊ、給水バルブ２−ｊ、排水バルブ３−ｊより計測された各水田４−ｊの水位、給水量、排水量、及び気象ロボット１０７より計測された当該水田地域の雨量を、所定周期ごとに受信している。処理プログラム２０３では、１日１回乃至数回（例えば３時間ごと）、上記受信データを用いて、各水田４−ｊの現在の浸透量の推定を行う。土壌中の浸透量は、Darcyの法則「浸透速度は動水勾配に正比例する」に支配され、水田４−ｊにおける浸透量（浸透速度）Ｖｊは、閉鎖浸透時（負圧の土層が存在しない浸透状態）、
【００５４】
【数２】
【００５５】
として記述できる。ここで、土壌条件の違いにより、浸透には閉鎖浸透と開放浸透の２つの状態が存在しうるが、開放浸透時でも（数２）と同様の式が成り立ち、開放浸透時に成り立つ式の方が未知変数も少なく単純であるため、より複雑な閉鎖浸透時に的を絞って説明を進める。
【００５６】
ここで、各水田４−ｊの水の移動に関連する物理量を図５に示す記号で表すことにする。通常、同一水田地域内の各水田においては、ほぼ同一の土壌・土層を有しており、またほほ同一の地下水位を有していることから、上記（数２）における合成透水係数Ｋ、全土層厚Ｌ、地下水位Ｈ(t）は各水田４−ｊにおいて共通であると仮定できる。また同一水田地域内の各水田においては、ほほ同一の気候、気温、湿度、日射、風速等）条件下にあることから、各水田４−ｊの（単位面積当りの）蒸発散量Ｅは共通であると仮定できる。上述のことから、現在までの所定期間△Ｔ（前回浸透量推定実行時〜今回浸透量推定実行時）中の各水田４−ｊにおける水収支式は、
【００５７】
【数３】
【００５８】
として記述できる。
【００５９】
連立方程式（数３）において、未知変数は、Ｋ、Ｌ、Ｈ（t）、Ｈ（ｔ−△Ｔ）、Ｅ（t）の５個であり、その他の変数は、コンピュータ１０４が計測している物理量から一意かつ容易に算出できる。従って、未知変数の個数と同じ５枚以上の水田があれば連立方程式（数３）は求解可能となり、本実施形態では水管理の対象水田を３０枚としているため十分求解可能となる。
【００６０】
処理プログラム２０３では、所定周期ごとに計測している各種物理量から、連立方程式（数３）を構築し、適当な求解手法に基づいて連立方程式（数３）の求解を行い、各水田４−ｊの浸透量を推定する。但し、浸透量推定期間△Ｔ中の水田水位、及び給排水量が完全に一致する水田が複数枚存在する場合は、該水田に対応する（数３）の方程式が一致してしまい、連立方程式（数３）が不定になり求解できない。
【００６１】
この場合、処理２０３プログラムでは、上記水田に対して、後述する目標水田水位の許容範囲内に水田水位が収まるよう給水又は排水を行い、連立方程式（数３）が不定になるのを防止する。上述のようにして、ＣＰＵ２２１において、処理プログラム２０３により、各水田４−ｊの浸透量の推定が行われる。
【００６２】
ところで、処理プログラム２０３では、各水田４−ｊの給水量及び排水量を計測するものと仮定しているが、以下の方法を用いることにより上記給排水量は計測不要となり、更に上記給排水量は推定可能とさえなる。即ち、連立方程式（数３）を構成する水収支式を、給排水を行っていない水田に対応する水収支式（給排水量の項＝０）に限定して該連立方程式（数３）の求解を行い、各水田４−ｊの浸透量、合成透水係数Ｋ、全土層厚Ｌ、地下水位（t）、蒸発散量Ｅを推定する。上記推定値を、給排水を行った水田に対応する水収支式（数３）に代入することにより、給排水量の値が推定可能となる。なお、給排水を行っていない水田の存在は、各水田４−ｊへの給排水の順序を適当に割り振ることにより容易に実現可能となる。
【００６３】
ここで、本発明の大きな特徴である、上述の水田群の浸透量推定方法の基本的考え方について説明しておく。従来、１水田における水収支式に基づいて浸透量を推定することは、方程式の個数より未知数（蒸発散量、地下水位、透水係数）の個数の方が多いため、求解不可能であった。しかしながら本発明では、対象水田を、同一地域内に存在する複数枚の水田までに拡大することにより、その求解を可能にしている。即ち、対象水田を複数枚に増加させることで上記方程式の個数を増加させ、更に、対象水田か何一地域内に存在することから、各水田は同一土壌条件下、同一気象条件下にあると仮定し、上記未知変数（蒸発散量、地下水位、透水係数）は共通の値をとるとして、対象水田の増加に伴う未知変数の増加を防止している。これにより、方程式の個数が未知数の個数を上回るようになり、その求解を可能としている。
【００６４】
次に、目標水位決定処理プログラム２０４によって実現される、各水田４−ｊの目標浸透量に対応する目標水田水位を決定する処理（４）について、図２、図５を用いて説明する。処理プログラム２０４では、処理プログラム２０３の実行周期（浸透量推定周期）と同じタイミングで、処理プログラム２０２にて算出した各水田４−ｊの目標浸透量、及び許容範囲に対応する、各水田４−ｊの目標水位、及び許容範囲を決定する。処理プログラム２０３の実行結果から、図５に示す当該水田地域の土層の合成透水係数Ｋ、全土層厚Ｌ、地下水位Ｈ（t）は全て算出されているため、上記値を（数２）に代入すれば、（数２）は水田水位と水田浸透量との１対１の関係を示す式となる。従って、上記各水田４−ｊにおける３つの浸透量（目標浸透量、上限浸透量、下限浸透量）を上記関係式に代入することにより、各水田４−ｊの目標水位、及び許容範囲を決定する。決定された目標水位、及び許容範囲は、処理２０４の今回実行時から次回実行時まで有効とする。上述のようにして、ＣＰＵ２２１において、処理プログラム２０４により、各水田４−ｊの目標未位、及びその許容範囲の決定が行われる。
【００６５】
次に、給排水バルブ制御処理プログラム２０５によって実現される、各水田４−ｊの水位を上記目標水田水位の所定の近傍内に収めるよう、給水バルブ２−ｊ、排水バルブ３−ｊの自動制御を行う処理（５）について、図２、図６を用いて説明する。コンピュータ１０４では、水位計１−ｊより計測された各水田４−ｊの水位を、所定周期ごとに受信している。処理プログラム２０５では、上記受信データを用いて、上記水位計測周期と同じタイミングで、各給水バルブ２−ｊ、排水バルブ３−ｊの開閉状態を決定し、各バルブに対して開閉信号の送信を行う。
【００６６】
図６は、処理プログラム２０５が実行するバルブ制御処理の概要を示すフローチャートである。はじめに、各水位計１−ｊより計測した最新の水田水位を読み込み（ステップ６０１）、また各水田４−ｊの目標水位、及びその許容範囲を、処理プログラム２０４において決定された最新値に更新する（ステップ６０２）。
【００６７】
次に処理プログラム２０５では、各水田４−ｊの水位が目標水位の所定の近傍Ｉ内（目標水位±許容範囲）に存在するように、以下に示す給排水バルブの開放／閉鎖条件に基づいて、給水バルブ２−ｊ、排水バルブ３−ｊの操作を行う。
【００６８】
・給水バルブ開放条件（閉鎖から開放に移行する条件）：水田水位＜目標水位−許容範囲
・給水バルブ閉鎖条件（開放から閉鎖に移行する条件）：水田水位＞目標水位
・排水バルブ開放条件：水田水位＞目標水位＋許容範囲
・排水バルブ閉鎖条件：水田水位＜目標水位
即ち、水田水位が近傍Ｉの下限値（目標水位−許容範囲）を下回っていればデータ伝送路１０５を介して、給水バルブ２−ｊに開放信号を送信し、給水を実行させる（ステップ６０３、６０４）。給水バルブ２−ｊでは、コンピュータ１０４からのバルブ開放信号を受信し、バルブを開けて給水パイプライン１０８から水田４−ｊ内に水を供給する。これにより水田４−ｊの水位が漸次上昇し、水田水位が目標水位を上回るようになればデータ伝送路１０５を介して、給水バルブ２−ｊに閉鎖信号を送信し、給水を停止させる（ステップ６０５、６０６）。これにより水田４−ｊの水位は目標水位付近で停止することになる。
【００６９】
逆に、水田水位が近傍Ｉの上限値（目標水位＋許容範囲）を上回っていればデータ伝送路１０５を介して、排水バルブ３−ｊに開放信号を送信し、排水を実行させる（ステップ６０７、６０８）。排水バルブ３−ｊでは、コンピュータ１０４からのバルブ開放信号を受信し、バルブを開けて水田４−ｊから排水パイプライン１０９内に水を排水する。これにより水田４−ｊの水位が漸次下降し、水田水位が目標水位を下回るようになれば、データ伝送路１０５を介して、排水バルブ３−ｊに閉鎖信号を送信し、排水を停止させる（ステップ６０９、６１０）。これにより水田４−ｊの水位は目標水位付近で停止することになる。上述のようにして、ＣＰＵ２２１において、処理プログラム２０５により、水田水位が目標水位の所定の近傍Ｉ内に存在するように、給水バルブ２−ｊ、及び排水バルブ３−ｊの自動制御が行われる。
【００７０】
上述のように、本発明の実施形態によれば、水稲収量に大きな影響を及ぼす物理量として浸透量に着目し、同一地域内の水田群では、各水田とも同一土壌下、同一気候条件下にあるという物理現象を利用して、水田群の水位トレンド、給排水、雨量情報から、水田の土壌・土層、地下水位等の状態に依存することなく水田の浸透量を推定し、各水田の浸透量が稲の生育に応じた適正浸透量なるよう水田水位を適切に管理することが可能である。
【００７１】
以上説明したように本発明の実施の形態によれば各水田において成り立つ水収支式に基づいて、地下水位、透水係数、蒸発散量等を計測不要な未知変数として浸透量推定計算を行っているため、水田の土壌・土質、地下水位、蒸発散量を把握するための特殊な計器を用いることなく、また昼夜の別に限定されることなく、各水田の現在の浸透量を推定できる。
【００７２】
また、稲の生育段階推定結果と、稲の生育段階に対応した目標浸透量テーブルに基づいて、各水田の適正浸透量の目標値を決定し、上記現推定浸透量との比較を行っているため、各水田において現在の稲の生育に応じた適正浸透が行われているか否かを判断できる。
【００７３】
適正浸透が行われていない水田に対しては、上記浸透量推定計算において、浸透に影響を及ぼす物理量（地下水位、透水係数）をも同時に把握しているため、適正浸透量を実現するのに過不足ない目標水田水位を算出でき、更に水田水位が上記目標水田水位に等しくなるよう給排水バルブの操作を行っているため、各水田の浸透量が適正浸透量になるような水管理を行うことができる。
【００７４】
且つ、上記浸透量推定周期、即ち水管理周期を任意の周期（例えば３時間程度の短期間）に設定できるため、各水田の浸透量を稲の生育期間中常に適正浸透量に維持することができ、より緻密な水管理を行うことができる。
【００７５】
即ち、上述の効果を総合すれば、本発明によれば、稲の植付から刈り取りまでの全生育期間に対して、各水田の浸透量を、現在の稲の生育に応じた適正浸透量に絶え間なく維持することができるため、水稲収量の最大化を図ることができる。
【００７６】
また、水田の土壌・土層の状態に応じて水田からの浸透特性は大きく違ってくるため、全国には漏水田や排水不良田など様々な水田が存在するが、本発明では、上記浸透量推定計算において、土壌・土層の状態を表すパラメータである透水係数をオンライン推定し、該透水係数に基づいて、水田の適正浸透量管理を行っているため、任意の水田地域に対して適用可能な、汎用性の高い水管理システムを構築できる。
【００７７】
【発明の効果】
以上の構成により、本発明においては、特殊な計器を用いなくとも、各水田の水の浸透量を把握することが可能になる。
【図面の簡単な説明】
【図１】本発明の水管理システムの実施形態を示すシステム構成図である。
【図２】コンピュータ内部のハード構成図及び各種処理プログラムのブロック図である。
【図３】平均気温、日照時間と稲の生育速度との関係を示す説明図である。
【図４】稲の生育段階に応じた目標水深テーブルの一例を示す説明図である。
【図５】１つの水田における水収支の関係を示す説明図である。
【図６】給排水バルブ制御の概要を示すフローチャートである。
【符号の説明】
１−ｊ…水位計、２−ｊ…給水バルブ、３−ｊ…排水バルブ、４−ｊ…水田、１０４…コンピュータ、１０５…データ伝送路、１０６…気象情報提供システム、１０７…気象ロボット、１０８…給水パイプライン、１０９…排水パイプライン。
【発明の属する技術分野】
本発明は、複数の水田から構成された水田地域において、収量の最大化を図るべく、水田からの地下浸透量を適正浸透量に保つように、給排水バルブ操作、揚水ポンプ操作その他の水管理作業を行うための水管理方法に係り、特に（ほぼ）同一土壌・土層、同一気候を有する水田地域に対して好適な、複数水田の水管理方法に関する。
【０００２】
【従来の技術】
従来、給排水パイプライン、給排水路の整備された水田地域では、稲の生育段階に応じた好ましい水位管理を行うため、各水田に１つの水位計を設置しておき、水位計測値が所定の目標水位に近づくように給排水バルブの操作を行っている。例えば、特開平９−６５７７６によれば、気象観測データ、気象予報データ、稲の生育データとに基づいて、目標水位テーブル及び該目標水位に対応した水位許容範囲テーブルを予め用意しておき、水位センサを介して水田の水位を検知しながら、該水田水位が目標水位に対して上記水位許容範囲内となるように、灌水装置の操作を行っている。
【０００３】
【発明が解決しようとする課題】
上記従来技術によれば、田植えから落水までの間、稲の生育段階に応じた適切な水田水位となるように、水田水位を自動的に変化させながら灌水を行うことができ、農家の人の灌水作業の省力化を図ることができる。また水位管理には、長年の経験が必要であったが、経験の少ない人でも適切な水田水位となるように水田の灌水を行うことができる。
【０００４】
ところで、「新編灌漑排水上巻（養賢堂）」によれば、水田水の地下浸透は、水に含まれる各種養分・酸素の根への供給、及び土壌の酸化促進等の役割を持ち、適正浸透量と呼ばれる水稲生育に好ましい浸透量が存在する。実際、ある地方の浸透量と水稲収量との関係の調査結果によれば、最低収量は減水深（蒸発散量＋浸透量）の増加とともに減少するが、最高収量は、減水深が２０〜３０mm/dayのときピークを持ち、それより多くても少なくても減少することが示されている。
【０００５】
但し、適正浸透量は水稲生育期間中一定ではなく、時期及び各種条件によって変化する性質を持つ。生育初期から最高分けつ中子期までは、土壌中の養分濃度も高く、土壌の還元による障害も少ないため、浸透量の増加は肥料分の流亡につながり、浸透は水稲収量に対してマイナスに働く。中干し後になると、夏期の高温のために有機物の分解が進み、土壌は著しい還元状態となり、各種の有害物質を生じるため、適度の透水と間断灌漑による土壌の酸化が必要となる。
【０００６】
即ち、植付より中干しまでは浸透を抑え、中干し後から落水までは透水と間断灌漑を行う、といったように、稲の生育段階に応じて適切な浸透量となるように目標浸透量管理を行うことが、高収量を得るために必要不可欠となる。
【０００７】
しかるに、従来技術では、浅水、深水、間断潅漑等のように水位を制御量とした目標水位管理を行っているのみで、水田からの浸透量を全く考慮してしない。浸透量は水田水位のみならず地下水位等にも影響されるため、偶然に水田浸透量が適正浸透量になる可能性はありうるが、従来技術では適正浸透量となるよう水田浸透量を能動的に操作できず、積極的な高収量化を図ることができない。
【０００８】
浸透量を制御量とした適正浸透量管理を行うためには、水田の浸透量推定を行い、推定浸透量が適正浸透量になるよう水位管理（給排水バルブ制御）を行うことか必要不可欠となる。しかるに、現在の技術レベルでは、浸透量を精度よく推定することは容易ではない。例えば、浸透量推定のためによく用いられている方法として、水田において成り立つ水収支式：水田貯留変化量＝給水量−排水量＋雨量−蒸発散量−浸透量、に対して、水田貯留変化量、給水量、排水量、雨量を測定し、蒸発散量を何らかの方法により推定することにより、残された浸透量を算出する方法がある。
【０００９】
しかしながら実際は、上記蒸発散量の推定が容易でないため、上記水収支式に基づく方法は有効とは言えない。傾度法、熱収支法等を用いれば蒸発散量推定可能であるが、特殊な測定器を必要とするので一般水田の水管理に用いることは現実的ではない。夜間は蒸発散量が０になることを利用すれば、上記方法でも浸透量の推定は可能となるが、昼間の浸透量推定には利用できず一般的ではない。
【００１０】
また、仮に水田の浸透量を推定できたとしても、推定浸透量が適正浸透量になるように水管理を行うことは容易ではない。浸透量は水田水位のみならず、地下水位、土壌・土層の状態等にも大きく左右されるため、適正浸透量を実現するための水田水位を決定することが困難なためである。
【００１１】
本発明の目的は、水稲収量に大きな影響を及ぼす物理量として浸透量に着目し、水田群の水位トレンド、給排水・雨量情報のみから、水田の土壌・土層、地下水位等の状態に依存することなく水田の浸透量を推定し、各水田の浸透量が稲の生育に応じた適正浸透量となるよう水田水位を適切に管理することができる、水田水管理の技術を提供することにある。
【００１２】
【課題を解決するための手段】
上記目的達成のため、次の手段を設ける。
（手段１）
複数区画からなる水田に対する水管理を行う水田群の水管理システムにおいて、上記水田の存在する水田地域の雨量を把握する手段と、上記各水田の水位を計測する水位計と、各水田に給水および排水を行う給排水装置と、各水田の浸透量の目標値およびその許容範囲を決定する目標浸透量決定手段と、上記把握された雨量および各水田の水位に基づいて、所定アルゴリズムに従って各水田の浸透量を推定する浸透量推定手段と、上記推定された浸透量の目標値の許容範囲に収まるように前記給排水装置を動作させる制御手段とを有することを特徴とする。
【００１３】
（手段２）
複数区画からなる水田に対する水管理を行う水田群の水管理システムにおいて、上記水田の存在する水田地域の雨量を把握する手段と、上記各水田の水位を計測する水位計と、各水田に給水および排水を行う給排水装置と、各水田の浸透量の目標値およびその許容範囲を決定する目標浸透量決定手段と、上記把握された雨量および各水田の水位に基づいて、所定アルゴリズムに従って各水田の浸透量を推定する浸透量推定手段と、上記推定された浸透量の目標値の許容範囲に収まるように各水田の目標水位及びその許容範囲を決定する目標水位決定手段とを有することを特徴とする。
【００１４】
（手段３）
本発明の水田の水管理システムでは、手段１、２における所定のアルゴリズムは、各水田に対して成り立つ水収支式において、蒸発散量、及び地下水位、透水係数等の浸透に影響を及ぼす物理量が上記各式とも共通の値を取ると仮定して、上記水収支式を連立して解くことにより各水田の浸透量を算出すること、を特徴とする。
【００１５】
（手段４）
本発明の水田の水管理システムでは、手段１、２における所定の目標量及びその許容範囲は、稲の生育段階又は日付と、該当希稲種の目標浸透量及びその許容範囲との対応関係を記録したテーブルに基づいて決定すること、を特徴とする。
【００１６】
（手段５）
本発明の水田の水管理システムでは、水田地域の気温、水温、湿度、日射豊、日照時間等の気象情報を把握するための手段を有し、手段１、２における所定の目標量及びその許容範囲は、稲の生育段階又は日付と、該当稲種の目標浸透量及びその許容範囲との対応関係を記録したテーブルに基づいて決定し、上記目標浸透量又はその許容範囲は、上記把握された気象情報に基づいて所定量だけ修正すること、を特徴とする。
【００１７】
（手段６）
本発明の水田の水管理システムでは、手段２における目標水位決定手段は、各水田に対して成り立つ水収支式において、蒸発散量、及び地下水位、透水係数等の浸透に影響を及ぼす物理量が上記各式とも共通の値を取ると仮定し、上記水収支式を連立して解くことにより上記物理量の値を算出し、上記算出された物理量の値及び前記決定された浸透量の目標値及び浸透量の上下限値を上記各水田ごとの水収支式に代入して解くことにより、前記各水田の目標水位及びその許容範囲を決定すること、を特徴とする。
【００１８】
（手段７）
本発明の水田の水管理システムでは、複数区画からなる水田地域において、上記水田地域の雨量を把握するための手段と、上記各水田の水位を計測するための水位計と、上記把握された雨量及び各水田の水位に基づいて、所定のアルゴリズムに従って各水田の水収支に関連する物理量を推定する物理量推定手段と、を有することを特徴とする。
【００１９】
（手段８）
本発明の水田の水管理システムでは、手段７における水収支に関連する物理量は、浸透量、地下水位、透水係数、その他浸透に影響を及ぼす物理量、蒸発散量、給水量、排水量のいずれかであること、を特徴とする。
【００２０】
（手段９）
本発明の水田の水管理システムでは、手段７における所定のアルゴリズムは、各水田に対して成り立つ水収支式において、蒸発散量、及び地下水位、透水係数等の浸透に影響を及ぼす物理量が上記各式とも共通の値を取ると仮定して、上記水収支式を連立立して解くことにより各水田の任意の期間中の水の移動に関連する物理量を算出すること、を特徴とする。
【００２１】
（手段１０）
本発明の水田の水管理システムでは、複数区画からなる水田地域において、上記水田地域の雨量を把握するための手段と、上記各水田の水位を計測するための水位計と、各水田に給水又は排水を行うための給排水装置と、各水田の水位変化量の目標値及びその許容範囲を決定する目標水位変化量決定手段と、上記把握された雨量及び各水田の水位に基づいて、所定のアルゴリズムに従って降雨の影響又は給排水の影響又は蒸発散の影響を取り除いた各水田の水位変化量を推定する水位変化量推定手段と、上記推定された水位変化量が上記水位変化量の目標値の許容範囲内に収まるように各水田の給排水装置を動作させる制御手段と、を有することを特徴とする。
【００２２】
（手段１１）
本発明の水田の水管理システムでは、複数区画からなる水田地域において、上記水田地域の雨量を把握するための手段と、上記客水田の水位を計測するための水位計と、各水田に給水又は排水を行うための給排水装置と、各水田の水位変化量の目標位及びその許容範囲を決定する目標水位変化量決定手段と、上記把握された雨量及び各水田の水位に基づいて、所定のアルゴリズムに従って降雨の影響又は給排水の影響又は蒸発散の影響を取り除いた各水田の水位変化量を推定する水位変化量推定手段と、上記推定された水位変化量が上記水位変化量の目標値の許容範囲内に収まるように各水田の目標水位及びその許容範囲を決定する目標水位決定手段と、を有することを特徴とする。
【００２３】
（手段１２）
本発明の複数区画からなる水田地域における水田の水管理を行うためのコンピュータ読み取り可能な媒体では、上記水画地域の雨量を把握する機能と、上記客水田に設置された水位計からの水位計測値を把握する機能と、各水田の浸透量の目標値及びその許容範囲を決定する目標浸透量決定機能と、上記把握された雨量及び各水田の水位に基づいて、所定のアルゴリズムに従って各水田の浸透量を推定する浸透量推定技能と、上記推定された浸透量が上記浸透量の目標値の許容範囲内に収まるように各水田に設置された給排水装置を動作させる制御機能と、を実行させるためのプログラムを記録したことを特徴とする。
【００２４】
（手段１３）
本発明の複数区画からなる水田地域における水田の水管理を行うためのコンピュータ読み取り可能な媒体では、上記水田地域の雨量を把握する機能と、上記客水田に設置された水位計からの水位計測値を把握する機能と、各水田の浸透量の目標値及びその許容範囲を決定する目標浸透量決定機能と、上記把握された雨量及び各水田の水位に基づいて、所定のアルゴリズムに従って各水田の浸透量を推定する浸透量推定機能と、上記推定された浸透量が上記浸透量の目標値の許容範囲内に収まるように各水田の目標水位及びその許容範囲を決定する目標水位決定機能と、を実行させるためのプログラムを記録したことを特徴とする。
【００２５】
（手段１４）
本発明の複数区画からなる水田地域における水田の水管理を行うためのコンピュータ読み取り可能な媒体では、手段１２、１３における所定のアルゴリズムは、各水田に対して成り立つ水収支式において、蒸発散量、及び地下水位、透水係数等の浸透に影響を及ぼす物理量が上記各式とも共通の値を取ると仮定して、上記水収支式を連立して解くことにより各水田の浸透量を算出すること、を特徴とする。
【００２６】
（手段１５）
本発明の複数区画からなる水田地域における水田の水管理を行うためのコンピュータ読み取り可能な媒体では、手段１２、１３における所定の目標量及びその許容範囲は、稲の生育段階又は日付と、該当稲種の目標浸透量及びその許容範囲との対応関係を記録したテーブルに基づいて決定すること、を特徴とするとする。
【００２７】
（手段１６）
本発明の複数区画からなる水田地域における水田の水管理を行うためのコンピュータ読み取り可能な媒体では、上記水田地域の気温、水温、湿度、日射量、日照時間等の気象情報を把握する機能を実行させるためのプログラムを記録し、手段１２、１３における所定の目標量及びその許容範囲は、稲の生育段階又は日付と、該当稲種の目標浸透量及びその許容範囲との対応関係を記録したテーブルに基づいて決定し、上記目標浸透量又はその許容範囲は、上記把握された気象情報に基づいて所定量だけ修正すること、を特徴とする。
【００２８】
（手段１７）
本発明の複数区画からなる水田地域における水田の水管理を行うためのコンピュータ読み取り可能な媒体では、手段１３における目標水位決定機能は、各水田に対して成り立つ水収支式において、蒸発散量、及び地下水位、透水係数等の浸透に影響饗を及ぼす物理量が上記各式とも共通の値を取ると仮定し、上記水収支式を連立して解くことにより上記物理量の値を算出し、上記算出された物理量の値及び前記決定された浸透量の目標値及び浸透量の上下限値を上記各水田ごとの水収支式に代入して解くことにより、前記各水田の目標水位及びその許容範囲を決定すること、を特徴とする。
【００２９】
（手段１８）
本発明の複数区画からなる水田地域における水田の水管理を行うためのコンピュータ読み取り可能な媒体では、上記水田地域の雨量を把握する機能と、上記各水田に設置された水位計からの水位計測値を把握する機能と、上記把握された雨量及び各水田の水位に基づいて、所定のアルゴリズムに従って各水田の水収支に関連する物理量を推定する物理量推定機能と、を実行させるためのプログラムを記録したことを特徴とする。
【００３０】
（手段１９）
本発明の複数区画からなる水田地域における水田の水管理を行うためのコンピュータ読み取り可能な媒体では、手段１８における水収支に関連する物理量は、浸透量、地下水位、透水係数、その他浸透に影響を及ぼす物理量、蒸発散量、給水量、排水量のいずれかであることを特徴とするコンピュータ読み取り可能な媒体。
【００３１】
（手段２０）
本発明の複数区画からなる水田地域における水田の水管理を行うためのコンピュータ読み取り可能な媒体では、手段１８における所定のアルゴリズムは、各水田に対して成り立つ水収支式において、蒸発散量、及び地下水位、透水係数等の浸透に影響を及ぼす物理量が上記各式とも共通の値を取ると仮定して、上記水収支式を連立して解くことにより各水田の任意の期間中の水の移動に関連する物理量を算出すること、を特徴とする。
【００３２】
（手段２１）
本発明の複数区画からなる水田地域における水田の水管理を行うためのコンピュータ読み取り可能な媒体では、上記水田地域の雨量を把握する機能と、上記各水田に設置された水位計からの水位計測値を把握する機能と、各水田の水位変化量の目標値及びその許容範囲を決定する目標水位変化量決定機能と、上記把握された雨量及び各水田の水位に基づいて、所定のアルゴリズムに従って降雨の影響又は給排水の影響又は蒸発散の影響を取り除いた各水田の水位変化量を推定する水位変化量推定機能と、上記推定された水位変化量が上記水位変化量の目標値の許容範囲内に収まるように各水田に設置された給排水装置を動作させる制御機能と、を実行させるためのプログラムを記録したことを特徴とする。
【００３３】
（手段２２）
本発明の複数区画からなる水田地域における水田の水管理を行うためのコンピュータ読み取り可能な媒体では、上記水田地域の雨量を把握する機能と、上記各水田に設置された水位計からの水位計測値を把握する機能と、各水田の水位変化量の目標値及びその許容範囲を決定する目標水位変化量決定機能と、上記把握された雨量及び各水田の水位に基づいて、所定のアルゴリズムに従って降雨の影響又は給排水の影響又は蒸発散の影響を取り除いた各水田の水位変化量を推定する水位変化量推定機能と、上記推定された水位変化量が上記水位変化量の目標値の許容範囲内に収まるように各水田の目標水位及びその許容範囲を決定する目標水位決定機能と、を実行させるためのプログラムを記録したことを特徴とする。
【００３４】
【発明の実施の形態】
以下、本発明を図面に示す実施形態によって詳細に説明する。
図１は、水田群における各水田の水位を管理する水管理システムの、実施形態を示すシステム構成図である。本実施形態の水管理システムは、水位計１−ｊ（ｊ＝１、…、３０）、給水バルブ２−ｊ、排水バルブ３−ｊ、コンピュータ１０４、データ伝送路１０５、気象情報提供システム１０６、気象ロボット１０７、給水パイプライン１０８、排水パイプライン１０９から構成される。
【００３５】
水管理の対象となる水田地域は、農水省が推進している大区画水田モデル地区と同様に、３０枚の１ｈａ規模の水田４−ｊによって構成される。各水田４−ｊ内には、水位計１−ｊ、給水バルブ２−ｊ、排水バルブ３−ｊが１つづつ設置されており、所定計測周期ごとに各水田４−ｊの水位、給水量、排水量をデータ伝送路１０５を介してコンピュータ１０４に送信している。ここで、データ伝送路１０５は、無線又は有線のどちらでも構わない。上記水田地域には、気象ロボット１０７が設置されており、上記水位計測周期と同じタイミングで、上記水田地域における気温、湿度、日照時間、日射量、風向・風速、雨量を、データ伝送路１０５を介してコンピュータ１０４に送信している。
【００３６】
更に、上記気象データ（気温、湿度、日照時間、日射量、風向・風速、雨量）の所定周期、所定時刻先までの予報値を、日本気象協会、又は民間団体の保有する気象情報提供システム１０６から、所定タイミングごとに、データ伝送路１０５を介してコンピュータ１０４に送信している。
【００３７】
コンピュータ１０４は、図２に示すように、ＣＰＵ２２1、ＲＡＭ２２２、ハードディスク装置２２３、フロッピーディスク装置２２４、ＭＯディスク装置２２５、ＤＶＤ装置２２６、ＣＤ−ＲＯＭ装置２２７、ＲＡＭカード装置２２８、入力装置２２９、表示装置２３０によって構成される。ＣＰＵ２２１は、水管理システム全体の動作を制御して、水稲収量に大きな影響を及ぼす物理量である浸透量を制御量として、給水バルブ２−ｊ、排水バルブ３−ｊの操作を行い、各水田４−ｊの水管理を行う中央処理装置である。ＲＡＭ２２２は、各種処理プログラムや各種データをロードする記憶装置である。
【００３８】
ハードディスク装置２２３は、各種処理プログラム２０１〜２０５、及び各種データベース２１１、２１２を磁気ディスクに格納する記憶装置である。フロッピーディスク装置２２４は、各種処理プログラム２０１〜２０５、及び各種データベース２１１、２１２が記録されたフロッピーディスクの読み書きを行う装置である。ＭＯディスク装置２２５は、各種処理プログラム２０１〜２０５、及び各種データベース２１１、２１２が記録されたＭＯディスク（光磁気ディスク）の読み書きを行う装置である。
【００３９】
ＤＶＤ装置２２６は、各種処理プログラム２０１〜２０５、及び各種データベース２１１、２１２が記録されたＤＶＤ（デジタルビデオディスク）の読み書きを行う装置である。ＣＤ−ＲＯＭ装置２２７は、各種処理プログラム２０１〜２０５、及び各種データベース２１１、２１２が記録されたＣＤ−ＲＯＭ（コンパクトディスク）の読み書きを行う装置である。ＲＡＭカード装置２２７は、各種処理プログラム２０１〜２０５、及び各種データベース２１１、２１２が記録されたＲＡＭカードの読み書きを行う装置である。
【００４０】
入力装置２２９は、水田浸透量推定プログラム２０３実行に必要となる各水田４−ｊの物理形状データや、稲の生育段階推定処理プログラム２０１実行に必要となる各水田４−ｊの稲の品種や田植え日等を入力するための装置である。表示装置２３０は、各種処理プログラム２０１〜２０５によって算出された各水田４−ｊの目標浸透量、水田浸透量、目標水位、バルブ開閉状態／及び稲の生育段階等、並びに気象情報提供システム１０６から送信された上記水田地域の気象データ（気温、湿度、日照時間、日射量、風向・風速、雨量）の予報値、並びに気象ロボット１０７、水位計１−ｊ、給水バルブ２−ｊ、排水バルブ３−ｊから送信された各水田４−ｊの水位、給水量、及び排水量等の表示を行う装置である。
【００４１】
ここで、上記各種処理プログラム２０１〜２０５、及び各種データベース２１１、２１２は上記記憶媒体のいずれに記録していてもかまわないため、以下、記録媒体の読み書きを行う装置がハードディスク装置２２３のみの場合を例にとって説明する。コンピュータ１０４では、各水田４−ｊの浸透量が稲の生育に応じた適正浸透量となるように、水田水位の管理を行う。上記水田水位管理は、ハードディスク装置２２３に記憶された各種処理プログラム２０１〜２０５、及び各種データベース２１１、２１２をＲＡＭ２２２にロード・記憶し、ＣＰＵ２２１によって、各種処理プログラム２０１〜２０５に対応する以下の処理：
（１）各水田４−ｊの現在の稲の生育段階を推定する処理
（２）各水田４−ｊの上記生育段階に応じた適切な目標浸透量を決定する処理
（３）各水田４−ｊの浸透量を推定する処理
（４）各水田４−ｊの上記目標浸透量に対応、する目標水田水位を決定する処理
（５）各水田４−ｊの水位を上記目標水田水位の所定の近傍内に収めるよう、給水バルブ２−ｊ、排水バルブ３−ｊの自動制御を行う処理
を順次実行していくことにより実現される。
【００４２】
以下、ＣＰＵ２２１において実行される上記処理（１）〜（５）の実行方法について説明する。はじめに、生育段階推定処理プログラム２０１によって実現される、各水田４−ｊの現在の稲の生育段階を推定する処理（１）について、図２、図３を用いて説明する。コンピュータ１０４では、気象ロボット１０７より計測された気温日照時間のデータを、所定周期ごとに受信している。処理プログラム２０１では、１日１回日没後、上記所定周期ごとの気温、日照時間データから当日の昼間時の平均気温、日照時間を算出する。例えば、昼間時の平均気温は、昼間時における上記所定周期ごとの気温の平均を取ることにより算出できる。また、昼間時の日照時間は、昼間時における上記所定周期ごとの日照時間を積算することにより算出できる。
【００４３】
生育モデルデータベース２１１では、稲の品種別の生育段階推定モデルを配慮している。上記生育段階推定モデルとしては、国内において広く用いられているＤＶＩモデルを採用する。ＤＶＩモデルによれば、田植え後ｉ日目の生育値ＤＶＩは、出芽日からその日までの生育速度ＤＶＲを積算して、
【００４４】
【数１】
と記述される。ここで、図３に示すように、生育速度ＤＶＲは気温が高いほど、また短日条件ほど大きくなることが知られている。
【００４６】
処理プログラム２０１では、入力装置２２９より入力された各水田４−ｊの稲の品種に基づいて、上記データベース２１１から対応する生育段階推定モデルを選択する。更に、入力装置２２９より入力された各水田４−ｊの稲の田植え日、及び処理プログラム２０１にて上記算出した田植え日から当日までの昼間時における平均気温、日照時間を（数１）に代入することにより、各水田４−ｊの当日までの稲の生育段階（ＤＶＩ値）を算出する。上述のように処理プログラム２０１において、自動的に稲の生育段階を算出してもよいが、もちろん農家において稲の生育状態を観察し、生育段階（ＤＶＩ値）を処理２０１に代わって入力装置２２９からコンピュータ１０４に入力してもよい。上述のようにして、ＣＰＵ２２１において、処理プログラム２０１により、各水田４−ｊの当日までの稲の生育段階の推定が行われる。
【００４７】
次に、目標浸透量決定処理理プログラム２０２によって実現される、現在の稲の生育段階に応じた適切な目標浸透量を決定する処理（２）について、図２、図４を用いて説明する。目標浸透量データベース２１２では、稲の生育段階（ＤＶＩ値）に対応した、稲の品種別・地方別の、田植え日から収穫までの、目標浸透量テーブル、及び当該目標浸透量データに対する許容範囲テーブルを記憶している。
【００４８】
図４に、目標浸透量テーブル、及び許容範囲テーブルの一例を示す。田植え直後では、土壌の養分濃度も高く、土壌の還元障害も少ないため、浸透量を低めに抑え、養分の流亡を防ぐ。有効分けつ期では地温上昇のため浸透を増やし、無効分けつ期では、分けつを抑制するため浸透を減らし、中干しを行う。
【００４９】
幼穂形成期、穂ばらみ期を通して豊熟期の前半までは、夏期の高温のために有機物の分解が盛んで土壌は著しい還元状態となり、各種の有害物質を生じるため、適度の浸透と土壌の酸化を目的とした間断灌漑を行う。収穫期では、出穂２５〜３０日後を目安に落水を行い、浸透を０にし、土壌を固め刈り取りに備える。
【００５０】
上記目標浸透量テーブル、及び許容範囲テーブルは、過去の稲の生育データや水位・浸透量データ、気象データ等に基づいて、稲の生育段階（ＤＶＩ値）から目標浸透量、及び許容範囲が特定できるよう「低温時には目標浸透量を所定量だけ増加させる」といったノウハウを組み込んでおく。また、上記テーブルは固定ではなく、農家のオリジナリティを組み込めるよう、いつでも変更可能なフレキシブルなテーブルとする。
【００５１】
処理プログラム２０２では、処理プログラム２０１の実行周期と同じタイミングで、以下の処理を行う。処理プログラム２０２では、入力装置２２９より入力された各水田４−ｊの稲の品種・地方に基づいて、上記データベース２１２から対応する目標浸透量テーブル、及び許容範囲テーブルを選択する。更に、上記処理プログラム２０１より算出された各水田４−ｊの稲の生育段階と、上記選択された各水田４−ｊの目標浸透量テーブル、及び許容範囲テーブルとを比較して、各水田４−ｊの適切な目標浸透量、及び許容範囲を決定する。決定された目標浸透量、及び許容範囲は、処理プログラム２０２の今回実行時から次回実行時まで有効とする。
【００５２】
また、コンピュータ１０４では、気象情報提供システム１０６より、所定周期、所定時刻先までの気象予報データを、所定タイミングごとに受信している。上記目標浸透量テーブルに、例えば上記の「低温時には目標浸透量を所定量だけ増加させる」といったノウハウが対応している場合は、処理プログラム２０２では、決定された目標浸透量の有効期間（当日の所定時刻〜次の日の該所定時刻）における気温予報データを考慮し、決定された目標浸透量の修正を行う。上述のようにして、ＣＰＵ２２１において、処理プログラム２０２により、各水田４−ｊの目標浸透量、及びその許容範囲の決定が行われる。
【００５３】
次に、水田浸透量推定処理プログラム２０３によって実現される、各水田４−ｊの浸透量を推定する処理（３）について、図２、図５を用いて説明する。コンピュータ１０４では、水位計１−ｊ、給水バルブ２−ｊ、排水バルブ３−ｊより計測された各水田４−ｊの水位、給水量、排水量、及び気象ロボット１０７より計測された当該水田地域の雨量を、所定周期ごとに受信している。処理プログラム２０３では、１日１回乃至数回（例えば３時間ごと）、上記受信データを用いて、各水田４−ｊの現在の浸透量の推定を行う。土壌中の浸透量は、Darcyの法則「浸透速度は動水勾配に正比例する」に支配され、水田４−ｊにおける浸透量（浸透速度）Ｖｊは、閉鎖浸透時（負圧の土層が存在しない浸透状態）、
【００５４】
【数２】
として記述できる。ここで、土壌条件の違いにより、浸透には閉鎖浸透と開放浸透の２つの状態が存在しうるが、開放浸透時でも（数２）と同様の式が成り立ち、開放浸透時に成り立つ式の方が未知変数も少なく単純であるため、より複雑な閉鎖浸透時に的を絞って説明を進める。
【００５６】
ここで、各水田４−ｊの水の移動に関連する物理量を図５に示す記号で表すことにする。通常、同一水田地域内の各水田においては、ほぼ同一の土壌・土層を有しており、またほほ同一の地下水位を有していることから、上記（数２）における合成透水係数Ｋ、全土層厚Ｌ、地下水位Ｈ(t）は各水田４−ｊにおいて共通であると仮定できる。また同一水田地域内の各水田においては、ほほ同一の気候、気温、湿度、日射、風速等）条件下にあることから、各水田４−ｊの（単位面積当りの）蒸発散量Ｅは共通であると仮定できる。上述のことから、現在までの所定期間△Ｔ（前回浸透量推定実行時〜今回浸透量推定実行時）中の各水田４−ｊにおける水収支式は、
【００５７】
【数３】
として記述できる。
【００５９】
連立方程式（数３）において、未知変数は、Ｋ、Ｌ、Ｈ（t）、Ｈ（ｔ−△Ｔ）、Ｅ（t）の５個であり、その他の変数は、コンピュータ１０４が計測している物理量から一意かつ容易に算出できる。従って、未知変数の個数と同じ５枚以上の水田があれば連立方程式（数３）は求解可能となり、本実施形態では水管理の対象水田を３０枚としているため十分求解可能となる。
【００６０】
処理プログラム２０３では、所定周期ごとに計測している各種物理量から、連立方程式（数３）を構築し、適当な求解手法に基づいて連立方程式（数３）の求解を行い、各水田４−ｊの浸透量を推定する。但し、浸透量推定期間△Ｔ中の水田水位、及び給排水量が完全に一致する水田が複数枚存在する場合は、該水田に対応する（数３）の方程式が一致してしまい、連立方程式（数３）が不定になり求解できない。
【００６１】
この場合、処理２０３プログラムでは、上記水田に対して、後述する目標水田水位の許容範囲内に水田水位が収まるよう給水又は排水を行い、連立方程式（数３）が不定になるのを防止する。上述のようにして、ＣＰＵ２２１において、処理プログラム２０３により、各水田４−ｊの浸透量の推定が行われる。
【００６２】
ところで、処理プログラム２０３では、各水田４−ｊの給水量及び排水量を計測するものと仮定しているが、以下の方法を用いることにより上記給排水量は計測不要となり、更に上記給排水量は推定可能とさえなる。即ち、連立方程式（数３）を構成する水収支式を、給排水を行っていない水田に対応する水収支式（給排水量の項＝０）に限定して該連立方程式（数３）の求解を行い、各水田４−ｊの浸透量、合成透水係数Ｋ、全土層厚Ｌ、地下水位（t）、蒸発散量Ｅを推定する。上記推定値を、給排水を行った水田に対応する水収支式（数３）に代入することにより、給排水量の値が推定可能となる。なお、給排水を行っていない水田の存在は、各水田４−ｊへの給排水の順序を適当に割り振ることにより容易に実現可能となる。
【００６３】
ここで、本発明の大きな特徴である、上述の水田群の浸透量推定方法の基本的考え方について説明しておく。従来、１水田における水収支式に基づいて浸透量を推定することは、方程式の個数より未知数（蒸発散量、地下水位、透水係数）の個数の方が多いため、求解不可能であった。しかしながら本発明では、対象水田を、同一地域内に存在する複数枚の水田までに拡大することにより、その求解を可能にしている。即ち、対象水田を複数枚に増加させることで上記方程式の個数を増加させ、更に、対象水田か何一地域内に存在することから、各水田は同一土壌条件下、同一気象条件下にあると仮定し、上記未知変数（蒸発散量、地下水位、透水係数）は共通の値をとるとして、対象水田の増加に伴う未知変数の増加を防止している。これにより、方程式の個数が未知数の個数を上回るようになり、その求解を可能としている。
【００６４】
次に、目標水位決定処理プログラム２０４によって実現される、各水田４−ｊの目標浸透量に対応する目標水田水位を決定する処理（４）について、図２、図５を用いて説明する。処理プログラム２０４では、処理プログラム２０３の実行周期（浸透量推定周期）と同じタイミングで、処理プログラム２０２にて算出した各水田４−ｊの目標浸透量、及び許容範囲に対応する、各水田４−ｊの目標水位、及び許容範囲を決定する。処理プログラム２０３の実行結果から、図５に示す当該水田地域の土層の合成透水係数Ｋ、全土層厚Ｌ、地下水位Ｈ（t）は全て算出されているため、上記値を（数２）に代入すれば、（数２）は水田水位と水田浸透量との１対１の関係を示す式となる。従って、上記各水田４−ｊにおける３つの浸透量（目標浸透量、上限浸透量、下限浸透量）を上記関係式に代入することにより、各水田４−ｊの目標水位、及び許容範囲を決定する。決定された目標水位、及び許容範囲は、処理２０４の今回実行時から次回実行時まで有効とする。上述のようにして、ＣＰＵ２２１において、処理プログラム２０４により、各水田４−ｊの目標未位、及びその許容範囲の決定が行われる。
【００６５】
次に、給排水バルブ制御処理プログラム２０５によって実現される、各水田４−ｊの水位を上記目標水田水位の所定の近傍内に収めるよう、給水バルブ２−ｊ、排水バルブ３−ｊの自動制御を行う処理（５）について、図２、図６を用いて説明する。コンピュータ１０４では、水位計１−ｊより計測された各水田４−ｊの水位を、所定周期ごとに受信している。処理プログラム２０５では、上記受信データを用いて、上記水位計測周期と同じタイミングで、各給水バルブ２−ｊ、排水バルブ３−ｊの開閉状態を決定し、各バルブに対して開閉信号の送信を行う。
【００６６】
図６は、処理プログラム２０５が実行するバルブ制御処理の概要を示すフローチャートである。はじめに、各水位計１−ｊより計測した最新の水田水位を読み込み（ステップ６０１）、また各水田４−ｊの目標水位、及びその許容範囲を、処理プログラム２０４において決定された最新値に更新する（ステップ６０２）。
【００６７】
次に処理プログラム２０５では、各水田４−ｊの水位が目標水位の所定の近傍Ｉ内（目標水位±許容範囲）に存在するように、以下に示す給排水バルブの開放／閉鎖条件に基づいて、給水バルブ２−ｊ、排水バルブ３−ｊの操作を行う。
【００６８】
・給水バルブ開放条件（閉鎖から開放に移行する条件）：水田水位＜目標水位−許容範囲
・給水バルブ閉鎖条件（開放から閉鎖に移行する条件）：水田水位＞目標水位
・排水バルブ開放条件：水田水位＞目標水位＋許容範囲
・排水バルブ閉鎖条件：水田水位＜目標水位
即ち、水田水位が近傍Ｉの下限値（目標水位−許容範囲）を下回っていればデータ伝送路１０５を介して、給水バルブ２−ｊに開放信号を送信し、給水を実行させる（ステップ６０３、６０４）。給水バルブ２−ｊでは、コンピュータ１０４からのバルブ開放信号を受信し、バルブを開けて給水パイプライン１０８から水田４−ｊ内に水を供給する。これにより水田４−ｊの水位が漸次上昇し、水田水位が目標水位を上回るようになればデータ伝送路１０５を介して、給水バルブ２−ｊに閉鎖信号を送信し、給水を停止させる（ステップ６０５、６０６）。これにより水田４−ｊの水位は目標水位付近で停止することになる。
【００６９】
逆に、水田水位が近傍Ｉの上限値（目標水位＋許容範囲）を上回っていればデータ伝送路１０５を介して、排水バルブ３−ｊに開放信号を送信し、排水を実行させる（ステップ６０７、６０８）。排水バルブ３−ｊでは、コンピュータ１０４からのバルブ開放信号を受信し、バルブを開けて水田４−ｊから排水パイプライン１０９内に水を排水する。これにより水田４−ｊの水位が漸次下降し、水田水位が目標水位を下回るようになれば、データ伝送路１０５を介して、排水バルブ３−ｊに閉鎖信号を送信し、排水を停止させる（ステップ６０９、６１０）。これにより水田４−ｊの水位は目標水位付近で停止することになる。上述のようにして、ＣＰＵ２２１において、処理プログラム２０５により、水田水位が目標水位の所定の近傍Ｉ内に存在するように、給水バルブ２−ｊ、及び排水バルブ３−ｊの自動制御が行われる。
【００７０】
上述のように、本発明の実施形態によれば、水稲収量に大きな影響を及ぼす物理量として浸透量に着目し、同一地域内の水田群では、各水田とも同一土壌下、同一気候条件下にあるという物理現象を利用して、水田群の水位トレンド、給排水、雨量情報から、水田の土壌・土層、地下水位等の状態に依存することなく水田の浸透量を推定し、各水田の浸透量が稲の生育に応じた適正浸透量なるよう水田水位を適切に管理することが可能である。
【００７１】
以上説明したように本発明の実施の形態によれば各水田において成り立つ水収支式に基づいて、地下水位、透水係数、蒸発散量等を計測不要な未知変数として浸透量推定計算を行っているため、水田の土壌・土質、地下水位、蒸発散量を把握するための特殊な計器を用いることなく、また昼夜の別に限定されることなく、各水田の現在の浸透量を推定できる。
【００７２】
また、稲の生育段階推定結果と、稲の生育段階に対応した目標浸透量テーブルに基づいて、各水田の適正浸透量の目標値を決定し、上記現推定浸透量との比較を行っているため、各水田において現在の稲の生育に応じた適正浸透が行われているか否かを判断できる。
【００７３】
適正浸透が行われていない水田に対しては、上記浸透量推定計算において、浸透に影響を及ぼす物理量（地下水位、透水係数）をも同時に把握しているため、適正浸透量を実現するのに過不足ない目標水田水位を算出でき、更に水田水位が上記目標水田水位に等しくなるよう給排水バルブの操作を行っているため、各水田の浸透量が適正浸透量になるような水管理を行うことができる。
【００７４】
且つ、上記浸透量推定周期、即ち水管理周期を任意の周期（例えば３時間程度の短期間）に設定できるため、各水田の浸透量を稲の生育期間中常に適正浸透量に維持することができ、より緻密な水管理を行うことができる。
【００７５】
即ち、上述の効果を総合すれば、本発明によれば、稲の植付から刈り取りまでの全生育期間に対して、各水田の浸透量を、現在の稲の生育に応じた適正浸透量に絶え間なく維持することができるため、水稲収量の最大化を図ることができる。
【００７６】
また、水田の土壌・土層の状態に応じて水田からの浸透特性は大きく違ってくるため、全国には漏水田や排水不良田など様々な水田が存在するが、本発明では、上記浸透量推定計算において、土壌・土層の状態を表すパラメータである透水係数をオンライン推定し、該透水係数に基づいて、水田の適正浸透量管理を行っているため、任意の水田地域に対して適用可能な、汎用性の高い水管理システムを構築できる。
【００７７】
【発明の効果】
以上の構成により、本発明においては、特殊な計器を用いなくとも、各水田の水の浸透量を把握することが可能になる。
【図面の簡単な説明】
【図１】本発明の水管理システムの実施形態を示すシステム構成図である。
【図２】コンピュータ内部のハード構成図及び各種処理プログラムのブロック図である。
【図３】平均気温、日照時間と稲の生育速度との関係を示す説明図である。
【図４】稲の生育段階に応じた目標水深テーブルの一例を示す説明図である。
【図５】１つの水田における水収支の関係を示す説明図である。
【図６】給排水バルブ制御の概要を示すフローチャートである。
【符号の説明】
１−ｊ…水位計、２−ｊ…給水バルブ、３−ｊ…排水バルブ、４−ｊ…水田、１０４…コンピュータ、１０５…データ伝送路、１０６…気象情報提供システム、１０７…気象ロボット、１０８…給水パイプライン、１０９…排水パイプライン。
Claims (6)
- 複数の水田を有する水田地域において，前記各水田の浸透量を制御する水田群の水管理システムであって，
前記各水田に，水位を測定する水位計と，給水装置と，排水装置と，を備え，
当該水田地域について測定された雨量を取得する手段と，
予め求めた各水田の目標浸透量を満たす，目標水位を求める推定手段と，
前記水位計により測定した水位が前記目標水位となるように，前記各水田の給水装置または排水装置を動作させる制御手段と，を備え，
前記推定手段は，浸透量の推定式
ことを特徴とする水田群の水管理システム。 - 請求項1に記載の水田群の水管理システムであって，
前記給水装置および排水装置の各々は，前記各水田への給水量および排水量を測定する手段を備え，
前記推定手段は，前記測定した給水量および排水量に基づき，前記目標水位を求める
ことを特徴とする水田群の水管理システム。 - 請求項1または2に記載の水田群の水管理システムであって，
さらに，
前記各水田における稲の生育段階または日付と適正浸透量の目標値との対応関係を記憶する手段と，
前記記憶された対応関係と，前記各水田の稲の生育段階または日付と，から，各水田の前記目標浸透量を決定する手段と，を備える
ことを特徴とする水田群の水管理システム。 - 請求項1ないし3いずれか一に記載の水田群の水管理システムであって，
前記制御手段は，前記推定手段により前記目標水位が求まらなかった場合に，前記各水田の水位が異なるように前記給水装置または排水装置を動作させ，
前記推定手段は，前記各水田の水位が異なっているときの水位測定値hj(t)に基づき，前記目標水位を求める
ことを特徴とする水田群の水管理システム。 - 各々が，水位を測定する水位計と，給水装置と，排水装置と，を備える複数の水田を有する水田地域において，コンピュータによる前記各水田の浸透量を制御する水田群の水管理方法であって，
当該水田地域について測定された雨量を取得する工程と，
予め求めた各水田の目標浸透量を満たす，目標水位を求める推定工程と，
前記水位計により測定した水位が前記目標水位となるように，前記各水田の給水装置または排水装置を動作させる制御工程と，を備え，
前記推定工程は，浸透量の推定式
ことを特徴とする水田群の水管理方法。 - 各々が，水位を測定する水位計と，給水装置と，排水装置と，を備える複数の水田を有する水田地域において，コンピュータに前記各水田の浸透量を制御させる，水田群の水管理プログラムを記録したコンピュータ読み取り可能な媒体であって，
前記プログラムは，前記コンピュータに，
当該水田地域について測定された雨量を取得する工程と，
予め求めた各水田の目標浸透量を満たす，目標水位を求める推定工程と，
前記水位計により測定した水位が前記目標水位となるように，前記各水田の給水装置または排水装置を動作させる制御工程と，を行わせるものであり，
前記推定工程は，浸透量の推定式
ことを特徴とする水田群の水管理プログラムを記録したコンピュータ読み取り可能な媒体。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP15416698A JP4006098B2 (ja) | 1998-06-03 | 1998-06-03 | 水田群の水管理システム，水田群の水管理方法，および，水田群の水管理を行うためのプログラムを記録したコンピュータ読み取り可能な媒体 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP15416698A JP4006098B2 (ja) | 1998-06-03 | 1998-06-03 | 水田群の水管理システム，水田群の水管理方法，および，水田群の水管理を行うためのプログラムを記録したコンピュータ読み取り可能な媒体 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPH11346581A JPH11346581A (ja) | 1999-12-21 |
| JP4006098B2 true JP4006098B2 (ja) | 2007-11-14 |
Family
ID=15578283
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP15416698A Expired - Fee Related JP4006098B2 (ja) | 1998-06-03 | 1998-06-03 | 水田群の水管理システム，水田群の水管理方法，および，水田群の水管理を行うためのプログラムを記録したコンピュータ読み取り可能な媒体 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP4006098B2 (ja) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2021087388A (ja) * | 2019-12-04 | 2021-06-10 | Kddi株式会社 | 水門制御装置、水門、水門制御方法及び水門制御システム |
| WO2021157033A1 (ja) * | 2020-02-06 | 2021-08-12 | 株式会社ナイルワークス | 水管理システム、水管理サーバ及び水管理方法 |
Families Citing this family (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| KR100449442B1 (ko) * | 2002-06-11 | 2004-10-01 | 주식회사 한성전자산업개발 | 농업용수 관리시스템 |
| JP4320463B2 (ja) * | 2003-03-19 | 2009-08-26 | ヤンマー農機株式会社 | 水田管理システム |
| KR100813881B1 (ko) | 2006-08-28 | 2008-03-18 | 전자부품연구원 | 센서 네트워크를 이용한 농경지 수위 조절 방법 및 시스템 |
| JP6980246B2 (ja) * | 2017-03-24 | 2021-12-15 | 独立行政法人国立高等専門学校機構 | プログラム及び潅漑制御装置 |
| CN110604045A (zh) * | 2019-10-08 | 2019-12-24 | 北京享云智汇科技有限公司 | 一种农田灌溉排水控制系统 |
| JP2021078461A (ja) * | 2019-11-21 | 2021-05-27 | 株式会社ナイルワークス | 植物の生育管理システム、植物の生育管理方法、及び植物の生育管理プログラム |
| CN114698431B (zh) * | 2022-05-07 | 2023-05-05 | 浙江农林大学 | 一种沼液定期定量还田智能化控制系统 |
| JP7425846B1 (ja) * | 2022-09-30 | 2024-01-31 | 株式会社インターネットイニシアティブ | 水位管理システムおよび水位管理方法 |
-
1998
- 1998-06-03 JP JP15416698A patent/JP4006098B2/ja not_active Expired - Fee Related
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2021087388A (ja) * | 2019-12-04 | 2021-06-10 | Kddi株式会社 | 水門制御装置、水門、水門制御方法及び水門制御システム |
| WO2021157033A1 (ja) * | 2020-02-06 | 2021-08-12 | 株式会社ナイルワークス | 水管理システム、水管理サーバ及び水管理方法 |
| JP7412021B2 (ja) | 2020-02-06 | 2024-01-12 | 株式会社ナイルワークス | 水管理システム、水管理サーバ及び水管理方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| JPH11346581A (ja) | 1999-12-21 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| George et al. | Development and testing of an irrigation scheduling model | |
| EP3179319B1 (en) | Method for irrigation planning and system for its implementation | |
| Wiyo et al. | Effect of tied-ridging on soil water status of a maize crop under Malawi conditions | |
| Panigrahi et al. | Field test of a soil water balance simulation model | |
| Smith et al. | Estimating irrigation duration for high performance furrow irrigation on cracking clay soils | |
| Huang et al. | Water–yield relationships and optimal water management for winter wheat in the Loess Plateau of China | |
| Palomo et al. | Determining water consumption in olive orchards using the water balance approach | |
| JP4006098B2 (ja) | 水田群の水管理システム，水田群の水管理方法，および，水田群の水管理を行うためのプログラムを記録したコンピュータ読み取り可能な媒体 | |
| Reddy | A simple method of estimating the soil water balance | |
| CN115455086A (zh) | 一种灌区农业灌溉用水量的估算方法及系统 | |
| De Silva et al. | Groundwater recharge estimation using improved soil moisture balance methodology for a tropical climate with distinct dry seasons | |
| Rose et al. | A simulation model of growth-field environment relationships for Townsville Stylo (Stylosanthes humilis HBK) pasture | |
| Reyes-Esteves et al. | Modeling approaches for determining appropriate depth of subsurface drip irrigation tubing in alfalfa | |
| Ojha et al. | Evaluation of a nonlinear root-water uptake model | |
| KR102461749B1 (ko) | 관수 결정 장치 및 방법 | |
| Pandey et al. | Lysimeter based crop coefficients for estimation of crop evapotranspiration of black gram (Vigna Mungo L.) in sub-humid region | |
| Coolong et al. | The relationship between soil water potential, environmental factors, and plant moisture status for poblano pepper grown using tensiometer-scheduled irrigation | |
| CN110432046B (zh) | 一种温室内的智能灌溉系统 | |
| Bautista-Capetillo et al. | Using thermal units for crop coefficient estimation and irrigation scheduling improves yield and water productivity of corn (Zea mays L.) | |
| KC | Irrigation scheduling: a soft adaptor to weather uncertainties and irrigation efficiency improvement initiatives | |
| Flores-Cayuela et al. | Verifiable water use inventory using ICTs in industrial agriculture | |
| WO2020047587A1 (en) | System and method for sensor-based auto-calibration of soil-moisture levels | |
| Ravelo et al. | Incorporating crop needs into drainage system design | |
| Amini et al. | Determination of water requirement and crop coefficient for strawberry using lysimeter experiment in a semi-arid climate | |
| Schepers et al. | Integrated water and nitrogen management |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20050324 |
|
| A711 | Notification of change in applicant |
Free format text: JAPANESE INTERMEDIATE CODE: A712 Effective date: 20061110 |
|
| A521 | Written amendment |
Free format text: JAPANESE INTERMEDIATE CODE: A821 Effective date: 20061110 |
|
| A521 | Written amendment |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20061227 |
|
| RD02 | Notification of acceptance of power of attorney |
Free format text: JAPANESE INTERMEDIATE CODE: A7422 Effective date: 20061227 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20070403 |
|
| A521 | Written amendment |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20070427 |
|
| TRDD | Decision of grant or rejection written | ||
| A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20070807 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20070827 |
|
| R150 | Certificate of patent or registration of utility model |
Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100831 Year of fee payment: 3 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100831 Year of fee payment: 3 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110831 Year of fee payment: 4 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120831 Year of fee payment: 5 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120831 Year of fee payment: 5 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130831 Year of fee payment: 6 |
|
| S533 | Written request for registration of change of name |
Free format text: JAPANESE INTERMEDIATE CODE: R313533 |
|
| S533 | Written request for registration of change of name |
Free format text: JAPANESE INTERMEDIATE CODE: R313533 |
|
| R350 | Written notification of registration of transfer |
Free format text: JAPANESE INTERMEDIATE CODE: R350 |
|
| LAPS | Cancellation because of no payment of annual fees |