Method for estimating optimal efficiency point parameters and performance curve in axial-flow PAT power generation mode

This application is a national stage application of PCT/CN2018/090714 filed Jun. 11, 2018, which claims the benefit of priority from Chinese Application No. 201810037215.9 filed Jan. 15, 2018, the disclosures of both applications being incorporated herein by reference in their entirety.

The present invention belongs to the field of micro-hydropower technology, and specifically relates to a method for estimating the optimal efficiency point parameters and performance curve in an axial-flow Pump as turbine (PAT) power generation mode based on related parameters.

Based on data of the international energy agency, the total primary energy consumption is increasing rapidly in recent decades, and the gap of energy supply and demand is increasing day by day. For national and social sustainable development, it is very necessary to develop and utilize the renewable energy resource, where the waterpower resource as a green and environmentally friendly energy source still has a huge technical development potential in the world. It can be seen from the hydropower development process: the current development level of large and medium-sized hydropower stations is high and tend to be stable in developed nations and partial newly-emerging nations, and development of the hydropower resource is gradually transformed from the large-medium-sized hydropower stations to small-micro hydropower stations.

A micro hydropower station is a hydropower system with an installed capacity less than 100 kW, and can be flexibly established in any water system with surplus energy according to the demand. PAT is a currently most-effective manner to reduce all costs and shorten the investment payback period, which has attracted the attention of many international scholars. Compared with the traditional turbine, the PAT has the characteristics of low price, high standardization and serialization, simple and convenient installation, and easy accessibility of spare parts. These characteristics reduce the comprehensive cost of the PAT and make up for the disadvantage of low efficiency of PAT compared with the traditional turbine to a certain extent. However, manufacturers generally do not test performance of the PAT in a power generation mode, so an optimal efficiency point and a characteristic curve of the PAT are unknown, which causes model selection difficulty during actual application.

Currently, the prediction of hydropower characteristic of an axial-flow PAT power generation mode mainly depends on model test or three-dimensional CFD numerical simulation. These two methods have high reliability and accuracy, but their disadvantages are also obvious. The model test needs a large amount of manpower and resources, and the debugging and installing period of a test platform and test equipment is long. The numerical simulation has high requirements on theory and experience of engineers, and takes a long time. Therefore, these two methods are difficult to be popularized and applied in engineering practice.

In conclusion, there is a lack of a simple, convenient and feasible method for estimating the optimal efficiency point parameters and performance curve in an axial-flow PAT power generation mode to meet the needs of unit selection in engineering practice.

In order to solve the above problem, the present invention is to provide a method for estimating optimal efficiency point parameters and performance curve in an axial-flow PAT power generation mode, which is simple, convenient and feasible, has a short calculation period, and also has great prediction accuracy in a normal operation range of a PAT.

To achieve the above objective, the present invention adopts the following solutions:

Optimal efficiency point parameters estimation method

The present invention provides a method for estimating optimal efficiency point parameters in an axial-flow PAT power generation mode, where the method includes the following steps: I1, calculating an axial velocity ν ml,t of an optimal efficiency point in a power generation mode by utilizing the following formula:

where ν ml,t is an axial component of an absolute velocity (m/s), Q pb is an flow rate of an optimal efficiency point in a pumping mode (m 3 /s), R is a radius of an airfoil position (m), n is a rotational velocity of an impeller (r/min), D is impeller diameter (m), d h is a diameter of a hub (m), β e is an impeller airfoil setting angle)(°, and β e,g is a guide blade setting angle (°);

I2, calculating a flow rate Q tb of the optimal efficiency point in the power generation mode based on the axial velocity and flow area of the optimal efficiency point:

where Q tb is the flow rate of the optimal efficiency point in the power generation mode (m 3 /s);

I3, calculating a theoretical hydraulic head H t,th of an axial-flow PAT:

where u l,t a circumferential velocity (m/s), ν ul,t is a circumferential component of the absolute velocity (m/s), and g is a gravitational acceleration (m/s 2 );

I4, segmenting a whole flow channel of the axial-flow PAT according to geometrical characteristics, and respectively calculating a frictional hydraulic head loss and a local hydraulic head loss of each segment according to formulas listed in Table 1:

where λ f,e is a hydraulic frictional resistance coefficient of the elbow pipe, θ e is a bend angle of the elbow pipe (°), r e is ea radius of the elbow center line (m), r 0,e is a radius of the elbow pipe (m), ν e is an average flow velocity of the elbow pipe (m/s), f f,s is a hydraulic frictional resistance coefficient, θ s is a convergence angle (°), A s1 and A s2 are a flow area of an inlet of the reductive pipe and a flow area of an outlet of the reductive pipe (m 2 ) respectively, ν s1 and ν s2 respectively are an average flow velocity of the inlet of the reductive pipe and an average flow velocity of the outlet of the reductive pipe (m/s), ξ s is a local loss coefficient and is determined by the size of the convergence angle, h f is a frictional hydraulic head loss (friction loss) in an impeller flow channel (m), Z is the number of impeller blades, λ is a hydraulic frictional resistance coefficient, l b is a length of the impeller flow channel and can be calculated approximately according to an airfoil chord length (m), D hyd is an equivalent hydraulic diameter of the impeller flow channel (m), is an average relative velocity of the impeller flow channel:

w _ = ( w 1 ′ , p 2 + w 2 ′ , p 2 ) 2 ,

w 1′,p is a relative velocity of a blade inlet (m/s), w 2′,p is a relative velocity of a blade outlet (m/s), h shock is a fluid shock loss of the blade inlet (m), w i is a fluid shock velocity of the blade inlet, that is a component of the relative velocity in an airfoil vertical direction (m/s), h wake is a wake loss at the blade outlet (m), ξ is an excretion coefficient at the impeller outlet, subscripts g of Z g , λ g , l b,g , D hyd,g and w g represents the guide blade zone, the meaning of each symbol is consistent with the symbols of the impeller zone, λ f,d is a hydraulic frictional resistance coefficient of the diffuser, θ d is a diffusion angle of the diffuser (°), A d1 and A d2 respectively are a flow area of a diffuser inlet and a flow area of a diffuser outlet (m 2 ), ν d1 and ν d1 respectively are an average flow velocity of the diffuser inlet and an average flow velocity of the diffuser outlet (m/s), k is a diffuser coefficient and is determined by the size of the diffuser;

I5, calculating an output power P tb of the optimal efficiency point in the power generation mode based on the theoretical hydraulic head and the flow of the optimal efficiency point in the power generation mode:

P tb =ρgQ tb H t,th ,

I6, calculating an actual hydraulic head H tb of the axial-flow PAT at the optimal efficiency point in the power generation mode:

H tb =H t,th +Σh t,loss ,

where Σh t,loss is a sum of all hydraulic head losses; and

I7, calculating an optimal efficiency η tb based on a ratio of the theoretical hydraulic head to a total hydraulic head:

Performance Curve Estimation Method

The present invention further provides a method for estimating a performance curve in an axial-flow PAT power generation mode, where the performance curve is estimated based on the above described estimation method of an optimal efficiency point parameter in an axial-flow PAT power generation mode, and the method includes the following steps:

II1, calculating a normalized flow-hydraulic head curve based on similarity hypothesis of an axial-flow PAT normalized performance curve:

h t =2.55 q t 2 −1.30 q t −0.25,

where h t is a normalized hydraulic head, and calculated according to

h t = H t H tb ,

H t is a hydraulic head, H tb is an optimal efficiency point hydraulic head, and q t is a normalized flow rate; according to

q t = Q t Q tb ,

Q t is flow rate, and Q tb is the optimal efficiency point flow rate; and a normalized parameter is defined as be a ratio of an actual parameter value to an optimal efficiency point parameter value;

II2, calculating a normalized hydraulic head-output power curve:

p t =1.27 h t −0.27,

where p t is a normalized output power, and according to

p t = P t P tb ,

P t is the output power, and P tb is an optimal efficiency point output power; and

II3, calculating a hydraulic head-efficiency curve:

where η tb is optimal efficiency.

The present invention obtains all performance curves, including a head-flow curve, an output power-head curve and an efficiency-head curve, of a specific axial-flow PAT by defining the normalized parameter and performing a reversal treatment on a normalized curve.

In the calculation of the parameters of the optimal efficiency point, the flow rate can be calculated by using the binary theory to analyze and solve the fluid motion, that is, a stream surface of 50% blade-height position is selected to be analyzed; or the fluid movement is analyzed and solved by utilizing a quasi-ternary theory, that is, stream surfaces of N blade-height positions are selected to be analyzed and solved, and the solved N axial velocities conduct weighted average on the overflow area according to their radius positions in order to obtain an average axial velocity.

In the calculation of hydraulic head, a friction coefficient of the frictional hydraulic head loss is calculated according to the Haaland formula, where the wall roughness is selected according to pipe conditions. To an actual PAT unit, the roughness of its blade surfaces and guide blade surfaces generally can be determined according to industrial standards, however, there is no uniform processing standard for the roughness of other parts of the supporting pipeline system, so determination and trial are performed in advance according to the material of the pipeline and the utilization condition, correction is performed by using an optimal efficiency point parameter in a pumping model, and actual roughness is obtained based on multiple trials.

Function and Effect of the Present Invention

Firstly, the optimal efficiency point parameter is the most important fundamental data in all kinds of PAT applications, and is the basis of unit selection in the planning and designing stage. The estimation method provided by the present invention can directly predict each energy parameter of a PAT unit at an operating point. In the existing method, the energy parameters of the optimal efficiency point can be further obtained based on the obtained complete performance curve.

Secondly, the estimation method provided by the invention is based on similarity of the normalized performance curve of the axial-flow PAT, when predicting multiple operating conditions, the method is fast, simple and easy to use, and can flexibly select a predicted operating condition range according to hydrological characteristics.

Thirdly, when the feasibility is also considered, the estimation method provided by the present invention has a great compatibility and can be applied to an axial-flow PAT with the specific velocity 120-300 (the specific velocity is the specific velocity of the optimal efficiency point of the pumping mode, and the formula is

n ⁢ Q H 3 / 4 ) .

The present invention has high adaptability and can provide reliable accuracy in an allowed axial-flow PAT operating range.

FIG. 1 is a velocity triangle of an inlet and an outlet of an optimal efficiency point in a pumping mode, u represents a circumferential velocity, w represents a relative velocity, ν represents an absolute velocity, A represents an airfoil setting angle, a subscript 1 represents an inlet while a subscript 2 represents an outlet, and p represents an operating condition of the pumping mode;

FIG. 2 is a velocity triangle of an inlet and an outlet of an optimal efficiency point in a power generation mode, u represents a circumferential velocity, w represents a relative velocity, ν represents an absolute velocity, A represents an airfoil setting angle, a subscript 1 represents an inlet while a subscript 2 represents an outlet, and t represents an operating condition of the pumping mode;

FIG. 3 is a schematic structural diagram of a flow channel of an axial-flow PAT unit;

FIG. 4 is a schematic diagram of a segment of an axial-flow PAT unit, where (A), (B) and (D) are side views, and (C) is an axial view; that is, an impeller zone is divided into several blade channels according to the number of blades (there are three blades in the figure, so there are three blade channels), parameter calculation of its hydraulic diameter is conducted according to a single blade channel, the final friction loss of the single blade channel is multiplied by the number of the blades so as to obtain a total friction loss of the impeller zone, and calculation method of a friction loss in the guide vane zone is similar;

FIG. 5 is a schematic diagram showing decomposition of shock velocity of a blade inlet, that is, a relative velocity w of a fluid is orthogonally decomposed into two components in parallel with a wing chord and in a direction vertical to the wing chord, where a component wi in the direction vertical to the wing chord is the shock velocity;

FIG. 6 is a three-dimensional model of an axial-flow PAT; and

FIG. 7 , FIG. 8 and FIG. 9 are comparison diagrams of a performance curve obtained by the estimation method provided by the present invention and original data, where FIG. 7 is a comparison diagram of a hydraulic head-flow curve, FIG. 8 is a comparison diagram of a hydraulic head-output power curve, and FIG. 9 is a comparison diagram of a hydraulic head-power curve.

The following describes specific embodiments of a method for estimating an optimal efficiency point parameter and performance curve in an axial-flow PAT power generation mode according to the present invention in detail with reference to the accompanying drawings.

As shown in FIG. 6 , the rotational velocity of an impeller of an axial-flow PAT n is equal to 1450 r/min, the diameter of the impeller D is equal to 0.30 m, the diameter of the hub d h is equal to 0.108 m, an airfoil setting angle of the impeller A is equal to 23°, a setting angle of a guide blade β e.g is equal to 80°, and a radius R of 50% blade-height airfoil is equal to 0.105 m. Energy parameters of an optimal efficiency point in a pumping mode are: the hydraulic head is 3.33 m, the flow rate is 0.326 m 3 /s, the shaft power is 12.96 kW and the hydraulic efficiency is 82.09%. Energy parameters of an optimal efficiency point in a power generation mode are: the hydraulic head is 4.47 m, the flow rate is 0.457 m 3 /s, the output power is 16.04 kW and the hydraulic efficiency is 80.09%.

Prediction results by using an optimal efficiency point estimation method provided by the present invention are: the flow rate is 0.471 m 3 /s (corresponding to the steps I1 and I2), the hydraulic head is 4.29 m (corresponding to the steps I3, I4 and I6), the output power is 15.49 kW (corresponding to the step I5), and the hydraulic efficiency is 78.11% (corresponding to the step I7). Prediction error of each parameter is as follows: the prediction error of the flow rate is 3.06%, the prediction error of the hydraulic head is 4.03%, the prediction error of the hydraulic efficiency is 1.98%, and the prediction error of the output power is 3.43%.

The method provided by the present invention is used for predicting a performance curve of a certain axial-flow PAT (corresponding to the steps II1 to II3), and compares the obtained performance curve with the original data, as shown in FIG. 7 , and it can be seen that the prediction results of the present invention is consistent with the trend of the original data and is in good agreement. Within the allowable operation range, the average error of the hydraulic head is 5.72%, the average error of the output power is 4.68% and the average error of the hydraulic efficiency is 2.42%.

The above embodiments merely illustrate the technical solution of the present invention. The method for estimating an optimal efficiency point parameter and performance curve in an axial-flow PAT power generation mode associated with the present invention is not only limited to contents described in the above embodiments, but also is subject to a scope defined in the claims. Any modifications, supplementations or equivalent replacements made by a person skilled in the part based on the embodiments shall fall within the protection scope of the present invention.