Integral Sliding Mode Control for Maximum Power Point Tracking in DFIG Based Floating Offshore Wind Turbine and Power to Gas

This paper proposes a current decoupling controller for a Doubly-fed Induction Generator (DFIG) based on floating offshore wind turbine and power to gas. The proposed controller realizes Maximum Power Point Tracking (MPPT) through integral sliding mode compensation. By using the internal model control strategy, an open-loop controller is designed to ensure that the system has good dynamic performance. Furthermore, using the integral Sliding Mode Control (SMC) strategy, a compensator is designed to eliminate the parameter perturbation and external disturbance of the open-loop control. The parameters of the designed controller are designed through Grey Wolf Optimization (GWO). Simulation results show that the proposed control strategy has better response speed and smaller steady-state error than the traditional control strategy. This research is expected to be applied to the field of hydrogen production by floating offshore wind power.


Introduction
Since the industrial revolution, energy has become the fundamental power for human production and development. On the one hand, natural reserves of non-renewable resources such as oil and coal are limited. After being consumed by humans for centuries, such resources are gradually depleted. On the other hand, after the combustion of various fossil fuels, a large amount of toxic gas (including nitrogen oxides and sulfur oxides) and particulate matter will be emitted, which will increase the risk of respiratory diseases [1]. Under such circumstances, various renewable energy utilization schemes such as wind energy and solar energy have been proposed [2]. Due to its unique advantages, Doubly-fed Induction Generators (DFIG) occupy a large proportion in the field of wind power generation. The performance of the drive control system is directly related to the operating performance of the generator [3]. Therefore, the generator control strategy is the key to affecting generator performance.
Including doubly-fed induction generators, the induction motors are multi-variable and strongly coupled systems. The vector control is the strategy, which transforms the three-phase asynchronous motor system into the DC motor system through coordinate transformation, thereby achieving independent control of torque and flux components [4]. The vector control will overlap the coordinate system with a certain vector of the motor.
They are the available vectors, such as stator flux, stator voltage, air gap flux, rotor flux, etc. [5]. Generally, a DFIG wind power system using vector control is designed as a double closed loop structure. The two PI controllers are used to control the speed and reactive power separately in the outer loop. At the same time, the two PI controllers are used to control the d-axis current component and the q-axis current component, respectively, in the inner loop [6]. In [7], the author proposed a robust proportional-integral (PI) controller with its parameters designed by constrained population extremal optimization for the load frequency control problem. Zeng proposed a hybrid method combined adaptive control, neural network, and PID and proved its effectiveness for a nonlinear system in [8]. For this control structure, some studies have adjusted the parameters of the PI controller by applying a heuristic algorithm to obtain better power capture performance [9,10]. Qi investigated the event-triggered L ∞ control for network-based switched linear systems with transmission delay in [11] and proposed H ∞ filtering for networked switched systems in [12].
With the development of floating offshore wind turbines over the world, their reliability and dynamic performance requirements have gradually increased [13,14]. Although vector control theory can be used to convert a doubly-fed floating offshore wind turbine system into a linear system. However, in fact, the doubly-fed wind turbine system is still nonlinear, and it often contains a large amount of time delay and coupling components in the actual working environment [15][16][17]. With the increase of coupling effect, its torque may produce transient distortion, which affects the dynamic performance of the system, and then causes the maximum power point to not coincide with the actual working state [18]. Therefore, it has been difficult for the simple PI control system to meet the increasingly high work requirements, and more advanced controllers need to be designed according to modern control technology. It is worth mentioning that floating offshore wind energy develops rapidly in developing countries in recent years, such as China.
Floating offshore wind power generation technology provides new ideas for solving the problem of wind curtailment. It is of great significance for solving the problem of floating offshore wind power consumption, and development of decentralized wind power generation technology, realizing multi-channel and efficient use of renewable energy [19]. Wind power generation hydrogen is the method by using electrolysis of water energy storage. On the one hand, it can integrate hydrogen as a clean and high-energy fuel into the existing gas supply network and realize the complementary conversion of electricity to gas. On the other hand, it can directly use hydrogen energy in high-efficiency and clean technologies such as fuel cells. For example, Figure 1 shows the description of floating offshore wind turbines and power to gas [20]. Usually, the floating offshore wind turbines are located in the deep sea and have the advantages of abundant wind resources.
Sliding mode control is a nonlinear control method, which is characterized by discontinuous control. It can make the system switch from one continuous structure to another continuous structure in the state space, so this is a variable structure control method. The control design of sliding mode control will change the dynamic characteristics of the system so that the state trajectory moves along the boundary of the control structure. This boundary is called the sliding surface. The sliding surface can be designed as required, and the movement of the system trajectory on the sliding surface is not affected by the parameters of the controlled object and external interference. Therefore, sliding mode control has the advantages of fast response and insensitivity to external interference [21].
In recent years, with the development of computer technology and high-power electronic switching elements, the realization of sliding mode control has become easier. For such complex nonlinear systems as wind power systems, it is suitable to use sliding mode control theory to design controllers that meet performance requirements. It can improve the dynamic performance of wind turbines by using the sliding mode controller instead of the PI controller for speed control [22,23]. In Ref. [24], SMC is used for direct power control of a doubly-fed induction generator wind turbine operating in unbalanced grid voltage conditions. In Ref. [25], an improved strategy of SMC is used for the controller design of a variable speed wind turbine with a fractional order doubly fed induction generator. In Ref. [26], a combination of sliding mode control and phase-locked loop estimator is used for the sensorless vector control strategy of a rotor-tied doubly fed induction generator system. In Ref. [27], an adaptive gain second-order sliding mode direct power control strategy is proposed, which can work normally in a DFIG wind turbine system with unbalanced grid voltage. In view of the rich characteristics of floating offshore wind energy resources, this study focuses on high-power floating offshore wind power generation systems. Based on the above studies, in order to solve the problem of insufficient dynamic performance and robustness of traditional dual closed-loop PI controllers, this paper proposes a new MPPT control strategy. According to the principle of internal model control, an internal model open-loop controller is designed. In addition, in order to improve the controller's ability to resist external disturbances, the sliding mode theory is used to compensate the errors of the internal-mode open-loop controller. The novelty of this study is the use of a new sliding mode controller structure. In the proposed control method, a sliding mode controller is used to compensate the output error instead of directly outputting the control signal. This is different from most DFIG wind turbine systems based on sliding mode control. In order to verify the accuracy and reliability of the proposed controller, the proposed controller was simulated and verified in Matlab software.
The structure of the article is as follows: Section 2 describes the DFIG floating offshore wind turbine system model. Section 3 introduces the design of the proposed controller and gives a detailed proof of system stability. In Section 4, the proposed control method is simulated, and the simulation results are analyzed. Finally, it is summarized in Section 5.

DFIG Floating Offshore Wind Turbine System Modeling
The doubly-fed floating offshore wind turbine system is shown in Figure 2. The generator used in this power generation system is a Doubly-Fed Induction Generator [28]. The floating offshore wind turbine is connected to the generator rotor through a gearbox. At the same time, the rotor is excited by a power converter, and the generator stator is connected to the grid. A variable-speed constant-frequency power generation system using a doubly-fed floating offshore wind turbine is one of the popular forms of current wind energy conversion systems. Decoupling control of active and reactive power in floating offshore wind turbine systems can be achieved by applying dual PWM converters. This kind of wind energy conversion system has the regulation capability similar to that of a synchronous motor, and can regulate the voltage and frequency of the power grid. In addition, the power converter is located on the rotor side, and the power of the rotor side circuit depends on the generator slip power. Therefore, its actual power is only about 30% of the generator capacity, which greatly reduces the size and cost of the power converter. Due to such advantages of the doubly-fed floating offshore wind turbine system, the industry has conducted in-depth research on the converter control strategy of this system.

Floating Offshore Wind Turbine Model
In the doubly-fed floating offshore wind turbine system, the floating offshore wind turbine is used to convert wind energy into mechanical energy, and the generator is used to convert mechanical energy into electrical energy. The mechanical energy is transmitted between the floating offshore wind turbine and the generator through a transmission structure such as a transmission shaft and a gearbox. The total amount of mechanical energy output by a floating offshore wind turbine is mainly determined by the wind speed and the aerodynamic characteristics of the floating offshore wind turbine. At the same time, the mechanical energy characteristics delivered to the generator are also closely related to the mechanical characteristics of the transmission system. The wind energy captured by an floating offshore wind turbine is usually described by Equation (1) [29]: where P m is the output power of the floating offshore wind turbine, C p is the power coefficient, and it is related to the tip speed ratio λ and pitch angle β, ρ is the air density, A is the area swept by the pales of the turbine, and v is the wind speed. The power coefficient can be written as Equation (2) [30]: The definition of tip-speed-ratio λ is as the follows. The ratio of the linear velocity of the floating offshore wind turbine blade tip to the wind speed is λ, that is, λ = ωR/v, ω is the angular speed of the floating offshore wind turbine, R is the blade length, and c i (i = 1, 2, 3, 4, 5, 6) are constants related to the characteristics of the floating offshore wind turbine. The parameters used in this study are c 1 = 0.22, c 2 = 116, c 3 = 0.4, c 4 = 5, c 5 = 12.5, and c 6 = 0. The relationship between power coefficient, tip speed ratio, and pitch angle is shown in Figure 3. It can be known from the above analysis that there is an optimal tip speed ratio λ opt , which makes the captured wind energy the largest by the floating offshore wind turbine at this tip speed ratio. Therefore, the reference value of the rotational speed of the floating offshore wind turbine can be given by Equation (3):

