An electrical automation regulation method and system

Abstract

The present invention relates to the field of adaptive control technology, and specifically to an electrical automation regulation method and system, which includes the following steps: obtaining three-phase current signals, collecting single-cycle high-frequency sampling, extracting instantaneous extreme points, recording time coordinates and amplitudes, calculating time intervals and change rates, screening the distribution of extreme points, and obtaining harmonic distortion characteristic parameters. In the present invention, based on high-frequency sampling, instantaneous extreme points are extracted, time intervals and change rates are calculated, the description of harmonic distortion characteristics is strengthened, the extreme value growth rate and phase drift rate are used to screen the periods with prominent fluctuation amplitudes, the prediction accuracy of harmonic change trends is improved, through the calculation of zero-sequence components and phase offset amounts, combined with the periodic change characteristics, precise separation of harmonic sources is achieved, the ability to identify the causes of harmonics is enhanced, the compensation current is calculated based on the harmonic sources, the injection timing is optimized, the compensation parameters are adjusted before the harmonic peak, and the compensation response speed and matching accuracy are improved.

Description

Technical Field

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

[0002] The technical field of adaptive control includes technical methods for adjusting control parameters or strategies during the operation of a system to enable the control system to maintain desired performance under different working conditions in the face of environmental changes, system dynamic characteristic changes, etc. 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 applied in industrial automation, robot control, power systems, aerospace, etc. 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 internal system parameters, improving control accuracy and stability.

[0003] Among them, the electrical automation regulation method refers to a technical method in an electrical control system that dynamically adjusts control parameters by real-time monitoring the system operation state according to set regulation rules to ensure that the system operates according to the expected requirements. This patent theme involves real-time regulation strategies based on feedback signals, including analyzing the input signals of electrical equipment using adaptive control algorithms and dynamically adjusting voltage, current, or power factor by calculating specific regulation parameters. In addition, this method also includes selecting appropriate control modes according to the system operation state and maintaining stable operation of the system under different load or environmental changes by adjusting the control parameters of the drive circuit, power electronic device, or programmable controller.

[0004] The prior art is based on fixed-time window mean analysis or Fourier transform, making it difficult to accurately identify instantaneous extreme points, with limited accuracy in harmonic distortion feature extraction and insufficient ability to identify sudden distortions. The identification of harmonic sources relies on static parameter matching, making it difficult to adapt to dynamic working conditions and prone to misjudgment when the load complexity increases. The compensation control uses fixed compensation coefficients, lacking adaptability to real-time working conditions, with inflexible adjustment of current amplitude and phase, and limited harmonic suppression effect, affecting the stable operation of the system. Summary of the Invention

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

[0006] To achieve the above purpose, the present invention adopts the following technical solution: An electrical automation regulation method, comprising the following steps:

[0007] S1: Obtain three-phase current signals, collect single-cycle high-frequency sampling, extract instantaneous extreme points, record the time coordinates and amplitudes, calculate the time intervals and change rates, screen the distribution of extreme points, and obtain harmonic distortion characteristic parameters;

[0008] S2: Based on the harmonic distortion characteristic parameters, call the instantaneous extreme values, calculate the phase drift rate and the extreme value growth rate, analyze the distortion stability, screen the time periods with prominent fluctuation amplitudes, determine the dominant harmonic components, measure the peak time and the distortion amplitude, construct a trend change curve, and calculate the harmonic distortion trend prediction value;

[0009] S3: Call the harmonic distortion trend prediction value, extract the A, B, and C phase currents in the three-phase current signals, calculate the zero-sequence component, obtain the phase offset and the harmonic distortion amplitude, compare the periodic changes, determine the harmonic source, and obtain the harmonic source separation result;

[0010] S4: Based on the harmonic source separation result, calculate the A, B, and C phase compensation currents, compare the multi-phase effects, calculate the injection time points of the harmonic compensation current, and set the compensation current parameters before the harmonic peak arrives to obtain the asymmetric compensation parameters.

[0011] As a further solution of the present invention, the harmonic distortion characteristic parameters include time intervals, change rates, and the distribution of extreme points. The harmonic distortion trend prediction value includes phase drift rate, extreme value growth rate, distortion stability, dominant harmonic components, 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 the A-phase compensation current, B-phase compensation current, C-phase compensation current, and the injection time point of the harmonic compensation current.

[0012] As a further solution of the present invention, the steps for obtaining the harmonic distortion characteristic parameters are specifically as follows:

[0013] S101: Obtain three-phase current signals, perform single-cycle high-frequency sampling, record the current amplitudes corresponding to multiple time points, extract the instantaneous extreme points at multiple moments, store the time coordinates and amplitudes of the extreme points, and generate an instantaneous extreme value data set;

[0014] S102: Based on the instantaneous extreme value data set, calculate the time intervals and amplitude change rates between adjacent extreme points, establish a time interval sequence and a change rate sequence, screen the outliers in the time interval sequence, eliminate the time interval data exceeding the set threshold, and use the formula:

[0015] ;

[0016] Calculate and generate the distribution deviation of the extreme point change rate;

[0017] Where, Represents the distribution deviation of the change rate of extreme points, Represents the change rate of the th extreme point, Represents the average value of the change rate, Represents the time interval of the th extreme point, Represents the average value of the time interval,

[0018] S103: Invoke the distribution deviation of the change rate of extreme points, screen abnormal data points exceeding 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.

[0019] As a further solution of the present invention, the steps for obtaining the harmonic distortion trend prediction value are specifically as follows:

[0020] S201: Based on the harmonic distortion characteristic parameters, invoke the instantaneous extreme values, compare the extreme value change amplitudes in different time periods, calculate the extreme value growth rate, and screen according to the set extreme value growth rate reference value, extract the time periods with prominent growth rates, and obtain the extreme value growth rate screening results;

[0021] S202: Invoke the extreme value growth rate screening results, calculate the corresponding phase drift rate for the screened time periods, judge the fluctuation of the phase drift rate, screen the time intervals according to the fluctuation amplitude, and obtain the phase drift rate fluctuation interval;

[0022] S203: Invoke the phase drift rate fluctuation interval, measure the peak time and distortion amplitude of the dominant harmonic component within the fluctuation interval, and construct a trend change curve, using the formula:

[0023] ;

[0024] Calculate to obtain the harmonic distortion trend prediction value;

[0025] Wherein, Represents the harmonic distortion trend prediction value, Represents the distortion amplitude of the th dominant harmonic component, Represents the attenuation factor of the th dominant harmonic component, Represents the time variable, Represents the The phase offset of the dominant harmonic component, represents the total number of dominant harmonic components.

[0026] As a further solution of the present invention, the step of obtaining the harmonic source separation result is specifically as follows:

[0027] S301: Call the predicted value of the harmonic distortion trend, extract the currents of phases A, B, and C in the three-phase current signal, calculate the root mean square value of the polyphase current signal in multiple cycles, obtain the in-cycle fluctuation characteristics of the currents of phases A, B, and C, and at the same time perform a difference operation based on the root mean square value of the polyphase current and the zero-sequence component, calculate the current offset of the zero-sequence component, and obtain the zero-sequence component offset rate;

[0028] S302: Based on the zero-sequence component offset rate, calculate the phase offset of the currents of phases A, B, and C relative to the reference phase, select the currents of phases A and B for difference calculation, construct the ratio relationship of the zero-sequence component, and use the formula:

[0029] ;

[0030] Calculate the phase change trend of the currents of phases A and B, obtain the phase difference parameter of the currents of phases A and B, judge the comparison relationship between the phase change and the zero-sequence component, screen the key phase distortion characteristics, and obtain the main phase distortion amount;

[0031] Wherein, represents the phase offset, represents the th cycle current signals of phases A and B, represents the current signal of 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, represents the arctangent function operation for calculating the phase angle corresponding to the current signal ratio;

[0032] S303: Call the main phase distortion amount, calculate the corresponding harmonic distortion amplitude, and at the same time, based on the change trend of the harmonic distortion amplitude in multiple cycles, select the harmonic amplitude extreme points for fluctuation range calculation, and compare the periodic change characteristics to judge the type and contribution degree of the harmonic source, and obtain the harmonic source separation result.

[0033] As a further solution of the present invention, the step of obtaining the asymmetric compensation parameter is specifically as follows:

[0034] S401: Based on the harmonic source separation result, calculate the compensation currents of phases A, B, and C, extract the polyphase harmonic current components, call the harmonic current component data, calculate the peak amplitude of the harmonic current of each phase, and at the same time, based on the peak amplitude data, calculate the phase change trend. 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;

[0035] S402: Call the three-phase harmonic phase offset to calculate the time offset value of the polyphase compensation current, call the three-phase compensation current time offset data, compare the harmonic compensation influence, calculate the timing difference of the compensation current, obtain the optimal injection time point of the compensation current, and at the same time, call the harmonic amplitude data to calculate the compensation amplitude and establish the optimal compensation current parameters;

[0036] S403: Based on the optimal compensation current parameters, set the compensation current adjustment parameters before the harmonic peak arrives, calculate the parameters of the polyphase compensation current, and use the formula:

[0037] ;

[0038] Obtain the asymmetric compensation parameters;

[0039] where, represents the overall parameter of the compensation current, represents the th-order harmonic voltage, represents the phase angle of the th-order harmonic voltage, represents the phase adjustment parameter, represents the total number of harmonic orders,

[0040] As a further solution of the present invention, the method further includes:

[0041] S5: Call the asymmetric compensation parameters, collect the data before compensation, calculate the multi-band priorities, compare the amplitude changes, adjust the current amplitude and phase, and record the adaptive compensation execution data;

[0042] The adaptive compensation execution data includes the data before compensation, the multi-band priorities, the amplitude changes, the current amplitude adjustment, and the current phase adjustment.

[0043] As a further solution of the present invention, the step of obtaining the adaptive compensation execution data is specifically:

[0044] S501: Call the asymmetric compensation parameters, collect the data before compensation, calculate the signal amplitudes and spectral characteristics of multiple bands, obtain the signal change rates of different bands, and at the same time, based on the change rates, screen the bands, calculate the weights of the bands, and sort them to obtain the multi-band priority sorting;

[0045] S502: Based on the multi - band priority sorting, calculate the change in the signal amplitude of the selected band, call the initial current amplitude and phase parameters, calculate the amplitude change ratio of the multi - band signal, and compare it with the set reference change threshold. Use the formula:

[0046] ;

[0047] Calculate the current adjustment factor, adjust the current amplitude and phase of multiple bands, and obtain the current correction parameter;

[0048] Wherein, represents the current adjustment factor, represents the current signal amplitude of the th band, represents the reference signal amplitude, represents the weight of the th band, represents the number of selected bands, 、 represent the operations of capturing the maximum and minimum values among all selected bands;

[0049] S503: Call the current correction parameter, adjust the actually applied current amplitude and phase, and record the adjusted signal response in real - time. Calculate the difference in the signal offset before and after adjustment to obtain the adaptive compensation execution data.

[0050] An electrical automation regulation system, which is used to execute the above - mentioned electrical automation regulation method. The system includes:

[0051] The harmonic distortion feature extraction module acquires three - phase current signals, extracts single - cycle high - frequency sampling data, filters the instantaneous extreme points in the sampling points, records the time coordinates and amplitudes of the instantaneous extreme points, calculates the time interval between consecutive instantaneous extreme points, measures the extreme value change rate, filters out the extreme points with prominent fluctuation amplitudes of the change rate, analyzes the distribution density of the extreme points, extracts the 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 the harmonic distortion feature parameters;

[0052] The phase offset calculation module, based on the harmonic distortion feature parameters, calls the instantaneous extreme point data, calculates the phase drift rate of each harmonic, filters out the harmonic components with prominent changes in the phase drift rate, measures the extreme value growth rate, analyzes the stability of the harmonic distortion, extracts the time periods with prominent distortion amplitudes, calculates the peak time points and the corresponding distortion amplitudes, and establishes a phase offset trend curve to obtain the harmonic distortion trend prediction value;

[0053] The harmonic source separation module calls the predicted value of the harmonic distortion trend, extracts the phase A, B, and C 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 change rate of the distortion amplitude within a period, screens out the harmonic components with prominent changes in the distortion amplitude within the period, compares the distortion distribution of each phase harmonic component, analyzes the harmonic source, and obtains the harmonic source separation result;

[0054] Based on the harmonic source separation result, the compensation current calculation module calculates the harmonic compensation currents for phases A, B, and C, screens out 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, and measures the differences in the compensation currents of different phases to obtain the asymmetric compensation parameters;

[0055] 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 multi-frequency band priority, compares the change in the current amplitude before and after compensation, adjusts the amplitude and phase of the compensation current, records the compensation execution data, and obtains the adaptive compensation execution data.

[0056] Compared with the prior art, the advantages and positive effects of the present invention are as follows:

[0057] In the present invention, based on high-frequency sampling, the instantaneous extreme points are accurately extracted, the time interval and change rate are calculated, and the harmonic distortion characteristics are strengthened. By using the extreme value growth rate and phase drift rate, the time periods with prominent fluctuation amplitudes are screened out to improve the prediction accuracy of the harmonic change trend. Through the calculation of the zero-sequence component and phase offset, combined with the periodic change characteristics, the harmonic source is accurately separated, and the ability to identify the cause of harmonics 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-frequency band priority strategy ensures the compensation optimization of harmonics with different frequencies, reduces system harmonic interference, and enhances the stability of power quality. Description of the Drawings

[0058] Figure 1 is a schematic diagram of the working process of the present invention;

[0059] Figure 2 is a flowchart of the steps for obtaining the harmonic distortion characteristic parameters of the present invention;

[0060] Figure 3 is a flowchart of the steps for obtaining the predicted value of the harmonic distortion trend of the present invention;

[0061] Figure 4 is a flowchart of the steps for obtaining the harmonic source separation result of the present invention;

[0062] Figure 5 is a flowchart of the steps for obtaining the asymmetric compensation parameters of the present invention;

[0063] Figure 6 This is the flowchart of the acquisition steps for the adaptive compensation execution data of the present invention. Detailed implementation manners

[0064] In order to make the objectives, technical solutions and advantages of the present invention clearer and more understandable, the present invention will be 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 used to limit the present invention.

[0065] In the description of the present invention, it should be understood that the orientation or positional relationship indicated by the terms "length", "width", "upper", "lower", "front", "rear", "left", "right", "vertical", "horizontal", "top", "bottom", "inner", "outer", etc. is based on the orientation or positional relationship shown in the accompanying drawings, and is 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 should not be construed as a limitation of the present invention. In addition, in the description of the present invention, "a plurality of" means two or more, unless otherwise specifically defined.

[0066] Embodiment 1: Please refer to Figure 1 , the present invention provides a technical solution: an electrical automation regulation method, including the following steps:

[0067] S1: Obtain three-phase current signals, collect single-cycle high-frequency sampling, extract instantaneous extreme points, record the time coordinates and amplitudes, calculate the time intervals and change rates, screen the distribution of extreme points, and obtain harmonic distortion characteristic parameters;

[0068] S2: Based on the harmonic distortion characteristic parameters, call the instantaneous extreme values, calculate the phase drift rate and the extreme value growth rate, analyze the distortion stability, screen the time periods with prominent fluctuation amplitudes, determine the dominant harmonic components, measure the peak time and the distortion amplitude, construct a trend change curve, and calculate the harmonic distortion trend prediction value;

[0069] S3: Call the harmonic distortion trend prediction value, extract the A, B, and C phase currents in the three-phase current signals, calculate the zero-sequence component, obtain the phase offset and the harmonic distortion amplitude, compare the periodic changes, determine the harmonic source, and obtain the harmonic source separation result;

[0070] S4: Based on the harmonic source separation result, calculate the compensation currents of the A, B, and C phases, compare the multi-phase effects, calculate the injection time points of the harmonic compensation currents, set the compensation current parameters before the harmonic peak arrives, and obtain the asymmetric compensation parameters;

[0071] S5: Invoke the asymmetric compensation parameters, collect the data before compensation, calculate the multi-band priorities, compare the amplitude changes, adjust the current amplitude and phase, and record the adaptive compensation execution data.

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

[0073] Please refer to Figure 2 , and the steps for obtaining the harmonic distortion characteristic parameters are specifically as follows:

[0074] S101: Obtain the three-phase current signals, perform single-cycle high-frequency sampling, record the current amplitudes corresponding to multiple time points, extract the instantaneous extreme points at multiple moments, store the time coordinates and amplitudes of the extreme points, and generate an instantaneous extreme value dataset;

[0075] Obtain the three-phase current signals, and use a high-frequency sampling device to perform single-cycle acquisition of the current signals. The specific operations include performing at least thousands of equally spaced samplings within a complete power frequency cycle (20 ms for a 50 Hz system and 16.67 ms for a 60 Hz system), recording the current amplitude at each sampling point. The sampling interval needs to satisfy the Nyquist sampling theorem, that is, at least twice the highest frequency of the harmonic components. For example, for sampling of a 2 kHz harmonic component, a sampling rate of at least 4 kHz is required. For each sampling point, determine whether it is a local extreme point according to the positive and negative change trends of the current signal. The specific method is as follows: Calculate the change trend of the current values at three adjacent time points. If the current value at the current point is greater than the values at the previous and next points, it is determined as a local maximum; if the current value at the current point is less than the values at the previous and next points, it is determined as a local minimum. The time coordinates and amplitude data of the extreme points are stored in the dataset to generate an instantaneous extreme value dataset, as shown in Table 1.

[0076] Table 1 Example of the instantaneous extreme value dataset

[0077]

[0078] As shown in Table 1, the sampling point number represents the extreme points that appear in the time series. Its time coordinates are recorded by the high-frequency sampling device, and the amplitude is the corresponding instantaneous current value. This dataset provides a basis for subsequent analysis.

[0079] 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.

[0080] 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:

[0081] ;

[0082] Calculate and generate the distribution deviation of the rate of change of extreme value points;

[0083] 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;

[0084] 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:

[0085] ;

[0086] in, is the mean of the rate of change of all extreme points, and the calculation formula is:

[0087] ;

[0088] Assume that the aforementioned extreme point change rate sequence is A / ms, then its mean value is:

[0089] ;

[0090] Calculate the standard deviation:

[0091] ;

[0092] ;

[0093] Calculate the numerator part of the value:

[0094] ;

[0095] ;

[0096] Calculate the denominator part of the value:

[0097] ;

[0098] where is the mean of the time intervals. Assuming , then:

[0099] ;

[0100] ;

[0101] Final calculation:

[0102] ;

[0103] This result indicates that the distribution deviation value of the extreme point change rate has been calculated. The larger the value, the more discrete the distribution of the extreme point change rate. This value will be used for subsequent screening of abnormal data points to determine whether there is an abnormal harmonic distortion situation.

[0104] S103: Invoke the distribution deviation of the extreme point change rate, screen abnormal data points exceeding 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 meeting the set distribution characteristics, obtain the extreme point distribution interval parameters, and generate harmonic distortion characteristic parameters.

[0105] If it exceeds the set threshold (such as 6.0), it is determined that the current data is abnormal. Screen the abnormal data points and count the proportion of abnormal data points. For example, if the total number of extreme points in the dataset is 100 and the number of abnormal data points is 7, then the abnormal distribution rate of data points:

[0106] ;

[0107] ​Combine the time interval screening results, extract the extreme points that meet the set distribution characteristics, that is, when the change rate and the time interval both meet the set thresholds, retain the corresponding extreme points, and finally obtain the extreme point distribution interval parameters. For example, if the set extreme point change rate distribution range is [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.

[0108] This result indicates 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 a significant harmonic interference effect in the power grid.

[0109] Please refer to Figure 3 , and the steps for obtaining the predicted value of the harmonic distortion trend are specifically as follows:

[0110] S201: Based on the harmonic distortion characteristic parameters, call the instantaneous extreme values, compare the extreme value change amplitudes in different time periods, calculate the extreme value growth rate, and screen according to the set extreme value growth rate reference value to extract the time periods with prominent growth rates and obtain the extreme value growth rate screening results;

[0111] First, call the instantaneous extreme values, 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 between adjacent time points, set a time window to statistically analyze these changes to obtain the change trend of the characteristic parameters, and 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:

[0112] ;

[0113] Among them, is the extreme value growth rate of the th harmonic component, and are the harmonic extremes of adjacent time points respectively, is the time interval. By calculating the extreme value growth rates within multiple time windows and further screening according to the set extreme value growth rate reference value, the determination of the reference value can be based on the statistical analysis of historical data. For example, statistically analyze the mean and standard deviation of the growth rates within the past 30 time windows, then the set reference value is:

[0114] ;

[0115] Then screen out the growth rates that exceed Time periods, and these time periods are the time periods with prominent extreme growth rates. As shown in Table 2, a set of example data is listed.

[0116] Table 2 Extreme growth rates

[0117] ;

[0118] As shown in Table 2, based on the set reference value, time points 1 and 3 with relatively high extreme growth rates are screened out, and the screening result of the extreme growth rate is obtained.

[0119] S202: Invoke the screening result of the extreme growth rate, calculate the corresponding phase drift rate for the screened time period, judge the fluctuation condition of the phase drift rate, screen the time interval according to the fluctuation amplitude, and obtain the phase drift rate fluctuation interval;

[0120] It is necessary to obtain the phase change condition of the dominant harmonic component at adjacent time points. By calculating the phase drift rate formula:

[0121] ;

[0122] Among them, and respectively represent the phase angles at time and moments, represents the phase drift rate. Then judge the fluctuation condition of the phase drift rate, set the phase drift rate fluctuation reference value, similar to the reference value setting of the extreme growth rate mentioned above, and set it based on the statistical data:

[0123] ;

[0124] Screen out the time intervals with larger fluctuation amplitudes. For example, if the drift rates of the third harmonic within a certain period are 0.1, 0.2, 0.15, 0.4, 0.05 respectively, the average value is 0.18, and the standard deviation is 0.12, then the reference value is 0.42, and the time periods with drift rates higher than this value are screened out, and finally the phase drift rate fluctuation interval is obtained.

[0125] S203: Invoke the phase drift rate fluctuation interval, measure the peak moment and distortion amplitude of the dominant harmonic component within the fluctuation interval, and construct a trend change curve, using the formula:

[0126] ;

[0127] Calculate to obtain the harmonic distortion trend prediction value;

[0128] Among them, represents the harmonic distortion trend prediction value, represents the The distortion amplitude of the dominant harmonic component represents the attenuation factor of the th dominant harmonic component, represents the time variable, represents the drift factor of the th dominant harmonic component, represents the phase offset of the th dominant harmonic component,

[0129] Formula:

[0130] ;

[0131] 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 etc., is the time variable, and the calculation step size is set , is the drift factor, measured by statistical regression method , is the phase offset, and the value is taken according to the drift rate calculated above

[0132] Assume that the distortion amplitudes of the dominant harmonic components measured within a certain time period are respectively 、 , then calculate the predicted values:

[0133] ;

[0134] ;

[0135] ;

[0136] ;

[0137] The result shows that the predicted value of the harmonic distortion trend is 2.72V, which can be used for subsequent harmonic distortion trend evaluation

[0138] Please refer to Figure 4 , and the specific steps for obtaining the harmonic source separation result are as follows:

[0139] S301: Invoke the predicted value of harmonic distortion trend, extract the currents of phases A, B, and C in the three-phase current signal, calculate the root mean square (RMS) values of the poly-phase current signal over multiple cycles, obtain the in-cycle fluctuation characteristics of the currents of phases A, B, and C. At the same time, perform a difference operation based on the RMS value of the poly-phase current and the zero-sequence component, calculate the current offset of the zero-sequence component, and obtain the zero-sequence component offset rate.

[0140] First, it is necessary to obtain the currents of phases A, B, and C in the three-phase current signal, and perform per-cycle sampling on each phase current signal. Assume that the number of sampling points in one cycle is , then within the range of cycles, accumulate all the sampling values of each phase current. Next, calculate the root mean square (RMS) value for each phase current, that is:

[0141] ;

[0142] where represents phases A, B, and C, is the current value at the th sampling point. After obtaining the RMS values of the three-phase currents of A, B, and C through the above calculations, analyze the fluctuation of their values within the cycle range, that is, calculate the change in the RMS value per cycle, which represents the periodic fluctuation amplitude of the current. Assume that the RMS value of cycle is , then the fluctuation amplitude is calculated as follows:

[0143] ;

[0144] For the zero-sequence component, its calculation is based on the sum of the instantaneous values of the three-phase currents:

[0145] ;

[0146] Subsequently, calculate the RMS value of the zero-sequence component within the cycle and perform a difference operation:

[0147] ;

[0148] To further analyze the change trend of the zero-sequence component, it is necessary to calculate its relative offset rate. Assume that the calculation method of the zero-sequence component offset rate is as follows:

[0149] ;

[0150] This offset rate is used to judge the current fluctuation of the zero-sequence component. If exceeds the set threshold (for example, 0.05), it indicates that the fluctuation of the zero-sequence component within this cycle is significant. As shown in Table 1, the RMS values of the three-phase currents and the zero-sequence component offset rate in a certain measurement are listed.

[0151] Table 3 Root mean square values of three-phase currents and offset rates of zero-sequence components

[0152]

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

[0154] S302: Based on the offset rate of the zero-sequence component, calculate the phase offset of the A, B, and C phase currents relative to the reference phase. Select the A and B phase currents for difference calculation, construct the ratio relationship of the zero-sequence component, and use the formula:

[0155] ;

[0156] Calculate the phase change trends of the A and B phase currents, obtain the phase difference parameters of the A and B phase currents, judge the comparison relationship between the phase change and the zero-sequence component, screen the key phase distortion characteristics, and obtain the main phase distortion amount;

[0157] Among them, represents the phase offset, represents the current signals of the A and B phases in the th cycle, represents the current signal of the zero-sequence component, represents the current difference between the A and B phases, represents the maximum current during the measurement, represents the number of cycles selected during the calculation process, represents the arctangent function operation, which is used to calculate the phase angle corresponding to the current signal ratio;

[0158] Based on the offset rate of the zero-sequence component calculated in the previous paragraph , calculate the phase offset of the A, B, and C phase currents relative to the reference phase. Assuming that the A phase is the reference phase, calculate the phase offsets of the B and C phases relative to the A phase. Assume that within the cycle , the reference phase of the A phase current is , and the instantaneous phase of the B phase current is , then the calculation of its phase offset is as follows:

[0159] ;

[0160] Adopt the differential calculation method to calculate the difference ratio of the A and B phase currents, and construct the ratio relationship of the zero-sequence component:

[0161] ;

[0162] Among them, and are the current values of phase A and phase B in the th cycle respectively, is the zero-sequence component current signal, indicating the current difference between phase A and phase B, and the calculation method is as follows:

[0163] ;

[0164] ;

[0165] Suppose cycles, and the measured phase A and phase B currents are as follows:

[0166] ;

[0167] Calculate the difference:

[0168] ;

[0169] Assume the zero-sequence component current , the maximum current , calculate the phase offset:

[0170] ;

[0171] ;

[0172] ;

[0173] This calculation shows that the phase offset of the phase B current relative to the phase A current is approximately 54.7°, indicating a significant phase difference, which can be used to judge the source of harmonic distortion.

[0174] S303: Call the main phase distortion amount, calculate the corresponding harmonic distortion amplitude, and at the same time, based on the change trend of the harmonic distortion amplitude within multiple cycles, select the extreme points of the harmonic amplitude for calculating the fluctuation range, and compare the periodic change characteristics to judge the type and contribution degree of the harmonic source, and obtain the harmonic source separation result.

[0175] Use Fourier transform to calculate the harmonic components, obtain the amplitudes of each harmonic, and draw the change curve of the harmonic distortion amplitude within multiple cycles. Through the peak detection method, select the extreme points to calculate the fluctuation range:

[0176] ;

[0177] Among them, and are the maximum and minimum harmonic amplitudes respectively. Combining the fluctuation characteristics of multiple cycles, calculate the periodic change trend of the harmonic amplitude. Suppose the 5th harmonic amplitudes of a certain measurement are:

[0178] ;

[0179] Then the fluctuation range:

[0180] ;

[0181] If the harmonic source determination threshold is set to 1.0, then the contribution of this harmonic source is relatively large, and its source and influence can be further analyzed.

[0182] Please refer to Figure 5 , and the specific steps for obtaining the asymmetric compensation parameters are as follows:

[0183] 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 each phase harmonic current, and at the same time, based on the peak amplitude data, calculate the phase change trend. 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;

[0184] 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. Among them, the calculation method of the peak amplitude is to select the maximum positive and negative values of the harmonic current in each period, take the absolute value and calculate the mean value. For example, if the maximum value of the harmonic current of phase A in a certain period is 3.2 A and the minimum value is -3.0 A, then the peak amplitude is calculated as , based on the peak amplitude data in all periods, calculate the average peak amplitude of each phase, and at the same time call the three-phase harmonic current phase data, calculate the phase change trend, use the difference to calculate the phase change rate between adjacent moments, and compare the phase change data of each phase to obtain the phase offset value between the three phases. Among them, 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 , and finally obtain the three-phase harmonic phase offset. Assume that the harmonic current phases of phases A, B, and C are respectively , , , then the phase offsets are respectively , , and finally form a data set of the three-phase harmonic phase offset.

[0185] S402: Call the three-phase harmonic phase offset, calculate the time offset value of the multi-phase compensation current, call the three-phase compensation current time offset data, compare the harmonic compensation influence, calculate the time series difference of the compensation current, obtain the optimal injection time point of the compensation current, and at the same time call the harmonic amplitude data, calculate the compensation amplitude, and establish the optimal compensation current parameters;

[0186] 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 period. For example, if the power frequency period is 20 ms, then the phase offset The corresponding time offset calculation is , similarly The corresponding time offset value is . Combining the time offset data of the three-phase compensation current, comparing the harmonic compensation influence, calculate the timing difference of the compensation current. Among them, the calculation of the timing difference is based on the time alignment degree between the phase of the compensation current and the phase of the actual harmonic current. If the phase of the compensation current 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, obtain the optimal injection time point of the compensation current. At the same time, call the harmonic amplitude data to calculate the compensation amplitude. The calculation method is as follows: Let the amplitude of the harmonic current be , and the amplitude of the compensation current be , then , where is the compensation coefficient. Assuming , the compensation coefficient , then the amplitude of the compensation current . Finally, establish the optimal compensation current parameters.

[0187] S403: Based on the optimal compensation current parameters, set the compensation current adjustment parameters before the harmonic peak arrives, and calculate the parameters of the multi-phase compensation current. Use the formula:

[0188] ;

[0189] Obtain the asymmetric compensation parameters;

[0190] Among them, represents the overall parameter of the compensation current, represents the th-order harmonic voltage, represents the phase angle of the th-order harmonic voltage, represents the phase adjustment parameter, represents the total number of harmonic orders, is the imaginary unit.

[0191] Formula:

[0192] ;

[0193] Among them, represents the overall parameter of the compensation current, represents the th-order harmonic voltage, represents the phase angle of the th-order harmonic voltage, represents the phase adjustment parameter, represents the total number of harmonic orders, is the imaginary unit. Based on the known data for calculation, assuming the total number of harmonic orders , the harmonic voltage amplitudes of each order are respectively , , , and their phase angles are respectively , , , and the phase adjustment parameters are respectively , , , substitute into the formula:

[0194] ;

[0195] Further calculation:

[0196] ;

[0197] Convert the exponential form to the rectangular coordinate system:

[0198] ;

[0199] Calculate three numerical values:

[0200] ;

[0201] ;

[0202] ;

[0203] Calculate the modulus value:

[0204] ;

[0205] Finally, obtain the asymmetric compensation parameter .

[0206] Table 4 Compensation current calculation parameters

[0207]

[0208] As shown in Table 4, the voltage values and their phase angle information of different harmonic orders, as well as the phase adjustment parameters, are listed. The finally calculated compensation current parameter is , and its modulus value is . This result indicates that the calculated asymmetric compensation parameter meets the harmonic current compensation requirements and can be used in the actual compensation control link.

[0209] Please refer to Figure 6 , the specific steps for obtaining the adaptive compensation execution data are as follows:

[0210] S501: Call the asymmetric compensation parameters, collect the data before compensation, calculate the signal amplitudes and spectral characteristics of multiple frequency bands, obtain the signal change rates of different frequency bands, and at the same time screen the frequency bands based on the change rates, calculate the weights of the frequency bands, and sort them to obtain the multi-frequency band priority ranking;

[0211] First, set the signal sampling points of multiple frequency bands, scan in the range of 1 MHz to 100 MHz with a step frequency of 1 kHz, and collect the amplitude data of the signals , and store them in the matrix , and record the initial value of the reference signal amplitude at the same time, for example, set the reference amplitude to 1 V. Then, calculate the signal amplitude change rate of each frequency band, and use the method to obtain the signal change amount of each frequency band, where if V, it is determined that there is a significant change in this frequency band, otherwise the data of this frequency band is ignored. Subsequently, screen the frequency bands with significant changes and calculate the weights of each frequency band , and the weight calculation uses the normalization method, that is:

[0212] ;

[0213] 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 change amplitudes have higher priorities, and finally form a multi-frequency band priority ranking list.

[0214] Table 5 Calculation Table of Multi-Frequency Band Signal Change Rates

[0215]

[0216] As shown in Table 5, the signal change rates of different frequency bands are obtained through calculation. Among them, the change amount of the 20 MHz frequency band is the largest, so it ranks highest in the priority ranking.

[0217] S502: Based on the multi-frequency band priority ranking, calculate the signal amplitude change of the selected frequency band, call the initial current amplitude and phase parameters, calculate the amplitude change ratio of the multi-frequency band signals, and compare it with the set reference change threshold, using the formula:

[0218] ;

[0219] Calculate the current adjustment factor, adjust the current amplitude and phase of the multi-frequency band, and obtain the current correction parameter;

[0220] Among them, represents the current adjustment factor, represents the The current signal amplitude of the frequency band, represents the reference signal amplitude, represents the weight of the frequency band, represents the number of selected frequency bands, 、 indicates the operation of capturing the maximum and minimum values among all selected frequency bands;

[0221] Select the first frequency bands and calculate the signal amplitude change ratio. First, call the initial current amplitude A and the phase parameter to calculate the current amplitude adjustment factor for the selected frequency bands. The calculation method is as follows:

[0222] 1. Obtain the amplitude change range of the selected frequency bands, that is, and , for example, for the selected frequency bands {5MHz, 10MHz, 20MHz}, the maximum change is 0.20V and the minimum change is 0.08V;

[0223] 2. Calculate the signal change amplitude difference:

[0224] ;

[0225] 3. Calculate the weighted average weight:

[0226] ;

[0227] 4. Calculate the current adjustment factor:

[0228] ;

[0229] The calculation result indicates that the current amplitude and phase need to be adjusted to compensate for the signal change.

[0230] S503: Call the current correction parameter to adjust the actually applied current amplitude and phase, and record the adjusted signal response in real time. Calculate the difference in signal offset before and after adjustment to obtain the adaptive compensation execution data.

[0231] Call the calculated current correction parameter A to adjust the actually applied current, that is, the new current amplitude:

[0232] ;

[0233] At the same time, the phase adjustment amount is determined according to the correction parameter, and the new phase is set to , where depends on the phase compensation model of the specific system and is set to , namely . Record the adjusted signal response in real time, calculate the signal offset before and after adjustment, and finally obtain the adaptive compensation execution data.

[0234] An electrical automation regulation system, which is used to execute the above electrical automation regulation method. The system includes:

[0235] The harmonic distortion feature extraction module obtains the three-phase current signal, extracts the single-cycle high-frequency sampling data, screens the instantaneous extreme points in the sampling points, records the time coordinates and amplitudes of the instantaneous extreme points, calculates the time intervals of consecutive instantaneous extreme points, measures the extreme value change rate, screens out the extreme points with prominent fluctuation amplitudes of the change rate, analyzes the distribution density of the extreme points, extracts the distortion frequency characteristics, calculates the relative amplitudes of the fundamental wave and each harmonic, measures the total harmonic distortion rate, establishes the distortion trend curve, and obtains the harmonic distortion characteristic parameters;

[0236] The phase offset calculation module, based on the harmonic distortion characteristic parameters, calls the instantaneous extreme point data, calculates the phase drift rate of each harmonic, screens out the harmonic components with prominent changes in the phase drift rate, measures the extreme value growth rate, analyzes the stability of the harmonic distortion, extracts the time periods with prominent distortion amplitudes, calculates the peak time points and the corresponding distortion amplitudes, establishes the phase offset trend curve, and obtains the harmonic distortion trend prediction value;

[0237] 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 change rate of the distortion amplitude within a period, screens out the harmonic components with prominent changes in the distortion amplitude within a period, compares the distortion distribution of each phase harmonic component, analyzes the harmonic source, and obtains the harmonic source separation result;

[0238] The compensation current calculation module, based on the harmonic source separation result, calculates the A, B, and C phase harmonic compensation currents, screens out the harmonic components with prominent multi-phase influence, calculates the injection time points of the compensation current, measures the time interval before the harmonic peak arrives, calculates the compensation current parameters, measures the differences of the compensation currents in different phases, and obtains the asymmetric compensation parameters;

[0239] 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 multi-frequency band priority, compares the changes in the current amplitudes before and after compensation, adjusts the amplitude and phase of the compensation current, records the compensation execution data, and obtains the adaptive compensation execution data.

[0240] The above are only the preferred embodiments of the present invention, and are not intended to limit the present invention in other forms. Any person skilled in the relevant art may use the technical content disclosed above to make changes or modifications into equivalent embodiments with equivalent changes and apply them to other fields. However, as long as it does not depart from the technical solution content of the present invention, any simple modification, equivalent change and modification made to the above embodiments based on the technical essence of the present invention still fall 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: 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.

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 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.

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, 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.

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 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 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, 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.

6. The electrical automation adjustment method according to claim 5, characterized in that: The steps of 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 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.

7. The electrical automation adjustment method according to claim 6, characterized in that: 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.

8. The electrical automation adjustment method according to claim 7, characterized in that: 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.

9. An electrical automation regulation system, characterized in that: According to the electrical automation regulation method according to any one of claims 1 to 8, 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.

Images