CN110318744A - 一种用于预测页岩气资源的方法 - Google Patents

一种用于预测页岩气资源的方法 Download PDF

Info

Publication number
CN110318744A
CN110318744A CN201810289872.2A CN201810289872A CN110318744A CN 110318744 A CN110318744 A CN 110318744A CN 201810289872 A CN201810289872 A CN 201810289872A CN 110318744 A CN110318744 A CN 110318744A
Authority
CN
China
Prior art keywords
point
equal
interpolation
less
shale
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.)
Granted
Application number
CN201810289872.2A
Other languages
English (en)
Other versions
CN110318744B (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.)
China Petroleum and Chemical Corp
Sinopec Exploration and Production Research Institute
China Petrochemical Corp
Original Assignee
Sinopec Exploration and Production Research Institute
China Petrochemical Corp
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 Sinopec Exploration and Production Research Institute, China Petrochemical Corp filed Critical Sinopec Exploration and Production Research Institute
Priority to CN201810289872.2A priority Critical patent/CN110318744B/zh
Publication of CN110318744A publication Critical patent/CN110318744A/zh
Application granted granted Critical
Publication of CN110318744B publication Critical patent/CN110318744B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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
    • E21B49/00—Testing the nature of borehole walls; Formation testing; Methods or apparatus for obtaining samples of soil or well fluids, specially adapted to earth drilling or wells

Landscapes

  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Geology (AREA)
  • Mining & Mineral Resources (AREA)
  • Physics & Mathematics (AREA)
  • Environmental & Geological Engineering (AREA)
  • Fluid Mechanics (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Geochemistry & Mineralogy (AREA)
  • Investigating Or Analyzing Materials Using Thermal Means (AREA)

Abstract

本发明提出了一种用于预测页岩气资源的方法，该方法利用页岩气富集主控因素等地质资料，按照页岩气富集条件评价和经济可采分级评价标准，优选页岩气资源禀赋较高的有利区块以及具有商业开采价值的页岩气储层。在此基础上针对优选的富集区分别对页岩厚度、含气性等参数进行不确定性评价，从而准确预测出页岩气资源量。

Description

一种用于预测页岩气资源的方法
技术领域
本发明涉及地质勘探领域，特别涉及页岩油气勘探领域。
背景技术
为了准确、科学、快速评价非常规油气勘探项目，探寻页岩油气勘探开发接替区块，需要建立非常规油气勘探项目特有的评价方法和技术，展开适应性配套技术攻关，满足页岩油气多阶段的勘探开发格局的需求。
通常，在决策进入钻井勘探阶段之前的页岩气地质资源量评价，需要通过区域地质调查、钻探取样分析、页岩气成藏地质条件评价，搞清页岩气含气系统是否普遍存在，优选出页岩气勘探目的层与目标区，评价页岩气资源勘探的潜力与规模；在此基础上，通过选区评价工作，圈定有利区以及具有商业开采价值的页岩气“甜点”。但针对页岩储层非均质性强和地质参数(富有机质页岩厚度、含气性等参数)的变化，传统的页岩气有利区选区仍采用统一标准，无法实现油气富集区和经济可采区的分级评价；资源评价方法直接采用单一参数利用体积法和地质参数的单一均值直接计算页岩气地质资源量，利用这种方法会导致页岩气地质资源量评价误差较大，忽视非均质性、采用单一评价参数给出确定性评价结果会导致有利区选区和资源评价结果可信度降低。
现有的页岩油气资源评价方法主要存在以下问题：(1)无法实现页岩气富集有利区和经济可采区的分级评价并在此基础上进行资源量评价；(2)目前广泛应用的体积法预测结果多为确定性数据，缺少不确定评价结果。
发明内容
为解决上述技术问题，本发明提出了一种用于预测页岩气资源的方法，其包括以下步骤：
步骤S100：对待评价区域内的多个处进行测井以获得每个井位处目的储层的有机碳含量、烃源岩热成熟度指标、页岩厚度、页岩密度、页岩含气量、埋深；
对所有页岩密度求平均以获得平均页岩密度；
步骤S200：对待评价区域进行均匀的网格划分，待评价区域内设置有井位的单元格为已知单元格，其他为未知单元格；
步骤S300：将每个已知单元格中的一个井位处的有机碳含量、烃源岩热成熟度指标、页岩厚度、页岩含气量、埋深的值依次分别作为目的储层在该已知单元格范围内有机碳含量、烃源岩热成熟度指标、页岩厚度、页岩含气量、埋深的均值；
步骤S400：随机构建模拟所有未知单元格的地质属性的顺序，按顺序依次对每个未知单元格的地质属性进行模拟，
模拟第一个未知单元格的地质属性则进行步骤S410，
模拟第二个未知单元格至最后一个未知单元格的地质属性时，每对一个未知单元格的地质属性进行模拟则进行步骤S420；
步骤S410：根据每个已知单元格的有机碳含量的均值、烃源岩热成熟度指标的均值、页岩厚度的均值、页岩含气量的均值、埋深的均值、各个单元格的位置和当前待计算的未知单元格的位置采用克里金插值法计算出该未知单元格中有机碳含量插值、烃源岩热成熟度指标插值、页岩厚度插值、页岩含气量插值、埋深插值；
步骤S420：根据每个已知单元格的有机碳含量的均值、烃源岩热成熟度指标的均值、页岩厚度的均值、页岩含气量的均值、埋深的均值，已经计算出的未知单元格中有机碳含量插值、烃源岩热成熟度指标插值、页岩厚度插值、页岩含气量插值、埋深插值，以及各个单元格的位置和当前待计算的未知单元格的位置采用克里金插值法计算出该未知单元格中有机碳含量插值、烃源岩热成熟度指标插值、页岩厚度插值、页岩含气量插值、埋深插值；
步骤S500：对每个已知单元格的有机碳含量的均值进行评价获得有机碳评分值，对每个未知单元格的有机碳含量插值进行评价获得有机碳评分值，对每个已知单元格的烃源岩热成熟度指标的均值进行评价获得热成熟度评分值，对每个未知单元格的烃源岩热成熟度指标插值进行评价获得热成熟度评分值，对每个已知单元格的页岩厚度的均值进行评价获得厚度评分值，对每个未知单元格的页岩厚度插值进行评价获得厚度评分值，对每个已知单元格的埋深的均值进行评价获得经济可采性评分，对每个未知单元格的埋深插值进行评价获得经济可采性评分；
步骤S600：将每个已知单元格的有机碳评分值、热成熟度评分值、厚度评分值相乘获得每个已知单元格的页岩气富集概率；
将每个未知单元格的有机碳评分值、热成熟度评分值、厚度评分值相乘获得每个未知单元格的页岩气富集概率；
步骤S700：将有机碳评分值、热成熟度评分值、厚度评分值、经济可采性评分均大于或等于0.5分的单元格的集合划为有利区；
步骤S800：计算出有利区内每个已知单元格对应区域中的天然气地质资源量和有利区内每个未知单元格对应区域中的天然气地质资源量，
将有利区内每个已知单元格对应区域中的天然气地质资源量与有利区内每个未知单元格对应区域中的天然气地质资源量进行累加获得累积天然气地质资源量；
步骤S900：重复步骤S400至步骤S800预设次数，以获得预设个数的累积天然气地质资源量；
步骤S1000：以累积天然气地质资源量数值从大到小的顺序对所有累积天然气地质资源量从序号1开始进行连续编号，按照算式(1)计算每个序号对应的超概率，
超概率＝(序号/N+(序号-1)/N)/2 (1)
其中，N是累积天然气地质资源量的个数，
将每个累积天然气地质资源量及其对应的超概率以坐标点的形式标记在以横轴表示累积天然气地质资源量数值、纵轴表示超概率的直角坐标系上，将直角坐标系上相邻的两个点用直线段连接起来得到超概率曲线，将超概率曲线上超概率等于百分之五十所对应的累积天然气地质资源量作为待评价区域的天然气地质资源量的最终估值。
在一个具体的实施例中，在步骤S500中，
已知单元格的有机碳评分值评分标准为：当有机碳含量的均值小于0.4％时获得0分，当有机碳含量的均值小于0.5％大于等于0.4％时随机获得0分和0.5分中的一个分值，当有机碳含量的均值大于等于0.5％小于0.9％时获得0.5分，当有机碳含量的均值大于等于0.9％小于1％时随机获得0.5分和0.8分中的一个分值，当有机碳含量的均值大于等于1％小于1.4％时获得0.8分，当有机碳含量的均值大于等于1.4％小于1.5％时随机获得0.8分和1分中的一个分值，当有机碳含量的均值大于等于1.5％时获得1分；
未知单元格的有机碳评分值评分标准为：当有机碳含量插值小于0.4％时获得0分，当有机碳含量插值小于0.5％大于等于0.4％时随机获得0分和0.5分中的一个分值，当机碳含量插值大于等于0.5％小于0.9％时获得0.5分，当有机碳含量插值大于等于0.9％小于1％时随机获得0.5分和0.8分中的一个分值，当机碳含量插值大于等于1％小于1.4％时获得0.8分，当机碳含量插值大于等于1.4％小于1.5％时随机获得0.8分和1分中的一个分值，当有机碳含量插值大于等于1.5％时获得1分；
已知单元格的热成熟度评分值评分标准为：当烃源岩热成熟度指标的均值小于0.6％时获得0分，当烃源岩热成熟度指标的均值小于0.7％大于等于0.6％时随机获得0分和0.5分中的一个分值，当烃源岩热成熟度指标的均值大于等于0.7％小于1％时获得0.5分，当烃源岩热成熟度指标的均值大于等于1％小于1.1％时随机获得0.5分和0.8分中的一个分值，当烃源岩热成熟度指标的均值大于等于1.1％小于1.4％时获得0.8分，当有烃源岩热成熟度指标的均值大于等于1.4％小于1.5％时随机获得0.8分和1分中的一个分值，当烃源岩热成熟度指标的均值大于等于1.5％时获得1分；
未知单元格的热成熟度评分值评分标准为：当烃源岩热成熟度指标插值小于0.6％时获得0分，当烃源岩热成熟度指标插值小于0.7％大于等于0.6％时随机获得0分和0.5分中的一个分值，当烃源岩热成熟度指标插值大于等于0.7％小于1％时获得0.5分，当烃源岩热成熟度指标插值大于等于1％小于1.1％时随机获得0.5分和0.8分中的一个分值，当烃源岩热成熟度指标插值大于等于1.1％小于1.4％时获得0.8分，当有烃源岩热成熟度指标插值大于等于1.4％小于1.5％时随机获得0.8分和1分中的一个分值，当烃源岩热成熟度指标插值大于等于1.5％时获得1分；
已知单元格的厚度评分值评分标准为：当页岩厚度的均值小于5m时获得0分，当页岩厚度的均值小于10m大于等于5m时随机获得0分和0.5分中的一个分值，当页岩厚度的均值大于等于10m小于15m时获得0.5分，当页岩厚度的均值大于等于15m小于20m时随机获得0.5分和0.8分中的一个分值，当页岩厚度的均值大于等于20m小于25m时获得0.8分，当有页岩厚度的均值大于等于30m小于40m时随机获得0.8分和1分中的一个分值，当页岩厚度的均值大于等于40m时获得1分；
未知单元格的厚度评分值评分标准为：当页岩厚度插值小于5m时获得0分，当页岩厚度插值小于10m大于等于5m时随机获得0分和0.5分中的一个分值，当页岩厚度插值大于等于10m小于15m时获得0.5分，当页岩厚度插值大于等于15m小于20m时随机获得0.5分和0.8分中的一个分值，当页岩厚度插值大于等于20m小于25m时获得0.8分，当有页岩厚度插值大于等于30m小于40m时随机获得0.8分和1分中的一个分值，当页岩厚度插值大于等于40m时获得1分；
已知单元格的经济可采性评分标准为：当埋深的均值小于1000m时获得0分，当埋深的均值小于1500m大于等于1000m时随机获得0分和1分中的一个分值，当埋深的均值大于等于1500m小于3000m时获得0.5分，当埋深的均值大于等于3000m小于3500m时随机获得1分和0.5分中的一个分值，当埋深的均值大于等于3500m小于4000m时获得0.5分，当有埋深的均值大于等于4000m小于4500m时随机获得0.5分和0.2分中的一个分值，当埋深的均值大于等于4500m时获得0.2分；
未知单元格的经济可采性评分标准为：当埋深插值小于1000m时获得0分，当埋深插值小于1500m大于等于1000m时随机获得0分和1分中的一个分值，当埋深插值大于等于1500m小于3000m时获得0.5分，当埋深插值大于等于3000m小于3500m时随机获得1分和0.5分中的一个分值，当埋深插值大于等于3500m小于4000m时获得0.5分，当有埋深插值大于等于4000m小于4500m时随机获得0.5分和0.2分中的一个分值，当埋深插值大于等于4500m时获得0.2分。
在一个具体的实施例中，在步骤S800中，
采用算式(2)计算有利区内每个已知单元格对应区域中的天然气地质资源量，采用算式(3)计算有利区内每个未知单元格对应区域中的天然气地质资源量；
Q1＝Pg×Pe×S×H×ρ×G÷100 (2)
其中，
Q1：已知单元格对应区域内的天然气地质资源量，m3；
Pg：已知单元格的页岩气富集概率，无量纲；
Pe：已知单元格的经济可采性评分，无量纲；
S：已知单元格的面积，km2；
H：已知单元格的页岩厚度的均值，m；
ρ：平均页岩密度，t/m3；
G：已知单元格的页岩含气量的均值，m3/t；
Q2＝Pg×Pe×S×H×ρ×G÷100 (3)
其中，
Q2：未知单元格对应区域内的天然气地质资源量，m3；
Pg：未知单元格的页岩气富集概率，无量纲；
Pe：未知单元格的经济可采性评分，无量纲；
S：未知单元格的面积，km2；
H：未知单元格的页岩厚度插值，m；
ρ：平均页岩密度，t/m3；
G：未知单元格的页岩含气量插值，m3/t。
在一个具体的实施例中，在步骤S900中，预设次数大于或等于1000。
本发明记载的方法利用页岩气富集主控因素等地质资料，按照页岩气富集条件评价和经济可采分级评价标准，优选出页岩气资源禀赋较高的有利区块以及具有商业开采价值的页岩气储层。在此基础上针对优选的富集区分别对页岩厚度、含气性等参数进行不确定性评价，从而准确预测出页岩气资源量。
附图说明
在下文中将基于实施例并参考附图来对本发明进行更详细的描述。其中：
图1为本发明的一个实施例中有利区的分布图；
图2为本发明的一个实施例中的超概率曲线分布图。
在附图中，相同的部件使用相同的附图标记。附图并未按照实际的比例绘制。
具体实施方式
下面将结合附图对本发明作进一步说明。
本实施例中的一种用于预测页岩气资源的方法包括以下步骤：
步骤S100：对待评价区域内的多个处进行测井以获得每个井位处目的储层的有机碳含量、烃源岩热成熟度指标、页岩厚度、页岩密度、页岩含气量、埋深；
对所有页岩密度求平均以获得平均页岩密度。
通过测井来获取每个井位处的有机碳含量、烃源岩热成熟度指标、页岩厚度、页岩密度、页岩含气量、埋深为本领域所公知的技术，在此不再赘述。
步骤S200：对待评价区域进行均匀的网格划分，待评价区域内设置有井位的单元格为已知单元格，其他为未知单元格。
步骤S300：将每个已知单元格中的一个井位处的有机碳含量、烃源岩热成熟度指标、页岩厚度、页岩含气量、埋深的值依次分别作为目的储层在该已知单元格范围内有机碳含量、烃源岩热成熟度指标、页岩厚度、页岩含气量、埋深的均值。
步骤S400：随机构建模拟所有未知单元格的地质属性的顺序，按顺序依次对每个未知单元格的地质属性进行模拟，
模拟第一个未知单元格的地质属性则进行步骤S410，
模拟第二个未知单元格至最后一个未知单元格的地质属性时，每对一个未知单元格的地质属性进行模拟则进行步骤S420。
步骤S410：根据每个已知单元格的有机碳含量的均值、烃源岩热成熟度指标的均值、页岩厚度的均值、页岩含气量的均值、埋深的均值、各个单元格的位置和当前待计算的未知单元格的位置采用克里金插值法计算出该未知单元格中有机碳含量插值、烃源岩热成熟度指标插值、页岩厚度插值、页岩含气量插值、埋深插值。
步骤S420：根据每个已知单元格的有机碳含量的均值、烃源岩热成熟度指标的均值、页岩厚度的均值、页岩含气量的均值、埋深的均值，已经计算出的未知单元格中有机碳含量插值、烃源岩热成熟度指标插值、页岩厚度插值、页岩含气量插值、埋深插值，以及各个单元格的位置和当前待计算的未知单元格的位置采用克里金插值法计算出该未知单元格中有机碳含量插值、烃源岩热成熟度指标插值、页岩厚度插值、页岩含气量插值、埋深插值。
步骤S500：对每个已知单元格的有机碳含量的均值进行评价获得有机碳评分值，当有机碳含量的均值小于0.4％时获得0分，当有机碳含量的均值小于0.5％大于等于0.4％时随机获得0分和0.5分中的一个分值，当有机碳含量的均值大于等于0.5％小于0.9％时获得0.5分，当有机碳含量的均值大于等于0.9％小于1％时随机获得0.5分和0.8分中的一个分值，当有机碳含量的均值大于等于1％小于1.4％时获得0.8分，当有机碳含量的均值大于等于1.4％小于1.5％时随机获得0.8分和1分中的一个分值，当有机碳含量的均值大于等于1.5％时获得1分；
对每个未知单元格的有机碳含量插值进行评价获得有机碳评分值，当有机碳含量插值小于0.4％时获得0分，当有机碳含量插值小于0.5％大于等于0.4％时随机获得0分和0.5分中的一个分值，当机碳含量插值大于等于0.5％小于0.9％时获得0.5分，当有机碳含量插值大于等于0.9％小于1％时随机获得0.5分和0.8分中的一个分值，当机碳含量插值大于等于1％小于1.4％时获得0.8分，当机碳含量插值大于等于1.4％小于1.5％时随机获得0.8分和1分中的一个分值，当有机碳含量插值大于等于1.5％时获得1分；
对每个已知单元格的烃源岩热成熟度指标的均值进行评价获得热成熟度评分值，当烃源岩热成熟度指标的均值小于0.6％时获得0分，当烃源岩热成熟度指标的均值小于0.7％大于等于0.6％时随机获得0分和0.5分中的一个分值，当烃源岩热成熟度指标的均值大于等于0.7％小于1％时获得0.5分，当烃源岩热成熟度指标的均值大于等于1％小于1.1％时随机获得0.5分和0.8分中的一个分值，当烃源岩热成熟度指标的均值大于等于1.1％小于1.4％时获得0.8分，当有烃源岩热成熟度指标的均值大于等于1.4％小于1.5％时随机获得0.8分和1分中的一个分值，当烃源岩热成熟度指标的均值大于等于1.5％时获得1分；
对每个未知单元格的烃源岩热成熟度指标插值进行评价获得热成熟度评分值，当烃源岩热成熟度指标插值小于0.6％时获得0分，当烃源岩热成熟度指标插值小于0.7％大于等于0.6％时随机获得0分和0.5分中的一个分值，当烃源岩热成熟度指标插值大于等于0.7％小于1％时获得0.5分，当烃源岩热成熟度指标插值大于等于1％小于1.1％时随机获得0.5分和0.8分中的一个分值，当烃源岩热成熟度指标插值大于等于1.1％小于1.4％时获得0.8分，当有烃源岩热成熟度指标插值大于等于1.4％小于1.5％时随机获得0.8分和1分中的一个分值，当烃源岩热成熟度指标插值大于等于1.5％时获得1分；
对每个已知单元格的页岩厚度的均值进行评价获得厚度评分值，当页岩厚度的均值小于5m时获得0分，当页岩厚度的均值小于10m大于等于5m时随机获得0分和0.5分中的一个分值，当页岩厚度的均值大于等于10m小于15m时获得0.5分，当页岩厚度的均值大于等于15m小于20m时随机获得0.5分和0.8分中的一个分值，当页岩厚度的均值大于等于20m小于25m时获得0.8分，当有页岩厚度的均值大于等于30m小于40m时随机获得0.8分和1分中的一个分值，当页岩厚度的均值大于等于40m时获得1分；
对每个未知单元格的页岩厚度插值进行评价获得厚度评分值，当页岩厚度插值小于5m时获得0分，当页岩厚度插值小于10m大于等于5m时随机获得0分和0.5分中的一个分值，当页岩厚度插值大于等于10m小于15m时获得0.5分，当页岩厚度插值大于等于15m小于20m时随机获得0.5分和0.8分中的一个分值，当页岩厚度插值大于等于20m小于25m时获得0.8分，当有页岩厚度插值大于等于30m小于40m时随机获得0.8分和1分中的一个分值，当页岩厚度插值大于等于40m时获得1分；
对每个已知单元格的埋深的均值进行评价获得经济可采性评分，当埋深的均值小于1000m时获得0分，当埋深的均值小于1500m大于等于1000m时随机获得0分和1分中的一个分值，当埋深的均值大于等于1500m小于3000m时获得0.5分，当埋深的均值大于等于3000m小于3500m时随机获得1分和0.5分中的一个分值，当埋深的均值大于等于3500m小于4000m时获得0.5分，当有埋深的均值大于等于4000m小于4500m时随机获得0.5分和0.2分中的一个分值，当埋深的均值大于等于4500m时获得0.2分；
对每个未知单元格的埋深插值进行评价获得经济可采性评分，当埋深插值小于1000m时获得0分，当埋深插值小于1500m大于等于1000m时随机获得0分和1分中的一个分值，当埋深插值大于等于1500m小于3000m时获得0.5分，当埋深插值大于等于3000m小于3500m时随机获得1分和0.5分中的一个分值，当埋深插值大于等于3500m小于4000m时获得0.5分，当有埋深插值大于等于4000m小于4500m时随机获得0.5分和0.2分中的一个分值，当埋深插值大于等于4500m时获得0.2分。
步骤S600：将每个已知单元格的有机碳评分值、热成熟度评分值、厚度评分值相乘获得每个已知单元格的页岩气富集概率；
将每个未知单元格的有机碳评分值、热成熟度评分值、厚度评分值相乘获得每个未知单元格的页岩气富集概率。
步骤S700：将有机碳评分值、热成熟度评分值、厚度评分值、经济可采性评分均大于或等于0.5分的单元格的集合划为有利区(如图1所示)。
上述单元格可以是已知单元格也可以是未知单元格，只要单元格所对应的有机碳评分值、热成熟度评分值、厚度评分值、经济可采性评分均大于0.5即可被划到有利区内。
步骤S800：采用算式(1)计算有利区内每个已知单元格对应区域中的天然气地质资源量，采用算式(2)计算有利区内每个未知单元格对应区域中的天然气地质资源量；
Q1＝Pg×Pe×S×H×ρ×G÷100 (1)
其中，
Q1：已知单元格对应区域内的天然气地质资源量，m3；
Pg：已知单元格的页岩气富集概率，无量纲；
Pe：已知单元格的经济可采性评分，无量纲；
S：已知单元格的面积，km2；
H：已知单元格的页岩厚度的均值，m；
ρ：平均页岩密度，t/m3；
G：已知单元格的页岩含气量的均值，m3/t；
Q2＝Pg×Pe×S×H×ρ×G÷100 (2)
其中，
Q2：未知单元格对应区域内的天然气地质资源量，m3；
Pg：未知单元格的页岩气富集概率，无量纲；
Pe：未知单元格的经济可采性评分，无量纲；
S：未知单元格的面积，km2；
H：未知单元格的页岩厚度插值，m；
ρ：平均页岩密度，t/m3；
G：未知单元格的页岩含气量插值，m3/t。
将有利区内每个已知单元格对应区域中的天然气地质资源量与有利区内每个未知单元格对应区域中的天然气地质资源量进行累加获得累积天然气地质资源量。
步骤S900：重复步骤S400至步骤S800预设次数，以获得预设个数的累积天然气地质资源量。
优选地，预设次数大于或等于1000。
步骤S1000：以累积天然气地质资源量数值从大到小的顺序对所有累积天然气地质资源量从序号1开始进行连续编号，按照算式(3)计算每个序号对应的超概率，
超概率＝(序号/N+(序号-1)/N)/2 (3)
其中，N是累积天然气地质资源量的个数，
将每个累积天然气地质资源量及其对应的超概率以坐标点的形式标记在以横轴表示累积天然气地质资源量数值、纵轴表示超概率的直角坐标系上，将直角坐标系上相邻的两个点用直线段连接起来得到如图2所示的超概率曲线，将超概率曲线上超概率等于百分之五十所对应的累积天然气地质资源量(即图2中的点P50的纵坐标的数值)作为待评价区域的天然气地质资源量的最终估值。
虽然已经参考优选实施例对本发明进行了描述，但在不脱离本发明的范围的情况下，可以对其进行各种改进并且可以用等效物替换其中的部件。尤其是，只要不存在结构冲突，各个实施例中所提到的各项技术特征均可以任意方式组合起来。本发明并不局限于文中公开的特定实施例，而是包括落入权利要求的范围内的所有技术方案。

Claims (4)

1.一种用于预测页岩气资源的方法，包括：
步骤S100：对待评价区域内的多个处进行测井以获得每个井位处目的储层的有机碳含量、烃源岩热成熟度指标、页岩厚度、页岩密度、页岩含气量、埋深；
对所有页岩密度求平均以获得平均页岩密度；
步骤S200：对待评价区域进行均匀的网格划分，待评价区域内设置有井位的单元格为已知单元格，其他为未知单元格；
步骤S300：将每个已知单元格中的一个井位处的有机碳含量、烃源岩热成熟度指标、页岩厚度、页岩含气量、埋深的值依次分别作为目的储层在该已知单元格范围内有机碳含量、烃源岩热成熟度指标、页岩厚度、页岩含气量、埋深的均值；
步骤S400：随机构建模拟所有未知单元格的地质属性的顺序，按顺序依次对每个未知单元格的地质属性进行模拟，
模拟第一个未知单元格的地质属性则进行步骤S410，
模拟第二个未知单元格至最后一个未知单元格的地质属性时，每对一个未知单元格的地质属性进行模拟则进行步骤S420；
步骤S410：根据每个已知单元格的有机碳含量的均值、烃源岩热成熟度指标的均值、页岩厚度的均值、页岩含气量的均值、埋深的均值、各个单元格的位置和当前待计算的未知单元格的位置采用克里金插值法计算出该未知单元格中有机碳含量插值、烃源岩热成熟度指标插值、页岩厚度插值、页岩含气量插值、埋深插值；
步骤S420：根据每个已知单元格的有机碳含量的均值、烃源岩热成熟度指标的均值、页岩厚度的均值、页岩含气量的均值、埋深的均值，已经计算出的未知单元格中有机碳含量插值、烃源岩热成熟度指标插值、页岩厚度插值、页岩含气量插值、埋深插值，以及各个单元格的位置和当前待计算的未知单元格的位置采用克里金插值法计算出该未知单元格中有机碳含量插值、烃源岩热成熟度指标插值、页岩厚度插值、页岩含气量插值、埋深插值；
步骤S500：对每个已知单元格的有机碳含量的均值进行评价获得有机碳评分值，对每个未知单元格的有机碳含量插值进行评价获得有机碳评分值，对每个已知单元格的烃源岩热成熟度指标的均值进行评价获得热成熟度评分值，对每个未知单元格的烃源岩热成熟度指标插值进行评价获得热成熟度评分值，对每个已知单元格的页岩厚度的均值进行评价获得厚度评分值，对每个未知单元格的页岩厚度插值进行评价获得厚度评分值，对每个已知单元格的埋深的均值进行评价获得经济可采性评分，对每个未知单元格的埋深插值进行评价获得经济可采性评分；
步骤S600：将每个已知单元格的有机碳评分值、热成熟度评分值、厚度评分值相乘获得每个已知单元格的页岩气富集概率；
将每个未知单元格的有机碳评分值、热成熟度评分值、厚度评分值相乘获得每个未知单元格的页岩气富集概率；
步骤S700：将有机碳评分值、热成熟度评分值、厚度评分值、经济可采性评分均大于或等于0.5分的单元格的集合划为有利区；
步骤S800：计算出有利区内每个已知单元格对应区域中的天然气地质资源量和有利区内每个未知单元格对应区域中的天然气地质资源量，
将有利区内每个已知单元格对应区域中的天然气地质资源量与有利区内每个未知单元格对应区域中的天然气地质资源量进行累加获得累积天然气地质资源量；
步骤S900：重复步骤S400至步骤S800预设次数，以获得预设个数的累积天然气地质资源量；
步骤S1000：以累积天然气地质资源量数值从大到小的顺序对所有累积天然气地质资源量从序号1开始进行连续编号，按照算式(1)计算每个序号对应的超概率，
超概率＝(序号/N+(序号-1)/N)/2 (1)
其中，N是累积天然气地质资源量的个数，
将每个累积天然气地质资源量及其对应的超概率以坐标点的形式标记在以横轴表示累积天然气地质资源量数值、纵轴表示超概率的直角坐标系上，将直角坐标系上相邻的两个点用直线段连接起来得到超概率曲线，将超概率曲线上超概率等于百分之五十所对应的累积天然气地质资源量作为待评价区域的天然气地质资源量的最终估值。
2.根据权利要求1所述的方法，其特征在于，在步骤S500中，
已知单元格的有机碳评分值评分标准为：当有机碳含量的均值小于0.4％时获得0分，当有机碳含量的均值小于0.5％且大于等于0.4％时随机获得0分和0.5分中的一个分值，当有机碳含量的均值大于等于0.5％且小于0.9％时获得0.5分，当有机碳含量的均值大于等于0.9％且小于1％时随机获得0.5分和0.8分中的一个分值，当有机碳含量的均值大于等于1％且小于1.4％时获得0.8分，当有机碳含量的均值大于等于1.4％且小于1.5％时随机获得0.8分和1分中的一个分值，当有机碳含量的均值大于等于1.5％时获得1分；
未知单元格的有机碳评分值评分标准为：当有机碳含量插值小于0.4％时获得0分，当有机碳含量插值小于0.5％且大于等于0.4％时随机获得0分和0.5分中的一个分值，当机碳含量插值大于等于0.5％且小于0.9％时获得0.5分，当有机碳含量插值大于等于0.9％且小于1％时随机获得0.5分和0.8分中的一个分值，当机碳含量插值大于等于1％且小于1.4％时获得0.8分，当机碳含量插值大于等于1.4％且小于1.5％时随机获得0.8分和1分中的一个分值，当有机碳含量插值大于等于1.5％时获得1分；
已知单元格的热成熟度评分值评分标准为：当烃源岩热成熟度指标的均值小于0.6％时获得0分，当烃源岩热成熟度指标的均值小于0.7％且大于等于0.6％时随机获得0分和0.5分中的一个分值，当烃源岩热成熟度指标的均值大于等于0.7％且小于1％时获得0.5分，当烃源岩热成熟度指标的均值大于等于1％且小于1.1％时随机获得0.5分和0.8分中的一个分值，当烃源岩热成熟度指标的均值大于等于1.1％且小于1.4％时获得0.8分，当有烃源岩热成熟度指标的均值大于等于1.4％且小于1.5％时随机获得0.8分和1分中的一个分值，当烃源岩热成熟度指标的均值大于等于1.5％时获得1分；
未知单元格的热成熟度评分值评分标准为：当烃源岩热成熟度指标插值小于0.6％时获得0分，当烃源岩热成熟度指标插值小于0.7％且大于等于0.6％时随机获得0分和0.5分中的一个分值，当烃源岩热成熟度指标插值大于等于0.7％且小于1％时获得0.5分，当烃源岩热成熟度指标插值大于等于1％且小于1.1％时随机获得0.5分和0.8分中的一个分值，当烃源岩热成熟度指标插值大于等于1.1％且小于1.4％时获得0.8分，当有烃源岩热成熟度指标插值大于等于1.4％且小于1.5％时随机获得0.8分和1分中的一个分值，当烃源岩热成熟度指标插值大于等于1.5％时获得1分；
已知单元格的厚度评分值评分标准为：当页岩厚度的均值小于5m时获得0分，当页岩厚度的均值小于10m且大于等于5m时随机获得0分和0.5分中的一个分值，当页岩厚度的均值大于等于10m且小于15m时获得0.5分，当页岩厚度的均值大于等于15m且小于20m时随机获得0.5分和0.8分中的一个分值，当页岩厚度的均值大于等于20m且小于25m时获得0.8分，当有页岩厚度的均值大于等于30m且小于40m时随机获得0.8分和1分中的一个分值，当页岩厚度的均值大于等于40m时获得1分；
未知单元格的厚度评分值评分标准为：当页岩厚度插值小于5m时获得0分，当页岩厚度插值小于10m且大于等于5m时随机获得0分和0.5分中的一个分值，当页岩厚度插值大于等于10m且小于15m时获得0.5分，当页岩厚度插值大于等于15m且小于20m时随机获得0.5分和0.8分中的一个分值，当页岩厚度插值大于等于20m且小于25m时获得0.8分，当有页岩厚度插值大于等于30m且小于40m时随机获得0.8分和1分中的一个分值，当页岩厚度插值大于等于40m时获得1分；
已知单元格的经济可采性评分标准为：当埋深的均值小于1000m时获得0分，当埋深的均值小于1500m且大于等于1000m时随机获得0分和1分中的一个分值，当埋深的均值大于等于1500m且小于3000m时获得0.5分，当埋深的均值大于等于3000m且小于3500m时随机获得1分和0.5分中的一个分值，当埋深的均值大于等于3500m且小于4000m时获得0.5分，当有埋深的均值大于等于4000m且小于4500m时随机获得0.5分和0.2分中的一个分值，当埋深的均值大于等于4500m时获得0.2分；
未知单元格的经济可采性评分标准为：当埋深插值小于1000m时获得0分，当埋深插值小于1500m且大于等于1000m时随机获得0分和1分中的一个分值，当埋深插值大于等于1500m且小于3000m时获得0.5分，当埋深插值大于等于3000m且小于3500m时随机获得1分和0.5分中的一个分值，当埋深插值大于等于3500m且小于4000m时获得0.5分，当有埋深插值大于等于4000m且小于4500m时随机获得0.5分和0.2分中的一个分值，当埋深插值大于等于4500m时获得0.2分。
3.根据权利要求1所述的方法，其特征在于，在步骤S800中，
采用算式(2)计算有利区内每个已知单元格对应区域中的天然气地质资源量，采用算式(3)计算有利区内每个未知单元格对应区域中的天然气地质资源量；
Q1＝Pg×Pe×S×H×ρ×G÷100 (2)
其中，
Q1：已知单元格对应区域内的天然气地质资源量，m3；
Pg：已知单元格的页岩气富集概率，无量纲；
Pe：已知单元格的经济可采性评分，无量纲；
S：已知单元格的面积，km2；
H：已知单元格的页岩厚度的均值，m；
ρ：平均页岩密度，t/m3；
G：已知单元格的页岩含气量的均值，m3/t；
Q2＝Pg×Pe×S×H×ρ×G÷100 (3)
其中，
Q2：未知单元格对应区域内的天然气地质资源量，m3；
Pg：未知单元格的页岩气富集概率，无量纲；
Pe：未知单元格的经济可采性评分，无量纲；
S：未知单元格的面积，km2；
H：未知单元格的页岩厚度插值，m；
ρ：平均页岩密度，t/m3；
G：未知单元格的页岩含气量插值，m3/t。
4.根据权利要求1所述的方法，其特征在于，在步骤S900中，预设次数大于或等于1000。
CN201810289872.2A 2018-03-30 2018-03-30 一种用于预测页岩气资源的方法 Active CN110318744B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810289872.2A CN110318744B (zh) 2018-03-30 2018-03-30 一种用于预测页岩气资源的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810289872.2A CN110318744B (zh) 2018-03-30 2018-03-30 一种用于预测页岩气资源的方法

Publications (2)

Publication Number Publication Date
CN110318744A true CN110318744A (zh) 2019-10-11
CN110318744B CN110318744B (zh) 2022-01-21

Family

ID=68112234

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810289872.2A Active CN110318744B (zh) 2018-03-30 2018-03-30 一种用于预测页岩气资源的方法

Country Status (1)

Country Link
CN (1) CN110318744B (zh)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110930020A (zh) * 2019-11-20 2020-03-27 中国地质大学（北京） 非常规油气资源的经济可采资源量确定方法
CN110991933A (zh) * 2019-12-19 2020-04-10 西南石油大学 一种山地页岩气资源评价方法及系统
CN112282742A (zh) * 2020-10-22 2021-01-29 中国石油大学(华东) 一种页岩油优质储层的预测方法
WO2022012368A1 (zh) * 2020-07-15 2022-01-20 中国石油天然气股份有限公司 一种页岩原位转化可采油气资源量的预测方法及装置
CN117216669A (zh) * 2023-11-09 2023-12-12 中国地质大学（北京） 一种页岩储层分类评价图板的建立方法及应用

Citations (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20130262069A1 (en) * 2012-03-29 2013-10-03 Platte River Associates, Inc. Targeted site selection within shale gas basins
US20130346040A1 (en) * 2012-06-20 2013-12-26 Schlumberger Technology Corporation Shale gas production forecasting
CN104453873A (zh) * 2014-10-23 2015-03-25 中国石油天然气股份有限公司 页岩油气经济有效层段的评价方法
CN104750884A (zh) * 2013-12-26 2015-07-01 中国石油化工股份有限公司 基于多因素非线性回归的页岩油气富集指数定量评价方法
CN104866940A (zh) * 2014-02-26 2015-08-26 中国石油化工股份有限公司 一种地质风险评价随机模拟方法
CN105134195A (zh) * 2015-09-02 2015-12-09 中国石油天然气股份有限公司 一种基于测井资料的页岩气储层品质评价方法
CN105138848A (zh) * 2015-09-07 2015-12-09 中国地质大学（北京） 一种页岩气概率选区方法
CN105528522A (zh) * 2015-12-15 2016-04-27 中国石油天然气股份有限公司 一种基于空间网格计算连续性油气藏资源量的方法及装置
CN105697003A (zh) * 2016-03-17 2016-06-22 成都创源油气技术开发有限公司 页岩气储量评价方法
CN105804730A (zh) * 2015-01-19 2016-07-27 国际商业机器公司 使用历史井数据的资源识别的方法和系统
CN106951686A (zh) * 2017-02-28 2017-07-14 中国石油大学(北京) 页岩气选区评价方法及装置
CN107153755A (zh) * 2016-03-03 2017-09-12 中国石油化工股份有限公司 一种页岩气井数值模拟的求解方法
CN107589466A (zh) * 2016-07-06 2018-01-16 中国石油化工股份有限公司 一种用于评价待发现资源量空间分布特征的方法
CN107842359A (zh) * 2016-09-19 2018-03-27 中国石油化工股份有限公司 单井可采储量评价方法及装置

Patent Citations (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20130262069A1 (en) * 2012-03-29 2013-10-03 Platte River Associates, Inc. Targeted site selection within shale gas basins
US20130346040A1 (en) * 2012-06-20 2013-12-26 Schlumberger Technology Corporation Shale gas production forecasting
CN104750884A (zh) * 2013-12-26 2015-07-01 中国石油化工股份有限公司 基于多因素非线性回归的页岩油气富集指数定量评价方法
CN104866940A (zh) * 2014-02-26 2015-08-26 中国石油化工股份有限公司 一种地质风险评价随机模拟方法
CN104453873A (zh) * 2014-10-23 2015-03-25 中国石油天然气股份有限公司 页岩油气经济有效层段的评价方法
CN105804730A (zh) * 2015-01-19 2016-07-27 国际商业机器公司 使用历史井数据的资源识别的方法和系统
CN105134195A (zh) * 2015-09-02 2015-12-09 中国石油天然气股份有限公司 一种基于测井资料的页岩气储层品质评价方法
CN105138848A (zh) * 2015-09-07 2015-12-09 中国地质大学（北京） 一种页岩气概率选区方法
CN105528522A (zh) * 2015-12-15 2016-04-27 中国石油天然气股份有限公司 一种基于空间网格计算连续性油气藏资源量的方法及装置
CN107153755A (zh) * 2016-03-03 2017-09-12 中国石油化工股份有限公司 一种页岩气井数值模拟的求解方法
CN105697003A (zh) * 2016-03-17 2016-06-22 成都创源油气技术开发有限公司 页岩气储量评价方法
CN107589466A (zh) * 2016-07-06 2018-01-16 中国石油化工股份有限公司 一种用于评价待发现资源量空间分布特征的方法
CN107842359A (zh) * 2016-09-19 2018-03-27 中国石油化工股份有限公司 单井可采储量评价方法及装置
CN106951686A (zh) * 2017-02-28 2017-07-14 中国石油大学(北京) 页岩气选区评价方法及装置

Non-Patent Citations (22)

* Cited by examiner, † Cited by third party
Title
JAMES W.SCHMOKER: ""resource-assessment prespectives for unconventional gas systems"", 《AAPG》 *
NINA E.TRICHE,ETC: ""shales gas volumetrics of unconventional resource plays in canning basin, western australia"", 《SOCIETY OF PETROLEUM ENGINEERS》 *
YUANYIN ZHANG, ZHIJUN JIN, YEQUAN CHEN, XIWU LIU: ""Prestack density inversion based on integrated norm regularization in shale reservoir"", 《SOCIETY OF PETROLEUM ENGINEERS》 *
刘欢，等: ""页岩气地质选区及资源潜力评价方法"", 《西南石油大学学报（自然科学版）》 *
刘超英: ""页岩气勘探选区评价方法探讨"", 《石油实验地质》 *
吴博，等: ""圈闭资源量计算中体积法的软件实现"", 《江汉石油职工大学学报》 *
孙金凤，等: ""页岩气区块优选的模糊动态多属性决策模型研究"", 《运筹与管理》 *
朱华，等: ""页岩气资源评价方法体系及其应用——以川西坳陷为例"", 《天然气工业》 *
柳波，等: ""松辽盆地北部青山口组页岩油形成地质条件及勘探潜力"", 《石油与天然气地质》 *
步少峰，等: ""湘中地区海相页岩气资源潜力评价方法及参数选取"", 《成都理工大学学报（自然科学版）》 *
游声刚: ""重庆地区页岩气有利区优选及开发策略"", 《中国博士学位论文全文数据库 (工程科技Ⅰ辑)》 *
王文广，等: ""东濮凹陷北部沙河街组页岩油有利区优选"", 《大庆石油地质与开发》 *
甘辉: ""长宁地区龙马溪组页岩气资源潜力分析"", 《中国优秀硕士学位论文全文数据库基础科学辑》 *
盛秀杰，等: ""不同地质条件约束下的油气资源量概率组合加和方法及应用"", 《地球科学》 *
盛秀杰，等: ""融合不同地质场景假设的油气区带及圈闭资源量计算方法"", 《天然气地球科学》 *
石石，等: ""基于概率法的油气储量不确定性分析"", 《天然气勘探与开发》 *
谷宁，等: ""基于概率分析的非常规油气选区方法"", 《天然气与石油》 *
贾成业，等: ""概率法在油气储量计算中的应用"", 《天然气工业》 *
邹才能，等: ""常规-非常规油气"有序聚集"理论认识及实践意义"", 《石油勘探与开发》 *
郭秀英，等: ""页岩气选区评价指标筛选及其权重确定方法——以四川盆地海相页岩为例"", 《地质勘探》 *
郭秋麟，等: ""致密油资源评价方法、软件及关键技术"", 《天然气地球科学》 *
金之钧，等: ""油气资源评价方法的基本原则"", 《石油学报》 *

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110930020A (zh) * 2019-11-20 2020-03-27 中国地质大学（北京） 非常规油气资源的经济可采资源量确定方法
CN110930020B (zh) * 2019-11-20 2022-02-11 中国地质大学（北京） 非常规油气资源的经济可采资源量确定方法
CN110991933A (zh) * 2019-12-19 2020-04-10 西南石油大学 一种山地页岩气资源评价方法及系统
WO2022012368A1 (zh) * 2020-07-15 2022-01-20 中国石油天然气股份有限公司 一种页岩原位转化可采油气资源量的预测方法及装置
CN112282742A (zh) * 2020-10-22 2021-01-29 中国石油大学(华东) 一种页岩油优质储层的预测方法
CN117216669A (zh) * 2023-11-09 2023-12-12 中国地质大学（北京） 一种页岩储层分类评价图板的建立方法及应用
CN117216669B (zh) * 2023-11-09 2024-03-22 贵州能源产业研究院有限公司 一种页岩储层分类评价图板的建立方法及应用

Also Published As

Publication number Publication date
CN110318744B (zh) 2022-01-21

Similar Documents

Publication Publication Date Title
CN110318744A (zh) 一种用于预测页岩气资源的方法
CN105719063B (zh) 一种页岩气储量品质分类的综合评价方法
CN104040377A (zh) 用于碳酸盐中的岩石物理学岩石定型的集成工作流程或方法
Milad et al. Impact of lithofacies variations and structural changes on natural fracture distributions
CN111695303B (zh) 一种煤层顶板砂岩含水层充水强度评价方法
CN107092995A (zh) 地热资源评价方法和装置
CN114254960A (zh) 一种低勘探程度区油气资源的评价方法
KR102051225B1 (ko) 셰일층의 파쇄효율 및 가스부존 유망성 예측 방법
CN114638300A (zh) 一种识别页岩油气藏“甜点”的方法、装置及存储介质
CN105549115A (zh) 一种利用录、测、试资料定量描述油气分布的方法
Voegele et al. Microstratigraphic analysis of fossil distribution in the lower Hornerstown and upper Navesink formations at the Edelman Fossil Park, NJ
CN107103377A (zh) 油气区带勘探经济评价方法及装置
Gunnarsson 3D modeling in Petrel of geological CO2 storage site
CN110930020B (zh) 非常规油气资源的经济可采资源量确定方法
CN112147713B (zh) 泥页岩总有机碳含量分段预测方法
CN107060744A (zh) 一种测井地质评价系统
Ali et al. Determination of reservoir hydraulic flow units and permeability estimation using flow zone indicator method
Ungar et al. Fluid Identification from Mud Gas in the Overburden–A Case Study for the Snorre Field
Amanipoor Providing a subsurface reservoir quality maps in oil fields by geostatistical methods
CN112989257B (zh) 海相页岩油气储层的产气量测定方法
CN115422789B (zh) 一种考虑全过程优化断块油藏水驱采收率预测方法及系统
CN110929365A (zh) 一种油砂资源量计算系统
Poursamad et al. Reservoir Quality Evaluation of Sarvak Formation in Gachsaran Oil Field, SW of Iran
Wawrzyniak-Guz et al. Classification of Polish shale gas formations from Baltic Basin, Poland based on well logging data by statistical methods
Pandey et al. Mineral reserve estimation of a limestone mine by inverse square distance method using datamine software

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