DFIG Model
The mathematical model of an induction motor has the characteristics of multivariable, strong coupling, nonlinearity, and time-varying. It is very inconvenient to analyze it. The original model can be transformed by vector control theory to obtain the state equation of DFIG in the synchronous rotating coordinate system [28]. Its dynamic characteristics can be described by Equations (4) and (5): where u sd and u sq are voltage components of the stator in d-axis and q-axis, u rd and u rq are voltage components of the rotor in d-axis and q-axis, i sd and i sq are current components of the stator in d-axis and q-axis, i rd and i rq are current components of the rotor in d-axis and q-axis, R s and R r are the resistances of the stator and rotor, ψ sd and ψ sq are the flux components of the stator in d-axis and q-axis, and ψ rd and ψ rq are the flux components of the rotor in d-axis and q-axis. ω 1 is the synchronous electromagnetic angular velocity, ω r is the rotor electromagnetic angular velocity, p is the differential operator, T e is the electromagnetic torque, n p is the number of pole pairs, T L is the external load, J is the moment of inertia, and D is the damping coefficient. The relationship between flux and current is shown in Equation (6): where L s is the stator inductance, L r is the rotor inductance, and L m is the mutual inductance.
In stator flux orientation vector control, the direction of the stator flux is coincident with the d-axis of the synchronous rotation coordinate system. Assuming that the DFIG winding is directly connected to an infinitely large power grid, parameters such as the voltage amplitude, phase, and angular frequency of the grid can be considered as the constant, and the stator flux can be regarded as a fixed value. If the stator resistance is ignored, ψ sd = ψ s and ψ sq = 0. One has the following results in Equations (7) and (8): Substituting Equations (7) and (8) into (4), we know that where σ is the leakage coefficient, and its definition is shown in Equation (10): It can be known from the above analysis that the decoupling of active power and reactive power of the stator has been realized by vector control. Generator stator active power can be controlled by the rotor current q-axis component, and stator reactive power can be controlled by the rotor current d-axis component. Considering that the rotor voltage balance is not completely decoupled, the feedforward term is defined as shown in Equation (11): where u rd and u rq are voltage feedforward components of the rotor in d-axis and q-axis. The d-axis and q-axis current component can be decoupled, and the entire system can be regarded as a first-order inertial system. The transfer function of the system is shown in Equation (12): where u rd0 and u rq0 are voltage components of the rotor in d-axis and q-axis after decoupling. Equation (12) shows that, if the grid voltage is constant, the electromagnetic torque of DFIG is proportional to the current component of the rotor in q-axis, that is, the electromagnetic torque can be controlled by the rotor q-axis current.
In summary, the DFIG equivalent model under stator flux orientation vector control is shown in Figure 4. By controlling the q-axis component of the rotor current, the active power of the DFIG can be controlled. By controlling the d-axis component of the rotor current, the DFIG reactive power can be controlled. Therefore, controlling the d-axis and q-axis components of the DFIG rotor current can control the output power of the generator. Since the electromagnetic torque can also be controlled by the rotor current, the DFIG speed control can be further realized.

Internal Model Controller
Internal model control is based on the process mathematical model for the controller design. It was proposed and applied to wind turbine in 2018 [31]. This control method has a good effect on solving tracking problems and anti-jamming problems, and it also has robustness to model uncertainty problems. Its basic control structure is shown in Figure 5.
In Figure 5, where G(s) is the system model transfer function, M(s) is the internal model transfer function, C(s) is the controller transfer function, F(s) is the internal model control transfer function, R(s) is the reference input signal, E(s) is the error signal, U(s) is the control system output signal, D(s) is the interference signal, and Y(s) is the system output signal. According to Figure 5, the system input-output transfer function is shown in Equation (13).
It is assumed that the internal model can be fully described by the system model, that is, M(s) = G(s). When C(s) = M −1 (s), there is Y(s) = R(s). In this ideal state, it is not necessary to adjust any control parameters to make the output of the system fully track the reference input. In addition, since the system output does not include interference signal components, the controller can be considered to be able to overcome any interference signals.
However, in an actual system, the above ideal controller is difficult to implement for the following reasons: (i) When the controlled system G(s) contains time-delay components, the controller C(s) will contain lead components, which does not conform to physical reality and is difficult to achieve; (ii) When the right half-plane zeros exist in the controlled system G(s), the controller C(s) will contain the right half-plane poles, and the controller will be unstable at this time, which will affect the stability of the entire system; (iii) When the order of the denominator polynomial in the controlled system G(s) is higher than the numerator, the controller C(s) will include a differentiator. Since the differentiator is extremely sensitive to signal noise, it is not suitable for practical use; (iv) When there is a model error, that is, M(s) = G(s), the ideal controller will not guarantee the stability of the system.
Due to the above reasons, there are two steps in the designing of the controller. First, the controlled system G(s) is been decomposed. The smallest phase part has been found that does not include the time delay components and the zero point in the right half plane. The minimum phase part is the internal mode M(s) of the control system. Then, the controller is constructed by inverting the internal model and adding a filter, which is shown in Equation (14): where L(s) is a low-pass filter, as shown in Equation (15): where α is a modulation coefficient. Combined with Equation (12), the internal model controller can be obtained as shown in formula (16): The controller proposed according to the internal model theory is not difficult from the PI controller in form, but the internal model controller has only one parameter that needs to be adjusted. If the internal model parameters are completely consistent with the actual system parameters, the actual input signal can be replaced by the internal model's response to the output signal. Giving the definition of the system output estimateŶ(s), which is shown in Equations (17) and (18) An open-loop internal model controller is obtained, which is shown in Figure 6.

