CN106651060B - 一种定量评价断控天然气藏富集程度的方法 - Google Patents
一种定量评价断控天然气藏富集程度的方法 Download PDFInfo
- Publication number
- CN106651060B CN106651060B CN201710039408.3A CN201710039408A CN106651060B CN 106651060 B CN106651060 B CN 106651060B CN 201710039408 A CN201710039408 A CN 201710039408A CN 106651060 B CN106651060 B CN 106651060B
- Authority
- CN
- China
- Prior art keywords
- natural gas
- enrichment degree
- pool
- class
- factor
- 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
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 title claims abstract description 178
- 239000003345 natural gas Substances 0.000 title claims abstract description 89
- 238000000034 method Methods 0.000 title claims abstract description 43
- 239000007789 gas Substances 0.000 claims abstract description 41
- 238000004321 preservation Methods 0.000 claims abstract description 35
- 230000006835 compression Effects 0.000 claims abstract description 32
- 238000007906 compression Methods 0.000 claims abstract description 32
- 238000013508 migration Methods 0.000 claims abstract description 25
- 230000005012 migration Effects 0.000 claims abstract description 25
- 238000003860 storage Methods 0.000 claims abstract description 24
- 238000012360 testing method Methods 0.000 claims abstract description 14
- 238000012417 linear regression Methods 0.000 claims abstract description 9
- 239000011435 rock Substances 0.000 claims description 30
- 238000003325 tomography Methods 0.000 claims description 21
- 239000004215 Carbon black (E152) Substances 0.000 claims description 20
- 229930195733 hydrocarbon Natural products 0.000 claims description 20
- 150000002430 hydrocarbons Chemical group 0.000 claims description 20
- 238000002347 injection Methods 0.000 claims description 14
- 239000007924 injection Substances 0.000 claims description 14
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 14
- 238000004458 analytical method Methods 0.000 claims description 6
- 238000012935 Averaging Methods 0.000 claims description 4
- 230000000694 effects Effects 0.000 description 10
- 238000006073 displacement reaction Methods 0.000 description 5
- 239000012530 fluid Substances 0.000 description 5
- 230000015572 biosynthetic process Effects 0.000 description 4
- 239000000470 constituent Substances 0.000 description 4
- 230000001276 controlling effect Effects 0.000 description 4
- 238000005520 cutting process Methods 0.000 description 4
- 238000004451 qualitative analysis Methods 0.000 description 4
- 238000004445 quantitative analysis Methods 0.000 description 4
- 150000003839 salts Chemical class 0.000 description 4
- 239000004575 stone Substances 0.000 description 4
- 238000011156 evaluation Methods 0.000 description 3
- 230000008901 benefit Effects 0.000 description 2
- 239000003245 coal Substances 0.000 description 2
- 238000007796 conventional method Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- -1 life Chemical compound 0.000 description 2
- 238000000691 measurement method Methods 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 241001269238 Data Species 0.000 description 1
- 230000035508 accumulation Effects 0.000 description 1
- 238000009825 accumulation Methods 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 230000002776 aggregation Effects 0.000 description 1
- 238000004220 aggregation Methods 0.000 description 1
- 238000012512 characterization method Methods 0.000 description 1
- 238000002224 dissection Methods 0.000 description 1
- 238000013210 evaluation model Methods 0.000 description 1
- 230000002349 favourable effect Effects 0.000 description 1
- 238000005429 filling process Methods 0.000 description 1
- 230000013632 homeostatic process Effects 0.000 description 1
- 238000000265 homogenisation Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000006386 neutralization reaction Methods 0.000 description 1
- 238000011017 operating method Methods 0.000 description 1
- 230000008092 positive effect Effects 0.000 description 1
- 238000011158 quantitative evaluation Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 239000000243 solution Substances 0.000 description 1
- 238000006467 substitution reaction 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/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
-
- 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)
- Business, Economics & Management (AREA)
- Human Resources & Organizations (AREA)
- Economics (AREA)
- Strategic Management (AREA)
- Physics & Mathematics (AREA)
- Tourism & Hospitality (AREA)
- Theoretical Computer Science (AREA)
- Marketing (AREA)
- General Physics & Mathematics (AREA)
- General Business, Economics & Management (AREA)
- Quality & Reliability (AREA)
- Animal Husbandry (AREA)
- Game Theory and Decision Science (AREA)
- Operations Research (AREA)
- Entrepreneurship & Innovation (AREA)
- Life Sciences & Earth Sciences (AREA)
- Agronomy & Crop Science (AREA)
- Development Economics (AREA)
- Marine Sciences & Fisheries (AREA)
- Mining & Mineral Resources (AREA)
- Health & Medical Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明提出一种定量评价断控天然气藏富集程度的方法，属于气田勘探技术领域，包括以下步骤：选取与断控天然气藏富集程度相关的六个地质因素：源岩‑断层接触处平均生气强度、源岩‑断层接触面积、成藏期源‑储剩余压力差、天然气纵向运移距离、断层倾角余弦值、保存条件参数；选取多个已探明气藏，获取其对应的六个因素和富集程度参数，绘制散点图，通过拟合函数获得单因素与富集程度的关系式；将获得的关系式通过多元线性回归建立断控天然气藏富集程度的多因素预测模型；初选勘探潜力区，获取其对应的六个因素数据，结合多因素预测模型进行勘探潜力区的选择。本发明建立了多因素综合预测模型，能够有效提高断控天然气藏富集程度的定量评价精度。
Description
技术领域
本发明属于气田勘探技术领域，尤其涉及一种定量评价断控天然气藏富集程度的方法。
背景技术
断控天然气藏是指以断层为主要输导通道的烃源岩层系外圈闭中聚集形成的下生上储式天然气藏，目前中国东部盆地发现的气藏多数属于此类。评价断控天然气藏富集程度的传统方法，主要是基于气田勘探实践结果，对生、运、储、聚、散等成藏要素与天然气富集之间的关系进行剖析，从而对断控天然气藏富集程度及控制因素进行宏观定性分析。然而，随着气田勘探程度的不断提高，传统的定性分析已不能满足现今天然气藏精细勘探的要求，因此，迫切需要能够定量评价断控天然气藏富集程度的方法。
近十年来，越来越多的学者致力于天然气藏富集程度的单一影响因素的半定量和定量化分析。例如：孙明亮分别对盖层厚度和盖层排替压力两个因素进行了单一因素的定量分析，公开了一组形成天然气藏的定量指标，具体的：形成中效天然气藏需要盖层厚度大于40m、盖层排替压力大于15MPa，而形成高效天然气藏则需要厚度超过100m的直接盖层、盖层排替压力不应小于20MPa。此外，柳广弟根据对典型气藏的解剖和大量气田的统计，指出天然气的成藏过程主要受气藏成藏期的源-储剩余压力差、输导体系的类型和输导效率、盖层厚度和排替压力等几方面因素的控制，并对这些因素分别进行了单一因素的定量分析，得到的指标为：形成中、高效气藏需要成藏期的源-储剩余压力差大于25MPa，输导体系为汇聚型，盖层厚度大于40m，排替压力大于15MPa。
然而，现有技术中的这些半定量和定量化分析仍局限于单一因素的定量化分析，由于断控天然气藏所处的地质环境复杂，其富集程度是由多种成藏要素共同控制的结果，且不同因素的影响贡献存在差异，因而这种单一因素的定量评价精度低，不能有效指导断控天然气藏的勘探。因此，如何提供一种多因素综合预测模型，提高断控天然气藏富集程度的定量评价精度，是当前急需解决的一项技术难题。
发明内容
本发明针对上述的现有断控天然气藏富集程度的定量评价方法局限于单一因素、评价精度低的技术问题，提出一种定量评价断控天然气藏富集程度的方法，建立了多因素综合预测模型，能够有效提高断控天然气藏富集程度的定量评价精度。
为了达到上述目的，本发明采用的技术方案为：
一种定量评价断控天然气藏富集程度的方法，包括以下步骤：
选取影响因素：选取与断控天然气藏富集程度相关的六个地质因素，分别为源岩-断层接触处平均生气强度X1、源岩-断层接触面积X2、成藏期源-储剩余压力差X3、天然气纵向运移距离X4、断层倾角余弦值X5、保存条件参数X6；
分析单因素相关性：选取多个属于同一构造背景、相同类型的已探明断控天然气藏为样本气藏，分别获取多个样本气藏对应的源岩-断层接触处平均生气强度X1、源岩-断层接触面积X2、成藏期源-储剩余压力差X3、天然气纵向运移距离X4、断层倾角余弦值X5、保存条件参数X6以及样本气藏的富集程度参数的数据，分别绘制六个因素与富集程度的散点图，通过拟合函数获得单因素与富集程度的关系式：Fi＝D(Xi)，i＝1,2,……,6；
建立预测模型：将获得的六个单因素与富集程度的关系式通过多元线性回归法建立断控天然气藏富集程度的多因素预测模型:
其中，ai为权重系数，通过多元线性回归法获得；
选择勘探潜力区：初选勘探潜力区，获取初选潜力区对应的源岩-断层接触处平均生气强度X1、源岩-断层接触面积X2、成藏期源-储剩余压力差X3、天然气纵向运移距离X4、断层倾角余弦值X5、保存条件参数X6的数据，代入所述多因素预测模型计算得到每个初选潜力区的富集程度，选择其中富集程度高的地区作为重点勘探潜力区。
作为优选，所述的富集程度参数为储量。
作为优选，所述源岩-断层接触处平均生气强度X1通过将烃源岩层与断层的接触处分成若干小段，由每一小段处的生气强度相加求平均获得。
作为优选，所述保存条件参数X6采用断接厚度参数表征，通过对断接厚度分类赋值的方法获得。
与现有技术相比，本发明的优点和积极效果在于：
1、本发明充分考虑了供气量、成藏动力、断裂输导分流作用、封盖保存作用等多种成藏要素对断控天然气藏富集程度的影响，并具体选择了源岩-断层接触处平均生气强度、源岩-断层接触面积、成藏期源-储剩余压力差、天然气纵向运移距离、断层倾角余弦值、保存条件参数六个因素，建立了多因素综合预测模型，并基于预测模型进行有利勘探区优选，能够有效提高断控天然气藏富集程度的定量评价精度；
2、本发明提供的定量评价断控天然气藏富集程度的方法，可以为断控天然气藏的勘探部署提供指导，降低勘探风险，具有显著的经济效益。
附图说明
图1为断控天然气藏所处地质环境的纵向地质剖面示意图；
图2为标注有生气强度的烃源岩层的水平地质剖面示意图；
图3为本发明实施例所提供的定量评价断控天然气藏富集程度的方法的流程图；
图4为本发明实施例1中煤成气探明储量与源岩-断层接触处平均生气强度X1的关系图；
图5为本发明实施例1中煤成气探明储量与源岩-断层接触面积X2的关系图；
图6为本发明实施例1中煤成气探明储量与成藏期源-储剩余压力差X3的关系图；
图7为本发明实施例1中煤成气探明储量与天然气纵向运移距离X4的关系图；
图8为本发明实施例1中煤成气探明储量与断层倾角余弦值X5的关系图；
图9为本发明实施例1中煤成气探明储量与保存条件参数X6的关系图；
图10为本发明实施例1所提供的多因素预测模型的计算值与实际值之间的关系图；
以上各图中：1、烃源岩层；2、断层；3、储集层；4、盖层；5、断控天然气藏。
具体实施方式
下面将对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
本发明实施例提供了一种定量评价断控天然气藏富集程度的方法，包括以下步骤：
S1选取影响因素：选取与断控天然气藏富集程度相关的六个地质因素，分别为源岩-断层接触处平均生气强度X1、源岩-断层接触面积X2、成藏期源-储剩余压力差X3、天然气纵向运移距离X4、断层倾角余弦值X5、保存条件参数X6。
在本步骤中，影响因素的选取是基于天然气聚散动平衡理论，即任何气藏的形成均是天然气的“聚”和“散”之间的动态平衡过程，“聚”大于“散”才能形成有效的天然气聚集。对于断控天然气藏而言，“聚”的量主要取决于供气量和成藏动力，供气量越多、成藏动力越大，天然气充注进入气藏的量越多；“散”的量则主要取决于运移过程中断裂输导分流作用和封盖保存作用，断裂输导分流作用越弱、封盖保存作用越强，天然气的散失量越少。因此，断控天然气藏的富集程度主要受供气量、成藏动力、断裂输导分流作用、封盖保存作用四个成藏要素的共同控制。
如图1和图2所示，供气量由烃源岩层1与断层2接触处烃源岩的生气强度与接触面积共同控制，生气强度越大、接触面积越大，则气源供给越充足，因此，选取了源岩-断层接触处平均生气强度X1和源岩-断层接触面积X2两个影响因素。成藏动力主要取决于成藏期烃源岩层1的剩余压力与储集层3中圈闭储层的孔隙流体剩余压力的差值，此压力差越大，成藏动力越大，因此，选取了成藏期源-储剩余压力差X3作为第三个影响因素。断裂输导分流作用由输导距离及输导速度共同控制，其中，输导距离取决于天然气纵向运移距离，天然气纵向运移距离越远，输导分流作用越强；而输导速度受断层2的倾角大小的影响，断层2的倾角越小，浮力沿断层方向分解力越小，输导速度越慢，天然气沿断层纵向输导能力越弱、输导分流作用越强。因而，天然气纵向运移距离越短、断层2的倾角越大，断裂输导分流作用越弱，因此，选取了天然气纵向运移距离X4、断层倾角余弦值X5分别作为第四个和第五个影响因素。封盖保存作用主要取决于盖层4中泥岩(或盐膏岩)的保存条件参数，例如厚度参数，厚度越大，盖层4封盖天然气的能力越强，封盖保存作用越强，因此，选取了保存条件参数X6作为第六个影响因素。
S2分析单因素相关性：选取多个属于同一构造背景、相同类型的已探明断控天然气藏为样本气藏，分别获取多个样本气藏对应的源岩-断层接触处平均生气强度X1、源岩-断层接触面积X2、成藏期源-储剩余压力差X3、天然气纵向运移距离X4、断层倾角余弦值X5、保存条件参数X6以及样本气藏的富集程度参数的数据，分别绘制六个因素与富集程度的散点图，通过拟合函数获得单因素与富集程度的关系式：Fi＝D(Xi)，i＝1,2,……,6。
在本步骤中，源岩-断层接触处平均生气强度X1是通过测得的生气强度进一步求平均获得，其单位为m3/km2，生气强度的测量方法采用本领域的常规方法，例如：生烃动力学法、热模拟实验法和盆地模拟法等；可以理解的是，本实施例并不局限于上述所列举的生气强度的测量方法，还可以为本领域技术人员所熟知的其它合理方法。
如图1和图2所示，烃源岩层1与断层2的纵向接触长度如图1中标注的L1所示,烃源岩层1与断层2的平面接触长度如图2中标注的L2所示,源岩-断层接触面积X2即为源岩-断层纵向接触长度L1与平面接触长度L2的乘积，其单位为km2，其中，L1和L2可通过地震解释的平剖面数据结合微分叠加的方法获得，本领域技术人员还可以采用其他合理的方法定量表征源岩-断层接触面积X2。
成藏期源-储剩余压力差X3，即成藏期烃源岩层1的剩余压力与储集层3中圈闭储层的孔隙流体剩余压力的差值，其单位为MPa，其中，成藏期烃源岩层1的剩余压力由排烃压力替代表征，圈闭储层的孔隙流体剩余压力可以通过流体包裹体分析和欠压实等多种方法综合确定。可以理解的是，本领域技术人员还可以采用其他合理的方法获取成藏期源-储剩余压力差X3。
天然气纵向运移距离X4为烃源岩层1的上边界与储集层3的上边界之间的距离，如图2中标注的H所示，其单位为km，可通过地震解释的平剖面数据结合微分叠加的方法获得，本领域技术人员还可以采用其他合理的方法获取此参数。
断层2与水平面的夹角即为断层倾角，如图2中标注的θ所示，断层倾角余弦值X5即为断层倾角θ的余弦值，其中，断层倾角θ的角度通过桌面直尺量角器在地震剖面上读取。当然，本领域技术人员还可以采用其他合理的方法来测量断层倾角θ的角度。
保存条件参数X6为盖层4中泥岩(或盐膏岩)的保存条件参数，可选厚度参数，本领域技术人员还可以选取其他合理的保存条件参数。
所述气藏的富集程度参数可采用储量来表征，本领域技术人员还可以采用其他合理的富集程度参数。
此外，还需要说明的是，根据六个单因素与富集程度关系式的相关系数的大小，可分析富集程度与各因素之间的关系密切程度，相关系数的值越大，则富集程度与该因素的关系越密切，影响越大，由此可排除相关系数很低的因素。
S3建立预测模型：将获得的六个单因素与富集程度的关系式通过多元线性回归法建立断控天然气藏富集程度的多因素预测模型:
其中，ai为权重系数，通过多元线性回归法获得。
在本步骤中，需要说明的是，权重系数ai为通过多元线性回归法获得的六个单因素均一化变量的最优组合回归系数。此外，通过所得多因素预测模型相关系数的大小，可以评价模型的合理性和可靠性。
S4选择勘探潜力区：初选勘探潜力区，获取初选潜力区对应的源岩-断层接触处平均生气强度X1、源岩-断层接触面积X2、成藏期源-储剩余压力差X3、天然气纵向运移距离X4、断层倾角余弦值X5、保存条件参数X6的数据，代入所述多因素预测模型计算得到每个初选潜力区的富集程度，选择其中富集程度高的地区作为重点勘探潜力区。
在本步骤中，需要说明的是，勘探潜力区的初选可采用本领域技术人员所熟知的评价断控天然气藏富集程度的定性分析方法，选择成藏条件有利的地区为勘探潜力区。进一步通过将初选潜力区对应的六个因素的数据代入所述多因素预测模型中，可计算获得每个初选潜力区的富集程度，选择其中富集程度高的地区作为重点勘探潜力区。
本发明实施例所提供的定量评价断控天然气藏富集程度的方法的流程图如图3所示，本发明实施例针对不同地质背景的待勘探区，选择同一构造背景下、相同类型的已探明断控天然气藏为研究区，建立了六个单因素与富集程度的关系式，进一步通过多元线性回归法，建立了特定研究区天然气藏富集程度的定量表征模型——多因素预测模型，模型中不同因素对应的权重系数即可反映所有不同因素之间以及各因素与富集程度参数之间的模切程度的大小。同时，利用该模型可对同一研究区内待勘探区的天然气富集程度进行预测，从而优选出重点勘探潜力区。这种精细定量评价方式可以有效提高断控天然气藏富集程度的定量评价精度，提高勘探成功率，同时可降低断控天然气藏的勘探风险。
在一可选实施例中，所述的富集程度参数为储量，其单位为m3。由于储量是易测定且具有代表性的富集程度参数，因此，在本实施例中，采用储量表征断控天然气藏的富集程度，由此获得的多因素预测模型更具有代表性。可以理解的是，本领域技术人员还可以采用其他合理的富集程度参数。
在一可选实施例中，所述源岩-断层接触处平均生气强度X1通过将烃源岩层1与断层2的接触处分成若干小段，由每一小段处的生气强度相加求平均获得，采用此方法获得的数据更为准确。
在一可选实施例中，为了便于获取保存条件参数的数据，所述保存条件参数X6采用断接厚度参数表征，通过对断接厚度分类赋值的方法获得。在本实施例中，所述保存条件参数X6采用断接厚度参数表征，即盖层4的厚度与断层2的断距差值，如图2中标注的h所示。由于在一定厚度范围内，泥岩(或盐膏岩)的保存作用相差较小，因此，可以通过对断接厚度分类赋值法获取保存条件系数X6。具体的，按照断接厚度大小对盖层4进行分类，即A类：h>500m、B类：h＝500-300m、C类：h＝300-100m、D类：h＝100-40m、E类：h＝0-40m，分别对五类盖层的保存条件系数X6进行赋值，即A类：X6＝4、B类：X6＝3、C类：X6＝2、D类：X6＝1、E类：X6＝0。采用分类赋值法获取保存条件参数X6，不需要测定盖层中泥岩(或盐膏岩)的保存作用的具体数值，简化了操作步骤，降低了工作量、耗时短，且成本低。
为了更清楚详细地介绍本发明实施例所提供的定量评价断控天然气藏富集程度的方法，下面将结合具体实施例进行描述。
实施例1
以渤海湾盆地东濮凹陷煤成气藏为例，具体说明定量评价断控天然气藏富集程度的方法。渤海湾盆地东濮凹陷煤成气的烃源岩为下部石炭-二叠系煤系地层，主要富集层系为古近系，区内主干基底断裂为其主要的垂向输导通道，属于典型的断控天然气藏。具体操作步骤如下：
1、选取影响因素：选取与断控天然气藏富集程度相关的六个地质因素，分别为源岩-断层接触处平均生气强度X1、源岩-断层接触面积X2、成藏期源-储剩余压力差X3、天然气纵向运移距离X4、断层倾角余弦值X5、保存条件参数X6；
2、分析单因素相关性：选取文23、户部寨、马厂、白庙、卫城、濮城6个已知断控天然气藏，分别获取其对应的源岩-断层接触处平均生气强度X1、源岩-断层接触面积X2、成藏期源-储剩余压力差X3、天然气纵向运移距离X4、断层倾角余弦值X5、保存条件参数X6以及其煤成气探明储量的数据，分别绘制六个因素与富集程度的散点图，如图4-图9所示，通过拟合函数获得以下具体单因素与富集程度的关系式：
F1＝2.2012*EXP(0.0747*X1)，R2＝0.8846 (1)
F2＝3.6529*EXP(0.0521*X2)，R2＝0.5784 (2)
F3＝5.8874*EXP(0.0553*X3)，R2＝0.4334 (3)
F4＝-138.2*ln(X4)+282.09，R2＝0.2727 (4)
F5＝-430.9*X5+410.95，R2＝0.5466 (5)
F6＝6.5609*X6 2.0762，R2＝0.8567 (6)
其中，源岩-断层接触处平均生气强度X1通过将烃源岩层1与断层2的接触处分成若干小段，由每一小段处的生气强度相加求平均获得；源岩-断层接触面积X2为源岩-断层纵向接触长度L1与平面接触长度L2的乘积，其中，L1和L2通过地震解释的平剖面数据结合微分叠加的方法获得；成藏期源-储剩余压力差X3为成藏期烃源岩层1的剩余压力与储集层3中圈闭储层的孔隙流体剩余压力的差值，其中，成藏期烃源岩层1的剩余压力由排烃压力替代表征，圈闭储层的孔隙流体剩余压力通过流体包裹体分析和欠压实方法综合确定；天然气纵向运移距离X4通过地震解释的平剖面数据结合微分叠加的方法获得；断层倾角余弦值X5即为断层倾角θ的余弦值，断层倾角θ的角度通过桌面直尺量角器在地震剖面上读取；保存条件参数X6采用断接厚度参数表征，通过对断接厚度分类赋值的方法获得。
从式(1)-式(6)可以看出：渤海湾盆地东濮凹陷煤成气藏的煤成气探明储量与源岩-断层接触处平均生气强度X1、源岩-断层接触面积X2、成藏期源-储剩余压力差X3以及保存条件参数X6呈正相关关系，与天然气纵向运移距离X4及断层倾角余弦值X5呈负相关关系。由相关系数R2的大小可以看出各因素对煤成气探明储量影响的强弱顺序为：源岩-断层接触处平均生气强度X1>保存条件参数X6>源岩-断层接触面积X2>断层倾角余弦值X5>成藏期源-储剩余压力差X3>天然气纵向运移距离X4，由此表明：源岩生气强度和保存条件是东濮凹陷煤成气成藏的首要控制因素，而断裂输导分流作用对其控制作用相对较弱。
3、建立预测模型：将式(1)-式(6)通过多元线性回归法建立断控天然气藏富集程度的多因素预测模型为:
F＝0.894F1-1.519F2+0.121F3-1.474F4-0.336F5+0.797F6+55.484，R2＝0.9868(7)
从相关系数R2＝0.9868可以看出，拟合的公式具有较高的相关性。从图10中可以看出，通过多因素预测模型式(7)得到的计算值与实际值的误差相对较小，误差均在允许范围之内，保证该模型的建立合理、可靠。
4、选择勘探潜力区：通过对东濮凹陷煤成气待勘探地区的成藏要素与富集程度之间的关系进行定性分析，初步选取凹陷内的桥口、方里集、文中、文西、刘庄和徐集六个勘探潜力区，获取初选潜力区对应的源岩-断层接触处平均生气强度X1、源岩-断层接触面积X2、成藏期源-储剩余压力差X3、天然气纵向运移距离X4、断层倾角余弦值X5、保存条件参数X6的数据，代入式(7)中，计算得到六个地区的预测煤成气储量，如表1所示：
表1六个初选潜力区测定数据及预测煤成气探明储量的数据统计表
文中地区、方里集地区预测煤成气储量相对较高，分别达74.80亿方和55.14亿方，文西、刘庄、桥口地区次之，预测储量分别为34.57、27.08和15.84亿方，徐集地区最小，预测储量仅为8.84亿方，因此，优选文中和方里集两个预测储量较高的地区作为重点勘探潜力区。
Claims (2)
1.一种定量评价断控天然气藏富集程度的方法，其特征在于，包括以下步骤：
选取影响因素：选取与断控天然气藏富集程度相关的六个地质因素，分别为源岩-断层接触处平均生气强度X1、源岩-断层接触面积X2、成藏期源-储剩余压力差X3、天然气纵向运移距离X4、断层倾角余弦值X5、保存条件参数X6；
分析单因素相关性：选取多个属于同一构造背景、相同类型的已探明断控天然气藏为样本气藏，分别获取多个样本气藏对应的源岩-断层接触处平均生气强度X1、源岩-断层接触面积X2、成藏期源-储剩余压力差X3、天然气纵向运移距离X4、断层倾角余弦值X5、保存条件参数X6以及样本气藏的富集程度参数的数据，分别绘制六个因素与富集程度的散点图，通过拟合函数获得单因素与富集程度的关系式：Fi＝D(Xi)，i＝1,2,……,6；
建立预测模型：将获得的六个单因素与富集程度的关系式通过多元线性回归法建立断控天然气藏富集程度的多因素预测模型:
其中，ai为权重系数，通过多元线性回归法获得；
选择勘探潜力区：初选勘探潜力区，获取初选潜力区对应的源岩-断层接触处平均生气强度X1、源岩-断层接触面积X2、成藏期源-储剩余压力差X3、天然气纵向运移距离X4、断层倾角余弦值X5、保存条件参数X6的数据，代入所述多因素预测模型计算得到每个初选潜力区的富集程度，选择其中富集程度高的地区作为重点勘探潜力区；
所述富集程度参数为储量；
所述保存条件参数X6采用断接厚度参数表征，通过对断接厚度分类赋值的方法获得；赋值方法为：按照断接厚度大小对盖层进行分类，A类：h>500m、B类：h＝500-300m、C类：h＝300-100m、D类：h＝100-40m、E类：h＝0-40m，分别对五类盖层的保存条件系数X6进行赋值，A类：X6＝4、B类：X6＝3、C类：X6＝2、D类：X6＝1、E类：X6＝0。
2.根据权利要求1所述的定量评价断控天然气藏富集程度的方法，其特征在于：所述源岩-断层接触处平均生气强度X1通过将烃源岩层与断层的接触处分成若干小段，由每一小段处的生气强度相加求平均获得。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710039408.3A CN106651060B (zh) | 2017-01-19 | 2017-01-19 | 一种定量评价断控天然气藏富集程度的方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710039408.3A CN106651060B (zh) | 2017-01-19 | 2017-01-19 | 一种定量评价断控天然气藏富集程度的方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106651060A CN106651060A (zh) | 2017-05-10 |
| CN106651060B true CN106651060B (zh) | 2019-12-03 |
Family
ID=58841912
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710039408.3A Active CN106651060B (zh) | 2017-01-19 | 2017-01-19 | 一种定量评价断控天然气藏富集程度的方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106651060B (zh) |
Families Citing this family (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109064351B (zh) * | 2018-10-15 | 2021-12-10 | 中国石油大学(华东) | 一种潜山油气充注能力的定量评价方法 |
| CN109948904B (zh) * | 2019-02-21 | 2022-11-18 | 西南石油大学 | 一种求取含油气区带评价参数权重的方法 |
| CN112177605B (zh) * | 2020-09-16 | 2022-04-22 | 广州海洋地质调查局 | 一种海域天然气水合物有利聚集区主控因素的确定方法 |
| CN112330109A (zh) * | 2020-10-22 | 2021-02-05 | 中国石油大学(华东) | 油气勘探区带地质综合评价的量化表征方法、系统、设备 |
| CN114544917B (zh) * | 2020-11-24 | 2024-01-30 | 中国石油天然气股份有限公司 | 原油裂解气藏天然气失散量确定方法及装置 |
| CN117236036A (zh) * | 2023-09-25 | 2023-12-15 | 中国石油天然气集团有限公司 | 一种定量预测气藏中氦气含量的方法及装置 |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104750884A (zh) * | 2013-12-26 | 2015-07-01 | 中国石油化工股份有限公司 | 基于多因素非线性回归的页岩油气富集指数定量评价方法 |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2015001558A1 (en) * | 2013-07-01 | 2015-01-08 | Salespredict Sw Ltd. | System and method for predicting sales |
-
2017
- 2017-01-19 CN CN201710039408.3A patent/CN106651060B/zh active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104750884A (zh) * | 2013-12-26 | 2015-07-01 | 中国石油化工股份有限公司 | 基于多因素非线性回归的页岩油气富集指数定量评价方法 |
| CN104750884B (zh) * | 2013-12-26 | 2018-02-02 | 中国石油化工股份有限公司 | 基于多因素非线性回归的页岩油气富集指数定量评价方法 |
Non-Patent Citations (1)
| Title |
|---|
| 《页岩含气量主控因素及其定量模型》;李明明;《中国优秀硕士学位论文全文数据库 基础科学辑》;中国学术期刊（光盘版）电子杂志社;20160415(第4期);摘要，正文第47-67页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106651060A (zh) | 2017-05-10 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106651060B (zh) | 一种定量评价断控天然气藏富集程度的方法 | |
| CN110288258B (zh) | 一种高含水油藏剩余油挖潜方法 | |
| CN103806911B (zh) | 利用随钻钻井液录井资料进行储层流体类型的判别方法 | |
| Tianjian et al. | Classification and characterization of barrier-intercalation in sandy braided river reservoirs: Taking Hegli Oilfield of Muglad Basin in Sudan as an example | |
| CN103883322B (zh) | 页岩气储层勘探方法及装置 | |
| CN105201490B (zh) | 一种页岩层段岩石相分析方法 | |
| CN110412661A (zh) | 细粒岩油气藏甜点段优势段簇的评价方法及装置 | |
| CN111046585A (zh) | 基于多元线性回归分析的页岩气甜点预测方法 | |
| Li | Geophysical exploration technology: Applications in lithological and stratigraphic reservoirs | |
| CN108956952A (zh) | 盐间泥页岩层系油藏地质甜点的预测方法及勘探方法 | |
| CN105092820A (zh) | 一种生物碎屑灰岩储层非均质性研究方法 | |
| CN104991286A (zh) | 一种基于沉积模式的沉积相表征方法 | |
| Wang et al. | Quantitative evaluation of transport efficiency of fault-reservoir composite migration pathway systems in carbonate petroliferous basins | |
| CN107292435A (zh) | 一种基于含油高度定量预测的断块圈闭评价方法 | |
| Qi'an et al. | Accumulation and exploration of petroleum reservoirs in west slope of northern Songliao Basin, China | |
| CN104516025A (zh) | 碳酸盐储层物性随钻分类和评价方法 | |
| CN109653738A (zh) | 一种评价复杂改造区页岩气保存有效性的方法 | |
| CN106372751A (zh) | 一种相对古水深指标参数的测井预测方法 | |
| Förster et al. | Spatial analysis of temperature (BHT/DST) data and consequences for heat-flow determination in sedimentary basins | |
| Hemmings-Sykes | The influence of faulting on hydrocarbon migration in the Kupe area, south Taranaki Basin, New Zealand | |
| Feng et al. | Discovery of “Enveloping surface of oil and gas overpressure migration” in the Songliao Basin and its bearings on hydrocarbon migration and accumulation mechanisms | |
| CN110244356A (zh) | 一种构造破坏形成的后期油气藏的判识方法 | |
| Yu et al. | Examining and applying the theory of “exploring petroleum inside source kitchens” for continental shale oil: A case study from the Kong 2 member of the Cangdong sag in the Bohai Bay Basin, China | |
| CN110930020B (zh) | 非常规油气资源的经济可采资源量确定方法 | |
| Liang et al. | Opportunities and challenges for gas coproduction from coal measure gas reservoirs with coal‐shale‐tight sandstone layers: A review |
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 |