CN109858719A - Vineyard energy efficiency calculation method and device based on DEA and SFA - Google Patents
Vineyard energy efficiency calculation method and device based on DEA and SFA Download PDFInfo
- Publication number
- CN109858719A CN109858719A CN201711237222.5A CN201711237222A CN109858719A CN 109858719 A CN109858719 A CN 109858719A CN 201711237222 A CN201711237222 A CN 201711237222A CN 109858719 A CN109858719 A CN 109858719A
- Authority
- CN
- China
- Prior art keywords
- energy value
- value
- vineyard
- energy
- output
- 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.)
- Pending
Links
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The present invention provides a kind of vineyard energy efficiency calculation method based on DEA and SFA, estimates including the use of production function of the SFA model to vineyard, obtains the estimated value of each unknown parameter in production function；According to the estimated value of each unknown parameter, the vineyard in the energy value of the output index within a production cycle of energy value and the vineyard of each input-occupancy-output analysis in a production cycle, the correction value of the energy value of output index is calculated；The energy efficiency in vineyard is calculated using DEA model according to the energy value of each input-occupancy-output analysis of the correction value of the energy value of output index and vineyard within a production cycle.Vineyard energy efficiency calculation method provided by the invention based on DEA and SFA, production function is estimated using SFA model, then the correction value of the energy value of output index is calculated, DEA model is recycled to calculate the energy efficiency in vineyard, random noise is eliminated, the result obtained is stable, accurate and effective.
Description
Technical field
The present invention relates to environment efficiency assessment technology field more particularly to a kind of vineyard energy effects based on DEA and SFA
Rate calculation method and device.
Background technique
In recent years, there is the advantages that product market price is high, good in economic efficiency as the grape of industrial crops, in adjustment agriculture
Industry structure and aspect of increasing farmers' income have played important function；On the other hand, establishment planting technology, southern area are yielded two crops a year
Cultivation technique is captured, and has also laid solid foundation for grape industry development, the enthusiasm of each department Developing Grape production is significantly
It improves.But modern agriculture has the characteristics of heavy dependence fossil energy and other resource inputs, the life of this energy height consumption
Production mode can undoubtedly aggravate the burden of ecological environment；On the other hand, production investment unreasonable structure, energy utilization rate is low, resource
The problems such as waste, also seriously affects the sustainable development of grape industry.
For a long time, academia, China focuses on always wine-growing management for the research emphasis of grape industry and sets
Apply etc., it is less for the research in terms of grape Evaluation of Sustainability.In the prior art, pass through DEA (data
Envelopment analysis, DEA) come calculate vineyard energy efficiency specifically, firstly, determine vineyard production process
In good input-occupancy-output analysis and output index, obtain the data value of each input-occupancy-output analysis and output index, based on each input-occupancy-output analysis and
The data value of output index establishes vineyard energy efficiency evaluation model, according to vineyard energy efficiency evaluation model, calculates Portugal
Grape garden energy efficiency value.
The energy efficiency calculation method in vineyard in the prior art, using DEA model to the energy efficiency in vineyard into
Row calculates, and that there are calculated results is unstable for this method, and investment and output index are single, is not applied for different cultivation modes,
It is the disadvantages of being influenced by enchancement factor, lower so as to cause the accuracy and validity of calculated result.
Summary of the invention
(1) technical problems to be solved
The vineyard energy efficiency calculation method based on DEA and SFA that the object of the present invention is to provide a kind of, solves existing
The calculated result of calculation method in technology is influenced by enchancement factor, accuracy and the lower technical problem of validity.
(2) technical solution
In order to solve the above-mentioned technical problem, on the one hand, the present invention provides a kind of vineyard energy effect based on DEA and SFA
Rate calculation method, comprising:
S10, estimated using production function of the stochastic frontier analysis SFA model to the vineyard, obtain the life
Produce the estimated value of each unknown parameter in function；
S20, each investment according to the estimated value, the vineyard of each unknown parameter within a production cycle
The energy value of the output index of the energy value of index and the vineyard within a production cycle, calculates the output and refers to
The correction value of target energy value；
It is S30, every within a production cycle according to the correction value of the energy value of the output index and the vineyard
The energy value of a input-occupancy-output analysis calculates the energy efficiency in the vineyard using DEA DEA model.
Further, before step S10 further include:
The measurement of each input-occupancy-output analysis of the vineyard within a production cycle is determined according to the cultivation mode of grape
Value, the cultivation mode include at least outdoor cultivation and facility cultivation；
According to the survey of the output index of the measured value of each input-occupancy-output analysis and the vineyard within a production cycle
Magnitude calculates separately the energy value of each input-occupancy-output analysis and the energy value of the output index.
Further, the measurement of each input-occupancy-output analysis of the vineyard according to acquisition within a production cycle
The measured value of the output index of value and the vineyard within a production cycle, calculates separately each input-occupancy-output analysis
Energy value and the output index energy value specifically:
Ei=γi×Mi
Wherein, EiFor the energy value of i-th kind of investment or output index, MiFor i-th kind investment or output index measured value,
γiThe folding energy coefficient of i-th kind of investment or output index.
Further, the input-occupancy-output analysis includes at least:
Steel, iron wire, irrigation water, electricity, chemical fertilizer, pesticide, organic fertilizer, agricultural film and labour.
Further, the production function are as follows:
Lnq=β0+β1lnA+β2lnB+β3lnC+β4lnD+β5lnE+β6lnF+β7lnG+β8lnH+β9lnK+v+u
Wherein, q is the energy value of the output index, β0、β1、β2、β3、β4、β5、β6、β7、β8、β9It is unknown ginseng with u
Number, A are the energy value of the steel, and B is the energy value of the iron wire, and C is the energy value of the irrigation water, and D is the electricity
Energy value, E are the energy value of the chemical fertilizer, and F is the energy value of the pesticide, and G is the energy value of the organic fertilizer, and H is described
The energy value of agricultural film, K are the energy value of the labour, and v is random noise.
Further, the production function are as follows: step S20 specifically:
q1=q/exp (lnq- β '0-β'1lnA-β'2lnB-β'3lnC-β'4lnD-β'5lnE-β'6lnF-β'7lnG-β'8lnH-β'9lnK+u')
Wherein, q1For the correction value of the energy value of the output index, q is the energy value of the output index, β '0、β'1、
β'2、β'3、β'4、β'5、β'6、β'7、β'8、β'9It is the estimated value of unknown parameter with u', A is the energy value of the steel, and B is
The energy value of the iron wire, C are the energy value of the irrigation water, and D is the energy value of the electricity, and E is the energy value of the chemical fertilizer,
F is the energy value of the pesticide, and G is the energy value of the organic fertilizer, and H is the energy value of the agricultural film, and K is the labour's
Energy value.
On the other hand, the present invention provides a kind of vineyard energy efficiency computing device based on DEA and SFA, comprising:
Estimation module is obtained for being estimated using production function of the stochastic frontier analysis SFA model to the vineyard
Take the estimated value of each unknown parameter in the production function；
Correction module, for according to the estimated value of each unknown parameter, the vineyard within a production cycle
Each input-occupancy-output analysis output index within a production cycle of energy value and the vineyard energy value, calculate
The correction value of the energy value of the output index；
Computing module, for the correction value and the vineyard according to the energy value of the output index a production week
The energy value of each input-occupancy-output analysis in phase calculates the energy efficiency in the vineyard using DEA DEA model.
In another aspect, the present invention provide it is a kind of for calculating the electronic equipment of vineyard energy efficiency, comprising:
Memory and processor, the processor and the memory complete mutual communication by bus；It is described to deposit
Reservoir is stored with the program instruction that can be executed by the processor, and it is above-mentioned that the processor calls described program instruction to be able to carry out
Method.
Another aspect, the present invention provide a kind of computer program product, and the computer program product is non-including being stored in
Computer program in transitory computer readable storage medium, the computer program include program instruction, when described program refers to
When order is computer-executed, the computer is made to execute above-mentioned method.
Another aspect, the present invention provide a kind of computer readable storage medium, are stored thereon with computer program, the meter
Calculation machine program realizes above-mentioned method when being executed by processor.
(3) beneficial effect
Vineyard energy efficiency calculation method provided by the invention based on DEA and SFA, using SFA model to production letter
Number is estimated, the correction value of the energy value of output index is then calculated, and DEA model is recycled to calculate the energy effect in vineyard
Rate eliminates random noise, suitable for the calculating of the energy efficiency under different cultivation modes, improves the accuracy of calculated result
And validity.
Detailed description of the invention
Fig. 1 is the vineyard energy efficiency calculation method schematic diagram based on DEA and SFA according to the embodiment of the present invention；
Fig. 2 is the vineyard energy efficiency computing device schematic diagram based on DEA and SFA according to the embodiment of the present invention；
Fig. 3 is provided in an embodiment of the present invention for calculating the structural schematic diagram of the electronic equipment of vineyard energy efficiency.
Specific embodiment
In order to keep the purposes, technical schemes and advantages of the embodiment of the present invention clearer, implement below in conjunction with the present invention
Attached drawing in example, technical scheme in the embodiment of the invention is clearly and completely described, it is clear that described embodiment
It is a part of the embodiment of the present invention, instead of all the embodiments.Based on the embodiment of the present invention, those of ordinary skill in the art
Every other embodiment obtained without making creative work, shall fall within the protection scope of the present invention.
Embodiment 1:
Fig. 1 is the vineyard energy efficiency calculation method schematic diagram based on DEA and SFA according to the embodiment of the present invention, such as
Shown in Fig. 1, the embodiment of the present invention provides a kind of vineyard energy efficiency calculation method based on DEA and SFA, comprising:
S10, estimated using production function of the stochastic frontier analysis SFA model to the vineyard, obtain the life
Produce the estimated value of each unknown parameter in function；
S20, each investment according to the estimated value, the vineyard of each unknown parameter within a production cycle
The energy value of the output index of the energy value of index and the vineyard within a production cycle, calculates the output and refers to
The correction value of target energy value；
It is S30, every within a production cycle according to the correction value of the energy value of the output index and the vineyard
The energy value of a input-occupancy-output analysis calculates the energy efficiency in the vineyard using DEA DEA model.
Further, before step S10 further include:
The measurement of each input-occupancy-output analysis of the vineyard within a production cycle is determined according to the cultivation mode of grape
Value, the cultivation mode include at least outdoor cultivation and facility cultivation；
According to the survey of the output index of the measured value of each input-occupancy-output analysis and the vineyard within a production cycle
Magnitude calculates separately the energy value of each input-occupancy-output analysis and the energy value of the output index.
Further, the measurement of each input-occupancy-output analysis of the vineyard according to acquisition within a production cycle
The measured value of the output index of value and the vineyard within a production cycle, calculates separately each input-occupancy-output analysis
Energy value and the output index energy value specifically:
Ei=γi×Mi
Wherein, EiFor the energy value of i-th kind of investment or output index, MiFor i-th kind investment or output index measured value,
γiThe folding energy coefficient of i-th kind of investment or output index.
Further, the input-occupancy-output analysis includes at least:
Steel, iron wire, irrigation water, electricity, chemical fertilizer, pesticide, organic fertilizer, agricultural film and labour.
Further, the production function are as follows:
Lnq=β0+β1lnA+β2lnB+β3lnC+β4lnD+β5lnE+β6lnF+β7lnG+β8lnH+β9lnK+v+u
Wherein, q is the energy value of the output index, β0、β1、β2、β3、β4、β5、β6、β7、β8、β9It is unknown ginseng with u
Number, A are the energy value of the steel, and B is the energy value of the iron wire, and C is the energy value of the irrigation water, and D is the electricity
Energy value, E are the energy value of the chemical fertilizer, and F is the energy value of the pesticide, and G is the energy value of the organic fertilizer, and H is described
The energy value of agricultural film, K are the energy value of the labour, and v is random noise.
Further, the production function are as follows: step S20 specifically:
q1=q/exp (lnq- β '0-β'1lnA-β'2lnB-β'3lnC-β'4lnD-β'5lnE-β'6lnF-β'7lnG-β'8lnH-β'9lnK+u')
Wherein, q1For the correction value of the energy value of the output index, q is the energy value of the output index, β '0、β'1、
β'2、β'3、β'4、β'5、β'6、β'7、β'8、β'9It is the estimated value of unknown parameter with u', A is the energy value of the steel, and B is
The energy value of the iron wire, C are the energy value of the irrigation water, and D is the energy value of the electricity, and E is the energy value of the chemical fertilizer,
F is the energy value of the pesticide, and G is the energy value of the organic fertilizer, and H is the energy value of the agricultural film, and K is the labour's
Energy value.
Specifically, the embodiment of the present invention is carried out by taking Chinese facility grape creation data in 2015 as an example based on DEA and SFA
Energy efficiency assessment, data source is the questionnaire data of national grape industrial system, obtains 81 vineyaries altogether at one
Input-occupancy-output analysis and output achievement data in production cycle.
Firstly, determining the measurement of each input-occupancy-output analysis of the vineyard within a production cycle according to the cultivation mode of grape
Value, cultivation mode include at least outdoor cultivation and facility cultivation, wherein facility cultivation includes at least Forcing cultivation, delayed cultivation
And rain shelter cultivation.The energy investment of different cultivation modes is different, in practical application, can determine have according to specific cultivation mode
The input-occupancy-output analysis of body.Chinese Grape by investigating acquisition on the spot produces actual conditions, obtains with access pertinent literature, books
Information combines, and the embodiment of the present invention is comprehensive to have chosen steel (Kg/ mus), iron wire (Kg/ mus), fill according to facility cultivation mode
Irrigate water (cubic meter/mu), electric (degree/mu), chemical fertilizer (Kg/ mus), pesticide (Kg/ mus), organic fertilizer (Kg/ mus), agricultural film (Kg/ mus),
Labour's (man day/mu) is used as input-occupancy-output analysis, has chosen (Kg/ mus) of vintage as output index.
It should be noted that the measured value of each input-occupancy-output analysis of the vineyard within a production cycle is by converting it
Data afterwards, vineyard is in the orchard establishment period that the measured value of each input-occupancy-output analysis in a production cycle includes after conversion
The input-occupancy-output analysis disposably put into measured value, that is, including vineyard orchard establishment period investment steel, iron wire, labour, agriculture
Chemical fertilizer, pesticide, organic fertilizer, agricultural film, irrigation water, the electricity, labour that the input of film and vineyard are put within a production cycle
The input of power, and according to vineyard average life, the production cycle in such as 20 years and vineyard, such as 1 year, by orchard establishment
Each index input in period is converted, and is merged in the input with vineyard in a production cycle, is obtained vineyard and exist
Input-occupancy-output analysis and output achievement data in one production cycle, totally 9 input-occupancy-output analysis, including steel, iron wire, irrigation water, electricity,
Chemical fertilizer, pesticide, organic fertilizer, agricultural film and labour, there are one output index, the yield in vineyard.
Complete, detailed energy investment and output index system are constructed, the calculated result of energy efficiency can be made more quasi-
Really, effectively.
Then, the measured value of each input-occupancy-output analysis of the vineyard of acquisition within a production cycle and vineyard exist
The measured value of output index in one production cycle.
Then, domestic and foreign literature is consulted, and domestic grape industry is combined to plant actual conditions, compiles to obtain vineyard
The folding energy coefficient of input and output index passes through the measurement of folding energy coefficient value and the previous step input-occupancy-output analysis obtained and output index
Value Data calculates separately the energy value of each input-occupancy-output analysis and the energy value of the output index.Calculation formula is as follows:
Ei=γi×Mi
Wherein, EiFor the energy value of i-th kind of investment or output index, MiFor i-th kind investment or output index measured value,
γiThe folding energy coefficient of i-th kind of investment or output index.
Then, estimated using production function of the stochastic frontier analysis SFA model to the vineyard, obtain the life
Produce the estimated value of each unknown parameter in function.That is, saving as TXT after the energy Value Data of investment and output index is taken logarithm
File imports in Frontier4.1 software, calculates the estimated value of every unknown parameter in production function.Investment after taking logarithm
It is as shown in table 1 with the energy Value Data of output index.
Table 1 takes the energy Value Data of the investment after logarithm and output index
The production function in vineyard are as follows:
Lnq=β0+β1lnA+β2lnB+β3lnC+β4lnD+β5lnE+β6lnF+β7lnG+β8lnH+β9lnK+v+u
Wherein, q is the energy value of the output index, and β 0, β 1, β 2, β 3, β 4, β 5, β 6, β 7, β 8, β 9 and u are unknown
Parameter, wherein u is unknown random error.A is the energy value of the steel, and B is the energy value of the iron wire, and C is the filling
The energy value of water is irrigate, D is the energy value of the electricity, and E is the energy value of the chemical fertilizer, and F is the energy value of the pesticide, and G is institute
The energy value of organic fertilizer is stated, H is the energy value of the agricultural film, and K is the energy value of the labour, and v is random noise.
The ins document parameter setting of Frontier4.1 is as follows:
1 1=ERROR COMPONENTS MODEL, 2=TE EFFECTS MODEL
test.dta DATA FILE NAME
test.out OUTPUT FILE NAME
1 1=PRODUCTION FUNCTION, 2=COST FUNCTION
y LOGGED DEPENDENT VARIABLE(Y/N)
81 NUMBER OF CROSS-SECTIONS
1 NUMBER OF TIME PERIODS
81 NUMBER OF OBSERVATIONS IN TOTAL
9 NUMBER OF REGRESSOR VARIABLES(Xs)
n MU(Y/N)[OR DELTA0(Y/N)IF USING TE EFFECTS MODEL]
n ETA(Y/N)[OR NUMBER OF TE EFFECTS REGRESSORS(Zs)]
n STARTING VALUES(Y/N)
The result of the estimated value of the unknown parameter obtained after operation is as shown in table 2 and table 3.
The estimated value of 2 unknown parameter β of table
| Parameter beta | Numerical value |
| β0 | 8.6414766 |
| β1 | 0.00111164 |
| β2 | -0.0296534 |
| β3 | -0.0470274 |
| β4 | 0.04891689 |
| β5 | -0.0479722 |
| β6 | 0.13232002 |
| β7 | 0.02741187 |
| β8 | 0.11638066 |
| β9 | 0.06255305 |
The estimated value of 3 unknown parameter u of table
Then, each input-occupancy-output analysis according to the estimated value of each unknown parameter, vineyard within a production cycle
The energy value of the output index of energy value and vineyard within a production cycle calculates repairing for the energy value of output index
Positive value.Specific calculation formula is as follows:
q1=q/exp (lnq- β '0-β'1lnA-β'2lnB-β'3lnC-β'4lnD-β'5lnE-β'6lnF-β'7lnG-β'8lnH-β'9lnK+u')
Wherein, q1For the correction value of the energy value of the output index, q is the energy value of the output index, β '0、β'1、
β'2、β'3、β'4、β'5、β'6、β'7、β'8、β'9It is the estimated value of unknown parameter with u', A is the energy value of the steel, and B is
The energy value of the iron wire, C are the energy value of the irrigation water, and D is the energy value of the electricity, and E is the energy value of the chemical fertilizer,
F is the energy value of the pesticide, and G is the energy value of the organic fertilizer, and H is the energy value of the agricultural film, and K is the labour's
Energy value.
By calculating, the correction value of the energy value of obtained output index is as shown in table 4.
The correction value (unit: MJ/ha) of the energy value of the energy value and output index of 4 input-occupancy-output analysis of table
Finally, being referred to according to each investment of the correction value of the energy value of output index and vineyard within a production cycle
Target energy value calculates the energy efficiency in vineyard using DEA DEA model.That is, the data in table 4 are inputted
The measuring and calculating of DEA method is carried out in software DEAP2.1, parameter setting is as follows:
By the energy efficiency TE1 that vineyard is calculated.In addition, it is calculated to sample DEA model in the prior art
The energy efficiency TE2 in vineyard samples the energy efficiency TE3 in the calculated vineyard of SFA model in the prior art, three kinds of sides
The calculated energy efficiency result of case and comparison are as shown in table 5.
5 energy efficiency calculated result of table and compare
And sufficiently small disturbing factor is generated using Matlab software, to the solution of the present invention (vehicle economy A&SFA), DEA
And the anti-interference ability of SFA method has carried out stability test, the results are shown in Table 6.
6 stability test result of table
By analysis, in the energy efficiency measuring and calculating in above 81 vineyaries, from the point of view of the average value of energy efficiency, TE1>
TE2>TE3, wherein DEA&SFA model and DEA model result are increasingly similar, and having 31 vineyard energy efficiency values to reach 1 (has
Effect), in three kinds of models, DEA&SFA model has effectively eliminated the influence of enchancement factor bring.
From, as can be seen that SFA model is best in three kinds of models, DEA&SFA model is secondly, DEA in stability test result
Model is worst.
In general, DEA&SFA model can eliminate random error that SFA model can not be eliminated and DEA model and have
Better stability can more accurate, efficiently calculate the energy efficiency in facility grape garden.
Vineyard energy efficiency calculation method provided by the invention based on DEA and SFA, using SFA model to production letter
Number is estimated, the correction value of the energy value of output index is then calculated, and DEA model is recycled to calculate the energy effect in vineyard
Rate eliminates random noise, suitable for the calculating of the energy efficiency under different cultivation modes, improves the accuracy of calculated result
And validity.
Embodiment 2:
Fig. 2 is the vineyard energy efficiency computing device schematic diagram based on DEA and SFA according to the embodiment of the present invention, such as
Shown in Fig. 2, the embodiment of the present invention provides a kind of vineyard energy efficiency computing device based on DEA and SFA, for realizing above-mentioned
Method as described in the examples specifically includes estimation module 10, correction module 20 and computing module 30, wherein
Estimation module 10 for being estimated using production function of the stochastic frontier analysis SFA model to the vineyard,
Obtain the estimated value of each unknown parameter in the production function；
Correction module 20 is used for the estimated value according to each unknown parameter, the vineyard within a production cycle
Each input-occupancy-output analysis output index within a production cycle of energy value and the vineyard energy value, calculate
The correction value of the energy value of the output index；
Computing module 30 is used to be produced according to the correction value and the vineyard of the energy value of the output index at one
The energy value of each input-occupancy-output analysis in period calculates the energy effect in the vineyard using DEA DEA model
Rate.
The method and above-mentioned reality calculated using energy efficiency of the device described in the embodiment of the present invention to vineyard
Apply identical in example, details are not described herein again.
Vineyard energy efficiency computing device provided by the invention based on DEA and SFA, using SFA model to production letter
Number is estimated, the correction value of the energy value of output index is then calculated, and DEA model is recycled to calculate the energy effect in vineyard
Rate eliminates random noise, improves the accuracy and validity of calculated result.
Embodiment 3:
Fig. 3 is provided in an embodiment of the present invention for calculating the structural schematic diagram of the electronic equipment of vineyard energy efficiency,
As shown in figure 3, the equipment includes: processor (processor) 801, memory (memory) 802 and bus 803；
Wherein, processor 801 and memory 802 complete mutual communication by the bus 803；
Processor 801 is used to call the program instruction in memory 802, to execute provided by above-mentioned each method embodiment
Method, for example,
S10, estimated using production function of the stochastic frontier analysis SFA model to the vineyard, obtain the life
Produce the estimated value of each unknown parameter in function；
S20, each investment according to the estimated value, the vineyard of each unknown parameter within a production cycle
The energy value of the output index of the energy value of index and the vineyard within a production cycle, calculates the output and refers to
The correction value of target energy value；
It is S30, every within a production cycle according to the correction value of the energy value of the output index and the vineyard
The energy value of a input-occupancy-output analysis calculates the energy efficiency in the vineyard using DEA DEA model.
Embodiment 4:
The embodiment of the present invention discloses a kind of computer program product, and the computer program product is non-transient including being stored in
Computer program on computer readable storage medium, the computer program include program instruction, when described program instructs quilt
When computer executes, computer is able to carry out method provided by above-mentioned each method embodiment, for example,
S10, estimated using production function of the stochastic frontier analysis SFA model to the vineyard, obtain the life
Produce the estimated value of each unknown parameter in function；
S20, each investment according to the estimated value, the vineyard of each unknown parameter within a production cycle
The energy value of the output index of the energy value of index and the vineyard within a production cycle, calculates the output and refers to
The correction value of target energy value；
It is S30, every within a production cycle according to the correction value of the energy value of the output index and the vineyard
The energy value of a input-occupancy-output analysis calculates the energy efficiency in the vineyard using DEA DEA model.
Embodiment 5:
The embodiment of the present invention provides a kind of non-transient computer readable storage medium, the non-transient computer readable storage
Medium storing computer instruction, the computer instruction make the computer execute side provided by above-mentioned each method embodiment
Method, for example,
S10, estimated using production function of the stochastic frontier analysis SFA model to the vineyard, obtain the life
Produce the estimated value of each unknown parameter in function；
S20, each investment according to the estimated value, the vineyard of each unknown parameter within a production cycle
The energy value of the output index of the energy value of index and the vineyard within a production cycle, calculates the output and refers to
The correction value of target energy value；
It is S30, every within a production cycle according to the correction value of the energy value of the output index and the vineyard
The energy value of a input-occupancy-output analysis calculates the energy efficiency in the vineyard using DEA DEA model.
Those of ordinary skill in the art will appreciate that: realize that all or part of the steps of above method embodiment can pass through
The relevant hardware of program instruction is completed, and program above-mentioned can be stored in a computer readable storage medium, the program
When being executed, step including the steps of the foregoing method embodiments is executed；And storage medium above-mentioned includes: ROM, RAM, magnetic disk or light
The various media that can store program code such as disk.
The embodiments such as device and equipment described above are only schematical, wherein described be used as separate part description
Unit may or may not be physically separated, component shown as a unit may or may not be
Physical unit, it can it is in one place, or may be distributed over multiple network units.It can be according to the actual needs
Some or all of the modules therein is selected to achieve the purpose of the solution of this embodiment.Those of ordinary skill in the art are not paying
In the case where creative labor, it can understand and implement.
Through the above description of the embodiments, those skilled in the art can be understood that each embodiment can
It realizes by means of software and necessary general hardware platform, naturally it is also possible to pass through hardware.Based on this understanding, on
Stating technical solution, substantially the part that contributes to existing technology can be embodied in the form of software products in other words, should
Computer software product may be stored in a computer readable storage medium, such as ROM/RAM, magnetic disk, CD, including several fingers
It enables and using so that a computer equipment (can be personal computer, server or the network equipment etc.) executes each implementation
Method described in certain parts of example or embodiment.
Finally, it should be noted that the above embodiments are merely illustrative of the technical solutions of the present invention, rather than its limitations；Although
Present invention has been described in detail with reference to the aforementioned embodiments, those skilled in the art should understand that: it still may be used
To modify the technical solutions described in the foregoing embodiments or equivalent replacement of some of the technical features；
And these are modified or replaceed, technical solution of various embodiments of the present invention that it does not separate the essence of the corresponding technical solution spirit and
Range.
Claims (10)
1. a kind of vineyard energy efficiency calculation method based on DEA and SFA characterized by comprising
S10, estimated using production function of the stochastic frontier analysis SFA model to the vineyard, obtain the production letter
The estimated value of each unknown parameter in number；
S20, each input-occupancy-output analysis according to the estimated value, the vineyard of each unknown parameter within a production cycle
Output index within a production cycle of energy value and the vineyard energy value, calculate the output index
The correction value of energy value；
S30, each throwing according to the correction value and the vineyard of the energy value of the output index within a production cycle
The energy value for entering index calculates the energy efficiency in the vineyard using DEA DEA model.
2. the method according to claim 1, wherein before step S10 further include:
The measured value of each input-occupancy-output analysis of the vineyard within a production cycle, institute are determined according to the cultivation mode of grape
Cultivation mode is stated including at least outdoor cultivation and facility cultivation；
According to the measurement of the output index of the measured value of each input-occupancy-output analysis and the vineyard within a production cycle
Value, calculates separately the energy value of each input-occupancy-output analysis and the energy value of the output index.
3. according to the method described in claim 2, it is characterized in that, the vineyard according to acquisition is a production week
The measured value of output index of the measured value and the vineyard of each input-occupancy-output analysis in phase within a production cycle,
Calculate separately the energy value of each input-occupancy-output analysis and the energy value of the output index specifically:
Ei=γi×Mi
Wherein, EiFor the energy value of i-th kind of investment or output index, MiFor the measured value of i-th kind of investment or output index, γiThe
The folding energy coefficient of i kind investment or output index.
4. according to the method described in claim 3, it is characterized in that, the input-occupancy-output analysis includes at least:
Steel, iron wire, irrigation water, electricity, chemical fertilizer, pesticide, organic fertilizer, agricultural film and labour.
5. according to the method described in claim 4, it is characterized in that, the production function are as follows:
Lnq=β0+β1lnA+β2lnB+β3lnC+β4lnD+β5lnE+β6lnF+β7lnG+β8lnH+β9lnK+v+u
Wherein, q is the energy value of the output index, β0、β1、β2、β3、β4、β5、β6、β7、β8、β9It is unknown parameter with u, A is
The energy value of the steel, B are the energy value of the iron wire, and C is the energy value of the irrigation water, and D is the energy value of the electricity,
E is the energy value of the chemical fertilizer, and F is the energy value of the pesticide, and G is the energy value of the organic fertilizer, and H is the energy of the agricultural film
Magnitude, K are the energy value of the labour, and v is random noise.
6. according to the method described in claim 5, it is characterized in that, the production function are as follows: step S20 specifically:
q1=q/exp (lnq- β '0-β'1lnA-β'2lnB-β'3lnC-β'4lnD-β'5lnE-β'6lnF-β'7lnG-β'8lnH-β'9lnK+u')
Wherein, q1For the correction value of the energy value of the output index, q is the energy value of the output index, β '0、β'1、β'2、
β'3、β'4、β'5、β'6、β'7、β'8、β'9It is the estimated value of unknown parameter with u', A is the energy value of the steel, and B is described
The energy value of iron wire, C are the energy value of the irrigation water, and D is the energy value of the electricity, and E is the energy value of the chemical fertilizer, and F is
The energy value of the pesticide, G are the energy value of the organic fertilizer, and H is the energy value of the agricultural film, and K is the energy of the labour
Magnitude.
7. a kind of vineyard energy efficiency computing device based on DEA and SFA characterized by comprising
Estimation module obtains institute for estimating using production function of the stochastic frontier analysis SFA model to the vineyard
State the estimated value of each unknown parameter in production function；
Correction module, for every within a production cycle according to the estimated value of each unknown parameter, the vineyard
The energy value of the output index of the energy value of a input-occupancy-output analysis and the vineyard within a production cycle, described in calculating
The correction value of the energy value of output index；
Computing module, for the correction value and the vineyard according to the energy value of the output index within a production cycle
The energy value of each input-occupancy-output analysis calculate the energy efficiency in the vineyard using DEA DEA model.
8. a kind of for calculating the electronic equipment of vineyard energy efficiency characterized by comprising
Memory and processor, the processor and the memory complete mutual communication by bus；The memory
It is stored with the program instruction that can be executed by the processor, the processor calls described program instruction to be able to carry out right such as and wants
Seek 1 to 6 any method.
9. a kind of computer program product, which is characterized in that the computer program product includes being stored in non-transient computer
Computer program on readable storage medium storing program for executing, the computer program include program instruction, when described program is instructed by computer
When execution, the computer is made to execute the method as described in claim 1 to 6 is any.
10. a kind of computer readable storage medium, is stored thereon with computer program, which is characterized in that the computer program
The method as described in claim 1 to 6 is any is realized when being executed by processor.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711237222.5A CN109858719A (en) | 2017-11-30 | 2017-11-30 | Vineyard energy efficiency calculation method and device based on DEA and SFA |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711237222.5A CN109858719A (en) | 2017-11-30 | 2017-11-30 | Vineyard energy efficiency calculation method and device based on DEA and SFA |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN109858719A true CN109858719A (en) | 2019-06-07 |
Family
ID=66887921
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201711237222.5A Pending CN109858719A (en) | 2017-11-30 | 2017-11-30 | Vineyard energy efficiency calculation method and device based on DEA and SFA |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109858719A (en) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113869688A (en) * | 2021-09-18 | 2021-12-31 | 深圳市裕展精密科技有限公司 | Equipment efficiency prediction method, electronic device, production scheduling method and device |
| CN118228889A (en) * | 2024-05-11 | 2024-06-21 | 广州华夏汇海科技有限公司 | Method and system for generating college entrance examination football score based on neural network |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US9076118B1 (en) * | 2015-01-23 | 2015-07-07 | Iteris, Inc. | Harvest advisory modeling using field-level analysis of weather conditions, observations and user input of harvest condition states, wherein a predicted harvest condition includes an estimation of standing crop dry-down rates, and an estimation of fuel costs |
| CN106651184A (en) * | 2016-12-26 | 2017-05-10 | 中国农业大学 | Vineyard chemical fertilizer utilization efficiency measurement and calculation method and device |
-
2017
- 2017-11-30 CN CN201711237222.5A patent/CN109858719A/en active Pending
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US9076118B1 (en) * | 2015-01-23 | 2015-07-07 | Iteris, Inc. | Harvest advisory modeling using field-level analysis of weather conditions, observations and user input of harvest condition states, wherein a predicted harvest condition includes an estimation of standing crop dry-down rates, and an estimation of fuel costs |
| CN106651184A (en) * | 2016-12-26 | 2017-05-10 | 中国农业大学 | Vineyard chemical fertilizer utilization efficiency measurement and calculation method and device |
Non-Patent Citations (4)
| Title |
|---|
| ALIREZA KHOSHROO等: ""一种可提高葡萄生产能源效率的非参数数据包络分析方法"", 《ENERGY》 * |
| 戴碧艳等: ""我国葡萄生产能源效率评价——基于数据包络与随机前沿的方法融合及验证"", 《中国农业大学学报》 * |
| 李双杰等: ""随机前沿分析与数据包络分析方法的评析与比较"", 《统计与决策》 * |
| 陈锡康主编: "《农业投入产出技术与模型》", 31 October 1987 * |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113869688A (en) * | 2021-09-18 | 2021-12-31 | 深圳市裕展精密科技有限公司 | Equipment efficiency prediction method, electronic device, production scheduling method and device |
| CN118228889A (en) * | 2024-05-11 | 2024-06-21 | 广州华夏汇海科技有限公司 | Method and system for generating college entrance examination football score based on neural network |
| CN118228889B (en) * | 2024-05-11 | 2024-08-20 | 广州华夏汇海科技有限公司 | Method and system for generating college entrance examination football score based on neural network |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Perry | Accounting for water use: Terminology and implications for saving water and increasing production | |
| Feuerbacher et al. | An analytical framework to estimate the economics and adoption potential of dual land-use systems: The case of agrivoltaics | |
| Zhao et al. | Evapotranspiration estimation methods in hydrological models | |
| CN112348292B (en) | Short-term wind power prediction method and system based on deep learning network | |
| Barradas et al. | A decision support system-fertigation simulator (DSS-FS) for design and optimization of sprinkler and drip irrigation systems | |
| Mushtaq et al. | Climate change and water security: estimating the greenhouse gas costs of achieving water security through investments in modern irrigation technology | |
| Pan et al. | Spatial-temporal dynamics of grain yield and the potential driving factors at the county level in China | |
| Liu et al. | The calculation of productivity factor for ecological footprints in China: A methodological note | |
| CN108510107B (en) | Agricultural planting guidance method, electronic equipment and storage medium | |
| Ma et al. | Will the decline of efficiency in China's agriculture come to an end? An analysis based on opening and convergence | |
| CN106600037B (en) | Multi-parameter auxiliary load prediction method based on principal component analysis | |
| Huang et al. | Assessing grain crop water productivity of China using a hydro-model-coupled-statistics approach: Part I: Method development and validation | |
| CN114881544B (en) | Reservoir flow regulation and control method and device, electronic equipment and storage medium | |
| CN109858719A (en) | Vineyard energy efficiency calculation method and device based on DEA and SFA | |
| Irauschek et al. | Evaluating five forest models using multi-decadal inventory data from mountain forests | |
| CN116187823A (en) | Regional carbon emission calculation method and device based on electric power data | |
| CN112715322B (en) | Method and device for obtaining agricultural irrigation water | |
| CN117236650B (en) | Water-fertilizer integrated intelligent pump house control method | |
| CN113255999A (en) | Multi-model-based crop yield assessment method, device, equipment and medium | |
| Teixeira et al. | Modelling seasonality of dry matter partitioning and root maintenance respiration in lucerne (Medicago sativa L.) crops | |
| Verburg et al. | Why benefits from controlled release fertilisers can be lower than expected on some soils | |
| CN112541611A (en) | Method and system for predicting pollution discharge amount of rain farming ground source | |
| Subbotin et al. | Factor analysis of environmental impact of manure utilization | |
| CN110782112A (en) | Method and system for estimating greenhouse gas emission reduction potential in crop production | |
| Bulatewicz et al. | Calibration of a crop model to irrigated water use using a genetic algorithm |
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 | ||
| RJ01 | Rejection of invention patent application after publication | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20190607 |