JP3602951B2 - 環境状態測定方法及び装置 - Google Patents
環境状態測定方法及び装置 Download PDFInfo
- Publication number
- JP3602951B2 JP3602951B2 JP34528397A JP34528397A JP3602951B2 JP 3602951 B2 JP3602951 B2 JP 3602951B2 JP 34528397 A JP34528397 A JP 34528397A JP 34528397 A JP34528397 A JP 34528397A JP 3602951 B2 JP3602951 B2 JP 3602951B2
- Authority
- JP
- Japan
- Prior art keywords
- sound wave
- temperature
- measured
- space
- wall
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims description 43
- 230000007613 environmental effect Effects 0.000 title claims description 39
- 238000009826 distribution Methods 0.000 claims description 214
- 238000005259 measurement Methods 0.000 claims description 68
- 230000002093 peripheral effect Effects 0.000 claims description 17
- 238000005192 partition Methods 0.000 claims description 6
- 238000000638 solvent extraction Methods 0.000 claims description 4
- 238000000691 measurement method Methods 0.000 claims description 2
- 238000001514 detection method Methods 0.000 description 18
- 230000000694 effects Effects 0.000 description 16
- 238000010586 diagram Methods 0.000 description 15
- 230000005540 biological transmission Effects 0.000 description 14
- 230000001902 propagating effect Effects 0.000 description 6
- 238000004378 air conditioning Methods 0.000 description 4
- 238000006243 chemical reaction Methods 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 230000006870 function Effects 0.000 description 3
- 238000009529 body temperature measurement Methods 0.000 description 2
- 238000002591 computed tomography Methods 0.000 description 2
- 230000007423 decrease Effects 0.000 description 2
- 230000014509 gene expression Effects 0.000 description 2
- 238000010521 absorption reaction Methods 0.000 description 1
- 230000036760 body temperature Effects 0.000 description 1
- 239000003814 drug Substances 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01K—MEASURING TEMPERATURE; MEASURING QUANTITY OF HEAT; THERMALLY-SENSITIVE ELEMENTS NOT OTHERWISE PROVIDED FOR
- G01K11/00—Measuring temperature based upon physical or chemical changes not covered by groups G01K3/00, G01K5/00, G01K7/00 or G01K9/00
- G01K11/22—Measuring temperature based upon physical or chemical changes not covered by groups G01K3/00, G01K5/00, G01K7/00 or G01K9/00 using measurement of acoustic effects
- G01K11/24—Measuring temperature based upon physical or chemical changes not covered by groups G01K3/00, G01K5/00, G01K7/00 or G01K9/00 using measurement of acoustic effects of the velocity of propagation of sound
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2291/00—Indexing codes associated with group G01N29/00
- G01N2291/02—Indexing codes associated with the analysed material
- G01N2291/028—Material parameters
- G01N2291/02845—Humidity, wetness
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2291/00—Indexing codes associated with group G01N29/00
- G01N2291/02—Indexing codes associated with the analysed material
- G01N2291/028—Material parameters
- G01N2291/02881—Temperature
Landscapes
- Physics & Mathematics (AREA)
- Acoustics & Sound (AREA)
- General Physics & Mathematics (AREA)
- Measuring Temperature Or Quantity Of Heat (AREA)
- Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
- Air Conditioning Control Device (AREA)
Description
【０００１】
【発明の属する技術分野】
本発明は環境状態測定方法及び環境状態測定装置に係り、特に、空間内の温度分布や湿度分布等の環境状態を測定する環境状態測定方法、及び該環境状態測定方法を適用可能な環境状態測定装置に関する。
【０００２】
【従来の技術】
被空調空間内の温度や湿度を所望の値に制御するためには、被空調空間内の温度、湿度を測定する必要があるが、被空調空間内の各箇所における温度、湿度が一定であることは稀であり、特にドーム球場や劇場等の閉鎖された大空間では、空間内の各箇所における温度、湿度は大きくばらついていることが多い。このため、例えばドーム球場や劇場等の大空間において快適性と省エネルギーの双方を満足する空気調和を行う場合や、室内の温度や湿度を精密に所望の値に制御する場合には、被空調空間内の温度や湿度の分布を測定する必要がある。
【０００３】
しかし、温度センサや湿度センサによって温度や湿度の分布を直接測定しようとすると、センサを、被空調空間の中央部付近を含む被空調空間内にマトリクス状に多数配設する必要があるが、コストが嵩み景観が損なわれると共に被測定空間の利用も制限されるので、上記のように多数のセンサを配設することは現実的ではない。またドーム球場や劇場等では、被測定空間の周縁部（例えば壁体等）に取付けたセンサにより周縁部の温度等を測定すると共に、風船等でセンサを吊り上げることで被測定空間の中央部付近の温度等を測定することが考えられるが、この方法では、被測定空間内で例えば野球等のイベントが実施されている最中に温度や湿度の分布を測定することは不可能である。
【０００４】
一方、医学の分野では、人体に対して多方向からＸ線を照射し、人体を透過したＸ線量を高感度のセンサで測定し、測定されたＸ線量に基づき人体の横断面に沿った各位置における体内組織のＸ線吸収量をコンピュータにより演算して２次元画像として再構成し、人体の横断面像としてモニタに表示するコンピュータ断層撮影（ＣＴ：Ｃｏｍｐｕｔｅｄ Ｔｏｍｏｇｒａｐｈｙ）技術が広く利用されている。また、Ｘ線に代えて超音波を用いることも知られている。
【０００５】
上記のＣＴ技術によれば、測定対象（この場合は人体）内部の状態を非接触で測定することができるので、このＣＴ技術を利用して被測定空間内の温度分布を求めることが提案されている。具体的には、Ｘ線に代えて、媒質の温度によって伝搬時間（速度）が変化する音波を用い、例えば被測定空間を区画している壁体の多数箇所にスピーカ、マイクロフォン、温度センサを設け、スピーカ及びマイクロフォンによって前記多数箇所の間の音波の伝搬時間を各々測定し、被測定空間内の多数箇所の間の音波の伝搬時間から、被測定空間内の温度分布を演算によって求めている。
【０００６】
【発明が解決しようとする課題】
しかしながら、上記のＣＴ技術を利用した温度分布の測定は、音波の伝搬時間のみから温度分布を求めているので、実際の温度分布との誤差が大きい、という問題があった。実際の温度分布との誤差を小さくするためには、例えば被測定空間の周縁部に温度センサを設け、被測定空間内の多数箇所の間の音波の伝搬時間に加えて、温度センサによって測定した被測定空間の周縁部の温度も用いて被測定空間の温度分布を求めることが考えられる。
【０００７】
しかし、被測定空間が壁体によって区画されている場合、該壁体の温度は被測定空間内の温度と必ずしも一致しておらず、被測定空間内のうち壁体近傍の温度は、壁体の温度の影響を受けて大きく変化していることが一般的である。このため、被測定空間の壁体の温度を用いて被測定空間内の温度分布を求めたとしても、被測定空間内の温度分布の測定精度の向上には繋がらない、という問題があった。
【０００８】
また、上記のＣＴ技術を利用した温度分布の測定では、音波が直線的に伝搬するものとして温度分布を求めているが、音波はＸ線とは異なり直線的に伝搬するとは限らず、伝搬経路上に存在している空気の温度の影響を受けて音波の伝搬経路自体も変化し、これに伴って、スピーカ及びマイクロフォンによって測定される音波の伝搬時間も変化する。従って、被測定空間内の温度分布に依存する音波の伝搬経路の変化も、被測定空間内の温度分布の測定精度低下の一因となっていた。
【０００９】
更に、湿度分布については測定方法が確立されておらず、被測定空間内の中央部付近に湿度センサ等を設けて中央部付近の湿度を測定することなく、被測定空間内の湿度分布を求めることは不可能であった。
【００１０】
本発明は上記事実を考慮して成されたもので、被測定空間内の温度分布を高精度に求めることができる環境状態測定方法及び環境状態測定装置を得ることが目的である。
【００１１】
また本発明は、被測定空間内の湿度分布を、被測定空間の中央部付近に湿度検出手段を設けることなく求めることができる環境状態測定方法及び環境状態測定装置を得ることが目的である。
【００１２】
【課題を解決するための手段】
上記目的を達成するために請求項１記載の発明に係る環境状態測定方法は、被測定空間の周縁部の互いに異なる複数箇所の間の音波の伝搬時間を各々測定すると共に、前記被測定空間内でかつ前記被測定空間を区画する壁体から所定距離離れた位置の温度を測定し、前記複数箇所の間の音波の伝搬時間、及び前記壁体から所定距離離れた位置の温度に基づいて、前記被測定空間内の温度分布を演算により求める。
【００１３】
請求項１記載の発明では、被測定空間の周縁部の互いに異なる複数箇所の間の音波の伝搬時間を各々測定すると共に、被測定空間内でかつ被測定空間を区画する壁体から所定距離離れた位置の温度を測定している。なお、複数箇所の間の音波の伝搬時間は、例えば複数箇所の各々にスピーカ等の音波発生手段と、マイクロフォン等の音波検知手段を設け、或る箇所で音波を発生させてから他の箇所で音波が検知される迄の時間（音波の伝搬時間）を測定することを、前記複数箇所の間で各々行うことで得ることができる。
【００１４】
また、壁体から所定距離離れた位置の温度は、壁体から所定距離離れた位置に配設した温度センサ等の温度測定手段によって測定することができる。この所定距離としては、被測定空間を区画する壁体の温度の影響が十分に小さくなるか、又は前記影響が殆ど無視できる大きさとなる距離を適用することができ、例えば壁体からの距離と温度との関係を測定して定めることができる。これにより、壁体の温度の影響が非常に小さく、又は排除された被測定空間内の温度を得ることができる。
【００１５】
また、請求項１の発明では、上記のようにして測定した複数箇所の間の音波の伝搬時間、及び壁体から所定距離離れた位置の温度に基づいて、被測定空間内の温度分布を演算によって求めている。これは、例えば被測定空間を仮想的に多数の領域（ボクセル）に区切り、或る一対の箇所の間を音波が伝搬した際に前記一対の箇所の間に存在する各領域を音波が横切る距離を求め、前記一対の箇所の間の音波の伝搬時間と、前記各領域を音波が横切る際の伝搬速度との関係を、前記各領域を音波が横切る距離を用いて数式で表すことを前記複数箇所の間について各々行って連立方程式を立て、測定した複数箇所の間の音波の伝搬時間を連立方程式に代入すると共に、前記壁体から所定距離離れた位置に対応する領域については、測定した温度に基づき音波の伝搬速度を演算して連立方程式に代入した後に、伝搬速度が未知の領域について音波が横切る際の伝搬速度を演算し、伝搬速度から各領域の温度を演算することによって実現できる。
【００１６】
このように、請求項１の発明では、壁体から所定距離離れた位置の温度を測定し、複数箇所の間の音波の伝搬時間に加えて、壁体から所定距離離れた位置の温度の実測値を用いて被測定空間内の温度分布を演算しており、被測定空間内の温度分布の演算における未知数が減少すると共に、壁体の温度の影響を非常に小さく又は排除することができるので、被測定空間内の温度分布を精度良く求めることができる。
【００１７】
請求項２記載の発明に係る環境状態測定方法は、被測定空間を区画する壁体と、前記被測定空間内でかつ前記壁体から所定距離離れた位置と、の温度差を予め測定しておき、前記被測定空間の周縁部の互いに異なる複数箇所の間の音波の伝搬時間を各々測定すると共に、前記壁体の温度を測定し、測定した壁体の温度及び前記予め測定した温度差に基づいて、前記壁体から所定距離離れた位置の温度を推定し、前記複数箇所の間の音波の伝搬時間、及び前記壁体から所定距離離れた位置の温度に基づいて、前記被測定空間内の温度分布を演算により求める。
【００１８】
請求項２記載の発明では、被測定空間を区画する壁体と、被測定空間内でかつ壁体から所定距離離れた位置と、の温度差を予め測定しておき、壁体の温度を測定し、測定した壁体の温度及び予め測定した温度差に基づいて、壁体から所定距離離れた位置の温度を推定している。これにより、壁体の温度の影響が非常に小さく、又は排除された被測定空間内の温度を得ることができる。
【００１９】
そして請求項２の発明では、被測定空間の周縁部の互いに異なる複数箇所の間の音波の伝搬時間を各々測定し、複数箇所の間の音波の伝搬時間、及び壁体から所定距離離れた位置の温度に基づいて、被測定空間内の温度分布を演算により求めており、請求項１の発明と同様に、被測定空間内の温度分布の演算における未知数が減少すると共に、壁体の温度の影響を小さく又は排除することができるので、被測定空間内の温度分布を精度良く求めることができる。
【００２０】
また、請求項２の発明では、壁体から所定距離離れた位置の温度として推定値を用いているので、請求項１の発明と比較すると温度分布の演算精度は若干低下するものの、壁体から所定距離離れた位置に温度検出手段を設ける必要がないので、温度測定に際しての制約が少なくなる。
【００２１】
請求項３記載の発明に係る環境状態測定方法は、被測定空間の周縁部の互いに異なる複数箇所の間の音波の伝搬時間を各々測定し、前記複数箇所の間の音波の伝搬時間に基づいて、前記複数箇所の間の音波の伝搬経路を直線と仮定して前記被測定空間内の温度分布を演算により求め、求めた温度分布に基づき前記複数箇所の間の音波の伝搬経路を補正し、前記複数箇所の間の音波の伝搬時間、及び前記補正によって得られた複数箇所の間の音波の伝搬経路に基づいて、前記被測定空間内の温度分布を演算により再度求める。
【００２２】
被測定空間内における音波の伝搬経路は、被測定空間内の温度分布に依存して変化するが、請求項３記載の発明では、被測定空間の周縁部の互いに異なる複数箇所の間の音波の伝搬時間を各々測定し、複数箇所の間の音波の伝搬経路を直線と仮定して被測定空間内の温度分布を演算により求め、求めた温度分布に基づき複数箇所の間の音波の伝搬経路を補正している。これにより、実際の音波の伝搬経路により近い伝搬経路を得ることができる。
【００２３】
そして、請求項３の発明では、測定した複数箇所の間の音波の伝搬時間及び補正によって得られた複数箇所の間の音波の伝搬経路に基づいて、被測定空間内の温度分布を演算により再度求めているので、被測定空間内の温度分布を精度良く求めることができる。なお請求項３の発明において、音波の伝搬経路を補正し、被測定空間内の温度分布を再度求めることを、複数回繰り返すようにすれば、被測定空間内の温度を更に精度良く求めることができる。
【００２４】
ところで、本願発明者等は被測定空間内の湿度分布を測定するために、湿度に依存して変化する物理量について検討を行った。その結果、空気中を伝搬する音波の減衰率が、空気の温度、湿度、及び音波の周波数によって変化することを見い出し、被測定空間内の音波の減衰率の分布、被測定空間内の温度分布、及び音波の周波数から、被測定空間内の湿度分布を求めることができる、との知見を得た。
【００２５】
このため、請求項４記載の発明に係る環境状態測定方法は、被測定空間の周縁部の互いに異なる複数箇所の間の音波の伝搬時間及び減衰量を各々測定し、前記複数箇所の間の音波の伝搬時間に基づいて、前記被測定空間内の温度分布を演算により求めると共に、前記複数箇所の間の音波の減衰量に基づいて、前記被測定空間内の音波の減衰率の分布を演算によって求め、前記被測定空間内の音波の減衰率の分布、前記被測定空間内の温度分布、及び前記減衰量の測定に用いた音波の周波数に基づいて、前記被測定空間内の湿度分布を演算によって求める。
【００２６】
請求項４記載の発明では、被測定空間の周縁部の互いに異なる複数箇所の間の音波の伝搬時間及び減衰量を各々測定している。なお、複数箇所の間の音波の減衰量は、音波の伝搬時間の測定と同様に、例えば複数箇所の各々にスピーカ等の音波発生手段と、マイクロフォン等の音波検知手段を設け、或る箇所で発生させた所定の音量の音波が、他の箇所でどの程度の音量で検知されるかを測定し、音波発生時の音量、音波検知時の音量、音波発生手段及び音波検知手段の効率に基づいて音波の減衰量を求めることを、前記複数箇所の間で各々行うことによって得ることができる。
【００２７】
また請求項４の発明では、複数箇所の間の音波の伝搬時間に基づいて被測定空間内の温度分布を演算により求めると共に、複数箇所の間の音波の減衰量に基づいて被測定空間内の音波の減衰率の分布を演算によって求めている。被測定空間内の音波の減衰率の分布は、先に述べた被測定空間内の温度分布と同様に、例えば被測定空間を仮想的に多数の領域（ボクセル）に区切り、或る一対の箇所の間を音波が伝搬した際に前記一対の箇所の間に存在する各領域を音波が横切る距離を求め、前記一対の箇所の間を音波が伝搬した際の減衰量と、前記各領域における音波の減衰率との関係を、前記各領域を音波が横切る距離を用いて数式で表すことを前記複数箇所の間について各々行って連立方程式を立て、測定した複数箇所の間の音波の減衰量を連立方程式に代入した後に、各領域における音波の減衰率を演算することによって求めることができる。
【００２８】
そして請求項４の発明では、被測定空間内の音波の減衰率の分布、被測定空間内の温度分布、及び減衰量の測定に用いた音波の周波数に基づいて、被測定空間内の湿度分布を演算によって求めている。従って請求項４の発明によれば、従来は測定方法が確立されていなかった被測定空間内の湿度分布を求めることを、被測定空間の中央部付近に湿度センサ等の湿度検出手段を設けることなく実現することができる。
【００２９】
なお、請求項３又は請求項４の発明における被測定空間内の温度分布の演算に際しては、請求項５に記載したように、被測定空間内でかつ被測定空間を区画する壁体から所定距離離れた位置の温度を求め、被測定空間内の温度分布を、前記壁体から所定距離離れた位置の温度も用いて演算により求めることが好ましい。壁体から所定距離離れた位置の温度は、請求項１のように測定によって得るようにしてもよいし、請求項２のように壁体と壁体から所定距離離れた位置の温度差を予め測定しておき、壁体の温度を測定し、測定した壁体の温度及び予め測定した温度差に基づいて推定するようにしてもよい。これにより、被測定空間内の温度分布をより高精度に求めることができる。また請求項４の発明では、高精度に求めた温度分布に基づいて、被測定空間内の湿度分布もより高精度に求めることができる。
【００３０】
請求項６記載の発明は、請求項４の発明において、被測定空間内でかつ前記被測定空間を区画する壁体から所定距離離れた位置の湿度を求め、前記被測定空間内の湿度分布を、前記壁体から所定距離離れた位置の湿度も用いて演算により求めることを特徴としている。
【００３１】
前述のように、被測定空間内のうち壁体の温度は被測定空間内の温度と必ずしも一致しておらず、被測定空間内の壁体近傍では、壁体の温度の影響を受けて温度が大きく変化していることが一般的であるが、これに伴って、前記壁体近傍では湿度も大きく変化していることが多い。これに対し請求項６の発明では、被測定空間内でかつ被測定空間を区画する壁体から所定距離離れた位置の湿度を求め、被測定空間内の湿度分布を、壁体から所定距離離れた位置の湿度も用いて演算により求めており、被測定空間内の湿度分布の演算における未知数が減少すると共に、壁体の温度による影響を小さく又は排除することができるので、被測定空間内の湿度分布をより精度良く求めることができる。
【００３２】
なお、壁体から所定距離離れた位置の湿度は、請求項１と同様に測定によって得るようにしてもよいし、請求項２と同様に、壁体と壁体から所定距離離れた位置の湿度差を予め測定しておき、壁体の湿度を測定し、測定した壁体の湿度及び予め測定した湿度差に基づいて推定するようにしてもよい。
【００３３】
また本願発明者等は、空気中を伝搬する音波の減衰率と、空気の温度、湿度、及び音波の周波数との関係について考察し、空気の温度及び湿度が或る範囲内の場合には、湿度の変化に拘らず音波の減衰率が変化しない不感領域があること、及びこの不感領域が生ずるときの空気の温度及び湿度は、音波の周波数によって異なっている、との知見を得た。
【００３４】
上記に基づき請求項７記載の発明は、請求項４の発明において、前記複数箇所の間の音波の減衰量を測定し、前記被測定空間内の音波の減衰率の分布を演算によって求めることを、複数種の周波数の音波について各々行い、前記複数種の周波数の音波について各々求めた前記被測定空間内の音波の減衰率の分布、前記被測定空間内の温度分布、及び前記複数種の周波数に基づいて、前記被測定空間内の湿度分布を演算によって求めることを特徴としている。
【００３５】
請求項７記載の発明では、前記複数箇所の間の音波の減衰量を測定し、前記被測定空間内の音波の減衰率の分布を演算によって求めることを、複数種の周波数の音波について各々行い、複数種の周波数の音波について各々求めた前記被測定空間内の音波の減衰率の分布を用いて被測定空間内の湿度分布を演算によって求めている。これにより、被測定空間内の環境条件（温度及び湿度）が、特定の周波数の音波では湿度の変化に拘らず減衰率が変化しない条件（前記特定の周波数の音波では湿度の変化に対して減衰率が不感領域となる条件）であったとしても、別の周波数の音波の減衰率の分布から被測定空間内の湿度分布を求めることができる。
【００３６】
ところで、音波の減衰量の測定に際し、音波の減衰量が非常に大きい場合には、発生させる音波の音量を大きくしたり、高感度の音波検知手段によって音波を検知する必要が有り、音波の減衰量が非常に小さい場合には、高感度の音波検知手段によって音波の音量を精密に検知する必要が有る。そして、或る二点間を伝搬する音波の減衰量は、前記二点間の距離によって大きく変化する。
【００３７】
このため、請求項８記載の発明は、請求項４の発明において、前記複数箇所の間の音波の減衰量の測定に用いる音波の周波数を、被測定空間の大きさに応じて変更することを特徴としている。
【００３８】
音波の減衰率は音波の周波数によって大きく変化し、温度及び湿度が一定であっても音波の周波数が高周波になるに従って減衰率は高くなる。請求項８の発明では、音波の減衰量の測定に用いる音波の周波数を、被測定空間の大きさに応じて変更しているので、被測定空間の大きさ、すなわち複数箇所の間の距離に拘らず、複数箇所の間を伝搬する音波の減衰量が、測定に好適な値となるように調整することができる。これにより、複数箇所の間の音波の減衰量を効率良く測定することができる。
【００３９】
なお、請求項８の発明における音波の周波数の変更は、具体的には、被測定空間の大きさが大きくなるに従って、すなわち複数箇所の間の距離が大きくなるに従って、周波数が低くなるように変更することができる。
【００４０】
請求項９記載の発明に係る環境状態測定装置は、被測定空間の周縁部の互いに異なる複数箇所の間の音波の伝搬時間を各々測定する伝搬時間測定手段と、前記被測定空間内でかつ前記被測定空間を区画する壁体から所定距離離れた位置の温度を測定する温度測定手段と、前記伝搬時間測定手段によって測定された前記複数箇所の間の音波の伝搬時間、及び前記温度測定手段によって測定された前記壁体から所定距離離れた位置の温度に基づいて、前記被測定空間内の温度分布を演算により求める温度分布演算手段と、を含んで構成しているので、請求項１の発明と同様に、被測定空間内の温度分布を高精度に求めることができる。
【００４１】
請求項１０記載の発明に係る環境状態測定装置は、予め測定された、被測定空間を区画する壁体と、前記被測定空間内でかつ前記壁体から所定距離離れた位置と、の温度差を記憶する記憶手段と、前記被測定空間の周縁部の互いに異なる複数箇所の間の音波の伝搬時間を各々測定する伝搬時間測定手段と、前記壁体の温度を測定する壁体温度測定手段と、前記壁体温度測定手段によって測定された壁体の温度及び前記記憶手段に記憶されている温度差に基づいて、前記壁体から所定距離離れた位置の温度を推定する温度推定手段と、前記伝搬時間測定手段によって測定された前記複数箇所の間の音波の伝搬時間、及び前記温度推定手段によって推定された前記壁体から所定距離離れた位置の温度に基づいて、前記被測定空間内の温度分布を演算により求める温度分布演算手段と、を含んで構成しているので、請求項２の発明と同様に、被測定空間内の温度分布を高精度に求めることができる。
【００４２】
請求項１１記載の発明に係る環境状態測定装置は、被測定空間の周縁部の互いに異なる複数箇所の間の音波の伝搬時間を各々測定する伝搬時間測定手段と、前記伝搬時間測定手段によって測定された前記複数箇所の間の音波の伝搬時間に基づいて、前記複数箇所の間の音波の伝搬経路を直線と仮定して前記被測定空間内の温度分布を演算により求め、求めた温度分布に基づき前記複数箇所の間の音波の伝搬経路を補正し、前記複数箇所の間の音波の伝搬時間、及び前記補正によって得られた複数箇所の間の音波の伝搬経路に基づいて、前記被測定空間内の温度分布を演算により再度求める温度分布演算手段と、を含んで構成しているので、請求項３の発明と同様に、被測定空間内の温度分布を高精度に求めることができる。
【００４３】
請求項１２記載の発明に係る環境状態測定装置は、被測定空間の周縁部の互いに異なる複数箇所の間の音波の伝搬時間及び減衰量を各々測定する伝搬時間・減衰量測定手段と、前記伝搬時間・減衰量測定手段によって測定された前記複数箇所の間の音波の伝搬時間に基づいて、前記被測定空間内の温度分布を演算により求める温度分布演算手段と、前記伝搬時間・減衰量測定手段によって測定された前記複数箇所の間の音波の減衰量に基づいて前記被測定空間内の音波の減衰率の分布を演算によって求め、前記被測定空間内の音波の減衰率の分布及び前記温度分布演算手段によって演算された前記被測定空間内の温度分布に基づいて、前記被測定空間内の湿度分布を演算によって求める湿度分布演算手段と、を含んで構成しているので、請求項４の発明と同様に、被測定空間の中央部付近に湿度センサ等の湿度検出手段を設けることなく、被測定空間内の湿度分布を求めることができる。
【００４４】
【発明の実施の形態】
以下、図面を参照して本発明の実施形態の一例を詳細に説明する。
【００４５】
〔第１実施形態〕
図１には、本第１実施形態に係る環境状態測定装置１０が示されている。環境状態測定装置１０は、多数のノード１２Ａ、１２Ｂ、１２Ｃ、…と、ホストコンピュータ１４を備えており、これらが伝送媒体１６を介して互いに接続されて構成されている。ホストコンピュータ１４は、ＣＰＵ１４Ａ、ＲＯＭ１４Ｂ、ＲＡＭ１４Ｃ、入出力ポート１４Ｄを備え、これらがバス１４Ｅを介して互いに接続されて構成されている。入出力ポート１４Ｄはネットワーク伝送部４９を介して伝送媒体１６に接続されている。また、入出力ポート１４Ｄには、各種の情報を表示するためのディスプレイ１８と、オペレータが各種のデータやコマンド等を入力するためのキーボード２０が接続されている。
【００４６】
一方、多数のノード１２は、被測定空間の周縁部に、周縁部の全周に亘って略一定の間隔で配置される。例として図２には、円筒状の壁体２２によって区画され、上方が屋根（図示省略）によって閉塞されたドーム状の被測定空間２４（例えばドーム球場等）内の環境状態を測定する場合の各ノード１２の配置が示されており、各ノード１２は壁体２２に埋設されている。なお、ノード１２の数及び間隔は図２に示した例に限定されるものではない。また、各ノード１２を識別するために、各ノードには便宜的にノード番号が付与されている。
【００４７】
ノード１２Ａ、１２Ｂ、１２Ｃ、…は各々同一の構成であるので、以下では図１を参照し、ノード１２Ａの構成について説明する。ノード１２Ａは、各々被測定空間２４内に向けて配置された、音波発生手段としてのスピーカ３０、及び音波検知手段としてのマイクロフォン３２を備えている。スピーカ３０は、増幅器３４、Ｄ／Ａ変換器３６を介してデータ処理部３８に接続されており、マイクロフォン３２は、増幅器４０、Ａ／Ｄ変換器４２を介してデータ処理部３８に接続されている。
【００４８】
データ処理部３８はマイクロコンピュータ等を含んで構成され、ネットワーク伝送部４８を介して伝送媒体１６に接続されており、伝送媒体１６及びネットワーク伝送部４８を介してホストコンピュータ１４から、周波数ｆの音波の発生が指示されると、指示された周波数ｆの音波を指示された音量で発生させるためのデータを、所定のタイミングでＤ／Ａ変換器３６に出力する。データ処理部３８から出力されたデータは、Ｄ／Ａ変換器３６でアナログの電気信号に変換され、増幅器３４で増幅された後にスピーカ３０に供給される。これにより、スピーカ３０からは、周波数ｆの音波が所定の時刻に指示された音量でスピーカ３０から発せられる。
【００４９】
また、マイクロフォン３２から出力される電気信号は、増幅器４０で増幅されＡ／Ｄ変換器４２でディジタルデータに変換されてデータ処理部３８に入力される。データ処理部３８は、Ａ／Ｄ変換器４２を介して入力されるデータに基づいて、マイクロフォン３２から出力される電気信号のレベルの変化を監視し、閾値以上の音量の音波が到来したか否か判断する。そして、閾値以上の音量の音波が到来したと判断した場合には、到来した音波の音量及び音波の到来時刻（検知時刻）をメモリ等に記憶する。そして、ネットワーク伝送部４８及び伝送媒体１６を介してホストコンピュータ１４から音波の検知時刻及び音量の問い合わせがあった場合には、記憶している音波の検知時刻及び音量を表すデータをホストコンピュータ１４に送信する。
【００５０】
また、ノード１２は熱電対等から成る温度センサ４４を備えている。被測定空間２４内の温度は壁体２２の温度と異なっており、壁体２２の近傍では、壁体２２の壁面に近づくに従って壁体２２の温度の影響をより強く受けるので、図３（Ｂ）に示すように、被測定空間２４内の温度は壁体２２の近傍で大きく変化している。このため、本実施形態では、壁体２２の壁面と、壁体２２の温度の影響が略０となる位置との距離Ｌ１を予め求め、図３（Ａ）に示すように、温度センサ４４が、壁体２２の壁面から被測定空間２４の中央部に向けて距離Ｌ１だけ隔てた位置の温度を測定するように構成されている。
【００５１】
温度センサ４４は請求項９に記載の温度測定手段に対応しており、温度センサ４４を駆動するドライバ４６を介してデータ処理部３８に接続されている。データ処理部３８は、伝送媒体１６及びネットワーク伝送部４８を介してホストコンピュータ１４から温度の測定が指示されると、壁体２２の壁面から距離Ｌだけ隔てた位置の温度を温度センサ４４によって測定し、測定結果を表す温度データをホストコンピュータ１４に送信する。
【００５２】
次に本第１実施形態の作用として、まず図４のフローチャートを参照し、被測定空間２４内の環境状態（温度分布及び湿度分布）を求める場合にホストコンピュータ１４で実行される伝搬時間・減衰量測定処理について説明する。
【００５３】
ステップ１００では、被測定空間２４の大きさに基づき、測定に用いる音波の周波数を複数種設定する。具体的には、複数種の周波数（ｆ１ 、ｆ２ 、…）として、スピーカ３０が音波を発生可能でマイクロフォン３２が音波を検知可能な周波数帯域内の値を設定すると共に、被測定空間２４の大きさが大きくなるに従って（すなわち被測定空間１２の中央部を挟んで対向するノード１２間の距離が大きくなるに従って）、前記複数種の周波数が全体的に低周波側に偏倚するように設定する。上記の周波数の設定は請求項８の発明に対応している。
【００５４】
次のステップ１０２では、ステップ１００で設定した複数種の周波数（ｆ１ 、ｆ２ 、…）の中から、測定未実行の周波数ｆを選択する。ステップ１０４ではカウンタｍに１を代入し、ステップ１０６では、ノード番号ｍのノード１２（以下、単にノードｍという）に対し、スピーカ３０から所定の時刻に周波数ｆの音波を所定の音量で発生させるよう指示する。これにより、ノードｍのスピーカ３０からは、所定の時刻に周波数ｆの音波が所定の音量で発生され、発生された音波は他のノード１２のマイクロフォン３２で各々検知され、音波の検知時刻及び音量が各々メモリ等に記憶される。
【００５５】
ところで本実施形態では、音波を発生させるノードｍに対し、残りのノードのうちノードｍから比較的離れた位置に存在している所定数のノード１２のみを測定対象のノードとしている（図２及び図６には、各ノードから音波を発生させるときの測定対象のノードを線で結んで示す（この線は音波の伝搬時間及び総減衰量の測定対象区間を表している））。このため、ステップ１０８では、ノードｍに対応する測定対象の所定数のノード１２のノード番号（ｎ１ 、ｎ２ 、…）を各々判断し、次のステップ１１０では、判断したノード番号に基づいて、測定対象の所定数のノード１２に対して音波の検知時刻及び音量を各々問い合わせ、測定対象の所定数のノード１２から送信された音波の検知時刻及び音量を表すデータを各々取り込む。
【００５６】
ステップ１１２では、ノードｍから音波が発生された時刻及び発生された音波の音量、測定対象の所定数のノード１２の各々における音波の検知時刻及び音量に基づいて、ノードｍのスピーカ３０から測定対象のノードｎのマイクロフォン３２への音波（周波数ｆの音波）の伝搬時間ｔｍｎ〔秒〕、ノードｍのスピーカ３０から測定対象のノードｎのマイクロフォン３２へ周波数ｆの音波が伝搬する際の総減衰量Ｒｍｎｆ 〔ｄＢ〕を、測定対象の所定数のノード１２について各々演算し、周波数ｆと対応させて記憶する。
【００５７】
次のステップ１１４では、カウンタｍの値が最終ノード番号に一致したか否か判定する。判定が否定された場合にはステップ１１６へ移行し、カウンタｍの値を１だけインクリメントしてステップ１０６に戻る。これにより、ステップ１１４の判定が肯定される迄の間は、音波を発するノード１２を順に切り換えてステップ１０６〜１１２の処理が繰り返され、例として図２に示す全ての測定対象区間について、音波の検知時刻及び音量の測定、音波の伝搬時間ｔｍｎ及び総減衰量Ｒｍｎｆ の演算が行われることになる。
【００５８】
上記のステップ１０６〜１１６は、各ノードのスピーカ３０、マイクロフォン３２、データ処理部３８と共に、請求項９乃至請求項１１に記載の伝搬時間測定手段、請求項１２に記載の伝搬時間・減衰量測定手段に対応している。
【００５９】
なお、上記では全てのノード１２から各々音波を発生させ、測定対象の所定数のノードから音波の検知時刻及び音量を毎回取り込んで音波の伝搬時間ｔｍｎ及び総減衰量Ｒｍｎｆ を演算しており、例として図２に示す全ての測定対象区間の各々に対し、２方向（一方のノードから他方のノードへ向かう方向、及び他方のノードから一方のノードに向かう方向）について音波の検知時刻及び音量の測定、音波の伝搬時間ｔｍｎ及び総減衰量Ｒｍｎｆ の演算を行っているが、同一の区間の前記２方向についての音波の伝搬時間及び総減衰量は一般に等しいので、各区間に対し何れか一方の方向についてのみ、音波の検知時刻及び音量の測定、音波の伝搬時間ｔｍｎ及び総減衰量Ｒｍｎｆ の演算を行うようにしてもよい。
【００６０】
これにより、図２の例では１個のノードにはマイクロフォン３２を設ける必要がなくなり、別の５個のノードにはスピーカ３０を設ける必要がなくなるので、構成を簡単にすることができると共に、図４に示した伝搬時間・減衰量測定処理に要する時間も短縮することができる。
【００６１】
ステップ１１４の判定が肯定されるとステップ１１８に移行し、先のステップ１００で設定した複数種の周波数の各々について、全ての測定対象区間に対する音波の検知時刻及び音量の測定、及び音波の伝搬時間ｔｍｎ及び総減衰量Ｒｍｎｆ の演算を行ったか否か判定する。判定が否定された場合にはステップ１０２に戻り、複数種の周波数の中から測定未実行の別の周波数ｆを選択し、各ノードから選択した周波数ｆの音波を発生させて、全ての測定対象区間に対する音波の検知時刻及び音量の測定、及び音波の伝搬時間ｔｍｎ及び総減衰量Ｒｍｎｆ の演算を行う。これにより、全ての測定対象区間における音波の伝搬時間及び総減衰量が、ステップ１００で選択した複数種の周波数の音波について各々求められることになる。
【００６２】
なお、音速は温度にのみ依存し、同一の区間を伝搬する音波の伝搬時間は、音波の周波数に拘らず等しいので、音波の伝搬時間は単一の周波数の音波についてのみ測定・演算を行い、音波の総減衰量は複数種の周波数について測定・演算を行うようにしてもよい。
【００６３】
全ての測定対象区間における音波の伝搬時間及び総減衰量が、複数種の周波数について各々求められると、ステップ１１８の判定が肯定されてステップ１２０へ移行する。ステップ１２０では、全てのノード１２に対し温度センサ４４による温度の測定を指示する。そして次のステップ１２２では、全てのノード１２からの温度データの取込みを行う。
【００６４】
上記のようにして伝搬時間・減衰量測定処理が行われると、続いてホストコンピュータ１４では温度・湿度分布演算処理が実行される。この温度・湿度分布演算処理について、図５のフローチャートを参照して説明する。
【００６５】
ステップ１５０では、複数種の周波数について各々測定・演算した音波の伝搬時間のうち、特定の周波数について測定・演算した全ての測定対象区間における音波の伝搬時間ｔｍｎを取り込む。
【００６６】
ところで、本実施形態では図６に示すように、被測定空間２４内を仮想的に多数の矩形状の領域（ボクセル）に区切り、各ボクセル毎に温度及び湿度を演算する。また本実施形態では、各測定対象区間における音波の伝搬経路を各々直線と仮定したときの、音波の伝搬経路が該伝搬経路上に存在する各ボクセルを横切っている距離ｌｍｎｘｙ〔ｍ〕（図６参照、但しｍは測定対象区間の一端の音波発生側のノードのノード番号、ｎは測定対象区間の他端の音波検知側のノードのノード番号、ｘ及びｙは各ボクセルを識別する符号である）が、各測定対象区間について予め演算されてＲＯＭ１４Ｂ等に記憶されている。
【００６７】
次のステップ１５２では、音波の伝搬経路が該伝搬経路上に存在する各ボクセルを横切っている距離ｌｍｎｘｙをＲＯＭ１４Ｂから取り込む。なお、距離ｌｍｎｘｙはボクセルの大きさ及び位置、ノードｍ及びノードｎの位置から幾何学的に算出可能である。
【００６８】
次のステップ１５４では、ステップ１５２で取り込んだ距離ｌｍｎｘｙを用い、全ての測定対象区間について、音波の伝搬時間ｔｍｎと、音波の伝搬経路上に存在する各ボクセルを音波が通過する際の音速ｖｘｙ〔ｍ／秒〕との関係を表す方程式を各々生成することにより、被測定空間２４内の各ボクセルの温度Ｔｘｙを各々求めるための連立方程式を生成する。以下、この連立方程式について説明する。
【００６９】
音速ｖは次の（１）式に示すように、音波が伝搬する媒体（空気）の温度Ｔにのみ依存することが知られている。
【００７０】
ｖ＝ ３３１．４５ ＋ ０．６０７・Ｔ 〔ｍ／秒〕 …（１）
ノードｍのスピーカ３０からノードｎのマイクロフォン３２に至る測定対象区間を音波が伝搬する際に、音波の伝搬経路上に存在するボクセルｘｙを音波が通過する時間ｔｍｎｘｙ〔秒〕は次の（２）式により表され、この（２）式を（１）式に代入すると次の（３）式が得られる。
【００７１】
ｔｍｎｘｙ＝ｌｍｎｘｙ／ｖｘｙ …（２）
ｔｍｎｘｙ＝ｌｍｎｘｙ／（３３１．４５＋ ０．６０７・Ｔｘｙ） …（３）
但し、Ｔｘｙはボクセルｘｙ内の温度〔℃〕である。ノードｍのスピーカ３０からノードｎのマイクロフォン３２に至る音波の伝搬経路は多数のボクセルを横切っており、スピーカｍから前記経路を伝搬してマイクロフォンｎに至る音波の伝搬時間ｔｍｎは、前記多数のボクセルの各々における音波の通過時間の総和であるので、
上記の（４）式で表すことができる。ここで、
Ｑｘｙ＝１／ｖｘｙ＝１／（３３１．４５＋ ０．６０７・Ｔｘｙ） …（５）
とすれば、先の（４）式から、ノードｍのスピーカ３０からノードｎのマイクロフォン３２に至る音波の伝搬経路が横切っている各ボクセル内の温度の関数であるＱｘｙ（各ボクセルを音波が通過する際の音速ｖｘｙの逆数）の一次多項式である（６）式が得られる。
【００７２】
ｔｍｎ＝Σ（ｌｍｎｘｙ・Ｑｘｙ） …（６）
各測定対象区間における音波の伝搬時間ｔｍｎ及び距離ｌｍｎｘｙは既知であるので、ステップ１５４では、各測定対象区間における音波の伝搬時間ｔｍｎ及び距離ｌｍｎｘｙを（６）式に各々代入することにより、測定対象区間の数と同数の１次多項式（連立方程式）を生成する。
【００７３】
なお、上記の連立方程式を解くことで、各ボクセルについて変数Ｑｘｙの値を得ることができ、変数Ｑｘｙを（５）式に代入することで各ボクセル内の温度Ｔｘｙを逆算により求めることができるが、被測定空間２４内の周縁部付近に位置している一部のボクセルについては、温度センサ４４によって温度が実測されている。このためステップ１５６では、温度が既知のボクセルについて、各ノードから取込んだ温度データが表す温度Ｔｘｙを（５）式に代入して変数Ｑｘｙの値を演算し、演算したＱｘｙの値を連立方程式に代入する。
【００７４】
そして次のステップ１５８では、温度が既知のボクセルについて変数Ｑｘｙの値を代入した連立方程式を解いて、温度が未知のボクセルの変数Ｑｘｙの値を求め、求めた変数Ｑｘｙの値を（５）式に代入することで各ボクセル内の温度Ｔｘｙを逆算によって求める。これにより、被測定空間２４内の各ボクセル毎の温度、すなわち被測定空間２４内の温度分布が演算により求まることになる。
【００７５】
上記では温度センサ４４によって壁体２２の壁面から距離Ｌ１だけ隔てた位置（壁体２２の温度の影響が略０となる位置）の温度を測定し、測定した温度を連立方程式に代入した後に、連立方程式を解いて各ボクセルの温度を求めており、被測定空間２４の周縁部の温度の実測値を用いて温度分布を求めていると共に、前記実測値は壁体２２の温度の影響が略０となる位置で測定した値であるので、被測定空間２４内の温度分布を精度良く求めることができる。
【００７６】
ステップ１６０では、先のステップ１５８で求めた被測定空間２４内の温度分布に基づき、各測定対象区間における音波の伝搬経路を各々補正する。この音波の伝搬経路の補正は、例えば以下のように行うことができる。
【００７７】
例として図７（Ａ）に示すように被測定空間２４内の温度分布が求まり、測定対象区間の１つであるノードｍとノードｎとの間における音波の伝搬経路を補正する場合、補正前の伝搬経路（図７（Ａ）では直線で示す）に直交する方向に沿った温度変化の勾配を、補正前の伝搬経路の一端から他端に亘って演算する。なお、図７（Ｂ）には、伝搬経路に直交する方向に沿った温度変化の勾配を、勾配が大きくなるに従って長さを長くした太線で示している。
【００７８】
次に、伝搬経路上の各部分における曲率半径を、温度変化の勾配の大きさに反比例するように（温度変化の勾配が大きくなるに従って曲率半径が小さくなるように）決定し、伝搬経路に直交する方向に沿った温度変化における高温側が凸となるように、決定した曲率半径に従って伝搬経路を湾曲させる。図７（Ｃ）に示しているノードｍとノードｎを結ぶ曲線は、決定した曲率半径に従って湾曲させた後の伝搬経路を示している。
【００７９】
続いて、湾曲させた伝搬経路に対し、その両側に、伝搬経路から等しい間隔を隔てて伝搬経路に平行な一対の仮想線（図７（Ｃ）に示す仮想線５０Ａ、５０Ｂ参照）を設定すると共に、伝搬経路及び一対の仮想線を、伝搬経路に沿って所定間隔毎に、伝搬経路に直交する仮想的な区分線（図７（Ｃ）に示す区分線５２参照）によって複数の区間に区切り、各区間において、一対の仮想線上を伝搬する音波の伝搬時間が互いに等しくなるように、伝搬経路の湾曲度合いを更に補正する。
【００８０】
上記の補正処理を、全ての測定対象区間における音波の伝搬経路に対して各々行うと、次のステップ１６２へ移行し、ステップ１６０における伝搬経路の補正において、伝搬経路に対する補正量（例えば補正前の伝搬経路に対する補正後の伝搬経路の距離の最大値や、補正前の伝搬経路と補正後の伝搬経路とで囲まれた部分の面積等）が所定値以上となった伝搬経路が有るか否か判定する。判定が肯定された場合には、ステップ１６４へ移行する。
【００８１】
図８（Ａ）と（Ｂ）を比較しても明らかなように、音波の伝搬経路を補正した場合、伝搬経路上に存在するボクセルが変化したり、伝搬経路が各ボクセルを横切っている距離が変化する。このためステップ１６４では、先のステップ１６０の補正を行った後の全ての測定対象区間における音波の伝搬経路が、該伝搬経路上に存在する各ボクセルを横切っている距離ｌｍｎｘｙを各々演算し、ステップ１５４に戻る。
【００８２】
これにより、ステップ１５４以降では、補正後の音波の伝搬経路に基づいて被測定空間２４内の温度分布が再度演算されるので、例として図８（Ｃ）にも示すように一部のボクセルの温度Ｔｘｙが変化し、各測定対象区間における音波の伝搬経路を直線と仮定して被測定空間２４内の温度分布を演算したことによって生じた誤差が小さくされる。そして、ステップ１５４〜１６４はステップ１６２の判定が肯定されている間繰り返されるので、被測定空間２４内の温度分布を高精度に求めることができる。
【００８３】
なお、上述したステップ１５０〜１６４は、請求項９乃至請求項１２に記載の温度分布演算手段に対応しており、より詳しくは、ステップ１５２、及びステップ１５４〜１５８の処理を第１回目に実行することは請求項１１に記載の温度分布演算手段における「複数箇所の間の音波の伝搬経路を直線と仮定して被測定空間内の温度分布を演算により求め」ることに対応しており、ステップ１６０は、請求項１１に記載の温度分布演算手段における「温度分布に基づき複数箇所の間の音波の伝搬経路を補正」することに対応しており、ステップ１６２の判定が肯定されることにより、ステップ１６４、及びステップ１５４以降の処理を再度実行することは、請求項１１に記載の温度分布演算手段における「複数箇所の間の音波の伝搬時間、及び補正によって得られた複数箇所の間の音波の伝搬経路に基づいて、被測定空間内の温度分布を演算により再度求める」ことに対応している。
【００８４】
ステップ１６２の判定が否定されるとステップ１６６へ移行し、複数種の周波数について各々測定・演算した音波の総減衰量のうち、或る周波数ｆについて測定・演算した全ての測定対象区間における音波の総減衰量Ｒｍｎｆ を取り込む。次のステップ１６８では、ステップ１６６で取り込んだ総減衰量Ｒｍｎｆ を用い、全ての測定対象区間について、音波の総減衰量Ｒｍｎｆ と、音波の伝搬経路上に存在する各ボクセルの周波数ｆの音波についての減衰率Ｒｘｙｆ との関係を表す方程式を各々生成し、被測定空間２４内の各ボクセルの周波数ｆの音波についての減衰率Ｒｘｙｆ を各々求めるための連立方程式を生成する。以下、この連立方程式について説明する。
【００８５】
ノードｍのスピーカ３０からノードｎのマイクロフォン３２に至る測定対象区間を周波数ｆの音波が伝搬する際に、音波の伝搬経路上に存在するボクセルｘｙにおける減衰量Ｒｍｎｘｙｆ は、ボクセルｘｙを周波数ｆの音波が通過する際の単位距離当りの減衰量（減衰率）をＲｘｙｆ 〔ｄＢ／ｍ〕とすると、
Ｒｍｎｘｙｆ ＝ｌｍｎｘｙ・Ｒｘｙｆ …（７）
上記の（７）式で表すことができる。従って、ノードｍのスピーカ３０からノードｎのマイクロフォン３２に至る測定対象区間を周波数ｆの音波が伝搬する際の総減衰量Ｒｍｎｆ 〔ｄＢ〕は、各ボクセルにおける周波数ｆの音波についての減衰率Ｒｘｙｆ の一次多項式である次の（８）式で表すことができる。
【００８６】
但し、Φｍｎｓｐはノードｍのスピーカ３０の指向特性や周波数特性等を考慮した電気−音響変換特性、Φｍｎｍｉｃ はノードｎのマイクロフォン３２の指向特性や周波数特性等を考慮した音響−電気変換特性である。
【００８７】
各測定対象区間における周波数ｆの音波の総減衰量Ｒｍｎｆ 及び距離ｌｍｎｘｙは既知であり、本実施形態では、変換特性Φｍｎｓｐ及びΦｍｎｍｉｃ は、スピーカ３０やマイクロフォン３２の指向特性や周波数特性に基づき予め算出されてＲＯＭ１４Ｂに記憶されている。従って、ステップ１６８では、各測定対象区間における音波の総減衰量Ｒｍｎｆ 、及び補正後の音波の伝搬経路から先のステップ１６４で演算した距離ｌｍｎｘｙを（８）式に各々代入することにより、測定対象区間の数と同数の１次多項式（連立方程式）を生成する。
【００８８】
次のステップ１７０では、上記の連立方程式を解いて各ボクセルにおける周波数ｆの音波についての減衰率Ｒｘｙｆ を求める。これにより、周波数ｆの音波についての被測定空間２４内の減衰率の分布が演算により求まることになる。次のステップ１７２では複数種の周波数について被測定空間２４内の減衰率の分布を求めたか否か判定する。判定が否定された場合にはステップ１７２に戻り、減衰率の分布を求めていない他の周波数ｆについて測定・演算した音波の総減衰量Ｒｍｎｆ を取り込み、上記と同様にして被測定空間２４内の減衰率の分布を求める。
【００８９】
ステップ１６６〜１７２を複数回繰り返し、複数種の周波数について被測定空間２４内の減衰率の分布を全て求めると、ステップ１７２の判定が肯定されてステップ１７４へ移行する。ステップ１７４では、被測定空間２４内の温度分布、及び複数種の周波数について各々求めた被測定空間２４内の減衰率の分布に基づいて、各ボクセルの湿度を各々演算する。
【００９０】
具体的には、図９〜図１３から明らかなように、空気中を通過する音波の単位距離当りの減衰率（距離減衰率）は、気温及び湿度によって変化すると共に、音波の周波数によっても変化する。また、空気の温度及び音波の周波数によっては、湿度の変化に拘らず減衰率が変化しない不感領域が生ずる。例えば図１０に示す気温が０℃の場合には、周波数 ２５０〔ＨＺ 〕の音波について湿度の変化に拘らず減衰率が殆ど変化しない不感領域が生じ、同様に、図１１に示す気温が１０℃の場合には周波数 ２５０〔ＨＺ 〕及び ５００〔ＨＺ 〕の音波、図１２に示す気温が２０℃の場合には周波数２〔ｋＨＺ 〕の音波、図１３に示す気温が３０℃の場合には周波数４〔ｋＨＺ 〕の音波について湿度の変化に拘らず減衰率が殆ど変化しない不感領域が生じている。
【００９１】
本実施形態では、上記の不感領域の影響で湿度が不定となることを回避するために複数種の周波数について減衰率の分布を求めている。このため、湿度の演算にあたっては、温度、周波数及び減衰率をパラメータ（入力）とし、これらのパラメータに合致する湿度を出力する３次元のルックアップテーブルを作成するか、或いは次の（９）式のように、温度、周波数及び減衰率と湿度との関係を関数式として定めておく。
【００９２】
Ｈｘｙ＝ｆｕｎｃｔｉｏｎ（Ｔｘｙ，ｆ，Ｒｘｙｆ ） …（９）
但し、Ｈｘｙはボクセルｘｙ内の相対湿度〔％〕である。
【００９３】
そして、先に求めた被測定空間２４内の温度分布、先に求めた複数種の周波数についての被測定空間２４内の減衰率の分布から、前述のルックアップテーブル又は関数式を用いて、各ボクセル毎に相対湿度Ｈｘｙを求める。これにより、被測定空間内の湿度分布が求まることになる。
【００９４】
なお、不感領域に関しては、温度、減衰率、及び周波数の各値の組み合わせから、不感領域に相当する条件か否かを判断し、不感領域に相当する条件であると判断したデータを用いないようにするか、或いは入力が不感領域に相当する条件のときには相対湿度データが出力されないようにルックアップテーブルを作成しておけば、不感領域の影響を受けることなく、被測定空間２４内の湿度分布を精度良く求めることができる。
【００９５】
上述したステップ１６６〜１７４は請求項１２に記載の湿度分布演算手段に対応しており、より詳しくは、ステップ１６６〜１７２は、前記湿度分布演算手段における「複数箇所の間の音波の減衰量に基づいて被測定空間内の音波の減衰率の分布を演算によって求め」ることに対応しており、ステップ１７４は、前記湿度分布演算手段における「被測定空間内の音波の減衰率の分布、被測定空間内の温度分布、及び減衰量の測定に用いた音波の周波数に基づいて、被測定空間内の湿度分布を演算によって求める」ことに対応している。また、ステップ１７４において複数種の周波数についての被測定空間２４内の減衰率の分布を用いることは、請求項７記載の発明に対応している。
【００９６】
次のステップ１７６では、上記により得られた被測定空間２４内の温度分布及び湿度分布を、例として図１４に示すコンターマップ等の形態でディスプレイ１８に表示し、処理を終了する。なお、ディスプレイ１８への表示に代えて、被空調空間２４内の空気調和を行う空調装置に対し、被測定空間２４内の温度分布及び湿度分布を表すデータを出力するようにしてもよい。これにより、空調装置が被測定空間２４に対し、快適性と省エネルギーの双方を満足する空気調和を行うことが可能となる。
【００９７】
〔第２実施形態〕
次に本発明の第２実施形態について説明する。なお、第１実施形態と同一の部分には同一の符号を付し、説明を省略する。図１５には本第２実施形態に係る環境状態測定装置６０が示されている。
【００９８】
この環境状態測定装置６０では、各ノード１２が、壁体２２の壁面近傍の位置の湿度を検出する湿度センサ６２を備えている。湿度センサ６２はドライバ６４を介してデータ処理部３８に接続されている。データ処理部３８は、伝送媒体１６及びネットワーク伝送部４８を介してホストコンピュータ１４から湿度の測定が指示されると、壁体２２の壁面近傍の位置の湿度を湿度センサ６２によって測定し、測定結果を表す湿度データをホストコンピュータ１４に送信する。
【００９９】
また、本第２実施形態に係る温度センサ４４は、請求項１０に記載の壁体温度測定手段に対応しており、図１６（Ａ）に示すように、壁体２２の壁面近傍の温度（より詳しくは壁体２２の壁面から被測定空間２４の中央部に向けて距離Ｌ２（Ｌ２≪Ｌ１）だけ隔てた位置の温度）を測定するように構成されている。
【０１００】
本第２実施形態では、温度センサ４４によって測定される壁体２２の壁面から距離Ｌ２だけ隔てた位置の温度と、壁体２２の温度の影響が略０となる位置（例えば前記壁面から距離Ｌ１以上離れた位置）の温度と、の温度差ΔＴが予め測定され、ホストコンピュータ１４のＲＯＭ１４Ｂ（請求項１０に記載の記憶手段に相当）に記憶されている。また、湿度に関しても、湿度センサ６２によって測定される壁体２２の壁面近傍の湿度と、壁体２２の温度の影響が略０となる位置（例えば前記壁面から距離Ｌ１以上離れた位置）の湿度と、の温度差ΔＨが予め測定され、ホストコンピュータ１４のＲＯＭ１４Ｂに記憶されている。
【０１０１】
次に本第２実施形態の作用として、まず図１７のフローチャートを参照し、本第２実施形態に係る伝搬時間・減衰量測定処理について、第１実施形態に係る伝搬時間・減衰量測定処理（図４）と異なる部分についてのみ説明する。
【０１０２】
本第２実施形態に係る伝搬時間・減衰量測定処理では、ステップ１００〜１１８において、複数種の周波数について、全ての測定対象区間に対する音波の検知時刻及び音量の測定、及び音波の伝搬時間ｔｍｎ及び総減衰量Ｒｍｎｆ の演算を行った後に、ステップ１２１において、全てのノード１２に対し温度センサ４４による温度の測定及び湿度センサ６２による湿度の測定を指示する。そして次のステップ１２３では、全てのノード１２からの温度データ及び湿度データの取込みを行う。
【０１０３】
次に図１８のフローチャートを参照し、本第２実施形態に係る温度・湿度分布演算処理について、第１実施形態に係る温度・湿度分布演算処理（図５）と異なる部分についてのみ説明する。
【０１０４】
本第２実施形態に係る温度・湿度分布演算処理では、ステップ１５４で被測定空間２４内の各ボクセルの温度Ｔｘｙを各々求めるための連立方程式を生成した後に、次のステップ１５５において、各ノードから取込んだ温度データを、ＲＯＭ１４Ｂに記憶されている温度差ΔＴにより補正する。このステップ１５５は、請求項１０に記載の温度推定手段に対応している。これにより、壁体２２の温度の影響が略０となるように補正した温度データを得ることができる。そして、次のステップ１５７では、補正後の温度データが表す温度Ｔｘｙを（５）式に代入して変数Ｑｘｙの値を演算し、演算したＱｘｙの値を連立方程式に代入する。
【０１０５】
本第２実施形態では、壁体２２の壁面近傍の温度を測定し、測定によって得られた温度データを、予め測定した温度差ΔＴにより壁体２２の温度の影響が略０となるように補正した後に連立方程式に代入しているので、第１実施形態と同様に被測定空間２４内の温度分布を精度良く求めることができる。また、本第２実施形態によれば、温度センサ４４を壁体２２の壁面から大きく突出させる必要がないので、見栄えが向上する。
【０１０６】
また、本第２実施形態に係る温度・湿度分布演算処理では、ステップ１６６〜１７２で複数種の周波数について被測定空間２４内の減衰率の分布を求めた後に、ステップ１７３において、各ノードから取込んだ湿度データを、ＲＯＭ１４Ｂに記憶されている湿度差ΔＨにより補正する。これにより、壁体２２の温度の影響が略０となるように補正した湿度データを得ることができる。そして、次のステップ１７５では、被測定空間２４内の温度分布、及び複数種の周波数について各々求めた被測定空間２４内の減衰率の分布に基づいて、各ボクセルの湿度を各々演算するが、補正後の湿度データにより相対湿度が既知である壁体２２の壁面近傍のボクセルｘｙについては、補正後の湿度データを相対湿度Ｈｘｙとして設定する。これにより被測定空間２４内の湿度分布が求まる。
【０１０７】
なお、上記のように、補正後の湿度データも用いて被測定空間２４内の湿度分布を求めることは、請求項６の発明に対応している。
【０１０８】
このように、本第２実施形態では、壁体２２の壁面近傍の湿度を測定し、測定によって得られた湿度データを、予め測定した湿度差ΔＨにより壁体２２の温度の影響が略０となるように補正し、補正後の湿度データも用いて被測定空間２４内の湿度分布を求めているので、被測定空間２４内の湿度分布をより精度良く求めることができる。
【０１０９】
なお、第２実施形態では湿度センサ６２によって壁体２２の壁面近傍の湿度を測定し、予め測定した湿度差ΔＨにより壁体２２の温度の影響が略０となるように湿度データを補正していたが、これに限定されるものではなく、第１実施形態で説明した温度センサ４４と同様に、湿度センサ６２によって壁体２２の温度の影響が略０となる位置（例えば壁面から距離Ｌ１だけ隔てた位置）の湿度を測定し、測定した湿度を補正することなく用いて被測定空間２４内の湿度分布を求めるようにしてもよい。
【０１１０】
また、上記では本発明を、円筒状の壁体によって区画され、上方が屋根によって閉塞されたドーム球場等のドーム状の被測定空間２４の温度分布、湿度分布の測定に適用した場合を説明したが、被測定空間の形状や用途等は上記に限定されるものではなく、本発明は、例えば劇場、病院等の空間において快適性と省エネルギーの双方を満足する空気調和を行う場合や、或いはクリーンルーム内の温度及び湿度を精密に所望の値に制御する等の場合に適用可能であることは言うまでもない。
【０１１１】
【発明の効果】
以上説明したように請求項１及び請求項９記載の発明は、被測定空間の周縁部の互いに異なる複数箇所の間の音波の伝搬時間を各々測定すると共に、被測定空間を区画する壁体から所定距離離れた位置の温度を測定し、音波の伝搬時間及び壁体から所定距離離れた位置の温度に基づいて、被測定空間内の温度分布を演算により求めるようにしたので、被測定空間内の温度分布を高精度に求めることができる、という優れた効果を有する。
【０１１２】
請求項２及び請求項１０記載の発明は、被測定空間を区画する壁体と該壁体から所定距離離れた位置との温度差を予め測定し、被測定空間の周縁部の互いに異なる複数箇所の間の音波の伝搬時間を各々測定すると共に、壁体の温度を測定し、測定した壁体の温度及び前記温度差に基づいて壁体から所定距離離れた位置の温度を推定し、音波の伝搬時間及び壁体から所定距離離れた位置の温度に基づいて、被測定空間内の温度分布を演算により求めるようにしたので、被測定空間内の温度分布を高精度に求めることができる、という優れた効果を有する。
【０１１３】
請求項３及び請求項１１記載の発明は、被測定空間の周縁部の互いに異なる複数箇所の間の音波の伝搬時間を各々測定し、複数箇所の間の音波の伝搬時間に基づいて音波の伝搬経路を直線と仮定して被測定空間内の温度分布を演算により求め、求めた温度分布に基づき音波の伝搬経路を補正し、音波の伝搬時間及び補正によって得られた音波の伝搬経路に基づいて、被測定空間内の温度分布を演算により再度求めるようにしたので、被測定空間内の温度分布を高精度に求めることができる、という優れた効果を有する。
【０１１４】
請求項４及び請求項１２記載の発明は、被測定空間の周縁部の互いに異なる複数箇所の間の音波の伝搬時間及び減衰量を各々測定し、音波の伝搬時間に基づいて被測定空間内の温度分布を演算により求めると共に、音波の減衰量に基づいて被測定空間内の音波の減衰率の分布を演算によって求め、音波の減衰率の分布、温度分布、及び減衰量の測定に用いた音波の周波数に基づいて、被測定空間内の湿度分布を演算によって求めるようにしたので、被測定空間の中央部付近に湿度センサ等を設けることなく、被測定空間内の湿度分布を求めることができる、という優れた効果を有する。
【０１１５】
請求項５記載の発明は、請求項３又は請求項４の発明において、被測定空間を区画する壁体から所定距離離れた位置の温度を求め、壁体から所定距離離れた位置の温度も用いて被測定空間内の温度分布を求めるようにしたので、上記効果に加え、被測定空間内の温度分布をより高精度に求めることができる、という効果を有する。
【０１１６】
請求項６記載の発明は、請求項４の発明において、被測定空間を区画する壁体から所定距離離れた位置の湿度を求め、壁体から所定距離離れた位置の湿度も用いて被測定空間内の湿度分布を求めるようにしたので、上記効果に加え、被測定空間内の湿度分布をより精度良く求めることができる、という効果を有する。
【０１１７】
請求項７記載の発明は、請求項４の発明において、複数箇所の間の音波の減衰量を測定して被測定空間内の音波の減衰率の分布を求めることを、複数種の周波数の音波について各々行い、複数種の周波数の音波について各々求めた被測定空間内の音波の減衰率の分布、被測定空間内の温度分布、及び複数種の周波数に基づいて、被測定空間内の湿度分布を演算によって求めるようにしたので、上記効果に加え、被測定空間内の温度及び湿度が、特定の周波数の音波では減衰率が不感領域となる条件であったとしても、被測定空間内の湿度分布を求めることができる、という効果を有する。
【０１１８】
請求項８記載の発明は、請求項４の発明において、音波の減衰量の測定に用いる音波の周波数を、被測定空間の大きさに応じて変更するようにしたので、上記効果に加え、音波の減衰量を効率良く測定することができる、という効果を有する。
【図面の簡単な説明】
【図１】第１実施形態に係る環境状態測定装置の概略構成を示すブロック図である。
【図２】ノードの配置及び測定対象区間を示す被測定空間の平面図である。
【図３】（Ａ）は第１実施形態に係る温度センサの概略構成図、（Ｂ）は壁面からの距離と被測定空間内の温度との関係を示す線図である。
【図４】第１実施形態に係る伝搬時間・減衰量測定処理を示すフローチャートである。
【図５】第１実施形態に係る温度・湿度分布演算処理を示すフローチャートである。
【図６】被測定空間内のボクセル、音波の伝搬経路がボクセル内を横切っている距離を示す概念図である。
【図７】（Ａ）乃至（Ｃ）は、被測定空間内の温度分布に基づく音波の伝搬経路の補正を説明するための概念図である。
【図８】（Ａ）乃至（Ｃ）は、補正後の音波の伝搬経路に基づく被測定空間内の温度分布の再演算を説明するための概念図である。
【図９】気温が−１０℃の条件での、相対湿度と、各種周波数の音波の距離減衰率と、の関係を示す線図である。
【図１０】気温が０℃の条件での、相対湿度と、各種周波数の音波の距離減衰率と、の関係を示す線図である。
【図１１】気温が＋１０℃の条件での、相対湿度と、各種周波数の音波の距離減衰率と、の関係を示す線図である。
【図１２】気温が＋２０℃の条件での、相対湿度と、各種周波数の音波の距離減衰率と、の関係を示す線図である。
【図１３】気温が＋３０℃の条件での、相対湿度と、各種周波数の音波の距離減衰率と、の関係を示す線図である。
【図１４】被測定空間内の温度分布及び湿度分布の表示例を示すイメージ図である。
【図１５】第２実施形態に係る環境状態測定装置の概略構成を示すブロック図である。
【図１６】（Ａ）は第２実施形態に係る温度センサの概略構成図、（Ｂ）は壁面からの距離と被測定空間内の温度との関係を示す線図である。
【図１７】第２実施形態に係る伝搬時間・減衰量測定処理を示すフローチャートである。
【図１８】第２実施形態に係る温度・湿度分布演算処理を示すフローチャートである。
【符号の説明】
１０ 環境状態測定装置
１２ ノード
１４ ホストコンピュータ
２２ 壁体
２４ 被測定空間
３０ スピーカ
３２ マイクロフォン
３８ データ処理部
４４ 温度センサ
６０ 環境状態測定装置
６２ 湿度センサ
【発明の属する技術分野】
本発明は環境状態測定方法及び環境状態測定装置に係り、特に、空間内の温度分布や湿度分布等の環境状態を測定する環境状態測定方法、及び該環境状態測定方法を適用可能な環境状態測定装置に関する。
【０００２】
【従来の技術】
被空調空間内の温度や湿度を所望の値に制御するためには、被空調空間内の温度、湿度を測定する必要があるが、被空調空間内の各箇所における温度、湿度が一定であることは稀であり、特にドーム球場や劇場等の閉鎖された大空間では、空間内の各箇所における温度、湿度は大きくばらついていることが多い。このため、例えばドーム球場や劇場等の大空間において快適性と省エネルギーの双方を満足する空気調和を行う場合や、室内の温度や湿度を精密に所望の値に制御する場合には、被空調空間内の温度や湿度の分布を測定する必要がある。
【０００３】
しかし、温度センサや湿度センサによって温度や湿度の分布を直接測定しようとすると、センサを、被空調空間の中央部付近を含む被空調空間内にマトリクス状に多数配設する必要があるが、コストが嵩み景観が損なわれると共に被測定空間の利用も制限されるので、上記のように多数のセンサを配設することは現実的ではない。またドーム球場や劇場等では、被測定空間の周縁部（例えば壁体等）に取付けたセンサにより周縁部の温度等を測定すると共に、風船等でセンサを吊り上げることで被測定空間の中央部付近の温度等を測定することが考えられるが、この方法では、被測定空間内で例えば野球等のイベントが実施されている最中に温度や湿度の分布を測定することは不可能である。
【０００４】
一方、医学の分野では、人体に対して多方向からＸ線を照射し、人体を透過したＸ線量を高感度のセンサで測定し、測定されたＸ線量に基づき人体の横断面に沿った各位置における体内組織のＸ線吸収量をコンピュータにより演算して２次元画像として再構成し、人体の横断面像としてモニタに表示するコンピュータ断層撮影（ＣＴ：Ｃｏｍｐｕｔｅｄ Ｔｏｍｏｇｒａｐｈｙ）技術が広く利用されている。また、Ｘ線に代えて超音波を用いることも知られている。
【０００５】
上記のＣＴ技術によれば、測定対象（この場合は人体）内部の状態を非接触で測定することができるので、このＣＴ技術を利用して被測定空間内の温度分布を求めることが提案されている。具体的には、Ｘ線に代えて、媒質の温度によって伝搬時間（速度）が変化する音波を用い、例えば被測定空間を区画している壁体の多数箇所にスピーカ、マイクロフォン、温度センサを設け、スピーカ及びマイクロフォンによって前記多数箇所の間の音波の伝搬時間を各々測定し、被測定空間内の多数箇所の間の音波の伝搬時間から、被測定空間内の温度分布を演算によって求めている。
【０００６】
【発明が解決しようとする課題】
しかしながら、上記のＣＴ技術を利用した温度分布の測定は、音波の伝搬時間のみから温度分布を求めているので、実際の温度分布との誤差が大きい、という問題があった。実際の温度分布との誤差を小さくするためには、例えば被測定空間の周縁部に温度センサを設け、被測定空間内の多数箇所の間の音波の伝搬時間に加えて、温度センサによって測定した被測定空間の周縁部の温度も用いて被測定空間の温度分布を求めることが考えられる。
【０００７】
しかし、被測定空間が壁体によって区画されている場合、該壁体の温度は被測定空間内の温度と必ずしも一致しておらず、被測定空間内のうち壁体近傍の温度は、壁体の温度の影響を受けて大きく変化していることが一般的である。このため、被測定空間の壁体の温度を用いて被測定空間内の温度分布を求めたとしても、被測定空間内の温度分布の測定精度の向上には繋がらない、という問題があった。
【０００８】
また、上記のＣＴ技術を利用した温度分布の測定では、音波が直線的に伝搬するものとして温度分布を求めているが、音波はＸ線とは異なり直線的に伝搬するとは限らず、伝搬経路上に存在している空気の温度の影響を受けて音波の伝搬経路自体も変化し、これに伴って、スピーカ及びマイクロフォンによって測定される音波の伝搬時間も変化する。従って、被測定空間内の温度分布に依存する音波の伝搬経路の変化も、被測定空間内の温度分布の測定精度低下の一因となっていた。
【０００９】
更に、湿度分布については測定方法が確立されておらず、被測定空間内の中央部付近に湿度センサ等を設けて中央部付近の湿度を測定することなく、被測定空間内の湿度分布を求めることは不可能であった。
【００１０】
本発明は上記事実を考慮して成されたもので、被測定空間内の温度分布を高精度に求めることができる環境状態測定方法及び環境状態測定装置を得ることが目的である。
【００１１】
また本発明は、被測定空間内の湿度分布を、被測定空間の中央部付近に湿度検出手段を設けることなく求めることができる環境状態測定方法及び環境状態測定装置を得ることが目的である。
【００１２】
【課題を解決するための手段】
上記目的を達成するために請求項１記載の発明に係る環境状態測定方法は、被測定空間の周縁部の互いに異なる複数箇所の間の音波の伝搬時間を各々測定すると共に、前記被測定空間内でかつ前記被測定空間を区画する壁体から所定距離離れた位置の温度を測定し、前記複数箇所の間の音波の伝搬時間、及び前記壁体から所定距離離れた位置の温度に基づいて、前記被測定空間内の温度分布を演算により求める。
【００１３】
請求項１記載の発明では、被測定空間の周縁部の互いに異なる複数箇所の間の音波の伝搬時間を各々測定すると共に、被測定空間内でかつ被測定空間を区画する壁体から所定距離離れた位置の温度を測定している。なお、複数箇所の間の音波の伝搬時間は、例えば複数箇所の各々にスピーカ等の音波発生手段と、マイクロフォン等の音波検知手段を設け、或る箇所で音波を発生させてから他の箇所で音波が検知される迄の時間（音波の伝搬時間）を測定することを、前記複数箇所の間で各々行うことで得ることができる。
【００１４】
また、壁体から所定距離離れた位置の温度は、壁体から所定距離離れた位置に配設した温度センサ等の温度測定手段によって測定することができる。この所定距離としては、被測定空間を区画する壁体の温度の影響が十分に小さくなるか、又は前記影響が殆ど無視できる大きさとなる距離を適用することができ、例えば壁体からの距離と温度との関係を測定して定めることができる。これにより、壁体の温度の影響が非常に小さく、又は排除された被測定空間内の温度を得ることができる。
【００１５】
また、請求項１の発明では、上記のようにして測定した複数箇所の間の音波の伝搬時間、及び壁体から所定距離離れた位置の温度に基づいて、被測定空間内の温度分布を演算によって求めている。これは、例えば被測定空間を仮想的に多数の領域（ボクセル）に区切り、或る一対の箇所の間を音波が伝搬した際に前記一対の箇所の間に存在する各領域を音波が横切る距離を求め、前記一対の箇所の間の音波の伝搬時間と、前記各領域を音波が横切る際の伝搬速度との関係を、前記各領域を音波が横切る距離を用いて数式で表すことを前記複数箇所の間について各々行って連立方程式を立て、測定した複数箇所の間の音波の伝搬時間を連立方程式に代入すると共に、前記壁体から所定距離離れた位置に対応する領域については、測定した温度に基づき音波の伝搬速度を演算して連立方程式に代入した後に、伝搬速度が未知の領域について音波が横切る際の伝搬速度を演算し、伝搬速度から各領域の温度を演算することによって実現できる。
【００１６】
このように、請求項１の発明では、壁体から所定距離離れた位置の温度を測定し、複数箇所の間の音波の伝搬時間に加えて、壁体から所定距離離れた位置の温度の実測値を用いて被測定空間内の温度分布を演算しており、被測定空間内の温度分布の演算における未知数が減少すると共に、壁体の温度の影響を非常に小さく又は排除することができるので、被測定空間内の温度分布を精度良く求めることができる。
【００１７】
請求項２記載の発明に係る環境状態測定方法は、被測定空間を区画する壁体と、前記被測定空間内でかつ前記壁体から所定距離離れた位置と、の温度差を予め測定しておき、前記被測定空間の周縁部の互いに異なる複数箇所の間の音波の伝搬時間を各々測定すると共に、前記壁体の温度を測定し、測定した壁体の温度及び前記予め測定した温度差に基づいて、前記壁体から所定距離離れた位置の温度を推定し、前記複数箇所の間の音波の伝搬時間、及び前記壁体から所定距離離れた位置の温度に基づいて、前記被測定空間内の温度分布を演算により求める。
【００１８】
請求項２記載の発明では、被測定空間を区画する壁体と、被測定空間内でかつ壁体から所定距離離れた位置と、の温度差を予め測定しておき、壁体の温度を測定し、測定した壁体の温度及び予め測定した温度差に基づいて、壁体から所定距離離れた位置の温度を推定している。これにより、壁体の温度の影響が非常に小さく、又は排除された被測定空間内の温度を得ることができる。
【００１９】
そして請求項２の発明では、被測定空間の周縁部の互いに異なる複数箇所の間の音波の伝搬時間を各々測定し、複数箇所の間の音波の伝搬時間、及び壁体から所定距離離れた位置の温度に基づいて、被測定空間内の温度分布を演算により求めており、請求項１の発明と同様に、被測定空間内の温度分布の演算における未知数が減少すると共に、壁体の温度の影響を小さく又は排除することができるので、被測定空間内の温度分布を精度良く求めることができる。
【００２０】
また、請求項２の発明では、壁体から所定距離離れた位置の温度として推定値を用いているので、請求項１の発明と比較すると温度分布の演算精度は若干低下するものの、壁体から所定距離離れた位置に温度検出手段を設ける必要がないので、温度測定に際しての制約が少なくなる。
【００２１】
請求項３記載の発明に係る環境状態測定方法は、被測定空間の周縁部の互いに異なる複数箇所の間の音波の伝搬時間を各々測定し、前記複数箇所の間の音波の伝搬時間に基づいて、前記複数箇所の間の音波の伝搬経路を直線と仮定して前記被測定空間内の温度分布を演算により求め、求めた温度分布に基づき前記複数箇所の間の音波の伝搬経路を補正し、前記複数箇所の間の音波の伝搬時間、及び前記補正によって得られた複数箇所の間の音波の伝搬経路に基づいて、前記被測定空間内の温度分布を演算により再度求める。
【００２２】
被測定空間内における音波の伝搬経路は、被測定空間内の温度分布に依存して変化するが、請求項３記載の発明では、被測定空間の周縁部の互いに異なる複数箇所の間の音波の伝搬時間を各々測定し、複数箇所の間の音波の伝搬経路を直線と仮定して被測定空間内の温度分布を演算により求め、求めた温度分布に基づき複数箇所の間の音波の伝搬経路を補正している。これにより、実際の音波の伝搬経路により近い伝搬経路を得ることができる。
【００２３】
そして、請求項３の発明では、測定した複数箇所の間の音波の伝搬時間及び補正によって得られた複数箇所の間の音波の伝搬経路に基づいて、被測定空間内の温度分布を演算により再度求めているので、被測定空間内の温度分布を精度良く求めることができる。なお請求項３の発明において、音波の伝搬経路を補正し、被測定空間内の温度分布を再度求めることを、複数回繰り返すようにすれば、被測定空間内の温度を更に精度良く求めることができる。
【００２４】
ところで、本願発明者等は被測定空間内の湿度分布を測定するために、湿度に依存して変化する物理量について検討を行った。その結果、空気中を伝搬する音波の減衰率が、空気の温度、湿度、及び音波の周波数によって変化することを見い出し、被測定空間内の音波の減衰率の分布、被測定空間内の温度分布、及び音波の周波数から、被測定空間内の湿度分布を求めることができる、との知見を得た。
【００２５】
このため、請求項４記載の発明に係る環境状態測定方法は、被測定空間の周縁部の互いに異なる複数箇所の間の音波の伝搬時間及び減衰量を各々測定し、前記複数箇所の間の音波の伝搬時間に基づいて、前記被測定空間内の温度分布を演算により求めると共に、前記複数箇所の間の音波の減衰量に基づいて、前記被測定空間内の音波の減衰率の分布を演算によって求め、前記被測定空間内の音波の減衰率の分布、前記被測定空間内の温度分布、及び前記減衰量の測定に用いた音波の周波数に基づいて、前記被測定空間内の湿度分布を演算によって求める。
【００２６】
請求項４記載の発明では、被測定空間の周縁部の互いに異なる複数箇所の間の音波の伝搬時間及び減衰量を各々測定している。なお、複数箇所の間の音波の減衰量は、音波の伝搬時間の測定と同様に、例えば複数箇所の各々にスピーカ等の音波発生手段と、マイクロフォン等の音波検知手段を設け、或る箇所で発生させた所定の音量の音波が、他の箇所でどの程度の音量で検知されるかを測定し、音波発生時の音量、音波検知時の音量、音波発生手段及び音波検知手段の効率に基づいて音波の減衰量を求めることを、前記複数箇所の間で各々行うことによって得ることができる。
【００２７】
また請求項４の発明では、複数箇所の間の音波の伝搬時間に基づいて被測定空間内の温度分布を演算により求めると共に、複数箇所の間の音波の減衰量に基づいて被測定空間内の音波の減衰率の分布を演算によって求めている。被測定空間内の音波の減衰率の分布は、先に述べた被測定空間内の温度分布と同様に、例えば被測定空間を仮想的に多数の領域（ボクセル）に区切り、或る一対の箇所の間を音波が伝搬した際に前記一対の箇所の間に存在する各領域を音波が横切る距離を求め、前記一対の箇所の間を音波が伝搬した際の減衰量と、前記各領域における音波の減衰率との関係を、前記各領域を音波が横切る距離を用いて数式で表すことを前記複数箇所の間について各々行って連立方程式を立て、測定した複数箇所の間の音波の減衰量を連立方程式に代入した後に、各領域における音波の減衰率を演算することによって求めることができる。
【００２８】
そして請求項４の発明では、被測定空間内の音波の減衰率の分布、被測定空間内の温度分布、及び減衰量の測定に用いた音波の周波数に基づいて、被測定空間内の湿度分布を演算によって求めている。従って請求項４の発明によれば、従来は測定方法が確立されていなかった被測定空間内の湿度分布を求めることを、被測定空間の中央部付近に湿度センサ等の湿度検出手段を設けることなく実現することができる。
【００２９】
なお、請求項３又は請求項４の発明における被測定空間内の温度分布の演算に際しては、請求項５に記載したように、被測定空間内でかつ被測定空間を区画する壁体から所定距離離れた位置の温度を求め、被測定空間内の温度分布を、前記壁体から所定距離離れた位置の温度も用いて演算により求めることが好ましい。壁体から所定距離離れた位置の温度は、請求項１のように測定によって得るようにしてもよいし、請求項２のように壁体と壁体から所定距離離れた位置の温度差を予め測定しておき、壁体の温度を測定し、測定した壁体の温度及び予め測定した温度差に基づいて推定するようにしてもよい。これにより、被測定空間内の温度分布をより高精度に求めることができる。また請求項４の発明では、高精度に求めた温度分布に基づいて、被測定空間内の湿度分布もより高精度に求めることができる。
【００３０】
請求項６記載の発明は、請求項４の発明において、被測定空間内でかつ前記被測定空間を区画する壁体から所定距離離れた位置の湿度を求め、前記被測定空間内の湿度分布を、前記壁体から所定距離離れた位置の湿度も用いて演算により求めることを特徴としている。
【００３１】
前述のように、被測定空間内のうち壁体の温度は被測定空間内の温度と必ずしも一致しておらず、被測定空間内の壁体近傍では、壁体の温度の影響を受けて温度が大きく変化していることが一般的であるが、これに伴って、前記壁体近傍では湿度も大きく変化していることが多い。これに対し請求項６の発明では、被測定空間内でかつ被測定空間を区画する壁体から所定距離離れた位置の湿度を求め、被測定空間内の湿度分布を、壁体から所定距離離れた位置の湿度も用いて演算により求めており、被測定空間内の湿度分布の演算における未知数が減少すると共に、壁体の温度による影響を小さく又は排除することができるので、被測定空間内の湿度分布をより精度良く求めることができる。
【００３２】
なお、壁体から所定距離離れた位置の湿度は、請求項１と同様に測定によって得るようにしてもよいし、請求項２と同様に、壁体と壁体から所定距離離れた位置の湿度差を予め測定しておき、壁体の湿度を測定し、測定した壁体の湿度及び予め測定した湿度差に基づいて推定するようにしてもよい。
【００３３】
また本願発明者等は、空気中を伝搬する音波の減衰率と、空気の温度、湿度、及び音波の周波数との関係について考察し、空気の温度及び湿度が或る範囲内の場合には、湿度の変化に拘らず音波の減衰率が変化しない不感領域があること、及びこの不感領域が生ずるときの空気の温度及び湿度は、音波の周波数によって異なっている、との知見を得た。
【００３４】
上記に基づき請求項７記載の発明は、請求項４の発明において、前記複数箇所の間の音波の減衰量を測定し、前記被測定空間内の音波の減衰率の分布を演算によって求めることを、複数種の周波数の音波について各々行い、前記複数種の周波数の音波について各々求めた前記被測定空間内の音波の減衰率の分布、前記被測定空間内の温度分布、及び前記複数種の周波数に基づいて、前記被測定空間内の湿度分布を演算によって求めることを特徴としている。
【００３５】
請求項７記載の発明では、前記複数箇所の間の音波の減衰量を測定し、前記被測定空間内の音波の減衰率の分布を演算によって求めることを、複数種の周波数の音波について各々行い、複数種の周波数の音波について各々求めた前記被測定空間内の音波の減衰率の分布を用いて被測定空間内の湿度分布を演算によって求めている。これにより、被測定空間内の環境条件（温度及び湿度）が、特定の周波数の音波では湿度の変化に拘らず減衰率が変化しない条件（前記特定の周波数の音波では湿度の変化に対して減衰率が不感領域となる条件）であったとしても、別の周波数の音波の減衰率の分布から被測定空間内の湿度分布を求めることができる。
【００３６】
ところで、音波の減衰量の測定に際し、音波の減衰量が非常に大きい場合には、発生させる音波の音量を大きくしたり、高感度の音波検知手段によって音波を検知する必要が有り、音波の減衰量が非常に小さい場合には、高感度の音波検知手段によって音波の音量を精密に検知する必要が有る。そして、或る二点間を伝搬する音波の減衰量は、前記二点間の距離によって大きく変化する。
【００３７】
このため、請求項８記載の発明は、請求項４の発明において、前記複数箇所の間の音波の減衰量の測定に用いる音波の周波数を、被測定空間の大きさに応じて変更することを特徴としている。
【００３８】
音波の減衰率は音波の周波数によって大きく変化し、温度及び湿度が一定であっても音波の周波数が高周波になるに従って減衰率は高くなる。請求項８の発明では、音波の減衰量の測定に用いる音波の周波数を、被測定空間の大きさに応じて変更しているので、被測定空間の大きさ、すなわち複数箇所の間の距離に拘らず、複数箇所の間を伝搬する音波の減衰量が、測定に好適な値となるように調整することができる。これにより、複数箇所の間の音波の減衰量を効率良く測定することができる。
【００３９】
なお、請求項８の発明における音波の周波数の変更は、具体的には、被測定空間の大きさが大きくなるに従って、すなわち複数箇所の間の距離が大きくなるに従って、周波数が低くなるように変更することができる。
【００４０】
請求項９記載の発明に係る環境状態測定装置は、被測定空間の周縁部の互いに異なる複数箇所の間の音波の伝搬時間を各々測定する伝搬時間測定手段と、前記被測定空間内でかつ前記被測定空間を区画する壁体から所定距離離れた位置の温度を測定する温度測定手段と、前記伝搬時間測定手段によって測定された前記複数箇所の間の音波の伝搬時間、及び前記温度測定手段によって測定された前記壁体から所定距離離れた位置の温度に基づいて、前記被測定空間内の温度分布を演算により求める温度分布演算手段と、を含んで構成しているので、請求項１の発明と同様に、被測定空間内の温度分布を高精度に求めることができる。
【００４１】
請求項１０記載の発明に係る環境状態測定装置は、予め測定された、被測定空間を区画する壁体と、前記被測定空間内でかつ前記壁体から所定距離離れた位置と、の温度差を記憶する記憶手段と、前記被測定空間の周縁部の互いに異なる複数箇所の間の音波の伝搬時間を各々測定する伝搬時間測定手段と、前記壁体の温度を測定する壁体温度測定手段と、前記壁体温度測定手段によって測定された壁体の温度及び前記記憶手段に記憶されている温度差に基づいて、前記壁体から所定距離離れた位置の温度を推定する温度推定手段と、前記伝搬時間測定手段によって測定された前記複数箇所の間の音波の伝搬時間、及び前記温度推定手段によって推定された前記壁体から所定距離離れた位置の温度に基づいて、前記被測定空間内の温度分布を演算により求める温度分布演算手段と、を含んで構成しているので、請求項２の発明と同様に、被測定空間内の温度分布を高精度に求めることができる。
【００４２】
請求項１１記載の発明に係る環境状態測定装置は、被測定空間の周縁部の互いに異なる複数箇所の間の音波の伝搬時間を各々測定する伝搬時間測定手段と、前記伝搬時間測定手段によって測定された前記複数箇所の間の音波の伝搬時間に基づいて、前記複数箇所の間の音波の伝搬経路を直線と仮定して前記被測定空間内の温度分布を演算により求め、求めた温度分布に基づき前記複数箇所の間の音波の伝搬経路を補正し、前記複数箇所の間の音波の伝搬時間、及び前記補正によって得られた複数箇所の間の音波の伝搬経路に基づいて、前記被測定空間内の温度分布を演算により再度求める温度分布演算手段と、を含んで構成しているので、請求項３の発明と同様に、被測定空間内の温度分布を高精度に求めることができる。
【００４３】
請求項１２記載の発明に係る環境状態測定装置は、被測定空間の周縁部の互いに異なる複数箇所の間の音波の伝搬時間及び減衰量を各々測定する伝搬時間・減衰量測定手段と、前記伝搬時間・減衰量測定手段によって測定された前記複数箇所の間の音波の伝搬時間に基づいて、前記被測定空間内の温度分布を演算により求める温度分布演算手段と、前記伝搬時間・減衰量測定手段によって測定された前記複数箇所の間の音波の減衰量に基づいて前記被測定空間内の音波の減衰率の分布を演算によって求め、前記被測定空間内の音波の減衰率の分布及び前記温度分布演算手段によって演算された前記被測定空間内の温度分布に基づいて、前記被測定空間内の湿度分布を演算によって求める湿度分布演算手段と、を含んで構成しているので、請求項４の発明と同様に、被測定空間の中央部付近に湿度センサ等の湿度検出手段を設けることなく、被測定空間内の湿度分布を求めることができる。
【００４４】
【発明の実施の形態】
以下、図面を参照して本発明の実施形態の一例を詳細に説明する。
【００４５】
〔第１実施形態〕
図１には、本第１実施形態に係る環境状態測定装置１０が示されている。環境状態測定装置１０は、多数のノード１２Ａ、１２Ｂ、１２Ｃ、…と、ホストコンピュータ１４を備えており、これらが伝送媒体１６を介して互いに接続されて構成されている。ホストコンピュータ１４は、ＣＰＵ１４Ａ、ＲＯＭ１４Ｂ、ＲＡＭ１４Ｃ、入出力ポート１４Ｄを備え、これらがバス１４Ｅを介して互いに接続されて構成されている。入出力ポート１４Ｄはネットワーク伝送部４９を介して伝送媒体１６に接続されている。また、入出力ポート１４Ｄには、各種の情報を表示するためのディスプレイ１８と、オペレータが各種のデータやコマンド等を入力するためのキーボード２０が接続されている。
【００４６】
一方、多数のノード１２は、被測定空間の周縁部に、周縁部の全周に亘って略一定の間隔で配置される。例として図２には、円筒状の壁体２２によって区画され、上方が屋根（図示省略）によって閉塞されたドーム状の被測定空間２４（例えばドーム球場等）内の環境状態を測定する場合の各ノード１２の配置が示されており、各ノード１２は壁体２２に埋設されている。なお、ノード１２の数及び間隔は図２に示した例に限定されるものではない。また、各ノード１２を識別するために、各ノードには便宜的にノード番号が付与されている。
【００４７】
ノード１２Ａ、１２Ｂ、１２Ｃ、…は各々同一の構成であるので、以下では図１を参照し、ノード１２Ａの構成について説明する。ノード１２Ａは、各々被測定空間２４内に向けて配置された、音波発生手段としてのスピーカ３０、及び音波検知手段としてのマイクロフォン３２を備えている。スピーカ３０は、増幅器３４、Ｄ／Ａ変換器３６を介してデータ処理部３８に接続されており、マイクロフォン３２は、増幅器４０、Ａ／Ｄ変換器４２を介してデータ処理部３８に接続されている。
【００４８】
データ処理部３８はマイクロコンピュータ等を含んで構成され、ネットワーク伝送部４８を介して伝送媒体１６に接続されており、伝送媒体１６及びネットワーク伝送部４８を介してホストコンピュータ１４から、周波数ｆの音波の発生が指示されると、指示された周波数ｆの音波を指示された音量で発生させるためのデータを、所定のタイミングでＤ／Ａ変換器３６に出力する。データ処理部３８から出力されたデータは、Ｄ／Ａ変換器３６でアナログの電気信号に変換され、増幅器３４で増幅された後にスピーカ３０に供給される。これにより、スピーカ３０からは、周波数ｆの音波が所定の時刻に指示された音量でスピーカ３０から発せられる。
【００４９】
また、マイクロフォン３２から出力される電気信号は、増幅器４０で増幅されＡ／Ｄ変換器４２でディジタルデータに変換されてデータ処理部３８に入力される。データ処理部３８は、Ａ／Ｄ変換器４２を介して入力されるデータに基づいて、マイクロフォン３２から出力される電気信号のレベルの変化を監視し、閾値以上の音量の音波が到来したか否か判断する。そして、閾値以上の音量の音波が到来したと判断した場合には、到来した音波の音量及び音波の到来時刻（検知時刻）をメモリ等に記憶する。そして、ネットワーク伝送部４８及び伝送媒体１６を介してホストコンピュータ１４から音波の検知時刻及び音量の問い合わせがあった場合には、記憶している音波の検知時刻及び音量を表すデータをホストコンピュータ１４に送信する。
【００５０】
また、ノード１２は熱電対等から成る温度センサ４４を備えている。被測定空間２４内の温度は壁体２２の温度と異なっており、壁体２２の近傍では、壁体２２の壁面に近づくに従って壁体２２の温度の影響をより強く受けるので、図３（Ｂ）に示すように、被測定空間２４内の温度は壁体２２の近傍で大きく変化している。このため、本実施形態では、壁体２２の壁面と、壁体２２の温度の影響が略０となる位置との距離Ｌ１を予め求め、図３（Ａ）に示すように、温度センサ４４が、壁体２２の壁面から被測定空間２４の中央部に向けて距離Ｌ１だけ隔てた位置の温度を測定するように構成されている。
【００５１】
温度センサ４４は請求項９に記載の温度測定手段に対応しており、温度センサ４４を駆動するドライバ４６を介してデータ処理部３８に接続されている。データ処理部３８は、伝送媒体１６及びネットワーク伝送部４８を介してホストコンピュータ１４から温度の測定が指示されると、壁体２２の壁面から距離Ｌだけ隔てた位置の温度を温度センサ４４によって測定し、測定結果を表す温度データをホストコンピュータ１４に送信する。
【００５２】
次に本第１実施形態の作用として、まず図４のフローチャートを参照し、被測定空間２４内の環境状態（温度分布及び湿度分布）を求める場合にホストコンピュータ１４で実行される伝搬時間・減衰量測定処理について説明する。
【００５３】
ステップ１００では、被測定空間２４の大きさに基づき、測定に用いる音波の周波数を複数種設定する。具体的には、複数種の周波数（ｆ１ 、ｆ２ 、…）として、スピーカ３０が音波を発生可能でマイクロフォン３２が音波を検知可能な周波数帯域内の値を設定すると共に、被測定空間２４の大きさが大きくなるに従って（すなわち被測定空間１２の中央部を挟んで対向するノード１２間の距離が大きくなるに従って）、前記複数種の周波数が全体的に低周波側に偏倚するように設定する。上記の周波数の設定は請求項８の発明に対応している。
【００５４】
次のステップ１０２では、ステップ１００で設定した複数種の周波数（ｆ１ 、ｆ２ 、…）の中から、測定未実行の周波数ｆを選択する。ステップ１０４ではカウンタｍに１を代入し、ステップ１０６では、ノード番号ｍのノード１２（以下、単にノードｍという）に対し、スピーカ３０から所定の時刻に周波数ｆの音波を所定の音量で発生させるよう指示する。これにより、ノードｍのスピーカ３０からは、所定の時刻に周波数ｆの音波が所定の音量で発生され、発生された音波は他のノード１２のマイクロフォン３２で各々検知され、音波の検知時刻及び音量が各々メモリ等に記憶される。
【００５５】
ところで本実施形態では、音波を発生させるノードｍに対し、残りのノードのうちノードｍから比較的離れた位置に存在している所定数のノード１２のみを測定対象のノードとしている（図２及び図６には、各ノードから音波を発生させるときの測定対象のノードを線で結んで示す（この線は音波の伝搬時間及び総減衰量の測定対象区間を表している））。このため、ステップ１０８では、ノードｍに対応する測定対象の所定数のノード１２のノード番号（ｎ１ 、ｎ２ 、…）を各々判断し、次のステップ１１０では、判断したノード番号に基づいて、測定対象の所定数のノード１２に対して音波の検知時刻及び音量を各々問い合わせ、測定対象の所定数のノード１２から送信された音波の検知時刻及び音量を表すデータを各々取り込む。
【００５６】
ステップ１１２では、ノードｍから音波が発生された時刻及び発生された音波の音量、測定対象の所定数のノード１２の各々における音波の検知時刻及び音量に基づいて、ノードｍのスピーカ３０から測定対象のノードｎのマイクロフォン３２への音波（周波数ｆの音波）の伝搬時間ｔｍｎ〔秒〕、ノードｍのスピーカ３０から測定対象のノードｎのマイクロフォン３２へ周波数ｆの音波が伝搬する際の総減衰量Ｒｍｎｆ 〔ｄＢ〕を、測定対象の所定数のノード１２について各々演算し、周波数ｆと対応させて記憶する。
【００５７】
次のステップ１１４では、カウンタｍの値が最終ノード番号に一致したか否か判定する。判定が否定された場合にはステップ１１６へ移行し、カウンタｍの値を１だけインクリメントしてステップ１０６に戻る。これにより、ステップ１１４の判定が肯定される迄の間は、音波を発するノード１２を順に切り換えてステップ１０６〜１１２の処理が繰り返され、例として図２に示す全ての測定対象区間について、音波の検知時刻及び音量の測定、音波の伝搬時間ｔｍｎ及び総減衰量Ｒｍｎｆ の演算が行われることになる。
【００５８】
上記のステップ１０６〜１１６は、各ノードのスピーカ３０、マイクロフォン３２、データ処理部３８と共に、請求項９乃至請求項１１に記載の伝搬時間測定手段、請求項１２に記載の伝搬時間・減衰量測定手段に対応している。
【００５９】
なお、上記では全てのノード１２から各々音波を発生させ、測定対象の所定数のノードから音波の検知時刻及び音量を毎回取り込んで音波の伝搬時間ｔｍｎ及び総減衰量Ｒｍｎｆ を演算しており、例として図２に示す全ての測定対象区間の各々に対し、２方向（一方のノードから他方のノードへ向かう方向、及び他方のノードから一方のノードに向かう方向）について音波の検知時刻及び音量の測定、音波の伝搬時間ｔｍｎ及び総減衰量Ｒｍｎｆ の演算を行っているが、同一の区間の前記２方向についての音波の伝搬時間及び総減衰量は一般に等しいので、各区間に対し何れか一方の方向についてのみ、音波の検知時刻及び音量の測定、音波の伝搬時間ｔｍｎ及び総減衰量Ｒｍｎｆ の演算を行うようにしてもよい。
【００６０】
これにより、図２の例では１個のノードにはマイクロフォン３２を設ける必要がなくなり、別の５個のノードにはスピーカ３０を設ける必要がなくなるので、構成を簡単にすることができると共に、図４に示した伝搬時間・減衰量測定処理に要する時間も短縮することができる。
【００６１】
ステップ１１４の判定が肯定されるとステップ１１８に移行し、先のステップ１００で設定した複数種の周波数の各々について、全ての測定対象区間に対する音波の検知時刻及び音量の測定、及び音波の伝搬時間ｔｍｎ及び総減衰量Ｒｍｎｆ の演算を行ったか否か判定する。判定が否定された場合にはステップ１０２に戻り、複数種の周波数の中から測定未実行の別の周波数ｆを選択し、各ノードから選択した周波数ｆの音波を発生させて、全ての測定対象区間に対する音波の検知時刻及び音量の測定、及び音波の伝搬時間ｔｍｎ及び総減衰量Ｒｍｎｆ の演算を行う。これにより、全ての測定対象区間における音波の伝搬時間及び総減衰量が、ステップ１００で選択した複数種の周波数の音波について各々求められることになる。
【００６２】
なお、音速は温度にのみ依存し、同一の区間を伝搬する音波の伝搬時間は、音波の周波数に拘らず等しいので、音波の伝搬時間は単一の周波数の音波についてのみ測定・演算を行い、音波の総減衰量は複数種の周波数について測定・演算を行うようにしてもよい。
【００６３】
全ての測定対象区間における音波の伝搬時間及び総減衰量が、複数種の周波数について各々求められると、ステップ１１８の判定が肯定されてステップ１２０へ移行する。ステップ１２０では、全てのノード１２に対し温度センサ４４による温度の測定を指示する。そして次のステップ１２２では、全てのノード１２からの温度データの取込みを行う。
【００６４】
上記のようにして伝搬時間・減衰量測定処理が行われると、続いてホストコンピュータ１４では温度・湿度分布演算処理が実行される。この温度・湿度分布演算処理について、図５のフローチャートを参照して説明する。
【００６５】
ステップ１５０では、複数種の周波数について各々測定・演算した音波の伝搬時間のうち、特定の周波数について測定・演算した全ての測定対象区間における音波の伝搬時間ｔｍｎを取り込む。
【００６６】
ところで、本実施形態では図６に示すように、被測定空間２４内を仮想的に多数の矩形状の領域（ボクセル）に区切り、各ボクセル毎に温度及び湿度を演算する。また本実施形態では、各測定対象区間における音波の伝搬経路を各々直線と仮定したときの、音波の伝搬経路が該伝搬経路上に存在する各ボクセルを横切っている距離ｌｍｎｘｙ〔ｍ〕（図６参照、但しｍは測定対象区間の一端の音波発生側のノードのノード番号、ｎは測定対象区間の他端の音波検知側のノードのノード番号、ｘ及びｙは各ボクセルを識別する符号である）が、各測定対象区間について予め演算されてＲＯＭ１４Ｂ等に記憶されている。
【００６７】
次のステップ１５２では、音波の伝搬経路が該伝搬経路上に存在する各ボクセルを横切っている距離ｌｍｎｘｙをＲＯＭ１４Ｂから取り込む。なお、距離ｌｍｎｘｙはボクセルの大きさ及び位置、ノードｍ及びノードｎの位置から幾何学的に算出可能である。
【００６８】
次のステップ１５４では、ステップ１５２で取り込んだ距離ｌｍｎｘｙを用い、全ての測定対象区間について、音波の伝搬時間ｔｍｎと、音波の伝搬経路上に存在する各ボクセルを音波が通過する際の音速ｖｘｙ〔ｍ／秒〕との関係を表す方程式を各々生成することにより、被測定空間２４内の各ボクセルの温度Ｔｘｙを各々求めるための連立方程式を生成する。以下、この連立方程式について説明する。
【００６９】
音速ｖは次の（１）式に示すように、音波が伝搬する媒体（空気）の温度Ｔにのみ依存することが知られている。
【００７０】
ｖ＝ ３３１．４５ ＋ ０．６０７・Ｔ 〔ｍ／秒〕 …（１）
ノードｍのスピーカ３０からノードｎのマイクロフォン３２に至る測定対象区間を音波が伝搬する際に、音波の伝搬経路上に存在するボクセルｘｙを音波が通過する時間ｔｍｎｘｙ〔秒〕は次の（２）式により表され、この（２）式を（１）式に代入すると次の（３）式が得られる。
【００７１】
ｔｍｎｘｙ＝ｌｍｎｘｙ／ｖｘｙ …（２）
ｔｍｎｘｙ＝ｌｍｎｘｙ／（３３１．４５＋ ０．６０７・Ｔｘｙ） …（３）
但し、Ｔｘｙはボクセルｘｙ内の温度〔℃〕である。ノードｍのスピーカ３０からノードｎのマイクロフォン３２に至る音波の伝搬経路は多数のボクセルを横切っており、スピーカｍから前記経路を伝搬してマイクロフォンｎに至る音波の伝搬時間ｔｍｎは、前記多数のボクセルの各々における音波の通過時間の総和であるので、
Ｑｘｙ＝１／ｖｘｙ＝１／（３３１．４５＋ ０．６０７・Ｔｘｙ） …（５）
とすれば、先の（４）式から、ノードｍのスピーカ３０からノードｎのマイクロフォン３２に至る音波の伝搬経路が横切っている各ボクセル内の温度の関数であるＱｘｙ（各ボクセルを音波が通過する際の音速ｖｘｙの逆数）の一次多項式である（６）式が得られる。
【００７２】
ｔｍｎ＝Σ（ｌｍｎｘｙ・Ｑｘｙ） …（６）
各測定対象区間における音波の伝搬時間ｔｍｎ及び距離ｌｍｎｘｙは既知であるので、ステップ１５４では、各測定対象区間における音波の伝搬時間ｔｍｎ及び距離ｌｍｎｘｙを（６）式に各々代入することにより、測定対象区間の数と同数の１次多項式（連立方程式）を生成する。
【００７３】
なお、上記の連立方程式を解くことで、各ボクセルについて変数Ｑｘｙの値を得ることができ、変数Ｑｘｙを（５）式に代入することで各ボクセル内の温度Ｔｘｙを逆算により求めることができるが、被測定空間２４内の周縁部付近に位置している一部のボクセルについては、温度センサ４４によって温度が実測されている。このためステップ１５６では、温度が既知のボクセルについて、各ノードから取込んだ温度データが表す温度Ｔｘｙを（５）式に代入して変数Ｑｘｙの値を演算し、演算したＱｘｙの値を連立方程式に代入する。
【００７４】
そして次のステップ１５８では、温度が既知のボクセルについて変数Ｑｘｙの値を代入した連立方程式を解いて、温度が未知のボクセルの変数Ｑｘｙの値を求め、求めた変数Ｑｘｙの値を（５）式に代入することで各ボクセル内の温度Ｔｘｙを逆算によって求める。これにより、被測定空間２４内の各ボクセル毎の温度、すなわち被測定空間２４内の温度分布が演算により求まることになる。
【００７５】
上記では温度センサ４４によって壁体２２の壁面から距離Ｌ１だけ隔てた位置（壁体２２の温度の影響が略０となる位置）の温度を測定し、測定した温度を連立方程式に代入した後に、連立方程式を解いて各ボクセルの温度を求めており、被測定空間２４の周縁部の温度の実測値を用いて温度分布を求めていると共に、前記実測値は壁体２２の温度の影響が略０となる位置で測定した値であるので、被測定空間２４内の温度分布を精度良く求めることができる。
【００７６】
ステップ１６０では、先のステップ１５８で求めた被測定空間２４内の温度分布に基づき、各測定対象区間における音波の伝搬経路を各々補正する。この音波の伝搬経路の補正は、例えば以下のように行うことができる。
【００７７】
例として図７（Ａ）に示すように被測定空間２４内の温度分布が求まり、測定対象区間の１つであるノードｍとノードｎとの間における音波の伝搬経路を補正する場合、補正前の伝搬経路（図７（Ａ）では直線で示す）に直交する方向に沿った温度変化の勾配を、補正前の伝搬経路の一端から他端に亘って演算する。なお、図７（Ｂ）には、伝搬経路に直交する方向に沿った温度変化の勾配を、勾配が大きくなるに従って長さを長くした太線で示している。
【００７８】
次に、伝搬経路上の各部分における曲率半径を、温度変化の勾配の大きさに反比例するように（温度変化の勾配が大きくなるに従って曲率半径が小さくなるように）決定し、伝搬経路に直交する方向に沿った温度変化における高温側が凸となるように、決定した曲率半径に従って伝搬経路を湾曲させる。図７（Ｃ）に示しているノードｍとノードｎを結ぶ曲線は、決定した曲率半径に従って湾曲させた後の伝搬経路を示している。
【００７９】
続いて、湾曲させた伝搬経路に対し、その両側に、伝搬経路から等しい間隔を隔てて伝搬経路に平行な一対の仮想線（図７（Ｃ）に示す仮想線５０Ａ、５０Ｂ参照）を設定すると共に、伝搬経路及び一対の仮想線を、伝搬経路に沿って所定間隔毎に、伝搬経路に直交する仮想的な区分線（図７（Ｃ）に示す区分線５２参照）によって複数の区間に区切り、各区間において、一対の仮想線上を伝搬する音波の伝搬時間が互いに等しくなるように、伝搬経路の湾曲度合いを更に補正する。
【００８０】
上記の補正処理を、全ての測定対象区間における音波の伝搬経路に対して各々行うと、次のステップ１６２へ移行し、ステップ１６０における伝搬経路の補正において、伝搬経路に対する補正量（例えば補正前の伝搬経路に対する補正後の伝搬経路の距離の最大値や、補正前の伝搬経路と補正後の伝搬経路とで囲まれた部分の面積等）が所定値以上となった伝搬経路が有るか否か判定する。判定が肯定された場合には、ステップ１６４へ移行する。
【００８１】
図８（Ａ）と（Ｂ）を比較しても明らかなように、音波の伝搬経路を補正した場合、伝搬経路上に存在するボクセルが変化したり、伝搬経路が各ボクセルを横切っている距離が変化する。このためステップ１６４では、先のステップ１６０の補正を行った後の全ての測定対象区間における音波の伝搬経路が、該伝搬経路上に存在する各ボクセルを横切っている距離ｌｍｎｘｙを各々演算し、ステップ１５４に戻る。
【００８２】
これにより、ステップ１５４以降では、補正後の音波の伝搬経路に基づいて被測定空間２４内の温度分布が再度演算されるので、例として図８（Ｃ）にも示すように一部のボクセルの温度Ｔｘｙが変化し、各測定対象区間における音波の伝搬経路を直線と仮定して被測定空間２４内の温度分布を演算したことによって生じた誤差が小さくされる。そして、ステップ１５４〜１６４はステップ１６２の判定が肯定されている間繰り返されるので、被測定空間２４内の温度分布を高精度に求めることができる。
【００８３】
なお、上述したステップ１５０〜１６４は、請求項９乃至請求項１２に記載の温度分布演算手段に対応しており、より詳しくは、ステップ１５２、及びステップ１５４〜１５８の処理を第１回目に実行することは請求項１１に記載の温度分布演算手段における「複数箇所の間の音波の伝搬経路を直線と仮定して被測定空間内の温度分布を演算により求め」ることに対応しており、ステップ１６０は、請求項１１に記載の温度分布演算手段における「温度分布に基づき複数箇所の間の音波の伝搬経路を補正」することに対応しており、ステップ１６２の判定が肯定されることにより、ステップ１６４、及びステップ１５４以降の処理を再度実行することは、請求項１１に記載の温度分布演算手段における「複数箇所の間の音波の伝搬時間、及び補正によって得られた複数箇所の間の音波の伝搬経路に基づいて、被測定空間内の温度分布を演算により再度求める」ことに対応している。
【００８４】
ステップ１６２の判定が否定されるとステップ１６６へ移行し、複数種の周波数について各々測定・演算した音波の総減衰量のうち、或る周波数ｆについて測定・演算した全ての測定対象区間における音波の総減衰量Ｒｍｎｆ を取り込む。次のステップ１６８では、ステップ１６６で取り込んだ総減衰量Ｒｍｎｆ を用い、全ての測定対象区間について、音波の総減衰量Ｒｍｎｆ と、音波の伝搬経路上に存在する各ボクセルの周波数ｆの音波についての減衰率Ｒｘｙｆ との関係を表す方程式を各々生成し、被測定空間２４内の各ボクセルの周波数ｆの音波についての減衰率Ｒｘｙｆ を各々求めるための連立方程式を生成する。以下、この連立方程式について説明する。
【００８５】
ノードｍのスピーカ３０からノードｎのマイクロフォン３２に至る測定対象区間を周波数ｆの音波が伝搬する際に、音波の伝搬経路上に存在するボクセルｘｙにおける減衰量Ｒｍｎｘｙｆ は、ボクセルｘｙを周波数ｆの音波が通過する際の単位距離当りの減衰量（減衰率）をＲｘｙｆ 〔ｄＢ／ｍ〕とすると、
Ｒｍｎｘｙｆ ＝ｌｍｎｘｙ・Ｒｘｙｆ …（７）
上記の（７）式で表すことができる。従って、ノードｍのスピーカ３０からノードｎのマイクロフォン３２に至る測定対象区間を周波数ｆの音波が伝搬する際の総減衰量Ｒｍｎｆ 〔ｄＢ〕は、各ボクセルにおける周波数ｆの音波についての減衰率Ｒｘｙｆ の一次多項式である次の（８）式で表すことができる。
【００８６】
【００８７】
各測定対象区間における周波数ｆの音波の総減衰量Ｒｍｎｆ 及び距離ｌｍｎｘｙは既知であり、本実施形態では、変換特性Φｍｎｓｐ及びΦｍｎｍｉｃ は、スピーカ３０やマイクロフォン３２の指向特性や周波数特性に基づき予め算出されてＲＯＭ１４Ｂに記憶されている。従って、ステップ１６８では、各測定対象区間における音波の総減衰量Ｒｍｎｆ 、及び補正後の音波の伝搬経路から先のステップ１６４で演算した距離ｌｍｎｘｙを（８）式に各々代入することにより、測定対象区間の数と同数の１次多項式（連立方程式）を生成する。
【００８８】
次のステップ１７０では、上記の連立方程式を解いて各ボクセルにおける周波数ｆの音波についての減衰率Ｒｘｙｆ を求める。これにより、周波数ｆの音波についての被測定空間２４内の減衰率の分布が演算により求まることになる。次のステップ１７２では複数種の周波数について被測定空間２４内の減衰率の分布を求めたか否か判定する。判定が否定された場合にはステップ１７２に戻り、減衰率の分布を求めていない他の周波数ｆについて測定・演算した音波の総減衰量Ｒｍｎｆ を取り込み、上記と同様にして被測定空間２４内の減衰率の分布を求める。
【００８９】
ステップ１６６〜１７２を複数回繰り返し、複数種の周波数について被測定空間２４内の減衰率の分布を全て求めると、ステップ１７２の判定が肯定されてステップ１７４へ移行する。ステップ１７４では、被測定空間２４内の温度分布、及び複数種の周波数について各々求めた被測定空間２４内の減衰率の分布に基づいて、各ボクセルの湿度を各々演算する。
【００９０】
具体的には、図９〜図１３から明らかなように、空気中を通過する音波の単位距離当りの減衰率（距離減衰率）は、気温及び湿度によって変化すると共に、音波の周波数によっても変化する。また、空気の温度及び音波の周波数によっては、湿度の変化に拘らず減衰率が変化しない不感領域が生ずる。例えば図１０に示す気温が０℃の場合には、周波数 ２５０〔ＨＺ 〕の音波について湿度の変化に拘らず減衰率が殆ど変化しない不感領域が生じ、同様に、図１１に示す気温が１０℃の場合には周波数 ２５０〔ＨＺ 〕及び ５００〔ＨＺ 〕の音波、図１２に示す気温が２０℃の場合には周波数２〔ｋＨＺ 〕の音波、図１３に示す気温が３０℃の場合には周波数４〔ｋＨＺ 〕の音波について湿度の変化に拘らず減衰率が殆ど変化しない不感領域が生じている。
【００９１】
本実施形態では、上記の不感領域の影響で湿度が不定となることを回避するために複数種の周波数について減衰率の分布を求めている。このため、湿度の演算にあたっては、温度、周波数及び減衰率をパラメータ（入力）とし、これらのパラメータに合致する湿度を出力する３次元のルックアップテーブルを作成するか、或いは次の（９）式のように、温度、周波数及び減衰率と湿度との関係を関数式として定めておく。
【００９２】
Ｈｘｙ＝ｆｕｎｃｔｉｏｎ（Ｔｘｙ，ｆ，Ｒｘｙｆ ） …（９）
但し、Ｈｘｙはボクセルｘｙ内の相対湿度〔％〕である。
【００９３】
そして、先に求めた被測定空間２４内の温度分布、先に求めた複数種の周波数についての被測定空間２４内の減衰率の分布から、前述のルックアップテーブル又は関数式を用いて、各ボクセル毎に相対湿度Ｈｘｙを求める。これにより、被測定空間内の湿度分布が求まることになる。
【００９４】
なお、不感領域に関しては、温度、減衰率、及び周波数の各値の組み合わせから、不感領域に相当する条件か否かを判断し、不感領域に相当する条件であると判断したデータを用いないようにするか、或いは入力が不感領域に相当する条件のときには相対湿度データが出力されないようにルックアップテーブルを作成しておけば、不感領域の影響を受けることなく、被測定空間２４内の湿度分布を精度良く求めることができる。
【００９５】
上述したステップ１６６〜１７４は請求項１２に記載の湿度分布演算手段に対応しており、より詳しくは、ステップ１６６〜１７２は、前記湿度分布演算手段における「複数箇所の間の音波の減衰量に基づいて被測定空間内の音波の減衰率の分布を演算によって求め」ることに対応しており、ステップ１７４は、前記湿度分布演算手段における「被測定空間内の音波の減衰率の分布、被測定空間内の温度分布、及び減衰量の測定に用いた音波の周波数に基づいて、被測定空間内の湿度分布を演算によって求める」ことに対応している。また、ステップ１７４において複数種の周波数についての被測定空間２４内の減衰率の分布を用いることは、請求項７記載の発明に対応している。
【００９６】
次のステップ１７６では、上記により得られた被測定空間２４内の温度分布及び湿度分布を、例として図１４に示すコンターマップ等の形態でディスプレイ１８に表示し、処理を終了する。なお、ディスプレイ１８への表示に代えて、被空調空間２４内の空気調和を行う空調装置に対し、被測定空間２４内の温度分布及び湿度分布を表すデータを出力するようにしてもよい。これにより、空調装置が被測定空間２４に対し、快適性と省エネルギーの双方を満足する空気調和を行うことが可能となる。
【００９７】
〔第２実施形態〕
次に本発明の第２実施形態について説明する。なお、第１実施形態と同一の部分には同一の符号を付し、説明を省略する。図１５には本第２実施形態に係る環境状態測定装置６０が示されている。
【００９８】
この環境状態測定装置６０では、各ノード１２が、壁体２２の壁面近傍の位置の湿度を検出する湿度センサ６２を備えている。湿度センサ６２はドライバ６４を介してデータ処理部３８に接続されている。データ処理部３８は、伝送媒体１６及びネットワーク伝送部４８を介してホストコンピュータ１４から湿度の測定が指示されると、壁体２２の壁面近傍の位置の湿度を湿度センサ６２によって測定し、測定結果を表す湿度データをホストコンピュータ１４に送信する。
【００９９】
また、本第２実施形態に係る温度センサ４４は、請求項１０に記載の壁体温度測定手段に対応しており、図１６（Ａ）に示すように、壁体２２の壁面近傍の温度（より詳しくは壁体２２の壁面から被測定空間２４の中央部に向けて距離Ｌ２（Ｌ２≪Ｌ１）だけ隔てた位置の温度）を測定するように構成されている。
【０１００】
本第２実施形態では、温度センサ４４によって測定される壁体２２の壁面から距離Ｌ２だけ隔てた位置の温度と、壁体２２の温度の影響が略０となる位置（例えば前記壁面から距離Ｌ１以上離れた位置）の温度と、の温度差ΔＴが予め測定され、ホストコンピュータ１４のＲＯＭ１４Ｂ（請求項１０に記載の記憶手段に相当）に記憶されている。また、湿度に関しても、湿度センサ６２によって測定される壁体２２の壁面近傍の湿度と、壁体２２の温度の影響が略０となる位置（例えば前記壁面から距離Ｌ１以上離れた位置）の湿度と、の温度差ΔＨが予め測定され、ホストコンピュータ１４のＲＯＭ１４Ｂに記憶されている。
【０１０１】
次に本第２実施形態の作用として、まず図１７のフローチャートを参照し、本第２実施形態に係る伝搬時間・減衰量測定処理について、第１実施形態に係る伝搬時間・減衰量測定処理（図４）と異なる部分についてのみ説明する。
【０１０２】
本第２実施形態に係る伝搬時間・減衰量測定処理では、ステップ１００〜１１８において、複数種の周波数について、全ての測定対象区間に対する音波の検知時刻及び音量の測定、及び音波の伝搬時間ｔｍｎ及び総減衰量Ｒｍｎｆ の演算を行った後に、ステップ１２１において、全てのノード１２に対し温度センサ４４による温度の測定及び湿度センサ６２による湿度の測定を指示する。そして次のステップ１２３では、全てのノード１２からの温度データ及び湿度データの取込みを行う。
【０１０３】
次に図１８のフローチャートを参照し、本第２実施形態に係る温度・湿度分布演算処理について、第１実施形態に係る温度・湿度分布演算処理（図５）と異なる部分についてのみ説明する。
【０１０４】
本第２実施形態に係る温度・湿度分布演算処理では、ステップ１５４で被測定空間２４内の各ボクセルの温度Ｔｘｙを各々求めるための連立方程式を生成した後に、次のステップ１５５において、各ノードから取込んだ温度データを、ＲＯＭ１４Ｂに記憶されている温度差ΔＴにより補正する。このステップ１５５は、請求項１０に記載の温度推定手段に対応している。これにより、壁体２２の温度の影響が略０となるように補正した温度データを得ることができる。そして、次のステップ１５７では、補正後の温度データが表す温度Ｔｘｙを（５）式に代入して変数Ｑｘｙの値を演算し、演算したＱｘｙの値を連立方程式に代入する。
【０１０５】
本第２実施形態では、壁体２２の壁面近傍の温度を測定し、測定によって得られた温度データを、予め測定した温度差ΔＴにより壁体２２の温度の影響が略０となるように補正した後に連立方程式に代入しているので、第１実施形態と同様に被測定空間２４内の温度分布を精度良く求めることができる。また、本第２実施形態によれば、温度センサ４４を壁体２２の壁面から大きく突出させる必要がないので、見栄えが向上する。
【０１０６】
また、本第２実施形態に係る温度・湿度分布演算処理では、ステップ１６６〜１７２で複数種の周波数について被測定空間２４内の減衰率の分布を求めた後に、ステップ１７３において、各ノードから取込んだ湿度データを、ＲＯＭ１４Ｂに記憶されている湿度差ΔＨにより補正する。これにより、壁体２２の温度の影響が略０となるように補正した湿度データを得ることができる。そして、次のステップ１７５では、被測定空間２４内の温度分布、及び複数種の周波数について各々求めた被測定空間２４内の減衰率の分布に基づいて、各ボクセルの湿度を各々演算するが、補正後の湿度データにより相対湿度が既知である壁体２２の壁面近傍のボクセルｘｙについては、補正後の湿度データを相対湿度Ｈｘｙとして設定する。これにより被測定空間２４内の湿度分布が求まる。
【０１０７】
なお、上記のように、補正後の湿度データも用いて被測定空間２４内の湿度分布を求めることは、請求項６の発明に対応している。
【０１０８】
このように、本第２実施形態では、壁体２２の壁面近傍の湿度を測定し、測定によって得られた湿度データを、予め測定した湿度差ΔＨにより壁体２２の温度の影響が略０となるように補正し、補正後の湿度データも用いて被測定空間２４内の湿度分布を求めているので、被測定空間２４内の湿度分布をより精度良く求めることができる。
【０１０９】
なお、第２実施形態では湿度センサ６２によって壁体２２の壁面近傍の湿度を測定し、予め測定した湿度差ΔＨにより壁体２２の温度の影響が略０となるように湿度データを補正していたが、これに限定されるものではなく、第１実施形態で説明した温度センサ４４と同様に、湿度センサ６２によって壁体２２の温度の影響が略０となる位置（例えば壁面から距離Ｌ１だけ隔てた位置）の湿度を測定し、測定した湿度を補正することなく用いて被測定空間２４内の湿度分布を求めるようにしてもよい。
【０１１０】
また、上記では本発明を、円筒状の壁体によって区画され、上方が屋根によって閉塞されたドーム球場等のドーム状の被測定空間２４の温度分布、湿度分布の測定に適用した場合を説明したが、被測定空間の形状や用途等は上記に限定されるものではなく、本発明は、例えば劇場、病院等の空間において快適性と省エネルギーの双方を満足する空気調和を行う場合や、或いはクリーンルーム内の温度及び湿度を精密に所望の値に制御する等の場合に適用可能であることは言うまでもない。
【０１１１】
【発明の効果】
以上説明したように請求項１及び請求項９記載の発明は、被測定空間の周縁部の互いに異なる複数箇所の間の音波の伝搬時間を各々測定すると共に、被測定空間を区画する壁体から所定距離離れた位置の温度を測定し、音波の伝搬時間及び壁体から所定距離離れた位置の温度に基づいて、被測定空間内の温度分布を演算により求めるようにしたので、被測定空間内の温度分布を高精度に求めることができる、という優れた効果を有する。
【０１１２】
請求項２及び請求項１０記載の発明は、被測定空間を区画する壁体と該壁体から所定距離離れた位置との温度差を予め測定し、被測定空間の周縁部の互いに異なる複数箇所の間の音波の伝搬時間を各々測定すると共に、壁体の温度を測定し、測定した壁体の温度及び前記温度差に基づいて壁体から所定距離離れた位置の温度を推定し、音波の伝搬時間及び壁体から所定距離離れた位置の温度に基づいて、被測定空間内の温度分布を演算により求めるようにしたので、被測定空間内の温度分布を高精度に求めることができる、という優れた効果を有する。
【０１１３】
請求項３及び請求項１１記載の発明は、被測定空間の周縁部の互いに異なる複数箇所の間の音波の伝搬時間を各々測定し、複数箇所の間の音波の伝搬時間に基づいて音波の伝搬経路を直線と仮定して被測定空間内の温度分布を演算により求め、求めた温度分布に基づき音波の伝搬経路を補正し、音波の伝搬時間及び補正によって得られた音波の伝搬経路に基づいて、被測定空間内の温度分布を演算により再度求めるようにしたので、被測定空間内の温度分布を高精度に求めることができる、という優れた効果を有する。
【０１１４】
請求項４及び請求項１２記載の発明は、被測定空間の周縁部の互いに異なる複数箇所の間の音波の伝搬時間及び減衰量を各々測定し、音波の伝搬時間に基づいて被測定空間内の温度分布を演算により求めると共に、音波の減衰量に基づいて被測定空間内の音波の減衰率の分布を演算によって求め、音波の減衰率の分布、温度分布、及び減衰量の測定に用いた音波の周波数に基づいて、被測定空間内の湿度分布を演算によって求めるようにしたので、被測定空間の中央部付近に湿度センサ等を設けることなく、被測定空間内の湿度分布を求めることができる、という優れた効果を有する。
【０１１５】
請求項５記載の発明は、請求項３又は請求項４の発明において、被測定空間を区画する壁体から所定距離離れた位置の温度を求め、壁体から所定距離離れた位置の温度も用いて被測定空間内の温度分布を求めるようにしたので、上記効果に加え、被測定空間内の温度分布をより高精度に求めることができる、という効果を有する。
【０１１６】
請求項６記載の発明は、請求項４の発明において、被測定空間を区画する壁体から所定距離離れた位置の湿度を求め、壁体から所定距離離れた位置の湿度も用いて被測定空間内の湿度分布を求めるようにしたので、上記効果に加え、被測定空間内の湿度分布をより精度良く求めることができる、という効果を有する。
【０１１７】
請求項７記載の発明は、請求項４の発明において、複数箇所の間の音波の減衰量を測定して被測定空間内の音波の減衰率の分布を求めることを、複数種の周波数の音波について各々行い、複数種の周波数の音波について各々求めた被測定空間内の音波の減衰率の分布、被測定空間内の温度分布、及び複数種の周波数に基づいて、被測定空間内の湿度分布を演算によって求めるようにしたので、上記効果に加え、被測定空間内の温度及び湿度が、特定の周波数の音波では減衰率が不感領域となる条件であったとしても、被測定空間内の湿度分布を求めることができる、という効果を有する。
【０１１８】
請求項８記載の発明は、請求項４の発明において、音波の減衰量の測定に用いる音波の周波数を、被測定空間の大きさに応じて変更するようにしたので、上記効果に加え、音波の減衰量を効率良く測定することができる、という効果を有する。
【図面の簡単な説明】
【図１】第１実施形態に係る環境状態測定装置の概略構成を示すブロック図である。
【図２】ノードの配置及び測定対象区間を示す被測定空間の平面図である。
【図３】（Ａ）は第１実施形態に係る温度センサの概略構成図、（Ｂ）は壁面からの距離と被測定空間内の温度との関係を示す線図である。
【図４】第１実施形態に係る伝搬時間・減衰量測定処理を示すフローチャートである。
【図５】第１実施形態に係る温度・湿度分布演算処理を示すフローチャートである。
【図６】被測定空間内のボクセル、音波の伝搬経路がボクセル内を横切っている距離を示す概念図である。
【図７】（Ａ）乃至（Ｃ）は、被測定空間内の温度分布に基づく音波の伝搬経路の補正を説明するための概念図である。
【図８】（Ａ）乃至（Ｃ）は、補正後の音波の伝搬経路に基づく被測定空間内の温度分布の再演算を説明するための概念図である。
【図９】気温が−１０℃の条件での、相対湿度と、各種周波数の音波の距離減衰率と、の関係を示す線図である。
【図１０】気温が０℃の条件での、相対湿度と、各種周波数の音波の距離減衰率と、の関係を示す線図である。
【図１１】気温が＋１０℃の条件での、相対湿度と、各種周波数の音波の距離減衰率と、の関係を示す線図である。
【図１２】気温が＋２０℃の条件での、相対湿度と、各種周波数の音波の距離減衰率と、の関係を示す線図である。
【図１３】気温が＋３０℃の条件での、相対湿度と、各種周波数の音波の距離減衰率と、の関係を示す線図である。
【図１４】被測定空間内の温度分布及び湿度分布の表示例を示すイメージ図である。
【図１５】第２実施形態に係る環境状態測定装置の概略構成を示すブロック図である。
【図１６】（Ａ）は第２実施形態に係る温度センサの概略構成図、（Ｂ）は壁面からの距離と被測定空間内の温度との関係を示す線図である。
【図１７】第２実施形態に係る伝搬時間・減衰量測定処理を示すフローチャートである。
【図１８】第２実施形態に係る温度・湿度分布演算処理を示すフローチャートである。
【符号の説明】
１０ 環境状態測定装置
１２ ノード
１４ ホストコンピュータ
２２ 壁体
２４ 被測定空間
３０ スピーカ
３２ マイクロフォン
３８ データ処理部
４４ 温度センサ
６０ 環境状態測定装置
６２ 湿度センサ
Claims (12)
- 被測定空間の周縁部の互いに異なる複数箇所の間の音波の伝搬時間を各々測定すると共に、
前記被測定空間内でかつ前記被測定空間を区画する壁体から所定距離離れた位置の温度を測定し、
前記複数箇所の間の音波の伝搬時間、及び前記壁体から所定距離離れた位置の温度に基づいて、前記被測定空間内の温度分布を演算により求める
環境状態測定方法。 - 被測定空間を区画する壁体と、前記被測定空間内でかつ前記壁体から所定距離離れた位置と、の温度差を予め測定しておき、
前記被測定空間の周縁部の互いに異なる複数箇所の間の音波の伝搬時間を各々測定すると共に、
前記壁体の温度を測定し、測定した壁体の温度及び前記予め測定した温度差に基づいて、前記壁体から所定距離離れた位置の温度を推定し、
前記複数箇所の間の音波の伝搬時間、及び前記壁体から所定距離離れた位置の温度に基づいて、前記被測定空間内の温度分布を演算により求める
環境状態測定方法。 - 被測定空間の周縁部の互いに異なる複数箇所の間の音波の伝搬時間を各々測定し、
前記複数箇所の間の音波の伝搬時間に基づいて、前記複数箇所の間の音波の伝搬経路を直線と仮定して前記被測定空間内の温度分布を演算により求め、
求めた温度分布に基づき前記複数箇所の間の音波の伝搬経路を補正し、
前記複数箇所の間の音波の伝搬時間、及び前記補正によって得られた複数箇所の間の音波の伝搬経路に基づいて、前記被測定空間内の温度分布を演算により再度求める
環境状態測定方法。 - 被測定空間の周縁部の互いに異なる複数箇所の間の音波の伝搬時間及び減衰量を各々測定し、
前記複数箇所の間の音波の伝搬時間に基づいて、前記被測定空間内の温度分布を演算により求めると共に、
前記複数箇所の間の音波の減衰量に基づいて、前記被測定空間内の音波の減衰率の分布を演算によって求め、
前記被測定空間内の音波の減衰率の分布、前記被測定空間内の温度分布、及び前記減衰量の測定に用いた音波の周波数に基づいて、前記被測定空間内の湿度分布を演算によって求める
環境状態測定方法。 - 前記被測定空間内でかつ前記被測定空間を区画する壁体から所定距離離れた位置の温度を求め、
前記被測定空間内の温度分布を、前記壁体から所定距離離れた位置の温度も用いて演算により求める
ことを特徴とする請求項３又は請求項４記載の環境状態測定方法。 - 被測定空間内でかつ前記被測定空間を区画する壁体から所定距離離れた位置の湿度を求め、
前記被測定空間内の湿度分布を、前記壁体から所定距離離れた位置の湿度も用いて演算により求める
ことを特徴とする請求項４記載の環境状態測定方法。 - 前記複数箇所の間の音波の減衰量を測定し、前記被測定空間内の音波の減衰率の分布を演算によって求めることを、複数種の周波数の音波について各々行い、
前記複数種の周波数の音波について各々求めた前記被測定空間内の音波の減衰率の分布、前記被測定空間内の温度分布、及び前記複数種の周波数に基づいて、前記被測定空間内の湿度分布を演算によって求める
ことを特徴とする請求項４記載の環境状態測定方法。 - 前記複数箇所の間の音波の減衰量の測定に用いる音波の周波数を、被測定空間の大きさに応じて変更することを特徴とする請求項４記載の環境状態測定方法。
- 被測定空間の周縁部の互いに異なる複数箇所の間の音波の伝搬時間を各々測定する伝搬時間測定手段と、
前記被測定空間内でかつ前記被測定空間を区画する壁体から所定距離離れた位置の温度を測定する温度測定手段と、
前記伝搬時間測定手段によって測定された前記複数箇所の間の音波の伝搬時間、及び前記温度測定手段によって測定された前記壁体から所定距離離れた位置の温度に基づいて、前記被測定空間内の温度分布を演算により求める温度分布演算手段と、
を含む環境状態測定装置。 - 予め測定された、被測定空間を区画する壁体と、前記被測定空間内でかつ前記壁体から所定距離離れた位置と、の温度差を記憶する記憶手段と、
前記被測定空間の周縁部の互いに異なる複数箇所の間の音波の伝搬時間を各々測定する伝搬時間測定手段と、
前記壁体の温度を測定する壁体温度測定手段と、
前記壁体温度測定手段によって測定された壁体の温度及び前記記憶手段に記憶されている温度差に基づいて、前記壁体から所定距離離れた位置の温度を推定する温度推定手段と、
前記伝搬時間測定手段によって測定された前記複数箇所の間の音波の伝搬時間、及び前記温度推定手段によって推定された前記壁体から所定距離離れた位置の温度に基づいて、前記被測定空間内の温度分布を演算により求める温度分布演算手段と、
を含む環境状態測定装置。 - 被測定空間の周縁部の互いに異なる複数箇所の間の音波の伝搬時間を各々測定する伝搬時間測定手段と、
前記伝搬時間測定手段によって測定された前記複数箇所の間の音波の伝搬時間に基づいて、前記複数箇所の間の音波の伝搬経路を直線と仮定して前記被測定空間内の温度分布を演算により求め、求めた温度分布に基づき前記複数箇所の間の音波の伝搬経路を補正し、前記複数箇所の間の音波の伝搬時間、及び前記補正によって得られた複数箇所の間の音波の伝搬経路に基づいて、前記被測定空間内の温度分布を演算により再度求める温度分布演算手段と、
を含む環境状態測定装置。 - 被測定空間の周縁部の互いに異なる複数箇所の間の音波の伝搬時間及び減衰量を各々測定する伝搬時間・減衰量測定手段と、
前記伝搬時間・減衰量測定手段によって測定された前記複数箇所の間の音波の伝搬時間に基づいて、前記被測定空間内の温度分布を演算により求める温度分布演算手段と、
前記伝搬時間・減衰量測定手段によって測定された前記複数箇所の間の音波の減衰量に基づいて前記被測定空間内の音波の減衰率の分布を演算によって求め、前記被測定空間内の音波の減衰率の分布、前記温度分布演算手段によって演算された前記被測定空間内の温度分布、及び前記伝搬時間・減衰量測定手段が前記減衰量の測定に用いた音波の周波数に基づいて、前記被測定空間内の湿度分布を演算によって求める湿度分布演算手段と、
を含む環境状態測定装置。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP34528397A JP3602951B2 (ja) | 1997-12-15 | 1997-12-15 | 環境状態測定方法及び装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP34528397A JP3602951B2 (ja) | 1997-12-15 | 1997-12-15 | 環境状態測定方法及び装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPH11173925A JPH11173925A (ja) | 1999-07-02 |
| JP3602951B2 true JP3602951B2 (ja) | 2004-12-15 |
Family
ID=18375551
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP34528397A Expired - Fee Related JP3602951B2 (ja) | 1997-12-15 | 1997-12-15 | 環境状態測定方法及び装置 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP3602951B2 (ja) |
Families Citing this family (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| KR100306121B1 (ko) * | 1999-10-21 | 2001-11-14 | 윤덕용 | 음파를 이용한 온도측정시 센서위치의 최적 결정방법 |
| JP4783873B2 (ja) * | 2000-08-11 | 2011-09-28 | 有限会社 サイプラス | 環境状態測定装置及び方法 |
| KR100623634B1 (ko) * | 2003-09-22 | 2006-09-13 | 김형윤 | 구조물의 건전상태 감시방법 |
| WO2021221117A1 (ja) | 2020-04-30 | 2021-11-04 | ダイキン工業株式会社 | 環境検出システム |
| CN116324354A (zh) * | 2020-10-29 | 2023-06-23 | 大金工业株式会社 | 环境状况测量装置 |
| CN112815487B (zh) * | 2021-01-18 | 2022-06-24 | 深圳市晨北科技有限公司 | 加湿器的控制方法及装置、设备及存储介质 |
Family Cites Families (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS55109836U (ja) * | 1979-01-30 | 1980-08-01 | ||
| JPS59162430A (ja) * | 1983-03-07 | 1984-09-13 | Mitsubishi Electric Corp | 超音波による空間温度分布の測定方法およびその表示装置 |
| JPS61265540A (ja) * | 1985-05-20 | 1986-11-25 | Tokyo Electric Power Co Inc:The | 燃焼ガス体温度測定方法 |
| JPH01107942U (ja) * | 1988-01-14 | 1989-07-20 | ||
| JPH01219664A (ja) * | 1988-02-29 | 1989-09-01 | Takasago Thermal Eng Co Ltd | 湿度測定法 |
| JP2524806B2 (ja) * | 1988-04-05 | 1996-08-14 | ニッタン株式会社 | 火災監視装置 |
| JP2686298B2 (ja) * | 1988-12-27 | 1997-12-08 | バブコツク日立株式会社 | ガス温度計測装置と、これを用いたガス分析装置 |
| JPH0418236U (ja) * | 1990-06-07 | 1992-02-14 | ||
| JP3095637B2 (ja) * | 1994-09-28 | 2000-10-10 | 松下電工株式会社 | 超音波湿度センサ |
-
1997
- 1997-12-15 JP JP34528397A patent/JP3602951B2/ja not_active Expired - Fee Related
Also Published As
| Publication number | Publication date |
|---|---|
| JPH11173925A (ja) | 1999-07-02 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP4269623B2 (ja) | 血流可視化診断装置 | |
| US6230042B1 (en) | Method and arrangement for determining the location of a catheter within an animal body | |
| JP5314322B2 (ja) | ボリュメトリックフローを計測するための方法及び装置 | |
| JPS62164400A (ja) | 電子消音システム | |
| EP1223864A1 (en) | A method and apparatus for improving the accuracy with which the speed of a fluid is measured | |
| JP3602951B2 (ja) | 環境状態測定方法及び装置 | |
| CN108646221A (zh) | 一种基于声源阵列的空间麦克风定位方法 | |
| JP2009207899A (ja) | 超音波映像を処理する超音波システム及び方法 | |
| JP4783873B2 (ja) | 環境状態測定装置及び方法 | |
| JPH0595946A (ja) | 超音波の音速測定法及び音速測定手段を備えた超音波診断装置 | |
| Henry et al. | Localization of adventitious respiratory sounds | |
| JP3699269B2 (ja) | 環境状態測定方法及び装置 | |
| US4082088A (en) | Apparatus for determination of respiratory passageway resistance | |
| JP2020189094A (ja) | Ｌａｔ体積マップ | |
| CN113167660A (zh) | 温度测量装置、测量装置、声波接收装置以及程序 | |
| JP4558904B2 (ja) | 画像処理装置及び記憶媒体 | |
| JP2711208B2 (ja) | 組織中音速測定方法 | |
| JP2004257877A (ja) | 音源検出方法、音源検出装置及びロボット | |
| JPH08159864A (ja) | 音響特性測定装置 | |
| CN110584710B (zh) | 一种血流成像方法及装置、存储介质 | |
| JP3977559B2 (ja) | 伝播時間測定方法、環境状態測定方法及び装置 | |
| EP4014857A1 (en) | Probe-cavity motion modeling | |
| CN108760899B (zh) | 一种超声换能器模型制作方法和装置 | |
| JP5298388B2 (ja) | 超音波を用いた温度測定方法及び温度測定装置 | |
| JP3006270B2 (ja) | 超音波位置測定装置 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20040809 |
|
| TRDD | Decision of grant or rejection written | ||
| A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20040921 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20040927 |
|
| R150 | Certificate of patent or registration of utility model |
Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20081001 Year of fee payment: 4 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20091001 Year of fee payment: 5 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20101001 Year of fee payment: 6 |
|
| LAPS | Cancellation because of no payment of annual fees |