CN107942386B - 一种计算槽波地震数据激发延迟时间的方法 - Google Patents
一种计算槽波地震数据激发延迟时间的方法 Download PDFInfo
- Publication number
- CN107942386B CN107942386B CN201711132141.9A CN201711132141A CN107942386B CN 107942386 B CN107942386 B CN 107942386B CN 201711132141 A CN201711132141 A CN 201711132141A CN 107942386 B CN107942386 B CN 107942386B
- Authority
- CN
- China
- Prior art keywords
- wave
- delay time
- seismic data
- data
- refracted
- 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
- 238000000034 method Methods 0.000 title claims abstract description 42
- 230000005284 excitation Effects 0.000 claims abstract description 44
- 238000012937 correction Methods 0.000 claims abstract description 12
- 239000011435 rock Substances 0.000 claims description 22
- 230000005540 biological transmission Effects 0.000 claims description 11
- 238000004364 calculation method Methods 0.000 claims description 4
- 238000005516 engineering process Methods 0.000 claims description 4
- 238000000926 separation method Methods 0.000 claims description 4
- 238000004458 analytical method Methods 0.000 abstract description 4
- 238000012545 processing Methods 0.000 abstract description 3
- 239000003245 coal Substances 0.000 description 40
- 230000003111 delayed effect Effects 0.000 description 6
- 239000012141 concentrate Substances 0.000 description 5
- 239000002360 explosive Substances 0.000 description 4
- 238000010586 diagram Methods 0.000 description 3
- 239000011159 matrix material Substances 0.000 description 3
- 230000003044 adaptive effect Effects 0.000 description 2
- 238000007664 blowing Methods 0.000 description 2
- 238000010276 construction Methods 0.000 description 2
- 238000004880 explosion Methods 0.000 description 2
- 238000000605 extraction Methods 0.000 description 2
- 230000010287 polarization Effects 0.000 description 2
- 238000003969 polarography Methods 0.000 description 2
- 230000000644 propagated effect Effects 0.000 description 2
- 238000003860 storage Methods 0.000 description 2
- 238000009412 basement excavation Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000033228 biological regulation Effects 0.000 description 1
- 239000002817 coal dust Substances 0.000 description 1
- 230000006854 communication Effects 0.000 description 1
- 230000001934 delay Effects 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 238000005553 drilling Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000005611 electricity Effects 0.000 description 1
- 239000004744 fabric Substances 0.000 description 1
- 230000000977 initiatory effect Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- VIKNJXKGJWUCNN-XGXHKTLJSA-N norethisterone Chemical compound O=C1CC[C@@H]2[C@H]3CC[C@](C)([C@](CC4)(O)C#C)[C@@H]4[C@@H]3CCC2=C1 VIKNJXKGJWUCNN-XGXHKTLJSA-N 0.000 description 1
- 239000002245 particle Substances 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000003325 tomography Methods 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/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/30—Analysis
- G01V1/308—Time lapse or 4D effects, e.g. production related effects to the formation
-
- 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/36—Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/61—Analysis by combining or comparing a seismic data set with other data
- G01V2210/612—Previously recorded data, e.g. time-lapse or 4D
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Life Sciences & Earth Sciences (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Fluid Mechanics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明涉及槽波地震数据处理领域,具体涉及一种计算槽波地震数据激发延迟时间的方法。本发明通过分析槽波地震数据中折射波的传播特征,将折射波拾取出来并选取相应的数据点,根据选取的数据点从而计算得到震源激发延迟时间。本发明可以准确得到震源的激发延迟时间,为延迟校正提供了依据。
Description
技术领域
本发明涉及槽波地震数据处理领域,具体涉及一种计算槽波地震数据激发延迟时间的方法。
背景技术
自从德国DMT公司2010年将SUMMIT槽波地震仪引入中国后,槽波勘探技术再次引起煤矿广泛关注,在煤层厚度和断层构造等方面取得很好的探测效果。地震数据采集是槽波勘探最基础,也是最重要的环节,采集数据时,需要用***引爆安置在巷道煤层钻孔中的***,***产生振动,形成地震波在煤层中传播,被煤层中检波器接收并记录下来。在***瞬时激发引爆***的同时,槽波地震仪立即开始启动,记录下来的信号能够准确反映地震波在煤层中传播的时间;当***并不是瞬时起爆,有一定延迟时间,则槽波地震仪记录的信号传播时间比真实时间要长,这会对后期的资料处理造成错误结果。因此,当所用***不是瞬时起爆,必须确定激发延迟时间,以便对地震数据进行延迟校正。
发明内容
本发明的目的是提供一种计算槽波地震数据激发延迟时间的方法,用以解决现有技术无法计算震源激发延迟时间导致不能对地震数据进行延迟校正的问题。
为实现上述目的,本发明的方案包括一种计算槽波地震数据激发延迟时间的方法,包括以下步骤:
采集槽波地震数据并拾取所述槽波地震数据中的初始折射波数据;
利用坐标旋转和波场分离技术对所述初始折射波数据进行拾取,获得精确折射波数据;
在所述精确折射波数据中拾取设定数量的数据点,每个所述数据点包括折射波的旅行时和对应的炮检距;
对所述数据点进行拟合计算得到围岩速度和截距;
根据所述围岩速度和截距计算得到延迟时间。
进一步的,所述延迟时间ΔT的计算公式为:
其中,Vw是围岩速度,C是截距,n为所述设定数量,Ti、Di和Wi分别为第i个所述数据点对应的旅行时、炮检距和数据权值。
进一步的,所述计算围岩速度和截距的公式为:
其中,和分别为n个所述数据点对应的n个炮检距数据的平均值和n个旅行时的平均值。
进一步的,对所述围岩速度和截距进行处理得到折射波时距曲线,计算公式为:
其中,TC为所述折射波时距曲线,D为炮检距,C为截距。
进一步的,使用透射法或者折射法获取所述槽波地震数据。
进一步的,将所述延迟时间与不同类型的标准延期时间比较,得到对应的校正时间。
进一步的,所述标准延期时间分为一段0毫秒、二段25毫秒、三段50毫秒、四段75毫秒和五段100毫秒共五个类型。
本发明的有益效果是:通过槽波地震数据中折射波的传播特征,将折射波拾取出来并选取相应的数据点,根据选取的数据点从而计算得到震源激发延迟时间。本发明可以准确得到震源的激发延迟时间,为延迟校正提供了依据。
同时本发明设定了不同类型的标准延期时间,将延迟时间与标准延期时间进行比较,进而得到对应的校正时间。
附图说明
图1是煤层激发产生的地震波类型;
图2是典型槽波地震炮集记录图;
图3是煤层折射波传播剖面示意图;
图4是透射法煤层折射波平面示意图;
图5是反射法煤层折射波平面示意图;
图6是波场分离后的折射波信号图;
图7是瞬时激发与延时激发的折射波时距曲线图;
图8是本发明实施例的槽波数据采集观测***图;
图9是本发明实施例的地震槽波炮集中拾取的折射波数据图。
具体实施方式
下面结合附图对本发明做进一步详细的说明。
《煤矿安全规程》中明确规定,为防止煤矿瓦斯和煤尘***危险,在矿井采掘工作面放炮,必须使用煤矿许用瞬发电***或煤矿许用毫秒延期电***起爆。瞬发电***在电能作用下,从通电到起爆时间不大于130ms,一般为40-70ms,瞬时起爆的均一性取决于电***内的全电阻和桥丝电阻;毫秒延期电***的迟发间隔时间通过调整***内缓燃剂的化学成份和数量得到。
增加***激发延迟时间,使延时小于点燃瓦斯的***时间,可以防止放炮引起的瓦斯事故。因此,目前国内煤矿多数采用毫秒延期电***,共有五种类型,分别为一段、二段、三段、四段和五段,对应延迟时间为0、25、50、75和100毫秒。并且规定煤矿用电***总延期时间不能大于130ms(即不能使用大于五段的电***)。在进行槽波地震数据采集时,如果混合使用不同段的延期电***,将会出现激发延迟现象,必须进行延迟校正。
如图1所示,包括震源和检波器,在***引爆煤层中***后,由于煤层和围岩的波阻抗夹层组合,在煤层中产生多种类型地震波,包括槽波、直达波、透射波、滑行波和折射波等。煤矿生产上主要利用槽波,采用透射法或反射法的观测***进行数据采集,通过分析槽波频散或者包络特征进行地质勘探。对于其它类型地震波均为干扰波,研究程度相对较弱。
检波器采集得到的典型槽波地震炮集记录如图2所示,其中左半边地震道数据为检波器记录的X分量,右半边地震道数据为对应检波器接收到的Y分量,炮集记录中能够明显识别出三个波组,依据地震波传播原理,可知折射波传播时间最短,检波器最先接收到的初至波就是是折射波,对应波组①,而槽波转播速度最慢,时间最长,最后被记录,也就是图中波组③。
如图3所示,由于折射波传播速度最快,地震信号最先起跳,从地震记录中识别最容易,本次发明主要利用折射波的传播规律进行震源延迟时间的计算。当***引爆***,激发产生地震波,入射角度达到临界角后,地震波将会沿着煤层和围岩的分界面传播,形成滑行波,根据惠更斯原理,滑行波传播过程中,将会形成折射波返回煤层,被煤层中检波器所接收,所以,折射波在煤层中传播路径为图中所示OABR;图中还包括煤厚h、围岩速度VW、煤层速度VC、炮检距D和临界角θ。
如图4所示为透射法采集地震数据的示意图。透射法采集方式是将炮点和检波器布置在工作面不同巷道内,接收工作面内的透射地震信号,常用于对煤层厚度进行探测。采集得到的炮集地震记录中,由于震源和检波器不在一条直线上,槽波地震炮集记录中(如图2所示)横轴并不表示炮检距D,而是震源O在检波器排列线上的投影O’与检波器R之间的距离x。
因此,槽波地震炮集记录中折射波的传播时间与传播距离,也就是时距关系为:
公式中VW、VC分别为围岩速度和煤层速度,h为煤层厚度,L为工作面宽度,D为震源到检波器距离,C是由煤层厚度、煤层速度和围岩速度确定的常数,由于震源和检波器排列不在一条直线上,槽波地震炮集记录中的折射波时距曲线和地面地震勘探中折射波时距曲线特征不同,表现为双曲线形态,双曲线的顶点即最小值位于震源在检波器排列线上的投影处,随着x增大,双曲线变陡,围岩速度越大,曲线越缓。
图5所示是反射法采集地震数据的示意图。反射法采集方式是将炮点和检波器布置在矿井同一巷道内,排列方式与常规地面地震勘探数据采集方式类似,接收来自工作面内的反射地震信号,适用于对煤层内地质构造(断层、陷落柱等)进行探测,对应折射波时距方程为:
公式中变量含义与折射波时距公式(1)中变量一致。可以看出,由于反射法中震源和检波器排列在一条直线上,槽波地震炮集记录中折射波时距曲线为直线,与常规地面勘探中折射波时距曲线特征一致。
采集的槽波地震数据中,包含多种类型地震波,为提高分析精度,需要准确提取出其中的折射波信号。提取槽波地震数据中的折射波共包含两个步骤,一是进行坐标旋转;二是进行波场分离。以目前煤矿使用较多的德国SUMMIT槽波地震仪为例进行说明,该仪器采用水平双分量检波器,垂直于煤壁水平安置在巷道煤层中。由于各检波器X分量和Y分量接收到的信号传播方向不统一,在两个分量上均记录有折射波信号(如图2),这就需要进行坐标旋转,将X和Y分量校正到水平平行波传播方向和水平垂直波传播方向。为获得坐标旋转角度θ,依据纵波质点振动方向与波传播方向一致,也就是纵波在波传播方向能量最大,在波传播垂直方向能量最小的特点,可求得旋转角度θ为:
公式中xi、yi分别为槽波地震炮集中折射波数据点在X分量和Y分量上记录的数据,n为拾取的折射波信号数据点个数。
首先进行第一次折射波旅行时拾取,分别在X分量和Y分量上确定折射波旅行时的起始和终止位置,依据公式(3)计算需要旋转的角度。经过坐标旋转后,原X分量存储水平平行分量,记录水平面内波振动方向与传播方向相同的地震信号;原Y分量存储水平垂直分量,记录水平面内波振动方向和传播方向垂直的地震信号。
然后采用基于S变换的自适应协方差矩阵极化分析方法进行折射波提取,通过得到的矩阵特征值和特征向量计算极化滤波器中椭圆率和方位角参数,由于噪音等干扰,设定椭圆率小于某一特定值(如0.2)时近似认为线性极化,方位角在0°到45°之间,近似认为平行波传播方向,分离得到的信号近似认为是线性极化,且振动方向与波传播方向平行,即折射波,结果如图6所示。由于折射纵波(波组①)振动方向与传播方向一致,在水平平行分量上,折射纵波能量最强,同相轴清晰明显,而在水平垂直分量上,初至折射纵波(波组①)能量最弱,几乎看不到。这与折射波的传播规律一致,说明该方法提取的折射波准确可信。
通过槽波地震数据中折射波的时距特征,可以看出不管是采用透射法还是反射法观测***,折射波时距曲线为一条不过原点的直线,其斜率为围岩速度VW,截距为由煤层厚度、煤层速度和围岩速度确定的常数C。只要知道VW和C就可确定折射波的时距曲线特征。但不同的采集方式得到的地震炮集记录中,横轴含义不同。在反射法中,横轴与常规地面地震勘探得到的地震炮集记录含义相同,折射波表现为一条直线;在透射法中,横轴为震源在检波器排列方向上的投影,炮集记录中折射波表现为一条双曲线。
在提取的水平平行分量地震炮集记录上,进行第二次折射波旅行时拾取(图6中红色短线),假定拾取了n组折射波旅行时数据(T1,D1),(T2,D2),…,(Tn,Dn)。通过误差最小平方和确定最佳匹配函数,利用最小二乘法拟合槽波地震炮集中折射波时距曲线,经过推导具体公式为:
公式中Ti和Di分别为炮集中拾取折射波的旅行时和炮检距,和为n个旅行时和炮检距数据的平均值。由此可得到围岩速度VW和截距C,即可确定对应折射波的时距曲线。
当采用毫秒延期电***起爆***时,从槽波地震仪器开始记录到***起爆有一定时间延迟。由上述折射波时距特征可以看出,***延时激发采集的地震炮集记录中,折射波时距曲线的斜率为围岩速度,保持不变,只是传播时间滞后一段时间。因此,震源延时激发产生的折射波与瞬时激发产生的折射波在时距曲线上表现为截距不同的平行直线,如图7所示。
因此,毫秒延期电***激发延迟时间为折射波时距曲线与瞬时激发折射波时距曲线之间的水平距离。只要确定瞬时激发对应的折射波时距曲线,那么***延时激发的延时时间ΔT为:
公式中Wi为数据权值,根据拾取初至波旅行时数据点的可信度人为设定,数值在0到1之间,数据点越可信,数值越接近1。Ti和Di分别为延时激发采集炮集中地震道拾取折射波的旅行时和炮检距,VW和C为瞬时激发炮集中折射波时距曲线拟合得到的围岩速度和截距,n为拾取初至波数据点的个数。
考虑到拾取折射波旅行时数据点的人工误差,以及仪器计时、***质量等误差,对计算得到***延时结果进行校正,如表1所示,当延时ΔT小于25ms时,认为采用一段延期***激发,不进行校正;当延时ΔT大于25ms且小于50ms时,认为采用两段延期***激发,延时校正为25ms;当延时ΔT大于50ms且小于75ms时,认为采用三段延期***激发,延时校正为50ms;当延时ΔT大于75ms且小于100ms时,认为采用四段延期***激发,延时校正为75ms;当延时ΔT大于100ms,认为采用五段延期***激发,延时校正为100ms。如果延时ΔT大于130ms时,则认为计算结果有问题。
表1 震源激发延时校正
下面给出一个具体的实施例来说明本发明所述方法。
某煤矿工作面进行槽波地震数据采集,设置震源48炮,检波器47道,采用透射法进行勘探,观测***排列见图8。在施工过程中,由于放炮人员没有严格遵循施工要求,以及客观条件限制等因素,震源激发混合使用了不同段的***,造成采集数据存在激发延迟现象,需要确定炮集记录对应的延迟时间。
在采集的槽波地震炮集记录中,首先在X分量和Y分量上进行第一次折射波旅行时数据拾取,拾取时不需要精准确定折射波的起始和终止点位置,只要包含折射波即可,依据公式(3)获得旋转角度,进行坐标旋转,将地震波信号校正到平行波传播分量和垂直分量上。
然后利用基于S变换的自适应协方差矩阵极化分析方法提取折射波,在水平平行分量上进行第二次折射波旅行时提取。拾取过程中,只需要拾取折射波最清晰,不存在争议的地震道,需要注意透射法采集炮集记录中折射波旅行时表现为双曲线特征。在时间炮检距图中显示拾取的折射波旅行时数据,如图9所示,相同颜色的数据点为同一炮拾取的折射波旅行时。
从数据分布特征可以明显看出,槽波炮集中部分炮集记录明显存在激发延迟现象。越靠近左侧则激发延时时间越短,因此一般选取最左侧的数据点作为瞬时激发的(延时为0ms)基准,本实施例中从最左侧折射波旅行时数据中,选择位于中间位置的第3炮激发的折射波旅行时作为基准,认为该炮为瞬时激发或一段***激发(延时为0ms),依据公式(4)和(5),拟合出折射波时距曲线对应的围岩速度和截距,并相应得到煤矿现用五种类型毫秒延期电***激发方式中另外四种激发方式所对应的折射波时距曲线,如图9中直线所示,分别对应延迟时间为0ms、25ms、50ms、75ms和100ms。通过公式(6)分别计算各炮的激发延迟时间,经过计算发现第12、13、15、16、28炮的延时ΔT分别为77ms、82ms、79ms、80ms、83ms。依据表1校正数据,最终认定第12、13、15、16、28炮采用的是四段延期***激发,需要进行75ms延迟校正处理。
以上给出了本发明涉及的具体实施方式,通过分析槽波地震数据中折射波的传播特征,将折射波拾取出来并选取相应的数据点,根据选取的数据点从而计算得到震源激发延迟时间,同时设定了不同类型的标准延期时间,将延迟时间与标准延期时间进行比较,进而得到对应的校正时间。
但本发明不局限于所描述的实施方式,例如采集槽波地震数据时方法的具体选取,或者拾取折射波时的具体手段的改变,这样形成的技术方案是对上述实施例进行微调形成的,这种技术方案仍落入本发明的保护范围内。
Claims (6)
1.一种计算槽波地震数据激发延迟时间的方法,其特征在于,包括以下步骤:
采集槽波地震数据并拾取所述槽波地震数据中的初始折射波数据;
利用坐标旋转和波场分离技术对所述初始折射波数据进行拾取,获得精确折射波数据;
在所述精确折射波数据中拾取设定数量的数据点,每个所述数据点包括折射波的旅行时和对应的炮检距;
对所述数据点进行拟合计算得到围岩速度和截距;
根据所述围岩速度和截距计算得到延迟时间;
所述延迟时间ΔT的计算公式为:
其中,Vw是围岩速度,C是截距,n为所述设定数量,Ti、Di和Wi分别为第i个所述数据点对应的旅行时、炮检距和数据权值。
2.根据权利要求1所述一种计算槽波地震数据激发延迟时间的方法,其特征在于:所述计算围岩速度和截距的公式为:
其中,和分别为n个所述数据点对应的n个炮检距数据的平均值和n个旅行时的平均值。
3.根据权利要求2所述一种计算槽波地震数据激发延迟时间的方法,其特征在于:对所述围岩速度和截距进行处理得到折射波时距曲线,计算公式为:
其中,TC为所述折射波时距曲线,D为炮检距,C为截距。
4.根据权利要求3所述一种计算槽波地震数据激发延迟时间的方法,其特征在于:使用透射法或者反射法获取所述槽波地震数据。
5.根据权利要求1所述一种计算槽波地震数据激发延迟时间的方法,其特征在于:将所述延迟时间与不同类型的标准延期时间比较,得到对应的校正时间。
6.根据权利要求5所述一种计算槽波地震数据激发延迟时间的方法,其特征在于:所述标准延期时间分为一段0毫秒、二段25毫秒、三段50毫秒、四段75毫秒和五段100毫秒共五个类型。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711132141.9A CN107942386B (zh) | 2017-11-15 | 2017-11-15 | 一种计算槽波地震数据激发延迟时间的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711132141.9A CN107942386B (zh) | 2017-11-15 | 2017-11-15 | 一种计算槽波地震数据激发延迟时间的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107942386A CN107942386A (zh) | 2018-04-20 |
CN107942386B true CN107942386B (zh) | 2019-05-07 |
Family
ID=61932349
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711132141.9A Active CN107942386B (zh) | 2017-11-15 | 2017-11-15 | 一种计算槽波地震数据激发延迟时间的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107942386B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112230285B (zh) * | 2019-07-15 | 2023-09-26 | 中国石油天然气集团有限公司 | 地震数据激发延迟时间的校正方法及校正装置 |
CN111208557A (zh) * | 2020-03-03 | 2020-05-29 | 何良 | 井下分布式槽波地震仪智能化数据回收和校正方法 |
CN111650636A (zh) * | 2020-04-28 | 2020-09-11 | 徐州天浩勘探有限公司 | 一种煤层底板隐伏构造探测方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO1992012444A1 (en) * | 1991-01-02 | 1992-07-23 | Western Atlas International, Inc. | Method for determining geometry of subsurface features while drilling |
CN102788991A (zh) * | 2012-07-16 | 2012-11-21 | 中煤科工集团西安研究院 | 基于z分量的透射槽波高密度快速探测方法 |
CN105785440A (zh) * | 2016-02-29 | 2016-07-20 | 河南理工大学 | 一种矿井槽波双分量地震信号频散曲线提取方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US10353092B2 (en) * | 2015-12-10 | 2019-07-16 | Pgs Geophysical As | Velocity model update with an inversion gradient |
-
2017
- 2017-11-15 CN CN201711132141.9A patent/CN107942386B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO1992012444A1 (en) * | 1991-01-02 | 1992-07-23 | Western Atlas International, Inc. | Method for determining geometry of subsurface features while drilling |
CN102788991A (zh) * | 2012-07-16 | 2012-11-21 | 中煤科工集团西安研究院 | 基于z分量的透射槽波高密度快速探测方法 |
CN105785440A (zh) * | 2016-02-29 | 2016-07-20 | 河南理工大学 | 一种矿井槽波双分量地震信号频散曲线提取方法 |
Non-Patent Citations (2)
Title |
---|
槽波地震反射法在断裂构造探测中的应用;姚小帅 等;《中州煤炭》;20151231(第9期);第101-104页 * |
透射法地震勘探中槽波几何运动学特征;杜艳艳 等;《地球物理学进展》;20171031;第32卷(第5期);第1978-1983页 * |
Also Published As
Publication number | Publication date |
---|---|
CN107942386A (zh) | 2018-04-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107942386B (zh) | 一种计算槽波地震数据激发延迟时间的方法 | |
Zhao et al. | Tunnel seismic tomography method for geological prediction and its application | |
CN104459797B (zh) | 一种井中微地震事件识别和拾取的方法 | |
CN107642114B (zh) | 桩基浇注前桩底隐患探查方法及其探查装置 | |
Malehmir et al. | Landstreamer seismics and physical property measurements in the Siilinjärvi open-pit apatite (phosphate) mine, central Finland | |
CN102262240B (zh) | 三维观测***满覆盖区域自动布设方法 | |
CN105093297B (zh) | 一种微地震定位精度的质量控制方法 | |
CN100349006C (zh) | 一种用波阻抗反演技术确定岩性和流体分布的方法 | |
CN109799529A (zh) | 一种基于互相关的横波***vsp裂缝预测方法 | |
Nakamura et al. | Preliminary results of a large offset seismic survey west of Hengchun Peninsula, southern Taiwan | |
CN105467445B (zh) | 一种三维地震观测***的建立方法及装置 | |
Ma | Multi-dipole observation system and study on the abnormal variation of the geoelectric field observed at Capital Circle area before the Wen’an MS 5.1 earthquake | |
Poulter | The Poulter seismic method of geophysical exploration | |
CN108732618A (zh) | 一种用于地震监视记录上快速识别横波的装置和方法 | |
Li | Viability of azimuthal variation in P-wave moveout for fracture detection | |
Bennett | 3-D seismic refraction for deep exploration targets | |
Hall et al. | Seismic velocities in the Lewisian metamorphic complex, northwest Britain—‘in situ’measurements | |
CN109521467A (zh) | 一种基于突出煤层巷道的超前探测方法 | |
Zhao et al. | Application of improved ant-tracking method in carbonate-fracture reservoir | |
Mintrop | On the stratification of the Earth’s crust according to seismic studies of a large explosion and of earthquakes | |
Kay et al. | Imaging the Moho and Vp/Vs ratio in the western Superior Archean craton with wide angle reflections | |
Maunde et al. | Seismic refraction investigation of fracture zones and bedrock configuration for geohydrologic and geotechnical studies in part of Nigeria’s Capital City, Abuja | |
Menke et al. | Observations of the transition from linear polarization to complex polarization in short-period compressional waves | |
Larroque et al. | How 3D seismic can help enhance mining | |
Corrigan et al. | Estimate of shear-wave anisotropy using multicomponent seismic data |
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 |