2. 上海波睿能源科技有限公司, 上海 201913;
3. Elsiwave Reservoir, 北京 100083;
4. 石油工业出版社, 北京 100011
2. Shanghai Borui Technology Co.Ltd., Shanghai 201913, China;
3. Elsiwave Reservoir, Beijing 100083, China;
4. Petroleum Industry Press, Beijing 100011, China
松辽盆地北部中浅层广泛发育河流相、三角洲相及湖相砂岩, 它们是大庆油田主要的储集层, 但其中有很大一部分砂岩层物性差、产能低, 应用常规开采手段难以获得有效动用[1]。
根据松辽盆地的地质特点及勘探开发实际, 本地区致密油定义[2]为:①储层地面空气渗透率小于或等于1mD(1mD≈0.987×10-3μm2), 孔隙度小于或等于12%;②赋存于油页岩或与其互层共生的致密砂岩储层中, 未经大规模、长距离运移; ③直井常规压裂产量低, 需采用水平井大规模体积压裂技术提高产量的石油资源。
为提高水平井和压裂设计水平, 大庆油田根据测井资料总结了致密油储层评价方法, 主要利用“两品质”指标评价致密油储层[2-3]。“两品质”包括“储层品质”(或称为“地质甜点”)和“工程品质”(或称为“工程甜点”)。大庆油田“工程品质”参数主要包括岩石脆性指数、最小水平主应力等。
国内外通常使用孔隙压力、岩石脆性、最小水平主应力、最大水平主应力与最小水平主应力之差(即水平应力差)、破裂压力等参数描述“工程甜点”。只有水平井同时钻遇“双甜点”即“地质甜点”和“工程甜点”, 且在“双甜点”段实施压裂, 水平井+大型水力压裂才可以获得比较高的油气产量。
前人对大庆地区致密油开发中影响水平井压裂效果的主要因素进行了研究, 与“工程甜点”有关的因素包括裂缝密度和裂缝长度。裂缝密度即裂缝的发育程度, 主要与岩石脆性和水平应力差有关。脆性越大, 越容易产生裂缝。水平应力差越小, 越容易产生复杂的裂缝系统[4-6]。龙虎泡油田龙26外扩区高台子组油层整体脆性不高, 很多压裂施工出现了压裂后岩石没有发生脆性破裂而是塑性形变现象。
2016年, 刘勇等[7]利用叠前反演的杨氏模量和泊松比计算了岩石脆性指数, 刘建伟等[8]利用杨氏模量、体积模量、剪切模量线性回归公式计算破裂压力; 2018年, 印兴耀等[9]介绍了一种利用地层曲率、杨氏模量和泊松比计算最大水平主应力的方法, 任岩等[10]利用测井资料计算岩石破裂压力、评价吉木萨尔致密油地层的可压裂性, 朱斗星等[11]利用叠后地震资料开展了“工程甜点”的预测。上述文献评价“工程甜点”时偏重于岩石脆性指数或破裂压力, 很少提及水平应力差。
本次研究重点是开展地震高分辨弹性参数反演, 获得杨氏模量、泊松比、纵波速度和密度4个弹性参数。利用这些弹性参数和本地区的地应力计算公式, 计算并预测岩石脆性、水平应力差等关键的“工程甜点”指标, 进而指导和优化水平井压裂设计。
1 技术方法 1.1 叠前反演 1.1.1 4个弹性参数反射率直接反演纵波速度、密度、杨氏模量和泊松比这4个关键的弹性参数是计算岩石脆性和水平主应力的基础数据, 因此必须开展叠前弹性参数反演。目前成熟稳定的叠前弹性参数反演技术是基于Zoeppritz方程的线性公式, 即:
$ \begin{array}{*{20}{l}} {{R_{{\rm{PP}}}}(\theta ) \approx \left( {\frac{1}{4}{{\sec }^2}\theta - 2k{{\sin }^2}\theta } \right) \cdot \frac{{\Delta E}}{E} + \left[ {\frac{1}{4}{{\sec }^2}\theta \cdot } \right.}\\ {\left. {\frac{{(2k - 3) \cdot {{(2k - 1)}^2}}}{{(4k - 3) \cdot k}} + 2k{{\sin }^2}\theta \cdot \frac{{1 - 2k}}{{3 - 4k}}} \right] \cdot \frac{{\Delta \sigma }}{\sigma } + }\\ {\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\left( {\frac{1}{2} - \frac{1}{4}{{\sec }^2}\theta } \right) \cdot \frac{{\Delta \rho }}{\rho }} \end{array} $ | (1) |
$ \begin{aligned} R_{\mathrm{PP}}(\theta) \approx & \frac{1}{2}\left(1-4 \frac{v_{\mathrm{s}}^{2}}{v_{\mathrm{P}}^{2}} \sin ^{2} \theta\right) \cdot \frac{\Delta \rho}{\rho}+\frac{1}{2 \cos ^{2} \theta}\cdot \\ & \frac{\Delta v_{\mathrm{p}}}{v_{\mathrm{p}}}-4 \frac{v_{\mathrm{s}}^{2}}{v_{\mathrm{p}}^{2}} \sin ^{2} \theta \cdot \frac{\Delta v_{\mathrm{s}}}{v_{\mathrm{s}}} \end{aligned} $ | (2) |
式中:R为反射系数; θ为入射角; vP为纵波速度; vS为横波速度; k为vS2/vP2; E为杨氏模量; σ为泊松比; ρ为密度。
公式(1)由宗兆云等给出[12], 可以直接反演杨氏模量、泊松比和密度的反射率, 而且目前只有这个公式才可以直接反演杨氏模量和泊松比这两个参数反射率。公式(2)是Aki & Rechards公式, 可以直接反演纵波速度、横波速度和密度反射率。
本次叠前弹性参数反演采用全道集参与的直接弹性反演参数方法[13], 即不是采用部分叠加而是采用全部道集数据, 也不是采用间接计算方法由纵波速度、横波速度和密度3个基本的参数计算其它弹性参数, 而是采用直接反演方法得到这3个参数, 方法本身在实现过程中加入了抗噪技术。
1.1.2 高分辨反演由于龙虎泡油田目的层高三段油层相对于地震数据分辨率较薄, 储层平均厚度4m左右, 目前地震数据的分辨率有限, 难以识别这么薄的油层, 因此需要开展高分辨率地震反演。
首先利用地震解释层位数据和测井数据在层位控制下由测井数据内插得到一个测井数据模型, 该模型包含低频和高频信息。需要建立纵波速度、密度、杨氏模量和泊松比4个弹性参数测井数据模型。本次高分辨率地震反演采用常规地震反演结果(中低频部分)与测井高频模型融合方法[14], 其中常规地震反演即是从上述测井数据模型中补充低频背景模型得到绝对值数据。
1.2 岩石脆性指数及地应力计算公式中相关系数的确定脆性和水平主应力计算需要4个关键的弹性参数:杨氏模量、泊松比、纵波速度和密度。
岩石脆性指数计算的方法有两类:一类是利用测井资料计算脆性矿物含量从而计算出脆性指数; 另一类是利用杨氏模量和泊松比计算岩石脆性指数[15]。本次研究采用后一种方法, 首先对杨氏模量和泊松比做归一化处理, 然后求二者的算术平均值得到岩石的脆性指数(BRIT-T):
$ B_{\mathrm{RIT}-E}=\frac{E-E_{\min }}{E_{\max }-E_{\min }} \times 100 \% $ | (3) |
$ B_{\mathrm{RIT}_{\sigma}}=\frac{\sigma-\sigma_{\min }}{\sigma_{\max }-\sigma_{\min }} \times 100 \% $ | (4) |
$ B_{\mathrm{RIT}-T}=\frac{B_{\mathrm{RTT}-E}+B_{\mathrm{RTT}-\sigma}}{2} $ | (5) |
式中:B为脆性指数, 即岩石中脆性矿物含量。根据研究区测井资料统计结果, 本次研究Emax、Emin、σmax、σmin分别取值为50GPa、15GPa、0.35、0.25。
计算水平应力差首先需要计算水平主应力, 水平主应力计算需要泊松比、声波速度和密度3个弹性参数。在本次研究中, 水平主应力计算模型使用黄荣樽等[16]改进的Anderson公式, 即在公式中增加一个构造系数修正项, 使得水平主应力不仅仅与上覆岩层压力有关, 也与构造应力有关, 即:
$ P_{\mathrm{h}}=\left(\frac{\sigma}{1-\sigma}+\xi_{1}\right)\left(P_{\mathrm{o}}-\alpha P_{\mathrm{P}}\right)+\alpha P_{\mathrm{P}} $ | (6) |
$ P_{\mathrm{H}}=\left(\frac{\sigma}{1-\sigma}+\xi_{2}\right)\left(P_{\mathrm{o}}-\alpha P_{\mathrm{p}}\right)+\alpha P_{P} $ | (7) |
式中:Ph为最小水平主应力; PH为最大水平主应力; Po为上覆岩层压力; PP为孔隙压力; ξ1, ξ2为构造应力系数; α为有效应力系数, 也称Boit系数。
根据实验室差应变法分析本地区地应力及岩石力学参数, 确定公式中两个重要参数ξ1和ξ2分别取值为0.35、0.53。α可由测井参数计算得到(即由测井体积模量K计算得到); 也可由泥质含量VCL得到。一般认为, 当VCL < 0.3时, α=0.8;当0.3≤VCL < 0.5时, α=0.65;当0.5≤VCL≤1时, α=0.5。本区VCL为0.35~0.40, 故本次研究α取值为0.65。
上覆岩层压力Po由叠前反演的密度数据体积分得到, 孔隙压力PP由叠前反演的纵波速度vP结合上覆岩层压力Po计算得到。泊松比σ通过叠前反演得到。
本次研究的技术思路及流程如图 1所示。
龙虎泡油田位于大庆长垣构造带西部, 中间夹齐家—古隆凹陷。龙虎泡油田除局部地区高零段、高一段含油外, 其它地区主力油层均集中于高三段和高四段油层, 尤其是高三段。总体上, 龙26试验区高台子组高三段油层平均渗透率为1.23mD, 平均孔隙度为13.9%, 平均厚度为3.9m, 为中低孔、偏低渗致密储层。该区高三段油层在压裂前普遍无自然产能, 必须采用水平井体积压裂增产措施才能获得有经济效益的产能。
2.2 脆性指数及应力差预测利用龙虎泡油田龙26试验区4口直井(钻达目的层高三段油组)测井资料建立测井弹性参数模型, 针对每个弹性参数反射率数据体, 分别开展高分辨反演, 得到4个弹性参数高分辨反演成果数据体, 结果如图 2所示。
利用公式(3)、公式(4)和公式(5)计算岩石脆性三维数据体, 利用公式(6)和公式(7)分别计算最小和最大水平主应力, 并得到水平应力差。
高岩石脆性指数是形成裂缝的必要条件。一般认为脆性破坏是指岩石在弹性变形阶段结束、没有或仅有很小的塑性变形时发生的破坏。如果岩石在构造应力作用下容易发生脆性破坏从而形成裂缝, 那么岩石的泊松比和杨氏模量要同时满足条件:泊松比<0.3[5]且杨氏模量E为30~80GPa, 即脆性指数为40%左右。
图 3是龙26试验区高三段主力油层脆性指数平面分布图, 图中大多数开发井轨迹钻遇地层的脆性指数都大于40%, 具有较高的脆性指数, 如平35、平38和平48井。图 4是平26井高三段脆性指数剖面, 钻遇储层段的岩石脆性指数都大于45%。
水平应力差越小, 在压裂时越容易形成更复杂的裂缝, 即产生的裂缝并不仅仅在某一个方向, 而是在更多的方向产生裂缝从而形成一个复杂的裂缝系统。这样复杂的裂缝系统更容易沟通众多微小孔隙, 从而扩大储集空间并增强渗透性。同时较小的水平应力差还有助于保持井身的稳定性, 不易导致井壁垮塌等钻井事故。由在四川盆地的实践经验可知, 水平应力差一般小于15MPa时比较有利。图 5、图 6分别为过平26井高三段主力油层最大、最小水平主应力剖面。图 7是龙26试验区高三段水平应力差平面分布图, 图中水平应力差都小于5MPa, 从水平应力差这个角度看, 龙26试验区高三段水平应力差完全满足压裂的要求。图 8是过平26井高三段水平应力差剖面, 可见, 钻遇高三段的水平应力差都小于5MPa。
今后在龙26试验区进行水平井压裂设计时, 应该优先在水平应力差小(<15MPa)、岩石脆性指数高(>40%)的水平段实施压裂, 这样才能获得更好的压裂效果。对于岩石脆性指数和水平应力差两个“工程甜点”指标, 岩石脆性指数更加重要, 因为在龙26试验区水平应力差都较小。
2.3 对井误差分析利用研究区平26、平14h和平48h这3口井资料对脆性指数和水平应力差两个指标参数的预测结果进行误差统计分析, 结果表明:脆性指数预测结果相对误差小于10%, 水平应力差预测结果相对误差小于12%, 预测精度满足生产需要。
2.4 压裂效果分析龙虎泡油田生产井如果不做压裂则普遍无自然产能。经过水平井体积压裂改造后, 龙26试验区3口水平井(平35、平38、平48)产能最高增加到8~12t/d, 平均产能2.0~3.7t/d, 增产效果明显(表 1)。这3口井的脆性指数都在40%~50%, 脆性指数较高。大庆油田在该区压裂实践经验表明, 大多数井压裂后增产效果较好的井段的脆性指数通常也都大于40%。
本文利用叠前反演得到的4个弹性参数(纵波速度、密度、杨氏模量和泊松比), 采用龙虎泡油田本地的岩石脆性指数计算公式和水平主应力计算公式, 计算得到岩石脆性指数、最大水平主应力、最小水平主应力和水平应力差, 油田实际数据计算的最大相对误差为12%, 证明本文方法可行。
龙虎泡油田主力构造区高台子高三段脆性指数为35%~50%, 水平应力差都小于5MPa。压裂生产实践表明, 压裂后增产效果明显的水平段脆性指数都大于40%且水平应力差较小, 通常小于5MPa。平48等3口生产井压裂后增产效果明显。
建议在水平井轨迹设计和部署时参考“工程甜点”预测成果以提高压裂成功率。本文研究思路和技术适用于类似地区致密储层脆性指数和水平应力差的预测研究。
[1] |
崔保文, 林铁峰, 董万百, 等. 松辽盆地北部致密油水平井技术及勘探实践[J]. 大庆石油地质与开发, 2014, 33(5): 17-22. CUI B W, LIN T F, DONG W B, et al. Horizontal well techniques and their exploration practices in the tight oil reservoir in north Songliao basin[J]. Petroleum Geology and Oilfield Development in Daqing, 2014, 33(5): 17-22. |
[2] |
王玉华, 蒙启安, 梁江平, 等. 松辽盆地北部致密油勘探[J]. 中国石油勘探, 2015, 20(4): 44-52. WANG Y H, MENG Q A, LIANG J P, et al. Tight oil exploration in northern Songliao basin[J]. China Petroleum Exploration, 2015, 20(4): 44-52. DOI:10.3969/j.issn.1672-7703.2015.04.005 |
[3] |
金成志. 松辽盆地北部陆相非均质致密油整体勘探开发"2+3"模式[J]. 大庆石油地质与开发, 2014, 33(5): 9-15. JIN C Z. Overall exploration and development "2+3"model of the continental heterogeneous tight oil in north Songliao basin[J]. Petroleum Geology and Oilfield Development in Daqing, 2014, 33(5): 9-15. DOI:10.3969/J.ISSN.1000-3754.2014.05.002 |
[4] |
朱维耀, 王明, 刘合, 等. 特低渗油层非达西渗流压裂水平井非定常渗流产能预测新方法[J]. 特种油气藏, 2010, 17(5): 63-66. ZHU W Y, WANG M, LIU H, et al. A new method for predicting unsteady production capacity of non-darcy flow fracturing horizontal wells in extra-low permeability reservoirs[J]. Special Oil & Gas Reservoirs, 2010, 17(5): 63-66. DOI:10.3969/j.issn.1006-6535.2010.05.017 |
[5] |
张年学, 盛祝平, 李晓, 等. 岩石泊松比与内摩擦角关系研究[J]. 岩石力学与工程学报, 2011, 30(增1): 2600-2609. ZHANG N X, SHENG Z P, LI X, et al. Study of relationship between Poisson's ratio and angle of internal friction for rocks[J]. Chinese Journal of Rock Mechanics and Engineering, 2011, 30(S1): 2600-2609. |
[6] |
崔明月, 刘玉章, 修乃领, 等. 形成复杂缝网体积的影响因素分析[J]. 石油钻采工艺, 2012, 36(2): 82-87. CUI M Y, LIU Y Z, XIU N L, et al. The analysis of influencing Factors of forming complex fracture network volume[J]. Oil Driling & Production Technology, 2012, 36(2): 82-87. |
[7] |
刘勇, 方伍宝, 李振春, 等. 基于叠前地震的脆性预测方法及应用研究[J]. 石油物探, 2016, 55(3): 425-432. LIU Y, FANG W B, LI Z C, et al. Brittleness prediction and application based on pre-stack seismic inversion[J]. Geophysical Prospecting for Petroleum, 2016, 55(3): 425-432. DOI:10.3969/j.issn.1000-1441.2016.03.013 |
[8] |
刘建伟, 张云银, 曾联波, 等. 非常规油藏地应力和应力甜点地球物理预测——渤南地区沙三下亚段页岩油藏勘探实例[J]. 石油地球物理勘探, 2016, 51(4): 792-800. LIU J W, ZHENG Y Y, ZENG L B, et al. Geophysical prediction of stress and stress desserts in unconventional reservoirs:An example in Bonan area[J]. Oil Geophysical Prospecting, 2016, 51(4): 792-800. |
[9] |
印兴耀, 马妮, 马正乾, 等. 地应力预测技术的研究现状与进展[J]. 石油物探, 2018, 57(4): 488-504. YIN X Y, MA N, MA Z Q, et al. Review of in-situ stress prediction technology[J]. Geophysical Prospecting for Petroleum, 2018, 57(4): 488-504. |
[10] |
任岩, 曹宏, 姚逢昌, 等. 吉木萨尔致密油储层脆性及可压裂性预测[J]. 石油地球物理勘探, 2018, 53(3): 511-519. REN Y, CAO H, YAO F C, et al. Brittleness and fracability prediction for tight oil reservoir in Jimsar Sag, Junggar Basin[J]. Oil Geophysical Prospecting, 2018, 53(3): 511-519. |
[11] |
朱斗星, 蒋立伟, 牛卫涛, 等. 页岩气地震地质工程一体化技术的应用[J]. 石油地球物理勘探, 2018, 53(S1): 249-255. ZHU D X, JIANG L W, NIU W T, et al. Seismic and geological integration applied in the shale gas exploration[J]. Oil Geophysical Prospecting, 2018, 53(S1): 249-255. |
[12] |
宗兆云, 印兴耀, 张峰, 等. 杨氏模量和泊松比反射系数近似方程及叠前地震反演[J]. 地球物理学报, 2012, 55(11): 3786-3794. ZONG Z Y, YIN X Y, ZHANG F, et al. Reflection coefficient equation and pre-stack seismic inversion with Young's modulus and Poisson's ratio[J]. Chinese Journal of Geophysics, 2012, 55(11): 3786-3794. DOI:10.6038/j.issn.0001-5733.2012.11.025 |
[13] |
杨建礼, 常新伟.提高非气藏"亮点"和非"亮点"气藏弹性预测识别精度的方法[R].第十七届中国科协年会, 广州, 2015 YANG J L, CHANG X W.The method of prediction accuracy improving of non-gas bright-spot and non-bright-spot gas[R].Essays of 17th Annual Meeting of China Association for Science and Technology, Guangzhou, 2015 |
[14] |
李庆忠. 论地震约束反演的策略[J]. 石油地球物理勘探, 1998, 33(4): 424-438. LI Q Z. On strategy of seismic restricted inversion[J]. Oil Geophysical Prospecting, 1998, 33(4): 424-438. |
[15] |
RICKMAN R, MULLEN M J, PETRE J E.A practical use of shale petrophysics for stimulation design optimization: All shale plays are not clones of the barnett shale[R].SPE Annual Technical Conference and Exhibition, Colorado USA, 2008
|
[16] |
黄荣樽, 庄锦江. 一种新的地层破裂压力预测方法[J]. 石油钻采工艺, 1986, 9(3): 1-14. HUANG R Z, ZHUANG J J. A new method for predicting formation fracture pressure[J]. Oil Driling & Production Technology, 1986, 9(3): 1-4. |