TKE/MKEを、4D Flow MRIのエネルギー読影として読む

まず一言で

結論: MKEは「主流がどれだけ勢いよく走るか」、TKEは「その流れがどれだけ崩れているか」を見る。

4D Flow MRIでは、速度ベクトルそのものだけでなく、速度から派生するエネルギー指標を読めます。代表が MKE と TKE です。MKEは秩序だった平均速度から計算され、TKEはボクセル内の速度ばらつきから計算されます。

つまり、同じ「高い値」でも意味が違います。MKEが高いなら主流の勢いが強い。TKEが高いなら、狭窄後、弁周囲、分岐、逆流、混合で流れが崩れている可能性があります。ただしTKE hot spotは、低SNR、部分容積、VENC不適合、ROI境界でも生じ得ます。

この記事の立ち位置

公開文献にもとづく教育用記事です。特定施設、個別案件、非公開実装情報、個人を特定し得る情報は含めていません。

式は暗記ではなく読影の地図

まず式を出します。ここで大事なのは計算そのものより、どの項が画像上のどの現象を表すかです。

MKE = 1/2 ρ |v|²
TKE ≈ 1/2 ρ (σ_x² + σ_y² + σ_z²)

ρ
血液密度。エネルギー密度へ換算する係数。文献ではおおむね 1060 kg/m³ 程度として扱われる。

v
平均速度ベクトル。そろって進む主流の速さと向きを表す。

|v|²
速度の二乗。速さが2倍になるとMKEは約4倍になる。

σ_x,y,z
3方向の速度標準偏差。ボクセル内で速度がどれだけばらつくかを表す。

MKEとTKEの数式項の意味 — 図2. MKEは平均速度 `v` の二乗、TKEは速度ばらつき `σ` の二乗から読む。式は「何を高値として見ているか」を分解する地図である。

さらに深く: なぜ二乗が出てくるのか

古典力学の運動エネルギーは 1/2 m v² です。4D Flow MRIでは、個々の粒子の質量ではなく、ボクセルやROI内の単位体積あたりで考えるため、質量密度 ρ が出てきます。

MKEは平均速度そのものの二乗です。一方TKEは、平均から外れた速度成分、つまり揺らぎ v' の二乗平均に対応します。直感的には、同じ平均速度でも、ボクセル内で速度がそろっていればTKEは低く、速度方向や大きさがばらつけばTKEが上がります。

Reynolds分解で見る: 平均とばらつきを分ける

TKEの考え方は、流体力学の Reynolds decomposition に近い読み方です。瞬間速度を平均成分と揺らぎ成分へ分けます。

v(t) = \bar{v} + v'(t)
k = 1/2 ρ \langle v'_x{}^2 + v'_y{}^2 + v'_z{}^2 \rangle

\bar{v}
時間またはアンサンブルで見た平均速度。MKEの主役。

v'(t)
平均からのずれ。TKEの主役。

k
乱流運動エネルギー密度。4D FlowではTKE mapとして扱われる。

読影の意味
平均流が強いのか、流れが崩れているのかを分ける。

数式が苦手なら、次のように読むと十分です。MKEは矢印がそろって太く走るイメージ、TKEは同じ場所で矢印が散らばるイメージです。

4D Flow MRIではσをどう測るのか

TKEは、単に速度mapを見て「乱れていそう」と言っているわけではありません。phase-contrast MRIでは、ボクセル内に異なる速度成分が混在すると、velocity encodingに応じて信号が減衰します。この減衰から、ボクセル内速度標準偏差 σ を推定する考え方が使われます。

S(k_v) ≈ S₀ exp(i k_v v) exp(-k_v² σ²/2)

exp(i k_v v)
平均速度 v による位相回転。MKEへつながる。

exp(-k_v²σ²/2)
速度ばらつきによる信号減衰。TKEへつながる。

Multi-VENC
複数の感度で信号を観測し、速度とばらつきを同時に推定しやすくする。

注意
低SNRや部分容積でも減衰・ばらつき様の所見が出るため、QCが必要。

ここで前の記事の Bayesian Unfolding と接続します。平均速度 v と速度分散 σ を同じ信号モデルから読むことで、MKEとTKEを同じ4D Flow datasetから考えられます。

J/m³ と mJ を混同しない

TKE/MKEで最も混乱しやすいのが単位です。ボクセルごとのエネルギー密度は J/m³、ROI内で体積積分した総量は mJ で報告されることがあります。

TKE/MKEのJ/m3とmJの違い — 図3. J/m³はボクセルごとのエネルギー密度、mJはROI全体を足し合わせた量。単位、ROI、空間分解能、時間フレームをそろえずに比較してはいけない。

報告形式	何を表すか	比較時の注意
J/m³	ボクセル単位のエネルギー密度	空間分解能、SNR、平滑化、VENC設定の影響を受ける。
mJ	ROI内を体積積分した総エネルギー	ROIサイズ、血管範囲、時間積分の定義で大きく変わる。
peak値	局所最大値	ノイズ、部分容積、境界設定に敏感。
平均値	ROI内の代表値	局所hot spotを薄める可能性がある。

臨床読影: 4象限で考える

MKE/TKEは、絶対値だけでなく組み合わせで読むと臨床的な意味が見えます。高速で整った流れ、高速で崩れた流れ、低速で混合する流れ、静かな低流量は、同じ「流れの異常」でも意味が違います。

MKE/TKEの4象限読影 — 図4. MKE/TKEの4象限。高MKE高TKEは高速jet後の乱れ、低MKE高TKEは低速混合や複雑分岐を疑う。hot spotは必ずQCと合わせる。

高MKE・低TKE

速いが整った流れ。弁通過後の整流域や正常大動脈の主流など、主流の勢いが主体。

高MKE・高TKE

高速jet後の乱れ。狭窄後流、弁周囲逆流、壁衝突後などで、エネルギー損失や組織負荷の議論へつながる。

低MKE・高TKE

平均流は弱いが、混合・逆流・分岐で速度がばらつく状態。低SNR由来の偽高値も疑う。

低MKE・低TKE

静かな低流量。循環自体が弱いのか、ROIが対象外なのか、時間フレームがずれていないかを確認する。

TKE hot spotをすぐ病態と決めない

TKEは「乱れ」を見たい指標ですが、測定・解析の要因でも上がります。特に局所hot spotは、病態のサインである前に、品質確認のトリガーとして扱うべきです。

TKE/MKEの失敗モードとQC — 図5. TKE hot spotの主な偽高値要因。低SNR、部分容積、VENC設定、ROI境界を確認してから臨床解釈へ進む。

アンチパターン

「TKEが高い = 乱流 = 重症」と短絡しない。速度画像、magnitude、位相折り返り、近接断面、ROI境界、時間フレームを見て、再現性を確認する。

Energy Loss、WSS、OSIとどう違うか

TKE/MKEは流れのエネルギーを読む入口です。ただし、Energy Loss、WSS、OSI、helicityとは同じものではありません。

指標	中心に見るもの	TKE/MKEとの関係
MKE	平均流の勢い	jetや高流量を読みやすいが、乱れそのものではない。
TKE	速度ばらつき・乱れ	狭窄後や混合の指標になり得るが、QC依存。
Energy Loss	散逸・圧損に近い概念	TKEと関連するが同義ではない。境界条件と計算モデルが重要。
WSS	壁面のせん断応力	壁近傍速度勾配に依存し、空間分解能とsegmentationに敏感。
OSI/RRT	せん断の方向変動・滞留	低WSSだけでは見えない時間方向の情報を足す。
Helicity/Vorticity	回転・渦・螺旋流	エネルギー量ではなく、流れの組織化を読む補助線。

さらに深く: TKEがEnergy Lossに見える理由と、同じでない理由

乱れが強い領域では、平均流のエネルギーが小さな速度揺らぎへ移り、最終的には粘性散逸へ向かいます。そのためTKE高値はエネルギー損失を連想させます。

ただしTKEはボクセル内速度ばらつきのエネルギー密度であり、圧損や散逸率そのものではありません。Energy Lossを議論するには、流量、圧力、粘性、境界条件、時間積分の前提が必要です。ここを混同すると、TKE mapだけで治療効果を断言する危険があります。

文献紹介: 式を支える論文、読影を支える図解

図6. TKE/MKEの文献地図。Dyverfeldtらは速度ばらつきのMRI計測、BinterらはMulti-VENCで速度と乱れを同時推定する流れを支える。

Dyverfeldtら 2006

phase-contrast MRIを拡張し、ボクセル内速度標準偏差と乱流強度を定量する基礎文献です。TKEを「見た目の乱れ」ではなく、MRI信号から推定する量として読む起点になります。

Binterら 2013

Bayesian multipoint velocity encodingにより、平均速度と乱流関連量を同時に扱う枠組みを示した重要文献です。前記事のBayesian Unfoldingと本記事のTKE/MKEをつなぎます。

4D Flow CMR consensus / phase-contrast CMR review

4D Flowの撮像、VENC、速度-to-noise、可視化、解析、報告の基本を整理する入口です。TKE/MKEを読む前に、速度map自体の品質と限界を押さえる必要があります。

Claim / Evidence / Limitation

Claim

MKEとTKEは、平均流と速度ばらつきを分けることで、4D Flow MRIのエネルギー読影を深くする。

Evidence

DyverfeldtらのTKE計測、BinterらのBayesian multipoint velocity encoding、PC-CMR/4D Flow consensus文献。

Limitation

TKE/MKEは撮像条件、VENC、SNR、部分容積、ROI、時間分解能に依存する。単独指標で病態を断言しない。

Zettelkastenへ戻す問い

MKEとTKEを、それぞれ「平均」と「ばらつき」の式から説明できるか。
J/m³とmJを混同せず、ROIとボクセルの違いを説明できるか。
TKE hot spotを病態と判断する前に、どのQC項目を確認するか。
Energy Loss、WSS、OSI、HelicityとTKE/MKEの違いを一文で言えるか。
次の記事では、Energy Lossを「圧損・散逸・TKEとの違い」から分解する。

次の記事への橋渡し

次はEnergy Lossを独立記事にするのが自然です。TKEが高い場所で本当に損失が大きいのか、圧較差や散逸率とどう違うのかを、数式と図解で分けていきます。