CN112904133B - 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 PDFInfo
- Publication number
- CN112904133B CN112904133B CN202110183805.4A CN202110183805A CN112904133B CN 112904133 B CN112904133 B CN 112904133B CN 202110183805 A CN202110183805 A CN 202110183805A CN 112904133 B CN112904133 B CN 112904133B
- 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.)
- Active
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]
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
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 the content of the first and second substances,
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:
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 as u (t), t e [ t ∈ [ [ t ]1，t2]Wherein t is1＜t2Defining the filter length of the mathematical form filter of the direct current voltage, and establishing the mathematical form filter 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 lengthA processing function M (u);
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 the content of the first and second substances,
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:
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:
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.
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
wherein the content of the first and second substances,
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]。
Wherein the moving average filter function n (u) is:
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 the filtered curve are eliminated, the countless small maximum values of the DC voltage waveform are eliminated, the shape is changed into countless step shapes connected by a plurality of transverse lines, the vicinity of an overshoot value is flattened, in order to find the maximum point, the filtering is carried out by adopting the 20ms sliding average filtering and then is moved leftwards by 10ms to offset the time delay, and the final filtered curve y (u) is formed.
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.
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.
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 processing the function according to the maximum value and the minimum value in the filter lengthAnd a time offset function, establishing a mathematical morphological filter processing function M (u) of the DC voltage.
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]Whereint1＜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 the content of the first and second substances,
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:
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 the content of the first and second substances,
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.
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 CN112904133A (en) | 2021-06-04 |
| CN112904133B true 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) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115079577B (en) * | 2022-07-22 | 2022-11-11 | 浙江中控技术股份有限公司 | Closed loop step test method and test device based on real-time control performance evaluation |
Citations (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| 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 |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6990422B2 (en) * | 1996-03-27 | 2006-01-24 | World Energy Labs (2), Inc. | Method of analyzing the time-varying electrical response of a stimulated target substance |
-
2021
- 2021-02-10 CN CN202110183805.4A patent/CN112904133B/en active Active
Patent Citations (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| 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)
| Title |
|---|
| Parameter estimation and controller design for dynamic systems from the step responses based on the Newton iteration;Ling Xu，et al.;《Nonlinear Dynamics》;Springer;20141117;第79卷(第3期);2155-2163 * |
| 特高压直流电流测量装置阶跃响应特性校验系统研究;邱进等;《电力系统保护与控制》;20210101;第49卷(第1期);141-147 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112904133A (en) | 2021-06-04 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| EP3304687B1 (en) | Automatic network device electrical phase identification | |
| CN109901385B (en) | Advanced observation method and device | |
| CN112904133B (en) | Step performance testing method and system of direct current control protection system | |
| US20160143421A1 (en) | Makeup application assistance device, makeup application assistance system, and makeup application assistance method | |
| CN106104610B (en) | Display situation analysis apparatus, display situation analysis method, and program recording medium | |
| CN110579682B (en) | Fault recording data transient homologous comparison method and device | |
| WO2001075660A1 (en) | Method, computer program, and system for automated real-time signal analysis for detection, quantification, and prediction of signal changes | |
| CN107368768B (en) | User identification method suitable for weighing scale and weighing scale | |
| CN107609458A (en) | Emotional feedback method and device based on expression recognition | |
| CN111366887A (en) | Mutual inductor transient sampling precision calibration method, device, terminal and medium | |
| EP2265965B1 (en) | Relative noise of a measured signal | |
| Garnier et al. | Latest developments for the Matlab CONTSID toolbox | |
| CN110728712B (en) | Furnace mouth positioning method and device based on machine vision, storage medium and electronic terminal | |
| 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 | |
| CN110149104A (en) | A kind of robot zero phase-shift Real-Time Filtering method | |
| CN115982269A (en) | Test data extraction method and system | |
| CN113820600B (en) | Alternating current resistance testing system, method and device and computer readable storage medium | |
| CN105550396B (en) | A kind of non parametric system identification method based on step response correlation analysis | |
| CN108542381A (en) | A kind of data processing method and device | |
| CN103401534B (en) | Based on the method for the final voltage fast prediction of RC circuit discharge and recharge | |
| Kotas | Projective filtering of time-aligned ECG beats for repolarization duration measurement | |
| Barbé et al. | An ARMA time series approach for analyzing long memory dynamics in measurements | |
| CN113473523B (en) | Wireless node proximity sensing method based on node neighbor relation and RSSI frequency distribution | |
| Allafi et al. | Identification of Fractional-Order Continuous-Time Hybrid Box-Jenkins Models Using Refined Instrumental Variable Continuous-Time Fractional-Order 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 |