Electrical automation adjusting method and system

Abstract

The invention relates to the technical field of adaptive control, in particular to an electrical automation adjustment method and system, and the method comprises the following steps: obtaining a three-phase current signal, collecting single-cycle high-frequency sampling, extracting an instantaneous extreme point, recording a time coordinate and an amplitude value, calculating a time interval and a change rate, screening the distribution of the extreme point, and obtaining a harmonic distortion characteristic parameter. According to the method, on the basis of high-frequency sampling, instantaneous extreme point extraction, time interval and change rate calculation, harmonic distortion feature description enhancement, extreme value growth rate and phase drift rate are utilized, time periods with prominent fluctuation amplitude are screened, harmonic change trend prediction precision is improved, and the harmonic change trend prediction precision is improved through zero-sequence component and phase offset calculation in combination with periodic change features. The harmonic source is accurately separated, the harmonic cause identification capability is improved, the compensation current is calculated based on the harmonic source, the injection opportunity is optimized, the compensation parameter is adjusted before the harmonic peak value, and the compensation response speed and the matching precision are improved.

Description

Technical Field

[0001] The present invention relates to the field of adaptive control technology, and in particular to an electrical automation regulation method and system. Background Art

[0002] The field of adaptive control technology includes technical methods that adjust control parameters or strategies to keep the control system's expected performance under different operating conditions in response to environmental changes and changes in the system's dynamic characteristics during system operation. This technical field mainly involves control strategies based on mathematical models, self-learning control algorithms, and real-time parameter adjustment methods. Adaptive control technology is widely used in industrial automation, robotic control, power systems, aerospace, and other fields. It usually includes technical means such as parameter adaptive control, model reference adaptive control, and intelligent adaptive control to ensure that the control system can adapt to external disturbances and changes in system internal parameters, and improve control accuracy and stability.

[0003] Among them, the electrical automation regulation method refers to a technical method in which the control parameters are dynamically adjusted according to the set regulation rules by real-time monitoring of the system operation status in the electrical control system to ensure that the system operates as expected. The subject of this patent involves a real-time regulation strategy based on feedback signals, including the use of adaptive control algorithms to analyze the input signals of electrical equipment and the dynamic adjustment of voltage, current or power factor by calculating specific regulation parameters. In addition, the method also includes selecting a suitable control mode according to the system operation status, and adjusting the control parameters of the drive circuit, power electronic device or programmable controller to enable the system to maintain stable operation under different load or environmental changes.

[0004] Existing technologies are based on fixed time window mean analysis or Fourier transform, which makes it difficult to accurately identify instantaneous extreme points, and the accuracy of harmonic distortion feature extraction is limited, and the ability to identify sudden distortion is insufficient. Harmonic source identification relies on static parameter matching, which is difficult to adapt to dynamic working conditions and prone to misjudgment when load complexity increases. Compensation control uses a fixed compensation coefficient, which lacks adaptability to real-time working conditions, lacks flexibility in current amplitude and phase adjustment, and has limited harmonic suppression effects, affecting the stable operation of the system. Summary of the invention

[0005] The purpose of the present invention is to solve the shortcomings of the prior art and to propose an electrical automation regulation method and system.

[0006] In order to achieve the above object, the present invention adopts the following technical solution: an electrical automation adjustment method, comprising the following steps: S1: Obtain three-phase current signals, collect single-cycle high-frequency samples, extract instantaneous extreme points, record time coordinates and amplitudes, calculate time intervals and change rates, screen extreme point distributions, and obtain harmonic distortion characteristic parameters; S2: Based on the harmonic distortion characteristic parameters, call the instantaneous extreme value, calculate the phase drift rate and the extreme value growth rate, analyze the distortion stability, select the period with prominent fluctuation amplitude, determine the dominant harmonic component, calculate the peak time and distortion amplitude, construct the trend change curve, and calculate the harmonic distortion trend prediction value; S3: calling the harmonic distortion trend prediction value, extracting the A, B, and C phase currents in the three-phase current signal, calculating the zero-sequence component, obtaining the phase offset and harmonic distortion amplitude, comparing the periodic changes, determining the harmonic source, and obtaining the harmonic source separation result; S4: Based on the harmonic source separation result, calculate the compensation currents of phases A, B, and C, compare the multi-phase impact, calculate the injection time point of the harmonic compensation current, set the compensation current parameters before the harmonic peak arrives, and obtain the asymmetric compensation parameters.

[0007] As a further scheme of the present invention, the harmonic distortion characteristic parameters include time interval, rate of change, and extreme point distribution; the harmonic distortion trend prediction values ​​include phase drift rate, extreme value growth rate, distortion stability, dominant harmonic component, peak time, distortion amplitude, and trend change curve; the harmonic source separation results include zero-sequence component, phase offset, harmonic distortion amplitude, periodic change, and harmonic source; the asymmetric compensation parameters include A-phase compensation current, B-phase compensation current, C-phase compensation current, and injection time point of harmonic compensation current.

[0008] As a further solution of the present invention, the step of acquiring the harmonic distortion characteristic parameter is specifically as follows: S101: Acquire three-phase current signals, perform single-cycle high-frequency sampling, record current amplitudes corresponding to multiple time points, extract instantaneous extreme value points at multiple times, store the time coordinates and amplitudes of the extreme value points, and generate an instantaneous extreme value data set; S102: Based on the instantaneous extreme value data set, calculate the time intervals and amplitude change rates of adjacent extreme value points, establish a time interval sequence and a change rate sequence, filter outliers in the time interval sequence, and remove time interval data that exceeds a set threshold, using the formula: ; Calculate and generate the distribution deviation of the rate of change of extreme value points; in, Represents the distribution deviation of the rate of change of extreme points, Representative The rate of change of extreme points, represents the average value of the rate of change, represents the standard deviation of the rate of change, Representative The time interval between extreme points, represents the mean of the time interval, Represents the total number of extreme points; S103: Call the extreme point change rate distribution deviation, filter out abnormal data points that exceed the set threshold, count the proportion of abnormal data points, calculate the abnormal distribution rate of data points, combine the time interval screening results, extract extreme points that meet the set distribution characteristics, obtain extreme point distribution interval parameters, and generate harmonic distortion characteristic parameters.

[0009] As a further solution of the present invention, the step of obtaining the harmonic distortion trend prediction value is specifically as follows: S201: based on the harmonic distortion characteristic parameters, calling the instantaneous extreme value, comparing the extreme value change amplitudes in the differentiated time periods, calculating the extreme value growth rate, and screening according to the set extreme value growth rate reference value, extracting the time period with outstanding growth rate, and obtaining the extreme value growth rate screening result; S202: calling the extreme value growth rate screening result, calculating the corresponding phase drift rate for the screened time period, determining the fluctuation of the phase drift rate, screening the time interval for the fluctuation amplitude, and obtaining the phase drift rate fluctuation interval; S203: calling the phase drift rate fluctuation interval, calculating the peak time and distortion amplitude of the dominant harmonic component in the fluctuation interval, and constructing a trend change curve, using the formula: ; Calculate and obtain the predicted value of harmonic distortion trend; in, represents the predicted value of harmonic distortion trend, Representative The distortion amplitude of the dominant harmonic component, Representative The attenuation factor of the dominant harmonic component, represents the time variable, Representative The drift factor of the dominant harmonic component, Representative The phase shift of the dominant harmonic component, Represents the total number of dominant harmonic components.

[0010] As a further solution of the present invention, the step of obtaining the harmonic source separation result is specifically: S301: calling the harmonic distortion trend prediction value, extracting the A, B, and C phase currents in the three-phase current signal, calculating the root mean square value of the multi-phase current signal in multiple cycles, obtaining the intra-cycle fluctuation characteristics of the A, B, and C phase currents, and performing differential operation based on the root mean square value of the multi-phase current and the zero-sequence component, calculating the current offset of the zero-sequence component, and obtaining the zero-sequence component offset rate; S302: Based on the zero-sequence component offset rate, the phase offsets of the A, B, and C phase currents relative to the reference phase are calculated, the A and B phase currents are selected for difference calculation, and the ratio relationship of the zero-sequence component is constructed using the formula: ; Calculate the phase change trend of A-phase current and B-phase current, obtain the phase difference parameter of A-phase current and B-phase current, determine the comparison relationship between phase change and zero-sequence component, screen the key phase distortion characteristics, and obtain the main phase distortion amount; in, represents the phase offset, Representative The current signal of phase A and phase B in one cycle, The current signal representing the zero-sequence component, Represents the current difference between phases A and B, represents the maximum current during the measurement, represents the number of cycles selected during the calculation process, Represents the inverse tangent function operation, which is used to calculate the phase angle corresponding to the ratio of the current signal; S303: calling the main phase distortion variable, calculating the corresponding harmonic distortion amplitude, and selecting the extreme value points of the harmonic amplitude to calculate the fluctuation range based on the variation trend of the harmonic distortion amplitude in multiple periods, and comparing the periodic variation characteristics to determine the type and contribution of the harmonic source, and obtain the harmonic source separation result.

[0011] As a further solution of the present invention, the step of obtaining the asymmetric compensation parameter is specifically as follows: S401: Based on the harmonic source separation result, calculate the compensation currents of phases A, B, and C, extract multi-phase harmonic current components, call the harmonic current component data, calculate the peak amplitude of each phase harmonic current, and calculate the phase change trend based on the peak amplitude data. For the three-phase harmonic current data, compare the phase deviation, calculate the phase offset value between the three phases, and obtain the three-phase harmonic phase offset; S402: calling the three-phase harmonic phase offset, calculating the time offset value of the multi-phase compensation current, calling the three-phase compensation current time offset data, comparing the harmonic compensation impact, calculating the timing difference of the compensation current, obtaining the optimal injection time point of the compensation current, and calling the harmonic amplitude data at the same time, calculating the compensation amplitude, and establishing the optimal compensation current parameters; S403: Based on the optimal compensation current parameters, compensation current adjustment parameters are set before the harmonic peak value is reached, and parameters of the multi-phase compensation current are calculated using the formula: ; Obtain asymmetric compensation parameters; in, Represents the overall parameter of the compensation current, Representative Order harmonic voltage, Representative The phase angle of the harmonic voltage of order, represents the phase adjustment parameter, represents the total number of harmonic orders, Is an imaginary unit.

[0012] As a further embodiment of the present invention, the method further comprises: S5: calling the asymmetric compensation parameters, collecting data before compensation, calculating multi-band priorities, comparing amplitude changes, adjusting current amplitude and phase, and recording adaptive compensation execution data; The adaptive compensation execution data includes pre-compensation data, multi-band priorities, amplitude changes, current amplitude adjustment, and current phase adjustment.

[0013] As a further solution of the present invention, the step of acquiring the adaptive compensation execution data is specifically as follows: S501: calling the asymmetric compensation parameter, collecting data before compensation, calculating the signal amplitudes and spectrum characteristics of multiple frequency bands, obtaining the signal change rates of differentiated frequency bands, and filtering the frequency bands based on the change rates, calculating the weights of the frequency bands, and sorting them to obtain the priority sorting of multiple frequency bands; S502: Based on the multi-band priority sorting, the signal amplitude change of the selected frequency band is calculated, the initial current amplitude and phase parameters are called, the amplitude change ratio of the multi-band signal is calculated, and compared with the set reference change threshold, using the formula: ; Calculate the current adjustment factor, adjust the current amplitude and phase of multiple frequency bands, and obtain the current correction parameters; in, represents the current adjustment factor, Representative The current signal amplitude of the frequency band, represents the reference signal amplitude, Representative The weight of the frequency band, Represents the number of frequency bands to be filtered, , represents the operation of capturing the maximum and minimum values ​​in all filtered frequency bands; S503: calling the current correction parameter, adjusting the actual applied current amplitude and phase, and recording the adjusted signal response in real time, calculating the difference of the signal offset before and after the adjustment, and obtaining adaptive compensation execution data.

[0014] An electrical automation regulation system, the electrical automation regulation system is used to execute the above electrical automation regulation method, the system comprises: The harmonic distortion feature extraction module obtains three-phase current signals, extracts single-cycle high-frequency sampling data, screens instantaneous extreme value points among the sampling points, records the time coordinates and amplitudes of the instantaneous extreme value points, calculates the time intervals between consecutive instantaneous extreme value points, measures the extreme value change rate, screens out extreme value points with prominent fluctuation amplitudes of the change rate, analyzes the distribution density of extreme value points, extracts distortion frequency characteristics, calculates the relative amplitudes of the fundamental wave and each harmonic, measures the total harmonic distortion rate, establishes a distortion trend curve, and obtains harmonic distortion feature parameters; The phase shift calculation module calls the instantaneous extreme point data based on the harmonic distortion characteristic parameters, calculates the phase drift rate of each harmonic, screens the harmonic components with prominent phase drift rate changes, measures the extreme value growth rate, analyzes the stability of harmonic distortion, extracts the time period with prominent distortion amplitude, calculates the peak time point and the corresponding distortion amplitude, establishes a phase shift trend curve, and obtains the harmonic distortion trend prediction value; The harmonic source separation module calls the harmonic distortion trend prediction value, extracts the A, B, and C phase currents in the three-phase current signal, calculates the zero-sequence component, measures the phase offset and harmonic distortion amplitude of each phase current, calculates the distortion amplitude change rate within the cycle, screens the harmonic components with prominent distortion amplitude changes within the cycle, compares the distortion distribution of the harmonic components of each phase, analyzes the harmonic source, and obtains the harmonic source separation result; The compensation current calculation module calculates the harmonic compensation currents of phases A, B, and C based on the harmonic source separation result, screens the harmonic components with prominent multi-phase influence, calculates the injection time point of the compensation current, measures the time interval before the harmonic peak arrives, calculates the compensation current parameters, measures the difference in compensation currents of different phases, and obtains asymmetric compensation parameters; The adaptive compensation control module calls the asymmetric compensation parameters, collects the three-phase current data before compensation, calculates the harmonic distortion amplitude of each frequency band before compensation, measures the priority of multiple frequency bands, compares the current amplitude changes before and after compensation, adjusts the compensation current amplitude and phase, records the compensation execution data, and obtains the adaptive compensation execution data.

[0015] Compared with the prior art, the advantages and positive effects of the present invention are: In the present invention, based on high-frequency sampling, instantaneous extreme points are accurately extracted, time intervals and change rates are calculated, and the characterization of harmonic distortion characteristics is strengthened. The extreme value growth rate and phase drift rate are used to screen time periods with prominent fluctuation amplitudes, thereby improving the accuracy of harmonic change trend prediction. By calculating the zero-sequence component and phase offset, combined with the periodic change characteristics, accurate separation of harmonic sources is achieved, and the ability to identify harmonic causes is improved. The compensation current is calculated based on the harmonic source, the injection timing is optimized, and the compensation parameters are adjusted before the harmonic peak to improve the compensation response speed and matching accuracy. The multi-band priority strategy ensures the optimization of compensation for harmonics of different frequencies, reduces system harmonic interference, and enhances the stability of power quality. BRIEF DESCRIPTION OF THE DRAWINGS

[0016] Figure 1 It is a schematic diagram of the workflow of the present invention; Figure 2 Flow chart of the steps for obtaining harmonic distortion characteristic parameters of the present invention; Figure 3 A flowchart of the steps for obtaining the harmonic distortion trend prediction value of the present invention; Figure 4 A flowchart of the steps for obtaining the harmonic source separation result of the present invention; Figure 5 is a flow chart of the steps for obtaining asymmetric compensation parameters of the present invention; Figure 6 This is a flow chart of the steps for acquiring data for adaptive compensation execution of the present invention. DETAILED DESCRIPTION

[0017] In order to make the purpose, technical solution and advantages of the present invention more clearly understood, the present invention is further described in detail below in conjunction with the accompanying drawings and embodiments. It should be understood that the specific embodiments described herein are only used to explain the present invention and are not intended to limit the present invention.

[0018] In the description of the present invention, it should be understood that the terms "length", "width", "up", "down", "front", "back", "left", "right", "vertical", "horizontal", "top", "bottom", "inside", "outside" and the like indicate positions or positional relationships based on the positions or positional relationships shown in the drawings, and are only for the convenience of describing the present invention and simplifying the description, rather than indicating or implying that the device or element referred to must have a specific orientation, be constructed and operated in a specific orientation, and therefore cannot be understood as limiting the present invention. In addition, in the description of the present invention, "multiple" means two or more, unless otherwise clearly and specifically defined.

[0019] Example 1: Please refer to Figure 1 The present invention provides a technical solution: an electrical automation adjustment method, comprising the following steps: S1: Obtain three-phase current signals, collect single-cycle high-frequency samples, extract instantaneous extreme points, record time coordinates and amplitudes, calculate time intervals and change rates, screen extreme point distributions, and obtain harmonic distortion characteristic parameters; S2: Based on the harmonic distortion characteristic parameters, call the instantaneous extreme value, calculate the phase drift rate and the extreme value growth rate, analyze the distortion stability, screen the period with prominent fluctuation amplitude, determine the dominant harmonic component, calculate the peak time and distortion amplitude, construct the trend change curve, and calculate the harmonic distortion trend prediction value; S3: Call the harmonic distortion trend prediction value, extract the A, B, and C phase currents in the three-phase current signal, calculate the zero-sequence component, obtain the phase offset and harmonic distortion amplitude, compare the periodic changes, determine the harmonic source, and obtain the harmonic source separation result; S4: Based on the harmonic source separation results, calculate the compensation currents of phases A, B, and C, compare the multi-phase impact, calculate the injection time point of the harmonic compensation current, set the compensation current parameters before the harmonic peak reaches, and obtain the asymmetric compensation parameters; S5: Call asymmetric compensation parameters, collect pre-compensation data, calculate multi-band priorities, compare amplitude changes, adjust current amplitude and phase, and record adaptive compensation execution data.

[0020] The characteristic parameters of harmonic distortion include time interval, rate of change, and extreme point distribution. The trend prediction values ​​of harmonic distortion include phase drift rate, extreme value growth rate, distortion stability, dominant harmonic component, peak time, distortion amplitude, and trend change curve. The results of harmonic source separation include zero-sequence component, phase offset, harmonic distortion amplitude, periodic change, and harmonic source. The asymmetric compensation parameters include A-phase compensation current, B-phase compensation current, C-phase compensation current, and injection time point of harmonic compensation current. The adaptive compensation execution data includes pre-compensation data, multi-band priority, amplitude change, current amplitude adjustment, and current phase adjustment.

[0021] See also Figure 2 , the specific steps for obtaining the harmonic distortion characteristic parameters are: S101: Acquire three-phase current signals, perform single-cycle high-frequency sampling, record current amplitudes corresponding to multiple time points, extract instantaneous extreme value points at multiple times, store the time coordinates and amplitudes of the extreme value points, and generate an instantaneous extreme value data set; The three-phase current signal is obtained, and a high-frequency sampling device is used to collect the current signal in a single cycle. The specific operation includes at least thousands of equal-interval samplings within a complete power frequency cycle (20ms for a 50Hz system and 16.67ms for a 60Hz system), and the current amplitude of each sampling point is recorded. The sampling interval must meet the Nyquist sampling theorem, that is, at least twice the highest frequency of the harmonic component. For example, for sampling a 2kHz harmonic component, a sampling rate of at least 4kHz is required. For each sampling point, determine whether it is a local extreme point based on the positive and negative change trend of the current signal. The specific method is: calculate the current value change trend of three adjacent time points. If the current value of the current point is greater than the two previous and next points, it is determined to be a local maximum value; if the current value of the current point is less than the two previous and next points, it is determined to be a local minimum value. The time coordinates and amplitude data of the extreme point are stored in the data set to generate an instantaneous extreme value data set, as shown in Table 1.

[0022] Table 1 Examples of instantaneous extreme value datasets

[0023] As shown in Table 1, the sampling point number represents the extreme point appearing in the time series. Its time coordinate is recorded by the high-frequency sampling device, and the amplitude is the corresponding instantaneous current value. This data set provides the basis for subsequent analysis.

[0024] The results show that the instantaneous extreme point data of the current signal was successfully obtained, and the matching relationship between the time coordinate and the amplitude was established. This data set provides a basis for subsequent harmonic distortion analysis and can be further used to calculate characteristic parameters such as the time interval and change rate of the extreme points.

[0025] S102: Based on the instantaneous extreme value data set, calculate the time intervals and amplitude change rates of adjacent extreme value points, establish time interval sequences and change rate sequences, filter outliers in the time interval sequences, and remove time interval data that exceeds the set threshold, using the formula: ; Calculate and generate the distribution deviation of the rate of change of extreme value points; in, Represents the distribution deviation of the rate of change of extreme points, Representative The rate of change of extreme points, represents the average value of the rate of change, represents the standard deviation of the rate of change, Representative The time interval between extreme points, represents the mean of the time interval, Represents the total number of extreme points; For example, the time interval between the second extreme point and the first extreme point in Table 1 is . Calculate the time intervals between all adjacent extreme points and store them in the time interval sequence. At the same time, calculate the amplitude change rate of adjacent extreme points, defined as For example, the amplitude change rate of the second extreme point in Table 1 is , repeat the calculation to obtain the complete change rate sequence. Filter the time interval sequence and remove abnormal values ​​that exceed the set threshold. For example, if the time interval threshold is set to [0.5ms, 1.5ms], if the time interval value exceeds this range, remove the data point and recalculate the time interval mean. Use the following formula to calculate the distribution deviation of the extreme point change rate: ; in, is the mean of the rate of change of all extreme points, and the calculation formula is: ; Assume that the aforementioned extreme point change rate sequence is A / ms, then its mean value is: ; Calculate the standard deviation: ; ; calculate The numerator part of the value: ; ; calculate The denominator part of the value: ; in is the mean of the time interval, assuming ,but: ; ; Final calculation: ; The result shows that the distribution deviation value of the extreme point change rate is calculated. The larger the value, the more discrete the distribution of the extreme point change rate. This value will be used for subsequent abnormal data point screening to determine whether there is any abnormal harmonic distortion.

[0026] S103: Call the extreme point change rate distribution deviation, filter out abnormal data points that exceed the set threshold, count the proportion of abnormal data points, calculate the abnormal distribution rate of data points, combine the time interval screening results, extract extreme points that meet the set distribution characteristics, obtain the extreme point distribution interval parameters, and generate harmonic distortion characteristic parameters.

[0027] If it exceeds the set threshold (such as 6.0), the current data is judged to be abnormal, and the abnormal data points are screened and the proportion of abnormal data points is counted. If the total number of extreme points in the data set is 100 and the number of abnormal data points is 7, then the data point abnormal distribution rate is: ; Combined with the time interval screening results, the extreme points that meet the set distribution characteristics are extracted. That is, when the change rate and the time interval both meet the set thresholds, the corresponding extreme points are retained, and finally the extreme point distribution interval parameters are obtained. For example, the extreme point change rate distribution range is set to [25A / ms, 40A / ms], and the time interval range is [0.5ms, 1.5ms]. The remaining data points after screening constitute the final harmonic distortion characteristic parameters.

[0028] The results show that the calculation of the abnormal distribution rate of data points is used to evaluate the fluctuation of the harmonic distortion signal and assist in determining which data points may belong to abnormal distortion data. The screened extreme point distribution interval parameters can be used for further harmonic feature extraction to determine whether there is significant harmonic interference in the power grid.

[0029] See also Figure 3 , the specific steps for obtaining the harmonic distortion trend prediction value are: S201: based on the harmonic distortion characteristic parameters, calling the instantaneous extreme value, comparing the extreme value change amplitudes in the differentiated time periods, calculating the extreme value growth rate, and screening according to the set extreme value growth rate reference value, extracting the time period with outstanding growth rate, and obtaining the extreme value growth rate screening result; First, call the instantaneous extreme value, that is, extract the peak points and valley points of each dominant harmonic component from the time series, use the numerical difference method to calculate the extreme value increase of adjacent time points, set the time window to count these changes, so as to obtain the change trend of the characteristic parameters, then compare the extreme value change amplitudes in different time periods, and calculate the extreme value growth rate. Specifically, within a certain time window, the extreme value growth rate of the dominant harmonic component is calculated as: ; in, For the The extreme value growth rate of the harmonic components, and are the harmonic extreme values ​​at adjacent time points, The extreme value growth rate in multiple time windows is calculated, and further screened according to the set extreme value growth rate benchmark value. The determination of the benchmark value can be based on the statistical analysis of historical data. For example, the average growth rate in the past 30 time windows is calculated. and standard deviation , then the reference value is set as: ; Then filter out the growth rate exceeding These time periods are the time periods with prominent extreme growth rates, as shown in Table 2, which lists a set of sample data.

[0030] Table 2 Extreme growth rate ; As shown in Table 2, based on the set benchmark value, time points 1 and 3 with higher extreme value growth rates are screened out, and the extreme value growth rate screening results are obtained.

[0031] S202: calling the extreme value growth rate screening result, calculating the corresponding phase drift rate for the screened time period, determining the fluctuation of the phase drift rate, screening the time interval for the fluctuation amplitude, and obtaining the phase drift rate fluctuation interval; It is necessary to obtain the phase change of the dominant harmonic component at adjacent time points, and calculate the phase drift rate formula: ; in, and Respectively indicate time and The phase angle at time, Represents the phase drift rate. Then, the fluctuation of the phase drift rate is determined and the phase drift rate fluctuation benchmark value is set. This is similar to the benchmark value setting of the extreme value growth rate mentioned above and is set based on statistical data: ; Filter out the time intervals with large fluctuations. For example, if the drift rate of the third harmonic within a certain period of time is They are 0.1, 0.2, 0.15, 0.4, and 0.05 respectively, with a mean of 0.18 and a standard deviation of 0.12. The benchmark value is 0.42, and the time period with a drift rate higher than this value is screened out to finally obtain the phase drift rate fluctuation range.

[0032] S203: Call the phase drift rate fluctuation interval, calculate the peak time and distortion amplitude of the dominant harmonic component in the fluctuation interval, and construct a trend change curve using the formula: ; Calculate and obtain the predicted value of harmonic distortion trend; in, represents the predicted value of harmonic distortion trend, Representative The distortion amplitude of the dominant harmonic component, Representative The attenuation factor of the dominant harmonic component, represents the time variable, Representative The drift factor of the dominant harmonic component, Representative The phase shift of the dominant harmonic component, Represents the total number of dominant harmonic components.

[0033] formula: ; In this formula, is the distortion amplitude of the dominant harmonic component, represents the attenuation factor, which can be determined experimentally, such as by exponential fitting wait, For time variables, set the calculation step size , is the drift factor, which is measured by statistical regression method. , is the phase offset, which is calculated based on the drift rate.

[0034] Assume that the distortion amplitudes of the dominant harmonic components measured in a certain period of time are , , then calculate the predicted value: ; ; ; ; This result shows that the calculated harmonic distortion trend prediction value is 2.72V, which can be used for subsequent harmonic distortion trend evaluation.

[0035] See also Figure 4 , the specific steps for obtaining the harmonic source separation results are: S301: calling the harmonic distortion trend prediction value, extracting the A, B, and C phase currents in the three-phase current signal, calculating the root mean square value of the multi-phase current signal in multiple cycles, obtaining the intra-cycle fluctuation characteristics of the A, B, and C phase currents, and performing differential operation based on the root mean square value of the multi-phase current and the zero-sequence component, calculating the current offset of the zero-sequence component, and obtaining the zero-sequence component offset rate; First, we need to obtain the A, B, and C phase currents in the three-phase current signal, and sample each phase current signal in cycles. Assume that the number of sampling points in one cycle is , then in Within the range of cycles, all sampling values ​​of each phase current are accumulated. Then, the root mean square value (RMS) of each phase current is calculated, that is: ; in, Represents phases A, B, and C. For the After obtaining the RMS values ​​of the three-phase currents A, B, and C through the above calculation, the fluctuations within the cycle range are analyzed, that is, the change in the RMS value within each cycle is calculated to represent the periodic fluctuation amplitude of the current. Assuming the cycle The RMS value is , the fluctuation range is calculated as follows: ; For the zero-sequence component, its calculation is based on the sum of the instantaneous values ​​of the three-phase currents: ; Then, the RMS value of the zero-sequence component within the cycle is calculated and a differential operation is performed: ; In order to further analyze the changing trend of the zero-sequence component, it is necessary to calculate its relative offset rate and set the zero-sequence component offset rate The calculation method is as follows: ; The offset rate is used to determine the current fluctuation of the zero-sequence component. If it exceeds the set threshold (for example, 0.05), it indicates that the zero-sequence component fluctuates significantly during the period. As shown in Table 1, the three-phase current RMS values ​​and zero-sequence component offset rates in a certain measurement are listed.

[0036] Table 3 Three-phase current RMS value and zero-sequence component offset rate

[0037] As shown in Table 3, the zero-sequence component offset rate of period 3 reaches 0.071, which is higher than the set threshold of 0.05. Therefore, there may be significant harmonic distortion fluctuations in this period.

[0038] S302: Based on the zero-sequence component offset rate, the phase offsets of the A, B, and C phase currents relative to the reference phase are calculated, the A and B phase currents are selected for difference calculation, and the ratio relationship of the zero-sequence component is constructed using the formula: ; Calculate the phase change trend of A-phase current and B-phase current, obtain the phase difference parameter of A-phase current and B-phase current, determine the comparison relationship between phase change and zero-sequence component, screen the key phase distortion characteristics, and obtain the main phase distortion amount; in, represents the phase offset, Representative The current signal of phase A and phase B in one cycle, The current signal representing the zero-sequence component, Represents the current difference between phases A and B, represents the maximum current during the measurement, represents the number of cycles selected during the calculation process, Represents the inverse tangent function operation, which is used to calculate the phase angle corresponding to the ratio of the current signal; Based on the zero-sequence component offset rate calculated in the previous paragraph , calculate the phase offset of the three-phase currents A, B, and C relative to the reference phase. Assume that phase A is the reference phase, and calculate the phase offset of phases B and C relative to phase A. Assume that in the cycle The reference phase of phase A current is , the instantaneous phase of phase B current is , then the phase offset is calculated as follows: ; The differential calculation method is used to calculate the difference ratio of the A and B phase currents, and the ratio relationship of the zero-sequence component is constructed: ; in, and Phase A and phase B are The current value of each cycle, is the zero-sequence component current signal, Indicates the current difference between phases A and B, calculated as follows: ; ; set up The measured currents of phases A and B are as follows: ; Calculate the difference: ; Assuming the zero-sequence current , maximum current , calculate the phase offset: ; ; ; This calculation shows that the phase shift of the B phase current relative to the A phase is about 54.7°, indicating a significant phase difference that can be used to determine the source of harmonic distortion.

[0039] S303: Call the main phase distortion variable, calculate the corresponding harmonic distortion amplitude, and based on the change trend of the harmonic distortion amplitude in multiple cycles, select the extreme value points of the harmonic amplitude to calculate the fluctuation range, and compare the periodic change characteristics to determine the type and contribution of the harmonic source, and obtain the harmonic source separation result.

[0040] Fourier transform is used to calculate the harmonic components, obtain the amplitude of each harmonic, and draw the harmonic distortion amplitude change curve within multiple cycles. Through the peak detection method, the extreme point is selected to calculate the fluctuation range: ; in, and are the maximum and minimum harmonic amplitudes respectively. Combined with the fluctuation characteristics of multiple cycles, the periodic variation trend of the harmonic amplitude is calculated. Suppose the 5th harmonic amplitude of a certain measurement is: ; The fluctuation range is: ; If the harmonic source determination threshold is set to 1.0, the harmonic source contributes more and its source and impact can be further analyzed.

[0041] See also Figure 5 , the specific steps for obtaining the asymmetric compensation parameters are: S401: Based on the harmonic source separation result, calculate the compensation current of phases A, B, and C, extract the multi-phase harmonic current components, call the harmonic current component data, calculate the peak amplitude of each phase harmonic current, and calculate the phase change trend based on the peak amplitude data. For the three-phase harmonic current data, compare the phase deviation, calculate the phase offset value between the three phases, and obtain the three-phase harmonic phase offset; Call the three-phase harmonic current time series data, analyze the time domain distribution of each phase harmonic current, and calculate the peak amplitude of each phase harmonic current. The peak amplitude is calculated by selecting the maximum positive and negative values ​​of the harmonic current in each cycle, taking their absolute values ​​for average calculation. For example, the maximum value of the harmonic current of phase A in a certain cycle is 3.2A, and the minimum value is -3.0A, then the peak amplitude is calculated as , based on the peak amplitude data in all cycles, calculate the average peak amplitude of each phase, and call the three-phase harmonic current phase data at the same time to calculate the phase change trend, use the difference to calculate the phase change rate at adjacent moments, and compare the phase change data of each phase to obtain the phase offset value between the three phases. The phase offset value between the three phases is calculated based on the phase difference between adjacent phases. For example, the phase difference between phase B and phase A is , the phase difference between phase C and phase A is Finally, the three-phase harmonic phase offset is obtained. Assuming that the harmonic current phases of phases A, B, and C are , , , then the phase offsets are , , and finally a data set of three-phase harmonic phase offsets is formed.

[0042] S402: calling the three-phase harmonic phase offset, calculating the time offset value of the multi-phase compensation current, calling the three-phase compensation current time offset data, comparing the harmonic compensation impact, calculating the timing difference of the compensation current, obtaining the optimal injection time point of the compensation current, and calling the harmonic amplitude data at the same time, calculating the compensation amplitude, and establishing the optimal compensation current parameters; Calculate the time offset value of the three-phase compensation current. The calculation method of the time offset value is based on the relationship between the phase offset and the power frequency cycle. For example, if the power frequency cycle is 20ms, the phase offset The corresponding time offset is calculated as , similarly The corresponding time offset value is , combined with the time offset data of the three-phase compensation current, compared with the harmonic compensation effect, the timing difference of the compensation current is calculated. The calculation of the timing difference is based on the time alignment degree of the compensation current phase and the actual harmonic current phase. If the compensation current phase lags behind the harmonic current, the compensation current needs to be injected in advance. If it is ahead, the compensation current needs to be injected later. Based on this adjustment strategy, the optimal injection time point of the compensation current is obtained. At the same time, the harmonic amplitude data is called to calculate the compensation amplitude. The calculation method is: Assume that the amplitude of the harmonic current is , the compensation current amplitude is ,but ,in is the compensation coefficient, assuming , compensation coefficient , then the compensation current amplitude , and finally establish the optimal compensation current parameters.

[0043] S403: Based on the optimal compensation current parameters, the compensation current adjustment parameters before the harmonic peak value is reached are set, and the parameters of the multi-phase compensation current are calculated using the formula: ; Obtain asymmetric compensation parameters; in, Represents the overall parameter of the compensation current, Representative Order harmonic voltage, Representative The phase angle of the harmonic voltage of order, represents the phase adjustment parameter, represents the total number of harmonic orders, Is an imaginary unit.

[0044] formula: ; in, Represents the overall parameter of the compensation current, Representative Order harmonic voltage, Representative The phase angle of the harmonic voltage of order, represents the phase adjustment parameter, represents the total number of harmonic orders, is an imaginary unit. Calculation is based on known data, assuming that the total number of harmonic orders is , the harmonic voltage amplitudes of each order are , , , and their phase angles are , , The phase adjustment parameters are , , , put it into the formula: ; Further calculation: ; Convert exponential form to rectangular coordinates: ; Calculate three values: ; ; ; Calculate the modulus: ; Finally, the asymmetric compensation parameters are obtained .

[0045] Table 4 Compensation current calculation parameters

[0046] As shown in Table 4, the voltage values ​​and phase angle information of different harmonic orders, as well as the phase adjustment parameters are listed. The final calculated compensation current parameters are: , whose modulus is ,The result shows that the calculated asymmetric compensation parameters meet ,the requirements of harmonic current compensation and can be used in the actual ,compensation control link.

[0047] See also Figure 6 , the specific steps for obtaining adaptive compensation execution data are: S501: calling asymmetric compensation parameters, collecting data before compensation, calculating signal amplitudes and spectrum characteristics of multiple frequency bands, obtaining signal change rates of differentiated frequency bands, and filtering frequency bands based on the change rates, calculating weights of frequency bands, and sorting them to obtain multi-band priority sorting; First, set the signal sampling points of multiple frequency bands, scan in the range of 1MHz to 100MHz with a step frequency of 1kHz, and collect the amplitude data of the signal , and stored in the matrix In the process, the reference signal amplitude is recorded at the same time The initial value of the signal amplitude is set to 1V. Then, the signal amplitude change rate of each frequency band is calculated using The signal variation of each frequency band is obtained by V determines that the frequency band has significant changes, otherwise the data of the frequency band is ignored. Subsequently, the frequency bands with significant changes are screened and the weights of each frequency band are calculated. , the weight calculation adopts the normalization method, that is: ; To ensure that the sum of the weights of all frequency bands is 1. Then sort the calculated weights in descending order so that the frequency bands with larger changes have higher priority, and finally form a multi-band priority sorting list.

[0048] Table 5 Multi-band signal change rate calculation table

[0049] As shown in Table 5, the signal change rates of different frequency bands are calculated, among which the 20MHz frequency band has the largest change, so it ranks highest in the priority sorting.

[0050] S502: Based on the priority sorting of multiple frequency bands, the signal amplitude change of the selected frequency band is calculated, the initial current amplitude and phase parameters are called, the amplitude change ratio of the multi-frequency band signal is calculated, and compared with the set reference change threshold, using the formula: ; Calculate the current adjustment factor, adjust the current amplitude and phase of multiple frequency bands, and obtain the current correction parameters; in, represents the current adjustment factor, Representative The current signal amplitude of the frequency band, represents the reference signal amplitude, Representative The weight of the frequency band, Represents the number of frequency bands to be filtered, , represents the operation of capturing the maximum and minimum values ​​in all filtered frequency bands; Before selection frequency bands and calculate the ratio of the signal amplitude change. First, call the initial current amplitude A and Phase Parameters , calculate the current amplitude adjustment factor for the selected frequency band The calculation method is as follows: 1. Get the amplitude variation range of the selected frequency band, that is and , for example, for the selected frequency band {5MHz, 10MHz, 20MHz}, the maximum change is 0.20V and the minimum change is 0.08V; 2. Calculate the difference in signal change amplitude: ; 3. Calculate the weighted average weight: ; 4. Calculate the current adjustment factor: ; The calculation results show that the current amplitude and phase need to be adjusted to compensate for the signal changes.

[0051] S503: calling the current correction parameter, adjusting the actual applied current amplitude and phase, and recording the adjusted signal response in real time, calculating the difference of the signal offset before and after the adjustment, and obtaining the adaptive compensation execution data.

[0052] Call the calculated current correction parameters A, adjust the actual applied current, that is, the new current amplitude: ; At the same time, the phase adjustment amount is determined according to the correction parameter, and the new phase is set to ,in It depends on the phase compensation model of the specific system. This setting is ,Right now The adjusted signal response is recorded in real time, the signal offset before and after the adjustment is calculated, and finally the adaptive compensation execution data is obtained.

[0053] An electrical automation regulation system, the electrical automation regulation system is used to execute the above electrical automation regulation method, the system comprises: The harmonic distortion feature extraction module obtains three-phase current signals, extracts single-cycle high-frequency sampling data, screens instantaneous extreme value points among the sampling points, records the time coordinates and amplitudes of the instantaneous extreme value points, calculates the time intervals between consecutive instantaneous extreme value points, measures the extreme value change rate, screens out extreme value points with prominent fluctuation amplitudes of the change rate, analyzes the distribution density of extreme value points, extracts distortion frequency characteristics, calculates the relative amplitudes of the fundamental wave and each harmonic, measures the total harmonic distortion rate, establishes a distortion trend curve, and obtains harmonic distortion feature parameters; The phase shift calculation module is based on the harmonic distortion characteristic parameters, calls the instantaneous extreme point data, calculates the phase drift rate of each harmonic, screens the harmonic components with prominent phase drift rate changes, measures the extreme value growth rate, analyzes the stability of harmonic distortion, extracts the time period with prominent distortion amplitude, calculates the peak time point and the corresponding distortion amplitude, establishes the phase shift trend curve, and obtains the harmonic distortion trend prediction value; The harmonic source separation module calls the harmonic distortion trend prediction value, extracts the A, B, and C phase currents in the three-phase current signal, calculates the zero-sequence component, measures the phase offset and harmonic distortion amplitude of each phase current, calculates the distortion amplitude change rate within the cycle, screens the harmonic components with prominent distortion amplitude changes within the cycle, compares the distortion distribution of each phase harmonic component, analyzes the harmonic source, and obtains the harmonic source separation result; The compensation current calculation module calculates the harmonic compensation currents of phases A, B, and C based on the results of harmonic source separation, screens the harmonic components with prominent multi-phase influence, calculates the injection time point of the compensation current, measures the time interval before the harmonic peak arrives, calculates the compensation current parameters, measures the difference in compensation currents of different phases, and obtains the asymmetric compensation parameters; The adaptive compensation control module calls the asymmetric compensation parameters, collects the three-phase current data before compensation, calculates the harmonic distortion amplitude of each frequency band before compensation, measures the priority of multiple frequency bands, compares the current amplitude changes before and after compensation, adjusts the compensation current amplitude and phase, records the compensation execution data, and obtains the adaptive compensation execution data.

[0054] The above are only preferred embodiments of the present invention and are not intended to limit the present invention in other forms. Any technician familiar with the profession may use the technical contents disclosed above to change or modify them into equivalent embodiments with equivalent changes and apply them to other fields. However, any simple modification, equivalent change and modification made to the above embodiments based on the technical essence of the present invention without departing from the technical solution of the present invention still falls within the protection scope of the technical solution of the present invention.

Claims

1. An electrical automation adjustment method, characterized in that: The following steps are involved: S1: Acquire three-phase current signals, collect single-cycle high-frequency samples, extract instantaneous extreme points, record time coordinates and amplitudes, calculate time intervals and change rates, screen extreme point distributions, and obtain harmonic distortion characteristic parameters; S2: Based on the harmonic distortion characteristic parameters, call the instantaneous extreme value, calculate the phase drift rate and the extreme value growth rate, analyze the distortion stability, screen the time period with prominent fluctuation amplitude, determine the dominant harmonic component, measure the peak time and distortion amplitude, construct a trend change curve, and calculate the harmonic distortion trend prediction value; S3: Calling the harmonic distortion trend prediction value, extracting the A, B, and C phase currents in the three-phase current signal, calculating the zero-sequence component, obtaining the phase offset and harmonic distortion amplitude, comparing the periodic changes, determining the harmonic source, and obtaining the harmonic source separation result; S4: Based on the harmonic source separation result, calculate the compensation currents of phases A, B, and C, compare the multi-phase impact, calculate the injection time point of the harmonic compensation current, set the compensation current parameters before the harmonic peak arrives, and obtain the asymmetric compensation parameters.

2. The electrical automation adjustment method according to claim 1, characterized in that: The harmonic distortion characteristic parameters include time interval, rate of change, and extreme point distribution; the harmonic distortion trend prediction value includes phase drift rate, extreme value growth rate, distortion stability, dominant harmonic component, peak time, distortion amplitude, and trend change curve; the harmonic source separation result includes zero-sequence component, phase offset, harmonic distortion amplitude, periodic change, and harmonic source; the asymmetric compensation parameters include A-phase compensation current, B-phase compensation current, C-phase compensation current, and injection time point of harmonic compensation current.

3. The electrical automation adjustment method according to claim 2, characterized in that: The steps for obtaining the harmonic distortion characteristic parameters are specifically as follows: S101: Acquire three-phase current signals, perform single-cycle high-frequency sampling, record current amplitudes corresponding to multiple time points, extract instantaneous extreme value points at multiple moments, store the time coordinates and amplitudes of the extreme value points, and generate an instantaneous extreme value data set; S102: Based on the instantaneous extreme value data set, calculate the time intervals and amplitude change rates of adjacent extreme value points, establish a time interval sequence and a change rate sequence, filter outliers in the time interval sequence, and remove time interval data that exceeds a set threshold, using the formula: ； Calculate and generate the distribution deviation of the rate of change of extreme points; in, Represents the distribution deviation of the rate of change of extreme points, Representative The rate of change of extreme points, represents the average value of the rate of change, represents the standard deviation of the rate of change, Representative The time interval between extreme points, represents the mean of the time interval, Represents the total number of extreme points; S103: Call the extreme point change rate distribution deviation, filter out abnormal data points that exceed the set threshold, count the proportion of abnormal data points, calculate the abnormal distribution rate of data points, combine the time interval screening results, extract extreme points that meet the set distribution characteristics, obtain extreme point distribution interval parameters, and generate harmonic distortion characteristic parameters.

4. The electrical automation adjustment method according to claim 3, characterized in that: The steps for obtaining the harmonic distortion trend prediction value are specifically as follows: S201: Based on the harmonic distortion characteristic parameters, call the instantaneous extreme value, compare the extreme value change amplitudes in the differentiated time periods, calculate the extreme value growth rate, and screen according to the set extreme value growth rate benchmark value to extract the time period with the outstanding growth rate to obtain the extreme value growth rate screening result; S202: Calling the extreme value growth rate screening result, calculating the corresponding phase drift rate for the screened time period, determining the fluctuation of the phase drift rate, screening the time interval based on the fluctuation amplitude, and obtaining the phase drift rate fluctuation interval; S203: Call the phase drift rate fluctuation range, calculate the peak time and distortion amplitude of the dominant harmonic component in the fluctuation range, and construct a trend change curve using the formula: ； Calculate and obtain the predicted value of harmonic distortion trend; in, represents the predicted value of harmonic distortion trend, Representative The distortion amplitude of the dominant harmonic component, Representative The attenuation factor of the dominant harmonic component, represents the time variable, Representative The drift factor of the dominant harmonic component, Representative The phase shift of the dominant harmonic component, Represents the total number of dominant harmonic components.

5. The electrical automation adjustment method according to claim 4, characterized in that: The steps for obtaining the harmonic source separation result are specifically as follows: S301: calling the harmonic distortion trend prediction value, extracting the A, B, and C phase currents from the three-phase current signal, calculating the root mean square value of the multi-phase current signal within multiple cycles, obtaining the intra-cycle fluctuation characteristics of the A, B, and C phase currents, and performing a differential operation based on the root mean square value of the multi-phase current and the zero-sequence component to calculate the current offset of the zero-sequence component and obtain the zero-sequence component offset rate; S302: Based on the zero-sequence component offset rate, calculate the phase offsets of the A, B, and C phase currents relative to the reference phase, select the A and B phase currents for difference calculation, and construct the ratio relationship of the zero-sequence component using the formula: ； Calculate the phase change trend of the A and B phase currents, obtain the phase difference parameters of the A and B phase currents, determine the comparison relationship between the phase change and the zero-sequence component, screen the key phase distortion characteristics, and obtain the main phase distortion; in, represents the phase offset, Representative The current signal of phase A and phase B in a cycle, The current signal representing the zero-sequence component, Represents the current difference between phases A and B, represents the maximum current during the measurement period, represents the number of cycles selected during the calculation process, Represents the inverse tangent function operation, which is used to calculate the phase angle corresponding to the ratio of current signals; S303: Call the main phase distortion value and calculate the corresponding harmonic distortion amplitude. At the same time, based on the change trend of the harmonic distortion amplitude in multiple cycles, select the extreme value point of the harmonic amplitude to calculate the fluctuation range, and compare the periodic change characteristics to determine the type and contribution of the harmonic source to obtain the harmonic source separation result.

6. The electrical automation adjustment method according to claim 5, characterized in that: The steps for obtaining the asymmetric compensation parameters are specifically as follows: S401: Based on the harmonic source separation result, calculate the compensation currents of phases A, B, and C, extract the multi-phase harmonic current components, call the harmonic current component data, calculate the peak amplitude of the harmonic current of each phase, and calculate the phase change trend based on the peak amplitude data. Compare the phase deviations of the three-phase harmonic current data, calculate the phase offset values ​​between the three phases, and obtain the three-phase harmonic phase offset; S402: Calling the three-phase harmonic phase offset, calculating the time offset value of the multi-phase compensation current, calling the three-phase compensation current time offset data, comparing the harmonic compensation impact, calculating the timing difference of the compensation current, obtaining the optimal injection time point of the compensation current, and calling the harmonic amplitude data at the same time, calculating the compensation amplitude, and establishing the optimal compensation current parameters; S403: Based on the optimal compensation current parameters, compensation current adjustment parameters are set before the harmonic peak reaches the target, and multi-phase compensation current parameters are calculated using the formula: ； Obtain asymmetric compensation parameters; in, Represents the overall parameter of the compensation current, Representative Order harmonic voltage, Representative The phase angle of the order harmonic voltage, represents the phase adjustment parameter, represents the total number of harmonic orders, Is an imaginary unit.

7. The electrical automation adjustment method according to claim 6, characterized in that: The method further comprises: S5: calling the asymmetric compensation parameters, collecting pre-compensation data, calculating multi-band priorities, comparing amplitude changes, adjusting current amplitude and phase, and recording adaptive compensation execution data; The adaptive compensation execution data includes pre-compensation data, multi-band priorities, amplitude change, current amplitude adjustment, and current phase adjustment.

8. The electrical automation adjustment method according to claim 7, characterized in that: The steps for acquiring the adaptive compensation execution data are specifically as follows: S501: Calling the asymmetric compensation parameters, collecting pre-compensation data, calculating the signal amplitudes and spectral characteristics of multiple frequency bands, obtaining signal change rates of differentiated frequency bands, filtering frequency bands based on the change rates, calculating weights of the frequency bands, and sorting the frequency bands to obtain a priority ranking of the multiple frequency bands; S502: Based on the multi-band priority ranking, the signal amplitude change of the selected frequency band is calculated, the initial current amplitude and phase parameters are called, the amplitude change ratio of the multi-band signal is calculated, and compared with the set reference change threshold, using the formula: ； Calculate the current adjustment factor, adjust the current amplitude and phase of multiple frequency bands, and obtain the current correction parameters; in, represents the current adjustment factor, Representative The current signal amplitude of the frequency band, represents the reference signal amplitude, Representative The weight of the frequency band, Represents the number of frequency bands to be filtered, 、 Represents the operation of capturing the maximum and minimum values ​​in all filtered frequency bands; S503: calling the current correction parameter, adjusting the actual applied current amplitude and phase, and recording the adjusted signal response in real time, calculating the difference between the signal offsets before and after the adjustment, and obtaining adaptive compensation execution data.

9. An electrical automation regulation system, characterized in that: According to any one of claims 1 to 8, the electrical automation regulation method comprises: The harmonic distortion feature extraction module acquires three-phase current signals, extracts single-cycle high-frequency sampling data, screens instantaneous extreme points among the sampling points, records the time coordinates and amplitudes of the instantaneous extreme points, calculates the time intervals between consecutive instantaneous extreme points, measures the extreme value change rate, screens out extreme points with prominent fluctuation amplitudes in the change rate, analyzes the distribution density of the extreme points, extracts distortion frequency characteristics, calculates the relative amplitudes of the fundamental wave and each harmonic, measures the total harmonic distortion rate, establishes a distortion trend curve, and obtains harmonic distortion feature parameters. The phase shift calculation module calls the instantaneous extreme point data based on the harmonic distortion characteristic parameters, calculates the phase drift rate of each harmonic, screens the harmonic components with prominent phase drift rate changes, measures the extreme value growth rate, analyzes the stability of the harmonic distortion, extracts the time period with prominent distortion amplitude, calculates the peak time point and the corresponding distortion amplitude, establishes a phase shift trend curve, and obtains the harmonic distortion trend prediction value; The harmonic source separation module calls the harmonic distortion trend prediction value, extracts the A, B, and C phase currents in the three-phase current signal, calculates the zero-sequence component, measures the phase offset and harmonic distortion amplitude of each phase current, calculates the distortion amplitude change rate within the cycle, screens the harmonic components with prominent distortion amplitude changes within the cycle, compares the distortion distribution of the harmonic components of each phase, analyzes the harmonic source, and obtains the harmonic source separation result; The compensation current calculation module calculates the harmonic compensation currents of phases A, B, and C based on the harmonic source separation results, screens the harmonic components with prominent multi-phase effects, calculates the injection time point of the compensation current, measures the time interval before the harmonic peak arrives, calculates the compensation current parameters, measures the difference in compensation currents of different phases, and obtains asymmetric compensation parameters; The adaptive compensation control module calls the asymmetric compensation parameters, collects the three-phase current data before compensation, calculates the harmonic distortion amplitude of each frequency band before compensation, measures the priority of multiple frequency bands, compares the current amplitude changes before and after compensation, adjusts the compensation current amplitude and phase, records the compensation execution data, and obtains the adaptive compensation execution data.

Images