JPS60167007A - プラントの故障原因推定方式 - Google Patents
プラントの故障原因推定方式Info
- Publication number
- JPS60167007A JPS60167007A JP59021694A JP2169484A JPS60167007A JP S60167007 A JPS60167007 A JP S60167007A JP 59021694 A JP59021694 A JP 59021694A JP 2169484 A JP2169484 A JP 2169484A JP S60167007 A JPS60167007 A JP S60167007A
- Authority
- JP
- Japan
- Prior art keywords
- failure
- cause
- matrix
- information
- propagation
- 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
Classifications
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B23/00—Testing or monitoring of control systems or parts thereof
- G05B23/02—Electric testing or monitoring
- G05B23/0205—Electric testing or monitoring by means of a monitoring system capable of detecting and responding to faults
- G05B23/0259—Electric testing or monitoring by means of a monitoring system capable of detecting and responding to faults characterized by the response to fault detection
- G05B23/0275—Fault isolation and identification, e.g. classify fault; estimate cause or root of failure
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Automation & Control Theory (AREA)
- Monitoring And Testing Of Nuclear Reactors (AREA)
- Testing And Monitoring For Control Systems (AREA)
Abstract
(57)【要約】本公報は電子出願前の出願データであるた
め要約のデータは記録されません。
め要約のデータは記録されません。
Description
【発明の詳細な説明】
〔発明の利用分野〕
本発明は、原子カプラント、化学プラントなど複数個の
機器から構成される系において、単一故障(故障原因が
1個)が発生した場合、故障原因となっている機器を、
限られた数のセンサの情報から推定するプラントの故障
原因推定方式に関する。
機器から構成される系において、単一故障(故障原因が
1個)が発生した場合、故障原因となっている機器を、
限られた数のセンサの情報から推定するプラントの故障
原因推定方式に関する。
プラントの故障原因推定方式として、本発明者はすでに
特願昭56−187878号、特願昭57−79954
号、特願昭58−87429号および特願昭58−19
0726号を提案したが、これらの方式では、プラント
を構成する機器の故障波及関係に条件付故障波及関係(
弁の開、閉などの条件が成立して初めて先へ影響が波及
する関係)がある場合には、事前に条件の成立、不成立
を仮定して、モデル化しておく必要があった口上のため
、事前の仮定と異なる状況下では、推定を誤る可能性が
あった。
特願昭56−187878号、特願昭57−79954
号、特願昭58−87429号および特願昭58−19
0726号を提案したが、これらの方式では、プラント
を構成する機器の故障波及関係に条件付故障波及関係(
弁の開、閉などの条件が成立して初めて先へ影響が波及
する関係)がある場合には、事前に条件の成立、不成立
を仮定して、モデル化しておく必要があった口上のため
、事前の仮定と異なる状況下では、推定を誤る可能性が
あった。
本発明の目的は、上記従来技術の問題点を解決するため
に、プランl−を構成する機器の故障波及関係に条件付
故障波及関係がある場合にも、真の故障原因を的確に推
定することのできる故障原因推定方式を提供することに
ある。
に、プランl−を構成する機器の故障波及関係に条件付
故障波及関係がある場合にも、真の故障原因を的確に推
定することのできる故障原因推定方式を提供することに
ある。
上記目的を達成するために、本発明は、複数個の機器か
らなる系において、機器毎には故障率を与え、かつ2機
器間には直接的な故障波及関係(波及方向、波及時間、
波及確率)を与えておき、センサから得られる異常情報
とそれらの検知時刻および正常情報と条件成立情報に基
づいて故障原因を推定している。本発明は特願昭58−
190726号で提案した方式に、条件付故障波及関係
を故障波及関係として扱うか、削除するかを条件成立情
報に基づいて判定する段階を加える事により故障原因推
定の精度を上げた点に特徴がある。
らなる系において、機器毎には故障率を与え、かつ2機
器間には直接的な故障波及関係(波及方向、波及時間、
波及確率)を与えておき、センサから得られる異常情報
とそれらの検知時刻および正常情報と条件成立情報に基
づいて故障原因を推定している。本発明は特願昭58−
190726号で提案した方式に、条件付故障波及関係
を故障波及関係として扱うか、削除するかを条件成立情
報に基づいて判定する段階を加える事により故障原因推
定の精度を上げた点に特徴がある。
以下、本発明の一実施例を第1図から第10図により詳
細に説明する。
細に説明する。
第1図は本発明による故障原因推定方式を実現するプラ
ント系の一実施例の構成を示すものである。
ント系の一実施例の構成を示すものである。
第1図において、プラント101は、複数個の構成機器
102と、その中のいくつかの機器の状態を検出するた
めのセンサ103とからなる。これらのセンサ103は
各構成機器102の動作状態、たとえば、流量、温度、
弁の状態などの信号104を検出し、検出信号105を
各センサ103に対応した異常検知器106からなる異
常検知装置107、または状態判定器109からなる条
件判定装置110に出力する。異常検知装置107には
、あらかじめ、各信号が正常かを判定する基準信号が記
憶されており、基準信号と検出信号105を比較し、各
検出信号105が正常か異常かの信号108を故障原因
推定装置112に出力する。条件判定装置110には、
あらかじめ、各信号に対応する条件が成立か不成立かを
判定する基準信号が記憶されており、基準信号と検出信
号105を比較し、各検出信号105に対応する条件が
成立か不成立かの信号111を故障原因推定装M112
に出力する。故障原因推定装置112では、一定の時間
毎に信号108をチェックし、正常信号から異常信号に
変化した場合にはその時刻を記憶する。オペレータコン
ソール113から演算指示信号114が出力されると。
102と、その中のいくつかの機器の状態を検出するた
めのセンサ103とからなる。これらのセンサ103は
各構成機器102の動作状態、たとえば、流量、温度、
弁の状態などの信号104を検出し、検出信号105を
各センサ103に対応した異常検知器106からなる異
常検知装置107、または状態判定器109からなる条
件判定装置110に出力する。異常検知装置107には
、あらかじめ、各信号が正常かを判定する基準信号が記
憶されており、基準信号と検出信号105を比較し、各
検出信号105が正常か異常かの信号108を故障原因
推定装置112に出力する。条件判定装置110には、
あらかじめ、各信号に対応する条件が成立か不成立かを
判定する基準信号が記憶されており、基準信号と検出信
号105を比較し、各検出信号105に対応する条件が
成立か不成立かの信号111を故障原因推定装M112
に出力する。故障原因推定装置112では、一定の時間
毎に信号108をチェックし、正常信号から異常信号に
変化した場合にはその時刻を記憶する。オペレータコン
ソール113から演算指示信号114が出力されると。
故障原因推定装置112では、演算指示信号114に応
じて、信号108と正常信号から異常信号に変化した時
刻、信号111、および初期データ入力装置115から
入力された故障波及方向。
じて、信号108と正常信号から異常信号に変化した時
刻、信号111、および初期データ入力装置115から
入力された故障波及方向。
波及時間、波及確率、各機器の故障率に対応した信号1
16とに基づいて、故障原因の推定または推定された故
障原因の優先順位付けを行って、その信号117を表示
装置118に出力表示する。
16とに基づいて、故障原因の推定または推定された故
障原因の優先順位付けを行って、その信号117を表示
装置118に出力表示する。
なお、初期データは一度入力しておけば、データの変更
を行わない限り、再入力の必要はない。
を行わない限り、再入力の必要はない。
第2図は、第1図の故障原因推定装置112での処理の
流れの一例を示すフローチャートである6以下、本発明
の実施例を、プラント構成機器の故障波及関係が第3図
のネットワークで表わされるプラントへの適用例を用い
て説明する。
流れの一例を示すフローチャートである6以下、本発明
の実施例を、プラント構成機器の故障波及関係が第3図
のネットワークで表わされるプラントへの適用例を用い
て説明する。
第3図のネットワークで節点はプラント構成機器を表わ
し、矢印の向きは故障の影響が直接波及する方向を表わ
す。四角で囲んだ節点はセンサによって、その機器が異
常がどうかを判定できる機器を表わす。また、三角は条
件を表わす。たとえば第3図では、機器4から3への波
及には条件aの成立が必要であることを示している。
し、矢印の向きは故障の影響が直接波及する方向を表わ
す。四角で囲んだ節点はセンサによって、その機器が異
常がどうかを判定できる機器を表わす。また、三角は条
件を表わす。たとえば第3図では、機器4から3への波
及には条件aの成立が必要であることを示している。
このような複数個(第3図では19個)の機器において
、2機器間で故障の影響が直接波及する場合には、波及
方向と波及時間を与え、この関係を第4図に示すように
行列表現する。この行列をA= (aij)、i、j=
1.=2.−19とすると、たとえば、行列Aの1行2
列目の100は、機器lから2へ故障の影響が直接波及
し、その時間が100秒であることを示す。空白部分は
直接波及がないことを意味し、演算上φとして扱う。行
列Aは、すべての条件が成立している場合の関係を表わ
す。
、2機器間で故障の影響が直接波及する場合には、波及
方向と波及時間を与え、この関係を第4図に示すように
行列表現する。この行列をA= (aij)、i、j=
1.=2.−19とすると、たとえば、行列Aの1行2
列目の100は、機器lから2へ故障の影響が直接波及
し、その時間が100秒であることを示す。空白部分は
直接波及がないことを意味し、演算上φとして扱う。行
列Aは、すべての条件が成立している場合の関係を表わ
す。
さらに、2機器間で故障が直接波及する場合の波及方向
と波及確率を与え、この関係を第5図に示すように行列
表現する。この行列をB = [b i、j)。
と波及確率を与え、この関係を第5図に示すように行列
表現する。この行列をB = [b i、j)。
1、j=l、2.・・・、19とすると、たとえば、行
列Bの1行2列目の1は、機器1から2へ故障の影響が
直接に波及する確率が1であることを示すことになる。
列Bの1行2列目の1は、機器1から2へ故障の影響が
直接に波及する確率が1であることを示すことになる。
空白部分は直接波及がないことを意味し、演算上はOと
して扱う。行列Bは、すべての条件が成立し、でいる場
合の関係を表わす。
して扱う。行列Bは、すべての条件が成立し、でいる場
合の関係を表わす。
第1図の初期データ入力袋W115から故障原因推定袋
W112へ入力する故障波及時間は行列Aの形で、波及
確率は行列Bの形で入力する。
W112へ入力する故障波及時間は行列Aの形で、波及
確率は行列Bの形で入力する。
いま、第3図のネットワークで故障波及関係が表わされ
るプラントにおいて、R常信号を出している機器3,9
,18,19、正常信号を出している機器7,11,1
5.17とする。条件aは不成立、条件すは成立とする
。また、ある基準時刻から異常検知時刻までの経過時間
をそれぞれt3=100秒、tg=200秒*t+a=
300秒*t+e=300秒とし、オペレータコンソー
ル113から演算開始の指示信号が故障原因推定装置1
12に入力されるまでの経過時間t0=300秒とする
。
るプラントにおいて、R常信号を出している機器3,9
,18,19、正常信号を出している機器7,11,1
5.17とする。条件aは不成立、条件すは成立とする
。また、ある基準時刻から異常検知時刻までの経過時間
をそれぞれt3=100秒、tg=200秒*t+a=
300秒*t+e=300秒とし、オペレータコンソー
ル113から演算開始の指示信号が故障原因推定装置1
12に入力されるまでの経過時間t0=300秒とする
。
まず、第2図のステップ119により、正常情報用故障
波及時間行列を計算する。ここでは、行列Bで1未満の
値の要素に対応する行列Aの要素の値と条件付故障波及
関係にある行列Aの要素の値をべてQに置き換える。そ
の結果を第6図に示す。この行列をC= (elj)
t 1 t j= 1 t 2e・・・、19とすると
、たとえば、行列Cの1行2列目の100は、機器1か
ら2へ故障の影響が100秒で必ず直接波及することを
示す。空白部分は必ずしも直接波及しないことを意味し
、演算上はωとして扱う。この行列Cを用いて、すべて
の機器から正常信号を出している機器までの最短故障波
及時間をめる。これは、既知の最短距離と最短路をめる
解法(たとえば、伊理正夫・古林隆:ネットワーク理論
、日科技連出版社(1976) 。
波及時間行列を計算する。ここでは、行列Bで1未満の
値の要素に対応する行列Aの要素の値と条件付故障波及
関係にある行列Aの要素の値をべてQに置き換える。そ
の結果を第6図に示す。この行列をC= (elj)
t 1 t j= 1 t 2e・・・、19とすると
、たとえば、行列Cの1行2列目の100は、機器1か
ら2へ故障の影響が100秒で必ず直接波及することを
示す。空白部分は必ずしも直接波及しないことを意味し
、演算上はωとして扱う。この行列Cを用いて、すべて
の機器から正常信号を出している機器までの最短故障波
及時間をめる。これは、既知の最短距離と最短路をめる
解法(たとえば、伊理正夫・古林隆:ネットワーク理論
、日科技連出版社(1976) 。
PP、34−37.を参照のこと)によりめることがで
きる。得られた波及時間を行列で表わすと第7図のよう
になる。この行列をD=(dij〕、i=1.2.・・
・、19、j=7.10,11,15゜17とすると、
たとえば1行列りの1行1列目の410は、機器1から
7へ故障の影響が波及するのに要する時間は高々410
秒であ゛ることを示す。
きる。得られた波及時間を行列で表わすと第7図のよう
になる。この行列をD=(dij〕、i=1.2.・・
・、19、j=7.10,11,15゜17とすると、
たとえば1行列りの1行1列目の410は、機器1から
7へ故障の影響が波及するのに要する時間は高々410
秒であ゛ることを示す。
空白部分は波及しない可能性があることを意味し、演算
上はωとして扱う。
上はωとして扱う。
次に、第2図のステップ120により、正常情報に基づ
く故障波及関係の更新を行う。ここでは。
く故障波及関係の更新を行う。ここでは。
行列Aについて、正常信号を出している機器に対応する
行と列、および条件が不成立の条件付故障波及関係に対
応する要素をすべてのにする。その結果を第8図に示す
。この行列をE ” (81J) +i、j=1.2.
・・・、19とする。この演算は、第3図のネットワー
ク上では、正常信号を出している機器の入出力技および
条材不成立の枝を削除することに相当する。
行と列、および条件が不成立の条件付故障波及関係に対
応する要素をすべてのにする。その結果を第8図に示す
。この行列をE ” (81J) +i、j=1.2.
・・・、19とする。この演算は、第3図のネットワー
ク上では、正常信号を出している機器の入出力技および
条材不成立の枝を削除することに相当する。
第2図のステップ121から125までは、前述の特願
昭58−190726号と同じであるので、以下のステ
ップは結果のみを示す。
昭58−190726号と同じであるので、以下のステ
ップは結果のみを示す。
第2図のステップ121により、異常情報用故障波及時
間行列と波及確率行列を計算する。第9図には、異常情
報用故障波及時間行列を、第10図には、波及確率行列
を示す。
間行列と波及確率行列を計算する。第9図には、異常情
報用故障波及時間行列を、第10図には、波及確率行列
を示す。
次に、第2図のステップ122により、異常情報に基づ
く故障原因の推定を行う。ここでは、機器1,2.3が
故障原因と推定される。
く故障原因の推定を行う。ここでは、機器1,2.3が
故障原因と推定される。
次に、第2図のステップ123により、故障率に基づく
推定故障原因のチェックを行う。本例では1機器1,2
.3のいずれも故障率は0でないとする。したがって、
推定故障原因はそのまま。
推定故障原因のチェックを行う。本例では1機器1,2
.3のいずれも故障率は0でないとする。したがって、
推定故障原因はそのまま。
次に、第2図のステップ124により、故障波及時間に
基づく推定故障原因のチェックを行う。
基づく推定故障原因のチェックを行う。
この結果、推定故障原因は機器2,3となる。
次に、第2図のステップ125により、推定故障原因の
優先度付けを行う。
優先度付けを行う。
上述した実施例によれば、
(1)複数個の機器から構成されるプラントを対象とし
、設置できるセンサの数が限られていても故障原因の推
定が可能、 (2)推定故障原因が複数個ある場合には、故障原因で
ある可能性の大きいものから順に表示可能。
、設置できるセンサの数が限られていても故障原因の推
定が可能、 (2)推定故障原因が複数個ある場合には、故障原因で
ある可能性の大きいものから順に表示可能。
(3)正常信号を出している機器の故障波及関係および
条件不成立の条件付故障波及関係を除去して、故障波及
時間を計算するため、正確な故障原因の推定が可能。
条件不成立の条件付故障波及関係を除去して、故障波及
時間を計算するため、正確な故障原因の推定が可能。
という効果がある。
以上説明したごとく本発明によれば、
(a)センサおよび異常検知装置が設置されていない機
器を含む系においても故障原因を推定可能、(b)推定
故障原因が複数個ある場合には、各構成機器の故障率、
故障波及確率を用いて、故障原因である可能性の優先度
付けが可能、 (c)正常信号を出している機器の故障波及関係および
条件不成立の条件付故障波及関係を除去して、故障波及
時間を計算するため、正確な故障原因の推定が可能、 となり、この結果、従来方式ではできなかった故障原因
推定が可能となり、故障発生時の対策を容易にできると
いう効果がある。
器を含む系においても故障原因を推定可能、(b)推定
故障原因が複数個ある場合には、各構成機器の故障率、
故障波及確率を用いて、故障原因である可能性の優先度
付けが可能、 (c)正常信号を出している機器の故障波及関係および
条件不成立の条件付故障波及関係を除去して、故障波及
時間を計算するため、正確な故障原因の推定が可能、 となり、この結果、従来方式ではできなかった故障原因
推定が可能となり、故障発生時の対策を容易にできると
いう効果がある。
第1図は本発明の故障原因推定方式によるプラント系の
一実施例の構成図、第2図は第1図の故障原因推定装置
での処理の一例を示すフローチャート第3図はプラント
の故障波及関係と設置センサの位置を示す図、第4図は
直接故障波及時間行列を示す図、第5図は直接故障波及
確率行列を示す図、第6図は正常情報用直接故障波及時
間行列を示す図、第7図は正常情報用故障波及時間行列
を示す図、第8図は異常情報用直接故障波及時間行列を
示す図、第9図は異常情報用故障波及時間行列を示す図
、第1O図は故障波及確率行列を第1図 第 2 巴 第 3 図 第 4 口 第 !5 国 第 6 口 第1頁の続き @発明者橋本 茂男 日立重大みか町5丁目2番1号 株式会社日立製作所大
みか工場内
一実施例の構成図、第2図は第1図の故障原因推定装置
での処理の一例を示すフローチャート第3図はプラント
の故障波及関係と設置センサの位置を示す図、第4図は
直接故障波及時間行列を示す図、第5図は直接故障波及
確率行列を示す図、第6図は正常情報用直接故障波及時
間行列を示す図、第7図は正常情報用故障波及時間行列
を示す図、第8図は異常情報用直接故障波及時間行列を
示す図、第9図は異常情報用故障波及時間行列を示す図
、第1O図は故障波及確率行列を第1図 第 2 巴 第 3 図 第 4 口 第 !5 国 第 6 口 第1頁の続き @発明者橋本 茂男 日立重大みか町5丁目2番1号 株式会社日立製作所大
みか工場内
Claims (1)
- 複数個の機器−からなる系における限られた機器に対応
するセンサの情報により故障原因となる機器を推定する
プラントの故障原因推定方式において、初期入力である
2機器間の故障波及方向、波及時間、波及確率および各
機器の故障率、センサからの入力である異常情報とそれ
らの検知時刻および正常情報と条件成立情報とに基づい
て、故障原因となる機器を、正常情報用故障波及時間行
列を計算するステップと、条件成立情報・正常情報に基
づいて故障波及関係を更新するステップと、異常情報用
故障波及時間行列・波及確率行列を計算するステップと
、異常情報に基づいて故障原因を推定するステップと、
機器の故障率に基づいて推定故障原因をチェックするス
テップと、故障波及時間に基づいて推定故障原因をチェ
ックするステップと、故障率と波及確率とに基づいて推
定故障原因を優先度付けするステップとを実行すること
により推定することを特徴とするプラントの故障原因推
定方式。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP59021694A JPS60167007A (ja) | 1984-02-10 | 1984-02-10 | プラントの故障原因推定方式 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP59021694A JPS60167007A (ja) | 1984-02-10 | 1984-02-10 | プラントの故障原因推定方式 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| JPS60167007A true JPS60167007A (ja) | 1985-08-30 |
Family
ID=12062170
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP59021694A Pending JPS60167007A (ja) | 1984-02-10 | 1984-02-10 | プラントの故障原因推定方式 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JPS60167007A (ja) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP0541326A3 (en) * | 1991-11-04 | 1993-06-09 | General Electric Company | Determination of control system status |
| CN117608274A (zh) * | 2024-01-23 | 2024-02-27 | 四川鼎鸿智电装备科技有限公司 | 一种工程机械故障诊断方法及系统、设备、介质 |
-
1984
- 1984-02-10 JP JP59021694A patent/JPS60167007A/ja active Pending
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP0541326A3 (en) * | 1991-11-04 | 1993-06-09 | General Electric Company | Determination of control system status |
| CN117608274A (zh) * | 2024-01-23 | 2024-02-27 | 四川鼎鸿智电装备科技有限公司 | 一种工程机械故障诊断方法及系统、设备、介质 |
| CN117608274B (zh) * | 2024-01-23 | 2024-04-02 | 四川鼎鸿智电装备科技有限公司 | 一种工程机械故障诊断方法及系统、设备、介质 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US10976729B2 (en) | Method of generating plant normal state prediction data and apparatus using the same | |
| US9268330B2 (en) | Method for fusing data from sensors using a consistency criterion | |
| US20040153815A1 (en) | Methodology for temporal fault event isolation and identification | |
| Fricker Jr | Directionally sensitive multivariate statistical process control methods with application to syndromic surveillance | |
| Li et al. | Historical-orbital-data-based method for monitoring the operational status of satellites in low Earth orbit | |
| JP2890815B2 (ja) | プラントの異常診断装置 | |
| JPS61245210A (ja) | プラントの故障原因推定方式 | |
| EP0697638B1 (fr) | Procédé de détection des séquences complètes et des séquences d'échec dans un système de reconnaissance de situations | |
| JPH04211829A (ja) | 複合型エキスパートシステム | |
| JPS6083108A (ja) | プラントの故障原因推定方式 | |
| JPS60173612A (ja) | 故障源推定方式 | |
| JPS61175709A (ja) | プラントの故障波及範囲予測方式 | |
| JPS59214913A (ja) | プラントの故障源推定方式 | |
| JPS61173313A (ja) | プラントの故障原因推定方式 | |
| JPS60167006A (ja) | プラントの故障波及予測・診断方式 | |
| JPS59214914A (ja) | 故障源推定方式 | |
| JPS60107111A (ja) | プラントの故障波及範囲予測方式 | |
| JPS6126104A (ja) | センサ配置決定装置 | |
| Patton et al. | A review of parity space approaches to fault diagnosis applicable toaerospace systems | |
| Alkaya et al. | Non-linear minimum variance estimation for fault detection systems | |
| Sepulveda | Optimizing a Bank of Kalman Filters for Navigation Integrity for using Parallel Computing and Efficient Software Design | |
| Deyst et al. | Application of likelihood ratio methods to failure detection and identification in the NASA F-8 DFBW aircraft | |
| JP3110060B2 (ja) | 異常事故原因診断方法 | |
| JPS5890122A (ja) | プラントの故障源推定方式 | |
| CN119811638A (zh) | 一种胃食管酸暴露程度分析方法、装置、设备及存储介质 |