发明名称 一种获取水泥基材料中氯离子扩散和结合参数的方法
摘要 本发明公开了一种获取水泥基材料中氯离子扩散和结合参数的方法,该方法只需要一条总氯离子浓度曲线,通过回归计算即可同时获取氯离子的有效扩散系数和吸附等温曲线参数。该方法包括如下步骤:通过非稳态扩散试验得到一条总氯离子浓度分布曲线;尝试不同的待求参数组合,计算获得理论总氯离子浓度分布曲线;用非线性的最小二乘法对理论曲线和试验曲线进行拟合回归,获得拟合度最高的一套参数值;计算各参数的置信区间;进一步计算得到吸附曲线,表观扩散系数随时间和空间分布曲线,总氯离子、结合氯离子、自由氯离子在时间和空间上的分布曲线。本发明方法不仅可以同时得到扩散参数和结合参数,还可以得到各参数的区间。
申请公布号 CN105973757A 申请公布日期 2016.09.28
申请号 CN201610274452.8 申请日期 2016.04.28
申请人 东南大学 发明人 万克树;司先军;徐自强
分类号 G01N13/00(2006.01)I;G01N15/00(2006.01)I;G01N17/00(2006.01)I 主分类号 G01N13/00(2006.01)I
代理机构 南京苏高专利商标事务所(普通合伙) 32204 代理人 李倩
主权项 一种获取水泥基材料中氯离子扩散和结合参数的方法,其特征在于:包括如下步骤:步骤1,通过氯离子的非稳态扩散试验,得到水泥基材料中总氯离子浓度随扩散深度x和/或扩散时间t的分布曲线;步骤2,考虑氯离子在水泥基材料中的扩散和结合,建立如下的传输方程:方程(1)是总氯离子C<sub>t</sub>(x,t)、自由氯离子C<sub>f</sub>(x,t)和结合氯离子C<sub>b</sub>(x,t)浓度的关系式,其中,m为自由氯离子与总氯离子单位换算的参数;C<sub>t</sub>(x,t)=C<sub>b</sub>(x,t)+m·C<sub>f</sub>(x,t)      (1);方程(2)是自由氯离子C<sub>f</sub>(x,t)和结合氯离子C<sub>b</sub>(x,t)浓度之间的函数关系,采用氯离子吸附等温线来表示,其中,α和β为结合参数:<maths num="0001"><math><![CDATA[<mrow><msub><mi>C</mi><mi>b</mi></msub><mo>=</mo><mfrac><mrow><msub><mi>&alpha;C</mi><mi>f</mi></msub></mrow><mrow><mn>1</mn><mo>+</mo><msub><mi>&beta;C</mi><mi>f</mi></msub></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow><mo>;</mo></mrow>]]></math><img file="FDA0000977204120000011.GIF" wi="1054" he="108" /></maths>方程(3)采用菲克第二定律及其边界条件和初值条件,其中,D<sub>a</sub>是水泥基材料孔溶液中氯离子的表观扩散系数,C<sub>0</sub>(mol/L)是浸泡溶液的浓度,L是无穷大;<maths num="0002"><math><![CDATA[<mrow><mfrac><mrow><mo>&part;</mo><msub><mi>C</mi><mi>f</mi></msub><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>t</mi><mo>)</mo></mrow></mrow><mrow><mo>&part;</mo><mi>t</mi></mrow></mfrac><mo>=</mo><mfrac><mo>&part;</mo><mrow><mo>&part;</mo><mi>x</mi></mrow></mfrac><mrow><mo>(</mo><msub><mi>D</mi><mi>a</mi></msub><mo>&CenterDot;</mo><mfrac><mrow><mo>&part;</mo><msub><mi>C</mi><mi>f</mi></msub><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>t</mi><mo>)</mo></mrow></mrow><mrow><mo>&part;</mo><mi>x</mi></mrow></mfrac><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000977204120000012.GIF" wi="554" he="103" /></maths>C<sub>f</sub>(x,0)=0        (3)C<sub>f</sub>(0,t)=c<sub>Q</sub>C<sub>f</sub>(L,0)=0;方程(4)是D<sub>a</sub>和水泥基材料孔溶液中氯离子有效扩散系数D<sub>e</sub>之间的关系式;其中,<img file="FDA0000977204120000013.GIF" wi="157" he="109" />是氯离子结合能力,可由方程(2)获得;<maths num="0003"><math><![CDATA[<mrow><msub><mi>D</mi><mi>a</mi></msub><mo>=</mo><mfrac><msub><mi>D</mi><mi>e</mi></msub><mrow><mn>1</mn><mo>+</mo><mfrac><mn>1</mn><mi>m</mi></mfrac><mo>&CenterDot;</mo><mfrac><mrow><mo>&part;</mo><msub><mi>C</mi><mi>b</mi></msub><mrow><mo>(</mo><mrow><mi>x</mi><mo>,</mo><mi>t</mi></mrow><mo>)</mo></mrow></mrow><mrow><mo>&part;</mo><msub><mi>C</mi><mi>f</mi></msub><mrow><mo>(</mo><mrow><mi>x</mi><mo>,</mo><mi>t</mi></mrow><mo>)</mo></mrow></mrow></mfrac></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow><mo>;</mo></mrow>]]></math><img file="FDA0000977204120000014.GIF" wi="1099" he="135" /></maths>方程(5)是方程(4)的解析解,<maths num="0004"><math><![CDATA[<mrow><mfrac><mrow><msubsup><mo>&Integral;</mo><msub><mi>C</mi><mn>0</mn></msub><msub><mi>C</mi><mi>f</mi></msub></msubsup><msub><mi>D</mi><mi>a</mi></msub><mrow><mo>(</mo><msub><mi>C</mi><mi>f</mi></msub><mo>)</mo></mrow><msub><mi>dC</mi><mi>f</mi></msub></mrow><mrow><msubsup><mo>&Integral;</mo><msub><mi>C</mi><mn>0</mn></msub><mn>0</mn></msubsup><msub><mi>D</mi><mi>a</mi></msub><mrow><mo>(</mo><msub><mi>C</mi><mi>f</mi></msub><mo>)</mo></mrow><msub><mi>dC</mi><mi>f</mi></msub></mrow></mfrac><mo>=</mo><mi>e</mi><mi>r</mi><mi>f</mi><mrow><mo>(</mo><mfrac><mi>x</mi><msqrt><mrow><mn>4</mn><msub><mi>D</mi><mi>a</mi></msub><mrow><mo>(</mo><msub><mi>C</mi><mi>f</mi></msub><mo>)</mo></mrow><mi>t</mi></mrow></msqrt></mfrac><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow><mo>;</mo></mrow>]]></math><img file="FDA0000977204120000021.GIF" wi="1174" he="199" /></maths>步骤3,在一定的参数范围内,选择任意一套有效扩散参数D<sub>e</sub>和结合参数,根据方程(2)~(5)可计算获取自由氯离子C<sub>f</sub>(x,t)浓度分布曲线,再代入方程(1)~(2)可获得总氯离子C<sub>t</sub>(x,t)浓度分布曲线,然后用非线性的最小二乘法与步骤1得到的实验总氯离子浓度分布曲线进行拟合回归,最终获得拟合度最高的一套参数值;步骤4,假定拟合度最高的参数值对应的是真实的参数值,将其代入扩散方程计算得到一条理论曲线,将理论曲线作为输入曲线,其最高拟合度将是100%,由步骤3得到的最高拟合度与100%拟合度之间的参数范围即为每个参数的置信区间;步骤5,将步骤3得到的拟合度最高的参数值再次代入方程(1)~(5),可计算得到等温吸附曲线,表观扩散系数随扩散时间和扩散深度的分布曲线以及总氯离子浓度、结合氯离子浓度和自由氯离子浓度随扩散时间和扩散深度的分布曲线。
地址 210096 江苏省南京市玄武区四牌楼2号