CN106124378A - 一种基于压汞法毛管压力曲线的复杂岩石分形维数分布谱构建方法 - Google Patents
一种基于压汞法毛管压力曲线的复杂岩石分形维数分布谱构建方法 Download PDFInfo
- Publication number
- CN106124378A CN106124378A CN201610436104.6A CN201610436104A CN106124378A CN 106124378 A CN106124378 A CN 106124378A CN 201610436104 A CN201610436104 A CN 201610436104A CN 106124378 A CN106124378 A CN 106124378A
- Authority
- CN
- China
- Prior art keywords
- fractal dimension
- sample
- formula
- matrix
- fractal
- 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.)
- Granted
Links
- QSHDDOUJBYECFT-UHFFFAOYSA-N mercury Chemical compound [Hg] QSHDDOUJBYECFT-UHFFFAOYSA-N 0.000 title claims abstract description 73
- 238000000034 method Methods 0.000 title claims abstract description 46
- 229910052753 mercury Inorganic materials 0.000 title claims abstract description 42
- 239000011435 rock Substances 0.000 title claims abstract description 34
- 238000002347 injection Methods 0.000 title claims abstract description 25
- 239000007924 injection Substances 0.000 title claims abstract description 25
- 238000010276 construction Methods 0.000 title claims abstract description 12
- 239000011148 porous material Substances 0.000 claims abstract description 15
- 238000000354 decomposition reaction Methods 0.000 claims abstract description 4
- 239000011159 matrix material Substances 0.000 claims description 37
- 238000005259 measurement Methods 0.000 claims description 10
- IJGRMHOSHXDMSA-UHFFFAOYSA-N Atomic nitrogen Chemical compound N#N IJGRMHOSHXDMSA-UHFFFAOYSA-N 0.000 claims description 6
- 238000002474 experimental method Methods 0.000 claims description 4
- 230000001186 cumulative effect Effects 0.000 claims description 3
- 238000001035 drying Methods 0.000 claims description 3
- 239000001307 helium Substances 0.000 claims description 3
- 229910052734 helium Inorganic materials 0.000 claims description 3
- SWQJXJOGLNCZEY-UHFFFAOYSA-N helium atom Chemical compound [He] SWQJXJOGLNCZEY-UHFFFAOYSA-N 0.000 claims description 3
- 239000000203 mixture Substances 0.000 claims description 3
- 229910052757 nitrogen Inorganic materials 0.000 claims description 3
- 238000005457 optimization Methods 0.000 claims description 3
- 230000017105 transposition Effects 0.000 claims description 3
- 230000000903 blocking effect Effects 0.000 claims 1
- 238000004458 analytical method Methods 0.000 abstract description 3
- 230000011218 segmentation Effects 0.000 abstract description 3
- 238000012512 characterization method Methods 0.000 abstract description 2
- BVKZGUZCCUSVTD-UHFFFAOYSA-L Carbonate Chemical compound [O-]C([O-])=O BVKZGUZCCUSVTD-UHFFFAOYSA-L 0.000 description 4
- 239000003208 petroleum Substances 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 241000208340 Araliaceae Species 0.000 description 1
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 description 1
- 235000003140 Panax quinquefolius Nutrition 0.000 description 1
- 238000000205 computational method Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 239000007789 gas Substances 0.000 description 1
- 235000008434 ginseng Nutrition 0.000 description 1
- 230000035699 permeability Effects 0.000 description 1
- 239000004576 sand Substances 0.000 description 1
- 239000000243 solution 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/088—Investigating volume, surface area, size or distribution of pores; Porosimetry
- G01N15/0886—Mercury porosimetry
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N1/00—Sampling; Preparing specimens for investigation
- G01N1/28—Preparing specimens for investigation including physical details of (bio-)chemical methods covered elsewhere, e.g. G01N33/50, C12Q
- G01N1/286—Preparing specimens for investigation including physical details of (bio-)chemical methods covered elsewhere, e.g. G01N33/50, C12Q involving mechanical work, e.g. chopping, disintegrating, compacting, homogenising
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N1/00—Sampling; Preparing specimens for investigation
- G01N1/28—Preparing specimens for investigation including physical details of (bio-)chemical methods covered elsewhere, e.g. G01N33/50, C12Q
- G01N1/44—Sample treatment involving radiation, e.g. heat
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N1/00—Sampling; Preparing specimens for investigation
- G01N1/28—Preparing specimens for investigation including physical details of (bio-)chemical methods covered elsewhere, e.g. G01N33/50, C12Q
- G01N1/286—Preparing specimens for investigation including physical details of (bio-)chemical methods covered elsewhere, e.g. G01N33/50, C12Q involving mechanical work, e.g. chopping, disintegrating, compacting, homogenising
- G01N2001/2873—Cutting or cleaving
Landscapes
- Chemical & Material Sciences (AREA)
- Biochemistry (AREA)
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Analytical Chemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Dispersion Chemistry (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
- Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)
Abstract
本发明公开了一种基于压汞法毛管压力曲线的复杂岩石分形维数分布谱构建方法,该方法基于进汞饱和度和进汞压力的幂指数关系构建了具有多个分形维数的累加公式,引用迭代奇异值分解算法求解得到多个分形维数对应的比例,从而得到分形维数分布谱。与传统的分形理论和模型相比,本发明可实现对具有分段分形特征的毛管压力曲线的细致刻画。传统的单重分形维数方法针对一块样品的毛管压力曲线只能得到一个分形维数;而通过本发明,针对一块样品的毛管压力曲线,可得到多个分形维数以及它们的比例,从而实现了分形维数分布谱的构建,为复杂岩石的孔隙结构分析和分形表征带来新的思路。
Description
技术领域
本发明涉及一种基于压汞法毛管压力曲线的复杂岩石分形维数分布谱构建方法,属于石油工程及岩石物理领域。
背景技术
分形理论自发展以来,已广泛应用于石油工程和岩石物理的各个领域,成为表征岩石孔隙结构的重要手段之一,取得了较好的应用效果。
毛管压力曲线是分形理论用于表征岩石孔隙结构的重要基础资料,以往通常根据饱和度和压力的幂指数关系构建简单的单重分形维数模型并用最小二乘法求解。随着研究的进一步深入,所面临的对象也日益复杂,从简单的高孔隙度、高渗透率砂岩向碳酸盐岩、火山岩和致密砂岩等复杂岩石发展,分形理论在应用中遇到一些瓶颈。传统的分形理论和分形维数计算方法很难将复杂岩石的毛管压力曲线进行很好的拟合,特别是对于孔隙结构极其复杂、非均质性极强的岩石而言,它们的毛管压力曲线呈明显的分段特征。为克服传统分形理论和模型在描述毛管压力曲线中的缺陷,非常有必要对分形理论进行深入分析,发展新的毛管压力曲线表征方法。
发明内容
基于上述技术问题,本发明提供一种基于压汞法毛管压力曲线的复杂岩石分形维数分布谱构建方法,为复杂岩石的孔隙结构评价和分形特征研究提供参考。
本发明所采用的技术解决方案是:
一种基于压汞法毛管压力曲线的复杂岩石分形维数分布谱构建方法,包括以下步骤:
a岩石预处理、基础数据采集和压汞实验
a1对岩石进行钻取、切割、打磨,得到呈柱塞样的样品,将样品放入烘箱,烘干后取出,用游标卡尺测量样品的长度、直径,计算样品的总体积;用电子天平测量样品的质量;
a2将样品放入孔隙度仪,采用氦气或氮气作为介质,应用波义耳定律测量得到样品的孔隙体积,计算样品孔隙度;
a3将样品放入压汞仪,采用恒速或者恒压法将汞注入样品,记录注入过程中进汞压力和汞体积的变化,进汞压力和汞体积采集的个数为m;
b将汞体积除以样品的孔隙体积,得到进汞饱和度;根据进汞饱和度和进汞压力的幂指数关系,构建多个分形维数的累加公式:
式(1)中:Dj为预设的第j个分形维数序列;aj为Dj所对应的比例;n是预设的分形维数的个数;Pc为进汞压力;S为进汞压力为Pc时的进汞饱和度;
c将式(1)进行变形,得到:
d将式(2)表示为最优化求解的矩阵形式,得到:
min||S-AX||2 (3)
式(3)中:S是由进汞饱和度形成的大小为m×1矩阵,已知量;A是由构成的大小为m×n矩阵,已知量;X是分形维数的比例构成的大小为n×1矩阵;m表示压汞法毛管压力曲线测量过程中进汞压力记录的点数;
e对式(3)中的A进行奇异值分解,得到:
A=UZVT (4)
式(4)中:U是大小为m×m的正交矩阵,V是大小为n×n的正交矩阵;VT是V的转置;Z是大小为m×n的对角矩阵;
f设置阈值q将对角矩阵进行截断,令对角矩阵Z中所有小于q的元素为0,并将Z改写为Zq,代入式(4),可得:
Aq=UZqVT (5)
g采用矩阵变换法,将式(5)代入式(3),得到X的初值:
式(6)中:是Zq的逆矩阵;
h应用迭代算法进行求解,得到最优的X
h1根据式(6)得到初值X1,设置最小平均相对误差为T和最大迭代次数为I,令i=1;
h2计算误差矩阵ΔSi=S-UZqVTXi和平均相对误差ΔTi=||ΔSi||/||S||;
h3如果ΔTi<T或者i=I,停止迭代,输出Xi;否则转入步骤h4;
h4计算和Xi=Xi+ΔXi,令i=i+1;并转入步骤h2;
k完成优化迭代步骤h后输出Xi,绘制分形维数分布谱。
优选的,预设的分形维数的个数n小于进汞压力记录的点数m。
优选的,预设的分形维数Dj最大值和最小值分别为3和2。
本发明的有益技术效果是:
本发明可实现对具有分段分形特征的毛管压力曲线的细致刻画,传统的单重分形维数方法针对一块样品的毛管压力曲线只能得到一个分形维数;而通过本发明,针对一块样品的毛管压力曲线,可得到多个分形维数以及它们的比例,从而实现了分形维数分布谱的构建,为复杂岩石的孔隙结构分析和分形表征带来新的思路。
附图说明
下面结合附图与具体实施方式对本发明作进一步说明:
图1为本发明提供的一种基于压汞法毛管压力曲线的复杂岩石分形维数分布谱构建方法的技术流程图;
图2是某复杂碳酸盐岩样品的压汞法毛管压力曲线图;
图3是应用本发明提供的方法对图2所述样品进行处理得到的分形维数分布谱。
具体实施方式
结合附图,一种基于压汞法毛管压力曲线的复杂岩石分形维数分布谱构建方法,包括以下步骤:
a岩石预处理、基础数据采集和压汞实验
a1对岩石进行钻取、切割、打磨成直径为2.54cm,长度为约3-5cm的柱塞样,将样品放入烘箱,选择温度为95℃并恒温24小时后,取出,用游标卡尺测量样品的长度、直径,计算样品的总体积;用电子天平测量样品的质量。
a2将样品放入气测法孔隙度仪,采用氦气或氮气作为介质,应用波义耳定律测量得到样品的孔隙体积,计算样品孔隙度。
a3将样品放入压汞仪,采用恒速或者恒压法将汞注入样品,记录注入过程中进汞压力和汞体积的变化,进汞压力和汞体积采集的个数为m。
b将汞体积除以样品的孔隙体积,得到进汞饱和度。根据进汞饱和度和进汞压力的幂指数关系,构建多个分形维数的累加公式:
式(1)中:Dj为预设的第j个分形维数序列;aj为Dj所对应的比例;n是预设的分形维数的个数;Pc为进汞压力;S为进汞压力为Pc时的进汞饱和度。
c将式(1)进行变形,得到:
d将式(2)表示为最优化求解的矩阵形式,得到:
min||S-AX||2 (3)
式(3)中:S是由进汞饱和度形成的大小为m×1矩阵,已知量;A是由构成的大小为m×n矩阵,已知量;X是分形维数的比例构成的大小为n×1矩阵;m表示压汞法毛管压力曲线测量过程中进汞压力记录的点数。
e对式(3)中的A进行奇异值分解,得到:
A=UZVT (4)
式中:U是大小为m×m的正交矩阵,V是大小为n×n的正交矩阵;VT是V的转置;Z是大小为m×n的对角矩阵。
f设置阈值q将对角矩阵进行截断,令对角矩阵Z中所有小于q的元素为0,并将Z改写为Zq,代入式(4),可得:
Aq=UZqVT (5)
g采用矩阵变换法,将式(5)代入式(3),得到X的初值为:
式中:是Zq的逆矩阵。
h应用迭代算法进行求解,得到最优的X
h1根据式(6)得到初值X1,设置最小平均相对误差为T和最大迭代次数为I,令i=1。
h2计算误差矩阵ΔSi=S-UZqVTXi和平均相对误差ΔTi=||ΔSi||/||S||。
h3如果ΔTi<T或者i=I,停止迭代,输出Xi;否则转入步骤h4。
h4计算和Xi=Xi+ΔXi,令i=i+1;并转入步骤h2。
k完成优化迭代步骤h后输出Xi,绘制分形维数分布谱。
上述步骤中:预设的分形维数的个数n小于进汞压力记录的点数m。
上述步骤中:预设的分形维数Dj最大值和最小值分别为3和2。
上述方法不仅仅限于压汞法毛管压力曲线的处理,还适用于半渗透率隔板法毛管压力曲线和离心法毛管压力曲线的处理。
图1是一种基于压汞法毛管压力曲线的岩石分形维数分布谱构建方法的技术流程图,主要包括岩石钻取、切割、打磨、烘干等预处理,岩石直径、长度、质量、孔隙度等基本参数测量,岩石压汞法毛管压力测量与数据整理,岩石分形维数的分布谱求解等四个部分,这四个部分缺一不可,且顺序不可颠倒。
图2是某碳酸盐岩样品实测的压汞法毛管压力曲线,压汞实验采用恒速法。需要说明的是,采用恒压法得到的毛管压力曲线,也可用本发明提供的方法进行处理。
图3是应用本发明所提供的方法对图2所述的某碳酸盐岩样品进行处理得到的分形维数的分布谱,从图中可知,该样品的分形维数主要在2到2.45之间,呈单峰分布。
上述方式中未述及的部分采取或借鉴已有技术即可实现。
需要说明的是,在本说明书的教导下,本领域技术人员所作出的任何等同替代方式,或明显变型方式,均应在本发明的保护范围之内。
Claims (3)
1.一种基于压汞法毛管压力曲线的复杂岩石分形维数分布谱构建方法,其特征在于包括以下步骤:
a岩石预处理、基础数据采集和压汞实验
a1对岩石进行钻取、切割、打磨,得到呈柱塞样的样品,将样品放入烘箱,烘干后取出,用游标卡尺测量样品的长度、直径,计算样品的总体积;用电子天平测量样品的质量;
a2将样品放入孔隙度仪,采用氦气或氮气作为介质,应用波义耳定律测量得到样品的孔隙体积,计算样品孔隙度;
a3将样品放入压汞仪,采用恒速或者恒压法将汞注入样品,记录注入过程中进汞压力和汞体积的变化,进汞压力和汞体积采集的个数为m;
b将汞体积除以样品的孔隙体积,得到进汞饱和度;根据进汞饱和度和进汞压力的幂指数关系,构建多个分形维数的累加公式:
式(1)中:Dj为预设的第j个分形维数序列;aj为Dj所对应的比例;n是预设的分形维数的个数;Pc为进汞压力;S为进汞压力为Pc时的进汞饱和度;
c将式(1)进行变形,得到:
d将式(2)表示为最优化求解的矩阵形式,得到:
min||S-AX||2 (3)
式(3)中:S是由进汞饱和度形成的大小为m×1矩阵,已知量;A是由构成的大小为m×n矩阵,已知量;X是分形维数的比例构成的大小为n×1矩阵;m表示压汞法毛管压力曲线测量过程中进汞压力记录的点数;
e对式(3)中的A进行奇异值分解,得到:
A=UZVT (4)
式(4)中:U是大小为m×m的正交矩阵,V是大小为n×n的正交矩阵;VT是V的转置;Z是大小为m×n的对角矩阵;
f设置阈值q将对角矩阵进行截断,令对角矩阵Z中所有小于q的元素为0,并将Z改写为Zq,代入式(4),可得截断后的矩阵为:
Aq=UZqVT (5)
g采用矩阵变换法,将式(5)代入式(3),得到X的初值为:
式(6)中:是Zq的逆矩阵;
h应用迭代算法进行求解,得到最优的X
h1根据式(6)得到初值X1,设置最小平均相对误差为T和最大迭代次数为I,令i=1;
h2计算误差矩阵ΔSi=S-UZqVTXi和平均相对误差ΔTi=||ΔSi||/||S||;
h3如果ΔTi<T或者i=I,停止迭代,输出Xi;否则转入步骤h4;
h4计算和Xi=Xi+ΔXi,令i=i+1;并转入步骤h2;
k完成优化迭代步骤h后输出Xi,绘制分形维数分布谱。
2.根据权利要求1所述的一种基于压汞法毛管压力曲线的复杂岩石分形维数分布谱构建方法,其特征在于:预设的分形维数的个数n小于进汞压力记录的点数m。
3.根据权利要求1所述的一种基于压汞法毛管压力曲线的复杂岩石分形维数分布谱构建方法,其特征在于:预设的分形维数Dj最大值和最小值分别为3和2。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610436104.6A CN106124378B (zh) | 2016-06-16 | 2016-06-16 | 一种基于压汞法毛管压力曲线的复杂岩石分形维数分布谱构建方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610436104.6A CN106124378B (zh) | 2016-06-16 | 2016-06-16 | 一种基于压汞法毛管压力曲线的复杂岩石分形维数分布谱构建方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106124378A true CN106124378A (zh) | 2016-11-16 |
CN106124378B CN106124378B (zh) | 2018-12-18 |
Family
ID=57470976
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610436104.6A Expired - Fee Related CN106124378B (zh) | 2016-06-16 | 2016-06-16 | 一种基于压汞法毛管压力曲线的复杂岩石分形维数分布谱构建方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106124378B (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106979918A (zh) * | 2017-04-20 | 2017-07-25 | 中国石油大学(北京) | 一种获取致密油藏岩心的液体渗透率的方法及装置 |
CN109100812A (zh) * | 2018-07-25 | 2018-12-28 | 中国石油大学(北京) | 基于核磁共振的岩石孔隙分形维数评价方法及装置 |
CN109238868A (zh) * | 2018-11-20 | 2019-01-18 | 山东农业大学 | 一种微破损现场检测混凝土结构抗压强度的测试方法 |
CN113075105A (zh) * | 2021-03-25 | 2021-07-06 | 河南理工大学 | 一种利用综合分形维Dp定量评价砂岩储层的方法 |
CN114018780A (zh) * | 2021-10-11 | 2022-02-08 | 北京市科学技术研究院分析测试研究所(北京市理化分析测试中心) | 一种陶瓷涂层孔隙和分形维数的测定方法及应用 |
CN115235977A (zh) * | 2022-07-22 | 2022-10-25 | 中国石油大学(华东) | 基于分形特征计算页岩多尺度孔径分布的方法及*** |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102735821A (zh) * | 2012-07-06 | 2012-10-17 | 中国石油大学(北京) | 测定岩石研磨性的方法及*** |
CN104573198A (zh) * | 2014-12-23 | 2015-04-29 | 长江大学 | 基于随机分形理论的数字岩心及孔隙网络模型重构方法 |
CN104712330A (zh) * | 2015-01-30 | 2015-06-17 | 中国地质大学(武汉) | 一种测井渗透率解释方法 |
CN104834807A (zh) * | 2015-03-18 | 2015-08-12 | 成都北方石油勘探开发技术有限公司 | 一种基于分形理论的应力敏感储层相对渗透率计算方法 |
CN105092820A (zh) * | 2015-08-06 | 2015-11-25 | 中国石油集团长城钻探工程有限公司地质研究院 | 一种生物碎屑灰岩储层非均质性研究方法 |
-
2016
- 2016-06-16 CN CN201610436104.6A patent/CN106124378B/zh not_active Expired - Fee Related
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102735821A (zh) * | 2012-07-06 | 2012-10-17 | 中国石油大学(北京) | 测定岩石研磨性的方法及*** |
CN104573198A (zh) * | 2014-12-23 | 2015-04-29 | 长江大学 | 基于随机分形理论的数字岩心及孔隙网络模型重构方法 |
CN104712330A (zh) * | 2015-01-30 | 2015-06-17 | 中国地质大学(武汉) | 一种测井渗透率解释方法 |
CN104834807A (zh) * | 2015-03-18 | 2015-08-12 | 成都北方石油勘探开发技术有限公司 | 一种基于分形理论的应力敏感储层相对渗透率计算方法 |
CN105092820A (zh) * | 2015-08-06 | 2015-11-25 | 中国石油集团长城钻探工程有限公司地质研究院 | 一种生物碎屑灰岩储层非均质性研究方法 |
Non-Patent Citations (2)
Title |
---|
刘义坤等: ""毛管压力曲线和分形理论在储层分类中的应用"", 《岩性油气藏》 * |
张宪国等: ""基于孔隙分形特征的低渗透储层孔隙结构评价"", 《岩性油气藏》 * |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106979918A (zh) * | 2017-04-20 | 2017-07-25 | 中国石油大学(北京) | 一种获取致密油藏岩心的液体渗透率的方法及装置 |
CN109100812A (zh) * | 2018-07-25 | 2018-12-28 | 中国石油大学(北京) | 基于核磁共振的岩石孔隙分形维数评价方法及装置 |
CN109100812B (zh) * | 2018-07-25 | 2020-06-05 | 中国石油大学(北京) | 基于核磁共振的岩石孔隙分形维数评价方法及装置 |
CN109238868A (zh) * | 2018-11-20 | 2019-01-18 | 山东农业大学 | 一种微破损现场检测混凝土结构抗压强度的测试方法 |
CN113075105A (zh) * | 2021-03-25 | 2021-07-06 | 河南理工大学 | 一种利用综合分形维Dp定量评价砂岩储层的方法 |
CN113075105B (zh) * | 2021-03-25 | 2023-06-23 | 河南理工大学 | 一种利用综合分形维Dp定量评价砂岩储层的方法 |
CN114018780A (zh) * | 2021-10-11 | 2022-02-08 | 北京市科学技术研究院分析测试研究所(北京市理化分析测试中心) | 一种陶瓷涂层孔隙和分形维数的测定方法及应用 |
CN115235977A (zh) * | 2022-07-22 | 2022-10-25 | 中国石油大学(华东) | 基于分形特征计算页岩多尺度孔径分布的方法及*** |
Also Published As
Publication number | Publication date |
---|---|
CN106124378B (zh) | 2018-12-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106124378A (zh) | 一种基于压汞法毛管压力曲线的复杂岩石分形维数分布谱构建方法 | |
CN109443867B (zh) | 一种对致密岩石的物性参数进行连续检测的方法 | |
CN104634804B (zh) | 一种利用核磁共振t2谱确定储层相对渗透率的方法 | |
EP3593159B1 (en) | Absolute porosity and pore size determination of pore types in media with varying pore sizes using nmr | |
CN109932297B (zh) | 一种致密砂岩储层渗透率的计算方法 | |
CN104278989B (zh) | 一种获取低孔低渗储层饱和度指数的方法 | |
CN107622139B (zh) | 裂缝渗透率的计算方法 | |
CN110296931B (zh) | 一种致密砂岩油水相对渗透率信息的表征方法及*** | |
CN109100812B (zh) | 基于核磁共振的岩石孔隙分形维数评价方法及装置 | |
CN105445441B (zh) | 储层含气饱和度确定方法及装置 | |
CN104374683A (zh) | 一种岩心孔隙压缩系数测试装置及其测试方法 | |
CN107688037B (zh) | 一种利用核磁测井t2分布确定井下岩石粒度曲线的方法 | |
CN102565858A (zh) | 一种多孔介质含水饱和度的计算方法 | |
CN114239367A (zh) | 一种室内岩心的数字化多相流固耦合渗流数值模拟方法 | |
CN110161071A (zh) | 一种评价致密储层可动流体最小喉道半径的方法 | |
CN114283254B (zh) | 基于核磁共振数据的岩心数字化孔隙网络模型构建方法 | |
CN102323199B (zh) | 一种孔隙度测定实验设计方法及其设计*** | |
CN112505084A (zh) | 一种注气提高页岩油可动性的评价模型、评价方法、应用 | |
Carles et al. | Low permeability measurements using steady-state and transient methods | |
CN108007954A (zh) | 一种致密岩心束缚水环境建立的实验装置及测试方法 | |
CA3107041A1 (en) | Assessment of inaccessible pore volume for polymer flooding | |
CN106323833A (zh) | 岩心孔隙度测量方法及装置 | |
CN103675722A (zh) | 岩石t2-g实验采集参数自动匹配方法 | |
CN104374679B (zh) | 一种确定储层渗透率贡献值的方法 | |
CN105138761B (zh) | 降雨条件下坡面地表糙率及土壤吸湿率的估算方法 |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20181218 Termination date: 20200616 |