CN112096348B - 一种天然气产层组非均质性评价方法 - Google Patents
一种天然气产层组非均质性评价方法 Download PDFInfo
- Publication number
- CN112096348B CN112096348B CN202011001164.8A CN202011001164A CN112096348B CN 112096348 B CN112096348 B CN 112096348B CN 202011001164 A CN202011001164 A CN 202011001164A CN 112096348 B CN112096348 B CN 112096348B
- Authority
- CN
- China
- Prior art keywords
- heterogeneity
- parameter
- reservoir
- evaluated
- group
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000011156 evaluation Methods 0.000 title claims abstract description 70
- 238000004519 manufacturing process Methods 0.000 title claims abstract description 64
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 title claims abstract description 48
- 239000003345 natural gas Substances 0.000 title claims abstract description 24
- 239000007789 gas Substances 0.000 claims abstract description 56
- 230000015572 biosynthetic process Effects 0.000 claims abstract description 9
- 239000003245 coal Substances 0.000 claims description 42
- 238000000034 method Methods 0.000 claims description 32
- 238000004458 analytical method Methods 0.000 claims description 10
- 239000005416 organic matter Substances 0.000 claims description 8
- 239000002734 clay mineral Substances 0.000 claims description 7
- 230000035699 permeability Effects 0.000 claims description 7
- 229910052500 inorganic mineral Inorganic materials 0.000 claims description 5
- 239000011707 mineral Substances 0.000 claims description 5
- 239000004576 sand Substances 0.000 claims 2
- 239000010410 layer Substances 0.000 description 62
- XQCFHQBGMWUEMY-ZPUQHVIOSA-N Nitrovin Chemical group C=1C=C([N+]([O-])=O)OC=1\C=C\C(=NNC(=N)N)\C=C\C1=CC=C([N+]([O-])=O)O1 XQCFHQBGMWUEMY-ZPUQHVIOSA-N 0.000 description 13
- 238000004364 calculation method Methods 0.000 description 9
- 230000006870 function Effects 0.000 description 9
- 230000008569 process Effects 0.000 description 7
- 238000011160 research Methods 0.000 description 6
- 239000004215 Carbon black (E152) Substances 0.000 description 5
- 229930195733 hydrocarbon Natural products 0.000 description 5
- 150000002430 hydrocarbons Chemical class 0.000 description 5
- 238000004587 chromatography analysis Methods 0.000 description 3
- 238000013507 mapping Methods 0.000 description 3
- 238000005065 mining Methods 0.000 description 3
- 239000010453 quartz Substances 0.000 description 3
- VYPSYNLAJGMNEJ-UHFFFAOYSA-N silicon dioxide Inorganic materials O=[Si]=O VYPSYNLAJGMNEJ-UHFFFAOYSA-N 0.000 description 3
- 238000010276 construction Methods 0.000 description 2
- 230000008021 deposition Effects 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 239000011435 rock Substances 0.000 description 2
- 238000003860 storage Methods 0.000 description 2
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 229910052799 carbon Inorganic materials 0.000 description 1
- 238000007405 data analysis Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 230000002349 favourable effect Effects 0.000 description 1
- 239000012530 fluid Substances 0.000 description 1
- 239000011229 interlayer Substances 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000035515 penetration Effects 0.000 description 1
- 230000000704 physical effect Effects 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 238000004445 quantitative analysis Methods 0.000 description 1
- 238000001179 sorption measurement Methods 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
Images
Classifications
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B43/00—Methods or apparatus for obtaining oil, gas, water, soluble or meltable materials or a slurry of minerals from wells
- E21B43/006—Production of coal-bed methane
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B47/00—Survey of boreholes or wells
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/02—Agriculture; Fishing; Forestry; Mining
Landscapes
- Engineering & Computer Science (AREA)
- Mining & Mineral Resources (AREA)
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- Geology (AREA)
- Business, Economics & Management (AREA)
- Geochemistry & Mineralogy (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Environmental & Geological Engineering (AREA)
- Fluid Mechanics (AREA)
- Agronomy & Crop Science (AREA)
- Economics (AREA)
- General Chemical & Material Sciences (AREA)
- Geophysics (AREA)
- Chemical & Material Sciences (AREA)
- Animal Husbandry (AREA)
- Marine Sciences & Fisheries (AREA)
- Chemical Kinetics & Catalysis (AREA)
- Health & Medical Sciences (AREA)
- Oil, Petroleum & Natural Gas (AREA)
- General Health & Medical Sciences (AREA)
- Human Resources & Organizations (AREA)
- Marketing (AREA)
- Primary Health Care (AREA)
- Strategic Management (AREA)
- Tourism & Hospitality (AREA)
- General Business, Economics & Management (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明提供的一种天然气产层组非均质性评价方法，应用于燃气开采技术领域，所述方法在确定待评价产层组所属的产层组类型之后，基于产层组类型表征的待评价产层组的储层构成，分别计算待评价产层组的各储层非均质性影响参数的参数值，以及各非均质性影响参数的权重系数；针对每一个非均质性影响参数，根据非均质性影响参数的参数值，计算用于表征非均质性影响参数与非均质评价结果的关联程度的隶属度，最终根据各储层非均质性影响参数的隶属度和权重系数，确定待评价产层组的非均质性评价结果，通过本方法，可以实现对包含两种及以上储层的天然气产层组进行非均质性评价，进而为煤系“三气”的整体开采提供参考依据。
Description
技术领域
本发明属于燃气开采技术领域，尤其涉及一种天然气产层组非均质性评价方法。
背景技术
在油气开采领域，非均质性主要指油气储层在经过漫长的沉积、成岩以及后期构造作用的综合影响后，在地层空间分布以及储层内部各种属性等方面所呈现出的差异。油气储层的非均质性越弱，开采难度越低，相反的，油气储层的非均质性越强，开采难度也就越高。因此，通过对油气储层进行非均质性评价，可以为油气开采提供参考依据，非均质性评价是油气开采领域中的重要工作。
以煤层气、页岩气和致密砂岩气共生叠置为特征的煤系“三气”，是近年来发现的一类重要的非常规天然气资源，根据现有气藏合层开采经验，把煤系天然气储层视作一个整体，综合勘探开发煤系煤层气、页岩气、致密砂岩气，不仅可以减少勘探开发成本，增大非常规天然气总储量和技术可采资源量，还可以提高气井使用效率和单井利润。
然而，现有技术中针对煤系天然气储层的非均质性评价方法，还大都以单一煤系气储层为评价对象，即只能分别对煤层气储层、页岩气储层、致密砂岩气储层的非均质性进行评价，对于包含两种及以上储层的煤系天然气产层组还缺乏行之有效的非均质性评价方法，难以对煤系“三气”的整体开采提供参考依据。
发明内容
有鉴于此，本发明的目的在于提供一种天然气产层组非均质性评价方法，对包含两种及以上储层的天然气产层组进行非均质性评价，为煤系“三气”的整体开采提供参考依据，具体方案如下：
本发明提供一种天然气产层组非均质性评价方法，包括：
确定待评价产层组所属的产层组类型；
基于所述产层组类型表征的所述待评价产层组的储层构成，分别计算所述待评价产层组的各非均质性影响参数的参数值，以及各所述非均质性影响参数的权重系数；
针对每一所述非均质性影响参数，根据所述非均质性影响参数的参数值，计算所述非均质性影响参数的隶属度，其中，所述隶属度用于表征所述非均质性影响参数与非均质评价结果的关联程度；
根据各所述非均质性影响参数的隶属度和权重系数，确定所述待评价产层组的非均质性评价结果。
可选的，所述根据各所述非均质性影响参数的隶属度和权重系数，确定所述待评价产层组的非均质性评价结果，包括：
针对每一所述非均质性影响参数，计算所述非均质性影响参数的隶属度与相应的权重系数的乘积，得到各所述非均质性影响参数的加权值；
计算各所述加权值之和，得到非均质性评价系数；
根据所述非均质性评价系数以及预设评价规则，确定所述待评价产层组的非均质性评价结果。
可选的，所述根据所述非均质性评价系数以及预设评价规则，确定所述待评价产层组的非均质性评价结果，包括：
若所述非均质性评价系数小于等于第一系数阈值，确定所述待评价产层组为强非均质性产层组；
若所述非均质性评价系数大于所述第一系数阈值，且小于第二系数阈值，确定所述待评价产层组为中等非均质性产层组；
若所述非均质性评价系数大于等于所述第二系数阈值，确定所述待评价产层组为弱非均质性产层组，其中，所述第一系数阈值小于所述第二系数阈值。
可选的，所述基于所述产层组类型表征的所述待评价产层组的各储层构成，分别计算所述待评价产层组的各非均质性影响参数的参数值，包括：
根据所述产层组类型表征的储层构成，确定所述待评价产层组包含的各储层；
分别获取各非均质性影响参数在各所述储层中的储层参数值和变异系数；
针对每一所述非均质性影响参数，根据所述非均质性影响参数的各所述储层参数值和各所述变异系数，计算所述非均质性影响参数在所述待评价产层组中的参数值。
可选的，所述根据所述非均质性影响参数的各所述储层参数值和各所述变异系数，计算所述非均质性影响参数在所述待评价产层组中的参数值，包括：
计算各所述变异系数之和，得到变异系数基值；
将所述非均质性影响参数的各所述储层参数值、各所述变异系数，以及所述变异系数基值代入下式，得到所述非均质性影响参数在所述待评价产层组中的参数值：
K＝KSHALE×C1/C+KCOAL×C2/C+KSAND×C3/C
上式中：K为所述非均质性影响参数在所述待评价产层组中的参数值；
KSHALE为所述非均质性影响参数在页岩储层中的参数值；
KCOAL为所述非均质性影响参数在煤层储层中的参数值；
KSAND为所述非均质性影响参数在致密砂岩储层中的参数值；
C为所述变异系数基值；
C1为所述非均质性影响参数在页岩储层中的变异系数；
C2为所述非均质性影响参数在煤层储层中的变异系数；
C3为所述非均质性影响参数在致密砂岩储层中的变异系数。
可选的，计算各所述非均质性影响参数的权重系数的过程，包括：
以非均质性评价结果为目标层、以各所述非均质性影响参数为对象层，以各所述非均质性影响参数与所述非均质性评价结果之间的关系作为准则，构建层次分析模型；
根据所述层次分析模型，计算各所述非均质性影响参数的权重系数。
可选的，所述根据所述非均质性影响参数的参数值，计算所述非均质性影响参数的隶属度，包括：
调用所述非均质性影响参数对应的隶属度函数；
将所述非均质性影响参数的参数值输入所述隶属度函数，得到所述非均质性影响参数的隶属度。
可选的，所述产层组类型至少包括：致密砂岩-页岩组合型、页岩-煤层组合型、致密砂岩-煤层组合型，以及致密砂岩-页岩-煤层组合型。
可选的，所述非均质性影响参数至少包括：厚度、含气量、渗透率、孔隙度、脆性矿物、黏土矿物、有机质含量、泊松比和弹性模量。
基于上述技术方案，本发明提供的天然气产层组非均质性评价方法，在确定待评价产层组所属的产层组类型之后，基于产层组类型表征的待评价产层组的储层构成，分别计算待评价产层组的各储层非均质性影响参数的参数值，以及各非均质性影响参数的权重系数；然后，针对每一个非均质性影响参数，根据非均质性影响参数的参数值，计算用于表征非均质性影响参数与非均质评价结果的关联程度的隶属度，最终根据各储层非均质性影响参数的隶属度和权重系数，确定待评价产层组的非均质性评价结果，通过本发明提供的产层组非均质性评价方法，可以实现对包含两种及以上储层的天然气产层组进行非均质性评价，进而为煤系“三气”的整体开采提供参考依据。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1是本发明实施例提供的一种天然气产层组非均质性评价方法流程图；
图2是本发明实施例提供的层析分析模型的结构框图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚，下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
参见图1，图1是本发明实施例提供的一种天然气产层组非均质性评价方法流程图，该方法可应用于电子设备，该电子设备可选如笔记本电脑、智能手机、PC(个人计算机)等具有数据分析能力的终端设备，显然，该电子设备在某些情况下也可选用网络侧的服务器实现；参照图1，本发明实施例提供的天然气产层组非均质性评价方法，可以包括：
S100，确定待评价产层组所属的产层组类型。
在实际应用中，气藏合层开采需要对研究区按照藏气压力系数划分为不同的含气系统，针对不同的含气系统采取不同的开采措施，比如，可以将压力系数比较接近1的地质区域划分为近正常压力系统，将压力系数小于1的地质区域划分为略微欠压系统。不同的含气系统之间相互独立，限制各含气系统间的流体交换及能量传递。需要说明的是，对于研究区含气系统的划分，可以按照现有技术中的方法实现，本发明对于含气系统的具体划分过程不做限定。
进一步的，还需要预先按照储层的不同，将属于同一含气系统的广大区域按照预设的划分规则划分为不同的产层组，当然，划分得到的各个产层组，相互之间也是相互独立的。
可选的，对于产层组的划分可以按照如下规则进行：
首先，考虑储层压裂时，产层不同岩性组合厚度及接触关系；其次，根据测井曲线考虑页岩层气测全烃含量，气测全烃含量较好的页岩储层将被划分进产层组；最后，考虑储层是符合单压还是合压条件，单压条件主要考虑煤层厚度大于10m、隔层厚度大于15m以及地层倾角大于45度，比如煤层厚度大于10米，相邻层只有煤层气存在，则大段煤层符合单压条件。
按照以上产层组划分原则，以致密砂岩层和煤层作为主要产气层，可以将同一含气系统内的广大区域内的产层组至少划分为四种类型：致密砂岩-页岩组合型、页岩-煤层组合型、致密砂岩-煤层组合型，以及致密砂岩-页岩-煤层组合型。比如，如果某一地质区域内，致密砂岩储层和页岩储层为主要的产气层，则该地质区域则可划分为致密砂岩-页岩组合型产层组，其他类型产层组的划分不再赘述。
在实际操作中，为便于对研究区内不同地质区域进行识别和标识，往往需要按照地理位置或其他相关条件，对研究区进行命名，比如，可以命名为山西组-太1段、山西组-太2段等，这种命名的方式可以参照现有技术或本领域中常用的方式实现，主要是为了简单、清晰的对各个产层组进行划分。
在实际应用中，本步骤确定待评价产层组所属的产层组类型，至少可以包括两种应用场景：其一，按照上述内容将研究区进行含气系统、产层组的划分，确定同一含气系统内，任一产层组所属的产层组类型。在这种情况下，在进行产层组的划分时，基本可以同步的确定待评价产层组所属的产层组类型。其二，如前所述，针对已经划分好的产层组，各个产层组是按照一定的命名规则予以命名的，本步骤中获取得到的是待评价产层组的标识信息，比如，产层组名称，或ID编码等，此种场景下，可以通过预设映射关系的方式，确定待评价产层组的产层组类型。当然，该预设映射关系中记录有待评价产层组标识信息与产层组类型的对应关系。比如，待评价产层组为山西组-太1段，通过查询预设映射关系，则可以确定山西组-太1段的产层组类型为致密砂岩-页岩组合型。
S110、基于产层组类型表征的待评价产层组的储层构成，分别计算待评价产层组的各非均质性影响参数的参数值，以及各非均质性影响参数的权重系数。
影响产层组非均质性的因素有很多，每一种因素对应一个非均质性影响参数，因此，在进行量化分析时，需要结合多个非均质性影响参数对待评价产层组的非均质性进行评价。
可选的，本发明实施例所提供的评价方法中所涉及的非均质性影响参数至少包括：厚度、含气量、渗透率、孔隙度、脆性矿物、黏土矿物、TOC、泊松比和弹性模量。下面对各个非均质性影响参数进行简要的介绍：
厚度是影响产层组非均质性的重要因素，产层组厚度的大小影响产层组的规模，产层组厚度越大，产层组规模越大，产层组采出气的可能性越大，是评价产层组非均质性的最基本的参数；
含气量反映产层组内各储层的含气性，是储层开采的物质基础，直接影响产层组开采的产量，具有经济价值的产层组必然具有一定的含气量；
渗透率是产层组内储层渗透性能的表征，反映储层的排采条件，是影响储层产气量高低的关键参数；
孔隙度提供储集空间,是影响储气量多少的重要参数；
在储层成分构成参数中，主要包括脆性矿物、黏土矿物和有机质含量。其中，有机质含量是评价有机质丰度的重要指标，有机质含量越高，生烃潜力就越大，越有利于烃类气体的吸附，而有机碳丰度是控制煤系烃源岩生烃能力的重要因素；
储层物质组成中的脆性矿物，比如石英，的含量影响产层组的可改造性；
黏土矿物的含量则可以体现沉积环境；
泊松比和弹性模量属于力学性能参数，力学性质的各向异性是评价储层非均质性的一个重要方面，力学参数作为保证油气井梳理开展勘探开发工作的基础数据，它反映了岩石受外力作用的物理性质，有利于解释地应力解释和评价井壁的稳定性。
如前所述，在现有应用中，更多的是对单一储层的分析，对于上述非均质性影响参数的计算也大都基于单一储层进行，还没有如何对产层组的上述各个非均质性影响参数进行计算的方法，因此，本发明实施例提供一种计算产层组各个非均质性影响参数的参数值的方法。
可选的，在采用上述任一种方法确定待评价产层组所属的产层组类型后，即可根据产层组类型所表征的储层构成，明确确定待评价产层组具体包含的各个储层。比如，如果待评价产层组为页岩-煤层组合型，即表明待评价产层组是由页岩储层和煤层储层构成，也就是说，产层组类型用于表征产层组的储层构成，在确定待评价产层组所属的产层组类型后，即可直接确定待评价产层组所包含的具体储层。
可以想到的是，待评价产层组的每一种储层，均可以用上述各非均质性影响参数进行分析，换句话说，针对每一个非均质性影响参数而言，在待评价产层组的不同储层中，都对应一个参数值，在本方案中定义为储层参数值。为计算待评价产层组的各个非均质性影响参数的参数值，首先需要获取各非均质性影响参数在各储层中的储层参数值和变异系数。
需要说明的是，针对单一储层中，各个非均质性影响参数的储层参数值以及变异系数的计算均可以按照现有技术实现，本发明对于非均质性影响参数的在单一储层中的储层参数值以及相对应的变异系数的计算方法不做限定。
在得到上述各个非均质性影响参数在各个储层中的储层参数值，以及变异系数之后，即可分别计算各非均质性影响参数在待评价产层组中的参数值。由于各个非均质性影响参数的参数值的计算过程一致，此处以任一非均质性影响参数为例，对计算过程进行说明。
首先，需要计算各非均质性影响参数对应的全部变异系数之和，并将所得之和作为变异系数基值。
然后，针对每一个非均质性影响参数，将该非均质性影响参数的各储层参数值、各变异系数，以及前述计算得到的变异系数基值代入下式，得到非均质性影响参数在待评价产层组中的参数值：
K＝KSHALE×C1/C+KCOAL×C2/C+KSAND×C3/C
上式中：K为非均质性影响参数在待评价产层组中的参数值；
KSHALE为非均质性影响参数在页岩储层中的参数值；
KCOAL为非均质性影响参数在煤层储层中的参数值；
KSAND为非均质性影响参数在致密砂岩储层中的参数值；
C为变异系数基值；
C1为非均质性影响参数在页岩储层中的变异系数；
C2为非均质性影响参数在煤层储层中的变异系数；
C3为非均质性影响参数在致密砂岩储层中的变异系数。
需要说明的是，采用上述公式进行计算时，对于未包含的储层，其对应的计算项取0即可。比如，待评价产层组属于致密砂岩-页岩组合型，并不包含煤层储层，在采用详述公式计算时，KCOAL和C2取0即可。
可以想到的是，在实际应用中，上述各个非均质性影响参数对最终的非均质性评价结果的影响程度必然是不同的，因此，本发明实施例提供的评价方法引入权重系数，为每一个非均质性影响参数设置一个权重系数，通过权重系数表征各非均质性影响参数对最终评价结果的影响程度。
可选的，在本发明实施例中，各个非均质性影响参数的权重系数通过层次分析法确定得到。以非均质性评价结果为目标层、以各非均质性影响参数为对象层，以各非均质性影响参数与非均质性评价结果之间的关系作为准则，分为含气参数和产气参数，构建层次分析模型，进而根据层次分析模型，计算各非均质性影响参数的权重系数。作为一种可选的层次分析模型，参见图2，图2是本发明实施例提供的层析分析模型的结构框图。需要说明的是，对于构建层析分析模型的具体过程，可以参照现有技术中层次分析理论实现，本发明对于层次分析模型的具体构建过程不做限定。
可选的，参见表1，表1示出了一种应用场景中，各个非均质性影响参数的权重系数的计算结果。
S120、针对每一非均质性影响参数，根据非均质性影响参数的参数值，计算非均质性影响参数的隶属度。
在本发明实施例中，隶属度用于表征非均质性影响参数与非均质评价结果的关联程度。为量化计算各非均质性影响参数的隶属度，本发明实施例为每一非均质性影响参数设置一对应的隶属度函数。
具体的，
①含气量
②渗透率
③厚度
④孔隙度
⑤有机质含量
⑥石英含量
⑦黏土矿物
⑧弹性模量
⑨泊松比
在上述各式中，x表示相应的非均质性影响参数的参数值。
在经过前述步骤，计算得到各非均质性影响参数的参数值之后，调用上述各非均质性影响参数对应的隶属度函数，将各参数的参数值输入相应的隶属度函数中，即可计算得到各非均质性影响参数所对应的隶属度。
可选的，本发明实施例中的隶属度函数是基于模糊数学理论得到。首先，建立评价对象的因素集合U，在本实施例中，因素集合U＝{含气量，渗透率，厚度，孔隙度，有机质含量，石英，黏土矿物，泊松比，弹性模量}。
然后，建立评价对象的结语集合V，结语集合V是与因素集合U中评价因素相应的评价标准集合，具体到本实施例中，结语集合V＝{强非均质性，中等非均质性，弱非均质性}。
最后，建立模糊矩阵，即在得到各非均质性影响参数的参数值的情况下，从一个单因素指标出发建立隶属度函数，以确定因素集合中某因素对结语集合中各元素的隶属度，进而得到各非均质性影响因素的隶属度函数。
需要说明的是，在应用模糊数学进行隶属度函数的构造时，实现细节可以参照现有技术中模糊数学的相关方法实现，本发明对此不做具体限定。
S130、根据各非均质性影响参数的隶属度和权重系数，确定待评价产层组的非均质性评价结果。
经过前述步骤，已经得到各非均质性影响参数的隶属度和权重系数，基于此，针对每一非均质性影响参数，计算该非均质性影响参数的隶属度与相应的权重系数的乘积，可以得到各非均质性影响参数对应的加权值。
然后，计算各加权值之和，便可得到待评价产层组的非均质性评价系数。最终根据非均质性评价系数以及预设评价规则，确定待评价产层组的非均质性评价结果。
可选的，可以根据非均质性评价系数的取值范围以及前述结语集合中的非均质性评价结论，设置不同的阈值，结合待评价产层组的非均质性评价系数以及预设的阈值，确定最终的非均质性评价结果。
具体的，预设第一系数阈值和第二系数阈值，且第一系数阈值小于第二系数阈值。如果待评价产层组的非均质性评价系数小于等于第一系数阈值，确定待评价产层组为强非均质性产层组；如果待评价产层组的非均质性评价系数大于第一系数阈值，且小于第二系数阈值，确定待评价产层组为中等非均质性产层组；如果待评价产层组的非均质性评价系数大于等于第二系数阈值，确定待评价产层组为弱非均质性产层组。
综上所述，本发明提供的天然气产层组非均质性评价方法，在确定待评价产层组所属的产层组类型之后，基于产层组类型表征的待评价产层组的储层构成，分别计算待评价产层组的各非均质性影响参数的参数值，以及各非均质性影响参数的权重系数；然后，针对每一个非均质性影响参数，根据非均质性影响参数的参数值，计算用于表征非均质性影响参数与非均质评价结果的关联程度的隶属度，最终根据各非均质性影响参数的隶属度和权重系数，确定待评价产层组的非均质性评价结果，通过本发明提供的产层组非均质性评价方法，可以实现对包含两种及以上储层的天然气产层组进行非均质性评价，进而为煤系“三气”的整体开采提供参考依据。
本说明书中各个实施例采用递进的方式描述，每个实施例重点说明的都是与其他实施例的不同之处，各个实施例之间相同相似部分互相参见即可。对于实施例公开的装置而言，由于其与实施例公开的方法相对应，所以描述的比较简单，相关之处参见方法部分说明即可。
专业人员还可以进一步意识到，结合本文中所公开的实施例描述的各示例的单元及算法步骤，能够以电子硬件、计算机软件或者二者的结合来实现，为了清楚地说明硬件和软件的可互换性，在上述说明中已经按照功能一般性地描述了各示例的组成及步骤。这些功能究竟以硬件还是软件方式来执行，取决于技术方案的特定应用和设计约束条件。专业技术人员可以对每个特定的应用来使用不同方法来实现所描述的功能，但是这种实现不应认为超出本发明的范围。
结合本文中所公开的实施例描述的方法或算法的步骤可以直接用硬件、处理器执行的软件模块，或者二者的结合来实施。软件模块可以置于随机存储器(RAM)、内存、只读存储器(ROM)、电可编程ROM、电可擦除可编程ROM、寄存器、硬盘、可移动磁盘、CD-ROM、或技术领域内所公知的任意其它形式的存储介质中。
对所公开的实施例的上述说明，使本领域专业技术人员能够实现或使用本发明。对这些实施例的多种修改对本领域的专业技术人员来说将是显而易见的，本文中所定义的一般原理可以在不脱离本发明的核心思想或范围的情况下，在其它实施例中实现。因此，本发明将不会被限制于本文所示的这些实施例，而是要符合与本文所公开的原理和新颖特点相一致的最宽的范围。
Claims (7)
1.一种天然气产层组非均质性评价方法，其特征在于，包括：
确定待评价产层组所属的产层组类型；
基于所述产层组类型表征的所述待评价产层组的储层构成，分别计算所述待评价产层组的各非均质性影响参数的参数值，以及各所述非均质性影响参数的权重系数；
针对每一所述非均质性影响参数，根据所述非均质性影响参数的参数值，计算所述非均质性影响参数的隶属度，其中，所述隶属度用于表征所述非均质性影响参数与非均质评价结果的关联程度；
根据各所述非均质性影响参数的隶属度和权重系数，确定所述待评价产层组的非均质性评价结果；
其中，所述基于所述产层组类型表征的所述待评价产层组的储层构成，分别计算所述待评价产层组的各非均质性影响参数的参数值，包括：
根据所述产层组类型表征的储层构成，确定所述待评价产层组包含的各储层；
分别获取各非均质性影响参数在各所述储层中的储层参数值和变异系数；
针对每一所述非均质性影响参数，根据所述非均质性影响参数的各所述储层参数值和各所述变异系数，计算所述非均质性影响参数在所述待评价产层组中的参数值；
其中，所述根据所述非均质性影响参数的各所述储层参数值和各所述变异系数，计算所述非均质性影响参数在所述待评价产层组中的参数值，包括：
计算各所述变异系数之和，得到变异系数基值；
将所述非均质性影响参数的各所述储层参数值、各所述变异系数，以及所述变异系数基值代入下式，得到所述非均质性影响参数在所述待评价产层组中的参数值：
K＝KSHALE×C1/C+KCOAL×C2/C+KSAND×C3/C
上式中：K为所述非均质性影响参数在所述待评价产层组中的参数值；
KSHALE为所述非均质性影响参数在页岩储层中的参数值；
KCOAL为所述非均质性影响参数在煤层储层中的参数值；
KSAND为所述非均质性影响参数在致密砂岩储层中的参数值；
C为所述变异系数基值；
C1为所述非均质性影响参数在页岩储层中的变异系数；
C2为所述非均质性影响参数在煤层储层中的变异系数；
C3为所述非均质性影响参数在致密砂岩储层中的变异系数。
2.根据权利要求1所述的天然气产层组非均质性评价方法，其特征在于，所述根据各所述非均质性影响参数的隶属度和权重系数，确定所述待评价产层组的非均质性评价结果，包括：
针对每一所述非均质性影响参数，计算所述非均质性影响参数的隶属度与相应的权重系数的乘积，得到各所述非均质性影响参数的加权值；
计算各所述加权值之和，得到非均质性评价系数；
根据所述非均质性评价系数以及预设评价规则，确定所述待评价产层组的非均质性评价结果。
3.根据权利要求2所述的天然气产层组非均质性评价方法，其特征在于，所述根据所述非均质性评价系数以及预设评价规则，确定所述待评价产层组的非均质性评价结果，包括：
若所述非均质性评价系数小于等于第一系数阈值，确定所述待评价产层组为强非均质性产层组；
若所述非均质性评价系数大于所述第一系数阈值，且小于第二系数阈值，确定所述待评价产层组为中等非均质性产层组；
若所述非均质性评价系数大于等于所述第二系数阈值，确定所述待评价产层组为弱非均质性产层组，其中，所述第一系数阈值小于所述第二系数阈值。
4.根据权利要求1所述的天然气产层组非均质性评价方法，其特征在于，计算各所述非均质性影响参数的权重系数的过程，包括：
以非均质性评价结果为目标层、以各所述非均质性影响参数为对象层，以各所述非均质性影响参数与所述非均质性评价结果之间的关系作为准则，构建层次分析模型；
根据所述层次分析模型，计算各所述非均质性影响参数的权重系数。
5.根据权利要求1所述的天然气产层组非均质性评价方法，其特征在于，所述根据所述非均质性影响参数的参数值，计算所述非均质性影响参数的隶属度，包括：
调用所述非均质性影响参数对应的隶属度函数；
将所述非均质性影响参数的参数值输入所述隶属度函数，得到所述非均质性影响参数的隶属度。
6.根据权利要求1-5任一项所述的天然气产层组非均质性评价方法，其特征在于，所述产层组类型至少包括：致密砂岩-页岩组合型、页岩-煤层组合型、致密砂岩-煤层组合型，以及致密砂岩-页岩-煤层组合型。
7.根据权利要求1-5任一项所述的天然气产层组非均质性评价方法，其特征在于，所述非均质性影响参数至少包括：厚度、含气量、渗透率、孔隙度、脆性矿物、黏土矿物、有机质含量、泊松比和弹性模量。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011001164.8A CN112096348B (zh) | 2020-09-22 | 2020-09-22 | 一种天然气产层组非均质性评价方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011001164.8A CN112096348B (zh) | 2020-09-22 | 2020-09-22 | 一种天然气产层组非均质性评价方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112096348A CN112096348A (zh) | 2020-12-18 |
| CN112096348B true CN112096348B (zh) | 2023-01-03 |
Family
ID=73755951
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202011001164.8A Active CN112096348B (zh) | 2020-09-22 | 2020-09-22 | 一种天然气产层组非均质性评价方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112096348B (zh) |
Family Cites Families (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US10352099B2 (en) * | 2015-09-02 | 2019-07-16 | Exxonmobil Upstream Research Company | Methods for drilling a wellbore within a subsurface region and drilling assemblies that include and/or utilize the methods |
| CN106651656A (zh) * | 2016-12-28 | 2017-05-10 | 国网江西省电力公司经济技术研究院 | 基于改进的隶属度函数的供电可靠性模糊综合评价方法 |
| CN106894800B (zh) * | 2017-02-23 | 2019-05-21 | 中国海洋石油集团有限公司 | 一种适用于海上稠油油藏的调剖选井决策方法 |
| CN109492860A (zh) * | 2018-09-26 | 2019-03-19 | 中国石油天然气股份有限公司 | 一种基于气井产能的致密砂岩储层多参数定量评价方法 |
| CN109184677A (zh) * | 2018-09-26 | 2019-01-11 | 西南石油大学 | 用于非均质互层砂体的储层评价方法 |
-
2020
- 2020-09-22 CN CN202011001164.8A patent/CN112096348B/zh active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN112096348A (zh) | 2020-12-18 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| KR102170765B1 (ko) | 딥러닝을 이용한 셰일가스 생산량 예측모델 생성 방법 | |
| CN105134191B (zh) | 致密油油井储量的评价方法 | |
| CN103775057A (zh) | 一种致密油气藏有效储层的判识方法及装置 | |
| Park et al. | Improved decision making with new efficient workflows for well placement optimization | |
| WO2022186182A1 (ja) | 予測装置、予測方法、及び、記録媒体 | |
| CN113361111B (zh) | 一种低渗透油藏储层特征模型的划分方法 | |
| Kalantari-Dahaghi et al. | New Insight into Integrated Reservoir Management using Top-Down, Intelligent Reservoir Modeling Technique; Application to a Giant and Complex Oil Field in the Middle East | |
| CN111101914B (zh) | 水平井压裂段簇优选方法及设备 | |
| CN112096348B (zh) | 一种天然气产层组非均质性评价方法 | |
| CN113792499A (zh) | 基于沉积特征的松散承压含水层富水性动态确定方法 | |
| Nnamdi | Conceptual Reservoir Development and Short-Term Forecasting Using Material Balance Based Integrated Asset Models and Neural Network Proxy Models | |
| Saputra | Physics-guided data-driven production forecasting in shales | |
| CN111691878A (zh) | 用于页岩的岩相品质分级评价方法及装置 | |
| CN111287739B (zh) | 一种基于地层原油粘度的剩余油分布预测方法 | |
| CN110400093B (zh) | 页岩油工业价值等级评价方法、装置及电子设备 | |
| CN107060744A (zh) | 一种测井地质评价系统 | |
| Zeng et al. | Optimal selection of stimulation wells using a fuzzy multicriteria methodology | |
| CN111898267A (zh) | 一种基于沉积地质条件的冲击地压危险性评价方法 | |
| Aslam et al. | Capacitance Resistance Clustered Model for Mature Peripheral Waterflood Performance Prediction & Optimization | |
| Al-Sane et al. | Intelligent Well Drainage and Interference Density Analysis Tool for Optimizing Infill Drilling | |
| Aslam et al. | An Integrated Multi-Disciplinary Modeling and Simulation Approach to History Matching and Forecasting Shale Wells-An Eagle Ford Case Study | |
| US20240093592A1 (en) | Quantification of pore-filling dolomite and calcite cement in carbonate reservoirs in post hydrocarbon charge stage | |
| CN105447291A (zh) | 一种获取可采资源量分布模型的方法 | |
| CN115295092A (zh) | 一种砾岩油藏储层分类方法 | |
| Lu et al. | A new production forecasting method of the multifractured horizontal wells based on cluster analysis |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |