基于页岩油水两相渗流特性的油井产能模拟研究.pdf

doi:10.11911/syztjs.2023084引用格式：孙鑫，刘礼军，侯树刚，等.基于页岩油水两相渗流特性的油井产能模拟研究 J.石油钻探技术，2023,51（5）：167-172.SUNXin,LIULijun,HOUShugang,etal.Numericalsimulationofshaleoilwellproductivitybasedonshaleoil-watertwo-phaseflowcharacteristicsJ.PetroleumDrillingTechniques，2023,51(5)：167-172.基于页岩油水两相渗流特性的油井产能模拟研究孙鑫1,2，刘礼军3，侯树刚1，戴彩丽2，杜焕福1，王春伟1（1.中石化经纬有限公司地质测控技术研究院,山东青岛266003；2.中国石油大学（华东）石油工程学院,山东青岛266580；3.成都理工大学能源学院,四川成都610059）摘要:页岩孔隙结构及固液相互作用复杂，其微观渗流特性加大了页岩油产能预测的难度。为准确评价体积压裂后多尺度孔隙结构发育的页岩油藏产能，基于页岩储层油水两相相渗计算方法和嵌入式离散裂缝模型，考虑页岩真实孔隙结构作用下的微观油水两相渗流特性，形成了考虑页岩体积压裂页岩油藏产能的数值模拟方法。基于页岩储层孔径分布计算油水相渗曲线，结合页岩油藏压裂/生产流程，开展了页岩油藏压裂液空间分布以及油井产能评价模拟分析。结果表明，不同孔径分布下的页岩油水两相相渗曲线存在差异，压裂液主要分布在压裂裂缝、与其相连的天然裂缝以及其周边基质中，在闷井过程中裂缝内压裂液逐渐渗吸进入基质并置换基质中原油，经体积压裂可实现改造区域的整体动用。研究结果可以从微观油水两相渗流特性与宏观产能评价角度为页岩油藏高效开发提供技术支撑。关键词:页岩油藏；水力压裂；油水相渗；压裂液分布；油井产能；数值模拟中图分类号:TE32+8文献标志码:A文章编号:10010890（2023）05016706Numerical Simulation of Shale Oil Well Productivity Based onShale Oil-Water Two-Phase Flow CharacteristicsSUN Xin1,2,LIU Lijun3,HOU Shugang1,DAI Caili2,DU Huanfu1,WANG Chunwei1(1.Geosteering&Logging Research Institute,Sinopec Matrix Corporation,Qingdao,Shandong,266003,China;2.School ofPetroleum Engineering,China University of Petroleum(East China),Qingdao,Shandong,266580,China;3.College of Energy,Chengdu University of Technology,Chengdu,Sichuan,610059,China)Abstract:Theporestructureofshaleiscomplex,andsolid-liquidinteractionoccurs.Inaddition,itsmicroscopicflow characteristics increase the difficulty of shale oil productivity prediction.In order to accurately evaluate theproductivityofshaleoilreservoirswithmulti-scaleporestructuresaftervolumefracturing,themicroscopicmulti-phaseflowcharacteristicsundertheactionofrealporestructureofshalewereconsideredbasedontheoil-watertwo-phaserelativepermeabilitycalculationmethodandembeddeddiscretefracturemodel(EDFM)ofshalereservoirs.Asaresult,anumericalsimulationmethodforshaleoilreservoirproductivityconsideringshalevolumefracturingwasdeveloped.Theoil-watertwo-phaserelativepermeabilitycurvewascalculatedbasedontheporesizedistributionofshalereservoirs,andcombinedwiththefracturing/productionprocessofshalereservoirs,thespatialdistributionoffracturingfluidinshaleoilreservoirsandtheproductivityevaluationofoilwellweresimulatedandanalyzed.Theresultsshowthattherearedifferencesintheoil-watertwo-phaserelativepermeabilitycurvesofshaleunderdifferentporesizedistributions.Fracturingfluidsaremainlydistributedinfracturingfractures,naturalfracturesconnectedwiththem,andthesurroundingmatrix.Duringtheprocessofshut-in,thefracturingfluidinthefractureisgraduallyimbibedintothematrix,displacingthecrudeoilinthematrixandrealizingthewholeutilizationofthestimulatedareaby收稿日期:2023-05-17；改回日期:2023-08-01。作者简介:孙鑫（1993），男，山东昌邑人，2016 年毕业于中国石油大学（华东）石油工程专业，2021 年获中国石油大学（华东）油气田开发工程专业博士学位，博士后，主要从事非常规储层评价与产能预测技术研究工作。E-mail：。通信作者:刘礼军，。基金项目：四川省自然科学基金项目“基于流固耦合数值模拟的陆相页岩凝析气藏合理开发方式探索”(编号：2022NSFSC1077)、青岛市博士后项目“东营凹陷页岩油测录井产能预测方法研究”(编号：QDBSH20220201023)、中石化经纬有限公司博士后研究项目“东营凹陷页岩油测录井产能预测方法研究”(编号：JWBH2203)联合资助。第51卷第5期石油钻探技术Vol.51No.52023年9月PETROLEUMDRILLINGTECHNIQUESSep.,2023volume fracturing.The research results can provide technical support for the efficient development of shale oilreservoirsfromtheperspectiveofmicroscopicoil-watertwo-phaseflowcharacteristicsandmacroscopicproductivityevaluation.Key words:shaleoilreservoir;hydraulicfracturing;oil-waterrelativepermeability;distributionoffracturingfluid;wellproductivity;numericalsimulation我国油气对外依存度持续上升，常规油气田进入开发后期，页岩油气等非常规油气资源的重要性日益显著12。2022 年，中国石油的页岩油产量突破 300104t，页岩油作为油气资源的后起之秀，其规模化开发正加速推进3。作为实现页岩油藏规模效益开发的关键技术，水平井和体积压裂技术受页岩复杂孔隙结构和固液相互作用影响45，目前在页岩油渗流机理、压裂液返排规律及产能预测等方面仍面临诸多问题与挑战。页岩油藏压裂过程中，高压泵注的压裂液使主裂缝和次生裂缝延伸，形成复杂人工裂缝网络；同时，压裂液会通过渗吸作用进入并滞留在基质中，返排过程中少量排出，影响页岩油后续产能68。实践表明，不同页岩油藏储层的压裂液返排率与产能差异较大，其返排特征主要受页岩储层中油水两相渗流特性影响，由页岩孔隙结构特征控制。页岩油储层孔隙结构通常具有纳米级孔隙发育、孔径分布范围广的特点，前人基于页岩孔隙内流体流动规律、孔径分布等提出了多种两相相对渗透率计算方法。WangJinxun 等人9将多孔介质概念化为不同尺寸管道串并联的毛细管模型，考虑孔隙尺寸分布及孔隙形状，推导了储层相渗的经典计算方法。LiRan 等人10综合考虑单孔内两相流动特征和孔隙分形结构，建立了页岩两相相渗计算方法，并分析了孔隙尺寸和结构对相渗特征的影响。SuYuliang 等人11考虑页岩有机和无机孔隙内油水赋存特征，建立了页岩油水两相相渗计算方法。数值模拟是常用的油藏产能评价手段，但体积压裂后的页岩油藏多尺度孔隙和裂缝储渗空间发育1213，传统双重介质模型、局部网格加密等模拟方法具有计算量大、难以处理复杂结构裂缝等局限，同时，页岩油流动受多种机理影响1417，数值模拟难度大。嵌入式离散裂缝模型是将复杂裂缝几何形态直接嵌入正交背景网格中，简化了裂缝的几何剖分过程，极大地降低了计算量和复杂度18。国内外学者综合页岩油渗流机理和嵌入式离散裂缝模型，开展了页岩油井产能模拟分析1920，但目前尚未实现微观页岩油水相渗计算与宏观页岩油井产能的耦合。因此，笔者结合页岩油藏相渗计算方法、嵌入式离散裂缝模型和油水两相渗流数学模型，提出了考虑页岩孔隙结构作用下油水两相渗流特性的页岩油井产能数值模拟方法，分析了体积压裂后页岩油藏压裂液空间分布特征和油井产能，实现了页岩微观油水两相渗流特性与宏观油井产能的一体化评价。1数学模型及求解 1.1页岩油藏相渗计算方法基于毛细管相渗计算模型，结合实际页岩孔隙形状和孔径分布，建立页岩油藏油水两相相渗计算方程。考虑页岩储层中复杂的孔隙形状，采用三角形毛细管模型表征页岩油藏储层。根据三角形毛细管中油水分布状态，单个毛细管中的油水两相流动规律可表示为：Qo=r4eff8oLp（1）Qw=(14G)r4dwLp（2）式中：Qo和 Qw分别为毛细管中油相和水相的流量，m3/s；o和 w分别为油相和水相黏度，Pas；L 为毛细管长度，m；p 为施加在毛细管上的压差，Pa；为水与孔隙壁面相互作用的无因次阻力系数（用于表征孔隙表面性质对流体流动的影响）；reff为毛细管有效油相半径，m；G 为三角形毛细管的形状因子；rd为油水稳定状态下的界面曲率半径，m。reff=21vuutGP2(14G)r2d+1rin（3）rd=P0.5G+G（4）式中：P 为三角形毛细管截面周长，m；rin为毛细管内切圆半径，m。结合单个毛细管中油水流动规律和页岩孔径分布，可得页岩储层油水相对渗透率计算公式：168石油钻探技术2023年9月Kro,m=5mk=1fkr4eff,k6nk=1fkr2in,kAk（5）Krw,m=20mk=1fk(14G)r4d,k3nk=1fkr2in,kAk+nk=m+1fkr2in,kAknk=1fkr2in,kAk（6）Sw,m=mk=1fk(14G)r2d,k+14Gnk=m+1fkr2in,knk=1fkAk（7）式中：n 为不同尺寸孔隙总数，其中，1m 为中心含水的孔隙尺寸数量，m+1n 为边缘含水的孔隙尺寸数量；Ak为第 k 个尺寸的毛细管截面积，m2；fk为第k 个尺寸的毛细管所占比例。因此，已知页岩储层的孔径分布后，便可根据式（5）式（7）计算出页岩油水相对渗透率 1.2页岩油藏油水两相渗流数学模型考虑页岩油藏中的油水两相渗流过程，其基质和裂缝中油水两相流体质量守恒关系可统一表达为连续性方程。t(S)=(v)+q（8）式中：为 o 或 w，代表油相或水相；为孔隙度；为 相流体的密度，kg/m3；S为 相流体饱和度；q为 相流体的源汇项，kg/(m3s)；v为 相流体的渗流速度，m/s。考虑页岩中流体流动的最小启动压力梯度效应，用非线性渗流模型描述基质中油水两相流动21。裂缝中通常不存在启动压力梯度效应，因此采用常规达西定律描述裂缝内的油水两相流动过程：v=KKr11a+b?（9）=pgh（10）式中：K 为绝对渗透率，m2；Kr为 相流体的相对渗透率；b 为拟启动压力梯度的倒数，(Pa/m)1；a 为非线性渗流凹形曲线段的影响因子；为 相流体的流动势，Pa；p为 相流体压力，Pa；h 为深度，m。1.3页岩油藏渗流模型求解方法为了高效求解页岩油藏油水两相流体流动，基于嵌入式离散裂缝模型对体积压裂后页岩油藏中复杂裂缝进行几何离散和网格剖分（见图 1）。对于给定的体积压裂页岩油藏模型，采用结构化网格对基质区域进行剖分，将水力压裂缝和天然裂缝网络嵌入至剖分后的结构化网格中，利用结构化网格边界切割裂缝网络，形成离散裂缝网格单元，综合形成页岩油藏数值模拟的网格单元系统。基质基质基质基质单元裂缝单元基质基质基质裂缝裂缝裂缝裂缝裂缝图 1 嵌入式离散裂缝模型示意Fig.1 Embedded discrete fracture model基于网格单元系统，采用有限体积法对油水两相渗流模型进行数值离散，并推导得到离散方程的残差形式：Rt+1,i=ji()t+1ij+12Tt+1ij(t+1jt+1i)(1t+1ij)+(Vq)t+1i(VS)t+1i(VS)tit（11）ij+12式中：表示单元 i 和 j 界面上的加权平均；R,i为单元 i 中 相连续性方程的残差，kg/s；i为单元 i 的邻近单元集合；t+1 为当前时间步；t 为上一时间步；t 为当前时间步长，s；V 为单元体积，m3；为流度，定义为=Kr/，(Pas)1；Tij为单元 i 和 j 间的传导率，可分为基质和裂缝不同介质单元组合间的传导率2223；ij为启动压力梯度引起的附加阻力系数。ij=1a+b?t+1jt+1i?dij（12）式中：dij为单元 i 和 j 间的距离，m。采用牛顿-拉夫森方法求解离散的残差方程。cRt+1,i(xt+1k)xcxc,k+1=Rt+1,i(xt+1k)（13）xt+1k+1=xt+1k+xk+1（14）第51卷第5期孙鑫等.基于页岩油水两相渗流特性的油井产能模拟研究169式中：k 为迭代层次；c 为主变量向量元素；x 为主变量向量，选取油相压力和含水饱和度为主变量。在每个时间步中，采用上述求解格式进行迭代计算，并更新主变量至残差向量的范数小于设定的允许误差，进入下一个时间步进行计算。2页岩油藏压裂液分布及产能模拟分析 2.1页岩油藏油水相对渗透率计算选取单峰型孔径分布和双峰型孔径分布 2 种典型孔径分布页岩（见图 2），在孔隙形状参数相同的基础上，采用页岩油藏相渗计算方法计算油水相对渗透率，结果见图 3。由图 3 可知，相比于单峰型孔径分布，双峰型孔径分布的孔隙尺寸更大，油相的流动能力更强。因此，双峰型孔径分布的页岩储层油相相对渗透率更大，水相相对渗透率更小。00.010.020.030.040.050.060.07102030孔隙尺寸/nm单峰型双峰型405060孔径分布频率图 2 两种典型页岩孔径分布概率曲线Fig.2 Two typical probability curves for shale pore sizedistribution0.20.40.60.8水相相渗-单峰型油相相渗-单峰型水相相渗-双峰型油相相渗-双峰型1.0含水饱和度00.20.40.60.81.01.2相对渗透率图 3 两种孔径分布计算的油水两相相渗曲线Fig.3 Oil-water two-phase relative permeability curvescalculated with two pore size distributions 2.2页岩油藏压裂液分布特征分析为分析压裂过程中压裂液的分布特征，基于胜利油田某页岩油井地质及压裂设计资料，结合嵌入式离散裂缝模型，建立体积压裂页岩油藏模型（见图 4）。该页岩油藏基质孔隙度 7.0%，渗透率 0.5D；水力裂缝开度为 4mm，渗透率为 5D；天然裂缝开度为 0.3mm，渗透率 0.1D；初始油藏压力 40MPa，初始含水饱和度为 0.05，油相和水相的黏度分别为0.40 和 0.25mPas，压裂液注入量为 1.0104m3，压裂后闷井时间为 30d，油水相对渗透率曲线采用图 3中单峰型孔径分布的计算结果。注入压裂液后的页岩油藏基质、裂缝中的压力和含水饱和度分布模拟结果如图 5 所示，压裂结束闷井 30d 后的基质、裂缝中的压力和含水饱和度分布则如图 6 所示。水力裂缝天然裂缝370 m1 700 m图 4 体积压裂页岩油藏模型Fig.4 Shale oil reservoir model by volume fracturing（a）基质压力压力/MPa压力/MPa41 42 43 44 45 46 47 48 49 50 51 52 53 54 5541 42 43 44 45 46 47 48 49 50 51 52 53 54 550.052含水饱和度含水饱和度0.056 0.060 0.064 0.068 0.072 0.076 0.08000.10.20.30.40.50.60.70.80.91.0（b）基质含水饱和度（c）裂缝压力（d）裂缝含水饱和度图 5 压裂结束时页岩油藏压力和含水饱和度分布模拟结果Fig.5 Simulation results of pressure and water saturationdistribution in shale oil reservoir after fracturing（a）基质压力（b）基质含水饱和度（c）裂缝压力（d）裂缝含水饱和度压力/MPa 41 42 43 44 45 46 47 48 49 50 51 52 53 54 550.052含水饱和度0.056 0.060 0.064 0.068 0.072 0.076 0.080压力/MPa 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55含水饱和度 00.10.20.30.40.50.60.70.80.91.0图 6 闷井 30 d 后页岩油藏压力和含水饱和度分布模拟结果Fig.6 Simulation results of pressure and water saturationdistribution in shale oil reservoir after 30 days ofshut-in从图 5 和图 6 可以看出，压裂液主要进入压裂缝及其周边天然裂缝和基质中，引起近水力裂缝周170石油钻探技术2023年9月边区域压力升高，该区域裂缝内含水饱和度显著上升，近水力裂缝基质内含水饱和度有所提升。进入闷井阶段后，裂缝和基质中的压力逐渐向周围区域耗散，近水力裂缝高压区域内的压力逐渐降低。同时，裂缝内的压裂液在毛细管力作用下渗吸进入基质，裂缝内含水饱和度降低，对基质中的原油产生一定的渗吸置换作用。2.3页岩油藏压后产能分析在压裂液注入和闷井的模拟结果基础上，开展页岩油藏压后产能数值模拟，对页岩油藏衰竭开发动用范围和产油量进行评价。衰竭开发 1000d 后的储层基质和裂缝中的压力和含油饱和度分布如图 7 所示。经体积压裂后的页岩油藏裂缝网络发育，衰竭开发过程中油藏动用程度高，天然裂缝发育范围内基本可以动用开发。此外，衰竭开发后裂缝内含水量低，但基质内的含水饱和度分布与生产前差异不大，其原因在于压裂液在毛管力作用下滞留在基质中，生产压差难以克服毛管阻力，这也解释了实际页岩储层压裂后压裂液返排率低的现象。页岩油藏生产 1000d 的日产油量和累计产油量曲线如图 8 所示。体积压裂页岩油藏衰竭开发产量递减速度快，经 1000d 开发后油井累计产油量可达61145m3。此外，开发过程中累计产水量为 3335m3，忽略地层水产出，计算得出压裂液反排率为 33%，与现场实际基本符合。由此可见，受页岩油水两相渗流特性及毛管力作用影响，页岩储层压裂后压裂液返排率较低，但体积压裂后的页岩油藏动用程度较好。（a）基质压力（b）基质含水饱和度（c）裂缝压力（d）裂缝含水饱和度压力/MPa 1015202530350.052含水饱和度0.056 0.060 0.064 0.068 0.072 0.076 0.080压力/MPa 10152025303500.010.020.030.040.05含水饱和度图 7 生产 1 000 d 后页岩油藏压力和含水饱和度分布模拟结果Fig.7 Simulation results of pressure and water saturationdistribution in shale oil reservoir after 1000 days ofproduction3结论1）考虑页岩孔隙结构作用下油水两相渗流特性，压裂页岩油藏产能数值模拟方法可实现页岩油藏油水相对渗透率、压裂液分布和返排以及油井产能的全流程评价。2）基于页岩储层孔径分布以及孔隙结构参数，采用毛细管模型可得到页岩油藏油水相对渗透率，不同页岩孔径分布下油水相对渗透率存在较大差异。3）压裂液主要分布于压裂缝及其周边的天然裂缝和基质中，闷井阶段进入周边基质，基质毛管阻力作用导致压裂液返排率较低，但体积压裂后的页岩油藏动用范围和程度较好。参考文献References邹才能，马锋，潘松圻，等.全球页岩油形成分布潜力及中国陆相页岩油理论技术进展 J.地学前缘，2023，30（1）：128142.ZOUCaineng,MAFeng,PANSongqi,etal.Formationanddistri-butionpotentialofglobalshaleoilandthedevelopmentsofcontin-ental shale oil theory and technology in ChinaJ.Earth ScienceFrontiers,2023,30(1):128142.1杨雷，金之钧.全球页岩油发展及展望 J.中国石油勘探，2019，24（5）：553559.YANGLei,JIN Zhijun.Global shale oil development and pro-spectsJ.ChinaPetroleumExploration,2019,24(5):553559.2200400600时间/d（a）日产油量8001 00001234567200400600时间/d累计产油量/（104m3）02004006008001 0008001 000（b）累计产油量日产油量/m3图 8 生产 1 000 d 后页岩油藏日产油量和累积产油量曲线Fig.8 Daily oil production and cumulative oil productioncurves of shale oil reservoir after 1000 days of pro-duction第51卷第5期孙鑫等.基于页岩油水两相渗流特性的油井产能模拟研究171佚名.我国页岩油勘探开发在多领域获重要进展 J.天然气勘探与开发，2023，46（1）：118.Anon.ShaleoilexplorationanddevelopmentinChinahasmadeim-portantprogressinmanyfieldsJ.NaturalGasExplorationandDe-velopment,2023,46(1):118.3于学亮，胥云，翁定为，等.页岩油藏“密切割”体积改造产能影响因素分析 J.西南石油大学学报(自然科学版)，2020，42（3）：132143.YUXueliang,XUYun,WENGDingwei,etal.Factorsinfluencingtheproductivityofthemulti-fracturedshaleoilreservoirwithtight-erclustersJ.JournalofSouthwestPetroleumUniversity(Science&TechnologyEdition),2020,42(3):132143.4倪华峰，杨光，张延兵.长庆油田页岩油大井丛水平井钻井提速技术 J.石油钻探技术，2021，49（4）：2933.NIHuafeng,YANGGuang,ZHANGYanbing.ROPimprovementtechnologies for large-cluster horizontal shale oil wells in theChangqingOilfieldJ.PetroleumDrillingTechniques,2021,49(4):2933.5欧阳伟平，张冕，孙虎，等.页岩油水平井压裂渗吸驱油数值模拟研究 J.石油钻探技术，2021，49（4）：143149.OUYANGWeiping,ZHANGMian,SUNHu,etal.Numericalsim-ulation of oil displacement by fracturing imbibition in horizontalshale oil wellsJ.Petroleum Drilling Techniques,2021,49(4):143149.6王继伟，曲占庆，郭天魁，等.考虑压裂液渗吸的压后压裂液返排的数值模拟 J.深圳大学学报（理工版），2023，40（1）：5665.WANGJiwei,QUZhanqing,GUOTiankui,etal.Numericalsimu-lationofpost-fracturingflowbackconsideringfracturingfluidim-bibitionJ.JournalofShenzhenUniversity(Science&Engineering),2023,40(1):5665.7徐润滋，杨胜来，王吉涛，等.高温高压下陆相致密油藏非稳态压裂液渗吸机理研究 J.油气地质与采收率，2023，30（3）：94103.XURunzi,YANGShenglai,WANGJitao,etal.Studyonunsteadyimbibitionmechanismoffracturingfluidintightcontinentalreser-voirsunderhightemperatureandhighpressureJ.PetroleumGeo-logyandRecoveryEfficiency,2023,30(3):94103.8WANGJinxun,DONGMingzhe,YAOJun.CalculationofrelativepermeabilityinreservoirengineeringusinganinteractingtriangulartubebundlemodelJ.Particuology,2012,10(6):710721.9LIRan,CHENZhangxin,WUKeliu,etal.Afractalmodelforgas-waterrelativepermeabilitycurveinshalerocksJ.JournalofNatur-alGasScienceandEngineering,2020,81:103417.10SUYuliang,XUJilong,WANGWendong,etal.Relativepermeab-ilityestimationofoilwatertwo-phaseflowinshalereservoirJ.PetroleumScience,2022,19(3):11531164.11何利，肖阳，孙宜成，等.车 21 井区裂缝性油藏地质建模与工程设计一体化研究 J.特种油气藏，2021，28（5）：2329.HELi,XIAOYang,SUNYicheng,etal.Onintegrationofgeolo-gicalmodelingandengineeringdesignoffracturedoilreservoirsinWellBlockChe21J.SpecialOil&GasReservoirs,2021,28(5):2329.12孟勇，贾庆升，张潦源，等.东营凹陷页岩油储层层间干扰及裂缝扩展规律研究 J.石油钻探技术，2021，49（4）：130138.MENGYong,JIAQingsheng,ZHANGLiaoyuan,etal.Research13oninterlayerinterferenceandthefracturepropagationlawofshaleoilreservoirs in the Dongying SagJ.Petroleum Drilling Tech-niques,2021,49(4):130138.何吉祥，徐有杰，高阳，等.裂缝性致密油藏多级压裂水平井试井模型 J.断块油气田，2021，28（2）：241246.HEJixiang,XUYoujie,GAOYang,etal.Welltestmodelofmulti-stagefracturedhorizontalwellinfracturedtightreservoirsJ.Fault-BlockOil&GasField,2021,28(2):241246.14赵国翔，姚约东，王链，等.页岩油藏微尺度流动特征及应力敏感性分析 J.断块油气田，2021，28（2）：247252.ZHAOGuoxiang,YAOYuedong,WANGLian,etal.Microscaletransport behaviors of shale oil and stress sensitivity analysisJ.Fault-BlockOilandGasField,2021,28(2):247252.15李江，陈先超，高平，等.考虑应力敏感效应的裂缝性碳酸盐岩气井拟稳态产能预测方法 J.石油钻探技术，2021，49（3）：111116.LIJiang,CHENXianchao,GAOPing,etal.Apseudo-steady-stateproductivity prediction method for fractured carbonate gas wellsconsidering stress-sensitivity effectsJ.PetroleumDrilling Tech-niques,2021,49(3):111116.16李虹，于海洋，杨海烽，等.裂缝性非均质致密储层自适应应力敏感性研究 J.石油钻探技术，2022，50（3）：99105.LIHong,YUHaiyang,YANGHaifeng,etal.Adaptivestresssensit-ivitystudyoffracturedheterogeneoustightreservoirJ.PetroleumDrillingTechniques,2022,50(3):99105.17LEESH,LOUGHMF,JENSENCL.HierarchicalmodelingofflowinnaturallyfracturedformationswithmultiplelengthscalesJ.WaterResourcesResearch,2001,37(3):443455.18刘礼军，姚军，孙海，等.考虑启动压力梯度和应力敏感的页岩油井产能分析 J.石油钻探技术，2017，45（5）：8491.LIULijun,YAOJun,SUNHai,etal.Theeffectofthresholdpres-suregradientandstresssensitivityonshaleoilreservoirproductiv-ityJ.PetroleumDrillingTechniques,2017,45(5):8491.19BAIYuhu,LIULijun,FANWeipeng,etal.Coupledcompositionalflow and geomechanics modeling of fractured shale oil reservoirwithconfinedphasebehaviorJ.JournalofPetroleumScienceandEngineering,2021,196:107608.20杨清立，杨正明，王一飞，等.特低渗透油藏渗流理论研究 J.钻采工艺，2007，30（6）：5254.YANGQingli,YANGZhengming,WANGYifei,etal.Studyonflow theory in ultra-low permeability oil reservoirJ.Drilling&ProductionTechnology,2007,30(6):5254.21LIULijun,HUANGZhaoqin,YAOJun,etal.Simulatingtwo-phaseflowandgeomechanicaldeformationinfracturedkarstreservoirsbasedonacoupledhydro-mechanicalmodelJ.InternationalJourn-alofRockMechanicsandMiningSciences,2021,137:104543.22姚军，刘礼军，孙海，等.复杂裂缝性致密油藏注水吞吐数值模拟及机制分析 J.中国石油大学学报（自然科学版），2019，43（5）：108117.YAOJun,LIULijun,SUNHai,etal.Numericalsimulationandmechanismanalysisofwaterhuffandpuffprocessincomplexfrac-turedtightoilreservoirsJ.JournalofChinaUniversityofPetro-leum(EditionofNaturalScience),2019,43(5):108117.23编辑曹耐172石油钻探技术2023年9月