CN101276437A - Enterprise energy consumption process model building and emulation method - Google Patents

Enterprise energy consumption process model building and emulation method Download PDF

Info

Publication number
CN101276437A
CN101276437A CNA200810037456XA CN200810037456A CN101276437A CN 101276437 A CN101276437 A CN 101276437A CN A200810037456X A CNA200810037456X A CN A200810037456XA CN 200810037456 A CN200810037456 A CN 200810037456A CN 101276437 A CN101276437 A CN 101276437A
Authority
CN
China
Prior art keywords
fuzzy
energy consumption
transition
storehouse
eec
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.)
Granted
Application number
CNA200810037456XA
Other languages
Chinese (zh)
Other versions
CN100570634C (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.)
Individual
Original Assignee
Individual
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 Individual filed Critical Individual
Priority to CNB200810037456XA priority Critical patent/CN100570634C/en
Publication of CN101276437A publication Critical patent/CN101276437A/en
Application granted granted Critical
Publication of CN100570634C publication Critical patent/CN100570634C/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02PCLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
    • Y02P80/00Climate change mitigation technologies for sector-wide applications
    • Y02P80/10Efficient use of energy, e.g. using compressed air or pressurized fluid as energy carrier

Abstract

The invention particularly relates to a modelling and emulation method for enterprise energy consumption process, belonging to the technical field of energy consumption control, wherein the modelling method comprises using fuzzy variation of a fuzzy Petri net model to represent each energy consumption activity, using fuzzy warehouse to represent the warehouse of raw materials and finished products relevant to the energy consumption activities; the number of marks in the fuzzy warehouse represents the number of the raw materials and the finished products, the relevant supply and demand relationship is denoted by connection line between crunodes; the connection intensity of the connection line signifies the loss rate when conveying the raw materials and finished products from the start point to the end point. The emulation method comprises an initial network model and an excitation fuzzy variation which determines the enabling fuzzy variation set at the current moment, disposes the fuzzy variation at certain excitation moment into the future event table and finds out the minimum excitation time. The method of the invention can clearly and completely describe all kinds of interaction activities in the process of the energy consumption of the enterprise, and the invention has high computing efficiency, strong versatility, which can reinforce the management of energy consumption and save the energy for the enterprise.

Description

Enterprise energy consumption process model building and emulation mode
Technical field
The invention belongs to energy consumption control technology field, be specifically related to a kind of modeling and simulation method of enterprise energy consumption process.
Background technology
Enterprise is the energy resource consumption main body as the intensive terminal user of energy resource consumption.At energy-saving and cost-reducing developing goal, many big-and-middle-sized energy consumption enterprise implements serial of methods and measure.But in actual applications,, the concrete operating position of the energy and the benefit of generation thereof are lacked accurate cognition, influenced the estimation of cost of products and the effective control and the decision-making of enterprise energy consumption because energy consumption enterprise lacks effective efficiency evaluation measures.Therefore, be evaluated as target with enterprise's efficiency, set up the enterprise energy consumption system model, reflect energy consumption situation and influence factor in the enterprise production process comprehensively, exactly, and the energy uses and the dynamic behaviour that consumes in the simulation analysis enterprise production process, will have important Research Significance.It is not only enterprise and realizes that the energy resource consumption situation quantitatively and the important method of qualitative analysis, also is the effective means that enterprise carries out the efficiency assessment, the while also for enterprise implement comprehensively energy-saving and cost-reducing, efficiently use the development strategy of the energy that certain guidance is provided.
Current, both at home and abroad in the development of flowsheeting software, mainly be based on the aspect of technological design, and mostly be the tailored version simulation softward greatly, as, external Aspen Tech company, the petrochemical process simulation softward that SimSci company is released one after another, these softwares are aspect detachment process, and private processes aspect such as olefin polymerization process is ripe relatively, but because modeling technique is still immature, real unitized flowsheeting business software is actually rare, and it is more rare to be based on the software that enterprise's efficiency analyzes.Compare abroad, through years of researches and tachnical storage, China has obtained some achievements in modeling, emulation and the optimization of unit, but with compare the still bigger gap of existence abroad.The main problem that exists is that (1) integration capability is poor.Though each class model is independent to have advanced level technically, for many devices, many processes, multiple operation have serial, the integrated modeling and simulation of whole process of characteristic such as walk abreast, backflow a little less than; (2) model commonality is poor.The model of being developed is often at the specific device of independent enterprise and develop, and causes the model commonality developed poor, and the adaptability of applying is not strong.
At this situation, growth requirement towards high energy consumption industry energy conservation consumption reduction, angle from system, technology and methods such as system ensemble engineering, infotech, the energy and power engineering, environmental engineering, research is based on modeled enterprise energy consumption appraisal procedure, and set up enterprise energy consumption system integration model around this method, with energy consumption process, the material of enterprise move, digitizings such as resource distribution, complementary energy recycling; Based on institute's established model, utilization simulation analysis technology obtains the scientific and rational efficiency evaluation measures of high energy consumption enterprise simultaneously; For the direct energy consumption of business unit's product, indirectly energy consumption and fully enterprise energy consumption situation such as energy consumption forecast analysis is provided, and provide analysis platform for the efficient utilization of enterprise energy most optimum distribution of resources, balanced supply and demand of energy analysis and the energy.
Angle from the enterprise energy consumption method for establishing model, domestic and international currently used modeling method is based on mathematical model, and commonly used have power consumption statistic model, input-output model, linear programming model and system dynamics method and a System Discrimination modeling method etc.These methods have reflected the production scale and the energy consumption situation of enterprise to a certain extent.Especially input-output model is in the enterprise energy consumption cost analysis, the production of the various products that can quantitatively reflect enterprise comprehensively and produced consumes and constitutes, comprising to the consumption of self-made products, to the consumption of outsourcing raw material and fuel and to the consumption of work etc., this model is the good tool of cost accounting and analysis.But input-output model is to be based upon on the basis of historical data analysis to the accounting of energy consumption, when the industrial structure and technical factor change, static input-output model does not then have the ability of energy consumption forecast analysis, input-output model can not reflect process of manufacture intuitively simultaneously, can only reflect quantitatively that production consumption constitutes.Linear programming is an important branch of operational research, research be how some systems keep optimum duty under static state problem, as: the arrangement problem of the production schedule, the selection problem of technological process, produce blending problem, the coordination problem of operation and the assignment problem of Ministry Of Fuel And Power power resource up and down.Though this model analytically has stronger ability optimizing, but and be not suitable for analyzing whole enterprise dynamic, complicated production run.In addition, system dynamics is as the basic skills of socioeconomic systems analysis, can be used for the interaction of internal factor composition and interaction and the system and the external environmental factor of analytic system, and can determine to influence the key factor of system dynamics behavior and restriction relation each other thereof on this basis.This method is suitable for the qualitative analysis between the systematic influence factor having stronger advantage aspect the economic system prediction, and in the changeable energy consumption process of compositing factor, and be not suitable for doing quantitative analysis with this method.The System Discrimination modeling method need not the complicated mechanism of the process of understanding in depth as setting up one of important means of systematic procedure mathematical model, and the system action information that only need utilize the input and output data of system to be provided is set up the mathematical model of systematic procedure.This method is usually used in setting up the model of Business entity energy consumption unit, if but be used for the production of whole, complicated, dynamic enterprise, identification model lacks certain of overall importance and accuracy.Above-mentioned various energy consumption model modeling method, though advantage is respectively arranged, can not be from the overall situation, angle is set up for the enterprise energy consumption process and melted visual, mathematicization in multi-level, the integrated model of one intuitively.In addition, the model of being developed is developed at the specific device of independent enterprise often, and because enterprise's technological process kind is many, causes the model commonality developed poor, and the adaptability of applying is not strong.How fully to inherit on the basis of active cell Research on Equipment Model, To enterprises efficiency assessment set up be independent of specific can Source Type and with energy equipment, and embody that the energy resource consumption process is installed more, the enterprise energy consumption system model of parallel, the serial of many processes, multiple operation, characteristic such as backflow, be one of characteristic innovation of present patent application.
From the angle of flowsheeting method, the process analogy method that current flowsheeting software is adopted mainly contains: sequential modular approach (as Aspen Plus, PRO/II simulation system) and equation solving approach (as HYSYS, SPEEDUP simulation system).Wherein, sequential modular approach is hidden into each model of element equation in each module dispersedly, is that a kind of module of pursuing is carried out Calculation Method.For system architecture just cell by cell rearwards unit transmit, and not only do not had the quick-reading flow sheets situation that logistics but also incompetent stream unit from behind transmit to the unit of front conversely, adopt sequential modular approach to simulate, will be very smooth.But for the flowsheeting problem that has closed circuit, then need could use sequential modular approach by to after flow process piecemeal, cutting and the increase convergence piece, this way is comparatively burdensome, has strengthened the workload of analog computation.Equation solving approach is meant a huge system of equations as whole complicated flow system model of listing, just usually said simultaneous equations, the way of directly finding the solution.This method has been avoided the multi-level iterative computation in the sequential modular approach, so counting yield is higher.But the weak point of equation solving approach is: correctly set up difficulty of huge system of equations; Can not inherit a large amount of unit modules of having developed; Lack large-scale efficiently Nonlinear System of Equations derivation algorithm.
Summary of the invention
The objective of the invention is to propose the modeling and simulation method of the enterprise energy consumption process of a kind of counting yield height, highly versatile.
The enterprise energy consumption process model is the core of enterprise energy consumption system model, at this model, the method that the present invention proposes comprises: " based on the energy consumption process model building method of expansion Fuzzy Petri Net " and " based on the energy consumption process emulation mode of expansion Fuzzy Petri Net ", and particular content is as follows:
(1) based on the energy consumption process model building method of expanding Fuzzy Petri Net
The enterprise energy consumption process model mainly is meant the process logic of enterprise energy consumption and the interaction between the influence factor, comprises all energy consumption activities and the dependence between them of forming energy consumption process.Reflected how energy consumption system compositing factor (energy consumption equipment, energy transmission pipe network etc.) is worked in coordination with and finished Energy Consumption.The foundation of energy consumption process model should be independent of the specific energy and with can equipment, embodies that the energy resource consumption process is installed more, parallel, the serial of many processes, multiple operation, characteristic such as backflow.
At the characteristics of enterprise energy consumption process model, the enterprise energy consumption process model building method should be chosen a kind of processor-oriented modeling method, describes the various interbehaviors of energy consumption process clear, completely; Simultaneously, model should possess certain graphical expressive ability, can reflect production run intuitively, is easy to and user communication.
In the existing process modeling method, fuzzy Bi Teli (Petri) net is as a kind of mathematics and graphical the description and analysis tool of system, not only have and describe parallel, concurrent capacity, and when operation unlike common Petri nets state be discretize, but have successional.Therefore, the present invention selects the ultimate principle based on Fuzzy Petri Net for use, does corresponding expansion according to the characteristics of energy consumption process, thereby proposes " based on the energy consumption process model building method of expansion Fuzzy Petri Net ".
The substance of this method is: represent each energy consumption activity with the fuzzy transition of Fuzzy Petri Nets Model, warehouse with fuzzy storehouse corresponding raw material of represented energy consumption activity and finished product, fuzzy storehouse in reference numerals represent the quantity of raw material and finished product, corresponding relation between supply and demand is represented with the line between node, the proportion of goods damageds the when strength of joint of line is represented from origin-to-destination transportation raw material or finished product.After being familiar with the general energy consumption process of enterprise, using this method enterprise energy consumption process can be mapped as the form of Fuzzy Petri Net, and introduce corresponding transformation rule.
The ultimate principle of method is by providing to give a definition.
Definition 1: the definition of enterprise energy consumption process Fuzzy Petri Nets Model (being designated as EEC-FPN)
Be defined as one 9 tuple: EEC-FPN=(P, T, G, Y I, Y O, SP, EP, M 0(p), ∑).
Wherein, P={P 1, P 2..., P nBe fuzzy storehouse finite aggregate, be used to represent the warehouse of enterprise's raw material (comprising materials such as the energy and the non-energy) and finished product;
T={T 1, T 2..., T mBe fuzzy transition finite aggregate ( Be that P collection and T collection are non-intersect), be used to represent energy consumption activity, energy consumption equipment or the production and processing unit of enterprise;
G={G 1, G 2..., G kBe limited non-Buddhism collection, the control of it and transition interrelates, and can be used for representing branch in the energy transmission pipe network valve that converges and shed, and can control the different material stream that inputs or outputs in the production;
Y IBe Marked Fuzzy relation on P * T, library representation to the connection situation of transition and the rated input i on the connecting line k, strength of joint α kAnd corresponding input intensity computing function Y I(i k, α k), in the energy resource consumption process model, Function Y I(i k, α k) and α kCan adopt different definition as the case may be, for example: as strength of joint α kDuring expression " maximum transfer rate of an energy transmission pipe network " class physical significance, can make Y I(i k, α k)=min{i k, α k, as strength of joint α kDuring expression " loss rate " and so on implication, available Y I(i k, α k)=i k* α kEnergy consumption in the expression energy transmission course;
Y OBe Marked Fuzzy relation on T * P, expression be transitted towards the storehouse the connection situation and the rated output i on the connecting line j, strength of joint β jAnd corresponding mark incremental computations Function Y O(i j, β j), in the energy resource consumption process model, β jAnd Y O(i j, β j) definition, same α kAnd Y I(i k, α k) definition similar;
SP ⋐ P By initial fuzzy storehouse is collected, also be the beginning node of EEC-FPN network. EP ⋐ P For stopping that fuzzy storehouse collects also is the end node of EEC-FPN network.They are respectively applied for beginning and the end position of representing production run;
M 0(P) be defined on the P a value in [0, ∞) in the function of real number, the initial markers state during the operation beginning of library representation place is used to represent initial resource (material) to distribute;
∑=(E, Q Δ ∏) are the additional information set, and wherein E is the external event set; Q is external arithmetic control set, such as: Q can be a System Discrimination algorithm, can identify the mathematical model of certain unit energy consumption equipment in the production through this algorithm, then can calculate under the situation of the specified quantitative input energy energy consumption situation of this equipment and output products amount by this mathematical model; The external tool collection of Δ for supporting to produce; ∏ is personal information set related in the production run.
In the definition 2:EEC-FPN network fuzzy storehouse definition
Right ∀ P i ∈ P , P i=(M 0, M i(t), d i, h i), wherein,
Figure A20081003745600095
M 0Be the fuzzy P of storehouse institute iInitial marking, promptly initial resource is distributed, if P iBe continuous, M then 0Value is arithmetic number and zero, if P iDisperse, then M 0Value is positive integer and zero; M i(t) be the P of storehouse institute iIn the sign of moment t, the variation of interior resource of the fuzzy over time storehouse of its reflection institute or material; d iBe and the fuzzy storehouse P of institute iThe time limit that interrelates.
The definition of fuzzy transition in the definition 3:EEC-FPN network
Right ∀ T i ∈ T , T iBe 10 tuples: T i=(τ i, f i, q i, e i, R i, Δ i, PO i, G i, N i, h Ti, d Ti).
Wherein, τ iBe be defined on the T a value in [0, ∞) in the function of real number, expression transition T iIgniting (enabling) threshold, in the energy consumption process model, produce necessary minimum raw material number or the required energy quantity of energy consumption equipment work in order to be expressed as;
f iBe a map that is defined on the T, it is one to the transition map among the T and is defined in its each input intensity Y I(i k, α k) on a nonnegative function, be called fuzzy transition T iThe state transitions control function, work as f i〉=τ iThe time fuzzy transition have the ability of igniting (enabling);
q i∈ Q is fuzzy transition T iS operation control, can be used for representing the dynamics of energy consumption equipment or unit energy consumption activity passing through q iComputing can obtain the output products that the input energy of specified quantitative produced or the quantity of intermediate product behind certain physics, chemical reaction;
e i∈ E, e iBe and T iThe external event that interrelates, typical external event is as specific message issue, specific incident generation etc.;
R i=(ε, π i), ε represents T iExecutor's information, promptly
Figure A20081003745600102
And π i∈ ∏ ∪ null}, if ε=0, π i=null, otherwise π i∈ ∏;
Δ iThe ∈ Δ represents to carry out the external tool of this task; PO i=IPO i∪ OPO iBe fuzzy transition T iThe inputoutput data object set, IPO i ⋐ P Be T iThe input object collection, OPO i ⋐ P Be T iThe object output collection;
G i∈ G is T iDoor collection (make G i=IG iBe T iThe input gate collection, G i=OG iBe T iThe out gate collection), it with the management functions such as transmission, reception and route of message token or different material stream from fuzzy transition and fuzzy storehouse separate, make structure of models more clear, reasonable;
N i=(l, ∫), l represents this T iNested attribute,
Figure A20081003745600105
If l=1, ∫ ∈ EEC-FPN S, (EEC-FPN SBe the subnet of EEC-FPN);
For discrete T i, d TiBe T iThe time delay that causes, and for continuous T i, then definable blurs transition T iMaximum excitation speed be v Ti, v Ti=1/d Ti
Getting in touch between father's net and the subnet among the definition 4:EEC-FPN
Blur transition T if there is nested in the EEC-FPN model i, i.e. l (T i)=1, ∫ (T i)=EEC-FPN S, and T iThe inputoutput data object set be respectively IPO iAnd OPO i, SP (EEC-FPN then S)=IPO i, EP (EEC-FPN S)=OPO i
(2) based on the energy consumption process emulation mode of expanding Fuzzy Petri Net
Simulation algorithm has partly been described what use is made of institute established model and has been carried out emulation.The enterprise energy consumption process model is based on the expansion Fuzzy Petri Net, so the research of simulation algorithm should launch based on the analytical approach of expansion Fuzzy Petri Nets Model.The evolutionary process of energy consumption system is to advance by the state of discrete event and continuous dynamic process is common, therefore, must consider these two kinds of evolution modes simultaneously in emulation mechanisms.The emulation of energy consumption process is formed by the quantitative change and qualitative change mixed and alternate, in emulation mechanisms, each class variable and status information all by fuzzy storehouse in sign embody, and all kinds of Energy Consumption and dependent event are all finished by fuzzy exciting of transition.The excitation process of fuzzy transition is equal to the consumption of energy consumption activity to the energy, each fuzzy transition excite beginning all will consume its preposition fuzzy storehouse in relevant holder agree (energy), and in these fuzzy transition, produce new Tuo Ken, and this new Tuo Ken is positioned over its rearmounted bluring in the institute of storehouse.According to such thinking, select the step-length of time minimum in each incident (transition excite beginning and transition to excite end all can regard the incident that occurs in the emulation as) as the system emulation clock, realize the propelling of emulation.Concrete simulation algorithm is as follows:
1, initialization pessimistic concurrency control, and global simulation clock t and emulation concluding time tmax, determine each fuzzy storehouse in the model initial marking M0 (P).
2, determine whether current sign is target identification, if emulation finishes, otherwise, change step 3 over to.(perhaps judge the value of current t, if t 〉=t Max, emulation finishes, otherwise changes step 3) over to.
3, determine current t constantly enable fuzzy transition collection T E
3.1, make A={a 1, a 2..., a n} T=0, E={e 1, e 2..., e m} T=0.Here, A represent each fuzzy storehouse the sign situation, E represents the situation that enables of each fuzzy transition, its element is 0 or 1.As the fuzzy P of storehouse institute i(i=1,2 ..., identification number n) greater than 0 or its when being supplied to, a iBe 1, otherwise a iBe 0; As fuzzy transition T j(j=1,2 ..., when m) enabling, e jBe 1, otherwise e jBe 0;
3.2, for the fuzzy P of storehouse institute i, i carries out cycle criterion from 1 to n, if M i(t)>0, a then i=1; Otherwise a i=0;
3.3, for a that satisfies condition arbitrarily i=1 the fuzzy storehouse P of institute iCarry out cycle criterion, whether then to judge the time limit link with it, if, a i=1, otherwise make a i=0;
3.4, for fuzzy transition T j, j carries out cycle criterion from 1 to m, works as T jInitial conditions be " with " time, ∀ P i ∈ T j * , (a i=1) ∧ (f j〉=τ j), e then j=1, otherwise, e j=0; Work as T jInitial conditions be " or " time, ∃ P i ∈ T j * , (a i=1) ∧ (f j〉=τ j), e then j=1, otherwise, e j=0;
3.5, T E={ T j| e jThe fuzzy transition collection of=1} for enable this moment.
4, to enabling fuzzy transition collection T EIn each fuzzy transition T j, determine its excitation instant t respectively according to time delay j, and with t jAnd excite fuzzy transition T accordingly jPut into event table T in the future FIn.
5, the event table T in future of scan model F, find out minimum excitation instant min (t j) and excite fuzzy transition T accordingly j
6, advance global simulation clock t=min (t j), excite corresponding fuzzy transition T j, according to following formula calculate relevant fuzzy storehouse in sign, obtain the sign that next is netted constantly, change step 2 over to.
1. be discrete when bluring transition
Work as T jExcite the back at t+dt:
M i ( t + dt ) = M i ( t ) - I ( P i , T j ) , ∀ P i ∈ IOP j M i ( t ) + O ( P i , T j ) , ∀ P i ∈ OOP j , Here I (P i, T j) and O (P i, T j) be respectively and blur the P of storehouse institute iWith fuzzy transition T jBetween input, output intensity computing function; M i(t) the fuzzy P of storehouse institute of expression moment t iResource quantity.
2. be continuous when bluring transition
In moment t → t+dt, P iSign change:
M i ( t + dt ) = M i ( t ) - v j ( t ) · I ( P i , T j ) · dt , ∀ P i ∈ IOP j M i ( t ) + v j ( t ) · O ( P i , T j ) · dt , ∀ P i ∈ OOP j , Here I (P i, T j) and O (P i, T j) the fuzzy P of storehouse institute of expression equally respectively iWith fuzzy transition T jBetween input, output intensity computing function; M i(t) the fuzzy P of storehouse institute of expression moment t iResource quantity, and v j(t) be fuzzy transition T jExcitation rate.
The inventive method can be described the various interbehaviors of enterprise clear, completely, can react production run intuitively, be easy to and user communication, and, counting yield height, highly versatile, and can instruct enterprise to strengthen managing power consumption, energy savings.
Description of drawings
Fig. 1 is the EEC-FPN model diagram of embodiment.
Among the figure, P 1---the warehouse of unstripped gas A, P 2---the warehouse of unstripped gas B, P 3---by reactor powered can, P 4---the storage of mixer output, P 5---storage power supply state sign, P 6---deposit the input of the required chilled water of heat interchanger, P 7---the storage of reactor output, P 8---deposit the output of the required chilled water of heat interchanger, P 9---the output storage of heat interchanger, P 10---flash evaporator power can, P 11---the storage of products C, P 12---the warehouse of unreacting gas, P 13---inside reactor capacity, P 14---storage power down mode sign, P 15---the storage of waste gas D, T 1---mixer, T 2---abstract type activity, expression reactor working portion, T 3---heat interchanger, T 4---cooling activity, T 5---flash evaporator, T 6---dispenser, promptly complementary energy reclaims link, T 7---powered operation, T 8---the operation that has a power failure, G 1---input and door, the G of mixer 2---input and door, the G of reactor 3---input and door, the G of heat interchanger 4---output and door, the G of heat interchanger 5---input and door, the G of flash evaporator 6---output and door, the G of flash evaporator 7, G 8---the input and the door of confession, power failure logical operation.
Embodiment
With certain reaction under high pressure flow process is example, further specifies the inventive method.
Reaction under high pressure flow process: enter reactor after two kinds of different unstripped gas A, B and circular flow (with A for leading and containing a small amount of B and C) mix.Reactor carries out A+B → C+D reaction under the electrical heating situation, its working temperature should maintain 80 ~ 100 ℃, and when reactor volume reaches 60L, should stop charging.Enter flash evaporator after the cooling of reactor outlet stream stock-traders' know-how heat interchanger.Major product C flows out bottom flash evaporator, and unreacted A (and a spot of B and C) and waste gas D deliver to dispenser after the flash evaporator vapor phase exit is discharged, make waste gas D discharge, and unreacted A (and a spot of B and C) returns use.This production run is except that the consumption to raw material A, B, and the main usefulness of its production equipment can be electricity and water.
The EEC-FPN model of this reaction under high pressure flow process as shown in Figure 1.
In addition, T1 represents mixer and satisfies τ 1=2, f 1 = I ( P 1 , T 1 ) I ( P 2 , T 1 ) , q 1: O (P 3, T 1)=I (P 1, T 1)+I (P 2, T 1), R 1=0, IPO 1={ P 1, P 2, OPO 1={ P 3, l 1=0, h T1=1, V 1=1; T2 is abstract type activity, expression reactor working portion, and satisfy τ 2=0, f 2=(I (P 3, T 2)>0) ∧ (I (P 4, T 2)>0), q 2: O (P 7, T 2)=60, R 2=0, IPO 2={ P 3, P 4, P 5, OPO 2={ P 7, P 13, l 2=1, ∫ 2: slightly, h T2=1, V 2=1/3; M 0(P 1)=49, M 0(P 2)=30, M 0(P 2)=80, M 0(P 14)=1, M 0(P 6)=60, M 0(P 10)=50, other each fuzzy storehouse initial marking be 0.
Provide emulation below and advance 30s, the sampling period is when being 1s, and the number change of the part energy (or material) and equipment state switch instances are as shown in table 1 below:
Number change of each energy of table 1 (or material) and equipment state switch instances
Time Unstripped gas A Unstripped gas B Potpourri Electricity Power supply state Reactant Unreacting gas The reaction output Waste gas Power down mode
Initial value 49 30 0 80 0 0 0 0 0 1
1 47 29 3 80 1 0 0 0 0 0
2 47 29 3 80 1 0 0 0 0 0
3 45 28 2.94 79.99529 1 30.6 0 0 0 0
4 45 28 2.94 79.99529 1 30.6 0 0 0 0
5 43 27 5.94 79.99529 1 30.6 0 0 0 0
6 43 27 5.94 79.99529 1 30.6 0 0 0 0
7 41 26 5.880001 79.99057 1 61.2 0 0 0 0
8 41 26 5.880001 79.99057 0 1.200001 1 60 0 1
9 39.3 25 8.880001 79.99057 1 1.200001 0 60 0.7 0
10 39.3 25 8.880001 79.99057 1 1.200001 0 60 0.7 0
11 37.3 24 8.820001 79.98586 1 31.8 0 60 0.7 0
12 37.3 24 8.820001 79.98586 1 31.8 0 60 0.7 0
13 35.3 23 11.82 79.98586 1 31.8 0 60 0.7 0
14 35.3 23 11.82 79.98586 1 31.8 0 60 0.7 0
15 33.3 22 11.76 79.98114 1 62.4 0 60 0.7 0
16 33.3 22 11.76 79.98114 0 2.400002 1 120 0.7 1
17 31.6 21 14.76 79.98114 1 2.400002 0 120 1.4 0
18 31.6 21 14.76 79.98114 1 2.400002 0 120 1.4 0
19 29.6 20 14.7 79.97643 1 33 0 120 1.4 0
20 29.6 20 14.7 79.97643 1 33 0 120 1.4 0
21 27.6 19 17.7 79.97643 1 33 0 120 1.4 0
22 27.6 19 17.7 79.97643 1 33 0 120 1.4 0
23 25.6 18 17.64 79.97171 1 63.6 0 120 1.4 0
24 25.6 18 17.64 79.97171 0 3.599998 1 180 1.4 1
25 23.9 17 20.64 79.97171 1 3.599998 0 180 2.1 0
26 23.9 17 20.64 79.97171 1 3.599998 0 180 2.1 0
27 21.9 16 20.58 79.967 1 34.2 0 180 2.1 0
28 21.9 16 20.58 79.967 1 34.2 0 180 2.1 0
29 19.9 15 23.58 79.967 1 34.2 0 180 2.1 0
30 19.9 15 23.58 79.967 1 34.2 0 180 2.1 0

Claims (1)

1, a kind of enterprise energy consumption process model building and emulation mode is characterized in that comprising " based on the energy consumption process model building method of expansion Fuzzy Petri Net " and " based on the energy consumption process emulation mode of expansion Fuzzy Petri Net ", and particular content is as follows:
(1) based on the energy consumption process model building method of expanding Fuzzy Petri Net
The substance of this method is: represent each energy consumption activity with the fuzzy transition of enterprise energy consumption process Fuzzy Petri Nets Model, warehouse with fuzzy storehouse corresponding raw material of represented energy consumption activity and finished product, fuzzy storehouse in reference numerals represent the quantity of raw material and finished product, corresponding relation between supply and demand is represented with the line between node, the proportion of goods damageds the when strength of joint of line is represented from origin-to-destination transportation raw material or finished product; Wherein:
Definition 1: the enterprise energy consumption process Fuzzy Petri Nets Model, be designated as EEC-FPN, it is defined as:
One 9 tuple: EEC-FPN=(P, T, G, Y I, Y O, SP, EP, M 0(p), ∑);
Wherein, P={P 1, P 2..., P nBe fuzzy storehouse finite aggregate, be used to represent the warehouse of enterprise's raw material and finished product;
T={T 1, T 2..., T mBe the finite aggregate of fuzzy transition,
Figure A20081003745600021
Be that P collection and T collection are non-intersect, be used to represent energy consumption activity, energy consumption equipment or the production and processing unit of enterprise;
G={G 1, G 2..., G kBe limited non-Buddhism collection, the control of it and transition interrelates, and be used for representing the branch of the energy transmission pipe network valve that converges and shed, and control inputs or outputs the different material stream in the production;
Y IBe Marked Fuzzy relation on P * T, library representation to the connection situation of transition and the rated input i on the connecting line k, strength of joint α kAnd corresponding input intensity computing function Y I(i k, α k), in the energy resource consumption process model, Function Y I(i k, α k) and α kCan adopt different definition as the case may be: as strength of joint α kDuring expression " maximum transfer rate of an energy transmission pipe network " class physical significance, make Y I(i k, α k)=min{i k, α k, as strength of joint α kDuring expression " loss rate " and so on implication, use Y I(i k, α k)=i k* α kEnergy consumption in the expression energy transmission course;
Y OBe Marked Fuzzy relation on T * P, expression be transitted towards the storehouse the connection situation and the rated output i on the connecting line j, strength of joint β jAnd corresponding mark incremental computations Function Y O(i j, β j), in the energy resource consumption process model, β jAnd Y O(i j, β j) definition, same α kAnd Y I(i k, α k) definition similar;
SP ⋐ P For initial fuzzy storehouse collects is the beginning node of EEC-FPN network, EP ⋐ P For stopping that fuzzy storehouse collects is the end node of EEC-FPN network, and they are respectively applied for beginning and the end position of representing production run;
M 0(P) be defined on the P a value in [0, ∞) in the function of real number, the initial markers state during the operation beginning of library representation place is used to represent initial resource to distribute;
∑=(Δ ∏) is the additional information set for E, Q, and wherein E is the external event set, and Q is external arithmetic control set, and the external tool collection of Δ for supporting to produce, ∏ are personal information set related in the production run;
In the definition 2:EEC-FPN network fuzzy storehouse definition
Right ∀ P i ∈ P , P i=(M 0, M i(t), d i, h i), wherein,
Figure A20081003745600032
M 0Be the fuzzy P of storehouse institute iInitial marking, promptly initial resource is distributed, if P iBe continuous, M then 0Value is arithmetic number and zero, if P iDisperse, then M 0Value is positive integer and zero; M i(t) be the P of storehouse institute iIn the sign of moment t, the variation of interior resource of the fuzzy over time storehouse of its reflection institute or material; d iBe and the fuzzy storehouse P of institute iThe time limit that interrelates;
The definition of fuzzy transition in the definition 3:EEC-FPN network
Right ∀ T i ∈ T , T iBe 10 tuples: T i=(τ i, f i, q i, e i, R i, Δ i, PO i, G i, N i, h Ti, d Ti);
Wherein, τ iBe be defined on the T a value in [0, ∞) in the function of real number, expression transition T iThe igniting threshold, in the energy consumption process model, produce necessary minimum raw material number or the required energy quantity of energy consumption equipment work in order to be expressed as;
f iBe a map that is defined on the T, it is one to the transition map among the T and is defined in its each input intensity Y I(i k, α k) on a nonnegative function, be called fuzzy transition T iThe state transitions control function, work as f i〉=τ iThe time fuzzy transition have the ability of igniting;
q i∈ Q is fuzzy transition T iS operation control, be used to represent the dynamics of energy consumption equipment or unit energy consumption activity pass through q iComputing obtain the output products that the input energy of specified quantitative produced or the quantity of intermediate product behind certain physics, chemical reaction;
e i∈ E, e iBe and T iThe external event that interrelates;
R i=(ε, π i), ε represents T iExecutor's information, promptly
Figure A20081003745600034
And π i∈ ∏ ∪ null}, if ε=0, π i=null, otherwise π i∈ ∏;
Δ iThe ∈ Δ represents to carry out the external tool of this task; PO i=IPO i∪ OPO iBe fuzzy transition T iThe inputoutput data object set, IPO i ⋐ P Be T iThe input object collection, OPO i ⋐ P Be T iThe object output collection;
G i∈ G is T iDoor collection, make G i=IG iBe T iThe input gate collection, G i=OG iBe T iThe out gate collection, it with the management function of transmission, reception and the route of message token or different material stream from fuzzy transition and fuzzy storehouse separate, make structure of models more clear, reasonable;
N i=(l, ∫), l represents this T iNested attribute,
Figure A20081003745600037
If l=1, ∫ ∈ EEC-FPN S, EEC-FPN SIt is the subnet of EEC-FPN;
Figure A20081003745600041
For discrete T i, d TiBe T iThe time delay that causes, and for continuous T i, then definable blurs transition T iMaximum excitation speed be v Ti, v Ti=1/d Ti
Getting in touch between father's net and the subnet among the definition 4:EEC-FPN
Blur transition T if there is nested in the EEC-FPN model i, i.e. l (T i)=1, ∫ (T i)=EEC-FPN S, and T iThe inputoutput data object set be respectively IPO iAnd OPO i, SP (EEC-FPN then S)=IPO i, EP (EEC-FPN S)=OPO i
(2) based on the energy consumption process emulation mode of expanding Fuzzy Petri Net
Concrete steps are as follows:
Step 1, initialization pessimistic concurrency control, and global simulation clock t and emulation concluding time tmax, determine each fuzzy storehouse in the model initial marking M0 (P);
Step 2, determine whether current sign is target identification, if emulation finishes, otherwise, change step 3 over to.(perhaps judge the value of current t, if t 〉=t Max, emulation finishes, otherwise changes step 3 over to;
Step 3, determine current t constantly enable fuzzy transition collection T E:
3.1, make A={a 1, a 2..., a n} T=0, E={e 1, e 2..., e m} T=0; Here, A represent each fuzzy storehouse the sign situation, E represents the situation that enables of each fuzzy transition, its element is 0 or 1; As the fuzzy P of storehouse institute i, i=1,2 ..., n, identification number greater than 0 or its when being supplied to, a iBe 1, otherwise a iBe 0; As fuzzy transition T j(j=1,2 ..., when m) enabling, e jBe 1, otherwise e jBe 0;
3.2, for the fuzzy P of storehouse institute i, i carries out cycle criterion from 1 to n, if M i(t)>0, a then i=1; Otherwise a i=0;
3.3, for a that satisfies condition arbitrarily i=1 the fuzzy storehouse P of institute iCarry out cycle criterion, whether then to judge the time limit link with it, if, a i=1, otherwise make a i=0;
3.4, for fuzzy transition T j, j carries out cycle criterion from 1 to m, works as T jInitial conditions be " with " time, ∀ P i ∈ T j * , (a i=1) ∧ (f j〉=τ j), e then j=1, otherwise, e j=0; Work as T jInitial conditions be " or " time, ∃ P i ∈ T j * , (a i=1) ∧ (f j〉=τ j), e then j=1, otherwise, e j=0;
3.5, T E={ T j| e jThe fuzzy transition collection of=1} for enable this moment;
Step 4, to enabling fuzzy transition collection T EIn each fuzzy transition T j, determine its excitation instant t respectively according to time delay j, and with t jAnd excite fuzzy transition T accordingly jPut into event table T in the future FIn;
The event table T in future of step 5, scan model F, find out minimum excitation instant min (t j) and excite fuzzy transition T accordingly j
Step 6, propelling global simulation clock t=min (t j), excite corresponding fuzzy transition T j, according to following formula calculate relevant fuzzy storehouse in sign, obtain the sign that next is netted constantly, change step 2 over to:
1. be discrete when bluring transition
Work as T jExcite the back at t+dt:
M i ( t + dt ) = M i ( t ) - I ( P i , T j ) , ∀ P i ∈ IOP j M i ( t ) + O ( P i , T j ) , ∀ P i ∈ OOP j , Here I (P i, T j) and O (P i, T j) be respectively and blur the P of storehouse institute iWith fuzzy transition T jBetween input, output intensity computing function; M i(t) the fuzzy P of storehouse institute of expression moment t iResource quantity;
2. be continuous when bluring transition
In moment t → t+dt, P iSign change:
M i ( t + dt ) = M i ( t ) - v j ( t ) · I ( P i , T j ) · dt , ∀ P i ∈ IOP j M i ( t ) + v j ( t ) · O ( P i , T j ) · dt , ∀ P i ∈ OOP j , Here I (P i, T j) and O (P i, T j) the fuzzy P of storehouse institute of expression equally respectively iWith fuzzy transition T jBetween input, output intensity computing function; M i(t) the fuzzy P of storehouse institute of expression moment t iResource quantity, and v j(t) be fuzzy transition T jExcitation rate.
CNB200810037456XA 2008-05-15 2008-05-15 Enterprise energy consumption process model building and emulation mode Expired - Fee Related CN100570634C (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CNB200810037456XA CN100570634C (en) 2008-05-15 2008-05-15 Enterprise energy consumption process model building and emulation mode

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CNB200810037456XA CN100570634C (en) 2008-05-15 2008-05-15 Enterprise energy consumption process model building and emulation mode

Publications (2)

Publication Number Publication Date
CN101276437A true CN101276437A (en) 2008-10-01
CN100570634C CN100570634C (en) 2009-12-16

Family

ID=39995860

Family Applications (1)

Application Number Title Priority Date Filing Date
CNB200810037456XA Expired - Fee Related CN100570634C (en) 2008-05-15 2008-05-15 Enterprise energy consumption process model building and emulation mode

Country Status (1)

Country Link
CN (1) CN100570634C (en)

Cited By (19)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102629340A (en) * 2012-03-06 2012-08-08 常州节安得能源科技有限公司 Energy efficiency artificial intelligence analysis system
CN102855354A (en) * 2012-08-21 2013-01-02 南京南瑞继保电气有限公司 Industrial-enterprise-oriented online multi-dimensional energy consumption data statistics modeling method
CN105139161A (en) * 2015-10-14 2015-12-09 常州大学 Virtual enterprise modeling and scheduling method based on Petri network
CN105528470A (en) * 2015-04-29 2016-04-27 南京财经大学 Process industry energy consumption process modeling and simulation method taking uncertain factors into account
CN105700384A (en) * 2016-03-14 2016-06-22 东北大学 Method of building wastewater pollutant discharge Petri net simulation model
WO2016101624A1 (en) * 2014-12-26 2016-06-30 重庆大学 Machining workshop workpiece energy consumption quota formulation method
CN105759615A (en) * 2016-04-06 2016-07-13 浙江工业大学 Fault-tolerant flexible small piece assembly control method based on collaborative Petri network
CN105792270A (en) * 2014-12-24 2016-07-20 国家电网公司 Discrete event simulation method applied to WiMAX (World Interoperability for MicrowaveAccess) system and platform thereof
CN106168769A (en) * 2016-07-19 2016-11-30 同济大学 The modeling of a kind of multiple coupled hybrid flowsheet industrial process and emulation mode
CN108599834A (en) * 2018-03-07 2018-09-28 中国地质大学(武汉) A kind of satellite communication network link utilization analysis method and system
CN110750920A (en) * 2019-09-30 2020-02-04 华中科技大学 Fuzzy Petri network-based manufacturing system energy-saving operation decision method and equipment
CN111722599A (en) * 2020-05-07 2020-09-29 杭州电子科技大学 CPS modeling and analyzing method based on object-oriented generalized stochastic Petri network
CN112036053A (en) * 2020-09-04 2020-12-04 中国航空无线电电子研究所 Petri net-based airport surface sliding simulation and optimization method
CN113378416A (en) * 2021-08-16 2021-09-10 傲林科技有限公司 Event network simulation-based optimization method and device and computer equipment
CN113609697A (en) * 2021-08-16 2021-11-05 傲林科技有限公司 Event network-based analog simulation method and device and computer equipment
CN113610313A (en) * 2021-08-16 2021-11-05 傲林科技有限公司 Event network-based cost reduction optimization method, system and storage medium
CN113610419A (en) * 2021-08-16 2021-11-05 傲林科技有限公司 Event network-based in-business monitoring analysis method, system and storage medium
CN113609780A (en) * 2021-08-16 2021-11-05 傲林科技有限公司 Event network-based clock operation strategy control method and device and electronic equipment
CN113610315A (en) * 2021-08-16 2021-11-05 傲林科技有限公司 Event network-based prediction method and device and electronic equipment

Cited By (30)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102629340B (en) * 2012-03-06 2015-04-08 江苏节安得能源科技有限公司 Energy efficiency artificial intelligence analysis system
CN102629340A (en) * 2012-03-06 2012-08-08 常州节安得能源科技有限公司 Energy efficiency artificial intelligence analysis system
CN102855354A (en) * 2012-08-21 2013-01-02 南京南瑞继保电气有限公司 Industrial-enterprise-oriented online multi-dimensional energy consumption data statistics modeling method
CN102855354B (en) * 2012-08-21 2014-11-05 南京南瑞继保电气有限公司 Industrial-enterprise-oriented online multi-dimensional energy consumption data statistics modeling method
CN105792270A (en) * 2014-12-24 2016-07-20 国家电网公司 Discrete event simulation method applied to WiMAX (World Interoperability for MicrowaveAccess) system and platform thereof
WO2016101624A1 (en) * 2014-12-26 2016-06-30 重庆大学 Machining workshop workpiece energy consumption quota formulation method
CN105528470A (en) * 2015-04-29 2016-04-27 南京财经大学 Process industry energy consumption process modeling and simulation method taking uncertain factors into account
CN105528470B (en) * 2015-04-29 2018-09-18 南京财经大学 The lower flow industrial energy consumption process modeling and simulation method of uncertain factor effect
CN105139161B (en) * 2015-10-14 2018-10-16 常州大学 A kind of Modeling of Virtual Enterprise and dispatching method based on Petri network
CN105139161A (en) * 2015-10-14 2015-12-09 常州大学 Virtual enterprise modeling and scheduling method based on Petri network
CN105700384A (en) * 2016-03-14 2016-06-22 东北大学 Method of building wastewater pollutant discharge Petri net simulation model
CN105759615A (en) * 2016-04-06 2016-07-13 浙江工业大学 Fault-tolerant flexible small piece assembly control method based on collaborative Petri network
CN105759615B (en) * 2016-04-06 2018-09-07 浙江工业大学 A kind of fault tolerant flexibility smallclothes assembly control method based on the Petri network that can cooperate
CN106168769A (en) * 2016-07-19 2016-11-30 同济大学 The modeling of a kind of multiple coupled hybrid flowsheet industrial process and emulation mode
CN108599834A (en) * 2018-03-07 2018-09-28 中国地质大学(武汉) A kind of satellite communication network link utilization analysis method and system
CN110750920A (en) * 2019-09-30 2020-02-04 华中科技大学 Fuzzy Petri network-based manufacturing system energy-saving operation decision method and equipment
CN111722599B (en) * 2020-05-07 2021-10-29 杭州电子科技大学 CPS modeling and analyzing method based on object-oriented generalized stochastic Petri network
CN111722599A (en) * 2020-05-07 2020-09-29 杭州电子科技大学 CPS modeling and analyzing method based on object-oriented generalized stochastic Petri network
CN112036053A (en) * 2020-09-04 2020-12-04 中国航空无线电电子研究所 Petri net-based airport surface sliding simulation and optimization method
CN112036053B (en) * 2020-09-04 2024-05-03 中国航空无线电电子研究所 Airport scene sliding simulation and optimization method based on Petri network
CN113609780A (en) * 2021-08-16 2021-11-05 傲林科技有限公司 Event network-based clock operation strategy control method and device and electronic equipment
CN113610313A (en) * 2021-08-16 2021-11-05 傲林科技有限公司 Event network-based cost reduction optimization method, system and storage medium
CN113610419A (en) * 2021-08-16 2021-11-05 傲林科技有限公司 Event network-based in-business monitoring analysis method, system and storage medium
CN113609697A (en) * 2021-08-16 2021-11-05 傲林科技有限公司 Event network-based analog simulation method and device and computer equipment
CN113610315A (en) * 2021-08-16 2021-11-05 傲林科技有限公司 Event network-based prediction method and device and electronic equipment
CN113610315B (en) * 2021-08-16 2023-10-13 傲林科技有限公司 Event network-based prediction method and device and electronic equipment
CN113610419B (en) * 2021-08-16 2023-10-13 傲林科技有限公司 Event network-based in-event monitoring analysis method, system and storage medium
CN113609780B (en) * 2021-08-16 2023-11-14 傲林科技有限公司 Control method and device for clock running strategy based on event network and electronic equipment
CN113610313B (en) * 2021-08-16 2024-02-06 傲林科技有限公司 Cost reduction optimization method, system and storage medium based on event network
CN113378416A (en) * 2021-08-16 2021-09-10 傲林科技有限公司 Event network simulation-based optimization method and device and computer equipment

Also Published As

Publication number Publication date
CN100570634C (en) 2009-12-16

Similar Documents

Publication Publication Date Title
CN100570634C (en) Enterprise energy consumption process model building and emulation mode
Ma et al. Big data driven predictive production planning for energy-intensive manufacturing industries
Pan et al. Digital twin based real-time production logistics synchronization system in a multi-level computing architecture
Wu et al. Energy efficiency optimization in scheduling crude oil operations of refinery based on linear programming
Gupta et al. Branch and bound experiments in convex nonlinear integer programming
Diwekar Green process design, industrial ecology, and sustainability: A systems analysis perspective
Ho et al. Electricity cost minimisation for optimal makespan solution in flow shop scheduling under time-of-use tariffs
Jalving et al. Graph-based modeling and simulation of complex systems
Ooka et al. A review on optimization techniques for active thermal energy storage control
Guo et al. Reinforcement learning for disassembly system optimization problems: A survey
Zhang et al. A new unrelated parallel machine scheduling problem with tool changes to minimise the total energy consumption
CN106295878A (en) A kind of flexible job shop scheduling system based on Petri network Yu improved adaptive GA-IAGA
Ikeda et al. Application of differential evolution-based constrained optimization methods to district energy optimization and comparison with dynamic programming
Cheng et al. Bi-objective optimization of single-machine batch scheduling under time-of-use electricity prices
Thomas et al. Are intelligent manufacturing systems sustainable?
Zhou et al. Hybrid multi-objective opposite-learning evolutionary algorithm for integrated production and maintenance scheduling with energy consideration
CN101751613A (en) Information processing method for visual chemical-process operational energy efficiency mandatory administration system
Bekker Applying the cross-entropy method in multi-objective optimisation of dynamic stochastic systems
CN102651115B (en) Parallel asynchronous hybrid algorithm processing system and reservoir or Optimal Scheduling of Multi-reservoir System method
CN113848904A (en) Method for optimizing task allocation of multiple mobile robots based on punished energy consumption
Bao et al. A collaborative iterated greedy algorithm with reinforcement learning for energy-aware distributed blocking flow-shop scheduling
CN105528470A (en) Process industry energy consumption process modeling and simulation method taking uncertain factors into account
Zhang et al. A framework of energy-consumption driven discrete manufacturing system
Liu et al. Mixed-flow assembly line balancing with uncertain assembly times in remanufacturing
Zhang et al. Multiobjective particle swarm optimization with direction search and differential evolution for distributed flow-shop scheduling problem

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
C17 Cessation of patent right
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20091216

Termination date: 20120515