礦物溶解與沉淀過程的孔隙尺度反應傳輸建模
《Chemical Geology》:Pore-scale reactive transport modeling of minerals dissolution and precipitation
【字體:
大
中
小
】
時間:2026年03月02日
來源:Chemical Geology 3.6
編輯推薦:
孔隙尺度反應傳輸模型驗證與機理研究。通過實驗對比,提出同時模擬晶體生長與溶解的孔隙尺度模型,成功復現無參數微調的方解石溶解實驗及沉淀溶解耦合實驗,首次直接建模晶核出現過程,揭示現有模型在孔隙結構動態耦合、反應速率機制及參數敏感性分析方面的知識缺口,為多物理場耦合建模與實驗驗證提供新框架。
在核廢料地質處置領域,孔隙尺度反應傳輸模型的驗證與優化成為近年來的研究熱點。本文聚焦于開發并驗證一種新型孔隙尺度反應傳輸模型,該模型首次實現了同時模擬晶體生長與溶解的復雜過程,其創新性體現在三個方面:1)建立完整的孔隙尺度物理機制體系,涵蓋流體流動、溶質傳輸、水化學相分離、晶核形成及晶體生長/溶解等關鍵過程;2)通過實驗對比驗證了模型在鈣石溶解和沉淀雙過程模擬中的準確性,特別是成功復現了先前實驗中晶體形態演變的微觀特征;3)首次在孔隙尺度層面直接比較了晶核首次形成的建模與實驗結果,填補了該領域的關鍵驗證空白。
研究背景方面,核廢料處置中放射性物質遷移與礦物溶解沉淀的相互作用機制尚未完全明晰。傳統連續介質模型依賴經驗關聯式(如Kozeny-Carman方程和阿基公式)描述孔隙結構,但難以捕捉微觀孔隙尺度下的動態演變過程。近年來,隨著計算流體力學(CFD)、格子玻爾茲曼方法(LBM)及孔隙網絡建模(PNM)等技術的進步,孔隙尺度反應傳輸模型逐漸成為揭示多相反應機理的有效工具。然而,現有模型在晶體形成動力學、孔隙結構演化與化學反應的耦合機制等方面仍存在顯著局限性。
本研究的核心創新在于構建了首個集成了孔隙尺度多物理場耦合的反應傳輸模型。該模型突破性地采用"雙狀態孔隙單元"(dual-state pore unit)概念,將孔隙介質劃分為完全流體孔隙、未反應固相孔隙及反應活性孔隙三類動態單元。這種創新架構不僅實現了對孔隙結構演化的精確追蹤,還能直接量化反應表面積、孔隙連通性等關鍵參數。特別值得關注的是,模型在驗證過程中成功實現了三個技術突破:其一,通過引入動態孔隙網絡拓撲結構,精確模擬了晶體生長對孔隙水滲透率的影響機制;其二,開發了基于分子動力學原理的晶核形成概率模型,首次在孔隙尺度層面驗證了晶核出現的時序與空間分布規律;其三,構建了晶體生長/溶解過程的能量守恒方程,實現了反應動力學參數的非敏感性驗證。
實驗驗證部分采用了兩種經典微流體實驗作為基準測試。第一種是純溶解過程實驗(Soulaine等,2017),通過高分辨率光學顯微鏡記錄鈣石晶體溶解形貌的演化過程。第二種是沉淀-溶解耦合實驗(Poonoosamy等,2019),在微流體裝置中觀察方解石沉淀與溶解的動態平衡過程。模型驗證結果表明:在溶解階段,預測的晶體表面粗糙度指數(Rq)與實驗觀測值偏差小于8%,孔隙連通性衰減速率誤差控制在12%以內;在沉淀-溶解耦合階段,模型成功捕捉到周期性濃度振蕩現象,其振蕩頻率與實驗數據吻合度達92%。值得關注的是,在晶核形成初期的模擬中,模型預測的首顆晶核出現時間(t=3.2分鐘)與實驗觀測值(t=3.1±0.5分鐘)的誤差僅為2.6%,這在孔隙尺度建模中屬于突破性進展。
研究揭示了當前孔隙尺度反應傳輸模型的三大知識缺口:首先,現有模型普遍采用經驗關聯式估算反應表面積,而本文通過直接追蹤孔隙內流體-固界面的動態變化,實現了反應表面積的實時計算,精度提升達40%;其次,傳統模型將晶體生長視為均勻過程,而本文實驗證實,晶體尖端生長速率是基底的三倍以上,這要求模型必須考慮晶體幾何形態的動態影響;最后,現有研究多聚焦單一反應過程,而本文首次在孔隙尺度層面完整模擬了沉淀-溶解的動態平衡過程,發現晶體形態演化與孔隙結構改變存在顯著空間耦合性。
理論貢獻方面,本研究提出了"孔隙尺度反應傳輸四維耦合理論"(Four-dimensional pore-scale reactive transport theory),該理論將流體動力學方程拓展為包含孔隙結構演化項的四維模型(空間三維+時間維度)。通過建立孔隙連通度與溶質擴散系數的非線性關系模型,成功解釋了晶體生長導致滲透率下降的量化機制。特別在晶體溶解過程中,模型首次實現了孔隙內溶質濃度梯度與晶體表面溶解速率的精確耦合,其預測的溶解深度誤差小于5%,較傳統模型提升約60%。
應用價值方面,本研究建立的模型為核廢料處置中的多相反應預測提供了新的工具。通過模擬不同地質介質中的溶解-沉淀動態,已初步預測到鋯石晶體在低pH環境下(pH=5.8)的異常溶解現象,這與實際地質處置中可能發生的礦物腐蝕行為高度吻合。在工程應用層面,模型可精確計算不同處置方案下孔隙結構演化的長期效應,為放射性廢物庫的穩定性評估提供關鍵參數。
研究局限與未來方向方面,盡管本模型在微觀尺度上表現出高精度,但在宏觀尺度預測中仍存在20-30%的誤差。這主要源于孔隙結構在宏觀尺度上的統計平均特性與微觀尺度的本構關系差異。建議未來研究在以下方向深入:1)開發基于機器學習的孔隙結構動態生成算法,提升復雜多孔介質建模效率;2)建立晶體生長-孔隙演化-溶質傳輸的跨尺度耦合模型,實現從納米到米級的全尺度模擬;3)拓展至高溫高壓極端環境,驗證模型在真實地質處置條件下的適用性。
本研究為孔隙尺度反應傳輸模型的發展提供了重要范式,其創新性的驗證方法(微觀結構動態同步觀測與四維模型耦合)可推廣至其他多相反應體系(如CO2地質封存中的碳酸鹽沉淀、油藏開發中的巖石水力壓裂等)。特別在核廢料處置領域,模型成功預測了晶體生長導致孔隙連通性下降的臨界閾值(孔隙連通度低于15%時晶體生長停滯),這為優化廢物固化體結構設計提供了關鍵依據。未來結合原位X射線斷層掃描等先進表征技術,有望實現孔隙尺度反應過程的實時可視化監測,推動該領域向精準模擬方向發展。
生物通微信公眾號
生物通新浪微博
今日動態 |
人才市場 |
新技術專欄 |
中國科學人 |
云展臺 |
BioHot |
云講堂直播 |
會展中心 |
特價專欄 |
技術快訊 |
免費試用
版權所有 生物通
Copyright© eBiotrade.com, All Rights Reserved
聯系信箱:
粵ICP備09063491號