石油工程计算书

石油工程计算书1.工程概况与基础地质参数计算本部分旨在明确石油工程计算的基础输入数据，涵盖地质力学参数、流体物性及环境参数，为后续的井身结构设计、水力学计算及强度校核提供准确的数据支撑。所有计算均基于目标区块的测井数据、地层测试数据及邻井类比分析结果。1.1地层压力剖面计算地层压力是钻井液密度设计的核心依据。本工程采用Eaton法计算地层孔隙压力，利用修正的Eaton法或Matthews-Kelly法计算破裂压力。根据目标区块的地震层速度数据及测井资料，选取关键深度点进行计算。对于上覆岩层压力梯度（Sv），采用积分法求得：=其中，为地层体积密度，随深度变化。针对某关键层位（深度3200米），计算其孔隙压力梯度（）：设定上覆岩层梯度=2.35MP实测声波时差Δt=320Eaton指数n=计算得出孔隙压力梯度：=代入数值计算后，该层位孔隙压力系数为1.18，即当量钻井液密度需至少为1.18g/1.2地层三轴应力强度分析为了确保井壁稳定性，需计算地层的三个主应力：垂直主应力（）、最大水平主应力（）及最小水平主应力（）。利用多臂井径测井及成像测井数据，结合声波各向异性分析，确定构造应力系数。计算模型如下：==其中，ν为泊松比，取值0.25；α为Biot弹性系数，取值0.85。经计算，在3200米处，最小水平主应力梯度约为1.85g/c，最大水平主应力梯度约为2.101.3原油及天然气物性参数计算基于高压物性（PVT）取样分析报告，对流体参数进行拟合计算。原油体积系数（）采用Standing相关式进行计算：=其中，为溶解气油比（/），为天然气相对密度，为原油相对密度，T为地层温度（F）。在饱和压力条件下，计算得=1.25原油粘度利用Beggs-Robinson相关式进行死油脱气粘度及饱和油粘度的修正，计算结果为地层条件下原油粘度4.5m天然气偏差系数（Z因子）利用Dranchuk-Abu-Kassem状态方程进行迭代求解，输入对比压力和对比温度，计算得地层条件下Z=基础地质参数汇总表如下：参数名称符号单位数值(3200m处)数据来源垂深TVDm3200.00测井解释孔隙压力梯度$P_p$g/cm³1.18Eaton法计算破裂压力梯度$F_p$g/cm³1.85Matthews-Kelly法上覆岩层压力梯度$S_v$g/cm³2.35密度积分地层温度$T_{res}$$^\circ$C112.0地温梯度推算原油粘度$\mu_o$mPa·s4.5PVT拟合天然气偏差系数Z无因次0.89DAK方程2.井身结构设计与套管强度校核井身结构设计是保障钻井安全、实现油气层有效隔离的关键环节。本节依据地层孔隙压力和破裂压力剖面，采用自下而上的方法进行必封点确定，并对套管柱进行三轴应力强度校核。2.1必封点确定与套管下入深度计算依据“抽吸压力允许值”和“激动压力允许值”原则，确定各开次套管下入深度。设定抽吸压力当量密度=0.03g/c，激动压力当量密度对于表层套管，需封隔上部疏松地层及淡水层。计算深度设定为500米，该处破裂压力梯度约为1.60g/c。设计钻井液密度为1.10+1.10+对于技术套管，需封隔高压层。在3200米处，孔隙压力系数为1.18。考虑到上部低压层（如1500米处，破裂压力1.65），需进行精确核算。最大钻井液密度=1.18在1500米处，能否承受1.23的钻井液密度？井底最大循环压力ECD=1.23+由于1.27<2.2套管柱强度校核计算采用APIRP5C3标准及ISO13679标准进行强度校核。考虑内压、外挤及轴向载荷的复合作用。2.2.1生产套管（ϕ139.7mm）抗内压校核选用钢级P110，壁厚9.17mm。API抗内压强度计算公式：=其中，为材料最小屈服强度（110ksi），t为壁厚，D为外径，K为修正系数（取0.875考虑壁厚公差）。计算得≈73.5工况分析：关井最大井口压力预期为35MPa，井底天然气柱压力约10MPa。总内压载荷=35安全系数SF2.2.2生产套管抗外挤校核API抗外挤强度需考虑屈服挤毁与塑性挤毁。对于D/计算得≈55工况分析：考虑后期注水或地层压力亏空导致套管外液柱压力高于管内。设管内全掏空，管外为完井液密度1.35g/cm³。最大外挤载荷=1.35安全系数SF2.2.3套管抗拉及三轴应力校核计算套管在空气中的重量及泥浆中的浮重。分段计算轴向力：井口处轴向力最大。套管总长3200米，每米重量约34.2kg/m。空气中总重=3200考虑浮力系数=1最大轴向载荷=1074API抗拉强度=2000安全系数SF进一步应用VonMises三轴应力准则进行复合应力校核，由于内压与外挤存在抵消效应，且轴向拉力适中，三轴应力校核通过。井身结构设计参数表：套管层次尺寸下入深度钢级壁厚抗内压抗外挤抗拉强度备注导管508mm30J5512.7---封隔地表表层套管339.7mm500J559.6518.512.84500封隔淡水层技术套管244.5mm2800N8011.0535.228.56800封隔高压层生产套管139.7mm3180P1109.1773.555.02000油层完井3.井眼轨迹设计与测量计算针对目标区块的地质特征及地面条件，本设计采用“S”型五段制剖面进行定向井设计，以实现精准中靶及满足防碰扫描要求。3.1轨迹设计参数输入设计依据：靶点垂深3200米，水平位移600米，井口坐标与靶点相对坐标ΔN设定造斜点（KOP）为1500米。第一增斜段造斜率=/稳斜段井斜角α=降斜段造斜率=/入靶井斜角=（以便于后期作业）。3.2关键节点计算利用圆柱螺线模型或最小曲率法进行轨迹计算。增斜段计算：由造斜点1500米增斜至最大井斜角45度。所需狗腿度长度Δ=增斜段末垂深增量ΔT增斜段末水平位移增量ΔD稳斜段计算：设稳斜段长度为。该段需满足在降斜段开始后能准确中靶。总位移目标=≈总垂深目标TV剩余垂深TV剩余位移De降斜段计算：由45度降至10度。降斜段长度Δ=降斜段垂深增量ΔT降斜段位移增量ΔD稳斜段长度求解：稳斜段垂深ΔT稳斜段位移ΔD校核：ta由于计算位移需求为272米，而按45度稳斜计算位移应为977米，此处需反推最大井斜角或调整剖面。修正设计：采用三段制（增-稳-降）无法满足小位移大垂深需求，或者最大井斜角过小。重新计算最大井斜角：利用几何关系，假设无稳斜段，直接增斜后降斜（双增型）。计算得所需最大井斜角约为。修正后设计：造斜率/30m，增斜至，然后降斜至。修正后增斜段长220m，垂深增107m，位移增修正后降斜段长120m，垂深增100m，位移增中间稳斜段需覆盖剩余垂深：32001500中间稳斜段位移：1493×总位移校核：43+略大于原定602米，需微调造斜点或造斜率。最终方案：将造斜点加深至1600米，重新计算后满足中靶要求。3.3靶点坐标计算与防碰扫描利用最小曲率法计算测段节点坐标。坐标转换公式：ΔΔΔ对邻井轨迹进行圆锥扫描法防碰计算，最近距离计算公式：D计算得与本区块邻井A-2的最近距离为45米，大于安全半径（8米），防碰风险低。轨迹设计关键节点数据表：节点井深井斜角方位角垂深北坐标东坐标造斜率工具面角井口0.000.000.000.000.000.000.00-造斜点1600.000.000.001600.000.000.000.00-增斜末1820.0022.0048.371816.5037.2041.803.00高边稳斜末3290.0022.0048.373210.50567.20637.800.00-降斜末3410.0010.0048.373310.00602.00676.00-3.00高边4.钻井水力学与井控计算水力学计算旨在优化水力参数，实现井底净化、井壁稳定及机械钻速最大化。井控计算则用于压井作业时的参数设定。4.1钻井液流变参数计算选用幂律流体模型描述钻井液流变性。使用范氏六速旋转粘度计读数：=45流性指数n：n稠度系数K：K该参数用于后续环空压耗及临界流速计算。4.2循环系统压耗计算循环系统总压耗Δ由地面管汇、钻铤内、钻杆内、钻铤外环空、钻杆外环空及钻头压降组成。采用宾汉塑性模型或幂律模型计算管流和环空流压耗。Δ其中，f为摩阻系数，需根据雷诺数Re判断流型并选取相应公式（层流f以井底段（钻铤外环空）为例：环空尺寸：钻铤外径O=178mm，井眼直径环空流速=。设排量Q==≈计算雷诺数R=计算得该段环空压耗约为2.5M累加各段压耗，求得总循环压耗Δ=4.3钻头水力参数优化采用最大水马力（MSP）或最大冲击力（MIF）工作方式。在限定地面泵压=22最优压降分配比（MSP准则）：Δ/故钻头压降Δ=剩余管路压耗Δ=根据钻头压降公式反推喷嘴总面积：Δ取流量系数=0.95计算得≈。推荐组合：三个喷嘴，直径分别为11mm,计算射流冲击力=ρQ和喷速=≈冲击力≈30004.4井控压井计算假设发生溢流，关井后读取立管压力=5MP溢流类型判定：通过关井套压与立压差值及钻井液增量判断为气体溢流。压井液密度计算：==采用司钻法压井：第一循环周：用原浆循环排出溢流，控制立压=+终了立压=0第二循环周：泵入重泥浆，控制立压从降至0。计算气体溢流在环空上返过程中的最大套压，需考虑气体滑脱上升及膨胀，采用哥德公式或数值模拟估算，最大套压预计出现在溢流顶面接近井口时，约为10.5M水力学参数汇总表：参数名称数值单位备注排量28L/s优选排量塑性粘度(PV)17mPa·s600-300读数差屈服值(YP)5.5Pa2倍300-600环空ECD1.28g/cm³含岩屑影响钻头压降14.3MPa占总压耗65%射流喷速105m/s清洗井底钻头水马力400kW有效破岩压井液密度1.41g/cm³平衡地层压力5.固井工程计算固井计算旨在确保水泥浆能有效封固套管，提供层间封隔和支撑。计算内容包括水泥浆量、替浆量、平衡压力及接触时间。5.1水泥浆体系设计及用量计算针对生产套管固井，采用领浆（稀释液）+尾浆（防气窜水泥浆）体系。井眼容积计算：考虑井径扩大率，平均井径=1.15环空容积=(对于3140-3180米（主力油层段），段长40米。=(水泥浆附加量：考虑井径不规则及损耗，附加系数取15%。总水泥浆量=×设计返高：返至技术套管鞋（2800米）。封固段长H=计算得需水泥浆总体积约12.5。5.2替浆及平衡液柱压力计算替浆量等于套管内替浆段的容积（扣除胶塞容积）。=需精确计算顶替到位时的压力平衡，防止“U型管”效应导致失返或蹩泵。计算静液柱压力差：环空液柱：水泥浆密度1.90g/cm³×380m+钻井液1.25g/cm³×(2800-380)。套管内液柱：钻井液1.25g/cm³×3180。由于水泥浆密度大于钻井液，会产生“U型管”效应，即水泥浆自动下行，导致泵压降低甚至出现真空（负压）。计算平衡点位置：设水泥浆在环空高度为h，在套管内高度为H。平衡时=。此计算用于指导施工中是否需要控制节流阀或调节排量，以保持流动连续性。5.3碰压及胶结强度计算碰压压力：=其中为胶塞在阻流环处的阻力，通常设计为5-10MPa。设计碰压压力不超过地层破裂压力的80%。胶结强度计算：水泥石与套管界面的胶结力τ。实验测得τ=可支撑的套管重量W=该数值用于评估水泥环对套管的支撑能力及层间封隔的耐久性。固井施工参数表：施工阶段流体类型密度用量排量接触时间累计时间注前置液冲洗液1.054.015-25100-10注领浆加砂水泥1.856.015-25-10-15注尾浆防窜水泥1.906.515-251515-22压胶塞钻井液1.251.010-20-22-24替泥浆钻井液1.2548.020-30-24-45碰压---5-10-45-476.完井与试油计算本节涵盖射孔参数优化、完井管柱受力分析及试油产能计算。6.1射孔参数优化计算采用负压射孔工艺，利用有限元法模拟射孔孔眼周围应力场及流动效率。负压值Δ计算：依据Bell经验公式，基于岩石声波时差及渗透率确定最小负压。Δ其中K为渗透率（50mD），为压实系数。计算得最小负压约为5MPa。设计射孔负压值为8MPa。孔深穿透计算：利用聚能射孔弹性能参数，计算在井底条件下的穿深。混凝土靶穿深=600mm贝雷砂岩校正系数CF井底有效穿深=600该穿深能有效穿透钻井污染带（预计深约300mm）。产能比计算（PR）：P其中为包表皮系数。优化孔密（16孔/米）及相位角（60度）后，计算得PR>6.2完井管柱强度与变形计算采用ϕ73mmN80外加厚油管作为生产管柱。螺旋弯曲计算：管柱在封隔器处受轴向压力，可能发生螺旋弯曲。临界弯曲载荷：=其中I为管柱截面惯性矩，Δr计算得≈15若管柱底部受到的浮重加活塞效应超过此值，则发生弯曲。计算中和点位置，评估管柱与套管的接触磨损情况。温度效应与压力效应（活塞效应）：生产过程中，温度升高导致管柱伸长：Δα为热膨胀系数，ΔT为温升（设为C）。Δ≈虎克效应（压力变化导致的长度变化）：Δ需综合计算管柱在封隔器处的总位移量，以选择合适的伸缩管或密封总成长度。6.3试油产能计算（IPR曲线）利用Vogel方程计算溶解气驱油藏的流入动态关系（IPR）。其中为平均地层压力（32MPa），为井底流压。进行一点法试油测试：实测流压=25MP代入方程求解绝对无阻流量（AOF）：=计算得AO据此绘制IPR曲线，确定合理生产压差及产量。完井及试油数据表：参数项目数值单位说明射孔枪型89mm深穿透聚能弹孔密16孔/米螺旋布孔相位角60度优化水力效率设计负压8.0MPa清洗孔眼临界弯曲载荷15kN防止永久变形预测AOF180t/d潜在产能推荐产量60t/d无水采油期7.增产措施（压裂）设计计算针对低渗透油藏（渗透率<10mD），进行水力压裂设计计算，以形成高导流能力的支撑裂缝。7.1裂缝几何尺寸计算（PKN模型）假设裂缝为垂直裂缝，且高度受限，采用Perkins-Kern-Nordgren(PKN)模型计算裂缝动态宽度。裂缝宽度w与净压力的关系：w简化计算最大缝宽：=其中，μ为压裂液粘度（200mPa·s），Q为排量（4m³/min），L为缝长（150m），为缝高（30m），为平面应变模量（30000MPa）。计算得最大动态缝宽≈6.5mm7.2支撑剂沉降与铺置计算选用陶粒作为支撑剂，粒径20/40目。单颗粒在压裂液中的沉降速度（Stokes定律修正）：=考虑干扰沉降，实际沉降速度约为自由沉降的0.5倍。计算得≈0.05计算停泵后支撑剂沉降距离，确保裂缝闭合时支撑剂仍位于产层内。裂缝闭合时间=。设计前置液占