[GRASS-SVN] r66495 - grass-addons/grass7/raster/r.green/r.green.hydro/r.green.hydro.technical
svn_grass at osgeo.org
svn_grass at osgeo.org
Wed Oct 14 05:44:49 PDT 2015
Author: Giulia
Date: 2015-10-14 05:44:49 -0700 (Wed, 14 Oct 2015)
New Revision: 66495
Modified:
grass-addons/grass7/raster/r.green/r.green.hydro/r.green.hydro.technical/r.green.hydro.technical.html
grass-addons/grass7/raster/r.green/r.green.hydro/r.green.hydro.technical/r.green.hydro.technical.py
Log:
r.green: changes in gui and in the manual
Modified: grass-addons/grass7/raster/r.green/r.green.hydro/r.green.hydro.technical/r.green.hydro.technical.html
===================================================================
--- grass-addons/grass7/raster/r.green/r.green.hydro/r.green.hydro.technical/r.green.hydro.technical.html 2015-10-14 12:39:59 UTC (rev 66494)
+++ grass-addons/grass7/raster/r.green/r.green.hydro/r.green.hydro.technical/r.green.hydro.technical.html 2015-10-14 12:44:49 UTC (rev 66495)
@@ -4,7 +4,7 @@
<h2>NOTES</h2>
-Firstly, the module computes the <b>structure</b>. As the current potential concerns especially small hydropower (inferior to 20 MW), the structure suggested is the one for small hydropower detailed in the picture below. It is composed of an intake (A) which diverts water from the river. This water is conveyed into a derivation channel (B) with a very low slope and arrives in a forebay tank (C) which regulates the fluctuation of discharge. Finally, the penstock (D) conveys the water with the highest possible head to the turbine-alternator group (E) which produces electricity. The water is then released in the river (restitution F). We use the following vocabulary: the structure of the plant corresponds to the part with the derivation channel, the forebay tank and the penstock, whereas the segment of the plant corresponds to the part of the river (water not diverted) between the intake (A) and the restitution (F).<br><br>
+Firstly, the module computes the <b>structure</b>. As the current potential concerns especially small hydropower (inferior to 20 MW), the structure suggested is the one for small hydropower detailed in the picture below. It is composed of an intake (A) which diverts water from the river. This water is conveyed into a derivation channel (B) with a very low slope and arrives in a forebay tank (C) which regulates the fluctuation of discharge. Finally, the penstock (D) conveys the water with the highest possible head to the turbine-alternator group (E) which produces electricity. The water is then released into the river (restitution F). We use the following vocabulary: the structure of the plant corresponds to the part with the derivation channel, the forebay tank and the penstock, whereas the segment of the plant corresponds to the part of the river (water not diverted) between the intake (A) and the restitution (F).<br><br>
<center>
<img src="r_green_hydro_technical_picstruct.png" alt="structure"><br>
@@ -13,9 +13,9 @@
The power is maximized for the highest head in the penstock so the derivation channel is computed along the same quote (the low slope is neglected here) until the point which maximizes the head along the penstock. The structure is computed for both sides of the river in order to determine which one produces the most power.<br><br>
-Using the computed structure, the module calculates the <b>head losses</b> :<br><br>
+Using the computed structure, the module calculates the <b>head losses</b>:<br><br>
<blockquote>- in the derivation channel <br>
-<blockquote>There are regular losses calculated thanks to Manning's formula : <i>
+<blockquote>There are regular losses calculated thanks to Manning's formula: <i>
Δh<SUB>deriv</SUB>=L*(Q/(ks*A*Rh<SUP>2/3</SUP>))<SUP>2</SUP></i><br>
<blockquote>where Rh is the hydraulic radius (m),<br>
A the cross sectional area of flow (m<SUP>2</SUP>),<br>
@@ -24,29 +24,29 @@
ks the Strickler coefficient (m<SUP>1/3</SUP>/s), we consider steel as default parameter, with ks=75 m<SUP>1/3</SUP>/s.<br><br></blockquote></blockquote>
- in the forebay tank<br>
-<blockquote>There are singular losses caused by the change of section in the forebay tank and to the change of direction in the penstock (steep slope).<br>
+<blockquote>There are singular losses caused by the change of section in the forebay tank and the change of direction in the penstock (steep slope).<br>
-In any case, singular head losses are expressed like this : <i>Δh<SUB>sing</SUB>=K*V<SUP>2</SUP>/(2g)</i><br>
+In any case, singular head losses are expressed like this: <i>Δh<SUB>sing</SUB>=K*V<SUP>2</SUP>/(2g)</i><br>
<blockquote>where V is the velocity (m/s), <br>
g the gravity term (9,81 m/s<sup>2</sup>), <br>
K is a coefficient determined according to the kind of singularity.</blockquote>
-In our case, the singular losses are the sum of the ones for these three phenomena :<br>
-<blockquote>- enlargement at the entrance of the forebay tank : K<SUB>1</SUB>=1 and V=1 m/s<br>
-- narrowing at the exit of the forebay tank : K<SUB>2</SUB>=0.5 and V=4Q/(πD<SUP>2</SUP>) m/s<br>
-- bend at the beginning of the penstock : K<SUB>3</SUB>=(gross head/L)<SUP>2</SUP>+2*sin(ASIN(gross head/L)/2)<SUP>4</SUP> and V=4Q/(πD<SUP>2</SUP>) m/s
+In our case, the singular losses are the sum of the ones for these three phenomena:<br>
+<blockquote>- enlargement at the entrance of the forebay tank: K<SUB>1</SUB>=1 and V=1 m/s<br>
+- narrowing at the exit of the forebay tank: K<SUB>2</SUB>=0.5 and V=4Q/(πD<SUP>2</SUP>) m/s<br>
+- bend at the beginning of the penstock: K<SUB>3</SUB>=(gross head/L)<SUP>2</SUP>+2*sin(ASIN(gross head/L)/2)<SUP>4</SUP> and V=4Q/(πD<SUP>2</SUP>) m/s
<br><br></blockquote></blockquote>
- in the penstock<br>
-<blockquote>There are regular losses calculated thanks to this formula : <i>Δh<SUB>pen</SUB>=(f*8*L*Q<SUP>2</SUP>)/(π<SUP>2</SUP>*D<SUP>5</SUP>*g)</i><br>
-<blockquote>where L is penstock length (m), <br>
+<blockquote>There are regular losses calculated thanks to this formula: <i>Δh<SUB>pen</SUB>=(f*8*L*Q<SUP>2</SUP>)/(π<SUP>2</SUP>*D<SUP>5</SUP>*g)</i><br>
+<blockquote>where L is the penstock length (m), <br>
D is the penstock diameter (m), <br>
Q is the discharge (m<SUP>3</SUP>/s), <br>
-f is Darcy-Weisbach friction coefficient, which can be determined by Colebrooke formula. We consider steel by default with absolute roughness of ε = 0,015 mm.<br><br></blockquote></blockquote></blockquote>
+f is the Darcy-Weisbach friction coefficient, which can be determined by Colebrooke formula. We consider steel by default with absolute roughness of ε = 0,015 mm.<br><br></blockquote></blockquote></blockquote>
-Then, the module chooses the <b>turbine</b> which is the most accurate for each plant. The data of possible turbines are gathered in the folder turbine. For each turbine there is a file text with the ranges of discharge and head required to use it. There is also the efficiency in function of QW/Q_design. The turbine is designed to work at Q_design and QW corresponds to the real discharge which flows is the turbine. As we don't consider the duration curves but only the mean annual discharge, we assume that QW=Q_design.<br><br>
+Then, the module chooses the <b>turbine</b> which is the most accurate for each plant. The data of possible turbines of which you have to enter the path into the turbine_folder=string field, are gathered in the folder turbine. For each turbine there is a text file with the ranges of discharge and head required to use it. There is also the efficiency in function of QW/Q_design. The turbine is designed to work at Q_design and QW corresponds to the real discharge which flows in the turbine. As we don't consider the duration curves but only the mean annual discharge, we assume that QW=Q_design.<br><br>
-There is already a large choice of turbines available in turbine folder but the user is encouraged to create a file text with the turbines he wants to consider. The file has to have the same scheme with the correct information at the corresponding lines :<br>
+If you want to create an additional text file for another turbine model, the file has to have this scheme with the correct information at the corresponding lines:<br>
<blockquote><pre>TURBINE ALPHA_C
Name of the turbine Value of alpha_c
Q_MIN Q_MAX
@@ -56,9 +56,11 @@
QW/Q_design ETA
Coordinates of the curve efficiency=f(QW/Q_design)<br><br></pre></blockquote>
+In the turbine folder there is already the text file called list with a large choice of turbines available. You have to enter the path of the list file into the turbine_list=string field of the GUI. But the user can also create his own text file that has to have the same structure with a list of the names of the turbines he has selected and wants to be considered. <br><br>
+
To choose the turbine, the module first selects the turbines with ranges of discharge and head containing the values of the corresponding potential plant. Among these turbines, it chooses the one which has the best efficiency for QW=Q_design.<br><br>
-Thus the efficiency of the turbine is found. The global efficiency also includes the <b>efficiencies</b> of the shaft, the alternator and the transformer which can be chosen or are respectively equal to 1, 0.96 and 0.99 by default.<br><br>
+Thus the efficiency of the turbine is found. The global efficiency also includes the <b>efficiencies</b> of the shaft, the alternator and the transformer which can be chosen or they are respectively equal to 1, 0.96 and 0.99 by default.<br><br>
Finally, the corrected value of power which can be exploited is calculated.<br>
It corresponds to <i>P=η*ρ*g*Q*Δh<SUB>net</SUB></i><br>
@@ -68,15 +70,15 @@
Q the discharge (m<SUP>3</SUP>/s),<br>
Δh<SUB>net</SUB> the net head, that means the gross head minus head losses<br><br></blockquote>
-The output map of the module is the one with the structure for each plant, including in the <b>table</b> the data of :<br>
+The output map of the module is the one with the structure for each plant, including in the <b>table</b> the data of:<br>
<blockquote>- discharge (m<sup>3</sup>/s)<br>
- gross head (m)<br>
-- kind of the channel : derivation (conduct) or penstock<br>
+- kind of the channel: derivation (conduct) or penstock<br>
- side of the river (option0 or option1)<br>
- diameter of the channel (m)<br>
- losses in the channel (m)<br><br>
-Moreover, only in the penstock's line of the structure, there are :<br>
+Moreover, only in the penstock's line of the structure, there are:<br>
- singular losses (m) in the forebay tank between the derivation channel and the penstock<br>
- the total losses (m) for each structure, which are the sum of the regular losses in the derivation channel and in the penstock and the singular losses in the forebay tank<br>
- net head (m), which is the gross head minus the total losses<br>
@@ -85,22 +87,22 @@
- kind of the selected turbine (turbine)<br>
- power (kW) which can be exploited considering the technical constrains<br>
- global efficiency (power/hyd_power)<br>
-- max_power : yes or no, yes for the side (option1 or option0) which produces the most power</blockquote>
+- max_power: yes or no, yes for the side (option1 or option0) which produces the most power</blockquote>
<h2>EXAMPLE</h2>
This example is based on the case-study of Gesso and Vermenagna valleys in the Natural Park of the Maritime Alps, Piedmont, Italy.<br><br>
-Here is the input vector map potentialplants with the intakes and restitutions (in red) computed by r.green.hydro.recommended. The vector map with the segments of river is also visibile in blue on this picture. These potential plants have a maximum length of 800m and a distance of 800m between them.
+Here is the input vector map potentialplants with the intakes and restitutions (in red) computed by r.green.hydro.recommended. The vector map with the segments of river is also visible in blue in this picture. These potential plants have a maximum length of 800 m and a distance of 800 m between them.
<center>
<img src="r_green_hydro_technical_input.png" alt="input"><br>
Potential intakes and restitutions
</center><br><br>
-The following command computes the structure of the potential plants for each side of the river and include the corrected power on the output map table :<br>
-<div class="code"><pre>r.green.hydro.technical plant=potentialplants elevation=elevation structplants=techplants</pre></div><br>
-
+The following command that you can either put in the command console or the GUI of r.green.hydro.technical, computes the structure of the potential plants for each side of the river and includes the corrected power on the output map table:<br>
+<div class="code"><pre>r.green.hydro.technical plant=potentialplants elevation=elevation output_struct=techplants output_plant=segmentplants
+turbine_folder=/pathtothefileoftheturbine_folder turbine_list=/pathtothefileoftheturbine_list</pre></div><br>
The result is shown in the following vector map called techplants. The table of this map is completed as explained in the end of the NOTES part.<br><br>
<center>
Modified: grass-addons/grass7/raster/r.green/r.green.hydro/r.green.hydro.technical/r.green.hydro.technical.py
===================================================================
--- grass-addons/grass7/raster/r.green/r.green.hydro/r.green.hydro.technical/r.green.hydro.technical.py 2015-10-14 12:39:59 UTC (rev 66494)
+++ grass-addons/grass7/raster/r.green/r.green.hydro/r.green.hydro.technical/r.green.hydro.technical.py 2015-10-14 12:44:49 UTC (rev 66495)
@@ -143,16 +143,16 @@
#%option
#% key: turbine_folder
#% type: string
-#% description: Path to the folder contening the text file with info about all kind of turbines
-#% required: no
+#% description: Path to the folder containing the text file with info about all kind of turbines
+#% required: yes
#% answer:
#% guisection: Turbine
#%end
#%option
#% key: turbine_list
#% type: string
-#% description: Path to the text file contening the list of the turbines considered
-#% required: no
+#% description: Path to the text file containing the list of the turbines considered
+#% required: yes
#% answer:
#% guisection: Turbine
#%end
More information about the grass-commit
mailing list