An Adaptive Grid Voltage/Frequency Tracking Method Based on SOGIs on a Shipboard PV–Diesel-Battery Hybrid Power System

: This paper addresses the unbalanced voltage, subharmonic/dc-offset voltage, and low-frequency (LF)/high-frequency (HF) harmonics of a grid voltage tracking method based on second-order generalized integrators (SOGIs) in high voltage/frequency swing on a shipboard photovoltaic (PV)-diesel-battery hybrid power system. To perform this work, a kind of shipboard PV–diesel-battery hybrid power system structure was ﬁrst analyzed, emphasizing both the active and reactive power (PQ) control strategy and the sensitivity of the phase-locked loop (PLL) that is crucial to the vessel’s electrical networks. Then, the effect of grid voltage harmonics in SOGIs and of voltage/frequency swing on SOGI frequency-locked loop (SOGI-FLL) was studied. Meanwhile, aiming to the adverse power qualities of a shipboard power system (SPS), a SOGI-based structure with preﬁlter, a dc-offset block, and a positive sequence extractor (SOGI-FDE) was proposed. Finally, to overcome all of the vessel’s grid problems, a new SOGI-based voltage tracking structure, SOGI-FDE-FLL, consisting of SOGI-FDE and SOGI-FLL, was proposed to achieve accurate grid voltage tracking rapidly. This proposed schematic was used as an adaptive grid voltage tracking method to a three-phase three-wire shipboard PV–diesel-battery hybrid power system. Experimental results were obtained validating this proposal.


Introduction
Nowadays, diesel-engine synchronous generation is a primary shipboard power system (SPS). However, with the tendency of reducing the use of traditional fossil fuels and limiting the ship large-tonnage, which is increasing quickly with the size of ships and their associated electrical power equipment, the development of new energy supplies instead of diesel is an urgent demand with regard to the quick development of shipboard electric power systems [1,2]. Photovoltaic (PV) generation utilizing solar energy is one of the best choices for power supply on ships due to its excellent economy and high efficiency. However, the features of PV such as intermittency and limited development area, has limited solar-powered ships working within a short navigation range with rather high costs [3]. Hence, the PV-diesel-battery hybrid power system is a promising power supply system in comparison with traditional power generation systems and can enhance fuel economy and the continuity of power supply, reduce marine pollution, and improve grid power quality, system stability, and flexibility [4][5][6].
Onboard AC electrical power system may be categorized as an islanded AC microgrid [7], which causes numerous power quality issues. One prominent problem is the voltage/frequency swing during

Power System Description
The adaptive grid voltage tracking method described here for applying a photovoltaic battery to SPS revealed superior performance in adverse grid conditions. A shipboard PV-diesel-battery hybrid power system was chosen to introduce the proposed method in a practical environment. Figure 1 depicts the configuration of a shipboard PV-diesel-battery hybrid power system, which was composed of three kinds of power supply (a PV unit, a diesel-engine synchronous generator, and a battery unit), onboard electrical loads, and a PV-battery inverter. The PV unit included a group of PV cells and a dc-dc converter with one-way power flow for maximum power point tracking (MPPT) and dc output voltage stability. The battery unit consisted of a quantity of lead-acid or lithium battery cells, with parallel and series connection, and one bidirectional dc-dc converter as a charging or discharging circuit. The power flow of this battery converter was provided by a power management system (PMS) based on the battery's state-of-charge (SOC), the surrounding temperature, the load power, etc. The above two dc-dc converters connected the inverter to the AC bus to which the diesel-engine synchronous generator was connected. These three units provide electricity together to the onboard electrical loads, which differ from a terrestrial power-electronics-based distributed power generation system (DPGS). For the limitation of deck area, the electric power capacity of the PV unit, compared to the diesel generator unit, was lower. Meanwhile, considering the vessel's safety on the ocean as a moving platform, the PV unit equipment must satisfy the rules and requirements of navigation safety.

Shipboard PV-Diesel-Battery Hybrid Power System Configuration
In this case, the diesel-engine synchronous generator worked in standalone mode as a slack bus where the diesel governor and synchronous excitation system maintained the grid voltage frequency and amplitude of the shipboard PV-diesel-battery hybrid power system. The PV-battery inverter acts as a PQ bus to inject the desired active and reactive powers given by the PMS to the grid.
Energies 2018, 11, x FOR PEER REVIEW 4 of 20 and amplitude of the shipboard PV-diesel-battery hybrid power system. The PV-battery inverter acts as a PQ bus to inject the desired active and reactive powers given by the PMS to the grid.

Figure 1.
Schematic of the shipboard PV-diesel-battery hybrid power system. In this power system structure, the PV-battery inverter is the key role for the interface of PV and battery to the grid, which was composed of a three-phase full bridge circuit and inductor and capacitor (LC) output filter. The extra line inductance Lg connecting the PV-battery inverter and AC grid reduced the active and reactive power flows coupling. The active and reactive power outputs of the inverter are given by (2) where P and Q are the active and reactive power outputs;  is the angle between the diesel-engine generator and the PV-battery inverter; X is the line impedance; Vg and Vo are the output voltages of the diesel-engine generator and PV-battery inverter, respectively. P and Q can be controlled by V and  from Equations (1) and (2).

PV-Battery Inverter Control Strategy
The control part of the PV-battery inverter consisted of inner voltage and current control loops, and outer P and Q droop control loops, assigning power by mimicking the static droop feature of a synchronous generator. The voltage and current control loops were to maintain the frequency and amplitude of the output voltage to the reference, which was regulated by the P and Q droop loops. Furthermore, the P and Q droop loops may adjust the R/X ratio of the output impedance "seen" by the inverter, which was composed of the power calculation, the power droop control, and the reference voltage generator. The inverter active power P and reactive power Q was obtained by calculating the inverter output current and voltage. The difference in active and reactive power (ΔP and ΔQ) was found via comparison with the corresponding reference active power P* and reactive power Q* from the PMS by the given optimization algorithm. Then, the angular frequency and voltage differences (Δω and ΔV) were obtained as the active power and reactive power PI controller outputs with ΔP and ΔQ as the relevant inputs. The inverter output voltage angular frequency reference ω* and amplitude reference V* were derived by tracking the output voltage angular frequency ωg and amplitude Vg by the phase-locked loop (PLL) and minus Δω and ΔV. ω* and V* had droop characteristics, which can be given by  Figure 1. Schematic of the shipboard PV-diesel-battery hybrid power system. In this power system structure, the PV-battery inverter is the key role for the interface of PV and battery to the grid, which was composed of a three-phase full bridge circuit and inductor and capacitor (LC) output filter. The extra line inductance L g connecting the PV-battery inverter and AC grid reduced the active and reactive power flows coupling. The active and reactive power outputs of the inverter are given by where P and Q are the active and reactive power outputs; ϕ is the angle between the diesel-engine generator and the PV-battery inverter; X is the line impedance; V g and V o are the output voltages of the diesel-engine generator and PV-battery inverter, respectively. P and Q can be controlled by V and ϕ from Equations (1) and (2).

