本发明公开一种低渗透多孔介质材料的流固耦合模型构建方法，包括如下步骤：建立低渗透岩石渗流场对应力场影响的数学模型；建立低渗透岩石孔隙度与绝对渗透率的关系，在此基础上建立低渗透岩石应力场对渗流场影响的数学模型；将低渗透岩石渗流场对应力场影响的数学模型和应力场对渗流场影响的数学模型进行联立，建立低渗透岩石流固耦合关系。本发明在研究应力场对渗流场影响的数学模型时，采用气体作为渗流介质测得了低渗透岩石的绝对渗透率，并采用函数拟合消除了低渗透岩石气体滑脱效应的影响，使得建立的应力场对渗流场影响的数学模型较为精确，最终建立的低渗透岩石流固耦合关系准确可靠；而且，该方法不必进行原位试验，节省了物力和财力。
1.一种低渗透多孔介质材料的流固耦合模型构建方法，其特征在于，包括如下步骤：步骤1，建立低渗透岩石渗流场对应力场影响的数学模型；
式(1)中，f为渗流体积力，rw是水的比重，H为水头，fx和fy是渗流体积力在x和y方向的分量，Jx和Jy是水力梯度在x和y方向的分量；
其次将单元的渗流体积力转换为等效节点荷载，如下式(2)和式(3)：式(2)中，{Fs}为渗流体积力产生的等效节点力，[N]为单元形状函数，Ω为已知力的边界条件；式(3)中，{ΔFs}为渗流体积力增量产生的等效节点力；
最后建立低渗透岩石渗流场对应力场影响的数学模型，如下式(4)：式中，kliq为低渗透岩石的绝对渗透率，{ΔF}为节点荷载增量，{δ}为位移，{δ0}为初始位移，{Δδ}为位移增量，σij为应力，nj为等效节点，Su为已知的位移边界条件，Sa为确定的应力边界条件；
步骤2，建立低渗透岩石孔隙度与绝对渗透率的关系，在此基础上建立低渗透岩石应力场对渗流场影响的数学模型；所述低渗透岩石应力场对渗流场影响的数学模型如下式(5)：式(5)中，φ为孔隙度，kliq(φ)为低渗透岩石绝对渗透率与孔隙度的函数关系，h为试样高度，H为水头，q1为流量；Γ1为有确定水头的边界条件，Γ2为有确定流量的边界条件，n2为Γ2边界条件边界的法线；Γ3为混合边界条件，即 α1和β1是常量，n3为Γ3边界条件边界的法线；
步骤3，将低渗透岩石渗流场对应力场影响的数学模型和应力场对渗流场影响的数学模型进行联立，建立低渗透岩石流固耦合关系。
2.根据权利要求1所述的低渗透多孔介质材料的流固耦合模型构建方法，其特征在于，步骤2中，所述建立低渗透岩石孔隙度与绝对渗透率的关系包括下述步骤：步骤21，测量低渗透岩石在不同围压下的渗透率、确定相应的绝对渗透率；
步骤23，分别确定低渗透岩石绝对渗透率与围压的函数关系、以及低渗透岩石孔隙度与围压的函数关系；
步骤24，通过公式推导建立低渗透岩石孔隙度与绝对渗透率的关系。
3.根据权利要求2所述的低渗透多孔介质材料的流固耦合模型构建方法，其特征在于，步骤21中，采用下式(6)对不同围压下低渗透岩石气测渗透率与孔隙压力倒数进行拟合，得到不同围压下低渗透岩石消除滑脱效应的绝对渗透率：式(6)中，kliq为绝对渗透率，k为气测渗透率，b为滑脱因子，a为二次滑脱因子，pc为孔隙压力。
4.根据权利要求3所述的低渗透多孔介质材料的流固耦合模型构建方法，其特征在于，步骤23中，所述低渗透岩石绝对渗透率与围压的函数关系如下式(7)：式(7)中，d1和h1为常数，p为围压。
5.根据权利要求4所述的低渗透多孔介质材料的流固耦合模型构建方法，其特征在于，步骤23中，所述低渗透岩石孔隙度与围压的函数关系如下式(8)：φ＝φ0e-ip (8)；
6.根据权利要求5所述的低渗透多孔介质材料的流固耦合模型构建方法，其特征在于，步骤24中，结合公式(7)和公式(8)，确定所述低渗透岩石孔隙度与绝对渗透率的关系，如下式(10)：
7.根据权利要求1所述的低渗透多孔介质材料的流固耦合模型构建方法，其特征在于，步骤3中，最终建立的低渗透岩石的流固耦合关系如下式(13)：式(13)中，{f}为渗流场水头分布函数。
一种低渗透多孔介质材料的流固耦合模型构建方法
技术领域
[0001] 本发明涉及一种流固耦合模型构建方法，特别涉及一种低渗透多孔介质材料的流固耦合模型构建方法。
背景技术
[0002] 低渗透岩石是许多土木工程项目的主要介质，如煤矿开采，CO2储存，石油燃料开采和核废料处理。低渗透岩石的力学性质对工程安全有着显著的影响，许多研究人员已经做了大量关于低渗透岩石力学性质的研究，这些研究普遍没有考虑液体渗流对低渗透岩石力学性质的影响；但是考虑流固耦合的低渗透岩石的力学响应与不考虑渗流的低渗透岩石的力学响应有较大不同。因此研究考虑流固耦合情况下的低渗透岩石的力学行为较为迫切。
[0003] 构建低渗透岩石解析模型是一种有效研究低渗透岩石的力学响应的方法，但是构建低渗透岩石的流固耦合解析模型是较为困难的。因为构成低渗透岩石流固耦合解析模型的关键参数-低渗透岩石的液体渗透率(绝对渗透率)是很难通过液体作为介质去测量的，而通过气体测量的低渗透岩石的绝对渗透率由于滑脱效应的存在导致测量结果偏大。因此，如何选择一种合适的方法去构建低渗透岩石的流固耦合解析模型需要进一步的研究。
发明内容
[0004] 发明目的：本发明的目的针对现有技术中存在的问题，提供一种低渗透多孔介质材料的流固耦合模型构建方法。
[0005] 技术方案：本发明所述的一种低渗透多孔介质材料的流固耦合模型构建方法，包括如下步骤：
[0006] 步骤1，建立低渗透岩石渗流场对应力场影响的数学模型；
[0007] 步骤2，建立低渗透岩石孔隙度与绝对渗透率的关系，在此基础上建立低渗透岩石应力场对渗流场影响的数学模型；
[0008] 步骤3，将低渗透岩石渗流场对应力场影响的数学模型和应力场对渗流场影响的数学模型进行联立，建立低渗透岩石流固耦合关系。
[0009] 上述步骤1中，低渗透岩石的渗流场是通过改变低渗透岩石的渗流体积力去影响应力场的分布，故低渗透岩石渗流场对应力场影响的数学模型通过下述步骤建立：
[0010] 在二维渗流环境中，低渗透岩石的渗流体积力与水力梯度成正比，两者关系如下式(1)：
[0011]
[0012] 式(1)中，f为渗流体积力，rw是水的比重，H为水头，fx和fy是渗流体积力在x和y方向的分量，Jx和Jy是水力梯度在x和y方向的分量；
[0013] 其次将单元的渗流体积力转换为等效节点荷载，如下式(2)和式(3)：
[0014]
[0015]
[0016] 式(2)中，{Fs}为渗流体积力产生的等效节点力，[N]为单元形状函数，Ω为已知力的边界条件；式(3)中，{ΔFs}为渗流体积力增量产生的等效节点力；
[0017] 在此基础上，低渗透岩石渗流场对应力场影响的数学模型得到了，如下式(4)：
[0018]
[0019] 式中，kliq为低渗透岩石的绝对渗透率，{ΔF}为节点荷载增量，{δ}为位移，{δ0}为初始位移，{Δδ}为位移增量，σij为应力，nj为等效节点，Su为已知的位移边界条件，Sa为确定的应力边界条件。
[0020] 低渗透岩石应力场对渗流场的影响主要通过影响低渗透岩石的孔隙度从而影响其渗透率来实现。因此，步骤2中，考虑应力场影响的渗流场数学模型如公式(5)：
[0021]
[0022] 式(5)中，φ为孔隙度，kliq(φ)为低渗透岩石绝对渗透率与孔隙度的函数关系，h为试样高度，H为水头，q1为流量；Γ1为有确定水头的边界条件，Γ2为有确定流量的边界条件，n2为Γ2边界条件边界的法线；Γ3为混合边界条件，即 α1和β1是常量，n3为Γ3边界条件边界的法线。
[0023] 因此，确定应力场对渗流场影响的关键是确定低渗透岩石渗透率与孔隙度之间的关系。步骤2中，建立低渗透岩石孔隙度与绝对渗透率的关系包括下述步骤：
[0024] 步骤21，测量低渗透岩石在不同围压下的渗透率、确定相应的绝对渗透率；
[0025] 步骤22，测量低渗透岩石在对应围压下的孔隙度；
[0026] 步骤23，分别确定低渗透岩石绝对渗透率与围压的函数关系、以及低渗透岩石孔隙度与围压的函数关系；
[0027] 步骤24，通过公式推导建立低渗透岩石孔隙度与绝对渗透率的关系。
[0028] 具体的，步骤21中，采用下式(6)对不同围压下低渗透岩石气测渗透率与孔隙压力倒数进行拟合，得到不同围压下低渗透岩石消除滑脱效应的绝对渗透率：
[0029]
[0030] 式(6)中，kliq为绝对渗透率，k为气测渗透率，b为滑脱因子，a为二次滑脱因子，pc为孔隙压力。
[0031] 由此，步骤23中，低渗透岩石绝对渗透率与围压的函数关系如下式(7)：
[0032]
[0033] 式(7)中，d1和h1为常数，p为围压。
[0034] 上述步骤23中，低渗透岩石孔隙度与围压的函数关系如下式(8)：
[0035] φ＝φ0e-ip (8)；
[0036] 式(8)中，φ0为初始孔隙度，i为常数。
[0037] 步骤24中，结合公式(7)和公式(8)，确定低渗透岩石孔隙度与绝对渗透率的关系，如下式(10)：
[0038]
[0039] 最终步骤3中，将步骤1和步骤2所得数学模型联系，建立的低渗透岩石的流固耦合关系如下式(13)：
[0040]
[0041] 式(13)中，{f}为渗流场水头分布函数。
[0042] 有益效果：与现有技术相比，本发明的显著优点在于：(1)本发明通过将低渗透岩石渗流场对应力场影响的数学模型和应力场对渗流场影响的数学模型联立、得到低渗透岩石的流固耦合关系，其中，研究应力场对渗流场影响的数学模型时，采用气体作为渗流介质测得了低渗透岩石的绝对渗透率，并采用函数拟合消除了低渗透岩石气体滑脱效应的影响，使得建立的应力场对渗流场影响的数学模型较为精确，最终建立的低渗透岩石流固耦合关系准确可靠；(2)本发明通过室内试验和数学公式推导来确定低渗透岩石流固耦合的关系，不必进行原位试验，方法简便易行，节省了大量物力和财力。
具体实施方式
[0043] 下面对本发明的技术方案作进一步说明。
[0044] 以确定某地下水封洞库低渗透岩石流固耦合关系为例，对本发明的一种低渗透多孔介质材料的流固耦合模型构建方法进行说明。
[0045] 步骤1建立低渗透岩石渗流场对应力场影响的数学模型
[0046] 低渗透岩石的渗流场是通过改变低渗透岩石的渗流体积力去影响应力场的分布，在二维渗流环境中，低渗透岩石的渗流体积力与水力梯度成正比，如公式(1)所示：
[0047]
[0048] 式中，f为渗流体积力，rw是水的比重，H为水头，fx和fy是渗流体积力在x和y方向的分量，Jx和Jy是水力梯度在x和y方向的分量；
[0049] 以下公式可以转换单元的渗流体积力到等效节点荷载，如公式(2)和公式(3)所示：
[0050]
[0051] 式中，{Fs}为渗流体积力产生的等效节点力，[N]是单元形状函数，Ω为已知力的边界条件。
[0052]
[0053] 式中，{ΔFs}为渗流体积力增量产生的等效节点力。
[0054] 因此考虑渗流场影响的应力场的数学模型得到了，如公式(4)所示。
[0055]
[0056] 式中，kliq为低渗透岩石的绝对渗透率，{ΔF}为节点荷载增量，{δ}为位移，{δ0}为初始位移，{Δδ}为位移增量，σij为应力，nj为等效节点，Su为已知的位移边界条件，Sa为确定的应力边界条件。
[0057] 步骤2：所述建立低渗透岩石应力场对渗流场影响的数学模型：
[0058] 应力场对低渗透岩石渗流场的影响主要通过影响低渗透岩石的孔隙度从而影响其渗透率实现。因此，考虑应力场影响的渗流场二维模型如公式(5)所示：
[0059]
[0060] 式中，φ为孔隙度，kliq(φ)为低渗透岩石绝对渗透率与孔隙度的函数关系，h为试样高度，H为水头，q1为流量；Γ1为有确定水头的边界条件，Γ2为有确定流量的边界条件，n2为Γ2边界条件边界的法线；Γ3为混合边界条件，即 α1和β1是常量，n3为Γ3边界条件边界的法线。
[0061] 因此，确定应力场对渗流场影响的关键是确定低渗透岩石渗透率与孔隙度之间的关系。
[0062] 步骤3：确定低渗透岩石渗透率与围压的函数关系
[0063] 采用公式(6)对不同围压下低渗透岩石气测渗透率与孔隙压力倒数进行拟合，得到不同围压下低渗透岩石消除滑脱效应的绝对渗透率。
[0064]
[0065] 式中，kliq为绝对渗透率，k为气测渗透率，b为滑脱因子，a为二次滑脱因子，pc为孔隙压力。
[0066] 通过拟合得到的低渗透岩石绝对渗透率如表1。
[0067] 表1低渗透岩石绝对渗透率
[0068]
[0069] 通过函数拟合，得到低渗透岩石绝对渗透率与围压成指数函数关系如式(7)，拟合参数如表2。
[0070]
[0071] 式中d1和h1是常数，p为围压。
[0072] 表2拟合参数
[0073]
[0074] 步骤4：确定低渗透岩石孔隙度与围压的函数关系
[0075] 通过函数拟合低渗透岩石孔隙度与围压成指数函数关系，如式(8)；拟合参数i为0.01。
[0076] φ＝φ0e-ip (8)
[0077] 式中，φ0是初始孔隙度，i是常数。
[0078] 步骤5：通过公式推导建立低渗透岩石孔隙度与渗透率的关系。
[0079] 结合公式(7)和公式(8)，得到公式(9)。
[0080]
[0081] 化简之后，低渗透岩石绝对渗透率与孔隙度关系曲线如公式(10)。
[0082]
[0083] 因此，低渗透岩石绝对渗透率与孔隙度成幂函数关系，将参数的值代入式(10)、式(11)得到。
[0084] kliq＝0.313φ2 (11)
[0085] 步骤6：将低渗透岩石渗流场对应力场影响的数学模型和应力场对渗流场影响的数学模型进行联立，建立低渗透岩石流固耦合关系。
[0086] 公式(11)被代入公式(5)，公式(12)得到。
[0087]
[0088] 因此低渗透岩石流固耦合的关系式得到通过结合公式(12)和公式(4)，得到公式(13)。
[0089]
[0090] 式中，{f}是渗流场水头分布函数。
