JPH06129941A - パイプラインの監視方法及びその装置 - Google Patents

パイプラインの監視方法及びその装置

Info

Publication number
JPH06129941A
JPH06129941A JP5145959A JP14595993A JPH06129941A JP H06129941 A JPH06129941 A JP H06129941A JP 5145959 A JP5145959 A JP 5145959A JP 14595993 A JP14595993 A JP 14595993A JP H06129941 A JPH06129941 A JP H06129941A
Authority
JP
Japan
Prior art keywords
pressure
time
value
straight line
wave
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.)
Granted
Application number
JP5145959A
Other languages
English (en)
Other versions
JP2565086B2 (ja
Inventor
Deii Uorutaazu Guretsugu
グレッグ・ディー・ウォルターズ
Bii Jiyaadeii Jieemusu
ジェームス・ビー・ジャーディー
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
JFE Engineering Corp
Original Assignee
NKK Corp
Nippon Kokan Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by NKK Corp, Nippon Kokan Ltd filed Critical NKK Corp
Publication of JPH06129941A publication Critical patent/JPH06129941A/ja
Application granted granted Critical
Publication of JP2565086B2 publication Critical patent/JP2565086B2/ja
Anticipated expiration legal-status Critical
Expired - Fee Related legal-status Critical Current

Links

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
    • G01M3/00—Investigating fluid-tightness of structures
    • G01M3/02—Investigating fluid-tightness of structures by using fluid or vacuum
    • G01M3/26—Investigating fluid-tightness of structures by using fluid or vacuum by measuring rate of loss or gain of fluid, e.g. by pressure-responsive devices, by flow detectors
    • G01M3/28—Investigating fluid-tightness of structures by using fluid or vacuum by measuring rate of loss or gain of fluid, e.g. by pressure-responsive devices, by flow detectors for pipes, cables or tubes; for pipe joints or seals; for valves ; for welds
    • G01M3/2807—Investigating fluid-tightness of structures by using fluid or vacuum by measuring rate of loss or gain of fluid, e.g. by pressure-responsive devices, by flow detectors for pipes, cables or tubes; for pipe joints or seals; for valves ; for welds for pipes
    • G01M3/2815—Investigating fluid-tightness of structures by using fluid or vacuum by measuring rate of loss or gain of fluid, e.g. by pressure-responsive devices, by flow detectors for pipes, cables or tubes; for pipe joints or seals; for valves ; for welds for pipes using pressure measurements
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
    • G01M3/00—Investigating fluid-tightness of structures
    • G01M3/02—Investigating fluid-tightness of structures by using fluid or vacuum
    • G01M3/04—Investigating fluid-tightness of structures by using fluid or vacuum by detecting the presence of fluid at the leakage point
    • G01M3/24—Investigating fluid-tightness of structures by using fluid or vacuum by detecting the presence of fluid at the leakage point using infrasonic, sonic, or ultrasonic vibrations
    • G01M3/243—Investigating fluid-tightness of structures by using fluid or vacuum by detecting the presence of fluid at the leakage point using infrasonic, sonic, or ultrasonic vibrations for pipes

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Pipeline Systems (AREA)
  • Testing Or Calibration Of Command Recording Devices (AREA)
  • Measurement Of Unknown Time Intervals (AREA)
  • Examining Or Testing Airtightness (AREA)
  • Measuring Fluid Pressure (AREA)

Abstract

(57)【要約】 【目的】 パイプラインに圧力変動を伴う事象が発生し
た場合に、たとえＳ／Ｎ比が低い圧力信号であったとし
ても、時間遅れなく、前記事象に起因する圧力先頭波の
到達時刻及び強度を正しく検出する。 【構成】 パイプライン中の流体に圧力変動を発生させ
る事象に起因して、流体を媒体として伝播する圧力先頭
波の到達時刻を検出する圧力先頭波の到達時刻検出方法
において、パイプライン中の流体圧力に関係する特性値
を測定し、測定された特性値を電気信号に変換し、所定
時間に亘る電気信号を記憶保持し、記憶保持された所定
時間の電気信号から統計的手法を用いて圧力先頭波の到
達時刻を決定する。

Description

【発明の詳細な説明】
【０００１】
【産業上の利用分野】本発明は、パイプラインにおける
漏洩現象や弁の開閉操作等の各種の事象の発生を監視す
るパイプラインの監視方法及びその装置に係わり、特
に、前記事象に起因して発生しパイプライン内の流体中
を伝播する圧力先頭波の到達時刻を正確に検出する圧力
先頭波の到達時刻検出方法と圧力先頭波の到達時刻検出
装置に関する。
【０００２】
【従来の技術】天然ガス等の気体や石油等の液体を輸送
するパイプラインにおいては、弁操作、ポンプの起動／
停止等の通常の運転状態における各事象、又は漏洩等の
事象が発生すると、その事象の発生した地点でパイプラ
イン内を通流する前述した各流体の圧力が急激に変化す
る。この圧力の瞬時の変化により、圧力波が発生し、発
生地点から上流側と下流側に向かって伝播していく。こ
の伝播速度はこのパイプライン中を流れる流体を伝搬す
る音の速度にほぼ等しい。
【０００３】パイプライン中を通流する流体の圧力を検
出する圧力検出器を取付け、この圧力検出器の出力信号
を常時監視し、適切な処理を行えば、前記各種の事象に
起因して発生する圧力波、すなわち圧力先頭波を検出で
きることになる。これにより、例えばパイプライン上に
互いに離間した２か所に圧力検出器を取付ける事によっ
て、各圧力検出器が各圧力先頭波を検出した時刻の時間
差でもって、漏洩等の事象の発生位置を推定したり、弁
操作やポンプの作動操作が正常に行われたかどうかを確
認できる。
【０００４】この事象の発生位置を求めるための一つの
技術がオッテンシュタインに付与された米国特許第 3,8
51,521号に開示されている。また、欧州特許出願第 0
042212 号の公開公報には、圧力先頭波の強度から事象
の発生位置を推定する技術が開示されている。
【０００５】
【発明が解決しようとする課題】一般に、パイプライン
中を通流する流体の圧力は絶えず変動している。その結
果、例えば図２に示すように、圧力検出器で検出された
圧力信号に高いレベルの雑音が混入して、Ｓ／Ｎ比が大
幅に低下し、圧力波の先端（即ち、圧力先頭波）を正確
に検出するための妨げとなる。
【０００６】なお、圧力先頭波の到達の検出を意図する
先行技術の例が、コビントンに付与された米国特許第
4,144,748号，4,091,658 号と 4,012,944号や、ファー
マに付与された米国特許第4,796,446 号に開示されてい
る。しかし、この開示された技術によると、検出された
圧力信号に含まれる雑音成分を除去するために、圧力先
頭波に続く圧力信号を著しいフィルタリングを介して観
察し、圧力先頭波通過後の圧力変化の傾向を検出するよ
うにしている。したがって、この技術は圧力先頭波が通
過（到達）したことを確認する技術であり、圧力先頭波
の到達時刻そのものを検出するものではない。
【０００７】コビントンの積分器やファーマのフィルタ
のような単一のフィルタを使えば、ある程度の雑音を除
去でき、圧力変化の傾向を把握するのに役立つ。しか
し、雑音を充分除去するためには、ローパスフィルタの
遮断周波数を例えば５Hz程度に設定する必要がある。し
かし、このように遮断周波数を低い値に設定すると、本
来の検出しようとする圧力先頭波の変化特性がフィルタ
の時定数に影響される。その結果、検出された圧力先頭
波の到達時刻が実際の到達時刻より遅れてしまい、正確
な到達時刻を検出できない問題がある。
【０００８】本発明はこのような事情に鑑みてなされた
ものであり、クロック信号等に同期して得られる各時刻
の圧力値のうちの所定個数の圧力値を統計処理し、この
統計的な変化量でもって圧力先頭波の到達を判定する事
によって、たとえＳ／Ｎ比が低い圧力信号であったとし
ても、時間遅れなく、正しい圧力先頭波の到達時刻を検
出できる圧力先頭波の到達時刻検出方法及び圧力先頭波
の到達時刻検出装置を提供することを目的とする。
【０００９】
【課題を解決するための手段】上記課題を解消するため
に、請求項１の発明は、パイプライン中の流体に圧力変
動を発生させる事象に起因して流体を媒体として伝播す
る圧力先頭波の到達時刻を検出する圧力先頭波の到達時
刻検出方法である。
【００１０】そして、この方法は、パイプライン中の流
体圧力に関係する特性値を測定する工程と、測定された
特性値を電気信号に変換する工程と、所定時間に亘る電
気信号を記憶保持する工程と、記憶保持された所定時間
の電気信号から統計的手法を用いて前記圧力先頭波の到
達時刻を決定する工程とで構成されている。
【００１１】また、請求項２の発明は、パイプライン中
の流体に圧力変動を発生させる事象に起因して流体を媒
体として伝播する圧力先頭波の到達時刻をパイプライン
上の所定位置で検出する圧力先頭波の到達時刻検出方法
である。
【００１２】そして、この方法は、パイプライン上の所
定位置で流体の圧力値を所定時間間隔で測定して、各時
刻における各圧力値からなる圧力値群を求める工程と、
圧力値群に含まれる複数の圧力値のうちから連続する所
定個数の圧力値からなるサブ圧力値群を指定する工程
と、各圧力値がそれぞれ最新の圧力値として得られたと
き、自己が最新値として所属するサブ圧力値群に含まれ
る全部の圧力値から所定位置における圧力の時間変化特
性を直線近似する工程と、時間変化特性が所定の閥値を
越えた時刻を特定する工程と、特定された時刻以前の所
定期間における各圧力値の時間変化特性を示す参照直線
を求める工程と、特定された時刻以降の所定期間におけ
る各圧力値の時間変化特性の最大勾配直線を求める工程
と、参照直線と最大勾配直線との交点を求める工程と、
交点に対応する時刻を圧力先頭波の到達時刻として出力
する工程とで構成されている。
【００１３】また、請求項３の発明は、上記請求項２の
発明に対して、最大勾配直線を求めた以降に得られた圧
力の時間変化特性が所定の変化閥値未満になった時点に
おける各圧力値の最新直線を求める工程と、最新直線と
参照直線との間隔で示される圧力値変化量を前記圧力先
頭波の強度として出力する工程とを加えたものである。
【００１４】さらに、請求項４の発明は、パイプライン
中の流体に圧力変動を発生させる事象に起因して流体を
媒体として伝播する圧力先頭波の到達時刻をパイプライ
ン上の所定位置で検出する圧力先頭波の到達時刻検出装
置である。
【００１５】そして、この装置は、パイプライン上の所
定位置で流体の圧力値を所定時間間隔で測定して、各時
刻における各圧力値からなる圧力値群を求める手段と、
圧力値群に含まれる複数の圧力値のうちから連続する所
定個数の圧力値からなるサブ圧力値群を指定する手段
と、各圧力値がそれぞれ最新の圧力値として得られたと
き、自己が最新値として所属するサブ圧力値群に含まれ
る全部の圧力値から所定位置における圧力の時間変化率
を直線近似する手段と、得られた時間変化特性と所定の
閥値とを比較することによって、圧力先頭波の到達時刻
を特定する手段とで構成されている。
【００１６】また、請求項５の発明は、請求項４の装置
における圧力先頭波の到達時刻を特定する手段を、圧力
の時間変化特性が所定の閥値を越えた時刻を特定する手
段と、特定された時刻以前の所定期間における各圧力値
の時間変化特性を示す参照直線を求める手段と、特定さ
れた時刻以降の所定期間における各圧力値の時間変化特
性の最大勾配直線を求める手段と、参照直線と最大勾配
直線との交点を求める手段と、交点に対応する時刻を圧
力先頭波の到達時刻として出力する手段とで構成してい
る、また、請求項６の発明は、上述した請求項４の装置
に対して、所定位置における圧力の時間変化特性から圧
力先頭波の強度を求める手段を加えたものである。
【００１７】さらに、請求項７の発明は、パイプライン
中の流体に圧力変動を発生させる事象に起因して流体を
媒体として伝播する圧力先頭波の到達時刻をパイプライ
ン上の所定位置で検出する圧力先頭波の到達時刻検出装
置である。
【００１８】そして、この装置は、パイプライン上の所
定位置における流体の圧力に関する特性値をアナログ信
号として検出する検出器と、標準時刻情報によりクロッ
ク信号を出力する基準時刻発生部と、検出器から出力さ
れる特性値のアナログ信号を、クロック信号を付加しＡ
／Ｄ変換して各時刻におけるデジタル信号として出力す
るＡ／Ｄ変換器と、Ａ／Ｄ変換器から出力される最新の
所定時間内に含まれる特性値の時間的変化を示す最適直
線を求める最適直線算出回路と、最適直線の勾配が所定
の閥値を越えたとき圧力先頭波検出信号を出力する比較
器と、圧力先頭波検出信号と基準時刻発生部から出力さ
れるクロック信号とから圧力先頭波の到達時刻を求める
到達時刻算出回路とで構成されている。
【００１９】
【作用】このように構成された圧力先頭波の到達時刻検
出方法及び到達時刻検出装置の基本動作を説明する。パ
イプライン内に圧力変動を伴う事象が発生した場合、圧
力波がパイプライン内の流体中を伝播していくが、この
圧力波（即ち圧力先頭波）はパイプライン内の流体の圧
力に関係する特性値を測定する事によって検出される。
測定された圧力に関係する特性値は電気信号に変換され
る。そして、測定された各圧力に関係する特性値は所定
時間に亘り記憶保持される。したがって、常時最新の所
定時間分の圧力に関係する特性値が記憶保持される。そ
して、この最新の所定時間分の圧力圧力に関係する特性
値から統計的手法を用いて圧力先頭波の到達時刻を決定
している。すなわち、常に最新の統計量が得られ、例え
ばこの統計量の変化度合い等によって圧力先頭波の到達
時刻が一義的に定まる。
【００２０】また、請求項２においては、上記手法にお
ける好ましいより具体的な手法を示している。すなわ
ち、パイプライン上の所定位置で流体の圧力値を所定時
間間隔で測定して、各時刻における各圧力値からなる圧
力値群を定義し、さらに、この圧力値群に含まれる複数
の圧力値のうちから連続する所定個数の圧力値からなる
サブ圧力値群を定義している。
【００２１】そして、このサブ圧力値群に所属する所定
個数の圧力値から所定位置（測定位置）における圧力の
時間変化特性が直線近似される。すなわち、最新の所定
個数の圧力値から測定位置における圧力の時間的変化を
示す直線が統計的に求められる。そして、この直線は、
新規に圧力値が測定される毎に、新たな直線に更新され
る。
【００２２】圧力の時間的変化特性を示すこの直線の勾
配が所定の閥値を越えると、圧力先頭波の影響による圧
力変動が生じ始めたと判断する。そして、この時点以前
の所定個数の圧力値から前記圧力変動が生じていない通
常の状態での各圧力値の時間変化特性を示す参照直線を
求める。
【００２３】次に、閥値を越えた後の所定期間における
各圧力値の時間変化特性の最大勾配直線を求める。そし
て、参照直線と最大勾配直線との交点に対応する時刻が
圧力先頭波の到達時刻となる。
【００２４】さらに、最大勾配直線を求めた以降に得ら
れた各圧力値の時間変化特性が所定の閥値未満になった
時点における各圧力値の最新直線を求め、この最新直線
と前記参照直線との間隔で示される圧力値変化量を圧力
先頭波の強度とする。
【００２５】すなわち、圧力先頭波が測定位置（所定位
置）を通過すると、この測定位置（所定位置）における
圧力値は圧力先頭波が到達する前に比較して大きく変化
するが、圧力先頭波が通過した後は、圧力値の時間変化
特性は元の通常状態時に近くなる。したがって、圧力値
の時間変化特性が所定の閥値未満に低下した時点で、こ
の時点における最新直線と先の参照直線との間隔が圧力
先頭波の圧力変化幅、すなわち圧力先頭波の強度に対応
する値となる。
【００２６】また、請求項７の発明においては、パイプ
ライン上の所定位置に配置された検出器でパイプライン
内の流体の圧力がアナログ信号として検出され、Ａ／Ｄ
変換器でもって、基準時刻発生部から出力されるクロッ
ク信号を付加して、デジタルの圧力値に変換される。こ
のＡ／Ｄ変換器から出力される最新の所定時間内に含ま
れる各時刻における各圧力値の時間的変化を示す最新直
線が最適直線算出回路で算出される。
【００２７】そして、比較器において最新直線の勾配が
所定の閥値を越えたと判定された場合に、圧力先頭波が
到達したと判断して、先頭波検出信号が出力される。さ
らに、到達時刻算出回路でもって、先頭波検出信号と基
準時刻情報から圧力先頭波の到達時刻が算出される。
【００２８】
【実施例】以下本発明の一実施例を図面を用いて説明す
る。図１は実施例の圧力先頭波の到達時刻検出方法を適
用した圧力先頭波の到達時刻検出装置の概略構成を示す
ブロック図である。
【００２９】図１に示すように、パイプライン２の所定
位置に、このパイプライン２内部の流体圧力をアナログ
電気信号に変換する圧力検出器１が取付けられている。
この圧力検出器１は、パイプの管壁の変形を測定する歪
み計、またはマイクロフォン、その他の圧力の過渡的事
象に感応する検出器であれば、どんな検出器であっても
よい。
【００３０】この実施例においては、最大流体圧力及び
雑音を処理するに十分な機能範囲を備え、更に、パイプ
ライン内における圧力先頭波の減衰に鑑み、測定が望ま
れる長距離のパイプラインにおける微小な圧力先頭波の
強度も分析するに充分な機能を備える高品質の工業用圧
力検出器が使用されている。
【００３１】また、その周波数応答は、過渡的事象の圧
力先頭波の形状を正確に測定するに充分なものでなけれ
ばならない。例えば、ガソリン輸送のパイプラインの場
合、少なくとも0.03 psiの分解能を持ち、また少なくと
も50 Hz の周波数応答を持ち、かつ 0〜 2000 psi の圧
力測定レンジを有する圧力検出器であれば、0.1 psiの
過渡的圧力変化の検出が可能である。
【００３２】一般に、圧力検出器１からの出力信号は、
流体の運動の乱れや検出器に起因する雑音を含んでい
る。この雑音は、Ａ／Ｄ変換器３のサンプリング周波数
の半分（ナイキスト周波数）より高い周波数を含む場合
もある。また、パイプライン２の圧力値の絶対値に相当
する信号のオフセット成分（直流成分）もある。
【００３３】これらの圧力信号の構成成分がデジタル信
号演算処理部５における圧力先頭波検出の処理動作を妨
げる。これらの不要成分をアナログ信号処理部６内に設
けられたフィルタを通過させて、できるだけ排除する必
要がある。このフィルタのパラメーターとしては、ノイ
ズの特性次第であるが、およそ0.1 から50 Hz の周波数
通過帯域のものが通常使用される。
【００３４】なお、このフィルタは、先行技術の説明で
述べたフィルタとはその機能が異なる。すなわち、先行
技術のフィルタは、遮断周波数がかなり低い値に設定さ
れたローパスフィルタである。したがって、圧力先頭波
の強度と到達時刻を正確に決定するため必要となる圧力
先頭波に関する情報を切り捨てることになる。
【００３５】本発明では、このような極めて有用な圧力
先頭波の情報を保持する広い周波数の通過帯域を持つフ
ィルタを使用する。勿論、一部の雑音は、このようなフ
ィルタで除去されずにそのまま通過する。しかし、本発
明では、以下に述べる統計的な演算の手法を用いて、こ
のような雑音を排除できる。
【００３６】実際、この発明の手法を採用することによ
って、パイプラインの中の例えば図２に示したpeak to
peakの雑音が、漏洩等を示す圧力先頭波の強度に近似し
たり又はこの強度を超えて、すなわちＳ／Ｎ比が１を下
回る程度である場合でも、圧力先頭波を正確に測定でき
る。
【００３７】圧力検出器１からの圧力値を示すアナログ
電気信号は、アナログ信号処理部６内のフィルタに導か
れる。Ａ／Ｄ変換器３は、フィルタから出力されたアナ
ログ信号を受信し、圧力先頭波の先端の形状を分析する
に充分なサンプリング周波数でデジタル信号に変換す
る。サンプリング周波数としては、上昇または下降して
いる圧力先頭波内に含まれる２乃至３個の圧力サンプル
値が採取できれば、通常はそれで充分である。このサン
プリング周波数は、パイプライン２の種類によるが、10
〜100 Hzの間で適当な値に設定される。この実施例にお
いては、２５Hzのサンプリング周波数が使用されてい
る。
【００３８】基準時刻発生部４により各圧力サンプル値
採取の正確な時刻が定まる。これは、圧力先頭波の到達
時刻の正確な決定に必要である。この実施例では、衛星
通信（ＧＰＳ）を一定周期で受信して、基準時刻を得
る。なお、衛星通信（ＧＰＳ）を利用する他にも、ＪＪ
Ｙ（標準時刻電波）のラジオ放送を受信したり、原子の
半減期を利用した原子時計等を利用することができる。
【００３９】デジタル信号演算処理部５は、Ａ／Ｄ変換
器３から出力された各圧力サンプル値を受取り、圧力先
頭波検出のアルゴリズムを実行して、到達時刻と強度と
を出力する。この実施例においては、デジタル信号演算
処理部５はマイクロコンピュータで構成されている。
【００４０】図４に示す圧力先頭波検出アルゴリズム
は、パイプライン２中を伝搬する圧力先頭波が流体圧力
の統計的に意味のある変化を示すことを前提として作動
する。すなわち、あるデータ群が他のデータ群に含まれ
るか否かを統計的に決定する手法である。
【００４１】この手法においては、最も最新の所定個数
の圧力サンプル値の統計的な代表値が母集団の統計的代
表値に比較して異なっている場合は、この所定個数の圧
力サンプル値は、異なる母集団に所属すると判断する。
したがって、パイプライン２の運転状態に変化があった
ことを意味する。
【００４２】なお、図４のフローチャートに関する以下
の議論は、負の圧力先頭波の検出について説明するが、
正の圧力先頭波の検出においても全く同様の手順で実行
できる。
【００４３】図４において、先ず、Ａ／Ｄ変換器３のサ
ンプリング周波数で定まる所定周期（サンプリング周
期）に相当する時間が経過して、Ａ／Ｄ変換器３から出
力される１個の圧力サンプル値を読取る（ステップ
１）。
【００４４】ステップ２でその読取った圧力サンプル値
が、例えば図３に示すように、次々に生じる最新のＮP
個の圧力サンプル値を常に含むメモリに蓄積される。こ
の実施例では、メモリに常時蓄積されている最新の圧力
サンプル値の個数ＮP は５００個である。したがって、
この個数ＮP は、１秒当たり２５サンプルの採取速度
（すなわち、２５Hzのサンプリング周波数）において
は、２０秒間の時間間隔、又は時間窓を表わす。
【００４５】各圧力サンプル値がメモリに蓄積される
と、ステップ３で、時間と圧力サンプル値とを直線的な
関数関係に置き換えるような最小自乗法を用いて、最新
のＮC個の圧力サンプル値に最もよく適合する最新直線
を算定する。但し、ここでは、個数ＮC は個数ＮP より
小さいものとする。この実施例では、ＮC は１２５個で
あり、５秒間の時間窓を表わす。
【００４６】この算出された最新直線の勾配はステップ
４で最新のＮh 個の勾配を含む勾配メモリに蓄積され
る。なお、この実施例では、個数Ｎh はＮP に等しい
（Ｎh ＝ＮP ）。この勾配メモリに記憶されているＮh
個全部の勾配の平均値Ｍh 及び標準偏差σh がステップ
５で計算される。
【００４７】ステップ６で最も最新に得られた勾配と平
均値Ｍh との差（勾配−Ｍh ）の絶対値が、標準偏差σ
h のＫS 倍を超えた場合、この最新の勾配は、Ｎh 個の
勾配群とは異なる母集団に含まれると推定される。すな
わち、この勾配は、通常の圧力状態とは異なる圧力先頭
波を意味すると判断される。なお、係数ＫS は、パイプ
ライン固有の雑音環境で最もよく機能するように調整す
る必要がある。この実施例では、この係数ＫS は約５〜
１０である。
【００４８】一旦、勾配と平均値Ｍh との差の絶対値が
この閾値（ＫS ・σh ）を超えると、ステップ７で過渡
的変化事象が既に検出されていたかどうかを判断するた
めに、二者択一式の偏差フラグの設定状態を調べる。偏
差フラグがまだ設定されていなければ、該当偏差フラグ
を設定する（ステップ８）。その後、この圧力低下の前
に生じた圧力状況を知る必要がある。そのため、この閾
値（ＫS ・σh ）を超えたとき、図３に示す参照直線が
Ｎr 個の各圧力サンプル値を用いて算出される。そし
て、算出された参照直線は一旦記憶される（ステップ
９）。
【００４９】この個数Ｎr は、監視対象パイプラインの
条件により変わるが、例えば、個数Ｎr は個数ＮC に等
しくてもよい。但し、この参照直線を求めるためのＮr
個の圧力サンプル値に対応する時間窓は、圧力低下部分
の採取時点を含んではならない。
【００５０】すなわち、Ｎi 個目の圧力サンプル値を採
取した時点で、このＮi 個目の圧力サンプル値を含むＮ
C 個の圧力サンプル値を代表する勾配と平均値Ｍh との
差の絶対値が前記閥値（ＫS ・σh ）を越えた場合に
は、前記参照直線がＮi 個目の圧力サンプル値を含む圧
力低下部分の影響を受けないようにするため、その参照
直線の時間窓はＮi-Nr個目からＮi 個目までとすること
はできない。
【００５１】この結果、具体的には、監視対象パイプラ
インの条件次第で変わるが、本発明によると、この参照
直線の時間窓はＮi-Roの点で終了するのが望ましい。な
お、ＲO は、Ｎi 個目の圧力サンプル値から、それ以前
へＲO 個遡ることを意味する。監視対象パイプラインの
条件により変わるが、ＲO は１０程度であるのが望まし
い。即ち、Ｎr 個の圧力サンプル値は、採取時点Ｎ
i-Ro-Nr から採取時点Ｎi-Roまでに採取されたものであ
る。この算出された参照直線の勾配ＳR 及び切片ＯＦＦ
R は前述したように後で使用するために蓄積される。
【００５２】ステップ６において、勾配の平均値Ｍh と
の差（勾配−Ｍh ）の絶対値が標準偏差σh のＫS 倍か
らなる閾値（ＫS ・σh ）を越えても、圧力サンプル値
の読取りは継続される。すなわち、ステップ１からステ
ップ７を繰返すことによって最適直線が算出される。
【００５３】そして、最も最新に得られたＮC 個の圧力
サンプル値を代表する直線の勾配と平均値Ｍh との差の
絶対値がこの閾値（ＫS ・σh ）を下回る値に入るまで
戻ったとき、ステップ６からステップ１０へ進む。そし
て、偏差フラグが設定されていた場合は、ステップ１１
において、過渡的変化事象がなお観察されるか否かを調
べる。具体的には、最新の直線の勾配と、その時点での
Ｍh 個の勾配群の平均値Ｍh との差（勾配−Ｍh ）の絶
対値が、閥値（Ｋe ・σh ）未満に低下しているか否か
を調べる。
【００５４】低下した場合は、過渡的変化事象が終了し
たと判断できる。この閥値（Ｋe ・σh ）は、この圧力
低下の終了の検出を最適に行うように調整される。な
お、この実施例においては、Ｋe は１に等しい。
【００５５】次に、発見されたこの圧力低下現象がパイ
プライン運転状態の著しい変化によるものであるかどう
かを評価する必要がある。具体的には、圧力先頭波の強
度をある閾値と比較する。この強度を求めるために、図
３に示すように、参照直線が時間軸の正方向に延ばさ
れ、この参照直線と最新のＮC 個の圧力サンプル値を代
表する勾配を持つ最新直線とのｙ軸方向の差異がＮC 個
の圧力サンプル値の中の最も古い採取点においてステッ
プ１２で計算される。
【００５６】ステップ１２は、偏差フラグがステップ７
で設定された後に１回だけ、即ち、ステップ１１が「Ｙ
ＥＳ」を示した直後にのみ実行される。この様な場合、
偏差フラグは、ステップ１７で解除される。偏差フラグ
がステップ７で再度設定されない限り、ステップ１０で
「ＮＯ」を示すことにより、ステップ１２に到達するこ
とはない。なお、ステップ１２に相当する回路が実施さ
れる場合の詳細は、後述する図６，図７に示す。
【００５７】ステップ１２で算出された図３におけるｙ
軸方向の差異の絶対値が、ステップ１３において、予め
定められた閾値ＫP を越えている場合には、この差異
は、ノイズではなく、過渡的変化事象に起因する圧力低
下によるものと判断できる。そして、この差異の絶対値
が過渡的変化事象による圧力先頭波の強度と見なされ
る。
【００５８】なお、差異の絶対値が閾値ＫP より小さい
場合には、ノイズに起因するものと判断され、無視され
る。この閾値ＫP は過渡的変化事象に対する本発明の検
出シテスムの感度を決めるために調整される。この実施
例においては、この閾値ＫPは約0.1 psi である。
【００５９】ステップ１３において閾値ＫP を超える変
化量（図３のｙ軸方向の差異）が検出されると、次のス
テップでは、その変化の正確な開始時点、即ち、圧力先
頭波の到達時刻が検出される。先ず、最大勾配直線がス
テップ１４で設定される。この最大勾配直線は、圧力が
実際に低下している時のＮM 個の圧力サンプル値を使用
して計算される全ての直線のうち最大勾配を持っている
直線である。
【００６０】かくして、ＮM 個の圧力サンプル値は、ス
テップ６で最初に閥値（ＫS ・σh）を越えた採取時点
と、ステップ１１が成立（ＹＥＳとなる）した時の直線
を決定するのに使用された圧力サンプル値の最も古い採
取時点との間に含まれる圧力サンプル値から得られる。
個数ＮM （この実施例では５である）は、個数ＮC より
相当小さい個数である。この最大勾配直線と前述の参照
直線との交差点に対応する時刻がステップ１５で算出さ
れる。この算出された時刻が圧力先頭波の到達時刻とし
てステップ１６において、先に求めた強度（変化量）と
共に出力される。
【００６１】なお、この採取時点とステップ１５に相当
する回路は、図６及び図７に記載されている。また、ス
テップ１３において、算出された変化量（強度）が閾値
ＫP を超えていない場合であっても、次の圧力先頭波の
到達を検出できるようにシステムをリセットするため
に、ステップ１７において偏差フラグが解除される。
【００６２】なお、上述した各定数ＮP ，ＮC ，Ｎh ，
Ｎr ，ＮM ，ＫS ，Ｋe ，ＫP の各値は、監視対象パイ
プラインの固有のノイズ環境下においても圧力先頭波の
到達時刻及び強度検出をより正確に実行するために調整
される。この調整は、熟練した運転技術者により経験的
に行われる。また、これは、運転技術者のパイプライン
全般の専門的知識、且つ、この発明の原理を具現化する
システムが設置された監視対象パイプラインにおいて、
このシステムの運用をすることで得られる知見により行
われる。
【００６３】これらの各値は、固定することもでき、ま
た、実際のシステム運用状態に合わせて、継続的に変更
することができる。例えば、本発明を使用しての運用経
験を積むことによって、この発明の使用結果が蓄積さ
れ、分析されるにつれて、このシステムは、前述の各定
数を誤警報を発しないぎりぎりの感度の高い水準に調整
等、自動的にシステム自体を最適に調節することができ
ると考えらる。
【００６４】図５は本発明の他の実施例に係わる圧力先
頭波の到達時刻検出装置の概略構成を示すブロック図で
ある。図１の実施例と同一部分には同一符号が付してあ
る。したがって、重複する部分の説明は省略する。この
実施例では、図１の実施例におけるアナログ処理部６が
除去されて、デジタル演算処理部５ａで全てが処理され
ている。
【００６５】パイプライン２内の流体の圧力は圧力検出
器１で検出されてアナログ信号に変換されてＡ／Ｄ変換
器３へ入力される。Ａ／Ｄ変換器３には基準時刻発生部
４内の一定周期で標準時刻を受領するクロック信号発振
器からのクロック信号が印加されている。したがって、
Ａ／Ｄ変換器３により各圧力信号を採取した時刻が基準
時刻により正確に決められる。Ａ／Ｄ変換器３の圧力信
号の取込み速度、すなわちサンプリング周波数は、サン
プリング定理により、圧力検出器１の出力信号中の検出
しようとするノイズを含む情報の周波数の少なくとも２
倍以上に設定される。
【００６６】デシタル演算処理部５ａで全てを処理する
実施例においては、サンプリング周波数は100 〜1000Hz
程度に設定する。ノイズを除去するためのデジタル・フ
ィルタはデジタル演算処理部５ａ内に組込まれている。
このデジタル・フィルタを実現するための詳細技術は既
に既知であるのでここでは説明を省略する。この技術に
関して、参考書として、オッペンハイム、アラン、Ｖ．
とシェイファー、ロナルド Ｗ．著、「離散的時系列信
号処理」、プレンティス ホール、エングルェウッド
クリフス社、Ｎ．Ｊ．，１９８９年等がある。
【００６７】図６は本発明の他の実施例に係わる圧力先
頭波の到達時間検出装置の概略構成を示すブロック図で
ある。パイプライン２の流体の圧力を検出する圧力検出
器１、Ａ／Ｄ変換器３、アナログ信号処理部６，基準時
刻発生部４は、図１に示す実施例と同一番号が付されて
いるので、重複する部分の説明は省略する。この実施例
は、図１の実施例におけるデジタル演算処理部５をハー
ド回路で構成したものである。
【００６８】Ａ／Ｄ変換器３から出力された各時刻にお
ける各圧力サンプル値は例えばシフトレジスタから構成
されたメモリ３０へ入力される。このメモリ３０は、例
えば５００個のレジスタからなり、最新の圧力サンプル
値が１個入力すると最も古い１個が掃出される。また、
全部でＮP （＝５００）個の圧力サンプル値を蓄積可能
である、このレジスタ３０のクロック端子には基準時刻
発生部４内の一定周期で基準時刻を受領するクロック信
号発振器からのクロック信号が入力されている。
【００６９】メモリ３０を構成するＮP （＝５００）個
の各レジスタから各圧力サンプル値が並列に出力され
る。すなわち、クロック信号が入力される毎にメモリ３
０からその時点における最新の５００個の各圧力サンプ
ル値が同時に出力される。メモリ３０から出力されたＮ
P （＝５００）個の圧力サンプル値は、最新直線算出回
路３２，参照直線算出回路３４，最大勾配直線算出回路
３６へ入力される。
【００７０】この３つの各最適直線を求める各回路３
２，３４，３６は同一構成であり２個の入力信号のみが
異なる。１番目の入力信号は、前記ＮP （＝５００）個
の圧力サンプル値のうち、これら各々の最適直線の算出
に用いる複数の圧力サンプル値のうちの先頭（時間軸で
表現すると最も新しい時点）の圧力サンプル値に対応す
る採取時刻の開始位置を示す信号である。２番目の信号
は、前記開始位置から何個の圧力サンプル値まで時間的
に遡って用いるかを示す、圧力サンプル値の個数Ｎc ，
Ｎr ，ＮM のいずれかを示す信号である。
【００７１】図７において、圧力サンプル値個数を示す
カウンタ４４は、外部からリセット信号が入力されると
計数値が０にクリアされ、クロック信号発生回路４５か
らクロック信号が入力される毎に、計数値を１ずつカウ
ントアップする。計数値は加算器４６でもって開始位置
が加算されて、マルチプレクサ回路４０へ印加される。
【００７２】したがって、マルチプレクサ回路４０へ入
力される信号は、圧力サンプル値の最適直線の算出に使
用する圧力サンプル値の先頭の採取時点から１採取時点
ずつ加算されていく信号となる。そして、カウンタ４４
は、外部から入力された圧力サンプル値個数ｎ（＝Ｎ
c ，Ｎr ，ＮM ）に達すると、計数値の計数動作を停止
する。
【００７３】マルチプレクサ回路４０にはメモリ３０か
ら出力されたＮP （＝５００）個の各圧力サンプル値が
入力されている。マルチプレクサ回路４０はＮP （＝５
００）個の圧力サンプル値のうち、加算器４６から出力
される信号の示す採取時点に対応する１個の圧力サンプ
ル値を次の算術ロジック・ユニット（以下ＡＬＵと略記
する）４２のｙ端子へ送出する。
【００７４】この場合、加算器４６の値が小さい程、シ
フトレジスタに後から入力された圧力サンプル値を示
す。したがって、加算器４６の値が０の場合は、メモリ
３０としてのシフトレジスタに最後に入力された最新の
圧力サンプル値がＡＬＵ４２へ送出される。
【００７５】したがって、カウンタ４４の計数値が０か
ら指定された圧力サンプル値個数に達すると、メモリ３
０に記憶されているＮP （＝５００）個の圧力サンプル
値のうち、外部から指定された採取時間開始位置から、
同じく指定された圧力サンプル値個数ｎ（＝Ｎc ，Ｎ
r ，ＮM ）までの各圧力サンプル値がＡＬＵ４２へ入力
されたことになる。
【００７６】また、ＡＬＵ４２のｘ入力端子には、カウ
ンタ４４から出力された計数値ｘが入力される。また、
ｎ端子には、外部から指定された圧力サンプル値個数ｎ
（＝Ｎc ，Ｎr ，ＮM ）が入力される。さらに、ＡＬＵ
４２には前記クロック信号発生回路４５からクロック信
号が入力され、さらに、リセット信号も入力される。
【００７７】ＡＬＵ４２のｘ入力端子に入力される各計
数値は各圧力サンプル値の各採取時刻位置を示している
ので、ＡＬＵ４２のｙ入力端子に入力される各圧力サン
プル値とｘ入力端子に入力される各採取時刻との関係
は、図３に示すように、ｙ軸を圧力サンプル値とし、ｘ
軸を時間軸とする（ｘ−ｙ）座標で表示することが可能
である。
【００７８】このＡＬＵ４２は一種のマイクロコンピュ
ータで構成されており、入力されたｎ個の各採取時刻ｘ
における各圧力サンプル値ｙとの直線関係を、例えば最
小自乗法を用いて求める。この最適直線はｙ＝Ｓｘ＋Ｏ
ffの形式で求まり、Ｓが最適直線の勾配であり、Ｏffが
ｙ軸を横切る位置を示す切片である。この勾配Ｓ及び切
片ＯFFは(1)(2)式で求まる。
【００７９】 勾配＝（ｎΣｘｙ−ΣｘΣｙ）／［ｎΣｘ2 −（Σｘ）2 ］ …(1) 切片＝（Σｘ2 Σｙ−ΣｘΣｘｙ）／［ｎΣｘ2 −（Σｘ）2 ］ …
（２） ＡＬＵ４２は算出した最適直線の勾配Ｓと切片Ｏｆｆを
出力する。適当な時間に、以下に説明するように、カウ
ンタ４４とＡＬＵ４２はリセットされる。ＡＬＵ４２は
クロック信号発生回路４５から出力されるクロック信号
に同期して作動する。
【００８０】与えられた圧力サンプル値の採取時点Ｎi
に対するｘの値は、第３図に示されるように過去に向か
って延びるｘ軸に対応する。最新の採取点は、ｘ＝０に
おけるものであり、先行する採取点は、ｘの正の数値側
におけるものである。最新の圧力サンプル値がシフトレ
ジスタからなるメモリ３０に入力される毎にｘ＝０とな
るレジスタ番地に蓄積される。クロック信号が入力する
毎に、各レジスタに記憶されている各圧力サンプル値は
１つずつシフトされる。すなわち、最も古い採取時刻の
圧力サンプル値が掃出される。
【００８１】図７に示した構造を有する各直線算出回路
３２，３４，３６にて算出された３つの最適直線は、そ
れぞれ最新直線、参照直線、最大勾配直線として出力さ
れる。各直線算出回路３２，３４，３６毎に各３つの最
適直線は、互いに異なるｘ軸座標系（直線の出発点でｘ
＝０）において計算され、表示される。３つの最適直線
は、次の方程式により定義される。
【００８２】 最新直線 ｙ＝ＳC ｘC ＋ＯffC …(3) 参照直線 ｙ＝ＳR ｘR ＋ＯffR …(4) 最大勾配直線 ｙ＝ＳM ｘM ＋ＯffM …(5) 但し、ＳとＯffは、それぞれ、勾配と切片を示す。
【００８３】図４のフローチャートのステップ１１の条
件を満たした時点において、すなわち、圧力先頭波が圧
力検出器１の設置位置を通過した場合に、上述した(3)
(4)(5) 式を用いて圧力先頭波の到達時刻ＸTDA と強度
Ａm を算出する必要がある。よって、上述した(3)(4)
(5) 式は次の関係を有する。
【００８４】 ｘc ＝ｘ（最新の直線は常に０から出発するから） …(6) ｘR ＝ｘ − ＰR …(7) ｘM ＝ｘ − ＰM …(8) 但し、図４に示すように、ＰR とＰM は、参照直線と最
大勾配直線の図中左方向へ延長する場合の開始位置を示
す採取時点である。
【００８５】方程式(3)(4)(5) に上記(6)(7)(8) を代入
して、上述した過渡的変化事象の終了時刻における各最
適直線は次の(10)(11)(12)式となる。 最新直線 ｙ＝ＳC ｘ＋ ＯffC …(10) 参照直線 ｙ＝Ｓ（ｘ−ＰR ）＋ＯffR …(11) 最大勾配直線 ｙ＝Ｓ（ｘ−ＰM ）＋ＯffM …(12) 図４に示すフローチャートのステップ１２で求めること
になる圧力先頭波の強度Ａm は、参照直線のｘ＝０方向
に延長された直線部分と最新直線のｘ＝ＮC の時点にお
けるｙ軸方向の差異として次のように計算される。
【００８６】 強度Ａm ＝［ＳR （ＮC −ＰR ）＋ＯffR ］−［ＳC ＮC ＋ＯffC ］ …(13) ＝ＳR ＮC −ＳR ＰR ＋ＯffR −ＳC ＮC −ＯffC …(14) 圧力先頭波の到達時刻ＸTDA は、参照直線と最大勾配直
線との交差点のｘ座標と考えられる。これは方程式(11)
(12)から次の(15)(16)式で算出される。
【００８７】 ＳR （ＸTDA −ＰR ）＋ＯFFR ＝ＳM （ＸTDA −ＰM ）＋ＯffM …(15) この式を解いて、到達時刻ＸTDA は、次の通りである。 ＸTDA ＝（ＳM ＰM −ＳR ＰR ＋ＯffR −ＯffM ）／（ＳM −ＳR ） …(16) さて、図６に戻る。最新直線算出回路３２は、上述した
ように動作し、ｘ＝０の出発採取時点の信号と例えば１
２５等のＮC の採取時点（圧力サンプル値数）の番号の
信号を受信する。その結果、最新直線算出回路３２は求
めた最適直線の勾配ＳC と切片ＯffC を出力する。
【００８８】同様に、参照直線回路３４は、例えば１０
等の最初の採取点Ｒ0 の信号を受信する。換言すれば、
上に述べた通り、出発点となる採取点は、図４のフロー
チャートにおけるステップ６において、（勾配−Ｍｈ）
の絶対値が所定の閾値（ＫS・σh ）を越えた時点（ｘ
＝０）から遡って１０番目の採取時点（Ｒ0 ＝１０，図
３においてはｘ軸上のＰR 点）が採取開始時点となる。
【００８９】そして、参照直線を算出するために用いる
圧力サンプル値個数はこの採取開始時点に先行するＮr
個である。その結果、参照直線算出回路３４は求めた参
照直線の勾配ＳR と切片ＯffR を出力する。
【００９０】最大勾配直線算出回路３６は最大勾配位置
カウンタ５０からの開始信号を受信する。また、この最
大勾配直線算出回路３６は例えば５個等の採取点数（計
算に使用する圧力サンプル値個数）ＮM を受信する。そ
して、この最大勾配直線算出回路３６は最大勾配直線を
算出して、勾配ＳM と切片ＯFFM を出力する。
【００９１】次に、図６に示す回路の各部の動作を順を
追って説明する。クロック信号発生回路５２からのクロ
ック信号に同期して、Ａ／Ｄ変換器３が圧力検出器１か
ら出力されたアナログ信号を各採取時点に対応するデジ
タルの各圧力サンプル値に変換して、シフトレジスタか
らなるメモリ３０へ格納していく。最新直線算出回路３
２はクロック信号に同期して、新たに算出した最新直線
の勾配ＳC と切片ＯFFC を出力する。
【００９２】勾配メモリ５４は、メモリ３０と同様にシ
フトレジスタで構成されており、最新直線算出回路３２
から順次出力される各勾配ＳC を受信し、例えば５００
個等の最新の勾配群Ｎh 個を記憶保持する。マイクロコ
ンピュータから構成されたＡＬＵ５６は、勾配メモリ５
４に蓄積されている５００個の勾配ＳC から、平均値Ｍ
h と標準偏差σh を計算する。
【００９３】また、図６の回路図において、次のよう
に、図４のフローチャートにおけるステップ６に相当す
る作業を行う。減算器５８は、最新直線算出回路３２か
ら出力された最新直線の勾配ＳC からＡＬＵ５６により
提供される勾配の平均値Ｍh を差し引く。そして、差の
絶対値｜ＳC −Ｍh ｜が絶対値算出回路６０で計算さ
れ、比較器６２の一方の入力端子へ入力される。
【００９４】この比較器６２の他方の入力端子には、乗
算器６４にて算出された係数ＫS とＡＬＵ５６から得た
標準偏差σh との乗算値（ＫS ・σh ）が入力されてい
る。比較器６２は両者の比較を行う（ステップ６）。差
の絶対値｜ＳC −Ｍh ｜が乗算値（ＫS ・σh ）を超え
ると、比較器６２はハイ（Ｈ）レベル信号を出力する。
その結果、アンドゲート６６の出力がハイ（Ｈ）レベル
となる。
【００９５】これは、初期状態においては偏差フラグ回
路７０に偏差フラグが設定されておらず、偏差フラグ回
路７０からの出力がロー（Ｌ）レベルであるので、イン
バータ６８からの出力はハイ（Ｈ）レベルであり、した
がって、比較器６２の出力がハイ（Ｈ）レベルに変化す
ると、アンドゲート６６の出力もハイ（Ｈ）レベルに変
化するからである。
【００９６】したがって、これらの回路が図４における
ステップ６の［ＹＥＳ］に対応する。アンドゲート６６
の出力がハイ（Ｈ）レベルに変化すると、例えばＲＳフ
リップフロップで構成された偏差フラグ回路７０におけ
る偏差フラグが設定される。
【００９７】偏差フラグ回路７０に偏差フラグが設定さ
れると、インバータ６８を介してアンドゲート７２の一
方の入力端子がロー（Ｌ）レベルに変化する。その結
果、参照直線算出回路３４のリセット端子には、クロッ
ク信号発生回路５２から新たなクロック信号は入力され
ない。
【００９８】したがって、偏差フラグが設定されている
限り、参照直線算出回路３４は新たな参照直線の算出処
理を停止し、偏差フラグが設定された時点における参照
直線の勾配ＳR と切片ＯFFR を継続して出力する。
【００９９】また、アンドゲート６６が成立すると、参
照直線位置カウンタ１０２がリセットされ、圧力サンプ
ル値の採取時点ｘの計数値を一旦０に戻す。また、図４
のフローチャートにおけるステップ１１に相当する手順
は、図６の回路において、次のように実施される。
【０１００】絶対値算出回路６０から出力される絶対値
｜勾配ＳC −Ｍh ｜は比較器７４の一方の入力端子へ入
力される。この比較器７４の他方の入力端子には、乗算
器７６にて算出された係数Ｋe とＡＬＵ５６から得た標
準偏差σh との乗算値（Ｋe・σh ）が入力されてい
る。比較器７４は両者の比較を行う（ステップ１１）。
【０１０１】乗算値（Ｋe ・σh ）が差の絶対値｜ＳC
−Ｍh ｜を超えると、比較器７４はハイ（Ｈ）レベル信
号を出力する。この時、偏差フラグ回路７０からの出力
は偏差フラグが既に設定されているために、ハイ（Ｈ）
レベルとなっているので、アンドゲート７８の出力がハ
イ（Ｈ）レベルとなる。アンドゲート７８の出力がハイ
（Ｈ）レベルとなると、偏差フラグ回路７０の偏差フラ
グがリセット（解除）される。
【０１０２】また、アンドゲート７８のハイレベル信号
は遅延回路８１を介してアンドゲート８０の一方の入力
端子へ入力される。アンドゲート８０の他方の入力端子
には比較器８２の出力信号が入力されている。比較器８
２は、絶対値算出回路８３から出力される前述した(14)
式で示す強度Ａm の絶対値｜Ａm ｜（変化量）と予め設
定された閥値ＫP とを比較し、 強度Ａm ＝ＳR ＮC −ＳR ＰR ＋ＯffR −ＳC ＮC −ＯffC …(14) この強度が閥値を越えた場合には、ハイ（Ｈ）レベル信
号をアンドゲート８０へ送出する。アンドゲート８０が
成立すると、ラッチ回路９８へラッチ信号が入力され
る。
【０１０３】次に、絶対値算出回路８３に入力される強
度Ａm の算出手順を説明する。乗算器８４は、最新直線
算出回路３２から出力される勾配ＳC と予め設定されて
いる圧力サンプル値個数ＮC との乗算値ＳC ・ＮC を算
出する。乗算値ＳC ・ＮC は次の加算器８６で、最新直
線算出回路３２から出力される切片ＯFFC に加算され
る。
【０１０４】また、別の乗算器９０は参照直線算出回路
３４から出力される勾配ＳR に同じく圧力サンプル値数
ＮC を乗算して、乗算値ＳR ・ＮC を次の減算器８８へ
送出する。減算器８８は、乗算器９０からの乗算値ＳR
・ＮC と加算器８６からの加算値との差 ＳR ＮC −（ＳC ＮC ＋ＯFFC ） を算出して、次の加算器９６へ送出する。
【０１０５】また、別の乗算器９４は、参照直線位置カ
ウンタ１０２から出力される計数値ＰR に、参照直線算
出回路３４から出力される勾配ＳR を乗算し、乗算値Ｓ
R ・ＰR を次の減算器９２へ送出する。この減算器９２
は、参照直線算出回路３４の切片ＯffR から乗算値ＳR
・ＰR を減算する。減算値を次の各加算器９６，１１０
へ送出する。加算器９６は、前述した(14)式で示す強度
Ａm を算出する。
【０１０６】加算器９６から出力された圧力先頭波の強
度Ａm はラッチ回路９８に入力される。アンドゲート７
８と比較器８２が共にハイレベル信号を出力するとアン
ドゲート８０が成立し、ラッチ回路９８が、前記圧力先
頭波の強度Ａm 及び加算器１１２から出力された圧力先
頭波の到達時刻ＴTDA をラッチして、出力する。
【０１０７】次に、圧力先頭波の到達時刻ＴTDA の算出
手順を説明する。最大勾配直線算出回路３６は、ｘ＝Ｎ
C とｘ＝ＰR との間に存在するＮM 個（例えば５個）の
採取時点の各圧力サンプル値を用いて、図３に示す最大
勾配直線を算出して、その勾配ＳM 及び切片ＯffM を出
力する。
【０１０８】最大勾配直線位置カウンタ５０は、図６に
示すように、アンドゲート７８のハイレベル出力に応動
してリセットされる。アンドゲート７８からハイレベル
信号が出力されると、最大勾配直線を設定するための各
圧力サンプル値の採取時点が図３におけるｘ＝ＮC の点
から始まる。クロック信号発生回路５１から出力される
クロック信号に同期して、最大勾配直線位置カウンタ５
０は入力されている初期値ｘ＝ＮC から計数値を増加し
ていく。
【０１０９】最大勾配直線算出回路３６の内部クロック
信号発生回路４５の周波数は、全ての計算がクロック信
号発生回路５１の１パルス分の周期内に十分終了できる
程に十分に高いものである。同様に、クロック信号発生
回路５１の周波数は、クロック信号発生回路５２のクロ
ック信号に同期して最大勾配直線位置カウンタ５０の作
動開始から終了までに亘り、全ての圧力サンプル値の採
取時点に対する計算を完全に行う程に十分に高いもので
ある。
【０１１０】最大勾配直線位置カウンタ５０の作動範囲
の最後の採取時点は、参照直線位置カウンタ１０２によ
って設定されている。参照直線位置カウンタ１０２はア
ンドゲート６６が成立したことにより０にリセットされ
る。すなわち、ステップ６の条件が成立している。図３
においては、ｘ＝ＰR −Ｒ0 の時である。参照直線位置
カウンタ１０２は時間の経過とともに、すなわち、ｘｙ
座標の原点が右に動くにつれて、前記ｘ（＝ＰR −Ｒ
0 ）に動いた分だけの個数を加えていく。
【０１１１】アンドゲート７８が成立したとき（ステッ
プ１１の条件が適合したとき）、参照直線位置カウンタ
１０２は、参照直線が計算され始めた採取時点から演算
処理されてきた採取時点の正確な番号を登録する（即
ち、これは、１回目のステップ６に対応する）。
【０１１２】最大勾配直線算出回路３６にて計算され
た、最大勾配直線の勾配ＳM と切片ＯFFM とはこの最大
勾配直線のｘ軸座標上の位置計算値と一緒にラッチ回路
１００に入力される。最大勾配直線に対する位置計算値
は、最大勾配直線位置カウンタ５０から出力される。ラ
ッチ回路１００は、入力に際して、最新の勾配を、既に
入力されている勾配と比較する。そして、これらの２つ
の直線の勾配のうちで急な勾配を持つ方の直線に対応す
る勾配、切片、時刻位置データを蓄積する。
【０１１３】この手順は、最終点（即ち、例えばｘ＝Ｐ
R −Ｒ0 ）のＸ軸座標の数値に達するまで繰り返され
る。したがって、最終的に、このラッチ回路１００にラ
ッチされて出力される最大勾配は、ＳM であり、これに
対応する切片と位置計算値は、それぞれ、ＯFFM とＰM
である。
【０１１４】乗算器１０４、減算器１０６、１０８及び
加算器１１０、除算器１１２は、ラッチ回路９８の入力
時に、圧力先頭波の到達時刻ＸTDA の数値を、前述した
式(16)に基づいて算出する。算出された圧力先頭波の到
達時刻ＸTDA は、前述したように、アンドゲート８０
が、ラッチ回路９８にラッチ信号を入力した時点で、ラ
ッチされ、出力される。
【０１１５】なお、遅延回路８１は、アンドゲート７８
が成立したときに使用され、ラッチ回路１００の入力時
に発信される信号は、最大勾配直線位置カウンタ５０と
最大勾配直線算出回路３６が再設定された後の時点にの
み（即ち、かかる遅延の後に）、ラッチ回路９８が作動
した場合に有用である。すなわち、アンドゲート７４が
成立した後に、最大勾配直線位置カウンタ５０と最大勾
配直線算出回路３６が正常に動作が終了するまでの時間
を稼ぐ回路である。
【０１１６】本発明をもってすれば、簡単な設置形態
で、圧力先頭波の到達時刻の検出が可能であり、この到
達時刻情報と、漏洩の位置及び規模を発見するのに役立
つ強度情報とともに、警報を発することに利用できる。
もっと高度な適用形態では、パイプライン上の幾つかの
地点のそれぞれから時刻、強度及びその他の情報を中央
ステーションに送ることを含む。この中央ステーション
では、運転技術者、或いは、マスターコンピュータによ
り前記情報を利用して漏洩あるいはその他の事象（例え
ば弁の開閉等）の発生位置を特定することができる。
【０１１７】なお、弁の開閉等を監視する利用形態の場
合、圧力検出器１の設置地点は、監視対象とする弁から
数マイル離れたパイプライン上の都合のよい場所に設定
することができる。弁が作動の指令を受ける度に、適当
な強度（極性を含む）と到達時刻の過渡的圧力事象を、
本発明によって検出することができ、単純な論理的比較
によって、弁に対する指令が正しく行われていることを
確認できる。
【０１１８】ここに述べたのは、圧力先頭波の到達時刻
と強度が、一旦、多くのノイズを含むパイプラインの圧
力信号から抽出されると、いろいろな面に適用できる多
くの実施例のうちのごく僅かのものである。
【０１１９】また、前記の如く、実施例について詳しく
述べてきたが、これに対する種々の改良がこの技術分野
に熟練している通常の知識を有する者によれば、容易に
行うことができる。例えば、直線以外に曲線を使用して
圧力・時間の関係を近似することもできる。即ち、近似
のために、「最小自乗法」以外の何らかの基準に基づい
て、曲線による適合手順を使用することもできようが、
これらの修正やその他の修正は、本発明の範疇に属する
ものである。
【０１２０】
【発明の効果】以上説明したように本発明の圧力先頭波
の到達時刻検出方法及び圧力先頭波の到達時刻検出装置
によれば、クロック信号等に同期して所定時間毎に得ら
れる各時刻の圧力サンプル値のうちの常時最新の所定個
数の圧力サンプル値を処理し、これにより圧力先頭波の
到達時刻及び圧力先頭波の強度を算出することができ
る。したがって、たとえＳ／Ｎ比が低い圧力信号であっ
たとしても、時間遅れなく、正しい圧力先頭波の到達時
刻及び強度を検出できる。
【図面の簡単な説明】
【図１】 本発明の一実施例に係わる圧力先頭波の到達
時刻検出装置の概略構成を示すブロック図。
【図２】 典型的なパイプライン上の１箇所の地点にお
ける圧力波形の時間経過図。
【図３】 同実施例装置を用いて圧力先頭波の到達時刻
と強度を決定するために使用される諸直線を示す説明
図。
【図４】 同実施例装置の動作を示す流れ図。
【図５】 本発明の他の実施例に係わる圧力先頭波の到
達時刻検出装置の概略構成を示すブロック図。
【図６】 本発明のさらに別の実施例に係わる圧力先頭
波の到達時刻検出装置の概略構成を示すブロック図。
【図７】 同実施例装置の要部を取出して示すブロック
図。
【符号の説明】
１…圧力検出器、２…パイプライン、３…Ａ／Ｄ変換
器、４…基準時刻発生部、５，５ａ…デジタル演算処理
部、６…アナログ信号処理部、３０…メモリ、３２…最
新直線算出回路、３４…参照直線算出回路、３６…最大
勾配直線算出回路、５０…最大勾配直線位置カウンタ、
４５，５１，５２…クロック信号発生回路、５４…勾配
メモリ、７０…偏差フラグ回路、９８，１００…ラッチ
回路、１０２…参照直線位置カウンタ。

Claims (7)

    【特許請求の範囲】
  1. 【請求項１】 パイプライン中の流体に圧力変動を発生
    させる事象に起因して前記流体を媒体として伝播する圧
    力先頭波の到達時刻を検出する方法において、 前記パイプライン中の流体圧力に関係する特性値を測定
    する工程と、 この測定された特性値を電気信号に変換する工程と、 所定時間に亘る前記電気信号を記憶保持する工程と、 この記憶保持された所定時間の電気信号から統計的手法
    を用いて前記圧力先頭波の到達時刻を決定する工程とを
    備えた圧力先頭波の到達時刻検出方法。
  2. 【請求項２】 パイプライン中の流体に圧力変動を発生
    させる事象に起因して前記流体を媒体として伝播する圧
    力先頭波の到達時刻を前記パイプライン上の所定位置で
    検出する方法において、 前記パイプライン上の所定位置で前記流体の圧力値を所
    定時間間隔で測定して、各時刻における各圧力値からな
    る圧力値群を求める工程と、 前記圧力値群に含まれる複数の圧力値のうちから連続す
    る所定個数の圧力値からなるサブ圧力値群を指定する工
    程と、 前記各圧力値がそれぞれ最新の圧力値として得られたと
    き、自己が最新値として所属するサブ圧力値群に含まれ
    る全部の圧力値から前記所定位置における圧力の時間変
    化特性を直線近似する工程と、 この時間変化特性が所定の閥値を越えた時刻を特定する
    工程と、 この特定された時刻以前の所定期間における各圧力値の
    時間変化特性を示す参照直線を求める工程と、 前記特定された時刻以降の所定期間における各圧力値の
    時間変化特性の最大勾配直線を求める工程と、 前記参照直線と最大勾配直線との交点を求める工程と、 この交点に対応する時刻を前記圧力先頭波の到達時刻と
    して出力する工程とを備えた圧力先頭波の到達時刻検出
    方法。
  3. 【請求項３】 前記最大勾配直線を求めた以降に得られ
    た圧力の時間変化特性が所定の閥値未満になった時点に
    おける各圧力値の最新直線を求める工程と、この最新直
    線と前記参照直線との間隔で示される圧力値変化量を前
    記圧力先頭波の強度として出力する工程とを備えた請求
    項２記載の圧力先頭波の到達時刻検出方法。
  4. 【請求項４】 パイプライン中の流体に圧力変動を発生
    させる事象に起因して前記流体を媒体として伝播する圧
    力先頭波の到達時刻を前記パイプライン上の所定位置で
    検出する圧力先頭波の到達時刻検出装置において、 前記パイプライン上の所定位置で前記流体の圧力値を所
    定時間間隔で測定して、各時刻における各圧力値からな
    る圧力値群を求める手段と、 前記圧力値群に含まれる複数の圧力値のうちから連続す
    る所定個数の圧力値からなるサブ圧力値群を指定する手
    段と、 前記各圧力値がそれぞれ最新の圧力値として得られたと
    き、自己が最新値として所属するサブ圧力値群に含まれ
    る全部の圧力値から前記所定位置における圧力の時間変
    化特性を直線近似する手段と、 この得られた時間変化特性と所定の閥値とを比較するこ
    とによって、前記圧力先頭波の到達時刻を特定する手段
    とを備えた圧力先頭波の到達時刻検出装置。
  5. 【請求項５】 前記圧力先頭波の到達時刻を特定する手
    段は、 前記圧力の時間変化特性が所定の閥値を越えた時刻を特
    定する手段と、 この特定された時刻以前の所定期間における各圧力値の
    時間変化特性を示す参照直線を求める手段と、 前記特定された時刻以降の所定期間における各圧力値の
    時間変化特性の最大勾配直線を求める手段と、 前記参照直線と最大勾配直線との交点を求める手段と、 この交点に対応する時刻を前記圧力先頭波の到達時刻と
    して出力する手段とを備えた請求項４記載の圧力先頭波
    の到達時刻検出装置。
  6. 【請求項６】 前記所定位置における圧力の時間変化特
    性から前記圧力先頭波の強度を求める手段を備えた請求
    項４記載の圧力先頭波の到達時刻検出装置。
  7. 【請求項７】 パイプライン中の流体に圧力変動を発生
    させる事象に起因して前記流体を媒体として伝播する圧
    力先頭波の到達時刻を前記パイプライン上の所定位置で
    検出する圧力先頭波の到達時刻検出装置において、 前記パイプライン上の所定位置における前記流体の圧力
    に関する特性値をアナログ信号として検出する検出器
    と、 標準時刻情報によりクロック信号を出力する基準時刻発
    生部と、 前記検出器から出力される特性値のアナログ信号を、前
    記クロック信号を付加しＡ／Ｄ変換して各時刻における
    デジタル信号として出力するＡ／Ｄ変換器と、 このＡ／Ｄ変換器から出力される最新の所定時間内に含
    まれる特性値の時間的変化を示す最適直線を求める最適
    直線算出回路と、 前記最適直線の勾配が所定の閥値を越えたとき圧力先頭
    波検出信号を出力する比較器と、 この圧力先頭波検出信号と前記基準時刻発生部から出力
    されるクロック信号とから前記圧力先頭波の到達時刻を
    求める到達時刻算出回路とを備えた圧力先頭波の到達時
    刻検出装置。
JP5145959A 1992-10-16 1993-06-17 パイプラインの監視方法及びその装置 Expired - Fee Related JP2565086B2 (ja)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
US962524 1992-10-16
US07/962,524 US5388445A (en) 1992-10-16 1992-10-16 Method for determining arrival and amplitude of a wave front and apparatus therefor

Publications (2)

Publication Number Publication Date
JPH06129941A true JPH06129941A (ja) 1994-05-13
JP2565086B2 JP2565086B2 (ja) 1996-12-18

Family

ID=25506010

Family Applications (1)

Application Number Title Priority Date Filing Date
JP5145959A Expired - Fee Related JP2565086B2 (ja) 1992-10-16 1993-06-17 パイプラインの監視方法及びその装置

Country Status (2)

Country Link
US (1) US5388445A (ja)
JP (1) JP2565086B2 (ja)

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2000266626A (ja) * 1999-03-16 2000-09-29 Techno Craft Co Ltd 同時多点測定装置
JP2008064494A (ja) * 2006-09-05 2008-03-21 Tsurumi Seiki:Kk 漏水判定装置、漏水判定方法、圧力波センサ
JP6115001B1 (ja) * 2016-04-06 2017-04-19 光俊 菅原 センサ信号検出装置
JP2019012050A (ja) * 2017-07-03 2019-01-24 国立研究開発法人農業・食品産業技術総合研究機構 埋設管水路の漏水原因モニタリング装置、及び埋設管水路の漏水原因モニタリング方法
WO2021106966A1 (ja) * 2019-11-28 2021-06-03 株式会社フジキン 流体制御機器の動作情報収集システム、流体制御機器、流体制御機器の動作情報収集方法、及びコンピュータプログラム
US20220275610A1 (en) * 2019-08-02 2022-09-01 The University Of Adelaide Method and system to monitor pipeline condition

Families Citing this family (22)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5948969A (en) * 1997-10-20 1999-09-07 Vista Research, Inc. Methods for measuring the flow rate due to a leak in a pressurized pipe system
US6289723B1 (en) * 1999-03-04 2001-09-18 Robert L. Leon Detecting seal leaks in installed valves
US6389881B1 (en) 1999-05-27 2002-05-21 Acoustic Systems, Inc. Method and apparatus for pattern match filtering for real time acoustic pipeline leak detection and location
AUPS236502A0 (en) * 2002-05-16 2002-06-13 Peterson, Paul Method of monitoring measurements
DE10242162A1 (de) * 2002-09-10 2004-03-18 Cegelec Anlagen- Und Automatisierungstechnik Gmbh & Co. Kg Verfahren zur Leckprüfung von Rohren und Rohrsystemen
AU2003902318A0 (en) * 2003-05-14 2003-05-29 Vision Fire And Security Pty Ltd Improved Sensing Apparatus And Method
US7219553B1 (en) * 2003-09-11 2007-05-22 Loren Worthington Dynamic transient pressure detection system
FR2883979B1 (fr) * 2005-03-29 2007-05-11 Areva T & D Sa Procede de detection de la position d'un front d'onde dans un signal recu par un detecteur
GB2423562B (en) * 2005-05-27 2007-01-17 Brinker Technology Ltd Determining leak location and size in pipes
US8365831B2 (en) * 2007-12-18 2013-02-05 Exxonmobil Upstream Research Company Determining connectivity architecture in 2-D and 3-D heterogeneous data
CA2728561A1 (en) * 2008-06-27 2009-12-30 Exxonmobil Research And Engineering Company A method and apparatus for real time enhancing of the operation of a fluid transport pipeline
US10352505B2 (en) * 2008-06-27 2019-07-16 Exxonmobil Research And Engineering Company Method and apparatus for real time enhancing of the operation of a fluid transport pipeline
FR2951561B1 (fr) 2009-10-20 2011-12-09 Areva T & D Sas Procede de detection de la position d'un front d'onde correspondant a un evenement dans un signal recu par un detecteur
RU2565112C2 (ru) * 2013-07-18 2015-10-20 Федеральное государственное бюджетное образовательное учреждение высшего образования "Тюменский государственный нефтегазовый университет" (ТюмГНГУ) Способ диагностики герметичности магистрального трубопровода
US10036684B2 (en) * 2013-11-12 2018-07-31 Nec Corporation Leakage determination system and leakage determination method
LT2963405T (lt) * 2014-07-03 2020-05-25 Kamstrup A/S Komunalinio tinklo slėgio jutiklio įtaisas
US10481036B2 (en) 2015-04-29 2019-11-19 Medeng Research Institute Ltd. Pipeline leak detection system
CN104930355B (zh) * 2015-06-09 2017-10-13 段焕丰 适用于城市供水管道系统的在线无损检测方法及其装置
WO2018052675A1 (en) 2016-09-19 2018-03-22 Exxonmobil Research And Engineering Company A method and apparatus for real time enhancing of the operation of a fluid transport pipeline
PL3607290T3 (pl) 2017-04-05 2022-08-22 Tenova Goodfellow Inc. Sposób i urządzenie do wykrywania akustycznego wycieków płynów
RU2700491C1 (ru) * 2019-01-10 2019-09-17 Федеральное государственное унитарное предприятие "Российский Федеральный ядерный центр - Всероссийский научно-исследовательский институт экспериментальной физики" (ФГУП "РФЯЦ-ВНИИЭФ") Способ обнаружения нештатной ситуации на многониточном магистральном трубопроводе
RU2758876C1 (ru) * 2020-12-29 2021-11-02 Акционерное Общество "Атомэнергопроект" Способ определения уровня утечки через негерметичный затвор шарового крана запорно-регулирующей арматуры в рабочем режиме и устройство для его осуществления

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS502990A (ja) * 1973-04-09 1975-01-13
JPS5839975A (ja) * 1981-09-02 1983-03-08 Mitsubishi Electric Corp 磁気異常自動検出器
JPS5934165A (ja) * 1982-08-20 1984-02-24 Fujitsu Ltd 入力デ−タの立上り検出方式
JPS6042632A (ja) * 1983-08-05 1985-03-06 エレクトリシテイ　デ　フランス　セルビスナシヨナル 音響信号の急激的な発生を検出する装置
JPS63276182A (ja) * 1987-05-07 1988-11-14 Ricoh Co Ltd 図形デ−タの近似処理方式
JPH0472537A (ja) * 1989-09-19 1992-03-06 Tokyo Gas Co Ltd 配管の異常監視装置

Family Cites Families (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US3903729A (en) * 1970-12-30 1975-09-09 Taft Broadcasting Corp Method and apparatus for detecting a break or other occurrence in a pipeline containing gas under pressure
US3851521A (en) * 1973-01-19 1974-12-03 M & J Valve Co System and method for locating breaks in liquid pipelines
US4144743A (en) * 1974-12-09 1979-03-20 Shafer Valve Company Electronic fluid pipeline leak detector and method
US4091658A (en) * 1974-12-09 1978-05-30 Shafer Valve Company Electronic fluid pipeline leak detector and method
US4012944A (en) * 1974-12-09 1977-03-22 Shafer Valve Company Electronic fluid pipeline leak detector and method
US4198855A (en) * 1978-03-06 1980-04-22 The Japan Steel Works Ltd. Line break detector
EP0042212A1 (en) * 1980-06-16 1981-12-23 Imperial Chemical Industries Plc Pipeline leak location
JPS58100730A (ja) * 1981-12-11 1983-06-15 Sumitomo Metal Ind Ltd パイプラインの漏洩検知方法およびその装置
EP0186478A3 (en) * 1984-12-25 1987-09-23 Nippon Kokan Kabushiki Kaisha Method and apparatus for detecting leaks in a gas pipe line
US4796466A (en) * 1987-02-17 1989-01-10 Ed Farmer System for monitoring pipelines
JPS649237A (en) * 1987-03-24 1989-01-12 Kanegafuchi Chemical Ind Heat-resistant polyimide film
JPH02147932A (ja) * 1988-11-30 1990-06-06 Nkk Corp パイプラインの漏洩検知方法
JPH02147931A (ja) * 1988-11-30 1990-06-06 Nkk Corp パイプラインの漏洩検知方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS502990A (ja) * 1973-04-09 1975-01-13
JPS5839975A (ja) * 1981-09-02 1983-03-08 Mitsubishi Electric Corp 磁気異常自動検出器
JPS5934165A (ja) * 1982-08-20 1984-02-24 Fujitsu Ltd 入力デ−タの立上り検出方式
JPS6042632A (ja) * 1983-08-05 1985-03-06 エレクトリシテイ　デ　フランス　セルビスナシヨナル 音響信号の急激的な発生を検出する装置
JPS63276182A (ja) * 1987-05-07 1988-11-14 Ricoh Co Ltd 図形デ−タの近似処理方式
JPH0472537A (ja) * 1989-09-19 1992-03-06 Tokyo Gas Co Ltd 配管の異常監視装置

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2000266626A (ja) * 1999-03-16 2000-09-29 Techno Craft Co Ltd 同時多点測定装置
JP2008064494A (ja) * 2006-09-05 2008-03-21 Tsurumi Seiki:Kk 漏水判定装置、漏水判定方法、圧力波センサ
JP6115001B1 (ja) * 2016-04-06 2017-04-19 光俊 菅原 センサ信号検出装置
WO2017175318A1 (ja) * 2016-04-06 2017-10-12 光俊 菅原 センサ信号検出装置
JP2019012050A (ja) * 2017-07-03 2019-01-24 国立研究開発法人農業・食品産業技術総合研究機構 埋設管水路の漏水原因モニタリング装置、及び埋設管水路の漏水原因モニタリング方法
US20220275610A1 (en) * 2019-08-02 2022-09-01 The University Of Adelaide Method and system to monitor pipeline condition
WO2021106966A1 (ja) * 2019-11-28 2021-06-03 株式会社フジキン 流体制御機器の動作情報収集システム、流体制御機器、流体制御機器の動作情報収集方法、及びコンピュータプログラム

Also Published As

Publication number Publication date
US5388445A (en) 1995-02-14
JP2565086B2 (ja) 1996-12-18

Similar Documents

Publication Publication Date Title
JP2565086B2 (ja) パイプラインの監視方法及びその装置
US5272646A (en) Method for locating leaks in a fluid pipeline and apparatus therefore
JP2560978B2 (ja) パイプラインの監視方法及びその監視装置
RU2386992C2 (ru) Датчик технологического параметра с возможностью диагностики
US5932812A (en) Method and devices use in flow measurement
KR101226774B1 (ko) 도압관의 막힘 진단 장치 및 막힘 진단 방법
CN101506629A (zh) 流量测量诊断
WO1997048974A1 (en) Impulse line blockage detector systems and methods
US8453517B2 (en) Pressure guiding tube blockage diagnosing device and blockage diagnosing method
JP2575810B2 (ja) 弁漏洩監視装置
US6813588B1 (en) Control system and method for detecting plugging in differential pressure cells
GB2408100A (en) Acoustic detection of foam formation in pipe flows
JP2765456B2 (ja) パイプラインの漏洩検知方法
EP0300460B1 (en) Apparatus for detecting presence/absence of water leakage from water pipe
CN1109890C (zh) 粒子计数器的估计方法和实施所说的方法的装置
JP3358787B2 (ja) ガス漏洩検知装置、及びこれに用いられるガス漏洩検知方法
JPH07140033A (ja) パイプラインの漏洩検知方法
JP4019419B2 (ja) 導圧管詰まり検出器およびそれを内蔵した差圧・圧力伝送器
JPS6235050B2 (ja)
US5016186A (en) Method of detecting noise disappearance and detecting device therefor
CN220102890U (zh) 一种油气管道泄漏监测系统
JPH01199131A (ja) 漏水検出装置
JPH0774769B2 (ja) ガスパイプラインの漏洩検知方法
JPH0142367B2 (ja)
JPS58167938A (ja) パイプラインの漏洩発生検出装置

Legal Events

Date Code Title Description
FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20081003

Year of fee payment: 12

S531 Written request for registration of change of domicile

Free format text: JAPANESE INTERMEDIATE CODE: R313531

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20081003

Year of fee payment: 12

R350 Written notification of registration of transfer

Free format text: JAPANESE INTERMEDIATE CODE: R350

LAPS Cancellation because of no payment of annual fees