CN104166159B - 四维微地震监测的裂缝形态处理方法和*** - Google Patents
四维微地震监测的裂缝形态处理方法和*** Download PDFInfo
- Publication number
- CN104166159B CN104166159B CN201410335280.1A CN201410335280A CN104166159B CN 104166159 B CN104166159 B CN 104166159B CN 201410335280 A CN201410335280 A CN 201410335280A CN 104166159 B CN104166159 B CN 104166159B
- Authority
- CN
- China
- Prior art keywords
- fracture
- dimensional
- energy
- seismic
- source
- 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.)
- Expired - Fee Related
Links
- 238000012544 monitoring process Methods 0.000 title claims abstract description 43
- 238000000034 method Methods 0.000 claims abstract description 26
- 238000004364 calculation method Methods 0.000 claims abstract description 17
- 238000010276 construction Methods 0.000 claims abstract description 14
- 238000003384 imaging method Methods 0.000 claims abstract description 14
- 206010017076 Fracture Diseases 0.000 claims description 92
- 208000010392 Bone Fractures Diseases 0.000 claims description 90
- 238000003672 processing method Methods 0.000 claims description 24
- 239000011159 matrix material Substances 0.000 claims description 23
- 230000008569 process Effects 0.000 claims description 12
- 238000010606 normalization Methods 0.000 claims description 9
- 238000012545 processing Methods 0.000 claims description 8
- 238000011084 recovery Methods 0.000 claims description 7
- 239000011435 rock Substances 0.000 claims description 6
- 238000011161 development Methods 0.000 claims description 4
- 230000000694 effects Effects 0.000 claims description 4
- 239000012530 fluid Substances 0.000 claims description 3
- 238000012417 linear regression Methods 0.000 claims description 3
- 230000008859 change Effects 0.000 abstract description 9
- 238000005516 engineering process Methods 0.000 description 8
- 230000015572 biosynthetic process Effects 0.000 description 5
- 238000006243 chemical reaction Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 4
- 230000007547 defect Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000001228 spectrum Methods 0.000 description 2
- 238000010521 absorption reaction Methods 0.000 description 1
- 230000004075 alteration Effects 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 230000001186 cumulative effect Effects 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明公开了四维微地震监测的裂缝形态处理方法和***,其裂缝形态处理方法包括:A、选定需监测计算的空间三维地质体范围;B、按照预设分辨率将所述空间三维地质体划分为多个体积单元;C、建立施工区域地层的速度波模型,根据速度波模型进行射线追踪,计算各体积单元与检波器的信号传播时间和路径衰减;D、对目标数据体中的每个体积单元,根据检波器的接收信号恢复原始信号;E计算各体积单元的能量,形成四维成像结果并输出。本发明大大提高了裂缝形态描述的准确性,还能反映出各个时刻裂缝的变化过程。
Description
技术领域
本发明涉及微地震监测技术,特别涉及一种四维微地震监测的裂缝形态处理方法和***。
背景技术
微地震监测技术是一种裂缝监测方法,目前普遍应用于国内各大油气田的开发、规划。微地震事件发生在裂隙之类的断面上,地层内地应力呈各向异性分布。随着外力影响,原有裂缝会变形、扩展,从而使应力松弛,储藏能量的一部分以弹性波的形式释放出来产生小的地震,即微地震。微地震事件一般表现为清晰的脉冲。越弱的微地震事件,其频率越高,能量越小,破裂长度越短。在传播过程中,由于地层结构复杂,岩石介质吸收以及不同地质环境,都会影响微地震信号的能量。
现有的微地震监测解释技术默认探区地层是均匀地层,地下微震波是以直线方式传播至地表被检波器接收到的,在对微地震事件进行反演计算定位时存在较大误差,导致对微地震震源的定位不准确。由于地下裂缝形态是根据微地震事件反映出来的,所以对于最终的裂缝形态的描述也会存在较大误差。
因此,现有技术还有待于改进和发展。
发明内容
鉴于上述现有技术的不足之处,本发明的目的在于提供一种四维微地震监测的裂缝形态处理方法和***,以解决现有技术对裂缝形态的描述存在较大误差的问题。
为了达到上述目的,本发明采取了以下技术方案:
一种四维微地震监测的裂缝形态处理方法,其包括:
A、选定需监测计算的空间三维地质体范围;
B、按照预设分辨率将所述空间三维地质体划分为多个体积单元;
C、建立施工区域地层的速度波模型,根据速度波模型进行射线追踪,计算各体积单元与检波器的信号传播时间和路径衰减;
D、对目标数据体中的每个体积单元,根据检波器的接收信号恢复原始信号;
E、计算各体积单元的能量,形成四维成像结果并输出。
所述的四维微地震监测的裂缝形态处理方法中,在所述步骤D中,还包括:计算微震波震源座标;当有6个接收站接收到同一震源处的微震信号时,地面各接收站的P波到时时间与发震时刻的关系如下:
(1),
对公式(1)进行推导可得到公式(2):
(2),
其中,~表示任意六个接收站的P波到时,是发震时刻;(,,)…(,,)是各接收站的坐标,~的值默认为0;VP是P波速度,(,,)是微震震源的空间坐标;,,,是待求的未知数。
所述的四维微地震监测的裂缝形态处理方法中,所述公式(2)采用矩阵方式求解,将公式(2)按照[A]=[K][B]方式矩阵化,其中,矩阵[A]为:
,
矩阵[B]为:
,
矩阵[K]为:
。
所述的四维微地震监测的裂缝形态处理方法中,从目的层某点到地面的高度的公式如下:
,
微震波震源的空间位置为(,,,)。
所述的四维微地震监测的裂缝形态处理方法中,在所述步骤E中,地震波震源处破裂能量的计算公式:
,
其中,k是被扫描的目标体积中第k个震源点,M是最大检波器数量,,w是设定的时间窗口,fij是矢量迭加K点到所有采集站记录的信号振幅的平方,F是归一化因子。
一种用于实现所述四维微地震监测的裂缝形态处理方法的***,其包括:
扫描模块,用于选定需监测计算的空间三维地质体范围, 按照预设分辨率将所述空间三维地质体划分为多个体积单元;
计算模块,用于建立施工区域地层的速度模型,根据速度模型进行射线追踪,计算各体积单元与检波器的信号传播时间和路径衰减;
恢复模块,对目标数据体中的每个体积单元,根据检波器的接收信号恢复原始信号;
能量图像模块,用于计算各体积单元的能量,形成四维成像结果并输出。
所述的四维微地震监测的裂缝形态处理方法的***中,所述恢复模块还用于计算微震波震源座标;当有6个接收站接收到同一震源处的微震信号时,地面各接收站的P波到时时间与发震时刻的关系如下:
(1),
对公式(1)进行推导可得到公式(2):
(2),
其中,~表示任意六个接收站的P波到时,是发震时刻;(,,)…(,,)是各接收站的坐标,~的值默认为0;VP是P波速度,(,,)是微震震源的空间坐标;,,,是待求的未知数。
所述的四维微地震监测的裂缝形态处理方法的***中,所述恢复模块对所述公式(2)采用矩阵方式求解,将公式(2)按照[A]=[K][B]方式矩阵化,其中,矩阵[A]为:
,
矩阵[B]为:
,
矩阵[K]为:
。
所述的四维微地震监测的裂缝形态处理方法的***中,从目的层某点到地面的高度的公式如下:
,
微震波震源的空间位置为(,,,)。
所述的四维微地震监测的裂缝形态处理方法的***中,所述能量图像模块计算地震波震源处破裂能量采用的计算公式为:
,
其中,k是被扫描的目标体积中第k个震源点,M是最大检波器数量, w是设定的时间窗口,fij是矢量迭加K点到所有采集站记录的信号振幅的平方,F是归一化因子。
相较于现有技术,本发明提供的四维微地震监测的裂缝形态处理方法和***,选定需监测计算的空间三维地质体范围, 按照预设分辨率将所述空间三维地质体划分为多个体积单元,建立施工区域地层的速度波模型,根据微震波定位理论定位出微震波震源、即破裂点的坐标;通过记各破裂点的破裂能量,形成裂缝形态的各时段破裂能量图;大大提高了裂缝形态描述的准确性,还能反映出各个时刻裂缝的变化过程。
附图说明
图1为本发明四维微地震监测的裂缝形态处理方法的方法流程图。
图2为本发明四维微地震监测的裂缝形态处理方法中空间三维立方体的示意图。
图3为本发明四维微地震监测的裂缝形态处理方法中震源与检波器之间的射线示意图。
图4为本发明四维微地震监测的裂缝形态处理方法中各时段破裂能量示意图。
图5为本发明四维微地震监测的裂缝形态处理***的结构框图。
具体实施方式
本发明提供一种四维微地震监测的裂缝形态处理方法和***,在基于现有裂缝形态处理的基础上,增加速度模型理论,对所在探区进行分层,提高了微震源定位准确性,通过微地震事件反演定位,大大提高了裂缝形态描述的准确性,又可以反映出各个时刻裂缝的变化过程。为使本发明的目的、技术方案及效果更加清楚、明确,以下参照附图并举实施例对本发明进一步详细说明。应当理解,此处所描述的具体实施例仅用以解释本发明,并不用于限定本发明。
本实施例应用于油田目的井压裂监测,通过获取震源的几何位置来大致确定地下裂缝的形态,再获取震源的破裂能量来细化裂缝形态;请参阅图1,所述的四维微地震监测的裂缝形态处理方法包括以下步骤:
S100、选定需监测计算的空间三维地质体范围。默认为是由施工目的层段至地表的空间三维地质体。
S200、按照预设分辨率将所述空间三维地质体划分为多个体积单元。为后面速度波模型的建立以及通过微震信号反演得出地下能量场的分布变化做准备。
S300、建立施工区域地层的速度波模型,根据速度波模型进行射线追踪,计算各体积单元与检波器的信号传播时间和路径衰减。为后面地下微震源的定位以及破裂能量计算做准备。
S400、对目标数据体中的每个体积单元,根据检波器的接收信号恢复原始信号。根据地面检波器接收到的微震信号,结合S300计算出的各体积单元与检波器的信号传播时间和路径衰减,得出地下微震源的空间分布。
S500、计算各体积单元的能量,形成四维成像结果并输出。由于地下某一微震波传播至地面时会被多个检波器接收到,将各检波器接收到的微震信号进行反演得出某一点的微震信号,若微地震发生在这一点,则所得微地震信号会在此点叠加,若不是则所得微震信号会相互抵消。
其中,所述需要监测的空间三维地质体通常以井水平井压裂层段参考,如距水平段中心1000M(X轴)*1000M(Y轴)*3000M(Z轴),M表示米。按照预设分辨率(如5M*5M*5M)将所述空间三维地质体划分为多个体积单元101(也可叫网格节点),如图2中所示的各个方格。空间三维地质体的顶面为地表最高点海拔的平面。监测的目的层为对应Z轴坐标距顶面的深度值。
体积单元划分完成后,即可开始建立施工区域地层(即从目的层到地表的所有区域)的速度波模型。则在所述步骤S300中,速度波模型的建立方法包括:
步骤1、在有井的区块,收集井声波和VSP测井(专用于求取准确地震速度的一种测井,可获得准确的时深关系,进而得到准确的地震平均速度和地震层速度),再结合时深表资料计算出井孔层速度Vint(相当于纵波);以提高井的纵向分辨率,其中,井孔层速度Vint的公式如下:
,其中,
表示厚度; 表示通过该段厚度地层的双程时间,可用VSP资料直接求取;Hb、Ha分别为某一层段的底部深度值、顶部深度值;Tob、Ton为通过地层底部与顶部的反射双轨时间。
步骤2、在无井的区块,通过速度普建立速度模型(将地层分层,同一层的传播速度相同),并根据速度普-DIX公式(常规层速度求取公式)将速度普转换为层速度,层速度沿构造层序内插得到层速度数据体,然后将层速度数据体转换为平均速度体(相当于横波),其公式如下:
其中,、是相邻两地层的平均速度,t01、t02是从基准面至相邻两地层界面的反射双程时间。
本实施例中将地层分为多层,认为同一层中的传播速度相同。通过计算平均速度体来提高横向分辨率。
步骤3、井与速度谱协同构建高精度的平均速度体V3,把井的纵向分辨率和速度谱的横向分辨率结合起来,在构造层序的约束下,得到高精度的速度波模型。
现有的微地震监测技术都是默认微震波是以直线方式由震源处传播至地面检波器102的,而实际上由于地层结构复杂,地震波往往是以折射方式经地层的层层衰减才传至地面被检波器接收到的。本实施例通过引入速度波模型,在施工区域建立速度波模型,可以精确划分地层速度,对微震波震源进行反演定位的准确性大大提高。之后,再根据速度波模型进行射线追踪,计算各体积单元与检波器的信号(即微震波信号)传播时间和路径衰减,如图3所示。图3中仅示出微震波从震源处传播至地面某检波器的路线示意图。图3中需要监测的空间三维立方体10中、在各个层上分布了多个震源(以空间三维立方体10内部的圆点表示)。基于不同层地质结构不同,各震源处微震信号是以折射波形式传输至地面检波器的。由于大地坐标与数学坐标并不一致,在图3中,N表示正北方向、相当于X轴;E表示正东方向、相当于Y轴;D表示深度,相当于Z轴。通过确定震源映射在X、Y、Z轴上的坐标即可获得震源的几何位置。
所述步骤S400中,所述目标数据体即选定需监测计算的空间三维地质体范围。在该空间三维地质体范围内,有的体积单元有微震波震源产生,有的体积单元没有。通过对所有的体积单元进行计算,有微震波震源的先定位震源坐标、再计算震源处破裂能量来描述裂缝形态;无微震波震源的保持当前能量状况(即无能量反应)。将所有体积单元的计算结果按照时间顺序整合起来即可获得裂缝的变化过程。则在步骤S400中,还包括计算微震波震源座;本实施例采用微震波震源定位原理来判定微震波震源坐标。
假设有6个接收站接收到同一震源处信号,首先需要计算地面各接收站的P波到时时间与发震时刻的关系如下:
(1),
对上述公式(1)进行推导可得到公式(2),如下所示:
(2); 其中,~表示任意六个接收站的P波到时,是发震时刻;(,,)…(,,)是各接收站的坐标,一般~的值默认为0;VP是P波速度,(,,)是微震震源的空间坐标;,,,是待求的未知数。求解这一系列方程组(实际转化为矩阵求解),就可确定微震震源位置。在一段时间内,许许多多微震震源的空间分布可以描述裂缝实时轮廓,进而给出裂缝的方位、长度、高度、产状及地应力方向等地层参数。
上述公式(2)可写成标准系数矩阵形式,求解矩阵的方法很多,可以解出,,;再把,,代入公式(1)中就可以得出,就是相对压裂深度的裂缝高度,由于计算过程的累积作用,高度误差较大。
本实施例中,可将上述公式(2)按照[A]=[K][B]方式矩阵化,其中,矩阵[A]写为:
,
矩阵[B]写为:
,
矩阵[K]写为:
。
在具体实施时,也可以将上述公式(2)按照[K]=[A]×[B]-1 方式矩阵化,[B]-1是[B]的逆矩阵。此处对此不作详述。通过上述公式即可求得,,值。
本实施例中确定深度的方法是综合考虑各接收站走时的方法,即对以下任一组数据作为走时的函数进行线性回归,回归常数即为相对观测段(即从目的层某点到地面的高度,其公式如下:
。
依据上述计算过程即可确定微震波震源的空间位置(,,,)。
微震波震源的空间分布能反应出地下裂缝的大体形态,裂缝从无到有,至不断扩展的变化形态则可通过破裂能量的变化过程形象的描述。由于微震波发出至传播到地表检波器时,传输过程中存在一定能量损耗和信号衰减,还有其他外界环境干扰。因此,检波器接收到的信号包括了剩余的微震波和杂波,在步骤S400中,对目标数据体中的每个体积单元,根据接收信号来恢复原始信号(即微震波原始的信号能量),也即是反演出微震波震源能量大小。则在步骤S500中,微震波震源处破裂能量Sc(k)的计算公式为:
,
其中,k是被扫描的目标体积中第k个震源点,M是最大检波器数量, w是设定的时间窗口,fij是矢量迭加K点到所有采集站记录的信号振幅的平方,F是归一化因子。
对所有的接收站在时间窗口w内记录的f求和,并适当使用归一化因子(normalization_factor), 即得到k点的破裂能量Sc(k),扫描目标体积中所有的数值点,即获得一段时间内破裂能量的分布。破裂传播到一个台站的较大振幅(主要能量)集中在一个时间段内。如果地震波速度分布已知,只要正确的移动(即射线追踪)各接收站的记录并迭加其振幅,破裂点处的破裂能量Sc(k)将会增强。
通过对监测记录的能量数据进行反演叠加,得出地下能量场的分布趋势,从而立体的描述出裂缝的产状、方位、长度等信息(即四维成像结果)。如图4所示。从图4中可以看出,a时刻刚开始无地下能量场反应,之后在b时刻出现地下能量反应,c时刻地下能量反应随压裂施工进度的进行慢慢扩大,d时刻继续扩大;e时刻继续扩大且出现两个较高的能量区域,f时刻裂缝继续扩大且高能量区域也扩大、能量增加。图4清晰地反映出各个时刻地下能量的变化过程,亦反映出地下裂缝的变化过程,大大方便了微地震监测。
现有的微地震监测技术仅以所测得的微震事件(判定的地下微震波震源)的分布来判定地下裂缝的几何形态,但实际上由于地层结构的非均质性,同一区域的地下应力场分布并不均匀以及天然裂缝的存在,这就造成了压裂施工时,压裂段近端的岩石未破裂可远端的岩石优先破裂的现象产生。针对现有有微地震技术的这一缺漏,本实施例将微地震监测技术与影像成图技术结合起来,运用无源地震的微地震三分量数据,进行多波(纵波和横波)振幅属性分析,并采用相关体数据计算处理方法,得出监测期内各时间域三维空间体地下储层岩石破裂和高压流体活动释放的能量分布情况。
本发明实施例还相应提供一种用于实现所述四维微地震监测的裂缝形态处理方法的***,如图5所示,所述***包括:
扫描模块110,选定需监测计算的空间三维地质体范围,按照预设分辨率将所述空间三维地质体划分为多个体积单元。一般为监测施工层段至地表,按照(5M*5M*5M)的分辨率将空间三维地质体划分为多个体积单元。
计算模块120,用于建立施工区域地层的速度波模型,根据速度波模型进行射线追踪,计算各体积单元与检波器的信号传播时间和路径衰减。
恢复模块130,用于对目标数据体中的每个体积单元,根据检波器的接收信号恢复原始信号;
能量图像模块140,用于计算各体积单元的能量,形成四维成像结果。相当于反演出微震波震源能量大小;记录破裂点不同时段的破裂能量,形成裂缝形态的各时段破裂能量图,从而反应裂缝发育过程中各时刻地下能量场分布。
由于上述***中各个部分的功能已经在上述方法实施例中进行了详细介绍,此处不再赘述。
综上所述,本发明按照预设网格将需要监测的空间区域划分为多个体积单元,按照一定分辨率将地质体划分为多个体积单元;根据速度模型进行射线追踪,计算各体积单元与检波器的信号传播时间和路径衰减;对目标数据体中的每个体积单元,根据检波器的接收信号,恢复原始信号;计算各体积单元的能量,形成四维(XYZT)成像结果。
可以理解的是,对本领域普通技术人员来说,可以根据本发明的技术方案及其发明构思加以等同替换或改变,而所有这些改变或替换都应属于本发明所附的权利要求的保护范围。
Claims (6)
1.一种四维微地震监测的裂缝形态处理方法,其特征在于,包括:
A、选定需监测计算的空间三维地质体范围;
B、按照预设分辨率将所述空间三维地质体划分为多个体积单元;
C、建立施工区域地层的速度波模型,根据速度波模型进行射线追踪,计算各体积单元与检波器的信号传播时间和路径衰减;
D、对目标数据体中的每个体积单元,根据检波器的接收信号恢复原始信号;
E、计算各体积单元的能量,形成四维成像结果并输出;
所述步骤D中,对所有的体积单元进行计算,有微震波震源的先定位震源坐标、再计算震源处破裂能量来描述裂缝形态;无微震波震源的保持当前能量状况;
在所述步骤D中,还包括:计算微震波震源座标;当有6个接收站接收到同一震源处的微震信号时,地面各接收站的P波到时时间与发震时刻的关系如下:
(1);
其中,~表示任意六个接收站的P波到时,是发震时刻;(,,)…(,,)是各接收站的坐标,~的值默认为0;VP是P波速度,(,,)是微震震源的空间坐标;,,,是待求的未知数;
对任一组数据作为走时的函数进行线性回归,回归常数即为从目的层某点到地面的高度,其公式如下:
;
微震波震源的空间位置为(,,,);
在所述步骤E中,地震波震源处破裂能量的计算公式:
,
其中,k是被扫描的目标体积中第k个震源点,M是最大检波器数量, w是设定的时间窗口,fij是矢量迭加K点到所有采集站记录的信号振幅的平方,F是归一化因子;
对所有的接收站在时间窗口w内记录的f求和,并使用归一化因子, 即得到k点的破裂能量Sc(k),扫描目标体积中所有的数值点,即获得一段时间内破裂能量的分布;
运用无源地震的微地震三分量数据,进行纵波和横波振幅属性分析,并采用相关体数据计算处理方法,得出监测期内各时间域三维空间体地下储层岩石破裂和高压流体活动释放的能量分布情况;记录破裂点不同时段的破裂能量,形成裂缝形态的各时段破裂能量图,从而立体的描述出裂缝的产状、方位、长度的信息,反应裂缝发育过程中各时刻地下能量场分布。
2.根据权利要求1所述的四维微地震监测的裂缝形态处理方法,其特征在于,
对公式(1)进行推导可得到公式(2):
(2)。
3. 根据权利要求2所述的四维微地震监测的裂缝形态处理方法,其特征在于,所述公式(2)采用矩阵方式求解,将公式(2)按照[A]=[K][B]方式矩阵化,其中,矩阵[A]为:
,
矩阵[B]为:
,
矩阵[K]为:
。
4.一种用于实现权利要求1所述四维微地震监测的裂缝形态处理方法的***,其特征在于,包括:
扫描模块,用于选定需监测计算的空间三维地质体范围, 按照预设分辨率将所述空间三维地质体划分为多个体积单元;
计算模块,用于建立施工区域地层的速度模型,根据速度模型进行射线追踪,计算各体积单元与检波器的信号传播时间和路径衰减;
恢复模块,对目标数据体中的每个体积单元,根据检波器的接收信号恢复原始信号;对所有的体积单元进行计算,有微震波震源的先定位震源坐标、再计算震源处破裂能量来描述裂缝形态;无微震波震源的保持当前能量状况;
所述恢复模块还用于计算微震波震源座标;当有6个接收站接收到同一震源处的微震信号时,地面各接收站的P波到时时间与发震时刻的关系如下:
(1);
其中,~表示任意六个接收站的P波到时,是发震时刻;(,,)…(,,)是各接收站的坐标,~的值默认为0;VP是P波速度,(,,)是微震震源的空间坐标;,,,是待求的未知数;
对任一组数据作为走时的函数进行线性回归,回归常数即为从目的层某点到地面的高度,其公式如下:
;
微震波震源的空间位置为(,,,);
能量图像模块,用于计算各体积单元的能量,形成四维成像结果并输出;地震波震源处破裂能量的计算公式:
,
其中,k是被扫描的目标体积中第k个震源点,M是最大检波器数量, w是设定的时间窗口,fij是矢量迭加K点到所有采集站记录的信号振幅的平方,F是归一化因子;
对所有的接收站在时间窗口w内记录的f求和,并使用归一化因子, 即得到k点的破裂能量Sc(k),扫描目标体积中所有的数值点,即获得一段时间内破裂能量的分布;
运用无源地震的微地震三分量数据,进行纵波和横波振幅属性分析,并采用相关体数据计算处理方法,得出监测期内各时间域三维空间体地下储层岩石破裂和高压流体活动释放的能量分布情况;记录破裂点不同时段的破裂能量,形成裂缝形态的各时段破裂能量图,从而立体的描述出裂缝的产状、方位、长度的信息,反应裂缝发育过程中各时刻地下能量场分布。
5.根据权利要求4所述的四维微地震监测的裂缝形态处理方法的***,其特征在于,
对公式(1)进行推导可得到公式(2):
(2)。
6. 根据权利要求5所述的四维微地震监测的裂缝形态处理方法的***,其特征在于,所述恢复模块对所述公式(2)采用矩阵方式求解,将公式(2)按照[A]=[K][B]方式矩阵化,其中,矩阵[A]为:
,
矩阵[B]为:
,
矩阵[K]为:
。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410335280.1A CN104166159B (zh) | 2014-07-15 | 2014-07-15 | 四维微地震监测的裂缝形态处理方法和*** |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410335280.1A CN104166159B (zh) | 2014-07-15 | 2014-07-15 | 四维微地震监测的裂缝形态处理方法和*** |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104166159A CN104166159A (zh) | 2014-11-26 |
CN104166159B true CN104166159B (zh) | 2015-10-21 |
Family
ID=51910058
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410335280.1A Expired - Fee Related CN104166159B (zh) | 2014-07-15 | 2014-07-15 | 四维微地震监测的裂缝形态处理方法和*** |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104166159B (zh) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105891886A (zh) * | 2014-12-08 | 2016-08-24 | 北京恒欣泰克科技有限公司 | 一种基于三维地层模型的地质体自动分析方法 |
CN105929445B (zh) * | 2016-04-18 | 2018-05-18 | 北京京援伟达技术有限公司 | 微破裂向量扫描方法 |
CN106054242B (zh) * | 2016-05-04 | 2018-11-13 | 中国地质大学(北京) | 三维各向异性衰减介质波场模拟方法 |
CN107784383A (zh) * | 2016-08-31 | 2018-03-09 | 王淳讙 | 地层监测数据四维实境还原***及其方法 |
CN110414675B (zh) * | 2019-09-02 | 2022-05-27 | 中北大学 | 一种基于深度学习的地下浅层震源定位方法 |
CN112727423B (zh) * | 2020-12-30 | 2022-11-18 | 中国海洋石油集团有限公司 | 精细评价海上砂岩储层水力压裂裂缝高度的方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102262220A (zh) * | 2011-04-28 | 2011-11-30 | 中南大学 | 一种基于非线性拟合的微震源或声发射源的定位方法 |
CN102565855A (zh) * | 2012-01-02 | 2012-07-11 | 吉林大学 | 油田压裂地面微地震数据处理方法 |
CN102841373A (zh) * | 2012-08-23 | 2012-12-26 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 基于方位角约束的微地震事件定位方法 |
CN103336297A (zh) * | 2013-06-03 | 2013-10-02 | 北京京援伟达技术有限公司 | 微破裂向量扫描方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US10061046B2 (en) * | 2012-02-06 | 2018-08-28 | Ion Geophysical Corporation | Integrated passive and active seismic surveying using multiple arrays |
-
2014
- 2014-07-15 CN CN201410335280.1A patent/CN104166159B/zh not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102262220A (zh) * | 2011-04-28 | 2011-11-30 | 中南大学 | 一种基于非线性拟合的微震源或声发射源的定位方法 |
CN102565855A (zh) * | 2012-01-02 | 2012-07-11 | 吉林大学 | 油田压裂地面微地震数据处理方法 |
CN102841373A (zh) * | 2012-08-23 | 2012-12-26 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 基于方位角约束的微地震事件定位方法 |
CN103336297A (zh) * | 2013-06-03 | 2013-10-02 | 北京京援伟达技术有限公司 | 微破裂向量扫描方法 |
Also Published As
Publication number | Publication date |
---|---|
CN104166159A (zh) | 2014-11-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104166159B (zh) | 四维微地震监测的裂缝形态处理方法和*** | |
Toomey et al. | Microearthquakes beneath median valley of Mid‐Atlantic Ridge near 23° N: Tomography and tectonics | |
EP1946129B1 (en) | Continuous reservoir monitoring for fluid pathways using 3d microseismic data | |
Kent et al. | Distribution of magma beneath the East Pacific Rise between the Clipperton transform and the 9 17′ N Deval from forward modeling of common depth point data | |
AU667253B2 (en) | Method for simulating crosswell seismic data | |
EP2255222B1 (en) | Method for imaging the earth's subsurface using passive seismic sensing | |
CN104199090B (zh) | 一种地面监测微地震定位的速度模型构建与求解方法 | |
RU2010103987A (ru) | Способы и системы для обработки микросейсмических данных | |
US10670754B2 (en) | System and method for processing microseismic data | |
CN102213769A (zh) | 一种利用三维垂直地震剖面资料确定各向异性参数的方法 | |
CN101086535A (zh) | 表征油藏随时间演变的方法和程序 | |
CN109884710B (zh) | 针对激发井深设计的微测井层析成像方法 | |
CN105510880A (zh) | 一种基于双差法的微地震震源定位方法 | |
AU2014407527B2 (en) | Integrating vertical seismic profile data for microseismic anisotropy velocity analysis | |
GB2505042A (en) | Seismic trace attribute | |
CN109557582A (zh) | 一种二维多分量地震资料偏移成像方法及*** | |
CN103105622B (zh) | 基于数据库技术的同型波时差定位方法 | |
CN104199088B (zh) | 一种提取入射角道集的方法及*** | |
CN102877828A (zh) | 一种三维多井联合井地ct成像方法 | |
CN109100792A (zh) | 基于台站与三维地震联合采集资料的速度反演方法 | |
CN110967751B (zh) | 基于地面浅井监测的微地震事件的定位方法及存储介质 | |
CN102466818B (zh) | 一种利用井间地震数据对各向异性介质成像的方法 | |
Wang et al. | Reverse time migration imaging of tunnels via the finite element method using an unstructured mesh | |
CN110780345A (zh) | 一种隧道超前地震探测地震数据的三维速度分析方法 | |
CN112444885B (zh) | 基于走时相等的微测井和层析反演浅表层联合建模方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20151021 Termination date: 20190715 |
|
CF01 | Termination of patent right due to non-payment of annual fee |