CN109473981B - Network charge calculation method considering voltage out-of-limit index under electric power market environment - Google Patents

Network charge calculation method considering voltage out-of-limit index under electric power market environment Download PDF

Info

Publication number
CN109473981B
CN109473981B CN201811465947.4A CN201811465947A CN109473981B CN 109473981 B CN109473981 B CN 109473981B CN 201811465947 A CN201811465947 A CN 201811465947A CN 109473981 B CN109473981 B CN 109473981B
Authority
CN
China
Prior art keywords
bus
reactive
limit
load
voltage
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.)
Expired - Fee Related
Application number
CN201811465947.4A
Other languages
Chinese (zh)
Other versions
CN109473981A (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.)
State Grid Corp of China SGCC
State Grid Hebei Electric Power Co Ltd
Beijing King Star Hi Tech System Control Co Ltd
Original Assignee
State Grid Corp of China SGCC
State Grid Hebei Electric Power Co Ltd
Beijing King Star Hi Tech System Control 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 State Grid Corp of China SGCC, State Grid Hebei Electric Power Co Ltd, Beijing King Star Hi Tech System Control Co Ltd filed Critical State Grid Corp of China SGCC
Priority to CN201811465947.4A priority Critical patent/CN109473981B/en
Publication of CN109473981A publication Critical patent/CN109473981A/en
Application granted granted Critical
Publication of CN109473981B publication Critical patent/CN109473981B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00—Circuit arrangements for ac mains or ac distribution networks
    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
    • H02J2203/20—Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]

Landscapes

  • Engineering & Computer Science (AREA)
  • Power Engineering (AREA)
  • Supply And Distribution Of Alternating Current (AREA)

Abstract

The invention relates to a method for calculating a network charge by considering a voltage out-of-limit index in an electric power market environment, and belongs to the technical field of safety and reliability protection of an electric power system. According to the calculation method, when the voltage is not out of limit, the corresponding bus voltage out-of-limit index is 0, the power grid is indicated to operate within the preset electricity price cost, and when the bus voltage out-of-limit index is larger than 0, the power grid is indicated to exceed the preset cost range, and the operation cost of the power grid is possibly increased due to electric power market transaction. By the calculation method, the bus voltage out-of-limit indexes under the electric power market environment are distributed to market main bodies participating in transaction, such as power generation (power selling), power utilization (power purchasing) and the like, so that certain net-passing fees are charged for the market main bodies causing the voltage out-of-limit risks on the basis of the price of the power transmission and distribution network.

Description

