Digital Core Construction Methods for High Stress in Deep and Ultra-Deep Oil and Gas Reservoirs
-
摘要:
深层超深层油气藏由于埋藏深，其地应力达200 MPa，会显著改变储层岩石孔隙的微观结构。数字岩心是孔隙尺度数值模拟的重要载体，但是现有数字岩心重构方法是基于常温常压下岩心的扫描图像重构，不能反映高应力下的孔隙结构。为此，提出了一种基于离散元法考虑高应力影响的数字岩心重构方法。首先，采用分水岭算法分割CT图像，利用球面谐波分析方法建立轮廓数据库，并在PFC3D中建立Clump（团簇）模板库；然后，根据孔隙度和粒径分布使用模板库中的Clump建立离散元模型，并用两点相关和线性路径相关函数曲线评价模型的准确性；随后，标定颗粒间微观力学参数，并加载应力模拟得到不同应力下的数字岩心；最后，分析了不同应力下数字岩心的孔隙几何拓扑结构，计算孔隙度和渗透率。以Bentheim砂岩为例，构建了其不同应力下的数字岩心，研究结果表明，应力增大，导致孔隙和喉道半径缩小、喉道伸长、连通性变差、孔隙度和渗透率减小。研究结果为深层超深层油气藏孔隙尺度模拟提供了技术途径。
-
关键词:
- 油气藏 /
- 深层超深层；数字岩心重构 /
- 离散元法 /
- 孔隙几何拓扑结构 /
- 渗透率
Abstract:Deep and ultra-deep oil and gas reservoirs buried at significant depths and subjected to ground stresses of 200 MPa, undergo notable changes in the pore microstructure of reservoir rocks. Digital core modeling serves as a crucial tool for pore-scale numerical simulations. However, current digital core reconstruction methods are based on scanning image reconstruction under normal temperature and pressure conditions and thus they fail to reflect the pore structure under high pressure conditions. Therefore, a digital core reconstruction method based on the discrete element method (DEM) was proposed by considering the effect of high stress. Initially, the watershed algorithm was employed to segment computed tomography (CT) images, and the contour database was established by the spherical harmonic analysis method. The Clump template library was established in PFC3D. Then, according to porosity and particle size distribution, the Clump in template library was used to build a discrete element model. After, the accuracy of the model was evaluated via calculations of two-point correlation and linear path correlation function curves. Next, the micromechanical parameters between particles were calibrated, enabling simulation of the digital core under varying stress conditions. Finally, the pore geometry topology of the digital core under different stresses was analyzed, and porosity and permeability were calculated. Bentheim sandstone was taken as an example to construct digital cores under different stresses.. The research results show that high stress leads to reduced pore and throat radius, elongated throats, diminished connectivity, and lower porosity and permeability. The results provide technical support for pore-scale simulations of deep and ultra-deep oil and gas reservoirs..
-
油气井工程中，为了保持井筒稳定和地层–套管环空长效密封，通过注水泥固井封隔油、气、水层，为持续钻探、完井开采提供安全作业环境。固井工程是建井技术的核心环节，直接影响钻井安全、工程成本、油气井寿命甚至油气开采效率。中国陆上深层超深层地质条件复杂，尤其是塔里木盆地和四川盆地海相油气藏存在埋藏深（超8 000 m）、高温高压（温度超200 ℃）、富含酸性气体（如元坝气田H2S含量6.87%，CO2 含量15.51%）和分段改造作业压力高（深层页岩气压裂压力≥90 MPa）等问题，固井施工面临一系列技术难题[1-4]：
1）复杂地质环境及工况下的水泥环力学响应与损伤机制不清。对大型储层改造工况下的复杂交变动载，以及CO2和H2S等高酸性气体、流体冲蚀等复杂地层界面环境的影响机理研究不充分，导致水泥浆材料与体系的构建缺乏理论依据。
2）水泥浆功能材料缺乏，水泥环力化性能差，易脆性开裂、导致环空气窜。常规水泥石为脆性体，弹性模量高（高于10 GPa）、渗透率高（高于0.2 mD）和腐蚀速率高（高于1.0 mm/a），不能满足复杂油气井水泥环长效密封的要求。例如，涪陵页岩气田一期井口带压率高达70%。
3）尾管悬挂器承载不足，易腐蚀。常规尾管悬挂器承载低（低于1800 kN）、不耐腐蚀，深井超深井尾管悬挂器坐挂、丢手成功率低，常出现套管下放不到位的问题，难以满足超深井复杂环境需求。
4）固井设计与施工精细程度不足。常规固井关键参数设计依赖于经验，精准指导固井施工难度大，固井顶替效率差（低于65%），无法保证固井质量与施工安全。
针对上述技术难题，中国石化经过10余年技术攻关，在固井方法、材料与装置、工艺和技术等方面取得了突破，揭示了水泥环密封失效机理，研发了高性能固井水泥浆，研制了新型尾管悬挂器，提出了优化设计新方法。这些成果应用于现场后，解决了深层超深层油气井固井关键技术难题；同时，以“极端环境、绿色环保、智能跨越”为理念，持续进行传统固井技术升级迭代，有力支撑了我国深层超深层油气高效勘探开发。
1. 水泥环密封失效机理及控制方法
1.1 水泥石力学性能演化规律
常规水泥石单轴压缩、三轴压缩、直接/间接拉伸等力学试验所获取的抗压强度、弹性模量和抗拉强度等力学参数[5-6]，已无法满足深层超深层油气井水泥石力学、物化性能的定量设计的需要。为此，研制了全尺寸水泥石密封评价装置[7]，消除了尺寸效应，该装置耐温200 ℃、耐压150 MPa（见图1）；建立了水泥石密封能力定量评价方法，可模拟高温、高压、交变载荷和流体冲蚀等环境条件下的水泥石密封能力演化过程。
试验研究结果表明，高温高压下水泥石表现为高弹性模量（高于10 GPa）、脆性劈裂，在温压场耦合下产生显著的塑性变形，无明显宏观裂缝，但收缩显著。周期交变荷载作用下的水泥石应力–应变曲线如图2所示（σ为水泥石的应力，MPa；ε为水泥石的应变）：第1次加载卸载后应变不可完全恢复，产生残余应变；第2次加载卸载后继续产生残余应变；残余应变逐渐累积，多次加载卸载后残余应变率达0.60%～0.95%。研究表明，残余应变与循环次数、水泥石弹性模量及孔隙度呈正相关，主要由水泥石的孔壁坍塌和压密塑性变形导致。
![]()
图 2 交变荷载下水泥石应力–应变曲线Figure 2. Elastic-plastic constitutive model of set cement under alternate loads为此，建立了水泥石刚度退化模型、累积塑性应变与荷载水平及荷载周期的函数关系（见图2）；同时，推导了考虑损伤阈值的力学本构模型：
σ1={Eε1+μ(σ2+σ3)ε<εdEε1+μ(σ2+σ3)+(B1−Eε1){1−exp[−(αdI1+√J2−TdS0)n]}ε⩾ (1) 式中：σ1为水泥石最大主应力，MPa；σ2，σ3为水泥石次大和最小主应力，MPa；ε1为最大主应力方向的应变；εd为扩容屈服点处的应变；Td为水泥石损伤时强度，MPa；I1为球应变张量的第一不变量；J2为偏应力张量的第二不变量；αd，B1为水泥石达到扩容时的材料强度参数，由水泥石的黏聚力cd及内摩擦角φd确定；E为水泥石弹性模量，GPa；μ为水泥石泊松比；S0为水泥石宏观统计的平均强度参数；n为水泥石微元强度分布集中程度。
1.2 交变动载下水泥环密封失效机理及控制方法
揭示了页岩气井多级分段压裂动载下的水泥环密封失效机理[8-10]，即交变动载下套管–水泥环–地层界面变形不一致，当累积塑性应变超过界面极限胶结应变时，产生微环隙，导致密封失效、界面窜流；建立了套管–水泥环–围岩弹塑性有限元分步计算模型；另外，提出了“低弹模（5～7 GPa）、低渗透（渗透率0.001 mD）”防损伤水泥控制准则，创建了“压裂段数–泵压–弹性模量”耦合设计图版（见图3），指导了水泥浆新材料研发和新体系设计。
1.3 酸性环境下水泥环密封失效机理及控制方法
明确了高温条件下CO2和H2S耦合作用下水泥石结构疏松、膨胀开裂腐蚀机制[11-12]，即CO2、H2S基于酸碱反应腐蚀水泥石，消耗Ca(OH)2晶体，使水化硅酸钙凝胶（C−S−H）结晶化；H2S腐蚀生成二水石膏和钙矾石，导致水泥石膨胀开裂，CO2腐蚀生成Ca(HCO3)2，造成钙离子流失，强度和渗透率下降；CO2和H2S共同作用加剧了水化硅酸钙凝胶结晶化，存在成倍叠加效应，CO2是主导因素。
为了堵塞酸性气体渗入通道、减少碱性反应物质，提出了“低渗透（<0.001 mD）、低碱性晶体含量（<5%）”防腐蚀控制准则。水泥石中碱性晶体含量为5%时，腐蚀后的抗压强度为31.35 MPa；而碱性晶体含量为10%时，腐蚀后的抗压强度为28.14 MPa；随着碱性晶体含量降低，腐蚀的强度有所增加。同样，水泥石的渗透率越低，腐蚀后渗透率保持较低（见图4）。
![]()
图 4 水泥石碱性晶体含量和渗透率对腐蚀的影响规律Figure 4. Influence of alkali crystal content and permeability of set cement on corrosion建立了“降碱度、增惰性”水泥浆构建方法，通过添加硅铝网状聚合物活性材料或高价金属离子材料，消耗Ca(OH)2晶体，提高水化产物与酸性气体的反应惰性；同时，可添加胶乳类纳米弹性颗粒材料，提高水泥环的致密性。
1.4 地层–水泥界面劣化机理及控制方法
地层–水泥界面胶结容易失效，主要受井筒内钻井液滤饼、地层油气水、温度和压力等复杂地质工况的影响。
文献[13-14]发现了Ca(OH)2晶体富集与溶蚀是界面胶结劣化的主因：由于泌水效应，地层–水泥界面易富集大量Ca(OH)2晶体，在地层水环境下Ca(OH)2晶体溶蚀，导致界面胶结强度降低，易出现环空间隙，降低了防窜能力。为此，建立了以“界面等效渗透率”为目标函数、以“透射声波幅度”为主要变量的界面密封能力预测模型（见式（2）、式（3）），实现了对界面密封能力演化的定量评价与不同水泥浆体系的适应性分析；同时，提出了“晶体抑制、限位支撑”水泥浆控制方法，降低了界面溶蚀程度，增强了界面压实能力。
K_{\rm{eq}}={{\rm{e}}}^{-1.90-0.018u-1.6\times {10}^{-5}{u}^{2}} (2) u=\frac{{A}_{{\rm{E}}}}{{A}_{{\rm{EITC}}}}=\left(1-U\right)\times 95\% (3) 其中：Keq为界面等效渗透率，mD·mm2；u为透射首波信号幅度占比；AE为界面透射声幅值，mV；AEITC为仪器发射的声波总能量，mV；U为实测反射声波总能量占比。
2. 高性能水泥浆
2.1 高温高密度水泥浆
高温使水泥石强度大幅降低，同时降低水泥浆粘滞力，水泥浆体系发生沉降失稳，用常规铁矿粉加重的水泥浆存在“冷浆稠、高温稀”等问题。为了解决水泥流体和固体存在的问题，研制了高温强度衰退抑制剂，在高温245 ℃下28 d强度不衰退；开发了亚微米球形三氧化二铁加重材料，平均粒径160 nm，密度5.10～5.15 kg/L，该材料具有“滚动轴承”作用，能显著降低水泥浆体系的内摩擦力，冷浆稠度低，地面混配性能好，颗粒粒径小，不易发生沉降；研发了热增黏聚合物，高温条件下分子结构相态由网状转变为链状，水泥浆体系黏度增加，利于加重材料的均匀悬浮。基于紧密堆积理论，开发了高温高密度水泥浆，密度1.95～2.40 kg/L，耐温245 ℃，沉降密度差小于0.03 kg/L，48 h控压强度32 MPa。
2.2 高压防气窜水泥浆
针对高温高压气井渗透压差大、易窜流等问题，依据 “高致密、强胶结”防气窜准则，研制开发了纳米液硅和苯丙胶乳水泥浆体系[15-18]。
1）纳米液硅水泥浆。纳米液硅是一种纳米SiO2颗粒在水中分散悬浮的乳液，它以50～300 nm二氧化硅为原料，利用高能物理球磨分散、电位阻斥力等工艺制备而成[15]。乳液固相含量45%，球形颗粒，无定形态。纳米液硅在水泥浆中能改善水泥浆性能[16]：①限位支撑作用，降低水泥浆收缩率46%，界面胶结强度达2.6 MPa；②孔隙填充作用，水泥石渗透率0.007 mD；③增强作用，7 d强度53 MPa。纳米液硅水泥浆密度1.50～2.60 kg/L，耐温220 ℃，7 min静胶凝强度由48 Pa发展至240 Pa。
2）苯丙胶乳水泥浆。胶乳是一种高聚物粒子分散在水介质中的胶体分散体系，采用苯乙烯等为原料制备而成[17]，通过攻关研究，解决了乳液的高温稳定性。纳米苯丙乳液热裂解温度350 ℃。苯丙胶乳颗粒具有高弹体吸能作用，可将水泥石的弹性模量降至5.8 GPa；同时，它具有孔隙填充防气窜作用，室内试验表明，在水泥浆候凝“失重”时，即使气层压力高于液柱压力也未发生气窜 [18]。苯丙胶乳水泥浆密度1.85～2.60 kg/L，API滤失量32 mL，水泥石渗透率0.005 mD，耐温220 ℃。
2.3 防损伤自愈合水泥浆
自愈合乳液是在聚苯乙烯链段中加入疏水亲油链段，形成自修复分子主链和侧链，通过缠绕、自舒展实现膨胀。自愈合乳液固相含量44%，耐温160 ℃。自修复乳液在水泥石中以高弹乳膜存在，当水泥石受到载荷时，高弹颗粒发挥吸能缓冲作用，可降低水泥开裂、塑性变形和残余应变，遇烃后封堵水泥石裂隙，修复裂缝。乳液吸烃体积膨胀率达700%，自修复水泥浆耐温160 ℃，弹性模量为5.5 GPa，可修复宽度为350 μm的裂缝，可用于井深超4 500 m的页岩气井，能够解决深层页岩气井的井口带压难题[19]。
2.4 防腐增强型水泥浆
降碱增惰防腐材料以SiO2和Al2O3为原料，在激活剂作用下聚合而成。该材料可消耗Ca(OH)2碱性晶体，促进生成高强惰性产物5CaO·6SiO2·5H2O，降低水化产物与H2S和CO2的反应活性，降低水泥石的碱度和腐蚀速度，提高水泥石的致密性，使酸性环境下的水泥石强度稳定。防腐增强型水泥浆耐温220 ℃，渗透率小于0.001 mD，腐蚀速率小于0.1 mm/a，水泥石腐蚀后强度无明显衰退，确保了酸性气井水泥石的持久密封能力。
2.5 防漏泡沫水泥浆
在常规水泥浆或低密度水泥浆中加入高性能发泡剂、稳泡剂，利用自主研制的注氮装置及工艺，可形成泡沫水泥浆，其密度可通过调整注气量快速调节。泡沫水泥浆具有弹韧性强和可压缩等独特性能，密度为1.10～1.60 kg/L，水泥石抗压强度达16 MPa、弹性模量为4～8 GPa。在页岩油气井、超深井等的技术套管、生产套管固井中规模化应用[20-22]，大幅提高了固井质量，解决了固井低压漏失、浅层气窜与储层水泥环密封失效等技术难题（其中页岩气井水平段固井长度突破4 000 m）；同时，解决了川渝地区喀斯特地貌浅表层破碎带恶性漏、塌、垮并存的难题。
3. 超深井固井工具
3.1 高承载尾管悬挂器
塔里木盆地、四川盆地等地区多口超深井尾管长度超3 000 m、重量超3 000 kN，远高于常规尾管悬挂器的承载极限。常规尾管悬挂器是通过楔形卡瓦将尾管重力全部转化为对上层套管的外张力，从而实现尾管悬挂的，其承载方式为轴向–径向二维承载（见图5（a））。当尾管超重时，易使上层套管变形，导致尾管悬挂器坐挂失败。
![]()
图 5 常规承载机构与复合斜面重载机构对比Figure 5. Comparison between conventional bearing mechanism and compound-bevel heavy load mechanism高承载尾管悬挂器采用复合斜面重载机构，通过卡瓦侧面承载产生相互抵消的周向分力N，其承载方式为“轴向–径向–周向”三维承载（见图5（b）），从而显著降低对套管的外张力，承载能力较常规悬挂器提高1倍以上。
为此，研发了系列高承载尾管悬挂器产品[23]，其中ϕ365.1 mm×ϕ273.1 mm尾管悬挂器的最大承载能力达到3 400 kN，玉中2井创造了国内悬挂尾管最重纪录（3 100 kN）。
3.2 压力平衡式尾管悬挂器
深层超深层油气井具有裸眼段长、井眼复杂等特点，尾管下入遇阻和不到位的问题突出。为此，研发了系列压力平衡式尾管悬挂器[24-26]（见表1），其特有的双向自平衡液压驱动机构可消除大排量循环时的坐挂力，并且允许循环排量不受工具限制，从而防止下套管过程中尾管悬挂器提前坐挂，形成了下尾管中途循环解阻及尾管快速下入工艺和技术，工具耐温150 ℃，密封能力35 MPa。现场应用时，该压力平衡式悬挂器中途循环排量最高达2.6 m3/min、泵压达42 MPa，应用最大井深7 849 m，尾管最长3 509 m，解决了复杂油气井尾管下入慢和到位率低的难题。
表 1 压力平衡式尾管悬挂器的规格和性能参数Table 1. Specifications and performance parameters for pressure-balanced liner hangers序号 规格 最大外径/mm 最小内径/mm 承载能力/kN 坐挂压力/MPa 耐温能力/℃ 1 ϕ244.5 mm×ϕ177.8 mm 215.0 157.0 1 800 8～9 150 2 ϕ237.1 mm×ϕ193.7 mm 238.0 172.0 2 400 8～9 150 3 ϕ219.1 mm×ϕ168.3 mm 188.0 145.7 1 200 8～9 150 4 ϕ365.1 mm×ϕ273.1 mm 322.0 245.0 3 400 6～7 120 3.3 防腐蚀尾管悬挂器
针对尾管悬挂器在酸性环境下的腐蚀问题，根据尾管悬挂器结构特点和性能要求，建立了不同腐蚀环境下的零部件选材规范和评价准则，研制了防H2S、防CO2和防H2S/CO2等的系列防腐蚀尾管悬挂器[27]，满足了H2S、CO2含量大于20%的高酸性环境防腐蚀要求。为此，研发了耐H2S/CO2共同腐蚀的橡胶材料，在温度204 ℃、H2S分压3.5 MPa、CO2分压3.5 MPa条件下腐蚀老化后的关键性能保持率在70%以上。
3.4 封隔式分级注水泥器
针对常规分级注水泥器难以解决的超深井漏失问题，研发了封隔式分级注水泥器[28]，主要由水力扩张式封隔器和分级注水泥器组成（见图6）。一级注水泥后，憋压打开水力扩张封隔器的注液通道并填充胀封，封隔套管与井眼的环空，彻底封隔分级注水泥器下部漏层，从而降低二级固井漏失风险。同时，封隔器胀封后具有一定的环空封隔能力，并在分级箍循环孔关闭后能够保证管内密封能力不低于套管抗内压强度，提高了套管柱完整性。工具密封能力70 MPa，在10%井径扩大率下环空封隔能力大于35 MPa。
该封隔式分级注水泥器已经在塔里木盆地超深易漏失井应用了100余口井，二级固井无漏失率提高了90%。
4. 固井优化设计新方法
4.1 防气窜设计新方法
针对深层超深层油气井固井中的防气窜问题，开发了环空液柱压力精细控制技术，建立了环空液柱压力精细计算新模型[29-30]，预测精度较传统模型提高30%。实际应用时，可根据无量纲化实测井深、井径、地温梯度和水泥浆密度等参数，分段（段长1 m）拟合环空液柱压力。根据井底失重压力预测数值，建立了静胶凝过渡期环空压力实时调控技术，以确保地层流体受控，降低环空气窜风险。
4.2 顶替效率优化设计新方法
针对固井中的顶替效率问题，建立了复杂油气井顶替效率优化设计方法[31-32]，提出了真三维固井顶替数值模型，在国家天河一号超级计算中心开展了模拟研究，定量评价了不同工艺参数对顶替效率的影响机制。研究发现，直井段切力和黏度差、水平段密度差均与顶替效率正相关，并提出了黏度差、密度差、居中度、排量等关键参数及注水泥作业优化设计方法，应用结果表明，固井顶替效率可由65%提高至95%。
4.3 智能化设计新方法
针对智能化固井设计问题，基于计算模型，利用机器学习并融合标准规范及专家经验，开发了智能化固井设计软件，建立了云端化固井数据库，开发了智能化固井设计决策平台。该平台可根据现场井筒数据自动生成最优设计，同时利用机器学习对水泥浆配方大数据进行学习，基于水泥浆性能参数要求推荐水泥浆配方，可实现固井参数优化设计、过程模拟、实时监测、参数动态调整和质量分析的全过程控制。
深层超深层油气井固井技术在顺北油气田、元坝气田、川渝页岩气田等规模应用，固井质量合格率从55%提高至100%、优良率从23%提高至90%，页岩气井井口带压率从70.2%降低至7.6%，解决了高温高压气窜、酸性气田腐蚀、页岩气损伤等重大固井技术难题，支撑了深层超深层油气高效勘探开发。高端尾管悬挂器国内市场占有率70%。创造了水泥浆应用密度最高世界纪录（2.78 kg/L）和工具应用井深最深亚洲纪录（9 300 m）。
5. 发展建议
我国深层超深层油气资源丰富，是当前和未来油气勘探开发的重点和热点，2023年中国石化将在塔里木盆地钻探万米井，由于井更深、井下地质环境及工况更加复杂，固井技术将面临新的挑战。同时，深层页岩油气及干热岩等非常规资源的开发、国家对环保要求及数智化技术转型等也对固井技术提出新的要求或带来新的发展机遇。因此，在后续技术研究与发展过程中，应以“极端环境、绿色环保、智能跨越”为理念，在传统技术提升、环保及新能源技术拓展、智能化技术转型及基础理论研究推进等方面开展攻关研究，推动深层超深层油气井固井技术不断升级迭代。
1）传统优势技术方面，尽快攻克极端高温、极端低温、极端井型和极端应变等油气井固井技术，以满足万米井、超长水平段水平井、超大型压裂的固井要求。
2）环保及新能源技术方面，着力开发高应变、高酸性新材料及非硅酸盐基凝胶、地质聚合物等绿色环保材料，逐步替代传统水泥及助剂，以满足干热岩热储开采、高动载储气库、高酸性储氢库与CO2封存及国家对环保的要求。
3）智能化技术转型方面，建议研究大数据技术、智能化水泥浆和智能化固井监测仪器及装备，引领未来技术发展。
4）基础理论研究持续推进方面，以油气井全生命周期井筒完整性固井水泥密封为目标，进一步揭示水泥浆–水泥石力学化学演化规律，明确固井新材料新体系构建方向，提升水泥浆（石）性能，以满足更加复杂油气井长期密封的技术要求。
6. 结束语
2010年以来，中国石化立足自主创新，围绕“防窜、防腐、防损伤，保持井筒密封完整性”技术目标，经过持续攻关，在固井理论方法、水泥材料、工具装置和控制工艺方面取得了升级，解决了水泥环损伤、环空气窜、胶结质量差等重大技术难题，提高了深层超深层复杂油气井固井技术水平。后续技术研究过程中，应紧密结合国家能源发展战略，基于“极端环境、绿色环保、智能跨越”理念自主攻关，完善“三防一保”技术体系，攻关固井技术信息化与智能化，研发绿色、智能环保材料，持续提高复杂油气井的固井质量。
-
![]()
图 2 常温常压数字岩心离散元建模流程
Figure 2. Flowchart of digital core modeling at normal temperature and pressure by DEM
![]()
图 3 样品粒径及颗粒形状分析结果
Figure 3. Analysis results of particle size and particle shapes of the sample
![]()
图 4 球面描述符与SH阶数的关系曲线
Figure 4. Relationship between the spherical descriptor and SH degree
![]()
图 5 常温常压数字岩心离散元建模结果
Figure 5. Results of digital core modeling at normal temperature and pressure by DEM
![]()
图 7 重构岩心与真实岩心的两点相关和线性路径相关函数曲线对比结果
Figure 7. Comparison of two-point correlation and linear path correlation function curves between reconstructed and real cores
![]()
图 8 颗粒间微观力学参数标定模型及偏应力-应变曲线
Figure 8. Micromechanical parameter calibration model among particles and the deviational stress-strain curve
![]()
图 9 相同围压、不同轴压下的数字岩心孔隙几何拓扑结构
Figure 9. Pore geometry topology of digital core under same confining pressure and different axial pressures
![]()
图 10 相同轴向应力、不同围压下的数字岩心孔隙几何拓扑结构
Figure 10. Pore geometry topology of digital core under same axial stress and different confining pressures
![]()
图 11 不同应力组合作用下数字岩心孔隙度和渗透率的分布
Figure 11. Porosity and permeability distribution of digital core under different stress combinations
表 1 不同应力组合下的数字岩心构建结果
Table 1 Digital core construction results under different stress combinationss
水平应力 σz = 120 MPa σz = 150 MPa σz = 180 MPa σx = 50 MPa
σy = 70 MPa
σx = 70 MPa
σy = 90 MPa
σx = 90 MPa
σy = 110 MPa
下载: 导出CSV
-
[1] 贾承造,庞雄奇. 深层油气地质理论研究进展与主要发展方向[J]. 石油学报,2015,36(12):1457–1469. JIA Chengzao, PANG Xiongqi. Research processes and main development directions of deep hydrocarbon geological theories[J]. Acta Petrolei Sinica., 2015, 36(12): 1457–1469.
[2] 范家伟,袁野,李绍华,等. 塔里木盆地深层致密油藏地质工程一体化模拟技术[J]. 断块油气田,2022,29(2):194–198. FAN Jiawei, YUAN Ye, LI Shaohua, et al. Geology-engineering integrated simulation technology of deep tight oil reservoir in Tarim Basin[J]. Fault-Block Oil & Gas Field, 2022, 29(2): 194–198.
[3] 陈建勋. 深层高压碳酸盐岩气藏孔隙结构特征及衰竭开发规律[J]. 特种油气藏,2022,29(5):80–87. CHEN Jianxun. Pore structure characteristics and natural depletion law of deep high-pressure carbonate gas reservoirs[J]. Special Oil & Gas Reservoirs, 2022, 29(5): 80–87.
[4] 曾义金. 中国石化深层超深层油气井固井技术新进展与发展建议[J]. 石油钻探技术,2023,51(4):66–73. ZENG Yijin. Novel advancements and development suggestions of cementing technologies for deep and ultra-deep wells of Sinopec[J]. Petroleum Drilling Techniques, 2023, 51(4): 66–73.
[5] 王建云,韩涛,赵宽心,等. 塔深5井超深层钻井关键技术[J]. 石油钻探技术,2022,50(5):27–33. WANG Jianyun, HAN Tao, ZHAO Kuanxin, et al. Key drilling technologies for the ultra-deep Well Tashen 5[J]. Petroleum Drilling Techniques, 2022, 50(5): 27–33.
[6] 李虹,于海洋,杨海烽,等. 裂缝性非均质致密储层自适应应力敏感性研究[J]. 石油钻探技术,2022,50(3):99–105. LI Hong, YU Haiyang, YANG Haifeng, et al. Adaptive stress sensitivity study of fractured heterogeneous tight reservoir[J]. Petroleum Drilling Techniques, 2022, 50(3): 99–105.
[7] 孙强,孙志刚,张超. DLH油田低渗砂岩孔隙分形定量表征方法研究[J]. 西南石油大学学报(自然科学版),2023,45(1):105–116. SUN Qiang, SUN Zhigang, ZHANG Chao. A study on fractal quantitative characterization method of low permeability sandstone pore in DLH Oilfield[J]. Joumal of Southwest Petroleum University (Seience & Technology Edition), 2023, 45(1): 105–116.
[8] 姚军,黄朝琴,刘文政,等. 深层油气藏开发中的关键力学问题[J]. 中国科学:物理学 力学 天文学,2018,48(4):044701. YAO Jun, HUANG Zhaoqin, LIU Wenzheng, et al. Key mechanical problems in the development of deep oil and gas reservoirs[J]. SCIENTIA SINICA Physica, Mechanica & Astronomica, 2018, 48(4): 044701.
[9] 李荣强,高莹,杨永飞,等. 基于CT扫描的岩心压敏效应实验研究[J]. 石油钻探技术,2015,43(5):37–43. 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.
[10] 王晨晨,姚军,杨永飞,等. 基于CT扫描法构建数字岩心的分辨率选取研究[J]. 科学技术与工程,2013,13(4):1049–1052. WANG Chenchen, YAO Jun, YANG Yongfei, et al. Study on resolution selection for digital rock construction with CT scanning method[J]. Science Technology and Engineering, 2013, 13(4): 1049–1052.
[11] 刘洋,王春生,孙启冀,等. 低渗砂岩储层数字岩心构建及渗流模拟[J]. 断块油气田,2017,24(6):817–821. LIU Yang, WANG Chunsheng, SUN Qiji, et al. Digital core construction and seepage simulation of low permeability sandstone reservoir[J]. Fault-Block Oil & Gas Field, 2017, 24(6): 817–821.
[12] 姚军,宋文辉,李阳,等. 有机质孔隙对页岩气流动能力影响研究[J]. 中国科学:物理学 力学 天文学,2017,47(9):094702. doi: 10.1360/SSPMA2017-00040 YAO Jun, SONG Wenhui, LI Yang, et al. Study on the influence of organic pores on shale gas flow ability[J]. Scientia Sinaca: Physica, Mechanica & Astronomica, 2017, 47(9): 094702. doi: 10.1360/SSPMA2017-00040
[13] ØREN P E, BAKKE S. Process based reconstruction of sandstones and prediction of transport properties[J]. Transport in Porous Media, 2002, 46(2): 311–343.
[14] 姚军,赵秀才,衣艳静,等. 数字岩心技术现状及展望[J]. 油气地质与采收率,2005,12(6):52–54. YAO Jun, ZHAO Xiucai, YI Yanjing, et al. The current situation and prospect on digital core technology[J]. Petroleum Geology and Recovery Efficiency, 2005, 12(6): 52–54.
[15] YANG Yongfei, LIU Fugui, YAO Jun, et al. Multi-scale reconstruction of porous media from low-resolution core images using conditional generative adversarial networks[J]. Journal of Natural Gas Science and Engineering, 2022, 99: 104411. doi: 10.1016/j.jngse.2022.104411
[16] 杨永飞,刘夫贵,姚军,等. 基于生成对抗网络的页岩三维数字岩芯构建[J]. 西南石油大学学报(自然科学版),2021,43(5):73–83. YANG Yongfei, LIU Fugui, YAO Jun, et al. Reconstruction of 3D shale digital rock based on generative adversarial network[J]. Journal of Southwest Petroleum University (Science & Technology Edition), 2021, 43(5): 73–83.
[17] 陈林,黎棚武,张绍俊,等. 基于机器学习的岩芯渗透率及裂缝开度预测[J]. 西南石油大学学报(自然科学版),2023,45(4):155–163. CHEN Lin, LI Pengwu, ZHANG Shaojun, et al. Prediction method of core permeability and fracture aperture based on machine learning[J]. Journal of Southwest Petroleum University (Science & Technology Edition), 2023, 45(4): 155–163.
[18] SAENGER E H, LEBEDEV M, URIBE D, et al. Analysis of high-resolution X-ray computed tomography images of Bentheim sandstone under elevated confining pressures[J]. Geophysical Prospecting, 2016, 64(4): 848–859.
[19] ZHU Linqi, ZHANG Chong, ZHANG Chaomo, et al. Challenges and prospects of digital core-reconstruction researc[J]. Geofluids, 2019, 2019: 7814180.
[20] FAN Ming, MCCLURE J, HAN Yanhui, et al. Interaction between proppant compaction and single-/multiphase flows in a hydraulic fracture[J]. SPE Journal, 2018, 23(4): 1290–1303. doi: 10.2118/189985-PA
[21] CUNDALL P A, STRACK O D L. A discrete numerical model for granular assemblies[J]. Géotechnique, 1979, 29(1): 47–65. doi: 10.1680/geot.1979.29.1.47
[22] ZHENG Junxing, HRYCIW R D. Segmentation of contacting soil particles in images by modified watershed analysis[J]. Computers and Geotechnics, 2016, 73: 142–152. doi: 10.1016/j.compgeo.2015.11.025
[23] ZHAO Shiwei, ZHAO Jidong. A poly-superellipsoid-based approach on particle morphology for DEM modeling of granular media[J]. International Journal for Numerical and Analytical Methods in Geomechanics, 2019, 43(13): 2147–2169.
[24] ANDRADE J E, LIM K W, AVILA C F, et al. Granular element method for computational particle mechanics[J]. Computer Methods in Applied Mechanics and Engineering, 2012, 241(): 262–274.
[25] BOWMAN E T, SOGA K, DRUMMOND W. Particle shape characterisation using Fourier descriptor analysis[J]. Geotechnique, 2001, 51(6): 545–554.
[26] GARBOCZI E J. Three-dimensional mathematical analysis of particle shape using X-ray tomography and spherical harmonics: application to aggregates used in concrete[J]. Cement and concrete research, 2002, 32(10): 1621–1638.
[27] WEI Deheng, WANG Jianfeng, ZHAO Budi. A simple method for particle shape generation with spherical harmonics[J]. Powder Technology, 2018, 330: 284–291. doi: 10.1016/j.powtec.2018.02.006
[28] LI Zechuang, LIU Zhibin. Influence of particle shape on the macroscopic and mesolevel mechanical properties of slip zone soil based on 3D scanning and 3D DEM[J]. Advances in Materials Science and Engineering, 2021, 2021: 1–14.
[29] LAI Zhengshou, CHEN Qiushi. Characterization and discrete element simulation of grading and shape-dependent behavior of JSC-1A Martian regolith simulant[J]. Granular Matter, 2017, 19(4): 69.
[30] KHALEGHI K, TALMAN S, SHOKRI A R, et al. A coupled pore-scale modelling approach to capture macro-scale stress-dependent permeability of rocks[R]. URTEC 198264, 2019.
[31] KAWAMOTO R, ANDÒ E, VIGGIANI G, et al. Level set discrete element method for three-dimensional computations with triaxial case study[J]. Journal of the Mechanics and Physics of Solids, 2016, 91: 1–13.
[32] KLEIN E, BAUD P, REUSCHLÉ T et al. Mechanical behaviour and failure mode of bentheim sandstone under triaxial compression[J]. Physics and Chemistry of the Earth, Part A: Solid Earth and Geodesy, 2001, 26(1/2): 21–25.
[33] YANG Yongfei, ZHANG Wenjie, GAO Ying, et al. Influence of stress sensitivity on microscopic pore structure and fluid flow in porous media[J]. Journal of Natural Gas Science and Engineering, 2016, 36(part A): 20–31.
[34] LU Jun, YIN Guangzhi, DENG Bozhi, et al. Permeability characteristics of layered composite coal-rock under true triaxial stress conditions[J]. Journal of Natural Gas Science and Engineering, 2019, 66: 60–76. doi: 10.1016/j.jngse.2019.03.023
-
期刊类型引用(0)
其他类型引用(1)
计量
- 文章访问数: 225
- HTML全文浏览量: 47
- PDF下载量: 108
- 被引次数: 1
