CN111914209A - Drainage gas production effect fuzzy comprehensive evaluation method based on entropy method - Google Patents
Drainage gas production effect fuzzy comprehensive evaluation method based on entropy method Download PDFInfo
- Publication number
- CN111914209A CN111914209A CN202010399587.3A CN202010399587A CN111914209A CN 111914209 A CN111914209 A CN 111914209A CN 202010399587 A CN202010399587 A CN 202010399587A CN 111914209 A CN111914209 A CN 111914209A
- Authority
- CN
- China
- Prior art keywords
- evaluation
- gas production
- entropy
- index
- drainage
- 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.)
- Pending
Links
- 238000011156 evaluation Methods 0.000 title claims abstract description 161
- 238000004519 manufacturing process Methods 0.000 title claims abstract description 67
- 230000000694 effects Effects 0.000 title claims abstract description 45
- 238000000034 method Methods 0.000 title claims abstract description 38
- 239000006260 foam Substances 0.000 claims abstract description 35
- 238000011084 recovery Methods 0.000 claims abstract description 30
- 239000011159 matrix material Substances 0.000 claims abstract description 22
- 238000004458 analytical method Methods 0.000 claims abstract description 6
- 230000008859 change Effects 0.000 claims description 22
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 11
- 238000010606 normalization Methods 0.000 claims description 7
- 230000015572 biosynthetic process Effects 0.000 claims description 3
- 238000012545 processing Methods 0.000 claims description 3
- 238000003786 synthesis reaction Methods 0.000 claims description 2
- 230000005484 gravity Effects 0.000 claims 1
- 239000000243 solution Substances 0.000 description 8
- 238000004364 calculation method Methods 0.000 description 6
- 230000008569 process Effects 0.000 description 4
- 230000008901 benefit Effects 0.000 description 3
- 238000010276 construction Methods 0.000 description 3
- 238000013079 data visualisation Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 239000012530 fluid Substances 0.000 description 2
- 239000004088 foaming agent Substances 0.000 description 2
- 238000009825 accumulation Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 239000003795 chemical substances by application Substances 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000000265 homogenisation Methods 0.000 description 1
- 238000002347 injection Methods 0.000 description 1
- 239000007924 injection Substances 0.000 description 1
- 239000007788 liquid Substances 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
Images
Classifications
-
- 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/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
-
- 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
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Physics (AREA)
- Computational Mathematics (AREA)
- Mathematical Optimization (AREA)
- Data Mining & Analysis (AREA)
- Pure & Applied Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- Mathematical Analysis (AREA)
- Computing Systems (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Algebra (AREA)
- Geometry (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- Feedback Control In General (AREA)
Abstract
本发明公开了一种基于熵值法的排水采气效果模糊综合评价方法,针对泡沫排水采气井的特性,选取评价泡沫排水采气效果的指标,遵循定性与定量相结合的原则,构建出泡沫排水采气效果评价指标体系,考虑气井排水采气效果评价指标之间的相互关系,采用线性分析方法建立隶属度矩阵,采用熵值法确定指标权重,结合模糊综合评价法计算出气井泡沫排水采气效果的综合评价指数,利用本发明的方法可得到更为科学、准确和客观的评价信息。
The invention discloses a fuzzy comprehensive evaluation method of drainage gas production effect based on an entropy method. According to the characteristics of foam drainage gas production wells, an index for evaluating foam drainage gas production effect is selected, and a foam is constructed by following the principle of combining qualitative and quantitative. The evaluation index system of drainage gas recovery effect, considering the relationship between the evaluation indicators of gas well drainage gas recovery effect, adopts linear analysis method to establish membership degree matrix, adopts entropy value method to determine the index weight, and calculates the foam drainage recovery rate of gas well with fuzzy comprehensive evaluation method. The comprehensive evaluation index of the gas effect can be obtained by using the method of the present invention to obtain more scientific, accurate and objective evaluation information.
Description
技术领域technical field
本发明涉及排水采气效果评价技术领域,特别是涉及一种基于熵值法的排水采气效果模糊综合评价方法。The invention relates to the technical field of drainage gas recovery effect evaluation, in particular to a fuzzy comprehensive evaluation method for drainage gas recovery effect based on an entropy method.
背景技术Background technique
排水采气是解决气井井筒及井底附近地层积液过多或产水,并使气井恢复正常生产,是气井生产过程中一项提升生产效率的重要措施。泡沫排水采气因其具有设备简单、容易施工、见效快、成本低、不影响气井生产的优点,在全球的排水采气作业生产中得到广泛应用,因而系统地评价泡沫排水采气效果就显得尤为重要。Drainage and gas recovery is an important measure to improve the production efficiency in the process of gas well production to solve the problem of excessive liquid accumulation or water production near the wellbore and the bottom of the gas well, and to restore the gas well to normal production. Foam drainage gas recovery has been widely used in drainage gas recovery operations around the world because of its simple equipment, easy construction, quick effect, low cost, and no impact on gas well production. especially important.
然而,影响泡沫排水采气工艺效果的因素分析较多,从而导致泡沫排水采气领域并没有一套完整且获得普遍行业认可的评价体系,无法对不同气田的出水气井的泡沫排水效果进行对比评价,因此基于气井产气领域的生产技术管理需要,本文提出一种基于熵值法与模糊综合评价的排水采气综合评价方法。However, there are many factors that affect the effect of foam drainage gas recovery process, resulting in the absence of a complete and generally industry-recognized evaluation system in the field of foam drainage gas recovery, and it is impossible to compare and evaluate the foam drainage effect of water-producing gas wells in different gas fields. Therefore, based on the needs of production technology management in the field of gas production from gas wells, this paper proposes a comprehensive evaluation method for drainage and gas production based on entropy method and fuzzy comprehensive evaluation.
发明内容SUMMARY OF THE INVENTION
本发明的目的是针对现有技术中存在的技术缺陷,而提供一种基于熵值法的排水采气效果模糊综合评价方法。The purpose of the present invention is to provide a fuzzy comprehensive evaluation method of drainage gas recovery effect based on entropy method, aiming at the technical defects existing in the prior art.
为实现本发明的目的所采用的技术方案是:The technical scheme adopted for realizing the purpose of the present invention is:
一种基于熵值法的排水采气效果模糊综合评价方法,包括以下步骤:A fuzzy comprehensive evaluation method for drainage gas recovery effect based on entropy method, comprising the following steps:
步骤1,选取n个影响排水采气效果的评价指标Δui,i=1,2,…,n;Step 1, select n evaluation indexes Δu i that affect the effect of drainage and gas recovery, i=1, 2, ..., n;
步骤2,对应排水采气效果的综合评价指数D与评价级别的关系,将综合评价指数D设定为m个取值区间,vj为取值区间的中值,j=1,2,…,m,由此确定评判集V={v1,v2,v3,…vm};Step 2: Corresponding to the relationship between the comprehensive evaluation index D of the drainage and gas recovery effect and the evaluation level, the comprehensive evaluation index D is set as m value intervals, v j is the median value of the value interval, j=1, 2, ... , m, thus determine the judgment set V={v 1 , v 2 , v 3 ,...v m };
步骤3,采用线性分析方法建立隶属度矩阵R:Step 3, use the linear analysis method to establish the membership matrix R:
首先,将每个评价指标划分成m个评价区间,并设定每个评价指标的评价区间对应的级别界限值aj,j=1,2,…,m;First, each evaluation index is divided into m evaluation intervals, and the level limit value a j corresponding to the evaluation interval of each evaluation index is set, j=1, 2, ..., m;
然后,计算每一评价指标的评价区间相对于评判集的隶属度进而可求出隶属度构成的矩阵,也就是模糊关系矩阵R;Then, calculate the membership degree of the evaluation interval of each evaluation index relative to the evaluation set Then, the matrix composed of membership degrees, that is, the fuzzy relationship matrix R, can be obtained;
步骤4,采用熵值法确定各评价指标的权重W:Step 4, using the entropy method to determine the weight W of each evaluation index:
首先,对步骤3求算得到的隶属度转化为相对值,令然后进行归一化处理,得到tij;First, calculate the membership degree obtained in step 3 Converted to relative values, let Then perform normalization to obtain t ij ;
然后计算第j项评价区间下第i个影响因素tij占该评价区间的比重,以及第j项评价区间的熵值;进而计算各项评价指标的权值,再得到各项评价指标权重W;Then calculate the proportion of the ith influencing factor tij under the jth evaluation interval in the evaluation interval, and the entropy value of the jth evaluation interval; then calculate the weight of each evaluation index, and then obtain the weight of each evaluation index W ;
步骤5,采用(+,*)算子模糊合成步骤4得到的各项评价指标权重W与步骤3得到的模糊关系矩阵R,得到最终评价向量S,再由最终评价向量S计算气井泡沫排水采气效果的综合评价指数D。Step 5: Use the (+, *) operator to fuzzy synthesize the weights of each evaluation index obtained in step 4 and the fuzzy relationship matrix R obtained in step 3 to obtain the final evaluation vector S, and then calculate the gas well foam drainage recovery from the final evaluation vector S. The comprehensive evaluation index D of the gas effect.
在上述技术方案中,所述步骤1中,n=4,所述评价指标分别为日产气变化率Δqg、日产水变化率Δqw、油套压差变化率ΔP以及日泡沫排水采气作业成本ΔC。In the above technical solution, in the step 1, n=4, and the evaluation indicators are the daily gas production change rate Δq g , the daily water production change rate Δq w , the oil jacket pressure difference change rate ΔP and the daily foam drainage gas production operation. Cost ΔC.
在上述技术方案中,所述步骤2中,m=5,所述评判集V={v1,v2,v3,v4,v5},评价级别为很好、好、较好、一般和差,v1=90,v2=70,v3=50,v4=30,v5=10。In the above technical solution, in the
在上述技术方案中,所述步骤3中,In the above technical solution, in the step 3,
其中,求出构成的矩阵得到模糊关系矩阵R。Among them, find The formed matrix gets the fuzzy relation matrix R.
在上述技术方案中,所述步骤3中,归一化处理公式如下:In the above technical solution, in the step 3, the normalization processing formula is as follows:
正向指标:Positive indicators:
负向指标:Negative indicators:
则r′ij或r″ij为第i个影响因素的第j个评价区间的数值(i=1,2,…,n;j=1,2,…,m),归一化后的数据记为tij。Then r′ ij or r″ ij is the value of the j-th evaluation interval of the i-th influencing factor (i=1, 2,...,n; j=1, 2,...,m), the normalized data Denoted as t ij .
在上述技术方案中,所述步骤4中,tij占该评价区间的比重第j项评价区间的熵值 In the above technical solution, in the step 4, tij accounts for the proportion of the evaluation interval The entropy value of the jth evaluation interval
在上述技术方案中,所述步骤4中,计算信息熵冗余度:dj=l-ej,然后计算各项指标的权值计算得出权重,确定权重W:W=(w1,w2,…,wn)式中权重和为 1。In the above technical solution, in the step 4, the information entropy redundancy is calculated: d j =le j , and then the weights of each index are calculated The weight is obtained by calculation, and the weight W is determined: W=(w 1 , w 2 , . . . , wn ) where the weight sum is 1.
在上述技术方案中,所述步骤5中,S=W·R=(s1,s2,...,sm)。In the above technical solution, in the step 5, S=W·R=(s 1 , s 2 , . . . , s m ).
在上述技术方案中,所述步骤5中,D=S*V=s1*v1+s2*v2+…+sm*vm。In the above technical solution, in the step 5, D=S*V=s 1 *v 1 +s 2 *v 2 +...+s m *v m .
与现有技术相比,本发明的有益效果是:Compared with the prior art, the beneficial effects of the present invention are:
1.求权重时用熵值法可以弥补主成分用方差贡献率做权重的不足,解决实际存在的等级界限模糊性,避免了人为主观因素的干扰。并且此方法考虑到了数据的离散程度,可以根据指标对综合评价的影响程度进行指标权重的调整。1. The entropy method can make up for the insufficiency of using the variance contribution rate as the weight of the principal component when calculating the weight, solve the ambiguity of the actual level boundary, and avoid the interference of human subjective factors. And this method takes into account the discrete degree of data, and can adjust the index weight according to the influence degree of the index on the comprehensive evaluation.
2.从评价结果反映问题的多少来看,两种方法结合的模型既可以得出具有实际意义的新指标,又可以反映出每一个主成分的主要指标,还可以对各项指标的综合情况进行排序,使复杂问题简化,同时得到更为科学、准确和客观的评价信息。2. Judging from how much the evaluation results reflect the problem, the model combined with the two methods can not only draw new indicators with practical significance, but also reflect the main indicators of each principal component, and can also evaluate the comprehensive situation of each indicator. Sorting can simplify complex problems and obtain more scientific, accurate and objective evaluation information.
附图说明Description of drawings
图1为本发明涉及的基于熵值法的排水采气效果模糊综合评价方法流程图;Fig. 1 is the flow chart of the fuzzy comprehensive evaluation method of drainage gas recovery effect based on entropy value method involved in the present invention;
图2为本发明涉及的泡沫排水采气作业前后产气量对比图;Fig. 2 is the comparison chart of gas production before and after the foam drainage gas production operation involved in the present invention;
图3为本发明涉及的泡沫排水采气作业前后产水量对比图;Fig. 3 is a comparison diagram of water production before and after the foam drainage gas production operation involved in the present invention;
图4为本发明涉及的泡沫排水采气作业前后油套压差值对比图;Fig. 4 is a comparison diagram of oil jacket pressure difference before and after the foam drainage gas production operation involved in the present invention;
图5为本发明涉及的基于熵值法的泡沫排水采气作业综合评价指数和油套压差变化率曲线图。FIG. 5 is a graph showing the comprehensive evaluation index of the foam drainage gas recovery operation and the change rate of the oil jacket pressure difference based on the entropy method according to the present invention.
具体实施方式Detailed ways
以下结合具体实施例对本发明作进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。The present invention will be further described in detail below in conjunction with specific embodiments. It should be understood that the specific embodiments described herein are only used to explain the present invention, but not to limit the present invention.
实施例1Example 1
如附图1所示,一种基于熵值法的排水采气效果模糊综合评价方法,包括以下步骤:As shown in accompanying drawing 1, a kind of fuzzy comprehensive evaluation method of drainage gas recovery effect based on entropy value method, comprises the following steps:
步骤1,针对泡沫排水采气井的特性,选取评价泡沫排水采气效果的评价指标;Step 1, according to the characteristics of the foam drainage gas production well, select the evaluation index for evaluating the foam drainage gas production effect;
步骤2,遵循定性与定量相结合的原则,构建出泡沫排水采气效果评价指标体系:
对应综合评价指数与评价级别的关系,设定m个取值区间vj,j=1,2,…,m;由此确定评判集V={v1,v2,v3,…vm};Corresponding to the relationship between the comprehensive evaluation index and the evaluation level, set m value intervals v j , j = 1, 2, ..., m; thus determine the evaluation set V = {v 1 , v 2 , v 3 , ... v m };
步骤3,考虑气井排水采气效果评价指标之间的相互关系,采用线性分析方法建立隶属度矩阵,得到模糊关系矩阵R:Step 3: Considering the relationship between the evaluation indicators of gas well drainage and gas recovery, a linear analysis method is used to establish a membership matrix, and a fuzzy relationship matrix R is obtained:
3.1,将每个评价指标划分成m个评价区间,并设定每个评价指标的评价区间对应的级别界限值aj,j=1,2,…,m;3.1. Divide each evaluation index into m evaluation intervals, and set the level limit values a j , j=1, 2, ..., m corresponding to the evaluation interval of each evaluation index;
3.2,计算每一评价指标的评价区间相对于评判集的隶属度进而可求出隶属度构成的矩阵,也就是模糊关系矩阵R;3.2, Calculate the membership degree of the evaluation interval of each evaluation index relative to the evaluation set Then, the matrix composed of membership degrees, that is, the fuzzy relationship matrix R, can be obtained;
步骤4,采用熵值法确定各评价指标的权重:Step 4, using the entropy method to determine the weight of each evaluation index:
4.1,对步骤3求算得到的隶属度转化为相对值,令然后进行归一化处理,得到tij;4.1, the membership degree calculated in step 3 Converted to relative values, let Then perform normalization to obtain t ij ;
4.2,计算第j项评价区间下第i个影响因素tij占该评价区间的比重,以及第j项评价区间的熵值;进而计算各项评价指标的权值,再得到各项评价指标权重W;4.2, calculate the proportion of the i-th influencing factor t ij under the j-th evaluation interval in the evaluation interval, and the entropy value of the j-th evaluation interval; then calculate the weight of each evaluation index, and then obtain the weight of each evaluation index W;
步骤5,结合模糊综合评价法计算出气井泡沫排水采气效果的综合评价指数。Step 5: Combine the fuzzy comprehensive evaluation method to calculate the comprehensive evaluation index of the foam drainage effect of the gas well.
实施例2Example 2
参照图1所示的方法流程图,基于熵值法的排水采气效果模糊综合评价方法过程如下:Referring to the method flow chart shown in Figure 1, the process of the fuzzy comprehensive evaluation method for drainage and gas recovery based on the entropy method is as follows:
步骤1,确定排水采气评价指标,通过对影响泡沫排水采气效果与效益的各项因素的分析,选取影响排水采气效果的四个重要指标:日产气变化率(Δqg)、日产水变化率(Δqw)、油套压差变化率(ΔP)以及日泡沫排水采气作业成本(ΔC)作为评价指标。将这四个评价指标组成模糊集合U,将此模糊集合为排水采气效果评价目标。U的表达式如下:Step 1: Determine the evaluation index of drainage and gas recovery. Through the analysis of various factors affecting the effect and benefit of foam drainage and gas recovery, four important indicators affecting the drainage and gas recovery effect are selected: daily gas production change rate (Δq g ), daily water production The rate of change (Δq w ), the rate of change of oil jacket pressure difference (ΔP) and the daily foam drainage gas production cost (ΔC) were used as evaluation indicators. These four evaluation indicators are formed into a fuzzy set U, and this fuzzy set is the evaluation target of drainage gas production effect. The expression for U is as follows:
U={u1,u2,u3,u4}={Δqg,Δqw,ΔP,ΔC}.U={u 1 ,u 2 ,u 3 ,u 4 }={Δq g ,Δq w ,ΔP,ΔC}.
步骤2,确定评判集V,为了对应综合评价指数与评价级别的关系,将综合评价指数D 设定为m个取值区间,vj为取值区间的中值,在本实施例中
V={v1,v2,v3,v4,v5}={很好,好,较好,一般,差}V={v 1 , v 2 , v 3 , v 4 , v 5 }={very good, good, better, fair, poor}
也就是j=1,2,3,4,5。That is, j=1, 2, 3, 4, 5.
为了对应综合评价指数与评价效果的关系,将最后计算出的综合评价指数归到该评判集,本实施例中设定综合评价指数的取值区间如下表表1所示:In order to correspond to the relationship between the comprehensive evaluation index and the evaluation effect, the final calculated comprehensive evaluation index is classified into the evaluation set. In this embodiment, the value interval of the comprehensive evaluation index is set as shown in Table 1 below:
表1取值区间Table 1 Value range
vj为对应综合评价指数取值区间的中值。v j is the median value of the corresponding comprehensive evaluation index value interval.
也就是,当D大于等于80时小于100时,评价级别为很好,当D大于等于60小于80时,评价级别为好,当D大于等于40小于60时,评价级别为较好,当D大于等于20小于 40时,评价级别为一般,当D大于等于0小于20时,评价级别为差。That is, when D is greater than or equal to 80 and less than 100, the evaluation level is very good, when D is greater than or equal to 60 and less than 80, the evaluation level is good, when D is greater than or equal to 40 and less than 60, the evaluation level is good, and when D is greater than or equal to 40 and less than 60, the evaluation level is good. When D is greater than or equal to 20 and less than 40, the evaluation level is general, and when D is greater than or equal to 0 and less than 20, the evaluation level is poor.
步骤3,为了考虑气井排水采气效果评价的评价指标之间的相互关系,采用线性分析方法建立隶属度矩阵,为了划分评价指标相对于评判集的隶属度,需要将每个评价指标的不同级别划分区间,具体划分形式如表2所示。Step 3, in order to consider the relationship between the evaluation indexes of the gas well drainage and gas recovery effect evaluation, a linear analysis method is used to establish a membership degree matrix. The specific division forms are shown in Table 2.
表2评价指标划分出的评价区间Table 2 Evaluation interval divided by evaluation index
举例说明,在表2中,日产气变化率大于等于0.5时,其评价级别为很好,日产气变化率大于等于0.2小于0.5时,其评价级别为好,日产气变化率大于等于0.1小于0.2时,其评价级别为较好,日产气变化率大于等于0.05小于0.1时,其评价级别为一般,日产气变化率大于等于0小于0.05时,其评价级别为差。For example, in Table 2, when the daily gas production change rate is greater than or equal to 0.5, the evaluation level is very good, when the daily gas production change rate is greater than or equal to 0.2 and less than 0.5, the evaluation level is good, and the daily gas production change rate is greater than or equal to 0.1 and less than 0.2 When the daily gas production change rate is greater than or equal to 0.05 and less than 0.1, the evaluation level is fair, and when the daily gas production change rate is greater than or equal to 0 and less than 0.05, the evaluation level is poor.
然后根据评价指标划分区间计算隶属度,隶属度的计算如下:Then, the membership degree is calculated by dividing the interval according to the evaluation index. The membership degree is calculated as follows:
设若将评价指标Δui划分m个级别,设定aj(j=1,2,....,m)为评价指标对应的级别界限值,评价指标对应评价级别的计算公式(即隶属度的计算公式)如下:If the evaluation index Δu i is divided into m levels, set a j (j=1, 2, . . . , m) as the level limit value corresponding to the evaluation index, and the calculation formula of the evaluation index corresponding to the evaluation level (that is, the degree of membership) The calculation formula of ) is as follows:
根据Δui与ai的关系和对应的关系式求解的值,得出各评价指标相对于各评价区间的隶属度,进而可求出构成的矩阵,也就是模糊关系矩阵R。According to the relationship between Δu i and a i and the corresponding relationship to solve , the membership degree of each evaluation index relative to each evaluation interval can be obtained, and then the degree of membership of each evaluation index can be obtained. The matrix formed is the fuzzy relation matrix R.
在本实施例中,i=1,2,3,4,也就是n=4,各评价指标划分为了5个级别,也就是 m=5,则对于u1(日产气变化率)而言,a1=0.5,a2=0.2,a3=0.1,a4=0.05,a5=0,再由隶属度的计算公式,计算得到同理计算u2(日产水变化率)、u3(油套压差变化率)、u4日加注量成本变化率中各隶属度,得到模糊关系矩阵R 如下:In this embodiment, i=1, 2, 3, 4, that is, n=4, and each evaluation index is divided into 5 levels, that is, m=5, then for u 1 (daily gas production rate of change), a 1 = 0.5, a 2 = 0.2, a 3 = 0.1, a 4 = 0.05, a 5 = 0, and then by the calculation formula of the membership degree, we can get Similarly, calculate the membership degrees in u 2 (daily water production rate of change), u 3 (oil jacket pressure difference rate of change), and u 4 daily injection cost rate of change, and obtain the fuzzy relationship matrix R as follows:
步骤4,确定影响因素权重,选取n个影响因素(本实施例中选取了四个因素,即n=4), m个评价区间(本实施例中m=5),则为第i个影响因素的第j个评价区间的数值(i=1, 2,...,n;j=1,2,...,m)。Step 4: Determine the weight of the influencing factors, select n influencing factors (four factors are selected in this embodiment, that is, n=4), and m evaluation intervals (m=5 in this embodiment), then is the value of the j-th evaluation interval of the i-th influencing factor (i=1, 2,...,n; j=1, 2,...,m).
指标的归一化处理:异质指标同质化,由于各项指标的计量单位并不统一,首先要进行标准化处理,即把指标的绝对值转化为相对值,即令从而解决各项不同质指标值的同质化问题。指标归一化式子如下:Normalization of indicators: Homogenization of heterogeneous indicators. Since the measurement units of various indicators are not uniform, standardization must be carried out first, that is, the absolute value of the indicator is converted into a relative value, that is, the So as to solve the problem of homogeneity of different index values. The index normalization formula is as follows:
正向指标:Positive indicators:
负向指标:Negative indicators:
则r′ij或r″ij为第i个影响因素的第j个评价区间的数值(i=1,2,...,n;j=1,2,...,m),归一化后的数据记为tij。Then r′ ij or r″ ij is the value of the j-th evaluation interval of the i-th influencing factor (i=1, 2,...,n; j=1, 2,...,m), normalized The transformed data is denoted as t ij .
计算第j项评价区间下第i个影响因素(tij)占该评价区间的比重:Calculate the proportion of the i-th influencing factor (t ij ) under the j-th evaluation interval in the evaluation interval:
其中i=1,...,n,j=1,...,m。where i=1,...,n,j=1,...,m.
计算第j项评价区间的熵值:Calculate the entropy value of the jth evaluation interval:
其中k=1/ln(n),满足ej≥0。Where k=1/ln(n), e j ≥ 0 is satisfied.
计算信息熵冗余度:Calculate the information entropy redundancy:
dj=1-ej d j =1-e j
计算各项指标的权值;Calculate the weight of each indicator;
计算得出权重,确定权重W:Calculate the weight and determine the weight W:
W=(w1,w2,…,wn)W=(w 1 , w 2 , . . . , wn )
式中权重和为1。where the weight sum is 1.
本实施例中 In this example
步骤5,计算综合评价指数,采用(+,*)算子模糊合成权重W与矩阵R得到最终评价向量S:Step 5: Calculate the comprehensive evaluation index, and use the (+, *) operator fuzzy synthesis weight W and matrix R to obtain the final evaluation vector S:
S=W·R=(s1,s2,...,sm)S=W·R=(s 1 , s 2 , ..., s m )
本实施例中, In this embodiment,
s1=0*0.30403643+0*0.53612817+0.073*0.15329247+0*0.00654293=0.01119035031s 1 =0*0.30403643+0*0.53612817+0.073*0.15329247+0*0.00654293=0.01119035031
s2=0.14*0.30403643+0.14*0.53612817+0.927*0.15329247+0*0.00654293=0.25972516369s 2 =0.14*0.30403643+0.14*0.53612817+0.927*0.15329247+0*0.00654293=0.25972516369
s3=0.86*0.30403643+0.86*0.53612817+0*0.15329247+0*0.00654293=0.722541556s 3 =0.86*0.30403643+0.86*0.53612817+0*0.15329247+0*0.00654293=0.722541556
s4=0*0.30403643+0*0.1378+0*0.53612817+0*0.15329247+0.47*0.00654293=0.0030751771s 4 =0*0.30403643+0*0.1378+0*0.53612817+0*0.15329247+0.47*0.00654293=0.0030751771
s5=0*0.30403643+0*0.1378+0*0.53612817+0*0.15329247+0.53*0.00654293=0.0034677529s 5 =0*0.30403643+0*0.1378+0*0.53612817+0*0.15329247+0.53*0.00654293=0.0034677529
S=(s1,s2,s3,s4,s5)S=(s 1 ,s 2 ,s 3 ,s 4 ,s 5 )
综合评价指数计算公式为:D=S*VThe calculation formula of the comprehensive evaluation index is: D=S*V
由表2计算D=S*VCalculate D=S*V from Table 2
D=S*V=s1*90+s2*70+s3*50+s4*30+s5*10=55.4419036282D=S*V=s 1 *90+s 2 *70+s 3 *50+s 4 *30+s 5 *10=55.4419036282
实施例3Example 3
利用实施例2的方法对中国西南油气田某地区的气井进行泡排效果评价,该井从2012 年4月18日开始进行泡沫排水采气工艺,截止到至2012年8月22日实施工艺天数是126天,故选择近等量的未泡排昨天天数进行对比,时间是2011年12月1日至2012年4月17 日。在现场工艺实施的过程中泡排剂的用量和周期是根据现场数据不断调整的。Utilize the method of
通过Python处理泡沫排水采气作业前后数据,用matplotlib工具库实现数据可视化,作出附图2泡沫排水采气作业前后产气量对比图、附图3泡沫排水采气作业前后产水量对比图以及附图3泡沫排水采气作业前后油套压差值对比图。The data before and after the foam drainage gas production operation is processed through Python, and the data visualization is realized by the matplotlib tool library, and the comparison chart of the gas production before and after the foam drainage gas production operation in Figure 2, and Figure 3 the comparison chart of the water production volume before and after the foam drainage gas production operation and the attached drawings are made 3. Comparison chart of oil jacket pressure difference before and after foam drainage gas production operation.
整合数据并结合本发明的计算公式计算得出基于熵值法的综合评价指数,整合油套压差变化率的数据,通过基于Python的matplotlib工具库实现数据可视化,如附图5所示,并分析综合评价指数与油套压差变化率的关系,从而更好的分析评价排水采气效果。Integrate the data and calculate the comprehensive evaluation index based on the entropy method in combination with the calculation formula of the present invention, integrate the data of the oil jacket pressure difference rate of change, and realize data visualization through the Python-based matplotlib tool library, as shown in Figure 5, and The relationship between the comprehensive evaluation index and the change rate of the oil jacket pressure difference is analyzed, so as to better analyze and evaluate the drainage gas recovery effect.
如图2-4所示,随着泡排剂的不断添加使得油套压差不断减小,产水量和产气量在2012 年5月16日开始逐渐增大,随后保持比较稳定的生产,油套压差在6月25日至7月10日减小至14.56和14.78MPa的最小值区间,井筒内的积液基本被排出,在此随后油套压差值稳定在一个区间内,地层中的产出水量相对比较稳定,表示此时泡沫排水采气生产作业制度较合理达到排出井底积液实现稳定生产的目的。As shown in Figure 2-4, with the continuous addition of the foam discharge agent, the pressure difference of the oil jacket continued to decrease, and the water production and gas production began to increase gradually on May 16, 2012, and then maintained a relatively stable production. From June 25 to July 10, the casing pressure difference decreased to a minimum value of 14.56 and 14.78MPa, and the accumulated fluid in the wellbore was basically discharged. After that, the oil casing pressure difference stabilized within an interval, and the formation The produced water volume is relatively stable, indicating that the foam drainage gas production operation system is more reasonable at this time to achieve the purpose of draining the bottom hole fluid and achieving stable production.
从气井生产效果来看,泡沫排水采气作业在2012年5月16日开始产气量和产水量开始增大,达到了施工作业目的,但是基于熵值法的泡沫排水采气作业综合评价指数只有40,属于效果一般的阶段,因为考虑到日排水采气作业成本的话,并未达到较好的经济效益。当油套压差变化率达到30%时,泡沫排水采气作业综合评价指数才达到60以上,施工作业效果达到好的阶段,具有相对较好的经济效益。Judging from the production effect of gas wells, the gas production and water production of foam drainage gas production began to increase on May 16, 2012, which achieved the purpose of the construction operation. However, the comprehensive evaluation index of foam drainage gas production based on the entropy method is only 40. It belongs to the stage with general effect, because considering the cost of daily drainage and gas production, it has not achieved good economic benefits. When the change rate of the oil jacket pressure difference reaches 30%, the comprehensive evaluation index of the foam drainage gas recovery operation reaches above 60, the construction operation effect reaches a good stage, and it has relatively good economic benefits.
从图5泡沫排水采气作业综合评价指数曲线与传统的生产参数曲线对比来看,本文提出的评价方法在结合排水采气作业成本的基础上对泡排效果评价参数更加全面,表现更加丰富,不同作业效果阶段表现更加明显,更加有利于泡沫排水采气的设计和起泡剂的优选。From the comparison between the comprehensive evaluation index curve of foam drainage and gas production operation and the traditional production parameter curve in Figure 5, the evaluation method proposed in this paper is more comprehensive and has richer performance on the evaluation parameters of foam drainage effect on the basis of the cost of drainage gas production operation. The performance of different operation effect stages is more obvious, which is more conducive to the design of foam drainage and gas recovery and the selection of foaming agent.
以上所述仅是本发明的优选实施方式,应当指出的是,对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以做出若干改进和润饰,这些改进和润饰也应视为本发明的保护范围。The above are only the preferred embodiments of the present invention. It should be noted that, for those skilled in the art, without departing from the principles of the present invention, several improvements and modifications can be made. These improvements and Retouching should also be regarded as the protection scope of the present invention.
Claims (10)
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010399587.3A CN111914209A (en) | 2020-05-12 | 2020-05-12 | Drainage gas production effect fuzzy comprehensive evaluation method based on entropy method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010399587.3A CN111914209A (en) | 2020-05-12 | 2020-05-12 | Drainage gas production effect fuzzy comprehensive evaluation method based on entropy method |
Publications (1)
Publication Number | Publication Date |
---|---|
CN111914209A true CN111914209A (en) | 2020-11-10 |
Family
ID=73237582
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010399587.3A Pending CN111914209A (en) | 2020-05-12 | 2020-05-12 | Drainage gas production effect fuzzy comprehensive evaluation method based on entropy method |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111914209A (en) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115859223A (en) * | 2023-02-27 | 2023-03-28 | 四川省计算机研究院 | Multi-source data industry fusion analysis method and system |
CN117391551A (en) * | 2023-12-12 | 2024-01-12 | 天津朔程科技有限公司 | Gas well drainage gas production evaluation method |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104794361A (en) * | 2015-05-05 | 2015-07-22 | 中国石油大学(华东) | Comprehensive evaluation method for water flooding oil reservoir development effect |
CN105469196A (en) * | 2015-11-18 | 2016-04-06 | 山东科技大学 | Comprehensive evaluation method and comprehensive evaluation system for evaluating mine construction project process |
CN106055871A (en) * | 2016-05-20 | 2016-10-26 | 中国石油天然气股份有限公司 | Evaluation method for gas recovery effect of gas well foam drainage |
CN109559017A (en) * | 2018-10-30 | 2019-04-02 | 中国石油天然气股份有限公司 | Technical index limit before foam drainage gas production implementation of gas well in Qinghai gas area and determination of technical index limit |
-
2020
- 2020-05-12 CN CN202010399587.3A patent/CN111914209A/en active Pending
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104794361A (en) * | 2015-05-05 | 2015-07-22 | 中国石油大学(华东) | Comprehensive evaluation method for water flooding oil reservoir development effect |
CN105469196A (en) * | 2015-11-18 | 2016-04-06 | 山东科技大学 | Comprehensive evaluation method and comprehensive evaluation system for evaluating mine construction project process |
CN106055871A (en) * | 2016-05-20 | 2016-10-26 | 中国石油天然气股份有限公司 | Evaluation method for gas recovery effect of gas well foam drainage |
CN109559017A (en) * | 2018-10-30 | 2019-04-02 | 中国石油天然气股份有限公司 | Technical index limit before foam drainage gas production implementation of gas well in Qinghai gas area and determination of technical index limit |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115859223A (en) * | 2023-02-27 | 2023-03-28 | 四川省计算机研究院 | Multi-source data industry fusion analysis method and system |
CN117391551A (en) * | 2023-12-12 | 2024-01-12 | 天津朔程科技有限公司 | Gas well drainage gas production evaluation method |
CN117391551B (en) * | 2023-12-12 | 2024-03-22 | 天津朔程科技有限公司 | Gas well drainage gas production evaluation method |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103729550B (en) | Multiple-model integration Flood Forecasting Method based on propagation time cluster analysis | |
CN106651656A (en) | Power supply reliability fuzzy comprehensive evaluation method based on improved membership degree function | |
CN109635435A (en) | A kind of natural river course stage discharge relation based on bayesian theory determines method | |
CN106021961A (en) | Urban standard cyclic working condition constructing method based on genetic algorithm optimization | |
CN111914209A (en) | Drainage gas production effect fuzzy comprehensive evaluation method based on entropy method | |
CN116658155A (en) | Shale gas well yield prediction method | |
CN109657327B (en) | The evaluation method of pump installation outlet passage comprehensive performance | |
CN111784536A (en) | A method for estimating water level overrun time based on measured water level changes in open channels | |
CN110130884A (en) | Method for determining residual potential of gas well | |
CN110084426A (en) | A kind of wind power interval prediction method based on industry and mining city | |
CN106223394B (en) | The mono- water reservoir system water resource optimal allocation method of single pumping plant-in library is directly mended under the conditions of insufficient irrigation | |
CN116776717A (en) | Multi-objective dynamic optimization method of drilling parameters based on improved NSGA-III algorithm | |
CN109635372A (en) | Design flood method based on the Bayesian model for improving prior probability | |
CN109559017A (en) | Technical index limit before foam drainage gas production implementation of gas well in Qinghai gas area and determination of technical index limit | |
CN111598429A (en) | Method for evaluating construction operation effect of foaming agent | |
CN113445988A (en) | Method for evaluating productivity of gas well of low-permeability carbonate rock gas reservoir | |
CN113139692A (en) | Waste slag dam deformation prediction method based on DPSO-ANFIS | |
CN111340293A (en) | A regional building energy consumption prediction method based on energy distribution pattern recognition | |
CN112253102A (en) | Method and device for determining air discharge pressure of oil well casing | |
CN117287192A (en) | Multi-layer oil and gas well yield prediction and reserve estimation method | |
CN117520931A (en) | Neural network-based optimization method for parameters of oil pumping well | |
CN107103412B (en) | A river water resource asset liability system and its calculation method | |
CN109919491A (en) | A joint calculation method of water and heat balance | |
CN110765660B (en) | Rapid determination method of movable reserves of low-permeability gas reservoirs under different gas production rates | |
CN115012925A (en) | Experimental determination method for vertical gas well shaft flow pattern under high pressure condition |
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 |