模型說明
腹膜透析藥水的脫水(ultrafiltration, UF)與溶質清除,在生物物理學上主要由三個經典數學框架描述:Kedem-Katchalsky 不可逆熱力學方程、Pyle-Popovich 雙室動力學模型,以及 Rippe 的三孔模型(three-pore model)。本模擬器採用簡化的宏觀模型,參數以三孔模型族群曲線(Mujais & Vonesh)與 PET 分型範圍(Twardowski)擬合,再以台灣腎臟醫學會病例集的實際病人資料檢核。數值用來教學判讀方向與大小,不取代 PET、24 小時收集與臨床判斷。
1. 水分傳輸:三孔模型與修正 Starling 方程
三孔模型(Rippe 1991)把腹膜微血管壁視為有三種孔徑的半透膜。下標 B=血漿、D=透析液,k=孔徑(C 超小孔、S 小孔、L 大孔),i=溶質。
(1)腹腔液量變化(淨脫水率)
dVD/dt = Jv,C + Jv,S + Jv,L − L
(2)各孔徑的跨膜水通量(修正 Starling 方程)
Jv,k = αk · LpS · [ ΔP − σk,prot · Δπprot + Σi σk,i · RT · (CD,i − CB,i) ]
ΔP Pcap − IPP ≈ 17 − (1 + 3·VD) mmHg(VD 以公升計;2 L 時約 7 mmHg,留置量越大越小)
Δπprot 血漿蛋白膠體滲透壓,約 22–26 mmHg;σprot:超小孔 1、小孔約 0.99、大孔約 0.10
RT 37 °C 約 19.3 mmHg/(mmol/L)
σ 超小孔對所有溶質 σ=1;小孔對葡萄糖約 0.024–0.031、整體約 0.04
LpS 成人約 0.07–0.085 mL/min/mmHg(Rippe 1991 用 0.082)
L 腹腔淋巴回流:Rippe 1991 的直接淋巴流約 0.3 mL/min;臨床示蹤劑量到的「有效吸收」約 1–1.5 mL/min,含組織間隙吸收
膠體滲透壓項不能省略:少了它,等張液留置不會出現淨吸收,也算不出長留置的負脫水。超小孔的 σ=1,它看到的是全部溶質的滲透壓差;3.86%(4.25% dextrose)高張液留置第 1 小時,約一半脫水經超小孔,1.36% 或留置後期比例明顯較低。
2. 溶質傳輸:Patlak 方程與 Pyle-Popovich 模型
(3)跨膜溶質通量(Patlak/Kedem-Katchalsky)
Js,k = PSk · (CB − CD) + Jv,k · (1 − σk) · [ CB − f · (CB − CD) ]
f = 1/Pe − 1/(ePe − 1),Pe = Jv,k · (1 − σk) ÷ PSk (Pe → 0 時 f → 0.5)
與 PS·ΔC·Pe/(ePe−1) + Jv(1−σ)CB 完全等價
(4)Pyle-Popovich 腹腔溶質質量平衡(宏觀雙室模型)
d(VD CD)/dt = MTAC · (CB − CD) + QU · S · C̄ − L · CD
C̄ = CB − f · (CB − CD),S = 1 − σ(篩選係數),QU=經膜脫水率
3. 臨床估算:Garred/Waniewski MTAC、OCG 與自由水
(5)單次留置估算 MTAC(Waniewski 修正 Garred 公式)
MTAC = (V̄ ÷ t) · ln[ V01−F · (CB − CD,0) ÷ ( Vt1−F · (CB − CD,t) ) ]
Waniewski 建議 F=0.5;F=1 時體積項消失,還原為 Garred 公式
MTAC = (V̄ ÷ t) · ln[ (CB − CD,0) ÷ (CB − CD,t) ]
(6)葡萄糖滲透傳導度(OCG)與自由水(FWT)
OCG = LpS · σglu Rippe 參數 0.082 × 0.043 ≈ 3.5 μL/min/mmHg
FWT = UF60 − ΔMNa ÷ PNa (3.86% 葡萄糖、60 分鐘;ΔMNa 為移除的鈉量)
ISPD 2021 判讀:3.86% 葡萄糖 4 小時淨脫水 < 400 mL(或 2.27% < 100 mL)應懷疑脫水不足;3.86% 留置 1 小時鈉下降 ≤ 5 mmol/L 或鈉篩比 ≤ 0.03,提示超小孔或 OCG 功能低下。間質纖維化會同時降低 LpS 與 OCG(Morelle 2015)。
4. 本模擬器採用的簡化
葡萄糖藥水 dV/dt = Kuf × (C葡萄糖,透析液 − C葡萄糖,血漿) − L
Kuf 相當於公式(2)中 LpS × σ × RT 的合併值,三個孔道不分開計算
L 取臨床「有效吸收」1.3–1.6 mL/min(已包含膠體滲透壓與靜水壓的淨效果)
愛多尼爾 dV/dt = A × e−t/τ − L愛多尼爾 (以經驗曲線取代膠體滲透、分子量級分與澱粉酶水解)
葡萄糖吸收 dM葡萄糖/dt = −MTAC葡萄糖 × (C透析液 − C血漿)
溶質清除 公式(4):MTAC × (C血漿 − C透析液) + 對流項
脫水時對流項 = S × C血漿 × J(相當於 f=0);回吸時 = C透析液 × J
D/P = 透析液濃度 ÷ 血漿濃度;D/D₀ = 留置 t 小時葡萄糖濃度 ÷ 剛注入時濃度
每分鐘逐步計算;殘留液 200 mL;S:尿素 0.7、肌酸酐 0.6;血漿葡萄糖 6 mmol/L(糖尿病個案 10)
腹腔容量(腹內壓) 公式(2)的 ΔP = 17 − IPP;一般腹腔 IPP ≈ 1 + 3·V,2–3 L 時仍有 5–8 mmHg 的正向靜水壓
腹腔接近容量時壓力-容積曲線變陡,IPP 迅速逼近微血管壓,ΔP 趨近 0 再轉負
模擬器以容量係數 f 近似這段陡升:
Vcap = 坐姿注入上限 × 1.35(一般病人)
Vcap = 坐姿注入上限(多囊腎、疝氣等腹腔空間受限:注入上限就是剩餘空間)
f = 1 當 V腹腔 ≤ 0.8·Vcap
f = (Vcap − V腹腔) ÷ (0.2·Vcap) 當 0.8·Vcap < V腹腔 ≤ Vcap
dV/dt = f × 滲透流入 − 回吸(超過 Vcap 時回吸再加 1%/min × 超出量)
一般腹腔在 IPP ≈ 1 + 3·V 的線性段,每多 1 L 液量,脫水率只降約 LpS × 3 ≈ 0.2 mL/min,所以一般體型灌 2 L 時不受影響,族群曲線比對不變。腹腔空間小的病人(多囊腎)壓力-容積曲線很早就進入陡升段,愛多尼爾因為膠體滲透壓持續整個留置,最容易頂到上限:脫水量受剩餘空間限制,而不是藥水本身。注入量滑桿最多只能拉到注入上限,灌滿時脫水幾乎為零。陡升段的係數 1.35 與 0.8 為教學用假設值,尚無文獻直接量測。
因為三個孔道合併成一個滲透係數,模擬器呈現得出「葡萄糖先升後降、愛多尼爾持續緩升」的曲線,但無法呈現鈉篩效應與超小孔、小孔各自的貢獻。
5. 不同藥水的動力學比較
愛多尼爾的脫水量:Rippe & Levin 2000 預測未曾使用者 12 小時約 600 mL、長期使用者(血中已有代謝物)約 400 mL,並顯示脫水隨有效腹膜面積增加;Akonur 2015 預測長期使用者 9/14 小時約 385/460–509 mL,慢到快傳輸者之間差距小;Stachowska-Pietka 2024 的模擬也只見到傳輸型之間的小差異。本模擬器依 Mujais & Vonesh 族群曲線(14 小時約 710–810 mL,接近未曾使用者),14 小時低傳輸略高於高傳輸;長期使用者的實際脫水可能較低。
6. 各傳輸型參數與 PET 對照
括號內為 Twardowski 標準 PET 的 4 小時分型範圍。
7. 與族群淨脫水曲線的比對
比對值為 Mujais & Vonesh 2 L 注入淨脫水曲線(數位化)讀出的峰值、達峰與歸零時間;「歸零」為淨脫水回到 0 的時間。綠色為差距在容許範圍內,橘色為偏差較大。
8. 以病例集實際資料檢核
9. 臨床計算式
每週 Kt/V(腹膜)= 每日引流尿素量 ÷ 血漿尿素 × 7 ÷ V
每週 Kt/V(殘腎)= U/P 尿素 × 每日尿量 × 7 ÷ V;U/P 尿素預設 2.8
V(尿素分布容積,Watson)
男 = 2.447 − 0.09516 × 年齡 + 0.1074 × 身高 + 0.3362 × 體重
女 = −2.097 + 0.1069 × 身高 + 0.2466 × 體重
體表面積(Mosteller)= √(身高 × 體重 ÷ 3600)
每週肌酸酐清除 = (腹膜肌酸酐清除 + 殘腎尿素與肌酸酐清除平均) × 7 × 1.73 ÷ 體表面積
注入量:建議 1.25 L/m²、上限 1.5 L/m²(體表面積經驗法則);多囊腎約 1,000–1,200 mL(臨床經驗值)
尿素氮生成 UNA(g/日)= (蛋白質攝取 − 20.1) ÷ 7.5(由 Bergström PNA 公式反推)
穩定期 BUN(mg/dL)= UNA × 100 ÷ 每日總尿素清除(L)
U/P 尿素 2.8 取自病例集兩位病人實測(2.75、2.99)。第二關 10 位個案的傳輸型比例(高 1、中高 4、中低 4、低 1)參考台北榮總 898 位新進腹膜透析病人的分布(高 13%、中高 44.8%、中低 35.7%、低 6.5%)。
10. 限制
- 三孔合併為單一滲透係數,未模擬鈉篩效應、無溶質水傳輸與腹膜纖維化對滲透傳導度的影響;淋巴回流使用臨床有效吸收值,未拆分直接淋巴流與組織吸收。
- 愛多尼爾以經驗曲線表示,未拆分子量級分與澱粉酶水解(原始研究用 5–8 個級分),也未區分未曾使用者與長期使用者。
- 脫水曲線依 2 L 標準注入量的族群平均擬合;個別病人的淋巴回吸、導管引流與腹內壓差異很大。
- 血漿尿素、肌酸酐與葡萄糖視為固定值,沒有模擬一天內的濃度變化。
- APD 以每次灌入、引流共 25 分鐘估算,未模擬潮式透析與引流不完全。
- 殘腎清除以固定 U/P 比推估,實際應以 24 小時尿液收集計算;BUN 估算未計入透析液蛋白流失。
參考文獻
理論基礎
- Kedem O, Katchalsky A. Thermodynamic analysis of the permeability of biological membranes to non-electrolytes. Biochim Biophys Acta. 1958;27(2):229-246. doi:10.1016/0006-3002(58)90330-5
- Popovich RP, Moncrief JW, Nolph KD, Ghods AJ, Twardowski ZJ, Pyle WK. Continuous ambulatory peritoneal dialysis. Ann Intern Med. 1978;88(4):449-456. doi:10.7326/0003-4819-88-4-449
- Garred LJ, Canaud B, Farrell PC. A simple kinetic model for assessing peritoneal mass transfer in chronic ambulatory peritoneal dialysis. ASAIO J. 1983;6:131-137.(未收錄於 PubMed)
- Flessner MF, Dedrick RL, Schultz JS. A distributed model of peritoneal-plasma transport: theoretical considerations. Am J Physiol. 1984;246(4 Pt 2):R597-R607. doi:10.1152/ajpregu.1984.246.4.R597
- Rippe B, Stelin G. Simulations of peritoneal solute transport during CAPD. Application of two-pore formalism. Kidney Int. 1989;35(5):1234-1244. doi:10.1038/ki.1989.115
- Rippe B, Stelin G, Haraldsson B. Computer simulations of peritoneal fluid transport in CAPD. Kidney Int. 1991;40(2):315-325. doi:10.1038/ki.1991.216
- Waniewski J, Werynski A, Heimbürger O, Lindholm B. Simple models for description of small-solute transport in peritoneal dialysis. Blood Purif. 1991;9(3):129-141. doi:10.1159/000170009
- Vonesh EF, Lysaght MJ, Moran J, Farrell P. Kinetic modeling as a prescription aid in peritoneal dialysis. Blood Purif. 1991;9(5-6):246-270. doi:10.1159/000170024
- Vonesh EF, Story KO, O'Neill WT. A multinational clinical validation study of PD ADEQUEST 2.0. Perit Dial Int. 1999;19(6):556-571. PMID: 10641777
- Ni J, Verbavatz JM, Rippe A, et al. Aquaporin-1 plays an essential role in water permeability and ultrafiltration during peritoneal dialysis. Kidney Int. 2006;69(9):1518-1525. doi:10.1038/sj.ki.5000285
- Morelle J, Sow A, Hautem N, et al. Interstitial fibrosis restricts osmotic water transport in encapsulating peritoneal sclerosis. J Am Soc Nephrol. 2015;26(10):2521-2533. doi:10.1681/ASN.2014090939
腹膜功能評估與處方
- Twardowski ZJ, Nolph KD, Khanna R, et al. Peritoneal equilibration test. Perit Dial Bull. 1987;7(3):138-147.(未收錄於 PubMed)
- Mujais S, Vonesh E. Profiling of peritoneal ultrafiltration. Kidney Int Suppl. 2002;(81):S17-S22. doi:10.1046/j.1523-1755.62.s81.4.x
- Morelle J, Stachowska-Pietka J, Öberg C, et al. ISPD recommendations for the evaluation of peritoneal membrane dysfunction in adults. Perit Dial Int. 2021;41(4):352-372. doi:10.1177/0896860820982218
- Brown EA, Blake PG, Boudville N, et al. ISPD practice recommendations: prescribing high-quality goal-directed peritoneal dialysis. Perit Dial Int. 2020;40(3):244-253. doi:10.1177/0896860819895364
- van Biesen W, Heimburger O, Krediet R, et al. Evaluation of peritoneal membrane characteristics: clinical advice for prescription management by the ERBP working group. Nephrol Dial Transplant. 2010;25(7):2052-2062. doi:10.1093/ndt/gfq100
- Chou YH, Chen YT, Chen JY, et al. Baseline peritoneal membrane transport characteristics are associated with peritonitis risk in incident peritoneal dialysis patients. Membranes (Basel). 2022;12(3):276. doi:10.3390/membranes12030276
愛多尼爾
- Mistry CD, Gokal R, Peers E. A randomized multicenter clinical trial comparing isosmolar icodextrin with hyperosmolar glucose solutions in CAPD. Kidney Int. 1994;46(2):496-503. doi:10.1038/ki.1994.300
- Rippe B, Levin L. Computer simulations of ultrafiltration profiles for an icodextrin-based peritoneal fluid in CAPD. Kidney Int. 2000;57(6):2546-2556. doi:10.1046/j.1523-1755.2000.00114.x
- Akonur A, Holmes CJ, Leypoldt JK. Predicting the peritoneal absorption of icodextrin in rats and humans including the effect of α-amylase activity in dialysate. Perit Dial Int. 2015;35(3):288-296. doi:10.3747/pdi.2012.00247
- Stachowska-Pietka J, Waniewski J, Olszowska A, et al. Modelling of icodextrin hydrolysis and kinetics during peritoneal dialysis. Sci Rep. 2023;13:6526. doi:10.1038/s41598-023-33480-w
- Stachowska-Pietka J, Waniewski J, Olszowska A, et al. Can one long peritoneal dwell with icodextrin replace two short dwells with glucose? Front Physiol. 2024;15:1339762. doi:10.3389/fphys.2024.1339762
臨床計算與資料來源
- Watson PE, Watson ID, Batt RD. Total body water volumes for adult males and females estimated from simple anthropometric measurements. Am J Clin Nutr. 1980;33(1):27-39. doi:10.1093/ajcn/33.1.27
- Mosteller RD. Simplified calculation of body-surface area. N Engl J Med. 1987;317(17):1098. doi:10.1056/NEJM198710223171717
- Bergström J, Heimbürger O, Lindholm B. Calculation of the protein equivalent of total nitrogen appearance from urea appearance. Which formulas should be used? Perit Dial Int. 1998;18(5):467-473. PMID: 9848623
- Blake PG, Bargman JM, Brimble KS, et al. Clinical practice guidelines and recommendations on peritoneal dialysis adequacy 2011. Perit Dial Int. 2011;31(2):218-239. doi:10.3747/pdi.2011.00026
- 台灣腎臟醫學會。腹膜透析病例集 2021(TSN PD Casebook 2021)。個案 PD01 之 PET 與 24 小時收集資料。