CN112904133A - Step performance testing method and system of direct current control protection system - Google Patents

Step performance testing method and system of direct current control protection system Download PDF

Info

Publication number
CN112904133A
CN112904133A CN202110183805.4A CN202110183805A CN112904133A CN 112904133 A CN112904133 A CN 112904133A CN 202110183805 A CN202110183805 A CN 202110183805A CN 112904133 A CN112904133 A CN 112904133A
Authority
CN
China
Prior art keywords
value
time
point
curve
direct 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.)
Granted
Application number
CN202110183805.4A
Other languages
Chinese (zh)
Other versions
CN112904133B (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.)
China South Power Grid International Co ltd
Electric Power Dispatch Control Center of Guangdong Power Grid Co Ltd
Original Assignee
China South Power Grid International Co ltd
Electric Power Dispatch Control Center of Guangdong Power Grid 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 China South Power Grid International Co ltd, Electric Power Dispatch Control Center of Guangdong Power Grid Co Ltd filed Critical China South Power Grid International Co ltd
Priority to CN202110183805.4A priority Critical patent/CN112904133B/en
Publication of CN112904133A publication Critical patent/CN112904133A/en
Application granted granted Critical
Publication of CN112904133B publication Critical patent/CN112904133B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R31/00—Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
    • 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
    • H02J1/00—Circuit arrangements for dc mains or dc 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
    • H02J3/00—Circuit arrangements for ac mains or ac distribution networks
    • H02J3/36—Arrangements for transfer of electric power between ac networks via a high-tension dc link
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E60/00—Enabling technologies; Technologies with a potential or indirect contribution to GHG emissions mitigation
    • Y02E60/60—Arrangements for transfer of electric power between AC networks or generators via a high voltage DC link [HVCD]

Landscapes

  • Engineering & Computer Science (AREA)
  • Power Engineering (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Measurement Of Current Or Voltage (AREA)
  • Emergency Protection Circuit Devices (AREA)

Abstract

The application discloses a step performance test method and a system of a direct current control protection system, firstly, a method of combining a mathematical morphology filter with a sliding average filter is provided, harmonic waves and interference of direct current voltage are filtered, the complete inhibition of the harmonic waves and the interference during step response of the electrical quantity of the direct current voltage in a direct current transmission system is realized, meanwhile, deviation of measurement time is not caused, a smooth waveform trend curve of the direct current voltage is further obtained, a foundation is laid for finding out the overshoot points of the electrical quantity such as the direct current voltage of the direct current transmission system and the like, further, response time is identified, then, the overshoot points are identified according to the direct current voltage curve without the harmonic waves and the interference, whether the curve has the overshoot points or not is judged, finally, the response time with the overshoot points and the response time without the overshoot points are respectively calculated, and the automatic batch judgment of the direct current transmission system direct current voltage step response test result of a computer is realized, greatly improving the efficiency and accuracy.

Description

Step performance testing method and system of direct current control protection system
Technical Field
The application relates to the technical field of electric power, in particular to a step performance testing method and system of a direct current control protection system.
Background
The core of a direct current transmission system is a control protection system, in order to test the dynamic response performance of direct current transmission disturbed, an obvious step impact is usually given to direct current, voltage and power instruction values, the dynamic response characteristics of the direct current system after the direct current transmission disturbed are observed, the two key indexes comprise recovery time and overshoot, and whether the two indexes meet the design requirements or not is determined.
At present, the response time and overshoot of the step response are strictly defined, and the actual waveform has harmonic waves and interference due to the influence of sampling and system operating characteristics, so that the step response time cannot be intelligently identified through a computer. Therefore, the conventional test for the dynamic response performance of the disturbed direct current transmission can only judge the step response time and the overshoot by manual observation, but has limited precision, complex test and large workload.
Disclosure of Invention
The embodiment of the application provides a method and a system for testing the step performance of a direct current control protection system, which are used for solving the technical problems of large error and large workload caused by the fact that the judgment of step response time and overshoot in direct current transmission can only be realized through manual observation in the prior art.
In view of this, a first aspect of the present application provides a method for testing step performance of a dc control protection system, where the method includes:
s1, presetting the original direct current voltage as u (t), and t e [ t ∈1，t2]Wherein t is1＜t2Defining the filtering length delta t of the mathematical form filter of the direct-current voltage, and establishing a mathematical form filter processing function M (u) of the direct-current voltage according to a maximum value processing function, a minimum value processing function and a time offset function in the filtering length;
s2, after setting the offset time, the filtering length and the filtering frequency value of the time offset function, establishing a sliding average filtering function N (u) of the direct current voltage according to the offset time, the filtering length and the filtering frequency value;
s3, establishing a filter function y (u) of the dc voltage according to the moving average filter function and the mathematical morphology filter processing function, wherein y (u) is N (M (u (t));
s4, acquiring an initial value vb of the curve according to a curve starting point tb, a step starting point t0 and a curve ending point te of the curve of the filter function, and a first time td1, a second time td2 and a third time td3 corresponding to a first ending value ve, a second ending value ve90 and a third ending value ve 110;
s5, judging whether the curve has overshoot points according to the first end value ve, the second end value ve90 and the third end value ve110, if yes, calculating response time according to the step starting point and the fourth end value ve130, and if not, calculating the response time according to the step starting point and the second end value.
Optionally, step S1 specifically includes:
s01, presetting the original direct current voltage as u (t), and t e [ t ∈1，t2]Wherein t is1＜t2And defining the filter length of the mathematical form filter of the direct current voltage as delta t, and satisfying that delta t is less than t2-t1And the time offset function is O (u, t)o) Wherein, toIs the offset time;
s02, respectively establishing a maximum value processing function f (u) and a minimum value processing function g (u) of the mathematical morphology filter according to the maximum value and the minimum value of the direct current voltage in the filtering length;
s03, establishing the mathematical morphology filter processing function as M (u) according to the maximum processing function, the minimum processing function and the time offset function;
wherein,
f(u)＝f(u(t))＝max[t-Δt，t]u (t), and finally t e [ t ]1+Δt，t2]；
g(u)＝g(u(t))＝min[t-Δt，t]u (t), and finally t e [ t ]1+Δt，t2]；
M(u)=M(u(t))=0.5*O(f(g(g(f(u(t)))))+g(f(f(g(u(t)))))，to) Finally t ∈ [ t ]1+4Δt-to，t2-to]。
Optionally, step S4 specifically includes:
s04, defining the curve starting point of the filter function curve of the direct current voltage as tb, the step starting point as t0, the first overshoot point as t1, the curve ending point as te, the initial value as vb, the first ending value as ve, the second ending value as ve90, the third ending value as ve110 and the fourth ending value as ve 130;
s05, setting an average value of a first interval between the curve starting point and the step starting point as the initial value, where the first interval is [ tb, tb +0.02 × N ], and N is a maximum positive integer such that tb +0.02 × N < t 0;
s06, setting the average value of a second interval as the curve termination point, wherein the second interval is [ te-0.02, te ];
s07, regarding a time continuously lower than the second end value as the first time td1, the second end value is ve90 ═ ve-0.1(ve-vb), regarding a time corresponding to a first time greater than the third end value as the second time td2, the third end value ve110 ═ ve +0.1(ve-vb), and regarding a time corresponding to a first time smaller than the second end value as the third time td3, in a direction from the end point of the curve to the start point of the curve.
Optionally, step S5 specifically includes:
s08, traversing the derivative value of each point on the curve to the right by taking the first moment as a starting point, judging whether the derivative value is not larger than the first point corresponding to a zero value, if so, executing a step S09, and if not, setting the response time to tr as td1-t 0;
s09, using the first point as the first overshoot point, when the value of the first overshoot point is greater than the fourth end value, setting the response time tr to max (td2, td3) -t0, when the value of the first overshoot point is less than the fourth end value, traversing the corresponding value on the curve to the right, setting the corresponding time of the first point which is less than the third end value as a fourth time td4, when the fourth time is equal to the second time, setting the response time tr to td3-t0, and when the fourth time is less than the second time, setting the response time tr to max (td2, td3) -t 0.
Optionally, the time offset function is:
O(u，to)＝O(u(t)，to)＝u(t-to) Finally t ∈ [ t ]1-to，t2-to]；
In the formula, toIs offset by a time to the left.
Optionally, the moving average filter function n (u) is:
Figure BDA0002942813610000031
a second aspect of the present application provides a system for testing step performance of a dc control protection system, the system comprising:
a first modeling unit for presetting an original DC voltage tou(t)，t∈[t1，t2]Wherein t is1＜t2Defining the filtering length of a mathematical form filter of the direct current voltage, and establishing a mathematical form filter processing function M (u) of the direct current voltage according to a maximum value processing function, a minimum value processing function and a time offset function in the filtering length;
the second modeling unit is used for establishing a sliding average filtering function N (u) of the direct-current voltage according to the offset time, the filtering length and the filtering frequency value after setting the offset time, the filtering length and the filtering frequency value of the time offset function;
a third modeling unit, configured to establish a filter function y (u) of the dc voltage according to the moving average filter function and the mathematical morphology filter processing function, where y (u) ═ N (M (u (t));
the overshoot point analysis unit is used for acquiring an initial value vb of the curve according to a curve starting point tb, a step starting point t0 and a curve ending point te of the curve of the filter function, and a first time td1, a second time td2 and a third time td3 which correspond to the first ending value ve, the second ending value ve90 and the third ending value ve 110;
and the response time calculating unit is used for judging whether the curve has an overshoot point according to the first termination value ve, the second termination value ve90 and the third termination value ve110, if so, calculating the response time according to the step starting point and the fourth termination value ve130, and otherwise, calculating the response time according to the step starting point and the second termination value.
Optionally, the first modeling unit is specifically configured to:
presetting the original direct current voltage as u (t), t e [ t ∈ [ [ t ]1，t2]Wherein t is1＜t2And defining the filter length of the mathematical form filter of the direct current voltage as delta t < t2-t1And the time offset function is O (u, t)o) Wherein, toIs the offset time;
respectively establishing a maximum value processing function f (u) and a minimum value processing function g (u) of the mathematical form filter according to the maximum value and the minimum value of the direct-current voltage in the filtering length;
establishing a mathematical form filter processing function M (u) according to the maximum processing function, the minimum processing function and the time offset function;
wherein,
f(u)＝f(u(t))＝max[t-Δt，t]u (t), and finally t e [ t ]1+Δt，t2]；
g(u)＝g(u(t))＝min[t-Δt，t]u (t), and finally t e [ t ]1+Δt，t2]；
M(u)=M(u(t))=0.5*O(f(g(g(f(u(t)))))+g(f(f(g(u(t)))))，to) Finally t ∈ [ t ]1+4Δt-to，t2-to]。
Optionally, the overshoot point analysis unit is specifically configured to:
defining the curve starting point of a filter function curve of the direct-current voltage as tb, the step starting point as t0, the first overshoot point as t1, the curve ending point as te, the initial value as vb, the first ending value as ve, the second ending value as ve90, the third ending value as ve110 and the fourth ending value as ve 130;
setting an average value of a first interval between the curve starting point and the step starting point as the initial value, wherein the first interval is [ tb, tb +0.02 × N ], and N is a maximum positive integer such that tb +0.02 × N < t 0;
setting the average value of a second interval as the curve termination point, wherein the second interval is [ te-0.02, te ];
in the direction from the curve ending point to the curve starting point, a time continuously lower than the second ending value is taken as the first time td1, the second ending value is ve90 ═ ve-0.1(ve-vb), a time corresponding to the first time being greater than the third ending value is taken as the second time td2, the third ending value is ve110 ═ ve +0.1(ve-vb), and a time corresponding to the first time being less than the second ending value is taken as the third time td 3.
Optionally, the response time calculating unit is specifically configured to:
traversing the derivative value of each point on the curve to the right by taking the first moment as a starting point, judging whether the derivative value is not greater than the first point corresponding to a zero value, if so, executing a step S09, otherwise, setting the response time to tr to td1-t 0;
and taking the first point as the first overshoot point, when the value of the first overshoot point is greater than the fourth end value, setting the response time tr to be max (td2, td3) -t0, when the value of the first overshoot point is less than the fourth end value, traversing the value corresponding to each point on the curve to the right, setting the time corresponding to the first point which is less than the third end value as a fourth time td4, when the fourth time is equal to the second time, setting the response time tr to be td3-t0, and when the fourth time is less than the second time, setting the response time tr to be max (td2, td3) -t 0.
According to the technical scheme, the embodiment of the application has the following advantages:
the embodiment of the application provides a step performance test method of a direct current control protection system, firstly, a method of combining a mathematical morphology filter with a sliding average filter is provided, harmonic waves and interference of direct current voltage are filtered, complete suppression of the harmonic waves and the interference during step response of the electrical quantity of the direct current voltage in a direct current transmission system is realized, meanwhile, deviation of measurement time is not caused, a smooth waveform trend curve of the direct current voltage is further obtained, a foundation is laid for finding out overshoot points of the electrical quantity such as the direct current voltage of the direct current transmission system and the like, further response time is identified, then overshoot points are identified according to the direct current voltage curve without the harmonic waves and the interference, whether the curve has the overshoot points or not is judged, finally, response time with the overshoot points and without the overshoot points is respectively calculated, and the result of a computer automatically judging the step response test of the direct current voltage of the direct current transmission system in batches is realized, greatly improving the efficiency and accuracy. Therefore, the technical problems that the error is large and the workload is large due to the fact that the judgment of the step response time and the overshoot in the direct current transmission can only be achieved through manual observation in the prior testing technology are solved.
Drawings
Fig. 1 is a schematic flowchart of a first embodiment of a method for testing step performance of a dc control protection system provided in an embodiment of the present application;
fig. 2 is a schematic flowchart of a second embodiment of a step performance testing method for a dc control protection system provided in the embodiment of the present application;
fig. 3 is a schematic structural diagram of an embodiment of a step performance testing system of a dc control protection system provided in the embodiment of the present application;
FIGS. 4a, 4b and 4c are schematic diagrams of the definition for the last step response;
FIG. 5 is a schematic diagram of a waveform of an extreme 1 DC voltage in a manually triggered DC voltage step in a certain extra-high voltage DC system;
fig. 6 is a schematic waveform diagram of the dc voltage signal after performing sliding filtering and left shifting according to the present application.
Detailed Description
In order to make the technical solutions of the present application better understood, the technical solutions in the embodiments of the present application will be clearly and completely described below with reference to the drawings in the embodiments of the present application, and it is obvious that the described embodiments are only a part of the embodiments of the present application, and not all of the embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present application.
Referring to fig. 4a, 4b and 4c, the upper step response definition is as follows:
(1) if the measured value can be maintained within a range of ± 10% of the variation of the setting value with the error being the error with the new setting value as the center line after the first overshoot, the response time 'trl' can be defined as the time required from the moment of the variation of the setting value until the measured value reaches 90% of the variation of the setting value.
The conditions that should be satisfied with 'trl' are: the first overshoot does not exceed 30% of the variation of the setting value, and the measurement value is finally stabilized to the new setting value after the second overshoot.
(2) If the measured value exceeds the range of ± 10% of the setting value variation amount again after the first overshoot, the response time 'tr 2' is defined as: and the time from the moment of change of the setting value to the time required for the measured value to enter and be kept within the range of +/-10% of the variation of the setting value by taking the new setting value as a central line and the error.
(3) If there is no overshoot, i.e., an over-damped system, the response time 'tr 3' is defined as: the time required from the moment of the change of the setting value to the time when the measured value reaches 90% of the change amount of the setting value.
The lower step is similar to the upper step, and the only difference is that the step direction is opposite, which is not described herein again.
The following is a description of the prior art identification of step responses by human or computer:
referring to fig. 5, fig. 5 is a waveform of the pole 1 dc voltage in an extra-high voltage dc system when the dc voltage step is triggered manually.
As can be seen from fig. 5, the dc voltage is finally stabilized at 420kV after being overshot for about 2 times from 0V, and the step response time needs to be calculated in order to determine whether the step response performance of the dc control system meets the design requirement. The common method is manual identification, and human eyes can directly and quickly find the real moment of the first overshoot point based on the channel 1 direct current voltage, but the error is large. If a computer is adopted for automatic identification, a complete method is not available at present, and only a basic idea is provided, namely, a low-pass filter (such as moving average filtering) is adopted for filtering to obtain a direct-current voltage filtered by a channel 2, and then identification is carried out, but the method has great difficulty in actual operation and is embodied as follows:
1. a second overshoot point larger than the first overshoot point may exist, so that the global maximum value cannot be directly found;
2. because the amplified channel 2 curve always has micro interference and fluctuation and a large number of local maximum values, a real first overshoot point cannot be found. And because of the principle of low-pass filtering, these local maxima can't be eliminated with the increase of the filtering intensity, exist all the time, and simultaneously increasing the filtering intensity will lead to the curve delay after filtering to increase, influence the moment identification of overshoot point, lead to unable from the step starting point to look for the first maximum.
Referring to fig. 1, a first embodiment of a method for testing step performance of a dc control protection system provided in the embodiment of the present application includes:
step 101, presetting an original direct current voltage as u (t), and t e [ t ∈ [ [ t ]1，t2]Wherein t is1＜t2Defining the filter length delta t of the mathematical form filter of the direct current voltage, and establishing a mathematical form filter processing function M (u) of the direct current voltage according to a maximum value processing function, a minimum value processing function and a time offset function in the filter length.
Step 102, after setting the offset time, the filter length and the filter frequency value of the time offset function, establishing a sliding average filter function n (u) of the direct current voltage according to the offset time, the filter length and the filter frequency value.
Step 103, a filter function y (u) of the direct current voltage is established according to the moving average filter function and the mathematical form filter processing function, wherein y (u) is equal to N (M (u (t)).
Step 104, obtaining an initial value vb of the curve according to a curve starting point tb, a step starting point t0 and a curve ending point te of the curve of the filter function, and a first time td1, a second time td2 and a third time td3 corresponding to the first ending value ve, the second ending value ve90 and the third ending value ve 110.
And 105, judging whether the curve has an overshoot point according to the first termination value ve, the second termination value ve90 and the third termination value ve110, if so, calculating response time according to the step starting point and the fourth termination value ve130, and otherwise, calculating response time according to the step starting point and the second termination value.
The first embodiment of the application provides a step performance testing method of a direct current control protection system, which firstly provides a method of combining a mathematical morphology filter with a sliding average filter to filter harmonic waves and interference of direct current voltage, so as to completely inhibit the harmonic waves and the interference when the electrical quantity of the direct current voltage in a direct current transmission system responds to the step, and simultaneously, the deviation of measurement time is not caused, and further, a smooth waveform trend curve of the direct current voltage is obtained, a foundation is laid for finding out the overshoot points of the electrical quantities such as the direct current voltage of the direct current transmission system and the like, and further identifying the response time, then, the overshoot points are identified according to the direct current voltage curve without the harmonic waves and the interference, whether the curve has the overshoot points is judged, finally, the response time with the overshoot points and the response time without the overshoot points are respectively calculated, and the result of the step response test of the direct current voltage of the direct current transmission system is automatically judged in batches by a computer, greatly improving the efficiency and accuracy. Therefore, the technical problems that the error is large and the workload is large due to the fact that the judgment of the step response time and the overshoot in the direct current transmission can only be achieved through manual observation in the prior testing technology are solved.
The above is a first embodiment of a method for testing the step performance of the dc control protection system provided in the embodiment of the present application, and the following is a second embodiment of a method for testing the step performance of the dc control protection system provided in the embodiment of the present application.
Referring to fig. 2, a second embodiment of a method for testing step performance of a dc control protection system provided in the embodiment of the present application includes:
step 201, presetting an original direct current voltage as u (t), and t e [ t ∈ [ [ t ]1，t2]Wherein t is1＜t2And the filter length of the mathematical form filter of the direct current voltage is defined as delta t, and the condition that delta t is less than t is satisfied2-t1And the time offset function is O (u, t)o) Wherein, toIs the offset time.
Wherein the time-left shift function is:
O(u，to)＝O(u(t)，to)＝u(t-to) Finally t ∈ [ t ]1-to，t2-to]；
In the formula, toIs the offset time.
Step 202, respectively establishing a maximum value processing function f (u) and a minimum value processing function g (u) of the mathematical form filter according to the maximum value and the minimum value of the direct current voltage in the filtering length.
It should be noted that the maximum value can be defined as max[t-Δt，t]u(t)，max[t-Δt，t]u (t) represents [ t- Δ t, t]The maximum value of u (t) in the interval, and the minimum value can be defined as min[t-Δt，t]u(t)，min[t-Δt，t]u (t) represents [ t- Δ t, t]Minimum value of u (t) within the interval
Step 203, establishing a mathematical form filter processing function M (u) according to the maximum value processing function, the minimum value processing function and the time offset function;
wherein,
f(u)＝f(u(t))＝max[t-Δt，t]u (t), and finally t e [ t ]1+Δt，t2]；
g(u)＝g(u(t))＝min[t-Δt，t]u (t), and finally t e [ t ]1+Δt，t2]；
M(u)=M(u(t))=0.5*O(f(g(g(f(u(t)))))+g(f(f(g(u(t)))))，to) Finally t ∈ [ t ]1+4Δt-to，t2-to]。
Step 204, after setting the offset time, the filter length and the filter frequency value of the time offset function, establishing a sliding average filter function n (u) of the dc voltage according to the offset time, the filter length and the filter frequency value.
Wherein the moving average filter function n (u) is:
Figure BDA0002942813610000091
it should be noted that in the embodiment of the present application, Δ t is preferably 0.02, and the offset time t is preferably seto0.02 was taken. The filtering is carried out through the sliding average filtering function, all harmonic waves and interference higher than 50Hz of a curve after filtering are eliminated, the countless small maximum value of the direct current voltage waveform is eliminated, the shape of countless steps connected by a plurality of transverse lines is changed, the vicinity of an overshoot value is flattened, in order to find the maximum point, the filtering is carried out by adopting the sliding average filtering of 20ms and then leftwards moving by 10ms to offset time delay, and the final filtering curve is formedLine y (u).
Step 205, a filter function y (u) of the dc voltage is established according to the moving average filter function and the mathematical morphology filter processing function, where y (u) is equal to N (M (u (t)).
Please refer to fig. 6, fig. 6 is a schematic diagram of a waveform of a dc voltage signal after performing sliding filtering and left shifting according to the present application, and it should be noted that a curve of the original dc voltage signal after performing sliding filtering and left shifting according to a filtering function of the dc voltage established in the present application becomes smoother, so that harmonics and interferences of the dc voltage are filtered, and the harmonics and interferences are completely suppressed when an electrical quantity step response of the dc voltage in the dc power transmission system is realized, thereby obtaining a smooth waveform trend curve of the dc voltage.
Step 206, defining a curve starting point tb, a step starting point t0, a first overshoot point t1, a curve ending point te, an initial value vb, a first ending value ve, a second ending value ve90, a third ending value ve110, and a fourth ending value ve130 of the filter function curve of the direct current voltage.
It should be noted that the variables defined in this application: the curve starting point tb, the step starting point t0 and the curve ending point te are all known quantities, and the first overshoot point t1, the initial value vb, the second ending value ve90, the third ending value ve110 and the fourth ending value ve130 are all unknown quantities.
For the definitions of the second end value ve90, the third end value ve110, and the fourth end value ve130, the skilled person in the application may set the values according to practical situations, for example, the fourth end value ve120, and the fourth end value ve130 is expressed by the above definition of the last step: the first overshoot must not exceed 30% of the setting value variation.
Similarly, for the analysis of the step-down response, the skilled person can set the analysis according to the specific situation, and the step-down and step-up principles are similar.
Step 207, setting an average value of a first interval between the curve starting point and the step starting point as an initial value, wherein the first interval is [ tb, tb +0.02 × N ], and N is a maximum positive integer which enables tb +0.02 × N to be less than t 0;
and step 208, setting the average value of a second interval as a curve termination point, wherein the second interval is [ te-0.02, te ].
In step 209, in the direction from the curve ending point to the curve starting point, the time continuously lower than the second ending value is taken as the first time td1, the second ending value is ve90 ═ ve-0.1(ve-vb), the time corresponding to the first time being greater than the third ending value is taken as the second time td2, the third ending value ve110 is ve +0.1(ve-vb), and the time corresponding to the first time being less than the second ending value is taken as the third time td 3.
And step 210, traversing the derivative value of each point on the curve to the right by taking the first moment as a starting point, judging whether the derivative value is not greater than the first point corresponding to a zero value, if so, executing step S09, and if not, setting the response time to tr to td1-t 0.
And step 211, taking the first point as a first overshoot point, when the value of the first overshoot point is greater than a fourth end value, the response time tr is max (td2, td3) -t0, when the value of the first overshoot point is less than the fourth end value, traversing the value corresponding to each point on the curve to the right, setting the time corresponding to the first point which is less than the third end value as a fourth time td4, when the fourth time is equal to the second time, the response time tr is td3-t0, and when the fourth time is less than the second time, the response time tr is max (td2, td3) -t 0.
It can be understood that the core algorithm of the step performance testing method of the direct current control protection system provided by the application is mathematical morphology filtering, sliding average filtering and curve forward moving algorithm, harmonic waves and interference are filtered, the time of the filtered curve is consistent with that of the original curve, and the first overshoot point is found out to calculate the response time.
The second embodiment of the method for testing the step performance of the dc control protection system provided in the embodiment of the present application is as follows.
Referring to fig. 3, an embodiment of a step performance testing system of a dc control protection system according to the present application includes:
a first modeling unit 301, configured to preset an original dc voltage as u (t), and t e [ t ∈ [ ]1，t2]Wherein t is1＜t2Defining the filter length of the mathematical form filter of the direct current voltage, and establishing a mathematical form filter processing function M (u) of the direct current voltage according to a maximum value processing function, a minimum value processing function and a time offset function in the filter length.
The second modeling unit 302 is configured to set the offset time, the filter length, and the filter frequency value of the time offset function, and then establish a sliding average filter function n (u) of the dc voltage according to the offset time, the filter length, and the filter frequency value.
The third modeling unit 303 is configured to establish a filter function y (u) of the dc voltage according to the moving average filter function and the mathematical morphology filter processing function, where y (u) is equal to N (M (u (t)).
The overshoot point analysis unit 304 is configured to obtain an initial value vb of the curve according to a curve starting point tb, a step starting point t0, and a curve ending point te of the curve of the filter function, and a first time td1, a second time td2, and a third time td3 corresponding to the first ending value ve, the second ending value ve90, and the third ending value ve 110.
The response time calculating unit 305 is configured to determine whether an overshoot point exists in the curve according to the first end value ve, the second end value ve90, and the third end value ve110, if yes, calculate a response time according to the step start point and the fourth end value ve130, otherwise, calculate a response time according to the step start point and the second end value.
It is clear to those skilled in the art that, for convenience and brevity of description, the specific working processes of the above-described systems and units may refer to the corresponding processes in the foregoing method embodiments, and are not described herein again.
The terms "first," "second," "third," "fourth," and the like in the description of the present application and in the above-described drawings are used for distinguishing between similar elements and not necessarily for describing a particular sequential or chronological order. It is to be understood that the data so used is interchangeable under appropriate circumstances such that the embodiments of the application described herein are, for example, capable of operation in sequences other than those illustrated or otherwise described herein. Furthermore, the terms "comprises," "comprising," and "having," and any variations thereof, are intended to cover a non-exclusive inclusion, such that a process, method, system, article, or apparatus that comprises a list of steps or elements is not necessarily limited to those steps or elements expressly listed, but may include other steps or elements not expressly listed or inherent to such process, method, article, or apparatus.
It should be understood that in the present application, "at least one" means one or more, "a plurality" means two or more. "and/or" for describing an association relationship of associated objects, indicating that there may be three relationships, e.g., "a and/or B" may indicate: only A, only B and both A and B are present, wherein A and B may be singular or plural. The character "/" generally indicates that the former and latter associated objects are in an "or" relationship. "at least one of the following" or similar expressions refer to any combination of these items, including any combination of single item(s) or plural items. For example, at least one (one) of a, b, or c, may represent: a, b, c, "a and b", "a and c", "b and c", or "a and b and c", wherein a, b, c may be single or plural.
In the several embodiments provided in the present application, it should be understood that the disclosed system, apparatus and method may be implemented in other manners. For example, the above-described apparatus embodiments are merely illustrative, and for example, the division of the units is only one logical division, and other divisions may be realized in practice, for example, a plurality of units or components may be combined or integrated into another system, or some features may be omitted, or not executed. In addition, the shown or discussed mutual coupling or direct coupling or communication connection may be an indirect coupling or communication connection through some interfaces, devices or units, and may be in an electrical, mechanical or other form.
The units described as separate parts may or may not be physically separate, and parts displayed as units may or may not be physical units, may be located in one place, or may be distributed on a plurality of network units. Some or all of the units can be selected according to actual needs to achieve the purpose of the solution of the embodiment.
In addition, functional units in the embodiments of the present application may be integrated into one processing unit, or each unit may exist alone physically, or two or more units are integrated into one unit. The integrated unit can be realized in a form of hardware, and can also be realized in a form of a software functional unit.
The integrated unit, if implemented in the form of a software functional unit and sold or used as a stand-alone product, may be stored in a computer readable storage medium. Based on such understanding, the technical solution of the present application may be substantially implemented or contributed to by the prior art, or all or part of the technical solution may be embodied in a software product, which is stored in a storage medium and includes instructions for causing a computer device (which may be a personal computer, a server, or a network device) to execute all or part of the steps of the method according to the embodiments of the present application. And the aforementioned storage medium includes: a U disk, a removable hard disk, a Read-only Memory (ROM), a Random Access Memory (RAM), a magnetic disk or an optical disk, and other various media capable of storing program codes.
The above embodiments are only used for illustrating the technical solutions of the present application, and not for limiting the same; although the present application has been described in detail with reference to the foregoing embodiments, it should be understood by those of ordinary skill in the art that: the technical solutions described in the foregoing embodiments may still be modified, or some technical features may be equivalently replaced; and such modifications or substitutions do not depart from the spirit and scope of the corresponding technical solutions in the embodiments of the present application.

Claims (10)

1. A step performance test method of a direct current control protection system is characterized by comprising the following steps:
s1, presetting the original direct current voltage as u (t), and t e [ t ∈1，t2]Wherein t is1＜t2Defining the filtering length delta t of the mathematical form filter of the direct-current voltage, and establishing a mathematical form filter processing function M (u) of the direct-current voltage according to a maximum value processing function, a minimum value processing function and a time offset function in the filtering length;
s2, after setting the offset time, the filtering length and the filtering frequency value of the time offset function, establishing a sliding average filtering function N (u) of the direct current voltage according to the offset time, the filtering length and the filtering frequency value;
s3, establishing a filter function y (u) of the dc voltage according to the moving average filter function and the mathematical morphology filter processing function, wherein y (u) is N (M (u (t));
s4, acquiring an initial value vb of the curve according to a curve starting point tb, a step starting point t0 and a curve ending point te of the curve of the filter function, and a first time td1, a second time td2 and a third time td3 corresponding to a first ending value ve, a second ending value ve90 and a third ending value ve 110;
s5, judging whether the curve has overshoot points according to the first end value ve, the second end value ve90 and the third end value ve110, if yes, calculating response time according to the step starting point and the fourth end value ve130, and if not, calculating the response time according to the step starting point and the second end value.
2. The method for testing step performance of a dc control protection system according to claim 1, wherein step S1 specifically includes:
s01, presetting the original direct current voltage as u (t), and t e [ t ∈1，t2]Wherein t is1＜t2And defining the filter length of the mathematical form filter of the direct current voltage as delta t, and satisfying that delta t is less than t2-t1And the time offset function is O (u, t)o) Wherein, toIs the offset time;
s02, respectively establishing a maximum value processing function f (u) and a minimum value processing function g (u) of the mathematical morphology filter according to the maximum value and the minimum value of the direct current voltage in the filtering length;
s03, establishing the mathematical morphology filter processing function as M (u) according to the maximum processing function, the minimum processing function and the time offset function;
wherein,
f(u)＝f(u(t))＝max[t-Δt，t]u (t), and finally t e [ t ]1+Δt，t2]；
g(u)＝g(u(t))＝min[t-Δt，t]u (t), and finally t e [ t ]1+Δt，t2]；
M(u)＝M(u(t))＝0.5*O(f(g(g(f(u(t)))))+g(f(f(g(u(t)))))，to) Finally t ∈ [ t ]1+4Δt-to，t2-to]。
3. The step performance testing method of the dc control protection system according to claim 2, wherein step S4 specifically includes:
s04, defining the curve starting point of the filter function curve of the direct current voltage as tb, the step starting point as t0, the first overshoot point as t1, the curve ending point as te, the initial value as vb, the first ending value as ve, the second ending value as ve90, the third ending value as ve110 and the fourth ending value as ve 130;
s05, setting an average value of a first interval between the curve starting point and the step starting point as the initial value, where the first interval is [ tb, tb +0.02 × N ], and N is a maximum positive integer such that tb +0.02 × N < t 0;
s06, setting the average value of a second interval as the curve termination point, wherein the second interval is [ te-0.02, te ];
s07, regarding a time continuously lower than the second end value as the first time td1, the second end value is ve90 ═ ve-0.1(ve-vb), regarding a time corresponding to a first time greater than the third end value as the second time td2, the third end value ve110 ═ ve +0.1(ve-vb), and regarding a time corresponding to a first time smaller than the second end value as the third time td3, in a direction from the end point of the curve to the start point of the curve.
4. The step performance testing method of the dc control protection system according to claim 3, wherein step S5 specifically includes:
s08, traversing the derivative value of each point on the curve to the right by taking the first moment as a starting point, judging whether the derivative value is not larger than the first point corresponding to a zero value, if so, executing a step S09, and if not, setting the response time tr to td1-t 0;
s09, using the first point as the first overshoot point, when the value of the first overshoot point is greater than the fourth end value, setting the response time tr to max (td2, td3) -t0, when the value of the first overshoot point is less than the fourth end value, traversing the corresponding value on the curve to the right, setting the corresponding time of the first point which is less than the third end value as a fourth time td4, when the fourth time is equal to the second time, setting the response time tr to td3-t0, and when the fourth time is less than the second time, setting the response time tr to max (td2, td3) -t 0.
5. The step performance testing method of the dc control protection system according to claim 2, wherein the time offset function is:
O(u，to)＝O(u(t)，to)＝u(t-to) Finally t ∈ [ t ]1-to，t2-to]；
In the formula, toIs offset by a time to the left.
6. The method for testing the step performance of the DC control protection system according to claim 5, wherein the moving average filter function N (u) is:
Figure FDA0002942813600000031
the final t e [ t ∈ ]1+0.01，t2-0.01]。
7. A step performance test system of a direct current control protection system is characterized by comprising:
a first modeling unit for presetting an original DC voltage as u (t), t e [ t ∈ [ [ t ]1，t2]Wherein t is1＜t2Defining the filtering length of a mathematical form filter of the direct current voltage, and establishing a mathematical form filter processing function M (u) of the direct current voltage according to a maximum value processing function, a minimum value processing function and a time offset function in the filtering length;
the second modeling unit is used for establishing a sliding average filtering function N (u) of the direct-current voltage according to the offset time, the filtering length and the filtering frequency value after setting the offset time, the filtering length and the filtering frequency value of the time offset function;
a third modeling unit, configured to establish a filter function y (u) of the dc voltage according to the moving average filter function and the mathematical morphology filter processing function, where y (u) ═ N (M (u (t));
the overshoot point analysis unit is used for acquiring an initial value vb of the curve according to a curve starting point tb, a step starting point t0 and a curve ending point te of the curve of the filter function, and a first time td1, a second time td2 and a third time td3 which correspond to the first ending value ve, the second ending value ve90 and the third ending value ve 110;
and the response time calculating unit is used for judging whether the curve has an overshoot point according to the first termination value ve, the second termination value ve90 and the third termination value ve110, if so, calculating the response time according to the step starting point and the fourth termination value ve130, and otherwise, calculating the response time according to the step starting point and the second termination value.
8. The system for testing step performance of a dc control protection system according to claim 7, wherein the first modeling unit is specifically configured to:
presetting the original direct current voltage as u (t), t e [ t ∈ [ [ t ]1，t2]Wherein t is1＜t2And defining the filter length of the mathematical form filter of the direct current voltage as delta t < t2-t1And the time offset function is O (u, t)o) Wherein, toIs the offset time;
respectively establishing a maximum value processing function f (u) and a minimum value processing function g (u) of the mathematical form filter according to the maximum value and the minimum value of the direct-current voltage in the filtering length;
establishing a mathematical form filter processing function M (u) according to the maximum processing function, the minimum processing function and the time offset function;
wherein,
f(u)＝f(u(t))＝max[t-Δt，t]u (t), and finally t e [ t ]1+Δt，t2]；
g(u)＝g(u(t))＝min[t-Δt，t]u (t), and finally t e [ t ]1+Δt，t2]；
M(u)＝M(u(t))＝0.5*O(f(g(g(f(u(t)))))+g(f(f(g(u(t)))))，to) Finally t ∈ [ t ]1+4Δt-to，t2-to]。
9. The system for testing step performance of a dc control protection system according to claim 8, wherein the overshoot point analyzing unit is specifically configured to:
defining the curve starting point of a filter function curve of the direct-current voltage as tb, the step starting point as t0, the first overshoot point as t1, the curve ending point as te, the initial value as vb, the first ending value as ve, the second ending value as ve90, the third ending value as ve110 and the fourth ending value as ve 130;
setting an average value of a first interval between the curve starting point and the step starting point as the initial value, wherein the first interval is [ tb, tb +0.02 × N ], and N is a maximum positive integer such that tb +0.02 × N < t 0;
setting the average value of a second interval as the curve termination point, wherein the second interval is [ te-0.02, te ];
in the direction from the curve ending point to the curve starting point, a time continuously lower than the second ending value is taken as the first time td1, the second ending value is ve90 ═ ve-0.1(ve-vb), a time corresponding to the first time being greater than the third ending value is taken as the second time td2, the third ending value is ve110 ═ ve +0.1(ve-vb), and a time corresponding to the first time being less than the second ending value is taken as the third time td 3.
10. The system for testing step performance of a dc control protection system according to claim 9, wherein the response time calculating unit is specifically configured to:
traversing the derivative value of each point on the curve to the right by taking the first moment as a starting point, judging whether the derivative value is not greater than the first point corresponding to a zero value, if so, executing a step S09, otherwise, setting the response time to tr to td1-t 0;
and taking the first point as the first overshoot point, when the value of the first overshoot point is greater than the fourth end value, setting the response time tr to be max (td2, td3) -t0, when the value of the first overshoot point is less than the fourth end value, traversing the value corresponding to each point on the curve to the right, setting the time corresponding to the first point which is less than the third end value as a fourth time td4, when the fourth time is equal to the second time, setting the response time tr to be td3-t0, and when the fourth time is less than the second time, setting the response time tr to be max (td2, td3) -t 0.
CN202110183805.4A 2021-02-10 2021-02-10 Step performance testing method and system of direct current control protection system Active CN112904133B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110183805.4A CN112904133B (en) 2021-02-10 2021-02-10 Step performance testing method and system of direct current control protection system

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110183805.4A CN112904133B (en) 2021-02-10 2021-02-10 Step performance testing method and system of direct current control protection system

Publications (2)

Publication Number Publication Date
CN112904133A true CN112904133A (en) 2021-06-04
CN112904133B CN112904133B (en) 2022-02-25

Family

ID=76123505

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110183805.4A Active CN112904133B (en) 2021-02-10 2021-02-10 Step performance testing method and system of direct current control protection system

Country Status (1)

Country Link
CN (1) CN112904133B (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115079577A (en) * 2022-07-22 2022-09-20 浙江中控技术股份有限公司 Closed loop step test method and test device based on real-time control performance evaluation

Citations (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20040128088A1 (en) * 1996-03-27 2004-07-01 Laletin William H. Method of analyzing the time-varying electrical response of a stimulated target substance
CN101498761A (en) * 2008-02-02 2009-08-05 北京芯慧同用微电子技术有限责任公司 Test method for step response performance of phase-locked loop system
CN103684193A (en) * 2014-01-07 2014-03-26 南京埃斯顿自动化股份有限公司 Parameter setting method of alternating current servo system controller
CN103762595A (en) * 2014-01-24 2014-04-30 阳光电源股份有限公司 Frequency locking control method and system for suppressing power harmonics and amplitude disturbance
CN104035051A (en) * 2014-06-16 2014-09-10 中国科学院长春光学精密机械与物理研究所 Input voltage step response testing method of DC (Direct Current) power converter
KR20180032971A (en) * 2016-09-23 2018-04-02 한국전력공사 Apparatus and method for eliminating harmonics of hvdc system
CN109327152A (en) * 2017-07-28 2019-02-12 南京理工大学 Grid-connected current ring critical damping parameter determination method comprising digital control delay
CN111024186A (en) * 2019-12-31 2020-04-17 上海市计量测试技术研究院 Microwave level meter step response time measuring device and method
CN111276993A (en) * 2020-01-22 2020-06-12 华南理工大学 Controller parameter setting method, medium and equipment applied to high-voltage direct-current transmission
CN111650512A (en) * 2020-05-22 2020-09-11 国网天津市电力公司电力科学研究院 Follow-up load disturbance test method for automatic power generation control of thermal generator set
CN111679236A (en) * 2020-05-11 2020-09-18 国网江苏省电力有限公司营销服务中心 Direct current transient state step response delay test method, system and device

Patent Citations (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20040128088A1 (en) * 1996-03-27 2004-07-01 Laletin William H. Method of analyzing the time-varying electrical response of a stimulated target substance
CN101498761A (en) * 2008-02-02 2009-08-05 北京芯慧同用微电子技术有限责任公司 Test method for step response performance of phase-locked loop system
CN103684193A (en) * 2014-01-07 2014-03-26 南京埃斯顿自动化股份有限公司 Parameter setting method of alternating current servo system controller
CN103762595A (en) * 2014-01-24 2014-04-30 阳光电源股份有限公司 Frequency locking control method and system for suppressing power harmonics and amplitude disturbance
CN104035051A (en) * 2014-06-16 2014-09-10 中国科学院长春光学精密机械与物理研究所 Input voltage step response testing method of DC (Direct Current) power converter
KR20180032971A (en) * 2016-09-23 2018-04-02 한국전력공사 Apparatus and method for eliminating harmonics of hvdc system
CN109327152A (en) * 2017-07-28 2019-02-12 南京理工大学 Grid-connected current ring critical damping parameter determination method comprising digital control delay
CN111024186A (en) * 2019-12-31 2020-04-17 上海市计量测试技术研究院 Microwave level meter step response time measuring device and method
CN111276993A (en) * 2020-01-22 2020-06-12 华南理工大学 Controller parameter setting method, medium and equipment applied to high-voltage direct-current transmission
CN111679236A (en) * 2020-05-11 2020-09-18 国网江苏省电力有限公司营销服务中心 Direct current transient state step response delay test method, system and device
CN111650512A (en) * 2020-05-22 2020-09-11 国网天津市电力公司电力科学研究院 Follow-up load disturbance test method for automatic power generation control of thermal generator set

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
LING XU，ET AL.: "Parameter estimation and controller design for dynamic systems from the step responses based on the Newton iteration", 《NONLINEAR DYNAMICS》 *
邱进等: "特高压直流电流测量装置阶跃响应特性校验系统研究", 《电力系统保护与控制》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115079577A (en) * 2022-07-22 2022-09-20 浙江中控技术股份有限公司 Closed loop step test method and test device based on real-time control performance evaluation
CN115079577B (en) * 2022-07-22 2022-11-11 浙江中控技术股份有限公司 Closed loop step test method and test device based on real-time control performance evaluation

Also Published As

Publication number Publication date
CN112904133B (en) 2022-02-25

Similar Documents

Publication Publication Date Title
EP3304687B1 (en) Automatic network device electrical phase identification
CN110579682B (en) Fault recording data transient homologous comparison method and device
CN110687468B (en) Method and device for estimating state of charge of battery
CN112904133B (en) Step performance testing method and system of direct current control protection system
CN106645934B (en) Electricity consumption behavior based on dynamic grid outlier is opposed electricity-stealing diagnostic method and device
CN106405230B (en) Frequency measurement method and device
EP3172639A1 (en) System and method for estimating a move using object measurements
KR20170085369A (en) Method and apparatus for managing of battery capacity
CN110568373A (en) Lithium battery health state evaluation method, system, terminal and storage medium
US20200166577A1 (en) Method and Device for Monitoring a Stable Convergence Behavior of a Kalman Filter
CN109933264A (en) Graph data display methods and device
CN113821938B (en) Short-term prediction method and device for metering error state of mutual inductor
CN109061482B (en) Battery health degree prediction method and device
CN114441984A (en) Lithium battery health state estimation method
CN112043258A (en) Dynamic heart rate prediction method, device, equipment and readable storage medium
CN110149104B (en) Zero-phase-shift real-time filtering method for robot
WO2020147870A1 (en) Electric drive unit and method for temperature calculation in an electric drive unit
CN118364233B (en) Waveform segmentation processing device
Chen et al. Refined instrumental variable identification of continuous-time OE and BJ models from irregularly sampled data
CN115834859A (en) Four-eye stereoscopic vision camera 5G software calibration method and system
CN114519890B (en) Micro expression peak automatic detection method, system, equipment and storage medium
CN108542381A (en) A kind of data processing method and device
JP2016023818A (en) Method for detecting operation of air conditioner and system for detecting operation of air conditioner
Duhé et al. Recursive system identification for coefficient estimation of continuous-time fractional order systems
CN112415392B (en) Method for determining forgetting factor, electronic equipment, storage medium and device

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