CN111898785B - 基于页岩的断裂韧性空间分布特征预测方法及系统 - Google Patents
基于页岩的断裂韧性空间分布特征预测方法及系统 Download PDFInfo
- Publication number
- CN111898785B CN111898785B CN201910367277.0A CN201910367277A CN111898785B CN 111898785 B CN111898785 B CN 111898785B CN 201910367277 A CN201910367277 A CN 201910367277A CN 111898785 B CN111898785 B CN 111898785B
- Authority
- CN
- China
- Prior art keywords
- fracture toughness
- shale
- parameters
- longitudinal
- transverse 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 23
- 238000004364 calculation method Methods 0.000 claims abstract description 77
- 230000000704 physical effect Effects 0.000 claims abstract description 24
- 230000014509 gene expression Effects 0.000 claims abstract description 22
- 238000000611 regression analysis Methods 0.000 claims abstract description 20
- 238000004590 computer program Methods 0.000 claims description 12
- 238000012512 characterization method Methods 0.000 claims description 11
- 238000005259 measurement Methods 0.000 claims description 5
- 238000007656 fracture toughness test Methods 0.000 claims description 4
- 235000021185 dessert Nutrition 0.000 description 11
- 238000011161 development Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000011156 evaluation Methods 0.000 description 2
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 229910052799 carbon Inorganic materials 0.000 description 1
- 238000011545 laboratory measurement Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000011435 rock Substances 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression 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
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Business, Economics & Management (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Human Resources & Organizations (AREA)
- Economics (AREA)
- Data Mining & Analysis (AREA)
- Strategic Management (AREA)
- Tourism & Hospitality (AREA)
- Mathematical Optimization (AREA)
- Operations Research (AREA)
- Marketing (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Physics (AREA)
- Computational Mathematics (AREA)
- General Business, Economics & Management (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Engineering & Computer Science (AREA)
- Algebra (AREA)
- Animal Husbandry (AREA)
- Marine Sciences & Fisheries (AREA)
- Mining & Mineral Resources (AREA)
- Software Systems (AREA)
- Health & Medical Sciences (AREA)
- Databases & Information Systems (AREA)
- General Health & Medical Sciences (AREA)
- Agronomy & Crop Science (AREA)
- Probability & Statistics with Applications (AREA)
- Primary Health Care (AREA)
- Evolutionary Biology (AREA)
- Bioinformatics & Computational Biology (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Development Economics (AREA)
- Game Theory and Decision Science (AREA)
- Entrepreneurship & Innovation (AREA)
- Quality & Reliability (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
本发明提供一种基于页岩的断裂韧性空间分布特征预测方法及系统。该基于页岩的断裂韧性空间分布特征预测方法包括：获取页岩的基础物性参数、页岩的纵横波属性参数和页岩的断裂韧性值；根据基础物性参数、纵横波属性参数和断裂韧性值建立多组断裂韧性值与几何参数，纵横波属性参数的关系式；对多组断裂韧性值与几何参数，纵横波属性参数的关系式进行多元回归分析，得到断裂韧性计算模型；输入测井参数至断裂韧性计算模型中，得到测井的断裂韧性值；输入地震数据体至断裂韧性计算模型中，得到地震的断裂韧性值；根据测井的断裂韧性值和地震的断裂韧性值得到准确的断裂韧性空间分布特征，进一步确定页岩气的分布范围。
Description
技术领域
本发明涉及页岩气技术领域，具体地，涉及一种基于页岩的断裂韧性空间分布特征预测方法及系统。
背景技术
页岩气甜点关键参数预测包括TOC(Total Organic Carbon，总有机碳)含量、含气性、孔隙度、脆性指数等参数，根据甜点参数及地质解释可以共同确定页岩气的分布，脆性指数则作为后期压裂开发页岩气的重要参数。但是脆性指数这个参数不具备实际的物理含义，也无法通过对岩石等材料进行实验室测得，在页岩气脆性指数预测中，一般的做法是根据其他弹性参数计算拟合而得到。然而与脆性指数相对应的材料学中的一个参数——断裂韧性参数则具备实际的物理含义，因此断裂韧性参数给页岩气的脆性和可压性评价带来了实际的意义，也有望实现脆性和可压性评价的准确预测。现有技术无法得到准确的断裂韧性空间分布特征，因此也无法确定页岩气的分布范围。
发明内容
本发明实施例的主要目的在于提供一种基于页岩的断裂韧性空间分布特征预测方法及系统，以得到准确的断裂韧性空间分布特征，实现对页岩气甜点参数的有效预测，进一步确定页岩气的分布范围。
为了实现上述目的，本发明实施例提供一种基于页岩的断裂韧性空间分布特征预测方法，包括：
获取页岩的基础物性参数、页岩的纵横波属性参数和页岩的断裂韧性值；
根据基础物性参数、纵横波属性参数和断裂韧性值建立多组断裂韧性值与几何参数，纵横波属性参数的关系式；
对多组断裂韧性值与几何参数，纵横波属性参数的关系式进行多元回归分析，得到断裂韧性计算模型；
输入测井参数至断裂韧性计算模型中，得到测井的断裂韧性值；
输入地震数据体至断裂韧性计算模型中，得到地震的断裂韧性值；
根据测井的断裂韧性值和地震的断裂韧性值得到断裂韧性空间分布特征。
本发明实施例还提供一种基于页岩的断裂韧性空间分布特征预测系统，包括：
获取单元，用于获取页岩的基础物性参数、页岩的纵横波属性参数和页岩的断裂韧性值；
关系式建立单元，用于根据基础物性参数、纵横波属性参数和断裂韧性值建立多组断裂韧性值与几何参数，纵横波属性参数的关系式；
断裂韧性计算模型单元，用于对多组断裂韧性值与几何参数，纵横波属性参数的关系式进行多元回归分析，得到断裂韧性计算模型；
第一计算单元，用于输入测井参数至断裂韧性计算模型中，得到测井的断裂韧性值；
第二计算单元，用于输入地震数据体至断裂韧性计算模型中，得到地震的断裂韧性值；
断裂韧性空间分布特征单元，用于根据测井的断裂韧性值和地震的断裂韧性值得到断裂韧性空间分布特征。
本发明实施例还提供一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，处理器执行计算机程序时实现以下步骤：
获取页岩的基础物性参数、页岩的纵横波属性参数和页岩的断裂韧性值；
根据基础物性参数、纵横波属性参数和断裂韧性值建立多组断裂韧性值与几何参数，纵横波属性参数的关系式；
对多组断裂韧性值与几何参数，纵横波属性参数的关系式进行多元回归分析，得到断裂韧性计算模型；
输入测井参数至断裂韧性计算模型中，得到测井的断裂韧性值；
输入地震数据体至断裂韧性计算模型中，得到地震的断裂韧性值；
根据测井的断裂韧性值和地震的断裂韧性值得到断裂韧性空间分布特征。
本发明实施例还提供一种计算机可读存储介质，其上存储有计算机程序，计算机程序被处理器执行时实现以下步骤：
获取页岩的基础物性参数、页岩的纵横波属性参数和页岩的断裂韧性值；
根据基础物性参数、纵横波属性参数和断裂韧性值建立多组断裂韧性值与几何参数，纵横波属性参数的关系式；
对多组断裂韧性值与几何参数，纵横波属性参数的关系式进行多元回归分析，得到断裂韧性计算模型；
输入测井参数至断裂韧性计算模型中，得到测井的断裂韧性值；
输入地震数据体至断裂韧性计算模型中，得到地震的断裂韧性值；
根据测井的断裂韧性值和地震的断裂韧性值得到断裂韧性空间分布特征。
本发明实施例的基于页岩的断裂韧性空间分布特征预测方法及系统，先根据基础物性参数、纵横波属性参数和断裂韧性值建立多组断裂韧性值与几何参数，纵横波属性参数的关系式，再对多组断裂韧性值与几何参数，纵横波属性参数的关系式进行多元回归分析，得到断裂韧性计算模型，然后输入测井参数至断裂韧性计算模型中，得到测井的断裂韧性值，输入地震数据体至断裂韧性计算模型中，得到地震的断裂韧性值，最后根据测井的断裂韧性值和地震的断裂韧性值得到准确的断裂韧性空间分布特征，实现对页岩气甜点参数的有效预测，进一步确定页岩气的分布范围。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1是本发明实施例中基于页岩的断裂韧性空间分布特征预测方法的流程图；
图2是本发明实施例中基于页岩的断裂韧性空间分布特征预测系统的结构框图。
具体实施方式
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
本领域技术技术人员知道，本发明的实施方式可以实现为一种系统、装置、设备、方法或计算机程序产品。因此，本公开可以具体实现为以下形式，即：完全的硬件、完全的软件(包括固件、驻留软件、微代码等)，或者硬件和软件结合的形式。
鉴于现有技术无法得到准确的断裂韧性空间分布特征和页岩气的分布范围，本发明实施例提供一种基于页岩的断裂韧性空间分布特征预测方法，以得到准确的断裂韧性空间分布特征，实现对页岩气甜点参数的有效预测，进一步确定页岩气的分布范围。以下结合附图对本发明进行详细说明。
图1是本发明实施例中基于页岩的断裂韧性空间分布特征预测方法的流程图。如图1所示，基于页岩的断裂韧性空间分布特征预测方法包括：
S101：获取页岩的基础物性参数、页岩的纵横波属性参数和页岩的断裂韧性值。
S102：根据基础物性参数、纵横波属性参数和断裂韧性值建立多组断裂韧性值与几何参数，纵横波属性参数的关系式。
S103：对多组断裂韧性值与几何参数，纵横波属性参数的关系式进行多元回归分析，得到断裂韧性计算模型。
具体实施时，可以分别计算根据每组关系式得到的断裂韧性值与实际断裂韧性值之间的相关系数。最大的相关系数对应的关系式为断裂韧性计算模型。
S104：输入测井参数至断裂韧性计算模型中，得到测井的断裂韧性值。
其中，测井参数包括偶极子声波测井数据或全波列声波测井数据，如测井的横波速度和测井的纵波速度。
S105：输入地震数据体至断裂韧性计算模型中，得到地震的断裂韧性值。
执行S105之前，可以先获取地震数据，再对地震数据进行地震反演以得到地震数据体。地震数据体包括地震的横波速度和地震的纵波速度。
S106：根据测井的断裂韧性值和地震的断裂韧性值得到断裂韧性空间分布特征。
S106包括：先根据测井的断裂韧性值和地震的断裂韧性值得到断裂韧性数据体，再根据断裂韧性数据体获得断裂韧性空间分布特征。
图1所示的基于页岩的断裂韧性空间分布特征预测方法的执行主体可以为计算机。由图1所示的流程可知，本发明实施例的基于页岩的断裂韧性空间分布特征预测方法先根据基础物性参数、纵横波属性参数和断裂韧性值建立多组断裂韧性值与几何参数，纵横波属性参数的关系式，再对多组断裂韧性值与几何参数，纵横波属性参数的关系式进行多元回归分析，得到断裂韧性计算模型，然后输入测井参数至断裂韧性计算模型中，得到测井的断裂韧性值，输入地震数据体至断裂韧性计算模型中，得到地震的断裂韧性值，最后根据测井的断裂韧性值和地震的断裂韧性值得到准确的断裂韧性空间分布特征，实现对页岩气甜点参数的有效预测，进一步确定页岩气的分布范围。
一实施例中，断裂韧性计算模型如下：
KIC＝4.282×10-5Vp+6.256×10-4Vs-0.8091；
其中，KIC为断裂韧性值，Vp为纵波速度，Vs为横波速度。
一实施例中，执行S101之前，还包括：测量页岩式样的几何参数，根据页岩式样的几何参数获取页岩的基础物性参数；对页岩式样进行声波测量，获取页岩的纵横波属性参数，对页岩式样进行断裂韧性测试，获取页岩的断裂韧性值。
具体实施时，可以根据从矿场取回的页岩样品制备人字形切槽巴西圆盘试样作为页岩式样，在测量页岩式样的几何参数后获取页岩的基础物性参数。其中，基础物性参数包括体积和密度，纵横波属性参数包括横波速度、纵波速度、衰减系数和波阻抗。
本发明的其中一个具体实施例如下：
1、根据从矿场取回的页岩样品制备人字形切槽巴西圆盘试样作为页岩式样。
2、测量页岩式样的几何参数，并根据页岩式样的几何参数获取页岩的基础物性参数。
3、对页岩式样进行声波测量，获取页岩的纵横波属性参数；对页岩式样进行断裂韧性测试，获取页岩的断裂韧性值。
4、根据基础物性参数、纵横波属性参数和断裂韧性值建立多组断裂韧性值与几何参数，纵横波属性参数的关系式。
5、对多组断裂韧性值与几何参数，纵横波属性参数的关系式进行多元回归分析，得到断裂韧性计算模型。例如，可以分别计算根据每组关系式得到的断裂韧性值与实际断裂韧性值之间的相关系数，最大的相关系数对应的关系式即为断裂韧性计算模型。
6、获取地震数据，对地震数据进行地震反演以得到地震数据体。
7、输入测井参数至断裂韧性计算模型中，得到测井的断裂韧性值；输入地震数据体至断裂韧性计算模型中，得到地震的断裂韧性值。
8、先根据测井的断裂韧性值和地震的断裂韧性值得到断裂韧性数据体，再根据断裂韧性数据体获得断裂韧性空间分布特征。
本发明可以应用在实际的页岩气地震工区的甜点预测中。从单井预测结果可以看到，本发明预测得到的断裂韧性值和脆性指数相对呈负相关关系，与实验室测量及模拟分析结果保持一致。在储层最好的发育层段可以看出断裂韧性值偏低，而脆性指数偏高的特征，符合地质认识。另外从井震联合断裂韧性预测结果也可以看到，本发明预测得到的断裂韧性值的特征与脆性指数特征呈相对负相关性，符合该区域的页岩地质认识。
综上，本发明实施例的基于页岩的断裂韧性空间分布特征预测方法先根据基础物性参数、纵横波属性参数和断裂韧性值建立多组断裂韧性值与几何参数，纵横波属性参数的关系式，再对多组断裂韧性值与几何参数，纵横波属性参数的关系式进行多元回归分析，得到断裂韧性计算模型，然后输入测井参数至断裂韧性计算模型中，得到测井的断裂韧性值，输入地震数据体至断裂韧性计算模型中，得到地震的断裂韧性值，最后根据测井的断裂韧性值和地震的断裂韧性值得到准确的断裂韧性空间分布特征，实现对页岩气甜点参数的有效预测，进一步确定页岩气的分布范围。
基于同一发明构思，本发明实施例还提供了一种基于页岩的断裂韧性空间分布特征预测系统，由于该系统解决问题的原理与基于页岩的断裂韧性空间分布特征预测方法相似，因此该系统的实施可以参见方法的实施，重复之处不再赘述。
图2是本发明实施例中基于页岩的断裂韧性空间分布特征预测系统的结构框图。如图2所示，基于页岩的断裂韧性空间分布特征预测系统包括：
获取单元，用于获取页岩的基础物性参数、页岩的纵横波属性参数和页岩的断裂韧性值；
关系式建立单元，用于根据基础物性参数、纵横波属性参数和断裂韧性值建立多组断裂韧性值与几何参数，纵横波属性参数的关系式；
断裂韧性计算模型单元，用于对多组断裂韧性值与几何参数，纵横波属性参数的关系式进行多元回归分析，得到断裂韧性计算模型；
第一计算单元，用于输入测井参数至断裂韧性计算模型中，得到测井的断裂韧性值；
第二计算单元，用于输入地震数据体至断裂韧性计算模型中，得到地震的断裂韧性值；
断裂韧性空间分布特征单元，用于根据测井的断裂韧性值和地震的断裂韧性值得到断裂韧性空间分布特征。
在其中一种实施例中，基础物性参数包括体积和密度；
纵横波属性参数包括横波速度、纵波速度、衰减系数和波阻抗；
测井参数包括测井的横波速度和测井的纵波速度；
地震数据体包括地震的横波速度和地震的纵波速度。
在其中一种实施例中，断裂韧性计算模型如下：
KIC＝4.282×10-5Vp+6.256×10-4Vs-0.8091；
其中，KIC为断裂韧性值，Vp为纵波速度，Vs为横波速度。
在其中一种实施例中，获取单元具体用于：测量页岩式样的几何参数；
根据页岩式样的几何参数获取页岩的基础物性参数；
对页岩式样进行声波测量，获取页岩的纵横波属性参数；
对页岩式样进行断裂韧性测试，获取页岩的断裂韧性值。
综上，本发明实施例的基于页岩的断裂韧性空间分布特征预测系统先根据基础物性参数、纵横波属性参数和断裂韧性值建立多组断裂韧性值与几何参数，纵横波属性参数的关系式，再对多组断裂韧性值与几何参数，纵横波属性参数的关系式进行多元回归分析，得到断裂韧性计算模型，然后输入测井参数至断裂韧性计算模型中，得到测井的断裂韧性值，输入地震数据体至断裂韧性计算模型中，得到地震的断裂韧性值，最后根据测井的断裂韧性值和地震的断裂韧性值得到准确的断裂韧性空间分布特征，实现对页岩气甜点参数的有效预测，进一步确定页岩气的分布范围。
本发明实施例还提供了一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，处理器执行计算机程序时实现以下步骤：
获取页岩的基础物性参数、页岩的纵横波属性参数和页岩的断裂韧性值；
根据基础物性参数、纵横波属性参数和断裂韧性值建立多组断裂韧性值与几何参数，纵横波属性参数的关系式；
对多组断裂韧性值与几何参数，纵横波属性参数的关系式进行多元回归分析，得到断裂韧性计算模型；
输入测井参数至断裂韧性计算模型中，得到测井的断裂韧性值；
输入地震数据体至断裂韧性计算模型中，得到地震的断裂韧性值；
根据测井的断裂韧性值和地震的断裂韧性值得到断裂韧性空间分布特征。
综上，本发明实施例的计算机设备先根据基础物性参数、纵横波属性参数和断裂韧性值建立多组断裂韧性值与几何参数，纵横波属性参数的关系式，再对多组断裂韧性值与几何参数，纵横波属性参数的关系式进行多元回归分析，得到断裂韧性计算模型，然后输入测井参数至断裂韧性计算模型中，得到测井的断裂韧性值，输入地震数据体至断裂韧性计算模型中，得到地震的断裂韧性值，最后根据测井的断裂韧性值和地震的断裂韧性值得到准确的断裂韧性空间分布特征，实现对页岩气甜点参数的有效预测，进一步确定页岩气的分布范围。
本发明实施例还提供了一种计算机可读存储介质，其上存储有计算机程序，计算机程序被处理器执行时实现以下步骤：
获取页岩的基础物性参数、页岩的纵横波属性参数和页岩的断裂韧性值；
根据基础物性参数、纵横波属性参数和断裂韧性值建立多组断裂韧性值与几何参数，纵横波属性参数的关系式；
对多组断裂韧性值与几何参数，纵横波属性参数的关系式进行多元回归分析，得到断裂韧性计算模型；
输入测井参数至断裂韧性计算模型中，得到测井的断裂韧性值；
输入地震数据体至断裂韧性计算模型中，得到地震的断裂韧性值；
根据测井的断裂韧性值和地震的断裂韧性值得到断裂韧性空间分布特征。
综上，本发明实施例的计算机可读存储介质先根据基础物性参数、纵横波属性参数和断裂韧性值建立多组断裂韧性值与几何参数，纵横波属性参数的关系式，再对多组断裂韧性值与几何参数，纵横波属性参数的关系式进行多元回归分析，得到断裂韧性计算模型，然后输入测井参数至断裂韧性计算模型中，得到测井的断裂韧性值，输入地震数据体至断裂韧性计算模型中，得到地震的断裂韧性值，最后根据测井的断裂韧性值和地震的断裂韧性值得到准确的断裂韧性空间分布特征，实现对页岩气甜点参数的有效预测，进一步确定页岩气的分布范围。
以上所述的具体实施例，对本发明的目的、技术方案和有益效果进行了进一步详细说明，所应理解的是，以上所述仅为本发明的具体实施例而已，并不用于限定本发明的保护范围，凡在本发明的精神和原则之内，所做的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。
Claims (8)
1.一种基于页岩的断裂韧性空间分布特征预测方法，其特征在于，包括：
获取页岩的基础物性参数、页岩的纵横波属性参数和页岩的断裂韧性值；
根据所述基础物性参数、所述纵横波属性参数和所述断裂韧性值建立多组断裂韧性值与几何参数，纵横波属性参数的关系式；
对所述多组断裂韧性值与几何参数，纵横波属性参数的关系式进行多元回归分析，得到断裂韧性计算模型；
输入测井参数至所述断裂韧性计算模型中，得到测井的断裂韧性值；
输入地震数据体至所述断裂韧性计算模型中，得到地震的断裂韧性值；
根据所述测井的断裂韧性值和所述地震的断裂韧性值得到断裂韧性空间分布特征；
所述断裂韧性计算模型如下：
KIC＝4.282′10-5Vp+6.256′10-4Vs-0.8091；
其中，KIC为断裂韧性值，Vp为纵波速度，Vs为横波速度。
2.根据权利要求1所述的基于页岩的断裂韧性空间分布特征预测方法，其特征在于，
所述基础物性参数包括体积和密度；
所述纵横波属性参数包括横波速度、纵波速度、衰减系数和波阻抗；
所述测井参数包括测井的横波速度和测井的纵波速度；
所述地震数据体包括地震的横波速度和地震的纵波速度。
3.根据权利要求1所述的基于页岩的断裂韧性空间分布特征预测方法，其特征在于，获取页岩的基础物性参数之前，还包括：
测量页岩式样的几何参数；
根据所述页岩式样的几何参数获取页岩的基础物性参数；
对所述页岩式样进行声波测量，获取页岩的纵横波属性参数；
对所述页岩式样进行断裂韧性测试，获取页岩的断裂韧性值。
4.一种基于页岩的断裂韧性空间分布特征预测系统，其特征在于，包括：
获取单元，用于获取页岩的基础物性参数、页岩的纵横波属性参数和页岩的断裂韧性值；
关系式建立单元，用于根据所述基础物性参数、所述纵横波属性参数和所述断裂韧性值建立多组断裂韧性值与几何参数，纵横波属性参数的关系式；
断裂韧性计算模型单元，用于对所述多组断裂韧性值与几何参数，纵横波属性参数的关系式进行多元回归分析，得到断裂韧性计算模型；
第一计算单元，用于输入测井参数至所述断裂韧性计算模型中，得到测井的断裂韧性值；
第二计算单元，用于输入地震数据体至所述断裂韧性计算模型中，得到地震的断裂韧性值；
断裂韧性空间分布特征单元，用于根据所述测井的断裂韧性值和所述地震的断裂韧性值得到断裂韧性空间分布特征；
所述断裂韧性计算模型如下：
KIC＝4.282′10-5Vp+6.256′10-4Vs-0.8091；
其中，KIC为断裂韧性值，Vp为纵波速度，Vs为横波速度。
5.根据权利要求4所述的基于页岩的断裂韧性空间分布特征预测系统，其特征在于，
所述基础物性参数包括体积和密度；
所述纵横波属性参数包括横波速度、纵波速度、衰减系数和波阻抗；
所述测井参数包括测井的横波速度和测井的纵波速度；
所述地震数据体包括地震的横波速度和地震的纵波速度。
6.根据权利要求4所述的基于页岩的断裂韧性空间分布特征预测系统，其特征在于，所述获取单元具体用于：测量页岩式样的几何参数；
根据所述页岩式样的几何参数获取页岩的基础物性参数；
对所述页岩式样进行声波测量，获取页岩的纵横波属性参数；
对所述页岩式样进行断裂韧性测试，获取页岩的断裂韧性值。
7.一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，其特征在于，所述处理器执行所述计算机程序时实现以下步骤：
获取页岩的基础物性参数、页岩的纵横波属性参数和页岩的断裂韧性值；
根据所述基础物性参数、所述纵横波属性参数和所述断裂韧性值建立多组断裂韧性值与几何参数，纵横波属性参数的关系式；
对所述多组断裂韧性值与几何参数，纵横波属性参数的关系式进行多元回归分析，得到断裂韧性计算模型；
输入测井参数至所述断裂韧性计算模型中，得到测井的断裂韧性值；
输入地震数据体至所述断裂韧性计算模型中，得到地震的断裂韧性值；
根据所述测井的断裂韧性值和所述地震的断裂韧性值得到断裂韧性空间分布特征；
所述断裂韧性计算模型如下：
KIC＝4.282′10-5Vp+6.256′10-4Vs-0.8091；
其中，KIC为断裂韧性值，Vp为纵波速度，Vs为横波速度。
8.一种计算机可读存储介质，其上存储有计算机程序，其特征在于，所述计算机程序被处理器执行时实现以下步骤：
获取页岩的基础物性参数、页岩的纵横波属性参数和页岩的断裂韧性值；
根据所述基础物性参数、所述纵横波属性参数和所述断裂韧性值建立多组断裂韧性值与几何参数，纵横波属性参数的关系式；
对所述多组断裂韧性值与几何参数，纵横波属性参数的关系式进行多元回归分析，得到断裂韧性计算模型；
输入测井参数至所述断裂韧性计算模型中，得到测井的断裂韧性值；
输入地震数据体至所述断裂韧性计算模型中，得到地震的断裂韧性值；
根据所述测井的断裂韧性值和所述地震的断裂韧性值得到断裂韧性空间分布特征；
所述断裂韧性计算模型如下：
KIC＝4.282′10-5Vp+6.256′10-4Vs-0.8091；
其中，KIC为断裂韧性值，Vp为纵波速度，Vs为横波速度。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910367277.0A CN111898785B (zh) | 2019-05-05 | 2019-05-05 | 基于页岩的断裂韧性空间分布特征预测方法及系统 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910367277.0A CN111898785B (zh) | 2019-05-05 | 2019-05-05 | 基于页岩的断裂韧性空间分布特征预测方法及系统 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN111898785A CN111898785A (zh) | 2020-11-06 |
| CN111898785B true CN111898785B (zh) | 2024-08-02 |
Family
ID=73169179
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910367277.0A Active CN111898785B (zh) | 2019-05-05 | 2019-05-05 | 基于页岩的断裂韧性空间分布特征预测方法及系统 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN111898785B (zh) |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109577970A (zh) * | 2018-12-11 | 2019-04-05 | 西南石油大学 | 一种页岩储层ⅰ型断裂韧性的测井评价方法 |
Family Cites Families (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP4494884B2 (ja) * | 2004-06-30 | 2010-06-30 | 三菱スペース・ソフトウエア株式会社 | 地震発生情報急報装置及び地震発生情報急報方法 |
| CN104977618B (zh) * | 2014-04-09 | 2017-09-22 | 中国石油集团东方地球物理勘探有限责任公司 | 一种评价页岩气储层及寻找甜点区的方法 |
| CN105221141B (zh) * | 2014-06-23 | 2016-08-31 | 中国石油化工股份有限公司 | 一种泥页岩脆性指数预测方法 |
| WO2017020267A1 (zh) * | 2015-08-05 | 2017-02-09 | 深圳朝伟达科技有限公司 | 一种水力压裂微地震事件的解释方法 |
| CN106204303A (zh) * | 2016-07-08 | 2016-12-07 | 西安石油大学 | 一种基于权重分配的页岩气储层可压性评价方法 |
| CN107063905A (zh) * | 2017-03-29 | 2017-08-18 | 西南石油大学 | 一种基于划痕实验的页岩断裂韧性计算方法 |
| CN109543352B (zh) * | 2018-12-29 | 2023-05-09 | 西南石油大学 | 页岩断裂韧性预测方法、装置及电子设备 |
-
2019
- 2019-05-05 CN CN201910367277.0A patent/CN111898785B/zh active Active
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109577970A (zh) * | 2018-12-11 | 2019-04-05 | 西南石油大学 | 一种页岩储层ⅰ型断裂韧性的测井评价方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN111898785A (zh) | 2020-11-06 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US11098565B2 (en) | Method for estimating permeability of fractured rock formations from induced slow fluid pressure waves | |
| CN105221133B (zh) | 一种基于测井多参数确定烃源岩有机碳含量的方法和装置 | |
| CN104453874B (zh) | 一种基于核磁共振的砂砾岩储层含油饱和度的计算方法 | |
| CN105277982B (zh) | 一种泥页岩总有机碳含量地震预测方法 | |
| Socco et al. | Feasibility of ultrasonic tomography for nondestructive testing of decay on living trees | |
| Bonal et al. | Wavelet analysis of bender element signals | |
| CN108562938B (zh) | 一种消除频散效应的方法、装置及系统 | |
| Liu et al. | An integrated 3D fracture network reconstruction method based on microseismic events | |
| CN110284873A (zh) | 一种油井储集性质的检测方法及检测装置 | |
| CN103513270A (zh) | 一种基于岩石声学特性的气层识别评价方法及装置 | |
| CN111381292B (zh) | 一种预测砂岩含烃储层的测井解释方法与装置 | |
| CN103091342B (zh) | 一种对岩芯样品进行ct扫描分析处理的方法 | |
| CN106355571A (zh) | 一种白云岩储层质量的确定方法及装置 | |
| CN113158315A (zh) | 一种基于静力触探数据的岩土体参数三维非平稳条件随机场建模方法 | |
| CN107861162B (zh) | 基于微电极测井资料的天然裂缝识别方法及系统 | |
| Lv et al. | The effects of sorbing and non-sorbing gases on ultrasonic wave propagation in fractured coal | |
| CN110344822B (zh) | 含钙质致密储层含水饱和度的确定方法 | |
| CN113216945A (zh) | 一种致密砂岩储层渗透性定量评价方法 | |
| CN109901238A (zh) | 一种基于应力差电阻率实验的高应力地层电阻率校正方法 | |
| CN111898785B (zh) | 基于页岩的断裂韧性空间分布特征预测方法及系统 | |
| Zhao | Investigating brittle rock failure and associated seismicity using laboratory experiments and numerical simulations | |
| CN117094222A (zh) | 基于神经网络的岩土体物性参数推演物理力学指标方法 | |
| Govindaraj | A new approach for estimation of properties of metamorphic rocks | |
| CN108490490B (zh) | 裂缝带表征方法、装置、电子设备及计算机存储介质 | |
| Song et al. | Uniaxial compressive strength estimation based on the primary wave velocity in coal: considering scale effect and anisotropy |
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 |