岩石由不同尺寸和形状的矿物颗粒胶结而成，其内部包含着随机分布的微裂纹和孔隙，而不同矿物颗粒之间的胶结程度也各有差异。由于岩石内部缺陷分布的演化具有随机性，遵循统计学规律，可以通过一定的概率分布进行描述。定义损伤为在某一应力作用下已经破坏的微元数目n与岩石初始状态下的微元总数目N之比。损伤变量即为微元破坏的概率，表示为已损伤的微元数目n与总微元数目N之比。

$D=\frac{n}{N}$(1)

在外部荷载作用下，岩石微元的破坏表现出随机性，并且根据Krajicinovic模型的描述，损伤变量可被定义为微元破坏概率小于等于 $F$的累积概率 $P$。如果我们设定微元破坏的概率密度函数为 $f\left(x\right)$，那么 $P\left(F\right)$可以被表述为 $f\left(x\right)$的累积分布函数，即：

$D_M=P\left(F\right)={\displaystyle {\int }_0^{F}f\left(x\right)\text{d}x}$(2)

研究表明，当建立岩石损伤统计本构模型时，基于岩石微元强度服从Weibull分布或正态分布的假定，来研究岩石破坏过程中的损伤演化特性，能够达到较好的研究效果，相对于正态分布，Weibull分布是一种较不保守的分布类型，它在一定程度上更真实地反映了随机变量的分布。鉴于此，本文假定岩石微元强度服从Weibull分布，其概率密度函数为

$\varphi \left(F\right)=\frac{m}{F_0}{\left(\frac{F}{F_0}\right)}^{m-1}\exp \left[-{\left(\frac{F}{F_0}\right)}^m\right]$(3)

式中， $F$表示微元强度， $m$和 $F_0$为Weibull分布的参数。

由式(2)得岩石在荷载作用下的损伤变量为

$D_M={\displaystyle \underset{0}{\overset{F}{\int }}\varphi \left(x\right)\text{d}x}=1-\exp \left[-{\left(\frac{F}{F_0}\right)}^m\right]$ (4)

当岩石处在水溶液环境中，岩石的组成成分中有部分会与水以及其中的氢离子或者其他离子反应生成新的化学物质，在此过程中，岩石组分溶解，孔隙率增大，岩石的结构也发生变化，导致岩石的强度以及弹性模量均发生变化。事实上，岩石的结构和组成极其复杂，理想情况下可以认为是由可溶性胶结性和不溶性材料制成的孔隙度岩石。水化作用会导致可溶性胶结材料的损失，减少不溶性颗粒之间的接触，最终影响岩石的宏观变形性能。因此，本文认为砂岩中可溶性胶结材料的损失是导致水化学腐蚀作用下岩石力学性能恶化的根本原因。根据李宁 [28] 之前的研究，一个基于微观结构的化学损伤变量可以将可溶性胶结材料的摩尔数变化来表征岩石水化学损伤：

$D_{CH}=\frac{\Delta S^d}{{\left(S^d\right)}_0}={\left[\frac{\Delta {\omega }^d}{{\left({\omega }^d\right)}_0}\right]}^{\frac{2}{3}}=1-\frac{S^d}{{\left(S^d\right)}_0}={\left[1-\frac{{\omega }^d}{{\left({\omega }^d\right)}_0}\right]}^{\frac{2}{3}}$ (5)

其中， ${\left(S^d\right)}_0$和 ${\left({\omega }^d\right)}_0$分别为初始无腐蚀时有效承载面积和岩石材料的总摩尔数； $\Delta S^d$和 $\Delta {\omega }^d$分别为化学腐蚀过程中有效承载面积的减小值和可溶胶结材料摩尔数的衰减值， $S^d$和 ${\omega }^d$分别为经过水化学腐蚀后的有效承载面积和可溶材料的摩尔数。考虑到岩石中还有不可溶材料，即不与水和酸性溶液反应的材料，这里的 ${\left({\omega }^d\right)}_0$表示岩石试样中所有材料的总摩尔数。

$\Delta {\omega }^d$可根据以下公式(10)~(12)得到。 ${\left({\omega }^d\right)}_0$可以根据岩石的矿物组成计算出来。

水与岩石之间存在着复杂而多样的相互作用，这种相互作用在地质历史上持续了很长时间，其机理涉及到多个方面。其中一个重要的方面是水和岩石中可溶性胶结物之间的化学反应，这种反应会影响岩石的力学性质和耐久性。水和岩石中可溶性胶结物之间的化学反应主要包括两个过程：(1) 溶解过程，即水分子与可溶性胶结物分子之间发生静电力或氢键作用，使得可溶性胶结物分子从固态转变为液态或气态，并进入到水中。(2) 反应过程，即水中的溶质与岩石中的可溶性胶结物分子之间发生化学键断裂或形成，产生新的物质，并改变原有物质的结构和组成。由于水在岩体内部具有一定的流动性，因此水和岩石之间的化学反应不仅发生在表面接触处，还会发生在内部孔隙处。这样就会导致岩石内部孔隙度、渗透率、强度、刚度等参数发生变化，从而造成岩石的化学损伤。为了描述和预测岩石在不同环境条件下(如酸度、温度、压力等)由于水和岩石之间的化学反应而引起的力学性质变化，需要建立合理的化学损伤本构模型。岩石与水之间的化学反应在这一复杂相互作用过程中起着决定性的作用。考虑到实验室中相对静止的浸泡条件和自然岩体中水流缓慢的情况，所提出的化学损伤模型忽略了水传输引起的化学损伤，仅将溶质溶解和化学反应作为化学损伤的主要原因。

岩石中可溶性胶结物与酸性溶液中的氢离子之间的化学反应是导致岩石化学损伤的主要原因之一。这种反应会降低砂岩的强度和刚度，增加其孔隙度和渗透率。反应速率取决于酸的类型和浓度、温度、压力和砂岩的矿物组成，岩石中的可溶性胶结物和氢离子反应方程式如下：

$m_i{A}_i+n_i{\text{H}}^{+}\to p_i{B}_i+q_i{C}_i$ (6)

其中， $A_i$和 ${\text{H}}^{+}$为反应物， $B_i$和 $C_i$代表生成物， $m_i$、 $n_i$、 $p_i$和 $q_i$为化学方程的平衡系数。

为了预测与 ${\text{H}}^{+}$的化学反应造成的损伤，应知道可溶性水泥的反应物损失。根据化学动力学理论，用 ${\text{H}}^{+}$的浓度变化速率表示的各化学反应的速度如下：

$v^i\left({\text{H}}^{+}\right)=k{c}^{x_i}\left({\text{H}}^{+}\right)c^{y_i}\left(A_i\right)$(7)

其中， $v^i\left({\text{H}}^{+}\right)$为 ${\text{H}}^{+}$与 $A_i$的反应速率， $c^{x_i}\left({\text{H}}^{+}\right)$和 $c^{y_i}\left(A_i\right)$分别表示 ${\text{H}}^{+}$和反应物 $A_i$的瞬时浓度。比例常数k为速率参数， $x_i$和 $y_i$分别为反应物的浓度参数。

反应物 $A_i$难溶于水，其在水溶液中的浓度是一个常数。由于存在多重反应，反应机理复杂， ${\text{H}}^{+}$的反应速率可以用来综合反映所有反应速率。 ${\text{H}}^{+}$的反应速率简化如下：

$v_{{\text{H}}^{+}}=\lambda c^x\left({\text{H}}^{+}\right)$(8)

即：

$-\frac{\text{d}c\left({\text{H}}^{+}\right)}{\text{d}t}=\lambda c^x\left({\text{H}}^{+}\right)$ (9)

由式(9)，可得：

$C_{{\text{H}}^{+}}\left(t\right)={\left({\left(C_{{\text{H}}^{+}}^0\right)}^{1-x}-\left(1-x\right)\lambda t\right)}^{\frac{1}{1-x}}$(10)

其中， $C_{{\text{H}}^{+}}^0$为 ${\text{H}}^{+}$的初始浓度， $C_{{\text{H}}^{+}}\left(t\right)$为时刻 $t$时的 ${\text{H}}^{+}$的浓度。

反映溶液pH值与 ${\text{H}}^{+}$浓度关系的公式如下：

$\text{pH}=-\lg c\left({\text{H}}^{+}\right)$(11)

为了预测岩石和 ${\text{H}}^{+}$的化学反应造成的损伤，应知道岩石可溶性化合物中反应物的损失。由式(6)与各反应物含量的关系，可以推断出与 ${\text{H}}^{+}$反应一定时间引起的岩石可溶性化合物 $D_i$摩尔数反应的量如下：

$\Delta {\omega }_1^d={\displaystyle \underset{i=1}{\overset{n}{\sum }}\frac{m_i{\alpha }_i}{n_i}\left(C_{{\text{H}}^{+}}\left(t\right)-C_{{\text{H}}^{+}}^0\right)V_{cs}}$ (12)

其中， $V_{cs}$表示化学溶液的体积。 ${\alpha }_i$表示所有能与 ${\text{H}}^{+}$反应的反应物中反应物 $A_i$的初始摩尔分数。

另一方面，化学溶液中的水可以与岩石中的相关物质发生反应。由于水量充足，反应在短时间内完成，时间效应可以忽略。因此，可溶性化合物的摩尔数损失可以定量地表示如下：

$\Delta {\omega }_2^d={\displaystyle \underset{i=1}{\overset{m}{\sum }}\frac{m_{Di}}{M_i}}$ (13)

其中， $m_{Di}$和 $M_i$分别为可溶性化合物 $D_i$的质量和摩尔质量。

由于任何物质与水反应实际上都是和水中水分子电离出来的 ${\text{H}}^{+}$以及 ${\text{OH}}^{-}$发生反应，岩石化合物不管是和酸性 ${\text{H}}^{+}$反应还是和水发生反应，本质上都是与 ${\text{H}}^{+}$发生反应，因此岩石化合物可溶性化合物的摩尔数损失总量应表示为：

(14)

结合式(5)、(12)~(14)，可以得到化学损伤的表达式如下：

$D_c={\left[\frac{{\displaystyle \underset{i=1}{\overset{n}{\sum }}\frac{m_i{\alpha }_i}{n_i}\left({10}^{-{\text{pH}}_0}-{10}^{-\text{pH}\left(t\right)}\right)V_{cs}}+{\displaystyle \underset{i=1}{\overset{m}{\sum }}\frac{m_{Di}}{M_i}}}{{\left({\omega }^d\right)}_0}\right]}^{\frac{2}{3}}$ (15)

岩石经循环的加热和冷却作用后，由于热应力的反复作用使得岩石在受荷载之前就存在微元破坏，从而产生热损伤。高温作用下，岩石内产生的热膨胀力以及结构热应力，使其内部细观结构(如微孔隙、微裂纹等)发生不可逆的变化，从而导致岩石宏观力学性质产生劣化。根据损伤力学理论，岩石热处理后，其损伤变量可用岩石物理力学参数的变化来表征。通过大量的文献及试验数据可以得出温度对岩石具有软化的作用，即温度升高，岩石的弹性模量减小，因此可采用弹性模量来定义热损伤变量：

$D_T=1-\frac{E_T}{E_0}$ (16)

式中， $E_T$为高温下岩石的弹性模量， $E_0$为常温下岩石的弹性模量。

由于岩石在温度、化学腐蚀和荷载作用下会表现出不同的损伤特性，因此，围压–化学腐蚀–温度耦合作用下岩石的总损伤并不是各种单一损伤变量的叠加。用缺陷材料的占比来分别表示岩石在三种单一状态下的损伤，缺陷材料的占比及缺陷材料体积占材料总体积的比值，假设初始状态未损伤岩石，高温损伤岩石，化学和高温耦合损伤岩石和围压–化学腐蚀–温度三种应力状态耦合作用后岩石的有效体积分别为 $V_0$， $V_1$， $V_2$和 $V_3$。则有：

$D_T=1-\frac{V_1}{V_0}$(17)

$D_C=1-\frac{V_2}{V_1}$(18)

$D_M=1-\frac{V_3}{V_2}$(19)

联立式可得围压–化学腐蚀–温度耦合作用下岩石的总损伤变量：

$D_H=D_T+D_C+D_M-D_T{D}_C-D_T{D}_M-D_C{D}_M+D_T{D}_C{D}_M$(20)

式中， $D_T$为温度损伤变量， $D_C$为化学损伤变量， $D_M$为力损伤变量， $D_H$为围压–化学腐蚀–温度耦合作用下岩石的总损伤变量。

真三轴压缩下的岩石损伤本构关系为：

${\sigma }_1=E_{CT}{\epsilon }_1\left(1-D_T\right)\left(1-D_C\right)\left(1-D_M\right)+\nu \left({\sigma }_2+{\sigma }_3\right)$(21)

假三轴压缩下的岩石损伤本构关系为：

${\sigma }_1=E_{CT}{\epsilon }_1\left(1-D_T\right)\left(1-D_C\right)\left(1-D_M\right)+2\nu {\sigma }_3$(22)

式中， $E_T$为高温处理后岩石的弹性模量， $\nu$为岩石的泊松比。

根据有效应力概念和Lemaitre应变等价性原理，名义应力下含损伤岩石产生的应变等价于有效应力下无损伤岩石产生的有效应变(其中名义应力即为试验测得的应力)。则可建立岩石损伤本构关系如下：

${\sigma }_i^{*}=\frac{{\sigma }_i}{1-D};i=1,2,3$(23)

式中， ${\sigma }_i^{*}$为岩石的有效应力， ${\sigma }_i$为岩石的名义应力， $D$为岩石的损伤变量。

在对损伤变量修正以后，有效应力 ${\sigma }_i^{*}$和名义应力 ${\sigma }_i$关系式变为：

${\sigma }_i^{*}=\frac{{\sigma }_i}{1-\lambda D}=\frac{{\sigma }_i}{1-D\sqrt{\frac{{\sigma }_r}{{\sigma }_p}}};i=1,2,3$(24)

由于修正系数只是用于修正力损伤变量，对温度损伤变量和化学损伤变量没有影响，则修正后的围压–化学腐蚀–温度耦合作用下岩石的总损伤变量耦合公式变为：

$D_H=D_T+D_C+\lambda D_M-D_T{D}_C-D_T\lambda D_M-D_C\lambda D_M+D_T{D}_C\lambda D_M$ (25)

将公式(4)、(15)和(16)代入式(25)，岩石的损伤总变量为：

$D_H=1-\frac{E_T}{E_0}\left[1-{\left[\frac{A}{{\left({\omega }^d\right)}_0}\right]}^{\frac{2}{3}}\right]\exp \left[-{\left(\frac{F}{F_0}\right)}^m\right]$ (26)

其中

$A={\displaystyle \underset{i=1}{\overset{n}{\sum }}\frac{m_i{\alpha }_i}{n_i}\left({10}^{-{\text{pH}}_0}-{10}^{-\text{pH}\left(t\right)}\right)V_{cs}}+{\displaystyle \underset{i=1}{\overset{m}{\sum }}\frac{m_{Di}}{M_i}}$(27)

岩石的三轴压缩试验，由于不是轴向压力和围压同时加载，而是先施加围压，待围压达到预定值并稳定后再施加轴向压力。因此在施加轴向压力之前，岩石在围压作用下已经产生了初始轴向应变，所以实际的轴向应变应为试验测量应变值与初始轴向应变之和。

根据Lemaitre应变等价性原理和有效应力的概念，名义应力下含损伤岩石产生的应变等价于有效应力下无损伤岩石产生的有效应变(其中名义应力即为试验测得的应力)。则可建立岩石损伤本构关系如下

${\sigma }_i^{*}=\frac{{\sigma }_i}{1-D_S};i=1,2,3$ (28)

式中， ${\sigma }_i^{*}$为岩石的有效应力， ${\sigma }_i^{}$为岩石的名义应力， $D_S$为岩石的损伤变量。则根据广义胡克定律可知

${\epsilon }_i^{*}=\left[{\sigma }_i^{*}-v\left({\sigma }_j^{*}+{\sigma }_k^{*}\right)\right]/E_S$(29)

式中： $\left(i,j,k\right)$为 $\left(1,2,3\right)$， $E_S$为除力场外的所有场叠加作用下岩石的弹性模量，v为泊松比， ${\epsilon }_i^{*}$为与有效应力 ${\sigma }_i^{*}$所对应的有效应变。

根据变形协调条件得

${\epsilon }_i^{*}={\epsilon }_i$(30)

由式(20)，(21)和(22)可得岩石损伤本构关系

${\sigma }_1=E_S{\epsilon }_1\left(1-D_S\right)+v\left({\sigma }_2+{\sigma }_3\right)$(31)

由岩石的三轴压缩试验可知，三轴压缩试验的加载顺序是通过液压油先施加围压，待围压稳定至预定值时，再加载轴向压力。因此在施加轴向压力之前，岩石会在围压的作用下产生初始轴向应变，所以试验中记录的轴向偏应力 ${\sigma }_{1t}$实际上为轴向应力 ${\sigma }_1$和围压 ${\sigma }_3$的差值，即

${\sigma }_{1t}={\sigma }_1-{\sigma }_3$(32)

在围压作用下，岩石产生的初始轴向应变为

${\epsilon }_{10}=\frac{1-2v}{E_T}\sigma {}_3$(33)

由于试验结果中不包含轴向压力施加前产生的初始轴向应变，故实际轴向应变 ${\epsilon }_1$为试验中测量的轴向应变 ${\epsilon }_{1t}$与初始轴向应变 ${\epsilon }_{10}$之和，即

${\epsilon }_1={\epsilon }_{1t}+{\epsilon }_{10}$(34)

修正后的化学和围压耦合作用下岩石损伤本构关系为：

$\begin{array}{c}{\sigma }_{1t}=\left[E_{CT}{\epsilon }_{1t}+\left(1-2v\right){\sigma }_3\right]\left(1-D_H\right)+\left(2\nu -1\right){\sigma }_3\\ =\left[E_{CT}{\epsilon }_{1t}+\left(1-2v\right){\sigma }_3\right]\left(1-D_M\right)\left(1-D_C\right)\left(1-D_T\right)+\left(2\nu -1\right){\sigma }_3\end{array}$ (35)

将式(26)代入(35)得：

${\sigma }_{1t}\text{=}\left[E_{CT}{\epsilon }_{1t}+\left(1-2v\right){\sigma }_3\right]\frac{E_T}{E_0}\left[1-{\left[\frac{A}{{\left({\omega }^d\right)}_0}\right]}^{\frac{2}{3}}\right]\exp \left[-{\left(\frac{F}{F_0}\right)}^m\right]+\left(2\nu -1\right){\sigma }_3$ (36)

其中

$A={\displaystyle \underset{i=1}{\overset{n}{\sum }}\frac{m_i{\alpha }_i}{n_i}\left({10}^{-{\text{pH}}_0}-{10}^{-\text{pH}\left(t\right)}\right)V_{cs}}+{\displaystyle \underset{i=1}{\overset{m}{\sum }}\frac{m_{Di}}{M_i}}$(37)

$F=\frac{E_{CT}{\epsilon }_{1t}+{\sigma }_3+c\cot \phi }{{\sigma }_3+c\cot \phi }$(38)

从上述模型可以看出，该模型中需要确定的参数有 $E_T$和 $m,F_0$，对于参数 $E_T$可采用2次多项式拟合弹性模量随温度变化的趋势来确定，而对于参数 $m,F_0$可采用岩石应力应变曲线的峰值点 $\left({\epsilon }_p,{\sigma }_p\right)$来确定，即极值法。

在上文建立的损伤演化方程和本构关系中，包含2个Weibull分布参数 $m$和 $F_0$作为损伤模型的待求参量，本文采用岩石应力应变曲线的峰值点 $\left({\epsilon }_p,{\sigma }_p\right)$来确定模型的参数，即极值法，具体求解过程如下。

峰值应力 ${\sigma }_p$与所对应的峰值应变 ${\epsilon }_p$满足下列两个几何条件

$\epsilon ={\epsilon }_p,{\sigma }_1={\sigma }_p$ (39)

$\epsilon ={\epsilon }_p,\frac{\text{d}{\sigma }_1}{\text{d}{\epsilon }_1}=0$ (40)

将式(39)和(40)代入(36)中，可得到两个关于Weibull分布参数 $m$和 $F_0$的关系式，并最终求出模型参数 $m$和 $F_0$。

首先将式(39)代入式(36)中可得

$\frac{E_T}{E_0}\left[1-{\left[\frac{A}{{\left({\omega }^d\right)}_0}\right]}^{\frac{2}{3}}\right]\exp \left[-{\left(\frac{F}{F_0}\right)}^m\right]=\frac{{\sigma }_p+\left(1-2v\right){\sigma }_3}{E_{CT}{\epsilon }_p+\left(1-2v\right){\sigma }_3}$ (41)

其中

$F_{sc}=\frac{E_{CT}{\epsilon }_{1t}+{\sigma }_3+c\cot \phi }{{\sigma }_3+c\cot \phi }$ (42)

$A={\displaystyle \underset{i=1}{\overset{n}{\sum }}\frac{m_i{\alpha }_i}{n_i}\left({10}^{-{\text{pH}}_0}-{10}^{-\text{pH}\left(t\right)}\right)V_{cs}}+{\displaystyle \underset{i=1}{\overset{m}{\sum }}\frac{m_{Di}}{M_i}}$ (43)

对式(36)取偏导得：

${\frac{\partial {\sigma }_{1t}}{\partial {\epsilon }_{1t}}|}_{\begin{array}{l}{\sigma }_{1t}={\sigma }_p\\ {\epsilon }_{1t}={\epsilon }_p\end{array}}=\frac{E_T}{E_0}\left[1-{\left[\frac{A}{{\left({\omega }^d\right)}_0}\right]}^{\frac{2}{3}}\right]\exp \left[-{\left(\frac{F}{F_0}\right)}^m\right]\left\{E_{CT}-\frac{m\left[E_{CT}{\epsilon }_p+\left(1-2v\right){\sigma }_3\right]}{F_{sc}}{\left(\frac{F_{sc}}{F_0}\right)}^m\frac{\partial F_{sc}}{\partial {\epsilon }_{1t}}\right\}$ (44)

将式(40)代入(44)得：

$E_{CT}-\frac{m\left[E_{CT}{\epsilon }_p+\left(1-2v\right){\sigma }_3\right]}{F_{sc}}{\left(\frac{F_{sc}}{F_0}\right)}^m\frac{\partial F_{sc}}{\partial {\epsilon }_{1t}}=0$ (45)

对式(42)求偏导得：

${\frac{\partial F_{sc}}{\partial {\epsilon }_{1t}}|}_{\begin{array}{l}{\sigma }_{1t}={\sigma }_p\\ {\epsilon }_{1t}={\epsilon }_p\end{array}}=\frac{E_{CT}}{{\sigma }_3+c\cot \phi }$ (46)

将式(46)代入式(45)得：

${\left(\frac{F_{\text{sc}}}{F_0}\right)}^m=\frac{E_{CT}{\epsilon }_p+{\sigma }_3+c\cot \phi }{m\left[E_{CT}{\epsilon }_p+\left(1-2v\right){\sigma }_3\right]}$ (47)

联立式(41)和(47)可得

$m=\frac{1}{\ln \left(B\right)}\frac{E_{CT}{\epsilon }_p+{\sigma }_3+c\cot \phi }{\left[E_{CT}{\epsilon }_p+\left(1-2v\right){\sigma }_3\right]}$ (48)

$F_0=\frac{E_{CT}{\epsilon }_p+{\sigma }_3+c\cot \phi }{{\sigma }_3+c\cot \phi }{\left(\frac{m\left[E_{CT}{\epsilon }_p+\left(1-2v\right){\sigma }_3\right]}{E_{CT}{\epsilon }_p+{\sigma }_3+c\cot \phi }\right)}^{\frac{1}{m}}$ (49)

其中

$B=\frac{E_{CT}{\epsilon }_p+\left(1-2v\right){\sigma }_3}{{\sigma }_p+\left(1-2v\right){\sigma }_3}\frac{E_T}{E_0}\left[1-{\left[\frac{A}{{\left({\omega }^d\right)}_0}\right]}^{\frac{2}{3}}\right]$(50)

为了验证本文所建立的岩石损伤演化方程、统计损伤本构模型和参数确定方法的合理性及准确性，本文通过引用Chen [29] 的相关研究成果来验证。

在该试验过程中，水化学处理组砂岩试样分别经过pH为1，3，5的硝酸溶液处理，并设置了经过pH为7的蒸馏水处理的对照组。高温处理组加热升温速率为10℃/min，恒温处理4小时，处理温度分别为25℃、300℃、600℃、900℃。以及高温和水化学耦合处理组，先经过水化学处理再高温处理。

该实验中选取砂岩试样的元素含量如 表1 所示。

温度–水化学–围压耦合作用涉及岩石中氧化物与水溶液反应，温度–水化学–围压耦合作用涉及岩石中氧化物与水溶液反应，砂岩中粘土矿物与酸性溶液之间的化学反应方程主要包括：

(1) 氧化物与酸反应：

$\{\begin{array}{l}{\text{NaAlSi}}_3{\text{O}}_8+4{\text{H}}^{+}+4{\text{H}}_{\text{2}}\text{O}\to 3{\text{H}}_{\text{4}}{\text{SiO}}_{\text{4}}+{\text{Na}}^{+}+{\text{Al}}^{3+}\\ {\text{KAlSi}}_{\text{3}}{\text{O}}_{\text{8}}+4{\text{H}}^{+}+4{\text{H}}_{\text{2}}\text{O}\to 3{\text{H}}_{\text{4}}{\text{SiO}}_{\text{4}}+{\text{K}}^{+}+{\text{Al}}^{3+}\\ {\text{CaAl}}_{\text{2}}{\text{Si}}_{\text{2}}{\text{O}}_{\text{8}}+8{\text{H}}^{+}\to 2{\text{H}}_{\text{4}}{\text{SiO}}_{\text{4}}+{\text{Ca}}^{2+}+2{\text{Al}}^{3+}\\ {\text{KAl}}_{\text{3}}{\text{Si}}_{\text{3}}{\text{O}}_{\text{10}}{\left(\text{OH}\right)}_{\text{2}}+10{\text{H}}^{+}\to 3{\text{H}}_{\text{4}}{\text{SiO}}_{\text{4}}+{\text{K}}^{+}+3{\text{Al}}^{3+}\\ {\text{CaCO}}_3+2{\text{H}}^{+}\to {\text{Ca}}^{2+}+{\text{H}}_{\text{2}}\text{O}+{\text{CO}}_2\uparrow \\ \text{MgO}+2{\text{H}}^{+}\to {\text{Mg}}^{2+}+{\text{H}}_{\text{2}}\text{O}\\ {\text{Al}}_{\text{2}}{\text{O}}_3+6{\text{H}}^{+}\to 2{\text{Al}}^{3+}+3{\text{H}}_{\text{2}}\text{O}\\ {\text{Fe}}_{\text{2}}{\text{O}}_3+6{\text{H}}^{+}\to 2{\text{Fe}}^{3+}+3{\text{H}}_{\text{2}}\text{O}\\ \text{CaO}+2{\text{H}}^{+}\to {\text{Ca}}^{2+}+{\text{H}}_{\text{2}}\text{O}\\ {\text{Na}}_{\text{2}}\text{O}+2{\text{H}}^{+}\to 2{\text{Na}}^{+}+{\text{H}}_{\text{2}}\text{O}\end{array}$ (51)

(2) 氧化物与水的反应：

$\{\begin{array}{l}\text{CaO}+{\text{H}}_{\text{2}}\text{O}\to \text{Ca}{\left(\text{OH}\right)}_2\\ \text{MgO}+{\text{H}}_{\text{2}}\text{O}\to \text{Mg}{\left(\text{OH}\right)}_2\end{array}$ (52)

在酸性环境中，砂岩中的物质与 ${\text{H}}^{+}$反应后只能置换出金属离子，不能将 $\text{Si}$变成离子融入水溶液使之脱离砂岩，因此，通过计算反应物和所涉及的产物的摩尔量，就可以确定已经丢失的溶解混凝土的摩尔量。

根据公式(12)和(13)，计算得出不同pH值水化学环境下砂岩可溶性胶结材料反应量，见 表1 。再根据公式(48)和(49)计算得出不同条件下威布尔分布系数 $m$和 $F_0$，见 表2 。

续表