# Optimal Sizing of the Stand-Alone Photovoltaic System for a Solar-Powered Translational Sprinkler Irrigation Machine considering the Loss of Power Supply Probability

Authors: Kenan Liu; Bugong Sun; Xiaoyang Gao; Yang Zhang; Wei Sun; Quan Feng; Wanxia Yang
Journal: Mathematical Problems in Engineering
 (2022)

---

## Abstract

The translational sprinkler irrigation machine is widely used because of its high degree of automation, less manual investment, and convenient movement. However, when using the sprinkler irrigation machine, the power supply is difficult to be guaranteed in some remote power shortage areas. Solar energy has become one of the best choices for the power source of sprinkler irrigation machines. An important issue is the optimal sizing of the stand-alone photovoltaic system for a solar-powered translational sprinkler irrigation machine. This work conducts the optimal sizing of the stand-alone photovoltaic system for a solar-powered translational sprinkler irrigation machine considering the loss of power supply probability. Firstly, the self-developed translational sprinkler irrigation machine is introduced. The load power, which includes the driving power of the translational sprinkler irrigation machine, water intake pressure driving power, and the loss power of the controller and sensors, is considered in the calculation process. Subsequently, the photovoltaic generator model and the battery storage model are established. Then, the stand-alone photovoltaic system is optimized considering the loss of power supply probability (LPSP) and the life cycle cost (LCC). After the solar irradiance and ambient temperature of a typical sunny day in the test area are given, the optimal combination of the PV module and battery is obtained. The optimal sizing result is verified using the PV power, state of charge (SOC), and load power. The presented optimal sizing method is also compared with an existing method. Besides, the operation test of the sprinkler irrigation machine is also carried out to verify the optimal sizing result. The optimal sizing result is proved to be effective and applicable.


---

## Body

## 1. Introduction

The development of irrigation machinery is one of the guarantee factors to improve crop yield. How to save water and energy for irrigation has become an inevitable trend in modern agricultural development. The translational sprinkler irrigation machine is widely used because of its high degree of automation, less manual investment, and convenient movement. However, the power supply is difficult to guarantee when using the sprinkler irrigation machine in some remote power shortage areas. These power shortage areas cannot be irrigated in time. Solar energy is clean and environment-friendly energy. With the maturity of photovoltaic technology, solar energy has become one of the best choices for sprinkler irrigation machine power.The solar-powered translational sprinkler irrigation machine transforms the traditional farmland irrigation mode into mobile sprinkler irrigation with high uniformity of irrigation. An important issue is the optimal sizing of the stand-alone photovoltaic system for a solar-powered translational sprinkler irrigation machine. At present, there are fewer studies on the matching of power demand and solar photovoltaic power during the operation of solar-powered sprinkler irrigation machines. The sprinkler irrigation machinery powered by solar is not mature, which limits its popularization and application in agriculture. At present, the research on solar-powered machinery mostly focuses on the research of solar vehicles, such as [1–3].Many efforts have been made regarding the optimal sizing of the stand-alone photovoltaic system. Benmouiza et al. [4] presented an optimization method for designing the stand-alone photovoltaic system. This work mainly compares the proposed method with the existing methods. However, the specific application of the method is not illustrated. Aziz et al. [5] optimized the stand-alone photovoltaic system by minimizing the loss of power supply probability. In this study, only the loss of power supply probability (LPSP) was considered without considering the stand-alone photovoltaic system's life cycle cost (LCC). Fezai et al. [6] proposed a methodology to study the load profile impact on the stand-alone photovoltaic system reliability and storage component. Senol et al. [7] studied the technical and economic feasibility of photovoltaic pumping of water in Turkey. This study only considered the life cycle cost (LCC), ignoring the loss of power supply probability (LPSP). Sallem et al. [8, 9] presented a new management approach that makes the decision on the optimum connection times of the elements of a photovoltaic water pumping installation. However, this paper does not optimize the configuration from the aspects of LCC and LPSP.At present, most of the research on the optimal sizing of the stand-alone photovoltaic system focuses on the power supply-demand for residential households or rural areas. Abbes et al. [10] proposed a new methodology with the aim to design an autonomous hybrid PV-wind-battery system. The electric demand of a residential house is taken as an example to verify the proposed method. Bataineh et al. [11] designed a stand-alone photovoltaic (PV) system to provide the required electricity for a single residential household in a rural area in Jordan. Koutroulis et al. [12] presented a methodology for optimal sizing of stand-alone photovoltaic/wind generator systems. The proposed method has been applied to design a power generation system that supplies a residential household. Okoye et al. [13] studied the optimal sizing of stand-alone photovoltaic systems in residential buildings. As mentioned above, there are a few studies on the walking power demand and photovoltaic power matching design for the solar-powered sprinkler irrigation machine. Therefore, this work aims to optimize the stand-alone photovoltaic system for a solar-powered translational sprinkler irrigation machine considering the loss of power supply probability.The rest of this paper is organized as follows. In Section2, the stand-alone photovoltaic system configuration model is introduced. In Section 3, the stand-alone photovoltaic system is optimized considering the loss of power supply probability. In Section 4, the optimal combination of PV module and battery are obtained, and relevant discussions are carried out. Conclusions are drawn in Section 5.


## 2. Stand-Alone Photovoltaic System Configuration Model

### 2.1. The Self-Developed Translational Sprinkler Irrigation Machine

This study uses the self-developed translational sprinkler irrigation machine as the research objective, and its structural diagram is shown in Figure1.Figure 1
Structural diagrams of the developed translational sprinkler irrigation machine.Figure1 mainly introduces the main structure of the self-developed translational sprinkler irrigation machine. The self-developed translational sprinkler irrigation machine comprises the solar photovoltaic power supply, mechanical structure, and control system parts. As shown in Figure 1, the self-developed sprinkler irrigation machine mainly includes the solar photovoltaic module, battery, main water pipe, water pump, steering motor, high-clearance chassis, control cabinet, walking motor reducer, nozzle, and so on. The detailed information of the developed translational sprinkler irrigation machine is described in [14].


### 2.2. Driving Power of the Translational Sprinkler Irrigation Machine

The driving power of the translational sprinkler irrigation machinePd can be given by(1)Pd=mgfvmaxηα,where m is the maximum mass of the translational sprinkler irrigation machine after filling water into the water tank, g is the acceleration of gravity, f is the rolling resistance coefficient, vmax is the maximum travelling speed of the translational sprinkler irrigation machine, a is the safety factor, and η is the mechanical transmission efficiency. For this study, the parameters m, g, f, vmax, α, and η are chosen as 3500 kg, 9.8 m/s2, 0.3, 0.028 m/s, 1.2, and 0.8, respectively.


### 2.3. Water Intake Pressure Driving Power

The translational sprinkler irrigation machine is mainly used for crop irrigation, pesticide spraying, and fertilization. According to the inflow pressure and flow requirements of the unit design, the water intake pressure driving power required by sprinkler irrigation unit can be expressed as follows:(2)Pp=ρgQH3600ηPα,where ρ is the water density, Q is the flow of the pumping unit, H is the water-lift of the self-priming pump of the translational sprinkler irrigation machine, α is the safety factor, and ηp is the efficiency of the pump pressurization system, which includes water pump efficiency and DC brushless motor efficiency. For this study, the parameters ρ, Q, H, α, and ηp are chosen as 1000 kg/m3, 12 m3/h, 10 m, 1.2, and 0.54, respectively.


### 2.4. The Loss Power of the Controller and Sensors

The loss power of the controller and sensors of the unitPc is also taken into consideration. In this study, Pc is chosen as 100 W.The total load powerPl can be represented as follows:(3)Pl=Pd+Pp+Pc.


### 2.5. Photovoltaic Generator Model

The output electric power of the photovoltaic generator can be written as follows:(4)Ppv=Npv×ηpv×Apv×Ir,where Npv is the number of PV modules, ηpv is the power conversion efficiency, Apv is the surface area of PV panels of a single module, and Ir is the solar radiation. The power conversion efficiency ηpv can be expressed by the following equation:(5)ηpv=ηr×ηpc×1−βgetc×Tc−NOCT,where ηr is the reference module efficiency, ηpc is the power condition efficiency, βgetc is the generator efficiency temperature coefficient, Tc is the cell temperature, and NOCT is the normal operating cell temperature when cells operate under standard operating conditions. Tc can be estimated as follows:(6)TC=30+0.0175×Ir−300+1.14×Ta−25,where Ta is ambient temperature (°C).


### 2.6. Battery Storage Model

In this study, the lead-acid battery is used in the stand-alone photovoltaic system due to low cost, high reliability, convenient maintenance, and so on. The state of charge (SOC) is determined by using the difference between the generated PV power and the load power. When the generated PV power is greater than the load power, the energy produced by the PV system can satisfy the demand of the load. Meanwhile, the extra energy produced by the PV system will be stored in the battery, which corresponds to the charging process. On the contrary, when the generated PV power is less than the load power, the demand of load cannot be satisfied. The stored energy will be released to guarantee the demand of load, which corresponds to the discharging process.The state of charge (SOC) is used to reflect the battery’s remaining power in real-time. It can be defined by using the ratio of the battery's remaining capacity to the rated capacity at a certain time.When charging, the SOC of the battery at the timet is expressed as follows:(7)SOCt=SOCt−Δt+ΔEstoreηinNbErate.When discharging, the SOC of the battery at the timet is expressed as follows:(8)SOCt=SOCt−Δt−ΔEstoreηoutNbErate,where SOCt and SOCt−Δt are the states of charge of the battery at the time t and t−Δt, respectively, ΔEstore is the amount of energy charged or discharged by the battery during period Δt, ηin and ηout are the charging and discharging efficiency of the battery system, respectively, Nb is the number of batteries, and ΔEstore is the rated capacity of the battery.When the PV power is greater than the load power, the extra PV power will be stored in the battery. The theoretical charging capacity of the battery within timet and t+Δt is expressed as follows:(9)ΔEstore=Ppvt−PltηinvΔt.When the PV power is less than the load power, and the sum of the PV power and the battery discharge power is greater than the load power, the battery needs to release power to meet the load demand. The theoretical discharging capacity of the battery within timet and t+Δt is expressed as follows:(10)ΔEstore=Pltηinv−PpvtΔt.The state of charge of the battery should not exceed the upper limitSOCmax and not be lower than the lower limit SOCmin, which is expressed as follows:(11)SOCmin≤SOCt≤SOCmax.WhenSOCt+Δt<SOCmin, the actual discharge capacity of the battery during time t and t+Δt can be expressed as follows:(12)ΔEdischarge=NbErateSOCt−SOCmin.When the sum of the PV power and the battery discharge power is less than the load power, the loss of power supply (LPS) will happen. The capacity of the LPS, i.e.,QLPSt, can be calculated as follows:(13)QLPSt=Plt−Ppvt+Pdischarge×ηinvΔt,where Pdischarge is the discharge power of the battery.WhenSOCt+Δt>SOCmax, the actual charge capacity of the battery during time t and t+Δt can be expressed as follows:(14)ΔEcharge=NbErateSOCmax−SOCt.In this study, it is assumed that the battery is full charged. That is, the initial value of SOC is chosen as 100%.


## 2.1. The Self-Developed Translational Sprinkler Irrigation Machine

This study uses the self-developed translational sprinkler irrigation machine as the research objective, and its structural diagram is shown in Figure1.Figure 1
Structural diagrams of the developed translational sprinkler irrigation machine.Figure1 mainly introduces the main structure of the self-developed translational sprinkler irrigation machine. The self-developed translational sprinkler irrigation machine comprises the solar photovoltaic power supply, mechanical structure, and control system parts. As shown in Figure 1, the self-developed sprinkler irrigation machine mainly includes the solar photovoltaic module, battery, main water pipe, water pump, steering motor, high-clearance chassis, control cabinet, walking motor reducer, nozzle, and so on. The detailed information of the developed translational sprinkler irrigation machine is described in [14].


## 2.2. Driving Power of the Translational Sprinkler Irrigation Machine

The driving power of the translational sprinkler irrigation machinePd can be given by(1)Pd=mgfvmaxηα,where m is the maximum mass of the translational sprinkler irrigation machine after filling water into the water tank, g is the acceleration of gravity, f is the rolling resistance coefficient, vmax is the maximum travelling speed of the translational sprinkler irrigation machine, a is the safety factor, and η is the mechanical transmission efficiency. For this study, the parameters m, g, f, vmax, α, and η are chosen as 3500 kg, 9.8 m/s2, 0.3, 0.028 m/s, 1.2, and 0.8, respectively.


## 2.3. Water Intake Pressure Driving Power

The translational sprinkler irrigation machine is mainly used for crop irrigation, pesticide spraying, and fertilization. According to the inflow pressure and flow requirements of the unit design, the water intake pressure driving power required by sprinkler irrigation unit can be expressed as follows:(2)Pp=ρgQH3600ηPα,where ρ is the water density, Q is the flow of the pumping unit, H is the water-lift of the self-priming pump of the translational sprinkler irrigation machine, α is the safety factor, and ηp is the efficiency of the pump pressurization system, which includes water pump efficiency and DC brushless motor efficiency. For this study, the parameters ρ, Q, H, α, and ηp are chosen as 1000 kg/m3, 12 m3/h, 10 m, 1.2, and 0.54, respectively.


## 2.4. The Loss Power of the Controller and Sensors

The loss power of the controller and sensors of the unitPc is also taken into consideration. In this study, Pc is chosen as 100 W.The total load powerPl can be represented as follows:(3)Pl=Pd+Pp+Pc.


## 2.5. Photovoltaic Generator Model

The output electric power of the photovoltaic generator can be written as follows:(4)Ppv=Npv×ηpv×Apv×Ir,where Npv is the number of PV modules, ηpv is the power conversion efficiency, Apv is the surface area of PV panels of a single module, and Ir is the solar radiation. The power conversion efficiency ηpv can be expressed by the following equation:(5)ηpv=ηr×ηpc×1−βgetc×Tc−NOCT,where ηr is the reference module efficiency, ηpc is the power condition efficiency, βgetc is the generator efficiency temperature coefficient, Tc is the cell temperature, and NOCT is the normal operating cell temperature when cells operate under standard operating conditions. Tc can be estimated as follows:(6)TC=30+0.0175×Ir−300+1.14×Ta−25,where Ta is ambient temperature (°C).


## 2.6. Battery Storage Model

In this study, the lead-acid battery is used in the stand-alone photovoltaic system due to low cost, high reliability, convenient maintenance, and so on. The state of charge (SOC) is determined by using the difference between the generated PV power and the load power. When the generated PV power is greater than the load power, the energy produced by the PV system can satisfy the demand of the load. Meanwhile, the extra energy produced by the PV system will be stored in the battery, which corresponds to the charging process. On the contrary, when the generated PV power is less than the load power, the demand of load cannot be satisfied. The stored energy will be released to guarantee the demand of load, which corresponds to the discharging process.The state of charge (SOC) is used to reflect the battery’s remaining power in real-time. It can be defined by using the ratio of the battery's remaining capacity to the rated capacity at a certain time.When charging, the SOC of the battery at the timet is expressed as follows:(7)SOCt=SOCt−Δt+ΔEstoreηinNbErate.When discharging, the SOC of the battery at the timet is expressed as follows:(8)SOCt=SOCt−Δt−ΔEstoreηoutNbErate,where SOCt and SOCt−Δt are the states of charge of the battery at the time t and t−Δt, respectively, ΔEstore is the amount of energy charged or discharged by the battery during period Δt, ηin and ηout are the charging and discharging efficiency of the battery system, respectively, Nb is the number of batteries, and ΔEstore is the rated capacity of the battery.When the PV power is greater than the load power, the extra PV power will be stored in the battery. The theoretical charging capacity of the battery within timet and t+Δt is expressed as follows:(9)ΔEstore=Ppvt−PltηinvΔt.When the PV power is less than the load power, and the sum of the PV power and the battery discharge power is greater than the load power, the battery needs to release power to meet the load demand. The theoretical discharging capacity of the battery within timet and t+Δt is expressed as follows:(10)ΔEstore=Pltηinv−PpvtΔt.The state of charge of the battery should not exceed the upper limitSOCmax and not be lower than the lower limit SOCmin, which is expressed as follows:(11)SOCmin≤SOCt≤SOCmax.WhenSOCt+Δt<SOCmin, the actual discharge capacity of the battery during time t and t+Δt can be expressed as follows:(12)ΔEdischarge=NbErateSOCt−SOCmin.When the sum of the PV power and the battery discharge power is less than the load power, the loss of power supply (LPS) will happen. The capacity of the LPS, i.e.,QLPSt, can be calculated as follows:(13)QLPSt=Plt−Ppvt+Pdischarge×ηinvΔt,where Pdischarge is the discharge power of the battery.WhenSOCt+Δt>SOCmax, the actual charge capacity of the battery during time t and t+Δt can be expressed as follows:(14)ΔEcharge=NbErateSOCmax−SOCt.In this study, it is assumed that the battery is full charged. That is, the initial value of SOC is chosen as 100%.


## 3. Optimal Sizing of the Stand-Alone Photovoltaic System

### 3.1. Objective Function

Optimal sizing of the stand-alone photovoltaic system aims to minimize the life cycle cost (LCC) of the system on the premise of meeting the load power demand and power supply reliability. In this study, the LCC includes equipment costCE, the present worth of equipment replacement cost after n years CEnPW, installation cost CInst, and the present worth of maintenance cost CMPW. Therefore, the objective function of the LCC for the stand-alone photovoltaic system can be expressed as(15)minLCC=CE+CEnPW+CInst+CMPW,where the equipment cost CE mainly contains the cost of the PV modules, batteries, and controller, which are represented by CPV, CB, and Cc, respectively. Accordingly, the equipment cost CE can be represented as(16)CE=CPV+CB+Cc.In our study, the cost of each photovoltaic panel is 1000 RMB, the cost of each battery is 600 RMB, and the cost of the controller is 1600 RMB. Then, equation (16) can be rewritten as(17)CE=1000×Npv+600×Nb+1600.In this work, the operating life of the photovoltaic system N is set as 20 years, namely,N = 20. The life of the battery is assumed to be five years, and the life of the PV modules and controller is assumed to be 20 years. Therefore, only the present worth of the battery replacement cost after n years i.e., CBnPW, should be considered in the operating life of the photovoltaic system, i.e., CEnPW  =  CBnPW. The present worth of the battery replacement cost after n years can be calculated as follows:(18)CBnPW=CB1+i1+dn,where i and d represent the inflation rate and discount rate, respectively. During calculation, the inflation rate i and discount rate d are taken as 3% and 5%, respectively. Since the operating life of the photovoltaic system is 20 years and the life of the battery is assumed to be five years, n is chosen as 5, 10, and 15. Substituting the value of the parameters into equation (18), the following equation can be obtained:(19)CBnPW=600×Nb×1+3%1+5%5+1+3%1+5%10+1+3%1+5%15=1489.7×Nb,The installation cost is assumed 10% of the PV modules cost, which can be expressed as follows:(20)CInst=10%×1000×Npv=100×Npv.The present worth of maintenance cost CMPW can be acquired as follows:(21)CMPW=CMPY×1+i1+d×1−1+i/1+dN1−1+i/1+d,where CMPY is the maintenance cost per year, which is assumed 2% of the PV modules cost. Substituting the parameter value into equation (21), we can get(22)CMPY=2%×1000×Npv×1+3%1+5%×1−1+3%/1+5%201−1+3%/1+5%=328.87×Npv.Combining equations (17), (19), (20), and (22), the life cycle cost of the system LCC can be expressed as follows:(23)LCC=1428.87×Npv+2089.7×Nb+1600.


### 3.2. Optimization Constraints

#### 3.2.1. Power Supply Reliability Constraint

The loss of power supply probability (LPSP) represents the probability that the photovoltaic system cannot meet the load power demand within the evaluation periodT. The LPSP can be expressed as the ratio of load power shortage to total load demand in the time period. In this study, LPSP is taken as the power supply reliability index, which can be calculated as follows:(24)LPSPt=∑t=t0t0+kΔtPlt−PPVt+Pdischarget×ηinvΔt∑t=t0t0+kΔtPltΔt,where k is the time series for evaluation. When LPSP = 0, it indicates that the power supply of the system can be satisfied; when LPSP = 1, it indicates that the power supply of the system cannot be satisfied. Therefore, the constraint corresponding to the power supply reliability is expressed as(25)g1=maxLPSPt−LPSPmax≤0.


#### 3.2.2. Battery Power Constraint

To prevent the service life of the battery from being affected by over charge or over discharge, the state of charge of the battery should meet the following constraint:(26)SOCmin≤SOCt≤SOCmax.In this work, the SOCmin and SOCmax are chosen as the 20% and 80%, respectively.To better demonstrate the optimization process and algorithm logic, the corresponding flow chart is presented, as shown in Figure2.Figure 2
Flow chart of the optimization process and algorithm logic.


## 3.1. Objective Function

Optimal sizing of the stand-alone photovoltaic system aims to minimize the life cycle cost (LCC) of the system on the premise of meeting the load power demand and power supply reliability. In this study, the LCC includes equipment costCE, the present worth of equipment replacement cost after n years CEnPW, installation cost CInst, and the present worth of maintenance cost CMPW. Therefore, the objective function of the LCC for the stand-alone photovoltaic system can be expressed as(15)minLCC=CE+CEnPW+CInst+CMPW,where the equipment cost CE mainly contains the cost of the PV modules, batteries, and controller, which are represented by CPV, CB, and Cc, respectively. Accordingly, the equipment cost CE can be represented as(16)CE=CPV+CB+Cc.In our study, the cost of each photovoltaic panel is 1000 RMB, the cost of each battery is 600 RMB, and the cost of the controller is 1600 RMB. Then, equation (16) can be rewritten as(17)CE=1000×Npv+600×Nb+1600.In this work, the operating life of the photovoltaic system N is set as 20 years, namely,N = 20. The life of the battery is assumed to be five years, and the life of the PV modules and controller is assumed to be 20 years. Therefore, only the present worth of the battery replacement cost after n years i.e., CBnPW, should be considered in the operating life of the photovoltaic system, i.e., CEnPW  =  CBnPW. The present worth of the battery replacement cost after n years can be calculated as follows:(18)CBnPW=CB1+i1+dn,where i and d represent the inflation rate and discount rate, respectively. During calculation, the inflation rate i and discount rate d are taken as 3% and 5%, respectively. Since the operating life of the photovoltaic system is 20 years and the life of the battery is assumed to be five years, n is chosen as 5, 10, and 15. Substituting the value of the parameters into equation (18), the following equation can be obtained:(19)CBnPW=600×Nb×1+3%1+5%5+1+3%1+5%10+1+3%1+5%15=1489.7×Nb,The installation cost is assumed 10% of the PV modules cost, which can be expressed as follows:(20)CInst=10%×1000×Npv=100×Npv.The present worth of maintenance cost CMPW can be acquired as follows:(21)CMPW=CMPY×1+i1+d×1−1+i/1+dN1−1+i/1+d,where CMPY is the maintenance cost per year, which is assumed 2% of the PV modules cost. Substituting the parameter value into equation (21), we can get(22)CMPY=2%×1000×Npv×1+3%1+5%×1−1+3%/1+5%201−1+3%/1+5%=328.87×Npv.Combining equations (17), (19), (20), and (22), the life cycle cost of the system LCC can be expressed as follows:(23)LCC=1428.87×Npv+2089.7×Nb+1600.


## 3.2. Optimization Constraints

### 3.2.1. Power Supply Reliability Constraint

The loss of power supply probability (LPSP) represents the probability that the photovoltaic system cannot meet the load power demand within the evaluation periodT. The LPSP can be expressed as the ratio of load power shortage to total load demand in the time period. In this study, LPSP is taken as the power supply reliability index, which can be calculated as follows:(24)LPSPt=∑t=t0t0+kΔtPlt−PPVt+Pdischarget×ηinvΔt∑t=t0t0+kΔtPltΔt,where k is the time series for evaluation. When LPSP = 0, it indicates that the power supply of the system can be satisfied; when LPSP = 1, it indicates that the power supply of the system cannot be satisfied. Therefore, the constraint corresponding to the power supply reliability is expressed as(25)g1=maxLPSPt−LPSPmax≤0.


### 3.2.2. Battery Power Constraint

To prevent the service life of the battery from being affected by over charge or over discharge, the state of charge of the battery should meet the following constraint:(26)SOCmin≤SOCt≤SOCmax.In this work, the SOCmin and SOCmax are chosen as the 20% and 80%, respectively.To better demonstrate the optimization process and algorithm logic, the corresponding flow chart is presented, as shown in Figure2.Figure 2
Flow chart of the optimization process and algorithm logic.


## 3.2.1. Power Supply Reliability Constraint

The loss of power supply probability (LPSP) represents the probability that the photovoltaic system cannot meet the load power demand within the evaluation periodT. The LPSP can be expressed as the ratio of load power shortage to total load demand in the time period. In this study, LPSP is taken as the power supply reliability index, which can be calculated as follows:(24)LPSPt=∑t=t0t0+kΔtPlt−PPVt+Pdischarget×ηinvΔt∑t=t0t0+kΔtPltΔt,where k is the time series for evaluation. When LPSP = 0, it indicates that the power supply of the system can be satisfied; when LPSP = 1, it indicates that the power supply of the system cannot be satisfied. Therefore, the constraint corresponding to the power supply reliability is expressed as(25)g1=maxLPSPt−LPSPmax≤0.


## 3.2.2. Battery Power Constraint

To prevent the service life of the battery from being affected by over charge or over discharge, the state of charge of the battery should meet the following constraint:(26)SOCmin≤SOCt≤SOCmax.In this work, the SOCmin and SOCmax are chosen as the 20% and 80%, respectively.To better demonstrate the optimization process and algorithm logic, the corresponding flow chart is presented, as shown in Figure2.Figure 2
Flow chart of the optimization process and algorithm logic.


## 4. Results and Discussion

### 4.1. Determination of the Optimal Combination of PV Module and Battery

Solar irradiance and ambient temperature are important input parameters of the optimization model. The AV6592 type portable tester (voltage accuracy 0.01 V, current accuracy 0.001 A, and power range 0.1∼500 W) was used to monitor the solar irradiance and ambient temperature in the test area (east longitude 108.07°, northern latitude 34.28°, and height above sea level 521 m). The single-crystal photovoltaic panel type CS5M32-260 and the valve regulated fully sealed lead-acid battery typed 190H52 (capacity 200 A·h, rated voltage 12 V) are used in this work. The monitoring time is from April 29 to May 29, 2018, a total of 31 days. During the days of monitoring, the translational sprinkler irrigation machine works 8 hours a day, from 9 a.m. to 5 p.m. The solar irradiance and ambient temperature within 248 hours are shown in Figure3.Figure 3
Solar irradiance and ambient temperature during the total monitoring time. (a) Solar irradiance. (b) Ambient temperature.
(a)(b)As shown in Figure3, the solar irradiance and ambient temperature within the total monitoring time show randomness and irregular variation due to the complex climatic conditions. In this study, the developed translational sprinkler irrigation machine has a lightweight and small- and medium-size structure powered by the stand-alone photovoltaic system. Suppose the sprinkler irrigation machine is guaranteed to operate continuously under non-sunny weather conditions. In that case, more photovoltaic modules and batteries will be required, which may exceed the bearing range of the design structure. Therefore, the developed translational sprinkler irrigation machine has the best applicability under sunny weather conditions. In addition, sprinkler irrigation is more necessary in sunny weather conditions. The solar irradiance and ambient temperature of a typical sunny day in the test area are shown in Figure 4.Figure 4
Solar irradiance and ambient temperature of a typical sunny day in the test area. (a) Solar irradiance. (b) Ambient temperature.
(a)(b)In Figure4(b), the red curve is the fitting of the ambient temperature for a typical sunny day. It can be seen from Figure 4 that under the typical sunny weather condition in the test area the solar irradiance and ambient temperature reach the maximum between 13:00 and 14:00. The other main input parameters are listed in Table 1.Table 1
The other main input parameters.
ParametersValueηr0.12ηpc0.9βgetc0.0045NOCT45Apv (m2)1.5ηin0.9ηout0.85Erate (A·h)200ηinv0.95As mentioned above, when the photovoltaic (PV) system cannot meet the load power demand, the loss of power supply will happen. The loss of power supply probability (LPSP) varies with the combinations of PV module and battery. Considering the small- and medium-size structure of the developed translational sprinkler irrigation machine, the number of photovoltaic modules and batteries should not be too large. In this study, the number of photovoltaic modules ranges from 1 to 25, and the number of batteries also ranges from 1 to 25. That is, 1 ≤Npv ≤ 25, and 1 ≤ Nb ≤ 25. The relationship among the number of photovoltaic modules, the number of batteries, and the LPSPs is shown in Figure 5.Figure 5
Relationship among the number of photovoltaic modules, the number of batteries, and the LPSPs.In Figure5, the curves are the boundaries between the desirable and undesirable LPSPs. It can be seen from Figure 5 that the larger the LPSP, the more combinations of PV module and battery can be selected. The purpose of optimal sizing is to minimize the life cycle cost (LCC) on the premise that LPSP meets the requirements. The life cycle cost (LCC) comparison for different combinations of photovoltaic module and battery is shown in Table 2.Table 2
LCC comparison for different combinations of photovoltaic module and battery.
No.Life cycle cost (LCC)Loss of power supply probability (LPSP)Number of PV modules (Npv) (%)Number of batteries (Nb)114245.4176.4734217763.9870.6845321282.5531.8456424801.1224.4367528319.6922.4678631838.2612.6589735356.834.74910838875.4001011942393.97011121045912.54012131149431.11013141252949.68014151356468.25015161459986.82016171528641.320162As shown in Table2, with the increase of the number of PV modules and batteries, the LPSPs decrease, but the corresponding LCC increases. In addition, when Npv = 16 and Nb = 2, it can also make LPSP = 0. The relationship between LPSP and LCC is shown in Figure 6.Figure 6
Relationship between LPSP and LCC.As can be seen from Figure6, we mainly divide the parameter plane into three regions, i.e., regions A, B, and C. The LPSP corresponding to region A is small, but the smaller LPSP is based on the sacrifice of LCC. As for region B, the corresponding LPSP is greater than that in region A, while the corresponding LCC is not reduced much. The LPSP corresponding to region C is larger than that in regions A and B, but the corresponding LCC is smaller than that in regions A and B. In this study, to ensure the continuous operation of the sprinkler irrigation machine, the maximum LPSP is set to 0; that is, there is no power loss during the operation of the sprinkler irrigation machine under typical sunny weather conditions.After comparative analysis, the combination of PV module and battery withNpv = 16 and Nb = 2 is proved the optimal one, the LCC of which is minimum and whose LPSP also meets the requirements (LPSP = 0).


### 4.2. Comparison with the Existing Method

In this section, by combining the data listed in Table2, the optimization sizing method presented in this manuscript is compared with the existing method reported in [7]. In [7], the technical and economic feasibility of photovoltaic pumping of water is studied. However, the loss of power supply probability (LPSP) is not considered in [7].According to Table2, the optimal combination of PV module and battery determined by the proposed method corresponds to No. 15 with Npv = 16 and Nb = 2, and the minimum LCC is 28641.32 RMB. Besides, the sprinkler irrigation machine can operate continuously without loss of power supply. If the loss of power supply probability (LPSP) is not considered, the optimal combination of PV module and battery obtained using the method presented in [7] may correspond to No. 1 with Npv = 3 and Nb = 4, and the minimum LCC is 14245.41 RMB. Although the LCC obtained using the method presented in [7] is relatively smaller, the corresponding LPSP is as high as 76.47%. Consequently, the machine cannot work continuously, and the effective working time is very short.


### 4.3. Operation Test of Sprinkler Irrigation Machine

To verify the optimal sizing results, the solar irradiance and ambient temperature on June 5, 2018 (typical sunny day) are observed, as shown in Figure7. Meanwhile, the sprinkler irrigation operation is also carried out using the developed sprinkler irrigation machine in the test area.Figure 7
Solar irradiance and ambient temperature of the test day. (a) Solar irradiance. (b) Ambient temperature.
(a)(b)It can be seen from Figure7 that the overall change trend of solar irradiance and ambient temperature on the test day is consistent with that on the typical sunny day described in Figure 4. The solar irradiance increases first and then decreases and reaches the maximum between 13:00 and 14:00. The ambient temperature rises first and then decreases slightly.After substituting the optimal combination of PV module and battery (Npv = 16 and Nb = 2) into the PV system model and battery storage model, the PV power, state of charge (SOC), and load power at each time can be obtained, as shown in Figure 8. Note that the factors affecting the operation of sprinkler irrigation machines are complicated. The load power at each time is assumed to be equal to the full load power in the calculation process to ensure the normal operation of the sprinkler irrigation machine.Figure 8
The PV power, state of charge (SOC), and load power at different time points.As shown in Figure8, the PV power is greater than the load power at most time points. In Figure 8, between 9:00 and 10:00, the PV power is less than the load power due to the low irradiance and ambient temperature. Consequently, the battery needs to be discharged to meet the load demand. Accordingly, the state of charge (SOC) of the battery also shows a downward trend. From 10:00 to 16:30, the PV power is greater than the load power, which can meet the demand of the load. During this time period, the SOC curve fluctuates with time. The SOC value rises because the PV power is greater than the load power, and the SOC is calculated according to equation (7). The decrease of SOC value is because SOC is greater than SOCmax, and SOC is calculated according to equation (14). From 16:30 to 17:00, the solar irradiance and ambient temperature decrease, which makes the PV power lower than the load power, and then the SOC also shows a downward trend.In order to further verify the practicability of the self-developed translational sprinkler irrigation machine and the feasibility of the optimal sizing result, the sprinkler irrigation experiment using the developed sprinkler irrigation machine was carried out in the experimental site, as shown in Figure9.Figure 9
Sprinkler irrigation experiment site using the developed sprinkler irrigation machine.The developed sprinkler irrigation machine, which installed 16 photovoltaic modules and two batteries, runs from 9:00 to 17:00 for a total of 8 hours, and the maximum running speed is 0.028 m/s. During operation, the flow of the pumping unit is 12 m3/h, and the water-lift of the self-priming pump of the translational sprinkler irrigation machine is 10 m. Within 8 hours of operation, the sprinkler irrigation machine operated continuously without loss of power supply, which further verified the effectiveness of the optimal sizing results.


## 4.1. Determination of the Optimal Combination of PV Module and Battery

Solar irradiance and ambient temperature are important input parameters of the optimization model. The AV6592 type portable tester (voltage accuracy 0.01 V, current accuracy 0.001 A, and power range 0.1∼500 W) was used to monitor the solar irradiance and ambient temperature in the test area (east longitude 108.07°, northern latitude 34.28°, and height above sea level 521 m). The single-crystal photovoltaic panel type CS5M32-260 and the valve regulated fully sealed lead-acid battery typed 190H52 (capacity 200 A·h, rated voltage 12 V) are used in this work. The monitoring time is from April 29 to May 29, 2018, a total of 31 days. During the days of monitoring, the translational sprinkler irrigation machine works 8 hours a day, from 9 a.m. to 5 p.m. The solar irradiance and ambient temperature within 248 hours are shown in Figure3.Figure 3
Solar irradiance and ambient temperature during the total monitoring time. (a) Solar irradiance. (b) Ambient temperature.
(a)(b)As shown in Figure3, the solar irradiance and ambient temperature within the total monitoring time show randomness and irregular variation due to the complex climatic conditions. In this study, the developed translational sprinkler irrigation machine has a lightweight and small- and medium-size structure powered by the stand-alone photovoltaic system. Suppose the sprinkler irrigation machine is guaranteed to operate continuously under non-sunny weather conditions. In that case, more photovoltaic modules and batteries will be required, which may exceed the bearing range of the design structure. Therefore, the developed translational sprinkler irrigation machine has the best applicability under sunny weather conditions. In addition, sprinkler irrigation is more necessary in sunny weather conditions. The solar irradiance and ambient temperature of a typical sunny day in the test area are shown in Figure 4.Figure 4
Solar irradiance and ambient temperature of a typical sunny day in the test area. (a) Solar irradiance. (b) Ambient temperature.
(a)(b)In Figure4(b), the red curve is the fitting of the ambient temperature for a typical sunny day. It can be seen from Figure 4 that under the typical sunny weather condition in the test area the solar irradiance and ambient temperature reach the maximum between 13:00 and 14:00. The other main input parameters are listed in Table 1.Table 1
The other main input parameters.
ParametersValueηr0.12ηpc0.9βgetc0.0045NOCT45Apv (m2)1.5ηin0.9ηout0.85Erate (A·h)200ηinv0.95As mentioned above, when the photovoltaic (PV) system cannot meet the load power demand, the loss of power supply will happen. The loss of power supply probability (LPSP) varies with the combinations of PV module and battery. Considering the small- and medium-size structure of the developed translational sprinkler irrigation machine, the number of photovoltaic modules and batteries should not be too large. In this study, the number of photovoltaic modules ranges from 1 to 25, and the number of batteries also ranges from 1 to 25. That is, 1 ≤Npv ≤ 25, and 1 ≤ Nb ≤ 25. The relationship among the number of photovoltaic modules, the number of batteries, and the LPSPs is shown in Figure 5.Figure 5
Relationship among the number of photovoltaic modules, the number of batteries, and the LPSPs.In Figure5, the curves are the boundaries between the desirable and undesirable LPSPs. It can be seen from Figure 5 that the larger the LPSP, the more combinations of PV module and battery can be selected. The purpose of optimal sizing is to minimize the life cycle cost (LCC) on the premise that LPSP meets the requirements. The life cycle cost (LCC) comparison for different combinations of photovoltaic module and battery is shown in Table 2.Table 2
LCC comparison for different combinations of photovoltaic module and battery.
No.Life cycle cost (LCC)Loss of power supply probability (LPSP)Number of PV modules (Npv) (%)Number of batteries (Nb)114245.4176.4734217763.9870.6845321282.5531.8456424801.1224.4367528319.6922.4678631838.2612.6589735356.834.74910838875.4001011942393.97011121045912.54012131149431.11013141252949.68014151356468.25015161459986.82016171528641.320162As shown in Table2, with the increase of the number of PV modules and batteries, the LPSPs decrease, but the corresponding LCC increases. In addition, when Npv = 16 and Nb = 2, it can also make LPSP = 0. The relationship between LPSP and LCC is shown in Figure 6.Figure 6
Relationship between LPSP and LCC.As can be seen from Figure6, we mainly divide the parameter plane into three regions, i.e., regions A, B, and C. The LPSP corresponding to region A is small, but the smaller LPSP is based on the sacrifice of LCC. As for region B, the corresponding LPSP is greater than that in region A, while the corresponding LCC is not reduced much. The LPSP corresponding to region C is larger than that in regions A and B, but the corresponding LCC is smaller than that in regions A and B. In this study, to ensure the continuous operation of the sprinkler irrigation machine, the maximum LPSP is set to 0; that is, there is no power loss during the operation of the sprinkler irrigation machine under typical sunny weather conditions.After comparative analysis, the combination of PV module and battery withNpv = 16 and Nb = 2 is proved the optimal one, the LCC of which is minimum and whose LPSP also meets the requirements (LPSP = 0).


## 4.2. Comparison with the Existing Method

In this section, by combining the data listed in Table2, the optimization sizing method presented in this manuscript is compared with the existing method reported in [7]. In [7], the technical and economic feasibility of photovoltaic pumping of water is studied. However, the loss of power supply probability (LPSP) is not considered in [7].According to Table2, the optimal combination of PV module and battery determined by the proposed method corresponds to No. 15 with Npv = 16 and Nb = 2, and the minimum LCC is 28641.32 RMB. Besides, the sprinkler irrigation machine can operate continuously without loss of power supply. If the loss of power supply probability (LPSP) is not considered, the optimal combination of PV module and battery obtained using the method presented in [7] may correspond to No. 1 with Npv = 3 and Nb = 4, and the minimum LCC is 14245.41 RMB. Although the LCC obtained using the method presented in [7] is relatively smaller, the corresponding LPSP is as high as 76.47%. Consequently, the machine cannot work continuously, and the effective working time is very short.


## 4.3. Operation Test of Sprinkler Irrigation Machine

To verify the optimal sizing results, the solar irradiance and ambient temperature on June 5, 2018 (typical sunny day) are observed, as shown in Figure7. Meanwhile, the sprinkler irrigation operation is also carried out using the developed sprinkler irrigation machine in the test area.Figure 7
Solar irradiance and ambient temperature of the test day. (a) Solar irradiance. (b) Ambient temperature.
(a)(b)It can be seen from Figure7 that the overall change trend of solar irradiance and ambient temperature on the test day is consistent with that on the typical sunny day described in Figure 4. The solar irradiance increases first and then decreases and reaches the maximum between 13:00 and 14:00. The ambient temperature rises first and then decreases slightly.After substituting the optimal combination of PV module and battery (Npv = 16 and Nb = 2) into the PV system model and battery storage model, the PV power, state of charge (SOC), and load power at each time can be obtained, as shown in Figure 8. Note that the factors affecting the operation of sprinkler irrigation machines are complicated. The load power at each time is assumed to be equal to the full load power in the calculation process to ensure the normal operation of the sprinkler irrigation machine.Figure 8
The PV power, state of charge (SOC), and load power at different time points.As shown in Figure8, the PV power is greater than the load power at most time points. In Figure 8, between 9:00 and 10:00, the PV power is less than the load power due to the low irradiance and ambient temperature. Consequently, the battery needs to be discharged to meet the load demand. Accordingly, the state of charge (SOC) of the battery also shows a downward trend. From 10:00 to 16:30, the PV power is greater than the load power, which can meet the demand of the load. During this time period, the SOC curve fluctuates with time. The SOC value rises because the PV power is greater than the load power, and the SOC is calculated according to equation (7). The decrease of SOC value is because SOC is greater than SOCmax, and SOC is calculated according to equation (14). From 16:30 to 17:00, the solar irradiance and ambient temperature decrease, which makes the PV power lower than the load power, and then the SOC also shows a downward trend.In order to further verify the practicability of the self-developed translational sprinkler irrigation machine and the feasibility of the optimal sizing result, the sprinkler irrigation experiment using the developed sprinkler irrigation machine was carried out in the experimental site, as shown in Figure9.Figure 9
Sprinkler irrigation experiment site using the developed sprinkler irrigation machine.The developed sprinkler irrigation machine, which installed 16 photovoltaic modules and two batteries, runs from 9:00 to 17:00 for a total of 8 hours, and the maximum running speed is 0.028 m/s. During operation, the flow of the pumping unit is 12 m3/h, and the water-lift of the self-priming pump of the translational sprinkler irrigation machine is 10 m. Within 8 hours of operation, the sprinkler irrigation machine operated continuously without loss of power supply, which further verified the effectiveness of the optimal sizing results.


## 5. Conclusions

This paper aims to optimize the stand-alone photovoltaic system for a solar-powered translational sprinkler irrigation machine. The following conclusions can be drawn:(1)
The load power, which includes the driving power of the translational sprinkler irrigation machine, water intake pressure driving power, and the loss power of the controller and sensors, is considered in the calculation process. Meanwhile, the photovoltaic generator model and the battery storage model are established.(2)
The stand-alone photovoltaic system is optimized considering the loss of power supply probability (LPSP) and the life cycle cost (LCC). The combination of PV module and battery withNpv = 16 and Nb = 2 is proved the optimal one.(3)
By analyzing the PV power, state of charge (SOC), and load power of a typical sunny day, the optimal sizing result is proven effective and applicable. Besides, the optimal sizing result is also verified by the operation test of the sprinkler irrigation machine.



---

*Source: 6133969-2022-01-19.xml*