CN113703047B - 一种预测水平地应力差的方法、装置及存储介质 - Google Patents
一种预测水平地应力差的方法、装置及存储介质 Download PDFInfo
- Publication number
- CN113703047B CN113703047B CN202111097949.4A CN202111097949A CN113703047B CN 113703047 B CN113703047 B CN 113703047B CN 202111097949 A CN202111097949 A CN 202111097949A CN 113703047 B CN113703047 B CN 113703047B
- Authority
- CN
- China
- Prior art keywords
- horizontal
- ground stress
- stress difference
- representing
- horizontal ground
- 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 59
- 238000003860 storage Methods 0.000 title claims abstract description 15
- 239000012530 fluid Substances 0.000 claims abstract description 13
- 238000010276 construction Methods 0.000 claims description 10
- 238000004590 computer program Methods 0.000 claims description 5
- 230000005484 gravity Effects 0.000 claims description 5
- 230000001133 acceleration Effects 0.000 claims description 4
- 238000005211 surface analysis Methods 0.000 claims description 4
- 230000010354 integration Effects 0.000 claims description 2
- 238000011161 development Methods 0.000 abstract description 7
- 102100023340 3-ketodihydrosphingosine reductase Human genes 0.000 abstract description 4
- 101001050680 Homo sapiens 3-ketodihydrosphingosine reductase Proteins 0.000 abstract description 4
- 238000012544 monitoring process Methods 0.000 description 6
- 238000004364 calculation method Methods 0.000 description 4
- 230000015572 biosynthetic process Effects 0.000 description 3
- 238000005553 drilling Methods 0.000 description 3
- 238000011156 evaluation Methods 0.000 description 3
- 230000003287 optical effect Effects 0.000 description 3
- 238000012545 processing Methods 0.000 description 3
- 238000005070 sampling Methods 0.000 description 3
- 230000005483 Hooke's law Effects 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 230000003068 static effect Effects 0.000 description 2
- 238000004422 calculation algorithm Methods 0.000 description 1
- 238000013016 damping Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 238000000691 measurement method Methods 0.000 description 1
- 239000011148 porous material Substances 0.000 description 1
- 239000004575 stone Substances 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/30—Analysis
- G01V1/306—Analysis for determining physical properties of the subsurface, e.g. impedance, porosity or attenuation profiles
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本公开涉及一种预测水平地应力差的方法、装置及存储介质，其中该方法包括：确定水平地应力差预测模型的输入参数的值，其中，水平地应力差预测模型是在组合弹簧地应力模型中引入各向异性引起的水平应力扰动得到的，以各向异性介质的法向柔度和杨氏模量指示水平应力扰动，输入参数包括：垂向地应力、流体压力、有效压力系数、泊松比、杨氏模量、法向柔度、水平最大主应变，以及水平最小主应变；使用水平地应力差预测模型，根据输入参数的值确定水平地应力差。实现了水平地应力差绝对值的预测，相比水平地应力差异系数(DHSR)，实用性和适用性更强，可有效指导页岩气勘探及开发。
Description
技术领域
本公开涉及油气勘探领域，尤其涉及一种预测水平地应力差的方法、装置及存储介质。
背景技术
地应力是油气勘探开发领域的重要研究内容，其重要作用贯穿于各个阶段，勘探阶段与构造及其形成演化密切相关，开发阶段是井筒稳定性的重要考虑因素，同时在生产阶段可以预防套损等。近年来页岩气勘探开发实践揭示，水平地应力差是决定压裂能否形成复杂缝网的关键参数，水平地应力差的评价及预测对于页岩气商业勘探开发至关重要。
相关技术中，主要的地应力研究方法主有测量法、测井计算法、数值模拟方法及地震预测方法。测量法数据较为准确，方法又包括声发射法及水利压裂法等；测井法主要是利用测井的声波、密度等，基于一定的理论进行计算，在实际资料中也有一定的应用，尤其黄荣樽教授(1983)考虑了构造应力对地应力的影响，根据不同的地区选取不同的构造应力值，建立了一种地应力计算模式，称为黄氏模型，后又有学者根据广义胡克定律推导了均匀各向同性线弹性介质的地应力计算模型，称为弹簧模型；但这两种方法只有在钻井的基础上进行，成本较高，尤其无法在钻前预知和评价。
为了能够对地下介质进行有效的钻前预测，数值模拟是一种有效的方法，通过建立区域力学地质模型，不断调整边界值迭代计算，模拟应力场分布规律，这种方法横向分辨率较低，且计算效率较低，无法满足大面积勘探三维区的评价需求。
基于地震资料的地应力预测方法中，常用的水平地应力差异系数(DHSR)预测方法是由Gray(2012)基于地应力的基础理论以及地震各向异性介质地震岩石物理理论，推导的地应力差系数与各种弹性参数、各向异性参数之间的定量关系，该方法基于水平介质横向位移为零的假设条件，只能实现水平地应力差异比(DHSR)的预测，无法得到水平地应力差的大小，影响了地应力的有效评价。
发明内容
为了解决上述技术问题或者至少部分地解决上述技术问题，本公开提供了一种预测水平地应力差的方法、装置及存储介质。
第一方面，本公开提供了一种预测水平地应力差的方法，包括：确定水平地应力差预测模型的输入参数的值，其中，水平地应力差预测模型是在组合弹簧地应力模型中引入各向异性引起的水平应力扰动得到的，以各向异性介质的法向柔度和杨氏模量指示水平应力扰动，输入参数包括：垂向地应力、流体压力、有效压力系数、泊松比、杨氏模量、法向柔度、水平最大主应变，以及水平最小主应变；使用水平地应力差预测模型，根据输入参数的值确定水平地应力差。
在一些实施例中，水平地应力差预测模型表示为：
其中，Δσ表示水平地应力差，σV表示垂向地应力，PP表示流体压力，α表示有效压力系数，ZN表示法向柔度，E表示杨氏模量，v表示泊松比，εH表示水平最大主应变；εh表示水平最小主应变。
在一些实施例中，确定法向柔度的值，包括：接收叠前地震数据，其中，叠前地震数据包括多个方位且每个方位包括多个入射角对应的数据；根据叠前地震数据进行弹性阻抗反演，得到方位弹性阻抗数据，其中，方位弹性阻抗数据包括多个方位且每个方位的多个入射角对应的弹性阻抗；使用方位弹性阻抗方程根据方位弹性阻抗数据确定拉梅系数、剪切模量和法向弱度的值，其中，方位弹性阻抗方程表征方位弹性阻抗与拉梅系数、剪切模量和法向弱度的确定性关系；根据拉梅系数、剪切模量和法向弱度的值确定法向柔度的值。
在一些实施例中，确定泊松比和杨氏模量的值，包括：根据拉梅系数和剪切模量确定泊松比和杨氏模量的值。
在一些实施例中，方位弹性阻抗方程表示为：
其中，
其中，g表示横波速度与纵波速度的比值的平方，λ、μ、ρ分别表示各向异性介质的拉梅系数、剪切模量与密度，ΔN、ΔT分别表示各向异性介质Schoenberg模型的法向弱度和切向弱度，Δ表示为上下介质参数的差异，θ表示地震波入射角，φ为方位角，AEI(θ，φ)表示入射角为θ方位角为φ时的弹性阻抗，EI0、λ0、μ0、ρ0分别表示介质的平均弹性阻抗、平均拉梅系数、平均剪切模量、平均密度。
在一些实施例中，确定水平最大主应变和水平最小主应变的值，包括：根据薄板理论和趋势面分析确定水平最大主应变和水平最小主应变的值，其中，水平最大主应变的值为构造趋势面的垂直构造方向的二阶导数，水平最小主应变的值为构造趋势面的平行构造方向的二阶导数。
在一些实施例中，确定垂向地应力的值，包括：根据密度、重力加速度在深度域积分确定垂向地应力的值。
第二方面，本公开提供了一种预测水平地应力差的装置，包括：第一确定模块，用于确定水平地应力差预测模型的输入参数的值，其中，水平地应力差预测模型是在组合弹簧地应力模型中引入各向异性引起的水平应力扰动得到的，以各向异性介质的法向柔度和杨氏模量指示水平应力扰动，输入参数包括：垂向地应力、流体压力、有效压力系数、泊松比、杨氏模量、法向柔度、水平最大主应变，以及水平最小主应变；第二确定模块，用于使用水平地应力差预测模型，根据输入参数的值确定水平地应力差。
第三方面，本公开提供了一种计算机设备，该计算机设备包括：存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序；计算机程序被处理器执行时实现本公开的任一方法的步骤。
第四方面，本公开提供了一种计算机可读存储介质，该计算机可读存储介质上存储有预测水平地应力差的程序，该预测水平地应力差的程序被处理器执行时实现本公开的任一预测水平地应力差的方法的步骤。
本公开实施例提供的上述技术方案与相关技术相比具有如下优点：本公开实施例提供的该方法，实现了水平地应力差绝对值的预测，相比水平地应力差异系数(DHSR)，实用性和适用性更强，可有效指导页岩气勘探及开发。
附图说明
此处的附图被并入说明书中并构成本说明书的一部分，示出了符合本公开的实施例，并与说明书一起用于解释本公开的原理。
为了更清楚地说明本公开实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，对于本领域普通技术人员而言，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。
图1为本公开实施例提供的预测水平地应力差的方法一种实施方式的流程图；
图2为本公开实施例提供的使用叠前地震数据确定输入参数一种实施方式的流程图；
图3为本公开实施例一个示例的水平地应力差与微地震监测事件的对比图；
图4为本公开实施例提供的预测水平地应力差的装置一种实施方式的结构框图；
图5为本公开实施例提供的计算机设备一种实施方式的硬件示意图。
具体实施方式
应当理解，此处所描述的具体实施例仅仅用以解释本公开，并不用于限定本公开。
在后续的描述中，使用用于表示元件的诸如“模块”、“部件”或“单元”的后缀仅为了有利于本公开的说明，其本身没有特定的意义。因此，“模块”、“部件”或“单元”可以混合地使用。
本公开实施例提供了一种预测水平地应力差的方法，如图1所示，该方法包括步骤S102至步骤S104。
步骤S102，确定水平地应力差预测模型的输入参数的值。
在本公开实施例中，水平地应力差预测模型是在组合弹簧地应力模型中引入各向异性引起的水平应力扰动得到的，以各向异性介质的法向柔度和杨氏模量指示水平应力扰动。在本公开中，裂隙、裂缝等导致各向异性，各向异性引起水平应力扰动。杨氏模量表征各向同性固有刚度属性，法向柔度表征各向异性对固有各向同性属性引起的扰动。在本公开实施例中，以各向异性介质的法向柔度和杨氏模量指示水平应力扰动。
在本公开实施例中，水平地应力差预测模型的输入参数包括：垂向地应力、流体压力(也称为孔隙压力)、有效压力系数(Biot系数)、泊松比、杨氏模量、法向柔度、水平最大主应变，以及水平最小主应变。
步骤S104，使用水平地应力差预测模型，根据输入参数的值确定水平地应力差。
水平地应力差预测模型
在一些实施例中，应用地震各向异性介质理论(包括本构方程、Thomsen各向异性参数、Schoenberg线性滑动理论)，薄板理论及组合弹簧地应力模型，建立水平地应力差预测模型，该水平地应力差预测模型表示为：
其中，Δσ表示水平地应力差，σV表示垂向地应力，PP表示流体压力，α表示有效压力系数，ZN表示法向柔度，E表示杨氏模量，v表示泊松比，εH表示水平最大主应变；εh表示水平最小主应变。
在该水平地应力差预测模型中，第一项可以理解为由于垂直地应力引起的水平地应力分量差异，其可以通过各向异性介质理论下的各向异性特征参数及弹性参数求取；第二项可以理解为薄板理论下由于构造作用引起的水地平应力差，可通过广义胡克定律求取。
在本公开实施例中，可以通过多种公知的方法确定本公开的上述输入参数，如测井评价技术，然而相关技术中的方法可操作性不强。因此，本公开实施例提供了一种使用叠前地震确定部分上述输入参数的方法，以实现使用叠前地震数据预测水平地应力差，可操作性更强。
法向柔度
在一些实施例中，使用叠前地震数据确定法向柔度的值，如图2所示，该方法包括步骤S202至步骤S208。
步骤S202，接收叠前地震数据。
在本公开实施例中，叠前地震数据包括多个方位且每个方位包括多个入射角对应的数据。
在本公开中，对于一定区域的地层，可获得多道地震数据。每道地震数据包括多个方位且每个方位包括多个入射角对应的数据。每道地震数据包括一位置沿地层深度方向上的多个采样点的数据。每个采样点包括多个方位且每个方位包括多个入射对应的数据。
步骤S204，根据叠前地震数据进行弹性阻抗反演，得到方位弹性阻抗数据。
在本公开实施例中，方位弹性阻抗数据包括多个方位且每个方位的多个入射角对应的弹性阻抗。地震数据对应的每个采样点包括不同方位及不同入射角的弹性阻抗。
在本公开实施例中，可采用公知的方法进行弹性阻抗反演，如基于贝叶斯的弹性阻抗反演方法，本公开对此不做赘述。
在本公开实施例中，在进行弹性阻抗反演之前，可对原始叠前地震数据进行预处理，例如开展部分叠加处理，本公开实施例对此不做限定。
步骤S206，使用方位弹性阻抗方程根据方位弹性阻抗数据确定拉梅系数、剪切模量和法向弱度的值。
在本公开实施例中，方位弹性阻抗方程表征方位弹性阻抗与拉梅系数、剪切模量和法向弱度的确定性关系。由此，实现先由叠前地震数据反演得到弹性阻抗，再由弹性阻抗确定拉梅系数、剪切模量和法向弱度的值。由于弹性阻抗反演稳定性高，使得确定输入参数的可靠性高。
在一些实施例中，方位弹性阻抗方程表示为：
其中，
其中，g表示横波速度与纵波速度的比值的平方，λ、μ、ρ分别表示各向异性介质的拉梅系数、剪切模量与密度，ΔN、ΔT分别表示各向异性介质Schoenberg模型的法向弱度和切向弱度，Δ表示为上下介质参数的差异，θ表示地震波入射角，φ为方位角，AEI(θ，φ)表示入射角为θ方位角为φ时的弹性阻抗，EI0、λ0、μ0、ρ0分别表示介质的平均弹性阻抗、平均拉梅系数、平均剪切模量、平均密度。
示例性的，在本公开的一些实施例中，基于上述方位弹性阻抗方程，加入测井约束，通过阻尼最小二乘算法求解弹性阻抗矩阵，开展各向异性参数和弹性参数反演，直接一次求取拉梅系数、剪切模量和法向弱度。
步骤S208，根据拉梅系数、剪切模量和法向弱度的值确定法向柔度的值。
在本公开实施例中，法向柔度的确定方式如下所示：
其中，ZN为法向柔度，λ为各向异性介质的拉梅系数，μ为各向异性介质的剪切模量，ΔN为各向异性介质的法向弱度。
在一些实施例中，如图2所示，还包括步骤S210，根据拉梅系数和剪切模量确定泊松比和杨氏模量的值。具体参见：
式中，λ为拉梅系数，μ为剪切模量，v为泊松比，E为杨氏模量。
通过本公开实施例，实现了地震数据对应的采样点对应的输入参数的确定。
水平最大、最小主应变
在一些实施例中，根据薄板理论和趋势面分析确定水平最大主应变和水平最小主应变的值，其中，水平最大主应变的值为构造趋势面的垂直构造方向的二阶导数，水平最小主应变的值为构造趋势面的平行构造方向的二阶导数。
在本公开实施例中，结合薄板理论，构造趋势面二阶导数可以表征弹性应变，则最大、最小主应变可表示为垂直构造方向和平行构造方向的二阶导数。
其中，z表示为构造趋势面的海拔；x垂直表示垂直构造方向；x平行表示平行构造方向。
垂向地应力
在一些实施例中，确定垂向地应力的值，包括：根据密度、重力加速度在深度域积分确定垂向地应力的值，本公开实施例对此不做赘述。
流体压力
在本公开实施例中，流体压力可通过公知的压力预测方法确定，例如，Eaton法、Fillippone法、Stone法等，本公开实施例对比不做赘述和限定。
在本公开实施例的一个优选示例中，优化组合弹簧地应力模型，提出了一种水平地应力差预测模型(如式(1所示))；利用各向异性介质方位AVO的反射系数近似方程，通过叠前地震数据反演各向异性参数和弹性参数(如式(2)所示)，并求取法向柔度(如式(3)所示)；基于薄板理论，用构造趋势面二阶导数表征应变，求取水平最大和最小主应变(如式(4)所示)；使用求出的参数利用水平地应力差预测模型实现水平地应力差预测。
其中，反演理论公式为叠前各向异性介质方位AVO的反射系数近似方程：
公式(6)中，g为横波速度与纵波速度的比值的平方，λ、μ、ρ分别表示各向异性介质的拉梅系数、剪切模量与密度，ΔN、ΔT分别表示各向异性介质Schoenberg模型的法向弱度和切向弱度，Δ表示为上下介质参数的差异，θ表示地震波入射角，φ为方位角。
由于弹性阻抗反演稳定性高，是实际生产中得到了广泛的应用，根据Connolly(1999)推导弹性阻抗以及Whitcombe(2002)标准化弹性阻抗的思想，因而可以将公式(6)推导为方位弹性阻抗方程(如式(2)所示)。
该示例，同时考虑了重力引起的水平应力分量差异和由于构造作用引起的水平地应力差异，且实现了水平地应力差绝对值的预测，相比水平地应力差异系数(DHSR)，实用性和适用性更强，可有效指导页岩气勘探及开发。
图3示出了基于该示例预测的水平地应力差与微地震监测事件的对比图，图3中，色标的灰度程度指示水平地应力差的大小，圆形气泡代表微地震监测事件。从图3中可以看出，A井东侧水平地应力差较低，易于压裂，东侧微地震监测事件数多、微地震能量强；其中5-7、13-14、18-22段压裂微地震监测事件数较高，水平地应力差预测低值井段与微地震监测事件基本一致，吻合率大于80％。
本公开实施例还提供了一种预测水平地应力差的装置，下面结合图4对本公开实施例的装置进行描述。
如图4所示，该装置包括：第一确定模块410，用于确定水平地应力差预测模型的输入参数的值。第二确定模块420，与第一确定模块410相连，用于使用水平地应力差预测模型，根据输入参数的值确定水平地应力差。
在本公开中，水平地应力差预测模型是在组合弹簧地应力模型中引入各向异性引起的应力扰动得到的，以各向异性介质的法向柔度和杨氏模量指示应力扰动，输入参数包括：垂向地应力、流体压力、有效压力系数、泊松比、杨氏模量、法向柔度、水平最大主应变，以及水平最小主应变。本公开的一个是示例性水平地应力差预测模型如本公开前述式(1)所示，在此不做赘述。
在一些实施例中，第一确定模块410确定法向柔度的值，具体包括：接收叠前地震数据，其中，叠前地震数据包括多个方位且每个方位包括多个入射角对应的数据；根据叠前地震数据进行弹性阻抗反演，得到方位弹性阻抗数据，其中，方位弹性阻抗数据包括多个方位且每个方位的多个入射角对应的弹性阻抗；使用方位弹性阻抗方程(示例性的，如本公开前述式(2)所示)根据方位弹性阻抗数据确定拉梅系数、剪切模量和法向弱度的值，其中，方位弹性阻抗方程表征方位弹性阻抗与拉梅系数、剪切模量和法向弱度的确定性关系；根据拉梅系数、剪切模量和法向弱度的值确定法向柔度的值。
在一些实施例中，第一确定模块410确定泊松比和杨氏模量的值，具体包括：根据拉梅系数和剪切模量确定泊松比和杨氏模量的值。
在一些实施例中，第一确定模块410确定水平最大主应变和水平最小主应变的值，具体包括：根据薄板理论和趋势面分析确定水平最大主应变和水平最小主应变的值，其中，水平最大主应变的值为构造趋势面的垂直构造方向的二阶导数(如本公开前述式(5)所示)，水平最小主应变的值为构造趋势面的平行构造方向的二阶导数(如本公开前述式(5)所示)。
在一些实施例中，第一确定模块410确定垂向地应力的值，包括：根据密度、重力加速度在深度域积分确定垂向地应力的值。
本公开实施例还提供了一种计算机设备。图5为本公开实施例提供的计算机设备一种实施方式的硬件结构示意图，如图5所示，本公开实施例的计算机设备10包括：至少包括但不限于：可通过系统总线相互通信连接的存储器11和处理器12。需要指出的是，图5仅示出了具有组件11-12的计算机设备10，但是应理解的是，并不要求实施所有示出的组件，可以替代的实施更多或者更少的组件。
本实施例中，存储器11(即可读存储介质)包括闪存、硬盘、多媒体卡、卡型存储器(例如，SD或DX存储器等)、随机访问存储器(RAM)、静态随机访问存储器(SRAM)、只读存储器(ROM)、电可擦除可编程只读存储器(EEPROM)、可编程只读存储器(PROM)、磁性存储器、磁盘、光盘等。在一些实施例中，存储器11可以是计算机设备10的内部存储单元，例如计算机设备10的硬盘或内存。在另一些实施例中，存储器11也可以是计算机设备10的外部存储设备，例如该计算机设备10上配备的插接式硬盘，智能存储卡(Smart Media Card,SMC)，安全数字(Secure Digital,SD)卡，闪存卡(Flash Card)等。当然，存储器11还可以既包括计算机设备10的内部存储单元也包括其外部存储设备。本实施例中，存储器11通常用于存储安装于计算机设备10的操作系统和各类软件。此外，存储器11还可以用于暂时地存储已经输出或者将要输出的各类数据。
处理器12在一些实施例中可以是中央处理器(Central Processing Unit，CPU)、控制器、微控制器、微处理器、或其他数据处理芯片。该处理器12通常用于控制计算机设备10的总体操作。本实施例中，处理器12用于运行存储器11中存储的程序代码或者处理数据，例如本公开实施例的任一或多个方法。
本实施例还提供一种计算机可读存储介质，如闪存、硬盘、多媒体卡、卡型存储器(例如，SD或DX存储器等)、随机访问存储器(RAM)、静态随机访问存储器(SRAM)、只读存储器(ROM)、电可擦除可编程只读存储器(EEPROM)、可编程只读存储器(PROM)、磁性存储器、磁盘、光盘、服务器、App应用商城等等，其上存储有计算机程序，程序被处理器执行时实现相应功能。本实施例的计算机可读存储介质用于存储本公开实施例的任一或多个的程序代码，被处理器执行时实现本公开实施例的任一或多个的方法。
需要说明的是，在本文中，术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、物品或者装置不仅包括那些要素，而且还包括没有明确列出的其他要素，或者是还包括为这种过程、方法、物品或者装置所固有的要素。在没有更多限制的情况下，由语句“包括一个……”限定的要素，并不排除在包括该要素的过程、方法、物品或者装置中还存在另外的相同要素。
上述本公开实施例序号仅仅为了描述，不代表实施例的优劣。
通过以上的实施方式的描述，本领域的技术人员可以清楚地了解到上述实施例方法可借助软件加必需的通用硬件平台的方式来实现，当然也可以通过硬件，但很多情况下前者是更佳的实施方式。基于这样的理解，本公开的技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来，该计算机软件产品存储在一个存储介质(如ROM/RAM、磁碟、光盘)中，包括若干指令用以使得一台终端(可以是手机，计算机，服务器，空调器，或者网络设备等)执行本公开各个实施例所述的方法。
上面结合附图对本公开的实施例进行了描述，但是本公开并不局限于上述的具体实施方式，上述的具体实施方式仅仅是示意性的，而不是限制性的，本领域的普通技术人员在本公开的启示下，在不脱离本公开宗旨和权利要求所保护的范围情况下，还可做出很多形式，这些均属于本公开的保护之内。
Claims (9)
1.一种预测水平地应力差的方法，其特征在于，包括：
确定水平地应力差预测模型的输入参数的值，其中，所述水平地应力差预测模型是在组合弹簧地应力模型中引入各向异性引起的水平应力扰动得到的，以各向异性介质的法向柔度和杨氏模量指示所述水平应力扰动，所述输入参数包括：垂向地应力、流体压力、有效压力系数、泊松比、杨氏模量、法向柔度、水平最大主应变，以及水平最小主应变；
使用所述水平地应力差预测模型，根据所述输入参数的值确定水平地应力差；
所述水平地应力差预测模型表示为：
其中，Δσ表示所述水平地应力差，σV表示所述垂向地应力，PP表示所述流体压力，α表示所述有效压力系数，ZN表示所述法向柔度，E表示所述杨氏模量，v表示所述泊松比，εH表示所述水平最大主应变；εh表示所述水平最小主应变。
2.根据权利要求1所述的方法，其特征在于，确定所述法向柔度的值，包括：
接收叠前地震数据，其中，所述叠前地震数据包括多个方位且每个方位包括多个入射角对应的数据；
根据所述叠前地震数据进行弹性阻抗反演，得到方位弹性阻抗数据，其中，所述方位弹性阻抗数据包括多个方位且每个方位的多个入射角对应的弹性阻抗；
使用方位弹性阻抗方程根据所述方位弹性阻抗数据确定拉梅系数、剪切模量和法向弱度的值，其中，所述方位弹性阻抗方程表征方位弹性阻抗与拉梅系数、剪切模量和法向弱度的确定性关系；
根据所述拉梅系数、所述剪切模量和所述法向弱度的值确定法向柔度的值。
3.根据权利要求2所述的方法，其特征在于，确定所述泊松比和所述杨氏模量的值，包括：根据所述拉梅系数和所述剪切模量确定所述泊松比和所述杨氏模量的值。
4.根据权利要求2所述的方法，其特征在于，所述方位弹性阻抗方程表示为：
其中，
其中，g表示横波速度与纵波速度的比值的平方，λ、μ、ρ分别表示各向异性介质的拉梅系数、剪切模量与密度，ΔN、ΔT分别表示各向异性介质Schoenberg模型的法向弱度和切向弱度，Δ表示为上下介质参数的差异，θ表示地震波入射角，φ为方位角，AEI(θ，φ)表示入射角为方位角为φ时的弹性阻抗，EI0、λ0、μ0、ρ0分别表示介质的平均弹性阻抗、平均拉梅系数、平均剪切模量、平均密度。
5.根据权利要求1所述的方法，其特征在于，确定所述水平最大主应变和所述水平最小主应变的值，包括：
根据薄板理论和趋势面分析确定所述水平最大主应变和所述水平最小主应变的值，其中，
所述水平最大主应变的值为构造趋势面的垂直构造方向的二阶导数；
所述水平最小主应变的值为所述构造趋势面的平行构造方向的二阶导数。
6.根据权利要求1所述的方法，其特征在于，确定所述垂向地应力的值，包括：根据密度、重力加速度在深度域积分确定所述垂向地应力的值。
7.一种预测水平地应力差的装置，其特征在于，包括：
第一确定模块，用于确定水平地应力差预测模型的输入参数的值，其中，所述水平地应力差预测模型是在组合弹簧地应力模型中引入各向异性引起的水平应力扰动得到的，以各向异性介质的法向柔度和杨氏模量指示所述水平应力扰动，所述输入参数包括：垂向地应力、流体压力、有效压力系数、泊松比、杨氏模量、法向柔度、水平最大主应变，以及水平最小主应变；
第二确定模块，用于使用所述水平地应力差预测模型，根据所述输入参数的值确定水平地应力差；
所述水平地应力差预测模型表示为：
其中，Δσ表示所述水平地应力差，σV表示所述垂向地应力，PP表示所述流体压力，α表示所述有效压力系数，ZN表示所述法向柔度，E表示所述杨氏模量，v表示所述泊松比，εH表示所述水平最大主应变；εh表示所述水平最小主应变。
8.一种计算机设备，其特征在于，所述计算机设备包括：
存储器、处理器及存储在所述存储器上并可在所述处理器上运行的计算机程序；
所述计算机程序被所述处理器执行时实现如权利要求1至6中任一项所述的方法的步骤。
9.一种计算机可读存储介质，其特征在于，所述计算机可读存储介质上存储有预测水平地应力差的程序，所述预测水平地应力差的程序被处理器执行时实现如权利要求1至6中任一项所述的预测水平地应力差的方法的步骤。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111097949.4A CN113703047B (zh) | 2021-09-18 | 2021-09-18 | 一种预测水平地应力差的方法、装置及存储介质 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111097949.4A CN113703047B (zh) | 2021-09-18 | 2021-09-18 | 一种预测水平地应力差的方法、装置及存储介质 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113703047A CN113703047A (zh) | 2021-11-26 |
| CN113703047B true CN113703047B (zh) | 2024-04-09 |
Family
ID=78661189
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202111097949.4A Active CN113703047B (zh) | 2021-09-18 | 2021-09-18 | 一种预测水平地应力差的方法、装置及存储介质 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113703047B (zh) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114371076B (zh) * | 2022-01-06 | 2024-08-20 | 上海电气集团股份有限公司 | 工件应力值的测试方法、系统、电子设备及存储介质 |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2011010989A1 (en) * | 2009-07-20 | 2011-01-27 | Services Petroliers Schlumberger | Estimating formation stresses using radial profiles of three shear moduli |
| CN109856674A (zh) * | 2019-03-04 | 2019-06-07 | 西南石油大学 | 工程甜点评测方法及装置 |
| CN110456414A (zh) * | 2019-08-21 | 2019-11-15 | 河南理工大学 | 矿区应力方向估算方法 |
| CN111460602A (zh) * | 2019-12-06 | 2020-07-28 | 西南石油大学 | 一种基于岩石物理建模的横观各向同性地层地应力预测方法 |
| CN113283108A (zh) * | 2021-06-10 | 2021-08-20 | 西安石油大学 | 一种定量评价页岩油储层可压裂性的方法及系统 |
| CN113341458A (zh) * | 2021-04-25 | 2021-09-03 | 中国石油大学(华东) | 基于构造应变的水平主应力确定方法、装置及存储介质 |
| CN113341465A (zh) * | 2021-06-11 | 2021-09-03 | 中国石油大学(北京) | 方位各向异性介质的地应力预测方法、装置、介质及设备 |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP2529254A4 (en) * | 2010-01-25 | 2017-08-02 | CGG Veritas Services (U.S.) Inc. | Methods and systems for estimating stress using seismic data |
-
2021
- 2021-09-18 CN CN202111097949.4A patent/CN113703047B/zh active Active
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2011010989A1 (en) * | 2009-07-20 | 2011-01-27 | Services Petroliers Schlumberger | Estimating formation stresses using radial profiles of three shear moduli |
| CN109856674A (zh) * | 2019-03-04 | 2019-06-07 | 西南石油大学 | 工程甜点评测方法及装置 |
| CN110456414A (zh) * | 2019-08-21 | 2019-11-15 | 河南理工大学 | 矿区应力方向估算方法 |
| CN111460602A (zh) * | 2019-12-06 | 2020-07-28 | 西南石油大学 | 一种基于岩石物理建模的横观各向同性地层地应力预测方法 |
| CN113341458A (zh) * | 2021-04-25 | 2021-09-03 | 中国石油大学(华东) | 基于构造应变的水平主应力确定方法、装置及存储介质 |
| CN113283108A (zh) * | 2021-06-10 | 2021-08-20 | 西安石油大学 | 一种定量评价页岩油储层可压裂性的方法及系统 |
| CN113341465A (zh) * | 2021-06-11 | 2021-09-03 | 中国石油大学(北京) | 方位各向异性介质的地应力预测方法、装置、介质及设备 |
Non-Patent Citations (3)
| Title |
|---|
| Estimating the transversely isotropic elastic intact rock properties for in situ stress measurement data reduction: a case study of the Olkiluoto mica gneiss, Finland;Hakala, Matti 等;International Journal of Rock Mechanics and Mining Sciences;第44卷(第1期);14-46 * |
| 叠前反演在龙虎泡致密油"工程甜点"预测中的应用;王勇 等;石油物探(第01期);118-125 * |
| 地应力预测技术的研究现状与进展;印兴耀 等;石油物探(第04期);5-21 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113703047A (zh) | 2021-11-26 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US20210277778A1 (en) | Method for evaluating brittleness of deep shale reservoir and computer readable storage medium | |
| US20210255359A1 (en) | Method for estimating rock brittleness from well-log data | |
| CN103336305B (zh) | 一种基于灰色理论划分致密砂岩储层岩石物理相的方法 | |
| CN104749622A (zh) | 一种基于岩石物理的泥页岩可压性定量表征方法 | |
| Fan et al. | Prediction of the horizontal stress of the tight sandstone formation in eastern Sulige of China | |
| CN113703047B (zh) | 一种预测水平地应力差的方法、装置及存储介质 | |
| CN104153768A (zh) | 一种评价花岗岩储层储集性能的方法 | |
| Cheng et al. | Identifying and estimating formation stress from borehole monopole and cross-dipole acoustic measurements | |
| CN110727031B (zh) | 一种基于三维叠前地震反演结果的地应力获取方法 | |
| CN104483706B (zh) | 一种基于煤岩力学参数的煤体结构测井定量识别方法 | |
| Falcon‐Suarez et al. | CO2‐brine substitution effects on ultrasonic wave propagation through sandstone with oblique fractures | |
| Hua et al. | Rock mass deformation modulus estimation models based on in situ tests | |
| CN113341458B (zh) | 基于构造应变的水平主应力确定方法、装置及存储介质 | |
| CN114063162B (zh) | 一种基于小样本机器学习SVM Vp/Vs预测方法 | |
| Zhang et al. | Natural fracture prediction in Keshen 2 ultra-deep tight gas reservoir based on R/S analysis, Kuqa Depression, Tarim Basin | |
| CN109190292A (zh) | 基于测井资料的含水层富水性预测方法 | |
| US20210381354A1 (en) | Geomechanical applications using biots coefficient based on advanced downhole logging measurements | |
| CN113514890B (zh) | 一种使用地震数据预测地层压力系数的方法、装置及设备 | |
| CN113777668A (zh) | 用于砂泥岩互层致密气藏的地应力计算方法及装置 | |
| Joshi et al. | Applicability of attenuation relations for regional studies | |
| Abdideh et al. | Analysis of deep stress field using well log and wellbore breakout data: a case study in Cretaceous oil reservoir, southwest Iran | |
| CN114441402A (zh) | 一种评价致密砂岩渗透性的方法 | |
| Vasco et al. | Seismicity and Stress Associated With a Fluid‐Driven Fracture: Estimating the Evolving Geometry | |
| Huang et al. | Development characteristics and prediction of fractures in deep-ultradeep tight sandstone reservoirs of Kuqa Depression: a case study of the Bashijiqike and Baxigai formations in Dabei 9 Block | |
| CN117077446B (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 | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |