CN107885893A - 描述二氧化碳驱储层非均质性的构建方法 - Google Patents

描述二氧化碳驱储层非均质性的构建方法 Download PDF

Info

Publication number
CN107885893A
CN107885893A CN201610867559.3A CN201610867559A CN107885893A CN 107885893 A CN107885893 A CN 107885893A CN 201610867559 A CN201610867559 A CN 201610867559A CN 107885893 A CN107885893 A CN 107885893A
Authority
CN
China
Prior art keywords
permeability
carbon dioxide
dioxide flooding
gas
grade
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
Application number
CN201610867559.3A
Other languages
English (en)
Other versions
CN107885893B (zh
Inventor
王杰
戴涛
于金彪
谭保国
赵莹莹
史华
宋勇
宋桂茹
张波
汪勇
陈燕虎
易红霞
苏海波
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China Petroleum and Chemical Corp
Exploration and Development Research Institute of Sinopec Henan Oilfield Branch Co
Exploration and Development Research Institute of Sinopec Shengli Oilfield Co
Original Assignee
China Petroleum and Chemical Corp
Exploration and Development Research Institute of Sinopec Henan Oilfield Branch Co
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by China Petroleum and Chemical Corp, Exploration and Development Research Institute of Sinopec Henan Oilfield Branch Co filed Critical China Petroleum and Chemical Corp
Priority to CN201610867559.3A priority Critical patent/CN107885893B/zh
Publication of CN107885893A publication Critical patent/CN107885893A/zh
Application granted granted Critical
Publication of CN107885893B publication Critical patent/CN107885893B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • EFIXED CONSTRUCTIONS
    • E21EARTH OR ROCK DRILLING; MINING
    • E21BEARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
    • E21B43/00Methods or apparatus for obtaining oil, gas, water, soluble or meltable materials or a slurry of minerals from wells
    • E21B43/16Enhanced recovery methods for obtaining hydrocarbons
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q50/00Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/02Agriculture; Fishing; Forestry; Mining
    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02ATECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A10/00TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE at coastal zones; at river basins
    • Y02A10/40Controlling or monitoring, e.g. of flood or hurricane; Forecasting, e.g. risk assessment or mapping

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Theoretical Computer Science (AREA)
  • Mining & Mineral Resources (AREA)
  • Business, Economics & Management (AREA)
  • General Physics & Mathematics (AREA)
  • Geology (AREA)
  • Economics (AREA)
  • Strategic Management (AREA)
  • Animal Husbandry (AREA)
  • Marine Sciences & Fisheries (AREA)
  • Geochemistry & Mineralogy (AREA)
  • Health & Medical Sciences (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • General Health & Medical Sciences (AREA)
  • Human Resources & Organizations (AREA)
  • Marketing (AREA)
  • Primary Health Care (AREA)
  • Agronomy & Crop Science (AREA)
  • Tourism & Hospitality (AREA)
  • General Business, Economics & Management (AREA)
  • Fluid Mechanics (AREA)
  • Environmental & Geological Engineering (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明提供一种描述二氧化碳驱储层非均质性的构建方法,包括:步骤1,计算区块的渗透率级差,并根据研究区的渗透率级差,确定合理的渗透率级差界限;步骤2,对某一研究区内的渗透率解释结果进行处理,计算网格模型的渗透率级差;步骤3,利用数值模拟手段研究不同方法对二氧化碳驱见气时间的影响程度;步骤4,分析二氧化碳驱中原油物理化学性质对气油比的影响;步骤5,利用数值模拟手段分析不同渗透率值模型的流量差;步骤6,利用实际模型的二氧化碳驱气油比拟合对比校正渗透率平均值。该构建方法根据数值模拟中所能提供的参数,建立了渗透率曲线离散化方法,为二氧化碳驱数值模拟精度的提高以及其驱替机理的定量表征提供了数值模拟依据。

Description

描述二氧化碳驱储层非均质性的构建方法
技术领域
本发明涉及油藏数值模拟理论研究、技术应用及油气田开发提高油藏采收率领域,特别是涉及到一种描述二氧化碳驱储层非均质性的构建方法。
背景技术
随着中石化绿色环保理念的深入、二氧化碳驱技术代表了三次采油的前沿方向,市场需求量大,经济效益显著,既可以提高原油采收率又可以降低二氧化碳排放。低渗透油藏存在注水开发困难、水驱效果差、采收率低。胜利油田小于5毫达西的特低渗油藏无法实现注水开发,二氧化碳驱具有很大的优势。二氧化碳驱的优势在于降低原油粘度、膨胀原油体积系数、蒸发原油中间烃组分、降低界面张力、溶解气驱、增加注入率。但是同时,二氧化碳相对于原油存在密度和粘度差异,导致高渗透层的早见气、采收程度差。当储层存在纵向非均质性时,高渗透层的层厚会加剧这种影响,数值模拟计算的储层见气时间比较短、采出程度较实际偏高。原因在于储层建模时,对渗透率曲线离散化时,降低了储层的非均质性,渗透率更加平均、二氧化碳驱推进更加均匀,和矿场实际不符合。曲线离散化的常规方法是算术平均、几何平均、调和平均、均方根平均等方法,无论哪种方法都是整个储层参与属性的平均,降低了气体推进的速度。目前国内外研究人员对于气驱推进速度的分析、对比和优化,缺乏有针对性的描述。为此我们发明了一种新的描述二氧化碳驱储层非均质性的构建方法,解决了以上技术问题。
发明内容
本发明的目的是提供一种利用数值模拟手段定义储层渗透率率级差对见气时间和采出程度的影响程度,对渗透率曲线离散化进行定量校正的描述二氧化碳驱储层非均质性的构建方法。
本发明的目的可通过如下技术措施来实现:描述二氧化碳驱储层非均质性的构建方法,该描述二氧化碳驱储层非均质性的构建方法包括:步骤1,计算区块的渗透率级差,并根据研究区的渗透率级差,确定合理的渗透率级差界限;步骤2,对某一研究区内的渗透率解释结果进行处理,计算网格模型的渗透率级差;步骤3,利用数值模拟手段研究不同方法对二氧化碳驱见气时间的影响程度;步骤4,分析二氧化碳驱中原油物理化学性质对气油比的影响;步骤5,利用数值模拟手段分析不同渗透率值模型的流量差;步骤6,利用实际模型的二氧化碳驱气油比拟合对比校正渗透率平均值。
本发明的目的还可通过如下技术措施来实现:
在步骤1中,对某一研究区目的层的单井测井解释结果进行区间分布特征概率分析,计算区块的渗透率级差。
在步骤1中,利用概念模型针对研究区的渗透率级差,研究储层的渗透率级差对见气时间和采出程度的影响因素分析,建立纵向渗透率级差与见气时间和采出程度的相关曲线,确定合理的渗透率级差界限。
在步骤2中,对某一研究区内的渗透率解释结果进行网格离散化处理,分析网格数据分布特征,计算网格模型的渗透率级差。
在步骤3中,对比不同的曲线离散化方法包括算术平均、最大值、最小值、均方根平均、调和平均这些方法得到的网格的渗透率值,利用数值模拟手段研究不同离散化方法对二氧化碳驱见气时间的影响程度。
在步骤5中,利用数值模拟手段分析不同渗透率值模型的流量差,求解平均含气饱和度,计算边界流量、计算网格间的压差,对于渗透率小于网格块内渗透率最大值/渗透率级差界限时,进行人为的截断,剩余的渗透率值用RMS方法进行渗透率离散化处理。
在步骤6中,进行实际区块模拟,利用现场实测的气油比曲线进行拟合差异分析,当误差大于5%时,重新步骤5中的流量差计算,计算新的渗透率属性值,通过反复验证修改,直至计算的气油比曲线与实测气油比曲线误差小于5%时,模拟得到渗透率场为合理的渗透率属性场。
本发明中的描述二氧化碳驱储层非均质性的构建方法,更清楚的描述这种非均质性对数值模拟的结果影响,通过储层的渗透率级差对二氧化碳驱见气的影响因素进行分析,对储层渗透率曲线的离散化时进行物性下限的截断,消除低渗透层的气体推进速度的降低,加快高渗透层的气体推进速度,更加符合矿场实际。该方法利用数值模拟手段定义储层渗透率率级差对见气时间和采出程度的影响程度,对渗透率曲线离散化进行定量校正,准确预测储层的见气时间和采出程度。为气驱矿场实践提供有效的数据支持,为油藏工程师提供准确的数据指导。该方法根据数值模拟中所能提供的参数,建立了渗透率曲线离散化方法,为二氧化碳驱数值模拟精度的提高以及其驱替机理的定量表征提供了数值模拟依据。
附图说明
图1为本发明的描述二氧化碳驱储层非均质性的构建方法的一具体实施例的流程图;
图2为本发明的一具体实施例中渗透率级差概念模型设计有无隔层的对比图;
图3为本发明的一具体实施例中层间渗透率级差与无因次见气时间曲线;
图4为本发明的一具体实施例中储层空间渗透率级差分布图;
图5为本发明的一具体实施例中不同储层非均质性离散化方法结果对比曲线;
图6为本发明的一具体实施例中渗透率曲线离散化处理方法示意图;
图7为本发明的一具体实施例中粗细模型离散化气油比曲线对比曲线;
图8为本发明的一具体实施例中精确拟合单井的气油比曲线图版。
具体实施方式
为使本发明的上述特征和其他目的、特征和优点能更明显易懂,下文特举出较佳实施例,并配合附图所示,作详细说明如下。
如图1所示,图1为本发明的描述二氧化碳驱储层非均质性的构建方法的流程图。
在步骤101,对某一研究区目的层的单井测井解释结果进行区间分布特征概率分析,计算区块的渗透率级差。
在步骤102,利用概念模型针对研究区的渗透率级差,研究储层的渗透率级差对见气时间和采出程度的影响因素分析,建立纵向渗透率级差与见气时间和采出程度的相关曲线,确定合理的渗透率级差界限;
在步骤103,对某一研究区内的渗透率解释结果进行网格离散化处理,分析网格数据分布特征,计算网格模型的渗透率级差;
在步骤104,以单网格为例,对比不同的曲线离散化方法(算术平均、最大值、最小值、RMS均方根、调和平均等)得到的网格的渗透率值,利用数值模拟手段研究不同离散化方法对二氧化碳驱见气时间的影响程度;
在步骤105,二氧化碳驱过程中,二氧化碳在不同压力下的粘度和密度存在差异,分析二氧化碳驱中原油物理化学性质对气油比的影响分析。
在步骤106,利用数值模拟手段分析不同渗透率值模型的流量差,求解平均含气饱和度,计算边界流量、计算网格间的压差,对于渗透率小于(网格块内渗透率最大值/渗透率级差界限)值进行人为的截断,剩余的渗透率值用RMS方法进行渗透率离散化处理。
在步骤107,进行模型的二氧化碳驱气油比拟合对比分析。利用实际模型的二氧化碳驱气油比拟合对比校正渗透率平均值。进行实际区块模拟,利用现场实测的气油比曲线进行拟合差异分析,当误差大于5%时,重新步骤5中的流量差计算,计算新的渗透率属性值,通过反复验证修改,直至计算的气油比曲线与实测气油比曲线误差小于5%时,模拟得到渗透率场为合理的渗透率属性场。
以不同的渗透率级别为例,利用数值模拟手段测定渗透率级差对气油比影响曲线。图2是设计有无隔层的概念模型对比图,研究隔夹层和高低渗透层的厚度对见气时间的影响程度进行分析,去除这两个因素对见气时间的影响。根据本发明中所提出的方法计算出层间非均质级差与无因次见气时间曲线(如图3所示),研究影响因素影响程度,可以判断渗透率级差对见气时间的影响程度,对高89地区的地质模型的网格渗透率级差分布如图4所示。不同渗透率曲线离散化方法(算术平均、最大值、最小值、RMS均方根、调和平均等)计算的见气时间不同,无法通过单一的方法计算符合实际的气油比曲线(如图5所示),通过对渗透率曲线进行特殊的处理,利用数值模拟手段分析不同渗透率值模型的流量差,求解平均含气饱和度,计算边界流量、计算网格间的压差,对于渗透率小于(网格块内渗透率最大值/渗透率级差界限)值进行人为的截断,剩余的渗透率值应用RMS进行渗透率离散化处理(如图6所示),得到合适的相渗曲线级差以及网格渗透率值。图7是采用该方法确定的合适的概念模型的气油比拟合曲线。根据这些结果可以得到实际区块的单井的气油比曲线拟合图版如图8。

Claims (7)

1.描述二氧化碳驱储层非均质性的构建方法,其特征在于,该描述二氧化碳驱储层非均质性的构建方法包括:
步骤1,计算区块的渗透率级差,并根据研究区的渗透率级差,确定合理的渗透率级差界限;
步骤2,对某一研究区内的渗透率解释结果进行处理,计算网格模型的渗透率级差;
步骤3,利用数值模拟手段研究不同方法对二氧化碳驱见气时间的影响程度;
步骤4,分析二氧化碳驱中原油物理化学性质对气油比的影响;
步骤5,利用数值模拟手段分析不同渗透率值模型的流量差;
步骤6,利用实际模型的二氧化碳驱气油比拟合对比校正渗透率平均值。
2.根据权利要求1所述的描述二氧化碳驱储层非均质性的构建方法,其特征在于,在步骤1中,对某一研究区目的层的单井测井解释结果进行区间分布特征概率分析,计算区块的渗透率级差。
3.根据权利要求1所述的描述二氧化碳驱储层非均质性的构建方法,其特征在于,在步骤1中,利用概念模型针对研究区的渗透率级差,研究储层的渗透率级差对见气时间和采出程度的影响因素分析,建立纵向渗透率级差与见气时间和采出程度的相关曲线,确定合理的渗透率级差界限。
4.根据权利要求1所述的描述二氧化碳驱储层非均质性的构建方法,其特征在于,在步骤2中,对某一研究区内的渗透率解释结果进行网格离散化处理,分析网格数据分布特征,计算网格模型的渗透率级差。
5.根据权利要求1所述的描述二氧化碳驱储层非均质性的构建方法,其特征在于,在步骤3中,对比不同的曲线离散化方法包括算术平均、最大值、最小值、均方根平均、调和平均这些方法得到的网格的渗透率值,利用数值模拟手段研究不同离散化方法对二氧化碳驱见气时间的影响程度。
6.根据权利要求1所述的描述二氧化碳驱储层非均质性的构建方法,其特征在于,在步骤5中,利用数值模拟手段分析不同渗透率值模型的流量差,求解平均含气饱和度,计算边界流量、计算网格间的压差,对于渗透率小于网格块内渗透率最大值/渗透率级差界限时,进行人为的截断,剩余的渗透率值用RMS方法进行渗透率离散化处理。
7.根据权利要求1所述的描述二氧化碳驱储层非均质性的构建方法,其特征在于,在步骤6中,进行实际区块模拟,利用现场实测的气油比曲线进行拟合差异分析,当误差大于5%时,重新步骤5中的流量差计算,计算新的渗透率属性值,通过反复验证修改,直至计算的气油比曲线与实测气油比曲线误差小于5%时,模拟得到渗透率场为合理的渗透率属性场。
CN201610867559.3A 2016-09-29 2016-09-29 描述二氧化碳驱储层非均质性的构建方法 Active CN107885893B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610867559.3A CN107885893B (zh) 2016-09-29 2016-09-29 描述二氧化碳驱储层非均质性的构建方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610867559.3A CN107885893B (zh) 2016-09-29 2016-09-29 描述二氧化碳驱储层非均质性的构建方法

Publications (2)

Publication Number Publication Date
CN107885893A true CN107885893A (zh) 2018-04-06
CN107885893B CN107885893B (zh) 2021-06-01

Family

ID=61769836

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610867559.3A Active CN107885893B (zh) 2016-09-29 2016-09-29 描述二氧化碳驱储层非均质性的构建方法

Country Status (1)

Country Link
CN (1) CN107885893B (zh)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109635410A (zh) * 2018-12-06 2019-04-16 中国华能集团清洁能源技术研究院有限公司 一种基于孔隙系统的渗透率离散元模拟方法
CN109884726A (zh) * 2019-03-07 2019-06-14 中国石油大学(北京) 气驱油藏的见气时间预测方法和装置
CN111810142A (zh) * 2020-08-06 2020-10-23 西南石油大学 一种通过平板模型模拟井组开发过程中防气窜的方法
CN112800640A (zh) * 2019-11-14 2021-05-14 北京中科联华石油科学研究院 一种确定地质模型垂向网格精度的方法
CN113389531A (zh) * 2020-03-12 2021-09-14 中国石油化工股份有限公司 Co2注采方法、电子设备及介质

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101876241A (zh) * 2009-04-30 2010-11-03 中国石油天然气股份有限公司 一种提高正韵律厚油层水驱采收率的方法
US20120109603A1 (en) * 2009-06-22 2012-05-03 Ning Li Quantitative calculation method for oil (gas) saturation of fractured reservoir
CN102507412A (zh) * 2011-11-07 2012-06-20 中国石油集团川庆钻探工程有限公司 碳酸盐岩油藏等效模型的裂缝-基质渗透率级差判别方法
CN104453804A (zh) * 2014-10-17 2015-03-25 中国海洋石油总公司 一种气驱油藏开发动态监测与评价方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101876241A (zh) * 2009-04-30 2010-11-03 中国石油天然气股份有限公司 一种提高正韵律厚油层水驱采收率的方法
US20120109603A1 (en) * 2009-06-22 2012-05-03 Ning Li Quantitative calculation method for oil (gas) saturation of fractured reservoir
CN102507412A (zh) * 2011-11-07 2012-06-20 中国石油集团川庆钻探工程有限公司 碳酸盐岩油藏等效模型的裂缝-基质渗透率级差判别方法
CN104453804A (zh) * 2014-10-17 2015-03-25 中国海洋石油总公司 一种气驱油藏开发动态监测与评价方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
宋道万: "二氧化碳混相驱数值模拟结果的主要影响因素", 《油气地质与采收率》 *
王杰 等: "一种通过数值模拟手段划分CO2驱替相带的新方法——以高89块油藏为例", 《科技导报》 *

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109635410A (zh) * 2018-12-06 2019-04-16 中国华能集团清洁能源技术研究院有限公司 一种基于孔隙系统的渗透率离散元模拟方法
CN109635410B (zh) * 2018-12-06 2021-08-17 中国华能集团清洁能源技术研究院有限公司 一种基于孔隙系统的渗透率离散元模拟方法
CN109884726A (zh) * 2019-03-07 2019-06-14 中国石油大学(北京) 气驱油藏的见气时间预测方法和装置
CN109884726B (zh) * 2019-03-07 2020-08-04 中国石油大学(北京) 气驱油藏的见气时间预测方法和装置
CN112800640A (zh) * 2019-11-14 2021-05-14 北京中科联华石油科学研究院 一种确定地质模型垂向网格精度的方法
CN113389531A (zh) * 2020-03-12 2021-09-14 中国石油化工股份有限公司 Co2注采方法、电子设备及介质
CN111810142A (zh) * 2020-08-06 2020-10-23 西南石油大学 一种通过平板模型模拟井组开发过程中防气窜的方法

Also Published As

Publication number Publication date
CN107885893B (zh) 2021-06-01

Similar Documents

Publication Publication Date Title
CN104879103B (zh) 一种分层注水效果分析方法
CN107885893A (zh) 描述二氧化碳驱储层非均质性的构建方法
MacMinn et al. CO2 migration in saline aquifers. Part 2. Capillary and solubility trapping
CN103114830B (zh) 一种富气驱水/气交替注入方法
CN109948272A (zh) 基于井间连通性的调堵动态预测方法和系统
CN104616350A (zh) 缝洞型碳酸盐油藏三维物理模型建立方法
CN105740563A (zh) 一种成熟油田二次开发之优势通道识别技术
CN104632157A (zh) 低渗透油藏均衡驱替方法
CN105696986A (zh) 一种新型复合驱驱油实验/试验模拟方法
CN107292074B (zh) 一种井间连通性判断方法
CN107448197B (zh) 一种海相砂岩油藏窜流通道空间发育定量识别的方法
CN103939066A (zh) 一种一注多采井组定注水量确定油井产液量的方法
CN104912537B (zh) 一种巨厚潜山裂缝油藏的井网部署方法
CN104091069A (zh) 确定非均质储层各层位和位置驱油效率和波及系数的方法
CN104060985A (zh) 一种层状油藏调剖堵水堵剂进入深度测试方法及系统
CN112147684B (zh) 表征同沉积断层活动强度的方法和装置
CN105089659B (zh) 一种砾岩油藏渗流单元识别方法
CN106481315B (zh) 陆上砂岩油藏单井可采储量快速确定模型及建立方法
CN104712301A (zh) 一种确定火山岩储层压裂规模的方法
CN107288595A (zh) 一种注水利用率的评价方法
Wei et al. A comprehensive performance evaluation methodology for miscible gas flooding: A case study in a giant carbonate reservoir in Middle East
CN111209687A (zh) 一种用于判识co2驱气窜的方法
CN113032996B (zh) 一种用于低渗裂缝型油藏水平井的水窜通道识别方法
CN110513089A (zh) 一种海上稠油热采吞吐开发产能倍数确定方法
CN102606120B (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