CN106841002B - 确定岩石试验中产生滑脱效应的隐节理密度阈值的方法 - Google Patents
确定岩石试验中产生滑脱效应的隐节理密度阈值的方法 Download PDFInfo
- Publication number
- CN106841002B CN106841002B CN201710050242.5A CN201710050242A CN106841002B CN 106841002 B CN106841002 B CN 106841002B CN 201710050242 A CN201710050242 A CN 201710050242A CN 106841002 B CN106841002 B CN 106841002B
- Authority
- CN
- China
- Prior art keywords
- sample
- blind joint
- joint density
- pressure
- slippage effect
- 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
Links
- 239000011435 rock Substances 0.000 title claims abstract description 51
- 230000000694 effects Effects 0.000 title claims abstract description 38
- 238000012360 testing method Methods 0.000 title claims abstract description 25
- 238000000034 method Methods 0.000 title claims abstract description 21
- 230000035699 permeability Effects 0.000 claims abstract description 47
- 230000003204 osmotic effect Effects 0.000 claims abstract description 18
- 239000007789 gas Substances 0.000 claims description 41
- 239000002994 raw material Substances 0.000 claims description 12
- ATJFFYVFTNAWJD-UHFFFAOYSA-N Tin Chemical compound [Sn] ATJFFYVFTNAWJD-UHFFFAOYSA-N 0.000 claims description 6
- 239000011083 cement mortar Substances 0.000 claims description 6
- 238000002360 preparation method Methods 0.000 claims description 6
- 238000001764 infiltration Methods 0.000 claims description 5
- 230000008595 infiltration Effects 0.000 claims description 5
- 230000035515 penetration Effects 0.000 claims description 4
- 239000011261 inert gas Substances 0.000 claims description 3
- 238000012423 maintenance Methods 0.000 claims description 3
- 238000004519 manufacturing process Methods 0.000 claims description 3
- 239000000463 material Substances 0.000 claims description 3
- 238000012887 quadratic function Methods 0.000 claims description 3
- 238000002347 injection Methods 0.000 claims description 2
- 239000007924 injection Substances 0.000 claims description 2
- 238000000691 measurement method Methods 0.000 claims description 2
- 238000000465 moulding Methods 0.000 claims description 2
- 238000005259 measurement Methods 0.000 abstract description 3
- 239000000523 sample Substances 0.000 description 51
- 239000007788 liquid Substances 0.000 description 6
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 4
- 239000004568 cement Substances 0.000 description 2
- 238000002474 experimental method Methods 0.000 description 2
- 239000011148 porous material Substances 0.000 description 2
- 239000007787 solid Substances 0.000 description 2
- 239000000243 solution Substances 0.000 description 2
- 238000003860 storage Methods 0.000 description 2
- 239000002518 antifoaming agent Substances 0.000 description 1
- 238000009412 basement excavation Methods 0.000 description 1
- 239000003638 chemical reducing agent Substances 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000005065 mining Methods 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
- 230000003014 reinforcing effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 239000004576 sand Substances 0.000 description 1
- 238000009738 saturating Methods 0.000 description 1
- 238000010998 test method Methods 0.000 description 1
- 239000002699 waste material Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
- G01N15/08—Investigating permeability, pore-volume, or surface area of porous materials
- G01N15/082—Investigating permeability by forcing a fluid through a sample
- G01N15/0826—Investigating permeability by forcing a fluid through a sample and measuring fluid flow rate, i.e. permeation rate or pressure change
Landscapes
- Chemical & Material Sciences (AREA)
- Physics & Mathematics (AREA)
- Analytical Chemistry (AREA)
- Dispersion Chemistry (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Fluid Mechanics (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
本发明公开一种确定岩石试验中产生滑脱效应的隐节理密度阈值的方法,包括如下步骤:(1)制备不同隐节理密度的试样;(2)从隐节理密度最低的试样开始,测量该试样在相同围压不同渗压下的气体渗透率,并记录压力室的入口气体压力和出口气体压力、取两者的平均值作为平均压力;(3)绘制该试样的气体渗透率随平均压力倒数的变化曲线,通过函数对该曲线进行拟合;(4)若拟合相关系数大于0.90,则选择隐节理密度大于此试样的下一个试样,重复上述步骤;若相关系数小于0.90,则此试样的隐节理密度即为产生滑脱效应的隐节理密度阈值。本方法可将滑脱效应与隐节理密度直接联系起来,更直观、快捷的对不同隐节理密度试样气体渗透率测量结果进行评估。
Description
技术领域
本发明属于地下洞室、采矿、地下核废料储存、水利水电等岩石工程中以节理岩体为研究对象的工程领域,具体是涉及一种确定含隐节理的致密岩石试验中产生滑脱效应的隐节理密度阈值的方法。
背景技术
当液体在多孔介质中流动时,孔道中心液体的流速比孔壁处液体的流速明显高,孔壁处固液接触面的液体速度流动速度几乎为零。而当气体在多孔介质中流动时,除了气体分子之间发生碰撞外,还有气体分子和孔壁之间的碰撞(固气边界的联结力),当岩石孔隙半径接近气体分子自由程,气体分子和固体壁之间的碰撞频率增加,这样在边壁表面就产生了一个多余的流量,即滑脱流。导致气体测量的渗透率大于其绝对渗透率也大于液测渗透率,这种现象称为滑脱效应,又称林肯伯格效应(Klinkenberg effect)。由于滑脱效应的存在,使得气测渗透率大于试样本身固有渗透率而不能反映试样的真实渗透属性。特别是在孔隙度和渗透率极低的岩石内部,气体平均自由程和岩石内部平均孔隙半径接近时,管壁处的气体分子流动会更快,会产生明显的滑脱效应。
目前,地下水封洞库利用低渗透岩石的致密性来进行石油、石油液化气的储存。此类岩石一般有整体性好、裂隙少、弱透水性等特点,此时的气测渗透率的滑脱效应亦明显,而随着洞室的开挖与加固岩石的围压亦发生变化,导致岩石内部部分微裂隙扩展,当裂隙数量增大到一定程度时,渗透率也在显著增加,气渗中的滑脱效应也会随之减小到忽略不计。因此,研究岩石微裂纹的扩展中,产生滑脱效应的隐节理密度阈值对确定气测渗透率是否考虑滑脱效应具有非常重要的工程意义和科学意义。
发明内容
发明目的:本发明的针对现有技术中由于滑脱效应存在导致的气测渗透率不能反映试样的真实渗透的问题,提供一种确定岩石试验中产生滑脱效应的隐节理密度阈值的方法。
技术方案:本发明的确定岩石试验中产生滑脱效应的隐节理密度阈值的方法,包括如下步骤:
步骤1,制备不同隐节理密度的致密岩石试样;隐节理密度具体为节理法线方向上单位长度内的节理条数;
具体的,致密岩石试样的制备方法可包括下述步骤:
步骤11,依据天然致密岩石的强度及渗透特性,配制水泥砂浆作为制作试样的原料;
步骤12,将配制好的原料注入模具中,注至模具一半的位置,垂直放入用于模拟实际节理的锡片,然后继续注入原料,直至试样完全成形;
步骤13,养护、至原料固化成型,然后脱模、取出试样养护,试样制备完成。
优选的,步骤11中,控制水泥砂浆的配比,使得配制的原料制成不含隐节理的完整样时,完整样的渗透率量级与实际天然岩石的渗透率量级相同。实际天然岩石的渗透率量级为10-14m2~10-22m2。
步骤2,从隐节理密度最低的试样开始,测量该试样在相同围压不同渗压下的气体渗透率,并记录每一渗压下压力室的入口气体压力和出口气体压力、取两者的平均值作为平均压力;
具体的,试样在相同围压不同渗压下的气体渗透率测量方法包括下述步骤:
步骤21,将试样放入压力室中,加围压至与工程实际符合;
步骤22,待围压稳定后,设定不同的渗压等级,从低到高依次测量试样的气体渗透率。
优选的,利用致密岩石惰性气体测试系统测试试样的气体渗透率。
进一步的,不同渗压取自0.2~2.5MPa的压力范围。
步骤3,绘制该试样的气体渗透率随平均压力倒数的变化曲线,并通过函数对该曲线进行拟合;较优的,通过二次函数对试样的气体渗透率随平均压力倒数的变化曲线进行拟合。
步骤4,若拟合得到的相关系数大于0.90,则选择隐节理密度大于此试样的下一个试样,重复步骤1~4;若相关系数小于0.90,则此试样的隐节理密度即为产生滑脱效应的隐节理密度阈值。
有益效果:与现有技术相比,本发明的显著优点在于:本发明将滑脱效应与隐节理密度直接联系起来,通过测量不同隐节理密度试样的气体渗透率,并对其结果进行评估,更直观、快捷的确定出隐节理密度的阈值,从而可有效判断实际工程中是否需要考虑气体渗透率的滑脱效应,为工程中预测含隐节理致密岩石的渗透率提供方法和科学依据,对实际工程有极大的运用价值。
附图说明
图1为本发明确定岩石试验中产生滑脱效应的隐节理密度阈值的方法流程图;
图2为隐节理密度为20条时,岩样的渗透率与平均压力倒数的函数关系拟合曲线;
图3为隐节理密度为40条时,岩样的渗透率与平均压力倒数的函数关系拟合曲线;
图4为隐节理密度为60条时,岩样的渗透率与平均压力倒数的函数关系拟合曲线;
图5为隐节理密度为80条时,岩样的渗透率与平均压力倒数的函数关系拟合曲线。
具体实施方式
下面结合附图及具体实施例对本发明的技术方案作进一步说明。
本发明的确定岩石试验中产生滑脱效应的隐节理密度阈值的方法通过运用物模实验的手段,制备出不同隐节理密度的试样模型,测出不同渗压下试样的气体渗透率,再对其关系进行函数拟合进而判断是否存在滑脱效应,反复多次实验过程进而确定产生滑脱效应的隐节理密度阈值,为工程中预测含隐节理致密岩石的渗透率提供方法和科学依据,对实际工程有极大的运用价值。
如图1,本发明的确定在岩样试验中产生滑脱效应的隐节理密度阈值的方法,包括以下步骤:
步骤1,制备不同隐节理密度的致密岩石试样;隐节理密度具体为节理法线方向上单位长度内的节理条数,分别制备隐节理密度为20条、40条、60条和80条的试样;
具体的,致密岩石试样的制备方法包括如下步骤:
(1)选用0.1mm的锡片模拟实际隐节理,将锡片剪成边长10mm的正方形;
(2)依据天然致密岩石砂岩的强度及渗透特性,配制水泥砂浆作为制作岩样的原料,要求控制水泥砂浆的配比,使得配制的原料制成不含隐节理的完整样时,完整样的渗透率量级与实际天然岩石的渗透率量级(10-14m2~10-22m2)保持一致;
本次试验经多次测量标定,最终确定用于制备试样原料的组分配比,其中,水泥、砂、水等的配比为1:0.5:0.35,另外加入水泥质量0.2%的消泡剂和减水剂,依此配比作为制作试样的原料。
(3)将原料由注射器抽取后,注射到模具中。先注射到模具内一半的位置,根据隐节理密度垂直放入对应数量剪好的锡片,使裂隙面与渗流方向一致;当锡片放好后继续注射原料至试样完全成形,制成试样为直径50mm、高50mm的圆柱样。
(4)养护30h,然后脱模,取出试样养护,试样制备完成。
步骤2,从隐节理密度最低的试样开始,将试样在致密岩石惰性气体测试系统中进行不同渗压下的气体渗透率测量,具体测量步骤包括:
(1)首先将隐节理密度为20条的试样放入压力室,加围压至与工程实际较符合大小,本次试验围压加至25MPa;
(2)待围压稳定后,设定不同的渗压等级,由低到高,依次测量该试样在不同渗压下的气体渗透率;渗压的具体取值可在0.2~2.5MPa的范围内,本次实验选取在渗压0.2MPa、0.4MPa、0.6MPa、0.8MPa和1.0MPa下测试气体渗透率;
本步骤中,测量气体渗透率的同时,记录每一渗压下压力室的入口气体压力P1和出口气体压力P0,设置平均压力Pm,Pm=(P1+P0)/2;
步骤3,根据获得的试验数据绘制该试样的气体渗透率随平均压力倒数的变化曲线,并用二次函数对渗透率与平均压力倒数的函数关系进行拟合,如图2;
步骤4,若拟合得到的相关系数大于0.90,则选择隐节理密度大于此试样的下一个试样,重复步骤1~4;若拟合得到的相关系数小于0.90,则此试样的隐节理密度即为产生滑脱效应的隐节理密度阈值。
通过对隐节理密度分别为20条、40条、60条、80条试样的试验,拟合曲线如图2~5,当隐节理密度为80条时,得到的相关系数为0.86,小于0.90,因此得到该类试样在此围压等级下产生滑脱效应的隐节理密度阈值为80条。说明小于此隐节理密度的岩石气测渗透率需要考虑其滑脱效应的影响,气测渗透率需要修正,否则过高的估计了其渗透性能的大小。
Claims (8)
1.一种确定岩石试验中产生滑脱效应的隐节理密度阈值的方法,其特征在于,包括如下步骤:
步骤1,制备不同隐节理密度的致密岩石试样;
步骤2,从隐节理密度最低的试样开始,测量该试样在相同围压不同渗压下的气体渗透率,并记录每一渗压下压力室的入口气体压力和出口气体压力、取两者的平均值作为平均压力;
步骤3,绘制该试样的气体渗透率随平均压力倒数的变化曲线,并通过函数对该曲线进行拟合;
步骤4,若拟合得到的相关系数大于0.90,则选择隐节理密度大于此试样的下一个试样,重复步骤1~4;若相关系数小于0.90,则此试样的隐节理密度即为产生滑脱效应的隐节理密度阈值。
2.根据权利要求1所述的确定岩石试验中产生滑脱效应的隐节理密度阈值的方法,其特征在于,步骤1中,所述致密岩石试样的制备方法包括下述步骤:
步骤11,依据天然致密岩石的强度及渗透特性,配制水泥砂浆作为制作试样的原料;
步骤12,将配制好的原料注入模具中,注至模具一半的位置,垂直放入用于模拟实际节理的锡片,然后继续注入原料,直至试样完全成形;
步骤13,养护、至原料固化成型,然后脱模、取出试样养护,试样制备完成。
3.根据权利要求2所述的确定岩石试验中产生滑脱效应的隐节理密度阈值的方法,其特征在于,步骤11中,控制水泥砂浆的配比,使得配制的原料制成不含隐节理的完整样时,完整样的渗透率量级与实际天然岩石的渗透率量级相同。
4.根据权利要求2所述的确定岩石试验中产生滑脱效应的隐节理密度阈值的方法,其特征在于,所述实际天然岩石的渗透率量级为10-14m2~10-22m2。
5.根据权利要求1所述的确定岩石试验中产生滑脱效应的隐节理密度阈值的方法,其特征在于,步骤2中,利用致密岩石惰性气体测试系统测试试样的气体渗透率。
6.根据权利要求1所述的确定岩石试验中产生滑脱效应的隐节理密度阈值的方法,其特征在于,步骤2中,所述试样在相同围压不同渗压下的气体渗透率测量方法包括:
步骤21,将试样放入压力室中,加围压至与工程实际符合;
步骤22,待围压稳定后,设定不同的渗压等级,从低到高依次测量试样的气体渗透率。
7.根据权利要求1所述的确定岩石试验中产生滑脱效应的隐节理密度阈值的方法,其特征在于,步骤2中,所述不同渗压取自0.2~2.5MPa的压力范围。
8.根据权利要求1所述的确定岩石试验中产生滑脱效应的隐节理密度阈值的方法,其特征在于,步骤3中,通过二次函数对试样的气体渗透率随平均压力倒数的变化曲线进行拟合。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710050242.5A CN106841002B (zh) | 2017-01-23 | 2017-01-23 | 确定岩石试验中产生滑脱效应的隐节理密度阈值的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710050242.5A CN106841002B (zh) | 2017-01-23 | 2017-01-23 | 确定岩石试验中产生滑脱效应的隐节理密度阈值的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106841002A CN106841002A (zh) | 2017-06-13 |
CN106841002B true CN106841002B (zh) | 2019-05-28 |
Family
ID=59120670
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710050242.5A Active CN106841002B (zh) | 2017-01-23 | 2017-01-23 | 确定岩石试验中产生滑脱效应的隐节理密度阈值的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106841002B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108760598B (zh) * | 2018-04-12 | 2019-12-10 | 河海大学 | 确定充填节理岩石压密所需最小循环加卸载次数的方法 |
CN108956412B (zh) * | 2018-04-12 | 2020-10-02 | 河海大学 | 一种确定高海拔寒区岩石冻融损伤程度的方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103063557A (zh) * | 2012-12-31 | 2013-04-24 | 河海大学 | 一种岩石气体渗透测定装置及测定方法 |
CN103163057A (zh) * | 2013-03-18 | 2013-06-19 | 河海大学 | 一种致密岩石材料气体渗透率测试装置及测算方法 |
CN104713812A (zh) * | 2015-04-01 | 2015-06-17 | 西南石油大学 | 一种基于岩心气测渗透率测量装置的校正方法 |
CN105651570A (zh) * | 2015-12-29 | 2016-06-08 | 河海大学 | 一种含无充填隐节理的柱状节理岩体试样的制备方法 |
CN105651572A (zh) * | 2015-12-29 | 2016-06-08 | 河海大学 | 一种含充填隐节理的柱状节理岩体试样的制备方法 |
-
2017
- 2017-01-23 CN CN201710050242.5A patent/CN106841002B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103063557A (zh) * | 2012-12-31 | 2013-04-24 | 河海大学 | 一种岩石气体渗透测定装置及测定方法 |
CN103163057A (zh) * | 2013-03-18 | 2013-06-19 | 河海大学 | 一种致密岩石材料气体渗透率测试装置及测算方法 |
CN104713812A (zh) * | 2015-04-01 | 2015-06-17 | 西南石油大学 | 一种基于岩心气测渗透率测量装置的校正方法 |
CN105651570A (zh) * | 2015-12-29 | 2016-06-08 | 河海大学 | 一种含无充填隐节理的柱状节理岩体试样的制备方法 |
CN105651572A (zh) * | 2015-12-29 | 2016-06-08 | 河海大学 | 一种含充填隐节理的柱状节理岩体试样的制备方法 |
Non-Patent Citations (1)
Title |
---|
致密岩石气体渗流滑脱效应试验研究;王环玲等;《岩土工程学报》;20160531;第38卷(第5期);第777-784页 |
Also Published As
Publication number | Publication date |
---|---|
CN106841002A (zh) | 2017-06-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106884635B (zh) | 一种低、特低渗透油藏co2驱最小混相压力的测定方法 | |
CN105067781B (zh) | 泡沫驱油评价装置及其评价方法 | |
CN103267722B (zh) | 一种承压渗透注浆加固试验装置及方法 | |
CN107063963B (zh) | 一种致密储层微裂缝扩展及渗流特征的测试装置和方法 | |
CN103471976B (zh) | 一种测量含水合物多孔沉积物渗透率的装置 | |
CN103513019B (zh) | 固井水泥浆气窜模拟评价仪 | |
CN101354334B (zh) | 基于瞬态压力脉冲法的原位小型渗透系数测量系统 | |
CN106121601A (zh) | 一种泡沫驱油物理模拟装置及方法 | |
CN110242290B (zh) | 一种页岩岩芯层理缝流动模拟评价装置及方法 | |
CN103969160B (zh) | 高温高压泡沫液体动滤失检测系统及其检测方法 | |
CN114354809B (zh) | 一种二氧化碳脉冲驱替置换甲烷的实验系统及实验评价方法 | |
CN202814960U (zh) | 裂缝性底水气藏水侵动态物理模拟实验装置 | |
CN205538580U (zh) | 裂隙介质系统渗透张量室内测定装置 | |
CN106841002B (zh) | 确定岩石试验中产生滑脱效应的隐节理密度阈值的方法 | |
CN111878075B (zh) | 一种倾斜油藏气水协同注入分区采出程度的测试方法 | |
CN111006952A (zh) | 高压渗透注浆加固裂隙岩石试样的实验测试装置及注浆方法 | |
CN105510207A (zh) | 确定裂隙岩体不同方向渗透率对孔隙度敏感程度的方法 | |
CN114152510A (zh) | 用于富水破碎岩层动水注浆加固模型试验装置及试验方法 | |
CN105298488A (zh) | 非连续充填方式下导流能力测试方法 | |
CN108828190B (zh) | 一种裂缝性致密砂岩油气藏的裂缝模拟方法 | |
CN104695952B (zh) | 非常规气藏地层参数的解释方法及解释系统 | |
CN107907464B (zh) | 一种压裂用渗透石水泥浆性能测定装置及方法 | |
Guo et al. | Water invasion and remaining gas distribution in carbonate gas reservoirs using core displacement and NMR | |
CN105675441A (zh) | 一种径向流条件下不同水合物饱和度时的气水相对渗透率测定方法 | |
CN217542768U (zh) | 一种可视化泡沫驱替实验装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |