CN108873063B - 一种微地震矩张量反演的方法及装置 - Google Patents
一种微地震矩张量反演的方法及装置 Download PDFInfo
- Publication number
- CN108873063B CN108873063B CN201710343149.3A CN201710343149A CN108873063B CN 108873063 B CN108873063 B CN 108873063B CN 201710343149 A CN201710343149 A CN 201710343149A CN 108873063 B CN108873063 B CN 108873063B
- Authority
- CN
- China
- Prior art keywords
- waveform
- frequency
- residual
- simulated
- frequency domain
- 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
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. analysis, for interpretation, for correction
Abstract
本发明实施例提供了一种微地震矩张量反演的方法及装置，其中，该方法包括：获取频率域的观测微地震波形，在观测微地震波形中选取信噪比大于预设值的波形频段，观测微地震波形是通过传感器记录得到的；获取频率域的模拟微地震波形，模拟微地震波形是根据格林函数计算得到的；计算所选的波形频段的振幅谱和模拟微地震波形的振幅谱的频率域振幅谱残差，计算所选的波形频段的相位谱和模拟微地震波形的相位谱的频率域相位谱残差，计算频率域振幅谱残差和频率域相位谱残差的联合残差；根据联合残差最小值对应的模拟微地震波形的相位谱和模拟微地震波形的振幅谱，反演微地震矩张量。本申请提高了反演结果的可靠性。
Description
技术领域
本发明涉及微地震监测技术领域，特别涉及一种微地震矩张量反演的方法及装置。
背景技术
微地震监测是了解井下水力压裂效果的重要手段，震源机制的矩张量反演是了解压裂过程中裂缝破裂情况的主要工具。采用广义反透射系数方法正演理论地震记录，用矩张量描述震源属性，求解观测记录和矩张量的线性方程组，反演出震源机制解，并得到水力压裂所产生裂缝的方位和倾角等参数，为压裂裂缝的解释和压裂效果评价提供依据。
震源参数包括震源断层面的走向、倾向和倾角，震源断层两盘错动的方向、幅度，震源断层面的长度、宽度、断层破裂的扩展速度、震源主应力状态以及错动时释放的应力等。震源机制和震源参数的资料对微地震机理分析至关重要。
矩张量反演现有技术中有许多种方法，通常采用在时域里进行矩张量反演，要求获取高质量的三分量记录，但是由于噪音、不精确的速度模型和不合理的检波器布设等因素的影响，大量微地震虽然能被检波器记录到，却由于信噪比低不够计算地震的震源机制解，使得反演结果可靠性低。
发明内容
本发明实施例提供了一种微地震矩张量反演的方法，以解决现有技术中在时域里进行矩张量反演时存在的反演结果可靠性低的技术问题。该方法包括：获取频率域的观测微地震波形，在所述观测微地震波形中选取信噪比大于预设值的波形频段，所述观测微地震波形是通过传感器记录得到的；获取频率域的模拟微地震波形，所述模拟微地震波形是根据格林函数模拟得到的；计算所选的波形频段的振幅谱和所述模拟微地震波形的振幅谱的频率域振幅谱残差，计算所选的波形频段的相位谱和所述模拟微地震波形的相位谱的频率域相位谱残差，计算所述频率域振幅谱残差和所述频率域相位谱残差的联合残差，并计算所述联合残差的最小值；根据所述联合残差最小值对应的所述模拟微地震波形的相位谱和所述模拟微地震波形的振幅谱，反演地震矩张量。
在一个实施例中，所述预设值为3。
在一个实施例中，还包括：通过快速傅立叶变换获取所述频率域的观测微地震波形和所述频率域的模拟微地震波形。
在一个实施例中，通过以下公式计算所述频率域振幅谱残差和所述频率域相位谱残差的联合残差：ε(f)＝1-(1-εph(f))(1-εamp(f))；其中，ε(f)是联合残差，εph(f)是所述频率域相位谱残差；是所述模拟微地震波形的频率f在x点的相位；是所选的波形频段的频率f在x点的相位；εamp(f)是所述频率域振幅谱残差；是所述模拟微地震波形的频率f在x点的振幅；是所选的波形频段的频率f在x点的振幅。
在一个实施例中，通过以下公式根据所述联合残差最小值对应的所述模拟微地震波形的相位谱和所述模拟微地震波形的振幅谱，反演微地震矩张量：其中，u(x,f)是频率域里频率f在点x所产生的位移；Mij(f)是频率f的微地震矩张量；Gi,j(f)是频率f包括震源至接收点间的传播效应的格林函数；i是微地震矩张量矩阵的行；j是微地震矩张量矩阵的列。
本发明实施例还提供了一种微地震矩张量反演的装置，以解决现有技术中在时域里进行矩张量反演时存在的反演结果可靠性低的技术问题。该装置包括：第一波形获取模块，用于获取频率域的观测微地震波形，在所述观测微地震波形中选取信噪比大于预设值的波形频段，所述观测微地震波形是通过传感器记录得到的；第二波形获取模块，用于获取频率域的模拟微地震波形，所述模拟微地震波形是根据格林函数模拟得到的；计算模块，用于计算所选的波形频段的振幅谱和所述模拟微地震波形的振幅谱的频率域振幅谱残差，计算所选的波形频段的相位谱和所述模拟微地震波形的相位谱的频率域相位谱残差，计算所述频率域振幅谱残差和所述频率域相位谱残差的联合残差，并计算所述联合残差的最小值；反演模块，用于根据所述联合残差最小值对应的所述模拟微地震波形的相位谱和所述模拟微地震波形的振幅谱，反演微地震矩张量。
在一个实施例中，所述预设值为3。
在一个实施例中，还包括：波形变换模块，用于通过快速傅立叶变换获取所述频率域的观测微地震波形和所述频率域的模拟微地震波形。
在一个实施例中，所述计算模块具体通过以下公式计算所述频率域振幅谱残差和所述频率域相位谱残差的联合残差：ε(f)＝1-(1-εph(f))(1-εamp(f))；其中，ε(f)是联合残差，εph(f)是所述频率域相位谱残差；是所述模拟微地震波形的频率f在x点的相位；是所选的波形频段的频率f在x点的相位；εamp(f)是所述频率域振幅谱残差；是所述模拟微地震波形的频率f在x点的振幅；是所选的波形频段的频率f在x点的振幅。
在一个实施例中，所述反演模块具体通过以下公式根据所述联合残差最小值对应的所述模拟微地震波形的相位谱和所述模拟微地震波形的振幅谱，反演微地震矩张量：其中，u(x,f)是频率域里频率f在点x所产生的位移；Mij(f)是频率f的微地震矩张量；Gi,j(f)是频率f包括震源至接收点间的传播效应的格林函数；i是微地震矩张量矩阵的行；j是微地震矩张量矩阵的列。
在本发明实施例中，提出采用频率域的观测微地震波形和模拟微地震波形进行微地震矩张量反演，具体采用频率域振幅谱残差和频率域相位谱残差的联合残差进行反演；同时，在反演前在观测微地震波形中选取信噪比大于预设值的波形频段，以便基于选取的信噪比高的波形频段进行反演，摒弃信噪比低的频带。使得反演结果的稳健性对噪音干扰不敏感，提高了反演结果的可靠性。
附图说明
此处所说明的附图用来提供对本发明的进一步理解，构成本申请的一部分，并不构成对本发明的限定。在附图中：
图1是本发明实施例提供的一种微地震矩张量反演的方法的流程图；
图2是本发明实施例提供的一种二阶矩张量的分量示意图；
图3是本发明实施例提供的一种具体的微地震矩张量反演方法的流程图；
图4是本发明实施例提供的一种微地震矩张量反演的装置的结构框图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚明白，下面结合实施方式和附图，对本发明做进一步详细说明。在此，本发明的示意性实施方式及其说明用于解释本发明，但并不作为对本发明的限定。
在本发明实施例中，提供了一种微地震矩张量反演的方法，如图1所示，该方法包括：
步骤101：获取频率域的观测微地震波形，在所述观测微地震波形中选取信噪比大于预设值的波形频段，所述观测微地震波形是通过传感器记录得到的；
步骤102：获取频率域的模拟微地震波形，所述模拟微地震波形是根据格林函数计算得到的；
步骤103：计算所选的波形频段的振幅谱和所述模拟微地震波形的振幅谱的频率域振幅谱残差，计算所选的波形频段的相位谱和所述模拟微地震波形的相位谱的频率域相位谱残差，计算所述频率域振幅谱残差和所述频率域相位谱残差的联合残差，并计算所述联合残差的最小值；
步骤104：根据所述联合残差最小值对应的所述模拟微地震波形的相位谱和所述模拟微地震波形的振幅谱，反演微地震矩张量。
由图1所示的流程可知，在本发明实施例中，提出采用频率域的观测微地震波形和模拟微地震波形进行微地震矩张量反演，具体采用频率域振幅谱残差和频率域相位谱残差的联合残差进行反演；同时，在反演前在观测微地震波形中选取信噪比大于预设值的波形频段，以便基于选取的信噪比高的波形频段进行反演，摒弃信噪比低的频带。使得反演结果的稳健性对噪音干扰不敏感，提高了反演结果的可靠性。
具体实施时，在所述观测微地震波形中选取信噪比大于预设值的波形频段时，所述预设值可以为3。即基于信噪比大于3的波形频段进行微地震矩张量反演。
具体实施时，上述微地震矩张量反演的方法，还包括：通过快速傅立叶变换获取所述频率域的观测微地震波形和所述频率域的模拟微地震波形。
具体的，可以通过传感器观测记录时域的观测微地震波形，然后，通过快速傅立叶变换将时域的观测微地震波形变换为频率域的观测微地震波形。同时，可以先根据震源至接收点间的传播效应的格林函数模拟计算时域的模拟微地震波形，然后采用快速傅立叶变换将时域的模拟微地震波形变换为频率域的模拟微地震波形。
具体实施时，可以通过以下公式计算所述频率域振幅谱残差和所述频率域相位谱残差的联合残差：
ε(f)＝1-(1-εph(f))(1-εamp(f)) (1)
其中，
ε(f)是联合残差，εph(f)是所述频率域相位谱残差；是所述模拟微地震波形的频率f在x点的相位；是所选的波形频段的频率f在x点的相位；εamp(f)是所述频率域振幅谱残差；是所述模拟微地震波形的频率f在x点的振幅；是所选的波形频段的频率f在x点的振幅。
上述联合残差可采用L1范型或L2范型均可。
具体实施时，在计算上述联合残差时通过对上述公式(1)求偏导数来计算联合残差的最小值，计算得到联合残差的最小值后，获得联合残差最小值对应的所选的波形频段的相位谱、所选的波形频段的振幅谱、所述模拟微地震波形的相位谱以及所述模拟微地震波形的振幅谱，进而通过以下公式根据所述联合残差最小值对应的所选的波形频段的相位谱、所选的波形频段的振幅谱、所述模拟微地震波形的相位谱以及所述模拟微地震波形的振幅谱，反演微地震矩张量：
其中，u(x,f)是频率域里频率f在点x所产生的位移；Mij(f)是频率f的微地震矩张量；Gi,j(f)是频率f包括震源至接收点间的传播效应的格林函数；i是微地震矩张量矩阵的行；j是微地震矩张量矩阵的列。
上述公式(4)可以采用最小二乘法进行求解。
具体的，
mij为常数，代表二阶矩张量Mij的分量。若i＝j，表明力和力臂在同一方向，为无矩单力偶；若i≠j，表明力作用于i方向，力臂在j方向，为一个力矩为mij的单力偶，具体的，如图2所示。
现有技术中，则采用以下公式(5)在时域进行微地震矩张量反演：
u(x,t)＝Mij[Gi,j*s(t)] (5)
u(x,t)为时间域里时间t在点x所产生的位移；Mij为微地震矩张量；Gi,j为包括震源至接收点间的传播效应的格林函数；s(t)为震源的时间函数。
以下具体描述上述微地震矩张量反演的方法的工作流程，如图3所示，该方法包括：
1、获取微地震波的原始资料，即获取传感器观测得到的时域的观测微地震波形，同时，根据格林函数模拟时域的模拟微地震波形，对时域的观测微地震波形和时域的模拟微地震波形进行预处理；
2、通过快速傅立叶变换将时域的观测微地震波形变换为频率域的观测微地震波形，将时域的模拟微地震波形变换为频率域的模拟微地震波形；
3、在频率域的观测微地震波形中选取信噪比高(例如，信噪比大于3)的频带；
4、计算频率域的模拟微地震波形的振幅谱和相位谱，计算所选频段波形的振幅谱和相位谱，计算频率域的模拟微地震波形的振幅谱和所选频段波形的振幅谱的频率域振幅谱残差，计算模拟微地震波形的相位谱和所选频段波形的相位谱的频率域相位谱残差，计算频率域振幅谱残差和频率域相位谱残差的联合残差的最小值；可以采用公式(1)、(2)以及(3)；
5、根据所述联合残差最小值对应的所述模拟微地震波形的相位谱以及所述模拟微地震波形的振幅谱，反演微地震矩张量，具体可以采用公式(4)，直至小于某一限定的法制，计算结束。
上述微地震矩张量反演的方法的关键在于，根据波形记录找到频率域信噪比高的频带进行计算，摈弃信噪比低的频带，从而提高反演可靠性。
基于同一发明构思，本发明实施例中还提供了一种微地震矩张量反演的装置，如下面的实施例所述。由于微地震矩张量反演的装置解决问题的原理与微地震矩张量反演的方法相似，因此微地震矩张量反演的装置的实施可以参见微地震矩张量反演的方法的实施，重复之处不再赘述。以下所使用的，术语“单元”或者“模块”可以实现预定功能的软件和/或硬件的组合。尽管以下实施例所描述的装置较佳地以软件来实现，但是硬件，或者软件和硬件的组合的实现也是可能并被构想的。
图4是本发明实施例的微地震矩张量反演的装置的一种结构框图，如图4所示，该装置包括：
第一波形获取模块401，用于获取频率域的观测微地震波形，在所述观测微地震波形中选取信噪比大于预设值的波形频段，所述观测微地震波形是通过传感器记录得到的；
第二波形获取模块402，用于获取频率域的模拟微地震波形，所述模拟微地震波形是根据格林函数计算得到的；
计算模块403，用于计算所选的波形频段的振幅谱和所述模拟微地震波形的振幅谱的频率域振幅谱残差，计算所选的波形频段的相位谱和所述模拟微地震波形的相位谱的频率域相位谱残差，计算所述频率域振幅谱残差和所述频率域相位谱残差的联合残差，并计算所述联合残差的最小值；
反演模块404，用于根据所述联合残差最小值对应的所述模拟微地震波形的相位谱和所述模拟微地震波形的振幅谱，反演微地震矩张量。
在一个实施例中，所述预设值为3。
在一个实施例中，还包括：波形变换模块，用于通过快速傅立叶变换获取所述频率域的观测微地震波形和所述频率域的模拟微地震波形。
在一个实施例中，所述计算模块具体通过以下公式计算所述频率域振幅谱残差和所述频率域相位谱残差的联合残差：ε(f)＝1-(1-εph(f))(1-εamp(f))；其中，ε(f)是联合残差，εph(f)是所述频率域相位谱残差；是所述模拟微地震波形的频率f在x点的相位；是所选的波形频段的频率f在x点的相位；εamp(f)是所述频率域振幅谱残差；是所述模拟微地震波形的频率f在x点的振幅；是所选的波形频段的频率f在x点的振幅。
在一个实施例中，所述反演模块具体通过以下公式根据所述联合残差最小值对应的所述模拟微地震波形的相位谱和所述模拟微地震波形的振幅谱，反演微地震矩张量：其中，u(x,f)是频率域里频率f在点x所产生的位移；Mij(f)是频率f的微地震矩张量；Gi,j(f)是频率f包括震源至接收点间的传播效应的格林函数；i是微地震矩张量矩阵的行；j是微地震矩张量矩阵的列。
在本发明实施例中，提出采用频率域的观测微地震波形和模拟微地震波形进行微地震矩张量反演，具体采用频率域振幅谱残差和频率域相位谱残差的联合残差进行反演；同时，在反演前在观测微地震波形中选取信噪比大于预设值的波形频段，以便基于选取的信噪比高的波形频段进行反演，摒弃信噪比低的频带。使得反演结果的稳健性对噪音干扰不敏感，提高了反演结果的可靠性。
显然，本领域的技术人员应该明白，上述的本发明实施例的各模块或各步骤可以用通用的计算装置来实现，它们可以集中在单个的计算装置上，或者分布在多个计算装置所组成的网络上，可选地，它们可以用计算装置可执行的程序代码来实现，从而，可以将它们存储在存储装置中由计算装置来执行，并且在某些情况下，可以以不同于此处的顺序执行所示出或描述的步骤，或者将它们分别制作成各个集成电路模块，或者将它们中的多个模块或步骤制作成单个集成电路模块来实现。这样，本发明实施例不限制于任何特定的硬件和软件结合。
以上所述仅为本发明的优选实施例而已，并不用于限制本发明，对于本领域的技术人员来说，本发明实施例可以有各种更改和变化。凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。
Claims (8)
1.一种微地震矩张量反演的方法，其特征在于，包括：
获取频率域的观测微地震波形，在所述观测微地震波形中选取信噪比大于预设值的波形频段，所述观测微地震波形是通过传感器记录得到的；
获取频率域的模拟微地震波形，所述模拟微地震波形是根据格林函数计算得到的；
计算所选的波形频段的振幅谱和所述模拟微地震波形的振幅谱的频率域振幅谱残差，计算所选的波形频段的相位谱和所述模拟微地震波形的相位谱的频率域相位谱残差，计算所述频率域振幅谱残差和所述频率域相位谱残差的联合残差，并计算所述联合残差的最小值；
根据所述联合残差最小值对应的所述模拟微地震波形的相位谱和所述模拟微地震波形的振幅谱，反演微地震矩张量；
通过以下公式根据所述联合残差最小值对应的所述模拟微地震波形的相位谱和所述模拟微地震波形的振幅谱，反演微地震矩张量：
其中，u(x,f)是频率域里频率f在点x所产生的位移；Mij(f)是频率f的微地震矩张量；Gi,j(f)是频率f包括震源至接收点间的传播效应的格林函数；i是微地震矩张量矩阵的行；j是微地震矩张量矩阵的列；是所述模拟微地震波形的频率f在x点的振幅；是所述模拟微地震波形的频率f在x点的相位。
2.如权利要求1所述的微地震矩张量反演的方法，其特征在于，所述预设值为3。
3.如权利要求1所述的微地震矩张量反演的方法，其特征在于，还包括：
通过快速傅立叶变换获取所述频率域的观测微地震波形和所述频率域的模拟微地震波形。
4.如权利要求1至3中任一项所述的微地震矩张量反演的方法，其特征在于，通过以下公式计算所述频率域振幅谱残差和所述频率域相位谱残差的联合残差：
ε(f)＝1-(1-εph(f))(1-εamp(f))；
其中，ε(f)是联合残差，εph(f)是所述频率域相位谱残差；是所述模拟微地震波形的频率f在x点的相位；是所选的波形频段的频率f在x点的相位；εamp(f)是所述频率域振幅谱残差；是所述模拟微地震波形的频率f在x点的振幅；是所选的波形频段的频率f在x点的振幅。
5.一种微地震矩张量反演的装置，其特征在于，包括：
第一波形获取模块，用于获取频率域的观测微地震波形，在所述观测微地震波形中选取信噪比大于预设值的波形频段，所述观测微地震波形是通过传感器记录得到的；
第二波形获取模块，用于获取频率域的模拟微地震波形，所述模拟微地震波形是根据格林函数计算得到的；
计算模块，用于计算所选的波形频段的振幅谱和所述模拟微地震波形的振幅谱的频率域振幅谱残差，计算所选的波形频段的相位谱和所述模拟微地震波形的相位谱的频率域相位谱残差，计算所述频率域振幅谱残差和所述频率域相位谱残差的联合残差，并计算所述联合残差的最小值；
反演模块，用于根据所述联合残差最小值对应的所述模拟微地震波形的相位谱和所述模拟微地震波形的振幅谱，反演微地震矩张量；
所述反演模块具体通过以下公式根据所述联合残差最小值对应的所述模拟微地震波形的相位谱和所述模拟微地震波形的振幅谱，反演微地震矩张量：
其中，u(x,f)是频率域里频率f在点x所产生的位移；Mij(f)是频率f的微地震矩张量；Gi,j(f)是频率f包括震源至接收点间的传播效应的格林函数；i是微地震矩张量矩阵的行；j是微地震矩张量矩阵的列；是所述模拟微地震波形的频率f在x点的振幅；是所述模拟微地震波形的频率f在x点的相位。
6.如权利要求5所述的微地震矩张量反演的装置，其特征在于，所述预设值为3。
7.如权利要求5所述的微地震矩张量反演的装置，其特征在于，还包括：
波形变换模块，用于通过快速傅立叶变换获取所述频率域的观测微地震波形和所述频率域的模拟微地震波形。
8.如权利要求5至7中任一项所述的微地震矩张量反演的装置，其特征在于，所述计算模块具体通过以下公式计算所述频率域振幅谱残差和所述频率域相位谱残差的联合残差：
ε(f)＝1-(1-εph(f))(1-εamp(f))；
其中，ε(f)是联合残差，εph(f)是所述频率域相位谱残差；是所述模拟微地震波形的频率f在x点的相位；是所选的波形频段的频率f在x点的相位；εamp(f)是所述频率域振幅谱残差；是所述模拟微地震波形的频率f在x点的振幅；是所选的波形频段的频率f在x点的振幅。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710343149.3A CN108873063B (zh) | 2017-05-16 | 2017-05-16 | 一种微地震矩张量反演的方法及装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710343149.3A CN108873063B (zh) | 2017-05-16 | 2017-05-16 | 一种微地震矩张量反演的方法及装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108873063A CN108873063A (zh) | 2018-11-23 |
| CN108873063B true CN108873063B (zh) | 2019-12-13 |
Family
ID=64320699
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710343149.3A Active CN108873063B (zh) | 2017-05-16 | 2017-05-16 | 一种微地震矩张量反演的方法及装置 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108873063B (zh) |
Families Citing this family (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110907538A (zh) * | 2019-11-29 | 2020-03-24 | 神华准格尔能源有限责任公司 | 岩石声发射的混合矩张量反演计算方法、存储介质和系统 |
| CN111175815B (zh) * | 2020-01-06 | 2022-04-15 | 中国石油化工股份有限公司 | 油藏改造微地震监测裂缝震源机制求解方法及系统 |
| CN111462106A (zh) * | 2020-04-09 | 2020-07-28 | 中山易美杰智能科技有限公司 | 一种深度学习图像识别输入所用张量的生成方法及其应用 |
| CN112748465B (zh) * | 2020-12-30 | 2021-12-10 | 中国矿业大学(北京) | 基于岩石特征的震源机制反演方法及装置 |
| CN113050158B (zh) * | 2021-03-19 | 2022-05-31 | 中国科学院武汉岩土力学研究所 | 近场微震信号波形的分析方法、装置、设备及存储介质 |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US10302786B2 (en) * | 2013-12-30 | 2019-05-28 | Cgg Services Sas | Methods and systems of determining a fault plane of a microseismic event |
| CN105467444B (zh) * | 2015-12-10 | 2017-11-21 | 中国石油天然气集团公司 | 一种弹性波全波形反演方法及装置 |
-
2017
- 2017-05-16 CN CN201710343149.3A patent/CN108873063B/zh active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN108873063A (zh) | 2018-11-23 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108873063B (zh) | 一种微地震矩张量反演的方法及装置 | |
| Pioldi et al. | Output-only modal dynamic identification of frames by a refined FDD algorithm at seismic input and high damping | |
| Ji et al. | Source description of the 1999 Hector Mine, California, earthquake, part I: Wavelet domain inversion theory and resolution analysis | |
| Boore et al. | Processing of strong-motion accelerograms: needs, options and consequences | |
| CN111221037B (zh) | 解耦弹性逆时偏移成像方法和装置 | |
| RU2010103987A (ru) | Способы и системы для обработки микросейсмических данных | |
| KR20110057124A (ko) | 지진 표면파들의 파형들을 사용하는 토양 특성들의 추정 | |
| CN105588883B (zh) | 三维岩石力学参数获取方法和系统 | |
| US20150355356A1 (en) | Method for Temporal Dispersion Correction for Seismic Simulation, RTM and FWI | |
| Pioldi et al. | Earthquake‐induced structural response output‐only identification by two different Operational Modal Analysis techniques | |
| CN103852785B (zh) | 地层各向异性的评价方法 | |
| WO2014084751A1 (en) | A method for processing acoustic waveforms | |
| US20070064531A1 (en) | Non-linear seismic trace matching to well logs | |
| Zhu et al. | Removing mass loading effects of multi-transducers using Sherman-Morrison-Woodbury formula in modal test | |
| CN113341455A (zh) | 一种粘滞各向异性介质地震波数值模拟方法、装置及设备 | |
| KR101820850B1 (ko) | 직접 파형 역산의 반복 적용을 이용한 탄성파 영상화 장치 및 방법 | |
| Ghaderi et al. | Practical high-fidelity frequency-domain force and location identification | |
| CN106574980A (zh) | 用于地下地质体的岩石性质估计的系统和方法 | |
| CN114114403A (zh) | 一种基于分数阶拉氏算子的各向异性衰减介质模拟方法 | |
| Beresnev | Simulation of near-fault high-frequency ground motions from the representation theorem | |
| Aimar et al. | Novel techniques for in-situ estimation of shear-wave velocity and damping ratio through MASW testing Part I: A beamforming procedure for extracting Rayleigh-wave phase velocity and phase attenuation | |
| Yablokov et al. | Uncertainty quantification of multimodal surface wave inversion using artificial neural networks | |
| CN109425892B (zh) | 地震子波的估计方法及系统 | |
| Ramos-Martinez et al. | Source-parameter estimation by full waveform inversion in 3D heterogeneous, viscoelastic, anisotropic media | |
| CN112649848B (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 |