CN109113688B - 一种确定非纯co2近混相驱最低混相压力的方法 - Google Patents
一种确定非纯co2近混相驱最低混相压力的方法 Download PDFInfo
- Publication number
- CN109113688B CN109113688B CN201811066498.6A CN201811066498A CN109113688B CN 109113688 B CN109113688 B CN 109113688B CN 201811066498 A CN201811066498 A CN 201811066498A CN 109113688 B CN109113688 B CN 109113688B
- Authority
- CN
- China
- Prior art keywords
- oil
- pure
- mmp
- content
- injected gas
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 25
- 230000015572 biosynthetic process Effects 0.000 claims description 13
- 238000006073 displacement reaction Methods 0.000 claims description 7
- 239000003921 oil Substances 0.000 abstract description 31
- 238000002474 experimental method Methods 0.000 abstract description 11
- 210000005239 tubule Anatomy 0.000 abstract description 10
- 238000004364 calculation method Methods 0.000 abstract description 7
- 239000010779 crude oil Substances 0.000 abstract description 7
- 238000002347 injection Methods 0.000 abstract description 6
- 239000007924 injection Substances 0.000 abstract description 6
- 239000012535 impurity Substances 0.000 abstract description 4
- 230000001419 dependent effect Effects 0.000 abstract description 3
- 238000011084 recovery Methods 0.000 description 5
- 125000004432 carbon atom Chemical group C* 0.000 description 4
- 150000001875 compounds Chemical class 0.000 description 4
- 230000000694 effects Effects 0.000 description 2
- 239000000203 mixture Substances 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 102000002274 Matrix Metalloproteinases Human genes 0.000 description 1
- 108010000684 Matrix Metalloproteinases Proteins 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 230000009977 dual effect Effects 0.000 description 1
- 238000001704 evaporation Methods 0.000 description 1
- 239000012530 fluid Substances 0.000 description 1
- 238000012417 linear regression Methods 0.000 description 1
- 239000011148 porous material Substances 0.000 description 1
- 238000012216 screening Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 239000007858 starting material Substances 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
Classifications
-
- 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/164—Injecting CO2 or carbonated water
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
-
- 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)
- Theoretical Computer Science (AREA)
- Mining & Mineral Resources (AREA)
- Geology (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- General Physics & Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- Chemical Kinetics & Catalysis (AREA)
- Geometry (AREA)
- Environmental & Geological Engineering (AREA)
- Fluid Mechanics (AREA)
- Evolutionary Computation (AREA)
- Computer Hardware Design (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geochemistry & Mineralogy (AREA)
- Carbon And Carbon Compounds (AREA)
- Organic Low-Molecular-Weight Compounds And Preparation Thereof (AREA)
Abstract
本发明公开了一种确定非纯CO2近混相驱最低混相压力的方法。该方法确定的非纯CO2的经验公式拟合过程中对杂质气体均应予以考虑。此外,考虑到注入气中的C1含量对MMP影响与原油组分中已有C1组分含量密切相关,因此在拟合过程中,采用非线性拟合。同时,考虑到注入气对MMP的影响程度与温度密切相关,因此将温度也作为基础参变量。拟合过程中,将非纯CO2条件下MMP与纯CO2条件下的MMP差值,作为因变量,将注入气杂质气体的各组分摩尔含量,油组分中C1含量,油层温度作为参变量,进行多元非线性拟合,得到注入气为非纯CO2时最小混相压力。本发明充分考虑了原油和注入气情况,计算结果更接近细管实验结果。
Description
技术领域
本发明涉及一种确定非纯CO2近混相驱最低混相压力的方法,属于气驱提高采收率技术领域。
背景技术
气驱尤其是CO2驱已成为重要的提高油田采收率的技术手段,在国内外获得广泛应用。最小混相压力(MMP)作为CO2驱应用研究的重要基础指标,越来越受到重视。由于油藏和流体条件的复杂性以及气源纯度的影响,国内东部大部分油田难以实现混相驱,CO2提高采收率效果并不明显。1986年Zick等在实验和理论计算的基础上,首次提出了近混相驱的概念,认为在油气两相界面张力(IFT)值较低时采收率能够达到 95%以上,近混相驱替过程为凝析-蒸发气驱双重作用的结果,提出了“油气界面张力为 0并非必要”的观点,引起广泛关注。从最终采收率的角度看,黏度、界面张力和孔隙大小分布的关系比是否完全混相更重要,在适当降低注入压力的条件下,近混相驱能够接近或达到混相驱的实施效果。在非纯CO2驱MMP研究方面,考虑到细管实验工作量大和细管模拟方法的复杂性,国内外学者一般采用经验公式的方法对目标油藏实施注气可行性进行初步筛选,但目前提出的经验公式与细管实验结果偏差较大,因此,研究一种确定非纯CO2近混相驱压力区间的方法很有必要。
发明内容
本发明的目的是提供一种确定非纯CO2近混相驱最低混相压力的方法。
本发明提供的确定非纯CO2驱油时最小混相压力的方法,包括:
1)测定油藏各组份含量和油层温度,其中,油藏各组分含量的单位为mol%;油层温度记为TR,单位为℃;
2)将步骤1)所得油藏各组分含量和TR代入如下公式I中,得到非纯CO2驱油时最小混相压力,记为MMPIMP;
所述公式I中,MMPIMP是注入气为非纯CO2时最小混相压力,单位为MPa;
XinjC1是注入气中C1摩尔百分含量,mol%;
XinjN2是注入气中N2摩尔百分含量,mol%;
XinjH2S是注入气中H2S摩尔百分含量,mol%;
Xinjc2-c6是注入气中C2-C6摩尔百分含量,mol%;(所述C2-C6是指注入气中含有的碳原子总数为2-6的所有化合物)
vol是地层油中C1和N2摩尔百分含量,mol%;
TR是油层温度,单位为℃;
MMPPURE是相同驱油条件下注入气为纯CO2的最小混相压力,单位为MPa。
上述方法中,所述MMPPURE是相同条件下注入气为纯CO2的最小混相压力;当注入气为纯CO2时,以油层温度、原油组分中(C2-C4)的摩尔含量、(C5-C6)的摩尔含量、以及C7+的平均摩尔质量作为自变量,以MMP作为因变量,采用regress函数利用matlab进行线性回归,即可按照如下公式II计算得到MMPPURE:
所述公式II中,TR是步骤1)所得油层温度,单位为℃;
MWC7+是地层油中C7+摩尔质量,g/mol;(所述C7+是指地层油中含有的碳原子总数为不小于7的所有化合物)
XC1是地层油中C1+N2摩尔百分含量,mol%;
XC2-C4是地层油中C2-C4摩尔百分含量;(所述C2-C4是指地层油中含有的碳原子总数为2-4的所有化合物),mol%
XC5-C6是地层油中C5-C6摩尔百分含量(所述C5-C6是指地层油中含有的碳原子总数为5-6的所有化合物),mol%。
而当注入气为非纯CO2时,由于注入气所含杂质气体一般为N2、H2S、C1、C2- C6,因此,非纯CO2的经验公式拟合过程中对杂质气体均应予以考虑。此外,考虑到注入气中的C1含量对MMP影响与原油组分中已有C1组分含量密切相关,因此在拟合过程中,采用非线性拟合。同时,考虑到注入气对MMP的影响程度与温度密切相关,因此将温度也作为基础参变量。
拟合过程中,将非纯CO2条件下MMP与纯CO2条件下的MMP差值,作为因变量,将注入气杂质气体的各组分摩尔含量,油组分中C1含量,油层温度作为参变量,进行多元非线性拟合,最终得到所述公式I。
本发明的计算方法较已有经验公式具有如下优点:
本发明方法充分考虑了非纯CO2气体中杂质含量,温度等关键因素,利用该经验公式计算的MMP结果与细管实验结果较为接近,误差小。
具体实施方式
下面结合具体实施例对本发明作进一步阐述,但本发明并不限于以下实施例。所述方法如无特别说明均为常规方法。所述原材料如无特别说明均能从公开商业途径获得。
实施例1、渤海Q油田非纯CO2(含量80%和55%)驱油时最小混相压力MMP 计算
细管实验用注入气和原油组成如表2和表3。
细管实验采用CFS-100多功能综合驱替系统,细管参数如表4。
表2、注入气组成
表3、实验油样主要物性
表4、实验所用细管参数表
1、细管实验得到MMP:
注CO2纯度为80%的注入气与原油的MMP为32.96MPa,
注CO2纯度为55%的注入气与原油的MMP为37.96MPa。
2、本发明公式计算得到MMP如下表5:
表5、纯/非纯CO2驱经验公式计算的MMP与细管实验/数模结果对比
3、现有经验公式计算MMP结果
对于目标油藏来说,注入气中杂质气体的含量范围比较大,现有经验公式中,仅有Yuan(2004)公式可用于非纯CO2驱MMP的计算,表6结果表明:非纯条件下,使用该现有经验公式相对误差均在15%以上。
表6、现有经验公式计算的MMP与细管实验结果对比(非纯CO2)
由上可知,本发明方法由于充分考虑了非纯CO2气体中杂质含量,温度等关键因素,利用本发明提供的方法计算的MMP结果与细管实验结果较为接近,误差较现有经验公式的15%以上可减小到3%左右,误差明显减小。
Claims (1)
1.一种确定用非纯CO2驱油时最小混相压力的方法,包括:
1)测定油藏各组份含量和油层温度,其中,油藏各组分含量的单位为mol%;油层温度记为TR,单位为℃;
2)将步骤1)所得油藏各组分含量和TR代入如下公式I中,得到非纯CO2驱油时最小混相压力,记为MMPIMP;
所述公式I中,MMPIMP是注入气为非纯CO2时最小混相压力,单位为MPa;
XinjC1是注入气中C1摩尔百分含量;
XinjN2是注入气中N2摩尔百分含量;
XinjH2S是注入气中H2S摩尔百分含量;
Xinjc2-c6是注入气中C2-C6摩尔百分含量;
vol是地层油中C1和N2摩尔百分含量;
TR是油层温度,单位为℃;
MMPPURE是相同驱油条件下注入气为纯CO2的最小混相压力,单位为MPa;
所述MMPPURE按照如下公式II计算得到:
所述公式II中,TR是步骤1)所得油层温度,单位为℃;
MWC7+是地层油中C7+摩尔质量,g/mol;
XC1是地层油中C1+N2摩尔百分含量;
XC2-C4是地层油中C2-C4摩尔百分含量;
XC5-C6是地层油中C5-C6摩尔百分含量。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811066498.6A CN109113688B (zh) | 2018-09-10 | 2018-09-10 | 一种确定非纯co2近混相驱最低混相压力的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811066498.6A CN109113688B (zh) | 2018-09-10 | 2018-09-10 | 一种确定非纯co2近混相驱最低混相压力的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109113688A CN109113688A (zh) | 2019-01-01 |
CN109113688B true CN109113688B (zh) | 2021-03-02 |
Family
ID=64859258
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811066498.6A Active CN109113688B (zh) | 2018-09-10 | 2018-09-10 | 一种确定非纯co2近混相驱最低混相压力的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109113688B (zh) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110059435A (zh) * | 2019-04-27 | 2019-07-26 | 西南石油大学 | 一种非纯二氧化碳混相驱最小混相压力gwo-lssvm预测方法 |
CN112177576A (zh) * | 2019-06-12 | 2021-01-05 | 中国石油天然气股份有限公司 | 一种注高含硫化氢和二氧化碳气体的混相机理判别方法 |
CN110598340B (zh) * | 2019-09-18 | 2020-09-01 | 中国石油大学(北京) | 一种注气驱油实验流体的确定方法及装置 |
CN112647905B (zh) * | 2019-10-10 | 2023-12-22 | 中国石油化工股份有限公司 | 一种确定天然气驱注入气组分的方法及天然气驱油方法 |
CN111256900B (zh) * | 2020-02-20 | 2021-01-12 | 中国石油大学(北京) | 油气最小混相压力确定方法及装置 |
CN111861026B (zh) * | 2020-07-28 | 2022-06-07 | 中国石油化工股份有限公司 | 一种快速获取苏北陆相油田最小混相压力的方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4513821A (en) * | 1984-02-03 | 1985-04-30 | Mobil Oil Corporation | Lowering CO2 MMP and recovering oil using carbon dioxide |
CN107013192A (zh) * | 2017-06-05 | 2017-08-04 | 中海石油(中国)有限公司 | 一种气驱近混相驱压力区间的界定方法 |
CN108266165A (zh) * | 2018-01-15 | 2018-07-10 | 中国石油化工股份有限公司 | 低渗油藏co2驱最小混相压力计算方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7707013B2 (en) * | 2005-12-05 | 2010-04-27 | Shell Oil Company | Method for estimating minimum miscibility enrichment |
-
2018
- 2018-09-10 CN CN201811066498.6A patent/CN109113688B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4513821A (en) * | 1984-02-03 | 1985-04-30 | Mobil Oil Corporation | Lowering CO2 MMP and recovering oil using carbon dioxide |
CN107013192A (zh) * | 2017-06-05 | 2017-08-04 | 中海石油(中国)有限公司 | 一种气驱近混相驱压力区间的界定方法 |
CN108266165A (zh) * | 2018-01-15 | 2018-07-10 | 中国石油化工股份有限公司 | 低渗油藏co2驱最小混相压力计算方法 |
Non-Patent Citations (1)
Title |
---|
非纯CO2原油体系混相机理研究;孙蓉;《万方》;20081231;第35-39页 * |
Also Published As
Publication number | Publication date |
---|---|
CN109113688A (zh) | 2019-01-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109113688B (zh) | 一种确定非纯co2近混相驱最低混相压力的方法 | |
CN102329768B (zh) | 用于油藏残余油气化开采的菌群构建方法 | |
Zang et al. | CO2 removal from synthesized ternary gas mixtures used hydrate formation with sodium dodecyl sulfate (SDS) as additive | |
CN110699058A (zh) | 一种用于水驱稠油的渗透改性降粘驱油剂及其制备方法 | |
CN106590606B (zh) | 耐温抗盐低张力泡沫及其在泡沫封堵中的应用 | |
CN204830680U (zh) | 一种控制二氧化碳冻堵的天然气乙烷回收装置 | |
CN1837574A (zh) | 热敏水溶性凝胶堵水调剖封窜技术 | |
CN104449641A (zh) | 一种降低co2非混相驱最小混相压力的调节剂及其应用方法 | |
CN101538461B (zh) | 一种注水井用降压增注剂及制备方法 | |
CN108729893A (zh) | 一种提高稠油油藏采收率的泡沫复合冷采方法 | |
CN103468236A (zh) | 一种含有丁烷的压裂液及其制备方法 | |
CN106350035A (zh) | 长链叔胺作为水‑co2交替驱中水相原位增粘剂的应用 | |
CN106635060A (zh) | 一种具有近晶相—胆甾相转变的液晶组合物 | |
CN106590607A (zh) | 耐温抗盐低张力泡沫剂组合物及其制备方法 | |
CN107384349A (zh) | 一种低渗孔道型油藏深部双段塞调剖剂及其使用方法 | |
CA1092966A (en) | Method of oil recovery employing enriched gas drive with control of evolved gas | |
CN109054795A (zh) | 一种超低渗透油藏驱油用杂化含氟表面活性剂的制备方法 | |
CN105754565A (zh) | 一种有机微球与无机硅酸盐复配型稠油热采封窜剂及其制备方法 | |
CN115703862A (zh) | 一种稠化二氧化碳驱油高效增稠剂及其制备方法 | |
CN114753818B (zh) | 提高稠油吞吐井产量的方法 | |
CN110683968B (zh) | 一种稠油开采用泡沫剂及其制备方法 | |
CN107384351A (zh) | 一种低渗孔道型油藏双段塞复配调剖剂及其使用方法 | |
CN113882840A (zh) | 一种利用生纳材料进行稠油吞吐冷采的方法 | |
CN104211188A (zh) | 一种钡锶铁防垢剂及其制备工艺 | |
Wang et al. | Low-temperature oxidation catalytic technology of light crude oil |
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 |