CN101261291B - Wattless service pricing method based on fees susceptance - Google Patents

Abstract

The invention discloses a wattless service pricing method based on a fee susceptance, which includes the following steps: 1) the power supply watt, the load watt and the load idle work of a wattless fee paying subject are equivalently converted into a node injection current of a matching current of a steady state; the power idle work of the wattless fee charging subject is equivalently converted into a node parallel susceptance of a matching current of a steady state, that is, a fee paying susceptance; 2) the reactive component of each user on the fee susceptance is calculated; 3) the wattlessproduction cost weight of the fee susceptance is obtained by multiplying each reactive component of the fee susceptance by the wattless production cost of the wattless power unit relative to the feesusceptance; 4) the wattless service electricity price of each user relative to watt and idle work is calculated. The wattless service pricing method involved in the invention is suitable for any structure power net with advance and/or lagging non-power sources. By using the wattless service pricing method, not only the electricity grid is guaranteed to have the balance of income and outlay, but also the users of the electricity grid are provided with precise economic incentive information.

Description

Wattless service pricing method based on fees susceptance

Technical field

The present invention relates to a kind of electric power pricing method, the Wattless service pricing method in particularly a kind of electric energy transmitting process.

Background technology

The market-oriented reform of electric power was risen from early 1990s, have swept the globe at present.The electricity price of the electricity market of tending to act running, the broad applicability of its theoretical tightness, algorithm, excitation user make full use of existing resources effective etc., are to guarantee the electricity market key of healthy development.But because capacitive reactive power in opposition to each other and the idle complicacy that all will charge of perception and the unintentional nonlinearity of power expression own make existing idle service electricity price theory seriously lack with algorithm.

Wattless service pricing method mainly contains marginal costing and integrated cost method two big classes at present.Marginal costing contains the cost minimization model of node reactive balance constraint by foundation, and utilization optimization method solving model is again with the shadow price of each node reactive balance constraint of the model solution correspondence Reactive Power Price as this node.This is a kind of method with strict economic theory basis.It can provide good economic incentives information, promotes the optimization utilization of electric power resource, but can not realize the balance between revenue and expenditure of idle cost, and along with objective function in the model and constraint condition consider different, the Reactive Power Price of each node has different.

The integrated cost method is to give each bar circuit with the total idle absorption of costs of system earlier, distribute to each load by the composition (following the trail of by the following current to load) of circuit trend again, the idle cost of system of asking each node specific load to share promptly gets the payable idle service electricity price of this load.This method thinking is succinct, directly perceived, easy to understand and can guarantee balance between revenue and expenditure, but it is with the idle cost summation back reallocation of system, obscured the Different Effects of the characteristic of opposition naturally of the different location of reactive power source and capacitive thereof and perception to each node Reactive Power Price, the while, it can not provide effective economic incentives information to the user.

Therefore, tightly, extensively be suitable for, have the idle service pricing problem of good excitation information to require study.

Summary of the invention

The purpose of this invention is to provide a kind of user oriented idle service pricing new method--based on the Wattless service pricing method of fees susceptance, this new method not only is suitable for and contains the electrical network of capacitive reactive power power supply, also be suitable for the arbitrary structures electrical network that contains perceptual reactive power source or contain these two kinds of reactive power sources simultaneously, and it can be guaranteed the balance between revenue and expenditure of idle cost, can provide effective economic incentives information to the user again.

For achieving the above object, the technical solution adopted in the present invention is:

A kind of Wattless service pricing method based on fees susceptance may further comprise the steps:

1) power supply of idle paying main body is meritorious, load is meritorious and the equivalent respectively node injection current that becomes corresponding stable state trend to separate of reactive load; With the idle equivalent node shunt susceptance that becomes corresponding stable state trend to separate of idle charge-collector's power supply, i.e. fees susceptance;

2) calculate the idle component of each user on fees susceptance;

3) the idle cost of unit that each idle component on the fees susceptance be multiply by the reactive power source of this fees susceptance correspondence gets the idle cost component on the fees susceptance;

4) calculate each user with respect to meritorious and idle idle service electricity price.

Wherein, step 1) may further comprise the steps:

(i) the stable state trend of calculating electrical network is separated;

(ii) the user of i node is injected the active power P of electrical network _iAnd reactive power Q _iIn pairs equivalent respectively

The node injection current I that answers the stable state trend to separate _i ^pAnd I _i ^q:

${\stackrel{\‾}{I}}_i^p=P_i/{\hat{V}}_i\left(\mathrm{kA}\right)$

${\stackrel{\‾}{I}}_i^q=-j{Q}_i/{\hat{V}}_i\left(\mathrm{kA}\right)$

Wherein

Be i node voltage V _iConjugate;

(iii) with the idle charge-collector's of k node reactive power generation power Q _kThe node fees susceptance that the corresponding stable state trend of equivalent one-tenth is separated:

$B_k=Q_k/{\stackrel{\‾}{V}}_k^2\left(S\right)$

Wherein, V _kBe the k node voltage;

Step 2 wherein) may further comprise the steps:

(i) the meritorious equivalent electric current I of calculating user i _i ^pAt fees susceptance B _kOn the reactive power generation component Q that causes _Ik ^p:

$Q_{\mathrm{ik}}^p=B_k[\left({\stackrel{\‾}{V}}_k\right)\·\left({\stackrel{\‾}{Z}}_{\mathrm{ki}}{\stackrel{\‾}{I}}_i^p\right)]\left(\mathrm{MVAR}\right)$

Wherein " " represents dot-product, Z _KiElement for the capable i row of grid nodes impedance matrix k.

(ii) calculate the idle equivalent electric current I of user i _i ^qAt fees susceptance B _kOn the reactive power generation component Q that causes _Ik ^q:

$Q_{\mathrm{ik}}^q=B_k[\left({\stackrel{\‾}{V}}_k\right)\·\left({\stackrel{\‾}{Z}}_{\mathrm{ki}}{\stackrel{\‾}{I}}_i^q\right)]\left(\mathrm{MVAR}\right)$

Wherein " " represents dot-product, Z _KiElement for the capable i row of grid nodes impedance matrix k.

Step 3) may further comprise the steps:

(i) the meritorious equivalent electric current I of utilization user i _i ^pAt fees susceptance B _kOn the reactive power generation component Q that causes _Ik ^p, calculating corresponding reactive power generation cost component is γ _kQ _Ik ^p(unit); γ wherein _kThe idle cost of unit for k node reactive power source.

(ii) use the idle equivalent electric current I of user i _i ^qAt fees susceptance B _kOn the reactive power generation component Q that causes _Ik ^q, calculating corresponding reactive power generation cost component is γ _kQ _Ik ^q(unit); γ wherein _kThe idle cost of unit for k node reactive power source.

Described step 4) is: meritorious, the idle idle service electricity price λ that calculates transmission user i _i ^pAnd λ _i ^q:

${\mathrm{\λ}}_i^p=\underset{k}{\mathrm{\Σ}}{\mathrm{\γ}}_k{Q}_{\mathrm{ik}}^p/P_i=\underset{k}{\mathrm{\Σ}}{\mathrm{\γ}}_k{B}_k[\left({\stackrel{\‾}{V}}_k\right)\·({\stackrel{\‾}{Z}}_{\mathrm{ki}}/{\hat{V}}_i)]$ (unit/MWh)

${\mathrm{\λ}}_i^q=\underset{k}{\mathrm{\Σ}}{\mathrm{\γ}}_k{Q}_{\mathrm{ik}}^q/Q_i=\underset{k}{\mathrm{\Σ}}{\mathrm{\γ}}_k{B}_k[\left({\stackrel{\‾}{V}}_k\right)\·(-j{\stackrel{\‾}{Z}}_{\mathrm{ki}}/{\hat{V}}_i)]$ (unit/MVARh)

Adopt the beneficial effect that the present invention obtained to have:

1) utilization branch road dissipated power component lapses to algorithm, realized that power supply is meritorious, load is meritorious, the calculating of the idle component of reactive load (user) on (providing idle service) on the fees susceptance, different location and the different attribute (capacitive and perception) of having taken into account reactive power source to the Different Effects of the idle service electricity price of each node, make that Wattless service pricing method provided by the invention is more fair and rational.The present invention has considered the idle service that the meritorious transmission of power supply is accepted, and has quantitatively taken into account the cross influence (being the P-Q coupling) between meritorious and idle.The present invention had both guaranteed balance between revenue and expenditure, and economic incentives information accurately is provided again.Therefore, Wattless service pricing method provided by the invention is more effective, more adapts to the actual needs of electricity market.

(2) the present invention is applied widely, line mutual-ground capacitor need not be done approximate processing, not only is suitable for the electrical network that contains the capacitive reactive power power supply, also is suitable for the arbitrary structures electrical network that contains perceptual reactive power source or contain these two kinds of reactive power sources simultaneously.

(3) be easy to implement, computing formula of the present invention is simple, and the interative computation or the convergence that do not contain any uncertain number of times are judged computing, and computing velocity is fast.

The inventive method is separated based on trend, set up the equivalent electrical network that contains fees susceptance, use the algorithm that lapses to of branch road dissipated power component, realized that power supply is meritorious, load is meritorious, the idle component of reactive load (user) on fees susceptance (idle service is provided) calculated.The idle cost of unit in conjunction with each charge reactive power source accurately calculates the idle cost of serving component that each user deals with, and calculates meritorious, the idle idle service electricity price of each user based on this cost component.The model that provides is simple and clear, line mutual-ground capacitor need not be done power supply or migration process, and a demand is separated simple sinusoidal steady-state circuit.The present invention has simultaneously considered between power supply and meritorious-idle cross influence, makes Reactive Power Pricing more fair and rational, and it can also provide effective economic incentives information to the user simultaneously.

In sum, the present invention has overcome the deficiencies in the prior art, compared to existing technology, has significant theory and technology advantage, has very high using value.

Description of drawings

Fig. 1 is simple 3 node power system equivalent electrical networks.

Embodiment

A kind of Wattless service pricing method based on fees susceptance, step is as follows:

With the node injection current that the power supply of idle paying main body is gained merit, load is meritorious, the corresponding stable state trend of the equivalent respectively one-tenth of reactive load (being called for short " user ") is separated; With the equivalent node shunt susceptance (being called for short " fees susceptance ") that becomes corresponding stable state trend to separate of idle charge-collector's power supply idle (be called for short " idle service ").

In the equivalent electrical network of electric system as shown in Figure 1, I _i ^pAnd I _i ^qThe user who is respectively the i node injects that the meritorious of electrical network and the corresponding electrical network stable state of reactive power trend separate waits the value node injection current.B _kWait value node fees susceptance (equivalent susceptance) for what the idle corresponding stable state trend that the charge reactive power source sent of k node was separated.If node is non-transformer but also zero load not only, then its injection current and fees susceptance do not exist.

(1) determines that with ordinary tides string routine or state estimation program the stable state trend of electrical network separates; Specific algorithm has: Gauss-Seidel procedure, the inferior method of newton-pressgang, P-Q decomposition method etc.Specific algorithm sees also following document: document [1], " Power System Analysis ", and He Yangzan, Wen Zengyin work, the Central China University of Science and Technology publishes, and 2002,52～72 pages.

(2) user of i node is injected the active power P of electrical network _iAnd reactive power Q _iThe corresponding stable state trend of equivalent respectively one-tenth is separated (i node voltage V _iConjugation use

Expression) node injection current I _i ^pAnd I _i ^qTheir concrete computing formula is:

${\stackrel{\‾}{I}}_i^p=P_i/{\hat{V}}_i\left(\mathrm{kA}\right)---\left(1\right)$

${\stackrel{\‾}{I}}_i^q=-j{Q}_i/{\hat{V}}_i\left(\mathrm{kA}\right)---\left(2\right)$

(3) with the idle charge-collector's of k node reactive power generation power Q _kThe corresponding stable state trend of equivalent one-tenth is separated (k node voltage usefulness

Expression) node fees susceptance B _kIts concrete computing formula is:

$B_k=Q_k/{\stackrel{\‾}{V}}_k^2\left(S\right)---\left(3\right)$

Wherein, V _kBe the k node voltage;

(4) the meritorious equivalent electric current I that lapses to algorithm computation user i of utilization branch road dissipated power component _i ^pAt fees susceptance B _kOn the reactive power generation component Q that causes _Ik ^pAnd cost component.Its concrete computing formula is:

$Q_{\mathrm{ik}}^p=B_k[\left({\stackrel{\‾}{V}}_k\right)\·\left({\stackrel{\‾}{Z}}_{\mathrm{ki}}{\stackrel{\‾}{I}}_i^p\right)]\left(\mathrm{MVAR}\right)---\left(4\right)$

Corresponding reactive power generation cost component is γ _kQ _Ik ^p(unit).Wherein " " represents dot-product, Z _KiBe the element of the capable i row of grid nodes impedance matrix k, γ _kThe idle cost of unit for k node reactive power source.This part is calculated and is related to list of references--document [2]: Jian-Chun Peng (Peng Jianchun), Hui Jiang (Jiang Hui) andYong-Hua Song.A Weakly Conditioned Imputation Of An Impedance BranchDissipation Power[J] .IEEE Trans.Power Syst., Nov., 2007,22 (4): 2124-2133

(5) the idle equivalent electric current I that lapses to algorithm computation user i of utilization branch road dissipated power component _i ^qAt fees susceptance B _kOn the reactive power generation component Q that causes _Ik ^qAnd cost component.Its concrete computing formula is:

$Q_{\mathrm{ik}}^q=B_k[\left({\stackrel{\‾}{V}}_k\right)\·\left({\stackrel{\‾}{Z}}_{\mathrm{ki}}{\stackrel{\‾}{I}}_i^q\right)]\left(\mathrm{MVAR}\right)---\left(5\right)$

Corresponding reactive power generation cost component is γ _kQ _Ik ^q(unit).

(6) the idle cost component of all fees susceptances adds up respectively and meritorious, idle divided by user i in meritorious, the idle electrical network that causes of each user i, promptly get the gaining merit of transmission user i, idle idle service electricity price λ _i ^pAnd λ _i ^qTheir concrete computing formula is:

${\mathrm{\λ}}_i^p=\underset{k}{\mathrm{\Σ}}{\mathrm{\γ}}_k{Q}_{\mathrm{ik}}^p/P_i=\underset{k}{\mathrm{\Σ}}{\mathrm{\γ}}_k{B}_k[\left({\stackrel{\‾}{V}}_k\right)\·({\stackrel{\‾}{Z}}_{\mathrm{ki}}/{\hat{V}}_i)]$ (unit/MWh) (6)

${\mathrm{\λ}}_i^q=\underset{k}{\mathrm{\Σ}}{\mathrm{\γ}}_k{Q}_{\mathrm{ik}}^q/Q_i=\underset{k}{\mathrm{\Σ}}{\mathrm{\γ}}_k{B}_k[\left({\stackrel{\‾}{V}}_k\right)\·(-j{\stackrel{\‾}{Z}}_{\mathrm{ki}}/{\hat{V}}_i)]$ (unit/MVARh) (7)

Claims (5)

1. Wattless service pricing method based on fees susceptance may further comprise the steps:

1) power supply of idle paying main body is meritorious, load is meritorious and the equivalent respectively node injection current that becomes corresponding stable state trend to separate of reactive load; With the idle equivalent node shunt susceptance that becomes corresponding stable state trend to separate of idle charge-collector's power supply, i.e. fees susceptance;

2) calculate the idle component of each user on fees susceptance;

3) the idle cost of unit that each idle component on the fees susceptance be multiply by the reactive power source of this fees susceptance correspondence gets the idle cost component on the fees susceptance;

4) calculate each user with respect to meritorious and idle idle service electricity price.

2. the Wattless service pricing method based on fees susceptance according to claim 1, step 1) may further comprise the steps:

(i) the stable state trend of calculating electrical network is separated;

(ii) the user of i node is injected the active power P of electrical network _iAnd reactive power Q _iThe node injection current I that the corresponding stable state trend of equivalent respectively one-tenth is separated _i ^pAnd I _i ^q:

${\stackrel{\‾}{I}}_i^p=P_i/{\hat{V}}_i---\left(\mathrm{kA}\right)$

${\stackrel{\‾}{I}}_i^q={-\mathrm{jQ}}_i/{\hat{V}}_i---\left(\mathrm{kA}\right)$

Wherein

Be i node voltage V _iConjugate;

(iii) with the idle charge-collector's of k node reactive power generation power Q _kThe node fees susceptance that the corresponding stable state trend of equivalent one-tenth is separated:

$B_k=Q_k/{\stackrel{\‾}{V}}_k^2$ (S)

Wherein, V _kBe the k node voltage;

3. the Wattless service pricing method based on fees susceptance according to claim 2, step 2) may further comprise the steps:

(i) the meritorious equivalent electric current I of calculating user i _i ^pAt fees susceptance B _kOn the reactive power generation component Q that causes _Ik ^p:

$Q_{\mathrm{ik}}^p=B_k[\left({\stackrel{\‾}{V}}_k\right)\·\left({\stackrel{\‾}{Z}}_{\mathrm{ki}}{\stackrel{\‾}{I}}_i^p\right)]---\left(\mathrm{MVAR}\right)$

Wherein " " represents dot-product, Z _KiElement for the capable i row of grid nodes impedance matrix k;

(ii) calculate the idle equivalent electric current I of user i _i ^qAt fees susceptance B _kOn the reactive power generation component Q that causes _Ik ^q:

$Q_{\mathrm{ik}}^q=B_k[\left({\stackrel{\‾}{V}}_k\right)\·\left({\stackrel{\‾}{Z}}_{\mathrm{ki}}{\stackrel{\‾}{I}}_i^q\right)]---\left(\mathrm{MVAR}\right)$

Wherein " " represents dot-product, Z _KiElement for the capable i row of grid nodes impedance matrix k.

4. the Wattless service pricing method based on fees susceptance according to claim 3, step 3) may further comprise the steps:

(i) the meritorious equivalent electric current I of utilization user i _i ^pAt fees susceptance B _kOn the reactive power generation component Q that causes _Ik ^p, calculating corresponding reactive power generation cost component is γ _kQ _Ik ^p(unit); γ wherein _kThe idle cost of unit for k node reactive power source;

(ii) use the idle equivalent electric current I of user i _i ^qAt fees susceptance B _kOn the reactive power generation component Q that causes _Ik ^q, calculating corresponding reactive power generation cost component is γ _kQ _Ik ^q(unit); γ wherein _kThe idle cost of unit for k node reactive power source.

5. the Wattless service pricing method based on fees susceptance according to claim 4, described step 4) is: meritorious, the idle idle service electricity price λ that calculates transmission user i _i ^pAnd λ _i ^q:

${\mathrm{\λ}}_i^p=\underset{k}{\mathrm{\Σ}}{\mathrm{\γ}}_k{Q}_{\mathrm{ik}}^p/P_i=\underset{k}{\mathrm{\Σ}}{\mathrm{\γ}}_k{B}_k[\left({\stackrel{\‾}{V}}_k\right)\·({\stackrel{\‾}{Z}}_{\mathrm{ki}}/{\hat{V}}_i)]$ (unit/MWh)

${\mathrm{\λ}}_i^q=\underset{k}{\mathrm{\Σ}}{\mathrm{\γ}}_k{Q}_{\mathrm{ik}}^q/Q_i=\underset{k}{\mathrm{\Σ}}{\mathrm{\γ}}_k{B}_k[\left({\stackrel{\‾}{V}}_k\right)\·(-j{\stackrel{\‾}{Z}}_{\mathrm{ki}}/{\hat{V}}_i)]$ (unit/MVARh)

Images