随着国内新发现油气日趋以低品位复杂油气藏为主，勘探开发难度日益增加，高效开发这类复杂油气藏急需新理念与新工艺。随着地质工程一体化创新发展，以北美Bakken、Eagle Ford为代表的致密油藏，借助地质工程一体化思路，采用长水平段和多级压裂等技术实现了原来认为没有效益的低品位资源规模化效益开发，为国内复杂难动用油气藏开发提供了新的途径^[1]。1994年Alderman^[2]率先在油气领域引入“一体化管理体系”的概念，随后在实践过程中“勘探开发一体化”不断完善，形成了地质、工程、经济、管理一体化模式，成为新油气资源开发的主要方法。国外斯伦贝谢、哈利伯顿等服务公司先后提出一体化项目管理（PM或IPM）的概念并实践，2002年哈利伯顿公司本着综合一体化服务高效共赢的理念重新定义了IPM作业服务模式，并在国外多个油气田勘探开发中进行了实践。

中国在20世纪90年代引入“一体化”理念，2015年吴奇等^[3]提出地质—工程一体化理念和“钻井品质”概念并在中国南方海相页岩气进行了实践，实现了中国南方海相页岩气的经济有效开发；2015年陈颖杰等^[4]利用开发地质、精细建模、井眼轨迹设计等一体化理念进行了页岩气地质工程一体化导向钻井探索。2016年5月31日，由中国石油学会石油工程专业委员会等单位组织召开的“第一届油气地质工程一体化论坛”，极大地推动了国内地质工程一体化的深入研究。2017年，胡文瑞院士^[1]在对国内外开发对比的基础上提出地质工程一体化是实现复杂油气藏效益勘探开发的必由之路。章敬等^[5]基于多学科协同，采用地质与工程技术相结合，管理与工程相结合，探索出工程技术服务一体化总包模式开发新疆玛湖致密油。王昕等^[6]提出以规模发现与效益动用为主线，聚焦从成藏研究到压裂参数优化的全过程地质工程一体化技术路线。鲜成钢等^[7]利用岩心、地震等数据，基于单井约束、迭代逼近的方法建立了三维三应力模型。常少英等^[8]将工程人员“嵌入”地质研究团队，实现了地质、钻井、油藏等学科协同工作，探索了剩余油挖潜研究中的地质工程一体化路线，实现了油田的高产稳产。2018年，张以明等^[9]根据超高温超深碳酸盐岩储层特点，建立了储层评价、改造工艺优选与设计等地质工程一体化体积改造技术。2019年，张合文等^[10]根据低渗透、薄层边际油藏难动用的问题，采用地质与钻完井、储层压裂改造和采油工程等的协同，实现了Zanazour油田水平井分段压裂有效开发。2020年，杨智等^[11]以致密油“甜点”识别为基础，根据储层精细刻画，分析了国内几种典型致密油油藏地质工程一体化的进展。冯张斌等^[12]在储层开发地质精细描述的基础上，进行了鄂尔多斯盆地延长组7段致密油地质工程一体化建模和油藏数值模拟方面的研究，重点探索了水平井井眼轨迹设计与钻井优化方面的工作。2021年，蒋廷学等^[13]按照地质工程一体化理念，采用压前“甜点”预测储层评价、压裂大数据协同优化、渗吸实验评价和压后评估与管理等技术建立了体积压裂地质工程一体化技术。

目前，国内地质工程一体化在多学科交叉融合方面进行了大量探索，为国内多个复杂油气藏经济有效开发提供了坚实的基础，但是在储层精细描述—地应力场三维模型—井网优化—开发效果评价等多学科协同提高超低渗油藏开发方面仍有提升的空间，储层非均质性、地应力场、井网和改造技术一体化的研究在鄂尔多斯盆地页岩油藏的适用性仍需要进一步探索。

姬塬油田位于鄂尔多斯盆地中西部，跨天环坳陷和伊陕斜坡两大构造单元，主要位于定边、吴起两县（图 1）。姬塬油田构造平缓，为西倾单斜，坡降为6m/1000m左右，局部发育利于油气聚集的微构造。鄂尔多斯盆地延长组自下而上分为长10至长1共10个油层组^[14-16]，主要生产层位是长6段及长7段。姬塬地区长9段和长7段的黑色页岩及暗色泥岩是主要烃源岩。长6段厚度在120m左右，自下而上可进一步划分为长6₃、长6₂和长6₁共3个亚段。长7段厚度在110m左右，也已实现了规模建产，自下而上可进一步划分为长7₃、长7₂和长7₁共3个亚段。

图 1

图 1 鄂尔多斯盆地姬塬油田位置图 Fig. 1 Structural location map of Jiyuan Oilfield in Ordos Basin

姬塬地区延长组发育一套典型的河流—湖泊沉积体系，姬塬油田长7段属三角洲前缘和半深湖—深湖沉积。长7段沉积期，鄂尔多斯湖盆鼎盛发育，半深湖—深湖范围为延长组沉积期最大^[17]。长7段单砂体厚度分布以2~4m为主，个别砂体厚度可达20m以上，平均砂体厚度不超过10m，以零星分布为主；砂地比为27.88%~34.86%，其中长7₂亚段最高，长7₁亚段次之，长7₃亚段砂体最不发育^[18]。储层岩性以细粒长石砂岩、岩屑长石砂岩为主，颗粒分选差，填隙物以绿泥石、高岭石和铁方解石为主，总含量平均为12.9%~14.7%^[19]。长7段孔隙类型以粒间孔和溶蚀孔为主，喉道类型以弯片状和片状为主，发育微米—纳米级孔喉，半径在0.1~1.0μm之间的亚微米级喉道占比超64%，中值半径为0.16μm。研究区孔隙度为5.00%~16.95%，平均为8.86%，低孔占比约为69%；渗透率为0.01~1.65mD，平均为0.2mD，根据行业标准《油气储层评价方法》（SY/T 6285—2011），储层物性以低孔—特低孔、超低渗为主，孔隙度与渗透率整体呈正线性相关^[20]。

岩心观察时发现一定数量裂缝，并且以高角度裂缝为主，裂缝倾角介于30°~90°，以60°~80°为主，裂缝走向介于25°~180°，以75°~80°和85°~90°两个范围分布最多，裂缝走向玫瑰花图表明裂缝方向主要为北东东67.5°、北东东80°和北东东88°三个方向^[21]。前期地质认识表明，长7段油藏为页岩油藏，单个油藏规模较小，没有统一油水界面，油水分布关系复杂；研究区长7₃亚段、长6₁亚段均有一定数量油藏发现。试油数据显示长6段油气更为富集，长6₁亚段平均含油饱和度最高，长6段各亚段孔隙度、渗透率、有效厚度和含油饱和度平面分布趋势基本相同，表明研究区为主要受岩性控制的油藏；油藏平面呈条带状分布，纵向呈透镜体状分布，局部发育一定数量边水或零星水体^[16]，少部分井为油水同产井或水井，研究区水型以CaCl₂型地层水为主，少量分布Na₂SO₄型和MgCl₂型地层水，平均矿化度为61g/L。安83井区长7段油藏自2008年5月投产以来，已经成为姬塬油田长7段油藏试验区块，为长7段油藏开发部署提供了良好的指导，本文以安83井区长7段油藏为例进行研究。

安83井区长7段储层处于辫状河三角洲前缘亚相，储层厚度变化大，发育水下分流河道、分流间湾等沉积微相，水下分流河道砂体厚度为5~25m，单个水下分流河道宽度为12~23m。姬塬油田安83井区长7段砂体平面延伸距离长，局部超16km，纵向多层叠置厚度可达近30m。

受沉积环境影响，长7段砂体普遍存在复杂的纵向、横向叠置关系，平面上砂体展布方向、高渗条带展布方向与最大水平主应力方向不一致，长7段砂体主要沿北东40°展布；渗透率等值线整体沿北东40°展布，但高渗条带以沿北东5°~北东10°展布为主；最大水平主应力以北东70°展布为主，局部存在最大水平主应力方向变化。用统一井排距和井排方向布井造成研究区储量控制程度低，适应性不足。目前姬塬油田采用了矩形、菱形、正方形和水平井网等不同井网形式、6种不同井排距，以菱形反九点井网为主要开发井网，井排距约为480m×200m，井排方向以北东60°~80°为主，初期井网按照统一井排距在砂体范围内均匀分布，后期部分井区逐步采用了水平井五点法井网，井距为300~600m，水平段长度为500~1000m。由于姬塬油田长7段砂体纵向多期叠置，横向存在多砂体叠置展布，造成现有开发井网控制程度不高。通过对安83井区部分井组计算表明，其水驱控制程度平均为58.85%，按照水驱储量控制程度评价标准可判断水驱控制程度差。注水开发区域能量补充不足，造成其产量递减率偏高，第一年递减率超过20%，水平井开发区年递减率甚至超过45%。

姬塬油田长7段油藏普遍存在纵向多期砂体叠置，各油层组砂体在油田不同区域呈现不同特征，油田西部区块砂体规模小，隔夹层发育，优质储层横向连通性较好，而中等储层和差储层横向连通性较差；油田东部区块砂体横向规模大且连片性好，纵向叠置厚度大；油田中部区块砂体纵向以薄砂体为主，具有隔夹层发育、横向连通性差、零散分布的特点。受沉积微相和水动力等影响，长7段各亚段水下分流河道等沉积微相以北西向带状分布为主，局部呈北东向带状分布，通常存在多条水下分流河道并在局部交叉汇集，孔隙度、渗透率平面上也以北西向带状分布为主，局部存在一定弯折或其他走向，沿砂体展布方向物性要好于垂直砂体展布方向。姬塬油田长7段油藏类型为岩性油藏，没有统一油水界面，但是在油藏中部区块砂体低部位、砂体边部及其他区块部分砂体局部低部位零散分布一定数量天然水体。由于研究区整体物性较差，必须经过压裂才能工业化开发，压裂过程中易与天然水体沟通，造成部分井投产后即产水，控水困难。此外，井网采用统一井排距和井排方向在砂体范围内均匀布置，井网设计时没有建立三维地应力场模型，根据现有动态和压裂监测表明部分井区井排方向与裂缝展布方向不一致，致使部分井组注水开发后过早见水。同时，沿最大水平主应力北东70°设置的井排方向与砂体展布方向不一致，没有建立起井网、储层改造压裂、砂体展布与渗透率展布之间的协同优化，导致部分注水井区出现地层压力保持水平低，油井暴性水淹，控水困难(图 2)。

图 2

图 2 安231-38井生产曲线图 Fig. 2 Production curve of Well An 231-38

地质工程一体化是对原有工作流程和结构的变革。实施地质工程一体化首先需要大油藏的概念，一切以油藏为中心，通过地质、油藏、工艺、经济等多学科技术协同，实现油藏最佳规模效益开发。根据地质工程一体化目标的需要，建议采用扁平化地质工程一体化组织结构（图 3）和多学科协同地质工程一体化油藏工程研究流程（图 4）。

图 3

图 3 地质工程一体化组织结构图 Fig. 3 Organization chart with geology and engineering integration

图 4

图 4 地质工程一体化油藏工程研究流程图 Fig. 4 Flow chart of oil reservoir engineering study with geology and engineering integration

地质工程一体化实践中采用扁平化管理，围绕油藏经济高效开发这一中心任务，坚持多学科协同解决问题，建立闭环地质工程与油藏开发一体化系统，通过对油藏进行持续跟踪与更新，持续优化井网与开发技术对策，始终以最优决策运用到油藏开发中。项目实施过程中，虽然不同阶段各学科工作量和工作重点有所不同，但均需要以决策目标为中心，整个流程中，围绕油藏模型，将地质研究作为基础环节，通过地震数据、测井数据、岩心分析等，地球物理工程师和开发地质工程师能够建立精细构造模型。同时，开发地质工程师和油藏工程师等协同，利用钻井数据、三轴应力实验结果、测井计算力学参数等建立单井岩石力学参数和地应力模型。此外，根据生产动态数据、动态监测等，油藏工程师、采油工程师等对储层进行认识。通过以上基础模型，在综合地质研究基础上建立研究区三维属性模型、三维力学模型及三维地应力场模型，基于此开展井网优化、井轨迹优化部署和其他工艺设计，为后续随钻跟踪提供支撑。基于三维属性模型，根据随钻监测与跟踪数据，动态调整预测砂体变化，指导随钻调整，预测钻井风险，切实提高钻井砂体钻遇率，降低钻井风险。结合新钻井随钻数据、录井和测井数据进行储层评价，对三维属性模型进行实时调整，根据调整后的三维属性模型优化完井和压裂方案。进而根据压裂监测、压后效果分析调整完善三维属性模型、三维力学模型和三维地应力场模型，根据生产动态在油藏工程分析的基础上，利用油藏数值模拟进行历史拟合，持续优化更新，对储层进行再认识。针对开发存在的问题，根据油藏工程分析提出增产增注措施，利用油藏数值模拟等方法进行优化，再次更新储层认识，形成持续更新、以油藏高效开发为核心的地质工程一体化油藏工程研究体系。

以三维属性模型为目标，通过地球物理、开发地质、油藏等学科协同，开展油藏精细描述工作。根据地震断层解释结果进行断层建模。首先，在对岩心观察及地震反演属性等分析的基础上，以标志层识别为主线，借助现代沉积学理论，充分利用地震层位解释、测井标志层识别和生产动态裂缝识别等，进行研究区地层对比划分，建立研究区各储层的构造地质模型。姬塬油田延长组通常自下而上分布有K0—K9共10个标志层，其中位于长7₃亚段底部湖相黑色页岩（张家滩页岩）底部的凝灰岩为K1标志层，电性具有高声波时差、高电阻率、高自然伽马和低密度的“三高一低”特征，长6段和长7段交界的K2标志层则体现出高自然伽马、高声波时差、高自然电位、高电阻率的特征（图 5）。其次，以细分亚段为单元，研究分层构造特征。参考区域沉积背景和早期沉积相研究资料，结合岩心观察，运用沉积学原理，通过油藏储层的地质沉积特征在各种测井信息上的具体分析，明确研究区沉积特征，建立起研究区单井测井相模式，绘制沉积微相剖面图和平面分布图。根据中国石油长庆油田公司勘探开发研究院研究表明，姬塬油田延长组物源方向主要为北东向和北西向，油田西部区块距离北西向物源近，而油田东部区块则主要受北东向物源控制。由于砂体受沉积相的控制，在沉积相研究的基础上可进一步落实油藏砂体分布规律，绘制砂体展布图。然后，在一次测井解释结果基础上，结合室内岩心分析结果，对实际开发油层性质校正并进行二次测井解释，编绘油藏剖面图。随后，开展储层物性分布特征及其非均质性、隔层的分布和油层连通状况、注采关系等研究，分类评价储层物性分布与油井产能分布特征。最后，应用油藏地质建模软件并采用多参数协同砂体约束建模技术，以尽可能进行精细的地质分层，建立研究区描述构造、储层物性、流体空间分布的静态三维属性模型，三维属性模型主要步骤建立流程见图 6。

图 5

图 5 研究区典型井地层划分测井响应图 Fig. 5 Stratigraphic division and logging response in typical wells in the study area

图 6

图 6 研究区三维属性模型建立流程图 Fig. 6 Flow chart of 3D attribute modeling in the study area

由于鄂尔多斯盆地储层相对稳定，地应力方向相对单一，长期开发中形成了以北东向为最大水平主应力方向的基本认识，各区块布井过程中通常会根据部分井测试结果确定的某一最大水平主应力方向进行均匀井网设计，而很少系统建立三维力学模型和三维地应力场模型。同时，姬塬油田开展的岩石力学参数室内实验和矿场测试较少，造成直接基于实验结果建立三维力学模型的复杂性和难度较高。为了建立可靠合理的三维力学模型和三维地应力场模型以指导油田开发，本文综合利用室内实验、测井数据、岩石力学理论及压裂资料等，采用随机建模方法建立了一套三维力学模型和三维地应力场模型，同时利用钻井、压裂等资料对模型不断进行更新和完善，持续提高模型的精度和可靠性。

首先，开展室内岩心岩石力学实验，确定区块岩石力学参数与地应力参数。其次，将钻井数据、固井数据和压裂数据等相结合，对已有的岩石力学参数和地应力参数计算模型进行修正，随后利用修正后的计算模型和各井测井资料计算得到单井岩石力学参数和单井地应力参数，建立单井剖面。接着，基于单井岩石力学参数和单井地应力参数，利用建模软件，采用随机建模的方法，建立三维力学模型和三维地应力场模型，并根据已有监测井对多个实现进行优选。最后，及时将新增数据和监测资料导入，持续更新、修正并检验模型。多方法融合的三维力学模型和三维地应力场模型建立流程见图 7。

图 7

图 7 多方法融合三维力学模型和三维地应力场模型建立流程图 Fig. 7 Flow chart of 3D mechanic modeling and 3D in-situ stress field modeling by fusing multi method

常规测井中声波时差测井常为纵波，而进行岩石力学参数计算时需要用到纵波和横波两种，因此首先通过实验室测定研究区纵波与横波关系。实验时用波长为毫米级的超声波来反演推断地震波在地质体中的传播特征，根据岩样声学特征测试纵波和横波实验速度数据，可以建立研究区纵波与横波转化关系式。根据实验结果可以得到其纵波与横波关系如图 8所示，该区纵波与横波转换方程为

图 8

图 8 鄂尔多斯盆地安83井区纵波与横波转换关系图 Fig. 8 P-S wave conversion diagram in An 83 well block in Ordos Basin

式中  ν_s——横波速度，m/s；

ν_p——纵波速度，m/s。

在此基础上可以利用测井资料计算岩石泊松比、杨氏模量、剪切模量、体积模量等弹性力学参数，相关模型可从相关资料获得。由于弹性力学参数模型计算时需要用到静态弹性模量等静态参数，而测井资料直接求取结果为动态弹性模量等动态参数，因此需要根据相关室内实验建立动态与静态弹性模量等岩石力学参数转换关系。

利用“TAW-1000深水孔隙压力伺服实验系统”对区块岩样进行岩石力学参数测定工作，实验内容包括密度、弹性模量、泊松比、单轴抗压强度、单轴抗拉强度、三轴抗压强度、内聚力和内摩擦角等。可以根据实验结果和测井计算结果回归得到修正后的各岩石力学动态、静态弹性模量和泊松比之间的转换方程。

如修正后动态、静态弹性模量的转换公式为

式中  E_s——岩石的静态弹性模量，MPa；

E_d——岩石的动态弹性模量，MPa。

同时也可以根据内聚力和内摩擦角等实验结果，得到修正后岩石的抗压强度计算公式、抗压强度与抗拉强度的关系式、内聚力计算模型和内摩擦角模型。进而编写软件，利用测井曲线计算得到研究区各井强度参数和弹性参数等岩石力学参数。在黄荣樽模型的基础上，采用改进后弹簧组合模型来计算地应力^[22]。基于以上模型，利用测井、压裂和室内实验结果等资料可以得到各井的单井岩石力学参数和地应力参数。

考虑到成熟的三维地质建模软件建立的模型不仅可以忠实于控制点的实测数据，而且还可以对井间数据实现可靠精度下的内插和外推，因此本文基于单井计算得到的地应力参数，尝试利用三维地质建模软件建立区域三维地应力场模型。目前三维建模方法可分为确定性建模方法和随机建模方法两大类^[23-24]。考虑到实测中发现不同区域地应力变化的复杂性，因此采用基于象元的随机建模方法建立三维地应力场模型。

在建立储层模型的时候，要充分利用前期地质研究的成果，如考虑主应力方向，以建立正确的变差函数模型，并选择合适的模拟方法。将计算得到的单井地应力（三应力）作为新的测井数据，进而进行随机模拟，选用高斯模型和序贯高斯模拟方法建立研究区三维地应力场模型。图 9、图 10为基于此方法建立的鄂尔多斯盆地姬塬油田安83井区最大水平主应力和最小水平主应力的三维模型。该区最大水平主应力方向主体为北东60°~北东80°，局部如安233-36、安231-38等井最大水平主应力方向为北东20°~北东50°。

图 9

图 9 安83井区最大水平主应力三维图 Fig. 9 3D diagram of the maximum horizontal principal stress in An 83 well block

图 10

图 10 安83井区最小水平主应力三维图 Fig. 10 3D diagram of the minimum horizontal principal stress in An 83 well block

鄂尔多斯盆地长6段、长7段油藏井网设计时，通常在砂体分布范围内以均匀分布井网为主，井排方向与最大水平主应力方向相同。以安83井区为例，该井区采用6套不同井排距开采，其中中部采用450m×140m矩形井网；北部采用350m×150m菱形井网；西部采用220m×220m正方形井网；中西部采用400m×140m矩形井网；东南部采用水平井井网；南部采用480m×200m菱形井网^[25]。井网主要井排方向为北东60°~北东80°，与最大水平主应力方向（裂缝发育方向）平行，生产动态反映出压力保持水平低、主向井易水淹、水驱控制程度低、个别侧向井水淹等问题。针对实际开发出现的问题，考虑到压裂裂缝走向与地应力的关系，以及裂缝与砂体沉积、渗透率分布的关系，提出了基于动态三维属性模型和三维地应力场模型，综合考虑裂缝和砂体展布方向的井网部署方案。

首先，根据砂体展布方向、渗透率平面分布特征和室内岩心渗透率实验结果等，确定沿砂体展布方向的渗透率与垂直砂体展布方向渗透率的极差倍数，将极差倍数作为确定井排距关系的重要依据，同时建立归一化渗透率分布场。接着，根据三维地应力场模型，刻画压裂裂缝分布规律，确定人工裂缝展布方向，重点关注与研究区北东60°~北东80°方向不一致的最大水平主应力分布区域。然后，对三维最大水平主应力进行归一化处理，建立归一化最大水平主应力分布场。随后，结合其他区块压裂裂缝与原始储层渗透率等效处理关系，将归一化后的渗透率分布场和归一化最大水平主应力分布场进行矢量加和，将加和后的矢量分布方向作为井排方向。最后，根据渗透率的极差倍数和已开发井网效果分析结果，采用菱形反九点井网进行井网部署，后期根据动态特征可将菱形反九点井网调整成矩形井网；对于滚动开发井，根据动态三维属性模型和三维地应力场模型采用相同流程确定井网。

地质工程一体化中的开发技术政策应该是整个一体化项目中的落脚点与新的起始点，围绕油藏高效开发的目标，充分借助地质工程一体化的理念，在精细描述和井网优化的基础上，先通过常规油藏工程分析，进一步明确现有开发效果及存在问题，再借助油藏数值模拟技术，实现地质与油藏工程和压裂的耦合模拟^[26]。

开发技术政策优化流程见图 11。首先，通过常规油藏工程方法适时开展动态分析；通过区块生产曲线、井组生产曲线、单井生产曲线及时了解各部分生产动态；通过水驱储量控制程度计算井网与地质因素的匹配关系；通过水驱特征曲线、含水率变化关系等明确水驱开发效果；通过递减规律分析明确区块开发状况，及时了解区块开发状态和存在问题；通过对常规油藏工程各指标计算，对于出现含水率突然快速增加等情况，及时分析重点井组、井开发状态和可能见水层位、方向等。接着，根据建立的三维属性模型和三维地应力场模型，结合地质认识和裂缝系统认识，利用油藏数值模拟方法，将包含刻画后裂缝系统的粗化三维属性模型导入数值模拟软件。根据研究区室内相渗及高压物性等实验、井史等建立油藏数值模型并根据动态认识和测试资料成果进行历史拟合，重点拟合中高含水井及含水暴性增加井，利用数值模拟方法确定地质裂缝等认识，定量确定研究区剩余油分布，明确剩余油富集规律。然后，根据修正后的油藏数值模型，对于新部署井或完善井网区域，进行井网优化部署模拟，确定合理井网，进而确定新布井开发层位、合理压裂改造参数；对于已开发井区，进行开发技术政策优化，进一步优化各井注入量、产量等开发参数；对于拟改造井或调整井，开展重复压裂、井网调整、层位调整等压裂参数优化，确定最佳改造井及改造措施。最后，将新的油藏数值模拟认识加载到三维属性模型中，持续更新模型，根据新钻井和生产动态新认识及时将新增数据和监测资料导入，持续更新、修正并检验三维属性模型和油藏数值模拟模型。

图 11

图 11 地质工程一体化开发技术政策优化流程图 Fig. 11 Flow chart of development countermeasure optimization with geology and engineering integration

以安83井区为例，区内共有21口油井、13口注水井，最早投产的井是安236-39井和安234-33井，均于2010年9月投产，截至目前，两口井平均日产油量小于1t。井区内含水率以40%~80%最多，区块整体含水率为70%左右，安234-37井等4口井含水率均超过90%，采用连通关系统计法估算水驱储量控制程度，计算得出安83井区水驱控制程度为58.85%，可判断水驱控制程度差。通过对含水率与采出程度回归可知，研究区含水率呈“S”形变化，表现为早期含水率快速增长，随着油田的开发，中期含水率缓慢上升，后期也有明显的小幅度上升，含水率与采出程度关系如下

式中  f_w——含水率；

R_K——可采储量采出程度。

根据Arps递减模型，安83井区月产量呈双曲递减，递减公式如下

式中  Q——月产量，t；

t——生产时间，月。

根据各井动态分析，特别是注水量波动和周边油井产量变化可知，安83井区在平面方向上存在方向性见水，见水方向有两个方向，主要方向是北偏东方向，小部分见水方向是北偏西方向，局部存在裂缝性水淹。

根据油藏工程动态分析认识，结合三维属性模型、三维地应力场模型、室内实验结果，建立安83井区油藏数值模型。对于安83井区长7段油藏设置北西方向35°为I轴，垂直此方向为J轴，网格步长沿I、J方向均为10m，纵向上划分为5个网格，网格数为N_I×N_J×N_K=206×161×5=165830，其中N_I、N_J、N_K分别为I、J、K轴方向的网格数。本文对安83井区长7段油藏34口油水井全部进行了拟合，采用的方式为：给定油井单井产液量、给定注水井注水量，重点拟合含水率等指标，拟合历史从2010年8月开始，模拟过程以1个月为一时间段。模拟区内34口井中31口井的拟合程度较高，拟合率达到91.2%。

安83井区长7段油藏剩余油分布主要受到渗透率非均质性、压裂裂缝与井网系统、注采不平衡、岩性的影响。模拟区目前剩余油分布见图 12。

图 12

图 12 安83井区模拟区目前剩余油分布图 Fig. 12 Distribution of remaining oil in the simulation area in An 83 well block

对于安83井区而言，目前各层动用情况不均衡，剩余油从纵向看主要集中在中下部；而从平面看，大部分井周边30~100m范围内的剩余油饱和度较低，含水率较高，但是不同井区差异较大，部分井区目前动用砂体周边基本水淹，如安233-35注水井周边剩余油饱和度较小，此外对于该区含水较高地带主要受区域内裂缝系统影响，整体沿裂缝方向含水突进比较明显。

为了进一步提高安83井区采收率，在井网优化的基础上，对该井区采取重复体积压裂增产措施。考虑到裂缝半长、主裂缝导流能力和次裂缝导流能力等多个因素均会影响压裂效果，因此采用正交实验的方式模拟多因素对压裂效果的影响，以便选择最佳的施工参数。采用L₉(3⁴)正交试验表（表 1），影响因素包括裂缝半长、缝宽、主裂缝导流能力和次裂缝导流能力4个，每种因素考虑3种水平。安83井区相关模拟结果见图 13。

表 1

表 1 安83井区注水开发下体积压裂优化正交方案设计参数表（L9(34)） Table 1 Parameters of volume fracturing optimization orthogonal scheme by using water injection development in An 83 well block(L9(34)) 因素 裂缝半长/m 缝宽/m 主裂缝导流能力/(D·cm) 次裂缝导流能力/(D·cm) 实验1 100 30 15 0.4 实验2 100 50 20 0.7 实验3 100 70 25 1 实验4 120 30 20 1 实验5 120 50 25 0.4 实验6 120 70 15 0.7 实验7 140 30 25 0.7 实验8 140 50 15 1 实验9 140 70 20 0.4

表 1 安83井区注水开发下体积压裂优化正交方案设计参数表（L9(34)） Table 1 Parameters of volume fracturing optimization orthogonal scheme by using water injection development in An 83 well block(L9(34))

图 13

图 13 注水开发时安83井区体积压裂各方案2015年末累计产油效应分析曲线 Fig. 13 Cumulative oil production analysis curve at the end of 2015 of each volume fracturing scheme during water injection development in An 83 well block

通过图 13看，采用目前井网注水开发时，对于体积压裂各方案模拟结果表明，如果仅仅考虑模拟期末累计产油量而不考虑其他因素，则最佳的组合应该是裂缝半长为120m、缝宽为30m、主裂缝导流能力为15D·cm、次裂缝导流能力为1D·cm。而从影响因素方面分析，主裂缝导流能力影响最大，次裂缝导流能力影响最小。可以看出，存在注水井时过长的裂缝可能会造成采油井过早见水，因此并不是主裂缝越长越好，适当的主裂缝长度尤为重要。此外由于裂缝并非直接与采油井处于同一连通方向，因此缝宽和主裂缝导流能力也不易过宽或过大，否则会增加油井过早见水的可能；而增加次裂缝导流能力则可扩大泄油面积，有利于提高采出程度。

根据优化流程，选择了安231-44等5口井进行重复压裂改造，改造前后效果见表 2。本次设计5口井进行施工，从施工前后对比看各井均取得了不错的效果，而与其他井相比也有更好的效果，目前本次优化设计井日产量较邻井增产30.46%。

表 2

表 2 安83井区增产措施施工前后效果对比表 Table 2 Comparison of production results before and after well stimulation in An 83 well block 井号 措施前 措施后 日产液/m3 日产油/t 含水率/% 日产液/m3 日产油/t 含水率/% 安231-44 0.68 0.47 30.93 5.31 1.28 75.1 安233-38 1.35 0.57 57.66 7.26 1.76 76.6 安234-31 — — — 5.05 1.93 60.9 安241-24 — — — 5.01 2.43 50.2 安277-28 1.54 0.02 96.01 4.6 1.71 61.1

表 2 安83井区增产措施施工前后效果对比表 Table 2 Comparison of production results before and after well stimulation in An 83 well block

（1）基于岩心、测井资料和压裂施工监测相结合的单井岩石力学参数及地应力计算方法，利用随机模拟方法建立储层三维地应力场，确定了安83井区最大水平主应力方向主要为北东60°~80°左右，而在安231-38等井附近最大水平主应力为北西30°左右，基于地应力分析设置人工裂缝后数值模拟结果与目前动态分析结果比较吻合。

（2）形成了精细地质描述—地应力分析—地质建模—油藏工程设计—数值模拟一体化技术方法，优化出安83井区注水开发时最佳的压裂工艺参数组合应该是裂缝半长为120m，缝宽为30m，主裂缝导流能力15D·cm，次裂缝导流能力为1D·cm，主裂缝导流能力影响最大，而次裂缝导流能力影响最小，优化后压裂井平均日产量较邻井增产30.46%。

（3）在动态地质模型的基础上，利用油藏数值模拟方法持续优化井网和开发技术政策，有效指导了安83井区等姬塬油田长7段高效开发。