特許 6971332 流体測定装置、流体測定方法、及びプログラム

(19)【発行国】日本国特許庁(JP)

(12)【公報種別】特許公報(B2)

(11)【特許番号】6971332

(24)【登録日】2021年11月4日

(45)【発行日】2021年11月24日

(54)【発明の名称】流体測定装置、流体測定方法、及びプログラム

(51)【国際特許分類】

   G01F 1/66 20060101AFI20211111BHJP

   G01F 1/00 20060101ALI20211111BHJP

   A61B 5/026 20060101ALI20211111BHJP

   A61B 5/0285 20060101ALI20211111BHJP

【ＦＩ】

   G01F1/66 103

   G01F1/00 J

   A61B5/026 120

   A61B5/0285 H

【請求項の数】14

【全頁数】17

(21)【出願番号】特願2019-567126(P2019-567126)

(86)(22)【出願日】2019年1月23日

(86)【国際出願番号】JP2019002148

(87)【国際公開番号】WO2019146663

(87)【国際公開日】20190801

【審査請求日】2020年7月7日

(31)【優先権主張番号】特願2018-8501(P2018-8501)

(32)【優先日】2018年1月23日

(33)【優先権主張国】JP

(73)【特許権者】

【識別番号】000006633

【氏名又は名称】京セラ株式会社

(74)【代理人】

【識別番号】100147485

【弁理士】

【氏名又は名称】杉村 憲司

(74)【代理人】

【識別番号】230118913

【弁護士】

【氏名又は名称】杉村 光嗣

(74)【代理人】

【識別番号】100132045

【弁理士】

【氏名又は名称】坪内 伸

(72)【発明者】

【氏名】長坂 優志

【審査官】 岡田 卓弥

(56)【参考文献】

【文献】 国際公開第２０１５／０３３４６９（ＷＯ，Ａ１）

【文献】 特開２０１５−１３２５８０（ＪＰ，Ａ）

【文献】 特開２０１８−９９２２（ＪＰ，Ａ）

【文献】 米国特許第６０４１６６２（ＵＳ，Ａ）

(58)【調査した分野】（Int.Cl.，ＤＢ名）

Ｇ０１Ｆ １／００− ９／０２

Ｇ０１Ｐ ５／００− ５／２６

Ａ６１Ｂ ５／０２− ５／０３

(57)【特許請求の範囲】

【請求項1】

流体を含む被照射物に向けて光を照射可能な発光部と、
前記流体で散乱した散乱光を受光可能な受光部と、
前記散乱光に基づいて周波数スペクトルを生成する生成部、及び前記周波数スペクトルの特徴成分に基づいて前記流体の流動状態を推定する推定部、を有する制御部と、
既知の流動状態の前記流体に基づく第２周波数スペクトルを取得する取得部と、を備え、
前記制御部は、測定対象の前記流体に基づく第１周波数スペクトルを前記生成部で生成し、前記第２周波数スペクトルを前記取得部により取得した後、前記第１周波数スペクトルの特徴成分と前記第２周波数スペクトルの特徴成分を前記推定部で対比することで、測定対象の前記流体の流動状態を推定する、流体測定装置。

【請求項2】

請求項１に記載の流体測定装置であって、
前記被照射物は、流体が流れる流路、をさらに含み、
前記生成部は、前記流体で散乱した散乱光と、前記流路によって散乱した散乱光とが干渉して生じた干渉光に基づいて、前記周波数スペクトルを生成する、流体測定装置。

【請求項3】

請求項１または２に記載の流体測定装置であって、
前記周波数スペクトルは、前記散乱光のドップラーシフトに基づく信号から生成される、流体測定装置。

【請求項4】

請求項１〜３のいずれかに記載の流体測定装置であって、
前記推定部は、前記周波数スペクトルにおける特定周波数の強度に基づく値及び特定区間の周波数の強度に基づく値の少なくとも一方を前記特徴成分として、前記流体の流動状態を推定する、流体測定装置。

【請求項5】

請求項４に記載の流体測定装置であって、
前記推定部は、少なくとも一つの前記特定周波数に対応する強度と、前記特定周波数以外の少なくとも一つの周波数に対応する強度との比を前記特徴成分として、前記流体の流動状態を推定する、流体測定装置。

【請求項6】

請求項４または５に記載の流体測定装置であって、
前記推定部は、前記周波数スペクトルにおける前記特定区間の強度の平均及び分散の少なくとも一方を前記特徴成分として、前記流体の流動状態を推定する、流体測定装置。

【請求項7】

請求項４〜６のいずれかに記載の流体測定装置であって、
前記推定部は、前記周波数スペクトルの前記特定区間の形状を前記特徴成分として、前記流体の流動状態を推定する、流体測定装置。

【請求項8】

請求項７に記載の流体測定装置であって、
前記推定部は、２直線近似または指数関数近似によって、前記第１周波数スペクトルの形状と第２周波数スペクトルの形状の類否を判断する、流体測定装置。

【請求項9】

請求項１〜８のいずれかに記載の流体測定装置であって、
前記推定部は、前記既知の流動状態と前記第２周波数スペクトルの前記特徴成分との関係を、学習技術に学習させた後、前記第１周波数スペクトルの特徴成分と前記関係に基づいて、測定対象の前記流体の流動状態を推定する、測定装置。

【請求項10】

請求項１〜９のいずれかに記載の流体測定装置であって、
前記推定部は、前記流体の流動状態として、前記流体の流速及び流量の少なくとも一方を推定する、流体測定装置。

【請求項11】

請求項１〜１０のいずれかに記載の流体測定装置であって、
前記発光部は、レーザ光を照射する、流体測定装置。

【請求項12】

流体を含む被照射物に向けて光を照射するステップと、
前記流体で散乱した散乱光を受光するステップと、
前記散乱光に基づいて第１周波数スペクトルを生成するステップと、
既知の流動状態の前記流体に基づく第２周波数スペクトルを取得するステップと、
前記第１周波数スペクトルの特徴成分に基づいて前記流体の流動状態を推定するステップと、を備え、
前記第１周波数スペクトルの特徴成分と前記第２周波数スペクトルの特徴成分を対比することで、測定対象の前記流体の流動状態を推定する、流体測定方法。

【請求項13】

コンピュータに、
流体を含む被照射物に向けて光を照射させるステップと、
前記流体で散乱した散乱光を受光させるステップと、
前記散乱光に基づいて第１周波数スペクトルを生成させるステップと、
既知の流動状態の前記流体に基づく第２周波数スペクトルを取得するステップと、
前記第１周波数スペクトルの特徴成分に基づいて前記流体の流動状態を推定させるステップと、を備え、
前記第１周波数スペクトルの特徴成分と前記第２周波数スペクトルの特徴成分を対比させることで、測定対象の前記流体の流動状態を推定させる、プログラム。

【請求項14】

請求項４に記載の流体測定装置であって、
前記推定部は、前記周波数スペクトルの前記特定区間において積分した強度の値を前記特徴成分として、前記流体の流動状態を推定する、流体測定装置。

【発明の詳細な説明】

【関連出願の相互参照】

【0001】

本出願は、２０１８年１月２３日に日本国に特許出願された特願２０１８−８５０１の優先権を主張するものであり、この先の出願の開示全体を、ここに参照のために取り込む。

【技術分野】

【0002】

本開示は、流体測定装置、流体測定方法、及びプログラムに関する。

【背景技術】

【0003】

従来、流れる物体の流量または流速を計測する装置が知られている。例えば、特許文献１には、光学的に流量を計測する装置が記載されている。

【先行技術文献】

【特許文献】

【0004】

【特許文献1】特開２０１７−１１３３２０号公報

【発明の概要】

【0005】

一実施形態に係る流体測定装置は、流体を含む被照射物に向けて光を照射可能な発光部と、流体で散乱した散乱光を受光可能な受光部と、散乱光に基づいて周波数スペクトルを生成する生成部、及び周波数スペクトルの特徴成分に基づいて流体の流動状態を推定する推定部、を有する制御部と、既知の流動状態の流体に基づく第２周波数スペクトルを取得する取得部と、を備える。そして、制御部は、測定対象の流体に基づく第１周波数スペクトルを生成部で生成し、第２周波数スペクトルを取得部により取得した後、第１周波数スペクトルの特徴成分と第２周波数スペクトルの特徴成分を推定部で対比することで、測定対象の流体の流動状態を推定することができる。

【0006】

一実施形態に係る流体測定方法は、流体を含む被照射物に向けて光を照射するステップと、流体で散乱した散乱光を受光するステップと、散乱光に基づいて第１周波数スペクトルを生成するステップと、既知の流動状態の流体に基づく第２周波数スペクトルを取得するステップと、第１周波数スペクトルの特徴成分に基づいて流体の流動状態を推定するステップと、を備える。そして、一実施形態に係る流体測定方法は、第１周波数スペクトルの特徴成分と第２周波数スペクトルの特徴成分を対比することで、測定対象の流体の流動状態を推定することができる。

【0007】

一実施形態に係るプログラムは、コンピュータに、流体を含む被照射物に向けて光を照射させるステップと、流体で散乱した散乱光を受光させるステップと、散乱光に基づいて第１周波数スペクトルを生成させるステップと、既知の流動状態の流体に基づく第２周波数スペクトルを取得するステップと、第１周波数スペクトルの特徴成分に基づいて流体の流動状態を推定させるステップと、を備える。そして、一実施形態に係るプログラムは、コンピュータに、第１周波数スペクトルの特徴成分と第２周波数スペクトルの特徴成分を対比させることで、測定対象の流体の流動状態を推定させることができる。

【図面の簡単な説明】

【0008】

【図1】一実施形態に係る流体測定装置の概略構成の一例を示すブロック図である。

【図2】一実施形態に係る流体測定装置による干渉光の検出を説明する図である。

【図3】一実施形態に係る流体測定装置が実行するフローチャートの例である。

【図4】一実施形態に係る流体測定の原理を説明するための図である。

【図5】一実施形態に係る流体測定の原理を説明するための図である。

【図6】一実施形態に係る流体測定の原理を説明するための図である。

【図7】一実施形態に係る流体測定の原理を説明するための図である。

【図8】一実施形態に係る流体測定の原理を説明するための図である。

【発明を実施するための形態】

【0009】

種々の条件下で流体の流動状態を測定することができれば、流体測定装置の利便性を高めることができる。本開示は、利便性の高い流体測定装置、流体測定方法、及びプログラムの提供に関する。一実施形態によれば、利便性の高い流体測定装置、流体測定方法、及びプログラムを提供することができる。以下、本開示の一実施形態を、図面を参照して説明する。まず、一実施形態に係る流体測定装置の構成について説明する。

【0010】

図１は、一実施形態に係る流体測定装置の概略構成例を示すブロック図である。

【0011】

図１は、一実施形態に係る流体測定装置が有する機能部を示すブロック図である。なお、図１は、流路Ｂを流れる流体Ａを模式的に示している。また、各機能部に電力を供給する電源、及び電源から各機能部に電力を供給する構成などは省略して示している。

【0012】

一実施形態に係る流体測定装置は、流れる物体（流体）の流れの状態（流動状態）を算出することができる。具体的には、流体測定装置は、流体の流量または流速を流体の流動状態として算出することができる。なお、流量は、単位時間あたりに流れる流体の体積または質量であり、流速は、単位時間あたりに流体が進む距離である。

【0013】

流体測定装置は、光のドップラー効果を利用して流体の流動状態を算出することができる。被照射物（流体及び流体が流れる流路など）に照射される光は、流体によって散乱し、周波数が流体の流動状態に応じてドップラー効果によってシフト（ドップラーシフト）する。そのため、ドップラー効果を利用すれば、流動状態を算出することができる。具体的には、流体測定装置は、測定対象の流体を含む被照射物に向けて発光部から光を照射し、当該被照射物で散乱された光を含む干渉光を受光部で受光することができる。そして、流体測定装置は、受光部の出力に基づいて、流体の流動状態を算出することができる。

【0014】

測定対象の流体は、光のドップラー効果を利用して流動状態を算出可能なものであればよい。具体的には、流体は、それ自体が光を散乱するもの、または光を散乱する物質（散乱物質）を流動させるものであればよい。流体は、例えば、水、血液、プリンター用のインク、または粉体を含む気体などであればよい。なお、散乱物質または粉体などが流体に追従して流動する場合、流体測定装置は、散乱物質または粉体などの流量または流速を、流体の流量または速度とみなすこともできる。すなわち、「流体の流量または流速」とは、「散乱物質または粉体などの流量または流速」と解釈することもできる。

【0015】

一実施形態に係る流体測定装置１は、流体測定装置１の制御を行う制御部１０を備える。一実施形態に係る流体測定装置１において、センサ部６０は、流路Ｂを流れる流体Ａに関する検出を行う。そして、流体測定装置１において、制御部１０は、センサ部６０によって検出された結果に基づいて、流路Ｂを流れる流体Ａの流動状態を推定することができる。

【0016】

センサ部６０の位置は、流路Ｂを流れる流体Ａの流動状態を推定できるように、流路Ｂに対して位置決めされればよい。センサ部６０は、発光部６２と、受光部６４とを備えている。

【0017】

発光部６２は、流路Ｂに向けて光を照射することができる。発光部６２は、例えばレーザ光などを照射することができる。発光部６２は、例えば、流体Ａ又は流体Ａに含まれる固体など、特定成分を検出可能な波長のレーザ光を、測定光として照射してよい。発光部６２は、例えば、任意の数のＬＤ（レーザダイオード：Laser Diode）により構成される。

【0018】

発光部６２は、制御部１０の駆動部５０によって駆動されればよい。駆動部５０は、任意のレーザ駆動回路などにより構成されればよい。なお、駆動部５０は、流体測定装置１の外部に設けてもよいし、センサ部６０に内蔵されてもよい。その結果、流体測定装置１は、設計の自由度を向上させることができる。

【0019】

受光部６４は、発光部６２から照射された光のうち、流体Ａで散乱した光を受光することができる。また、受光部６４は、発光部６２から照射された光のうち、流路Ｂで散乱した光を受光することができる。すなわち、受光部６４は、流体Ａで散乱した光と流路Ｂで散乱した光を含む干渉光を受光することができる。受光部６４は、例えば、任意の数のＰＤ（フォトダイオード：Photo Diode）により構成される。

【0020】

受光部６４によって受光された光に関する信号（光信号）は、制御部１０の生成部１２に送信される。生成部１２については後述する。受光部６４によって受光された光信号を生成部１２に送信して処理する際は、各種の増幅器及び／又はフィルタなどを用いてもよいが、これらは図示を省略してある。

【0021】

なお、センサ部６０は、図１に示すような構成に限定されない。例えば、センサ部６０は、発光部６２及び受光部６４の双方をワンパッケージにして含むものではなく、発光部６２と受光部６４とをそれぞれ別体として構成してもよい。その結果、流体測定装置１は、設計の自由度を向上させることができるため、利便性を向上させることができる。

【0022】

次に、流体測定装置１の制御部１０について説明する。

【0023】

流体測定装置１の制御部１０は、生成部１２と、取得部１４と、推定部１６とを備えている。また、制御部１０は、記憶部２０と、通信部３０と、表示部４０と、駆動部５０との少なくともいずれかを、適宜含んで構成してもよい。

【0024】

制御部１０は、生成部１２、取得部１４、及び推定部１６などの種々の機能を実行するための制御及び処理能力を提供するために、例えばＣＰＵ(Central Processing Unit)のような、少なくとも１つのプロセッサを含む。制御部１０は、生成部１２、取得部１４、及び推定部１６などの機能を、まとめて１つのプロセッサで実現してもよいし、いくつかのプロセッサで実現してもよいし、それぞれ個別のプロセッサで実現してもよい。また、プロセッサは、集積回路（ＩＣ：Integrated Circuit）またはディスクリート回路として実現されればよい。なお、プロセッサは、他の種々の既知の技術に基づいて実現されればよい。一実施形態において、制御部が実行する生成部１２、取得部１４、及び推定部１６などの機能は、例えばＣＰＵ及び当該ＣＰＵで実行されるプログラムとして構成されてもよい。

【0025】

生成部１２は、センサ部６０の受光部６４の出力に基づいて周波数スペクトルを生成し、推定部１６へ出力することができる。取得部１４は、記憶部２０に記憶した周波数スペクトルを取得し、推定部１６へ出力することができる。推定部１６は、生成部１２で生成した周波数スペクトルと、取得部１４で取得した周波数スペクトルに基づいて、流体Ａの流動状態を推定することができる。

【0026】

記憶部２０は、半導体メモリ又は磁気メモリ等で構成されてよい。記憶部２０は、各種情報及び実行されるプログラム等を記憶することができる。記憶部２０は、取得部１４及び／又は推定部１６のワークメモリとして機能してもよい。また、記憶部２０は、周波数スペクトルを記憶することができる。記憶部２０は、生成部１２が生成した周波数スペクトルを予め記憶しておいてもよいし、外部から通信などにより取得した周波数スペクトルを記憶してもよい。なお、記憶部２０は、各種のメモリカードなどであってもよい。

【0027】

通信部３０は、無線通信をはじめとする各種の通信機能を実現することができる。通信部３０は、例えばＬＴＥ（Long Term Evolution）等の種々の通信方式による通信を実現してよい。通信部３０は、例えばＩＴＵ−Ｔ（International Telecommunication Union Telecommunication Standardization Sector）において通信方式が標準化されたモデムを含んでよい。通信部３０は、例えばアンテナを介して、例えば外部サーバ又はクラウドサーバのような外部機器と、ネットワークを介して無線通信してよい。一実施形態において、通信部３０は、例えば外部サーバ又はクラウドサーバなどの外部のデータベースから、第２周波数スペクトルＳ２を受信してよい。また、このようにして通信部３０が受信した第２周波数スペクトルＳ２は、記憶部２０に記憶してもよい。

【0028】

表示部４０は、流体Ａの流動状態の測定結果などの種々の情報を各種の表示デバイスに表示させ、ユーザに通知することができる。表示デバイスは、液晶ディスプレイ（ＬＣＤ）、有機ＥＬディスプレイ、又は無機ＥＬディスプレイ等であればよい。表示部４０は、文字、図形、記号、又はグラフ等の画像を表示させてよい。また、表示部４０は、操作用オブジェクト等の画像を表示させてもよい。

【0029】

また、表示部４０がユーザに通知する情報は、必ずしもユーザに視覚的効果を与えるものに限定されない。例えば、表示部４０は、情報をユーザに伝えることができれば、各種の情報を音としてスピーカから出力させてもよい。

【0030】

制御部１０の記憶部２０、通信部３０、及び表示部４０は、それぞれ流体測定装置１に内蔵されてもよいし、流体測定装置１の外部に設けられてもよい。また、例えば表示部４０は、センサ部６０に内蔵されてもよい。その結果、流体測定装置は、設計の自由度を向上させることができる。

【0031】

なお、流体Ａが流れる流路Ｂは、各種の素材で構成されたチューブ状の部材であればよい。具体的には、流路Ｂは、発光部６２が照射する光の少なくとも一部を透過させる素材で構成されればよい。流路Ｂの素材は、例えば、プラスチック、塩化ビニール、またはガラスなどであればよい。また、流路Ｂは、例えば、血管などの人または動物の体液が流れる体組織であってもよい。なお、流路Ｂは、流体Ａの流動状態が適切に測定されるように、流体Ａが漏れ出したりしないように構成されればよい。

【0032】

次に、一実施形態に係る流体測定における散乱光の検出について説明する。図２は、一実施形態に係る流体測定における散乱光の検出を説明する図である。

【0033】

図２は、一実施形態に係る流体測定装置による干渉光の検出を説明する図である。なお、図２において、流体Ａは、散乱物質を含んでおり、説明の便宜のためいくつかを白の楕円で示してある。また、図２において、散乱物質は速度Ｖで右方向に流動している。

【0034】

発光部６２から流路Ｂに向けて照射される光には、入射光Ｌｅ１及び入射光Ｌｅ２が含まれる。入射光Ｌｅ１及び入射光Ｌｅ２は、発光部６２から照射される際、周波数ｆ０の光とする。なお、周波数ｆ０の入射光Ｌｅ１を、Ｌｅ１（ｆ０）と示し、周波数ｆ０の入射光Ｌｅ２を、Ｌｅ２（ｆ０）と示す。

【0035】

入射光Ｌｅ１は、静止した流路Ｂの表面において界面反射する。すなわち、入射光Ｌｅ１は、流体Ａによって散乱されず、流路Ｂの表面によって散乱される。入射光Ｌｅ１は、流路Ｂの表面によって散乱されて、散乱光Ｌｒ１になる。散乱光Ｌｒ１は、入射光Ｌｅ１が静止している流路Ｂの表面によって散乱されたものである。この時、流路Ｂは静止しており、ドップラー効果が生じないため、入射光Ｌｅ１の周波数ｆ０は変化しない。なお、周波数ｆ０の散乱光Ｌｒ１を、Ｌｒ１（ｆ０）と示す。

【0036】

入射光Ｌｅ２は、流路Ｂの表面において界面反射せずに、流路Ｂの表面を透過する。すなわち、入射光Ｌｅ２は、流路Ｂの表面によって散乱されずに、流体Ａによって散乱される。入射光Ｌｅ２は、流体Ａによって散乱されて散乱光Ｌｒ２になる。この時、入射光Ｌｅ２は、流体Ａによって散乱されているため、周波数ｆ０はドップラーシフトする。なお、周波数ｆ０が周波数Δｆだけ変化した散乱光Ｌｒ２を、Ｌｒ２（ｆ０＋Δｆ）と示す。

【0037】

受光部６４は、上述した散乱光Ｌｒ１及び散乱光Ｌｒ２を受光する。したがって、流体測定装置１は、受光部６４の出力を一要素として、流体Ａの流動状態を推定することができる。

【0038】

図３は、一実施形態に係る流体測定装置１の制御部１０が実行する流動状態推定のフローチャートの例である。

【0039】

まず、生成部１２は、受光部６４からの出力に基づいて、静止している流路Ｂからの散乱光Ｌｒ１と、測定対象の流体Ａからの乱光Ｌｒ２との干渉によって生じる干渉光のうなりの信号（ビート信号）を取得する（ステップＳ１）。ビート信号は、うなりの強度と時間の関係を示すものである。

【0040】

図４は、図３のステップＳ１において取得されるビート信号の一例を示す図である。図４の縦軸は信号出力の強度を表し、横軸は時間を表している。図４に示す例において、受光部６４からの信号出力は電圧値であり、信号出力Ｖｓとして電圧の単位（Ｖ）で示してある。なお、ビート信号は、信号出力の強度の時間変化を示すものであれば、電圧に限定されない。信号出力は、例えば、電流値または抵抗値などであってもよい。

【0041】

次に、生成部１２は、取得したビート信号に基づいて、測定対象の流体に基づく周波数スペクトル（第１周波数スペクトル）Ｓ１を生成する（ステップＳ２）。周波数スペクトルは、ビート信号に含まれる周波数ｆと周波数ごとの強度Ｐ（ｆ）との関係を示したスペクトルである。生成部１２は、ステップＳ１で取得したビート信号を高速フーリエ変換（ＦＦＴ：Fast Fourier Transform）することで、周波数スペクトルを生成することができる。

【0042】

図５は、周波数スペクトルの例を示す図である。図５の縦軸は任意単位（arbitrary unit）の強度Ｐ（ｆ）を表し、横軸は周波数ｆを表している。なお、図５では、流体Ａの流速が２ｍｍ／ｓ、４ｍｍ／ｓ、６ｍｍ／ｓ、８ｍｍ／ｓ、１０ｍｍ／ｓ、及び１４ｍｍ／ｓの場合における周波数スペクトルを例として示している。

【0043】

周波数スペクトルにおいて、流体Ａの流速が速くなるほど、高周波数側の周波数の強度が増加する。流体測定装置１は、このような周波数スペクトルの変化に基づいて、流体の流動状態を推定することができる。

【0044】

推定部１６は、周波数スペクトルＰ（ｆ）に対し、周波数の重み付けを行うことができる（式（１））。すなわち、各周波数ｆとその強度Ｐ（ｆ）の積を算出することができる。

【数1】

推定部１６は、上記式（１）を適当な周波数範囲で積分することにより、以下の積分値を算出することができる（式（２））。

【数2】

推定部１６は、上記式（２）のようにして得られた積分値に比例定数Ｋを乗じる。そして、推定部１６は、レーザ光のような光の受光強度に依存しなくなるように、受光信号の全パワー（Ｉの２乗）、すなわちＤＣ成分でこれを除することで規格化することにより、以下のような値を算出することができる（式（３））。

【数3】

【0045】

推定部１６は、算出した測定対象の流体の式（３）の値と、予め既知の流動状態の流体を測定して算出した式（３）の値を比較することで、測定対象の流体の流動状態を推定することができる。

【0046】

図６は、一実施形態に係る流体測定装置１が推定した流体の流量の例を示す図である。なお、図６は、例として、血管内の血液の流量、すなわち血流量の時間変化を示している。縦軸に示している血流量は、任意単位（arbitrary unit）である。

【0047】

ここで、測定環境等によっては、周波数スペクトルに意図しないノイズが含まれる場合があった。したがって、ノイズの影響を排除せずに測定を行った場合、測定精度が低減していた。これに対し、一実施形態に係る流体測定装置１は、周波数スペクトルＰ（ｆ）に含まれる流体の流動状態ごとに特徴的な成分（特徴成分）に基づいて、流体の流動状態を推定することができる。その結果、流体測定装置１は、種々のノイズの影響を低減することができるため、推定精度を向上させることができる。すなわち、流体測定装置１は、利便性を高めることができる。

【0048】

また、ドップラーシフトの原理を利用した流量計は、レーザ光の発光量のぶれがノイズとして検出されることがあった。したがって、周波数スペクトルで示される強度は、遅い流体の散乱光成分とレーザのぶれの成分との総和となる。すなわち、測定精度が軽減していた。これに対し、一実施形態に係る流体測定装置１は、ノイズの影響を低減させることができるため、推定精度を向上させることができる。

【0049】

また、ＦＦＴ解析を行う際の周波数軸の上限は、サンプリングレートに依存する。このため、ＦＦＴ解析の周波数の上限を超える周波数のドップラーシフト量は観測されなかった。すなわち、流体の流量または流速が大きくなるにつれて、測定精度が低減していた。これに対し、流体測定装置１は、種々の特徴成分に基づいて流体の流動状態を推定することができるため、推定精度を向上させることができる。

【0050】

また、流体の流量の大きさは、流路の断面積にも依存する。すなわち、同じ流量の流体でも、細い流路であれば速く流れ、太い流路であれば遅く流れることになる。このため、流量を測定する際、流路の断面積が不明であれば流量を測定することが困難であった。これに対し、一実施形態に係る流体測定装置１は、必ずしも流路の断面積の情報を用いることなく、流量を推定することができる。すなわち、流体測定装置１は、利便性を高めることができる。

【0051】

また、一実施形態に係る流体測定装置１は、非侵襲に測定することができ、利便性を向上させることができる。

【0052】

以下、一実施形態に係る流体測定装置１による流体測定について、さらに説明する。

【0053】

流体測定装置１の取得部１４は、ステップＳ２において取得した第１周波数スペクトルＳ１との対比の対象となる、既知の流動状態の流体に基づく周波数スペクトル（第２周波数スペクトル）Ｓ２を記憶部２０から取得することができる（ステップＳ３）。この場合、流体測定装置１は、予め既知の流動状態の流体の測定を行い、既知の流動状態と生成部１２で生成した第２周波数スペクトルＳ２との関係を記憶部２０に記憶させておけばよい。また、ステップＳ３において、流体測定装置１は、第２周波数スペクトルＳ２を、通信部３０を介して例えば外部サーバ若しくはクラウドサーバのような外部機器に記憶させてもよい。この場合、取得部１４は、通信部３０を介して、第２周波数スペクトルＳ２を取得すればよい。流体測定装置１は、予め記憶部２０等に第２周波数スペクトルを記憶させておくことで、実測時間を短縮することができる。すなわち、流体測定装置１は、利便性を高めることができる。

【0054】

なお、既知の流動状態は、任意の手法によって把握すればよい。既知の流動状態は、例えば、熱式又は渦電流式などの他の流量計、または設定された流量で流体を流動させることができるポンプなどを用いて把握すればよい。

【0055】

また、流体測定装置１は、周波数スペクトルとともに、周波数スペクトルを実現する各種の測定条件も含めて記憶しておいてもよい。測定条件とは、例えば、サンプリングレート、ＡＤＣ変換レートなどのアルゴリズムの条件、または気温、気圧などの測定環境の条件等であればよい。また、流体測定装置１は、ＬＤなどで構成される発光部６２の強度分散、発光部６２が照射する光の周波数、流路Ｂとセンサ部６０との位置関係、流路Ｂの断面積、流路Ｂを構成する素材・材質、想定される実験環境のノイズなども適宜記録してよい。また、流体測定装置１は、流体の種類、物理・化学特性なども記録してよい。その結果、流体測定装置１は、測定条件の差異に起因する測定誤差の影響を低減することができるため、推定精度を向上させることができる。

【0056】

また、記憶した第２周波数スペクトルＳ２及び当該スペクトルを実現する各種の変数等は、定期的又は不定期に、通信部３０等を介して外部サーバなどから情報の取得を行うことで適宜更新されてもよい。その結果、測定に適した条件を適宜選択することができるため、流体測定装置１は、利便性を向上させることができる。

【0057】

推定部１６は、第１周波数スペクトルＳ１の特徴成分Ｓ１ｃと、第２周波数スペクトルＳ２の特徴成分Ｓ２ｃとを対比する（ステップＳ４）。取得部１４は、第１周波数スペクトルＳ１の特徴成分Ｓ１ｃと対比するのに好適な特徴成分を有する第２周波数スペクトルＳ２を選定して取得することができる。

【0058】

推定部１６は、以下の式（４）の値が最小になるときの第２周波数スペクトルＳ２を採用して、第１周波数スペクトルＳ１と対比することで、流体Ａの流動状態を推定してよい。その結果、流体測定装置１は推定精度を向上させることができる。なお、以下の式（４）において、Ｆ及びｆの一方は第１周波数スペクトルＳ１に基づく値であり、他方は第２周波数スペクトルＳ２に基づく値としてよい。

【数4】

【0059】

また、推定部１６は、上記式（４）に周波数スペクトルの平均周波数の強度を変数として加えてもよい。平均周波数の強度は、以下の式（５）で算出すればよい。

【数5】

【0060】

推定部１６は、上記式（５）に基づいて、第１周波数スペクトルＳ１における平均周波数の強度ｆと、第２周波数スペクトルＳ２における平均周波数の強度Ｆを算出することができる。そして、推定部１６は、以下の式（６）のように、平均周波数の強度ｆと平均周波数の強度Ｆの平均二乗誤差を比較式に加えて、最も値が小さくなる第２周波数スペクトルＳ２対応する流動状態を、測定対象の流動状態として推定してもよい。その結果、流体測定装置１は推定精度を向上させることができる。

【数6】

【0061】

流体測定装置１の制御部１０は、上述のように記憶した周波数スペクトルから、特徴成分を抽出することができる。その結果、流体測定装置１は、測定したい流体ごとに特徴成分を取得することができるため、利便性を高めることができる。なお、記憶部２０は、抽出した特徴成分をさらに記憶してもよい。

【0062】

周波数スペクトルの特徴成分は、例えば次のようなものであればよい。
（１）周波数スペクトルにおける特定周波数の強度の値
制御部１０は、測定する流体ごとに、ノイズが少ない周波数の強度を特徴成分として抽出することができる。これによれば、流体測定装置１は、測定結果のノイズを低減させることができるため、好適な測定が可能となる。ノイズが少ない周波数の強度とは、例えば、周波数ｆとその強度Ｐ（ｆ）の積が最大となる周波数の強度等であればよい。
（２）周波数スペクトルにおける特定周波数の強度と、他の特定周波数の強度との比
これによれば、各周波数の強度に共通するノイズを低減することができる。すなわち、推定部１６は、推定精度を向上させることができるため、好適な測定が可能となる。
（３）周波数スペクトルの強度の平均
これによれば、推定部１６は、測定データを平滑化・規格化することができるため、流動状態ごとの強度の特徴を把握しやすくなる。したがって、推定部１６は、第１周波数スペクトルＳ１と第２周波数スペクトルＳ２との対比が容易となる。また、角周波数の強度に共通するノイズを低減することができるため、推定部１６は、推定精度を向上させることができる。
（４）周波数スペクトルの特定周波数区間において積分した強度の値
推定部１６は、ノイズの少ない周波数区間を選択して演算することができる。これによれば、流体測定装置１は、測定結果のノイズを低減させることができるため、好適な測定が可能となる。また、流動状態に応じて周波数スペクトルの強度は変化するため、これらを積分した値は、流動状態に応じて変化しやすく、流動状態ごとの特徴を反映しやすくなる。すなわち、推定部１６は、推定精度を向上させることができる。
（５）周波数スペクトルの特定周波数区間における強度の平均
これによれば、推定部１６は、測定データを平滑化・規格化することができるため、流動状態ごとの強度の特徴を把握しやすくなる。したがって、推定部１６は、第１周波数スペクトルＳ１と第２周波数スペクトルＳ２との対比が容易となる。また、推定部１６は、ノイズの少ない周波数区間を選択して演算することができるため、ノイズの影響を低減して演算することができる。すなわち、推定部１６は、推定精度を向上させることができる。
（６）周波数スペクトルの特定周波数区間における強度の分散
これによれば、推定部１６は、ノイズの影響を低減して演算することができるため、好適な測定が可能となる。
（７）特定時刻ｔ１から他の特定時刻ｔ２までの区間の周波数スペクトルにおける（２）の値の変化
これによれば、推定部１６は、さらに時間方向の変化に基づく推定が可能となる。すなわち、推定部１６は、流動状態の推定に用いる情報が増えるため、推定精度を向上させることができる。

【0063】

以上のような特徴成分は、第１周波数スペクトルＳ１と第２周波数スペクトルＳ２とが対比可能なパラメータであればよく、上記の数値に限られない。例えば、特徴成分は、パワースペクトルにおける特定区間の形状としてもよい。

【0064】

また、生成部１２は、ビート信号を細かい所定時間で区切り、それぞれの所定時間における周波数スペクトルを生成することができるため、さらに時間の経過を示した三次元の周波数スペクトルを用意してもよい。流体測定装置１は、このような三次元の周波数スペクトルから、特徴成分を抽出してもよい。これによれば、流動状態の変化に伴う特徴成分の推移の傾向から、流動状態の推定が可能となる。すなわち、パワースペクトルの特徴成分ごとの推移の傾向が必ずしも一致する必要はなく、パワースペクトルの全体の推移の傾向が同一又は類似と判断されれば、流動状態の推定が可能となる。したがって、流体測定装置１の有用性を向上可能である。また、推定に用いることができる情報量が多くなるため、流体測定装置１は推定精度を向上させることができる。

【0065】

上述のような特徴成分の抽出は、それぞれ個別に用いるのみならず、それらの２つ以上を適宜組み合わせて用いてもよい。また、例えば、特徴成分は、後述のステップＳ４において用いることができる。

【0066】

次に、推定部１６は、ステップＳ４の対比結果に基づいて、流体Ａの流動状態を推定し（ステップＳ５）、図３に示す動作を終了する。ステップＳ５において推定される流体Ａの流動状態は、例えば、流体Ａの流速及び流量の少なくとも一方としてよい。推定部１６は、第１周波数スペクトルＳ１の特徴成分Ｓ１ｃと同値または近似する第２周波数スペクトルＳ２の特徴成分Ｓ２ｃに対応する流動状態を、測定対象の流動状態と推定することができる。なお、推定部１６が、特徴成分Ｓ１ｃと特徴成分Ｓ２ｃが近似していると判断する範囲は、ユーザの任意で定めればよい。

【0067】

推定部１６は、周波数スペクトルの特定周波数区間の曲線の形状に基づいて、流動状態を推定してもよい。図７及び図８は、係る推定方法を説明するための図である。なお、周波数スペクトルの曲線は、直線も含む概念とする。以下、このような手法について説明する。なお、スペクトルの形状とは、例えば各周波数の強度の最大値同士を結んだ曲線、または強度の中間値同士を結んだ曲線などである。言い換えれば、各周波数の強度の最大値の集合、または中間値の集合などである。

【0068】

推定部１６は、第１周波数スペクトルＳ１に対して２直線近似を行うことで、流動状態を推定してもよい。まず、推定部１６は、任意の低周波数側の周波数から、高周波数側の方向へ直線近似を行う（例えば図７の区間１）。この直線近似は、近似度が劣化するまで行えばよい。ここで、近似度が劣化したか否かは、例えば、カイ２乗検定によって判断すればよい。具体的には、カイ２乗値をその自由度で除したものが、特定閾値を超えた場合、近似度が劣化したとみなすことができる。そして、推定部１６は、第１区間の終点と異なる任意の高周波数側の周波数から、低周波方向へ直線近似を行う（例えば図７の区間２）。この直線近似も、区間１と同様に近似度が劣化するまで行えばよい。そして、推定部１６は、２直線の交点（図７における点Ｘ）を特徴周波数ｆとして算出することができる。

【0069】

次に、推定部１６は、第２周波数スペクトルＳ２に対しても２直線近似を行う。近似の方法は、第１周波数スペクトルＳ１と同様であればよく、推定部１６は、第２周波数スペクトルＳ２において特定周波数Ｆを算出することができる。

【0070】

そして、推定部１６は、特定周波数ｆと最も近い特定周波数Ｆに対応する流動状態を、測定対象の流動状態として推定することができる。

【0071】

推定部１６は、指数関数で近似することで流動状態を算出してもよい。まず、推定部１６は、第１周波数スペクトルＳ１及び第２周波数スペクトルＳ２における特定の周波数区間（例えば２ｋＨｚ〜２０ｋＨｚ）に対して、指数関数で近似を行う（例えば図８における区間３）。次に、推定部１６は、第１周波数スペクトルＳ１を近似して指数関数Ｅ１を取得し、第２周波数スペクトルＳ２を近似して指数関数Ｅ２を取得する。

【0072】

そして、推定部１６は、指数関数Ｅ１のパラメータと最も合致度の高いパラメータを有する指数関数Ｅ２に対応する流動状態を、測定対象の流動状態として推定してもよい。

【0073】

ここで、上記で例示したように抽出した特徴成分は、例えばマシンラーニング又はディープラーニングのような、ＡＩ（Artificial Intelligence）に基づく各種の学習技術を用いて、さらに流動状態との重み付けされてもよい。すなわち、流体測定装置１は、各種の学習技術を用いて、流体の流動状態を算出してもよい。

【0074】

具体的には、推定部１６は、既知の流動状態と第２周波数スペクトルの特徴成分との関係を各種の学習技術に学習させることができる。そして、推定部１６は、第１周波数スペクトルの特徴成分と、各種の学習技術が重み付けを行なった流動状態と特徴成分の関係に基づいて、測定対象の流体の流動状態を算出することができる。これによれば、推定部１６は、第１周波数スペクトルＳ１と第２周波数スペクトルＳ２との対比の精度を向上させることができるため、流体測定装置１は、有用性を向上させることができる。

【0075】

なお、例えば、教師あり学習の場合において、学習データは上記に例示した特徴成分であり、教師データは周波数スペクトルに対応する流体の流動状態であればよい。また、流体測定装置１は、周波数スペクトルを学習データ、周波数スペクトルに対応する流動状態を教師データとして、上記の特徴成分以外の特徴成分に基づいて流動状態を推定してもよい。また、学習で得られた流動状態と特徴成分の関係は、記憶部２０などに適宜記憶されればよい。

【0076】

以上説明したように、特徴成分は、周波数スペクトルにおける少なくとも１つの特定周波数に対応する強度と、前記特定周波数以外の少なくとも１つの周波数に対応する強度との比をとしてもよい。また、特徴成分は、周波数スペクトルにおける特定区間に対応する強度の平均及び分散の少なくとも一方をとしてもよい。また、特徴成分は、周波数スペクトルにおける特定区間のスペクトルの形状としてもよい。

【0077】

なお、ステップ５において推定された流体Ａの流動状態に関する情報は、例えば表示部４０を介して各種の表示デバイスに表示してもよい。これにより、一実施形態に係る流体測定装置１のユーザは、流体Ａの流動状態に関する情報を視認することができる。

【0078】

以上説明したように、一実施形態に係る流体測定装置１によれば、流体の流動状態に係る種々の条件において、流体の流動状態を適切かつ精度よく推定することができる。また、流体測定装置１は、従来のレーザドップラー流量計に対して比較的広範な範囲の流動状態の測定が可能である。すなわち、一実施形態に係る流体測定装置１は、利便性を高めることができる。

【0079】

本開示を諸図面及び実施例に基づき説明してきたが、当業者であれば本開示に基づき種々の変形又は修正を行うことが容易であることに注意されたい。従って、これらの変形又は修正は本開示の範囲に含まれることに留意されたい。例えば、各機能部に含まれる機能などは論理的に矛盾しないように再配置可能である。複数の機能部等は、１つに組み合わせられたり、分割されたりしてよい。上述した本開示に係る各実施形態は、それぞれ説明した各実施形態に忠実に実施することに限定されるものではなく、適宜、各特徴を組み合わせたり、一部を省略したりして実施されうる。

【0080】

また、上述した実施形態は、流体測定装置１としての実施に限定されない。例えば、上述した実施形態は、流体測定装置１において実行されるような流体測定方法、及び流体測定装置１のような装置を制御するコンピュータに実行させるプログラムとして実施してもよい。

【0081】

一実施形態に係る流体測定方法は、流体を含む被照射物に向けて光を照射するステップと、流体で散乱した散乱光を受光するステップと、散乱光に基づいて周波数スペクトルを生成するステップと、周波数スペクトルの特徴成分に基づいて流体の流動状態を推定するステップと、を備える。そして、一実施形態に係る流体測定方法は、測定対象の流体に基づいて第１周波数スペクトル、及び既知の流動状態の流体に基づいて第２周波数スペクトルを生成した後、第１周波数スペクトルの特徴成分と第２周波数スペクトルの特徴成分を対比することで、測定対象の流体の流動状態を推定することができる。

【0082】

一実施形態に係るプログラムは、コンピュータに、流体を含む被照射物に向けて光を照射させるステップと、流体で散乱した散乱光を受光させるステップと、散乱光に基づいて周波数スペクトルを生成させるステップと、周波数スペクトルの特徴成分に基づいて流体の流動状態を推定させるステップと、を備える。そして、一実施形態に係るプログラムは、コンピュータに、測定対象の流体に基づいて第１周波数スペクトル、及び既知の流動状態の流体に基づいて第２周波数スペクトルを生成させた後、第１周波数スペクトルの特徴成分と第２周波数スペクトルの特徴成分を対比させることで、測定対象の流体の流動状態を推定させることができる。

【符号の説明】

【0083】

１ 流体測定装置
１０ 制御部
１２ 生成部
１４ 取得部
１６ 推定部
２０ 記憶部
３０ 通信部
４０ 表示部
５０ 駆動部
６０ センサ部
６２ 発光部
６４ 受光部

【図1】

【図2】

【図3】

【図4】

【図5】

【図6】

【図7】

【図8】