一种水合物地层固井过程高压气水反侵的判别方法技术

本发明专利技术公开了一种水合物地层固井过程高压气水反侵的判别方法，包括：建立水合物地层数值模型；将水合物地层数值模型进行网格划分，划分为若干单元；设计各项固井工艺参数的取值，采用该取值进行模拟试验，判别该固井工艺参数取值下是否发生反侵。本发明专利技术的优点是：能够高度还原水合物地层固井过程；能够判别水合物地层固井中高压气水是否反侵入水泥环。

【技术实现步骤摘要】
一种水合物地层固井过程高压气水反侵的判别方法
本专利技术涉及天然气水合物开发
，特别涉及一种天然气水合物地层固井过程高压气水是否反侵环空水泥浆的判别方法。
技术介绍
天然气水合物是由水分子和天然气分子(主要是CH4、C2H6、CO2等)在低温和高压条件下组成的似冰状非化学计量笼形化合物。由于具有资源储量大、主要燃烧产物环保无污染等优点，在能源、环境和全球气候变化中的重要意义使其成为当前地球科学和能源研究的一大热点。在深水油气钻遇水合物地层以及水合物储层开采过程中，固井是钻完井作业过程中一个重要环节，它包括下套管和注水泥，其主要目的是保护和支撑油气井内的套管，封隔油、气和水等地层。对于稳定井壁、固定套管和延长油气井生产寿命等方面具有重要意义。天然气水合物地层通常具有地层温度低和孔隙压力高的特点。同时，天然气水合物储层通常胶结性较差(甚至未胶结)、力学强度弱，且水合物所处的相平衡状态较为脆弱，轻微的温升就会导致水合物的分解并引发链锁反应。而在固井过程中，固井水泥浆的水化放热会引起近井壁一定范围内地层温度的升高，极易导致水合物的分解，从而形成局部高压游离气、水区域。而高压游离气和水在压差的驱动下易反向侵入水泥浆中，从而对水泥环的内部结构、力学强度、二界面胶结强度和储层封隔能力产生不利影响，甚至会导致二界面缩胀开裂、井壁垮塌、水气窜流和固井报废等严重钻井事故。因此，寻求一种针对该反向侵入是否发生的判别手段显得十分重要。由于水合物现场或室内试验成本高，技术难度大，耗时耗力，选取数值模拟方法对其进行研究更加合理。但水合物地层固井过程的复杂性为数值模拟带来了一些难点，主要包括：(1)在水泥浆初凝前，侵入地层的水泥浆作为移动热源，边侵入边继续发生水化反应放热。(2)水泥浆放热速率是不断变化的。由于数值模拟手段的局限性，目前尚未拥有模拟水合物地层固井过程的类似方法。
技术实现思路
本专利技术针对现有技术的缺陷，提供了一种水合物地层固井过程高压气水反侵的判别方法，解决了现有技术中存在的缺陷。为了实现以上专利技术目的，本专利技术采取的技术方案如下：一种水合物地层固井过程高压气水反侵的判别方法，包括以下步骤：步骤1，建立水合物地层数值模型，水合物地层数值模型内应包含地层、套管、套管内部空间、套管与地层间的环状空间；步骤2，将水合物地层数值模型进行网格划分，划分为若干单元，若干单元之间进行物质交换，物质交换或是能量交换，或是物质及能量交换；步骤3，设计各项固井工艺参数的取值，取值根据实际地层情况及施工工况共同确定，采用该组固井工艺参数取值进行模拟试验，根据模拟试验的结果，判别该固井工艺参数下是否发生反侵。进一步地，步骤1中水合物地层数值模型中包括：完整的地层物理性质参数，地层物理性质参数包括：地层孔隙压力、地层骨架密度、地层孔隙度、地层绝对渗透率、地层骨架比热容、地层导热系数、地层温度和地层孔隙水合物饱和度。进一步地，步骤2中，若干单元之间的物质交换，或是能量交换，或是物质及能量交换，原理如下：模型内的质量和能量守恒可通过有限体积积分差的方式来实现：式中：V,Vn——体积，子域n的体积[L3]；Mκ——组分κ的质量积累项[kgm-3]；A,Γn——表面积，子域n的表面积[L2]；Fκ——组分κ的达西通量向量[kgm-2s-1]；n——向内单位法向量；qκ——组分κ的源/汇项[kgm-3s-1]；t——时间[s]。β包含H(水合物相)，A(液相)，G(气相)，I(冰相)；κ包含w(水),m(天然气),h(水合物),i(水溶性抑制剂),θ(热量)。其中：①水合物相平衡曲线式中：Pe为水合物相平衡压力，MPa；Te为相平衡温度，K。②质量累计项Mκ式中：φ——孔隙度；Sβ——相β的饱和度；ρβ——相β的密度[kg/m3]；Xβκ——相β中组分κ的质量分数。③热量累积项Mθ式中：ρR——岩石密度[kgm-3]；CR——干岩的热容[Jkg-1K-1]；Uβ——相β的比内能[Jkg-1]；Δ()——当前时间步长括号内的数量变化；ΔUH——水合物解离/形成的比焓[Jkg-1]④流量项Fκ水、甲烷和抑制剂的质量流量取决于液相和气相的变化，由于两个固相(β≡H，I)不会对流体流量产生影响，所以在采用动力学模型时，水合物组分经过所有单元边界时的质量流量为0。液相：式中：k——岩土体绝对渗透率[m2]；krA——为液相的相对渗透率；μA——液相粘度[Pa·s]；PA——液相压力[Pa]；g——重力加速度矢量[m·s-2]。气相：式中：k0——高气压条件下的绝对渗透率[m2]；b——气体滑脱因子[Pa]；krG——气相渗透率；μG——气相粘度[Pa·s]；JκG——组分κ在气相中的扩散质量流量。热流量：式中：KR——岩土体热导率[W·m-1·K-1]；Kβ——相β的热导率[W·m-1·K-1]；hβ——相β的比焓[J·kg-1]；fσ——热辐射因子；σ0——Stefan-Boltzmann常数，其值为5.6687×10-8J·m-2·K-4。⑤源/汇项qκ地层流体开采或注入抑制剂：式中：qβ——相β的开采或注入速率[kg·m-3]。地层注入热量或抽取热量：对于方程(1)在空间上的离散而言，可采用有限体积积分差的方式来实现：①体积积分式中：M是体积归一化广延量；Mn是在体积Vn里的平均值。②表面积分式中：Fnm是F通过体积单元Vn和Vm之间界面Anm法向(向内分量的平均值)。在实际求解时，通过Newton-Raphson迭代法进行数值计算。进一步地，步骤3中，固井工艺参数主要包括：水泥浆水化放热速率、固井压差、保压时间、水泥浆密度、水泥初凝时间等。进一步地，步骤3中，模拟试验包括以下步骤：步骤31，在水合物地层数值模型中设置固井工艺参数，先通过环状空间的体积和水泥浆此时的放热速率来计算整个环空空间单元格的放热量，再将放热量设置在相应单元中；步骤32，模拟开始时，在环空空间单元处增加等同于固井压差的压力值，保压一段时间之后再减去这部分压力，固井压差及保压时间根据实际地层情况和施工工况确定；步骤33，整个模拟试验从施加压力开始，到卸去压力，再到水泥浆初凝为止。将整个模拟在时间上分为若干个段，在第一段模拟结束后，根据该段模拟结果得到水泥浆的侵入量，之后依次计算紧靠环空空间的若干地层单元中的孔隙体积，并根据实际工况或经验假定水泥浆侵入过程中会充填孔隙体积的百分数，得出水泥浆的侵入深度和在相应单元内的体积，由此配合水泥浆此时的放热速率计算出各单元内此刻的放热量，下一段模拟开始时将各单元的放热量设置完成后继续进行模拟，以此循环直到最后一段本文档来自技高网...

【技术保护点】
1.一种水合物地层固井过程高压气水反侵的判别方法，其特征在于，包括以下步骤：/n步骤1，建立水合物地层数值模型，水合物地层数值模型内应包含地层、套管、套管内部空间、套管与地层间的环状空间；/n步骤2，将水合物地层数值模型进行网格划分，划分为若干单元，若干单元之间可以进行物质交换，或是能量交换，或是物质及能量交换；/n步骤3，设计各项固井工艺参数的取值，取值根据实际地层情况及施工工况共同确定，采用该组固井工艺参数取值进行模拟试验，根据模拟试验的结果，判别该固井工艺参数下是否发生反侵。/n

【技术特征摘要】
1.一种水合物地层固井过程高压气水反侵的判别方法，其特征在于，包括以下步骤：
步骤1，建立水合物地层数值模型，水合物地层数值模型内应包含地层、套管、套管内部空间、套管与地层间的环状空间；
步骤2，将水合物地层数值模型进行网格划分，划分为若干单元，若干单元之间可以进行物质交换，或是能量交换，或是物质及能量交换；
步骤3，设计各项固井工艺参数的取值，取值根据实际地层情况及施工工况共同确定，采用该组固井工艺参数取值进行模拟试验，根据模拟试验的结果，判别该固井工艺参数下是否发生反侵。


2.根据权利要求1所述的判别方法，其特征在于：步骤1中水合物地层数值模型中包括：完整的地层物理性质参数，地层物理性质参数包括：地层孔隙压力、地层骨架密度、地层孔隙度、地层绝对渗透率、地层骨架比热容、地层导热系数、地层温度和地层孔隙水合物饱和度。


3.根据权利要求1所述的判别方法，其特征在于：步骤2中，若干单元之间的物质交换，或是能量交换，或是物质及能量交换，原理如下：
模型内的质量和能量守恒可通过有限体积积分差的方式来实现：



式中：V,Vn——体积，子域n的体积[L3]；Mκ——组分κ的质量积累项[kgm-3]；A,Γn——表面积，子域n的表面积[L2]；Fκ——组分κ的达西通量向量[kgm-2s-1]；n——向内单位法向量；qκ——组分κ的源/汇项[kgm-3s-1]；t——时间[s]；
β包含H(水合物相)，A(液相)，G(气相)，I(冰相)；
κ包含w(水),m(天然气),h(水合物),i(水溶性抑制剂),θ(热量)；
其中：
①水合物相平衡曲线






式中：Pe为水合物相平衡压力，MPa；Te为相平衡温度，K；
②质量累计项Mκ



式中：——孔隙度；Sβ——相β的饱和度；ρβ——相β的密度[kg/m3]；Xβκ——相β中组分κ的质量分数；
③热量累积项Mθ






式中：ρR——岩石密度[kgm-3]；CR——干岩的热容[Jkg-1K-1]；Uβ——相β的比内能[Jkg-1]；Δ()——当前时间步长括号内的数量变化；ΔUH——水合物解离/形成的比焓[Jkg-1]
④流量项Fκ
水、甲烷和抑制剂的质量流量取决于液相和气相的变化，由于两个固相(β≡H，I)不会对流体流量产生影响，所以在采用动力学模型时，水合物组分经过所有单元边界时的质量流量为0；
液相：



式中：k——岩土体绝对渗透率[m2]；krA——为液相的相对渗透率；μA——液相粘度[Pa·s]；PA——液相压力[Pa]；g——重力加速度矢量[m·s-2]；
气相：
...

【专利技术属性】
技术研发人员：郑明明，王晓宇，王凯，周珂锐，吴喆，刘天乐，郑少军，李丽霞，
申请(专利权)人：成都理工大学，
类型：发明
国别省市：四川;51