CN101276437A - Enterprise energy consumption process model building and emulation method - Google Patents
Enterprise energy consumption process model building and emulation method Download PDFInfo
- 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
Links
Images
Classifications
-
- Y—GENERAL 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
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02P—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
- Y02P80/00—Climate change mitigation technologies for sector-wide applications
- Y02P80/10—Efficient 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
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;
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=(M
0, M
i(t), d
i, h
i), wherein,
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
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
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,
Be T
iThe input object collection,
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,
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,
(a
i=1) ∧ (f
j〉=τ
j), e then
j=1, otherwise, e
j=0; Work as T
jInitial conditions be " or " time,
(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:
2. be continuous when bluring transition
In moment t → t+dt, P
iSign change:
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,
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,
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;
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=(M
0, M
i(t), d
i, h
i), wherein,
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
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
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,
Be T
iThe input object collection,
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,
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,
(a
i=1) ∧ (f
j〉=τ
j), e then
j=1, otherwise, e
j=0; Work as T
jInitial conditions be " or " time,
(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:
2. be continuous when bluring transition
In moment t → t+dt, P
iSign change:
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)
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 |
-
2008
- 2008-05-15 CN CNB200810037456XA patent/CN100570634C/en not_active Expired - Fee Related
Cited By (30)
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 |