2. 中国石油集团测井有限公司长庆分公司
2. Changqing Branch of CNPC Logging Co., Ltd.
越来越多的油气田采用水平井开发技术,水平井对于提高产能和油气田高效开发意义重大。由于水平井与直井中测井仪器响应不同,目前尚未形成较完善的水平井测井解释方法,地层各向异性是引起水平井测井响应复杂的因素之一[1-5],水平井周围会部署多口控制井,且控制井一般为直井,利用水平井与邻近控制井测井响应对比可以研究地层各向异性。水平井最常见的测井系列有电阻率测井、自然伽马测井、声波时差测井,由于水平井与邻近控制井电阻率测井系列往往不一致,地层各向异性对自然伽马测井影响很小,因此,水平井与邻近控制井电阻率和自然伽马对比实际意义不大,但是受地层沉积物颗粒大小、沉积韵律影响,声波时差在地层纵向和横向上差异较大[6-9],水平井与邻近控制井相同地层的声波时差对比,可更好地研究地层声波时差各向异性。本文在鄂尔多斯盆地陇东地区水平井与邻近控制井测井数据分析基础上,提出了基于测井曲线对比的水平井声波时差各向异性校正方法,这将有助于提高水平井储层参数测井计算精度。
1 储层地质特征及水平井部署概况鄂尔多斯盆地陇东地区水平井目的层为延长组长6X层。该层主要由深灰色细砂岩、粉砂岩、钙质砂岩、灰黑色泥岩、粉砂质泥岩组成,孔隙度为8%~12%,渗透率为0.1~0.3mD,长6X储层物性差,属于低孔、低渗致密油储层[10-12]。主力储层砂体向东北—西南方向延伸,砂体宽度为150~400m。陇东地区水平井水平段长度为800~1500m,水平段井斜角一般大于85°,水平段方位平均为162°,周围控制井一般有3~5口[13],控制井一般为直井。
2 声波时差各向异性校正方法由于复杂的测井环境、地层各向异性等因素影响,水平井测井响应远不同于直井[14-16],要正确评价水平井测井资料,就要充分利用控制井测井资料,找出在岩性接近条件下水平井与控制井之间差异,利用两者之间差异分析储层各向异性。岩层沉积时相同岩性地层在平面上总会延伸一定的距离,地层横向可对比性强,陇东地区水平井与邻近控制井目的层相同垂深处最小距离为200~400m,岩性在平面上变化不大[17-20],地层拉平条件下相同深度处控制井目的层的测井数值,可近似认为是仪器垂直条件下的响应,而此深度下水平井测井数值则是仪器水平条件下的响应。
水平井在目的层钻进时,水平段钻头并不是完全水平钻进,轨迹一般会沿着目的层上下、左右偏移,水平井与控制井相同深度平面上距离最近两点的测井数值最具有可对比性,因为这两点岩性最为接近。计算水平段上每一点与相同深度平面上多个控制井距离,取距离最小的控制井测量值与同深度点上水平段测量值进行对比,将水平段上每一测量点逐一与控制井对比,控制井越多,得到的对比数值越可靠。由于水平段延伸一定长度,测井数据量大,因此,在目的层相同深度条件下,通过距离最近的控制井与水平段声波时差对比,可以分析大量测量点垂直方向与水平方向声波时差差异,根据分析结果可进一步建立水平井声波时差校正方法,用校正后声波时差就能较为准确地计算出水平段物性参数。
图 1是地层拉平条件下水平井与两口控制井轨迹示意图,S井为水平井,A井、B井为两口控制井,控制井一般为直井,F上、F下为目的层上、下界面,S1、S2为目的层水平段上两点,通常要求水平段上各点井斜角大于85°,S1、S2垂深在控制井A井上分别对应A1、A2,在控制井B井上分别对应B1、B2,计算A1S1、B1S1之间距离,并比较A1S1、B1S1大小,当A1S1<B1S1时,就将控制井A1点测井数值与水平井S1点测井数值提取出来进行对比; 当A1S1>B1S1时,就将控制井B1点测井数值与水平井S1点测井数值提取出来进行对比; 当A1S1=B1S1时,S1点测井数值与A1、B1测井数值之一进行对比,同样也可以提取出与水平井S2点相对比的数值。水平段上每一点都可以提取出与之对比的控制井数值,水平井与控制井对比曲线实现过程具体如下:
(1) 设水平井井口坐标为(xS0, yS0, zS0)、第i个(i≥1)测量点斜深为dSi(m)、井斜角为αSi(°)、方位角为φSi(°),第i个测量点测井曲线值为λS(xSi, ySi, zSi)。依据以上已知数据,用公式(1)至公式(3)求出邻近测量点间隔Δl(m)及第i个测量点东西向位移xeSi(m)、南北向位移ynSi(m)。
$ \Delta l=d_{\mathrm{S} i}-d_{\mathrm{S}(i-1)} $ | (1) |
$ x e_{\mathrm{S} i}=x e_{\mathrm{S}(i-1)}+\Delta l \times \sin \left[\frac{\alpha_{\mathrm{S} i}+\alpha_{\mathrm{S}(i-1)}}{2}\right] \times \sin \left[\frac{\varphi_{\mathrm{S} i}+\varphi_{\mathrm{S}(i-1)}}{2}\right] $ | (2) |
$ y n_{\mathrm{S} i}=y n_{\mathrm{S}(i-1)}+\Delta l \times \sin \left[\frac{\alpha_{\mathrm{S} i}+\alpha_{\mathrm{S}(i-1)}}{2}\right] \times \cos \left[\frac{\varphi_{\mathrm{S} i}+\varphi_{\mathrm{S}(i-1)}}{2}\right] $ | (3) |
上式中当i=1时,
(2) 用公式(4)、公式(5)计算第i个测量点水平位移sSi(m)、水平方位角θSi(°)。
$ \theta_{\mathrm{S} i}=\frac{180}{\pi} \arctan \left(\frac{x e_{\mathrm{S} i}}{y n_{\mathrm{S} i}}\right) $ | (4) |
$s_{\mathrm{S} i}=\sqrt{x e_{\mathrm{S} i}^{2}+y n_{\mathrm{S} i}^{2}} $ | (5) |
(3) 用公式(6)至公式(8)计算第i个测量点坐标(xSi, ySi, zSi)。
$ x_{\mathrm{S} i}=x_{\mathrm{S} 0}+s_{\mathrm{S} i} \times \sin \theta_{\mathrm{S} i} $ | (6) |
$ y_{\mathrm{S} i}=y_{\mathrm{S} 0}+s_{\mathrm{S} i} \times \cos \theta_{\mathrm{S} i} $ | (7) |
$ z_{\mathrm{S} i}=z_{\mathrm{S}(i-1)}+\Delta l \times \cos \left[\frac{\alpha_{\mathrm{S} i}+\alpha_{\mathrm{S}(i-1)}}{2}\right] $ | (8) |
(4) 水平井曲线校深后,水平段上每一点斜深dSi可能对应一个垂深vSi(m),也可能多个斜深dSi对应一个垂深vSi,根据vSi最大值(vSmax)和最小值(vSmin)确定出水平段轨迹在目的层上、下移动量Δv(m),Δv=ΔvSmax-ΔvSmin。
(5)(xA0,yA0)为A控制井井口坐标,水平井与A井补心海拔差为ΔESA(m),A井与水平井补心海拔对齐后,A井每一个测量点深度为zAi+ΔESA,zAi(m)为A井测井记录深度。由于水平井与A井目的层高程差为ΔHSA(m),因此,A井与水平井目的层对齐后,A井每一个测量点深度从zAi+ΔESA变为zAi+ΔESA-ΔHSA,将zAi+ΔESA-ΔHSA记作zAci,经过海拔差和高程差校正后,控制井A井每一个测量点坐标为(xA0, yA0, zAci)。由水平段vSi最大值和最小值确定出对应A控制井zAci变化范围,在zAci变化范围内,控制井A井每一点测井曲线值记作λA(zAci); 同样,控制井B井经过海拔差和高程差校正后,每一个测量点坐标为(xB0, yB0, zBci),测井曲线值记作λB(zBci)。
(6) 提取出水平段上每一点在相同深度条件下对应最近控制井的测量值。例如,水平段上S1点坐标为(xS1, yS1, zS1),在控制井A井和B井上分别找出与zSi深度相同的zAci和zBci,在zSi深度平面上计算(xS1, yS1)与(xA0, yA0)之间距离,记作DA1(m),由公式(9)计算。
$ \mathrm{D}_{\mathrm{A} 1}=\sqrt{\left(x_{\mathrm{S} 1}-x_{\mathrm{A} 0}\right)^{2}+\left(y_{\mathrm{S} 1}-y_{\mathrm{A} 0}\right)^{2}} $ | (9) |
由公式(10)计算(xS1, yS1)与(xB0, yB0)之间距离,记作DB1(m)。
$\mathrm{D}_{\mathrm{B} 1}=\sqrt{\left(x_{\mathrm{S} 1}-x_{\mathrm{B} 0}\right)^{2}+\left(y_{\mathrm{S} 1}-y_{\mathrm{B} 0}\right)^{2}} $ | (10) |
比较DA1与DB1大小,当DA1<DB1时,将S1点原始测量值λS(xS1, yS1, zS1)与λA(zAci)对比; 当DA1 > DB1时,将λB(zBci)与S1点原始数值进行对比; DA1=DB1时,S1点原始数值与λA(zAci)或λB(zBci)任意一点对比。当有3口以上控制井时,计算相同深度平面上多个水平井与控制井间距离,提取距离最小控制井上同深度点的测量值与水平段上测量值进行对比。
用以上方法提取出水平段上每一点相同深度条件下对应最近控制井的数值,通过程序实现上述方法,就可完成水平井与控制井测量值快速对比。控制井越多,得到的对比数值越可靠。
3 校正方法应用实例为了精确对比水平井与控制井数值,建立储层各向异性校正公式,需要尽可能对比相同深度条件下,水平段与控制井距离较小范围内两点测井数值,为此设计一个距离阀值,这一阀值可根据砂体在平面上最小变化距离来确定。对于陇东地区长6X储层,地质分析认为在平面215m范围内地层岩性、物性变化不大,因此,计算相同深度条件下水平段上每一点与控制井距离,将距离小于或等于215m的所有点集中在一起分析,通过这些点水平段与控制井测量值对比,建立水平井声波时差校正方法。
C17-2A井是鄂尔多斯盆地陇东地区部署的一口水平井,目的层为延长组长6X层,水平段长度为1152m,井口坐标x: 4044317.910,y: 36486762.160,该水平井邻近部署了3口控制井,控制井目的层平均厚度为16.5m,用本文的计算方法提取出水平段上每一点对应相同垂深下最近控制井的测量值,筛选出与控制井距离小于或等于215m的数据点4685个,这些测量点在水平段上自然伽马平均值为85.1API,声波时差平均值为209.2μs/m,对应控制井自然伽马平均值为83.2API,声波时差平均值为232.2μs/m,研究区目的层水平段声波时差比垂直井小,平均偏小23μs/m,呈现出较强的各向异性,自然伽马在水平井与直井中差异不明显。为了分析岩性对于声波时差各向异性的影响,将统计的4685个对比数据依据自然伽马数值大小,按间隔10API分区间统计,从表 1中可以看出随着水平井自然伽马值升高,垂直方向与水平方向声波时差差值逐渐增大,声波各向异性逐渐增强,声波时差差值与自然伽马相关性较好。
当研究区目的层自然伽马值低于80API时,一般认为是泥质含量较低、物性较好的储层,表 1的数据表明,这些较好储层水平段声波时差比垂直井低5.8~17.3μs/m,平均为11.55μs/m,直接应用水平段声波时差计算的孔隙度将会减小3.8%。研究区目的层14块饱含流体垂直岩样纵波时差实验测量值平均为226μs/m,14块水平岩样纵波时差实验测量值平均为212μs/m,两者相差14μs/m,本文提供的方法对比结果与实验接近,可以应用到实际测井解释中。
用同样方法分析了研究区C23-1A井等5口水平井与其控制井声波时差差异,将每口井水平段上相同垂深下与控制井距离小于或等于215m的数据点集中统计,同样依据水平井自然伽马值大小,按每间隔10API分区间统计,通过回归分析得到研究区目的层垂直方向和水平方向声波时差差值ΔAC与水平井自然伽马GR关系式(图 2):
$ \Delta A C=0.697 G R-38.878 \quad\left(R^{2}=0.802\right) $ | (11) |
由公式(11)和实测水平井声波时差(AC水平井)可以得到校正后的声波时差(AC校正)。
$ A C_{\text {校正 }}=A C_{\text {水平井 }}+\Delta A C $ | (12) |
该地区已有直井条件下声波时差(AC)与孔隙度(ϕ)经验关系式如下:
$\phi=0.16 A C_{\text {校正 }}-26.29 $ | (13) |
根据公式(13)和校正后的声波时差计算得到水平段储层孔隙度。对研究区16口井声波时差做了校正,图 3是C17-2A井2700.0~2792.0m测井解释综合图,该井2704.0~2718.0m、2718.0~ 2777.3m、2777.3~2788.0m分别解释为差油层、油层、干层,应用本文提供的方法计算水平段声波时差校正量,得到了校正后的声波时差,应用校正后的声波时差计算2704.0~2718.0m、2718.0~ 2777.3m、2777.3~2788.0m孔隙度平均值分别为10.25%、11.02%、4.15%,应用实际测量声波时差计算的孔隙度分别为6.02%、8.79%、3.79%,计算结果普遍偏小。因此,水平井声波时差各向异性校正很有必要,这将有助于提高水平井储层参数测井计算精度。
(1) 声波时差是计算储层孔隙度的重要测井曲线,受地层沉积特征影响,声波时差各向异性特征显著,在水平井与邻近控制井相同地层声波时差对比分析基础上,提出了水平井声波时差各向异性校正方法和步骤,建立了鄂尔多斯盆地陇东地区长6X层水平井声波时差各向异性校正公式。
(2) 陇东地区长6X层水平井与控制井对比分析表明,水平段声波时差比邻近直井普遍偏低,邻近直井与水平段声波时差差值与水平井自然伽马关系密切,随着水平井自然伽马升高,声波时差差值逐渐增大,声波各向异性逐渐增强。通过建立的声波时差校正公式计算了目的层孔隙度,计算结果明显提高了水平井储层参数测井计算精度。
[1] |
覃世银, 管志宁, 王昌学, 等. 层状各向异性地层的识别与评价[J]. 测井技术, 2003, 27(3): 194-197. Qin Shiyin, Guan Zhining, Wang Changxue, et al. Identification and evaluation of layered anisotropic formations[J]. Well Logging Technology, 2003, 27(3): 194-197. DOI:10.3969/j.issn.1004-1338.2003.03.005 |
[2] |
孔繁达, 何涛, 姜洪福. 薄互层水平井测井电阻率的各向异性校正[J]. 北京大学学报(自然科学版), 2011, 47(2): 302-308. Kong Fanda, He Tao, Jiang Hongfu. Anisotropic correction for resistivity logging of horizontal well in thin interlayer reservoir[J]. Acta Scientiarum Naturalium Universitatis Pekinensis, 2011, 47(2): 302-308. |
[3] |
乔悦东, 孙建孟, 耿尊博. 斜井泥岩声波速度各向异性校正新方法研究[J]. 石油天然气学报, 2010, 32(5): 104-108. Qiao Yuedong, Sun Jianmeng, Geng Zunbo. New methods of shale acoustic velocity anisotropy correction in deviated wells[J]. Journal of Oil and Gas Technology, 2010, 32(5): 104-108. |
[4] |
刘国强. 非常规油气时代的测井采集技术挑战与对策[J]. 中国石油勘探, 2021, 26(5): 24-37. Liu Guoqiang. Challenges and countermeasures of well logging data acquisition technology in unconventional petroleum exploration and development[J]. China Petroleum Exploration, 2021, 26(5): 24-37. |
[5] |
邹才能, 赵群, 王红岩, 等. 非常规油气勘探开发理论技术助力我国油气增储上产[J]. 石油科技论坛, 2021, 40(3): 72-79. Zou Caineng, Zhao Qun, Wang Hongyan, et al. Theory and technology of unconventional oil and gas exploration and development helps China increase oil and gas reserves and production[J]. Petroleum Science and Technology Forum, 2021, 40(3): 72-79. |
[6] |
邓继新, 史哥, 刘瑞珣, 等. 泥岩、页岩声速各向异性及其影响因素分析[J]. 地球物理学报, 2004, 47(5): 862-868. Deng Jixin, Shi Ge, Liu Ruixun, et al. Analysis of the velocity anisotropy and its affection factors in shale and mudstone[J]. Journal of Geophysics, 2004, 47(5): 862-868. |
[7] |
周灿灿, 王昌学. 水平井测井解释技术综述[J]. 地球物理学进展, 2006, 21(1): 152-160. Zhou Cancan, Wang Changxue. Review of horizontal well logging interpretation technology[J]. Progress in Geophysics, 2006, 21(1): 152-160. |
[8] |
赵江青, 王成龙. 岩石各向异性在水平井测井解释中的应用[J]. 测井技术, 1998, 22(1): 36-41. Zhao Jiangqing, Wang Chenglong. Laboratory measurement of rock anisotropy and its application to horizontal well log interpretation[J]. Well Logging Technology, 1998, 22(1): 36-41. |
[9] |
时建超, 屈雪峰, 雷启鸿, 等. 致密油水平井声波时差测井影响因素分析及测井响应特征研究: 以鄂尔多斯盆地陇东地区长7储层为例[J]. 西北大学学报(自然科学版), 2017, 47(4): 586-592. Shi Jianchao, Qu Xuefeng, Lei Qihong, et al. The influencing factors and response characteristics of acoustic time logging in tight oil horizontal well: a case study of Chang 7 reservoir in Longdong area of Ordos Basin[J]. Journal of Northwest University(Natural Science Edition), 2017, 47(4): 586-592. |
[10] |
夏青松, 田景春, 张锦泉, 等. 鄂尔多斯盆地陇东地区三叠系延长组长6—长8储层评价及有利区带预测[J]. 油气地质与采收率, 2003, 10(4): 11-14. Xia Qingsong, Tian Jingchun, Zhang Jinquan, et al. Reservoir evaluation and favorable area prospecting in Chang 6 to Chang 8 of the Triassic period at Longdong region in Ordos Basin[J]. Petroleum Geology and Recovery Efficiency, 2003, 10(4): 11-14. |
[11] |
付金华, 刘显阳, 李士祥, 等. 鄂尔多斯盆地三叠系延长组长7段页岩油勘探发现与资源潜力[J]. 中国石油勘探, 2021, 26(5): 1-11. Fu Jinhua, Liu Xianyang, Li Shixiang, et al. Discovery and resource potential of shale oil of Chang 7 member, Triassic Yanchang Formation, Ordos Basin[J]. China Petroleum Exploration, 2021, 26(5): 1-11. |
[12] |
王志坤, 王多云, 郑希民, 等. 陕甘宁盆地陇东地区三叠系延长统长6—长8储层沉积特征及物性分析[J]. 天然气地球科学, 2003, 14(5): 380-386. Wang Zhikun, Wang Duoyun, Zheng Ximin, et al. Depositional characteristics and physical property analysis of the Chang 6-8 reservoir of Triassic Yanchang Formation in Longdong area Shanganning[J]. Natural Gas Geoscience, 2003, 14(5): 380-386. |
[13] |
莫邵元, 何顺利, 王帅, 等. 华庆超低渗透油藏压裂水平井适应性分析[J]. 科学技术与工程, 2013, 13(12): 3256-3259. Mo Shaoyuan, He Shunli, Wang Shuai, et al. Feasibility analysis of fractured-horizontal well development in Huaqing ultra-low permeability reservoir[J]. Science Technology and Engineering, 2013, 13(12): 3256-3259. |
[14] |
徐波, 汪忠浩, 伍东, 等. 基于常规电测井资料的各向异性储层水平井测井解释[J]. 测井技术, 2018, 42(6): 653-657. Xu Bo, Wang Zhonghao, Wu Dong, et al. Log interpretation of horizontal well in anisotropic reservoir based on wire line resistivity logging[J]. Well Logging Technology, 2018, 42(6): 653-657. |
[15] |
李贤胜, 刘向君, 梁利喜, 等. 层理性页岩声波各向异性校正方法研究[J]. 油气藏评价与开发, 2020, 10(5): 49-54. Li Xiansheng, Liu Xiangjun, Liang Lixi, et al. Correction methods for acoustic anisotropy of bedding shale[J]. Reservoir Evaluation and Development, 2020, 10(5): 49-54. |
[16] |
侯振坤, 杨春和, 郭印同, 等. 单轴压缩下龙马溪组页岩各向异性特征研究[J]. 岩土力学, 2015, 36(9): 2541-2550. Hou Zhenkun, Yang Chunhe, Guo Yintong, et al. Experimental study on anisotropic properties of Longmaxi Formation shale under uniaxial compression[J]. Rock and Soil Mechanics, 2015, 36(9): 2541-2550. |
[17] |
邓秀芹. 鄂尔多斯盆地三叠系延长组超低渗透大型岩性油藏成藏机理研究[D]. 西安: 西北大学, 2011. Deng Xiuqin. Accumulation mechanism research on ultra-low-permeability and large scale lithogical reservoirs of Triassic Yanchang Formation in Ordos Basin[D]. Xi'an: Northwest University, 2011. |
[18] |
付金华, 董国栋, 周新平, 等. 鄂尔多斯盆地油气地质研究进展与勘探技术[J]. 中国石油勘探, 2021, 26(3): 19-40. Fu Jinhua, Dong Guodong, Zhou Xinping, et al. Research progress of petroleum geology and exploration technology in Ordos Basin[J]. China Petroleum Exploration, 2021, 26(3): 19-40. |
[19] |
李忠兴, 李松泉, 廖广志, 等. 长庆油田超低渗透油藏持续有效开发重大试验攻关探索与实践[J]. 石油科技论坛, 2021, 40(4): 1-11. Li Zhongxing, Li Songquan, Liao Guangzhi, et al. Practice from significant experiment on effective development of Changqing ultra-low permeability reservoirs[J]. Petroleum Science and Technology Forum, 2021, 40(4): 1-11. |
[20] |
付锁堂, 付金华, 牛小兵, 等. 庆城油田成藏条件及勘探开发关键技术[J]. 石油学报, 2020, 41(7): 777-793. Fu Suotang, Fu Jinhua, Niu Xiaobing, et al. Accumulation conditions and key exploration and development technologies in Qingcheng Oilfield[J]. Acta Petrolei Sinica, 2020, 41(7): 777-793. |