Characterization Method and Application of Electrical Imaging Logging in Conglomerate Reservoir: A Case Study in Mahu Sag of Junggar Basin
-
摘要: 准噶尔盆地玛湖凹陷砾岩油藏电成像测井资料丰富,为了充分利用高分辨率电成像在岩石组分、非均质性以及储层结构表征方面的优势,开展了砾岩储层电成像测井资料处理方法及深化应用研究。依据砾岩储层电成像上岩石组分电阻率的差异,通过岩心标定确定了电阻率截止值,构建了砾石、砂质和泥质含量的计算方法及伪粒度概率累计曲线的粒度分析方法;通过数理统计,构建了分选系数、孔隙度、裂缝孔隙度、高分辨率电阻率的计算方法。基于对砾岩储层主控因素的认识,提出了利用岩性组分因子、电成像测井孔隙度、非均质性因子、裂缝因子等参数评价砾岩储层储集性能的储能指数,构建了基于电成像视地层水电阻率谱方差与储层指数识别储层流体性质的新方法。研究结果表明,电成像测井资料可以有效评价砾岩储层的物性、含油性,应用于油田勘探评价,能取得较好的效果,为砾岩储层电成像测井资料的深化应用提供了借鉴。Abstract: Because the conglomerate reservoir in Mahu Sag of Junggar Basin has rich electrical imaging logging data, the processing method and further applications of electrical imaging logging data for conglomerate reservoirs were studied to make full use of the advantage of high-resolution electrical imaging in rock composition, heterogeneity, and reservoir structure characterization. Considering the resistivity difference of rock compositions in electrical ima-ging of conglomerate reservoirs, the cut-off value of resistivity was determined through core calibration, and the calcu-lation method for the relative content of gravel, sand, and argillaceous parts as well as the grain size analysis method of cumulative pseudo-grain-size probability curves were constructed. Through mathematical statistics, the methods for calculating the sorting coefficient, porosity, fracture porosity, and high-precision resistivity were established. Based on the understanding of the main controlling factors of conglomerate reservoirs, it was proposed that the energy storage indexes for evaluating the reservoir performance of conglomerate reservoir could be used to develop a new method for identifying reservoir fluid properties by the variance of apparent formation water resistivity spectra and reservoir indexes on the basis of electrical imaging. The energy storage indexes mainly included the rock composition factor, porosity of electrical imaging logging, heterogeneity factor, and fracture factor, etc. The research results show that the electrical imaging logging data can effectively evaluate the physical properties and oil-bearing properties of conglo-merate reservoirs, with good application effects in evaluation of field exploration. It has provided a reference for further application of electrical imaging logging data in conglomerate reservoirs.
-
随着超深井、定向井、水平井的逐渐增多,复杂载荷作用下套管的应力计算问题越来越受到关注。目前的研究主要针对居中套管(包括偏磨)在均匀和非均匀地应力条件下的强度计算[1–9],但在非均匀地应力条件下无法求得不居中套管应力的解析解,W. J. Rodriguez[10],P. D. Pattillo[11]和A. Nabipour[12]等人利用有限元方法进行了研究,只有窦益华[13]讨论了解析求解方法,提出了解决此类问题的新思路。在实际工程中,由于所遇情况复杂多变,常常需要针对不同对象进行有限元建模和分析,所需周期长,很难满足快速、实时解决问题的要求。为此,笔者提出了一种基于支持向量机(support vector machine,SVM)的非均匀地应力条件下不居中套管最大von Mises应力的快速预测方法,并对ϕ311.1 mm垂直井眼中的套管进行了分析研究。
1. SVM原理及其Matlab实现
SVM是一种基于统计学习理论的机器学习算法,通过寻求结构风险最小化来实现实际风险的最小化,在有限信息条件下得到最优结果[14],在解决小样本、非线性和高维模式识别问题中表现出许多特有的优势。SVM可以作为一种广义的线性分类器,它的原理是利用非线性变换将输入空间变换到一个高维的特征空间,并在新的空间寻找最优线性分界面。线性可分的情况下,为确保经验风险最小,选取最优分界线H时不仅需要分类准确,还要使分类间隔M(H1和H2之间的距离)最大[14],如图1所示。
线性不可分情况下,SVM的主要思想是将输入向量
x 映射到一个高维的特征向量空间,用特征向量ϕ(x) 来代替,从而得到最优分类函数。由于在运算过程中,无论是目标函数还是决策函数都只涉及到训练样本之间的内积运算,而原空间的核函数K(xi,xj) 可以取代该内积运算,因此可通过选择合适的核函数完成从线性问题到非线性问题的推广[14]。支持向量机回归算法(support vector regression,SVR)是SVM的衍生算法,其本质是需要寻求一个最优超平面,使所有样本点与该超平面的总体偏差最小[15]。也就是说,对于一个容量为n的样本集:
T={(xi,yi),i=1,2,⋯,n},xi∈Rd,yi∈R (1) 式中:
i 为样本序号;d为xi 的维数。需要寻找一个最优回归函数f,使
f(xi) 尽可能接近对应的样本值yi 。通过构造Lagrange函数,将支持向量机回归问题转化为凸二次规划寻优的对偶问题[15]:
{max (2) \!\;\!{\text{其中}}\quad\qquad K\left( {{{{x}}_i}, {{{x}}_j}} \right) = \exp {\left( - \frac{{{{\left\| {{{{x}}_i} - {{{x}}_j}} \right\|}^2}}}{{2{\sigma ^2}}}\right)^{}}\quad\quad (3) 式中:C为惩罚因子,是一个大于0的常数;
a_i^{} ,a_i^ * ,a_j^{} 和a_j^ * 为Lagrange乘子;\varepsilon 为不敏感系数;K({{{x}}_i}, {{{x}}_j}) 为径向基核函数[15];\sigma 为核函数参数。若问题的解为
a_{i{\rm{o}}} 、a_{i{\rm{o}}}^ * ,则相应的回归函数为[15]:f({{x}}) = \sum\limits_{i = 1}^n {(a_{{i{\rm{o}}}}^ * - a_{i{\rm{o}}}^{})} K\left( {{{{x}}_i}, {{{x}}_j}} \right) + {b_{\rm{o}}} (4) 式中:
{b_{\rm{o}}} 为回归参数值。式(4)即为基于径向基核函数的
\varepsilon {\rm{ - SVR}} 模型。\varepsilon {\rm{ - SVR}} 模型的回归效果主要受核函数参数\sigma 、不敏感系数\varepsilon 及惩罚因子C 等参数的影响[16–17]。通过比较各种参数优选的方法,选用网络搜索法进行参数寻优[18],以训练样本均方根误差(root mean squared error,RMSE)最小作为寻优标准,即:\mathop {\min }\limits_{\sigma , C, \varepsilon } \sqrt {\frac{1}{n}\sum\limits_{i = 1}^n {{{\left( {f\left( {{{{x}}_i}} \right) - {y_i}} \right)}^2}} } (5) 式中:
f\left( {{{{x}}_i}} \right) 为预测值。基于Matlab软件,可以实现程序运算,具体运算框图见图2。
2. SVM实验样本构建
为了利用SVM方法实现非均匀地应力和套管不居中情况下套管最大应力的预测,必须建立实验样本。由于实际的“套管–水泥环–地层”系统及其受载十分复杂,这些样本无法用真正的实验方法构建,为此利用ANSYS有限元方法进行“实验样本”的构建。关于“套管–水泥环–地层”系统的有限元分析已十分成熟,也已经过实际数据的验证。这种有限元分析的一个特点是在弹性条件下计算套管的应力,最后根据应力大小和套管的强度来确定套管是否失效[19]。本文的“实验样本”也基于弹性条件进行构建。
影响“套管–水泥环–地层”系统套管应力的因素很多,包括套管的类型和尺寸、最大及最小水平地应力、钻井液密度、水泥环的弹性模量及泊松比、地层的弹性模量及泊松比、套管偏心距和偏心角。关于地应力的影响,有不同的观点,但对于刚钻穿盐膏层的井筒来说,由于应力释放需要较长的时间,因此在一定时间段内,套管将会承受非均匀地应力的作用。同时,研究表明,最大应力一般出现在偏心角为90°或270°的方向[20]。为了减少模型变量,提高预测效率,在构建“实验样本”时,将偏心角固定为90°。同时,当套管给定时,其类型和尺寸都将确定。因此,针对给定套管进行应力预测时,将主要包含最大地应力、最小地应力和水泥环特性等8个关键影响因素。由于实际情况非常复杂,很难确定上述影响因素的合理范围,因此采用了文献调研和实际调查的方法,确定了影响因素的近似取值范围,见表1[21–22]。
表 1 主要影响因素及取值范围Table 1. Main influencing factors and range of values影响因素 取值范围 最大水平主应力σH/MPa 80~135 最小水平主应力σh/MPa 30~80 钻井液密度ρf/(kg∙L–1) 1.15~2.05 水泥环的弹性模量Ec/GPa 10~60 水泥环的泊松比μc 0.15~0.35 地层的弹性模量Es/GPa 1~30 地层的泊松比μs 0.10~0.30 套管偏心距δ/mm 1.5~25.7 某典型的“套管–水泥环–地层”系统如图3所示,井眼直径为311.1 mm,套管内径为219.0 mm,套管外径为250.8 mm。假设套管在井眼中的偏心距为δ,偏心方位角为φ,最大水平主应力σH沿x方向,最小水平主应力σh沿y方向。对几何模型进行单元划分,单元类型选择二维四边形PLANE183单元,按内密外疏的方式划分网格。划分后的网格如图4所示,共计6 215单元,18 953节点。对其施加载荷及约束,偏心方位角固定为90°,分别计算不同偏心距时的套管von Mises应力分布及其最大值。
某一工况下对应的计算参数:偏心距为1.5 mm,σH 为92.0 MPa,σh为59.0 MPa,地层弹性模量和泊松比分别为30 GPa和0.27,水泥环的弹性模量和泊松比分别为35 GPa和0.17,套管的弹性模量和泊松比分别为210 GPa和0.30,钻井液密度为1.64 kg/L,对应井深为4 995.00 m,钻井液静液柱产生的内压约8.03 MPa。地层用以井筒为中心、边长为3.00 m的正方形代替,该工况下的套管内壁von Mises应力云图见图5。从图5可以看出,套管内壁沿y方向(即最小水平主应力方向)的von Mises应力最大达321 MPa。
类似地,可计算得到其他99个“实验样本”数据的最大von Mises应力,因样本数量较大,只给出了1~10号和91~100号样本计算结果(见表2),并将1~95号样本作为训练样本,96~100号样本作为测试样本。为了避免人为干扰,计算过程中,每一个因素的取值都采用随机方法。
表 2 SVM“实验样本”数据Table 2. Data of the SVM “experimental samples”序号 ρf/(kg∙L–1) Ec/GPa μc Es/GPa μs σH/MPa σh/MPa δ/mm σv/MPa 1 1.73 35.00 0.26 15.70 0.25 55.00 107.50 25.7 642.24 2 1.48 38.57 0.18 25.91 0.23 40.71 123.21 1.5 766.62 3 1.48 38.57 0.18 25.91 0.23 40.71 123.21 11.2 768.59 4 1.23 38.57 0.18 17.74 0.31 37.14 80.00 1.5 458.83 5 1.48 20.71 0.32 25.91 0.23 51.43 127.14 20.8 754.53 6 1.48 38.57 0.18 25.91 0.23 40.71 123.21 6.3 767.72 7 1.48 20.71 0.32 25.91 0.23 51.43 127.14 11.2 754.29 8 2.05 42.14 0.35 23.87 0.16 65.71 111.43 1.5 487.51 9 1.73 35.00 0.26 9.57 0.40 51.43 127.14 25.7 1 008.85 10 1.81 27.86 0.17 7.53 0.34 44.29 103.57 1.5 835.23 … … … … … … … … … … 91 1.89 24.29 0.24 21.83 0.10 33.57 99.64 16.0 550.76 92 1.15 56.43 0.29 11.61 0.12 76.43 115.36 25.7 957.12 93 1.73 35.00 0.26 9.57 0.40 51.43 127.14 6.3 1 006.05 94 1.48 38.57 0.18 25.91 0.23 40.71 123.21 16.0 769.33 95 1.23 17.14 0.15 30.00 0.36 65.71 111.43 6.3 468.30 96 1.73 10.00 0.30 21.83 0.10 55.00 107.50 25.7 586.98 97 1.48 60.00 0.21 7.53 0.34 62.14 131.07 16.0 1 192.18 98 1.48 60.00 0.21 7.53 0.34 62.14 131.07 20.8 1 195.77 99 1.40 38.57 0.18 13.66 0.19 37.14 80.00 20.8 531.50 100 1.97 24.29 0.24 5.49 0.27 76.43 115.36 25.7 862.73 3. 套管最大von Mises应力预测及精度分析
基于SVM的非均匀地应力条件下不居中套管最大应力预测,关键在于确定
\varepsilon {\rm{ - SVR}} 模型的3个核心参数,即\sigma 、\varepsilon 和C 的寻优,这需要在一定范围内对3个参数进行搜索。笔者选择\sigma 的范围为[0.1,3.0],搜索步长为0.01;\varepsilon 的范围为[0.01,3.05],搜索步长为0.05;C 的范围为[1,10],搜索步长为0.1。经参数寻优,最佳的\sigma 、\varepsilon 和C 值分别为2.01、0.01和3.00。基于这3个模型参数,预测得到表2中96~100号参数对应的最大应力值(见表3和图6)。为了比较,同时列出了对应的样本值。表 3 测试样本的预测结果Table 3. Predictive effect of test samples样本序号 模型参数 最大von Mises应力/MPa 绝对误差/MPa 相对误差,% 平均相对误差,% 样本值 预测值 96 \sigma =2.01
\varepsilon =0.01
C=3.00586.98 603.21 16.23 2.76 1.32 97 1 192.18 1 200.54 8.36 0.70 98 1 195.77 1 207.74 11.97 1.00 99 531.50 530.85 –0.65 –0.12 100 862.73 845.36 –17.37 –2.01 从图6可以看出,95个训练样本的预测结果和样本值基本重合,表明训练样本的拟合精度很高。5个测试样本的相对误差分别为2.76%、0.70%、1.00%、–0.12%和–2.01%,平均相对误差仅为1.32%。可见,利用SVM方法,可以快速预测非均匀地应力条件下不居中套管的最大von Mises应力。
4. 结论与建议
1)非均匀地应力条件下,目前主要采用有限元分析方法求解不居中套管的应力,缺少数值解,给套管强度校核带来了不便。利用基于SVM的智能预测方法可以解决这一难题,为实际套管应力预测提供了一种新的方法。
2)基于径向基函数的
\varepsilon - {\rm{SVR}} 模型的预测精度较高,测试样本预测值与“实验样本”“真实值”的平均相对误差仅为1.32%。3)在已知影响因素合理范围内,利用有限元分析方法可以获得大量的“实验样本”,不但成本低,而且速度快,为利用SVM方法解决套管应力预测难题提供了有效手段。
4)本文仅以一种套管为例进行了最大von Mises应力SVM预测,对于其他种类套管,需要重新进行“实验样本”计算和
\varepsilon - {\rm{SVR}} 模型关键参数寻优。在保证计算获得的“实验样本”准确可靠的前提下,可以采用该方法预测各种套管的最大von Mises应力。 -
-
[1] 袁阳,侯国强,马涛涛,等. 一种基于测井电成像图像的砂砾岩岩性自动识别方法[J]. 测控技术,2021,40(7):30–35. YUAN Yang, HOU Guoqiang, MA Taotao, et al. An automatic lithology recognition method of sand conglomerate based on electric imaging logging[J]. Measurement & Control Technology, 2021, 40(7): 30–35.
[2] 蔺建华,何雪莹,董丽新. 致密砂砾岩储层岩性识别技术应用[J]. 石油管材与仪器,2021,7(4):80–83. LIN Jianhua, HE Xueying, DONG Lixin. The application of the complex lithology identification technology on tight glutenite formation[J]. Petroleum Tubular Goods & Instruments, 2021, 7(4): 80–83.
[3] 罗兴平,庞旭,苏东旭,等. 电成像测井在复杂砂砾岩储集层岩性识别中的应用:以准噶尔盆地玛湖凹陷西斜坡百口泉组为例[J]. 新疆石油地质,2018,39(3):345–351. LUO Xingping, PANG Xu, SU Dongxu, et al. Recognition of complicated sandy conglomerate reservoir based on micro-resistivity imaging logging: a case study of Baikouquan Formation in western slope of Mahu Sag, Junggar Basin[J]. Xinjiang Petroleum Geology, 2018, 39(3): 345–351.
[4] 李茂兵. 电成像测井自动识别和定量评价研究[D]. 青岛: 中国石油大学(华东), 2010. LI Maobing. Automatic recognition and quantitative evaluation study on electric image logging[D]. Qingdao: China University of Petroleum(East China), 2010.
[5] 张涛. BP神经网络在测井解释中的应用研究[D]. 西安: 西北大学, 2010. ZHANG Tao. Application of BP neural network in logging interpretation[D]. Xi’an: Northwest University, 2010.
[6] 孙健,周魁,冉小丰,等. Bayes判别分析方法在岩性识别中的应用[J]. 石油天然气学报,2009,31(2):74–77. doi: 10.3969/j.issn.1000-9752.2009.02.017 SUN Jian, ZHOU Kui, RAN Xiaofeng, et al. Bayes discriminant analysis method in lithology recognition[J]. Journal of Oil and Gas Technology, 2009, 31(2): 74–77. doi: 10.3969/j.issn.1000-9752.2009.02.017
[7] 寻知锋,余继峰. 聚类和判别分析在测井岩性识别中的应用[J]. 山东科技大学学报(自然科学版),2008,27(5):10–13. XUN Zhifeng, YU Jifeng. The application of cluster and discriminant analyses in logging lithology recognition[J]. Journal of Shandong University of Science and Technology (Natural Science), 2008, 27(5): 10–13.
[8] 张治国. 人工神经网络及其在地学中的应用研究[D]. 长春: 吉林大学, 2006. ZHANG Zhiguo. Study on artificial neural networks and their applications in geoscience[D]. Changchun: Jilin University, 2006.
[9] 李灿元. 成像测井的图像自适应识别方法[D]. 哈尔滨: 哈尔滨工业大学, 2003. LI Canyuan. An adaptive image recognition method for imaging logging[D]. Harbin: Harbin Institute of Technology, 2003.
[10] 庄东志,周凤鸣,司兆伟,等. 南堡5号构造深层火山岩岩性测井识别技术[J]. 断块油气田,2011,18(6):743–745. ZHUANG Dongzhi, ZHOU Fengming, SI Zhaowei, et al. Logging identification of deep volcanic lithology in No.5 structure of Nanpu Sag [J]. Fault-Block Oil & Gas Field, 2011, 18(6): 743–745.
[11] 周云才, 刘瑞林. 图像识别方法在FMI图像处理中的应用[J]. 江汉石油学院学报, 2003, 25(增刊1): 50–51. ZHOU Yuncai, LIU Ruilin. Application of image recognition for FMI image processing[J]. Journal of Jianghan Petroleum Institute, 2003, 25(supplement 1): 50–51.
[12] 周萍. 王府断陷火石岭组火山岩岩性及岩相识别[J]. 断块油气田,2020,27(2):188–192. ZHOU Ping. Identification of volcanics reservoir lithology and lithofacies in Huoshiling Formation of Wangfu Fault Depression [J]. Fault-Block Oil & Gas Field, 2020, 27(2): 188–192.
[13] AKBAR M, CHAKRAVORTY S, RUSSELL S D, et al. Unconventional approach to resolving primary and secondary porosity in gulf carbonates from conventional logs and borehole images[R]. SPE 87297, 2000.
[14] 刘丹,潘保芝,房春慧,等. 利用成像测井评价储层孔隙空间的二维非均质性[J]. 世界地质,2014,33(3):640–646. doi: 10.3969/j.issn.1004-5589.2014.03.015 LIU Dan, PAN Baozhi, FANG Chunhui, et al. Evaluation of 2D heterogeneity of reservoir’s pore space based on imaging logging[J]. Global Geology, 2014, 33(3): 640–646. doi: 10.3969/j.issn.1004-5589.2014.03.015
[15] 侯振学,于雪娴,李东旭,等. 电成像测井处理新技术在储层评价方面的应用[J]. 地球物理学进展,2020,35(2):573–578. doi: 10.6038/pg2020DD0047 HOU Zhenxue, YU Xuexian, LI Dongxu, et al. Application of new processing technology of electrical imaging logging in reservoir evaluation[J]. Progress in Geophysics, 2020, 35(2): 573–578. doi: 10.6038/pg2020DD0047
[16] 梁劳劳. 致密砂岩储层裂缝最优识别方法及软件开发[D]. 成都: 成都理工大学, 2012. LIANG Laolao. Optimal identification method and software development of cracks in tight sandstone reservoir[D]. Chengdu: Chengdu University of Technology, 2012.
[17] 赫志兵. 辽河坳陷兴马潜山储层测井解释方法研究[D]. 北京: 中国地质大学(北京), 2012. HE Zhibing. Logging interpretation method of Xing-Ma buried hill reservoir in Liaohe Depression[D]. Beijing: China University of Geosciences(Beijing), 2012.
[18] 吴斌,唐洪,王素荣,等. 利用蚁群算法识别及预测碳酸盐岩裂缝的方法探讨[J]. 重庆大学学报(自然科学版),2012,35(9):131–138. WU Bin, TANG Hong, WANG Surong, et al. Recognition and prediction colony algorithm analysis of fracture by using ant in carbonate reservoir[J]. Journal of Chongqing University(Natural Science Edition), 2012, 35(9): 131–138.
[19] 石广仁. 支持向量机在裂缝预测及含气性评价应用中的优越性[J]. 石油勘探与开发,2008,35(5):588–594. doi: 10.3321/j.issn:1000-0747.2008.05.010 SHI Guangren. Superiorities of support vector machine in fracture prediction and gassiness evaluation[J]. Petroleum Exploration and Development, 2008, 35(5): 588–594. doi: 10.3321/j.issn:1000-0747.2008.05.010
[20] 刘伟,张晋言,张文姣,等. 基于电成像测井资料的砂砾岩储层有效性分类评价方法[J]. 石油钻探技术,2016,44(4):114–119. LIU Wei, ZHANG Jinyan, ZHANG Wenjiao, et al. An evaluation method for glutenite reservoir effectiveness classification based on electrical imaging logging data[J]. Petroleum Drilling Techniques, 2016, 44(4): 114–119.
[21] LAI Jin, PANG Xiaojiao, XIAO Qiyao, et al. Prediction of reservoir quality in carbonates via porosity spectrum from image logs[J]. Journal of Petroleum Science and Engineering, 2019, 173: 197–208. doi: 10.1016/j.petrol.2018.10.022
[22] 李宁,肖承文,伍丽红,等. 复杂碳酸盐岩储层测井评价:中国的创新与发展[J]. 测井技术,2014,38(1):1–10. LI Ning, XIAO Chengwen, WU Lihong, et al. The innovation and development of log evaluation for complex carbonate reservoir in China[J]. Well Logging Technology, 2014, 38(1): 1–10.
-
期刊类型引用(4)
1. 苗娟,黄兵,谢力,汤明. 高钢级套管段铣工具优化及性能评价. 特种油气藏. 2021(02): 163-170 . 百度学术
2. 席岩,李军,陶谦,曾夏茂,钟文力. 循环载荷作用下微环隙的产生及演变. 断块油气田. 2020(04): 522-527 . 百度学术
3. 冯翔宇,石茂林,马跃,宋学官,孙伟,田腾. 基于支持向量回归的定性-定量因子混合建模方法. 大连理工大学学报. 2020(06): 599-609 . 百度学术
4. 孙挺,赵颖,杨进,殷启帅,汪文星,陈缘. 基于支持向量机的钻井工况实时智能识别方法. 石油钻探技术. 2019(05): 28-33 . 本站查看
其他类型引用(1)