CN111781453B - Fault moment-based direct current system commutation failure risk assessment method - Google Patents
Fault moment-based direct current system commutation failure risk assessment method Download PDFInfo
- Publication number
- CN111781453B CN111781453B CN202010634607.0A CN202010634607A CN111781453B CN 111781453 B CN111781453 B CN 111781453B CN 202010634607 A CN202010634607 A CN 202010634607A CN 111781453 B CN111781453 B CN 111781453B
- Authority
- CN
- China
- Prior art keywords
- commutation
- harmonic
- voltage
- commutation failure
- time
- 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
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R19/00—Arrangements for measuring currents or voltages or for indicating presence or sign thereof
- G01R19/0084—Arrangements for measuring currents or voltages or for indicating presence or sign thereof measuring voltage only
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R29/00—Arrangements for measuring or indicating electric quantities not covered by groups G01R19/00 - G01R27/00
- G01R29/18—Indicating phase sequence; Indicating synchronism
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Mathematical Physics (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Theoretical Computer Science (AREA)
- Operations Research (AREA)
- Algebra (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Inverter Devices (AREA)
Abstract
The invention provides a direct current system commutation failure risk assessment method based on fault time aiming at an alternating current-direct current interconnection system. The method firstly analyzes the influence of each subharmonic on the commutation phase voltage-time area, and further calculates the additional commutation phase voltage-time area of the harmonic whole; from the practical engineering angle, defining the influence coefficient of n-th harmonic according to the characteristics that the higher the harmonic frequency is and the smaller the corresponding harmonic content is, and taking the influence coefficient of 2-th harmonic as the key index of the additional commutation voltage-time area; and finally, evaluating the risk of commutation failure through the mapping relation between the fault moment and the key indexes. The method can effectively evaluate the risk of the commutation failure, reflects the influence of the fault moment on the commutation failure of the direct-current system, and is beneficial to the safe, stable and efficient operation of the power grid.
Description
Technical Field
The invention relates to the field of situation perception of direct current commutation failure in the technical field of electric power, in particular to a fault moment-based direct current system commutation failure risk assessment method.
Background
The direct current transmission technology based on the thyristor is widely applied to practical engineering due to the advantages of large capacity and long-distance transmission. Typical multi-feed-in direct current systems are formed in load center areas such as east China and south China, wherein 11 loops of direct current are fed in by east China as far as 2018, and the intensive access of the direct current brings new challenges to the safe and stable operation of a power system while effectively relieving the power utilization pressure, wherein the failure of direct current commutation is an important challenge. The phase commutation failure can cause the reduction of direct current transmission power, and even cause continuous phase commutation failure to cause more serious faults such as direct current blocking and the like under the condition of serious faults.
The commutation bus voltage amplitude is generally adopted in engineering to evaluate the commutation failure risk, but the evaluation method treats the commutation voltage as a fundamental wave, and in fact, the commutation voltage contains a larger harmonic component during the fault. The fault moment can influence the harmonic initial phase, and further influence the commutation process. Therefore, the influence of the fault moment on the commutation failure of the direct current system is researched, and the method has an important function of making and taking control and protection measures in time. However, a practical quantitative evaluation method is still lacking.
Disclosure of Invention
In view of the above, the invention provides a fault-time-based direct current system commutation failure risk assessment method, which includes analyzing the influence of each subharmonic on the commutation phase voltage-time area, and further calculating the additional commutation voltage-time area of the whole harmonic; from the practical engineering angle, defining the influence coefficient of n-th harmonic according to the characteristics that the higher the harmonic frequency is and the smaller the corresponding harmonic content is, and taking the influence coefficient of 2-th harmonic as the key index of the additional commutation voltage-time area; and finally, evaluating the risk of commutation failure through the mapping relation between the fault moment and the key indexes.
In order to achieve the purpose, the invention provides the following technical scheme:
a failure-time-based commutation failure risk assessment method for a direct-current system comprises the following steps:
(1) analyzing the influence of each harmonic on the phase-voltage-time area of the phase-voltage exchange;
(2) calculating the additional commutation voltage-time area of the whole harmonic;
(3) extracting a key index of additional commutation voltage-time area for measuring commutation failure risk;
(4) and calculating the value of the key index according to the fault moment, and evaluating the commutation failure risk.
Further, the influence of each harmonic in step (1) on the voltage-time area of the commutation phase is specifically expressed by the following formula:
Wherein S is n Representing the commutation voltage-time area of the nth harmonic, E n 、Respectively representing the amplitude and phase of each harmonic, omega being 100 pi (rad/s) representing the angular velocity of the fundamental frequency of the system, and the upper and lower limits of integration being the commutation starting time t β And the commutation end time t γ ;
Change the integration interval from (t) β ,t γ ) Is transformed intoS n Represented by the formula:
wherein the content of the first and second substances,μ=ω(t γ -t β ) Indicating the commutation overlap angle.
Further, the failure of phase commutation in step (2) refers to a phenomenon that when two valves of the converter perform phase commutation, a valve which is out of conduction in the phase commutation process fails to recover blocking capability in time under the action of a reverse voltage, or the phase commutation process fails to end during the action of the reverse voltage, so that the valve which is turned off is turned on again under the action of the forward voltage.
Further, the commutation failure fault is identified by the following formula:
γ<γ min
wherein γ ═ ω (t) 0 -t β ) Indicating the valve extinction angle, t 0 Is the zero crossing time of the commutation voltage; gamma ray min Indicating the valve intrinsic limit extinction angle.
Further, the integral additional commutation voltage-time area of the harmonic in the step (2) is calculated by the following formula:
further, the key index of the additional commutation voltage-time area for measuring the commutation failure risk in the step (3) is extracted through the following process:
Define index n The influence coefficient of the nth harmonic component is calculated according to the following formula:
therefore, the temperature of the molten metal is controlled,index when the harmonic frequency n is small n The following approximation is made:
known from the above equation, index n Andapproximately in a trigonometric function distribution with the amplitude value of
Since the higher the number of harmonics, the fast attenuation of the harmonic content, Δ S is approximated by S 2 Determine, i.e. from index 2 And (6) determining.
Further, when index 2 At the minimum, the commutation failure risk was considered to be the highest.
Further, the step (4) is specifically realized by the following steps:
wherein the function represents the initial phase of the fault time t _ fault and the 2 nd harmonic componentFunctional relationship between;
(42) calculating a key index 2 :
(43) Assessing commutation failure risk:
index 2 if the harmonic component is more than 0, the fault harmonic component is beneficial to the commutation process, and the commutation failure risk is lower when the numerical value is larger;
index 2 < 0 indicates that the harmonic component of the fault is in phase commutationThe process is unfavorable, and the smaller the number, the higher the risk of commutation failure.
Further, in the step (42), when the dc system is in steady operation, the extinction angle γ ═ pi- ω t γ Pi/12, trigger advance angle beta pi- ω t β ≈0.211π,index 2 Written approximately as:
compared with the prior art, the invention has the following advantages and beneficial effects:
The invention provides a fault moment-based direct-current system commutation failure risk assessment method for an alternating-current and direct-current interconnection system, which can effectively assess the commutation failure risk, reflect the influence of the fault moment on the commutation failure of the direct-current system, and is beneficial to the safe, stable and efficient operation of a power grid.
Drawings
Fig. 1 is a schematic flow chart of a method for evaluating a commutation failure risk of a dc system based on a fault time according to the present invention.
Detailed Description
The technical solutions provided by the present invention will be described in detail below with reference to specific examples, and it should be understood that the following specific embodiments are only illustrative of the present invention and are not intended to limit the scope of the present invention.
As shown in fig. 1, a method for evaluating a risk of commutation failure of a dc system based on a fault time provided in an embodiment of the present invention includes:
step S1: analyzing the influence of each harmonic on the phase-voltage-time area of the phase-voltage switching, and concretely showing that:
wherein S is n Representing the commutation voltage-time area of the nth harmonic, E n 、Respectively representing the amplitude and phase of each harmonic, omega being 100 pi (rad/s) representing the angular velocity of the fundamental frequency of the system, and the upper and lower limits of integration being the commutation starting time t β And the commutation end time t γ ;
Change the integration interval from (t) β ,t γ ) Is transformed intoS n Can be represented by the following formula:
wherein the content of the first and second substances,μ=ω(t γ -t β ) Indicating the commutation overlap angle.
Step S2: calculating the integral additional commutation voltage-time area Delta S of the harmonic wave by the following formula:
step S3: extracting a key index of additional commutation voltage-time area for measuring commutation failure risk;
the failure fault of commutation refers to the phenomenon that when two valves of the converter carry out commutation, the valve which is out of conduction in the commutation process cannot restore the blocking capability in time under the action of reverse voltage, or the commutation process cannot be finished in the reverse voltage action period, so that the valve which is turned off is conducted again under the action of forward voltage. Specifically, the following equation can be used for discrimination:
γ<γ min
wherein γ ═ ω (t) 0 -t β ) Indicating the valve extinction angle, t 0 Is the zero crossing time of the commutation voltage; gamma ray min Indicating the inherent extreme extinction angle of the valve, is generally related to the inherent characteristics of the valve.
Extracting a key index of additional commutation voltage-time area for measuring commutation failure risk, wherein the specific implementation process comprises the following steps:
define index n The influence coefficient of the nth harmonic component can be calculated according to the following formula:
therefore, the temperature of the molten metal is controlled,index when the harmonic frequency n is small n The following approximation can be made:
as shown in the above formula, index n Andapproximately in a trigonometric function distribution with the amplitude value of
Since the higher the number of harmonics, the faster the harmonic content will decay, Δ S can be approximated as being represented by S 2 Determine, i.e. by index 2 And (6) determining. When index 2 At the minimum, the commutation failure risk was considered to be the highest.
Step S4: calculating the value of a key index according to the fault moment, and evaluating the commutation failure risk, wherein the specific implementation process comprises the following steps:
wherein the function represents the initial phase of the fault time t _ fault and the 2 nd harmonic componentFunctional relationship between; the mapping can be obtained by fitting off-line simulation data, and the result is an expression form of a linear function;
(42) calculating a key index 2 :
When a general direct current system operates in a steady state, the extinction angle gamma is pi- ω t γ Pi/12, trigger advance angle beta pi- ω t β 0.211 π, hence index 2 Can be written approximately as:
(43) assessing commutation failure risk:
index 2 if the harmonic component is more than 0, the fault harmonic component is beneficial to the commutation process, and the commutation failure risk is lower when the numerical value is larger;
index 2 and < 0 indicates that the fault harmonic component is unfavorable to the commutation process, and the smaller the value, the higher the risk of commutation failure.
The technical means disclosed in the invention scheme are not limited to the technical means disclosed in the above embodiments, but also include the technical scheme formed by any combination of the above technical features. It should be noted that those skilled in the art can make various improvements and modifications without departing from the principle of the present invention, and such improvements and modifications are also considered to be within the scope of the present invention.
Claims (7)
1. A failure-time-based commutation failure risk assessment method for a direct-current system is characterized by comprising the following steps:
(1) analyzing the influence of each harmonic on the phase voltage-time area of the phase-change voltage;
(2) calculating the additional commutation voltage-time area of the whole harmonic;
(3) extracting a key index of additional commutation voltage-time area for measuring commutation failure risk; the key index of the additional commutation voltage-time area for measuring the commutation failure risk is extracted by the following process:
define index n The influence coefficient of the n-th harmonic component is calculated according to the following formula:
therefore, the temperature of the molten metal is controlled,E n 、respectively representing the amplitude and phase of each harmonic, ω -100 π rad/s representing the angular velocity of the fundamental frequency of the system,μ=ω(t γ -t β ) Representing the commutation overlap angle, t β For the start of commutation, t γ Is the commutation end time;
since the higher the number of harmonics, the faster the harmonic content will decay, Δ S is approximated by S 2 Determine, i.e. by index 2 Determining;
(4) calculating the value of a key index according to the fault moment, and evaluating the commutation failure risk; the specific implementation process is as follows:
wherein the function represents the initial phase of the fault time t _ fault and the 2 nd harmonic component A functional relationship therebetween;
(42) calculating a key index 2 :
(43) Assessing commutation failure risk:
index 2 if the harmonic component is more than 0, the fault harmonic component is beneficial to the commutation process, and the commutation failure risk is lower when the numerical value is larger;
index 2 and < 0 indicates that the fault harmonic component is unfavorable to the commutation process, and the smaller the value, the higher the risk of commutation failure.
2. The fault moment-based direct current system commutation failure risk assessment method according to claim 1, wherein the influence of each harmonic in step (1) on the commutation voltage-time area is specifically expressed by the following formula:
wherein S is n Representing the commutation voltage-time area of the nth harmonic, E n 、Respectively representing the amplitude and the phase of each harmonic, wherein omega is 100 pi rad/s to represent the angular speed of the fundamental frequency of the system, and the upper limit and the lower limit of the integral are respectively the commutation starting time t β And the commutation end time t γ ;
Change the integration interval from (t) β ,t γ ) Is transformed intoS n Represented by the formula:
3. The method for evaluating the risk of commutation failure of the dc system according to claim 1, wherein the commutation failure is a phenomenon that when two valves of the converter perform commutation, a valve that is out of conduction in a commutation process fails to recover blocking capability in time under the action of a reverse voltage, or a commutation process fails to end during the action of a reverse voltage, so that a valve that is turned off is turned back on under the action of a forward voltage.
4. The fault-time-based commutation failure risk assessment method for a direct current system according to claim 1 or 3, wherein the commutation failure is determined by the following formula:
γ<γ min
wherein γ ═ ω (t) 0 -t β ) Indicating the valve extinction angle, t 0 Is the zero crossing time of the commutation voltage; gamma ray min Indicating the valve intrinsic limit extinction angle.
6. the method for evaluating commutation failure risk of a DC system according to claim 1, wherein when index is used as the index 2 At the minimum, the commutation failure risk was considered highest.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010634607.0A CN111781453B (en) | 2020-07-04 | 2020-07-04 | Fault moment-based direct current system commutation failure risk assessment method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010634607.0A CN111781453B (en) | 2020-07-04 | 2020-07-04 | Fault moment-based direct current system commutation failure risk assessment method |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111781453A CN111781453A (en) | 2020-10-16 |
CN111781453B true CN111781453B (en) | 2022-07-29 |
Family
ID=72757940
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010634607.0A Active CN111781453B (en) | 2020-07-04 | 2020-07-04 | Fault moment-based direct current system commutation failure risk assessment method |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111781453B (en) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112421666B (en) * | 2020-11-10 | 2022-09-09 | 国网新疆电力有限公司经济技术研究院 | High-voltage direct-current power transmission system subsequent commutation failure risk assessment and inhibition method |
Citations (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101701996A (en) * | 2009-11-26 | 2010-05-05 | 华南理工大学 | Method and device for determining commutation failure of high voltage direct current system |
CN103543357A (en) * | 2013-10-23 | 2014-01-29 | 华北电力大学 | Method for prejudging phase change failure of converter valve |
CN104614640A (en) * | 2014-12-26 | 2015-05-13 | 四川大学 | Method for detecting phase change failure of high-voltage direct-current power transmission system |
KR20150142763A (en) * | 2014-06-11 | 2015-12-23 | 한국전력공사 | Apparatus for predicting and controlling commutation failure of high voltage direct current system |
CN106033889A (en) * | 2015-03-20 | 2016-10-19 | 国家电网公司 | Method for determining risk degrees of mutual influence of multi-infeed alternating current-direct current system inversion stations |
CN106803153A (en) * | 2017-01-23 | 2017-06-06 | 中国电力科学研究院 | A kind of appraisal procedure and system of many feed-in ac and dc systemses commutation failure risks |
CN107134801A (en) * | 2017-05-12 | 2017-09-05 | 国网河南省电力公司电力科学研究院 | The direct current transportation commutation failure probability acquiring method of meter and commutation failure PREDICTIVE CONTROL |
CN108418238A (en) * | 2018-03-08 | 2018-08-17 | 国网上海市电力公司 | The method occurred based on harmonic synthesis commutation coefficient assessment commutation failure |
CN110011346A (en) * | 2019-03-19 | 2019-07-12 | 国家电网有限公司 | A kind of interactional commutation failure methods of risk assessment of consideration Inverter Station |
CN110323773A (en) * | 2019-07-03 | 2019-10-11 | 中国电力科学研究院有限公司 | A kind of commutation failure prediction technique and system considering multifactor impact |
CN110460082A (en) * | 2018-05-07 | 2019-11-15 | 中国电力科学研究院有限公司 | A kind of more feed-in commutation failure of high voltage direct current system risk method of discrimination and system |
CN110518622A (en) * | 2019-09-24 | 2019-11-29 | 华北电力大学 | A kind of commutation failure detects in advance and prevention method, device |
-
2020
- 2020-07-04 CN CN202010634607.0A patent/CN111781453B/en active Active
Patent Citations (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101701996A (en) * | 2009-11-26 | 2010-05-05 | 华南理工大学 | Method and device for determining commutation failure of high voltage direct current system |
CN103543357A (en) * | 2013-10-23 | 2014-01-29 | 华北电力大学 | Method for prejudging phase change failure of converter valve |
KR20150142763A (en) * | 2014-06-11 | 2015-12-23 | 한국전력공사 | Apparatus for predicting and controlling commutation failure of high voltage direct current system |
CN104614640A (en) * | 2014-12-26 | 2015-05-13 | 四川大学 | Method for detecting phase change failure of high-voltage direct-current power transmission system |
CN106033889A (en) * | 2015-03-20 | 2016-10-19 | 国家电网公司 | Method for determining risk degrees of mutual influence of multi-infeed alternating current-direct current system inversion stations |
CN106803153A (en) * | 2017-01-23 | 2017-06-06 | 中国电力科学研究院 | A kind of appraisal procedure and system of many feed-in ac and dc systemses commutation failure risks |
CN107134801A (en) * | 2017-05-12 | 2017-09-05 | 国网河南省电力公司电力科学研究院 | The direct current transportation commutation failure probability acquiring method of meter and commutation failure PREDICTIVE CONTROL |
CN108418238A (en) * | 2018-03-08 | 2018-08-17 | 国网上海市电力公司 | The method occurred based on harmonic synthesis commutation coefficient assessment commutation failure |
CN110460082A (en) * | 2018-05-07 | 2019-11-15 | 中国电力科学研究院有限公司 | A kind of more feed-in commutation failure of high voltage direct current system risk method of discrimination and system |
CN110011346A (en) * | 2019-03-19 | 2019-07-12 | 国家电网有限公司 | A kind of interactional commutation failure methods of risk assessment of consideration Inverter Station |
CN110323773A (en) * | 2019-07-03 | 2019-10-11 | 中国电力科学研究院有限公司 | A kind of commutation failure prediction technique and system considering multifactor impact |
CN110518622A (en) * | 2019-09-24 | 2019-11-29 | 华北电力大学 | A kind of commutation failure detects in advance and prevention method, device |
Non-Patent Citations (2)
Title |
---|
Evaluation of Commutation failure Risk for HVDC Caused by Harmonic Voltage;Jiahao Chen等;《2019 IEEE 3rd Conference on Energy Internet and Energy System Integration (EI2)》;20200409;全文 * |
基于改进换相面积的直流输电换相失败判别方法;刘济豪等;《华北电力大学学报(自然科学版)》;20140228;全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN111781453A (en) | 2020-10-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US10951037B2 (en) | Method and apparatus for time-domain droop control with integrated phasor current control | |
CN103078312B (en) | Phase-changing failure suppression method based on direct current prediction control | |
Xiao et al. | Efficient approach to quantify commutation failure immunity levels in multi‐infeed HVDC systems | |
CN103762581B (en) | A kind of defence method of one-tower double-circuit DC transmission system periodicity commutation failure | |
CN103078311B (en) | Direct current predicting and setting method for inhibiting commutation failure | |
CN107147120B (en) | RBF dual neural network self-adaptive sliding mode control method of active power filter | |
CN110429635B (en) | Commutation failure prediction method based on direct current finite time domain prediction | |
CN103543357B (en) | A kind of pre-judging method of converter valve commutation failure | |
Liu et al. | Extinction angle predictive control strategy for commutation failure mitigation in HVDC systems considering voltage distortion | |
CN111781453B (en) | Fault moment-based direct current system commutation failure risk assessment method | |
Zhu et al. | Prevention and mitigation of high‐voltage direct current commutation failures: a review and future directions | |
CN104218788A (en) | Network side power factor sliding-mode self-adaptive control method of matrix rectifier | |
CN107425533B (en) | Method for guiding synchronization inverter in grid-connected operation state when voltage is unbalanced | |
Selim et al. | Fast quasi‐static time‐series analysis and reactive power control of unbalanced distribution systems | |
Ouyang et al. | Current-limit control method to prevent subsequent commutation failure of LCC-HVDC based on adaptive trigger voltage | |
CN109066656A (en) | A kind of stability of power system judgment method based on single feed-in broad sense operation short-circuit ratio | |
CN116436327A (en) | Adaptive prediction control method and system for converter based on neural network | |
CN109802380B (en) | Low-voltage current limiting control method, system and device for high-voltage direct-current transmission | |
CN111697613A (en) | Method and system for calculating critical voltage of commutation failure of direct-current power transmission system | |
CN104734172A (en) | Self-adaption PI control method for improving HVDC system suppressing phase commutation failure capacity | |
Liu et al. | An enhanced commutation failure prevention control in LCC based HVDC systems | |
CN107332267A (en) | It is a kind of to avoid extra-high voltage direct-current power from continuing relatively low control method and system | |
CN116581763A (en) | Method for selectively switching working modes of photovoltaic inverter system in power distribution network | |
CN106099918A (en) | A kind of method of simulation calculation multi-infeed DC mains frequency emergency control policy | |
CN103474964B (en) | A kind of prevent inversion from overturning intelligent overcurrent protective device and control 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 |