Model Property Based Material Balance and Energy Conservation Analysis for Process Industry Energy Transfer Systems

Conventional historical data based material and energy balance analyses are static and isolated computations. Such methods cannot embody the cross-coupling effect of energy flow, material flow and information flow in the process industry; furthermore, they cannot easily realize the effective evaluation and comparison of different energy transfer processes by alternating the model module. In this paper, a novel method for material balance and energy conservation analysis of process industry energy transfer system is developed based on model property. Firstly, a reconfigurable energy transfer process model, which is independent of energy types and energy-consuming equipment, is presented from the viewpoint of the cross-coupling effect of energy flow, material flow and information flow. Thereafter the material balance determination is proposed based on both a dynamic incidence matrix and dynamic balance quantity. Moreover, the model-weighted conservation determination theorem is proved, and the energy efficiency analysis method is also discussed. Results confirmed the efficacy of the proposed methods, confirming its potential for use by process industry in energy efficiency analyses. OPEN ACCESS Energies 2015, 8 12284


Introduction
Model property analysis can be harnessed to reflect the characteristics of an existing system, and furthermore, such an approach can also be adopted in system design to help ensure the correct implementation of a new system [1].An energy transfer process model, as the key sub-model of a process industry energy transfer system, embodies both behaviour and structure characteristics [2].Analyzing the properties of the energy transfer process model can help evaluate and validate any practical system simulated by the model.
A process industry energy transfer system consists of energy-consuming units (for example, equipment) and the coupling relationship between units comprising the system [3].Different equipment exhibit different working characteristics, so the models for equipment are diverse.Research concerning models of energy-consuming equipment is significant and has been predominantly based on the system identification modeling method [4][5][6].However, most of this research only focuses on the energy saving of individual equipment; it neglects the issues of relationship and coordination between the different equipment that comprise an arbitrary system.Nevertheless, the coupling relationship between the different equipment determines the structure of the energy transfer system, and the upstream and downstream interaction; thus a process model reflecting the coupling relationship between the different constituent units is of crucial importance.
At present, common process modeling methods mainly include sequential modular approach, equation-oriented approach and simultaneous modular approach [7][8][9].The sequential modular approach can easily solve the sequential structure process but cannot be well utilized for the simulation of a process with one or more loops because the process partitioning and cutting means must be adopted initially; this will increase the workload and difficulty of the simulation calculation [10].The equation-oriented approach needs to establish the huge equations for a complex process [11].Indeed, computational complexity is the primary drawback for both the equation-oriented approach and the simultaneous modular approach [12].Moreover, the technological process of the process industry is complex and heterogeneous, and the material flow, energy flow and information flow are cross-coupled.How to promote the adaptability of the process model and meanwhile embody the coupling effect of these three different flows is an important issue.Although the various production processes differ with respect to energy type, energy-consuming equipment, materials, products, and so forth, they have in common the process characteristics of parallel, serial, bypass and reflux structures.Therefore, constructing a reconfigurable energy transfer process model, which not only can embody the coupling relationship between energy flow, material flow and information flow, but also can describe the interaction and logic relation between the key factors of the energy transfer process, is meaningful and indispensable for material flow and energy flow analyses of process industry energy transfer systems.
The energy flow, material flow and information flow all coexist in process industry energy transfer systems.Material flow reflects the processing and implementation procedure of the main material product, and energy flow drives the transformation of the material flow.So, the material flow and energy flow, which are the specific embodiments of the production process, determine the resource utilization and energy efficiency.Information flow signals and influences the behaviour of the energy flow and material flow, and also reflects the external environment information.The evolution state of the energy transfer system can be affected or changed by information regulating or external interventions such as the equipment startup/stop.Moeller et al. [13], using Petri nets, researched the simulation of material flow, focusing exclusively on pure material factors such as material input and output, but neglecting the cross-impact of energy, material, information, operator and so forth.Although the material flow and energy flow enter the process or equipment simultaneously, on the output side, there is the material flow along with energy outputs on account of the interaction of energy flow and material flow inside the equipment [14].Only by comprehensively considering the coupling and separation effect between energy and material to research the energy transfer process can the integration and coordination of energy flow and material flow be properly understood and realized.Other researchers have investigated the mass and energy flow models for steel plant based on historical data, and analyzed the influences of enterprise energy structures in different technological processes, for example in [15,16].Here, the researchers analyzed the energy consumption using static matching approaches; however such approaches cannot reflect the dynamic changes of energy.Furthermore, when the industry structure and technical factors are changed, methods based on historical data cannot be easily and effectively applied.As a result of being confined by isolated mass balance and thermal equilibrium, much conventional research on energy flow and material flow is static and insular.Research on process industry material and energy flow must be based on the concepts of energy transfer process networking and dynamic behaviour [2,17].Therefore, research presented in this paper is based on the model of energy transfer process, and the behaviour and structure properties of that model.
So far, taking advantage of graphic description, efficiently computation, as well as the concurrent, asynchronous and resource sharing modeling ability of Petri nets [18], the application of Petri nets and extended Petri nets in manufacture modeling and control is prevalent [19,20].As the process industry is a continuous process along with discrete events, we have explored the fuzzy Petri net with hybrid behavior states and applied it to construct a model of an energy-consuming network including energy and non-energy flow between the different workshops [21,22].This paper, in view of the characteristics of product heterogeneity, uncertainty and continuity of the process industry, firstly constructs the Reconfigurable Process Model of Process Industry Energy Transfer System (RPM_PIETS) from the perspective of comprehensively considering the inherent coupling effects of material, energy and information flows.Based on the semantics of the RPM_PIETS model, the dynamic and structure properties of the model can then be defined, and the material balance and energy conservation analysis of the process industry can be realized based on model properties.
Due to the RPM_PIETS model being developed by improving and extending fuzzy Petri nets in view of the particular modeling demands of process industry energy transfer processes, the analysis methods for RPM_PIETS are different from that of general Petri nets.At present, network analysis methods mainly include the algebraic analysis technique based on state equations, the diagram analysis method based on coverability tree, as well as the inductive analysis method based on decomposition and reduction [23,24].The algebraic analysis technique, firstly proposed by Peterson [25], describes the net system as an incidence matrix and builds the linear relationship with state reachability.Murata [26] showed the properties of Petri nets, especially structure properties such as conservativeness, repetitiveness and so on, based on the relevant results of algebraic analysis.But this method cannot describe the dynamic characteristics adequately.The diagram analysis method with reachability tree reveals the dynamic running mechanism of a system by the finite directed tree, just like the state machine, but the complexity of the states of the complex system is unacceptable [27,28].In order to avoid the analytical difficulty incurred by the state complexity of Petri nets, the inductive analysis method was presented according to the concepts of decomposition and reduction, and crucially, the original properties will still be retained during decomposition and reduction [29][30][31].
From the above discussion, the current model analysis methods have their own characteristics, but there are still certain application limitations and adaptability constraints.As the process industry energy transfer system has the distinct characteristics of continuity, uncertainty, concurrency and asynchronism, RPM_PIETS model is proposed with these characteristics being paramount.When analyzing the material balance and energy conservation of the process industry energy transfer process, this must be based on the semantics of the RPM_PIETS model.Because the normal incidence matrix cannot describe the dynamic features of the model, the respective concepts of static and dynamic incidence matrices are developed.In terms of dynamic and static incidence matrices, the RPM_PIETS model properties are researched, and the material balance analysis is conducted based on dynamic balance quantity.And then, based on an analysis of the model-weighted conservation, the energy conservation calculation of process industry energy transfer system is proposed.
The remainder of this paper is structured as follows: in Section 2, we introduce the formal specification of the reconfigurable process model of process industry energy transfer system after analyzing the key factors and their interactive relations which comprise process industry energy transfer process.In Section 3, we contribute the material balance analysis method based on model dynamic properties; related material balance and material overstock determination are put forward.Next, Section 4 presents the judgment theorem of model-weighted conservation, and describes the energy conservation and energy efficiency analysis for process industry energy transfer process.In Section 5, the validation experiment is presented after which the paper is concluded.

Key Factors and Interactive Relations
Different process industries differ in terms of technological processes, energy-consuming equipment and energy types; nonetheless, there still exists commonalities from the perspective of energy transfer processes.For example, all of the energy transfer processes not only consume the primary energy sources, but also consume the secondary energy sources and the energy-consumed medium, which are all energy carriers; additionally, the energy flow, material flow and information flow all coexist in the energy transfer process, and the enterprise personnel, management and other multiple factors, either individually or in combination, can affect the energy transfer process in the process industry.
To derive the common factors in energy transfer processes, it is necessary to extract the key factors in every aspect of energy transfer processes.The interior factors include energy intensity, energy transmission pipe loss rate, and so forth; the external factors that interact with the energy transfer process include external computation and control, personnel, and so forth.The key factors and interactive relations of the process industry energy transfer process are outlined in Figure 1.In the subsequent section, the reconfigurable process model of process industry energy transfer system is given in terms of Figure 1.The key factors of process industry energy transfer process.

Formal Specification of Reconfigurable Energy Transfer Process Model
Energy transfer behaviour in the process industry is normally continuous; however, discrete events also occur in production processes.The various energy and non-energy materials flow through the equipment with the characteristics of concurrency and asynchronism.Thus, the modeling method for process industry energy transfer process must be process-centric, capable of describing both certain and uncertain information, and support highly effective analyses.In order to realize the behaviour and structure description of an enterprise energy transfer process, a fuzzy Petri nets model oriented to enterprise energy consumption process was explored [21].Here, the Reconfigurable Process Model of Process Industry Energy Transfer System (RPM_PIETS) is further presented in view of energy flow, material flow and information flow.
Element 1: is a finite non-empty set of fuzzy places.Here, P C denotes the set of continuous fuzzy places, which indicates the storehouses of process industry energy and non-energy materials or the warehouses of products.Here, the symbol "◎" represents P C .P D , denoting the set of discrete fuzzy places, indicates the information set such as the message and equipment state.Here, the symbol "○" represents P D .For , the detailed definition for the fuzzy place Pi is shown in is a finite nonempty set of fuzzy transitions, φ P T   and φ P T   .Here, T C denotes the set of continuous fuzzy transitions, which represents the energy-consuming activities or equipment, and T C is symbolized by "□".T D denotes the set of discrete fuzzy transitions, which represents the logic control or message delivery in process industry energy transfer processes such as the start or stop operation for equipment, and T D is symbolized by "▌".
, the fuzzy transition Tj is defined further in Definition 3.
is related to the input-output control of the fuzzy transition, and represents the join or split valves in the energy transmission pipe net.It can control the flowing directions of different input-output material, energy and information flows, and reflect the logic behaviour of concurrence or asynchrony and so forth in the energy transfer process.
Here, it is symbolized by "Δ" in different directions in view of the input or output control gate.Element 4: is a weighted fuzzy relation about T P  .It shows the connectivity from fuzzy places to fuzzy transitions.For , the rated input from Pi to Tj is defined as , λ i j , the fuzzy link strength is defined as , ω i j and the input weight function can be as ) , ( is the set of auxiliary information.Here, E is the external event set, and Q is the set of external computing models and could reflect the working characteristic of the energy-consuming equipment.For example, Q is the set of mathematical models of energy-consuming activities, by which the product or semi-product output can be computed under a given energy input.Π is the set of relevant personnel information.
Definition 2: For , the fuzzy place Pi is defined as ) , ), ( ), ( ( 0 Where, M0(Pi) denotes the initial marking or resource allocation of fuzzy place Pi.
is a mapping that denotes the marking of fuzzy place Pi at time t, which reflects the quantitative variation of the marking of fuzzy place Pi with time.At any time, , where Where, τ : { 0 } indicates the enabling threshold of Tj. ; here j T * denotes the preposition fuzzy place set of Tj. j q is the mathematical or rule model for Tj and Q q j  .It reflects the dynamic characteristic of the energy-consuming activity or equipment.According to j q , the output of the product or semi-finished product of equipment can be computed after the relevant physical and chemical reaction.E e j  is the external event associated with Tj, such as special information release, the occurrence of particular events or other typical external events, and so forth.
(ε,π ) , where ε denotes the operating mode of Tj, 0, automatic mode ε 1, manual mode 2, man-machine interactive mode , and π Π j  associates with the operator of Tj, if ε 0  , π j null  .j G is the logic control gate of Tj and G G j  as mentioned in Definition 1 (where Gj = IGj and Gj = OGj respectively denote the input and output logic control gate of Tj).It fulfills the role of routing management for information flow and energy or non-energy flow, and reflects the logic behaviour of concurrence or asynchrony, etc.It is related to the input-output control of the fuzzy transition Tj, and represents the join or split valves in the energy transmission pipe net.
is the set of input and output fuzzy places of Tj. , where l denotes the hierarchy of Tj, (RPM_PIETSS is the subnet of RPM_PIETS).  .n j v is the rated firing speed which is restricted by the equipment production capacity and the material and energy transmission rate.l j v is the minimum firing speed of Tj.At any time t, the actual firing speed of Tj is , and is the expected value of the triangular fuzzy number ) ( ~t v j , and 0 λ 1   , λ can be valued according to the risk attitude of experts and λ 0.5  means risk neutrality.If comprehensively considering several uncertain factors, ) ( ~t v j can be dynamically computed based on the triangular fuzzy number multi-attribute decision making method [32].

Dynamic Incidence Matrix
An incidence matrix can reflect the net topology of the model in the algebraic analysis technique.For the RPM_PIETS model, due to the existence of the logic control gate as AND along with OR, and the dynamic changes of input and output intensity between energy supplier (fuzzy place) and energy-consuming equipment (fuzzy transition), the logical relation of material flows is different at different time and it changes dynamically.
Definition 4: Let , where , ( ) ( , ) represents the actual input intensity at time t from fuzzy place Pi to fuzzy transition Tj; its value will change according to the enabling condition of fuzzy transition Tj and the marking of fuzzy place Pi.As such, it will change dynamically., ( ) ( , ) represents the actual output intensity at time t from fuzzy transition Tj to fuzzy place Pi, and its value will also change according to the operating mechanism of Tj; here, is defined as the dynamic incidence matrix, and its matrix element is It can be seen from the above Definition 4 that the dynamic incidence matrix reflects the running logical relation between different model nodes at some point, and embodies the flow direction of energy flow and material flow, as well as the logical connections of parallel, series, convergence and divergence and so on.The value of the dynamic incidence matrix element is a positive real number or zero.

Material Balance Determination Based on Dynamic Balance Quantity
In RPM_PIETS model, the marking of the fuzzy place is a function of time and fuzzy transition firing speed, and all the enabled fuzzy transitions will be fired simultaneously.Therefore, the traditional reachability tree method is difficult to be used when analyzing of RPM_PIETS model.Based on the above definition of the dynamic incidence matrix, the dynamic balance quantity which can reflect the quantitative change of energy storehouse or production warehouse is defined below.
Definition 5: In view of the dynamic incidence matrix of RPM_PIETS , for any continuous fuzzy place Pi, let: where is the instantaneous firing speed of the continuous fuzzy transition.For any discrete fuzzy place Pi, let: where σ ( )  , we can find that the dynamic balance quantity can reflect the input and output quantitative balance of the various warehouses of energy, material or production in the whole net system at any time.) (t B i reflects the instantaneous quantitative change of the warehouse Pi, whereby we can get the following inference: If keeping ( ) 0 i B t  for a period of time, it means the material in the warehouse Pi will gradually diminish, representing a shortage in practice.
If keeping ( ) 0 i B t  for a period of time, it means the material in the warehouse Pi will gradually increase, representing an overstock in practice.
If keeping in a steady periodic variation around zero coordinates, the material in the warehouse Pi can be kept in balance, avoiding an inadequate or excessive material supply.

Material Overstock Analysis Based on Boundedness
Definition 6: For , the fuzzy place i P is bounded.
If the maximum capacity, Ki, is defined for the fuzzy place Pi in the model RPM_PIETS, it can ensure the boundedness of Pi during the running of the RPM_PIETS model under the enabling rule.
The model boundedness can help to analyze whether there is production or material overstock under the initial material allocation scheme, and test for bottlenecks in the production process.

Energy Carrier and Energy Conservation
There exists energy and non-energy materials, such as primary energy, secondary energy, energy-consumed medium and diverse raw materials and products and so on.Primary energy refers to the natural and unprocessed energy resources, such as coal, natural gas, and water amongst others.Secondary energy, in order to meet the requirements of production, is transformed directly or indirectly from primary energy, such as electricity, all kinds of oil products, coal gas and so forth.The energy-consumed medium specifies that kind of working material which is not being used as raw material and not entering the product, but which must consume energy when forming itself.Moreover, all kinds of finished products and semi-finished products are made from primary energy, secondary energy and energy-consumed medium.Whatever energy and non-energy materials, there is an equivalent energy value corresponding to each of them.The equivalent energy value of the secondary energy material refers to the amount of the primary energy consumption during producing a unit quantity of the secondary energy.The equivalent energy values of the energy-consumed medium and all kinds of productions also mean the amount of the primary energy consumption during producing a unit quantity of the energy-consumed medium or production.Therefore, each kind of energy can be converted into the equivalent tons of standard coal.
It can be seen that whatever primary energy, secondary energy, energy-consumed medium and diverse raw materials and products, they are all energy carriers and can be transformed or produced by consuming primary energy.So, the analysis for process industry energy transfer system should be executed based on the total energy consumption.For the overall energy transfer system, no matter where and what the energy flows and converts, there exists the characteristic of energy conservation in the whole system.This means all the energy consumed in the production process can be regarded as being transferred to the different productions.So, all of the inputted energy and materials for the equipment can be converted into the supplied energy; all the output products of the equipment can be converted into the effective energy; All the emission waste of the equipment can then be converted into the loss energy.
According to the above idea, the determination method of model-weighted conservation is presented in this section.Through the calculation of energy conservation weighted matrix, the comprehensive energy consumption and the energy efficiency of equipment and production in the process industry are analyzed.In this section, the definitions of the static incidence matrix and the model-weighted conservation are respectively given at first; one kind of model-weighted conservation calculation based on S-invariant is then proposed.

Static Incidence Matrix
The static topology structure of a process industry energy transfer process can be reflected by a static incidence matrix, and the system structure properties depend on the network topology structure.Thus, the definition of the static incidence matrix is also indispensable to the analysis of RPM_PIETS model.
Definition 7: Let denotes whether there exists a connection between fuzzy place Pi and fuzzy transition Tj.If the connection exists, denotes whether there exists a connection between fuzzy transition Tj and fuzzy place Pi.If the connection exists, C  can be referred to as the static incidence matrix of the RPM_PIETS model, and the matrix element . Here,  C is the output static incidence matrix, and  C is the input static incidence matrix.
Based on the static incidence matrix and the literature [29], the S-invariant can be defined as follows: Definition 8: If and only if and any element of vector I is a non-negative real number, this 1  n vector I is the S-invariant of a RPM_PIETS model.

Model-Weighted Conservation Determination
Definition 9: For a RPM_PIETS model, if there exists one |) | ( P n n  dimensional weighted vector , such that for any initial marking set 0 M of

RPM_PIETS and }
, there always exists is the dynamic incidence matrix during the running range dt .σ( ) t , the simultaneous component with the continuous fuzzy transition is zero, reflects the firing situation of every discrete fuzzy transition between the last time dt t  and the current time t.The component of v(u) represents the instantaneous firing speed of the associated continuous fuzzy transition, and the simultaneous component with the discrete fuzzy transition is zero.
From Definitions 4 and 7, it can be known that the static incidence matrix C and the dynamic incidence matrix W(t) are same type of matrix.At any time t, if there exists t X such that ) , we can derive Equation ( 5) in terms of Equation ( 4).
Here, if , and then Thereby, the solvable condition for equation is that the inverse of the square matrix T CC can be found, namely 0  T CC .
Due to I being the S-invariant of RPM_PIETS model and 0 Based on the above derivation and Definition 9 on model-weighted conservation, it can be concluded that the weighted conservation determination of RPM_PIETS model is just as the above Theorem 2.

Energy Efficiency Analysis based on Model-Weighted Conservation
According to Theorem 2, each element value of S-invariant is the weighted value of the relative fuzzy place.Each fuzzy place marking multiplies of the corresponding weighted value and then sums all of them to get the weighted sum which remains unchanged.In other words, for any reachable marking set . For a process industry energy transfer system, the characteristic of model-weighted conservation illustrates the energy conservation of the process industry energy transfer system.At any time, the amounts of all the resources and productions are multiplied by their corresponding weighted values, and their weighted sum is just the total energy value contained in the energy transfer system.
For energy-consuming equipment, multiplying the decreasing marking amounts of each input fuzzy place by its corresponding weighted value to get their weighted sum over a period of time, is just the supplied energy for this equipment.In addition, the increased marking of the fuzzy places related to the production warehouses are multiplied by their corresponding weighted value, and the weighted sum is the effective energy of this equipment.Furthermore, the fuzzy places related to the discharged waste are multiplied by their corresponding weighted value, and the weighted sum is the loss energy of this equipment.In other words, the effective energy of arbitrary equipment can be converted to the production in weighted value, and the loss energy of equipment can be converted to the discharge in weighted value.
Based on above analysis, the total energy contained in the process industry energy transfer system is Assuming that one terminal energy-consuming equipment is related to the fuzzy transition j T , the supplied energy of this equipment is: The effective energy of this equipment is: where e P represents the output production warehouse of the equipment.
The loss energy of the equipment is: where c P represents the discharging waste warehouse of the equipment.Hence, the energy usage ratio M t is the amount of the final product, the comprehensive energy consumption per unit product is just ) (e I , namely the weighted value of the fuzzy place e P .

The Design of the Experiment Platform
The experiment platform was developed using Visual C#.Net, the GDI+, COM, and XML parsing techniques were used also.This software platform, named PIETS_Mod & Sim (Modeling and Simulation Software for Process Industry Energy Transfer System), can support process model construction and simulation for process industry energy transfer systems, and the modeling principle is in accordance with the RPM_PIETS model.The simulation propulsion process was designed as shown in Figure 2.

Example Analysis
Taking a chemical reaction process for example, two kinds of material A and B are, at first, mixed through a mixer, and are then inputted to a reactor which operates in the temperature range of 80~100 °C.The reactor is electrically heated by the model (where, Q represents the electricity consumption, m and c respectively represent the mass and the specific heat capacity of the mixture, and t  represents the temperature change).The least feeding mixture which maintains the reactor working is 4.5L.For this production process, the main energy consumption is electricity except for the material A and B. The transmission loss rates are considered negligible, and there will be m j n j v v  without the limit of pipe transmission rate.Aiming to model property based material balance and energy conservation analysis, the effect of uncertain factors on the equipment firing speed is considered negligible too.The state place of power on M 0 (P 5 ) = 0; d 5 =0; K 5 = 1 The depository of the output reactant M 0 (P 6 ) = 0; d 6 =0; K 6 = +∞ P 7 The temperature change of reactor M 0 (P 7 ) = 80; d 7 = 0; K 6 = 100 P 8 The state place of power off M 0 (P 8 ) = 1; d 8 = 0; K 8 = 1 The discrete control for power on ) 1 ) ( ( ) 80 ) ( ( : The discrete control for power off ) 1 ) ( ( ) 100 ) ( ( : From Figure 3, in the RPM_PIETS model of the chemical reaction process, the discrete fuzzy places 5 P and 8 P can play the signal and control role for electricity supplying, and can ensure that the reactor functions properly in the temperature range 80~100 °C. ) ( 5 t M must be of value 1 if the working temperature of the reactor is below 80 °C, and 0 ) ( if the temperature is over 100 °C; otherwise ) ( 5 t M will be kept.The RPM_PIETS model presented in this paper uses the marking of fuzzy place to embody the quantitative variation of energy and material consumption, and the firing of fuzzy transition to represent the energy consuming activity or equipment.The firing of fuzzy transition depends not only on upstream and downstream material and energy flows, but also on the related discrete fuzzy places with information flow characteristics.According to Definition 7, the static incidence matrix of the RPM_PIETS model of the chemical reaction process can be computed by , where: The known standard coal coefficients of material A, material B and the electricity are ; thus the model-weighted conservation determination approach presented in this paper is validated.Additionally, by the simulation results, the comprehensive energy consumption per unit reactant can be computed: stays above zero longer than it stays below zero; therefore we still need to enhance the operation efficiency of the mixer to satisfy the consecutive material supply for the subsequent reactor.Comparing with the conventional methods, the model property based material balance and energy conservation analysis method developed in this paper has the following advantages:  Some conventional methods for material and energy flow balance analysis are based on the statistical analysis of historical data [33,34].Our model property based method can analyze material and energy flow during the dynamic simulation process of the model, so it is much more intuitive to realize dynamic analysis and find the potential material overstock and scarcity.Furthermore, it can be much easier to compare different technical processes by simply replacing constituent unit modules. Methods only focusing on separate balance analyses of special equipment or units are static matching computations through isolated mass balance and thermal equilibrium; this neglects the complicated correlation amongst constituent units.Our method based on the process model of energy transfer system fully considers the upstream and downstream interaction and the coupling relationship between the different units; thus, the method developed in this paper can help realize the global and dynamic analysis. In the case of known input-output models or model identification module of energy-consuming equipment (as j q related to j T mentioned in Definition 3), following the modeling principle derived in this paper, the process model of energy transfer system can be reconfigured by just resetting the relative production constraints.So, the proposed method can be applied to different industry production processes.

Conclusions
This paper, from the perspective of the cross-coupling effect of energy flow, material flow and information flow, further researches the Reconfigurable Process Model of Process Industry Energy Transfer System (RPM_PIETS).Considering the inherent dynamism and uncertainty of material, energy and information flows, the static and dynamic incidence matrices differentiating the traditional incidence matrix are presented.Thereafter the material balance and energy conservation analyses are primarily formulated.The main contribution of this paper can be summarized as follows:  The Reconfigurable Process Model of Process Industry Energy Transfer System (RPM_PIETS) is presented, which can embody the interaction of the key factors and the complex networking connection of the energy transfer process.It is independent of energy type and energy consuming equipment. Based on the dynamic property of RPM_PIETS, the material balance determination and material overstock analysis methods are proposed. In terms of the concept of energy carrier, the model-weighted conservation determination is advocated and the energy efficiency analysis is also given.
The presented analysis method based on process model property not only benefits the checking of the material balance and energy conservation for the actual system itself, but also improves the reliability design and analysis of process industry energy transfer systems.
In order to enable high integration of both energy flow and material flow, this paper analyzed the energy transfer process by focusing on energy conversion and consumption.The intensity of energy

Figure 1 .
Figure 1.The key factors of process industry energy transfer process.
can be defined differently in view of the actual requirement.For example, when , ω i j represents the maximum transfer rate from Pi to Tj, represents the transmission loss rate from Pi to Tj, consumption in the transmission pipe.Element 5: O(P, T) O(P, T) is a weighted fuzzy relation about T × P. It shows the connectivity from fuzzy transitions to fuzzy places.For i j P P T T     , , the rated output from Tj to Pi is defined as , λ i j , the fuzzy link strength is defined as , ω i j and the output weight function is ) of function O(Pi, Tj) is similar to I(Pi, Tj).Element 6: M0(P) M0(P) denotes the initial marking set of fuzzy places and is the set valued on ] , 0 [  .It can reflect the initial resource allocation in the energy transfer process.Element 7: ∑ For D j T T  , Tj d is the firing delay, which generally describes the firing duration time of the discrete transition.
instantaneous firing speed of Tj is affected by both the material and energy supplies and the uncertain factors, from the viewpoint of uncertain factors, the instantaneous firing speed of Tj is defined as triangular fuzzy number.Here, fuzzy number and denotes the instantaneous firing speed.m j v is the maximum firing speed of Tj under the ideal working condition and

Figure 4 .Figure 5 .
Figure 4.The resource profile of fuzzy place 4 P .
denotes the reserved marking of Pi that can't enable other fuzzy transitions at time t.denotes the non-reserved marking of Pi that can enable the related fuzzy transitions at time t.di is the time limit associated with the fuzzy place Pi.
i K RPM_PIETS model, it specifies the indispensable minimum energy or materials required to enable the equipment to begin working.fj is the control function for the state changing of Tj.If τ Here, P IPO j  is the input fuzzy places.If IGj is the AND gate, * denotes the preposition fuzzy place set of Tj);If IGj is the XOR gate, So we can regard that RPM_PIETS as weighted conservation.
Based on the model derived in this paper, the RPM_PIETS model of this chemical reaction process is illustrated in Figure3.Each element of this model corresponding to Definition 1 is given below.

Table 1 .
 , according to the technical principle of the chemical reaction process and the working mechanism of the equipment, all the model elements and their attribute restrictions corresponding to Definitions 2 and 3 are illustrated in more detail in Table1.Model elements and attribute restrictions.
j T 1.42, 0.81 and 0.123 respectively; the weighted values of both discrete place P5 and P8 which only reflect the state of electricity supplied are 0. In term of Definition 8,

Table 2 .
The simulation results of the example.