CN108171376B - 总有机碳的预测方法、装置、电子设备及存储介质 - Google Patents

总有机碳的预测方法、装置、电子设备及存储介质 Download PDF

Info

Publication number
CN108171376B
CN108171376B CN201711450913.3A CN201711450913A CN108171376B CN 108171376 B CN108171376 B CN 108171376B CN 201711450913 A CN201711450913 A CN 201711450913A CN 108171376 B CN108171376 B CN 108171376B
Authority
CN
China
Prior art keywords
impedance
curve
organic carbon
total organic
longitudinal wave
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
CN201711450913.3A
Other languages
English (en)
Other versions
CN108171376A (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
Exploration and Development Research Institute of Sinopec Jianghan Oilfield Co
Original Assignee
China Petroleum and Chemical Corp
Exploration and Development Research Institute of Sinopec Jianghan Oilfield Co
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 Petroleum and Chemical Corp, Exploration and Development Research Institute of Sinopec Jianghan Oilfield Co filed Critical China Petroleum and Chemical Corp
Priority to CN201711450913.3A priority Critical patent/CN108171376B/zh
Publication of CN108171376A publication Critical patent/CN108171376A/zh
Application granted granted Critical
Publication of CN108171376B publication Critical patent/CN108171376B/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/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/06—Energy or water supply
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02P—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
    • Y02P90/00—Enabling technologies with a potential contribution to greenhouse gas [GHG] emissions mitigation
    • Y02P90/80—Management or planning
    • Y02P90/84—Greenhouse gas [GHG] management systems

Landscapes

  • Business, Economics & Management (AREA)
  • Engineering & Computer Science (AREA)
  • Economics (AREA)
  • Human Resources & Organizations (AREA)
  • Strategic Management (AREA)
  • Theoretical Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Marketing (AREA)
  • General Physics & Mathematics (AREA)
  • General Business, Economics & Management (AREA)
  • Tourism & Hospitality (AREA)
  • Public Health (AREA)
  • General Health & Medical Sciences (AREA)
  • Primary Health Care (AREA)
  • Water Supply & Treatment (AREA)
  • Development Economics (AREA)
  • Game Theory and Decision Science (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Operations Research (AREA)
  • Quality & Reliability (AREA)
  • Investigating Or Analyzing Materials By The Use Of Electric Means (AREA)

Abstract

本发明提供了一种总有机碳的预测方法、装置、电子设备及存储介质，涉及油气勘探技术领域。该方法包括：基于待测页岩气储层的纵波阻抗曲线以及梯度阻抗曲线拟合出在多个角度的总有机碳指示因子曲线；获取多个角度的总有机碳指示因子曲线中与待测页岩气储层的实测总有机碳曲线的相似度最高的总有机碳指示因子曲线对应的第一角度；基于扩展弹性阻抗计算公式以及待测页岩气储层的叠前地震道集数据进行反演，得到纵波阻抗体以及梯度阻抗体；将第一角度、纵波阻抗体以及梯度阻抗体代入预设总有机碳指示因子体表达式进行计算从而获得待测页岩气储层的总有机碳的预测含量。该总有机碳的预测方法可以使总有机碳的预测含量的精确度提高。

Description

总有机碳的预测方法、装置、电子设备及存储介质
技术领域
本发明涉及油气勘探技术领域，具体而言，涉及一种总有机碳的预测方法、装置、电子设备及存储介质。
背景技术
页岩气储层的TOC即为页岩气储层总有机碳，总有机碳含量是评价页岩有机质丰度和生烃能力的主要指标。TOC的含量影响了页岩气储层的有机孔隙度和吸附气的含量，一般来说，页岩中TOC越高，有机孔隙度和吸附气的含量就越大，因而，精确的页岩气储层的TOC含量的预测对页岩气藏开发尤为重要。
目前用地震资料直接定量预测页岩气储层的TOC方法是以地震属性与TOC建立数学拟合关系进行TOC预测，而由于地震属性存在固有的多解性，导致该方法预测的TOC结果不精确。
发明内容
有鉴于此，本发明实施例提供了一种总有机碳的预测方法、装置、电子设备及存储介质，由于考虑到地震道集资料中的梯度变化因素，使得页岩气储层的TOC定量预测更加可信可靠。
为了实现上述目的，本发明采用的技术方案如下：
第一方面，本发明实施例提供了一种总有机碳的预测方法，所述方法包括：基于待测页岩气储层的纵波阻抗曲线以及梯度阻抗曲线拟合出在多个角度的总有机碳指示因子曲线，其中，所述多个角度包括0°～360°范围内的各个角度；获取所述多个角度的总有机碳指示因子曲线中与所述待测页岩气储层的实测总有机碳曲线的相似度最高的总有机碳指示因子曲线对应的第一角度；基于扩展弹性阻抗计算公式以及所述待测页岩气储层的叠前地震道集数据进行反演，得到纵波阻抗体以及梯度阻抗体；将所述第一角度、所述纵波阻抗体以及所述梯度阻抗体代入预设总有机碳指示因子体表达式进行计算从而获得所述待测页岩气储层的总有机碳的预测含量。
第二方面，本发明实施例提供了一种总有机碳的预测装置，所述装置包括第一拟合模块、角度获取模块、反演执行模块以及含量计算模块，其中，所述第一拟合模块用于基于待测页岩气储层的纵波阻抗曲线以及梯度阻抗曲线拟合出在多个角度的总有机碳指示因子曲线，其中，所述多个角度包括0°～360°范围内的各个角度；所述角度获取模块用于获取所述多个角度的总有机碳指示因子曲线中与所述待测页岩气储层的实测总有机碳曲线的相似度最高的总有机碳指示因子曲线对应的第一角度；所述反演执行模块用于基于扩展弹性阻抗计算公式以及所述待测页岩气储层的叠前地震道集数据进行反演，得到纵波阻抗体以及梯度阻抗体；所述含量计算模块用于将所述第一角度、所述纵波阻抗体以及所述梯度阻抗体代入预设总有机碳指示因子体表达式进行计算从而获得所述待测页岩气储层的总有机碳的预测含量。
第三方面，本发明实施例提供了一种电子设备，所述电子设备包括存储器和处理器，所述存储器存储有计算机指令，当所述计算机指令由所述处理器读取并执行时，使所述处理器执行上述第一方面提供的总有机碳的预测方法。
第四方面，本发明实施例提供了一种存储介质，所述存储介质中存储有计算机指令，其中，所述计算机指令在被读取并运行时执行上述第一方面提供的总有机碳的预测方法。
本发明实施例提供的总有机碳的预测方法、装置、电子设备及存储介质，该方法通过基于待测页岩气储层的纵波阻抗曲线以及梯度阻抗曲线拟合出在多个角度的总有机碳指示因子曲线，其中，多个角度包括0°～360°范围内的各个角度，然后获取多个角度的总有机碳指示因子曲线中与待测页岩气储层的实测总有机碳曲线的相似度最高的总有机碳指示因子曲线对应的第一角度，再基于扩展弹性阻抗计算公式以及待测页岩气储层的叠前地震道集数据进行反演，得到纵波阻抗体以及梯度阻抗体，最后将第一角度、纵波阻抗体以及梯度阻抗体代入预设总有机碳指示因子体表达式进行计算从而获得待测页岩气储层的总有机碳的预测含量。该总有机碳的预测方法考虑了梯度阻抗因素，从而使预测的待测页岩气储层的总有机碳含量的精确度较高。
为使本发明的上述目的、特征和优点能更明显易懂，下文特举较佳实施例，并配合所附附图，作详细说明如下。
附图说明
为使本发明实施例的目的、技术方案和优点更加清楚，下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
图1示出了本发明实施例提供的电子设备的方框示意图；
图2示出了本发明实施例提供的总有机碳的预测方法的流程图；
图3示出了本发明实施例提供的总有机碳的预测方法中步骤S120的流程图；
图4示出了本发明实施例提供的各个角度总有机碳指示因子曲线与实测总有机碳曲线相似性关系的示意图；
图5示出了本发明实施例提供的有机碳实测曲线与最佳的总有机碳指示因子曲线示意图；
图6示出了本发明实施例提供的总有机碳的预测方法中步骤S130的流程图；
图7示出了本发明实施例提供的反演得到的纵波阻抗剖面的示意图；
图8示出了本发明实施例提供的反演得到的梯度阻抗剖面的示意图；
图9示出了本发明实施例提供的总有机碳指示因子体的连井剖面图；
图10示出了本发明实施例提供的总有机碳的预测装置的模块图；
图11示出了本发明实施例提供的总有机碳的预测装置中角度获取模块的模块图；
图12示出了本发明实施例提供的总有机碳的预测装置中反演执行模块的模块图。
具体实施方式
下面将结合本发明实施例中附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。通常在此处附图中描述和示出的本发明实施例的组件可以以各种不同的配置来布置和设计。因此，以下对在附图中提供的本发明的实施例的详细描述并非旨在限制要求保护的本发明的范围，而是仅仅表示本发明的选定实施例。基于本发明的实施例，本领域技术人员在没有做出创造性劳动的前提下所获得的所有其他实施例，都属于本发明保护的范围。
应注意到：相似的标号和字母在下面的附图中表示类似项，因此，一旦某一项在一个附图中被定义，则在随后的附图中不需要对其进行进一步定义和解释。同时，在本发明的描述中，术语“第一”、“第二”等仅用于区分描述，而不能理解为指示或暗示相对重要性。
图1示出了一种可应用于本发明实施例中的电子设备的结构框图。如图1所示，电子设备100包括存储器102、存储控制器104，一个或多个(图中仅示出一个)处理器106、外设接口108、射频模块110、音频模块112、显示单元114等。这些组件通过一条或多条通讯总线/信号线116相互通讯。
存储器102可用于存储软件程序以及模块，如本发明实施例中的总有机碳的预测方法及装置对应的程序指令/模块，处理器106通过运行存储在存储器102内的软件程序以及模块，从而执行各种功能应用以及数据处理，如本发明实施例提供的总有机碳的预测方法。
存储器102可包括高速随机存储器，还可包括非易失性存储器，如一个或者多个磁性存储装置、闪存、或者其他非易失性固态存储器。处理器106以及其他可能的组件对存储器102的访问可在存储控制器104的控制下进行。
外设接口108将各种输入/输出装置耦合至处理器106以及存储器102。在一些实施例中，外设接口108，处理器106以及存储控制器104可以在单个芯片中实现。在其他一些实例中，他们可以分别由独立的芯片实现。
射频模块110用于接收以及发送电磁波，实现电磁波与电信号的相互转换，从而与通讯网络或者其他设备进行通讯。
音频模块112向用户提供音频接口，其可包括一个或多个麦克风、一个或者多个扬声器以及音频电路。
显示单元114在电子设备100与用户之间提供一个显示界面。具体地，显示单元114向用户显示视频输出，这些视频输出的内容可包括文字、图形、视频及其任意组合。
可以理解，图1所示的结构仅为示意，电子设备100还可包括比图1中所示更多或者更少的组件，或者具有与图1所示不同的配置。图1中所示的各组件可以采用硬件、软件或其组合实现。
第一实施例
如图2示出了本发明实施例提供的总有机碳的预测方法的流程图。请参见图2，该总有机碳的预测方法包括：
步骤S110：基于待测页岩气储层的纵波阻抗曲线以及梯度阻抗曲线拟合出在多个角度的总有机碳指示因子曲线，其中，所述多个角度包括0°～360°范围内的各个角度。
在对待测页岩气储层的总有机碳的含量进行预测时，由于需要考虑梯度阻抗的影响，因此可以先对用于预测总有机碳的总有机碳指示因子曲线进行拟合，并且使总有机碳指示因子曲线中包括梯度阻抗。
在本发明实施例中，在基于待测页岩气储层的纵波阻抗曲线以及梯度阻抗曲线拟合总有机碳指示因子曲线时之前，需要获得该待测页岩气层的纵波阻抗曲线以及梯度阻抗曲线。
因此，在本发明实施例中，步骤S110之前，该总有机碳的预测方法还可以包括：
获取所述待测页岩气储层的测井曲线；基于纵波阻抗计算公式以及所述测井曲线获得所述纵波阻抗曲线；基于梯度阻抗计算公式以及所述测井曲线获得所述梯度阻抗曲线。
在本发明实施例中，待测页岩气储层的测井曲线可以从该页岩气储层的测井数据中获取。该测井曲线可以是用户输入的数据，也可以是其他终端设备传输的数据。
当然，待测页岩气储层的测井曲线的具体获取方式在本发明实施例中并不作为限定。
在本发明实施例中，基于测井曲线中可以获得密度数据、纵波速度、横波速度等数据。
在本发明实施例中，纵波阻抗计算公式可以为：AI＝Vpρ，其中，Vp为纵波速度，ρ为密度，AI为纵波阻抗。
在本发明实施例中，梯度阻抗计算公式可以为：GI＝Vpρ-4kVs -8k，其中，Vp为纵波速度，Vs为纵波速度，ρ为密度，k为k＝(vs/vp)2。
从而，基于待测页岩气储层的测井曲线利用上述纵波阻抗计算公式以及梯度阻抗计算公式，可以获得该待测页岩气储层对应的纵波阻抗曲线以及该待测页岩气储层对应的梯度阻抗曲线。
在获得上述的该待测页岩气储层对应的纵波阻抗曲线以及该待测页岩气储层对应的梯度阻抗曲线之后，在纵波阻抗与梯度阻抗对应的空间中，通过坐标旋转，可以拟合出该待测页岩气储层对应的多个角度的总有机碳指示因子曲线，其表达式可以是：TOCN(λ)＝AIcosλ+GIsinλ，其中，TOCN为拟合的总有机碳指示因子，λ为旋转角度，λ的范围可以为0°～360°。
可以理解的是，上述表达式可以是在λ取不同的值时，对应了一个角度的总有机碳指示因子曲线。
例如，TOCN(30°)＝AIcos30°+GIsin30°表示旋转角度为30°的总有机碳指示因子曲线。
步骤S120：获取所述多个角度的总有机碳指示因子曲线中与所述待测页岩气储层的实测总有机碳曲线的相似度最高的总有机碳指示因子曲线对应的第一角度。
在本发明实施例中，可以基于该待测页岩气储层的实测总有机碳曲线，获取拟合的相似度最高的某个角度的总有机碳指示因子曲线，以使后续进行总有机碳预测的结果的准确度和可靠性高。
具体的，请参见图3，步骤S120可以包括：
步骤S121：获取所述待测页岩气储层的实测总有机碳曲线。
在本发明实施例中，待测页岩气储层的实测总有机碳曲线可以为用户输入，也可以从其他终端设备获取。
当然，在本发明实施例中，待测页岩气储层的实测总有机碳曲线的获取方式在本发明实施例中并不作为限定。
步骤S122：获取所述多个角度的总有机碳指示因子曲线与所述待测页岩气储层的实测总有机碳曲线的夹角余弦作为所述相似度。
在本发明实施例中，可以以实测总有机碳曲线与多个角度的总有机碳指示因子曲线的夹角余弦作为曲线之间的相似度。
具体的，相似度可以为：
Figure BDA0001528464840000081
其中，cosβ为上述的夹角余弦，TOCN为拟合的总有机碳指示因子，TOCO为实测总有机碳。
从而，可以获得多个角度中的每个角度对应的总有机碳指示因子曲线与实测总有机碳曲线的相似度。
步骤S123：获取所述夹角余弦最大时对应的总有机碳指示因子曲线对应的第一角度。
在步骤S122中获得的多个夹角余弦中，可以获得多个夹角余弦中最大的值对应的总有机碳指示因子曲线，并且获得该总有机碳指示因子曲线对应的角度，即旋转角度。
如图4所示，为拟合的各个角度总有机碳指示因子曲线与实测总有机碳曲线相似性关系图。图4中例举的示范井的总有机碳指示因子曲线在在纵波阻抗与梯度阻抗空间中经过115度旋转与实测的总有机碳曲线相似度最高，相似系数为0.828，这个角度的总有机碳指示因子曲线可以理解为最佳总有机碳指示因子曲线。
图5所示，为总有机碳实测曲线与最佳的总有机碳指示因子曲线示意图，可以得出，两条曲线相似性很高。
从而，可以获得多个角度的总有机碳指示因子曲线中与待测页岩气储层的实测总有机碳曲线的相似度最高的总有机碳指示因子曲线对应的第一角度θ。即，经过上述的相似运算获得相似度最高的总有机碳指示因子曲线及其对应的旋转角度θ。
步骤S130：基于扩展弹性阻抗计算公式以及所述待测页岩气储层的叠前地震道集数据进行反演，得到纵波阻抗体以及梯度阻抗体。
在本发明实施例中，需要获取待测页岩气层的三维空间中各个位置的总有机碳预测含量，即获得总有机碳指示因子三维数据，因此，需要获取用于计算的纵波阻抗体数据以及梯度阻抗体数据。
在本发明实施例中，请参见图6，基于扩展弹性阻抗计算公式以及所述待测页岩气储层的叠前地震道集数据进行反演，得到纵波阻抗体以及梯度阻抗体，可以包括：
步骤S131：定义地震波垂直入射对应的纵波阻抗为所述扩展弹性阻抗公式中常数，将所述扩展弹性阻抗公式转换为关于纵波阻抗和梯度阻抗的第一表达式。
在本发明实施例中，扩展弹性阻抗计算公式可以为：
Figure BDA0001528464840000101
其中，p＝(cosx+sinx)，q＝-Ksinx，r＝(cosx-4Ksinx)，EEI(x)为扩展弹性阻抗，α为纵波速度，β为横波速度，ρ为密度，x为角度，k＝(α/β)2，α0、β0、ρ0分别为垂直入射角为0°时对应的纵波速度、横波速度以及密度。
上述扩展弹性阻抗计算公式为Whitcombe进行修正后的归一化方程。
可以理解的是，由于引入了常数α0，β0，ρ0，弹性波阻抗就归一化到了声波阻抗的尺度上，当x＝0°，其值为声波阻抗AI，当x＝90°，其值为梯度阻抗GI。
定义垂直入射时对应的纵波阻抗为常数，即AI0＝α0ρ0。则上述扩展弹性计算公式可以转换为关于纵波阻抗和梯度阻抗的第一表达式：
Figure BDA0001528464840000102
其中，AI为纵波阻抗，GI为梯度阻抗，AI0垂直入射纵波阻抗值，GI0垂直入射梯度阻抗值，x为角度。
步骤S132：基于所述叠前地震道集数据以及所述扩展弹性阻抗公式计算获得扩展弹性阻抗数据。
利用第一表达式对纵波阻抗以及梯度阻抗进行反演时，需要获取扩展弹性阻抗数据。具体的，可以基于叠前地震道集数据中的纵波速度、横波速度、密度、入射角度等，采用上述扩展弹性阻抗公式计算，获得扩展弹性阻抗数据。
步骤S133：基于所述扩展弹性阻抗数据以及所述第一表达式对纵波阻抗以及梯度阻抗进行反演，得到所述纵波阻抗体以及所述梯度阻抗体。
在获得扩展弹性阻抗数据，以及上述的关于纵波阻抗和梯度阻抗的第一表达式之后，可以对纵波阻抗以及梯度阻抗进行反演，从而获得纵波阻抗体以及梯度阻抗体。
如图7所示为反演得到的纵波阻抗剖面，如图8所示为反演得到的梯度阻抗剖面。
步骤S140：将所述第一角度、所述纵波阻抗体以及所述梯度阻抗体代入预设总有机碳指示因子体表达式进行计算从而获得所述待测页岩气储层的总有机碳的预测含量。
基于步骤S120中获得的多个角度的总有机碳指示因子曲线中与待测页岩气储层的实测总有机碳曲线的相似度最高的总有机碳指示因子曲线对应的第一角度θ，可以获得预设总有机碳指示因子体表达式为：TOC＝AIcosθ+GIsinθ，其中，TOC为总有机碳指示因子体，AI为纵波阻抗体，GI为梯度阻抗体，θ为总有机碳指示因子的最佳旋转角度即第一角度θ。
从而将纵波阻抗体以及梯度阻抗体代入预设总有机碳指示因子体表达式可以获得总有机碳指示因子体。如图9所示为获得的总有机碳指示因子体的连井剖面图，基于总有机碳指示因子体可以获得总有机碳含量的定了预测，从而获得该待测页岩气储层的总有机碳的预测含量。
本发明第一实施例提供的总有机碳的预测方法，用叠前地震资料预测页岩气储层的TOC。该方法首先在纵波阻抗与梯度阻抗空间中，拟合出TOC的指示因子曲线，并经过相似运算得到最佳TOC的指示因子曲线及其对应的旋转角度；其次运用修改扩展弹性阻抗反演得到纵波阻抗和梯度阻抗；最后通过纵波阻抗和梯度阻抗实现页岩气储层的TOC预测。从而，该总有机碳的预测方法充分考虑了地震数据的梯度变化对应的梯度阻抗因素，最终准确地实现了页岩气储层的总有机碳的预测。
第二实施例
本发明第二实施例提供了一种总有机碳的预测装置100，请参见图10，所述总有机碳的预测装置100包括第一拟合模块110、角度获取模块120、反演执行模块130以及含量计算模块140。其中，所述第一拟合模块110用于基于待测页岩气储层的纵波阻抗曲线以及梯度阻抗曲线拟合出在多个角度的总有机碳指示因子曲线，其中，所述多个角度包括0°～360°范围内的各个角度；所述角度获取模块120用于获取所述多个角度的总有机碳指示因子曲线中与所述待测页岩气储层的实测总有机碳曲线的相似度最高的总有机碳指示因子曲线对应的第一角度；所述反演执行模块130用于基于扩展弹性阻抗计算公式以及所述待测页岩气储层的叠前地震道集数据进行反演，得到纵波阻抗体以及梯度阻抗体；所述含量计算模块140用于将所述第一角度、所述纵波阻抗体以及所述梯度阻抗体代入预设总有机碳指示因子体表达式进行计算从而获得所述待测页岩气储层的总有机碳的预测含量。
在本发明实施例中，所述总有机碳的预测装置100还包括测井曲线获取模块、第一曲线获取模块以及第二曲线获取模块。其中，所述测井曲线获取模块用于获取所述待测页岩气储层的测井曲线；所述第一曲线获取模块用于基于纵波阻抗计算公式以及所述测井曲线获得所述纵波阻抗曲线；所述第二曲线获取模块用于基于梯度阻抗计算公式以及所述测井曲线获得所述梯度阻抗曲线。
在本发明实施例中，请参见图11，所述角度获取模块120包括实测曲线获取单元121、相似度获取单元122以及角度获取执行单元123。其中，所述实测曲线获取单元121用于获取所述待测页岩气储层的实测总有机碳曲线；所述相似度获取单元122用于获取所述多个角度的总有机碳指示因子曲线与所述待测页岩气储层的实测总有机碳曲线的夹角余弦作为所述相似度；所述角度获取执行单元123获取所述夹角余弦最大时对应的总有机碳指示因子曲线对应的第一角度。
在本发明实施例中，请参见图12，所述反演执行模块130包括表达式获取单元131、阻抗数据计算单元132以及反演单元133。其中，所述表达式获取单元131用于定义地震波垂直入射对应的纵波阻抗为所述扩展弹性阻抗公式中常数，将所述扩展弹性阻抗公式转换为关于纵波阻抗和梯度阻抗的第一表达式；所述阻抗数据计算单元132用于基于所述叠前地震道集数据以及所述扩展弹性阻抗公式计算获得扩展弹性阻抗数据；所述反演单元133用于基于所述扩展弹性阻抗数据以及所述第一表达式对纵波阻抗以及梯度阻抗进行反演，得到所述纵波阻抗体以及所述梯度阻抗体。
第三实施例
本发明第三实施例提供了一种电子设备100，请参见图1，所述电子设备100包括存储器102和处理器106，所述存储器102存储有计算机指令，当所述计算机指令由所述处理器106读取并执行时，使所述处理器106执行如本发明第一实施例提供的总有机碳的预测方法。
第四实施例
本发明第四实施例提供了一种存储介质，所述存储介质中存储有计算机指令，其中，所述计算机指令在被读取并运行时执行如本发明第一实施例提供的总有机碳的预测方法。
综上所述，本发明实施例提供的总有机碳的预测方法、装置、电子设备及存储介质，该方法通过基于待测页岩气储层的纵波阻抗曲线以及梯度阻抗曲线拟合出在多个角度的总有机碳指示因子曲线，其中，多个角度包括0°～360°范围内的各个角度，然后获取多个角度的总有机碳指示因子曲线中与待测页岩气储层的实测总有机碳曲线的相似度最高的总有机碳指示因子曲线对应的第一角度，再基于扩展弹性阻抗计算公式以及待测页岩气储层的叠前地震道集数据进行反演，得到纵波阻抗体以及梯度阻抗体，最后将第一角度、纵波阻抗体以及梯度阻抗体代入预设总有机碳指示因子体表达式进行计算从而获得待测页岩气储层的总有机碳的预测含量。该总有机碳的预测方法考虑了梯度阻抗因素，从而使预测的待测页岩气储层的总有机碳含量的精确度较高。
需要说明的是，本说明书中的各个实施例均采用递进的方式描述，每个实施例重点说明的都是与其他实施例的不同之处，各个实施例之间相同相似的部分互相参见即可。对于装置类实施例而言，由于其与方法实施例基本相似，所以描述的比较简单，相关之处参见方法实施例的部分说明即可。
本申请所提供的几个实施例中，应该理解到，所揭露的装置和方法，也可以通过其它的方式实现。以上所描述的装置实施例仅仅是示意性的，例如，附图中的流程图和框图显示了根据本发明的多个实施例的装置、方法和计算机程序产品的可能实现的体系架构、功能和操作。在这点上，流程图或框图中的每个方框可以代表一个模块、程序段或代码的一部分，所述模块、程序段或代码的一部分包含一个或多个用于实现规定的逻辑功能的可执行指令。也应当注意，在有些作为替换的实现方式中，方框中所标注的功能也可以以不同于附图中所标注的顺序发生。例如，两个连续的方框实际上可以基本并行地执行，它们有时也可以按相反的顺序执行，这依所涉及的功能而定。也要注意的是，框图和/或流程图中的每个方框、以及框图和/或流程图中的方框的组合，可以用执行规定的功能或动作的专用的基于硬件的系统来实现，或者可以用专用硬件与计算机指令的组合来实现。
另外，在本发明各个实施例中的各功能模块可以集成在一起形成一个独立的部分，也可以是各个模块单独存在，也可以两个或两个以上模块集成形成一个独立的部分。
所述功能如果以软件功能模块的形式实现并作为独立的产品销售或使用时，可以存储在一个计算机可读取存储介质中。基于这样的理解，本发明的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的部分可以以软件产品的形式体现出来，该计算机软件产品存储在一个存储介质中，包括若干指令用以使得一台计算机设备(可以是个人计算机，服务器，或者网络设备等)执行本发明各个实施例所述方法的全部或部分步骤。而前述的存储介质包括：U盘、移动硬盘、只读存储器(ROM，Read-Only Memory)、随机存取存储器(RAM，Random Access Memory)、磁碟或者光盘等各种可以存储程序代码的介质。需要说明的是，在本文中，诸如第一和第二等之类的关系术语仅仅用来将一个实体或者操作与另一个实体或操作区分开来，而不一定要求或者暗示这些实体或操作之间存在任何这种实际的关系或者顺序。而且，术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、物品或者设备不仅包括那些要素，而且还包括没有明确列出的其他要素，或者是还包括为这种过程、方法、物品或者设备所固有的要素。在没有更多限制的情况下，由语句“包括一个……”限定的要素，并不排除在包括所述要素的过程、方法、物品或者设备中还存在另外的相同要素。
以上所述仅为本发明的优选实施例而已，并不用于限制本发明，对于本领域的技术人员来说，本发明可以有各种更改和变化。凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。应注意到：相似的标号和字母在下面的附图中表示类似项，因此，一旦某一项在一个附图中被定义，则在随后的附图中不需要对其进行进一步定义和解释。
以上所述，仅为本发明的具体实施方式，但本发明的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本发明揭露的技术范围内，可轻易想到变化或替换，都应涵盖在本发明的保护范围之内。因此，本发明的保护范围应所述以权利要求的保护范围为准。

Claims (10)

1.一种总有机碳的预测方法，其特征在于，所述方法包括：
基于待测页岩气储层的纵波阻抗曲线以及梯度阻抗曲线拟合出在多个角度的总有机碳指示因子曲线，其中，所述多个角度包括0°～360°范围内的各个角度；
获取所述多个角度的总有机碳指示因子曲线中与所述待测页岩气储层的实测总有机碳曲线的相似度最高的总有机碳指示因子曲线对应的第一角度；
基于扩展弹性阻抗计算公式以及所述待测页岩气储层的叠前地震道集数据进行反演，得到纵波阻抗体以及梯度阻抗体；
将所述第一角度、所述纵波阻抗体以及所述梯度阻抗体代入预设总有机碳指示因子体表达式进行计算从而获得所述待测页岩气储层的总有机碳的预测含量；
其中，所述扩展弹性阻抗计算公式为：
Figure FDA0002559361590000011
其中，p＝(cos x+sin x)，q＝-K sin x，r＝(cos x-4K sin x)，EEI(x)为扩展弹性阻抗数据，α为纵波速度，β为横波速度，ρ为密度，x为角度，K＝(α/β)2，α0、β0、ρ0分别为垂直入射角为0°时对应的纵波速度、横波速度以及密度。
2.根据权利要求1所述的方法，其特征在于，所述基于纵波阻抗曲线以及梯度阻抗曲线拟合出在多个角度的总有机碳指示因子曲线之前，所述方法还包括：
获取所述待测页岩气储层的测井曲线；
基于纵波阻抗计算公式以及所述测井曲线获得所述纵波阻抗曲线；
基于梯度阻抗计算公式以及所述测井曲线获得所述梯度阻抗曲线。
3.根据权利要求1所述的方法，其特征在于，所述获取所述多个角度的总有机碳指示因子曲线中与所述待测页岩气储层的实测总有机碳曲线的相似度最高的总有机碳指示因子曲线对应的第一角度，包括：
获取所述待测页岩气储层的实测总有机碳曲线；
获取所述多个角度的总有机碳指示因子曲线与所述待测页岩气储层的实测总有机碳曲线的夹角余弦作为所述相似度；
获取所述夹角余弦最大时对应的总有机碳指示因子曲线对应的第一角度。
4.根据权利要求1所述的方法，其特征在于，所述基于扩展弹性阻抗计算公式以及所述待测页岩气储层的叠前地震道集数据进行反演，得到纵波阻抗体以及梯度阻抗体，包括：
定义地震波垂直入射对应的纵波阻抗为所述扩展弹性阻抗计算公式中常数，将所述扩展弹性阻抗计算公式转换为关于纵波阻抗和梯度阻抗的第一表达式；
基于所述叠前地震道集数据以及所述扩展弹性阻抗计算公式计算获得所述扩展弹性阻抗数据；
基于所述扩展弹性阻抗数据以及所述第一表达式对纵波阻抗以及梯度阻抗进行反演，得到所述纵波阻抗体以及所述梯度阻抗体。
5.一种总有机碳的预测装置，其特征在于，所述装置包括第一拟合模块、角度获取模块、反演执行模块以及含量计算模块，其中，
所述第一拟合模块用于基于待测页岩气储层的纵波阻抗曲线以及梯度阻抗曲线拟合出在多个角度的总有机碳指示因子曲线，其中，所述多个角度包括0°～360°范围内的各个角度；
所述角度获取模块用于获取所述多个角度的总有机碳指示因子曲线中与所述待测页岩气储层的实测总有机碳曲线的相似度最高的总有机碳指示因子曲线对应的第一角度；
所述反演执行模块用于基于扩展弹性阻抗计算公式以及所述待测页岩气储层的叠前地震道集数据进行反演，得到纵波阻抗体以及梯度阻抗体；
所述含量计算模块用于将所述第一角度、所述纵波阻抗体以及所述梯度阻抗体代入预设总有机碳指示因子体表达式进行计算从而获得所述待测页岩气储层的总有机碳的预测含量；
其中，所述扩展弹性阻抗计算公式为：
Figure FDA0002559361590000031
其中，p＝(cos x+sin x)，q＝-K sin x，r＝(cos x-4K sin x)，EEI(x)为扩展弹性阻抗数据，α为纵波速度，β为横波速度，ρ为密度，x为角度，K＝(α/β)2，α0、β0、ρ0分别为垂直入射角为0°时对应的纵波速度、横波速度以及密度。
6.根据权利要求5所述的装置，其特征在于，所述装置还包括测井曲线获取模块、第一曲线获取模块以及第二曲线获取模块，其中，
所述测井曲线获取模块用于获取所述待测页岩气储层的测井曲线；
所述第一曲线获取模块用于基于纵波阻抗计算公式以及所述测井曲线获得所述纵波阻抗曲线；
所述第二曲线获取模块用于基于梯度阻抗计算公式以及所述测井曲线获得所述梯度阻抗曲线。
7.根据权利要求5所述的装置，其特征在于，所述角度获取模块包括实测曲线获取单元、相似度获取单元以及角度获取执行单元，其中，
所述实测曲线获取单元用于获取所述待测页岩气储层的实测总有机碳曲线；
所述相似度获取单元用于获取所述多个角度的总有机碳指示因子曲线与所述待测页岩气储层的实测总有机碳曲线的夹角余弦作为所述相似度；
所述角度获取执行单元获取所述夹角余弦最大时对应的总有机碳指示因子曲线对应的第一角度。
8.根据权利要求5所述的装置，其特征在于，所述反演执行模块包括表达式获取单元、阻抗数据计算单元以及反演单元，其中，
所述表达式获取单元用于定义地震波垂直入射对应的纵波阻抗为所述扩展弹性阻抗计算公式中常数，将所述扩展弹性阻抗计算公式转换为关于纵波阻抗和梯度阻抗的第一表达式；
所述阻抗数据计算单元用于基于所述叠前地震道集数据以及所述扩展弹性阻抗计算公式计算获得所述扩展弹性阻抗数据；
所述反演单元用于基于所述扩展弹性阻抗数据以及所述第一表达式对纵波阻抗以及梯度阻抗进行反演，得到所述纵波阻抗体以及所述梯度阻抗体。
9.一种电子设备，其特征在于，所述电子设备包括存储器和处理器，所述存储器存储有计算机指令，当所述计算机指令由所述处理器读取并执行时，使所述处理器执行如权利要求1-4中任一权项所述的方法。
10.一种存储介质，其特征在于，所述存储介质中存储有计算机指令，其中，所述计算机指令在被读取并运行时执行如权利要求1-4中任一权项所述的方法。
CN201711450913.3A 2017-12-27 2017-12-27 总有机碳的预测方法、装置、电子设备及存储介质 Active CN108171376B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201711450913.3A CN108171376B (zh) 2017-12-27 2017-12-27 总有机碳的预测方法、装置、电子设备及存储介质

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201711450913.3A CN108171376B (zh) 2017-12-27 2017-12-27 总有机碳的预测方法、装置、电子设备及存储介质

Publications (2)

Publication Number Publication Date
CN108171376A CN108171376A (zh) 2018-06-15
CN108171376B true CN108171376B (zh) 2020-09-18

Family

ID=62518852

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201711450913.3A Active CN108171376B (zh) 2017-12-27 2017-12-27 总有机碳的预测方法、装置、电子设备及存储介质

Country Status (1)

Country Link
CN (1) CN108171376B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US11880776B2 (en) 2021-11-25 2024-01-23 Institute Of Geology And Geophysics, Chinese Academy Of Sciences Graph neural network (GNN)-based prediction system for total organic carbon (TOC) in shale

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113740911B (zh) * 2021-09-06 2023-09-26 北京海润联创石油科技有限公司 一种基于坐标旋转波阻抗反演提高储层预测精度的方法

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104977618B (zh) * 2014-04-09 2017-09-22 中国石油集团东方地球物理勘探有限责任公司 一种评价页岩气储层及寻找甜点区的方法
CN105277982B (zh) * 2014-07-25 2016-08-31 中国石油化工股份有限公司 一种泥页岩总有机碳含量地震预测方法
CN104345345B (zh) * 2014-10-30 2017-09-01 中国石油天然气集团公司 一种页岩储层总有机碳toc含量预测方法
CN106324674B (zh) * 2016-08-23 2018-06-12 中国石油大学(华东) 一种页岩气toc叠前地震反演预测方法
CN106526669B (zh) * 2016-09-19 2018-10-23 中国石油化工股份有限公司 一种页岩油气藏的地震储层预测方法

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US11880776B2 (en) 2021-11-25 2024-01-23 Institute Of Geology And Geophysics, Chinese Academy Of Sciences Graph neural network (GNN)-based prediction system for total organic carbon (TOC) in shale

Also Published As

Publication number Publication date
CN108171376A (zh) 2018-06-15

Similar Documents

Publication Publication Date Title
CN107679358B (zh) 一种确定储层渗透率的方法及装置
Knetter et al. Four‐point discontinuity observations using Cluster magnetic field data: A statistical survey
CN106646611B (zh) 一种纵波和转换波地震数据时间匹配的方法及装置
CN107092031B (zh) 一种声波测井数据校正方法及装置
CN108171376B (zh) 总有机碳的预测方法、装置、电子设备及存储介质
CN110794458B (zh) 基于时频分析的含气性检测方法、装置及存储介质
CN110068502B (zh) 砾岩强度确定方法及装置
CN106324665A (zh) 反演裂缝密度的方法及系统
CN107829727A (zh) 一种确定地应力的方法及装置
Lyu et al. Efficiency of the spectral element method with very high polynomial degree to solve the elastic wave equation
Poulos et al. Evaluation of earthquake response spectra directionality using stochastic simulations
WO2022159862A1 (en) Two methods of determining permeabilities of naturally fractured rocks from laboratory measurements
Li et al. Joint inversion of formation radial shear-velocity profiles by dipole acoustic logging while drilling
CN106772599A (zh) 一种计算地层横波速度的方法及装置
Mignan et al. The spatial scale of detected seismicity
Andersen et al. Machine learning based prediction of porosity and water saturation from Varg field reservoir well logs
CN109444061B (zh) 过敏原溶液浓度检测方法、装置、设备和存储介质
Almeida et al. Improving parameter priors for data‐scarce estimation problems
US11327188B2 (en) Robust arrival picking of seismic vibratory waves
Chunxi et al. Formation permeability estimation using Stoneley waves from logging while drilling: theory, method, and application
Zeng et al. Quaternion-based anisotropic inversion for flexural waves in horizontal transverse isotropic formations with unmatched sources: A synthetic example
Xue et al. Application of variational mode decomposition–based Hilbert marginal differential cepstrum for hydrocarbon detection
Thombare et al. Determining P-wave 45 laboratory acoustic measurement using a vertically drilled plug
US7861432B1 (en) Method and apparatus for hole measurement
CN118094339B (zh) 地层温度预测方法、装置及计算设备

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
TA01 Transfer of patent application right

Effective date of registration: 20200414

Address after: 100027 Beijing, Chaoyangmen, North Street, No. 22, No.

Applicant after: China Petroleum & Chemical Corp.

Applicant after: EXPLORATION AND DEVELOPMENT RESEARCH INSTITUTE OF SINOPEC JIANGHAN OILFIELD Co.

Address before: 430000 18 University Garden Road, East Lake hi tech Development Zone, Wuhan, Hubei (China Science and Technology Park)

Applicant before: EXPLORATION AND DEVELOPMENT RESEARCH INSTITUTE OF SINOPEC JIANGHAN OILFIELD Co.

TA01 Transfer of patent application right
GR01 Patent grant
GR01 Patent grant