水平な不透水層上に集水埋渠を設け伏流水がその両側から流入する場合、CALCULATION OF WATER VOLUME FLOWING IN DRAIN. CASE:RIVER-BED WATER,WATER COLLECTING DRAIN.

Equation Q=\frac{kL(H^{2}-h^{2})}{R}000183Scilab 183sc01183sc02183sc03183sc04
 -->clear;
-->funcprot(0);
-->deff('[Q]=f(R,h)','Q=k/100*L*(H^2-h^2)/R');//地下水揚水量:Q(m^3/sec)
-->k=10^-1;//透水係数(cm/sec)
-->H=30;//(m)
-->L=20;//(m)
-->R=linspace(5,20);//(m)
-->h=linspace(0.5,3);//(m)
-->xset("window",1);
-->clf();
-->xgrid();
-->condition='k='+string(k)+',H='+string(H);
-->graphtitle='地下水揚水量Q(m^3/sec).' + condition;
-->xtitle(graphtitle,'R(m)','h(m)','Q(m^3/s)',boxed=1);
-->color(240,255,240);
-->idcolor=color(240,255,240);
-->fplot3d(R,h,f,flag=[idcolor , 2 , 4]);
-->xset("window",2);
-->clf();
-->xgrid();
-->xtitle(graphtitle,'R(m)','h(m)','Q(m^3/s)',boxed=1);
-->contour(R,h,f,30,flag=[2,2,4]);
-->f(5,3)
 ans  =
 
    3.564  
-->f(20,0.5)
 ans  =
 
    0.89975  
-->clear;
-->funcprot(0); 
-->deff('[Q]=f(R,H)','Q=k/100*L*(H^2-h^2)/R');
-->//地下水揚水量:Q(m^3/sec)
-->k=10^-1;//透水係数(cm/sec)
-->h=3;//(m)
-->L=20;//(m)
-->R=linspace(5,20);//(m)
-->H=linspace(4,10);//(m)
-->xset("window",3);
-->clf();
-->xgrid();
-->condition='k='+string(k)+',h='+string(h);
-->graphtitle='地下水揚水量Q(m^3/sec).' + condition;
-->xtitle(graphtitle,'R(m)','H(m)','Q(m^3/s)',boxed=1);
-->color(255,239,219);
-->idcolor=color(255,239,219);
-->fplot3d(R,H,f,flag=[idcolor , 2 , 4]);
-->xset("window",4);
-->clf();
-->xgrid();
-->xtitle(graphtitle,'R(m)','H(m)','Q(m^3/s)',boxed=1);
-->contour(R,H,f,30,flag=[2,2,4]);
-->f(15,7)
 ans  =
 
    0.0533333  
SciNotes
clear;
funcprot(0);
deff('[Q]=f(R,h)','Q=k/100*L*(H^2-h^2)/R');//地下水揚水量:Q(m^3/sec)
k=10^-1;//透水係数(cm/sec)
H=30;//(m)
L=20;//(m)
R=linspace(5,20);//(m)
h=linspace(0.5,3);//(m)
 
xset("window",1);
clf();
xgrid();
condition='k='+string(k)+',H='+string(H);
graphtitle='地下水揚水量Q(m^3/sec).' + condition;
xtitle(graphtitle,'R(m)','h(m)','Q(m^3/s)',boxed=1);
color(240,255,240);
idcolor=color(240,255,240);
fplot3d(R,h,f,flag=[idcolor , 2 , 4]);
 
xset("window",2);
clf();
xgrid();
xtitle(graphtitle,'R(m)','h(m)','Q(m^3/s)',boxed=1);
contour(R,h,f,30,flag=[2,2,4]);
 
f(5,3)
 
f(20,0.5)
   
clear;
funcprot(0); 
deff('[Q]=f(R,H)','Q=k/100*L*(H^2-h^2)/R');
//地下水揚水量:Q(m^3/sec)
k=10^-1;//透水係数(cm/sec)
h=3;//(m)
L=20;//(m)
R=linspace(5,20);//(m)
H=linspace(4,10);//(m)
 
xset("window",3);
clf();
xgrid();
condition='k='+string(k)+',h='+string(h);
graphtitle='地下水揚水量Q(m^3/sec).' + condition;
xtitle(graphtitle,'R(m)','H(m)','Q(m^3/s)',boxed=1);
color(255,239,219);
idcolor=color(255,239,219);
fplot3d(R,H,f,flag=[idcolor , 2 , 4]);
 
xset("window",4);
clf();
xgrid();
xtitle(graphtitle,'R(m)','H(m)','Q(m^3/s)',boxed=1);
contour(R,H,f,30,flag=[2,2,4]);
 
f(15,7)
参考文献 川島普、篠原紀、西川泰治(1999).『新訂衛生工学』.明現社.184pp.アプリケーション Scilab http://www.scilab.org/