An Effective Harmonic Suppression Method Based on Linear Active Disturbance Rejection Control for Copper Electrowinning Rectifier

Because the harmonics in the production process of copper electrowinning have an important impact on the electrical energy consumption, it is necessary to suppress the harmonics effectively. In this paper, a copper electrowinning rectifier with double inverse star circuit is selected as a study object in which a large number of harmonics mainly including the 5th, 7th, 11th, and 13th harmonics are generated and injected back into the power grid. ,e total harmonic distortion rate of the power grid is up to 29.19% before filtering.,erefore, a method combining the induction filtering method and the active filtering method is proposed to carry out comprehensive filtering. Simulation results demonstrate that the total harmonic distortion rate of the system decreases to 4.20%, which indicates that the proposed method can track the corresponding changes of harmonics when the load changes in real time and filter them out. In order to ensure and improve the effect of active filter, a current harmonic tracking control method based on linear active disturbance rejection control is proposed. Simulation results show that the total harmonic distortion rate decreases to 3.34%, which is also lower than that of hysteresis control. Compared with the conventional single filtering method, the new filtering method combining induction filtering with active filtering based on linear active disturbance rejection control in the copper electrowinning rectifier has obvious advantages.


Introduction
In recent years, copper and its alloy materials are widely used in various fields, and the demand of copper is increasing in many enterprises. As one of the basic technologies of copper production, copper electrolysis and electrowinning rectification technology has attracted more and more attention. With the continuous increase of the number and capacity of copper electrowinning rectifier units, the impact of the system on the power grid in industrial application is becoming more and more serious. A large number of harmonics are injected back into the power grid, which will directly affect the quality of copper production. It is very important to ensure the safe and stable operation and power quality of power grid system [1,2]. Many researchers and scholars at home and abroad have studied power quality problems, especially the research on current quality in power grid [3][4][5][6]. e use of a large number of power electronic converter equipment will directly produce harmonic and reactive power problems and affect the power quality of the grid [7][8][9]. e copper electrowinning rectifier is such a typical device. e rectifying device of the copper electrowinning device consists of voltage regulating transformer and double inverse star rectifier, which is a typical low-voltage highcurrent nonlinear load. In the process of operation, turning on and off semicontrolled thyristor can cause a distortion of the grid current and voltage. From the perspective of frequency domain, these waveforms contain not only the power frequency sinusoidal quantity but also some voltage or current components [2] which are positive integral multiples of power frequency. ese components are called harmonics whose existence will not only affect the production quality of the system but also affect the normal operation of other equipment in the same power supply system and even lead to safety accidents in the power system. erefore, it is an urgent problem to control the harmonics to ensure the power quality and build a clean electrical environment.
Yuan et al. [10] conducted the harmonic suppression characteristic analysis of a phase-shifting reactor in a rectifier system. Takeshita et al. [11] presented the input current waveform control of the rectifier circuit which realizes simultaneously the high input power factor and the harmonic suppression of the receiving-end voltage and the source current under the distorted receiving-end voltage. Méndez et al. [12] designed an active high-power filter for the harmonic suppression in coils powered separately by 12pulsed thyristor converters. Liu et al. [13] proposed a novel controllable inductive power filtering (CIPF) method to effectively eliminate the harmonics flowing through the transformer in the industrial DC power supply system. Liu et al. [14] proposed a novel power factor correction (PFC) controller based on proportional resonant for selective harmonic suppression of a UPS system. Ramesh and Habeebullah Sait [15] introduced harmonic removal in a switched capacitor multilevel inverter utilizing the artificial bee colony (ABC) method. In [16], an impedance reconstruction control method for the source PWM inverter was proposed, which improves the phase of the output sequence impedance of the source PWM inverter at high-frequency areas to effectively suppress the high-frequency oscillation of the island power system.
Wang et al. [17] proposed a frequency-domain harmonic model for 12-pulse series-connected thyristor-controlled rectifier under unbalanced supply voltage. Wiechmann et al. [18] proposed an optimized sequential control technique to improve converter's reactive power consumption presentation for copper electrowinning high-current rectifiers. In [19], a modified DPWM (MDPWM) scheme was proposed to regulate the neutral point voltage by using redundant clamping modes in high-frequency Vienna-type rectifiers. An et al. [20] simulated diode rectifier + FC filter schemebased large electrolytic aluminum filter device. Mao et al. [9] used the filter with single tuning in a high-power thyristor rectifier as a medium-frequency furnace and high-frequency furnace for metal smelting. Dai et al. [21] proposed a comprehensive treatment scheme of "automatic switch single tuning passive filter + high voltage side series filter inductor" in medium frequency smelting furnace. Ma et al. and Li et al. [22,23] applied the first-order linear active disturbance rejection controller (LADRC) to the current tracking control of a 10 kV line parallel hybrid active power filter, which can well deal with the contradiction between overshoot and rapidity, and the performance is better than the traditional PI controller, which verifies the effectiveness of the linear active disturbance rejection control strategy. In [24], aiming at the disadvantages of the grid side inverter of doubly fed induction generator, such as the complexity and poor stability, and the problems of the conventional PI closed-loop control, such as poor stability and slow response speed, an ADRC control strategy was proposed, in which a resonance link is added to the linear state error feedback. Instead of PI control of current loop, it can suppress the loworder harmonics of grid connected current, improve the response speed of the system, and ensure smooth operation when the inverter is connected to the grid.
Analyzing above relevant researches on harmonic suppression methods for rectifier indicates that there are few special effective researches on harmonic suppression of high-power copper electrowinning rectifier. e effect of the conventional harmonic suppression method is not ideal. It is very important to study an effective harmonic suppression technology. erefore, in this paper, an "induction filter-+ active filter" method is proposed to effectively suppress the harmonics. Moreover, the LADRC technology is added to the active filter current tracking control in active filter, which further improves the harmonic suppression effect.

Copper Electrowinning Process.
In copper electrowinning production, cathode copper is produced through the leaching-extraction-electrowinning process; in the electrolyte, the insoluble lead-alloy plate (Pb-Ca-Sn) is used as the anode, and the thin copper starting sheet is used as the cathode. In essence, copper electrowinning is the process where copper ion is reduced by DC and deposited on the cathode [25,26]. Main reactions are as follows.
e copper sulfate solution C u SO 4 is decomposed into a copper ion Cu 2+ and sulfate ion SO 2− 4 in the applied electric field, and the process is expressed as e copper ion obtains electrons on the cathode and is reduced to copper which is deposited on the cathode, and the process can be expressed as

Current Efficiency.
According to Faraday's law, the theory is that copper of 1.1864 × 10 − 3 kg should be deposited on the cathode by DC of 1 ampere-hour. In fact, compared with the theoretical value, the actual amount of produced copper is decreased by impurities, oxidation and dissolution of cathode sediment, and electrode short circuit and leakage loss. In production practice, current efficiency is proposed to evaluate effective utilization of current, and it is the percentage ratio between the actual production of copper and theoretical production of the copper and is expressed by [27] η � G q · I · t · n × 100%, where η represents current efficiency in %, G represents actual production of copper in kg during t, q represents electrochemical equivalent of copper in 1.1864 × 10 − 3 kg/Ah, I represents current intensity in A, t represents the time of electrowinning in hour, and n represents the number of electrowinning electrolytic tanks. Current efficiency is an important technical and economic index of electrowinning copper production. In production practice, due to different specific conditions, the current efficiency of copper electroproduction is different, and effective measures should be taken to improve current efficiency.

DC Power Consumption.
In copper electrowinning production, DC power consumption in kWh/kg is defined as the electricity used to produce copper of 1 kg in 1 hour, which is an important economic and technical index. DC power consumption is calculated by where W represents DC power consumption in kWh/kg, v represents electrolytic tank voltage in V, and η represents current efficiency in %. Formula (4) indicates that electrolytic tank voltage and current efficiency are the direct factors affecting energy consumption. In the industrial process, the voltage of conventional electrolytic tank voltage is in an approximate range of 1.8V ∼ 2V, and the electrolytic tank voltage will increase with the increase of current density; consequently, the power consumption will also increase. However, if lowcurrent-density electrowinning is used, the output will be reduced, although the power consumption can be reduced in the electrowinning process. erefore, the actual production situation of the factory should be taken into account when a current density is determined. At present, the current density of copper electrowinning used in the world is different, fluctuating in a range; the low one is approximately 100A/m 2 , and the high one is approximately 270A/m 2 .

Double Inverse Star Rectifying Circuit.
With the continuous development of electrolysis, electroplating, electrolytic degreasing, and other industries, the demand for low-voltage and high-current DC power supply is increasing. Double inverse star rectifying circuit stands out among all kinds of circuits for its superior performance. It is a rectifying circuit composed of two three-phase half-wave circuits and a balance reactor, which is depicted in Figure 1. Compared with the common six-phase half-wave rectifier circuit and three-phase six-pulse bridge rectifier circuit, the current effective value of the rectifier arm is smaller, and the capacity utilization rate is high. Adding the balance reactor can improve the reverse voltage drop, and the rectifier can parallel two output different instantaneous values of groups of three-phase half wave to supply power to the load. In the working process of double inverse star rectifier circuit, the rectifier is a nonlinear load, which will inevitably produce harmonics. e internal structure diagram of double inverse star transformer is illustrated in Figure 2. On the primary side of the transformer, the windings are connected according to the triangle structure, and the secondary side is composed of two-star structure windings with opposite directions. On the premise of meeting the safety requirements of industrial application, Figure 3 shows that the input three-phase grid voltages u a , u b , u c are changed into six voltages u a1 , u b1 , u c1 and u a2 , u b2 , u c2 and there is a difference of 60°between them.
According to the internal structure and voltage vector relationship of double inverse star rectifier transformer, the corresponding mathematical model can be established. Firstly, the input three-phase grid voltage is set as where U m is defined as the magnitude of the grid voltage. According to the voltage vector relationship shown in Figure 3, the other two voltage equations can be obtained as follows: where U n is the amplitude of secondary voltage of the transformer. Assuming that the ratio of variation is k, the relationship between U m and U n satisfies 3.2. Analysis of Characteristic Subharmonic Current. In order to analyze the types and contents of harmonics generated by the double inverse star rectifying transformer in the working process and effectively filter out the harmonics, it is assumed that both the grid voltage and the transformer are ideal. It can be seen from Figures 1 and 2 that the primary side current of the transformer can be represented by the secondary side current as e system input currents are obtained as By substituting (9) into (10), the expressions of system input current and secondary side current are obtained as According to (11), it is difficult to analyze the input current. In this study, the load will be simplified to a constant current. It is difficult to determine the instantaneous voltage difference of three-phase half-wave output I d , but its waveform is a triangular wave with three times the fundamental frequency. Since the balance reactor can make the rectifier system work independently with symmetry, the current can be supplied by two transformers equally, namely, In this paper, the switching function is introduced to study the current relationship between AC side and DC side directly. Taking the phase a 1 as an example, its switching function is S a1 � i a1 /i d1 whose Fourier series expansion is According to S a1 , Fourier series of phases b 1 and c 1 can be obtained, and then the expression of three-phase halfwave group I is expressed as On the basis of (14), the expression of three-phase halfwave group II can be obtained as e relationships between the switching function of each phase and its corresponding current input and output are as follows:

Mathematical Problems in Engineering
According to formulas (13)- (17), the Fourier series decomposition of each secondary side current is as follows: It can be seen from (12), (18), and (19) that taking phase a as an example, its current expression is It is obvious that the harmonic content in the power grid mainly contains (6k ± 1) times (k � 1, 2, 3, ...) harmonics. e electrowinning copper rectifier used in this paper is a double inverse star rectifying system, and its system parameters are summarized in Table 1.
According to Table 1, the parameters of the double inverse star rectifying simulation circuit are set. Figure 4 shows waveforms of output voltage and output current. Figure 5 shows waveforms of voltage and current measured on the power grid. Figure 6 shows the spectrum analysis diagram of harmonic measurement on the power grid.
Spectrum of harmonics in Figure 6 demonstrates that a large number of harmonics are generated and injected into the power grid in the rectifying process of double inverse star rectifier, and these harmonics are mainly the 5th, 7th, 11th, and 13th harmonics. Before filtering, the total measured distortion rate of the power grid is up to 29.19%.

Filtering Method Combining Active Filter with Inductive Filter
Because the induction filtering method can meet the requirements of the total harmonic distortion rate after filtering and the active power filtering method can meet the requirements of tracking and filtering the harmonic changes in real time after the load changes, the two filtering methods are combined to carry out comprehensive filtering for the copper electrodeposited rectifying system. A system diagram of combining the two filtering methods is shown in Figure 7, an induction filter is added to the transformer side, and the active power filter is added to the power grid side, which will make the electrowinning copper rectifier achieve a good filtering effect.
In the simulation process, the active filter is added at 0.02 s, the inductive filter is added at 0.06 s, and the SCR trigger angle is changed from 10°to 20°at 0.12 s, that is, when the load changes, the corresponding harmonic also changes.
e simulation results are shown in Figures 8  and 9. e corresponding harmonic spectrum is shown in Figure 9.
From the above simulation results, it can be seen that when the active filter is added at 0.02 s, the waveform has obvious changes and the filtering effect is good. When the induction filter is added at 0.06 s, the waveform is stable after one cycle, that is, the total harmonic distortion rate is 4.83% at 0.08 s, which basically meets the requirements of the power grid. When the load changes at 0.12 s, the total distortion rate rises to 8.07%, which indicates that the harmonic has changed. After a period of 0.16 s, the total distortion rate of the system decreases to 3.76%, which indicates that the filtering system can track the corresponding changes of harmonics when the load changes in real time and filter them out.

Construction of LADRC Applied in APF.
In this paper, the current tracking control method of instantaneous value comparison is adopted. Although it has an obvious filtering effect, its tracking effect is limited by the choice of hysteresis loop width. If it is too large, the tracking effect is not good. If it is too small, the switching frequency is too high which results in higher power loss. erefore, a new current tracking control method-linear active disturbance rejection control (LADRC)-is introduced. e LADRC technology is an improvement of linearization based on active disturbance rejection control (ADRC). e key part is the linearization of the extended state observer (ESO). Compared with ADRC, the parameter tuning of ACRC is reduced, and the control process is relatively simple, which improves its practical value in the control process.
In this paper, the LADRC technology is applied in current tracking control, which means that ESO and control are linear, and the differential tracker is not needed in   Mathematical Problems in Engineering practical application. e control block diagram of LADRC application in APF is shown in Figure 10, taking singlephase a as an example. Figure 11 shows the internal structure block diagram of LADRC.

Design of Linear Expansion State Observer (LESO).
In the entire LADRC system, the LESO is the most critical part. It can calculate each state quantity and disturbance value of the unknown model system according to the input data. For any given n-order unknown system, there is where y(t) represents the input signal of the system, the nonlinear dynamic characteristic of the system is f(•), w(t) represents the external disturbance of the system, u(t) represents the control input signal of the system, and b 0 represents the gain of the controller. In LACRC, f(•) is the expression of total disturbance. All state spaces of the system satisfy where, the state variable x n+1 , which is the total disturbance passing through the LESO, can be expressed as In this way, the original n-order system is expanded into a new n + 1-order linear system. e first-order linear ADRC is designed in this paper, and the state space of the system can be written as X is added to the system as an expanded state. For the convenience of design, the matrix form of (26) is where • is the first derivative of the disturbance. e LESO of the system is written as where L � β 1 β 2 is the gain of LESO.
If there is no steady-state error, the eigenvalues of the observation error characteristic matrix should be all in the left half of the plane, and the characteristic equation can be obtained as When the bandwidth of the observer is fixed, if the following conditions are satisfied: the effect of the controller is best. According to formulas (30) and (31), e specific form of LESO is Inputs of LESO are u(t) and y(t). e two system states z 1 and z 2 track y and expand the state x 2 , respectively.

