狭窄流の圧力損失推定を数値解析で比べる — どの式が真値に近いか
4D Flow MRI-Based Pressure Loss Estimation in Stenotic Flows: Evaluation Using Numerical Simulations
4D Flow から圧力損失を出す方法は複数ある。どれが真の非可逆圧較差に近いのか。
真値 6.76±0.54 mmHg に対し、簡易ベルヌーイは 8.74±0.67 mmHg と過大、拡張ベルヌーイは 6.57±0.53 mmHg で最も近い。圧ポアソン方程式はほぼゼロになる。TKE は分解能の影響をほとんど受けないが、粘性散逸は過小評価され分解能に依存する。
当社の視点
圧力損失の指標を選ぶときに、最初に見るべき比較表です。CFD で真値が分かる条件を作り、そこに MRI 計測を模擬して各手法を当てているため、手法そのものの優劣を切り分けられます。実務で効くのは2点。第一に、圧ポアソン方程式は非可逆圧較差をほぼ拾えないこと。第二に、TKE と粘性散逸は形状ごとの相関は強い（r²>0.93）のに、形状をまたぐと落ちる（r²=0.756 と 0.776）ことです。つまり同一症例の経時変化には使えても、症例間の絶対比較には向きません。
扱う範囲
-
01
検証の設計
狭窄度・狭窄後径・流量を変えた形状群について CFD で真値を作り、そこに空間分解能を変えた MRI 計測を模擬した。簡易ベルヌーイ、拡張ベルヌーイ、圧ポアソン方程式、TKE の体積積分、粘性散逸を、真の非可逆圧較差と突き合わせている。
-
02
ベルヌーイ2種
簡易ベルヌーイは真値 6.76±0.54 mmHg に対し 8.74±0.67 mmHg と過大に出る。狭窄後の圧力回復を無視するため。拡張ベルヌーイは 6.57±0.53 mmHg で最も近いが、真値が小さい領域では誤差が残る。
-
03
圧ポアソン方程式
この方程式による推定は常にほぼゼロになった。相対圧の分布は出せても、非可逆な圧力損失そのものは拾えていない。
-
04
TKE と粘性散逸
どちらも形状ごとに見れば真値と強く相関する（r²>0.93）。ただし全形状をまとめると TKE で r²=0.756、粘性散逸で r²=0.776 まで落ちる。TKE の推定は正確で分解能の影響が小さく、粘性散逸は全体に過小評価され分解能に依存した。
真の非可逆圧較差との比較
数値で確かめる
| 項目 | 値 | 補足 |
|---|---|---|
| 真値（TPGnet） | 6.76 ± 0.54 mmHg | CFD による基準値 |
| 簡易ベルヌーイ | 8.74 ± 0.67 mmHg | 圧力回復を無視するため過大 |
| 拡張ベルヌーイ | 6.57 ± 0.53 mmHg | 最も近い。ただし小さい真値で誤差が残る |
| 圧ポアソン方程式 | ほぼゼロ | 非可逆圧較差を拾えない |
| TKE との相関（形状ごと） | r² > 0.93 | 粘性散逸も同程度 |
| TKE との相関（全形状） | r² = 0.756 | 粘性散逸は r² = 0.776 |
| 分解能の影響 | TKE は小、粘性散逸は大 | 粘性散逸は全体に過小評価される |
いずれもこの論文が報告した値です。ノイズを含まない理想的な軸対称形状での結果。
実務チェック
- 圧較差を報告するとき、どの式で出したかを必ず書く。同じデータで 8.74 と 6.57 に分かれる。
- 非可逆圧較差を狙うなら圧ポアソン方程式を使わない。ほぼゼロになる。
- TKE や粘性散逸を症例間で比べない。形状をまたぐと相関が落ちる。
- 粘性散逸を使うなら空間分解能を明記する。値が分解能に依存する。
この研究の限界
- 理想化された軸対称形状での検証で、ノイズも含んでいません。
- 実際の患者形状や生体内での検証は行われていません。
- 参照となる圧測定は CFD の計算値であり、侵襲的な実測ではありません。
書誌情報
| 原題 | 4D Flow MRI-Based Pressure Loss Estimation in Stenotic Flows: Evaluation Using Numerical Simulations |
|---|---|
| 著者 | Casas B, Lantz J, Dyverfeldt P, Ebbers T |
| 掲載誌 | Magnetic Resonance in Medicine（2016） |
| DOI | 10.1002/mrm.25772 |
| PMID | 26016805 |