确定裂隙岩体渗透率及孔隙度对有效应力敏感程度的方法
【技术领域】
[0001 ]本发明涉及一种确定裂隙岩体对有效应力敏感程度的方法，具体涉及一种确定裂 隙岩体渗透率及孔隙度对有效应力敏感程度的方法。
【背景技术】
[0002] 经历了漫长的地质成生改造作用和表生演化，赋存于一定地质应力环境中的岩体 内部都广泛发育着节理、层理、接触带、剪切带、断层等结构面，因此在实际工程中遇到的岩 体大多都是由结构面和完整岩块组成的裂隙岩体，由于裂隙岩体中复杂的不连续的结构面 对岩体的渗透性质有控制作用，使得裂隙岩体的渗透性质比一般的均质、连续、各向同性材 料的渗透性质复杂的多，对工程安全有着显著的影响。
[0003] 在实际环境中影响岩体渗透率和孔隙结构的因素有很多，如岩体所处的应力环 境，岩体的埋深，岩体所处渗流场的影响等，因此很难对各因素的影响关系逐一描述清楚， 于是引入渗透率和孔隙度对有效应力的敏感程度，对影响岩体渗透率和孔隙度的因素进行 归一化处理，通过这一系数可以清晰的得到有效应力对渗透率和孔隙度的影响关系，敏感 程度值越大，表示岩体渗透率和孔隙度对有效应力的变化就越敏感，反之则敏感程度不高。
[0004] 现有关于岩体渗透率和孔隙度对应力敏感程度的研究如下：
[0005] 专利号为104977195A的中国发明专利"采用人造裂缝性碳酸盐岩岩心进行应力敏 感实验的方法"公开了一种测试不同充填状况裂隙岩体渗透率对应力敏感程度的方法，该 方法只能确定含充填的裂隙岩体的渗透率对应力的敏感程度，对于无充填的裂隙岩体不同 方向渗透率对应力的敏感程度没有涉及；
[0006] 专利号为104405374A的中国发明专利"一种致密气藏储层应力敏感程度的测量方 法"公开了测量一种致密气藏储层应力敏感程度的方法，该方法通过数值分析计算岩体的 渗透率对应力的敏感程度，可靠性较差；
[0007] 专利号为101038251的中国发明专利"全自动敏感综合测试仪"公布了一种自动测 试各向同性岩体应力敏感程度的方法，该方法仅能测量各向同性岩体渗透率对应力的敏感 程度，不能测量裂隙岩体不同方向渗透率对应力的敏感程度；
[0008] 专利号为104834807A的中国发明专利"一种基于分形理论的应力敏感储层相对渗 透率计算方法"公开了一种基于分形理论计算储层渗透率对应力敏感程度的方法，该方法 采用理论分析求得岩体渗透率对应力的敏感程度，采用了太多的假设，结果可靠性差；
[0009] 专利号为103528934A的中国发明专利"一种测量超低渗岩石渗透率应力敏感程度 的互相关技术"公开了一种测量超低渗透完整岩心渗透率对应力敏感程度的方法"，该方法 仅能测量超低渗透率的完整岩石渗透率对应力的敏感程度，不能测量渗透率较高的裂隙岩 体渗透率对应力的敏感程度；
[0010]专利号为102721636A的中国发明专利"一种基于尺度模型确定裂缝性低渗储层应 力敏感程度的方法"公开了一种考虑了尺度效应的完整岩体渗透率对应力敏感程度的测试 方法，该方法采用数值模拟的方法求得岩体渗透率对应力的敏感程度，采用了太多的假设， 结果可靠性差。
[0011]专利号为104089823A的中国专利申请"一种基于孔隙压缩实验确定岩石有效应力 系数的方法"，该方法通过液体在应力作用下的压缩测量岩石的孔隙度，在测量过程中液体 难免会在岩石孔隙结构中有残留，得到的岩石孔隙度对应力的敏感性不够精确。
[0012] 专利号为104133050A的中国发明专利"应力条件下多孔岩石有效应力系数与孔隙 率同时测试方法"，该方法利用在应力加载过程中的岩石应变数据、应力数据以及孔隙体积 变化数据通过理论推导得到岩石的孔隙率参数，运用了较多假设，结果难免会有误差，得到 的岩石孔隙度对应力的敏感性不够精确。
[0013] 综上所述，目前关于测试岩体渗透率和孔隙度对应力敏感程度的方法多是针对各 向同性岩体，对于各向异性岩体渗透率对有效应力敏感程度的测试方法未见有报道，但实 际工程中遇到的岩体内部大都广泛发育着节理、层理、接触带、剪切带、断层等结构面，岩体 表现出显著的渗透各向异性，其不同方向渗透率的大小与岩体所受的有效应力关系紧密， 对工程安全有着显著的影响，且目前多采用数值分析和理论推导的方法确定岩体渗透率和 孔隙度对应力的敏感程度，采用过多假设，可靠性较差，结果精确性不足。因此迫切需要一 种较为完备的确定渗透各向异性岩体不同方向渗透率和孔隙度对有效应力敏感程度的方 法。
【发明内容】
[0014] 发明目的:本发明的第一目的是提供一种确定裂隙岩体渗透率对有效应力敏感程 度的方法，本发明的第二目的是提供一种确定裂隙岩体孔隙度对有效应力敏感程度的方 法。
[0015] 技术方案:本发明涉及一种确定裂隙岩体渗透率对有效应力敏感程度的方法，包 括如下步骤：
[0016] 步骤1，模拟实际裂隙岩体，制备与水平方向呈不同倾角的相似材料；
[0017] 步骤2,选任一倾角相似材料逐级施加围压，每级围压稳定后测量其渗透率，确定 该倾角方向相似材料的渗透率与围压的函数关系式；
[0018] 步骤3:对该函数关系式求导，所得导函数表示该倾角方向相似材料的渗透率对有 效应力的变化率；
[0019] 步骤4,对该导函数进行标准化处理，得到该倾角方向单位应力变化引起的单位渗 透率的变化量，该变化量即为待测实际裂隙岩体该倾角方向的渗透率对有效应力变化的敏 感程度；
[0020] 步骤5,重复步骤2~4,得到待测实际裂隙岩体不同倾角方向的渗透率对有效应力 变化的敏感程度。
[0021] 通过本发明的方法得到不同倾角方向的渗透率对有效应力变化的敏感程度，充分 反映了各向异性岩体渗透率对有效应力的敏感程度，并且，标准化处理后能够得到单位应 力变化引起的单位渗透率变化量，可准确表征渗透率差异较大的渗透各向异性岩体不同方 向渗透率对有效应力敏感程度。
[0022] 具体的，步骤1中，制备与水平方向呈不同倾角的相似材料的方法包括：
[0023] (1)根据实际裂隙岩体的几何形态以及该岩体中裂隙的形态和分布，制作该裂隙 岩体的三维数字模型；
[0024] (2)将该三维数字模型输入3D打印机，旋转该三维数字模型，以实际裂隙岩体的粉 末为原料打印出与水平方向呈不同倾角的相似材料。
[0025] 上述倾角为0。、15。、30。、45。、60。、75。或90。。
[0026] 上述步骤2中，根据测得的每级围压下的渗透率，以围压为横坐标、渗透率为纵坐 标，作出该倾角方向围压与渗透率的散点图，通过迭代法得出与散点图拟合度最高的函数 关系式，即为该倾角方向相似材料的渗透率与围压的函数关系式。
[0027] 较优的，步骤2中，逐级施加310^、810^、151^?、2510^、3010^及3510^的围压。
[0028]上述步骤4中，对导函数进行标准化处理的方法为:将步骤3所得的导函数除以该 倾角相似材料的初始渗透率，该初始渗透率为初始围压下测得的相似材料的渗透率。
[0029] 本发明所述的一种确定裂隙岩体孔隙度对有效应力敏感程度的方法，包括如下步 骤：
[0030] 步骤1，模拟实际裂隙岩体，制备裂隙岩体的相似材料；
[0031] 步骤2,对该相似材料逐级施加围压，每级围压稳定后测量其孔隙度，确定孔隙度 与围压的函数关系式；
[0032]步骤3:对该函数关系式求导，所得导函数表示孔隙度对有效应力的变化率；<