CN113107599A - 采煤工作面内隐伏构造透射槽波相邻道振幅比成像方法 - Google Patents
采煤工作面内隐伏构造透射槽波相邻道振幅比成像方法 Download PDFInfo
- Publication number
- CN113107599A CN113107599A CN202110400859.1A CN202110400859A CN113107599A CN 113107599 A CN113107599 A CN 113107599A CN 202110400859 A CN202110400859 A CN 202110400859A CN 113107599 A CN113107599 A CN 113107599A
- Authority
- CN
- China
- Prior art keywords
- wave
- equal
- attenuation
- amplitude
- imaging
- 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
- 239000003245 coal Substances 0.000 title claims abstract description 109
- 238000003384 imaging method Methods 0.000 title claims abstract description 91
- 230000005540 biological transmission Effects 0.000 title claims abstract description 56
- 238000000034 method Methods 0.000 claims abstract description 47
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 4
- 230000005284 excitation Effects 0.000 claims description 16
- 238000005070 sampling Methods 0.000 claims description 6
- 238000005065 mining Methods 0.000 claims description 5
- 239000008186 active pharmaceutical agent Substances 0.000 claims description 3
- 238000001914 filtration Methods 0.000 claims description 3
- 238000006467 substitution reaction Methods 0.000 claims description 3
- 230000009466 transformation Effects 0.000 claims description 3
- 238000001514 detection method Methods 0.000 description 16
- 238000009792 diffusion process Methods 0.000 description 15
- 239000006185 dispersion Substances 0.000 description 15
- 238000010521 absorption reaction Methods 0.000 description 14
- 230000000694 effects Effects 0.000 description 8
- 230000001902 propagating effect Effects 0.000 description 7
- 238000010276 construction Methods 0.000 description 4
- 230000007547 defect Effects 0.000 description 4
- 235000019687 Lamb Nutrition 0.000 description 2
- 230000005856 abnormality Effects 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000000644 propagated effect Effects 0.000 description 2
- 230000005855 radiation Effects 0.000 description 2
- 239000000523 sample Substances 0.000 description 2
- 239000004576 sand Substances 0.000 description 2
- 238000007789 sealing Methods 0.000 description 2
- 238000000926 separation method Methods 0.000 description 2
- 230000002411 adverse Effects 0.000 description 1
- 230000002238 attenuated effect Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 238000003786 synthesis reaction Methods 0.000 description 1
- 238000003325 tomography Methods 0.000 description 1
Images
Classifications
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21F—SAFETY DEVICES, TRANSPORT, FILLING-UP, RESCUE, VENTILATION, OR DRAINING IN OR OF MINES OR TUNNELS
- E21F17/00—Methods or devices for use in mines or tunnels, not covered elsewhere
- E21F17/18—Special adaptations of signalling or alarm devices
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21C—MINING OR QUARRYING
- E21C41/00—Methods of underground or surface mining; Layouts therefor
- E21C41/16—Methods of underground mining; Layouts therefor
- E21C41/18—Methods of underground mining; Layouts therefor for brown or hard coal
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21F—SAFETY DEVICES, TRANSPORT, FILLING-UP, RESCUE, VENTILATION, OR DRAINING IN OR OF MINES OR TUNNELS
- E21F17/00—Methods or devices for use in mines or tunnels, not covered elsewhere
Landscapes
- Engineering & Computer Science (AREA)
- Mining & Mineral Resources (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geochemistry & Mineralogy (AREA)
- Geology (AREA)
- Remote Sensing (AREA)
- Geophysics And Detection Of Objects (AREA)
- Analysing Materials By The Use Of Radiation (AREA)
Abstract
Description
技术领域
本发明涉及煤矿开采技术领域,具体涉及一种采煤工作面内隐伏构造透射槽波衰减成像方法。
背景技术
在煤矿开采过程中,由于煤层地质条件复杂,煤层内部发育大量尺度近似煤层厚度的小规模地质构造如断层、陷落柱等,传统的地面探测方法分辨率不足,不能精确识别煤层内的小构造异常。采煤工作面内隐伏的未知小规模构造异常成为致灾隐患,严重影响采煤作业。透射槽波衰减成像是探测采煤工作面内隐伏构造的有效方法,具有探测距离大、构造成像精度高等优势。
当前应用透射槽波资料进行采煤工作面内隐伏构造衰减成像均采用常规振幅衰减成像方法。
设煤层的总衰减系数为α,槽波在煤层中传播,接收点接收到的槽波振幅A与槽波激发的初始振幅A0和传播距离L之间的关系如下:
A=A0exp(-αL)
则
对上式两边取自然对数得,
基于以上理论,目前透射槽波采煤工作面内隐伏构造常规振幅衰减成像方法的技术方案如下:
1).采用如图1所示的观测***,在采煤工作面一侧巷道煤层内布置M个炮点激发槽波,在对面巷道煤层内布置P个接收点接收每一个炮点激发的共M个共炮点槽波数据。
2).将需要构造成像的采煤工作面范围进行网格划分,划分成N个小网格,如图1所示。令每个小网格的衰减系数为αn,1≤n≤N且n为整数。图1所示的观测***共形成M×P条射线路径。令每一条射线路径Li,j在成像范围内每一个小网格上传播的长度为1≤i≤M且i为整数,1≤j≤P且j为整数。
3).在槽波数据中提取沿每条射线路径Li,j传播被检波器接收到的透射槽波振幅值。设第i个炮点激发的槽波初值振幅为Ai。第j个检波器接收到的第i个炮点激发沿传播路径Li,j传播的透射槽波振幅为Ai,j(i=1,2,…,M;j=1,2,…,P)。
令
列出线性方程:
5).应用代数重建方法、联合代数重建法或联合迭代重建法对上面方程组进行求解,得到采煤工作面成像范围内N个网格衰减系数αn(n=1,2,…,N)。
6).用αn(n=1,2,…,N)对采煤工作面成像范围进行衰减系数成像,根据αn异常确定采煤工作面内构造分布。
常规振幅衰减成像方法必须已知初始振幅A0,并且要求所有震源激发的透射槽波初始振幅必须保持一致。但是,由于槽波是典型的频散波,其振幅沿射线路径并不是纯衰减过程,在距离震源的一定距离内存在槽波的合成过程,在此过程中随着传播距离增大槽波振幅逐渐增大,因此透射槽波衰减规律不同于非频散或弱频散的体波衰减规律。理论上可以在炮点附近布置检波器来记录震源的初始振幅,但是在现场实际槽波探测中,增加炮点初始振幅检波器会增加施工工作量和风险,由于受煤层厚度和震源主频的影响槽波最大振幅位置与震源之间的距离不固定,难以确定记录最大初始槽波振幅的记录位置。并且由于距离炮点过近初始振幅过大往往超出检波器最大振幅量程,从而难以精确记录每个震源的初始振幅。同时由于各炮点激发药量、炮孔深度和炮泥封孔效果难以保持一致使得各震源激发的槽波初始振幅难以保持一致。并且常规振幅衰减成像方法求得衰减系数中难以消除煤层吸收衰减、波前扩散衰减、频散衰减和槽波向顶底板漏能衰减等非构造因素造成的振幅衰减对构造成像的影响。以上不利条件严重影响了常规振幅衰减成像方法在透射槽波衰减采煤工作面内隐伏构造成像中的应用效果。
上述透射槽波常规振幅衰减采煤工作面内隐伏构造成像方法存在以下不足:
1、由于强烈的频散特性,槽波振幅沿射线路径并不是纯衰减过程,在距离震源的一定距离内槽波振幅随传播距离不但不衰减反而增大,透射槽波衰减规律不同于非频散或弱频散的体波衰减规律,常规振幅衰减成像方法不适用于频散波振幅衰减成像;
2、常规振幅衰减成像方法需要精确求取震源激发的透射槽波初始振幅,而现场槽波探测中槽波初始振幅精确记录不可实现;
3、常规振幅衰减成像要求所有炮点激发的槽波初始振幅保持一致,但是现场探测时由于各炮点激发药量、炮孔深度和炮泥封孔效果不能保持一致使得各震源激发的槽波初始振幅难以保持一致;
4、常规振幅衰减成像方法求不能消除煤层吸收衰减、波前扩散衰减、频散衰减、槽波向顶底板漏能衰减等与构造无关的因素造成的振幅衰减对构造成像的影响。
同时,也有相关文献被公开,如蔺国华,透射槽波方法在中厚煤层隐伏小断层探测中的应用[J],中国煤炭地质,2018,30(1):113-117;王季,李建政,吴海,等,透射槽波能量衰减系数成像与陷落柱探测[J].煤炭科学技术,2015,43(1):108-111;段天柱,透射槽波层析成像在煤矿工作面隐伏构造探测中的应用[J].矿业安全与环保,2018,45(2):68-71等,以上文献均采用常规振幅衰减成像方法进行采煤工作面内隐伏构造成像。其成像过程中,初始振幅采用没有误差约束的估计方法得到;没有对不同炮点激发的槽波能量不同进行精确校正;没有采用相应的措施消除煤层吸收衰减、波前扩散衰减、频散衰减、槽波向顶底板漏能衰减等与构造无关的因素造成的振幅衰减的影响,因此难以保证成像精度。
发明内容
本发明的目的是为克服上述现有技术的不足,提供一种相邻地震道振幅比值构造衰减采煤工作面内隐伏构造透射槽波成像方法。
透射槽波在煤层中传播,造成其振幅衰减的主要因素有:煤层吸收衰减、波前扩散衰减、频散衰减、构造衰减和漏能衰减;因此,透射槽波振幅衰减是以上所有因素共同作用的总效应。设总衰减系数为α,则α可以等效为煤层吸收衰减系数αα、波前扩散衰减系数αg、频散衰减系数αd、构造衰减系数αs和漏能衰减系数αl之和,即:
α=αα+αg+αd+αs+αl
应用煤层对槽波的总衰减系数进行构造成像受到非构造衰减系数的影响较大。在非构造衰减系数相对于构造衰减系数较大的情况下,常规振幅衰减成像方法的技术方案难以对采煤工作面内隐伏构造进行精确成像。
透射槽波振幅A、槽波初始振幅A0和传播距离L与各吸收系数的关系表达式为:
A=A0exp(-(αα+αg+αd+αl+αs)L),
透射槽波探测在煤层内激发槽波,并在同一煤层内接收透射槽波,透射槽波仅在煤层内顺层传播。在煤层较均匀或非均质变化较弱的条件下,煤层吸收衰减系数αα、波前扩散衰减系数αg、频散衰减系数αd和漏能衰减系数αl在局部范围内变化不大各自近似相等,它们造成的透射槽波振幅衰减与射线传播距离L长度有关,随着射线传播距离L长度增大振幅衰减的越多。而构造衰减系数αs决定于射线是否穿过构造。
沿采煤工作面一侧煤层均匀密集布置炮点,在对侧均匀密集布置检波器如图1所示。设第i炮激发槽波的初始振幅为Ai,透射槽波沿煤层传播被相邻的第j-1,j两个接收点接收,两条相邻射线传播路径长度分别为Li,j-1,Li,j。透射槽波沿两条射线路径传播的总衰减系数分别为αi,j-1和αi,j。j-1和j接收点接收到的透射槽波振幅分别为Ai,j-1,Ai,j。
则
αi,j-1=ααi,j-1+αgi,j-1+αdi,j-1+αli,j-1+αsi,j-1
其中ααi,j-1、αgi,j-1、αdi,j-1、αli,j-1、αsi,j-1分比为透射槽波沿Li,j-1射线路径传播的煤层吸收衰减系数、波前扩散衰减系数、频散衰减系数、漏能衰减系数和构造衰减系数。
αi,j=ααi,j+αgi,j+αdi,j+αli,j+αsi,j
其中ααi,j、αgi,j、αdi,j、αli,j、αsi,j分别为透射槽波沿Li,j射线路径传播的煤层吸收衰减系数、波前扩散衰减系数、频散衰减系数、漏能衰减系数和构造衰减系数。
则可以得到透射槽波沿两条射线路径Li,j和Li,j-1传播的振幅Ai,j-1和Ai,j与初始振幅Ai和各自路径衰减系数的关系如下。
Ai,j-1=Aiexp(-(ααi,j-1+αgi,j-1+αdi,j-1+αli,j-1)Li,j-1-αsi,j-1Li,j-1)
Ai,j=Aiexp(-(ααi,j+αgi,j+αdi,j+αli,j)Li,j-αsi,jLi,j)
为了保证对采煤工作面内的构造进行精确探测,要求炮点和检波点间距小于要求探测构造最小尺寸的一半。为了保证相邻检波点距离相对于其射线路径长度足够小要求检波点间距小于最小炮检距的十分之一。在以上条件下两条相邻的共炮点射线路径近似重合长度近似相等,在微小的局部范围内除构造衰减系数外,其它衰减系数近似相等。而每条射线路径受是否穿过构造位置影响,其构造衰减系数差异较大。
令上面两式相比,
由以上分析可知,在煤层较均匀或非均质变化较弱的条件下,并且检波器均匀密集布置,相邻检波点距离相对于其射线路径长度足够小。透射槽波相邻地震道振幅比值主要与构造衰减系数有关,受非构造衰减因素影响较小。由此提出一种不需要求取槽波初始振幅、不需要各震源激发的槽波初始振幅保持一致,可以有效克服非构造衰减因影响的采煤工作面内隐伏构造成像方法。
采煤工作面内隐伏构造透射槽波相邻道振幅比成像方法,步骤如下:
S1.采用如图1所示的观测***,在采煤工作面一侧巷道煤层内均匀布置M个炮点激发槽波,在对面巷道煤层内均匀布置P个接收器,每个接收器均接收每一个炮点激发的共M个共炮点透射槽波数据,共形成M×P条射线路径,M和P均为正整数,炮点和检波器之间最小距离为L,需要探测的隐伏构造在平行工作面的平面内最小尺寸要求为D,炮点间距离DS,要求接收器间距为DR,要求同时要求
S2.将需要构造成像的采煤工作面范围进行网格划分,划分成N个边长相等的正方形小网格,N为正整数,正方形网格边长为Lg,为了保证构造成像精度要求令每个小网格的构造衰减系数为1≤n≤N且n为整数;令每一条射线路径Li,j在成像范围内每一个小网格上传播的长度为其中1≤i≤M且i为整数,1≤j≤P且j为整数;
S3.采用改进的S变换离散化算法提取槽波数据中沿每一条长度为Li,j的射线路径传播被检波器接收到的透射槽波振幅值Ai,j,具体步骤如下:第i炮激发沿长度为Li,j的射线路径传播被第j个检波器接收的地震记录信号为xi,j(t),其中t为时间采样点,0≤t≤T-1且t为整数,T为2的非负整数次幂;采用下面公式进行改进的S变换,变换为时间-频率域信号Si,j[t,f],其中f为频率采样点,0≤f≤T-1且f为整数,T为2的非负整数次幂;
在Si,j[t,f]中识别出能量最强的透射槽波埃里相(Airy相)的时间域范围[t1,t2]和频率域范围[f1,f2],其中t2>t1,1≤t1≤T-1,1≤t2≤T-1,f2>f1,1≤f1≤T-1,1≤f2≤T-1);对Si,j[t,f]在时间域和频率域进行[t1,t2]和[f1,f2]带通滤波滤除除透射槽波信号外的其它地震信号,然后进行反S变换得到第i炮激发沿长度为Li,j的射线路径传播被第j个检波器接收的透射槽波时域信号的的xci,j(t);在xci,j(t)中取得最大振幅作为透射槽波振幅值Ai,j;
S4.第i个炮点激发的槽波分别沿传播路径Li,j-1、Li,j传播的透射槽波振幅分别为Ai,j-1、Ai,j,其中1≤i≤M且i为整数,2≤j≤P且j为整数,传播路径Li,j-1、Li,j在成像范围内每一个小网格上的传播长度为
本发明的有益效果是:
1、本方法采用相邻道透射槽波振幅比值构建线性方程组求解构造衰减系数,不需要透射槽波的初始振幅。改进了常规振幅衰减成像方法求解衰减系数时必须需要透射槽波的初始振幅,而透射槽波初始振幅几乎不可能精确测得的缺点。本方法有效避免了常规振幅衰减成像方法由于初始振幅选择不当对工作面内隐伏构造成像的误差。
2、本方法在求解透射槽波构造衰减系数时,不要求所有炮点激发的槽波初始振幅必须保持一致的条件。本方法有效克服了常规振幅衰减成像方法由于初始振幅不一致对工作面内隐伏构造成像的误差。
3、常规振幅衰减成像方法求得的透射槽波衰减系数是非构造因素造成的透射槽波振幅衰减(煤层吸收衰减、波前扩散衰减、频散衰减、槽波向顶底板漏能衰减等)和构造引起的透射槽波振幅衰减的总效应。因此用常规振幅衰减成像方法求的衰减系数进行构造成像,不能消除非构造衰减因素的影响存在较大误差。本方法求得的构造衰减系数主要反演采煤工作面内隐伏的构造特征,可以有效减弱非构造衰减因素的影响。因此,用本方法求得的构造衰减系数进行构造成像精度更高。
4、常规振幅衰减成像方法进行采煤工作面内隐伏构造成像时,需要进行波前扩散校正。本方法进行采煤工作面内隐伏构造成像不需要波前扩散校正,简化了构造成像过程。
5、发明的本方法同样适用于板状构件缺陷兰姆波(Lamb Wave)衰减检测、板状建筑构件的破裂、损伤衰减检测和地表浅层面波衰减构造探测。
附图说明
图1是透射槽波探测示意图;
图2a是采煤工作面隐伏构造模型图;
图2b是常规振幅衰减成像方法成像(初始振幅0.0001);
图2c是常规振幅衰减成像方法成像(初始振幅0.0002);
图2d是常规振幅衰减成像方法成像(初始振幅0.0005);
图2e是常规振幅衰减成像方法成像(初始振幅0.001);
图2f是常规振幅衰减成像方法成像(初始振幅0.002);
图2g是常规振幅衰减成像方法成像(初始振幅0.005);
图2h是常规振幅衰减成像方法成像(初始振幅0.01);
图2i是常规振幅衰减成像方法成像(初始振幅0.1);
图2j是本方法衰减成像(不需要初始振幅)。
具体实施方式
下面结合附图和实施例对本发明进一步说明。
本说明书所附图式所绘示的结构、比例、大小等,均仅用以配合说明书所揭示的内容,以供熟悉此技术的人士了解与阅读,并非用以限定本发明可实施的限定条件,故不具技术上的实质意义,任何结构的修饰、比例关系的改变或大小的调整,在不影响本发明所能产生的功效及所能达成的目的下,均应仍落在本发明所揭示的技术内容涵盖的范围内。同时,本说明书中所引用的如“上”、“下”、“左”、“右”、“中间”及“一”等的用语,亦仅为便于叙述的明了,而非用以限定本发明可实施的范围,其相对关系的改变或调整,在无实质变更技术内容下,当亦视为本发明可实施的范畴。
透射槽波在煤层中传播,造成其振幅衰减的主要因素有:煤层吸收衰减、波前扩散衰减、频散衰减、构造衰减和漏能衰减;因此,透射槽波振幅衰减是以上所有因素共同作用的总效应。设总衰减系数为α,则α可以等效为煤层吸收衰减系数αα、波前扩散衰减系数αg、频散衰减系数αd、构造衰减系数αs和漏能衰减系数αl之和,即:
α=αα+αg+αd+αs+αl
应用煤层对槽波的总衰减系数进行构造成像受到非构造衰减系数的影响较大。在非构造衰减系数相对于构造衰减系数较大的情况下,常规振幅衰减成像方法的技术方案难以对采煤工作面内隐伏构造进行精确成像。
透射槽波振幅A、槽波初始振幅A0和传播距离L与各吸收系数的关系表达式为:
A=A0exp(-(αα+αg+αd+αl+αs)L),
透射槽波探测在煤层内激发槽波,并在同一煤层内接收透射槽波,透射槽波仅在煤层内顺层传播。在煤层较均匀或非均质变化较弱的条件下,煤层吸收衰减系数αα、波前扩散衰减系数αg、频散衰减系数αd和漏能衰减系数αl在局部范围内变化不大各自近似相等,它们造成的透射槽波振幅衰减与射线传播距离L长度有关,随着射线传播距离L长度增大振幅衰减的越多。而构造衰减系数αs决定于射线是否穿过构造,即受射线在介质中相对于构造的位置影响。
沿采煤工作面一侧煤层均匀密集布置炮点,在对侧均匀密集布置检波器如图1所示。设第i炮激发槽波的初始振幅为Ai,透射槽波沿煤层传播被相邻的第j-1,j两个接收点接收,两条相邻射线传播路径长度分别为Li,j-1,Li,j。透射槽波沿两条射线路径传播的总衰减系数分别为αi,j-1和αi,j。j-1和j接收点接收到的透射槽波振幅分别为Ai,j-1,Ai,j。
则
αi,j-1=ααi,j-1+αgi,j-1+αdi,j-1+αli,j-1+αsi,j-1
其中ααi,j-1、αgi,j-1、αdi,j-1、αli,j-1、αsi,j-1分别为透射槽波沿Li,j-1射线路径传播的煤层吸收衰减系数、波前扩散衰减系数、频散衰减系数、漏能衰减系数和构造衰减系数。
αi,j=ααi,j+αgi,j+αdi,j+αli,j+αsi,j
其中ααi,j、αgi,j、αdi,j、αli,j、αsi,j分比为透射槽波沿Li,j射线路径传播的煤层吸收衰减系数、波前扩散衰减系数、频散衰减系数、漏能衰减系数和构造衰减系数。
则可以得到透射槽波沿两条射线路径Li,j和Li,j-1传播的振幅Ai,j-1和Ai,j与初始振幅Ai和各自路径衰减系数的关系如下。
Ai,j-1=Aiexp(-(ααi,j-1+αgi,j-1+αdi,j-1+αli,j-1)Li,j-1-αsi,j-1Li,j-1)
Ai,j=Aiexp(-(ααi,j+αgi,j+αdi,j+αli,j)Li,j-αsi,jLi,j)
为了保证对采煤工作面内的构造进行精确探测,要求炮点和检波点间距小于要求探测构造最小尺寸的一半。为了保证相邻检波点距离相对于其射线路径长度足够小要求检波点间距小于最小炮检距的十分之一。在以上条件下两条相邻的共炮点射线路径近似重合长度近似相等,在微小的局部范围内除构造衰减系数外,其它衰减系数近似相等。而每条射线路径受是否穿过构造位置影响,其构造衰减系数差异较大。
令上面两式相比。
由以上分析可知,在煤层较均匀或非均质变化较弱的条件下,并且检波器均匀密集布置,相邻检波点距离相对于其射线路径长度足够小。透射槽波相邻地震道振幅比值主要与构造衰减系数有关,受非构造衰减因素影响较小。由此提出一种不需要求取槽波初始振幅、不需要各震源激发的槽波初始振幅保持一致,可以有效克服非构造衰减因影响的采煤工作面内隐伏构造成像方法。
采煤工作面内隐伏构造透射槽波相邻道振幅比成像方法,步骤如下:
S1.采用如图1所示的观测***,在采煤工作面一侧巷道煤层内均匀布置M个炮点激发槽波,在对面巷道煤层内均匀布置P个接收器,每个接收器均接收每一个炮点激发的共M个共炮点透射槽波数据,共形成M×P条射线路径,M和P均为正整数,炮点和检波器之间最小距离为L,需要探测的隐伏构造在平行工作面的平面内最小尺寸为D,炮点间距离DS,要求接收器间距为DR,要求同时要求
S2.将需要构造成像的采煤工作面范围进行网格划分,划分成N个边长相等的正方形小网格,N为正整数,正方形网格边长为Lg,为了保证构造成像精度要求令每个小网格的构造衰减系数为1≤n≤N且n为整数;令每一条射线路径Li,j在成像范围内每一个小网格上传播的长度为其中1≤i≤M且i为整数,1≤j≤P且j为整数;
S3.采用改进的S变换离散化算法提取槽波数据中沿每一条长度为Li,j的射线路径传播被检波器接收到的透射槽波振幅值Ai,j,具体步骤如下:第i炮激发沿长度为Li,j的射线路径传播被第j个检波器接收的地震记录信号为xi,j(t),其中t为时间采样点,0≤t≤T-1且t为整数,T为2的非负整数次幂;采用下面公式进行改进的S变换,变换为时间-频率域信号Si,j[t,f],其中f为频率采样点,0≤f≤T-1且f为整数,T为2的非负整数次幂;
在Si,j[t,f]中识别出能量最强的透射槽波埃里相(Airy相)的时间域范围[t1,t2]和频率域范围[f1,f2],其中t2>t1,1≤t1≤T-1,1≤t2≤T-1,f2>f1,1≤f1≤T-1,1≤f2≤T-1);对Si,j[t,f]在时间域和频率域进行[t1,t2]和[f1,f2]带通滤波滤除除透射槽波信号外的其它地震信号,然后进行反S变换得到第i炮激发沿长度为Li,j的射线路径传播被第j个检波器接收的透射槽波时域信号的的xci,j(t);在xci,j(t)中取得最大振幅作为透射槽波振幅值Ai,j;
S4.第i个炮点激发的槽波分别沿传播路径Li,j-1、Li,j传播的透射槽波振幅分别为Ai,j-1、Ai,j,其中1≤i≤M且i为整数,2≤j≤P且j为整数,传播路径Li,j-1、Li,j在成像范围内每一个小网格上的传播长度为
设计如图2a所示的采煤工作面数值模型,通过人工数值模拟透射槽波数据分别用常规振幅衰减成像方法和本方法进行构造成像。如图2b-图2j所示,透射槽波常规振幅衰减成像方法采煤工作面内构造成像效果依赖于初始振幅选择。而实际探测中,透射槽波的初始振幅未知,盲目选择不当的初始振幅严重影响构造成像效果。而本方法不需要透射槽波的初始振幅,在相同的槽波资料的条件下一次计算能得到高精度的采煤工作面内隐伏构造的精确成像。
上述虽然结合附图对本发明的具体实施方式进行了描述,但并非对本发明保护范围的限制,所属领域技术人员应该明白,在本发明的技术方案的基础上,本领域技术人员不需要付出创造性劳动即可做出的各种修改或变形仍在本发明的保护范围以内。
Claims (3)
1.一种采煤工作面内隐伏构造透射槽波相邻道振幅比成像方法,其特征是,步骤如下:
S1.在采煤工作面一侧巷道煤层内均匀布置M个炮点激发槽波,在对面巷道煤层内均匀布置P个接收器,每个接收器均接收每一个炮点激发的共M个共炮点透射槽波数据,共形成M×P条射线路径,M和P均为正整数,炮点和检波器之间最小距离为L,需要探测的隐伏构造在平行工作面的平面内最小尺寸为D,炮点间距离DS,接收器间距为DR;
S2.将需要构造成像的采煤工作面范围进行网格划分,划分成N个边长相等的正方形小网格,N为正整数,正方形网格边长为Lg,要求令每个小网格的构造衰减系数为1≤n≤N且n为整数,令每一条射线路径Li,j在成像范围内每一个小网格上传播的长度为其中1≤i≤M且i为整数,1≤j≤P且j为整数;
S3.采用改进的S变换离散化算法提取槽波数据中沿每一条长度为Li,j的射线路径传播被检波器接收到的透射槽波振幅值Ai,j;
S4.第i个炮点激发的槽波分别沿传播路径Li,j-1、Li,j传播的透射槽波振幅分别为Ai,j-1、Ai,j,其中1≤i≤M且i为整数,2≤j≤P且j为整数,传播路径Li,j-1、Li,j在成像范围内每一个小网格上的传播长度为
3.如权利要求1所述的采煤工作面内隐伏构造透射槽波相邻道振幅比成像方法,其特征是,所述步骤S3中,计算透射槽波振幅值Ai,j的具体步骤如下:
第i炮激发沿长度为Li,j的射线路径传播被第j个检波器接收的地震记录信号为xi,j(t),其中t为时间采样点,0≤t≤T-1且t为整数,T为2的非负整数次幂;采用下面公式进行改进的S变换,变换为时间-频率域信号Si,j[t,f],其中f为频率采样点,0≤f≤T-1且f为整数,T为2的非负整数次幂;
在Si,j[t,f]中识别出能量最强的透射槽波埃里相的时间域范围[t1,t2]和频率域范围[f1,f2],其中t2>t1,1≤t1≤T-1,1≤t2≤T-1,f2>f1,1≤f1≤T-1,1≤f2≤T-1);对Si,j[t,f]在时间域和频率域进行[t1,t2]和[f1,f2]带通滤波滤除除透射槽波信号外的其它地震信号,然后进行反S变换得到第i炮激发沿长度为Li,j的射线路径传播被第j个检波器接收的透射槽波时域信号的的xci,j(t);在xci,j(t)中取得最大振幅作为透射槽波振幅值Ai,j。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110400859.1A CN113107599B (zh) | 2021-04-14 | 2021-04-14 | 采煤工作面内隐伏构造透射槽波相邻道振幅比成像方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110400859.1A CN113107599B (zh) | 2021-04-14 | 2021-04-14 | 采煤工作面内隐伏构造透射槽波相邻道振幅比成像方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113107599A true CN113107599A (zh) | 2021-07-13 |
CN113107599B CN113107599B (zh) | 2022-07-15 |
Family
ID=76717000
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110400859.1A Active CN113107599B (zh) | 2021-04-14 | 2021-04-14 | 采煤工作面内隐伏构造透射槽波相邻道振幅比成像方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113107599B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114814946A (zh) * | 2022-04-24 | 2022-07-29 | 安徽理工大学 | 一种基于透射槽波相邻道质心频率的层析成像方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106610503A (zh) * | 2016-11-29 | 2017-05-03 | 北京中矿大地地球探测工程技术有限公司 | 一种煤矿采掘过程中的全方位槽波地震探测方法 |
CN107015273A (zh) * | 2017-04-10 | 2017-08-04 | 中煤科工集团西安研究院有限公司 | 一种基于频率域的煤层Love型槽波压缩方法及装置 |
CN107831530A (zh) * | 2017-10-25 | 2018-03-23 | 山东科技大学 | 厚煤层沿底板或顶板掘进煤巷反射槽波超前探测方法 |
CN111650636A (zh) * | 2020-04-28 | 2020-09-11 | 徐州天浩勘探有限公司 | 一种煤层底板隐伏构造探测方法 |
-
2021
- 2021-04-14 CN CN202110400859.1A patent/CN113107599B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106610503A (zh) * | 2016-11-29 | 2017-05-03 | 北京中矿大地地球探测工程技术有限公司 | 一种煤矿采掘过程中的全方位槽波地震探测方法 |
CN107015273A (zh) * | 2017-04-10 | 2017-08-04 | 中煤科工集团西安研究院有限公司 | 一种基于频率域的煤层Love型槽波压缩方法及装置 |
CN107831530A (zh) * | 2017-10-25 | 2018-03-23 | 山东科技大学 | 厚煤层沿底板或顶板掘进煤巷反射槽波超前探测方法 |
CN111650636A (zh) * | 2020-04-28 | 2020-09-11 | 徐州天浩勘探有限公司 | 一种煤层底板隐伏构造探测方法 |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114814946A (zh) * | 2022-04-24 | 2022-07-29 | 安徽理工大学 | 一种基于透射槽波相邻道质心频率的层析成像方法 |
CN114814946B (zh) * | 2022-04-24 | 2024-04-05 | 安徽理工大学 | 一种基于透射槽波相邻道质心频率的层析成像方法 |
Also Published As
Publication number | Publication date |
---|---|
CN113107599B (zh) | 2022-07-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102937721B (zh) | 利用初至波走时的有限频层析成像方法 | |
Lebedev et al. | Mapping the Moho with seismic surface waves: A review, resolution analysis, and recommended inversion strategies | |
Huthwaite | Evaluation of inversion approaches for guided wave thickness mapping | |
CN103630933B (zh) | 基于非线性优化的时空域交错网格有限差分方法和装置 | |
WO2017167191A1 (zh) | 地震数据处理方法和装置 | |
CN102636568B (zh) | 一种检测混凝土内部缺陷的有限元超声成像方法 | |
EA019484B1 (ru) | Способ получения сейсмических данных | |
CN105158796B (zh) | 确定toc含量的方法和装置 | |
Dalton et al. | The thermal structure of cratonic lithosphere from global Rayleigh wave attenuation | |
US11906472B2 (en) | Non-destructive concrete stress evaluation | |
CN104775811A (zh) | 一种地层各向异性信息的提取及校正方法和*** | |
CN113107599B (zh) | 采煤工作面内隐伏构造透射槽波相邻道振幅比成像方法 | |
Duret et al. | Near-surface velocity modeling using a combined inversion of surface and refracted P-waves | |
CN104990994A (zh) | 基于脉冲压缩技术的沉积物中声速和衰减的宽带测量方法 | |
Niu et al. | Integrated assessment of concrete structure using Bayesian theory and ultrasound tomography | |
Vilhelm et al. | Velocity dispersion in fractured rocks in a wide frequency range | |
CN114994752A (zh) | 一种弹性波地震频率谐振勘探方法 | |
CN105571986B (zh) | 一种基于散射能谱双能窗计算岩石密度的方法 | |
Zuo et al. | Quantitative mapping of thickness variations along a ray path using geometrical full waveform inversion and guided wave mode conversion | |
CN100495022C (zh) | 混凝土超声层析成像算法 | |
Mostafanejad et al. | Velocity structure of the northern Mississippi Embayment sediments, Part I: Teleseismic P‐wave spectral ratios analysis | |
Gaswint et al. | New results on angular reconstruction of deep pulser radio signals | |
CN115508879A (zh) | 一种海水放射性检测的伽马能谱全谱分析方法 | |
Meléndez et al. | Anisotropic P-wave travel-time tomography implementing Thomsen's weak approximation in TOMO3D | |
Zhang et al. | Illumination compensation for the adjoint-state first-arrival traveltime tomography in VTI media |
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 |