非常规天然气

基于四参数随机生长模型的页岩储层应力敏感分析

  • 徐加祥 , 1, 2, 3 ,
  • 杨立峰 2, 3 ,
  • 丁云宏 1 ,
  • 刘哲 2, 3 ,
  • 高睿 2, 3 ,
  • 王臻 2, 3
展开
  • 1. 中国石油勘探开发研究院, 北京 100083
  • 2. 中国石油勘探开发研究院压裂酸化技术服务中心, 河北 廊坊 065007
  • 3. 中国石油天然气集团公司油气藏改造重点实验室, 河北 廊坊 065007

徐加祥(1991-),男,山东泰安人,博士研究生,主要从事油气藏增产改造技术研究. E-mail:

收稿日期: 2019-01-08

  修回日期: 2019-04-16

  网络出版日期: 2019-10-14

基金资助

国家科技重大专项“致密油储层高效体积改造技术”(2016ZX05046-004)

Stress sensitivity analysis of the shale reservoir by the quartet structure generation set

  • Jia-xiang Xu , 1, 2, 3 ,
  • Li-feng Yang 2, 3 ,
  • Yun-hong Ding 1 ,
  • Zhe Liu 2, 3 ,
  • Rui Gao 2, 3 ,
  • Zhen Wang 2, 3
Expand
  • 1. PetroChina Research Institute of Petroleum Exploration & Development, Beijing 100083, China
  • 2. Fracturing and Acidification Service Center of Research Institution of Petroleum Exploration & Development, PetroChina, Langfang 065007, China
  • 3. CNPC Key Laboratory for Reformation of Oil and Gas Reservoirs, Langfang 065007, China

Received date: 2019-01-08

  Revised date: 2019-04-16

  Online published: 2019-10-14

本文亮点

页岩储层的应力敏感性是影响其后期开发效果的关键因素,从微观的角度深入认识其应力敏感机理及其影响因素对页岩气的开发具有重要意义。借助四参数随机生长模型构建了不同孔隙度和不同孔隙大小分布的岩心样本,利用弹性力学理论模拟了不同有效应力作用下各个岩心孔隙半径的分布变化及其对岩心固有渗透率的影响,深入分析了孔隙大小及其形状因子与上述两者之间的关系。结果表明,导致页岩应力敏感的直接原因是有效应力作用下孔隙面积的减小及孔隙位置的迁移。有效应力的增大使得各孔隙半径均有减小,孔隙半径的减小比例分别与孔隙初始面积和孔隙的形状因子呈正相关关系和负相关关系。在相同的孔隙度条件下,孔隙半径越均匀,平均孔隙半径越小,应力敏感性越强。有效应力的增加使得岩心固有渗透率呈指数型下降且孔隙度越小、固有渗透率越低的岩心,其应力敏感性越强,孔隙度对固有渗透率的影响大于孔隙半径均匀性的影响。

本文引用格式

徐加祥 , 杨立峰 , 丁云宏 , 刘哲 , 高睿 , 王臻 . 基于四参数随机生长模型的页岩储层应力敏感分析[J]. 天然气地球科学, 2019 , 30(9) : 1341 -1348 . DOI: 10.11764/j.issn.1672-1926.2019.04.014

Highlights

The stress sensitivity of permeability in the shale reservoir is the key factor influencing its later development effect. It is of great significance for shale gas development to deeply understand its stress sensitivity mechanism and its influencing factors from the microscopic point of view. In this paper, core samples with different porosity and pore size distribution are constructed by means of the quartet structure generation set. The variation of pore radius distribution and its influence on intrinsic permeability of cores under different effective stress are simulated by using elasticity theory. The relationship between pore size and shape factor and the above two parameters is deeply analyzed. The results show that the direct cause of shale stress sensitivity is the decrease of pore area and the migration of pore position under effective stress. With the increase of effective stress, the pore radius decreases. The proportion of reduction of pore radius is positively and negatively correlated with the initial pore area and the pore shape factor, respectively. Under the same porosity condition, the more uniform the pore radius, the stronger the stress sensitivity. With the increase of effective stress, the intrinsic permeability of core decreases exponentially, and the smaller the porosity and the lower the intrinsic permeability, the stronger the stress sensitivity of core, and the effect of porosity on intrinsic permeability is obviously greater than that of pore radius uniformity.

0 引言

目前,页岩气和致密气等非常规资源已成为我国天然气产量持续增长的重要支柱[1,2]。在实际开发过程中,游离气和吸附气的采出使得储层压力不断降低,基质岩石所受到的有效应力持续增加,岩石孔隙度和渗透率等物性参数随之变化,进而影响气井产能[3,4]。针对致密储层的应力敏感问题,国内外学者从不同角度作了大量研究。Chalmers等[5]和Bustin等[6]通过改变页岩上覆压力使其有效应力变化并采用脉冲法测量,结果表明页岩具有强应力敏感性。Kwon等[7]和郭为等[8]利用脉冲法研究了页岩内外压力变化对其应力敏感的影响。张睿等[9]也利用脉冲法测量了页岩岩心的应力敏感曲线并借助毛管理论推导了孔隙半径和渗透率与有效应力的关系。张睿等[10]以Walsh模型与幂律公式的内在联系为基础,归纳了5种不同孔隙类型所适用的应力敏感公式。张睿等[11]从渗透率应力敏感机理出发分别得到了针对毛细管、平板裂缝和双重孔隙介质渗透率应力敏感系数的通用表达式。肖文联等[12]根据实验数据和岩石微观特征,利用3种应力敏感系数,划分了致密砂岩应力敏感性的标准。雷刚等[13]根据Hertz接触变形法则,定量表征了非均质致密砂岩毛管孔隙渗透率应力敏感性。李传亮[14]由应力敏感指数与压缩系数之间的关系建立了两者之间的理论公式。
以上研究中实验方法的岩心获取存在诸多困难,无法定量表示孔隙结构,而理论研究往往局限于对应力敏感系数的刻画,对影响应力敏感的内在因素缺乏深入分析。本文借助数字岩心技术对不同孔隙结构的岩心进行重构,采用弹性力学原理对不同有效应力作用下孔隙半径及岩心固有渗透率的变化进行分析,深入揭示了孔隙大小及其形状因子对岩心应力敏感的影响。

1 岩心重构

数字岩心是研究孔隙介质微观渗流特征的重要技术手段,建立数字岩心主要有物理方法[15,16]和数学方法[17,18],其中物理方法虽然可以还原真实岩心特征但是费用昂贵且耗时较长,难以有针对性的获得具有不同储层参数的岩心。铸体薄片以及粒度分析资料等岩心二维信息相对比较容易获得,可以在此基础上利用数学方法重构不同性质的数字岩心。
四参数随机生长法(QSGS)是一种构建土体微观结构的方法,其中的四参数指初始分布概率P c、生长概率Pi 、概率密度p i nm和孔隙度φ。该方法物理意义明确,能够考虑岩石颗粒的随机分布及统计特性,通过调整参数可以实现页岩多种微观形貌特征的重构[19,20]
假设岩石完全饱和,其固相颗粒为生长相,孔隙为非生长相,初始条件下构造区全部为孔隙,固相颗粒的构造过程如下:
(1)在构造区域的网格点上随机布置固相颗粒的生长核,初始分布概率为P c,并确保P c不大于最终岩石中固相体积分数(P c≤1-φ)。
(2)在每个生长核的网格点上,按照给定不同方向上的生长概率Pi i代表生长方向,i=1,2…,8),向周围相邻网格点上生长,生长方向如图1所示。其中,1~4是4个主要生长方向,5~8为4个次要生长方向。为了获得各向均质的微观结构,设定主要生长方向和次要生长方向的生长比例为:P 1-4P 5-8=4∶1。
图1 四参数随机生长法固相生长方向示意

Fig.1 Diagram of solid growth directions in QSGS

(3)重复步骤(2),直至非生长相(孔隙)的体积分数达到给定的孔隙率φ。
本文所重构的岩心为各向同性的,如上文所述固相颗粒在主要生长方向和次要生长方向的生长比例为:P 1-4P 5-8=4∶1,即P 1-4=0.04,P 5-8=0.01。设置初始分布概率为0.002,分别构建孔隙度为1%、5%和10%的岩心,如图2(a)、图2(b)和图2(c)所示。然后固定孔隙度为10%,分别构造初始概率为0.000 2和0.02的岩心,如图2(d)和图2(e)所示,以获得孔隙度相同但孔隙半径分布不同的岩心。
图2 不同结构孔隙构造

Fig.2 Different pore structures in digital cores

识别并提取上述岩心中所有孔隙半径并统计其分布频率,如图3所示。可以看出5种岩心孔隙半径均符合对数正态分布,该分布特征与Morrow[21]得到的孔径分布结论一致,说明了所构建岩心的准确性。同时,可以发现当岩心孔隙度较小时,其内部孔隙的半径主要以小孔隙(0.2~0.5μm)为主,随着孔隙度的增大,岩心内部小孔隙比例有所降低,中等孔隙(0.7~1.2μm)比例有所增加。在相同孔隙度条件下,随着构造初始概率的减小,岩心中孔隙大小愈发不均匀,孔隙半径范围增大,最大孔隙半径接近3.5μm,该认识在图2中有比较直观的体现。
图3 不同结构岩心中孔隙半径分布

Fig.3 Distribution of pore diameter in different digital cores

2 岩心应力敏感模型

对于非常规气藏,在未开发的条件下,地层应力σ、有效应力σ´和孔隙压力P处于平衡状态,如图4所示。在开发过程中,由于采用衰竭式开发,游离气和吸附气被采出后,孔隙压力下降,导致基岩受到的有效应力增加,对于页岩这样的脆性材料,研究证明其在载荷作用下的应力—应变呈现弹性关系[22]。在有效应力作用下,页岩的应力—位移关系表示如下:
图4 岩心中有效应力示意

Fig.4 Effective stress in pore structure

ρ 2 u t 2 = · σ ' + F v
式中:ρ为岩石颗粒密度,kg/m3u为岩石颗粒位移,mm;σ´为有效应力张量; F v为单位体积所受外力,N。其应力—应变关系及应变—位移关系如式(2)和式(3)所示:
σ = C ε
ε = 1 2 ( u ) T + u
对于二维的应力—应变关系,其刚度矩阵为4阶张量,如下:
C = E 1 + υ 1 + 2 υ 1 - υ υ 0 0 υ 1 - υ 0 0 0 0 1 - 2 υ 2 0 0 0 0 1 - 2 υ 2
式中:ε为应变张量;E为杨氏模量,GPa;υ为泊松比。
岩心中孔隙形状一般是不规则的,为了研究其不规则程度对应力敏感性的影响,这里采用形状因子对孔隙形状的不规则性进行刻画:
G = A C 2
式中:A为单个孔隙的面积,μm2C为该孔隙的周长,μm。
页岩的固有渗透率与其孔隙度、孔隙半径以及孔隙的迂曲度等因素有关,一般采用Kozeny⁃Carman方程进行描述[23],表示如下:
k = C k c ϕ 3 1 - ϕ 2
其中, C k c = c 8 a v 2 τ
式中:φ为孔隙度,%;c为经验几何系数;a v为比表面积,1/μm;τ为孔隙迂曲度。

3 模拟结果

基于上述模型,设置岩石基质弹性模量为23GPa,泊松比为0.15,岩石密度为2 600kg/m3,在介质上下和左右方向分别设置最大和最小主应力,在岩石孔隙中设置油藏压力,限制研究介质的自由运动,分别对不同有效应力作用下页岩岩心中孔隙面积、孔隙周长以及孔隙渗透率进行模拟,并分析了其影响因素。

3.1 孔隙半径变化

图2中5种不同孔隙结构岩心在不同有效应力作用下孔隙半径的变化如图6所示。图2(c)中部分孔隙在有效应力为20MPa时,孔隙结构的变化如图5所示。
图6 不同孔隙结构岩心中孔隙半径随有效应力变化

Fig.6 Diameter variation of the pore in different cores under the effective stress

图5 有效应力作用下孔隙结构变化

Fig.5 Variation of the pore structure under the effective stress

图5中黑色线框为原始孔隙形状,蓝色区域为变形后孔隙部分,红色区域为基质部分。可以看出,岩心中孔隙在有效应力的作用下有大幅度减小,在这个过程中会使得孔隙原本狭小的区域闭合,连通的孔隙被截断成一个主要孔隙及分散在其周围的多个小孔隙。同时,基质变形还会使孔隙位置发生偏移,这一现象可能会在一定程度上增大岩石孔隙的迂曲度,进而增大其中流体的渗流阻力。
图6中孔隙半径的变化可以看出,随着有效应力的增加,孔隙半径的分布曲线逐渐向左偏移,说明小孔隙所占比例不断增加,大孔隙所占比例逐渐减小,该结论与李荣强等[24]的实验结果一致。分析其原因在于有效应力的增大使得孔隙空间不断压缩,原本体积较大的孔隙被挤压成小孔隙,且岩心孔隙度越大,变化趋势越明显[图6(c)],孔隙大小分布越均匀,平均孔隙度越小,其变化趋势越明显[图6(d)]。

3.2 孔隙面积变化比例与孔隙面积的关系

在有效应力作用下,岩心中各个孔隙面积的变化是导致其应力敏感性的直接原因。为了研究不同孔隙条件下影响孔隙结构变化的因素,将有效应力为20MPa时岩心中各孔隙面积的变化比例与其对应的初始孔隙面积及孔隙形状因子绘图,如图7图8所示。就单独某种孔隙结构的岩心而言,当孔隙初始面积较小时,孔隙面积的减小比例在较大范围内分布,但是随着孔隙初始面积的增大,孔隙面积的减小比例都比较大。在整个孔隙面积分布范围内,两者存在一定的正相关关系,即孔隙初始面积越大其变化比例越大。对比不同岩心可以看出,随着孔隙度的增加和孔隙大小越均匀,两者的正相关关系愈加明显,当孔隙度为10%[图7(c)]和孔隙大小很均匀[图7(d)]时,最大孔隙减小比例都接近90%。
图7 各孔隙面积变化比例与孔隙初始面积大小的关系

Fig.7 Effect of the initial pore area on its change ratio

图8 各孔隙面积变化比例与孔隙形状因子的关系

Fig.8 Effect of the pore shape factor on the change ratio of pore area

3.3 孔隙面积变化比例与孔隙形状因子的关系

对于孔隙形状因子而言,圆形和正方形结构的形状因子分别为0.079 6和0.062 5,三角形的形状因子范围为0~0.048 1,可见形状因子越小,孔隙形状越不规则。5种岩心中孔隙的形状因子大部分分布在0.05以下,并有少部分在0.05~0.06之间,可以看出其孔隙以三角形和正方形为主。同时,孔隙面积减小比例与其形状因子具有明显的负相关关系,即孔隙形状因子越小,孔隙越不规则越容易发生较大的形变,其应力敏感性越强,这种关系随着孔隙度的增加和孔隙大小变均匀都有一定的加强。

3.4 岩心固有渗透率变化

不同孔隙结构的岩心固有渗透率随有效应力的变化如图9图10所示。通过模拟发现,岩心渗透率随着孔隙度的减小会大幅度降低,孔隙度为10%的岩心初始渗透率为0.024×10-3μm2,当孔隙度减小至1%时,其初始固有渗透率仅为0.000 61×10-3μm2。随着有效应力的增加,岩心的固有渗透率在不断减小,但减小幅度逐渐减缓。初始渗透率较小的岩心,渗透率降低幅度越大,孔隙度为1%的岩心在有效应力为20MPa时,固有渗透率下降了58%。
图9 不同孔隙度条件下固有渗透率随有效应力的变化

Fig.9 Variation of permeability with the effective stress under different porosity

图10 不同孔径比例下固有渗透率随有效应力的变化

Fig.10 Variation of permeability with the effective stress under different radius ratio of pore

对于孔隙大小均匀性不同的岩心,在相同的孔隙度条件下,孔隙大小参差不齐的岩心由于大孔隙的存在使得其初始固有渗透率较大,约为0.11×10-3μm2。虽然在上述研究中证明有效应力作用下大孔隙变化比例较大,但依然具有相当导流能力,其渗透率降低幅度较小。但是与孔隙度相比,岩心中孔隙的均匀性对固有渗透的影响相对较小,三者的渗透率下降比例基本都为50%左右。

4 结论

