1、录 井 工 程第34卷第2期基于六臂地层倾角资料的井眼扩大系数对录井迟到时间的校正徐长凯刘晓亮王 俊(中国石油大庆钻探工程公司地质录井一公司)徐长凯工程师,1970年生,2014年毕业于东北石油大学石油工程专业,现在大庆钻探地质录井一公司资料采集二分公司从事综合录井管理工作。通信地址:163411 黑龙江省大庆市地质录井一公司资料采集二分公司。电话:(0459)5691861。Email:0引言伊拉克B油田位于伊拉克南部与伊朗交界处,构造上位于中东波斯湾盆地扎格罗斯山构造带西缘,为轴向近似南北向的长轴背斜1。该区域曾经历早白垩世被动大陆边缘与晚白垩世-新生代前陆盆地多期较强烈的南北向差异性构造
2、活动与近东西向的挤压2。受地质构造应力的剪切、拉伸作用以及钻头附近地应力释放等多重力学作用影响,Sulaiy地层“糖葫芦串”形的椭圆井眼非常发育,其椭圆长轴(最小水平主应力)发育方位几乎一致,集中在144324,椭圆短轴(最大水平主应力)方位集中在54234,井眼受地质构造应力作用明显。在录井资料采集过程中,鉴于椭圆井眼的发育程度对录井迟到时间的影响较明显,本文基于六臂地层倾角测井资料,采用拟合最佳逼近椭圆求解井眼扩大系数的方法,对迟到时间理论摘要伊拉克B油田的Sulaiy地层“糖葫芦串”形的椭圆井眼非常发育,对迟到时间影响较大。基于六臂地层倾角测井资料,采用最佳逼近的椭圆拟合井眼实际情况,并
3、引入井眼扩大系数克服井径扩大率表征裸眼段井眼规则度方面的不足,解决了井眼不规则而引起的迟到时间不精准问题。在实际生产中,应用文中迟到时间校正公式对伊拉克B油田5口井迟到时间进行校正,使录井深度(迟到井深)与电测深度绝对误差降至1.4 m以内,满足了B油田生产需求,效果良好。关键词迟到时间Sulaiy地层椭圆井眼六臂地层倾角最佳逼近椭圆拟合井径扩大率井眼扩大系数中图分类号:TE 132.1文献标识码:ADOI:10.3969/j.issn.16729803.2023.02.007Correction of mud logging lag time by borehole enlargement
4、coefficient based on sixarmdipmeter logXU Changkai,LIU Xiaoliang,WANG JunNo.1 Mud Logging Company of CNPC Daqing Drilling&Exploration Engineering Corporation,Daqing,Heilongjiang163411,ChinaAbstract:The sugarcoated haws on sticks string shaped elliptical hole of Sulaiy Formation in Iraq B oilfiled is
5、 verydeveloped and has a great impact on the lag time.Based on the sixarm dipmeter log,the best approximation ellipse is used tofit the actual situation of the borehole.And the introduction of borehole enlargement coefficient overcomes the shortcomings ofhole enlargement ratio in characterizing the
6、borehole regularity of the open hole interval,and solves the problem of inaccuratelag time caused by borehole irregularity.In actual production,the lag time correction formula in this paper was applied tocorrect the lag time of five wells in Iraq B oilfield,which can reduce the absolute error betwee
7、n mud logging depth(lag welldepth)and electric depth to within 1.4 m,and meet the production demand of B oilfield with good results.Key words:lag time,Sulaiy Formation,elliptical hole,sixarm stratigraphic dip,best approximation,ellipse fitting,holeenlargement ratio,borehole enlargement coefficient引用
8、:徐长凯,刘晓亮,王俊.基于六臂地层倾角资料的井眼扩大系数对迟到时间的校正J.录井工程,2023,34(2):3944.XU Changkai,LIU Xiaoliang,WANG Jun.Correction of mud logging lag time by borehole enlargement coefficient based onsixarm dipmeter logJ.Mud Logging Engineering,2023,34(2):3944.39录 井 工 程工 艺 技 术2023年6月公式进行修正,使之更符合实际情况,提高了录井资料采集精准度。1录井迟到时间理论计算公
9、式录井迟到时间是指在石油天然气钻井过程中,钻井液携带岩屑与油气水,沿着井筒与钻具所形成的环形空间从井底上返到井口取样位置所需的时间,其精准度直接影响上返岩屑、油气水、钻井液及气测等多项资料的采集质量,特别是在快速钻井条件下,迟到时间不准确容易导致对储层及储层流体性质的误判3。迟到时间理论计算公式如下:T迟V1+V2Q()d12-d22h1+()D2-d22h24Q(1)其中:V114()d12-d22h1V214()D2-d22h2式中:T迟为迟到时间理论值,min;V1为套管段环空体积,m3;V2为裸眼段环空体积,m3;Q为钻井液泵总排量,m3/min;D为钻头直径,m;d1为套管内径,m;
10、d2为钻具外径,m;h1为套管下深,m;h2为钻头在裸眼井段下深,m。由录井迟到时间的理论计算公式可知:在钻井液泵总排量、钻头直径、钻具外径、钻头下深等参数一定时,裸眼段环空体积是影响迟到时间精准度的最关键因素。在井筒形成过程中,钻头周围应力场的改变会导致井壁失稳,井眼通常呈现不规则形态(图1),而现场迟到时间的计算一般采用钻头尺寸代替井径,忽略了裸眼段井眼空间规则与否的影响,使得迟到时间理论计算值与实际情况误差较大,且随着井深的增加,误差会越来越大4。因此,利用井眼扩大系数量化裸眼井段井眼规则程度,并对理论迟到时间进行校正,可提高资料采集的精准度。2录井迟到时间的校正方法2.1六臂地层倾角测
11、井原理简述研究表明57:力学原因引起的井眼坍塌最终形成椭圆井眼,且椭圆面垂直于井眼轴线8,其长轴远大于钻头直径,短轴则接近钻头直径或者比钻头直径大的多9。伊拉克B油田作业现场采用六臂地层倾角测井获取椭圆井眼的数据,所用测井仪器具有6个彼此对称、互成60夹角的推靠臂10,确保6个极板与井壁有可靠的接触,以便获取更多的测井信息,具有判定井眼几何形状的能力11,对井眼的几何形状反映更为灵敏12。6条极板测量的井径曲线分别为C1、C2、C3、C4、C5、C6,3条相隔60的双井径曲线分别为C14、C25、C36(图2)。测井仪器在测量过程中会随着提升而不断旋转,当测遇崩落井段时,一对测臂将嵌入椭圆井眼
12、长轴方向的槽内,使仪器不再随着提升而旋转13,从而获得椭圆形井眼长轴所在的方位。当C14为最大时,长轴所在方位即为1号极板方位P1AZ;当C25为最大时,长轴所在方位P1AZ+60;当C36为最大时,长轴所在方位P1AZ+120。2.2求解最佳逼近井眼的拟合椭圆长轴与短轴由于六臂地层倾角资料没有能代表椭圆井眼长短轴的垂直双井径曲线,需要拟合一个最佳逼近的椭圆井眼轨迹14。由于研究区的井型均为直井,实际最大井底井斜均小于3,那么同一口井中不同井段所形成的椭圆井眼可认为是共中心点,即:每个裸眼段椭图1井眼及环空示意图26条可伸缩臂结构极板分布及井眼示意 40第34卷第2期徐长凯等:基于六臂地层倾角
13、资料的井眼扩大系数对迟到时间的校正圆井眼的中心点均位于钻具正中心。椭圆井眼轨迹拟合具体步骤如下。2.2.1建立标准拟合椭圆的参数方程标准拟合椭圆的参数方程为:xa2cosyb2sin(ab0)(2)式中:为椭圆上的任何一点与原点之间连线与X轴之间的夹角,();a为椭圆长轴;b为椭圆短轴;x、y分别为椭圆上横、纵坐标值。2.2.2建立直角坐标系求解椭圆长轴与短轴当C14最大时,以井眼中心O为原点,C14为X轴建立直角坐标系(图2),由椭圆的参数方程可知,6条井径在拟合椭圆上的坐标分别为:A(a2cos0,b2sin0),B(a2cos300,b2sin300),C(a2cos240,b2sin2
14、40),D(a2cos180,b2sin180),E(a2cos120,b2sin120),F(a2cos60,b2sin60)其 中,令 cos0 K1,sin0 K1,cos300 K2,sin300K2,cos240K3,sin240K3,cos180K4,sin180K4,cos120K5,sin120K5,cos60K6,sin60K6。令f(a,b)(a2K1-C12K1)2+(a2K2-C22K2)2+(a2K3-C32K3)2+(a2K4-C42K4)2+(a2K5-C52K5)2+(a2K6-C62K6)2+(b2K1-C12K1)2+(b2K2-C22K2)2+(b2K3-
15、C32K3)2+(b2K4-C42K4)2+(b2K5-C52K5)2+(b2K6-C62K6)2当所拟合的椭圆最佳逼近实际井眼轨迹时,则满足1415:f(a,b)a0f(a,b)b0(3)对公式(3)偏导求解椭圆的长轴与短轴分别为:aC1+C43+C2+C3+C5+C612bC2+C3+C5+C64同理,当C25最大时:aC2+C53+C1+C3+C4+C612b=C1+C3+C4+C64当C36最大时:aC3+C63+C1+C2+C4+C512bC1+C2+C4+C542.3求解井眼扩大系数并修正迟到时间理论公式椭圆井眼通常采用长轴a与钻头直径D表征井径扩大率(k)8:ka-DD100%(
16、4)平均井径扩大率(-k):-k1ni=1nki(5)校正后的裸眼段环空体积(V2)4:V2(1+-k)V214(1+-k)(D2-d22)h2(6)采用平均井径扩大率校正后迟到时间(T迟)计算公式:T迟V1+V2Q()d12-d22h1+()1+-k()D2-d22h24Q(7)由于椭圆井眼不同位置处井径扩大率各不相同,仅采用长轴以“以线带面”的方式无法精准反映复杂井眼扩大情况,利用井径扩大系数校正迟到时间的公式(7)在椭圆井眼中应用效果不佳。笔者采用裸眼段井眼面积与钻头面积关系表征井眼坍塌程度,并定名为井眼扩大系数(),即裸眼段井眼面积相对钻头面积的差值与钻头面积的比值:S井眼-S钻头S钻
17、头(8)式中:为井眼扩大系数,无量纲;S井眼为裸眼段井眼面积,m2;S钻头为钻头面积,m2。最佳逼近井眼的拟合椭圆接近实际井眼形状,可近似认为:S井眼S拟椭(9)式中:S拟椭为拟合椭圆面积,m2。41录 井 工 程工 艺 技 术2023年6月由公式(8)、(9)得:S拟椭-S钻头S钻头(a/2)(b/2)-(D/2)2(D/2)2ab-D2D2(10)统计并计算该区域椭圆井眼的平均井眼扩大系数():1ni=1ni(11)由公式(8)、(11)求解裸眼段平均井眼面积(-S井眼):-S井眼()1+S钻头(12)校正后的裸眼段环空体积:V2()-S井眼-S钻具h2h2()1+D2-d224(13)采
18、用平均井眼扩大系数校正后的迟到时间T迟校计算公式:T迟校h1(d12-d22)+h2(1+)D2-d224Q(14)3应用效果为了更加准确地反映施工区域地应力集中井段的椭圆井眼情况,笔者选取了伊拉克B油田北部、中部、南部具有代表性的FH12井等3口井的电测数据,对采用上述方法计算的井眼扩大系数进行了统计分析(表1),井眼扩大系数区间为(0.207,0.545),平均井眼扩大系数 为0.461。注:1 in25.4 mm井号FH12FH13FH11井深/m4453.64454.04491.44492.24493.04454.04455.04487.14487.24487.34452.84454.
19、74490.74491.34494.0钻头直径/in8.58.58.58.58.58.58.58.58.58.58.58.58.58.58.5六臂地层倾角测井响应值/inC114.0114.5010.8511.5713.8014.5015.2312.6412.4112.1915.4812.7011.7710.7711.60C414.3014.2311.5712.8614.7714.2315.5912.8513.1113.3015.1812.0512.1311.5614.96C28.598.538.368.468.388.538.408.488.488.508.448.648.418.328.3
20、5C58.758.788.768.818.718.788.638.608.648.658.648.748.688.748.74C38.518.498.448.448.428.498.568.368.378.408.378.698.698.429.03C68.648.688.448.608.778.688.488.648.959.058.718.438.668.489.13最佳逼近拟合椭圆长轴(a)/in12.3112.4510.3111.0012.3812.4511.0911.3411.3811.3813.0711.1210.8410.2711.79短轴(b)/in8.628.628.508.
21、588.578.628.528.528.618.658.548.628.618.498.81井眼扩大系数0.4690.4850.2150.3060.4690.4850.3070.3370.3560.3620.5450.3280.2910.2070.438表1采用本文方法计算的井眼扩大系数统计分析 42第34卷第2期徐长凯等:基于六臂地层倾角资料的井眼扩大系数对迟到时间的校正以B油田Sulaiy地层为例(表2):通过对比发现,采用井眼扩大系数校正后迟到时间精准度更高,录井深度与电测深度误差减少至1.4 m以内,能够满足椭圆井眼条件下的资料采集工作。4结论(1)在录井资料采集过程中,由于无法实时跟
22、踪椭圆井眼井径变化态势,通常忽略井眼井径变化的影响,直接采用迟到时间理论公式计算,导致最终结果误差较大。(2)伊拉克B油田Sulaiy地层发育“糖葫芦串”形椭圆井眼,可采用最佳逼近椭圆拟合井眼周轨迹以反映真实情况。在满足坐标差的平方和f(a,b)关于椭圆长轴(a)和短轴(b)的偏导数均为零的条件下,利用六臂地层倾角测井资料可求解拟合椭圆的长轴与短轴。(3)针对井径扩大率表征井眼存在“以线代面”的不足,本文引入井眼扩大系数。其采用裸眼段井眼面积相对钻头面积的差值与钻头面积的比值表征井眼坍塌程度,在描述井眼方面更贴近实际情况。由于最佳逼近井眼的拟合椭圆接近实际井眼形状,可使用拟合椭圆的面积代替实际
23、井眼的面积以获取井眼扩大系数。(4)通过对文中校正公式使用前、后对比验证:伊拉克 B油田 Sulaiy地层井眼平均扩大系数为 0.461,采用井眼扩大系数对迟到时间进行校正后录井深度与电测深度误差减少至1.4 m以内,可满足本区域生产实践。参考文献 1 付志方,孙红军,董立,等.差异构造演化控制倾斜油水界面之实例:伊朗雅达油田J.海相油气地质,2017,22(3):6772.FU Zhifang,SUN Hongjun,DONG Li,et al.Controleffects on the tilted oilwater contact caused by differentialstruct
24、ural evolution:A case of Yada Oilfield,IranJ.Marine Origin Petroleum Geology,2017,22(3):6772.2 杜洋,衣英杰,辛军,等.伊朗SA油田Sarvak油藏大幅度倾斜油水界面成因探讨J.石油实验地质,2015,37(2):187193.DU Yang,YI Yingjie,XIN Jun,et al.Genesis of largeamplitude tilting oilwater contact in Sarvak Formation inSouth Azadegan Oilfield,IranJ.Pet
25、roleum Geology&Experiment,2015,37(2):187193.3 张建斌,苑仁国.迟到时间对储集层流体性质评价的影响分析与对策J.录井工程,2020,31(1):4245,51.ZHANGJianbin,YUANRenguo.Analysisandcountermeasuresoftheimpactoflagtimeontheevaluation of reservoir fluid propertiesJ.Mud LoggingEngineering,2020,31(1):4245,51.4 沈书锋,郑玉朝,李伟.井径扩大率对岩屑迟到时间的影响校正J.录井工程,20
26、09,20(1):6162.SHEN Shufeng,ZHENG Yuchao,LI Wei.Correction ofthe effect of hole enlargement ratio on lag time of cuttingsJ.Mud Logging Engineering,2009,20(1):6162.5 王金凤,邓金根,李宾.井壁坍塌破坏过程的数值模拟及井径扩大率预测J.石油钻探技术,2000,28(6):1314.WANG Jinfeng,DENG Jingen,LI Bin.Numericalsimulation of borehole collapse proced
27、ure and prediction onwellhole enlargementJ.Petroleum Drilling Techniques,2000,28(6):1314.表2Sulaiy地层迟到时间理论计算公式与校正公式计算结果对比井号FH08FH09FH17FH20FH21FH22FH29迟到时间/min理论计算50.4550.4747.9150.2048.8949.0051.02校正计算-53.6056.1754.6954.5357.45录井深度/m4458.04460.04564.04570.04462.04437.04518.0电测深度/m4465.44468.04564.84
28、571.04462.94438.24519.4误差/m7.48.00.81.00.91.21.4 43录 井 工 程工 艺 技 术2023年6月 6 刘书杰,周建良,何保生,等.欠平衡钻井条件下的井径扩大率预测技术J.中国海上油气,2007,19(5):335337.LIU Shujie,ZHOU Jianliang,HE Baosheng,et al.Theprediction technology of hole enlargement ratio in thecondition of underbalanced drillingJ.China Offshore Oiland Gas,20
29、07,19(5):335337.7 李军,陈勉,张广清,等.易坍塌地层椭圆形井眼内套管应力的有限元分析J.石油大学学报(自然科学版),2004,28(2):4548.LI Jun,CHEN Mian,ZHANG Guangqing,et al.Analysis onstress distribution of casing in sloughing formation with finiteelement methodJ.Journal of China University of Petroleum(Edition of Natural Science),2004,28(2):4548.8
30、赵凯,韩继勇,许永华.定向井井径扩大率计算模型及影响因素分析J.断块油气田,2016,23(5):638642.ZHAOKai,HANJiyong,XUYonghua.Boreholediameter enlargement ratio calculation model of directionalwells and influence factors analysisJ.FaultBlock Oil&Gas Field,2016,23(5):638642.9 曲国胜,张华,叶洪.黄骅坳陷井壁崩落法现代构造应力场测量J.华北地震科学,1993,11(3):918.QU Guosheng,ZH
31、ANG Hua,YE Hong.The principalstress orientations by the borehole breakouts in HuanghuaDepressionJ.North China Earthquake Sciences,1993,11(3):918.10 高兴军,付鑫生,于兴河,等.六臂地层倾角处理中的球面映像技术分析J.石油地球物理勘探,2002,37(2):129133.GAOXingjun,FUXinsheng,YUXinghe,etal.Analysis of spherical projection technology in six armd
32、ipmeterprocessingJ.OilGeophysicalProspecting,2002,37(2):129133.11 张志伟.地层倾角测井资料在海拉尔盆地乌尔逊油田的应用J.国外测井技术,2010(175):3436.ZHANG Zhiwei.Application of diplog data in WuerxunOilfield,HailarBasinJ.WorldWellLoggingTechnology,2010(175):3436.12 程希,许杰,王骏,等.六臂地层倾角测井仪功能电路及液压测试器实现J.石油管材与仪器,2015,1(1):8486.CHENG Xi,X
33、U Jie,WANG Jun,et al.Design of thelogging function circuits and hydraulic test device of sixarmdipmeterJ.PetroleumTubularGoods&Instruments,2015,1(1):8486.13 印兴耀,马妮,马正乾.地应力预测技术的研究现状与进展J.石油物探,2018,57(4):488504.YIN Xingyao,MA Ni,MA Zhengqian,et al.Review ofin situstresspredictiontechnologyJ.Geophysical
34、Prospecting for Petroleum,2018,57(4):488504.14 魏阳庆,刘峰,孙彩云.利用六臂地层倾角测井资料分析地应力J.断块油气田,2005,12(5):7880.WEI Yangqing,LIU Feng,SUN Caiyun.Using sixarmdipmeter log to analyze ground stressJ.Fault BlockOil&Gas Field,2005,12(5):7880.15 李岳生,黄友谦.数值逼近M.北京:人民教育出版社,1978.LIYuesheng,HUANGYouqian.NumericalapproximationM.Beijing:People s Education Press,1978.(返修收稿日期20230409编辑王丙寅)44