昨天把 δ ≈ 2.2 毫米這個數字標了存疑,理由很直接:D 借用的是氯離子的組織擴散係數(拿 Graiver 量的 NaCl 值頂替,因為查不到 NO₂⁻ 本身的組織擴散係數),k 是從「培根裡 80–90% 的亞硝酸鹽在熱處理後測不到」反推出來的假設值,兩個都是替代品,不是亞硝酸根本身在生肉裡的實測數字。今天把兩個替代品都換成真的——查到 Nordin(1969)的完整動力學公式,也查到 NO₂⁻ 本身在豬肉組織裡的擴散係數。結果是:兩個數字各自錯了一個數量級以上,但 δ 本身幾乎沒有變。這不是「昨天矇對了」的好消息,是一個更值得記錄的壞消息——兩個獨立的大誤差方向相反、量級相近,在開根號的比值裡幾乎完全抵銷。今天要把這個抵銷拆開來看,順便揪出一個更根本的問題:δ 這個穩態解,可能從一開始就不該套在「冷藏醃漬」這個階段上。

📖 學

一、真正的 k:Nordin(1969)的動力學公式

Kolari 與 Aunan(1972,美國肉品協會基金會會議論文,PDF 可公開下載)整理了 Nordin(1969,Canadian Institute of Food Technology Journal 2(2):79–85)對火腿裡亞硝酸鹽消耗速率的系統性研究。方法是把不同溫度、不同 pH 下量到的殘留亞硝酸鹽濃度取對數對時間畫直線(一階動力學:log[NaNO₂] = a + bt),每一組條件得到一個半衰期,再把半衰期對溫度與 pH 做二次迴歸,擬合出:

文中直接給出兩個定錨點驗證這條式子:pH 固定在 6.2 時,4.4°C(40°F)下半衰期是 513 小時,104.4°C(220°F)下只剩 1.6 小時;溫度固定在 76.7°C(170°F)時,pH 從 6.7 降到 5.6,半衰期從 12 小時掉到 4.9 小時——換算成公式裡 0.35×(6.7−5.6)=0.385,跟實測的 log₁₀(12/4.9)=0.389 幾乎完全吻合。這條迴歸式不是空話,是有內部一致性的。

代入 4.4°C、pH 6.2(文獻本身的定錨點,剛好貼近冷藏醃漬溫度):t½ = 512.9 小時,換算一階速率常數 k = ln2/t½ = 3.75 × 10⁻⁷ /s。用這個 k 回頭算「24 小時內消耗掉多少亞硝酸鹽」:kt = 3.75×10⁻⁷ × 86400 = 0.0324,C/C₀ = e^(−0.0324) = 0.968——24 小時只消耗了 3.2%。昨天拿「85% 在 24 小時內被消耗掉」反推出來的 k 是 2.2×10⁻⁵ /s,是今天這個真值的 58.5 倍。

二、為什麼差了 58 倍:一份真的培根實測表

差異的來源不難找。Kolari & Aunan 那份 PDF 裡有一組 1972 年 AMI 與 Swift 公司合作的培根實測數據,配方 120 ppm(跟 FSIS 現行注射醃上限一模一樣)那一列是:注射後、尚未煙燻 96 ppm(已損失 20%,對應「加進生肉當下就損失 20–25%」的說法);進煙燻室加熱後掉到 46 ppm(這一步再損失 52 個百分點);之後冷藏 7 天(7.2°C)剩 20 ppm,同樣 7 天但存放在 26.7°C 剩 6 ppm;28 天後 7.2°C 剩 10 ppm,26.7°C 只剩 1 ppm。

這張表把損失分成兩段,而兩段的溫度差了近 100°C:冷藏醃漬階段(4°C 左右)損失緩慢,煙燻/加熱階段(通常到 60–80°C 以上)損失暴增。昨天引用 Cassens & Mandigo「10–20% 殘留」那個數字,講的其實是熱加工之後的殘留量,反推出來的 k 因此混進了煙燻/蒸煮那段的高速消耗,卻被拿去套在「生肉冷藏醃漬」這個溫度差了近 100°C 的階段上。Nordin 公式裡溫度項是指數的,58.5 倍的誤差完全在合理範圍內——不是運氣不好,是把兩個溫度差一百度的階段當成同一回事。

三、真正的 D:NO₂⁻ 在豬肉組織裡到底擴散多快

Gómez、Sanjuán、Bon、Arnau 與 Clemente(2015,Journal of Food Engineering 149:188–194;同組 2019 年後續建模論文重複引用同一組數字)直接用豬半膜肌(Semimembranosus)圓柱樣品,在 2°C、7°C、12°C 三個溫度下量 NaNO₂ 與水的有效擴散係數。NO₂⁻ 的擴散係數範圍是 1.02×10⁻¹² 到 4.58×10⁻¹² m²/s(隨溫度上升而增加),活化能 60.32 kJ/mol;同一份研究裡水分擴散的活化能只有 32.24 kJ/mol——鹽類擴散對溫度的敏感度幾乎是水分擴散的兩倍。他們也確認了肌肉的方向性:垂直肌纖維方向的 NO₂⁻ 擴散比平行方向快。

這組數字比昨天拿來頂替的氯離子組織係數(Graiver 2006,70 g/L 鹵水、4°C,1.10×10⁻¹⁰ m²/s)小了 24 到 108 倍。這代表 Fox(1980)當年量到「亞硝酸根系統性比氯離子慢」這個現象,有很大一部分本來就是純擴散變慢——不是反應項偽裝出來的假象,是 NO₂⁻ 這個離子本身在肌肉組織裡走得比 Cl⁻ 慢兩個數量級。

四、重算 δ:1.22 倍的淨誤差,是巧合抵銷不是方法對

把兩個真值配到同一個溫度上。用活化能 60.32 kJ/mol 把 Gómez 等人在 2°C 量到的下限(1.02×10⁻¹²)內插到 4.4°C,得到 D = 1.28×10⁻¹² m²/s。配上前面算的 k = 3.75×10⁻⁷ /s(4.4°C、pH 6.2):

若肉的實際 pH 偏酸一點(生鮮豬肉常見的 5.6,而非 Nordin 定錨用的 6.2),k 會變成 6.09×10⁻⁷ /s,δ 跟著縮到約 1.45 mm。所以今天用兩份完全獨立的真實文獻重算出來的穿透深度,落在 1.45–1.85 毫米這個範圍。

跟昨天「已標存疑」拼出來的 2.2 毫米相比,誤差只有 20% 上下——但這不是「昨天的方法站得住腳」。拆開看誤差來源:D 被高估了 r_D = 1.10×10⁻¹⁰ / 1.28×10⁻¹² ≈ 86 倍,k 被高估了 r_k = 2.2×10⁻⁵ / 3.75×10⁻⁷ ≈ 58.5 倍。δ 只吃這兩個誤差的比值:√(86/58.5) ≈ 1.22 倍。2.2 / 1.22 ≈ 1.8 mm,跟直接重算的 1.85 mm 幾乎對上。這是統計上運氣好,不是方法論正確——如果只有一邊算錯(比如 D 用對了、k 用錯了),δ 會直接偏移 √58.5 ≈ 7.6 倍;只有兩邊剛好都錯、而且錯的方向剛好把比值修正回來,才會出現昨天那種「歪打正著」的結果。

五、更根本的問題:1/k ≈ 31 天,穩態可能根本沒到

δ = √(D/k) 是穩態解,只有在時間遠大於特徵時間 1/k 時才成立——這一點跟這個系列一直在強調的「先檢查 Fo 再選工具」是同一件事。1/k = 1/3.75×10⁻⁷ ≈ 2.66×10⁶ 秒 ≈ 31 天。任何在冷藏溫度下操作的醃漬——不管是浸漬醃培根的幾天,還是撒鹽乾醃的一兩週——都遠短於 31 天。也就是說,在冷藏醃漬這個溫度下,反應慢到系統可能根本沒進入穩態,δ = √(D/k) 這個框架本身可能問錯了階段:那幾天裡亞硝酸鹽的消耗多半可以忽略,行為更接近昨天更早幾天算的那種純擴散,而不是今天在算的反應-擴散穩態。真正會讓 1/k 掉到分鐘到小時級、讓反應-擴散穩態有意義的,是後面煙燻或蒸煮那段高溫加熱期——Barbieri、Bergamaschi 與 Franceschini(2013,Meat Science 93(2):282–286)直接在 55–72°C 這個加熱溫度區間測了不同起始亞硝酸鹽量(40–150 mg/kg)下的動力學參數,那才是 δ 真正該被拿來描述的溫度區間。連帶一個有趣的對照:09-22 那篇算過乾式醃製整條後腿的均勻化時間是 56 天——量級上跟這裡的 31 天已經接近,代表長期乾式醃製確實可能部分踩進反應限制的區域,快速醃漬(培根注射、幾天浸漬)則完全踩不到。

六、k 是四條路徑加總的表觀常數,不是單一反應

