Discussion on a Formation Pore Pressure Detection Method for Carbonate Rocks Based on the Thin Plate Theory
-
摘要: 现有碳酸盐岩地层压力检测方法均存在不足,因此,为保证Y油田F地层碳酸盐岩地层的钻井安全,开展了碳酸盐岩地层压力的检测方法研究。基于薄板理论,考虑体积弹性模量的影响,结合碳酸盐岩特征,建立了构造挤压条件下的地层压力地质力学识别模型;通过分析F地层碳酸盐岩的地质构造、异常高压及测井响应特征分布规律,结合排除法研究了异常高压的成因机制;形成的碳酸盐岩地层压力检测方法在Y油田F层进行了实例应用。F层小断层发育,纵波速度在5 500 m/s左右;FU层和FL层的地层压力系数分别约为1.45和1.30;FU层岩石骨架变形量大于FL层,构造挤压是异常高压产生的主要成因机制。应用结果表明,该模型检测值与SFT实测值间的相对误差小于10%;地层压力随构造变形曲率、地层压力系数、弹性模量的增大而增大,随泊松比的增大而减小,且呈线性关系。研究认为,基于薄板理论的碳酸盐岩地层压力检测方法,能够比较准确地检测由构造挤压作用下的碳酸盐岩地层压力。Abstract: The current formation pore pressure detection methods for carbonate rock have shortcomings.In order to ensure drilling safety of carbonate rock in the F formation,Y Oilfield,the formation pore pressure detection method was studied.The geo-mechanical model of formation pore pressure was proposed under tectonic compression,based on the thin plate theory,considering the effect of bulk modulus and in accordance with carbonate rock characteristics.Based on the analysis of the geological structure,abnormal over-pressure and logging response of carbonate rock in F formation,combined with exclusive method,the generation mechanism of abnormal over-pressure was studied.When a formation pore pressure detection method was applied to F formation,Y Oilfield,the study showed that the minor faults were developed and the P-wave velocity is about 5 500 m/s in F formation.Formation pore pressure coefficient of FU formation and FL formation was about 1.45 and 1.30 respectively.The deformation of rock framework of FU formation was larger than that of FL formation,and the main generation mechanism of abnormal over-pressure was based on tectonic compression.The application results showed that the error of formation pore pressure detected value was less than 10% when compared with SFT measured value.Formation pore pressure increased with the increase of tectonic deformation curvature,formation pore pressure coefficient and elasticity modulus,and decreased with the increase of Poisson’s ratio,which was shown in the linear relationship.The study suggests that formation pore pressure detection method of carbonate rock based on the thin plate theory could be used to detect pore pressure of carbonate rock caused by tectonic compression.
-
分层注水是通过控制不同注水层的水嘴实施注水,以缓解层间矛盾,提高注水波及系数的一项技术措施[1-6]。目前,国内分层注水技术经过 60多年的发展,形成了以固定式、钢丝投捞式和电缆测调式为代表的分层注水技术,在支撑油田持续高产稳产、提高水驱采收率等方面发挥了重要作用[7]。
油田开发中后期,由于油藏储层的非均质性,油田开发效果较差,而降低高渗透率注水层的注水量、增加低渗透率注水层的注水量,可获得较佳的注水开发效果。因此,分层注水量的计算和调配是分层注水的关键[8]。目前,计算分层注水量的方法主要有统计学方法、数值模拟方法和油藏工程方法。统计学方法主要研究各个影响因素与注水量之间的关系,对样本的数量及质量要求较高,分层注水测试数据量较少的区块难以应用。数值模拟方法主要是通过建立分层注水数值模型,模拟计算不同注水量下油田开发效果及评价剩余油分布规律,确定分层注水的合理配注量[9-10],数值模拟方法考虑的因素较多,计算复杂,在现场应用难度大。油藏工程方法主要有注采比法、水驱前缘法、有效厚度法、吸水剖面系数法、静态地层系数法和动态综合多因素劈分系数法等 [11-13],这些方法计算较简单,但精度较低。张俊廷等人[14]采用灰色关联分析方法,研究了影响注水量劈分系数的主要因素及其在劈分系数中所占权重,建立了计算分层注水量的公式。郭哲源和李国庆等人[15-16]建立了基于BP神经网络的分层注水量预测模型。赵广渊等人[17]建立了考虑管流、嘴流和地层流动的注水量计算方法,但只是简单地采用有效注入压力和地层静压的压差表述地层流动。董驰等人[18]建立了综合考虑地质参数和原油黏度动态变化的渗流阻力计算模型,明确了特高含水期分层注水界限。
目前,已有计算分层注水量的油藏工程方法对于地层渗流考虑得较为简单,并且未考虑注入过程中压力及注入流量随时间变化的特征。同时,现场注水过程中,还需要考虑井口泵注最大压力及地层破裂压力等约束条件[19-20],以保证注水安全,但现有计算方法均未涉及该类约束条件。因此,笔者建立了考虑井筒管流–水嘴孔流–地层渗流的耦合模型,并考虑井口压力及地层压力的约束,构建了一种水嘴直径反演计算方法,通过模型正演及反演计算,分析了模型的计算精度、验证了模型计算结果的可靠性,为分层注水提供了一种新的、更准确的注水量计算及水嘴调配方法。
1. 分层注水模型的构建及求解
分层注水过程中,注入水从井口注入地层,主要包含了井筒管流、水嘴孔流和地层渗流3个流动阶段,而压力是注入水在各流动阶段的驱动力。首先分别建立井筒管流、水嘴孔流和地层渗流模型,然后考虑连接点处压力流量的耦合条件,最后建立分层注水井管流–嘴流–渗流耦合模型。
1.1 数学模型的建立
为方便数学模型的建立,作出以下假设:1)油藏有多个注水层,每层注水层厚度均匀,但厚度、渗透率、孔隙度等油藏物性参数可以不同,各层之间在地层中互不连通;2)各注水层通过配水器水嘴与井筒连通;3)油藏中只有水流动,其流动符合达西定律,油藏内不考虑重力及毛细管力;4)注入水在井筒中的流动视为管流;5)考虑表皮因子影响,不考虑井筒储集系数。分层注水物理模型如图1所示。
1)井筒管流模型。注入水在井筒的流动主要包括从井口到配水器水嘴前的流动,根据动量守恒定理可得注入水在井筒中的动量守恒方程。
dpdz=ρgsinθ−ρvdvdz−fρv22d (1) 式中:
p 为压力,Pa;z 为垂直深度,m;ρ 为流体密度,kg/m3;v 为注入水在井筒中的流动速度,m/s;f 为摩阻系数;d 为井筒直径,m;θ 为井筒倾斜角,(°)。由于分层注水的原因,摩阻需要分段计算,从井口到第1注水层的流量为总注入量,第1注水层到第2注水层间的流量为总流量减去第1注水层的注入量,以此类推,进行不同段的摩阻计算及井筒内水嘴前压力的计算。
2)水嘴孔流模型。由流体力学基本原理可知,流体经过标准节流装置(孔板、喷嘴、文丘里管等)时,节流前后的静压力差Δp(嘴损)与体积流量q的关系为:
Δp=8ρζπ 2dw2q2 (2) 式中:Δp为水嘴前后的压差,Pa;
q 为体积流量,m3/s;dw 为水嘴直径,m;ζ 为局部阻力系数。注水过程中,不同注水层配水器水嘴的直径不同,导致压差不同以及流量不同,因此,分层注水量的调整通常是通过调整水嘴来实现。
3)油层渗流模型。油藏内的渗流模型与储层类型等有关,笔者对常见的直井均质无限大模型进行分析,依据储层流体质量守恒方程及达西渗流定律得到储层渗流模型:
∇2p=ϕμCtK∂p∂t (3) 式中:
∇ 为哈密顿算子;ϕ 为孔隙度;μ为黏度,Pa·s;Ct为综合压缩系数,Pa−1;K为渗透率,m2;t 为时间,s。根据初始条件和边界条件,可以求解模型得到井底压力,模型求解方法较成熟,此处不再赘述。实际应用过程中,可以根据实际储层情况选取不同的渗流模型进行求解,如双重孔隙介质模型、径向复合模型及压裂井模型等;同样,不同注水层也可以采用不同的模型进行求解,主要目的是得到井壁处的压力。
基于上述3个独立的模型,考虑井筒管流计算得到的注水层处井筒压力与水嘴前的压力相等和水嘴后的压力与渗流模型计算出的井壁处压力相等2个连续条件,将3个模型耦合形成分层注水井的管流–嘴流–渗流模型。
1.2 模型求解方法
数学模型中,油管管流、水嘴孔流及地层渗流构成了油管–水嘴–地层的流动系统。采用全隐式方法耦合求解法很难求解该数学模型,其原因是油管、水嘴与地层流动系统的求解方程形式存在很大差异。同时,油管、水嘴与储层网格的体积相差3~4个数量级,求解的雅可比矩阵很难收敛,从而导致效率较低。因此,考虑对3个模型分别进行求解,再通过判断是否满足连续条件来确定求解结果是否正确,然后采用迭代的方式求得最终解。
模型求解过程中,需要注意的是,由于上文所构建的模型考虑到注水过程中井底压力是随时间持续变化的,导致分层注水量、井口压力也随时间持续变化。因此,在给定分层配注量的条件下,并不能通过给定一个固定的分层注水量,利用所建模型直接计算出配水器水嘴直径,而配水器水嘴直径通常是固定的。故模型的正演考虑从给定配水器水嘴直径出发,计算不同时刻的分层注水量、井底压力和井口压力。模型正演的计算流程如图2所示。
模型正演计算流程包含了2个循环和2个迭代过程:循环1为时间上的循环,计算不同时刻的分层注水量及压力,需要注意的是,应用渗流模型计算井壁处的压力时,由于不同时刻的流量不同,需采用叠加原理计算井壁处的压力;循环2为层段的循环,需要循环计算所有层段的流量和压力;迭代1为各层段注入量的迭代,主要是结合嘴流模型及渗流模型以及对应的连续条件迭代计算各层段的注水量;迭代2为井口压力的迭代,要满足所构建的各模型方程,井口压力也需满足对应的连续条件,通过判断各层段注入量之和是否为总注入量来迭代计算所需井口压力。
2. 考虑压力约束的反演模型
由于分层注入量随时间变化的特征,所构建的模型主要通过配水器水嘴直径计算分层注水量;而现场应用过程中,主要是通过配注量确定所需水嘴直径。因此,考虑采用最优化理论,确定各注水层配注量对应的最优水嘴直径。
注入量和水嘴直径的关系为:
q=f(d) (4) 式中:
q=(q1,q2,⋯,qn) 为所给注水时间最后时刻所有层的注入量,m3/d;d=(d1,d2,⋯,dn) 为所有层的水嘴直径,m。因此,最优化目标函数为:
g(d)=‖ (5) 式中:
Q = ({Q_1},{Q_2}, \cdots ,{Q_n}) 为所有层的配注量,m3/d。计算过程中,结合现场实际注水情况及储层特性,还需考虑2个约束条件:一是井口压力的约束,一般来说,现场通过柱塞泵将水注入井内,井口注入压力会存在一个最大值,因此在计算过程中计算出的井口压力需满足该约束条件;二是油藏压力的约束,注水过程中油层压力不能大于油层的破裂压力,通常还会给一个0.9的安全系数,因此,通过渗流模型计算出的井壁处压力需满足该约束条件。
由于井口压力及井底压力均属于模型计算过程中的中间参数,考虑采用最优化方法中的惩罚函数法引入压力约束条件,定义函数为:
\begin{split} F(d) =& g(d) + \sigma \sum\limits_{i = 1}^n {\max {{\left[ {0,({p_{{\text{wf}}i}} - 0.9{p_{{\text{wf}}\max }})} \right]}^2}} + \\ &\sigma \sum\limits_{j = 1}^m {\max {{\left[ {0,({p_{{\text{c}}j}} - {p_{{\text{c}}\max }})} \right]}^2}} \end{split} (6) 式中:
{p_{{\text{wf}}i}} 为各层井壁处压力,MPa;{p_{{\text{wf}}\max }} 为地层破裂压力,MPa;{p_{{\text{c}}j}} 为不同时刻井口压力,MPa;{p_{{\text{c}}\max }} 为最大井口压力,MPa;n为总注水层数量;m为总时间段数量;σ为一个较大的正数。当计算出的井底压力大于0.9倍破裂压力、计算出的井口压力大于最大井口压力时,式(6)右侧后2项为一个较大的数,无法满足最优化条件,而当计算出的井底压力小于0.9倍破裂压力、计算出的井口压力小于最大井口压力时,右侧后2项为0,此时考虑压力约束的最优化问题转化为:
\min \;\;\;F(d) (7) 通过给定目标配注量及约束条件,即可反演计算得到各注水层最优的水嘴直径,指导现场分层注水工作。
3. 计算结果分析
3.1 模型正演计算结果分析
给定一组参数,利用上文构建的模型计算不同注水时间下井口压力、各注水层压力及各注水层的注入量,计算结果如图3和图4所示。给定的参数:总注水量100 m3/d;注水时间20 d;井筒半径0.076 2 m;注入水的黏度0.6 mPa·s、体积系数1m3/m3、密度1 kg/L;分3层注水,各层孔隙度与有效厚度相同,渗透率不同,各注水层的参数见表1。
表 1 注水层的参数Table 1. Parameters of each water injection layer注水层序号 有效厚度/m 储层深度/m 渗透率/mD 孔隙度 水嘴直径/mm 第1层 15 2 000 2 0.12 2 第2层 15 2 050 3 0.12 3 第3层 15 2 110 5 0.12 5 由图3可知,随着注水时间增长,井口压力及各注水层井壁处压力均呈逐渐升高趋势,井口压力相对较低,各注水层井壁处压力相对较高,各注水层之间压差较小,第1注水层井壁处压力最低,第3注水层井壁处压力最高。由图4可知,随着注水时间增长,各注水层的注入流量也呈逐渐变化的趋势,其中第1、第2注水层的注入量随注水时间增长逐渐增大,第3注水层的注入量逐渐减小,各注水层的注入量总体均逐渐趋于平稳,最后时刻第1、第2和第3注水层的注入量分别为16.3,28.8和54.9 m3/d。
在此基础上,通过给定各注水层配水器水嘴的直径,研究其对各注水层注入量的影响。各注水层配水器水嘴直径的变化,笔者设计了4种变化方式,前3种变化方式是2个注水层配水器的水嘴直径不变,只是1个注水层配水器的水嘴直径变化,最后1种变化方式是3个注水层配水器的水嘴直径同时变化,具体变化情况见表2。
表 2 各注水层配水器水嘴直径的变化方式Table 2. Change modes of nozzle diameter of water distributors in each water injection layer变化方式 第1层水嘴直径/mm 第2层水嘴直径/mm 第3层水嘴直径/mm 方式1 2,3,4,5,6,7,8 3 5 方式2 2 2,3,4,5,6,7,8 5 方式3 2 3 2,3,4,5,6,7,8 方式4 2,3,4,5,6,7,8 3,4,5,6,7,8,9 5,6,7,8,9,10,11 利用上文所建模型计算各注水层配水器水嘴直径不同变化方式下,各注水层最后时刻的注入量,结果如图5—图8所示。
由图5—图7可知,当只改变某一注水层配水器水嘴直径时,随着该层配水器水嘴直径逐渐增大,该层注入量逐渐增大,但同时其他2层注入量逐渐减小。由图8可知,当配水器水嘴直径均逐渐增大时,配水器水嘴直径较小注水层的注入量逐渐增大,而配水器水嘴直径较大注水层的注入量逐渐减小,此时,由于层与层之间相互影响,各注水层的注入量与与各层配水器水嘴直径不再是简单的正相关关系。
3.2 模型反演计算结果分析
在正演的基础上,给定一组分层配注量进行反演计算,考虑增大低渗层配注量,给定3个注水层的目标配注量分别为30,20和50 m3/d,基于3.1节的反演计算结果,首先不考虑压力约束条件进行反演计算,结果如图9和图10所示。第1、第2和第3层注水层配水器的水嘴直径分别为1.7,1.2和2.0 mm。
从图10可以看出,3个注水层的反演注入量逐渐趋于配注量,相对误差小于0.1%,计算精度高,模型的拟合效果好。从图9可以看出,由于第1注水层的配注量较大,导致该层井壁处压力最大,井口压力也较高。考虑压力约束条件时,油层破裂压力为37.8 MPa,最大井口压力为20.0 MPa,在此基础上反演计算井口压力和各注水层井壁处的压力及注入量,结果如图11和图12所示。第1、第2和第3注水层配水器的水嘴直径分别为1.7,1.2和2.0 mm。
从图11可以看出,井口压力和各注水层井壁处的压力逐渐升高,但均满足压力约束条件,各注水层井壁处压力最高为34 MPa,低于0.9倍的破裂压力,井口压力最高为19.8 MPa,低于最大井口压力。从图12可以看出,第1、第2和第3注水层最后时刻的反演注入量分别为28.9,20.2和50.9 m3/d,反演出的各层注入量与各层的配注量存在一定的误差,分析其原因可能是压力约束造成的。因此,在实际计算过程中,应给定一个现场误差,分析其误差是否满足要求,若误差满足要求,则说明计算结果可行;若误差超出要求,则说明在所给定压力约束条件下,无法达到所需配注量要求,建议修改配注量或采取压裂酸化等措施对油层进行改造,增大其吸水能力。
4. 结 论
1)考虑井筒管流–水嘴孔流–地层渗流,建立了分层注水量计算模型,并考虑压力约束条件,构建了各注水层配水器水嘴直径反演计算模型,以指导现场分层配注及水嘴调配。
2)正演计算结果表明,模型能够精确计算各注水层随时间变化的注入量及井口压力、井底压力,并且其与配水器水嘴直径有关,单层配水器水嘴直径与该层注入量呈正相关关系,与其他层注入量呈反相关关系,当各注水层配水器水嘴直径均发生变化时,注入量变化规律与各层储层物性等相关。
3)反演计算结果表明,模型计算精度较高,压力约束条件能够较好地约束计算过程中的井口压力、井底压力,保证现场注水安全,同时计算出的各层注入量与配注量存在误差,需根据误差情况提出下一步分层注水调配建议。
-
[1] 马海.Fillippone地层压力预测方法的改进及应用[J].石油钻探技术,2012,40(6):56-61. Ma Hai.Pore pressure prediction with seismic interval velocity by the modified Fillippone method[J].Petroleum Drilling Techniques,2012,40(6):56-61. [2] Eaton B A.The equation for geopressure prediction from well logs[R].SPE 5544,1975.
[3] Han D H,Nur A,Morgan D.Effects of porosity and clay content on wave velocities in sandstone[J].Geophysics,1986,51(11):2093-2107.
[4] Li Qiuguo,Heliot Denis,Zhao Liangxiao,et al.Abnormal pressure detection and wellbore stability evaluation in carbonate formations of east Sichuan,China[R].IADC/SPE Drilling Conference,New Orleans,America,2000.
[5] 樊洪海.利用声速检测欠压实泥岩异常高压的简易方法与应用[J].石油钻探技术,2001, 29(5):9-11. Fan Honghai.A simple pore pressure estimation method for a disequilibrium compaction shale using sonic velocity[J].Petroleum Drilling Techniques,2001,29(5):9-11. [6] Ye Zhi,Fan Honghai,Liu Gang,et al.Estimating formation pore pressure in tectonic compression zones[J].Petroleum Science and Technology,2012,30(8):766-774.
[7] Atashbari V,Tingay M R.Pore pressure prediction in a carbonate reservoir[R].SPE 150836,2012.
[8] Bowers G L.Pore pressure estimation from velocity data:accounting for overpressure mechanisms besides undercompaction[R].IADC/SPE 27488,1994.
[9] Jin Yan,Chen Mian,Chen Kangping,et al.Pre-caspian basin wells in salt-gypsum beds require an optimized drilling fluid[J].Oil Gas Journal,2010,108(4):46-53.
[10] 李茂,何俊才,李玮.地层压力的多井对比综合预测技术[J].特种油气藏,2008,15(1):88-91. Li Mao,He Juncai,Li Wei.Comprehensive prediction of formation pressure by multiwell correlation[J].Special Oil and Gas Reservoirs,2008,15(1):88-91. [11] Weakley P R.Determination of formation pore pressures in carbonate environments from sonic logs[R].SPE 9009,1990.
[12] 杨进,刘书杰,石磊,等.挤压构造地层压力预测模型研究[J].石油学报,2009,30(5):764-768. Yang Jin,Liu Shujie,Shi Lei,et al.Research on prediction model for formation pressure in compression structure[J].Acta Petrolei Sinica,2009,30(5):764-768. [13] Timoshenko Sp,Krieger W.Theory of plates and shells[M] .2nd ed.New York:McGraw-Hill,1959.
[14] Li Shulin,George Jeff,Pardy Cary.Pore-Pressure and wellbore-stability prediction to increase drilling efficiency[J].Journal of Petroleum Technology,2012,64(2):99-101.
[15] 陈勉,金衍,张广清.石油工程岩石力学[M].北京:科学出版社,2008:34-40. Chen Mian,Jin Yan,Zhang Guangqing.Petroleum related rock mechanics[M].Beijing:Science Press,2008:34-40. [16] He Xilei,He Zhenhua,Wang Ruiliang,et al.Calculations of rock matrix modulus based on a linear regression relation[J].Applied Geophysics,2011,8(3):155-162.
[17] Biot M A.Generalized theory of acoustic propagation in porous dissipative media[J].The Journal of the Acoustic al Society of America,1962,34(5):1254-1262.
[18] Pride S R,Berryman J G,Harris J M.Seismic attenuation due to wave-induced flow[J].Journal of Geophysical Research:Solid Earth,2004,109(B1):1201-1221.
-
期刊类型引用(11)
1. 成海,张逸群,王印,刘超. 超深定向井地质工程力学耦合机理研究进展与认识. 油气藏评价与开发. 2024(04): 593-599+646 . 百度学术
2. 余金秋,孙夕平,于永才,李璇,李晓明,徐光成. 深层地层压力预测技术研究进展. 海相油气地质. 2024(04): 337-347 . 百度学术
3. 伍贤柱,胡旭光,韩烈祥,罗园,许期聪,庞平,李黔. 井控技术研究进展与展望. 天然气工业. 2022(02): 133-142 . 百度学术
4. 路保平,王志战,张元春. 碳酸盐岩孔隙压力预监测理论与方法进展. 石油学报. 2022(04): 571-580 . 百度学术
5. 邓虎,贾利春. 四川盆地深井超深井钻井关键技术与展望. 天然气工业. 2022(12): 82-94 . 百度学术
6. 蔡文军,邓金根,张建峰,王大勇,孙源秀,冯永存. 印尼A油田礁灰岩储层的三维孔隙压力预测及可视化. 科学技术与工程. 2021(21): 8863-8869 . 百度学术
7. 时梦璇,刘之的,杨学峰,杨珺茹,陈秀娟,刘洪壮,曹晋璐. 地层孔隙压力地球物理测井预测技术综述及展望. 地球物理学进展. 2020(05): 1845-1853 . 百度学术
8. 熊晓军,黄劲,陈容,LIAO Yiduo,袁野. 一种适用于欠压实成因的地层孔隙流体压力预测技术. 物探化探计算技术. 2019(01): 1-5 . 百度学术
9. 王明飞,苏克露,肖伟,陈超,刘晓晶,王静波. 泥页岩应力场约束的叠后地震裂缝预测技术——以焦石坝区块五峰组-龙马溪组一段为例. 石油与天然气地质. 2018(01): 198-206 . 百度学术
10. 路保平,鲍洪志,余夫. 基于流体声速的碳酸盐岩地层孔隙压力求取方法. 石油钻探技术. 2017(03): 1-7 . 本站查看
11. 豆宁辉,赵向阳,余夫,金衍,潘堤. 碳酸盐岩纵波速度响应特征试验研究. 科学技术与工程. 2016(16): 9-12+34 . 百度学术
其他类型引用(4)
计量
- 文章访问数: 4001
- HTML全文浏览量: 90
- PDF下载量: 3710
- 被引次数: 15