Design of Linear State Error Feedback (LSEF).
e main function of linear error feedback control rate is not only to preserve the advantages of integral link but also to avoid the error caused by integral transformation when designing parameters. erefore, the disturbance compensation value −z 2 /b 0 is required to replace the integral part of PID control in the design process. For the second-order and even higherorder system, although there is no integral term, it can still reflect the function of the integral. e LADRC system used in current tracking control is a first-order system in this paper, which does not need differential tracker. erefore, we only need to use the control rate for P.
where u 0 is the error feedback control quantity and k p represents the proportional constant. It can be seen from Figure 11 that the final error feedback control quantity is For the convenience of frequency-domain analysis, Figure 11 can be transformed into Figure 12. Figure 12 shows the transfer function in the frequency domain. e input signal of the controller is R(s), and Y(s) indicates the output signal. e disturbance signal is W(s), and P(s) indicates the controlled object model. e control signal is U(s), and H(s) and G(s) are the undetermined terms of the controller, respectively. e control signal can be expressed as

Convergence Analysis of Linear Extended State Observer.
Taking the controlled object represented by equation (22) as an example, for the n-order controlled object represented by equation (23), the LESO is represented by e parameter of the observer is L � [B 1 , B 2 , . . . , B n+1 ] which is adjusted by pole assignment method.
at is to say, the pole assignment method is used to design the convergence of the LESO.

Harmonic Suppression with LADRC.
On the basis of the filtering method combining active filter with inductive filter, the LADRC technology is applied in current tracking control in ARF, filtering effects of which are obtained by conducting simulation. Simulation results are obtained in Figures 13 and  14. e active filter is added at 0.02 s, the inductive filter is added at 0.06 s, and the SCR trigger angle changes from 10°t o 20°at 0.12 s, that is, the load changes.
Comparison of total distortion rates between hysteresis control and LADRC is summarized in Table 2.
It can be seen from Table 2 that obviously the total harmonic distortion rate under the control of LADRC is 3.99% at 0.08 s, which is lower than that under hysteresis control; at 0.12 s, the trigger angle of SCR changes from 10°to 20°, that is, when the load changes, the harmonic also changes accordingly. After one cycle, the total harmonic distortion rate returns to 3.41%, which is also lower than that of hysteresis control. erefore, the results verify that the harmonic suppression effect of LADCR is better than that of hysteresis tracking control.

Conclusions
Since the induction filtering method can meet the requirements of the total harmonic distortion rate after filtering and the active power filtering method can meet the requirements of tracking and filtering the harmonic changes in real time when the load changes, in order to effectively suppress a large number of harmonics mainly including the 5th, 7th, 11th, and 13th harmonics generated during operation of a copper electrowinning rectifier with double inverse star circuit, firstly, Firstly, this paper proposes a method combining the induction filtering method and the active filtering method to carry out comprehensive filtering. Simulation results indicate that the total distortion rate of the system decreases to 4.83% from 29.19% before filtering, which indicates that the proposed method can track the corresponding changes of harmonics when the load changes in real time and filter them out. Secondly, a current harmonic tracking control method based on LADRC is proposed to ensure and improve the effect of the active filter. Simulation results show that the total harmonic distortion rate can be decreased to 3.41%, which is also lower than that of hysteresis control. Compared with the conventional single filtering method, it is obvious that the proposed filtering method combining induction filtering with active filtering based on LADRC in the copper electrowinning rectifier has obvious advantages in setting parameter and debugging system.

Data Availability
e data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest
e authors declare that there are no conflicts of interest regarding the publication of this paper.