CN110907538A - 岩石声发射的混合矩张量反演计算方法、存储介质和系统 - Google Patents

岩石声发射的混合矩张量反演计算方法、存储介质和系统 Download PDF

Info

Publication number
CN110907538A
CN110907538A CN201911205651.3A CN201911205651A CN110907538A CN 110907538 A CN110907538 A CN 110907538A CN 201911205651 A CN201911205651 A CN 201911205651A CN 110907538 A CN110907538 A CN 110907538A
Authority
CN
China
Prior art keywords
seismic source
source mechanism
moment tensor
mechanism solution
acoustic emission
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.)
Pending
Application number
CN201911205651.3A
Other languages
English (en)
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.)
Shenhua Zhungeer Energy Co Ltd
Original Assignee
Shenhua Zhungeer Energy Co Ltd
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 Shenhua Zhungeer Energy Co Ltd filed Critical Shenhua Zhungeer Energy Co Ltd
Priority to CN201911205651.3A priority Critical patent/CN110907538A/zh
Publication of CN110907538A publication Critical patent/CN110907538A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N29/00—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
    • G01N29/14—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object using acoustic emission techniques
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N3/00—Investigating strength properties of solid materials by application of mechanical stress
    • G01N3/02—Details
    • G01N3/06—Special adaptations of indicating or recording means
    • 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/15—Correlation function computation including computation of convolution operations
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2203/00—Investigating strength properties of solid materials by application of mechanical stress
    • G01N2203/0058—Kind of property studied
    • G01N2203/006—Crack, flaws, fracture or rupture
    • G01N2203/0062—Crack or flaws
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2203/00—Investigating strength properties of solid materials by application of mechanical stress
    • G01N2203/02—Details not specific for a particular testing method
    • G01N2203/06—Indicating or recording means; Sensing means
    • G01N2203/0658—Indicating or recording means; Sensing means using acoustic or ultrasonic detectors
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2291/00—Indexing codes associated with group G01N29/00
    • G01N2291/02—Indexing codes associated with the analysed material
    • G01N2291/028—Material parameters
    • G01N2291/0289—Internal structure, e.g. defects, grain size, texture

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • Chemical & Material Sciences (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Theoretical Computer Science (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Health & Medical Sciences (AREA)
  • Mathematical Physics (AREA)
  • Analytical Chemistry (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Data Mining & Analysis (AREA)
  • Algebra (AREA)
  • Computing Systems (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Acoustics & Sound (AREA)
  • Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)

Abstract

本发明提出一种岩石声发射的混合矩张量反演计算方法、存储介质和系统，方法包括利用最小二乘法进行绝对矩张量反演计算，得到修正前的震源机制解；引入迭代参数，根据迭代参数计算权重修正值；通过权重修正值对修正前的震源机制解进行修正，得到修正后的震源机制解。本发明解决了由于岩石的的非均质性导致的P波传播分散的技术问题，为岩石内部微破裂的震源机制解提出一种较好的计算方法。

Description

岩石声发射的混合矩张量反演计算方法、存储介质和系统
技术领域
本发明涉及计算机领域，尤其涉及一种岩石声发射的混合矩张量反演计算方法、存储介质和系统。
背景技术
岩石材料是一种各向异性的脆性材料，一般在受到外部荷载的作用下岩石内部会发生脆性破坏，产生微破裂。随着荷载的改变，内部的微裂隙逐渐的萌生扩展，直至完全贯通。因此，通过声发射实验对岩石加载过程中产生的声发射信号进行监测对分析岩石的破坏状态具有重要意义。
一般来说，岩石微破裂引起的声发射信息可以用于分析岩石微破裂的性质(声发射震源机制解)、数量(声发射撞击数)、空间位置(声发射信号定位)。目前，针对声发射撞击数的监测以及微破裂空间位置的定位技术已经非常成熟了，但针对震源机制解的矩张量反演技术依旧不成熟。
发明内容
基于以上问题，本发明提出一种岩石声发射的混合矩张量反演计算方法、存储介质和系统，解决了由于岩石的的非均质性导致的P波传播分散的技术问题，为岩石内部微破裂的震源机制解提出一种较好的计算方法。
本发明提出一种岩石声发射的混合矩张量反演计算方法，包括：
利用最小二乘法进行绝对矩张量反演计算，得到修正前的震源机制解；
引入迭代参数，根据迭代参数计算权重修正值；
通过权重修正值对修正前的震源机制解进行修正，得到修正后的震源机制解。
此外，利用最小二乘法进行绝对矩张量反演计算，得到修正前的震源机制解的公式为：
Mi＝Gi -1ui obs，其中Gi为格林函数，ui为接收到的弹性波的波形位移振幅，obs代表观测值，Mi为震源机制解。
此外，引入迭代参数，根据迭代参数计算权重修正值至少包括：
引入第一迭代参数和第二迭代参数，第一迭代参数的计算公式为：
Figure BDA0002296873550000021
其中
Figure BDA0002296873550000022
为第一迭代系数，median为取中位值，uijkl th为弹性波的波形位移振幅理论值，uijkl obs为弹性波的波形位移振幅实际值，Neq是指所需计算的方程总数。
此外，计算第二迭代参数的公式为：
Figure BDA0002296873550000023
其中，ωIterNo为第二迭代参数，IterNo为自然数。
此外，根据迭代参数计算权重修正值的公式为：
Figure BDA0002296873550000024
其中，ajkl为权重修正值。
此外，通过权重修正值对修正前的震源机制解进行修正，得到修正后的震源机制解包括：先计算修正后的弹性波的波形位移振幅实际值，再计算修正后的震源机制解；
修正后的弹性波的波形位移振幅实际值的计算公式为：
uijkl new＝uijkl old+ajkl·uijkl old
其中，uijkl old为弹性波的波形位移振幅实际值，uijkl new为修正后的弹性波的波形位移振幅实际值。
此外，计算修正后的震源机制解的计算公式为：
Mi new＝Gi -1ui new其中，ui new为修正后的弹性波的波形位移振幅实际值，Mi new为修正后的震源机制解。
此外，利用公式ui th＝GiMi重新计算弹性波的波形位移振幅理论值，重复引入迭代参数，根据迭代参数计算权重修正值，通过权重修正值对修正前的震源机制解进行修正，得到修正后的震源机制解的步骤，通过多次修正得到多次修正后的震源机制解。
本发明提出一种存储介质，所述存储介质存储计算机指令，当计算机执行所述计算机指令时，用于执行上述任一项所述的岩石声发射的混合矩张量反演计算方法。
本发明提出一种岩石声发射的混合矩张量反演计算系统，包括至少一个处理器；以及，与所述至少一个处理器通信连接的存储器；其中，所述存储器存储有可被所述一个处理器执行的指令，所述指令被所述至少一个处理器执行，以使所述至少一个处理器能够：
利用最小二乘法进行绝对矩张量反演计算，得到修正前的震源机制解；
引入迭代参数，根据迭代参数计算权重修正值；
通过权重修正值对修正前的震源机制解进行修正，得到修正后的震源机制解。
通过采用上述技术方案，具有如下有益效果：
本发明解决了由于岩石的的非均质性导致的P波传播分散的技术问题，为岩石内部微破裂的震源机制解提出一种较好的计算方法。
附图说明
图1是本发明一个实施例提供的岩石声发射的混合矩张量反演计算方法的流程图；
图2是本发明一个实施例提供的岩石声发射的混合矩张量反演计算方法中的岩石三点弯曲声发射探头布置示意图；
图3是本发明一个实施例提供的岩石声发射的混合矩张量反演计算方法中的三点弯曲声发射事件震源定位结果示意图；
图4是本发明一个实施例提供的岩石声发射的混合矩张量反演计算方法中的三点弯曲声发射震源机制图。
具体实施方式
以下结合具体实施方案和附图对本发明进行进一步的详细描述。其只意在详细阐述本发明的具体实施方案，并不对本发明产生任何限制，本发明的保护范围以权利要求书为准。
参照图1，本发明提出一种岩石声发射的混合矩张量反演计算方法，包括：
步骤S001，利用最小二乘法进行绝对矩张量反演计算，得到修正前的震源机制解；
步骤S002，引入迭代参数，根据迭代参数计算权重修正值；
步骤S003，通过权重修正值对修正前的震源机制解进行修正，得到修正后的震源机制解。
现有技术中岩石声发射矩张量反演主要用的是绝对矩张量反演和相对矩张量反演。然而，绝对矩张量反演和相对矩张量反演都非常依靠声波传递的格林函数建立的准确性，尤其是在复杂条件(岩石内部不同矿物组成以及非连续结构的存在导致岩石的非均质性)下，P波和S波的到时拾取的准确性会对格林函数的建立造成极大的问题，而且信号和噪声的比例较小会导致矩张量计算结果失真。因此，本实施例在之前的绝对矩张量反演的基础上，结合相对矩张量计算的优势，提出了一种岩石声发射的混合矩张量反演计算方法。
岩石声发射的混合矩张量反演计算方法，是通过引入了迭代参数，对计算的震源机制解进行修正，可以有效地降低由于信号和噪声比例较低而导致的误差，并修正了由于岩石的的非均质性导致的P波传播分散的技术问题，为岩石内部微破裂的震源机制解提出一种较好的计算方法。
本发明提出一种岩石声发射的混合矩张量反演计算方法，包括：
矩张量反演的前提为点源假设，即：破裂处震源的尺度远小于震源和传感器之间的距离以及由震动产生的弹性波波长，此时震源可以被视为一个点源。震源的非弹性变形特征可用矩张量来描述。在传感器k接收到的弹性波的波形位移振幅uk为
Figure BDA0002296873550000051
式(1)中，*号表示卷积运算。
Figure BDA0002296873550000052
假设力作用的时间短暂，则式(1)可变为线性方程：
uk(x,t)＝Gki,jMij＝GM (2)
为了符合在一个集群(空间内)的地震事件，并区分观察位移与计算的矩张量值，将公式(2)的下标进行延长：
ui obs＝GiMi (3)
向量u的一般形式为uijkl，各字母的物理含义为：i是指事件数目，i＝(1,2,3,…,N)，N是指事件簇中的最大的数字；j是指传感器数目，最大值为p；k是指三轴传感器的分量，k＝(1,2,3)；l＝1是指P波。
格林函数Gi是对应于第i个事件的，是由下列元素组成的：
Gijkl1,Gijkl2,Gijkl3,…,Gijkl
第五个下标对应的是矩张量的六个独立分量。基于总的位移是等于位移分量的向量和来计算的，那么每个格林函数的分量可以按如下方式计算：
Figure BDA0002296873550000053
式(4)中，γi11＝(xi1-ξj1)/rij，γi12＝(xi2-ξj2)/rij和γi13＝(xi3-ξj3)/rij，其分别是指震源到接收器的单位方向向量；震源i的位置：xi＝(xi1,xi2,xi3)，传感器j的位置：ξj＝(ξj1,ξj2,ξj3)；ηk1，ηk2和ηk3是描述三轴传感器三个方向的余弦；ρ是值试件密度；α是指P波波速；r是指震源和传感器的距离，r＝|x-ξ|。
本实施例提出的方法是建立在绝对矩张量反演的基础上，利用最小二乘法进行绝对矩张量反演计算，得到修正前的震源机制解的公式为：
Mi＝Gi -1ui obs (5)，其中Gi为格林函数，ui为接收到的弹性波的波形位移振幅，obs代表观测值，Mi为震源机制解。
在进行最小二乘法进行绝对矩张量反演时，由于方程为超定方程(即方程数>未知数数量)，导致上式(5)计算得到的理论值与实际值并不相等，所以需要进行修正。
利用最小二乘法进行绝对矩张量反演计算，得到修正前的震源机制解的公式为：
Mi＝Gi -1ui obs，其中Gi为格林函数，ui为接收到的弹性波的波形位移振幅，obs代表观测值，Mi为震源机制解。
引入迭代参数，根据迭代参数计算权重修正值至少包括：
引入第一迭代参数和第二迭代参数，第一迭代参数的计算公式为：
Figure BDA0002296873550000061
其中
Figure BDA0002296873550000062
为第一迭代系数，median为取中位值，uijkl th为弹性波的波形位移振幅理论值，uijkl obs为弹性波的波形位移振幅实际值，Neq是指所需计算的方程总数。在没有传感器失效的情况下Neq＝N，N是指事件簇中的最大的数字。
计算第二迭代参数的公式为：
Figure BDA0002296873550000063
其中，ωIterNo为第二迭代参数，IterNo为自然数。引入第二迭代参数ωIterNo，这是为了在计算过程中引入迭代系数来提高计算的平滑度。
IterNo从1取到10，ωIterNo＝{0.10,0.12,0.15,0.19,0.25,0.32,0.40,0.50,0.63,0.79,1.0}。在迭代开始之后，迭代系数不断增大，并在最后阶段，其增加的比例明显增大。
根据迭代参数计算权重修正值的公式为：
Figure BDA0002296873550000071
其中，ajkl为权重修正值。
通过权重修正值对修正前的震源机制解进行修正，得到修正后的震源机制解包括：先计算修正后的弹性波的波形位移振幅实际值，再计算修正后的震源机制解；
修正后的弹性波的波形位移振幅实际值的计算公式为：
uijkl new＝uijkl old+ajkl·uijkl old (9)
其中，uijkl old为弹性波的波形位移振幅实际值，uijkl new为修正后的弹性波的波形位移振幅实际值。
根据在一定范围内的事件群中所有情况，计算得出修正权重ajkl，修正之前所得到的弹性波的波形位移振幅实际值。
计算修正后的震源机制解的计算公式为：
Mi new＝Gi -1ui new (10)，其中，ui new为修正后的弹性波的波形位移振幅实际值，Mi new为修正后的震源机制解。
利用公式ui th＝GiMi (11)重新计算弹性波的波形位移振幅理论值，重复引入迭代参数，根据迭代参数计算权重修正值，通过权重修正值对修正前的震源机制解进行修正，得到修正后的震源机制解的步骤，通过多次修正得到多次修正后的震源机制解。
本实施解决了现有技术中由于岩石的的非均质性导致的P波传播分散的技术问题，为岩石内部微破裂的震源机制解提出一种较好的计算方法。
在实际应用时，以岩石三点弯曲实验为例，对岩石破裂过程中进行声发射监测，如图2所示，在岩石表面布置八个压电陶瓷探头，在加载过程中监测岩石内部产生的声发射信号，然后对其进行信号依次进行如下处理。
首先，利用SwansongⅡ信号滤波法对岩石破裂过程中的声发射信号进行滤噪，其滤噪的主要原则基于声发射信号与岩石内部裂纹的相关性。通过该方法滤噪后，可以有效的去除监测信号中一些非相关的信号，提高定位与反演的精度。
第二步，在滤噪的基础上，对声发射信号进行P波到时拾取，这一步是为了得到P波的初到时间和初动振幅。这一步用AR-AIC到时拾取方法，主要原则是：在P波到达的时候，其振幅会发生较大的变化，根据振幅的期望和方差的变化来识别P波的到达时间。
第三步，利用Geiger定位技术，对内部破裂震源进行定位，其原理是依靠P波从震源位置到达不同探头时间的差值。如图3所示，是岩石在三点弯曲条件下的声发射事件定位图。
第四步，根据本发明提出一种岩石声发射的混合矩张量反演计算方法，进行混合矩张量反演计算，得到岩石内部破裂震源的时空演化特征，如图4所示。
本发明提出一种存储介质，所述存储介质存储计算机指令，当计算机执行所述计算机指令时，用于执行上述任一项所述的岩石声发射的混合矩张量反演计算方法。
本发明提出一种岩石声发射的混合矩张量反演计算系统，包括至少一个处理器；以及，与所述至少一个处理器通信连接的存储器；其中，所述存储器存储有可被所述一个处理器执行的指令，所述指令被所述至少一个处理器执行，以使所述至少一个处理器能够：
利用最小二乘法进行绝对矩张量反演计算，得到修正前的震源机制解；
引入迭代参数，根据迭代参数计算权重修正值；
通过权重修正值对修正前的震源机制解进行修正，得到修正后的震源机制解。
以上所述的仅是本发明的原理和较佳的实施例。应当指出，对于本领域的普通技术人员来说，在本发明原理的基础上，还可以做出若干其它变型，也应视为本发明的保护范围。

Claims (10)

1.一种岩石声发射的混合矩张量反演计算方法，其特征在于，包括：
利用最小二乘法进行绝对矩张量反演计算，得到修正前的震源机制解；
引入迭代参数，根据迭代参数计算权重修正值；
通过权重修正值对修正前的震源机制解进行修正，得到修正后的震源机制解。
2.根据权利要求1所述的岩石声发射的混合矩张量反演计算方法，其特征在于，
利用最小二乘法进行绝对矩张量反演计算，得到修正前的震源机制解的公式为：
Mi＝Gi -1ui obs，其中Gi为格林函数，ui为接收到的弹性波的波形位移振幅，obs代表观测值，Mi为震源机制解。
3.根据权利要求1所述的岩石声发射的混合矩张量反演计算方法，其特征在于，
引入迭代参数，根据迭代参数计算权重修正值至少包括：
引入第一迭代参数和第二迭代参数，第一迭代参数的计算公式为：
Figure FDA0002296873540000011
其中
Figure FDA0002296873540000012
为第一迭代系数，median为取中位值，uijkl th为弹性波的波形位移振幅理论值，uijkl obs为弹性波的波形位移振幅实际值，Neq是指所需计算的方程总数。
4.根据权利要求3所述的岩石声发射的混合矩张量反演计算方法，其特征在于，计算第二迭代参数的公式为：
Figure FDA0002296873540000013
其中，ωIterNo为第二迭代参数，IterNo为自然数。
5.根据权利要求4所述的岩石声发射的混合矩张量反演计算方法，其特征在于，根据迭代参数计算权重修正值的公式为：
Figure FDA0002296873540000021
其中，ajkl为权重修正值。
6.根据权利要求5所述的岩石声发射的混合矩张量反演计算方法，其特征在于，
通过权重修正值对修正前的震源机制解进行修正，得到修正后的震源机制解包括：先计算修正后的弹性波的波形位移振幅实际值，再计算修正后的震源机制解；
修正后的弹性波的波形位移振幅实际值的计算公式为：
uijkl new＝uijkl old+ajkl·uijkl old
其中，uijkl old为弹性波的波形位移振幅实际值，uijkl new为修正后的弹性波的波形位移振幅实际值。
7.根据权利要求6所述的岩石声发射的混合矩张量反演计算方法，其特征在于，
计算修正后的震源机制解的计算公式为：
Mi new＝Gi -1ui new其中，ui new为修正后的弹性波的波形位移振幅实际值，Mi new为修正后的震源机制解。
8.根据权利要求7所述的岩石声发射的混合矩张量反演计算方法，其特征在于，
利用公式ui th＝GiMi重新计算弹性波的波形位移振幅理论值，重复引入迭代参数，根据迭代参数计算权重修正值，通过权重修正值对修正前的震源机制解进行修正，得到修正后的震源机制解的步骤，通过多次修正得到多次修正后的震源机制解。
9.一种存储介质，其特征在于，所述存储介质存储计算机指令，当计算机执行所述计算机指令时，用于执行权利要求1至8中任一项所述的岩石声发射的混合矩张量反演计算方法。
10.一种岩石声发射的混合矩张量反演计算系统，其特征在于，包括至少一个处理器；以及，与所述至少一个处理器通信连接的存储器；其中，所述存储器存储有可被所述一个处理器执行的指令，所述指令被所述至少一个处理器执行，以使所述至少一个处理器能够：
利用最小二乘法进行绝对矩张量反演计算，得到修正前的震源机制解；
引入迭代参数，根据迭代参数计算权重修正值；
通过权重修正值对修正前的震源机制解进行修正，得到修正后的震源机制解。
CN201911205651.3A 2019-11-29 2019-11-29 岩石声发射的混合矩张量反演计算方法、存储介质和系统 Pending CN110907538A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201911205651.3A CN110907538A (zh) 2019-11-29 2019-11-29 岩石声发射的混合矩张量反演计算方法、存储介质和系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911205651.3A CN110907538A (zh) 2019-11-29 2019-11-29 岩石声发射的混合矩张量反演计算方法、存储介质和系统

Publications (1)

Publication Number Publication Date
CN110907538A true CN110907538A (zh) 2020-03-24

Family

ID=69820876

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911205651.3A Pending CN110907538A (zh) 2019-11-29 2019-11-29 岩石声发射的混合矩张量反演计算方法、存储介质和系统

Country Status (1)

Country Link
CN (1) CN110907538A (zh)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111734490A (zh) * 2020-05-11 2020-10-02 中南大学 确定巷道支护区域和支护方式的方法及系统
CN112034425A (zh) * 2020-09-15 2020-12-04 中南大学 未知波速的声发射源线性更正定位方法、系统及存储介质
CN113218766A (zh) * 2021-05-27 2021-08-06 江西理工大学 一种基于矩张量分析的岩石起裂应力与损伤应力辨识方法
CN113866278A (zh) * 2021-09-26 2021-12-31 江西理工大学 基于优势声发射震源能量特征的岩石长期强度确定方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105334266A (zh) * 2015-11-11 2016-02-17 华北理工大学 一种岩石声发射源定位方法
CN105891884A (zh) * 2016-06-28 2016-08-24 中国石油大学(北京) 一种微地震震源机制反演方法及装置
CN108195941A (zh) * 2017-12-29 2018-06-22 郑州智谷工业技术有限公司 一种岩石压裂声发射在线监测系统
CN108873063A (zh) * 2017-05-16 2018-11-23 大庆市弘洋石油科技开发有限公司 一种微地震矩张量反演的方法及装置
CN110297040A (zh) * 2019-06-10 2019-10-01 华北水利水电大学 一种基于声发射上升波段采集参数的声发射类型判断方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105334266A (zh) * 2015-11-11 2016-02-17 华北理工大学 一种岩石声发射源定位方法
CN105891884A (zh) * 2016-06-28 2016-08-24 中国石油大学(北京) 一种微地震震源机制反演方法及装置
CN108873063A (zh) * 2017-05-16 2018-11-23 大庆市弘洋石油科技开发有限公司 一种微地震矩张量反演的方法及装置
CN108195941A (zh) * 2017-12-29 2018-06-22 郑州智谷工业技术有限公司 一种岩石压裂声发射在线监测系统
CN110297040A (zh) * 2019-06-10 2019-10-01 华北水利水电大学 一种基于声发射上升波段采集参数的声发射类型判断方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
张凌凡等: "岩石断裂混合矩张量反演与数值分析", 《矿业科学学报》 *

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111734490A (zh) * 2020-05-11 2020-10-02 中南大学 确定巷道支护区域和支护方式的方法及系统
CN111734490B (zh) * 2020-05-11 2022-03-29 中南大学 确定巷道支护区域和支护方式的方法及系统
CN112034425A (zh) * 2020-09-15 2020-12-04 中南大学 未知波速的声发射源线性更正定位方法、系统及存储介质
CN112034425B (zh) * 2020-09-15 2023-03-21 中南大学 未知波速的声发射源线性更正定位方法、系统及存储介质
CN113218766A (zh) * 2021-05-27 2021-08-06 江西理工大学 一种基于矩张量分析的岩石起裂应力与损伤应力辨识方法
CN113866278A (zh) * 2021-09-26 2021-12-31 江西理工大学 基于优势声发射震源能量特征的岩石长期强度确定方法
CN113866278B (zh) * 2021-09-26 2023-12-29 江西理工大学 基于优势声发射震源能量特征的岩石长期强度确定方法

Similar Documents

Publication Publication Date Title
CN110907538A (zh) 岩石声发射的混合矩张量反演计算方法、存储介质和系统
JP2008134182A5 (zh)
CN110132308B (zh) 一种基于姿态确定的usbl安装误差角标定方法
WO2009105964A1 (zh) 一种三维地震记录的叠前似二维变换的方法
Zhu et al. Calculation of dynamic response sensitivity to substructural damage identification under moving load
CN108873063B (zh) 一种微地震矩张量反演的方法及装置
CN108717201A (zh) 一种隧道围岩微震源定位方法
JP3837099B2 (ja) 構造物の損傷推定システムおよびプログラム
CN110703313B (zh) 考虑传感器感度的声发射事件震级获取方法、系统及可读存储介质
CN111400898B (zh) 基于主模态法和应变的阵列天线振动变形预测方法及设备
Zhu et al. Removing mass loading effects of multi-transducers using Sherman-Morrison-Woodbury formula in modal test
CN105204065A (zh) 一种拾取初至波的方法和装置
CN111781639B (zh) 针对obs多分量数据的炮检互易弹性波全波形反演方法
CN110954952B (zh) 一种微震信号初至波类型判别和波速校正方法
CN115508908A (zh) 一种地震面波走时和重力异常联合反演方法与系统
CN112130204B (zh) 一种适用于深埋隧洞围岩波速的现场测试方法
CN115453622A (zh) 一种优化加速的最小二乘逆时偏移成像方法
CN109061723B (zh) 一种隧洞岩爆孕育过程的微震源高精度定位方法及系统
CN112578435A (zh) 岩石超声测试初至拾取方法及系统
CN111736208A (zh) 变权重联合P波和S波初至数据的微震事件Bayes定位方法、系统及介质
Zhong et al. A composite beam integrating an in-situ FPCB sensor membrane with PVDF arrays for modal curvature measurement
CN109975869B (zh) 一种沿地层走向光滑约束的反射波波形反演方法
CN110045014B (zh) 基于贝叶斯学习的Lamb波频散消除方法及其系统
CN109188527B (zh) 滩浅海地区快速建立三维近海底速度模型的方法
CN114964259A (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
RJ01 Rejection of invention patent application after publication

Application publication date: 20200324

RJ01 Rejection of invention patent application after publication