CN107657365B - 地质资源开采价值评价方法及装置 - Google Patents

地质资源开采价值评价方法及装置 Download PDF

Info

Publication number
CN107657365B
CN107657365B CN201710804295.1A CN201710804295A CN107657365B CN 107657365 B CN107657365 B CN 107657365B CN 201710804295 A CN201710804295 A CN 201710804295A CN 107657365 B CN107657365 B CN 107657365B
Authority
CN
China
Prior art keywords
geological
value
trap
representing
parameter
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
Application number
CN201710804295.1A
Other languages
English (en)
Other versions
CN107657365A (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.)
Petrochina Co Ltd
Original Assignee
Petrochina Co Ltd
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 Petrochina Co Ltd filed Critical Petrochina Co Ltd
Priority to CN201710804295.1A priority Critical patent/CN107657365B/zh
Publication of CN107657365A publication Critical patent/CN107657365A/zh
Priority to US16/039,823 priority patent/US11037084B2/en
Application granted granted Critical
Publication of CN107657365B publication Critical patent/CN107657365B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

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/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/02—Agriculture; Fishing; Forestry; Mining
    • 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/06315—Needs-based resource requirements planning or analysis

Landscapes

  • Business, Economics & Management (AREA)
  • Human Resources & Organizations (AREA)
  • Engineering & Computer Science (AREA)
  • Strategic Management (AREA)
  • Economics (AREA)
  • Entrepreneurship & Innovation (AREA)
  • General Business, Economics & Management (AREA)
  • Educational Administration (AREA)
  • Development Economics (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Marketing (AREA)
  • Physics & Mathematics (AREA)
  • Tourism & Hospitality (AREA)
  • Health & Medical Sciences (AREA)
  • Agronomy & Crop Science (AREA)
  • Mining & Mineral Resources (AREA)
  • Animal Husbandry (AREA)
  • General Health & Medical Sciences (AREA)
  • Marine Sciences & Fisheries (AREA)
  • Primary Health Care (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Quality & Reliability (AREA)
  • Operations Research (AREA)
  • Game Theory and Decision Science (AREA)
  • Geophysics And Detection Of Objects (AREA)
  • Sampling And Sample Adjustment (AREA)
  • Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)

Abstract

本发明的目的是提供一种地质资源开采价值评价方法及装置，所述方法包括：获取地质构造在多个指定维度上的地质参数；其中，所述地质参数用于表征所述地质构造位于所述指定维度的状态；将所述地质参数分别与每个所述指定维度的阈值进行匹配，得到每个所述指定维度的表征值；多个所述指定维度的所述表征值，形成表征向量；基于所述表征向量进行评价运算，得到所述地质构造中地质资源的评价值。通过本发明的实施方式充分考虑地质资源的阈值对地质资源开采价值的作用，可以对地质资源开采价值进行评价。

Description

地质资源开采价值评价方法及装置
技术领域
本发明涉及地质资源勘探领域，特别涉及一种地质资源开采价值评价方法及装置。
背景技术
地质资源开采包括煤炭、石油、天然气、金属等，但由于地质资源开采的成本比较高，所以对一个区域的开采价值的评价尤其重要。现有技术有的评价方法基于专家经验，有的评价方法对地质资源开采条件中的短板没有统筹考虑，都存在对地质资源开采价值评价准确性不足的问题。例如：在对于圈闭的有效性评价中，所述圈闭有效性评价是地质资源开采价值评价的一种评价。圈闭是油气聚集和存储的主要场所，也是油气勘探和研究的主要对象，而油气勘探中所钻探的圈闭目标是有效圈闭。有效圈闭是指能够获得工业油气流的圈闭。目前，现有的油气圈闭评价技术主要包括：基于圈闭识别和描述的评价技术，侧重于圈闭形态的识别和描述；基于专家知识的圈闭评价技术，主要依靠人为经验对圈闭进行定性综合评价；基于圈闭地质参数赋权重值的评价技术，主要依靠人为对圈闭参数进行赋值半定量评价；圈闭潜在资源量评价技术，主要根据圈闭潜在资源量大小对圈闭进行排队优选和评价；与已知油气藏类比方法，包括多维决策法、人工神经网络、蒙特卡洛、最优化线性理论等圈闭评价技术，主要依据被评价圈闭参数与已知油气藏的比较确定其可钻探性。但是，对圈闭的有效性评价人为比重较大，不能得到特别客观准确的评价。
因此亟需一种地质资源开采价值评价方法获得更好的评价值。
发明内容
本说明书实施方式的目的是提供一种地质资源开采价值评价方法、装置、电子设备及计算机存储介质，能够评价圈闭的有效性。
为实现上述目的，本说明书实施方式提供一种地质资源开采价值评价方法，所述方法包括：获取地质构造在多个指定维度上的地质参数；其中，所述地质参数用于表征所述地质构造位于所述指定维度的状态；将所述地质参数分别与每个所述指定维度的阈值进行匹配，得到每个所述指定维度的表征值；多个所述指定维度的所述表征值，形成表征向量；其中，所述阈值用于表示有价值的地质资源对应的指定维度上的地质参数最小值，所述地质参数最小值对应的地质资源与所述地质资源相同种类且已开采；基于所述表征向量进行评价运算，得到所述地质构造中地质资源的评价值。
为实现上述目的，本说明书实施方式提供一种地质资源开采价值评价装置，所述装置包括：获取单元，获取地质构造在多个指定维度上的地质参数；其中，所述地质参数用于表征所述地质构造位于所述指定维度的状态；匹配单元，将所述地质参数分别与每个所述指定维度的阈值进行匹配，得到每个所述指定维度的表征值；多个所述指定维度的所述表征值，形成表征向量；其中，所述阈值用于表示有价值的地质资源对应的指定维度上的地质参数最小值，所述地质参数最小值对应的地质资源与所述地质资源相同种类且已开采；评价单元，基于所述表征向量进行评价运算，得到所述地质构造中地质资源的评价值。
为实现上述目的，本说明书实施方式提供一种电子设备，包括：输入设备，处理器；输入设备，所述输入设备用于获取所述地质构造在多个指定维度上的地质参数；其中，所述地质参数用于表征所述地质构造位于所述指定维度的状态；处理器，所述处理器用于将所述地质参数分别与每个所述指定维度的阈值进行匹配，得到每个所述指定维度的表征值；多个所述指定维度的所述表征值，形成表征向量；其中，所述阈值用于表示有价值的地质资源对应的指定维度上的地质参数最小值，所述地质参数最小值对应的地质资源与所述地质资源相同种类且已开采；基于所述表征向量进行评价运算，得到所述地质构造中地质资源的评价值。
为实现上述目的，本说明书实施方式提供一种计算机存储介质，所述计算机存储介质存储有计算机程序指令，在所述计算机程序指令被执行时实现：获取所述地质构造在多个指定维度上的地质参数；其中，所述地质参数用于表征所述地质构造位于所述指定维度的状态；将所述地质参数分别与每个所述指定维度的阈值进行匹配，得到每个所述指定维度的表征值；多个所述指定维度的所述表征值，形成表征向量；其中，所述阈值用于表示有价值的地质资源对应的指定维度上的地质参数最小值，所述地质参数最小值对应的地质资源与所述地质资源相同种类且已开采；基于所述表征向量进行评价运算，得到所述地质构造中地质资源的评价值。
由以上本申请实施方式提供的技术方案可见，本申请实施方式通过获取所述地质构造在多个指定维度上的地质参数；将所述地质参数分别与每个所述指定维度的阈值进行匹配，得到每个所述指定维度的表征值；基于所述表征向量进行评价运算，得到所述地质构造中地质资源的评价值。以上做法考虑多个维度，尤其考虑了多个维度的地质参数与其对应的阈值，所述阈值用于表示有价值的地质资源对应的指定维度上的地质参数最小值，所述地质参数最小值对应的地质资源与所述地质资源相同种类且已开采，可以实现对地质资源开采价值的定量评价。
附图说明
为了更清楚地说明本说明书实施方式或现有技术中的技术方案，下面将对实施方式或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本说明书中记载的一些实施方式，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。
图1是本说明书提供的一种地质资源开采价值评价方法的流程图；
图2为本说明书提供的一个圈闭及上覆膏盐岩盖层分布地震剖面特征图；
图3为本说明书提供的一个圈闭特征图；
图4为本说明书提供的一个圈闭等高线地形图；
图5为本说明书提供的一个圈闭膏盐岩盖层厚度分布特征图；
图6为本说明书提供的一个烃源岩岩心镜质体反射率Ro分布特征图；
图7为本说明书提供的一个烃源岩生气丰度分布特征图；
图8为本说明书提供的一个沟通圈闭与烃源岩的断层特征图；
图9为本说明书提供的天然气有效圈闭面积阈值(下限值)特征图；
图10为本说明书提供的工业气流层段储层孔隙度阈值(下限值)特征图；
图11为本说明书提供的天然气有效圈闭的有效储层厚度阈值(下限值)特征图；
图12为本说明书提供的一个目的层孔隙度归一化量化数据特征图；
图13为本说明书提供的一个油气圈闭有效性评价示意图；
图14为本说明书提供的圈闭有效性评价效果图；
图15是本说明书提供的一种电子设备的示意图。
具体实施方式
下面将结合本说明书实施方式中的附图，对本说明书实施方式中的技术方案进行清楚、完整地描述，显然，所描述的实施方式仅仅是本说明书一部分实施方式，而不是全部的实施方式。基于本说明书中的实施方式，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施方式，都应当属于本申请保护的范围。
请参阅图1。本说明书提供的一种地质资源开采价值评价方法的流程图。所述评价方法可以包括以下步骤。
在本实施方式中，执行所述地质资源开采价值评价的客体可以是具有逻辑运算功能的电子设备。所述电子设备可以是服务器和客户端。所述客户端可以为台式电脑、平板电脑、笔记本电脑、工作站等。当然，客户端并不限于上述具有一定实体的电子设备，其还可以为运行于上述电子设备中的软体。还可以是一种通过程序开发形成的程序软件，该程序软件可以运行于上述电子设备中。
在本实施方式中，所述地质资源包括煤、石油、天然气、地热、金属矿产(如铁、锰、铜)、非金属矿产等。所述地质资源处于地质构造中。具体的，例如对于油气的地质资源开采价值评价即圈闭的有效性评价，其他运用此方法进行地质资源开采价值的评价也属于本说明书保护范围。
在本实施方式中，有效圈闭为聚集油气并形成具有工业价值油气藏的圈闭，圈闭有效性定量评价就是评价圈闭是否具有工业价值，有效圈闭的控制因素很多，只有将控制圈闭形成具有工业价值的因素全部考虑才能提高评价符合率。
步骤S10：获取所述地质构造在多个指定维度上的地质参数；其中，所述地质参数用于表征所述地质构造位于所述指定维度的状态。
在本实施方式中，所述地质构造可以是指地质体(如沉积岩体，侵入岩体或矿体等)本身所具有的形态特征，具体的，是指可以储存或运移地质资源的一种地质，包括该地区的结构、本身性质等。在一个具体实施方式中，所述地质构造为圈闭，需要评价该区域的油气开采价值即就是评价所述圈闭的有效性。
在本实施方式中，所述多个指定维度可以是与所述地质资源种类相关的维度。具体的，例如评价煤矿资源时，所述维度可以包括煤层倾角，煤层厚度，煤层范围，煤层深度等多个维度参数；评价石油资源时，所述维度可以包括孔隙度、圈闭面积、盖层厚度、有效储层厚度、输导体系、距生烃中心距离。
在本实施方式中，所述地质参数可以是指所述地质构造中与地质资源开采价值相关的参数；所述地质参数用于表征所述地质构造位于所述指定维度的状态，是指所述地质构造在该所述维度的参数值。具体的，在圈闭有效性评价的场景中，所述地质参数是指表示该维度上状态，包括孔隙度圈闭参数、圈闭面积圈闭参数、盖层厚度圈闭参数、有效储层厚度圈闭参数、输导体系圈闭参数、距生烃中心距离圈闭参数。其中，所述圈闭参数为所述地质参数在圈闭评价时的表述方式。
在本实施方式中，获取所述地质构造在多个指定维度上的地质参数可以通过输入设备的输入，所述输入设备可以是键盘、鼠标、语音输入设备等；也可以通过外接存储设备的输入，所述外接存储设备可以是U盘、机械硬盘等；也可以通过网络接收，例如因特网、局域网；也可以通过读取本地数据等。
步骤S12：将所述地质参数分别与每个所述指定维度的阈值进行匹配，得到每个所述指定维度的表征值；多个所述指定维度的所述表征值，形成表征向量；其中，所述阈值用于表示与所述地质资源相同种类的已开采的有价值的地质资源对应的指定维度上的地质参数最小值。
在本实施方式中，所述表征值是指在所述每一个维度上所述地质参数对地质资源开采价值的贡献值。具体的，在一个圈闭有效性评价的实施例中，所述表征值是指在所述每一个维度上所述圈闭参数对形成工业油气藏的贡献值，用以评价该维度，体现该维度所述阈值的作用和所述地质参数的作用。
在本实施方式中，所述阈值是指具有地质资源开发价值的控因参数的阈值，所述阈值用于表示与所述地质资源相同种类的已开采的有价值的地质资源对应的指定维度上的地质参数最小值。具体的，在圈闭有效性评价的实施方式中，所述阈值是指具有工业价值圈闭控因参数的阈值，当被评价圈闭的某一圈闭参数值低于对应阈值时，该圈闭不能获得工业油气流，只有当圈闭的所有控因参数均在对应阈值之上时，圈闭才为有效圈闭。具体的，所述阈值包括与圈闭参数对应的孔隙度阈值、圈闭面积阈值、盖层厚度阈值、有效储层厚度阈值、输导体系阈值、距生烃中心距离阈值；在该场景中本发明通过考虑至少此6中阈值，克服阈值考虑不足的问题。
在本实施方式中，多个所述指定维度的所述表征值，形成表征向量；所述表征向量是所述表征值的载体，也是所述表征值的集合，用以获得地质资源开采评价值做准备。所述表征向量中的所述表征值具有相同的权重。
在本实施方式中，将所述地质参数分别与每个所述指定维度的阈值进行匹配包括通过所述表征值，将不同维度的数据统一成无量纲参数，不同维度对应的表征值具有相同的权重，归一量化后，便于统筹运算获得圈闭评价值。
步骤S14：基于所述表征向量进行评价运算，得到所述地质构造中地质资源的评价值。
在本实施方式中，所述地质资源开采价值评价值用于评价所述地质资源开采价值，所述地质资源开采价值评价值包括第一特征，所述第一特征包括，在地质资源开采价值评价值处于指定范围情况下，所述地质资源具有开采价值。在本实施方式中，所述圈闭评价值是指反映圈闭有效概率的数值。
在本实施方式中，基于所述表征向量进行评价运算，是指对表征向量中相同权重的表征值进行运算，用以得到一个所述地质构造中地质资源的评价值。所述运算方式可以包括算术平均、几个平均、调和平均、加权平均等算法。当然，所属领域技术人员在本申请技术精髓启示下，还可能采用其它算法得到所述地质构造中地质资源的评价值，但只要其实现的功能和效果与本申请方式相同或相似，均应涵盖于本申请保护范围内。
在一个实施方式中，所述多个维度可以包括孔隙度、圈闭面积、盖层厚度、有效储层厚度、输导体系、距生烃中心距离。所述多个指定维度与圈闭有效性相关；本发明通过考虑至少此6个维度克服有效圈闭考虑的维度不够全面的不足。
在本实施方式中，所述圈闭是一种能阻止油气继续运移并能在其中聚集的场所，所述圈闭面积圈闭参数表示所述圈闭的覆盖面积。所述圈闭面积阈值是指形成具有工业油气藏的最小圈闭面积。
在本实施方式中，获取所述圈闭面积圈闭参数，通过采集研究区及目的层位的测井、地震资料，通过测井标定地震，确定储层发育层位，进行地震解释并成图，确定圈闭面积、闭合高度和溢出点等圈闭参数，请参阅图2、3、4，在一个具体的实施例中，该地质构造中圈闭溢出点海拔-5725m，圈闭高点海拔-5105m，圈闭闭合高度620m，图4中圈闭范围标定的区域就是圈闭面积，获得此实施例的圈闭面积62.4km2。
在本实施方式中，所述盖层为位于储集层之上能封隔储集层使其中流体免于向上逸散的保护层。所述盖层厚度圈闭参数为封堵圈闭全部充满天然气的覆盖层厚度，所述盖层厚度阈值是指盖层可以封堵圈闭全部充满天然气的最小值，圈闭覆盖层封堵条件是有效的。
在本实施方式中，采集研究区及目的层位的测井资料和岩心物性分析资料，通过岩心物性分析数据标定测井资料，通过测井解释确定孔隙度大于形成工业价值气层储层层段孔隙度阈值3.5％的有效储层厚度。在一个具体实施例中，获得该地质构造中圈闭有效储层厚度平均值为94.2m。
在本实施方式中，所述盖层厚度圈闭参数的确定，通过采集目的层位的测井资料、地震资料，通过测井标定地震，确定圈闭上覆盖层顶底界面，请参阅图2，通过地震解释确定圈闭上覆盖层厚度，请参阅图5。在一个具体实施例中，所述地质构造区域为克深1-克深2，所述圈闭上覆盖层厚度大于400m。
在本实施方式中，所述孔隙度石孔隙体积与岩石总体积之比，反映地层储集流体的能力。所述孔隙度阈值是指形成具有工业气流储层的最小孔隙度。
在本实施方式中，所述输导体系是指油气从烃源岩运移到圈闭过程中所经历的所有路径网及其相关围岩，包括连通砂体、断层、不整合及其组合。
在本实施方式中，所述距生烃中心距离圈闭参数是指反映被评价圈位置与油气运移范围关系的参数。所述被评价圈位置与有效生烃范围越近，越能形成有效圈闭。所述距生烃中心距离阈值是指当圈闭处于油气运移最远边界位置时，所述距生烃中心距离圈闭参数的参数值。
在一个场景示例中，获得所述输导体系数据、所述距生烃中心距离维度数据，可以包括以下步骤。
步骤S20：获得目的层的有效生烃范围
在本实施方式中，有效生烃范围是能够形成具有工业价值气藏的区域范围。
在一个场景示例中，采集研究区及目的层位同一油气系统中烃源岩岩心分析总有机碳(TOC，Total Organic Carbon)、镜质体反射率(Ro)分析资料，具体的，可以通过沉积岩中总有机碳的测定(GB/T 19145-2003)、沉积岩中镜质体反射率测定方法(SY/T 5124-2012)进行测定，确定所述地质构造区圈闭天然气来源于侏罗系和三叠系。采集研究区及目的层位的测井资料、地震资料，通过地震解释确定烃源岩厚度，通过测井解释确定TOC分布。请参阅图6、图7，通过盆地模拟确定Ro分布、生气强度分布。随着TOC和Ro的增大，生气强度会增大，根据对已发现工业油气藏的含油气盆地，采集272个已开发气田及所在区域的生气强度资料，当生气强度大于10×108m3/km2时，能够形成具有工业价值气藏，从而根据生气强度大于10×108m3/km2时，确定有效生烃范围。
步骤S22：获得所述输导体系数据。
在本实施方式中，所述输导体系对圈闭油气成藏的贡献采用如下方法确定，当无有效输导体系时为0.1，因此所述输导体系阈值取0.1，当断层或砂体或不整合组合的间接输导体系时为0.5，当砂体为直接输导体系时为0.55，当不整合为输导体系时为0.6，当烃源岩与圈闭储层直接接触为输导体系时为0.8，当断层为直接输导体系时为0.9。请参阅图8，在一个具体实施例中，所述圈闭断层为直接输导体系，其值为0.9。
步骤S24：获得所述距生烃中心距离数据。
在本实施方式中，当被评价圈闭位于有效生烃范围内时的距离参数为0.9，当被评价圈闭位于油气运移最远边界范围处的距离参数为0.1，当被评价圈闭位于油气运移最远边界范围之外时的距离参数为0。具体的，采集已发现油气藏对应圈闭、有效生烃范围、油气运移最远边界位置资料，获取圈闭与有效生烃中心边界的油气运移距离参数，当被评价圈闭位于有效生烃范围到油气运移范围边界时，被评价圈闭与生烃中心的距离参数根据如下函数获得。
Figure BDA0001402320740000071
其中，L表示被评价圈闭的油气运移距离参数，单位km；LS表示被评价圈闭距有效生烃范围边界的距离，单位km；Lmax表示有效生烃中心边界距油气最远运移范围边界距离，单位km。目的层的烃源岩有效疏导体系外边界即为所述油气最远运移范围边界。在一个具体实施例中，所述被评价圈闭的位置，可以是所述圈闭面积的几何中心，当圈闭位于有效生烃中心范围内，此时距生烃中心距离圈闭参数为0.9。
下面结合表1具体说明。
表1某圈闭评价相关参数
参数名称 阈值 圈闭参数
圈闭面积(km<sup>2</sup>) 5 62.4
盖层厚度(m) 0.45 400
孔隙度(％) 3.5 6.6
有效储层厚度(m) 10 94.2
输导体系 0.1 0.9
距生烃中心距离 0.1 0.9
在一个实施方式中，所述圈闭面积阈值为5km2、所述盖层厚度阈值为0.45m、所述孔隙度阈值为3.5％、所述有效储层厚度阈值为10m、所述输导体系阈值为0.1、所述距生烃中心距离阈值为0.1。所述阈值可以运用于其他被评价圈闭，具有通用性。
在本实施方式中，获取所述圈闭面积阈值，通过采集已发现并有效开发的360个圈闭气田，获取的最小圈闭面积为5km2，另外从中随机抽取49个，当单个圈闭天然气可采储量大于20×108m3时为有效圈闭，获取的圈闭面积阈值为5km2，请参阅图9。因此，研究区形成工业价值气藏的圈闭面积阈值采用5km2，此时，所述研究区气柱高度(圈闭闭合高度)为32m。
在本实施方式中，获取所述盖层厚度阈值，可以包括以下步骤。
步骤S30，根据如下函数获取封堵天然气的有效盖层厚度。
Hs＝(g×Hts(ρw-ρgsc/Bgi)/145)/Ts
其中，Hs表示有效盖层厚度，单位m；Hts表示圈闭闭合高度，单位m；ρw表示地层条件水密度，单位g/cm3；ρgsc表示地面条件天然气密度，单位g/cm3；Bgi表示天然气体积系数；Ts表示膏盐岩突破压力梯度差，单位MPa/m；g表示重力加速度，单位m/s2。
在一个具体的实施场景中，圈闭闭合高度为620m；圈闭上覆盖层为膏盐岩，采集被评价圈闭上覆盖层膏盐岩3块岩心样品，开展模拟地层温压条件的突破压力实验，具体的，可以通过岩石气体突破压力测定方法(SY/T 5748-2013)进行测定，3块样品的突破压力梯度差分别为3.5MPa/m、3.8MPa/m、4.2MPa/m，出于保守性考虑，取最小的突破压力梯度差3.5MPa/m作为圈闭上覆盖层的突破压力梯度差。被评价圈闭的地温计算采用该圈闭中部埋藏深度，圈闭高点海拔-5105m，海平面之上地层厚度1671m，计算得到圈闭中部埋藏深度7086m，被评价圈闭所在地区的地温梯度2.188℃/100m，地表平均温度取18℃，计算得到平均地层温度446K；圈闭中部地层压力借用邻区已发现气藏的克拉205井压力梯度1.655MPa/100m，计算得到被评价圈闭中部埋深压力117.30MPa。采集被评价圈闭附近同层位天然气样品，通过PVT实验(压力、体积、温度Pressure-Volume-Temperature)获得天然气偏差系数1.718，具体的，可以按天然气压缩因子的计算标准(GB/T 17747.1-2011、GB/T 17747.2-2011、GB/T 17747.3-2011)，测量获得天然气偏差系数。地面标准温度取291K，地面标准压力取0.101MPa。所述天然气体积系数根据如下函数计算。
Figure BDA0001402320740000091
其中，Bgi表示天然气体积系数，无量纲；Psc表示地面标准压力，单位MPa；Tsc表示地面标准温度，单位K；Zi表示天然气偏差系数，无量纲；Pi表示圈闭中部压力，单位MPa。以此获得天然气体积系数0.00223。采集被评价圈闭地层水，在圈闭中部深度温压条件下，地层水密度为1.067g/cm3。采集被评价圈闭同地区、同层位天然气样品，通过实验测定地面条件下的天然气密度为0.000768g/cm3。以上数据也可以通过采集被评价圈闭或研究区的油气分析资料和突破压力实验资料进行获取。
将以上参数代入函数Hs＝(g×Hts(ρw-ρgsc/Bgi)/145)/Ts中，获得圈闭全部充满天然气时的有效膏盐岩厚度为8.66m，也就是该实施例中有效盖层厚度为8.66m。当膏盐岩盖层厚度大于该值时，盖层可以封堵圈闭全部充满天然气；被评价圈闭上覆膏盐岩厚度大于400m，圈闭上覆盖层封堵条件是有效的。
步骤S32，根据类比原则获取研究区盖层厚度阈值。
具体的，在步骤S30的实施场景中，圈闭闭合高度为620m，对应有效膏盐岩厚度为8.66m，所述研究区形成圈闭面积阈值5km2对应的圈闭中的气柱高度为32m，也是圈闭闭合高度为32m。因此，对应有效封堵32m气柱高度的膏盐岩盖层厚度阈值为0.45m，即形成具有工业价值气藏圈闭的膏盐岩盖层厚度阈值为0.45m。
在另一个场景示例中，所述盖层厚度阈值也可以通过直接对研究区取样测量。具体过程同步骤S30，研究对象直接为研究区，此处只进行简要说明，不再详细复述；获得的盖层厚度阈值与所述类比原则获取研究区盖层厚度阈值差异不大，具体的获得的盖层厚度阈值可以采用研究区为对象的获取结果，也可以采用类比原则获取的结果，也可以是两者平均值。
具体的，根据如下函数，获取以研究区为对象的封堵天然气的盖层厚度阈值。
Hs＝(g×Hts(ρw-ρgsc/Bgi)/145)/Ts
其中，Hs表示有效盖层厚度，单位m；Hts表示研究区圈闭闭合高度，单位m；ρw表示地层条件水密度，单位g/cm3；ρgsc表示地面条件天然气密度，单位g/cm3；Bgi表示天然气体积系数；Ts表示膏盐岩突破压力梯度差，单位MPa/m；g表示重力加速度，单位m/s2。
模拟地层温压条件的突破压力实验，获得膏盐岩突破压力梯度差；采集的天然气样品，通过实验获得天然气偏差系数，根据如下函数计算天然气体积系数。
Figure BDA0001402320740000101
其中，Bgi表示天然气体积系数，无量纲；Psc表示地面标准压力，单位MPa；Tsc表示地面标准温度，单位K；Zi表示天然气偏差系数，无量纲；Pi表示圈闭中部压力，单位MPa。
通过采集研究区圈闭地层水，获得地层水密度；采集研究区圈闭同地区、同层位天然气样品，获得天然气密度。将上述数据带入Hs＝(g×Hts(ρw-ρgsc/Bgi)/145)/Ts，获得研究区有效盖层厚度，即为盖层厚度阈值。
在本实施方式中，采集研究区临近区域、同目的层岩心物性分析资料、试气资料和测井资料，岩心分析孔隙度标定，通过测井资料解释获得工业气流层段储层孔隙度阈值为3.5％，请参阅图10。采集目的层的气层岩心分析孔隙度资料、测井资料，岩心分析孔隙度标定测井，通过测井资料解释，获得被评价圈闭储层孔隙度大于孔隙度阈值3.5％的储层平均孔隙度为6.6％，如图10所示的测井解释(气层)。若没有该项资料，也可以通过采集被评价圈闭目的层段储层岩心样品，按覆压下岩石孔隙和渗透率测定方法(SY/T 6385-1999)测量样品的孔隙度，获取被评价圈闭储层孔隙度大于孔隙度阈值的孔隙度平均值，该值可以采用加权得到。
在本实施方式中，根据360个已开发气田的最小储层厚度，获得形成具有工业气流储层的最小孔隙度对应储层厚度10m，即有效储层厚度阈值为10m，请参阅图11。在一个场景中，圈闭有效储层厚度平均值为94.2m，大于有效储层厚度阈值10m，就被评价圈闭有效储层厚度而言，该圈闭是有效的。
在本实施方式中，根据当无有效输导体系时所述疏导体系为0.1，获得所述输导体系阈值为0.1。
在本实施方式中，当被评价圈闭位于油气运移最远边界范围处的距离参数为0.1，当被评价圈闭位于油气运移最远边界范围之外时的距离参数为0，因此所述距生烃中心距离阈值为0.1。
在一个实施方式中，根据如下函数将所述地质参数分别与每个所述指定维度的阈值进行匹配，得到每个所述指定维度的表征值。
Figure BDA0001402320740000111
其中，WS表示所述表征值，其值范围为0～1；aS为可变系数，与所述地质构造中该关键参数的阈值有关；Smin_P表示所述地质构造中的地质参数与阈值对成藏概率的贡献值，所述Smin_P满足Smin_P＝Slim_p·S/Slim。S表示所述地质构造中对应维度的地质参数值，Slim表示所述地质构造中与表征值对应维度的阈值。具体的，在圈闭有效性评价的场景中，Smin_P表示所述地质构造中的圈闭参数与阈值对成藏概率的贡献值，所述Smin_P满足Smin_P＝Slim_p·S/Slim，S表示对应维度的圈闭参数值，Slim表示与表征值对应维度的阈值。
在一个实施方式中，可以根据如下函数获得所述可变系数aS：
Figure BDA0001402320740000112
其中，Slim_P表示阈值控制地质资源成藏概率值，对应的表征值为WS_50％；Smid_P表示中值控制地质资源成藏概率值，对应的表征值为WS_90％；aS_50％表示在圈闭参数等于阈值时对应的可变系数aS，aS_90％表示在圈闭参数等于参数中值时对应的可变系数aS；所述aS_50％满足aS_50％＝[ln(1-WS_50％)]/Slim_P，所述aS_90％满足aS_90％＝[ln(1-WS_90％)]/Smid_P。具体的，在圈闭有效性评价的实施例中，Slim_P表示阈值控制油气成藏概率值，对应的表征值为WS_50％；Smid_P表示中值控制油气成藏概率值，对应的表征值为WS_90％；aS_50％表示在圈闭参数等于阈值时对应的可变系数aS，aS_90％表示在圈闭参数等于参数中值时对应的可变系数aS；所述aS_50％满足aS_50％＝[ln(1-WS_50％)]/Slim_P，所述aS_90％满足aS_90％＝[ln(1-WS_90％)]/Smid_P。
在本实施方式中，针对圈闭有效性评价的场景。根据对大量已探明油气藏的圈闭参数分析可知，当圈闭参数均为形成工业油气藏对应参数的阈值时，圈闭能够形成工业油气藏的概率约为50％；当圈闭参数均为形成工业油气藏对应参数的中值时，圈闭能够形成工业油气藏的概率约为90％。具体的，一个维度上阈值控制油气成藏概率为50％，即Slim_P＝0.5，当该维度圈闭参数值等于该维度的阈值时对应的表征值为0.5，即WS_50％＝0.5；当圈闭参数的工业油气藏中的中值控制油气成藏概率为90％，即Smid_P＝0.9，该维度圈闭参数等于中值时对应的表征值为0.9，即WS_90％＝0.9；阈值、中值对应的可变系数分别为aS_50％、aS_90％。
下面结合图12具体说明
Figure BDA0001402320740000113
在一个本发明实施例中，孔隙度阈值为3.5％，当孔隙度圈闭参数为获得工业气流层段储层孔隙度阈值3.5％时，对应的表征值为0.5,；当孔隙度为2％时，对应的表征值为0.195；当孔隙度为5％时，对应的表征值为0.763。可见，孔隙度表征值明显放大了获得工业气流层段储层孔隙度阈值附近的差值，提高了敏感度。通过WS方法可以有效放大圈闭成藏参数在阈值附近的精度，大大提高了圈闭评价的准确性。具体的，当圈闭参数处于阈值附近时，所述圈闭参数只要略微波动，所得到的该维度对应的表征值变化幅度较大，即圈闭参数在阈值附近时，表征值的斜率比较大。充分考虑了阈值圈闭有效性的贡献。
在一个实施方式中，根据如下函数进行评价运算，得到所述地质构造中地质资源的评价值。
Figure BDA0001402320740000121
其中，Zs表示地质资源开采价值评价值；Ws_i表示第i个维度的表征值；n表示维度个数。具体的，在圈闭有效性评价的场景中，Zs表示圈闭评价值；Ws_i表示第i个维度的表征值；n表示维度个数。
在本实施方式中，所述评价值Zs其值为0～1，该值越大，指示的地质资源可开采价值的概率越大；反之，则指示的开采价值概率就越小。具体的，当Zs大于或等于0.5时，所述地质资源具有开采价值。具体的，在圈闭有效性评价的场景中，圈闭评价值Zs其值为0～1，该值越大，指示的圈闭有效性概率越大；反之，则指示的圈闭有效性概率就越小。具体的，当Zs大于或等于0.5时，所述圈闭有效。
在本实施方式中，地质资源开采价值评价还可以根据下述方式进行评价，当地质参数小于对应的阈值时，所述地质资源不具有开采价值。
下面结合图13具体说明本说明书一个具体的场景示例，本说明书实施方式提供的地质资源开采价值评价方法，可以通过运行于电子设备中的软件实现。该地质资源开采价值评价方法对某一个地理位置的地质资源进行评价，提供作业依据。
在本场景示例中，所述地质资源可以是油气，所述地质构造可以是圈闭。获取所述地质构造在多个指定维度上的地质参数；例如，在圈闭有效性评价中，获取圈闭面积参数、盖层厚度参数、孔隙度参数、有效储层厚度参数、输导体系参数、距生烃中心距离参数。所述地质参数可以实验由来，也可以是测井资料中获取等；可以通过外接存储设备、网络传输、输入设备输入等方式获得。在该示例中，所述地质资源开采价值评价具体的是一个圈闭有效性评价方法。所述评价方法包括：采集研究区及目的层段测井资料、地震资料，获取储层发育层位、圈闭面积、圈闭闭合高度和圈闭溢出点；采集已有效开发油气田的圈闭面积资料，获取有效圈闭面积阈值(下限值)。采集研究区及目的层段油气分析资料和盖层突破压力实验资料，获取圈闭盖层厚度阈值(下限值)；采集研究区及目的层段测井资料、地震资料，获取圈闭盖层厚度及分布。采集研究区及邻近区域目的层段岩心物性分析资料、试油试气资料、测井资料，获取工业油气层段储层孔隙度阈值(下限值)，获取被评价圈闭储层孔隙度大于孔隙度阈值(下限值)的孔隙度平均值。采集研究区及目的层段测井资料，获取被评价圈闭储层孔隙度大于孔隙度阈值(下限值)的有效储层厚度；采集已有效开发油气田的有效储层厚度资料，获取有效储层厚度阈值(下限值)。采集已有效开发油气田所在区域烃源岩的生气强度、获得形成具有工业价值油气圈闭所在区域的生气强度阈值(下限值)；采集研究区及目的层段的烃源岩岩心总有机碳含量(TOC)分析资料、地球物理资料，获取TOC分布、烃源岩厚度分布；采集研究区及目的层段的烃源岩岩心镜质体反射率(Ro)分析资料，获取Ro分布、生气强度分布、有效生气范围。采集研究区及目的层段的地震资料，获取沟通圈闭与烃源岩的油气输导体系；采集研究区目的层的有效生烃中心位置、有效生烃范围、被评价圈闭位置，获取油气运移距离参数。
在本场景示例中，将所述地质参数分别与每个所述指定维度的阈值进行匹配，得到每个所述指定维度的表征值；其中，所述阈值用于表示有价值的地质资源对应的指定维度上的地质参数最小值，所述地质参数最小值对应的地质资源与所述地质资源相同种类且已开采；所述阈值可以通过实验获得，也可以通过历史资料获得。将所述地质参数与对应的所述阈值匹配包括，在对应维度上的地质参数与阈值进行匹配，通过计算获得该维度的表征值，用以评价该维度，体现该维度所述阈值的作用和所述地质参数的作用。例如，所述地质参数除以所述阈值获得的值也是一种表征值的获取方式，优选的，可以通过函数
Figure BDA0001402320740000131
获得，因为利用该运算规则，所述地质参数在所述阈值附近变化时，函数斜率比较大，可以体现阈值的木桶效应作用。计算方式可以使电子设备中预设的指令等。通过不同维度上所述地质参数与所述阈值的匹配，获得对应的表征值，不同维度的表征值，形成用来表示该地区地质结构的表征向量。具体的，对被评价圈闭的圈闭面积、储层孔隙度、盖层厚度、储层厚度、输导体系和运移距离进行归一量化处理，获取对应参数的表征值，形成表征向量。
在本场景示例中，基于所述表征向量进行评价运算，得到所述地质构造中地质资源的评价值。具体的，可以根据如下函数确定地质资源开采价值评价值：
Figure BDA0001402320740000132
其中，Zs表示地质资源开采价值评价值；Ws_i表示第i个维度的表征值；n表示维度个数。当圈闭评价结果值大于等于0.5时为有效圈闭，当圈闭评价结果值小于0.5时为无效圈闭。下面结合表1，具体说明。所述圈闭面积阈值为5km2、所述盖层厚度阈值为0.45m、所述孔隙度阈值为3.5％、所述有效储层厚度阈值为10m、所述输导体系阈值为0.1、所述距生烃中心距离阈值为0.1。所述圈闭面积圈闭参数为62.4km2、所述盖层厚度圈闭参数为400m、所述孔隙度圈闭参数为6.6％、所述有效储层厚度圈闭参数为94.2m、所述输导体系圈闭参数为0.9、所述距生烃中心距离圈闭参数为0.9。根据所述
Figure BDA0001402320740000141
获得所述圈闭面积的表征值为1，所述盖层厚度的表征值为1，所述孔隙度的表征值为0.920，所述有效储层厚度的表征值为1，所述疏导体系的表征值为1，所述距生烃中心距离的表征值为1。根据函数
Figure BDA0001402320740000142
获得该圈闭评价值Zs约为0.983，Zs大于或等于0.5，所以圈闭有效，且圈闭Zs的值为0～1，Zs越大，该圈闭为圈闭有效概率越大。
同理，当所述地质资源开采价值评价方法可以在煤矿资源开采价值的评价。获取煤矿资源开采价值的多个维度参数，所述参数包括地质参数和阈值。所述维度包括煤层倾角，煤层厚度，煤层范围，煤层深度等多个维度参数，所述评价值的确定方法同圈闭有效性评价值的确定方法此处不再复述。当地质资源为金属矿产时，所述多个维度可以包括矿床大小，高度，开采深度等。所述评价值的确定方法同圈闭有效性评价值的确定方法此处不再复述
为了清楚的说明本说明书实施方式的有益效果，下面结合附图14进行说明。利用本发明对已发现的13个工业气藏圈闭进行评价，其中，8个工业气藏圈闭钻探，未获得工业价值的5个圈闭，圈闭评价结果符合率较好。
本说明书实施方式中还提供了一种地质资源开采价值评价装置，如下面的实施方式所述。由于一种地质资源开采价值评价装置解决问题的原理与一种地质资源开采价值评价方法相似，因此一种地质资源开采价值评价的实施可以参见一种地质资源开采价值评价方法的实施，重复之处不再赘述。以下所使用的，术语“单元”或者“模块”可以实现预定功能的软件和/或硬件的组合。尽管以下实施例所描述的装置较佳地以软件来实现，但是硬件，或者软件和硬件的组合的实现也是可能并被构想的。该装置具体可以包括：获取单元、匹配单元、评价单元。下面对该结构进行具体说明。
获取单元，获取所述地质构造在多个指定维度上的地质参数；其中，所述地质参数用于表征所述地质构造位于所述指定维度的状态。
匹配单元，将所述地质参数分别与每个所述指定维度的阈值进行匹配，得到每个所述指定维度的表征值；多个所述指定维度的所述表征值，形成表征向量；其中，所述阈值用于表示与所述地质资源相同种类的已开采的有价值的地质资源对应的指定维度上的地质参数最小值。
在本实施方式中，所述表征值是指在所述每一个维度上所述地质参数对地质资源开采价值的贡献值。具体的，在一个圈闭有效性评价的实施例中，所述表征值是指在所述每一个维度上所述圈闭参数对形成工业油气藏的贡献值。通过所述表征值，将不同维度的数据统一成无量纲参数，不同维度对应的表征值具有相同的权重，归一量化后，便于统筹运算获得圈闭评价值。
评价单元，所述评价单元用于基于所述表征向量进行评价运算，得到所述地质构造中地质资源的评价值。
在一个实施方式中，所述多个维度包括孔隙度、圈闭面积、盖层厚度、有效储层厚度、输导体系、距生烃中心距离。所述多个维度与圈闭有效性相关；本发明通过考虑至少此6个维度的圈闭参数及阈值，克服有效圈闭考虑的维度不够全面以及阈值考虑不周到的不足。
在本实施方式中，阈值是指具有地质资源开发价值的控因参数的阈值。具体的，在圈闭有效性评价的场景中，所述阈值是指具有工业价值圈闭控因参数的阈值，当被评价圈闭的某一圈闭参数值低于对应阈值时，该圈闭不能获得工业油气流，只有当圈闭的所有控因参数均在对应阈值之上时，圈闭才为有效圈闭。具体的，包括与圈闭参数对应的孔隙度阈值、圈闭面积阈值、盖层厚度阈值、有效储层厚度阈值、输导体系阈值、距生烃中心距离阈值。
在一个实施方式中，根据如下函数将所述地质参数分别与每个所述指定维度的阈值进行匹配，得到每个所述指定维度的表征值。
Figure BDA0001402320740000151
其中，WS表示所述表征值，其值范围为0～1；aS为可变系数，与所述地质构造中该关键参数的阈值有关；Smin_P表示所述地质参数与阈值对成藏概率的贡献值，所述Smin_P满足Smin_P＝Slim_p·S/Slim。S表示对应维度的地质参数值，Slim表示与表征值对应维度的阈值。具体的，在圈闭有效性评价的场景中，Smin_P表示圈闭参数与阈值对成藏概率的贡献值，所述Smin_P满足Smin_P＝Slim_p·S/Slim，S表示对应维度的圈闭参数值，Slim表示与表征值对应维度的阈值。
在一个实施方式中，根据如下函数进行评价运算，得到所述地质构造中地质资源的评价值。
Figure BDA0001402320740000161
其中，Zs表示地质资源开采价值评价值；Ws_i表示第i个维度的表征值；n表示维度个数。具体的，在圈闭有效性评价的场景中，Zs表示圈闭评价值；Ws_i表示第i个维度的表征值；n表示维度个数。
具体的，在一个实施方式中，评价值Zs其值为0～1，该值越大，指示的地质资源可开采价值的概率越大；反之，则指示的开采价值概率就越小。具体的，当Zs大于或等于0.5时，所述地质资源具有开采价值。具体的，在圈闭有效性评价的场景中，圈闭评价值Zs其值为0～1，该值越大，指示的圈闭有效性概率越大；反之，则指示的圈闭有效性概率就越小。具体的，当Zs大于或等于0.5时，所述圈闭有效。
请参阅图15。本说明书实施方式中还提供了一种电子设备，包括：输入设备，处理器；所述输入设备用于获取所述地质构造在多个指定维度上的地质参数；其中，所述地质参数用于表征所述地质构造位于所述指定维度的状态；所述处理器用于将所述地质参数分别与每个所述指定维度的阈值进行匹配，得到每个所述指定维度的表征值；多个所述指定维度的所述表征值，形成表征向量；其中，所述阈值用于表示有价值的地质资源对应的指定维度上的地质参数最小值，所述地质参数最小值对应的地质资源与所述地质资源相同种类且已开采；基于所述表征向量进行评价运算，得到所述地质构造中地质资源的评价值。
所述输入设备是用户和计算机系统之间进行信息交换的主要装置之一。所述输入设备包括键盘、鼠标、摄像头、扫描仪、光笔、手写输入板、语音输入装置等；输入设备用于把原始数据和处理这些数的程序输入到计算机中。所述输入设备还可以获取接收其他模块、单元、设备传输过来的数据。
所述处理器可以按任何适当的方式实现。例如，处理器可以采取例如微处理器或处理器以及存储可由该(微)处理器执行的计算机可读程序代码(例如软件或固件)的计算机可读介质、逻辑门、开关、专用集成电路(Application Specific Integrated Circuit，ASIC)、可编程逻辑控制器和嵌入微控制器的形式等等。
在本实施方式中，该电子设备具体实现的功能和效果，可以与其它实施方式对照解释，在此不再赘述。
本说明书实施方式中还提供了一种计算机存储介质，所述计算机存储介质存储有计算机程序指令，在所述计算机程序指令被执行时实现：获取所述地质构造在多个指定维度上的地质参数；其中，所述地质参数用于表征所述地质构造位于所述指定维度的状态；将所述地质参数分别与每个所述指定维度的阈值进行匹配，得到每个所述指定维度的表征值；多个所述指定维度的所述表征值，形成表征向量；其中，所述阈值用于表示有价值的地质资源对应的指定维度上的地质参数最小值，所述地质参数最小值对应的地质资源与所述地质资源相同种类且已开采；基于所述表征向量进行评价运算，得到所述地质构造中地质资源的评价值。
在本实施方式中，所述存储器包括但不限于随机存取存储器(Random AccessMemory,RAM)、只读存储器(Read-Only Memory,ROM)、缓存(Cache)、硬盘(Hard DiskDrive,HDD)或者存储卡(Memory Card)。所述存储器可以用于存储计算机程序指令。网络通信单元可以是依照通信协议规定的标准设置的，用于进行网络连接通信的接口。
在本实施方式中，该计算机存储介质存储的程序指令具体实现的功能和效果，可以与其它实施方式对照解释，在此不再赘述。
本说明书实施方式提出地质资源开采价值评价最本质的特点是，考虑与该地质资源相关的多个维度，并充分考虑阈值对地质资源开采价值的影响，将所述多个维度的地质参数进行统一量化处理，制定了一种地质资源开采价值定量评价标准，提供了一种地质资源开采价值评价方法，尤其提供了一种圈闭有效性定量评价的新方法。
尽管本申请内容中提到一种地质资源开采价值评价方法、装置、电子设备或计算机存储介质，但是，本申请并不局限于必须是行业标准或实施例所描述的情况等，某些行业标准或者使用自定义方式或实施例描述的实施基础上略加修改后的实施方案也可以实现上述实施例相同、等同或相近、或变形后可预料的实施效果。应用这些修改或变形后的数据获取、处理、输出、判断方式等的实施例，仍然可以属于本申请的可选实施方案范围之内。
虽然本申请提供了如实施例或流程图所述的方法操作步骤，但基于常规或者无创造性的手段可以包括更多或者更少的操作步骤。实施例中列举的步骤顺序仅仅为众多步骤执行顺序中的一种方式，不代表唯一的执行顺序。在实际中的装置或客户端产品执行时，可以按照实施例或者附图所示的方法顺序执行或者并行执行(例如并行处理器或者多线程处理的环境，甚至为分布式数据处理环境)。术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、产品或者设备不仅包括那些要素，而且还包括没有明确列出的其他要素，或者是还包括为这种过程、方法、产品或者设备所固有的要素。在没有更多限制的情况下，并不排除在包括所述要素的过程、方法、产品或者设备中还存在另外的相同或等同要素。
上述实施例阐明的装置或模块等，具体可以由计算机芯片或实体实现，或者由具有某种功能的产品来实现。为了描述的方便，描述以上装置时以功能分为各种模块分别描述。当然，在实施本申请时可以把各模块的功能在同一个或多个软件和/或硬件中实现，也可以将实现同一功能的模块由多个子模块的组合实现等。以上所描述的装置实施例仅仅是示意性的，例如，所述模块的划分，仅仅为一种逻辑功能划分，实际实现时可以有另外的划分方式，例如多个模块或组件可以结合或者可以集成到另一个系统，或一些特征可以忽略，或不执行。
本领域技术人员也知道，除了以纯计算机可读程序代码方式实现控制器以外，完全可以通过将方法步骤进行逻辑编程来使得控制器以逻辑门、开关、专用集成电路、可编程逻辑控制器和嵌入微控制器等的形式来实现相同功能。因此这种控制器可以被认为是一种硬件部件，而对其内部包括的用于实现各种功能的装置也可以视为硬件部件内的结构。或者甚至，可以将用于实现各种功能的装置视为既可以是实现方法的软件模块又可以是硬件部件内的结构。
本申请可以在由计算机执行的计算机可执行指令的一般上下文中描述，例如程序模块。一般地，程序模块包括执行特定任务或实现特定抽象数据类型的例程、程序、对象、组件、数据结构、类等等。也可以在分布式计算环境中实践本申请，在这些分布式计算环境中，由通过通信网络而被连接的远程处理设备来执行任务。在分布式计算环境中，程序模块可以位于包括存储设备在内的本地和远程计算机存储介质中。
通过以上的实施方式的描述可知，本领域的技术人员可以清楚地了解到本申请可借助软件加必需的通用硬件平台的方式来实现。基于这样的理解，本申请的技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来，该计算机软件产品可以存储在存储介质中，如ROM/RAM、磁碟、光盘等，包括若干指令用以使得一台计算机设备(可以是个人计算机，移动终端，服务器，或者网络设备等)执行本申请各个实施例或者实施例的某些部分所述的方法。
本说明书中的各个实施例采用递进的方式描述，各个实施例之间相同或相似的部分互相参见即可，每个实施例重点说明的都是与其他实施例的不同之处。本申请可用于众多通用或专用的计算机系统环境或配置中。例如：个人计算机、服务器计算机、手持设备或便携式设备、平板型设备、多处理器系统、基于微处理器的系统、置顶盒、可编程的电子设备、网络PC、小型计算机、大型计算机、包括以上任何系统或设备的分布式计算环境等等。
虽然通过实施例描绘了本申请，本领域普通技术人员知道，本申请有许多变形和变化而不脱离本申请的精神，希望所附的权利要求包括这些变形和变化而不脱离本申请。

Claims (7)

1.一种地质资源开采价值评价方法，其特征在于，所述方法包括：
获取地质构造在多个指定维度上的地质参数；其中，所述地质参数用于表征所述地质构造位于所述指定维度的状态；所述多个指定维度包括：圈闭面积、盖层厚度、孔隙度、有效储层厚度、输导体系、距生烃中心距离；
将所述地质参数分别与每个所述指定维度的阈值进行匹配，得到每个所述指定维度的表征值，包括：根据如下函数将所述地质参数分别与每个所述指定维度的阈值进行匹配：
Figure FDA0002443753650000011
其中，WS表示所述表征值；aS为可变系数；Smin_P表示任一指定维度上的地质参数值对成藏概率的贡献值，所述Smin_P满足Smin_P＝Slim_p·S/Slim，S表示所述指定维度上的地质参数值，Slim表示所述指定维度上的阈值，Slim_p表示所述阈值控制地质资源成藏概率值；所述可变系数aS根据如下函数获得：
Figure FDA0002443753650000012
其中，Slim_P对应的表征值为WS_50％；Smid_P表示阈值控制成藏概率值，对应的表征值为WS_90％；aS_50％表示在圈闭参数等于阈值时对应的可变系数aS，aS_90％表示在圈闭参数等于参数中值时对应的可变系数aS；所述aS_50％满足aS_50％＝[ln(1-WS_50％)]/Slim_P，所述aS_90％满足aS_90％＝[ln(1-WS_90％)]/Smid_P；多个所述指定维度的所述表征值，形成表征向量；其中，所述阈值用于表示有价值的地质资源对应的指定维度上的地质参数最小值，所述地质参数最小值对应的地质资源与所述地质资源相同种类且已开采；
基于所述表征向量进行评价运算，得到所述地质构造中地质资源的评价值。
2.如权利要求1所述的方法，其特征在于，所述圈闭面积的阈值为5km2、所述盖层厚度的阈值为0.45m、所述孔隙度的阈值为3.5％、所述有效储层厚度的阈值为10m、所述输导体系的阈值为0.1、所述距生烃中心距离的阈值为0.1。
3.如权利要求1所述的方法，其特征在于，根据如下函数获得所述距生烃中心距离的地质参数：
Figure FDA0002443753650000013
其中，L表示被评价圈闭的油气运移距离参数，单位km；LS表示被评价圈闭距有效生烃范围边界的距离，单位km；Lmax表示有效生烃中心边界距油气最远运移范围边界距离，单位km。
4.如权利要求1所述的方法，其特征在于，根据如下函数确定地质资源开采价值评价值：
Figure FDA0002443753650000021
其中，Zs表示地质资源开采价值评价值；Ws_i表示第i个维度的表征值；n表示维度个数。
5.一种地质资源开采价值评价装置，其特征在于，所述装置包括：
获取单元，获取地质构造在多个指定维度上的地质参数；其中，所述地质参数用于表征所述地质构造位于所述指定维度的状态；所述多个指定维度包括：圈闭面积、盖层厚度、孔隙度、有效储层厚度、输导体系、距生烃中心距离；
匹配单元，将所述地质参数分别与每个所述指定维度的阈值进行匹配，得到每个所述指定维度的表征值，包括：根据如下函数将所述地质参数分别与每个所述指定维度的阈值进行匹配：
Figure FDA0002443753650000022
其中，WS表示所述表征值；aS为可变系数；Smin_P表示任一指定维度上的地质参数值对成藏概率的贡献值，所述Smin_P满足Smin_P＝Slim_p·S/Slim，S表示所述指定维度上的地质参数值，Slim表示所述指定维度上的阈值，Slim_p表示所述阈值控制地质资源成藏概率值；所述可变系数aS根据如下函数获得：
Figure FDA0002443753650000023
其中，Slim_P对应的表征值为WS_50％；Smid_P表示阈值控制成藏概率值，对应的表征值为WS_90％；aS_50％表示在圈闭参数等于阈值时对应的可变系数aS，aS_90％表示在圈闭参数等于参数中值时对应的可变系数aS；所述aS_50％满足aS_50％＝[ln(1-WS_50％)]/Slim_P，所述aS_90％满足aS_90％＝[ln(1-WS_90％)]/Smid_P；多个所述指定维度的所述表征值，形成表征向量；其中，所述阈值用于表示与所述地质资源相同种类的已开采的有价值的地质资源对应的指定维度上的地质参数最小值；
评价单元，基于所述表征向量进行评价运算，得到所述地质构造中地质资源的评价值。
6.一种电子设备，其特征在于，包括：输入设备，处理器；
所述输入设备用于获取地质构造在多个指定维度上的地质参数；其中，所述地质参数用于表征所述地质构造位于所述指定维度的状态；所述多个指定维度包括：圈闭面积、盖层厚度、孔隙度、有效储层厚度、输导体系、距生烃中心距离；
所述处理器用于将所述地质参数分别与每个所述指定维度的阈值进行匹配，得到每个所述指定维度的表征值，包括：根据如下函数将所述地质参数分别与每个所述指定维度的阈值进行匹配：
Figure FDA0002443753650000031
其中，WS表示所述表征值；aS为可变系数；Smin_P表示任一指定维度上的地质参数值对成藏概率的贡献值，所述Smin_P满足Smin_P＝Slim_p·S/Slim，S表示所述指定维度上的地质参数值，Slim表示所述指定维度上的阈值，Slim_p表示所述阈值控制地质资源成藏概率值；所述可变系数aS根据如下函数获得：
Figure FDA0002443753650000032
其中，Slim_P对应的表征值为WS_50％；Smid_P表示阈值控制成藏概率值，对应的表征值为WS_90％；aS_50％表示在圈闭参数等于阈值时对应的可变系数aS，aS_90％表示在圈闭参数等于参数中值时对应的可变系数aS；所述aS_50％满足aS_50％＝[ln(1-WS_50％)]/Slim_P，所述aS_90％满足aS_90％＝[ln(1-WS_90％)]/Smid_P；多个所述指定维度的所述表征值，形成表征向量；其中，所述阈值用于表示与所述地质资源相同种类的已开采的有价值的地质资源对应的指定维度上的地质参数最小值；基于所述表征向量进行评价运算，得到所述地质构造中地质资源的评价值。
7.一种计算机存储介质，其特征在于，所述计算机存储介质存储有计算机程序指令，在所述计算机程序指令被执行时实现：
获取地质构造在多个指定维度上的地质参数；其中，所述地质参数用于表征所述地质构造位于所述指定维度的状态；所述多个指定维度包括：圈闭面积、盖层厚度、孔隙度、有效储层厚度、输导体系、距生烃中心距离；
将所述地质参数分别与每个所述指定维度的阈值进行匹配，得到每个所述指定维度的表征值，包括：根据如下函数将所述地质参数分别与每个所述指定维度的阈值进行匹配：
Figure FDA0002443753650000041
其中，WS表示所述表征值；aS为可变系数；Smin_P表示任一指定维度上的地质参数值对成藏概率的贡献值，所述Smin_P满足Smin_P＝Slim_p·S/Slim，S表示所述指定维度上的地质参数值，Slim表示所述指定维度上的阈值，Slim_p表示所述阈值控制地质资源成藏概率值；所述可变系数aS根据如下函数获得：
Figure FDA0002443753650000042
其中，Slim_P对应的表征值为WS_50％；Smid_P表示阈值控制成藏概率值，对应的表征值为WS_90％；aS_50％表示在圈闭参数等于阈值时对应的可变系数aS，aS_90％表示在圈闭参数等于参数中值时对应的可变系数aS；所述aS_50％满足aS_50％＝[ln(1-WS_50％)]/Slim_P，所述aS_90％满足aS_90％＝[ln(1-WS_90％)]/Smid_P；多个所述指定维度的所述表征值，形成表征向量；其中，所述阈值用于表示与所述地质资源相同种类的已开采的有价值的地质资源对应的指定维度上的地质参数最小值；
基于所述表征向量进行评价运算，得到所述地质构造中地质资源的评价值。
CN201710804295.1A 2017-09-08 2017-09-08 地质资源开采价值评价方法及装置 Active CN107657365B (zh)

Priority Applications (2)

Application Number Priority Date Filing Date Title
CN201710804295.1A CN107657365B (zh) 2017-09-08 2017-09-08 地质资源开采价值评价方法及装置
US16/039,823 US11037084B2 (en) 2017-09-08 2018-07-19 Method and apparatus for evaluating exploitation value of a geological resource

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710804295.1A CN107657365B (zh) 2017-09-08 2017-09-08 地质资源开采价值评价方法及装置

Publications (2)

Publication Number Publication Date
CN107657365A CN107657365A (zh) 2018-02-02
CN107657365B true CN107657365B (zh) 2020-08-11

Family

ID=61129285

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710804295.1A Active CN107657365B (zh) 2017-09-08 2017-09-08 地质资源开采价值评价方法及装置

Country Status (2)

Country Link
US (1) US11037084B2 (zh)
CN (1) CN107657365B (zh)

Families Citing this family (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107807407B (zh) * 2017-09-30 2019-10-11 中国石油天然气股份有限公司 一种油气区带有效性评价方法和装置
CN111767624A (zh) * 2019-03-15 2020-10-13 中国石油化工股份有限公司 一种基于地质帕累托参数先验求解的资源评价方法
CN110176070A (zh) * 2019-05-07 2019-08-27 中山大学 一种基于钻孔数据进行机器学习的地层序列模拟方法
CN110470762B (zh) * 2019-08-22 2021-11-02 陕西延长石油(集团)有限责任公司研究院 一种提高芳烃参数定量评价烃源岩成熟度准确度的方法
CN110826936B (zh) * 2019-11-22 2022-09-20 中国地质大学（北京） 一种页岩油气资源分级评价方法
CN113496301B (zh) * 2020-04-01 2024-03-29 中国石油天然气股份有限公司 油气田资产评价方法及装置
EP3926368B1 (en) 2020-06-15 2024-01-10 TotalEnergies OneTech Analyzing a hydrocarbon trap
CN114167515B (zh) * 2020-09-11 2023-03-28 中国石油化工股份有限公司 一种岩性圈闭有效性识别的方法
CN112785448B (zh) * 2020-11-24 2024-03-05 中国石油天然气股份有限公司 一种油气资源区带评价方法
CN112989698A (zh) * 2021-03-11 2021-06-18 上海城勘信息科技有限公司 一种地层分类方法、装置、电子设备及存储介质
CN114047547B (zh) * 2021-11-25 2023-09-22 中国石油大学(华东) 膏盐岩类盖层封闭能力的定量评价方法
CN115099578B (zh) * 2022-06-01 2023-05-26 中国地质调查局油气资源调查中心 一种含油气盆地边缘地区的油气资源评价方法
CN115146976A (zh) * 2022-07-11 2022-10-04 中国地质大学（北京） 影响待勘测区块评价的主要地质参数的选取方法及设备
CN115354992A (zh) * 2022-08-31 2022-11-18 成都理工大学 一种基于岩性组合含气性特征的煤系气储层评价方法
CN116341925B (zh) * 2023-02-23 2024-06-18 成都理工大学 致密砂岩优质储层连续性评价系统、方法、设备及终端

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105888659A (zh) * 2016-06-06 2016-08-24 中国石油大学(北京) 确定岩性油气藏成藏概率的方法及装置
CN106022946A (zh) * 2016-06-06 2016-10-12 中国石油大学(北京) 确定岩性地层圈闭油气成藏概率的方法及装置
CN106886058A (zh) * 2017-04-17 2017-06-23 中国石油大学(北京) 致密油藏有效储层物性下限值的确定方法和装置
CN106991279A (zh) * 2017-03-29 2017-07-28 中国石油大学(北京) 一种致密砂岩气成藏概率的确定方法和装置

Family Cites Families (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB9813760D0 (en) 1998-06-25 1998-08-26 Geco Prakla Uk Ltd Seismic data signal processing method
US7467044B2 (en) * 2007-01-15 2008-12-16 Chevron U.S.A. Inc Method and system for assessing exploration prospect risk and uncertainty
CN102865081B (zh) 2012-04-28 2015-07-15 中国神华能源股份有限公司 一种保水开采方法
CN103901502B (zh) 2014-03-11 2016-06-08 中国海洋石油总公司 一种基于复合砂体构型图版的海上油田储层预测方法
CN105443119A (zh) 2014-08-19 2016-03-30 中国石油化工股份有限公司 一种获取油气藏储量规模阈值的方法
CN104573339B (zh) 2014-12-24 2018-03-09 中国石油大学(北京) 页岩气储层的地质参数确定方法和装置
CN105488583B (zh) 2015-11-23 2019-11-08 中国石油天然气股份有限公司 预测致密油待评价区域可采储量的方法及装置
CN105697003A (zh) 2016-03-17 2016-06-22 成都创源油气技术开发有限公司 页岩气储量评价方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105888659A (zh) * 2016-06-06 2016-08-24 中国石油大学(北京) 确定岩性油气藏成藏概率的方法及装置
CN106022946A (zh) * 2016-06-06 2016-10-12 中国石油大学(北京) 确定岩性地层圈闭油气成藏概率的方法及装置
CN106991279A (zh) * 2017-03-29 2017-07-28 中国石油大学(北京) 一种致密砂岩气成藏概率的确定方法和装置
CN106886058A (zh) * 2017-04-17 2017-06-23 中国石油大学(北京) 致密油藏有效储层物性下限值的确定方法和装置

Also Published As

Publication number Publication date
CN107657365A (zh) 2018-02-02
US11037084B2 (en) 2021-06-15
US20190080273A1 (en) 2019-03-14

Similar Documents

Publication Publication Date Title
CN107657365B (zh) 地质资源开采价值评价方法及装置
RU2669948C2 (ru) Оптимизация многоступенчатого проекта нефтяного месторождения в условиях неопределенности
CN106285660B (zh) 一种多层砂岩油藏低阻油层识别方法及装置
US8332155B2 (en) System and method for hydrocarbon gas pay zone characterization in a subterranean reservoir
Wang et al. Improved methods for determining effective sandstone reservoirs and evaluating hydrocarbon enrichment in petroliferous basins
Elkatatny A self-adaptive artificial neural network technique to predict total organic carbon (TOC) based on well logs
Bhatti et al. Permeability prediction using hydraulic flow units and electrofacies analysis
Liu et al. Permeability predictions for tight sandstone reservoir using explainable machine learning and particle swarm optimization
Wang et al. Determination of total organic carbon content in shale formations with regression analysis
CN104239743B (zh) 确定岩性油气藏成藏概率的方法和装置
J Al-Mudhafar Integrating lithofacies and well logging data into smooth generalized additive model for improved permeability estimation: Zubair formation, South Rumaila oil field
Tariq et al. Data-driven approaches to predict thermal maturity indices of organic matter using artificial neural networks
Obidegwu et al. Seismic assisted history matching using binary image matching
US11649721B2 (en) Hydrocarbon evaluation systems
Siddig et al. [Retracted] Utilization of Artificial Neural Network in Predicting the Total Organic Carbon in Devonian Shale Using the Conventional Well Logs and the Spectral Gamma Ray
Yasin et al. Fault and fracture network characterization using seismic data: a study based on neural network models assessment
Siddig et al. Application of various machine learning techniques in predicting total organic carbon from well logs
Ouenes et al. Estimation of propped volume permeability using strain from geomechanical modeling of interacting hydraulic and natural fractures-application to the eagle ford
CN108227033B (zh) 砂岩储层的孔隙度的预测方法、装置及计算机存储介质
Hudson et al. Unsupervised machine learning for detecting soil layer boundaries from cone penetration test data
Jayne et al. Permeability correlation structure of the Columbia River Plateau and implications for fluid system architecture in continental large igneous provinces
Zheng et al. Study on reservoir properties and critical depth in deep coal seams in Qinshui Basin, China
CN106285661B (zh) 一种基于判定指数的低阻油层识别方法及装置
Wang et al. Quantitative evaluation of unconsolidated sandstone heavy oil reservoirs based on machine learning
Lomask et al. A seismic to simulation unconventional workflow using automated fault-detection attributes

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