WO2013122610A1 - Maximum power point tracking (mppt) - Google Patents

Maximum power point tracking (mppt) Download PDF

Info

Publication number
WO2013122610A1
WO2013122610A1 PCT/US2012/025739 US2012025739W WO2013122610A1 WO 2013122610 A1 WO2013122610 A1 WO 2013122610A1 US 2012025739 W US2012025739 W US 2012025739W WO 2013122610 A1 WO2013122610 A1 WO 2013122610A1
Authority
WO
WIPO (PCT)
Prior art keywords
power
array
voltage
predictive model
arrays
Prior art date
Application number
PCT/US2012/025739
Other languages
French (fr)
Inventor
Djordje Garabandic
Original Assignee
Xantrex Technology, Inc.
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Xantrex Technology, Inc. filed Critical Xantrex Technology, Inc.
Priority to CN201280072419.7A priority Critical patent/CN104220951B/en
Priority to EP12868625.0A priority patent/EP2815286A4/en
Priority to AU2012370020A priority patent/AU2012370020A1/en
Priority to US14/379,030 priority patent/US20150013748A1/en
Priority to PCT/US2012/025739 priority patent/WO2013122610A1/en
Publication of WO2013122610A1 publication Critical patent/WO2013122610A1/en

Links

Classifications

    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02SGENERATION OF ELECTRIC POWER BY CONVERSION OF INFRARED RADIATION, VISIBLE LIGHT OR ULTRAVIOLET LIGHT, e.g. USING PHOTOVOLTAIC [PV] MODULES
    • H02S50/00Monitoring or testing of PV systems, e.g. load balancing or fault identification
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R21/00Arrangements for measuring electric power or power factor
    • G01R21/133Arrangements for measuring electric power or power factor by using digital technique
    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05FSYSTEMS FOR REGULATING ELECTRIC OR MAGNETIC VARIABLES
    • G05F1/00Automatic systems in which deviations of an electric quantity from one or more predetermined values are detected at the output of the system and fed back to a device within the system to restore the detected quantity to its predetermined value or values, i.e. retroactive systems
    • G05F1/66Regulating electric power
    • G05F1/67Regulating electric power to the maximum power available from a generator, e.g. from solar cell
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02EREDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E10/00Energy generation through renewable energy sources
    • Y02E10/50Photovoltaic [PV] energy
    • Y02E10/56Power conversion systems, e.g. maximum power point trackers

Definitions

  • the process of converting renewable energy into power is generally complex, as many external factors play a role in how that energy is harvested and utilized.
  • conditions are analyzed to try to maximize the harvested energy.
  • a solar panel implemented using a photovoltaic (PV) array(s) can collect more energy when positioned perpendicularly to the sun on a cloudless day.
  • PV photovoltaic
  • maximum power point tracking procedures may be utilized.
  • the maximum power point tracking (MPPT) procedures may determine the supply and demand of the PV array(s) and the load (e.g., appliances, machines, battery), and may be used to identify the optimal or near optimal voltage/current of the PV array at which power extraction may be maximized (or nearly maximized).
  • MPP maximum power point
  • a method for power tracking includes measuring a plurality of samples of power produced by a photovoltaic (PV) array over a first interval of time, determining, based on the measured plurality of samples of the power of the PV array, a non-linear predictive model of a behavior of the power produced by the PV array, and performing a first adjustment of the PV array's voltage at a second time instant subsequent to an end of the first interval of time.
  • PV photovoltaic
  • the method further includes measuring at a third time instant another sample of the power of the PV array with the adjusted voltage, and determining a power difference between the power of the PV array with the adjusted voltage at the third time instant and a computed power level of the PV array at the third time instant determined from the non-linear predictive model.
  • Embodiments of the method may include at least some of the features described in the present disclosure, including one or more of the following features.
  • the method may further include performing another adjustment of the voltage of the PV array based on the determined power difference between the power of the PV array with the adjusted voltage at the third time instant and the computed power level of the PV array at the third time instant determined from the non-linear predictive model.
  • Performing the other adjustment of the voltage of the PV array may include increasing the voltage of the PV array in response to a determination that the power difference is greater than ⁇ when the first adjustment of the voltage of the PV array was a voltage increase, wherein ⁇ is a pre-determined value smaller than power produced by the PV array.
  • Performing the other adjustment of the voltage of the PV array may include decreasing the voltage of the PV array in response to a determination that the power difference is less than ⁇ when the first adjustment of the voltage of the PV array was a voltage increase, wherein ⁇ is a pre-determined value smaller than power produced by the PV array.
  • Determining the non-linear predictive model of the behavior of the power produced by the PV array may include deriving the non-linear predictive model using one or more of, for example, interpolation polynomials, a least square technique, and/or trigonometric interpolation.
  • Measuring the plurality of samples of the power of the PV array over the first interval of time may include measuring the plurality of samples at irregular time periods during the first interval of time.
  • the non-linear predictive model of the behavior of the power of the PV array may include a non-linear predictive model of the behavior of the power of the PV array resulting from insolation changes.
  • the insolation changes may include one or more of, for example, a temperature change, a cloud cover change, and/or a sun position change.
  • a system in some embodiments, includes one or more photovoltaic (PV) arrays configured to convert solar irradiation into current, a sampling subsystem to measure one or more samples of power produced by the one or more PV arrays, and a controller.
  • PV photovoltaic
  • the controller is configured to determine, based on a plurality of samples of the power of the one or more PV array measured over a first time interval, a non-linear predictive model of a behavior of the power produced by the one or more PV arrays, cause a first adjustment of the one or more PV arrays' voltage at a second time instant subsequent to an end of the first interval of time, and determine a power difference between another sample of the power of the one or more PV arrays with the adjusted voltage measured at a third time instant and a computed power level of the one or more PV arrays at the third time instant determined from the non-linear predictive model.
  • Embodiments of the system may include at least some of the features described in the present disclosure, including at least some of the features described above in relation to the method, as well as one or more of the following features.
  • the controller may further be configured to cause another adjustment of the voltage of the one or more PV arrays based on the determined power difference between the power of the one or more PV arrays with the adjusted voltage at the third time instant and the computed power level of the one or more PV arrays at the third time instant determined from the non-linear predictive model.
  • the controller configured to cause the other adjustment of the voltage of the one or more PV arrays may be configured to cause an increase of the voltage of the one or more PV arrays in response to a determination that the power difference is greater than ⁇ when the first adjustment of the voltage of the one or more PV arrays was a voltage increase, wherein ⁇ is a pre-determined value smaller than power produced by the one or more PV arrays.
  • the controller configured to cause the other adjustment of the voltage of the one or more PV arrays may be configured to cause a decrease of the voltage of the PV array in response to a determination that the power difference is less than ⁇ when the first adjustment of the voltage of the one or more PV arrays was a voltage increase, wherein ⁇ is a predetermined value smaller than power produced by the one or more PV arrays.
  • the controller configured to determine the non-linear predictive model of the behavior of the power produced by the one or more PV arrays may be configured to derive the non- linear predictive model using one or more of: interpolation polynomials, a least square technique, and trigonometric interpolation.
  • the controller may include one or more processors.
  • the sampling subsystem configured to measure the one or more samples of the power produced by the one or more PV arrays may be configured to measure the plurality of samples at irregular time periods during the first interval of time.
  • a method for power tracking includes measuring a plurality of samples of power produced by a photovoltaic (PV) array over a first interval of time, determining, based on the measured plurality of samples of the power of the PV array, a non-linear predictive model of a behavior of the power produced by the PV array, and performing a first adjustment of the PV array's current at a second time instant subsequent to an end of the first interval of time.
  • PV photovoltaic
  • the method further includes measuring at a third time instant another sample of the power of the PV array with the adjusted current, and determining a power difference between the power of the PV array with the adjusted current at the third time instant and a computed power level of the PV array at the third time instant determined from the non-linear predictive model.
  • Embodiments of the method may include at least some of the features described in the present disclosure, including at least some of the features described above in relation to the first method and the system, as well as one or more of the following features.
  • the method may further include performing another adjustment of the current of the PV array based on the determined power difference between the power of the PV array with the adjusted current at the third time instant and the computed power level of the PV array at the third time instant determined from the non-linear predictive model.
  • a method for power tracking includes measuring a plurality of samples of power produced by a photovoltaic (PV) array over a first interval of time, and determining, based on the measured plurality of samples of the power of the PV array, a linear predictive model of a behavior of the power produced by the PV array, the linear predictive model derived using a least square technique.
  • PV photovoltaic
  • the method further includes performing a first adjustment of the PV array's voltage at a second time instant subsequent to an end of the first interval of time, measuring at a third time instant another sample of the power of the PV array with the adjusted voltage, and determining a power difference between the power of the PV array with the adjusted voltage at the third time instant and a computed power level of the PV array at the third time instant determined from the linear predictive model derived using the least square technique.
  • Embodiments of the method may include at least some of the features described in the present disclosure, including at least some of the features described above in relation to the methods and the system.
  • a list of "at least one of A, B, and C” includes any of the combinations A or B or C or AB or AC or BC and/or ABC (i.e., A and B and C).
  • a list of "at least one of A, B, and C” may also include AA, AAB, AAA, BB, etc.
  • FIG. 1 A is a schematic diagram of a solar panel system to harvest solar energy.
  • FIG. IB is an example graph of power output behavior for a PV array, such as the PV array shown in FIG. 1A, as a function of the array's voltage
  • FIG. 2 is a graph of the power behavior of a PV array over a period of time as a result of insolation changes and/or voltage perturbation.
  • FIG. 3 is a flow chart of a power tracking procedure that takes into account the effect of insolation changes on the power tracking functionality performed by, for example, the system of FIG. 1A.
  • FIG. 4 is a schematic diagram of an example system implementation to perform a tracking procedure such as the procedure of FIG. 3.
  • FIG. 5 is a schematic diagram of a generic computing system.
  • FIG. 6 is a schematic diagram of an example implementation to determine the power difference between the power measured following perturbation and the value computed using a predictive model based on a second order Lagrange polynomial.
  • FIG. 7 is a graph illustrating the performance of a tracking procedure.
  • a solar power system includes a device that converts an electrical direct current produced in a PV array into an alternating current.
  • the alternating current is then provided as an input into a computer-based system, which is configured, among other things, to track the MPP to enable efficiently harvesting power from the PV array.
  • the computer-based system may further include a clock that is utilized to time intervals of the signal sampling period, intervals of the perturb cycle period, a perturb point, and an observe point.
  • the alternating current is sampled at intervals, which are less than the interval of the perturb cycle, in order to generate a predictive model of power behavior of the array from which a trajectory of the predicted behavior of the PV array is generated.
  • the difference between the actual power of the PV array and the predicted power at a subsequent time instant is determined.
  • the MPP can be properly tracked in dynamic insolation conditions by allowing the direction of the next perturbation cycle to be correctly chosen based on the determined power difference.
  • Other system implementations are possible.
  • methods, systems, and other implementations including a method for power tracking that includes measuring a plurality of samples of power produced by a photovoltaic (PV) array over a first interval of time, determining, based on the measured plurality of samples of the power of the PV array, a non-linear predictive model of a behavior of the power produced by the PV array, and performing a first adjustment of the PV array's voltage at a second time instant subsequent to an end of the first interval of time.
  • PV photovoltaic
  • the method further includes measuring another sample of the power of the PV array with the adjusted voltage, and determining a power difference between the power of the PV array with the adjusted voltage at a third time instant and a computed power level of the PV array at the third time instant determined from the non-linear predictive model.
  • the system 100 includes a plurality of photovoltaic (PV) modules (or cells), such as the PV module 104, arranged as a PV array 103.
  • PV photovoltaic
  • each PV module may comprise a connected series of photovoltaic cells.
  • Surfaces of the PV cells which, in some implementations, may be configured similarly to photodiodes, are irradiated with radiation 101 from the sun 102 to excite electrons and create a potential between two electrodes.
  • the energy produced by the cells creates a direct current (DC).
  • DC direct current
  • the array 103 is connected to an inverter 1 10 (also referred to as a converter) configured to convert DC to AC.
  • Solar power may be measured in insolation, which is the amount of power (kilowatt hours) collected over a known surface area (square meters) and over a given time period (day).
  • insolation changes with changes to environmental conditions, such as temperature, cloud cover and position of the sun, and thus, the power produced by the PV modules fluctuates (generally in a non-linear manner) when changes to the environment in which the PV array 103 is deployed occur.
  • the PV array 103 may be mounted on a moveable support 105, which can be adjusted to position the planar surface of the array at an angle with the sun that seeks to maximize the radiation received from the sun.
  • the optimal voltage or current of the PV array varies.
  • the operating points of the PV array may be adjusted so as to identify the optimal (or near optimal) operating point for the PV array for the present load and/or environmental conditions.
  • MPPT maximum power point tracking
  • the MPPT procedure is configured to determine, based on the supply and demand of the PV array and the load, an optimal voltage (when the PV array is configured to operate with a constant voltage and a variable current that varies with changing environmental and/or load conditions), and/or the optimal current (when the PV array is configured to operate with a constant current and a variable voltage level that varies with changing environmental and/or load conditions) to enable optimal, e.g., maximum or near maximum, power to be provided by the PV array.
  • MPPT maximum power point tracking
  • IB showing an example graph 150 of power output behavior for a PV array as a function of the array's voltage
  • the point where the difference between the supply and demand is zero, or approximately zero corresponds to the maximum power point (MPP) of the PV array where the PV array is supplying power to the load substantially efficiently, and little or no additional power is being lost or required.
  • MPP maximum power point
  • the point on the graph at which maximum power is delivered corresponds to a point where the derivative of the graph, dP/dV, is equal to 0.
  • the derivative dP/dV (as V increases) is not equal to 0, and the PV array is thus not operating at an optimal state.
  • a section 154 of the graph 150 corresponds to a derivative dP/dV being smaller than 0, and may be indicative that voltage levels are too high for given load conditions.
  • points on the graph 150 where the derivative dP/dV is greater than 0 (dP/dV > 0) as V increases, as depicted in the section 156 of the graph, may represent voltage levels that are too low for existing load conditions.
  • a procedure such as the Maximum Power Point Tracking (MPPT) method, may be implemented to iteratively adjust the voltage level of the PV array until optimal or near optimal power delivery by the PV array for the given load and/or environmental conditions is established.
  • MPPT Maximum Power Point Tracking
  • an MPPT-based procedure may be implemented as a "perturb and observe” (P&0)-based procedure, in which the voltage of the PV array is periodically perturbed (e.g., the voltage produced by the PV array is adjusted), and the change in the current, and thus in the power, provided by the PV array as a result of the change in voltage, is observed.
  • P&0 perturb and observe
  • This procedure can be described according to the following equation:
  • Ppv (T) - Pp V (0) x (1)
  • Ppy the observed power of the photovoltaic array
  • T the period of the perturb interval (in this example)
  • the system 100 depicted in FIG. 1A may include a controller 120 (e.g., a processor-based controller, a DSP processor, an application-specific integrated circuit, and/or any combination thereof) configured to determine whether non-optimal conditions of the systems and devices exist (e.g., by tracking the behavior of the power produced as a function of changing voltage and/or current), and to controllably adjust the voltage/current to converge to an optimal or near optimal power behavior for the PV array.
  • a controller 120 e.g., a processor-based controller, a DSP processor, an application-specific integrated circuit, and/or any combination thereof
  • the controller is configured to cause an adjustment of the voltage of PV array, and to repeat these operations if the voltage adjustment did not result in establishment of optimal or near- optimal conditions (e.g., MPP was not established).
  • the controller 120 may be electrically coupled to, for example, the output of the inverter 110 (e.g., to enable determining power values) and/or to the PV array 103 (e.g., to control the array's voltage and/or current).
  • MPPT Procedures implemented in accordance with, for example, Equation 1 above, measure the power (e.g., multiplying the array's voltage by the array's current) at one time instant (e.g., when the system was perturbed by, for example, causing the voltage of the array to increase or decrease), and at a subsequent time instant, T (the time period between the time instants should be sufficiently long to enable the power behavior of the array to settle down).
  • the computed power difference according to Equation 1 does not indicate what contributed to the power measured at the subsequent time instant. Particularly, if during the period of measurement environmental and load conditions were unchanged, then it can be assumed that the power measured at time instant T resulted substantially from the adjustment made to the voltage of the PV array. However, if there were changes to the environmental conditions in which the PV array is operating (temperature, cloud coverage, sun brightness, etc.) then these environmental changes may have contributed, at least to some extent, to the measured change in the power supplied by the PV array (for the purpose of explanation, it will be assumed that load conditions are not changed during the PV array power measurement period).
  • power optimization procedures for PV arrays that take into account the effect/contribution of insolation changes may be implemented.
  • computation of the power difference may be performed in accordance with a model that incorporates the effect/contribution of insolation changes to changes in the power supplied/produced by a PV array.
  • FIG. 2 a graph 200 of the power behavior of a PV array over a period of time as a result of insolation changes and voltage perturbation is shown. As illustrated in FIG.
  • the power is shown to vary due to, in the example of the graph 200, insolation changes (e.g., changes in the temperature, cloud coverage, etc.).
  • insolation changes e.g., changes in the temperature, cloud coverage, etc.
  • the power behavior of the PV array resulting from insolation changes may be sampled, and based on the sample points recorded, a model representative of power behavior as a result of insolation changes is determined.
  • the derived model can therefore be used to predict the power behavior of the PV array, as a result of insolation changes, at a time subsequent to the end of the sampling period.
  • the dashed projection 204 is a predicted extension of sampled portion of the graph 200 prior to its perturbation, and therefore does not represent the behavior of the power level of the PV array resulting from perturbation of the PV array.
  • the PV array system is perturbed (e.g., its voltage level is altered), resulting in power behavior represented, in the example of FIG. 2, by the solid line portion 210 of the graph 200.
  • the solid line portion 210 represents the behavior of the power of the PV array resulting from the effect of both the perturbation and the insolation changes.
  • the model derived from the samples can provide reasonable short term predictions of the power behavior from insolation changes (or any other factor which is to be modeled), but the accuracy and validity of the predicted behavior diminishes over time. As can be seen from FIG.
  • the effect the perturbation has on the power behavior of the PV array can be determined by identifying, using the derived model (represented by the dashed portion 204 depicted in FIG. 2), the effect insolation changes are predicted to have on the power behavior, and subtracting from the overall power behavior (represented by the solid line section 210 of the graph 200) the predicted effect/contribution of the insolation changes on the power behavior. Accordingly, in some embodiments, the effect of perturbation can be represented as:
  • T is a time instant, subsequent to the end of the sampling period and to the time instant when the PV array was perturbed, when the power supplied by the PV array is sampled.
  • the resultant value will therefore provide a more accurate indication of dP/dV, and will thus provide a more accurate and reliable tracking of the power behavior of the PV array so as to more quickly achieve an optimal, or near optimal, voltage level (or current level, in situations where the PV array has a near constant current and fluctuating voltage).
  • the model derived to predict the power behavior of the PV array as a result of insolation changes may be a non-linear model, e.g., high order models such as, for example, second and/or third order models.
  • non-linear models provide a more accurate and robust representation (e.g., as compared to representations provided, for example, by linear models) of a rate of change (in this case, the PV array's power rate of change), and therefore may enable more accurate tracking of the PV array's MPP.
  • non-linear model enables non-constant sampling (i.e., using varying sampling intervals in which the sampling intervals between successive samples are not necessarily the same) to obtain the samples that are used to derive the non-linear representation corresponding to the samples.
  • derivation of non-linear predictive models to represent the behavior of PV array power resulting from insolation changes (or resulting from other factors) may be based on one or more of, for example, interpolation polynomials, a least square technique, and/or trigonometric interpolation.
  • a flow chart of a power tracking procedure 300 that takes into account the effect of insolation changes on the power tracking functionality performed by, for example, the system 100, is shown.
  • Implementation of the procedure 300 may be based, for example, on configuring the controller 120 depicted in FIG. 1 to perform the operations described herein.
  • a plurality of samples of power produced by the PV array over a first interval of time is measured 310.
  • Sample measurements can be performed, for example, by separately measuring the voltage and current produced by the PV array (e.g., using voltmeter, ammeters, etc.) and multiplying the two to obtain an instant power value supplied by the PV array.
  • the measured electrical attributes may be performed with respect to the output of the inverter (e.g., the inverter 1 10) of the PV array.
  • Measuring power samples may also be performed, for example, by dividing the identification period (also called the "sampling period) into a number of intervals. In some embodiments, and as shown in FIG. 2, the identification period may be divided to substantially equal length intervals of duration T s .
  • the sampling period does not have to be constant, and the intervals between successive samples may vary over the duration of the identification period.
  • the measured samples may be obtained by using a sample-and-hold circuit to determine the voltage level and/or current level at the measurement instants, and using, for example, an analog-to-digital converter (ADC) to digitize the measured samples.
  • ADC analog-to-digital converter
  • Such sample-and-hold circuit(s) and/or an ADC may be part of the example controller 120 of FIG. 1.
  • a non-linear predictive model of a behavior of the power produced by the PV array may be determined 320.
  • a non-linear model based on interpolation polynomials may be derived.
  • Various forms may be used to represent interpolating polynomials.
  • an interpolating polynomial function may be represented as a linear combination of orthogonal polynomials.
  • an interpolating polynomial may be represented as combination of Lagrangian polynomials.
  • a combination of Lagrangian polynomials may be represented as:
  • the sampling interval is substantially constant, a non-constant (irregular) sampling interval may be used instead.
  • Equation (5) can thus be utilized to compute the sampled signal power at any number of intervals " «" to collect data points which predict the trajectory of the PV power signal.
  • the integration of the Equation (5) provides the predictive PV power signal function Ppv modei(t) having n+1 samples over a period T s represented as: Ppv_modei ((n+l)T s )
  • Equation (6) each € (n) can be derived from equation (5), with each selection of "n” producing "n+1" data points.
  • FIG. 6 is a schematic diagram of an example DSP implementation 600 to determine the power difference between the power (or a discretized value thereof) measured following perturbation (i.e., P(k)) and the value computed using a second order Lagrange polynomial used to predict the power behavior of a PV array resulting, for example, from insolation changes.
  • the power differential can be expressed as:
  • the implementation 600 includes a summing operator 610 that adds the sample 602, corresponding to P(k), to the negative of the preceding sample 604 (corresponding to P(k-l/3) that was amplified by an amplifier 612 with a gain of 3.
  • the resultant sum of the summing operation at 610 is added, using the summing operator 616, to a resultant amplified valued of the sample 606, corresponding to P(k-2/3).
  • the sample 606 is amplified using an amplifier 614 that also has a gain of 3.
  • the result of the summing operator 616 is added to the negative of the sample 608, corresponding to P(k-1).
  • L(4) y 0 (4-l)(4-2)(4-3)/(-6) + yi 4(4-2)(4-3)/2 + y 2 4(4-l)(4-3)/(-2)
  • derivation of another type of non-linear predictive model to represent the power behavior due to, for example, insolation changes may be based on the Least Square Method.
  • the resulting signal model is generally a polynomial function with a function order of m where m ⁇ n +1 .
  • the parameter m may be selected based on the desired modeling accuracy.
  • the implementation of the least square method to facilitate derivation of a model based on the samples collected can be, in some embodiments, in the form of recursive calculations.
  • PV PV model can be represented as:
  • PVpv model ao + & ⁇ t where ao and ai are obtained from:
  • PVp V model may therefore be expressed as: pv pv model (yi - yo) 1 12 + 5/6 yo + 1/3 yi - l/6y 2 [0071]
  • derivation of another type of nonlinear predictive model to represent the power behavior due to, for example, insolation changes may be based on trigonometric interpolation.
  • This derivation technique is suitable, for example, in situations where the model is to be represented as a periodic function.
  • trigonometric interpolation for n+l data points, namely (xo,yo), ⁇ ⁇ ⁇ , (xn,y «).
  • a function F(x) may be represented as follows:
  • F(x) yo [sin 1 ⁇ 2 (x-xi) sin 1 ⁇ 2 (x-x 2 ) ... sin 1 ⁇ 2 (x-x n )] / [sin 1 ⁇ 2 (xo-xi) sin 1 ⁇ 2 (xo-x 2 ) .... sin 1 ⁇ 2 (x 0 -x n )] + yi [sin 1 ⁇ 2 (x-x 0 ) sin 1 ⁇ 2 (x-x 2 ) ... sin 1 ⁇ 2 (x-x n )] / [sin 1 ⁇ 2 (xi-xo) sin 1 ⁇ 2 (xi-x 2 ) ....
  • the sampling points were assumed to be taken at equal duration intervals.
  • derivations to obtain predictive model representations for the power behavior may be performed with sample points taken at non-constant (e.g., irregular) intervals.
  • other techniques/procedures to determine non-linear predictive model to represent the power behavior of the PV array may be used in addition to or instead of the any of the above illustrated techniques/procedures.
  • a first adjustment is performed 330 of the PV array's voltage at a second time instant subsequent to an end of the first interval of time (during which the power samples were measured/obtained).
  • the voltage adjustment performed i.e., perturbation of the system, may be to increase or decrease the voltage of the PV array by some pre-determined amount, or by an amount that depends on the present voltage level (as may have been measured, for example, when the last sample was obtained).
  • the determination of whether to initially increase or decrease the voltage may be based on some pre-determined set of rules (e.g., that the voltage is to initially be increased, that a decision as to initially increase or decrease the voltage may be based on the present voltage level or on the difference between power levels of the last two samples obtained, etc.)
  • the PV array may be configured to output a substantially constant current, and to have varying voltage that depends, for example, on load changes, insolation changes, etc. Under those circumstances, perturbation of the system may be performed by making a current adjustment to the PV array, and tracking the maximum power point resulting from such current perturbations.
  • T another sample of the power of the PV array with an adjusted electrical attribute (e.g., adjusted voltage, in the example of FIG. 2) is measured 340 (i.e., the observation of the system following the perturbation of the system).
  • the measurement of the power level of the PV array with the adjusted electrical attribute may be performed using a sample-and-hold circuit to determine the voltage level and/or current level at the instant of measurements, and by using, for example, an analog-to-digital converter to digitize the measured sample.
  • the time duration between "Tl" and "T” should be sufficiently long to enable the system to respond to the adjustment performed.
  • a power difference is determined 350 between the power of the PV array with the adjusted voltage at the third time instant and a computed power level of the PV array at the third time instant determined from using the non-linear predictive model.
  • a determination is made of Ppv(T) - Ppv modd(T).
  • the power difference therefore provides a value indicative of the effect the perturbation had on the power being tracked by attempting to remove the contribution that insolation changes (computed using the determined non-linear model), or any other factor represented by this or some other non-linear predictive model, had on the power being tracked.
  • the determined power difference may be utilized to track the MPP of the PV array and control the output of the PV array in relation to the load.
  • the determined power difference can be utilized to determine the next direction of perturbation (e.g., adjustment of an attribute of the PV array, be it the voltage or current) based on the sign of a result of the power difference, and to use the resultant sign (positive "+", or negative "-") to control the direction of the perturbation during the next cycle.
  • the procedures implemented to track the maximum power point of a PV array may include performing another adjustment of the voltage of the PV array based on the determined power difference between the power of the PV array with the adjusted voltage at the third time instant and the computed power level of the PV array at the third time instant as determined from the non-linear predictive model.
  • is a pre-determined value that is generally smaller than the power level of the PV array with the adjusted voltage ( ⁇ may be equal to 0 in some embodiments)
  • the preceding adjustment of the voltage (or another electrical attribute) of the PV array was a voltage increase then the next voltage adjustment for the PV array would be to increase the voltage.
  • the computed power difference is determined to be greater than ⁇ (e.g., greater than 0), this may indicate that the power is being tracked in the right direction.
  • the next adjustment should also be a voltage increase (e.g., by some pre-determined amount, by an amount that is based on, for example, the present voltage or power level, etc.) If the previous adjustment was a voltage decrease, then the next adjustment should also be a voltage decrease.
  • the next adjustment of the PV array's electrical attribute(s) may be to decrease the voltage.
  • a computed power difference that is less than ⁇ e.g., less than 0
  • the voltage may have to be decreased on the next adjustment so as to reverse to tracking direction.
  • the next adjustment may be a voltage increase.
  • the process of adjusting the electrical attribute(s) of the PV array may be repeated until, for example, an absolute value of the determined power difference is less than some pre-determined value, indicating that the PV array may have substantially or approximately reached its maximum power point.
  • FIG. 4 a schematic diagram of an example system implementation 400 to perform a tracking procedure, such as that performed by the procedure 300 of FIG. 3, is shown.
  • the system 400 may be implemented, at least in part, as the controller 120 of FIG. 1A.
  • the system 400 may include, for example, a sampling and PV power calculation module 410 (also referred to as a signal identification unit/module) to measure one or more samples of a power level produced by a PV array (such as the PV array 103), and may include sample-and-hold circuits and analog-to-digital converters.
  • the Data Sampling and PV Power Calculation module 410 samples the voltage and current of the PV generator, computes power and provides the input data samples to the Power Increment Calculation module 420.
  • the input to the Power Increment Calculation module 420 may be implemented, in some embodiments, to receive three (3) samples of the power signal taken before a perturbation (e.g., P(k-l/3), P(k-2/3) and P(k-l)), and one sample of the power signal taken after the perturbation P(k). Other implementations (e.g., using different sample sequences) may be used.
  • the Power Increment Calculation module 420 of the example system 400 uses the "before perturbation samples," e.g., P(k-l/3), P(k-2/3) and P(k-l), to construct, for example, a non- linear predictive model, and to determine a predicted (expected) power after the perturbation model, i.e., P mo dei(kT).
  • the example system 400 also includes a Perturb Direction Computation module 430 which maintains a perturbation history (e.g., previous perturbation direction). The module determines, based on the power difference dP(kT) and the previous perturbation direction, what the next perturbation direction should be. For example if the previous PV voltage perturbation was in a positive direction (direction(k-l)>0) and the power difference (also referred to as "power increment") is positive (i.e., dP(k)>0) then the next perturbation will be in a positive direction too.
  • a Perturb Direction Computation module 430 which maintains a perturbation history (e.g., previous perturbation direction). The module determines, based on the power difference dP(kT) and the previous perturbation direction, what the next perturbation direction should be. For example if the previous PV voltage perturbation was in a positive direction (direction(k-l)>0) and the power difference (also referred to as "power increment”) is positive (i.
  • the example system 400 also includes a PV Voltage Reference Calculation module 440 configured to receive the direction command (e.g., from a module/unit such as the Perturb Direction Computation module 430 of the system 400) and to compute a PV voltage reference as follows:
  • Vref(k) Vref(k-1) + V pe rturbe step
  • Vref(k) Vref(k-1) - Vperturbe step
  • Vperturbe step represents a perturbation step voltage, which may be a pre-determined step voltage, or may be an adjustable dynamically derived value.
  • Vref(k) may also be compared to an associated minimum and maximum limits and clipped to those limiting values. The minimum and maximum limits may be based on the PV generator, Inverter, and/or power converter specifications.
  • the computed Vref(k) is provided to an inverter module configured, among other things, to adjust the current of the PV generator such that the PV voltage converges to the reference voltage level.
  • a processor-based computing system which, for example, may be used to implement, at least partly, the controller 120 of FIG. 1.
  • FIG. 5 a schematic diagram of a generic computing system 500 is shown.
  • the computing system 500 includes a processor-based device 10 such as a personal computer, a specialized computing device, and so forth, that typically includes a central processor unit 512.
  • the system may include a main memory, cache memory and bus interface circuits (not shown).
  • the processor-based device 510 may include a mass storage element 514, such as a hard drive and/or a flash drive associated with the computer system.
  • the computing system 500 may further include a keyboard, or keypad, 516, and a monitor 520, e.g., a CRT (cathode ray tube) or LCD (liquid crystal display) monitor, that may be placed where a user can access them.
  • a monitor 520 e.g., a CRT (cathode ray tube) or LCD (liquid crystal display) monitor, that may be placed where a user can access them.
  • the processor-based device 510 is configured to facilitate, for example, the implementation of the tracking procedure, including facilitating measurement of sample points of the power output generated using a PV array, determining a non-linear predictive model representative of the power behavior of the PV array system resulting from such factors as insolation changes, determining whether power tracking is proceeding in the right direction, etc.
  • the storage device 514 may thus include a computer program product that when executed on the processor-based device 510 causes the processor-based device to perform operations to facilitate the implementation of the above-described procedures.
  • the processor-based device may further include peripheral devices to enable input/output functionality.
  • Such peripheral devices may include, for example, a CD-ROM drive and/or flash drive, a network connection (implemented using, for example, a USB port, a wireless transceiver, etc.), for downloading related content to the connected system.
  • Such peripheral devices may also be used for downloading software containing computer instructions to enable general operation of the respective system/device.
  • special purpose logic circuitry e.g., an FPGA (field programmable gate array), a DSP processor, or an ASIC (application-specific integrated circuit) may be used in the implementation of the system 500.
  • a DSP processor may be used to implement the non-linear predictive model used to compute the predicted power behavior of the PV array due to, for example, insolation changes (or some other factor).
  • An example of such an implementation is illustrated in FIG. 6 (as was more particularly described above).
  • Other modules that may be included with the processor-based device 510 are speakers, a sound card, a pointing device, e.g., a mouse or a trackball, by which the user can provide input to the computing system 500.
  • the processor-based device 510 may include an operating system.
  • Computer programs also known as programs, software, software applications or code
  • programs include machine instructions for a programmable processor, and may be implemented in a high-level procedural and/or object-oriented programming language, and/or in
  • machine-readable medium refers to any non-transitory computer program product, apparatus and/or device (e.g., magnetic discs, optical disks, memory, Programmable Logic Devices (PLDs)) used to provide machine instructions and/or data to a programmable processor, including a non-transitory machine- readable medium that receives machine instructions as a machine-readable signal.
  • PLDs Programmable Logic Devices
  • FIG. 7 is a graph illustrating the performance of a tracking procedure, such as the tracking procedure 300 described in relation to FIG. 3.
  • PV generator transients are simulated to include fast transients such as, for example, moving clouds, and slow transients, such as, for example, thermal effects.
  • the generator simulated can be one that can produce up to 340kW with a maximum power point voltage of around 520V.
  • the upper waveforms (in the upper box) of FIG. 7 shows the behavior of Vmppt(t) (the smooth line) and actual Vpv(t) in Volts.
  • Vmppt(t) is the "ideal" position of the maximum power point voltage and Vpv(t) is the achieved PV voltage.
  • the time axis is provided in seconds.
  • the lower waveform (lower box) of FIG. 7 shows the harvested power from the PV generator in Watts.
  • the graph shows that the PV power has some large and relatively fast transients. Those transients may be the result of cloud movements. As can be seen from the top graph, in such dynamic conditions the difference between the maximum power point voltage and the actual PV voltage are within +/-5V.

Landscapes

  • Engineering & Computer Science (AREA)
  • Power Engineering (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Sustainable Development (AREA)
  • Sustainable Energy (AREA)
  • Electromagnetism (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Automation & Control Theory (AREA)
  • Control Of Electrical Variables (AREA)
  • Supply And Distribution Of Alternating Current (AREA)

Abstract

Disclosed are methods, systems, and other implementations, including a method that includes measuring a plurality of samples of power produced by a photovoltaic (PV) array over a first interval of time, determining, based on the measured plurality of samples a non-linear predictive model of a behavior of the power produced by the PV array, and performing a first adjustment of the PV array's voltage at a second time instant subsequent to an end of the first interval of time. The method further includes measuring at a third time instant another sample of the power of the PV array with the adjusted voltage, and determining a power difference between the power of the PV array with the adjusted voltage at the third time instant and a computed power level of the PV array at the third time instant determined from the non-linear predictive model.

Description

MAXIMUM POWER POINT TRACKING (MPPT)
BACKGROUND
[0001] The process of converting renewable energy into power is generally complex, as many external factors play a role in how that energy is harvested and utilized. Typically, conditions are analyzed to try to maximize the harvested energy. For example, a solar panel, implemented using a photovoltaic (PV) array(s) can collect more energy when positioned perpendicularly to the sun on a cloudless day. As environmental conditions (e.g., cloud coverage, temperature, angle of irradiation of the sun with the solar panel, etc.) change, power extraction efficiently also changes.
[0002] Additionally, power efficiency provided by solar panels generally also varies as load conditions change. Accordingly, in order to properly extract the energy from solar panels in accordance with load conditions, maximum power point tracking procedures may be utilized. The maximum power point tracking (MPPT) procedures may determine the supply and demand of the PV array(s) and the load (e.g., appliances, machines, battery), and may be used to identify the optimal or near optimal voltage/current of the PV array at which power extraction may be maximized (or nearly maximized).
[0003] Rapid shifts in insolation conditions (e.g., changing cloud coverage, drop in temperature, changing angle of the sun with respect to the solar panels) may affect MPPT results, possibly causing erroneous MPPT results.
SUMMARY
[0004] Disclosed herein are systems, methods and other implementations that provide efficient realization of maximum power point (MPP) determination that can be achieved in rapidly changing insolation conditions (e.g., environmental conditions, such as cloud coverage, temperature, position of the sun in the sky, etc.)
[0005] Thus, in some embodiments, a method for power tracking is disclosed. The method includes measuring a plurality of samples of power produced by a photovoltaic (PV) array over a first interval of time, determining, based on the measured plurality of samples of the power of the PV array, a non-linear predictive model of a behavior of the power produced by the PV array, and performing a first adjustment of the PV array's voltage at a second time instant subsequent to an end of the first interval of time. The method further includes measuring at a third time instant another sample of the power of the PV array with the adjusted voltage, and determining a power difference between the power of the PV array with the adjusted voltage at the third time instant and a computed power level of the PV array at the third time instant determined from the non-linear predictive model.
[0006] Embodiments of the method may include at least some of the features described in the present disclosure, including one or more of the following features.
[0007] The method may further include performing another adjustment of the voltage of the PV array based on the determined power difference between the power of the PV array with the adjusted voltage at the third time instant and the computed power level of the PV array at the third time instant determined from the non-linear predictive model.
[0008] Performing the other adjustment of the voltage of the PV array may include increasing the voltage of the PV array in response to a determination that the power difference is greater than ε when the first adjustment of the voltage of the PV array was a voltage increase, wherein ε is a pre-determined value smaller than power produced by the PV array.
[0009] Performing the other adjustment of the voltage of the PV array may include decreasing the voltage of the PV array in response to a determination that the power difference is less than ε when the first adjustment of the voltage of the PV array was a voltage increase, wherein ε is a pre-determined value smaller than power produced by the PV array.
[0010] Determining the non-linear predictive model of the behavior of the power produced by the PV array may include deriving the non-linear predictive model using one or more of, for example, interpolation polynomials, a least square technique, and/or trigonometric interpolation.
[0011] Measuring the plurality of samples of the power of the PV array over the first interval of time may include measuring the plurality of samples at irregular time periods during the first interval of time.
[0012] The non-linear predictive model of the behavior of the power of the PV array may include a non-linear predictive model of the behavior of the power of the PV array resulting from insolation changes. [0013] The insolation changes may include one or more of, for example, a temperature change, a cloud cover change, and/or a sun position change.
[0014] In some embodiments, a system is disclosed. The system includes one or more photovoltaic (PV) arrays configured to convert solar irradiation into current, a sampling subsystem to measure one or more samples of power produced by the one or more PV arrays, and a controller. The controller is configured to determine, based on a plurality of samples of the power of the one or more PV array measured over a first time interval, a non-linear predictive model of a behavior of the power produced by the one or more PV arrays, cause a first adjustment of the one or more PV arrays' voltage at a second time instant subsequent to an end of the first interval of time, and determine a power difference between another sample of the power of the one or more PV arrays with the adjusted voltage measured at a third time instant and a computed power level of the one or more PV arrays at the third time instant determined from the non-linear predictive model.
[0015] Embodiments of the system may include at least some of the features described in the present disclosure, including at least some of the features described above in relation to the method, as well as one or more of the following features.
[0016] The controller may further be configured to cause another adjustment of the voltage of the one or more PV arrays based on the determined power difference between the power of the one or more PV arrays with the adjusted voltage at the third time instant and the computed power level of the one or more PV arrays at the third time instant determined from the non-linear predictive model.
[0017] The controller configured to cause the other adjustment of the voltage of the one or more PV arrays may be configured to cause an increase of the voltage of the one or more PV arrays in response to a determination that the power difference is greater than ε when the first adjustment of the voltage of the one or more PV arrays was a voltage increase, wherein ε is a pre-determined value smaller than power produced by the one or more PV arrays.
[0018] The controller configured to cause the other adjustment of the voltage of the one or more PV arrays may be configured to cause a decrease of the voltage of the PV array in response to a determination that the power difference is less than ε when the first adjustment of the voltage of the one or more PV arrays was a voltage increase, wherein ε is a predetermined value smaller than power produced by the one or more PV arrays.
[0019] The controller configured to determine the non-linear predictive model of the behavior of the power produced by the one or more PV arrays may be configured to derive the non- linear predictive model using one or more of: interpolation polynomials, a least square technique, and trigonometric interpolation.
[0020] The controller may include one or more processors.
[0021] The sampling subsystem configured to measure the one or more samples of the power produced by the one or more PV arrays may be configured to measure the plurality of samples at irregular time periods during the first interval of time.
[0022] In some embodiments, a method for power tracking is disclosed. The method includes measuring a plurality of samples of power produced by a photovoltaic (PV) array over a first interval of time, determining, based on the measured plurality of samples of the power of the PV array, a non-linear predictive model of a behavior of the power produced by the PV array, and performing a first adjustment of the PV array's current at a second time instant subsequent to an end of the first interval of time. The method further includes measuring at a third time instant another sample of the power of the PV array with the adjusted current, and determining a power difference between the power of the PV array with the adjusted current at the third time instant and a computed power level of the PV array at the third time instant determined from the non-linear predictive model.
[0023] Embodiments of the method may include at least some of the features described in the present disclosure, including at least some of the features described above in relation to the first method and the system, as well as one or more of the following features.
[0024] The method may further include performing another adjustment of the current of the PV array based on the determined power difference between the power of the PV array with the adjusted current at the third time instant and the computed power level of the PV array at the third time instant determined from the non-linear predictive model.
[0025] In some embodiments, a method for power tracking is disclosed. The method includes measuring a plurality of samples of power produced by a photovoltaic (PV) array over a first interval of time, and determining, based on the measured plurality of samples of the power of the PV array, a linear predictive model of a behavior of the power produced by the PV array, the linear predictive model derived using a least square technique. The method further includes performing a first adjustment of the PV array's voltage at a second time instant subsequent to an end of the first interval of time, measuring at a third time instant another sample of the power of the PV array with the adjusted voltage, and determining a power difference between the power of the PV array with the adjusted voltage at the third time instant and a computed power level of the PV array at the third time instant determined from the linear predictive model derived using the least square technique.
[0026] Embodiments of the method may include at least some of the features described in the present disclosure, including at least some of the features described above in relation to the methods and the system.
[0027] Unless defined otherwise, all technical and scientific terms used herein have the same meaning as commonly or conventionally understood. As used herein, the articles "a" and "an" refer to one or to more than one (i.e., to at least one) of the grammatical object of the article. By way of example, "an element" means one element or more than one element. "About" and/or "approximately" as used herein when referring to a measurable value such as an amount, a temporal duration, and the like, is meant to encompass variations of ±20% or ±10%, ±5%), or +0.1% from the specified value, as such variations are appropriate to in the context of the systems, devices, methods, and other implementations described herein.
[0028] As used herein, including in the claims, "and" as used in a list of items prefaced by "at least one of or "one or more of indicates that any combination of the listed items may be used. For example, a list of "at least one of A, B, and C" includes any of the combinations A or B or C or AB or AC or BC and/or ABC (i.e., A and B and C). Furthermore, to the extent more than one occurrence or use of the items A, B, or C is possible, multiple uses of A, B, and/or C may form part of the contemplated combinations. For example, a list of "at least one of A, B, and C" may also include AA, AAB, AAA, BB, etc.
[0029] Other and further objects, features, aspects, and advantages of the present disclosure will become better understood with the following detailed description of the accompanying drawings.
BRIEF DESCRIPTION OF THE DRAWINGS
[0030] FIG. 1 A is a schematic diagram of a solar panel system to harvest solar energy.
[0031] FIG. IB is an example graph of power output behavior for a PV array, such as the PV array shown in FIG. 1A, as a function of the array's voltage
[0032] FIG. 2 is a graph of the power behavior of a PV array over a period of time as a result of insolation changes and/or voltage perturbation. [0033] FIG. 3 is a flow chart of a power tracking procedure that takes into account the effect of insolation changes on the power tracking functionality performed by, for example, the system of FIG. 1A.
[0034] FIG. 4 is a schematic diagram of an example system implementation to perform a tracking procedure such as the procedure of FIG. 3.
[0035] FIG. 5 is a schematic diagram of a generic computing system.
[0036] FIG. 6 is a schematic diagram of an example implementation to determine the power difference between the power measured following perturbation and the value computed using a predictive model based on a second order Lagrange polynomial.
[0037] FIG. 7 is a graph illustrating the performance of a tracking procedure.
[0038] Like reference symbols in the various drawings indicate like elements.
DETAILED DESCRIPTION
[0039] Disclosed herein are techniques, systems, methods, devices, and other implementations, for estimating a maximum power point of a photovoltaic (PV) array. For example, in some embodiments, a solar power system includes a device that converts an electrical direct current produced in a PV array into an alternating current. The alternating current is then provided as an input into a computer-based system, which is configured, among other things, to track the MPP to enable efficiently harvesting power from the PV array. The computer-based system may further include a clock that is utilized to time intervals of the signal sampling period, intervals of the perturb cycle period, a perturb point, and an observe point. The alternating current is sampled at intervals, which are less than the interval of the perturb cycle, in order to generate a predictive model of power behavior of the array from which a trajectory of the predicted behavior of the PV array is generated. When the actual voltage of the array is perturbed, the difference between the actual power of the PV array and the predicted power at a subsequent time instant is determined. The MPP can be properly tracked in dynamic insolation conditions by allowing the direction of the next perturbation cycle to be correctly chosen based on the determined power difference. Other system implementations are possible. [0040] In some embodiments, methods, systems, and other implementations are provided, including a method for power tracking that includes measuring a plurality of samples of power produced by a photovoltaic (PV) array over a first interval of time, determining, based on the measured plurality of samples of the power of the PV array, a non-linear predictive model of a behavior of the power produced by the PV array, and performing a first adjustment of the PV array's voltage at a second time instant subsequent to an end of the first interval of time. The method further includes measuring another sample of the power of the PV array with the adjusted voltage, and determining a power difference between the power of the PV array with the adjusted voltage at a third time instant and a computed power level of the PV array at the third time instant determined from the non-linear predictive model.
[0041] With reference to FIG. 1A, a schematic diagram of a solar panel system 100 to harvest solar energy is shown. The system 100 includes a plurality of photovoltaic (PV) modules (or cells), such as the PV module 104, arranged as a PV array 103. In some implementations, each PV module may comprise a connected series of photovoltaic cells. Surfaces of the PV cells, which, in some implementations, may be configured similarly to photodiodes, are irradiated with radiation 101 from the sun 102 to excite electrons and create a potential between two electrodes. In some embodiments, the energy produced by the cells creates a direct current (DC). To convert this DC electricity into usable power (e.g., in the form of alternating current), the array 103 is connected to an inverter 1 10 (also referred to as a converter) configured to convert DC to AC.
[0042] Solar power may be measured in insolation, which is the amount of power (kilowatt hours) collected over a known surface area (square meters) and over a given time period (day). The insolation changes with changes to environmental conditions, such as temperature, cloud cover and position of the sun, and thus, the power produced by the PV modules fluctuates (generally in a non-linear manner) when changes to the environment in which the PV array 103 is deployed occur. In some embodiments, the PV array 103 may be mounted on a moveable support 105, which can be adjusted to position the planar surface of the array at an angle with the sun that seeks to maximize the radiation received from the sun.
[0043] In order to extract energy from the solar panels according to the load demand, in some embodiments, techniques and procedures to identify optimal, or near optimal, operating points of the PV array are performed. For example, as load and environmental conditions affecting performance of the PV array change, the optimal voltage or current of the PV array varies. However, because PV arrays are generally configured to operate at a constant operating point (e.g., at a constant voltage), fluctuating load and/or environmental conditions may render operation of the PV array, configured to operate at a set voltage or current, as non-optimal. Accordingly, in some implementations, as load and environmental conditions change, the operating points of the PV array may be adjusted so as to identify the optimal (or near optimal) operating point for the PV array for the present load and/or environmental conditions.
[0044] In some embodiments, techniques/procedures such as a maximum power point tracking (MPPT) procedure may be utilized. The MPPT procedure is configured to determine, based on the supply and demand of the PV array and the load, an optimal voltage (when the PV array is configured to operate with a constant voltage and a variable current that varies with changing environmental and/or load conditions), and/or the optimal current (when the PV array is configured to operate with a constant current and a variable voltage level that varies with changing environmental and/or load conditions) to enable optimal, e.g., maximum or near maximum, power to be provided by the PV array. For example, and with reference to FIG. IB showing an example graph 150 of power output behavior for a PV array as a function of the array's voltage, the point where the difference between the supply and demand is zero, or approximately zero (at a point 152 on the graph 150) corresponds to the maximum power point (MPP) of the PV array where the PV array is supplying power to the load substantially efficiently, and little or no additional power is being lost or required. The point on the graph at which maximum power is delivered corresponds to a point where the derivative of the graph, dP/dV, is equal to 0. At other points on the graph, the derivative dP/dV (as V increases) is not equal to 0, and the PV array is thus not operating at an optimal state. For example, a section 154 of the graph 150, corresponds to a derivative dP/dV being smaller than 0, and may be indicative that voltage levels are too high for given load conditions. In another example, points on the graph 150 where the derivative dP/dV is greater than 0 (dP/dV > 0) as V increases, as depicted in the section 156 of the graph, may represent voltage levels that are too low for existing load conditions.
[0045] When a PV array is not operating at optimal or near optimal conditions, e.g., when the derivative dP/dV of the power vs. voltage curve is non-zero, a procedure, such as the Maximum Power Point Tracking (MPPT) method, may be implemented to iteratively adjust the voltage level of the PV array until optimal or near optimal power delivery by the PV array for the given load and/or environmental conditions is established. [0046] For example, in some embodiments, an MPPT-based procedure may be implemented as a "perturb and observe" (P&0)-based procedure, in which the voltage of the PV array is periodically perturbed (e.g., the voltage produced by the PV array is adjusted), and the change in the current, and thus in the power, provided by the PV array as a result of the change in voltage, is observed. The difference between the observed power value and an earlier power value can thus provide a value representative of dP/dV. This procedure can be described according to the following equation:
Ppv (T) - PpV (0) = x (1) where, Ppy is the observed power of the photovoltaic array, T is the period of the perturb interval (in this example), and x is the difference in power over one period (T). If the change in the current (or power) is zero, i.e., x = 0, the MPP is being realized. However, if the change in the power is such that the power increased or decreased, the MPP has shifted in accordance with the observed change. For example, if dP/dV > 0, it can be concluded that the perturbation caused the array's operating power to move towards MPP, whereas if dP/dV < 0, it can be concluded that the perturbation caused the array's operating power to move away from the MPP.
[0047] Accordingly, to implement MPPT procedures, the system 100 depicted in FIG. 1A may include a controller 120 (e.g., a processor-based controller, a DSP processor, an application-specific integrated circuit, and/or any combination thereof) configured to determine whether non-optimal conditions of the systems and devices exist (e.g., by tracking the behavior of the power produced as a function of changing voltage and/or current), and to controllably adjust the voltage/current to converge to an optimal or near optimal power behavior for the PV array. For example, in situations where a controller coupled to a PV array indicates non-optimal power conditions, e.g., by determining that dP/dV is non-zero, the controller is configured to cause an adjustment of the voltage of PV array, and to repeat these operations if the voltage adjustment did not result in establishment of optimal or near- optimal conditions (e.g., MPP was not established). As shown in FIG. 1A, the controller 120 may be electrically coupled to, for example, the output of the inverter 110 (e.g., to enable determining power values) and/or to the PV array 103 (e.g., to control the array's voltage and/or current).
[0048] MPPT Procedures implemented in accordance with, for example, Equation 1 above, measure the power (e.g., multiplying the array's voltage by the array's current) at one time instant (e.g., when the system was perturbed by, for example, causing the voltage of the array to increase or decrease), and at a subsequent time instant, T (the time period between the time instants should be sufficiently long to enable the power behavior of the array to settle down).
[0049] The computed power difference according to Equation 1 does not indicate what contributed to the power measured at the subsequent time instant. Particularly, if during the period of measurement environmental and load conditions were unchanged, then it can be assumed that the power measured at time instant T resulted substantially from the adjustment made to the voltage of the PV array. However, if there were changes to the environmental conditions in which the PV array is operating (temperature, cloud coverage, sun brightness, etc.) then these environmental changes may have contributed, at least to some extent, to the measured change in the power supplied by the PV array (for the purpose of explanation, it will be assumed that load conditions are not changed during the PV array power measurement period). It is therefore difficult to determine if the "power delta" sign, based on which the direction of tracking is determined, was in response to either the voltage perturbation, insolation changes, or both. Accordingly, under these circumstances, a computation of power difference in accordance with Equation 1 may provide a value that would not be an accurate indication of whether further adjustments to the voltage of the PV array are necessary, and/or in what direction (e.g., voltage increase, voltage decrease) the adjustment should be made. Consequently, power optimization procedures, such as the MPPT procedure, may be tracking in the wrong direction.
[0050] Therefore, in some embodiments, power optimization procedures for PV arrays that take into account the effect/contribution of insolation changes may be implemented. Particularly, in some implementations, computation of the power difference may be performed in accordance with a model that incorporates the effect/contribution of insolation changes to changes in the power supplied/produced by a PV array. With reference to FIG. 2, a graph 200 of the power behavior of a PV array over a period of time as a result of insolation changes and voltage perturbation is shown. As illustrated in FIG. 2, during the period 0 < t < Tl (i.e., the period preceding the instant of perturbation), the power is shown to vary due to, in the example of the graph 200, insolation changes (e.g., changes in the temperature, cloud coverage, etc.). During that period, also referred to as the "identification period" or "sampling period," the power behavior of the PV array resulting from insolation changes may be sampled, and based on the sample points recorded, a model representative of power behavior as a result of insolation changes is determined. The derived model can therefore be used to predict the power behavior of the PV array, as a result of insolation changes, at a time subsequent to the end of the sampling period. FIG. 2 shows a dashed projection 204 corresponding to the power behavior of the PV array predicted from the model derived from the samples collected during the period 0 < t < Tl . The dashed projection 204 is a predicted extension of sampled portion of the graph 200 prior to its perturbation, and therefore does not represent the behavior of the power level of the PV array resulting from perturbation of the PV array.
[0051] As further shown in FIG. 2, at a time subsequent to the sampling period, the PV array system is perturbed (e.g., its voltage level is altered), resulting in power behavior represented, in the example of FIG. 2, by the solid line portion 210 of the graph 200. The solid line portion 210 represents the behavior of the power of the PV array resulting from the effect of both the perturbation and the insolation changes. Generally, the model derived from the samples can provide reasonable short term predictions of the power behavior from insolation changes (or any other factor which is to be modeled), but the accuracy and validity of the predicted behavior diminishes over time. As can be seen from FIG. 2, the effect the perturbation has on the power behavior of the PV array can be determined by identifying, using the derived model (represented by the dashed portion 204 depicted in FIG. 2), the effect insolation changes are predicted to have on the power behavior, and subtracting from the overall power behavior (represented by the solid line section 210 of the graph 200) the predicted effect/contribution of the insolation changes on the power behavior. Accordingly, in some embodiments, the effect of perturbation can be represented as:
PPV(T) - PpV_Model(T) (2) where T is a time instant, subsequent to the end of the sampling period and to the time instant when the PV array was perturbed, when the power supplied by the PV array is sampled. The difference between the actual power supplied by PV array (such actual power being affected by, for example, the perturbation and insolation changes) and the power level predicted using the mathematical model modeling the power behavior of the PV array as a result of insolation changes, provides a value representative of effect of perturbation on the power supplied by the PV array. The resultant value will therefore provide a more accurate indication of dP/dV, and will thus provide a more accurate and reliable tracking of the power behavior of the PV array so as to more quickly achieve an optimal, or near optimal, voltage level (or current level, in situations where the PV array has a near constant current and fluctuating voltage).
[0052] As will be described in greater details below, in some implementations, the model derived to predict the power behavior of the PV array as a result of insolation changes (at least for short term prediction) may be a non-linear model, e.g., high order models such as, for example, second and/or third order models. Use of non-linear models provide a more accurate and robust representation (e.g., as compared to representations provided, for example, by linear models) of a rate of change (in this case, the PV array's power rate of change), and therefore may enable more accurate tracking of the PV array's MPP. Additionally, use of a non-linear model enables non-constant sampling (i.e., using varying sampling intervals in which the sampling intervals between successive samples are not necessarily the same) to obtain the samples that are used to derive the non-linear representation corresponding to the samples. In some embodiments, and as will be described below in greater detail, derivation of non-linear predictive models to represent the behavior of PV array power resulting from insolation changes (or resulting from other factors) may be based on one or more of, for example, interpolation polynomials, a least square technique, and/or trigonometric interpolation.
[0053] With reference to FIG. 3, a flow chart of a power tracking procedure 300 that takes into account the effect of insolation changes on the power tracking functionality performed by, for example, the system 100, is shown. Implementation of the procedure 300 may be based, for example, on configuring the controller 120 depicted in FIG. 1 to perform the operations described herein. As noted, to determine the effect insolation changes (or other factors) have on the power behavior of a PV array, such as the PV array 103 depicted in FIG. 1 , a plurality of samples of power produced by the PV array over a first interval of time is measured 310. Sample measurements can be performed, for example, by separately measuring the voltage and current produced by the PV array (e.g., using voltmeter, ammeters, etc.) and multiplying the two to obtain an instant power value supplied by the PV array. In some embodiments, the measured electrical attributes may be performed with respect to the output of the inverter (e.g., the inverter 1 10) of the PV array. Measuring power samples may also be performed, for example, by dividing the identification period (also called the "sampling period) into a number of intervals. In some embodiments, and as shown in FIG. 2, the identification period may be divided to substantially equal length intervals of duration Ts. As noted, when non-linear models are derived to represent the behavior of power supplied by the PV array, the sampling period does not have to be constant, and the intervals between successive samples may vary over the duration of the identification period. In some implementations, the measured samples may be obtained by using a sample-and-hold circuit to determine the voltage level and/or current level at the measurement instants, and using, for example, an analog-to-digital converter (ADC) to digitize the measured samples. Such sample-and-hold circuit(s) and/or an ADC may be part of the example controller 120 of FIG. 1.
[0054] Based on the measured plurality of samples of the power level of the PV array, a non-linear predictive model of a behavior of the power produced by the PV array may be determined 320. As noted, in some embodiments, a non-linear model based on interpolation polynomials may be derived. Various forms may be used to represent interpolating polynomials. For example, an interpolating polynomial function may be represented as a linear combination of orthogonal polynomials. In some implementations, an interpolating polynomial may be represented as combination of Lagrangian polynomials. Particularly, a combination of Lagrangian polynomials may be represented as:
L(x) =A 7=0∑yj ¾ (t) (3)
[0055] In order to determine the number of samples from a given set of, for example, n+1 data points at the intervals ηΊ at which the signal is sampled, the data points can be represented as j =/Ts and y = Ppv (fTs), where j = 0.... n. It is to be noted that although in this particular example the sampling interval is substantially constant, a non-constant (irregular) sampling interval may be used instead. Utilizing the aforementioned representation in equation (3), the polynomials can be expanded as the following mathematical function: lj (t) = (0< i < n, i≠j) Π (t-x;) / (xj -xd =
(t-X0)/(x/- Xo)...(t-X/_l)/(Xj -X/_l) · (t-X/+1)/(Xj -x/+i). · · (4)
[0056] When the system is running on a clocked time interval, it is assumed that the time interval at which the signal is sampled is a unity time interval, e.g., Ts = 1 , 2, 3, etc. Equation (4) can then be reduced to: tj (t) = (0≤i < n, i≠j) U (t- / (/- =
(t-0)/(/--0)... (t-(/-- 1 ))/(/-(/·- 1 )) · (t-(/+l))/(/ -(/+; )· · . (5)
[0057] Equation (5) can thus be utilized to compute the sampled signal power at any number of intervals "«" to collect data points which predict the trajectory of the PV power signal. The integration of the Equation (5) provides the predictive PV power signal function Ppv modei(t) having n+1 samples over a period Ts represented as: Ppv_modei ((n+l)Ts) | =;) = L(n+1) = (6) y0£o(n+l) + y^n+l).... y(n)(n+l)£(n) (n+1)
In Equation (6), each€(n) can be derived from equation (5), with each selection of "n" producing "n+1" data points.
[0058] As an example, consider a situation where n = 2. In this case, there are 3 sampling points, giving rise to the following expressions: io(t) = (t-l)/(0-l) (t-2)/(0-2) = (t-l)(t-2)/2
£i(t) = (t-Oy(l-O) (t-2)/(l-2) = -t(t-2) l2(t) = (t-0)/(2-0) · (t-l)/(2-l) = t(t-l)/2
Ppv_modd(t) = y0£o(t) + yi*i(t) + y22(t)
[0059] With the above expression for Ppv modei(t), the value of the function for t = (n+l)Ts may be computed as follows:
Ppv_modei((n+l)Ts) = (n+\) \n=2 = L(3)
L(3) = y0€o(3) + yi€i(3) + y2£2(3) = y0(3-l)(3-2)/2 - 3yi(3-2) + y2(3-0)(3-l)/2 = yo - 3yi + 3y2
Therefore, PPV modei(3Ts) = PPV(0) - 3PPV(TS) + 3PPV(2TS)
[0060] FIG. 6 is a schematic diagram of an example DSP implementation 600 to determine the power difference between the power (or a discretized value thereof) measured following perturbation (i.e., P(k)) and the value computed using a second order Lagrange polynomial used to predict the power behavior of a PV array resulting, for example, from insolation changes. The power differential can be expressed as:
P(k) - PPV_model(k),
Substituting the expression for the 2nd order Lagrangian PPV model derived above provides:
P(k) - [P(k-1) - 3P(k-2/3) + 3P(k-l/3)], where the samples are a discretized distance of 1/3 from each other. Accordingly: dP(k) = P(k) - P(k-1) + 3P(k-2/3) - 3P(k-l/3)
[0061] Thus, to obtain output value 620, corresponding to dP(k), the implementation 600 includes a summing operator 610 that adds the sample 602, corresponding to P(k), to the negative of the preceding sample 604 (corresponding to P(k-l/3) that was amplified by an amplifier 612 with a gain of 3. The resultant sum of the summing operation at 610 is added, using the summing operator 616, to a resultant amplified valued of the sample 606, corresponding to P(k-2/3). The sample 606 is amplified using an amplifier 614 that also has a gain of 3. Finally, the result of the summing operator 616 is added to the negative of the sample 608, corresponding to P(k-1). It is to be noted that the implementation of FIG. 6 is just one example of a realization a predictive model and for determining the power difference according to, in this example, dP(k) = P(k) - P(k-1) + 3P(k-2/3) - 3P(k-l/3), and that other implementations may be used instead of, or in addition to, the DSP-based realization depicted in FIG. 6.
[0062] As another example, consider the situation where n = 3. In this case, there are 4 sampling points, giving rise to the following expressions:
*o(t) = (t-l)/(0-l) (t-2)/(0-2) · (t-3)/(0-3) = (t-l)(t-2)(t-3)/(-6)
£ι(ί) = (t-0)/(l-0) (t-2)/(l-2) · (t-3)/(l-3) = t(t-2)(t-3)/2 i2(t) = (t-0)/(2-0) (t-l)/(2-l) · (t-3)/(2-3) = t(t-l)(t-3)/(-2)
£3(t) = (t-0)/(3-0) (t-l)/(3-l) · (t-2)/(3-2) = t(t-l)(t-2)/6
[0063] With above expressions, Ppv modei(t) for t = (n+l)Ts may be computed as follows:
PpV_model((n+l)Ts) |„=3 = L(4)
L(4) = y0(4-l)(4-2)(4-3)/(-6) + yi4(4-2)(4-3)/2 + y24(4-l)(4-3)/(-2)
+ y34(4-l)(4-2)/6
= -yo + 4yi -6y2 + 4y3
[0064] Therefore, the predicted power level at 4Ts, as computed based on the above derived expression, is: Ppv_modei(4Ts) = -PPV(0) + 4PPV(TS) - 6PPV(2TS) + 4PpV(3Ts)
[0065] In some embodiments, derivation of another type of non-linear predictive model to represent the power behavior due to, for example, insolation changes, may be based on the Least Square Method. Using this technique, the resulting signal model is generally a polynomial function with a function order of m where m < n +1 . The parameter m may be selected based on the desired modeling accuracy. The implementation of the least square method to facilitate derivation of a model based on the samples collected can be, in some embodiments, in the form of recursive calculations.
[0066] Derivation of a predictive model based on a least squares procedure may be implemented as follows. Given a set of n+l data points, i.e., (xo,yo), · · · , (xn,yn), where xj = jTs and yj- = PpvOTs), and where j = 0 . . . n, the objective is to determine coefficients of the following polynomial:
Ppv_modei(t) = y(t) = a0 +ait + a2t2 + ... + amtm , such that the criterion function φ is minimized, where, for example, φ may be: φ =∑ni=o [yi - (a0 +aiXj + a2x; 2 + ... + Ά„Χ )]2
[0067] The coefficients obtained as a solution of the following system of linear equations: na0 + ai∑ i=o X; + a2∑ i=o x i + ... + am1=0 x i -∑ i=o y
„ i=0 Xi + a!∑ i=o x i + a2∑ i=0 X; + . .. + ami=0 Xi -∑ i=o i¼
„ i=o x i + ai∑ i=0 Xj + a21=o Xi + ... + am∑ i=o x; -∑ i=o x; y
[0068] As an example, consider the situation where n=2, in which case there are 3 data sets, namely, (xo,yo), (xi,yi), and (x2,y2). Let Ts =1 , in which case xo = 0, xi = 1 , and x2 = 2. In this example, m =1 , which means that in this case, when a statistical procedure, such as the least square procedure, is used to derive the predictive model, the predictive model may be represented as a first order (linear) function. Thus, PVPV model can be represented as:
PVpv model = ao + &\t where ao and ai are obtained from:
3a0 + a^xo + xi + x2) = y0 + yi + y2. ao (xO + xl + x2) + ai(xo2 + xi2 + x2 2) = xoyo + xiyi + x2y2
[0069] When substituting the values of xo = 0, xi = I , and x2 = 2, the following set of equations are obtained:
3a0 + 3ai = y0 + yi + y2; and
3ao + 5ai = yi + 2y2 and ao and ai can be determined to be: ai = (Υ2 - yo)/2, and ao = 5/6 y0 + 1/3 yi - l/6y2 [0070] PVpV model may therefore be expressed as: pvpv model (yi - yo) 1 12 + 5/6 yo + 1/3 yi - l/6y2 [0071] Thus, at the time t = (n+l)Ts |n=2,Ts=i = 3, PVpv model is determined to be:
PVpv model ((n+l)Ts |n=2,Ts=l) = PVpv model (3)
= 3/2 (y2 - yo) + 5/6 y0 + 1/3 yi - l/6y2
= -2/3 yo + 1/3 yi + 4/3 y2
[0072] In yet another example, in some embodiments, derivation of another type of nonlinear predictive model to represent the power behavior due to, for example, insolation changes, may be based on trigonometric interpolation. This derivation technique is suitable, for example, in situations where the model is to be represented as a periodic function. Consider, for example, the use of trigonometric interpolation for n+l data points, namely (xo,yo), · · · , (xn,y«). A function F(x) may be represented as follows:
F(x) = yo [sin ½ (x-xi) sin ½ (x-x2) ... sin ½ (x-xn)] / [sin ½ (xo-xi) sin ½ (xo-x2) .... sin ½ (x0-xn)] + yi [sin ½ (x-x0) sin ½ (x-x2) ... sin ½ (x-xn)] / [sin ½ (xi-xo) sin ½ (xi-x2) .... sin ½ (xi-xn)] + + yn [sin ½ (x-xo) sin ½ (x-x2) sin ½ (x-xn_i)] / [sin ½ (xn-xo) sin ½ (xn-X2) sin ½ (xn-xn- i)]
[0073] As an example, consider the situation where n = 3. In this case, there are 4 sampling points which, for the purpose of illustration, are sampled at Xj = jTs, where j = 0, 1, 2, 3. With four (4) sample points, the predictive model, PPV _modei(t may be expressed as follows:
PpV_model(t) = F(t) = yo [sin ½ (t-1) sin ½ (t-2) sin ½ (t-3)] / [sin ½ (0-1) sin ½ (0-2) sin ½ (0-3)] + yi [sin ½ (t-0) sin ½ (t-2) sin ½ (t-3)] / [sin ½ (1-0) sin ½ (1-2) sin ½ (1-3)] + y2 [sin ½ (t-0) sin ½ (t-1) sin ½ (t-3)] / [sin ½ (2-0) sin ½ (2-1) sin ½ (2-3)] + y3 [sin ½ (t-0) sin ½ (t-1) sin ½ (t-2)] / [sin ½ (3-0) sin ½ (3-1) sin ½ (3-2)]
[0074] For ease of illustration, in the above model derivation examples, the sampling points were assumed to be taken at equal duration intervals. However, as noted, derivations to obtain predictive model representations for the power behavior may be performed with sample points taken at non-constant (e.g., irregular) intervals. Additionally, although three different derivation techniques were illustrated, other techniques/procedures to determine non-linear predictive model to represent the power behavior of the PV array may be used in addition to or instead of the any of the above illustrated techniques/procedures.
[0075] Turning back to FIG. 3, after measuring/obtaining the samples, and determining, based on the samples obtained, a non-linear model to represent the power behavior of the PV array, tracking operations to determine the maximum power point using the derived nonlinear model may be performed. Thus, a first adjustment is performed 330 of the PV array's voltage at a second time instant subsequent to an end of the first interval of time (during which the power samples were measured/obtained). For example, the voltage adjustment performed, i.e., perturbation of the system, may be to increase or decrease the voltage of the PV array by some pre-determined amount, or by an amount that depends on the present voltage level (as may have been measured, for example, when the last sample was obtained). In some embodiments, the determination of whether to initially increase or decrease the voltage (i.e., before any power difference computation is performed) may be based on some pre-determined set of rules (e.g., that the voltage is to initially be increased, that a decision as to initially increase or decrease the voltage may be based on the present voltage level or on the difference between power levels of the last two samples obtained, etc.) In the example of FIG. 2, voltage adjustment, or perturbation, occurs at the time instant Tl = nTs (the time instant Tl is also referred to as the "point of perturbation"). As noted, in some embodiments, the PV array may be configured to output a substantially constant current, and to have varying voltage that depends, for example, on load changes, insolation changes, etc. Under those circumstances, perturbation of the system may be performed by making a current adjustment to the PV array, and tracking the maximum power point resulting from such current perturbations.
[0076] At a third time instant T, where T > Tl, another sample of the power of the PV array with an adjusted electrical attribute (e.g., adjusted voltage, in the example of FIG. 2) is measured 340 (i.e., the observation of the system following the perturbation of the system). The measurement of the power level of the PV array with the adjusted electrical attribute may be performed using a sample-and-hold circuit to determine the voltage level and/or current level at the instant of measurements, and by using, for example, an analog-to-digital converter to digitize the measured sample. Generally, the time duration between "Tl" and "T" should be sufficiently long to enable the system to respond to the adjustment performed.
[0077] Having obtained the other sample, a power difference is determined 350 between the power of the PV array with the adjusted voltage at the third time instant and a computed power level of the PV array at the third time instant determined from using the non-linear predictive model. Thus, in some implementations, a determination is made of Ppv(T) - Ppv modd(T). The power difference therefore provides a value indicative of the effect the perturbation had on the power being tracked by attempting to remove the contribution that insolation changes (computed using the determined non-linear model), or any other factor represented by this or some other non-linear predictive model, had on the power being tracked.
[0078] The determined power difference may be utilized to track the MPP of the PV array and control the output of the PV array in relation to the load. For example, the determined power difference can be utilized to determine the next direction of perturbation (e.g., adjustment of an attribute of the PV array, be it the voltage or current) based on the sign of a result of the power difference, and to use the resultant sign (positive "+", or negative "-") to control the direction of the perturbation during the next cycle. Thus, in some embodiments, the procedures implemented to track the maximum power point of a PV array may include performing another adjustment of the voltage of the PV array based on the determined power difference between the power of the PV array with the adjusted voltage at the third time instant and the computed power level of the PV array at the third time instant as determined from the non-linear predictive model.
[0079] For example, in response to a determination that the power difference is greater than ε, where ε is a pre-determined value that is generally smaller than the power level of the PV array with the adjusted voltage (ε may be equal to 0 in some embodiments), if the preceding adjustment of the voltage (or another electrical attribute) of the PV array was a voltage increase then the next voltage adjustment for the PV array would be to increase the voltage. Particularly, if the computed power difference is determined to be greater than ε (e.g., greater than 0), this may indicate that the power is being tracked in the right direction. Therefore, if the previous adjustment was to increase the voltage, then the next adjustment should also be a voltage increase (e.g., by some pre-determined amount, by an amount that is based on, for example, the present voltage or power level, etc.) If the previous adjustment was a voltage decrease, then the next adjustment should also be a voltage decrease.
[0080] Conversely, if it is determined that the power difference is less than ε, and the preceding adjustment of the PV array was a voltage increase, then the next adjustment of the PV array's electrical attribute(s) may be to decrease the voltage. Particularly, because a computed power difference that is less than ε (e.g., less than 0) may indicate that the maximum power point is being tracked in the wrong direction (i.e., moving away from the maximum power point), if the previous adjustment was to increase the voltage of the PV array, then the voltage may have to be decreased on the next adjustment so as to reverse to tracking direction. Similarly, if it is determined that the power difference is less than ε, and the previous adjustment was a voltage decrease, then the next adjustment may be a voltage increase.
[0081] In some embodiments, the process of adjusting the electrical attribute(s) of the PV array may be repeated until, for example, an absolute value of the determined power difference is less than some pre-determined value, indicating that the PV array may have substantially or approximately reached its maximum power point.
[0082] While the implementations discussed above refer to the perturb and observe technique, similar implementations may also be realized in relation to other types of tracking techniques, including, for example, the incremental conductance procedure, the constant voltage procedure, etc. For example, in implementations based on incremental conductance procedure, the incremental conductance dl/dV is determined so as to compute the sign of dP/dV. With this procedure, the system is determined to have reached its MPP when dl/dV = -i/v. If the system has not reached its MPP yet, the system continues to be perturbed.
[0083] With reference to FIG. 4, a schematic diagram of an example system implementation 400 to perform a tracking procedure, such as that performed by the procedure 300 of FIG. 3, is shown. In some embodiments, the system 400 may be implemented, at least in part, as the controller 120 of FIG. 1A. The system 400 may include, for example, a sampling and PV power calculation module 410 (also referred to as a signal identification unit/module) to measure one or more samples of a power level produced by a PV array (such as the PV array 103), and may include sample-and-hold circuits and analog-to-digital converters. The Data Sampling and PV Power Calculation module 410 samples the voltage and current of the PV generator, computes power and provides the input data samples to the Power Increment Calculation module 420. The input to the Power Increment Calculation module 420 may be implemented, in some embodiments, to receive three (3) samples of the power signal taken before a perturbation (e.g., P(k-l/3), P(k-2/3) and P(k-l)), and one sample of the power signal taken after the perturbation P(k). Other implementations (e.g., using different sample sequences) may be used. The Power Increment Calculation module 420 of the example system 400 uses the "before perturbation samples," e.g., P(k-l/3), P(k-2/3) and P(k-l), to construct, for example, a non- linear predictive model, and to determine a predicted (expected) power after the perturbation model, i.e., Pmodei(kT). As described herein, the predictive model may be derived based on interpolation polynomials (e.g., to derive a Lagrangian polynomial based model), a least square technique, a trigonometric interpolation technique, etc. Having computed a predicted power level using the derived model, a power difference may be computed, e.g., based on dP(kT) = P(kT) - Pmodei(kT).
[0084] The example system 400 also includes a Perturb Direction Computation module 430 which maintains a perturbation history (e.g., previous perturbation direction). The module determines, based on the power difference dP(kT) and the previous perturbation direction, what the next perturbation direction should be. For example if the previous PV voltage perturbation was in a positive direction (direction(k-l)>0) and the power difference (also referred to as "power increment") is positive (i.e., dP(k)>0) then the next perturbation will be in a positive direction too. However if the power difference is negative (i.e., dP(k)<0) than the perturbation direction at t=kT should have the opposite sign to the direction in t=(k- 1)T. [0085] As further shown in FIG. 4, the example system 400 also includes a PV Voltage Reference Calculation module 440 configured to receive the direction command (e.g., from a module/unit such as the Perturb Direction Computation module 430 of the system 400) and to compute a PV voltage reference as follows:
If direction(k)>0 than: Vref(k) = Vref(k-1) + Vperturbe step
If direction(k)<0 than: Vref(k) = Vref(k-1) - Vperturbe step where Vperturbe step represents a perturbation step voltage, which may be a pre-determined step voltage, or may be an adjustable dynamically derived value. In some embodiments, Vref(k) may also be compared to an associated minimum and maximum limits and clipped to those limiting values. The minimum and maximum limits may be based on the PV generator, Inverter, and/or power converter specifications. The computed Vref(k) is provided to an inverter module configured, among other things, to adjust the current of the PV generator such that the PV voltage converges to the reference voltage level.
[0086] Performing the tracking and control procedures described herein may be facilitated, in some implementations, by a processor-based computing system which, for example, may be used to implement, at least partly, the controller 120 of FIG. 1. With reference to FIG. 5, a schematic diagram of a generic computing system 500 is shown. The computing system 500 includes a processor-based device 10 such as a personal computer, a specialized computing device, and so forth, that typically includes a central processor unit 512. In addition to the CPU 512, the system may include a main memory, cache memory and bus interface circuits (not shown). The processor-based device 510 may include a mass storage element 514, such as a hard drive and/or a flash drive associated with the computer system. The computing system 500 may further include a keyboard, or keypad, 516, and a monitor 520, e.g., a CRT (cathode ray tube) or LCD (liquid crystal display) monitor, that may be placed where a user can access them.
[0087] The processor-based device 510 is configured to facilitate, for example, the implementation of the tracking procedure, including facilitating measurement of sample points of the power output generated using a PV array, determining a non-linear predictive model representative of the power behavior of the PV array system resulting from such factors as insolation changes, determining whether power tracking is proceeding in the right direction, etc. The storage device 514 may thus include a computer program product that when executed on the processor-based device 510 causes the processor-based device to perform operations to facilitate the implementation of the above-described procedures. The processor-based device may further include peripheral devices to enable input/output functionality. Such peripheral devices may include, for example, a CD-ROM drive and/or flash drive, a network connection (implemented using, for example, a USB port, a wireless transceiver, etc.), for downloading related content to the connected system. Such peripheral devices may also be used for downloading software containing computer instructions to enable general operation of the respective system/device. Alternatively and/or additionally, in some embodiments, special purpose logic circuitry, e.g., an FPGA (field programmable gate array), a DSP processor, or an ASIC (application-specific integrated circuit) may be used in the implementation of the system 500. For example, in some implementations, a DSP processor may be used to implement the non-linear predictive model used to compute the predicted power behavior of the PV array due to, for example, insolation changes (or some other factor). An example of such an implementation is illustrated in FIG. 6 (as was more particularly described above). Other modules that may be included with the processor-based device 510 are speakers, a sound card, a pointing device, e.g., a mouse or a trackball, by which the user can provide input to the computing system 500. The processor-based device 510 may include an operating system.
[0088] Computer programs (also known as programs, software, software applications or code) include machine instructions for a programmable processor, and may be implemented in a high-level procedural and/or object-oriented programming language, and/or in
assembly/machine language. As used herein, the term "machine-readable medium" refers to any non-transitory computer program product, apparatus and/or device (e.g., magnetic discs, optical disks, memory, Programmable Logic Devices (PLDs)) used to provide machine instructions and/or data to a programmable processor, including a non-transitory machine- readable medium that receives machine instructions as a machine-readable signal.
[0089] FIG. 7 is a graph illustrating the performance of a tracking procedure, such as the tracking procedure 300 described in relation to FIG. 3. PV generator transients are simulated to include fast transients such as, for example, moving clouds, and slow transients, such as, for example, thermal effects. The generator simulated can be one that can produce up to 340kW with a maximum power point voltage of around 520V. The upper waveforms (in the upper box) of FIG. 7 shows the behavior of Vmppt(t) (the smooth line) and actual Vpv(t) in Volts. Vmppt(t) is the "ideal" position of the maximum power point voltage and Vpv(t) is the achieved PV voltage. The time axis is provided in seconds. A tracking achieves good performance if the achieved PV voltage closely matches the maximum power point voltage. The lower waveform (lower box) of FIG. 7 shows the harvested power from the PV generator in Watts. The graph shows that the PV power has some large and relatively fast transients. Those transients may be the result of cloud movements. As can be seen from the top graph, in such dynamic conditions the difference between the maximum power point voltage and the actual PV voltage are within +/-5V.
[0090] Although solar energy is discussed in relation to the above implementations, the system, methods and all other implementations discussed above can also be used and applied in relation to other types of generators and for other forms of energy, such as energy harvested in wind generators and water pumps.
[0091] While particular embodiments have been disclosed herein in detail, this has been done by way of example for purposes of illustration only, and is not intended to be limiting with respect to the scope of the appended claims, which follow. In particular, it is contemplated that various substitutions, alterations, and modifications may be made without departing from the spirit and scope of the invention as defined by the claims. Other aspects, advantages, and modifications are considered to be within the scope of the following claims. The claims presented are representative of the embodiments and features disclosed herein. Other unclaimed embodiments and features are also contemplated. Accordingly, other embodiments are within the scope of the following claims.

Claims

WHAT IS CLAIMED IS :
1. A method for power tracking, the method comprising:
measuring a plurality of samples of power produced by a photovoltaic (PV) array over a first interval of time;
determining, based on the measured plurality of samples of the power of the PV array, a non-linear predictive model of a behavior of the power produced by the PV array;
performing a first adjustment of the PV array's voltage at a second time instant subsequent to an end of the first interval of time;
measuring at a third time instant another sample of the power of the PV array with the adjusted voltage; and
determining a power difference between the power of the PV array with the adjusted voltage at the third time instant and a computed power level of the PV array at the third time instant determined from the non-linear predictive model.
2. The method of claim 1 , further comprising:
performing another adjustment of the voltage of the PV array based on the determined power difference between the power of the PV array with the adjusted voltage at the third time instant and the computed power level of the PV array at the third time instant determined from the non-linear predictive model.
3. The method of claim 2, wherein performing the other adjustment of the voltage of the PV array comprises:
increasing the voltage of the PV array in response to a determination that the power difference is greater than ε when the first adjustment of the voltage of the PV array was a voltage increase, wherein ε is a pre-determined value smaller than power produced by the PV array.
4. The method of claim 2, wherein performing the other adjustment of the voltage of the PV array comprises:
decreasing the voltage of the PV array in response to a determination that the power difference is less than e when the first adjustment of the voltage of the PV array was a voltage increase, wherein ε is a pre-determined value smaller than power produced by the PV array.
5. The method of claim 1 , wherein determining the non-linear predictive model of the behavior of the power produced by the PV array comprises:
deriving the non-linear predictive model using one or more of: interpolation polynomials, a least square technique, and trigonometric interpolation.
6. The method of claim 1 , wherein measuring the plurality of samples the power of the array over the first interval of time comprises:
measuring the plurality of samples at irregular time periods during the first interval of time.
7. The method of claim 1 , wherein the non-linear predictive model of the behavior of the power of the PV array includes a non-linear predictive model of the behavior of the power of the PV array resulting from insolation changes.
8. The method of claim 7, wherein the insolation changes include one or more of: a temperature change, a cloud cover change, and a sun position change.
9. A system comprising:
one or more photovoltaic (PV) arrays configured to convert solar irradiation into current; a sampling subsystem to measure one or more samples of power produced by the one or more PV arrays; and
a controller configured to:
determine, based on a plurality of samples of the power of the one or more PV array measured over a first time interval, a non-linear predictive model of a behavior of the power produced by the one or more PV arrays,
cause a first adjustment of the one or more PV arrays' voltage at a second time instant subsequent to an end of the first time interval, and
determine a power difference between another sample of the power of the one or more PV arrays with the adjusted voltage measured at a third time instant and a computed power level of the one or more PV arrays at the third time instant determined from the non-linear predictive model.
10. The system of claim 9, wherein the controller is further configured to:
cause another adjustment of the voltage of the one or more PV arrays based on the determined power difference between the power of the one or more PV arrays with the adjusted voltage at the third time instant and the computed power level of the one or more PV arrays at the third time instant determined from the non-linear predictive model.
1 1. The system of claim 10, wherein the controller configured to cause the other adjustment of the voltage of the one or more PV arrays is configured to:
cause an increase of the voltage of the one or more PV arrays in response to a determination that the power difference is greater than ε when the first adjustment of the voltage of the one or more PV arrays was a voltage increase, wherein ε is a pre-determined value smaller than power produced by the one or more PV arrays.
12. The system of claim 10, wherein the controller configured to cause the other adjustment of the voltage of the one or more PV arrays is configured to:
cause a decrease of the voltage of the PV array in response to a determination that the power difference is less than ε when the first adjustment of the voltage of the one or more PV arrays was a voltage increase, wherein ε is a pre-determined value smaller than power produced by the one or more PV arrays.
13. The system of claim 9, wherein the controller configured to determine the nonlinear predictive model of the behavior of the power produced by the one or more PV arrays is configured to:
derive the non-linear predictive model using one or more of: interpolation polynomials, a least square technique, and trigonometric interpolation.
14. The system of claim 9, wherein the controller comprises one or more processors.
15. The system of claim 9, wherein the sampling subsystem configured to measure one or more samples of the power produced by the one or more PV arrays is configured to: measure the plurality of samples at irregular time periods during the first interval of time.
16. The system of claim 9, wherein the non-linear predictive model of the behavior of the power of the one or more PV arrays includes a non-linear predictive model of the behavior of the power of the one or more PV arrays resulting from insolation changes.
17. The system of claim 16, wherein the insolation changes include one or more of: a temperature change, a cloud cover change, and a sun position change.
18. A method for power tracking, the method comprising:
measuring a plurality of samples of power produced by a photovoltaic (PV) array over a first interval of time;
determining, based on the measured plurality of samples of the power of the PV array, a no n- linear predictive model of a behavior of the power of the PV array;
performing a first adjustment of the PV array's current at a second time instant subsequent to an end of the first interval of time;
measuring at a third time instant another sample of the power of the PV array with the adjusted current; and
determining a power difference between the power of the PV array with the adjusted current at the third time instant and a computed power level of the PV array at the third time instant determined from the non-linear predictive model.
19. The method of claim 18, further comprising:
performing another adjustment of the current of the PV array based on the determined power difference between the power of the PV array with the adjusted current at the third time instant and the computed power level of the PV array at the third time instant determined from the non-linear predictive model.
20. The method of claim 18, wherein the non-linear predictive model of the behavior of the power of the PV array includes a non-linear predictive model of the behavior of the power of the PV array resulting from insolation changes.
21. A method for power tracking, the method comprising: measuring a plurality of samples of power produced by a photovoltaic (PV) array over a first interval of time;
determining, based on the plurality of measured samples of the power of the PV array, a linear predictive model of a behavior of the power produced by the PV array, the linear predictive model derived using a least square technique;
performing a first adjustment of the PV array's voltage at a second time instant subsequent to an end of the first interval of time;
measuring at a third time instant another sample of the power of the PV array with the adjusted voltage; and
determining a power difference between the power of the PV array with the adjusted voltage at the third time instant and a computed power level of the PV array at the third time instant determined from the linear predictive model derived using the least square technique.
PCT/US2012/025739 2012-02-17 2012-02-17 Maximum power point tracking (mppt) WO2013122610A1 (en)

Priority Applications (5)

Application Number Priority Date Filing Date Title
CN201280072419.7A CN104220951B (en) 2012-02-17 2012-02-17 MPPT maximum power point tracking (MPPT)
EP12868625.0A EP2815286A4 (en) 2012-02-17 2012-02-17 Maximum power point tracking (mppt)
AU2012370020A AU2012370020A1 (en) 2012-02-17 2012-02-17 Maximum power point tracking (MPPT)
US14/379,030 US20150013748A1 (en) 2012-02-17 2012-02-17 Maximum power point tracking (mppt)
PCT/US2012/025739 WO2013122610A1 (en) 2012-02-17 2012-02-17 Maximum power point tracking (mppt)

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
PCT/US2012/025739 WO2013122610A1 (en) 2012-02-17 2012-02-17 Maximum power point tracking (mppt)

Publications (1)

Publication Number Publication Date
WO2013122610A1 true WO2013122610A1 (en) 2013-08-22

Family

ID=48984577

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/US2012/025739 WO2013122610A1 (en) 2012-02-17 2012-02-17 Maximum power point tracking (mppt)

Country Status (5)

Country Link
US (1) US20150013748A1 (en)
EP (1) EP2815286A4 (en)
CN (1) CN104220951B (en)
AU (1) AU2012370020A1 (en)
WO (1) WO2013122610A1 (en)

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103592992A (en) * 2013-11-18 2014-02-19 国家电网公司 System and method for rapid optimizing of maximum power point of photovoltaic array under shadow condition
CN104020814A (en) * 2014-06-26 2014-09-03 苏州青云能源科技有限公司 Near-linear maximum power point tracking method of photovoltaic power generation system
CN104035478A (en) * 2014-06-26 2014-09-10 苏州青云能源科技有限公司 Maximum power point tracing method for independent photovoltaic power generation system
WO2016040813A1 (en) * 2014-09-12 2016-03-17 Carnegie Mellon University Systems, methods, and software for planning, simulating, and operating electrical power systems
US20160248374A1 (en) * 2015-02-23 2016-08-25 Lsis Co., Ltd. Photovoltaic system
TWI553440B (en) * 2015-02-26 2016-10-11 國立中山大學 Maximum power point tracking method for photovoltaic generation
DE102015116455A1 (en) 2015-09-29 2017-03-30 Endress+Hauser Flowtec Ag Energy supply of measuring instruments from energy sources with poor performance
CN109814652A (en) * 2019-03-19 2019-05-28 武汉理工大学 A kind of thermoelectric power generation MPPT method based on non-equal part polychotomy

Families Citing this family (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR101965820B1 (en) * 2015-01-30 2019-04-04 엘에스산전 주식회사 Photovoltaic data collecting device
CN104950981B (en) * 2015-05-12 2017-03-01 江苏固德威电源科技股份有限公司 The multichannel MPPT method for tracing of Current Sensorless
CN105116957B (en) * 2015-07-22 2017-03-15 厦门科华恒盛股份有限公司 A kind of maximum power point of photovoltaic power generation system tracking
MY194129A (en) * 2015-11-19 2022-11-14 Sigmagen Inc Multi-modal maximum power point tracking optimzation solar photovoltaic system
CN105425894B (en) * 2015-12-01 2017-03-29 国网甘肃省电力公司电力科学研究院 A kind of photovoltaic system MPPT maximum power point tracking optimization method
CN106292829A (en) * 2016-08-31 2017-01-04 杭州康奋威科技股份有限公司 Solar energy photovoltaic system MPPT control method
US11101653B2 (en) * 2018-10-09 2021-08-24 S&C Electric Company Solar generation estimation

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20070290636A1 (en) * 2006-06-03 2007-12-20 Adensis Gmbh Photovoltaic system and method for operating a photovoltaic system
US20100117623A1 (en) * 2008-11-11 2010-05-13 Fife John M System and method of determining maximum power point tracking for a solar power inverter
US20100176771A1 (en) * 2009-01-15 2010-07-15 Xantrex Technology Inc. Perturb voltage as a decreasing non-linear function of converter power
US20110184583A1 (en) * 2010-01-22 2011-07-28 General Electric Company Model-based power estimation of photovoltaic power generation system
US20120013312A1 (en) * 2010-07-15 2012-01-19 Chih-Yuan Hsieh Power Control Device and Method thereof

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101651436B (en) * 2009-09-16 2013-04-03 阳光电源股份有限公司 High-precision maximum power point tracing method
US8455806B2 (en) * 2010-01-18 2013-06-04 Sunpower Corporation Photovoltaic assembly for use in diffuse weather conditions and related methods
JP5581965B2 (en) * 2010-01-19 2014-09-03 オムロン株式会社 MPPT controller, solar cell control device, photovoltaic power generation system, MPPT control program, and MPPT controller control method
CN102193561A (en) * 2011-04-17 2011-09-21 北京诚创科星科技有限公司 Maximum power point output tracking method for photovoltaic cell
CN102355003B (en) * 2011-09-13 2013-10-16 辽宁力迅风电控制系统有限公司 Control method and device for single-phase grid-connected photovoltaic power generation system

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20070290636A1 (en) * 2006-06-03 2007-12-20 Adensis Gmbh Photovoltaic system and method for operating a photovoltaic system
US20100117623A1 (en) * 2008-11-11 2010-05-13 Fife John M System and method of determining maximum power point tracking for a solar power inverter
US20100176771A1 (en) * 2009-01-15 2010-07-15 Xantrex Technology Inc. Perturb voltage as a decreasing non-linear function of converter power
US20110184583A1 (en) * 2010-01-22 2011-07-28 General Electric Company Model-based power estimation of photovoltaic power generation system
US20120013312A1 (en) * 2010-07-15 2012-01-19 Chih-Yuan Hsieh Power Control Device and Method thereof

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
ESRAM ET AL., COMPARISON OF PHOTOVOLTAIC ARRAY MAXIMUM POWER POINT TRACKING TECHNIQUE
R.ALONSO ET AL., A NEW METHOD SOLVING PROBLEMS WITH LOCAL MAXIMUM POWER PPOINTS
SALAS V ET AL., REVIEW OF THE MAXIMUM POWER POINT TRACKING ALGORITHMS FOR STAND-ALONE PHOTOVOLTAIC SYSTEMS
See also references of EP2815286A4 *

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103592992A (en) * 2013-11-18 2014-02-19 国家电网公司 System and method for rapid optimizing of maximum power point of photovoltaic array under shadow condition
CN103592992B (en) * 2013-11-18 2016-02-10 国家电网公司 The quick optimization system of maximum power point of photovoltaic array and method under a kind of shadowed condition
CN104020814A (en) * 2014-06-26 2014-09-03 苏州青云能源科技有限公司 Near-linear maximum power point tracking method of photovoltaic power generation system
CN104035478A (en) * 2014-06-26 2014-09-10 苏州青云能源科技有限公司 Maximum power point tracing method for independent photovoltaic power generation system
WO2016040813A1 (en) * 2014-09-12 2016-03-17 Carnegie Mellon University Systems, methods, and software for planning, simulating, and operating electrical power systems
US20160248374A1 (en) * 2015-02-23 2016-08-25 Lsis Co., Ltd. Photovoltaic system
TWI553440B (en) * 2015-02-26 2016-10-11 國立中山大學 Maximum power point tracking method for photovoltaic generation
DE102015116455A1 (en) 2015-09-29 2017-03-30 Endress+Hauser Flowtec Ag Energy supply of measuring instruments from energy sources with poor performance
CN109814652A (en) * 2019-03-19 2019-05-28 武汉理工大学 A kind of thermoelectric power generation MPPT method based on non-equal part polychotomy
CN109814652B (en) * 2019-03-19 2020-04-21 武汉理工大学 Thermoelectric power generation MPPT method based on non-equal multi-division method

Also Published As

Publication number Publication date
AU2012370020A1 (en) 2014-09-18
CN104220951A (en) 2014-12-17
CN104220951B (en) 2016-10-05
US20150013748A1 (en) 2015-01-15
EP2815286A1 (en) 2014-12-24
EP2815286A4 (en) 2015-12-02

Similar Documents

Publication Publication Date Title
WO2013122610A1 (en) Maximum power point tracking (mppt)
Fathabadi Novel solar powered electric vehicle charging station with the capability of vehicle-to-grid
Al Nabulsi et al. Efficiency optimization of a DSP-based standalone PV system using fuzzy logic and dual-MPPT control
US9354652B2 (en) Maximum power point tracking method and device, and photovoltaic power generation system
Li et al. An improved perturbation and observation maximum power point tracking algorithm based on a PV module four-parameter model for higher efficiency
CN104778352B (en) The seven parameter photovoltaic cell output characteristic modeling methods based on STFT operators
KR100989441B1 (en) Control system of solar cell generation using adaptive fuzzy model and method of tracking maximum power point using the same
Zhao et al. A novel maximum power point tracking strategy based on optimal voltage control for photovoltaic systems under variable environmental conditions
Tripathy et al. Photovoltaic system using Lambert W function-based technique
Yin et al. Multi-step depth model predictive control for photovoltaic power systems based on maximum power point tracking techniques
JP6287306B2 (en) Power estimation apparatus, power estimation method, and program
Khemliche et al. Bond graph modeling and optimization of photovoltaic pumping system: Simulation and experimental results
CN103729013B (en) The method and apparatus of tracking photovoltaic system maximum power point
CN113033136B (en) Simplified photovoltaic cell physical parameter extraction optimization method and system
CN106055019B (en) The photovoltaic generating system of MPPT maximum power point tracking based on online soft sensor model
Xu et al. Novel linear iteration maximum power point tracking algorithm for photovoltaic power generation
Raj et al. A two stage successive estimation based maximum power point tracking technique for photovoltaic modules
JP5903341B2 (en) Power generation control device, solar power generation system, and power generation control method
CN103995559B (en) A kind ofly determine voltage MPPT control method and system based on environment parameter model
Xiaoping et al. A novel maximum power point tracking in partially shaded PV systems using a hybrid method
Ammar et al. Evaluation of the performance of a FONN-based MPPT control for a photovoltaic watering system
KR101573277B1 (en) Solar generating apparatus having a function of dynamic mppt and the method thereof
CN111198592B (en) Photovoltaic maximum power point tracking method and device
CN105159388B (en) A kind of method of the MPPT maximum power point tracking in photovoltaic microgrid system
CN116048182A (en) Photovoltaic module flexible power point tracking method, device, equipment, medium and system

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 12868625

Country of ref document: EP

Kind code of ref document: A1

WWE Wipo information: entry into national phase

Ref document number: 14379030

Country of ref document: US

NENP Non-entry into the national phase

Ref country code: DE

WWE Wipo information: entry into national phase

Ref document number: 2012868625

Country of ref document: EP

ENP Entry into the national phase

Ref document number: 2012370020

Country of ref document: AU

Date of ref document: 20120217

Kind code of ref document: A