JP2021012078A - 温度分布推定装置及び温度分布推定方法 - Google Patents
温度分布推定装置及び温度分布推定方法 Download PDFInfo
- Publication number
- JP2021012078A JP2021012078A JP2019125910A JP2019125910A JP2021012078A JP 2021012078 A JP2021012078 A JP 2021012078A JP 2019125910 A JP2019125910 A JP 2019125910A JP 2019125910 A JP2019125910 A JP 2019125910A JP 2021012078 A JP2021012078 A JP 2021012078A
- Authority
- JP
- Japan
- Prior art keywords
- temperature distribution
- temperature
- semivariogram
- distribution estimation
- sample
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims description 50
- 238000005259 measurement Methods 0.000 claims abstract description 85
- 238000004364 calculation method Methods 0.000 claims description 14
- 238000012800 visualization Methods 0.000 claims description 12
- 239000011159 matrix material Substances 0.000 claims description 8
- 238000002360 preparation method Methods 0.000 claims description 5
- 239000006185 dispersion Substances 0.000 claims 2
- 238000004458 analytical method Methods 0.000 description 18
- 230000006870 function Effects 0.000 description 12
- 238000012795 verification Methods 0.000 description 11
- 238000010586 diagram Methods 0.000 description 10
- 239000012530 fluid Substances 0.000 description 10
- 238000000611 regression analysis Methods 0.000 description 5
- 238000007796 conventional method Methods 0.000 description 4
- 238000011156 evaluation Methods 0.000 description 4
- 238000012545 processing Methods 0.000 description 4
- 238000004891 communication Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 2
- 238000010438 heat treatment Methods 0.000 description 2
- 230000010365 information processing Effects 0.000 description 2
- 238000007476 Maximum Likelihood Methods 0.000 description 1
- ATJFFYVFTNAWJD-UHFFFAOYSA-N Tin Chemical compound [Sn] ATJFFYVFTNAWJD-UHFFFAOYSA-N 0.000 description 1
- 238000004378 air conditioning Methods 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 230000006866 deterioration Effects 0.000 description 1
- 238000005290 field theory Methods 0.000 description 1
- 230000020169 heat generation Effects 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 230000002194 synthesizing effect Effects 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Landscapes
- Measuring Temperature Or Quantity Of Heat (AREA)
- Air Conditioning Control Device (AREA)
Abstract
Description
所定空間の複数の計測位置でそれぞれ計測された複数の計測温度を含む温度データ群を取得する取得部と、
前記温度データ群に基づいて、前記所定空間の鉛直方向及び水平方向の間に異方性を有する空間温度特性を推定する空間温度特性推定部と、
前記複数の計測位置から前記所定空間の任意の推定位置までの距離と、前記空間温度特性とに応じた重み係数を用いて前記複数の計測温度に対する空間補間を行って、前記推定位置の温度を算定することにより、前記所定空間の温度分布を推定する温度分布推定部とを備え、
前記温度分布推定部は、
前記推定位置が前記計測位置に対して前記鉛直方向に位置するほど当該計測位置の前記重み係数を小さく、前記推定位置が前記計測位置に対して前記水平方向に位置するほど当該計測位置の前記重み係数を大きくして、前記空間補間を行う。
所定空間の複数の計測位置でそれぞれ計測された複数の計測温度を含む温度データ群を取得する取得工程と、
前記温度データ群に基づいて、前記所定空間における鉛直方向及び水平方向の空間温度特性を推定する空間温度特性推定工程と、
前記複数の計測位置から前記所定空間の任意の推定位置までの距離と、前記空間温度特性とに応じた重み係数を用いて前記複数の計測温度に対する空間補間を行って、前記推定位置の温度を算定することにより、前記所定空間の温度分布を推定する温度分布推定工程とを備え、
前記温度分布推定工程は、
前記推定位置が前記計測位置に対して前記鉛直方向に位置するほど当該計測位置の前記重み係数を小さく、前記推定位置が前記計測位置に対して前記水平方向に位置するほど当該計測位置の前記重み係数を大きくして、前記空間補間を行う。
図2は、本発明の実施形態に係る温度分布推定装置3を示すブロック図である。温度分布推定装置3は、例えば、汎用のコンピュータ等の情報処理装置で構成されている。すなわち、温度分布推定装置3は、キーボード、タッチパネル等により構成される入力部30と、HDD、メモリ等により構成される記憶部31と、CPU等のプロセッサにより構成される制御部32と、ネットワーク5との通信インターフェースである通信部33と、ディスプレイ、スピーカ等により構成される出力部34とを備える。なお、入力部30や出力部34は省略してもよく、その場合には、管理者端末4を、入力部30や出力部34として機能させればよい。
所定空間における任意のパラメータ分布を把握したい場合、当該パラメータを計測可能な複数のセンサを所定空間内に配置し、当該複数のセンサにより計測された離散的な計測値の間を内挿(点と点の間の連続値を予測)することにより、所定空間のパラメータ分布を推定する方法がある。このようなパラメータ分布を推定する方法として、様々なノンパラメトリック手法(例えば、多項式補間、逆距離荷重補間(IDW)、スプライン法、Natural Neighbor法、TIN法等)が存在するが、本実施形態では、その1つであるクリギング(Kriging)を用いる。以下に、一般的なクリギングによる空間補間手法について説明する。
次に、本実施形態に係る温度分布推定装置3による温度分布推定処理の推定性能を検証するため、流体解析による解析結果(検証用温度分布データ)と、温度分布推定結果との比較を行った。以下、流体解析の解析条件と、温度分布推定結果の検証条件及び検証結果について説明する。
図8は、流体解析の解析条件を示し、(a)は解析モデル、(b)は境界条件を示す図である。
図9は、検証用温度分布データと、温度分布推定結果とを比較するときの各評価位置を示し、(a)は平面図、(b)はA−A線断面図である。図10は、温度データ群の各計測位置と、温度分布推定処理の各ケースを示す図である。
図11は、検証用温度分布データ、温度分布推定結果及び平均絶対誤差を示す図である。
以上、本発明の実施形態について説明したが、本発明は上記実施形態に限定されるものではなく、本発明の技術的思想を逸脱しない範囲で適宜変更可能である。
4…管理者端末、5…ネットワーク、10…空間、
11V…鉛直方向の標本範囲、11H…水平方向の標本範囲、
30…入力部、31…記憶部、32…制御部、33…通信部、34…出力部、
310…温度分布推定プログラム、
320…取得部、321…空間温度特性推定部321、
321a…標本セミバリオグラム作成部、321b…異方性パラメータ算定部、
321c…理論セミバリオグラム推定部、
322…温度分布推定部、323…可視化情報出力部
Claims (8)
- 所定空間の複数の計測位置でそれぞれ計測された複数の計測温度を含む温度データ群を取得する取得部と、
前記温度データ群に基づいて、前記所定空間の鉛直方向及び水平方向の間に異方性を有する空間温度特性を推定する空間温度特性推定部と、
前記複数の計測位置から前記所定空間の任意の推定位置までの距離と、前記空間温度特性とに応じた重み係数を用いて前記複数の計測温度に対する空間補間を行って、前記推定位置の温度を算定することにより、前記所定空間の温度分布を推定する温度分布推定部とを備え、
前記温度分布推定部は、
前記推定位置が前記計測位置に対して前記鉛直方向に位置するほど当該計測位置の前記重み係数を小さく、前記推定位置が前記計測位置に対して前記水平方向に位置するほど当該計測位置の前記重み係数を大きくして、前記空間補間を行う、
ことを特徴とする温度分布推定装置。 - 前記空間温度特性推定部は、
前記温度データ群に基づいて、鉛直方向の標本セミバリオグラム及び水平方向の標本セミバリオグラムを作成する標本セミバリオグラム作成部と、
前記鉛直方向の標本セミバリオグラム及び前記水平方向の標本セミバリオグラムに基づいて、前記異方性を表す異方性パラメータを算定する異方性パラメータ算定部と、
前記鉛直方向の標本セミバリオグラム及び前記水平方向の標本セミバリオグラムと、前記異方性パラメータとに基づいて、前記空間温度特性として理論セミバリオグラムを推定する理論セミバリオグラム推定部とを備え、
前記温度分布推定部は、
前記理論セミバリオグラムに基づいて、前記空間補間を行う
ことを特徴とする請求項1に記載の温度分布推定装置。 - 前記異方性パラメータは、
前記理論セミバリオグラムが所定の値となる点の集合を楕円として表すためのパラメータであって、
前記楕円の長軸を前記水平方向に一致させる回転角度と、
前記楕円の長軸に対する前記楕円の短軸の比率である異方性比とを含む、
ことを特徴とする請求項2に記載の温度分布推定装置。 - 前記異方性パラメータ算定部は、
前記複数の計測温度に基づいて分散値を算出し、前記分散値でバリオグラムモデルのシルを固定した状態で、前記鉛直方向の標本セミバリオグラムに前記バリオグラムモデルをフィッティングさせたときの前記鉛直方向のレンジと、前記水平方向の標本セミバリオグラムに前記バリオグラムモデルをフィッティングさせたときの前記水平方向のレンジとをそれぞれ算定し、
前記水平方向のレンジに対する前記鉛直方向のレンジの比率に基づいて、前記異方性比を算定する、
ことを特徴とする請求項3に記載の温度分布推定装置。 - 前記理論セミバリオグラム推定部は、
前記異方性比を加味して、前記鉛直方向の標本セミバリオグラム及び前記水平方向の標本セミバリオグラムにバリオグラムモデルをフィッティングさせることで、等方的な理論バリオグラムを推定し、
前記等方的な理論バリオグラムと、1及び前記異方性比を要素とする対角行列とに基づいて、前記理論セミバリオグラムを推定する、
ことを特徴とする請求項4に記載の温度分布推定装置。 - 前記標本セミバリオグラム作成部は、
前記複数の計測位置をそれぞれ基準位置とするとき、当該基準位置からみて前記鉛直方向に対して所定の角度の範囲内に含まれる前記計測位置に基づいて、前記鉛直方向の標本セミバリオグラムを作成するとともに、当該基準位置からみて前記水平方向に対して所定の角度の範囲内に含まれる前記計測位置に基づいて、前記水平方向の標本セミバリオグラムを作成する、
ことを特徴とする請求項2乃至請求項5のいずれか一項に記載の温度分布推定装置。 - 前記温度分布推定部により推定された前記温度分布を可視化するための可視化情報を出力する可視化情報出力部をさらに備える、
ことを特徴とする請求項1乃至請求項6のいずれか一項に記載の温度分布推定装置。 - 所定空間の複数の計測位置でそれぞれ計測された複数の計測温度を含む温度データ群を取得する取得工程と、
前記温度データ群に基づいて、前記所定空間における鉛直方向及び水平方向の空間温度特性を推定する空間温度特性推定工程と、
前記複数の計測位置から前記所定空間の任意の推定位置までの距離と、前記空間温度特性とに応じた重み係数を用いて前記複数の計測温度に対する空間補間を行って、前記推定位置の温度を算定することにより、前記所定空間の温度分布を推定する温度分布推定工程とを備え、
前記温度分布推定工程は、
前記推定位置が前記計測位置に対して前記鉛直方向に位置するほど当該計測位置の前記重み係数を小さく、前記推定位置が前記計測位置に対して前記水平方向に位置するほど当該計測位置の前記重み係数を大きくして、前記空間補間を行う、
ことを特徴とする温度分布推定方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2019125910A JP7316855B2 (ja) | 2019-07-05 | 2019-07-05 | 温度分布推定装置及び温度分布推定方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2019125910A JP7316855B2 (ja) | 2019-07-05 | 2019-07-05 | 温度分布推定装置及び温度分布推定方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2021012078A true JP2021012078A (ja) | 2021-02-04 |
JP7316855B2 JP7316855B2 (ja) | 2023-07-28 |
Family
ID=74227483
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2019125910A Active JP7316855B2 (ja) | 2019-07-05 | 2019-07-05 | 温度分布推定装置及び温度分布推定方法 |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP7316855B2 (ja) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113884219A (zh) * | 2021-12-06 | 2022-01-04 | 中国空气动力研究与发展中心低速空气动力研究所 | 一种表面温度的测量装置和测量方法 |
WO2023199869A1 (ja) * | 2022-04-12 | 2023-10-19 | パナソニックIpマネジメント株式会社 | 推定方法、プログラム及び推定システム |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH04205570A (ja) * | 1990-11-30 | 1992-07-27 | Hitachi Plant Eng & Constr Co Ltd | 三次元温度分布表示方法 |
JP2008107963A (ja) * | 2006-10-24 | 2008-05-08 | Hitachi Ltd | 3次元データ推定方法およびプログラム |
JP2008241433A (ja) * | 2007-03-27 | 2008-10-09 | Nec Corp | 観測データ同化方法 |
-
2019
- 2019-07-05 JP JP2019125910A patent/JP7316855B2/ja active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH04205570A (ja) * | 1990-11-30 | 1992-07-27 | Hitachi Plant Eng & Constr Co Ltd | 三次元温度分布表示方法 |
JP2008107963A (ja) * | 2006-10-24 | 2008-05-08 | Hitachi Ltd | 3次元データ推定方法およびプログラム |
JP2008241433A (ja) * | 2007-03-27 | 2008-10-09 | Nec Corp | 観測データ同化方法 |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113884219A (zh) * | 2021-12-06 | 2022-01-04 | 中国空气动力研究与发展中心低速空气动力研究所 | 一种表面温度的测量装置和测量方法 |
WO2023199869A1 (ja) * | 2022-04-12 | 2023-10-19 | パナソニックIpマネジメント株式会社 | 推定方法、プログラム及び推定システム |
Also Published As
Publication number | Publication date |
---|---|
JP7316855B2 (ja) | 2023-07-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP2010078447A (ja) | Q値解析システム、q値解析方法及びq値解析プログラム | |
KR101912195B1 (ko) | 목표의 실시간 위치 측정과 지도 구축 방법 및 장치 | |
Wunsch et al. | Practical global oceanic state estimation | |
WO2019181032A1 (en) | Wind flow sensing system and method for determining velocity fields of wind flow | |
JP6001740B1 (ja) | 構造物の伝達関数の高精度評価・地震応答予測・劣化診断システム及びその方法 | |
Afrasiabi et al. | Thermal simulation in multiphase incompressible flows using coupled meshfree and particle level set methods | |
JP7316855B2 (ja) | 温度分布推定装置及び温度分布推定方法 | |
JP6159594B2 (ja) | 観測値処理装置、観測値処理方法及び観測値処理プログラム | |
JP4866896B2 (ja) | 気象予測装置、気象予測方法および気象予測プログラム | |
van Ophem et al. | Model based virtual intensity measurements for exterior vibro-acoustic radiation | |
Fan et al. | Propagation of vertical and horizontal source data errors into a TIN with linear interpolation | |
Lee et al. | System identification of spatial distribution of structural parameters using modified transitional Markov Chain Monte Carlo method | |
Ubertini et al. | Computer simulation of stochastic wind velocity fields for structural response analysis: comparisons and applications | |
Albuquerque et al. | Residual least-squares error estimate for unstructured h-adaptive meshes | |
JP2018194497A (ja) | 気象予測システム、気象予測方法、気象予測プログラム | |
EP4127457A1 (en) | System and method for wind flow turbulence measurement by lidar in a complex terrain | |
JP2018136671A (ja) | 内部波パラメータ推定装置および内部波パラメータ推定方法 | |
JP5981901B2 (ja) | 気流推定方法およびその装置 | |
JP2021096084A (ja) | 測定機器 | |
Yu et al. | Defect identification in heat transfer problems using boundary data | |
Gupta et al. | A reduced modal parameter based algorithm to estimate excitation forces from optimally placed accelerometers | |
WO2022004777A1 (ja) | 温度分布学習装置 | |
JP6097200B2 (ja) | データ推定方法およびその装置 | |
JP5562895B2 (ja) | 日射量推定装置、日射量推定方法及びプログラム | |
JP5881629B2 (ja) | 情報処理装置及び発生源推定方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20220627 |
|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20230315 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20230412 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20230530 |
|
TRDD | Decision of grant or rejection written | ||
A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20230705 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20230718 |
|
R150 | Certificate of patent or registration of utility model |
Ref document number: 7316855 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |