CN110419415B - Rainfall forecast-based large irrigation area paddy field irrigation plan optimization method - Google Patents

Rainfall forecast-based large irrigation area paddy field irrigation plan optimization method Download PDF

Info

Publication number
CN110419415B
CN110419415B CN201910352946.7A CN201910352946A CN110419415B CN 110419415 B CN110419415 B CN 110419415B CN 201910352946 A CN201910352946 A CN 201910352946A CN 110419415 B CN110419415 B CN 110419415B
Authority
CN
China
Prior art keywords
irrigation
day
area
water layer
group
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.)
Active
Application number
CN201910352946.7A
Other languages
Chinese (zh)
Other versions
CN110419415A (en
Inventor
程吉林
黄慧雯
汪靓
陈兴
蒋晓红
袁承斌
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Yangzhou University
Original Assignee
Yangzhou University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Yangzhou University filed Critical Yangzhou University
Priority to CN201910352946.7A priority Critical patent/CN110419415B/en
Publication of CN110419415A publication Critical patent/CN110419415A/en
Application granted granted Critical
Publication of CN110419415B publication Critical patent/CN110419415B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • A—HUMAN NECESSITIES
    • A01—AGRICULTURE; FORESTRY; ANIMAL HUSBANDRY; HUNTING; TRAPPING; FISHING
    • A01G—HORTICULTURE; CULTIVATION OF VEGETABLES, FLOWERS, RICE, FRUIT, VINES, HOPS OR SEAWEED; FORESTRY; WATERING
    • A01G25/00—Watering gardens, fields, sports grounds or the like
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00—Administration; Management
    • G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/02—Agriculture; Fishing; Forestry; Mining

Landscapes

  • Engineering & Computer Science (AREA)
  • Business, Economics & Management (AREA)
  • Strategic Management (AREA)
  • Human Resources & Organizations (AREA)
  • Economics (AREA)
  • Theoretical Computer Science (AREA)
  • General Business, Economics & Management (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Physics & Mathematics (AREA)
  • Tourism & Hospitality (AREA)
  • Marketing (AREA)
  • Health & Medical Sciences (AREA)
  • Environmental Sciences (AREA)
  • Animal Husbandry (AREA)
  • General Health & Medical Sciences (AREA)
  • Mining & Mineral Resources (AREA)
  • Marine Sciences & Fisheries (AREA)
  • Agronomy & Crop Science (AREA)
  • Water Supply & Treatment (AREA)
  • Primary Health Care (AREA)
  • Development Economics (AREA)
  • Game Theory and Decision Science (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Operations Research (AREA)
  • Quality & Reliability (AREA)
  • Cultivation Of Plants (AREA)

Abstract

The invention relates to a large irrigation area paddy field irrigation plan optimization method based on rainfall forecast, which comprises the following steps: (1) determining parameters of the required irrigation requirements; (2) measuring the real-time depth of the water layer of the rice field in the irrigation area one day before the initial day of irrigation, and determining the initial depth of the irrigation water layer; (3) dividing the irrigation area into irrigation rotation groups according to the distribution of main canal gates, and determining the time d required by each group for irrigationiAnd the time required to complete all irrigation ds(ii) a (4) Obtaining d required for optimal irrigationsForecasting data of daily average precipitation in the intra-day irrigation area; (5) calculating the irrigation required d through the average precipitation forecast data of the irrigation area, the initial field water layer depth of irrigation and the field leakagesThe depth of the water layer in the field day by day within the day; (6) and (3) completing the irrigation plan of each rotation irrigation group according to the field water layer depth and the irrigation requirements day by day: (7) and (4) repeating the steps (1) to (6) until the whole growth period of the rice is finished. By the method and the device, irrigation water can be optimized and reduced, and the utilization efficiency of water resources is improved.

Description

Rainfall forecast-based large irrigation area paddy field irrigation plan optimization method
Technical Field
The invention relates to a large irrigation area paddy field irrigation plan optimization method based on rainfall forecast, and belongs to the field of agricultural water and soil engineering.
Background
With the rapid development of economy and society in China and the continuous growth of population, the pressure of rice production and agricultural water resource utilization in China is increasing day by day. At present, irrigation water in the rice planting process is a main component of the whole agricultural water, and the optimization of the irrigation water of the rice has very important significance. At present, rice irrigation plans are made more extensively, although specific irrigation time of each rotation irrigation group in a large irrigation area is sequential, irrigation water consumption is fixed and unified, and natural rainfall resources are not fully considered and utilized. Such irrigation programs result in a large amount of water waste in the irrigation of large irrigation areas.
Disclosure of Invention
In order to overcome the defects of the prior art, an irrigation plan of a large-scale irrigation area rice field is scientifically and reasonably made, natural rainfall is fully utilized for irrigation, and waste of water resources is reduced; the invention provides a large irrigation area paddy field irrigation plan optimization method based on rainfall forecast. The invention is easy to be popularized in the practice of rice field irrigation in various large irrigation areas.
The invention aims to realize the purpose, and the method for optimizing the irrigation plan of the rice field in the large irrigation area based on rainfall forecast comprises the following steps:
(1) determining parameters of the required irrigation requirements according to the growth period, the irrigation mode, the irrigation time interval and the irrigation area location of the rice;
(2) measuring the real-time depth of the water layer of the rice field in the irrigation area in the day before the initial day of irrigation, and determining the initial depth of irrigation;
(3) dividing the irrigation area into irrigation rotation groups according to the distribution of main canal gates, and determining the time d required by each group for irrigationiAnd the time required to complete all irrigation ds(d s7 or less);
(4) obtaining optimal irrigation dsForecasting data of daily average precipitation in the intra-day irrigation area;
(5) calculating irrigation d by using a water balance equation according to the average precipitation forecast data of the irrigation area, the initial irrigation field water layer depth, the field leakage and the crop water demandsThe water layer depth of the field is gradually increased day by day in the irrigation area in the day;
(6) and finishing the irrigation plan of each rotation irrigation group according to the daily field water layer depth and the irrigation requirements.
(7) And (4) repeating the steps (1) to (6) until the whole growth period of the rice is finished.
2. Preferably, the step (1) specifically includes the steps of:
a. determining irrigation area sites, irrigation time periods, irrigation modes and growth periods of rice in the area during irrigation, wherein the irrigation area sites need to optimize irrigation water quantity;
b. determining the upper limit h of the suitable water layer of the irrigation area according to the informationmaxSuitable lower limit of water layer hminLower limit of soil moisture content thetaminDry volume weight of soil
Figure BDA0002044538220000021
Field seepage DP, water surface evaporation Z and alpha values, and field water holding capacity theta0The depth of the active layer S of the root system of the rice;
the α values are given by the following table:
Figure BDA0002044538220000022
c. according to the growth period of the rice, converting the water surface evaporation Z into crop water demand ET by using an alpha value method, namely the following formula;
ET＝αZ
d. according to the lower limit theta of the water content of the soilminDry volume weight of soil
Figure BDA0002044538220000023
The lower limit h of the (underground) wet water layer is calculated by the following formulasmax；
Figure BDA0002044538220000024
Where ρ is the density of water, ρ is 1g/cm3。
Preferably, the step (2) specifically includes the steps of:
a. determining the initial day of irrigation;
b. actually measuring the depth H of the water layer in the irrigation area one day before the initial day of irrigation0；
4. Preferably, the step (3) specifically includes the steps of:
a. determining the area of an irrigation area;
b. dividing the irrigation area into s rotation irrigation groups according to the area of the irrigation area and the distribution condition of main canal gates, wherein the area of each group is not more than 20 ten thousand mu;
c. determining the days d required by each group according to the area of the rotation irrigation groupi(ii) a The area of the irrigation block is less than or equal to 5 ten thousand mu, diEqual to 1; the area of the irrigation area group is more than 5 ten thousand mu and less than or equal to 10 ten thousand mu, diEqual to 2; the area of the irrigation area group is more than 10 ten thousand mu and less than or equal to 20 ten thousand mu, diEqual to 3;
e. calculating the total irrigation time d of the irrigation areas：
Figure BDA0002044538220000025
dsIt is required to be 7 or less, and if it is more than 7, it is not suitable for the method of the present invention.
f. Determining the sequence and date of irrigation of each rotation irrigation group;
preferably, the step (4) specifically includes the steps of:
a. acquiring all m weather forecast sites d in a circular area with the center of an irrigation area as the center of a circle and the radius of 100 kilometerssForecast data y of precipitation amount of each day in the dayi,j i＝1,2,...,ds；j＝1,2,...,m；
b. Calculating the distance r from the center of the irrigation area to m weather forecast stationsmAnd a total distance rs：
Figure BDA0002044538220000031
c. Calculating the weight omega of the rainfall of each meteorological site in the average rainfall forecast data of the irrigation districtj,j＝1,2,...,m：
Figure BDA0002044538220000032
d. Calculating irrigation area dsAverage precipitation forecast data T of each day in the dayi,i＝1,2,...,ds：
Figure BDA0002044538220000033
Preferably, the step (5) specifically includes the steps of:
a. the depth H of the water layer in the initial field of irrigation is predicted according to the average precipitation of the irrigation area0Field leakage DP and crop water demand ET, calculating irrigation d by using the following formulasDaily field water layer depth Hi: when no surface water layer is present, the depth of the water layer refers to the distance from the surface of the ground water layer to the surface of the paddy field.
Hi＝Hi-1+Ti-ET-DP,i＝1,2,...,ds(ii) a Surface water layers are present on days i, i-1;
or
Hi＝-Hi-1+Ti-ET-DP,i＝1,2,...,ds(ii) a A surface water layer exists on the ith day, and no surface water layer exists on the ith-1 day;
or
Hi＝|Hi-1+Ti-ET-DP|,i＝1,2,...,ds(ii) a No surface water layer on day i, and surface water layer on day i-1;
or
Hi＝|-Hi-1+Ti-ET-DP|,i＝1,2,...,ds(ii) a No surface water layer exists on days i, i-1;
b. when not irrigating, the depth of the field water layer of each rotation irrigation group is equal.
Preferably, the step (6) specifically includes the steps of:
a. comparing the depth H of the water layer in the field of the ith day of the rotation irrigation group when the surface water layer existsiAnd hmaxThe size of (d); if H isiGreater than hmaxThe rotation irrigation group needs to drain water on the ith day with the water discharge hp＝Hi-hmax(ii) a If H isiH is less than or equal tomaxAnd rotation irrigation planWhen the irrigation of the group is turned to the irrigation, the group is irrigated to the water layer of the field to reach h on the ith daymaxThen, the group is excluded from the irrigation plan and is not irrigated any more; otherwise the group does not need irrigation nor drainage on day i;
b. comparing the water layer depth H of the ith day without surface water layeriAnd hsmaxThe size of (d); if H isiIs greater than or equal to hsmaxThe group must start irrigation on day i to the upper limit of the suitable water layer hmax(ii) a If H isiLess than hsmaxAnd the irrigation of the group is not turned to in the rotation irrigation plan, the group does not need irrigation on the ith day; if H isiLess than hsmaxBut the irrigation of the group is turned to in the rotation irrigation plan, the group is irrigated to reach the depth h of the field water layer on the ith daymaxThen, the group is excluded from the irrigation plan and is not irrigated any more;
e. repeating the steps a and b until all the rotation irrigation groups complete the irrigation plan.
The method is advanced and scientific, and discloses a method for optimizing the irrigation water quantity of a large-scale irrigation area paddy field based on rainfall forecast, which comprises the following steps: (1) determining parameters of the required irrigation requirements according to the growth period, the irrigation mode, the irrigation time interval and the irrigation area location of the rice; (2) measuring the real-time depth of the water layer of the rice field in the irrigation area one day before the initial day of irrigation, and determining the initial depth of the irrigation water layer; (3) dividing the irrigation area into irrigation rotation groups according to the distribution of main canal gates, and determining the time d required by each group for irrigationiAnd the time required to complete all irrigation ds(ds7 or less); (4) obtaining d required for optimal irrigationsForecasting data of daily average precipitation in the intra-day irrigation area; (5) calculating the irrigation required d through the average precipitation forecast data of the irrigation area, the initial field water layer depth of irrigation and the field leakagesThe depth of the water layer in the field day by day within the day; (6) and (3) completing the irrigation plan of each rotation irrigation group according to the field water layer depth and the irrigation requirements day by day: (7) and (4) repeating the steps (1) to (6) until the whole growth period of the rice is finished. The invention utilizes the full utilization of rainfall information to adjust the irrigation water quantity of the rice field, can optimize and reduce the irrigation water and improve the utilization efficiency of water resources.
Has the advantages that: the invention applies the high-efficiency modern rainfall forecasting technology, scientifically and reasonably utilizes water resources formed by natural rainfall and the arrangement of wheel irrigation in large irrigation areas, reduces the water demand in the rice planting process and saves the water for irrigating rice fields. The method is easy to popularize and apply in the rice irrigation process of various large irrigation areas.
Drawings
FIG. 1 is a flow chart of the present invention;
FIG. 2 is the division of the rotation irrigation group of the large irrigation area.
Detailed Description
The invention will be further explained with reference to the attached drawings and the practice of irrigation planning of certain large-scale paddy rice irrigation areas in the east of China:
(1) according to the steps of the flow chart shown in FIG. 1:
a. the irrigation plan is formulated for a large-scale gravity irrigation area of a certain plain in the east of China, the irrigation period is 7 days/7 to 13 days/7 months, the irrigation mode is shallow irrigation, and the tillering stage is in the rice growth period;
b. determining the upper limit h of the suitable water layer of the irrigation area according to the informationmax20mm, suitable lower limit of water layer hmin0.0mm, lower limit of soil moisture content thetamin0.45 dry volume weight of soil
Figure BDA0002044538220000051
Field seepage DP of 2.2mm/d, water surface evaporation Z of 5.96mm/d and alpha of 1.04, and field water holding capacity theta0The depth S of the root system moving layer of the rice is 200mm when being equal to 0.4;
c. according to the growth period of the rice, converting the water surface evaporation Z into crop water demand ET by using an alpha value method, namely the following formula;
ET＝αZ
ET is 6.2mm/d
d. According to the lower limit theta of the water content of the soilminDry volume weight of soil
Figure BDA0002044538220000052
The lower limit h of the (underground) wet water layer is calculated by the following formulasmax；
Figure BDA0002044538220000053
Where h issmax＝11.6mm
(2) According to the steps of the specification and the flow chart:
a. the initial day of the irrigation is 7 months and 7 days;
b. the depth of the water layer on the day before the irrigation is 20mm, namely the depth of the initial water layer for irrigation is H0＝20mm。
(3) According to the steps of the specification and the flow chart:
a. determining the area of an irrigation area; the area of the irrigation area is 34.0 ten thousand mu.
b. Dividing the irrigation area into s rotation irrigation groups according to the area of the irrigation area and the distribution condition of main canal gates, wherein the area of each group is not more than 20 ten thousand mu;
the irrigation divides the whole irrigation area into 3 irrigation groups, namely s-3, and is shown in fig. 2 in detail.
c. Determining the days d required by each group according to the area of the rotation irrigation groupi(ii) a The area of the irrigation block is less than or equal to 5 ten thousand mu, diEqual to 1; the area of the irrigation area group is more than 5 ten thousand mu and less than or equal to 10 ten thousand mu, diEqual to 2; the area of the irrigation area group is more than 10 ten thousand mu and less than or equal to 20 ten thousand mu, diEqual to 3;
the irrigation area of the rotation irrigation group I is controlled to be 10.7 ten thousand mu, and the irrigation days are 2 days; the irrigation area of the rotation irrigation group II is controlled to be 7.1 ten thousand mu, and the irrigation days are 2 days; and the irrigation area of the rotation irrigation group III is controlled to be 16.2 ten thousand mu, and the irrigation days are 3 days.
e. Calculating the total irrigation time d of the irrigation areas：
Figure BDA0002044538220000054
This time ds＝7
f. Determining the sequence and date of irrigation of each rotation irrigation group;
the irrigation time of the rotation irrigation group I is determined to be 7-8 days in the current time, the irrigation time of the rotation irrigation group II is determined to be 7-9-7-10 days in the current time, and the irrigation time of the rotation irrigation group I is determined to be 7-11-7-13 days in the current time.
(4) According to the steps of the specification and the flow chart:
a. acquiring all m weather forecast sites d in a circular area with the center of an irrigation area as the center of a circle and the radius of 100 kilometerssForecast data y of precipitation amount of each day in the dayi,j i＝1,2,...,ds；j＝1,2,...,m；
The number of the weather forecast stations is 3 in the circular area with the radius of 100 kilometers, namely m is 3, and the rainfall forecast value of each weather station within 7 days is shown in a table.
TABLE 1 weather station 7 days precipitation forecast (mm)
Weather station 7 month and 7 days 7 month and 8 days 7 month and 9 days 7 month and 10 days 7 month and 11 days 7 month and 12 days 7 month and 13 days
1 0 1.2 3.0 0 0 1.8 2.0
2 0 2.1 0 0 0 2.0 2.0
3 0 3.2 1.0 0 0 1.1 2.0
b. Calculating the distance r from the center of the irrigation area to m weather forecast stationsmAnd a total distance rs：
Figure BDA0002044538220000061
The distances from the 3 weather forecast stations to the center of the irrigation area are 20 kilometers, 30 kilometers and 50 kilometers respectively.
c. Calculating the weight omega of the rainfall of each meteorological site in the average rainfall forecast data of the irrigation districtj,j＝1,2,...,m：
Figure BDA0002044538220000062
This time IIIThe weight of each weather station is omega1＝0.2，ω2＝0.3，ω3＝0.5，
d. Calculating irrigation area dsAverage precipitation forecast data T of each day in the dayi,i＝1,2,...,ds：
Figure BDA0002044538220000063
According to the calculation, the average precipitation forecast data of the irrigation at this time is as follows:
T1＝0.0mm,T2＝2.5mm,T3＝1.1mm,T4＝0.0mm,T5＝0.0mm,T6＝1.5mm,T7＝2.0mm
(5) according to the steps of the specification and the flow chart:
a. the depth H of the water layer in the initial field of irrigation is predicted according to the average precipitation of the irrigation area0Field leakage DP and crop water demand ET, calculating irrigation d by using the following formulasDaily field water layer depth Hi：
Hi＝Hi-1+Ti-ET-DP,i＝1,2,...,ds(ii) a Surface water layers are present on days i, i-1;
or
Hi＝-Hi-1+Ti-ET-DP,i＝1,2,...,ds(ii) a A surface water layer exists on the ith day, and no surface water layer exists on the ith-1 day;
or
Hi＝|Hi-1+Ti-ET-DP|,i＝1,2,...,ds(ii) a No surface water layer on day i, and surface water layer on day i-1;
or
Hi＝|-Hi-1+Ti-ET-DP|,i＝1,2,...,ds(ii) a No surface water layer exists on days i, i-1;
according to the above formula, the water depth change in 7 days of the current irrigation plan is calculated as shown in table 2.
TABLE 2 Water layer depth (mm)
Figure BDA0002044538220000071
b. When not irrigating, the depth of the field water layer of each rotation irrigation group is equal.
(6) According to the steps of the specification and the flow chart:
a. comparing the depth H of the water layer in the field of the ith day of the rotation irrigation group when the surface water layer existsiAnd hmaxThe size of (d); if H isiGreater than hmaxThe rotation irrigation group needs to drain water on the ith day with the water discharge hp＝Hi-hmax(ii) a If H isiH is less than or equal tomaxAnd when the irrigation of the group is turned to the rotation irrigation plan, irrigating the group on the ith day until the water layer in the field reaches hmaxThen, the group is excluded from the irrigation plan and is not irrigated any more; otherwise the group does not need irrigation nor drainage on day i;
b. comparing the water layer depth H of the ith day without surface water layeriAnd hsmaxThe size of (d); if H isiIs greater than or equal to hsmaxThe group must be irrigated on day i to an upper limit of a suitable water layer hmax(ii) a If H issiLess than hsmaxAnd the irrigation of the group is not turned to in the rotation irrigation plan, the group does not need irrigation on the ith day; if H issiLess than hsmaxAnd the irrigation of the group is turned to in the rotation irrigation plan, the group is irrigated on the ith day until the water layer in the field reaches hmaxThen, the group is excluded from the irrigation plan and is not irrigated any more;
e. repeating the steps a and b until all the rotation irrigation groups complete the irrigation plan.
By comparison, the depth of the irrigation water layer in no one day exceeds h in the irrigation planmaxI.e. forced drainage is not required a day. According to the rotation irrigation plan, the rotation irrigation group I needs to irrigate in 7 months and 8 days until the depth of a surface water layer reaches 20 millimeters, and the needed irrigation water is 14.3 millimeters; the rotation irrigation group II needs to irrigate water in 7 months and 10 days until the depth of a surface water layer reaches 20 millimeters, and the needed irrigation water is 30 millimeters; the rotation irrigation group III needs to be irrigated to the depth of a surface water layer of 20mm in 7 months and 11 days, and the irrigation water needs to be 38.4 mm. Up to this irrigation areaThe irrigation plan is completed, and no further irrigation or drainage is needed from 7-month-12 days to 7-month-13 days.
(7) According to the steps of the specification and the flow chart:
and (4) repeating the steps (1) to (6) until the whole growth period of the rice is finished.
Repeating the processes, and continuing to execute the irrigation plan in the rice growth period; until the rice is harvested.

Claims (1)

1. A large irrigation area paddy field irrigation plan optimization method based on rainfall forecast is characterized by comprising the following steps:
(1) determining parameters of the required irrigation requirements according to the growth period, the irrigation mode, the irrigation time period and the irrigation area location of the rice;
(2) measuring the real-time depth of the water layer of the rice field in the irrigation area one day before the initial day of irrigation, and determining the initial depth of the irrigation water layer;
(3) dividing the irrigation area into irrigation rotation groups according to the distribution of main canal gates, and determining the time d required by each group for irrigationiAnd the time required to complete all irrigation dsWherein d issLess than or equal to 7;
(4) obtaining optimal irrigation dsForecasting data of daily average precipitation in the intra-day irrigation area;
(5) calculating irrigation d by using a water balance equation according to the average precipitation forecast data of the irrigation area, the initial irrigation field water layer depth, the field leakage and the crop water demandsThe water layer depth of the field is gradually increased day by day in the irrigation area in the day;
(6) completing the irrigation plan of each rotation irrigation group according to the field water layer depth and the irrigation requirements day by day;
(7) repeating the steps (1) to (6) until the whole growth period of the rice is finished;
the step (1) specifically comprises the following steps:
a. determining irrigation area sites, irrigation time periods, irrigation modes and growth periods of rice in the area during irrigation, wherein the irrigation area sites and the irrigation time periods need to be optimized;
b. determining the upper limit h of the suitable water layer of the irrigation area according to the information in the step amaxSuitable lower limit of water layer hminLower limit of soil moisture content thetaminDry volume weight of soil
Figure FDA0003178869140000011
Field seepage DP, water surface evaporation Z and alpha values, and field water holding capacity theta0The depth of the active layer S of the root system of the rice;
the α values are given by the following table:
Figure FDA0003178869140000012
c. according to the growth period of the rice, the water surface evaporation Z is converted into the crop water demand ET by an alpha value method, namely the following formula:
ET＝αZ；
d. according to the lower limit theta of the water content of the soilminDry volume weight of soil
Figure FDA0003178869140000013
The lower limit h of the (underground) wet water layer is calculated by the following formulasmax：
Figure FDA0003178869140000014
Where ρ is the density of water, ρ is 1g/cm3；
The step (2) specifically comprises the following steps:
a. determining the initial day of irrigation;
b. actually measuring the depth H of the water layer in the irrigation area one day before the initial day of irrigation0；
The step (3) specifically comprises the following steps:
a. determining the area of an irrigation area;
b. dividing the irrigation area into s rotation irrigation groups according to the area of the irrigation area and the distribution condition of main canal gates, wherein the area of each group is not more than 20 ten thousand mu;
c. determining the days d required by each group according to the area of the rotation irrigation groupi(ii) a The area of the irrigation block is less than or equal to 5 ten thousand mu, diEqual to 1; the area of the irrigation area group is more than 5 ten thousand mu and less than or equal to 10 ten thousand mu, diEqual to 2; the area of the irrigation area group is more than 10 ten thousand mu and less than or equal to 20 ten thousand mu, diEqual to 3;
e. calculating the total irrigation time d of the irrigation areas：
Figure FDA0003178869140000021
dsIt is required to be 7 or less.
f. Determining the sequence and date of irrigation of each rotation irrigation group;
the step (4) specifically comprises the following steps:
a. acquiring all m weather forecast sites d in a circular area with the center of an irrigation area as the center of a circle and the radius of 100 kilometerssForecast data y of precipitation amount of each day in the dayi,j i＝1,2,...,ds；j＝1,2,...,m；
b. Calculating the distance r from the center of the irrigation area to m weather forecast stationsmAnd a total distance rs：
Figure FDA0003178869140000022
c. Calculating the weight omega of the rainfall of each meteorological site in the average rainfall forecast data of the irrigation districtj,j＝1,2,...,m：
Figure FDA0003178869140000023
d. Calculating irrigation area dsAverage precipitation forecast data T of each day in the dayi,i＝1,2,...,ds：
Figure FDA0003178869140000024
The step (5) specifically comprises the following steps:
a. the depth H of the water layer in the initial field of irrigation is predicted according to the average precipitation of the irrigation area0Field leakage DP and crop water demand ET, calculating irrigation d by using the following formulasDaily field water layer depth HiWhen no surface water layer exists, the depth of the water layer refers to the distance from the surface of the underground water layer to the surface of the rice field;
Hi＝Hi-1+Ti-ET-DP,i＝1,2,...,ds(ii) a Surface water layers are present on days i, i-1;
or
Hi＝-Hi-1+Ti-ET-DP,i＝1,2,...,ds(ii) a A surface water layer exists on the ith day, and no surface water layer exists on the ith-1 day;
or
Hi＝|Hi-1+Ti-ET-DP|,i＝1,2,...,ds(ii) a No surface water layer on day i, and surface water layer on day i-1;
or
Hi＝|-Hi-1+Ti-ET-DP|,i＝1,2,...,ds(ii) a No surface water layer exists on days i, i-1;
b. when not irrigating, the depth of the field water layer of each rotation irrigation group is equal;
the step (6) specifically comprises the following steps:
a. comparing the depth H of the water layer in the field of the ith day of the rotation irrigation group when the surface water layer existsiAnd hmaxThe size of (d); if H isiGreater than hmaxThe rotation irrigation group needs to drain water on the ith day with the water discharge hp＝Hi-hmax(ii) a If H isiH is less than or equal tomaxAnd when the irrigation of the group is turned to the rotation irrigation plan, the group is irrigated to the field water layer h on the ith daymaxThen, the group is excluded from the irrigation plan and is not irrigated any more; otherwise the group does not need irrigation nor drainage on day i;
b. comparing the water layer depth H of the ith day without surface water layeriAnd hsmaxThe size of (d); if H isiIs greater than or equal to hsmaxThe group must start irrigation on day i to the upper limit of the suitable water layer hmax(ii) a If H isiLess than hsmaxAnd the irrigation of the group is not turned to in the rotation irrigation plan, the group does not need irrigation on the ith day; if H isiLess than hsmaxBut the irrigation of the group is turned to in the rotation irrigation plan, the group is irrigated to reach the depth h of the field water layer on the ith daymaxThen, the group is excluded from the irrigation plan and is not irrigated any more;
e. repeating the steps a and b until all the rotation irrigation groups complete the irrigation plan.
CN201910352946.7A 2019-04-29 2019-04-29 Rainfall forecast-based large irrigation area paddy field irrigation plan optimization method Active CN110419415B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910352946.7A CN110419415B (en) 2019-04-29 2019-04-29 Rainfall forecast-based large irrigation area paddy field irrigation plan optimization method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910352946.7A CN110419415B (en) 2019-04-29 2019-04-29 Rainfall forecast-based large irrigation area paddy field irrigation plan optimization method

Publications (2)

Publication Number Publication Date
CN110419415A CN110419415A (en) 2019-11-08
CN110419415B true CN110419415B (en) 2021-10-12

Family

ID=68408410

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910352946.7A Active CN110419415B (en) 2019-04-29 2019-04-29 Rainfall forecast-based large irrigation area paddy field irrigation plan optimization method

Country Status (1)

Country Link
CN (1) CN110419415B (en)

Citations (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103477927A (en) * 2013-09-24 2014-01-01 曾敏 Rice cultivation water management method for moderate and mild arsenic pollution paddy field
CN104855195A (en) * 2015-04-27 2015-08-26 扬州大学 Rice water-saving irrigation method
CN106845808A (en) * 2017-01-10 2017-06-13 北京师范大学 Intelligently decision-making technique and system are arranged in filling in irrigated area rice field based on remotely-sensed data inverting
CN106934534A (en) * 2017-02-28 2017-07-07 武汉大学 Rattan knot melon Irrigation Project Design water effective utilization coefficients computational methods long
CN107122901A (en) * 2017-04-26 2017-09-01 中国水利水电科学研究院 A kind of determination method of the crucial growing stage of region of no relief field-crop draining
CN107330804A (en) * 2017-07-06 2017-11-07 贵州省水利科学研究院 A kind of wisdom water conservancy management and control cloud platform and method
WO2018009482A1 (en) * 2016-07-05 2018-01-11 Sostena, Inc. System and method for crop management
CN108401854A (en) * 2018-01-30 2018-08-17 华南农业大学 A kind of navel orange water-saving irrigation method based on soil moisture detection
CN108575673A (en) * 2018-01-17 2018-09-28 中国农业科学院农业资源与农业区划研究所 Drought-hit area crop irrigation fertilizing method based on same day Weather Forecast Information and system
CN108876005A (en) * 2018-05-07 2018-11-23 中国农业科学院农田灌溉研究所 Irrigation in winter wheat forecasting procedure based on Weather information
CN108920429A (en) * 2018-06-12 2018-11-30 河海大学 A kind of abnormal data analysis method of Water level trend monitoring
CN109447426A (en) * 2018-10-12 2019-03-08 河海大学 Response analysis method of the crop structure based on crop water mechanism to changing environment

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
SG10201704222VA (en) * 2017-05-24 2018-12-28 Remote Grid Pte Ltd System, method and apparatus for management of agricultural resource

Patent Citations (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103477927A (en) * 2013-09-24 2014-01-01 曾敏 Rice cultivation water management method for moderate and mild arsenic pollution paddy field
CN104855195A (en) * 2015-04-27 2015-08-26 扬州大学 Rice water-saving irrigation method
WO2018009482A1 (en) * 2016-07-05 2018-01-11 Sostena, Inc. System and method for crop management
CN106845808A (en) * 2017-01-10 2017-06-13 北京师范大学 Intelligently decision-making technique and system are arranged in filling in irrigated area rice field based on remotely-sensed data inverting
CN106934534A (en) * 2017-02-28 2017-07-07 武汉大学 Rattan knot melon Irrigation Project Design water effective utilization coefficients computational methods long
CN107122901A (en) * 2017-04-26 2017-09-01 中国水利水电科学研究院 A kind of determination method of the crucial growing stage of region of no relief field-crop draining
CN107330804A (en) * 2017-07-06 2017-11-07 贵州省水利科学研究院 A kind of wisdom water conservancy management and control cloud platform and method
CN108575673A (en) * 2018-01-17 2018-09-28 中国农业科学院农业资源与农业区划研究所 Drought-hit area crop irrigation fertilizing method based on same day Weather Forecast Information and system
CN108401854A (en) * 2018-01-30 2018-08-17 华南农业大学 A kind of navel orange water-saving irrigation method based on soil moisture detection
CN108876005A (en) * 2018-05-07 2018-11-23 中国农业科学院农田灌溉研究所 Irrigation in winter wheat forecasting procedure based on Weather information
CN108920429A (en) * 2018-06-12 2018-11-30 河海大学 A kind of abnormal data analysis method of Water level trend monitoring
CN109447426A (en) * 2018-10-12 2019-03-08 河海大学 Response analysis method of the crop structure based on crop water mechanism to changing environment

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
基于RZWQM模型的石羊河流域春小麦灌溉制度优化;周始威等;《农业工程学报》;20160331(第6期);第121-129页 *
轮灌分组灌溉优化模型与二维编码的遗传算法实现;李彬等;《节水灌溉》;20060331(第3期);第5-6页 *

Also Published As

Publication number Publication date
CN110419415A (en) 2019-11-08

Similar Documents

Publication Publication Date Title
CN1324949C (en) Insufficient irrigation forecast and control method
CN113407897A (en) Design method of distributed water circulation model based on multi-source mutual-aid water supply mode
CN109452146B (en) Winter wheat water-saving irrigation decision method, control device and control system
CN105528734A (en) Water volume distribution method based on multilevel water resource management
CN105850674B (en) A kind of rice field oxygenation fills row's analog control system and method
CN105494033B (en) A kind of intelligent water-saving irrigation method based on crop demand
CN111280019A (en) Soil moisture digital prediction and irrigation early warning method
CN110288169B (en) Partitioned optimization method for crop irrigation system of underground water shallow-buried area
CN112819332A (en) Water distribution method and device based on full-channel transmission and distribution and computer equipment
CN112700035A (en) Optimization method for regional scale crop zoning water and fertilizer management mode
CN114781217B (en) Hydraulic design method, device and equipment for movable drip irrigation system of circular sprinkler
Liu et al. Paddy water managements for diffuse nitrogen and phosphorus pollution control in China: A comprehensive review and emerging prospects
CN110419415B (en) Rainfall forecast-based large irrigation area paddy field irrigation plan optimization method
CN103374904B (en) Farmland discharging-irrigating patten transformation saves land water-saving method
CN105340658A (en) Cultivating device for measuring water consumption of water-saving and drought-resistant rice and use method thereof
CN111587751A (en) Annual water management method for mechanized planting of rice-ratoon rice in winter paddy field
CN115222152B (en) Rotation irrigation system optimization method for improving field drip irrigation uniformity
CN112369176B (en) Water and fertilizer integrated fertilization method based on crop root growth distribution rule
CN110050666B (en) Rainfall forecast-based irrigation optimization method for small electromechanical rice irrigation areas
CN101904292A (en) Soilless cultivation method using waste compost as substrate
CN113283717B (en) Evaluation method for field irrigation water-saving technology
CN110199605B (en) Construction method of ridge seepage-proofing structure for water retention of paddy field
CN108112432B (en) Method for planting rice in sandy soil
CN207011349U (en) A kind of blueberry is taken shelter from rain, irrigates integrated canopy body structure
CN111582766A (en) System and method for alternately irrigating plants according to environmental parameters

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant