A Method and System for Suppressing DC Ice-Melting Current Oscillation

Abstract

The present invention relates to the field of power electronics technology, and particularly to a method and system for suppressing DC ice melting current oscillation, including: obtaining an equivalent inductance value and an equivalent capacitance value according to the operation environment data of the DC ice melting system; when it is evaluated based on the equivalent inductance value and the equivalent capacitance value that the current DC ice melting system state does not meet the condition for suppressing current oscillation, generating an optimized valve group control strategy based on the equivalent inductance value and the operation environment data; during the process of controlling the DC ice melting system to perform ice melting operation through the optimized valve group control strategy, extracting current oscillation characteristic data according to the current waveform data, and determining the ice melting efficiency according to the current oscillation characteristic data; when the ice melting efficiency does not meet the standard, optimizing the valve group control optimization strategy by using the current oscillation characteristic data and a multi-objective optimization algorithm, and outputting an optimal ice melting control strategy. The present invention dynamically optimizes the valve group control strategy, taking into account both the suppression of current oscillation and the ice melting efficiency, and ensuring the stable operation of the DC ice melting system.

Description

Technical Field

[0001] The present invention relates to the technical field of power electronics, and particularly to a method and system for suppressing DC ice melting current oscillation. Background Art

[0002] In a mobile DC ice melting system, the turn-on process of the valve group current is prone to causing current oscillation. This phenomenon is mainly caused by the interaction between the internal inductance and capacitance of the system. At the moment when the valve group is turned on, due to the interaction between the internal inductance and capacitance, there is often a sharp rise in current and high-frequency oscillation. This non-ideal current change not only exacerbates the internal stress level of the system, but also may lead to equipment damage or unsatisfactory ice melting effect. To solve this problem, the prior art usually adopts the method of increasing the system inductance to limit the current rise rate, thereby suppressing the generation of oscillation. This will cause the system response speed to slow down. In ice melting operations that require rapid response, the slowdown of the response speed affects the efficiency and real-time performance of ice melting, making the prior art face many challenges in dealing with current oscillation problems.

[0003] Mobile DC ice melting systems usually need to operate in extreme environmental conditions such as extremely cold climates or high-altitude areas. These special environments have an inestimable impact on the electrical parameters of the equipment, making the characteristics of current oscillation more difficult to predict and control. Environmental parameter fluctuations such as temperature, altitude, and humidity changes not only directly change the equivalent inductance and capacitance values of the system, but also indirectly affect the ice melting efficiency, making the suppression of current oscillation more difficult. In addition, to effectively suppress current oscillation, it is necessary to precisely control the turn-on strategy of the valve group and the current path, but this often conflicts with the goal of pursuing the best ice melting efficiency. When real-time monitoring the system response and ice melting effect, the selection of sampling frequency and accuracy also becomes a thorny problem, which is directly related to the accuracy of data analysis and the effectiveness of control strategies. At the same time, in the process of optimizing the layout of heating elements to improve the ice melting efficiency, the optimization of the layout will introduce new oscillation factors. For example, the distribution and connection method of heating elements will change the equivalent inductance and capacitance of the system, thereby affecting the current oscillation characteristics. Therefore, how to avoid introducing new oscillation factors while optimizing the layout of heating elements has become a difficult problem in the prior art.

[0004] In summary, the prior art faces multiple challenges in suppressing current oscillation and improving ice melting efficiency. How to seek the best balance between suppressing oscillation and improving ice melting efficiency has become an urgent problem to be solved in mobile DC ice melting technology. Summary of the Invention

[0005] To solve the above technical problems, the present invention provides a method and system for suppressing DC ice melting current oscillation.

[0006] In a first aspect, the present invention provides a method for suppressing DC ice melting current oscillation, the method comprising the following steps:

[0007] Obtain in real time the operating environment data of the DC ice melting system, and obtain the equivalent inductance value and equivalent capacitance value of the DC ice melting system under the current environmental conditions according to the operating environment data;

[0008] Calculate the theoretical current rising rate according to the equivalent inductance value and the equivalent capacitance value, and evaluate whether the current state of the current DC ice melting system meets the condition for suppressing current oscillation according to the theoretical current rising rate;

[0009] When the current state of the current DC ice melting system does not meet the condition for suppressing current oscillation, optimize the opening time, opening duration and current path selection logic of the ice melting valve group based on the equivalent inductance value and the operating environment data, and generate an optimized valve group control strategy;

[0010] During the process of controlling the DC ice melting system to perform ice melting operation through the optimized valve group control strategy, obtain the current waveform data of the DC ice melting system in real time;

[0011] Extract the current oscillation characteristic data from the current waveform data by using the Fourier transform algorithm, and determine the ice melting efficiency according to the current oscillation characteristic data;

[0012] When it is detected that the ice melting efficiency does not meet the standard, optimize the valve group control optimization strategy by using the current oscillation characteristic data and the multi-objective optimization algorithm, and output the optimal ice melting control strategy.

[0013] In a further embodiment, the step of obtaining the equivalent inductance value and equivalent capacitance value of the DC ice melting system under the current environmental conditions according to the operating environment data includes:

[0014] Perform standardization processing on the operating environment data to obtain the preprocessed operating environment data, and construct an equivalent circuit model based on the actual circuit topology of the DC ice melting system through the electromagnetic field theory;

[0015] Input the preprocessed operating environment data into the equivalent circuit model, and use the least squares method to solve the equivalent inductance value and equivalent capacitance value under the current environmental conditions.

[0016] In a further embodiment, the step of calculating the theoretical current rising rate according to the equivalent inductance value and the equivalent capacitance value includes:

[0017] Based on the equivalent inductance value and the equivalent capacitance value, use Kirchhoff's voltage law to construct a differential equation of the current changing with time in the equivalent circuit model;

[0018] Set the initial current condition at the moment when the valve group is turned on to zero, solve the differential equation using the fourth-order Runge-Kutta method, and extract the current time-domain waveform after the valve group is turned on.

[0019] Calculate the instantaneous change rate when the current starts to rise after the valve group is turned on according to the current time-domain waveform, and obtain the theoretical rise rate.

[0020] In a further embodiment, the step of evaluating whether the current DC ice melting system state meets the condition for suppressing current oscillation according to the theoretical current rise rate includes:

[0021] Compare the theoretical current rise rate with a preset current rise rate threshold. If the theoretical current rise rate is less than or equal to the current rise rate threshold, it is determined that the current DC ice melting system state meets the condition for suppressing current oscillation, and the DC ice melting system is controlled to continue the ice melting operation;

[0022] If the theoretical current rise rate is greater than the current rise rate threshold, it is determined that the current DC ice melting system state does not meet the condition for suppressing current oscillation.

[0023] In a further embodiment, the step of optimizing the opening time, opening duration, and current path selection logic of the ice melting valve group based on the equivalent inductance value and operation environment data to generate a valve group control optimization strategy includes:

[0024] Use a sliding window to detect outliers in the equivalent inductance value and the operation environment data, and smooth the outliers to obtain a standardized equivalent inductance value and standardized operation environment data;

[0025] Based on the standardized equivalent inductance value and standardized operation environment data, use the association rule mining algorithm to analyze the correlation between the equivalent inductance value and environmental factors, and obtain the equivalent inductance distribution characteristics under the current operation environment conditions; the equivalent inductance distribution characteristics include the equivalent inductance change rate;

[0026] Analyze the equivalent inductance value distribution characteristics under the current operation environment conditions, and identify the oscillation influence index of the standardized equivalent inductance value on the amplitude of the ice melting current oscillation;

[0027] According to the equivalent inductance distribution characteristics and the oscillation influence index, use a two-dimensional fuzzy controller to determine the opening time correction amount and the opening duration correction coefficient;

[0028] Based on the opening time correction amount and the opening duration correction coefficient, use the centroid method to defuzzify to obtain the opening time of the ice melting valve group and the opening duration of the ice melting valve group;

[0029] Taking the impedance value between nodes as the weight, generate a power grid structure diagram according to the actual circuit topology of the DC ice melting system;

[0030] Under the opening time and the opening duration of the ice melting valve group, use the graph theory algorithm to search for candidate current paths from the power source to the ice melting point in the power grid structure diagram;

[0031] Calculate the impedance matching degree of each candidate current path, and select the candidate current path with the smallest impedance matching degree as the current optimal current path;

[0032] Integrate the opening time of the ice melting valve group, the opening duration of the ice melting valve group and the current optimal current path to form a valve group control optimization strategy.

[0033] In a further embodiment, the steps of using the Fourier transform algorithm to extract current oscillation characteristic data from the current waveform data and determining the ice melting efficiency according to the current oscillation characteristic data include:

[0034] Use the Fourier transform algorithm to perform frequency domain transformation on the current waveform data to obtain current waveform spectrum data;

[0035] Identify the main frequency components and oscillation intensity of the current oscillation from the current waveform spectrum data to obtain current oscillation characteristic data;

[0036] According to the relationship between the historical current oscillation characteristic data and the historical ice melting efficiency, use the random forest regression algorithm to construct an ice melting efficiency evaluation model;

[0037] Input the current oscillation characteristic data into the ice melting efficiency evaluation model to obtain the ice melting efficiency under the current operating environment conditions.

[0038] In a further embodiment, the current oscillation characteristic data includes oscillation frequency data and oscillation amplitude data; the oscillation intensity is the sum of the squares of the amplitudes of all frequency components;

[0039] Among them, the oscillation frequency data includes the main frequency components of the current oscillation; the oscillation amplitude data includes the amplitude of the fundamental wave component and the amplitudes of each harmonic component.

[0040] In a further embodiment, the steps of optimizing the valve group control optimization strategy by using the current oscillation characteristic data and the multi-objective optimization algorithm and outputting the optimal ice melting control strategy when it is detected that the ice melting efficiency does not meet the standard include:

[0041] When it is detected that the ice melting efficiency does not meet the standard, analyze the correlation relationship between different valve group control parameters and the current oscillation and the ice melting efficiency according to the current oscillation characteristic data, and determine the optimization variable range of each valve group control parameter;

[0042] Randomly generate individuals containing several groups of valve group control parameter combinations within the optimization variable range to form an initial population;

[0043] Simulate the ice melting process for each individual in the population, and evaluate the total harmonic distortion rate index and the ice melting efficiency index of each individual during the ice melting process through the valve group control parameters of each individual;

[0044] Calculate the fitness value of each individual according to the total harmonic distortion rate index and the ice melting efficiency index, and select the individuals that meet the fitness threshold from the current population as the parent individuals by using the roulette wheel selection method according to the fitness value;

[0045] Perform crossover operations on the continuous variables in the parent individuals by using the arithmetic crossover method, and perform crossover operations on the discrete variables in the parent individuals by using the single-point crossover method to generate new individuals; the continuous variables include the opening time and the opening duration; the discrete variables include the current path selection;

[0046] Perform mutation operations on the new individuals generated by crossover, and repeat the above selection, crossover, and mutation processes until the fitness converges to a preset convergence threshold;

[0047] After the iteration process ends, select the individual with the highest fitness value from the final population to obtain the optimal ice melting control strategy.

[0048] In a further implementation, the process of judging whether the ice melting efficiency meets the standard is specifically as follows:

[0049] Compare the ice melting efficiency with a preset ice melting efficiency threshold. If the ice melting efficiency is lower than the ice melting efficiency threshold, it is determined that the ice melting efficiency does not meet the standard; if the ice melting efficiency is higher than or equal to the ice melting efficiency threshold, continue to maintain the current valve group control optimization strategy.

[0050] In a second aspect, the present invention provides a DC ice melting current oscillation suppression system, and the system includes:

[0051] A data acquisition module, configured to acquire the operation environment data of the DC ice melting system in real time, and obtain the equivalent inductance value and the equivalent capacitance value of the DC ice melting system under the current environmental conditions according to the operation environment data;

[0052] An oscillation evaluation module, configured to calculate the theoretical current rise rate according to the equivalent inductance value and the equivalent capacitance value, and evaluate whether the current state of the current DC ice melting system meets the condition for suppressing current oscillation according to the theoretical current rise rate;

[0053] A strategy generation module, configured to optimize the opening time, the opening duration, and the current path selection logic of the ice melting valve group based on the equivalent inductance value and the operation environment data when the current state of the DC ice melting system does not meet the condition for suppressing current oscillation, and generate a valve group control optimization strategy;

[0054] The ice melting operation module is used to obtain the current waveform data of the DC ice melting system in real time during the process of controlling the DC ice melting system to perform ice melting operations through the valve group control optimization strategy;

[0055] The oscillation analysis module is used to extract current oscillation characteristic data from the current waveform data by using the Fourier transform algorithm, and determine the ice melting efficiency according to the current oscillation characteristic data;

[0056] The strategy optimization module is used to optimize the valve group control optimization strategy by using the current oscillation characteristic data and the multi-objective optimization algorithm when it is detected that the ice melting efficiency does not meet the standard, and output the optimal ice melting control strategy.

[0057] The present invention provides a method and system for suppressing DC ice melting current oscillation. The method obtains the operating environment data of the DC ice melting system in real time, and obtains the equivalent inductance value and equivalent capacitance value of the DC ice melting system under the current environmental conditions according to the operating environment data; calculates the theoretical current rise rate according to the equivalent inductance value and equivalent capacitance value, and evaluates whether the current state of the current DC ice melting system meets the condition for suppressing current oscillation according to the theoretical current rise rate; when the current state of the current DC ice melting system does not meet the condition for suppressing current oscillation, optimizes the opening time, opening duration and current path selection logic of the ice melting valve group based on the equivalent inductance value and operating environment data to generate a valve group control optimization strategy; obtains the current waveform data of the DC ice melting system in real time during the process of controlling the DC ice melting system to perform ice melting operations through the valve group control optimization strategy; extracts current oscillation characteristic data from the current waveform data by using the Fourier transform algorithm, and determines the ice melting efficiency according to the current oscillation characteristic data; when it is detected that the ice melting efficiency does not meet the standard, optimizes the valve group control optimization strategy by using the current oscillation characteristic data and the multi-objective optimization algorithm, and outputs the optimal ice melting control strategy. Compared with the prior art, this method dynamically optimizes the valve group control strategy by real-time monitoring of the operating environment data and current waveform data, ensuring that the DC ice melting system achieves the optimal ice melting efficiency while suppressing current oscillation, and guaranteeing the stable operation of the DC ice melting system. Description of the Drawings

[0058] Figure 1 is a schematic flowchart of the method for suppressing DC ice melting current oscillation provided by an embodiment of the present invention;

[0059] Figure 2 is a block diagram of the system for suppressing DC ice melting current oscillation provided by an embodiment of the present invention. Detailed Embodiments

[0060] The following specifically illustrates the embodiments of the present invention in conjunction with the accompanying drawings. The given examples are only for illustrative purposes and should not be construed as limiting the present invention. The included drawings are for reference and explanation only and do not constitute a limitation on the scope of patent protection of the present invention, because many changes can be made to the present invention without departing from the spirit and scope of the present invention.

[0061] Referring to Figure 1 , an embodiment of the present invention provides a method for suppressing direct current ice melting current oscillation. As Figure 1 shown, the method includes the following steps:

[0062] S1. Obtain the operating environment data of the direct current ice melting system in real time, and obtain the equivalent inductance value and equivalent capacitance value of the direct current ice melting system under the current environmental conditions according to the operating environment data.

[0063] In some embodiments, the step of obtaining the equivalent inductance value and equivalent capacitance value of the direct current ice melting system under the current environmental conditions according to the operating environment data includes:

[0064] Perform standardization processing on the operating environment data to obtain preprocessed operating environment data, and based on the actual circuit topology of the direct current ice melting system, construct an equivalent circuit model through electromagnetic field theory;

[0065] Input the preprocessed operating environment data into the equivalent circuit model, and use the least squares method to solve the equivalent inductance value and equivalent capacitance value under the current environmental conditions.

[0066] Specifically, in this embodiment, various types of sensors are arranged around the wires and key components of the direct current ice melting system, such as temperature sensors, humidity sensors, wind speed sensors, and air pressure sensors, as well as electrical parameter sensors such as voltage sensors and current sensors. Then, a sensor network is used to collect the operating environment data of the direct current ice melting system in real time. The operating environment data includes information such as temperature, altitude, and humidity. In this embodiment, the collected operating environment data is cleaned to remove outliers and noise, and the cleaned data is standardized to convert the data into a unified scale for subsequent analysis and processing. Then, in this embodiment, the circuit topology of the system is analyzed according to the actual wiring diagram and circuit principle of the direct current ice melting system. Based on the actual circuit topology of the direct current ice melting system, an equivalent circuit model is constructed using electromagnetic field theory. The model includes key circuit elements such as inductance and capacitance to reflect the electrical characteristics of the system under different environmental conditions. Input the preprocessed operating environment data into the equivalent circuit model, and by minimizing the gap between the actual measured value and the model predicted value, use methods such as the least squares method to solve the equivalent inductance value and equivalent capacitance value under the current environmental conditions.

[0067] S2. Calculate the theoretical current rise rate based on the equivalent inductance value and the equivalent capacitance value, and evaluate whether the current oscillation suppression condition is satisfied according to the current state of the current DC de-icing system based on the theoretical current rise rate.

[0068] In some embodiments, the step of calculating the theoretical current rise rate based on the equivalent inductance value and the equivalent capacitance value includes:

[0069] Based on the equivalent inductance value and the equivalent capacitance value, use Kirchhoff's voltage law to construct a differential equation for the current change with time in the equivalent circuit model;

[0070] Set the initial current condition at the moment of valve group opening to zero, use the fourth-order Runge-Kutta method to solve the differential equation, and extract the current time-domain waveform after the valve group is opened;

[0071] Calculate the instantaneous change rate when the current starts to rise after the valve group is opened according to the current time-domain waveform to obtain the theoretical rise rate.

[0072] Specifically, in the equivalent circuit model of the DC de-icing system, in this embodiment, the voltage equation is established by Kirchhoff's voltage law (KVL). For a circuit including inductance, capacitance, and resistance, the KVL equation is expressed as:

[0073]

[0074] where is the externally applied voltage; L is the equivalent inductance value; is the current; R is the total resistance in the equivalent circuit model; Q is the charge on the capacitor; C is the equivalent capacitance value; t is the time.

[0075] Since the circuit is in an open state before the valve group is turned on, in this embodiment, it is assumed that the current is zero at the moment when the valve group is turned on. Then, the fourth-order Runge-Kutta method (RK4) is used as the numerical method for solving the differential equation. According to the fourth-order Runge-Kutta method, the time domain is divided into multiple small intervals. For each small interval, the fourth-order Runge-Kutta method is used to calculate the approximate value of the current at the end of the small interval. This process is repeated until the entire required time range is covered, so as to obtain the complete waveform of the current changing with time. The calculated current values are arranged in chronological order to form the current time-domain waveform after the valve group is turned on. In this embodiment, the point at which the current starts to increase significantly after the valve group is turned on is determined by observing the position where the current time-domain waveform starts to increase from a value close to zero. Using a numerical differentiation method (such as the two-point difference method), the instantaneous change rate of the current at this point, that is, the theoretical current rise rate, is calculated. At the same time, in this embodiment, a current rise rate threshold is set according to the system design and safety standards. This current rise rate threshold should ensure that based on the maximum current change rate that the system can withstand without causing unwanted oscillations or damage. In this embodiment, the calculated theoretical current rise rate is compared with the preset current rise rate threshold. If the theoretical rise rate is less than or equal to the current rise rate threshold, it is considered that the current state of the DC ice melting system meets the conditions for suppressing current oscillations, and the system has good stability, and the DC ice melting system can be controlled to continue the ice melting operation; otherwise, if the theoretical rise rate is greater than the current rise rate threshold, it is considered that the current state of the DC ice melting system does not meet the conditions for suppressing current oscillations, there is a risk of current oscillations, and system parameters need to be adjusted or other measures need to be taken to reduce the risk of current oscillations in order to reduce the current rise rate and ensure the stable operation of the system.

[0076] S3. When the current state of the DC ice melting system does not meet the conditions for suppressing current oscillations, optimize the opening time, opening duration, and current path selection logic of the ice melting valve group based on the equivalent inductance value and operating environment data, and generate a valve group control optimization strategy.

[0077] In some embodiments, the steps of optimizing the opening time, opening duration, and current path selection logic of the ice melting valve group based on the equivalent inductance value and operating environment data, and generating a valve group control optimization strategy include:

[0078] Use a sliding window to detect the outliers in the equivalent inductance value and the operating environment data, and smooth the outliers to obtain the standardized equivalent inductance value and the standardized operating environment data;

[0079] Based on the standardized equivalent inductance value and the standardized operating environment data, use the association rule mining algorithm to analyze the correlation between the equivalent inductance value and environmental factors, and obtain the equivalent inductance distribution characteristics under the current operating environment conditions; the equivalent inductance distribution characteristics include the equivalent inductance change rate;

[0080] Analyze the distribution characteristics of the equivalent inductance value under the current operating environment conditions, and identify the oscillation influence index of the standardized equivalent inductance value on the amplitude of the ice melting current oscillation;

[0081] According to the equivalent inductance distribution characteristics and the oscillation influence index, use a two-dimensional fuzzy controller to determine the correction amount of the turn-on timing and the correction coefficient of the turn-on duration;

[0082] Based on the correction amount of the turn-on timing and the correction coefficient of the turn-on duration, use the centroid method to defuzzify to obtain the turn-on timing of the ice melting valve group and the turn-on duration of the ice melting valve group;

[0083] Taking the impedance value between nodes as the weight, generate a power grid structure diagram according to the actual circuit topology of the DC ice melting system;

[0084] At the turn-on timing of the ice melting valve group and the turn-on duration of the ice melting valve group, use graph theory algorithms to search for candidate current paths from the power source to the ice melting point in the power grid structure diagram;

[0085] Calculate the impedance matching degree of each candidate current path, and select the candidate current path with the smallest impedance matching degree as the current optimal current path;

[0086] Integrate the turn-on timing of the ice melting valve group, the turn-on duration of the ice melting valve group and the current optimal current path to form an optimized valve group control strategy.

[0087] Specifically, in this embodiment, a sliding window with a fixed length is set, and the size of the sliding window is adjusted according to the frequency and stability of data changes. For example, a window containing the latest N data points can be selected. Within the sliding window, this embodiment uses statistical methods such as Z-score or IQR to calculate the statistical characteristics of data points. The statistical characteristics include the mean and standard deviation. For each data point within the window, determine whether it exceeds the standard deviation of several times the mean. The data points outside the range are outlier points. This embodiment uses methods such as linear interpolation or average value replacement to smooth the detected outlier points to eliminate the influence of abnormal data on subsequent analysis, so as to obtain the standardized equivalent inductance value and the standardized operating environment data, and discretize the standardized equivalent inductance value and the standardized operating environment data, convert continuous numerical values into discrete intervals or categories, and use association rule mining algorithms such as FP-Growth to analyze the correlation between the equivalent inductance value and environmental factors to obtain the equivalent inductance distribution characteristics under the current operating environment conditions. The equivalent inductance distribution characteristics include the equivalent inductance change rate, etc.

[0088] Meanwhile, in this embodiment, the oscillation influence index is defined as the correlation index between the standardized equivalent inductance value and the oscillation amplitude of the ice melting current. The oscillation influence index is used to measure the influence degree of the change in the equivalent inductance value on the current oscillation. In this embodiment, regression analysis is adopted to establish a non - linear regression model between the standardized equivalent inductance value and the current oscillation amplitude, and according to the established non - linear regression model, the oscillation influence index of the standardized equivalent inductance value on the ice melting current oscillation amplitude under the current operating environment conditions is calculated. Specifically, in this embodiment, the regression coefficient of the non - linear regression model is used to determine the amplitude of the corresponding change in the current oscillation amplitude when the equivalent inductance value changes by one unit, so as to obtain the oscillation influence index. Then, in this embodiment, a two - dimensional fuzzy controller with the equivalent inductance distribution characteristics (such as the equivalent inductance change rate) and the oscillation influence index as input variables, and the opening time correction amount and the opening duration correction coefficient as output variables is established. The fuzzy sets and membership functions of the input and output variables are defined. Specifically, the equivalent inductance change rate is divided into three fuzzy subsets: fast, medium, and slow, the oscillation influence index is divided into three fuzzy subsets: high, medium, and low, and the opening time correction amount and the opening duration correction coefficient are respectively divided into fuzzy subsets such as positive, zero, and negative, and the corresponding membership functions are determined. Fuzzy control rules are formulated based on expert experience and system analysis. If the equivalent inductance change rate is fast and the oscillation influence index is high, the opening time correction amount should be advanced and the opening duration correction coefficient should be increased; if the equivalent inductance change rate is slow and the oscillation influence index is low, the opening time correction amount can be delayed and the opening duration correction coefficient can be decreased. According to the input equivalent inductance distribution characteristics and the oscillation influence index, the corresponding fuzzy output is determined by using the Mamdani fuzzy inference method, and the center - of - gravity method is used to defuzzify the output value of the two - dimensional fuzzy controller. In this embodiment, the initially set opening time is added to the opening time correction amount to obtain the optimized ice melting valve group opening time. For example, if the initial opening time is , and the opening time correction amount is , then the optimized opening time is ; meanwhile, the initially set opening duration is multiplied by the opening duration correction coefficient to obtain the optimized ice melting valve group opening duration. For example, if the initial opening duration is , and the opening duration correction coefficient is k, then the optimized opening duration is .

[0089] For the current path selection logic, in this embodiment, according to the wiring diagram and circuit principle of the DC ice melting system, the connection relationship and impedance parameters between each node are determined, and a power grid structure diagram is constructed with nodes as vertices and electrical connections between nodes as edges. In the power grid structure diagram, the weight value of each edge is set as the impedance value between the corresponding nodes. If there is no direct connection between nodes, the weight is set to infinity or zero. Then, at the determined opening time and opening duration of the ice melting valve group, this embodiment uses graph theory algorithms such as the Dijkstra algorithm to search for multiple short paths from the power supply node to the ice melting point node in ascending order of impedance value according to the weights (impedance values) of each edge in the power grid structure diagram, forming candidate current paths. At the same time, for each candidate current path, this embodiment calculates the impedance matching degree of each candidate current path according to the relative error between the total impedance of the candidate current path (the sum of the impedances of each segment in the path) and the ideal impedance of the system. Among all candidate current paths, this embodiment selects the candidate current path with the smallest impedance matching degree as the current optimal current path. This path can make the ice melting current have less loss during transmission and better meet the impedance matching requirements of the system, which is beneficial to improving the ice melting efficiency and suppressing current oscillation. This embodiment integrates the opening time of the ice melting valve group, the opening duration of the ice melting valve group, and the current optimal current path to form a valve group control optimization strategy, and outputs the valve group control optimization strategy to the control system of the DC ice melting system. The control system adjusts the opening operation of the ice melting valve group and the current path selection according to this strategy to implement the optimized ice melting control process. Through the above steps, this embodiment can dynamically adjust the ice melting control parameters and paths according to different operating environment conditions, effectively suppress current oscillation, improve the operating stability and ice melting efficiency of the DC ice melting system, and ensure the safe and reliable operation of the power system.

[0090] S4. During the process of controlling the DC ice melting system to perform ice melting operations through the valve group control optimization strategy, the current waveform data of the DC ice melting system is acquired in real time.

[0091] S5. Use the Fourier transform algorithm to extract the current oscillation characteristic data from the current waveform data, and determine the ice melting efficiency according to the current oscillation characteristic data.

[0092] In some embodiments, the steps of using the Fourier transform algorithm to extract the current oscillation characteristic data from the current waveform data and determining the ice melting efficiency according to the current oscillation characteristic data include:

[0093] Perform frequency domain transformation on the current waveform data using the Fourier transform algorithm to obtain the current waveform spectrum data;

[0094] Identify the main frequency components and oscillation intensity of the current oscillation from the current waveform spectrum data to obtain the current oscillation characteristic data;

[0095] According to the relationship between the historical current oscillation characteristic data and the historical ice melting efficiency, a random forest regression algorithm is used to construct an ice melting efficiency evaluation model;

[0096] Input the current oscillation characteristic data into the ice melting efficiency evaluation model to obtain the ice melting efficiency under the current operating environment conditions.

[0097] Specifically, in this embodiment, the collected current waveform data is preprocessed. The preprocessing includes operations such as denoising and removing the DC component to improve the accuracy of the Fourier transform. Then, the preprocessed current waveform data is transformed from the time domain to the frequency domain using the Fourier transform algorithm to obtain the current waveform spectrum data. In this embodiment, by analyzing the current waveform spectrum data, frequency components with larger amplitudes are searched for, and the main frequency components of the current oscillation are identified, that is, the frequency points with amplitudes significantly greater than the background noise. These frequency components reflect the main oscillation modes in the current. Then, in this embodiment, the sum of the squares of the amplitudes of all frequency components is calculated as a measure of the oscillation intensity. The identified main frequency components and their corresponding amplitude, phase and other information are organized into current oscillation characteristic data. The current oscillation characteristic data includes oscillation frequency data and oscillation amplitude data; the oscillation intensity is the sum of the squares of the amplitudes of all frequency components; the oscillation frequency data includes the specific frequency values of the main frequency components of the current oscillation; the oscillation amplitude data includes the fundamental wave component amplitude and the amplitudes of each harmonic component.

[0098] Meanwhile, this embodiment collects the historical operation data of the DC ice melting system over a period of time. The historical operation data includes historical current oscillation characteristic data and corresponding historical ice melting efficiency data. The historical current oscillation characteristic data includes oscillation frequency data, oscillation amplitude data, etc. under different operation environmental conditions and at different ice melting stages. The historical ice melting efficiency data can be obtained by methods such as actually measuring the change in the ice layer thickness on the wire surface and the change in the ice layer weight before and after ice melting. This embodiment selects characteristic variables with strong correlation with the historical ice melting efficiency from the historical current oscillation characteristic data as the model input, such as the amplitude of the main frequency component, the total oscillation intensity, etc., and uses the corresponding historical ice melting efficiency data as the output label of the model to construct a training sample set, and trains the ice melting efficiency evaluation model using the random forest regression algorithm. This embodiment uses the training sample set to train the random forest regression model. During the training process, the random forest algorithm extracts multiple sub-sample sets from the training sample set through the Bootstrap sampling method, and each sub-sample set is used to train a decision tree. During the construction of each decision tree, the node splitting is performed by randomly selecting a feature subset to increase the diversity of the model. The current current oscillation characteristic data obtained through the above steps is input into the trained ice melting efficiency evaluation model. The model predicts the input data through multiple internal decision trees. Each decision tree outputs an ice melting efficiency prediction value, and the random forest regression model then averages or weights and averages these prediction values to obtain the final ice melting efficiency prediction result.

[0099] S6. When it is detected that the ice melting efficiency does not meet the standard, use the current oscillation characteristic data and the multi-objective optimization algorithm to optimize the valve group control optimization strategy, and output the optimal ice melting control strategy.

[0100] In some embodiments, the step of using the current oscillation characteristic data and the multi-objective optimization algorithm to optimize the valve group control optimization strategy and output the optimal ice melting control strategy when it is detected that the ice melting efficiency does not meet the standard includes:

[0101] When it is detected that the ice melting efficiency does not meet the standard, analyze the correlation between different valve group control parameters and current oscillation and ice melting efficiency according to the current oscillation characteristic data, and determine the optimization variable range of each valve group control parameter;

[0102] Randomly generate individuals containing several groups of valve group control parameter combinations within the optimization variable range to form an initial population;

[0103] Perform simulation calculations on the ice melting process for each individual in the population, and evaluate the total harmonic distortion rate index and ice melting efficiency index of each individual during the ice melting process through the valve group control parameters of each individual;

[0104] Calculate the fitness value of each individual according to the total harmonic distortion rate index and the ice melting efficiency index, and select the individuals that meet the fitness threshold from the current population as parental individuals by using the roulette wheel selection method according to the fitness value;

[0105] Perform crossover operations on the continuous variables in the parental individuals by using the arithmetic crossover method, and perform crossover operations on the discrete variables in the parental individuals by using the single-point crossover method to generate new individuals; the continuous variables include the turning-on time and the turning-on duration; the discrete variables include the current path selection;

[0106] Perform mutation operations on the new individuals generated by the crossover, and repeat the above selection, crossover, and mutation processes until the fitness converges to a preset convergence threshold;

[0107] After the iteration process ends, select the individual with the highest fitness value from the final population to obtain the optimal ice melting control strategy.

[0108] Specifically, in this embodiment, the ice melting efficiency threshold is set according to the actual operation experience and design requirements of the DC ice melting system. The real-time obtained ice melting efficiency is compared with the preset ice melting efficiency threshold. If the ice melting efficiency is lower than the ice melting efficiency threshold, it is determined that the ice melting efficiency does not meet the standard, and the valve group control strategy needs to be further optimized; if the ice melting efficiency is higher than or equal to the ice melting efficiency threshold, it is considered that the current strategy is effective, and the current valve group control optimization strategy is continued. When it is detected that the ice melting efficiency does not meet the standard, in this embodiment, the current oscillation characteristic data (such as oscillation frequency, oscillation amplitude, oscillation intensity) and the ice melting efficiency data are used as the reference sequence, and different valve group control parameters (turn-on timing, turn-on duration, current path selection, etc.) are used as the comparison sequence. The length of each comparison sequence is the same as that of the reference sequence. For each data point, this embodiment calculates the absolute difference between the reference sequence and the comparison sequence at this point, finds the maximum and minimum values among all the absolute differences, and then calculates the grey correlation coefficient using the grey correlation degree algorithm according to the maximum and minimum values of the absolute differences. The average value of the grey correlation coefficients of each comparison sequence is obtained to get the grey correlation degree between the comparison sequence and the reference sequence. The grey correlation degrees between different valve group control parameters and current oscillation and ice melting efficiency are compared. The larger the correlation degree, the stronger the correlation between the control parameter and current oscillation or ice melting efficiency. For example, if the grey correlation degree of the turn-on timing is large, it indicates that the turn-on timing has a more significant impact on current oscillation and ice melting efficiency; conversely, if the grey correlation degree of the current path selection is small, its impact is relatively small. The key valve group control parameters that have a greater impact on current oscillation and ice melting efficiency are determined according to the correlation degree. Usually, the parameters with a correlation degree higher than the preset correlation degree threshold are selected as the key valve group control parameters. For the determined key valve group control parameters, this embodiment analyzes their change trends in historical data and the specific relationships with current oscillation and ice melting efficiency. For example, by plotting the scatter plot or line chart between the parameter and the index, observe how the current oscillation and ice melting efficiency change when the parameter increases or decreases. According to the correlation analysis results, the optimization variable range is determined for each key valve group control parameter to ensure that the parameter adjustment does not damage the system stability and can effectively improve the ice melting efficiency.

[0109] In this embodiment, several groups of valve bank control parameter combinations are randomly generated within the optimized variable range. Each combination constitutes an individual. The generated individuals are grouped together to form an initial population. The ice melting process is simulated and calculated for each individual in the population. During the simulation process, the valve bank control parameters of the individual are used. Through the simulation calculation, two optimization objectives, namely the total harmonic distortion rate index (THD) and the ice melting efficiency index, of each individual during the ice melting process are evaluated. The total harmonic distortion rate index reflects the magnitude of the harmonic content in the current waveform. According to the total harmonic distortion rate index and the ice melting efficiency index, the fitness value is calculated for each individual. The roulette wheel selection method is used to select individuals that meet the fitness threshold from the current population as parental individuals. The individuals with higher fitness values have a greater probability of being selected.

[0110] In this embodiment, different crossover methods are used to perform crossover operations on continuous variables and discrete variables to generate new individuals. In this embodiment, arithmetic crossover operations are performed on the continuous variables (such as turn-on timing and turn-on duration) in the parental individuals, that is, the corresponding variable values of two parental individuals are taken and linearly combined according to a certain ratio to generate new continuous variable values. For discrete variables (such as current path selection), single-point crossover operations are performed on the discrete variables (such as current path selection) in the parental individuals, and the values of some discrete variables are exchanged. At a randomly selected position, the gene segments of two parental individuals are exchanged to generate a new current path combination. The crossed continuous variables and discrete variables are combined together to form a new individual. Mutation operations are performed on the new individuals generated by crossover, that is, for continuous variables, a small value is randomly perturbed within its adjustment range; for discrete variables, another feasible current path is randomly selected to increase the diversity of the population. The above selection, crossover, and mutation processes are repeated to form a new generation of population. It is checked whether the fitness values of the new generation of population converge to a preset convergence threshold. If they converge, the iterative process ends; if not, the iteration continues. After the iterative process ends, the individual with the highest fitness value is selected from the final population, and the valve bank control parameters of the optimal individual are output as the optimal ice melting control strategy to guide the actual operation of the DC ice melting system. Through the above steps, this embodiment can optimize the valve bank control optimization strategy and output the optimal ice melting control strategy that can improve the ice melting efficiency and meet the harmonic distortion requirements when it is detected that the ice melting efficiency does not meet the standard, thereby improving the operation efficiency and stability of the DC ice melting system.

[0111] An embodiment of the present invention provides a method for suppressing DC ice melting current oscillation. The method includes: obtaining in real time the operation environment data of the DC ice melting system, and obtaining the equivalent inductance value and equivalent capacitance value of the DC ice melting system under the current environmental conditions according to the operation environment data; calculating the theoretical current rising rate according to the equivalent inductance value and equivalent capacitance value, and evaluating whether the current state of the current DC ice melting system meets the condition for suppressing current oscillation according to the theoretical current rising rate; when the current state of the current DC ice melting system does not meet the condition for suppressing current oscillation, optimizing the opening time, opening duration and current path selection logic of the ice melting valve group based on the equivalent inductance value and operation environment data, and generating a valve group control optimization strategy; during the process of controlling the DC ice melting system to perform ice melting operation through the valve group control optimization strategy, obtaining in real time the current waveform data of the DC ice melting system; extracting current oscillation characteristic data from the current waveform data by using the Fourier transform algorithm, and determining the ice melting efficiency according to the current oscillation characteristic data; when it is detected that the ice melting efficiency does not meet the standard, optimizing the valve group control optimization strategy by using the current oscillation characteristic data and multi-objective optimization algorithm, and outputting an optimal ice melting control strategy. Compared with the prior art, the method dynamically optimizes the valve group control strategy by real-time monitoring of the operation environment data and current waveform data, ensures that the DC ice melting system reaches the optimal ice melting efficiency while suppressing current oscillation, and guarantees the stable operation of the DC ice melting system.

[0112] It should be noted that the magnitudes of the serial numbers of the above processes do not mean the order of execution. The order of execution of each process should be determined according to its function and internal logic, and should not constitute any limitation to the implementation process of the embodiments of the present application.

[0113] In one embodiment, as Figure 2 shown, an embodiment of the present invention provides a DC ice melting current oscillation suppression system, and the system includes:

[0114] A data acquisition module 101, configured to obtain in real time the operation environment data of the DC ice melting system, and obtain the equivalent inductance value and equivalent capacitance value of the DC ice melting system under the current environmental conditions according to the operation environment data;

[0115] An oscillation evaluation module 102, configured to calculate the theoretical current rising rate according to the equivalent inductance value and the equivalent capacitance value, and evaluate whether the current state of the current DC ice melting system meets the condition for suppressing current oscillation according to the theoretical current rising rate;

[0116] A strategy generation module 103, configured to optimize the opening time, opening duration and current path selection logic of the ice melting valve group based on the equivalent inductance value and operation environment data when the current state of the current DC ice melting system does not meet the condition for suppressing current oscillation, and generate a valve group control optimization strategy;

[0117] The ice melting operation module 104 is used to obtain the current waveform data of the DC ice melting system in real time during the process of controlling and optimizing the DC ice melting system to perform ice melting operations through the valve group control optimization strategy;

[0118] The oscillation analysis module 105 is used to extract the current oscillation characteristic data from the current waveform data by using the Fourier transform algorithm, and determine the ice melting efficiency according to the current oscillation characteristic data;

[0119] The strategy optimization module 106 is used to optimize the valve group control optimization strategy by using the current oscillation characteristic data and the multi-objective optimization algorithm when it is detected that the ice melting efficiency does not meet the standard, and output the optimal ice melting control strategy.

[0120] For the specific limitations of a DC ice melting current oscillation suppression system, reference can be made to the above limitations on a DC ice melting current oscillation suppression method, which will not be elaborated here. Those of ordinary skill in the art can realize that, combined with the various modules and steps described in the embodiments disclosed in the present application, they can be implemented in hardware, software, or a combination of both. Whether these functions are executed in a hardware or software manner depends on the specific application and design constraints of the technical solution. Professional technicians can use different methods to implement the described functions for each specific application, but such implementation should not be considered to exceed the scope of the present application.

[0121] The embodiment of the present invention provides a DC ice melting current oscillation suppression system. The system obtains the operating environment data of the DC ice melting system in real time through the data acquisition module, and obtains the equivalent inductance value and equivalent capacitance value of the DC ice melting system under the current environmental conditions according to the operating environment data; the oscillation evaluation module calculates the theoretical current rise rate according to the equivalent inductance value and equivalent capacitance value, and evaluates whether the current state of the current DC ice melting system meets the conditions for suppressing current oscillation; the strategy generation module optimizes the opening time, opening duration and current path selection logic of the ice melting valve group based on the equivalent inductance value and the operating environment data when the current state of the DC ice melting system does not meet the conditions for suppressing current oscillation, and generates a valve group control optimization strategy; the ice melting operation module obtains the current waveform data of the DC ice melting system in real time during the process of controlling and optimizing the DC ice melting system to perform ice melting operations through the valve group control optimization strategy; the oscillation analysis module extracts the current oscillation characteristic data from the current waveform data by using the Fourier transform algorithm, and determines the ice melting efficiency according to the current oscillation characteristic data; the strategy optimization module optimizes the valve group control optimization strategy by using the current oscillation characteristic data and the multi-objective optimization algorithm when it is detected that the ice melting efficiency does not meet the standard, and outputs the optimal ice melting control strategy. Compared with the prior art, the system dynamically optimizes the valve group control strategy by real-time monitoring of the operating environment data and current waveform data, ensuring that the DC ice melting system achieves the optimal ice melting efficiency while suppressing current oscillation, and guaranteeing the stable operation of the DC ice melting system.

[0122] The above-described embodiments merely represent several preferred embodiments of the present application. The description thereof is relatively specific and detailed, but it should not be construed as a limitation on the scope of the invention patent. It should be noted that for those of ordinary skill in the art, without departing from the technical principle of the present invention, several improvements and substitutions can be made, and these improvements and substitutions should also be regarded as the protection scope of the present application. Therefore, the protection scope of the patent of the present application shall be subject to the protection scope of the claims.

Claims

1. A method for suppressing DC ice melting current oscillation, characterized in that It includes the following steps: Obtain the operation environment data of the DC de-icing system in real time, and obtain the equivalent inductance value and equivalent capacitance value of the DC de-icing system under the current environmental conditions according to the operation environment data; Calculate the theoretical current rise rate according to the equivalent inductance value and the equivalent capacitance value, and evaluate whether the current state of the current DC de-icing system meets the condition of suppressing current oscillation according to the theoretical current rise rate; When the current state of the current DC de-icing system does not meet the condition of suppressing current oscillation, optimize the opening time, opening duration and current path selection logic of the de-icing valve group based on the equivalent inductance value and the operation environment data, and generate a valve group control optimization strategy; During the process of controlling the DC de-icing system to perform de-icing operations through the valve group control optimization strategy, obtain the current waveform data of the DC de-icing system in real time; Use the Fourier transform algorithm to extract the current oscillation characteristic data from the current waveform data, and determine the de-icing efficiency according to the current oscillation characteristic data; When it is detected that the de-icing efficiency does not meet the standard, optimize the valve group control optimization strategy by using the current oscillation characteristic data and the multi-objective optimization algorithm, and output the optimal de-icing control strategy; Among them, the step of obtaining the equivalent inductance value and equivalent capacitance value of the DC de-icing system under the current environmental conditions according to the operation environment data includes: Perform standardization processing on the operation environment data to obtain the preprocessed operation environment data, and based on the actual circuit topology of the DC de-icing system, construct an equivalent circuit model through electromagnetic field theory; Input the preprocessed operation environment data into the equivalent circuit model, and use the least squares method to solve the equivalent inductance value and equivalent capacitance value under the current environmental conditions; The step of calculating the theoretical current rise rate according to the equivalent inductance value and the equivalent capacitance value includes: Based on the equivalent inductance value and the equivalent capacitance value, use Kirchhoff's voltage law to construct a differential equation for the change of current with time in the equivalent circuit model; Set the initial current condition at the moment of valve group opening to zero, use the fourth-order Runge-Kutta method to solve the differential equation, and extract the current time-domain waveform after the valve group is opened; Calculate the instantaneous change rate when the current starts to rise after the valve group is opened according to the current time-domain waveform to obtain the theoretical rise rate.

2. The DC ice-melting current oscillation suppression method according to claim 1, wherein The step of evaluating whether the current state of the current DC de-icing system meets the condition of suppressing current oscillation according to the theoretical current rise rate includes: Compare the theoretical current rise rate with a preset current rise rate threshold. If the theoretical current rise rate is less than or equal to the current rise rate threshold, it is determined that the current state of the DC de-icing system meets the condition of suppressing current oscillation, and control the DC de-icing system to continue de-icing operations; If the theoretical current rise rate is greater than the current rise rate threshold, it is determined that the current state of the current DC de-icing system does not meet the condition of suppressing current oscillation.

3. A method for suppressing DC ice melting current oscillation according to claim 1, characterized in that, The step of optimizing the opening time, opening duration and current path selection logic of the de-icing valve group based on the equivalent inductance value and the operation environment data, and generating a valve group control optimization strategy includes: Detect outliers in the equivalent inductance value and the operating environment data using a sliding window, and smooth the outliers to obtain a standardized equivalent inductance value and standardized operating environment data; Based on the standardized equivalent inductance value and the standardized operating environment data, use an association rule mining algorithm to analyze the correlation between the equivalent inductance value and environmental factors, and obtain the equivalent inductance distribution characteristics under the current operating environment conditions; the equivalent inductance distribution characteristics include the equivalent inductance change rate; Analyze the equivalent inductance value distribution characteristics under the current operating environment conditions, and identify the oscillation influence index of the standardized equivalent inductance value on the amplitude of the ice melting current oscillation; According to the equivalent inductance distribution characteristics and the oscillation influence index, use a two-dimensional fuzzy controller to determine the opening timing correction amount and the opening duration correction coefficient; Based on the opening timing correction amount and the opening duration correction coefficient, use the centroid method to defuzzify to obtain the ice melting valve group opening timing and the ice melting valve group opening duration; Taking the impedance value between nodes as the weight, generate a power grid structure diagram according to the actual circuit topology of the DC ice melting system; At the ice melting valve group opening timing and the ice melting valve group opening duration, use a graph theory algorithm to search for candidate current paths from the power source to the ice melting point in the power grid structure diagram; Calculate the impedance matching degree of each candidate current path, and select the candidate current path with the minimum impedance matching degree as the current optimal current path; Integrate the ice melting valve group opening timing, the ice melting valve group opening duration and the current optimal current path to form a valve group control optimization strategy.

4. The DC ice melting current oscillation suppression method according to claim 1, wherein The steps of extracting current oscillation characteristic data from the current waveform data using the Fourier transform algorithm and determining the ice melting efficiency according to the current oscillation characteristic data include: Use the Fourier transform algorithm to perform a frequency domain transformation on the current waveform data to obtain current waveform spectrum data; Identify the main current oscillation frequency components and the oscillation intensity from the current waveform spectrum data to obtain current oscillation characteristic data; According to the relationship between historical current oscillation characteristic data and historical ice melting efficiency, use a random forest regression algorithm to construct an ice melting efficiency evaluation model; Input the current oscillation characteristic data into the ice melting efficiency evaluation model to obtain the ice melting efficiency under the current operating environment conditions.

5. The DC ice melting current oscillation suppression method according to claim 4, characterized in that: The current oscillation characteristic data includes oscillation frequency data and oscillation amplitude data; the oscillation intensity is the sum of the squares of the amplitudes of all frequency components; Among them, the oscillation frequency data includes the main current oscillation frequency components; the oscillation amplitude data includes the fundamental wave component amplitude and the amplitudes of each harmonic component.

6. The method for suppressing the direct current ice melting current oscillation according to claim 1, wherein The steps of optimizing the valve group control optimization strategy using the current oscillation characteristic data and a multi-objective optimization algorithm and outputting an optimal ice melting control strategy when it is detected that the ice melting efficiency does not meet the standard include: When it is detected that the ice melting efficiency does not meet the standard, analyze the correlation between different valve group control parameters and current oscillation and ice melting efficiency according to the current oscillation characteristic data, and determine the optimization variable range of each valve group control parameter; Randomly generate individuals containing several groups of valve group control parameter combinations within the optimization variable range to form an initial population; Simulate the ice melting process for each individual in the population, and evaluate the total harmonic distortion rate index and the ice melting efficiency index of each individual during the ice melting process through the valve group control parameters of each individual; Calculate the fitness value of each individual according to the total harmonic distortion rate index and the ice melting efficiency index, and select individuals that meet the fitness threshold from the current population as parental individuals using the roulette wheel selection method according to the fitness value; Perform crossover operations on the continuous variables in the parental individuals using the arithmetic crossover method, and perform crossover operations on the discrete variables in the parental individuals using the single-point crossover method to generate new individuals; the continuous variables include the turn-on timing and the turn-on duration; the discrete variables include the current path selection; Perform mutation operations on the newly generated individuals after crossover, and repeat the above selection, crossover, and mutation processes until the fitness converges to a preset convergence threshold; After the iteration process ends, select the individual with the highest fitness value from the final population to obtain the optimal ice melting control strategy.

7. A method for suppressing DC ice melting current oscillation as described in claim 1, characterized in that The specific process for judging whether the ice melting efficiency meets the standard is as follows: Compare the ice melting efficiency with a preset ice melting efficiency threshold. If the ice melting efficiency is lower than the ice melting efficiency threshold, it is determined that the ice melting efficiency does not meet the standard; if the ice melting efficiency is higher than or equal to the ice melting efficiency threshold, continue to maintain the current valve group control optimization strategy.

8. A DC ice melting current oscillation suppression system, characterized in that, The system includes: A data acquisition module for acquiring the operating environment data of the DC ice melting system in real time, and obtaining the equivalent inductance value and the equivalent capacitance value of the DC ice melting system under the current environmental conditions according to the operating environment data; An oscillation evaluation module for calculating the theoretical current rise rate according to the equivalent inductance value and the equivalent capacitance value, and evaluating whether the current state of the current DC ice melting system meets the condition for suppressing current oscillation according to the theoretical current rise rate; A strategy generation module for optimizing the turn-on timing, turn-on duration of the ice melting valve group, and the current path selection logic based on the equivalent inductance value and the operating environment data when the current state of the DC ice melting system does not meet the condition for suppressing current oscillation, and generating a valve group control optimization strategy; An ice melting operation module for acquiring the current waveform data of the DC ice melting system in real time during the process of controlling the DC ice melting system to perform ice melting operations through the valve group control optimization strategy; An oscillation analysis module for extracting current oscillation characteristic data from the current waveform data using the Fourier transform algorithm, and determining the ice melting efficiency according to the current oscillation characteristic data; A strategy optimization module for optimizing the valve group control optimization strategy using the current oscillation characteristic data and a multi-objective optimization algorithm when it is detected that the ice melting efficiency does not meet the standard, and outputting the optimal ice melting control strategy; Among them, obtaining the equivalent inductance value and the equivalent capacitance value of the DC ice melting system under the current environmental conditions according to the operating environment data specifically includes: Perform standardization processing on the operating environment data to obtain preprocessed operating environment data, and construct an equivalent circuit model based on the actual circuit topology of the DC ice melting system through electromagnetic field theory; Input the preprocessed operating environment data into the equivalent circuit model, and use the least squares method to solve for the equivalent inductance value and equivalent capacitance value under the current environmental conditions; The calculating the theoretical current rise rate according to the equivalent inductance value and the equivalent capacitance value specifically includes: Based on the equivalent inductance value and the equivalent capacitance value, use Kirchhoff's voltage law to construct a differential equation for the current change with time in the equivalent circuit model; Set the initial current condition at the moment of valve group opening to zero, use the fourth-order Runge-Kutta method to solve the differential equation, and extract the current time-domain waveform after the valve group is opened; Calculate the instantaneous change rate when the current starts to rise after the valve group is opened according to the current time-domain waveform to obtain the theoretical rise rate.

Images