1.本发明涉及油气田开发评价领域，更具体地说，涉及一种水平井分段压裂温度分布预测方法。
背景技术：2.在过去几十年中，人们设计了监测和诊断裂缝效果和油井性能的技术，如非放射性示踪剂和生产测井。其技术具有局限性，基于现有技术的裂缝模拟分析具有很大的不确定性，近年来，光纤分布式温度传感器(dts)在裂缝处理中得到了广泛的应用，能实时监测压裂过程中的动态温度场。温度数据与压裂诊断的其它井信息一起使用，可以提供识别裂缝位置和流速剖面的信息。
3.近年来，光纤分布式温度传感器(dts)越来越多地用于水平井压裂动态监测，拟解决在水平井分段压裂过程中普遍面临的人工裂缝起裂位置不明、压裂液去向未知、裂缝扩展形态不清、压裂效果难以评价等技术难题，而温度预测模型是基于dts监测进行压裂诊断的基础，但目前定量预测水平井压裂过程中的温度分布仍是一项巨大挑战。
技术实现要素：4.本发明要解决的技术问题在于，供一种水平井分段压裂温度分布预测方法，通过使用温度预测模型，模拟温度分布来识别裂缝数量，并评价压裂井性能。
5.本发明解决其技术问题所采用的技术方案是：构造一种水平井分段压裂温度分布预测方法，包括以下步骤：
6.s1、建立水平井井筒流动模型：
[0007][0008]
式中，t为时间；r为井筒内径；γ为井筒打开程度；ρ
l
为井筒流入流体密度；v
l
井筒中注入液速度；ρ
wb
为井筒中的流体密度；
[0009]
s2、设定井筒流动模型的初始和边界条件：
[0010]
vi＝v|
t＝0
＝0
ꢀꢀꢀ
(2)
[0011][0012]
s3、建立水平井井筒热模型，瞬态井筒热模型的最终形式为：
[0013][0014]
式中，c
p
为比热容；θ为水平井筒倾角；tr表示储层温度；由储层热模型计算，t
l
表示注入液温度；g为重力加速度；u
t
为综合传热系数；t为温度；
[0015]
s4、设定水平井井筒热模型的初始和边界条件：
[0016]
t|
t＝0
＝t
wf
,
ꢀꢀꢀ
(5)
[0017][0018]
s5、对裂缝几何形态进行如下假设：(1)发生在裂缝处的流体滤失与裂缝面呈线性垂直；(2)注入过程中裂缝具有无限的导流能力；(3)压裂液不可压缩；并对裂缝扩展过程进行了建模，裂缝的流动模型为：
[0019][0020]
式中，τ(x)是裂缝延伸长度x的时间之后，流体滤失开始；xf是裂缝半长，m；w、h分别为任意微元段的宽和高；c是考虑流体压缩性效应cc、流体粘度效应cv和造壁性效应cw的总流体滤失系数；
[0021]
s6、设定裂缝流动模型的初始和边界条件为：
[0022]
xf|
t＝0
＝0
ꢀꢀꢀ
(8)
[0023][0024]
s7、建立裂缝热学模型，计算裂缝扩展和关井过程中裂缝内压裂液温度分布；能量守恒的概念与体积微分方程上的裂缝流动模型相似：
[0025][0026]
s8、设定裂缝热学模型的初始和边界条件为:
[0027]
t
l
|
t＝0
＝tiꢀꢀꢀ
(11)
[0028][0029]
式中，ti表示初始储层温度；t为注入流体温度；vi为x＝r和γ＝1时的流体速度；由井筒模型计算得出，在储层温度模型中，裂缝面边界与储层模型耦合；
[0030]
s9、压裂液从裂缝面向储层滤失的距离可通过将滤失速度与滤失时间相结合来计算，储层流动模型为：
[0031][0032]
s10、建立了储层热学模型，并与井筒和裂缝模型相结合，根据能量平衡方程推导局部地层能量守恒方程，储层热学模型为：
[0033][0034]
式(14)为局部地层能量守恒方程；对于被压裂液侵入的网格块，其平均有效热容和有效热导率不是恒定的；上述方程左端为随时间t在微元段中的能量积累：右端第一项是由于流体滤失进入地层而产生的热转换；第二和第三项是二维流动系统中的热传导；其中，ke是平均有效储层导热系数；tr表示地层岩石温度；
[0035]
s11、设定储层流动模型和储层热学模型的初始和边界条件为：
[0036]
tr|
t＝0
＝tiꢀꢀꢀ
(18)
[0037][0038]
s12、将上述步骤中的井筒、裂缝和储层模型耦合并结合边界条件进行求解，获取实时模拟井筒温度分布。
[0039]
按上述方案，通过dts测量压裂注液过程中动态温度数据，从而进行裂缝诊断，实现裂缝位置识别、判断裂缝半长大小，并对压裂段内形成的裂缝数量和位置进行准确诊断。
[0040]
实施本发明的水平井分段压裂温度分布预测方法，具有以下有益效果：
[0041]
1、本发明通过dts测量压裂注液过程中动态温度数据，从而进行裂缝诊断，如裂缝位置识别、裂缝半长大小；
[0042]
2、通过本发明可对压裂段内形成的裂缝数量和位置进行准确诊断；
[0043]
3、本发明在考虑多种微量热效应的基础之上，实现了水平井分段压裂注液过程中的温度分布动态模拟预测。
附图说明
[0044]
下面将结合附图及实施例对本发明作进一步说明，附图中：
[0045]
图1是本发明水平井分段压裂温度分布预测方法的技术路线图；
[0046]
图2是本发明水平井分段压裂温度分布预测模型的求解步骤示意图；
[0047]
图3是裂缝温度分布示意图；
[0048]
图4是储层温度分布示意图；
[0049]
图5是压裂段内分别形成1条、2条裂缝时的温度分布云图(从左至右)。
具体实施方式
[0050]
为了对本发明的技术特征、目的和效果有更加清楚的理解，现对照附图详细说明本发明的具体实施方式。
[0051]
本发明提供了一种水平井分段压裂温度分布预测方法，如图1所示，该方法包括下列步骤：
[0052]
s1、建立水平井井筒流动模型：
[0053][0054]
式中，t为时间，s；r为井筒内径，m；γ为井筒打开程度；ρl为井筒流入流体密度，kg/m3；vl井筒中注入液速度，m/s；ρ
wb
为井筒中的流体密度，kg/m3。
[0055]
s2、设定井筒流动模型的初始和边界条件：
[0056]
vi＝v|
t＝0
＝0
ꢀꢀꢀ
(2)
[0057][0058]
s3、建立水平井井筒热模型，瞬态井筒热模型的最终形式为：
[0059][0060]
式中，c
p
为比热容，j/(kg
·
k)；θ为水平井筒倾角，
°
；tr表示储层温度，℃；由储层热模型计算，t
l
表示注入液温度，℃；g为重力加速度，m/s2；u
t
为综合传热系数，w/(m2·
k)；t为温度，℃。
[0061]
s4、设定水平井井筒热模型的初始和边界条件：
[0062]
t|
t＝0
＝t
wf
,
ꢀꢀꢀ
(5)
[0063][0064]
s5、为了简单起见，对裂缝几何形态进行如下假设：(1)发生在裂缝处的流体滤失与裂缝面呈线性垂直；(2)注入过程中裂缝具有无限的导流能力；(3)压裂液不可压缩。并对裂缝扩展过程进行了建模，裂缝的流动模型为：
[0065][0066]
式中，τ(x)是裂缝延伸长度x的时间之后，流体滤失开始。xf是裂缝半长，m；w、h分别为任意微元段的宽和高，m；c是考虑流体压缩性效应cc、流体粘度效应cv和造壁性效应cw的总流体滤失系数。
[0067]
s6、设定裂缝流动模型的初始和边界条件为：
[0068]
xf|
t＝0
＝0
ꢀꢀꢀ
(8)
[0069][0070]
s7、建立了裂缝热学模型，计算了裂缝扩展和关井过程中裂缝内压裂液温度分布。能量守恒的概念与体积微分方程上的裂缝流动模型相似。
[0071][0072]
s8、设定裂缝热学模型的初始和边界条件为:
[0073]
t
l
|
t＝0
＝tiꢀꢀꢀ
(11)
[0074][0075]
式中，ti表示初始储层温度，℃；t为注入流体温度，℃。vi为x＝r和γ＝1时的流体速度，m/s；由井筒模型计算得出，在储层温度模型中，裂缝面边界与储层模型耦合。
[0076]
s9、压裂液从裂缝面向储层滤失的距离可通过将滤失速度与滤失时间相结合来计算，因此，储层流动模型为：
[0077]
[0078]
s10、为了研究压裂过程中的储层温度，建立了储层热学模型，并与井筒和裂缝模型相结合。根据能量平衡方程推导了局部地层能量守恒方程，储层热学模型为：
[0079][0080]
式(14)为局部地层能量守恒方程。对于被压裂液侵入的网格块，其平均有效热容和有效热导率不是恒定的。上述方程左端为随时间t在微元段中的能量积累：右端第一项是由于流体滤失进入地层而产生的热转换。第二和第三项是二维流动系统中的热传导。其中，ke是平均有效储层导热系数，(j/(m
·s·
k)；tr表示地层岩石温度，℃。
[0081]
s11、设定储层流动模型和储层热学模型的初始和边界条件为：
[0082]
tr|
t＝0
＝tiꢀꢀꢀ
(18)
[0083][0084]
s12、为了实时模拟井筒温度分布，将上述步骤中的井筒、裂缝和储层模型耦合并结合边界条件进行求解。
[0085]
通过模拟一口水平井分段压裂过程中的温度动态，分析水平井分段压裂过程中温度分布特征如图3和4所示，分别为该水平井压裂60分钟时储层和裂缝中的温度分布，验证了本文建立温度模型的可靠性。
[0086]
温度预测模型可模拟分段多簇压裂时的温度分布，考虑到多簇裂缝同时扩展。为了研究注入流量对水平井分段多簇压裂温度分布的影响，将总注入量分配到压裂段的每个射孔簇，模拟结果如图5所示。
[0087]
在理想条件下，段内形成单裂缝时，比段内形成多条裂缝时的温度响应更明显，而段内多簇起裂时，凡是有压裂液注入的位置，储层温度分布上都会形成一个“尖锥”，而“尖锥”的位置，就指示着裂缝起裂位置，对应地，井筒温度剖面中的温降清晰地指示了裂缝起裂的位置，据此，可以通过水平井压裂注液过程中的井筒温度剖面监测，对压裂段内形成的裂缝数量及进行准确诊断。
[0088]
上面结合附图对本发明的实施例进行了描述，但是本发明并不局限于上述的具体实施方式，上述的具体实施方式仅仅是示意性的，而不是限制性的，本领域的普通技术人员在本发明的启示下，在不脱离本发明宗旨和权利要求所保护的范围情况下，还可做出很多形式，这些均属于本发明的保护之内。
技术特征：1.一种水平井分段压裂温度分布预测方法，其特征在于，包括以下步骤：s1、建立水平井井筒流动模型：式中，t为时间；r为井筒内径；γ为井筒打开程度；ρ
l
为井筒流入流体密度；v
l
井筒中注入液速度；ρ
wb
为井筒中的流体密度；s2、设定井筒流动模型的初始和边界条件：v
i
＝v|
t＝0
＝0
ꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(2)s3、建立水平井井筒热模型，瞬态井筒热模型的最终形式为：式中，c
p
为比热容；θ为水平井筒倾角；t
r
表示储层温度；由储层热模型计算，t
l
表示注入液温度；g为重力加速度；u
t
为综合传热系数；t为温度；s4、设定水平井井筒热模型的初始和边界条件：t|
t＝0
＝t
wf
,
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(5)s5、对裂缝几何形态进行如下假设：(1)发生在裂缝处的流体滤失与裂缝面呈线性垂直；(2)注入过程中裂缝具有无限的导流能力；(3)压裂液不可压缩；并对裂缝扩展过程进行了建模，裂缝的流动模型为：式中，τ(x)是裂缝延伸长度x的时间之后，流体滤失开始；x
f
是裂缝半长，m；w、h分别为任意微元段的宽和高；c是考虑流体压缩性效应c
c
、流体粘度效应c
v
和造壁性效应c
w
的总流体滤失系数；s6、设定裂缝流动模型的初始和边界条件为：x
f
|
t＝0
＝0
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(8)
s7、建立裂缝热学模型，计算裂缝扩展和关井过程中裂缝内压裂液温度分布；能量守恒的概念与体积微分方程上的裂缝流动模型相似：s8、设定裂缝热学模型的初始和边界条件为:t
l
|
t＝0
＝t
i
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(11)式中，t
i
表示初始储层温度；t为注入流体温度；v
i
为x＝r和γ＝1时的流体速度；由井筒模型计算得出，在储层温度模型中，裂缝面边界与储层模型耦合；s9、压裂液从裂缝面向储层滤失的距离可通过将滤失速度与滤失时间相结合来计算，储层流动模型为：s10、建立了储层热学模型，并与井筒和裂缝模型相结合，根据能量平衡方程推导局部地层能量守恒方程，储层热学模型为：式(14)为局部地层能量守恒方程；对于被压裂液侵入的网格块，其平均有效热容和有效热导率不是恒定的；上述方程左端为随时间t在微元段中的能量积累：右端第一项是由于流体滤失进入地层而产生的热转换；第二和第三项是二维流动系统中的热传导；其中，k
e
是平均有效储层导热系数；t
r
表示地层岩石温度；s11、设定储层流动模型和储层热学模型的初始和边界条件为：t
r
|
t＝0
＝t
i
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(18)
s12、将上述步骤中的井筒、裂缝和储层模型耦合并结合边界条件进行求解，获取实时模拟井筒温度分布。2.根据权利要求1所述的水平井分段压裂温度分布预测方法，其特征在于，通过dts测量压裂注液过程中动态温度数据，从而进行裂缝诊断，实现裂缝位置识别、判断裂缝半长大小，并对压裂段内形成的裂缝数量和位置进行准确诊断。
技术总结本发明涉及一种水平井分段压裂温度分布预测方法，包括以下步骤：由质量、能量守恒建立水平井段井筒模型、裂缝模型和储层热学模型，对上述模型进行求解并进行模型耦合，运用真实数据来验证模型的输出结果，研究注入流量对水平井分段多簇压裂温度分布的影响，通过水平井压裂注液过程中的井筒温度剖面监测，对压裂段内形成的裂缝数量及进行准确诊断。本发明通过DTS测量压裂注液过程中动态温度数据，从而进行裂缝诊断，如裂缝位置识别、裂缝半长大小；通过本发明可对压裂段内形成的裂缝数量和位置进行准确诊断；本发明在考虑多种微量热效应的基础之上，实现了水平井分段压裂注液过程中的温度分布动态模拟预测。温度分布动态模拟预测。温度分布动态模拟预测。
技术研发人员：高云伟 张国锋 胡光 邹顺良 孙文常
受保护的技术使用者：中石化江汉石油工程有限公司 中石化江汉石油工程有限公司页岩气开采技术服务公司
技术研发日：2022.06.30
技术公布日：2022/11/1
