A stable constraint-based camera siting method

Abstract

The application provides a method for site selection of a phase modifier under synchronous stability constraint, comprising the following steps: constructing an index matrix of weighted active phase sensitivity of nodes based on active phase sensitivity of a power system and active power variation of nodes; selecting a maximum value element in the index of weighted active phase sensitivity of nodes, and determining a critical value of synchronous instability according to the size of elements in the column where the maximum value element is located in the index matrix; and determining the site selection range of the phase modifier under the synchronous stability constraint according to the critical value of synchronous instability in the node corresponding to the column where the maximum value element is located. The method overcomes the deficiency that the site selection of the phase modifier only considers small disturbance stability and overvoltage constraint in the current research results, and can provide a reference for site selection and capacity determination of the phase modifier in the sending end power grid under the constraint of synchronous stability.

Description

Technical Field

[0001] This invention belongs to the field of power systems, and particularly relates to a method for selecting a camera under stability constraints. Background Technology

[0002] In recent years, to further achieve the "carbon peak and carbon neutrality" strategy, the smart grid sector has vigorously developed new energy sources, planning to build large-scale wind and solar power bases in desert, Gobi, and arid regions, and accelerating the construction of a new energy supply and consumption system. However, the gradual phase-out of traditional energy sources must be based on the safety and reliability of new energy sources. Currently, large-scale new energy bases are prone to risks such as broadband oscillations and temporary overvoltages when transmitting energy over long distances. To improve the safety and reliability of new energy bases in the sending-end grid, the common practice is to install synchronous condensers near the new energy source to enhance the transmission capacity of new energy.

[0003] Centralized synchronous condensers are typically located at the converter station in the sending-end system to suppress transient overvoltages under typical DC faults. However, renewable energy sources generally use a step-by-step voltage boosting method to collect voltage at the converter bus for centralized transmission, resulting in significant electrical distances from the rectifier station. Considering the principle of local reactive power compensation, distributed synchronous condensers of a certain capacity are needed to effectively mitigate the amplifying effect of renewable energy on transient overvoltages in the sending-end system while simultaneously enhancing the voltage support strength of the sending-end grid. Therefore, coordinating the deployment of distributed synchronous condensers near renewable energy bases with centralized synchronous condensers near converter stations is a pressing issue that needs to be addressed.

[0004] Site selection and capacity optimization are two important aspects of synchronous condenser (SCDC) optimization. Site selection mainly determines the proposed installation node of the SCDC through relevant indicators, while capacity optimization is mainly based on economic efficiency and related constraints. As a special type of synchronous machine, current research on SCDC site selection neglects the synchronization stability issues that arise after the SCDC is connected, which can easily lead to unreasonable site selection. Summary of the Invention

[0005] To address the shortcomings and deficiencies of existing technologies in synchronous condenser location research, which neglects the synchronization and stability issues after synchronous condensers are connected and is prone to unreasonable location selection, this invention proposes a synchronous condenser location method under stability constraints. Based on the active phase sensitivity and the node active power change quantity, the node active weighted phase sensitivity index can be used to screen the most severely faulty nodes, thereby providing a method for determining the location range of synchronous condensers.

[0006] This invention proposes a camera location selection method under stability constraints, comprising:

[0007] Based on the active phase sensitivity of the power system and the change in active power at nodes, an index matrix of weighted active phase sensitivity at nodes is constructed.

[0008] Select the extreme value element in the active weighted phase sensitivity index of the node, and determine the critical value of synchronization instability based on the size of the element in the column where the extreme value element is located in the index matrix.

[0009] In the nodes corresponding to the columns containing the extreme values, the location range of the synchronous condenser under the synchronous stability constraint is determined based on the critical value of synchronous instability.

[0010] Optionally, the active power weighted phase sensitivity index of the node is:

[0011] J p =J -1 ΔP L ;

[0012] Among them, J p J is a matrix composed of active power weighted phase sensitivity indices, and J is a matrix composed of active power phase sensitivity indices; ΔP L It is a matrix composed of changes in the active power of nodes.

[0013] Optionally, the matrix formed by the active phase sensitivity is:

[0014]

[0015] In the formula, diag represents the diagonalized matrix; It is the derived matrix of the modified nodal admittance matrix; U L It is a column vector of the voltage magnitudes at each node in steady state;

[0016] in, n is the number of equipment nodes in the power system, and m is the number of passive nodes in the power system. It is the modified nodal admittance matrix The element in δ ii It is the voltage phase at node i, δ ij It is the voltage phase difference between node i and node j.

[0017] Optionally, the column vector of the voltage amplitude of each node in steady state is determined by the steady-state power flow value of the power system without a synchronous condenser.

[0018] Optionally, the matrix formed by the changes in active power at the nodes is:

[0019] ΔP L =P0-diag(U pxf )I vxf -diag(U pyf )I vyf ;

[0020] Wherein, P0 is the column vector of active power injected into each node of the power system without a synchronous condenser, determined by the steady-state value of the system power flow; Upxf U pyf These are the x and y components of the voltage at each node in the global coordinate system after a fault; I vxf I vyf These represent the x and y components of the current at each node in the global coordinate system after a fault.

[0021] Optionally, the extreme value element is the maximum value element of the active weighted phase sensitivity of the nodes in the index matrix.

[0022] Optionally, determining the critical value for synchronization instability based on the size of the elements in the column containing the extreme value element in the index matrix includes:

[0023] Determine the minimum element in the column containing the maximum value element, and take the node corresponding to the row number of the minimum element in the index matrix as the center position of the camera condenser location;

[0024] Analyze the power angle characteristics of other nodes after connecting the synchronous condenser when a specified fault occurs in the power system, assuming that a synchronous condenser is installed at the central location. Combine the equal area rule to determine the critical node for synchronous stability constraint.

[0025] Determine the element corresponding to the row number of the critical node in the column containing the maximum value element, and use the active weighted phase sensitivity of the node corresponding to the element as the critical value for synchronous instability.

[0026] Optionally, determining the location range of the synchronous condenser under synchronization stability constraints based on the synchronization instability critical value in the node corresponding to the column containing the extreme value element includes:

[0027] The nodes in the column containing the extreme value element that are lower than the synchronous instability threshold are used as the location of the synchronous condenser under the synchronous stability constraint.

[0028] The beneficial effects of the technical solution provided by this invention are:

[0029] This invention overcomes the shortcomings of current research results, which only consider small disturbance stability and overvoltage constraints in the selection of synchronous condensers. It can provide a reference for the selection and capacity determination of synchronous condensers in the sending-end power grid under the constraint of synchronous stability. Attached Figure Description

[0030] To more clearly illustrate the technical solution of the present invention, the accompanying drawings used in the description of the embodiments will be briefly introduced below. Obviously, the accompanying drawings described below are only some embodiments of the present invention. For those skilled in the art, other drawings can be obtained based on these drawings without creative effort.

[0031] Figure 1This is a flowchart illustrating a stable constraint-based camera location selection method proposed in an embodiment of the present invention.

[0032] Figure 2 This is a power system topology diagram proposed in an embodiment of the present invention;

[0033] Figure 3 A time-domain simulation diagram of a fault with a phase shifter installed on node 6, which has the lowest fault sensitivity in the region, as proposed in this embodiment of the invention.

[0034] Figure 4 A time-domain simulation diagram of a fault with a synchronous condenser installed on node 34, which has the highest fault sensitivity in the region, as proposed in this embodiment of the invention.

[0035] Figure 5 A time-domain simulation diagram of a synchronous instability fault occurring when a synchronous condenser is added to the fault node 15 of the infinite tie line according to an embodiment of the present invention.

[0036] Figure 6 The time-domain simulation diagram shows the failure of synchronous instability when a synchronous condenser is installed on the fault node 18 of the infinite tie line proposed in this embodiment of the invention. Detailed Implementation

[0037] To make the objectives, technical solutions, and advantages of the embodiments of the present invention clearer, the technical solutions of the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only some embodiments of the present invention, and not all embodiments. Based on the embodiments of the present invention, all other embodiments obtained by those skilled in the art without creative effort are within the scope of protection of the present invention.

[0038] The terms "first," "second," "third," "fourth," etc. (if present) in the specification, claims, and accompanying drawings of this invention are used to distinguish similar objects and are not necessarily used to describe a specific order or sequence. It should be understood that such data can be interchanged where appropriate so that embodiments of the invention described herein can be implemented in orders other than those illustrated or described herein.

[0039] It should be understood that in the various embodiments of the present invention, the sequence number of each process does not imply the order of execution. The execution order of each process should be determined by its function and internal logic, and should not constitute any limitation on the implementation process of the embodiments of the present invention.

[0040] It should be understood that in this invention, "comprising" and "having," and any variations thereof, are intended to cover non-exclusive inclusion, for example, a process, method, system, product, or device that includes a series of steps or units is not necessarily limited to those steps or units that are explicitly listed, but may include other steps or units that are not explicitly listed or that are inherent to such process, method, product, or device.

[0041] It should be understood that in this invention, "multiple" refers to two or more. "And / or" is merely a description of the relationship between related objects, indicating that three relationships can exist. For example, "and / or B" can represent: A existing alone, A and B existing simultaneously, and B existing alone. The character " / " generally indicates that the preceding and following related objects are in an "or" relationship. "Contains A, B, and C", "Contains A, B, and C" means that all three A, B, and C are contained; "Contains A, B, or C" means that one of A, B, and C is contained; "Contains A, B, and / or C" means that any one, two, or three of A, B, and C are contained.

[0042] It should be understood that in this invention, "B corresponding to A", "B corresponding to A", "A and B correspond", or "B and A correspond" means that B is associated with A, and B can be determined based on A. Determining B based on A does not mean determining B solely based on A; B can also be determined based on A and / or other information. Matching A and B is defined as a similarity between A and B that is greater than or equal to a preset threshold.

[0043] Depending on the context, "if" as used here can be interpreted as "when," "when," "in response to determination," or "in response to detection."

[0044] The technical solution of the present invention will be described in detail below with reference to specific embodiments. These specific embodiments can be combined with each other, and the same or similar concepts or processes may not be described again in some embodiments.

[0045] Example:

[0046] like Figure 1 As shown, this embodiment proposes a camera location selection method under stability constraints, including:

[0047] S1: Based on the active phase sensitivity of the power system and the change in active power at nodes, construct an index matrix for the weighted active phase sensitivity of nodes.

[0048] S2: Select the extreme value element in the active weighted phase sensitivity index of the node, and determine the critical value of synchronization instability based on the size of the element in the column where the extreme value element is located in the index matrix.

[0049] S3: In the nodes corresponding to the columns containing the extreme values, determine the location range of the synchronous condenser under the synchronous stability constraint based on the critical value of synchronous instability.

[0050] This embodiment addresses the issue that existing studies on synchronous condenser (SCE) location neglect the synchronization and stability problems that arise after SCE is connected, which can easily lead to unreasonable location selection. It proposes a node active power weighted phase sensitivity index based on active phase sensitivity and node active power change. This index can be used to screen the nodes with the most severe faults, thus providing a method for determining the location range of SCEs. This can provide a reference for the location and capacity determination of SCEs in the sending-end power grid under the constraint of synchronization and stability.

[0051] In this embodiment, the active power weighted phase sensitivity index of the node is first calculated as follows:

[0052] J p =J -1 ΔP L ;

[0053] Among them, J p J is a matrix composed of active power weighted phase sensitivity indices, and J is a matrix composed of active power phase sensitivity indices; ΔP L It is a matrix composed of changes in the active power of nodes.

[0054] Specifically, the matrix formed by the active phase sensitivity is as follows:

[0055]

[0056] In the formula, diag represents the diagonalized matrix; It is the derived matrix of the modified nodal admittance matrix; U L It is a column vector of the voltage magnitudes at each node in steady state;

[0057] in, n is the number of equipment nodes in the power system, and m is the number of passive nodes in the power system. It is the modified nodal admittance matrix The element in δ ii It is the voltage phase at node i, δ ij It is the voltage phase difference between node i and node j.

[0058] In this embodiment, the nodal admittance matrix is ​​corrected. The significance is that it equates the new energy source to a model of a constant current source with parallel impedance, where the parallel resistance is -SBk and the parallel reactance is 1 / Sk. The nodal admittance matrix of the original system is then corrected using this parallel impedance. Here, S is the capacity of the new energy equipment; k is the dynamic reactive power ratio coefficient of the new energy source; and the expression for B is:

[0059]

[0060] Among them, I max Indicates the current limiting of new energy sources; I q0 This represents the reactive current output of the new energy source in steady state.

[0061] U L The expression is:

[0062] U L =[U1,U2,...,U n ] T ;

[0063] The voltage amplitude is also determined by the steady-state power flow value of the system without a phase shifter.

[0064] In this embodiment, the matrix formed by the changes in the active power of the nodes is as follows:

[0065] ΔP L =P0-diag(U pxf )I vxf -diag(U pyf )I vyf ;

[0066] Wherein, P0 is the column vector of active power injected into each node of the power system without a synchronous condenser, determined by the steady-state value of the system power flow; U pxf U pyf These are the x and y components of the voltage at each node in the global coordinate system after a fault; I vxf I vyf These represent the x and y components of the current at each node in the global coordinate system after a fault.

[0067] U pf The formula is:

[0068]

[0069] in, It is the modified nodal impedance matrix, defined as the modified nodal admittance matrix. The inverse matrix, with elements of ΔI s ′ um It is a column vector of current mutations at each node of the system, including the faulty node.

[0070] ΔI s ′ um The calculation method is as follows:

[0071] (1) During a fault, the current surge of the equivalent constant current source of the new energy source satisfies the following equation:

[0072]

[0073] Among them, I d0 It is the active current output of the new energy source in steady state; δ p0 Let A be the steady-state voltage phase angle in the global coordinate system (without a phase shifter); the expression for A is:

[0074]

[0075] (2) There are sudden changes in current at the new energy equipment nodes and fault nodes in the system, and the following relationship is satisfied:

[0076]

[0077] Where, ΔI j It is a column vector consisting of the current surges of the equivalent constant current source of the new energy in (1) during the fault period; the subscript f represents the fault node; ΔI f It is the sudden change in current at the fault node to be determined; It is the voltage at the fault node.

[0078] (3) After one iteration, the voltage expressions for each node are:

[0079]

[0080] Where, ΔI sum It is the column vector composed of the current mutation of new energy sources and fault nodes in (1) and (2).

[0081] (4) Considering current limiting, the dynamic reactive power ratio coefficient of new energy sources is corrected as follows:

[0082]

[0083] Based on this, the corrected current surge amounts for each new energy node are:

[0084]

[0085] Where, the subscript i represents the i-th new energy source, and the bolded variables are all column vectors composed of new energy parameters, and the meaning of the elements in the vector is the same as in (1); δ pfi for The column vector formed by the voltage phase angles of each node.

[0086] (5)ΔI s ′ um It is a column vector composed of the current mutations of new energy sources and fault nodes in (4) and (2).

[0087] I vf The formula is:

[0088]

[0089] Where * represents the Hadamard product.

[0090] Based on the index matrix of active power weighted phase sensitivity of the above nodes, this embodiment uses the extreme values ​​of the elements in the index matrix to analyze the location of the synchronous condenser with better synchronization stability. Specifically, the extreme value element is the maximum value element of the active power weighted phase sensitivity of the nodes in the index matrix.

[0091] The step of determining the critical value of synchronous instability based on the size of the element in the column where the maximum value element is located in the index matrix includes: determining the minimum element in the column where the maximum value element is located, and taking the node corresponding to the row number of the minimum element in the index matrix as the center position of the synchronous condenser location.

[0092] Analyze the power angle characteristics of other nodes after connecting the synchronous condenser when a specified fault occurs in the power system, assuming a synchronous condenser is installed at the central location. Combine this with the equal area rule to determine the critical node for synchronous stability constraints.

[0093] Determine the element corresponding to the row number of the critical node in the column containing the maximum value element, and use the active weighted phase sensitivity of the node corresponding to the element as the critical value for synchronous instability.

[0094] As described above, the method for determining the location range of the synchronous condenser based on the aforementioned index threshold values ​​in this embodiment is as follows:

[0095] (1) Select J p The largest element in the array, whose corresponding column number indicates the node where the most serious synchronization stability problem occurred;

[0096] (2)J p The row number corresponding to the smallest element in the column corresponding to the largest element represents the node with the lowest sensitivity under the most severe fault condition, and the node with the best synchronization stability when equipped with a synchronous condenser is also represented by the row number corresponding to the largest element.

[0097] (3) To characterize the steady-state and fault-period power angle characteristics of the synchronous condenser after connecting to different nodes when the most severe fault occurs, and to find the node that just does not experience synchronous instability using the equal area rule, J p The index corresponding to the row number of the node in the column with the largest element is the critical value of the node's active weighted phase sensitivity. Nodes in the column with elements below this critical value constitute the location range of the synchronous condenser.

[0098] The following is based on Figure 2 Taking the power system topology shown as an example, the specific implementation process of the synchronous condenser location method described in this embodiment is as follows:

[0099] Consider a multi-unit system comprising six new energy bases (1-6). These six bases are connected to nodes 1-39, forming the power system's network topology. Node 31 is an infinite bus node. Each new energy unit has a rated capacity of 1.5MW. After connecting to the collection bus, it is connected to the main grid via a 38 / 230kV step-up transformer. The node for the synchronous condenser is yet to be determined. The system topology is as follows: Figure 2 As shown, the network parameters are those of the classic IEEE 39 system. Among them, the AC system capacity baseline value S... B It is 230MVA, and the voltage reference value U B The voltage is 230kV. Simulations will be conducted based on two fault types: the most severe regional fault and an infinite tie-line fault, taking into account the active power flow during a fault after the synchronous condenser is connected.

[0100] (1) Fault within the area

[0101] For a 6-unit 39-node system without a synchronous condenser, the per-unit values ​​of the power output of the 6 new energy bases under the system capacity baseline are set as 0.5+j0.15, 0.5+j0.15, 2+j0.5, 1+j0.25, 1+j0.25, and 1.5+j0.45, respectively, and the initial dynamic reactive power ratio coefficient of each new energy source is set to 2.

[0102] The active power weighted phase sensitivity matrix J of each node can be obtained through calculation. p Assume J p The node corresponding to the largest element is (33, 20), and the row number corresponding to the smallest element in the column corresponding to the largest element is 6. Since the infinite node 31 was eliminated during the calculation, the node with the most severe three-phase short circuit in the system is 20. When a metallic three-phase short circuit fault occurs at node 20, the relative phase change at node 6 is the smallest, while the relative phase change at node 34 is the largest. Therefore, when analyzing faults within the region later, the fault node is node 20.

[0103] First, based on the active weighted phase sensitivity matrix J of each node... p Regarding the element size, the synchronous condenser (PCC) is installed at node 6 with a capacity of 100 MVar, i.e., two 50 MVar PCCs connected in parallel. Since node 6 is directly connected to the infinite node, node 6 has the lowest power-weighted phase sensitivity, which is consistent with physical laws, thus verifying the effectiveness of the power-weighted phase sensitivity matrix for each node. When the PCC is installed at node 6, the average phase value of each renewable energy PCC is first calculated based on the power flow results. The phase value of the PCC at renewable energy base 3 is selected as representative, and then the power angle characteristics of the PCC are calculated. Without considering the damping effect, the acceleration / deceleration area of ​​the PCC can always remain equal during and after a fault, and the PCC will not experience synchronization instability.

[0104] In the time-domain simulation, a three-phase short-circuit ground fault is set at node 20, with a fault duration of 0.1s. The power angle time-domain simulation waveform when the synchronous condenser is installed at node 6 is as follows: Figure 3 As shown, its power angle first decelerates and then accelerates, without synchronous instability, which verifies the effectiveness of the proposed theory.

[0105] Next, the synchronous condenser was installed at node 34, because this node has the highest active power output of the new energy base 3, and its active power weighted phase sensitivity is also the highest. This indicates that the active power weighted phase sensitivity of each node is determined by both the equipment and the network. Observation shows that when node 20 fails, the active power output of the new energy base 3 drops to almost zero. At this time, the power fed back to the synchronous condenser is relatively small, and the risk of synchronization instability of the synchronous condenser is also relatively small. Figure 4 The time-domain simulation waveform of the synchronous condenser's power angle verifies the analysis results. The synchronous condenser accelerates slightly during the fault and immediately returns to its original operating state after the fault is cleared.

[0106] Therefore, when analyzing the synchronization stability of the synchronous condenser, faults within the region are not the most serious type of fault.

[0107] (2) Infinite connection line fault

[0108] Considering a three-phase short-circuit ground fault at node 6, the active power of each renewable energy source during the low-voltage ride-through period will be unable to be transmitted to the infinite grid. Because the transmission channel is blocked, the active power of each renewable energy source will be fed back to the synchronous condenser. After entering the low-voltage ride-through period, the equivalent active power output of each renewable energy source will change with the power angle of the synchronous condenser.

[0109] Further, based on the magnitude of the active power weighted sensitivity index of each node under the fault of node 6, the access nodes of the synchronous condenser are screened one by one. When the synchronous condenser is connected to node 15, and the fault duration is also set to 0.1s, the synchronous condenser experiences transient power angle instability, such as... Figure 5 As shown. When the synchronous condenser is connected to node 18 and the fault duration is set to 0.1s, the synchronous condenser will just avoid transient power angle instability, as... Figure 6 As shown in the figure. This further indicates that when a synchronous condenser is connected to the sending-end grid, if the transmission channel is blocked during a renewable energy fault, and the active power output of the faulted renewable energy is high, the synchronous condenser faces the risk of synchronous instability, which requires special attention during the grid planning and operation phase. At this time, the critical active power weighted sensitivity value of the node is 0.2290. When the active power weighted sensitivity value of the node is less than 0.2290, it indicates that synchronous instability will not occur when the synchronous condenser is connected, thus determining the installation range of the synchronous condenser under synchronous stability constraints.

[0110] The serial numbers in the above embodiments are for descriptive purposes only and do not represent the order in which the components are assembled or used.

[0111] The above description is merely an embodiment of the present invention and is not intended to limit the present invention. Any modifications, equivalent substitutions, improvements, etc., made within the spirit and principles of the present invention should be included within the protection scope of the present invention.

Claims

1. A method for adjusting camera location under stability constraints, characterized in that, The method for selecting the location of the camera condenser includes: Based on the active phase sensitivity of the power system and the change in active power at nodes, an index matrix of weighted active phase sensitivity at nodes is constructed. Select the extreme value element in the active weighted phase sensitivity index of the node, and determine the critical value of synchronization instability based on the size of the element in the column where the extreme value element is located in the index matrix. In the nodes corresponding to the columns containing the extreme value elements, the location range of the synchronous condenser under the synchronous stability constraint is determined based on the synchronous instability critical value. The active weighted phase sensitivity index of the node is: ； Among them, J p J is a matrix composed of active power weighted phase sensitivity indices, and J is a matrix composed of active power phase sensitivity indices; ∆P L It is a matrix composed of changes in the active power of nodes; The matrix formed by the active phase sensitivity is: ； In the formula, diag represents the diagonalized matrix; It is the derived matrix of the modified nodal admittance matrix; U L It is a column vector of the voltage magnitudes at each node in steady state; in, Where n is the number of equipment nodes in the power system, and m is the number of passive nodes in the power system. It is the modified nodal admittance matrix The elements in It is the voltage phase at node i, δ ij It is the voltage phase difference between node i and node j.

2. The method for camera location selection under stability constraints according to claim 1, characterized in that, The column vector of the voltage amplitude of each node in steady state is determined by the steady-state power flow value of the power system without a synchronous condenser.

3. The method for adjusting camera location under stability constraints according to claim 1, characterized in that, The matrix formed by the changes in active power at the nodes is as follows: ； Wherein, P0 is the column vector of active power injected into each node of the power system without a synchronous condenser, determined by the steady-state value of the system power flow; U pxf U pyf These are the x and y components of the voltage at each node in the global coordinate system after a fault; I vxf I vyf These represent the x and y components of the current at each node in the global coordinate system after a fault.

4. The method for camera location selection under stability constraints according to claim 1, characterized in that, The extreme value element is the maximum value element of the active weighted phase sensitivity of the nodes in the index matrix.

5. The camera location selection method under stability constraints according to claim 4, characterized in that, The step of determining the critical value for synchronous instability based on the size of the elements in the column containing the extreme value element in the index matrix includes: Determine the minimum element in the column containing the maximum value element, and take the node corresponding to the row number of the minimum element in the index matrix as the center position of the camera condenser location; Analyze the power angle characteristics of other nodes after connecting the synchronous condenser when a specified fault occurs in the power system, assuming a synchronous condenser is installed at the central location. Combine this with the equal area rule to determine the critical node for synchronous stability constraints. Determine the element corresponding to the row number of the critical node in the column containing the maximum value element, and use the active weighted phase sensitivity of the node corresponding to the element as the critical value for synchronous instability.

6. The method for adjusting camera location under stability constraints according to claim 1, characterized in that, The step of determining the location range of the synchronous condenser under synchronous stability constraints based on the synchronous instability critical value in the node corresponding to the column containing the extreme value element includes: The nodes in the column containing the extreme value element that are lower than the synchronous instability threshold are used as the location of the synchronous condenser under the synchronous stability constraint.

Images