CN110807573B - 一种海相碳酸盐岩油气有利勘探区的确定方法及系统 - Google Patents

一种海相碳酸盐岩油气有利勘探区的确定方法及系统 Download PDF

Info

Publication number
CN110807573B
CN110807573B CN201910942667.6A CN201910942667A CN110807573B CN 110807573 B CN110807573 B CN 110807573B CN 201910942667 A CN201910942667 A CN 201910942667A CN 110807573 B CN110807573 B CN 110807573B
Authority
CN
China
Prior art keywords
dredging
fault
reservoir
efficiency
determining
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
CN201910942667.6A
Other languages
English (en)
Other versions
CN110807573A (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 University of Petroleum Beijing
Original Assignee
China University of Petroleum Beijing
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 China University of Petroleum Beijing filed Critical China University of Petroleum Beijing
Priority to CN201910942667.6A priority Critical patent/CN110807573B/zh
Publication of CN110807573A publication Critical patent/CN110807573A/zh
Application granted granted Critical
Publication of CN110807573B publication Critical patent/CN110807573B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION 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/00Administration; Management
    • G06Q10/06Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
    • G06Q10/063Operations research, analysis or management
    • G06Q10/0639Performance analysis of employees; Performance analysis of enterprise or organisation operations
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V9/00Prospecting or detecting by methods not provided for in groups G01V1/00 - G01V8/00
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION 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/00Administration; Management
    • G06Q10/04Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION 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/00Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/02Agriculture; Fishing; Forestry; Mining

Landscapes

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

Abstract

本发明提供了一种海相碳酸盐岩油气有利勘探区的确定方法及系统,所述方法包括:获取断层发育的海相碳酸盐岩断层基础地质资料;根据所述断层基础地质资料确定断层垂向疏导体系疏导效率和储层横向疏导体系疏导效率,确定海相碳酸盐岩目标工区所有探井的探井单井相对产能;根据探井单井相对产能、所述断层垂向疏导体系疏导效率和储层横向疏导体系疏导效率确定断层‑储层复合疏导体系疏导效率;根据所述断层‑储层复合疏导体系疏导效率确定海相碳酸盐岩目标工区油气有利勘探区,本发明可提高海相碳酸盐岩油气勘探油气有利勘探区确定的准确性。

Description

一种海相碳酸盐岩油气有利勘探区的确定方法及系统
技术领域
本发明涉及确定油气有利勘探区技术领域,尤其涉及一种海相碳酸盐岩油气有利勘探区的确定方法及系统。
背景技术
油气运移是油气成藏研究最重要的一环,油气运移通道构成了疏导体系,疏导体系的疏导效率决定了油气成藏的难易程度及规模。疏导体系一般包括断层垂向运移通道组形成的垂向疏导体系,相对高孔渗储层横向运移通道形成的横向疏导体系,以及两者联合形成的复合疏导体系。海相碳酸盐岩地层油气成藏规律十分复杂,储层广泛发育裂缝与溶孔,油气既可通过断层疏导体系垂向运移,也可通过相对高孔渗储层疏导体系横向运移,断层-储层复合疏导体系联合控制着油气运移的有利区。定量表征断层-储层复合疏导效率对于预测海相碳酸盐岩有利勘探区具有重要意义。
前人分别对砂岩断层、砂体输导体系的输导效率进行过一些研究。其中,对断层垂向输导体系的疏导效率提出了相应的计算模型,砂体横向输导体系的输导效率研究主要停留在定性层面,目前,复合疏导体系疏导效率的研究比较欠缺,尚未出现一个可靠的计算模式来定量表征复合疏导体系疏导效率。在海相碳酸盐岩地层,广泛发育裂缝与溶孔,油气既可通过断层垂向疏导体系运移,也可通过相对高孔渗储层横向疏导体系运移,裂缝和相对高孔渗储层组成的复合疏导体系控制着油气的富集成藏,油气藏分布规律十分复杂,油气藏有利勘探区预测很困难。基于断层-储层复合疏导体系疏导效率,可以指明油气有利勘探区,能有效地指导海相碳酸盐岩油气的勘探。
发明内容
本发明的一个目的在于提供一种海相碳酸盐岩油气有利勘探区的确定方法,可提高海相碳酸盐岩油气勘探油气有利勘探区确定的准确性。本发明的另一个目的在于提供一种海相碳酸盐岩油气有利勘探区的确定系统。
为了达到以上目的,本发明一方面公开了一种海相碳酸盐岩油气有利勘探区的确定方法,包括:
获取断层发育的海相碳酸盐岩断层基础地质资料;
根据所述断层基础地质资料确定断层垂向疏导体系疏导效率和储层横向疏导体系疏导效率,确定海相碳酸盐岩目标工区所有探井的探井单井相对产能;
根据探井单井相对产能、所述断层垂向疏导体系疏导效率和储层横向疏导体系疏导效率确定断层-储层复合疏导体系疏导效率;
根据所述断层-储层复合疏导体系疏导效率确定海相碳酸盐岩目标工区油气有利勘探区。
优选的,所述根据所述断层基础地质资料确定断层垂向疏导体系疏导效率具体包括:
根据断层疏导系数确定模型和所述断层基础地质资料得到目标工区的断层垂向疏导体系疏导系数;
将断层垂向疏导体系疏导系数归一化处理,得到断层垂向疏导体系疏导效率。
优选的,根据所述断层基础地质资料确定储层横向疏导体系疏导效率具体包括:
根据储层疏导系数确定模型和所述断层基础地质资料得到目标工区的储层横向疏导体系疏导系数;
将储层横向疏导体系疏导系数归一化处理,得到储层横向疏导体系疏导效率。
优选的,所述确定海相碳酸盐岩目标工区所有探井的探井单井相对产能具体包括:
获取目标工区所有探井的单井产能;
将目标工区所有探井的单井产能归一化处理,得到探井单井相对产能。
优选的,所述断层疏导系数确定模型的计算公式为:
Figure GDA0003494814630000021
Figure GDA0003494814630000022
其中,Tf为断层垂向疏导体系输导系数,无量纲;Ci为水力半径,km;Li为第i条断层延伸长度,km;Hi为断层的断距,km;Ti为断层宽度,Ti=0.62×Hi×0.875i,km;θi为i条断层的倾角,度;S为研究区面积,km2;n为断层条数。
优选的,所述将断层垂向疏导体系疏导系数归一化处理,得到断层垂向疏导体系疏导效率通过以下公式实现:
Ef=Tfi/Tfmax
其中,Ef为断裂垂向疏导体系输导效率,0-1之间;Tfi为第i条断裂垂向输导体系疏导系数,无量纲;Tfmax为最大的断裂垂向疏导体系输导系数,无量纲。
优选的,所述储层疏导系数确定模型的计算公式为:
Tr=Hr×Kr
其中,Tr为碳酸盐岩储层横向疏导体系输导系数,无量纲;Hr为有效碳酸盐岩储层厚度,m;Kr为储层有效孔隙度,%。
优选的,所述将储层横向疏导体系疏导系数归一化处理,得到储层横向疏导体系疏导效率通过以下公式得到:
Er=Tri/Trmax
其中,Er为储层横向疏导体系输导效率,0-1之间;Tri为第i个储层横向疏导体系输导系数,无量纲;Trmax为最大储层疏导体系输导系数,无量纲。
优选的,所述根据探井单井相对产能、所述断层垂向疏导体系疏导效率和储层横向疏导体系疏导效率确定断层-储层复合疏导体系疏导效率通过以下公式实现:
Figure GDA0003494814630000031
其中,Ice为断层-储层复合疏导体系疏导效率,0~1之间;Ef为断层垂向疏导体系疏导效率,0~1之间;Er为白云岩储层横向疏导体系疏导效率,0~1之间。
本发明还公开了一种海相碳酸盐岩油气有利勘探区的确定系统,包括:
资料获取单元,用于获取断层发育的海相碳酸盐岩断层基础地质资料;
参数确定单元,用于根据所述断层基础地质资料确定断层垂向疏导体系疏导效率和储层横向疏导体系疏导效率,确定海相碳酸盐岩目标工区所有探井的探井单井相对产能;
复合疏导体系疏导效率确定单元,根据探井单井相对产能、所述断层垂向疏导体系疏导效率和储层横向疏导体系疏导效率确定断层-储层复合疏导体系疏导效率;
有利区确定单元,用于根据所述断层-储层复合疏导体系疏导效率确定海相碳酸盐岩目标工区油气有利勘探区。
本发明利用统计归纳法,公式法,数值分析法,分别计算断层垂向疏导体系疏导系数和疏导效率、储层横向疏导体系疏导系数和疏导效率。根据两种疏导体系疏导效率和油气藏指标(油气藏个数、单井产量、探井成功率等)之间的定量关系,回归分析得到碳酸盐岩断层-储层复合疏导体系疏导效率。复合疏导体系疏导效率和油气藏产能具有正相关关系,疏导效率越高,油气聚集量越大,形成的油气规模越大,是越有利于油气勘探的区块。该发明以地质基本原理出发,运用数学理论模型,综合运用统计法、归纳法、公式法,定量表征碳酸盐岩断层-储层复合疏导体系疏导效率,对海相碳酸盐岩油气勘探具有很大的指导意义。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1示出研究区位置及断层发育特征。
图2示出本发明一种海相碳酸盐岩油气有利勘探区的确定方法具体实施例的流程图之一。
图3示出本发明一种海相碳酸盐岩油气有利勘探区的确定方法具体实施例的流程图之二。
图4示出本发明一种海相碳酸盐岩油气有利勘探区的确定方法具体实施例的流程图之三。
图5示出本发明一种海相碳酸盐岩油气有利勘探区的确定方法具体实施例的流程图之四。
图6示出川中地区断层垂向疏导体系疏导系数分布。
图7示出川中地区断层垂向疏导体系疏导效率分布。
图8示出川中地区寒武系龙王庙组白云岩储层横向疏导体系疏导系数分布。
图9示出川中地区寒武系龙王庙组白云岩储层横向疏导体系疏导效率分布。
图10示出川中地区龙王庙组探井相对产能和断层垂向疏导体系疏导效率的关系。
图11示出川中地区龙王庙组探井相对产能和白云岩储层横向疏导体系疏导效率的关系。
图12示出川中地区龙王庙组白云岩断层-储层复合疏导体系疏导效率分布。
图13示出可靠性检验。
图14示出探井产能和断层-储层复合疏导体系疏导效率相关性。
图15示出本发明一种海相碳酸盐岩油气有利勘探区的确定系统具体实施例的结构图。
图16示出适于用来实现本发明实施例的计算机设备的结构示意图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
如图1所示,碳酸盐岩油气运移疏导体系很复杂,既有断层形成的垂向疏导体系,也有高孔渗储层形成的横向疏导体系,两种疏导体系交织形成的复合疏导体系疏导效率很难被定量表征,本发明根据盆地内断层发育特征、储层分布特征、探井产能数据,提出了一种定量表征碳酸盐岩断层-储层复合疏导体系疏导效率进行海相碳酸盐岩油气有利勘探区的确定方法。
根据本发明的一个方面,本实施例公开了一种海相碳酸盐岩油气有利勘探区的确定方法。如图2所示,本实施例中,所述方法包括:
S100:获取断层发育的海相碳酸盐岩断层基础地质资料。
S200:根据所述断层基础地质资料确定断层垂向疏导体系疏导效率和储层横向疏导体系疏导效率,确定海相碳酸盐岩目标工区所有探井的探井单井相对产能。海相碳酸盐岩油气研究区存在断层形成的垂向疏导体系和高孔渗储层形成的横向疏导体系,两种疏导体系交汇,形成网络状的纵横相交的断层-储层复合疏导体系。
S300:根据探井单井相对产能、所述断层垂向疏导体系疏导效率和储层横向疏导体系疏导效率确定断层-储层复合疏导体系疏导效率。
S400:根据所述断层-储层复合疏导体系疏导效率确定海相碳酸盐岩目标工区油气有利勘探区。其中,有利勘探区可表征储层资源,根据断层-储层复合疏导体系疏导效率确定海相碳酸盐岩目标工区油气有利勘探区可采用现有方法进行确定,例如可根据断层-储层复合疏导体系疏导效率的数值是否到达预定阈值来确定有利勘探区的范围,当然也可采用现有的其他方式进行有利勘探区的确定。
本发明利用统计归纳法,公式法,数值分析法,分别计算断层垂向疏导体系疏导系数和疏导效率、储层横向疏导体系疏导系数和疏导效率。根据两种疏导体系疏导效率和油气藏指标(油气藏个数、单井产量、探井成功率等)之间的定量关系,回归分析得到碳酸盐岩断层-储层复合疏导体系疏导效率。复合疏导体系疏导效率和油气藏产能具有正相关关系,疏导效率越高,油气聚集量越大,形成的油气规模越大,是越有利于油气勘探的区块。该发明以地质基本原理出发,运用数学理论模型,综合运用统计法、归纳法、公式法,定量表征碳酸盐岩断层-储层复合疏导体系疏导效率,对海相碳酸盐岩油气勘探具有很大的指导意义。
在优选的实施方式中,如图3所示,所述S200中根据所述断层基础地质资料确定断层垂向疏导体系疏导效率具体包括:
S211:根据断层疏导系数确定模型和所述断层基础地质资料得到目标工区的断层垂向疏导体系疏导系数。
S212:将断层垂向疏导体系疏导系数归一化处理,得到断层垂向疏导体系疏导效率。
具体的,断层对油气的垂向疏导系数可以用断裂带宽度、断裂密度、延伸长度和断层倾角等断层参数定量表征。将断层疏导系数归一化处理,即可以得到断层疏导效率。
在优选的实施方式中,根据孙永河2007年在《中国石油大学学报(自然科学版)》提出来的计算公式,结合断层延伸长度、断距、倾角,可建立所述断层疏导系数确定模型,其计算公式为:
Figure GDA0003494814630000061
Figure GDA0003494814630000062
其中,Tf为断层垂向疏导体系输导系数,无量纲;Ci为水力半径,km;Li为第i条断层延伸长度,km;Hi为断层的断距,km;Ti为断层宽度,Ti=0.62×Hi×0.875i,km;θi为i条断层的倾角,度;S为研究区面积,km2;n为断层条数。
在优选的实施方式中,所述将断层垂向疏导体系疏导系数归一化处理,得到断层垂向疏导体系疏导效率可通过以下公式实现:
Ef=Tfi/Tfmax
其中,Ef为断裂垂向疏导体系输导效率,0-1之间;Tfi为第i条断裂垂向输导体系疏导系数,无量纲;Tfmax为最大的断裂垂向疏导体系输导系数,无量纲。
在优选的实施方式中,如图4所示,根据S200中所述断层基础地质资料确定储层横向疏导体系疏导效率具体包括:
S221:根据储层疏导系数确定模型和所述断层基础地质资料得到目标工区的储层横向疏导体系疏导系数。
S222:将储层横向疏导体系疏导系数归一化处理,得到储层横向疏导体系疏导效率。
具体的,储层对油气的横向疏导系数可以有效储层的厚度和有效孔隙度等参数定量表征。将储层疏导系数归一化处理,即可以得到储层疏导效率。
在优选的实施方式中,所述储层疏导系数确定模型的计算公式为:
Tr=Hr×Kr
其中,Tr为碳酸盐岩储层横向疏导体系输导系数,无量纲;Hr为有效碳酸盐岩储层厚度,m;Kr为储层有效孔隙度,%。
在优选的实施方式中,所述将储层横向疏导体系疏导系数归一化处理,得到储层横向疏导体系疏导效率通过以下公式得到:
Er=Tri/Trmax
其中,Er为储层横向疏导体系输导效率,0-1之间;Tri为第i个储层横向疏导体系输导系数,无量纲;Trmax为最大储层疏导体系输导系数,无量纲。
在优选的实施方式中,如图5所示,所述S200中确定海相碳酸盐岩目标工区所有探井的探井单井相对产能具体包括:
S231:获取目标工区所有探井的单井产能。
S232:将目标工区所有探井的单井产能归一化处理,得到探井单井相对产能。
具体的,将收集到的探井单井产能数据,通过归一化处理,得到单井相对产能。根据单井相对产能和断裂垂向疏导体系输导效率,储层横向疏导体系输导效率之间的关系,回归分析,建立复合疏导体系疏导效率。
在优选的实施方式中,所述将目标工区所有探井的单井产能归一化处理,得到探井单井相对产能通过以下公式实现:
Pg=Pgi/Pgmax
其中,Pg为单井相对产能,无量纲,0-1之间;Pgi为i口探井日产气量,m3/d;Pgmax为最大探井日产气量,m3/d。
在优选的实施方式中,所述根据探井单井相对产能、所述断层垂向疏导体系疏导效率和储层横向疏导体系疏导效率确定断层-储层复合疏导体系疏导效率通过以下公式实现:
Figure GDA0003494814630000081
其中,Ice为断层-储层复合疏导体系疏导效率,0~1之间;Ef为断层垂向疏导体系疏导效率,0~1之间;Er为白云岩储层横向疏导体系疏导效率,0~1之间。
其中,在确定断层-储层复合疏导体系疏导效率时,可先建立断层-储层复合疏导体系疏导效率定量表征模型,其表达式可为:
Ice=af(Ef)+bf(Er)+c
其中,Ice为断层-储层复合疏导体系疏导效率,无量纲,0-1之间;f(Ef)、f(Er)是Pg关于Ef,Er的函数;a,b,c为常数。
将探井单井产能和断层垂向疏导体系输导效率,储层横向疏导体系输导效率进行回归分析,可求得断层-储层复合疏导体系疏导效率。
下面通过一个具体例子来对本发明作进一步的说明,四川盆地是中国最主要的含油气盆地之一,油气资源十分丰富,盆地面积约为19×104km3。本研究区川中地区主要位于四川省资阳市安岳县、遂宁市境内,资阳市乐至县以东,重庆市潼南县以西,面积约2330km2,如图1所示。2013年,在四川盆地川中地区发现了著名的安岳震旦系-寒武系特大型气田,其中,寒武系龙王庙组天然气资源量高达1.2×1012m3,天然气探明储量4404×108m3。本研究目的层为寒武系龙王庙组,其岩性为颗粒白云岩夹少量砂泥岩。龙王庙组与下部地层筇竹寺组、上部地层高台组形成良好的下生中储上盖生储盖组合关系。在龙王庙组之下发育烃源岩,即泥质粉砂岩、砂质炭质页岩、页岩发育的筇竹寺组;在龙王庙组之上发育区域盖层,即页岩、泥质白云质灰岩及石膏发育的高台组。川中地区寒武系发育14条主要断层,全部为正断层,断层延伸长度在5-10km,断距为20-40m。龙王庙组白云岩有效储层厚度分布在10-60m,储层有效孔隙度主要在2.5%-7%。川中地区断层和有效储层联合控制着油气成藏,其疏导效率影响油气的富集,如何定量评价断层-储层复合疏导体系疏导效率对于指导川中地区寒武系龙王庙下一步油气勘探具有重要的意义。
本发明分别计算已发现油气藏(以探井产能作为指标)和断层垂向疏导体系疏导效率、储层横向疏导体系疏导效率之间的定量关系,最终确定川中地区寒武系龙王庙组断层-储层复合疏导体系疏导效率。
确定影响川中地区寒武系龙王庙油气富集成藏的主控因素。通过地质剖析,本研究确定了断层,白云岩有效储层形成的疏导系统控制着该区的油气富集。断层垂向体系疏导能力越强,白云岩储层横向体系疏导能力越强,油气越易富集,两者组成的网络状复合疏导体系疏导能力控制着油气最富集的地区。
计算断层垂向疏导体系疏导系数(Tf)和疏导效率(Ef)。过程如下:
计算断层垂向疏导体系疏导系数(Tf)。影响断层输导性能的因素主要有断裂带宽度、断裂密度、延伸长度、断层倾角以及断层的密度等。断裂带越宽,延伸越长,断距越大,断裂密度越大,断层疏导性能越好。根据已有断层疏导系数计算模型,计算川中寒武系断层疏导体系疏导系数。图6示出了川中地区断层垂向疏导体系疏导系数的分布。
计算断层垂向疏导体系疏导效率(Ef)。对断层垂向疏导体系输导系数进行归一化处理,将得到的所有断层疏导系数除以最大的疏导系数即可得到疏导效率的值,其值分布在0-1之间。图7示出了川中地区断层垂向疏导体系疏导效率分布。
计算白云岩储层横向疏导体系疏导系数(Tr)和疏导效率(Er)。过程如下:
计算白云岩储层横向疏导体系疏导系数(Tr)。相对高孔渗白云岩储层即是储存油气的空间,也是油气运移的通道。根据有效白云岩储层的厚度和有效孔隙度的乘积来定量表征白云岩储层对油气的横向运移能力,即白云岩储层疏导系数Tr。图8示出了川中地区寒武系龙王庙组白云岩储层横向疏导体系疏导系数分布。
计算白云岩储层横向疏导体系疏导效率(Er)。对白云岩储层横向疏导体系输导系数进行归一化处理,将得到的所有白云岩储层疏导系数除以最大的疏导系数即可得到疏导效率的值,其值分布在0-1之间。图9示出了川中地区寒武系龙王庙组白云岩储层横向疏导体系疏导效率分布;
计算断层-储层复合疏导体系疏导效率(Ice)。将收集到的探井单井产能数据,通过归一化处理,即将每口已钻探井的日产气量除以最大的探井的日产气量得到单井相对产能。然后建立单井相对产能和断裂垂向疏导体系输导效率,储层横向疏导体系输导效率之间的关系,最后二元回归分析,建立复合疏导体系疏导效率计算公式。过程如下:
对川中寒武系龙王庙17口已钻探井产能进行归一化处理,得到单井相对产能(Pg)。
根据单井相对产能(Pg)和断层垂向疏导体系疏导效率(Ef)之间的分布关系,如图10所示,确定两者的函数关系:
Figure GDA0003494814630000102
其中,Pg为单井相对产能,0~1之间;Ef为断层垂向疏导体系疏导效率,0~1之间。
根据单井相对产能(Pg)和白云岩储层横向疏导体系疏导效率(Er)之间的分布关系,如图11所示,确定两者的函数关系:
Figure GDA0003494814630000103
其中,Pg为单井相对产能,0~1之间;Er为白云岩储层横向疏导体系疏导效率,0~1之间。
将探井单井产能(Pg)和断裂垂向疏导体系输导效率(Ef),储层横向疏导体系输导效率进行回归分析(Er),建立断层-储层复合疏导体系疏导效率的定量表征模型(Ice),计算公式是:
Figure GDA0003494814630000101
其中,Ice为断层-储层复合疏导体系疏导效率,0~1之间;Ef为断层垂向疏导体系疏导效率,0~1之间;Er为白云岩储层横向疏导体系疏导效率,0~1之间。
根据断层-储层复合疏导体系疏导效率计算公式,计算研究区川中地区龙王庙组复合疏导体系疏导效率。图12示出了川中地区龙王庙组白云岩断层-储层复合疏导体系疏导效率分布。
根据断层-储层复合疏导体系疏导效率计算公式,预测的探井产能并与已钻17口探井产能做对比,结果如图13所示,预测的探井产能与实际探井产能数值匹配度很高,说明该模型及计算公式具有很高的精确度和可靠性。统计复合疏导体系疏导效率和探井产能之间的关系如图14所示,两者具有正相关关系,复合疏导体系疏导效率越高,探井产能越大,反映了油气越易富集成藏。基于复合疏导体系疏导效率的大小,可以指出川中地区寒武系龙王庙组下一步油气有利勘探区。
本发明实施例通过定量表征断层-储层复合疏导体系疏导效率,解决了四川盆地川中地区寒武系龙王庙组确定下一步油气有利勘探区带选择的难题,探井结果验证该方法具有很高的精确度和可靠性。本方法为海相碳酸盐岩油气勘探提供了一种可行的技术方法,降低了油气勘探的风险,指明了油气的勘探方向,具有广泛的适用性。
基于相同原理,本实施例还公开了一种海相碳酸盐岩油气有利勘探区的确定系统。如图15所示,本实施例中,所述系统包括资料获取单元11、参数确定单元12、复合疏导体系疏导效率确定单元13和有利区确定单元14。
其中,资料获取单元11用于获取断层发育的海相碳酸盐岩断层基础地质资料。
参数确定单元12用于根据所述断层基础地质资料确定断层垂向疏导体系疏导效率和储层横向疏导体系疏导效率,确定海相碳酸盐岩目标工区所有探井的探井单井相对产能。
复合疏导体系疏导效率确定单元13根据探井单井相对产能、所述断层垂向疏导体系疏导效率和储层横向疏导体系疏导效率确定断层-储层复合疏导体系疏导效率。
有利区确定单元14用于根据所述断层-储层复合疏导体系疏导效率确定海相碳酸盐岩目标工区油气有利勘探区。
由于该系统解决问题的原理与以上方法类似,因此本系统的实施可以参见方法的实施,在此不再赘述。
上述实施例阐明的系统、装置、模块或单元,具体可以由计算机芯片或实体实现,或者由具有某种功能的产品来实现。一种典型的实现设备为计算机设备,具体的,计算机设备例如可以为个人计算机、膝上型计算机、蜂窝电话、相机电话、智能电话、个人数字助理、媒体播放器、导航设备、电子邮件设备、游戏控制台、平板计算机、可穿戴设备或者这些设备中的任何设备的组合。
在一个典型的实例中计算机设备具体包括存储器、处理器以及存储在存储器上并可在处理器上运行的计算机程序,所述处理器执行所述程序时实现如上所述的方法,或者,所述处理器执行所述程序时实现如上所述的方法。
下面参考图16,其示出了适于用来实现本申请实施例的计算机设备600的结构示意图。
如图16所示,计算机设备600包括中央处理单元(CPU)601,其可以根据存储在只读存储器(ROM)602中的程序或者从存储部分608加载到随机访问存储器(RAM))603中的程序而执行各种适当的工作和处理。在RAM603中,还存储有系统600操作所需的各种程序和数据。CPU601、ROM602、以及RAM603通过总线604彼此相连。输入/输出(I/O)接口605也连接至总线604。
以下部件连接至I/O接口605:包括键盘、鼠标等的输入部分606;包括诸如阴极射线管(CRT)、液晶反馈器(LCD)等以及扬声器等的输出部分607;包括硬盘等的存储部分608;以及包括诸如LAN卡,调制解调器等的网络接口卡的通信部分609。通信部分609经由诸如因特网的网络执行通信处理。驱动器610也根据需要连接至I/O接口605。可拆卸介质611,诸如磁盘、光盘、磁光盘、半导体存储器等等,根据需要安装在驱动器610上,以便于从其上读出的计算机程序根据需要被安装如存储部分608。
特别地,根据本发明的实施例,上文参考流程图描述的过程可以被实现为计算机软件程序。例如,本发明的实施例包括一种计算机程序产品,其包括有形地包含在机器可读介质上的计算机程序,所述计算机程序包括用于执行流程图所示的方法的程序代码。在这样的实施例中,该计算机程序可以通过通信部分609从网络上被下载和安装,和/或从可拆卸介质611被安装。
计算机可读介质包括永久性和非永久性、可移动和非可移动媒体可以由任何方法或技术来实现信息存储。信息可以是计算机可读指令、数据结构、程序的模块或其他数据。计算机的存储介质的例子包括,但不限于相变内存(PRAM)、静态随机存取存储器(SRAM)、动态随机存取存储器(DRAM)、其他类型的随机存取存储器(RAM)、只读存储器(ROM)、电可擦除可编程只读存储器(EEPROM)、快闪记忆体或其他内存技术、只读光盘只读存储器(CD-ROM)、数字多功能光盘(DVD)或其他光学存储、磁盒式磁带,磁带磁磁盘存储或其他磁性存储设备或任何其他非传输介质,可用于存储可以被计算设备访问的信息。按照本文中的界定,计算机可读介质不包括暂存电脑可读媒体(transitory media),如调制的数据信号和载波。
为了描述的方便,描述以上装置时以功能分为各种单元分别描述。当然,在实施本申请时可以把各单元的功能在同一个或多个软件和/或硬件中实现。
本发明是参照根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器,使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中,使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品,该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上,使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理,从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
还需要说明的是,术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含,从而使得包括一系列要素的过程、方法、商品或者设备不仅包括那些要素,而且还包括没有明确列出的其他要素,或者是还包括为这种过程、方法、商品或者设备所固有的要素。在没有更多限制的情况下,由语句“包括一个……”限定的要素,并不排除在包括所述要素的过程、方法、商品或者设备中还存在另外的相同要素。
本领域技术人员应明白,本申请的实施例可提供为方法、系统或计算机程序产品。因此,本申请可采用完全硬件实施例、完全软件实施例或结合软件和硬件方面的实施例的形式。而且,本申请可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本申请可以在由计算机执行的计算机可执行指令的一般上下文中描述,例如程序模块。一般地,程序模块包括执行特定任务或实现特定抽象数据类型的例程、程序、对象、组件、数据结构等等。也可以在分布式计算环境中实践本申请,在这些分布式计算环境中,由通过通信网络而被连接的远程处理设备来执行任务。在分布式计算环境中,程序模块可以位于包括存储设备在内的本地和远程计算机存储介质中。
本说明书中的各个实施例均采用递进的方式描述,各个实施例之间相同相似的部分互相参见即可,每个实施例重点说明的都是与其他实施例的不同之处。尤其,对于系统实施例而言,由于其基本相似于方法实施例,所以描述的比较简单,相关之处参见方法实施例的部分说明即可。
以上所述仅为本申请的实施例而已,并不用于限制本申请。对于本领域技术人员来说,本申请可以有各种更改和变化。凡在本申请的精神和原理之内所作的任何修改、等同替换、改进等,均应包含在本申请的权利要求范围之内。

Claims (9)

1.一种海相碳酸盐岩油气有利勘探区的确定方法,其特征在于,包括:
获取断层发育的海相碳酸盐岩断层基础地质资料;
根据所述断层基础地质资料确定断层垂向疏导体系疏导效率和储层横向疏导体系疏导效率,确定海相碳酸盐岩目标工区所有探井的探井单井相对产能;
根据探井单井相对产能、所述断层垂向疏导体系疏导效率和储层横向疏导体系疏导效率确定断层-储层复合疏导体系疏导效率;
根据所述断层-储层复合疏导体系疏导效率确定海相碳酸盐岩目标工区油气有利勘探区;
所述根据探井单井相对产能、所述断层垂向疏导体系疏导效率和储层横向疏导体系疏导效率确定断层-储层复合疏导体系疏导效率通过以下公式实现:
Figure FDA0003494814620000011
其中,Ice为断层-储层复合疏导体系疏导效率,0~1之间;Ef为断层垂向疏导体系疏导效率,0~1之间;Er为白云岩储层横向疏导体系疏导效率,0~1之间;
其中,在确定断层-储层复合疏导体系疏导效率时,先建立断层-储层复合疏导体系疏导效率定量表征模型,其表达式可为:
Ice=af(Ef)+bf(Er)+c
其中,Ice为断层-储层复合疏导体系疏导效率,无量纲,0-1之间;f(Ef)、f(Er)是Pg关于Ef,Er的函数;a,b,c为常数;
Figure FDA0003494814620000012
其中,Pg为单井相对产能,0~1之间;
Figure FDA0003494814620000013
其中,Pg为单井相对产能,0~1之间。
2.根据权利要求1所述的海相碳酸盐岩油气有利勘探区的确定方法,其特征在于,所述根据所述断层基础地质资料确定断层垂向疏导体系疏导效率具体包括:
根据断层疏导系数确定模型和所述断层基础地质资料得到目标工区的断层垂向疏导体系疏导系数;
将断层垂向疏导体系疏导系数归一化处理,得到断层垂向疏导体系疏导效率。
3.根据权利要求1所述的海相碳酸盐岩油气有利勘探区的确定方法,其特征在于,根据所述断层基础地质资料确定储层横向疏导体系疏导效率具体包括:
根据储层疏导系数确定模型和所述断层基础地质资料得到目标工区的储层横向疏导体系疏导系数;
将储层横向疏导体系疏导系数归一化处理,得到储层横向疏导体系疏导效率。
4.根据权利要求1所述的海相碳酸盐岩油气有利勘探区的确定方法,其特征在于,所述确定海相碳酸盐岩目标工区所有探井的探井单井相对产能具体包括:
获取目标工区所有探井的单井产能;
将目标工区所有探井的单井产能归一化处理,得到探井单井相对产能。
5.根据权利要求2所述的海相碳酸盐岩油气有利勘探区的确定方法,其特征在于,所述断层疏导系数确定模型的计算公式为:
Figure FDA0003494814620000021
Figure FDA0003494814620000022
其中,Tf为断层垂向疏导体系输导系数,无量纲;Ci为水力半径,km;Li为第i条断层延伸长度,km;Hi为断层的断距,km;Ti为断层宽度,Ti=0.62×Hi×0.875i,km;θi为i条断层的倾角,度;S为研究区面积,km2;n为断层条数。
6.根据权利要求2所述的海相碳酸盐岩油气有利勘探区的确定方法,其特征在于,所述将断层垂向疏导体系疏导系数归一化处理,得到断层垂向疏导体系疏导效率通过以下公式实现:
Ef=Tfi/Tfmax
其中,Ef为断裂垂向疏导体系输导效率,0-1之间;Tfi为第i条断裂垂向输导体系疏导系数,无量纲;Tfmax为最大的断裂垂向疏导体系输导系数,无量纲。
7.根据权利要求3所述的海相碳酸盐岩油气有利勘探区的确定方法,其特征在于,所述储层疏导系数确定模型的计算公式为:
Tr=Hr×Kr
其中,Tr为碳酸盐岩储层横向疏导体系输导系数,无量纲;Hr为有效碳酸盐岩储层厚度,m;Kr为储层有效孔隙度,%。
8.根据权利要求3所述的海相碳酸盐岩油气有利勘探区的确定方法,其特征在于,所述将储层横向疏导体系疏导系数归一化处理,得到储层横向疏导体系疏导效率通过以下公式得到:
Er=Tri/Trmax
其中,Er为储层横向疏导体系输导效率,0-1之间;Tri为第i个储层横向疏导体系输导系数,无量纲;Trmax为最大储层疏导体系输导系数,无量纲。
9.一种海相碳酸盐岩油气有利勘探区的确定系统,其特征在于,包括:
资料获取单元,用于获取断层发育的海相碳酸盐岩断层基础地质资料;
参数确定单元,用于根据所述断层基础地质资料确定断层垂向疏导体系疏导效率和储层横向疏导体系疏导效率,确定海相碳酸盐岩目标工区所有探井的探井单井相对产能;
复合疏导体系疏导效率确定单元,根据探井单井相对产能、所述断层垂向疏导体系疏导效率和储层横向疏导体系疏导效率确定断层-储层复合疏导体系疏导效率;
有利区确定单元,用于根据所述断层-储层复合疏导体系疏导效率确定海相碳酸盐岩目标工区油气有利勘探区;
所述根据探井单井相对产能、所述断层垂向疏导体系疏导效率和储层横向疏导体系疏导效率确定断层-储层复合疏导体系疏导效率通过以下公式实现:
Figure FDA0003494814620000031
其中,Ice为断层-储层复合疏导体系疏导效率,0~1之间;Ef为断层垂向疏导体系疏导效率,0~1之间;Er为白云岩储层横向疏导体系疏导效率,0~1之间;
其中,在确定断层-储层复合疏导体系疏导效率时,先建立断层-储层复合疏导体系疏导效率定量表征模型,其表达式可为:
Ice=af(Ef)+bf(Er)+c
其中,Ice为断层-储层复合疏导体系疏导效率,无量纲,0-1之间;f(Ef)、f(Er)是Pg关于Ef,Er的函数;a,b,c为常数;
Figure FDA0003494814620000032
其中,Pg为单井相对产能,0~1之间;
Pg=0.0034×e5.6052×E r
其中,Pg为单井相对产能,0~1之间。
CN201910942667.6A 2019-09-30 2019-09-30 一种海相碳酸盐岩油气有利勘探区的确定方法及系统 Active CN110807573B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910942667.6A CN110807573B (zh) 2019-09-30 2019-09-30 一种海相碳酸盐岩油气有利勘探区的确定方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910942667.6A CN110807573B (zh) 2019-09-30 2019-09-30 一种海相碳酸盐岩油气有利勘探区的确定方法及系统

Publications (2)

Publication Number Publication Date
CN110807573A CN110807573A (zh) 2020-02-18
CN110807573B true CN110807573B (zh) 2022-03-29

Family

ID=69488193

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910942667.6A Active CN110807573B (zh) 2019-09-30 2019-09-30 一种海相碳酸盐岩油气有利勘探区的确定方法及系统

Country Status (1)

Country Link
CN (1) CN110807573B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112965137B (zh) * 2021-02-05 2023-01-10 中国石油天然气股份有限公司 海相碳酸盐岩地层中下切谷的识别方法及装置
CN114427458A (zh) * 2021-12-17 2022-05-03 中国石油化工股份有限公司 一种海相碳酸盐岩断控缝洞型油气藏勘探方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103728671A (zh) * 2014-01-15 2014-04-16 中国石油大学(北京) 确定碳酸盐岩潜山含油气性的方法和装置
CN107237627A (zh) * 2017-05-20 2017-10-10 陕西延长石油(集团)有限责任公司研究院 一种岩溶型碳酸盐岩储层布井方法
CN110259429A (zh) * 2019-06-06 2019-09-20 中国石油天然气股份有限公司 一种多级动态暂堵体积酸压方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107807407B (zh) * 2017-09-30 2019-10-11 中国石油天然气股份有限公司 一种油气区带有效性评价方法和装置

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103728671A (zh) * 2014-01-15 2014-04-16 中国石油大学(北京) 确定碳酸盐岩潜山含油气性的方法和装置
CN107237627A (zh) * 2017-05-20 2017-10-10 陕西延长石油(集团)有限责任公司研究院 一种岩溶型碳酸盐岩储层布井方法
CN110259429A (zh) * 2019-06-06 2019-09-20 中国石油天然气股份有限公司 一种多级动态暂堵体积酸压方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
库车坳陷北带断裂输导效率及其物理模拟实验研究;孙永河 等;《中国石油大学学报(自然科学版)》;20071231;第31卷(第6期);135-140,151 *
断裂-砂体复合输导体系输导效率的定量评价——以库车前陆盆地北带白垩系为例;尚晓庆 等;《石油天然气学报》;20120331;第34卷(第3期);38-43 *
长岭断陷深层致密气藏综合研究及潜力评价;张祎;《中国博士学位论文全文数据库基础科学辑》;20150215(第02期);A011-87 *

Also Published As

Publication number Publication date
CN110807573A (zh) 2020-02-18

Similar Documents

Publication Publication Date Title
CN107657365B (zh) 地质资源开采价值评价方法及装置
CN103645516B (zh) 基于岩石物理相控油气作用确定油气产能的方法
CN105528522A (zh) 一种基于空间网格计算连续性油气藏资源量的方法及装置
CN111695303B (zh) 一种煤层顶板砂岩含水层充水强度评价方法
CN110807573B (zh) 一种海相碳酸盐岩油气有利勘探区的确定方法及系统
CN110807541B (zh) 一种气藏开发井型的选择方法
Soleymanzadeh et al. Effect of overburden pressure on determination of reservoir rock types using RQI/FZI, FZI* and Winland methods in carbonate rocks
CN104834006B (zh) 结构成因沉积参数图的处理方法和装置
Gupta et al. Haynesville shale: predicting long-term production and residual analysis to identify well interference and fracture hits
CN107725044B (zh) 基于阵列感应、侧向测井的砂岩含气储层产水率预测的方法
Popa et al. Optimizing horizontal well placement through stratigraphic performance prediction using artificial intelligence
Jayne et al. Permeability correlation structure of the Columbia River Plateau and implications for fluid system architecture in continental large igneous provinces
Wang et al. Quantitative evaluation of unconsolidated sandstone heavy oil reservoirs based on machine learning
CN112949901A (zh) 页岩油藏水平井准自然能量开发方法及系统
CN114460639A (zh) 页岩油储层渗透率的预测方法及装置
Clark The prediction of hydraulic fracture azimuth through geological, core, and analytical studies
Wang et al. Prediction of stability coefficient of open-pit mine slope based on artificial intelligence deep learning algorithm
CN113914853B (zh) 深水沉积体系砂泥岩薄互层砂岩储层含水饱和度确定方法
CN106227959B (zh) 一种基于四图叠合法预测岩性油藏有利区的方法及装置
CN112505761B (zh) 一种储层含气性检测方法和装置
CN107942400B (zh) 一种预测钻井风险的方法及装置
CN113946931A (zh) 一种测井遇阻遇卡的主导因素确定方法及装置
CN110930020A (zh) 非常规油气资源的经济可采资源量确定方法
CN114622878B (zh) 特低-超低渗透储层注水开发适应性判别方法和装置
CN111781661A (zh) 少井地区沉积微相平面展布预测方法及装置

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant