CN113050159B - 一种煤岩水力压裂裂缝微震定位及扩展机理监测方法 - Google Patents
一种煤岩水力压裂裂缝微震定位及扩展机理监测方法 Download PDFInfo
- Publication number
- CN113050159B CN113050159B CN202110307395.XA CN202110307395A CN113050159B CN 113050159 B CN113050159 B CN 113050159B CN 202110307395 A CN202110307395 A CN 202110307395A CN 113050159 B CN113050159 B CN 113050159B
- Authority
- CN
- China
- Prior art keywords
- microseismic
- waveform
- event
- hydraulic fracturing
- fracture
- 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
- 239000011435 rock Substances 0.000 title claims abstract description 71
- 239000003245 coal Substances 0.000 title claims abstract description 66
- 238000000034 method Methods 0.000 title claims abstract description 66
- 238000012544 monitoring process Methods 0.000 title claims abstract description 41
- 230000007246 mechanism Effects 0.000 title claims abstract description 31
- 238000012216 screening Methods 0.000 claims abstract description 70
- 230000033001 locomotion Effects 0.000 claims abstract description 15
- 230000008569 process Effects 0.000 claims abstract description 15
- 230000035945 sensitivity Effects 0.000 claims abstract description 8
- 230000006870 function Effects 0.000 claims description 27
- 238000001914 filtration Methods 0.000 claims description 17
- 230000006378 damage Effects 0.000 claims description 13
- 125000004432 carbon atom Chemical group C* 0.000 claims description 12
- 238000005553 drilling Methods 0.000 claims description 9
- 238000002347 injection Methods 0.000 claims description 5
- 239000007924 injection Substances 0.000 claims description 5
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 5
- 238000004364 calculation method Methods 0.000 claims description 4
- 238000010276 construction Methods 0.000 claims description 4
- 230000035515 penetration Effects 0.000 claims description 4
- 238000012935 Averaging Methods 0.000 claims description 3
- 238000002474 experimental method Methods 0.000 claims description 3
- 230000000630 rising effect Effects 0.000 claims description 3
- 238000005070 sampling Methods 0.000 claims description 3
- 230000009466 transformation Effects 0.000 claims description 3
- 230000005484 gravity Effects 0.000 claims description 2
- 238000009434 installation Methods 0.000 claims description 2
- 238000011156 evaluation Methods 0.000 abstract description 4
- 239000007789 gas Substances 0.000 description 14
- 238000005065 mining Methods 0.000 description 13
- 238000010586 diagram Methods 0.000 description 10
- 230000000694 effects Effects 0.000 description 10
- 238000005422 blasting Methods 0.000 description 5
- 238000005516 engineering process Methods 0.000 description 5
- 239000012530 fluid Substances 0.000 description 5
- 238000000605 extraction Methods 0.000 description 4
- 230000035699 permeability Effects 0.000 description 4
- 238000011161 development Methods 0.000 description 3
- 230000006698 induction Effects 0.000 description 3
- 238000012986 modification Methods 0.000 description 3
- 230000004048 modification Effects 0.000 description 3
- 239000000243 solution Substances 0.000 description 3
- 230000009286 beneficial effect Effects 0.000 description 2
- 230000008878 coupling Effects 0.000 description 2
- 238000010168 coupling process Methods 0.000 description 2
- 238000005859 coupling reaction Methods 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 230000008595 infiltration Effects 0.000 description 2
- 238000001764 infiltration Methods 0.000 description 2
- 230000000977 initiatory effect Effects 0.000 description 2
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 230000009471 action Effects 0.000 description 1
- 238000009412 basement excavation Methods 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 239000007788 liquid Substances 0.000 description 1
- 230000004807 localization Effects 0.000 description 1
- 238000012806 monitoring device Methods 0.000 description 1
- 239000003345 natural gas Substances 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
- 239000003209 petroleum derivative Substances 0.000 description 1
- 230000002265 prevention Effects 0.000 description 1
- 238000011084 recovery Methods 0.000 description 1
- 230000005641 tunneling Effects 0.000 description 1
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/288—Event detection in seismic signals, e.g. microseismics
-
- 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
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- Acoustics & Sound (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
本发明公开了一种煤岩水力压裂裂缝微震定位及扩展机理监测方法,首先根据实际情况布置微震传感器,采集压裂过程中的原始波形信号;在此基础上提出可更新长度的事件筛选窗口,利用长短时窗法扫描各通道数据,识别有效微震事件;通过定义波形敏感值筛选出煤岩体水力压裂诱发的有效微震事件,然后利用局部AIC法获得微震事件的高精度到时,并计算出破裂源的位置,同时根据P波幅值及初动信息反演震源机制;最后将微震定位和震源机制信息结合起来,可以揭示水力压裂裂缝的空间分布形态和破裂机制,确定煤层水力压裂影响范围,这对于压裂过程中参数实时动态调整、压裂半径确定、储层评价具有指导性意义。
Description
技术领域
本发明涉及一种煤岩体压裂监测评价方法,具体涉及一种煤岩水力压裂裂缝微震定位及扩展机理监测方法,属于岩体破裂及工程地球物理监测技术领域。
背景技术
随着煤炭资源开采深度和强度的增加,瓦斯灾害,特别是煤与瓦斯耦合的动力灾害愈来愈严重;而瓦斯(煤层气)虽是井下重大危险源,但它也是一种宝贵的清洁能源;对低透气性煤层进行卸压增透,提高储层的渗透率和抽采瓦斯率是深部瓦斯灾害治理与瓦斯资源高效开发利用的重要途径;水力压裂是石油、页岩气等油气资源开发的核心技术,已经广泛应用于低透气性储层渗透率的改造;近年来,水力压裂技术在煤矿井下煤层卸压增透防治瓦斯灾害及提高瓦斯抽采效率等方面得到了较为成功的应用;然而煤矿井下煤岩体水力压裂裂缝监测及压裂效果评价方面一直没有有效的技术手段,严重制约着煤岩体水力压裂理论、技术及应用进一步发展;微震定位监测技术通过采集煤岩体破裂产生的弹性波信号,根据震源定位结果对煤岩体破坏进行监测评价;目前,微震监测技术主要用来对矿井冲击地压、岩爆等动力灾害进行监测预警;作为一种能够对煤岩体破裂进行有效监测的手段,在煤岩体水力压裂监测评价方面的研究和应用较少;煤矿井下煤岩体水力压裂作业的同时,在压裂区域附近同时存在着煤层工作面回采、巷道掘进、***等采掘活动,水力压裂及上述采掘活动都会造成煤岩体损伤破裂,从而诱发微震波形信号;此外,在实验室煤岩试样水力压裂微震监测实验研究过程中,煤岩试样在受水压力作用破坏的同时,其外加载荷也会对煤岩试样造成破坏,产生微震波形信号;若不加区别的直接利用监测的微震波形信号对水力压裂裂缝形态和扩展过程进行监测,可能会造成与真实的水力压裂裂缝存在较大偏差;因此,如何区分水力压裂诱发和采掘等其它活动产生的微震信号,对于煤岩体水力压裂裂缝微震定位监测结果真实可靠性至关重要;近年来,一些学者和现场工程师也开始关注煤岩体水力压裂微震监测,提出了一些水力压裂微震监测装置与方法,然而现有的技术并没有涉及到上述提到的问题,也没有考虑到精细化处理微震的每一个环节;实际上实验室与现场传感器布置位置、有效微震事件的准确识别、高精度到时信息与震源定位结果都会影响煤岩体水力压裂裂缝微震真实可靠性;更为重要的是现有的各类压裂监测方法主要采用微震震源定位结果对水力压裂进行监测,较少涉及同时采用震源定位及震源机制反演结果联合表征水力压裂裂缝时空分布及扩展机理。目前这些问题仍然存在,并没被很好的解决。
目前迫切需要一种煤岩水力压裂裂缝微震定位及扩展机理监测方法。
发明内容
针对上述存在的技术不足,本发明提供一种煤岩水力压裂裂缝微震定位及扩展机理监测方法。
本发明采用如下技术方案:
本发明提供一种煤岩水力压裂裂缝微震定位及扩展机理监测方法,具体包括以下步骤:
(a)根据煤岩体压裂的作业地点确定需要监测的重点区域;根据井下巷道的几何形状、施工现场环境和岩体的性质,布置X个微震传感器,记做Ci(i=1、2、3...X),所布置的微震传感器能够监测整个压裂区域;
(c)对微震传感器采集的原始波形信号进行滤波预处理,滤波的频率范围是由水力压裂诱发微震信号的频率来确定;所述原始波形信号滤波结束后,基于背景噪音振幅的平均值对微震传感器Ci(i=1、2、3...X)采集的波形数据进行扫描,当任意微震传感器Cj采集波形信号的幅值超过时,将该微震传感器波形信号幅值超过的采样点记为起跳点Sj,起跳点所对应的时间点为Tj,然后从起跳点Sj前置一段时间Tb,后置一段时间Ta,将时间区间[Tb,Ta]定义为事件筛选窗口Wj;
(d)采用长短时窗法对事件筛选窗口Wj(j=1、2、3...)中所有微震传感器采集的波形信号进行扫描,若微震传感器Ci(i=1、2、3...X)中波形信号的特征函数值大于触发阈值λ时,判定微震传感器Ci内存在有效微震波形;当事件筛选窗口Wj内的有效微震波形数B≥m且m≥4时,确定该事件筛选窗口中存在有效微震事件,同时存储该有效微震事件的波形数据,否则,丢弃该事件筛选窗口内的波形数据;
(e)顺序扫描时刻Tj+Tb之后的波形信号,统计得到所有的有效微震事件数为E;
(f)利用波形敏感值对E个有效微震事件进行筛选;
上式中,表示第j个事件筛选窗口内第i个微震传感器采集的有效微震波形数,代表第j个事件筛选窗口内第i个微震传感器采集到的第k个有效微震波形的频率;当波形平均频率时,判定为优选信号,输出结果为1,反之输出结果为0;并对最终输出结果赋予权重w1;
上式中,代表第j个事件筛选窗口内第i个微震传感器采集到的第k个有效微震波形最大幅值和起跳点幅值的差值,表示第j个事件筛选窗口内第i个微震传感器采集到的第k个有效微震波形的上升时间;当平均幅值增长速率介于m3~m4,判断为优选信号,输出结果为1,反之输出结果为0,并对最终输出结果赋予权重w3;
上式中,表示第j个事件筛选窗口内第i个微震传感器采集的第k个有效微震波形的衰减时间,表示第j个事件筛选窗口内第i个微震传感器采集的第k个有效微震波形的持续时间;当平均衰减比介于m5~m6,判断为优选信号,输出结果为1,反之输出结果为0;并对最终输出结果赋予权重w4;
(5)计算波形敏感参数总权重,波形敏感参数总权重=(1|0)*w1+(1|0)*w2+(1|0)*w3+(1|0)*w4,当波形敏感参数总权重>w,判定这个事件为煤岩体水力压裂诱发产生的有效微震事件,按照时间顺序,依次对E个有效微震事件进行筛选,最终得到N个煤岩体水力压裂诱发产生的有效微震事件;
(g)利用局部赤迟信息准则法(AIC)拾取煤岩体水力压裂诱发微震事件的高精度P波到时、P波初至振幅及初动信息;
(h)在已知高精度P波到时的基础上,构建目标函数;引入迭代法,在三维空间坐标系中构建一个四面体,并以四面体的四个顶点作为破裂源,计算目标函数的最大值和最小值,所述最大值和最小值计算完毕后进行比较,抛弃使得目标函数取得最大值的点,同时按照变换规则寻求新的顶点构建新的四面体;重复上述步骤直至求得目标函数的最小值满足迭代结束的条件,即为定位事件点;
(i)在微震事件高精度定位的基础上,利用微震传感器采集的P波初至振幅信息结合定位事件点的坐标、微震传感器坐标反演震源矩张量,并在主轴坐标系下将震源矩张量分解为同性部分ISO、纯剪破坏部分DC、其他部分CLVD,同时根据P波初动信息绘出震源沙滩球,表示定位事件点的空间方位信息(走向、倾向、断层夹角);
(j)将定位事件点的坐标和空间方位信息结合起来,准确绘出裂缝的空间形态,分析裂缝的延伸与贯通,判断起裂时的破坏类型,裂缝延伸时的破坏类型,延伸的裂缝遇到天然节面和层理时的破坏类型。
优选地,所述步骤(a)中微震传感器放置的数量至少为5个,对于室内水力压裂实验,微震传感器应至少布置在样品的两个表面上;对于现场煤岩体水力压裂,微震传感器应采用钻孔安装方式布置在稳定性较好的岩层中。
优选地,所述步骤(c)中,对微震传感器采集的原始波形信号进行滤波预处理时,采用带通滤波器,滤波器的带宽为f1~f2Hz,仅允许高于较低的频率f1且低于较高的频率f2的特定频段的波形信号通过。
优选地,所述步骤(d)中,波形信号的特征函数值依照如下公式计算:
式中:STA代表短时窗、LTA代表长时窗,nl、ns分别为长时窗长度、短时窗长度,Xi为第i个采样点对应的波形信号的幅值。
优选地,所述步骤(f)中,波形敏感值包括波形平均频率、波形偏离值、平均幅值增长速率、平均衰减比。
优选地,所述步骤(h)中,所述目标函数是利用每个微震传感器采集走时与理论走时之差的绝对值来获得,具体采用以下公式:
上式中,f(x)为目标函数,t为发震时刻,ti为第i个微震传感器的采集到时,υp为P波波速,经过多次迭代,当目标函数的最小值小于预设的范围,迭代结束。
优选地,所述步骤(i)、(j)中,所述P波初至振幅信息反演水力压裂裂缝的破裂类型,需已知6个或6个以上微震传感器采集到的P波初至振幅,有效微震事件的破裂类型是根据纯剪切破坏部分MDC占震源矩张量比重PDC来确定,其中,PDC>60为剪切破裂,40<PDC<60为混合破裂,PDC>40为张拉破裂。所述P波初动信息是指P波从定位事件点传至地面时其初始震动的方向,即初动方向;将每个传感器采集的P波作适当的校正之后,把初动方向标到一个以震源为中心,且半径很小的球面上,即震源沙滩球上。当记录的数据足够多且传感器数量较多时,总可以找到两个互相垂直的节面,把震源沙滩球4等分;根据震源沙滩球上两个节面的方位就可以确定定位事件点的空间方位信息。
本发明的有益效果在于:
1、本发明考虑到现场压裂与实验室压裂微震传感器布置的区别,针对实际情况提出合理的传感器布置方式,提高波形信息采集完整性。
2、与传统的微震事件检测相比,本发明提出可更新长度的事件筛选窗口,避免固定时间窗口微震事件拾取精度低、事件丢失率高的问题;同时由于事件筛选窗口的存在,改善了赤迟信息准则法(AIC)在低信噪比的情况下拾取P波到时不唯一的情况,帮助相关从业人员获得高精度的P波到时。
3、本发明利用煤岩体压裂诱发的微震信号具有相对低频、持续时间偏长、起跳不明显等独特的特征,根据波形敏感参数(波形平均频率、波形偏离值、平均幅值增长速率、平均衰减比),可以区别水力压裂诱发和采掘等其它活动产生的微震信号,极大提高了煤岩体水力压裂裂缝微震定位监测结果真实可靠性,同时也降低时间成本与计算工作量,为后续震源机制反演提供便利。
4、引入了迭代法对震源的位置进行求解,避免通过增加微震传感器个数求解微震位置所带来的成本,提高了经济效益,与传统的水力压裂定位监测方法相比,引入迭代法可更为精确的求解震源的位置,由于迭代算法总是寻找目标函数的最优解,故求解的破裂源位置唯一,不会出现不收敛的情况。
5、本发明将微震定位结果与震源机制参数信息结合起来,能够表征煤岩体压裂裂缝扩展形态、空间分布特征和破裂机理。这对确定破裂半径、破裂面积、水力压裂模型的建立、压裂参数动态调整、现场指导等方面具有重要意义,同时也能够为石油气、天然气、页岩气等致密油气开采监测提供参考和对比,此外也可以拓展到地铁、井下巷道、桥隧开挖等岩体工程造成的裂缝评价与监测。
6、本发明是完整的一套煤岩水力压裂裂缝微震定位及扩展机理监测方法,包括微震传感器的布置、滤波处理、微震事件识别、高精度到时获取、震源定位与震源机制反演,这可以为相关从业人员进行微震数据处理提供经验借鉴。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本发明方法流程图;
图2为本发明含有背景噪音的原始波形信号示意图;
图3(a)为煤岩体水力压裂过程中采集到的第j(j=3)个事件筛选窗,第i(i=9)个微震传感器采集的有效微震波形与事件筛选窗口示意图;
图3(b)为煤岩体水力压裂过程中采集到的第j(j=3)个事件筛选窗口,所有微震传感器采集的有效微震波形与事件筛选窗口示意图;
图4波形特征参数示意图;
图5为煤岩体穿层钻孔水力压裂示意图;
图6为煤岩体水力压裂裂缝微震定位及扩展机理表征图;
图中:1-有效微震波形;2-长时窗;3-短时窗;4-起跳点;5-事件筛选窗口;6-曲率半径;7-包络线;8-煤层;9-泥岩;10-泥岩与砂岩;11-巷道;12-水力压力钻孔;13-封孔器;14-微震传感器;15-定位事件点;16-水力压裂裂缝;17-天然弱面或节理;18-震源沙滩球;19-破裂半径。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
如图1-6所示,本实例提供一种煤岩水力压裂裂缝微震定位及扩展机理监测方法,具体包括以下步骤:
(a)根据煤岩体压裂的作业地点确定需要监测的重点区域;根据井下巷道的几何形状、施工现场环境和岩体的性质,布置X个微震传感器,记做Ci(i=1、2、3...X),所布置的微震传感器能够监测整个压裂区域;
(c)对微震传感器采集的原始波形信号进行滤波预处理,滤波的频率范围是由水力压裂诱发微震信号的频率来确定;所述原始波形信号滤波结束后,基于背景噪音振幅的平均值对微震传感器Ci(i=1、2、3...X)采集的波形数据进行扫描,当任意微震传感器Cj采集波形信号的幅值超过时,将该微震传感器波形信号幅值超过的采样点记为起跳点Sj,起跳点所对应的时间点为Tj,然后从起跳点Sj前置一段时间Tb,后置一段时间Ta,将时间区间[Tb,Ta]定义为事件筛选窗口Wj;
(d)采用长短时窗法对事件筛选窗口Wj(j=1、2、3...)中所有微震传感器采集的波形信号进行扫描,若微震传感器Ci(i=1、2、3...X)内波形信号的特征函数值大于触发阈值λ时,判定微震传感器Ci内存在有效微震波形;当事件筛选窗口Wj内的有效微震波形数B≥m且m≥4时,确定该事件筛选窗口中存在有效微震事件,同时存储该有效微震事件的波形数据,否则,丢弃该事件筛选窗口内的波形数据;
(e)顺序扫描时刻Tj+Tb之后的波形信号,统计得到所有的有效微震事件数为E;
(f)利用波形敏感值对E个有效微震事件进行筛选;
上式中,表示第j个事件筛选窗口内第i个微震传感器采集的有效微震波形数,代表第j个事件筛选窗口内第i个微震传感器采集到的第k个有效微震波形的频率;当波形平均频率时,判定为优选信号,输出结果为1,反之输出结果为0;并对最终输出结果赋予权重w1;
上式中,代表第j个事件筛选窗口内第i个微震传感器采集到的第k个有效微震波形最大幅值和起跳点幅值的差值,表示第j个事件筛选窗口内第i个微震传感器采集到的第k个有效微震波形的上升时间;当平均幅值增长速率介于m3~m4,判断为优选信号,输出结果为1,反之输出结果为0,并对最终输出结果赋予权重w3;
上式中,表示第j个事件筛选窗口内第i个微震传感器采集的第k个有效微震波形的衰减时间,表示第j个事件筛选窗口内第i个微震传感器采集的第k个有效微震波形的持续时间;当平均衰减比介于m5~m6,判断为优选信号,输出结果为1,反之输出结果为0;并对最终输出结果赋予权重w4;
(5)计算波形敏感参数总权重,波形敏感参数总权重=(1|0)*w1+(1|0)*w2+(1|0)*w3+(1|0)*w4,当波形敏感参数总权重>w,判定这个事件为煤岩体水力压裂诱发产生的有效微震事件,按照时间顺序,依次对E个有效微震事件进行筛选,最终得到N个煤岩体水力压裂诱发产生的有效微震事件;
(g)利用局部赤迟信息准则法(AIC)拾取煤岩体水力压裂诱发微震事件的高精度P波到时、P波初至振幅及初动信息;
(h)在已知高精度P波到时的基础上,构建目标函数;引入迭代法,在三维空间坐标系中构建一个四面体,并以四面体的四个顶点作为破裂源,计算目标函数的最大值和最小值,所述最大值和最小值值计算完毕后进行比较,抛弃使得目标函数取得最大值的点,同时按照变换规则寻求新的顶点构建新的四面体;重复上述步骤直至目标函数的最小值满足迭代结束的条件,即为定位事件点;
(i)在微震事件高精度定位的基础上,利用采集的P波初至振幅信息结合定位事件点的坐标、微震传感器坐标反演震源矩张量,并在主轴坐标系下将震源矩张量分解为同性部分ISO、纯剪切破坏部分DC、其他部分CLVD,同时根据P波初动信息绘出震源沙滩球,表示破裂点的空间方位信息(走向、倾向);
(j)将定位事件点的坐标和空间方位信息结合起来,准确绘出裂缝的空间形态,分析裂缝的延伸与贯通,判断起裂时的破坏类型,裂缝延伸时的破坏类型,延伸的裂缝遇到天然节面和层理时的破坏类型。
如图3(a)、图3(b)是煤岩体底板穿层钻孔水力压裂过程中监测到的第j(j=3)个事件筛选窗口,第i(i=9)个微震传感器采集的有效微震波形。其中samples STA=20ms,samples LTA=50ms,采样率100HZ/s。由图3(a)和(b)可知波形起跳点对应的时间是1360ms,事件筛选窗口为[1210ms,1560ms],前置时间为150ms,后置时间为200ms;从图3(b)可知第j(j=3)个事件筛选窗口共采集到8个有效微震波形(1#~8#)。
图4是波形特征参数示意图。由图可知该有效微震波形的最大振幅为Amax,该有效微震波形起跳点对应振幅为A1,该有效微震波形的上升时间为trt,该有效微震波形的峰后衰减时间为ta,该有效微震波形的持续时间为ts,有效波形包络线曲率半径为ρ。
图6是煤岩体水力压裂裂缝微震定位及扩展机理表征图,由图可知水力压裂裂缝分布在以钻孔为轴线,破裂半径为400m的椭圆体范围内。
根据井下巷道的几何形状、施工现场环境以及岩体的性质,本次煤矿井下煤层底板穿层钻孔水力压裂微震监测共布置12个微震传感器14,记做Ci(i=1、2...12),如图5所示。
在开启高压注水泵的同时打开微震监测***,实时不间断的采集水力压裂全过程的原始波形信号,同时打开高压液体泵和微震监测***是为了保证信号不丢失,只要开始压裂作业就会有信号产生。煤岩体水力压裂产生的微震信号不同于煤层工作面回采、巷道掘进、***等采掘活动诱发的微震信号。煤岩体水力压裂微震信号是流体渗入与岩石破裂共同耦合产生,而上述采掘活动诱发的微震信号不涉及流体渗入;
原始波形信号滤波结束后,基于背景噪音振幅的平均值对微震传感器Ci(i=1、2、3...12)14采集的波形数据进行扫描,如图3(a)和(b)所示,扫描时发现微震传感器C914采集波形信号的幅值超过因此将该微震传感器14波形信号幅值超过的采样点记为起跳点S14,该起跳点S14对应时间为1360ms,然后从起跳点S14前置时间Tb=150ms,后置时间Ta=200ms,将时间区间[1210ms,1560ms]定义为事件筛选窗口W35;接着采用长短时窗法对事件筛选窗口W3(j=1、2、3...)5中所有微震传感器14的波形信号进行扫描,发现有8个有效微震波形1的波形特征函数值大于触发阈值λ=1.4,故第j(j=3)个事件筛选窗中存在一个有效微震事件。波形特征函数值按照以下公式计算:
按照上述方法扫描所有微震传感器14采集的波形信号,统计得到E=6988个有效微震事件。
为了能够区分煤岩体水力压裂诱发和煤层工作面回采、巷道掘进、***等采掘活动产生的有效微震事件,利用波形敏感值(波形平均频率、波形偏离值、平均幅值增长速率、平均衰减比)对E=6988个有效微震事件进行筛选;
根据以下公式计算得到第Wj个事件筛选窗口5内所有微震传感器14采集的有效微震波形1的波形平均频率
根据波形敏感参数(波形平均频率、波形偏离值、平均幅值增长速率、平均衰减比)可求出第j(j=3)个事件筛选窗5,有效微震事件的波形敏感参数总权重为:1*0.30+1*0.28+0*0.25+0*0.18=0.58,由于该波形敏感参数总权重>0.55,所以判定第j(j=3)个事件筛选窗5采集的有效微震事件是煤岩体底板穿层钻孔水力压裂诱发产生的有效微震事件而并非煤层工作面回采、巷道掘进、***等采掘活动诱发的产生的有效微震事件。
如图3(b),利用上述方法对第j(j=5)个事件筛选窗口5进行扫描并利用波形敏感参数进行判断,第j(j=5)个事件筛选窗口5内有效微震事件的波形敏感参数总权重为:0*0.30+0*0.28+1*0.25+1*0.18=0.43,由于该波形敏感参数总权重<0.55,所以判定该有效微震事件是煤层工作面回采、巷道掘进、***等采掘活动诱发产生。
接着利用波形敏感参数对所有微震传感器14采集的波形数据进行判别,最终得到4870个煤岩体底板穿层钻孔水力压裂诱发产生的有效微震事件;并采用局部AIC法拾取这些微震事件的高精度P波到时、P波初至振幅及初动信息。
在已知高精度P波到时的基础上,构建目标函数。引入迭代法,反复迭代,求得2178个定位事件点15,如图5所示。
由图6可知,根据定位事件点15可以粗略的绘出水力压裂裂缝的空间分布形态,判断出压裂的半径为400m,即水力压裂裂缝16分布在以水力压裂钻孔12为轴线,破裂半径19为400m的椭圆体范围内,仅利用定位事件点15分布情况来评价水力压裂的过程不具有说服力。震源机制蕴含着丰富的空间参数信息,因此利用P波初至振幅对2178个定位事件15进行震源矩张量反演,并把震源矩张量分解为同性部分ISO、纯剪切破坏部分DC、其他部分CLVD,统计剪切破坏部分占震源矩张量的比重,同时用P波初动信息绘出震源沙滩球18,从震源沙滩球18中可以获取定位事件点15的空间方位信息(走向、倾向、断层滑移夹角)。
如图6所示,将定位事件点15和震源沙滩球18所表示的空间方位信息结合起来,能够准确的绘出水力压裂裂缝16的空间形态;在煤岩体底板穿层钻孔水力压裂过程中,水力压裂裂缝16的扩展机制以剪切型占优,其中剪切破裂源1787个,占总事件数的81.6%;张拉破裂源266个,占事件总数的12.67%;混合破裂源125个,占事件总数的5.73%。由图6可知,裂缝起裂时发生的破裂类型为剪切破裂,这是由于压裂液的持续注入,在封孔器13中形成憋压环境,当封孔器13内部压裂液的压力大于周围应力环境时,就会将煤层8撑开裂缝;裂缝延伸时发生的破坏类型主要为剪切破裂但同时也出现了混合破裂与张拉破裂,出现混合破裂可能由于裂缝延伸时遇到了天然的弱面或节理17,延伸的裂缝直接贯穿弱面,产生混和破裂;在天然的弱面或节理17附近,定位事件点15颇多,这是因为天然的弱面或节理17本身就存在很多的微小裂纹,在压裂液作用下,更容易受到影响而贯通。
将震源定位和震源机制信息结合起来分析水力压裂裂缝延伸和扩展机理,这对确定破裂半径、破裂面积、水力压裂模型的建立、压裂参数动态调整、现场指导等方面具有重要意义,有益于压裂工程师评价水力压裂的过程。
显然,本领域的技术人员可以对本发明进行各种改动和变型而不脱离本发明的精神和范围。这样,倘若本发明的这些修改和变型属于本发明权利要求及其等同技术的范围之内,则本发明也意图包含这些改动和变型在内。
Claims (7)
1.一种煤岩水力压裂裂缝微震定位及扩展机理监测方法,其特征在于:具体包括以下步骤:
(a)根据煤岩体压裂的作业地点确定需要监测的重点区域;根据井下巷道的几何形状、施工现场环境和岩体的性质,布置X个微震传感器,记做Ci,i=1、2、3...X,所布置的微震传感器能够监测整个压裂区域;
(c)对微震传感器采集的原始波形信号进行滤波预处理,滤波的频率范围是由水力压裂诱发波形信号的频率来确定;所述原始波形信号滤波结束后,基于背景噪音振幅的平均值对微震传感器Ci采集的波形数据进行扫描,当任意微震传感器Cj采集波形信号的幅值超过时,将该微震传感器波形信号幅值超过的采样点记为起跳点Sj,起跳点所对应的时间点为Tj,然后从起跳点Sj前置一段时间Tb,后置一段时间Ta,将时间区间[Tb,Ta]定义为事件筛选窗口Wj;
(d)采用长短时窗法对事件筛选窗口Wj中所有微震传感器采集的波形信号进行扫描,j=1、2、3...,若微震传感器Ci中波形信号的特征函数值大于触发阈值λ时,判定微震传感器Ci内存在有效微震波形;当事件筛选窗口Wj内的有效微震波形数B≥m且m≥4时,确定该事件筛选窗口中存在有效微震事件,同时存储该有效微震事件的波形数据,否则,丢弃该事件筛选窗口内有效微震事件的波形数据;
(e)顺序扫描时刻Tj+Tb之后的微震波形信号,统计得到所有的有效微震事件数为E;
(f)利用波形敏感值对E个有效微震事件进行筛选;
上式中,表示第j个事件筛选窗口内第i个微震传感器采集的有效微震波形数,代表第j个事件筛选窗口内第i个微震传感器采集到的第k个有效微震波形的频率;当波形平均频率时,判定为优选信号,输出结果为1,反之输出结果为0;并对最终输出结果赋予权重w1;
上式中,代表第j个事件筛选窗口内第i个微震传感器采集到的第k个有效微震波形最大幅值和起跳点幅值的差值,表示第j个事件筛选窗口内第i个微震传感器采集到的第k个有效微震波形的上升时间;当平均幅值增长速率介于m3~m4,判断为优选信号,输出结果为1,反之输出结果为0,并对最终输出结果赋予权重w3;
上式中,表示第j个事件筛选窗口内第i个微震传感器采集的第k个有效微震波形的衰减时间,表示第j个事件筛选窗口内第i个微震传感器采集的第k个有效微震波形的持续时间;当平均衰减比介于m5~m6,判断为优选信号,输出结果为1,反之输出结果为0;并对最终输出结果赋予权重w4;
(5)计算波形敏感参数总权重,波形敏感参数总权重=(1|0)*w1+(1|0)*w2+(1|0)*w3+(1|0)*w4,当波形敏感参数总权重>w,判定这个事件为煤岩体水力压裂诱发产生的有效微震事件,按照时间顺序,依次对E个有效微震事件进行筛选,最终得到N个煤岩体水力压裂诱发产生的有效微震事件;
(g)利用局部赤迟信息准则法拾取煤岩体水力压裂诱发微震事件的高精度P波到时、P波初至振幅及初动信息;
(h)在已知高精度P波到时的基础上,构建目标函数;引入迭代法,在三维空间坐标系中构建一个四面体,并以四面体的四个顶点作为破裂源,计算目标函数的最大值和最小值,所述最大值和最小值计算完毕后进行比较,抛弃使得目标函数取得最大值的点,同时按照变换规则寻求新的顶点构建新的四面体;重复上述步骤直至求得目标函数的最小值满足迭代结束的条件,即为定位事件点;
(i)在微震事件高精度定位的基础上,利用采集的P波初至振幅信息结合定位事件点的坐标、微震传感器坐标反演出震源矩张量,并在主轴坐标系下将震源矩张量分解为同性部分ISO、纯剪切破坏部分DC、其他部分CLVD,同时根据P波初动信息绘出震源沙滩球,表示定位事件点的空间方位信息;
(j)将定位事件点的空间方位信息和坐标结合起来,准确绘出裂缝的空间形态,分析裂缝的延伸与贯通,判断起裂时的破坏类型,裂缝延伸时的破坏类型,延伸的裂缝遇到天然节面和层理时的破坏类型。
2.权利要求1所述的一种煤岩水力压裂裂缝微震定位及扩展机理监测方法,其特征在于:所述步骤(a)中微震传感器放置的数量至少为5个,对于室内水力压裂实验,微震传感器应至少布置在样品的两个表面上;对于现场煤岩体水力压裂,微震传感器应采用钻孔安装方式布置在稳定性较好的岩层中。
3.如权利要求1所述的一种煤岩水力压裂裂缝微震定位及扩展机理监测方法,其特征在于:所述步骤(c)中,对微震传感器采集的原始波形信号进行滤波预处理时,采用带通滤波器,滤波器的带宽为f1~f2Hz,仅允许高于较低的频率f1且低于较高的频率f2的特定频段的波形信号通过。
5.如权利要求1所述的一种煤岩水力压裂裂缝微震定位及扩展机理监测方法,其特征在于:所述步骤(f)中,波形敏感值包括波形平均频率、波形偏离值、平均幅值增长速率、平均衰减比。
7.如权利要求1所述的一种煤岩水力压裂裂缝微震定位及扩展机理监测方法,其特征在于:所述步骤(i)、(j)中,所述P波初至振幅信息反演水力压裂裂缝的破裂类型,需已知6个或6个以上微震传感器采集到的P波初至振幅,有效微震事件的破裂类型是根据纯剪切破坏部分MDC占震源矩张量比重PDC来确定,其中,PDC>60为剪切破裂,40<PDC<60为混合破裂,PDC>40为张拉破裂;所述P波初动信息是指P波从定位事件点传至地面时其初始震动的方向,即初动方向;将每个传感器采集的P波作适当的校正之后,把初动方向标到一个以震源为中心,且半径很小的球面上,即震源沙滩球上;当记录的数据足够多且传感器数量较多时,总可以找到两个互相垂直的节面,把震源沙滩球4等分;根据震源沙滩球上两个节面的方位就可以确定定位事件点的空间方位信息。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110307395.XA CN113050159B (zh) | 2021-03-23 | 2021-03-23 | 一种煤岩水力压裂裂缝微震定位及扩展机理监测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110307395.XA CN113050159B (zh) | 2021-03-23 | 2021-03-23 | 一种煤岩水力压裂裂缝微震定位及扩展机理监测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113050159A CN113050159A (zh) | 2021-06-29 |
CN113050159B true CN113050159B (zh) | 2021-11-16 |
Family
ID=76514368
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110307395.XA Active CN113050159B (zh) | 2021-03-23 | 2021-03-23 | 一种煤岩水力压裂裂缝微震定位及扩展机理监测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113050159B (zh) |
Families Citing this family (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113514878B (zh) * | 2021-07-09 | 2022-06-24 | 长沙矿山研究院有限责任公司 | 一种矿山微震波形类型自动识别方法 |
CN114137605A (zh) * | 2021-10-21 | 2022-03-04 | 陕西延长石油矿业有限责任公司 | 一种煤矿岩体中岩桥的智能定位及特征识别方法 |
CN114137600B (zh) * | 2021-11-26 | 2022-05-17 | 中国矿业大学 | 一种利用微震监测数据反演岩石破裂机理及失稳预测方法 |
CN114325823B (zh) * | 2021-12-30 | 2022-07-15 | 广西大学 | 一种基于nb-iot的岩体破裂失稳微震信号无线监测方法与装置 |
CN114417612B (zh) * | 2022-01-20 | 2024-03-29 | 北方工业大学 | 基于矩张量反演的采场微震震源机制求解方法 |
CN114412573B (zh) * | 2022-01-20 | 2023-08-08 | 中国矿业大学 | 一种矿井井壁破裂监测预警***及方法 |
CN114966849B (zh) * | 2022-04-27 | 2024-06-07 | 东北大学 | 基于微震或声发射及震源机制约束的岩体裂隙表征方法 |
CN115263427A (zh) * | 2022-05-25 | 2022-11-01 | 中煤科工开采研究院有限公司 | 覆岩破坏高度确定方法、装置、电子设备及存储介质 |
CN114966853B (zh) * | 2022-05-27 | 2023-03-31 | 中国矿业大学 | 基于微震监测信号确定冲击地点围岩运动参数极值的方法 |
CN116540299B (zh) * | 2023-07-05 | 2023-09-26 | 煤炭科学研究总院有限公司 | 一种用于煤矿场景的基于微震能量积聚倾向性的预警方法 |
CN117235428B (zh) * | 2023-11-16 | 2024-01-26 | 北京科技大学 | 煤岩损伤过程关键震源判识方法 |
CN117289344B (zh) * | 2023-11-24 | 2024-01-30 | 北京科技大学 | 基于震源空间分布的煤岩失稳破坏快速判识方法 |
CN118091754A (zh) * | 2023-11-30 | 2024-05-28 | 中国矿业大学 | 一种干热岩开采诱发地震机理分析方法和*** |
Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103336297A (zh) * | 2013-06-03 | 2013-10-02 | 北京京援伟达技术有限公司 | 微破裂向量扫描方法 |
CN103513272A (zh) * | 2012-06-19 | 2014-01-15 | 中国石油化工股份有限公司 | 一种微地震模拟监测方法 |
CN105528510A (zh) * | 2014-09-30 | 2016-04-27 | 中国石油化工股份有限公司 | 水力压裂微地震震源机制的反演方法及*** |
CN105629295A (zh) * | 2015-12-29 | 2016-06-01 | 四川圣诺油气工程技术服务有限公司 | 一种页岩气体积压裂微地震监测方法 |
WO2016134210A1 (en) * | 2015-02-20 | 2016-08-25 | Schlumberger Technology Corporation | Microseismic sensitivity analysis and scenario modelling |
CN106199703A (zh) * | 2016-08-26 | 2016-12-07 | 中国矿业大学 | 一种微震震源自动定位及可靠性综合评价方法 |
CN106353792A (zh) * | 2015-07-17 | 2017-01-25 | 中国石油化工股份有限公司 | 一种适用于水力压裂微震震源定位的方法 |
CN107728200A (zh) * | 2017-09-29 | 2018-02-23 | 中国石油化工股份有限公司 | 地面微地震压裂裂缝动态展布实时监测方法 |
CN111025392A (zh) * | 2019-12-27 | 2020-04-17 | 中国矿业大学 | 一种利用微震信号的煤岩体压裂裂缝实时快速监测评价方法 |
WO2020142074A1 (en) * | 2018-12-31 | 2020-07-09 | Halliburton Energy Services, Inc. | Distributed acoustic sensing: locating of microseismic events using travel time information with heterogeneous anisotropic velocity model |
-
2021
- 2021-03-23 CN CN202110307395.XA patent/CN113050159B/zh active Active
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103513272A (zh) * | 2012-06-19 | 2014-01-15 | 中国石油化工股份有限公司 | 一种微地震模拟监测方法 |
CN103336297A (zh) * | 2013-06-03 | 2013-10-02 | 北京京援伟达技术有限公司 | 微破裂向量扫描方法 |
CN105528510A (zh) * | 2014-09-30 | 2016-04-27 | 中国石油化工股份有限公司 | 水力压裂微地震震源机制的反演方法及*** |
WO2016134210A1 (en) * | 2015-02-20 | 2016-08-25 | Schlumberger Technology Corporation | Microseismic sensitivity analysis and scenario modelling |
CN106353792A (zh) * | 2015-07-17 | 2017-01-25 | 中国石油化工股份有限公司 | 一种适用于水力压裂微震震源定位的方法 |
CN105629295A (zh) * | 2015-12-29 | 2016-06-01 | 四川圣诺油气工程技术服务有限公司 | 一种页岩气体积压裂微地震监测方法 |
CN106199703A (zh) * | 2016-08-26 | 2016-12-07 | 中国矿业大学 | 一种微震震源自动定位及可靠性综合评价方法 |
CN107728200A (zh) * | 2017-09-29 | 2018-02-23 | 中国石油化工股份有限公司 | 地面微地震压裂裂缝动态展布实时监测方法 |
WO2020142074A1 (en) * | 2018-12-31 | 2020-07-09 | Halliburton Energy Services, Inc. | Distributed acoustic sensing: locating of microseismic events using travel time information with heterogeneous anisotropic velocity model |
CN111025392A (zh) * | 2019-12-27 | 2020-04-17 | 中国矿业大学 | 一种利用微震信号的煤岩体压裂裂缝实时快速监测评价方法 |
Non-Patent Citations (2)
Title |
---|
千秋煤矿微震震源参数特征以及震源机制分析;陈栋 等;《煤炭学报》;20190731;第44卷(第7期);第2011-2019页 * |
微震监测技术及其在煤矿的应用现状与展望;李楠 等;《煤炭学报》;20170731;第42卷;第83-96页 * |
Also Published As
Publication number | Publication date |
---|---|
CN113050159A (zh) | 2021-06-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN113050159B (zh) | 一种煤岩水力压裂裂缝微震定位及扩展机理监测方法 | |
CN106353792A (zh) | 一种适用于水力压裂微震震源定位的方法 | |
CN107479098B (zh) | 一种水力压裂过程中同井微地震监测工艺 | |
CN109736796B (zh) | 一种加深炮孔超前探测机构及其预报方法 | |
CN106248672B (zh) | 一种基于dic技术的现场孔内岩体裂纹扩展模式识别方法及*** | |
CN101968550A (zh) | 基于阵列光纤传感器的岩层识别装置及方法 | |
CN112731525B (zh) | 微震与电磁辐射同步监测的巷道围岩稳定性智能预测方法 | |
CN111025392B (zh) | 一种利用微震信号的煤岩体压裂裂缝实时快速监测评价方法 | |
CN111077576A (zh) | 地面压裂的监测方法 | |
RU2467171C1 (ru) | Способ диагностики опасных ситуаций при подземной добыче каменного угля и методика прогноза параметров зон трещиноватости, образованной гидроразрывом пласта | |
CN113153430B (zh) | 巷道围岩破坏声发射定位与波速成像监测及灾变预警方法 | |
CN112965136A (zh) | 一种富水岩溶隧道的多手段超前探测方法 | |
CN113360592A (zh) | 基于微震监测的n00工法煤层顶板突水危险性预警方法 | |
CN205743991U (zh) | 一种煤矿冲击地压联合监测预警*** | |
CN111708079A (zh) | 基于tsp的隧道有害气体综合超前预报方法 | |
CN111852566A (zh) | 煤矿井下水力压裂裂缝扩展规律的微震监测方法 | |
Butcher et al. | Evaluating rock mass disturbance within open-pit excavations using seismic methods: A case study from the Hinkley Point C nuclear power station | |
CN107942372A (zh) | 井间及井地联合地震波ct成像方法及装置 | |
CN110259442B (zh) | 一种煤系地层水力压裂破裂层位识别方法 | |
Igonin et al. | Seismic anisotropy reveals stress changes around a fault as it is activated by hydraulic fracturing | |
CN201837728U (zh) | 基于阵列光纤传感器的岩层识别装置 | |
CN106032750B (zh) | 一种基于钻探能谱的地质编录仪 | |
CN104932006A (zh) | 隧道地质超前预报中人工地震波的机械能自动化无损激发装置 | |
CN114658475A (zh) | 一种基于采动裂隙场量化优选的瓦斯精准抽采方法 | |
CN112748476B (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 |