塔里木盆地西部A区块油气资源丰富，油气藏具有超深、超高压等特点，储层平均埋藏深度超过7 300 m，勘探开发难度大。实钻资料证实，该区块钻井过程中发生过漏失的井占比69.5%，是对井漏进行大数据分析与钻前诊断研究的理想区块。井漏不仅会使钻井成本增加、建井周期增长，而且会诱发缩径、井壁失稳和地层破坏等一系列井下故障或复杂情况^[1-3]。

近年来，井漏影响因素分析为钻前漏失诊断提供了大量思路：S. Al Menhali等人^[4]认为，易发生井漏的主要是高渗透性弱固结地层、裂缝性碳酸盐岩储层和低含水地层；A. T. Al Hameedi等人^[5]利用机器学习算法分析了钻井参数对漏失速度的敏感性，回归了包含塑性黏度、屈服强度、当量循环密度、钻井液相对密度、喷嘴总流动面积、钻压、钻速等7个参数的漏失速度计算公式；H. B. Wibowo等人^[6]进行了包括地层孔隙压力、地层破裂压力、最小水平应力、最大水平应力、垂直应力、上覆岩层压力和剪切破坏梯度的地下压力分析，将井筒压力大于地层最小水平应力作为井漏的必要条件。

塔里木盆地西部A区块内走滑断裂带发育规模大，各类天然裂缝十分发育且分布复杂^[7]，且断裂带附近地层承压能力低，容易发生漏失。为此，笔者针对该区块钻井过程中出现的漏失问题，提出利用地层破裂压力与地层孔隙压力当量密度的重合程度来表征地层裂缝发育程度，并在此基础上建立了结合裂缝发育程度、裂缝宽度的井漏风险系数计算方法，综合分析了A区块易漏地层的裂缝发育程度、裂缝宽度不确定性范围及井漏风险系数，实现了钻前井漏风险诊断，为该区块钻井过程中防漏技术措施的制定提供了依据，为易漏地层提供了风险提示。

1.   A区块地质概况

1.1   地层构造特征

A区块位于塔里木盆地东北坳陷区西部，横跨阿瓦提断陷的北斜坡、顺托果勒低隆起西北部与沙雅隆起的西南部倾末端，自下而上可划分为奥陶系—志留系、泥盆系—二叠系、中生界、新生界等4套构造层。奥陶系地层发育齐全，自下而上发育有O_1-2y，O₂yj，O₃q，O₃l和O₃s等5个地层；二叠系仅发育P₂层，志留系有S₁t和S₁k等2个地层；O_1-2y和O₂yj发育碳酸盐岩缝洞型断溶体油藏^[8]，是主要的目的层段。

1.2   断裂体系特征

A区块内断裂体系主要由多组北西向、北东向、近南北向断穿基底的走滑断裂构成，断裂发育规模大，断穿层位多^[9]。断裂带及断裂带附近的碳酸盐岩受应力作用改造破碎，岩体周围以盖层滑脱为主，平面上断裂延伸距离较短、垂直断距较小，地层变形较弱，后期活动弱^[10-12]。区块内分布有18条断裂带，依据断裂等级可分为主干一级断裂（8条）、主干二级断裂（10条）2种^[13]，其中，主干一级断裂基底断裂特征明显，变形强，断距大；主干二级断裂基底断裂特征不明显，变形较弱^[9]。

1.3   井漏分析

实际钻井资料证实，2015年以来，A区块北东向走滑断裂主干一级断裂带上的13口井和主干二级断裂带上的12口井在钻井过程中均发生过井漏，在地震剖面上呈“断裂+串珠”反射特征。对A区块各地层发生井漏的次数和占比进行了统计，结果如图1所示。

图  1  A区块各地层井漏次数统计

Figure  1.  Statistics of the lost circulation of each layer in Block A

下载: 全尺寸图片 幻灯片

从图1可以看出，P₂，S₁t，S₁k，O₃s，O₂yj和O_1-2y等6个地层的漏失次数占比高达93.66%。因此，选取这6个地层作为A区块井漏分析与诊断的研究对象。

2.   基于自然语言处理的钻井漏失信息挖掘

自然语言处理方法（natural language processing，NLP）^[14]利用词法分析器和解析器，将文本转换为标记，并用解析器验证标记序列。文本信息检索方法（IR）^[15]采用高度可扩展统计方法搜索大量文本，而NLP的分析器是将正则表达式和Backus-Naur规范作为输入，生成确定词法或解析决策的代码和查找表，方便通过编程语言实现。基于自然语言处理算法的井漏信息挖掘技术，可以快速从各类钻井报表中提取与井漏相关的关键信息。首先分析各类钻井报表，总结其中涉及到井漏描述以及漏失的位置。首先，基于正则表达式筛选井漏信息，得到关键信息并提取展示；然后，将离散信息进行聚合关联，设计数据库表结构，为构建大规模的井漏案例库奠定基础。其主要逻辑顺序如下：

1）对各类钻井报表进行差异化分析和关键词总结，提取各类钻井报表中的差异化部分和“锚点”（见图2）。“锚点”指的是钻井报表中可作为定位关键信息的内容。

图  2  信息提取

Figure  2.  Information extraction

下载: 全尺寸图片 幻灯片

2）根据提取内容，设计基于Python的提取模板，再将提取的内容根据正则表达式进行过滤筛选。

3）将收集的离散信息进行聚合关联，如提取到某天某时间段发生什么故障或井下复杂情况，并匹配出当时的钻具组合和钻井液等信息。

4）将提取出来的信息存储到漏失数据库中，做永久化存储。漏失数据库主要包含静态数据、动态数据和结论数据。静态数据包括案例编号、井深、区块、井漏时作业工况、漏失层岩性、漏失层裂缝发育情况、漏失速度、漏失类型、钻井液性能、漏失通道尺寸、孔隙压力、破裂压力；动态数据包括钻井液池中钻井液的体积、出口流量、入口流量、立管压力、大钩载荷、扭矩、钻速、转速、泵压；结论数据包括井漏原因和处理措施。

3.   基于大数据分析的漏失分析与诊断

3.1   井漏分析

3.1.1   漏失数据处理

针对A区块P₂，S₁t，S₁k，O₃s，O₂yj和O_1-2y等6个易漏地层，利用大数据分析技术，汇总分析了该区块59口完钻井的井漏信息，包括了漏失次数、漏失率、钻井液密度、平均漏失量、平均损失周期、最大漏失速度、平均漏失速度，结果见表1。由表1可知，P₂，S₁t，S₁k，O₃s，O₂yj和O_1-2y等6个易漏地层都属于平均漏失量大于100 m³的大型漏失，漏失速度也都属于10～30 m³/h的中速漏失。

表  1  A区块井漏信息统计结果

Table  1.  Information statistics of lost circulation in Block A

地层	漏失次数	漏失率，%	钻井液密度/（kg·L−1）	平均漏失量/m3	平均损失周期/h	最大漏失速度/（m3·h−1）	平均漏失速度/（m3·h−1）
P2	48	37.29	1.22～1.42	176.942	90.637	27.115	25.818
S1t	62	20.34	1.27～1.40	115.674	69.352	27.098	14.952
S1k	31	18.64	1.27～1.44	221.714	237.366	25.204	18.590
O3s	32	20.34	1.29～1.87	136.866	37.015	12.352	9.411
O2yj	27	35.59	1.16～1.72	320.344	104.641	19.975	14.558
O1-2y	22	33.90	1.10～2.00	296.999	129.677	9.162	11.416

下载: 导出CSV 

| 显示表格

P₂为火成岩地层，岩性主要为英安岩、凝灰岩、玄武岩，岩性硬、脆，在构造运动和风化作用下形成了气孔、孔隙、裂缝，构成了主要漏失通道，局部裂缝、缝洞发育，地层承压能力低^[16]，该地层共漏失48次，漏失率高达37.29%，最大漏失速度和平均漏失速度也都位列第一。S₁t和S₁k含有断层和裂缝，地层薄弱，承压能力低，其中S₁t漏失了62次，漏失率为20.34%，最大漏失速度位列第二，平均漏失速度位列第三；S₁k漏失了31次，漏失率为18.64%，最大漏失速度位列第三，平均漏失速度位列第二，平均损失周期最长。A区块由于受到的应力相对较小，主要断裂发育在O₃s之下，产状近直立，平面呈“X”分支断裂数量较多^[17]，地层变形较弱，断距较小，后期活动弱，O₃s漏失了32次，漏失率为20.34%，平均漏失速度最小；O₂yj漏失了27次，漏失率为35.59%，平均漏失量最大；O_1-2y漏失了22次，漏失率为33.90%，平均漏失量和平均损失周期都位列第二。

3.1.2   漏失影响因素分析

利用大数据分析技术，统计了A区块59口完钻井6个易漏地层的钻井数据、钻井液数据、井漏数据和地层数据。钻井数据包括事故编号、井号、井深、开次、钻头直径；钻井液数据包括其密度、漏斗黏度、屈服应力、流性指数、稠度系数；井漏数据包括漏失开始时间、漏失结束时间、漏失量、漏失速度、漏失时间、漏失类型、漏失原因；地层数据包括地层代号、实际井筒压力、孔隙度、地层压力和破裂压力。地层数据中的实际井筒压力是根据钻井日报中的钻井液密度计算得到的，预测地层压力和预测破裂压力是根据钻前工程设计中的预测地层压力当量密度和预测破裂压力当量密度计算得到的，实际地层压力和实际破裂压力是根据钻井地层破裂压力试验和地层压漏试验得到的。

提取井漏信息中的压力信息，分析破裂压力、地层压力、井筒压力与每次井漏中漏失量的相关性，结果如图3所示（图3中，一个红色半透明圆点代表一次井漏的漏失量，圆点越大，漏失量越大）。从图3可以看出：井筒压力高于地层压力时，由于存在天然裂缝，容易出现漏失量较小的天然裂缝性井漏；井筒压力位于破裂压力与地层压力重合范围时，容易出现漏失量中等的渗透性井漏；井筒压力高于破裂压力时，容易出现漏失量较大的压裂井漏。将破裂压力与地层压力的当量密度重合程度类比为裂缝发育程度，其与裂缝宽度及井筒压力与地层压力、破裂压力之间的差值确定为井漏风险系数的3大影响因素。

图  3  破裂压力、地层压力、井筒压力与漏失量的相关性

Figure  3.  Correlation among fracture pressure, formation pressure, wellbore pressure, and leakage

下载: 全尺寸图片 幻灯片

3.2   裂缝发育程度分析

3.2.1   实际地层压力当量密度分布情况

地层压力是地层孔隙中地层水所产生的液柱压力，A区块的实际地层压力大部分是在钻井过程中通过dc指数法利用机械钻速与压差的负相关原理检测得到，其余则是利用地层破裂压力试验、地层压漏试验、当量密度预测法、关井套压折算法、试油静压、射孔测试等得到，当量密度预测法多用在O₂yj和O_1-2y地层，关井套压折算法多用在O_1-2y地层。

基于大数据分析理论，分析了6个易漏地层钻井过程中监测到的实际地层压力数据，并将其转换为地层压力当量密度，以消除井深的影响。以位于1号断裂带上的23口井为例（见图4），6个易漏地层压力当量密度不确定性分布范围的重合度较大。

图  4  地层压力当量密度分布

Figure  4.  Equivalent density distribution of the formation pressure

下载: 全尺寸图片 幻灯片

3.2.2   实际破裂压力当量密度分布情况

破裂压力是使地层产生水力裂缝或使原有裂缝张开时的井底压力。利用地层压漏试验测量A区块的实际破裂压力。一般在二开、三开钻进套管鞋以下第1个砂层3～5 m时进行地层压漏试验，将不再随注入量增大而增大或随注入量增大而减小时的压力作为地层破裂压力，或将达到井口承压设备中的最小额定工作压力或套管承受压力达到套管最小抗压强度80%时的压力作为地层破裂压力。

基于大数据分析理论，分析了A区块59口井P₂，S₁t，S₁k，O₃s，O₂yj和O_1-2y等6个易漏地层在钻井过程中监测到的实际破裂压力，并将其转换为破裂压力当量密度，以消除井深的影响。以位于1号断裂带上的23口井为例，综合分析6个易漏地层，发现破裂压力当量密度不确定性的重合度较大（见图5）。由于大部分已完钻井只进行了2次地层压漏试验，实际破裂压力当量密度数据较少，不确定性分布范围对比地层压力当量密度较差。

图  5  破裂压力当量密度分布

Figure  5.  Equivalent density distribution of the fracture pressure

下载: 全尺寸图片 幻灯片

3.2.3   压力当量密度不确定性分布范围分析

由于地层破裂压力试验和地层压漏试验次数较少，且无法得到井漏时的地层孔隙压力和破裂压力，所以汇总各个地层的实际地层压力当量密度和破裂压力当量密度进行分布范围分析。为了消除超出物理边界的值，将其拟合为正态分布，将正态分布中95.449 974%（μ−2δ～μ+2δ）的概率范围类比为各地层实际地层压力和破裂压力当量密度的不确定性分布范围，以1号断裂带上23口井为例进行统计分析，将各地层的实际压力当量密度先分箱处理，再计算每个分箱的分布频率，结果见图6、图7。图6、图7中，蓝色曲线为高斯核密度估计曲线，横轴上方蓝色小细条为实际压力当量密度数据分布，橙色曲线为正态分布曲线，可直观地评估其与实际压力当量密度数据的对应关系。

图  6  各地层压力当量密度不确定性分布

Figure  6.  Uncertainty distribution of equivalent density of the formation pressure in each layer

下载: 全尺寸图片 幻灯片

图  7  各地层破裂压力当量密度不确定性分布

Figure  7.  Uncertainty distribution of equivalent density of the fracture pressure in each layer

下载: 全尺寸图片 幻灯片

表2和表3为1号断裂带上23口井6个易漏地层实际地层压力、破裂压力当量密度不确定性分布统计结果。由表2和表3可知：地层埋藏深度越大，地层压力当量密度均值越小，其中S₁t的数据较少导致均值偏小；破裂压力当量密度的均值和离散程度呈现先增大后减小的趋势，均值峰值在厚度最大的O₃s，离散程度峰值在厚度较小、深度较大的O₂yj；除O₂yj外，地层压力当量密度的离散程度普遍比破裂压力当量密度高，破裂压力当量密度由于数据较少导致其不确定分布幅度小，分布较集中。

表  2  各地层实际地层压力当量密度分布统计结果

Table  2.  Statistics of actual equivalent density distribution of the formation pressure in each layer

地层	均值/（kg·L−1）	标准差/（kg·L−1）	分布范围下限/（kg·L−1）	分布范围上限/（kg·L−1）
P2	1.311692221	0.321 164 305	0.669 363 610	1.954 020 831
S1t	1.151 520 157	0.265 860 717	0.619 798 723	1.683 241 592
S1k	1.279 467 735	0.181 530 373	0.916 406 989	1.642 528 481
O3s	1.224 136 368	0.288 242 768	0.647 650 831	1.800 621 904
O2yj	1.126 302 960	0.128 635 852	0.869 031 257	1.383 574 664
O1-2y	1.076 552 962	0.249 590 264	0.577 372 434	1.575 733 490

下载: 导出CSV 

| 显示表格

表  3  各地层实际破裂压力当量密度分布统计结果

Table  3.  Statistics of actual equivalent density distribution of the fracture pressure in each layer

