CN104133245A - 一种地震资料的静校正方法及系统 - Google Patents

一种地震资料的静校正方法及系统 Download PDF

Info

Publication number
CN104133245A
CN104133245A CN201410351186.5A CN201410351186A CN104133245A CN 104133245 A CN104133245 A CN 104133245A CN 201410351186 A CN201410351186 A CN 201410351186A CN 104133245 A CN104133245 A CN 104133245A
Authority
CN
China
Prior art keywords
shot point
geophone station
velocity
data
experiment
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
Application number
CN201410351186.5A
Other languages
English (en)
Other versions
CN104133245B (zh
Inventor
冯发全
张宇生
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China National Petroleum Corp
BGP Inc
Original Assignee
China National Petroleum Corp
BGP Inc
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by China National Petroleum Corp, BGP Inc filed Critical China National Petroleum Corp
Priority to CN201410351186.5A priority Critical patent/CN104133245B/zh
Publication of CN104133245A publication Critical patent/CN104133245A/zh
Application granted granted Critical
Publication of CN104133245B publication Critical patent/CN104133245B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Geophysics And Detection Of Objects (AREA)

Abstract

本发明公开了一种地震资料的静校正方法及系统,该方法包括:根据实验地区的地震数据获取初至时间,在实验地区的CMP道集内进行折射分层,通过计算得到每一炮点或检波点的折射速度、延迟时;建立实验地区的初始模型,根据初至时间层析反演得到地下速度空间模型;选择某一部分炮点或检波点进行微测井调查,获取某一部分炮点或检波点的高速层埋深数据,并在地下速度空间模型上进行深度标定,获得对应的层析反演的地层速度;对其它炮点或检波点进行空间内插,得到每一炮点或检波点在高速层埋深处对应的层析反演的地层速度,再逐点提取每一炮点或检波点的高速层埋深数据;计算得到每一炮点或检波点的静校正量,通过静校正量对地震资料进行静校正处理。

Description

一种地震资料的静校正方法及系统
技术领域
本发明是关于地球物理勘探领域,涉及一种地震资料的静校正方法及系统。
背景技术
静校正是对地震资料所作的校正,用于补偿由高程、风化层厚度以及风化层速度产生的影响,把资料校到一个指定的基准面上。其目的通常是获得在一个平面上进行采集,且没有风化层或低速介质存在时的反射波到达时间。在地震勘探区域,一般要进行表层调查,获取有关表层的地球物理参数,直接采用控制点内插或者采用大炮初至信息并结合表层调查信息来反演建立表层模型,通过表层模型计算基准面静校正来消除近地表的影响从而提高地震反射波的成像质量,是一种地震勘探中的主要数据处理方法。
在现有的静校正应用中,例如《折射层析反演静校正在伊朗卡山地区的应用》提到:折射层析反演静校正能够较好地解决复杂近地表引起的静校正量问题。该技术最新的算法利用非显式的射线追踪方法求解近地表模型,有更强的适应性和灵活性,减少了复杂地表和地质条件下速度/深度模型的不确定性,提高了模型的精度。在伊朗卡山地区三维地震资料处理过程中,利用该项技术很好地解决了近地表引起的静校正问题。
另外,如《山前冲积扇表层静校正对策》提到:由于山前地表冲积扇的厚度和速度变化较大,利用传统的表层调查和静校正方法计算无法彻底解决“中长波长静校正量”问题。具体表现为小折射排列长度太小(当低速带厚度较大时,无法追踪到高速层顶界面)、小折射资料的解释算法不合适及微测井深度明显不够。通过对雁木西地区多条过井测线进行系统分析(内容包括采集方法、处理技术、速度建场及成图方法等),从采集、处理和解释入手对表层结构进行研究,采用大排列小观测折射及微测井观测,层析成像等技术,并利用初至波形成了一套完善的表层静校正技术,先后在雁木西和胜北地区应用,取得了很好的效果。
上述方法都采用表层调查信息(小折射、微测井)和大炮初至信息反演建立表层模型,能解决许多复杂地区的静校正问题,但在低降速带巨厚、表层结构复杂多变的地区,上述建模方法存在一定的局部不确定性,即对模型的描述和控制精度不够,有可能存在中长波长静校正的风险,也就是无法建立准确的表层模型用于解决地震资料处理中的静校正问题。
发明内容
对于上述地震资料处理中的静校正问题,本发明针对低降速带巨厚并且表层结构复杂多变地区,将微测井或其它方法得到的高速层埋深的深度数据、结合大炮初至采用折射反演方法得到的延迟时和折射速度数据、以及层析反演得到的空间速度分布数据,通过深度标定的方法将这三种方法的优势综合在一起联合建立表层模型,提出了一种地震资料的静校正方法及系统,用于对地震资料进行静校正处理,从而克服了地震勘探中复杂地区表层模型的精度较低,中长波长静校正可能存在的风险,对模型的描述困难且存在较大的误差等问题。
为达到上述目的,本发明提出了一种地震资料的静校正方法,包括以下步骤:步骤1,利用地震数据采集装置获取实验地区的地震数据;步骤2,根据所述地震数据,获取每一炮点的一定偏移距范围内的初至时间,以及获取所述实验地区的CMP道集,在所述CMP道集内进行折射分层,并通过计算得到每一炮点或检波点的折射速度以及所述每一炮点或检波点的延迟时;步骤3,设定所述实验地区的面元网格,并根据所述地震数据建立所述实验地区的初始模型,再根据所述初至时间对所述初始模型进行层析反演,得到所述面元网格一定深度范围内的地下速度空间模型;步骤4,根据所述实验地区地表和表层岩性分布特点,在所述实验地区内选择某一部分炮点或检波点进行微测井调查,获取所述某一部分炮点或检波点的高速层埋深数据,并在所述地下速度空间模型上进行深度标定,获得所述某一部分炮点或检波点在高速层埋深处对应的层析反演的地层速度;步骤5,根据所述某一部分炮点或检波点在所述高速层埋深处对应的层析反演的地层速度,对所述实验地区内除某一部分炮点或检波点外的其它炮点或检波点进行空间内插,得到所述实验地区内每一炮点或检波点在高速层埋深处对应的层析反演的地层速度,再根据所述实验地区内每一炮点或检波点在高速层埋深处对应的层析反演的地层速度逐点提取所述每一炮点或检波点的高速层埋深数据;步骤6,根据所述每一炮点或检波点的高速层埋深数据、所述每一炮点或检波点的折射速度、所述每一炮点或检波点的延迟时计算得到所述每一炮点或检波点的表层平均速度;步骤7,设定水平基准面高程和填充速度,并根据所述每一炮点或检波点的表层平均速度、所述每一炮点或检波点的高速层埋深数据计算得到所述每一炮点或检波点的静校正量,通过所述静校正量对所述地震资料进行静校正处理。
为达到上述目的,本发明还提出了一种地震资料的静校正系统,所述系统包括:地震数据采集模块,用于获取实验地区的地震数据;参数获取模块,用于根据所述地震数据,获取每一炮点的一定偏移距范围内的初至时间,以及获取所述实验地区的CMP道集,在所述CMP道集内进行折射分层,并通过计算得到每一炮点或检波点的折射速度以及所述每一炮点或检波点的延迟时;模型处理模块,用于设定所述实验地区的面元网格,并根据所述地震数据建立所述实验地区的初始模型,再根据所述初至时间对所述初始模型进行层析反演,得到所述面元网格一定深度范围内的地下速度空间模型;微测井及深度标定处理模块,用于根据所述实验地区地表和表层岩性分布特点,在所述实验地区内选择某一部分炮点或检波点进行微测井调查,获取所述某一部分炮点或检波点的高速层埋深数据,并在所述地下速度空间模型上进行深度标定,获得所述某一部分炮点或检波点在高速层埋深处对应的层析反演的地层速度;空间内插处理模块,用于根据所述某一部分炮点或检波点在所述高速层埋深处对应的层析反演的地层速度,对所述实验地区内除某一部分炮点或检波点外的其它炮点或检波点进行空间内插,得到所述实验地区内每一炮点或检波点在高速层埋深处对应的层析反演的地层速度,再根据所述实验地区内每一炮点或检波点在高速层埋深处对应的层析反演的地层速度逐点提取所述每一炮点或检波点的高速层埋深数据;表层平均速度计算模块,用于根据所述每一炮点或检波点的高速层埋深数据、所述每一炮点或检波点的折射速度、所述每一炮点或检波点的延迟时计算得到所述每一炮点或检波点的表层平均速度;静校正处理模块,用于设定水平基准面高程和填充速度,并根据所述每一炮点或检波点的表层平均速度、所述每一炮点或检波点的高速层埋深数据计算得到所述每一炮点或检波点的静校正量,通过所述静校正量对所述地震资料进行静校正处理。
本发明的地震资料的静校正方法及系统用于对地震资料进行静校正处理,进一步提高了对实验地区表层模型的描述能力和表层建模精度,消除表层介质对地震波传播的影响,更好满足复杂地区地震成像对静校正精度的高要求。
附图说明
此处所说明的附图用来提供对本发明的进一步理解,构成本申请的一部分,并不构成对本发明的限定。在附图中:
图1为本发明一实施例的地震资料的静校正方法流程图。
图2为本发明一实施例的地震资料的静校正系统示意图。
图3A是某地区某二维测线应用常规方法建立的表层模型剖面示意图。
图3B是与图3A相同二维测线采用本发明的地震资料的静校正方法建模的表层模型剖面示意图。
图4A是根据图3A模型计算的基准面静校正量曲线示意图。
图4B是根据图3B模型计算的基准面静校正量曲线示意图。
图5A是某地区某二维测线应用常规建模方法的静校正剖面示意图。
图5B是某地区某二维测线采用本发明的地震资料的静校正方法建模的静校正剖面示意图。
具体实施方式
以下配合图式及本发明的较佳实施例,进一步阐述本发明为达成预定发明目的所采取的技术手段。
图1为本发明一实施例的地震资料的静校正方法流程图。如图1所示,该方法包括以下步骤:
步骤S101,利用地震数据采集装置获取实验地区的地震数据。
其中,本实施例的步骤方法主要在进行过野外采集作业后的实验地区实施。
在本步骤中,获取的地震数据可以包括:根据野外地震资料采集并整理生成的SPS文件,如炮点文件、检波点文件、关系文件等等,这些文件是石油地球物理勘探中通用的标准数据文件,对数据内容及格式有严格的规定和要求,再根据这些文件可以建立一数据库,以供后续步骤使用。
步骤S102,根据地震数据,获取每一炮点的一定偏移距范围内的初至时间;另外,根据地震数据,还可获取实验地区的CMP道集,在CMP道集内进行折射分层,并通过计算得到每一炮点或检波点的折射速度以及每一炮点或检波点的延迟时。
在本步骤中,一定偏移距范围内是指根据不同地区地震资料的具体情况,选择某一折射初至的地震道的分布范围,对该范围内的地震道拾取初至时间。初至时间指地震波最先到达某一地震道的时间。
在本步骤中,CMP道集可由软件通过数据库(由前面的SPS文件生成)自动计算得到,然后在CMP道集内进行折射分层后,可以通过采用互换速度分析方法计算得到每一炮点或检波点的折射速度,在计算出折射速度的基础上,再计算延迟时;例如,再通过扩展广义互换算法(EGRM)或高斯-赛德尔方法计算每一炮点或检波点的延迟时。
步骤S103,设定实验地区的面元网格,并根据地震数据建立实验地区的初始模型,再根据初至时间对初始模型进行层析反演,得到面元网格一定深度范围内的地下速度空间模型。
其中,设定实验地区的面元表格包括:
在二维地震勘探时,设定测线网格,以单条测线建立实验地区的初始模型。
在三维地震勘探时,设定测网网格。建立整个实验地区的初始模型。
在本实施例中,初始模型也就是表层模型,是在数据库的基础上建立,网格的大小考虑了道距(三维还考虑接收线距)及表层模型的深度。
在本步骤中,一定深度范围内要根据工区的具体情况来决定,一般位于该工区最低表层界面下面的200米左右。根据经验或者试验确定。
步骤S104,根据实验地区地表和表层岩性分布特点,在实验地区内选择某一部分炮点或检波点进行微测井调查,获取某一部分炮点或检波点的高速层埋深数据,并在地下速度空间模型上进行深度标定,获得某一部分炮点或检波点在高速层埋深处对应的层析反演的地层速度。
由于微测井测量是一种较准确的表层结构调查方法,不受复杂地表条件限制,适用于复杂地区观测,但其生产成本高,不能保证逐点连续观测,因而不能直接获得连续的近地表模型。因此,在本步骤中,根据实验地区地表和表层岩性分布特点,在实验地区内选择某一部分炮点或检波点进行微测井调查。
另外,在本步骤中,采用微测井调查时可以结合其他资料,如大钻取芯资料、声波测井资料等,这些是工区内由第三方完成的,如果有,可以参考其结果。
步骤S105,根据某一部分炮点或检波点在高速层埋深处对应的层析反演的地层速度,对实验地区内除某一部分炮点或检波点外的其它炮点或检波点进行空间内插,得到实验地区内每一炮点或检波点在高速层埋深处对应的层析反演的地层速度,再根据实验地区内每一炮点或检波点在高速层埋深处对应的层析反演的地层速度逐点提取每一炮点或检波点的高速层埋深数据。
在本步骤中,空间内插包括:
当进行二维地震勘探时,对实验地区内除某一部分炮点或检波点外的其它炮点或检波点进行线性内插,也就得到了其他炮点或检波点在高速层埋深处对应的层析反演的地层速度,从而得到了实验地区每一炮点或检波点在高速层埋深处对应的层析反演的地层速度。
当进行三维地震勘探时,是对实验地区内除某一部分炮点或检波点外的其它炮点或检波点进行平面内插。
在本步骤中,最后得到每一炮点或检波点的高速层埋深数据,实际上也就是建立了整个实验地区(二维测线或三维侧网)的等效表层模型,通过该表层模型,结合下面步骤S106,进一步计算出每一炮点或检波点的表层平均速度。
步骤S106,根据每一炮点或检波点的高速层埋深数据、每一炮点或检波点的折射速度、每一炮点或检波点的延迟时计算得到每一炮点或检波点的表层平均速度。
其中,计算每一炮点或检波点的表层平均速度利用的公式为:
V W = ( Z - H G ) V R τ 2 V R 2 + ( Z - H G ) 2 ; - - - ( 1 )
其中,VW为每一炮点或检波点的表层平均速度;
VR为步骤S102获得的每一炮点或检波点的折射速度;
τ为步骤S102获得的每一炮点或检波点的延迟时;
Z为每一炮点或检波点的地表高程,从步骤S101获取的地震数据中直接得到,主要是由炮点文件和检波点文件提供;
HG为步骤S105获得的每一炮点或检波点的高速层埋深数据。
步骤S107,设定水平基准面高程和填充速度,并根据每一炮点或检波点的表层平均速度、每一炮点或检波点的高速层埋深数据计算得到每一炮点或检波点的静校正量,通过静校正量对地震资料进行静校正处理。
其中,计算每一炮点或检波点的静校正量利用的公式为:
TS r = E D - H G V C - Z - H G V W + h s V 0 ; - - - ( 2 )
TR j = E D - H G V C - Z - H G V W ; - - - ( 3 )
其中,TSr为每一炮点的静校正量;
TRj为每一检波点的静校正量;
ED为设定的水平基准面高程;
VC为设定的填充速度;
Z为每一炮点或检波点的地表高程,从步骤S101获取的地震数据中直接得到,主要是由炮点文件和检波点文件提供;
HG为步骤S105获得的每一炮点或检波点的高速层埋深数据;
VW为公式(1)得到的每一炮点或检波点的表层平均速度;
hs为每一炮点的井深,从步骤S101获取的地震数据中直接得到;
V0为每一炮点的表层速度,直接根据实验地区设定。一般某一实验地区可以根据地区情况或结合地区资料直接设定一个地区经验值。
在前述公式(2)中,求出的量为每一炮点的井口时间。
基于同一发明构思,本发明实施例中还提供了一种地震资料的静校正系统,如下面的实施例所述。
由于地震资料的静校正系统解决问题的原理与地震资料的静校正方法相似,因此地震资料的静校正系统的实施可以参见地震资料的静校正方法的实施,重复之处不再赘述。以下所使用的,术语“单元”或者“模块”可以实现预定功能的软件和/或硬件的组合。尽管以下实施例所描述的装置较佳地以软件来实现,但是硬件,或者软件和硬件的组合的实现也是可能并被构想的。
图2为本发明一实施例的地震资料的静校正系统示意图。如图2所示,包括:
地震数据采集模块11,用于获取实验地区的地震数据。
参数获取模块12,用于根据地震数据,获取每一炮点的一定偏移距范围内的初至时间,以及获取实验地区的CMP道集,在CMP道集内进行折射分层,并通过计算得到每一炮点或检波点的折射速度以及每一炮点或检波点的延迟时。
模型处理模块13,用于设定实验地区的面元网格,建立实验地区的初始模型,再根据初至时间对初始模型进行层析反演,得到面元网格一定深度范围内的地下速度空间模型。
微测井及深度标定处理模块14,用于根据实验地区地表和表层岩性分布特点,在实验地区内选择某一部分炮点或检波点进行微测井调查,获取某一部分炮点或检波点的高速层埋深数据,并在地下速度空间模型上进行深度标定,获得某一部分炮点或检波点在高速层埋深处对应的层析反演的地层速度。
空间内插处理模块15,用于根据某一部分炮点或检波点在高速层埋深处对应的层析反演的地层速度,对实验地区内除某一部分炮点或检波点外的其它炮点或检波点进行空间内插,得到实验地区内每一炮点或检波点在高速层埋深处对应的层析反演的地层速度,再根据实验地区内每一炮点或检波点在高速层埋深处对应的层析反演的地层速度逐点提取每一炮点或检波点的高速层埋深数据。
表层平均速度计算模块16,用于根据每一炮点或检波点的高速层埋深数据、每一炮点或检波点的折射速度、每一炮点或检波点的延迟时计算得到每一炮点或检波点的表层平均速度。
静校正处理模块17,用于设定水平基准面高程和填充速度,并根据每一炮点或检波点的表层平均速度、每一炮点或检波点的高速层埋深数据计算得到每一炮点或检波点的静校正量,通过静校正量对地震资料进行静校正处理。
在本实施例中,地震数据采集模块11,用于获取实验地区的地震数据,还包括:炮点文件、检波点文件、关系文件在内的SPS文件数据。
在本实施例中,参数获取模块12还包括:根据地震数据,获取实验地区的CMP道集,并在CMP道集内进行折射分层,采用互换速度分析方法计算得到每一炮点或检波点的折射速度,再通过扩展广义互换算法或高斯-赛德尔方法计算每一炮点或检波点的延迟时。
在本实施例中,模型处理模块13,用于设定实验地区的面元网格,建立实验地区的初始模型,再根据初至时间对初始模型进行层析反演,得到面元网格一定深度范围内的地下速度空间模型,还包括:
当进行二维地震勘探时,设定测线网格,根据地震数据以单条测线建立实验地区的表层模型;
当进行三维地震勘探时,设定测网网格,根据地震数据建立整个实验地区的表层模型。
在本实施例中,空间内插处理模块15,还包括:
当进行二维地震勘探时,根据某一部分炮点或检波点在高速层埋深处对应的层析反演的地层速度,对实验地区内除某一部分炮点或检波点外的其它炮点或检波点进行线性内插,得到实验地区内每一炮点或检波点在高速层埋深处对应的层析反演的地层速度。
当进行三维地震勘探时,根据某一部分炮点或检波点在高速层埋深处对应的层析反演的地层速度,对实验地区内除某一部分炮点或检波点外的其它炮点或检波点进行平面内插,得到实验地区内每一炮点或检波点在高速层埋深处对应的层析反演的地层速度。
在本实施例中,表层平均速度计算模块16计算每一炮点或检波点的表层平均速度利用的公式为:
V W = ( Z - H G ) V R τ 2 V R 2 + ( Z - H G ) 2 ; - - - ( 1 )
其中,VW为每一炮点或检波点的表层平均速度;
VR为参数获取模块12获得的每一炮点或检波点的折射速度;
τ为参数获取模块12获得的每一炮点或检波点的延迟时;
Z为每一炮点或检波点的地表高程,从地震数据采集模块11获取的地震数据中直接得到;
HG为空间内插处理模块15获得的每一炮点或检波点的高速层埋深数据。
在本实施例中,静校正处理模块17计算每一炮点或检波点的静校正量利用的公式为:
TS r = E D - H G V C - Z - H G V W + h s V 0 ; - - - ( 2 )
TR j = E D - H G V C - Z - H G V W ; - - - ( 3 )
其中,TSr为每一炮点的静校正量;
TRj为每一检波点的静校正量;
ED为设定的水平基准面高程;
VC为设定的填充速度;
Z为每一炮点或检波点的地表高程,从地震数据采集模块11获取的地震数据中直接得到;
HG为空间内插处理模块15获得的每一炮点或检波点的高速层埋深数据;
VW为表层平均速度计算模块16得到的每一炮点或检波点的表层平均速度;
hs为每一炮点的井深,从地震数据采集模块11获取的地震数据中直接得到;
V0为每一炮点的表层速度,直接根据实验地区设定。
为了对上述弹性波正演模拟方法及系统进行更为清楚的解释,下面结合一个具体的实施例来进行说明,然而值得注意的是该实施例仅是为了更好地说明本发明,并不构成对本发明不当的限定。
结合步骤S101,利用野外地震资料采集并经整理生成的SPS文件(炮点、检波点和关系文件)建立数据库。
结合步骤S102,根据数据库,拾取全部生产炮记录一定偏移距范围(偏移距3000米)内的初至时间。
再结合步骤S102,在CMP道集内进行折射分层,然后采用互换速度分析方法计算折射速度,采用高斯—赛德尔方法计算延迟时,得到每一个检波点和炮点的折射速度和延迟时。
结合步骤S103,定义测线网格(网格大小为20m×2m),建立初始模型(速度梯度模型,初始速度400m/s,速度梯度6,模型底界为海拔200米),根据初至时间数据进行层析反演,得到该测线的地下速度空间模型。
结合步骤S104及步骤S105,根据实验地区地表和表层岩性分布特点,选择某一部分炮点或检波点进行微测井调查,获取高速层埋深数据,并在地下速度空间模型上进行深度标定,获得某一部分炮点或检波点在高速层埋深处对应的层析反演的地层速度。再依此数据对实验地区内其它炮点或检波点进行空间内插,得到实验地区内所有炮点或检波点在高速层埋深处对应的层析反演的地层速度,最后根据实验地区内每一炮点或检波点在高速层埋深处对应的层析反演的地层速度逐点提取每一炮点或检波点的高速层埋深数据。
结合步骤S106及S107,计算得到每一炮点或检波点的表层平均速度,定义全区统一的水平基准面高程ED(1400米)和填充速度VC(4000m/s),最后计算出实验地区内每一炮点或检波点的基准面静校正量,根据炮点基准面静校正量和检波点基准面静校正量对地震资料的处理。
在一具体实施例中,图3A是某地区某二维测线应用常规方法建立的表层模型剖面示意图。图3B是与图3A相同二维测线采用本发明的地震资料的静校正方法建模的表层模型剖面示意图。
结合图3A和图3B表层模型来看,常规方法建立的表层模型厚度较深,而采用本方法建立的表层模型厚度较小,与实际模型更为接近。
图4A是根据图3A模型计算的基准面静校正量曲线示意图。图4B是根据图3B模型计算的基准面静校正量曲线示意图。
结合4A及图4B可以看出,静校正量整体变化趋势一致,都在数值上有一定差异,长波长有一定变化;表明采用深度标定的综合建模方法建立的表层模型更为准确,基准面静校正精度高,长波长和短波长静校正都解决的好,成像质量高,构造准确。
图5A是某地区某二维测线应用常规建模方法的静校正剖面示意图。图5B是某地区某二维测线采用本发明的地震资料的静校正方法建模的静校正剖面示意图。
再结合图5A及图5B来看地震资料的剖面效果,对比两图后,图5A为在应用常规静校正的剖面上,存在一定的长波长静校正,构造不准确,成像质量局部地区较差;而如图5B所示,采用深度标定的综合建模静校正的静校正剖面上,构造准确,成像质量明显提高,浅、中、深反射都得到加强,信噪比也得到提高,整个剖面的成像质量变好。
本发明的地震资料的静校正方法及系统用于对地震资料进行静校正处理,进一步提高了对实验地区表层模型的描述能力和表层建模精度,消除表层介质对地震波传播的影响,更好满足复杂地区地震成像对静校正精度的高要求。
以上所述的具体实施例,对本发明的目的、技术方案和有益效果进行了进一步详细说明,所应理解的是,以上所述仅为本发明的具体实施例而已,并不用于限定本发明的保护范围,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。

Claims (14)

1.一种地震资料的静校正方法,其特征在于,包括以下步骤:
步骤1,利用地震数据采集装置获取实验地区的地震数据;
步骤2,根据所述地震数据,获取每一炮点的一定偏移距范围内的初至时间,以及获取所述实验地区的CMP道集,在所述CMP道集内进行折射分层,并通过计算得到每一炮点或检波点的折射速度以及所述每一炮点或检波点的延迟时;
步骤3,设定所述实验地区的面元网格,并根据所述地震数据建立所述实验地区的初始模型,再根据所述初至时间对所述初始模型进行层析反演,得到所述面元网格一定深度范围内的地下速度空间模型;
步骤4,根据所述实验地区地表和表层岩性分布特点,在所述实验地区内选择某一部分炮点或检波点进行微测井调查,获取所述某一部分炮点或检波点的高速层埋深数据,并在所述地下速度空间模型上进行深度标定,获得所述某一部分炮点或检波点在高速层埋深处对应的层析反演的地层速度;
步骤5,根据所述某一部分炮点或检波点在所述高速层埋深处对应的层析反演的地层速度,对所述实验地区内除某一部分炮点或检波点外的其它炮点或检波点进行空间内插,得到所述实验地区内每一炮点或检波点在高速层埋深处对应的层析反演的地层速度,再根据所述实验地区内每一炮点或检波点在高速层埋深处对应的层析反演的地层速度逐点提取所述每一炮点或检波点的高速层埋深数据;
步骤6,根据所述每一炮点或检波点的高速层埋深数据、所述每一炮点或检波点的折射速度、所述每一炮点或检波点的延迟时计算得到所述每一炮点或检波点的表层平均速度;
步骤7,设定水平基准面高程和填充速度,并根据所述每一炮点或检波点的表层平均速度、所述每一炮点或检波点的高速层埋深数据计算得到所述每一炮点或检波点的静校正量,通过所述静校正量对所述地震资料进行静校正处理。
2.根据权利要求1所述的地震资料的静校正方法,其特征在于,所述步骤1,利用地震数据采集装置获取实验地区的地震数据还包括:炮点文件、检波点文件、关系文件在内的SPS文件数据。
3.根据权利要求1所述的地震资料的静校正方法,其特征在于,所述步骤2中,根据所述地震数据,获取所述实验地区的CMP道集,在所述CMP道集内进行折射分层,并通过计算得到每一炮点或检波点的折射速度以及所述每一炮点或检波点的延迟时包括:
根据所述地震数据,获取所述实验地区的CMP道集,并在所述CMP道集内进行折射分层,采用互换速度分析方法计算得到每一炮点或检波点的折射速度,再通过扩展广义互换算法或高斯-赛德尔方法计算所述每一炮点或检波点的延迟时。
4.根据权利要求1所述的地震资料的静校正方法,其特征在于,所述步骤3中,设定所述实验地区的面元网格,并根据所述地震数据建立所述实验地区的初始模型还包括:
当进行二维地震勘探时,设定测线网格,根据所述地震数据以单条测线建立所述实验地区的初始模型;
当进行三维地震勘探时,设定测网网格,根据所述地震数据建立整个所述实验地区的初始模型。
5.根据权利要求4所述的地震资料的静校正方法,其特征在于,所述步骤5中,所述空间内插还包括:
当进行二维地震勘探时,根据所述某一部分炮点或检波点在所述高速层埋深处对应的层析反演的地层速度,对所述实验地区内除某一部分炮点或检波点外的其它炮点或检波点进行线性内插,得到所述实验地区内每一炮点或检波点在高速层埋深处对应的层析反演的地层速度;
当进行三维地震勘探时,根据所述某一部分炮点或检波点在所述高速层埋深处对应的层析反演的地层速度,对所述实验地区内除某一部分炮点或检波点外的其它炮点或检波点进行平面内插,得到所述实验地区内每一炮点或检波点在高速层埋深处对应的层析反演的地层速度。
6.根据权利要求5所述的地震资料的静校正方法,其特征在于,所述步骤6中,根据所述每一炮点或检波点的高速层埋深数据、所述每一炮点或检波点的折射速度、所述每一炮点或检波点的延迟时计算得到所述每一炮点或检波点的表层平均速度利用的公式为:
V W = ( Z - H G ) V R τ 2 V R 2 + ( Z - H G ) 2 ; - - - ( 1 )
其中,VW为所述每一炮点或检波点的表层平均速度;
VR为所述步骤2获得的所述每一炮点或检波点的折射速度;
τ为所述步骤2获得的所述每一炮点或检波点的延迟时;
Z为所述每一炮点或检波点的地表高程,从所述步骤1获取的地震数据中直接得到;
HG为所述步骤5获得的所述每一炮点或检波点的高速层埋深数据。
7.根据权利要求6所述的地震资料的静校正方法,其特征在于,所述步骤8中,设定水平基准面高程和填充速度,并根据所述每一炮点或检波点的表层平均速度、所述每一炮点或检波点的高速层埋深数据计算得到所述每一炮点或检波点的静校正量利用的公式为:
TS r = E D - H G V C - Z - H G V W + h s V 0 ; - - - ( 2 )
TR j = E D - H G V C - Z - H G V W ; - - - ( 3 )
其中,TSr为所述每一炮点的静校正量;
TRj为每一所述检波点的静校正量;
ED为设定的所述水平基准面高程;
VC为设定的填充速度;
Z为所述每一炮点或检波点的地表高程,从所述步骤1获取的地震数据中直接得到;
HG为所述步骤5获得的所述每一炮点或检波点的高速层埋深数据;
VW为公式(1)得到的所述每一炮点或检波点的表层平均速度;
hs为所述每一炮点的井深,从所述步骤1获取的地震数据中直接得到;
V0为所述每一炮点的表层速度,直接根据所述实验地区设定。
8.一种地震资料的静校正系统,其特征在于,所述系统包括:
地震数据采集模块,用于获取实验地区的地震数据;
参数获取模块,用于根据所述地震数据,获取每一炮点的一定偏移距范围内的初至时间,以及获取所述实验地区的CMP道集,在所述CMP道集内进行折射分层,并通过计算得到每一炮点或检波点的折射速度以及所述每一炮点或检波点的延迟时;
模型处理模块,用于设定所述实验地区的面元网格,并根据所述地震数据建立所述实验地区的初始模型,再根据所述初至时间对所述初始模型进行层析反演,得到所述面元网格一定深度范围内的地下速度空间模型;
微测井及深度标定处理模块,用于根据所述实验地区地表和表层岩性分布特点,在所述实验地区内选择某一部分炮点或检波点进行微测井调查,获取所述某一部分炮点或检波点的高速层埋深数据,并在所述地下速度空间模型上进行深度标定,获得所述某一部分炮点或检波点在高速层埋深处对应的层析反演的地层速度;
空间内插处理模块,用于根据所述某一部分炮点或检波点在所述高速层埋深处对应的层析反演的地层速度,对所述实验地区内除某一部分炮点或检波点外的其它炮点或检波点进行空间内插,得到所述实验地区内每一炮点或检波点在高速层埋深处对应的层析反演的地层速度,再根据所述实验地区内每一炮点或检波点在高速层埋深处对应的层析反演的地层速度逐点提取所述每一炮点或检波点的高速层埋深数据;
表层平均速度计算模块,用于根据所述每一炮点或检波点的高速层埋深数据、所述每一炮点或检波点的折射速度、所述每一炮点或检波点的延迟时计算得到所述每一炮点或检波点的表层平均速度;
静校正处理模块,用于设定水平基准面高程和填充速度,并根据所述每一炮点或检波点的表层平均速度、所述每一炮点或检波点的高速层埋深数据计算得到所述每一炮点或检波点的静校正量,通过所述静校正量对所述地震资料进行静校正处理。
9.根据权利要求8所述的地震资料的静校正系统,其特征在于,所述地震数据采集模块,用于获取实验地区的地震数据,还包括:炮点文件、检波点文件、关系文件在内的SPS文件数据。
10.根据权利要求8所述的地震资料的静校正系统,其特征在于,所述参数获取模块,用于根据所述地震数据,获取每一炮点的一定偏移距范围内的初至时间,以及获取所述实验地区的CMP道集,在所述CMP道集内进行折射分层,并通过计算得到每一炮点或检波点的折射速度以及所述每一炮点或检波点的延迟时,还包括:
根据所述地震数据,获取所述实验地区的CMP道集,并在所述CMP道集内进行折射分层,采用互换速度分析方法计算得到每一炮点或检波点的折射速度,再通过扩展广义互换算法或高斯-赛德尔方法计算所述每一炮点或检波点的延迟时。
11.根据权利要求8所述的地震资料的静校正系统,其特征在于,所述模型处理模块,用于设定所述实验地区的面元网格,并根据所述地震数据建立所述实验地区的初始模型,再根据所述初至时间对所述初始模型进行层析反演,得到所述面元网格一定深度范围内的地下速度空间模型,还包括:
当进行二维地震勘探时,设定测线网格,根据所述地震数据以单条测线建立所述实验地区的初始模型;
当进行三维地震勘探时,设定测网网格,根据所述地震数据建立整个所述实验地区的初始模型。
12.根据权利要求11所述的地震资料的静校正系统,其特征在于,所述空间内插处理模块,用于根据所述某一部分炮点或检波点在所述高速层埋深处对应的层析反演的地层速度,对所述实验地区内除某一部分炮点或检波点外的其它炮点或检波点进行空间内插,得到所述实验地区内每一炮点或检波点在高速层埋深处对应的层析反演的地层速度,再根据所述实验地区内每一炮点或检波点在高速层埋深处对应的层析反演的地层速度逐点提取所述每一炮点或检波点的高速层埋深数据,还包括:
当进行二维地震勘探时,根据所述某一部分炮点或检波点在所述高速层埋深处对应的层析反演的地层速度,对所述实验地区内除某一部分炮点或检波点外的其它炮点或检波点进行线性内插,得到所述实验地区内每一炮点或检波点在高速层埋深处对应的层析反演的地层速度;
当进行三维地震勘探时,根据所述某一部分炮点或检波点在所述高速层埋深处对应的层析反演的地层速度,对所述实验地区内除某一部分炮点或检波点外的其它炮点或检波点进行平面内插,得到所述实验地区内每一炮点或检波点在高速层埋深处对应的层析反演的地层速度。
13.根据权利要求12所述的地震资料的静校正系统,其特征在于,所述表层平均速度计算模块,用于根据所述每一炮点或检波点的高速层埋深数据、所述每一炮点或检波点的折射速度、所述每一炮点或检波点的延迟时计算得到所述每一炮点或检波点的表层平均速度利用的公式为:
V W = ( Z - H G ) V R τ 2 V R 2 + ( Z - H G ) 2 ; - - - ( 1 )
其中,VW为所述每一炮点或检波点的表层平均速度;
VR为所述参数获取模块获得的所述每一炮点或检波点的折射速度;
τ为所述参数获取模块获得的所述每一炮点或检波点的延迟时;
Z为所述每一炮点或检波点的地表高程,从所述地震数据采集模块获取的地震数据中直接得到;
HG为所述空间内插处理模块获得的所述每一炮点或检波点的高速层埋深数据。
14.根据权利要求13所述的地震资料的静校正系统,其特征在于,所述静校正处理模块,用于设定水平基准面高程和填充速度,并根据所述每一炮点或检波点的表层平均速度、所述每一炮点或检波点的高速层埋深数据计算得到所述每一炮点或检波点的静校正量利用的公式为:
TS r = E D - H G V C - Z - H G V W + h s V 0 ; - - - ( 2 )
TR j = E D - H G V C - Z - H G V W ; - - - ( 3 )
其中,TSr为所述每一炮点的静校正量;
TRj为每一所述检波点的静校正量;
ED为设定的所述水平基准面高程;
VC为设定的填充速度;
Z为所述每一炮点或检波点的地表高程,从所述地震数据采集模块获取的地震数据中直接得到;
HG为所述空间内插处理模块获得的所述每一炮点或检波点的高速层埋深数据;
VW为所述表层平均速度计算模块得到的所述每一炮点或检波点的表层平均速度;
hs为所述每一炮点的井深,从所述地震数据采集模块获取的地震数据中直接得到;
V0为所述每一炮点的表层速度,直接根据所述实验地区设定。
CN201410351186.5A 2014-07-22 2014-07-22 一种地震资料的静校正方法及系统 Active CN104133245B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410351186.5A CN104133245B (zh) 2014-07-22 2014-07-22 一种地震资料的静校正方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410351186.5A CN104133245B (zh) 2014-07-22 2014-07-22 一种地震资料的静校正方法及系统

Publications (2)

Publication Number Publication Date
CN104133245A true CN104133245A (zh) 2014-11-05
CN104133245B CN104133245B (zh) 2018-01-02

Family

ID=51805988

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410351186.5A Active CN104133245B (zh) 2014-07-22 2014-07-22 一种地震资料的静校正方法及系统

Country Status (1)

Country Link
CN (1) CN104133245B (zh)

Cited By (18)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104570091A (zh) * 2014-12-15 2015-04-29 中国石油天然气集团公司 一种获得初至波射线的方法和装置
CN105866839A (zh) * 2016-06-28 2016-08-17 中国矿业大学(北京) 基于共姿态道集的静校正方法及装置
CN106338770A (zh) * 2015-07-10 2017-01-18 中国石油化工股份有限公司 一种炮检点数据互查方法和系统
CN107656314A (zh) * 2017-10-20 2018-02-02 中国石油天然气集团公司 道集记录处理方法、装置及计算机存储介质
CN108051858A (zh) * 2017-12-13 2018-05-18 中国矿业大学(北京) 优化的地震勘探方法及装置
CN109143365A (zh) * 2017-06-27 2019-01-04 中国石油化工股份有限公司 山前带地震静校正方法及计算机可读存储介质
CN109387871A (zh) * 2017-08-09 2019-02-26 中国石油化工股份有限公司 P-sv转换波组合静校正方法及系统
CN109884707A (zh) * 2019-03-20 2019-06-14 中国石油化工股份有限公司 近地表分层时深曲线静校正方法
CN111060966A (zh) * 2019-12-31 2020-04-24 恒泰艾普(北京)能源科技研究院有限公司 一种用于静校正技术的初至拾取方法
RU2722861C1 (ru) * 2019-07-09 2020-06-04 Общество с ограниченной ответственностью "СИГМА-ГЕО" Способ расчета статических поправок
CN111399037A (zh) * 2019-01-02 2020-07-10 中国石油天然气集团有限公司 高速顶界面提取的方法和装置
CN111965702A (zh) * 2020-07-21 2020-11-20 中国石油天然气集团有限公司 一种基于地震初至的剩余静校正量处理方法及系统
CN112162322A (zh) * 2020-09-28 2021-01-01 中海油田服务股份有限公司 一种地震勘探采集系统的同步性质控方法和装置
CN112241022A (zh) * 2019-07-16 2021-01-19 中国石油天然气集团有限公司 基于射线密度生成层析反演模型速度界面的方法和装置
CN112379434A (zh) * 2020-10-30 2021-02-19 中国石油天然气集团有限公司 适用于沙漠区的层析反演静校正的方法及装置
CN112444879A (zh) * 2019-08-30 2021-03-05 中国石油化工股份有限公司 地震勘探数据的处理方法及存储介质
CN112859172A (zh) * 2021-03-05 2021-05-28 中国石油天然气集团有限公司 纵横波微测井资料处理方法及装置
CN113093279A (zh) * 2020-01-08 2021-07-09 中国石油天然气集团有限公司 转换波静校正方法及装置

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6757616B1 (en) * 2003-02-04 2004-06-29 Charles Wayne Emmons Model-free refraction statics solution
US20110010098A1 (en) * 2008-02-25 2011-01-13 China National Petroleum Corporation Method of pre-stack two-dimension-like transformation of three-dimensional seismic record
CN102590864A (zh) * 2011-12-31 2012-07-18 中国石油集团西北地质研究所 两步法层析反演近地表建模方法
CN102819042A (zh) * 2011-06-10 2012-12-12 中国石油天然气集团公司 一种抽稀地震数据的近地表层析反演静校正方法
CN102944896A (zh) * 2012-11-09 2013-02-27 中国石油集团川庆钻探工程有限公司地球物理勘探公司 表层调查数据的模型法静校正方法
CN103869368A (zh) * 2012-12-10 2014-06-18 中国石油天然气集团公司 一种无表层调查资料约束的大炮初至综合建模静校正方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6757616B1 (en) * 2003-02-04 2004-06-29 Charles Wayne Emmons Model-free refraction statics solution
US20110010098A1 (en) * 2008-02-25 2011-01-13 China National Petroleum Corporation Method of pre-stack two-dimension-like transformation of three-dimensional seismic record
CN102819042A (zh) * 2011-06-10 2012-12-12 中国石油天然气集团公司 一种抽稀地震数据的近地表层析反演静校正方法
CN102590864A (zh) * 2011-12-31 2012-07-18 中国石油集团西北地质研究所 两步法层析反演近地表建模方法
CN102944896A (zh) * 2012-11-09 2013-02-27 中国石油集团川庆钻探工程有限公司地球物理勘探公司 表层调查数据的模型法静校正方法
CN103869368A (zh) * 2012-12-10 2014-06-18 中国石油天然气集团公司 一种无表层调查资料约束的大炮初至综合建模静校正方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
徐凯军 等: "复杂山地近地表结构综合调查及建模技术研究", 《物探化探计算技术》 *
李培明 等: "模型约束的三维初至折射静校正", 《石油地球物理勘探》 *

Cited By (29)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104570091A (zh) * 2014-12-15 2015-04-29 中国石油天然气集团公司 一种获得初至波射线的方法和装置
CN104570091B (zh) * 2014-12-15 2018-02-02 中国石油天然气集团公司 一种获得初至波射线的方法和装置
CN106338770A (zh) * 2015-07-10 2017-01-18 中国石油化工股份有限公司 一种炮检点数据互查方法和系统
CN106338770B (zh) * 2015-07-10 2018-09-18 中国石油化工股份有限公司 一种炮检点数据互查方法和系统
CN105866839A (zh) * 2016-06-28 2016-08-17 中国矿业大学(北京) 基于共姿态道集的静校正方法及装置
CN109143365A (zh) * 2017-06-27 2019-01-04 中国石油化工股份有限公司 山前带地震静校正方法及计算机可读存储介质
CN109143365B (zh) * 2017-06-27 2020-08-25 中国石油化工股份有限公司 山前带地震静校正方法及计算机可读存储介质
CN109387871A (zh) * 2017-08-09 2019-02-26 中国石油化工股份有限公司 P-sv转换波组合静校正方法及系统
CN107656314A (zh) * 2017-10-20 2018-02-02 中国石油天然气集团公司 道集记录处理方法、装置及计算机存储介质
CN107656314B (zh) * 2017-10-20 2020-03-10 中国石油天然气集团公司 道集记录处理方法、装置及计算机存储介质
CN108051858A (zh) * 2017-12-13 2018-05-18 中国矿业大学(北京) 优化的地震勘探方法及装置
CN108051858B (zh) * 2017-12-13 2018-11-13 中国矿业大学(北京) 优化的地震勘探方法及装置
CN111399037A (zh) * 2019-01-02 2020-07-10 中国石油天然气集团有限公司 高速顶界面提取的方法和装置
CN111399037B (zh) * 2019-01-02 2022-12-02 中国石油天然气集团有限公司 高速顶界面提取的方法和装置
CN109884707A (zh) * 2019-03-20 2019-06-14 中国石油化工股份有限公司 近地表分层时深曲线静校正方法
CN109884707B (zh) * 2019-03-20 2021-09-07 中国石油化工股份有限公司 近地表分层时深曲线静校正方法
RU2722861C1 (ru) * 2019-07-09 2020-06-04 Общество с ограниченной ответственностью "СИГМА-ГЕО" Способ расчета статических поправок
CN112241022A (zh) * 2019-07-16 2021-01-19 中国石油天然气集团有限公司 基于射线密度生成层析反演模型速度界面的方法和装置
CN112444879A (zh) * 2019-08-30 2021-03-05 中国石油化工股份有限公司 地震勘探数据的处理方法及存储介质
CN112444879B (zh) * 2019-08-30 2022-10-14 中国石油化工股份有限公司 地震勘探数据的处理方法及存储介质
CN111060966A (zh) * 2019-12-31 2020-04-24 恒泰艾普(北京)能源科技研究院有限公司 一种用于静校正技术的初至拾取方法
CN113093279B (zh) * 2020-01-08 2024-05-28 中国石油天然气集团有限公司 转换波静校正方法及装置
CN113093279A (zh) * 2020-01-08 2021-07-09 中国石油天然气集团有限公司 转换波静校正方法及装置
CN111965702A (zh) * 2020-07-21 2020-11-20 中国石油天然气集团有限公司 一种基于地震初至的剩余静校正量处理方法及系统
CN111965702B (zh) * 2020-07-21 2023-08-22 中国石油天然气集团有限公司 一种基于地震初至的剩余静校正量处理方法及系统
CN112162322A (zh) * 2020-09-28 2021-01-01 中海油田服务股份有限公司 一种地震勘探采集系统的同步性质控方法和装置
CN112379434B (zh) * 2020-10-30 2024-04-30 中国石油天然气集团有限公司 适用于沙漠区的层析反演静校正的方法及装置
CN112379434A (zh) * 2020-10-30 2021-02-19 中国石油天然气集团有限公司 适用于沙漠区的层析反演静校正的方法及装置
CN112859172A (zh) * 2021-03-05 2021-05-28 中国石油天然气集团有限公司 纵横波微测井资料处理方法及装置

Also Published As

Publication number Publication date
CN104133245B (zh) 2018-01-02

Similar Documents

Publication Publication Date Title
CN104133245A (zh) 一种地震资料的静校正方法及系统
CN101980054B (zh) 一种在高密度地震静校正处理中建立近地表速度模型的方法
CN102967882B (zh) 地层的层速度模型的建模方法
CN103869368B (zh) 一种无表层调查资料约束的大炮初至综合建模静校正方法
CN108064348A (zh) 一种基于两点射线追踪的地震走时层析反演方法
US10795053B2 (en) Systems and methods of multi-scale meshing for geologic time modeling
Arnulf et al. Structure, seismicity, and accretionary processes at the hot spot‐influenced Axial Seamount on the Juan de Fuca Ridge
CN105093319B (zh) 基于三维地震数据的地面微地震静校正方法
CN102879820B (zh) 基于三角网格的三维表层模型构建方法
CN104459782B (zh) 使用薄层的沿层层析反演地震速度建模方法以及建模单元
CN109884710B (zh) 针对激发井深设计的微测井层析成像方法
CN102879819B (zh) 保持地震波场运动学特征的地震数据处理方法及装置
CN103412333B (zh) 一种静校正基准面确定方法
CN104749617A (zh) 一种多尺度裂缝储层正演模型建立方法
CN105445789A (zh) 基于多次反射折射波约束的三维菲涅尔体旅行时层析成像方法
CN104375172B (zh) 火山岩下伏地层的构造形态正确成像的方法
CN104570068A (zh) 一种基于三维gis的观测系统crp属性分析方法
CN107703540B (zh) 一种微地震定位及层析成像方法
CN106291725A (zh) 一种快速反演地下地质体空间位置的方法
CN102877828A (zh) 一种三维多井联合井地ct成像方法
CN104166159A (zh) 四维微地震监测的裂缝形态处理方法和系统
CN109100792A (zh) 基于台站与三维地震联合采集资料的速度反演方法
CN104237946B (zh) 基于井控的单层反射纵波和反射转换横波的振幅匹配方法
CN104570070B (zh) 一种建立二维近地表地质模型的方法及设备
CN107356972A (zh) 一种各向异性介质的成像方法

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant