CN116629496A - 一种油藏窜流通道定量识别与分级方法 - Google Patents
一种油藏窜流通道定量识别与分级方法 Download PDFInfo
- Publication number
- CN116629496A CN116629496A CN202210123676.4A CN202210123676A CN116629496A CN 116629496 A CN116629496 A CN 116629496A CN 202210123676 A CN202210123676 A CN 202210123676A CN 116629496 A CN116629496 A CN 116629496A
- Authority
- CN
- China
- Prior art keywords
- evaluation
- longitudinal
- plane
- flow channel
- channel
- 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
- 238000000034 method Methods 0.000 title claims abstract description 50
- 230000005465 channeling Effects 0.000 title claims abstract description 48
- 238000011156 evaluation Methods 0.000 claims abstract description 193
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 70
- 238000004519 manufacturing process Methods 0.000 claims abstract description 31
- 238000011161 development Methods 0.000 claims abstract description 27
- 239000011159 matrix material Substances 0.000 claims description 36
- 238000002347 injection Methods 0.000 claims description 30
- 239000007924 injection Substances 0.000 claims description 30
- 238000010521 absorption reaction Methods 0.000 claims description 17
- 239000007788 liquid Substances 0.000 claims description 17
- 238000004088 simulation Methods 0.000 claims description 17
- 239000003129 oil well Substances 0.000 claims description 9
- 230000015572 biosynthetic process Effects 0.000 claims description 7
- 238000003786 synthesis reaction Methods 0.000 claims description 7
- 238000010276 construction Methods 0.000 claims description 6
- 230000002349 favourable effect Effects 0.000 abstract description 2
- 238000010586 diagram Methods 0.000 description 10
- 241000102542 Kara Species 0.000 description 6
- 230000035699 permeability Effects 0.000 description 4
- 230000008021 deposition Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000007547 defect Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000012407 engineering method Methods 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 239000000700 radioactive tracer Substances 0.000 description 1
- 238000011084 recovery Methods 0.000 description 1
- 239000000243 solution Substances 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
Classifications
-
- 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
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0631—Resource planning, allocation, distributing or scheduling for enterprises or organisations
- G06Q10/06311—Scheduling, planning or task assignment for a person or group
-
- 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
- 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
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0639—Performance analysis of employees; Performance analysis of enterprise or organisation operations
- G06Q10/06393—Score-carding, benchmarking or key performance indicator [KPI] analysis
-
- 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/06—Energy or water supply
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A10/00—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE at coastal zones; at river basins
- Y02A10/40—Controlling or monitoring, e.g. of flood or hurricane; Forecasting, e.g. risk assessment or mapping
Landscapes
- Business, Economics & Management (AREA)
- Engineering & Computer Science (AREA)
- Human Resources & Organizations (AREA)
- Physics & Mathematics (AREA)
- Economics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Strategic Management (AREA)
- Entrepreneurship & Innovation (AREA)
- Educational Administration (AREA)
- Marketing (AREA)
- Mathematical Physics (AREA)
- Development Economics (AREA)
- Tourism & Hospitality (AREA)
- General Business, Economics & Management (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Health & Medical Sciences (AREA)
- Data Mining & Analysis (AREA)
- Pure & Applied Mathematics (AREA)
- Quality & Reliability (AREA)
- Operations Research (AREA)
- Computational Mathematics (AREA)
- Game Theory and Decision Science (AREA)
- General Health & Medical Sciences (AREA)
- Computing Systems (AREA)
- Primary Health Care (AREA)
- Water Supply & Treatment (AREA)
- Algebra (AREA)
- Public Health (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种油藏窜流通道定量识别与分级方法，包括如下步骤：构建窜流通道平面评价图版；构建窜流通道纵向评价图版；基于窜流通道平面评价图版和纵向评价图版，利用模糊综合评判对窜流通道进行综合评价，得到窜流通道的分级结果。本发明先利用图版法对窜流通道在平面和纵向两个维度进行分级，然后在图版识别结果的基础上利用模糊综合评价算法进行综合评价与分级，得到窜流通道的分级结果。本发明适用于因窜流通道发育而导致高含水的油藏。本发明充分利用窜流通道生产特征，并且获取资料简单，识别和分级方法一体化，识别和分级结果方便有效，有利于现场的推广应用，为窜流通道的差异化处理提供基础。
Description
技术领域
本发明属于油气开发技术领域，具体涉及一种油藏窜流通道定量识别与分级方法。
背景技术
定量描述窜流通道是高含水老油田后期综合治理的关键技术之一，油井高含水往往是老油田重要的生产特征。油井高含水多是由于储层内部发育原生或次生的窜流通道，注入水沿着窜流通道向前快速突进，导致低效或无效注水循环，水驱波及系数和采收率低。为此需要明确地下储层窜流通道的分布位置和发育程度，并对发育程度进行分级，从而针对不同级别的窜流通道开展差异化的对策。
窜流通道的评价包括窜流通道的识别和窜流通道的分级。窜流通道的识别方法包括生产测井资料法、取芯井资料法、示踪剂监测资料法、试井资料法和油藏动静态资料法等方法。窜流通道的分级方法主要为油藏工程理论与数学算法结合的方法，其中包括模糊综合评判法。但是以上窜流通道的识别方法和分级方法都存在无法定量识别、资料获取困难、缺乏识别与分级方法一体化等缺点。
发明内容
本发明的目的在于提供一种油藏窜流通道定量识别与分级方法，以解决上述背景技术中提出的问题。
为实现上述目的，本发明提供如下技术方案：一种油藏窜流通道定量识别与分级方法，包括如下步骤：构建窜流通道平面评价图版；构建窜流通道纵向评价图版；基于窜流通道平面评价图版和纵向评价图版，利用模糊综合评判对窜流通道进行综合评价，得到窜流通道的分级结果。
进一步，所述构建窜流通道平面评价图版包括以下步骤：对窜流通道在平面评价图版进行识别；确定窜流通道在平面评价图版的分级界限；利用聚类算法对窜流通道的平面非均质情况分级，构建窜流通道平面评价图版，得到窜流通道在平面上的识别结果。
进一步，所述对窜流通道在平面评价图版进行识别包括：定义注水井方向流量比例和注水井方向油井在该层的归一化含水率指标，将所述指标作为窜流通道平面评价图版的横纵坐标值，实现窜流通道在平面评价图版上的识别。
进一步，所述确定窜流通道在平面评价图版的分级界限包括：建立反九点井网的数值模拟概念模型，设置不同窜流通道发育强度和发育规模，确定窜流通道在平面评价图版的分级界限。
进一步，所述利用聚类算法对窜流通道的平面非均质情况分级，构建窜流通道平面评价图版包括：以窜流通道平面评价初始图版上的直线为基础和约束，将实际数据点代入窜流通道平面评价初始图版，利用聚类算法将全区油水井间窜流通道的平面非均质情况进行分级，构建窜流通道平面评价图版，得到窜流通道在平面上的识别结果。
进一步，所述构建窜流通道纵向评价图版包括以下步骤：对窜流通道在纵向评价图版进行识别；确定窜流通道在纵向评价图版的分级界限；利用聚类算法对窜流通道的纵向非均质情况分级，构建窜流通道纵向评价图版，得到窜流通道在纵向上的识别结果。
进一步，所述对窜流通道在纵向评价图版进行识别包括：定义吸水或产液强度均质性系数和吸水或产液量非均质性系数指标，将所述指标作为窜流通道纵向评价图版的横纵坐标值，实现窜流通道在纵向评价图版上的识别。
进一步，所述确定窜流通道在纵向评价图版的分级界限包括：建立一注一采的数值模拟概念模型，设置不同的窜流通道发育强度和发育规模，确定窜流通道在纵向评价图版的分级界限。
进一步，所述利用聚类算法对窜流通道的纵向非均质情况分级，构建窜流通道纵向评价图版包括：以窜流通道纵向评价初始图版上的直线为基础和约束，将实际数据点代入窜流通道纵向评价初始图版，利用聚类算法将全区油水井间窜流通道的纵向非均质情况进行分级，构建窜流通道纵向评价图版，得到窜流通道在纵向上的识别结果。
进一步，所述得到窜流通道的分级结果包括以下步骤：确定因素集；确定评语集；确定因素集中各因素的权重；基于窜流通道平面评价图版和纵向评价图版确定隶属函数；基于隶属函数确定模糊综合评价矩阵；基于各因数权重和模糊综合评价矩阵综合判断，得到对窜流通道综合评价结果。
进一步，所述因素集为模糊综合评价指标，包括窜流通道平面评价结果和纵向评价结果。
进一步，所述评语集由窜流通道的平面和纵向非均质情况分级情况确定。
进一步，所述基于窜流通道平面评价图版和纵向评价图版确定隶属函数包括：将窜流通道平面评价图版和纵向评价图版中的每个区间的数据点进行聚类，得到聚类中心；通过各级非均质性聚类中心和对应的截距，得到模糊综合评价平面和纵向各级非均质性的评价标准；利用指派法确定隶属函数的模糊分布类型。
进一步，所述基于隶属函数确定模糊综合评价矩阵包括：
通过窜流通道平面和纵向评价图版的截距代入相应的非均质级别隶属函数得到隶属度，从而得到平面和纵向非均质性的模糊综合判断矩阵R，表达式如下所示：
其中R为模糊综合判断矩阵，r1n为平面非均质性对第n个评语的隶属度，r2n为纵向非均质性对第n个评语的隶属度，n为平面和纵向非均质性级别的数量。
进一步，所述基于各因数权重和模糊综合评价矩阵综合判断，得到对窜流通道综合评价结果，通过加入平面和纵向非均质性权重进行矩阵合成运算，得到综合评判矩阵B＝AR＝[b1,b2,b3,b4]，其中B为综合评判矩阵，A为指标权重，R为模糊综合判断矩阵，AR为矩阵合成运算，b1、b2、b3、b4为评价对象对各个评语的隶属度，取b1、b2、b3、b4中数值最大的评级作为综合评判结果。
本发明的技术效果和优点：该油藏窜流通道定量识别与分级方法先利用图版法对窜流通道在平面和纵向两个维度进行分级，然后在图版识别结果的基础上利用模糊综合评价算法进行综合评价与分级，得到窜流通道的分级结果。本发明适用于因窜流通道发育而导致高含水的油藏，包括长期水驱导致次生窜流通道发育油藏，因地质沉积、沉积微相、储层构型等原因导致的原生窜流通道发育油藏、裂缝性油藏等。本发明充分利用窜流通道生产特征，并且获取资料简单，识别和分级方法一体化，识别和分级结果方便有效，有利于现场的推广应用，为窜流通道的差异化处理提供基础。
附图说明
图1为本发明的油藏窜流通道定量识别与分级方法的流程示意图；
图2为本发明的反九点井网数值模拟概念模型示意图；
图3为本发明的基于反九点井网数值模拟概念模型确定的窜流通道平面评价初始图版示意图；
图4为本发明的一注一采的数值模拟概念模型示意图；
图5为本发明的基于一注一采数值模拟概念模型的窜流通道纵向评价初始图版示意图；
图6为本发明的卡拉姆卡斯油田J2C层窜流通道平面评价图版示意图；
图7为本发明的卡拉姆卡斯油田J2C层窜流通道纵向评价图版示意图；
图8为本发明的卡拉姆卡斯油田J2C层窜流通道平面非均质性各级隶属函数示意图；
图9为本发明的卡拉姆卡斯油田J2C层窜流通道纵向非均质性各级隶属函数示意图。
具体实施方式
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
本发明提供了如图1所示的一种油藏窜流通道定量识别与分级方法的流程示意图。
所述定量识别与分级方法包含如下步骤：构建窜流通道平面评价图版；构建窜流通道纵向评价图版；基于窜流通道平面评价图版和纵向评价图版，利用模糊综合评判对窜流通道进行综合评价，得到窜流通道的分级结果。
所述构建窜流通道平面评价图版包括以下步骤：对窜流通道在平面评价图版进行识别；确定窜流通道在平面评价图版的分级界限；利用聚类算法对窜流通道的平面非均质情况分级，构建窜流通道平面评价图版，得到窜流通道在平面上的识别结果。
所述对窜流通道在平面评价图版进行识别包括：定义注水井各方向流量比例和注水井方向油井在该层的归一化含水率指标，将所述指标作为窜流通道平面评价图版的横纵坐标值，实现窜流通道在平面评价图版上的识别。
在平面评价图版上窜流通道具有注水井流量分配比例大、油井饱和度低和含水率高的特点。注水井各方向流量比例和注水井各方向油井在该层的归一化含水率定义如下：
其中，fwmax为注水井各方向最大含水率，fw为注水井各方向含水率，fwmin为注水井各方向最小含水率，注水井各方向流量比例数值在0-1之间，数值越大，非均质性越高；注水井各方向油井含水率数值在0-1之间，数值越小，非均质性越高。
所述确定窜流通道在平面评价图版的分级界限包括：建立反九点井网的数值模拟概念模型，设置不同窜流通道发育强度和发育规模，确定窜流通道在平面评价图版的分级界限。
所述窜流通道发育强度用反九点井网的数值模拟概念模型中窜流通道和非窜流通道渗透率的级差表征。所述窜流通道的发育规模用窜流通道和非窜流通道的厚度比例表征，如表1所示。通过数值模拟模型得出结果，再计算相应的注水井各方向流量比例与各方向油井含水率，将结果数据画在平面评价图版上并连点成线得到窜流通道平面评价初始图版，如图3所示。
表1反九点井网概念模型参数
所述利用聚类算法对窜流通道的平面非均质情况分级，构建窜流通道平面评价图版包括：以窜流通道平面评价初始图版上的直线为基础和约束，将实际数据点代入窜流通道平面评价初始图版，利用聚类算法将全区油水井间窜流通道的平面非均质情况进行分级，构建窜流通道平面评价图版，得到窜流通道在平面上的识别结果。
在窜流通道平面评价初始图版中，以渗透率级差为1作为平面非均质性的界限，当数据点落在图3中级差为1直线的左上方时，不存在平面非均质性(非均质程度为四级)，当数据点落在图3中级差为1右下方时，注水井和生产井之间存在平面非均质性，且数据从左上角到右下角非均质性逐渐增强。
所述构建窜流通道纵向评价图版包括以下步骤：对窜流通道在纵向评价图版进行识别；确定窜流通道在纵向评价图版的分级界限；利用聚类算法对窜流通道的纵向非均质情况分级，构建窜流通道纵向评价图版，得到窜流通道在纵向上的识别结果。
所述对窜流通道在纵向评价图版进行识别包括：定义吸水或产液强度均质性系数和吸水或产液量非均质性系数指标，将所述指标作为窜流通道纵向评价图版的横纵坐标值，实现窜流通道在纵向评价图版上的识别。
在纵向评价图版上窜流通道具有吸水/产液强度高、吸水/产液量大的特点。吸水/产液强度均质系数和吸水/产液量非均质性系数定义如下：
其中，吸水/产液强度均质性系数在0-1之间，数值越大，非均质性越高；吸水/产液量非均质性系数在0-1之间，数值越大，非均质性越高。
所述确定窜流通道在纵向评价图版的分级界限包括：建立一注一采的数值模拟概念模型，设置不同的窜流通道发育强度和发育规模，确定窜流通道在纵向评价图版的分级界限。
所述窜流通道发育强度用一注一采的数值模拟概念模型上下两层渗透率的级差表征。所述窜流通道的发育规模用上下两层的厚度比例表征，如表2所示。通过数值模拟模型得出结果，再计算得出相应的吸水/产液强度均质性系数和吸水/产液量非均质系数，将结果数据画在纵向评价图版上并连点成线得到窜流通道纵向评价初始图版，如图5所示。
表2一注一采概念模型参数
所述利用聚类算法对窜流通道的纵向非均质情况分级，构建窜流通道纵向评价图版包括：以窜流通道纵向评价初始图版上的直线为基础和约束，将实际数据点代入窜流通道纵向评价初始图版，利用聚类算法将全区油水井间窜流通道的纵向非均质情况进行分级，构建窜流通道纵向评价图版，得到窜流通道在纵向上的识别结果。
在窜流通道纵向评价初始图版中，以渗透率级差为1作为纵向非均质性的界限，当数值落在图5中级差为1直线的左上方时，不存在纵向非均质性，当数值落在图5中级差为1直线的右下方时，注水井和生产井之间存在纵向非均质性，且数据从左上角到右下角非均质性逐渐增强。
所述得到窜流通道的分级结果包括以下步骤：确定因素集；确定评语集；确定因素集中各因素的权重；基于窜流通道平面评价图版和纵向评价图版确定隶属函数；基于隶属函数确定模糊综合评价矩阵；基于各因数权重和模糊综合评价矩阵综合判断，得到对窜流通道综合评价结果。
所述因素集为模糊综合评价指标，包括窜流通道平面评价结果和纵向评价结果，所述因素集表达式为U＝{平面非均质性，纵向非均质性}。
所述评语集由窜流通道的平面和纵向非均质情况分级情况确定。
所述基于窜流通道平面评价图版和纵向评价图版确定隶属函数包括：将窜流通道平面评价图版和纵向评价图版中的每个区间的数据点进行聚类，得到聚类中心；通过各级非均质性聚类中心和对应的截距，得到模糊综合评价平面和纵向各级非均质性的评价标准；利用指派法确定隶属函数的模糊分布类型。
常规模糊综合评价中，指标界限λi(i＝1,2,3…)为一维值，本发明中两个指标确定隶属度，相当于模糊综合评判方法在二维的拓展，界限值变成一条直线。在窜流通道平面评价图版和纵向评价图版中，各个分级界限直线斜率大致相同，这样通过求直线在纵轴上的截距将二维问题转换成一维问题。
所述基于隶属函数确定模糊综合评价矩阵包括：
通过窜流通道平面和纵向评价图版的截距代入相应的非均质级别隶属函数得到隶属度，从而得到平面和纵向非均质性的模糊综合判断矩阵R，表达式如下所示：
其中R为模糊综合判断矩阵，r1n为平面非均质性对第n个评语的隶属度，r2n为纵向非均质性对第n个评语的隶属度，n为平面和纵向非均质性级别的数量。
所述基于各因数权重和模糊综合评价矩阵综合判断，得到对窜流通道综合评价结果，通过加入平面和纵向非均质性权重进行矩阵合成运算，得到综合评判矩阵B＝AR＝[b1,b2,b3,b4]，其中B为综合评判矩阵，A为指标权重，R为模糊综合判断矩阵，AR为矩阵合成运算，b1、b2、b3、b4为评价对象对各个评语的隶属度，取b1、b2、b3、b4中数值最大的评级作为综合评判结果。
实施例
本实施例为卡拉姆卡斯油田J2C层油藏窜流通道定量识别与分级方法，包括如下步骤：
步骤1、构建卡拉姆卡斯油田J2C层窜流通道平面评价图版，包括以下步骤：
1、定义注水井各方向流量比例和注水井各方向油井在该层的归一化含水率两个指标，将所述两个指标作为窜流通道平面评价图版的横纵坐标值；
2、建立反九点井网的数值模拟概念模型，如图2所示，设置不同窜流通道发育强度和发育规模，得到平面评价初始图版，如图3所示；
3、以初始图版上的直线为基础和约束，将卡拉姆卡斯油田J2C层油水井实际数据点代入初始图版，利用聚类算法将全区油水井间窜流通道的平面非均质情况划分为4级，构建窜流通道平面评价图版，如图6所示，其中五角星为每个分级区域的聚类中心。
步骤2、构建卡拉姆卡斯油田J2C层窜流通道纵向评价图版，包括以下步骤：
1、定义吸水/产液强度均质性系数和吸水/产液量非均质性系数两个指标，将所述两个指标作为窜流通道纵向评价图版的横纵坐标值；
2、建立一注一采的数值模拟概念模型，如图4所示，设置不同的窜流通道发育强度和发育规模，如表2所示，得到纵向评价初始图版，如图5所示；
3、以初始图版上的直线为基础和约束，将卡拉姆卡斯油田J2C层油水井实际数据点代入初始图版，利用聚类算法将全区油水井间窜流通道的纵向非均质情况划分为4级，构建窜流通道纵向评价图版，如图7所示，其中五角星为每个分级区域的聚类中心。
步骤3、基于模糊综合评判对窜流通道综合评价，包括以下步骤：
1、因素集计为U＝{平面非均质性，纵向非均质性}。
2、评语集计为V＝{一级，二级，三级，四级}。
3、因素集中各因素的权重计为A＝{0.5,0.5}。
4、根据卡拉姆卡斯油田J2C层窜流通道平面和纵向评价图版求出窜流通道各级非均质性聚类中心坐标及截距，如表3所示，得到平面和纵向非均质性评价标准，如表4所示。
5、隶属函数采用梯形分布隶属函数，根据平面得到平面各级非均质性隶属函数，如图8所示，根据纵向非均质性评价标准得到纵向各级非均质性隶属函数，如图9所示。
6、以卡拉姆卡斯油田J2C层1028井为例，得出卡拉姆卡斯油田J2C层窜流通道平面图板上的坐标为(0.0345,0.0684)，代入平面图版分级界限直线方程，得到过点(0.0345,0.0684)且与分级界限直线平行的直线方程为f＝1.629q+0.0122，可知该直线在纵轴上的截距为0.0122；该井在纵向图版上的坐标为(0.5400,0.8072)，代入纵向图版分级界限直线方程，得到过点(0.5400,0.8072)且与分级界限直线平行的直线方程为H＝1.061p+0.2343，可知该直线在纵轴上的截距为0.2343。
7、代入相应隶属函数得到模糊综合评判矩阵为
8、加入平面和纵向非均质性权重进行矩阵合成运算，得到综合评判矩阵B＝AR＝[0,0.39185,0.10815,0.5]，取数值最大的评级(四级)作为综合评判结果，即卡拉姆卡斯油田J2C层1028井与注水井间发育四级非均质性窜流通道。
表3卡拉姆卡斯油田J2C层窜流通道平面和纵向图版各级非均质性聚类中心坐标及截距
表4卡拉姆卡斯油田J2C层窜流通道平面和纵向各级非均质评价标准
工作原理：该油藏窜流通道定量识别与分级方法基于图版法和模糊综合评判法，将反映窜流通道生产特征的油/水井产液/吸水剖面数据、注水井流量数据、含水率数据等数据充分利用，并利用油藏工程方法、油藏数值模拟方法和数学算法，构建窜流通道平面和纵向评价图版，进行窜流通道在平面和纵向上的识别，最后利用模糊综合评判算法将平面和纵向上的识别结果按照特定的权重综合，得到最终的窜流通道的分级结果。
最后应说明的是：以上所述仅为本发明的优选实施例而已，并不用于限制本发明，尽管参照前述实施例对本发明进行了详细的说明，对于本领域的技术人员来说，其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分技术特征进行等同替换,凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。
Claims (15)
1.一种油藏窜流通道定量识别与分级方法，其特征在于：所述定量识别与分级方法包括如下步骤：
构建窜流通道平面评价图版；
构建窜流通道纵向评价图版；
基于窜流通道平面评价图版和纵向评价图版，利用模糊综合评判对窜流通道进行综合评价，得到窜流通道的分级结果。
2.根据权利要求1所述的油藏窜流通道定量识别与分级方法，其特征在于：
所述构建窜流通道平面评价图版包括以下步骤：
对窜流通道在平面评价图版进行识别；
确定窜流通道在平面评价图版的分级界限；
利用聚类算法对窜流通道的平面非均质情况分级，构建窜流通道平面评价图版，得到窜流通道在平面上的识别结果。
3.根据权利要求2所述的油藏窜流通道定量识别与分级方法，其特征在于：
所述对窜流通道在平面评价图版进行识别包括：定义注水井方向流量比例和注水井方向油井在该层的归一化含水率指标，将所述指标作为窜流通道平面评价图版的横纵坐标值，实现窜流通道在平面评价图版上的识别。
4.根据权利要求2所述的油藏窜流通道定量识别与分级方法，其特征在于：
所述确定窜流通道在平面评价图版的分级界限包括：建立反九点井网的数值模拟概念模型，设置不同窜流通道发育强度和发育规模，确定窜流通道在平面评价图版的分级界限。
5.根据权利要求2所述的油藏窜流通道定量识别与分级方法，其特征在于：
所述利用聚类算法对窜流通道的平面非均质情况分级，构建窜流通道平面评价图版包括：以窜流通道平面评价初始图版上的直线为基础和约束，将实际数据点代入窜流通道平面评价初始图版，利用聚类算法将全区油水井间窜流通道的平面非均质情况进行分级，构建窜流通道平面评价图版，得到窜流通道在平面上的识别结果。
6.根据权利要求1所述的油藏窜流通道定量识别与分级方法，其特征在于：
所述构建窜流通道纵向评价图版包括以下步骤：
对窜流通道在纵向评价图版进行识别；
确定窜流通道在纵向评价图版的分级界限；
利用聚类算法对窜流通道的纵向非均质情况分级，构建窜流通道纵向评价图版，得到窜流通道在纵向上的识别结果。
7.根据权利要求6所述的油藏窜流通道定量识别与分级方法，其特征在于：
所述对窜流通道在纵向评价图版进行识别包括：定义吸水或产液强度均质性系数和吸水或产液量非均质性系数指标，将所述指标作为窜流通道纵向评价图版的横纵坐标值，实现窜流通道在纵向评价图版上的识别。
8.根据权利要求6所述的油藏窜流通道定量识别与分级方法，其特征在于：
所述确定窜流通道在纵向评价图版的分级界限包括：建立一注一采的数值模拟概念模型，设置不同的窜流通道发育强度和发育规模，确定窜流通道在纵向评价图版的分级界限。
9.根据权利要求6所述的油藏窜流通道定量识别与分级方法，其特征在于：
所述利用聚类算法对窜流通道的纵向非均质情况分级，构建窜流通道纵向评价图版包括：以窜流通道纵向评价初始图版上的直线为基础和约束，将实际数据点代入窜流通道纵向评价初始图版，利用聚类算法将全区油水井间窜流通道的纵向非均质情况进行分级，构建窜流通道纵向评价图版，得到窜流通道在纵向上的识别结果。
10.根据权利要求1所述的油藏窜流通道定量识别与分级方法，其特征在于：
所述得到窜流通道的分级结果包括以下步骤：
确定因素集；
确定评语集；
确定因素集中各因素的权重；
基于窜流通道平面评价图版和纵向评价图版确定隶属函数；
基于隶属函数确定模糊综合评价矩阵；
基于各因数权重和模糊综合评价矩阵综合判断，得到对窜流通道综合评价结果。
11.根据权利要求10所述的油藏窜流通道定量识别与分级方法，其特征在于：
所述因素集为模糊综合评价指标，包括窜流通道平面评价结果和纵向评价结果。
12.根据权利要求10或11所述的油藏窜流通道定量识别与分级方法，其特征在于：
所述评语集由窜流通道的平面和纵向非均质情况分级情况确定。
13.根据权利要求10所述的油藏窜流通道定量识别与分级方法，其特征在于：
所述基于窜流通道平面评价图版和纵向评价图版确定隶属函数包括：将窜流通道平面评价图版和纵向评价图版中的每个区间的数据点进行聚类，得到聚类中心；通过各级非均质性聚类中心和对应的截距，得到模糊综合评价平面和纵向各级非均质性的评价标准；利用指派法确定隶属函数的模糊分布类型。
14.根据权利要求10所述的油藏窜流通道定量识别与分级方法，其特征在于：
所述基于隶属函数确定模糊综合评价矩阵包括：
通过窜流通道平面和纵向评价图版的截距代入相应的非均质级别隶属函数得到隶属度，从而得到平面和纵向非均质性的模糊综合判断矩阵R，表达式如下所示：
其中R为模糊综合判断矩阵，r1n为平面非均质性对第n个评语的隶属度，r2n为纵向非均质性对第n个评语的隶属度，n为平面和纵向非均质性级别的数量。
15.根据权利要求10所述的油藏窜流通道定量识别与分级方法，其特征在于：
所述基于各因数权重和模糊综合评价矩阵综合判断，得到对窜流通道综合评价结果，通过加入平面和纵向非均质性权重进行矩阵合成运算，得到综合评判矩阵B＝AR＝[b1,b2,b3,b4]，其中B为综合评判矩阵，A为指标权重，R为模糊综合判断矩阵，AR为矩阵合成运算，b1、b2、b3、b4为评价对象对各个评语的隶属度，取b1、b2、b3、b4中数值最大的评级作为综合评判结果。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210123676.4A CN116629496A (zh) | 2022-02-10 | 2022-02-10 | 一种油藏窜流通道定量识别与分级方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210123676.4A CN116629496A (zh) | 2022-02-10 | 2022-02-10 | 一种油藏窜流通道定量识别与分级方法 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN116629496A true CN116629496A (zh) | 2023-08-22 |
Family
ID=87601334
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202210123676.4A Pending CN116629496A (zh) | 2022-02-10 | 2022-02-10 | 一种油藏窜流通道定量识别与分级方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN116629496A (zh) |
-
2022
- 2022-02-10 CN CN202210123676.4A patent/CN116629496A/zh active Pending
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108009716B (zh) | 一种水平井体积压裂效果影响因素多层次评价方法 | |
| CN110288258B (zh) | 一种高含水油藏剩余油挖潜方法 | |
| CN108446797B (zh) | 一种致密油藏水平井体积压裂初期产能预测方法 | |
| CN111425193B (zh) | 一种基于聚类分析测井岩石物理相划分的储层可压性评价方法 | |
| CN105654236B (zh) | 一种地下水型饮用水水源地污染风险评价方法 | |
| CN108960651A (zh) | 一种致密油气藏多级压裂水平井完井效率的综合评价方法 | |
| CN107291667B (zh) | 一种井间连通程度确定方法及系统 | |
| CN111173507B (zh) | 高含水油田剩余油预测方法 | |
| CN110821486B (zh) | 一种储层优势通道物性参数计算方法 | |
| CN110096815B (zh) | 一种流场表征的方法 | |
| CN113570226A (zh) | 断层破碎带隧道突涌水灾害发生概率等级评价方法 | |
| CN110533237A (zh) | 一种砂岩储层含油气产能预测方法 | |
| CN106569272A (zh) | 一种基于数据性质空间升维的地震属性融合方法 | |
| CN116384618A (zh) | 一种大井距油藏调驱选井方法和系统 | |
| CN114595628A (zh) | 一种水平井的体积压裂差异化改造方法 | |
| CN111832962B (zh) | 一种油田探明储量品质快速评价图版的建立方法 | |
| CN112560310A (zh) | 一种多层油藏渗透率动态变化定量表征的方法 | |
| CN116629496A (zh) | 一种油藏窜流通道定量识别与分级方法 | |
| CN111155980A (zh) | 一种水流优势通道识别方法及装置 | |
| CN115705452A (zh) | 整装砂岩油藏开发中后期的新型采收率预测方法 | |
| CN112983406B (zh) | 一种天然气水合物储层参数指标评价方法 | |
| CN114066166A (zh) | 一种沙砾岩低渗透油藏co2混相驱筛选方法 | |
| CN115144900A (zh) | 一种基于地震多属性聚类融合的碳酸盐岩储层厚度识别方法、系统、装置及存储介质 | |
| Zhang et al. | Check for | |
| CN115688345A (zh) | 非均相驱后剩余油定量预测方法 |
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 |