地层	均值/（kg·L−1）	标准差/（kg·L−1）	分布范围下限/（kg·L−1）	分布范围上限/（kg·L−1）
P2	1.762 290 617	0.083 384 714	1.595 521 189	1.929 060 044
S1t	1.763 333 333	0.099 247 166	1.564 839 001	1.961 827 666
S1k	1.780 773 826	0.127 570 246	1.525 633 335	2.035 914 317
O3s	1.952 118 014	0.220 338 077	1.511 441 860	2.392 794 168
O2yj	1.900 000 000	0.241 660 919	1.416 678 161	2.383 321 839
O1-2y	1.724 581 876	0.115 840 777	1.492 900 323	1.956 263 429

下载: 导出CSV 

| 显示表格

3.2.4   裂缝发育程度分析

预先了解区块内同一断裂带的天然裂缝发育特性，可以提高漏失量的预测准确度和目的储层的整体生产能力^[18–20]。通过分析漏失影响因素，提出利用地层压力与破裂压力当量密度不确定性分布范围的重合程度来表征地层的裂缝发育程度。通过计算各地层的裂缝发育系数来表征地层的裂缝发育程度：分别将某一地层的地层压力和破裂压力的当量密度范围映射到（−5，5）区间内，利用双曲正切函数将裂缝发育程度系数集中在（0，1）区间内，则裂缝发育系数为：

式中：Y为裂缝发育系数；${F}_{{\rm{pu}}}$为某一地层内的地层压力当量密度不确定性分布范围上限，kg/L；${F}_{{\rm{fd}}}$为某一地层内的破裂压力当量密度不确定性分布范围下限，kg/L。

以1号断裂带上的23口井为例，将表2中的地层压力当量密度和表3中的破裂压力当量密度不确定性分布范围代入式（1），计算得6个易漏地层的裂缝发育系数，发现：O₃s的裂缝发育程度最高，发育系数为0.486 475；O₂yj的裂缝发育程度最低，发育系数为−0.020 127；其他4个地层的裂缝发育程度由低到高依次为P₂（0.011 860）、S₁t（0.013 689）、O₁₋₂y（0.025 093）和S₁k（0.243 685）。

3.3   裂缝宽度分析

自然裂缝大小的估计是否可靠对井漏处理至关重要^[21]，也可以用其表征自然裂缝储集层，为完井设计提供依据。R. Majidi等人^[22–23]将裂缝中钻井液的损失假设为2个平行圆盘之间流体的径向层流流动，建立了基于屈服–幂律流变学理论的模型。利用裂缝中瞬态径向流的质量和线性动量的守恒原理，假设裂缝宽度恒定，利用数值方法可以求得该模型闭合形式的近似解，获得流动方向上任意2点之间压降的闭合控制方程：

式中：K为钻井液的稠度系数，Pa·sⁿ；q为漏失速度，m³/h；n为钻井液的流性指数；w为地层裂缝的宽度，m；r为裂缝到井的径向距离，m；${\tau }_{y}$为钻井液的屈服应力，Pa。

式（2）等号右侧的第1项是由流体的剪切减薄/增稠效应产生，径向发散流的速度降低会引起钻井液（非牛顿流体）的剪切速率和表观黏度显著变化，特别是在井漏初期剪切变稀大大提高了压降损失率。式（2）等号右侧的第二项由流体的屈服应力产生，表明具有屈服应力的流体只有压力梯度达到特定值时才能流动，而最终压降损失也由流体屈服应力控制。

由于A区块大多深井存在裂缝，且多为垂直缝或斜裂缝，所以采用考虑裂缝宽度、深度、长度的漏失控制方程^[24]。假设流体损失期间过平衡压力恒定，施加在整个钻井液液柱中，从井筒半径${r}_{\mathrm{w}}$积分到裂缝中钻井液侵入深度${r}_{\mathrm{f}}$，整个钻井液液柱的局部压降为：

式中：p_w为井筒压力，MPa；p_p为钻井液入侵前沿压力（地层压力），MPa；${r}_{\mathrm{f}}$为钻井液侵入前沿的径向长度，m；${r}_{\mathrm{w}}$为井筒半径，m；h为裂缝长度，m。

