1.   研究区概况

顺北地区位于塔里木盆地中部，构造应力复杂，断裂样式多样(图 1)，目的层走滑断裂有利于断控储层发育，但非目的层段走滑断裂及其破碎带不利于井壁稳定，影响高效成井^[1-8]。志留系地层岩性复杂，主要包括依木干他乌组、塔塔埃尔塔格组和柯坪塔格组。其中下志留统柯坪塔格组根据岩性、电性特征分析，自上而下分为3个岩性段，柯上段、柯中段和柯下段。志留系与上覆泥盆系平行不整合接触，与下覆奥陶系呈角度不整合接触(表 1)^[9]。顺北油田钻井钻至上古生界复杂断裂断距较大、裂缝发育部位时，易发生严重漏失。研究区内SHB52井志留系严重漏失2 500 m³，无法继续钻井，造成井口移位。为进一步降低漏失风险，亟需开展研究区断裂空间发育特征及主应力分析，预测断裂发育强度，评价裂缝发育规律，优选非目的层段裂缝发育较弱区域，部署设计井位，降低志留系严重漏失风险，为高效成井提供条件。

图  1  塔里木盆地顺北地区构造位置及断裂分布

Figure  1.  Structural location and fracture distribution of Shunbei area in Tarim Basin

下载: 全尺寸图片 幻灯片

表  1  塔中顺北地区志留系相关地层格架及地震波组综合划分

Table  1.  Comprehensive division of Silurian-related stratigraphic framework and seismic wave groups in Shunbei area, central Tarim Basin

下载: 导出CSV 

| 显示表格

2.   断裂解析与发育特征

2.1   断裂增强解释性处理

通过研究区三维地震资料志留系频谱分析，利用带通滤波解释性处理技术，优选走滑雁列断裂敏感频段，进一步进行倾角导向滤波处理，提高走滑雁列断裂成像精度。志留系—奥陶系桑塔木组主频为22 Hz，带宽为6~38 Hz。开展低频段、中频段、高频段带通滤波测试，其中带通滤波参数频段滤去6 Hz以下及22 Hz以上信息，保留10~15 Hz低频信息，倾角导向滤波敏感平滑参数线方向7线、道方向7道、纵向15 ms，基于走滑雁列断裂敏感频段，结合倾角导向滤波解释性处理，志留系断裂断点、断面特征较断裂增强处理前资料更加清晰(图 2)，结合区域已钻井标定，进行了引层综合对比分析及解析，确定了志留系、奥陶系各地层相应地震波组位置，同时结合志留系断裂发育区相轴错断位置，通过断面、断点及断裂发育继承性关系分析对比，进行了志留系断裂精细解释，志留系发育多条断裂，不同位置断裂断至位置存在差异。

图  2  塔中顺北地区志留系突出断面特征的解释性增强处理

a.原始地震资料；b.断裂增强解释性处理；c.精细断裂解释；d.柯坪塔格组顶界面断裂平面展布。

Figure  2.  Enhanced processing of prominent fracture characteristics of Silurian in Shunbei area, central Tarim Basin

下载: 全尺寸图片 幻灯片

2.2   断裂及发育特征

2.2.1   平面特征

不同位置走滑断裂对储层及钻井影响存在差异，走滑断裂利于储层发育^[10-11]，但非目的层位置对钻井而言，走滑断裂破碎带不利于井壁稳定。走滑断裂不同位置产生不同应力场环境。左阶左行式、右阶右行式表现出伸展应力场环境；右阶左行式、左阶右行式表现出挤压应力场环境^[12-13]。基于地震反射标志层界面志留系顶(T₆⁰)、志留系柯坪塔格组顶(T₆³)、志留系柯坪塔格组内幕(T₆⁵)、志留系底(T₇⁰)及其敏感断裂检测相干属性，志留系断裂右阶雁列式展布规律明显(图 3a)，雁列断裂对应性较好，断裂强度大。

图  3  塔中顺北地区志留系雁列断裂空间精细解释及分布

a.志留系断裂检测相干属性及雁列断裂特征；b.剖面断裂特征；c.“丝带效应”特征；d.志留系断裂空间特征。

Figure  3.  Detailed spatial interpretation and distribution of en echelon fractures in Silurian section, Shunbei area, central Tarim Basin

下载: 全尺寸图片 幻灯片

2.2.2   剖面特征

顺北地区断裂带加里东中期断裂样式主要为走滑挤压段正花状构造、平移段直立断层样式、走滑拉分段下掉正断层负花状样式^[14-15]。结合志留系地震资料断裂响应平面、剖面特征，基于断裂空间闭合解释，主要为拉张伸展应力作用下拉分段的负花状正断层特征(图 3b), 断裂空间特征表现为同一条断裂在不同位置，其走向、倾向会发生微弱变化，具有明显的“丝带效应”(图 3c-d)。

2.3   断裂解析

顺北油田志留系整体以雁列式负花状正断层组合样式为主，呈右阶雁列式展布，走滑伸展作用强烈^[16-17]。结合志留系漏失位置及剖面、平面地震特征，志留系雁列断裂地震响应特征主要为地震反射同相轴牵引挠曲弯折或错断。顺北地区志留系断裂结合断裂生长指数分析，志留系克孜尔塔格组、依木干他乌组断裂生长指数大于1(图 4)，表明雁列断裂单元内厚度较断裂带外厚度大，判断志留系断裂期次主要为加里东晚期到海西早期。其中断层生长指数具体为同一地层断裂上盘地层厚度除以断裂下盘相同地层厚度，断层生长指数大于1，代表地层沉积加厚。同时结合顺北地区下古生界中、下奥陶统目的层鹰山组、蓬莱坝组及寒武系发育刚性灰岩、白云岩地层及上奥陶统桑塔木组—志留系发育塑性碎屑岩地层，刚性—塑形叠合地层结构在走滑应力作用下控制形成了顺北走滑断裂带“下伏陡直走滑段与上覆雁列正断层”分层变形^[18-21](图 5)。其中下伏陡直走滑段断裂倾角为80°~90°，在部分区域寒武系膏盐岩段发育逆冲滑脱断裂，上覆雁列正断层倾角为50°~55°。

图  4  塔中顺北地区志留系各组地层变化

a.不同层断裂相干检测属性及剖面位置(a1.克孜尔塔格组顶；a2.柯坪塔格组顶)；b.断裂解释及地层变化规律；c.志留系不同地层生长指数。

Figure  4.  Stratigraphic variation of Silurian strata in Shunbei area, central Tarim Basin

下载: 全尺寸图片 幻灯片

图  5  塔中顺北地区走滑断裂带特征

剖面位置图见图 1的a—a’。

Figure  5.  Characteristics of strike-slip fault zone in Shunbei area, central Tarim Basin

下载: 全尺寸图片 幻灯片

3.   应力扰动特征分析

3.1   现今最大主应力方向确定

塔里木盆地顺北地区不同地层应力场分布特征具有微弱差异。利用实验测试和测井资料计算可以确定研究区目的层的岩石力学参数和关键井的地应力大小。水平主应力差异、裂缝开启程度、断裂强度等影响横波波速各向异性。水平最大主应力方向或断裂走向与快横波方位角对应，偶极阵列声波测井数据分离后可得快、慢横波速度及方位。在顺北地区志留系砂泥岩地层中，横波传播速度沿最大水平主应力方向快，沿最小水平主应力方向速度慢。基于顺北地区SHZ1井等已钻井志留系偶极阵列声波测井数据波速各向异性方向分析，确定研究区区域现今地应力场的方向为北东方向54°左右(图 6)。

图  6  塔中顺北地区SHZ1井基于偶极阵列声波测井数据志留系波速各向异性方向分析

Figure  6.  Anisotropic directions of Silurian velocity based on dipole array acoustic data from well SHZ1, Shunbei area, central Tarim Basin

下载: 全尺寸图片 幻灯片

3.2   应力扰动特征分析

3.2.1   地质模型建立

走滑断裂位置不同会产生不同性质的应力场环境，平面上不同构造单元力学性质有所差异。左阶左行式和右阶右行式在连接带表现出伸展应力场的环境。而右阶左行式和左阶右行式在连接带表现出聚敛应力场的环境。在三维地震资料解释基础上，利用实验测试和测井资料计算确定研究区目的层的岩石力学参数和关键井的地应力大小，结合井壁坍塌、椭圆井眼、诱导缝和GPS数据等被用来确定地应力的方向，将研究区中解释好的断裂带进行数字化，以顺北地区顺北4号断裂带为例，将数字化的断层结合现今地貌及断裂体系图建立地质模型(图 7)。一般断裂带较正常沉积的地层强度有弱化。根据顺北地区测井资料计算的岩石力学参数结果，参考邻井SHB5-1X井的地应力梯度：上覆岩层压力梯度为0.025 6 MPa/m，最大水平地应力梯度为0.023 4 MPa/m，最小水平地应力梯度为0.018 7 MPa/m。SHB51X井的地应力梯度：上覆岩层压力梯度为0.024 6 MPa/m，最大水平地应力梯度为0.022 5 MPa/m，最小水平地应力梯度为0.018 8 MPa/m。地应力方向参考塔里木盆地现今应力场方向，模拟时采用北东向40°，不同区块的应力场模拟分别赋以不同的力学参数。

图  7  塔中顺北4号断裂带志留系三维地质模型

Figure  7.  Three-dimensional geological model of Silurian in Shunbei No.4 fault belt, central Tarim Basin

下载: 全尺寸图片 幻灯片

3.2.2   应力扰动特征

对塔里木盆地顺北4号断裂带开展了迭代边界有限元法数值模拟，过程为基于4号断裂带构造断裂精细解析，利用三维地震数据解释的地层起伏、构造高低，建立地质模型。进一步利用已钻井资料、区域应力特征，结合力学模型，求解运算得到现今三轴主应力分布规律。相应位移、应变、应力等是根据地质体离散后利用离散单元的弹性模量、泊松比等岩石力学参数计算所得。根据不同岩石的应变特性，选用不同的网格单元类型进行地质模型的划分以得到有限元模型，地质模型网格化后，进行网格检查，其中纵横比检查、平衡度偏差检查、最大扭角检查等均在可接受范围，网格质量良好。在建立地质模型及确定各项岩石力学参数的基础上，即可按照现今应力场的方向对模型进行加载求解(本次研究边界力只包括水平构造力)。根据顺北地区现今应力状态，最终确定在模型北东—南西及北西—南东边界施加挤压应力。应力场模拟结果揭示了研究区现今地应力的平面分布规律，包括最大(S₁)、最小(S₃)和中间主应力(S₂)，并将应力场模拟结果与实际测井资料解释结果进行对比分析(图 8)。3个主应力在断裂端部集中，受断裂产状和断距的影响，平面上应力扰动带宽度不同，同一断裂的上下盘应力扰动范围呈现出非对称分布。应力扰动的不同进一步增加了井筒周围应力分布(井周应力场)的复杂性，可能会导致断裂两盘的漏失情况存在差异^[22-26]。进一步验证了地质模型和模拟结果的合理性，指导了该地区油气勘探开发。

图  8  塔中顺北地区SHB4号走滑断裂带志留系迭代边界元法应力场模拟结果

Figure  8.  Stress field simulation of Silurian in SHB4 strike-slip fault zone in Shunbei area of central Tarim Basin using iterative boundary element method

下载: 全尺寸图片 幻灯片

3.3   裂缝发育系数与漏失关系

裂缝发育强度影响泥浆漏失程度。岩石内部受应力作用发生张性、剪切破裂。利用Griffith、Mohr-Coulomb相应破裂准则分别预测张性破裂、剪切破裂，综合分析顺北地区志留系地层裂缝的平面分布情况，顺北地区4号断裂带志留系的裂缝发育系数(I_F)分布在0.61~2.69。地层I_F为1时，为岩石产生破裂临界值。根据顺北地区的I_F值分布范围，I_F值越大，岩石越容易产生破裂。针对顺北地区4号断裂带上的已钻井，统计裂缝发育系数，分析与漏失量的关系，总体表现出裂缝发育程度越高，漏失量逐渐增大。而SHB43井漏失量较大，裂缝发育系数较小，但应力差小，为46 MPa，较其他应力差小2~12 MPa，揭示该井漏失主要受控于地应力。

4.   井位优选

顺北油田雁列断裂空间差异性强，断裂相邻、相向而倾两条断裂所围限区域为雁列单元，雁列单元边界断裂强度大，其他伴生断裂断裂强度小^[27-30]；志留系雁列断裂丝带效应明显，断裂强度大，控制范围内裂缝较发育，断裂旁强能量异常区域预测裂缝较发育；结合空间立体雕刻，避开断裂解析及应力扰动分析成果中认为志留系断裂强度大、倾角大的雁列单元边界断裂、断裂应力集中部位及强能量异常断裂、裂缝发育区，结合目的层靶点空间特征，优选设计井轨迹。顺北41X井优选井口后，井深结构从前期六开优化为四开，志留系未发生严重漏失(志留系仅漏失14 m³)，节约钻井周期63 d(提速28%)。

5.   结论

(1) 顺北地区志留系克孜尔塔格组、依木干他乌组断裂生长指数大于1，揭示断裂单元内厚度较断裂带外厚度大，判定志留系断裂期次主要为加里东晚期到海西早期活动。结合剖面断裂样式及平面分布特征，断裂展布规律主要为雁列式负花状正断层组合排列样式，呈右阶雁列式展布。走滑伸展作用强烈，延伸长度较小，部分成对存在，地层破碎严重。

(2) 顺北地区志留系主应力方向为北东向54°。迭代边界元数值模拟明确研究区现今地应力的平面分布规律，包括最大、最小和中间主应力，3个主应力在断裂端部集中，受断裂产状和断距的影响，平面上应力扰动带宽度不同，同一断裂的上下盘应力扰动范围呈现出非对称分布。

(3) 避开断裂解析及应力扰动分析成果中认为志留系断裂强度大、倾角大的雁列单元边界断裂、断裂应力集中部位及强能量异常断裂、裂缝发育区，结合目的层靶点空间特征，优选设计井轨迹。顺北41斜井优选井口后，井深结构从前期六开优化为四开，志留系未发生严重漏失(志留系仅漏失14 m³)，节约钻井周期63 d(提速28%)。