本发明属于稠油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.00292隔夹层位置0.01733油层有效厚度/m0.08464沉积韵律0.001735储层非均质性(洛伦兹系数Vk)0.09616投产含水/f0.01547产水段长度/m0.01548产水段位置0.02509产水段非均质程度(渗透率极差Kv)0.012310水平段位置0.021111吞吐时机(含水率)/f0.031712周期注气量/sm30.072113注气速度/(sm3/d)0.013514焖井时间/d0.039415开井后采液速度/(rm3/d)0.145116平行井距离/m0.145117平行井采液速度/(rm3/d)0.161418吞吐轮次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吞吐评价具有重要的现实意义。当前第1页1 2 3