▶ 東京電力株式会社の特許一覧 ▶ 国立大学法人神戸大学の特許一覧
特許7246634流動様式判別装置、流動様式判別システムおよび流動様式判別方法
(19)【発行国】日本国特許庁(JP)
(12)【公報種別】特許公報(B2)
(11)【特許番号】
(24)【登録日】2023-03-17
(45)【発行日】2023-03-28
(54)【発明の名称】流動様式判別装置、流動様式判別システムおよび流動様式判別方法
(51)【国際特許分類】
G01N 29/48 20060101AFI20230320BHJP
G01N 29/36 20060101ALI20230320BHJP
G01N 29/024 20060101ALI20230320BHJP
G01F 1/66 20220101ALI20230320BHJP
G01F 1/00 20220101ALI20230320BHJP
【ＦＩ】
G01N29/48
G01N29/36
G01N29/024
G01F1/66 101
G01F1/66 Z
G01F1/00 Z
【請求項の数】
11
(21)【出願番号】P 2019058120
(22)【出願日】2019-03-26
(65)【公開番号】P2020159809
(43)【公開日】2020-10-01
【審査請求日】2022-03-18
【新規性喪失の例外の表示】特許法第３０条第２項適用 １．刊行物 ・発行日 平成３０年１１月２９日 日本機械学会第９６期流体工学部門講演会 講演論文集 ２．学会発表（ポスター） ・開催日 平成３０年１１月２９日～３０日（発表日：平成３０年１１月２９日） ・集会名 日本機械学会第９６期流体工学部門講演会 ・開催場所 蓬▲らい▼殿（北海道室蘭市宮の森町１－１）
(73)【特許権者】
【識別番号】000003687
【氏名又は名称】東京電力ホールディングス株式会社
(73)【特許権者】
【識別番号】504150450
【氏名又は名称】国立大学法人神戸大学
(74)【代理人】
【識別番号】100149548
【弁理士】
【氏名又は名称】松沼 泰史
(74)【代理人】
【識別番号】100175824
【弁理士】
【氏名又は名称】小林 淳一
(74)【代理人】
【識別番号】100140774
【弁理士】
【氏名又は名称】大浪 一徳
(74)【代理人】
【識別番号】100179833
【弁理士】
【氏名又は名称】松本 将尚
(74)【代理人】
【識別番号】100114937
【弁理士】
【氏名又は名称】松本 裕幸
(72)【発明者】
【氏名】梅沢 修一
(72)【発明者】
【氏名】杉田 勝彦
(72)【発明者】
【氏名】横坂 雅樹
(72)【発明者】
【氏名】村川 英樹
(72)【発明者】
【氏名】杉本 勝美
(72)【発明者】
【氏名】市村 修平
【審査官】清水 靖記
(56)【参考文献】
【文献】特開２００８－９６１７８（ＪＰ，Ａ）
【文献】特開２０１２－６３１８７（ＪＰ，Ａ）
【文献】特開平５－３０５０８７（ＪＰ，Ａ）
【文献】特開平１１－３０９１４７（ＪＰ，Ａ）
【文献】米国特許第６３８６０１８（ＵＳ，Ｂ１）
【文献】村川英樹 市村修平 島田昌輝 杉本勝美 梅沢修一 杉田勝彦 ，クランプオン式超音波流量計による湿り蒸気計測のための透過波形評価手法，日本機械学会動力・エネルギー技術シンポジウム講演論文集（ＣＤ－ＲＯＭ），2019年，A221, 1-5
(58)【調査した分野】(Int.Cl.，ＤＢ名)
Ｇ０１Ｎ ２９／００ － Ｇ０１Ｎ ２９／５２
Ｇ０１Ｂ １７／００ － Ｇ０１Ｂ １７／０８
Ｇ０１Ｈ １／００ － Ｇ０１Ｈ １７／００
Ｇ０１Ｎ １１／００ － Ｇ０１Ｎ １１／１６
Ｇ０１Ｆ １／００ － Ｇ０１Ｆ ９／０２
ＪＳＴＰｌｕｓ／ＪＭＥＤＰｌｕｓ／ＪＳＴ７５８０（ＪＤｒｅａｍＩＩＩ）
Ｓｃｏｐｕｓ
(57)【特許請求の範囲】
【請求項1】
配管内を流れる流体の流動様式を判別する流動様式判別装置であって、超音波を送信する第１超音波トランスデューサと、
前記第１超音波トランスデューサから送信された前記超音波を受信する第２超音波トランスデューサと、
前記超音波を受信した前記第２超音波トランスデューサによって出力される受信信号に基づいて、前記流体の流動様式を判別する演算部とを備え、
前記第１超音波トランスデューサは、前記超音波を複数回送信し、
前記第２超音波トランスデューサは、前記第１超音波トランスデューサから複数回送信された前記超音波を受信して、複数回分の前記受信信号を出力し、
前記演算部は、前記第１超音波トランスデューサによる前記超音波の送信時刻以降の期間である送信後期間中の各時刻における複数回分の前記受信信号の振幅の標準偏差、分散、および、標準偏差から得られる値のいずれかである標準偏差同等値に基づいて、前記流体の流動様式を判別し、
前記送信後期間は、前記第１超音波トランスデューサから送信され、前記配管の壁部内を伝わったガイド波としての前記超音波が前記第２超音波トランスデューサによって受信される期間である、
流動様式判別装置。
【請求項2】
前記演算部は、前記送信後期間中の各時刻における複数回分の前記受信信号の振幅の標準偏差同等値に基づいて、前記流体の流動様式が、乾き蒸気、層状流、波状流、波状環状遷移および環状噴霧流のいずれであるかを判別する、請求項１に記載の流動様式判別装置。
【請求項3】
前記流体の流動様式が層状流である場合における前記送信後期間中の各時刻における複数回分の前記受信信号の振幅の標準偏差同等値は、前記流体の流動様式が乾き蒸気である場合における前記送信後期間中の各時刻における複数回分の前記受信信号の振幅の標準偏差同等値よりも大きく、前記流体の流動様式が波状流である場合における前記送信後期間中の各時刻における複数回分の前記受信信号の振幅の標準偏差同等値は、前記流体の流動様式が層状流である場合における前記送信後期間中の各時刻における複数回分の前記受信信号の振幅の標準偏差同等値よりも大きく、
前記流体の流動様式が波状環状遷移である場合における前記送信後期間中の各時刻における複数回分の前記受信信号の振幅の標準偏差同等値は、前記流体の流動様式が波状流である場合における前記送信後期間中の各時刻における複数回分の前記受信信号の振幅の標準偏差同等値よりも大きく、
前記流体の流動様式が環状噴霧流である場合における前記送信後期間中の各時刻における複数回分の前記受信信号の振幅の標準偏差同等値は、前記流体の流動様式が波状環状遷移である場合における前記送信後期間中の各時刻における複数回分の前記受信信号の振幅の標準偏差同等値よりも大きい、
請求項２に記載の流動様式判別装置。
【請求項4】
前記第１超音波トランスデューサは、第１送信部と第１受信部とを有し、前記第２超音波トランスデューサは、第２送信部と第２受信部とを有し、前記第１超音波トランスデューサよりも前記流体の流れの下流側に配置されており、
前記第２受信部は、前記第１送信部から複数回送信された前記超音波を受信して、複数回分の前記受信信号を出力し、
前記第１受信部は、前記第２送信部から複数回送信された前記超音波を受信して、複数回分の前記受信信号を出力し、
前記演算部は、
前記第２受信部による前記第１送信部から送信された前記超音波の受信期間を含む期間中の各時刻における前記第２受信部の複数回分の前記受信信号の振幅の標準偏差同等値を算出し、
前記第１受信部による前記第２送信部から送信された前記超音波の受信期間を含む期間中の各時刻における前記第１受信部の複数回分の前記受信信号の振幅の標準偏差同等値を算出し、
前記第２受信部の複数回分の前記受信信号の振幅の標準偏差同等値と、前記第１受信部の複数回分の前記受信信号の振幅の標準偏差同等値とに基づいて、第１超音波の伝播時間と第２超音波の伝播時間との差分である伝播時間差を算出するか、あるいは、前記第２受信部の前記受信信号と、前記第１受信部の前記受信信号とに基づいて、前記第１超音波の伝播時間と前記第２超音波の伝播時間との差分である伝播時間差を算出することによって、前記流体の流量を算出し、
前記第２受信部による前記第１送信部から送信された前記超音波の受信期間は、前記第１送信部から送信され、前記配管の流路内を透過した目的信号としての前記超音波が前記第２受信部によって受信される期間であり、
前記第１受信部による前記第２送信部から送信された前記超音波の受信期間は、前記第２送信部から送信され、前記配管の前記流路内を透過した前記目的信号としての前記超音波が前記第１受信部によって受信される期間である、
請求項２または請求項３に記載の流動様式判別装置。
【請求項5】
前記演算部は、予め算出された前記流体の流動様式と前記配管内の気相容積流束と前記配管内の液相容積流束との関係と、
前記演算部によって判別された前記流体の流動様式と、
前記演算部によって算出された前記流体の流量とに基づいて、
前記流体の湿り度を算出する、
請求項４に記載の流動様式判別装置。
【請求項6】
配管内を流れる流体の流動様式を判別する流動様式判別システムであって、前記流動様式判別システムは、超音波流量計と、演算部とを備え、
前記超音波流量計は、
超音波を送信する第１超音波トランスデューサと、
前記第１超音波トランスデューサから送信された前記超音波を受信する第２超音波トランスデューサとを備え、
前記第１超音波トランスデューサは、前記超音波を複数回送信し、
前記第２超音波トランスデューサは、前記第１超音波トランスデューサから複数回送信された前記超音波を受信して、複数回分の受信信号を出力し、
前記演算部は、前記第１超音波トランスデューサによる前記超音波の送信時刻以降の期間である送信後期間中の各時刻における複数回分の前記受信信号の振幅の標準偏差、分散、および、標準偏差から得られる値のいずれかである標準偏差同等値に基づいて、前記流体の流動様式を判別し、
前記送信後期間は、前記第１超音波トランスデューサから送信され、前記配管の壁部内を伝わったガイド波としての前記超音波が前記第２超音波トランスデューサによって受信される期間である、
流動様式判別システム。
【請求項7】
前記演算部は、前記送信後期間中の各時刻における複数回分の前記受信信号の振幅の標準偏差同等値に基づいて、前記流体の流動様式が、乾き蒸気、層状流、波状流、波状環状遷移および環状噴霧流のいずれであるかを判別する、請求項６に記載の流動様式判別システム。
【請求項8】
前記流体の流動様式が層状流である場合における前記送信後期間中の各時刻における複数回分の前記受信信号の振幅の標準偏差同等値は、前記流体の流動様式が乾き蒸気である場合における前記送信後期間中の各時刻における複数回分の前記受信信号の振幅の標準偏差同等値よりも大きく、前記流体の流動様式が波状流である場合における前記送信後期間中の各時刻における複数回分の前記受信信号の振幅の標準偏差同等値は、前記流体の流動様式が層状流である場合における前記送信後期間中の各時刻における複数回分の前記受信信号の振幅の標準偏差同等値よりも大きく、
前記流体の流動様式が波状環状遷移である場合における前記送信後期間中の各時刻における複数回分の前記受信信号の振幅の標準偏差同等値は、前記流体の流動様式が波状流である場合における前記送信後期間中の各時刻における複数回分の前記受信信号の振幅の標準偏差同等値よりも大きく、
前記流体の流動様式が環状噴霧流である場合における前記送信後期間中の各時刻における複数回分の前記受信信号の振幅の標準偏差同等値は、前記流体の流動様式が波状環状遷移である場合における前記送信後期間中の各時刻における複数回分の前記受信信号の振幅の標準偏差同等値よりも大きい、
請求項７に記載の流動様式判別システム。
【請求項9】
前記第１超音波トランスデューサは、第１送信部と第１受信部とを有し、前記第２超音波トランスデューサは、第２送信部と第２受信部とを有し、前記第１超音波トランスデューサよりも前記流体の流れの下流側に配置されており、
前記第２受信部は、前記第１送信部から複数回送信された前記超音波を受信して、複数回分の前記受信信号を出力し、
前記第１受信部は、前記第２送信部から複数回送信された前記超音波を受信して、複数回分の前記受信信号を出力し、
前記演算部は、
前記第２受信部による前記第１送信部から送信された前記超音波の受信期間を含む期間中の各時刻における前記第２受信部の複数回分の前記受信信号の振幅の標準偏差同等値を算出し、
前記第１受信部による前記第２送信部から送信された前記超音波の受信期間を含む期間中の各時刻における前記第１受信部の複数回分の前記受信信号の振幅の標準偏差同等値を算出し、
前記第２受信部の複数回分の前記受信信号の振幅の標準偏差同等値と、前記第１受信部の複数回分の前記受信信号の振幅の標準偏差同等値とに基づいて、第１超音波の伝播時間と第２超音波の伝播時間との差分である伝播時間差を算出するか、あるいは、前記第２受信部の前記受信信号と、前記第１受信部の前記受信信号とに基づいて、前記第１超音波の伝播時間と前記第２超音波の伝播時間との差分である伝播時間差を算出することによって、前記流体の流量を算出し、
前記第２受信部による前記第１送信部から送信された前記超音波の受信期間は、前記第１送信部から送信され、前記配管の流路内を透過した目的信号としての前記超音波が前記第２受信部によって受信される期間であり、
前記第１受信部による前記第２送信部から送信された前記超音波の受信期間は、前記第２送信部から送信され、前記配管の前記流路内を透過した前記目的信号としての前記超音波が前記第１受信部によって受信される期間である、
請求項７または請求項８に記載の流動様式判別システム。
【請求項10】
前記演算部は、予め算出された前記流体の流動様式と前記配管内の気相容積流束と前記配管内の液相容積流束との関係と、
前記演算部によって判別された前記流体の流動様式と、
前記演算部によって算出された前記流体の流量とに基づいて、
前記流体の湿り度を算出する、
請求項９に記載の流動様式判別システム。
【請求項11】
超音波を送信する第１超音波トランスデューサと、前記第１超音波トランスデューサから送信された前記超音波を受信する第２超音波トランスデューサと
を用いて配管内を流れる流体の流動様式を判別する流動様式判別方法であって、
前記超音波を受信した前記第２超音波トランスデューサによって出力される受信信号に基づいて、前記流体の流動様式を判別する演算ステップと、
前記第１超音波トランスデューサが前記超音波を複数回送信する超音波送信ステップと、
前記第２超音波トランスデューサが前記第１超音波トランスデューサから複数回送信された前記超音波を受信して、複数回分の前記受信信号を出力する受信信号出力ステップとを備え、
前記演算ステップでは、前記第１超音波トランスデューサによる前記超音波の送信時刻以降の期間である送信後期間中の各時刻における複数回分の前記受信信号の振幅の標準偏差、分散、および、標準偏差から得られる値のいずれかである標準偏差同等値に基づいて、前記流体の流動様式を判別し、
前記送信後期間は、前記第１超音波トランスデューサから送信され、前記配管の壁部内を伝わったガイド波としての前記超音波が前記第２超音波トランスデューサによって受信される期間である、
流動様式判別方法。
【発明の詳細な説明】
【技術分野】
【0001】
本発明は、流動様式判別装置、流動様式判別システムおよび流動様式判別方法に関する。
【背景技術】
【0002】
従来から、上流側トランスデューサと下流側トランスデューサとを用いて管路の内部を流れる流体の流量を測定する超音波流量計が知られている（例えば特許文献１参照）。
特許文献１に記載された超音波流量計では、上流側トランスデューサによって送信された超音波が、下流側トランスデューサにおいて受信され、制御部が、上流から下流への超音波の伝播時間を算出する。また、下流側トランスデューサによって送信された超音波が、上流側トランスデューサにおいて受信され、制御部が、下流から上流への超音波の伝播時間を算出する。更に、制御部が、下流から上流への超音波の伝播時間から、上流から下流への超音波の伝播時間を減算して伝播時間差を算出し、伝搬時間差を用いて流体の流速および流量を演算する。
特許文献１に記載された超音波流量計では、上流側トランスデューサによって送信された超音波が、下流側トランスデューサにおいて受信され、制御部が、上流から下流への超音波の伝播時間を算出する。また、下流側トランスデューサによって送信された超音波が、上流側トランスデューサにおいて受信され、制御部が、下流から上流への超音波の伝播時間を算出する。更に、制御部が、下流から上流への超音波の伝播時間から、上流から下流への超音波の伝播時間を減算して伝播時間差を算出し、伝搬時間差を用いて流体の流速および流量を演算する。
【0003】
ところで、特許文献１に記載された技術では、管路の内部を流れる流体の流量が測定されるものの、管路の内部を流れる流体の流動様式は判別されない。
従来の一般的な技術においても、抜管する必要なく、例えば蒸気配管などのような配管内を流れる流体の流動様式を判別することは困難であり、提唱されている方式も存在しなかった。
従来の一般的な技術においても、抜管する必要なく、例えば蒸気配管などのような配管内を流れる流体の流動様式を判別することは困難であり、提唱されている方式も存在しなかった。
【先行技術文献】
【特許文献】
【0004】
【文献】特開２０１５－２２２１８８号公報
【発明の概要】
【発明が解決しようとする課題】
【0005】
本発明者等は、鋭意研究において、一方の超音波トランスデューサが超音波を複数回送信し、他方の超音波トランスデューサが一方の超音波トランスデューサから複数回送信された超音波を受信して複数回分の受信信号を出力し、更に、一方の超音波トランスデューサから送信されて配管の壁部内を伝わったガイド波としての超音波が他方の超音波トランスデューサによって受信される期間中の各時刻における複数回分の受信信号の振幅の標準偏差、分散、および、標準偏差から得られる値のいずれかの時間波形を比較することによって、配管内を流れる流体の流動様式を判別できることを見い出したのである。
つまり、本発明は、配管内を流れる流体の流動様式を判別することができる流動様式判別装置、流動様式判別システムおよび流動様式判別方法を提供することを目的とする。
つまり、本発明は、配管内を流れる流体の流動様式を判別することができる流動様式判別装置、流動様式判別システムおよび流動様式判別方法を提供することを目的とする。
【課題を解決するための手段】
【0006】
本発明の一態様は、配管内を流れる流体の流動様式を判別する流動様式判別装置であって、超音波を送信する第１超音波トランスデューサと、前記第１超音波トランスデューサから送信された前記超音波を受信する第２超音波トランスデューサと、前記超音波を受信した前記第２超音波トランスデューサによって出力される受信信号に基づいて、前記流体の流動様式を判別する演算部とを備え、前記第１超音波トランスデューサは、前記超音波を複数回送信し、前記第２超音波トランスデューサは、前記第１超音波トランスデューサから複数回送信された前記超音波を受信して、複数回分の前記受信信号を出力し、前記演算部は、前記第１超音波トランスデューサによる前記超音波の送信時刻以降の期間である送信後期間中の各時刻における複数回分の前記受信信号の振幅の標準偏差、分散、および、標準偏差から得られる値のいずれかである標準偏差同等値に基づいて、前記流体の流動様式を判別し、前記送信後期間は、前記第１超音波トランスデューサから送信され、前記配管の壁部内を伝わったガイド波としての前記超音波が前記第２超音波トランスデューサによって受信される期間である、流動様式判別装置である。
【0007】
本発明の一態様の流動様式判別装置では、前記演算部は、前記送信後期間中の各時刻における複数回分の前記受信信号の振幅の標準偏差同等値に基づいて、前記流体の流動様式が、乾き蒸気、層状流、波状流、波状環状遷移および環状噴霧流のいずれであるかを判別してもよい。
【0008】
本発明の一態様の流動様式判別装置では、前記流体の流動様式が層状流である場合における前記送信後期間中の各時刻における複数回分の前記受信信号の振幅の標準偏差同等値は、前記流体の流動様式が乾き蒸気である場合における前記送信後期間中の各時刻における複数回分の前記受信信号の振幅の標準偏差同等値よりも大きく、前記流体の流動様式が波状流である場合における前記送信後期間中の各時刻における複数回分の前記受信信号の振幅の標準偏差同等値は、前記流体の流動様式が層状流である場合における前記送信後期間中の各時刻における複数回分の前記受信信号の振幅の標準偏差同等値よりも大きく、前記流体の流動様式が波状環状遷移である場合における前記送信後期間中の各時刻における複数回分の前記受信信号の振幅の標準偏差同等値は、前記流体の流動様式が波状流である場合における前記送信後期間中の各時刻における複数回分の前記受信信号の振幅の標準偏差同等値よりも大きく、前記流体の流動様式が環状噴霧流である場合における前記送信後期間中の各時刻における複数回分の前記受信信号の振幅の標準偏差同等値は、前記流体の流動様式が波状環状遷移である場合における前記送信後期間中の各時刻における複数回分の前記受信信号の振幅の標準偏差同等値よりも大きくてもよい。
【0009】
本発明の一態様の流動様式判別装置では、前記第１超音波トランスデューサは、第１送信部と第１受信部とを有し、前記第２超音波トランスデューサは、第２送信部と第２受信部とを有し、前記第１超音波トランスデューサよりも前記流体の流れの下流側に配置されており、前記第２受信部は、前記第１送信部から複数回送信された前記超音波を受信して、複数回分の前記受信信号を出力し、前記第１受信部は、前記第２送信部から複数回送信された前記超音波を受信して、複数回分の前記受信信号を出力し、前記演算部は、前記第２受信部による前記第１送信部から送信された前記超音波の受信期間を含む期間中の各時刻における前記第２受信部の複数回分の前記受信信号の振幅の標準偏差同等値を算出し、前記第１受信部による前記第２送信部から送信された前記超音波の受信期間を含む期間中の各時刻における前記第１受信部の複数回分の前記受信信号の振幅の標準偏差同等値を算出し、前記第２受信部の複数回分の前記受信信号の振幅の標準偏差同等値と、前記第１受信部の複数回分の前記受信信号の振幅の標準偏差同等値とに基づいて、第１超音波の伝播時間と第２超音波の伝播時間との差分である伝播時間差を算出するか、あるいは、前記第２受信部の前記受信信号と、前記第１受信部の前記受信信号とに基づいて、前記第１超音波の伝播時間と前記第２超音波の伝播時間との差分である伝播時間差を算出することによって、前記流体の流量を算出し、前記第２受信部による前記第１送信部から送信された前記超音波の受信期間は、前記第１送信部から送信され、前記配管の流路内を透過した目的信号としての前記超音波が前記第２受信部によって受信される期間であり、前記第１受信部による前記第２送信部から送信された前記超音波の受信期間は、前記第２送信部から送信され、前記配管の前記流路内を透過した前記目的信号としての前記超音波が前記第１受信部によって受信される期間であってもよい。
【0010】
本発明の一態様の流動様式判別装置では、前記演算部は、予め算出された前記流体の流動様式と前記配管内の気相容積流束と前記配管内の液相容積流束との関係と、前記演算部によって判別された前記流体の流動様式と、前記演算部によって算出された前記流体の流量とに基づいて、前記流体の湿り度を算出してもよい。
【0011】
本発明の一態様は、配管内を流れる流体の流動様式を判別する流動様式判別システムであって、前記流動様式判別システムは、超音波流量計と、演算部とを備え、前記超音波流量計は、超音波を送信する第１超音波トランスデューサと、前記第１超音波トランスデューサから送信された前記超音波を受信する第２超音波トランスデューサとを備え、前記第１超音波トランスデューサは、前記超音波を複数回送信し、前記第２超音波トランスデューサは、前記第１超音波トランスデューサから複数回送信された前記超音波を受信して、複数回分の受信信号を出力し、前記演算部は、前記第１超音波トランスデューサによる前記超音波の送信時刻以降の期間である送信後期間中の各時刻における複数回分の前記受信信号の振幅の標準偏差、分散、および、標準偏差から得られる値のいずれかである標準偏差同等値に基づいて、前記流体の流動様式を判別し、前記送信後期間は、前記第１超音波トランスデューサから送信され、前記配管の壁部内を伝わったガイド波としての前記超音波が前記第２超音波トランスデューサによって受信される期間である、流動様式判別システムである。
【0012】
本発明の一態様の流動様式判別システムでは、前記演算部は、前記送信後期間中の各時刻における複数回分の前記受信信号の振幅の標準偏差同等値に基づいて、前記流体の流動様式が、乾き蒸気、層状流、波状流、波状環状遷移および環状噴霧流のいずれであるかを判別してもよい。
【0013】
本発明の一態様の流動様式判別システムでは、前記流体の流動様式が層状流である場合における前記送信後期間中の各時刻における複数回分の前記受信信号の振幅の標準偏差同等値は、前記流体の流動様式が乾き蒸気である場合における前記送信後期間中の各時刻における複数回分の前記受信信号の振幅の標準偏差同等値よりも大きく、前記流体の流動様式が波状流である場合における前記送信後期間中の各時刻における複数回分の前記受信信号の振幅の標準偏差同等値は、前記流体の流動様式が層状流である場合における前記送信後期間中の各時刻における複数回分の前記受信信号の振幅の標準偏差同等値よりも大きく、前記流体の流動様式が波状環状遷移である場合における前記送信後期間中の各時刻における複数回分の前記受信信号の振幅の標準偏差同等値は、前記流体の流動様式が波状流である場合における前記送信後期間中の各時刻における複数回分の前記受信信号の振幅の標準偏差同等値よりも大きく、前記流体の流動様式が環状噴霧流である場合における前記送信後期間中の各時刻における複数回分の前記受信信号の振幅の標準偏差同等値は、前記流体の流動様式が波状環状遷移である場合における前記送信後期間中の各時刻における複数回分の前記受信信号の振幅の標準偏差同等値よりも大きくてもよい。
【0014】
本発明の一態様の流動様式判別システムでは、前記第１超音波トランスデューサは、第１送信部と第１受信部とを有し、前記第２超音波トランスデューサは、第２送信部と第２受信部とを有し、前記第１超音波トランスデューサよりも前記流体の流れの下流側に配置されており、前記第２受信部は、前記第１送信部から複数回送信された前記超音波を受信して、複数回分の前記受信信号を出力し、前記第１受信部は、前記第２送信部から複数回送信された前記超音波を受信して、複数回分の前記受信信号を出力し、前記演算部は、前記第２受信部による前記第１送信部から送信された前記超音波の受信期間を含む期間中の各時刻における前記第２受信部の複数回分の前記受信信号の振幅の標準偏差同等値を算出し、前記第１受信部による前記第２送信部から送信された前記超音波の受信期間を含む期間中の各時刻における前記第１受信部の複数回分の前記受信信号の振幅の標準偏差同等値を算出し、前記第２受信部の複数回分の前記受信信号の振幅の標準偏差同等値と、前記第１受信部の複数回分の前記受信信号の振幅の標準偏差同等値とに基づいて、第１超音波の伝播時間と第２超音波の伝播時間との差分である伝播時間差を算出するか、あるいは、前記第２受信部の前記受信信号と、前記第１受信部の前記受信信号とに基づいて、前記第１超音波の伝播時間と前記第２超音波の伝播時間との差分である伝播時間差を算出することによって、前記流体の流量を算出し、前記第２受信部による前記第１送信部から送信された前記超音波の受信期間は、前記第１送信部から送信され、前記配管の流路内を透過した目的信号としての前記超音波が前記第２受信部によって受信される期間であり、前記第１受信部による前記第２送信部から送信された前記超音波の受信期間は、前記第２送信部から送信され、前記配管の前記流路内を透過した前記目的信号としての前記超音波が前記第１受信部によって受信される期間であってもよい。
【0015】
本発明の一態様の流動様式判別システムでは、前記演算部は、予め算出された前記流体の流動様式と前記配管内の気相容積流束と前記配管内の液相容積流束との関係と、前記演算部によって判別された前記流体の流動様式と、前記演算部によって算出された前記流体の流量とに基づいて、前記流体の湿り度を算出してもよい。
【0016】
本発明の一態様は、超音波を送信する第１超音波トランスデューサと、前記第１超音波トランスデューサから送信された前記超音波を受信する第２超音波トランスデューサとを用いて配管内を流れる流体の流動様式を判別する流動様式判別方法であって、前記超音波を受信した前記第２超音波トランスデューサによって出力される受信信号に基づいて、前記流体の流動様式を判別する演算ステップと、前記第１超音波トランスデューサが前記超音波を複数回送信する超音波送信ステップと、前記第２超音波トランスデューサが前記第１超音波トランスデューサから複数回送信された前記超音波を受信して、複数回分の前記受信信号を出力する受信信号出力ステップとを備え、前記演算ステップでは、前記第１超音波トランスデューサによる前記超音波の送信時刻以降の期間である送信後期間中の各時刻における複数回分の前記受信信号の振幅の標準偏差、分散、および、標準偏差から得られる値のいずれかである標準偏差同等値に基づいて、前記流体の流動様式を判別し、前記送信後期間は、前記第１超音波トランスデューサから送信され、前記配管の壁部内を伝わったガイド波としての前記超音波が前記第２超音波トランスデューサによって受信される期間である、流動様式判別方法である。
【発明の効果】
【0017】
本発明によれば、配管内を流れる流体の流動様式を判別することができる流動様式判別装置、流動様式判別システムおよび流動様式判別方法を提供することができる。
【図面の簡単な説明】
【0018】
【図1】第１実施形態の流動様式判別装置の一例などを示す図である。
【図2】配管の流路内の流体の流れが水を含まない空気の単相流である場合に第２超音波トランスデューサの受信部によって出力される受信信号の時間波形、および、第１超音波トランスデューサの受信部によって出力される受信信号の時間波形の一例を示す図である。
【図3】配管の流路内の流体の流れが水を含まない空気の単相流である場合における第２超音波トランスデューサの受信部の受信信号の振幅の標準偏差同等値の時間波形と、第１超音波トランスデューサの受信部の受信信号の振幅の標準偏差同等値の時間波形とを図２に追記した図である。
【図4】配管の流路内を流れる流体の流動様式が層状流である場合に第２超音波トランスデューサの受信部によって出力される受信信号の時間波形、第１超音波トランスデューサの受信部によって出力される受信信号の時間波形、第２超音波トランスデューサの受信部の受信信号の振幅の標準偏差同等値の時間波形、および、第１超音波トランスデューサの受信部の受信信号の振幅の標準偏差同等値の時間波形の一例を示す図である。
【図5】演算部の流動様式線図取得部によって取得される流動様式線図の一例を示す図である。
【図6】配管の流路内を流れる流体の流動様式が波状流である場合に第２超音波トランスデューサの受信部によって出力される受信信号の時間波形、第１超音波トランスデューサの受信部によって出力される受信信号の時間波形、第２超音波トランスデューサの受信部の受信信号の振幅の標準偏差同等値の時間波形、および、第１超音波トランスデューサの受信部の受信信号の振幅の標準偏差同等値の時間波形の一例を示す図である。
【図7】配管の流路内を流れる流体の流動様式が波状環状遷移である場合に第２超音波トランスデューサの受信部によって出力される受信信号の時間波形、第１超音波トランスデューサの受信部によって出力される受信信号の時間波形、第２超音波トランスデューサの受信部の受信信号の振幅の標準偏差同等値の時間波形、および、第１超音波トランスデューサの受信部の受信信号の振幅の標準偏差同等値の時間波形の一例を示す図である。
【図8】配管の流路内を流れる流体の流動様式が環状噴霧流である場合に第２超音波トランスデューサの受信部によって出力される受信信号の時間波形、第１超音波トランスデューサの受信部によって出力される受信信号の時間波形、第２超音波トランスデューサの受信部の受信信号の振幅の標準偏差同等値の時間波形、および、第１超音波トランスデューサの受信部の受信信号の振幅の標準偏差同等値の時間波形の一例を示す図である。
【図9】第１実施形態の流動様式判別装置において実行される処理の一例を説明するためのフローチャートである。
【図10】第１実施形態の流動様式判別装置において実行される他の処理の一例を説明するためのフローチャートである。
【図11】第２実施形態の流動様式判別システムの一例を示す図である。
【図12】第２実施形態の流動様式判別システムにおいて実行される処理の一例を説明するためのフローチャートである。
【図13】第２実施形態の流動様式判別システムにおいて実行される他の処理の一例を説明するためのフローチャートである。
【発明を実施するための形態】
【0019】
以下、図面を参照し、本発明の流動様式判別装置、流動様式判別システムおよび流動様式判別方法の実施形態について説明する。
【0020】
［第１実施形態］
図１は第１実施形態の流動様式判別装置１の一例などを示す図である。詳細には、図１（Ａ）は第１実施形態の流動様式判別装置１の一例を示しており、図１（Ｂ）は第１実施形態の流動様式判別装置１が適用される配管Ａと第１実施形態の流動様式判別装置１の第１超音波トランスデューサ１１および第２超音波トランスデューサ１２との位置関係の一例を示している。
図１に示す例では、流動様式判別装置１が、配管Ａの流路Ａ２内を流れる流体の流動様式を判別する。流動様式判別装置１は、第１超音波トランスデューサ１１と、第２超音波トランスデューサ１２と、演算部１３と、制御部１４とを備えている。
第１超音波トランスデューサ１１は、送信部１１Ａと、受信部１１Ｂとを備えており、第２超音波トランスデューサ１２は、送信部１２Ａと、受信部１２Ｂとを備えている。第１超音波トランスデューサ１１の送信部１１Ａは、超音波を送信し、第２超音波トランスデューサ１２の受信部１２Ｂは、第１超音波トランスデューサ１１の送信部１１Ａから送信された超音波などを受信する。詳細には、第２超音波トランスデューサ１２の受信部１２Ｂは、受信した第１超音波トランスデューサ１１の送信部１１Ａからの超音波などを受信信号（電圧信号）に変換して出力する。
第２超音波トランスデューサ１２の送信部１２Ａは、超音波を送信し、第１超音波トランスデューサ１１の受信部１１Ｂは、第２超音波トランスデューサ１２の送信部１２Ａから送信された超音波などを受信する。詳細には、第１超音波トランスデューサ１１の受信部１１Ｂは、受信した第２超音波トランスデューサ１２の送信部１２Ａからの超音波などを受信信号（電圧信号）に変換して出力する。
第２超音波トランスデューサ１２は、第１超音波トランスデューサ１１よりも、配管Ａの流路Ａ２内を流れる流体の流れ（図１（Ａ）の右向き）の下流側（図１（Ａ）の右側）に配置されている。
図１は第１実施形態の流動様式判別装置１の一例などを示す図である。詳細には、図１（Ａ）は第１実施形態の流動様式判別装置１の一例を示しており、図１（Ｂ）は第１実施形態の流動様式判別装置１が適用される配管Ａと第１実施形態の流動様式判別装置１の第１超音波トランスデューサ１１および第２超音波トランスデューサ１２との位置関係の一例を示している。
図１に示す例では、流動様式判別装置１が、配管Ａの流路Ａ２内を流れる流体の流動様式を判別する。流動様式判別装置１は、第１超音波トランスデューサ１１と、第２超音波トランスデューサ１２と、演算部１３と、制御部１４とを備えている。
第１超音波トランスデューサ１１は、送信部１１Ａと、受信部１１Ｂとを備えており、第２超音波トランスデューサ１２は、送信部１２Ａと、受信部１２Ｂとを備えている。第１超音波トランスデューサ１１の送信部１１Ａは、超音波を送信し、第２超音波トランスデューサ１２の受信部１２Ｂは、第１超音波トランスデューサ１１の送信部１１Ａから送信された超音波などを受信する。詳細には、第２超音波トランスデューサ１２の受信部１２Ｂは、受信した第１超音波トランスデューサ１１の送信部１１Ａからの超音波などを受信信号（電圧信号）に変換して出力する。
第２超音波トランスデューサ１２の送信部１２Ａは、超音波を送信し、第１超音波トランスデューサ１１の受信部１１Ｂは、第２超音波トランスデューサ１２の送信部１２Ａから送信された超音波などを受信する。詳細には、第１超音波トランスデューサ１１の受信部１１Ｂは、受信した第２超音波トランスデューサ１２の送信部１２Ａからの超音波などを受信信号（電圧信号）に変換して出力する。
第２超音波トランスデューサ１２は、第１超音波トランスデューサ１１よりも、配管Ａの流路Ａ２内を流れる流体の流れ（図１（Ａ）の右向き）の下流側（図１（Ａ）の右側）に配置されている。
【0021】
図１（Ｂ）に示す例では、第１超音波トランスデューサ１１の送信部１１Ａから送信されて配管Ａの流路Ａ２内を透過した超音波が、配管Ａの壁部Ａ１の内側表面で反射することなく第２超音波トランスデューサ１２の受信部１２Ｂによって受信されるように、かつ、第２超音波トランスデューサ１２の送信部１２Ａから送信されて配管Ａの流路Ａ２内を透過した超音波が、配管Ａの壁部Ａ１の内側表面で反射することなく第１超音波トランスデューサ１１の受信部１１Ｂによって受信されるように、第１超音波トランスデューサ１１と第２超音波トランスデューサ１２とが配置され、いわゆるＺ法の流体の流量計測が行われる。
他の例では、第１超音波トランスデューサ１１の送信部１１Ａから送信されて配管Ａの流路Ａ２内を透過した超音波が、配管Ａの壁部Ａ１の内側表面で１回反射した後に配管Ａの流路Ａ２内を透過して第２超音波トランスデューサ１２の受信部１２Ｂによって受信されるように、かつ、第２超音波トランスデューサ１２の送信部１２Ａから送信されて配管Ａの流路Ａ２内を透過した超音波が、配管Ａの壁部Ａ１の内側表面で１回反射した後に配管Ａの流路Ａ２内を透過して第１超音波トランスデューサ１１の受信部１１Ｂによって受信されるように、第１超音波トランスデューサ１１と第２超音波トランスデューサ１２とが配置され、いわゆるＶ法の流体の流量計測が行われてもよい。
更に他の例では、第１超音波トランスデューサ１１の送信部１１Ａから送信されて配管Ａの流路Ａ２内を透過した超音波が、配管Ａの壁部Ａ１の内側表面で２回反射した後に配管Ａの流路Ａ２内を透過して第２超音波トランスデューサ１２の受信部１２Ｂによって受信されるように、かつ、第２超音波トランスデューサ１２の送信部１２Ａから送信されて配管Ａの流路Ａ２内を透過した超音波が、配管Ａの壁部Ａ１の内側表面で２回反射した後に配管Ａの流路Ａ２内を透過して第１超音波トランスデューサ１１の受信部１１Ｂによって受信されるように、第１超音波トランスデューサ１１と第２超音波トランスデューサ１２とが配置され、いわゆるＮ法の流体の流量計測が行われてもよい。
他の例では、第１超音波トランスデューサ１１の送信部１１Ａから送信されて配管Ａの流路Ａ２内を透過した超音波が、配管Ａの壁部Ａ１の内側表面で１回反射した後に配管Ａの流路Ａ２内を透過して第２超音波トランスデューサ１２の受信部１２Ｂによって受信されるように、かつ、第２超音波トランスデューサ１２の送信部１２Ａから送信されて配管Ａの流路Ａ２内を透過した超音波が、配管Ａの壁部Ａ１の内側表面で１回反射した後に配管Ａの流路Ａ２内を透過して第１超音波トランスデューサ１１の受信部１１Ｂによって受信されるように、第１超音波トランスデューサ１１と第２超音波トランスデューサ１２とが配置され、いわゆるＶ法の流体の流量計測が行われてもよい。
更に他の例では、第１超音波トランスデューサ１１の送信部１１Ａから送信されて配管Ａの流路Ａ２内を透過した超音波が、配管Ａの壁部Ａ１の内側表面で２回反射した後に配管Ａの流路Ａ２内を透過して第２超音波トランスデューサ１２の受信部１２Ｂによって受信されるように、かつ、第２超音波トランスデューサ１２の送信部１２Ａから送信されて配管Ａの流路Ａ２内を透過した超音波が、配管Ａの壁部Ａ１の内側表面で２回反射した後に配管Ａの流路Ａ２内を透過して第１超音波トランスデューサ１１の受信部１１Ｂによって受信されるように、第１超音波トランスデューサ１１と第２超音波トランスデューサ１２とが配置され、いわゆるＮ法の流体の流量計測が行われてもよい。
【0022】
また、図１（Ｂ）に示す例では、第１超音波トランスデューサ１１および第２超音波トランスデューサ１２が、配管Ａの壁部Ａ１の外側に配置され、配管Ａの流路Ａ２内を流れる流体に接触しないが、他の例では、第１超音波トランスデューサ１１および第２超音波トランスデューサ１２が、配管Ａの流路Ａ２内を流れる流体に接触するように、配管Ａに対して配置されていてもよい。
【0023】
また、図１（Ｂ）に示す例では、ダンピング材Ｂが、配管Ａの壁部Ａ１の外側に配置されているが、他の例では、ダンピング材Ｂが、配管Ａの壁部Ａ１の外側に配置されていなくてもよい。
【0024】
図１に示す例では、演算部１３が、第１超音波トランスデューサ１１の送信部１１Ａから送信された超音波などを受信した第２超音波トランスデューサ１２の受信部１２Ｂによって出力される受信信号に基づいて、配管Ａの流路Ａ２内を流れる流体の流動様式を判別する等の演算を行う。制御部１４は、第１超音波トランスデューサ１１の制御および第２超音波トランスデューサ１２の制御を実行する。
他の例では、演算部１３が、第２超音波トランスデューサ１２の送信部１２Ａから送信された超音波などを受信した第１超音波トランスデューサ１１の受信部１１Ｂによって出力される受信信号に基づいて、配管Ａの流路Ａ２内を流れる流体の流動様式を判別してもよい。
他の例では、演算部１３が、第２超音波トランスデューサ１２の送信部１２Ａから送信された超音波などを受信した第１超音波トランスデューサ１１の受信部１１Ｂによって出力される受信信号に基づいて、配管Ａの流路Ａ２内を流れる流体の流動様式を判別してもよい。
【0025】
図１に示す例では、配管Ａの流路Ａ２内を流れる流体の流動様式を判別するために、第１超音波トランスデューサ１１の送信部１１Ａが超音波を複数回（例えば２０００回以上）送信し、第２超音波トランスデューサ１２の受信部１２Ｂは、第１超音波トランスデューサ１１の送信部１１Ａから複数回送信された超音波を受信して、複数回分の受信信号（電圧信号）を出力する。
第２超音波トランスデューサ１２の受信部１２Ｂが出力する複数回分の受信信号（電圧信号）は、配管Ａの流路Ａ２内を流れる流体の流量の算出にも用いられる。
更に、配管Ａの流路Ａ２内を流れる流体の流量を算出するために、第２超音波トランスデューサ１２の送信部１２Ａも超音波を複数回（例えば２０００回以上）送信し、第１超音波トランスデューサ１１の受信部１１Ｂが、第２超音波トランスデューサ１２の送信部１２Ａから複数回送信された超音波を受信して、複数回分の受信信号（電圧信号）を出力する。
つまり、演算部１３は、第１超音波トランスデューサ１１の送信部１１Ａから送信された超音波などを受信した第２超音波トランスデューサ１２の受信部１２Ｂによって出力される受信信号と、第２超音波トランスデューサ１２の送信部１２Ａから送信された超音波などを受信した第１超音波トランスデューサ１１の受信部１１Ｂによって出力される受信信号とに基づいて、配管Ａの流路Ａ２内を流れる流体の流量を算出する。
詳細には、配管Ａの流路Ａ２内を流れる流体の流動様式の判別および流量の算出を行うために、演算部１３は、第１受信信号取得部１３Ａと、第２受信信号取得部１３Ｂと、第１標準偏差算出部１３Ｃと、第２標準偏差算出部１３Ｄと、流動様式判別部１３Ｅと、第１超音波受信時刻算出部１３Ｆと、第２超音波受信時刻算出部１３Ｇと、第１超音波伝播時間算出部１３Ｈと、第２超音波伝播時間算出部１３Ｊと、伝播時間差算出部１３Ｋと、流量算出部１３Ｌと、流動様式線図取得部１３Ｍと、湿り度算出部１３Ｎとを備えている。
図１に示す例では、演算部１３が、第１超音波受信時刻算出部１３Ｆと、第２超音波受信時刻算出部１３Ｇと、第１超音波伝播時間算出部１３Ｈと、第２超音波伝播時間算出部１３Ｊとを備えているが、他の例では、演算部１３が、第１超音波受信時刻算出部１３Ｆと、第２超音波受信時刻算出部１３Ｇと、第１超音波伝播時間算出部１３Ｈと、第２超音波伝播時間算出部１３Ｊとを備えていなくてもよい。
第２超音波トランスデューサ１２の受信部１２Ｂが出力する複数回分の受信信号（電圧信号）は、配管Ａの流路Ａ２内を流れる流体の流量の算出にも用いられる。
更に、配管Ａの流路Ａ２内を流れる流体の流量を算出するために、第２超音波トランスデューサ１２の送信部１２Ａも超音波を複数回（例えば２０００回以上）送信し、第１超音波トランスデューサ１１の受信部１１Ｂが、第２超音波トランスデューサ１２の送信部１２Ａから複数回送信された超音波を受信して、複数回分の受信信号（電圧信号）を出力する。
つまり、演算部１３は、第１超音波トランスデューサ１１の送信部１１Ａから送信された超音波などを受信した第２超音波トランスデューサ１２の受信部１２Ｂによって出力される受信信号と、第２超音波トランスデューサ１２の送信部１２Ａから送信された超音波などを受信した第１超音波トランスデューサ１１の受信部１１Ｂによって出力される受信信号とに基づいて、配管Ａの流路Ａ２内を流れる流体の流量を算出する。
詳細には、配管Ａの流路Ａ２内を流れる流体の流動様式の判別および流量の算出を行うために、演算部１３は、第１受信信号取得部１３Ａと、第２受信信号取得部１３Ｂと、第１標準偏差算出部１３Ｃと、第２標準偏差算出部１３Ｄと、流動様式判別部１３Ｅと、第１超音波受信時刻算出部１３Ｆと、第２超音波受信時刻算出部１３Ｇと、第１超音波伝播時間算出部１３Ｈと、第２超音波伝播時間算出部１３Ｊと、伝播時間差算出部１３Ｋと、流量算出部１３Ｌと、流動様式線図取得部１３Ｍと、湿り度算出部１３Ｎとを備えている。
図１に示す例では、演算部１３が、第１超音波受信時刻算出部１３Ｆと、第２超音波受信時刻算出部１３Ｇと、第１超音波伝播時間算出部１３Ｈと、第２超音波伝播時間算出部１３Ｊとを備えているが、他の例では、演算部１３が、第１超音波受信時刻算出部１３Ｆと、第２超音波受信時刻算出部１３Ｇと、第１超音波伝播時間算出部１３Ｈと、第２超音波伝播時間算出部１３Ｊとを備えていなくてもよい。
【0026】
第１受信信号取得部１３Ａは、第２超音波トランスデューサ１２の受信部１２Ｂによって出力される受信信号（つまり、第２超音波トランスデューサ１２の受信部１２Ｂが、第１超音波トランスデューサ１１の送信部１１Ａから送信された超音波などを受信して電圧信号に変換したもの）（例えば図２および図３に「時間波形Ｗ１」で示す受信信号）を取得する。詳細には、第１受信信号取得部１３Ａは、第２超音波トランスデューサ１２の受信部１２Ｂによって出力された複数回分の受信信号（電圧信号）を取得する。
第２受信信号取得部１３Ｂは、第１超音波トランスデューサ１１の受信部１１Ｂによって出力される受信信号（つまり、第１超音波トランスデューサ１１の受信部１１Ｂが、第２超音波トランスデューサ１２の送信部１２Ａから送信された超音波などを受信して電圧信号に変換したもの）（例えば図２および図３に「時間波形Ｗ２」で示す受信信号）を取得する。詳細には、第２受信信号取得部１３Ｂは、第１超音波トランスデューサ１１の受信部１１Ｂによって出力された複数回分の受信信号（電圧信号）を取得する。
第２受信信号取得部１３Ｂは、第１超音波トランスデューサ１１の受信部１１Ｂによって出力される受信信号（つまり、第１超音波トランスデューサ１１の受信部１１Ｂが、第２超音波トランスデューサ１２の送信部１２Ａから送信された超音波などを受信して電圧信号に変換したもの）（例えば図２および図３に「時間波形Ｗ２」で示す受信信号）を取得する。詳細には、第２受信信号取得部１３Ｂは、第１超音波トランスデューサ１１の受信部１１Ｂによって出力された複数回分の受信信号（電圧信号）を取得する。
【0027】
第１標準偏差算出部１３Ｃは、第１超音波トランスデューサ１１の送信部１１Ａが超音波を送信した送信時刻Ｔ０（図２（Ａ）および図３（Ａ）参照）から所定時間経過後の時刻に第２超音波トランスデューサ１２の受信部１２Ｂによって出力される受信信号の振幅（図２（Ａ）および図２（Ｂ）の縦軸）の標準偏差（詳細には、複数回分の受信信号の振幅の標準偏差、分散、および、標準偏差から得られる値のいずれかである標準偏差同等値）（例えば図３（Ａ）および図３（Ｂ）に「時間波形Ｗ３」で示す標準偏差）を算出する。
また、第１標準偏差算出部１３Ｃは、第２超音波トランスデューサ１２の受信部１２Ｂによって出力される受信信号の振幅の標準偏差同等値の時間波形Ｗ３（図３（Ｂ）参照）を描くことができる程度の時間間隔で、複数の時刻における受信信号の振幅の標準偏差同等値を算出する。
また、第１標準偏差算出部１３Ｃは、第２超音波トランスデューサ１２の受信部１２Ｂによって出力される受信信号の振幅の標準偏差同等値の時間波形Ｗ３（図３（Ｂ）参照）を描くことができる程度の時間間隔で、複数の時刻における受信信号の振幅の標準偏差同等値を算出する。
【0028】
詳細には、第１標準偏差算出部１３Ｃは、第１超音波トランスデューサ１１の送信部１１Ａから送信され、配管Ａの壁部Ａ１内を伝わったガイド波としての超音波を、第２超音波トランスデューサ１２の受信部１２Ｂが受信する受信期間（詳細には、第２超音波トランスデューサ１２の受信部１２Ｂが受信すると予測される受信期間）ＴＡ～ＴＢ（図２（Ａ）および図３（Ａ）参照）を含む期間（つまり、受信期間ＴＡ～ＴＢよりも長い期間）Ｔ０～Ｔ３（図２（Ａ）および図３（Ａ）参照）中の各時刻における複数回分の受信信号の標準偏差同等値を算出する。
すなわち、第１標準偏差算出部１３Ｃは、第２超音波トランスデューサ１２の受信部１２Ｂによる第１超音波トランスデューサ１１の送信部１１Ａから送信された超音波の受信期間ＴＡ～ＴＢを含む期間Ｔ０～Ｔ３中の各時刻における複数回分の第２超音波トランスデューサ１２の受信部１２Ｂの受信信号の振幅の標準偏差同等値を算出する。
その結果、第２超音波トランスデューサ１２の受信部１２Ｂによる第１超音波トランスデューサ１１の送信部１１Ａから送信された超音波の受信期間ＴＡ～ＴＢを含む期間Ｔ０～Ｔ３中における第２超音波トランスデューサ１２の受信部１２Ｂの受信信号の振幅の標準偏差同等値の時間波形Ｗ３（図３（Ａ）および図３（Ｂ）参照）が得られる。
つまり、受信期間ＴＡ～ＴＢは、第１超音波トランスデューサ１１の送信部１１Ａによる超音波の送信時刻Ｔ０以降の期間である送信後期間に相当し、第１超音波トランスデューサ１１の送信部１１Ａから送信されて配管Ａの壁部Ａ１内を伝わったガイド波としての超音波が第２超音波トランスデューサ１２の受信部１２Ｂによって受信される期間である。
すなわち、第１標準偏差算出部１３Ｃは、第２超音波トランスデューサ１２の受信部１２Ｂによる第１超音波トランスデューサ１１の送信部１１Ａから送信された超音波の受信期間ＴＡ～ＴＢを含む期間Ｔ０～Ｔ３中の各時刻における複数回分の第２超音波トランスデューサ１２の受信部１２Ｂの受信信号の振幅の標準偏差同等値を算出する。
その結果、第２超音波トランスデューサ１２の受信部１２Ｂによる第１超音波トランスデューサ１１の送信部１１Ａから送信された超音波の受信期間ＴＡ～ＴＢを含む期間Ｔ０～Ｔ３中における第２超音波トランスデューサ１２の受信部１２Ｂの受信信号の振幅の標準偏差同等値の時間波形Ｗ３（図３（Ａ）および図３（Ｂ）参照）が得られる。
つまり、受信期間ＴＡ～ＴＢは、第１超音波トランスデューサ１１の送信部１１Ａによる超音波の送信時刻Ｔ０以降の期間である送信後期間に相当し、第１超音波トランスデューサ１１の送信部１１Ａから送信されて配管Ａの壁部Ａ１内を伝わったガイド波としての超音波が第２超音波トランスデューサ１２の受信部１２Ｂによって受信される期間である。
【0029】
また、第１標準偏差算出部１３Ｃは、第１超音波トランスデューサ１１の送信部１１Ａから送信され、配管Ａの流路Ａ２内を透過した目的信号としての超音波を、第２超音波トランスデューサ１２の受信部１２Ｂが受信する受信期間Ｔ４Ａ～Ｔ４Ｂ（図３（Ｂ）参照）中の各時刻における複数回分の受信信号の標準偏差同等値も算出する。
すなわち、第１標準偏差算出部１３Ｃは、第２超音波トランスデューサ１２の受信部１２Ｂによる第１超音波トランスデューサ１１の送信部１１Ａから送信された超音波の受信期間Ｔ４Ａ～Ｔ４Ｂ中の各時刻における複数回分の第２超音波トランスデューサ１２の受信部１２Ｂの受信信号の振幅の標準偏差同等値も算出する。
その結果、第２超音波トランスデューサ１２の受信部１２Ｂによる第１超音波トランスデューサ１１の送信部１１Ａから送信された超音波の受信期間Ｔ４Ａ～Ｔ４Ｂ中における第２超音波トランスデューサ１２の受信部１２Ｂの受信信号の振幅の標準偏差同等値の時間波形Ｗ３（図３（Ｂ）参照）も得られる。
すなわち、第１標準偏差算出部１３Ｃは、第２超音波トランスデューサ１２の受信部１２Ｂによる第１超音波トランスデューサ１１の送信部１１Ａから送信された超音波の受信期間Ｔ４Ａ～Ｔ４Ｂ中の各時刻における複数回分の第２超音波トランスデューサ１２の受信部１２Ｂの受信信号の振幅の標準偏差同等値も算出する。
その結果、第２超音波トランスデューサ１２の受信部１２Ｂによる第１超音波トランスデューサ１１の送信部１１Ａから送信された超音波の受信期間Ｔ４Ａ～Ｔ４Ｂ中における第２超音波トランスデューサ１２の受信部１２Ｂの受信信号の振幅の標準偏差同等値の時間波形Ｗ３（図３（Ｂ）参照）も得られる。
【0030】
第２標準偏差算出部１３Ｄは、第２超音波トランスデューサ１２の送信部１２Ａが超音波を送信した送信時刻Ｔ０（図２（Ａ）および図３（Ａ）参照）から所定時間経過後の時刻に第１超音波トランスデューサ１１の受信部１１Ｂによって出力される受信信号の振幅（図２（Ａ）および図２（Ｂ）の縦軸）の標準偏差（詳細には、複数回分の受信信号の振幅の標準偏差、分散、および、標準偏差から得られる値のいずれかである標準偏差同等値）（例えば図３（Ａ）および図３（Ｂ）に「時間波形Ｗ４」で示す標準偏差）を算出する。
また、第２標準偏差算出部１３Ｄは、第１超音波トランスデューサ１１の受信部１１Ｂによって出力される受信信号の振幅の標準偏差同等値の時間波形Ｗ４（図３（Ｂ）参照）を描くことができる程度の時間間隔で、複数の時刻における受信信号の振幅の標準偏差同等値を算出する。
詳細には、第２標準偏差算出部１３Ｄは、第２超音波トランスデューサ１２の送信部１２Ａから送信され、配管Ａの流路Ａ２内を透過した目的信号としての超音波を、第１超音波トランスデューサ１１の受信部１１Ｂが受信する受信期間（詳細には、第１超音波トランスデューサ１１の受信部１１Ｂが受信すると予測される受信期間）Ｔ５Ａ～Ｔ５Ｂ（図３（Ｂ）参照）を含む期間（つまり、受信期間Ｔ５Ａ～Ｔ５Ｂよりも長い期間）Ｔ０～Ｔ３（図２（Ａ）および図３（Ａ）参照）中の各時刻における複数回分の受信信号の標準偏差同等値を算出する。
すなわち、第２標準偏差算出部１３Ｄは、第１超音波トランスデューサ１１の受信部１１Ｂによる第２超音波トランスデューサ１２の送信部１２Ａから送信された超音波の受信期間Ｔ５Ａ～Ｔ５Ｂを含む期間Ｔ０～Ｔ３中の各時刻における複数回分の第１超音波トランスデューサ１１の受信部１１Ｂの受信信号の振幅の標準偏差同等値を算出する。
その結果、第１超音波トランスデューサ１１の受信部１１Ｂによる第２超音波トランスデューサ１２の送信部１２Ａから送信された超音波の受信期間Ｔ５Ａ～Ｔ５Ｂを含む期間Ｔ０～Ｔ３中における第１超音波トランスデューサ１１の受信部１１Ｂの受信信号の振幅の標準偏差同等値の時間波形Ｗ４（図３（Ａ）および図３（Ｂ）参照）が得られる。
また、第２標準偏差算出部１３Ｄは、第１超音波トランスデューサ１１の受信部１１Ｂによって出力される受信信号の振幅の標準偏差同等値の時間波形Ｗ４（図３（Ｂ）参照）を描くことができる程度の時間間隔で、複数の時刻における受信信号の振幅の標準偏差同等値を算出する。
詳細には、第２標準偏差算出部１３Ｄは、第２超音波トランスデューサ１２の送信部１２Ａから送信され、配管Ａの流路Ａ２内を透過した目的信号としての超音波を、第１超音波トランスデューサ１１の受信部１１Ｂが受信する受信期間（詳細には、第１超音波トランスデューサ１１の受信部１１Ｂが受信すると予測される受信期間）Ｔ５Ａ～Ｔ５Ｂ（図３（Ｂ）参照）を含む期間（つまり、受信期間Ｔ５Ａ～Ｔ５Ｂよりも長い期間）Ｔ０～Ｔ３（図２（Ａ）および図３（Ａ）参照）中の各時刻における複数回分の受信信号の標準偏差同等値を算出する。
すなわち、第２標準偏差算出部１３Ｄは、第１超音波トランスデューサ１１の受信部１１Ｂによる第２超音波トランスデューサ１２の送信部１２Ａから送信された超音波の受信期間Ｔ５Ａ～Ｔ５Ｂを含む期間Ｔ０～Ｔ３中の各時刻における複数回分の第１超音波トランスデューサ１１の受信部１１Ｂの受信信号の振幅の標準偏差同等値を算出する。
その結果、第１超音波トランスデューサ１１の受信部１１Ｂによる第２超音波トランスデューサ１２の送信部１２Ａから送信された超音波の受信期間Ｔ５Ａ～Ｔ５Ｂを含む期間Ｔ０～Ｔ３中における第１超音波トランスデューサ１１の受信部１１Ｂの受信信号の振幅の標準偏差同等値の時間波形Ｗ４（図３（Ａ）および図３（Ｂ）参照）が得られる。
【0031】
流動様式判別部１３Ｅは、受信期間ＴＡ～ＴＢ中の各時刻における複数回分の第２超音波トランスデューサ１２の受信部１２Ｂの受信信号の振幅の標準偏差同等値（例えば図３（Ａ）、図４（Ａ）、図６（Ａ）、図７（Ａ）および図８（Ａ）に「時間波形Ｗ３」で示す標準偏差）に基づいて、配管Ａの流路Ａ２内を流れる流体の流動様式を判別する。
詳細には、流動様式判別部１３Ｅは、受信期間ＴＡ～ＴＢ中の各時刻における複数回分の第２超音波トランスデューサ１２の受信部１２Ｂの受信信号の振幅の標準偏差同等値に基づいて、配管Ａの流路Ａ２内を流れる流体の流動様式が、例えば、乾き蒸気、層状流、波状流、波状環状遷移および環状噴霧流のいずれであるかを判別する。
詳細には、流動様式判別部１３Ｅは、受信期間ＴＡ～ＴＢ中の各時刻における複数回分の第２超音波トランスデューサ１２の受信部１２Ｂの受信信号の振幅の標準偏差同等値に基づいて、配管Ａの流路Ａ２内を流れる流体の流動様式が、例えば、乾き蒸気、層状流、波状流、波状環状遷移および環状噴霧流のいずれであるかを判別する。
【0032】
第１超音波受信時刻算出部１３Ｆは、第１標準偏差算出部１３Ｃによって算出された第２超音波トランスデューサ１２の受信部１２Ｂによる第１超音波トランスデューサ１１の送信部１１Ａから送信された超音波の受信期間Ｔ４Ａ～Ｔ４Ｂ（図３（Ｂ）参照）を含む期間Ｔ０～Ｔ３（図３（Ａ）参照）中の各時刻における複数回分の第２超音波トランスデューサ１２の受信部１２Ｂの受信信号の振幅の標準偏差同等値（例えば図３（Ｂ）に「時間波形Ｗ３」で示す標準偏差）に基づいて、第２超音波トランスデューサ１２の受信部１２Ｂによる第１超音波トランスデューサ１１の送信部１１Ａから送信された超音波の受信時刻（例えば時刻Ｔ４Ａ（図３（Ｂ）参照））を算出する。
詳細には、第１超音波受信時刻算出部１３Ｆは、第２超音波トランスデューサ１２の受信部１２Ｂによる第１超音波トランスデューサ１１の送信部１１Ａから送信された超音波の受信期間Ｔ４Ａ～Ｔ４Ｂを含む期間Ｔ０～Ｔ３中における第２超音波トランスデューサ１２の受信部１２Ｂの受信信号の振幅の標準偏差同等値の時間波形Ｗ３（図３（Ａ）および図３（Ｂ）参照）に基づいて、第２超音波トランスデューサ１２の受信部１２Ｂによる第１超音波トランスデューサ１１の送信部１１Ａから送信された超音波の受信時刻（例えば時刻Ｔ４Ａ）を算出する。
他の例（第１超音波トランスデューサ１１の送信部１１Ａから送信された超音波の受信期間Ｔ４Ａ～Ｔ４Ｂ中における第２超音波トランスデューサ１２の受信部１２Ｂの受信信号が、ノイズに埋もれていない例）では、第１超音波受信時刻算出部１３Ｆが、第２超音波トランスデューサ１２の受信部１２Ｂによる第１超音波トランスデューサ１１の送信部１１Ａから送信された超音波の受信期間Ｔ４Ａ～Ｔ４Ｂ中における第２超音波トランスデューサ１２の受信部１２Ｂの受信信号に基づいて、第２超音波トランスデューサ１２の受信部１２Ｂによる第１超音波トランスデューサ１１の送信部１１Ａから送信された超音波の受信時刻（例えば時刻Ｔ４Ａ）を算出してもよい。
詳細には、第１超音波受信時刻算出部１３Ｆは、第２超音波トランスデューサ１２の受信部１２Ｂによる第１超音波トランスデューサ１１の送信部１１Ａから送信された超音波の受信期間Ｔ４Ａ～Ｔ４Ｂを含む期間Ｔ０～Ｔ３中における第２超音波トランスデューサ１２の受信部１２Ｂの受信信号の振幅の標準偏差同等値の時間波形Ｗ３（図３（Ａ）および図３（Ｂ）参照）に基づいて、第２超音波トランスデューサ１２の受信部１２Ｂによる第１超音波トランスデューサ１１の送信部１１Ａから送信された超音波の受信時刻（例えば時刻Ｔ４Ａ）を算出する。
他の例（第１超音波トランスデューサ１１の送信部１１Ａから送信された超音波の受信期間Ｔ４Ａ～Ｔ４Ｂ中における第２超音波トランスデューサ１２の受信部１２Ｂの受信信号が、ノイズに埋もれていない例）では、第１超音波受信時刻算出部１３Ｆが、第２超音波トランスデューサ１２の受信部１２Ｂによる第１超音波トランスデューサ１１の送信部１１Ａから送信された超音波の受信期間Ｔ４Ａ～Ｔ４Ｂ中における第２超音波トランスデューサ１２の受信部１２Ｂの受信信号に基づいて、第２超音波トランスデューサ１２の受信部１２Ｂによる第１超音波トランスデューサ１１の送信部１１Ａから送信された超音波の受信時刻（例えば時刻Ｔ４Ａ）を算出してもよい。
【0033】
第２超音波受信時刻算出部１３Ｇは、第２標準偏差算出部１３Ｄによって算出された第１超音波トランスデューサ１１の受信部１１Ｂによる第２超音波トランスデューサ１２の送信部１２Ａから送信された超音波の受信期間Ｔ５Ａ～Ｔ５Ｂ（図３（Ｂ）参照）を含む期間Ｔ０～Ｔ３（図３（Ａ）参照）中の各時刻における複数回分の第１超音波トランスデューサ１１の受信部１１Ｂの受信信号の振幅の標準偏差同等値（例えば図３（Ｂ）に「時間波形Ｗ４」で示す標準偏差）に基づいて、第１超音波トランスデューサ１１の受信部１１Ｂによる第２超音波トランスデューサ１２の送信部１２Ａから送信された超音波の受信時刻（例えば時刻Ｔ５Ａ（図３（Ｂ）参照））を算出する。
詳細には、第２超音波受信時刻算出部１３Ｇは、第１超音波トランスデューサ１１の受信部１１Ｂによる第２超音波トランスデューサ１２の送信部１２Ａから送信された超音波の受信期間Ｔ５Ａ～Ｔ５Ｂを含む期間Ｔ０～Ｔ３中における第１超音波トランスデューサ１１の受信部１１Ｂの受信信号の振幅の標準偏差同等値の時間波形Ｗ４（図３（Ａ）および図３（Ｂ）参照）に基づいて、第１超音波トランスデューサ１１の受信部１１Ｂによる第２超音波トランスデューサ１２の送信部１２Ａから送信された超音波の受信時刻（例えば時刻Ｔ５Ａ）を算出する。
他の例（第２超音波トランスデューサ１２の送信部１２Ａから送信された超音波の受信期間Ｔ５Ａ～Ｔ５Ｂ中における第１超音波トランスデューサ１１の受信部１１Ｂの受信信号が、ノイズに埋もれていない例）では、第２超音波受信時刻算出部１３Ｇが、第１超音波トランスデューサ１１の受信部１１Ｂによる第２超音波トランスデューサ１２の送信部１２Ａから送信された超音波の受信期間Ｔ５Ａ～Ｔ５Ｂ中における第１超音波トランスデューサ１１の受信部１１Ｂの受信信号に基づいて、第１超音波トランスデューサ１１の受信部１１Ｂによる第２超音波トランスデューサ１２の送信部１２Ａから送信された超音波の受信時刻（例えば時刻Ｔ５Ａ）を算出してもよい。
詳細には、第２超音波受信時刻算出部１３Ｇは、第１超音波トランスデューサ１１の受信部１１Ｂによる第２超音波トランスデューサ１２の送信部１２Ａから送信された超音波の受信期間Ｔ５Ａ～Ｔ５Ｂを含む期間Ｔ０～Ｔ３中における第１超音波トランスデューサ１１の受信部１１Ｂの受信信号の振幅の標準偏差同等値の時間波形Ｗ４（図３（Ａ）および図３（Ｂ）参照）に基づいて、第１超音波トランスデューサ１１の受信部１１Ｂによる第２超音波トランスデューサ１２の送信部１２Ａから送信された超音波の受信時刻（例えば時刻Ｔ５Ａ）を算出する。
他の例（第２超音波トランスデューサ１２の送信部１２Ａから送信された超音波の受信期間Ｔ５Ａ～Ｔ５Ｂ中における第１超音波トランスデューサ１１の受信部１１Ｂの受信信号が、ノイズに埋もれていない例）では、第２超音波受信時刻算出部１３Ｇが、第１超音波トランスデューサ１１の受信部１１Ｂによる第２超音波トランスデューサ１２の送信部１２Ａから送信された超音波の受信期間Ｔ５Ａ～Ｔ５Ｂ中における第１超音波トランスデューサ１１の受信部１１Ｂの受信信号に基づいて、第１超音波トランスデューサ１１の受信部１１Ｂによる第２超音波トランスデューサ１２の送信部１２Ａから送信された超音波の受信時刻（例えば時刻Ｔ５Ａ）を算出してもよい。
【0034】
第１超音波伝播時間算出部１３Ｈは、第１超音波受信時刻算出部１３Ｆによって算出された第２超音波トランスデューサ１２の受信部１２Ｂによる第１超音波トランスデューサ１１の送信部１１Ａから送信された超音波の受信時刻（例えば時刻Ｔ４Ａ（図３（Ｂ）参照））と、第１超音波トランスデューサ１１の送信部１１Ａによる超音波の送信時刻Ｔ０（図２（Ａ）および図３（Ａ）参照）とに基づいて、第１超音波トランスデューサ１１の送信部１１Ａから送信された超音波の伝播時間（例えば伝播時間（Ｔ４Ａ－Ｔ０））を算出する。
第２超音波伝播時間算出部１３Ｊは、第２超音波受信時刻算出部１３Ｇによって算出された第１超音波トランスデューサ１１の受信部１１Ｂによる第２超音波トランスデューサ１２の送信部１２Ａから送信された超音波の受信時刻（例えば時刻Ｔ５Ａ（図３（Ｂ）参照））と、第２超音波トランスデューサ１２の送信部１２Ａによる超音波の送信時刻Ｔ０（図２（Ａ）および図３（Ａ）参照）とに基づいて、第２超音波トランスデューサ１２の送信部１２Ａから送信された超音波の伝播時間（例えば伝播時間（Ｔ５Ａ－Ｔ０））を算出する。
第２超音波伝播時間算出部１３Ｊは、第２超音波受信時刻算出部１３Ｇによって算出された第１超音波トランスデューサ１１の受信部１１Ｂによる第２超音波トランスデューサ１２の送信部１２Ａから送信された超音波の受信時刻（例えば時刻Ｔ５Ａ（図３（Ｂ）参照））と、第２超音波トランスデューサ１２の送信部１２Ａによる超音波の送信時刻Ｔ０（図２（Ａ）および図３（Ａ）参照）とに基づいて、第２超音波トランスデューサ１２の送信部１２Ａから送信された超音波の伝播時間（例えば伝播時間（Ｔ５Ａ－Ｔ０））を算出する。
【0035】
伝播時間差算出部１３Ｋは、第１超音波伝播時間算出部１３Ｈによって算出された第１超音波トランスデューサ１１の送信部１１Ａから送信された超音波の伝播時間（例えば伝播時間（Ｔ４Ａ－Ｔ０））と、第２超音波伝播時間算出部１３Ｊによって算出された第２超音波トランスデューサ１２の送信部１２Ａから送信された超音波の伝播時間（例えば伝播時間（Ｔ５Ａ－Ｔ０））との差分である伝播時間差（例えば伝播時間差（Ｔ５Ａ－Ｔ４Ａ））を算出する。
つまり、伝播時間差算出部１３Ｋは、第２超音波トランスデューサ１２の受信部１２Ｂの複数回分の受信信号の振幅の標準偏差（例えば図３（Ｂ）に「時間波形Ｗ３」で示す標準偏差）と、第１超音波トランスデューサ１１の受信部１１Ｂの複数回分の受信信号の振幅の標準偏差（例えば図３（Ｂ）に「時間波形Ｗ４」で示す標準偏差）とに基づいて、第１超音波トランスデューサ１１の送信部１１Ａから送信された超音波の伝播時間と第２超音波トランスデューサ１２の送信部１２Ａから送信された超音波の伝播時間との差分である伝播時間差（Ｔ５Ａ－Ｔ４Ａ）を算出する。
他の例では、伝播時間差算出部１３Ｋが、第２超音波トランスデューサ１２の受信部１２Ｂによる第１超音波トランスデューサ１１の送信部１１Ａからの超音波の受信時刻（例えば時刻Ｔ４Ａ）と、第１超音波トランスデューサ１１の受信部１１Ｂによる第２超音波トランスデューサ１２の送信部１２Ａからの超音波の受信時刻（例えば時刻Ｔ５Ａ）とを用いることなく、伝播時間差（例えば伝播時間差（Ｔ５Ａ－Ｔ４Ａ））を算出してもよい。
具体的には、伝播時間差算出部１３Ｋは、第２超音波トランスデューサ１２の受信部１２Ｂによって出力される受信信号の振幅の標準偏差同等値の時間波形と、第１超音波トランスデューサ１１の受信部１１Ｂによって出力される受信信号の振幅の標準偏差同等値の時間波形との相関処理（例えば特開２０１６－１８０６７９号公報に記載されているような公知の相関処理）を実行することによって、伝播時間差（例えば伝播時間差（Ｔ５Ａ－Ｔ４Ａ））を算出することができる。
更に他の例（第１超音波トランスデューサ１１の送信部１１Ａから送信された超音波の受信期間Ｔ４Ａ～Ｔ４Ｂ中における第２超音波トランスデューサ１２の受信部１２Ｂの受信信号が、ノイズに埋もれておらず、第２超音波トランスデューサ１２の送信部１２Ａから送信された超音波の受信期間Ｔ５Ａ～Ｔ５Ｂ中における第１超音波トランスデューサ１１の受信部１１Ｂの受信信号が、ノイズに埋もれていない例）では、伝播時間差算出部１３Ｋが、第２超音波トランスデューサ１２の受信部１２Ｂによって出力される受信信号の時間波形と、第１超音波トランスデューサ１１の受信部１１Ｂによって出力される受信信号の時間波形とに基づいて、伝播時間差（例えば伝播時間差（Ｔ５Ａ－Ｔ４Ａ））を算出してもよい。
つまり、伝播時間差算出部１３Ｋは、第２超音波トランスデューサ１２の受信部１２Ｂの複数回分の受信信号の振幅の標準偏差（例えば図３（Ｂ）に「時間波形Ｗ３」で示す標準偏差）と、第１超音波トランスデューサ１１の受信部１１Ｂの複数回分の受信信号の振幅の標準偏差（例えば図３（Ｂ）に「時間波形Ｗ４」で示す標準偏差）とに基づいて、第１超音波トランスデューサ１１の送信部１１Ａから送信された超音波の伝播時間と第２超音波トランスデューサ１２の送信部１２Ａから送信された超音波の伝播時間との差分である伝播時間差（Ｔ５Ａ－Ｔ４Ａ）を算出する。
他の例では、伝播時間差算出部１３Ｋが、第２超音波トランスデューサ１２の受信部１２Ｂによる第１超音波トランスデューサ１１の送信部１１Ａからの超音波の受信時刻（例えば時刻Ｔ４Ａ）と、第１超音波トランスデューサ１１の受信部１１Ｂによる第２超音波トランスデューサ１２の送信部１２Ａからの超音波の受信時刻（例えば時刻Ｔ５Ａ）とを用いることなく、伝播時間差（例えば伝播時間差（Ｔ５Ａ－Ｔ４Ａ））を算出してもよい。
具体的には、伝播時間差算出部１３Ｋは、第２超音波トランスデューサ１２の受信部１２Ｂによって出力される受信信号の振幅の標準偏差同等値の時間波形と、第１超音波トランスデューサ１１の受信部１１Ｂによって出力される受信信号の振幅の標準偏差同等値の時間波形との相関処理（例えば特開２０１６－１８０６７９号公報に記載されているような公知の相関処理）を実行することによって、伝播時間差（例えば伝播時間差（Ｔ５Ａ－Ｔ４Ａ））を算出することができる。
更に他の例（第１超音波トランスデューサ１１の送信部１１Ａから送信された超音波の受信期間Ｔ４Ａ～Ｔ４Ｂ中における第２超音波トランスデューサ１２の受信部１２Ｂの受信信号が、ノイズに埋もれておらず、第２超音波トランスデューサ１２の送信部１２Ａから送信された超音波の受信期間Ｔ５Ａ～Ｔ５Ｂ中における第１超音波トランスデューサ１１の受信部１１Ｂの受信信号が、ノイズに埋もれていない例）では、伝播時間差算出部１３Ｋが、第２超音波トランスデューサ１２の受信部１２Ｂによって出力される受信信号の時間波形と、第１超音波トランスデューサ１１の受信部１１Ｂによって出力される受信信号の時間波形とに基づいて、伝播時間差（例えば伝播時間差（Ｔ５Ａ－Ｔ４Ａ））を算出してもよい。
【0036】
流量算出部１３Ｌは、伝播時間差算出部１３Ｋによって算出された伝播時間差（例えば伝播時間差（Ｔ５Ａ－Ｔ４Ａ））に基づいて、配管Ａの流路Ａ２内を流れる流体の流量を算出する。詳細には、流量算出部１３Ｌは、伝播時間差算出部１３Ｋによって算出された伝播時間差、音速、配管Ａの内径、角度θ（図１（Ｂ）参照）、公知の流量算出式などに基づいて、配管Ａの流路Ａ２内を流れる流体の流量を算出する。
【0037】
流動様式線図取得部１３Ｍは、配管Ａの流路Ａ２内を流れる流体の流動様式と配管Ａの流路Ａ２内の気相容積流束ＪＧと配管Ａの流路Ａ２内の液相容積流束ＪＬとの関係を示す流動様式線図（例えば図５参照）を取得する。流動様式線図は、第１実施形態の流動様式判別装置１によって、あるいは、流動様式判別装置１の外部において予め作成（算出）される。
湿り度算出部１３Ｎは、流動様式線図取得部１３Ｍによって取得された流動様式線図と、流動様式判別部１３Ｅによって判別された配管Ａの流路Ａ２内を流れる流体の流動様式と、流量算出部１３Ｌによって算出された配管Ａの流路Ａ２内を流れる流体の流量とに基づいて、配管Ａの流路Ａ２内を流れる流体の湿り度［％］を算出する。
湿り度算出部１３Ｎは、流動様式線図取得部１３Ｍによって取得された流動様式線図と、流動様式判別部１３Ｅによって判別された配管Ａの流路Ａ２内を流れる流体の流動様式と、流量算出部１３Ｌによって算出された配管Ａの流路Ａ２内を流れる流体の流量とに基づいて、配管Ａの流路Ａ２内を流れる流体の湿り度［％］を算出する。
【0038】
図２は配管Ａの流路Ａ２内の流体の流れが水を含まない空気の単相流である場合（配管Ａの流路Ａ２内を流れる流体の流動様式が乾き蒸気である場合）に第１超音波トランスデューサ１１の送信部１１Ａから送信された超音波などを受信した第２超音波トランスデューサ１２の受信部１２Ｂによって出力される受信信号の時間波形Ｗ１、および、第２超音波トランスデューサ１２の送信部１２Ａから送信された超音波などを受信した第１超音波トランスデューサ１１の受信部１１Ｂによって出力される受信信号の時間波形Ｗ２の一例を示す図である。詳細には、図２（Ａ）は時間波形Ｗ１および時間波形Ｗ２の全体を示しており、図２（Ｂ）は図２（Ａ）の一部を拡大して示している。図２（Ａ）および図２（Ｂ）の縦軸は電圧［Ｖ］を示しており、図２（Ａ）および図２（Ｂ）の横軸は時刻［μｓ］を示している。
図３は配管Ａの流路Ａ２内の流体の流れが水を含まない空気の単相流である場合（配管Ａの流路Ａ２内を流れる流体の流動様式が乾き蒸気である場合）に演算部１３の第１標準偏差算出部１３Ｃによって算出された第２超音波トランスデューサ１２の受信部１２Ｂの受信信号の振幅の標準偏差同等値の時間波形Ｗ３と、演算部１３の第２標準偏差算出部１３Ｄによって算出された第１超音波トランスデューサ１１の受信部１１Ｂの受信信号の振幅の標準偏差同等値の時間波形Ｗ４とを図２に追記した図である。図３（Ａ）および図３（Ｂ）の縦軸は電圧［Ｖ］および標準偏差［Ｖ］を示しており、図３（Ａ）および図３（Ｂ）の横軸は時刻［μｓ］を示している。
図３は配管Ａの流路Ａ２内の流体の流れが水を含まない空気の単相流である場合（配管Ａの流路Ａ２内を流れる流体の流動様式が乾き蒸気である場合）に演算部１３の第１標準偏差算出部１３Ｃによって算出された第２超音波トランスデューサ１２の受信部１２Ｂの受信信号の振幅の標準偏差同等値の時間波形Ｗ３と、演算部１３の第２標準偏差算出部１３Ｄによって算出された第１超音波トランスデューサ１１の受信部１１Ｂの受信信号の振幅の標準偏差同等値の時間波形Ｗ４とを図２に追記した図である。図３（Ａ）および図３（Ｂ）の縦軸は電圧［Ｖ］および標準偏差［Ｖ］を示しており、図３（Ａ）および図３（Ｂ）の横軸は時刻［μｓ］を示している。
【0039】
図２および図３において、「単相流：ｕ＝２９．１ｍ／ｓ」は、配管Ａの流路Ａ２内の流体の流れが、水を含まない空気の単相流であって、配管Ａの流路Ａ２内における空気の平均流速が２９．１［ｍ／ｓ］であることを示している。
図２（Ａ）および図２（Ｂ）の時間波形Ｗ１は、第１超音波トランスデューサ１１の送信部１１Ａが超音波を例えば５０００回送信したときに第２超音波トランスデューサ１２の受信部１２Ｂによって出力される受信信号の時間波形の平均を示している。また、第１超音波トランスデューサ１１の送信部１１Ａによる超音波の送信時刻が、図２（Ａ）の横軸の時刻Ｔ０に設定されている。
図２（Ａ）および図２（Ｂ）の時間波形Ｗ２は、第２超音波トランスデューサ１２の送信部１２Ａが超音波を例えば５０００回送信したときに第１超音波トランスデューサ１１の受信部１１Ｂによって出力される受信信号の時間波形の平均を示している。また、第２超音波トランスデューサ１２の送信部１２Ａによる超音波の送信時刻が、図２（Ａ）の横軸の時刻Ｔ０に設定されている。
図２（Ａ）および図２（Ｂ）の時間波形Ｗ１は、第１超音波トランスデューサ１１の送信部１１Ａが超音波を例えば５０００回送信したときに第２超音波トランスデューサ１２の受信部１２Ｂによって出力される受信信号の時間波形の平均を示している。また、第１超音波トランスデューサ１１の送信部１１Ａによる超音波の送信時刻が、図２（Ａ）の横軸の時刻Ｔ０に設定されている。
図２（Ａ）および図２（Ｂ）の時間波形Ｗ２は、第２超音波トランスデューサ１２の送信部１２Ａが超音波を例えば５０００回送信したときに第１超音波トランスデューサ１１の受信部１１Ｂによって出力される受信信号の時間波形の平均を示している。また、第２超音波トランスデューサ１２の送信部１２Ａによる超音波の送信時刻が、図２（Ａ）の横軸の時刻Ｔ０に設定されている。
【0040】
図２および図３に示す例では、演算部１３の第１受信信号取得部１３Ａが、第２超音波トランスデューサ１２の受信部１２Ｂによって出力される受信信号（例えば図２（Ａ）および図２（Ｂ）に時間波形Ｗ１で示す受信信号）を取得する。
詳細には、第１受信信号取得部１３Ａは、第１超音波トランスデューサ１１の送信部１１Ａから例えば５０００回送信された超音波などを受信した第２超音波トランスデューサ１２の受信部１２Ｂによって出力された例えば５０００回分の受信信号の時間波形の平均（図２および図３に「Ｗ１」で示す時間波形）を取得する。
演算部１３の第１標準偏差算出部１３Ｃは、第２超音波トランスデューサ１２の受信部１２Ｂによって出力された例えば５０００回分の受信信号のそれぞれの時間波形（図示せず）の振幅の標準偏差同等値を算出する。
図３（Ａ）および図３（Ｂ）の時間波形Ｗ３は、第１標準偏差算出部１３Ｃによって算出された第２超音波トランスデューサ１２の受信部１２Ｂの例えば５０００回分の受信信号の振幅の標準偏差の時間波形を示している。
つまり、図３（Ａ）および図３（Ｂ）の時間波形Ｗ３は、第２超音波トランスデューサ１２の受信部１２Ｂによって出力された例えば５０００回分の受信信号のそれぞれの時間波形（図示せず）の振幅の標準偏差の時間波形を示している。
詳細には、第１受信信号取得部１３Ａは、第１超音波トランスデューサ１１の送信部１１Ａから例えば５０００回送信された超音波などを受信した第２超音波トランスデューサ１２の受信部１２Ｂによって出力された例えば５０００回分の受信信号の時間波形の平均（図２および図３に「Ｗ１」で示す時間波形）を取得する。
演算部１３の第１標準偏差算出部１３Ｃは、第２超音波トランスデューサ１２の受信部１２Ｂによって出力された例えば５０００回分の受信信号のそれぞれの時間波形（図示せず）の振幅の標準偏差同等値を算出する。
図３（Ａ）および図３（Ｂ）の時間波形Ｗ３は、第１標準偏差算出部１３Ｃによって算出された第２超音波トランスデューサ１２の受信部１２Ｂの例えば５０００回分の受信信号の振幅の標準偏差の時間波形を示している。
つまり、図３（Ａ）および図３（Ｂ）の時間波形Ｗ３は、第２超音波トランスデューサ１２の受信部１２Ｂによって出力された例えば５０００回分の受信信号のそれぞれの時間波形（図示せず）の振幅の標準偏差の時間波形を示している。
【0041】
図２および図３に示す例では、受信期間ＴＡ～ＴＢ中の各時刻における複数回分の第２超音波トランスデューサ１２の受信部１２Ｂの受信信号の振幅の標準偏差同等値（図３（Ａ）に時間波形Ｗ３で示す標準偏差同等値）がほぼ一定値である。そのため、演算部１３の流動様式判別部１３Ｅは、受信期間ＴＡ～ＴＢ中の各時刻における複数回分の第２超音波トランスデューサ１２の受信部１２Ｂの受信信号の振幅の標準偏差同等値の時間波形Ｗ３（図３（Ａ）参照）に基づいて、配管Ａの流路Ａ２内の流体の流れが、水を含まない空気の単相流である（配管Ａの流路Ａ２内を流れる流体の流動様式が乾き蒸気である）と判別する。
【0042】
更に、図２および図３に示す例では、演算部１３の第１超音波受信時刻算出部１３Ｆが、図３（Ａ）および図３（Ｂ）に示す第２超音波トランスデューサ１２の受信部１２Ｂの複数回分の受信信号の振幅の標準偏差同等値の時間波形Ｗ３に基づいて、第２超音波トランスデューサ１２の受信部１２Ｂによる第１超音波トランスデューサ１１の送信部１１Ａから送信された超音波の受信時刻（例えば時刻Ｔ４Ａ（図３（Ｂ）参照））を算出する。
【0043】
また、図２および図３に示す例では、演算部１３の第２受信信号取得部１３Ｂが、第１超音波トランスデューサ１１の受信部１１Ｂによって出力される受信信号（例えば図２（Ａ）および図２（Ｂ）に時間波形Ｗ２で示す受信信号）を取得する。
詳細には、第２受信信号取得部１３Ｂは、第２超音波トランスデューサ１２の送信部１２Ａから例えば５０００回送信された超音波などを受信した第１超音波トランスデューサ１１の受信部１１Ｂによって出力された例えば５０００回分の受信信号の時間波形の平均（図２および図３に「Ｗ２」で示す時間波形）を取得する。
演算部１３の第２標準偏差算出部１３Ｄは、第１超音波トランスデューサ１１の受信部１１Ｂによって出力された例えば５０００回分の受信信号のそれぞれの時間波形（図示せず）の振幅の標準偏差同等値を算出する。
図３（Ａ）および図３（Ｂ）の時間波形Ｗ４は、第２標準偏差算出部１３Ｄによって算出された第１超音波トランスデューサ１１の受信部１１Ｂの例えば５０００回分の受信信号の振幅の標準偏差の時間波形を示している。
つまり、図３（Ａ）および図３（Ｂ）の時間波形Ｗ４は、第１超音波トランスデューサ１１の受信部１１Ｂによって出力された例えば５０００回分の受信信号のそれぞれの時間波形（図示せず）の振幅の標準偏差の時間波形を示している。
詳細には、第２受信信号取得部１３Ｂは、第２超音波トランスデューサ１２の送信部１２Ａから例えば５０００回送信された超音波などを受信した第１超音波トランスデューサ１１の受信部１１Ｂによって出力された例えば５０００回分の受信信号の時間波形の平均（図２および図３に「Ｗ２」で示す時間波形）を取得する。
演算部１３の第２標準偏差算出部１３Ｄは、第１超音波トランスデューサ１１の受信部１１Ｂによって出力された例えば５０００回分の受信信号のそれぞれの時間波形（図示せず）の振幅の標準偏差同等値を算出する。
図３（Ａ）および図３（Ｂ）の時間波形Ｗ４は、第２標準偏差算出部１３Ｄによって算出された第１超音波トランスデューサ１１の受信部１１Ｂの例えば５０００回分の受信信号の振幅の標準偏差の時間波形を示している。
つまり、図３（Ａ）および図３（Ｂ）の時間波形Ｗ４は、第１超音波トランスデューサ１１の受信部１１Ｂによって出力された例えば５０００回分の受信信号のそれぞれの時間波形（図示せず）の振幅の標準偏差の時間波形を示している。
【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】
【0086】
図５および図８に示す例では、演算部１３の湿り度算出部１３Ｎが、流動様式線図取得部１３Ｍによって取得された流動様式線図（図５参照）と、流動様式判別部１３Ｅによって判別された配管Ａの流路Ａ２内を流れる流体の流動様式「環状噴霧流」と、流量算出部１３Ｌによって算出された配管Ａの流路Ａ２内を流れる流体の流量とに基づいて、配管Ａの流路Ａ２内を流れる流体の湿り度［％］を算出する。
詳細には、湿り度算出部１３Ｎは、流量算出部１３Ｌによって算出された配管Ａ内を流れる流体の流量を、配管Ａの流路Ａ２内の気相容積流束ＪＧに変換する。
また、湿り度算出部１３Ｎは、変換によって得られた配管Ａの流路Ａ２内の気相容積流束ＪＧの値を、流動様式線図（図５参照）の横軸に当てはめる。
更に、湿り度算出部１３Ｎは、流動様式線図の横軸の気相容積流束ＪＧの値に対応するプロット「×」の縦軸の値を、配管Ａの流路Ａ２内の液相容積流束ＪＬの値として算出する。また、湿り度算出部１３Ｎは、配管Ａの流路Ａ２内の液相容積流束ＪＬの値を、配管Ａの流路Ａ２内を流れる流体の湿り度［％］に変換する。
詳細には、湿り度算出部１３Ｎは、流量算出部１３Ｌによって算出された配管Ａ内を流れる流体の流量を、配管Ａの流路Ａ２内の気相容積流束ＪＧに変換する。
また、湿り度算出部１３Ｎは、変換によって得られた配管Ａの流路Ａ２内の気相容積流束ＪＧの値を、流動様式線図（図５参照）の横軸に当てはめる。
更に、湿り度算出部１３Ｎは、流動様式線図の横軸の気相容積流束ＪＧの値に対応するプロット「×」の縦軸の値を、配管Ａの流路Ａ２内の液相容積流束ＪＬの値として算出する。また、湿り度算出部１３Ｎは、配管Ａの流路Ａ２内の液相容積流束ＪＬの値を、配管Ａの流路Ａ２内を流れる流体の湿り度［％］に変換する。
【0087】
図９は第１実施形態の流動様式判別装置１において実行される処理の一例を説明するためのフローチャートである。
図９（Ａ）に示す例では、ステップＳ１１Ａにおいて、第１超音波トランスデューサ１１の送信部１１Ａが、超音波を送信する。
次いで、ステップＳ１１Ｂでは、第２超音波トランスデューサ１２の受信部１２Ｂが、ステップＳ１１Ａにおいて第１超音波トランスデューサ１１の送信部１１Ａから送信された超音波などを受信する。
次いで、ステップＳ１１Ｃでは、第２超音波トランスデューサ１２の受信部１２Ｂが、ステップＳ１１Ｂにおいて受信した第１超音波トランスデューサ１１の送信部１１Ａからの超音波などを受信信号（電圧信号）に変換して出力する。
次いで、ステップＳ１１Ｄでは、演算部１３の第１受信信号取得部１３Ａが、ステップＳ１１Ｃにおいて出力された第２超音波トランスデューサ１２の受信部１２Ｂの受信信号（電圧信号）を取得する。
次いで、ステップＳ１１Ｅでは、例えば制御部１４が、第１超音波トランスデューサ１１の送信部１１Ａによる超音波の送信回数が予め設定された回数に到達したか否かを判定する。第１超音波トランスデューサ１１の送信部１１Ａによる超音波の送信回数が予め設定された回数に到達した場合には、図９（Ａ）に示すルーチンを終了する。一方、第１超音波トランスデューサ１１の送信部１１Ａによる超音波の送信回数が予め設定された回数に到達していない場合には、ステップＳ１１Ａに戻る。
図９（Ａ）に示す例では、ステップＳ１１Ａにおいて、第１超音波トランスデューサ１１の送信部１１Ａが、超音波を送信する。
次いで、ステップＳ１１Ｂでは、第２超音波トランスデューサ１２の受信部１２Ｂが、ステップＳ１１Ａにおいて第１超音波トランスデューサ１１の送信部１１Ａから送信された超音波などを受信する。
次いで、ステップＳ１１Ｃでは、第２超音波トランスデューサ１２の受信部１２Ｂが、ステップＳ１１Ｂにおいて受信した第１超音波トランスデューサ１１の送信部１１Ａからの超音波などを受信信号（電圧信号）に変換して出力する。
次いで、ステップＳ１１Ｄでは、演算部１３の第１受信信号取得部１３Ａが、ステップＳ１１Ｃにおいて出力された第２超音波トランスデューサ１２の受信部１２Ｂの受信信号（電圧信号）を取得する。
次いで、ステップＳ１１Ｅでは、例えば制御部１４が、第１超音波トランスデューサ１１の送信部１１Ａによる超音波の送信回数が予め設定された回数に到達したか否かを判定する。第１超音波トランスデューサ１１の送信部１１Ａによる超音波の送信回数が予め設定された回数に到達した場合には、図９（Ａ）に示すルーチンを終了する。一方、第１超音波トランスデューサ１１の送信部１１Ａによる超音波の送信回数が予め設定された回数に到達していない場合には、ステップＳ１１Ａに戻る。
【0088】
図９（Ｂ）に示すステップＳ１２は、図９（Ａ）に示すルーチンが終了した場合に実行される。
図９（Ｂ）に示す例では、ステップＳ１２Ａにおいて、演算部１３の第１標準偏差算出部１３Ｃが、第２超音波トランスデューサ１２の受信部１２Ｂによる第１超音波トランスデューサ１１の送信部１１Ａから送信された超音波の受信期間ＴＡ～ＴＢを含む期間Ｔ０～Ｔ３中の各時刻における上述した回数分の第２超音波トランスデューサ１２の受信部１２Ｂの受信信号の振幅の標準偏差同等値を算出する。
次いで、ステップＳ１２Ｂでは、演算部１３の流動様式判別部１３Ｅは、受信期間ＴＡ～ＴＢ中の各時刻における複数回分の第２超音波トランスデューサ１２の受信部１２Ｂの受信信号の振幅の標準偏差同等値（例えば図３（Ａ）、図４（Ａ）、図６（Ａ）、図７（Ａ）および図８（Ａ）に「時間波形Ｗ３」で示す標準偏差）に基づいて、配管Ａの流路Ａ２内を流れる流体の流動様式を判別する。
図９（Ｂ）に示す例では、ステップＳ１２Ａにおいて、演算部１３の第１標準偏差算出部１３Ｃが、第２超音波トランスデューサ１２の受信部１２Ｂによる第１超音波トランスデューサ１１の送信部１１Ａから送信された超音波の受信期間ＴＡ～ＴＢを含む期間Ｔ０～Ｔ３中の各時刻における上述した回数分の第２超音波トランスデューサ１２の受信部１２Ｂの受信信号の振幅の標準偏差同等値を算出する。
次いで、ステップＳ１２Ｂでは、演算部１３の流動様式判別部１３Ｅは、受信期間ＴＡ～ＴＢ中の各時刻における複数回分の第２超音波トランスデューサ１２の受信部１２Ｂの受信信号の振幅の標準偏差同等値（例えば図３（Ａ）、図４（Ａ）、図６（Ａ）、図７（Ａ）および図８（Ａ）に「時間波形Ｗ３」で示す標準偏差）に基づいて、配管Ａの流路Ａ２内を流れる流体の流動様式を判別する。
【0089】
図１０は第１実施形態の流動様式判別装置１において実行される他の処理の一例を説明するためのフローチャートである。
図１０（Ａ）に示す例では、ステップＳ１３Ａにおいて、第２超音波トランスデューサ１２の送信部１２Ａが、超音波を送信する。
次いで、ステップＳ１３Ｂでは、第１超音波トランスデューサ１１の受信部１１Ｂが、ステップＳ１３Ａにおいて第２超音波トランスデューサ１２の送信部１２Ａから送信された超音波などを受信する。
次いで、ステップＳ１３Ｃでは、第１超音波トランスデューサ１１の受信部１１Ｂが、ステップＳ１３Ｂにおいて受信した第２超音波トランスデューサ１２の送信部１２Ａからの超音波などを受信信号（電圧信号）に変換して出力する。
次いで、ステップＳ１３Ｄでは、演算部１３の第２受信信号取得部１３Ｂが、ステップＳ１３Ｃにおいて出力された第１超音波トランスデューサ１１の受信部１１Ｂの受信信号（電圧信号）を取得する。
次いで、ステップＳ１３Ｅでは、例えば制御部１４が、第２超音波トランスデューサ１２の送信部１２Ａによる超音波の送信回数が予め設定された回数に到達したか否かを判定する。第２超音波トランスデューサ１２の送信部１２Ａによる超音波の送信回数が予め設定された回数に到達した場合には、図１０（Ａ）に示すルーチンを終了する。一方、第２超音波トランスデューサ１２の送信部１２Ａによる超音波の送信回数が予め設定された回数に到達していない場合には、ステップＳ１３Ａに戻る。
図１０（Ａ）に示す例では、ステップＳ１３Ａにおいて、第２超音波トランスデューサ１２の送信部１２Ａが、超音波を送信する。
次いで、ステップＳ１３Ｂでは、第１超音波トランスデューサ１１の受信部１１Ｂが、ステップＳ１３Ａにおいて第２超音波トランスデューサ１２の送信部１２Ａから送信された超音波などを受信する。
次いで、ステップＳ１３Ｃでは、第１超音波トランスデューサ１１の受信部１１Ｂが、ステップＳ１３Ｂにおいて受信した第２超音波トランスデューサ１２の送信部１２Ａからの超音波などを受信信号（電圧信号）に変換して出力する。
次いで、ステップＳ１３Ｄでは、演算部１３の第２受信信号取得部１３Ｂが、ステップＳ１３Ｃにおいて出力された第１超音波トランスデューサ１１の受信部１１Ｂの受信信号（電圧信号）を取得する。
次いで、ステップＳ１３Ｅでは、例えば制御部１４が、第２超音波トランスデューサ１２の送信部１２Ａによる超音波の送信回数が予め設定された回数に到達したか否かを判定する。第２超音波トランスデューサ１２の送信部１２Ａによる超音波の送信回数が予め設定された回数に到達した場合には、図１０（Ａ）に示すルーチンを終了する。一方、第２超音波トランスデューサ１２の送信部１２Ａによる超音波の送信回数が予め設定された回数に到達していない場合には、ステップＳ１３Ａに戻る。
【0090】
図１０（Ｂ）に示すステップＳ１４は、図９（Ａ）に示すルーチンが終了し、かつ、図１０（Ａ）に示すルーチンが終了した場合に実行される。
図１０（Ｂ）に示す例では、ステップＳ１４Ａにおいて、演算部１３の第２標準偏差算出部１３Ｄが、第１超音波トランスデューサ１１の受信部１１Ｂによる第２超音波トランスデューサ１２の送信部１２Ａから送信された超音波の受信期間Ｔ５Ａ～Ｔ５Ｂを含む期間Ｔ０～Ｔ３中の各時刻における上述した回数分の第１超音波トランスデューサ１１の受信部１１Ｂの受信信号の振幅の標準偏差同等値を算出する。
図１０（Ｂ）に示す例では、ステップＳ１４Ａにおいて、演算部１３の第２標準偏差算出部１３Ｄが、第１超音波トランスデューサ１１の受信部１１Ｂによる第２超音波トランスデューサ１２の送信部１２Ａから送信された超音波の受信期間Ｔ５Ａ～Ｔ５Ｂを含む期間Ｔ０～Ｔ３中の各時刻における上述した回数分の第１超音波トランスデューサ１１の受信部１１Ｂの受信信号の振幅の標準偏差同等値を算出する。
【0091】
次いで、ステップＳ１４Ｂでは、演算部１３の第１超音波受信時刻算出部１３Ｆが、図９（Ｂ）のステップＳ１２Ａにおいて算出された第２超音波トランスデューサ１２の受信部１２Ｂの受信信号の振幅の標準偏差同等値に基づいて、第２超音波トランスデューサ１２の受信部１２Ｂによる第１超音波トランスデューサ１１の送信部１１Ａから送信された超音波の受信時刻（例えば時刻Ｔ４Ａ）を算出する。
また、ステップＳ１４Ｃでは、演算部１３の第２超音波受信時刻算出部１３Ｇが、ステップＳ１４Ａにおいて算出された第１超音波トランスデューサ１１の受信部１１Ｂの受信信号の振幅の標準偏差同等値に基づいて、第１超音波トランスデューサ１１の受信部１１Ｂによる第２超音波トランスデューサ１２の送信部１２Ａから送信された超音波の受信時刻（例えば時刻Ｔ５Ａ）を算出する。
また、ステップＳ１４Ｃでは、演算部１３の第２超音波受信時刻算出部１３Ｇが、ステップＳ１４Ａにおいて算出された第１超音波トランスデューサ１１の受信部１１Ｂの受信信号の振幅の標準偏差同等値に基づいて、第１超音波トランスデューサ１１の受信部１１Ｂによる第２超音波トランスデューサ１２の送信部１２Ａから送信された超音波の受信時刻（例えば時刻Ｔ５Ａ）を算出する。
【0092】
次いで、ステップＳ１４Ｄでは、演算部１３の第１超音波伝播時間算出部１３Ｈが、ステップＳ１４Ｂにおいて算出された第２超音波トランスデューサ１２の受信部１２Ｂによる第１超音波トランスデューサ１１の送信部１１Ａから送信された超音波の受信時刻（例えば時刻Ｔ４Ａ）と、第１超音波トランスデューサ１１の送信部１１Ａによる超音波の送信時刻Ｔ０とに基づいて、第１超音波トランスデューサ１１の送信部１１Ａから送信された超音波の伝播時間（例えば伝播時間（Ｔ４Ａ－Ｔ０））を算出する。
また、ステップＳ１４Ｅでは、演算部１３の第２超音波伝播時間算出部１３Ｊが、ステップＳ１４Ｃにおいて算出された第１超音波トランスデューサ１１の受信部１１Ｂによる第２超音波トランスデューサ１２の送信部１２Ａから送信された超音波の受信時刻（例えば時刻Ｔ５Ａ）と、第２超音波トランスデューサ１２の送信部１２Ａによる超音波の送信時刻Ｔ０とに基づいて、第２超音波トランスデューサ１２の送信部１２Ａから送信された超音波の伝播時間（例えば伝播時間（Ｔ５Ａ－Ｔ０））を算出する。
また、ステップＳ１４Ｅでは、演算部１３の第２超音波伝播時間算出部１３Ｊが、ステップＳ１４Ｃにおいて算出された第１超音波トランスデューサ１１の受信部１１Ｂによる第２超音波トランスデューサ１２の送信部１２Ａから送信された超音波の受信時刻（例えば時刻Ｔ５Ａ）と、第２超音波トランスデューサ１２の送信部１２Ａによる超音波の送信時刻Ｔ０とに基づいて、第２超音波トランスデューサ１２の送信部１２Ａから送信された超音波の伝播時間（例えば伝播時間（Ｔ５Ａ－Ｔ０））を算出する。
【0093】
次いで、ステップＳ１４Ｆでは、演算部１３の伝播時間差算出部１３Ｋが、ステップＳ１４Ｄにおいて算出された第１超音波トランスデューサ１１の送信部１１Ａから送信された超音波の伝播時間（例えば伝播時間（Ｔ４Ａ－Ｔ０））と、ステップＳ１４Ｅにおいて算出された第２超音波トランスデューサ１２の送信部１２Ａから送信された超音波の伝播時間（例えば伝播時間（Ｔ５Ａ－Ｔ０））との差分である伝播時間差（例えば伝播時間差（Ｔ５Ａ－Ｔ４Ａ））を算出する。
他の例（第１超音波トランスデューサ１１の送信部１１Ａから送信された超音波の受信期間Ｔ４Ａ～Ｔ４Ｂ中における第２超音波トランスデューサ１２の受信部１２Ｂの受信信号が、ノイズに埋もれていない例）では、ステップＳ１４Ｆにおいて、演算部１３の伝播時間差算出部１３Ｋが、第２超音波トランスデューサ１２の受信部１２Ｂによって出力される受信信号の時間波形と、第１超音波トランスデューサ１１の受信部１１Ｂによって出力される受信信号の時間波形とに基づいて、伝播時間差（例えば伝播時間差（Ｔ５Ａ－Ｔ４Ａ））を算出してもよい。
図１０に示す例では、次いで、ステップＳ１４Ｇでは、演算部１３の流量算出部１３Ｌが、ステップＳ１４Ｆにおいて算出された伝播時間差（例えば伝播時間差（Ｔ５Ａ－Ｔ４Ａ））に基づいて、配管Ａの流路Ａ２内を流れる流体の流量を算出する。
他の例（第１超音波トランスデューサ１１の送信部１１Ａから送信された超音波の受信期間Ｔ４Ａ～Ｔ４Ｂ中における第２超音波トランスデューサ１２の受信部１２Ｂの受信信号が、ノイズに埋もれていない例）では、ステップＳ１４Ｆにおいて、演算部１３の伝播時間差算出部１３Ｋが、第２超音波トランスデューサ１２の受信部１２Ｂによって出力される受信信号の時間波形と、第１超音波トランスデューサ１１の受信部１１Ｂによって出力される受信信号の時間波形とに基づいて、伝播時間差（例えば伝播時間差（Ｔ５Ａ－Ｔ４Ａ））を算出してもよい。
図１０に示す例では、次いで、ステップＳ１４Ｇでは、演算部１３の流量算出部１３Ｌが、ステップＳ１４Ｆにおいて算出された伝播時間差（例えば伝播時間差（Ｔ５Ａ－Ｔ４Ａ））に基づいて、配管Ａの流路Ａ２内を流れる流体の流量を算出する。
【0094】
次いで、ステップＳ１４Ｈでは、演算部１３の流動様式線図取得部１３Ｍが、流動様式線図（例えば図５参照）を取得する。
次いで、ステップＳ１４Ｊでは、演算部１３の湿り度算出部１３Ｎが、ステップＳ１４Ｈにおいて取得された流動様式線図と、図９（Ｂ）のステップＳ１２Ｂにおいて判別された配管Ａの流路Ａ２内を流れる流体の流動様式と、ステップＳ１４Ｇにおいて算出された配管Ａの流路Ａ２内を流れる流体の流量とに基づいて、配管Ａの流路Ａ２内を流れる流体の湿り度［％］を算出する。
次いで、ステップＳ１４Ｊでは、演算部１３の湿り度算出部１３Ｎが、ステップＳ１４Ｈにおいて取得された流動様式線図と、図９（Ｂ）のステップＳ１２Ｂにおいて判別された配管Ａの流路Ａ２内を流れる流体の流動様式と、ステップＳ１４Ｇにおいて算出された配管Ａの流路Ａ２内を流れる流体の流量とに基づいて、配管Ａの流路Ａ２内を流れる流体の湿り度［％］を算出する。
【0095】
上述したように第１実施形態の流動様式判別装置１では、演算部１３の流動様式判別部１３Ｅが、受信期間ＴＡ～ＴＢ中の各時刻における複数回分の第２超音波トランスデューサ１２の受信部１２Ｂの受信信号の振幅の標準偏差同等値に基づいて、配管Ａの流路Ａ２内を流れる流体の流動様式を判別する。
そのため、第１実施形態の流動様式判別装置１では、配管Ａの流路Ａ２内を流れる流体の流動様式を、抜管する必要なく判別することができる。
そのため、第１実施形態の流動様式判別装置１では、配管Ａの流路Ａ２内を流れる流体の流動様式を、抜管する必要なく判別することができる。
【0096】
詳細には、第１実施形態の流動様式判別装置１では、演算部１３の流動様式判別部１３Ｅが、受信期間ＴＡ～ＴＢ中の各時刻における複数回分の第２超音波トランスデューサ１２の受信部１２Ｂの受信信号の振幅の標準偏差同等値の時間波形Ｗ３（図３（Ａ）、図４（Ａ）、図６（Ａ）、図７（Ａ）および図８（Ａ）参照）の形状に基づいて、配管Ａの流路Ａ２内を流れる流体の流動様式が、乾き蒸気、層状流、波状流、波状環状遷移および環状噴霧流のいずれであるかを判別する。
そのため、第１実施形態の流動様式判別装置１では、抜管する必要なく、配管Ａの流路Ａ２内を流れる流体の流動様式を詳細に判別することができる。
そのため、第１実施形態の流動様式判別装置１では、抜管する必要なく、配管Ａの流路Ａ２内を流れる流体の流動様式を詳細に判別することができる。
【0097】
また、第１実施形態の流動様式判別装置１では、演算部１３の湿り度算出部１３Ｎが、流動様式線図（図５参照）と、配管Ａの流路Ａ２内を流れる流体の流動様式と、配管Ａの流路Ａ２内を流れる流体の流量とに基づいて、配管Ａの流路Ａ２内を流れる流体の湿り度を算出する。
そのため、第１実施形態の流動様式判別装置１では、抜管する必要なく、配管Ａの流路Ａ２内を流れる流体の湿り度を得ることができる。
そのため、第１実施形態の流動様式判別装置１では、抜管する必要なく、配管Ａの流路Ａ２内を流れる流体の湿り度を得ることができる。
【0098】
＜適用例＞
エネルギーの有効利用の観点から、工場などの既設配管において蒸気流量を計測することが重要とされる。その手法の一つとして、第１実施形態の流動様式判別装置１を適用したクランプオン式超音波流量計が考えられる。
クランプオン式超音波流量計を用いて蒸気流量を計測する場合、蒸気と配管の音響インピーダンスの違いによる透過強度の低下、高温によるウェッジ材・センサ感度の低下、蒸気の湿りによる液滴・液膜の影響により、超音波透過強度が低下する。
第１実施形態の流動様式判別装置１を適用したクランプオン式超音波流量計では、低信号ノイズ比条件下においても、配管Ａ内を流れる流体の流量を適切に算出することができる。
第１実施形態の流動様式判別装置１を適用することにより、蒸気流量計測の実現によるエネルギーマネジメントへの応用の他、従来、伝播時間差式流量計では困難であった条件（低圧気体、天然ガスなどの混相流）において、計測を可能とする可能性がある。
蒸気を主力熱源とする工場などでは、省エネルギー、エネルギー診断から最適なエネルギーシステムの提案、設計などを行うＥＳＰ（エネルギーサービスプロバイダ）の観点から、蒸気流量計測を行うことが望まれる。第１実施形態の流動様式判別装置１を適用して蒸気流量計測を行うことによって、蒸気管を抜管する必要なく、蒸気の流量を把握することができる。第１実施形態の流動様式判別装置１を適用して蒸気流量計測を行うことにより、超音波トランスデューサによって受信される雑音が大きい場合であっても、適切に蒸気の流量を計測することができ、省エネルギーに資することができる。また、第１実施形態の流動様式判別装置１によって配管Ａの流路Ａ２内を流れる流体の流動様式、湿り度を把握することにより、配管Ａの熱ロスを解析し、適切な省エネルギー対策を施すことができる。
エネルギーの有効利用の観点から、工場などの既設配管において蒸気流量を計測することが重要とされる。その手法の一つとして、第１実施形態の流動様式判別装置１を適用したクランプオン式超音波流量計が考えられる。
クランプオン式超音波流量計を用いて蒸気流量を計測する場合、蒸気と配管の音響インピーダンスの違いによる透過強度の低下、高温によるウェッジ材・センサ感度の低下、蒸気の湿りによる液滴・液膜の影響により、超音波透過強度が低下する。
第１実施形態の流動様式判別装置１を適用したクランプオン式超音波流量計では、低信号ノイズ比条件下においても、配管Ａ内を流れる流体の流量を適切に算出することができる。
第１実施形態の流動様式判別装置１を適用することにより、蒸気流量計測の実現によるエネルギーマネジメントへの応用の他、従来、伝播時間差式流量計では困難であった条件（低圧気体、天然ガスなどの混相流）において、計測を可能とする可能性がある。
蒸気を主力熱源とする工場などでは、省エネルギー、エネルギー診断から最適なエネルギーシステムの提案、設計などを行うＥＳＰ（エネルギーサービスプロバイダ）の観点から、蒸気流量計測を行うことが望まれる。第１実施形態の流動様式判別装置１を適用して蒸気流量計測を行うことによって、蒸気管を抜管する必要なく、蒸気の流量を把握することができる。第１実施形態の流動様式判別装置１を適用して蒸気流量計測を行うことにより、超音波トランスデューサによって受信される雑音が大きい場合であっても、適切に蒸気の流量を計測することができ、省エネルギーに資することができる。また、第１実施形態の流動様式判別装置１によって配管Ａの流路Ａ２内を流れる流体の流動様式、湿り度を把握することにより、配管Ａの熱ロスを解析し、適切な省エネルギー対策を施すことができる。
【0099】
［第２実施形態］
以下、本発明の流動様式判別装置、流動様式判別システムおよび流動様式判別方法の第２実施形態について説明する。
第２実施形態の流動様式判別システムＳは、後述する点を除き、上述した第１実施形態の流動様式判別装置１と同様に構成されている。従って、第２実施形態の流動様式判別システムＳによれば、後述する点を除き、上述した第１実施形態の流動様式判別装置１と同様の効果を奏することができる。
以下、本発明の流動様式判別装置、流動様式判別システムおよび流動様式判別方法の第２実施形態について説明する。
第２実施形態の流動様式判別システムＳは、後述する点を除き、上述した第１実施形態の流動様式判別装置１と同様に構成されている。従って、第２実施形態の流動様式判別システムＳによれば、後述する点を除き、上述した第１実施形態の流動様式判別装置１と同様の効果を奏することができる。
【0100】
図１１は第２実施形態の流動様式判別システムＳの一例を示す図である。
図１１に示す例では、流動様式判別システムＳが、配管Ａ（図１（Ｂ）参照）の流路Ａ２内を流れる流体の流動様式を判別する。流動様式判別システムＳは、超音波流量計Ｓ１と、演算部Ｓ２とを備えている。演算部Ｓ２は、例えば配管Ａから離れた位置に配置される。
図１に示す例では、流動様式判別装置１が、第１超音波トランスデューサ１１と、第２超音波トランスデューサ１２と、演算部１３と、制御部１４とを備えているが、図１１に示す例では、超音波流量計Ｓ１が、第１超音波トランスデューサ１１と、第２超音波トランスデューサ１２と、制御部１４と、記憶部１５とを備えている。
記憶部１５は、第１超音波トランスデューサ１１の送信部１１Ａによる超音波の送信時刻、第２超音波トランスデューサ１２の送信部１２Ａによる超音波の送信時刻、第２超音波トランスデューサ１２の受信部１２Ｂによって出力される受信信号、第１超音波トランスデューサ１１の受信部１１Ｂによって出力される受信信号などを記憶する。
図１１に示す例では、流動様式判別システムＳが、配管Ａ（図１（Ｂ）参照）の流路Ａ２内を流れる流体の流動様式を判別する。流動様式判別システムＳは、超音波流量計Ｓ１と、演算部Ｓ２とを備えている。演算部Ｓ２は、例えば配管Ａから離れた位置に配置される。
図１に示す例では、流動様式判別装置１が、第１超音波トランスデューサ１１と、第２超音波トランスデューサ１２と、演算部１３と、制御部１４とを備えているが、図１１に示す例では、超音波流量計Ｓ１が、第１超音波トランスデューサ１１と、第２超音波トランスデューサ１２と、制御部１４と、記憶部１５とを備えている。
記憶部１５は、第１超音波トランスデューサ１１の送信部１１Ａによる超音波の送信時刻、第２超音波トランスデューサ１２の送信部１２Ａによる超音波の送信時刻、第２超音波トランスデューサ１２の受信部１２Ｂによって出力される受信信号、第１超音波トランスデューサ１１の受信部１１Ｂによって出力される受信信号などを記憶する。
【0101】
図１１に示す例では、演算部Ｓ２が、第１超音波トランスデューサ１１の送信部１１Ａから送信された超音波などを受信した第２超音波トランスデューサ１２の受信部１２Ｂによって出力される受信信号に基づいて、配管Ａの流路Ａ２内を流れる流体の流動様式を判別する等の演算を行う。
詳細には、図１１に示す例では、配管Ａの流路Ａ２内を流れる流体の流動様式を判別するために、第１超音波トランスデューサ１１の送信部１１Ａが超音波を複数回（例えば２０００回以上）送信し、第２超音波トランスデューサ１２の受信部１２Ｂは、第１超音波トランスデューサ１１の送信部１１Ａから複数回送信された超音波を受信して、複数回分の受信信号（電圧信号）を出力する。
第２超音波トランスデューサ１２の受信部１２Ｂが出力する複数回分の受信信号（電圧信号）は、配管Ａの流路Ａ２内を流れる流体の流量の算出にも用いられる。
更に、配管Ａの流路Ａ２内を流れる流体の流量を算出するために、第２超音波トランスデューサ１２の送信部１２Ａも超音波を複数回（例えば２０００回以上）送信し、第１超音波トランスデューサ１１の受信部１１Ｂが、第２超音波トランスデューサ１２の送信部１２Ａから複数回送信された超音波を受信して、複数回分の受信信号（電圧信号）を出力する。
つまり、演算部Ｓ２は、第１超音波トランスデューサ１１の送信部１１Ａから送信された超音波などを受信した第２超音波トランスデューサ１２の受信部１２Ｂによって出力される受信信号と、第２超音波トランスデューサ１２の送信部１２Ａから送信された超音波などを受信した第１超音波トランスデューサ１１の受信部１１Ｂによって出力される受信信号とに基づいて、配管Ａの流路Ａ２内を流れる流体の流量を算出する。
詳細には、配管Ａの流路Ａ２内を流れる流体の流動様式の判別および流量の算出を行うために、演算部Ｓ２は、第１受信信号取得部２Ａと、第２受信信号取得部２Ｂと、第１標準偏差算出部２Ｃと、第２標準偏差算出部２Ｄと、流動様式判別部２Ｅと、第１超音波受信時刻算出部２Ｆと、第２超音波受信時刻算出部２Ｇと、第１超音波伝播時間算出部２Ｈと、第２超音波伝播時間算出部２Ｊと、伝播時間差算出部２Ｋと、流量算出部２Ｌと、流動様式線図取得部２Ｍと、湿り度算出部２Ｎとを備えている。
図１１に示す例では、演算部Ｓ２が、第１超音波受信時刻算出部２Ｆと、第２超音波受信時刻算出部２Ｇと、第１超音波伝播時間算出部２Ｈと、第２超音波伝播時間算出部２Ｊとを備えているが、他の例では、演算部Ｓ２が、第１超音波受信時刻算出部２Ｆと、第２超音波受信時刻算出部２Ｇと、第１超音波伝播時間算出部２Ｈと、第２超音波伝播時間算出部２Ｊとを備えていなくてもよい。
詳細には、図１１に示す例では、配管Ａの流路Ａ２内を流れる流体の流動様式を判別するために、第１超音波トランスデューサ１１の送信部１１Ａが超音波を複数回（例えば２０００回以上）送信し、第２超音波トランスデューサ１２の受信部１２Ｂは、第１超音波トランスデューサ１１の送信部１１Ａから複数回送信された超音波を受信して、複数回分の受信信号（電圧信号）を出力する。
第２超音波トランスデューサ１２の受信部１２Ｂが出力する複数回分の受信信号（電圧信号）は、配管Ａの流路Ａ２内を流れる流体の流量の算出にも用いられる。
更に、配管Ａの流路Ａ２内を流れる流体の流量を算出するために、第２超音波トランスデューサ１２の送信部１２Ａも超音波を複数回（例えば２０００回以上）送信し、第１超音波トランスデューサ１１の受信部１１Ｂが、第２超音波トランスデューサ１２の送信部１２Ａから複数回送信された超音波を受信して、複数回分の受信信号（電圧信号）を出力する。
つまり、演算部Ｓ２は、第１超音波トランスデューサ１１の送信部１１Ａから送信された超音波などを受信した第２超音波トランスデューサ１２の受信部１２Ｂによって出力される受信信号と、第２超音波トランスデューサ１２の送信部１２Ａから送信された超音波などを受信した第１超音波トランスデューサ１１の受信部１１Ｂによって出力される受信信号とに基づいて、配管Ａの流路Ａ２内を流れる流体の流量を算出する。
詳細には、配管Ａの流路Ａ２内を流れる流体の流動様式の判別および流量の算出を行うために、演算部Ｓ２は、第１受信信号取得部２Ａと、第２受信信号取得部２Ｂと、第１標準偏差算出部２Ｃと、第２標準偏差算出部２Ｄと、流動様式判別部２Ｅと、第１超音波受信時刻算出部２Ｆと、第２超音波受信時刻算出部２Ｇと、第１超音波伝播時間算出部２Ｈと、第２超音波伝播時間算出部２Ｊと、伝播時間差算出部２Ｋと、流量算出部２Ｌと、流動様式線図取得部２Ｍと、湿り度算出部２Ｎとを備えている。
図１１に示す例では、演算部Ｓ２が、第１超音波受信時刻算出部２Ｆと、第２超音波受信時刻算出部２Ｇと、第１超音波伝播時間算出部２Ｈと、第２超音波伝播時間算出部２Ｊとを備えているが、他の例では、演算部Ｓ２が、第１超音波受信時刻算出部２Ｆと、第２超音波受信時刻算出部２Ｇと、第１超音波伝播時間算出部２Ｈと、第２超音波伝播時間算出部２Ｊとを備えていなくてもよい。
【0102】
第１受信信号取得部２Ａは、図１に示す第１受信信号取得部１３Ａと同様の機能を有する。第２受信信号取得部２Ｂは、図１に示す第２受信信号取得部１３Ｂと同様の機能を有する。第１標準偏差算出部２Ｃは、図１に示す第１標準偏差算出部１３Ｃと同様の機能を有する。第２標準偏差算出部２Ｄは、図１に示す第２標準偏差算出部１３Ｄと同様の機能を有する。流動様式判別部２Ｅは、図１に示す流動様式判別部１３Ｅと同様の機能を有する。第１超音波受信時刻算出部２Ｆは、図１に示す第１超音波受信時刻算出部１３Ｆと同様の機能を有する。第２超音波受信時刻算出部２Ｇは、図１に示す第２超音波受信時刻算出部１３Ｇと同様の機能を有する。第１超音波伝播時間算出部２Ｈは、図１に示す第１超音波伝播時間算出部１３Ｈと同様の機能を有する。第２超音波伝播時間算出部２Ｊは、図１に示す第２超音波伝播時間算出部１３Ｊと同様の機能を有する。伝播時間差算出部２Ｋは、図１に示す伝播時間差算出部１３Ｋと同様の機能を有する。流量算出部２Ｌは、図１に示す流量算出部１３Ｌと同様の機能を有する。流動様式線図取得部２Ｍは、図１に示す流動様式線図取得部１３Ｍと同様の機能を有する。湿り度算出部２Ｎは、図１に示す湿り度算出部１３Ｎと同様の機能を有する。
他の例では、伝播時間差算出部２Ｋが、第２超音波トランスデューサ１２の受信部１２Ｂによる第１超音波トランスデューサ１１の送信部１１Ａからの超音波の受信時刻（例えば時刻Ｔ４Ａ）と、第１超音波トランスデューサ１１の受信部１１Ｂによる第２超音波トランスデューサ１２の送信部１２Ａからの超音波の受信時刻（例えば時刻Ｔ５Ａ）とを用いることなく、伝播時間差（例えば伝播時間差（Ｔ５Ａ－Ｔ４Ａ））を算出してもよい。
具体的には、伝播時間差算出部２Ｋは、第２超音波トランスデューサ１２の受信部１２Ｂによって出力される受信信号の振幅の標準偏差同等値の時間波形と、第１超音波トランスデューサ１１の受信部１１Ｂによって出力される受信信号の振幅の標準偏差同等値の時間波形との相関処理（例えば特開２０１６－１８０６７９号公報に記載されているような公知の相関処理）を実行することによって、伝播時間差（例えば伝播時間差（Ｔ５Ａ－Ｔ４Ａ））を算出することができる。
更に他の例（第１超音波トランスデューサ１１の送信部１１Ａから送信された超音波の受信期間Ｔ４Ａ～Ｔ４Ｂ中における第２超音波トランスデューサ１２の受信部１２Ｂの受信信号が、ノイズに埋もれておらず、第２超音波トランスデューサ１２の送信部１２Ａから送信された超音波の受信期間Ｔ５Ａ～Ｔ５Ｂ中における第１超音波トランスデューサ１１の受信部１１Ｂの受信信号が、ノイズに埋もれていない例）では、伝播時間差算出部２Ｋが、第２超音波トランスデューサ１２の受信部１２Ｂによって出力される受信信号の時間波形と、第１超音波トランスデューサ１１の受信部１１Ｂによって出力される受信信号の時間波形とに基づいて、伝播時間差（例えば伝播時間差（Ｔ５Ａ－Ｔ４Ａ））を算出してもよい。
他の例では、伝播時間差算出部２Ｋが、第２超音波トランスデューサ１２の受信部１２Ｂによる第１超音波トランスデューサ１１の送信部１１Ａからの超音波の受信時刻（例えば時刻Ｔ４Ａ）と、第１超音波トランスデューサ１１の受信部１１Ｂによる第２超音波トランスデューサ１２の送信部１２Ａからの超音波の受信時刻（例えば時刻Ｔ５Ａ）とを用いることなく、伝播時間差（例えば伝播時間差（Ｔ５Ａ－Ｔ４Ａ））を算出してもよい。
具体的には、伝播時間差算出部２Ｋは、第２超音波トランスデューサ１２の受信部１２Ｂによって出力される受信信号の振幅の標準偏差同等値の時間波形と、第１超音波トランスデューサ１１の受信部１１Ｂによって出力される受信信号の振幅の標準偏差同等値の時間波形との相関処理（例えば特開２０１６－１８０６７９号公報に記載されているような公知の相関処理）を実行することによって、伝播時間差（例えば伝播時間差（Ｔ５Ａ－Ｔ４Ａ））を算出することができる。
更に他の例（第１超音波トランスデューサ１１の送信部１１Ａから送信された超音波の受信期間Ｔ４Ａ～Ｔ４Ｂ中における第２超音波トランスデューサ１２の受信部１２Ｂの受信信号が、ノイズに埋もれておらず、第２超音波トランスデューサ１２の送信部１２Ａから送信された超音波の受信期間Ｔ５Ａ～Ｔ５Ｂ中における第１超音波トランスデューサ１１の受信部１１Ｂの受信信号が、ノイズに埋もれていない例）では、伝播時間差算出部２Ｋが、第２超音波トランスデューサ１２の受信部１２Ｂによって出力される受信信号の時間波形と、第１超音波トランスデューサ１１の受信部１１Ｂによって出力される受信信号の時間波形とに基づいて、伝播時間差（例えば伝播時間差（Ｔ５Ａ－Ｔ４Ａ））を算出してもよい。
【0103】
図１２は第２実施形態の流動様式判別システムＳにおいて実行される処理の一例を説明するためのフローチャートである。
図１２（Ａ）に示す例では、ステップＳ２１Ａにおいて、第１超音波トランスデューサ１１の送信部１１Ａが、超音波を送信する。
次いで、ステップＳ２１Ｂでは、第２超音波トランスデューサ１２の受信部１２Ｂが、ステップＳ２１Ａにおいて第１超音波トランスデューサ１１の送信部１１Ａから送信された超音波などを受信する。
次いで、ステップＳ２１Ｃでは、第２超音波トランスデューサ１２の受信部１２Ｂが、ステップＳ２１Ｂにおいて受信した第１超音波トランスデューサ１１の送信部１１Ａからの超音波などを受信信号（電圧信号）に変換して出力する。
次いで、ステップＳ２１Ｄでは、記憶部１５が、ステップＳ２１Ｃにおいて出力された第２超音波トランスデューサ１２の受信部１２Ｂの受信信号（電圧信号）を記憶する。
次いで、ステップＳ２１Ｅでは、例えば制御部１４が、第１超音波トランスデューサ１１の送信部１１Ａによる超音波の送信回数が予め設定された回数に到達したか否かを判定する。第１超音波トランスデューサ１１の送信部１１Ａによる超音波の送信回数が予め設定された回数に到達した場合には、図１２（Ａ）に示すルーチンを終了する。一方、第１超音波トランスデューサ１１の送信部１１Ａによる超音波の送信回数が予め設定された回数に到達していない場合には、ステップＳ２１Ａに戻る。
図１２（Ａ）に示す例では、ステップＳ２１Ａにおいて、第１超音波トランスデューサ１１の送信部１１Ａが、超音波を送信する。
次いで、ステップＳ２１Ｂでは、第２超音波トランスデューサ１２の受信部１２Ｂが、ステップＳ２１Ａにおいて第１超音波トランスデューサ１１の送信部１１Ａから送信された超音波などを受信する。
次いで、ステップＳ２１Ｃでは、第２超音波トランスデューサ１２の受信部１２Ｂが、ステップＳ２１Ｂにおいて受信した第１超音波トランスデューサ１１の送信部１１Ａからの超音波などを受信信号（電圧信号）に変換して出力する。
次いで、ステップＳ２１Ｄでは、記憶部１５が、ステップＳ２１Ｃにおいて出力された第２超音波トランスデューサ１２の受信部１２Ｂの受信信号（電圧信号）を記憶する。
次いで、ステップＳ２１Ｅでは、例えば制御部１４が、第１超音波トランスデューサ１１の送信部１１Ａによる超音波の送信回数が予め設定された回数に到達したか否かを判定する。第１超音波トランスデューサ１１の送信部１１Ａによる超音波の送信回数が予め設定された回数に到達した場合には、図１２（Ａ）に示すルーチンを終了する。一方、第１超音波トランスデューサ１１の送信部１１Ａによる超音波の送信回数が予め設定された回数に到達していない場合には、ステップＳ２１Ａに戻る。
【0104】
図１２（Ｂ）に示すステップＳ２２は、図１２（Ａ）に示すルーチンが終了した場合に実行される。
図１２（Ｂ）に示す例では、ステップＳ２２Ｘにおいて、演算部Ｓ２の第１受信信号取得部２Ａが、ステップＳ２１Ｄにおいて記憶された第２超音波トランスデューサ１２の受信部１２Ｂの受信信号（電圧信号）を取得する。
次いで、ステップＳ２２Ａでは、演算部Ｓ２の第１標準偏差算出部２Ｃが、第２超音波トランスデューサ１２の受信部１２Ｂによる第１超音波トランスデューサ１１の送信部１１Ａから送信された超音波の受信期間ＴＡ～ＴＢを含む期間Ｔ０～Ｔ３中の各時刻における上述した回数分の第２超音波トランスデューサ１２の受信部１２Ｂの受信信号の振幅の標準偏差同等値を算出する。
次いで、ステップＳ２２Ｂでは、演算部Ｓ２の流動様式判別部２Ｅは、受信期間ＴＡ～ＴＢ中の各時刻における複数回分の第２超音波トランスデューサ１２の受信部１２Ｂの受信信号の振幅の標準偏差同等値に基づいて、配管Ａの流路Ａ２内を流れる流体の流動様式を判別する。
図１２（Ｂ）に示す例では、ステップＳ２２Ｘにおいて、演算部Ｓ２の第１受信信号取得部２Ａが、ステップＳ２１Ｄにおいて記憶された第２超音波トランスデューサ１２の受信部１２Ｂの受信信号（電圧信号）を取得する。
次いで、ステップＳ２２Ａでは、演算部Ｓ２の第１標準偏差算出部２Ｃが、第２超音波トランスデューサ１２の受信部１２Ｂによる第１超音波トランスデューサ１１の送信部１１Ａから送信された超音波の受信期間ＴＡ～ＴＢを含む期間Ｔ０～Ｔ３中の各時刻における上述した回数分の第２超音波トランスデューサ１２の受信部１２Ｂの受信信号の振幅の標準偏差同等値を算出する。
次いで、ステップＳ２２Ｂでは、演算部Ｓ２の流動様式判別部２Ｅは、受信期間ＴＡ～ＴＢ中の各時刻における複数回分の第２超音波トランスデューサ１２の受信部１２Ｂの受信信号の振幅の標準偏差同等値に基づいて、配管Ａの流路Ａ２内を流れる流体の流動様式を判別する。
【0105】
図１３は第２実施形態の流動様式判別システムＳにおいて実行される他の処理の一例を説明するためのフローチャートである。
図１３（Ａ）に示す例では、ステップＳ２３Ａにおいて、第２超音波トランスデューサ１２の送信部１２Ａが、超音波を送信する。
次いで、ステップＳ２３Ｂでは、第１超音波トランスデューサ１１の受信部１１Ｂが、ステップＳ２３Ａにおいて第２超音波トランスデューサ１２の送信部１２Ａから送信された超音波などを受信する。
次いで、ステップＳ２３Ｃでは、第１超音波トランスデューサ１１の受信部１１Ｂが、ステップＳ２３Ｂにおいて受信した第２超音波トランスデューサ１２の送信部１２Ａからの超音波などを受信信号（電圧信号）に変換して出力する。
次いで、ステップＳ２３Ｄでは、記憶部１５が、ステップＳ２３Ｃにおいて出力された第１超音波トランスデューサ１１の受信部１１Ｂの受信信号（電圧信号）を記憶する。
次いで、ステップＳ２３Ｅでは、例えば制御部１４が、第２超音波トランスデューサ１２の送信部１２Ａによる超音波の送信回数が予め設定された回数に到達したか否かを判定する。第２超音波トランスデューサ１２の送信部１２Ａによる超音波の送信回数が予め設定された回数に到達した場合には、図１３（Ａ）に示すルーチンを終了する。一方、第２超音波トランスデューサ１２の送信部１２Ａによる超音波の送信回数が予め設定された回数に到達していない場合には、ステップＳ２３Ａに戻る。
図１３（Ａ）に示す例では、ステップＳ２３Ａにおいて、第２超音波トランスデューサ１２の送信部１２Ａが、超音波を送信する。
次いで、ステップＳ２３Ｂでは、第１超音波トランスデューサ１１の受信部１１Ｂが、ステップＳ２３Ａにおいて第２超音波トランスデューサ１２の送信部１２Ａから送信された超音波などを受信する。
次いで、ステップＳ２３Ｃでは、第１超音波トランスデューサ１１の受信部１１Ｂが、ステップＳ２３Ｂにおいて受信した第２超音波トランスデューサ１２の送信部１２Ａからの超音波などを受信信号（電圧信号）に変換して出力する。
次いで、ステップＳ２３Ｄでは、記憶部１５が、ステップＳ２３Ｃにおいて出力された第１超音波トランスデューサ１１の受信部１１Ｂの受信信号（電圧信号）を記憶する。
次いで、ステップＳ２３Ｅでは、例えば制御部１４が、第２超音波トランスデューサ１２の送信部１２Ａによる超音波の送信回数が予め設定された回数に到達したか否かを判定する。第２超音波トランスデューサ１２の送信部１２Ａによる超音波の送信回数が予め設定された回数に到達した場合には、図１３（Ａ）に示すルーチンを終了する。一方、第２超音波トランスデューサ１２の送信部１２Ａによる超音波の送信回数が予め設定された回数に到達していない場合には、ステップＳ２３Ａに戻る。
【0106】
図１３（Ｂ）に示すステップＳ２４は、図１２（Ａ）に示すルーチンが終了し、かつ、図１３（Ａ）に示すルーチンが終了した場合に実行される。
図１３（Ｂ）に示す例では、ステップＳ２４Ｘにおいて、演算部Ｓ２の第２受信信号取得部２Ｂが、ステップＳ２３Ｄにおいて記憶された第１超音波トランスデューサ１１の受信部１１Ｂの受信信号（電圧信号）を取得する。
次いで、ステップＳ２４Ａでは、演算部Ｓ２の第２標準偏差算出部２Ｄが、第１超音波トランスデューサ１１の受信部１１Ｂによる第２超音波トランスデューサ１２の送信部１２Ａから送信された超音波の受信期間Ｔ５Ａ～Ｔ５Ｂを含む期間Ｔ０～Ｔ３中の各時刻における上述した回数分の第１超音波トランスデューサ１１の受信部１１Ｂの受信信号の振幅の標準偏差同等値を算出する。
図１３（Ｂ）に示す例では、ステップＳ２４Ｘにおいて、演算部Ｓ２の第２受信信号取得部２Ｂが、ステップＳ２３Ｄにおいて記憶された第１超音波トランスデューサ１１の受信部１１Ｂの受信信号（電圧信号）を取得する。
次いで、ステップＳ２４Ａでは、演算部Ｓ２の第２標準偏差算出部２Ｄが、第１超音波トランスデューサ１１の受信部１１Ｂによる第２超音波トランスデューサ１２の送信部１２Ａから送信された超音波の受信期間Ｔ５Ａ～Ｔ５Ｂを含む期間Ｔ０～Ｔ３中の各時刻における上述した回数分の第１超音波トランスデューサ１１の受信部１１Ｂの受信信号の振幅の標準偏差同等値を算出する。
【0107】
次いで、ステップＳ２４Ｂでは、演算部Ｓ２の第１超音波受信時刻算出部２Ｆが、図１２（Ｂ）のステップＳ２２Ａにおいて算出された第２超音波トランスデューサ１２の受信部１２Ｂの受信信号の振幅の標準偏差同等値に基づいて、第２超音波トランスデューサ１２の受信部１２Ｂによる第１超音波トランスデューサ１１の送信部１１Ａから送信された超音波の受信時刻（例えば時刻Ｔ４Ａ）を算出する。
また、ステップＳ２４Ｃでは、演算部Ｓ２の第２超音波受信時刻算出部２Ｇが、ステップＳ２４Ａにおいて算出された第１超音波トランスデューサ１１の受信部１１Ｂの受信信号の振幅の標準偏差同等値に基づいて、第１超音波トランスデューサ１１の受信部１１Ｂによる第２超音波トランスデューサ１２の送信部１２Ａから送信された超音波の受信時刻（例えば時刻Ｔ５Ａ）を算出する。
また、ステップＳ２４Ｃでは、演算部Ｓ２の第２超音波受信時刻算出部２Ｇが、ステップＳ２４Ａにおいて算出された第１超音波トランスデューサ１１の受信部１１Ｂの受信信号の振幅の標準偏差同等値に基づいて、第１超音波トランスデューサ１１の受信部１１Ｂによる第２超音波トランスデューサ１２の送信部１２Ａから送信された超音波の受信時刻（例えば時刻Ｔ５Ａ）を算出する。
【0108】
次いで、ステップＳ２４Ｄでは、演算部Ｓ２の第１超音波伝播時間算出部２Ｈが、ステップＳ２４Ｂにおいて算出された第２超音波トランスデューサ１２の受信部１２Ｂによる第１超音波トランスデューサ１１の送信部１１Ａから送信された超音波の受信時刻（例えば時刻Ｔ４Ａ）と、第１超音波トランスデューサ１１の送信部１１Ａによる超音波の送信時刻Ｔ０とに基づいて、第１超音波トランスデューサ１１の送信部１１Ａから送信された超音波の伝播時間（例えば伝播時間（Ｔ４Ａ－Ｔ０））を算出する。
また、ステップＳ２４Ｅでは、演算部Ｓ２の第２超音波伝播時間算出部２Ｊが、ステップＳ２４Ｃにおいて算出された第１超音波トランスデューサ１１の受信部１１Ｂによる第２超音波トランスデューサ１２の送信部１２Ａから送信された超音波の受信時刻（例えば時刻Ｔ５Ａ）と、第２超音波トランスデューサ１２の送信部１２Ａによる超音波の送信時刻Ｔ０とに基づいて、第２超音波トランスデューサ１２の送信部１２Ａから送信された超音波の伝播時間（例えば伝播時間（Ｔ５Ａ－Ｔ０））を算出する。
また、ステップＳ２４Ｅでは、演算部Ｓ２の第２超音波伝播時間算出部２Ｊが、ステップＳ２４Ｃにおいて算出された第１超音波トランスデューサ１１の受信部１１Ｂによる第２超音波トランスデューサ１２の送信部１２Ａから送信された超音波の受信時刻（例えば時刻Ｔ５Ａ）と、第２超音波トランスデューサ１２の送信部１２Ａによる超音波の送信時刻Ｔ０とに基づいて、第２超音波トランスデューサ１２の送信部１２Ａから送信された超音波の伝播時間（例えば伝播時間（Ｔ５Ａ－Ｔ０））を算出する。
【0109】
次いで、ステップＳ２４Ｆでは、演算部Ｓ２の伝播時間差算出部２Ｋが、ステップＳ２４Ｄにおいて算出された第１超音波トランスデューサ１１の送信部１１Ａから送信された超音波の伝播時間（例えば伝播時間（Ｔ４Ａ－Ｔ０））と、ステップＳ２４Ｅにおいて算出された第２超音波トランスデューサ１２の送信部１２Ａから送信された超音波の伝播時間（例えば伝播時間（Ｔ５Ａ－Ｔ０））との差分である伝播時間差（例えば伝播時間差（Ｔ５Ａ－Ｔ４Ａ））を算出する。
他の例（第１超音波トランスデューサ１１の送信部１１Ａから送信された超音波の受信期間Ｔ４Ａ～Ｔ４Ｂ中における第２超音波トランスデューサ１２の受信部１２Ｂの受信信号が、ノイズに埋もれていない例）では、ステップＳ２４Ｆにおいて、演算部Ｓ２の伝播時間差算出部２Ｋが、第２超音波トランスデューサ１２の受信部１２Ｂによって出力される受信信号の時間波形と、第１超音波トランスデューサ１１の受信部１１Ｂによって出力される受信信号の時間波形とに基づいて、伝播時間差（例えば伝播時間差（Ｔ５Ａ－Ｔ４Ａ））を算出してもよい。
図１３に示す例では、次いで、ステップＳ２４Ｇでは、演算部Ｓ２の流量算出部２Ｌが、ステップＳ２４Ｆにおいて算出された伝播時間差（例えば伝播時間差（Ｔ５Ａ－Ｔ４Ａ））に基づいて、配管Ａの流路Ａ２内を流れる流体の流量を算出する。
他の例（第１超音波トランスデューサ１１の送信部１１Ａから送信された超音波の受信期間Ｔ４Ａ～Ｔ４Ｂ中における第２超音波トランスデューサ１２の受信部１２Ｂの受信信号が、ノイズに埋もれていない例）では、ステップＳ２４Ｆにおいて、演算部Ｓ２の伝播時間差算出部２Ｋが、第２超音波トランスデューサ１２の受信部１２Ｂによって出力される受信信号の時間波形と、第１超音波トランスデューサ１１の受信部１１Ｂによって出力される受信信号の時間波形とに基づいて、伝播時間差（例えば伝播時間差（Ｔ５Ａ－Ｔ４Ａ））を算出してもよい。
図１３に示す例では、次いで、ステップＳ２４Ｇでは、演算部Ｓ２の流量算出部２Ｌが、ステップＳ２４Ｆにおいて算出された伝播時間差（例えば伝播時間差（Ｔ５Ａ－Ｔ４Ａ））に基づいて、配管Ａの流路Ａ２内を流れる流体の流量を算出する。
【0110】
次いで、ステップＳ２４Ｈでは、演算部Ｓ２の流動様式線図取得部２Ｍが、流動様式線図（例えば図５参照）を取得する。
次いで、ステップＳ２４Ｊでは、演算部Ｓ２の湿り度算出部２Ｎが、ステップＳ２４Ｈにおいて取得された流動様式線図と、図１２（Ｂ）のステップＳ２２Ｂにおいて判別された配管Ａの流路Ａ２内を流れる流体の流動様式と、ステップＳ２４Ｇにおいて算出された配管Ａの流路Ａ２内を流れる流体の流量とに基づいて、配管Ａの流路Ａ２内を流れる流体の湿り度［％］を算出する。
次いで、ステップＳ２４Ｊでは、演算部Ｓ２の湿り度算出部２Ｎが、ステップＳ２４Ｈにおいて取得された流動様式線図と、図１２（Ｂ）のステップＳ２２Ｂにおいて判別された配管Ａの流路Ａ２内を流れる流体の流動様式と、ステップＳ２４Ｇにおいて算出された配管Ａの流路Ａ２内を流れる流体の流量とに基づいて、配管Ａの流路Ａ２内を流れる流体の湿り度［％］を算出する。
【0111】
以上、本発明の実施形態を、図面を参照して詳述してきたが、具体的な構成はこれらの実施形態および各例に限られるものではなく、本発明の趣旨を逸脱しない範囲で適宜変更を加えることができる。上述した各実施形態および各例に記載の構成を組み合わせてもよい。
【符号の説明】
【0112】
１…流動様式判別装置、１１…第１超音波トランスデューサ、１１Ａ…送信部、１１Ｂ…受信部、１２…第２超音波トランスデューサ、１２Ａ…送信部、１２Ｂ…受信部、１３…演算部、１３Ａ…第１受信信号取得部、１３Ｂ…第２受信信号取得部、１３Ｃ…第１標準偏差算出部、１３Ｄ…第２標準偏差算出部、１３Ｅ…流動様式判別部、１３Ｆ…第１超音波受信時刻算出部、１３Ｇ…第２超音波受信時刻算出部、１３Ｈ…第１超音波伝播時間算出部、１３Ｊ…第２超音波伝播時間算出部、１３Ｋ…伝播時間差算出部、１３Ｌ…流量算出部、１３Ｍ…流動様式線図取得部、１３Ｎ…湿り度算出部、１４…制御部、１５…記憶部、Ｓ…流動様式判別システム、Ｓ１…超音波流量計、Ｓ２…演算部、２Ａ…第１受信信号取得部、２Ｂ…第２受信信号取得部、２Ｃ…第１標準偏差算出部、２Ｄ…第２標準偏差算出部、２Ｅ…流動様式判別部、２Ｆ…第１超音波受信時刻算出部、２Ｇ…第２超音波受信時刻算出部、２Ｈ…第１超音波伝播時間算出部、２Ｊ…第２超音波伝播時間算出部、２Ｋ…伝播時間差算出部、２Ｌ…流量算出部、２Ｍ…流動様式線図取得部、２Ｎ…湿り度算出部、Ａ…配管、Ａ１…壁部、Ａ２…流路、Ｂ…ダンピング材
【図1】
【図2】
【図3】
【図4】
【図5】
【図6】
【図7】
【図8】
【図9】
【図10】
【図11】
【図12】
【図13】