▶ 横河電機株式会社の特許一覧 ▶ 国立大学法人東京海洋大学の特許一覧
(19)【発行国】日本国特許庁(JP)
(12)【公報種別】特許公報(B2)
(11)【特許番号】
(24)【登録日】2023-12-12
(45)【発行日】2023-12-20
(54)【発明の名称】推定システム、推定装置及び推定方法
(51)【国際特許分類】
G01N 25/18 20060101AFI20231213BHJP
G01N 25/72 20060101ALI20231213BHJP
【ＦＩ】
G01N25/18 C
G01N25/72 D
【請求項の数】
7
(21)【出願番号】P 2020063026
(22)【出願日】2020-03-31
(65)【公開番号】P2021162434
(43)【公開日】2021-10-11
【審査請求日】2023-03-16
(73)【特許権者】
【識別番号】000006507
【氏名又は名称】横河電機株式会社
(73)【特許権者】
【識別番号】504196300
【氏名又は名称】国立大学法人東京海洋大学
(74)【代理人】
【識別番号】100106909
【弁理士】
【氏名又は名称】棚井 澄雄
(74)【代理人】
【識別番号】100146835
【弁理士】
【氏名又は名称】佐伯 義文
(74)【代理人】
【識別番号】100167553
【弁理士】
【氏名又は名称】高橋 久典
(74)【代理人】
【識別番号】100181124
【弁理士】
【氏名又は名称】沖田 壮男
(72)【発明者】
【氏名】伊藤 俊平
(72)【発明者】
【氏名】田中 仁章
(72)【発明者】
【氏名】波津久 達也
(72)【発明者】
【氏名】井原 智則
(72)【発明者】
【氏名】盛田 元彰
【審査官】野田 華代
(56)【参考文献】
【文献】国際公開第２０１９／２０２９８１（ＷＯ，Ａ１）
【文献】特開平８－１５１８９（ＪＰ，Ａ）
【文献】特表平７－５００６７３（ＪＰ，Ａ）
【文献】特開２０００－８８５６１（ＪＰ，Ａ）
【文献】特開昭４８－２９４８４（ＪＰ，Ａ）
【文献】米国特許出願公開第２０１１／３０８５４８（ＵＳ，Ａ１）
【文献】特開昭５４－１０７６２（ＪＰ，Ａ）
【文献】欧州特許出願公開第４１６３５８７（ＥＰ，Ａ１）
(58)【調査した分野】(Int.Cl.，ＤＢ名)
Ｇ０１Ｎ ２５／００－２５／７２
(57)【特許請求の範囲】
【請求項1】
流体が流れる配管の外表面温度の測定値を取得する取得部と、前記配管の内表面における析出物の除去処理の実行時刻を推測する推測部と、
前記析出物の厚さの推定に用いられる推定値を、前記析出物の除去処理の実行時刻の推測結果から所定時間内に得られた前記配管の外表面温度の測定値を用いて補正する補正部と、
補正された前記推定値と前記配管の外表面温度の測定値とに基づいて前記析出物の厚さを推定する推定部と、
を備える推定システム。
【請求項2】
前記補正部は、前記析出物の厚さの推定結果が予め規定された閾値以下になるように前記推定値を補正する、請求項１に記載の推定システム。
【請求項3】
前記補正部は、前記析出物の熱伝導率の前記推定値と、前記配管の内側における流体の温度の前記推定値と、前記配管の外側における流体の温度の前記推定値と、前記配管の内側における熱伝達率の前記推定値と、前記配管の外側における熱伝達率の前記推定値とのうちの少なくとも一つを補正する、請求項１に記載の推定システム。
【請求項4】
前記取得部は、更に、前記配管の外表面における熱流束の測定値と、前記配管の外表面を覆う断熱材の内部温度の測定値とのうちの少なくとも一つを取得し、前記推定部は、前記配管の外表面における熱流束の測定値と前記断熱材の内部温度の測定値とのうちの少なくとも一つと、前記配管の外表面温度の測定値と、補正された前記推定値とに基づいて、前記析出物の厚さを推定する、
請求項１又は請求項２に記載の推定システム。
【請求項5】
前記推測部は、前記析出物の除去処理を検出または予測した時刻から前記析出物の除去処理の実行時刻を推測し、前記補正部は、前記推定値を、前記析出物の除去処理の実行時刻の推測結果から所定時間内に得られた前記配管の外表面温度の測定値を用いて補正する、
請求項１から請求項３のいずれか一項に記載の推定システム。
【請求項6】
流体が流れる配管の内表面における析出物の厚さの推定に用いられる推定値を、前記析出物の除去処理の実行時刻の推測結果から所定時間内に得られた前記配管の外表面温度の測定値を用いて補正する補正部と、補正された前記推定値と前記配管の外表面温度の測定値とに基づいて前記析出物の厚さを推定する推定部と、
を備える推定装置。
【請求項7】
推定システムが実行する推定方法であって、流体が流れる配管の外表面温度の測定値を取得する取得ステップと、
前記配管の内表面における析出物の除去処理の実行時刻を推測する推測ステップと、
前記析出物の厚さの推定に用いられる推定値を、前記析出物の除去処理の実行時刻の推測結果から所定時間内に得られた前記配管の外表面温度の測定値を用いて補正する補正ステップと、
補正された前記推定値と前記配管の外表面温度の測定値とに基づいて前記析出物の厚さを推定する推定ステップと、
を含む推定方法。
【発明の詳細な説明】
【技術分野】
【0001】
本発明は、推定システム、推定装置及び推定方法に関する。
【背景技術】
【0002】
流体が流れる配管の内表面（流路壁）に、その流体に含まれている物質の析出物が付着することが、従来から知られている。例えば、オイル及びガスのパイプラインでは、温度及び圧力の条件に応じて、ハイドレート、ワックス、アスファルテン又はスケールといった析出物が付着することが知られている。このような析出物は、例えば、ピグと呼ばれる器具をパイプラインの内側に通過させることにより、或いは、析出物の発生を阻害する阻害剤をパイプラインに投入することにより、パイプラインの内表面から除去される。従来では、配管の内表面に付着した析出物の厚さが推定され、その厚さが閾値を超えた場合に、上述した除去処理が実行されていた。以下の特許文献１～４には、配管の内表面に付着した析出物の厚さの推定方法が、それぞれ開示されている。
【先行技術文献】
【特許文献】
【0003】
【文献】米国特許第６６４４８４８号明細書
【文献】米国特許第８９６０３０５号明細書
【文献】特開２０１６－１６６７８１号公報
【文献】米国特許出願公開第２００４／００５９５０５号明細書
【発明の概要】
【発明が解決しようとする課題】
【0004】
ところで、配管の内表面に付着した析出物の厚さを推定するために、配管の外表面温度、配管を流れる流体の温度、配管の外側及び内側における熱伝達率、及び、析出物の熱伝導率等の各物理量が必要になる場合がある。これらの物理量の中には、パイプライン等の現場において測定することができるものと、測定することができないものとがある。そこで、析出物の厚さを推定するために、測定することができない物理量の代わりに推定値が用いられる場合がある。しかしながら、推定値と真値との乖離が大きい場合には、配管の内表面における析出物の厚さの推定誤差を所定範囲に収めることができないという問題がある。
【0005】
また、特にパイプラインなど自然環境下に設置される配管の場合、周囲の環境変化（配管が埋もれるなど）が起こりやすく、推定値と真値との乖離が大きくなりやすいという問題がある。
【0006】
本発明は上記事情に鑑みてなされたものであり、配管の内表面における析出物の厚さの推定誤差を所定範囲に収めることが可能である推定システム、推定装置及び推定方法を提供することを目的としている。
【課題を解決するための手段】
【0007】
上記課題を解決するために、本発明の一態様による推定システムは、流体が流れる配管（２００）の外表面温度の測定値を取得する取得部（４０）と、前記配管（２００）の内表面における析出物の除去処理の実行時刻を推測する推測部（３）と、前記析出物の厚さの推定に用いられる推定値を、前記析出物の除去処理の実行時刻の推測結果から所定時間内に得られた前記配管の外表面温度の測定値を用いて補正する補正部（４２）と、補正された前記推定値と前記配管（２００）の外表面温度の測定値とに基づいて前記析出物の厚さを推定する推定部（４３）と、を備える。
【0008】
また、本発明の一態様による推定システムは、前記補正部（４２）が、前記析出物の厚さの推定結果が予め規定された閾値以下になるように前記推定値を補正する。
【0009】
また、本発明の一態様による推定システムは、前記補正部（４２）が、前記析出物の熱伝導率の前記推定値と、前記配管（２００）の内側における流体の温度の前記推定値と、前記配管（２００）の外側における流体の温度の前記推定値と、前記配管（２００）の内側における熱伝達率の前記推定値と、前記配管（２００）の外側における熱伝達率の前記推定値とのうちの少なくとも一つを補正する。
【0010】
また、本発明の一態様による推定システムは、前記取得部（４０）が、更に、前記配管（２００）の外表面における熱流束の測定値と、前記配管（２００）の外表面を覆う断熱材の内部温度の測定値とのうちの少なくとも一つを取得し、前記推定部（４３）が、前記配管（２００）の外表面における熱流束の測定値と前記断熱材の内部温度の測定値とのうちの少なくとも一つと、前記配管（２００）の外表面温度の測定値と、補正された前記推定値とに基づいて、前記析出物の厚さを推定する。
【0011】
また、本発明の一態様による推定システムは、前記推測部は、前記析出物の除去処理を検出または予測した時刻から前記析出物の除去処理の実行時刻を推測し、前記補正部（４２）が、前記推定値を、前記析出物の除去処理の実行時刻の推測結果から所定時間内に得られた前記配管の外表面温度の測定値を用いて補正する。
【0012】
本発明の一態様による推定装置（４ａ）は、流体が流れる配管（２００）の内表面における析出物の厚さの推定に用いられる推定値を、前記析出物の除去処理の実行時刻の推測結果から所定時間内に得られた前記配管の外表面温度の測定値を用いて補正する補正部（４２）と、補正された前記推定値と前記配管（２００）の外表面温度の測定値とに基づいて前記析出物の厚さを推定する推定部（４３）と、を備える。
【0013】
本発明の一態様による推定方法は、推定システムが実行する推定方法であって、流体が流れる配管（２００）の外表面温度の測定値を取得する取得ステップと、前記配管（２００）の内表面における析出物の除去処理の実行時刻を推測する推測ステップと、前記析出物の厚さの推定に用いられる推定値を、前記析出物の除去処理の実行時刻の推測結果から所定時間内に得られた前記配管の外表面温度の測定値を用いて補正する補正ステップと、補正された前記推定値と前記配管（２００）の外表面温度の測定値とに基づいて前記析出物の厚さを推定する推定ステップと、を含む。
【発明の効果】
【0014】
本発明によれば、配管の内表面における析出物の厚さの推定誤差を所定範囲に収めることが可能である。
【図面の簡単な説明】
【0015】
【図1】第１実施形態における、推定システムの構成例を示す図である。
【図2】第１実施形態における、センサの配置例を示す断面図である。
【図3】第１実施形態における、内表面に析出物が付着している配管の例を示す断面図である。
【図4】第１実施形態における、配管の外表面温度の例を示す図である。
【図5】第１実施形態における、厚さの推定結果の例を示す図である。
【図6】第１実施形態における、推定システムの動作例を示すフローチャートである。
【図7】第１実施形態の第２変形例における、センサの配置例を示す断面図である。
【図8】第２実施形態における、推定システムの構成例を示す図である。
【発明を実施するための形態】
【0016】
以下、図面を参照して本発明の実施形態による推定システム、推定装置及び推定方法について詳細に説明する。以下では、まず本発明の実施形態の概要について説明し、続いて本発明の各実施形態の詳細について説明する。
【0017】
〔概要〕
本発明の実施形態は、配管の内表面における析出物の厚さの推定誤差を所定範囲に収めることができるようにするものである。例えば、オイル及びガスのパイプラインの内表面に付着するハイドレート、ワックス、アスファルテン又はスケールといった析出物の厚さの推定誤差を、所定範囲に収めることができるようにするものである。
本発明の実施形態は、配管の内表面における析出物の厚さの推定誤差を所定範囲に収めることができるようにするものである。例えば、オイル及びガスのパイプラインの内表面に付着するハイドレート、ワックス、アスファルテン又はスケールといった析出物の厚さの推定誤差を、所定範囲に収めることができるようにするものである。
【0018】
特許文献３に開示された装置は、配管の外表面に設置された４個の温度センサを用いて、析出物の厚さを推定する。しかしながら、配管の外側及び内側における各流体の温度と、流体及び配管の間の熱の伝わり方を表すパラメータ（例えば、熱伝達率）とが精度よく推定又は測定されていなければ、配管の内表面における析出物の厚さの推定誤差を所定範囲に収めることができない。すなわち、配管の外側及び内側における各流体の温度の推定誤差が大きい場合には、流体の温度の推定値に基づいて推定された厚さの推定結果を所定範囲に収めることができない。
【0019】
また、配管の外側及び内側における流体の温度が時間経過に応じて変化するので、配管の外側及び内側における流体の圧力、流速及び流動は時間経過に応じて変化する。更に、配管の一部が海底に埋もれる場合がある。これらの要因により、過去の測定値に基づいて推定された推定値が推定誤差を有する場合がある。推定値が推定誤差を有する場合、推定誤差を有する推定値に基づいて推定された厚さの推定結果を所定範囲に収めることができない。
【0020】
これらのように特許文献１～４では、析出物の厚さの推定に用いられる推定値が補正されることがないので、配管の内表面における析出物の厚さの推定誤差は蓄積されたままである。
【0021】
本発明の実施形態では、推定システムは、配管の内表面における析出物の厚さの推定に用いられる推定値（熱等の物理量に関する推定値）を、析出物の除去処理の実行時刻（例えば、析出物を除去する器具が配管の内側を通過した通過時刻、あるいは、析出物を除去する器具が配管の内側を通過するであろう又は通過したであろう予測時刻）から所定時間内に得られた配管の外表面温度の測定値を用いて補正する。ここで、析出物の除去処理の実行時刻では析出物の厚さが閾値以下（例えば、０）になっているので、推定システムは、予め規定された閾値以下に厚さの推定結果がなるように推定値を補正する。推定値が補正されることによって、析出物の厚さの推定誤差は蓄積されなくなる。その結果、配管の内表面における析出物の厚さの推定誤差を所定範囲に収めることが可能である。なお、この所定時間は、析出物が除去された後に、センサの測定値に所定の影響を与える程度に再び付着し始めるまでの時間よりも短いことが望ましい。このような所定時間は、例えば析出物が付着するまでに要する時間を予め実験などで得ておくことで定められてもよい。
【0022】
〔第１実施形態〕
〈推定システム〉
図１は、第１実施形態における、推定システムの構成例を示す図である。本実施形態の推定システム１ａは、配管の外表面温度の測定値と、熱等の物理量に関する推定値とに基づいて、配管の内表面における析出物の厚さを推定するシステムである。
〈推定システム〉
図１は、第１実施形態における、推定システムの構成例を示す図である。本実施形態の推定システム１ａは、配管の外表面温度の測定値と、熱等の物理量に関する推定値とに基づいて、配管の内表面における析出物の厚さを推定するシステムである。
【0023】
推定システム１ａは、複数のセンサ２と、検出部（推測部）３と、推定装置４ａとを備える。推定装置４ａは、通信回線１００を経由してセンサ２及び検出部３と通信可能である。通信回線１００は、インターネット、ＷＡＮ（Wide Area Network）等の通信回線である。通信回線１００は、有線回線又は無線回線のいずれでもよく、有線回線及び無線回線の両方でもよい。
【0024】
配管２００は、流体を送る配管（流体が流れる配管）である。以下、配管２００を用いて送られる流体は、一例としてオイル及びガスである。オイル及びガスは、坑井から生産される。生産されたオイル及びガスは、配管２００の内側に流される。オイル及びガスの温度は、オイル及びガスが坑井から生産された直後において測定される。
【0025】
配管２００は例えば地上又は海中に設置されており、配管２００の一部又は全部は地中又は海底に埋設されていてもよい。配管２００の内表面には、析出物２０１が付着する。析出物２０１は、配管２００を用いて送られるオイル又はガス等に含まれている物質であり、例えば、ハイドレート、ワックス若しくはアスファルテン等の有機物、又はスケール等の無機物である。
【0026】
器具３００は、析出物２０１を除去する器具であり、例えば、ピグである。器具３００は、配管２００の内側を管軸方向に通過するという除去処理を実行することによって、析出物２０１を配管２００の内表面から除去する。
【0027】
センサ２は、温度センサであり、例えば、熱電対、測温抵抗体、ＤＴＳ（Distributed Temperature Sensor）又はサーモカメラである。第１実施形態においては、配管２００の外表面の温度（温度分布）を測定するために、複数のセンサ２は、配管２００の外周方向に所定の間隔で、配管２００の外表面に設置される。ここでセンサ２がサーモカメラであって、１台のサーモカメラで配管２００の外表面における所望の範囲の温度（温度分布）を測定することができる場合には、センサ２（サーモカメラ）は１台であってもよい。センサ２は、設置された位置において、配管２００の外表面温度を測定する。また、センサ２は、配管２００の外表面温度を測定できれば、配管外表面に限らず、配管に埋め込んで設置してもよい。センサ２は、配管２００の外表面温度の測定値を、通信回線１００を経由して推定装置４ａに出力する。
【0028】
図２は、第１実施形態における、センサの配置例を示す断面図である。なお、図２に示す断面図は、配管２００の半径方向の断面図である。図２（ａ），（ｂ）に示す例は何れも、４個のセンサ２（センサ２－１～２－４）が、配管２００の外表面に、配管２００の外周方向に沿って等間隔（９０度の間隔）で設置されている。これによって、センサ２－１～２－４は、配管２００の外周方向における配管２００の外表面の温度分布を測定することが可能である。
【0029】
ここで、図２（ｂ）に示す例では、配管２００の外表面に、センサ２－１～２－４とともに、配管２００の外表面を覆う断熱材２３０が設けられているが、図２（ａ）に示す例では、配管２００の外表面に断熱材２３０は設けられていない。断熱材２３０は、配管２００の外側からセンサ２への熱の影響を低減するために設けられる。断熱材２３０を設けることで、例えばセンサ２－１～２－４に対する海水の熱の影響を低減することができる。
【0030】
検出部３は、器具３００が所定の位置（例えば、検出部３の設置位置）を通過したことを検出する。例えば、検出部３は、配管２００の内側を管軸方向に移動する器具３００が発する振動や磁界を検出し、その検出結果に基づいて、器具３００が検出部３の設置位置を通過したこと（通過時刻）を検出（推測）する。或いは検出部３は、超音波若しくは放射線の反射又は透過を検出し、その検出結果に基づいて、検出部３の設置位置を器具３００が通過したことを検出（推測）する。
【0031】
検出部３は、器具３００が通過したことを検出（推測）することによって、その検出結果（検出時刻）を、配管２００の内表面における析出物２０１の除去処理の実行時刻（タイミング）として推測する。検出部３は、析出物２０１の除去処理の実行時刻の推測結果を、通信回線１００を経由して推定装置４ａに出力する。
【0032】
〈推定装置〉
推定装置４ａは、取得部４０と、記憶部４１と、補正部４２と、推定部４３と、判定部４４と、出力部４５とを備える。取得部４０は、複数のセンサ２が出力した測定値（配管２００の外表面温度の測定値）、検出部３が出力した推測結果（析出物２０１の除去処理の実行時刻の推測結果）を取得する。取得部４０は、上記の測定値及び推測結果を、通信回線１００を経由して取得する。取得部４０は、取得した測定値及び推測結果を記憶部４１に記憶させる。判定部４４は、器具３００の通過を検出部３が検出（推測）したか否かを判定する。すなわち、判定部４４は、析出物２０１の除去処理が実行されたか否かを判定する。
推定装置４ａは、取得部４０と、記憶部４１と、補正部４２と、推定部４３と、判定部４４と、出力部４５とを備える。取得部４０は、複数のセンサ２が出力した測定値（配管２００の外表面温度の測定値）、検出部３が出力した推測結果（析出物２０１の除去処理の実行時刻の推測結果）を取得する。取得部４０は、上記の測定値及び推測結果を、通信回線１００を経由して取得する。取得部４０は、取得した測定値及び推測結果を記憶部４１に記憶させる。判定部４４は、器具３００の通過を検出部３が検出（推測）したか否かを判定する。すなわち、判定部４４は、析出物２０１の除去処理が実行されたか否かを判定する。
【0033】
記憶部４１は、不揮発性の記録媒体（非一時的な記録媒体）を有しており、例えば、プログラム、モデル、パラメータ等を記憶する。モデルは、外表面温度の測定値と析出物２０１の厚さとの関係を表すものである。このモデルは、例えば関数又は変換テーブルの形式で表される。モデルが関数の形式で表される場合、パラメータは、モデル「関数ｆ（ｐ）」の引数「ｐ」である。
【0034】
このようなパラメータとしては、例えば、センサ２－１～２－４のセンサ２ごとの外表面温度の測定値、析出物２０１の熱伝導率の推定値、配管２００の内側における熱伝達率の推定値、配管２００の外側における熱伝達率の推定値が挙げられる。各パラメータは、例えば、過去の測定値と、経験値と、フローシミュレータによる推定値の解析結果とのうちの少なくとも一つに基づいて予め決定される。推定値の解析処理と測定とは、リアルタイムで繰り返し実行されてもよい。
【0035】
仮に、析出物２０１の厚さの推定に用いられる推定値が真値に対して大きな誤差を有している場合、器具３００が通過した直後において、実際には析出物２０１が除去されているにもかかわらず、析出物２０１の厚さの推定結果が０にならない（例えば、析出物２０１の厚さの推定結果が、０から大きく外れた値となる）。
【0036】
そこで、補正部４２は、析出物２０１の厚さの推定に用いられる推定値を、析出物２０１の除去処理の実行時刻から所定時間内に得られた測定値を用いて補正する。所定時間は、例えば、器具３００の通過に起因する配管２００の外表面温度（センサ２の測定値）の変化が収束し、配管２００の外表面温度（センサ２の測定値）が安定するのに要する時間であり、例えば、１～数分間、１～数時間、又は、１～数日間とすることができる。補正部４２は、析出物２０１の厚さの推定結果が閾値以下（例えば、０）になるように、析出物２０１の厚さの推定に用いられる推定値（複数の推定値が用いられる場合には、それら複数の推定値のうち少なくとも一つの推定値）を補正する。なお、閾値は予め規定される。
【0037】
析出物２０１の厚さの推定に用いられる推定値とは、配管２００周辺の環境における熱等の物理量に関する推定値であって、例えば、析出物２０１の熱伝導率の推定値と、配管２００の内側におけるオイル及びガスの温度の推定値と、配管２００の外側における流体（例えば、海水）の温度の推定値と、配管２００の内側（例えば、析出物、オイル及びガス）における熱伝達率の推定値と、配管２００の外側（例えば、海水、海底）における熱伝達率の推定値とのうちの少なくとも一つである。
【0038】
推定部４３は、補正された推定値と配管２００の外表面温度の測定値とに基づいて、析出物２０１の厚さを推定する。配管の外表面温度の測定値に基づいて析出物の厚さを推定する方法は、国際公開第２０１９／２０２９８１号（以下「参考文献」という。）に開示されている。なお、配管２００の外表面温度の測定値に基づいて推定部４３が析出物２０１の厚さを推定する方法は、上記の参考文献に開示された方法に限定されるものではない。
【0039】
推定部４３は、析出物２０１の厚さの推定結果を出力部４５に出力する。推定部４３は、析出物２０１の厚さの推定結果を記憶部４１に記録してもよく、また、析出物２０１の厚さの推定結果を記憶部４１に記録すると共に出力部４５に出力してもよい。推定部４３は、推定処理を終了するか否かを表す指示信号を、外部端末（不図示）から取得部４０を経由して取得してもよい。なお、推定部４３で行われる処理の詳細は、後述する。
【0040】
出力部４５は、データを外部装置（不図示）に出力する機能部である。例えば、出力部４５は、析出物２０１の厚さの推定結果を、外部装置に出力する。外部装置は、例えば、液晶ディスプレイ又は有機ＥＬディスプレイ等の表示デバイスである。表示デバイスは、析出物２０１の厚さの推定結果を表す画像（例えば、グラフ、断面図）を表示する。オペレータは、表示デバイスに表示された厚さの推定結果を確認することができる。
【0041】
上述した取得部４０、補正部４２、推定部４３、判定部４４及び出力部４５の一部又は全部は、ＣＰＵ（Central Processing Unit）等のプロセッサが、記憶部４１に記憶されたプログラムを実行することにより、ソフトウェア的に実現される。つまり、取得部４０、補正部４２、推定部４３、判定部４４及び出力部４５の一部又は全部の機能は、ソフトウェアとハードウェア資源とが協働することによって実現される。
【0042】
上記のプログラムは、コンピュータ読み取り可能な記録媒体に記録されてもよい。コンピュータ読み取り可能な記録媒体とは、例えばフレキシブルディスク、光磁気ディスク、ＲＯＭ（Read Only Memory）、ＣＤ－ＲＯＭ（Compact Disc Read Only Memory）等の可搬媒体、コンピュータシステムに内蔵されるハードディスク等の記憶装置などの非一時的な記録媒体である。プログラムは、通信回線１００を経由して、取得部４０によって受信されてもよい。
【0043】
なお、上述した取得部４０、補正部４２、推定部４３、判定部４４及び出力部４５の一部又は全部は、ハードウェアによって実現されてもよい。例えば、ＬＳＩ（Large Scale Integration circuit）、ＡＳＩＣ（Application Specific Integrated Circuit）、ＰＬＤ（Programmable Logic Device）又はＦＰＧＡ（Field Programmable Gate Array）等を用いた電子回路（electronic circuit又はcircuitry）によって実現されても良い。
【0044】
〈推定部で行われる処理〉
推定部４３は、少なくとも１個のセンサ２の外表面温度の測定値に基づいて、所定の位置（例えば、そのセンサ２の設置位置）における析出物２０１の厚さを推定する。例えば、推定部４３は、センサ２－１～２－４の各設置位置における外表面温度の測定値（温度分布）に基づいて、析出物２０１の形状（厚さ分布）を式（１）のように推定する。
推定部４３は、少なくとも１個のセンサ２の外表面温度の測定値に基づいて、所定の位置（例えば、そのセンサ２の設置位置）における析出物２０１の厚さを推定する。例えば、推定部４３は、センサ２－１～２－４の各設置位置における外表面温度の測定値（温度分布）に基づいて、析出物２０１の形状（厚さ分布）を式（１）のように推定する。
【0045】
ｔｈ＝ｆ（ｐ１，ｐ２，ｐ３，ｐ４） …（１）
【0046】
ここで、「ｔｈ」は、少なくとも１個のセンサ２の設置位置における析出物２０１の厚さの推定結果を表す。「ｆ」は、予め規定された関数を表す。パラメータ「ｐ１」は、少なくとも１個のセンサ２の外表面温度の測定値を表す。パラメータ「ｐ２」は、析出物２０１の熱伝導率の推定値を表す。パラメータ「ｐ３」は、配管２００の内側における流体（オイル及びガス）の温度の推定値又は測定値を表す。オイル及びガスの温度は、例えば、オイル及びガスが坑井から生産された直後において測定された温度に基づいて推定される。パラメータ「ｐ４」は、配管２００の外側における流体の温度の推定値又は測定値を表す。配管２００が地上に設置されている場合、配管２００の外側における流体は大気である。配管２００が海底に設置されている場合、配管２００の外側における流体は海水である。
【0047】
また、配管２００の内側に混相流が形成されている場合（熱の伝わり方に偏りがある場合）でも析出物２０１の厚さを推定できるように、記憶部４１は、配管２００の内側における流体の熱伝達率の推定値を、変数として予め記憶してもよい。推定部４３は、配管２００の内側における熱伝達率の推定値と、配管２００の外側における熱伝達率の推定値とに基づいて、析出物２０１の厚さを式（２）のように推定してもよい。
【0048】
ｔｈ＝ｆ（ｐ１，ｐ２，ｐ３，ｐ４，ｐ５） …（２）
【0049】
ここで、パラメータ「ｐ５」は、配管２００の内側における熱伝達率の推定値を表す。これにより推定部４３は、配管２００の内側に混相流が形成されている場合でも、配管２００の内側における熱の伝わり方を考慮して、析出物２０１の厚さを推定することができる。
【0050】
また、配管２００の一部が海底に埋設されている又は海流が激しい場合でも析出物２０１の厚さを推定できるように、記憶部４１は、配管２００の外側における熱伝達率の推定値を、変数として記憶してもよい。推定部４３は、配管２００の内側における熱伝達率の推定値と、配管２００の外側における熱伝達率の推定値とに基づいて、析出物２０１の厚さを式（３）のように推定してもよい。
【0051】
ｔｈ＝ｆ（ｐ１，ｐ２，ｐ３，ｐ４，ｐ５，ｐ６） …（３）
【0052】
ここで、パラメータ「ｐ６」は、配管２００の外側における熱伝達率の推定値を表す。これにより推定部４３は、配管２００の一部が例えば海底に埋設されている場合でも、配管２００の外側における熱の伝わり方を考慮して、析出物２０１の厚さを推定することができる。
【0053】
推定部４３は、配管２００の外側における熱伝達率の推定値に基づいて、析出物２０１の厚さを式（４）のように推定してもよい。
【0054】
ｔｈ＝ｆ（ｐ１，ｐ２，ｐ３，ｐ４，ｐ６） …（４）
【0055】
これにより推定部４３は、配管２００の内側における熱伝達率を推定せずに、配管２００の外側における熱の伝わり方を考慮して、析出物２０１の厚さを推定することができる。
【0056】
〈配管に付着した析出物の例〉
図３は、第１実施形態における、内表面に析出物が付着している配管の例を示す断面図である。なお、図３に示す断面図は、図２と同様に、配管２００の半径方向の断面図である。図３に例示する析出物２０１では、センサ２－１の設置位置における厚さが「ｔｈ１」であり、センサ２－２の設置位置における厚さが「ｔｈ２」であり、センサ２－３の設置位置における厚さが「ｔｈ３」であり、センサ２－４の設置位置における厚さが「ｔｈ４」である。
図３は、第１実施形態における、内表面に析出物が付着している配管の例を示す断面図である。なお、図３に示す断面図は、図２と同様に、配管２００の半径方向の断面図である。図３に例示する析出物２０１では、センサ２－１の設置位置における厚さが「ｔｈ１」であり、センサ２－２の設置位置における厚さが「ｔｈ２」であり、センサ２－３の設置位置における厚さが「ｔｈ３」であり、センサ２－４の設置位置における厚さが「ｔｈ４」である。
【0057】
図３では、配管２００の内側における一例として、オイル２１０及びガス２２０の混相流が形成されている。すなわち、配管２００の内側の天井付近にガス２２０が集まっている。配管２００の内側に混相流が形成されている場合、配管２００の内側において熱の伝わり方に偏りがあるので、配管２００の外周方向の位置に応じて熱の伝わり方が変化する。このようなことに起因して、厚さ「ｔｈ１」～「ｔｈ４」のように、析出物２０１の厚さが不均一になることがある。図３では、一例として、厚さ「ｔｈ４」が最も厚く、厚さ「ｔｈ１」が最も薄い。
【0058】
なお、配管２００の一部が例えば海底に埋設されている場合、埋設されていることの影響によっても、配管２００の外周方向の位置に応じて熱の伝わり方が変化する。また、配管２００の外側において海流が激しい場合、海流の影響によっても、配管２００の外周方向の位置に応じて熱の伝わり方が変化する。
【0059】
図４は、第１実施形態における、配管の外表面温度の例を示す図である。横軸は、配管２００の外周方向の位置（センサ２－１～２－４の各設置位置）を示す。縦軸は、図３に示された配管２００の外表面温度の測定値「ｔｐ」を示す。図４に示す例では、センサ２－１による測定値が「ｔｐ１」であり、センサ２－２による測定値が「ｔｐ２」であり、センサ２－３による測定値が「ｔｐ３」であり、センサ２－４による測定値が「ｔｐ４」である。図４に示す例では、センサ２－２による測定値「ｔｐ２」が最も低い。これは、図３に示された配管２００の内側の天井付近に集まっているガス２２０の熱伝導率が、オイル２１０の熱伝導率よりも低いためである。また、図４に示す例では、センサ２－４による測定値と、センサ２－３による測定値と、センサ２－２による測定値との順に、配管２００の外表面温度が高い。
【0060】
図５は、第１実施形態における、厚さの推定結果の例を示す図である。横軸は、配管２００の外周方向の位置（センサ２－１～２－４の各設置位置）を示す。縦軸は、図３に例示された配管２００における析出物２０１の厚さの推定結果を示す。なお、図５に示す測定結果は、推定部４３で、上述した処理が行われて得られたものである。
【0061】
推定部４３は、図４の縦軸に例示された外表面温度の測定値を用いて、センサ２－１～２－４のセンサ２ごとに析出物２０１の厚さを推定する。例えば、推定部４３は、外表面温度の測定値と式（１）～（４）のいずれかとを用いて、センサ２ごとに析出物２０１の厚さを推定する。なお、図５の縦軸に示された厚さの推定結果は一例である。
【0062】
〈推定システムの動作〉
図６は、第１実施形態における、推定システムの動作例を示すフローチャートである。図６に示すフローチャートの処理（推定値の補正フロー）は、通常の動作フロー（センサ２による検出結果（温度）に基づく析出物２０１の厚さを推定するフロー）中のいずれのタイミングであっても、優先して行われ得る。
図６は、第１実施形態における、推定システムの動作例を示すフローチャートである。図６に示すフローチャートの処理（推定値の補正フロー）は、通常の動作フロー（センサ２による検出結果（温度）に基づく析出物２０１の厚さを推定するフロー）中のいずれのタイミングであっても、優先して行われ得る。
【0063】
図６に示すフローチャートの処理が開始されると、まず、判定部４４は、器具３００の通過を検出部３が検出（推測）できたか否かを判定する。すなわち、判定部４４は、析出物２０１の除去処理が実行されたか否かを判定する（ステップＳ１０１）。例えば、判定部４４は、析出物２０１の除去処理が実行されたことを示す通知を受信したか否かに基づいて、析出物２０１の除去処理が実行されたか否かを判定する。
【0064】
器具３００の通過を検出部３が検出（推測）できていないと判定部４４が判定した場合（ステップＳ１０１：ＮＯ）、判定部４４は、配管２００内の析出物２０１の除去処理がなされていないと判定する（ステップＳ１０５）。推定システム１ａは、図６に示された推定処理を終了する。
【0065】
これに対し、器具３００の通過を検出部３が検出（推測）できたと判定部４４が判定した場合（ステップＳ１０１：ＹＥＳ）、判定部４４は、配管２００内の析出物２０１の除去処理がなされたと判定する（ステップＳ１０２）。器具３００の通過時刻から所定時間内にセンサ２は外表面温度を測定する（ステップＳ１０３）。このように、析出物２０１の厚さを推定するための外表面温度の測定とは別に、推定値の補正のためにセンサ２は外表面温度を測定する。もちろん、これに限定されず、例えば、除去処理の実行時刻が推測できた後の初回の外表面温度の測定結果（析出物２０１の厚さを推定のための測定結果）を用いて、推定値の補正が実行されてもよい。
【0066】
続いて、補正部４２は、析出物２０１の厚さの推定に用いられる各推定値を、除去処理後の外表面温度の測定値を用いて補正する（ステップＳ１０４）。具体的には、補正部４２は、析出物２０１の厚さの推定結果が閾値以下（例えば、０）になるように、析出物２０１の厚さの推定に用いられる推定値（複数の推定値を用いる場合には、少なくとも一つの推定値）を補正する。
【0067】
以上説明した通り、本実施形態では、補正部４２が、析出物２０１の厚さの推定に用いられる少なくとも一つの推定値を、除去処理の実行時刻としての検出部３の検出結果（検出時刻）から所定時間内に得られた配管２００の外表面温度の測定値を用いて補正するようにしている。ここで、補正部４２は、析出物２０１の厚さの推定結果が予め規定された閾値以下（例えば、０）になるように推定値を補正してもよい。推定値が補正されることによって、析出物２０１の厚さの推定誤差は蓄積されなくなるため、配管２００の内表面における析出物２０１の厚さの推定誤差を所定範囲に収めることが可能である。また、推定部４３は、配管２００の外表面温度の測定値をセンサ２ごとに取得するので、析出物２０１の形状（厚さ分布）を推定することも可能である。
【0068】
なお、センサ２－１～２－４の配置は、図２に示す配置以外の配置であっても良い。例えば、センサ２－１～２－４は、配管２００の外周方向に不均等な間隔で、配管２００の外表面に設置されてもよい。例えば、センサ２－１～２－４は、図３に例示されているように析出物２０１が厚くなり易い位置（例えば、配管２００の底部）の近傍において、密に設置されてもよく、厚くなり難い位置（例えば、配管２００の頂部）の近傍において、疎に設置されてもよい。また、外周方向に設置された複数のセンサ２（センサ２－１～２－４）の組は、配管２００の管軸方向に、予め定められた間隔で配管２００の外表面に設置されてもよい。これによって、配管２００の管軸方向についても、配管２００の内表面における析出物２０１の厚さの推定誤差を所定範囲に収めることが可能である。
【0069】
また、配管２００の内側に器具３００を通過させることによって析出物２０１を除去する代わりに、析出物２０１の発生を阻害する阻害剤を配管２００に投入することによって析出物２０１を除去するようにしてもよい。この場合、検出部３としては、配管２００内の阻害剤の有無を検出（または推測）可能なセンサ、または、阻害剤を投入する機器等の投入タイミングを検出（または推測）可能なセンサ等を利用することができる。配管２００の内側から析出物２０１を除去する方法には、上述したものに限られず、どのような方法が適用されてもよい。
【0070】
〈第１変形例〉
上述した実施形態では、センサ２－１～２－４が配管２００の外表面温度を測定するためにのみ用いられていた。本変形例は、図２（ａ）に示すセンサ２－１～２－４を、配管２００の外表面温度に加えて、配管２００の外表面の熱流束（例えば、配管２００の外表面と略直交する方向（法線方向）の熱流束）を測定するために用いるものである。センサ２－１～２－４は、その設置位置において、配管２００の外表面温度の測定値と、熱流束の測定値とを、通信回線１００を経由して推定部４３に送信する。
上述した実施形態では、センサ２－１～２－４が配管２００の外表面温度を測定するためにのみ用いられていた。本変形例は、図２（ａ）に示すセンサ２－１～２－４を、配管２００の外表面温度に加えて、配管２００の外表面の熱流束（例えば、配管２００の外表面と略直交する方向（法線方向）の熱流束）を測定するために用いるものである。センサ２－１～２－４は、その設置位置において、配管２００の外表面温度の測定値と、熱流束の測定値とを、通信回線１００を経由して推定部４３に送信する。
【0071】
本変形例において、推定部４３は、少なくとも１個のセンサ２の設置位置における外表面温度の測定値と、熱流束の測定値とに基づいて、そのセンサ２の設置位置における析出物２０１の厚さを式（５）のように推定する。
【0072】
ｔｈ＝ｆ（ｐ１，ｐ２，ｐ３，ｐ４，ｐ５，ｐ６，ｐ７） …（５）
【0073】
ここで、パラメータ「ｐ７」は、センサ２－１～２－４のうちの少なくとも１個のセンサ２による、配管２００の外表面における熱流束の測定値を表す。なお、式（５）において、パラメータ「ｐ５」及びパラメータ「ｐ６」は必須ではなく、無くてもよい。
【0074】
配管２００の外側における熱の環境と配管２００の内側における熱の環境とが共に変化する場合でも、推定部４３は、配管２００の外表面における熱流束の測定値を用いて、配管２００の外側の環境変化の影響と配管２００の内側の環境変化の影響とを分離する。これによって、熱流束の測定値を用いて配管２００の外側及び内側における環境の影響を補正して、配管２００の内表面における析出物２０１の厚さの推定誤差を所定範囲に収めることが可能である。
【0075】
推定部４３は、配管２００の外側及び内側における熱に関する推定結果を、出力部４５を用いてオペレータに通知してもよい。これによって、推定部４３は、例えば、パイプラインの管理において有用な情報を、オペレータに通知することが可能になる。パイプラインの管理において有用な情報とは、例えば、配管２００の外側における大気又は海水の温度変化と、配管２００の内側におけるオイル及びガスの各温度と、配管２００の内側における混相流の流動様式である。
【0076】
〈第２変形例〉
図７は、第１実施形態の第２変形例における、配管構成例を示す断面図である。なお、図７に示す断面図は、配管２００の半径方向の断面図である。本変形例は、図２（ｂ）に示す通り、配管２００の外表面に断熱材２３０が設けられている場合においても、配管２００の外表面温度に加えて、配管２００の外表面の熱流束の測定を可能にするものである。
図７は、第１実施形態の第２変形例における、配管構成例を示す断面図である。なお、図７に示す断面図は、配管２００の半径方向の断面図である。本変形例は、図２（ｂ）に示す通り、配管２００の外表面に断熱材２３０が設けられている場合においても、配管２００の外表面温度に加えて、配管２００の外表面の熱流束の測定を可能にするものである。
【0077】
図７に示す通り、本変形例では、センサ２－１～２－４に加えて、センサ２－５～２－８が設けられている。センサ２－５～２－８は、断熱材２３０内であって、センサ２－１～２－４から配管２００の半径方向に所定距離だけそれぞれ離れた位置に、配管２００の外周方向に沿って等間隔（等角度間隔であって、例えば９０度の間隔）で設置されている。センサ２－１～２－４とセンサ２－５～２－８とは、それぞれ対応して設けられている。
【0078】
センサ２－５～２－８は、設置位置において、配管２００の外表面温度に応じた断熱材２３０の内部温度を測定する。センサ２－５～２－８は、設置位置における断熱材２３０の内部温度の測定値を、通信回線１００を経由して推定部４３に送信する。
【0079】
推定部４３は、センサ２－１～２－４のうちの少なくとも１個のセンサの設置位置における外表面温度の測定値と、センサ２－５～２－８のうちの対応する少なくとも１個のセンサの設置位置における断熱材２３０の内部温度の測定値とに基づいて、センサ２－１～２－４の設置位置における析出物２０１の厚さを式（６）のように推定する。
【0080】
ｔｈ＝ｆ（ｐ１，ｐ２，ｐ３，ｐ４，ｐ５，ｐ６，ｐ８） …（６）
【0081】
ここで、パラメータ「ｐ８」は、センサ２－５～２－８のうちの少なくとも１個のセンサによる、断熱材２３０の内部温度の測定値を表す。なお、式（６）において、パラメータ「ｐ５」及びパラメータ「ｐ６」は必須ではなく、無くてもよい。
【0082】
本変形例においても、第１変形例と同様に、配管２００の外側及び内側における環境の影響を補正して、配管２００の内表面における析出物２０１の厚さの推定誤差を所定範囲に収めることが可能である。なお、本変形例においても、推定部４３は、配管２００の外側及び内側における熱に関する推定結果を、出力部４５を用いてオペレータに通知してもよい。
【0083】
〔第２実施形態〕
【0084】
〈推定システム〉
図８は、第２実施形態における、推定システムの構成例を示す図である。本実施形態の推定システム１ｂは、図１に示す推定システム１ａの推定装置４ａを推定装置４ｂに代えた構成である。図８に示す推定システム１ｂは、図１に示す推定システム１ａと同様に、配管の内表面における析出物の厚さを推定するシステムである。但し、図８に示す推定システム１ｂは、配管２００の所定位置を器具３００が通過すること又は通過したことを予測（推測）することによって、その予測結果（予測時刻）を、析出物の除去処理が実行される実行時刻（配管２００における所定の位置を器具３００が通過する通過時刻）として推測する点が、図１に示す推定システム１ａとは異なる。配管２００の所定位置とは、例えば、センサ２の設置位置である。
図８は、第２実施形態における、推定システムの構成例を示す図である。本実施形態の推定システム１ｂは、図１に示す推定システム１ａの推定装置４ａを推定装置４ｂに代えた構成である。図８に示す推定システム１ｂは、図１に示す推定システム１ａと同様に、配管の内表面における析出物の厚さを推定するシステムである。但し、図８に示す推定システム１ｂは、配管２００の所定位置を器具３００が通過すること又は通過したことを予測（推測）することによって、その予測結果（予測時刻）を、析出物の除去処理が実行される実行時刻（配管２００における所定の位置を器具３００が通過する通過時刻）として推測する点が、図１に示す推定システム１ａとは異なる。配管２００の所定位置とは、例えば、センサ２の設置位置である。
【0085】
〈推定装置〉
推定装置４ｂは、図１に示す推定装置４ａに予測部（推測部）４６を追加した構成である。予測部４６は、器具３００（ピグ）がピグ発射装置（pig launcher）（不図示）から発射された発射時刻の情報と、器具３００がピグ回収装置（pig receiver）（不図示）において回収された回収時刻の情報とのうちの少なくとも一つを取得する。なお、推定装置４ｂは、上記の発射時刻または回収時刻の情報を、ピグから直接、あるいは、ピグ発射装置又はピグ回収装置のうちの少なくとも一方を制御する制御端末（不図示）から取得する。
推定装置４ｂは、図１に示す推定装置４ａに予測部（推測部）４６を追加した構成である。予測部４６は、器具３００（ピグ）がピグ発射装置（pig launcher）（不図示）から発射された発射時刻の情報と、器具３００がピグ回収装置（pig receiver）（不図示）において回収された回収時刻の情報とのうちの少なくとも一つを取得する。なお、推定装置４ｂは、上記の発射時刻または回収時刻の情報を、ピグから直接、あるいは、ピグ発射装置又はピグ回収装置のうちの少なくとも一方を制御する制御端末（不図示）から取得する。
【0086】
予測部４６は、ピグが発射された位置とピグが回収された位置とのうちの少なくとも一つの情報を、記憶部４１から取得する。予測部４６は、ピグが発射された時刻及びピグが発射された位置とピグが回収された時刻及びピグが回収された位置とのうちの少なくとも一つと、所定位置と、器具３００の移動速度とに基づいて、所定位置を器具３００が通過すること又は通過したこと（通過時刻）を予測（推測）する。予測部４６は、器具３００が通過すること又は通過したことを予測（推測）することによって、その予測結果（予測時刻）を、配管２００のうち表面における析出物２０１の除去処理の実行時刻（タイミング）として推測する。予測部４６は、所定位置を器具３００が通過する通過時刻の予測結果、すなわち、析出物２０１の除去処理の実行時刻の推測結果を、記憶部４１に記憶させる。判定部４４は、器具３００の通過を予測部４６が予測（推測）したか否かを判定する。すなわち、判定部４４は、析出物２０１の除去処理が実行されたか否かを判定する。
【0087】
補正部４２は、析出物２０１の厚さの推定に用いられる推定値を、所定位置を器具３００が通過するであろう又は通過したであろう予測時刻から所定時間内（例えば、１秒以内）に得られた配管２００の外表面温度の測定値を用いて補正する。もちろん、検出部３の設置位置を器具３００が通過する通過時刻を予測する際には、上述したようにピグの発射または回収タイミングと通過推測位置とピグ移動速度とに基づいて通過時刻を予測するものに限定されず、予め記憶された移動時間に基づいて通過時刻を予測してもよい。具体的には、ピグがピグ発射装置を発射してから通過推測位置を通過するまでの移動時間、または、ピグが通過推測位置を通過してからピグ回収装置に回収されるまでの移動時間を、予め経験値または計算値で記憶部４１に記憶しておき、この記憶部４１に記憶された経験値または計算値に基づいて通過時刻を予測（推測）することができる。
【0088】
以上説明した通り、本実施形態では、補正部４２が、析出物２０１の厚さの推定に用いられる少なくとも一つの推定値を、除去処理の実行時刻としての予測部４６の予測結果（予測時刻）から所定時間内に得られた配管２００の外表面温度の測定値を用いて補正するようにしている。ここで、補正部４２は、析出物２０１の厚さの推定結果が予め規定された閾値以下（例えば、０）になるように推定値を補正してもよい。推定値が補正されることによって、析出物２０１の厚さの推定誤差は蓄積されなくなるため、配管２００の内表面における析出物２０１の厚さの推定誤差を所定範囲に収めることが可能である。よって、第1実施形態の検出部３を備えていない場合でも、本実施形態では配管２００の内表面における析出物２０１の厚さの推定誤差を所定範囲に収めることが可能である。また、たとえ第1実施形態の検出部３を備えているにも関わらず、第１実施形態の検出部３が器具３００の通過を検出することによって析出物２０１の除去処理の実行時刻を推測できない場合であっても、本実施形態では配管２００の内表面における析出物２０１の厚さの推定誤差を所定範囲に収めることが可能である。なお、第１実施形態の検出部３と第２実施形態の予測部４６とを推定システム１ｂが備えることで、予測部４６によるピグの通過予測と、検出部３によるピグの通過検出とのうちの少なくとも一方又は両方が実行されてもよい。
【0089】
以上、この発明の実施形態について図面を参照して詳述してきたが、具体的な構成はこの実施形態に限られるものではなく、この発明の要旨を逸脱しない範囲の設計等も含まれる。
【0090】
例えば、予測部４６は、析出物２０１の厚さ推定の結果が急激に変化（減少）した際に、器具３００が通過したと推定し（みなし）、推定値の補正フローが実行されてもよい。又は、予測部４６は、センサ２の測定値（外表面温度）が急激に変化した際に、器具３００が通過したと推定し（みなし）、推定値の補正フローが実行されてもよい。
【符号の説明】
【0091】
１ａ，１ｂ…推定システム、２…センサ、３…検出部、４ａ，４ｂ…推定装置、４０…取得部、４１…記憶部、４２…補正部、４３…推定部、４４…判定部、４５…出力部、４６…予測部、１００…通信回線、２００…配管、２０１…析出物、２１０…オイル、２２０…ガス、２３０…断熱材、３００…器具
【図1】
【図2】
【図3】
【図4】
【図5】
【図6】
【図7】
【図8】