CN114219674B - Method and system for determining maximum allowable breeding scale of livestock and poultry in drainage basin - Google Patents

Method and system for determining maximum allowable breeding scale of livestock and poultry in drainage basin Download PDF

Info

Publication number
CN114219674B
CN114219674B CN202111538185.8A CN202111538185A CN114219674B CN 114219674 B CN114219674 B CN 114219674B CN 202111538185 A CN202111538185 A CN 202111538185A CN 114219674 B CN114219674 B CN 114219674B
Authority
CN
China
Prior art keywords
phosphorus
target
basin
load
watershed
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
CN202111538185.8A
Other languages
Chinese (zh)
Other versions
CN114219674A (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.)
Heibei Agricultural University
Institute of Agricultural Resources and Regional Planning of CAAS
Original Assignee
Heibei Agricultural University
Institute of Agricultural Resources and Regional Planning of CAAS
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 Heibei Agricultural University, Institute of Agricultural Resources and Regional Planning of CAAS filed Critical Heibei Agricultural University
Priority to CN202111538185.8A priority Critical patent/CN114219674B/en
Publication of CN114219674A publication Critical patent/CN114219674A/en
Application granted granted Critical
Publication of CN114219674B publication Critical patent/CN114219674B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION 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/00Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/02Agriculture; Fishing; Forestry; Mining
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/18Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION 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/00Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/10Services
    • G06Q50/26Government or public services

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Business, Economics & Management (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Data Mining & Analysis (AREA)
  • Tourism & Hospitality (AREA)
  • Pure & Applied Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Primary Health Care (AREA)
  • Strategic Management (AREA)
  • Human Resources & Organizations (AREA)
  • General Business, Economics & Management (AREA)
  • General Health & Medical Sciences (AREA)
  • Economics (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Marketing (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Analysis (AREA)
  • Computational Mathematics (AREA)
  • Probability & Statistics with Applications (AREA)
  • Software Systems (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Educational Administration (AREA)
  • Evolutionary Biology (AREA)
  • Marine Sciences & Fisheries (AREA)
  • Animal Husbandry (AREA)
  • Operations Research (AREA)
  • Mining & Mineral Resources (AREA)
  • Development Economics (AREA)
  • Agronomy & Crop Science (AREA)
  • Algebra (AREA)
  • Databases & Information Systems (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • General Engineering & Computer Science (AREA)
  • Fodder In General (AREA)
  • Farming Of Fish And Shellfish (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)
  • Housing For Livestock And Birds (AREA)

Abstract

The invention provides a method and a system for determining the maximum allowable breeding scale of livestock and poultry in a drainage basin, wherein the allowable load of phosphorus in a drainage basin farmland is determined according to the allowable load of phosphorus in surface water of the drainage basin, the current loads of life sources of the drainage basin and point source phosphorus; determining the minimum load of phosphorus of the drainage farmland according to the crop planting area, the target yield, the phosphorus loss coefficient of the lowest unit yield and the phosphorus attenuation coefficient of the drainage river network; judging whether the minimum load of the phosphorus of the drainage basin farmland is less than or equal to the allowable load of the phosphorus of the drainage basin farmland; if yes, determining the phosphorus demand of the crops based on the target yield, the planting area and the phosphorus content in the harvested products; determining the equivalent weight of the livestock and poultry breeding pigs according to the phosphorus demand of various crops, the proportion of the yield contributed by fertilization to the target yield, the optimal organic-inorganic ratio, the seasonal utilization rate of organic fertilizer application and the phosphorus emission coefficient of unit pig equivalent weight; if not, adjusting the planting area of the crops with high phosphorus emission until the conditions are met. The invention gives consideration to the growth requirements of crops and water environment protection.

Description

Method and system for determining maximum allowable breeding scale of livestock and poultry in drainage basin
Technical Field
The invention relates to the field of agricultural planting, in particular to a method and a system for determining the maximum allowable breeding scale of livestock and poultry in a drainage basin.
Background
The first and second national pollution sources general survey results show that the livestock and poultry cultivation is an important source of water pollution, and effective prevention and control of the pollution of the livestock and poultry cultivation is one of the main problems in the field of ecological environment, wherein the determination of reasonable cultivation scale is one of effective means for standardizing cultivation production and preventing and controlling cultivation pollution. At present, a method for determining the bearing capacity of livestock and poultry manure in farmland soil from the perspective of the safe digestion capacity of the farmland soil exists, but a method for determining the breeding scale of the livestock and poultry in the watershed, which gives consideration to agricultural production and water environment protection, does not exist.
Disclosure of Invention
The invention aims to provide a method and a system for determining the maximum allowable breeding scale of the livestock and poultry in the drainage basin, which not only meet the growth requirements of crops, but also realize water environment protection.
In order to achieve the purpose, the invention provides the following scheme:
the invention provides a method for determining the maximum allowable breeding scale of livestock and poultry in a drainage basin, which comprises the following steps:
s1, determining the phosphorus allowable load of the target basin farmland according to the phosphorus allowable load of the surface water of the target basin, the current phosphorus load of the life source of the target basin and the current phosphorus load of the point source of the target basin;
s2, determining the minimum phosphorus load of the farmland in the target drainage basin according to the planting area of various crops in the target drainage basin, the target yield of various crops in the target drainage basin, the minimum phosphorus loss coefficient per unit yield of various crops in the target drainage basin and the phosphorus attenuation coefficient of the river network in the target drainage basin; the target yield of each crop is the average yield of each crop in nearly three years;
s3, judging whether the minimum phosphorus load of the target watershed farmland is less than or equal to the allowable phosphorus load of the target watershed farmland;
s31, if yes, determining the phosphorus demand of each crop in the target watershed based on the target yield of each crop in the target watershed, the phosphorus content of each crop harvest in the target watershed and the planting area of each crop in the target watershed;
s32, determining the equivalent weight of the livestock and poultry breeding pigs in the target watershed according to the phosphorus demand of various crops in the target watershed, the proportion of the yield contributed by fertilization to the target yield, the optimal organic-inorganic ratio, the in-season utilization rate of organic fertilizer application and the phosphorus emission coefficient of unit pig equivalent weight; the optimal organic-inorganic ratio is the ratio of organic fertilizer and inorganic fertilizer with minimum phosphorus loss in unit yield of crops; the proportion of the yield contributed by fertilization to the target yield is the proportion of the yield contributed by fertilization to the target yield; the current season utilization rate of the organic fertilizer application is the proportion of the amount of phosphorus absorbed by crops in the current season applied organic fertilizer to the amount of phosphorus in the applied organic fertilizer;
and S33, if not, adjusting the planting area of the crop with higher phosphorus load in the farmland in the S2 until the lowest phosphorus load of the target watershed farmland is less than or equal to the allowable phosphorus load of the target watershed farmland, and returning to execute S31-S32.
Optionally, the following formula is adopted to calculate the phosphorus allowable load of the surface water of the target basin:
F t =(QC t -QC 0 )/1000000;
wherein, F t Allowing the phosphorus to be loaded for surface water of the drainage basin, wherein Q is the perennial average flow of the surface water of the drainage basin, C t Is the phosphorus concentration of the target surface water quality of the drainage basin C 0 The concentration of phosphorus in the background water quality in the basin; the concentration of the phosphorus in the background water quality is the concentration of the phosphorus in the river water quality in an area which is less affected by artificial activities;
calculating the current status load of the phosphorus of the living source of the target basin by adopting the following formula:
F h =N h ×E h ×(1-R)/1000;
wherein, F h The current state load of the phosphorus source for watershed life,N h For the current population of the drainage basin, E h The phosphorus emission coefficient is the river basin living source, and R is the attenuation coefficient of river basin network phosphorus;
calculating the current load of the point source phosphorus element in the target drainage basin by adopting the following formula:
F p =W p ×C p ×(1-R)/1000000;
wherein, F p Is the current load of phosphorus element in a drainage basin point source, W p Is the drainage of sewage from drainage basin point source C p The concentration of phosphorus in sewage discharged for a drainage basin point source, and R is the attenuation coefficient of phosphorus in drainage basin river network;
determining the phosphorus allowable load of the farmland in the target basin by adopting the following formula:
Fc=F t -F h -F p
optionally, the minimum phosphorus load of the target watershed farmland is determined by adopting the following formula:
F c0 =∑(A i ×Y i ×L i )×(1-R);
wherein, F c0 Is the minimum load of phosphorus in the farmland in the drainage basin, A i Is the planting area of the i-th class crop in the drainage basin, Y i Target yield for class i crops in the watershed, L i The phosphorus loss coefficient of the i-th class of crops in the watershed is the lowest unit yield, and R is the attenuation coefficient of the phosphorus in the river network in the watershed.
Optionally, the phosphorus requirement of each type of crop in the target watershed is determined using the following formula:
D i =Y i ×C i ×A i
wherein D is i Phosphorus requirement for class i crops in the watershed, Y i Target yield for class i crops in the watershed, C i The content of phosphorus in the harvest of the i-th crop in the watershed, A i The planting area of the ith crop in the drainage basin.
Optionally, the target basin livestock breeding pig equivalent is calculated by adopting the following formula:
Figure BDA0003413179760000031
wherein N is l Pig equivalent for livestock and poultry breeding in target drainage basin, D i Phosphorus requirement for watershed i-th crops, P i The ratio of the yield contributed by fertilization to the target yield for the class i crop of the drainage basin, O i The optimal proportion of organic and inorganic substances for the ith class of crops in the drainage basin, U i In-season utilization of organic fertilizer applications for class i crops in watersheds, E h The phosphorus emission coefficient of river basin unit pig equivalent is shown, and R is the attenuation coefficient of river basin network phosphorus.
In order to achieve the purpose, the invention also provides the following scheme:
a system for determining the maximum allowable breeding scale of livestock and poultry in a river basin, which comprises:
the farmland phosphorus allowable load determining module is used for determining the phosphorus allowable load of the target basin farmland according to the phosphorus allowable load of the surface water of the target basin, the current state load of the phosphorus of the life source of the target basin and the current state load of the phosphorus of the point source of the target basin;
the farmland phosphorus minimum load determining module is used for determining the minimum load of phosphorus of a target basin according to the planting area of various crops in the target basin, the target yield of various crops in the target basin, the phosphorus loss coefficient of minimum unit yield of various crops in the target basin and the phosphorus attenuation coefficient of a target basin network; the target yield of each crop is the average yield of each crop in nearly three years;
the livestock breeding pig equivalent determining module is used for determining the livestock breeding pig equivalent of a target watershed, and specifically comprises the following steps: judging whether the minimum phosphorus load of the target watershed farmland is less than or equal to the allowable phosphorus load of the target watershed farmland, if so, determining the phosphorus demand of various crops in the target watershed based on the target yield of various crops in the target watershed, the phosphorus content in various crop harvests in the target watershed and the planting area of various crops in the target watershed;
determining the equivalent weight of the livestock and poultry breeding pigs in the target watershed according to the phosphorus demand of various crops in the target watershed, the proportion of the yield contributed by fertilization to the target yield, the optimal organic-inorganic ratio, the in-season utilization rate of organic fertilizer application and the phosphorus emission coefficient of unit pig equivalent weight; the optimal organic-inorganic ratio is the ratio of organic fertilizer and inorganic fertilizer with minimum phosphorus loss in unit yield of crops; the proportion of the yield contributed by fertilization to the target yield is the proportion of the yield contributed by fertilization to the target yield; the current season utilization rate of the organic fertilizer application is the proportion of the amount of phosphorus absorbed by crops in the current season applied organic fertilizer to the amount of phosphorus in the applied organic fertilizer;
if not, adjusting the planting area of the crops with higher phosphorus load in the farmland until the lowest phosphorus load of the target watershed farmland is less than or equal to the allowable phosphorus load of the target watershed farmland, and calculating the equivalent of the target watershed livestock and poultry breeding pigs.
Optionally, the farmland phosphorus allowable load determination module comprises:
the watershed surface water phosphorus allowable load determining unit calculates the target watershed surface water phosphorus allowable load by adopting the following formula:
F t =(QC t -QC 0 )/1000000;
wherein, F t Allowing the phosphorus to be loaded for surface water of the drainage basin, wherein Q is the perennial average flow of the surface water of the drainage basin, C t Is the phosphorus concentration of the target surface water quality of the drainage basin C 0 The concentration of phosphorus in the background water quality in the basin; the concentration of the phosphorus in the background water quality is the concentration of the phosphorus in the river water quality in an area which is less affected by artificial activities;
the watershed life source phosphorus status load determining unit calculates the target watershed life source phosphorus status load by adopting the following formula:
F h =N h ×E h ×(1-R)/1000;
wherein, F h Is the current status load of the river basin living source phosphorus, N h For the current population of the drainage basin, E h The phosphorus emission coefficient is the river basin living source, and R is the attenuation coefficient of river basin network phosphorus;
the current load determining unit of the phosphorus element of the basin point source adopts the following formula to calculate the current load of the phosphorus element of the target basin point source:
F p =W p ×C p ×(1-R)/1000000;
wherein, F p Is the current load of phosphorus element in a drainage basin point source, W p Is the drainage of sewage from drainage basin point source C p The concentration of phosphorus in sewage discharged by a drainage basin point source is shown, and R is the attenuation coefficient of the drainage basin river network phosphorus;
the watershed farmland phosphorus allowable load determining unit determines the target watershed farmland phosphorus allowable load by adopting the following formula:
Fc=F t -F h -F p
optionally, the farmland minimum phosphorus load determining module determines the target watershed farmland minimum phosphorus load by using the following formula:
F c0 =∑(A i ×Y i ×L i )×(1-R);
wherein, F c0 Is the minimum load of phosphorus in the farmland in the drainage basin, A i Is the planting area of the i-th class crop in the drainage basin, Y i Target yield for class i crops in the watershed, L i The phosphorus loss coefficient of the i-th class of crops in the watershed is the lowest unit yield, and R is the attenuation coefficient of the phosphorus in the river network in the watershed.
Optionally, the phosphorus requirement determining module of the crops calculates the phosphorus requirements of various crops in the target watershed by using the following formula:
D i =Y i ×C i ×A i
wherein D is i Phosphorus requirement for class i crops in the watershed, Y i Target yield for class i crops in the watershed, C i The content of phosphorus in the harvest of class i crops in the watershed, A i The planting area of the ith crop in the drainage basin.
Optionally, the equivalent weight determining module calculates the equivalent weight of the livestock-poultry-breeding pigs in the target drainage basin by using the following formula:
Figure BDA0003413179760000051
wherein N is l Pig equivalent for livestock and poultry breeding in target drainage basin, D i Phosphorus requirement for watershed i-th crops, P i The ratio of the yield contributed by fertilization to the target yield for the class i crop of the drainage basin, O i The optimal proportion of organic and inorganic substances for the ith class of crops in the drainage basin, U i In-season utilization of organic fertilizer applications for class i crops in watersheds, E h The phosphorus emission coefficient of the river basin unit pig equivalent is shown, and R is the attenuation coefficient of river basin river network phosphorus.
According to the specific embodiment provided by the invention, the invention discloses the following technical effects:
the invention provides a method and a system for determining the maximum allowable breeding scale of livestock and poultry in a drainage basin, wherein the method comprises the following steps: determining the phosphorus allowable load of the target watershed farmland according to the phosphorus allowable load of the surface water of the target watershed, the current-state load of the phosphorus of the life source of the target watershed and the current-state load of the phosphorus of the point source of the target watershed; determining the minimum phosphorus load of a target basin farmland according to the planting area of various crops in the target basin, the target yield of various crops in the target basin, the minimum unit yield phosphorus loss coefficient of various crops in the target basin and the target basin river network phosphorus attenuation coefficient; judging whether the minimum phosphorus load of the target watershed farmland is less than or equal to the allowable phosphorus load of the target watershed farmland; if so, determining the phosphorus demand of various crops in the target watershed based on the target yield of various crops in the target watershed, the phosphorus content in the harvested products of various crops in the target watershed and the planting area of various crops in the target watershed; determining the equivalent weight of the livestock and poultry breeding pigs in the target drainage basin according to the phosphorus demand of various crops in the target drainage basin, the proportion of the yield contributed by fertilization to the target yield, the optimal organic-inorganic ratio, the in-season utilization rate of organic fertilizer application and the phosphorus emission coefficient of unit pig equivalent weight; if not, adjusting the planting area of the crops with higher phosphorus load in the farmland until the minimum phosphorus load of the target watershed farmland is less than or equal to the allowable phosphorus load of the target watershed farmland, and returning to calculate the equivalent of the livestock and poultry breeding pigs. The invention not only meets the growth requirement of crops, but also realizes the water environment protection.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings needed to be used in the embodiments will be briefly described below, and it is obvious that the drawings in the following description are only some embodiments of the present invention, and it is obvious for those skilled in the art to obtain other drawings without inventive exercise.
FIG. 1 is a flow chart of a method for determining the maximum allowable breeding scale of the livestock and poultry in the drainage basin of the invention;
fig. 2 is a schematic block structure diagram of the system for determining the maximum allowable breeding scale of the livestock and poultry in the drainage basin of the invention.
Description of the symbols:
a farmland phosphorus allowable load determining module-1, a farmland phosphorus minimum load determining module-2 and a livestock breeding pig equivalent determining module-3.
Detailed Description
The technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, and not all of the embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
The invention aims to provide a method and a system for determining the maximum allowable breeding scale of the livestock and poultry in the drainage basin, which not only meet the growth requirements of crops, but also realize water environment protection.
In order to make the aforementioned objects, features and advantages of the present invention more comprehensible, the present invention is described in detail with reference to the accompanying drawings and the detailed description thereof.
As shown in figure 1, the method for determining the maximum allowable breeding scale of the livestock and poultry in the drainage basin comprises the following steps:
and S1, determining the phosphorus allowable load of the farmland in the target basin according to the phosphorus allowable load of the surface water of the target basin, the current phosphorus load of the life source of the target basin and the current phosphorus load of the point source of the target basin.
S2, determining the minimum phosphorus load of the farmland in the target watershed according to the planting area of various crops in the target watershed, the target yield of various crops in the target watershed, the minimum unit yield phosphorus loss coefficient of various crops in the target watershed and the target watershed river network phosphorus attenuation coefficient; the target yield of each crop is the average yield of each crop in nearly three years.
And S3, judging whether the minimum load of the phosphorus of the target basin is less than or equal to the allowable load of the phosphorus of the target basin.
And S31, if yes, determining the phosphorus demand of the crops in the target watershed based on the target yield of the crops in the target watershed, the phosphorus content of the crops in the target watershed and the planting area of the crops in the target watershed.
S32, determining the equivalent weight of the livestock and poultry breeding pigs in the target watershed according to the phosphorus demand of various crops in the target watershed, the proportion of the yield contributed by fertilization to the target yield, the optimal organic-inorganic ratio, the in-season utilization rate of organic fertilizer application and the phosphorus emission coefficient of unit pig equivalent weight; the optimal organic-inorganic ratio is the ratio of organic fertilizer and inorganic fertilizer with minimum phosphorus loss in unit yield of crops; the proportion of the yield contributed by fertilization to the target yield is the proportion of the yield contributed by fertilization to the target yield; the current season utilization rate of the organic fertilizer is the proportion of the amount of phosphorus absorbed by crops in the organic fertilizer applied in the current season to the amount of phosphorus in the organic fertilizer applied.
And S33, if not, adjusting the planting area of the crop with higher phosphorus load in the farmland in the S2 until the lowest phosphorus load of the target watershed farmland is less than or equal to the allowable phosphorus load of the target watershed farmland, and returning to execute S31-S32.
Specifically, the allowable load of watershed water pollutants is determined based on the water environment protection target.
Further, the allowable load of the phosphorus element in the surface water of the target basin is calculated by adopting the following formula:
F t =(QC t -QC 0 )/1000000;
wherein, F t The allowable load (t/a) of phosphorus for surface water of the basin, and Q is the annual average flow (m) of the surface water of the basin 3 /a),C t The concentration (mg/L) of phosphorus in the surface water of the drainage basin is the target water quality C 0 The concentration (mg/L) of phosphorus in the background water quality in the river basin is the concentration of phosphorus in the river water quality in an area which is less affected by artificial activities.
And calculating the current state load of the phosphorus of the target basin living source by adopting the following formula.
F h =N h ×E h ×(1-R)/1000;
Wherein, F h The current state load (t/a), N, of the phosphorus source for river basin life h Number of population for current drainage basin (p), E h The phosphorus emission coefficient (kg/p/a) of the watershed living source is obtained through a first/second pollution source general survey coefficient, and R is the attenuation coefficient (dimensionless) of the watershed river network phosphorus and is obtained from a literature result.
Calculating the current load of the point source phosphorus in the target basin by adopting the following formula:
F p =W p ×C p ×(1-R)/1000000;
wherein, F p Is the current load (t/a) of the phosphorus element of a drainage basin point source, W p Is the drainage (m) of point source sewage in drainage basin 3 /a),C p The concentration (mg/L) of phosphorus in sewage discharged from a drainage basin point source is shown, and R is the attenuation coefficient (dimensionless) of river network phosphorus in a drainage basin, and is obtained from the literature results.
Determining the phosphorus allowable load of the target watershed farmland by adopting the following formula:
Fc=F t -F h -F p
further, the minimum phosphorus load of the target watershed farmland is determined by the following formula:
F c0 =∑(A i ×Y i ×L i )×(1-R);
wherein, F c0 Is minimum load (t/a) of phosphorus in the farmland in the watershed, A i Is the area (ha), Y of the class i crop in the drainage basin i Target yield (t/ha/a), L for class i crops in the watershed i For class i in the basinThe phosphorus loss coefficient (kg/kg) of the minimum unit yield of the material, R is the attenuation coefficient (dimensionless) of the phosphorus in the river network of the river basin, and is obtained from the literature results.
Further, the phosphorus demand of various crops in the target watershed is determined by the following formula:
D i =Y i ×C i ×A i
wherein D is i The phosphorus requirement (t/a), Y of the i-th crop in the drainage basin i Target yield (t/ha/a), C for class i crops in the watershed i The content (kg/kg) of phosphorus in the harvest of the i-th class crop in the watershed, A i Is the area (h) of the class i crop in the drainage basin a )。
Further, the equivalent weight of the target basin livestock and poultry breeding pig is calculated by adopting the following formula:
Figure BDA0003413179760000091
wherein N is l Equivalent (h) and D) of pigs for livestock and poultry breeding in target drainage basin i Phosphorus requirement (t/a), P for class i crops in the watershed i The proportion (%) of the yield contributed by fertilization to the target yield for the class i crop of the watershed, obtained from the results of the literature meta-studies, O i The optimal organic-inorganic ratio O is determined for the i-th class crop in the drainage basin, specifically, on the premise of meeting the phosphorus demand of the crop, by the meta-extraction of the literature and the organic-inorganic ratio with the minimum phosphorus loss in the unit yield of the crop i Specifically, the percentage of organic fertilizer nutrients in the phosphorus demand of crops is expressed by (%); u shape i The current season utilization rate (%) for the application of the organic fertilizer of the ith crop in the watershed is obtained through the literature meta-research result; e h Obtaining a phosphorus emission coefficient (kg/h/a) of a drainage basin unit pig equivalent from a first/second pollution source general survey coefficient; r is the attenuation coefficient (dimensionless) of the river network phosphorus in the river basin, and is obtained from the literature results.
In order to achieve the above object, as shown in fig. 2, the present invention further provides a system for determining a maximum allowable breeding scale of livestock and poultry in a drainage basin, the system comprising: the system comprises a farmland phosphorus allowable load determining module 1, a farmland phosphorus minimum load determining module 2 and a livestock breeding pig equivalent determining module 3.
The farmland phosphorus allowable load determining module 1 is used for determining the phosphorus allowable load of the target watershed farmland according to the phosphorus allowable load of the surface water of the target watershed, the current state load of the target watershed life source phosphorus and the current state load of the target watershed point source phosphorus.
The farmland phosphorus minimum load determining module 2 is used for determining the minimum load of phosphorus of a farmland in a target basin according to the planting area of various crops in the target basin, the target yield of various crops in the target basin, the phosphorus loss coefficient of minimum unit yield of various crops in the target basin and the phosphorus attenuation coefficient of river network in the target basin; the target yield of each crop is the average yield of each crop in nearly three years.
The livestock breeding pig equivalent determining module 3 is used for determining the equivalent of the livestock breeding pigs in the target drainage basin, and specifically comprises the following steps: and judging whether the minimum load of the phosphorus of the farmland in the target watershed is less than or equal to the allowable load of the phosphorus of the farmland in the target watershed, if so, determining the phosphorus demand of various crops in the target watershed based on the target yield of various crops in the target watershed, the phosphorus content in the crops harvested in the target watershed and the planting area of various crops in the target watershed.
Determining the equivalent weight of the livestock and poultry breeding pigs in the target watershed according to the phosphorus demand of various crops in the target watershed, the proportion of the yield contributed by fertilization to the target yield, the optimal organic-inorganic ratio, the in-season utilization rate of organic fertilizer application and the phosphorus emission coefficient of unit pig equivalent weight; the optimal organic-inorganic ratio is the ratio of organic fertilizer and inorganic fertilizer with minimum phosphorus loss in unit yield of crops. The proportion of the yield contributed by fertilization to the target yield is the proportion of the yield contributed by fertilization to the target yield; the current season utilization rate of the organic fertilizer is the proportion of the amount of phosphorus absorbed by crops in the organic fertilizer applied in the current season to the amount of phosphorus in the organic fertilizer applied.
If not, adjusting the planting area of the crops with higher phosphorus load in the farmland until the lowest phosphorus load of the target watershed farmland is less than or equal to the allowable phosphorus load of the target watershed farmland, and calculating the equivalent of the target watershed livestock and poultry breeding pigs.
Further, the farmland phosphorus allowable load determination module comprises:
the watershed surface water phosphorus allowable load determining unit calculates the target watershed surface water phosphorus allowable load by adopting the following formula:
F t =(QC t -QC 0 )/1000000;
F t allowing the phosphorus to be loaded for surface water of the drainage basin, wherein Q is the perennial average flow of the surface water of the drainage basin, C t Is the phosphorus concentration of the target surface water quality of the drainage basin C 0 The concentration of phosphorus in the background water quality in the basin; the concentration of the phosphorus in the background water is the concentration of the phosphorus in the river water in an area which is less affected by artificial activities.
The current state load determining unit of the river basin living source phosphorus calculates the current state load of the target river basin living source phosphorus by adopting the following formula:
F h =N h ×E h ×(1-R)/1000;
wherein, F h Is the current status load of the river basin living source phosphorus, N h For the current population of the drainage basin, E h The coefficient of phosphorus emission of the watershed living source is shown, and R is the attenuation coefficient of the watershed river network phosphorus.
The current load determining unit of the phosphorus element of the basin point source adopts the following formula to calculate the current load of the phosphorus element of the target basin point source:
F p =W p ×C p ×(1-R)/1000000;
wherein, F p Is the current load of phosphorus element in a drainage basin point source, W p Is the drainage of sewage from drainage basin point source C p The concentration of phosphorus in sewage discharged by a drainage basin point source is shown, and R is the attenuation coefficient of the drainage basin river network phosphorus.
The watershed farmland phosphorus allowable load determining unit determines the target watershed farmland phosphorus allowable load by adopting the following formula:
Fc=F t -F h -F p
wherein, F c To the eyesStandard basin farmland phosphorus tolerance load, F t For watershed surface water phosphorus tolerance, F h Load of the current status of the river basin biogenic phosphorus, F p The load is the current load of the point source phosphorus in the drainage basin.
Further, the farmland phosphorus minimum load determining module determines the farmland phosphorus minimum load of the target basin by adopting the following formula:
F c0 =∑(A i ×Y i ×L i )×(1-R);
wherein, F c0 Is the minimum load of phosphorus in the farmland in the drainage basin, A i Is the planting area of the i-th class crop in the drainage basin, Y i Target yield for class i crops in the watershed, L i The phosphorus loss coefficient of the i-th class of crops in the watershed is the lowest unit yield, and R is the attenuation coefficient of the phosphorus in the river network in the watershed.
Further, the phosphorus demand determination module of the crops adopts the following formula to calculate the phosphorus demand of various crops in the target watershed:
D i =Y i ×C i ×A i
wherein D is i Phosphorus requirement for class i crops in the watershed, Y i Target yield for class i crops in the watershed, C i The content of phosphorus in the harvest of class i crops in the watershed, A i The planting area of the ith crop in the drainage basin.
Further, the equivalent weight determining module for the livestock and poultry breeding pigs calculates the equivalent weight of the livestock and poultry breeding pigs in the target drainage basin by adopting the following formula:
Figure BDA0003413179760000111
wherein N is l Pig equivalent for livestock and poultry breeding in target drainage basin, D i Phosphorus requirement for watershed i-th crops, P i The ratio of the yield contributed by fertilization to the target yield for the class i crop of the drainage basin, O i The optimal proportion of organic and inorganic substances for the ith class of crops in the drainage basin, U i Is the in-season utilization rate of the application of the organic fertilizer of the ith crop in the watershed,E h the phosphorus emission coefficient of river basin unit pig equivalent is shown, and R is the attenuation coefficient of river basin network phosphorus.
Compared with the traditional method for determining the livestock and poultry breeding scale only from the perspective of agricultural production, the method for determining the livestock and poultry breeding scale in the drainage basin in the planting and breeding complete combination mode is provided in the aspect of agricultural production and drainage basin water environment comprehensively, and the method is a technology for realizing agricultural sustainable development considering both agricultural production and environmental protection.
In the present specification, the embodiments are described in a progressive manner, each embodiment focuses on differences from other embodiments, and the same and similar parts among the embodiments are referred to each other. For the system disclosed by the embodiment, the description is relatively simple because the system corresponds to the method disclosed by the embodiment, and the relevant points can be referred to the method part for description.
The principles and embodiments of the present invention have been described herein using specific examples, which are provided only to help understand the method and the core concept of the present invention; meanwhile, for a person skilled in the art, according to the idea of the present invention, the specific embodiments and the application range may be changed. In view of the above, the present disclosure should not be construed as limiting the invention.

Claims (2)

1. A method for determining the maximum allowable breeding scale of livestock and poultry in a drainage basin is characterized by comprising the following steps:
s1, determining the phosphorus allowable load of the target basin farmland according to the phosphorus allowable load of the surface water of the target basin, the current state load of the phosphorus of the life source of the target basin and the current state load of the point source of the phosphorus of the target basin;
s2, determining the minimum phosphorus load of the farmland in the target watershed according to the planting area of various crops in the target watershed, the target yield of various crops in the target watershed, the minimum unit yield phosphorus loss coefficient of various crops in the target watershed and the target watershed river network phosphorus attenuation coefficient; the target yield of each crop is the average yield of each crop in nearly three years;
s3, judging whether the minimum phosphorus load of the target watershed farmland is less than or equal to the allowable phosphorus load of the target watershed farmland;
s31, if yes, determining the phosphorus demand of each crop in the target watershed based on the target yield of each crop in the target watershed, the phosphorus content of each crop harvest in the target watershed and the planting area of each crop in the target watershed;
s32, determining the equivalent weight of the livestock and poultry breeding pigs in the target watershed according to the phosphorus demand of various crops in the target watershed, the proportion of the yield contributed by fertilization to the target yield, the optimal organic-inorganic ratio, the in-season utilization rate of organic fertilizer application and the phosphorus emission coefficient of unit pig equivalent weight; the optimal organic-inorganic ratio is the ratio of organic fertilizer and inorganic fertilizer with minimum phosphorus loss in unit yield of crops; the proportion of the yield contributed by fertilization to the target yield is the proportion of the yield contributed by fertilization to the target yield; the current season utilization rate of the organic fertilizer application is the proportion of the amount of phosphorus absorbed by crops in the current season applied organic fertilizer to the amount of phosphorus in the applied organic fertilizer;
s33, if not, adjusting the planting area of the crop with higher phosphorus load in the farmland in the S2 until the lowest phosphorus load of the target watershed farmland is less than or equal to the allowable phosphorus load of the target watershed farmland, and returning to execute S31-S32;
calculating the allowable load of the phosphorus in the surface water of the target basin by adopting the following formula:
F t =(QC t -QC 0 )/1000000;
wherein, F t Allowing the phosphorus element to load for surface water in the drainage basin, Q is the annual average flow rate of the surface water in the drainage basin, and C t Is the phosphorus concentration of the target surface water quality of the drainage basin C 0 The concentration of phosphorus in the background water quality in the basin; the concentration of the phosphorus in the background water quality is the concentration of the phosphorus in the river water quality in an area which is less affected by artificial activities;
calculating the current status load of the phosphorus of the living source of the target basin by adopting the following formula:
F h =N h ×E h ×(1-R)/1000;
wherein, F h For generating in the drainage basinActive source phosphorus present load, N h For the current population of the drainage basin, E h The phosphorus emission coefficient is the river basin living source, and R is the attenuation coefficient of river basin network phosphorus;
calculating the current load of the point source phosphorus element in the target drainage basin by adopting the following formula:
F p =W p ×C p ×(1-R)/1000000;
wherein, F p Is the current load of phosphorus element in a drainage basin point source, W p Is the drainage of sewage from drainage basin point source C p The concentration of phosphorus in sewage discharged by a drainage basin point source is shown, and R is the attenuation coefficient of the drainage basin river network phosphorus;
determining the phosphorus allowable load of the target watershed farmland by adopting the following formula:
Fc=F t -F h -F p
determining the minimum phosphorus load of the target watershed farmland by adopting the following formula:
F c0 =∑(A i ×Y i ×L i )×(1-R);
wherein, F c0 Is the minimum load of phosphorus in the farmland in the drainage basin, A i Is the planting area of the i-th class crop in the drainage basin, Y i Target yield for class i crops in the watershed, L i The phosphorus loss coefficient of the i-th class of crops in the drainage basin per unit yield is shown, and R is the attenuation coefficient of the phosphorus in river network in the drainage basin;
determining the phosphorus demand of various crops in the target watershed by adopting the following formula:
D i =Y i ×C i ×A i
wherein D is i Phosphorus requirement, Y, for class i crops in the watershed i Target yield for class i crops in the watershed, C i The content of phosphorus in the harvest of class i crops in the watershed, A i The planting area of the ith crop in the drainage basin is shown;
calculating the equivalent of the target basin livestock breeding pig by adopting the following formula:
Figure FDA0003779637720000021
wherein N is l Pig equivalent for livestock and poultry breeding in target drainage basin, D i Phosphorus requirement, P, for class i crops in the watershed i The ratio of the yield contributed by fertilization to the target yield for the class i crop of the drainage basin, O i The optimal proportion of organic and inorganic substances for the ith class of crops in the drainage basin, U i In-season utilization of organic fertilizer applications for class i crops in watersheds, E h The phosphorus emission coefficient of river basin unit pig equivalent is shown, and R is the attenuation coefficient of river basin network phosphorus.
2. A system for determining the maximum allowable breeding scale of livestock and poultry in a drainage basin is characterized by comprising:
the farmland phosphorus allowable load determining module is used for determining the phosphorus allowable load of the target basin farmland according to the phosphorus allowable load of the surface water of the target basin, the current state load of the phosphorus of the life source of the target basin and the current state load of the phosphorus of the point source of the target basin;
the farmland phosphorus minimum load determining module is used for determining the minimum load of phosphorus of a target basin according to the planting area of various crops in the target basin, the target yield of various crops in the target basin, the phosphorus loss coefficient of minimum unit yield of various crops in the target basin and the phosphorus attenuation coefficient of a target basin network; the target yield of each crop is the average yield of each crop in nearly three years;
the livestock and poultry breeding pig equivalent determining module is used for determining the target basin livestock and poultry breeding pig equivalent, and specifically comprises the following steps: judging whether the minimum phosphorus load of the target watershed farmland is less than or equal to the allowable phosphorus load of the target watershed farmland, if so, determining the phosphorus demand of various crops in the target watershed based on the target yield of various crops in the target watershed, the phosphorus content in various crop harvests in the target watershed and the planting area of various crops in the target watershed;
determining the equivalent weight of the livestock and poultry breeding pigs in the target watershed according to the phosphorus demand of various crops in the target watershed, the proportion of the yield contributed by fertilization to the target yield, the optimal organic-inorganic ratio, the in-season utilization rate of organic fertilizer application and the phosphorus emission coefficient of unit pig equivalent weight; the optimal organic-inorganic ratio is the ratio of organic fertilizer and inorganic fertilizer with minimum phosphorus loss in unit yield of crops; the proportion of the yield contributed by fertilization to the target yield is the proportion of the yield contributed by fertilization to the target yield; the current season utilization rate of the organic fertilizer application is the proportion of the amount of phosphorus absorbed by crops in the current season applied organic fertilizer to the amount of phosphorus in the applied organic fertilizer;
if not, adjusting the planting area of crops with higher phosphorus load of the farmland until the minimum phosphorus load of the target watershed farmland is less than or equal to the allowable phosphorus load of the target watershed farmland, and calculating the equivalent weight of the livestock and poultry breeding pigs in the target watershed;
the farmland phosphorus allowable load determining module comprises:
the watershed surface water phosphorus allowable load determining unit calculates the target watershed surface water phosphorus allowable load by adopting the following formula:
F t =(QC t -QC 0 )/1000000;
wherein, F t Allowing the phosphorus element to load for surface water in the drainage basin, Q is the annual average flow rate of the surface water in the drainage basin, and C t Is the phosphorus concentration of the target surface water quality of the drainage basin C 0 The concentration of phosphorus in the background water quality in the basin; the concentration of the phosphorus in the background water quality is the concentration of the phosphorus in the river water quality in an area which is slightly influenced by artificial activities;
the current state load determining unit of the river basin living source phosphorus calculates the current state load of the target river basin living source phosphorus by adopting the following formula:
F h =N h ×E h ×(1-R)/1000;
wherein, F h Is the current status load of the river basin living source phosphorus, N h For the current population of the drainage basin, E h The phosphorus emission coefficient is the river basin living source, and R is the attenuation coefficient of river basin network phosphorus;
the current load determining unit of the phosphorus element of the basin point source adopts the following formula to calculate the current load of the phosphorus element of the target basin point source:
F p =W p ×C p ×(1-R)/1000000;
wherein, F p Is the current load of phosphorus element in a drainage basin point source, W p Is the drainage of sewage from drainage basin point source C p The concentration of phosphorus in sewage discharged by a drainage basin point source is shown, and R is the attenuation coefficient of the drainage basin river network phosphorus;
the watershed farmland phosphorus allowable load determining unit determines the target watershed farmland phosphorus allowable load by adopting the following formula:
Fc=F t -F h -F p
the farmland phosphorus minimum load determining module determines the target watershed farmland phosphorus minimum load by adopting the following formula:
F c0 =∑(A i ×Y i ×L i )×(1-R);
wherein, F c0 Is the minimum load of phosphorus in the farmland in the drainage basin, A i Is the planting area of the i-th class crop in the drainage basin, Y i Target yield for class i crops in the watershed, L i The phosphorus loss coefficient of the i-th class of crops in the drainage basin per unit yield is shown, and R is the attenuation coefficient of the phosphorus in river network in the drainage basin;
the phosphorus demand determination module of the crops adopts the following formula to calculate the phosphorus demand of various crops in the target watershed:
D i =Y i ×C i ×A i
wherein D is i Phosphorus requirement, Y, for class i crops in the watershed i Target yield for class i crops in the watershed, C i The content of phosphorus in the harvest of class i crops in the watershed, A i The planting area of the ith crop in the drainage basin is shown;
the equivalent weight determining module for the livestock and poultry breeding pigs calculates the equivalent weight of the livestock and poultry breeding pigs in the target drainage basin by adopting the following formula:
Figure FDA0003779637720000051
wherein N is l Pig equivalent for livestock and poultry breeding in target drainage basin, D i Phosphorus requirement for class i crops in watershedsQuantity, P i The ratio of the yield contributed by fertilization to the target yield for the class i crop of the drainage basin, O i The optimal proportion of organic and inorganic substances for the ith class of crops in the drainage basin, U i In-season utilization of organic fertilizer applications for class i crops in watersheds, E h The phosphorus emission coefficient of river basin unit pig equivalent is shown, and R is the attenuation coefficient of river basin network phosphorus.
CN202111538185.8A 2021-12-15 2021-12-15 Method and system for determining maximum allowable breeding scale of livestock and poultry in drainage basin Active CN114219674B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111538185.8A CN114219674B (en) 2021-12-15 2021-12-15 Method and system for determining maximum allowable breeding scale of livestock and poultry in drainage basin

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111538185.8A CN114219674B (en) 2021-12-15 2021-12-15 Method and system for determining maximum allowable breeding scale of livestock and poultry in drainage basin

Publications (2)

Publication Number Publication Date
CN114219674A CN114219674A (en) 2022-03-22
CN114219674B true CN114219674B (en) 2022-09-06

Family

ID=80702612

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111538185.8A Active CN114219674B (en) 2021-12-15 2021-12-15 Method and system for determining maximum allowable breeding scale of livestock and poultry in drainage basin

Country Status (1)

Country Link
CN (1) CN114219674B (en)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116485574B (en) * 2023-04-24 2024-01-26 河北农业大学 River basin planting structure determining method, system and equipment based on target water quality requirement

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112946239A (en) * 2021-02-01 2021-06-11 中国农业科学院农业资源与农业区划研究所 Method for measuring and calculating maximum manure amount allowed to be applied to soil based on soil phosphorus environment threshold
CN113158591A (en) * 2021-03-17 2021-07-23 南京信息职业技术学院 Method for determining utilization bearing capacity of river basin land development

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104951986B (en) * 2015-06-03 2018-07-06 中国农业科学院农业资源与农业区划研究所 Basin agricultural area source pollutants enter lake load estimate method
CN106503447B (en) * 2016-10-31 2020-01-14 中科宇图科技股份有限公司 Rapid processing method for non-point source pollution discharge list
CN108409047B (en) * 2018-04-16 2020-05-05 中国科学院南京地理与湖泊研究所 Method for systematically controlling rapid proliferation of blue algae cells in lake spring
CN109081530A (en) * 2018-10-20 2018-12-25 湖南艾布鲁环保科技股份有限公司 A kind of livestock and poultry feces recycle processing system and method
CN110046836A (en) * 2019-05-13 2019-07-23 生态环境部华南环境科学研究所 Region cultivates overall control line accounting method

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112946239A (en) * 2021-02-01 2021-06-11 中国农业科学院农业资源与农业区划研究所 Method for measuring and calculating maximum manure amount allowed to be applied to soil based on soil phosphorus environment threshold
CN113158591A (en) * 2021-03-17 2021-07-23 南京信息职业技术学院 Method for determining utilization bearing capacity of river basin land development

Also Published As

Publication number Publication date
CN114219674A (en) 2022-03-22

Similar Documents

Publication Publication Date Title
CN105532342A (en) Organic ecological farming system for rice field
CN113487079A (en) Method and device for low-carbon layout of urban land utilization scale structure
CN111955387A (en) Polder area paddy field coupling pond and ditch three-level wetland system
uddin Mahtab et al. Population and agricultural land use: towards a sustainable food production system in Bangladesh
CN114219674B (en) Method and system for determining maximum allowable breeding scale of livestock and poultry in drainage basin
Zhang et al. Modern greenhouse culture of juvenile soft-shelled turtle, Pelodiscus sinensis
Jørgensen et al. Nitrogen distribution as affected by stocking density in a combined production system of energy crops and free-range pigs
CN113627092A (en) Method for simulating time-space process of migration and transformation of pollutants in lake wetland
CN113463560A (en) Polder area water circulation paddy field non-point source pollution treatment system
Jeevandas et al. Concerns of groundwater depletion and irrigation efficiency in Punjab agriculture: A micro-level study
Srivastava et al. Integrated farming approach for runoff recycling systems in humid plateau areas of eastern India
CN1915006A (en) Ladder type well field installation, and method of use
CN110009193B (en) Reservoir ecological scheduling method considering suitable growth of downstream agricultural crops
CN107117711A (en) Non-point Source Pollution in Irrigation District emission reduction increases income pattern
HOU et al. “Coastal Grass Belt” as Paradigm for Grass-based Livestock Husbandry around Bohai Bay
Lamm Subsurface drip irrigation and possibilities in alfalfa
UPADHYAYA et al. Land and Water Management Strategies for Improving Agricultural Productivity of Farmers: Land and Water Management Strategies for Improving Agricultural Productivity of Farmers
CN114757405A (en) Irrigated area water resource balance optimization configuration method
Behailu et al. Water harvesting in northern Ethiopia: Environmental, health and socio-economic impacts
Machiwal et al. Land and water resource management in the Thar Desert region
Huang et al. Optimizing field management to promote the ecologicalization of agriculture in loess hilly-gully region, China
Kumar et al. Declining water table scenario in Haryana-A review
CN105777426A (en) Ecotype soil modifier for garden
Quach et al. Assessing the Feasibility of Microgrid Supported Open Hydroponics (MSOHCC) for A Resilient Fresh Food Supply in SIDS
Porter et al. Deconstructing agronomic resource use efficiencies to increase food production

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