以1号断裂带上的23口井为例，利用式（3）定量分析这23口井的井漏信息，计算得到6个易漏地层的裂缝宽度分布情况，如图8所示。计算分析用到的漏失速度为井漏初期且未泵入堵漏浆时的瞬时漏失速度^[25]，地层压力为图7中的地层压力不确定性分布范围。由图8可知，6个易漏地层的裂缝宽度都比较集中，且都有明显离群点出现，S₁t地层的离群点最少，O_1-2y地层的离群点最多，S₁t，O₃s和O₂yj等3个地层的裂缝宽度分布较为集中。

图  8  裂缝宽度分布

Figure  8.  Fracture width distribution

下载: 全尺寸图片 幻灯片

汇总各地层的所有裂缝宽度^[26]，进行分布范围分析，并将其进行正态分布拟合，各地层裂缝宽度的不确定性分布范围统计结果如图9、表4所示。基于裂缝宽度的物理性质，将正态分布的起始裂缝宽度值都设定为0。6个易漏地层裂缝宽度的最大值和上限值相差较大，其他参数均相差不多。其中O₂yj裂缝宽度的最大值最大，O_1-2y的最大值最小；P₂裂缝宽度的上限值最大，O_1-2y的上限值最小。

图  9  裂缝宽度不确定性分布

Figure  9.  Uncertainty distribution of the fracture width

下载: 全尺寸图片 幻灯片

表  4  裂缝宽度不确定性分布统计

Table  4.  Statistics of uncertainty distribution of fracture width

地层	最大裂缝宽度/m	最小裂缝宽度/m	裂缝宽度均值/m	裂缝宽度标准差/m	裂缝宽度上限95%/m
P2	0.058 90	0.001 6	0.016 70	0.012 70	0.042 00
S1t	0.024 30	0.000 8	0.006 70	0.003 90	0.014 50
S1k	0.020 60	0.001 3	0.009 10	0.006 80	0.022 70
O3s	0.057 90	0.000 4	0.009 80	0.001 09	0.031 60
O2yj	0.131 60	0.001 1	0.008 30	0.013 20	0.034 70
O1-2y	0.005 71	0.001 2	0.001 86	0.014 50	0.004 75

下载: 导出CSV 

| 显示表格

3.4   井漏风险诊断

井漏的风险系数受井筒压力与地层压力及破裂压力的差值、裂缝宽度及裂缝发育程度的影响，结合这3个影响因素的不确定分布范围，提出了井漏风险系数不确定分布范围的计算公式：

式中：K为井漏风险系数。

针对所有记录了泵入堵漏浆前瞬时漏失速度的井漏，利用式（4）计算得到井漏风险系数，并分析井漏风险系数与瞬时漏失速度的相关性，结果如图10所示。

图  10  井漏风险系数与瞬时漏失速度的相关性

Figure  10.  Correlation analysis of lost circulation risk coefficient and instantaneous leakage velocity

下载: 全尺寸图片 幻灯片

从图10可以看出，井漏风险系数与瞬时漏失速度具有统计学上的关联性，符合正相关关系，表明利用井漏风险系数可以较好地评价井漏的严重程度，实现钻前漏失诊断。

3.5   钻前井漏风险诊断实例分析

以1号断裂带A区块的8井为例进行钻前井漏风险诊断。由8井的《钻井工程设计》得知：该井位于塔里木盆地顺托果勒低隆北缘，目的层位为奥陶系中—下统鹰山组（O_1-2y），提取地层孔隙压力、破裂压力的预测当量密度和设计钻井液密度^[27]（见表5），预测该井地层的孔隙压力破裂压力（见表6）。

表  5  A区块8井的预测当量密度

Table  5.  Predicted equivalent density for eight wells in Block A

层位	井深/m	预测地层孔隙压力当量 密度/（kg·L−1）	预测地层破裂压力当量 密度/（kg·L−1）	设计钻井液密度/（kg·L−1）
P2	4 507～4 991	1.17～1.20	1.89～1.96	1.22～1.26
S1t	5 393～6 161	1.20～1.24	1.89～1.96	1.26～1.30
S1k	6 161～6 541	1.20～1.24	1.89～1.96	1.26～1.30
O3s	6 541～7 415	1.20～1.24	1.89～1.96	1.26～1.30
O2yj	7 450～7 572	1.15～1.18	1.65～1.85	1.21～1.29
O1-2y	7 572～7 814	1.15～1.18	1.65～1.85	1.21～1.29

下载: 导出CSV 

| 显示表格

表  6  A区块8井的预测压力

Table  6.  Predicted pressure for eight wells in Block A

层位	井深/m	预测地层孔隙压力/MPa	预测地层破裂压力/MPa	设计井筒压力/MPa
P2	4 507～4 991	51.677～58.694	83.479～95.867	53.886～61.628
S1t	5 393～6 161	63.422～74.868	99.889～118.340	66.592～78.491
S1k	6 161～6 541	72.453～79.486	114.114～125.639	76.076～83.332
O3s	6 541～7 415	76.922～90.107	121.152～142.427	80.768～94.467
O2yj	7 450～7 572	83.962～87.563	120.466～137.280	88.342～95.725
O1-2y	7 572～7 814	85.336～90.361	122.439～141.668	89.789～98.785

下载: 导出CSV 

| 显示表格

8井预测压力与实测压力不确定性的分布情况如图11所示（图11中，绿色填充区域为实测地层压力分布范围，粉色填充区域为实测破裂压力分布范围，粉色与绿色中间重合的灰绿色为破裂压力和地层压力的重合范围）。从图11可以看出：预测地层孔隙压力都在实测地层压力分布范围内；S₁t的预测破裂压力比实测破裂压力高，P₂和S₁k的预测破裂压力与实测破裂压力接近。

图  11  预测压力与实测压力不确定性分布范围对比

Figure  11.  Comparison of uncertainty distribution ranges of predicted pressure and measured pressure

下载: 全尺寸图片 幻灯片

将设计井筒压力、预测破裂压力、地层裂缝宽度和地层裂缝发育系数代入式（4），计算P₂，S₁t，S₁k，O₃s，O₂yj和O_1-2y的井漏风险系数和预测漏失速度，结果见表7。从表7可以看出：P₂井漏风险系数为0.40～0.77，最可能发生漏失，预测的漏失速度为43.27～82.98 m³/s；其他5个地层发生漏失的概率都较小，按井漏风险系数和预测漏失速度从大到小依次为O₃s，S₁k，O_1-2y，S₁t和O₂yj。

表  7  钻前井漏风险诊断结果

Table  7.  Pre-drilling lost circulation risk diagnosis

地层	厚度/m	井漏风险系数	预测漏失速度/（m3·s−1）
P2	484	0.245 0～0.475 8	26.464～51.390
S1t	768	0.000 2～0.000 3	0.020～0.029
S1k	380	0.004 3～0.004 4	0.463～0.470
O3s	874	0.011 5	1.242～1.246
O2yj	122	0～0.000 1	0.006～0.012
O1-2y	240	0.001 0～0.001 1	0.114

下载: 导出CSV 

| 显示表格

4.   结　论

1）针对塔里木盆地A区块井漏问题，利用自然语言处理技术，挖掘了A区块所有完钻井的井漏信息，提出用地层破裂压力与地层压力不确定范围的重合程度来表征地层裂缝发育程度的方法。

2）利用大数据分析技术，分析了A区块易漏地层实际地层压力和实际破裂压力的不确定性分布情况，通过处理漏失数据和分析影响漏失的因素，将裂缝发育程度、裂缝宽度和井筒压力与破裂压力之间的差值确定为影响井漏风险系数的3大因素。

3）漏失影响因素的不确定性分析结果表明：井漏风险系数与瞬时漏失速度具有统计学上的关联性，符合正相关关系；利用井漏风险系数可以较好地评价井漏的严重程度，实现钻前漏失诊断。