【技术实现步骤摘要】
本专利技术属于石油勘探开发领域，尤其是涉及一种中高孔渗疏松砂岩油藏速敏伤害计算方法。
技术介绍
1、储层速敏性是指在试油、采油、注水等作业过程中，当流体在储层中流动时，由于流体流动速度变化引起地层微粒运移，堵塞孔隙喉道，造成储层渗透率变化的现象。储层速敏影响因素较多，大体可分为内因和外因两类，内因包括：岩石矿物成分、孔隙结构、胶结程度等，外因包括：油水渗流速度、外来固相颗粒等。目前国内外对储层速敏性的研究手段主要为室内岩心实验，一般认为：低渗储层孔吼较小，流体速度越大，孔壁表面的矿物颗粒越容易被剥蚀，从而在孔吼处堆积，造成储层渗透率下降；中高渗储层孔吼较大，流体携带的矿物颗粒不易在孔道中堆积，渗透率变化不明显，甚至会出现流速越大渗透率越大的现象。
2、目前，虽然中高孔渗疏松砂岩油藏的注采井网相对完善，地层能量充足，但是定向井在高含水期开发过程中表现出产液量逐渐下降、井底流压逐渐下降的现象，相关测试显示表皮系数较大，储层伤害严重。经统计发现，定向井初期液量越高，后期液量递减趋势越明显，即初期渗流速度越大，后期储层伤害越大，速敏越严重。海上油田开展储层伤害评价的主要方法是压力恢复测试，但该方法费用较高，且需经常关停测试井，对油田开发过程中的成本控制和产量形势都造成了一定挑战。因此有必要在微观尺度的岩心实验基础上，开展宏观矿场尺度的储层速敏渗流机理研究，为中高孔渗疏松砂岩油藏后续高效开发提供指导。
技术实现思路
1、有鉴于此，本专利技术旨在提出一种中高孔渗疏松砂岩油藏速敏伤害
2、为达到上述目的，本专利技术的技术方案是这样实现的：一种中高孔渗疏松砂岩油藏速敏伤害计算方法，包括以下步骤：
3、s1，首先引入速敏伤害系数的概念，将渗流速度与渗透率关联；
4、s2，计算目前含水率下油水两相有效粘度；
5、s3，计算地层导压系数；
6、s4，建立定液量条件下单相流不稳定渗流公式；
7、s5，建立定流压条件下单相流不稳定渗流公式；
8、s6，采用中心显示格式离散化压力控制方程；
9、s7，将目标油田实际参数代入上述所有公式中，与无速敏条件进行对比，计算速敏伤害导致的附加表皮。
10、进一步的，步骤s1中，速敏伤害系数的取值范围为0-1，渗流速度与渗透率关系式如式(1)所示，
11、
12、式中，k为渗透率，μm2；ki为原始渗透率，μm2；αk为速敏伤害系数，无量纲；vw为井壁处渗流速度，m/d；v为任意半径处的渗流速度，m/d，计算公式如下：
13、
14、式中，μr为地层流体有效粘度，mpa·s；p为地层压力；r为渗流半径，m；q为日产液量，m3/d；h为油层厚度，m。
15、进一步的，步骤s2中，地层流体有效粘度的计算如式(3)所示，
16、μr＝1/(kro/μo+krw/μw) (3)
17、式中，kro为油相相对渗透率，无量纲；krw为水相相对渗透率，无量纲；μo为油相粘度，mpa.s；μw为水相粘度，mpa.s；μr为地层流体有效粘度，mpa·s。
18、进一步的，油相相对渗透率kro和水相相对渗透率krw分别根据相渗曲线获得。
19、进一步的，步骤s3中，地层导压系数的计算如式(4)所示，
20、
21、式中，η为地层导压系数，m2/d；k为渗透率，μm2；μr为地层流体有效粘度，mpa·s；ct为综合压缩系数，mpa-1；综合压缩系数ct根据实验参数获取。
22、进一步的，定液量条件下单相流不稳定渗流公式如式(5)所示，
23、
24、定液量条件单相流不稳定渗流公式如式(6)所示，
25、
26、式中，p为地层压力，mpa；r为渗流半径，m；rw为井筒半径，m；μr为地层流体有效粘度，mpa·s；pi为地层压力，mpa；t为生产时间，d；q为产液量，m3/d；k为渗透率，μm2；h为储层厚度，m；pwf为井底流压，mpa。
27、进一步的，所述步骤s4中，产液量q为稳定生产阶段的平均值，压力导数取井筒半径处的值；步骤s5中，井底流压pwf为泵吸入口压力变化较缓慢阶段时的平均值。
28、进一步的，步骤s6具体为，将式(5)和式(6)中的控制方程进行离散化，离散方程采用中心显示格式，如下式所示：
29、
30、式中，为下一时间步当前空间步地层压力，mpa；为当前时间步及下一空间步地层压力，mpa；为当前时间步及当前空间步地层压力，mpa；为当前时间步及前一空间步地层压力，mpa；为当前时间步地层导压系数，m2/d；δr为空间步长，m；δt为时间步长，d；r为渗流半径，m；
31、空间步长δr不超过井筒半径，时间步长δt不超过10分钟。
32、进一步的，步骤s7中，计算速敏伤害导致的附加表皮，指的是将有速敏条件下的比采液指数与无速敏条件下的比采液指数进行对比，即可得到疏松砂岩油藏在速敏条件下的储层伤害程度，计算公式如下式所示，
33、
34、式中，s为表皮系数，无量纲；jl为有速敏情况下的比采液指数，m3/d/mpa/m；jl’为无速敏情况下的比采液指数，m3/d/mpa/m；re为定向井泄油半径，m；rw为井筒半径，m。
35、进一步的，将计算的速敏伤害理论值与压力恢复测试实际值进行对比，判断符合程度。
36、相对于现有技术，本专利技术所述的中高孔渗疏松砂岩油藏速敏伤害计算方法具有以下优势：能够在不关停油井的情况下，仅利用实际生产动态参数，就可以评估任意时刻的储层伤害程度，同时节约油田测试成本；本专利技术所述方法的计算结果与现场实际压力恢复测试结果较为接近，证明了该方法的可靠性，可应用于海上类似油田储层伤害评价。
37、(1)本专利技术所述方法可实现油井生产全过程模拟；
38、(2)本专利技术所述方法井底流压计算结果与实际井底流压拟合程度较高；
39、(3)本专利技术所述方法计算表皮与实际测试表皮接近，误差在15％以内；
本文档来自技高网...【技术保护点】
1.一种中高孔渗疏松砂岩油藏速敏伤害计算方法，其特征在于，包括以下步骤：
2.根据权利要求1所述的中高孔渗疏松砂岩油藏速敏伤害计算方法，其特征在于：步骤S1中，速敏伤害系数的取值范围为0-1，渗流速度与渗透率关系式如式(1)所示，
3.根据权利要求2所述的中高孔渗疏松砂岩油藏速敏伤害计算方法，其特征在于：步骤S2中，地层流体有效粘度的计算如式(3)所示，
4.根据权利要求3所述的中高孔渗疏松砂岩油藏速敏伤害计算方法，其特征在于：油相相对渗透率kro和水相相对渗透率krw分别根据相渗曲线获得。
5.根据权利要求3所述的中高孔渗疏松砂岩油藏速敏伤害计算方法，其特征在于：步骤S3中，地层导压系数的计算如式(4)所示，
6.根据权利要求4或5所述的中高孔渗疏松砂岩油藏速敏伤害计算方法，其特征在于：定液量条件下单相流不稳定渗流公式如式(5)所示，
7.根据权利要求6所述的中高孔渗疏松砂岩油藏速敏伤害计算方法，其特征在于：所述步骤S4中，产液量Q为稳定生产阶段的平均值，压力导数取井筒半径处的值；步骤S5中，井底流压Pwf
8.根据权利要求7所述的中高孔渗疏松砂岩油藏速敏伤害计算方法，其特征在于：步骤S6具体为，将式(5)和式(6)中的控制方程进行离散化，离散方程采用中心显示格式，如下式所示：
9.根据权利要求8所述的中高孔渗疏松砂岩油藏速敏伤害计算方法，其特征在于：步骤S7中，计算速敏伤害导致的附加表皮，指的是将有速敏条件下的比采液指数与无速敏条件下的比采液指数进行对比，即可得到疏松砂岩油藏在速敏条件下的储层伤害程度，计算公式如下式所示，
10.根据权利要求9所述的中高孔渗疏松砂岩油藏速敏伤害计算方法，其特征在于：将计算的速敏伤害理论值与压力恢复测试实际值进行对比，判断符合程度。
...【技术特征摘要】
1.一种中高孔渗疏松砂岩油藏速敏伤害计算方法，其特征在于，包括以下步骤：
2.根据权利要求1所述的中高孔渗疏松砂岩油藏速敏伤害计算方法，其特征在于：步骤s1中，速敏伤害系数的取值范围为0-1，渗流速度与渗透率关系式如式(1)所示，
3.根据权利要求2所述的中高孔渗疏松砂岩油藏速敏伤害计算方法，其特征在于：步骤s2中，地层流体有效粘度的计算如式(3)所示，
4.根据权利要求3所述的中高孔渗疏松砂岩油藏速敏伤害计算方法，其特征在于：油相相对渗透率kro和水相相对渗透率krw分别根据相渗曲线获得。
5.根据权利要求3所述的中高孔渗疏松砂岩油藏速敏伤害计算方法，其特征在于：步骤s3中，地层导压系数的计算如式(4)所示，
6.根据权利要求4或5所述的中高孔渗疏松砂岩油藏速敏伤害计算方法，其特征在于：定液量条件下单相流不稳定渗流公式如式(5)所示，
...
【专利技术属性】
技术研发人员：于登飞，王月杰，郑举，李彪，张俊，黄琴，李景玲，龙明，白冰，张博文，
申请(专利权)人：中海石油中国有限公司天津分公司，
类型：发明
国别省市：
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。