哲学者の朝食はタバコとコーヒー

*フロイト(病気)とボルツマン(気象)の描く千佳画*//

rain.18.04.02.01.ー降水と大気は量子効果を持つ(4)

容易に想像つくことだが、降水の日数と大気の日数は、全体では1:2であるが、降水が自由度1、大氣が自由度2の降水量量分布関数であるから、降水量に対しては、傾きが違っている。従って、ある降水量以上の降水と大気ではその比率は1:2とはなっていない。
降水と大気の降水量が同じだから、ある降水量以上を分布関数に対して積分すれその比率は容易に求まる。問題は検証だ。大気の降水量をどうやって計るかで、ある。計れなければ、ある降水量以上の大気の個数が分からない。

いま、いったん降った雨をすべて貯留しよう。それを一定量ずつ毎日放流しよう。放流量をzとする。雨が集中すれば貯留地の水量は増し、晴れが集中すれば貯留池の水量は減っていく。雨の日数と晴れの日数の比、P(z)とする、によって、その貯留量は決定され、ある比率で最大になるだろう。最大の貯留容量をVとする。雨も晴れも集中の度合いには、全体で1:2になるから、上限があり、貯留量Vが最大になる比率P(z)は、放流量zによって決まってくることになる。
降水の降水量分布関数を自由度1のボルツマン分布、大氣の降水量分布関数を自由度2のボルツマン分布、指数分布と一般には呼ばれる、とする。ここで古典統計のボルツマン分布としたのは、量子統計の量子効果を考えないためであるし、地域の観測所の降水量分布は古典統計を示すためである。放流量zでは、それ以上の降水は貯留され、また雨の日と晴れの日に放流される。分布関数と置いたから、降水と大気の降水量はランダムに出現し、また、降水と大気の日数の比率P(z)は全体も部分も同じになっていると考えられる。
放流量zを変えていけば、最大貯留量Vも比率P(z)も変わって来るだろう。実際、1年間の降水を時系列に従って、出し入れ計算すると、放流量zと最大貯留容量Vの間に、右下がりのQ-V曲線が得られる。ここで、Qは年間降水量である。このQ-V曲線を大気と降水の分布関数から理論的に求めることができれば、降水と大気の降水量xに対する比率P(z)が求まることになる。降水量x1m^2に降った降水に対してであるが、放流量zは例えば、1haに降った降水に対してになる。

最大貯留量をV、放流量z以上の貯留量をV_1、雨の日の放流量をV_2、晴れの日の放流量をV_3とすれば、最大貯留量Vは、
         V = V_1 -V_2-V_3
放流量zに対して雨の日に従う分布関数をf(z)_1、自由度1のボルツマン分布である、晴れの日の分布関数をf(z)_2、自由度2のボルツマン分布である、とする。添え字は自由度である。雨の日数をA_1 、晴れの日数をA_2とし、分布関数の平均をそれぞれδ_1 /2δ_2とすると、
         V_1 = A_1\int_z^∞ zf(z)_1 dz
         V_2 = zA_1\int_z^∞ f(z)_1 dz
         V_3 = zA_2\int_z^∞ f(z)_2 dz

         f(z)_1 = \displaystyle\frac{1}{\sqrt{π}}
                                                               \displaystyle\frac{1}{\sqrt{δ_1}}z^{-1/2} e^{-z/δ_1}
         f(z)_2 =\displaystyle\frac{1}{δ_2} e^{-z/δ_2}
         δ_1 = 2*δ_2

雨の日の日数A_1と晴れの日の日数A_2の比をP(z)とおこう。
         P(z)=\displaystyle\frac{A_2}{A_1}
雨の年間降水量Qを基準としよう。
         Q = \displaystyle\frac{δ_1}{2}A_1
         Q =\langle δ_2\rangle\displaystyle\frac{A_2}{2}

         \langle δ_2\rangle =
                                              \displaystyle\frac{\int_z^∞ zf(z)_2 dz}{\int_z^∞ f(z)_2 dz}
従って、
         P(z) = \displaystyle\frac{δ_1}{\langle δ_2\rangle}
これを0から∞まで積分すると2である。

降水量xから放流量zを求めよう。

埋立地B m^2の最終処分場を考えよう。埋立地に降った雨はゴミの中を通って出てくる。浸出水という。今降水量をxとし、月別、地域別の進出係数をq_iとすれば、浸出水量zBを埋め立て地面積とすれば、
         z_i = xq_i B
これを次の自由度1のボルツマン分布、
         f(x) =\displaystyle\frac{1}{\sqrt{π}\sqrt{λ}}x^{-1/2} e^{-x/λ}
に、代入すると、
         f(z_i) = \displaystyle\frac{1}{\sqrt{π}}
                                                               \displaystyle\frac{\sqrt{q_i B λ}}{\sum {q_i Bλ}}
                                                                     z_i ^{-1/2} e^{-z_i/q_i Bλ}
平均は、
         \langle z\rangle =\sum \int_0^∞ z_i f(z_i)dz
                                                              =\displaystyle\frac{1}{2}
                                                                           \displaystyle\frac{\sum (q_i Bλ)^2}{\sum q_i Bλ}
従って、\langle z\rangle = δ_1 /2より、
         f(z)_1 = \displaystyle\frac{1}{\sqrt{π}\sqrt{δ_1}}z^{-1/2}e^{-z/δ_1}
         δ_1 = \displaystyle\frac{\sum q_i^2}{\sum q_i }Bλ
同じように、
         f(z)_2 =\displaystyle\frac{1}{δ_2}e^{-z/δ_2}
         δ_2 = δ_1 /2

検証しよう。。

検証は豪雨地域の、高知県南国市御免観測所の過去15年間について行った。表6.と図6.に過去15年間の日降水量分布の平均値と理論値を示した。
理論値と観測値がよく一致している。日降水量分布は自由度1のボルツマン分布と置いてよいことが分かる。

f:id:silverforce69:20180403232153j:plain

f:id:silverforce69:20180403232240j:plain

表7.と図7-1.にQ-V曲線の、1年ごとに時系列に従って出し入れ計算し、15年間の平均値を取ったものと理論値を示した。
理論値と観測値がよく一致している。日降水量分布は、降水は自由度1の、大気は自由度2のボルツマン分布と置いてよいこと、それらの平均降水量は等しいこと、が分かる。

f:id:silverforce69:20180403233234j:plain

f:id:silverforce69:20180403233258j:plain

図7-2.と図7-3.に、V_1V_2V_3P(z)を示した。V_1V_2P(z)は緩やかに減少、または増加するが、V_3には最大値があることが分かる。そのため、V_3の最大値以上では放流量zを増加しても、調整池容量Vはあまり減少しなくなる。また、Q-V曲線の理論値と観測値が、この最大値付近で若干ズレており、理論値が若干小さく計算されていることが分かる。その他、降水量が大きくなれば、雨の日の割合が多くなっていくとが分かる。

f:id:silverforce69:20180403233321j:plain

f:id:silverforce69:20180403233339j:plain