CN116842789B - 一种考虑生产时间窗的加密井井壁稳定分析方法 - Google Patents
一种考虑生产时间窗的加密井井壁稳定分析方法 Download PDFInfo
- Publication number
- CN116842789B CN116842789B CN202310751169.XA CN202310751169A CN116842789B CN 116842789 B CN116842789 B CN 116842789B CN 202310751169 A CN202310751169 A CN 202310751169A CN 116842789 B CN116842789 B CN 116842789B
- Authority
- CN
- China
- Prior art keywords
- well
- pressure
- parameters
- stress
- encrypted
- 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.)
- Active
Links
- 238000004519 manufacturing process Methods 0.000 title claims abstract description 45
- 238000004458 analytical method Methods 0.000 title claims abstract description 30
- 239000011435 rock Substances 0.000 claims abstract description 67
- 238000005553 drilling Methods 0.000 claims abstract description 36
- 238000004088 simulation Methods 0.000 claims abstract description 34
- 230000008859 change Effects 0.000 claims abstract description 30
- 239000011148 porous material Substances 0.000 claims abstract description 28
- 238000012937 correction Methods 0.000 claims abstract description 12
- 238000002474 experimental method Methods 0.000 claims abstract description 7
- 238000013461 design Methods 0.000 claims abstract description 6
- 238000000034 method Methods 0.000 claims description 61
- 238000004364 calculation method Methods 0.000 claims description 24
- 230000015572 biosynthetic process Effects 0.000 claims description 18
- 230000008878 coupling Effects 0.000 claims description 3
- 238000010168 coupling process Methods 0.000 claims description 3
- 238000005859 coupling reaction Methods 0.000 claims description 3
- 238000005065 mining Methods 0.000 claims description 3
- 239000007787 solid Substances 0.000 claims description 3
- 230000007423 decrease Effects 0.000 claims description 2
- 238000009826 distribution Methods 0.000 description 17
- 230000008569 process Effects 0.000 description 14
- 238000010276 construction Methods 0.000 description 7
- 238000005516 engineering process Methods 0.000 description 7
- 238000004422 calculation algorithm Methods 0.000 description 4
- 238000011161 development Methods 0.000 description 4
- 230000000694 effects Effects 0.000 description 4
- 238000012545 processing Methods 0.000 description 4
- 238000012360 testing method Methods 0.000 description 4
- 238000007405 data analysis Methods 0.000 description 3
- 230000003993 interaction Effects 0.000 description 3
- 238000012986 modification Methods 0.000 description 3
- 230000004048 modification Effects 0.000 description 3
- 238000005457 optimization Methods 0.000 description 3
- 238000007781 pre-processing Methods 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 239000000243 solution Substances 0.000 description 3
- 230000009466 transformation Effects 0.000 description 3
- 238000013528 artificial neural network Methods 0.000 description 2
- 238000004140 cleaning Methods 0.000 description 2
- 230000007613 environmental effect Effects 0.000 description 2
- 239000011159 matrix material Substances 0.000 description 2
- 230000035699 permeability Effects 0.000 description 2
- 238000003908 quality control method Methods 0.000 description 2
- 230000000452 restraining effect Effects 0.000 description 2
- 238000010008 shearing Methods 0.000 description 2
- 230000006641 stabilisation Effects 0.000 description 2
- 238000011105 stabilization Methods 0.000 description 2
- 230000003068 static effect Effects 0.000 description 2
- 238000007619 statistical method Methods 0.000 description 2
- 230000001360 synchronised effect Effects 0.000 description 2
- 230000001133 acceleration Effects 0.000 description 1
- 230000003044 adaptive effect Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000009933 burial Methods 0.000 description 1
- 230000001808 coupling effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000011439 discrete element method Methods 0.000 description 1
- 230000009977 dual effect Effects 0.000 description 1
- 239000012530 fluid Substances 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 238000011065 in-situ storage Methods 0.000 description 1
- 238000002347 injection Methods 0.000 description 1
- 239000007924 injection Substances 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 239000003129 oil well Substances 0.000 description 1
- 239000002245 particle Substances 0.000 description 1
- 230000002093 peripheral effect Effects 0.000 description 1
- 230000000704 physical effect Effects 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000007670 refining Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 239000002002 slurry Substances 0.000 description 1
- 230000000087 stabilizing effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B49/00—Testing the nature of borehole walls; Formation testing; Methods or apparatus for obtaining samples of soil or well fluids, specially adapted to earth drilling or wells
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/08—Fluids
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Fluid Mechanics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Geometry (AREA)
- Evolutionary Computation (AREA)
- Mining & Mineral Resources (AREA)
- General Engineering & Computer Science (AREA)
- Computer Hardware Design (AREA)
- Geology (AREA)
- Computing Systems (AREA)
- Mathematical Physics (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Environmental & Geological Engineering (AREA)
- Mathematical Analysis (AREA)
- Algebra (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geochemistry & Mineralogy (AREA)
- Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)
Abstract
本发明公开了一种考虑生产时间窗的加密井井壁稳定分析方法,其技术方案是:包括将初步井轨迹加载至储层三维地质力学模型中;其次,开展父井的水力压裂数值模拟,得到水力裂缝的实际扩展形态;之后,模拟不同生产时间条件的地层压力参数,并通过历史拟合对地层压力参数进行反演校正;对目的储层的岩心进行采集,实验得到不同孔隙压力条件下对应的力学参数变化规律,构建不同孔隙压力条件下的岩石力学参数变化方程;最后,导入加密井井轨迹到动态地质力学场模型中,根据优化后的安全钻井泥浆窗口指导加密井钻井方案设计。本发明能够提高加密井钻井的井壁稳定性,从而减少钻井事故的发生,提高加密井钻井效率,保证后续加密井压裂改造顺利完成。
Description
技术领域
本发明涉及油气钻井工程技术领域,特别涉及一种考虑生产时间窗的加密井井壁稳定分析方法。
背景技术
随着页岩气等非常规油气的开发持续进行,体积压裂井逐渐暴露出缝控储层有限,产气量下降速度过快等问题。考虑到开发初期井网较大(400m左右),井间剩余储层较高,因此,对页岩等非常规油气储层进行加密井体积压裂可以有效增加整个储层的井间改造程度,实现井间剩余储量的突破。加密井进行钻井设计时需要考虑老井生产的影响。主要原因在于,老井的衰竭生产将导致井周附近出现应力衰竭和反转,同时生产的强度和生产时间制度会造成应力场也存在较大差异,进而会影响后续的加密井井壁稳定性。目前在进行加密井钻井设计时,主要考虑增加加密井后期压裂的预期改造和提产效果,较少从井壁稳定性角度对加密井钻井进行优化设计。
井壁稳定研究是目前油气工业界中重点关注的热点方向。目前多数的井壁稳定主要通过解析解进行破裂和坍塌压力的计算,不少学者也采用数值模拟方法进行井壁稳定分析,常见的数值模拟方法主要包括有限元法和离散元法等。由于不少非常规储层具有埋深大、断裂裂缝发育、地温高、地应力强以及地层压力大等问题,容易导致在钻井过程中常发生严重的钻井事故。对于密集式布缝的井平台而言,由于油气长期开采产生的父井地层压力持续下降将对地层的应力产生非常重要的影响,从而会改变加密井子井的井壁稳定性。但是,目前研究加密井井壁稳定问题相关成果较少。搞清加密井动态地质力学变化规律,构建基于动态地质力学模型的加密井井壁稳定模型,将对于减少加密井钻井时的井壁失稳事故的发生极为重要。
发明内容
本发明的目的就是针对现有技术存在的上述缺陷,提供一种考虑生产时间窗的加密井井壁稳定分析方法。
本发明提到的一种考虑生产时间窗的加密井井壁稳定分析方法,其技术方案是:包括以下步骤:
(1)建立储层三维地质力学模型,并将初步井轨迹加载至所述储层三维地质力学模型中;
(2)其次,利用地质力学模型参数,开展父井的水力压裂数值模拟,得到水力裂缝的实际扩展形态;之后,通过油藏数值模拟软件模拟不同生产时间条件的地层压力参数,并通过历史拟合对地层压力参数进行反演校正;
(3)对目的储层的岩心进行采集,通过卸载条件下岩石力学实验得到不同孔隙压力条件下对应的力学参数变化规律,构建不同孔隙压力条件下的岩石力学参数变化方程。以地层压力参数为输入参数,再次导入地质力学模型中进行动态岩石力学参数和地应力参数计算,得到不同时间窗条件的动态地质力学场;
(4)导入加密井井轨迹到动态地质力学场模型中,并同步采集不同模拟动态地应力场井轨迹条件下的应力数据和更新岩石力学参数导入热流固耦合模型中进行破裂压力和坍塌压力的计算,结合更新地层孔隙压力数据形成加密井钻井的泥浆密度窗口;
(5)根据优化后的安全钻井泥浆窗口指导加密井钻井方案设计。
优选的,地质力学模型的建模方法需要包括构造建模、属性建模、岩石力学建模和地应力建模。
优选的,利用地质力学模型的参数,开展父井的水力压裂数值模拟,得到水力裂缝的实际扩展形态;同时,基于父井压裂得到的水力压裂裂缝形态,应当进一步进行数值模拟,分析不同开采时间条件下地层压力,并实现历史拟合。
优选的,由于地层衰竭产生的孔隙压力变化,并由于孔隙压力变化导致的岩石力学参数变化需要通过室内真三轴卸载实验完成,同时需要根据室内实验回归数据得到孔隙压力场变化条件下的岩石力学参数下降回归方程。
优选的,通过历史拟合和数值模拟得到的地层压力数据应当进一步导入到地应力场分析软件中来计算不同时刻条件下的动态地应力变化规律。
优选的,在动态的地质力学模型中提取加密井井身轨迹附近的相关岩石力学参数和地应力参数,相关岩石力学参数应当至少包括杨氏模量、泊松比、抗压强度、内摩擦角、内聚力等,相关地应力参数则主要包括水平最大主应力、水平最小主应力和垂向应力的所述地应力;以及相关井身结构参数、如井周角和井斜角、以及岩石孔隙度、以及Biot系数。
优选的,需要将不同时间窗条件下的参数导入至井壁稳定力学模型中对地层的孔隙压力P p (t)、破裂压力P f (t)和坍塌压力P c (t)进行计算,由此得到不同时间窗条件下的加密井钻井泥浆密度窗口。
优选的,根据安全钻井泥浆窗口压力值可能最小值对应数据应当为加密井井周不同时刻条件下井周孔隙压力值或者坍塌压力值,该值需要通过不同时刻的地层压力值或者坍塌压力计算公式获取。
优选的,根据安全钻井泥浆窗口压力值最大值对应数据应当为加密井井周不同时刻条件下井周破裂压力值。
与现有技术相比,本发明的有益效果具体如下:
(1)本发明提出了地质力学和油藏数值模拟来分析加密井井身规律的地质力学参数变化,通过地质力学和油藏数值模拟器实现力学参数和渗流参数的同步交互,构建的模型考虑了父井生产导致的加密井井周动态地质力学特征变化;
(2)本发明考虑了父井生产导致的岩石力学场和应力场的同步变化,并利用提取的岩石力学和应力场计算结果作为钻井泥浆密度窗口中破裂压力和坍塌压力数值的计算输入参数,同时坍塌压力和破裂压力中相关的关键岩石力学参数,如杨氏模量、内聚力和内摩擦角等实验结果同样需要室内试验回归校正后输入到破裂压力和坍塌压力计算模型中,从而通过实验和数值模拟双重技术手段实现了动态地质力学参数的不同更新,较原有的纯数值模拟方法具有更高的精度;
(3)本发明旨在创造一种考虑生产时间窗影响下的加密井井壁稳定计算模型,而传统考虑常规静态地质力学特征的井壁研究模型相对较多,并采用静态原始状态下的力学和应力参数作为输入,无法有效指导加密井的钻井井壁失稳分析。
附图说明
图1 为本发明的整体流程示意图;
图2 为本发明构建的不同条件下三维地应力模型数据图;
图3为本发明根据构建模型提取的不同时刻的沿井孔隙压力场数据图;
图4为本发明根据构建模型提取的不同时刻的生产井地应力场云图;
图5为本发明不同生产时间窗条件下计算的破裂压力分析数据图;
图6为本发明不同生产时间窗条件下计算的坍塌压力分析数据图;
图7为本发明不同生产时间窗条件下计算的破裂压力云图。
具体实施方式
以下结合附图对本发明的优选实施例进行说明,应当理解,此处所描述的优选实施例仅用于说明和解释本发明,并不用于限定本发明。
实施例1,本发明提到的一种考虑生产时间窗的加密井井壁稳定分析方法,
基于体元构造框架建模技术建立准确的区块构造模型,通过测、录井和地震数据,并结合小层厚度标定,开展三维精细地质建模;
基于建立的油藏构造模型,运用地震资料、钻井岩心等多种手段和数据源进行综合分析,以确定不同岩性的空间分布和岩石物理参数的变化规律,并结合相控技术建立油藏地质属性模型;
利用神经网络结合地质约束条件与多点统计分析方法,通过学习地质属性的空间分布和变异规律与克里金插值模拟岩性的空间变异性,生成符合实际情况的岩性分布;
收集测井、岩心、室内试验等数据,以岩性分析为约束结果,通过数值模拟方法得到关键油藏属性的参数分布特征,主要包括孔隙度、渗透率、含气饱和度等;
收集测井、岩心、地震等数据,包括储层横纵波、密度测井数据与地震横波阻抗、拉梅阻抗、剪切阻抗等。对收集到的数据进行数据清洗、异常值处理、缺失值填充等预处理,使数据符合建模要求。对预处理后的数据进行分析,建立储层属性与岩石力学属性之间的关系,如密度、杨氏模量、泊松比等。根据数据分析结果,建立三维岩石力学模型,包括储层模型、岩石力学属性模型、地震模型等。建立模型的过程需要考虑储层的空间结构、岩石力学属性的分布规律、地震波传播的特性等因素。利用建立的模型,对储层属性进行三维反演,得到储层密度、杨氏模量、泊松比等岩石力学属性。最后,利用测井解释的岩石力学结果和室内岩石力学实验成果进行环境约束和质量控制,保证反演结果的准确性和可靠性;
结合地震解释成果与有限元方法,将地下结构与地质构造导入有限元模型,利用网格交互算法将地层属性同步导入有限元计算平台,通过粒子群优化算法进行有限元数值计算,迭代求解得到三维地应力模型,同时通过收集地应力实验结果和现场小压分析等实验结果对应力大小进行约束,通过古地磁和波向异性结果,以及成像测井解释结果对应力方向进行约束;
将初步井轨迹加载至所述储层三维地质力学模型中,获取初步井轨迹与储层三维地质力学解释成果;所述在加密井井身轨迹附近提取井眼数据主要包括:井眼半径、井斜角、方位角等,在所述加密井井身附近提取的应力参数主要包括水平最大主应力、水平最小主应力和垂向应力的所述地应力;在所述加密井井身附近提取的岩石力学参数主要包括岩石抗压强度、岩石抗剪强度、内聚力、内摩擦角等参数。
通过构建的三维储层属性模型、岩石力学模型和地应力模型开展父井的水力压裂数值模拟,得到水力裂缝的实际扩展形态。结合多种数据和技术手段,包括地质学、岩石力学等方面的知识,并在模拟过程中对数据进行充分的处理和分析,以保证水力压裂模型的可靠性和准确性。同时,利用井间微震数据反演对水力压裂数值模拟形态进行校正和更新;
通过油藏数值模拟软件模拟不同生产时间条件的地层压力参数等属性数据,通过历史拟合对地层压力等参数进行反演校正,得到不同开采时间窗条件下的孔隙压力场结果;
对目的储层的岩心进行采集,通过卸载条件下岩石力学实验得到不同孔隙压力条件下对应的力学参数变化规律,构建不同孔隙压力条件下的岩石力学参数变化方程。以地层压力参数为输入参数,再次导入地质力学模型中进行动态岩石力学参数和地应力参数计算,得到不同时间窗条件的动态地质力学场;
导入加密井井轨迹到动态地质力学场模型中,并同步采集不同模拟动态地应力场沿井轨迹条件下的井周更新岩石力学参数和更新地应力数据,并导入力学计算模型中进行破裂压力计算,其中破裂压力计算模型主要采用Stephen破裂力学模型进行破裂压力计算;
导入加密井井轨迹到动态地质力学场模型中,并同步采集不同模拟动态地应力场沿井轨迹条件下的井周更新岩石力学参数和更新地应力数据,并导入力学计算模型中进行坍塌压力计算,其中坍塌压力计算模型主要采用Mogi-Column坍塌力学模型进行坍塌压力计算;
本发明提到的一种考虑生产时间窗的加密井井壁稳定分析方法,得到所述安全钻井泥浆窗口包括:最小安全钻井泥浆窗口压力值通过不同时刻条件下的孔隙压力值和坍塌压力值进行确定;根据所述最大安全钻井泥浆窗口压力值通过不同时刻条件下的破裂压力值进行确定。
实施例2,本发明提到的一种考虑生产时间窗的加密井井壁稳定分析方法,包括以下步骤:
根据目标储层钻井资料、断层解释结果和临层对比数据,基于体元构造框架建模技术建立准确的区块构造模型与小层发育构造框架。通过测、录井和地震数据,并结合导眼井分层数据进行各小层的厚度标定,同时利用井间数据修正与校核,实现包含小层的区块三维精细地质建模;
基于建立的油藏构造模型,运用地震资料、钻井岩心等多种手段和数据源进行综合分析,以确定不同岩性的空间分布和岩石物理参数的变化规律,通过数值模拟建立油藏属性模型。由于地震资料、钻井岩心等多种手段和数据源获取的地质信息数据格式不同,需要对数据进行校正、插值等预处理,确保数据的一致性和准确性,以便更好地进行后续的分析和建模。再进行综合分析,确定不同岩性的空间分布和岩石物理参数的变化规律,并结合相控技术更加准确地确定岩石属性和构造特征,建立油藏地质属性模型。
通过地震资料、钻井岩心、地质测量等多种数据来建立准确的地质属性模型。对数据进行分析和处理后,可以得到地质属性的空间分布和变异规律,从而对随机模拟结果进行约束。在使用地震资料进行井间岩性随机模拟时,可以根据地震速度和密度等数据,推导出地质属性的分布规律将这些规律应用到井间岩性随机模拟中,以获得更符合实际情况的结果。利用克里金插值方法来描述岩性的空间变异性,并对随机模拟结果进行控制。利用神经网络结合地质约束条件与多点统计分析方法,通过学习地质属性的空间分布和变异规律与克里金插值模拟岩性的空间变异性,生成符合实际情况的岩性分布。
收集测井、岩心、室内试验等数据,以岩性分析为约束结果,通过数值模拟方法得到关键油藏属性的参数分布特征,主要包括孔隙度、渗透率、含气饱和度等;
收集测井、岩心、地震等数据,包括储层横纵波、密度测井数据与地震横波阻抗、拉梅阻抗、剪切阻抗等。对收集到的数据进行数据清洗、异常值处理、缺失值填充等预处理,使数据符合建模要求。对预处理后的数据进行分析,建立储层属性与物理属性之间的关系,如密度、杨氏模量、泊松比等。通过数据分析,可以了解储层的物理特征,为建立模型提供依据。根据数据分析结果,建立三维岩石力学模型,包括储层模型、岩石力学属性模型、地震模型等。建立模型的过程需要考虑储层的空间结构、岩石力学属性的分布规律、地震波传播的特性等因素。利用建立的模型,对储层属性进行三维反演,得到储层密度、杨氏模量、泊松比等岩石力学属性。最后,利用测井解释的岩石力学结果和室内岩石力学实验成果进行环境约束和质量控制,保证反演结果的准确性和可靠性。利用测井数据计算杨氏模量的公式为:
(1)
泊松比的计算公式为:
(2)
密度的计算公式为:
(3)
垂向应力计算公式:
(4)
式中,E为杨氏模量,GPa;t s、t p为横波、纵波时间差,无因次;υ为泊松比,无因次;ρ为地层密度,g/m3;ρ m 为岩石基质密度;ρ f 为钻井液密度,g/m3;ΔT为测井曲线测量的时间(或深度)差,无因次;ΔT m为岩石基质的时间(或深度)差,无因次;ρ i 为密度,g/m3;ΔD i 为深度采样间隔,m。
开展地应力场的构建,通过单井地应力纵向剖面数据,以及古地磁、波速各向异性、声发射和差应变等相关实验数据对构建的地应力场进行空间反演,得到三维应力场分布特征;在三维岩石力学参数场的基础上,结合地应力分析测试资料,运用有限元方法模拟三维地应力场分布特征,同时运用断层局部应力修正技术对断层附近区域的应力值进行局部修正,同时对最大水平主应力方向进行偏转分析,一般遵循最大主应力呈垂直于断层走向偏转,最小主应力呈平行于断层走向偏转的相关标准。
(5)
式中,σ H ,σ h 分别为最大和最小水平主应力,MPa;υ为岩石泊松比,无因次,σ v为垂向应力,MPa;E为岩石杨氏模量,MPa;β 1 ,β 2 分别为最大和最小水平主应力方向的构造应力系数,无因次;α为比奥特系数,无因次;P p 为地层孔隙压力,MPa。
水平地应力的范围可以通过孔隙压力与上覆岩层压力同时进行约束得到,即形成应力多边形。通常情况下,在地下任何位置的水平应力大小都小于应力多边形的边界值,这是由于地下断层与天然裂缝的存在,其产生的摩擦力对水平应力的大小起到了一定的控制作用。
根据应力多边形理论,走滑断层的最大主应力上限为:
(6)
最小水平主应力的下限为:
(7)
式中:S v 为垂向应力,MPa;P p 为孔隙压力,MPa;S H 为最大水平有效应力,MPa;S h 为最小水平有效应力,MPa;q f 为内摩擦角计算的参数。
结合地震驱动技术与有限元方法,将地下结构与地质构造导入有限元模型,利用网格交互算法将地层属性导入地应力数值模拟平台,通过有限元模拟计算复杂地应力场,迭代求解得到三维地应力模型,同时也可以进行参数反演和模型优化,提高模拟结果的准确性和可靠性;
进行有限元分析时,需要进行网格调整与优化。有限元计算时网格规模的大小对计算结果和计算效率都有很大的影响,采用基于误差控制的自适应网格技术,通过评估每个网格单元的误差大小,对不同区域的网格进行细化或粗化,有效地控制网格规模。并采用基于三角剖分的网格生成算法,通过将模型划分为多个三角形来生成网格,有效地提高模拟计算的效率和准确性;
将初步井轨迹加载至上述建立的储层三维地质力学模型中,获取初步井轨迹与储层三维地质力学模型关键井眼参数、岩石力学参数和地应力参数;
通过构建的三维储层属性模型、岩石力学模型和地应力模型开展父井的水力压裂数值模拟,得到水力裂缝的实际扩展形态。结合多种数据和技术手段,包括地质学、地球物理学、岩石力学等方面的知识,并在模拟过程中对数据进行充分的处理和分析,以保证水力压裂模型的可靠性和准确性。利用井间微震数据反演对水力压裂数值模拟形态进行校正和更新;
通过油藏数值模拟软件模拟不同生产时间条件的地层压力参数,并通过历史拟合对地层压力参数进行反演校正;以地质属性与生产数据作为输入参数,
以油井的生产压力和注入压力为边界条件开展油藏数值模拟,得到不同时刻的储层压力场。将模拟计算得到的地层压力曲线与实际生产数据中的地层压力曲线进行比较,进行历史拟合。根据历史拟合得到的校正系数,对模拟计算得到的地层压力曲线进行反演校正,得到更加准确的地层压力参数P p (t)。压力衰竭地层井壁孔隙压力公式
(8)
式中,p p (t)为不同生产时间窗下的地层孔隙压力,g/m3;
对目的储层的岩心进行采集,通过卸载条件下岩石力学实验得到不同孔隙压力条件下对应的力学参数变化规律,构建不同孔隙压力条件下的岩石力学参数变化方程,在此基础上推导致密储层开发过程中热流固耦合作用本构方程,建立关键耦合参数。
以历史拟合校正后的地层压力参数为输入参数,再次导入地质力学模型中进行动态岩石力学参数和地应力参数计算,利用有限元平台进行动态地应力场建模求解计算,分析不同时刻条件下的动态地应力变化规律,形成不同开采时间窗条件的动态地质力学场;
定义5个坐标体系,分别是大地坐标系,原地应力坐标系,井轴坐标系,井轴柱坐标系,以及弱面坐标系表示。根据坐标转转换,得到井周应力场方程为:
(9)
(10)
式中,、/>、/>分别为直角坐标下X、Y、Z轴方向上地应力分量,MPa;σ H为最大水平地应力,MPa;σ h为最小水平地应力,MPa;σ v为上覆岩层地应力,MPa;P m为井筒压力,MPa;α为井斜角,(°);γ为方位角,(°);/>为圆周角,(°);σ r为井周地层的径向应力分量,MPa;σ θ为井周地层的环向应力分量,MPa;σ z为井周地层的轴向应力分量,MPa;/>、/>、为井周地层的剪切应力分量,MPa。
根据压力衰竭更新地应力和岩石力学参数,加密井井壁破裂压力公式:
(11)
式中,p f 为不同生产时间窗下的破裂压力,g/m3;St为不同生产时间窗下的岩石抗拉强度,MPa。
对于井壁坍塌,井壁处三个应力状态为:
(12)
式中,σ 1为最大主应力,MPa;σ 2为中间主应力,MPa;σ 3为最小主应力,MPa;;δ为渗流效应系数,无因次;ϕ为孔隙度,%。
根据MG-C准则计算坍塌压力:
(13)
式中,p c 为不同生产时间窗下的坍塌压力,g/m3;
通过牛顿迭代法求解,得到安全泥浆密度窗口上、下限:
(14)
(15)
式中,g为重力加速度,9.8 m·s-2;TVD为垂深,m。
地层压力衰竭后必然导致坍塌压力和破裂压力变化,导致安全泥浆密度窗口发生改变。将不同时间窗条件下的参数导入地层的破裂压力和坍塌压力进行计算,由于地应力随时间动态变化,直接体现在坍塌压力和破裂压力随地应力动态变化而动态变化,由此得到不同生产时间窗条件下的泥浆密度窗口。
结合更新地层孔隙压力和坍塌压力中最小压力数据形成加密井钻井的泥浆密度最小密度窗口,根据更新地层破裂压力数据形成加密井钻井的泥浆密度最大密度窗口。根据所述优化后的安全钻井泥浆窗口指导加密井钻井方案设计。
以上所述,仅是本发明的部分较佳实施例,任何熟悉本领域的技术人员均可能利用上述阐述的技术方案加以修改或将其修改为等同的技术方案。因此,依据本发明的技术方案所进行的相应简单修改或等同变换,尽属于本发明要求保护的范围。
Claims (9)
1.一种考虑生产时间窗的加密井井壁稳定分析方法,其特征是:包括以下步骤:
(1)建立储层三维地质力学模型,并将初步井轨迹加载至所述储层三维地质力学模型中;
(2)其次,利用地质力学模型参数,开展父井的水力压裂数值模拟,得到水力裂缝的实际扩展形态;之后,通过油藏数值模拟软件模拟不同生产时间条件的地层压力参数,并通过历史拟合对地层压力参数进行反演校正;
(3)对目的储层的岩心进行采集,通过卸载条件下岩石力学实验得到不同孔隙压力条件下对应的力学参数变化规律,构建不同孔隙压力条件下的岩石力学参数变化方程;以地层压力参数为输入参数,再次导入地质力学模型中进行动态岩石力学参数和地应力参数计算,得到不同时间窗条件的动态地质力学场;
(4)导入加密井井轨迹到动态地质力学场模型中,并同步采集不同模拟动态地应力场井轨迹条件下的应力数据和更新岩石力学参数导入热流固耦合模型中进行破裂压力和坍塌压力的计算,结合更新地层孔隙压力数据形成加密井钻井的泥浆密度窗口;
(5)根据优化后的安全钻井泥浆窗口指导加密井钻井方案设计。
2.根据权利要求1所述的考虑生产时间窗的加密井井壁稳定分析方法,其特征是:地质力学模型的建模方法需要包括构造建模、属性建模、岩石力学建模和地应力建模。
3.根据权利要求2所述的考虑生产时间窗的加密井井壁稳定分析方法,其特征是:利用地质力学模型的参数,开展父井的水力压裂数值模拟,得到水力裂缝的实际扩展形态;同时,基于父井压裂得到的水力压裂裂缝形态,应当进一步进行数值模拟,分析不同开采时间条件下地层压力,并实现历史拟合。
4.根据权利要求3所述的考虑生产时间窗的加密井井壁稳定分析方法,其特征是:由于地层衰竭产生的孔隙压力变化,并由于孔隙压力变化导致的岩石力学参数变化需要通过室内真三轴卸载实验完成,同时需要根据室内实验回归数据得到孔隙压力场变化条件下的岩石力学参数下降回归方程。
5.根据权利要求4所述的考虑生产时间窗的加密井井壁稳定分析方法,其特征是:通过历史拟合和数值模拟得到的地层压力数据应当进一步导入到地应力场分析软件中来计算不同时刻条件下的动态地应力变化规律。
6.根据权利要求5所述的考虑生产时间窗的加密井井壁稳定分析方法,其特征是:在动态的地质力学模型中提取加密井井身轨迹附近的相关岩石力学参数和地应力参数,相关岩石力学参数应当至少包括杨氏模量、泊松比、抗压强度、内摩擦角、内聚力等,相关地应力参数则主要包括水平最大主应力、水平最小主应力和垂向应力的所述地应力;以及相关井身结构参数、如井周角和井斜角、以及岩石孔隙度、以及Biot系数。
7.根据权利要求6所述的考虑生产时间窗的加密井井壁稳定分析方法,其特征是:需要将不同时间窗条件下的参数导入至井壁稳定力学模型中对地层的孔隙压力P p (t)、破裂压力P f (t)和坍塌压力P c (t)进行计算,由此得到不同时间窗条件下的加密井钻井泥浆密度窗口。
8.根据权利要求7所述的考虑生产时间窗的加密井井壁稳定分析方法,其特征是:根据安全钻井泥浆窗口压力值可能最小值对应数据应当为加密井井周不同时刻条件下井周孔隙压力值或者坍塌压力值,该井周孔隙压力值或者坍塌压力值需要通过不同时刻的地层压力值或者坍塌压力计算公式获取。
9.根据权利要求8所述的考虑生产时间窗的加密井井壁稳定分析方法,其特征是:根据安全钻井泥浆窗口压力值最大值对应数据应当为加密井井周不同时刻条件下井周破裂压力值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310751169.XA CN116842789B (zh) | 2023-06-25 | 2023-06-25 | 一种考虑生产时间窗的加密井井壁稳定分析方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310751169.XA CN116842789B (zh) | 2023-06-25 | 2023-06-25 | 一种考虑生产时间窗的加密井井壁稳定分析方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN116842789A CN116842789A (zh) | 2023-10-03 |
CN116842789B true CN116842789B (zh) | 2024-04-05 |
Family
ID=88164580
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202310751169.XA Active CN116842789B (zh) | 2023-06-25 | 2023-06-25 | 一种考虑生产时间窗的加密井井壁稳定分析方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN116842789B (zh) |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110454127A (zh) * | 2019-07-04 | 2019-11-15 | 成都理工大学 | 一种非常规油气藏加密井有利加密时间窗口确定方法 |
CN110704888A (zh) * | 2019-07-04 | 2020-01-17 | 成都理工大学 | 一种非常规油气藏加密井体积压裂施工参数优化设计方法 |
CN110705000A (zh) * | 2019-07-04 | 2020-01-17 | 成都理工大学 | 非常规储层加密井压裂动态微地震事件屏障区确定方法 |
CN113821953A (zh) * | 2021-09-22 | 2021-12-21 | 中国石油大学(华东) | 基于动态地应力和剩余油特征的暂堵重复压裂设计方法 |
CN113901681A (zh) * | 2021-09-22 | 2022-01-07 | 中国石油大学(华东) | 一种全寿命周期页岩气储层双甜点三维可压性评估方法 |
CN115906419A (zh) * | 2022-10-26 | 2023-04-04 | 安东石油技术(集团)有限公司 | 超深储层钻井井壁稳定的方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US11009620B2 (en) * | 2019-07-04 | 2021-05-18 | Chengdu University Of Technology | Method for determining favorable time window of infill well in unconventional oil and gas reservoir |
-
2023
- 2023-06-25 CN CN202310751169.XA patent/CN116842789B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110454127A (zh) * | 2019-07-04 | 2019-11-15 | 成都理工大学 | 一种非常规油气藏加密井有利加密时间窗口确定方法 |
CN110704888A (zh) * | 2019-07-04 | 2020-01-17 | 成都理工大学 | 一种非常规油气藏加密井体积压裂施工参数优化设计方法 |
CN110705000A (zh) * | 2019-07-04 | 2020-01-17 | 成都理工大学 | 非常规储层加密井压裂动态微地震事件屏障区确定方法 |
CN113821953A (zh) * | 2021-09-22 | 2021-12-21 | 中国石油大学(华东) | 基于动态地应力和剩余油特征的暂堵重复压裂设计方法 |
CN113901681A (zh) * | 2021-09-22 | 2022-01-07 | 中国石油大学(华东) | 一种全寿命周期页岩气储层双甜点三维可压性评估方法 |
CN115906419A (zh) * | 2022-10-26 | 2023-04-04 | 安东石油技术(集团)有限公司 | 超深储层钻井井壁稳定的方法 |
Non-Patent Citations (4)
Title |
---|
Experimental study on hydraulic fracture propagation behavior on dolomite rocks with open-hole completion;shixian等;《JOURNAL OF PETROLEUM SCIENCE AND ENGINEERING》;20220930;全文 * |
四维动态地应力建模方法及其在寿阳区块煤层气开发中的应用;张滨海;陈峥嵘;艾传志;李莹莹;彭成勇;朱海燕;唐煊赫;;中国海上油气;20180822(第04期);全文 * |
基于扩展有限元的水平井改进拉链式压裂数值模拟;冯其红等;《中国石油大学学报( 自然科学版)》;20190430;第43卷(第2期);全文 * |
声波全波列测井资料的处理方法及应用;江明, 李兆阳;石油仪器;20051030(第05期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN116842789A (zh) | 2023-10-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US10101498B2 (en) | Well survivability in multidimensional geomechanical space | |
EP3036398B1 (en) | Formation stability modeling | |
CN113821953B (zh) | 基于动态地应力和剩余油特征的暂堵重复压裂设计方法 | |
EP2497900B1 (en) | Modeling hydraulic fractures | |
US10810331B2 (en) | System for predicting induced seismicity potential resulting from injection of fluids in naturally fractured reservoirs | |
WO2020080973A1 (ru) | Способ и система комбинированного сопровождения процесса бурения скважины | |
CN106014399B (zh) | 一种非均质地层高精度三维地应力模型建立方法 | |
Liu et al. | Asymmetric propagation mechanism of hydraulic fracture networks in continental reservoirs | |
CN106126866A (zh) | 基于地质力学模型的油气井防砂方法及装置 | |
CN105386756A (zh) | 一种应用应变量计算脆性地层孔隙度的方法 | |
CN117313472A (zh) | 一种缝洞体碳酸盐岩储层的重复压裂参数优化设计方法 | |
Bui et al. | A Coupled Geomechanics-Reservoir Simulation Workflow to Estimate the Optimal Well-Spacing in the Wolfcamp Formation in Lea County | |
CN115324559B (zh) | 压裂诱发油气套管变形的多因素综合预测和预防方法 | |
CN110705168A (zh) | 构造应力场的模拟方法 | |
CN112412434B (zh) | 一种改进的疏松砂岩地应力计算方法 | |
CN116842789B (zh) | 一种考虑生产时间窗的加密井井壁稳定分析方法 | |
CN113356843B (zh) | 针对地层的井壁稳定性分析方法、装置、介质及设备 | |
Wang et al. | Application of anisotropic wellbore stability model and unconventional fracture model for lateral landing and wellbore trajectory optimization: a case study of shale gas in Jingmen Area, China | |
CN115324556A (zh) | 一种压裂诱发油气套管变形风险级别的综合预测方法 | |
CN109339776B (zh) | 一种各向异性地层地应力方位测量方法 | |
CN113187463B (zh) | 一种基于地层超压单成因贡献率的随钻孔隙压力预测方法 | |
Zhang et al. | Hydromechanical Modeling of Nonplanar Three-Dimensional Fracture Propagation Using an Iteratively Coupled Approach | |
CN118187842A (zh) | 一种考虑加密井欠平衡钻井井壁坍塌压力计算方法 | |
Ibrahim et al. | A Novel Approach Linking Structural Restoration and Multiphase Fluid-Flow Modeling Using 3D Geomechanics in a North Kuwait Giant Carbonate Reservoir | |
Liu | Geomechanical Reservoir Characteriazation and Applications |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |