CN110596455B - Power frequency electrical parameter extraction method, system and computer readable storage medium - Google Patents
Power frequency electrical parameter extraction method, system and computer readable storage medium Download PDFInfo
- Publication number
- CN110596455B CN110596455B CN201910905574.6A CN201910905574A CN110596455B CN 110596455 B CN110596455 B CN 110596455B CN 201910905574 A CN201910905574 A CN 201910905574A CN 110596455 B CN110596455 B CN 110596455B
- Authority
- CN
- China
- Prior art keywords
- phase
- axis
- voltage
- components
- current
- 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
Links
- 238000000605 extraction Methods 0.000 title claims abstract description 16
- 238000000034 method Methods 0.000 claims abstract description 35
- 230000008859 change Effects 0.000 claims abstract description 19
- 238000001914 filtration Methods 0.000 claims description 17
- 238000012545 processing Methods 0.000 claims description 13
- 238000006243 chemical reaction Methods 0.000 claims description 12
- 238000004590 computer program Methods 0.000 claims description 12
- 230000000737 periodic effect Effects 0.000 claims description 12
- 230000009466 transformation Effects 0.000 claims description 10
- 238000004364 calculation method Methods 0.000 claims description 9
- 230000008569 process Effects 0.000 claims description 8
- 230000001360 synchronised effect Effects 0.000 abstract description 5
- 238000010586 diagram Methods 0.000 description 10
- 230000006870 function Effects 0.000 description 5
- 238000001514 detection method Methods 0.000 description 3
- 230000009286 beneficial effect Effects 0.000 description 2
- 230000005540 biological transmission Effects 0.000 description 2
- 238000005070 sampling Methods 0.000 description 2
- 230000001052 transient effect Effects 0.000 description 2
- RZVAJINKPMORJF-UHFFFAOYSA-N Acetaminophen Chemical compound CC(=O)NC1=CC=C(O)C=C1 RZVAJINKPMORJF-UHFFFAOYSA-N 0.000 description 1
- 230000002159 abnormal effect Effects 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 230000002238 attenuated effect Effects 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000011084 recovery Methods 0.000 description 1
- 230000004044 response Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R19/00—Arrangements for measuring currents or voltages or for indicating presence or sign thereof
- G01R19/06—Measuring real component; Measuring reactive component
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R21/00—Arrangements for measuring electric power or power factor
- G01R21/127—Arrangements for measuring electric power or power factor by using pulse modulation
- G01R21/1271—Measuring real or reactive component, measuring apparent energy
Abstract
The invention discloses a power frequency electrical parameter extraction method, a system and a computer readable storage medium, comprising the following steps: acquiring three-phase voltage signals, three-phase current signals and phase information of three-phase voltage on an extra-high voltage direct current inversion end alternating current line natural coordinate system at the phase change failure moment; converting three-phase voltage signals and three-phase current signals on a natural coordinate system to a dq0 coordinate system based on a synchronous rotation coordinate change method to obtain voltage components and current components corresponding to d-axis and q-axis respectively; and then obtaining the direct current components corresponding to the d axis and the q axis, and calculating the apparent power, the active power and the reactive power corresponding to each axis by utilizing the direct current components corresponding to the d axis and the q axis. The method is convenient for rapidly extracting the power frequency component when the commutation fails.
Description
Technical Field
The invention belongs to the technical field of power grid data processing, and particularly relates to a power frequency electric parameter extraction method and system suitable for extra-high voltage direct current commutation failure and a computer readable storage medium.
Background
The distance between a large energy base and the load center in the middle east reaches 1000-. Therefore, China vigorously develops the ultra-high voltage alternating current and direct current hybrid power grid suitable for long-distance and large-capacity power transmission. However, with the rapid development of extra-high voltage alternating current and direct current, especially the rapid increase of the scale of extra-high voltage direct current transmission, the operation characteristics of the power grid are changed deeply. On one hand, the contradiction of 'strong direct current and weak alternating current' is prominent, and common disturbance or conventional operation such as alternating current line fault, reclosing, main transformer idle charging and the like can cause near-region direct current commutation failure, so that cascading failure of an alternating current and direct current system is caused, and huge transient energy impact is caused.
Meanwhile, the mutual influence characteristics of the alternating-current and direct-current hybrid power grid are obvious, and the dynamic response process of the direct-current control system in the fault period has discreteness and nonlinearity, so that the complexity and the variability of the fault transient process are further increased, and the simple fault global trend is presented to cause serious influence on the correct actions of the traditional power system relay protection device and the safety and stability control device. The extraction of power frequency of traditional power grid relay protection devices and various safety and stability control devices generally adopts full-wave or half-wave FFT (fast Fourier transform) technology, and the method at least needs 10mS or more to obtain correct power frequency electric parameters.
Disclosure of Invention
The invention aims to provide a power frequency electrical parameter extraction method, a system and a readable storage medium suitable for extra-high voltage direct current commutation failure, which are used for improving the extraction speed of power frequency point parameters, the method transforms abc three-phase sampling value data in a natural coordinate system to a dq0 coordinate system in real time based on synchronous rotating coordinate transformation (abc-dq0), then obtains direct current components of all phases in a rotating coordinate system, wherein the direct current components reflect the power frequency components in the abc coordinate system at the commutation failure moment, and then calculates the apparent power, the active power and the reactive power by utilizing the direct current components.
The invention provides a power frequency electrical parameter extraction method suitable for extra-high voltage direct current commutation failure, which comprises the following steps:
s1: acquiring three-phase voltage signals, three-phase current signals and phase information of three-phase voltage of an extra-high voltage direct current inversion end alternating current line on a natural coordinate system at the phase change failure moment;
s2: converting the three-phase voltage signals and the three-phase current signals on a natural coordinate system into a dq0 coordinate system by utilizing the phase information of the three-phase voltage to obtain voltage components and current components corresponding to d-axis and q-axis respectively;
s3: calculating any one or more of apparent power, active power and reactive power corresponding to each step by using the voltage components and the direct current components corresponding to the d axis and the q axis in the step S2;
and the direct current component of the d axis represents active components in the power frequency current and voltage, and the direct current component of the q axis represents reactive components in the power frequency current and voltage.
Further preferably, in step S3, each corresponding apparent power, active power, and reactive power is calculated according to the following formulas:
In the formula, S, P, Q pointsApparent power, active power and reactive power corresponding to other lambda,
lambda corresponds to the active component in the power frequency current and voltage in the direct current component of d-axis,
respectively represents the reactive components in the working frequency current and voltage in the direct current component of q axis corresponding to lambda, j is a complex sign, and the lambda phase is a phase a, a phase b or a phase c.
More preferably, in step S2, the three-phase voltage signals and the three-phase current signals are coordinate-converted according to the following formulas to obtain voltage components and current components corresponding to the d-axis and the q-axis:
In the formula iλp(t)、iλq(t) represents the currents i of the phases λλ(t) a d-axis current component and a q-axis current component after transformation to the dq0 coordinate system; u. ofλp(t)、uλq(t) each represents a voltage u of a lambda phaseλ(T) d-axis voltage component, q-axis voltage component, T after conversion to dq0 coordinate systemabc-dqIs a varying coordinate between the natural coordinate system and the dq0 coordinate systemλAnd phase information indicating a phase of a phase voltage, wherein the phase is a phase, b phase or c phase.
Preferably, the three-phase voltage signals and the three-phase current signals of the extra-high voltage direct current inverter end alternating current line at the phase conversion failure moment are expressed as follows:
in the formula iλ(t) is the current of lambda phase of the extra-high voltage DC inversion terminal AC line at the time of phase change failure, Iλ0Amplitude of non-periodic DC component corresponding to λ, IλpAnd IλqAmplitude values of power frequency cosine coordinate axis and power frequency sine coordinate axis components corresponding to the lambda are respectively set; i isλnThe amplitude of the nth higher harmonic component corresponding to λ; u. ofa(t) the voltage of lambda phase of the extra-high voltage DC inversion end AC line at the phase change failure moment, Uλ0Representing a non-periodic DC component, U, in the line voltage corresponding to λλpIndicating the component of the line voltage corresponding to λ that coincides with the fundamental active reference axis (coinciding with the phase instant of the phase-locked loop output, referred to in this patent as the active component), UλqRepresenting the component of the line voltage that is orthogonal to the fundamental active reference (kept orthogonal to the phase instants of the phase-locked loop output, referred to in this patent as the reactive component), UλnRepresenting the amplitude of the nth higher harmonic component in the line voltage corresponding to the lambda;
n represents the harmonic number, ω represents the angular frequency, t represents time,
indicating the initial phase of the nth harmonic component.
Further preferably, the voltage components and current components of d-axis and q-axis corresponding to the a-phase, b-phase, or c-phase are as follows:
further preferably, the process of acquiring the dc component is:
and filtering the voltage components of the d axis and the q axis, the alternating current components of the current components and the high-order harmonic components in the step S2 to obtain direct current components corresponding to the d axis and the q axis respectively.
Further preferably, the voltage components of the d-axis and q-axis, the alternating current components of the current components, and the harmonic components in step S2 are filtered out by using a low-pass filter.
In another aspect, the present invention provides a system based on the above method, including: the device comprises a signal acquisition module, a conversion module and a power calculation module;
the signal acquisition module is used for acquiring three-phase voltage signals, three-phase current signals and phase information of three-phase voltage of an extra-high voltage direct current inversion end alternating current line on a natural coordinate system at the phase change failure moment;
the conversion module is used for converting the three-phase voltage signals and the three-phase current signals on the natural coordinate system into a dq0 coordinate system by utilizing the phase information of the three-phase voltage to obtain voltage components and current components corresponding to the d axis and the q axis;
and the power calculation module is used for calculating any one or more of apparent power, active power and reactive power corresponding to each step by using the voltage components and the direct current components corresponding to the d axis and the q axis respectively.
Preferably, the system further comprises a filtering module, wherein the filtering module is used for filtering voltage components of the d axis and the q axis, alternating current components of current components and higher harmonic components to obtain direct current components corresponding to the d axis and the q axis.
Further preferably, the filtering module is a low-pass filter.
Furthermore, the present invention also provides a computer-readable storage medium comprising computer program instructions which, when executed by a processing terminal, cause the processing terminal to perform the above-mentioned method.
Advantageous effects
After the ultrahigh voltage direct current is accessed, once the receiving end has phase commutation failure or continuous phase commutation failure, the active and reactive fluctuation is very rapid in a power frequency period, the measuring method provided by the patent can effectively realize active power, reactive power and apparent power in the scene, and improve the working performance of various relay protection devices and safety and stability control devices in the ultrahigh voltage direct current receiving end converter station.
The invention combines synchronous rotation coordinate transformation and digital low-pass filtering, realizes the rapid extraction of direct current components for a single-phase system for the first time, and particularly provides a method for rapidly calculating power based on the direct current components, so that the instantaneous detection of active power and reactive power of fundamental wave signals is realized, the limitation that the time consumption is long when the existing FFT (full wave/half wave) transformation is used for obtaining the power frequency components of electric parameters is overcome, and the rapid extraction of the power frequency components is convenient for the failure of phase change.
Drawings
FIG. 1 is a schematic diagram of the implementation process of the method of the present invention in a single-phase power-frequency current component;
FIG. 2 is a schematic diagram of the implementation process of the method of the present invention in a single-phase power frequency voltage component;
fig. 3 is a comparison between the instantaneous electrical parameter extraction method and the conventional method, and curves 1, 2, and 3 correspond to A, B, C three phases respectively.
Detailed Description
The present invention will be further described with reference to the following examples.
When the extra-high voltage alternating current and direct current system has serious faults, the direct current system enters an abnormal operation state. From the perspective of harmonic components, the direct current system can be equivalent to a time-varying harmonic current source for an alternating current system, and a large amount of harmonics and inter-harmonics are injected into an adjacent alternating current line; from the view of power frequency components, equivalent power frequency current injected into the alternating current power grid by the direct current system is ahead of voltage, so that equivalent power frequency variable impedance sensed by the receiving-end alternating current power grid is capacitive; from the view of the non-periodic component, during the serious fault and fault recovery of the direct current system, the non-periodic component is increased rapidly, the attenuation time constant of the extra-high voltage system power grid is large, and the non-periodic component is attenuated slowly under some faults. The traditional electric parameter acquisition and calculation method based on the cycle signals can be stabilized only by 1-2 cycles, and is not beneficial to the rapid action of a protection and stability control device in the ultrahigh voltage direct current receiving end converter station. The power frequency electrical parameter extraction method suitable for the extra-high voltage direct current commutation failure realizes instantaneous power detection and is suitable for the requirements of extra-high voltage direct current inverter commutation failure time on various electrical parameter detection speeds.
As shown in fig. 1 and fig. 2, the specific implementation steps of the power frequency electrical parameter extraction method provided by the embodiment of the present invention, which is suitable for the extra-high voltage direct current commutation failure, are as follows:
step 1: and acquiring three-phase voltage signals, three-phase current signals and phase information of the three-phase voltage of the extra-high voltage direct current inverter end alternating current circuit at the phase change failure moment.
The embodiment of the invention can collect the three-phase voltage u of the alternating current circuit at the extra-high voltage direct current inversion end in real timea，ub，ucThree-phase current signal ia，ib，icAnd using a single-phase-locked loop to respectively align the three-phase voltages ua，ub，ucThe signal is subjected to phase tracking to obtain phase information theta of each phase voltage in real timea、θbAnd thetac. Once the commutation failure occurs, the corresponding phase voltage u is obtaineda，ub，ucThree-phase current signal ia，ib，icAnd phase information theta of three-phase voltagesa、θbAnd thetac。
In addition, the invention also provides phase information theta of three-phase voltagea、θbAnd thetacThe sine and cosine functions are respectively calculated and are expressed as: [ cos θ ]a sinθa]T，[cosθb sinθb]TAnd [ cos θ ]c sinθc]T。
Because the line current contains a large amount of non-periodic direct current components, alternating current power frequency components and alternating current higher harmonic components at the time of phase commutation failure, and similarly, the line voltage is also the same, the three-phase current and the three-phase voltage at the time of phase commutation failure can be expressed as follows (because the three phases a, b and c are similar, the embodiment of the invention uses λ to express the three phases a, b and c):
in the formula iλ(t) is the current of lambda phase of the extra-high voltage DC inversion terminal AC line at the time of phase change failure, Iλ0Amplitude of non-periodic DC component corresponding to λ, IλpAnd IλqAmplitude values of power frequency cosine coordinate axis and power frequency sine coordinate axis components corresponding to the lambda are respectively set; i isλnThe amplitude of the nth higher harmonic component corresponding to λ; u. ofa(t) the voltage of lambda phase of the extra-high voltage DC inversion end AC line at the phase change failure moment, Uλ0Representing a non-periodic DC component, U, in the line voltage corresponding to λλpIndicating the component of the line voltage corresponding to λ that coincides with the fundamental active reference axis (coinciding with the phase instant of the phase-locked loop output, referred to in this patent as the active component), UλqRepresenting the component of the line voltage that is orthogonal to the fundamental active reference (kept orthogonal to the phase instants of the phase-locked loop output, referred to in this patent as the reactive component), UλnRepresenting the amplitude of the nth higher harmonic component in the line voltage corresponding to the lambda; n represents the harmonic number, ω represents the angular frequency, t represents time,
indicating the initial phase of the nth harmonic component.
Based on equation 1, the embodiment of the present invention assumes phase a as an example, and the phase a current and the voltage at the phase commutation failure time are represented as:
step 2: converting the three-phase voltage signals and the three-phase current signals in the step 1 on the natural coordinate system into a dq0 coordinate system based on a synchronous rotating coordinate change method to obtain voltage components and current components corresponding to d-axis and q-axis respectively; the voltage components and the current components of the d axis and the q axis respectively comprise direct current components, alternating current components and higher harmonic components.
And performing coordinate conversion on each phase according to the following formula to obtain voltage components and current components corresponding to the d axis and the q axis:
in the formula iλp(t)、iλq(t) represents the currents i of the phases λλ(t) a d-axis current component and a q-axis current component after transformation to the dq0 coordinate system; u. ofλp(t)、uλq(t) each represents a voltage u of a lambda phaseλ(T) d-axis voltage component, q-axis voltage component, T after conversion to dq0 coordinate systemabc-dqIs a varying coordinate between the natural coordinate system and the dq0 coordinate systemλAnd phase information indicating a phase of a phase voltage, wherein the phase is a phase, b phase or c phase.
Similarly, taking phase a as an example, the voltage components and current components corresponding to d-axis and q-axis are:
wherein, taking d axis as an example, iap(t) the transformation process is as follows:
further transformed by equation (11):
for a phase i, the same appliesaq(t)、uap(t)、uaq(t) making a change to exist:
from the above formula, the component i on the d-axis after coordinate transformationap(t) contains, a direct current componentAlternating current component and higher harmonic component, wherein the amplitude of the direct current component is the amplitude I of the component on the working frequency cosine coordinate axis under the abc coordinate systemapThe amplitude I of the non-periodic direct current component of the alternating current component in the abc coordinate systema0And the amplitude I of the second harmonic componenta2The amplitudes of the higher harmonic components are constituted by the amplitudes of the alternating current components in all abc coordinate systems. Similarly, the component i in the q-axisaq(t) comprises a direct current component, an alternating current component and a higher harmonic component, wherein the amplitude of the direct current component is the amplitude I of the component on the working frequency sine coordinate axis under the abc coordinate systemaqThe amplitude I of the non-periodic direct current component of the alternating current component in the abc coordinate systema0And the amplitude I of the second harmonic componenta2The amplitudes of the higher harmonic components are constituted by the amplitudes of the alternating current components in all abc coordinate systems.
Similarly, the phase a is inferred to the phase b and the phase c, and the same method is adopted for variation respectively to obtain:
similarly, as can be seen from the above formula, the components on the d-axis and q-axis after coordinate transformation include a dc component, an ac component, and a harmonic component.
And step 3: and (3) processing and quickly filtering voltage components of the d axis and the q axis, alternating current components in current components and higher harmonic components in the step (2) by using a low-pass filter (LPF) to obtain direct current components corresponding to the d axis and the q axis. Since the voltage component and the current component of the d-axis and the q-axis are 50Hz, the cut-off frequency of the low-pass filter is selected to be not lower than 50 Hz.
For example, for the phase a, the reference phase adopted is the real-time phase of the phase voltage a, so the d-axis direct current component represents the active component in the power frequency current, and the q-axis represents the reactive component in the power frequency current, and therefore, the d-axis and q-axis direct current components after filtering can be expressed in a matrix form as follows:
by the same token, from phase a to phases b and c
Respectively showing the amplitude of active components in the power frequency current and voltage in the direct current component corresponding to the d axis after the filtering in the step 3,
respectively showing the amplitudes of the idle components in the working frequency current and voltage in the direct current component of the q axis corresponding to lambda after filtering in the step 3.
And 4, step 4: and 3, calculating the apparent power, the active power and the reactive power corresponding to each step by using the direct current components corresponding to each step on the d axis and the q axis.
Wherein, every corresponding apparent power, active power and reactive power are calculated according to the following formulas:
In the formula, S, P, Q λ respectively corresponds to the apparent power, the active power, and the reactive power, j is a complex symbol, such as phase a, and it can be known from the basic formula of the apparent power calculation that at any time of the sampling point, the formula of the power calculation is:
the description of the active power P and the reactive power Q of the line based on the synchronous rotating coordinate transformation is respectively as follows:
based on the method, the invention provides a power frequency electrical parameter extraction system, which comprises the following steps: the device comprises a signal acquisition module, a conversion module, a filtering module and a power calculation module;
the signal acquisition module is used for acquiring three-phase voltage signals, three-phase current signals and phase information of three-phase voltage of an extra-high voltage direct current inversion end alternating current line on a natural coordinate system at the phase change failure moment;
the conversion module is used for converting the three-phase voltage signals and the three-phase current signals on the natural coordinate system into a dq0 coordinate system by utilizing the phase information of the three-phase voltage to obtain voltage components and current components corresponding to the d axis and the q axis;
and the filtering module is used for filtering voltage components of the d axis and the q axis, alternating current components in the current components and higher harmonic components to obtain direct current components corresponding to the d axis and the q axis. In this embodiment, the filtering module selects a low-pass filter.
And the power calculation module is used for calculating any one or more of the apparent power, the active power and the reactive power corresponding to each direct current component by using the direct current components corresponding to each direct current component on the d axis and the q axis.
It should be understood that the functional unit modules in the embodiments of the present invention may be integrated into one processing unit, or each unit module may exist alone physically, or two or more unit modules are integrated into one unit module, and may be implemented in the form of hardware or software.
The embodiment of the invention also provides a computer readable storage medium, which comprises computer program instructions, and when the computer program instructions are executed by a processing terminal, the processing terminal executes the power frequency electrical parameter extraction method. The beneficial effects are referred to in the method part, and are not described in detail herein.
As will be appreciated by one skilled in the art, embodiments of the present application may be provided as a method, system, or computer program product. Accordingly, the present application may take the form of an entirely hardware embodiment, an entirely software embodiment or an embodiment combining software and hardware aspects. Furthermore, the present application may take the form of a computer program product embodied on one or more computer-usable storage media (including, but not limited to, disk storage, CD-ROM, optical storage, and the like) having computer-usable program code embodied therein.
The present application is described with reference to flowchart illustrations and/or block diagrams of methods, apparatus (systems), and computer program products according to embodiments of the application. It will be understood that each flow and/or block of the flow diagrams and/or block diagrams, and combinations of flows and/or blocks in the flow diagrams and/or block diagrams, can be implemented by computer program instructions. These computer program instructions may be provided to a processor of a general purpose computer, special purpose computer, embedded processor, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, create means for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be stored in a computer-readable memory that can direct a computer or other programmable data processing apparatus to function in a particular manner, such that the instructions stored in the computer-readable memory produce an article of manufacture including instruction means which implement the function specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be loaded onto a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer implemented process such that the instructions which execute on the computer or other programmable apparatus provide steps for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
According to the above description, the method for rapidly extracting the power frequency component based on the multi-synchronous rotation coordinate transformation and the digital low-pass filtering provided by the invention is suitable for rapidly decoupling and extracting the power frequency electrical parameters (voltage, current, active power and reactive power) at the phase commutation failure moment of the extra-high voltage direct current inverter.
Finally, it should be noted that: the above embodiments are only for illustrating the technical solutions of the present invention and not for limiting the same, and although the present invention is described in detail with reference to the above embodiments, those of ordinary skill in the art should understand that: modifications and equivalents may be made to the embodiments of the invention without departing from the spirit and scope of the invention, which is to be covered by the claims.
Claims (9)
1. A power frequency electrical parameter extraction method is characterized in that: the method comprises the following steps:
s1: acquiring three-phase voltage signals, three-phase current signals and phase information of three-phase voltage of an extra-high voltage direct current inversion end alternating current line on a natural coordinate system at the phase change failure moment;
s2: converting the three-phase voltage signals and the three-phase current signals on a natural coordinate system into a dq0 coordinate system by utilizing the phase information of the three-phase voltage to obtain voltage components and current components corresponding to d-axis and q-axis respectively;
s3: calculating any one or more of apparent power, active power and reactive power corresponding to each step by using the voltage components and the direct current components corresponding to the d axis and the q axis in the step S2;
the direct current component of the d axis represents active components in power frequency current and voltage, and the direct current component of the q axis represents reactive components in the power frequency current and voltage;
the method comprises the following steps of firstly, obtaining a three-phase voltage signal and a three-phase current signal of an extra-high voltage direct current inversion end alternating current circuit at the phase change failure moment, wherein the three-phase voltage signal and the three-phase current signal of the extra-high voltage direct current inversion end alternating current circuit at the phase change failure:
in the formula iλ(t) is the current of lambda phase of the extra-high voltage DC inversion terminal AC line at the time of phase change failure, Iλ0Amplitude of non-periodic DC component corresponding to λ, IλpAnd IλqAmplitude values of power frequency cosine coordinate axis and power frequency sine coordinate axis components corresponding to the lambda are respectively set; i isλnThe amplitude of the nth higher harmonic component corresponding to λ; u. ofa(t) the voltage of lambda phase of the extra-high voltage DC inversion end AC line at the phase change failure moment, Uλ0Representing a non-periodic DC component, U, in the line voltage corresponding to λλpRepresenting the component of the line voltage corresponding to λ that coincides with the fundamental active reference axis, UλqRepresenting the component of the line voltage that is orthogonal to the fundamental active reference, UλnRepresenting the amplitude of the nth higher harmonic component in the line voltage corresponding to the lambda;
n represents the harmonic number, ω represents the angular frequency, t represents time,
indicating the initial phase of the nth harmonic component.
2. The method of claim 1, wherein: in step S3, the apparent power, the active power, and the reactive power corresponding to each other are calculated according to the following formulas:
In the formula, S, P, Q represents the apparent power, the active power and the reactive power corresponding to lambda;
respectively represents the active components of working frequency current and voltage in the direct current component of d axis corresponding to lambda,
respectively represents the reactive components in the working frequency current and voltage in the direct current component of q axis corresponding to lambda, j is a complex sign, and the lambda phase is a phase a, a phase b or a phase c.
3. The method of claim 1, wherein: in step S2, coordinate conversion is performed on the three-phase voltage signal and the three-phase current signal according to the following formulas to obtain voltage components and current components corresponding to the d axis and the q axis:
In the formula iλp(t)、iλq(t) represents the currents i of the phases λλ(t) a d-axis current component and a q-axis current component after transformation to the dq0 coordinate system; u. ofλp(t)、uλq(t) each represents a voltage u of a lambda phaseλ(T) d-axis voltage component, q-axis voltage component, T after conversion to dq0 coordinate systemabc-dqIs a varying coordinate between the natural coordinate system and the dq0 coordinate systemλAnd phase information indicating a phase of a phase voltage, wherein the phase is a phase, b phase or c phase.
4. The method of claim 3, wherein: voltage components and current components of d-axis and q-axis corresponding to the phase a, the phase b or the phase c are shown as follows:
5. the method of claim 1, wherein: the process of acquiring the direct current component comprises the following steps:
and filtering the voltage components of the d axis and the q axis, the alternating current components of the current components and the high-order harmonic components in the step S2 to obtain direct current components corresponding to the d axis and the q axis respectively.
6. The method of claim 5, wherein: in step S2, the voltage components of the d axis and the q axis, the alternating current components of the current components, and the higher harmonic components are filtered out by using a low-pass filter.
7. A system based on the method of any one of claims 1-6, characterized by: the method comprises the following steps: the device comprises a signal acquisition module, a conversion module and a power calculation module;
the signal acquisition module is used for acquiring three-phase voltage signals, three-phase current signals and phase information of three-phase voltage of an extra-high voltage direct current inversion end alternating current line on a natural coordinate system at the phase change failure moment;
the conversion module is used for converting the three-phase voltage signals and the three-phase current signals on the natural coordinate system into a dq0 coordinate system by utilizing the phase information of the three-phase voltage to obtain voltage components and current components corresponding to the d axis and the q axis;
and the power calculation module is used for calculating any one or more of apparent power, active power and reactive power corresponding to each step by using the voltage components and the direct current components corresponding to the d axis and the q axis respectively.
8. The system of claim 7, wherein: the filter module is used for filtering voltage components of the d axis and the q axis, alternating current components in current components and higher harmonic components to obtain direct current components corresponding to the d axis and the q axis.
9. A computer readable storage medium comprising computer program instructions characterized in that: the computer program instructions, when executed by a processing terminal, cause the processing terminal to perform the method of any of claims 1 to 6.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910905574.6A CN110596455B (en) | 2019-09-24 | 2019-09-24 | Power frequency electrical parameter extraction method, system and computer readable storage medium |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910905574.6A CN110596455B (en) | 2019-09-24 | 2019-09-24 | Power frequency electrical parameter extraction method, system and computer readable storage medium |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110596455A CN110596455A (en) | 2019-12-20 |
| CN110596455B true CN110596455B (en) | 2021-07-13 |
Family
ID=68862816
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910905574.6A Active CN110596455B (en) | 2019-09-24 | 2019-09-24 | Power frequency electrical parameter extraction method, system and computer readable storage medium |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110596455B (en) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112600470B (en) * | 2020-12-16 | 2022-07-05 | 华北电力大学 | Power tracking excitation control method and system for double-shaft excitation generator |
| CN116223900B (en) * | 2022-12-30 | 2024-01-09 | 广州视骁科技有限公司 | Method, device, electronic equipment and storage medium for determining alternating current circuit power |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102222922A (en) * | 2011-06-10 | 2011-10-19 | 湖南大学 | STATCOM (static synchronous compensator) control system adopting droop control strategy and control method thereof |
| CN102957153A (en) * | 2012-07-31 | 2013-03-06 | 上海交通大学 | Direct-current neutral point balance controlling method of three-phase four-wire system three-bridge arm static synchronous compensator |
| CN103050988A (en) * | 2013-01-21 | 2013-04-17 | 四川极度电控系统制造有限责任公司 | Design method of converter station controller of flexible direct-current transmission system |
| JP2013192398A (en) * | 2012-03-14 | 2013-09-26 | Nissan Motor Co Ltd | Controller of motor and control method of motor |
| KR20170037748A (en) * | 2015-09-25 | 2017-04-05 | 전남대학교산학협력단 | Instantaneous reactive power and active power measurement techniques using all pass filter |
-
2019
- 2019-09-24 CN CN201910905574.6A patent/CN110596455B/en active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102222922A (en) * | 2011-06-10 | 2011-10-19 | 湖南大学 | STATCOM (static synchronous compensator) control system adopting droop control strategy and control method thereof |
| JP2013192398A (en) * | 2012-03-14 | 2013-09-26 | Nissan Motor Co Ltd | Controller of motor and control method of motor |
| CN102957153A (en) * | 2012-07-31 | 2013-03-06 | 上海交通大学 | Direct-current neutral point balance controlling method of three-phase four-wire system three-bridge arm static synchronous compensator |
| CN103050988A (en) * | 2013-01-21 | 2013-04-17 | 四川极度电控系统制造有限责任公司 | Design method of converter station controller of flexible direct-current transmission system |
| KR20170037748A (en) * | 2015-09-25 | 2017-04-05 | 전남대학교산학협력단 | Instantaneous reactive power and active power measurement techniques using all pass filter |
Non-Patent Citations (1)
| Title |
|---|
| "混合多馈入HVDC换相失败时控制与保护策略的研究";刘泽锴;《中国优秀硕士学位论文全文数据库工程科技Ⅱ辑》;20160515;第18-20页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN110596455A (en) | 2019-12-20 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN101673952B (en) | Precise phase locking method based on cross decoupling self-adaptive complex filter | |
| CN108599165B (en) | Three-level active filter based on composite control strategy, application system and method thereof | |
| CN107632200B (en) | Harmonic detection method based on discrete Fourier transform | |
| CN109412191B (en) | Phase locking method, device and equipment for high-voltage direct-current power transmission system | |
| CN104333244B (en) | Positive sequence component-based three-phase inverter control method and device | |
| CN104659813B (en) | A kind of multi-inverter parallel control method of quick harmonic circulating current suppression | |
| CN106655277B (en) | Improved phase-locked loop method for permanent magnet synchronous motor | |
| CN110596455B (en) | Power frequency electrical parameter extraction method, system and computer readable storage medium | |
| CN106410858A (en) | Software digital phase-locking method based on dual dq coordination conversion | |
| CN104502705A (en) | Non-phase-locked-loop rotating vector detection method suitable for power grid voltage distortion and imbalance | |
| CN104597367A (en) | Transducer drive induction motor stator turn-to-turn short circuit fault diagnosis method | |
| CN109600064B (en) | Method and system for calculating alternating current asymmetric fault main loop of modular multilevel converter | |
| CN105071390A (en) | Control method of H-bridge three-level active power filter and system | |
| CN104393809B (en) | Pumped storage group low-speed position detection method applicable to SCR static frequency converter | |
| CN113241748B (en) | Method and system for suppressing transient overvoltage of power electronic converter accessing weak power grid | |
| CN103117562A (en) | Control method of high-voltage cascade energy feedback converter power module | |
| CN110323767B (en) | Power distribution station area unbalance management method, system and storage medium | |
| CN102623996B (en) | Active power filter closed loop control method based on decoupling resonance regulator array | |
| CN104280593A (en) | Quick harmonic and reactive current detection method based on instantaneous reactive power theory | |
| CN107994565B (en) | Simulation method and system of unified power flow controller | |
| CN113541185B (en) | Dq conversion wind power converter grid voltage fault detection method | |
| CN112904073B (en) | Method for estimating capacitance value of bus capacitor of driving system of permanent magnet compressor without electrolytic capacitor | |
| CN106771507B (en) | Reactive current rapid detection method based on voltage reference split-phase synchronization | |
| CN113625066B (en) | Distribution transformer phase unbalance detection method, system, device and storage medium | |
| JP2013116019A (en) | Independent operation detection device, system interconnection inverter system, and independent operation detection method |
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 |