CN105606513B - 确定裂隙岩体渗透率及孔隙度对有效应力敏感程度的方法 - Google Patents

确定裂隙岩体渗透率及孔隙度对有效应力敏感程度的方法 Download PDF

Info

Publication number
CN105606513B
CN105606513B CN201610064446.XA CN201610064446A CN105606513B CN 105606513 B CN105606513 B CN 105606513B CN 201610064446 A CN201610064446 A CN 201610064446A CN 105606513 B CN105606513 B CN 105606513B
Authority
CN
China
Prior art keywords
permeability
porosity
effective stress
crack rock
confining pressure
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201610064446.XA
Other languages
English (en)
Other versions
CN105606513A (zh
Inventor
巢志明
王环玲
徐卫亚
吉华
贾朝军
闫龙
林志南
杨兰兰
赵恺
夏季
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Hohai University HHU
Original Assignee
Hohai University HHU
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Hohai University HHU filed Critical Hohai University HHU
Priority to CN201610064446.XA priority Critical patent/CN105606513B/zh
Publication of CN105606513A publication Critical patent/CN105606513A/zh
Application granted granted Critical
Publication of CN105606513B publication Critical patent/CN105606513B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N15/00Investigating characteristics of particles; Investigating permeability, pore-volume, or surface-area of porous materials
    • G01N15/08Investigating permeability, pore-volume, or surface area of porous materials
    • G01N15/082Investigating permeability by forcing a fluid through a sample
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N1/00Sampling; Preparing specimens for investigation
    • G01N1/28Preparing specimens for investigation including physical details of (bio-)chemical methods covered elsewhere, e.g. G01N33/50, C12Q

Abstract

本发明公开一种确定裂隙岩体渗透率及孔隙度对有效应力敏感程度的方法,其中,确定渗透率对有效应力敏感程度的方法包括:模拟实际裂隙岩体,制备不同倾角相似材料;选任一倾角相似材料逐级施加围压,每级围压稳定后测量其渗透率,确定该倾角方向的渗透率与围压的函数关系式;对函数关系式求导,所得导函数表示该倾角方向相似材料的渗透率对有效应力的变化率;对导函数进行标准化处理,得到该倾角方向单位应力变化引起的单位渗透率的变化量,即为该倾角方向的渗透率对有效应力变化的敏感程度;重复上述,得到不同倾角方向的渗透率对有效应力变化的敏感程度。本发明可准确表征渗透率差异较大的岩体不同方向渗透率对有效应力的敏感程度。

Description

确定裂隙岩体渗透率及孔隙度对有效应力敏感程度的方法
技术领域
[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:对该函数关系式求导,所得导函数表示孔隙度对有效应力的变化率;
[0033] 步骤4,对该导函数进行标准化处理,得到单位应力变化引起的单位孔隙度的变化 量,该变化量即为待测实际裂隙岩体孔隙度对有效应力变化的敏感程度。
[0034] 通过本发明的方法得到孔隙度对有效应力变化的敏感程度,标准化处理后能够得 到单位应力变化引起的单位孔隙度变化量,可准确表征孔隙度差异较大的裂隙岩体孔隙度 对有效应力敏感程度。
[0035] 具体的,步骤2中,根据测得的每级围压下的孔隙度,以围压为横坐标、孔隙度为纵 坐标,作出围压与孔隙度的散点图,通过迭代法得出与散点图拟合度最高的函数关系式,即 为孔隙度与围压的函数关系式。
[0036] 上述步骤2中,逐级施加310^、810^、151^»、2510^、3010^及3510^的围压。
[0037] 上述步骤4中,对导函数进行标准化处理的方法为:将步骤3所得的导函数除以初 始孔隙度,该初始孔隙度为初始围压下测得的相似材料的孔隙度。
[0038] 有益效果:与现有技术相比,本发明的显著优点在于:(1)本发明通过制备相似材 料来模拟实际裂隙岩体,不必进行原位试验,通过室内模型试验即可确定渗透各向异性裂 隙岩体不同方向渗透率对有效应力的敏感程度,节省了大量物力和财力;同时,以数学分析 方法定量确定裂隙岩体不同方向渗透率对有效应力的变化率,并对该变化率表达式进行标 准化处理,得到不同倾角方向上单位有效应力变化引起的单位渗透率的变化量,可以准确 表征渗透率差异较大的裂隙岩体不同方向渗透率对有效应力的敏感程度;(2)本发明以数 学分析方法定量确定裂隙岩体孔隙度对有效应力的变化率,并对该变化率表达式进行标准 化处理,得到单位有效应力变化引起的单位孔隙度的变化量,可准确表征孔隙度大小差异 较大的裂隙岩体孔隙度对有效应力敏感程度;
[0039] (3)本发明的相似材料可通过3D打印方法制备,可以精确反应各裂隙的空间形态 和分布,所得相似材料能够充分模拟实际裂隙岩体的性质,将其用于定量分析渗透各向异 性岩体的渗透率及孔隙度对有效应力的敏感系数时,实验结果具有代表性。
附图说明
[0040] 图1为实施例中测得的不同倾角方向相似材料的渗透率与围压的散点图;
[0041] 图2为实施例中测得的相似材料的孔隙度与围压的散点图。
具体实施方式
[0042] 下面结合附图对本发明的技术方案作进一步说明。
[0043] 本发明的一种确定裂隙岩体渗透率对有效应力敏感程度的方法,包括如下步骤:
[0044] 步骤1,模拟实际裂隙岩体,制备与水平方向呈不同倾角的相似材料;
[0045] 倾角可为0°、15°、30°、45°、60°、75°或90°;具体的,制备与水平方向呈不同倾角的 相似材料的方法包括如下两种:
[0046] 第一种,根据实际裂隙岩体的几何形态以及该岩体中裂隙的形态和分布,分别制 作用于制备不同倾角相似材料的模具;根据实际裂隙岩体的力学性能配置水泥砂浆,倒入 模具,养护、脱模,得到不同倾角的相似材料。
[0047] 第二种,根据实际裂隙岩体的几何形态以及该岩体中裂隙的形态和分布,制作该 裂隙岩体的三维数字模型;将该三维数字模型输入3D打印机,旋转该三维数字模型,以实际 裂隙岩体的粉末为原料打印出与水平方向呈不同倾角的相似材料。通过3D打印方法制备相 似材料,可以精确反应各裂隙的空间形态和分布,所得相似材料能够充分模拟实际裂隙岩 体的性质,将其用于定量分析渗透各向异性时,实验结果具有代表性。
[0048] 步骤2,选任一倾角相似材料逐级施加围压,每级围压稳定后测量其渗透率,确定 该倾角方向相似材料的渗透率与围压的函数关系式;
[0049] 可逐级施加310^、810^、151^»、2510^、3010^及3510^的围压。
[0050] 根据测得的每级围压下的渗透率,以围压为横坐标、渗透率为纵坐标,在坐标轴内 作出该倾角方向围压与渗透率的散点图,通过迭代法得出与散点图拟合度最高的函数关系 式,即为该倾角方向相似材料的渗透率与围压的函数关系式。
[0051] 步骤3:对该函数关系式求导,所得导函数表示该倾角方向相似材料的渗透率对有 效应力的变化率;
[0052] 步骤4,对该导函数进行标准化处理,得到该倾角方向单位应力变化引起的单位渗 透率的变化量,该变化量即为待测实际裂隙岩体该倾角方向的渗透率对有效应力变化的敏 感程度;
[0053] 初始围压下测得的相似材料的渗透率为初始渗透率,将导函数除以该初始渗透 率,作为裂隙岩体渗透率对有效应力敏感程度的定量表达式,该表达式表示单位应力变化 引起的单位渗透率的变化量,表征了渗透率对有效应力变化的敏感程度。
[0054] 步骤5,重复步骤2〜4,得到待测实际裂隙岩体不同倾角方向的渗透率对有效应力 变化的敏感程度。
[0055] 通过制备相似材料来模拟实际裂隙岩体,不必进行原位试验,通过室内模型试验 即可确定渗透各向异性裂隙岩体不同方向渗透率对有效应力的敏感程度,节省了大量物力 和财力;同时,以数学分析方法定量确定裂隙岩体不同方向渗透率对有效应力的变化率,并 对该变化率表达式进行标准化处理,得到不同倾角方向上单位有效应力变化引起的单位渗 透率的变化量,该变化量可以准确表征渗透率差异较大的裂隙岩体不同方向渗透率对有效 应力的敏感程度;
[0056] 本发明的一种确定裂隙岩体孔隙度对有效应力敏感程度的方法,包括如下步骤:
[0057] 步骤1,模拟实际裂隙岩体,制备裂隙岩体的相似材料;
[0058] 步骤2,对该相似材料逐级施加围压,每级围压稳定后测量其孔隙度,确定孔隙度 与围压的函数关系式;
[0059] 可逐级施加3MPa、8MPa、15Mpa、25MPa、30MPa及35MPa的围压。根据测得的每级围压 下的孔隙度,以围压为横坐标、孔隙度为纵坐标,在坐标轴内作出围压与孔隙度的散点图, 通过迭代法得出与散点图拟合度最高的函数关系式,即为孔隙度与围压的函数关系式。
[0060] 步骤3:对该函数关系式求导,所得导函数表示孔隙度对有效应力的变化率;
[0061] 步骤4,对该导函数进行标准化处理,得到单位应力变化引起的单位孔隙度的变化 量,该变化量即为待测实际裂隙岩体孔隙度对有效应力变化的敏感程度。
[0062] 初始围压下测得的相似材料的孔隙度为初始孔隙度,将导函数除以该初始孔隙 度,作为裂隙岩体孔隙度对有效应力敏感程度的定量表达式,该表达式表示单位应力变化 引起的单位孔隙度的变化量,表征了孔隙度对有效应力变化的敏感程度。
[0063] 以数学分析方法定量确定裂隙岩体孔隙度对有效应力的变化率,并对该变化率表 达式进行标准化处理,得到单位有效应力变化引起的单位孔隙度的变化量,可通过该变化 量准确表征孔隙度差异较大的裂隙岩体孔隙度对有效应力敏感程度。
[0064] 实施例
[0065] 以定量确定某水电站河床坝基处渗透各向异性岩体不同方向渗透率以及孔隙度 对应力敏感程度为例,本发明的确定裂隙岩体渗透率、孔隙度对有效应力敏感程度的方法 包括如下步骤:
[0066] (1)根据实际渗透各向异性裂隙岩体的体积及岩体中结构面的空间形态、数量及 分布情况,制作该裂隙岩体的三维数字模型,将三维数字模型输入3D打印机,旋转该三维数 字模型,以待测实际裂隙岩体的粉末为原材料,分别打印出与水平方向成0°、15°、30°、45°、 60 °、75 °及90 °倾角的裂隙岩体相似材料。
[0067] (2)对不同倾角相似材料分别逐级施加3MPa、8MPa、15Mpa、25MPa、30MPa及35MPa的 围压,利用岩石气体渗透率测试装置测量不同围压下不同倾角相似材料的渗透率,以围压 为横坐标,渗透率为纵坐标,做出不同倾角方向相似材料的渗透率与围压的散点图,如图1, 然后用迭代法分别计算出与各倾角方向相似材料的渗透率与围压的散点图拟合度最高的 函数关系式,作为裂隙岩体不同方向渗透率与围压的函数关系式,如表1。
[0068] 表1裂隙岩体不同方向渗透率与围压的函数关系式
[0069]
Figure CN105606513BD00081
Figure CN105606513BD00091
[0070] 同理,利用岩石气体渗透率测试装置测量不同围压下相似材料的孔隙度,以围压 为横坐标,孔隙度为纵坐标,做出相似材料的孔隙度与围压的散点图;由于不同倾角方向相 似材料孔隙度的大小以及随围压的变化规律接近一致,因此以〇°倾角相似材料在不同围压 下的孔隙度代表相似材料的孔隙度,并做出〇°倾角方向相似材料的孔隙度与围压的散点 图,如图2,并计算出与该散点图拟合度最高的函数关系式,作为裂隙岩体孔隙度与围压的 函数关系式y = 3 · 799x_Q ·026。
[0071] (3)对不同倾角方向相似材料的渗透率与围压的函数关系式进行求导,得到裂隙 岩体不同方向渗透率对围压的变化率的表达式,如表2。
[0072] 表2裂隙岩体不同方向渗透率对围压的变化率的表达式
[0073]
Figure CN105606513BD00092
Τ0074Ϊ~同理,对孔隙度与围压的函数关系式进行求导,得到裂隙岩体孔隙度对围压的变 化率的表达式,y = 0. 〇98x_1 ·026。
[0075] (4)对裂隙岩体不同方向渗透率对围压的变化率的表达式进行标准化处理,初始 围压3MPa时0°、15°、30°、45°、60°、75°及90°倾角方向处的初始渗透率分别为5.11\10_17、 5 · 07 X 10_16、7 · 38 X 10_16、5 · 23 X 10_16、2 · 05 X 10_15、1 · 33 X 10_16、5 · 87 X 10_16,以对应方向 裂隙岩体的渗透率对围压的变化率的表达式除以该方向的初始渗透率,并取绝对值,得到 裂隙渗透各向异性岩体渗透率对有效应力敏感程度的定量表达式,如表3,与倾角对应的表 达式分别表示该方向上单位应力变化引起的单位渗透率的变化量,该变化量表征了对应方 向渗透率对有效应力变化的敏感程度。
[0076] 表3裂隙岩体不同方向渗透率对有效应力敏感程度的定量表达式
[0077]
Figure CN105606513BD00093
[0078]
Figure CN105606513BD00101
[0079] 同理,对裂隙岩体不同方向孔隙度对围压的变化率的表达式进行标准化处理,初 始围压3MPa对应的初始孔隙度为5.531,以裂隙岩体的孔隙度对围压的变化率的表达式除 以该方向的初始孔隙度,并取绝对值,得到裂隙岩体孔隙度对有效应力敏感程度的定量表 达式,7 = 0.01791+°26,此表达式表示单位应力变化引起的单位孔隙度的变化量,该变化量 表征了孔隙度对有效应力变化的敏感程度。

Claims (10)

  1. CN 105606513 B 权利要求书 1/2页 1 · 一种确定裂隙岩体渗透率对有效应力敏感程度的方法,其特征在于,包括如下步骤: 步骤1,模拟实际裂隙岩体,制备与水平方向呈不同倾角的相似材料; 步骤2,选任一倾角相似材料逐级施加围压,每级围压稳定后测量其渗透率,确定该倾 角方向相似材料的渗透率与围压的函数关系式; 步骤3,对该函数关系式求导,所得导函数表示该倾角方向相似材料的渗透率对有效应 力的变化率; 步骤4,对该导函数进行标准化处理,得到该倾角方向单位应力变化引起的单位渗透率 的变化量,该变化量即为待测实际裂隙岩体该倾角方向的渗透率对有效应力变化的敏感程 度; 步骤5,重复步骤2〜4,得到待测实际裂隙岩体不同倾角方向的渗透率对有效应力变化 的敏感程度。
  2. 2. 根据权利要求1所述的确定裂隙岩体渗透率对有效应力敏感程度的方法,其特征在 于,步骤1中,所述制备与水平方向呈不同倾角的相似材料的方法包括: (1) 根据实际裂隙岩体的几何形态以及该岩体中裂隙的形态和分布,制作该裂隙岩体 的三维数字模型; (2) 将该三维数字模型输入3D打印机,旋转该三维数字模型,以实际裂隙岩体的粉末为 原料打印出与水平方向呈不同倾角的相似材料。
  3. 3. 根据权利要求1或2所述的确定裂隙岩体渗透率对有效应力敏感程度的方法,其特征 在于,所述倾角为0。、15°、30°、45°、60°、75°或90°。 4 ·根据权利要求1所述的确定裂隙岩体渗透率对有效应力敏感程度的方法,其特征在 于,步骤2中,根据测得的每级围压下的渗透率,以围压为横坐标、渗透率为纵坐标,作出该 倾角方向围压与渗透率的散点图,通过迭代法得出与散点图拟合度最高的函数关系式,即 为该倾角方向相似材料的渗透率与围压的函数关系式。
  4. 5. 根据权利要求1所述的确定裂隙岩体渗透率对有效应力敏感程度的方法,其特征在 于,步骤2中,逐级施加3MPa、8MPa、15Mpa、25MPa、30MPa及35MPa的围压。
  5. 6. 根据权利要求1所述的确定裂隙岩体渗透率对有效应力敏感程度的方法,其特征在 于,步骤4中,所述对导函数进行标准化处理的方法为:将步骤3所得的导函数除以该倾角方 向相似材料的初始渗透率,该初始渗透率为初始围压下测得的相似材料的渗透率。
  6. 7. —种确定裂隙岩体孔隙度对有效应力敏感程度的方法,其特征在于,包括如下步骤: 步骤1,模拟实际裂隙岩体,制备裂隙岩体的相似材料; 步骤2,对该相似材料逐级施加围压,每级围压稳定后测量其孔隙度,确定孔隙度与围 压的函数关系式; 步骤3:对该函数关系式求导,所得导函数表示孔隙度对有效应力的变化率; 步骤4,对该导函数进行标准化处理,得到单位应力变化引起的单位孔隙度的变化量, 该变化量即为待测实际裂隙岩体孔隙度对有效应力变化的敏感程度。
  7. 8. 根据权利要求7所述的确定裂隙岩体孔隙度对有效应力敏感程度的方法,其特征在 于,步骤2中,根据测得的每级围压下的孔隙度,以围压为横坐标、孔隙度为纵坐标,作出围 压与孔隙度的散点图,通过迭代法得出与散点图拟合度最高的函数关系式,即为孔隙度与 围压的函数关系式。 2 CN 105606513 B 权利要求书 2/2页 9 ·根据权利要求7所述的确定裂隙岩体孔隙度对有效应力敏感程度的方法,其特征在 于,步骤2中,逐级施力PMPa、8MPa、15Mpa、25MPa、3OMPa及35MPa的围压。
  8. 10.根据权利要求7所述的确定裂隙岩体孔隙度对有效应力敏感程度的方法,其特征在 于,步骤4中,所述对导函数进行标准化处理的方法为:将步骤3所得的导函数除以初始孔隙 度,该初始孔隙度为初始围压下测得的相似材料的孔隙度。 3
CN201610064446.XA 2016-01-29 2016-01-29 确定裂隙岩体渗透率及孔隙度对有效应力敏感程度的方法 Active CN105606513B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610064446.XA CN105606513B (zh) 2016-01-29 2016-01-29 确定裂隙岩体渗透率及孔隙度对有效应力敏感程度的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610064446.XA CN105606513B (zh) 2016-01-29 2016-01-29 确定裂隙岩体渗透率及孔隙度对有效应力敏感程度的方法

Publications (2)

Publication Number Publication Date
CN105606513A CN105606513A (zh) 2016-05-25
CN105606513B true CN105606513B (zh) 2018-07-24

Family

ID=55986622

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610064446.XA Active CN105606513B (zh) 2016-01-29 2016-01-29 确定裂隙岩体渗透率及孔隙度对有效应力敏感程度的方法

Country Status (1)

Country Link
CN (1) CN105606513B (zh)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107084886B (zh) * 2017-06-01 2020-01-03 河海大学 一种确定岩石有效应力系数的方法
CN108760598B (zh) * 2018-04-12 2019-12-10 河海大学 确定充填节理岩石压密所需最小循环加卸载次数的方法
CN109063334B (zh) * 2018-08-03 2020-10-02 河海大学 一种低渗透多孔介质材料的流固耦合模型构建方法
CN110243745A (zh) * 2019-06-18 2019-09-17 西安石油大学 评价增压过程中不同尺度孔隙对总孔隙压缩性影响的方法
CN110470581B (zh) * 2019-08-02 2021-05-11 中国石油大学(北京) 确定储层应力敏感程度的方法、装置及存储介质

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103207137A (zh) * 2013-04-08 2013-07-17 河海大学 全自动测定动态围压下致密岩石渗透率和孔隙度的装置
CN103528934A (zh) * 2013-10-29 2014-01-22 中国石油大学(北京) 一种测量超低渗岩石渗透率应力敏感性的互相关技术
CN103926267A (zh) * 2014-04-28 2014-07-16 西安石油大学 一种定量评价应力敏感过程中孔喉变化程度的方法
CN104515696A (zh) * 2014-12-09 2015-04-15 河海大学 应用3d打印技术的柱状节理岩体相似材料试样的制备方法
CN105158039A (zh) * 2015-09-08 2015-12-16 河海大学 一种用于裂隙岩体渗流试验仿真裂隙制作方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR101248531B1 (ko) * 2011-12-05 2013-04-03 한국지질자원연구원 이산화탄소 지중저장 매질의 공극률과 투과율 측정 장치 및 방법

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103207137A (zh) * 2013-04-08 2013-07-17 河海大学 全自动测定动态围压下致密岩石渗透率和孔隙度的装置
CN103528934A (zh) * 2013-10-29 2014-01-22 中国石油大学(北京) 一种测量超低渗岩石渗透率应力敏感性的互相关技术
CN103926267A (zh) * 2014-04-28 2014-07-16 西安石油大学 一种定量评价应力敏感过程中孔喉变化程度的方法
CN104515696A (zh) * 2014-12-09 2015-04-15 河海大学 应用3d打印技术的柱状节理岩体相似材料试样的制备方法
CN105158039A (zh) * 2015-09-08 2015-12-16 河海大学 一种用于裂隙岩体渗流试验仿真裂隙制作方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
Stress-dependence of the permeability and porosity of sandstone and shale from TCDP Hole-A;Jia-Jyun Dong,et al.;《International Journal of Rock Mechanics & Mining Sciences》;20101231;第47卷;第1141–1157页 *
低渗透岩石渗透率与孔隙率演化规律的气渗试验研究;王环玲等;《水利学报》;20150228;第46卷(第2期);第208-215页 *

Also Published As

Publication number Publication date
CN105606513A (zh) 2016-05-25

Similar Documents

Publication Publication Date Title
CN105606513B (zh) 确定裂隙岩体渗透率及孔隙度对有效应力敏感程度的方法
CN105510207B (zh) 确定裂隙岩体不同方向渗透率对孔隙度敏感程度的方法
CN107727679B (zh) 一种表征深层碳酸盐岩岩石物理学特征方法
CN108763711B (zh) 一种基于岩心扫描图像分块数值模拟的渗透率预测方法
CN106018740B (zh) 孔压静力触探标定罐系统
CN108106938B (zh) 一种实验确定酸液对致密碳酸盐岩杨氏模量影响的方法
CN108369177A (zh) 用于确定地下地层的气体渗透率的方法和装置
WO2016171971A1 (en) Nuclear magnetic resonance gas isotherm technique to evaluate reservoir rock wettability
CN108590601B (zh) 一种注水扩容膨胀施工参数优选的实验方法
CN105241802A (zh) 柱状节理岩体渗透张量测试方法及各向异性定量评估方法
CN109883889B (zh) 模拟co2在致密基质-裂缝扩散的实验装置及前缘预测方法
Hou et al. Experiments on foam texture under high pressure in porous media
CN107121354A (zh) 一种测试水下岩溶注浆材料抗分散性质的装置及其使用方法
CN110927035A (zh) 一种低渗致密砂岩束缚水饱和度计算方法
Zhang et al. A fluid‐solid coupling method for the simulation of gas transport in porous coal and rock media
CN105844011A (zh) 一种基于毛管模型的渗透率计算方法
CN105717021B (zh) 一种定量确定裂隙岩体渗透各向异性的方法
CN101377079A (zh) 一种室内测定基床系数的方法
US10732086B2 (en) Device and method for measuring magnitude of seepage force and its influence on effective stress of formation
CN109406362B (zh) 一种气水相对渗透率的确定方法
CN106248324B (zh) 长大软弱结构面强度参数的确定方法
CN110020397A (zh) 一种计算储层启动压力的方法
CN109083630B (zh) 一种用于评价钻井液封堵性能的方法
CN203929583U (zh) 一种含气页岩孔隙度和吸附参数的测试装置
Santana et al. Workflow on Incorporating Thick-Walled Cylinder Test Results in Finite Element Models of near Wellbore for Sanding Prediction Studies

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant