CN110543745B - 一种热采井水泥环封固完整性的评价方法 - Google Patents
一种热采井水泥环封固完整性的评价方法 Download PDFInfo
- Publication number
- CN110543745B CN110543745B CN201910930508.4A CN201910930508A CN110543745B CN 110543745 B CN110543745 B CN 110543745B CN 201910930508 A CN201910930508 A CN 201910930508A CN 110543745 B CN110543745 B CN 110543745B
- Authority
- CN
- China
- Prior art keywords
- cement
- temperature
- casing
- stratum
- cement sheath
- 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
- 239000004568 cement Substances 0.000 title claims abstract description 154
- 238000004519 manufacturing process Methods 0.000 title claims abstract description 28
- 238000000034 method Methods 0.000 title claims abstract description 24
- 238000007789 sealing Methods 0.000 title claims abstract description 21
- 238000010793 Steam injection (oil industry) Methods 0.000 claims abstract description 31
- 238000010276 construction Methods 0.000 claims abstract description 24
- 230000008878 coupling Effects 0.000 claims abstract description 16
- 238000010168 coupling process Methods 0.000 claims abstract description 16
- 238000005859 coupling reaction Methods 0.000 claims abstract description 16
- 239000007787 solid Substances 0.000 claims abstract description 16
- 230000008859 change Effects 0.000 claims abstract description 8
- 238000002474 experimental method Methods 0.000 claims abstract description 4
- 239000011435 rock Substances 0.000 claims description 29
- 238000004458 analytical method Methods 0.000 claims description 28
- 238000004364 calculation method Methods 0.000 claims description 19
- 230000035699 permeability Effects 0.000 claims description 16
- 238000012360 testing method Methods 0.000 claims description 10
- 239000004575 stone Substances 0.000 claims description 9
- 230000003993 interaction Effects 0.000 claims description 8
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 8
- 238000012423 maintenance Methods 0.000 claims description 7
- 238000009826 distribution Methods 0.000 claims description 6
- 239000011148 porous material Substances 0.000 claims description 6
- 238000005336 cracking Methods 0.000 claims description 5
- 230000001052 transient effect Effects 0.000 claims description 5
- 239000012530 fluid Substances 0.000 claims description 4
- 239000007788 liquid Substances 0.000 claims description 4
- 239000000243 solution Substances 0.000 claims description 4
- 150000001875 compounds Chemical class 0.000 claims description 3
- 238000004856 soil analysis Methods 0.000 claims description 3
- 230000005540 biological transmission Effects 0.000 claims description 2
- 238000007906 compression Methods 0.000 claims description 2
- 230000006835 compression Effects 0.000 claims description 2
- 238000012850 discrimination method Methods 0.000 claims description 2
- 230000005484 gravity Effects 0.000 claims description 2
- 238000012546 transfer Methods 0.000 claims description 2
- 230000007704 transition Effects 0.000 claims description 2
- 230000008016 vaporization Effects 0.000 claims description 2
- 238000009834 vaporization Methods 0.000 claims description 2
- 238000010008 shearing Methods 0.000 claims 2
- 238000011084 recovery Methods 0.000 abstract description 8
- 238000011156 evaluation Methods 0.000 abstract description 5
- 230000015556 catabolic process Effects 0.000 abstract description 3
- 238000006731 degradation reaction Methods 0.000 abstract description 3
- 239000007789 gas Substances 0.000 description 18
- 230000015572 biosynthetic process Effects 0.000 description 11
- 239000003921 oil Substances 0.000 description 11
- 238000011161 development Methods 0.000 description 5
- 239000000295 fuel oil Substances 0.000 description 4
- 238000002347 injection Methods 0.000 description 4
- 239000007924 injection Substances 0.000 description 4
- 239000011229 interlayer Substances 0.000 description 4
- VYPSYNLAJGMNEJ-UHFFFAOYSA-N Silicium dioxide Chemical compound O=[Si]=O VYPSYNLAJGMNEJ-UHFFFAOYSA-N 0.000 description 3
- 230000008901 benefit Effects 0.000 description 3
- 239000010410 layer Substances 0.000 description 3
- 230000007774 longterm Effects 0.000 description 3
- UFHFLCQGNIYNRP-UHFFFAOYSA-N Hydrogen Chemical compound [H][H] UFHFLCQGNIYNRP-UHFFFAOYSA-N 0.000 description 2
- 230000005465 channeling Effects 0.000 description 2
- 239000001257 hydrogen Substances 0.000 description 2
- 229910052739 hydrogen Inorganic materials 0.000 description 2
- 238000011065 in-situ storage Methods 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 description 2
- 238000012856 packing Methods 0.000 description 2
- 238000006467 substitution reaction Methods 0.000 description 2
- 238000010795 Steam Flooding Methods 0.000 description 1
- 239000002253 acid Substances 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 238000007792 addition Methods 0.000 description 1
- 230000005587 bubbling Effects 0.000 description 1
- OSMSIOKMMFKNIL-UHFFFAOYSA-N calcium;silicon Chemical compound [Ca]=[Si] OSMSIOKMMFKNIL-UHFFFAOYSA-N 0.000 description 1
- 239000010779 crude oil Substances 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000009472 formulation Methods 0.000 description 1
- 238000009533 lab test Methods 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 238000012821 model calculation Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000003345 natural gas Substances 0.000 description 1
- 239000003129 oil well Substances 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
- 239000003209 petroleum derivative Substances 0.000 description 1
- 238000012805 post-processing Methods 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 239000004576 sand Substances 0.000 description 1
- 239000000377 silicon dioxide Substances 0.000 description 1
- 239000002002 slurry Substances 0.000 description 1
- 239000002689 soil Substances 0.000 description 1
- 238000003860 storage Methods 0.000 description 1
Images
Classifications
-
- 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
- E21B33/00—Sealing or packing boreholes or wells
- E21B33/10—Sealing or packing boreholes or wells in the borehole
- E21B33/13—Methods or devices for cementing, for plugging holes, crevices or the like
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E10/00—Energy generation through renewable energy sources
- Y02E10/10—Geothermal energy
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Geology (AREA)
- Mining & Mineral Resources (AREA)
- Physics & Mathematics (AREA)
- Environmental & Geological Engineering (AREA)
- Fluid Mechanics (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geochemistry & Mineralogy (AREA)
- Earth Drilling (AREA)
- Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)
Abstract
本发明提供了一种热采井水泥环封固完整性的评价方法。该评价方法包括：确定注蒸汽热采井的热采施工参数；井筒内环境变化分析；实验确定温度对水泥石性能参数的影响；利用ABAQUS建立热传导模型；通过ABAQUS数值接口程序生成ODB文件；建立套管‑水泥环‑地层组合体的流固耦合模型；将热传导模型的ODB文件导入流固耦合模型中；对模型求解并分析数值结果。本发明的评价方法将水泥石强度衰退与高温高压的恶劣环境导致的水泥石破坏结合起来，对于评价热采井水泥环封固完整性具有非常重要的意义。
Description
技术领域
本发明涉及一种评价热采井水泥环固井封固完整性的方法，属于石油天然气勘探开发领域。
背景技术
固井时是向井内下入套管，并向井眼和套管之间的环形空间注入水泥的施工作业。固井后，套管-水泥环-地层组合体会胶结在一起，良好的胶结起到保护套管、支撑地层、封固油气层的作用，对安全高效生产有着重要的影响。但是，由于固井时的操作失误或是油气生产开发过程中的恶劣环境，导致水泥环封固失效情况的发生。如何提高套管与井壁间水泥环的密封质量、提高水泥环的强度时固井工程长期关注的问题。固井水泥环的力学性能，影响着固井的质量和油气井的寿命。现今全球剩余的石油资源中，稠油资源远远超过常规原油，稠油将是21世纪的重要能源之一。目前国内外稠油开采主要方式为蒸汽吞吐和蒸汽驱，热力采油时蒸汽温度高达380℃。在热采条件下，水泥石在高温下将不可避免的发生抗压强度衰退、渗透率增加等现象，水泥石的稳定性遭到破坏，井筒层间封隔效能失效，井口冒汽冒泡严重，导致层间气窜和环空带压，套损严重，大大缩短稠油井的生产寿命，影响稠油热采井开采效率，造成严重的经济损失，同时，酸性气田还可能对环境造成污染，给公众带来危害。根据辽河油田、克拉玛依等国内稠油生产区块现场数据统计，稠油井套损率高达50％，因水泥环层间封隔失效问题导致停产的井每年以10％的比例逐年上升，井口带压及井场周围冒汽冒泡、层间互窜严重。目前，对于常规地层油气井固井封固完整性的分析较多，而对于注汽井水泥环封固完整性的分析评价较少，这就对注气开发井建设和生产过程中风险的识别和控制造成阻碍，对于保证注气开发井固井作业安全及油气井长期高效开采，保证注气开发井固井水泥环的长期封隔稳定性存在很大的局限性。
发明内容
针对背景技术中存在的问题，本发明的目的是为了提供一种热采井高温高压环境下水泥环封固完整性的评价方法，该方法不仅对温度造成的水泥石性能下降进行分析，而且对外力作用下(温度、压力、注入流体等)水泥石本体失效和界面失效进行分析评价，从而建立了一套完整的力学体系。
为实现上述目的，本发明采用以下技术方案：
一种热采井水泥环封固完整性的评价方法，包括以下步骤：
1)获取热采井施工参数；
2)对水泥石进行不同温度条件养护并测试其性能，包括抗压强度、抗拉强度、弹性模量、泊松比；
3)基于1)、2)的数据，建立套管-水泥环-地层组合体的热传导模型；
4)将热传导模型计算得到ODB文件；
5)基于1)、2)的数据，建立套管-水泥环-地层组合体的热流固耦合模型；
6)对导入ODB文件的热流固耦合模型求解并分析数值结果，分析内容包括水泥石的温度分布、水泥石的周向应力、剪切应力以及一二界面是否出现微环隙。
其中步骤1)中各种参数可以由施工现场获取，并代入步骤3)和5)。步骤2)中可以通过实验室对水泥石测试获得，并代入步骤3)和5)。
优选的，所述步骤1)中热采井施工参数包括井口注蒸汽参数、井身轨迹参数、套管的属性参数、地层的属性参数、水泥石的属性参数以及目标层位载荷；所述井口注蒸汽参数包括注蒸汽压力、温度、速率、注蒸汽时间；所述井身轨迹参数包括套管尺寸、套管壁厚、环空尺寸、所钻井眼尺寸、目标层位深度；所述套管的属性参数包括密度、弹性模量、泊松比、热膨胀系数、热传导率；所述地层的属性参数包括密度、弹性模量、泊松比、热传导率、热膨胀系数、孔隙度、渗透率；所述水泥石的属性参数包括密度、热传导率、热膨胀系数、孔隙度、渗透率；所述目标层位载荷包括温度、压力；
优选的，所述步骤1)中目标层位载荷计算方法如下：
1)井筒内温度变化：将油管内的蒸汽温度看作与地面温度一致，即
T＝Tsur(t)(0＜t≤tend) (1)
其中，t为时间；T为蒸汽温度，℃；Tsur(t)为地面温度，℃；
地面温度指蒸汽地面注入时的温度。
2)井筒内压力降计算公式为：
其中，p-平均压力，pa；z-井深，m；d-输汽管线的内径，m；Ap为油管截面积，m2；is-蒸汽质量流量，kg/h；qg-蒸汽体积流量，m3/h。p、is、qg均为现场施工确定。
ρm为密度，kg/m3，其计算式为：
式中，M为气体分子量，无因次；R为通用气体常数，一般为R＝8.314J/(mol-1·℃-1)；γ为气体比重，无因次；Tf为气体温度，℃。
vm为速度，m/s，其计算式为：
式中，q为注蒸汽流量，m3/d，现场施工确定。
气体压缩因子Z的计算式为：
A1、A2、A3、A4、A5、A6、A7、A8分别取值为
A1＝0.31506237、A2＝-1.0467099、A3＝-0.57832729、A4＝0.53530771、A5＝-0.61232932、A6＝-0.1048813、A7＝0.68157001、A8＝0.68446549
式中，A1-A8的取值采用现今普遍认同的Dranchuk、Purris和Robinson的方法，通过经验公式确定；Tpr为井下某一位置处的温度，℃，由Tf迭代求得，或可简化为与Tf相等。
fm为湿蒸汽流体的摩擦阻力系数，其计算式为：
式中，e为管柱内壁绝对粗糙度，普通油管一般取4.57×10-5；
Re为气相雷诺数，无因次，计算式为
μm＝μs 1-H·μw H (8)
Vm＝XVs+(1-X)Vw (10)
式中，VT为套管体积；μs和μw分别为气态和液态水的粘度；Vw和Vs分别为注入的液态水和气态水的体积；Vm和H为中间过渡量；
X为干度，计算式为：
式中，
C1＝is(hs-hw) (12)
式中，Q为井筒热损失，W；hs为汽化潜热焓，hw为水的显热焓；采用W.S.Tortlke及S.M.Farouq Ail的计算方法；hm为总热焓，KJ/kg，计算方法为：
hm＝Xhs+(1-X)hw (15)
优选的，所述步骤2)中，首先利用常温养护箱(50℃)养护两天，转入高温养护釜养护两天，养护温度区间为50-380℃，使用压力机进行抗拉强度测试，利用巴西劈裂实验破坏载荷计算抗拉强度的方法为：
其中，σt为抗拉强度；D为试件直径；L为试件厚度；P为破坏载荷。
优选的，所述步骤3)、步骤5)中，模型建立不使用Tie约束，即允许组合体脱开。
优选的，所述步骤3)中，热传导模型的建立步骤为①建立部件；②施加属性；③装配部件；④施加分析步；⑤添加相互作用；⑥施加载荷；⑦划分网格；⑧提交作业。
其中，步骤①建立部件：建立三个不同的部件，分别为套管、水泥环和地层，地层厚度为套管-水泥环组合体的10倍。建立部件所需的套管尺寸、套管壁厚、环空尺寸、所钻井眼尺寸由现场井眼轨迹参数确定。
步骤②施加属性：套管的属性包括套管密度、套管弹性模量、套管泊松比、套管热膨胀系数、套管热传导率由施工现场给出。地层的属性包括岩石的密度、岩石弹性模量、岩石泊松比、岩石的热传导率、岩石热膨胀系数、岩石孔隙度、岩石渗透率，水泥石密度、水泥石的热传导率、水泥石的热膨胀系数、水泥石孔隙度、水泥石渗透率均由施工现场地层报告确定参数，现常用的1000米地层模型参数的数据如表1所示。水泥环的弹性模量、泊松比与温度变化的关系较大，因此由步骤2)给出，设置弹性参数时，选择使用与温度相关的数据选项，输入步骤2)得到的水泥石弹性模量和泊松比数据以及养护温度数值。
步骤③装配部件：选择创建部件选项，进行套管、水泥环、地层的非独立实例的创建。
步骤④施加分析步：施加瞬态热传递分析步。
步骤⑤添加相互作用：在套管与水泥环、水泥环与地层间添加热传导属性。
步骤⑥施加载荷：在套管内壁添加温度边界条件，并添加预定义温度条件。
步骤⑦划分网格：选择对象为部件，网格类型为DC3D8。
⑧提交作业。
这里所说的热传导模型的定义为：用软件建立不同材质的部件，对套管内壁施加温度，利用热能由高温传向低温的特性，使温度在不同材质中传播。
优选的，所述步骤5)中，热流固耦合模型的建立步骤为①建立部件；②施加属性；③装配部件；④施加分析步；⑤添加相互作用；⑥施加载荷；⑦划分网格；⑧提交作业。
步骤①建立部件：建立三个不同的部件，分别为套管、水泥环和地层，地层厚度为套管-水泥环组合体的10倍。建立部件所需的套管尺寸、套管壁厚、环空尺寸、所钻井眼尺寸由现场井眼轨迹参数确定。
步骤②施加属性：套管的属性包括套管密度、套管弹性模量、套管泊松比、套管热膨胀系数、套管热传导率由施工现场给出。地层的属性包括岩石的密度、岩石弹性模量、岩石泊松比、岩石的热传导率、岩石热膨胀系数、岩石孔隙度、岩石渗透率，水泥石密度、水泥石的热传导率、水泥石的热膨胀系数、水泥石孔隙度、水泥石渗透率均由施工现场地层报告确定参数，现常用的1000米地层模型参数的数据如表1所示。水泥环的弹性模量、泊松比与温度变化的关系较大，因此由步骤2)给出，设置弹性参数时，选择使用与温度相关的数据选项，输入步骤2)得到的水泥石弹性模量和泊松比数据以及养护温度数值。
步骤③装配部件：选择创建部件选项，进行套管、水泥环、地层的非独立实例的创建。
步骤④施加分析步：施加分析步2个，一为地应力平衡分析步，一为瞬态土分析步；在套管内壁施加载荷，计算方法如式1-16所示，大小如图8、图10、图12所示，由实际施工数据确定。
这里所说的地应力平衡分析步的目的是自动平衡地应力，瞬态土分析步的目的是把水泥石-地层看作可渗透的岩石。
步骤⑤添加相互作用：在套管与水泥环、水泥环与地层间添加热传导属性。
步骤⑥施加载荷：在套管内边界施加目标层位载荷，在地层外边界施加孔隙压力为0的外边界，并施加预定义孔隙比和预定义温度。
步骤⑦划分网格：选择对象为部件，套管网格类型为C3D8R，水泥石和地层网格类型为C3D8P。
步骤⑧提交作业。
这里所说的热流固耦合模型的定义为：建立一个新的流固耦合模型，该模型存在流体和固体间的相互作用，把上述热传导模型中的结果，即部件的每个节点上的温度导入新建立的模型中，使温度对流固耦合模型施加作用。
表1常用套管-水泥环-地层参数表
优选的，所述步骤6)中，定义剪切失效因子和径向开裂失效因子分别为ηc和ηt，公式分别为：
当ηc和ηt超过1时，水泥环分别发生剪切失效和径向开裂失效；
对于微环隙的判别方法，使用
一界面ro＝rif-roc＞0
二界面ri＝ric-rot＞0 (19)
式中，σi(i＝1,2,3)是水泥环中的最大、中间及最小主应力，MPa，步骤6中计算结果；是水泥石的内摩擦角，rad；σH为抗拉强度，步骤6中计算结果；σmax为实验结果，如图4所示，与水泥石养护温度有关；rif为套管中心至水泥环一界面的距离；roc为套管中心至套管外壁的距离；ric为套管中心至地层内边界的距离；rot为套管中心至水泥环外边界的距离，rif、roc、ric、rot结果均由步骤6)计算得出。
本发明同现有技术相比具有以下优点：
1)通过Mogi-Coulomb破坏准则将水泥石强度衰退与高温高压下水泥石破坏结合起来；
2)建立了热-流-固耦合模型，并能够正确的进行运算且得出相应的研究结果；3)本发明为一套完整的研究热采井固井封固完整性的方法，实现了只需利用现场注蒸汽参数，即可得到水泥环失效评价的方法。
附图说明
本发明的上述和/或附加的方面和优点从结合下面附图对实施例的描述中将变得明显和容易理解，其中：
图1为本发明的技术路线示意图。
图2为热-流-固耦合模型的技术路线示意图。
图3为实施例中的不同温度下水泥的抗压强度。
图4为实施例中的不同温度下水泥的抗拉强度。
图5为实施例中的不同温度下水泥的弹性模量。
图6为实施例中的不同温度下水泥的泊松比。
图7为实施例中注蒸汽速率与井筒干度的关系。
图8为实施例中注蒸汽速率与井筒压力的关系。
图9为实施例中注蒸汽压力与井筒干度的关系。
图10为实施例中注蒸汽压力与井筒压力的关系。
图11为实施例中注蒸汽温度与井筒干度的关系。
图12为实施例中注蒸汽温度与井筒压力的关系。
具体实施方式
为了对本发明的技术特征、目的和有益效果有更加清楚的理解，现对本发明的技术方案进行以下详细说明，但不能理解为对本发明的可实施范围的限定。
该实施例从施工现场获取热采井施工参数，包括井口注蒸汽参数(注蒸汽压力、温度、速率、注蒸汽时间)、井身轨迹参数(套管尺寸、套管壁厚、环空尺寸、所钻井眼尺寸、目标层位深度)、套管的属性参数(密度、弹性模量、泊松比、热膨胀系数、热传导率)、地层的属性参数(的密度、弹性模量、泊松比、热传导率、热膨胀系数、孔隙度、渗透率)、水泥石的属性参数(密度、热传导率、热膨胀系数、孔隙度、渗透率)。把获取的参数(或给定的数据，如表1)导入部件或属性，利用公式1-6计算目标层位载荷(包括温度、压力)，利用步骤2)获取水泥石弹性模量、泊松比参数。
将公式1-6进行求解计算，得到不同注蒸汽温度、不同注蒸汽压力、不同注蒸汽速率情况下井筒中压力和干度的变化，如图7-图12。
该实施例中采用的是目前常用的热采井固井水泥浆配方(嘉华G级水泥+硅砂)，钙硅比为1:1，水固比为0.44。使用50mm×50mm×50mm的标准模具对水泥石进行常温养护(50℃)，养护时间为2d，其目的是模拟固井时地层温度不高的情况。
待水泥石水化完全后，使用高温养护釜对水泥石进行高温养护2d，模拟热采井施工情况。养护温度分别为50℃、75℃、90℃、120℃、150℃、190℃、230℃、265℃、300℃、340℃、380℃。
养护完成后，使用取心机对水泥块进行取心处理，取心直径为25mm，取心长度为50mm。
对养护后的水泥块使用压力机测试其抗压强度，取心后利用巴西劈裂实验计算抗拉强度，计算公式如式7所示。利用三轴测试机测试不同养护温度条件下的弹性模量、泊松比，结果如图3-图6所示。
图3中，分段对不同温度条件下水泥的抗压强度进行拟合，得到抗压强度σc与温度T的关系为：
图4中，分段对不同温度条件下水泥的抗拉强度进行拟合，得到抗拉强度σH与温度T的关系为：
热传导模型建立：建立三个不同的部件，分别为套管、水泥环、地层，部件大小由现场获得的井身结构参数确定，包括套管尺寸、套管壁厚、环空尺寸、所钻井眼尺寸，地层厚度为套管-水泥环组合体的10倍，将上述得到的水泥石的弹性模量、泊松比与温度的关系数值以及现场得到的地层、套管和其他水泥环参数(或给定的数据，如表1)导入ABAQUS软件设为水泥石参数，设置水泥石弹性参数时，选择使用与温度相关的数据选项，输入步骤2)得到的弹性模量和泊松比数据以及养护温度数值。模拟高温水泥石强度衰退情况，建立热传导分析步，建立边界条件，施加温度边界条件，大小为步骤1)中施工现场给定的注蒸汽温度数据。分析步时间视注蒸汽时间而定，注蒸汽温度和压力施加部位为套管内壁，大小视模型计算结果而定。预定义温度设为地层实际温度。套管-水泥环-地层的单元均为DC3D8。
提交作业并运行完成后，在工作路径下找到热传导模型的ODB文件，此文件记录了套管-水泥环-地层组合体中不同时间不同单元和节点的温度分布，为后续热流固耦合模型的建立提供条件。
热流固耦合模型建立：同样的，建立三个不同的部件，分别为套管、水泥环、地层，部件大小由现场获得的井身结构参数确定，包括套管尺寸、套管壁厚、环空尺寸、所钻井眼尺寸，地层厚度为套管-水泥环组合体的10倍，将上述得到的水泥石的弹性模量、泊松比与温度的关系数值以及现场得到的地层、套管和其他水泥环参数(或给定的数据，如表1)导入ABAQUS软件设为水泥石参数，设置水泥石弹性参数时，选择使用与温度相关的数据选项，输入步骤2)得到的弹性模量和泊松比数据以及养护温度数值(式11、式12或图5、图6)。模拟高温水泥石强度衰退情况，建立两个不同的分析步，第一个分析步用于自动平衡地应力。第二个分析步为土体分析步，用于施加流固耦合模型，场输出选中NT(节点温度)、TEMP(单元温度)、COORD(节点坐标)等需要考虑的值。在第二个分析步中建立压力载荷，载荷大小为步骤1)中计算的目标层位压力数据，将地层外边界定义为孔隙压力边界，孔隙压力为0。在预定义场中，在第二个分析步中施加预定义温度场，温度场的区域选为全部部件，将ODB文件导入，此为热流固耦合模型最关键的一步。最后根据实际地层参数设置初始地应力并采用自动平衡地应力的方式平衡应力。套管单元为C3D8R，水泥环、地层的单元均为C3D8P，进行求解计算。
自动平衡地应力的操作方法是：在施加分析步模块时，选择地应力分析步，选择增量步的类型为自动，矩阵存储选择为非对称，直接迭代的方法，关闭几何非线性选项。
改变施工条件，如不同注蒸汽压力、不同注蒸汽时间、不同注蒸汽温度、不同注蒸汽速率；改变地层参数，如不同地层弹性模量、不同地层泊松比、不同地层渗透系数；考虑井底复杂工况，如不同水泥环厚度、不同套管偏心度。再次进行提交作业，以达到不同工况时热采井水泥环封固完整性的分析评价，需要注意的是，对提交的JOB命名应不同，以便于后处理的操作。
导出沿水泥环一界面至二界面的温度分布，导出一界面及二界面随时间变化的温度分布情况。通过ABAQUS后处理导出最大主应力、中间主应力、最小主应力分别为σ1、σ2、σ3，σc为不同温度下的抗压强度，如图3所示，剪切失效因子ηc计算公式为式(8)，其中，σc根据式11以及温度分布情况确定。导出周向应力S22为σH，σmax为不同温度下的抗拉强度，如图4所示。径向开裂失效因子ηt计算公式为式9，其中，σH根据式12以及温度分布情况确定。导出套管-水泥环界面和水泥环-地层界面的距离分别为ro和ri，界面失效的计算式为式10。当ηc和ηt超过1时，水泥环分别发生剪切失效和径向开裂失效；当ro和ri大于0时，水泥环发生界面失效。
此外，需要说明的是，本发明所属技术领域的技术人员可以对所描述的具体实施例做各种各样的修改或补充或采用类似的方式替代。应当理解，所属技术领域的普通技术人员在没有作出创造性劳动的前提下做出的其他实施例或者对部分技术特征进行等同替换，而不脱离本发明技术方案的精神，其均应涵盖在本发明请求保护的范围当中。
Claims (7)
1.一种热采井水泥环封固完整性的评价方法，其特征在于：包括以下步骤：
1)获取热采井施工参数；
2)对水泥石进行不同温度条件养护并测试其性能，包括抗压强度、抗拉强度、弹性模量、泊松比；
3)基于1)、2)的数据，建立套管-水泥环-地层组合体的热传导模型；
4)将热传导模型计算得到ODB文件；
5)基于1)、2)的数据，建立套管-水泥环-地层组合体的热流固耦合模型；
6)对导入ODB文件的热流固耦合模型求解并分析数值结果，分析内容包括水泥石的温度分布、水泥石的周向应力、剪切应力以及一二界面是否出现微环隙；
所述步骤3)中，热传导模型的建立步骤为：
步骤①建立部件：建立三个不同的部件，分别为套管、水泥环和地层，地层厚度为套管-水泥环组合体的10倍；
步骤②施加属性：套管的属性包括套管密度、套管弹性模量、套管泊松比、套管热膨胀系数、套管热传导率由施工现场给出；地层的属性包括岩石的密度、岩石弹性模量、岩石泊松比、岩石的热传导率、岩石热膨胀系数、岩石孔隙度、岩石渗透率，水泥石密度、水泥石的热传导率、水泥石的热膨胀系数、水泥石孔隙度、水泥石渗透率均由施工现场地层报告确定参数；水泥环的弹性模量、泊松比与温度变化的关系较大，由步骤2)给出，设置弹性参数时，选择使用与温度相关的数据选项，输入步骤2)得到的水泥石弹性模量和泊松比数据以及养护温度数值；
步骤③装配部件：选择创建部件选项，进行套管、水泥环、地层的非独立实例的创建；
步骤④施加分析步：施加瞬态热传递分析步；
步骤⑤添加相互作用：在套管与水泥环、水泥环与地层间添加热传导属性；
步骤⑥施加载荷：在套管内壁添加温度边界条件，并添加预定义温度条件；
步骤⑦划分网格：选择对象为部件，网格类型为DC3D8；
步骤⑧提交作业。
2.如权利要求1所述的一种热采井水泥环封固完整性的评价方法，其特征在于：所述步骤1)中热采井施工参数包括井口注蒸汽参数、井身轨迹参数、套管的属性参数、地层的属性参数、水泥石的属性参数以及目标层位载荷；所述井口注蒸汽参数包括注蒸汽压力、温度、速率、注蒸汽时间；所述井身轨迹参数包括套管尺寸、套管壁厚、环空尺寸、所钻井眼尺寸、目标层位深度；所述套管的属性参数包括密度、弹性模量、泊松比、热膨胀系数、热传导率；所述地层的属性参数包括密度、弹性模量、泊松比、热传导率、热膨胀系数、孔隙度、渗透率；所述水泥石的属性参数包括密度、热传导率、热膨胀系数、孔隙度、渗透率；所述目标层位载荷包括温度、压力。
3.如权利要求2所述的一种热采井水泥环封固完整性的评价方法，其特征在于：所述步骤1)中目标层位载荷计算方法如下：
1)井筒内温度变化：将油管内的蒸汽温度看作与地面温度一致，即
T＝Tsur(t)(0＜t≤tend) (1)
其中，t为时间；T为蒸汽温度，℃；Tsur(t)为地面温度，℃；
2)井筒内压力降计算公式为：
其中，p-平均压力，pa；z-井深，m；d-输汽管线的内径，m；Ap为油管截面积，m2；is-蒸汽质量流量，kg/h；qg-蒸汽体积流量，m3/h；p、is、qg均为现场施工确定；
ρm为密度，kg/m3，其计算式为：
式中，M为气体分子量，无因次；R为通用气体常数，一般为R＝8.314J/(mol-1·℃-1)；γ为气体比重，无因次；Tf为气体温度，℃；
vm为速度，m/s，其计算式为：
式中，q为注蒸汽流量，m3/d，现场施工确定；
气体压缩因子Z的计算式为：
A1、A2、A3、A4、A5、A6、A7、A8分别取值为
A1＝0.31506237、A2＝-1.0467099、A3＝-0.57832729、A4＝0.53530771、A5＝-0.61232932、A6＝-0.1048813、A7＝0.68157001、A8＝0.68446549
式中，A1-A8的取值采用现今普遍认同的Dranchuk、Purris和Robinson的方法，通过经验公式确定；Tpr为井下某一位置处的温度，℃，由Tf迭代求得，或可简化为与Tf相等；ρpr为拟对比密度；
fm为湿蒸汽流体的摩擦阻力系数，其计算式为：
式中，e为管柱内壁绝对粗糙度，普通油管一般取4.57×10-5；
Re为气相雷诺数，无因次，计算式为
μm＝μs 1-H·μw H (8)
Vm＝XVs+(1-X)Vw (10)
式中，VT为套管体积；μs和μw分别为气态和液态水的粘度；Vw和Vs分别为注入的液态水和气态水的体积；Vm和H为中间过渡量；
X为干度，计算式为：
式中，
C1＝is(hs-hw) (12)
式中，Q为井筒热损失，W；hs为汽化潜热焓，hw为水的显热焓；采用W.S.Tortlke及S.M.FarouqAil的计算方法；hm为总热焓，KJ/kg，计算方法为：
hm＝xhs+(1-x)hw (15)。
4.如权利要求1所述的一种热采井水泥环封固完整性的评价方法，其特征在于：所述步骤2)中，首先利用常温养护箱50℃养护两天，转入高温养护釜养护两天，养护温度区间为50-380℃，使用压力机进行抗拉强度测试，利用巴西劈裂实验破坏载荷计算抗拉强度的方法为：
其中，σt为抗拉强度；D为试件直径；L为试件厚度；P为破坏载荷。
5.如权利要求1所述的一种热采井水泥环封固完整性的评价方法，其特征在于：所述步骤3)、步骤5)中，模型建立不使用Tie约束，即允许组合体脱开。
6.如权利要求1所述的一种热采井水泥环封固完整性的评价方法，其特征在于：所述步骤5)中，热流固耦合模型的建立步骤为：
步骤①建立部件：建立三个不同的部件，分别为套管、水泥环和地层，地层厚度为套管-水泥环组合体的10倍；
步骤②施加属性：套管的属性包括套管密度、套管弹性模量、套管泊松比、套管热膨胀系数、套管热传导率由施工现场给出；地层的属性包括岩石的密度、岩石弹性模量、岩石泊松比、岩石的热传导率、岩石热膨胀系数、岩石孔隙度、岩石渗透率，水泥石密度、水泥石的热传导率、水泥石的热膨胀系数、水泥石孔隙度、水泥石渗透率均由施工现场地层报告确定参数；水泥环的弹性模量、泊松比与温度变化的关系较大，由步骤2)给出，设置弹性参数时，选择使用与温度相关的数据选项，输入步骤2)得到的水泥石弹性模量和泊松比数据以及养护温度数值；
步骤③装配部件：选择创建部件选项，进行套管、水泥环、地层的非独立实例的创建；
步骤④施加分析步：施加分析步2个，一为地应力平衡分析步，一为瞬态土分析步；在套管内壁施加载荷；
步骤⑤添加相互作用：在套管与水泥环、水泥环与地层间添加热传导属性；
步骤⑥施加载荷：在套管内边界施加目标层位载荷，在地层外边界施加孔隙压力为0的外边界，并施加预定义孔隙比和预定义温度；
步骤⑦划分网格：选择对象为部件，套管网格类型为C3D8R，水泥石和地层网格类型为C3D8P；
步骤⑧提交作业。
7.如权利要求1所述的一种热采井水泥环封固完整性的评价方法，其特征在于：所述步骤6)中，定义剪切失效因子和径向开裂失效因子分别为ηc和ηt，公式分别为：
当ηc和ηt超过1时，水泥环分别发生剪切失效和径向开裂失效；
对于微环隙的判别方法，使用
一界面ro＝rif-roc＞0
二界面ri＝ric-rot＞0 (19)
式中，σi(i＝1,2,3)是水泥环中的最大、中间及最小主应力，MPa，步骤6中计算结果；
是水泥石的内摩擦角，rad；σH为抗拉强度，步骤6中计算结果；σmax为实验结果，与水泥石养护温度有关；rif为套管中心至水泥环一界面的距离；roc为套管中心至套管外壁的距离；ric为套管中心至地层内边界的距离；rot为套管中心至水泥环外边界的距离，rif、roc、ric、rot结果均由步骤6计算得出。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910930508.4A CN110543745B (zh) | 2019-09-29 | 2019-09-29 | 一种热采井水泥环封固完整性的评价方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910930508.4A CN110543745B (zh) | 2019-09-29 | 2019-09-29 | 一种热采井水泥环封固完整性的评价方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110543745A CN110543745A (zh) | 2019-12-06 |
| CN110543745B true CN110543745B (zh) | 2023-04-07 |
Family
ID=68714911
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910930508.4A Active CN110543745B (zh) | 2019-09-29 | 2019-09-29 | 一种热采井水泥环封固完整性的评价方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110543745B (zh) |
Families Citing this family (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111859701A (zh) * | 2020-07-29 | 2020-10-30 | 西南石油大学 | 一种固井水泥浆塑性体积收缩与孔隙压降评价方法 |
| CN112163385B (zh) * | 2020-10-30 | 2023-09-22 | 云翼超算（北京）软件科技有限公司 | 求解强热流固耦合问题的并行无网格方法及*** |
| CN113688533B (zh) * | 2021-09-01 | 2022-08-26 | 西南石油大学 | 一种预测固井注水泥环空当量循环密度的计算方法 |
| CN115992669A (zh) * | 2021-10-20 | 2023-04-21 | 中国石油天然气股份有限公司 | 一种热采井水泥浆配方开发方法 |
| CN117235839B (zh) * | 2022-09-16 | 2024-04-19 | 中国葛洲坝集团三峡建设工程有限公司 | 基于自平衡应力场分析的井筒水泥环安全载荷计算方法 |
| CN117371822B (zh) * | 2023-12-04 | 2024-02-23 | 中国石油大学(华东) | 保障深水钻遇水合物稳定的固井调控、评价指标建立方法 |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104405366A (zh) * | 2014-10-22 | 2015-03-11 | 西南石油大学 | 一种高温高压固井水泥环力学完整性测试装置及方法 |
| CN104500034A (zh) * | 2014-12-29 | 2015-04-08 | 西南石油大学 | 一种评价压力变化对水泥环完整性影响的装置及方法 |
| CN105089639A (zh) * | 2014-04-23 | 2015-11-25 | 中国石油化工股份有限公司 | 一种煤层气井井筒流动动态预测方法 |
| CN106771096A (zh) * | 2016-12-09 | 2017-05-31 | 西南石油大学 | 一种固井水泥环封隔能力动态测试装置及实验方法 |
| CN109538185A (zh) * | 2018-10-30 | 2019-03-29 | 中国海洋石油集团有限公司 | 一种耦合温度场条件下的多层套管井井筒完整性分析模型 |
-
2019
- 2019-09-29 CN CN201910930508.4A patent/CN110543745B/zh active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105089639A (zh) * | 2014-04-23 | 2015-11-25 | 中国石油化工股份有限公司 | 一种煤层气井井筒流动动态预测方法 |
| CN104405366A (zh) * | 2014-10-22 | 2015-03-11 | 西南石油大学 | 一种高温高压固井水泥环力学完整性测试装置及方法 |
| CN104500034A (zh) * | 2014-12-29 | 2015-04-08 | 西南石油大学 | 一种评价压力变化对水泥环完整性影响的装置及方法 |
| CN106771096A (zh) * | 2016-12-09 | 2017-05-31 | 西南石油大学 | 一种固井水泥环封隔能力动态测试装置及实验方法 |
| CN109538185A (zh) * | 2018-10-30 | 2019-03-29 | 中国海洋石油集团有限公司 | 一种耦合温度场条件下的多层套管井井筒完整性分析模型 |
Non-Patent Citations (2)
| Title |
|---|
| 基于可靠性理论的复杂深井井筒***完整性研究;赵新波;《中国优秀博硕士学位论文全文数据库(硕士) 工程科技I辑》;20190715(第07期);全文 * |
| 考虑热固耦合作用的高温高压井井筒完整性分析;赵新波等;《工程科学学报》;20160201;第38卷(第1期);全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN110543745A (zh) | 2019-12-06 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110543745B (zh) | 一种热采井水泥环封固完整性的评价方法 | |
| Zhang et al. | Effect of thermal expansion annulus pressure on cement sheath mechanical integrity in HPHT gas wells | |
| Therond et al. | Large scale testing & modelling for cement zonal isolation of water injection wells | |
| Williams et al. | Flexible, expanding cement system (FECS) successfully provides zonal isolation across Marcellus shale gas trends | |
| McDaniel et al. | Cement sheath durability: increasing cement sheath integrity to reduce gas migration in the marcellus shale play | |
| CN109372466B (zh) | 利用天然地温场实现液-固-液相态转化的暂堵转向方法 | |
| Barreda et al. | Impact of cyclic pressure loading on well integrity in multi-stage hydraulic fracturing | |
| Yan et al. | Numerical investigation of debonding extent development of cementing interfaces during hydraulic fracturing through perforation cluster | |
| CN113756744A (zh) | 一种深水高温高压井套管-水泥环-地层***密封失效的数值模拟方法 | |
| Orlic et al. | Numerical investigations of cement interface debonding for assessing well integrity risks | |
| Han et al. | Study on the whole life cycle integrity of cement interface in heavy oil thermal recovery well under circulating high temperature condition | |
| CN111859701A (zh) | 一种固井水泥浆塑性体积收缩与孔隙压降评价方法 | |
| Wang et al. | Study on sealing integrity of cement sheath in ultra deep well reservoir reconstruction | |
| Tong et al. | Application of multiphysics coupling FEM on open wellbore shrinkage and casing remaining strength in an incomplete borehole in deep salt formation | |
| Miller et al. | Advancements in annular-pressure-buildup mitigation for thunder horse wells | |
| Khandka | Leakage behind casing | |
| Wang et al. | Study on simulation of cement plug—formation interface stripping failure and main influencing factors | |
| Li et al. | A Practical Approach to Maximise UGS Storage Efficiency in Order to Secure Energy Supply | |
| Khodami et al. | Numerical Investigation of the Impact of Geomechanical Parameters of Formations on Well Integrity of One of the Iranian Oil Fields | |
| Petty et al. | Life cycle modeling of wellbore cement systems used for enhanced geothermal system development | |
| Li et al. | Cement sheath integrity in anisotropic shale formations: a numerical investigation | |
| WU et al. | Optimum drilling design of high deviated directional wells and its application | |
| Johnson et al. | Energy Transition by Employing a Self-Healing, Reduced Carbon Dioxide Footprint Sealant in a Strategic Underground Gas Storage Project | |
| Speer et al. | Cemented Wellbore Experiments Reveals That Cement Interface Debonding Modeling Overestimate Potential Leakage Rates | |
| Krusche et al. | Application of engineered cementing solution to solve long-term cement integrity issues in Tunisia |
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 |