CN106802985B - 基于含水率的稠油边底水油藏多轮co2吞吐的评价方法 - Google Patents
基于含水率的稠油边底水油藏多轮co2吞吐的评价方法 Download PDFInfo
- Publication number
- CN106802985B CN106802985B CN201710003679.3A CN201710003679A CN106802985B CN 106802985 B CN106802985 B CN 106802985B CN 201710003679 A CN201710003679 A CN 201710003679A CN 106802985 B CN106802985 B CN 106802985B
- Authority
- CN
- China
- Prior art keywords
- factor
- membership
- formula
- calculated
- degree
- 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
- 238000011156 evaluation Methods 0.000 title claims abstract description 112
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 title claims abstract description 45
- 239000011159 matrix material Substances 0.000 claims abstract description 36
- 239000000243 solution Substances 0.000 claims description 38
- 238000004088 simulation Methods 0.000 claims description 21
- 239000013598 vector Substances 0.000 claims description 20
- 238000002347 injection Methods 0.000 claims description 18
- 239000007924 injection Substances 0.000 claims description 18
- 238000004519 manufacturing process Methods 0.000 claims description 16
- 239000012530 fluid Substances 0.000 claims description 15
- 239000012634 fragment Substances 0.000 claims description 12
- 238000011065 in-situ storage Methods 0.000 claims description 12
- 238000011160 research Methods 0.000 claims description 11
- 239000011229 interlayer Substances 0.000 claims description 9
- 230000015572 biosynthetic process Effects 0.000 claims description 6
- 238000013461 design Methods 0.000 claims description 6
- 239000000203 mixture Substances 0.000 claims description 6
- 230000033764 rhythmic process Effects 0.000 claims description 6
- 239000007788 liquid Substances 0.000 claims description 3
- 230000035699 permeability Effects 0.000 claims description 3
- 241000208340 Araliaceae Species 0.000 claims 1
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 claims 1
- 235000003140 Panax quinquefolius Nutrition 0.000 claims 1
- 235000008434 ginseng Nutrition 0.000 claims 1
- 238000004364 calculation method Methods 0.000 abstract description 8
- 230000010534 mechanism of action Effects 0.000 abstract description 5
- 238000005259 measurement Methods 0.000 abstract 1
- 239000003921 oil Substances 0.000 description 20
- 238000000034 method Methods 0.000 description 15
- 238000011161 development Methods 0.000 description 5
- 230000000694 effects Effects 0.000 description 5
- 238000011084 recovery Methods 0.000 description 3
- 239000010779 crude oil Substances 0.000 description 2
- HJUFTIJOISQSKQ-UHFFFAOYSA-N fenoxycarb Chemical compound C1=CC(OCCNC(=O)OCC)=CC=C1OC1=CC=CC=C1 HJUFTIJOISQSKQ-UHFFFAOYSA-N 0.000 description 2
- 238000011835 investigation Methods 0.000 description 2
- 230000003068 static effect Effects 0.000 description 2
- 239000003643 water by type Substances 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 235000013399 edible fruits Nutrition 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 230000004044 response Effects 0.000 description 1
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/164—Injecting CO2 or carbonated water
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION 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/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/02—Agriculture; Fishing; Forestry; Mining
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Mining & Mineral Resources (AREA)
- Theoretical Computer Science (AREA)
- Geology (AREA)
- Business, Economics & Management (AREA)
- General Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Animal Husbandry (AREA)
- Marine Sciences & Fisheries (AREA)
- General Engineering & Computer Science (AREA)
- Geometry (AREA)
- Evolutionary Computation (AREA)
- Economics (AREA)
- General Health & Medical Sciences (AREA)
- Human Resources & Organizations (AREA)
- Marketing (AREA)
- Agronomy & Crop Science (AREA)
- Strategic Management (AREA)
- Tourism & Hospitality (AREA)
- General Business, Economics & Management (AREA)
- Chemical & Material Sciences (AREA)
- Chemical Kinetics & Catalysis (AREA)
- Environmental & Geological Engineering (AREA)
- Fluid Mechanics (AREA)
- Computer Hardware Design (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geochemistry & Mineralogy (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明提供了一种基于含水率的稠油边底水油藏多轮CO2吞吐的评价方法。该评价方法为：测定实际油藏各因素的实际参数，以此为依据建立油藏模型，模拟得到吞吐后开井含水率；以吞吐后开井含水率为评价关联对象，计算各因素的修正权重并建立隶属度计算方法，得到模拟方案和实际方案的评价矩阵，再结合各因素的修正权重，计算各模拟方案和实际方案所有因素的综合得分的平均值和Cj并排序；对比和Cj的大小，给出评价结论；若则评价结果为较好；若则评价结果为较差。该评价方法能够综合反映CO2吞吐在控水和增油两个方面的作用机制，对指导稠油边底水油藏多轮CO2吞吐评价具有重要的意义。
Description
技术领域
本发明属于稠油CO2吞吐效果评价技术领域，涉及一种基于含水率的稠油边底水油藏多轮CO2吞吐的评价方法。
背景技术
稠油边底水油藏在开发过程中，原油粘度高，含水上升快，原油采出程度低；通过CO2吞吐，能够起到降粘、增油、控水的作用。
传统的开发效果评价主要以增油量为评价目标，但是对于稠油边底水油藏，CO2吞吐后的控水作用表现得尤为显著：一方面，井筒附近的含水率明显降低；另一方面，含水率下降提高了油水流度比，促进了原油的采出。
因此，建立一种基于吞吐后开井含水率的效果评价方法，能够综合反映CO2吞吐在控水和增油两个方面的作用机制，对指导稠油边底水油藏多轮CO2吞吐评价具有重要的现实意义。
常规评价方法通过模糊综合评价，建立各影响因素与评价目标的隶属度关系。然而，该方法在实际应用过程中，常常受到现场资料不全，或者评价井数少等原因的制约，使得权重以及评价目标与各因素隶属度关系的建立，由于缺乏数据支持，受到过多人为因素的干扰，从而导致评价方法不科学，评价结果不合理。
发明内容
针对现场资料不全带来的权重难赋值、评价矩阵难计算、基于吞吐后开井含水率效果的评价结果不合理等问题，本发明的目的在于提供一种基于吞吐后开井含水率的稠油边底水多轮CO2吞吐的评价方法，能够综合反映CO2吞吐在控水和增油两个方面的作用机制，对指导稠油边底水油藏多轮CO2吞吐评价具有现实意义。
本发明的目的通过以下技术方案得以实现：
本发明提供一种基于吞吐后开井含水率的稠油边底水多轮CO2吞吐的评价方法，其包括以下步骤：
步骤一，测定实际油藏各因素的实际参数，以此为依据建立油藏模型，模拟得到吞吐后开井含水率；
步骤二，以吞吐后开井含水率为评价关联对象，计算各因素的修正权重；
步骤三，借助油藏数值模拟技术建立各因素隶属度计算方法，分别得到模拟方案和实际方案的评价矩阵；
步骤四，根据模拟方案评价矩阵和各因素的修正权重，计算各模拟方案所有因素的综合得分的平均值
步骤五，根据实际方案评价矩阵和各因素的修正权重，计算实际方案所有因素的综合得分Cj并排序；
步骤六，对比和Cj的大小，给出评价结论；若则评价结果为较好；若则评价结果为较差。
上述评级方法中，所述吞吐后开井含水率指吞吐措施结束后，刚开井时测得的单井含水率。
上述评价方法中，所述隶属度属于模糊评价函数里的概念，模糊综合评价是对受多种因素影响的事物做出全面评价的一种十分有效的多因素决策方法，其特点是评价结果不是绝对地肯定或否定，而是以一个模糊集合来表示。若对论域(研究的范围) U中的任一元素x，都有一个数A(x)∈0，1与之对应，则称A为U上的模糊集， A(x)称为x对A的隶属度。当x在U中变动时，A(x)就是一个函数，称为A 的隶属函数。隶属度A(x)越接近于1，表示x属于A的程度越高，A(x)越接近于0表示x属于A的程度越低。用取值于区间0，1的隶属函数A(x)表征x属于A 的程度高低。
上述的评价方法中，优选地，各因素的修正权重是采用式(I)计算得到的：
其中，ai为第i项因素的经验权重；Ai为第i项因素的修正权重；i＝1，2，……， n(n≤N)；N为方案评价所需的因素总数；n为现场资料实际提供的因素总数。
上述的评价方法中，优选地，借助油藏数值模拟技术建立各因素隶属度计算方法，分别得到模拟方案和实际方案的评价矩阵的步骤包括：
根据测得的实际油藏各因素的实际参数，对各因素设计合理的参数变化范围(即根据油田实际情况中各因素参数分布范围，并适当扩大来设计，扩大范围一般为± 20％)；
借助油藏数值模拟技术，预测各因素参数变化时对应的吞吐后开井含水率；所述吞吐后开井含水率指吞吐措施结束后，刚开井时测得的单井含水率；
依据吞吐后开井含水率的大小，建立各因素与隶属度的函数关系式；
依据各因素对应的隶属度，得到模拟方案评价矩阵；
依据建立的隶属度计算方法，计算实际方案各因素取值对应的隶属度，得到实际方案评价矩阵。
上述评价方法中，所述油藏数值模拟技术是油气田开发生产预测研究中常用的一直方法，它是借助ECLIPSE或者CMG等数值模拟专业软件，以实际油藏特征参数为依据，建立油藏模型，模拟油田开采，从而预测油藏采收率及开发动态数据的一种方法。
上述的评价方法中，优选地，所述各因素包括隔夹层位置、沉积韵律、出水段长度、出水段位置、水平段位置、平行井距离、吞吐轮次、地层倾角、油层有效厚度、地层非均质性、投产含水、出水段非均质程度、吞吐时机、周期注气量、注气速度、焖井时间、开井后采液速度和平行井产液速度中的多种。
上述的评价方法，优选地，所述各因素与隶属度的函数关系式分别如下：
利用式(1)得到隔夹层位置对应的隶属度y1a、y1b和y1c：
(注：y1a＝1(无夹层)表示当无夹层时，隔夹层位置对应的隶属度为1；其他的关系式的含义表示与此类似。)
利用式(2)得到沉积韵律对应的隶属度y2a、y2b和y2c：
利用式(3)得到出水段长度对应的隶属度y3a、y3b和y3c：
利用式(4)得到出水段位置对应的隶属度y4a、y4b和y4c：
其中，A点为水平井的趾端、B点为水平井的跟端；
利用式(5)得到水平段位置对应的隶属度y5a、y5b和y5c：
利用式(6)得到平行井距离对应的隶属度y6a、y6b和y6c：
其中：“相距”是指相邻两个平行水平井之间的垂向距离；
利用式(7)得到吞吐轮次对应的隶属度y7：
利用式(8)计算地层倾角对应的隶属度y8a、y8b和y8c：
其中，θ为地层倾角，°；
利用式(9)计算油层有效厚度对应隶属度y9a、y9b和y9c：
其中，h为油层有效厚度，m；
利用式(10)计算地层非均质性对应的隶属度y10a、y10b和y10c：
其中，Vk为洛伦兹系数，用于表示地层非均质性；
利用式(11)计算投产含水对应的隶属度y11a、y11b和y11c：
其中，fw0为投产含水；
利用式(12)出水段非均质程度对应的隶属度y12a、y12b和y12c：
其中，Kv为渗透率极差，用于表示出水段非均质程度；
利用式(13)计算吞吐时机对应的隶属度y13a、y13b和y13c：
其中，f为吞吐时机；
利用式(14)计算周期注气量对应的隶属度y14a、y14b和y14c：
其中，Qig为周期注气量，sm3；
利用式(15)计算注气速度对应的隶属度y15a、y15b和y15c：
其中，vig为注气速度，sm3/d；
利用式(16)计算焖井时间对应的隶属度y16a、y16b和y16c：
其中，t为焖井时间，d；
利用式(17)计算开井后采液速度对应的隶属度y17a、y17b和y17c：
其中，vl为开井后采液速度，rm3/d；
利用式(18)计算平行井采液速度对应的隶属度y18a、y18b和y18c：
其中，ql为平行井采液速度，rm3/d。
上述的评价方法中，优选地，根据模拟方案评价矩阵和各因素的修正权重，计算各模拟方案所有因素的综合得分的平均值(作为评价界限)的步骤包括：
首先将计算得到的各因素的修正权重组成修正权重向量，计算修正权重向量与模拟方案评价矩阵的乘积，得到模拟方案的得分向量；利用式(II)计算模拟方案所有因素的综合得分cj：
其中，Ai为第i项因素的修正权重；i＝1，2，......，n(n≤N)；N为方案评价所需因素总数；n为现场资料实际提供的因素总数；bji为第j项模拟方案的第i项因素对应的隶属度；j＝1，2，......，m；m为模拟方案总数；cj为第j项模拟方案所有因素的综合得分；
然后利用式(III)计算各模拟方案所有因素的综合得分的平均值作为评价界限的参考值：
其中，为各模拟方案所有因素的综合得分的平均值；cj为第j项模拟方案所有因素的综合得分；j＝1，2，……，m；m为模拟方案总数。
上述的评价方法中，优选地，根据实际方案评价矩阵和各因素的修正权重，计算实际方案所有因素的综合得分并排序的步骤包括：
首先将计算得到的各因素的修正权重组成修正权重向量，计算修正权重向量与实际方案评价矩阵的乘积，得到实际方案的得分向量；利用式(IV)计算实际方案所有因素的综合得分Cj：
其中，Ai为第i项因素的修正权重；i＝1，2，……，n(n≤N)；N为方案评价所需因素总数；n为现场资料实际提供的因素总数；Bji为第j项实际方案的第i项因素对应的隶属度；j＝1，2，……，M；M为实际方案总数；Cj为第j项实际方案所有因素综合得分；
然后对实际方案所有因素综合得分Cj大小进行排序。
本发明提供的基于吞吐后开井含水率的稠油边底水多轮CO2吞吐评价方法通过油藏数值模拟穷举法涵盖了评价目标的考察范围和可能数据参数，建立了科学的隶属度计算方法，得到了基于吞吐后开井含水率的评价矩阵，提出了合理的评价界限；能够解决现场资料不全带来的权重难赋值、评价矩阵难计算、评价界限缺乏依据、基于吞吐后开井含水率效果的评价结果不合理等问题，能够综合反映CO2吞吐在控水和增油两个方面的作用机制，对指导稠油边底水油藏多轮CO2吞吐评价具有重要的现实意义。
具体实施方式
为了对本发明的技术特征、目的和有益效果有更加清楚的理解，现对本发明的技术方案进行以下详细说明，但不能理解为对本发明的可实施范围的限定。
实施例
本实施例提供一种基于吞吐后开井含水率的稠油边底水多轮CO2吞吐的评价方法，对动、静态资料较全的37口水平井进行基于吞吐后开井含水率的稠油边底水多轮CO2吞吐评价，包括以下步骤：
步骤一，测定实际油藏各因素的实际参数，以此为依据建立油藏模型，模拟得到吞吐后开井含水率；
步骤二，以吞吐后开井含水率为评价关联对象，采用式(I)计算各因素的修正权重Ai：
其中，ai为第i项因素的经验权重；Ai为第i项因素的修正权重；i＝1，2，……， n(n≤N)；N为方案评价所需的因素总数；n为现场资料实际提供的因素总数。
表1为根据实际水平井现场动、静态资料，计算并修正后得到的各因素的修正权重。
表1
| 序号 | 影响因素名称/单位 | 修正权重 |
| 1 | 地层倾角/° | 0.0029 |
| 2 | 隔夹层位置 | 0.0173 |
| 3 | 油层有效厚度/m | 0.0846 |
| 4 | 沉积韵律 | 0.00173 |
| 5 | 储层非均质性(洛伦兹系数V<sub>k</sub>) | 0.0961 |
| 6 | 投产含水/f | 0.0154 |
| 7 | 产水段长度/m | 0.0154 |
| 8 | 产水段位置 | 0.0250 |
| 9 | 产水段非均质程度(渗透率极差K<sub>v</sub>) | 0.0123 |
| 10 | 水平段位置 | 0.0211 |
| 11 | 吞吐时机(含水率)/f | 0.0317 |
| 12 | 周期注气量/sm<sup>3</sup> | 0.0721 |
| 13 | 注气速度/(sm<sup>3</sup>/d) | 0.0135 |
| 14 | 焖井时间/d | 0.0394 |
| 15 | 开井后采液速度/(rm<sup>3</sup>/d) | 0.1451 |
| 16 | 平行井距离/m | 0.1451 |
| 17 | 平行井采液速度/(rm<sup>3</sup>/d) | 0.1614 |
| 18 | 吞吐轮次 | 0.096l |
步骤三，借助油藏数值模拟技术建立各因素隶属度计算方法，分别得到模拟方案和实际方案的评价矩阵；具体步骤为：
(1)根据测得的实际油藏各因素的实际参数，对各因素设计合理的参数变化范围(根据油田实际情况中各因素参数分布范围，并扩大±20％来设计)；表2为根据油田实际现场资料设计的各因素合理参数变化范围。
表2
(2)借助油藏数值模拟技术，预测各因素参数变化时对应的吞吐后开井含水率；借助ECLIPSE或者CMG等数值模拟专业软件，以实际油藏特征参数为依据，建立油藏模型，模拟油田开采，从而预测油藏采收率及开发动态数据，预测各因素参数变化时对应的吞吐后开井含水率。
(3)依据吞吐后开井含水率的大小，建立各因素与隶属度的函数关系式；
建立的所述各因素与隶属度的函数关系式分别如下：
利用式(1)得到隔夹层位置对应的隶属度y1a、y1b和y1c：
利用式(2)得到沉积韵律对应的隶属度y2a、y2b和y2c：
利用式(3)得到出水段长度对应的隶属度y3a、y3b和y3c：
利用式(4)得到出水段位置对应的隶属度y4a、y4b和y4c：
其中，A点为水平井的趾端、B点为水平井的跟端；
利用式(5)得到水平段位置对应的隶属度y5a、y5b和y5c：
利用式(6)得到平行井距离对应的隶属度y6a、y6b和y6c：
利用式(7)得到吞吐轮次对应的隶属度y7：
利用式(8)计算地层倾角对应的隶属度y8a、y8b和y8c：
其中，θ为地层倾角，°；
利用式(9)计算油层有效厚度对应隶属度y9a、y9b和y9c：
其中，h为油层有效厚度，m；
利用式(10)计算地层非均质性对应的隶属度y10a、y10b和y10c：
其中，Vk为洛伦兹系数；
利用式(11)计算投产含水对应的隶属度y11a、y11b和y11c：
其中，fw0为投产含水；
利用式(12)出水段非均质程度对应的隶属度y12a、y12b和y12c：
其中，Kv为渗透率极差；
利用式(13)计算吞吐时机对应的隶属度y13a、y13b和y13c：
其中，f为吞吐时机；
利用式(14)计算周期注气量对应的隶属度y14a、y14b和y14c：
其中，Qig为周期注气量，sm3；
利用式(15)计算注气速度对应的隶属度y15a、y15b和y15c：
其中，vig为注气速度，sm3/d；
利用式(16)计算焖井时间对应的隶属度y16a、y16b和y16c：
其中，t为焖井时间，d；
利用式(17)计算开井后采液速度对应的隶属度y17a、y17b和y17c：
其中，vl为开井后采液速度，rm3/d；
利用式(18)计算平行井采液速度对应的隶属度y18a、y18b和y18c：
其中，ql为平行井采液速度，rm3/d。
由式(1)至式(18)计算得到各因素对应的不同参数的隶属度。
(4)依据各因素对应的隶属度，得到模拟方案评价矩阵；
(5)依据建立的隶属度计算方法，计算实际方案各因素取值对应的隶属度，得到实际方案评价矩阵。
步骤四，根据模拟方案评价矩阵和各因素的修正权重，计算各模拟方案所有因素的综合得分的平均值；其具体方法为：
首先将计算得到的各因素的修正权重组成修正权重向量，计算修正权重向量与模拟方案评价矩阵的乘积，得到模拟方案的得分向量；利用式(II)计算模拟方案所有因素的综合得分cj：
其中，Ai为第i项因素的修正权重；i＝1，2，......，n(n≤N)；N为方案评价所需因素总数；n为现场资料实际提供的因素总数；bji为第j项模拟方案的第i项因素对应的隶属度；j＝1，2，......，m；m为模拟方案总数；cj为第j项模拟方案所有因素的综合得分；
然后利用式(III)计算各模拟方案所有因素的综合得分的平均值作为评价界限的参考值：
其中，为各模拟方案所有因素的综合得分的平均值；cj为第j项模拟方案所有因素的综合得分；j＝1，2，……，m；m为模拟方案总数。
步骤五，根据实际方案评价矩阵和各因素的修正权重，计算实际方案所有因素的综合得分Cj并排序；其具体方法为：
首先将计算得到的各因素的修正权重组成修正权重向量，计算修正权重向量与实际方案评价矩阵的乘积，得到实际方案的得分向量；利用式(IV)计算实际方案所有因素的综合得分Cj：
其中，Ai为第i项因素的修正权重；i＝1，2，……，n(n≤N)；N为方案评价所需因素总数；n为现场资料实际提供的因素总数；Bji为第j项实际方案的第i项因素对应的隶属度；j＝1，2，……，M；M为实际方案总数；Cj为第j项实际方案所有因素综合得分；
然后对实际方案所有因素综合得分大小进行排序。
计算结果如表3和表4所示。表3为依据各因素对应的隶属度，得到模拟方案评价矩阵，以及修正权重向量与模拟方案评价矩阵的乘积后计算得到的各模拟方案所有因素的综合得分cj和各模拟方案所有因素的综合得分的平均值(评价界限)。
表4为依据所建立的隶属度计算方法，计算实际方案各因素取值对应的隶属度所得到的实际方案评价矩阵，以及修正权重向量与实际方案评价矩阵的乘积后计算得到的实际方案所有因素的综合得分Cj和排序结果。
步骤六，对比和Cj的大小，给出评价结论；若则评价结果为较好；若则评价结果为较差，由表4对比结果可知，该油田实际的37口稠油边底水油藏多轮CO2吞吐评价结果为：排序24号之前的井吞吐效果较好，排序24号以后的井的吞吐效果较差。
综上所述，本发明提供的基于吞吐后开井含水率的稠油边底水多轮CO2吞吐评价方法通过油藏数值模拟穷举法涵盖了评价目标的考察范围和可能数据参数，建立了科学的隶属度计算方法，得到了基于吞吐后开井含水率的评价矩阵，提出了合理的评价界限；能够解决现场资料不全带来的权重难赋值、评价矩阵难计算、评价界限缺乏依据、基于吞吐后开井含水率效果的评价结果不合理等问题，能够综合反映CO2吞吐在控水和增油两个方面的作用机制，对指导稠油边底水油藏多轮CO2吞吐评价具有重要的现实意义。
Claims (2)
1.一种基于吞吐后开井含水率的稠油边底水多轮CO2吞吐的评价方法，其包括以下步骤：
步骤一，测定实际油藏各因素的实际参数，以此为依据建立油藏模型，模拟得到吞吐后开井含水率；
步骤二，以吞吐后开井含水率为评价关联对象，计算各因素的修正权重；
步骤三，借助油藏数值模拟技术建立各因素与隶属度的函数关系式，分别得到模拟方案和实际方案的评价矩阵；
步骤四，根据模拟方案评价矩阵和各因素的修正权重，计算各模拟方案所有因素的综合得分的平均值
步骤五，根据实际方案评价矩阵和各因素的修正权重，计算第j项实际方案所有因素的综合得分Cj并排序；
步骤六，对比和Cj的大小，给出评价结论；若则评价结果为较好；若则评价结果为较差；
借助油藏数值模拟技术建立各因素与隶属度的函数关系式，分别得到模拟方案和实际方案的评价矩阵的步骤包括：
根据测得的实际油藏各因素的实际参数，对各因素设计参数变化范围；
借助油藏数值模拟技术，预测各因素参数变化时对应的吞吐后开井含水率；
依据吞吐后开井含水率的大小，建立各因素与隶属度的函数关系式；
依据各因素对应的隶属度，得到模拟方案评价矩阵；
依据建立的各因素与隶属度的函数关系式，计算实际方案各因素取值对应的隶属度，得到实际方案评价矩阵；
所述各因素包括隔夹层位置、沉积韵律、出水段长度、出水段位置、水平段位置、平行井距离、吞吐轮次、地层倾角、油层有效厚度、地层非均质性、投产含水、出水段非均质程度、吞吐时机、周期注气量、注气速度、焖井时间、开井后采液速度和平行井产液速度中的多种；
所述各因素与隶属度的函数关系式分别如下：
利用式(1)得到隔夹层位置对应的隶属度y1a、y1b和y1c：
利用式(2)得到沉积韵律对应的隶属度y2a、y2b和y2c：
利用式(3)得到出水段长度对应的隶属度y3a、y3b和y3c：
利用式(4)得到出水段位置对应的隶属度y4a、y4b和y4c：
其中，A点为水平井的趾端、B点为水平井的跟端；
利用式(5)得到水平段位置对应的隶属度y5a、y5b和y5c：
利用式(6)得到平行井距离对应的隶属度y6a、y6b和y6c：
利用式(7)得到吞吐轮次对应的隶属度y7：
利用式(8)计算地层倾角对应的隶属度y8a、y8b和y8c：
第1轮：y8a＝-0.0065θ2-0.0582θ+1.233
第2轮：y8b＝-0.0076θ2-0.242θ+1.657
第3轮：y8c＝-0.0065θ2-0.0582θ+1.233 式(8)
其中，θ为地层倾角，°；
利用式(9)计算油层有效厚度对应隶属度y9a、y9b和y9c：
第1轮：y9a＝-0.0826h+1.091
第2轮：y9b＝-0.0826h+1.091
第3轮：y9c＝-0.0826h+1.091 式(9)
其中，h为油层有效厚度，m；
利用式(10)计算地层非均质性对应的隶属度y10a、y10b和y10c：
其中，Vk为洛伦兹系数；
利用式(11)计算投产含水对应的隶属度y11a、y11b和y11c：
第1轮：y11a＝-fw0+1
第2轮：y11b＝-fw0+1
第3轮：y11c＝-fw0+1 式(11)
其中，fw0为投产含水；
利用式(12)计算出水段非均质程度对应的隶属度y12a、y12b和y12c：
其中，Kv为渗透率极差；
利用式(13)计算吞吐时机对应的隶属度y13a、y13b和y13c：
其中，f为吞吐时机；
利用式(14)计算周期注气量对应的隶属度y14a、y14b和y14c：
其中，Qig为周期注气量，sm3；
利用式(15)计算注气速度对应的隶属度y15a、y15b和y15c：
其中，vig为注气速度，sm3/d；
利用式(16)计算焖井时间对应的隶属度y16a、y16b和y16c：
其中，t为焖井时间，d；
利用式(17)计算开井后采液速度对应的隶属度y17a、y17b和y17c：
其中，vl为开井后采液速度，rm3/d；
利用式(18)计算平行井采液速度对应的隶属度y18a、y18b和y18c：
其中，ql为平行井采液速度，rm3/d；
根据模拟方案评价矩阵和各因素的修正权重，计算各模拟方案所有因素的综合得分的平均值的步骤包括：
首先将计算得到的各因素的修正权重组成修正权重向量，计算修正权重向量与模拟方案评价矩阵的乘积，得到模拟方案的得分向量；利用式(II)计算第j项模拟方案所有因素的综合得分cj：
其中，Ai为第i项因素的修正权重；i＝1，2，……，n(n≤N)；N为方案评价所需因素总数；n为现场资料实际提供的因素总数；bji为第j项模拟方案的第i项因素对应的隶属度；j＝1，2，……，m；m为模拟方案总数；cj为第j项模拟方案所有因素的综合得分；
然后利用式(III)计算各模拟方案所有因素的综合得分的平均值作为评价界限的参考值：
其中，为各模拟方案所有因素的综合得分的平均值；cj为第j项模拟方案所有因素的综合得分；j＝1，2，……，m；m为模拟方案总数；
根据实际方案评价矩阵和各因素的修正权重，计算实际方案所有因素的综合得分并排序的步骤包括：
首先将计算得到的各因素的修正权重组成修正权重向量，计算修正权重向量与实际方案评价矩阵的乘积，得到实际方案的得分向量；利用式(IV)计算第j项实际方案所有因素的综合得分Cj：
其中，Ai为第i项因素的修正权重；i＝1，2，……，n(n≤N)；N为方案评价所需因素总数；n为现场资料实际提供的因素总数；Bji为第j项实际方案的第i项因素对应的隶属度；j＝1，2，……，M；M为实际方案总数；Cj为第j项实际方案所有因素的综合得分；
然后对第j项实际方案所有因素的综合得分按Cj大小进行排序。
2.根据权利要求1所述的评价方法，其特征在于，各因素的修正权重是采用式(I)计算得到的：
其中，ai为第i项因素的经验权重；Ai为第i项因素的修正权重；i＝1，2，……，n(n≤N)；N为方案评价所需的因素总数；n为现场资料实际提供的因素总数。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710003679.3A CN106802985B (zh) | 2017-01-03 | 2017-01-03 | 基于含水率的稠油边底水油藏多轮co2吞吐的评价方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710003679.3A CN106802985B (zh) | 2017-01-03 | 2017-01-03 | 基于含水率的稠油边底水油藏多轮co2吞吐的评价方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106802985A CN106802985A (zh) | 2017-06-06 |
| CN106802985B true CN106802985B (zh) | 2019-09-24 |
Family
ID=58985325
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710003679.3A Active CN106802985B (zh) | 2017-01-03 | 2017-01-03 | 基于含水率的稠油边底水油藏多轮co2吞吐的评价方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106802985B (zh) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109372484B (zh) * | 2017-08-11 | 2021-05-25 | 中国石油化工股份有限公司 | 致密油藏co2吞吐有效作用半径的实验预测方法及*** |
| CN111946323B (zh) * | 2019-04-30 | 2023-10-31 | 中国石油天然气股份有限公司 | 进行增能吞吐试验的油井选择方法及装置 |
| CN111677486A (zh) * | 2020-05-27 | 2020-09-18 | 中国石油大学(北京) | 致密油二氧化碳吞吐模拟方法、装置及存储介质 |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105781505A (zh) * | 2016-03-22 | 2016-07-20 | 中国石油大学(北京) | 一种边水油藏水平井组注气吞吐的物理模拟方法 |
-
2017
- 2017-01-03 CN CN201710003679.3A patent/CN106802985B/zh active Active
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105781505A (zh) * | 2016-03-22 | 2016-07-20 | 中国石油大学(北京) | 一种边水油藏水平井组注气吞吐的物理模拟方法 |
Non-Patent Citations (3)
| Title |
|---|
| G104油藏二氧化碳吞吐物模参数优选及数值模拟研究;李亚辉;《中国优秀硕士学位论文全文数据库工程科技Ⅰ辑》;20150315(第03期);全文 * |
| 复杂断块油藏CO2吞吐试验及效果分析;谈士海，周正平，刘伟，等;《石油钻采工艺》;20020825;第24卷(第4期);全文 * |
| 榆树林油田CO2吞吐采油目标井优选方法研究;战菲，宋考平，秦亮，等;《数学的实践与认识》;20011108;第40卷(第21期);全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106802985A (zh) | 2017-06-06 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106802985B (zh) | 基于含水率的稠油边底水油藏多轮co2吞吐的评价方法 | |
| CN105587297B (zh) | 复杂断块油藏仿强边水驱技术适应性定量评价方法 | |
| US20190251460A1 (en) | Method for predicting oil and gas reservoir production | |
| CN104091069B (zh) | 确定非均质储层各层位和位置驱油效率和波及系数的方法 | |
| Zendehboudi et al. | Estimation of breakthrough time for water coning in fractured systems: Experimental study and connectionist modeling | |
| CN105740563B (zh) | 一种成熟油田二次开发之优势通道识别方法 | |
| CN108240208A (zh) | 一种油田水驱分类井组开发效果对标方法 | |
| CN105045941A (zh) | 基于无迹卡尔曼滤波的抽油机参数优化方法 | |
| CN107085942A (zh) | 一种基于狼群算法的交通流预测方法、装置及*** | |
| CN107730151A (zh) | 一种基于概念性水文模型的流域设计洪水推求方法 | |
| CN106909731B (zh) | 一种膨胀性土壤产汇流过程模拟方法 | |
| CN109034647A (zh) | 一种致密油气藏体积压裂水平井重复压裂选井的方法 | |
| CN107038516A (zh) | 一种中渗复杂断块油藏水驱开发效果定量评价方法 | |
| CN108509601A (zh) | 一种基于大数据分析的食品风味评定方法 | |
| CN114021466A (zh) | 基于返排数据和机器学习的页岩气有效缝网体积预测方法 | |
| CN109063403B (zh) | 一种滑溜水压裂优化设计方法 | |
| CN104680025B (zh) | 基于遗传算法极限学习机的抽油机参数优化方法 | |
| CN106802986B (zh) | 一种稠油边底水油藏多轮co2吞吐评价方法 | |
| CN106372320B (zh) | 一种采用亚滤波尺度模型对公路隧道湍流进行大涡模拟的方法 | |
| CN109469478A (zh) | 基于改进qfd的高含水期整装油田多级开发潜力评价方法 | |
| CN106920173B (zh) | 一种稠油边底水油藏多轮co2吞吐的选井方法 | |
| CN106055871A (zh) | 一种气井泡沫排水采气效果评价方法 | |
| CN104680023B (zh) | 基于多目标决策的抽油机参数优化方法 | |
| CN104268662A (zh) | 一种基于分步优化分位数回归的沉降预测方法 | |
| CN104809363B (zh) | 基于bk近似动态规划的三元复合驱优化方法 |
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 |