本文以四参数随机生长法构造的数字岩心为基础,利用弹性力学理论分析了不同孔隙结构岩心的应力敏感性及其影响因素,得到如下结论:
(1)岩心基质在有效应力作用下使得孔隙面积减小以及孔隙位置迁移,并由此导致的渗流通道减小及迂曲度增大是导致页岩应力敏感的直接原因。
(2)有效应力作用下岩心中孔隙的缩小程度受孔隙的初始面积及其形状因子的影响,并与孔隙初始面积呈正相关关系,与孔隙形状因子呈负相关关系,相关程度随孔隙度增大和孔隙均匀性增强而增大。
(3)页岩岩心的固有渗透率随有效应力的增加而呈放缓性减小,固有渗透率的减小程度随岩心孔隙度的减小和孔隙大小均匀性的增强而增大。
1
Zou Caineng , Zhang Guosheng , Yang Zhi ,et al .Geological concepts,characteristics,resource potential and key techniques of unconventional hydrocarbon: On unconventional petroleum geology[J]. Petroleum Exploration and Development, 2013, 40(4): 385-399, 454.

邹才能, 张国生, 杨智, 等 . 非常规油气概念、特征、潜力及技术——兼论非常规油气地质学[J]. 石油勘探与开发, 2013, 40(4): 385-399, 454.

2
Zou Caineng , Zhu Rukai , Wu Songtao , et al . Types, characteristics, genesis and prospects of conventional and unconventional hydrocarbon accumulations: Taking tight oil and tight gas in China as an instance[J]. Acta Petrolei Sinica, 2012, 33(2):173-187.

邹才能, 朱如凯, 吴松涛, 等 . 常规与非常规油气聚集类型、特征、机理及展望——以中国致密油和致密气为例[J]. 石油学报, 2012,33(2):173-187.

3
Zhu Weiyao , Ma Dongxu , Zhu Huayin , et al . Stress sensitivity of shale gas reservoir and its influence on productivity[J]. Natural Gas Geoscience, 2016, 27(5): 892-897.

朱维耀, 马东旭, 朱华银, 等 . 页岩储层应力敏感性及其对产能影响[J]. 天然气地球科学2016, 27(5): 892-897.

4
Tian Leng , Li Hongfan , Ma Jixiang ,et al .Multi-stage and multi-layer percolation model of tight gas reservoir based on pressure gradient and stress sensitivity[J]. Natural Gas Geoscience, 2017, 28(12): 1898-1907.

田冷, 李鸿范, 马继翔, 等 . 基于启动压力梯度与应力敏感的致密气藏多层多级渗流模型[J]. 天然气地球科学, 2017, 28(12): 1898-1907.

5
Chalmers G R L , Ross D J K , Bustin R M . Geological controls on matrix permeability of Devonian Gas Shale in the Horn River and Liard basins,northeastern British Columbia,Canada[J].International Journal of Coal Geology,2012,103:120-131.

6
Bustin A M M , Bustin R M , Cui Xiaojun . Importance of fabric on the production of gas shale[C]//SPE Unconventional Reservoirs Conference. Colorado: Society of Petroleum Engineers, 2008:1-29.

7
Kwon O , Kronenberg A K , Gangi A F , et al . Permeability of Wilcox shale and its effective pressure law[J]. Journal of Geophysical Research:Solid Earth,2001,106(B9):19339-19353.

8
Guo Wei , Xiong Wei , Gao Shusheng . Experimental study on stress sensitivity of shale gas reservoir[J]. Special Oil and Gas Reservoirs, 2012, 19(1): 95-97.

郭为, 熊伟, 高树生 . 页岩气藏应力敏感效应实验研究[J]. 特种油气藏, 2012, 19(1): 95-97.

9
Zhang Rui , Ning Zhengfu , Yang Feng , et al . Experimental study on microscopic pore structure controls on shale permeability under compaction process[J]. Natural Gas Geoscience, 2014, 25(8): 1284-1289.

张睿, 宁正福, 杨峰, 等 . 微观孔隙结构对页岩应力敏感影响的实验研究[J]. 天然气地球科学, 2014, 25(8): 1284-1289.

10
Zhang Rui , Ning Zhengfu , Zhang Haishan ,et al . New insights and discussion on stress sensitivity of fractured tight reservoir[J]. Natural Gas Geoscience, 2016, 27(5): 918-923.

张睿, 宁正福, 张海山, 等 . 裂缝性致密储层应力敏感机理新认识[J]. 天然气地球科学, 2016, 27(5): 918-923.

11
Zhang Rui , Ning Zhengfu , Yang Feng ,et al . Shale stress sensitivity experiment and mechanism[J].Acta Petrolei Sinica,2015, 36(2): 224-231, 237.

张睿, 宁正福, 杨峰, 等 . 页岩应力敏感实验与机理[J]. 石油学报, 2015, 36(2): 224-231, 237.

12
Xiao Wenlian , Li Tao , Li Min , et al . Evaluation of the stress sensitivity in tight reservoirs[J]. Petroleum Exploration and Development, 2016, 43(1): 107-114.

肖文联, 李滔, 李闽, 等 . 致密储集层应力敏感性评价[J]. 石油勘探与开发, 2016, 43(1): 107-114.

13
Lei Gang , Wang Hao , Dong Pingchuan ,et al . Quantitative analysis on stress sensitivity of heterogeneous tight sandstone[J]. Petroleum Geology and Recovery Efficiency, 2015, 22(3): 90-94.

雷刚, 王昊, 董平川, 等 . 非均质致密砂岩应力敏感性的定量表征[J]. 油气地质与采收率, 2015, 22(3): 90-94.

14
Li Chuanliang . A theoretical formula of stress sensitivity index with compressibility of rock[J]. Lithologic Reservoirs, 2007, 19(4): 95-98.

李传亮 . 岩石应力敏感指数与压缩系数之间的关系式[J]. 岩性油气藏, 2007, 19(4): 95-98.

15
Lymberopoulos D P , Payatakes A C . Derivation of topological, geometrical, and correlational properties of porous media from pore-chart analysis of serial section data[J]. Journal of Colloid and Interface Science, 1992, 150(1): 61-80.

16
Varslot T , Ghous A , Latham S , et al . Pore scale characterization of carbonates at multiple scale: integration of Micro CT, BSEM and FIBSEM[J].Petrophysics,2010,51(6): 379-387.

17
Zhu Yihua , Tao Guo . Sequential indicator simulation technique and its application in 3D digital core modeling[J]. Well Logging Technology, 2007, 31(2): 112-115.

朱益华, 陶果 . 顺序指示模拟技术及其在3D数字岩心建模中的应用[J]. 测井技术, 2007, 31(2): 112-115.

18
Zhao Xiucai , Yao Jun , Tao Jun , et al . A method of constructing digital core by simulated annealing algorithm[J]. Applied Mathematics A Journal of Chinese Universities: Serial A, 2007, 22(2): 127-133.

赵秀才, 姚军, 陶军, 等 . 基于模拟退火算法的数字岩心建模方法[J]. 高校应用数学学报: A辑, 2007, 22(2): 127-133.

19
Li Renmin , Liu Songyu , Fang Lei , et al . Micro-structure of clay generated by quartet structure generation set[J]. Journal of Zhejiang University: Engineering Science, 2010, 44(10): 1897-1901.

李仁民, 刘松玉, 方磊, 等 . 采用随机生长四参数生成法构造黏土微观结构[J]. 浙江大学学报: 工学版, 2010, 44(10): 1897-1901.

20
Shen Linfang , Wang Zhiliang , Li Shaojun . Numerical simulation for seepage field of soil based on mesoscopic structure reconfiguration technology[J]. Rock and Soil Mechanics, 2015, 36(11): 3307-3314.

申林方,王志良,李邵军 . 基于土体细观结构重构技术的渗流场数值模拟[J]. 岩土力学, 2015, 36(11): 3307-3314.

21
Morrow N R . Interfacial Phenomena in Petroleum Recovery[M]. New York: Marcel Dekker Inc, 1991.

22
Wei Yuanlong , Yang Chunhe , Guo Yintong , et al . Experimental research on deformation and fracture characteristics of shale under cyclic loading[J]. Chinese Journal of Geotechnical Engineering, 2015, 37(12): 2262-2271.

魏元龙, 杨春和, 郭印同, 等 . 三轴循环荷载下页岩变形及破坏特征试验研究[J]. 岩土工程学报, 2015, 37(12): 2262-2271.

23
Costa A .Permeability-porosity relationship:A reexamination of the Kozeny-Carman equation based on a fractal pore-space geometry assumption[J].Geophysical Research Letters, 2006, 33(2): L02318.

24
Li Rongqiang , Gao Ying , Yang Yongfei ,et al . Experimental study on the pressure sensitive effects of cores based on CT scanning[J]. Petroleum Drilling Techniques, 2015, 43(5): 37-43.

李荣强, 高莹, 杨永飞, 等 . 基于CT扫描的岩心压敏效应实验研究[J]. 石油钻探技术, 2015, 43(5): 37-43.

文章导航

/