Voltage source type disturbance generation method and system based on disturbance sequence linear interval fitting

Abstract

The invention discloses a voltage source type disturbance generation method and system based on disturbance sequence linear interval fitting, and the method comprises the steps: obtaining different types of voltage amplitude disturbance or frequency amplitude disturbance, and constructing a time-disturbance value sequence corresponding to each type of voltage amplitude disturbance or frequency amplitude disturbance; dividing each time-disturbance value sequence according to the size of a specified interval to obtain time sequences and disturbance value sequences of different intervals, and constructing the time sequences and disturbance value sequences of different intervals into corresponding parameter blocks; traversing the parameter blocks of different voltage amplitude perturbations or frequency amplitude perturbations to perform linear region fitting of the perturbation values, constructing voltage source type perturbation according to the fitting result of the perturbation values, and then selecting the phase and the perturbation type of the voltage source type perturbation to obtain final voltage source perturbation and outputting the final voltage source perturbation. According to the invention, the generation of multi-type voltage and frequency disturbance is realized under a unified framework, and the integrity of a test scene and the test efficiency are remarkably improved.

Description

Technical Field

[0001] The present invention relates to a voltage regulation and frequency modulation test control technology, and in particular to a voltage source type disturbance generation method and system based on disturbance sequence linear interval fitting. Background Art

[0002] With the large-scale integration of a high proportion of new energy and energy storage systems and a high proportion of power electronic equipment in the power system, power electronic equipment lacks rotating equipment such as mechanical rotors, the system's rotational inertia level is significantly reduced, and the voltage and frequency dynamic characteristics of the equipment show strong nonlinearity and multi-time scale coupling characteristics. The current power system dominated by traditional synchronous machines faces severe safety and stability challenges.

[0003] As voltage-source devices, grid-connected renewable energy and energy storage systems use control strategies to simulate generator voltage control and inertia characteristics, effectively mitigating the impact of the current "double-high" (high-voltage, high-voltage, high-inertia) on the power grid. However, their dynamic characteristics are significantly affected by their control structure and parameters. The industry currently requires laboratory dynamic response testing to assess their grid access qualifications.

[0004] Existing laboratory testing standards for grid-connected converters require the conduct of multi-type, continuous, and extreme operating condition disturbance excitation tests. However, these tests often use step-by-step, discrete disturbance excitations, and each disturbance excitation requires separate settings for duration and amplitude. For single or less variable primary frequency modulation, inertia response, and high / low voltage ride-through disturbance scenarios (such as frequency steps or voltage steps), the disturbance parameters are generally set manually for testing. If complex disturbance conditions with long durations and continuous multiple disturbances need to be simulated, manually setting the disturbance parameters is a huge workload. Moreover, for changing disturbance types, the time parameters and disturbance amplitudes at each time point need to be modified, making it difficult to reuse the original disturbance parameters. This results in long parameter modification times, low efficiency, and the risk of errors. Summary of the Invention

[0005] The technical problem to be solved by the present invention is as follows: In response to the above-mentioned problems of the prior art, a voltage source disturbance generation method and system based on linear interval fitting of the disturbance sequence are provided, which realizes the generation of multiple types of voltage and frequency disturbances under a unified framework, significantly improving the integrity of the test scenario and the test efficiency.

[0006] In order to solve the above technical problems, the technical solution adopted by the present invention is: A method for generating a voltage source type disturbance based on linear interval fitting of a disturbance sequence comprises the following steps: S101) acquiring different types of voltage amplitude disturbances or frequency amplitude disturbances, and constructing a time-disturbance value sequence corresponding to each voltage amplitude disturbance or frequency amplitude disturbance, wherein the time-disturbance value sequence is a two-dimensional array including time and disturbance value; S102) Dividing each time-disturbance value sequence into specified interval sizes to obtain time series and disturbance value sequences in different intervals, and constructing the time series and disturbance value sequences in different intervals into corresponding parameter blocks; S103) traversing parameter blocks of different voltage amplitude disturbances or frequency amplitude disturbances to perform linear region fitting of disturbance values, constructing a voltage source type disturbance based on the fitting results of the disturbance values, and then selecting the phase and disturbance type of the voltage source type disturbance to obtain the final voltage source disturbance and output it.

[0007] Furthermore, in step S101, when constructing the time-disturbance value sequence corresponding to each voltage amplitude disturbance or frequency amplitude disturbance, the time and disturbance value of the disturbance characteristic point in each voltage amplitude disturbance or frequency amplitude disturbance are recorded respectively to obtain the time series and disturbance value sequence of the disturbance characteristic point in each voltage amplitude disturbance or frequency amplitude disturbance.

[0008] Furthermore, in step S102, when constructing the time series and disturbance value sequences of different intervals into corresponding parameter blocks, the following steps are included: S201) Recording the interval index of each time-disturbance value sequence, and recording the boundary information of the time series and disturbance value sequence of each interval; S202) Dividing the first array and the second array into array segments corresponding to each time-disturbance value sequence; S203) The current interval index of each time-disturbance value sequence is stored in the corresponding array segment in the first array, the boundary information of the time series and the boundary information of the disturbance value sequence of the current interval of each time-disturbance value sequence is stored in the corresponding array segment in the second array, and each array segment in the second array is used as a parameter block.

[0009] Furthermore, in step S103, when traversing parameter blocks of different voltage amplitude disturbances or frequency amplitude disturbances to perform linear region fitting of disturbance values, the following steps are included: S301) traverse all parameter blocks; S302) If the simulation time of the current parameter block does not exceed the end boundary of the time series in the current parameter block and does not reach the total time of the time-disturbance value sequence corresponding to the current parameter block, then the current fitting result of the disturbance value is calculated based on the boundary information of the time series in the current parameter block, the boundary information of the disturbance value sequence in the current parameter block, and the simulation time of the current parameter block; S303) If the simulation time of the current parameter block exceeds the end boundary of the time series in the current parameter block and does not reach the total time of the time-disturbance value sequence corresponding to the current parameter block, then the current interval index of the corresponding array segment in the first array is updated to the index of the next interval, and then the boundary information of the time series and the boundary information of the disturbance value sequence of the corresponding interval are obtained according to the index of the next interval, and the boundary information of the time series and the boundary information of the disturbance value sequence in the current parameter block are updated. Finally, the current fitting result of the disturbance value is calculated according to the boundary information of the time series in the current parameter block, the boundary information of the disturbance value sequence in the current parameter block, and the simulation time of the current parameter block; S304) Jump to step S301 until the preset number of disturbance rotations is reached.

[0010] Furthermore, each array segment in the first array further includes an operation flag of the corresponding time-disturbance value sequence, and before step S301, the method further includes: detecting a disturbance start instruction, and if the disturbance start instruction is detected, setting the operation flag in each array segment in the first array to a first value; Before step S302, the method further includes: checking the running flag in the array segment corresponding to the current parameter block in the first array, and if the running flag is a first value, executing step S302; if the running flag is a second value, using the preset initial value as the current fitting result of the disturbance value; Before step S304 , the method further includes: if the simulation time of the current parameter block reaches the total time of the time-disturbance value sequence corresponding to the current parameter block, setting the running flag in the array segment corresponding to the current parameter block in the first array to the second value.

[0011] Furthermore, when the current fitting result of the disturbance value is calculated based on the boundary information of the time series in the current parameter block, the boundary information of the disturbance value sequence in the current parameter block, and the simulation time of the current parameter block, the calculation formula is as follows:

[0012] in, and They are the left boundary time and right boundary time in the boundary information of the time series, and are the left boundary value and the right boundary value in the boundary information of the perturbation value sequence, is the simulation time of the current parameter block.

[0013] Furthermore, the step of detecting the disturbance start instruction specifically includes: If a rising edge of a pulse is detected, the values ​​of the left and right edges are recorded. If no rising edge of a pulse is detected, the values ​​of the left and right edges are set to zero. Input the values ​​of the left and right edges into the judgment calculation expression. If the Boolean value of the judgment calculation expression is 1, the disturbance start instruction is received. The judgment calculation expression is as follows:

[0014] in, and They are the values ​​of the left and right edges of the pulse rising edge respectively.

[0015] Furthermore, when constructing a voltage source type disturbance based on the fitting result of the disturbance value, it includes: If the disturbance value is a disturbance value of the voltage amplitude disturbance, discrete integration is performed on the specified frequency value to obtain the three-phase phase, and the fitting result of the disturbance value and the three-phase phase are substituted into the voltage source type disturbance construction formula to obtain the voltage source type disturbance; If the disturbance value is a disturbance value of frequency amplitude disturbance, discrete integration is performed on the fitting result of the disturbance value to obtain the three-phase phase, and the specified voltage amplitude and three-phase phase are substituted into the voltage source type disturbance construction formula to obtain the voltage source type disturbance.

[0016] Furthermore, the voltage source type disturbance construction formula is as follows:

[0017]

[0018]

[0019]

[0020]

[0021]

[0022] in, 、 、 are the phases of the three phases A, B, and C respectively. 、 、 is the frequency amplitude of the three phases A, B and C, 、 、 is the voltage amplitude of phases A, B, and C, 、 、 A, B, C are three-phase AC voltage sources, is the sampling time interval.

[0023] The present invention also proposes a voltage source type disturbance generation system based on linear interval fitting of a disturbance sequence, comprising a processor and a computer storage medium, wherein the computer storage medium stores a computer program, and the computer program is executed by the processor to implement the steps of the voltage source type disturbance generation method based on linear interval fitting of a disturbance sequence.

[0024] Compared with the prior art, the advantages of the present invention are: The present invention constructs a time-disturbance value sequence based on the disturbance characteristics, divides the time-disturbance value sequence into specified interval sizes, and then performs real-time linear fitting output for each interval of different disturbance sequences, thereby realizing the generation of various types of disturbances. This reduces the workload of setting disturbance values ​​section by section when simulating complex disturbance conditions with long-term, continuous multiple disturbances, and significantly improves the integrity of the test scenario and the test efficiency. BRIEF DESCRIPTION OF THE DRAWINGS

[0025] Figure 1 The following is a brief flow chart of a method according to an embodiment of the present invention.

[0026] Figure 2 This is a detailed flowchart of performing linear region fitting of disturbance values ​​by traversing parameter blocks of different voltage amplitude disturbances or frequency amplitude disturbances in the method of an embodiment of the present invention.

[0027] Figure 3 Schematic diagram of the principle of switching the disturbance fitting interval.

[0028] Figure 4 Detailed flowchart of selecting the phase and disturbance type of the voltage source disturbance in the method according to an embodiment of the present invention.

[0029] Figure 5 Schematic diagram of various disturbances generated in the method according to an embodiment of the present invention.

[0030] Figure 6 Schematic diagram of the output of a three-phase voltage source after a disturbance signal is superimposed in the method of an embodiment of the present invention. DETAILED DESCRIPTION

[0031] The present invention will be further described below in conjunction with the accompanying drawings and specific preferred embodiments, but the scope of protection of the present invention is not limited thereby.

[0032] Example 1 This embodiment proposes a voltage source type disturbance generation method based on linear interval fitting of disturbance sequence. This method is oriented to new energy and energy storage systems. By constructing a time-disturbance value sequence based on the voltage source type voltage amplitude and frequency amplitude disturbance characteristics, the disturbance sequence interval is linearly fitted in real time, and voltage source type voltage and frequency disturbances are constructed based on the fitting results. The output voltage source type disturbance is selected according to the disturbance type and phase, and the generation of multiple types of voltage and frequency disturbances is achieved under a unified framework, which significantly improves the integrity of the test scenario and the test efficiency. This method not only provides an efficient tool for laboratory dynamic characteristic testing of converter systems of new energy and energy storage systems, but also provides a reusable generation framework for the generation of multiple types of voltage and frequency disturbances, effectively helping to build a laboratory test system for grid-type converters under new power systems.

[0033] like Figure 1 As shown, the method of this embodiment includes the following steps: S101) acquiring different types of voltage amplitude disturbances or frequency amplitude disturbances, and constructing a time-disturbance value sequence corresponding to each voltage amplitude disturbance or frequency amplitude disturbance, wherein the time-disturbance value sequence is a two-dimensional array including time and disturbance value; S102) Dividing each time-disturbance value sequence into specified interval sizes to obtain time series and disturbance value sequences in different intervals, and constructing the time series and disturbance value sequences in different intervals into corresponding parameter blocks; S103) traversing parameter blocks of different voltage amplitude disturbances or frequency amplitude disturbances to perform linear region fitting of disturbance values, constructing a voltage source type disturbance based on the fitting results of the disturbance values, and then selecting the phase and disturbance type of the voltage source type disturbance to obtain the final voltage source disturbance and output it.

[0034] In step S101 of this embodiment, based on the test standard requirements, the tested voltage and frequency disturbances are equivalently converted into a time-disturbance value sequence. When constructing the time-disturbance value sequence corresponding to each voltage amplitude disturbance or frequency amplitude disturbance, the characteristics of voltage amplitude disturbances or frequency amplitude disturbances such as high / low voltage ride-through, primary frequency modulation, inertia characteristics, and wide frequency scanning are specifically extracted, and the time and disturbance value of the disturbance feature points in each voltage amplitude disturbance or frequency amplitude disturbance are recorded respectively, and the time series and disturbance value sequence of the disturbance feature points in each voltage amplitude disturbance or frequency amplitude disturbance are obtained. The two-dimensional array of the time series and disturbance value sequence is used as the time-disturbance value sequence after the equivalent conversion of each voltage amplitude disturbance or frequency amplitude disturbance. The time-disturbance value sequence corresponding to the frequency amplitude disturbance can be expressed as: , the time-disturbance value sequence corresponding to the voltage amplitude disturbance can be expressed as: .

[0035] For example, the time-disturbance value sequence of the HVRT voltage disturbance is expressed as a time-voltage amplitude sequence, as shown in Table 1 below.

[0036] Table 1 Time-voltage amplitude sequence corresponding to HVRT voltage disturbance

[0037] The time-disturbance value sequence of a single FM frequency disturbance is expressed as a time-frequency amplitude sequence, as shown in Table 2 below.

[0038] Table 2 Time-frequency amplitude sequence corresponding to a single FM frequency disturbance

[0039] The voltage and frequency disturbances of the voltage source type are described uniformly through sequences. Other types of disturbances can be written in the "time-amplitude" form of the above example.

[0040] In this embodiment, step S102 constructs a block data structure to independently store the time-disturbance value sequence according to the time series and disturbance value sequence of multiple intervals. Variables such as the disturbance experience time, the current disturbance sequence interval, and the current disturbance sequence interval boundary are also constructed for the time-disturbance value sequence and the time series and disturbance value sequence of each interval, and these variables are also stored in the corresponding data blocks. Specifically, when constructing the time series and disturbance value sequences of different intervals into corresponding parameter blocks, the following steps are included: S201) Recording the interval index of each time-disturbance value sequence, and recording the boundary information of the time series and disturbance value sequence of each interval; S202) Dividing the first array and the second array into array segments corresponding to each time-disturbance value sequence; S203) The current interval index of each time-disturbance value sequence is stored in the corresponding array segment in the first array, the boundary information of the time series and the boundary information of the disturbance value sequence of the current interval of each time-disturbance value sequence is stored in the corresponding array segment in the second array, and each array segment in the second array is used as a parameter block.

[0041] In this embodiment, the first array is an integer array. In addition to storing the current interval index in each array segment, it also stores information such as the row and column characteristics, disturbance type, disturbance phase, and disturbance quantity of the corresponding time-disturbance value sequence, as shown in Table 3. Array indexes 0 to 9 of the integer array record common disturbance information. Within each array segment, segments 10*N to 10*(N+1)-1 record the row and column numbers of the Nth time-disturbance value sequence, as well as the current interval index and run flag, where N starts at 1. The time-disturbance value sequence is a two-dimensional array containing only two columns: time and disturbance value, and the number of rows is the number of disturbance feature points.

[0042] Table 3 Integer array storage information

[0043] In this embodiment, the second array uses a floating-point array. In addition to storing the boundary information of the time series and the boundary information of the disturbance value sequence corresponding to the current interval of the time-disturbance value sequence in each array segment, it also stores the time series data and disturbance value sequence data of the current interval, the calculation step size, the calculation time, the start instruction detection interval edge information, etc., as shown in Table 4. Array indexes 0 to 99 of the floating-point array are used to store program execution common information and external start instruction edge detection data. In each array segment, segments 100*N to 100*(N+1)-1 record the Nth time-disturbance value sequence point information, of which segments 100*N to 100*N+9 record the information of the current interval of the Nth time-disturbance value sequence, and segments 100*N+10 to 100*(N+1)-1 record the Nth time-disturbance value sequence point information, where N starts at 1.

[0044] Table 4 Floating-point array storage information

[0045] In this embodiment, after performing real-time linear fitting of the disturbance value on the stored disturbance sequence and outputting the fitting value in step S103, voltage source type voltage and frequency disturbances are constructed based on the fitting results, and the voltage source type disturbance is selected and output according to the disturbance type and phase.

[0046] In step S103, when traversing the parameter blocks of different voltage amplitude disturbances or frequency amplitude disturbances to perform linear region fitting of disturbance values, as shown in FIG. Figure 2 As shown, first, the external disturbance start command is detected. After the disturbance start command is detected, the current fitting interval of the disturbance sequence is determined. Based on the disturbance experience time, the multi-type disturbance is linearly fitted in real time and the fitting value is output. Specifically, the following steps are included: S300) detecting a disturbance start instruction. If a disturbance start instruction is detected, setting a running flag in each array segment in the first array to a first value, which in this embodiment is 1, and initializing the simulation time of all parameter blocks to 0; The external start command is generally a short step pulse with a value of 0 or 1. By detecting the rising edge of the pulse, it can be determined whether the start command has been received. The steps for detecting the disturbance start command include: If a rising edge of a pulse is detected, the values ​​of the left and right edges are recorded. If no rising edge of a pulse is detected, the values ​​of the left and right edges are set to zero. Input the values ​​of the left and right edges into the judgment calculation expression. If the Boolean value of the judgment calculation expression is 1, the disturbance start instruction is received. The judgment calculation expression is as follows: (1) in, and They are the values ​​of the left and right edges of the pulse rising edge respectively.

[0047] In each simulation step, write the new values ​​of Edge_pre and Edge_now into Table 4, then read the values ​​of Edge_pre and Edge_now in Table 4 and substitute them into formula (1) to determine whether the disturbance start instruction is detected in each step. S301) traversing all parameter blocks, thereby implementing linear fitting of multiple types of disturbance sequences one by one according to the storage order of the disturbance sequences in Table 4; S301a) Checking the running flag in the array segment corresponding to the current parameter block in the first array, and if it is the first value, executing step S302; if it is the second value, using the preset initial value as the current fitting result of the disturbance value. In this embodiment, the second value is 0; S302) If the simulation time of the current parameter block does not exceed the end boundary of the time series in the current parameter block and does not reach the total time of the time-disturbance value sequence corresponding to the current parameter block, then the current fitting result of the disturbance value is calculated based on the boundary information of the time series in the current parameter block, the boundary information of the disturbance value sequence in the current parameter block, and the simulation time of the current parameter block; In this embodiment, the calculation time Tsim in Table 4 is initialized to 0 with the receipt of the external start command as the starting time. The value of Tsim is updated at the beginning of each simulation step and serves as the simulation time at the beginning of the simulation step for the parameter block participating in the linear fitting (i.e., the parameter block whose corresponding operation flag matrixN_flag in Table 3 is 1, N≥1). The starting time of each simulation step is: Tsim=Tsim+Tstep. When the simulation step start time, i.e., the simulation time tseqnow at the beginning of a simulation step = (Tsim+Tstep), is within the current time series interval, linear fitting is performed. This embodiment stipulates the following for the linear fitting of the current disturbance type in the current disturbance interval: When the total start flag total_start_flag in Table 3 is 1, the multi-type disturbance sequence array is traversed. If the start flag matrixN_flag of the Nth group disturbance sequence array is 1, the disturbance sequence time has not reached the total sequence time tend, and has not reached the right end time Valnext of the current interval, then the linear region fitting is carried out and output. The output expression is as shown in formula (2). The current fitting result of the disturbance value is calculated based on the boundary information of the time series in the current parameter block, the boundary information of the disturbance value sequence in the current parameter block, and the simulation time of the current parameter block. The calculation formula is as follows:

[0048] in, and They are the left boundary time and right boundary time in the boundary information of the time series, and are the left boundary value and the right boundary value in the boundary information of the perturbation value sequence, is the simulation time of the current parameter block.

[0049] S303) If the simulation time of the current parameter block exceeds the end boundary of the time series in the current parameter block and does not reach the total time of the time-disturbance value sequence corresponding to the current parameter block, it is necessary to first update the time series interval boundary, and then perform linear fitting based on the updated boundary information. like Figure 3 As shown in FIG, when the simulation time tseqnow=(Tsim+Tstep) at the beginning of a simulation step exceeds the boundary of the current time series interval, the current interval index (matrixN_sectionIndex, N≥1) of the corresponding array segment in the first array is updated to the index of the next interval, and then the boundary information of the time series and the boundary information of the disturbance value sequence of the corresponding interval are obtained according to the index of the next interval, and the boundary information of the time series and the boundary information of the disturbance value sequence in the current parameter block are updated. Finally, according to formula (2), the current fitting result of the disturbance value is calculated based on the boundary information of the time series in the current parameter block, the boundary information of the disturbance value sequence in the current parameter block and the simulation time of the current parameter block.

[0050] S303a) If the simulation time of the current parameter block reaches the total time of the time-disturbance value sequence corresponding to the current parameter block, the running flag in the array segment corresponding to the current parameter block in the first array is set to the second value, that is, the running flag corresponding to the current parameter block is set to 0, and the simulation time of the current parameter block is also initialized to 0. The parameter block with the running flag set to 0 does not participate in linear fitting in subsequent iterations until the next disturbance start instruction is detected in the subsequent iteration, and the running flag of the parameter block is set to 1.

[0051] S304) Update the simulation time of the current parameter block and jump to step S301 until the preset number of perturbation rotations is reached. In this embodiment, the number of perturbation rotations is specifically the number of multi-type perturbation rotations set by Paranum in Table 3, so that the linear fitting of the current perturbation interval can be performed in each rotation.

[0052] In step S103, when constructing a voltage source type disturbance based on the fitting result of the disturbance value, since the linear fitting result of the disturbance sequence is the voltage amplitude and the frequency amplitude, and the discrete integral of the frequency amplitude over time is the phase, the voltage source disturbance is constructed based on the sine wave element. The specific steps include: If the disturbance value is a disturbance value of the voltage amplitude disturbance, discrete integration is performed on a specified frequency value (such as 50 Hz) to obtain the three-phase phase, and the fitting result of the disturbance value and the three-phase phase are substituted into the voltage source type disturbance construction formula to obtain the voltage source type disturbance; If the disturbance value is a disturbance value of frequency amplitude disturbance, discrete integration is performed on the fitting result of the disturbance value to obtain the three-phase phase, and the specified voltage amplitude (such as per-unit value 1) and the three-phase phase are substituted into the voltage source type disturbance construction formula to obtain the voltage source type disturbance.

[0053] In this embodiment, the voltage source type disturbance construction formula is as follows: (3) (4) (5) (6) (7) (8) in, 、 、 are the phases of the three phases A, B, and C respectively. 、 、 is the frequency amplitude of the three phases A, B and C, 、 、 is the voltage amplitude of phases A, B, and C, 、 、 A, B, C are three-phase AC voltage sources, is the sampling time interval.

[0054] The phase angles of each phase of the three-phase AC voltage source are calculated by equations (3) to (5), and the three-phase AC voltage source is synthesized by equations (6) to (8). For voltage amplitude disturbance, the voltage amplitude is 、 、 This is the result of linear fitting, the frequency is 50Hz, the phase 、 、 It is realized by frequency discrete summation, but the initial phase of the three phases needs to be set in advance. For frequency amplitude disturbance, the frequency amplitude 、 、 This is the result of linear fitting, the voltage amplitude 、 、 Use per-unit value 1, phase 、 、 It is achieved through discrete frequency summation, but the initial phases of the three phases need to be set in advance.

[0055] Since the voltage source type disturbance constructed by fitting is not divided into disturbance phase and disturbance type, in step S103, the phase and disturbance type of the voltage source type disturbance need to be manually selected. The specific process is as follows: Figure 4 As shown: First, the disturbance phase is selected. For three-phase symmetrical disturbances, this can be further divided into high / low voltage ride-through disturbances, primary frequency modulation disturbances, inertia characteristic disturbances, and wideband scanning disturbances. For two-phase or single-phase asymmetrical disturbances, only high / low voltage ride-through disturbances are available. After selecting the phase and disturbance type for voltage source disturbances, the desired voltage source disturbance is finally output.

[0056] Step S103 completes the linear fitting of multiple types of disturbance rotation in one simulation step through the above detailed content, and realizes the generation of single disturbance, such as Figure 4 As shown, when the loop end condition is not met, such as the number of disturbances does not meet the standard, it is necessary to generate disturbances for the next simulation step. The calculation time Tsim is superimposed on the simulation step Tstep as the starting simulation time for the next round of simulation steps. The next multi-type disturbance linear interval fitting generation can be started from step S300).

[0057] In this embodiment, the multi-type disturbances finally obtained in step S103 are as follows: Figure 5 As shown, Figure 5 Each part from top to bottom shows the output after fitting the voltage amplitude disturbance sequence in high / low voltage ride-through, the output after fitting the frequency amplitude disturbance sequence of primary frequency modulation, the output after fitting the frequency amplitude disturbance sequence of inertia characteristic, and the output after fitting the frequency amplitude disturbance sequence of wide frequency sweep. The voltage source type disturbance constructed by the fitting results is shown in the figure below. Figure 6 As shown, Figure 6 shows the output of the synthetic three-phase voltage source disturbance, Figure 6 Each part of Figure 5 The various parts in correspondence.

[0058] Example 2 This embodiment proposes a voltage source type disturbance generation system based on linear interval fitting of a disturbance sequence, including a processor and a computer storage medium, wherein the computer storage medium stores a computer program, and the computer program is executed by the processor to implement the steps of the voltage source type disturbance generation method based on linear interval fitting of a disturbance sequence described in Example 1.

[0059] In summary, the present invention constructs voltage and frequency disturbances based on real-time linear fitting of disturbance sequences, and also selects output voltage source disturbances according to disturbance type and phase. The following beneficial effects are achieved: 1. A sequence representation method for voltage and frequency disturbances of voltage sources in renewable energy and energy storage systems is proposed. The disturbance characteristics are described using time-disturbance values. The disturbance sequence points correspond to the disturbance characteristic points. This method is clear in meaning and easy to operate, providing a unified model for constructing various types of voltage and frequency disturbances.

[0060] 2. A data structure for storing perturbation sequence features is constructed, using an array to independently store individual time series and perturbation value sequences, along with other characteristic variables. This data structure can be scaled to the number of perturbation sequences, providing a template for large-scale perturbation processing.

[0061] 3. A method for real-time linear fitting and generation of stored disturbance sequences is proposed. It can identify external start-up instructions in real time and output real-time linear fitting for disturbance sequence intervals, which reduces the workload of manually setting disturbance values ​​segment by segment and provides a basis for the rapid generation of long-term, multi-segment disturbances.

[0062] 4. A modular perturbation sequence selection and generation architecture has been established, which has the ability to generate voltage and frequency perturbations of voltage source type in any combination.

[0063] The above description is merely a preferred embodiment of the present invention. The scope of protection of the present invention is not limited to the above embodiment. All technical solutions based on the concept of the present invention are within the scope of protection of the present invention. It should be noted that for those skilled in the art, various improvements and modifications that do not depart from the principles of the present invention should also be considered within the scope of protection of the present invention.

Claims

1. A voltage source disturbance generation method based on linear interval fitting of a disturbance sequence, characterized in that: The following steps are involved: S101) acquiring different types of voltage amplitude disturbances or frequency amplitude disturbances, and constructing a time-disturbance value sequence corresponding to each voltage amplitude disturbance or frequency amplitude disturbance, wherein the time-disturbance value sequence is a two-dimensional array including time and disturbance value; S102) Dividing each time-disturbance value sequence into specified interval sizes to obtain time series and disturbance value sequences in different intervals, and constructing the time series and disturbance value sequences in different intervals into corresponding parameter blocks; S103) traversing parameter blocks of different voltage amplitude disturbances or frequency amplitude disturbances to perform linear region fitting of disturbance values, constructing a voltage source type disturbance based on the fitting results of the disturbance values, and then selecting the phase and disturbance type of the voltage source type disturbance to obtain the final voltage source disturbance and output it.

2. The voltage source type disturbance generation method based on the linear interval fitting of the disturbance sequence according to claim 1 is characterized in that: In step S101, when constructing the time-disturbance value sequence corresponding to each voltage amplitude disturbance or frequency amplitude disturbance, specifically, the time and disturbance value of the disturbance characteristic point in each voltage amplitude disturbance or frequency amplitude disturbance are recorded respectively to obtain the time series and disturbance value sequence of the disturbance characteristic point in each voltage amplitude disturbance or frequency amplitude disturbance.

3. The voltage source type disturbance generation method based on the linear interval fitting of the disturbance sequence according to claim 1 is characterized in that: In step S102, when constructing the time series and disturbance value sequences of different intervals into corresponding parameter blocks, the following steps are included: S201) Recording the interval index of each time-disturbance value sequence, and recording the boundary information of the time series and disturbance value sequence of each interval; S202) Dividing the first array and the second array into array segments corresponding to each time-disturbance value sequence; S203) The current interval index of each time-disturbance value sequence is stored in the corresponding array segment in the first array, the boundary information of the time series and the boundary information of the disturbance value sequence of the current interval of each time-disturbance value sequence is stored in the corresponding array segment in the second array, and each array segment in the second array is used as a parameter block.

4. The voltage source type disturbance generation method based on linear interval fitting of disturbance sequence according to claim 3, characterized in that: In step S103, when traversing parameter blocks of different voltage amplitude disturbances or frequency amplitude disturbances to perform linear region fitting of disturbance values, the following steps are included: S301) traverse all parameter blocks; S302) If the simulation time of the current parameter block does not exceed the end boundary of the time series in the current parameter block and does not reach the total time of the time-disturbance value sequence corresponding to the current parameter block, then the current fitting result of the disturbance value is calculated based on the boundary information of the time series in the current parameter block, the boundary information of the disturbance value sequence in the current parameter block, and the simulation time of the current parameter block; S303) If the simulation time of the current parameter block exceeds the end boundary of the time series in the current parameter block and does not reach the total time of the time-disturbance value sequence corresponding to the current parameter block, then the current interval index of the corresponding array segment in the first array is updated to the index of the next interval, and then the boundary information of the time series and the boundary information of the disturbance value sequence of the corresponding interval are obtained according to the index of the next interval, and the boundary information of the time series and the boundary information of the disturbance value sequence in the current parameter block are updated. Finally, the current fitting result of the disturbance value is calculated according to the boundary information of the time series in the current parameter block, the boundary information of the disturbance value sequence in the current parameter block, and the simulation time of the current parameter block; S304) Jump to step S301 until the preset number of disturbance rotations is reached.

5. The voltage source type disturbance generation method based on disturbance sequence linear interval fitting according to claim 4 is characterized in that: Each array segment in the first array further includes an operation flag of the corresponding time-disturbance value sequence. Before step S301, the method further includes: detecting a disturbance start instruction, and if the disturbance start instruction is detected, setting the operation flag in each array segment in the first array to a first value; Before step S302, the method further includes: checking the running flag in the array segment corresponding to the current parameter block in the first array, and if the running flag is a first value, executing step S302; if the running flag is a second value, using the preset initial value as the current fitting result of the disturbance value; Before step S304 , the method further includes: if the simulation time of the current parameter block reaches the total time of the time-disturbance value sequence corresponding to the current parameter block, setting the running flag in the array segment corresponding to the current parameter block in the first array to the second value.

6. The voltage source type disturbance generation method based on disturbance sequence linear interval fitting according to claim 4, characterized in that: When the current fitting result of the disturbance value is calculated based on the boundary information of the time series in the current parameter block, the boundary information of the disturbance value sequence in the current parameter block, and the simulation time of the current parameter block, the calculation formula is as follows: in, and They are the left boundary time and right boundary time in the boundary information of the time series, and are the left boundary value and the right boundary value in the boundary information of the perturbation value sequence, is the simulation time of the current parameter block.

7. The voltage source type disturbance generation method based on disturbance sequence linear interval fitting according to claim 5, characterized in that: The steps of detecting the disturbance start instruction specifically include: If a rising edge of a pulse is detected, the values ​​of the left and right edges are recorded. If no rising edge of a pulse is detected, the values ​​of the left and right edges are set to zero. Substitute the values ​​of the left edge and the right edge into the judgment calculation expression. If the Boolean value of the judgment calculation expression is 1, the disturbance start instruction is detected. The judgment calculation expression is as follows: in, and They are the values ​​of the left and right edges of the pulse rising edge respectively.

8. The voltage source type disturbance generation method based on disturbance sequence linear interval fitting according to claim 1, characterized in that: When constructing a voltage source disturbance based on the fitting results of the disturbance value, it includes: If the disturbance value is a disturbance value of the voltage amplitude disturbance, discrete integration is performed on the specified frequency value to obtain the three-phase phase, and the fitting result of the disturbance value and the three-phase phase are substituted into the voltage source type disturbance construction formula to obtain the voltage source type disturbance; If the disturbance value is a disturbance value of frequency amplitude disturbance, discrete integration is performed on the fitting result of the disturbance value to obtain the three-phase phase, and the specified voltage amplitude and three-phase phase are substituted into the voltage source type disturbance construction formula to obtain the voltage source type disturbance.

9. The voltage source type disturbance generation method based on disturbance sequence linear interval fitting according to claim 8, characterized in that: The voltage source type disturbance construction formula is as follows: in, 、 、 are the phases of the three phases A, B, and C respectively. 、 、 is the frequency amplitude of the three phases A, B and C, 、 、 is the voltage amplitude of phases A, B, and C, 、 、 A, B, C are three-phase AC voltage sources, is the sampling time interval.

10. A voltage source disturbance generation system based on linear interval fitting of disturbance sequence, characterized in that: The method comprises a processor and a computer storage medium, wherein the computer storage medium stores a computer program, and the computer program is executed by the processor to implement the steps of the voltage source type disturbance generation method based on linear interval fitting of the disturbance sequence according to any one of claims 1 to 9.

Images