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 PDFInfo
- 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
Links
- 238000003973 irrigation Methods 0.000 title claims abstract description 271
- 230000002262 irrigation Effects 0.000 title claims abstract description 263
- 238000000034 method Methods 0.000 title claims abstract description 19
- 238000005457 optimization Methods 0.000 title claims abstract description 7
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 99
- 241000209094 Oryza Species 0.000 claims abstract description 36
- 235000007164 Oryza sativa Nutrition 0.000 claims abstract description 36
- 235000009566 rice Nutrition 0.000 claims abstract description 36
- 238000001556 precipitation Methods 0.000 claims abstract description 19
- 239000003621 irrigation water Substances 0.000 claims abstract description 15
- 238000009826 distribution Methods 0.000 claims abstract description 7
- 239000002352 surface water Substances 0.000 claims description 29
- 239000002689 soil Substances 0.000 claims description 13
- 230000008020 evaporation Effects 0.000 claims description 6
- 238000001704 evaporation Methods 0.000 claims description 6
- 239000005457 ice water Substances 0.000 claims description 3
- 239000002699 waste material Substances 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000005484 gravity Effects 0.000 description 1
- 239000003673 groundwater Substances 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
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
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 soilField 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:
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 soilThe lower limit h of the (underground) wet water layer is calculated by the following formulasmax;
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:
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:
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:
d. Calculating irrigation area dsAverage precipitation forecast data T of each day in the dayi,i=1,2,...,ds:
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 soilField 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 soilThe lower limit h of the (underground) wet water layer is calculated by the following formulasmax;
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:
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)
|
7 month and 7 |
7 month and 8 |
7 month and 9 |
7 month and 10 |
7 month and 11 |
7 month and 12 |
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:
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:
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:
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)
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 soilField 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:
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 soilThe lower limit h of the (underground) wet water layer is calculated by the following formulasmax:
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:
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:
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:
d. Calculating irrigation area dsAverage precipitation forecast data T of each day in the dayi,i=1,2,...,ds:
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.
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)
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)
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 |
-
2019
- 2019-04-29 CN CN201910352946.7A patent/CN110419415B/en active Active
Patent Citations (12)
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)
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 |