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 PDF

Info

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
Application number
CN202010399587.3A
Other languages
Chinese (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.)
Zhuojia Technology Tangshan Co ltd
North China University of Science and Technology
Original Assignee
Zhuojia Technology Tangshan Co ltd
North China University of Science and Technology
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 Zhuojia Technology Tangshan Co ltd, North China University of Science and Technology filed Critical Zhuojia Technology Tangshan Co ltd
Priority to CN202010399587.3A priority Critical patent/CN111914209A/en
Publication of CN111914209A publication Critical patent/CN111914209A/en
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/16Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design 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

本发明公开了一种基于熵值法的排水采气效果模糊综合评价方法,针对泡沫排水采气井的特性,选取评价泡沫排水采气效果的指标,遵循定性与定量相结合的原则,构建出泡沫排水采气效果评价指标体系,考虑气井排水采气效果评价指标之间的相互关系,采用线性分析方法建立隶属度矩阵,采用熵值法确定指标权重,结合模糊综合评价法计算出气井泡沫排水采气效果的综合评价指数,利用本发明的方法可得到更为科学、准确和客观的评价信息。

Figure 202010399587

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.

Figure 202010399587

Description

一种基于熵值法的排水采气效果模糊综合评价方法A fuzzy comprehensive evaluation method of drainage gas recovery effect based on entropy method

技术领域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;

然后,计算每一评价指标的评价区间相对于评判集的隶属度

Figure BDA0002488920730000041
进而可求出隶属度构成的矩阵,也就是模糊关系矩阵R;Then, calculate the membership degree of the evaluation interval of each evaluation index relative to the evaluation set
Figure BDA0002488920730000041
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求算得到的隶属度

Figure BDA0002488920730000021
转化为相对值,令
Figure BDA0002488920730000022
然后进行归一化处理,得到tij;First, calculate the membership degree obtained in step 3
Figure BDA0002488920730000021
Converted to relative values, let
Figure BDA0002488920730000022
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 step 2, m=5, the evaluation set V={v 1 , v 2 , v 3 , v 4 , v 5 }, and the evaluation level is very good, good, good, General sum and difference, v 1 =90, v 2 =70, v 3 =50, v 4 =30, v 5 =10.

在上述技术方案中,所述步骤3中,In the above technical solution, in the step 3,

Figure BDA0002488920730000023
Figure BDA0002488920730000023

其中,求出

Figure BDA0002488920730000024
构成的矩阵得到模糊关系矩阵R。Among them, find
Figure BDA0002488920730000024
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:

Figure BDA0002488920730000025
Figure BDA0002488920730000025

负向指标:Negative indicators:

Figure BDA0002488920730000031
Figure BDA0002488920730000031

则r′ij或r″ij为第i个影响因素的第j个评价区间的数值(i=1,2,…,n;j=1,2,…,m),归一化后的数据记为tijThen 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占该评价区间的比重

Figure BDA0002488920730000032
第j项评价区间的熵值
Figure BDA0002488920730000033
In the above technical solution, in the step 4, tij accounts for the proportion of the evaluation interval
Figure BDA0002488920730000032
The entropy value of the jth evaluation interval
Figure BDA0002488920730000033

在上述技术方案中,所述步骤4中,计算信息熵冗余度:dj=l-ej,然后计算各项指标的权值

Figure BDA0002488920730000034
计算得出权重,确定权重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
Figure BDA0002488920730000034
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*vmIn 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,遵循定性与定量相结合的原则,构建出泡沫排水采气效果评价指标体系:Step 2, following the principle of combining qualitative and quantitative, construct the evaluation index system of foam drainage gas recovery effect:

对应综合评价指数与评价级别的关系,设定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,计算每一评价指标的评价区间相对于评判集的隶属度

Figure BDA0002488920730000041
进而可求出隶属度构成的矩阵,也就是模糊关系矩阵R;3.2, Calculate the membership degree of the evaluation interval of each evaluation index relative to the evaluation set
Figure BDA0002488920730000041
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求算得到的隶属度

Figure BDA0002488920730000042
转化为相对值,令
Figure BDA0002488920730000043
然后进行归一化处理,得到tij;4.1, the membership degree calculated in step 3
Figure BDA0002488920730000042
Converted to relative values, let
Figure BDA0002488920730000043
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}.

Figure BDA0002488920730000051
Figure BDA0002488920730000051

Figure BDA0002488920730000052
Figure BDA0002488920730000052

Figure BDA0002488920730000053
Figure BDA0002488920730000053

Figure BDA0002488920730000054
Figure BDA0002488920730000054

步骤2,确定评判集V,为了对应综合评价指数与评价级别的关系,将综合评价指数D 设定为m个取值区间,vj为取值区间的中值,在本实施例中Step 2, determine the evaluation set V, in order to correspond to the relationship between the comprehensive evaluation index and the evaluation level, the comprehensive evaluation index D is set to m value intervals, v j is the median value of the value interval, in this embodiment

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

Figure BDA0002488920730000055
Figure BDA0002488920730000055

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

Figure BDA0002488920730000061
Figure BDA0002488920730000061

举例说明,在表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:

Figure BDA0002488920730000062
Figure BDA0002488920730000062

根据Δui与ai的关系和对应的关系式求解

Figure BDA0002488920730000063
的值,得出各评价指标相对于各评价区间的隶属度,进而可求出
Figure BDA0002488920730000064
构成的矩阵,也就是模糊关系矩阵R。According to the relationship between Δu i and a i and the corresponding relationship to solve
Figure BDA0002488920730000063
, 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.
Figure BDA0002488920730000064
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,再由隶属度的计算公式,计算得到

Figure BDA0002488920730000065
同理计算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
Figure BDA0002488920730000065
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:

Figure BDA0002488920730000071
Figure BDA0002488920730000071

步骤4,确定影响因素权重,选取n个影响因素(本实施例中选取了四个因素,即n=4), m个评价区间(本实施例中m=5),则

Figure BDA0002488920730000072
为第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
Figure BDA0002488920730000072
is the value of the j-th evaluation interval of the i-th influencing factor (i=1, 2,...,n; j=1, 2,...,m).

指标的归一化处理:异质指标同质化,由于各项指标的计量单位并不统一,首先要进行标准化处理,即把指标的绝对值转化为相对值,即令

Figure BDA0002488920730000073
从而解决各项不同质指标值的同质化问题。指标归一化式子如下: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
Figure BDA0002488920730000073
So as to solve the problem of homogeneity of different index values. The index normalization formula is as follows:

正向指标:Positive indicators:

Figure BDA0002488920730000074
Figure BDA0002488920730000074

负向指标:Negative indicators:

Figure BDA0002488920730000075
Figure BDA0002488920730000075

则r′ij或r″ij为第i个影响因素的第j个评价区间的数值(i=1,2,...,n;j=1,2,...,m),归一化后的数据记为tijThen 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:

Figure BDA0002488920730000076
Figure BDA0002488920730000076

其中i=1,...,n,j=1,...,m。where i=1,...,n,j=1,...,m.

计算第j项评价区间的熵值:Calculate the entropy value of the jth evaluation interval:

Figure BDA0002488920730000077
Figure BDA0002488920730000077

其中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;

Figure BDA0002488920730000081
Figure BDA0002488920730000081

计算得出权重,确定权重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.

本实施例中

Figure BDA0002488920730000082
In this example
Figure BDA0002488920730000082

步骤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 )

本实施例中,

Figure BDA0002488920730000083
In this embodiment,
Figure BDA0002488920730000083

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 embodiment 2 to carry out bubble discharge effect evaluation to a gas well in a certain area of China's southwest oil and gas field, this well starts to carry out foam drainage gas recovery process from April 18, 2012, and the number of days to implement the process as of August 22, 2012 is 126 days, so choose a nearly equal amount of unfoamed row yesterday for comparison, the time is from December 1, 2011 to April 17, 2012. During the implementation of the on-site process, the dosage and cycle of the foaming agent are continuously adjusted according to the on-site data.

通过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)

1. A drainage gas production effect fuzzy comprehensive evaluation method based on an entropy method is characterized by comprising the following steps:
step 1, selecting n evaluation indexes delta u influencing the drainage and gas production effectsi,i=1,2,…,n;
Step 2, corresponding to the relation between the comprehensive evaluation index D of the drainage gas production effect and the evaluation level, setting the comprehensive evaluation index D as m value intervals, vjThe median of the interval, j is 1, 2, …, m, from which the evaluation set V is determined1,v2,v3,…vm};
Step 3, establishing a membership matrix R by adopting a linear analysis method:
first, each evaluation index is divided into m evaluation sections, and a level limit value a corresponding to the evaluation section of each evaluation index is setj,j=1,2,…,m;
Then, the membership degree of the evaluation interval of each evaluation index relative to the evaluation set is calculated
Figure FDA0002488920720000013
Further, a matrix formed by membership degrees, namely a fuzzy relation matrix R can be solved;
step 4, determining the weight W of each evaluation index by adopting an entropy method:
first, the membership degree calculated in step 3 is calculated
Figure FDA0002488920720000012
Is converted into a relative value, such that
Figure FDA0002488920720000011
Then normalization processing is carried out to obtain tij
Then, the ith influence factor t under the jth evaluation interval is calculatedijThe proportion of the evaluation interval and the entropy value of the jth evaluation interval; further calculating the weight of each evaluation index, and then obtaining the weight W of each evaluation index;
and 5, carrying out fuzzy synthesis on each evaluation index weight W obtained in the step 4 and the fuzzy relation matrix R obtained in the step 3 by using (+,) operators to obtain a final evaluation vector S, and calculating a comprehensive evaluation index D of the gas recovery effect of the gas well foam drainage by using the final evaluation vector S.
2. The drainage gas production effect fuzzy comprehensive evaluation method based on the entropy method is characterized in that n is 4 in the step 1, and m is 5 in the step 2.
3. The entropy method-based drainage gas production effect fuzzy comprehensive evaluation method, according to claim 2, wherein the evaluation indexes are daily gas production change rates Δ qg. Rate of change of daily output water Δ qwThe oil jacket differential pressure change rate delta P and the daily foam drainage gas production operation cost delta C.
4. The drainage gas production effect fuzzy comprehensive evaluation method based on the entropy method as claimed in claim 2, wherein the evaluation set V ═ { V ═ V1,v2,v3,v4,v5The evaluation grade is good, general and poor; v. of1=90,v2=70,v3=50,v4=30,v5=10。
5. The drainage gas production effect fuzzy comprehensive evaluation method based on the entropy method as claimed in claim 1, wherein in the step 3,
Figure FDA0002488920720000021
wherein, find
Figure FDA0002488920720000022
The constructed matrix yields a fuzzy relation matrix R.
6. The drainage gas production effect fuzzy comprehensive evaluation method based on the entropy method as claimed in claim 1, wherein in the step 3, the normalization processing formula is as follows:
the forward direction index is as follows:
Figure FDA0002488920720000023
negative direction index:
Figure FDA0002488920720000024
then r'ijOr r ″)ijIs the ith influence factorThe value of the j-th evaluation interval of the element (i: 1, 2, …, n; j: 1, 2, …, m) is normalized and the normalized data is denoted as tij
7. The drainage gas production effect fuzzy comprehensive evaluation method based on the entropy method as claimed in claim 6, wherein in the step 4, t isijThe specific gravity of the evaluation interval
Figure FDA0002488920720000025
Entropy of j-th evaluation interval
Figure FDA0002488920720000026
8. The drainage gas production effect fuzzy comprehensive evaluation method based on the entropy method as claimed in claim 7, wherein in the step 4, the information entropy redundancy is calculated: dj=1-ejThen calculating the weight of each index
Figure FDA0002488920720000031
Calculating to obtain a weight, determining the weight W: w ═ W1,w2,…,wn) Wherein the sum of the weights is 1.
9. The method for fuzzy comprehensive evaluation of drainage gas production effect based on entropy method according to claim 8, wherein in the step 5, S-W-R (S)1,s2,...,sm)。
10. The entropy-method-based fuzzy comprehensive evaluation method for drainage gas production effects, according to claim 9, wherein in the step 5, D ═ S ═ V ═ S1*v1+s2*v2+…+sm*vm
CN202010399587.3A 2020-05-12 2020-05-12 Drainage gas production effect fuzzy comprehensive evaluation method based on entropy method Pending CN111914209A (en)

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)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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

Patent Citations (4)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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