An Instantly Dissolving Fracturing Fluid Technology Using Hydrophobic Associating Polymers Based on Swelling-Curing Mechanisms
-
摘要:
针对疏水缔合聚合物配制压裂液时分散性能差、溶胀时间长和不适合现场连续混配的问题,综合机械与化学方法,研究了疏水缔合聚合物速溶压裂液技术。通过优化助溶剂加量,减小配液用水的溶度参数,利用负熵因素加快聚合物溶胀速率;同时,设计了同轴双向搅拌装置,提高了固液混合物紊流程度,缩短了聚合物熟化时间。试验结果表明:溶度参数与溶胀胶团直径是影响聚合物溶胀与熟化的关键因素;4%助溶剂可使疏水缔合聚合物的溶胀率在4 min内达到92%;搅拌速率对疏水缔合聚合物熟化时间的影响较小,搅拌方式是影响疏水缔合聚合物熟化速度的主要因素,同轴双向搅拌6 min时的溶胀率达到86%;形成的疏水缔合聚合物速溶压裂液在90 ℃条件下的抗剪切性能大于180 mPa·s,满足现场施工要求。疏水缔合聚合物速溶压裂液技术实现了聚合物压裂液连续混配,降低了聚合物压裂液残液、余液对环境的影响,为大规模体积压裂工厂化作业提供了技术支撑。
Abstract:Due to hydrophobic associating polymers fracturing fluid have problems, such as poor dispersing performance, a long swelling time and which was unsuitable for continuous mixing at the wellsite. Thus, an instantly dissolving fracturing fluid technology using hydrophobic associating polymer has been developed by comprehensively combining mechanical and chemical methods. The solubility parameter for fluid preparing was reduced through the development of a cosolvent, and the swelling time of polymer was accelerated by utilizing the negative entropy factor. The coaxial bi-directional stirring device was designed to improve the turbulent flow of solid-liquid mixture and to shorten the curing time of the polymer. Experiments showed that the solubility parameter and the swelling micelle diameter were the key factors affecting the swelling and curing of polymer. The swelling ratio of hydrophobic associating polymer reached 92% in 4 min by adding 4% cosolvent. We found that the stirring rate has little effect on the curing time of such polymer, and further, the stirring method is the main factor affecting the curing speed of polymer. The swelling ratio of the coaxial bi-directional stirring method reached 86% in 6 min, and the developed instant dissolving hydrophobic associating polymer fracturing fluid system have the shear resistance of greater than 180 mPa·s at 90 °C, which satisfies the field requirements. The new fracturing fluid technology could achieve the continuous mixing of polymer fracturing fluid, and reduce the influence of residual fluids on the environment, which provided technical support for the implementation of large-scale volumetric fracturing in factory mode.
-
近年来,页岩油气等非常规油气资源的勘探开发逐渐成为当前研究的热门领域之一,其规模化开采得益于储层改造技术的发展,压裂裂缝质量决定着压后产量的高低,而支撑剂的运移与铺置效果是评价人工裂缝有效性的关键[1–6]。粗糙壁面阻碍缝内支撑剂的高效运移与铺置[7–11],国内外学者从物理试验和数值模拟方面开展了大量研究,Wen Qingzhi等人[12–13]通过平板试验模拟了支撑剂在不同复杂程度裂缝中的运移规律,分析了不同敏感因素对支撑剂运移和沉积的影响;郭天魁等人[14]研制了大尺寸压裂裂缝模拟装置,还原了支撑剂在裂缝中的真实运移情况;U. A. Inyang等人[15–16]通过室内试验模拟了不同压裂工艺和施工参数下支撑剂的运移和展布情况。但是,由于试验仪器尺度、材质、参数和试验目的的差异较大,未形成统一的测试标准和流程,难以系统指导压裂设计。因此,基于流体动力学和流固耦合的数值模拟技术成为研究支撑剂运移与展布的另一重要手段。现阶段数值模拟方法主要为欧拉−欧拉法和欧拉−拉格朗日法,前者将固相和液相处理成连续介质,采用欧拉法描述流体与固体颗粒的运动过程[17–19],该方法计算量小,常用于模拟低浓度固体颗粒的运移过程,但该方法无法考虑每个颗粒的运动与作用,而欧拉−拉格朗日法能较好地解决上述问题,模型收敛性较好,计算速度快[20–22]。然而,当前关于粗糙裂缝对流体及颗粒流动行为的研究较少,且缺乏定量描述复杂结构面流动机制的研究。
为此,针对不同粗糙度裂缝壁面编写了自相关高斯曲面代码,建立了不同粗糙度裂缝模型,采用CFD-DEM数值模拟方法,模拟了颗粒与流体在粗糙裂缝中的运移和粗糙裂缝内支撑剂的铺置过程,分析了裂缝内支撑剂颗粒滞留位置、堆积结构以及覆盖面积等展布规律,以期指导压裂参数优化。
1. 粗糙裂缝面模型
1.1 模拟方法
采用CFD-DEM双向耦合方法研究粗糙裂缝壁面间的颗粒−流体混合物流动问题,其耦合过程及过程中的固−液运动如图1所示[23]。通过计算同一时刻流域的压力和速度场,将流体压力和动量转移到颗粒运动,在循环结束前,在计算时间t上增加一个流场时间增量Δt,形成新的计算时间t'继续计算,直至计算时间与设置总时间ttotal相等时,循环结束。
运动连续性和动量守恒方程为[24]:
\frac{\partial }{{\partial t}}({\varepsilon _{\text{g}}}\rho {}_{\text{g}}) + \nabla \cdot ({\varepsilon _{\text{g}}}\rho {}_{\text{g}}{\boldsymbol{u}_{\text{g}}}) = 0 (1) \begin{split} \frac{\partial }{\partial t}\left({\varepsilon }_{\text{g}}{\rho }_{\text{g}}{\boldsymbol{u}}_{\text{g}}\right)+\nabla \cdot \left({\varepsilon }_{\text{g}}{\rho }_{\text{g}}{\boldsymbol{u}}_{\text{g}}\right)=&-{\varepsilon }_{\text{g}}\nabla p+\nabla \cdot \left({\varepsilon }_{\text{g}}{\boldsymbol{\tau} }_{\text{g}}\right)+ \\& {\varepsilon }_{\text{g}}{\rho }_{\text{g}}{g}-{\boldsymbol{F}}_{\text{g}-\text{p}} \end{split} (2) 考虑颗粒与颗粒、颗粒与壁面之间的碰撞,得到离散项控制方程:
{m_i}\frac{{{\text{d}}{{{u}}_i}}}{{{\text{d}}t}} = {m_i}g + {f_{{\text{p}} - {\text{g}},i}} + \sum\nolimits_{j = 1}^{{k_i}} {{f_{{\text{c}},ij}}} (3) {I_i}\frac{{{\text{d}}{w_i}}}{{{\text{d}}t}} = \sum\nolimits_{j = 1}^{{k_i}} {{T_{ij}}} (4) 颗粒−流体的曳力模型可表示为[25]:
{F_{\text{d}}} = \frac{{\beta ({u_{\text{g}}} - {u_{\text{p}}})}}{{{\rho _{\text{g}}}}} (5) \beta =\left\{\begin{array}{l}150\dfrac{(1-{\varepsilon }_{\text{g}})2{\mu }_{\text{g}}}{{\varepsilon }_{\text{g}}{d}_{\text{p}}{}^{2}}+1.75\dfrac{(1-{\varepsilon }_{\text{g}}){\rho }_{\text{g}}}{{d}_{\text{p}}}\left|{u}_{\text{g}}-{u}_{\text{p}}\right|\\ \dfrac{3}{4}{C}_{\text{D}}\dfrac{(1-{\varepsilon }_{\text{g}}){\rho }_{\text{g}}}{{d}_{\text{p}}}\left|{u}_{\text{g}}-{u}_{\text{p}}\right|{\varepsilon }_{\text{g}}{}^{-2.65}\\ {\varepsilon }_{\text{g}}>0.8\end{array}\right. (6) C_{\text{D}}=\left\{\begin{array}{ll}{\scriptstyle \dfrac{24}{Re}}(1+0.15Re^{0.687})&Re\le 1\;000\\ 0.44 & Re>1\;000\end{array}\right. (7) 式中:Fd为流体施加的曳力,N;β为曳力方程;CD为曳力系数;ρg为连续相流体密度,kg/m3;ug为流体运动速度,m/s;p为流体压力,Pa;εg为流体的体积分数;τg为流场内流体黏性应力张量;g为重力加速度,m/s2;mi为颗粒质量,kg;fp-g是颗粒相与流体相间作用力,N;fc为颗粒间接触力,N;Ii为颗粒转动惯量距,kg·m2;wi为颗粒的旋转角速度,rad/s;Tij为接触颗粒切向力力矩,N·m;up为颗粒流场速度;m/s;dp为颗粒直径,m;μg为流体黏度,Pa·s; Re为雷诺数;i,j分别为第i、j个颗粒。
1.2 模型建立
采用粗糙度系数CR描述结构面粗糙程度[26],CR与分形维数D之间的经验关系式为[27]:
\mathrm{ }C\mathrm{_R}=85.267\ 1(D-1)^{0.567\ 9} (8) 利用分形插值模拟岩石裂缝表面实际粗糙形态,对粗糙面尺寸进行高斯自相关函数变换:
C_{\text{h}}(u,v)=\sigma_{\text{h}}^2\text{exp}\left[-\left(\frac{u}{l_{\mathrm{c}1}}\right)^2-\left(\frac{v}{l_{\mathrm{c}2}}\right)^2\right] (9) 式中:CR为粗糙度系数;Ch(u,υ)为自相关函数;σh2为总体粗糙度的平方;u,υ分别为x,y方向上的滞后;lc1和 lc2是x,y方向上的相关长度系数。
通过自相关高斯函数得到不同粗糙度的点云数据,在此基础上建立三维曲面,通过盒维数法计算出处理后的二值图像,得到该曲面的分形维数D(见图2),代入式(8),求出曲面的粗糙度系数CR。
1.3 边界条件
假设模型边界条件与实际裂缝中流动方式一致,采用平面注入混相方式,定义近井端左侧为入口,右侧为出口,y轴负方向为重力方向。DEM采用不同平均粒径支撑剂颗粒,静摩擦因数为0.50,滚动摩擦因数为0.40,颗粒质量流速为0.02 kg/s,边界条件如图3所示。
采用结构化网格划分方法监测不同网格数量下同一高度沿y轴不同位置的湍流动能及绝对压力,将网格数量优化为10 500 个,如图4所示。
1.4 模型验证
为验证模型的准确性,采用文献[28]中的试验参数,建立380 mm×190 mm×2 mm平板模型,监测4 s内颗粒平均沉降速度。模型结果与试验结果间的平均误差为8.48 %,匹配性较好。误差成因主要在于部分颗粒进入分支裂缝,沉降速度偏低(见图5)。粗糙的壁面使颗粒运动受阻,颗粒运动不规则程度增加,沉降悬浮时间更长(见图6)。
2. 粗糙裂缝内支撑剂运移规律分析
2.1 支撑剂运移与沉降过程
支撑剂的沉降过程受多种因素影响,包括裂缝壁面粗糙度、支撑剂粒径、粒径组合方式等。颗粒在重力作用下沉降到裂缝底部,形成支撑剂床(见图6中蓝色区域)。随时间推移支撑剂床不断变高,其顶部上方空间变小,颗粒受到湍流作用影响,随着壁面粗糙度增大,缝道内流体湍流动能增大,扰动增加造成流道内支撑剂最大速度随之变大。当沉降时间为15 s时,颗粒沉降、摩擦和湍流效应达到平衡,此时支撑剂床高度为53 mm。此后已沉降支撑剂基本不再流动,继续泵入携砂液会造成砂堵。
2.2 壁面粗糙度
支撑剂在裂缝中的沉积高度和覆盖面积可反映支撑剂堵塞和裂缝的支撑效果,支撑剂悬浮区的形状也可表征支撑剂向深处的运移距离。在高粗糙度下,颗粒在缝口处的堆积高度更高,离缝口越远,颗粒沉积高度的下降越明显(见图7)。另外,压裂裂缝壁面越粗糙,支撑剂在裂缝中运移得越远,分布范围也会越广,模拟结果与Guo Tiankui等人[29]的研究结果一致。当裂缝粗糙度过大时(CR=76),支撑剂在裂缝曲率较大的区域形成卡槽和堵塞。
2.3 支撑剂粒径
为提高模型计算效率,采用支撑剂粒径与缝宽比(
\frac{{{d}}}{{{w}}} )来表征不同粒径通过缝宽的能力。选取CR=28、缝宽为1 mm作为初始条件,其余边界条件不变,分析不同\frac{{{d}}}{{{w}}} 条件下颗粒体积分数分布情况(见图8)。当\frac{{{d}}}{{{w}}} =0.3时,10 s后部分支撑剂发生沉降,形成多段凸起砂床,此时已沉降支撑剂在携砂液扰动下仍然具有流动性。当\frac{{{d}}}{{{w}}} =0.4时,支撑剂悬浮比例相对较大,在入口形成一定沉降,改善近井段充填效果。当\frac{{{d}}}{{{w}}} =0.8时,颗粒在缝道内的流动能形成固定通道(红色虚线),大尺寸支撑剂会提高近井端的支撑效果,但运移效率明显降低,造成入口短时间内堵塞,不利于支撑剂的高效铺置。通过监测得到相同质量流量下近井端已沉降支撑剂覆盖面积分别为41.7,83.3,175.0,195.0和215.0 cm2,覆盖面积显著增加,表明小粒径支撑剂携砂效率高,大粒径支撑剂铺置效果好。2.4 混合粒径组合
在边界条件一致的前提下,分析相同质量流量下不同粒径颗粒注入顺序对支撑剂铺置的影响。不同粒径颗粒在裂缝内的速度分布如图9所示。由图9可以看出:随着颗粒床高度增加,顶部颗粒速度较其他地方更大,颗粒堆积形态不规则,在裂缝曲面曲率较大的地方会形成间隙或堵塞;3种粒径支撑剂混合同时注入时,粒径0.4 mm支撑剂(黑色颗粒)主要呈悬浮态;粒径0.6 mm支撑剂(黄色颗粒)在沉积区底部呈层叠状分布;粒径0.8 mm支撑剂(红色颗粒)会导致近井端砂堤高度过高,后续支撑剂运移距离受限;“先小后大”的注入方式可在保证支撑剂整体运移效率的同时,增加后期井筒附近的支撑效果;“先大后小”的注入方式则受裂缝粗糙程度的影响,大粒径有加剧裂缝局部堵塞的风险。因此,对于岩性复杂或粗糙程度比较大的储层,在压裂过程中可先泵注加入小粒径支撑剂的携砂液,再泵注加入大粒径支撑剂的携砂液,以降低砂堵风险,提高输砂效率。
3. 结论与认识
1)粗糙度越大,流域内湍流动能和壁面剪切力在壁面处的变化越明显,这会影响支撑剂的运移和铺置。当粗糙度较大时,支撑剂易在缝内形成卡槽,且支撑剂在缝内的运移和铺置更不规则。
2)当粒径与缝宽比较小时,支撑剂在缝内不能形成稳定的支撑剂床,小粒径支撑剂在流体的作用下能运移得更远,大粒径支撑剂受重力影响短时间内在近井端沉积,运移距离短。
3)泵注混合粒径支撑剂时,大粒径支撑剂在缝道曲率较大区域形成间隙或卡槽效应,限制小颗粒支撑剂向深处运移,宜采用“先小后大”的泵注顺序。
-
表 1 压裂液破胶试验结果
Table 1 Experimental results of fracturing fluid
序号 胶囊破胶剂加量,% 过硫酸钠加量,% 不同破胶时间下的破胶状态 0.5 h 1.0 h 1.5 h 2.0 h 2.5 h 3.0 h 4.0 h 1 0.2 冻胶 冻胶 冻胶 冻胶 冻胶 冻胶 冻胶 2 0.2(碾碎) 冻胶 拉丝 轻微拉丝 破胶 破胶 破胶 破胶 3 0.8 冻胶 拉丝 轻微拉丝 轻微拉丝 基本破胶 基本破胶 基本破胶 4 1.0 冻胶 拉丝 轻微拉丝 基本破胶 破胶 破胶 破胶 5 1.5 冻胶 拉丝 基本破胶 破胶 破胶 破胶 破胶 注:胶囊破胶剂为过硫酸钠外层包裹一层胶囊的高温井用破胶剂,其主要破胶成分为过硫酸钠。 -
[1] 孙蒙. 疏水缔合物稠化剂的合成及评价[D]. 大庆: 东北石油大学, 2014. SUN Meng. The synthesis and evaluation of hydrophobic association thickener[D]. Daqing: Northeast Petroleum University, 2014.
[2] 肖兵, 范明福, 王延平, 等. 低摩阻超高温压裂液研究及应用[J]. 断块油气田, 2018, 25(4): 533–536. XIAO Bing, FAN Mingfu, WANG Yanping, et al. Study and application of low-friction high-temperature fracturing fluid[J]. Fault-Block Oil & Gas Field, 2018, 25(4): 533–536.
[3] 崔玉淼, 章敬, 李佳琦, 等. 玛湖凹陷白云质岩致密油藏改造技术研究与应用[J]. 新疆石油地质, 2014, 35(1): 101–104. CUI Yumiao, ZHANG Jing, LI Jiaqi, et al. Fracturing technology research and application on tight dolomite reservoir in Mahu Sag in Junggar Basin[J]. Xinjiang Petroleum Geology, 2014, 35(1): 101–104.
[4] 郭兴. 低粘压裂液体系及低粘压裂液用疏水缔合聚合物的制备与性能研究[D]. 西安: 陕西科技大学, 2017. GUO Xing. Preparation and property of low viscosity fracturing fluid and hydrophobically associating polymer for low viscosity fracturing fluid[D]. Xi’an: Shaanxi University of Science & Technology, 2017.
[5] 叶登胜, 王素兵, 蔡远红, 等. 连续混配压裂液及连续混配工艺应用实践[J]. 天然气工业, 2013, 33(10): 47–51. doi: 10.3787/j.issn.1000-0976.2013.10.007 YE Dengsheng, WANG Subing, CAI Yuanhong, et al. Application of continuously mixing fracturing fluid and such flow process[J]. Natural Gas Industry, 2013, 33(10): 47–51. doi: 10.3787/j.issn.1000-0976.2013.10.007
[6] 杜涛, 姚奕明, 蒋廷学, 等. 合成聚合物压裂液最新研究及应用进展[J]. 精细石油化工进展, 2016, 17(1): 1–5. doi: 10.3969/j.issn.1009-8348.2016.01.001 DU Tao, YAO Yiming, JIANG Tingxue, et al. Recent progress of research on synthetic polymer fracturing fluids and their application[J]. Advance in Fine Petrochemicals, 2016, 17(1): 1–5. doi: 10.3969/j.issn.1009-8348.2016.01.001
[7] YALKOWSKY S H, VALVANI S C. Solubility and partitioning I: solubility of nonelectrolytes in water[J]. Journal of Pharmaceutical Sciences, 1980, 69(8): 912–922. doi: 10.1002/jps.2600690814
[8] ASMUSSEN F, UEBERREITER K. Velocity of dissolution of polymers. Part II[J]. Journal of Polymer Science Part A: Polymer Chemistry, 1962, 57(165): 199–208.
[9] 唐康泰, 张自成. 极性聚合物聚丙烯酰胺的溶解Ⅰ[J]. 油田化学, 1988, 5(4): 249–257. TANG Kangtai, ZHANG Zicheng. Dissolution of polar polymer polyacrylamide: I[J]. Oilfield Chemistry, 1988, 5(4): 249–257.
[10] EZZELL S A, McCORMICK C L. Water-soluble copolymers. 39. synthesis and solution properties of associative acrylamido copoly-mers with pyrenesulfonamide fluorescence labels[J]. Macromole-cules, 1992, 25(7): 1881. doi: 10.1021/ma00033a007
[11] 钟传蓉. 疏水缔合丙烯酰胺共聚物的合成与性能及在溶液中的结构形态[D]. 成都: 四川大学, 2004. ZHONG Chuanrong. Study on the synthesis and properties of hydrophobically modified Polyacrylamides and its associating morphology in solution [D]. Chengdu: Sichuan University, 2004.
[12] 毛金成, 杨小江, 宋志峰, 等. 耐高温清洁压裂液体系HT-160的研制及性能评价[J]. 石油钻探技术, 2017, 45(6): 105–109. MAO Jincheng, YANG Xiaojiang, SONG Zhifeng, et al. Development and performance evaluation of hightemperature resi-stant clean fracturing fluid system HT-160[J]. Petroleum Drilling Techniques, 2017, 45(6): 105–109.
[13] 唐康泰, 张自成. 极性聚合物聚丙烯酰胺的溶解Ⅱ[J]. 油田化学, 1989, 6(1): 6–11. TANG Kangtai, ZHANG Zicheng. Dissolution of polar polymer polyacrylamide: Ⅱ[J]. Oilfield Chemistry, 1989, 6(1): 6–11.
[14] HILDEBRAND J H. Solubility of non-electrolytes[M]. 3rd ed. New York: Reinhold, 1936.
[15] 王刚, 孙刚, 刘洪兵. 超高分子量聚合物性能评价方法研究[J]. 大庆石油地质与开发, 2001, 20(2): 101–104. doi: 10.3969/j.issn.1000-3754.2001.02.033 WANG Gang, SUN Gang, LIU Hongbing. Extra high molecular polymer evaluation method[J]. Petroleum Geology & Oilfield Development in Daqing, 2001, 20(2): 101–104. doi: 10.3969/j.issn.1000-3754.2001.02.033
[16] 吴宗泽, 王序云, 高志. 机械设计[M]. 北京: 高等教育出版社, 2002: 38-55. WU Zongze, WANG Xuyun, GAO Zhi. Mechanical Design[M]. Beijing: Higher Education Press, 2002: 38-55.
[17] 刘景亮. 粘稠体系下同心双轴搅拌器气液分散和混合特性的研究[D]. 杭州: 浙江大学, 2015. LIU Jingliang. Gas dispersion and mixing characteristics of co-axial mixers in viscous systems[D]. Hangzhou: Zhejiang University, 2015.
-
期刊类型引用(20)
1. 马金龙. 浅层枯竭型气藏储气库钻井关键技术. 石油工业技术监督. 2024(02): 32-36 . 百度学术
2. 张思远,杨佳坤,宋丽娜,周栋梁,胡俊,许锋,施玉霞,胥洪成,裴根,范家屹. 储气库多周期运行井控参数预测模型. 油气藏评价与开发. 2024(05): 788-794 . 百度学术
3. 高小芃,赵虎,司西强. 近油基钻井液在文23储气库水平井的应用. 石油化工应用. 2024(10): 18-22 . 百度学术
4. 何斌斌,柳华杰,郑若臣,步玉环,殷慧,霍美桦,张军义,马小龙. 盐岩地层矿物离子对水泥浆性能的影响规律. 中国石油大学学报(自然科学版). 2024(06): 85-94 . 百度学术
5. 张强. 文23储气库储层段钻井液及储层保护技术. 断块油气田. 2023(03): 517-522 . 百度学术
6. 廖权文,胡建均,史怀忠,宋恒宇. 文23储气库钻井工程关键技术. 石油钻采工艺. 2023(02): 160-166 . 百度学术
7. 陈芳,马平平,杨立军,刘文超. 温西超低压储气库钻完井工程技术优化. 石油钻采工艺. 2023(02): 167-172 . 百度学术
8. 张赫. 永久式封隔器在储气库老井应用中的实践及认识. 石油和化工设备. 2023(08): 19-21 . 百度学术
9. 谢关宝,滕春鸣,柳华杰. 盐岩蠕变对水泥环气密封完整性影响规律研究. 石油钻探技术. 2022(02): 78-84 . 本站查看
10. 刘自广. 文23枯竭砂岩型储气库微泡钻井液技术. 钻探工程. 2022(02): 117-122 . 百度学术
11. 黄家根,段忠国,田霄,赵俊杰. 苏东39-61储气库钻井技术难点及对策. 石化技术. 2022(06): 54-56 . 百度学术
12. 路保平. 中国石化石油工程技术新进展与发展建议. 石油钻探技术. 2021(01): 1-10 . 本站查看
13. 曾大乾,张广权,张俊法,糜利栋,杨小松,贾跃玮,王丹丹,范照伟. 中石化地下储气库建设成就与发展展望. 天然气工业. 2021(09): 125-134 . 百度学术
14. 何祖清,何同,伊伟锴,孙鹏,艾爽. 中国石化枯竭气藏型储气库注采技术及发展建议. 地质与勘探. 2020(03): 605-613 . 百度学术
15. 袁光杰,张弘,金根泰,夏焱. 我国地下储气库钻井完井技术现状与发展建议. 石油钻探技术. 2020(03): 1-7 . 本站查看
16. 齐桂雪. 枯竭气藏型储气库储层应力敏感性实验研究. 石油地质与工程. 2020(03): 76-80 . 百度学术
17. 曾大乾,张俊法,张广权,糜利栋. 中石化地下储气库建库关键技术研究进展. 天然气工业. 2020(06): 115-123 . 百度学术
18. 秦文政,张彦华,高飞,孙超,党文辉. 呼图壁储气库水平井关键技术. 中国石油和化工标准与质量. 2020(15): 238-239+241 . 百度学术
19. 刘宏伟. 文23储气库老井安全封堵工艺设计及应用. 石油石化绿色低碳. 2020(05): 56-62 . 百度学术
20. 何祖清,伊伟锴,李军,孙鹏,蒋记伟,刘鹏林. 文23储气库注采管柱封隔器密封数值模拟研究. 石油管材与仪器. 2020(06): 47-51 . 百度学术
其他类型引用(4)