CN109683197B - 水平裂缝预测方法、装置和设备 - Google Patents
水平裂缝预测方法、装置和设备 Download PDFInfo
- Publication number
- CN109683197B CN109683197B CN201811502830.9A CN201811502830A CN109683197B CN 109683197 B CN109683197 B CN 109683197B CN 201811502830 A CN201811502830 A CN 201811502830A CN 109683197 B CN109683197 B CN 109683197B
- Authority
- CN
- China
- Prior art keywords
- shear wave
- incident angle
- angle direction
- velocity
- fracture
- 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 35
- 238000011161 development Methods 0.000 claims abstract description 34
- 238000004364 calculation method Methods 0.000 claims description 17
- 102100031418 EF-hand domain-containing protein D2 Human genes 0.000 claims description 12
- 101000802344 Homo sapiens Zinc finger SWIM domain-containing protein 7 Proteins 0.000 claims description 12
- 238000012545 processing Methods 0.000 claims description 8
- 239000003245 coal Substances 0.000 abstract description 9
- 206010017076 Fracture Diseases 0.000 description 102
- 208000010392 Bone Fractures Diseases 0.000 description 90
- 238000004458 analytical method Methods 0.000 description 5
- 238000010586 diagram Methods 0.000 description 5
- 230000000694 effects Effects 0.000 description 5
- 101100188686 Danio rerio opn1sw2 gene Proteins 0.000 description 4
- 102100031414 EF-hand domain-containing protein D1 Human genes 0.000 description 4
- 101150096151 EFHD1 gene Proteins 0.000 description 4
- 101150089053 SWS2 gene Proteins 0.000 description 4
- 230000006870 function Effects 0.000 description 4
- 238000007781 pre-processing Methods 0.000 description 4
- 238000004891 communication Methods 0.000 description 3
- 238000004590 computer program Methods 0.000 description 3
- 238000012937 correction Methods 0.000 description 3
- 239000011435 rock Substances 0.000 description 3
- 230000015572 biosynthetic process Effects 0.000 description 2
- 238000001514 detection method Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 239000004215 Carbon black (E152) Substances 0.000 description 1
- 235000015076 Shorea robusta Nutrition 0.000 description 1
- 244000166071 Shorea robusta Species 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 239000012530 fluid Substances 0.000 description 1
- 229930195733 hydrocarbon Natural products 0.000 description 1
- 150000002430 hydrocarbons Chemical class 0.000 description 1
- 229910052500 inorganic mineral Inorganic materials 0.000 description 1
- 239000011707 mineral Substances 0.000 description 1
Images
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/303—Analysis for determining velocity profiles or travel times
-
- 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
-
- 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/284—Application of the shear wave component and/or several components of the seismic signal
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/62—Physical property of subsurface
- G01V2210/622—Velocity, density or impedance
- G01V2210/6222—Velocity; travel time
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/64—Geostructures, e.g. in 3D data cubes
- G01V2210/646—Fractures
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
本发明提供了一种水平裂缝预测方法、装置和设备，其中，该方法包括：获取目标裂缝处第一入射角方向以及第二入射角方向的纵波速度和横波速度，其中，所述第一入射角小于所述第二入射角；计算所述第一入射方向的第一纵横波速度比和第一横波分裂系数；计算所述第二入射方向的第二纵横波速度比；在所述第一横波分裂系数小于第一预设值，在所述第一纵横波速度比与所述第二纵横波速度比之间的比值小于第二预设值的情况下，确定所述目标裂缝为水平裂缝。在本发明实施例中，通过计算第一入射角方向的纵横波速度比以及第一横波分裂系数、第二入射角方向的纵横波速度比，确定目标裂缝是否为水平裂缝，为煤层气、页岩气的勘探开发提供有效的数据。
Description
技术领域
本发明涉及地球物理勘探与开发技术领域，特别涉及一种水平裂缝预测方法、装置和设备。
背景技术
裂缝广泛分布于地下岩石中，对地震波传播有非常重要的影响，对地下流体储存运移也有着决定性的作用。因此，地层裂缝特征描述是裂缝性油气藏的一项重要研究内容。岩石在受到构造运动的作用会产生高角度裂缝，也可称作是垂直裂缝，这是在常规的油气勘探中主要研究的对象。此外，在页岩和煤层中，由于矿物的定向排列、地应力作用等因素的影响，发育着水平分布的天然裂缝，也可称作是水平裂缝。其中，水平裂缝为页岩气、煤层气提供了储集空间，对储层的含气量起到非常重要的作用，因此，水平裂缝预测对页岩气、煤层气等勘探开发具有重要的意义。
现有的裂缝预测方法如：基于横波分裂的地层裂缝预测方法对裂缝进行预测时只能对垂直裂缝进行有效的预测，而无法识别水平裂缝。
针对如何预测水平裂缝，目前尚未提出有效的解决方案。
发明内容
本发明实施例提供了一种水平裂缝预测方法，以解决现有技术中无法有效预测水平裂缝的问题。
本发明实施例提供了一种水平裂缝预测方法，包括：获取目标裂缝处第一入射角方向的纵波速度和横波速度、第二入射角方向的纵波速度和横波速度，其中，所述第一入射角小于所述第二入射角；根据所述第一入射角方向的纵波速度和横波速度计算所述第一入射方向的第一纵横波速度比和第一横波分裂系数；根据所述第二入射角方向的纵波速度和横波速度计算所述第二入射方向的第二纵横波速度比；在所述第一横波分裂系数小于第一预设值，在所述第一纵横波速度比与所述第二纵横波速度比之间的比值小于第二预设值的情况下，确定所述目标裂缝为水平裂缝。
在一个实施例中，根据所述第一入射角方向的纵波速度和横波速度，计算所述第一入射方向的第一横波分裂系数，包括：根据所述第一入射角方向的横波速度，获取所述第一入射角方向的快横波速度和慢横波速度；根据所述第一入射角方向的快横波速度、慢横波速度，计算所述第一横波分裂系数。
在一个实施例中，按照以下公式，根据所述第一入射角方向的快横波速度、慢横波速度，计算所述第一横波分裂系数：
SWS1＝(Vs11-Vs21)/Vs21
其中，Vs11为所述第一入射角方向的快横波速度；Vs21为所述第一入射角方向的慢横波速度；SWS1为所述第一横波分裂系数。
在一个实施例中，在确定所述目标裂缝为水平裂缝之后，还包括：获取所述第二入射角方向的快横波速度和慢横波速度；根据所述第二入射角方向的快横波速度、慢横波速度计算所述第二入射方向的第二横波分裂系数；将所述第二横波分裂系数作为所述水平裂缝的裂缝密度；根据所述裂缝密度，确定所述目标裂缝的发育程度。
在一个实施例中，所述第一入射角大于等于0°小于等于20°；所述第二入射角大于等于70°小于等于90°。
本发明实施例还提供一种水平裂缝预测装置，包括：获取模块，用于获取目标裂缝处第一入射角方向的纵波速度和横波速度和第二入射角方向的纵波速度和横波速度，其中，所述第一入射角小于所述第二入射角；第一计算模块，用于根据所述第一入射角方向的纵波速度和横波速度计算所述第一入射方向的第一纵横波速度比和第一横波分裂系数；第二计算模块，用于根据所述第二入射角方向的纵波速度和横波速度计算所述第二入射方向的第二纵横波速度比；处理模块，用于在所述第一横波分裂系数小于第一预设值，在所述第一纵横波速度比与所述第二纵横波速度比之间的比值小于第二预设值的情况下，确定所述目标裂缝为水平裂缝。
在一个实施例中，所述第一计算模块包括：第一获取单元，用于根据所述第一入射角方向的横波速度，获取所述第一入射角方向的快横波速度和慢横波速度；第一计算单元，用于根据所述第一入射角方向的快横波速度、慢横波速度，计算所述第一横波分裂系数。
在一个实施例中，所述第一计算单元按照以下公式计算所述第一横波分裂系数：
SWS1＝(Vs11-Vs21)/Vs21
其中，Vs11为所述第一入射角方向的快横波速度；Vs21为所述第一入射角方向的慢横波速度；SWS1为所述第一横波分裂系数。
在一个实施例中，还包括：第二获取单元，用于在确定所述目标裂缝为水平裂缝之后，获取所述第二入射角方向的快横波速度和慢横波速度；第二计算单元，用于根据所述第二入射角方向的快横波速度、慢横波速度计算所述第二入射方向的第二横波分裂系数；处理单元，用于将所述第二横波分裂系数作为所述水平裂缝的裂缝密度；确定单元，用于根据所述裂缝密度，确定所述目标裂缝的发育程度。
本发明实施例还提供一种水平裂缝预测设备，包括处理器以及用于存储处理器可执行指令的存储器，所述处理器执行所述指令时实现所述水平裂缝预测方法的步骤。
本发明实施例还提供一种计算机可读存储介质，其上存储有计算机指令，所述指令被执行时实现所述水平裂缝预测方法的步骤。
在本发明实施例中，提供了一种水平裂缝预测方法，通过计算第一入射角方向的第一纵横波速度比和第一横波分裂系数，以及第二入射角方向的第二纵横波速度比，并根据计算得到的第一横波分裂系数、第一纵横波速度比以及第二纵横波速度比确定目标裂缝是否为水平裂缝，为煤层气、页岩气的勘探开发提供有效的数据，进而有效提高了煤层气、页岩气的产量。
附图说明
此处所说明的附图用来提供对本发明的进一步理解，构成本申请的一部分，并不构成对本发明的限定。在附图中：
图1是根据本发明实施例提供的水平裂缝预测方法步骤示意图；
图2是根据本发明具体实施例提供的水平裂缝预测方法示意图；
图3是根据本发明实施例提供的水平裂缝预测装置的示意图；
图4是根据本发明实施例提供的水平裂缝预测电子设备的示意图。
具体实施方式
下面将参考若干示例性实施方式来描述本发明的原理和精神。应当理解，给出这些实施方式仅仅是为了使本领域技术人员能够更好地理解进而实现本发明，而并非以任何方式限制本发明的范围。相反，提供这些实施方式是为了使本申请公开更加透彻和完整，并且能够将本公开的范围完整地传达给本领域的技术人员。
本领域的技术人员知道，本发明的实施方式可以实现为一种系统、装置设备、方法或计算机程序产品。因此，本申请公开可以具体实现为以下形式，即：完全的硬件、完全的软件(包括固件、驻留软件、微代码等)，或者硬件和软件结合的形式。
考虑到现有的裂缝预测方法如：基于横波分裂的地层裂缝预测方法中激发的地震波其传播路径与水平裂缝平面接近垂直，表现不出水平裂缝处的各向异性特征，因此在使用该方法对裂缝进行预测时只能对垂直裂缝进行有效的预测，而无法有效预测水平裂缝。
基于上述无法有效预测水平裂缝的问题，本发明实施例提供了一种水平裂缝预测方法，如图1所示，可以包括以下步骤：
步骤S101：获取目标裂缝处第一入射角方向的纵波速度和横波速度、第二入射角方向的纵波速度和横波速度，其中，第一入射角小于第二入射角。
多波地震勘探是通过综合利用纵波、横波(或转换波)等多种地震波对地层进行勘探的方法，因此到多波地震资料中携带了比单波(纵波)分量更丰富的地层岩性信息，尤其还携带了地层裂缝信息。因此可以通过对获取的目标裂缝处的多波地震数据进行预处理和速度分析，得到第一入射角方向的纵波速度和横波速度、第二入射角方向的纵波速度和横波速度。其中，预处理包括但不限于：去噪、振幅补偿、增益；速度分析包括但不限于：通过速度扫描获得地震波动校正速度。
上述第一入射角小于第二入射角，具体的，上述第一入射角大于等于0°小于等于20°，第二入射角大于等于70°小于等于90°。其中，上述入射角可以为地震波传播方向与裂缝平面法线方向之间的夹角。进一步的，在地震勘探中，众多炮点(即激发地震波的震源所处位置)和接收点(即接收地震波的检波器所处位置)组成了观测系统。因此，存在不同角度的地震波入射-反射路径，上述地震波传播方向可以通过观测系统中的炮-检空间位置确定。
上述入射角的范围为0°(即垂直于裂缝面)到90°(即平行于裂缝面)，并且在多波地震数据中，各个不同角度的地震波入射-反射路径之间可以选取适当的间隔，保证在0°至20°之间以及70°至90°之间，有且仅有一个方向的多波地震数据，即上述第一入射角以及第二入射角是唯一的。
步骤S102：根据第一入射角方向的纵波速度和横波速度计算第一入射方向的第一纵横波速度比和第一横波分裂系数。
考虑到当地震横波在含有裂缝的各向异性介质中传播时，会产生横波分裂现象，形成快横波和慢横波，而且快、慢横波的差异程度与地层裂缝的发育程度成比例，横波分裂现象主要体现了裂缝信息。因此，可以获取目标裂缝处多波地震数据中第一入射角方向上的快横波速度和慢横波速度，并根据第一入射角方向上的快横波速度、慢横波速度，计算第一横波分裂系数。其中上述第一横波分裂系数通过以下公式得到：
SWS1＝(Vs11-Vs21)/Vs21
其中，Vs11为第一入射角方向的快横波速度；Vs21为第一入射角方向的慢横波速度；SWS1为第一横波分裂系数。
进一步的，上述第一纵横波速度比为目标裂缝处第一入射角方向的纵波速度与横波速度相除的值。
步骤S103：根据第二入射角方向的纵波速度和横波速度计算第二入射方向的第二纵横波速度比。
上述第一纵横波速度比为目标裂缝处第一入射角方向的纵波速度与横波速度相除的值。
步骤S104：在第一横波分裂系数小于第一预设值，在第一纵横波速度比与第二纵横波速度比之间的比值小于第二预设值的情况下，确定目标裂缝为水平裂缝。
其中，上述第一预设值可以为0.02，上述第二预设值可以为0.8。当第一横波分裂系数小于0.02，并且第一纵横波速度比与第二纵横波速度比之间的比值小于0.8时，可以确定目标裂缝为水平裂缝。可以理解的是，上述第一预设值以及上述第二预设值也可以根据实际情况设定为其它合理值，本申请对此不作限定。
裂缝密度是表示裂缝发育程度的一个概念值，用于反映裂缝发育程度。根据裂缝密度可以将目标裂缝划分为裂缝发育带、次发育带或者不发育带，在一个实施例中，可以将第二横波分裂系数作为目标裂缝的裂缝密度，以确定目标裂缝的发育程度。
具体的，在确定目标裂缝为水平裂缝后，可以获取目标裂缝处多波地震数据中第二入射角方向上的快横波速度和慢横波速度，并根据第二入射角方向上的快横波速度、慢横波速度，计算第二横波分裂系数。其中上述第二横波分裂系数通过以下公式得到：
SWS2＝(Vs12-Vs22)/Vs22
其中，Vs12为第二入射角方向的快横波速度；Vs22为第二入射角方向的慢横波速度；SWS2为第二横波分裂系数。
下面结合一个具体实施例对上述方法进行说明，然而，值得注意的是，该具体实施例仅是为了更好地说明本申请，并不构成对本申请的不当限定。
本发明实施提供了一种水平裂缝预测方法，如图2所示，可以包括：
步骤201，采集目标裂缝处的多波地震数据，对多波地震数据进行预处理、速度分析，得到不同入射角方向的纵波速度Vp、横波速度Vs、快横波速度Vs1以及慢横波速度Vs2。
其中，上述预处理包括去噪、振幅补偿、增益等；上述速度分析具体是通过对多波地震数据进行速度扫描获得地震波动校正速度，分析可以得到：不同入射角方向的纵波速度Vp、横波速度Vs、快横波速度Vs1以及慢横波速度Vs2。
步骤202，计算不同方向上纵横波速度比和横波分裂系数。
其中，上述纵横波速度比为纵波速度Vp与横波速度Vs相除的值。
考虑到当地震横波在含有裂缝的各向异性介质中传播时，会产生横波分裂现象，形成快横波和慢横波，而且快、慢横波的差异程度与地层裂缝的发育程度成比例，横波分裂现象主要体现了裂缝信息。因此，可以根据快横波速度Vs1、慢横波速度Vs2，计算横波分裂系数(SWS，shear-wave splitting)。其中上述横波分裂系数SWS通过以下公式得到：
SWS＝(Vs1-Vs2)/Vs2
步骤203，比较较小入射角方向与较大入射角方向纵横波速度比和横波分裂系数，并根据比较结果进行水平裂缝的预测。
其中，上述入射角为地震波传播方向与裂缝平面法线方向之间的夹角，其范围为0°(即垂直于裂缝面)到90°(即平行于裂缝面)。在地震勘探中，众多炮点(即激发地震波的震源所处位置)和接收点(即接收地震波的检波器所处位置)组成了观测系统。因此存在不同角度的地震波入射-反射路径，入射角可以通过观测系统中的炮-检空间位置确定。其中，较小入射角方向为0°～20°范围内的入射角，较大入射角方向为70°～90°范围内的入射角。
水平裂缝预测的判断依据为：在较小入射角横波分裂系小于0.02，较小入射角方向的纵横波速度比与较大入射角方向的纵横波速度比之间的比值小于0.8的情况下，确定目标裂缝为水平裂缝。
进一步的，裂缝密度是表示裂缝发育程度的一个概念值，用于反映裂缝发育程度。根据裂缝密度可以将目标裂缝划分为裂缝发育带、次发育带或者不发育带，在一个具体实施例中，水平裂缝发育程度的评估方式为：将较大入射角方向上的横波分裂系数作为目标裂缝的裂缝密度。
基于同一发明构思，本发明实施例中还提供了一种水平裂缝预测装置，如下面的实施例所述。由于水平裂缝预测装置解决问题的原理与水平裂缝预测方法相似，因此水平裂缝预测装置的实施可以参见水平裂缝预测方法的实施，重复之处不再赘述。以下所使用的，术语“单元”或者“模块”可以实现预定功能的软件和/或硬件的组合。尽管以下实施例所描述的装置较佳地以软件来实现，但是硬件，或者软件和硬件的组合的实现也是可能并被构想的。图3是本发明实施例的水平裂缝预测装置的一种结构框图，如图3所示，包括：获取模块301、第一计算模块302、第二计算模块303以及处理模块，下面对该结构进行说明。
获取模块301，用于获取目标裂缝处第一入射角方向的纵波速度和横波速度和第二入射角方向的纵波速度和横波速度，其中，第一入射角小于第二入射角。
多波地震勘探是通过综合利用纵波、横波(或转换波)等多种地震波对地层进行勘探的方法，因此到多波地震资料中携带了比单波(纵波)分量更丰富的地层岩性信息，尤其还携带了地层裂缝信息。因此可以通过对获取的目标裂缝处的多波地震数据进行预处理和速度分析，得到第一入射角方向的纵波速度和横波速度、第二入射角方向的纵波速度和横波速度。其中，预处理包括但不限于：去噪、振幅补偿、增益；速度分析包括但不限于：通过速度扫描获得地震波动校正速度。
上述第一入射角小于第二入射角，具体的，上述第一入射角可以大于等于0°小于等于20°，第二入射角可以大于等于70°小于等于90°。其中，上述入射角可以为地震波传播方向与裂缝平面法线方向之间的夹角。进一步的，在地震勘探中，众多炮点(即激发地震波的震源所处位置)和接收点(即接收地震波的检波器所处位置)组成了观测系统。因此，存在不同角度的地震波入射-反射路径，上述地震波传播方向可以通过观测系统中的炮-检空间位置确定。
上述入射角的范围为0°(即垂直于裂缝面)到90°(即平行于裂缝面)，并且在多波地震数据中，各个不同角度的地震波入射-反射路径之间可以选取适当的间隔，保证在0°至20°之间以及70°至90°之间，有且仅有一个方向的多波地震数据，即上述第一入射角以及第二入射角是唯一的。
第一计算模块302，用于根据第一入射角方向的纵波速度和横波速度计算第一入射方向的第一纵横波速度比和第一横波分裂系数。
考虑到当地震横波在含有裂缝的各向异性介质中传播时，会产生横波分裂现象，形成快横波和慢横波，而且快、慢横波的差异程度与地层裂缝的发育程度成比例，横波分裂现象主要体现了裂缝信息。在一个实施例中可以包括第一获取单元，用于获取目标裂缝处多波地震数据中第一入射角方向上的快横波速度和慢横波速度；第一计算单元，用于根据第一入射角方向上的快横波速度、慢横波速度，计算第一横波分裂系数。其中上述第一横波分裂系数通过以下公式得到：
SWS1＝(Vs11-Vs21)/Vs21
其中，Vs11为第一入射角方向的快横波速度；Vs21为第一入射角方向的慢横波速度；SWS1为第一横波分裂系数。
进一步的，上述第一纵横波速度比为目标裂缝处第一入射角方向的纵波速度与横波速度相除的值。
第二计算模块303，用于根据第二入射角方向的纵波速度和横波速度计算第二入射方向的第二纵横波速度比。
上述第一纵横波速度比可以为目标裂缝处第一入射角方向的纵波速度与横波速度相除的值。
处理模块304，用于在第一横波分裂系数小于第一预设值，在第一纵横波速度比与第二纵横波速度比之间的比值小于第二预设值的情况下，确定目标裂缝为水平裂缝。
其中，上述第一预设值可以为0.02，上述第二预设值可以为0.8。当第一横波分裂系数小于0.02，并且第一纵横波速度比与第二纵横波速度比之间的比值小于0.8时，可以确定目标裂缝为水平裂缝。可以理解的是，上述第一预设值以及上述第二预设值也可以根据实际情况设定为其它合理值，本申请对此不作限定。
裂缝密度是表示裂缝发育程度的一个概念值，用于反映裂缝发育程度。根据裂缝密度可以将目标裂缝划分为裂缝发育带、次发育带或者不发育带，在一个实施例中，可以将第二横波分裂系数作为目标裂缝的裂缝密度，以确定目标裂缝的发育程度。
具体的，在确定目标裂缝为水平裂缝后，还可以包括第二获取单元，用于获取目标裂缝处多波地震数据中第二入射角方向上的快横波速度和慢横波速度；第二计算单元，用于根据第二入射角方向上的快横波速度、慢横波速度，计算第二横波分裂系数。其中上述第二横波分裂系数通过以下公式得到：
SWS2＝(Vs12-Vs22)/Vs22
其中，Vs12为第二入射角方向的快横波速度；Vs22为第二入射角方向的慢横波速度；SWS2为第二横波分裂系数。
在计算得到第二横波分裂系数后，在处理单元中将第二横波分裂系数作为目标水平裂缝的裂缝密度，通过确定单元根据目标裂缝处的裂缝密度确定目标裂缝的发育程度。
从以上的描述中，可以看出，本发明实施例实现了如下技术效果：利用地震各向异性岩石物理特征，通过计算第一入射角方向的纵横波速度比以及第一横波分裂系数、第二入射角方向的纵横波速度比，确定目标裂缝是否为水平裂缝，并通过横波分裂系数进一步确定水平裂缝的发育程度，可以有效的预测煤层、页岩等储层中的水平裂缝及其发育程度，为煤层气、页岩气的勘探开发提供了有效数据，进而有效提高了煤层气、页岩气的产量。
本申请实施方式还提供了一种电子设备，具体可以参阅图4所示的基于本申请实施例提供的水平裂缝预测方法的电子设备组成结构示意图，所述电子设备具体可以包括输入设备41、处理器42、存储器43。其中，所述输入设备41具体可以用于输入目标裂缝处第一入射角方向的纵波速度和横波速度、第二入射角方向的纵波速度和横波速度，其中，所述第一入射角小于所述第二入射角。所述处理器42具体可以用于根据所述第一入射角方向的纵波速度和横波速度计算所述第一入射方向的第一纵横波速度比和第一横波分裂系数；根据所述第二入射角方向的纵波速度和横波速度计算所述第二入射方向的第二纵横波速度比；在所述第一横波分裂系数小于第一预设值，在所述第一纵横波速度比与所述第二纵横波速度比之间的比值小于第二预设值的情况下，确定所述目标裂缝为水平裂缝。所述存储器43具体可以用于存储第一入射角方向的纵波速度、横波速度、快横波速度、慢横波速度以及第二入射角方向的纵波速度、横波速度、快横波速度、慢横波速度等参数。
在本实施方式中，所述输入设备具体可以是用户和计算机系统之间进行信息交换的主要装置之一。所述输入设备可以包括键盘、鼠标、摄像头、扫描仪、光笔、手写输入板、语音输入装置等；输入设备用于把原始数据和处理这些数的程序输入到计算机中。所述输入设备还可以获取接收其他模块、单元、设备传输过来的数据。所述处理器可以按任何适当的方式实现。例如，处理器可以采取例如微处理器或处理器以及存储可由该(微)处理器执行的计算机可读程序代码(例如软件或固件)的计算机可读介质、逻辑门、开关、专用集成电路(Application Specific Integrated Circuit，ASIC)、可编程逻辑控制器和嵌入微控制器的形式等等。所述存储器具体可以是现代信息技术中用于保存信息的记忆设备。所述存储器可以包括多个层次，在数字系统中，只要能保存二进制数据的都可以是存储器；在集成电路中，一个没有实物形式的具有存储功能的电路也叫存储器，如RAM、FIFO等；在系统中，具有实物形式的存储设备也叫存储器，如内存条、TF卡等。
在本实施方式中，该电子设备具体实现的功能和效果，可以与其它实施方式对照解释，在此不再赘述。
本申请实施方式中还提供了一种基于水平裂缝预测方法的计算机存储介质，所述计算机存储介质存储有计算机程序指令，在所述计算机程序指令被执行时实现：确定目标裂缝是否为水平裂缝，如果确定为水平裂缝，进一步确认水平裂缝的发育程度。
在本实施方式中，上述存储介质包括但不限于随机存取存储器(Random AccessMemory,RAM)、只读存储器(Read-Only Memory,ROM)、缓存(Cache)、硬盘(Hard DiskDrive,HDD)或者存储卡(Memory Card)。所述存储器可以用于存储计算机程序指令。网络通信单元可以是依照通信协议规定的标准设置的，用于进行网络连接通信的接口。
在本实施方式中，该计算机存储介质存储的程序指令具体实现的功能和效果，可以与其它实施方式对照解释，在此不再赘述。
显然，本领域的技术人员应该明白，上述的本发明实施例的各模块或各步骤可以用通用的计算装置来实现，它们可以集中在单个的计算装置上，或者分布在多个计算装置所组成的网络上，可选地，它们可以用计算装置可执行的程序代码来实现，从而，可以将它们存储在存储装置中由计算装置来执行，并且在某些情况下，可以以不同于此处的顺序执行所示出或描述的步骤，或者将它们分别制作成各个集成电路模块，或者将它们中的多个模块或步骤制作成单个集成电路模块来实现。这样，本发明实施例不限制于任何特定的硬件和软件结合。
应该理解，以上描述是为了进行图示说明而不是为了进行限制。通过阅读上述描述，在所提供的示例之外的许多实施方式和许多应用对本领域技术人员来说都将是显而易见的。因此，本申请的范围不应该参照上述描述来确定，而是应该参照前述权利要求以及这些权利要求所拥有的等价物的全部范围来确定。
以上所述仅为本发明的优选实施例而已，并不用于限制本发明，对于本领域的技术人员来说，本发明实施例可以有各种更改和变化。凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。
Claims (11)
1.一种水平裂缝预测方法，其特征在于，包括：
获取目标裂缝处第一入射角方向的纵波速度和横波速度、第二入射角方向的纵波速度和横波速度，其中，所述第一入射角小于所述第二入射角；
根据所述第一入射角方向的纵波速度和横波速度计算所述第一入射角方向的第一纵横波速度比和第一横波分裂系数；
根据所述第二入射角方向的纵波速度和横波速度计算所述第二入射角方向的第二纵横波速度比；
在所述第一横波分裂系数小于第一预设值，在所述第一纵横波速度比与所述第二纵横波速度比之间的比值小于第二预设值的情况下，确定所述目标裂缝为水平裂缝。
2.根据权利要求1所述的方法，其特征在于，根据所述第一入射角方向的纵波速度和横波速度，计算所述第一入射角方向的第一横波分裂系数，包括：
根据所述第一入射角方向的横波速度，获取所述第一入射角方向的快横波速度和慢横波速度；
根据所述第一入射角方向的快横波速度、慢横波速度，计算所述第一横波分裂系数。
3.根据权利要求2所述的方法，其特征在于，按照以下公式，根据所述第一入射角方向的快横波速度、慢横波速度，计算所述第一横波分裂系数：
SWS1＝(Vs11-Vs21)/Vs21
其中，Vs11为所述第一入射角方向的快横波速度；Vs21为所述第一入射角方向的慢横波速度；SWS1为所述第一横波分裂系数。
4.根据权利要求1所述的方法，其特征在于，在确定所述目标裂缝为水平裂缝之后，还包括：
获取所述第二入射角方向的快横波速度和慢横波速度；
根据所述第二入射角方向的快横波速度、慢横波速度计算所述第二入射角方向的第二横波分裂系数；
将所述第二横波分裂系数作为所述水平裂缝的裂缝密度；
根据所述裂缝密度，确定所述目标裂缝的发育程度。
5.根据权利要求1所述的方法，其特征在于，所述第一入射角大于等于0°小于等于20°；所述第二入射角大于等于70°小于等于90°。
6.一种水平裂缝预测装置，其特征在于，包括：
获取模块，用于获取目标裂缝处第一入射角方向的纵波速度和横波速度和第二入射角方向的纵波速度和横波速度，其中，所述第一入射角小于所述第二入射角；
第一计算模块，用于根据所述第一入射角方向的纵波速度和横波速度计算所述第一入射角方向的第一纵横波速度比和第一横波分裂系数；
第二计算模块，用于根据所述第二入射角方向的纵波速度和横波速度计算所述第二入射角方向的第二纵横波速度比；
处理模块，用于在所述第一横波分裂系数小于第一预设值，在所述第一纵横波速度比与所述第二纵横波速度比之间的比值小于第二预设值的情况下，确定所述目标裂缝为水平裂缝。
7.根据权利要求6所述的装置，其特征在于，所述第一计算模块包括：
第一获取单元，用于根据所述第一入射角方向的横波速度，获取所述第一入射角方向的快横波速度和慢横波速度；
第一计算单元，用于根据所述第一入射角方向的快横波速度、慢横波速度，计算所述第一横波分裂系数。
8.根据权利要求7所述的装置，其特征在于，所述第一计算单元按照以下公式计算所述第一横波分裂系数：
SWS1＝(Vs11-Vs21)/Vs21
其中，Vs11为所述第一入射角方向的快横波速度；Vs21为所述第一入射角方向的慢横波速度；SWS1为所述第一横波分裂系数。
9.根据权利要求6所述的装置，其特征在于，还包括：
第二获取单元，用于在确定所述目标裂缝为水平裂缝之后，获取所述第二入射角方向的快横波速度和慢横波速度；
第二计算单元，用于根据所述第二入射角方向的快横波速度、慢横波速度计算所述第二入射角方向的第二横波分裂系数；
处理单元，用于将所述第二横波分裂系数作为所述水平裂缝的裂缝密度；
确定单元，用于根据所述裂缝密度，确定所述目标裂缝的发育程度。
10.一种水平裂缝预测设备，包括处理器以及用于存储处理器可执行指令的存储器，所述处理器执行所述指令时实现权利要求1至5中任一项所述方法的步骤。
11.一种计算机可读存储介质，其上存储有计算机指令，所述指令被执行时实现权利要求1至5中任一项所述方法的步骤。
Priority Applications (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811502830.9A CN109683197B (zh) | 2018-12-10 | 2018-12-10 | 水平裂缝预测方法、装置和设备 |
| US16/536,402 US10677947B1 (en) | 2018-12-10 | 2019-08-09 | Horizontal fracture prediction method, device and equipment |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811502830.9A CN109683197B (zh) | 2018-12-10 | 2018-12-10 | 水平裂缝预测方法、装置和设备 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109683197A CN109683197A (zh) | 2019-04-26 |
| CN109683197B true CN109683197B (zh) | 2020-05-01 |
Family
ID=66187342
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811502830.9A Active CN109683197B (zh) | 2018-12-10 | 2018-12-10 | 水平裂缝预测方法、装置和设备 |
Country Status (2)
| Country | Link |
|---|---|
| US (1) | US10677947B1 (zh) |
| CN (1) | CN109683197B (zh) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112363224A (zh) * | 2020-10-29 | 2021-02-12 | 中国石油天然气集团有限公司 | 一种横波地震勘探的应用潜力判断方法及装置 |
| CN115903039B (zh) * | 2022-12-26 | 2023-10-24 | 成都理工大学 | 一种基于进退法的横波分裂分析方法 |
| CN117572495B (zh) * | 2023-10-24 | 2024-05-17 | 成都理工大学 | 一种裂缝尺度定量预测方法 |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6556921B1 (en) * | 1999-11-08 | 2003-04-29 | Board Of Regents For The University Of Texas System | Determining vertical fractures in a stratum using scattered vertical and horizontal shear modes |
| CN103076633B (zh) * | 2012-09-20 | 2014-10-15 | 中国石油天然气集团公司 | 一种检测油气储层中裂缝走向方位的方法及装置 |
| CN104407378A (zh) * | 2014-11-25 | 2015-03-11 | 中国石油天然气股份有限公司 | 一种各向异性参数反演方法及装置 |
| CN108459346A (zh) * | 2017-02-20 | 2018-08-28 | 中国石油化工股份有限公司 | 一种页岩水平层理缝密度地震预测方法 |
| CN108957548A (zh) * | 2018-05-23 | 2018-12-07 | 中国地质调查局油气资源调查中心 | 一种多波多分量联合观测地震页岩气富集区预测技术 |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2016159987A1 (en) * | 2015-03-31 | 2016-10-06 | Halliburton Energy Services, Inc. | Simplified geomechanical model of stresses on an orthorhombic media |
| GB2555365B (en) * | 2015-07-28 | 2021-07-28 | Geoquest Systems Bv | Seismic constrained discrete fracture network |
-
2018
- 2018-12-10 CN CN201811502830.9A patent/CN109683197B/zh active Active
-
2019
- 2019-08-09 US US16/536,402 patent/US10677947B1/en active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6556921B1 (en) * | 1999-11-08 | 2003-04-29 | Board Of Regents For The University Of Texas System | Determining vertical fractures in a stratum using scattered vertical and horizontal shear modes |
| CN103076633B (zh) * | 2012-09-20 | 2014-10-15 | 中国石油天然气集团公司 | 一种检测油气储层中裂缝走向方位的方法及装置 |
| CN104407378A (zh) * | 2014-11-25 | 2015-03-11 | 中国石油天然气股份有限公司 | 一种各向异性参数反演方法及装置 |
| CN108459346A (zh) * | 2017-02-20 | 2018-08-28 | 中国石油化工股份有限公司 | 一种页岩水平层理缝密度地震预测方法 |
| CN108957548A (zh) * | 2018-05-23 | 2018-12-07 | 中国地质调查局油气资源调查中心 | 一种多波多分量联合观测地震页岩气富集区预测技术 |
Non-Patent Citations (1)
| Title |
|---|
| 利用三分量VSP资料预测裂缝方向;刘俊胜;《胜利油田职工大学学报》;20060630;第20卷(第3期);全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109683197A (zh) | 2019-04-26 |
| US20200183033A1 (en) | 2020-06-11 |
| US10677947B1 (en) | 2020-06-09 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Valoroso et al. | Radiography of a normal fault system by 64,000 high‐precision earthquake locations: The 2009 L'Aquila (central Italy) case study | |
| US9659252B2 (en) | Method to characterize heterogeneous anisotropic media | |
| CN110286421B (zh) | 一种致密砂岩储层天然裂缝建模方法 | |
| Bird et al. | Plate-tectonic analysis of shallow seismicity: Apparent boundary width, beta, corner magnitude, coupled lithosphere thickness, and coupling in seven tectonic settings | |
| Rohrbaugh Jr et al. | Estimating fracture trace intensity, density, and mean length using circular scan lines and windows | |
| CN109683197B (zh) | 水平裂缝预测方法、装置和设备 | |
| Böse et al. | PreSEIS: A neural network-based approach to earthquake early warning for finite faults | |
| EP2598915B1 (en) | Inverting geophysical data for geological parameters or lithology | |
| Seher et al. | Crustal velocity structure of the Lucky Strike segment of the Mid‐Atlantic Ridge at 37 N from seismic refraction measurements | |
| Lepillier et al. | Variational Phase‐field modeling of hydraulic fracture interaction with natural fractures and application to enhanced geothermal systems | |
| WO2016108896A1 (en) | Seismic elastic wave simulation for tilted transversely isotropic media using adaptive lebedev staggered grid | |
| CN103176211B (zh) | 基于多敏感弹性参数的含气储层预测方法及装置 | |
| WO2016163984A1 (en) | Fracture-size-correlated aperture mapping for localized porosity and permeability determination | |
| Far et al. | Fracture parameter inversion for Marcellus Shale | |
| Langet et al. | Joint focal mechanism inversion using downhole and surface monitoring at the Decatur, Illinois, CO2 injection site | |
| SA521421457B1 (ar) | تحديد خصائص السمات الترسيبية من خلال تصنيف زلزالي يستند إلى الجيولوجيا | |
| Darijani et al. | Clustering and constrained inversion of seismic refraction and gravity data for overburden stripping: Application to uranium exploration in the Athabasca Basin, Canada | |
| CN115373024A (zh) | 基于地层记录沉降反演被动陆缘地壳结构的方法及装置 | |
| CN112272783A (zh) | 使用折射走时层析成像产生针对地下结构的速度模型 | |
| Olsen et al. | Constraints of crustal heterogeneity and Q (f) from regional (< 4 Hz) wave propagation for the 2009 North Korea nuclear test | |
| Kaven et al. | Geometry of crustal faults: Identification from seismicity and implications for slip and stress transfer models | |
| CN108387711B (zh) | Toc表征参数、脆性表征参数的确定方法和装置 | |
| CN107817524B (zh) | 三维地震层析成像的方法和装置 | |
| CN113960667A (zh) | 围岩质量预测方法、预测系统、计算机设备、介质及终端 | |
| US10454713B2 (en) | Domain decomposition using a multi-dimensional spacepartitioning tree |
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 |