CN110988988A - 基于垂直裂缝介质的地震波波场模拟方法及装置 - Google Patents

基于垂直裂缝介质的地震波波场模拟方法及装置 Download PDF

Info

Publication number
CN110988988A
CN110988988A CN201911169482.2A CN201911169482A CN110988988A CN 110988988 A CN110988988 A CN 110988988A CN 201911169482 A CN201911169482 A CN 201911169482A CN 110988988 A CN110988988 A CN 110988988A
Authority
CN
China
Prior art keywords
value
expression
seismic wave
grid
spatial grid
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.)
Granted
Application number
CN201911169482.2A
Other languages
English (en)
Other versions
CN110988988B (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 Mining and Technology Beijing CUMTB
Original Assignee
China University of Mining and Technology Beijing CUMTB
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 Mining and Technology Beijing CUMTB filed Critical China University of Mining and Technology Beijing CUMTB
Priority to CN201911169482.2A priority Critical patent/CN110988988B/zh
Publication of CN110988988A publication Critical patent/CN110988988A/zh
Application granted granted Critical
Publication of CN110988988B publication Critical patent/CN110988988B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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/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/30—Analysis
    • G01V1/307—Analysis for determining seismic attributes, e.g. amplitude, instantaneous phase or frequency, reflection strength or polarity

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

基于垂直裂缝介质的地震波波场模拟方法及装置
技术领域
本申请涉及地震波技术领域，具体而言，涉及一种基于垂直裂缝介质的地震波波场模拟方法及装置。
背景技术
地震波正演模拟是模拟地震波在地球介质中的传播过程，以研究地震波传播特性与地下介质参数的关系。其中，数值模拟可以用于产生波传播的可视化，从而便于直观的展示地下构造的真实形态。油气勘探实践表明，地质层裂缝是油气的储存空间或运移通道，地质层裂缝介质对地震波的影响的研究越来越受到关注。实际地质层中的裂缝形成受多种因素控制。由于上覆地层的压力使得水平或低陡倾角裂缝存在较少(即：水平裂缝)，大多数为高陡倾角裂缝(即：垂直裂缝)，因此，模拟垂直裂缝对地震波传播的影响更具意义。
申请内容
鉴于此，本申请实施例的目的在于提供一种基于垂直裂缝介质的地震波波场模拟方法及装置，以模拟地质层的垂直裂缝对地震波传播的影响。
第一方面，本申请实施例提供一种基于垂直裂缝介质的地震波波场模拟方法，所述方法包括：在已知地震波信号、震源坐标、所述地震波信号的传播速度、垂直裂缝坐标的条件下，获取用于模拟所述垂直裂缝对所述地震波信号传播的影响的数值模拟表达式，其中，所述数值模拟表达式包括：
Figure BDA0002286401520000021
Figure BDA0002286401520000022
Figure BDA0002286401520000023
Figure BDA0002286401520000031
以及
Figure BDA0002286401520000032
获取所述地质密度的值、所述拉梅系数的值、所述切向屈服度的值、所述法向屈服度的值、所述采样时刻的值、所述时间采样率的值、所述空间网格模型的大小，以及所述空间网格的长度值和宽度值；将所述地质密度的值、所述拉梅系数的值、所述切向屈服度的值、所述法向屈服度的值、所述采样时刻的值、所述时间采样率的值、所述空间网格模型中的每个空间网格点的坐标位置，以及所述空间网格的长度值和宽度值分别输入所述数值模拟表达式，得到用于模拟所述垂直裂缝对所述地震波信号传播的影响的速度值和应力值。
在上述实现过程中，在已知地震波信号、震源坐标、所述地震波信号的传播速度、垂直裂缝坐标的条件下，利用用于模拟所述垂直裂缝对所述地震波信号传播的影响的数值模拟表达式来获取用于模拟所述垂直裂缝对所述地震波信号传播的影响的速度值和应力值，继而能够根据所述速度值和所述应力值来模拟地质层的垂直裂缝对地震波传播的影响。
基于第一方面，在一种可能的设计中，获取用于模拟所述垂直裂缝对所述地震波信号传播的影响的数值模拟表达式，包括：基于线性滑动理论，获取表征所述垂直裂缝位置处的边界条件；其中，所述边界条件为：
Figure BDA0002286401520000041
基于所述边界条件和虚拟网格，获取用于表征所述地震波信号在虚网格点处速度的第一表达式；获取表征应力和所述地震波信号在所述地质中的传播速度之间的关系的有限差分表达式；基于所述虚拟网格和所述旋转交错网格对所述有限差分表达式进行离散处理，得到离散化的差分表达式；基于所述第一表达式和所述离散化的差分表达式，获取所述数值模拟表达式。
在上述实现过程中，由于对所述有限差分表达式进行离散处理的前提条件是连续，由于线性滑动理论可知，在垂直裂缝处速度是不连续的，因此，基于虚拟网格和旋转交错网格对所述有限差分表达式进行离散处理，得到离散化的差分表达式，并利用虚网格来对表征垂直裂缝位置处的边界条件处理，以得到表征所述地震波信号在虚网格点处速度的第一表达式，以显示的描述裂缝处的速度的不连续性，最后基于所述第一表达式和所述离散化的差分表达式得到所述数值模拟表达式，以解决上述问题。
基于第一方面，在一种可能的设计中，基于所述虚拟网格和所述旋转交错网格对所述有限差分表达式进行离散处理，得到离散化的差分表达式，包括：利用所述虚拟网格和所述旋转交错网格对所述有限差分表达式中的速度和应力对所述第一差分方向和所述第二差分方向的偏导分别进行离散化处理，得到第一离散化的表达式；利用旋转交错网格对所述有限差分表达式中的速度和应力对时间的偏导分别进行离散化处理，得到第二离散化的表达式；利用所述第一离散化的表达式和所述第二离散化表达式对所述有限差分表达式进行离散处理，得到所述离散化的差分表达式。
在上述实现过程中，为了方便求解所述有限差分表达式中的速度值和应力值，因此，通过上述方式对所述有限差分表达式进行离散处理，以方便求解所述有限差分表达式中的速度值和应力值。
基于第一方面，在一种可能的设计中，获取表征应力和所述地震波信号在所述地质中的传播速度之间的关系的有限差分表达式，包括：获取表征所述应力和所述地震波信号在所述地质中的传播速度之间的关系的初始有限差分表达式；利用所述第一差分方向和所述第二差分方向对所述初始有限差分表达式中的差分方向进行处理，得到所述有限差分表达式。
在上述实现过程中，为了方便计算地震波信号沿x轴和沿z轴的速度，以及沿x轴和z轴的应力，因此，利用旋转交错网格来改变所述初始有限差分表达式中的差分方向，继而方便计算地震波信号沿x轴和沿z轴的速度，以及沿x轴和z轴的应力。
第二方面，本申请实施例提供一种基于垂直裂缝介质的地震波波场模拟方法，所述方法包括：在已知地震波信号、震源坐标、所述地震波信号的传播速度、垂直裂缝坐标的条件下，获取用于模拟所述垂直裂缝对所述地震波信号传播的影响的新的数值模拟表达式，其中，所述新的数值模拟表达式包括：
Figure BDA0002286401520000061
Figure BDA0002286401520000062
Figure BDA0002286401520000063
Figure BDA0002286401520000071
Figure BDA0002286401520000072
Figure BDA0002286401520000081
以及
Figure BDA0002286401520000091
获取所述地质密度的值、所述拉梅系数的值、所述切向屈服度的值、所述法向屈服度的值、所述采样时刻的值、所述时间采样率的值、所述空间网格模型的大小，以及所述空间网格的长度值和宽度值；将所述地质密度的值、所述拉梅系数的值、所述切向屈服度的值、所述法向屈服度的值、所述采样时刻的值、所述时间采样率的值、所述空间网格模型中的每个空间网格点的坐标位置，以及所述空间网格的长度值和宽度值分别输入所述新的数值模拟表达式，得到用于模拟所述垂直裂缝对所述地震波信号传播的影响的速度值和应力值。
在上述实现过程中，在已知地震波信号、震源坐标、所述地震波信号的传播速度、垂直裂缝坐标的条件下，利用用于模拟所述垂直裂缝对所述地震波信号传播的影响的新的数值模拟表达式来获取用于模拟所述垂直裂缝对所述地震波信号传播的影响的速度值和应力值，继而能够根据所述速度值和所述应力值来模拟地质层的垂直裂缝对地震波传播的影响。
第三方面，本申请实施例提供一种基于垂直裂缝介质的地震波波场模拟装置，所述装置包括：第一模拟表达式获取单元，用于在已知地震波信号、震源坐标、所述地震波信号的传播速度、垂直裂缝坐标的条件下，获取用于模拟所述垂直裂缝对所述地震波信号传播的影响的数值模拟表达式；第一参数值获取单元，获取所述地质密度的值、所述拉梅系数的值、所述切向屈服度的值、所述法向屈服度的值、所述采样时刻的值、所述时间采样率的值、所述空间网格模型的大小，以及所述空间网格的长度值和宽度值；数值模拟单元，用于将所述地质密度的值、所述拉梅系数的值、所述切向屈服度的值、所述法向屈服度的值、所述采样时刻的值、所述时间采样率的值、所述空间网格模型中的每个空间网格点的坐标位置，以及所述空间网格的长度值和宽度值分别输入所述数值模拟表达式，得到用于模拟所述垂直裂缝对所述地震波信号传播的影响的速度值和应力值。
第四方面，本申请实施例提供一种基于垂直裂缝介质的地震波波场模拟装置，所述装置包括：第二模拟表达式获取单元，用于在已知地震波信号、震源坐标、所述地震波信号的传播速度、垂直裂缝坐标的条件下，获取用于模拟所述垂直裂缝对所述地震波信号传播的影响的新的数值模拟表达式，其中，所述新的数值模拟表达式包括；第二参数值获取单元，用于获取所述地质密度的值、所述拉梅系数的值、所述切向屈服度的值、所述法向屈服度的值、所述采样时刻的值、所述时间采样率的值、所述空间网格模型的大小，以及所述空间网格的长度值和宽度值；第二数值模拟单元，用于将所述地质密度的值、所述拉梅系数的值、所述切向屈服度的值、所述法向屈服度的值、所述采样时刻的值、所述时间采样率的值、所述空间网格模型中的每个空间网格点的坐标位置，以及所述空间网格的长度值和宽度值分别输入所述新的数值模拟表达式，得到用于模拟所述垂直裂缝对所述地震波信号传播的影响的速度值和应力值。
第五方面，本申请实施例提供一种电子设备，包括处理器以及与所述处理器连接的存储器，所述存储器内存储计算机程序，当所述计算机程序被所述处理器执行时，使得所述电子设备执行第一方面和第二方面所述的方法。
第六方面，本申请实施例提供一种存储介质，所述存储介质中存储有计算机程序，当所述计算机程序在计算机上运行时，使得所述计算机执行第一方面和第二方面所述的方法。
本申请的其他特征和优点将在随后的说明书阐述，并且，部分地从说明书中变得显而易见，或者通过实施本申请实施例了解。本申请的目的和其他优点可通过在所写的说明书、权利要求书、以及附图中所特别指出的结构来实现和获得。
附图说明
为了更清楚地说明本申请实施例的技术方案，下面将对实施例中所需要使用的附图作简单地介绍，应当理解，以下附图仅示出了本申请的某些实施例，因此不应被看作是对范围的限定，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他相关的附图。
图1为本申请实施例提供的基于垂直裂缝介质的地震波波场模拟方法的流程示意图。
图2为本申请实施例提供的地震波和垂直裂缝的场景图。
图3为本申请实施例提高的旋转交错网格的网格示意图。
图4为本申请实施例提供的虚拟网格示意图。
图5为本申请实施例提供的表征虚网格和实网格重叠的示意图。
图6为本申请实施例提供的另一种基于垂直裂缝介质的地震波波场模拟方法的流程图。
图7为本申请实施例提供的网格模型中的网格点的示意图。
图8为本申请实施例提供的基于垂直裂缝介质的地震波波场模拟装置的结构示意图。
图9为本申请实施例提供的另一种基于垂直裂缝介质的地震波波场模拟装置的结构示意图。
图10为本申请实施例提供的地震波波场的一种传播示意图。
图11为本申请实施例提供的地震波波场的另一种传播示意图。
图12为本申请实施例提供的电子设备的结构示意图。
具体实施方式
下面将结合本申请实施例中附图，对本申请实施例中的技术方案进行描述。
应注意到：相似的标号和字母在下面的附图中表示类似项，因此，一旦某一项在一个附图中被定义，则在随后的附图中不需要对其进行进一步定义和解释。同时，在本申请的描述中，术语“第一”、“第二”等仅用于区分描述，而不能理解为指示或暗示相对重要性。
请参照图1，图1为本申请实施例提供的一种基于垂直裂缝介质的地震波波场模拟方法，所述方法包括步骤：S100、S200和S300。
S100：在已知地震波信号、震源坐标、所述地震波信号的传播速度、垂直裂缝坐标的条件下，获取用于模拟所述垂直裂缝对所述地震波信号传播的影响的数值模拟表达式，其中，所述数值模拟表达式包括：
Figure BDA0002286401520000131
Figure BDA0002286401520000132
Figure BDA0002286401520000133
Figure BDA0002286401520000141
以及
Figure BDA0002286401520000142
其中，t表示采样时刻，(m，n)表示空间网格点在空间网格模型中的坐标位置，Δt表示所述时间采样率，Δx表示所述空间网格模型中的空间网格的长度，Δz表示所述空间网格的宽度，NT表示所述垂直裂缝的切向屈服度、NS表示所述垂直裂缝的法向屈服度，λ、μ表示拉梅系数，ρ表示地质密度，
Figure BDA0002286401520000151
表示所述地震波信号在采样时刻为t+1/2、空间网格点的坐标位置为(m，n)的沿x轴方向的速度，
Figure BDA0002286401520000152
表示所述地震波信号在采样时刻为t-1/2、空间网格点的坐标位置为(m，n)的沿z轴方向的速度，
Figure BDA0002286401520000153
表示所述地震波信号在采样时刻为t+1、空间网格点的坐标位置为(m+1/2，n+1/2)的沿x轴方向的应力，
Figure BDA0002286401520000154
表示所述地震波信号在采样时刻为t、网格点的坐标位置为(m+1/2，n+1/2)的沿z轴方向的应力，
Figure BDA0002286401520000155
表示所述地震波信号在采样时刻为t、网格点的坐标位置为(m+1/2，n+1/2)的沿xz轴方向的应力。
S200：获取所述地质密度的值、所述拉梅系数的值、所述切向屈服度的值、所述法向屈服度的值、所述采样时刻的值、所述时间采样率的值、所述空间网格模型的大小，以及所述空间网格的长度值和宽度值；
S300：将所述地质密度的值、所述拉梅系数的值、所述切向屈服度的值、所述法向屈服度的值、所述采样时刻的值、所述时间采样率的值、所述空间网格模型中的每个空间网格点的坐标位置，以及所述空间网格的长度值和宽度值分别输入所述数值模拟表达式，得到用于模拟所述垂直裂缝对所述地震波信号传播的影响的速度值和应力值。
下面对上述方法进行详细阐述。
在本实施例中，所述地震波信号为雷克子波，在其他实施例中所述地震波信号可以为高斯子波等信号波、瑞克子波等，请参考图2，在本实施例中，建立包含x轴和z轴的二维坐标系，获取所述地震波信号开始传输的位置在所述二维坐标系中的坐标位置，即震源坐标(图中的点)，以及获取垂直裂缝(图中的线)在所述二维坐标系中的坐标位置，在本实施例中，所述震源坐标为(x＝750m，z＝750m)，所述垂直裂缝的坐标为(x＝1000m)，所述震源坐标和所述垂直裂缝的坐标根据实际需求设定，在此不做限制；在本实施例中，所述地震波信号的传播速度包括地震波的横波速度和地震波的纵波速度，所述地震波的横波速度为1400m/s，地震波的纵波速度为2800m/s，在其他实施例中，所述地震波信号的传播速度根据实际需求设定，在此不做限制。
作为一种实施方式，S100包括：获取表征所述应力和所述地震波信号在所述地质中的传播速度之间的关系的初始有限差分表达式；其中，所述初始有限差分表达式为：
Figure BDA0002286401520000161
其中，vx表示在地质中传播的地震波信号沿x轴方向的速度分量，vz表示表示在地质中传播的地震波信号沿z轴方向的速度分量，τxx表示沿x轴方向的应力分量，τxz表示沿xz轴方向的应力分量，
Figure BDA0002286401520000162
表示vx对时间t的偏导、
Figure BDA0002286401520000163
表示τxx对t的偏导，
Figure BDA0002286401520000164
表示应力τxz对x的偏导。由于获取所述初始有限差分表达式的方式为本领域技术人员所熟知的内容，因此在此不再赘述。
为了求解所述初始有限差分表达式中的速度和应力的值，因此，在获取到所述初始有限差分表达式之后，利用第一差分方向和第二差分方向对所述初始有限差分表达式中的差分方向进行处理，获得表征应力和所述地震波信号在所述地质中的传播速度之间的关系的有限差分表达式；其中，所述有限差分表达式为：
Figure BDA0002286401520000171
其中，
Figure BDA0002286401520000172
Figure BDA0002286401520000173
为旋转交错网格的第一差分方向；
Figure BDA0002286401520000174
为所述旋转交错网格的第二差分方向。
在实际实施过程中，获得表征应力和所述地震波信号在所述地质中的传播速度之间的关系的有限差分表达式的步骤包括：A、B和C。
A：利用旋转交错网格定义所述第一差分方向和所述第二差分方向。请参照图3，下面对所述旋转交错网格的第一差分方向和所述旋转交错网格的差分方向进行介绍：
本申请通过在预先建立的二维坐标系中建立网格模型，所述网格模型中有多个网格，定义所述第一差分方向
Figure BDA0002286401520000175
和所述第二差分方向
Figure BDA0002286401520000176
分别为沿着所述网格模型中的网格的两条对角线的上的方向，在所述网格模型大小相同的条件下，所述网格模型的长度和宽度的值越小，网格的数量越多，相反，网格的数量越少，网格的数量越多，因此，能够获取到的速度值和应力值越多，更有利于研究垂直裂缝对地震波的影响。
作为一种实施方式，第一差分方向也可以为图3中的
Figure BDA0002286401520000177
方向，第二差分方向也可以为图3中的
Figure BDA0002286401520000178
方向。
B：在定义了所述第一差分方向和所述第二差分方向之后，利用所述第一差分方向和所述第二方向的线性组合来重新表示所述初始有限差分方程中的差分方向，得到新的差分方向表达式
Figure BDA0002286401520000181
C：利用所述新的差分方向表达式替换所述初始有限差分表示式中的差分方向，获得表征应力和所述地震波信号在所述地质中的传播速度之间的关系的有限差分表达式。
为了求解所述有限差分表达式中的速度和应力值，因此，本申请需要利用差分算子代替微分算子，把计算区域离散化，以求取该区域每个网格节点在不同时刻的速度值和应力值。
下面对所述有限差分表达的离散过程进行详细阐述。
请参照图4，图4中的垂直裂缝的左侧被设置为实网格，右侧被设置为虚网格，图4中的
Figure BDA0002286401520000182
为所述地震波信号在虚网格点坐标位置为(m+1，n)的速度；(v)m,n为所述地震波信号在实网格点坐标位置为(m，n)的速度。
利用所述虚拟网格和所述旋转交错网格对所述有限差分表达式中的速度和应力对所述第一差分方向和所述第二差分方向的偏导分别进行离散化处理，得到第一离散化的表达式。
其中，所述第一离散化的表达式包括：
Figure BDA0002286401520000183
其中，
Figure BDA0002286401520000184
表示网格点位置为(m，n+1/2)处的速度(vi)m,n+1/2对
Figure BDA0002286401520000185
的偏导，
Figure BDA0002286401520000186
表示网格点位置为(m，n+1/2)处的速度(vi)m,n+1/2对
Figure BDA0002286401520000195
的偏导，
Figure BDA0002286401520000191
表示网格点位置为(m，n+1/2)处的应力(τij)m,n+1/2对
Figure BDA0002286401520000196
的偏导，其中，vi包括vx和vz两种情况。
利用旋转交错网格对所述有限差分表达式中的速度和应力对时间的偏导分别进行离散化处理，得到第二离散化的表达式。其中，所述第二离散化的表达式包括：
Figure BDA0002286401520000192
其中，i和j均可以为x和z。例如，τij可以表示τxz也可以表示τxx或τzz，
Figure BDA0002286401520000193
表示表示网格点位置为(m,n)处的应力(τij)m,n在第t时刻对时间的偏导。
在获取到所述第一离散化表达式和所述第二离散化表达式之后，将所述第一离散化的表达式和所述第二离散化表达式代入所述有限差分表达式中以进行离散处理，得到离散化的差分表达式。
其中，所述离散化的差分表达式包括：
Figure BDA0002286401520000194
Figure BDA0002286401520000201
其中，Δt为时间采样率。其中，时间采用率Δt越小，模拟精度越高，计算复杂的也越大。
基于线性滑动理论，获取表征所述垂直裂缝位置处的边界条件；其中，所述边界条件为：
Figure BDA0002286401520000202
和
Figure BDA0002286401520000203
其中，“+”表示所述垂直裂缝的右侧，“-”表示所述垂直裂缝的左侧。
其中，NT表示垂直裂缝的切向屈服度，Ns表示垂直裂缝的法向屈服度，vx +表示地震波信号在垂直裂缝的左侧沿x轴方向的传播速度，vx -表示地震波信号在垂直裂缝的右侧沿x轴方向的传播速度，vz +表示地震波信号在垂直裂缝的左侧沿z轴方向的传播速度，vz -表示地震波信号在垂直裂缝的右侧沿x轴方向的传播速度，
Figure BDA0002286401520000204
表示在垂直裂缝的左侧沿x轴方向的应力τxx +对时间t的偏导；
Figure BDA0002286401520000205
表示在垂直裂缝的右侧沿xz轴方向的应力τxz -对时间t的偏导。
基于所述边界条件和虚拟网格，获取用于表征所述地震波信号在虚网格点处速度的第一表达式，其中，所述第一表达式包括：
Figure BDA0002286401520000211
其中，
Figure BDA0002286401520000212
为所述地震波信号在虚网格点坐标位置为(m，n)的沿x轴方向的速度。
其中，获取用于表征所述地震波信号在虚网格点处速度的第一表达式的步骤可以按照如下方式实施：
请参照图5，由于实网格和虚网格重叠在一起，因此，网格模型中的每个网格点的速度即可以用虚网格处的速度表示又可以用实网格处的速度表示，其中，图5中的裂缝的右侧用“+”表示，左侧用“-”表示，当需要计算裂缝处(即x方向)的速度时，从右侧看，速度从裂缝(图5中的波浪线)右侧传递到左侧时，由于地震波信号需要穿过裂缝，所以速度不连续，因此需要利用到虚网格处的速度，可以看作是从实网格到虚网格，反方向类似，继而根据中心差分准则，得到切向的速度差分表达式
Figure BDA0002286401520000213
根据算术平均法则，得到速度平均表达式
Figure BDA0002286401520000214
在需要计算裂缝左右两侧(即z方向)的速度时，由于z方向没有裂缝，因此，地震波信号在z方向的速度为连续的，不需要利用虚网格，继而根据中心差分准则，得到法向的速度差分表达式
Figure BDA0002286401520000215
由于所述初始有限差分表达式中的切向和法向的速度-应力表达式为：
Figure BDA0002286401520000221
继而通过将所述切向的速度差分表达式、所述速度平均表达式、法向的速度差分表达式以及所述速度-应力表达式，分别代入所述边界条件中经过简单的化简可得表征所述地震波信号在虚网格点处速度的第一表达式。
在获取到了所述第一表达式和所述离散化的差分表达式之后，将所述第一表达式代入所述离散化的差分表达式中进行化简处理，获取所述数值模拟表达式。
S200：获取所述地质密度的值、所述拉梅系数的值、所述切向屈服度的值、所述法向屈服度的值、所述采样时刻的值、所述时间采样率的值、所述空间网格模型的大小，以及所述空间网格的长度值和宽度值。
其中，所述地质密度的值ρ＝2600kg/m3，所述地质密度的值根据实际情况设定，所述拉梅系数的值λ＝0.2、μ＝0.8，拉梅系数的取值为经验值，在本实施例中，所述切向屈服度的值NT＝1.27*10-9，所述法向屈服度的值Ns＝2.69-9，NT和Ns的取值为经验值，在其他实施例中，NT和Ns的取值也可以为其他值，其中，所述空间网格模型的大小为1500m*1500m，所述网格模型中的每个空间网格的长度值为5m，宽度值为5m，其中，网格的长度值和宽度值不一定相等，在其他实施例中，所述网格模型的大小以及网格的长度值和宽度值也可以设置为其他值，例如，网格长度值为6m，宽度值为5m；在本实施例中，采样时刻t＝40ms，采样时刻也可以为20ms，30ms等，时间采用率Δt＝0.0001S，在其他实施例中，时间采样率的值根据实际需求设置，其中，时间采用率Δt越小，模拟精度越高，计算复杂的也越大。
S300：将所述地质密度的值、所述拉梅系数的值、所述切向屈服度的值、所述法向屈服度的值、所述采样时刻的值、所述时间采样率的值、所述空间网格模型中的每个空间网格点的坐标位置，以及所述空间网格的长度值和宽度值分别输入所述数值模拟表达式，得到用于模拟所述垂直裂缝对所述地震波信号传播的影响的速度值和应力值。
请参照图6，本申请实施例提供了另一种基于垂直裂缝介质的地震波波场模拟方法，所述方法包括：S400、S500和S600。
S400：在已知地震波信号、震源坐标、所述地震波信号的传播速度、垂直裂缝坐标的条件下，获取用于模拟所述垂直裂缝对所述地震波信号传播的影响的新的数值模拟表达式，其中，所述新的数值模拟表达式包括：
Figure BDA0002286401520000231
Figure BDA0002286401520000232
Figure BDA0002286401520000233
Figure BDA0002286401520000241
Figure BDA0002286401520000242
Figure BDA0002286401520000251
以及
Figure BDA0002286401520000252
Figure BDA0002286401520000261
其中，t表示采样时刻，(m，n)表示空间网格点在空间网格模型中的坐标位置，Δt表示所述时间采样率，Δx表示所述空间网格模型中的空间网格的长度，Δz表示所述空间网格的宽度，NT表示所述垂直裂缝的切向屈服度、NS表示所述垂直裂缝的法向屈服度，λ、μ表示拉梅系数，ρ表示地质密度，
Figure BDA0002286401520000262
表示所述地震波信号在采样时刻为t+1/2、空间网格点的坐标位置为(m，n)的沿x轴方向的速度，
Figure BDA0002286401520000263
表示所述地震波信号在采样时刻为t-1/2、空间网格点的坐标位置为(m，n)的沿z轴方向的速度，
Figure BDA0002286401520000264
表示所述地震波信号在采样时刻为t+1、空间网格点的坐标位置为(m+1/2，n+1/2)的沿x轴方向的应力，
Figure BDA0002286401520000265
表示所述地震波信号在采样时刻为t、网格点的坐标位置为(m+1/2，n+1/2)的沿z轴方向的应力，
Figure BDA0002286401520000266
表示所述地震波信号在采样时刻为t、网格点的坐标位置为(m+1/2，n+1/2)的沿xz轴方向的应力；C1和C2为有限差分系数。
S500：获取所述地质密度的值、所述拉梅系数的值、所述切向屈服度的值、所述法向屈服度的值、所述采样时刻的值、所述时间采样率的值、所述空间网格模型的大小，以及所述空间网格的长度值和宽度值。
S600：将所述地质密度的值、所述拉梅系数的值、所述切向屈服度的值、所述法向屈服度的值、所述采样时刻的值、所述时间采样率的值、所述空间网格模型中的每个空间网格点的坐标位置，以及所述空间网格的长度值和宽度值分别输入所述新的数值模拟表达式，得到用于模拟所述垂直裂缝对所述地震波信号传播的影响的速度值和应力值。
利用有限差分进行数值模拟的原理是利用差分算子代替微分算子，把计算区域离散化，求取该区域每个网格节点在不同时刻的波场值。但是，当空间网格较大、震源主频过高时，会引起波形的畸变，使波场发生较严重的数值频散。有限差分计算中数值频散压制思路主要为提高有限差分阶数，阶数越高数值模拟精度越高，但提高差分阶数会增加计算量。其主要原理是数值离散时，使用泰勒展开，根据泰勒展开，使用的阶数越高，精度越高。
因此，为了抑制数值频散，以进一步提高数值模拟精度，因此，S400中获取所述新的数值模拟表达式的方式为：与S100不同的是，在获取到所述第一离散化表达式和所述第二离散化表达式之后，将所述第一离散化的表达式和所述第二离散化表达式带入所述有限差分表达式中以进行离散处理的过程中，在本申请实施例中，利用旋转交错网格的空间四阶来表示所述有限差分表达式中速度和应力对于第一差分方向和第二差分方向的偏导，请参照图7，在需要计算黑色方块处的速度时，利用黑色圆点处的速度来计算黑色方块处的速度，图7中的网格点包括虚拟网格点和实网格点。
在需要计算黑色方块处的速度时，若采用旋转交错网格的空间二阶的表示方式，则需要利用黑色方框附近的四个黑点的速度v(m，n)、v(m+1，n)、v(m，n+1)和v(m+1，n+1)，继而利用v(m+1，n+1)-v(m，n)代表
Figure BDA0002286401520000271
方向的二阶空间差分，另外两个相减代表
Figure BDA0002286401520000272
方向的二阶空间差分。若采用旋转交错网格的空间四阶的表示方式，则还需要利用到外面的四个黑点处的速度v(m+2，n+2)、v(m-1，n-1)、v(m+2，n-1)和v(m-1，n+2)，继而可以利用(v(m+1，n+1)-v(m，n))+(v(m+2，n+2)-v(m-1，n-1))代表
Figure BDA0002286401520000283
方向的四阶空间差分，利用上述方式对所述有限差分表达式进行离散化处理，继而得到新的离散化的差分表达式：
Figure BDA0002286401520000281
表1.几种有限差分系数
Figure BDA0002286401520000282
从表1可知，由于本申请是利用空间四阶来表示所述有限差分表达式中速度和应力对于第一差分方向和第二差分方向的偏导，因此，C1＝1.125，C2＝-0.0416666。
通过将S100中获取到所述地震波信号在虚网格点处速度的第一表达式代入所述新的离散化的差分表达式，得到所述新的数值模拟表达式。其中，S500和S600的具体实施方式与S200和S300相同，因此在此不再赘述。
请参照图8，本申请实施例提供一种基于垂直裂缝介质的地震波波场模拟装置，所述装置包括：第一模拟表达式获取单元410，用于在已知地震波信号、震源坐标、所述地震波信号的传播速度、垂直裂缝坐标的条件下，获取用于模拟所述垂直裂缝对所述地震波信号传播的影响的数值模拟表达式。第一参数值获取单元420，获取所述地质密度的值、所述拉梅系数的值、所述切向屈服度的值、所述法向屈服度的值、所述采样时刻的值、所述时间采样率的值、所述空间网格模型的大小，以及所述空间网格的长度值和宽度值。数值模拟单元430，用于将所述地质密度的值、所述拉梅系数的值、所述切向屈服度的值、所述法向屈服度的值、所述采样时刻的值、所述时间采样率的值、所述空间网格模型中的每个空间网格点的坐标位置，以及所述空间网格的长度值和宽度值分别输入所述数值模拟表达式，得到用于模拟所述垂直裂缝对所述地震波信号传播的影响的速度值和应力值。
作为一种实施方式，所述第一模拟表达式获取单元410，包括：边界条件获取单元，用于基于线性滑动理论，获取表征所述垂直裂缝位置处的边界条件；虚拟网格点处的速度获取单元，用于基于所述边界条件和虚拟网格，获取用于表征所述地震波信号在虚网格点处速度的第一表达式；差分表达式获取单元，用于获取表征应力和所述地震波信号在所述地质中的传播速度之间的关系的有限差分表达式；离散处理单元，用于基于所述虚拟网格和所述旋转交错网格对所述有限差分表达式进行离散处理，得到离散化的差分表达式；模拟表达式获取子单元，用于基于所述第一表达式和所述离散化的差分表达式，获取所述数值模拟表达式。
作为一种实施方式，所述离散处理单元，还用于利用所述虚拟网格和所述旋转交错网格对所述有限差分表达式中的速度和应力对所述第一差分方向和所述第二差分方向的偏导分别进行离散化处理，得到第一离散化的表达式；利用旋转交错网格对所述有限差分表达式中的速度和应力对时间的偏导分别进行离散化处理，得到第二离散化的表达式；以及利用所述第一离散化的表达式和所述第二离散化表达式对所述有限差分表达式进行离散处理，得到所述离散化的差分表达式。
作为一种实施方式，所述差分表达式获取单元，还用于获取表征所述应力和所述地震波信号在所述地质中的传播速度之间的关系的初始有限差分表达式；以及利用所述第一差分方向和所述第二差分方向对所述初始有限差分表达式中的差分方向进行处理，得到所述有限差分表达式。
请参照图9，本申请实施例提供另一种基于垂直裂缝介质的地震波波场模拟装置，所述装置包括：第二模拟表达式获取单元510，用于在已知地震波信号、震源坐标、所述地震波信号的传播速度、垂直裂缝坐标的条件下，获取用于模拟所述垂直裂缝对所述地震波信号传播的影响的新的数值模拟表达式。第二参数值获取单元520，用于获取所述地质密度的值、所述拉梅系数的值、所述切向屈服度的值、所述法向屈服度的值、所述采样时刻的值、所述时间采样率的值、所述空间网格模型的大小，以及所述空间网格的长度值和宽度值；第二数值模拟单元530，用于将所述地质密度的值、所述拉梅系数的值、所述切向屈服度的值、所述法向屈服度的值、所述采样时刻的值、所述时间采样率的值、所述空间网格模型中的每个空间网格点的坐标位置，以及所述空间网格的长度值和宽度值分别输入所述新的数值模拟表达式，得到用于模拟所述垂直裂缝对所述地震波信号传播的影响的速度值和应力值。
本实施例对的各功能单元实现各自功能的过程，请参见上述图1-图7所示实施例中描述的内容，此处不再赘述。
请参照图10和图11，图10为利用本申请所提供的方法在仿真实验中得到的地震波的垂直分量的传播示意图，图11为利用本申请所提供的方法在仿真实验中得到的地震波的水平分量的传播示意图，其中，仿真参数为：地震波采用主频为f＝30hz的瑞克子波；该地震波的纵波波速为2800m/s；该地震波的横波速度为1400m/s；地质密度ρ＝2600kg/m3。模型的空间网格参数Δx＝Δz＝5m，时间步长是Δt＝0.0001s；切向屈服度NT＝1.27×10-9m/Pa，法向屈服度Ns＝2.69×10-9m/Pa。
请参照图12，图为本申请实施例提供的一种电子设备100的结构示意图，电子设备100可以包括：存储器102、处理器101和通信总线，通信总线用于实现这些组件的连接通信。电子设备100可以为电脑、智能手机、个人数字助理(personal digital assistant，PDA)等。
存储器102用于存储所述地质密度的值、所述拉梅系数的值、所述切向屈服度的值、所述法向屈服度的值、所述采样时刻的值、所述时间采样率的值、所述空间网格模型的大小，以及所述空间网格的长度值和宽度值，所述数值模拟表达式，所述新的数值模拟表达式以及本申请实施例提供的基于垂直裂缝介质的地震波波场模拟方法及装置对应的计算指令，其中，存储器102可以是，但不限于，随机存取存储器(Random Access Memory，RAM)，只读存储器(Read Only Memory，ROM)等。
其中，处理器101用于读取并运行存储于存储器102中的计算机指令时，执行本申请实施例所提供的基于垂直裂缝介质的地震波波场模拟方法的步骤。处理器101可能是一种集成电路芯片，具有信号的处理能力。上述的处理器101可以是通用处理器，包括中央处理器(Central Processing Unit，CPU)、网络处理器(Network Processor，NP)等。
此外，本申请实施例还提供了一种存储介质，在该存储介质中存储有计算机程序，当所述计算机程序在计算机上运行时，使得所述计算机执行本申请任一项实施方式所提供的方法。
综上所述，本申请各实施例提出的基于垂直裂缝介质的地震波波场模拟方法及装置，在已知地震波信号、震源坐标、所述地震波信号的传播速度、垂直裂缝坐标的条件下，利用用于模拟所述垂直裂缝对所述地震波信号传播的影响的数值模拟表达式来获取用于模拟所述垂直裂缝对所述地震波信号传播的影响的速度值和应力值，继而能够根据所述速度值和所述应力值来模拟地质层的垂直裂缝对地震波传播的影响。
在本申请所提供的实施例中，应该理解到，所揭露的装置和方法，也可以通过其它的方式实现。以上所描述的装置实施例仅仅是示意性的，例如，附图中的流程图和框图显示了根据本申请的多个实施例的装置和方法。在这点上，流程图或框图中的每个方框可以代表一个模块、程序段或代码的一部分，所述模块、程序段或代码的一部分包含一个或多个用于实现规定的逻辑功能的可执行指令。也应当注意，在有些作为替换的实现方式中，方框中所标注的功能也可以以不同于附图中所标注的顺序发生。例如，两个连续的方框实际上可以基本并行地执行，它们有时也可以按相反的顺序执行，这依所涉及的功能而定。也要注意的是，框图和/或流程图中的每个方框、以及框图和/或流程图中的方框的组合，可以用执行规定的功能或动作的专用的基于硬件的装置来实现，或者可以用专用硬件与计算机指令的组合来实现。
以上所述仅为本申请的优选实施例而已，并不用于限制本申请，对于本领域的技术人员来说，本申请可以有各种更改和变化。凡在本申请的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本申请的保护范围之内。

Claims (10)

1.一种基于垂直裂缝介质的地震波波场模拟方法，其特征在于，所述方法包括：
在已知地震波信号、震源坐标、所述地震波信号的传播速度、垂直裂缝坐标的条件下，获取用于模拟所述垂直裂缝对所述地震波信号传播的影响的数值模拟表达式，其中，所述数值模拟表达式包括：
Figure FDA0002286401510000011
Figure FDA0002286401510000012
Figure FDA0002286401510000013
Figure FDA0002286401510000021
以及
Figure FDA0002286401510000022
其中，t表示采样时刻，(m，n)表示空间网格点在空间网格模型中的坐标位置，Δt表示所述时间采样率，Δx表示所述空间网格模型中的空间网格的长度，Δz表示所述空间网格的宽度，NT表示所述垂直裂缝的切向屈服度、NS表示所述垂直裂缝的法向屈服度，λ、μ表示拉梅系数，ρ表示地质密度，
Figure FDA0002286401510000031
表示所述地震波信号在采样时刻为t+1/2、空间网格点的坐标位置为(m，n)的沿x轴方向的速度，
Figure FDA0002286401510000032
表示所述地震波信号在采样时刻为t-1/2、空间网格点的坐标位置为(m，n)的沿z轴方向的速度，
Figure FDA0002286401510000033
表示所述地震波信号在采样时刻为t+1、空间网格点的坐标位置为(m+1/2，n+1/2)的沿x轴方向的应力，
Figure FDA0002286401510000034
表示所述地震波信号在采样时刻为t、网格点的坐标位置为(m+1/2，n+1/2)的沿z轴方向的应力，
Figure FDA0002286401510000035
表示所述地震波信号在采样时刻为t、网格点的坐标位置为(m+1/2，n+1/2)的沿xz轴方向的应力；
获取所述地质密度的值、所述拉梅系数的值、所述切向屈服度的值、所述法向屈服度的值、所述采样时刻的值、所述时间采样率的值、所述空间网格模型的大小，以及所述空间网格的长度值和宽度值；
将所述地质密度的值、所述拉梅系数的值、所述切向屈服度的值、所述法向屈服度的值、所述采样时刻的值、所述时间采样率的值、所述空间网格模型中的每个空间网格点的坐标位置，以及所述空间网格的长度值和宽度值分别输入所述数值模拟表达式，得到用于模拟所述垂直裂缝对所述地震波信号传播的影响的速度值和应力值。
2.根据权利要求1所述的方法，其特征在于，获取用于模拟所述垂直裂缝对所述地震波信号传播的影响的数值模拟表达式，包括：
基于线性滑动理论，获取表征所述垂直裂缝位置处的边界条件；其中，所述边界条件为：
Figure FDA0002286401510000036
和
Figure FDA0002286401510000037
其中，“+”表示所述垂直裂缝的右侧，“-”表示所述垂直裂缝的左侧；
基于所述边界条件和虚拟网格，获取用于表征所述地震波信号在虚网格点处速度的第一表达式，其中，所述第一表达式包括：
Figure FDA0002286401510000041
其中，
Figure FDA0002286401510000042
为所述地震波信号在虚网格点坐标位置为(m，n)的沿x轴方向的速度；
获取表征应力和所述地震波信号在所述地质中的传播速度之间的关系的有限差分表达式；其中，所述有限差分表达式为：
Figure FDA0002286401510000043
其中，
Figure FDA0002286401510000044
Figure FDA0002286401510000045
为旋转交错网格的第一差分方向；
Figure FDA0002286401510000046
为所述旋转交错网格的第二差分方向；
基于所述虚拟网格和所述旋转交错网格对所述有限差分表达式进行离散处理，得到离散化的差分表达式；其中，所述离散化的差分表达式包括：
Figure FDA0002286401510000051
基于所述第一表达式和所述离散化的差分表达式，获取所述数值模拟表达式。
3.根据权利要求2所述的方法，其特征在于，基于所述虚拟网格和所述旋转交错网格对所述有限差分表达式进行离散处理，得到离散化的差分表达式，包括：
利用所述虚拟网格和所述旋转交错网格对所述有限差分表达式中的速度和应力对所述第一差分方向和所述第二差分方向的偏导分别进行离散化处理，得到第一离散化的表达式；
利用旋转交错网格对所述有限差分表达式中的速度和应力对时间的偏导分别进行离散化处理，得到第二离散化的表达式；
利用所述第一离散化的表达式和所述第二离散化表达式对所述有限差分表达式进行离散处理，得到所述离散化的差分表达式。
4.根据权利要求2所述的方法，其特征在于，获取表征应力和所述地震波信号在所述地质中的传播速度之间的关系的有限差分表达式，包括：
获取表征所述应力和所述地震波信号在所述地质中的传播速度之间的关系的初始有限差分表达式；其中，所述初始有限差分表达式为：
Figure FDA0002286401510000061
利用所述第一差分方向和所述第二差分方向对所述初始有限差分表达式中的差分方向进行处理，得到所述有限差分表达式。
5.一种基于垂直裂缝介质的地震波波场模拟方法，其特征在于，所述方法包括：
在已知地震波信号、震源坐标、所述地震波信号的传播速度、垂直裂缝坐标的条件下，获取用于模拟所述垂直裂缝对所述地震波信号传播的影响的新的数值模拟表达式，其中，所述新的数值模拟表达式包括：
Figure FDA0002286401510000062
Figure FDA0002286401510000063
Figure FDA0002286401510000071
Figure FDA0002286401510000081
以及
Figure FDA0002286401510000091
其中，t表示采样时刻，(m，n)表示空间网格点在空间网格模型中的坐标位置，Δt表示所述时间采样率，Δx表示所述空间网格模型中的空间网格的长度，Δz表示所述空间网格的宽度，NT表示所述垂直裂缝的切向屈服度、NS表示所述垂直裂缝的法向屈服度，λ、μ表示拉梅系数，ρ表示地质密度，
Figure FDA0002286401510000092
表示所述地震波信号在采样时刻为t+1/2、空间网格点的坐标位置为(m，n)的沿x轴方向的速度，
Figure FDA0002286401510000101
表示所述地震波信号在采样时刻为t-1/2、空间网格点的坐标位置为(m，n)的沿z轴方向的速度，
Figure FDA0002286401510000102
表示所述地震波信号在采样时刻为t+1、空间网格点的坐标位置为(m+1/2，n+1/2)的沿x轴方向的应力，
Figure FDA0002286401510000103
表示所述地震波信号在采样时刻为t、网格点的坐标位置为(m+1/2，n+1/2)的沿z轴方向的应力，
Figure FDA0002286401510000104
表示所述地震波信号在采样时刻为t、网格点的坐标位置为(m+1/2，n+1/2)的沿xz轴方向的应力；C1和C2为有限差分系数；
获取所述地质密度的值、所述拉梅系数的值、所述切向屈服度的值、所述法向屈服度的值、所述采样时刻的值、所述时间采样率的值、所述空间网格模型的大小，以及所述空间网格的长度值和宽度值；
将所述地质密度的值、所述拉梅系数的值、所述切向屈服度的值、所述法向屈服度的值、所述采样时刻的值、所述时间采样率的值、所述空间网格模型中的每个空间网格点的坐标位置，以及所述空间网格的长度值和宽度值分别输入所述新的数值模拟表达式，得到用于模拟所述垂直裂缝对所述地震波信号传播的影响的速度值和应力值。
6.一种基于垂直裂缝介质的地震波波场模拟装置，所述装置包括：
第一模拟表达式获取单元，用于在已知地震波信号、震源坐标、所述地震波信号的传播速度、垂直裂缝坐标的条件下，获取用于模拟所述垂直裂缝对所述地震波信号传播的影响的数值模拟表达式，其中，所述数值模拟表达式包括：
Figure FDA0002286401510000105
Figure FDA0002286401510000111
Figure FDA0002286401510000112
Figure FDA0002286401510000121
以及
Figure FDA0002286401510000122
其中，t表示采样时刻，(m，n)表示空间网格点在空间网格模型中的坐标位置，Δt表示所述时间采样率，Δx表示所述空间网格模型中的空间网格的长度，Δz表示所述空间网格的宽度，NT表示所述垂直裂缝的切向屈服度、NS表示所述垂直裂缝的法向屈服度，λ、μ表示拉梅系数，ρ表示地质密度，
Figure FDA0002286401510000131
表示所述地震波信号在采样时刻为t+1/2、空间网格点的坐标位置为(m，n)的沿x轴方向的速度，
Figure FDA0002286401510000132
表示所述地震波信号在采样时刻为t-1/2、空间网格点的坐标位置为(m，n)的沿z轴方向的速度，
Figure FDA0002286401510000133
表示所述地震波信号在采样时刻为t+1、空间网格点的坐标位置为(m+1/2，n+1/2)的沿x轴方向的应力，
Figure FDA0002286401510000134
表示所述地震波信号在采样时刻为t、网格点的坐标位置为(m+1/2，n+1/2)的沿z轴方向的应力，
Figure FDA0002286401510000135
表示所述地震波信号在采样时刻为t、网格点的坐标位置为(m+1/2，n+1/2)的沿xz轴方向的应力；
第一参数值获取单元，获取所述地质密度的值、所述拉梅系数的值、所述切向屈服度的值、所述法向屈服度的值、所述采样时刻的值、所述时间采样率的值、所述空间网格模型的大小，以及所述空间网格的长度值和宽度值；
数值模拟单元，用于将所述地质密度的值、所述拉梅系数的值、所述切向屈服度的值、所述法向屈服度的值、所述采样时刻的值、所述时间采样率的值、所述空间网格模型中的每个空间网格点的坐标位置，以及所述空间网格的长度值和宽度值分别输入所述数值模拟表达式，得到用于模拟所述垂直裂缝对所述地震波信号传播的影响的速度值和应力值。
7.根据权利要求6所述的装置，所述第一模拟表达式获取单元，包括：
边界条件获取单元，用于基于线性滑动理论，获取表征所述垂直裂缝位置处的边界条件；其中，所述边界条件为：
Figure FDA0002286401510000136
和
Figure FDA0002286401510000141
其中，“+”表示所述垂直裂缝的右侧，“-”表示所述垂直裂缝的左侧；
虚拟网格点处的速度获取单元，用于基于所述边界条件和虚拟网格，获取用于表征所述地震波信号在虚网格点处速度的第一表达式，其中，所述第一表达式包括：
Figure FDA0002286401510000142
其中，
Figure FDA0002286401510000143
为所述地震波信号在虚网格点坐标位置为(m，n)的沿x轴方向的速度；
差分表达式获取单元，用于获取表征应力和所述地震波信号在所述地质中的传播速度之间的关系的有限差分表达式；其中，所述有限差分表达式为：
Figure FDA0002286401510000144
其中，
Figure FDA0002286401510000145
Figure FDA0002286401510000146
为旋转交错网格的第一差分方向；
Figure FDA0002286401510000147
为所述旋转交错网格的第二差分方向；
离散处理单元，用于基于所述虚拟网格和所述旋转交错网格对所述有限差分表达式进行离散处理，得到离散化的差分表达式；其中，所述离散化的差分表达式包括：
Figure FDA0002286401510000151
模拟表达式获取子单元，用于基于所述第一表达式和所述离散化的差分表达式，获取所述数值模拟表达式。
8.一种基于垂直裂缝介质的地震波波场模拟装置，其特征在于，所述装置包括：
第二模拟表达式获取单元，用于在已知地震波信号、震源坐标、所述地震波信号的传播速度、垂直裂缝坐标的条件下，获取用于模拟所述垂直裂缝对所述地震波信号传播的影响的新的数值模拟表达式，其中，所述新的数值模拟表达式包括：
Figure FDA0002286401510000161
Figure FDA0002286401510000162
Figure FDA0002286401510000171
Figure FDA0002286401510000181
以及
Figure FDA0002286401510000191
其中，t表示采样时刻，(m，n)表示空间网格点在空间网格模型中的坐标位置，Δt表示所述时间采样率，Δx表示所述空间网格模型中的空间网格的长度，Δz表示所述空间网格的宽度，NT表示所述垂直裂缝的切向屈服度、NS表示所述垂直裂缝的法向屈服度，λ、μ表示拉梅系数，ρ表示地质密度，
Figure FDA0002286401510000192
表示所述地震波信号在采样时刻为t+1/2、空间网格点的坐标位置为(m，n)的沿x轴方向的速度，
Figure FDA0002286401510000201
表示所述地震波信号在采样时刻为t-1/2、空间网格点的坐标位置为(m，n)的沿z轴方向的速度，
Figure FDA0002286401510000202
表示所述地震波信号在采样时刻为t+1、空间网格点的坐标位置为(m+1/2，n+1/2)的沿x轴方向的应力，
Figure FDA0002286401510000203
表示所述地震波信号在采样时刻为t、网格点的坐标位置为(m+1/2，n+1/2)的沿z轴方向的应力，
Figure FDA0002286401510000204
表示所述地震波信号在采样时刻为t、网格点的坐标位置为(m+1/2，n+1/2)的沿xz轴方向的应力；C1和C2为有限差分系数；
第二参数值获取单元，用于获取所述地质密度的值、所述拉梅系数的值、所述切向屈服度的值、所述法向屈服度的值、所述采样时刻的值、所述时间采样率的值、所述空间网格模型的大小，以及所述空间网格的长度值和宽度值；
第二数值模拟单元，用于将所述地质密度的值、所述拉梅系数的值、所述切向屈服度的值、所述法向屈服度的值、所述采样时刻的值、所述时间采样率的值、所述空间网格模型中的每个空间网格点的坐标位置，以及所述空间网格的长度值和宽度值分别输入所述新的数值模拟表达式，得到用于模拟所述垂直裂缝对所述地震波信号传播的影响的速度值和应力值。
9.一种电子设备，其特征在于，包括存储器以及处理器，所述存储器中存储有计算机程序指令，所述计算机程序指令被所述处理器读取并运行时，执行如权利要求1-5中任一项所述的方法的步骤。
10.一种存储介质，其特征在于，所述存储介质上存储有计算机程序指令，所述计算机程序指令被计算机读取并运行时，执行如权利要求1-5中任一项所述的方法的步骤。
CN201911169482.2A 2019-11-25 2019-11-25 基于垂直裂缝介质的地震波波场模拟方法及装置 Active CN110988988B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201911169482.2A CN110988988B (zh) 2019-11-25 2019-11-25 基于垂直裂缝介质的地震波波场模拟方法及装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911169482.2A CN110988988B (zh) 2019-11-25 2019-11-25 基于垂直裂缝介质的地震波波场模拟方法及装置

Publications (2)

Publication Number Publication Date
CN110988988A true CN110988988A (zh) 2020-04-10
CN110988988B CN110988988B (zh) 2020-10-23

Family

ID=70086697

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911169482.2A Active CN110988988B (zh) 2019-11-25 2019-11-25 基于垂直裂缝介质的地震波波场模拟方法及装置

Country Status (1)

Country Link
CN (1) CN110988988B (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111538079A (zh) * 2020-05-27 2020-08-14 中国矿业大学(北京) 基于全波形反演技术确定地质裂缝柔度参数的方法及装置
CN113569187A (zh) * 2021-07-28 2021-10-29 中国矿业大学(北京) 一种二维水平裂缝的弹性波数值模拟方法及装置
CN113945994A (zh) * 2020-06-30 2022-01-18 中国石油化工股份有限公司 使用有限差分模型进行高速多源加载和波场检索的方法

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20150029822A1 (en) * 2013-07-28 2015-01-29 Robert H. Brune Method and Apparatus for Enhanced Monitoring of Induced Seismicity and Vibration Using Linear Low Frequency and Rotational Sensors
CN104749617A (zh) * 2013-12-26 2015-07-01 中国石油化工股份有限公司 一种多尺度裂缝储层正演模型建立方法
WO2016114896A1 (en) * 2015-01-14 2016-07-21 Chevron U.S.A. Inc. Hydraulic fracture permeability characterization from electromagnetic measurements
CN106597531A (zh) * 2015-10-16 2017-04-26 中国石油化工股份有限公司 含垂直裂缝的页岩的波场传播特征的正演模拟方法
CN108627871A (zh) * 2017-03-24 2018-10-09 中国石油化工股份有限公司 一种tti介质裂缝性质参数的反演方法
CN109444954A (zh) * 2018-12-26 2019-03-08 中国矿业大学(北京) 裂缝数值的模拟方法、装置、电子设备及存储介质
CN109490963A (zh) * 2017-09-13 2019-03-19 中国石油化工股份有限公司 裂缝储层岩石物理建模方法及系统

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20150029822A1 (en) * 2013-07-28 2015-01-29 Robert H. Brune Method and Apparatus for Enhanced Monitoring of Induced Seismicity and Vibration Using Linear Low Frequency and Rotational Sensors
CN104749617A (zh) * 2013-12-26 2015-07-01 中国石油化工股份有限公司 一种多尺度裂缝储层正演模型建立方法
WO2016114896A1 (en) * 2015-01-14 2016-07-21 Chevron U.S.A. Inc. Hydraulic fracture permeability characterization from electromagnetic measurements
CN106597531A (zh) * 2015-10-16 2017-04-26 中国石油化工股份有限公司 含垂直裂缝的页岩的波场传播特征的正演模拟方法
CN108627871A (zh) * 2017-03-24 2018-10-09 中国石油化工股份有限公司 一种tti介质裂缝性质参数的反演方法
CN109490963A (zh) * 2017-09-13 2019-03-19 中国石油化工股份有限公司 裂缝储层岩石物理建模方法及系统
CN109444954A (zh) * 2018-12-26 2019-03-08 中国矿业大学(北京) 裂缝数值的模拟方法、装置、电子设备及存储介质

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
全红娟 等: "地震波横波分裂现象的数值模拟特征", 《西北大学学报》 *
刘财 等: "裂缝诱导双相ＨＴ Ｉ介质地震波场错格伪谱法模拟与波场特征分析", 《地球物理学报》 *

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111538079A (zh) * 2020-05-27 2020-08-14 中国矿业大学(北京) 基于全波形反演技术确定地质裂缝柔度参数的方法及装置
CN111538079B (zh) * 2020-05-27 2020-12-25 中国矿业大学(北京) 基于全波形反演技术确定地质裂缝柔度参数的方法及装置
CN113945994A (zh) * 2020-06-30 2022-01-18 中国石油化工股份有限公司 使用有限差分模型进行高速多源加载和波场检索的方法
CN113945994B (zh) * 2020-06-30 2023-07-04 中国石油化工股份有限公司 使用有限差分模型进行高速多源加载和波场检索的方法
CN113569187A (zh) * 2021-07-28 2021-10-29 中国矿业大学(北京) 一种二维水平裂缝的弹性波数值模拟方法及装置
CN113569187B (zh) * 2021-07-28 2022-03-29 中国矿业大学(北京) 一种二维水平裂缝的弹性波数值模拟方法及装置

Also Published As

Publication number Publication date
CN110988988B (zh) 2020-10-23

Similar Documents

Publication Publication Date Title
CN110988988B (zh) 基于垂直裂缝介质的地震波波场模拟方法及装置
Tarrass et al. New curvilinear scheme for elastic wave propagation in presence of curved topography
US9341729B2 (en) Amplitude contrast seismic attribute
Ladopoulos Elastodynamics for Non-linear Seismic Wave Motion in Real-Time Expert Seismology
Marsden et al. A study of infrasound propagation based on high-order finite difference solutions of the Navier-Stokes equations
CN107765308B (zh) 基于褶积思想与精确震源的重构低频数据频域全波形反演方法
BRPI1101394B1 (pt) método para separar campos de pressão e velocidade vertical medidos em tiras rebocadas em campos de ondas ascendentes e descendentes e meio legível por computador
CN107894618B (zh) 一种基于模型平滑算法的全波形反演梯度预处理方法
EP2497043A1 (en) Seismic imaging systems and methods employing a 3d reverse time migration with tilted transverse isotropy
CN111538079B (zh) 基于全波形反演技术确定地质裂缝柔度参数的方法及装置
Sharifi et al. Novel permeability upscaling method using fast marching method
Guo et al. Simulation of deep-water waves based on JONSWAP spectrum and realization by MATLAB
WO2012128873A2 (en) System and method for seismic imaging with reduced computational cost
CN105425298A (zh) 一种消除有限差分正演过程中数值频散的方法和装置
CN109709602B (zh) 一种远探测声波偏移成像方法、装置及系统
CA2953100A1 (en) System and method for rock property estimation of subsurface geologic volumes
CN104391325A (zh) 不连续非均质地质体检测方法和装置
Hjeij et al. Comparing advanced discretization methods for complex hydrocarbon reservoirs
US20170206288A1 (en) Method and apparatus for constructing and using absorbing boundary conditions in numberical computations of physical applications
CN113311484B (zh) 利用全波形反演获取粘弹性介质弹性参数的方法及装置
CN113221409B (zh) 一种有限元和边界元耦合的声波二维数值模拟方法和装置
Chambers et al. A practical implementation of wave front construction for 3-D isotropic media
Cai et al. Numerical manifold method with local pixel representation of finite covers for two-dimensional problems having complex discontinuities
Goldthorpe et al. Unconventional modelling of faulted reservoirs: a case study
CN114063161A (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