[論文レビュー] Optimization of Closed-Loop Shallow Geothermal Systems Using Analytical Models

O. Heinzel, Smajil Halilovic+2 more|arXiv (Cornell University)|Mar 26, 2026

Geothermal Energy Systems and Applications被引用数 0

ひとこと要約

この論文は、Finite Line Source 土壌モデルと 1U ボアホール熱伝達モデルを組み合わせ、Modified Lloyd’s アルゴリズムを用いて非凸領域におけるBHEを最適配置し、その後流体温度制約の下で均一ボアホール長を最小化する高速解析フレームワークを提示します。

ABSTRACT

Closed-loop shallow geothermal systems are one of the key technologies for decarbonizing the residential heating and cooling sector. The primary type of these systems involves vertical borehole heat exchangers (BHEs). During the planning phase, it is essential to find the optimal design for these systems, including the depth and spatial arrangement of the BHEs. In this work, we have developed a novel approach to find the optimal design of BHE fields, taking into account constraints such as temperature limits of the heat carrier fluid. These limits correspond to the regulatory practices applied during the planning phase. The approach uses a finite line source model to simulate temperature changes in the ground in combination with an analytical model of heat transport within the boreholes. Our approach is demonstrated using realistic scenarios and is expected to improve current practice in the planning and design of BHE systems.

研究の動機と目的

• 1U ボアホールの20年間の性能を予測するための迅速な解析モデルを開発する。
• 非凸領域で最適なBHE座標を自動的に位置付ける堅牢なヒューリスティックを作成する。
• シミュレーション期間中に流体温度制限を守りつつ総BHE長を最小化する深度最適化と空間配置を統合する。

提案手法

• Method of Images を用いた有限線源とGreenの関数で土壌温度をモデル化する。
• 境界条件と解析的出口温度 Tout を用いた準定常の1U伝達でボアホール流体結合をモデル化する。
• 自己およびボアホール間項とFFTによって加速されるデュアルグリッド畳み込みで土壌-流体相互作用を O(Nt log Nt) で計算する。
• BHEフィールド設計を2段階の最適化として定式化する: (i) 非凸領域で Modified Lloyd’s CVT による inter-BHE および境界間距離の最大化、(ii) Tout を Tmin/Tmax の制約の下で時間を通じて満たすように均一ボアホール長 L を最小化。
• 空間分布を非凸 CVT で最適化するための Modified Lloyd’s アルゴリズムにより中心点を領域内へ再投影。
• Tout を境界内に保ちながら L を最小化するための Sequential Least Squares Programming (SLSQP) bound + nonlinear constrained 問題で深度最適化を解く。

実験結果

リサーチクエスチョン

• RQ1 BHEフィールドを非凸領域でどのように最適配置して距離を最大化し干渉を最小化できるか？
• RQ2 与えられた熱需要プロファイルで20年間の流体温度制約を満たす最小の均一ボアホール長 L はどれか？
• RQ3 提案されたCVTベースの配置はグリッド、長方形領域、L字型領域など、異なる特性ジオメトリでどのように性能を発揮するか？

主な発見

Geometry	Optimal Length (m)
Grid (5x5)	91.47
Rect. (Medium)	91.44
Rect. (Small)	101.86
Rect. (Large)	84.07
L-Shape	80.72

• Modified Lloyd’s CVT 法は非 convex領域におけるBHE配置を近似最適グリッドに近づけ、最近傍距離の平均差は対称ケースで約2 m。
• 最適ボアホール長は特性ジオメトリによって異なる：Grid (5x5) 91.47 m; Rectangular Medium 91.44 m; Rectangular Small 101.86 m; Rectangular Large 84.07 m; L-shaped domain 80.72 m。
• より広い利用可能エリアはボアホール間隔を広げ、干渉を低減するため L を小さくできる。一方で小さなエリアは温度制約を満たすためにボアを深くする必要がある。
• このフレームワークでは20年の期間にわたり Tout を Tmin/Tmax 内に保ち、Tmax が冷却優勢プロファイルで制約要因となることが示され、領域認識深度設計の必要性を示唆する。