公開対象 / 安全性
本記事は MRI / 4D Flow / 統計 の教育目的で作成された一般化された解説です。特定患者、施設、製品仕様、非公開 WIP の代替ではありません。
まず一言で
WSS (Wall Shear Stress, 壁面せん断応力) は「血液が血管壁に接する直前で、どれだけ強く壁を擦っているか」を表す物理量である。単位は Pa (= N/m²)、生理的には大動脈で 0.1–1.5 Pa 程度のオーダー。
この一見地味な量が重要なのは、血管内皮細胞が WSS の 大きさ・方向・時間変動 を物理シグナルとして感知し、遺伝子発現・血管リモデリング・動脈硬化発生 を駆動するからである (Caro et al. 1971; Ku et al. 1985)。流速や圧較差を読むだけでは捉えられない「壁の経験する力」を、4D Flow MRI は速度ベクトル場から逆算する。
ただし WSS は速度勾配 (velocity gradient) に依存するため、空間分解能、Segmentation、ノイズ、平滑化 に極端に敏感で、典型的に 30–50% の系統的過小評価 が知られる。「絶対値は信用しないが、群間差・部位間差は読める」というのが現状のコンセンサスである。
---
1. WSS の物理的定義
1.1 連続体力学からの一次原理
血液をニュートン流体と仮定したときの WSS の定義は、壁面で血液要素にはたらく粘性応力テンソルの 接線成分 である。最も単純な形では:
τ_w = μ · (∂v_t / ∂n)|_wallここで:
τ_w: wall shear stress [Pa]μ: 粘性係数 (dynamic viscosity, 動的粘度とも) [Pa·s]、血液で約 3–4 mPa·sv_t: 壁面接線方向の血流速度 [m/s]∂/∂n: 壁面に垂直な方向 (wall-normal direction) の偏微分。慣習としてτ_wはその大きさで議論する
直感的には「壁から少し離れたところでどれだけ速く流れていて、それが壁にどれだけ近いところで急に止まるか (=速度勾配)」に粘性係数を掛けたものが WSS である。「壁に近いほど血液は遅い (no-slip)、遠いほど速い」という流体力学の基本仮定 (粘着条件) があるから、壁の直近に必ず速度勾配が立つ。WSS はその勾配の急峻さを定量化している。
1.2 単純化された解析解 — Poiseuille 流
直径 D の円管を通る定常 Poiseuille 流では、平均流速 V_mean から WSS が解析的に求まる:
τ_w = 8 · μ · V_mean / D例: 大動脈 (D = 30 mm, V_mean = 0.3 m/s, μ = 3.5 mPa·s) で τ_w ≈ 0.28 Pa。これは生理的値と一致する。
ただし、in vivo の血管はパルス流 (Womersley 流)、湾曲、分岐、テーパー、伸縮を伴い、Poiseuille の仮定は破綻する。これが「式は単純だが、実測は難しい」と言われる根本理由である。
1.3 ベクトル形式 — 三次元 WSS
実際の生体血管は曲がっており、壁面は曲面である。3D での WSS ベクトルは粘性応力テンソル T を壁面法線 n̂ に作用させ、接線成分のみを取り出す:
WSS_vec = T · n̂ - (n̂ · T · n̂) · n̂
T_ij = μ · (∂v_i/∂x_j + ∂v_j/∂x_i)4D Flow MRI が直接測れるのは速度場 v(x,t) だけなので、ここから空間勾配を取って T_ij を組み立て、壁面で評価し、法線方向と接線方向に分解する、という手続きを踏む。各ステップが誤差源になる。
---
2. なぜ WSS が臨床的に重要か — Mechanotransduction
血管内皮細胞は単なる「血液と平滑筋の仕切り」ではなく、WSS を感知する 力学センサー である。1971 年に Caro らが「動脈硬化 は 低 WSS 領域に好発する」という観察を報告して以来、半世紀の研究で以下のメカニズムが確立した:
- 高い層流 WSS (>1.5 Pa): NO 産生↑、抗炎症・抗血栓表現型、血管保護
- 低い WSS (<0.4 Pa): 内皮機能不全、接着分子発現↑、LDL 透過性↑、動脈硬化促進
- 振動性・乱流性 WSS: NF-κB 経路活性化、内皮細胞のターンオーバー↑
このため WSS は「測定可能な mechanotransduction 入力」として、以下の臨床問題に直結する:
- 二尖弁大動脈症 (BAV): 異常 WSS pattern が大動脈拡張を駆動するという仮説 (Mahadevia 2014, Guzzardi 2015)
- 大動脈解離リスク評価: 高 WSS 部位と内膜傷害の関連
- 脳動脈瘤: 高 WSS が瘤形成、低 WSS が瘤破裂と関連 (議論継続中)
- 冠動脈プラーク進展: 低 WSS 部位への lipid deposition
---
3. Newtonian 仮定と血液の非 Newtonian 性
§1 の式は「血液 = ニュートン流体 (せん断速度に依らず粘性が一定)」を仮定している。実際の血液は赤血球とその凝集 (rouleaux 形成) のため shear-thinning、すなわち低せん断速度で粘性が高く、高せん断速度で低くなる。この振る舞いは Carreau-Yasuda モデル でよく記述される:
μ(γ̇) = μ_∞ + (μ_0 - μ_∞) · [1 + (λ·γ̇)^a]^((n-1)/a)
γ̇ : せん断速度 [1/s]
μ_0 : ゼロせん断粘性 (高粘性極限, 約 50 mPa·s)
μ_∞ : 無限大せん断粘性 (低粘性極限, 約 3.5 mPa·s)[要原典確認: 上記パラメータ典型値は Cho & Kensey 1991、Boyd 2007 等で議論。一次PDF未確認の数値あり]
実用上は、大動脈や主要動脈のように γ̇ が高い (>100/s) 領域では Newtonian 近似が妥当 (μ ≈ μ_∞)。一方、小血管・低流速領域・瘤内部・心室内血流ではせん断速度が下がり、Newtonian 仮定で WSS を出すと体系的にずれる。研究レポートで「Newtonian 仮定」と明示されているか、Carreau モデルか、を必ず確認することが第一歩である。
---
4. 4D Flow MRI で WSS をどう推定するか
4D Flow MRI が測れるのは「voxel ごとの 3 方向速度ベクトル時系列」である。WSS を出すには次の 4 ステップを踏む:
4.1 Segmentation (壁の位置決め)
時間平均 magnitude / PC angiogram (PC-MRA) から血管 lumen を 3D segmentation する。出力は closed surface mesh (例: 三角形メッシュ)。これが「壁」である。
落とし穴: segmentation の閾値・smoothing が WSS に直接影響する。同じデータでも観察者間 variability が WSS 値で 10–20% 変動する報告がある (Potters 2015 系)。
4.2 Wall normal n̂ の決定
各 mesh 頂点で外向き法線ベクトルを計算する。この n̂ 方向の速度勾配を取るので、n̂ の精度がそのまま WSS の精度に伝播する。
4.3 速度勾配の計算
代表的な手法:
1. Direct gradient (Stalder 2008 — 4D Flow WSS の事実上の標準): 壁面近傍の速度を 2 次多項式 (Taylor 展開) で内挿し、壁面で ∂v/∂n を解析微分する。
2. B-spline smoothing (Potters 2015): 速度場を 3D B-spline で表現し、壁面で勾配評価。ノイズ耐性は高いが over-smoothing で WSS が下がる傾向。
3. Direct centered difference: 最も単純だがノイズに極端に弱く、研究的にはほぼ用いられない。
4.4 Tangential 成分の抽出
得られた速度勾配ベクトルから法線成分を引き、接線成分のみを残す。これが WSS ベクトル WSS_vec。スカラー指標としては:
- |WSS| : 大きさ (instantaneous)
- TAWSS : Time-Averaged WSS magnitude over cardiac cycle
- OSI : Oscillatory Shear Index (§6 で詳述)
- RRT : Relative Residence Time (§6 で詳述)
---
5. 系統的過小評価の正体
4D Flow MRI による WSS は、侵襲的なホットワイヤー流速計や CFD 比較で 30–50% 過小評価 されることが繰り返し報告されている (Petersson 2012, Cibis 2014 系)。原因は単一ではなく、複数の物理的・解析的因子の組み合わせである。
5.1 空間分解能依存性
WSS は壁面 1–2 mm の極薄境界層 (boundary layer) の速度勾配を読む。典型的 4D Flow voxel size は 2–3 mm 等方性で、boundary layer を 十分にサンプリングできない。voxel が大きいほど壁近傍の高速度勾配が「鈍化」され、WSS は系統的に下がる。Petersson 2012 は同じファントムで voxel size を変えながら WSS を測定し、voxel が 1.5 mm → 3 mm で WSS が 40% 以上低下することを示した。
5.2 Segmentation 依存性
壁の位置がわずかにズレるだけで、勾配を評価する点がずれて WSS が変わる。特に「壁を内側にずらす (under-segment)」と、より速度の高い点で勾配を取ることになり WSS が過大評価される一方、外側にずらす (over-segment) と低めに出る。実臨床では PC-MRA のコントラストが低く、観察者依存が大きい。
5.3 Partial volume effect
壁付近 voxel は血液信号と血管壁信号が混在する。位相は signal-weighted で平均化されるため、壁直近の速度は壁信号で「希釈」され、見かけの速度勾配が緩やかになる。これは voxel size を縮めない限り原理的に解決しない。
5.4 Smoothing と SNR のトレードオフ
ノイズ抑制のために空間 / 時間 smoothing をかけると、勾配が鈍り WSS が下がる。一方 smoothing なしだとノイズが勾配に乗り、空間的に荒い WSS map になる。Stalder 2008 protocol は「2 次多項式内挿 + 中程度 smoothing」のバランスで標準化を試みたが、ベンダー間・ソフト間で実装が異なる。
5.5 Wall motion
血管壁は心拍で 5–15% 径変化する。多くの WSS 解析は「壁を時間固定」で扱う。Aorta のように径変化が大きい血管では、これが追加誤差を生む。動的 segmentation は計算コストが高く、研究レベルに留まる。
---
6. 派生指標への橋 — OSI と RRT
WSS は時間変動するベクトル場なので、心拍 1 サイクルの平均だけでは「振動性」が捉えられない。次の 2 指標が併用される:
6.1 Oscillatory Shear Index (OSI)
He & Ku (1996) の定義:
OSI = (1/2) · [ 1 - |∫₀ᵀ WSS_vec dt| / ∫₀ᵀ |WSS_vec| dt ]- 値域: 0 (常に同方向) – 0.5 (完全に振動して時間平均がゼロ)
- 高 OSI 領域 = 「方向が揃わずぐらぐら」 = 内皮にとって不利な mechanical 環境
- 動脈硬化好発部位 (頚動脈分岐部、冠動脈分岐部後壁) と高 OSI 領域が一致する
直感: TAWSS が「平均してどれくらい強く擦っているか」、OSI は「方向が揃っているか、振り回されているか」。両者を組み合わせて初めて内皮環境が記述できる。
6.2 Relative Residence Time (RRT)
Himburg et al. (2004) の定義:
RRT = 1 / [ (1 - 2·OSI) · TAWSS ]- 単位は逆 Pa
- 「血液粒子が壁付近に滞留しがちか」の代理指標
- 高 RRT = 低 TAWSS かつ高 OSI 領域 = 動脈硬化リスク領域
OSI と TAWSS を 1 つにまとめた summary metric として、特に冠動脈・頚動脈研究で多用される。
---
7. 臨床応用と現状エビデンス
7.1 二尖弁大動脈 (BAV) と aortopathy
最も研究蓄積が厚い領域。BAV 患者では、弁形態 (right-left fusion vs right-non fusion) に応じて、上行大動脈の特定 region で 異常な高 WSS jet が生じる (Mahadevia et al. Circulation 2014)。Guzzardi 2015 は 4D Flow MRI で評価した高 WSS region と、実際に手術切除した大動脈組織の elastin degradation / matrix metalloproteinase 発現が一致することを示した。
→ 4D Flow WSS は「生検なしに局所組織リモデリングを推定できる」可能性。ただし WSS の絶対値ではなく patient-specific の異常 region 同定 が臨床的価値の中心。
7.2 大動脈瘤・解離
低 WSS 領域 = 血液滞留・血栓形成リスク、高 WSS 領域 = 内膜傷害リスク、という両端のロジックがある。van Ooij 2015 は atlas-based WSS mapping を提案し、健常者集団の正常 WSS map に対して、患者の局所異常を統計的に検出する手法を確立した。
7.3 脳動脈瘤
「高 WSS は瘤の発生・増大」「低 WSS は瘤の破裂」の二相仮説がある。CFD 主体だが 4D Flow MRI でも検証され、まだ「破裂予測指標として使える」段階には至っていない (2025 時点)。
7.4 冠動脈
Coronary 4D Flow は呼吸・心拍の二重 motion + 細い血管径で、WSS 推定はまだ研究段階。Coronary CTA + CFD の方が現実的。
---
8. 限界・落とし穴 (チェックリスト)
WSS 値を比較・解釈する際、以下を必ず確認する:
- vendor / sequence: PC-MRI / 4D Flow / WIP の implementation が異なれば WSS が直接比較できない
- voxel size: < 2 mm 等方性が望ましい (Petersson 系列)
- VENC: 過大な VENC は SNR を落とし、勾配ノイズを増やす
- 粘性モデル: Newtonian か Carreau-Yasuda か
- segmentation method: 観察者間 variability の報告があるか
- WSS algorithm: Stalder 2008 か、B-spline か、ベンダー独自か
- smoothing parameter: protocol 間で固定されているか
- 時間分解能: 心拍 1 周期で 20 時相以上が望ましい
- wall motion: 静止壁仮定か、動的 segmentation か
- 比較対象: 絶対値ではなく群間差・部位間差・atlas との偏差で議論しているか
「絶対値は ±50% 動く可能性がある」という前提を共有した上で、relative comparison として使うことが現時点のベストプラクティスである。
---
9. ZK化と次の記事への橋
この記事は次の atomic note に分解する (詳細は zk_notes/):
- WSS は粘性 × 壁法線速度勾配
- 4D Flow WSS は boundary layer サンプリング不足で 30–50% 過小評価
- Segmentation 依存性が観察者間 variability の主因
- 低 WSS 仮説 = atherogenesis 駆動 (Caro 1971 系)
- OSI は WSS 方向の振動性を 0–0.5 で表現
- RRT は OSI と TAWSS を統合した residence proxy
- BAV aortopathy では絶対値より局所異常 region 同定が臨床価値
連続記事の予定:
- OSI / RRT 単独記事: 振動性指標の数式・物理・臨床
- Helicity / Vorticity: 流れの「ねじれ」を読む
- Multi-VENC 完全解説: WSS 推定の前提となる velocity 精度
- Energy Loss との関係: 既出記事 260514_energy_loss_4dflow と内皮環境指標としての位置づけ比較
---
公開対象 / 安全性
- 施設名、患者情報、内部契約情報、メール、個人情報は含まない
- 一次文献の DOI / PubMed は
references.mdに集約 - 画像生成プロンプトは
image_prompts.jsonに分離 (Codex 側で GPT Image2 実行) - HTML 化、Cloudflare deploy、Discord 通知は Codex 担当