CN117687077A - 利用das和检波器混合阵列监测微震的方法和系统 - Google Patents
利用das和检波器混合阵列监测微震的方法和系统 Download PDFInfo
- Publication number
- CN117687077A CN117687077A CN202311650833.8A CN202311650833A CN117687077A CN 117687077 A CN117687077 A CN 117687077A CN 202311650833 A CN202311650833 A CN 202311650833A CN 117687077 A CN117687077 A CN 117687077A
- Authority
- CN
- China
- Prior art keywords
- source
- detector
- monitoring
- detectors
- das
- 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
Links
- 238000012544 monitoring process Methods 0.000 title claims abstract description 45
- 238000000034 method Methods 0.000 title claims abstract description 33
- 239000013307 optical fiber Substances 0.000 claims abstract description 24
- 238000004364 calculation method Methods 0.000 claims abstract description 21
- 238000012937 correction Methods 0.000 claims description 7
- 238000005070 sampling Methods 0.000 claims description 6
- 230000015572 biosynthetic process Effects 0.000 claims description 5
- 238000004458 analytical method Methods 0.000 claims description 4
- 238000005259 measurement Methods 0.000 claims description 4
- 238000006243 chemical reaction Methods 0.000 claims description 2
- 230000010287 polarization Effects 0.000 abstract 1
- 230000006870 function Effects 0.000 description 7
- 230000000694 effects Effects 0.000 description 4
- 239000013598 vector Substances 0.000 description 4
- 238000010586 diagram Methods 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 230000008859 change Effects 0.000 description 2
- 230000001427 coherent effect 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
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 239000000835 fiber Substances 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 239000002245 particle Substances 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
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/16—Receiving elements for seismic signals; Arrangements or adaptations of receiving elements
- G01V1/18—Receiving elements, e.g. seismometer, geophone or torque detectors, for localised single point measurements
- G01V1/181—Geophones
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/16—Receiving elements for seismic signals; Arrangements or adaptations of receiving elements
- G01V1/20—Arrangements of receiving elements, e.g. geophone pattern
-
- 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/288—Event detection in seismic signals, e.g. microseismics
Landscapes
- Engineering & Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- Remote Sensing (AREA)
- Environmental & Geological Engineering (AREA)
- Acoustics & Sound (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Business, Economics & Management (AREA)
- Emergency Management (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明公开了一种利用DAS和检波器混合阵列监测微震的系统和方法,包括:使用分布式声波传感器测量的信号,使用波形相似法获取震源位置,得到震源距离和深度。沿光纤或在光纤尾端布设三分量检波器,计算偏振角,得到方位信息。根据光纤传感器和三分量检波器得到的结果进行综合计算,得到微震事件三维空间位置。本发明的优点是:定位精度高,具有高准确性和可靠性,能够实时监测,提高了监测的全面性和有效性,成本低。
Description
技术领域
本发明涉及微震监测技术领域,特别涉及一种利用分布式声波传感器和检波器混合阵列微震监测系统和方法。
背景技术
分布式声波传感(DAS,Distributed Acoustic Sensing)的工作原理:DAS是一种传感技术,它通过收集和分析地面或井下产生的声音,这些声波会被光纤捕捉到。利用相干光时域反射测量的原理,将相干短脉冲激光注入到光纤中,当有外界振动作用于光纤上时,由于弹光效应,会微小地改变纤芯内部结构,从而导致背向瑞利散射信号的变化,使得接收到的反射光强发生变化,通过检测井下事件前后的瑞利散射光信号的强度变化。
传统的微震监测分为地面监测和井下监测。在井下监测中,高灵敏度传感器部署在靠近震源的钻孔中,以便在靠近阵列的地方记录更多微震事件,以最大限度地减少信号衰减和背景噪声。这可以检测小震级微震活动,其信噪比足以从稀疏检波器阵列确定源位置。最终检测范围仅限于监测位置周围的区域。
利用传统检波器的监测数据,记录传感器和微震位置之间的方向为由直接P波和/或S波到达的粒子运动确定。
但是利用传统检波器具有以下缺点:
1、布设检波器数目少时,使用井中单一垂直检波器排列,定位误差大,垂直方向定位误差十几米,水平方位误差几十米,使用地面二维检波器排列,垂向方位定位误差大于40米,水平方位误差在10米左右。
2、布设检波器数量多,定位误差会降低,但布设检波器的成本会上升。
使用波形叠加法,对每个波形进行叠加,不依赖于P波或S波的初至。
该方法使用波形进行叠加。在数据处理环节,使用P波进行计算,在计算时叠加波形的振幅值。但是振幅值有正有负,如果简单叠加会出现抵消的情况,该方法在计算之前选择使用走时校正振幅的极性,提高计算结果的准确性。
但是,波形叠加法使用传统三分量检波器,布设成本较高。
使用普通波形叠加法进行波形的叠加,单个事件记录的幅度在不同检波器上具有不同的极性和振幅,如果不进行极性的校正,直接叠加会出现正振幅和负振幅的互相抵消,进行的改进方法是使用极性校正后的振幅进行叠加,但是会引入噪声,局部奇异值的振幅会影响波形的信噪比。具体来讲,只有少数甚至一个高振幅噪声迹线可能会导致高叠加振幅,表明出现虚假检测,即所谓的“误报”。
发明内容
本发明针对现有技术的缺陷,提供了一种利用DAS和检波器混合阵列监测微震的方法和系统。
为了实现以上发明目的,本发明采取的技术方案如下:
一种利用DAS和检波器混合阵列监测微震的系统,包括:分布式声波传感器和检波器;
混合阵列布设方式如下:
分布式声波传感器的光纤水平布设在地表上,光纤标距为3m,三分量检波器沿光纤等间隔布设在水平地表;
作为优选,混合阵列布设方式如下:
分布式声波传感器的光纤垂直布设,沿垂直井方向,在光纤尾端布设3-5个三分量检波器,在井下的光纤,标距为3m。
本发明还公开了一种利用DAS和检波器混合阵列监测微震的方法,包括以下步骤:
步骤一,通过三分量检波器中的两个水平分量X和Y进行震源方位角的计算。使用互相关算子计算X分量和震源之间的夹角θ,并将其转换为相对于北向的方位角φr。
通过检查第一个初至点的极性,消除可能存在的两个潜在的检波器趋势,以消除歧义。
在存在任意偏差的观测井中,通过考虑与偏差相切的线l来解决问题,这需要使用球坐标参数井的倾角θw和井的水平方向φw,以及检波器的坐标xr,yr,zr。
利用参数l的方向定义一个垂直于井的平面,并在该平面上引入伪x、y和z轴,建立一个相对于井在某一点的坐标系。
为了进行地震检波器方向的分析,必须将震源坐标投影到定义的平面上,通过计算源坐标和伪x和y轴的内积来实现,并定义了一个伪源-检波器方位角φs'。最后,将φs'加到θ上将得到检波器方向方位角,相对于伪y轴。
步骤二,使用分布式声波传感器监测数据,在垂直剖面上进行相似系数S(k)的计算,并建立震源扫描网格。震源扫描网格的基本思想是将空间划分成规则的网格,根据检测到的到时进行波形校正,并通过相似系数S(k)来确定震源位置。
根据三分量检波器计算的方位信息,得到最终的震源坐标。
进一步地,震源方位角的计算如下式:
是零滞后的互相关算子,使用互相关算子计算X分量和震源之间的夹角θ如
式2:
进一步地,所述转换为相对于北向的方位角φr,公式如下:
φr=φs+θ (3)
φs是垂直井时的震源方位角。对于垂直井,通过将θ添加到源-检波器方位角φs相对于正北的角度来实现,并且假设检波器的水平分量在x-y平面。因此,在解震源-检波器方位角时,需使用震源位置的x和y坐标,如下所示:
进一步地,所述偏差相切的线l,利用球坐标,用参数化的方式表达为:
式中,θw为井相对于垂直方向的倾斜角度。
φw为井相对于正x轴的水平方向的角度。
xr,yr,zr为检波器在三维空间中的坐标,分别表示x、y和z坐标。
进一步地,建立一个相对于井在某一点的坐标系,具体如下:
进一步地,所述计算相似系数S(k)的公式如下:
式中,fij是第i道第j个采样点,M是参与计算的地震道数目,N为时窗大小,k为时窗中心位置。
进一步地,所述震源扫描网格的基本思想是把空间划分成规则的网格,在震源扫描网格中,以每个划分的网格为基本单元,在基本单元网格中,根据检波器拾取的到时,以最早到的初至时间为标准,进行波形校正,把时差对每个波形进行整体偏移,调整时差之后,使用公式7的相似函数计算,在不同震源扫描网格中进行滑动网格计算,其相似函数取最大值时,网格中的位置是震源位置。通过比较计算出的S(k)值,可以判断震源位置。
进一步地,计算得到最终的震源坐标,公式如下:
式中,(x1,y1,z1)为某一特定垂直剖面中,震源位置坐标,(x2,y2,z2)是根据三分量检波器计算出的方位信息,计算得到的最终的震源坐标。
进一步地,在步骤二之后,计算出震源定位的均方根误差ΔTRMS和速度加权的均方根误差ERMS,评估定位结果震源位置的准确性和可靠性,具体计算公式如下:
震源定位测量的走时与理论走时之间的误差,计算方法如公式8;
定位误差是已知地层速度模型的情况下进行计算,计算方法如公式9,
式中,NP为P波在分布式声波传感器阵列中到达的数量,NS是S波在分布式声波传感器阵列中到达的数量。是各道P波的监测实际走时与理论走时的误差,/>是各道S波的监测实际走时与理论走时的误差。
式中,VP是地层中P波速度,VS是地层中S波速度。
与现有技术相比,本发明的优点在于:
1.高定位精度:通过使用三分量检波器计算震源方位角,并结合分布式声波传感器监测数据,在垂直剖面上计算相似系数来确定震源位置,可以实现高精度的微震监测和定位。
2.准确性和可靠性:通过计算均方根误差来评估监测结果的准确性和可靠性,从而提高监测数据的可信度和实用性。
3.实时监测:该方法可以实时监测并定位微震活动,有助于实时掌握掌握地下微震活动情况,分析压裂裂缝扩展情况,提供了重要的实时监测和预警手段。
4.综合利用多种数据:该方法结合了三分量检波器和分布式声波传感器监测数据,综合利用多种数据,提高了监测的全面性和有效性。
5.成本低:布设方式简单,易于实现,且应用环境广泛,没有限制。
附图说明
图1是本发明实施例合阵列微震监测系统技术路线图。
图2是本发明实施例混合阵列水平方式布设图。
图3是本发明实施例混合阵列垂直方式布设图。
图4是本发明实施例混合阵列水平布设时的震源扫描网格示意图。
图5是本发明实施例混合阵列垂直布设时的震源扫描网格示意图。
具体实施方式
为使本发明的目的、技术方案及优点更加清楚明白,以下根据附图并列举实施例,对本发明做进一步详细说明。
如图1所示,一种利用DAS和检波器混合阵列监测微震的系统,包括:分布式声波传感器和检波器;
混合阵列布设方式有两种:
如图2所示,分布式声波传感器的光纤水平布设在地表上,三分量检波器沿光纤均匀分布,例如2km的水平光纤,标距为3m(等价于3m一个检波器),三分量检波器每隔200m布设一个,主要是为了控制三分量检波器可以接受到震源信号,为定位微震事件提供方位角信息。
如图3所示,分布式声波传感器的光纤垂直布设,沿垂直井方向,在光纤尾端布设3-5个三分量检波器,在井下的光纤,标距仍然设置为3m(等价于3m一个检波器),三分量检波器用于为定位微震事件提供方位角信息。
一种混合阵列微震监测方法,包括下面步骤:
一、根据三分量检波器监测的数据,极化震源方位角:
X和Y是三分量检波器中的两个水平分量,可根据XY水平分量进行震源方位角的计算:
是零滞后的互相关算子,互相关的计算公式如式2,θ是X分量和震源之间的夹角。
互相关在本质上是两个函数做内积运算。即向量内积的连续形式。其在线性空间角度上的意义是:一个向量在另一个向量上的投影,内积结果越大,投影越大,两个向量间夹角越小,方向越一致,相似度越高。
如果是垂直井,这个角度可以转换为相对于北向的方位角φr:
φr=φs+θ (3)
φs是垂直井时的震源方位角。对于垂直井,可以通过将θ添加到源-检波器方位角(φs)相对于正北的角度来实现,并且假设检波器的水平分量在x-y平面。因此,在解震源-检波器方位角时,只需使用震源位置的x和y坐标即可,如下所示:
此外,由于方程1只能产生±90°之间的角度,将存在两个潜在的检波器趋势,它们相隔180度;通过检查第一个初至点的极性,可以消除这种歧义。如果观测井存在任意偏差,那么在井的任何点,特别是在检波器位置,需要考虑一条与偏差相切的线l。利用球坐标,用参数化的方式表达为:
θw(井的倾角):井相对于垂直方向的倾斜角度。
φw(井的水平方向):井相对于正x轴的水平方向的角度。
xr,yr,zr(检波器的坐标):检波器在三维空间中的坐标,分别表示x、y和z坐标。
使用参数l的方向(在方程5的第一项中的向量)来定义一个垂直于井的平面。在这个平面上,引入了伪x、y和z轴。这些轴的选择是通过参数l的方向来确定的,这样就能够建立一个相对于井在某一点的坐标系。
为了进行地震检波器方向的分析,必须将震源坐标投影到上述定义的平面上。这可以通过计算由源坐标xs、ys和zs形成的向量与方程6中定义的伪x和y轴的内积来实现。然后通过在方程4中用投影震源坐标xs'和ys'替换实际源坐标,定义了一个伪源-检波器方位角φs';最后,将φs'加到θ上将得到检波器方向方位角,相对于伪y轴。请注意,方程6将在井为垂直井的情况下(即θw=0°,φw被选择为-90°)正确生成真实坐标。
二、使用三分量检波器监测的数据可以计算出震源的方位信息,下面使用分布式声波传感器计算震源的深度和距离。
与简单的波形直接叠加方法相比,本发明还可以检测较低信噪比事件,因为它基于信号一致性而不是信噪比本身,使用分布式声波传感器进行波形相似叠加,亮度函数即相似系数S(k):
对于特定的时窗,fij是第i道第j个采样点,M是参与计算的地震道数目,N为时窗大小,k为时窗中心位置。若所有的地震道数值相同,则S(k)=1。S(k)的取值范围是0到1,但通常等于1的情况太过理想,可以认为整个时空中的亮点函数取最大值时是地震发生的位置和事件发生的时刻。
使用分布式声波传感器监测的数据,在垂直剖面进行相似系数的计算,在某一特定的垂直剖面下建立震源扫描网格,混合阵列有两种布设方式,震源扫描网格如图4、图5所示;
震源扫描网格的基本思想是把空间划分成规则的网格,在震源扫描网格中,以每个划分的网格为基本单元,在基本单元网格中,根据检波器拾取的到时,以最早到的初至时间为标准,进行波形校正,把时差对每个波形进行整体偏移,调整时差之后,使用公式7的相似函数计算,在不同震源扫描网格中进行滑动网格计算,其相似函数取最大值时,网格中的位置是震源位置。如果是真实震源,那么到时最早的检波器应该是距离震源位置最近的,如图5所示,假设时差是3ms,那么使用3ms的时差去校正其他检波器接收到的信号,其计算波形相似性的值就会比较大。如果不是真实震源位置,那么在做时差校正后,波形在时间上是有到时差的,其波形相似性的值就会比较小,通过比较计算出的S(k)值,可以判断震源位置。
根据压裂设计方案,震源扫描网格包含压裂区域,已知压裂区域的坐标范围,震源扫描网格的坐标范围已知,根据网格划分的大小,计算出的震源网格的坐标也是已知的。根据三分量检波器计算的方位信息,最终三维空间中的震源位置为:
(x1,y1,z1)为某一特定垂直剖面中,震源位置坐标,(x2,y2,z2)是根据三分量检波器计算出的方位信息,计算得到的最终的震源坐标。
三、震源定位误差的分析方法:
震源定位测量的走时与理论走时之间的误差,计算方法如公式8;
定位误差是已知地层速度模型的情况下进行计算,计算方法如公式9:
NP为P波在分布式声波传感器阵列中到达的数量,NS是S波在分布式声波传感器阵列中到达的数量。是各道P波的监测实际走时与理论走时的误差,/>是各道S波的监测实际走时与理论走时的误差。
VP是地层中P波速度,VS是地层中S波速度。
上述根据本发明的方法可在硬件、固件中实现,或者被实现为可存储在记录介质(诸如CD ROM、RAM、软盘、硬盘或磁光盘)中的软件或计算机代码,或者被实现通过网络下载的原始存储在远程记录介质或非暂时机器可读介质中并将被存储在本地记录介质中的计算机代码,从而在此描述的方法可被存储在使用通用计算机、专用处理器或者可编程或专用硬件(诸如ASIC或FPGA)的记录介质上的这样的软件处理。可以理解,计算机、处理器、微处理器控制器或可编程硬件包括可存储或接收软件或计算机代码的存储组件(例如,RAM、ROM、闪存等),当所述软件或计算机代码被计算机、处理器或硬件访问且执行时,实现在此描述的混合阵列微震监测方法。此外,当通用计算机访问用于实现在此示出的处理的代码时,代码的执行将通用计算机转换为用于执行在此示出的处理的专用计算机。
本领域的普通技术人员将会意识到,这里所述的实施例是为了帮助读者理解本发明的实施方法,应被理解为本发明的保护范围并不局限于这样的特别陈述和实施例。本领域的普通技术人员可以根据本发明公开的这些技术启示做出各种不脱离本发明实质的其它各种具体变形和组合,这些变形和组合仍然在本发明的保护范围内。
Claims (10)
1.一种利用DAS和检波器混合阵列监测微震的系统,其特征在于,包括:分布式声波传感器和检波器;
混合阵列布设方式如下:
分布式声波传感器的光纤水平布设在地表上,光纤标距为3m,三分量检波器沿光纤等间隔布设在水平地表。
2.根据权利要求1所述的一种利用DAS和检波器混合阵列监测微震的系统,其特征在于:所述混合阵列布设方式如下:
分布式声波传感器的光纤垂直布设,沿垂直井方向,在光纤尾端布设3-5个三分量检波器,在井下的光纤,标距为3m。
3.一种利用DAS和检波器混合阵列监测微震的方法,其特征在于:基于权利要求1或2所述的一种利用DAS和检波器混合阵列监测微震的系统实现;
所述利用DAS和检波器混合阵列监测微震的方法包括以下步骤:
步骤一,通过三分量检波器中的两个水平分量X和Y进行震源方位角的计算;使用互相关算子计算X分量和震源之间的夹角θ,并将其转换为相对于北向的方位角φr;
通过检查第一个初至点的极性,消除可能存在的两个潜在的检波器趋势,以消除歧义;
在存在任意偏差的观测井中,通过考虑与偏差相切的线l来解决问题,这需要使用球坐标参数井的倾角θw和井的水平方向φw,以及检波器的坐标xr,yr,zr;
利用参数l的方向定义一个垂直于井的平面,并在该平面上引入伪x、y和z轴,建立一个相对于井在某一点的坐标系;
为了进行地震检波器方向的分析,必须将震源坐标投影到定义的平面上,通过计算源坐标和伪x和y轴的内积来实现,并定义了一个伪源-检波器方位角φs';最后,将φs'加到θ上将得到检波器方向方位角,相对于伪y轴;
步骤二,使用分布式声波传感器监测数据,在垂直剖面上进行相似系数S(k)的计算,并建立震源扫描网格;震源扫描网格的基本思想是将空间划分成规则的网格,根据检测到的到时进行波形校正,并通过相似系数S(k)来确定震源位置;
根据三分量检波器计算的方位信息,得到最终的震源坐标。
4.根据权利要求3所述的一种利用DAS和检波器混合阵列监测微震的方法,其特征在于:震源方位角的计算如下式:
是零滞后的互相关算子,使用互相关算子计算X分量和震源之间的夹角θ如式2:
5.根据权利要求4所述的一种利用DAS和检波器混合阵列监测微震的方法,其特征在于:所述转换为相对于北向的方位角φr,公式如下:
φr=φs+θ (3)
φs是垂直井时的震源方位角;对于垂直井,通过将θ添加到源-检波器方位角φs相对于正北的角度来实现,并且假设检波器的水平分量在x-y平面;因此,在解震源-检波器方位角时,需使用震源位置的x和y坐标,如下所示:
6.根据权利要求5所述的一种利用DAS和检波器混合阵列监测微震的方法,其特征在于:所述偏差相切的线l,利用球坐标,用参数化的方式表达为:
式中,θw为井相对于垂直方向的倾斜角度;
φw为井相对于正x轴的水平方向的角度;
xr,yr,zr为检波器在三维空间中的坐标,分别表示x、y和z坐标;
建立一个相对于井在某一点的坐标系,具体如下:
7.根据权利要求6所述的一种利用DAS和检波器混合阵列监测微震的方法,其特征在于:所述计算相似系数S(k)的公式如下:
式中,fij是第i道第j个采样点,M是参与计算的地震道数目,N为时窗大小,k为时窗中心位置。
8.根据权利要求7所述的一种利用DAS和检波器混合阵列监测微震的方法,其特征在于:所述震源扫描网格的基本思想是把空间划分成规则的网格,在震源扫描网格中,以每个划分的网格为基本单元,在基本单元网格中,根据检波器拾取的到时,以最早到的初至时间为标准,进行波形校正,把时差对每个波形进行整体偏移,调整时差之后,使用公式7的相似函数计算,在不同震源扫描网格中进行滑动网格计算,其相似函数取最大值时,网格中的位置是震源位置;通过比较计算出的S(k)值,可以判断震源位置。
9.根据权利要求8所述的一种利用DAS和检波器混合阵列监测微震的方法,其特征在于:计算得到最终的震源坐标,公式如下:
式中,(x1,y1,z1)为某一特定垂直剖面中,震源位置坐标,(x2,y2,z2)是根据三分量检波器计算出的方位信息,计算得到的最终的震源坐标。
10.根据权利要求9所述的一种利用DAS和检波器混合阵列监测微震的方法,其特征在于:在步骤二之后,计算出震源定位的均方根误差ΔTRMS和速度加权的均方根误差ERMS,评估定位结果震源位置的准确性和可靠性,具体计算公式如下:
震源定位测量的走时与理论走时之间的误差,计算方法如公式8;
定位误差是已知地层速度模型的情况下进行计算,计算方法如公式9,
式中,NP为P波在分布式声波传感器阵列中到达的数量,NS是S波在分布式声波传感器阵列中到达的数量;ΔTi P是各道P波的监测实际走时与理论走时的误差,ΔTi S是各道S波的监测实际走时与理论走时的误差;
式中,VP是地层中P波速度,VS是地层中S波速度。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202311650833.8A CN117687077B (zh) | 2023-12-05 | 2023-12-05 | 利用das和检波器混合阵列监测微震的方法和系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202311650833.8A CN117687077B (zh) | 2023-12-05 | 2023-12-05 | 利用das和检波器混合阵列监测微震的方法和系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN117687077A true CN117687077A (zh) | 2024-03-12 |
CN117687077B CN117687077B (zh) | 2024-05-03 |
Family
ID=90134378
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202311650833.8A Active CN117687077B (zh) | 2023-12-05 | 2023-12-05 | 利用das和检波器混合阵列监测微震的方法和系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN117687077B (zh) |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20150135819A1 (en) * | 2012-06-11 | 2015-05-21 | Kobold Services Inc. | Microseismic monitoring with fiber-optic noise mapping |
CN113484912A (zh) * | 2021-08-02 | 2021-10-08 | 中油奥博(成都)科技有限公司 | 页岩油气光纤智能地球物理数据采集系统及采集方法 |
CN113775330A (zh) * | 2021-09-23 | 2021-12-10 | 中油奥博(成都)科技有限公司 | 高温地热田综合地球物理勘探系统及地热甜点区评价方法 |
CN113791443A (zh) * | 2021-09-13 | 2021-12-14 | 中油奥博(成都)科技有限公司 | 分布式光纤地震数据采集系统及其数据收集转换方法 |
CN115079250A (zh) * | 2022-01-18 | 2022-09-20 | 中油奥博(成都)科技有限公司 | 基于光纤传感技术的co2封存地选址、安全监测系统及方法 |
WO2022203945A1 (en) * | 2021-03-24 | 2022-09-29 | Schlumberger Technology Corporation | True particle velocity wavefield processing in fiber optics – particle motion sensor hybrid array |
CN115685339A (zh) * | 2022-09-23 | 2023-02-03 | 中油奥博(成都)科技有限公司 | 三分量分布式光纤声波传感阵列声波测井装置及测量方法 |
-
2023
- 2023-12-05 CN CN202311650833.8A patent/CN117687077B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20150135819A1 (en) * | 2012-06-11 | 2015-05-21 | Kobold Services Inc. | Microseismic monitoring with fiber-optic noise mapping |
WO2022203945A1 (en) * | 2021-03-24 | 2022-09-29 | Schlumberger Technology Corporation | True particle velocity wavefield processing in fiber optics – particle motion sensor hybrid array |
CN113484912A (zh) * | 2021-08-02 | 2021-10-08 | 中油奥博(成都)科技有限公司 | 页岩油气光纤智能地球物理数据采集系统及采集方法 |
CN113791443A (zh) * | 2021-09-13 | 2021-12-14 | 中油奥博(成都)科技有限公司 | 分布式光纤地震数据采集系统及其数据收集转换方法 |
CN113775330A (zh) * | 2021-09-23 | 2021-12-10 | 中油奥博(成都)科技有限公司 | 高温地热田综合地球物理勘探系统及地热甜点区评价方法 |
CN115079250A (zh) * | 2022-01-18 | 2022-09-20 | 中油奥博(成都)科技有限公司 | 基于光纤传感技术的co2封存地选址、安全监测系统及方法 |
CN115685339A (zh) * | 2022-09-23 | 2023-02-03 | 中油奥博(成都)科技有限公司 | 三分量分布式光纤声波传感阵列声波测井装置及测量方法 |
Non-Patent Citations (1)
Title |
---|
ZHENDONG ZHANG 等: "Characterizing Microearthquakes Induced by Hydraulic Fracturing With Hybrid Borehole DAS and Three-Component Geophone Data", IEEE TRANSACTIONS ON GEOSCIENCE AND REMOTE SENSING, vol. 61, 5 April 2023 (2023-04-05), pages 1 - 15 * |
Also Published As
Publication number | Publication date |
---|---|
CN117687077B (zh) | 2024-05-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
RU2494418C1 (ru) | Способ измерения координат микросейсмических источников и параметров механизмов их очагов в условиях сильных сейсмических помех (варианты) | |
Liu et al. | A new 3D observation system designed for a seismic ahead prospecting method in tunneling | |
CN105510880A (zh) | 一种基于双差法的微地震震源定位方法 | |
CN112213768B (zh) | 一种联合震源机制反演的地面微地震定位方法及系统 | |
CN105676279A (zh) | 一种采用同心圆等炮检距的地震反射数据采集方法 | |
CN112305591B (zh) | 隧道超前地质预报方法、计算机可读存储介质 | |
CN110632667B (zh) | 一种基于冲击波激震条件下的隐伏陷落柱超前探测方法 | |
Zhang et al. | Automated microseismic event location by amplitude stacking and semblance | |
CN111103620A (zh) | 一种岩巷超前探测三维偏移成像方法 | |
CN108693560B (zh) | 一种基于互相关道的散射波成像方法及系统 | |
CN111352153B (zh) | 一种基于瞬时相位互相关加权的微地震干涉定位方法 | |
CN110579799B (zh) | 一种等旅行时间间隔的地震采集观测方法及系统 | |
CN110967751B (zh) | 基于地面浅井监测的微地震事件的定位方法及存储介质 | |
CN117687077B (zh) | 利用das和检波器混合阵列监测微震的方法和系统 | |
US20170299743A1 (en) | Estimating an earth response | |
CN110780346A (zh) | 一种隧道超前探测复杂地震波场的分离方法 | |
CN107664771B (zh) | 一种基于相似性系数的微地震全波形定位方法 | |
CN113568041B (zh) | 时移地震三维拖缆采集数据的可重复性分析方法及系统 | |
CN115980851A (zh) | 复合震源参数反演的方法、计算机设备及可读存储介质 | |
CN112904348B (zh) | 一种三维探测方法、装置、设备和存储介质 | |
CN110579798B (zh) | 一种等反射角度间隔的地震采集观测方法及系统 | |
CN109375261B (zh) | 用于地表微地震监测的观测台站的传感器布局方法及系统 | |
CN110531416B (zh) | 一种基于三分量反射信号时频域极化参数的断层确定方法 | |
CN115184987B (zh) | 地层倾角信息的提取方法、装置、服务器及存储介质 | |
CN112444869B (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 |