CN108563863B - Energy consumption calculation and scheduling method for urban rail transit system - Google Patents

Energy consumption calculation and scheduling method for urban rail transit system Download PDF

Info

Publication number
CN108563863B
CN108563863B CN201810321041.9A CN201810321041A CN108563863B CN 108563863 B CN108563863 B CN 108563863B CN 201810321041 A CN201810321041 A CN 201810321041A CN 108563863 B CN108563863 B CN 108563863B
Authority
CN
China
Prior art keywords
attribute
node
urban rail
rail transit
transit system
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201810321041.9A
Other languages
Chinese (zh)
Other versions
CN108563863A (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.)
Beijing Jiaotong University
CRRC Changchun Railway Vehicles Co Ltd
Original Assignee
Beijing Jiaotong University
CRRC Changchun Railway Vehicles Co Ltd
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 Beijing Jiaotong University, CRRC Changchun Railway Vehicles Co Ltd filed Critical Beijing Jiaotong University
Priority to CN201810321041.9A priority Critical patent/CN108563863B/en
Publication of CN108563863A publication Critical patent/CN108563863A/en
Application granted granted Critical
Publication of CN108563863B publication Critical patent/CN108563863B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)
  • Feedback Control In General (AREA)

Abstract

The invention provides an energy consumption calculating and scheduling method of an urban rail transit system. The method comprises the following steps: the method comprises the following steps: constructing a multi-layer network model of an urban rail transit energy consumption system based on the composition of the urban rail transit system; calculating the topological attribute and the data attribute of the node as input quantities, wherein the topological attribute comprises the degree, the betweenness and the clustering coefficient of the node, and the data attribute comprises the energy consumption and the node strength of the node; fusing the topological attribute and the data attribute of the node by using an OWA operator to obtain the weight of the node in the urban rail transit system; based on an emergence theory, the energy efficiency emergence state of the urban rail system, namely the energy efficiency characterization parameter, is calculated by using Choquet integral combined with the weight of the node; and regulating and controlling the energy efficiency of the urban rail operation system based on the urban rail system energy efficiency characterization parameters. For the urban rail transit system with a complex structure, the method can calculate the energy efficiency of the urban rail system, and can reasonably regulate and control the energy efficiency of the urban rail operation system by combining the calculated energy efficiency.

Description

