CN112036033A - 一种定量确定蒸汽吞吐转蒸汽驱地层温度时机的方法 - Google Patents
一种定量确定蒸汽吞吐转蒸汽驱地层温度时机的方法 Download PDFInfo
- Publication number
- CN112036033A CN112036033A CN202010894337.7A CN202010894337A CN112036033A CN 112036033 A CN112036033 A CN 112036033A CN 202010894337 A CN202010894337 A CN 202010894337A CN 112036033 A CN112036033 A CN 112036033A
- Authority
- CN
- China
- Prior art keywords
- temperature
- injection
- steam
- production
- different
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 230000015572 biosynthetic process Effects 0.000 title claims abstract description 54
- 238000010795 Steam Flooding Methods 0.000 title claims abstract description 40
- 238000000034 method Methods 0.000 title claims abstract description 26
- 238000004519 manufacturing process Methods 0.000 claims abstract description 111
- 239000003921 oil Substances 0.000 claims abstract description 84
- 238000002347 injection Methods 0.000 claims abstract description 60
- 239000007924 injection Substances 0.000 claims abstract description 60
- 239000010779 crude oil Substances 0.000 claims abstract description 55
- 239000000295 fuel oil Substances 0.000 claims abstract description 26
- 238000010586 diagram Methods 0.000 claims abstract description 14
- 238000002474 experimental method Methods 0.000 claims abstract description 9
- 238000006073 displacement reaction Methods 0.000 claims description 34
- 238000010793 Steam injection (oil industry) Methods 0.000 claims description 16
- 230000000694 effects Effects 0.000 claims description 13
- 239000012530 fluid Substances 0.000 claims description 13
- 230000035699 permeability Effects 0.000 claims description 12
- 230000008859 change Effects 0.000 claims description 10
- 238000012360 testing method Methods 0.000 claims description 9
- 238000004364 calculation method Methods 0.000 claims description 8
- 230000000638 stimulation Effects 0.000 claims description 8
- 238000012546 transfer Methods 0.000 claims description 7
- 230000000977 initiatory effect Effects 0.000 claims description 5
- 238000000518 rheometry Methods 0.000 claims description 5
- 238000007872 degassing Methods 0.000 claims description 4
- 239000011148 porous material Substances 0.000 claims description 4
- 238000011084 recovery Methods 0.000 description 21
- 238000004891 communication Methods 0.000 description 10
- 238000010438 heat treatment Methods 0.000 description 9
- 238000011161 development Methods 0.000 description 8
- 230000007423 decrease Effects 0.000 description 5
- 238000011160 research Methods 0.000 description 4
- 230000008569 process Effects 0.000 description 3
- IJGRMHOSHXDMSA-UHFFFAOYSA-N Atomic nitrogen Chemical compound N#N IJGRMHOSHXDMSA-UHFFFAOYSA-N 0.000 description 2
- CURLTUGMZLYLDI-UHFFFAOYSA-N Carbon dioxide Chemical compound O=C=O CURLTUGMZLYLDI-UHFFFAOYSA-N 0.000 description 2
- 241000336847 Luda Species 0.000 description 2
- 230000002411 adverse Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 2
- 230000005465 channeling Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000011002 quantification Methods 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- DIWRORZWFLOCLC-UHFFFAOYSA-N Lorazepam Chemical compound C12=CC(Cl)=CC=C2NC(=O)C(O)N=C1C1=CC=CC=C1Cl DIWRORZWFLOCLC-UHFFFAOYSA-N 0.000 description 1
- 238000010796 Steam-assisted gravity drainage Methods 0.000 description 1
- 229910002092 carbon dioxide Inorganic materials 0.000 description 1
- 239000001569 carbon dioxide Substances 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000008021 deposition Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 230000004069 differentiation Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 229910052757 nitrogen Inorganic materials 0.000 description 1
- 230000002028 premature Effects 0.000 description 1
- 238000013139 quantization Methods 0.000 description 1
- 239000004576 sand Substances 0.000 description 1
- 238000002791 soaking Methods 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B43/00—Methods or apparatus for obtaining oil, gas, water, soluble or meltable materials or a slurry of minerals from wells
- E21B43/16—Enhanced recovery methods for obtaining hydrocarbons
- E21B43/24—Enhanced recovery methods for obtaining hydrocarbons using heat, e.g. steam injection
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/08—Thermal analysis or thermal optimisation
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A10/00—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE at coastal zones; at river basins
- Y02A10/40—Controlling or monitoring, e.g. of flood or hurricane; Forecasting, e.g. risk assessment or mapping
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Physics (AREA)
- Geology (AREA)
- Mathematical Analysis (AREA)
- Mining & Mineral Resources (AREA)
- Pure & Applied Mathematics (AREA)
- Computational Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Mathematical Optimization (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Engineering & Computer Science (AREA)
- Databases & Information Systems (AREA)
- Algebra (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- Software Systems (AREA)
- Environmental & Geological Engineering (AREA)
- Fluid Mechanics (AREA)
- Computer Hardware Design (AREA)
- Operations Research (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geochemistry & Mineralogy (AREA)
- Geophysics And Detection Of Objects (AREA)
- Fats And Perfumes (AREA)
- Production Of Liquid Hydrocarbon Mixture For Refining Petroleum (AREA)
Abstract
本发明公开了一种定量确定蒸汽吞吐转蒸汽驱地层温度时机的方法:收集目标油藏基础数据资料,包括:原始油藏温度、地层条件下原油粘度、地面脱气原油的粘温曲线;开展稠油样品流变性实验;回归不同温度下的流变性方程;确定稠油样品流变性拐点温度;确定拐点温度对应的原油粘度;绘制多轮蒸汽吞吐后地层温度分布图;确定启动压力梯度与流度关系;确定注采井间的启动压力梯度;绘制不同注采井距对应的转驱温度图版;判断是否达到转驱的温度时机。本发明可对不同稠油油藏、不同注采条件下的蒸汽吞吐转蒸汽驱地层温度时机进行定量预测,能解决目前行业中关于稠油油藏多轮吞吐后转蒸汽驱的温度时机定性描述“注采井间形成热连通”缺乏可操作性的难题。
Description
技术领域
本发明属于石油开采技术领域,更具体地说,是涉及稠油油藏一种定量确定蒸汽吞吐转蒸汽驱地层温度时机的方法。
背景技术
热力采油技术(如蒸汽吞吐、蒸汽驱、蒸汽辅助重力泄油等)是一项可以大幅度提高稠油油田采收率的有效技术,是目前陆上稠油开发的有效手段。通常蒸汽吞吐的加热区域较小,井间动用程度低,导致蒸汽吞吐采收率仅为20%左右;而蒸汽吞吐后转蒸汽驱可以增加加热区域,提高井间的动用程度,采收率可达到40%左右,是蒸汽吞吐后进一步提高后续采收率的重要方式。
研究和实践均表明,蒸汽驱阶段的效果取决于在蒸汽吞吐阶段注汽井与生产井之间是否实现热连通,即井间的地层是否被注入蒸汽加热。若在蒸汽吞吐阶段注汽井与生产井之间没有实现热连通,即井间地层温度没有升高,仍然为原始地层温度,则对蒸汽驱产生不利影响;而蒸汽吞吐过度,会造成注汽井与生产井之间形成汽窜通道,导致蒸汽驱过早汽窜和平面驱替不均匀,从而影响蒸汽驱效果。因此,定量确定蒸汽吞吐转蒸汽驱的温度时机,对提高蒸汽吞吐后转蒸汽驱的效果至关重要。
关于稠油油藏多轮蒸汽吞吐后的转驱温度时机,目前行业的表述为“注采井间形成热连通”。在实际应用中,存在诸多问题,如:(1)“热连通”是一个定性判断,注采井间的温度上升到多少适合转蒸汽驱需要定量确定,即转驱温度T转驱需要定量化。陆地油田稠油热采大都采用小井距(70m-120m)开发模式,多轮次蒸汽吞吐后加热半径可以超过井距中心位置(RA+RB≥L),但中心位置的温度缺乏定量判断(图1)。(2)对于不同粘度的稠油油藏,转驱温度T转驱应该差异化。原油粘度不同,转驱时要求注采井间达到的温度条件应当有所不同。海上稠油热采由于经济性要求更高,通常采用大井距(200m-300m)开发模式。研究和实践表明,经过多轮蒸汽吞吐的加热半径仅为50m左右,按照目前的行业定性描述,海上稠油始终无法达到“注采井间形成热连通”的条件((RA+RB<L,图2),而这与油田实际不符。可见,目前关于转蒸汽驱温度时机的定性描述,缺乏定量化和差异化,判断模糊,在实际油田的研究中难以操作,对蒸汽吞吐后转蒸汽驱的研究带来了不利影响。
发明内容
本发明的目的是为了克服现有技术中的不足,提出一种定量确定蒸汽吞吐转蒸汽驱地层温度时机的方法,以此解决目前行业中关于稠油油藏多轮吞吐后转蒸汽驱的温度时机定性描述“注采井间形成热连通”缺乏可操作性,难满足矿场实际需要的难题。
本发明的目的是通过以下技术方案实现的。
本发明定量确定蒸汽吞吐转蒸汽驱地层温度时机的方法,包括以下步骤:
第一步,收集目标油藏基础数据资料,包括:原始油藏温度、地层条件下原油粘度、地面脱气原油的粘温曲线;
第二步,开展稠油样品流变性实验
提取目标油藏的稠油样品,脱气脱水后测试至少7个温度条件下的原油流变性,每个温度下测试剪切速率从1s-1到50s-1,间隔2s-1,共25个剪切速率下的原油粘度,同时根据公式(A-1)计算对应的屈服值τ;
τ=μoν/1000 (A-1)
式中,τ表示屈服值,Pa;μo表示原油粘度,mPa·s;v表示剪切速率,1s-1;
第三步,回归不同温度下的流变性方程
在直角坐标系中,以剪切速率v为横坐标,屈服值τ为纵坐标,分别回归第二步中7个温度条件的流变性方程;特点是不同温度条件下,流变性方程的形式是相同的,均为τ=Av+B的线性形式,只是不同温度条件下的斜率A值和常数项B值不同;
第四步,确定稠油样品流变性拐点温度
观察第三步得到的不同温度下的流变性方程,随着温度的增加,常数项B逐渐降低;当温度为T0时,常数项B约等于零,此时对应的温度T0即为稠油样品从非牛顿流体转变为牛顿流体的拐点温度T拐点;
第五步,确定拐点温度对应的原油粘度
在直角坐标上,以温度lgT为横坐标,原油粘度lgμo为纵坐标绘制散点图,回归得到(A-2-1)形式的粘温关系式,其中C和D为回归得到的系数,不同的稠油样品回归得到不同的C和D值;
lgμo=C+Dlg T (A-2-1)
式中,μo为原油粘度,mPa·s;T为温度,℃;得到原油粘度计算公式:
μo=10CTD (A-2-2)
第六步,绘制多轮蒸汽吞吐后地层温度分布图
在直角坐标系中,以距离注蒸汽井的距离为横坐标,以网格平均温度为纵坐标,绘制蒸汽吞吐N年后的网格平均温度随注蒸汽井的距离变化曲线,根据曲线确定注采井中间位置地层温度;如果注采井中间位置地层温度高于或等于拐点温度T中心≥T拐点,不考虑启动压力梯度的影响,认为达到转驱温度条件转蒸汽驱;如果注采井中间位置地层温度低于拐点温度T中心<T拐点,考虑启动压力梯度对能否建立有效注采***的影响;
第七步,确定启动压力梯度与流度关系
开展启动压力梯度实验,回归启动压力梯度与流度关系式。启动压力梯度与流体的性质、介质的表面作用和孔隙结构密切相关,不同油藏的启动压力梯度具有统一的公式样式(A-3-1):
G=E(K/μo)F (A-3-1)
式中,G表示启动压力梯度,MPa/m;K表示地层渗透率,10-3μm2,μo表示地层条件下的原油粘度,mPa·s;E和F为回归得到的系数,不同的油藏可回归得到不同的值;
第八步,确定注采井间的驱替压力梯度
根据最大注采压差ΔP,注采井间的井距L,确定注采井间中心位置的驱替压力梯度G驱替=ΔP/L;根据第七步中公式(A-3-1)计算启动压力梯度G启动;比较驱替压力梯度G驱替和启动压力梯度G启动的大小;
第九步,绘制不同注采井距对应的转驱温度图版
联立方程(A-2-2)和(A-3-1),得到:
对于等式右边,将地层渗透率K和第五步、第七步回归得到的系数C、D、E、F值代入,右边只剩下一个变量温度T;对于等式左边,在不同最大注采压差ΔP的情况下,分别代入不同注采井间的井距L,通过(A-4-1)计算对应的温度T;在直角坐标系中,以注采井间的井距L为横坐标,对应的温度T为纵坐标,绘制不同注采压差下不同注采井距的井间转驱温度图版;
第十步,判断是否达到转驱的温度时机
对于目标油田,根据最大注采压差ΔP,注采井间的井距L,查阅第九步的图版,确定转驱注采井中心位置所需的温度T转驱,与注采井中心位置地层温度T地层进行对比,当T地层≥T转驱时,可转蒸汽驱,当T地层<T转驱时,则不能转蒸汽驱。
与现有技术相比,本发明的技术方案所带来的有益效果是:
本发明从地层渗透率、地层原油粘度、油品的粘温特征等基本数据出发,研究了不同油品流变性拐点温度、启动压力梯度等影响热采过程渗流的关键参数,提出了只要注采井间的驱替压力梯度G驱替高于启动压力梯度就可以进行转驱,可定量得到不同稠油油藏蒸汽吞吐后转驱时注采井间的温度时机。
本发明提供的技术方法,改变了目前行业现行的多轮吞吐转驱时机——“注采井间达到热连通”这一非定量、非差异化表述。通过该方法,可以定量得到不同稠油油藏蒸汽吞吐后转驱时注采井间的温度时机。
附图说明
图1陆地油田小井距蒸汽吞吐后转蒸汽驱热采示意图
图2海上油田大井距蒸汽吞吐后转蒸汽驱热采示意图
图3为确定稠油油藏热采多轮吞吐后转驱地层温度时机方法的主要步骤
图4为LD27-2-6井在温度为50℃条件下的屈服值随剪切速率变化曲线
图5为LD27-2-6井在温度为75℃条件下的屈服值随剪切速率变化曲线
图6为LD27-2-6井原油屈服值随温度变化曲线
图7为LD27-2油田原油粘度随温度的变化曲线
图8为LD27-2油田蒸汽吞吐网格平均温度随距注汽井距离的分布曲线
图9为LD27-2油田不同注采压差下不同注采井距的井间转驱温度图版
图10为LD27-2油田蒸汽吞吐15轮次后的温度场图(左:目前井网;右:中间加密1口井进行蒸汽吞吐)
图11为LD27-2油田蒸汽吞吐后转蒸汽驱的温度场图(左:目前井网;右:中间加密1口井转注)。
具体实施方式
下面结合附图对本发明作进一步的描述。
本发明定量确定蒸汽吞吐转蒸汽驱地层温度时机的方法,具体为一种定量确定稠油油藏蒸汽吞吐(包括多元热流体吞吐、氮气辅助蒸汽吞吐、二氧化碳辅助蒸汽吞吐等)转蒸汽驱地层温度时机的方法,基于流变性拐点温度、启动压力梯度等稠油渗流特征参数,定量确定了不同稠油油藏蒸汽吞吐后转蒸汽驱时注采井间所需要达到的地层温度值,具有可量化、针对性强的特点,主要包括以下步骤,如图3所示。
第一步,收集目标油藏基础数据资料,
收集目标油藏的原始油藏温度、地层条件下原油粘度、地面脱气原油的粘温曲线等基础数据。
第二步,开展稠油样品流变性实验
提取目标油藏的稠油样品,脱气脱水后测试至少7个温度条件下的原油流变性,每个温度下测试剪切速率从1s-1到50s-1,间隔2s-1,共25个剪切速率下的原油粘度,同时根据公式(A-1)计算对应的屈服值τ。
τ=μoν/1000 (A-1)
式中,τ表示屈服值,Pa;μo表示原油粘度,mPa·s;ν表示剪切速率,1s-1。如50℃测得的不同剪切速率下的原油粘度变化值,根据公式(A-1)计算屈服值τ,计算结果见表1。
表1 LD27-2-6井流变性实验数据
第三步,回归不同温度下的流变性方程
在直角坐标系中,以剪切速率ν为横坐标,屈服值τ为纵坐标,分别回归第二步中7个温度条件的流变性方程。该方程的特点是不同温度条件下,流变性方程的形式是相同的,均为τ=Av+B的线性形式,只是不同温度条件下的斜率A值和常数项B值不同。
图4和图5分别为50℃和75℃条件下的屈服值随剪切速率变化散点图,对50℃和75℃回归结果见表2中第1行和第4行,第三列R2为回归方程的关联系数。按照同样的方法,可以分别绘制60℃、70℃、80℃、90℃、100℃等温度下的散点图,并回归不同温度条件的流变性方程,结果见表2。
表2 LD27-2-6井不同温度条件下的流变方程
第四步,确定稠油样品流变性拐点温度
观察第三步得到的不同温度下的流变性方程。随着温度的增加,常数项B逐渐降低;当温度为T0时,常数项B约等于零,此时对应的温度T0即为稠油样品从非牛顿流体转变为牛顿流体的拐点温度T拐点。
由表2可知,随着温度的增加,常数项B逐渐降低;当温度为75℃时,常数项B趋于0.002(约等于零),温度高于75℃时,B值略有波动,但都有一个共同特征:约等于零。因此LD27-2-6井稠油样品的拐点温度为75℃;为更好识别拐点温度,可绘制的不同温度下的屈服值变化曲线(见图6),可明显看出75℃时屈服值趋于零。
第五步,确定拐点温度对应的原油粘度
在直角坐标上,以温度lgT为横坐标,原油粘度lgμo为纵坐标绘制散点图,可回归得到(A-2-1)形式的粘温关系式,其中C和D为回归得到的系数,不同的稠油样品回归得到不同的C和D值。
lgμo=C+Dlg T (A-2-1)
式中,μo为原油粘度,mPa·s;T为温度,℃。可得到原油粘度计算公式:
μo=10CTD (A-2-2)
如图7所示,根据LD 27-2油田的稠油粘温数据,回归得到的粘温关系式为(A-2-3),C和D值分别为9.1168和-3.5405。
lgμo=9.1168-3.5405lg T (A-2-3)
根据(A-2-2)可得到原油粘度计算公式:
μo=1.3085×109T-3.5405 (A-2-4)
第六步,绘制多轮蒸汽吞吐后地层温度分布图
在直角坐标系中,以距离注蒸汽井的距离为横坐标,以网格平均温度为纵坐标,绘制蒸汽吞吐N年后的网格平均温度随注蒸汽井的距离变化曲线,根据曲线确定注采井中间位置地层温度。如果蒸汽吞吐阶段设计为5年,则绘制吞吐5年后的曲线;如果蒸汽吞吐阶段设计为10年,则绘制吞吐10年后的曲线。在加热半径范围内,地层温度T高于原始油藏温度T原始。如果注采井中间位置地层温度高于或等于拐点温度T中心≥T拐点,可不考虑启动压力梯度的影响,可认为达到转驱温度条件转蒸汽驱;如果注采井中间位置地层温度低于拐点温度T中心<T拐点,需要考虑启动压力梯度对能否建立有效注采***的影响。
如图8所示,绘制的LD27-2油田不同轮次蒸汽吞吐网格平均温度随注蒸汽井的距离变化曲线,其特征在于,随着距离注汽井距离的增加,网格的温度逐渐降低,在50m左右的时候,网格的温度为原始油层的温度。因为注采井距中心位置为距离注汽井200m(距离注汽井距离L=200m),吞吐10年中心位置的地层温度T中心仍然为原始地层温度50℃,小于第四步确定的拐点温度75℃,因此在蒸汽吞吐阶段难以达到热连通,此时需要考虑启动压力梯度是否能够建立有效注采***。
第七步,确定启动压力梯度与流度关系
开展启动压力梯度实验,回归启动压力梯度与流度关系式。启动压力梯度与流体的性质、介质的表面作用和孔隙结构密切相关,不同油藏的启动压力梯度具有统一的公式样式(A-3-1),只是不同油藏的回归系数E和F有区别。
G=E(K/μo)F (A-3-1)
式中,G表示启动压力梯度,MPa/m;K表示地层渗透率,10-3μm2;μo表示地层条件下的原油粘度,mPa·s;E和F为回归得到的系数,不同的油藏可回归得到不同的值。
稠油和低渗油藏相似,在多孔介质条件下,为具有启动压力的非达西渗流规律,稀油低渗透油藏启动压力梯度主要由低渗透率造成,而稠油的启动压力梯度则是由于原油粘度高造成,随粘度增加,启动压力梯度升高。
例如,根据渤海海域3个稠油油田油样实验室数据回归得到的启动压力梯度关系式见(A-3-2),系数E为0.1213,系数F为-0.8688。
G=0.1213(K/μo)-0.8688 (A-3-2)
上述方程(A-3-2)的特征在于不同区块或者不同油藏需要根据各自的实验数据进行回归,得到不同的系数E和F的值,从而得到启动压力梯度和流度的关系式。
第八步,确定注采井间的驱替压力梯度
根据最大注采压差ΔP,注采井间的井距L,确定注采井间中心位置的驱替压力梯度G驱替=ΔP/L;根据第七步中公式(A-3-1)计算启动压力梯度G启动;比较驱替压力梯度G驱替和启动压力梯度G启动的大小。其中,启动压力梯度G启动为中心位置处原油流动所需要的最小压力梯度;驱替压力梯度G驱替为油田实际注采过程中的实际驱替压力梯度;只有当G驱替≥G启动时中心位置处的原油才能流动。
第九步,绘制不同注采井距对应的转驱温度图版
联立方程(A-2-2)和(A-3-1),得到:
对于等式右边,将地层渗透率K和第五步、第七步回归得到的系数C、D、E、F值代入,右边只剩下一个变量温度T;对于等式左边,在不同最大注采压差ΔP(如ΔP=3MPa、4MPa、5MPa、6MPa等)的情况下,分别代入不同注采井间的井距L(如L=60m、80m、100m……360m、380m、400m等),通过(A-4-1)计算对应的温度T。在直角坐标系中,以注采井间的井距L为横坐标,对应的温度T为纵坐标,绘制不同注采压差下不同注采井距的井间转驱温度图版。
如,LD27-2油田得到的C、D、E、F值分别为,0.9598、-3.5405、0.1213、-0.8688,地层渗透率K为3785×10-3μm2,取ΔP=3MPa,通过(A-4-1)分别计算不同注采井距L=60m、80m、100m……360m、380m、400m的值,在直角坐标系中,以注采井距L为横坐标,对应的温度T为纵坐标,绘制成一条曲线;同理,按照相同的办法,可绘制出ΔP=4MPa、5MPa、6MPa等不同注采压差情况下的曲线,最终形成如图9的不同注采压差下不同注采井距的井间转驱温度图版。
第十步,判断是否达到转驱的温度时机
对于目标油田,根据最大注采压差ΔP,注采井间的井距L,查阅第九步的图版,确定转驱注采井中心位置所需的温度T转驱,与注采井中心位置地层温度T地层进行对比,当T地层≥T转驱时,可转蒸汽驱,当T地层<T转驱时,则不能转蒸汽驱。
其主要原理为在满足地层压力下降到3~5MPa的前提下,判断蒸汽吞吐结束后两口井中间位置处的驱替压力梯度G驱替是否大于启动压力梯度G启动。如果G驱替≥G启动,说明在目前地层温度下注采井间可建立有效驱替***,可认为满足了“井间达到热连通”的条件,可以转蒸汽驱;如果G驱替<G启动,说明在此温度条件下井间原油流动存在较大的启动压力梯度,未能建立有效驱替***,可认为没有满足“井间达到热连通”的条件,不能转蒸汽驱。
实施例:
本发明在LD27-2油田吞吐后转蒸汽驱的方案设计中得到了应用,应用实例表明本发明提供的方法有利于定量判断不同稠油油藏的转驱温度时机。
LD27-2油田主要含油层系为明化镇组,属于受断层控制的低幅鼻状构造,以河道、滩坝型浅水三角洲沉积为主,储层横向变化较大,非均质性强,属于高孔高渗储层。地层原油黏度2336.00mPa·s,属于普Ⅱ类稠油。明化镇组试采井冷采效果差,水平井冷采初期产能20~30m3/d,预测采收率仅7.0%,为提高采收率编制了总体热采开发方案,并在主力砂体部署热采井7口。为确保试验效果、降低试验风险,优选明化镇组两口井(A22H、A23H)进行蒸汽吞吐先导试验。两口井于2013年12月投产,截至2019年12月,A22H和A23H井累产油9.65×104m3,井控储量采出程度10.7%。
以LD27-2油田为例,展示该发明的实施过程。实施过程包含以下10个步骤:
第一步,目标油藏基础资料收集。收集旅大27-2油田的原始油藏温度、地层条件下原油粘度、地面脱气原油的粘温曲线等基础数据。
第二步,开展稠油样品流变性实验。提取旅大27-2油田的油样,脱气脱水后测试50℃~100℃之间7个温度条件下的原油流变性,每个温度下测试剪切速率从1s-1到50s-1(间隔2s-1)共25个剪切速率下的原油粘度,如50℃测得的不同剪切速率下的原油粘度变化值见表1中的第2列;根据公式(A-1)计算屈服值τ,计算结果见表1中的第3列。
τ=μov/1000 (A-1)
上式中的τ表示屈服值,Pa;μo表示原油粘度,mPa·s;v表示剪切速率,1s-1。
第三步,回归不同温度下的流变方程。在直角坐标系中,以剪切速率v为横坐标(表1中第1列),屈服值τ为纵坐标(表1中第3列),分别回归第二步7个温度条件的流变方程。图4和图5分别为50℃和75℃条件下的屈服值随剪切速率变化散点图,对50℃和75℃回归结果见表2中第1行和第4行,第三例R2为回归方程的关联系数。按照同样的方法,可以分别绘制60℃、70℃、80℃、90℃、100℃等温度下的散点图,并回归不同温度条件的流变性方程,结果见表2。
第四步,确定稠油样品流变性拐点温度。根据第三步,观察不同温度下的流变方程(表2中第2列)。随着温度的增加,常数项B逐渐降低;当温度为T0时,常数项B约等于零,此时对应的温度T0即为稠油从非牛顿流体转变为牛顿流体的拐点温度T拐点。由表2可知,随着温度的增加,常数项B逐渐降低;当温度为75℃时,常数项B趋于0.002(约等于零),温度高于75℃时,B值略有波动,但都有一个共同特征:约等于零。因此LD27-2-6井稠油样品的拐点温度为75℃;为更好识别拐点温度,可绘制的不同温度下的屈服值变化曲线(见图6),可明显看出75℃时屈服值趋于零。
第五步,确定拐点温度对应的稠油粘度。在直角坐标上,以温度lgT为横坐标,原油粘度lgμo为纵坐标绘制散点图。可回归得到(A-2-1)形式的粘温关系式,其中C和D为回归得到的系数,不同的原油样品可回归得到不同的C和D值。
lgμo=C+Dlg T (A-2-1)式中,μo为原油粘度,mPa·s;T为温度,℃。可得到原油粘度计算公式:
μo=10CTD (A-2-2)
如图7所示,根据LD 27-2油田的稠油粘温数据,回归得到的粘温关系式为(A-2-3),C和D值分别为9.1168和-3.5405。
lgμo=9.1168-3.5405lg T (A-2-3)
可得到原油粘度计算公式:
μo=1.3085×109T-3.5405 (A-2-4)
第六步,绘制多轮蒸汽吞吐后地层温度分布图。根据数值模拟结果,在直角坐标系中,以距离注蒸汽井的距离为横坐标,以网格平均温度为纵坐标,绘制距离井筒不同距离处的地层温度。如图8所示,绘制的LD27-2油田不同轮次蒸汽吞吐网格平均温度随注蒸汽井的距离变化曲线,其特征在于,随着距离注汽井距离的增加,网格的温度逐渐降低,在50m左右的时候,网格的温度为原始油层的温度。因为注采井距中心位置为距离注汽井200m,吞吐10年中心位置的地层温度T中心仍然为原始地层温度50℃,小于第四步确定的拐点温度75℃,因此在蒸汽吞吐阶段难以达到热连通,此时需要考虑启动压力梯度是否能够建立有效注采***。
LD27-2油田目前2口蒸汽吞吐井加热半径小,第3轮焖井结束时,A22H井加热半径28m左右,A23H井加热半径35m左右。数值模拟预测该油田蒸汽吞吐开发15轮次已经达到蒸汽吞吐开发的极限,水平井蒸汽吞吐的加热半径为50m,见图10。
第七步,确定启动压力梯度与流度关系。开展启动压力梯度实验,回归启动压力梯度与流度关系式。启动压力梯度与流体的性质、介质的表面作用和孔隙结构密切相关,不同油藏的启动压力梯度具有统一的公式样式(A-3-1)。
G=E(K/μo)F (A-3-1)
上式中的G表示启动压力梯度,MPa/m;K表示地层渗透率10-3μm2,μo表示地层条件下的原油粘度,mPa·s;E和F为回归得到的系数,不同的油藏可回归得到不同的值。例如,根据渤海海域3个稠油油田油样实验室数据回归得到的启动压力梯度关系式见(A-3-2),系数E为0.1213,系数F为-0.8688。
G=0.1213(K/μo)-0.8688 (A-3-2)
上述方程(A-3-2)的特征在于不同区块或者不同油藏需要根据各自的实验数据进行回归,得到不同的系数E和F的值,从而得到启动压力梯度和流度的关系式。
第八步,确定注采井间的驱替压力梯度。根据最大生产压差ΔP,注采井间的井距L,确定注采井间中心位置的驱替压力梯度G驱替=ΔP/L。根据第七步(A-3-1)计算启动压力梯度G启动。根据LD27-2油田实际情况,最大生产压差ΔP为6MPa,注采井间的井距L为400m,确定注采井间中心位置的驱替压力梯度G驱替=ΔP/L=0.015MPa/m。旅大27-2油田流度K/μo为1.62,根据第七步(A-3-1)计算启动压力梯度G启动=0.08MPa/m。
第九步,绘制不同注采井距对应的转驱温度图版。联立方程(A-2-2)和(A-3-1),得到:
对于等式右边,将地层渗透率K和第五步和第七步回归得到的系数C、D、E、F值代入,LD27-2油田得到的C、D、E、F值分别为,0.9598、-3.5405、0.1213、-0.8688,地层渗透率K为3785×10-3μm2,取ΔP=3MPa,通过(A-4-1)分别计算不同注采井距L=60m、80m、100m……360m、380m、400m的值,在直角坐标系中,以注采井距L为横坐标,对应的温度T为纵坐标,绘制成一条曲线;同理,按照相同的办法,可绘制出ΔP=4MPa、5MPa、6MPa等不同注采压差情况下的曲线,最终形成如图9的不同注采压差下不同注采井距的井间转驱温度图版。
第十步,判断是否达到转驱的温度时机。对于LD27-2油田,根据最大注采压差ΔP,注采井间的井距L,查阅第九步的图版(图9),查阅转驱注采井中心位置所需的温度T转驱。然后对比注采井中心位置地层温度T地层,当T地层≥T转驱时,可转蒸汽驱;当T地层<T转驱时,则不能转蒸汽驱。根据LD27-2油田实际情况,最大注采压差ΔP=6MPa,目前井网的注采井距L=400m,查阅第九步的图版(图9),确定注采井间中心位置的转驱所需要的温度为T转驱=70℃,而在蒸汽吞吐结束后,中心位置的地层温度为原始油藏温度T地层=55℃,T地层<T转驱,显然不能满足转驱的温度条件。
为了进一步提高蒸汽吞吐后的采收率,转驱是必然的选择。为此,可以考虑在注采井间加密一口井,在注采压差仍然为ΔP=6MPa的情况下,加密井网的注采井距L=200m,查阅第九步的图版(图9),确定注采井间中心位置的转驱所需要的温度为T转驱=56℃,而在蒸汽吞吐结束后,中心位置的地层温度为原始油藏温度T地层=55℃,T地层≈T转驱,此时可以满足转驱的温度条件。
为进一步论证本发明的有效性,可进一步研究不同井网蒸汽吞吐、多轮吞吐后转驱的指标见表3,按照目前400m井距预测吞吐15轮次(方案1)的动用储量采收率仅为23.8%,主要原因在于吞吐开发方式动用的储层半径较小,而实际油田的井距较大。为此,如图11,在两口井中间加密一口井H1,如果H1吞吐(方案2),采收率可以从23.8%提高到27.0%,H1单井净增油2.94×104m3。如果进一步增加动用的储层半径,仍用目前400m井距的现井网转驱(方案3和方案4),采收率仅有25.7%~26.3%;而通过加密井H1转驱(方案5),采收率可以从吞吐27.0%提高到35.9%,H1单井净增油8.13×104m3。
表3 LD27-2稠油油田不同开发方式开发效果对比
尽管上面结合附图对本发明的功能及工作过程进行了描述,但本发明并不局限于上述的具体功能和工作过程,上述的具体实施方式仅仅是示意性的,而不是限制性的,本领域的普通技术人员在本发明的启示下,在不脱离本发明宗旨和权利要求所保护的范围情况下,还可以做出很多形式,这些均属于本发明的保护之内。
Claims (1)
1.一种定量确定蒸汽吞吐转蒸汽驱地层温度时机的方法,其特征在于,包括以下步骤:
第一步,收集目标油藏基础数据资料,包括:原始油藏温度、地层条件下原油粘度、地面脱气原油的粘温曲线;
第二步,开展稠油样品流变性实验
提取目标油藏的稠油样品,脱气脱水后测试至少7个温度条件下的原油流变性,每个温度下测试剪切速率从1s-1到50s-1,间隔2s-1,共25个剪切速率下的原油粘度,同时根据公式(A-1)计算对应的屈服值τ;
τ=μov/1000 (A-1)
式中,τ表示屈服值,Pa;μo表示原油粘度,mPa·s;v表示剪切速率,1s-1;
第三步,回归不同温度下的流变性方程
在直角坐标系中,以剪切速率v为横坐标,屈服值τ为纵坐标,分别回归第二步中7个温度条件的流变性方程;特点是不同温度条件下,流变性方程的形式是相同的,均为τ=Av+B的线性形式,只是不同温度条件下的斜率A值和常数项B值不同;
第四步,确定稠油样品流变性拐点温度
观察第三步得到的不同温度下的流变性方程,随着温度的增加,常数项B逐渐降低;当温度为T0时,常数项B约等于零,此时对应的温度T0即为稠油样品从非牛顿流体转变为牛顿流体的拐点温度T拐点;
第五步,确定拐点温度对应的原油粘度
在直角坐标上,以温度lgT为横坐标,原油粘度lgμo为纵坐标绘制散点图,回归得到(A-2-1)形式的粘温关系式,其中C和D为回归得到的系数,不同的稠油样品回归得到不同的C和D值;
lgμo=C+DlgT (A-2-1)
式中,μo为原油粘度,mPa·s;T为温度,℃;得到原油粘度计算公式:
μo=10CTD (A-2-2)
第六步,绘制多轮蒸汽吞吐后地层温度分布图
在直角坐标系中,以距离注蒸汽井的距离为横坐标,以网格平均温度为纵坐标,绘制蒸汽吞吐N年后的网格平均温度随注蒸汽井的距离变化曲线,根据曲线确定注采井中间位置地层温度;如果注采井中间位置地层温度高于或等于拐点温度T中心≥T拐点,不考虑启动压力梯度的影响,认为达到转驱温度条件转蒸汽驱;如果注采井中间位置地层温度低于拐点温度T中心<T拐点,考虑启动压力梯度对能否建立有效注采***的影响;
第七步,确定启动压力梯度与流度关系
开展启动压力梯度实验,回归启动压力梯度与流度关系式。启动压力梯度与流体的性质、介质的表面作用和孔隙结构密切相关,不同油藏的启动压力梯度具有统一的公式样式(A-3-1):
G=E(K/μo)F (A-3-1)
式中,G表示启动压力梯度,MPa/m;K表示地层渗透率,10-3μm2,μo表示地层条件下的原油粘度,mPa·s;E和F为回归得到的系数,不同的油藏可回归得到不同的值;
第八步,确定注采井间的驱替压力梯度
根据最大注采压差ΔP,注采井间的井距L,确定注采井间中心位置的驱替压力梯度G驱替=ΔP/L;根据第七步中公式(A-3-1)计算启动压力梯度G启动;比较驱替压力梯度G驱替和启动压力梯度G启动的大小;
第九步,绘制不同注采井距对应的转驱温度图版
联立方程(A-2-2)和(A-3-1),得到:
对于等式右边,将地层渗透率K和第五步、第七步回归得到的系数C、D、E、F值代入,右边只剩下一个变量温度T;对于等式左边,在不同最大注采压差ΔP的情况下,分别代入不同注采井间的井距L,通过(A-4-1)计算对应的温度T;在直角坐标系中,以注采井间的井距L为横坐标,对应的温度T为纵坐标,绘制不同注采压差下不同注采井距的井间转驱温度图版;
第十步,判断是否达到转驱的温度时机
对于目标油田,根据最大注采压差ΔP,注采井间的井距L,查阅第九步的图版,确定转驱注采井中心位置所需的温度T转驱,与注采井中心位置地层温度T地层进行对比,当T地层≥T转驱时,可转蒸汽驱,当T地层<T转驱时,则不能转蒸汽驱。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010894337.7A CN112036033B (zh) | 2020-08-31 | 2020-08-31 | 一种定量确定蒸汽吞吐转蒸汽驱地层温度时机的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010894337.7A CN112036033B (zh) | 2020-08-31 | 2020-08-31 | 一种定量确定蒸汽吞吐转蒸汽驱地层温度时机的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112036033A true CN112036033A (zh) | 2020-12-04 |
CN112036033B CN112036033B (zh) | 2022-09-09 |
Family
ID=73586317
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010894337.7A Active CN112036033B (zh) | 2020-08-31 | 2020-08-31 | 一种定量确定蒸汽吞吐转蒸汽驱地层温度时机的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112036033B (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112412412A (zh) * | 2020-12-10 | 2021-02-26 | 中国石油天然气股份有限公司 | 超稠油蒸汽驱转驱启动条件的确定方法及装置 |
CN112796745A (zh) * | 2021-02-18 | 2021-05-14 | 中海油田服务股份有限公司 | 一种蒸汽驱油藏生产动态数据的确定方法和装置 |
CN114018782A (zh) * | 2021-10-18 | 2022-02-08 | 中国海洋石油集团有限公司 | 确定稠油水驱后转热水驱拟启动压力梯度动态变化的方法 |
CN114136838A (zh) * | 2021-11-19 | 2022-03-04 | 中国海洋石油集团有限公司 | 确定海上稠油不同含水阶段注热水驱原油粘度界限的方法 |
CN115324545A (zh) * | 2022-08-22 | 2022-11-11 | 中国石油大学(北京) | 变压式蒸汽辅助重力泄油的稠油开采方法 |
CN115422859A (zh) * | 2022-11-07 | 2022-12-02 | 西南石油大学 | 一种定量评价厚层稠油注蒸汽吞吐纵向波及系数的方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101852074A (zh) * | 2010-05-27 | 2010-10-06 | 中国石油天然气股份有限公司 | 一种用于层状超稠油油藏的开采方法及*** |
CN106593397A (zh) * | 2015-10-14 | 2017-04-26 | 中国石油天然气股份有限公司 | 稠油油藏的开采方法 |
CN107664031A (zh) * | 2016-07-29 | 2018-02-06 | 中国石油化工股份有限公司 | 通过确定水平井蒸汽驱井网形式提高采收率的方法 |
CN107664030A (zh) * | 2016-07-29 | 2018-02-06 | 中国石油化工股份有限公司 | 超稠油油藏水平井蒸汽驱提高采收率的方法 |
-
2020
- 2020-08-31 CN CN202010894337.7A patent/CN112036033B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101852074A (zh) * | 2010-05-27 | 2010-10-06 | 中国石油天然气股份有限公司 | 一种用于层状超稠油油藏的开采方法及*** |
CN106593397A (zh) * | 2015-10-14 | 2017-04-26 | 中国石油天然气股份有限公司 | 稠油油藏的开采方法 |
CN107664031A (zh) * | 2016-07-29 | 2018-02-06 | 中国石油化工股份有限公司 | 通过确定水平井蒸汽驱井网形式提高采收率的方法 |
CN107664030A (zh) * | 2016-07-29 | 2018-02-06 | 中国石油化工股份有限公司 | 超稠油油藏水平井蒸汽驱提高采收率的方法 |
Non-Patent Citations (2)
Title |
---|
宋传真等: "低渗稠油油藏蒸汽吞吐合理开发井距的确定――以叙利亚O油田为例", 《断块油气田》 * |
李延杰等: "海上M稠油油田吞吐后续转驱开发方案研究", 《西南石油大学学报(自然科学版)》 * |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112412412A (zh) * | 2020-12-10 | 2021-02-26 | 中国石油天然气股份有限公司 | 超稠油蒸汽驱转驱启动条件的确定方法及装置 |
CN112412412B (zh) * | 2020-12-10 | 2022-10-04 | 中国石油天然气股份有限公司 | 超稠油蒸汽驱转驱启动条件的确定方法及装置 |
CN112796745A (zh) * | 2021-02-18 | 2021-05-14 | 中海油田服务股份有限公司 | 一种蒸汽驱油藏生产动态数据的确定方法和装置 |
CN114018782A (zh) * | 2021-10-18 | 2022-02-08 | 中国海洋石油集团有限公司 | 确定稠油水驱后转热水驱拟启动压力梯度动态变化的方法 |
CN114136838A (zh) * | 2021-11-19 | 2022-03-04 | 中国海洋石油集团有限公司 | 确定海上稠油不同含水阶段注热水驱原油粘度界限的方法 |
CN114136838B (zh) * | 2021-11-19 | 2023-11-17 | 中国海洋石油集团有限公司 | 确定海上稠油不同含水阶段注热水驱原油粘度界限的方法 |
CN115324545A (zh) * | 2022-08-22 | 2022-11-11 | 中国石油大学(北京) | 变压式蒸汽辅助重力泄油的稠油开采方法 |
CN115324545B (zh) * | 2022-08-22 | 2023-10-03 | 中国石油大学(北京) | 变压式蒸汽辅助重力泄油的稠油开采方法 |
CN115422859A (zh) * | 2022-11-07 | 2022-12-02 | 西南石油大学 | 一种定量评价厚层稠油注蒸汽吞吐纵向波及系数的方法 |
CN115422859B (zh) * | 2022-11-07 | 2023-01-24 | 西南石油大学 | 一种定量评价厚层稠油注蒸汽吞吐纵向波及系数的方法 |
Also Published As
Publication number | Publication date |
---|---|
CN112036033B (zh) | 2022-09-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112036033B (zh) | 一种定量确定蒸汽吞吐转蒸汽驱地层温度时机的方法 | |
CN106651610B (zh) | 一种浅层超低渗砂岩油藏注水开发动态分析方法 | |
CN108952660B (zh) | 一种模拟注水井水压驱动裂缝延伸动态的方法 | |
CN105243182B (zh) | 致密油压裂水平井的动态储量计算方法 | |
CN104989341B (zh) | 一种确定低渗透油藏有效驱替注采井距的方法 | |
CN109614736B (zh) | 一种海上稠油油田蒸汽吞吐开发定向井与直井产能倍数的确定方法 | |
CN108386171B (zh) | 深水浊积砂岩注水开发油藏油井见水后注水强度优化方法 | |
CN106150463B (zh) | 一种砾岩油藏聚合物驱注入压力升幅确定方法 | |
CN111520136B (zh) | 考虑注水启动压力梯度的堵塞器嘴后压力计算方法 | |
CN109667564A (zh) | 一种海上稠油油田蒸汽吞吐开发定向井产能的确定方法 | |
Butler et al. | Steam and gas push (SAGP)-4; recent theoretical developments and laboratory results using layered models | |
RU2675115C1 (ru) | Способ разработки высоковязкой нефти | |
CN113356814B (zh) | 一种高相变油水原位乳化液提高稠油采收率的方法 | |
CN112597644B (zh) | 一种考虑层间干扰的多层合采油藏开发全周期产能校正方法 | |
CN107355200B (zh) | 一种纳微米颗粒分散体系改善水驱选井方法 | |
CN113673096A (zh) | 一种解堵增注剂处理半径的计算方法 | |
CN112031719A (zh) | 一种基于流动系数下启动压力的油藏开发方式优选方法 | |
CN115045642A (zh) | 一种水井油井推拉联动的稠油开采方法 | |
CN114018782A (zh) | 确定稠油水驱后转热水驱拟启动压力梯度动态变化的方法 | |
Yin | Present situation and development trend of oil production technology in unstable waterflooding | |
CN115746808A (zh) | 调堵剂及其应用 | |
RU2779868C1 (ru) | Способ разработки залежи высоковязкой или битумной нефти с использованием парных горизонтальных скважин | |
RU2776549C1 (ru) | Способ разработки месторождения высоковязкой нефти или битума с контролем развития паровой камеры в наблюдательных скважинах | |
Zhu | Application Precision Control Fracturing in Water Flooding Development During Extra High Water Cut | |
CN111764879B (zh) | 一种二元复合驱产剂窜流的评价方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |