مقدمة: من الخيار المملح إلى الأمن المائي العالمي — العملاق الخفي للتناضح (الأسموزية)

عندما نرش ملح الطعام على شرائح الخيار الرقيقة، تتساقط قطرات الماء من الأنسجة النباتية خلال دقائق. وفي الطب السريري، لو حُقن مريض بالخطأ بماء مقطر نقي وريدياً بدلاً من المحلول الملحي الفيزيولوجي ($0.9\,\text{wt}\%$ $\text{NaCl}$)، لاندفع الماء بقوة هائلة إلى داخل كريات الدم الحمراء حتى تنفجر أغشيتها الخلوية في حالة قاتلة تُعرف باسم انحلال الدم (Hemolysis). وعلى العكس، فإن البحار الغريق الذي يشرب ماء البحر يصاب بجفاف تناقضي قاتل لأن كليتيه تضطران لسحب ماء من خلايا جسمه لطرد فائض ملح البحر أكثر من حجم الماء الذي شربه.

تخضع جميع هذه الظواهر لعملية التناضح (الأسموزية - Osmosis)، وتُقاس قوتها الدافعة الديناميكية الحرارية بواسطة الضغط الاسموزي ($\Pi$). وتبلغ قيمة الضغط الاسموزي لمياه البحر القياسية ($35\,\text{g/L}$) نحو $27\text{ إلى }28\,\text{atm}$ ($\approx 2.75\,\text{MPa}$ أو $27.5\,\text{bar}$)، وهو ما يعادل ميكانيكياً الضغط الهيدروستاتيكي على عمق $280\,\text{متراً}$ تحت سطح المحيط! وبتطبيق ضغوط هيدروليكية تتراوح بين $55\text{ و}80\,\text{bar}$ على أغشية بوليمرية فائقة الرقة، تعكس محطات التناضح العكسي لمياه البحر (SWRO) هذه العملية الطبيعية لإنتاج أكثر من 100 مليون متر مكعب من مياه الشرب العذبة يومياً.

الفصل 1: تفنيد خرافة «تدرج تركيز الماء» — الكمون الكيميائي ($\mu_w$) وإنتروبيا المزج

1.1 لماذا يفشل التفسير المدرسي القائل بـ «انتشار الماء وفق تركيزه الحجمي»؟

تذكر معظم الكتب المدرسية أن الماء ينتشر ببساطة من حيث يكون تركيزه الحجمي ($c_w$ بوحدة $\text{mol/L}$) مرتفعاً إلى حيث يكون منخفضاً. هذا التفسير خاطئ فيزيائياً وكيميائياً:

1. الكثافة المولية الحجمية للماء ($c_w$) لا تحدد اتجاه التدفق الاسموزي: تبعاً للحجم المولي الجزئي للمذاب ($\bar{V}_s$) والتقبض الكهربائي، يتدفق الماء حصراً وفق كسره المولي $x_w$ وفعاليته الديناميكية الحرارية $a_w = \gamma_w x_w$.
2. انتشار فيك البسيط لا يفسر التوقف التام وانعكاس التدفق بالضغط الهيدروستاتيكي ($\Delta P = \Pi$): بما أن الماء السائل غير قابل للانضغاط تقريباً ($\kappa_T \approx 4.5 \times 10^{-5}\,\text{bar}^{-1}$)، فإن تطبيق ضغط $27.5\,\text{bar}$ لا يزيد كثافته إلا بنسبة $0.12\%$. ومع ذلك، فإن $\Delta P = \Pi$ يوقف التدفق تماماً و$\Delta P > \Pi$ يعكسه.
3. التناضح عبر الأغشية المسامية هو تدفق حملي جماعي، والضغط الاسموزي $\Pi$ خاصية تجميعية (Colligative property) لا تعتمد على حجم المسام أو مادة الغشاء.

1.2 الكمون الكيميائي لغيبس ($\mu_w$) والميكانيكا الإحصائية لإنتروبيا المزج ($\Delta S_{\text{mix}}$)

وفقاً لديناميكا غيبس الحرارية، عند ثبات درجة الحرارة $T$ والضغط $P$، تنتقل أي مادة كيميائية متحركة تلقائياً من منطقة ذات كمون كيميائي $\mu_i$ أعلى إلى منطقة ذات كمون كيميائي أدنى:

وللماء في المحلول مقارنة بالماء النقي $\mu_w^*(T, P_0)$:

حيث $\bar{V}_w \approx 18.07\,\text{cm}^3/\text{mol}$ هو الحجم المولي الجزئي للماء، و$x_w = \frac{n_w}{n_w + \sum \nu_k n_k} = 1 - x_s < 1$ هو الكسر المولي للماء. وفي المحلول المثالي ($\Delta H_{\text{mix}} = 0, \gamma_w = 1$)، ينشأ الانخفاض $\Delta \mu_w = RT \ln x_w < 0$ حصراً من زيادة إنتروبيا المزج ($\Delta S_{\text{mix}}$):

ولأن الغشاء شبه المنفذ يمنع أيونات الملح من العبور إلى الماء النقي، تعبر جزيئات الماء الغشاء لتخفيف الملح المحجوز، مما يعظم الإنتروبيا التشكيلية للكون ($S = k_B \ln W$) ويقلل طاقة غيبس الحرة $G = H - TS$.

1.3 الآلية الحركية لديباي–كرامر: انخفاض ضغط المذيب الموضعي عند فوهات المسام

كما أثبت بيتر ديباي (1923) وجيرالد مانينغ وإريك كرامر، عندما تصطدم جسيمات المذاب المتحركة حرارياً بالحاجز التنافري للغشاء $U_s(z)$ وترتد، فإنها تسحب جزيئات الماء معها نحو داخل المحلول:

وهذا يولد انخفاضاً ميكانيكياً موضعياً في ضغط الماء ($P_0 - \Pi$) عند مداخل المسام من جهة المحلول، فيسحب الماء النقي (الموجود عند ضغط $P_0$) عبر الغشاء.

  flowchart TD
    subgraph S1["1. التناضح الأمامي (Forward Osmosis): ΔP = 0 < Π"]
        L1["ماء نقي (L)<br/>ضغط P₀، كسر مولي x_w = 1<br/>كمون كيميائي μ_w(L) = μ_w*"]
        M1["غشاء شبه منفذ<br/>ارتداد المذاب عن الحاجز:<br/>ضغط الماء الموضعي يهبط إلى P₀ - Π"]
        R1["محلول ملحي (R)<br/>ضغط P₀، كسر مولي x_w < 1<br/>μ_w(R) = μ_w* + RT ln(γ_w x_w) < μ_w*"]
        L1 -->|"تدفق تلقائي للماء (Δμ_w > 0)<br/>إنتروبيا المزج والشفط عند المسام"| M1 --> R1
    end

    subgraph S2["2. الاتزان الاسموزي: ΔP = Π"]
        L2["ماء نقي (L)<br/>ضغط P₀، μ_w(L) = μ_w*"]
        M2["غشاء شبه منفذ<br/>تساوي الضغط عند المسام (P₀ = P₀)<br/>التدفق الصافي J_w = 0"]
        R2["محلول ملحي (R)<br/>تطبيق ضغط خارجي P₀ + Π<br/>μ_w(R) = μ_w* + V_w Π + RT ln(γ_w x_w) = μ_w*"]
        L2 <-->|"اتزان ديناميكي حراري (Δμ_w = 0)<br/>Π = -(RT / V_w) ln(γ_w x_w) ≈ iMRT"| M2 <--> R2
    end

    subgraph S3["3. التناضح العكسي (تحلية RO): ΔP > Π (55–80 bar)"]
        L3["مياه عذبة منتجة (L)<br/>ضغط P₀، μ_w(L) = μ_w*"]
        M3["غشاء RO TFC من البولي أميد<br/>نقل الذوبان-الانتشار:<br/>التدفق J_w = A(ΔP - ΔΠ)"]
        R3["محلول ملحي عالي الضغط (R)<br/>ضغط المضخة P₀ + ΔP (ΔP > Π)<br/>μ_w(R) = μ_w* + V_w(ΔP - Π) > μ_w*"]
        R3 -->|"انعكاس Δμ_w: μ_w(R) > μ_w(L)<br/>عصر الماء النقي من المحلول الملحي!"| M3 --> L3
    end

الفصل 2: الاشتقاق المتدرج لقانون فانت هوف (1887) وحساب الضغط الاسموزي لمياه البحر

2.1 الاشتقاق الديناميكي الحراري الدقيق

صاغ ياكوبوس هنريكوس فانت هوف (الحائز على أول جائزة نوبل في الكيمياء عام 1901) قانون الضغط الاسموزي استناداً إلى قياسات فيلهلم بفيفر بأغشية فيروسيانيد النحاس ($\text{Cu}_2[\text{Fe(CN)}_6]$):

1. تساوي الكمون الكيميائي للماء عند الاتزان بين الماء النقي ($L$، ضغط $P_0, x_w=1$) والمحلول ($R$، ضغط $P_0 + \Pi, x_w$):

2. المعادلة الديناميكية الحرارية الدقيقة (بافتراض عدم انضغاطية الماء):

3. مفكوك متسلسلة تايلور للمحاليل المخففة ($\gamma_w \to 1$، و$x_s = 1 - x_w = \frac{i n_s}{n_w + i n_s} \ll 1$ مع معامل فانت هوف $i$):

4. التحويل إلى التركيز المولاري $M = n_s / V$ (حيث $V \approx n_w \bar{V}_w$):

2.2 الحساب الدقيق لمياه البحر القياسية ($35\,\text{g/L}$ $\text{NaCl}$)

عند درجة حرارة $T = 298.15\,\text{K}$ ($25\,^\circ\text{C}$):

• التركيز المولاري: $M = \frac{35\,\text{g/L}}{58.44\,\text{g/mol}} \approx 0.60\,\text{mol/L}$.
• معامل فانت هوف: لملح $\text{NaCl}$ ($\nu = 2$)، المعامل الاسموزي $\phi \approx 0.925$ يعطي $i = \phi \nu \approx 1.85$ (الأسمولارية $\approx 1.11\,\text{Osm/L}$).
• الضغط الاسموزي ($R = 0.082057\,\text{L}\cdot\text{atm}/(\text{K}\cdot\text{mol})$):

وهذا يعادل عمود ماء بارتفاع $h = \frac{\Pi}{\rho g} \approx 280.6\,\text{متراً}$!

الفصل 3: التناضح البيولوجي — قنوات الأكوابورين (AQP1)، توترية كريات الدم، معادلة ستارلينغ وضغط الامتلاء

• قنوات الماء أكوابورين (AQP1، بيتر أغري، جائزة نوبل 2003): مسام بشكل الساعة الرملية مع مرشح فراغي $\text{ar/R}$ بقطر $2.8\,\text{\AA}$ يمرر $3 \times 10^9$ جزيء ماء في الثانية في صف واحد. ويجبر نمطان مركزيان من NPA (أسباراجين-برولين-ألانين) ثنائي القطب الكهربائي للماء على الدوران بزاوية $180^\circ$، مما يقطع سلسلة الروابط الهيدروجينية ويمنع تماماً تسرب البروتونات ($\text{H}_3\text{O}^+$) عبر آلية غروتهوس.
• توترية كريات الدم الحمراء: تحافظ الكريات على شكلها الطبيعي المقعر الوجهين في الوسط متساوي التوتر (Isotonic) ($0.9\,\text{wt}\%$ $\text{NaCl}$، $\approx 290\,\text{mOsm/L}$)، وتنكمش مشوكة (Crenation) في الوسط مفرط التوتر (Hypertonic)، وتنتفخ حتى تنفجر (Hemolysis) في الوسط ناقص التوتر (Hypotonic).
• معادلة ستارلينغ (1896) في الشعيرات الدموية: $$J_v = L_p S \left[ (P_c - P_i) - \sigma (\pi_c - \pi_i) \right]$$ يضمن الضغط الغرواني (الأونكوتي $\pi_c \approx 25\text{–}28\,\text{mmHg}$) الناتج عن ألبومين المصل ($66.5\,\text{kDa}$) إعادة امتصاص السوائل في الطرف الوريدي؛ ويؤدي نقصه إلى الوذمة (Edema) والاستسقاء (Ascites).
• ضغط الامتلاء النباتي ($\Psi_w = \Psi_s + \Psi_p + \Psi_g + \Psi_m$): يتحمل الجدار الخلوي ضغط امتلاء داخلي $\Psi_p \approx +0.8\text{ إلى }+1.5\,\text{MPa}$ يحافظ على انتصاب النباتات العشبية ويمنع البلزمة (Plasmolysis).

الفصل 4: تحلية المياه بالتناضح العكسي (RO) — أغشية البولي أميد المركبة (TFC) ونموذج الذوبان-الانتشار

بعد الأغشية غير المتماثلة للوب وسوريراجان (1960)، ابتكر جون كادوت في أواخر السبعينيات الغشاء المركب ذا الغشاء الرقيق (TFC) من البولي أميد عبر البلمرة البينية بين ميتا-فينيلين ثنائي الأمين (MPD) وكلوريد التريميسويل (TMC)، مكوناً طبقة فعالة متشابكة بسمك $100\text{ إلى }200\,\text{nm}$.

ووفقاً لنموذج الذوبان-الانتشار (Wijmans & Baker, 1995):

• تدفق الماء: $J_w = A(\Delta P - \Delta \Pi)$ حيث $A = \frac{D_w K_w c_w \bar{V}_w}{\delta R T}$.
• تدفق الملح: $J_s = B(C_m - C_p)$ حيث $B = \frac{D_s K_s}{\delta}$ (مستقل تماماً عن الضغط $\Delta P$!). ولذلك تقترب نسبة رفض الملح $R_s = \frac{A(\Delta P - \Delta \Pi)}{A(\Delta P - \Delta \Pi) + B}$ من $100\%$ مع زيادة $\Delta P$.
• حاجز نزع التميه وشذوذ البورون: يُرفض أيون $\text{Na}^+$ المتميه ($r_{\text{hyd}} \approx 3.58\,\text{\AA}$ مقابل $1.38\,\text{\AA}$ لجزيء $\text{H}_2\text{O}$) بسبب حاجز طاقة نزع التميه والاستبعاد العزلي لبورن. أما حمض البوريك المتعادل $\text{B(OH)}_3$ ($\text{p}K_a \approx 9.2$) عند $\text{pH} \approx 8.0$ فهو غير مشحون ويعبر الغشاء بتكوين روابط هيدروجينية (مما يتطلب رفع $\text{pH} > 10$ بواسطة $\text{NaOH}$ في مرحلة RO ثانية لتحويله إلى أيون البورات المشحون $\text{B(OH)}_4^-$).

  flowchart TD
    subgraph SWRO["محطة حديثة لتحلية مياه البحر بالتناضح العكسي (SWRO) مع مبادل ضغط أيزوباري (PX)"]
        Feed["مأخذ مياه البحر<br/>ملوحة 35 g/L، Π₀ ≈ 27.5 bar<br/>التدفق Q_f (100%)"]
        Pre["معالجة أولية UF/MF + مانع ترسب<br/>+ إزالة الكلور بـ NaHSO₃ (Cl₂ = 0)"]
        Split["مقسم التدفق"]
        HPP["مضخة الضغط العالي<br/>تضغط 50% من Q_f إلى 65 bar"]
        PX["مبادل الضغط الأيزوباري (PX)<br/>ينقل الضغط بكفاءة 96–98%<br/>إلى الـ 50% الأخرى من Q_f!"]
        BP["مضخة معززة (+2–3 bar)"]
        RO["وحدة أغشية RO TFC البولي أميد<br/>الضغط ΔP ≈ 65 bar > Π_exit (55 bar)"]
        Perm["المياه العذبة المنتجة (Permeate)<br/>استخلاص Y = 50%، TDS < 150 mg/L<br/>W_min = 1.06 kWh/m³ (الفعلي ≈ 2.1 kWh/m³)"]
        Brine["الرجيع الملحي عالي الضغط<br/>50% Q_f، 70 g/L (Π ≈ 55 bar)، ~63 bar"]
        Out["تصريف المحلول الملحي (~1 bar)"]

        Feed --> Pre --> Split
        Split -->|"50% Q_f"| HPP --> RO
        Split -->|"50% Q_f (ضغط منخفض)"| PX --> BP --> RO
        RO -->|"J_w = A(ΔP - ΔΠ)"| Perm
        RO -->|"الرجيع الملحي (63 bar)"| Brine --> PX --> Out
    end

الفصل 5: الحد الأدنى الديناميكي الحراري للطاقة ($1.06\,\text{kWh/m}^3$)، مبادلات الضغط (PX)، استقطاب التركيز والتلوث الحيوي

يُحسب الحد الأدنى للشغل العكوس للفصل $w_{\min}(Y)$ عند نسبة استخلاص $Y = V_p / V_f$ بتكامل $\Pi(y) = \frac{\Pi_0}{1 - y}$:

وعند $Y = 50\%$ ($\Pi_0 = 2.75\,\text{MPa}$)، نجد $w_{\min}(0.50) = 3.812\,\text{MJ/m}^3 \approx 1.06\,\text{kWh/m}^3$، مقارنة بنحو $\sim 650\,\text{kWh/m}^3$ للحرارة الكامنة لتبخير الماء. وبفضل مبادلات الضغط الأيزوبارية (PX) التي تسترجع $96\text{–}98\%$ من طاقة المحلول الملحي، انخفض الاستهلاك الفعلي لوحدة RO إلى $2.0\text{–}2.5\,\text{kWh/m}^3$.

وتتمثل العقبتان الهندسيتان الرئيسيتان في استقطاب التركيز ($\frac{C_m - C_p}{C_b - C_p} = \exp(J_w / k_m)$) والتلوث الحيوي (Biofouling) المتفاقم بسبب حساسية روابط الأميد في البولي أميد للكلور الحر (كلورة النيتروجين وإعادة ترتيب أورتون)، مما يحتم إزالة الكلور ببيسلفيت الصوديوم ($\text{NaHSO}_3$) قبل دخول الماء إلى الأغشية.