Energy consumption calculation and scheduling method for urban rail transit system
Technical Field
The invention relates to the technical field of rail transit energy consumption calculation, in particular to an energy consumption calculation and scheduling method of an urban rail transit system.
Background
With the development of urban rail transit, the huge energy consumption problem of an urban rail transit system also gets extensive attention. In the existing research on the energy efficiency of the urban rail transit system, no learner utilizes a network model to research the urban rail transit energy efficiency from the perspective of the system, and a strategy and a basis are provided for improving the overall energy efficiency of the system.
Network science is a new cross-science that is specialized in studying qualitative and quantitative laws of complex systems in nature and society, and has been applied to various scientific and engineering fields. Through long-term development, network science has formed the subject architecture and theoretical system of the system. The focus of research on complex network models has gradually shifted in recent years from single networks to multi-layer networks. The multilayer network model is developed as a single-layer network model, and can be used for modeling and analyzing a network which cannot be described by the single-layer network and comprises heterogeneous nodes and heterogeneous edges, so that the multilayer network model is used as a relatively new research branch of network science, is also paid attention to and applied by students, and has very important significance for modeling the energy efficiency of the urban rail transit system by using the multilayer network model.
Disclosure of Invention
The embodiment of the invention provides an energy consumption calculating and scheduling method of an urban rail transit system, which aims to overcome the defects of the prior art.
An energy consumption calculation and scheduling method of an urban rail transit system comprises the following steps:
the method comprises the following steps that firstly, an energy consumption network model of the urban rail transit system is constructed according to the composition relation of the urban rail transit system and by combining energy consumption data;
step two, calculating topological attributes and data attributes of nodes in the urban rail transit system by analyzing operation energy consumption data of the urban rail transit system and combining the constructed energy efficiency network model of the urban rail transit system;
fusing the topological attribute and the data attribute of the node by using an OWA operator to obtain the importance of the node in the urban rail transit system;
calculating energy efficiency characterization parameters of the urban rail transit system by using the importance of the Choquet integral combined with the node based on the emerging theory;
and fifthly, regulating and controlling the energy efficiency of the urban rail transit system based on the energy efficiency characterization parameters of the urban rail transit system.
Further, the constructing an energy consumption network model of the urban rail transit system according to the composition relationship of the urban rail transit system and by combining energy consumption data includes:
1.1, constructing a layer in an energy consumption network model of the urban rail transit system, wherein the layer refers to a set of subsystems with the same structural level or functional level and interaction relations between the subsystems, all layers in the urban rail transit system are divided according to a division principle from an upper layer to a lower layer, and L is a set of layers in the energy consumption network model of the urban rail transit system;
1.2, abstracting a subsystem which is composed of a certain amount of parts in a specific relation and has a relatively independent function and energy consumption in the urban rail transit system as a node, wherein V is a set of nodes in an energy consumption network model of the urban rail transit system;
1.3, abstracting an interaction relation between subsystems and subsystem energy consumption in the urban rail transit system as a connecting edge, taking a Pearson coefficient between the subsystems and subsystem energy consumption data as the weight of the edge, wherein E is a set of edges in an energy consumption network model of the urban rail transit system, and a calculation formula of the weight w of the edge is as follows:
Figure BDA0001625275900000031
in the formula:
n is the sample size;
ai,bi-statistics of two variables respectively.
Further, the calculating of the topology attribute and the data attribute of the nodes in the urban rail transit system by analyzing the operation energy consumption data of the urban rail transit system and combining the constructed energy efficiency network model of the urban rail transit system includes:
2.1, the calculation formula of the centrality of the degree of the node i in the topological attribute is as follows:
Figure BDA0001625275900000032
in the formula:
Figure BDA0001625275900000033
node vLi,jCentrality of degrees of (d);
Figure BDA0001625275900000034
and node vLi,jThe number of directly connected edges;
l V l-the number of nodes in the node set V;
2.2, the calculation formula of the betweenness of the nodes i in the topological attribute is as follows:
Figure BDA0001625275900000041
in the formula:
Figure BDA0001625275900000042
-passing node between any pair of nodes
Figure BDA0001625275900000043
The shortest path number of (2);
JS-the sum of all shortest path numbers between node pairs;
2.3, the calculation formula of the clustering coefficient of the node i in the topological attribute is as follows:
Figure BDA0001625275900000044
in the formula:
Figure BDA0001625275900000045
-a node
Figure BDA0001625275900000046
The number of edges existing between the node and the neighbor node;
Figure BDA0001625275900000047
-a node
Figure BDA0001625275900000048
The total possible number of edges between the adjacent nodes;
2.4, acquiring the energy consumption of the node i in the data attribute from the energy consumption data;
2.5, the calculation formula of the node strength of the node i in the data attribute is as follows:
Figure BDA0001625275900000049
in the formula:
Figure BDA00016252759000000410
-by node
Figure BDA00016252759000000411
The weight of an edge that is a starting or ending point.
Further, the fusion of the topological attribute and the data attribute of the node by using the OWA operator to obtain the importance of the node in the urban rail transit system includes:
3.1, collecting attribute data of each node obtained by investigation and calculation, and constructing a node evaluation attribute matrix;
3.2, standardizing the attribute data in the node evaluation attribute matrix according to an attribute data standardization formula;
3.3, sorting the attribute data according to the size of the attribute data value of the normalized node;
and 3.4, clustering the attribute data values of the nodes by using an OWA operator to obtain the importance value of each node.
Further, the normalizing the attribute data in the node evaluation attribute matrix according to the attribute data normalization formula includes:
setting attribute types of attribute data in the node evaluation attribute matrix, wherein the attribute types comprise benefit type, cost type, fixed type, deviation type, interval type and deviation interval type, and the benefit type attribute means that the larger the value of the attribute data is, the better the attribute data is; the cost attribute means that the smaller the attribute data value, the better; fixed-type attributes are attributes whose data values approach a fixed value alphajThe better; the offset attribute means that the more the attribute data value is offset from a fixed value betajThe better; the interval type attribute means that the closer the attribute data value is to a fixed interval
Figure BDA0001625275900000051
The better; the off-interval type attribute means that the more the attribute data value is off the fixed interval
Figure BDA0001625275900000052
The better;
the attribute data is normalized as follows:
if the attribute value is benefit type, then order
Figure BDA0001625275900000053
In the formula:
maxaj-the largest attribute value in the jth column of attribute data;
if the attribute value is cost type, then order
Figure BDA0001625275900000054
In the formula:
minaj-the smallest attribute value in the jth column of attribute data;
if the attribute value is fixed, order
Figure BDA0001625275900000061
If the attribute value is of offset type, then order
Figure BDA0001625275900000062
If the attribute value is interval type, then order
Figure BDA0001625275900000063
If the attribute value is of the deviated interval type, then order
Figure BDA0001625275900000064
Further, the using the OWA operator to aggregate the attribute data values of the nodes to obtain the importance value of each node includes:
clustering the attribute data values of the nodes by using an OWA operator to obtain the importance value of each node
Figure BDA0001625275900000065
Importance value
Figure BDA0001625275900000066
The calculation formula of (a) is as follows:
Figure BDA0001625275900000067
wherein brijThe br is a sorting result for sorting the normalized attributes of the nodes according to the sequence of the attribute values from big to smallijWhere i refers to the number of nodes, j refers to the number of attributes, m refers to the total number of attribute data used for clustering, and γ is the addition of the OWA operatorAnd the weight vector is determined by using a minimum variance method to the weight vector of the OWA operator, and the calculation formula of gamma is as follows:
Figure BDA0001625275900000071
Figure BDA0001625275900000072
where Disp (γ) represents the degree to which the attribute is used equally, and the order (γ) represents the degree to which the fusion operation is similar to the or operation.
Further, based on the emerging theory, the energy efficiency characterization parameters of the urban rail transit system are calculated by using the Choquet integral in combination with the importance of the nodes, and the method comprises the following steps:
normalizing the calculated importance value of the node, and taking the normalized result of the importance value of the node as a Shapley value of the node, wherein the interval of the Shapley value is 0-1;
calculating fuzzy measure of each node set by establishing an optimization model taking Marchal entropy as a target function, wherein the optimization model taking Marchal entropy as the target function is as follows:
Figure BDA0001625275900000073
Figure BDA0001625275900000074
in the formula
Figure BDA0001625275900000075
Figure BDA0001625275900000076
The | S | is the potential of the attribute set S;
using Choquet integrationCalculating the energy efficiency emergence state of the node, wherein the Choquet integral is defined as: let gλFor a lambda blur measure defined on P (X), f is a non-negative real valued measurable function defined on X, then f is with respect to gλIntegral (c) of (i) integral of (i) fdgλThe definition is shown as the following formula:
Figure BDA0001625275900000081
in the formula:
i——f(xi) Transformation of the vector so that 0 ≦ f (x)1)≤L≤f(xn);
Xi=(x1,x2,L,xn) And f (x)0)=0。
Further, the energy efficiency characterization parameter based on the urban rail transit system regulates and controls the energy efficiency of the urban rail transit system, and includes:
by utilizing the property of Choquet integral, seeking a functional relation between the energy efficiency characterization parameters of the subsystems and the energy efficiency characterization parameters of the urban rail transit system, and determining the energy efficiency importance of each subsystem according to the functional relation;
according to the formula of Choquet integral, when gλ(Ai) Under the known condition, the energy efficiency characterization parameter ^ fd of the urban rail transit systemAnd the energy efficiency characterization parameter f (x) of the subsystemi) The functional relationship between the two is shown as the following formula:
Figure BDA0001625275900000082
determining an energy consumption threshold of a subsystem which can obtain the energy consumption data according to the maximum value and the minimum value of the investigated energy consumption data, and determining an energy efficiency characterization parameter threshold of a subsystem which can not obtain the energy consumption data according to the energy efficiency characterization parameter of the subsystem and the energy consumption threshold of the subsystem which can obtain the energy consumption data; and determining an energy efficiency regulation and control method and strategy of the urban rail transit system according to the energy efficiency importance of each subsystem and the energy efficiency characterization parameter threshold of the subsystem, wherein the energy consumption data is unavailable.
According to the technical scheme provided by the embodiment of the invention, the multi-layer network model is utilized to perform modeling analysis on the urban rail transit system, and the emergence theory and the fuzzy integral method are combined to provide the energy efficiency regulation and control method and strategy of the urban rail transit system based on the energy efficiency emergence state. The method can calculate the energy efficiency of the urban rail transit system and reasonably regulate and control the energy efficiency of the urban rail transit system.
Drawings
In order to more clearly illustrate the technical solutions of the embodiments of the present invention, the drawings needed to be used in the description of the embodiments will be briefly introduced below, and it is obvious that the drawings in the following description are only some embodiments of the present invention, and it is obvious for those skilled in the art to obtain other drawings based on these drawings without inventive labor.
Fig. 1 is a schematic diagram illustrating an implementation principle of an energy consumption calculation and scheduling method of an urban rail transit system according to an embodiment of the present invention;
fig. 2 is a processing flow chart of an energy consumption calculating and scheduling method of an urban rail transit system according to an embodiment of the present invention;
fig. 3 is a schematic diagram of an energy consumption multi-layer network model of an urban rail transit system according to an embodiment of the present invention;
fig. 4 is an abstract diagram of an energy efficiency network model of an urban rail transit system according to an embodiment of the present invention.
Detailed Description
For the convenience of understanding the embodiments of the present invention, the following description will be further explained by taking several specific embodiments as examples in conjunction with the drawings, and the embodiments are not to be construed as limiting the embodiments of the present invention.
Example one
The method considers the complex interaction relationship among urban rail transit energy consumption systems, utilizes a multilayer network model to perform modeling analysis on the urban rail transit system, calculates the energy efficiency emergence state of the urban rail transit system by combining an emergence theory and a fuzzy integral method on the basis, and finally provides the energy efficiency regulation and control method and strategy of the urban rail transit system based on the energy efficiency emergence state.
The schematic diagram of the implementation principle of the energy consumption calculation and scheduling method of the urban rail transit system provided by the embodiment of the invention is shown in fig. 1, and the specific processing flow is shown in fig. 2, and the method comprises the following processing steps:
step 1, constructing an energy consumption network model of the urban rail transit system by combining energy consumption data according to the composition relationship of the urban rail transit system, and specifically comprising the following steps:
1.1, constructing a layer in an energy consumption network model of the urban rail transit system, wherein the layer refers to subsystems with the same structural level or functional level and a set of interaction relations between the subsystems. All layers in the urban rail transit system are divided according to a division principle from an upper layer to a lower layer, and L is a set of layers in an energy consumption network model of the urban rail transit system. Fig. 3 is a schematic diagram of an energy consumption multi-layer network model of an urban rail transit system according to an embodiment of the present invention.
1.2, abstracting subsystems which are formed by a certain amount of parts in a specific relation and have relatively independent functions and energy consumption in the urban rail transit system into nodes. And V is a set of nodes in an energy consumption network model of the urban rail transit system.
1.3, abstracting an interaction relation between subsystems and subsystem energy consumption in the urban rail transit system as a connecting edge, taking a Pearson coefficient between the subsystems and subsystem energy consumption data as the weight of the edge, wherein E is a set of edges in an energy consumption network model of the urban rail transit system, and a calculation formula of the weight w of the edge is as follows:
Figure BDA0001625275900000111
in the formula:
n is the sample size;
ai,bi-statistics of two variables respectively.
And secondly, calculating the topological attribute and the data attribute of the node by analyzing the operation energy consumption data of the urban rail transit system and combining the constructed energy consumption network model of the urban rail transit system.
2.1, the calculation formula of the centrality of the degree of the node i in the topological attribute is as follows:
Figure BDA0001625275900000112
in the formula:
Figure BDA0001625275900000113
-a node
Figure BDA0001625275900000114
Centrality of degrees of (d);
Figure BDA0001625275900000121
-and node
Figure BDA0001625275900000122
The number of directly connected edges;
i V I is the number of nodes in the node set V.
2.2, the calculation formula of the betweenness of the nodes i in the topological attribute is as follows:
Figure BDA0001625275900000123
in the formula:
Figure BDA0001625275900000124
-passing node between any pair of nodesDot
Figure BDA0001625275900000125
The shortest path number of (2);
JS-the sum of all shortest path numbers between node pairs.
2.3, the calculation formula of the clustering coefficient of the node i in the topological attribute is as follows:
Figure BDA0001625275900000126
in the formula:
Figure BDA0001625275900000127
-a node
Figure BDA0001625275900000128
The number of edges existing between the node and the neighbor node;
Figure BDA0001625275900000129
-a node
Figure BDA00016252759000001210
The total number of possible edges between its neighboring nodes.
And 2.4, acquiring the energy consumption of the node i in the data attribute from the energy consumption data.
2.5, the calculation formula of the node strength of the node i in the data attribute is as follows:
Figure BDA00016252759000001211
in the formula:
Figure BDA00016252759000001212
-by node
Figure BDA00016252759000001213
The weight of an edge that is a starting or ending point.
And thirdly, fusing the topological attribute and the data attribute of the node by utilizing an OWA (Ordered Weighted arithmetic mean operator) operator to obtain the importance of the node in the urban rail transit system.
And 3.1, constructing an evaluation target attribute matrix. Is provided with
Figure BDA0001625275900000131
And | η | is the number of attribute sets, wherein | η | is the attribute set of the evaluation target. For any element O in the evaluation target OiI belongs to | O |, all available attribute ηjJ is determined by | η |, to obtain oiAbout ηjAttribute value of aij. Measuring all elements in the evaluation target to obtain attribute values of all elements with respect to eta, thereby forming an attribute matrix A of the evaluation target O with respect to an attribute set eta,
Figure BDA0001625275900000132
as shown in table 1.
TABLE 1 Attribute matrix
Figure BDA0001625275900000133
And 3.2, normalizing the attribute data. The attribute types are generally benefit type, cost type, fixed type, partial type, interval type, partial interval type and the like, wherein the benefit type attribute means that the larger the attribute data value is, the better the attribute data value is; the cost attribute means that the smaller the attribute data value, the better; fixed-type attributes are attributes whose data values approach a fixed value alphajThe better; the offset attribute means that the more the attribute data value is offset from a fixed value betajThe better; the interval type attribute means that the closer the attribute data value is to a fixed interval
Figure BDA0001625275900000134
The better; the off-partition type attribute means that the more the attribute data value is off the fixed partitionWorkshop
Figure BDA0001625275900000135
The better. In order to eliminate the influence of different dimensions on the evaluation result, the attribute data may be normalized as follows.
If the attribute value is benefit type, then order
Figure BDA0001625275900000136
In the formula:
maxaj-the largest attribute value in the jth column of attribute data.
If the attribute value is cost type, then order
Figure BDA0001625275900000141
In the formula:
minaj-the smallest attribute value in the jth column of attribute data.
If the attribute value is fixed, order
Figure BDA0001625275900000142
If the attribute value is of offset type, then order
Figure BDA0001625275900000143
If the attribute value is interval type, then order
Figure BDA0001625275900000144
If the attribute value is of the deviated interval type, then order
Figure BDA0001625275900000145
3.3, sorting the attribute data. Will evaluate the target viNormalized attribute arijSorting the attribute values according to the sequence of the attribute values from large to small, wherein the sorting result is brijListing the attributes of all evaluation targets to form a new attribute matrix NA, wherein the calculation formula of the NA is as follows:
Figure BDA0001625275900000146
and 3.4, calculating the comprehensive attribute value of the evaluation target. The integrated attribute value of the evaluation target (i.e., node) is the importance of the evaluation target. Clustering the attribute values of the evaluation targets by using an OWA operator to obtain the importance of each evaluation target
Figure BDA0001625275900000151
Evaluation target oiThe formula for calculating the importance value is as follows:
Figure BDA0001625275900000152
where γ is the weight vector of the OWA operator and m refers to the total number of attribute data used for clustering. For example, a node has 5 attributes, 3 a attributes, and 2B attributes. Therefore, | η | is 5, M is 3 when calculating the A attribute aggregation, and M is 2 when calculating the B attribute aggregation.
The invention chooses to determine the weight vector gamma of the OWA operator by using the Minimum Variance Method (MVM). The calculation formula for determining the weight vector γ of the OWA operator by the minimum variance method is as follows:
Figure BDA0001625275900000153
Figure BDA0001625275900000154
where Disp (γ) represents the degree to which the attribute is used equally, and the order (γ) represents the degree to which the fusion operation is similar to the or operation.
And fourthly, calculating the energy efficiency emergence state (namely the energy efficiency characterization parameters) of the urban rail transit system by using the Choquet integral combined with the weight of the node based on the emergence theory.
And 4.1, obtaining an abstract diagram of the energy efficiency network model of the urban rail transit system by combining the energy efficiency network model, wherein the abstract diagram is shown in figure 4.
4.2, calculating the node energy value available for the energy consumption data.
And calculating a node energy effective value obtained by the energy consumption data according to the energy consumption investigation data of the urban rail transit system by using an energy efficiency calculation formula, wherein the energy efficiency is equal to the rated energy consumption compared with the actual energy consumption.
And 4.3, calculating to obtain a Shapley value of the node according to the node weight.
The interval of the Shapley value is 0-1, and in order to meet the definition, the invention normalizes the calculated node weight, and takes the normalized result of the node weight as the Shapley value of the node.
4.4, calculating fuzzy measure of attribute set of urban rail transit system
And calculating the fuzzy measure of each attribute set, namely each node set, by establishing an optimization model taking Marchal entropy as a target function. An optimization model with Marchal entropy as an objective function is as follows.
Figure BDA0001625275900000161
Figure BDA0001625275900000162
In the formula
Figure BDA0001625275900000163
Figure BDA0001625275900000164
And | S | is the potential of the attribute set S.
And 4.5, calculating the energy efficiency emergence state (namely the energy efficiency characterization parameters) of the urban rail transit system.
And calculating the energy efficiency emergence state of the urban rail transit system by using the Choquet integral. The definition of Choquet integral is: let gλFor a lambda blur measure defined on P (X), f is a non-negative real valued measurable function defined on X, then f is with respect to gλIntegral (c) of (i) integral of (i) fdgλIs defined as the formula:
Figure BDA0001625275900000171
in the formula:
i——f(xi) Transformation of the vector so that 0 ≦ f (x)1)≤L≤f(xn);Xi=(x1,x2,L,xn) And f (x)0)=0。f(xi) The energy efficiency emergence state of the urban rail transit system is represented; gλ(Ai) Is the fuzzy measure for each node set; the Shapley values of the nodes are used to compute a fuzzy measure for each set of nodes using Marichal entropy.
And fifthly, regulating and controlling the energy efficiency of the urban rail transit system based on the energy efficiency characterization parameters of the urban rail transit system.
And 5.1, calculating the node energy efficiency importance. And (3) seeking a functional relation between the subsystem energy efficiency representation parameters and the system energy efficiency representation parameters by using the property of the Choquet integral and adopting reverse thinking, and determining the energy efficiency importance of each subsystem node according to the coefficient.
According to the formula of Choquet integral, when gλ(Ai) It is known that the original equation can be transformed into ^ fd regarding the system energy efficiency emergence stateAnd subsystem energy efficiency emergence state f (x)i) The functional relation between the two is shown in the formula.
Figure BDA0001625275900000181
And 5.2, generating an energy efficiency control strategy of the urban rail transit system. And generating a control strategy based on the subsystem energy efficiency threshold according to the energy efficiency importance of each subsystem node.
Determining an energy consumption (energy efficiency) threshold of a subsystem which can obtain energy consumption data according to the maximum value and the minimum value of the investigated energy consumption data, and determining an energy efficiency characterization parameter threshold of a subsystem which can not obtain the energy consumption data according to an energy efficiency characterization parameter of the subsystem and the energy consumption (energy efficiency) threshold of the subsystem which can obtain the energy consumption data; and determining an energy efficiency regulation and control method and strategy of the urban rail transit system according to the energy efficiency importance of each subsystem and the energy efficiency characterization parameter threshold of the subsystem, wherein the energy consumption data is unavailable.
Example two
According to the above process, the following description is made of the embodiment using the actual data:
step 1: building a node attribute matrix
According to the node attribute classification, an attribute matrix of the node topology attribute is constructed
Figure BDA0001625275900000182
Wherein etaAAnd | is the number of the topological attributes of the node. According to the information of the nodes in the urban rail system energy efficiency network model, the | V | is 25 and the | eta is knownAAnd | is 3. And (3) inputting the obtained topological attribute of the multilayer network model into an attribute matrix, wherein the attribute matrix A of the node topological attribute is shown as a formula (a).
According to the node attribute classification, an attribute matrix of the node data attribute is constructed
Figure BDA0001625275900000192
Wherein etaBAnd | is the number of the node data attributes. According to the information of the nodes in the urban rail system energy efficiency network model, the | V | is 25 and the | eta is knownBAnd | is 2. According to the method for acquiring data attributes of the multi-layer network model provided in 3.3.2, each of the urban rail energy efficiency network models is acquiredAnd (4) inputting the rated energy consumption of the node and the node strength into the attribute matrix, wherein the attribute matrix B of the node data attribute is shown as a formula (B).
Figure BDA0001625275900000191
Step 2: attribute data normalization and attribute data ordering
It can be found that the topological attribute or the data attribute of the node is an attribute belonging to the benefit type, i.e. the larger the attribute value is, the more important the node is. Therefore, the attribute matrix of the node topology attribute and the attribute matrix of the data attribute of the node are normalized by the normalization formula of the benefit type attribute.
And (4) based on the normalized attribute matrix, sequencing the attributes belonging to the same node according to the attribute values to form a new matrix NA shown as a formula (c) and NB shown as a formula (d).
Figure BDA0001625275900000201
And step 3: calculating importance (i.e. weight) of a node
Due to the fact that the actual energy consumption of the urban rail system is unadditive, only L exists in an energy efficiency network model of the urban rail system4There is an interaction between the nodes of the layers, i.e. there is only L4The nodes in the layer possess the property of node strength. Therefore, the node weights between different layers are not comparable, and the weights of the nodes of different layers are evaluated separately herein.
Specific numerical values in the urban rail system energy efficiency network model are substituted into a weighted vector calculation formula of an OWA operator, and a weighted vector gamma of the node topological attribute can be obtained through MATLAB solutionA(0.584,0.333, 0.083). Weighting vector gamma of node data attributesB=(0.75,0.25)。
Clustering attribute values by using an OWA operator to obtain the weight of each node
Figure BDA0001625275900000213
As shown in table 2.
Table 2: weights of network model nodes
Figure BDA0001625275900000211
And 4, step 4: urban rail system energy efficiency network model abstract diagram
The energy efficiency network model of the urban rail system can be abstracted as shown in FIG. 4:
and 5: node energy value calculation
By using an energy efficiency calculation formula and combining rated energy consumption data and actual energy consumption data of the subsystem nodes obtained through investigation, the energy efficiency of the subsystem nodes can be calculated and obtained as shown in table 3.
Table 3: energy efficiency of measurable subsystem of urban rail system energy consumption data
Figure BDA0001625275900000212
Step 6: node sharley value calculation
The sharley value of a node herein will be calculated by the weight of the node. The sharley values for the nodes are shown in table 4.
Table 4: shapley value of node in urban rail system energy efficiency model
Figure BDA0001625275900000221
And 7: node energy efficiency emergence state calculation
The above is all known information for computing system or subsystem energy efficiency. The importance of the attributes and attribute sets obtained by calculating the optimization model constructed by using the maximization of the Marichal entropy as the objective function is shown in Table 5.
Table 5: fuzzy measure of urban rail system attribute set
Figure BDA0001625275900000231
According to an example of a calculation process of energy efficiency inrush states in a small system, Choquet integration is used to integrate the energy efficiency inrush states of nodes from bottom to top, and energy efficiency inrush states of an urban rail system and subsystems are shown in table 6.
Table 6: energy efficiency emergence state of urban rail system
Figure BDA0001625275900000241
According to the method for identifying the key energy efficiency nodes of the urban rail system, the nodes in the energy efficiency network of the urban rail system are identified, and the energy efficiency importance of each layer of nodes is shown in table 7.
Table 7: node energy efficiency importance degree in urban rail system energy efficiency model
Figure BDA0001625275900000242
According to the energy efficiency calculation formula, the energy efficiency threshold results of the subsystems, which can be obtained by calculating the energy consumption data, are shown in table 8.
Table 8: energy consumption data measurable node energy efficiency threshold of urban rail system
Figure BDA0001625275900000243
According to the method for calculating the energy efficiency characterization parameters of the urban rail system, energy efficiency thresholds of systems and subsystems, which are unavailable in energy consumption data, in the urban rail system are calculated and shown in table 9.
Table 9: urban rail system energy efficiency threshold
Figure BDA0001625275900000244
According to the energy efficiency importance of each layer of node of the urban rail system, a calculation formula corresponding to the energy efficiency characterization parameters of the urban rail system and the subsystems is obtained and is shown in table 10.
Table 10: urban rail system energy efficiency representation parameter calculation formula
Figure BDA0001625275900000245
Figure BDA0001625275900000251
The energy efficiency control strategy of the urban rail system and the variation values of the control nodes on the energy efficiency characterization parameters of the urban rail system are shown in table 11. In the table
Figure BDA0001625275900000253
Representing the variation value of the energy efficiency representation parameter of the urban rail system; and delta theta represents the variation value of the subsystem energy efficiency characterization parameter.
TABLE 11 energy efficiency characterization parameter variation chart of urban rail system
Figure BDA0001625275900000252
Figure BDA0001625275900000261
In summary, the embodiment of the invention provides an energy efficiency regulation and control method and strategy of an urban rail transit system based on an energy efficiency emergence state by utilizing a multilayer network model to perform modeling analysis on the urban rail transit system and combining an emergence theory and a fuzzy integral method. The method can calculate the energy efficiency of the urban rail transit system and reasonably regulate and control the energy efficiency of the urban rail transit system. The method can be particularly applied to the following aspects:
1. the system researches the composition of an urban rail system, the composition of an urban rail energy consumption system and the energy consumption characteristics of the urban rail system.
2. An urban rail system energy efficiency network model is built based on the multilayer network model, characteristics of the model are analyzed, and an idea is provided for identification of key energy consumption nodes of the urban rail system.
3. Based on the emerging theory and an urban rail system energy efficiency network model, the energy efficiency emerging process from microcosmic to macroscopic in the urban rail system is researched, and the energy efficiency of the urban rail system is regulated and controlled from the system perspective.
Those of ordinary skill in the art will understand that: the figures are merely schematic representations of one embodiment, and the blocks or flow diagrams in the figures are not necessarily required to practice the present invention.
Those of ordinary skill in the art will understand that: modules in the devices in the embodiments may be distributed in the devices in the embodiments according to the description of the embodiments, or may be located in one or more devices different from the embodiments with corresponding changes. The modules of the above embodiments may be combined into one module, or further split into multiple sub-modules.
It will be understood by those skilled in the art that all or part of the processes of the methods of the embodiments described above can be implemented by a computer program, which can be stored in a computer-readable storage medium, and when executed, can include the processes of the embodiments of the methods described above. The storage medium may be a magnetic disk, an optical disk, a Read-Only Memory (ROM), a Random Access Memory (RAM), or the like.
The above description is only for the preferred embodiment of the present invention, but the scope of the present invention is not limited thereto, and any changes or substitutions that can be easily conceived by those skilled in the art within the technical scope of the present invention are included in the scope of the present invention. Therefore, the protection scope of the present invention shall be subject to the protection scope of the claims.

Claims (7)

1. An energy consumption calculation and scheduling method of an urban rail transit system is characterized by comprising the following steps:
the method comprises the following steps that firstly, an energy consumption network model of the urban rail transit system is constructed according to the composition relation of the urban rail transit system and by combining energy consumption data;
step two, calculating the topological attribute and the data attribute of the nodes in the urban rail transit system by analyzing the operation energy consumption data of the urban rail transit system and combining the constructed energy consumption network model of the urban rail transit system;
fusing the topological attribute and the data attribute of the node by using an OWA operator to obtain the importance of the node in the urban rail transit system;
calculating energy efficiency characterization parameters of the urban rail transit system by using the importance of the Choquet integral combined with the node based on the emerging theory;
fifthly, regulating and controlling the energy efficiency of the urban rail transit system based on the energy efficiency characterization parameters of the urban rail transit system;
the energy consumption network model of the urban rail transit system is constructed by combining energy consumption data according to the composition relationship of the urban rail transit system, and comprises the following steps:
1.1, constructing a layer in an energy consumption network model of the urban rail transit system, wherein the layer refers to a set of subsystems with the same structural level or functional level and interaction relations between the subsystems, all layers in the urban rail transit system are divided according to a division principle from an upper layer to a lower layer, and L is a set of layers in the energy consumption network model of the urban rail transit system;
1.2, abstracting a subsystem which is composed of a certain amount of parts in a specific relation and has a relatively independent function and energy consumption in the urban rail transit system as a node, wherein V is a set of nodes in an energy consumption network model of the urban rail transit system;
1.3, abstracting an interaction relation between subsystems and subsystem energy consumption in the urban rail transit system as a connecting edge, taking a Pearson coefficient between the subsystems and subsystem energy consumption data as the weight of the edge, wherein E is a set of edges in an energy consumption network model of the urban rail transit system, and a calculation formula of the weight w of the edge is as follows:
Figure FDA0002814853980000011
in the formula:
n is the sample size;
ai,bi-statistics of two variables respectively.
2. The method according to claim 1, wherein the calculating of the topological attribute and the data attribute of the node in the urban rail transit system by analyzing the operation energy consumption data of the urban rail transit system and combining the constructed energy consumption network model of the urban rail transit system comprises:
2.1, the calculation formula of the centrality of the degree of the node i in the topological attribute is as follows:
Figure FDA0002814853980000021
in the formula:
Figure FDA0002814853980000022
-a node
Figure FDA0002814853980000023
Centrality of degrees of (d);
Figure FDA0002814853980000024
-and node
Figure FDA0002814853980000025
The number of directly connected edges;
l V l-the number of nodes in the node set V;
2.2, the calculation formula of the betweenness of the nodes i in the topological attribute is as follows:
Figure FDA0002814853980000026
in the formula:
Figure FDA0002814853980000027
-passing node between any pair of nodes
Figure FDA0002814853980000028
The shortest path number of (2);
JS-the sum of all shortest path numbers between node pairs;
2.3, the calculation formula of the clustering coefficient of the node i in the topological attribute is as follows:
Figure FDA0002814853980000029
in the formula:
Figure FDA00028148539800000210
-a node
Figure FDA00028148539800000211
The number of edges existing between the node and the neighbor node;
Figure FDA00028148539800000212
-a node
Figure FDA00028148539800000213
The total possible number of edges between the adjacent nodes;
2.4, acquiring the energy consumption of the node i in the data attribute from the energy consumption data;
2.5, the calculation formula of the node strength of the node i in the data attribute is as follows:
Figure FDA00028148539800000214
in the formula:
Figure FDA00028148539800000215
-by node
Figure FDA00028148539800000216
The weight of an edge that is a starting or ending point.
3. The method according to claim 2, wherein the fusing the topological attribute and the data attribute of the node by using the OWA operator to obtain the importance of the node in the urban rail transit system comprises:
3.1, collecting attribute data of each node obtained by investigation and calculation, and constructing a node evaluation attribute matrix;
3.2, standardizing the attribute data in the node evaluation attribute matrix according to an attribute data standardization formula;
3.3, sorting the attribute data according to the size of the attribute data value of the normalized node;
and 3.4, clustering the attribute data values of the nodes by using an OWA operator to obtain the importance value of each node.
4. The method of claim 3, wherein the normalizing the attribute data in the node evaluation attribute matrix according to the attribute data normalization formula comprises:
setting attribute types of attribute data in the node evaluation attribute matrix, wherein the attribute types comprise benefit type, cost type, fixed type, deviation type, interval type and deviation interval type, and the benefit type attribute means that the larger the value of the attribute data is, the better the attribute data is; the cost attribute means that the smaller the attribute data value, the better; fixed form of genusThe attribute is that the closer the attribute data value is to a certain fixed value alpha, the better; the offset attribute means that the more the attribute data value is offset from a fixed value betajThe better; the interval type attribute means that the closer the attribute data value is to a fixed interval
Figure FDA0002814853980000031
The better; the off-interval type attribute means that the more the attribute data value is off the fixed interval
Figure FDA0002814853980000032
The better;
the attribute data is normalized as follows:
if the attribute value is benefit type, then order
Figure FDA0002814853980000033
In the formula:
max cj-the largest attribute value in the jth column of attribute data;
if the attribute value is cost type, then order
Figure FDA0002814853980000034
In the formula:
min cj-the smallest attribute value in the jth column of attribute data;
if the attribute value is fixed, order
Figure FDA0002814853980000035
If the attribute value is of offset type, then order
Figure FDA0002814853980000036
If the attribute value is interval type, then order
Figure FDA0002814853980000037
If the attribute value is of the deviated interval type, then order
Figure FDA0002814853980000038
5. The method of claim 4, wherein the clustering the attribute data values of the nodes using the OWA operator to obtain the importance value of each node comprises:
clustering the attribute data values of the nodes by using an OWA operator to obtain the importance value of each node
Figure FDA0002814853980000041
Importance value
Figure FDA0002814853980000042
The calculation formula of (a) is as follows:
Figure FDA0002814853980000043
wherein brijThe br is a sorting result for sorting the normalized attributes of the nodes according to the sequence of the attribute values from big to smallijWherein i refers to the number of nodes, j refers to the number of attributes, m refers to the total number of attribute data for aggregation, γ is the weighted vector of the OWA operator, and the weighted vector of the OWA operator is determined by the minimum variance method, and the calculation formula of γ is as follows:
Figure FDA0002814853980000044
Figure FDA0002814853980000045
where Disp (γ) represents the degree to which the attribute is used equally, and the order (γ) represents the degree to which the fusion operation is similar to the or operation.
6. The method according to claim 3, 4 or 5, wherein the calculating of the energy efficiency characterization parameters of the urban rail transit system by using Choquet integral in combination with the importance of the nodes based on the emerging theory comprises:
normalizing the calculated importance value of the node, and taking the normalized result of the importance value of the node as a Shapley value of the node, wherein the interval of the Shapley value is 0-1;
calculating fuzzy measure of each node set by establishing an optimization model taking Marchal entropy as a target function, wherein the optimization model taking Marchal entropy as the target function is as follows:
Figure FDA0002814853980000051
Figure FDA0002814853980000052
in the formula
Figure FDA0002814853980000053
Figure FDA0002814853980000054
The | S | is the potential of the attribute set S;
energy efficiency flooding state for nodes using Choquet integrationLine calculations, the definition of Choquet integral is: let gλFor a lambda blur measure defined on P (X), f is a non-negative real valued measurable function defined on X, then f is with respect to gλIntegral (c) of (i) integral of (i) fdgλThe definition is shown as the following formula:
Figure FDA0002814853980000055
in the formula:
i——f(xi) Transformation of the vector so that 0 ≦ f (x)1)≤…≤f(xn);Xi=(x1,x2,…,xn) And f (x)0)=0,f(xi) Typical is the energy efficiency emergence state of the urban rail transit system, gλ(Ai) Is the fuzzy measure for each node set, and the Shapley values of the nodes are used to calculate the fuzzy measure for each node set using Marichal entropy.
7. The method according to claim 6, wherein the controlling the energy efficiency of the urban rail transit system based on the energy efficiency characterization parameter of the urban rail transit system comprises:
by utilizing the property of Choquet integral, seeking a functional relation between the energy efficiency characterization parameters of the subsystems and the energy efficiency characterization parameters of the urban rail transit system, and determining the energy efficiency importance of each subsystem according to the functional relation;
according to the formula of Choquet integral, when gλ(Ai) Under the known condition, the energy efficiency characterization parameter ^ fd of the urban rail transit systemAnd the energy efficiency characterization parameter f (x) of the subsystemi) The functional relationship between the two is shown as the following formula:
Figure FDA0002814853980000061
determining an energy consumption threshold of a subsystem which can obtain the energy consumption data according to the maximum value and the minimum value of the investigated energy consumption data, and determining an energy efficiency characterization parameter threshold of a subsystem which can not obtain the energy consumption data according to the energy efficiency characterization parameter of the subsystem and the energy consumption threshold of the subsystem which can obtain the energy consumption data; and determining an energy efficiency regulation and control method and strategy of the urban rail transit system according to the energy efficiency importance of each subsystem and the energy efficiency characterization parameter threshold of the subsystem, wherein the energy consumption data is unavailable.
CN201810321041.9A 2018-04-11 2018-04-11 Energy consumption calculation and scheduling method for urban rail transit system Active CN108563863B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810321041.9A CN108563863B (en) 2018-04-11 2018-04-11 Energy consumption calculation and scheduling method for urban rail transit system

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810321041.9A CN108563863B (en) 2018-04-11 2018-04-11 Energy consumption calculation and scheduling method for urban rail transit system

Publications (2)

Publication Number Publication Date
CN108563863A CN108563863A (en) 2018-09-21
CN108563863B true CN108563863B (en) 2021-03-30

Family

ID=63534562

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810321041.9A Active CN108563863B (en) 2018-04-11 2018-04-11 Energy consumption calculation and scheduling method for urban rail transit system

Country Status (1)

Country Link
CN (1) CN108563863B (en)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110544010B (en) * 2019-07-30 2023-04-07 同济大学 Identification method of key elements influencing global efficiency emergence of rail transit system
CN110688722B (en) * 2019-10-17 2023-08-08 深制科技(苏州)有限公司 Automatic generation method of part attribute matrix based on deep learning
CN113450558B (en) * 2020-03-27 2022-11-08 上海仙豆智能机器人有限公司 Method, system and storage medium for identifying network key node
CN112200320B (en) * 2020-12-02 2021-03-02 成都数联铭品科技有限公司 Model interpretation method, system, equipment and storage medium based on cooperative game method
CN113361792B (en) * 2021-06-21 2022-10-21 吉林大学 Urban electric bus travel energy consumption estimation method based on multivariate nonlinear regression

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103390111A (en) * 2013-07-29 2013-11-13 北京交通大学 Calculation method oriented to single-train running energy consumption of urban rail transit
CN106169103A (en) * 2016-06-24 2016-11-30 北京市地铁运营有限公司地铁运营技术研发中心 A kind of urban track traffic existing line energy consumption Calculating model
CN107203834A (en) * 2017-04-26 2017-09-26 广西电网有限责任公司电力科学研究院 A kind of electricity consumption energy consumption index appraisal procedure towards track traffic

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103699785A (en) * 2013-12-13 2014-04-02 北京交通大学 Urban rail traffic system safety identifying method based on emergence and entropy
CN105809251B (en) * 2016-03-01 2018-01-30 北京交通大学 A kind of EMUs system core part discrimination method of more attribute fusions

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103390111A (en) * 2013-07-29 2013-11-13 北京交通大学 Calculation method oriented to single-train running energy consumption of urban rail transit
CN106169103A (en) * 2016-06-24 2016-11-30 北京市地铁运营有限公司地铁运营技术研发中心 A kind of urban track traffic existing line energy consumption Calculating model
CN107203834A (en) * 2017-04-26 2017-09-26 广西电网有限责任公司电力科学研究院 A kind of electricity consumption energy consumption index appraisal procedure towards track traffic

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
A complex network-based importance measure for mechatronics systems;Yanhui Wang 等;《Physica A:Statistical Mechanics and its Application》;20170131;第466卷(第15期);第180-198页 *

Also Published As

Publication number Publication date
CN108563863A (en) 2018-09-21

Similar Documents

Publication Publication Date Title
CN108563863B (en) Energy consumption calculation and scheduling method for urban rail transit system
Zhang et al. Consensus-based TOPSIS-Sort-B for multi-criteria sorting in the context of group decision-making
CN104951425B (en) A kind of cloud service performance self-adapting type of action system of selection based on deep learning
Xue et al. A survey of application and classification on teaching-learning-based optimization algorithm
Bernal et al. Financial market time series prediction with recurrent neural networks
Ren et al. Group decision making with hesitant fuzzy linguistic preference relations based on modified extent measurement
Moradzadeh et al. Heating and cooling loads forecasting for residential buildings based on hybrid machine learning applications: A comprehensive review and comparative analysis
CN110837915B (en) Low-voltage load point prediction and probability prediction method for power system based on hybrid integrated deep learning
Gautam et al. TOPSIS for multi criteria decision making in intuitionistic fuzzy environment
CN106327357A (en) Load identification method based on improved probabilistic neural network
Nannan et al. Adaptive online time series prediction based on a novel dynamic fuzzy cognitive map
Ye et al. Bridging the micro and macro: Calibration of agent-based model using mean-field dynamics
Li et al. A modular neural network-based population prediction strategy for evolutionary dynamic multi-objective optimization
Handoyo et al. The fuzzy inference system with rule bases generated by using the fuzzy C-means to predict regional minimum wage in Indonesia
Fang et al. Learning decomposed spatial relations for multi-variate time-series modeling
Chen et al. Dynamic multi-objective evolutionary algorithm with center point prediction strategy using ensemble Kalman filter
Singaravel et al. Explainable deep convolutional learning for intuitive model development by non–machine learning domain experts
Orouskhani et al. Fuzzy adaptive cat swarm algorithm and Borda method for solving dynamic multi‐objective problems
Jamil et al. Technological innovation capability evaluation of high-tech firms using conjunctive and disjunctive belief rule-based expert system: a comparative study
CN112465253A (en) Method and device for predicting links in urban road network
Araichi et al. Vine copulas and fuzzy inference to evaluate the solvency capital requirement of multivariate dependent risks
Zhang et al. Prediction study and application of wind power development based on filtering error threshold
Poczeta et al. Structure optimization and learning of fuzzy cognitive map with the use of evolutionary algorithm and graph theory metrics
Li et al. Reduction of the criteria system for identifying effective reservoirs in the joint operation of a flood control system
Rang et al. A Long Short-Term Memory Prediction-Based Dynamic Multi-Objective Evolutionary Optimization Algorithm

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
TA01 Transfer of patent application right
TA01 Transfer of patent application right

Effective date of registration: 20200828

Address after: 100044 Beijing city Haidian District Xizhimen Shangyuan Village No. 3

Applicant after: Beijing Jiaotong University

Applicant after: CRRC CHANGCHUN RAILWAY VEHICLES Co.,Ltd.

Address before: 100044 Beijing city Haidian District Xizhimen Shangyuan Village No. 3

Applicant before: Beijing Jiaotong University

GR01 Patent grant
GR01 Patent grant