井组注气提高采收率数值模拟研究

文72467井组注气提高采收率数值模拟研究摘要文南油田已动用的地质储量中，渗透率小于20103M2的低渗透储层的地质储量有1747104T，占全油田储量的2724。这部分储层注水困难、井况较差、且油井产量低。针对低渗透油藏，“注气开发”是国内外最为行之有效的开采方法；本文针对文南油田文72467注气井组的地质特征，利用多组分模型软件，完成了三维地质模型的建立，生产历史拟合和动态预测等。该油藏通过注天然气开采，能大幅度地提高该区块的原油采收率。关键词注气；提高采收率；数值模拟；低渗透；组分模型中图分类号TE文献标识码A1油田地质概况文南油田位于东濮凹陷中央隆起带文留构造南部次级地堑内，是一个异常高压、低渗、高饱和的复杂断块油气藏。主要含油层位为下第三系沙河街组沙二下、沙三上和沙三中亚段，共分28个砂层组。油藏埋藏深度22103800M之间。自1983年以来，陆续投入开发了8个断块区共26个开发单元，动用含油面积462KM2，动用石油地质储量6413104T，标定可采储量1694104T，标定采收率2642。文南油田已动用的地质储量中，渗透率大于100103M2的地质储量有1005104T占全油田储量的1567；大于50103M2的地质储量有1423104T，占全油田储量的215。渗透率在2050103M2的地质储量有2238104T，占全油田储量的349。渗透率小于20103M2的地质储量有1747104T，占全油田储量的2724。这部分储量注水难，见效难，井况差，油井产量低，同时这部分的潜力大。怎样开发这部分储量是个重大的科研课题。通过反复论证，最终确定以文72467井组作为作为其中一个先导试验井组。数值模拟研究的区块向东北方向延伸到文72453井，西南方向延伸到接近文72494井处，如图1所示。确定注气井为文72467井，注气层位沙三中79砂组，注气井段3534236846米，注气井段射开油层24层，厚度358米。图1文72467井组数值模拟研究的范围文72467井组位于文72块沙三中油藏的中部，井组主要开发层位为沙三中79砂组，其次为沙三中2、5、6、10砂组，原始地层压力为4766MPA，压力系数为138189之间，目前沙三中79砂组地层压力为4435MPA，井组包括一线井有文72461、文72463、文72471、文72483、文72453井，二线井包括文72491、文72473、文72495等井。2地层流体相态参数拟合PVT拟合的目的是用状态方程型相态分析软件对实验室所做的恒组分膨胀、多次接触脱气实验等数据进行计算拟合，得到可代表真实储层流体特性的状态方程参数1。这些参数是组分模型中油、气性质计算的基础，用于组分模拟研究27。本次研究利用多组分模型的相态模拟分析软件对文72467井组的高压PVT实验数据进行拟合计算。其井流物组份组成见表1，单次闪蒸实验拟合数据见表2。恒组成膨胀实验和拟合数据对比见表3、图2和图3。可以看出降压过程中体系和液相的相对体积变化拟合程度良好。表1文72467井原始井流物组份组成表组分名称摩尔组成组分名称摩尔组成N200011IC400092CO200024NC500127C106580IC500124C200693C600147C300371C701641NC400190表2文72467井单次闪蒸实验拟合数据参数实验拟合误差气油比M3/M34518451694000235地层原油密度G/CM305437055377185212地面脱气油密度G/CM308073080521258887饱和压力MPA4393439290000230000200400600801001201402025303540455055606570体体MPA体体体体体体体体体体体体图1恒组成膨胀实验体系相对体积拟合表3恒组成膨胀实验和拟合数据对比表体系相对体积液相相对体积压力（MPA）实验值拟合值实验值拟合值6535093370942510000100006276093940948010000100006080094400952310000100005884094890956910000100005688095470961610000100005492096010966610000100005296096640971910000100005099097270977410000100004803097940986310000100004707098690989410000100004511099480995910000100004393100001000010000100004388100031000399049971438410010100059791994743311004510037947596564219101301010890009142402510296102438223845836501071910563727675043438109831078868577068315511413111556343653628631197411642582860142561127191231254205472226213688132234910491320961440013884495545880002000400060008000100002025303540455055606570体体MPA体体体体体体体体体体体体图2恒组成膨胀实验液相相对体积拟合0000200400600801001201400510152025303540体体MPA体体体体体体体体体体体体图3等容衰竭实验气相压缩因子拟合等容衰竭实验和拟合数据对比见表4、图4和图5。可以看出衰竭过程中气相的压缩因子和液相的相对体积变化拟合程度良好。多次接触脱气实验和拟合数据对比见表5、图6和图7。可以看出脱气过程中体系的气油比和原油体积系数拟合程度良好。0002000400060008000100000510152025303540体体MPA体体体体体体体体体体体体图4等容衰竭实验液相相对体积拟合0001000020000300004000050000051015202530354045体体MPA体体体M3/M3体体体体体体图5多次接触脱气实验气油比拟合000050100150200250051015202530354045体体MPA体体体体体体体体体体体体图6多次接触脱气实验原油体积系数拟合表4等容衰竭实验和拟合数据对比表气相压缩因子液相相对体积压力（MPA）实验值拟合值实验值拟合值4393100201000010000038251029609123806582285323809612084727333735132652091550803567686743620590890307833632862353147108867079445964577419810900308322564653986591092240880854135094929409464092845251484431440961309564513046962表5多次接触脱气实验和拟合数据对比表气油比（M3/M3）原油体积系数压力（MPA）实验值拟合值实验值拟合值43934607480512436235734253141272271965179332382354246331786172426521929181631649155520591318130501541142214719138824145313129815835718137512305913523417131911682941591651127911193油藏模型及储量拟合根据地质研究成果，本次模拟研究的主要对象是砂三中7号层、砂三中8号层和砂三中9号层三个主力层位，其中砂三中7号层分为两个小层，砂三中8号层分为两个小层，砂三中9号层分为三个小层。为此，在数值模拟网格划分中，纵向上划分为七个网格，并形成一套变深度的网格系统。平面上网格方向基本与构造长轴保持一致，网格数目总计为6846721896。网格参数场是根据地质研究成果提供的地质模型参数经插值得到，其中有效渗透率参数场通过数值模拟历史拟合研究工作而被修正。模型中初始化形成的参数场是否能代表实际油气藏情况，储量是重要标志之一。多组分模型综合地质、测井和油藏工程分析资料，对某些不确定参数进行修正，计算的储量与原始地质储量等各项参数都有比较好的吻合，见表6。表6文72467井组数值模拟区域储量拟合原油地质储量参数层位（104T）（104M3）拟合储量（104M3）砂三中71603971986831986827砂三中8829681027721027719砂三中94079355053085053080合计65130080676380676264油田开发动态预测在完成油气藏地质模型的建立，流体PVT相态计算，油、气、水三相相对渗透率和毛细管压力的确定，储量拟合和历史拟合以后，就可以展开针对不同的开发方式，开发速度，开发井网和开发层系进行动态预测模拟工作。从数值模拟研究的任务出发，针对文南油田文72467注气井组，设计了20个注富化天然气驱开采方案；开始注气时间为2002年1月1日，动态预测至2015年12月31日。文72467注气井的日注气水平分别为2、3、4和5万方，注富化天然气的组成见表7。预测方案编号分别为F01F20，具体参数见表6，其模拟计算结果见表8。表7富化天然气的组成拟组分组成加入液化气摩尔百分数拟组分名称开封气液化气02468N2TOC190535004300905350887329869308851287833266CO2TOC2833201911008332085476876318978791942C3TONC4107207377001072025260397995433968878FC5TOFC600540669000054001867031940452205849C07TOC1700070000000007000069000670006600064C1800000000000000000000000000000000000对各方案的模拟计算结果进行综合分析之后，建议以方案F016的开发设计参数对文南油田文72467注气井组进行开发；在方案F016中，日注气量为5万方，注入气源为“开封气”。表8文72467井组数值模拟开发动态指标方案编号累积产油（104M3）采出程度（）增加富化度增产油量（104M3）增加注气量增产油量（104M3）F00121278242637483F00221374482649413009624F00321469722661218019148F00421556442671967027820F00521616432679403033819F00624040492979870276225F00724105822987968006533273134F00824188602998229014811271888F00924285073010187024458272863F01024397583024132035709278115F01126315293261837503705F01226453303278943013801507882F01326604513297686028922513479F01426680993307166036570512455F01526821773324615050648520534F01627740593438506646235F01728099523482995035893672504F01828270463504184052987680074F01928431633524161069104687519F020285998735450160859286983445数值模拟研究结论针对低渗透（和特低渗透）油藏，“注气开发”是国内外最为行之有效的开采方法；文南油田文72467注气井组油藏注富化天然气驱开采，能大幅度地提高该区块的原油采收率；数值模拟研究结论表明随着注气强度（或注气量）的增加，累积产油量和原油的采出程度均会有较大幅度的增加，即注气强度（或注气量）是最为敏感的控制参数；数值模拟研究结论还表明随着注入天然气富化度的增加，累积产油量和原油的采出程度只有较小幅度的增加，即注入天然气富化度不是敏感的控制参数。参考文献1张茂林，俞高明油气体系拟组分相平衡及物性参数计算J天然气工业，1988,8（4）26322张茂林，吴清松等毛细管压力影响油气体系相平衡的理论研究与应用J天然气工业，2000，20（5）53563张茂林，俞高明等凝析油气藏拟组分数值模拟方法J西南石油学院学报，1991，13（2）38474梅海燕，孔祥言等油气烃类体系气液固三相相平衡计算J天然气工业，2000，20（3）75785梅海燕，孔祥言等预测石蜡沉积的热力学模型J石油勘探与开发，2000，27（1）84866张茂林，孙良田等凝析油气藏K值多组分模型数值模拟方法J石油学报，1991，12（1）60667张茂林，孙良田等分数步长算法在K值多组分模型中的应用J天然气工业，1991，11（4）5156（编辑）文72467井组注气提高采收率数值模拟研究摘要文南油田已动用的地质储量中，渗透率小于20103M2的低渗透储层的地质储量有1747104T，占全油田储量的2724。这部分储层注水困难、井况较差、且油井产量低。针对低渗透油藏，“注气开发”是国内外最为行之有效的开采方法；本文针对文南油田文72467注气井组的地质特征，利用多组分模型软件，完成了三维地质模型的建立，生产历史拟合和动态预测等。该油藏通过注天然气开采，能大幅度地提高该区块的原油采收率。关键词注气；提高采收率；数值模拟；低渗透；组分模型THENUMERICALSIMULATIONRESEARCHOFENHANCINGOILRECOVERYBYGASINJECTIONONWEN72467WELLGROUPABSTRACTTHELOWPERMEABLESTRATUM,INWHICHTHEPERMEABILITYI