JPH07140031A - ガス漏洩源領域、ガス漏洩源位置及びガス漏洩量の推定方法 - Google Patents
ガス漏洩源領域、ガス漏洩源位置及びガス漏洩量の推定方法Info
- Publication number
- JPH07140031A JPH07140031A JP31426693A JP31426693A JPH07140031A JP H07140031 A JPH07140031 A JP H07140031A JP 31426693 A JP31426693 A JP 31426693A JP 31426693 A JP31426693 A JP 31426693A JP H07140031 A JPH07140031 A JP H07140031A
- Authority
- JP
- Japan
- Prior art keywords
- gas
- gas leakage
- point
- detectors
- leak
- 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.)
- Pending
Links
Landscapes
- Examining Or Testing Airtightness (AREA)
Abstract
(57)【要約】
【目的】 ３つのガス検知器からのガス濃度データだ
け、または３つのガス検知器からのガス濃度データと該
ガス検知器の近傍に設置された風速計からの風速データ
のみからガス漏洩源位置とガス漏洩量とを推定する。 【構成】 仮定された位置のガス漏洩源とガス検知器と
の間の位置関係を定め、ガス検知器に於けるガス濃度
と、気象条件によって定まるパラメータと、ガス漏洩点
とガス検知器との間の前記位置関係との間の関係を表現
するガス拡散式から得られるガス漏洩量を風速で除した
値がそれぞれ等しくなる仮定されたガス漏洩源の軌跡を
求め、該軌跡の曲線をガス漏洩の可能性のある領域とす
る。このような曲線を複数求めれば、ガスの漏洩点及び
漏洩量を推定できる。
け、または３つのガス検知器からのガス濃度データと該
ガス検知器の近傍に設置された風速計からの風速データ
のみからガス漏洩源位置とガス漏洩量とを推定する。 【構成】 仮定された位置のガス漏洩源とガス検知器と
の間の位置関係を定め、ガス検知器に於けるガス濃度
と、気象条件によって定まるパラメータと、ガス漏洩点
とガス検知器との間の前記位置関係との間の関係を表現
するガス拡散式から得られるガス漏洩量を風速で除した
値がそれぞれ等しくなる仮定されたガス漏洩源の軌跡を
求め、該軌跡の曲線をガス漏洩の可能性のある領域とす
る。このような曲線を複数求めれば、ガスの漏洩点及び
漏洩量を推定できる。
Description
【０００１】
【産業上の利用分野】本発明は、ガス漏洩に於けるガス
漏洩領域及びガス漏洩量を、ガス検知器からのガス濃度
データだけ、またはガス濃度データ及び風速計からの風
速データから推定する方法に関し、主として可燃性ガ
ス、毒性ガス、各種油及び有機溶剤の蒸気等のガス漏洩
の可能性がある各種屋外プラントに於けるガス漏洩の早
期対処に有用なガス漏洩領域及び漏洩量の推定方法に関
する。
漏洩領域及びガス漏洩量を、ガス検知器からのガス濃度
データだけ、またはガス濃度データ及び風速計からの風
速データから推定する方法に関し、主として可燃性ガ
ス、毒性ガス、各種油及び有機溶剤の蒸気等のガス漏洩
の可能性がある各種屋外プラントに於けるガス漏洩の早
期対処に有用なガス漏洩領域及び漏洩量の推定方法に関
する。
【０００２】
【従来の技術】従来、ガス漏洩領域の推定のための時間
が掛かりすぎ、また、推定される漏洩場所が広すぎた
り、誤差が大きいという問題点があった。また、迅速
に、かつ誤差の少ないガス漏洩領域やガス漏洩量の推定
には風向風速計を必要とした。この場合、よく知られて
いるガス拡散を表す坂上の式やサットンの式を用いてこ
れらの推定を行っていた。
が掛かりすぎ、また、推定される漏洩場所が広すぎた
り、誤差が大きいという問題点があった。また、迅速
に、かつ誤差の少ないガス漏洩領域やガス漏洩量の推定
には風向風速計を必要とした。この場合、よく知られて
いるガス拡散を表す坂上の式やサットンの式を用いてこ
れらの推定を行っていた。
【０００３】以下の説明は、漏洩ガスの拡散による漏洩
ガス濃度と、風速、ガス漏洩量、位置関係、気象条件に
よって定まるパラメータとの関係を表現する式について
はサットンの式（通商産業省立地公害局監修「新訂 公
害防止の技術と法規（大気編）」」４版、昭和５６年８
月発行）をベースにする。もちろん、サットンの式の代
わりに坂上の式または他の理論式、実験式等を用いても
よいし、場所によってこれらを使い分けてもよい。
ガス濃度と、風速、ガス漏洩量、位置関係、気象条件に
よって定まるパラメータとの関係を表現する式について
はサットンの式（通商産業省立地公害局監修「新訂 公
害防止の技術と法規（大気編）」」４版、昭和５６年８
月発行）をベースにする。もちろん、サットンの式の代
わりに坂上の式または他の理論式、実験式等を用いても
よいし、場所によってこれらを使い分けてもよい。
【０００４】サットンの拡散式は漏洩ガスが濃度１００
％のときの任意の３次元位置の漏洩ガス濃度を次式で表
している。
％のときの任意の３次元位置の漏洩ガス濃度を次式で表
している。
【０００５】
【数１】
【０００６】ガス漏洩点からガス検知器までの上部空間
に、架台等の上下方向の拡散を妨げる水平障害物がない
場合、
に、架台等の上下方向の拡散を妨げる水平障害物がない
場合、
【０００７】
【数２】
【０００８】となり、ガス漏洩点からガス検知器までの
上部空間に、架台等の上下方向の拡散を妨げる水平障害
物がある場合で、上下の水平障害物（地面を含む）によ
る漏洩ガスの反復反射がある場合、
上部空間に、架台等の上下方向の拡散を妨げる水平障害
物がある場合で、上下の水平障害物（地面を含む）によ
る漏洩ガスの反復反射がある場合、
【０００９】
【数３】
【００１０】となる。
【００１１】上記の式に含まれる記号をプラントを考慮
して説明すると、 Ｃ（ｘ，ｙ，ｚ）：ガス漏洩点の座標を（０，０，Ｈ）
とし、風向きの方向にｘ軸をとったときの任意の位置座
標（ｘ，ｙ，ｚ）に於けるガス濃度。 Ｑ（０，０，Ｈ）：単位時間当たりの漏洩ガス量。 Ｕ：風速。 Ｈ：漏洩点の地表からの高さである。架台等の床がある
場合は床面からの高さとする。ｚも同じ床面からの任意
の高さとなる。 Ｔ：ガス漏洩点からガス検知器までの上部空間に架台等
の上下方向の拡散を妨げる水平障害物がある場合、上下
の水平障害物（地面を含む）の間の距離。 ｍ：反射の回数で、通常ｍ＝−２〜２をとれば充分に良
い近似になる。 Ｃy，Ｃz，ｎ：Ｈと気象条件によって定まるサットンの
拡散パラメータ。上下に水平障害物（地面を含む）があ
る場合、この間の大気の温度勾配は小さいので中立層を
なすと考えられる。ｎは１／５〜１／２の数値をとり、
中立層では１／４である。また、高さ１０ｍ以下の中立
層の場合、Ｃy＝０．２１（ｍ1/8）、Ｃ z＝０．１２
（ｍ1/8）である。
して説明すると、 Ｃ（ｘ，ｙ，ｚ）：ガス漏洩点の座標を（０，０，Ｈ）
とし、風向きの方向にｘ軸をとったときの任意の位置座
標（ｘ，ｙ，ｚ）に於けるガス濃度。 Ｑ（０，０，Ｈ）：単位時間当たりの漏洩ガス量。 Ｕ：風速。 Ｈ：漏洩点の地表からの高さである。架台等の床がある
場合は床面からの高さとする。ｚも同じ床面からの任意
の高さとなる。 Ｔ：ガス漏洩点からガス検知器までの上部空間に架台等
の上下方向の拡散を妨げる水平障害物がある場合、上下
の水平障害物（地面を含む）の間の距離。 ｍ：反射の回数で、通常ｍ＝−２〜２をとれば充分に良
い近似になる。 Ｃy，Ｃz，ｎ：Ｈと気象条件によって定まるサットンの
拡散パラメータ。上下に水平障害物（地面を含む）があ
る場合、この間の大気の温度勾配は小さいので中立層を
なすと考えられる。ｎは１／５〜１／２の数値をとり、
中立層では１／４である。また、高さ１０ｍ以下の中立
層の場合、Ｃy＝０．２１（ｍ1/8）、Ｃ z＝０．１２
（ｍ1/8）である。
【００１２】漏洩点の推定に関しても上式から分かるよ
うに、単位時間当たりの漏洩ガス量が一定であるとき、
丁度ガス検知器が漏洩点の風下にあるように風が吹いて
いるならば（上式に於てｙ＝０の時）ガス濃度と風速を
乗じたガスフラックスの値が極大となるので、ガスフラ
ックスの値が極大となるような（風）向きの逆方向が漏
洩点の方向を示すこととなる。
うに、単位時間当たりの漏洩ガス量が一定であるとき、
丁度ガス検知器が漏洩点の風下にあるように風が吹いて
いるならば（上式に於てｙ＝０の時）ガス濃度と風速を
乗じたガスフラックスの値が極大となるので、ガスフラ
ックスの値が極大となるような（風）向きの逆方向が漏
洩点の方向を示すこととなる。
【００１３】本発明者は、このことを利用して一定地域
内にある少なくとも２点のガスフラックスの値の時系列
上の極大値の向きの逆方向に線を引き交点を求め、この
点を囲む領域をガス漏洩源の位置として推定する方法を
特許出願した（特願平５−２０４５０３号）。また、一
定地域内にある少なくとも２つのガス検出器の位置に於
けるガスフラックスの値の時系列上の極大値の向きの逆
方向に線を引き交点を求め、この点を利用して丁度ガス
漏洩源の方向から風が吹かない場合でも、より正確なガ
ス漏洩点と、この場合のガス漏洩量を推定する方法を見
出し特許出願した（特願平５−２５０１３１号）。更
に、１つのガス検出器の位置に於ける２つのガスフラッ
クスの値から、ガス漏洩領域とガス漏洩量を推定する方
法を特許出願した（平成５年１０月２７日提出の特許
願）。しかし、これらはいずれも風向風速データを必要
とした。
内にある少なくとも２点のガスフラックスの値の時系列
上の極大値の向きの逆方向に線を引き交点を求め、この
点を囲む領域をガス漏洩源の位置として推定する方法を
特許出願した（特願平５−２０４５０３号）。また、一
定地域内にある少なくとも２つのガス検出器の位置に於
けるガスフラックスの値の時系列上の極大値の向きの逆
方向に線を引き交点を求め、この点を利用して丁度ガス
漏洩源の方向から風が吹かない場合でも、より正確なガ
ス漏洩点と、この場合のガス漏洩量を推定する方法を見
出し特許出願した（特願平５−２５０１３１号）。更
に、１つのガス検出器の位置に於ける２つのガスフラッ
クスの値から、ガス漏洩領域とガス漏洩量を推定する方
法を特許出願した（平成５年１０月２７日提出の特許
願）。しかし、これらはいずれも風向風速データを必要
とした。
【００１４】
【発明が解決しようとする課題】本発明は上記したよう
な従来技術の問題点に鑑みなされたものであり、本発明
の主な目的は、３つのガス検知器からのガス濃度データ
だけ、または３つのガス検知器からのガス濃度データと
該ガス検知器の近傍に設置された風速計からの風速デー
タからサットンの拡散式等を活用してガス漏洩源位置と
ガス漏洩量とを推定するための方法を提供することにあ
る。
な従来技術の問題点に鑑みなされたものであり、本発明
の主な目的は、３つのガス検知器からのガス濃度データ
だけ、または３つのガス検知器からのガス濃度データと
該ガス検知器の近傍に設置された風速計からの風速デー
タからサットンの拡散式等を活用してガス漏洩源位置と
ガス漏洩量とを推定するための方法を提供することにあ
る。
【００１５】
【課題を解決するための手段】上記した課題は、一定地
域内にある３つのガス検知器を選定し、これらから得ら
れた略同一時刻の第１のガス濃度データからガス漏洩領
域を推定するための方法であって、前記ガス検知器の高
さ、またはその近傍の高さをもつ水平面上の点にガス漏
洩源があると仮定し、更に、仮定された前記ガス漏洩源
から、仮定した風向の方位上に引いた線にそれぞれ前記
ガス検知器から垂線を引いて前記ガス漏洩源と前記ガス
検知器との間の位置関係を定め、前記ガス検知器に於け
るガス濃度と、気象条件によって定まるパラメータと、
ガス漏洩点とガス検知器との間の前記位置関係との間の
関係を表現するガス拡散式から得られるガス漏洩量を風
速で除した値がそれぞれ等しくなる点の軌跡を求め、該
軌跡の曲線をガス漏洩の可能性のある領域とすることを
特徴とするガス漏洩領域の推定方法、或いは、上記方法
の結果及び一定地域内にある３つのガス検知器を選定
し、これらから得られた第２のガス濃度データからガス
漏洩源位置及びガス漏洩量を推定するための方法であっ
て、前記で選定したのと同じ３つのガス検知器からの略
同一な別の時刻の第２のガス濃度データ、または前記で
選定したのと同じ２つのガス検知器までを含んでもよい
任意の３つのガス検知器からの任意の略同一な時刻の第
２のガス濃度データに基づいて、前記で仮定したガス漏
洩源の高さの水平面上に於て、前記と同様の演算を行い
新たな軌跡の曲線を求め、前記の曲線との交点をガス漏
洩源位置とし、このガス漏洩源位置に於けるガス漏洩量
を求めることを特徴とするガス漏洩源位置及びガス漏洩
量の推定方法を提供することにより達成される。
域内にある３つのガス検知器を選定し、これらから得ら
れた略同一時刻の第１のガス濃度データからガス漏洩領
域を推定するための方法であって、前記ガス検知器の高
さ、またはその近傍の高さをもつ水平面上の点にガス漏
洩源があると仮定し、更に、仮定された前記ガス漏洩源
から、仮定した風向の方位上に引いた線にそれぞれ前記
ガス検知器から垂線を引いて前記ガス漏洩源と前記ガス
検知器との間の位置関係を定め、前記ガス検知器に於け
るガス濃度と、気象条件によって定まるパラメータと、
ガス漏洩点とガス検知器との間の前記位置関係との間の
関係を表現するガス拡散式から得られるガス漏洩量を風
速で除した値がそれぞれ等しくなる点の軌跡を求め、該
軌跡の曲線をガス漏洩の可能性のある領域とすることを
特徴とするガス漏洩領域の推定方法、或いは、上記方法
の結果及び一定地域内にある３つのガス検知器を選定
し、これらから得られた第２のガス濃度データからガス
漏洩源位置及びガス漏洩量を推定するための方法であっ
て、前記で選定したのと同じ３つのガス検知器からの略
同一な別の時刻の第２のガス濃度データ、または前記で
選定したのと同じ２つのガス検知器までを含んでもよい
任意の３つのガス検知器からの任意の略同一な時刻の第
２のガス濃度データに基づいて、前記で仮定したガス漏
洩源の高さの水平面上に於て、前記と同様の演算を行い
新たな軌跡の曲線を求め、前記の曲線との交点をガス漏
洩源位置とし、このガス漏洩源位置に於けるガス漏洩量
を求めることを特徴とするガス漏洩源位置及びガス漏洩
量の推定方法を提供することにより達成される。
【００１６】前記曲線の交点が２つある場合には、前記
で選定したのと同じ３つのガス検知器から得られた略同
一ではあるが更に別の時刻の第３のガス濃度データ、ま
たは前記で選定したのと同じ２つのガス検知器までを含
んでもよい任意の３つのガス検知器から得られた任意の
略同一な時刻の第３のガス濃度データに基づいて、前記
で仮定したガス漏洩源の高さの水平面上の点に於て、前
記と同様の演算を行いまた新たな軌跡の曲線を求め、こ
れらの曲線と前記の２つの交点と交わる方の交点をガス
漏洩源位置とし、このガス漏洩源位置に於けるガス漏洩
量を求めることができる。
で選定したのと同じ３つのガス検知器から得られた略同
一ではあるが更に別の時刻の第３のガス濃度データ、ま
たは前記で選定したのと同じ２つのガス検知器までを含
んでもよい任意の３つのガス検知器から得られた任意の
略同一な時刻の第３のガス濃度データに基づいて、前記
で仮定したガス漏洩源の高さの水平面上の点に於て、前
記と同様の演算を行いまた新たな軌跡の曲線を求め、こ
れらの曲線と前記の２つの交点と交わる方の交点をガス
漏洩源位置とし、このガス漏洩源位置に於けるガス漏洩
量を求めることができる。
【００１７】また、風速データのある場合に、前記第１
から第３までのガス濃度データに対応した３つのガス漏
洩量が一致すればガス漏洩源の高さは仮定した高さと同
じとし、一致しない場合はガス漏洩源の高さを仮定し直
し、前記の方法で得られたガス漏洩源の水平面位置付近
に於て新たなガス漏洩源位置を求め、このとき得られる
３つのガス漏洩量が一致する３次元座標を、推定される
ガス漏洩源位置とすることができる。
から第３までのガス濃度データに対応した３つのガス漏
洩量が一致すればガス漏洩源の高さは仮定した高さと同
じとし、一致しない場合はガス漏洩源の高さを仮定し直
し、前記の方法で得られたガス漏洩源の水平面位置付近
に於て新たなガス漏洩源位置を求め、このとき得られる
３つのガス漏洩量が一致する３次元座標を、推定される
ガス漏洩源位置とすることができる。
【００１８】また、前記したガス漏洩量を風速で除した
値がそれぞれ等しくなる点の軌跡を求める過程が、ガス
漏洩源があると仮定した水平面上の各点に於て、風向の
方位を変えて１つのガス検知器のガス漏洩量を風速で除
した値と他の２つのガス検知器のガス漏洩量を風速で除
した値との比の対数の絶対値の和が最小となる点を求
め、次に、仮定した点を移動して該和が限りなくゼロに
収斂する点を求め、これらのゼロに収斂する点を繋ぐこ
とによって該軌跡とし、該軌跡上の各点の風向をそれぞ
れの点がガス漏洩点であるときの推定される風向として
求めることからなるのが好ましい。
値がそれぞれ等しくなる点の軌跡を求める過程が、ガス
漏洩源があると仮定した水平面上の各点に於て、風向の
方位を変えて１つのガス検知器のガス漏洩量を風速で除
した値と他の２つのガス検知器のガス漏洩量を風速で除
した値との比の対数の絶対値の和が最小となる点を求
め、次に、仮定した点を移動して該和が限りなくゼロに
収斂する点を求め、これらのゼロに収斂する点を繋ぐこ
とによって該軌跡とし、該軌跡上の各点の風向をそれぞ
れの点がガス漏洩点であるときの推定される風向として
求めることからなるのが好ましい。
【００１９】
【作用】まず（１）式をＱ／Ｕについて解いた形に変形
すると
すると
【００２０】
【数４】
【００２１】となる。
【００２２】ここで、図１に示すようにｉ番目のガス検
知器Ｇi（Ｘi、Ｙi、Ｚi）に関するＮ番目のデータとし
て漏洩ガス濃度Ｃi,Nが得られたこととする。３つのガ
ス検知器に対しＮ番目の風速と風向の方位をそれぞれＵ
N、αNとする。図１に於いてはガス漏洩点Ｐ（Ｘ，Ｙ，
Ｈ）をガス検知器Ｇiを含む水平面またはその近傍の任
意の高さに投影した点Ｐ′（Ｘ，Ｙ，Ｚi）で表してあ
る。図１から分かるように点Ｐ′から任意の風向の方位
に対し線を引き各ガス検知器からこの線に垂線を引い
て、（１）式に示されるＸに相当するＬと、Ｙに相当す
るＷをもとめている。
知器Ｇi（Ｘi、Ｙi、Ｚi）に関するＮ番目のデータとし
て漏洩ガス濃度Ｃi,Nが得られたこととする。３つのガ
ス検知器に対しＮ番目の風速と風向の方位をそれぞれＵ
N、αNとする。図１に於いてはガス漏洩点Ｐ（Ｘ，Ｙ，
Ｈ）をガス検知器Ｇiを含む水平面またはその近傍の任
意の高さに投影した点Ｐ′（Ｘ，Ｙ，Ｚi）で表してあ
る。図１から分かるように点Ｐ′から任意の風向の方位
に対し線を引き各ガス検知器からこの線に垂線を引い
て、（１）式に示されるＸに相当するＬと、Ｙに相当す
るＷをもとめている。
【００２３】（４）式に於いてはｉ番目のガス検知器Ｇ
iに関するＮ番目のデータに対応して距離ＬはＬi,N、距
離ＷはＷi,Nと表し、ガス検知器Ｇiの高さをＺi、ガス
漏洩点Ｐの床面からの高さはＨとして表し、この時
（４）式によって得られるガス漏洩点をＱi,Nとしてい
る。今、ｉ，ｊ，ｋ番目３つのガス検知器Ｇi、Ｇj、Ｇ
kの漏洩ガス濃度が測定可能下限値以上となり、Ｃi,N、
Ｃj,N、Ｃk,Nの３つのデータが得られたこととする。こ
の場合ガス漏洩点Ｐ（Ｘ，Ｙ，Ｈ）の高さが分からない
ので、ガス検知器Ｇiを含む水平面にまたはその近傍の
任意の高さにあると仮定する。Ｎ番目のデータから
（４）式によって得られるＱi,N／ＵN、Ｑj,N／ＵN、Ｑ
k,N／ＵNはそれぞれ等しくならなければならないから、 （Ｑi,N／ＵN）／（Ｑj,N／ＵN）＝１ （５） （Ｑi,N／ＵN）／（Ｑk,N／ＵN）＝１ （６） となる 。
iに関するＮ番目のデータに対応して距離ＬはＬi,N、距
離ＷはＷi,Nと表し、ガス検知器Ｇiの高さをＺi、ガス
漏洩点Ｐの床面からの高さはＨとして表し、この時
（４）式によって得られるガス漏洩点をＱi,Nとしてい
る。今、ｉ，ｊ，ｋ番目３つのガス検知器Ｇi、Ｇj、Ｇ
kの漏洩ガス濃度が測定可能下限値以上となり、Ｃi,N、
Ｃj,N、Ｃk,Nの３つのデータが得られたこととする。こ
の場合ガス漏洩点Ｐ（Ｘ，Ｙ，Ｈ）の高さが分からない
ので、ガス検知器Ｇiを含む水平面にまたはその近傍の
任意の高さにあると仮定する。Ｎ番目のデータから
（４）式によって得られるＱi,N／ＵN、Ｑj,N／ＵN、Ｑ
k,N／ＵNはそれぞれ等しくならなければならないから、 （Ｑi,N／ＵN）／（Ｑj,N／ＵN）＝１ （５） （Ｑi,N／ＵN）／（Ｑk,N／ＵN）＝１ （６） となる 。
【００２４】（５）及び（６）式が同時に成立するため
の条件は
の条件は
【００２５】
【数５】
【００２６】のようにＡijkを定義した時にＡijk＝０と
置いて得られる。
置いて得られる。
【００２７】この場合、上記のように対数の絶対値の和
をとる代わりに、対数の２乗の和をとってもよい。
をとる代わりに、対数の２乗の和をとってもよい。
【００２８】この場合、仮定した点Ｐから風向の方位に
引いた線に３つのガス検知器から垂線を引いて得られる
３つの交点のうち外側の２つの交点で挟まれた間の点
は、サットンの式の適用できない範囲であることに注意
する必要がある。
引いた線に３つのガス検知器から垂線を引いて得られる
３つの交点のうち外側の２つの交点で挟まれた間の点
は、サットンの式の適用できない範囲であることに注意
する必要がある。
【００２９】ガス漏洩点があると仮定した高さの水平面
の任意の点で風向の方位を変えてＡijkが最小となる値
を求め、この値を高さとする２次曲面を鳥瞰すると、複
数の尾根と谷があり、この谷のところにＡijk＝０とな
る軌跡がある。この軌跡を求めるためには、対象とする
地域を例えば５ｍ間隔の格子に分割し、その格子の交点
の位置でＡijkを求めると、大まかな山や谷が見えるの
で、この谷の部分を近似的に求める軌跡としてもよい。
の任意の点で風向の方位を変えてＡijkが最小となる値
を求め、この値を高さとする２次曲面を鳥瞰すると、複
数の尾根と谷があり、この谷のところにＡijk＝０とな
る軌跡がある。この軌跡を求めるためには、対象とする
地域を例えば５ｍ間隔の格子に分割し、その格子の交点
の位置でＡijkを求めると、大まかな山や谷が見えるの
で、この谷の部分を近似的に求める軌跡としてもよい。
【００３０】更に精密に軌跡の位置を求めようとするな
らば、例えば谷の部分に於てＸを一定値に保ったままＹ
を変えてＡijkが限りなく０に収束する点を求め、各Ｘ
の値についてこれらの点を結べば正確な軌跡を求めるこ
とができる。この場合ＸとＹを入れ換えてもよい。
らば、例えば谷の部分に於てＸを一定値に保ったままＹ
を変えてＡijkが限りなく０に収束する点を求め、各Ｘ
の値についてこれらの点を結べば正確な軌跡を求めるこ
とができる。この場合ＸとＹを入れ換えてもよい。
【００３１】また、各ガス検知器を中心とする極座標に
於て動径一定の円周を等間隔に分割してＡijkを求め、
この円周上でＡijkが極小値をとる位置付近を更に細か
く分割してＡijkが限りなく０に収束する点を求めても
よい。この場合、Ａijk＝０となる軌跡は各ガス検知器
の極めて近い近傍から放射状に伸びる曲線であるので、
動径を増加させて前の方向の延長部分を中心にＡijkが
限りなく０に収束する点を求めることによって容易に軌
跡が得られる。このようにして得られる軌跡の曲線がガ
ス漏洩の可能性のある領域である。
於て動径一定の円周を等間隔に分割してＡijkを求め、
この円周上でＡijkが極小値をとる位置付近を更に細か
く分割してＡijkが限りなく０に収束する点を求めても
よい。この場合、Ａijk＝０となる軌跡は各ガス検知器
の極めて近い近傍から放射状に伸びる曲線であるので、
動径を増加させて前の方向の延長部分を中心にＡijkが
限りなく０に収束する点を求めることによって容易に軌
跡が得られる。このようにして得られる軌跡の曲線がガ
ス漏洩の可能性のある領域である。
【００３２】Ｎ＝１のデータに基づいた複数の軌跡の曲
線とＮ＝２のデータに基づいた複数の軌跡の曲線とは一
般に２つの交点をもち、この２つの交点のいずれかが真
のガス漏洩源である。次にＮ＝３のデータに基づいた複
数の軌跡の曲線を求めれば、２つの交点のいずれかの交
点と再び交わる方が真のガス漏洩源である。Ｎ＝１のデ
ータ及びＮ＝２のデータに基づいた２つの交点は３つの
ガス検知器を中央に配した対称的な位置にあり、現場で
はその風向きからいずれが真のガス漏洩源かを直ちに判
断できよう。また、２つの交点のうちの片方が対象とす
る地域から外れていることもある。このようにして２つ
の交点が算出された段階でそのうちのいずれが真のガス
漏洩源かを状況により判断することもできる。
線とＮ＝２のデータに基づいた複数の軌跡の曲線とは一
般に２つの交点をもち、この２つの交点のいずれかが真
のガス漏洩源である。次にＮ＝３のデータに基づいた複
数の軌跡の曲線を求めれば、２つの交点のいずれかの交
点と再び交わる方が真のガス漏洩源である。Ｎ＝１のデ
ータ及びＮ＝２のデータに基づいた２つの交点は３つの
ガス検知器を中央に配した対称的な位置にあり、現場で
はその風向きからいずれが真のガス漏洩源かを直ちに判
断できよう。また、２つの交点のうちの片方が対象とす
る地域から外れていることもある。このようにして２つ
の交点が算出された段階でそのうちのいずれが真のガス
漏洩源かを状況により判断することもできる。
【００３３】風速データのある場合にＮ＝１からＮ＝３
までのガス濃度データに対応したガス漏洩量を風速で除
した値にそれぞれ風速を乗じて得られる３つのガス漏洩
量を算出し、これらが一致すればガス漏洩源の高さは仮
定した高さに等しい。
までのガス濃度データに対応したガス漏洩量を風速で除
した値にそれぞれ風速を乗じて得られる３つのガス漏洩
量を算出し、これらが一致すればガス漏洩源の高さは仮
定した高さに等しい。
【００３４】一致しない場合はガス漏洩源の高さを仮定
し直し、前記の方法で既に得られたガス漏洩源の水平面
位置付近に於て新たなガス漏洩源位置を求め、この時得
られる３つのガス漏洩量が一致する３次元座標がガス漏
洩源位置である。
し直し、前記の方法で既に得られたガス漏洩源の水平面
位置付近に於て新たなガス漏洩源位置を求め、この時得
られる３つのガス漏洩量が一致する３次元座標がガス漏
洩源位置である。
【００３５】漏洩ガス濃度は風向きの変化、風速の変
化、大気の乱流による揺らぎ等により、時間と共に変化
する。特に風向きが変化するときにはどのガス検知器の
ガス濃度も変化が大きいので、データ処理によりガス検
知器の位置による検出時刻の差を考慮して、ガス濃度が
大きく変化した後の比較的安定な時期の漏洩ガス濃度を
平均して計算に使用するデータとするとよい。
化、大気の乱流による揺らぎ等により、時間と共に変化
する。特に風向きが変化するときにはどのガス検知器の
ガス濃度も変化が大きいので、データ処理によりガス検
知器の位置による検出時刻の差を考慮して、ガス濃度が
大きく変化した後の比較的安定な時期の漏洩ガス濃度を
平均して計算に使用するデータとするとよい。
【００３６】風速データのない場合には、風速を強風、
弱風、微風のように風速を幅をもって捉え、ガス漏洩量
を風速で除した値にこれを乗じれば、幅はあるが大凡の
ガス漏洩量とガス漏洩源の高さが把握できる。
弱風、微風のように風速を幅をもって捉え、ガス漏洩量
を風速で除した値にこれを乗じれば、幅はあるが大凡の
ガス漏洩量とガス漏洩源の高さが把握できる。
【００３７】このようにガス漏洩点及びガス漏洩量が分
かれば任意の風向風速に対して漏洩ガス濃度分布などを
計算でき安全対策に役立てることができる。
かれば任意の風向風速に対して漏洩ガス濃度分布などを
計算でき安全対策に役立てることができる。
【００３８】上記においてもし気象条件が不明の場合若
しくは確定できない場合は気象条件を変数としてサット
ンの拡散パラメーターを変えてガス漏洩領域、ガス漏洩
点及びガス漏洩量を求めてもよい。
しくは確定できない場合は気象条件を変数としてサット
ンの拡散パラメーターを変えてガス漏洩領域、ガス漏洩
点及びガス漏洩量を求めてもよい。
【００３９】実際のプラントでは機器等による気流の乱
れがあり、ここに示したような理想的な状態でのガス拡
散からはズレるが、先ず理想的な状態でのガス漏洩点及
びガス漏洩量を推定し、次に機器の周りの流線を想定す
る等ズレを補正することにより実際の状態に近付けるこ
とができる。
れがあり、ここに示したような理想的な状態でのガス拡
散からはズレるが、先ず理想的な状態でのガス漏洩点及
びガス漏洩量を推定し、次に機器の周りの流線を想定す
る等ズレを補正することにより実際の状態に近付けるこ
とができる。
【００４０】本発明の方法は固定式ガス検知器だけでな
く、ガス検知器間の距離を短縮して移動式またはポータ
ブルとして適用しても良い。
く、ガス検知器間の距離を短縮して移動式またはポータ
ブルとして適用しても良い。
【００４１】
【実施例】本発明の好適実施例を図２〜図５に基づいて
詳しく説明する。
詳しく説明する。
【００４２】本発明のガス漏洩点とガス漏洩量の推定方
法を実現するシステム例は図２に示すようにＭ（≧３）
個のガス検知器Ｇ（１）〜Ｇ（Ｍ）と、Ｎ（≧１）個の
風速計Ｂ（１）〜Ｂ（Ｎ）とを有し、データ処理装置３
に接続されたガス濃度データ収集記憶装置１及び風速デ
ータ収集記憶装置２にガス濃度データ及び風速データを
入力するようになっている。
法を実現するシステム例は図２に示すようにＭ（≧３）
個のガス検知器Ｇ（１）〜Ｇ（Ｍ）と、Ｎ（≧１）個の
風速計Ｂ（１）〜Ｂ（Ｎ）とを有し、データ処理装置３
に接続されたガス濃度データ収集記憶装置１及び風速デ
ータ収集記憶装置２にガス濃度データ及び風速データを
入力するようになっている。
【００４３】また、データ処理装置３にはガス濃度警報
設定値、ガス検知器及び風速計位置データ等を入力した
り、各種操作を行うためのキーボード４、ガス検知器Ｇ
（１）〜Ｇ（Ｍ）及び風速計Ｂ（１）〜Ｂ（Ｎ）の番号
とその位置を記憶する計器位置記憶装置５、ガス検知器
の番号毎のガス濃度、風速等を必要な時間間隔、例えば
１０秒毎に過去の一定時間例えば過去１時間、時系列的
に記憶する外部記憶装置６、警報器７、ＣＲＴ表示装置
８及びプリンタ９が接続されている。
設定値、ガス検知器及び風速計位置データ等を入力した
り、各種操作を行うためのキーボード４、ガス検知器Ｇ
（１）〜Ｇ（Ｍ）及び風速計Ｂ（１）〜Ｂ（Ｎ）の番号
とその位置を記憶する計器位置記憶装置５、ガス検知器
の番号毎のガス濃度、風速等を必要な時間間隔、例えば
１０秒毎に過去の一定時間例えば過去１時間、時系列的
に記憶する外部記憶装置６、警報器７、ＣＲＴ表示装置
８及びプリンタ９が接続されている。
【００４４】各ガス検知器Ｇ（１）〜Ｇ（Ｍ）に対し
て、それぞれの位置での風速の近似値を与える風速計と
してその近傍にある風速計Ｂ（１）〜Ｂ（Ｎ）が選択さ
れ、各ガス検知器からのガス濃度データと風速データと
が１対１に対応するようになっている。
て、それぞれの位置での風速の近似値を与える風速計と
してその近傍にある風速計Ｂ（１）〜Ｂ（Ｎ）が選択さ
れ、各ガス検知器からのガス濃度データと風速データと
が１対１に対応するようになっている。
【００４５】この場合１個のガス検知器に対し複数の風
速計を選びそれらの風速データをその位置関係により内
挿または外挿して用いてもよい。通常ＮはＭの１００分
の１以上、望ましくは２０分の１以上でＭを超えない数
である。また、ガス検知器と風速計とを一体化して同一
場所の空気の風速とガス濃度とを同時に計れるようにし
てもよく、この場合これらの個数Ｍ，ＮはＭ＝Ｎとな
る。
速計を選びそれらの風速データをその位置関係により内
挿または外挿して用いてもよい。通常ＮはＭの１００分
の１以上、望ましくは２０分の１以上でＭを超えない数
である。また、ガス検知器と風速計とを一体化して同一
場所の空気の風速とガス濃度とを同時に計れるようにし
てもよく、この場合これらの個数Ｍ，ＮはＭ＝Ｎとな
る。
【００４６】本システムに於いては常時送信されるガス
検知器と風速計の測定データはデジタル信号として変換
され、データ処理装置３にてデータの選択、記憶、演算
等がなされるようになっている。
検知器と風速計の測定データはデジタル信号として変換
され、データ処理装置３にてデータの選択、記憶、演算
等がなされるようになっている。
【００４７】ガス濃度データ収集記憶装置１は、ガス検
知器からのガス濃度データを過去のデータを更新しなが
ら一定時間間隔で時系列的に収集記憶するようになって
いる。同様に風速データ収集記憶装置２は、風速計から
の風速データを更新しながら収集記憶するようになって
いる。
知器からのガス濃度データを過去のデータを更新しなが
ら一定時間間隔で時系列的に収集記憶するようになって
いる。同様に風速データ収集記憶装置２は、風速計から
の風速データを更新しながら収集記憶するようになって
いる。
【００４８】データ処理装置３は、ガス濃度データをガ
ス濃度データ収集記憶装置１から、また、各ガス検知器
に対応する風速データを風速データ収集記憶装置２から
取り出し、データの時間的変化が大きいデータは捨て、
比較的に安定なデータの平均値を記憶するようにデータ
を処理できるようになっている。
ス濃度データ収集記憶装置１から、また、各ガス検知器
に対応する風速データを風速データ収集記憶装置２から
取り出し、データの時間的変化が大きいデータは捨て、
比較的に安定なデータの平均値を記憶するようにデータ
を処理できるようになっている。
【００４９】図３及び図４は本発明のＣＲＴ表示装置８
による表示例であり、プラントの中の一部である長方形
をなす架台部分を示し、架台に設置されている機器の表
示はここで省略している。各ガス検知器Ｇ（１）〜Ｇ
（１０）は２０ｍ間隔で架台外周部に配置されている。
この例では、やや遠い位置にあるが計器室の屋上にある
風向風速計からの風速データを用いるようになってい
る。
による表示例であり、プラントの中の一部である長方形
をなす架台部分を示し、架台に設置されている機器の表
示はここで省略している。各ガス検知器Ｇ（１）〜Ｇ
（１０）は２０ｍ間隔で架台外周部に配置されている。
この例では、やや遠い位置にあるが計器室の屋上にある
風向風速計からの風速データを用いるようになってい
る。
【００５０】ガス濃度が警報設定値を超えた場合には警
報を発し、対応するガス検知器に対する表示を例えばオ
レンジ色に変えて点滅させるようになっている。以下に
数値例によって示すと、図３及び図４に示す地域は高さ
５ｍの架台とその上４ｍの位置にある架台に挟まれた空
間であり、濃度１００％のエチレンが漏洩し、最初にガ
ス検知器Ｇ（１）、Ｇ（２）、Ｇ（４）がガス濃度を検
知した。
報を発し、対応するガス検知器に対する表示を例えばオ
レンジ色に変えて点滅させるようになっている。以下に
数値例によって示すと、図３及び図４に示す地域は高さ
５ｍの架台とその上４ｍの位置にある架台に挟まれた空
間であり、濃度１００％のエチレンが漏洩し、最初にガ
ス検知器Ｇ（１）、Ｇ（２）、Ｇ（４）がガス濃度を検
知した。
【００５１】ガス検知器の床からの高さは、０．５ｍで
あり、最初はガス漏洩点の高さが不明なのでとりあえず
ガス検知器の高さと同じとする。即ち、Ｔ＝４（ｍ），
Ｚ1〜Ｚ10＝０．５（ｍ），Ｈ＝０．５（ｍ），Ｃy＝
０．２１，Ｃz＝０．１２，ｎ＝１／４とする。それぞ
れ最初の２分間の間にほぼ安定した３〜６個のガス濃度
データが得られ、これらの平均値として以下のデータが
得られた。即ち Ｃ1,1＝１０００×１０-6（ｍ3／ｍ3） Ｃ2,1＝５００×１０-6（ｍ3／ｍ3） Ｃ4,1＝２００×１０-6（ｍ3／ｍ3）
あり、最初はガス漏洩点の高さが不明なのでとりあえず
ガス検知器の高さと同じとする。即ち、Ｔ＝４（ｍ），
Ｚ1〜Ｚ10＝０．５（ｍ），Ｈ＝０．５（ｍ），Ｃy＝
０．２１，Ｃz＝０．１２，ｎ＝１／４とする。それぞ
れ最初の２分間の間にほぼ安定した３〜６個のガス濃度
データが得られ、これらの平均値として以下のデータが
得られた。即ち Ｃ1,1＝１０００×１０-6（ｍ3／ｍ3） Ｃ2,1＝５００×１０-6（ｍ3／ｍ3） Ｃ4,1＝２００×１０-6（ｍ3／ｍ3）
【００５２】これらの３つのガス濃度のデータから作用
の項で説明した方法で求めた対象地域のガス漏洩領域を
図３に曲線ｑ２１１、ｑ４１１及びＧ（１）、Ｇ
（２）、Ｇ（４）から伸びる短い実線で示している。対
象地域外についても参考に示せば図５の実線で示すｑ１
１１、ｑ１１２、ｑ２１１、ｑ２１２、ｑ４１１、ｑ４
１２の６本の曲線が得られた。
の項で説明した方法で求めた対象地域のガス漏洩領域を
図３に曲線ｑ２１１、ｑ４１１及びＧ（１）、Ｇ
（２）、Ｇ（４）から伸びる短い実線で示している。対
象地域外についても参考に示せば図５の実線で示すｑ１
１１、ｑ１１２、ｑ２１１、ｑ２１２、ｑ４１１、ｑ４
１２の６本の曲線が得られた。
【００５３】次の２分間の間に、それぞれ３〜６個のほ
ぼ安定したガス濃度データが得られ、これらの平均値と
して以下のデータが得られた。即ち Ｃ1,1＝３１０３×１０-6（ｍ3／ｍ3） Ｃ2,2＝１１２×１０-6（ｍ3／ｍ3） Ｃ4,2＝１０２９×１０-6（ｍ3／ｍ3）
ぼ安定したガス濃度データが得られ、これらの平均値と
して以下のデータが得られた。即ち Ｃ1,1＝３１０３×１０-6（ｍ3／ｍ3） Ｃ2,2＝１１２×１０-6（ｍ3／ｍ3） Ｃ4,2＝１０２９×１０-6（ｍ3／ｍ3）
【００５４】これらの３つのガス濃度のデータから、作
用の項で説明した方法で求めた対象地域のガス漏洩領域
を図３に曲線ｑ２２１、ｑ４２１及びＧ（１）、Ｇ
（２）、Ｇ（４）から伸びる短い点線で示している。
用の項で説明した方法で求めた対象地域のガス漏洩領域
を図３に曲線ｑ２２１、ｑ４２１及びＧ（１）、Ｇ
（２）、Ｇ（４）から伸びる短い点線で示している。
【００５５】対象地域外についても参考に示せば、図５
の点線で示す６本の曲線が得られた。これらの実線と点
線で示される曲線は対象地域内に於て１ヶ所で交差し、
その座標はＧ（１）を原点とし、Ｇ（２）方向をＸ軸と
し、Ｇ（４）方向をＹ軸とした場合（８．０，５０．
０，０．５）となる。この点はガス検知器の高さにある
水平面上にガス漏洩源があると仮定した場合の近似的な
ガス漏洩点位置Ｐ′を示している。
の点線で示す６本の曲線が得られた。これらの実線と点
線で示される曲線は対象地域内に於て１ヶ所で交差し、
その座標はＧ（１）を原点とし、Ｇ（２）方向をＸ軸と
し、Ｇ（４）方向をＹ軸とした場合（８．０，５０．
０，０．５）となる。この点はガス検知器の高さにある
水平面上にガス漏洩源があると仮定した場合の近似的な
ガス漏洩点位置Ｐ′を示している。
【００５６】図３にはこのＰ′点のそれぞれの風向をｑ
線に対応させて実線、点線で中央部に示し、また、ガス
漏洩量を風速で除した値の平均値を上部に示した。この
場合風速データがある場合には風速で除した値に風速を
乗じてその平均値をガス漏洩量としてよい。対象地域外
をも含む地域では図５に示すようにもう１つの交点即ち
ｑ１１１、ｑ１２１の交点があり、対象地域が限定され
ていない場合には２つのガス漏洩の候補点があることに
なる。現場での風向が分かればいずれが真のガス漏洩点
かが判断できる。
線に対応させて実線、点線で中央部に示し、また、ガス
漏洩量を風速で除した値の平均値を上部に示した。この
場合風速データがある場合には風速で除した値に風速を
乗じてその平均値をガス漏洩量としてよい。対象地域外
をも含む地域では図５に示すようにもう１つの交点即ち
ｑ１１１、ｑ１２１の交点があり、対象地域が限定され
ていない場合には２つのガス漏洩の候補点があることに
なる。現場での風向が分かればいずれが真のガス漏洩点
かが判断できる。
【００５７】更に次の２分間の間にガス検知器Ｇ
（１）、Ｇ（４）、Ｇ（６）がガス濃度を検知し、それ
ぞれ３〜６個のほぼ安定したガス濃度データが得られ、
これらの平均値として以下のデータが得られた。即ち Ｃ1,3＝１９９７×１０-6（ｍ3／ｍ3） Ｃ4,3＝１３２５２×１０-6（ｍ3／ｍ3） Ｃ6,3＝１０２×１０-6（ｍ3／ｍ3）
（１）、Ｇ（４）、Ｇ（６）がガス濃度を検知し、それ
ぞれ３〜６個のほぼ安定したガス濃度データが得られ、
これらの平均値として以下のデータが得られた。即ち Ｃ1,3＝１９９７×１０-6（ｍ3／ｍ3） Ｃ4,3＝１３２５２×１０-6（ｍ3／ｍ3） Ｃ6,3＝１０２×１０-6（ｍ3／ｍ3）
【００５８】これらの３つのガス濃度のデータから作用
の項で説明した方法で求めた対象地域のガス漏洩領域を
図４に曲線ｑ６３１及びＧ（１）、Ｇ（６）から伸びる
短い破線で示している。
の項で説明した方法で求めた対象地域のガス漏洩領域を
図４に曲線ｑ６３１及びＧ（１）、Ｇ（６）から伸びる
短い破線で示している。
【００５９】曲線ｑ６３１はｑ４１１とｑ４２１の交点
である近似的なガス漏洩点位置Ｐ′ともほぼ１点で交差
しており、この点が真のガス漏洩点の近似的な座標を与
えることになる。対象地域外についても参考に示せば図
５の破線で示すＹ軸に対称な４本の曲線が得られた。図
５からも分かるように、３つの曲線がほぼ１点に交わる
のはＰ′点のみであり、この点が真のガス漏洩点である
ことは明かである。
である近似的なガス漏洩点位置Ｐ′ともほぼ１点で交差
しており、この点が真のガス漏洩点の近似的な座標を与
えることになる。対象地域外についても参考に示せば図
５の破線で示すＹ軸に対称な４本の曲線が得られた。図
５からも分かるように、３つの曲線がほぼ１点に交わる
のはＰ′点のみであり、この点が真のガス漏洩点である
ことは明かである。
【００６０】風速データがある場合には作用の項で説明
した方法によりガス漏洩点の高さも算出できる。即ち近
似的なガス漏洩点位置Ｐ′の前後左右に各１ｍの広さを
もつ水平面を想定し、０からＴの間の高さを０．５ｍ間
隔で変えて３つの曲線が１点に交わるＰ′点を求め、こ
の点のガス漏洩量を風速で除した値にそれぞれの風速を
乗じて一致する高さを求めた。
した方法によりガス漏洩点の高さも算出できる。即ち近
似的なガス漏洩点位置Ｐ′の前後左右に各１ｍの広さを
もつ水平面を想定し、０からＴの間の高さを０．５ｍ間
隔で変えて３つの曲線が１点に交わるＰ′点を求め、こ
の点のガス漏洩量を風速で除した値にそれぞれの風速を
乗じて一致する高さを求めた。
【００６１】それぞれの風速は Ｕ1＝１．０（ｍ／ｓ） Ｕ2＝０．９（ｍ／ｓ） Ｕ3＝１．１（ｍ／ｓ） であり、Ｈ＝１．０ｍに於て全てのガス漏洩量が等しく
なり、この時の３次元位置座標とガス漏洩量を図４の上
部に、またこの時の風向を中央部に示した。
なり、この時の３次元位置座標とガス漏洩量を図４の上
部に、またこの時の風向を中央部に示した。
【００６２】このようにして全てのガス漏洩量が等しく
ならない時には比較的一致している２つ以上の高さの計
算結果を内挿してより精密な高さを求めればよい。
ならない時には比較的一致している２つ以上の高さの計
算結果を内挿してより精密な高さを求めればよい。
【００６３】
【発明の効果】上記した説明から明らかなように、本発
明によるガス漏洩領域、ガス漏洩源及びガス漏洩量の推
定方法によれば、３つ以上のガス検知器からの濃度デー
タだけを処理してガス漏洩量を風速で除した値が互いに
一致する点をガス漏洩源位置として求め、更に近傍の風
速データがある場合にはガス漏洩量が全て一致する点を
３次元的ガス漏洩源位置として求めることによって、早
期にかつ正確にガス漏洩源位置とガス漏洩量を推定する
ことができる。
明によるガス漏洩領域、ガス漏洩源及びガス漏洩量の推
定方法によれば、３つ以上のガス検知器からの濃度デー
タだけを処理してガス漏洩量を風速で除した値が互いに
一致する点をガス漏洩源位置として求め、更に近傍の風
速データがある場合にはガス漏洩量が全て一致する点を
３次元的ガス漏洩源位置として求めることによって、早
期にかつ正確にガス漏洩源位置とガス漏洩量を推定する
ことができる。
【図１】ガス検知器とガス漏洩点の位置関係を示す図で
ある。
ある。
【図２】本発明が適用されたガス漏洩監視システムの構
成を示すブロック図である。
成を示すブロック図である。
【図３】ガス漏洩領域及びガス漏洩量推定に於けるＣＲ
Ｔ表示装置への表示例を示す。
Ｔ表示装置への表示例を示す。
【図４】ガス漏洩領域及びガス漏洩量推定に於けるＣＲ
Ｔ表示装置への表示例を示す。
Ｔ表示装置への表示例を示す。
【図５】対象とする地域外をも含む広範囲な地域に於け
るガス漏洩の可能性のある領域を示す図である。
るガス漏洩の可能性のある領域を示す図である。
１ ガス濃度データ収集記憶装置 ２ 風向風速データ収集記憶装置 ３ データ処理装置 ４ キーボード ５ 計器位置記憶装置 ６ 外部記憶装置 ７ 警報器 ８ ＣＲＴ表示装置 ９ プリンタ
Claims (5)
- 【請求項１】 一定地域内にある３つのガス検知器を
選定し、これらから得られた略同一時刻の第１のガス濃
度データからガス漏洩領域を推定するための方法であっ
て、 前記ガス検知器の高さ、またはその近傍の高さをもつ水
平面上の点にガス漏洩源があると仮定し、 更に、仮定された前記ガス漏洩源から、仮定した風向の
方位上に引いた線にそれぞれ前記ガス検知器から垂線を
引いて前記ガス漏洩源と前記ガス検知器との間の位置関
係を定め、 前記ガス検知器に於けるガス濃度と、気象条件によって
定まるパラメータと、ガス漏洩点とガス検知器との間の
前記位置関係との間の関係を表現するガス拡散式から得
られるガス漏洩量を風速で除した値がそれぞれ等しくな
る点の軌跡を求め、 該軌跡の曲線をガス漏洩の可能性のある領域とすること
を特徴とするガス漏洩領域の推定方法。 - 【請求項２】 前記請求項１の結果及び一定地域内に
ある３つのガス検知器を選定し、これらから得られた第
２のガス濃度データからガス漏洩源位置及びガス漏洩量
を推定するための方法であって、 前記で選定したのと同じ３つのガス検知器からの略同一
な別の時刻の第２のガス濃度データ、または前記で選定
したのと同じ２つのガス検知器までを含んでもよい任意
の３つのガス検知器からの任意の略同一な時刻の第２の
ガス濃度データに基づいて、 前記で仮定したガス漏洩源の高さの水平面上に於て、前
記と同様の演算を行い新たな軌跡の曲線を求め、前記の
曲線との交点をガス漏洩源位置とし、このガス漏洩源位
置に於けるガス漏洩量を求めることを特徴とするガス漏
洩源位置及びガス漏洩量の推定方法。 - 【請求項３】 前記曲線の交点が２つある場合、前記
で選定したのと同じ３つのガス検知器から得られた略同
一ではあるが更に別の時刻の第３のガス濃度データ、ま
たは前記で選定したのと同じ２つのガス検知器までを含
んでもよい任意の３つのガス検知器から得られた任意の
略同一な時刻の第３のガス濃度データに基づいて、前記
で仮定したガス漏洩源の高さの水平面上の点に於て、前
記と同様の演算を行いまた新たな軌跡の曲線を求め、こ
れらの曲線と前記の２つの交点と交わる方の交点をガス
漏洩源位置とし、このガス漏洩源位置に於けるガス漏洩
量を求めることを特徴とする請求項２に記載のガス漏洩
源位置及びガス漏洩量の推定方法。 - 【請求項４】 風速データのある場合に、前記第１か
ら第３までのガス濃度データに対応した３つのガス漏洩
量が一致すればガス漏洩源の高さは仮定した高さと同じ
とし、一致しない場合はガス漏洩源の高さを仮定し直
し、前記の方法で得られたガス漏洩源の水平面位置付近
に於て新たなガス漏洩源位置を求め、このとき得られる
３つのガス漏洩量が一致する３次元座標を、推定される
ガス漏洩源位置とすることを特徴とする請求項３に記載
のガス漏洩源位置及びガス漏洩量の推定方法。 - 【請求項５】 前記したガス漏洩量を風速で除した値
がそれぞれ等しくなる点の軌跡を求める過程が、 ガス漏洩源があると仮定した水平面上の各点に於て、風
向の方位を変えて１つのガス検知器のガス漏洩量を風速
で除した値と他の２つのガス検知器のガス漏洩量を風速
で除した値との比の対数の絶対値の和が最小となる点を
求め、 次に、仮定した点を移動して該和が限りなくゼロに収斂
する点を求め、 これらのゼロに収斂する点を繋ぐことによって該軌跡と
し、該軌跡上の各点の風向をそれぞれの点がガス漏洩点
であるときの推定される風向として求めることを特徴と
する請求項１乃至４のいずれかに記載の推定方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP31426693A JPH07140031A (ja) | 1993-11-18 | 1993-11-18 | ガス漏洩源領域、ガス漏洩源位置及びガス漏洩量の推定方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP31426693A JPH07140031A (ja) | 1993-11-18 | 1993-11-18 | ガス漏洩源領域、ガス漏洩源位置及びガス漏洩量の推定方法 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| JPH07140031A true JPH07140031A (ja) | 1995-06-02 |
Family
ID=18051292
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP31426693A Pending JPH07140031A (ja) | 1993-11-18 | 1993-11-18 | ガス漏洩源領域、ガス漏洩源位置及びガス漏洩量の推定方法 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JPH07140031A (ja) |
Cited By (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| GB2350432A (en) * | 1999-05-27 | 2000-11-29 | Bg Intellectual Pty Ltd | Apparatus for tracing a gas leak |
| CN114754934A (zh) * | 2021-01-08 | 2022-07-15 | 中国石油天然气集团有限公司 | 气体泄漏检测方法 |
| CN116296120A (zh) * | 2023-03-16 | 2023-06-23 | 中国科学院武汉岩土力学研究所 | 一种面向盐穴储氦库井场的氦气泄漏源定位方法 |
| CN116429328A (zh) * | 2022-12-12 | 2023-07-14 | 宁波永大建设集团有限公司 | 一种外墙窗户渗漏检测方法、系统、存储介质及智能终端 |
| CN117191928A (zh) * | 2023-08-08 | 2023-12-08 | 广州科易光电技术有限公司 | 气体泄漏检测方法、装置、电子设备及存储介质 |
| WO2024155857A1 (en) * | 2023-01-19 | 2024-07-25 | Schlumberger Technology Corporation | Leak emissions sensor systems and processes for installing and using same |
| US12254622B2 (en) | 2023-06-16 | 2025-03-18 | Schlumberger Technology Corporation | Computing emission rate from gas density images |
| US12292310B2 (en) | 2022-12-15 | 2025-05-06 | Schlumberger Technology Corporation | Machine learning based methane emissions monitoring |
| US12480924B2 (en) | 2022-08-03 | 2025-11-25 | Schlumberger Technology Corporation | Automated record quality determination and processing for pollutant emission quantification |
| US12555125B2 (en) | 2023-06-09 | 2026-02-17 | Schlumberger Technology Corporation | Emission detecting camera placement planning using 3D models |
| US12577871B2 (en) | 2022-08-03 | 2026-03-17 | Schlumberger Technology Corporation | Linear cut generation method for sensor inversion constraint imposition |
-
1993
- 1993-11-18 JP JP31426693A patent/JPH07140031A/ja active Pending
Cited By (15)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| GB2350432B (en) * | 1999-05-27 | 2002-12-18 | Bg Intellectual Pty Ltd | Method for tracing a gas leak |
| GB2350432A (en) * | 1999-05-27 | 2000-11-29 | Bg Intellectual Pty Ltd | Apparatus for tracing a gas leak |
| CN114754934B (zh) * | 2021-01-08 | 2024-05-28 | 中国石油天然气集团有限公司 | 气体泄漏检测方法 |
| CN114754934A (zh) * | 2021-01-08 | 2022-07-15 | 中国石油天然气集团有限公司 | 气体泄漏检测方法 |
| US12577871B2 (en) | 2022-08-03 | 2026-03-17 | Schlumberger Technology Corporation | Linear cut generation method for sensor inversion constraint imposition |
| US12480924B2 (en) | 2022-08-03 | 2025-11-25 | Schlumberger Technology Corporation | Automated record quality determination and processing for pollutant emission quantification |
| CN116429328B (zh) * | 2022-12-12 | 2023-09-12 | 宁波永大建设集团有限公司 | 一种外墙窗户渗漏检测方法、系统、存储介质及智能终端 |
| CN116429328A (zh) * | 2022-12-12 | 2023-07-14 | 宁波永大建设集团有限公司 | 一种外墙窗户渗漏检测方法、系统、存储介质及智能终端 |
| US12292310B2 (en) | 2022-12-15 | 2025-05-06 | Schlumberger Technology Corporation | Machine learning based methane emissions monitoring |
| WO2024155857A1 (en) * | 2023-01-19 | 2024-07-25 | Schlumberger Technology Corporation | Leak emissions sensor systems and processes for installing and using same |
| CN116296120A (zh) * | 2023-03-16 | 2023-06-23 | 中国科学院武汉岩土力学研究所 | 一种面向盐穴储氦库井场的氦气泄漏源定位方法 |
| CN116296120B (zh) * | 2023-03-16 | 2025-08-15 | 中国科学院武汉岩土力学研究所 | 一种面向盐穴储氦库井场的氦气泄漏源定位方法 |
| US12555125B2 (en) | 2023-06-09 | 2026-02-17 | Schlumberger Technology Corporation | Emission detecting camera placement planning using 3D models |
| US12254622B2 (en) | 2023-06-16 | 2025-03-18 | Schlumberger Technology Corporation | Computing emission rate from gas density images |
| CN117191928A (zh) * | 2023-08-08 | 2023-12-08 | 广州科易光电技术有限公司 | 气体泄漏检测方法、装置、电子设备及存储介质 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JPS5926064B2 (ja) | 輪郭画像の特徴抽出装置 | |
| JPH07140031A (ja) | ガス漏洩源領域、ガス漏洩源位置及びガス漏洩量の推定方法 | |
| CN110851952B (zh) | 利用固定位置传感器定位化学气体泄漏点的方法和设备 | |
| CN110823465B (zh) | 利用移动传感器定位化学气体泄漏点的方法和设备 | |
| AU2023202555A1 (en) | Method and system for locating and quantifying fugitive emission leaks | |
| JPH07198523A (ja) | ガス濃度データと風速データからのガス漏洩源位置とガス漏洩量の推定方法 | |
| CN114755367B (zh) | 一种环保污染监测方法、系统、设备及介质 | |
| JP2740718B2 (ja) | ガス、蒸気等の漏洩地点および漏洩量推定システム | |
| JP2767088B2 (ja) | 気体状物質の漏洩地点推定方法 | |
| CN118781186B (zh) | 用于地下管廊渗漏的快速定位方法及系统 | |
| US20040011124A1 (en) | Method for determining aerodynamic parameters and method for detecting failure of a probe used for determining aerodynamic parameters | |
| JPH07120344A (ja) | ガス漏洩領域及び漏洩量の推定方法 | |
| JPH0240967B2 (ja) | Gasuroeigennokenchihoho | |
| JPH0783786A (ja) | ガス漏洩点及び漏洩量の推定方法 | |
| JPH07198524A (ja) | ガス濃度データからのガス漏洩点位置とガス漏洩量の推定方法 | |
| JPH07333097A (ja) | ガス漏洩領域の推定方法 | |
| CN117077238A (zh) | 一种用于消防点位精准跟踪的方法及系统 | |
| US8190376B2 (en) | System and method for source identification for a chemical release | |
| CN111562795A (zh) | 变电站无人机巡检系统和方法 | |
| JPH0843239A (ja) | 障害物を考慮したガス漏洩領域及びガス漏洩量の推定方法 | |
| JPH0264437A (ja) | 粉粒体貯蔵場からの粉塵飛散量測定方法及びその装置、並びにそれを利用する粉塵飛散監視装置 | |
| CN114061501A (zh) | 一种架体立柱参数测量方法、装置、电子设备及存储介质 | |
| JP2731713B2 (ja) | 航跡追尾装置 | |
| JP2687466B2 (ja) | ガスの漏洩場所を推定する方法 | |
| JP2730876B2 (ja) | プラント内における微小漏洩ガス雰囲気での新規漏洩発見法 |