Network charge calculation method considering voltage out-of-limit index under electric power market environment
Technical Field
The invention relates to a method for calculating a network charge by considering a voltage out-of-limit index in an electric power market environment, and belongs to the technical field of safety and reliability protection of an electric power system.
Background
From the situation of the power market developed in China at present, the power market mainly comprises a medium-long term market, a spot market, a frequency modulation auxiliary service market and the like. In the time dimension, the transaction varieties can be divided into two levels of medium-long term (more than day) and spot goods (day ahead and day inside). The medium and long term transaction varieties comprise bilateral negotiation transactions (market electric quantity and base number contracts) which are automatically carried out outside the field and common decomposition curve centralized competition transactions, user-defined decomposition curve listing transactions and base number contract transfer centralized transactions which are carried out in the field through a transaction platform in a unified mode, the transaction results of the varieties are medium and long term price difference contracts, and price difference settlement is carried out according to the day-ahead market prices of the delivery date. The traditional mode of medium-long term contract trading is to stipulate a contract period, contract electric quantity (total electric quantity in the contract period) and trading price, wherein the specific daily trading electric quantity in the contract period is not definitely stipulated and is mainly arranged by a schedule part; in recent electric power market rules (southern (start in Guangdong) electric power spot market rules (survey manuals)), there is a new trading mode of spot medium-long term contracts, which adds the provision of a decomposition curve on the traditional basis, that is, the electric quantity in the contract period needs to be decomposed to the electric quantity at each time of day, and the decomposition mode includes decomposition and publication of a market operation main body and negotiation and decomposition of both market main bodies.
Spot-transaction varieties include day-ahead electric energy markets, real-time electric energy markets and frequency modulation auxiliary service markets. The day-ahead and real-time electric energy market adopts a full-electric-quantity declaration and centralized optimization clearing mode, and forms the most economic (lowest generating cost) market trading result according to the quoted price of the generating side on the premise of meeting the physical constraint of the operation of the electric power system and serves as the basis for executing the scheduling plan. The frequency modulation auxiliary service market determines the calling and compensating mode of frequency modulation resources in a marketization mode, and aims to ensure the real-time running frequency stability of the power system. The result of the market trading is also reflected in the power generation plans of different market entities.
From the above, it can be seen that the analysis of the response of the power market transaction to the operation of the power grid, whether the market transaction is a medium-long term (day or more) market transaction or a spot-market (day ahead, day within), has an influence on the operation of the power grid, and mainly reflects different power generation (power selling) and power utilization (power purchasing) plan curves on different time scales. In medium and long term trading, a daily power trading plan curve is generally formulated by market operation or a scheduling department according to the completion condition of contract electric quantity in the day ahead. This trade curve is superimposed with the original power generation plan curve as part of the day-ahead planning. For the current goods transaction before the day, a market operation department directly designates a current power transaction plan curve before the day generally and adds the current power transaction plan curve into a current power generation plan before the day; for the spot-goods transaction in the day, the market operation department directly gives the power curve of the transaction in the day, and the power grid dispatching department completes the transaction through active rolling dispatching or real-time control in the day. It can be seen that the impact of power market trading on grid operation is reflected in the power generation/utilization plan curves required for market trading generated in different time dimensions (day-ahead, day-in, real-time).
Compared with the traditional power grid in a non-power market environment, the power generation/utilization requirements required by market transactions of different time dimensions in the power market environment mainly comprise the following two aspects;
1) compared with the traditional three-public scheduling, under the power grid scheduling taking the market as the leading factor, the operation mode of the power grid is greatly changed, and the operation mode outside the pre-checking range or the operation mode close to the operation safety boundary of the power grid can occur. For example, after a new energy station participates in medium and long term or spot transaction, when the capacity of the power grid for receiving new energy is limited, the load rate balanced three-party dispatching is not taken as a dispatching target, but the new energy power generation main body participating in a market transaction plan is dispatched preferentially, so that the problems of partial delivery section out-of-limit, insufficient spinning equipment under the system and the like can be caused.
2) Under market-dominated grid scheduling, the uncertainty of grid operation is further enhanced. Due to the trading behavior of the day-ahead spot, day-in spot and ancillary services markets, new power generation, utilization needs may be added in the day-in and real-time dimensions, which may be beyond the scope of traditional day-ahead plan safety checks. Besides the load fluctuation and the intermittent fluctuation of new energy power generation, the fluctuation caused by market transaction is increased, and new requirements are provided for a power grid dispatching operation department.
In the electric power market environment, due to the day-ahead, day-in and spot-shipment trades participated in by the market main body, the traditional planning and dispatching mode of the power grid is influenced, how to quantitatively evaluate the influence is reflected to the corresponding market trading main body, and then the influence and the constraint are carried out on the factors through the price lever, so that the electric power market environment is an important problem of power grid dispatching.
In the current electric power market, the prices of transmission and distribution networks are mainly checked in advance, on the basis of providing assets and costs related to transmission and distribution services, according to the principle of permitting the costs plus reasonable benefits and considering cost factors such as depreciation fees and operation maintenance fees, checking is performed periodically and executed (national committee for modification issues about issuing notices of pricing methods (trial execution) for provincial-level power grid transmission and distribution electricity, issue modification prices (2016) 2711).
When the power grid normally operates and the voltage is not out of limit, the corresponding bus voltage out-of-limit index is 0, which indicates that the power grid operates within the preset electricity price cost, and when the bus voltage out-of-limit index is larger than 0, which indicates that the power grid exceeds the preset cost range, the power market transaction may cause the increase of the power grid operation cost. The invention mainly provides a method for calculating a network-passing fee by considering a voltage out-of-limit index in an electric power market environment, which can distribute the bus voltage out-of-limit index in the electric power market environment to market main bodies participating in trading, such as power generation (power selling), power utilization (power purchasing) and the like, so that extra network-passing fee is charged for trading main bodies causing voltage out-of-limit risks on the basis of the price of a power transmission and distribution network.
Grand bin, zhanberming, and yearly, proposed a quasi-steady-state sensitivity method in the quasi-steady-state sensitivity analysis method (the report of motor engineering in china, V19N 4, 1999, 4 months, pp.9-13), which is different from the conventional static sensitivity analysis method, takes into account the quasi-steady-state physical response of the power system, and takes into account the total change between the new and old steady states before and after the system control, thereby effectively improving the accuracy of the sensitivity analysis. The method is based on a PQ decoupling model of a power system, and when a generator is provided with an Automatic Voltage Regulator (AVR), the generator node can be regarded as a PV node; when the generator is provided with automatic reactive power regulation (AQR) or Automatic Power Factor Regulation (APFR), the generator node is considered to be a PQ node, which is the same as a common load node. In addition, the load voltage static characteristic is considered as a primary or secondary curve of the node voltage. The established power flow model naturally takes these quasi-steady-state physical responses into account, so that the sensitivity calculated based on the power flow model is the quasi-steady-state sensitivity. The quasi-steady state sensitivity method described above was used in the calculations herein.
Disclosure of Invention
The invention aims to provide a method for calculating a net-crossing fee by considering a voltage out-of-limit index in an electric power market environment, which considers the quasi-steady-state sensitivity of the operation of a power grid, namely the active power change of each main transformer in the power grid after unit active power is added to a certain bus, and the bus voltage out-of-limit index in the operation of the power grid in the electric power market environment is distributed to each market main body participating in electric power market transaction.
The invention provides a method for calculating a net charge by considering a voltage out-of-limit index under a power market environment, which comprises the following steps of:
(1) establishing a market main body model, wherein the market main body model mainly refers to main body objects participating in electric power market transaction, and the market main body model comprises a power generation object and a power utilization object:
(1-1) market subject model S of Power Generation targetiIs defined as:
Si＝{Unj,j＝1,..,J,Pun,j＞0}∪{Ldk,k＝1,..,K,Pld,k＜0}，
the power generation object comprises a generator model with positive active output and a load model with negative active load, wherein UnjAnd LdkCalculating generators and loads in the model for the power grid;
(1-2) market subject model B of electric objectiIs defined as:
Bi＝{Unj,j＝1,..,J,Pun,j＜0}∪{Ldk,k＝1,..,K,Pld,k＞0}，
the power utilization object comprises a generator model with negative active output and a load model with positive active load, wherein UnjAnd LdkRespectively calculating a generator and a load in the power grid calculation model;
(2) calculating the out-of-limit index of the bus voltage:
the voltage index mainly evaluates whether the bus voltage is out of limit, and the index calculation method for a single bus comprises the following steps:
Figure GDA0003473316050000041
wherein V is the current value of the bus voltage, Vmax、VminIs the voltage limit value, k is less than or equal to 1, and k is generally 0.9;
(3) the method for apportioning out-of-limit indexes of the bus voltage comprises the following specific steps:
(3-1) for the bus with higher voltage limit, the out-of-limit index is for the power generation main body SiThe apportionment calculation method is as follows:
Figure GDA0003473316050000042
in the above formula (2), omegaBFor a set including the current upper bound buses, each of the upper bound buses is denoted as b, RbThe out-of-limit index of the b-th bus bar,
Figure GDA0003473316050000043
for generating (selling) electricityiComprising a generator UnjThe contribution ratio of the bus voltage to the upper limit is calculated by the following method:
Figure GDA0003473316050000044
wherein b is an out-of-limit bus, SBG(b,uj) Generator Un for reactive voltage sensitivity calculationjReactive voltage sensitivity to the b-th bus, uj being the injected bus-in matrix S of the generator UnjBGColumn number in, Qun,jFor the current reactive power output of the generator, Ω z is the set of buses with reactive power injection greater than 0 in the region, QxThe reactive value of the generator or the load is injected into the bus;
in the formula (2), the reaction mixture is,
Figure GDA0003473316050000045
for the main body S of electricity generation (electricity sale)iContaining load LdkThe contribution ratio of the bus voltage to the upper limit is calculated by the following method:
Figure GDA0003473316050000046
wherein b is an out-of-limit bus, SBG(b,dk) Load Ld given for reactive voltage sensitivity calculationkReactive voltage sensitivity to the b-th bus, dkIs a load LdkIs in the matrix SBGColumn number in, Qld,jFor the current reactive power output of the load, Ω z is the set of busbars with reactive injection greater than 0 in the region, QxThe reactive value of the generator or the load is injected into the bus;
(3-2) for the bus with lower voltage limit, the out-of-limit index is for the main body S of power generation (power selling)iThe allocation method comprises the following steps:
Figure GDA0003473316050000051
in the above formula (5), omegaBFor the set of buses including the current lower limit, each of the lower limit buses is denoted as b, RbIs the out-of-limit index of the b-th bus,
Figure GDA0003473316050000052
for generating (selling) electricityiComprising a generator UnjThe contribution ratio to the lower limit of the bus voltage is calculated by:
Figure GDA0003473316050000053
wherein b is an out-of-limit bus, SBG(b,uj) Reactive voltage sensitivity of the generator Unj to the b-th bus given for reactive voltage sensitivity calculation uj is the injected bus-in matrix S of the generator UnjBGColumn number in, Qun,jFor the current reactive power output of the generator, Ω z is the set of buses with reactive power injection less than 0 in the region, QxReactive value of the generator or load which is absorbing reactive power from the bus;
in the above-mentioned formula (5),
Figure GDA0003473316050000054
for the main body S of electricity generation (electricity sale)iContaining load LdkThe contribution ratio to the lower limit of the bus voltage is calculated by:
Figure GDA0003473316050000055
wherein b is an out-of-limit bus, SBG(b,dk) Reactive voltage sensitivity, d, of load Unj to the b-th bus given for reactive voltage sensitivity calculationkIs a load LdkIs in the matrix SBGColumn number in, Qld,jFor the current reactive power output of the load, Ω z is the set of busbars with reactive injection less than 0 in the region, QxReactive value of the generator or load which is absorbing reactive power from the bus;
(3-3) bus with higher voltage limit, its out-of-limit index is for main body BiThe apportionment calculation method is as follows:
Figure GDA0003473316050000061
in the above formula (8), omegaBThe method comprises the steps that a set comprising current bus bars exceeding the upper limit is marked as b; rbThe out-of-limit index of the b-th bus bar,
Figure GDA0003473316050000062
for electricity consumption (purchasing electricity) main body BiComprising a generator UnjThe contribution ratio of the bus voltage to the upper limit is calculated by the following method:
Figure GDA0003473316050000063
wherein b is an out-of-limit bus, SBG(b,uj) Generator Un for reactive voltage sensitivity calculationjReactive voltage sensitivity to the b-th bus, uj being the injected bus-in matrix S of the generator UnjBGColumn number in, Qun,jFor the current reactive power output of the generator, Ω z is the set of buses with reactive power injection greater than 0 in the region, QxThe reactive value of the generator or the load is injected into the bus;
in the above-mentioned formula (8),
Figure GDA0003473316050000064
for the main body B of electricity consumption (purchasing electricity)iContaining load LdkThe contribution ratio of the bus voltage to the upper limit is calculated by the following method:
Figure GDA0003473316050000065
wherein b is an out-of-limit bus, SBG(b,dk) Load Ld given for reactive voltage sensitivity calculationkReactive voltage sensitivity to the b-th bus, dkIs a load LdkIs in the matrix SBGColumn number in, Qld,jFor the current reactive power output of the load, Ω z is the reactive power in the regionInjecting a bus set, Q, greater than 0xThe reactive value of the generator or the load is injected into the bus;
(3-4) for the bus with lower voltage limit, the out-of-limit index is for the electricity main body BiThe allocation method comprises the following steps:
Figure GDA0003473316050000071
in the above formula (11), omegaBFor the set of buses including the current lower limit, each of the lower limit buses is denoted as b, RbIs the out-of-limit index of the b-th bus,
Figure GDA0003473316050000072
is a power consumption main body BiComprising a generator UnjThe contribution ratio to the lower limit of the bus voltage is calculated by:
Figure GDA0003473316050000073
wherein b is an out-of-limit bus, SBG(b,uj) Reactive voltage sensitivity of the generator Unj to the b-th bus given for reactive voltage sensitivity calculation uj is the injected bus-in matrix S of the generator UnjBGColumn number in, Qun,jFor the current reactive power output of the generator, Ω z is the set of buses with reactive power injection less than 0 in the region, QxReactive value of the generator or load which is absorbing reactive power from the bus;
in the above-mentioned formula (12),
Figure GDA0003473316050000074
for the main body B of electricity consumption (purchasing electricity)iContaining load LdkThe contribution ratio to the lower limit of the bus voltage is calculated by:
Figure GDA0003473316050000075
whereinB is an out-of-limit bus, SBG(b,dk) Reactive voltage sensitivity, d, of load Unj to the b-th bus given for reactive voltage sensitivity calculationkIs a load LdkIs in the matrix SBGColumn number in, Qld,jFor the current reactive power output of the load, Ω z is the set of busbars with reactive injection less than 0 in the region, QxReactive value of the generator or load which is absorbing reactive power from the bus;
(4) correcting the network passing fee of the market subject:
correcting the basic net charge P of the market subject according to the calculated distribution result of the voltage out-of-limit indexiMultiplying the basic net charge by a weight coefficient Fi：
Figure GDA0003473316050000076
Pi＝Pi·Fi (15)
Wherein the content of the first and second substances,
Figure GDA0003473316050000081
the voltage out-of-limit index value allocated to the market body,
Figure GDA0003473316050000082
is the sum of all bus voltage out-of-limit indexes in the current area,
Figure GDA0003473316050000083
Ω z is the set of all bus out-of-limit indicators in the region.
The invention provides a method for calculating the net charge by considering the voltage out-of-limit index under the power market environment, which has the advantages that:
according to the method for calculating the net-crossing fee by considering the voltage out-of-limit index in the power market environment, on one hand, the risk condition that the bus voltage is out-of-limit in the power market environment can be evaluated, and the risk condition that the bus voltage is out-of-limit possibly appears in the current and future power grids can be given; on the other hand, the risk indicators can be distributed to the main bodies (selling electricity and purchasing electricity) participating in market trading, so that the analysis of trading cost and the formulation of electricity price can be guided.
Detailed Description
The invention provides a method for calculating a network charge by considering a voltage out-of-limit index under a power market environment, which comprises the following steps of:
(1) establishing a market main body model, wherein the market main body model mainly refers to main body objects participating in electric power market transaction, and the market main body model comprises a power generation (electricity selling) object and a power utilization (electricity purchasing) object:
(1-1) market subject model S of Power Generation (selling of Power)iIs defined as:
Si＝{Unj,j＝1,..,J,Pun,j＞0}∪{Ldk,k＝1,..,K,Pld,k＜0}，
the power generation (selling) object comprises a generator model with positive active output and a load model with negative active load, wherein UnjAnd LdkCalculating generators and loads in the model for the power grid;
(1-2) market subject model B of electricity consumption (purchase) objectiIs defined as:
Bi＝{Unj,j＝1,..,J,Pun,j＜0}∪{Ldk,k＝1,..,K,Pld,k＞0}，
the electricity consumption (purchase) objects comprise a generator model with negative active output and a load model with positive active load, wherein UnjAnd LdkRespectively calculating a generator and a load in the power grid calculation model;
(2) calculating the out-of-limit index of the bus voltage:
the voltage index mainly evaluates whether the bus voltage is out of limit, and the index calculation method for a single bus comprises the following steps:
Figure GDA0003473316050000091
wherein V is the current value of the bus voltage, Vmax、VminIs the voltage limit value, k is less than or equal to 1, and k is generally 0.9;
(3) the bus voltage out-of-limit index is apportioned, and the sensitivity of active and reactive injection to the bus voltage is mainly adopted for the apportionment of the voltage out-of-limit index. For the bus with the voltage higher than the upper limit, the risk index is allocated to the generator with the reactive output larger than 0 or the load with the reactive load smaller than 0 only according to the reactive voltage sensitivity; and (4) allocating the buses with lower voltage limits according to the active voltage sensitivity and the reactive voltage sensitivity, and allocating the risk indexes to the generators with reactive power less than 0 or the loads with reactive loads more than 0. And further merging and calculating indexes shared by the generator and the load to market main bodies such as power plants participating in transaction, and the like, and the specific steps are as follows:
(3-1) for the bus with the higher voltage limit, the out-of-limit index is for the main body S of power generation (power selling)iThe apportionment calculation method is as follows:
Figure GDA0003473316050000092
in the above formula (2), omegaBFor a set including the current upper bound buses, each of the upper bound buses is denoted as b, RbThe out-of-limit index of the b-th bus bar,
Figure GDA0003473316050000093
for generating (selling) electricityiComprising a generator UnjThe contribution ratio of the bus voltage to the upper limit is calculated by the following method:
Figure GDA0003473316050000094
wherein b is an out-of-limit bus, SBG(b,uj) Generator Un for reactive voltage sensitivity calculationjReactive voltage sensitivity to the b-th bus, uj being the injected bus-in matrix S of the generator UnjBGColumn number in, Qun,jFor the current reactive power output of the generator, Ω z is the set of buses with reactive power injection greater than 0 in the region, QxThe reactive value of the generator or the load is injected into the bus;
in the formula (2), the reaction mixture is,
Figure GDA0003473316050000101
for the main body S of electricity generation (electricity sale)iContaining load LdkThe contribution ratio of the bus voltage to the upper limit is calculated by the following method:
Figure GDA0003473316050000102
wherein b is an out-of-limit bus, SBG(b,dk) Load Ld given for reactive voltage sensitivity calculationkReactive voltage sensitivity to the b-th bus, dkIs a load LdkIs in the matrix SBGColumn number in, Qld,jFor the current reactive power output of the load, Ω z is the set of busbars with reactive injection greater than 0 in the region, QxThe reactive value of the generator or the load is injected into the bus;
(3-2) for the bus with lower voltage limit, the out-of-limit index is for the main body S of power generation (power selling)iThe allocation method comprises the following steps:
Figure GDA0003473316050000103
in the above formula (5), omegaBFor the set of buses including the current lower limit, each of the lower limit buses is denoted as b, RbIs the out-of-limit index of the b-th bus,
Figure GDA0003473316050000104
for generating (selling) electricityiComprising a generator UnjThe contribution ratio to the lower limit of the bus voltage is calculated by:
Figure GDA0003473316050000105
wherein b is an out-of-limit bus, SBG(b,uj) Given for reactive voltage sensitivity calculationReactive voltage sensitivity of generator Unj to b-th bus uj is the injected bus of generator Unj in matrix SBGColumn number in, Qun,jFor the current reactive power output of the generator, Ω z is the set of buses with reactive power injection less than 0 in the region, QxReactive value of the generator or load which is absorbing reactive power from the bus;
in the above-mentioned formula (5),
Figure GDA0003473316050000106
for the main body S of electricity generation (electricity sale)iContaining load LdkThe contribution ratio to the lower limit of the bus voltage is calculated by:
Figure GDA0003473316050000111
wherein b is an out-of-limit bus, SBG(b,dk) Reactive voltage sensitivity, d, of load Unj to the b-th bus given for reactive voltage sensitivity calculationkIs a load LdkIs in the matrix SBGColumn number in, Qld,jFor the current reactive power output of the load, Ω z is the set of busbars with reactive injection less than 0 in the region, QxReactive value of the generator or load which is absorbing reactive power from the bus;
(3-3) bus with higher voltage limit, its out-of-limit index is for main body BiThe apportionment calculation method is as follows:
Figure GDA0003473316050000112
in the above formula (8), omegaBThe method comprises the steps that a set comprising current bus bars exceeding the upper limit is marked as b; rbThe out-of-limit index of the b-th bus bar,
Figure GDA0003473316050000113
for electricity consumption (purchasing electricity) main body BiComprising a generator UnjThe contribution ratio of the bus voltage to the upper limit is calculated by the following method:
Figure GDA0003473316050000114
wherein b is an out-of-limit bus, SBG(b,uj) Generator Un for reactive voltage sensitivity calculationjReactive voltage sensitivity to the b-th bus, uj being the injected bus-in matrix S of the generator UnjBGColumn number in, Qun,jFor the current reactive power output of the generator, Ω z is the set of buses with reactive power injection greater than 0 in the region, QxThe reactive value of the generator or the load is injected into the bus;
in the above-mentioned formula (8),
Figure GDA0003473316050000115
for the main body B of electricity consumption (purchasing electricity)iContaining load LdkThe contribution ratio of the bus voltage to the upper limit is calculated by the following method:
Figure GDA0003473316050000116
wherein b is an out-of-limit bus, SBG(b,dk) Load Ld given for reactive voltage sensitivity calculationkReactive voltage sensitivity to the b-th bus, dkIs a load LdkIs in the matrix SBGColumn number in, Qld,jFor the current reactive power output of the load, Ω z is the set of busbars with reactive injection greater than 0 in the region, QxThe reactive value of the generator or the load is injected into the bus;
(3-4) for the bus with lower voltage limit, the out-of-limit index is for the main body BiThe allocation method comprises the following steps:
Figure GDA0003473316050000121
in the above formula (11), omegaBFor a set of busbars comprising a current lower bound, each lower boundIs denoted as b, RbIs the out-of-limit index of the b-th bus,
Figure GDA0003473316050000122
for electricity consumption (purchasing electricity) main body BiComprising a generator UnjThe contribution ratio to the lower limit of the bus voltage is calculated by:
Figure GDA0003473316050000123
wherein b is an out-of-limit bus, SBG(b,uj) Reactive voltage sensitivity of the generator Unj to the b-th bus given for reactive voltage sensitivity calculation uj is the injected bus-in matrix S of the generator UnjBGColumn number in, Qun,jFor the current reactive power output of the generator, Ω z is the set of buses with reactive power injection less than 0 in the region, QxReactive value of the generator or load which is absorbing reactive power from the bus;
in the above-mentioned formula (12),
Figure GDA0003473316050000124
for the main body B of electricity consumption (purchasing electricity)iContaining load LdkThe contribution ratio to the lower limit of the bus voltage is calculated by:
Figure GDA0003473316050000125
wherein b is an out-of-limit bus, SBG(b,dk) Reactive voltage sensitivity, d, of load Unj to the b-th bus given for reactive voltage sensitivity calculationkIs a load LdkIs in the matrix SBGColumn number in, Qld,jFor the current reactive power output of the load, Ω z is the set of busbars with reactive injection less than 0 in the region, QxReactive value of the generator or load which is absorbing reactive power from the bus;
(4) correcting the network passing fee of the market subject:
according to the calculated voltage out-of-limit indexThe result of the distribution is corrected, and the basic net charge P of the market body is correctediMultiplying the basic net charge by a weight coefficient Fi：
Figure GDA0003473316050000131
Pi＝Pi·Fi (15)
Wherein the content of the first and second substances,
Figure GDA0003473316050000132
the voltage out-of-limit index value allocated to the market body,
Figure GDA0003473316050000133
is the sum of all bus voltage out-of-limit indexes in the current area,
Figure GDA0003473316050000134
Ω z is the set of all bus out-of-limit indicators in the region.
When the power grid normally operates and the voltage is not out of limit, the voltage out-of-limit index is 0, which indicates that the power grid operates within the preset electricity price cost; when the voltage out-of-limit condition occurs, the voltage out-of-limit index is larger than 0, which indicates that the power grid exceeds the preset cost range, the power market transaction may cause the risk of power grid operation, and then the basic power grid fee is intervened or adjusted.
The invention is further described in detail below with reference to one embodiment:
a net-passing fee calculation method considering voltage out-of-limit indexes in an electric power market environment forms an active power flow directed graph by using 4 bus nodes and 5 branches, wherein an arc represents a circuit in a system, the direction of the arc represents the flow direction of power flow, the beginning end of the arc is called an upstream node of the arc, and the terminal end of the arc is called a downstream node of the arc. The method comprises the following steps:
step 1, establishing a market main body model, wherein the market main body model mainly refers to main body objects participating in electric power market transaction, and the market main body model comprises a power generation (electricity selling) object and a power utilization (electricity purchasing) object.
1)Market body model S of power generation (selling) objectiIs defined as:
Si＝{Unj,j＝1,..,J,Pun,j＞0}∪{Ldk,k＝1,..,K,Pld,k＜0}
the power generation (power selling) objects comprise a plurality of generator models with positive active output and a plurality of load models with negative active load. Wherein UnjAnd LdkThe generators and loads in the model are calculated for the grid.
2) Market main body model B of electricity consumption (electricity purchasing) objectiIs defined as:
Bi＝{Unj,j＝1,..,J,Pun,j＜0}∪{Ldk,k＝1,..,K,Pld,k＞0}
the electricity consumption (purchase) objects comprise a plurality of generator models with negative active output and a plurality of load models with positive active load. Wherein UnjAnd LdkThe generators and loads in the model are calculated for the grid.
Step 2, calculating out-of-limit indexes of bus voltage
The voltage index mainly evaluates whether the bus voltage is out of limit, and the index calculation method for a single bus comprises the following steps:
Figure GDA0003473316050000141
wherein V is the current value of the bus voltage, Vmax、VminThe voltage limit is set, k is less than or equal to 1, and k is generally 0.9.
Step 3, apportionment of out-of-limit indexes of bus voltage:
the distribution of the out-of-limit voltage index mainly adopts the sensitivity of active and reactive injection to the bus voltage. For the bus with the voltage higher than the upper limit, the risk index is allocated to the generator with the reactive output larger than 0 or the load with the reactive load smaller than 0 only according to the reactive voltage sensitivity; and (4) allocating the buses with lower voltage limits according to the active voltage sensitivity and the reactive voltage sensitivity, and allocating the risk indexes to the generators with reactive power less than 0 or the loads with reactive loads more than 0. And further merging and calculating indexes shared by the generator and the load to market main bodies such as power plants participating in transaction, and the like, and the specific steps are as follows:
1) for the bus with higher voltage limit, the out-of-limit index is for the main body S of power generation (power selling)iThe apportionment calculation method is as follows:
Figure GDA0003473316050000142
in the formula (2), omegaBThe method comprises the steps that a set comprising current bus bars exceeding the upper limit is marked as b; rbThe out-of-limit index of the b-th bus bar,
Figure GDA0003473316050000143
for generating (selling) electricityiComprising a generator UnjThe contribution ratio of the bus voltage to the upper limit is calculated by the following method:
Figure GDA0003473316050000144
wherein b is an out-of-limit bus, SBG(b,uj) Generator Un for reactive voltage sensitivity calculationjReactive voltage sensitivity to the b-th bus, uj being the injected bus-in matrix S of the generator UnjBGColumn number in, Qun,jFor the current reactive power output of the generator, Ω z is the set of buses with reactive power injection greater than 0 in the region, QxIs the reactive value of the generator or load that injects reactive power into the bus.
In the formula (2), the reaction mixture is,
Figure GDA0003473316050000151
for the main body S of electricity generation (electricity sale)iContaining load LdkThe contribution ratio of the bus voltage to the upper limit is calculated by the following method:
Figure GDA0003473316050000152
wherein b is an out-of-limit bus, SBG(b,dk) Load Ld given for reactive voltage sensitivity calculationkReactive voltage sensitivity to the b-th bus, dkIs a load LdkIs in the matrix SBGColumn number in, Qld,jFor the current reactive power output of the load, Ω z is the set of busbars with reactive injection greater than 0 in the region, QxIs the reactive value of the generator or load that injects reactive power into the bus.
2) For the bus with lower limit voltage, the out-of-limit index is for the main body S of power generation (power selling)iThe allocation method comprises the following steps:
Figure GDA0003473316050000153
in the formula (5), omegaBThe method comprises the steps that a set comprising buses with current lower limits is recorded as b; rbIs the out-of-limit index of the b-th bus,
Figure GDA0003473316050000154
for generating (selling) electricityiComprising a generator UnjThe contribution ratio to the lower limit of the bus voltage is calculated by:
Figure GDA0003473316050000155
wherein b is an out-of-limit bus, SBG(b,uj) Reactive voltage sensitivity of the generator Unj to the b-th bus given for reactive voltage sensitivity calculation uj is the injected bus-in matrix S of the generator UnjBGColumn number in, Qun,jFor the current reactive power output of the generator, Ω z is the set of buses with reactive power injection less than 0 in the region, QxIs the reactive value of the generator or load that absorbs the reactive power from the bus.
In the formula (5), the reaction mixture is,
Figure GDA0003473316050000161
for the main body S of electricity generation (electricity sale)iContaining load LdkThe contribution ratio to the lower limit of the bus voltage is calculated by:
Figure GDA0003473316050000162
wherein b is an out-of-limit bus, SBG(b,dk) Reactive voltage sensitivity, d, of load Unj to the b-th bus given for reactive voltage sensitivity calculationkIs a load LdkIs in the matrix SBGColumn number in, Qld,jFor the current reactive power output of the load, Ω z is the set of busbars with reactive injection less than 0 in the region, QxIs the reactive value of the generator or load that absorbs the reactive power from the bus.
3) For the bus with higher voltage limit, the out-of-limit index is for the main body BiThe apportionment calculation method is as follows:
Figure GDA0003473316050000163
in the formula (8), omegaBThe method comprises the steps that a set comprising current bus bars exceeding the upper limit is marked as b; rbThe out-of-limit index of the b-th bus bar,
Figure GDA0003473316050000164
for electricity consumption (purchasing electricity) main body BiComprising a generator UnjThe contribution ratio of the bus voltage to the upper limit is calculated by the following method:
Figure GDA0003473316050000165
wherein b is an out-of-limit bus, SBG(b,uj) Generator Un for reactive voltage sensitivity calculationjReactive voltage sensitivity to the b-th bus, uj being the injected bus-in matrix S of the generator UnjBGColumn number in, Qun,jTo send outThe current reactive output of the motor, omega z, is a bus set with reactive injection greater than 0 in the region, QxIs the reactive value of the generator or load that injects reactive power into the bus.
In the formula (8), the reaction mixture is,
Figure GDA0003473316050000166
for the main body B of electricity consumption (purchasing electricity)iContaining load LdkThe contribution ratio of the bus voltage to the upper limit is calculated by the following method:
Figure GDA0003473316050000171
wherein b is an out-of-limit bus, SBG(b,dk) Load Ld given for reactive voltage sensitivity calculationkReactive voltage sensitivity to the b-th bus, dkIs a load LdkIs in the matrix SBGColumn number in, Qld,jFor the current reactive power output of the load, Ω z is the set of busbars with reactive injection greater than 0 in the region, QxIs the reactive value of the generator or load that injects reactive power into the bus.
2) For the bus with lower limit of voltage, the out-of-limit index is for the main body BiThe allocation method comprises the following steps:
Figure GDA0003473316050000172
in the formula (10), omegaBThe method comprises the steps that a set comprising buses with current lower limits is recorded as b; rbIs the out-of-limit index of the b-th bus,
Figure GDA0003473316050000173
for electricity consumption (purchasing electricity) main body BiComprising a generator UnjThe contribution ratio to the lower limit of the bus voltage is calculated by:
Figure GDA0003473316050000174
wherein b is an out-of-limit bus, SBG(b,uj) Reactive voltage sensitivity of the generator Unj to the b-th bus given for reactive voltage sensitivity calculation uj is the injected bus-in matrix S of the generator UnjBGColumn number in, Qun,jFor the current reactive power output of the generator, Ω z is the set of buses with reactive power injection less than 0 in the region, QxIs the reactive value of the generator or load that absorbs the reactive power from the bus.
In the formula (10), the compound represented by the formula (10),
Figure GDA0003473316050000175
for the main body B of electricity consumption (purchasing electricity)iContaining load LdkThe contribution ratio to the lower limit of the bus voltage is calculated by:
Figure GDA0003473316050000176
wherein b is an out-of-limit bus, SBG(b,dk) Reactive voltage sensitivity, d, of load Unj to the b-th bus given for reactive voltage sensitivity calculationkIs a load LdkIs in the matrix SBGColumn number in the moment, Qld,jFor the current reactive power output of the load, Ω z is the set of busbars with reactive injection less than 0 in the region, QxIs the reactive value of the generator or load that absorbs the reactive power from the bus.
Step 4, correcting the network passing fee of the market subject
Correcting the basic net charge P of the market subject according to the calculated distribution result of the voltage out-of-limit indexiMultiplying the basic net charge by a weight coefficient Fi：
Figure GDA0003473316050000181
Pi＝Pi·Fi (14)
Wherein the content of the first and second substances,
Figure GDA0003473316050000182
the voltage out-of-limit index value allocated to the market body,
Figure GDA0003473316050000183
is the sum of all bus voltage out-of-limit indexes in the current area,
Figure GDA0003473316050000184
Ω z is the set of all bus out-of-limit indicators in the region.
When the power grid normally operates and the voltage is not out of limit, the voltage out-of-limit index is 0, which indicates that the power grid operates within the preset electricity price cost; when the voltage out-of-limit condition occurs, the voltage out-of-limit index is larger than 0, which indicates that the power grid exceeds the preset cost range, the power market transaction may cause the risk of power grid operation, and then the basic power grid fee is intervened or adjusted.
An embodiment of the method of the invention is described in detail below:
step 1, establishing a market main body model S of a power generation (power selling) object1And S2Wherein S is1Comprising a generator G1，S2Comprising a generator G2(ii) a Market main body model B for establishing electricity consumption (electricity purchasing) object3And B4Wherein B is3Including the load Ld3，B4Including the load Ld4。
Wherein, the electricity generation (selling electricity) market main body S1Active power injection to bus 1
Figure GDA0003473316050000185
Injected reactive power
Figure GDA0003473316050000186
Power generation (selling electricity) market main body S2Active power injection to bus 2
Figure GDA0003473316050000187
Absorbing reactive power
Figure GDA0003473316050000188
Electricity consumption (purchasing) market main body B3Absorb active power to the bus bar 3
Figure GDA0003473316050000189
Absorbing reactive power
Figure GDA00034733160500001810
Electricity consumption (purchasing) market main body B4Absorb active power to the bus 4
Figure GDA00034733160500001811
Injected reactive power
Figure GDA00034733160500001812
Step 2, calculating out-of-limit indexes of bus voltage
The control target coefficient k is set to 0.996.
Setting the current voltage value V of the bus 11115.35kV, upper limit of voltage V1 max117kV, lower voltage limit V1 min＝113kV；
Setting the current voltage value V of the bus 22114.83kV, upper limit of voltage V2 max117kV, lower voltage limit V2 min＝113kV；
Setting the current voltage value V of the bus 33112.52kV, upper voltage limit V3 max117kV, lower voltage limit V3 min＝113kV；
Setting the current voltage value V of the bus 44118.49kV, upper limit of voltage V4 max117kV, lower voltage limit V4 min＝113kV；
The out-of-limit indexes of the voltage of each bus node are respectively as follows:
R1＝0；R2＝0；
Figure GDA0003473316050000191
Figure GDA0003473316050000192
step 3, apportionment of out-of-limit indexes of bus voltage:
according to a quasi-steady-state sensitivity calculation method provided by Sunnmacro, Zenberming and phase-year Ded in a quasi-steady-state sensitivity analysis method (the Chinese Motor engineering project, V19N 4, 1999 4 months, pp.9-13), a reactive voltage sensitivity matrix S of each power generation main body and each load to a bus is obtainedBGThe following were used:
Figure GDA0003473316050000193
the out-of-limit index of the voltage with the upper limit is distributed to a power generation (power selling) main body S1And S2The method comprises the following steps:
since the upper limit index of the current bus 4 is not 0, the current bus needs to be allocated to a power generation (power selling) main body, and the calculation process is as follows:
for the main body S of electricity generation (selling electricity)1The calculation process of the upper limit-exceeding index of the apportionment bus 4 is as follows:
Figure GDA0003473316050000194
wherein Qun,1＝10MVar，SBG(b,uj)＝SBG(4,1) ═ 0.03, and reactive load Q is injected into the busL4Sensitivity value S for bus bar 4 of-10 MVarBGAnd (4,4) ═ 0.1, and the reactive loss of the line is ignored.
Due to the main body S of power generation (selling electricity)1Does not contain a load, therefore
Figure GDA0003473316050000201
Finally, S1The out-of-limit indexes of the apportioned out-of-voltage upper limit are as follows:
Figure GDA0003473316050000202
for the main body S of electricity generation (selling electricity)2The calculation process of the upper limit-exceeding index of the apportionment bus 4 is as follows:
due to the main body S of power generation (selling electricity)2The method does not include a unit with positive reactive power output and a load, so that the distributed out-of-limit index of the over-voltage upper limit
Figure GDA0003473316050000203
The lower limit of the voltage out-of-limit index is distributed to the main body S of power generation (power selling)1And S2The method comprises the following steps:
since the lower limit index of the current bus 3 is not 0, the current bus needs to be allocated to a power generation (power selling) main body, and the calculation process is as follows:
for the main body S of electricity generation (selling electricity)1The lower limit index calculation process of the apportionment bus 3 is as follows:
due to the main body S of power generation (selling electricity)1The method does not include a unit with negative reactive power output and a load, so that the distributed out-of-limit index of the lower limit of the over-voltage
Figure GDA0003473316050000204
For the main body S of electricity generation (selling electricity)2The lower limit index calculation process of the apportionment bus 3 is as follows:
Figure GDA0003473316050000205
wherein Qun,2＝-5MVar，SBG(b,uj)＝SBG(3,2) ═ 0.03, absorb the reactive load Q to the busL315MVar, sensitivity value S to bus 3BGAnd (3,3) is 0.1, and the reactive loss of the line is ignored.
Due to the main body S of power generation (selling electricity)2Does not contain a load, therefore
Figure GDA0003473316050000211
Finally, S2The out-of-limit indexes of the lower limit of the distributed over-voltage are as follows:
Figure GDA0003473316050000212
the out-of-limit voltage index of the upper limit is distributed to the main body B of electricity consumption (electricity purchasing)3And B4The method comprises the following steps:
since the upper limit index of the current bus 4 is not 0, the current bus needs to be allocated to the electricity consumption (electricity purchase) main body, and the calculation process is as follows:
for electricity consumption (purchase) main body B3The calculation process of the upper limit-exceeding index of the apportionment bus 4 is as follows:
because the electricity consumption (purchasing electricity) main body B3The device does not comprise a unit and does not comprise a load with reactive power less than 0, so that the distributed out-of-limit index of the over-voltage upper limit
Figure GDA0003473316050000213
For electricity consumption (purchase) main body B4The calculation process of the upper limit-exceeding index of the apportionment bus 4 is as follows:
Figure GDA0003473316050000214
wherein Qld,4＝-10MVar，SBG(b,dk)＝SBG(4,4) ═ 0.1, inject the idle unit Q to the generating line1Sensitivity S to bus bar 4 of 10MVarBG(b,uj)＝SBGAnd (4,1) ═ 0.03, and the reactive loss of the line is ignored.
Because the electricity consumption (purchasing electricity) main body B4In the middle of the system, no machine set is included, so
Figure GDA0003473316050000215
Finally, B4The out-of-limit indexes of the apportioned out-of-voltage upper limit are as follows:
Figure GDA0003473316050000216
the lower limit of the voltage out-of-limit index is distributed to the electricity consumption (electricity purchasing) body B3And B4The method comprises the following steps:
since the lower limit index of the current bus 3 is not 0, the current bus needs to be allocated to the electricity consumption (electricity purchase) main body, and the calculation process is as follows:
for electricity consumption (purchase) main body B3The out-of-limit index calculation process of the apportionment bus 3 is as follows:
Figure GDA0003473316050000221
wherein Qld,3＝15MVar，SBG(b,dk)＝SBG(3,3) ═ 0.1, unit Q absorbing reactive power to bus2Sensitivity to bus 3S of-5 MVarBG(b,uj)＝SBGAnd (3,2) ═ 0.03, and the reactive loss of the line is ignored.
Because the electricity consumption (purchasing electricity) main body B3In the middle of the system, no machine set is included, so
Figure GDA0003473316050000222
Finally, B3The out-of-limit indexes of the lower limit of the distributed over-voltage are as follows:
Figure GDA0003473316050000223
for electricity consumption (purchase) main body B4The lower limit index calculation process of the apportionment bus 3 is as follows:
because the electricity consumption (purchasing electricity) main body B4The method does not include a unit or a load with positive reactive power output, so that the distributed out-of-limit index of the lower limit of the over-voltage
Figure GDA0003473316050000224
Step 4, correcting the network passing fee of the market subject
The basic internet fee P of the market subject is set to be 0.18 yuan.
The running cost of the power grid after the voltage is out of limit is distributed to a power generation (power selling) main body S1And S2The above results are as follows:
for the main body S of electricity generation (selling electricity)1Corrected net charge P1 s:
Figure GDA0003473316050000225
For the main body S of electricity generation (selling electricity)2Corrected net charge
Figure GDA0003473316050000231
Figure GDA0003473316050000232
The running cost of the power grid after the voltage is out of limit is distributed to a main body B of power consumption (electricity purchasing)3And B4The above results are as follows:
for electricity consumption (purchase) main body B3Corrected net charge
Figure GDA0003473316050000233
Figure GDA0003473316050000234
For electricity consumption (purchase) main body B4Corrected net charge
Figure GDA0003473316050000235
Figure GDA0003473316050000236

Claims (1)

1. A method for calculating a net charge considering a voltage out-of-limit index in an electric power market environment is characterized by comprising the following steps:
(1) establishing a market main body model, wherein the market main body model mainly refers to a main body object participating in electric power market transaction and comprises an electricity selling main body and an electricity purchasing main body:
(1-1) Electricity selling body SiIs defined as:
Si＝{Unj,j＝1,..,J,PUn,j＞0}∪{Ldk,k＝1,..,K,PLd,k＜0}，
electricity selling main body SiThe method comprises a generator model with positive active output in a power grid calculation model and a load model with negative active load, wherein UnjExpressed as the jth generator with positive active power output, Ld, in the grid calculation modelkThe k-th load is expressed as the negative active load in the power grid calculation model;
(1-2) Electricity purchasing entity GiIs defined as:
Gi＝{Unj’,j’＝1,..,J’,PUn,j’＜0}∪{Ldk’,k’＝1,..,K’,PLd,k’＞0}，
electricity purchase main body GiThe method comprises a generator model with negative active output in a power grid calculation model and a load model with positive active load, wherein Unj’Expressed as the j' th generator with negative active power output, Ld in the power grid calculation modelk’The k' th load with positive active load in the power grid calculation model is represented;
(2) calculating a bus voltage out-of-limit index R:
the voltage out-of-limit index mainly evaluates whether the bus voltage is out of limit, and the index calculation method for a single bus comprises the following steps:
Figure FDA0003465964820000011
wherein V is the current value of the bus voltage, VmaxIs the upper limit value of the bus voltage, VminAs a lower limit value of the bus voltageTaking m as 0.9;
(3) the method for apportioning out-of-limit indexes of the bus voltage comprises the following specific steps:
(3-1) for the bus with the higher voltage limit, the out-of-limit index is for the electricity selling main body SiThe apportionment calculation method is as follows:
Figure FDA0003465964820000021
in the above formula (2), omegaBSHTo include the set of current upper bound buses, each upper bound bus is denoted as b, RbThe out-of-limit index of the b-th bus bar,
Figure FDA0003465964820000022
for selling electricity main body SiComprising a generator UnjThe contribution ratio of the bus voltage to the upper limit is calculated by the following method:
Figure FDA0003465964820000023
wherein S isBG(b,uj) Generator Un for reactive voltage sensitivity calculationjReactive voltage sensitivity to the b-th bus, ujFor a generator UnjIs in the matrix SBGColumn number in, QUn,jFor a generator UnjCurrent reactive output, ΩZSHFor a set of busbars with reactive injection greater than 0 in a zone, Qz,xThe reactive value of the generator or the load is injected into the bus;
in the formula (2), the reaction mixture is,
Figure FDA0003465964820000024
is the electricity selling main body SiContaining load LdkThe contribution ratio of the bus voltage to the upper limit is calculated by the following method:
Figure FDA0003465964820000025
wherein S isBG(b,dk) Load Ld given for reactive voltage sensitivity calculationkReactive voltage sensitivity to the b-th bus, dkIs a load LdkIs in the matrix SBGColumn number in, QLd,kIs a load LdkCurrent reactive output, ΩZSHFor a set of busbars with reactive injection greater than 0 in a zone, Qz,xThe reactive value of the generator or the load is injected into the bus;
(3-2) for the bus with lower voltage limit, the out-of-limit index is for the electricity selling main body SiThe allocation method comprises the following steps:
Figure FDA0003465964820000026
in the above formula (5), omegaBSLFor a set of busbars comprising a current lower bound, each lower bound busbar is denoted b', Rb’For the out-of-limit indicator of the b' th bus,
Figure FDA0003465964820000031
for selling electricity main body SiComprising a generator UnjThe contribution ratio to the lower limit of the bus voltage is calculated by:
Figure FDA0003465964820000032
wherein S isBG(b’,uj) Generator Un for reactive voltage sensitivity calculationjReactive voltage sensitivity to the b' th bus, ujFor a generator UnjIs in the matrix SBGColumn number in, QUn,jFor a generator UnjCurrent reactive output, ΩFSLFor a set of busbars with reactive injection less than 0 in a zone, Qr,xReactive value of the generator or load which is absorbing reactive power from the bus;
in the above-mentioned formula (5),
Figure FDA0003465964820000033
is the electricity selling main body SiContaining load LdkThe contribution ratio to the lower limit of the bus voltage is calculated by:
Figure FDA0003465964820000034
wherein S isBG(b’,dk) Load Ld given for reactive voltage sensitivity calculationkReactive voltage sensitivity to the b' th bus, dkIs a load LdkIs in the matrix SBGColumn number in, QLd,kIs a load LdkCurrent reactive output, ΩFSLFor a set of busbars with reactive injection less than 0 in a zone, Qr,xReactive value of the generator or load which is absorbing reactive power from the bus;
(3-3) for the bus with higher voltage limit, the out-of-limit index thereof is for the electricity purchasing subject GiThe apportionment calculation method is as follows:
Figure FDA0003465964820000035
in the above formula (8), omegaBSHTo include the set of current upper bound buses, each upper bound bus is denoted as b, RbThe out-of-limit index of the b-th bus bar,
Figure FDA0003465964820000036
for purchasing electricity main body GiComprising a generator Unj'The contribution ratio of the bus voltage to the upper limit is calculated by the following method:
Figure FDA0003465964820000041
wherein S isBG(b,uj') Generator Un for reactive voltage sensitivity calculationj'Reactive voltage sensitivity to the b-th bus, uj'For a generator Unj'Is in the matrix SBGColumn number in, QUn,j'For the current reactive power output of the generator, omegaZSHFor a set of busbars with reactive injection greater than 0 in a zone, Qz,xThe reactive value of the generator or the load is injected into the bus;
in the above-mentioned formula (8),
Figure FDA0003465964820000042
to the electricity purchasing main body GiContaining load Ldk'The contribution ratio of the bus voltage to the upper limit is calculated by the following method:
Figure FDA0003465964820000043
wherein S isBG(b,dk') Load Ld given for reactive voltage sensitivity calculationk'Reactive voltage sensitivity to the b-th bus, dk'Is a load Ldk'Is in the matrix SBGColumn number in, QLd,k'For the current reactive power output of the load, omegaZSHFor a set of busbars with reactive injection greater than 0 in a zone, Qz,xThe reactive value of the generator or the load is injected into the bus;
(3-4) for the bus with lower voltage limit, the out-of-limit index is for the electricity purchasing subject GiThe allocation method comprises the following steps:
Figure FDA0003465964820000044
in the above formula (11), omegaBSLFor a set of busbars comprising a current lower bound, each lower bound busbar is denoted b', Rb'For the out-of-limit indicator of the b' th bus,
Figure FDA0003465964820000045
for purchasing electricity main body GiComprising a generator Unj'The contribution ratio to the lower limit of the bus voltage is calculated by:
Figure FDA0003465964820000046
wherein S isBG(b',uj') Generator Un for reactive voltage sensitivity calculationj'Reactive voltage sensitivity to the b' th bus, uj'For a generator Unj'Is in the matrix SBGColumn number in, QUn,j'For a generator Unj'Current reactive output, ΩFSLFor a set of busbars with reactive injection less than 0 in a zone, Qr,xReactive value of the generator or load which is absorbing reactive power from the bus;
in the above-mentioned formula (11),
Figure FDA0003465964820000051
to the electricity purchasing main body GiContaining load Ldk'The contribution ratio to the lower limit of the bus voltage is calculated by:
Figure FDA0003465964820000052
wherein S isBG(b',dk') Load Ld given for reactive voltage sensitivity calculationkReactive voltage sensitivity to the b' th bus, dk'Is a load Ldk'Is in the matrix SBGColumn number in, QLd,k'Is a load Ldk'Current reactive output, ΩFSLFor a set of busbars with reactive injection less than 0 in a zone, Qr,xReactive value of the generator or load which is absorbing reactive power from the bus;
(4) correcting the network passing fee of the market subject:
according to the calculated distribution result of the voltage out-of-limit index, the city is correctedBasic passing fee P of field bodyiMultiplying the basic net charge by a weight coefficient Fi：
Figure FDA0003465964820000053
Pi＝Pi·Fi (15)
Wherein the content of the first and second substances,
Figure FDA0003465964820000054
a voltage out-of-limit index value apportioned to the market subject
Figure FDA0003465964820000055
Or
Figure FDA0003465964820000056
Figure FDA0003465964820000057
Is the sum of all bus voltage out-of-limit indexes in the current area,
Figure FDA0003465964820000058
Ω z is the set of all bus out-of-limit indicators in the region.
CN201811465947.4A 2018-12-03 2018-12-03 Network charge calculation method considering voltage out-of-limit index under electric power market environment Expired - Fee Related CN109473981B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811465947.4A CN109473981B (en) 2018-12-03 2018-12-03 Network charge calculation method considering voltage out-of-limit index under electric power market environment

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811465947.4A CN109473981B (en) 2018-12-03 2018-12-03 Network charge calculation method considering voltage out-of-limit index under electric power market environment

Publications (2)

Publication Number Publication Date
CN109473981A CN109473981A (en) 2019-03-15
CN109473981B true CN109473981B (en) 2022-03-15

Family

ID=65674862

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811465947.4A Expired - Fee Related CN109473981B (en) 2018-12-03 2018-12-03 Network charge calculation method considering voltage out-of-limit index under electric power market environment

Country Status (1)

Country Link
CN (1) CN109473981B (en)

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107181284A (en) * 2017-04-24 2017-09-19 中国电力科学研究院 The method of adjustment and device of a kind of out-of-limit electricity of circuit

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107181284A (en) * 2017-04-24 2017-09-19 中国电力科学研究院 The method of adjustment and device of a kind of out-of-limit electricity of circuit

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
Allocation of the Reactive Power Support Requirements in Multitransaction Networks;Gianfranco Chicco等;《IEEE Transactions on Power Systems》;20020807;第243-249页 *
计及电压越限投资的有源配电网节点电价计算方法;范绚然等;《电力建设》;20171031;第38卷(第10期);第123-128页 *

Also Published As

Publication number Publication date
CN109473981A (en) 2019-03-15

Similar Documents

Publication Publication Date Title
Yan et al. Low-carbon economic dispatch method for integrated energy system considering seasonal carbon flow dynamic balance
Jiang et al. A multiagent-based hierarchical energy management strategy for maximization of renewable energy consumption in interconnected multi-microgrids
López et al. V2G strategies for congestion management in microgrids with high penetration of electric vehicles
Wang et al. Transactive energy sharing in a microgrid via an enhanced distributed adaptive robust optimization approach
CN111079971A (en) Charging station pricing method considering vehicle, station and network
KR20210058633A (en) Power trading apparatus between prosumer and consumer and its method
Liu et al. Wind power and photovoltaic power: How to improve the accommodation capability of renewable electricity generation in China?
Zhang et al. Research on the optimal allocation method of source and storage capacity of integrated energy system considering integrated demand response
CN107742223B (en) Provincial power grid power transmission and distribution pricing method considering power grid characteristics
CN114881749A (en) Virtual power plant participation main and auxiliary market bidding method considering risk constraint
Fu et al. Optimal scheduling strategy of grid-connected microgrid with ladder-type carbon trading based on Stackelberg game
Kong et al. Independence enhancement of distributed generation systems by integrating shared energy storage system and energy community with internal market
Cheng et al. Optimized configuration of distributed power generation based on multi-stakeholder and energy storage synergy
CN109473981B (en) Network charge calculation method considering voltage out-of-limit index under electric power market environment
CN111160767A (en) Comprehensive energy service benefit evaluation method
CN114693346A (en) Electric power market clearing method, system, equipment and storage medium
Liang et al. An Optimization Method for the Distributed Collaborative Operation of Multilateral Entities Considering Dynamic Time-of-Use Electricity Price in Active Distribution Network
Zhang et al. Multi-energy complementary virtual power plant economic scheduling considering demand response
Liang et al. A non-cooperative game theory based energy trading strategy of multi-microgrids
Wu et al. Multi-alliance market subject auxiliary peak shaving strategy for new energy consumption
CN117808565B (en) Virtual power plant multi-time bidding method considering green evidence and carbon transaction
CN112270432B (en) Energy management method of comprehensive energy system considering multi-subject benefit balance
Jiang et al. Two-stage optimization model of electric vehicle participating in balancing market
Hu et al. Modeling and Strategic Design of Multi-Layered Aggregation and Control for Flexible Loads
Wang et al. Consider the Way of Gathering Network Tariff for High Frequency Distributed Electricity Trading in Distribution Networks

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
GR01 Patent grant
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20220315

CF01 Termination of patent right due to non-payment of annual fee