CN106680868B - 一种槽波地震多属性参数层析方法 - Google Patents
一种槽波地震多属性参数层析方法 Download PDFInfo
- Publication number
- CN106680868B CN106680868B CN201710008578.5A CN201710008578A CN106680868B CN 106680868 B CN106680868 B CN 106680868B CN 201710008578 A CN201710008578 A CN 201710008578A CN 106680868 B CN106680868 B CN 106680868B
- Authority
- CN
- China
- Prior art keywords
- wave
- frequency
- energy
- tomography
- slot
- 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
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. analysis, for interpretation, for correction
- G01V1/282—Application of seismic models, synthetic seismograms
-
- 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. analysis, for interpretation, for correction
- G01V1/30—Analysis
- G01V1/301—Analysis for determining seismic cross-sections or geostructures
-
- 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. analysis, for interpretation, for correction
- G01V1/30—Analysis
- G01V1/306—Analysis for determining physical properties of the subsurface, e.g. impedance, porosity or attenuation profiles
-
- 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. analysis, for interpretation, for correction
- G01V1/36—Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
- G01V1/362—Effecting static or dynamic corrections; Stacking
-
- 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/62—Physical property of subsurface
- G01V2210/622—Velocity, density or impedance
- G01V2210/6222—Velocity; travel time
-
- 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/62—Physical property of subsurface
- G01V2210/624—Reservoir parameters
-
- 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/64—Geostructures, e.g. in 3D data cubes
Abstract
本发明公开了一种槽波地震多属性参数层析方法,方法包括如下步骤:1)理论模型分析2)地震数据预处理;3)参数的拾取方法;4)能量、频率矫正算法;5)层析成像。本发明的优点是:运用槽波、P波、S波的时间反演速度,槽波、P波、S波的能量反演断层、槽波的频率反演煤厚,并综合对比分析,解决工作面内构造、煤厚变化、顶底板破碎带等问题效果显著。
Description
技术领域
本发明涉及一种煤矿井下槽波地震勘探方法,更具体的说,是涉及一种槽波地震多属性参数层析方法。
背景技术
目前,槽波地震数据层析成像主要运用槽波的时间、P波的时间参数进行速度成像,通过时间参数反算煤层中槽波或P波的速度分布情况,从而反演地质构造异等。此方法在构造较为简单的条件下可取得良好的效果,但是在构造稍复杂的情况下,例如大断层可阻断槽波或P波的传播,从而使检波器接收不到槽波或P波,导致无法用时间参数来反算速度。而且,速度参数并非槽波对异常体反应最敏感的参数。
发明内容
本发明针对现有技术的不足,而提供一种运用槽波、P波、S波的时间反演速度,槽波、P波、S波的能量反演断层、槽波的频率反演煤厚,并综合对比分析,用于指导煤矿工作面内构造、煤厚变化、顶底板破碎带等问题的一种槽波地震多属性参数层析方法。
本发明的一种槽波地震多属性参数层析方法,方法包括如下步骤:
1)理论模型分析;
依据工作面内顶板、煤层、夹矸、底板的密度和速度等参数建立地层模型,运用槽波频散方程计算频散曲线,从而确立槽波埃里相位的速度和频率参数,为槽波的处理提供技术参数;
2)地震数据预处理;
地震数据预处理包括常规预处理和特殊预处理;常规预处理和地面地震波预处理类似,包括建立观测系统、剔除坏道、极性反转、一维频率域滤波、频谱分析,能量分析等;特殊预处理主要针对槽波进行的处理,包括柱面扩散矫正、极化旋转、两分量的旋转与合成、频散分析、时频分析、速度分析等;
3)参数的拾取方法;
多属性参数层析技术主要运用槽波、P波、S波的时间参数来反演速度,运用槽波、P波、S波的能量反演断层,运用槽波的频率反演煤厚,先提取出这些参数数据;实测地震记录上往往同时存在P波、S波和槽波,采取开时窗对三种类型的波进行手动参数拾取,参数包括时间、能量、频率三类;时间为拾取波形的初至时间,能量为时窗内记录包络计算求和并进行地层衰减补偿后的真振幅能量值,频率值为时窗内的计算进行傅里叶变换求取的主频值;
4)能量、频率矫正算法;
槽波数据随炮检距的变化呈指数衰减,而且受到激发药量、检波器耦合等因素的影响,使各道的槽波能量呈现出数量级的差异,构造引起的槽波能量变化难以体现出来,也无法直接进行层析成像,采用能量、频率矫正算法,矫正掉炮检距、激发药量、检波器耦合等因素造成的槽波能量差异,矫正后的能量参数可直接用于层析成像;
5)层析成像;
层析成像采用直射线正演和SIRT反演算法;层析成像算法包括走时层析成像算法,能量、频率层析成像算法。
作为优选进一步的,所述步骤(4)中所述能量、频率矫正算法具体分如下步骤:
a)对每一道地震数据做振幅谱,从0-500hz,每间隔5Hz取一个频点,共101个;
b)取出所有道的第一个频点,按炮检距大小排序,对排序后的数据做线性回归分析,然后回归线拉水平,频点跟着回归线一起上提;
c)同理,取出所有道的第二个、第三个频点,直到做完所有频点;
d)对于每一道,分别取出它校正后的101个频点中最大值,作为该道校正后的能量值,即完成能量矫正,矫正后的数据在0~1之间;
e)频率矫正和能量矫正方法相同。
作为优选进一步的,所述步骤(5)中的走时层析成像算法,包括以下步骤:
a)建立工作面初始模型,划分计算网格;
b)直射线追踪,找出地震波从炮点到接收点的传播路径,并算出理论走时;
c)根据理论走时与实际拾取走时的差值来反演求解,对模型进行修正;
d)对修正后的模型进行插值处理;
e)对修正后的模型进行平滑;
f)重复上述b~e四个步骤,直到模型修正满足一定的要求。
作为优选进一步的,所述步骤(5)中能量和频率层析成像算法:
a)建立工作面初始模型,划分计算网格;
b)直射线追踪,能量和频率首先乘以对应的炮检距,然后按照射线长度分配到各个网格中;
c)根据理论能量、频率与实际拾取的能量、频率差值来反演求解,对模型进行修正;
d)对修正后的模型进行插值处理;
e)对修正后的模型进行平滑;
f)重复上述b~e四个步骤,直到模型修正满足一定的要求。
本发明的有益效果是:和现有技术相比,运用槽波的时间、P波的时间参数进行速度层析成像,在构造稍复杂的情况下,大断层可阻断槽波或P波的传播,从而使检波器接收不到槽波或P波,导致用时间参数来反算速度不准确。对于煤矿生产,构造是首要问题,速度参数并非对构造反应最敏感的参数,实践表明,对构造反应最为敏感的参数是槽波埃里相能量,对煤厚反应敏感的参数是槽波的频率,此外直达P波、S波的能量、频率、速度也有一定的反应。基于此,本发明实现多属性参数层析技术,即运用槽波、P波、S波的时间反演速度,槽波、P波、S波的能量反演断层、槽波的频率反演煤厚,并综合对比分析,解决工作面内构造、煤厚变化、顶底板破碎带等问题效果显著。
附图说明
图1为理论槽波频散曲线;
图2为预处理后单炮记录;
图3为槽波能量层析成像结果;
图4为槽波频率层析成像结果;
图5为CT成像结果的地质解释结果图。
具体实施方式
下面通过具体实施例,对本发明的技术方案作进一步具体的说明,但是本发明并不限于实施例。
实施例1某个煤矿工作面槽波勘探实例:
(1)地质任务;
煤层内的断层等构造给煤矿安全生产带来一定的隐患,因此需要探明工作面内的地质异常区域,为工作面回采提供地质保障。本次工作面槽波地震探测的地质任务如下:查明工作面内落差大于1/2煤厚的断层发育情况、查明工作面内长轴直径大于20m的陷落柱、对其他地质异常体进行解释。
(2)工作面概况;
煤矿工作面走向长1938m,倾斜长219.7m。该工作面主采煤层为3#煤层,煤层厚度一般为4.2m-5.0m,平均煤厚为4.6m,煤层结构简单,煤层倾角为1~8°,平均5°。工作面东侧为5103巷和5105巷,且这两条巷道正在掘进。北侧为5303工作面,正在掘进。
工作面由53032巷、切眼、53162巷、53162联络巷组成,槽波勘探区域为该4条巷道包围区域(停采线以外不进行探测,观测系统布置延伸至停采线以外50m)。
5316工作面53032巷长度1317m;5316切眼探测长度为273m,53162巷探测长度为1317m。探测巷道总长度为3180m,探测区域面积为359541m2。
(3)槽波特征的理论分析;
根据测区煤层及围岩的密度、P波和S波速度、煤层厚度及结构等资料,从理论上计算槽波的频散曲线及振幅分布曲线,进一步合成槽波记录,预测槽波尤其是埃里震相的速度、频散特征与振幅分布等,为设计和施工、资料处理及分析提供理论指导。
如图1所示,根据该矿实际资料确定该工作面为对称三层模型,参数如下:顶底板横波速度为1800m/s、密度为2.6g/cm3,煤层厚度为5m、密度为1.2g/cm3、横波速度为900m/s。考虑Love型槽波,计算结果如图1所示。从中可以看出本工作面槽波的埃里相速度约为800m/s,频率约为140hz,煤层中心能量最强。
(4)透射槽波的原始数据预处理;
如图2所示,包括建立观测系统、剔除坏道、极性反转、一维频率域滤波、频谱分析,能量分析、柱面扩散矫正、极化旋转、两分量的旋转与合成、频散分析、时频分析、速度分析等处理步骤。预处理后单炮记录。
(5)层析成像
选用槽波的能量、频率参数进行反演构造和煤厚变化,层析成像采用直射线正演和SIRT反演算法。成像结果见图3、图4。
(6)地质成果
从图3中可以看出颜色从浅色到深色代表衰减的程度,颜色越深表示槽波能量被吸收得越严重,该区域往往有断层或陷落柱发育区域。异常区具有一定的形态,不同的形态往往代表着不同的地质现象,例如断层一般为“条带状”异常,而陷落柱的异常则呈“片状”分布。衰减的大小则用于判定异常体的定量参数如断层落差等。需要说明的是,槽波资料的解释离不开已知地质资料如巷道揭露资料的标定,透射槽波往往只能判定断层或陷落柱存在的位置,对于断层的倾向和落差大小还是要与巷道揭露的异常体参数进行对比方可得出,本次解释中有的断层完全隐伏于工作面内,其倾向就无法确定,而落差只能根据其延展长度大致推断。
图5为CT成像结果的地质解释结果图,图中断层以CF表示,陷落柱以CX表示。
通过对所有异常体的分析解释,本次勘探共发现断层16条,除了CF5落差较大以外,其它断层都为落差2米左右,或1.5米左右的小断层。通过分析共发现陷落柱4个,分别为CX1、CX2、CX3、CX4,它们长轴分别约为209m、33m、68m和33m。异常体的信息情况见表1。在工作面西北角有一异常被定义为资料不可靠区域,未做地质解释。值得注意的是CF5断层落差较大,其异常的影响范围也比较大,推测其附近可能还发育其它断层或陷落柱等。
表1槽波勘探成果解释汇总表
以上所述仅为本发明的实施例,并非因此限制本发明的专利范围,凡是利用本发明说明书内容所作的等效结构或等效流程变换,或直接或间接运用在其他相关的技术领域,均同理包括在本发明的专利保护范围内。
Claims (3)
1.一种槽波地震多属性参数层析方法,其特征在于,方法包括如下步骤:
1)理论模型分析;
依据工作面内顶板、煤层、夹矸、底板的密度和速度参数建立地层模型,运用槽波频散方程计算频散曲线,从而确立槽波埃里相位的速度和频率参数,为槽波的处理提供技术参数;
2)地震数据预处理;
地震数据预处理包括常规预处理和特殊预处理;常规预处理和地面地震波预处理类似,包括建立观测系统、剔除坏道、极性反转、一维频率域滤波、频谱分析,能量分析;特殊预处理针对槽波进行的处理,包括柱面扩散矫正、极化旋转、两分量的旋转与合成、频散分析、时频分析、速度分析;
3)参数的拾取方法;
多属性参数层析技术运用槽波、P波、S波的时间参数来反演速度,运用槽波、P波、S波的能量反演断层,运用槽波的频率反演煤厚,先提取出这些参数数据;采取开时窗对三种类型的波进行手动参数拾取,参数包括时间、能量、频率三类;时间为拾取波形的初至时间,能量为时窗内记录包络计算求和并进行地层衰减补偿后的真振幅能量值,频率值为时窗内的计算进行傅里叶变换求取的主频值;
4)能量、频率矫正算法;
槽波数据随炮检距的变化呈指数衰减,而且受到激发药量、检波器耦合因素的影响,使各道的槽波能量呈现出数量级的差异,构造引起的槽波能量变化难以体现出来,也无法直接进行层析成像,采用能量、频率矫正算法,矫正掉炮检距、激发药量、检波器耦合因素造成的槽波能量差异,矫正后的能量参数可直接用于层析成像;
5)层析成像;
层析成像采用直射线正演和SIRT反演算法;层析成像算法包括走时层析成像算法,能量和频率层析成像算法;
所述步骤(4)中所述能量、频率矫正算法具体分如下步骤:
a)对每一道地震数据做振幅谱,从0-500hz,每间隔5Hz取一个频点,共101个;
b)取出所有道的第一个频点,按炮检距大小排序,对排序后的数据做线性回归分析,然后回归线拉水平,频点跟着回归线一起上提;
c)同理,取出所有道的第二个、第三个频点,直到做完所有频点;
d)对于每一道,分别取出它校正后的101个频点中最大值,作为该道校正后的能量值,即完成能量矫正,矫正后的数据在0~1之间;
e)频率矫正和能量矫正方法相同。
2.根据权利要求1所述的一种槽波地震多属性参数层析方法,其特征在于,所述步骤(5)中的走时层析成像算法,包括以下步骤:
a)建立工作面初始模型,划分计算网格;
b)直射线追踪,找出地震波从炮点到接收点的传播路径,并算出理论走时;
c)根据理论走时与实际拾取走时的差值来反演求解,对模型进行修正;
d)对修正后的模型进行插值处理;
e)对修正后的模型进行平滑;
f)重复上述b~e四个步骤,直到模型修正满足一定的要求。
3.根据权利要求1所述的一种槽波地震多属性参数层析方法,其特征在于,所述步骤(5)中能量和频率层析成像算法:
a)建立工作面初始模型,划分计算网格;
b)直射线追踪,能量和频率首先乘以对应的炮检距,然后按照射线长度分配到各个网格中;
c)根据理论能量、频率与实际拾取的能量、频率差值来反演求解,对模型进行修正;
d)对修正后的模型进行插值处理;
e)对修正后的模型进行平滑;
f)重复上述b~e四个步骤,直到模型修正满足一定的要求。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710008578.5A CN106680868B (zh) | 2017-01-05 | 2017-01-05 | 一种槽波地震多属性参数层析方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710008578.5A CN106680868B (zh) | 2017-01-05 | 2017-01-05 | 一种槽波地震多属性参数层析方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106680868A CN106680868A (zh) | 2017-05-17 |
CN106680868B true CN106680868B (zh) | 2019-01-15 |
Family
ID=58850245
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710008578.5A Expired - Fee Related CN106680868B (zh) | 2017-01-05 | 2017-01-05 | 一种槽波地震多属性参数层析方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106680868B (zh) |
Families Citing this family (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108931812B (zh) * | 2017-05-24 | 2019-11-01 | 河南理工大学 | 一种基于透射法槽波埃里相判识构造煤的方法 |
CN109557583B (zh) * | 2017-09-26 | 2020-12-01 | 中国石油化工股份有限公司 | 一种地震属性提取方法及系统 |
CN108761532B (zh) * | 2018-05-08 | 2021-02-19 | 河北煤炭科学研究院 | 煤矿工作面断层预测方法 |
CN109597124B (zh) * | 2018-12-29 | 2021-06-04 | 煤炭科学技术研究院有限公司 | 一种工作面应力集中区槽波地震精细探测方法 |
CN109991660B (zh) * | 2019-04-18 | 2020-07-10 | 中煤科工集团西安研究院有限公司 | 基于槽波地震和克里金插值的煤层厚度预测方法及装置 |
CN110531415B (zh) * | 2019-08-21 | 2020-10-30 | 徐州工程学院 | 一种利用围岩松动圈影响的三维小断层超前探测方法 |
CN110531419B (zh) * | 2019-08-21 | 2020-10-30 | 徐州工程学院 | 一种利用Love型面波的陷落柱超前探测方法 |
CN111650636A (zh) * | 2020-04-28 | 2020-09-11 | 徐州天浩勘探有限公司 | 一种煤层底板隐伏构造探测方法 |
CN111812714B (zh) * | 2020-06-08 | 2022-11-15 | 中煤科工集团西安研究院有限公司 | 基于折射纵波与高频槽波的煤层纵横波速度求取方法 |
CN112363210B (zh) * | 2020-11-13 | 2023-10-13 | 福州华虹智能科技股份有限公司 | 基于透射槽波波速和衰减系数联合反演的煤厚定量预测方法 |
CN112987088B (zh) * | 2021-02-22 | 2023-04-18 | 成都理工大学 | 一种渗流介质地震横波数值模拟和成像方法 |
CN113075730A (zh) * | 2021-03-23 | 2021-07-06 | 西安中地博睿探测科技有限公司 | 一种透射槽波速度ct成像中走时信息自动提取方法 |
CN116755154B (zh) * | 2023-08-14 | 2023-12-12 | 煤炭科学技术研究院有限公司 | 地质构造的识别方法、装置、电子设备和存储介质 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2009029860A1 (en) * | 2007-08-29 | 2009-03-05 | Baker Hughes Incorporated | Downhole measurements of mud acoustic velocity |
WO2010142659A2 (en) * | 2009-06-09 | 2010-12-16 | Shell Internationale Research Maatschappij B.V. | Method for stratigraphic analysis of seismic data |
CN102788991A (zh) * | 2012-07-16 | 2012-11-21 | 中煤科工集团西安研究院 | 基于z分量的透射槽波高密度快速探测方法 |
CN103837897A (zh) * | 2013-12-13 | 2014-06-04 | 山东科技大学 | 煤矿井下槽波精细探测方法 |
CN105137475A (zh) * | 2015-09-17 | 2015-12-09 | 中煤科工集团西安研究院有限公司 | 一种基于皮带机的煤矿工作面实时探测系统及方法 |
CN105353410A (zh) * | 2015-12-10 | 2016-02-24 | 北京中矿大地地球探测工程技术有限公司 | 一种煤矿井下二分量槽波地震勘探方法 |
CN105700010A (zh) * | 2016-03-02 | 2016-06-22 | 河北煤炭科学研究院 | 煤矿联合地震全息预测方法 |
-
2017
- 2017-01-05 CN CN201710008578.5A patent/CN106680868B/zh not_active Expired - Fee Related
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2009029860A1 (en) * | 2007-08-29 | 2009-03-05 | Baker Hughes Incorporated | Downhole measurements of mud acoustic velocity |
WO2010142659A2 (en) * | 2009-06-09 | 2010-12-16 | Shell Internationale Research Maatschappij B.V. | Method for stratigraphic analysis of seismic data |
CN102788991A (zh) * | 2012-07-16 | 2012-11-21 | 中煤科工集团西安研究院 | 基于z分量的透射槽波高密度快速探测方法 |
CN103837897A (zh) * | 2013-12-13 | 2014-06-04 | 山东科技大学 | 煤矿井下槽波精细探测方法 |
CN105137475A (zh) * | 2015-09-17 | 2015-12-09 | 中煤科工集团西安研究院有限公司 | 一种基于皮带机的煤矿工作面实时探测系统及方法 |
CN105353410A (zh) * | 2015-12-10 | 2016-02-24 | 北京中矿大地地球探测工程技术有限公司 | 一种煤矿井下二分量槽波地震勘探方法 |
CN105700010A (zh) * | 2016-03-02 | 2016-06-22 | 河北煤炭科学研究院 | 煤矿联合地震全息预测方法 |
Non-Patent Citations (4)
Title |
---|
"CT技术确定煤矿中小断层的应用研究";吴翔飞 等;《西部探矿工程》;20030731(第86期);第92-95页 |
"含夹矸煤层的Love型槽波频散特征";匡伟 等;《煤矿开采》;20151231;第20卷(第6期);第23-26、43页 |
"槽波与P波对煤矿冲击地压风险区联合探测应用";朱师波 等;《能源技术与管理》;20161031;第41卷(第5期);第14-16页 |
"槽波地震勘探在李雅庄煤矿构造探测中的应用";常英俊 等;《山西焦煤科技》;20160630(第6期);第27-30页 |
Also Published As
Publication number | Publication date |
---|---|
CN106680868A (zh) | 2017-05-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106680868B (zh) | 一种槽波地震多属性参数层析方法 | |
CN104570125B (zh) | 一种利用井数据提高成像速度模型精度的方法 | |
CN104635269B (zh) | 一种基于叠前方位角道集预测火成岩裂缝型储层的方法 | |
Malehmir et al. | 3D reflection seismic imaging for open-pit mine planning and deep exploration in the Kevitsa Ni-Cu-PGE deposit, northern Finland | |
CN101915939B (zh) | 一种面波压制方法 | |
CN105510880A (zh) | 一种基于双差法的微地震震源定位方法 | |
CN107045145B (zh) | 地震层序控制下的叠前振幅随偏移距变化检测缝洞方法 | |
CN103809216B (zh) | 一种电阻率数据与地震数据联合速度建场方法 | |
CN102053261A (zh) | 一种地震数据处理方法 | |
CN101738636B (zh) | 一种三维vsp高斯束法多波联合偏移成像方法 | |
CN102375154A (zh) | 基于宽方位三维地震的裂缝参数确定方法 | |
Ahmadi et al. | High-resolution 2D seismic imaging and forward modeling of a polymetallic sulfide deposit at Garpenberg, central Sweden | |
CN105093297B (zh) | 一种微地震定位精度的质量控制方法 | |
CN109765615A (zh) | 一种地层品质因子反演方法及装置 | |
CN105467445B (zh) | 一种三维地震观测系统的建立方法及装置 | |
CN105938203A (zh) | 一种储层特性的检测方法及装置 | |
CN107728214A (zh) | 一种裂缝预测方法 | |
CN105093314B (zh) | 一种测定微地震震源的方法 | |
Lü et al. | Crustal-scale structure and deformation of Lu-Zong ore district: Joint interpretation from integrated geophysical data | |
CN106443791B (zh) | 求取倾斜地层或各向异性地层横波剩余静校正量的方法 | |
CN104422955B (zh) | 一种利用旅行时变化量进行各向异性参数提取的方法 | |
CN107942372A (zh) | 井间及井地联合地震波ct成像方法及装置 | |
Wang et al. | Polarization migration of three-component reflected waves under small migration aperture condition | |
Mari et al. | Contribution of seismic and acoustic methods to reservoir model building | |
Greenwood et al. | Characterization and imaging of a hydrothermally active near-vertical fault zone in crystalline rocks based on hydrophone VSP 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 | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20190115 Termination date: 20200105 |