PV-Battery Inverter Control Strategy
The control part of the PV-battery inverter consisted of inner voltage and current control loops, and outer P and Q droop control loops, assigning power by mimicking the static droop feature of a synchronous generator. The voltage and current control loops were to maintain the frequency and amplitude of the output voltage to the reference, which was regulated by the P and Q droop loops. Furthermore, the P and Q droop loops may adjust the R/X ratio of the output impedance "seen" by the inverter, which was composed of the power calculation, the power droop control, and the reference voltage generator. The inverter active power P and reactive power Q was obtained by calculating the inverter output current and voltage. The difference in active and reactive power (∆P and ∆Q) was found via comparison with the corresponding reference active power P* and reactive power Q* from the PMS by the given optimization algorithm. Then, the angular frequency and voltage differences (∆ω and ∆V) were obtained as the active power and reactive power PI controller outputs with ∆P and ∆Q as the relevant inputs. The inverter output voltage angular frequency reference ω* and amplitude reference V* were derived by tracking the output voltage angular frequency ω g and amplitude V g by the phase-locked loop (PLL) and minus ∆ω and ∆V. ω* and V* had droop characteristics, which can be given by where k pp and k ip are the proportional and integral coefficients of active power PI controller; k pq and k iq are the proportional and integral coefficients of the reactive power PI controller, respectively. Both the active and reactive power PI controllers are essentially adaptive droop coefficients.
The inner voltage and current control loops adopted the αβ stationary reference frame and PR controller for tracking the sinusoidal voltage reference with zero static error. The inverter's output power can follow the power reference given by PMS by space vector pulse width modulation (SVPWM) and the power switches control.
In all, for this shipboard PV-diesel-battery hybrid power system with finite power capacity, the PLL is one of the most important units for the PV-battery inverter as a current source. However, for special and serious grid problems regarding the power quality of broad voltage/frequency swing in heavy load transients, harmonics, unbalanced voltage, and subharmonic/dc-offset voltage in the SPS, the PLL's accuracy and response directly influences the overall SPS's dynamic characteristic and stability. In this paper, the SOGI algorithm was chosen as the basic unit for grid voltage tracking. In the next section, the SOGI structure and its limitations are analyzed in detail. Figure 2a depicts the block diagram of a SOGI structure. The SOGI is a frequency-adjustable resonator that is damped by an output loop. This outer loop introduces the error e v between the resonator output v α and the input signal v in the resonator input through gain k.

SOGI
where kpp and kip are the proportional and integral coefficients of active power PI controller; kpq and kiq are the proportional and integral coefficients of the reactive power PI controller, respectively. Both the active and reactive power PI controllers are essentially adaptive droop coefficients.
The inner voltage and current control loops adopted the αβ stationary reference frame and PR controller for tracking the sinusoidal voltage reference with zero static error. The inverter's output power can follow the power reference given by PMS by space vector pulse width modulation (SVPWM) and the power switches control.
In all, for this shipboard PV-diesel-battery hybrid power system with finite power capacity, the PLL is one of the most important units for the PV-battery inverter as a current source. However, for special and serious grid problems regarding the power quality of broad voltage/frequency swing in heavy load transients, harmonics, unbalanced voltage, and subharmonic/dc-offset voltage in the SPS, the PLL's accuracy and response directly influences the overall SPS's dynamic characteristic and stability. In this paper, the SOGI algorithm was chosen as the basic unit for grid voltage tracking. In the next section, the SOGI structure and its limitations are analyzed in detail. Figure 2a depicts the block diagram of a SOGI structure. The SOGI is a frequency-adjustable resonator that is damped by an output loop. This outer loop introduces the error ev between the resonator output vα and the input signal v in the resonator input through gain k. The SOGI outputs are vα and vβ, which are, respectively, the in-phase and quadrature-phase with the input signal v. From the input-to-output point of view, the SOGI behaves as a narrow bandpass filter (BPF) for vα to be a specific frequency selector, and as a low-pass filter (LPF) for vβ to be a highfrequency harmonics suppressor. The corresponding transfer functions are

SOGI
where ωg is oscillation frequency; and k is a damping constant of the SOGI system which is in direct proportion to bandwidth and can be determined by the scope frequency of input voltage. The SOGI The SOGI outputs are v α and v β , which are, respectively, the in-phase and quadrature-phase with the input signal v. From the input-to-output point of view, the SOGI behaves as a narrow bandpass filter (BPF) for v α to be a specific frequency selector, and as a low-pass filter (LPF) for v β to be a high-frequency harmonics suppressor. The corresponding transfer functions are where ω g is oscillation frequency; and k is a damping constant of the SOGI system which is in direct proportion to bandwidth and can be determined by the scope frequency of input voltage. The SOGI is kept as an underdamping system for a better transient response, given k < 2. The bode magnitude plots of Equations (5) and (6) are depicted in Figure 2b for k = 1.414. According to Equations (5) and (6), given that the error band is 5%, the key parameters related to the SOGI time response and frequency selectivity are the settling time: and the SOGI BPF gain at a given harmonic h order of ω g where 0 < k/2 < 1 and h is an odd integer number, and h >> 1. Taking into account the typical h values 3, 5, 7, and so on, Equation (8) can be further simplified as From Equations (7)-(9), it can be seen that the SOGI has a tradeoff relationship between harmonic attenuation and settling time. The value of settling time should be decreased when the SOGI is used on the occasion of serious voltage/frequency variation.
However, when the frequency of the grid voltage swing as well as the in-phase and quadrature-phase voltage frequencies are unstable for the oscillation frequency ω g coming from outside, an adaptive frequency-locked loop (FLL) unit is necessary.

SOGI-FLL
Based on the SOGI structure shown in Figure 2, the error transfer function E(s) from input signal v to error e can be deduced: Figure 3a depicts the bode magnitude and phase plots of Q(s) and E(s). It can be seen that the relationship between the input signal frequency and ω g can be obtained by the sign of v β multiplied by v. The SOGI-FLL scheme is shown in Figure 3b, where ω o is preset to the grid rated frequency, and γ is the proportional rate of the frequency differential coefficient; voltage normalization can eliminate the impact of voltage amplitude fluctuation on the SOGI's in-phase and quadrature-phase voltage. is kept as an underdamping system for a better transient response, given k < 2. The bode magnitude plots of Equations (5) and (6) are depicted in Figure 2b for k = 1.414. According to Equations (5) and (6), given that the error band is 5%, the key parameters related to the SOGI time response and frequency selectivity are the settling time: and the SOGI BPF gain at a given harmonic h order of ωg where 0 < k/2 < 1 and h is an odd integer number, and h >> 1. Taking into account the typical h values 3, 5, 7, and so on, Equation (8) can be further simplified as From Equations (7)-(9), it can be seen that the SOGI has a tradeoff relationship between harmonic attenuation and settling time. The value of settling time should be decreased when the SOGI is used on the occasion of serious voltage/frequency variation.
However, when the frequency of the grid voltage swing as well as the in-phase and quadraturephase voltage frequencies are unstable for the oscillation frequency ωg coming from outside, an adaptive frequency-locked loop (FLL) unit is necessary.

SOGI-FLL
Based on the SOGI structure shown in Figure 2, the error transfer function E(s) from input signal v to error e can be deduced: It can be seen that the relationship between the input signal frequency and ωg can be obtained by the sign of vβ multiplied by v. The SOGI-FLL scheme is shown in Figure 3b, where ωo is preset to the grid rated frequency, and γ is the proportional rate of the frequency differential coefficient; voltage normalization can eliminate the impact of voltage amplitude fluctuation on the SOGI's in-phase and quadrature-phase voltage.
where φ = ∠D(jω), and φ ∈ [−π/2 ∼ π/2]. Hence, the steady v α and v β can be deduced as where V α and V β are the corresponding amplitude of v α and v β . e is defined as the frequency error, which is where ∆ω is the deviation of ω and ω g , namely The corresponding average value of e in the fundamental frequency can be implied which is a one-order integrator frequency without a steady-state error. Figure 4 depicts the steady equivalent frequency closed loop system. ω g can be regulated and go back to ω ultimately for the integrator term. The system transfer function can be expressed as where T is the time constant, and T = k · ω o /γ. The rise time t r of the frequency in step response is where Vα and Vβ are the corresponding amplitude of vα and vβ. e is defined as the frequency error, which is v β where Δω is the deviation of ω and ωg, namely The corresponding average value of e in the fundamental frequency can be implied e k     (17) which is a one-order integrator frequency without a steady-state error. Figure 4 depicts the steady equivalent frequency closed loop system. ωg can be regulated and go back to ω ultimately for the integrator term. The system transfer function can be expressed as 1 s1 where T is the time constant, and   o Tk . The rise time tr of the frequency in step response is

Limitations of SOGIs
Quickly obtaining the accuracy grid voltage amplitude and frequency in the PV-diesel-battery hybrid power system is key for the serious voltage/frequency swings in the process of ship load mutation. The SOGI has the characteristics of sinusoidal signal tracking for second-order resonant and active filtering for HF attenuation as well as adaptive lock frequency when connecting the FLL unit. However, a SOGI will run into problems when it is applied to an SPS for a PV-battery inverter PLL.

Limitations of SOGIs
Quickly obtaining the accuracy grid voltage amplitude and frequency in the PV-diesel-battery hybrid power system is key for the serious voltage/frequency swings in the process of ship load mutation. The SOGI has the characteristics of sinusoidal signal tracking for second-order resonant and active filtering for HF attenuation as well as adaptive lock frequency when connecting the FLL unit. However, a SOGI will run into problems when it is applied to an SPS for a PV-battery inverter PLL.

. Unbalanced Voltage Impact on SOGI Input Signal
For unequal three-phase heavy loads connected to the grid easily arousing unbalanced grid voltage in an SPS, the unbalanced grid voltage v gabc can be expressed as (20) where v gabc+ , v gabc− , and v gabc0 represent its positive sequence, negative sequence, and zero sequence components, respectively. Accurate grid voltage information cannot be obtained by SOGI on unbalanced voltage due to the unequal three-phase voltage amplitude. Hence, for SPS application, a positive sequence extractor (PSE) unit should be included at the front of the SOGI.

Subharmonic and dc-Offset Impact on SOGI Quadrature-Phase Output Signal
Assuming there are dc-offset components in the input signal v of the SOGI, namely where F is the ratio of subharmonic amplitude to fundamental voltage amplitude; ω L is the subharmonic angle frequency and, when ω L is lower or zero, namely ω L << ω g , and φ L is close to π/2, can be expressed as From Equations (11)- (14), v α and v β are approximately simplified to Figure 5 shows the SOGI response to a subharmonic and a dc-bias in k = 1.414. The subharmonic was 5 Hz at 20% of V, and the dc bias was equal to 20% of V. As seen in this figure, this subharmonic made v β produce a double sinusoidal envelope with a 5 Hz and a 0.56 amplitude difference; the dc bias brought a 0.28 dc bias to v β . The presence of subharmonics and dc bias directly translated to the output v β with the gain of kF and had a direct impact on the SOGI quadrature-phase output, but had almost no impact on in-phase output, which can produce problems and distortion in all SOGI-based structures including SOGI-FLL. This impact can also be illuminated by the quadrature-phase behavior of an LPF with gain k, as seen in Figure 2b.

. Unbalanced Voltage Impact on SOGI Input Signal
For unequal three-phase heavy loads connected to the grid easily arousing unbalanced grid voltage in an SPS, the unbalanced grid voltage vgabc can be expressed as where vgabc+, vgabc−, and vgabc0 represent its positive sequence, negative sequence, and zero sequence components, respectively. Accurate grid voltage information cannot be obtained by SOGI on unbalanced voltage due to the unequal three-phase voltage amplitude. Hence, for SPS application, a positive sequence extractor (PSE) unit should be included at the front of the SOGI.

Subharmonic and dc-Offset Impact on SOGI Quadrature-Phase Output Signal
Assuming there are dc-offset components in the input signal v of the SOGI, namely cos cos where F is the ratio of subharmonic amplitude to fundamental voltage amplitude; ωL is the subharmonic angle frequency and, when ωL is lower or zero, namely ωL << ωg, and φL is close to π/2, can be expressed as 22 arctan From Equations (11)-(14), vα and vβ are approximately simplified to Figure 5 shows the SOGI response to a subharmonic and a dc-bias in k = 1.414. The subharmonic was 5 Hz at 20% of V, and the dc bias was equal to 20% of V. As seen in this figure, this subharmonic made vβ produce a double sinusoidal envelope with a 5 Hz and a 0.56 amplitude difference; the dc bias brought a 0.28 dc bias to vβ. The presence of subharmonics and dc bias directly translated to the output vβ with the gain of kF and had a direct impact on the SOGI quadrature-phase output, but had almost no impact on in-phase output, which can produce problems and distortion in all SOGI-based structures including SOGI-FLL. This impact can also be illuminated by the quadrature-phase behavior of an LPF with gain k, as seen in Figure 2b.  Moreover, when the SOGI-FLL is applied to SPS (as shown in Figure 1 for PLL), the existing subharmonics and dc bias may arouse the oscillation of voltage in amplitude and frequency, which makes the voltage reference of the voltage and current close loops unstable, which then affects the stability of the whole shipboard PV-diesel-battery power system. This kind of oscillation is more serious with the dc offset rise. Hence, SOGI-FLL cannot be used directly, and a low-frequency attenuation unit is needed.

Low-Order Harmonic Impact on SOGI In-Phase Output Signal
The presence of low-order harmonics in the input signal has a direct impact on the SOGI in-phase output since this output behaves as a bandpass filter (BPF) with −20 decibels per decade slope at frequencies greater than the oscillation frequency (see Figure 2b). This attenuation rate obtains a better filtering performance for switching harmonics but is not sufficient for low-order harmonics. To evaluate this problem, Figure 6 shows the SOGI in-phase response to 5th and 7th harmonics that are both at 15% of the fundamental voltage. The 5th and 7th harmonics were partly translated to the output v α and aroused a higher (5.21%) total harmonics distortion (THD). When SOGI was directly used in the proposed SPS shown in Figure 1, the in-phase output signals with low-frequency harmonics could be directly delivered to the grid through the PV-battery inverter making grid-injected current distortion, which then decreased the power qualities. Hence, a sharp slope of attenuation to harmonics is necessary to deal with the low-order harmonics existing in SPS. Moreover, when the SOGI-FLL is applied to SPS (as shown in Figure 1 for PLL), the existing subharmonics and dc bias may arouse the oscillation of voltage in amplitude and frequency, which makes the voltage reference of the voltage and current close loops unstable, which then affects the stability of the whole shipboard PV-diesel-battery power system. This kind of oscillation is more serious with the dc offset rise. Hence, SOGI-FLL cannot be used directly, and a low-frequency attenuation unit is needed.

Low-Order Harmonic Impact on SOGI In-Phase Output Signal
The presence of low-order harmonics in the input signal has a direct impact on the SOGI inphase output since this output behaves as a bandpass filter (BPF) with −20 decibels per decade slope at frequencies greater than the oscillation frequency (see Figure 2b). This attenuation rate obtains a better filtering performance for switching harmonics but is not sufficient for low-order harmonics. To evaluate this problem, Figure 6 shows the SOGI in-phase response to 5th and 7th harmonics that are both at 15% of the fundamental voltage. The 5th and 7th harmonics were partly translated to the output vα and aroused a higher (5.21%) total harmonics distortion (THD). When SOGI was directly used in the proposed SPS shown in Figure 1, the in-phase output signals with low-frequency harmonics could be directly delivered to the grid through the PV-battery inverter making gridinjected current distortion, which then decreased the power qualities. Hence, a sharp slope of attenuation to harmonics is necessary to deal with the low-order harmonics existing in SPS. Figure 6. The SOGI vα/V and vβ/V response with the 5th and 7th harmonics both with 15% voltage amplitude, and the corresponding THD is 5.21% and 0.95%.

The Proposed Adaptive Voltage/Frequency Tracking Method Based on SOGIs
A SOGI was considered in this study for its BPF for in-phase output and its LPF for quadraturephase output. However, when applied to an SPS for voltage tracking, the primary issues were the unbalanced voltage impact on the input signal, the subharmonic and dc-offset impact on the quadrature-phase output signal, and the low-order harmonic impact on the SOGI in-phase output signal.

The Proposed SOGI-FDE Configuration
The three-phase three-wire system is adopted in most cases for high reliability on the shipboard power supply. Once the three-phase voltage is unbalanced, there are positive and negative sequence voltage terms as shown in Equation (20). Setting vgαβ as a three-phase grid voltage in the αβ-frame, the corresponding positive sequence term  where Tαβ is the coordinate transformation (CT) from the abc-frame to the αβ-frame, namely,

The Proposed Adaptive Voltage/Frequency Tracking Method Based on SOGIs
A SOGI was considered in this study for its BPF for in-phase output and its LPF for quadrature-phase output. However, when applied to an SPS for voltage tracking, the primary issues were the unbalanced voltage impact on the input signal, the subharmonic and dc-offset impact on the quadrature-phase output signal, and the low-order harmonic impact on the SOGI in-phase output signal.

The Proposed SOGI-FDE Configuration
The three-phase three-wire system is adopted in most cases for high reliability on the shipboard power supply. Once the three-phase voltage is unbalanced, there are positive and negative sequence voltage terms as shown in Equation (20). Setting v gαβ as a three-phase grid voltage in the αβ-frame, the corresponding positive sequence term v + αβ0 can then be expressed as where T αβ is the coordinate transformation (CT) from the abc-frame to the αβ-frame, namely, In this paper, to simplify the schematic of the PLL under adverse grid conditions in the SPS, a synthetic SOGI-based structure with prefilter, a dc-offset block, and a positive sequence extractor (SOGI-FDE) as depicted in Figure 7, is proposed. This structure is composed of a coordinate transformation (CT) T αβ for the grid voltage abc/αβ-frame transformation and three SOGI units. SOGI I and II, connected to the grid voltage in αand β-frame components as the in-phase outputs, constitute two active filters (AF) and a subharmonics/dc-bias eliminator for grid voltage in the stationary α-axis and β-axis components, respectively, for LF attenuation. Meanwhile, SOGI III, linked to SOGI II, positively rotates the β-axis component v gβ ' π/2 for the imaginary unit j generator of PSE in Equation (25). The PSE unit extracts the positive sequence α-axis component from the unbalanced three-phase grid voltage.
Energies 2018, 11, x FOR PEER REVIEW 10 of 20 αβ 2 3 1 3 1 3 In this paper, to simplify the schematic of the PLL under adverse grid conditions in the SPS, a synthetic SOGI-based structure with prefilter, a dc-offset block, and a positive sequence extractor (SOGI-FDE) as depicted in Figure 7, is proposed. This structure is composed of a coordinate transformation (CT) Tαβ for the grid voltage abc/αβ-frame transformation and three SOGI units. SOGI I and II, connected to the grid voltage in α-and β-frame components as the in-phase outputs, constitute two active filters (AF) and a subharmonics/dc-bias eliminator for grid voltage in the stationary α-axis and β-axis components, respectively, for LF attenuation. Meanwhile, SOGI III, linked to SOGI II, positively rotates the β-axis component vgβ' π/2 for the imaginary unit j generator of PSE in Equation (25). The PSE unit extracts the positive sequence α-axis component from the unbalanced three-phase grid voltage.

Proposed SOGI-FDE-FLL for Vessel's SPS in Adverse Grid Conditions
For the mentioned shipboard PV-diesel-battery hybrid power system in Section 2, the PQ control method was adopted for the PV-battery inverter due to the lower capacity of the photovoltaic battery when compared to a diesel engine on the vessel. The grid-voltage and frequency information are needed to generate a voltage reference for the voltage and current closed-loop control. However, for the vessel's adverse grid condition of serious voltage/frequency swing, harmonics, unbalanced voltage, subharmonics, and dc bias component, an adaptive grid voltage/frequency tracking method based on SOGIs to a PV-diesel-battery hybrid SPS was proposed, as depicted in Figure 8. This schematic is composed of SOGI-FDE and SOGI-FLL, named SOGI-FDE-FLL.
In this scheme, the FLL is used only in the α-frame to make the sequence detection adaptive to the grid frequency. The SOGI-FDE unit, for a pre-filter, a subharmonic and dc-offset eliminator and a positive sequence extractor, is connected directly to the vessel's grid to obtain the positive sequence fundamental voltage vα+ without the subharmonics and dc bias component. From the SOGIs' characteristics of Equations (5) and (6) and the abc-frame to αβ-frame transformation of Equations (25) and (26), the voltage amplitude, frequency, and phase of vα+ is the a-axis positive sequence voltage in the abc-frame. This was the key opinion of this paper to track vα+ to obtain accurate voltage and frequency estimated values instead of tracking grid voltage vgabc with more power quality problems directly.
The SOGI-FLL unit, a grid voltage/frequency estimator, is used to quickly obtain accurate information of grid voltage and frequency in their serious fluctuations. Vg, the estimated voltage amplitude, is calculated as the root square of its in-phase and quadrature-phase output, whereas the estimated frequency ωg is given by an adaptive frequency schematic of an equivalent frequency closed-loop control. Vg and ωg subtract the differences of voltage and frequency obtained from Q-V and P-f droop relationships, respectively, and then reference the voltage generator to form the threephase reference signal of the voltage/current closed-loop control system as shown in Figure 1.

Proposed SOGI-FDE-FLL for Vessel's SPS in Adverse Grid Conditions
For the mentioned shipboard PV-diesel-battery hybrid power system in Section 2, the PQ control method was adopted for the PV-battery inverter due to the lower capacity of the photovoltaic battery when compared to a diesel engine on the vessel. The grid-voltage and frequency information are needed to generate a voltage reference for the voltage and current closed-loop control. However, for the vessel's adverse grid condition of serious voltage/frequency swing, harmonics, unbalanced voltage, subharmonics, and dc bias component, an adaptive grid voltage/frequency tracking method based on SOGIs to a PV-diesel-battery hybrid SPS was proposed, as depicted in Figure 8. This schematic is composed of SOGI-FDE and SOGI-FLL, named SOGI-FDE-FLL.
In this scheme, the FLL is used only in the α-frame to make the sequence detection adaptive to the grid frequency. The SOGI-FDE unit, for a pre-filter, a subharmonic and dc-offset eliminator and a positive sequence extractor, is connected directly to the vessel's grid to obtain the positive sequence fundamental voltage v α+ without the subharmonics and dc bias component. From the SOGIs' characteristics of Equations (5) and (6) and the abc-frame to αβ-frame transformation of Equations (25) and (26), the voltage amplitude, frequency, and phase of v α+ is the a-axis positive sequence voltage in the abc-frame. This was the key opinion of this paper to track v α+ to obtain accurate voltage and frequency estimated values instead of tracking grid voltage v gabc with more power quality problems directly.
The SOGI-FLL unit, a grid voltage/frequency estimator, is used to quickly obtain accurate information of grid voltage and frequency in their serious fluctuations. V g , the estimated voltage amplitude, is calculated as the root square of its in-phase and quadrature-phase output, whereas the estimated frequency ω g is given by an adaptive frequency schematic of an equivalent frequency closed-loop control. V g and ω g subtract the differences of voltage and frequency obtained from Q-V and P-f droop relationships, respectively, and then reference the voltage generator to form the three-phase reference signal of the voltage/current closed-loop control system as shown in Figure 1. In the above proposed schematic, the SOGI-FDE unit is in front of the SOGI-FLL. vα+, the input signal of SOGI-FLL, is associated with the vgα and jvgβ transmitted from the grid α-frame component vgα0 by SOGI I and from the β-frame component vgβ0 by SOGIs II and III for PSE. Meanwhile, it is transmitted through SOGI IV to the in-phase output vα. Hence, the schematic can be considered as a double-input single-output system. If the damping constant of all the SOGIs are equal, the Laplace transfer function G1(s), the ratio of vg to vgα0, and G2(s), the ratio of vg to vgβ0, can be expressed as where D(s) and Q(s) are the SOGI's transfer functions of corresponding in-phase and quadraturephase output to input signals (see Equations (5) and (6)). Figure 9 depicts the D(s), G1(s), and G2(s) magnitude plot in LF harmonics with the same damping constants of k = 1.414 and oscillation frequency of ωg = 314 rad/s. All of these have the capability of dc bias attenuation. D(s) represents the 1st-order SOGI of in-phase output with a narrow BPF, but insufficient LF harmonics suppression at 20 dB per decade slope. G1(s) was obtained by the 2nd-order SOGI, where the LF harmonics obtain substantial attenuation at 40 dB per decade slope for an s-domain of 2 zeros and 2 pairs of left-half poles. The 5th harmonic was reduced to 3%, and 0.1-order subharmonics can be degraded to lower than 3% (see Figure 9 for illustration). G2(s) was composed of the 3rd-order SOGI with the same G1(s) in the subharmonic zone and a sharp attenuation to the low-order harmonics at 80 dB per decade slope decrease. Hence, besides the positive sequence extractor from the unbalanced voltage, this proposed schematic has good features of LF harmonics, sharply fading out, including subharmonics and low-order harmonics without a dc bias, which allows this schematic to overcome the limitations of the SOGI's limitations -mentioned in Section 3.3.

SOGI II
0.5 In the above proposed schematic, the SOGI-FDE unit is in front of the SOGI-FLL. v α+ , the input signal of SOGI-FLL, is associated with the v gα and jv gβ transmitted from the grid α-frame component v gα0 by SOGI I and from the β-frame component v gβ0 by SOGIs II and III for PSE. Meanwhile, it is transmitted through SOGI IV to the in-phase output v α . Hence, the schematic can be considered as a double-input single-output system. If the damping constant of all the SOGIs are equal, the Laplace transfer function G 1 (s), the ratio of v g to v gα0 , and G 2 (s), the ratio of v g to v gβ0 , can be expressed as

SOGI IV
where D(s) and Q(s) are the SOGI's transfer functions of corresponding in-phase and quadrature-phase output to input signals (see Equations (5) and (6)). Figure 9 depicts the D(s), G 1 (s), and G 2 (s) magnitude plot in LF harmonics with the same damping constants of k = 1.414 and oscillation frequency of ω g = 314 rad/s. All of these have the capability of dc bias attenuation. D(s) represents the 1st-order SOGI of in-phase output with a narrow BPF, but insufficient LF harmonics suppression at 20 dB per decade slope. G 1 (s) was obtained by the 2nd-order SOGI, where the LF harmonics obtain substantial attenuation at 40 dB per decade slope for an s-domain of 2 zeros and 2 pairs of left-half poles. The 5th harmonic was reduced to 3%, and 0.1-order subharmonics can be degraded to lower than 3% (see Figure 9 for illustration). G 2 (s) was composed of the 3rd-order SOGI with the same G 1 (s) in the subharmonic zone and a sharp attenuation to the low-order harmonics at 80 dB per decade slope decrease. Hence, besides the positive sequence extractor from the unbalanced voltage, this proposed schematic has good features of LF harmonics, sharply fading out, including subharmonics and low-order harmonics without a dc bias, which allows this schematic to overcome the limitations of the SOGI's limitations -mentioned in Section 3.3. Energies 2018, 11, x FOR PEER REVIEW 12 of 20 In all, this proposed voltage/frequency tracking method based on SOGIs has good features regarding subharmonics and dc bias rejection, as well as voltage/frequency capture in grid voltage serious oscillation. This method is thus suitable for a shipboard power system with new energy generation.

The Parameters Optimization Design of SOGI-FDE-FLL
There were four damping constants in four SOGIs and one error gain in FLL in SOGI-FDE-FLL. The values of those parameters have close relations to the characteristics of the structure. Figure 10 depicts the relationships of the ratios of in-phase voltage amplitude Vα and quadrature-phase voltage amplitude Vβ to input voltage amplitude, respectively, with their frequency ratio ω/ωg from 0.8 to 1.2 and damping constant k from 0.5 to 2, obtained from Equations (11)- (14) to emulate the grid voltage/frequency swing. The values of k have great effects on both the in-phase and quadrature-phase outputs, especially on the former when the input signal amplitude and frequency fluctuate. When ω is equal to the oscillation frequency ωg, Vα and Vβ are both equal to V, and at that time, from Equations (13) and (14), V can be deduced from

The Characteristics of Damping Constant in SOGI
which means ωg is kept equal to ω, so the estimated amplitude can be obtained from SOGI's in-phase and quadrature-phase signals. Moreover, from the input-phase output-voltage curve shown in Figure 10a, it can be seen that In all, this proposed voltage/frequency tracking method based on SOGIs has good features regarding subharmonics and dc bias rejection, as well as voltage/frequency capture in grid voltage serious oscillation. This method is thus suitable for a shipboard power system with new energy generation.

The Parameters Optimization Design of SOGI-FDE-FLL
There were four damping constants in four SOGIs and one error gain in FLL in SOGI-FDE-FLL. The values of those parameters have close relations to the characteristics of the structure. Figure 10 depicts the relationships of the ratios of in-phase voltage amplitude V α and quadrature-phase voltage amplitude V β to input voltage amplitude, respectively, with their frequency ratio ω/ω g from 0.8 to 1.2 and damping constant k from 0.5 to 2, obtained from Equations (11)- (14) to emulate the grid voltage/frequency swing. The values of k have great effects on both the in-phase and quadrature-phase outputs, especially on the former when the input signal amplitude and frequency fluctuate. When ω is equal to the oscillation frequency ω g , V α and V β are both equal to V, and at that time, from Equations (13) and (14), V can be deduced from

The Characteristics of Damping Constant in SOGI
which means ω g is kept equal to ω, so the estimated amplitude can be obtained from SOGI's in-phase and quadrature-phase signals. In all, this proposed voltage/frequency tracking method based on SOGIs has good features regarding subharmonics and dc bias rejection, as well as voltage/frequency capture in grid voltage serious oscillation. This method is thus suitable for a shipboard power system with new energy generation.

The Parameters Optimization Design of SOGI-FDE-FLL
There were four damping constants in four SOGIs and one error gain in FLL in SOGI-FDE-FLL. The values of those parameters have close relations to the characteristics of the structure. Figure 10 depicts the relationships of the ratios of in-phase voltage amplitude Vα and quadrature-phase voltage amplitude Vβ to input voltage amplitude, respectively, with their frequency ratio ω/ωg from 0.8 to 1.2 and damping constant k from 0.5 to 2, obtained from Equations (11)- (14) to emulate the grid voltage/frequency swing. The values of k have great effects on both the in-phase and quadrature-phase outputs, especially on the former when the input signal amplitude and frequency fluctuate. When ω is equal to the oscillation frequency ωg, Vα and Vβ are both equal to V, and at that time, from Equations (13) and (14), V can be deduced from

The Characteristics of Damping Constant in SOGI
which means ωg is kept equal to ω, so the estimated amplitude can be obtained from SOGI's in-phase and quadrature-phase signals. Moreover, from the input-phase output-voltage curve shown in Figure 10a, it can be seen that Moreover, from the input-phase output-voltage curve shown in Figure 10a, it can be seen that and V α is increasingly smaller with decreasing k. The higher k is, the smaller the value of V minus V α is. However, the quadrature-phase output voltage curve shown in Figure 10b had variation trends different from that shown in Figure 10a. When the input signal fluctuating with k was greater than 1, V β /V and ω / ω g had approximately linear relations and V β decreased as ω increased. The higher k is, the more this linearity improves. The damping constant can be determined by the voltage variance from Figure 9, and was replotted to 2D as shown in Figure 11. Compared to V β , V α had less voltage variance. Increasing k is of benefit for decreasing this variance. On the other hand, V β , as shown in Figure 10b, has relatively less variance when k is from 1.0 to 1.5 but tends to have proportional variance to the continuously increasing k. This linear relationship improves the stability in the transient response and can be expressed as Energies 2018, 11, x FOR PEER REVIEW 13 of 20 and Vα is increasingly smaller with decreasing k. The higher k is, the smaller the value of V minus Vα is. However, the quadrature-phase output voltage curve shown in Figure 10b had variation trends different from that shown in Figure 10a. When the input signal fluctuating with k was greater than 1, Vβ/V and ω/ωg had approximately linear relations and Vβ decreased as ω increased. The higher k is, the more this linearity improves. The damping constant can be determined by the voltage variance from Figure 9, and was replotted to 2D as shown in Figure 11. Compared to Vβ, Vα had less voltage variance. Increasing k is of benefit for decreasing this variance. On the other hand, Vβ, as shown in Figure 10b, has relatively less variance when k is from 1.0 to 1.5 but tends to have proportional variance to the continuously increasing k. This linear relationship improves the stability in the transient response and can be expressed as Above all, the frequency variation of the input signal is the key point of amplitude variance generation in the in-phase and quadrature-phase outputs, and the damping constant is the key parameter for the SOGI unit, which directly affects the two output signals. Increasing the damping constant obviously decreased the amplitude variance of the in-phase output and improved the quadrature-phase output linearity, which is better for the stability of the transient response and voltage/frequency swing.

Three Damping Constants in SOGI-FDE Unit
The damping constant of the SOGI is also in direct relation to the bandwidth of the in-phase signal passing through. The higher the damping constant, the wider the frequency band signal passes through with worse harmonic attenuation; the lower damping constant represents a narrower bandwidth and better frequency selection, but worse stability.
For the proposed schematics of Figure 8, the inputs of SOGI I and II were directly connected to the grid voltage in the αβ-frame. Hence, k1 and k2, the corresponding damping constants of SOGI I and II, are key to the transmitting voltage bandwidth, letting k1 = k2 = k. The Chinese shipboard power system standard stipulates that the fundamental frequency of the vessel's grid fo is (50 ± 5) Hz [10]. Hence, to ensure that the proposed PLL worked normally across a wide range, the frequency bandwidth ΔfBW of SOGI I and SOGI II should be set greater than 25 Hz, namely, Above all, the frequency variation of the input signal is the key point of amplitude variance generation in the in-phase and quadrature-phase outputs, and the damping constant is the key parameter for the SOGI unit, which directly affects the two output signals. Increasing the damping constant obviously decreased the amplitude variance of the in-phase output and improved the quadrature-phase output linearity, which is better for the stability of the transient response and voltage/frequency swing.

Three Damping Constants in SOGI-FDE Unit
The damping constant of the SOGI is also in direct relation to the bandwidth of the in-phase signal passing through. The higher the damping constant, the wider the frequency band signal passes through with worse harmonic attenuation; the lower damping constant represents a narrower bandwidth and better frequency selection, but worse stability.
For the proposed schematics of Figure 8, the inputs of SOGI I and II were directly connected to the grid voltage in the αβ-frame. Hence, k 1 and k 2 , the corresponding damping constants of SOGI I and II, are key to the transmitting voltage bandwidth, letting k 1 = k 2 = k. The Chinese shipboard power system standard stipulates that the fundamental frequency of the vessel's grid f o is (50 ± 5) Hz [10]. Hence, to ensure that the proposed PLL worked normally across a wide range, the frequency bandwidth ∆f BW of SOGI I and SOGI II should be set greater than 25 Hz, namely, Meanwhile, for the dynamic response requirement of the SOGI unit in the SPS, the settling time should be limited to one fundamental period, namely, from Equation (7). Hence, for the furthest minimum, the in-phase output variance of SOGI I and II, shown in Figure 10a, that integrated Equations (27) and (28), k 1 and k 2 was set to 1.6. Similarly, SOGI III generated the positive rotation factor j with the quadrature-phase output. To obtain the minimum error shown in Figure 10b, and the optimal dynamic response, k 3 was set to 1.2.

Damping Constant and Error Gain in SOGI-FLL Unit
For the bandwidth requirement, SOGI-FLL should have a narrow bandwidth when compared to SOGI I and II. Hence, the damping constant k 4 should be set smaller than k 1 and k 2 from Equation (27). Let k 4 equal the square root of 2 to obtain the best response as an optimum second order system.
In the mentioned shipboard PV-diesel-battery power system, grid voltage stability is regulated by the exciter system of the diesel-engine synchronous generator. However, due to the large inertia of the diesel engine, grid voltage was in LF fluctuation on a large-scale time span. Due to the SPS with the slack bus, the fast response of the PV-battery inverter did not benefit system stability. For a good tradeoff between the speed of adaptive frequency and the stability of the power system, the rising time of the frequency system (shown in Figure 4) was made greater than the time of the fundamental period. Hence, the error gain γ was namely deduced from Equation (19) as (34) Figure 12 depicts the voltage/frequency dynamic response of SOGI-FLL when k 4 was equal to 1.414. When the grid frequency jumped from 50 to 55 Hz at 0.2 s, from 55 to 45 Hz at 0.35 s, and back to 50 Hz at 0.5 s; γ was equal to ω 2 o /2, ω 2 o /π, and ω o 2 /4, respectively. The corresponding frequency dynamic response is shown in Figure 12a, where γ is the key parameter. Higher γ brought about a short settling time with higher overshoot (see the curve of γ = ω o 2 /4), whereas a lower γ obtained good stability without overshoot, but had a longer settling time, as depicted in the curve of γ = ω o 2 /2 with a 50 ms settling time. When γ was ω o 2 /π, it showed good dynamic characteristics without overshoot and good stability. The optimal value of γ is namely (35) Figure 12b shows the estimated amplitude and frequency transient response of SOGI-FLL when the input signal's amplitude jumped from 1 to 1.2 per unit, or pu, at 0.2 s, from 1.2 to 0.8 pu at 0.35 s, and back to 1 pu at 0.5 s in γ = ω 2 o /π. In this plot, we can see that V g was equal to V, and f g was constant regardless of whether V changed or not, which means that V cannot affect V g and f g in steady states, which made the grid voltage V and frequency f have good consistency when the voltage/frequency was in swing. Moreover, f g had a slight fluctuation of 1.5 Hz and less than 50 ms settling time, about two grid periods. V g could smoothly shift in one grid period without oscillation.

Experiment and Discussion
Experimental results were obtained for a shipboard PV-diesel-battery hybrid power system where a PV-battery inverter generated a sinusoidal current synchronous to the positive sequence grid voltage. These results validated the proposed SOGI-based voltage/frequency tracking scheme in a practical shipboard power system with adverse grid conditions.

Experimental Setup
To confirm the theoretical analysis presented, the experimental setup of the shipboard PVdiesel-battery hybrid power system considered in this work is shown in Figure 13. It included a dc programmable source configured to emulate the PV-battery (TC.P. 16.1000.400.S, Regatron, Rorschach, Switzerland), a 2.2 kVA Danfoss three-phase inverter, and a programmable three-phase voltage source (61511, Chroma, Taoyuan, Taiwan) to emulate the shipboard grid and obtain repetitive unbalance faults with harmonics and dc offset. The setup used a digital control platform based on dSPACE (1007, dSPACE, Paderborn, Germany) to realize the proposed SOGI-based voltage/frequency tracking schematic and droop control strategies. The setup parameters are shown in Table 1.

Experiment and Discussion
Experimental results were obtained for a shipboard PV-diesel-battery hybrid power system where a PV-battery inverter generated a sinusoidal current synchronous to the positive sequence grid voltage. These results validated the proposed SOGI-based voltage/frequency tracking scheme in a practical shipboard power system with adverse grid conditions.

Experimental Setup
To confirm the theoretical analysis presented, the experimental setup of the shipboard PV-diesel-battery hybrid power system considered in this work is shown in Figure 13. It included a dc programmable source configured to emulate the PV-battery (TC.P. 16.1000.400.S, Regatron, Rorschach, Switzerland), a 2.2 kVA Danfoss three-phase inverter, and a programmable three-phase voltage source (61511, Chroma, Taoyuan, Taiwan) to emulate the shipboard grid and obtain repetitive unbalance faults with harmonics and dc offset. The setup used a digital control platform based on dSPACE (1007, dSPACE, Paderborn, Germany) to realize the proposed SOGI-based voltage/frequency tracking schematic and droop control strategies. The setup parameters are shown in Table 1.

Experimental Results and Discussion
The pre-fault grid voltages were V a = 220∠0 • V, V b = 220∠−120 • V, and V c = 220∠120 • V. The SPS grid voltages under adverse conditions consisted of 90% positive and 20% negative sequence components, namely V a = 242∠0 • V, V b = 180.4∠−132 • V, and V c = 180.4∠132 • V with 22 V, −11 V, and −11 V dc bias values, respectively, including low-order harmonics content (5% 5th and 5% 7th harmonics) to emulate the most serious situations, which led to a measured voltage THD in the pre-fault of 6.43% and 10.10% for the α and β frames, respectively. The programmed voltage swing of ±20% and frequency oscillation of ±20% is shown in Figure 14 to emulate the adverse grid voltage of SPS in load transient for its power limitation.

Experimental Results and Discussion
The pre-fault grid voltages were Va = 220∠0° V, Vb = 220∠−120° V, and Vc = 220∠120° V. The SPS grid voltages under adverse conditions consisted of 90% positive and 20% negative sequence components, namely Va = 242∠0° V, Vb = 180.4∠−132° V, and Vc = 180.4∠132° V with 22 V, −11 V, and −11 V dc bias values, respectively, including low-order harmonics content (5% 5th and 5% 7th harmonics) to emulate the most serious situations, which led to a measured voltage THD in the prefault of 6.43% and 10.10% for the α and β frames, respectively. The programmed voltage swing of ±20% and frequency oscillation of ±20% is shown in Figure 14 to emulate the adverse grid voltage of SPS in load transient for its power limitation. The parameters of this proposed adaptive grid voltage/frequency tracking schematic based on SOGIs to the shipboard PV-diesel-battery hybrid power system, shown in Figure 8, were optimized and are listed in Table 2.  The parameters of this proposed adaptive grid voltage/frequency tracking schematic based on SOGIs to the shipboard PV-diesel-battery hybrid power system, shown in Figure 8, were optimized and are listed in Table 2.  Figure 15 shows the grid voltage amplitude and frequency tracking transient response in the emulated adverse grid. Figure 15a,b depict the dynamic response of the locked amplitude V g and frequency f g in amplitude transient and show the filtered positive sequence of grid voltage in the αβ-frame with good quality without dc offset and harmonics. Note that the proposed SOGI-FDE-PLL could filter dc bias and harmonics and extract a 0.9 pu positive sequence voltage. The estimated grid voltage amplitude V g presented a fast response in one fundamental period without overshoot and had smooth dynamic features, whereas f g had oscillation in 0.25 Hz and the settling time was 2~3 periods, which can be neglected for the SPS with a ±5% fluctuation permit. Figure 15c,d depict the dynamic response of V g and f g in the alternative frequencies of 40 and 60 Hz, as well as the smoothly positive sequence of the grid voltage. The locked frequency f g had a smooth response with a 0.5 ms settling time without any fluctuations that have adaptive features of auto-adjustment, whereas the locked voltage V g oscillated with an 0.2 pu voltage deviation in approximately two periods settling time. For the low capacity of the PV-battery inverter when compared to the diesel generation system, this voltage deviation in two periods had little effect on the grid voltage. Hence, the proposed SOGI-FDE-FLL can capture grid voltage amplitude and frequency accurately and quickly in the vessel's power quality problems of serious voltage/frequency swing, unbalanced voltage, subharmonic/dc-offset voltage, and LF/HF harmonics. periods, which can be neglected for the SPS with a ±5% fluctuation permit. Figures 15c and 15d depict the dynamic response of Vg and fg in the alternative frequencies of 40 and 60 Hz, as well as the smoothly positive sequence of the grid voltage. The locked frequency fg had a smooth response with a 0.5 ms settling time without any fluctuations that have adaptive features of auto-adjustment, whereas the locked voltage Vg oscillated with an 0.2 pu voltage deviation in approximately two periods settling time. For the low capacity of the PV-battery inverter when compared to the diesel generation system, this voltage deviation in two periods had little effect on the grid voltage. Hence, the proposed SOGI-FDE-FLL can capture grid voltage amplitude and frequency accurately and quickly in the vessel's power quality problems of serious voltage/frequency swing, unbalanced voltage, subharmonic/dc-offset voltage, and LF/HF harmonics.  The emulated shipboard PV-diesel-battery hybrid power system (see Figure 13) was set up to test the validity of the proposed method. Figure 16 depicts the details of the emulated grid voltage v gabc and the a-phase injection-current i oa of the PV-battery inverter in rated power at a steady state. The THD of this obtained current was 2.05% (<3%) with higher quality, although the shipboard grid had more power quality issues of unbalanced voltage, subharmonic/dc-offset voltage, and LF/HF harmonics. The emulated shipboard PV-diesel-battery hybrid power system (see Figure 13) was set up to test the validity of the proposed method. Figure 16 depicts the details of the emulated grid voltage vgabc and the a-phase injection-current ioa of the PV-battery inverter in rated power at a steady state. The THD of this obtained current was 2.05% (<3%) with higher quality, although the shipboard grid had more power quality issues of unbalanced voltage, subharmonic/dc-offset voltage, and LF/HF harmonics.
(a) (b) Figure 16. Measured grid voltage and injection currents of the PV-battery inverter in a steady state: (a) the emulated adverse grid voltage: 100 V/div, 10 ms/div; (b) the corresponding a-phase injection current ioa of the inverter: 2 A/div, 10 ms/div.

Conclusions
In this paper, a kind of shipboard PV-diesel-battery power system structure was chosen, where the diesel-engine synchronous generator worked in standalone mode as a slack bus to maintain the grid voltage frequency and amplitude. PQ control was adopted for the PV-battery inverter, and the integrated power droop control strategy was established. Based on this SPS, an adaptive grid voltage/frequency tracking configuration (SOGI-FDE-FLL) was proposed for unbalanced voltage, subharmonic/dc-offset voltage, and LF/HF harmonic conditions. An experimental setup was established. When the emulated grid had an unbalanced voltage with 0.2 pu negative component with a 20% voltage swing and frequency oscillation including the 0.05 pu dc bias, 5th and 7th LF harmonics separately, the positive sequence voltage amplitude and frequency could be obtained precisely with a 2~3-period step response. Moreover, the captured voltage amplitude has consistence to the grid voltage amplitude envelope curve. Hence, the experimental results testified that the proposed SOGI-FDE-FLL could capture the grid voltage amplitude and frequency accurately and quickly. This proposed grid voltage/frequency tracking method based on SOGIs has a simple structure, tracks transient and accurate grid voltage/frequency, and has good prospects for vessels applying photovoltaic power generation technology. (a) the emulated adverse grid voltage: 100 V/div, 10 ms/div; (b) the corresponding a-phase injection current i oa of the inverter: 2 A/div, 10 ms/div.

Conclusions
In this paper, a kind of shipboard PV-diesel-battery power system structure was chosen, where the diesel-engine synchronous generator worked in standalone mode as a slack bus to maintain the grid voltage frequency and amplitude. PQ control was adopted for the PV-battery inverter, and the integrated power droop control strategy was established. Based on this SPS, an adaptive grid voltage/frequency tracking configuration (SOGI-FDE-FLL) was proposed for unbalanced voltage, subharmonic/dc-offset voltage, and LF/HF harmonic conditions. An experimental setup was established. When the emulated grid had an unbalanced voltage with 0.2 pu negative component with a 20% voltage swing and frequency oscillation including the 0.05 pu dc bias, 5th and 7th LF harmonics separately, the positive sequence voltage amplitude and frequency could be obtained precisely with a 2~3-period step response. Moreover, the captured voltage amplitude has consistence to the grid voltage amplitude envelope curve. Hence, the experimental results testified that the proposed SOGI-FDE-FLL could capture the grid voltage amplitude and frequency accurately and quickly. This proposed grid voltage/frequency tracking method based on SOGIs has a simple structure, tracks transient and accurate grid voltage/frequency, and has good prospects for vessels applying photovoltaic power generation technology.