Integral Sliding Mode Compensation
Due to external disturbances in the actual system, the model parameters may not be accurate enough. The estimated current of the open-loop controller will deviate from the actual current. Therefore, in order to improve the robustness of the control system, a sliding mode control strategy is introduced. The tracking error function is defined as shown in Equation (19): The open-loop controller is designed based on the internal model control theory. It can achieve good performance when the system reaches the equilibrium position, but when the system state is unbalanced at the initial moment and there is external interference, the internal-mode open-loop controller cannot obtain good control performance. For such applications, the integral sliding mode theory can be introduced to design the integral sliding mode surface. According to the integral sliding mode theory, the sliding mode surface is usually defined according to Equation (20): where r is the system order, δ is a positive real number, and k i is a sliding gain. It can be known from Equation (18) that r = 1 and δ = 1. Integral sliding mode control adds an integral part to the basic sliding mode control method to improve the performance of sliding mode control. Sliding mode control is robust to disturbances only in the sliding mode phase, but it is not robust in the arrival phase. Since the integral sliding mode control includes an integral part, a sliding surface can be designed so that the system is in the sliding mode phase at the initial moment. Thus, the entire state trajectory of the system is on the sliding surface [32]. The designed sliding surface is shown in Equation (21): The current estimate satisfies Equation (22): Substituting Equation (22) into Equation (9), the derivative of the sliding mode surface is shown in Equation (23): The stability of sliding mode control can usually be judged by Lyapunov function, as shown in Equation (24): In order to ensure that Lyapunov function is negative definite, that is,V ≤ 0, the control law is designed as shown in Equation (25): Substituting the above Equation ( Therefore, as long as the constant > 0 is satisfied,V ≤ 0 can be ensured. Considering the impact of uncertainties and disturbances on the system, the sliding mode surface can be rewritten into the form of Equation (27): where τ d is the influence of the uncertainties and disturbances on the sliding surface. Equation (26) can be rewritten as:V = S T (− sign(S) + τ d ) (28) As long as > |τ d | is satisfied,V ≤ 0 can be ensured. This shows that the proposed control method can ensure stability under disturbances.

Remark 1.
According to the characteristics of sliding mode control, the state trajectory of the system will be limited to the sliding surface after the system state moves to the sliding surface. However, in practical systems, due to the inertia of the controlled system and the delay of the switching components, the switching frequency of sliding mode control cannot be infinite. This will cause the trajectory of the system state to not stay on the sliding surface, but to traverse back and forth on the boundary of the sliding surface or to perform periodic motion near the equilibrium point. This high-frequency oscillation phenomenon is called chattering. Chattering affects the control accuracy of the system and cause frequent switching of switching components, shortening the component operation cycle and increasing power consumption.

Remark 2.
For floating offshore wind turbine systems, if chattering occurs in the DFIG drive control system, it may cause speed fluctuations, torque ripples, and current harmonics. If chattering occurs in the voltage or flux observer, it will cause observation noise and reduce system stability. Although sliding mode control has the advantage of good robustness, to some extent, chattering phenomenon hinders its control advantage. Therefore, the suppression of chattering has become the primary issue in practical application research of sliding mode control.
In the sliding mode control system of this study, due to the introduction of integral sliding mode control, the initial state of the system can be located on the sliding surface. Therefore, the design of the approach law should mainly be considered as the suppression of chattering. Based on the above analysis, the control law is shown by Equation (29): tanh(γS(t)) where 0 < < 1, η > 0, and γ > 0, tanh(S) is a hyperbolic tangent function, and its expression is shown in Equation (30): The relationship between the sliding function and its derivative under different parameter values is shown in Figure 7.
It can be known from Figure 7 that the approaching speed of the system is related to the current value of the sliding function. The larger the value of the sliding function, the faster the approaching speed. When the state of the system is close to the sliding surface, that is, the sliding function approaches 0, and the properties of the law of approximation approximate the hyperbolic tangent function. At this time, the characteristics of the system are basically determined by parameter γ. The hyperbolic tangent function is a continuous odd function. When it is far from the zero point on the positive semi-axis, its function value is approximately 1. When it is close to the zero point, its function value will rapidly decrease to 0. This characteristic allows the approach law to retain the robustness of the sliding mode control while suppressing chattering. When the parameter is close to 1, the characteristic of the law of approximation over the entire state space is similar to the hyperbolic tangent function. That is, when the system state is far from the sliding surface, the approaching speed is approximately constant. When the system state is close to the sliding surface, the approaching speed is rapidly reduced to 0. When the parameter is close to 0, the approaching speed farther from the sliding surface approximates an exponential relationship with the value of the sliding function. Therefore, the value of affects the characteristics of the approaching law when it is far away from the sliding surface, that is, the reaching phase of sliding mode control. The parameter η affects the transition process of the two characteristics of the system state when it is far away from and close to the sliding surface. When the value of this parameter is larger, the approach speed will converge closer to the sliding surface.  (29) and (30), the control law of the improved sliding mode controller is shown in Equation (31):

Theorem 1. According to the analysis in Equations
Proof. To verify the stability of the proposed system, substituting the control law (31) into the Lyapunov function described in Equation (24), we can geṫ Since 0 < < 1, η > 0.
It can be concluded thatV The equal sign holds if and only if S = 0. Therefore, it can be proved that the proposed sliding mode controller meets the stability requirements.
Taking the current d-axis component as an example, the structure of the proposed internal model-integral sliding model controller(IM-ISMC) is shown in Figure 8. As we can see in the figure, in the block diagram of the internal model-integral sliding model controller, an improved approach law has been used on u i . Moreover, the tracking error calculation has been used on error e, and the final output is u rd . In summary, by using the IM-ISMC controller instead of the PI controller in vector control, an IM-ISMC-based floating offshore wind turbine system can be obtained, as shown in Figure 9. As we can see in the figure, an IM-ISMC controller has been used at the back of PI controller. Moreover, the coordinate transformation has been used on PWM and IGBT, and angle calculation has been used at the back of DFIG. The proposed control system is based on vector control theory. The relationship between the q-axis component of the rotor current and the speed can be better described by a linear system, and it can be controlled by a PI controller. The relationship between the d-axis component of the rotor current and the reactive power is similar. The current is decoupled by the IM-ISMC controller.
The proposed control strategy contains many undetermined parameters, which include modulation coefficient α, sliding gain k i , and related parameters of sliding mode control law k p , , η, and γ. The Grey Wolf Optimization (GWO) will be used to optimize these parameters.
Grey Wolf Optimization is a meta-heuristic algorithm proposed by Mirjalili of Griffith University [33]. This algorithm is inspired by the hunting behavior and cooperative behavior in the hunting process of gray wolves. The gray wolf family has a strict hierarchical management system. According to leadership, the wolf group can be divided into four levels: α, β, δ, and ω. After obtaining the position of the prey, the gray wolves will surround the prey. This process can be expressed by Equation (35): Among them, k is the number of iterations, X(k) is the position vector of the gray wolf after performing the kth iteration, X p (k) is the position vector of the prey after performing the kth iteration, A and C are coefficient vectors, and the expressions are shown in Equation (36): Among them, r 1 and r 2 are random numbers in the range of [0,1]. a is a variable that linearly decreases with the number of iterations, and its value gradually decreases from 2 to 0.
Through the above iterative process, the individual can be redirected to any position around the prey. However, this is not enough to reflect the group wisdom in the gray wolves. When gray wolves are hunting prey, the higher-level gray wolves play a key role in this process. After the encirclement process is completed, the gray wolves will hunt under the leadership of wolves with three levels of α, β, and δ. This process can be described by the following equations: When the process of encirclement and hunting has been completed by the gray wolves, they will attack the prey. This stage is manifested as the convergence of the solution in the process of solving the optimization problem.
In the offshore wind power generation system proposed in this paper, a set of parameters satisfying the definition are regarded as feasible solutions, and all combinations of parameters satisfying the definition constitute the solution space of the optimization problem. Using the gray wolf optimization algorithm to search in the solution space, the sub-optimal solution can be obtained in a limited time.
Usually, the basic method of power to gas from floating floating offshore wind power is using electrolysis of sea water to produce hydrogen. The principle of electrolyzing water to produce hydrogen is shown in Figure 10. It can be seen from the figure that the loss of electrons at the anode produces oxygen, and the electrons at the cathode produce hydrogen. The solution mixed in water increases the conductivity of the electrolyte. That is to say, the electrolysis of water is to immerse the cathode and anode in the electrolytic cell in water and apply direct current, respectively. The water will be decomposed and produce hydrogen at the cathode and oxygen at the anode. The main chemical reaction equations of the electrode reaction are as follows: In the water, The point of cathode, The point of anode, The chemical equation for the overall reaction is: The industrial experiment intends is to use a megawatt-class floating offshore wind turbine and hydrogen production equipment to verify the optimization effect of the control system. From the perspective of the system, concrete floating foundation, offshore wind turbine, hydrogen production equipment, the theory, method, and algorithm of this research can be verified step by step. At the same time, during the verification process, the floating offshore wind power sensing device is used to obtain the real-time status information of the intelligent control operation energy conversion online, and the controller data acquisition device is used to obtain the status information related to the wind power operation and the hydrogen production process. All relevant data are collected and aggregated through the distributed information network system. The non-grid-connected hydrogen production and logistics of floating offshore wind power to be used in industrial experiments are shown in Figure 11.

Simulation Results
To verify the effectiveness of the proposed controller, MATLAB/Simulink is used for simulation verification. The DFIG parameters used in the simulation are shown in Table 1. The values of the parameters in the table are selected with reference to real 2 MW wind turbines. The effect of the controller is evaluated by testing the response of the control system at random wind speeds. Considering that the maximum point tracking operation of the floating offshore wind turbine is between the cut-in wind speed and the rated wind speed, the wind speeds used in the simulation are in the interval from 3 m/s to 12 m/s. It is shown in Figure 12.  Figure 13 shows the variation curve of the rotor speed under the random wind speed when the proposed controller and PI controller are used. Figures 14-16 are the rotor speed in the intervals of 3 s, 7 s, 7 s, 11 s, and 11 s, 15 s, which are obtained by magnifying the three parts of Figure 13. Due to the mechanical inertia of the generator, there is a certain delay between the actual speed and the reference speed. Overall, the proposed controller has better dynamic response speed than the PI controller.     Figure 17 shows the corresponding reactive power curves of the two controllers. The reference value of reactive power is 0 W. The proposed controller can converge faster, but the PI controller will oscillate near the equilibrium position. Therefore, the proposed controller has better steady-state performance than the PI controller.  Figure 18 shows the simulation results of the power coefficient. The maximum power coefficient is 0.438 approximately. The proposed controller can adjust the rotor speed according to the change of wind speed, so the power coefficient of IM-ISMC controller is larger than PI controller in general. The purpose of the maximum power point tracking is to capture wind energy as much as possible, and the offshore wind energy power coefficient can be used as the basis for determining the amount of captured offshore wind energy. This means that floating offshore wind turbine systems can capture more offshore wind energy by using the proposed controller.  Figure 19 shows the different active power between the PI controller and the proposed controller. They are obtained by subtracting the active power value of the PI controller from the active power value of the proposed controller. At the beginning of the simulation, the power differential reaches its maximum value of 90 kW or so. In general, the proposed controller has higher active power than a PI controller at all times, especially in the time of the rotor increasing speed.

Conclusions
A current decoupling controller is presented for DFIG wind power systems on a floating offshore wind turbine to achieve MPPT. By using the internal model control strategy, an internal model open-loop controller is designed to simplify the parameter design. Moreover, the sliding mode control theory is used to achieve dynamic compensation based on the internal model open-loop controller, whereas the integral sliding mode controller is used to compensate the output error of the open-loop controller, and the parameters of proposed controller are designed through Gray Wolf Optimization. In addition, the performance of the controller is verified by testing the response of the controller under random wind speed in simulation software. Simulation results show that the proposed controller is superior to traditional PI controllers in terms of dynamic response and steadystate error. However, the proposed controller is not implemented on real floating offshore wind turbines. At the same time, how the controller parameters affect the control results has not been studied in depth. In future research, it is possible to build an experimental platform to conduct further research on the field of hydrogen production by floating offshore wind power.  Data Availability Statement: The study did not report any data.