(57)【要約】  血糖値測定方法であって、生体の皮膚から排出された皮膚ガスを検出し、前記皮膚ガスの検出値に基づき、前記生体の血糖値を算出する。
2016072513.tif 000002
【特許請求の範囲】
【請求項１】
血糖値測定方法であって、
生体の皮膚から排出された皮膚ガスを検出し、
前記皮膚ガスの検出値に基づき、前記生体の血糖値を算出する、
血糖値測定方法。
【請求項２】
前記皮膚ガスは、前記生体のグルコース濃度に関連する物質を含む、
請求項１に記載の血糖値測定方法。
【請求項３】
血糖値測定装置であって、
生体の皮膚から排出された皮膚ガスを検出する検出部と、
前記検出部により検出された前記皮膚ガスの検出値に基づき、前記生体の血糖値を算出する算出部と、
を備える血糖値測定装置。
【請求項４】
前記皮膚ガスに含まれる複数の検出対象物質のそれぞれに対応し、前記皮膚ガスを検出する複数の検出部を備え、
前記算出部は、前記複数の検出部により検出された前記複数の検出対象物質の検出値に対し、関数を用いて前記複数の検出部の各々に応じた算術処理を行い、前記算術処理により複数の処理値を得て、前記複数の処理値を加算して得られた値を前記生体の前記血糖値として算出する、
請求項３に記載の血糖値測定装置。
【請求項５】
前記算術処理は、第１算術処理と、第２算術処理を含み、
前記第１算術処理を行った後に、前記第２算術処理を行う、
請求項４に記載の血糖値測定装置。
【請求項６】
前記算術処理においては、前記複数の検出対象物質の前記検出値の各々に、重み付け係数を乗算する、
請求項４または５に記載の血糖値測定装置。
【請求項７】
前記算術処理においては、前記複数の検出対象物質の前記検出値を、前記複数の検出部のうち所定の検出部により検出された検出値によって除算する、
請求項４から６のうちいずれか１項に記載の血糖値測定装置。
【請求項８】
前記算術処理においては、前記複数の検出対象物質の前記検出値をそれぞれ二乗する、
請求項４から６のうちいずれか１項に記載の血糖値測定装置。
【請求項９】
前記算出部により算出された前記生体の前記血糖値と、前記検出部によりそれぞれ検出された前記複数の検出対象物質の前記検出値とに基づき、前記生体の健康状態を評価する評価部をさらに備える、
請求項３から８のうちいずれか１項に記載の血糖値測定装置。
【請求項１０】
前記生体において血管が集中する前記生体の皮膚から、前記皮膚ガスを採取する採取部と、
前記採取部から前記複数の検出部に至るまで設けられた複数の搬送経路と、
前記複数の搬送経路を介して、前記採取部により採取された前記皮膚ガスを、前記複数の検出部に搬送する搬送部と、
前記採取部により前記皮膚ガスが採取されない期間において、前記複数の検出部を洗浄するように前記搬送部を制御する搬送制御部と、をさらに備える、
請求項４から９のうちいずれか１項に記載の血糖値測定装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、血糖値測定方法、および血糖値測定装置に関する。
本願は、２０１４年１１月７日に日本に出願された特願２０１４−２２７５３９号に基づき優先権を主張し、その内容をここに援用する。
【背景技術】
【０００２】
生体中の血糖値の測定は、穿刺針によって生体中の血液や間質液等を採取し測定する手法が一般的である。しかしながら、上記手法の場合、穿刺針が生体内に侵襲するため、利用者の負担が大きくなる場合があった。【０００３】
これに対し、生体中の血糖値を非侵襲で測定する技術として、皮膚を透過する近赤外の光を利用した赤外分光法に関する技術が知られている（例えば、特許文献１、２参照）。
一方、光を利用せず生体中の物質を非侵襲で測定する技術として、皮膚から排出される皮膚ガスを分析することにより、生体情報を測定する技術が知られている（例えば、特許文献３〜６）。
【先行技術文献】
【特許文献】
【０００４】
【特許文献１】日本国特開２００３−０４２９４８号公報
【特許文献２】日本国特開２００８−２５６３９８号公報
【特許文献３】日本国特開２００９−１７５１１１号公報
【特許文献４】日本国特開２００６−２３４８４５号公報
【特許文献５】日本国特開２０１０−２８１６９８号公報
【特許文献６】日本国実用新案登録第３１７５４６０号公報
【発明の概要】
【発明が解決しようとする課題】
【０００５】
しかしながら、特許文献１に記載の技術では、検出対象ごとに対応した波長範囲の特定が困難な場合や、非検出対象によって測定時に照射する光が吸収され、その反射光の強度が低下することにより誤差が多くなる場合があった。また、特許文献２に記載の技術では、検出対象として例えばグルコースが、非検出対象であるマルトース等と類似または一致した特性を示すために、両対象（グルコースとマルトース等）を区別することが困難な場合があった。
また、特許文献３〜６に記載の技術には、生体中の血糖値を測定する具体的な方法について記載がされていない。【０００６】
本発明は、このような事情を考慮してなされたものであり、非侵襲で生体中の血糖値を精度良く測定することができる血糖値測定方法、および血糖値測定装置を提供することを目的の一つとする。
【課題を解決するための手段】
【０００７】
上述した課題を解決する本発明の第一態様に係る血糖値測定方法は、生体の皮膚から排出された皮膚ガスを検出し、前記皮膚ガスの検出値に基づき、前記生体の血糖値を算出する。【０００８】
前記皮膚ガスは、前記生体のグルコース濃度に関連する物質を含んでいてもよい。【０００９】
また、本発明の第二態様に係る血糖値測定装置は、生体の皮膚から排出された皮膚ガスを検出する検出部と、前記検出部により検出された前記皮膚ガスの検出値に基づき、前記生体の血糖値を算出する算出部とを備える。【００１０】
上記態様に係る血糖値測定装置は、前記皮膚ガスに含まれる複数の検出対象物質のそれぞれに対応し、前記皮膚ガスを検出する複数の検出部を備え、前記算出部は、前記複数の検出部により検出された前記複数の検出対象物質の検出値に対し、関数を用いて前記複数の検出部の各々に応じた算術処理を行い、前記算術処理により複数の処理値を得て、前記複数の処理値を加算して得られた値を前記生体の前記血糖値として算出してもよい。【００１１】
前記算術処理は、第１算術処理と、第２算術処理を含み、前記第１算術処理を行った後に、前記第２算術処理を行ってもよい。【００１２】
前記算術処理においては、前記複数の検出対象物質の前記検出値の各々に、重み付け係数を乗算してもよい。【００１３】
前記算術処理においては、前記複数の検出対象物質の前記検出値を、前記複数の検出部のうち所定の検出部により検出された検出値によって除算してもよい。【００１４】
前記算術処理においては、前記複数の検出対象物質の前記検出値をそれぞれ二乗してもよい。【００１５】
上記態様に係る血糖値測定装置は、前記算出部により算出された前記生体の前記血糖値と、前記検出部により検出された前記複数の検出対象物質の検出値とに基づき、前記生体の健康状態を評価する評価部をさらに備えていてもよい。【００１６】
上記態様に係る血糖値測定装置は、前記生体において血管が集中する前記生体の皮膚から、前記皮膚ガスを採取する採取部と、前記採取部から前記複数の検出部に至るまで設けられた複数の搬送経路と、前記複数の搬送経路を介して、前記採取部により採取された前記皮膚ガスを、前記複数の検出部に搬送する搬送部と、前記採取部により前記皮膚ガスが採取されない期間において、前記複数の検出部を洗浄するように前記搬送部を制御する搬送制御部と、をさらに備えていてもよい。
【発明の効果】
【００１７】
本発明の上記態様によれば、生体の皮膚から排出される皮膚ガスを検出し、検出した皮膚ガスのガス濃度に基づき、生体中の血糖値を算出することにより、非侵襲で生体中の血糖値を精度良く測定することができる。
【図面の簡単な説明】
【００１８】
【図１】第１実施形態に係る血糖値測定装置１の一例を示す概略図である。
【図２】第１実施形態に係る制御装置３００の機能構成例を示す図である。
【図３】第１実施形態に係る血糖値測定装置１における皮膚ガスの供給経路の一例を示す図である。
【図４Ａ】予め取得された実測血糖値と、当該実測血糖値から得られた推定式を用いて所定周期ごとにガス測定から算出された算出血糖値と、を同グラフ上にプロットした図である。
【図４Ｂ】予め取得された実測血糖値と、当該実測血糖値から得られた推定式を用いて複数の人について所定周期ごとにガス測定から算出された算出血糖値と、を同グラフ上にプロットした図である。
【図５】第１実施形態に係る血糖値測定装置１により実行される処理の流れの一例を示すフローチャートである。
【図６】利用者の健康状態を評価する際に参照される参照テーブルの一例を示す図である。
【図７】アルコールＲＯＨが検出された利用者の健康状態を評価する際に参照される参照テーブルの一例を示す図である。
【発明を実施するための形態】
【００１９】
以下、本発明の実施形態について図面を参照して説明する。
＜第１実施形態＞
図１は、第１実施形態に係る血糖値測定装置１の一例を示す概略図である。血糖値測定装置１は、皮膚ガスプローブ１００と、皮膚ガス検出ユニット２００と、制御装置３００とを備える。
【００２０】
皮膚ガスプローブ１００は、容器内部の気体を容器外部に流出させない材質（例えばテトラフロロポリエチレン）で構成された円筒状の容器である。皮膚ガスプローブ（採取部）１００は、例えば、容器内部に挿入された指Ｆの皮膚表面から排出されるガス（以下、「皮膚ガス」と称する）を採取する。
【００２１】
皮膚ガス検出ユニット２００は、第１ポンプ２０２と、排気口２０４と、吸気口２０６と、乾燥器２０８と、第２ポンプ２１０と、電磁バルブ（搬送部）Ｖ１〜Ｖ１４と、センサ（検出部）Ｓａ〜Ｓｇと、ガス供給パイプ（搬送経路）Ｐとを備える。
第１ポンプ２０２および第２ポンプ２１０は、後述する制御装置３００と通信可能なように有線または無線等によって接続されている。第１ポンプ２０２および第２ポンプ２１０は、制御装置３００から送信される制御信号に基づき運転または停止する。第１ポンプ２０２は、運転時、皮膚ガス検出ユニット２００のガス供給パイプＰ内部の皮膚ガスを、排気口２０４に排気する。また、第２ポンプ２１０は、運転時、吸気口２０６から乾燥器２０８を介して皮膚ガス検出ユニット２００外部の空気を吸気し、吸気した空気を電磁バルブＶ１側へ供給する。以下、乾燥器２０８を介して皮膚ガス検出ユニット２００に供給される空気を、「乾燥空気」と記載する。
【００２２】
乾燥器２０８は、例えば、シリカゲル等の乾燥剤を内部に有し、通過する気体の水分を乾燥剤に吸収させて乾燥させる。通過する気体の湿度を低く一定に保つことにより、皮膚より蒸泄される水分量を精度よく検知することができる。
【００２３】
電磁バルブＶ１〜Ｖ１４は、制御装置３００と通信可能なように有線または無線等によって接続されている。電磁バルブＶ１〜Ｖ１４は、皮膚ガスの通路を開閉することのできる可動機構（弁）を内部に有しており、制御装置３００から送信される制御信号に基づき、可動機構（弁）を稼働させ、気体の流れる方向を切替える。
【００２４】
センサＳａ〜Ｓｇは、制御装置３００と通信可能なように有線または無線等によって接続されている。センサＳａ〜Ｓｇは、皮膚ガスに含まれる物質のうち、生体中のグルコース濃度に関連する物質ごとに設けられ、グルコース濃度に関連する物質を検出する半導体ガスセンサである。生体中のグルコース濃度に関連する物質は、二酸化炭素ＣＯ２と、一酸化炭素ＣＯと、アルコールＲＯＨ（例えばエタノール）と、アセトアルデヒドＲＣＨＯと、酢酸ＲＣＯＯＨと、アセトンＲＣＯＲと、水Ｈ２Ｏとを含む。本実施形態において、グルコース濃度に関連する物質は、気体で存在しているとする。
【００２５】
センサＳａ〜Ｓｇは、例えば、気体中に含まれる特定の物質を検出する金属酸化物がコーティングされた半導体表面（以下、「ガス検出面」と称する）と、抵抗部（例えば白金コイル）と、抵抗部の抵抗値を測定する抵抗測定部とを有する。供給された皮膚ガスがセンサＳａ〜Ｓｇ内部を通過する際、ガスに含まれる特定の物質は、ガス検出面に吸着する。特定の物質がガス検出面に吸着することによって、抵抗部に流れる電流の電気伝導率が変化し、抵抗測定部が測定する抵抗値が変化する。これによって、センサＳａ〜Ｓｇは、グルコース濃度に関連する物質を検出することができる。なお、センサＳａ〜Ｓｇの数は、生体中のグルコース濃度に関連する物質の種類に応じて増減させてもよく、センサＳａ〜Ｓｇの７つの数に限られない。
【００２６】
また、本実施形態におけるセンサＳａ〜Ｓｇは、特定の物質のみを検出することができず、複数の物質に対して所定の検出感度を有している。センサＳａは、例えば二酸化炭素ＣＯ２に対して最も検出感度が高いガス検出面を有し、センサＳｂは、例えば一酸化炭素ＣＯに対して最も検出感度が高いガス検出面を有し、センサＳｃは、例えばアルコールＲＯＨに対して最も検出感度が高いガス検出面を有している。また、センサＳｄは、例えばアセトアルデヒドＲＣＨＯに対して最も検出感度が高いガス検出面を有し、センサＳｅは、例えば、酢酸ＲＣＯＯＨに対して最も検出感度が高いガス検出面を有している。また、センサＳｆは、例えばアセトンＲＣＯＲに対して最も検出感度が高いガス検出面を有し、センサＳｇは、例えば水Ｈ２Ｏに対して最も検出感度が高いガス検出面を有している。
【００２７】
これらセンサＳａ〜Ｓｇにおいて、皮膚ガス検出ユニット２００内における前段から後段へ向かう順番や配置位置等は、上述したものに限られない。センサＳａ〜Ｓｇの順番や配置位置等は、例えば、物質検出時にセンサ内部で消費される皮膚ガスの割合を考慮して、皮膚ガスの消費量が少ないセンサを前段に、皮膚ガスの消費量が多いセンサを後段にするようにしてもよい。これによって、血糖値測定装置１は、皮膚ガスプローブ１００により採取される皮膚ガスが少量であっても、血糖値の測定を精度良く行うことができる。
【００２８】
ガス供給パイプＰは、皮膚ガスおよび乾燥空気を、各電磁バルブＶ１〜Ｖ１４を介し、それぞれのセンサＳａ〜Ｓｇに供給するために設けられるガス供給路である。
【００２９】
制御装置３００は、皮膚ガス検出ユニット２００の各種機能部を制御するコンピュータ装置である。以下、制御装置３００について図を参照して説明する。
図２は、第１実施形態に係る制御装置３００の機能構成例を示す図である。制御装置３００は、通信部３０２と、表示部３０４と、制御部３１０と、記憶部３２０とを備える。
【００３０】
通信部３０２は、皮膚ガス検出ユニット２００の各種機能部と無線または有線によって通信を行うための通信インターフェースを含む。また、通信部３０２は、制御装置３００の外部の他装置とネットワークを介して通信を行ってもよい。ネットワークは、例えば、インターネットや電話回線等の通信回線を含む。
【００３１】
表示部３０４は、例えば、ＬＣＤ（Liquid  Crystal  Display）や有機ＥＬ（Electroluminescence）ディスプレイ等の表示装置を含む。表示部３０４は、後述する算出部３１６の算出結果や評価部３１８の評価結果等を含む情報を表示する。また、表示部３０４は、表示装置の表示面を検知面とする接触検知機構とを含むタッチパネルであってもよい。
この場合、表示部３０４は、利用者からの操作入力を受け付けるユーザインターフェースとして機能する。
【００３２】
制御部３１０は、例えば、ＣＰＵ（Central  Processing  Unit）等のプロセッサが記憶部３２０に格納されたプログラムを実行することにより機能するソフトウェア機能部である。また、これらの制御部３１０の各機能部のうち一部または全部は、ＬＳＩ（Large  Scale  Integration）やＡＳＩＣ（Application  Specific  Integrated  Circuit）等のハードウェア機能部であってもよい。制御部３１０は、ガス供給制御部（搬送制御部）３１２と、取得部３１４と、算出部３１６と、評価部３１８とを備える。
【００３３】
ガス供給制御部３１２は、所定周期（例えば１０分）ごとに、第２ポンプ２１０から供給される乾燥空気を、皮膚ガスプローブ１００側（供給経路Ａ１）、または電磁バルブＶ３側（供給経路Ｂ１）へ供給するように電磁バルブＶ１を制御する。また、ガス供給制御部３１２は、第２ポンプ２１０から供給される乾燥空気を、皮膚ガスプローブ１００側（供給経路Ａ１）へ供給するように電磁バルブＶ１を制御する際、皮膚ガスプローブ１００から供給される皮膚ガスを、電磁バルブＶ３側（供給経路Ａ３）に供給するように電磁バルブＶ２を制御する。また、ガス供給制御部３１２は、電磁バルブＶ２を介して供給された皮膚ガスを電磁バルブＶ４側へ供給するように電磁バルブＶ３を制御する。これによって、ガス供給制御部３１２は、皮膚ガスプローブ１００により採取された皮膚ガスを、電磁バルブＶ２、Ｖ３を介して、後段の各種センサＳａ〜Ｓｇに供給する。
【００３４】
一方、ガス供給制御部３１２は、第２ポンプ２１０から供給される乾燥空気を、電磁バルブＶ３側（供給経路Ｂ１）へ供給するように電磁バルブＶ１を制御する際、電磁バルブＶ１から供給された乾燥空気を、電磁バルブＶ２、Ｖ４側へ供給するように電磁バルブＶ３を制御する。また、ガス供給制御部３１２は、皮膚ガスプローブ１００から皮膚ガスが供給される供給経路Ａ２を遮断しつつ、電磁バルブＶ３を介して供給された乾燥空気を電磁バルブＶ１４側（供給経路Ｂ２）へ供給するように電磁バルブＶ２を制御する。これによって、ガス供給制御部３１２は、皮膚ガスプローブ１００により採取された皮膚ガスを各種センサＳａ〜Ｓｇに供給しない期間、乾燥空気を各種センサＳａ〜Ｓｇに供給する。
この結果、血糖値測定装置１は、血糖値を測定しない期間において、各種センサＳａ〜Ｓｇのガス検出面に常に乾燥した空気を供給することにより、各種センサＳａ〜Ｓｇを洗浄することができ、皮膚ガスの検出精度を向上させることができる。
【００３５】
また、ガス供給制御部３１２は、第２ポンプ２１０を運転するように制御している期間、第１ポンプ２０２も運転するように制御する。これによって、血糖値測定装置１は、各種センサＳａ〜Ｓｇに供給された皮膚ガスや乾燥空気を排気口２０４へ排気する。
【００３６】
取得部３１４は、皮膚ガスに含まれる物質のガス濃度（例えば単位は［ｐｐｍ］）に関するデータ（検出値）を、センサ（検出部）Ｓａ〜Ｓｇからそれぞれ取得する。以下、センサＳａ（第１検出部）から取得された物質のガス濃度に関するデータを、「センサデータａ（第１検出値）」と称し、センサＳｂ（第２検出部）から取得された物質のガス濃度に関するデータを、「センサデータｂ（第２検出値）」と称する。また、センサＳｃ（第３検出部）から取得された物質のガス濃度に関するデータを、「センサデータｃ（第３検出値）」と称し、センサＳｄ（第４検出部）から取得された物質のガス濃度に関するデータを、「センサデータｄ（第４検出値）」と称し、センサＳｅ（第５検出部）から取得された物質のガス濃度に関するデータを、「センサデータｅ（第５検出値）」と称する。また、センサＳｆ（第６検出部）から取得された物質のガス濃度に関するデータを、「センサデータｆ（第６検出値）」と称し、センサＳｇ（第７検出部）から取得された物質のガス濃度に関するデータを、「センサデータｇ（第７検出値）」と称する。
【００３７】
算出部３１６は、取得部３１４により取得された各種センサデータａ〜ｇ（第１〜第７検出値）の各々の値の特徴に基づき、生体中の血糖値としてグルコースＧｌｃの濃度を算出する。以下、図および数式を参照して、グルコースＧｌｃの濃度の具体的な算出方法について説明する。
【００３８】
図３は、第１実施形態に係る血糖値測定装置１における皮膚ガスの供給経路の一例を示す図である。
算出部３１６は、取得部３１４により取得された各種センサデータａ〜ｇの各々の値の特徴に基づき、生体中のグルコースＧｌｃの濃度を血糖値として算出する。そのため、まず算出部３１６は、各種センサデータａ〜ｇの各々の値の特徴に基づき、グルコース濃度に関連する物質ごとの濃度を算出する。例えば、二酸化炭素ＣＯ２の濃度［ＣＯ２］は、式（１）により算出され、一酸化炭素ＣＯの濃度［ＣＯ］は、式（２）により算出され、アルコールＲＯＨの濃度［ＲＯＨ］は、式（３）により算出される。式（１）〜（３）に示すａ〜ｇは、センサデータａ〜ｇ（第１〜第７検出値）にそれぞれ対応する。
【００３９】
［ＣＯ２］＝Ａ１・ａ＋Ｂ１・ｂ＋Ｃ１・ｃ＋Ｄ１・ｄ＋Ｅ１・ｅ＋Ｆ１・ｆ＋Ｇ１・ｇ                                                       （１）
【００４０】
［ＣＯ］＝Ａ２・ａ＋Ｂ２・ｂ＋Ｃ２・ｃ＋Ｄ２・ｄ＋Ｅ２・ｅ＋Ｆ２・ｆ＋Ｇ２・ｇ                                                          （２）
【００４１】
［ＲＯＨ］＝Ａ３・ａ＋Ｂ３・ｂ＋Ｃ３・ｃ＋Ｄ３・ｄ＋Ｅ３・ｅ＋Ｆ３・ｆ＋Ｇ３・ｇ                                                        （３）
【００４２】
式（１）〜（３）に示す各種センサデータａ〜ｇに乗算される係数Ａ１〜３〜Ｇ１〜３（重み付け係数）は、予め実験やシミュレーション等によって求めた理想的な数値と、各種センサデータａ〜ｇと、が合致する場合に演算結果としての相関が最も高くなるように決定される。なお、これら係数Ａ１〜３〜Ｇ１〜３は、実験やシミュレーション等による数値に応じて、正または負の値を取る。また、グルコース濃度に関連するその他物質の濃度算出式（アセトアルデヒドＲＣＨＯ濃度［ＲＣＨＯ］、酢酸ＲＣＯＯＨ濃度［ＲＣＯＯＨ］、アセトンＲＣＯＲ濃度［ＲＣＯＲ］、水Ｈ２Ｏ濃度［Ｈ２Ｏ］等の算出式）については省略する。また、省略したその他物質の濃度算出式中の各係数（例えば、［ＲＣＨＯ］に係るＡ４〜Ｇ４、［ＲＣＯＯＨ］に係るＡ５〜Ｇ５、［ＲＣＯＲ］に係るＡ６〜Ｇ６、［Ｈ２Ｏ］に係るＡ７〜Ｇ７のような重み付け係数）においても、上述した方法と同じように決定される。
このように、グルコース濃度に関連する物質の濃度（例えば、［ＣＯ２］、［ＣＯ］、［ＲＯＨ］、［ＲＣＨＯ］、［ＲＣＯＯＨ］、［ＲＣＯＲ］、［Ｈ２Ｏ］）は、各センサデータａ〜ｇのいずれのセンサで検出されたかによって変化する係数（重み付け係数）Ａ１〜７〜Ｇ１〜７が乗算されるような関数を用いた算術処理（重み付け、重み付け計算）が行われることにより決定され、処理値として得られる。
【００４３】
算出部３１６は、算出したグルコース濃度に関連する物質ごとの濃度に基づき、グルコースＧｌｃの濃度［Ｇｌｃ］を算出する。算出部３１６は、例えば、式（４）に基づきグルコースＧｌｃの濃度を算出する。以下に示すように、グルコース濃度に関連する物質の濃度（処理値、第１の処理値）に対して、係数（第２の重み付け係数）Ａ＃〜Ｇ＃が乗算されるように、関数を用いた算術処理（第２の算術処理、重み付け、重み付け計算）が行われることによりさらに処理値（第２の処理値、例えばＡ＃・［ＣＯ２］、Ｂ＃・［ＣＯ］、Ｃ＃・［ＲＯＨ］等）が得られる。第２の処理値の和としてグルコース濃度は決定される。
【００４４】
［Ｇｌｃ］＝Ａ＃・［ＣＯ２］＋Ｂ＃・［ＣＯ］＋Ｃ＃・［ＲＯＨ］＋Ｄ＃・［ＲＣＨＯ］＋Ｅ＃・［ＲＣＯＯＨ］＋Ｆ＃・［ＲＣＯＲ］＋Ｇ＃・［Ｈ２Ｏ］                                                          （４）
【００４５】
式（４）に示す係数Ａ＃〜Ｇ＃は、上述した濃度算出式中の各種センサデータａ〜ｇに乗算される係数Ａ１〜７〜Ｇ１〜７に基づき決定される。例えば、係数Ａ＃、Ｂ＃、Ｃ＃は、それぞれ式（５）〜（７）により算出される。
【００４６】
Ａ＃＝Ｋ（Ａ１＋Ａ２＋Ａ３＋Ａ４＋Ａ５＋Ａ６＋Ａ７）      （５）
【００４７】
Ｂ＃＝Ｌ（Ｂ１＋Ｂ２＋Ｂ３＋Ｂ４＋Ｂ５＋Ｂ６＋Ｂ７）      （６）
【００４８】
Ｃ＃＝Ｍ（Ｃ１＋Ｃ２＋Ｃ３＋Ｃ４＋Ｃ５＋Ｃ６＋Ｃ７）      （７）
【００４９】
ここで、上記式（５）〜（７）中の係数Ｋ、Ｌ、Ｍは、グルコース濃度［Ｇｌｃ］に関連する物質が、グルコースＧｌｃの濃度［Ｇｌｃ］に対して寄与する割合（以下、「寄与率」と称する）に応じて決定される。以下、グルコースＧｌｃの濃度［Ｇｌｃ］に対する物質ごとの寄与率の決定方法について、生体内で生成される二酸化炭素ＣＯ２、水Ｈ２Ｏ、およびアセトンＲＣＯＲを例にして説明する。なお、その他物質（一酸化炭素ＣＯ、アルコールＲＯＨ、アセトアルデヒドＲＣＨＯ、酢酸ＲＣＯＯＨ等）の寄与率の決定方法については省略する。
【００５０】
生体中の代謝活動により、１個のグルコース分子は、細胞気質内で２つのピルビン酸分子に変換される。これらピルビン酸分子は、細胞内のミトコンドリアのクエン酸回路に運ばれ、例えば、２個の二酸化炭素ＣＯ２分子、１０個のＮＡＤＨ（ニコチンアミドアデニンジヌクレオチド）分子、２個のＦＡＤＨ２（フラビンアデニンジヌクレオチド）分子、その他分子（例えばＧＴＰ（グアノシン三リン酸））等に変換される。クエン酸回路内で変換された１０個のＮＡＤＨ分子および２個のＦＡＤＨ２分子は、同じミトコンドリア内の電子伝達系に運ばれる。電子伝達系において、１０個のＮＡＤＨ分子は、１０個のＨ２Ｏ分子とその他分子とに変換され、２個のＦＡＤＨ２分子は、２個のＨ２Ｏ分子とその他分子とに変換される。すなわち、１個のグルコース分子は、２個の二酸化炭素ＣＯ２分子と、１２個のＨ２Ｏ分子とに変換される。これによって、１個あたりのグルコース分子から生成される二酸化炭素ＣＯ２の分子数とＨ２Ｏの分子数とに応じて、それぞれ二酸化炭素ＣＯ２の寄与率とＨ２Ｏの寄与率とを決定する。
【００５１】
また、健康状態における生体中のグルコースＧｌｃの濃度に比して低いグルコースＧｌｃの濃度を示す状態（例えば空腹や飢餓状態）の場合、生体中の脂肪は、所定の代謝経路によって、脂肪酸とグリセロールとを含む物質に変換される。各脂肪酸分子は、生体中のβ酸化と呼ばれる代謝経路により、１個のアセト酢酸分子（アセチルＣＯＡ）とその他分子とに変換される。その後同じ代謝経路において、変換された１個のアセト酢酸分子（アセチルＣＯＡ）は、１個のアセトンＲＣＯＲ分子と１個の二酸化炭素ＣＯ２分子とに変換される。すなわち、生体中において、グルコースＧｌｃの濃度の低下と、脂肪酸またはアセト酢酸（アセチルＣＯＡ）の分子量の増加と、の対応関係が予め取得されている場合、グルコースＧｌｃの濃度に対するアセトンＲＣＯＲの寄与率と二酸化炭素ＣＯ２の寄与率とを決定できる。この対応関係は、例えば、複数人から生体採血により取得した血液のデータに基づき算出される。なお、上述した生体中の代謝経路によって生成される各物質の分子数は一例であるため、分子数の変更に伴い寄与率を適宜変更してもよい。
【００５２】
評価部３１８は、予め生体採血により取得された血糖値（以下、「実測血糖値」と称する）を示すデータに基づき、算出部３１６により算出された血糖値（以下、「算出血糖値」と称する）の評価を行い、評価結果を表示部３０４に出力する。実測血糖値は、例えば、本実施形態における血糖値測定装置１を利用する利用者から予め取得しておくものとする。また、利用者の年齢や性別、疾患履歴等の医療データと近しい医療データを有する他の利用者がいる場合、その利用者から取得した実測血糖値を当該利用者の実測血糖値の代わりに用いてもよい。以下、図を参照して、具体的な血糖値の評価方法について説明する。
【００５３】
図４Ａは、予め取得された実測血糖値と、当該実測血糖値から得られた推定式を用いて所定周期ごとにガス測定から算出された算出血糖値と、を同グラフ上にプロットした図である。なお、図４Ａ中に示す１５個の黒点の群と１０個の白点の群は、各々別の日の測定データ群である。
図４Ｂは、図４Ａと同様に、予め取得された実測血糖値と、当該実測血糖値から得られた推定式を用いて所定周期ごとにガス測定から算出された算出血糖値と、を同グラフ上にプロットした図であり、図４Ａに示したデータ群とは異なる日の複数の人に関する測定データ群である。
図４Ａおよび図４Ｂにおいて、縦軸および横軸の単位は、［ｍｇ／ｄＬ］である。
【００５４】
まず、評価部３１８は、予め取得された実測血糖値を示すデータに対して最小二乗法等の統計処理を行い、実測血糖値の示す最適な近似曲線Ｌ（例えば一次関数）を導出する。
次に、評価部３１８は、所定周期ごとに算出された算出血糖値が、導出した近似曲線Ｌに対して所定程度以上の強い相関を示す場合、強い相関を示す算出血糖値を表示部３０４に出力する。強い相関とは、例えば、近似曲線Ｌと算出血糖値との相関係数が０．８以上を示す場合である。評価部３１８は、所定周期ごとに算出された算出血糖値が、導出した近似曲線Ｌに対して強い相関を示さない場合、エラーを示す情報等を表示部３０４に出力する。
【００５５】
記憶部３２０は、ＲＯＭ（Read  Only  Memory）やＲＡＭ（Random  Access  Memory）、ＨＤＤ（Hard  Disk  Drive）、ＳＳＤ（Solid  State  Drive）、ＥＥＰＲＯＭ（Electrically  Erasable  Programmable  Read-Only  Memory）、フラッシュメモリ等の記憶装置を含む。記憶部３２０は、実測血糖値、算出血糖値、各種センサデータａ〜ｇ等を記憶するように制御される。なお、記憶部３２０は、制御装置３００に内蔵されるものに代えて、ＮＡＳ（Network  Attached  Storage）装置などの外付け型の不揮発性記憶装置でもよい。
【００５６】
図５は、第１実施形態に係る血糖値測定装置１により実行される処理の流れの一例を示すフローチャートである。本フローチャートは、例えば、所定周期で繰り返し実行される。まず、血糖値測定装置１の運転は、利用者が血糖値を測定するか否かに応じて決定される（ステップＳ１００）。利用者が血糖値を測定しない場合（ステップＳ１００：Ｎｏ）、血糖値測定装置１は、各種センサＳａ〜Ｓｇのガス検出面に乾燥した空気を供給して各種センサＳａ〜Ｓｇを洗浄する（ステップＳ１０２）。利用者が血糖値を測定する場合（ステップＳ１００：Ｎｏ）、利用者は皮膚ガスプローブ１００に指を挿入する。次に、皮膚ガスプローブ１００は、挿入された指から排出される皮膚ガスを採取する（ステップＳ１０４）。
【００５７】
次に、ガス供給制御部３１２は、皮膚ガスプローブ１００により採取された皮膚ガスを、電磁バルブＶ２、Ｖ３を介して、後段の各種センサＳａ〜Ｓｇに供給する（ステップＳ１０６）。次に、センサＳａ〜Ｓｇは、皮膚ガスに含まれる物質のうち、生体中のグルコース濃度に関連する物質を検出する。次に、取得部３１４は、各センサＳａ〜Ｓｇにより検出された皮膚ガスに含まれる物質のガス濃度に関するデータを、各センサＳａ〜Ｓｇからそれぞれ取得する。
【００５８】
次に、算出部３１６は、取得部３１４により取得された各種センサデータａ〜ｇの各々の値の特徴に基づき、生体中のグルコースＧｌｃの濃度を血糖値として算出する（ステップＳ１０８）。次に、評価部３１８は、予め取得された実測血糖値を示すデータに基づき、算出部３１６により算出された算出血糖値の評価を行う（ステップＳ１１０）。次に、評価部３１８は、評価した評価結果を表示部３０４に出力する（ステップＳ１１２）。これによって、血糖値測定装置１は、本フローチャートの処理を終了する。
【００５９】
以上説明した第１の実施形態に係る血糖値測定装置１によれば、生体の皮膚から排出される皮膚ガスを検出し、検出した皮膚ガスのガス濃度に基づき、生体中の血糖値を算出することにより、非侵襲で生体中の血糖値を精度良く測定することができる。
【００６０】
また、第１の実施形態に係る血糖値測定装置１は、非侵襲であり、且つ生体の排出物である皮膚ガスを検出対象としているため、利用者の負担が小さい。この結果、血糖値測定装置１は、利用者の血糖値を常時測定することができる。特に、血糖値測定装置１は、皮膚ガスが利用者の意思とは関係なく排出されることから、睡眠時においても利用者の血糖値を測定することができる。
【００６１】
＜第２の実施形態＞
以下、第２の実施形態に係る血糖値測定装置１について説明する。第２の実施形態では、第１の実施形態との相違点として、評価部３１８が各種センサデータａ〜ｇと算出血糖値とに基づき、利用者の健康状態を評価する場合について説明する。なお、上述した実施形態と共通する機能等についての説明は省略する。
【００６２】
図６は、利用者の健康状態を評価する際に参照される参照テーブルの一例を示す図である。
評価部３１８は、算出部３１６により算出されたグルコース濃度に関連する物質ごとの濃度および算出血糖値と、利用者が当該血糖値測定装置１を利用する前に食事したか否かの情報とに基づき、利用者の健康状態を評価する。食事したか否かの情報は、例えば、表示部３０４によって利用者からの操作入力として受け付け取得される。以下に、一例として、利用者の健康状態を評価する方法を説明する。
【００６３】
評価部３１８は、例えば、グルコース濃度に関連する物質のうち、水Ｈ２Ｏの濃度が閾値よりも低く、且つアセトンＲＣＯＲの濃度が閾値よりも高い場合、ミトコンドリア内のクエン酸回路や電子伝達系等の機能に異常をきたしているとして、利用者の健康状態を「生体機能異常」として評価する。また、評価部３１８は、例えば、利用者が食事済みの場合において、水Ｈ２Ｏの濃度が閾値よりも低く、且つアセトンＲＣＯＲの濃度が閾値よりも低い場合、上記クエン酸回路や電子伝達系等の機能の異常に加え、更にβ酸化の代謝経路にも異常をきたしているとして、利用者の健康状態を「糖尿病」として評価する。これら評価結果によって、利用者は、自身の健康状態を把握することができる。
【００６４】
また、評価部３１８は、利用者の皮膚ガスからアルコールＲＯＨが検出された場合、以下の処理を行ってもよい。図７は、アルコールＲＯＨが検出された利用者の健康状態を評価する際に参照される参照テーブルの一例を示す図である。
【００６５】
評価部３１８は、例えば、水Ｈ２Ｏの濃度が閾値よりも高く、アルコールＲＯＨの濃度が閾値よりも高い場合、算出血糖値が閾値よりも低い場合であっても、利用者の健康状態を「健康」として評価する。これは、飲酒等によって過剰にアルコールＲＯＨが摂取された場合、利用者の体内では、アルコールを分解するために、体内に存在するグルコースＧｌｃが多量に代謝され、その結果、一時的に水Ｈ２Ｏの濃度が閾値よりも高くなるためである。なお、アルコールＲＯＨの代わりに、アルコールＲＯＨの分解過程に生成されるアセトアルデヒドＲＣＨＯや酢酸ＲＣＯＯＨ等の物質を評価する際の指標としてもよい。これによって、利用者は、自身の健康状態を把握することができる。
【００６６】
なお、第２の実施形態において説明した図６および図７に示す参照テーブル図は、一例を示す図であり、医療従事者や医療関係者等による適切な指示（ガイドライン）に基づき適宜変更されてよい。
【００６７】
＜第３の実施形態＞
以下、第３の実施形態に係る血糖値測定装置１について説明する。第３の実施形態では、第１および第２の実施形態との相違点として、算出部３１６がグルコースＧｌｃの濃度を算出する際に用いる他の算出式について説明する。なお、上述した実施形態と共通する機能等についての説明は省略する。
【００６８】
上述した実施形態の場合、血糖値測定装置１により測定される血糖値は、利用者の年齢や活動量等によって、測定結果に誤差を含む場合がある。これは、利用者の年齢や活動量等に応じた皮膚の活性状態が異なるためである。皮膚の活性状態は、特に皮膚に含まれる水分量によって決められる。また、生体中のグルコース濃度が同じでも、皮膚の活性が高くなると皮膚ガスの排出量が増大する。
【００６９】
そのため、算出部３１６は、例えば、式（８）〜（１０）に基づき、グルコース濃度に関連する物質ごとの濃度を算出する。例えば、二酸化炭素ＣＯ２の濃度［ＣＯ２］は、式（８）により算出され、一酸化炭素ＣＯの濃度［ＣＯ］は、式（９）により算出され、アルコールＲＯＨの濃度［ＲＯＨ］は、式（１０）により、処理値として算出される。式（８）〜（１０）に示すａ〜ｇは、センサデータ（検出値）ａ〜ｇにそれぞれ対応する。なお、式（８）〜（１０）に示す各種センサデータａ〜ｇに乗算される係数（重み付け係数）Ａ１〜３〜Ｇ１〜３は、上述した決定方法と同様である。また、グルコース濃度に関連するその他物質の濃度算出式（［ＲＣＨＯ］、［ＲＣＯＯＨ］、［ＲＣＯＲ］、［Ｈ２Ｏ］等の処理値）については省略する。また、省略したその他物質の濃度算出式中の各係数（例えば、Ａ４〜Ｇ４、Ａ５〜Ｇ５、Ａ６〜Ｇ６、Ａ７〜Ｇ７のような重み付け係数）においても、上述した方法と同じように決定される。
【００７０】
［ＣＯ２］＝Ａ１・ａ/ｇ＋Ｂ１・ｂ/ｇ＋Ｃ１・ｃ/ｇ＋Ｄ１・ｄ/ｇ＋Ｅ１・ｅ/ｇ＋Ｆ１・ｆ/ｇ＋Ｇ１    （８）
【００７１】
［ＣＯ］＝Ａ２・ａ/ｇ＋Ｂ２・ｂ/ｇ＋Ｃ２・ｃ/ｇ＋Ｄ２・ｄ/ｇ＋Ｅ２・ｅ/ｇ＋Ｆ２・ｆ/ｇ＋Ｇ２    （９）
【００７２】
［ＲＯＨ］＝Ａ３・ａ/ｇ＋Ｂ３・ｂ/ｇ＋Ｃ３・ｃ/ｇ＋Ｄ３・ｄ/ｇ＋Ｅ３・ｅ/ｇ＋Ｆ３・ｆ/ｇ＋Ｇ３    （１０）
【００７３】
式（８）〜（１０）において、グルコース濃度に関連する物質ごとの濃度を算出する際、各種センサデータ（検出値）ａ〜ｇは、センサデータｇ（除算用検出値）によってそれぞれ除算されるような関数を用いた算術処理（重み付け計算）が行われる。センサデータｇは、皮膚ガスに含まれる水Ｈ２Ｏを最も良く検出したデータであることから、上記除算演算により、その他センサデータａ〜ｆに含まれる水Ｈ２Ｏの濃度の影響を補正することができる。このような関数を用いた算術処理（重み付け、重み付け計算）が行われる結果、血糖値測定装置１は、皮膚の活性の影響を抑制し、利用者の年齢や活動量に関わらず、血糖値を精度良く測定することができる。
【００７４】
また、算出部３１６は、例えば、式（１１）〜（１３）に基づき、グルコース濃度に関連する物質ごとの濃度を算出してもよい。例えば、二酸化炭素ＣＯ２の濃度［ＣＯ２］は、式（１１）により算出され、一酸化炭素ＣＯの濃度［ＣＯ］は、式（１２）により算出され、アルコールＲＯＨの濃度［ＲＯＨ］は、式（１３）により、処理値として算出される。式（１１）〜（１３）に示すａ〜ｇは、センサデータａ〜ｇにそれぞれ対応する。なお、式（１１）〜（１３）に示す各種センサデータａ〜ｇに乗算される係数（重み付け係数）Ａ１〜３〜Ｇ１〜３は、上述した決定方法と同様である。また、グルコース濃度に関連するその他物質の濃度算出式（［ＲＣＨＯ］、［ＲＣＯＯＨ］、［ＲＣＯＲ］、［Ｈ２Ｏ］等の処理値）については省略する。また、省略したその他物質の濃度算出式中の各係数（例えば、Ａ４〜Ｇ４、Ａ５〜Ｇ５、Ａ６〜Ｇ６、Ａ７〜Ｇ７のような重み付け係数）においても、上述した方法と同じように決定される。
【００７５】
［ＣＯ２］＝Ａ１・ａ２＋Ｂ１・ｂ２＋Ｃ１・ｃ２＋Ｄ１・ｄ２＋Ｅ１・ｅ２＋Ｆ１・ｆ２＋Ｇ１・ｇ２    （１１）
【００７６】
［ＣＯ］＝Ａ２・ａ２＋Ｂ２・ｂ２＋Ｃ２・ｃ２＋Ｄ２・ｄ２＋Ｅ２・ｅ２＋Ｆ２・ｆ２＋Ｇ２・ｇ２      （１２）
【００７７】
［ＲＯＨ］＝Ａ３・ａ２＋Ｂ３・ｂ２＋Ｃ３・ｃ２＋Ｄ３・ｄ２＋Ｅ３・ｅ２＋Ｆ３・ｆ２＋Ｇ３・ｇ２    （１３）
【００７８】
式（１１）〜（１３）において、グルコース濃度に関連する物質ごとの濃度を算出する際、各種センサデータａ〜ｇ（検出値）は２乗され、対応する係数（重み付け係数）Ａ１〜３〜Ｇ１〜３がそれぞれ乗算されるような関数を用いた算術処理（重み付け計算）が行われる。このような関数を用いた算術処理（重み付け、重み付け計算）を行うことによって、各種センサデータに示されるグルコース濃度に関連する物質ごとの濃度は、より顕著なデータとして扱われる。例えば、センサデータａは、皮膚ガスに含まれる二酸化炭素ＣＯ２の特徴を最も良く検出したデータであることから、センサデータａに含まれる二酸化炭素ＣＯ２の濃度は、その他物質の濃度や外乱値より突出した値になり、演算過程において、その他物質の濃度や外乱値の影響を抑制することができる。この結果、血糖値測定装置１は、血糖値を精度良く測定することができる。
【００７９】
また、算出部３１６は、例えば、式（１４）〜（１６）に基づき、グルコース濃度に関連する物質ごとの濃度を算出してもよい。例えば、二酸化炭素ＣＯ２の濃度［ＣＯ２］は、式（１４）により算出され、一酸化炭素ＣＯの濃度［ＣＯ］は、式（１５）により算出され、アルコールＲＯＨの濃度［ＲＯＨ］は、式（１６）により、処理値として算出される。式（１４）〜（１６）に示すａ〜ｇは、センサデータ（検出値）ａ〜ｇにそれぞれ対応する。なお、式（１４）〜（１６）に示す各種センサデータａ〜ｇに乗算される係数（重み付け係数）Ａ１〜３〜Ｇ１〜３は、上述した決定方法と同様である。また、グルコースＧｌｃの濃度に関連するその他物質の濃度算出式（［ＲＣＨＯ］、［ＲＣＯＯＨ］、［ＲＣＯＲ］、［Ｈ２Ｏ］等の処理値）については省略する。また、省略したその他物質の濃度算出式中の各係数（例えば、Ａ４〜Ｇ４、Ａ５〜Ｇ５、Ａ６〜Ｇ６、Ａ７〜Ｇ７のような重み付け係数）においても、上述した方法と同じように決定される。
【００８０】
［ＣＯ２］＝Ａ１・ｆ（ａ）＋Ｂ１・ｆ（ｂ）＋Ｃ１・ｆ（ｃ）＋Ｄ１・ｆ（ｄ）＋Ｅ１・ｆ（ｅ）＋Ｆ１・ｆ（ｆ）＋Ｇ１・ｆ（ｇ）    （１４）
【００８１】
［ＣＯ］＝Ａ２・ｆ（ａ）＋Ｂ２・ｆ（ｂ）＋Ｃ２・ｆ（ｃ）＋Ｄ２・ｆ（ｄ）＋Ｅ２・ｆ（ｅ）＋Ｆ２・ｆ（ｆ）＋Ｇ２・ｆ（ｇ）    （１５）
【００８２】
［ＲＯＨ］＝Ａ３・ｆ（ａ）＋Ｂ３・ｆ（ｂ）＋Ｃ３・ｆ（ｃ）＋Ｄ３・ｆ（ｄ）＋Ｅ３・ｆ（ｅ）＋Ｆ３・ｆ（ｆ）＋Ｇ３・ｆ（ｇ）    （１６）
【００８３】
式（１４）〜（１６）において、ｆ（ａ）〜ｆ（ｇ）は、各種センサデータａ〜ｇに関する関数（重み付け関数）である。関数ｆ（ａ）〜ｆ（ｇ）には、対応した各種センサデータａ〜ｇ以外に、測定中に混入することが想定される外乱等の誤差を考慮した数式項を含まれていてもよい。また、関数ｆ（ａ）〜ｆ（ｇ）は、指数関数や対数関数等の非線形関数で表されていてもよい。このような関数を用いた算術処理（重み付け、重み付け計算）を行うことによって、血糖値測定装置１は、グルコースＧｌｃの濃度の算出時に生じる誤差の影響を抑制することができ、血糖値を精度良く測定することができる。
【００８４】
＜その他の実施形態＞
以下に、その他の実施形態（変形例）について記載する。
皮膚ガスプローブ１００は、指から皮膚ガスを採取するものに限定されず、例えば、手のひら、腕、腹等の血管が集中する皮膚から皮膚ガスを採取するプローブであってもよい。この場合、皮膚ガスプローブ１００は、効率良く皮膚ガスを採取するため、皮膚表面に張り付けるパッド式のプローブであると好適である。また、皮膚ガスプローブ１００は、皮膚ガスを採取するための十分な採取面積があれば、どのような形状でもよく、例えば、手首や腕等に巻きつけるバンド型であってもよい。また、血糖値測定装置１は、皮膚ガスプローブ１００、皮膚ガス検出ユニット２００、および制御装置３００が一体となった装置であってもよく、利用者が携行可能な所謂ウェアラブル端末装置であってもよい。
【００８５】
以上、実施形態および実施例を用いて本発明の血糖値測定装置について説明したが、本発明の技術範囲は上記実施形態および実施例に限定されるものではなく、本発明の趣旨を逸脱しない範囲において構成要素の組み合わせを変えたり、各構成要素に種々の変更を加えたり、削除したりすることが可能である。
【符号の説明】
【００８６】
１‥血糖値測定装置
１００…皮膚ガスプローブ（採取部）
２００…皮膚ガス検出ユニット
２０２…第１ポンプ
２０４…排気口
２０６…吸気口
２０８…乾燥器
２１０…第２ポンプ
Ｖ１〜Ｖ１４…電磁バルブ（搬送部）
Ｓａ〜Ｓｇ…センサ（検出部）
３００…制御装置
３０２…通信部
３０４…表示部
３１０…制御部
３１２…ガス供給制御部（搬送制御部）
３１４…取得部
３１６…算出部
３１８…評価部
３２０…記憶部
【図１】
2016072513.tif 000003
【図２】
2016072513.tif 000004
【図３】
2016072513.tif 000005
【図４Ａ】
2016072513.tif 000006
【図４Ｂ】
2016072513.tif 000007
【図５】
2016072513.tif 000008
【図６】
2016072513.tif 000009
【図７】
2016072513.tif 000010
【国際調査報告】
2016072513.tif 000011
2016072513.tif 000012
2016072513.tif 000013
2016072513.tif 000014