Kolari & Aunan 那份 PDF 裡有一節「Mechanism of Nitrite Losses」,列出至少四條會同時發生、彼此獨立的消耗路徑:跟蛋白質的巰基(SH)反應生成亞硝基硫醇 RS-NO(Mirna & Hoffmann 1959);在加熱下跟蛋白質/胜肽反應放出 N₂、NO、CO₂、NO₂ 等氣體(Lewis 1936;Greenwood 1940 的模型實驗甚至量到 NaCl 濃度從 0 拉到 20% 會讓損失從 68% 升到 98%);跟氧合肌紅蛋白反應生成硝酸根與變性肌紅蛋白(Greenberg et al. 1943);以及抗壞血酸鹽會急遽加速一氧化氮釋出(Ando et al. 1963:100 ppm 緩衝亞硝酸鹽溶液加了抗壞血酸鹽後瞬間掉到 26–35 ppm)。Nordin 那條 log-linear 公式量的不是任何一條路徑的速率常數,是把這些路徑加總後的表觀一階消耗率——這跟這個系列前幾天在鹽擴散那邊「把對流吸收進有效 D」是同一種統計陷阱的另一個版本:曲線擬合得夠好,不代表底層是單一基元反應。順帶一提,亞硝基肌紅蛋白(醃肉粉紅色的來源)生成之後也不是化學上穩定的終點——Møller 與 Skibsted(2004,Chemistry – A European Journal)量到它在水溶液模型系統裡會經過兩個連續的一階氧化步驟緩慢分解,活化焓分別在 88 與 121 kJ/mol 量級,溫度影響同樣劇烈(此為生理 pH 緩衝液模型系統的數字,套到真實肉品基質上是推論,已標存疑)。

七、回到 FSIS:δ 不吃濃度,只吃 D/k 的比值

今天複查 9 CFR 424.22(b)確認昨天引用的數字仍然成立:注射醃培根的亞硝酸鈉上限是 120 ppm ingoing(或等量 148 ppm 亞硝酸鉀),搭配 550 ppm 抗壞血酸鈉或異抗壞血酸鈉;浸漬醃同樣 120 ppm ingoing,以生肉去皮重回推;乾醃(dry cured bacon)要求「premeasured amount of cure mixture…completely covering the surfaces」,上限拉高到 200 ppm ingoing(或 246 ppm 亞硝酸鉀)。今天多查到一條昨天沒提到的細節:現行規範裡還有一條更低亞硝酸鹽的替代路徑——注射醃若搭配可發酵醣類與乳酸菌培養物,亞硝酸鈉可以降到 100 ppm ingoing,甚至 40–80 ppm 的範圍。

這條替代規則本身就是今天算的東西的直接應用。δ = √(D/k) 這個穿透深度不吃表面濃度 C_s,只吃 D 跟 k 的比值——加大配方裡的亞硝酸鹽濃度不會讓它擴散得更深,只會讓每一點的絕對濃度等比例升高,深度不變。這正是為什麼規範選擇的補救手段從來不是「放寬濃度上限」,而是改變遞送方式(注射,把它物理性地送進整塊肉的體積裡)或加另一個獨立的柵欄(可發酵醣類+乳酸菌帶來的酸化,補上被降低的亞硝酸鹽原本要頂的那部分抑菌力)——兩條路都繞開了「靠擴散把濃度送到深處」這件事,因為這件事本來就對加大劑量無感。

🧠 記

  • 今天最重要的更正:昨天的 k 測錯了溫度階段——用的是熱處理後殘留量反推的消耗率,套在冷藏醃漬這個溫度差近 100°C 的階段上。真正的冷藏醃漬期 k(Nordin 1969,4.4°C、pH 6.2)= 3.75×10⁻⁷ /s,24 小時只消耗 3.2% 的亞硝酸鹽,不是昨天假設的 85%。
  • 真正的 NO₂⁻ 組織擴散係數(Gómez et al. 2015/2019)是 1.02–4.58×10⁻¹² m²/s,比昨天借用的氯離子類比值小 24–108 倍,活化能 60.32 kJ/mol(水分擴散只有 32.24 kJ/mol)。這證實 Fox(1980)量到的「亞硝酸根比氯離子慢」有很大一部分本來就是純擴散變慢,不全是反應項的假象。
  • 用兩個真值重算 δ,落在 1.45–1.85 mm(視肉的 pH 5.6–6.2 而定),跟昨天的 2.2 mm 只差 20%——但這是巧合:D 被高估 86 倍、k 被高估 58.5 倍,兩者在開根號的比值裡幾乎抵銷(淨誤差只剩約 1.22 倍)。任何一邊單獨算錯,δ 都會偏移 7–9 倍,不會有這種歪打正著。
  • 更根本的問題:δ = √(D/k) 是穩態解,只有在 t ≫ 1/k 時才成立。冷藏溫度下 1/k ≈ 31 天,遠長於任何一支培根實際泡鹵水或撒鹽的時間——反應-擴散這個框架可能根本不該套在「冷藏醃漬」階段,那裡反應慢到近乎純擴散;真正需要它的是後面煙燻/蒸煮那幾十分鐘到幾小時,那裡 k 暴增、1/k 掉到分鐘級(Barbieri et al. 2013 在 55–72°C 測的正是這個階段)。
  • Nordin 的 k 公式是把至少四條獨立化學路徑(跟 SH 基反應成亞硝基硫醇、跟蛋白質/胜肽反應放氣體、跟氧合肌紅蛋白反應生成硝酸根、被抗壞血酸鹽加速消耗)加總後的表觀一階擬合,不是單一基元反應的速率常數——跟前幾天把對流吸收進「有效 D」是同一種統計陷阱。
  • FSIS 9 CFR 424.22(b)昨天引用的 120/200 ppm 上限複查仍然成立;新查到低亞硝酸鹽替代路徑(注射醃搭配可發酵醣類+乳酸菌可降到 100 ppm,甚至 40–80 ppm)——本質是拿另一個柵欄去補償被拿掉的那部分亞硝酸鹽。
  • δ 不吃表面濃度,只吃 D/k 的比值——想靠「多加亞硝酸鹽」讓它擴散得更深是無效的,穿透深度不會變,只有絕對濃度會變。這是規範選擇注射或強制均勻覆蓋、而非放寬濃度上限的物理理由。

✍️ 實踐

如果在家做需要亞硝酸鹽的醃肉(乾式培根、煙燻香腸等),不要用「加大劑量」去補償穿透不夠深的疑慮——今天算出來的東西直接說明加量無效,真正該做的是縮小肉的厚度、改用注射,或依規範作法把配方均勻覆蓋整個表面,而不是在濃度上加碼。

控制並記錄冷藏溫度仍然值得做,即使 δ 本身在冷藏溫度下可能還沒進入穩態:Nordin 公式與 AMI 那份實測培根表都指出,儲存溫度每差 20°C 左右,亞硝酸鹽消耗速率可以差三到十倍以上(表裡 7.2°C vs 26.7°C 那組數字就是活生生的對照)。維持穩定的低溫冷藏,除了食安的其他理由外,直接決定亞硝酸鹽在到達要保護的深度之前還剩多少。

想自己核對今天這幾個數字的人,Kolari & Aunan(1972)那份 PDF 是公開的,裡面完整重印了 Nordin 公式與 AMI 的培根實測表(注射前後、7 天/28 天在 7.2°C 與 26.7°C 下的殘留量),可以直接拿來對照其他二手引用是否一致——這比信任任何單一次轉述都可靠。

安全提醒照舊:亞硝酸鹽相關的操作涉及肉毒桿菌中毒風險,今天算的是機制,不是操作建議;實際配方劑量請以 FSIS 或具體法規指引與經驗豐富的醃肉資料為準,劑量計算錯誤或分佈不均可能造成嚴重食安風險。

🔗 延伸學習

💬 問 AI

我在看一篇食物科學筆記,昨天算出亞硝酸鹽在肌肉組織裡的反應-擴散穿透深度
δ = √(D/k) ≈ 2.2 mm,但 D 用的是氯離子的組織擴散係數類比、k 是從「熱處理後
殘留 10-20%」反推的假設值。今天查到兩個真值——Nordin (1969) 的動力學公式給出
冷藏溫度下 k ≈ 3.75×10⁻⁷ /s,Gómez et al. (2015) 給出 NO2- 本身的組織擴散係數
D ≈ 1.02-4.58×10⁻¹² m²/s——重算後 δ 落在 1.45-1.85 mm,跟昨天的猜測只差 20%,
但拆開來看是兩個各錯了一個數量級以上的誤差(D 高估 86 倍、k 高估 58.5 倍)在
開根號的比值裡幾乎抵銷。我還發現一個問題:δ 是穩態解,只在 t ≫ 1/k 時成立,
而冷藏溫度下 1/k ≈ 31 天,遠長於實際醃漬時間。
 
請幫我:
1. 檢查「δ 只在 t ≫ 1/k 時成立」這個判斷是否正確,如果冷藏醃漬階段確實還沒
   進入穩態,那在這個階段亞硝酸鹽的分布應該用什麼模型描述(純擴散?還是某種
   過渡態的反應-擴散解)?
2. 幫我評估「D 被高估 86 倍、k 被高估 58.5 倍、兩者幾乎抵銷」這個說法在統計上
   是否站得住腳,還是我把不確定的數字算得太精確、掩蓋了兩份文獻本身的實驗
   不確定度(Gómez 等人只在 2、7、12°C 三個溫度點量測,Nordin 的公式是 1969 年
   的手動迴歸)?
3. 如果我想更嚴謹地把「冷藏醃漬階段」跟「加熱處理階段」分開建模,分別用各自
   適用的擴散/反應機制銜接起來,這種分段模型在食品工程文獻裡有沒有更成熟的
   先例可以參考?