CN104656135B - 一种地震剪切波的准确判别方法和系统 - Google Patents

一种地震剪切波的准确判别方法和系统 Download PDF

Info

Publication number
CN104656135B
CN104656135B CN201510107000.6A CN201510107000A CN104656135B CN 104656135 B CN104656135 B CN 104656135B CN 201510107000 A CN201510107000 A CN 201510107000A CN 104656135 B CN104656135 B CN 104656135B
Authority
CN
China
Prior art keywords
data
wave
shearing wave
component
earthquake
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
Application number
CN201510107000.6A
Other languages
English (en)
Other versions
CN104656135A (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.)
PowerChina Beijing Engineering Corp Ltd
Original Assignee
PowerChina Beijing Engineering Corp Ltd
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 PowerChina Beijing Engineering Corp Ltd filed Critical PowerChina Beijing Engineering Corp Ltd
Priority to CN201510107000.6A priority Critical patent/CN104656135B/zh
Publication of CN104656135A publication Critical patent/CN104656135A/zh
Application granted granted Critical
Publication of CN104656135B publication Critical patent/CN104656135B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Geophysics And Detection Of Objects (AREA)

Abstract

本发明公开了一种地震剪切波的准确判别方法和系统,该方法包括步骤:S1,数据格式转换,具体为将不同型号仪器记录的不同格式的波形数据转换为预定的专有数据格式;S2,数据预处理,具体为对S1中格式转换后得到的数据进行分量数据重排和编辑文件头;S3,建立参数文件,具体为建立分析计算所需的解译参数文件;S4,计算分析,具体为根据剪切波的特征对比分析并正确拾取地震剪切波初至走时,自动计算各分量的测试结果;S5,成果取值和输出,具体为对S4得出的测试结果进行取值并输出成果。本发明可根据情况选择任一分量、不同分量的平均值、分段选取不同分量进行组合解译最终成果,保证了成果的可靠性和合理性。

Description

一种地震剪切波的准确判别方法和系统
技术领域
本发明属于工程物探领域,涉及一种工程物探钻孔测试数据分析方法,尤其涉及一种地震剪切波的准确判别方法。
背景技术
钻孔物探测试是所有勘察中常用的手段之一,根据工程特点可选用不同测试方法,在需要对地基进行评价时(特别是需要对有液化土判别的地基),大多需要进行地震剪切波测试。由于地震波列中包含有压缩波、剪切波、面波等多种波型,如何正确判别和提取剪切波不仅是测试成果可靠性的关键所在,而且也考验着工程人员的经验和水平。
钻孔地震剪切波测试中,通常采用三分量检波器接收、正反向激发,因此每个测点的波形记录信息量大。长期以来,工程物探技术人员都是凭借自己丰富的经验和理论基础来分析判别地震剪切波,不仅工作效率低,还易发生误判,直接影响了测试结果的准确性和可靠性。
发明内容
本发明的目的在于提供一种地震剪切波的准确判别方法,该方法能适时切换地震剪切波三分量波形,并将单一分量正反向激发的波形叠加分析,提高了钻孔地震剪切测试数据分析的工作效率和剪切波判别的准确性,从而解决现有技术中存在的前述问题。主要应用于水电工程、水利工程、地震安全评价、工业与民用建筑地基勘察等领域。
为了实现上述目的,本发明采用的技术方案如下:
一种地震剪切波的准确判别方法,包括以下步骤:
S1,数据格式转换,具体为将不同型号仪器记录的不同格式的波形数据转换为预定的专有数据格式;
S2,数据预处理,具体为对S1中格式转换后得到的数据进行分量数据重排和编辑文件头;
S3,建立参数文件,具体为建立分析计算所需的解译参数文件;
S4,计算分析,具体为根据剪切波的特征对比分析并正确拾取地震剪切波初至走时,自动计算各分量的测试结果;
S5,成果取值和输出,具体为对S4得出的测试结果进行取值并输出成果。
优选的,步骤S1中被转换的源数据包括SEG-2格式数据和/或CSP格式数据,转换后得到的目标数据是RPS记录格式数据。
优选的,步骤S1中产生波形数据的仪器包括:R24、S12、NZ24、Geopen、Geode、WZG系列仪器和/或ZDQ系列仪器和/或SWS系列仪器。
优选的,步骤S2包括以下步骤:
S21,数据编辑,即按钻孔深度将各个测点的地震剪切波数据以正向和反向分别进行多文件抽道、经道重排为XYZ三个分量的数据;
S22,编辑文件头,即编辑重排后数据文件头参数。
优选的,所述文件头参数包括:孔号、激发点偏距、接收点深度、系统延时。
优选的,步骤S3具体为:根据步骤S2所得到的各分量数据,选择数据建立解译参数文件。
优选的,步骤S4包括以下步骤:
S41,拾取地震剪切波初至走时,即根据剪切波的特征,在波列中对比分析、识别剪切波,并拾取地震剪切波的初至走时;
S42,计算剪切波速,即根据测点间距和初至走时自动计算各分量的测试结果;
将拾取地震剪切波的初至走时与计算结果动态关联,快速检查和修改地震剪切波初至走时和计算结果,直至得到合理的成果。
优选的,步骤S5具体为:在XYZ三个分量测试结果中,通过分析合理综合取值,确定最终成果并输出。
一种地震剪切波的准确判别系统,包括:
数据格式转换模块,用于将不同型号仪器记录的不同格式的波形数据转换为预定的专有数据格式;
数据预处理模块,用于对格式转换后得到的数据进行分量数据重排和编辑文件头;
参数文件建立模块,用于建立分析计算所需的解译参数文件;
计算分析模块,用于根据剪切波的特征对比分析并正确拾取地震剪切波初至走时,自动计算各分量的测试结果;
成果取值和输出模块,用于对得出的测试结果进行取值并输出成果。
优选的,被转换的源数据包括SEG-2格式数据和/或CSP格式数据,转换后得到的目标数据是RPS记录格式数据。
本发明的有益效果是:
本发明对于地震波的显示是将各测点的记录分别按正反向激发转换成X、Y、Z三个分量的波列,并可切换对比显示各个分量波形;单一分量叠加显示正反两个方向激发的波形,便于地震剪切波的准确识别和分析,灵活选择最优波形作为解译的基础数据;可根据情况选择任一分量、不同分量的平均值、分段选取不同分量进行组合解译最终成果,保证了成果的可靠性和合理性。
附图说明
图1是本发明实施例的步骤流程示意图;
图2是本发明实施例中的建立解译参数界面示意图;
图3是本发明实施例中按钻孔深度合并后X分量地震剪切波形图;
图4是本发明实施例综合取值后的成果示意图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施方式仅仅用以解释本发明,并不用于限定本发明。
如附图1所示,本发明公开了一种地震剪切波的准确判别方法,该方法包括以下步骤:
1、数据格式转换,用于将不同型号仪器记录的不同格式的波形数据转换为本系统专有数据格式;
2、数据预处理,用于对格式转换后的数据进行分量数据重排和编辑文件头;
3、建立参数文件,用于建立分析计算所需的解译参数文件;
4、计算分析,用于根据剪切波的特征对比分析并正确拾取地震剪切波初至走时,自动计算各分量的测试结果;
5、成果选择和输出,用于合理取值(测试结果)、输出成果。
按上述方案,数据格式转换A的源数据包括SEG-2(R24、S12、NZ24、Geopen、Geode、WZG系列和ZDQ系列仪器记录的数据)、CSP(SWS系列仪器)等数据格式,目标数据是本系统特定的RPS记录格式数据。源数据兼容了国内外多种仪器的数据格式(公开),广泛适用于各个行业的钻孔地震剪切波数据分析。
数据预处理,包括:数据编辑,即按钻孔深度将各个测点的地震剪切波数据以正向和反方向(激发方向)分别进行多文件抽道、经道重排为XYZ三个分量的数据;编辑文件头,即编辑重排后数据文件头参数(如:孔号、激发点偏距、接收点深度、系统延时等)。
建立参数文件,根据数据预处理B所得到的各分量数据,合理选择数据建立解译参数文件。
计算分析,包括:拾取地震剪切波初至走时,即根据剪切波的特征(正方向和反方向激发的剪切波反相),在波列中对比分析、识别剪切波,并拾取地震剪切波的初至走时;计算剪切波速,即根据测点间距和初至走时自动计算各分量的测试结果。拾取地震剪切波的初至走时与计算结果动态关联,快速检查和修改地震剪切波初至走时和计算结果,直至得到合理的成果。
成果选择及输出,是在XYZ三个分量测试结果中,通过分析合理综合取值,确定最终成果,并输出。
数据预处理和计算分析,对于地震波的显示是将各测点的记录分别按正反向激发转换成X、Y、Z三个分量的波列,并可切换对比显示各个分量波形;单一分量叠加显示正反两个方向激发的波形,便于地震剪切波的准确识别和分析,灵活选择最优波形作为解译的基础数据;可根据情况选择任一分量、不同分量的平均值、分段选取不同分量进行组合解译最终成果,保证了成果的可靠性和合理性。
以下结合具体实施例对本发明的技术方案进行详细说明:
本实施例是某工程场地一钻孔24m深度(覆盖层)内地震剪切波速测试,源数据是美国产R24地震仪记录的数据(SEG-2格式)。
本实施例包括:数据格式转换A、数据预处理B、建立参数文件C、计算分析D、成果选择和输出E五个过程步骤。
1、数据格式转换A:
本实施例是将美国产R24地震仪记录的数据(SEG-2格式)转换为本系统特定的RPS数据格式。
本过程是通过交互界面选择源数据(自动识别仪器型号),进行数据格式转换,保存为本系统特定的RPS数据格式(目标数据)。源数据有SEG-1、SEG-2、SEG-D、CSP等多种格式可选。
2、数据预处理B:
本实施例的主要步骤是将由各测点记录的XYZ三分量地震波数据经多个数据文件抽道后按钻孔深度重新排列成X、Y、Z三个分量的数据文件、并对文件头(如:孔号、激发点偏距、接收点深度、系统延时等)进行编辑、适当调节各分量正反向激发的地震波的震幅。
本过程包含有数据编辑(单个或多个数据文件抽道、数据合并、数据拆分、数据重采样、数据剪切、数据冲零、道序重排)、文件头编辑等多个功能模块,针对实际中不同的数据,可灵活选择相应的处理方式。
3、建立参数文件C:
本实施例是采用正反方向激发、XYZ三个分量地震波数据建立参数文件。
本过程是根据预处理后地震波数据的质量、激发方式等因素,建立解译参数文件。
4、计算分析D:
本实施例是在同时显示正反向激发的X(或Y或Z)分量地震波形窗口,根据正反向激发的剪切波相互反相的特征对比分析,准确判别地震剪切波,并拾取初至走时,同步计算出各测试段的剪切波速。
本过程通过切换显示出XYZ三个分量地震波形(正反方向激发叠加显示),拾取各分量地震剪切波初至走时,并同步计算出各分量的测试结果(地震剪切波速)。拾取剪切波初至走时和测试结果动态关联。
5、成果选择和输出E:
本实施例是选取各测点XYZ三个分量的平均值作为最终成果。
XYZ三个分量的测试结果因为受到地震波数据质量的影响,在不同的测试深度上的可靠度会有所不同,本过程就是通过综合分析,合理选择某一分量或选择不同分量可靠测试段组合、或选择各分量的平均值作为最终成果,并输出。
通过采用本发明公开的上述技术方案,得到了如下有益的效果:
本发明对于地震波的显示是将各测点的记录分别按正反向激发转换成X、Y、Z三个分量的波列,并可切换对比显示各个分量波形;单一分量叠加显示正反两个方向激发的波形,便于地震剪切波的准确识别和分析,灵活选择最优波形作为解译的基础数据;可根据情况选择任一分量、不同分量的平均值、分段选取不同分量进行组合解译最终成果,保证了成果的可靠性和合理性。
以上所述仅是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以做出若干改进和润饰,这些改进和润饰也应视本发明的保护范围。

Claims (2)

1.一种地震剪切波的准确判别方法,其特征在于,包括以下步骤:
S1,数据格式转换,具体为将不同型号仪器记录的不同格式的波形数据转换为预定的专有数据格式;步骤S1中被转换的源数据包括SEG-2格式数据和/或CSP格式数据,转换后得到的目标数据是RPS记录格式数据;步骤S1中产生波形数据的仪器包括:R24、S12、NZ24、Geopen、Geode、WZG系列仪器和/或ZDQ系列仪器和/或SWS系列仪器;
S2,数据预处理,具体为对S1中格式转换后得到的数据进行分量数据重排和编辑文件头;步骤S2包括以下步骤:
S21,数据编辑,即按钻孔深度将各个测点的地震剪切波数据以正向和反向分别进行多文件抽道、经道重排为XYZ三个分量的数据;
S22,编辑文件头,即编辑重排后数据文件头参数,适当调节各分量正反向激发的地震波的震幅;所述文件头参数包括:孔号、激发点偏距、接收点深度、系统延时;
S3,建立参数文件,具体为建立分析计算所需的解译参数文件;步骤S3具体为:根据步骤S2所得到的各分量数据,选择数据建立解译参数文件;采用正反方向激发、XYZ三个分量地震波数据建立参数文件;根据预处理后地震波数据的质量、激发方式因素,建立解译参数文件;
S4,计算分析,具体为根据剪切波的特征对比分析并正确拾取地震剪切波初至走时,自动计算各分量的测试结果;步骤S4包括以下步骤:
S41,拾取地震剪切波初至走时,即根据剪切波的特征,在波列中对比分析、识别剪切波,并拾取地震剪切波的初至走时;在同时显示正反向激发的X(或Y或Z)分量地震波形窗口,根据正反向激发的剪切波相互反相的特征对比分析,准确判别地震剪切波,并拾取初至走时,同步计算出各测试段的剪切波速;
S42,计算剪切波速,即根据测点间距和初至走时自动计算各分量的测试结果;将拾取地震剪切波的初至走时与计算结果动态关联,快速检查和修改地震剪切波初至走时和计算结果,直至得到合理的成果;通过切换显示出XYZ三个分量地震波形,正反方向激发叠加显示,拾取各分量地震剪切波初至走时,并同步计算出各分量的测试结果地震剪切波速;拾取剪切波初至走时和测试结果动态关联;
S5,成果取值和输出,具体为对S4得出的测试结果进行取值并输出成果;步骤S5具体为:在XYZ三个分量测试结果中,通过分析合理综合取值,确定最终成果并输出;XYZ三个分量的测试结果因为受到地震波数据质量的影响,在不同的测试深度上的可靠度会有所不同,本方法就是通过综合分析,合理选择某一分量或选择不同分量可靠测试段组合、或选择各分量的平均值作为最终成果,并输出。
2.一种地震剪切波的准确判别系统,其特征在于,包括:
数据格式转换模块,用于将不同型号仪器记录的不同格式的波形数据转换为预定的专有数据格式;被转换的源数据包括SEG-2格式数据和/或CSP格式数据,转换后得到的目标数据是RPS记录格式数据;产生波形数据的仪器包括:R24、S12、NZ24、Geopen、Geode、WZG系列仪器和/或ZDQ系列仪器和/或SWS系列仪器;
数据预处理模块,用于对格式转换后得到的数据进行分量数据重排和编辑文件头;具体用于:数据编辑,即按钻孔深度将各个测点的地震剪切波数据以正向和反向分别进行多文件抽道、经道重排为XYZ三个分量的数据;编辑文件头,即编辑重排后数据文件头参数,适当调节各分量正反向激发的地震波的震幅;所述文件头参数包括:孔号、激发点偏距、接收点深度、系统延时;
参数文件建立模块,用于建立分析计算所需的解译参数文件;具体用于根据所述数据预处理模块所得到的各分量数据,选择数据建立解译参数文件;采用正反方向激发、XYZ三个分量地震波数据建立参数文件;根据预处理后地震波数据的质量、激发方式因素,建立解译参数文件;
计算分析模块,用于根据剪切波的特征对比分析并正确拾取地震剪切波初至走时,自动计算各分量的测试结果;具体用于:拾取地震剪切波初至走时,即根据剪切波的特征,在波列中对比分析、识别剪切波,并拾取地震剪切波的初至走时;在同时显示正反向激发的X(或Y或Z)分量地震波形窗口,根据正反向激发的剪切波相互反相的特征对比分析,准确判别地震剪切波,并拾取初至走时,同步计算出各测试段的剪切波速;计算剪切波速,即根据测点间距和初至走时自动计算各分量的测试结果;将拾取地震剪切波的初至走时与计算结果动态关联,快速检查和修改地震剪切波初至走时和计算结果,直至得到合理的成果;通过切换显示出XYZ三个分量地震波形,正反方向激发叠加显示,拾取各分量地震剪切波初至走时,并同步计算出各分量的测试结果地震剪切波速;拾取剪切波初至走时和测试结果动态关联;
成果取值和输出模块,用于对得出的测试结果进行取值并输出成果;具体用于:在XYZ三个分量测试结果中,通过分析合理综合取值,确定最终成果并输出;XYZ三个分量的测试结果因为受到地震波数据质量的影响,在不同的测试深度上的可靠度会有所不同,本方法就是通过综合分析,合理选择某一分量或选择不同分量可靠测试段组合、或选择各分量的平均值作为最终成果,并输出。
CN201510107000.6A 2015-03-11 2015-03-11 一种地震剪切波的准确判别方法和系统 Active CN104656135B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510107000.6A CN104656135B (zh) 2015-03-11 2015-03-11 一种地震剪切波的准确判别方法和系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510107000.6A CN104656135B (zh) 2015-03-11 2015-03-11 一种地震剪切波的准确判别方法和系统

Publications (2)

Publication Number Publication Date
CN104656135A CN104656135A (zh) 2015-05-27
CN104656135B true CN104656135B (zh) 2017-09-01

Family

ID=53247482

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510107000.6A Active CN104656135B (zh) 2015-03-11 2015-03-11 一种地震剪切波的准确判别方法和系统

Country Status (1)

Country Link
CN (1) CN104656135B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115267883B (zh) * 2022-06-23 2023-08-29 华南理工大学 地震响应预测模型训练及预测方法、系统、设备及介质

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102830170A (zh) * 2012-07-23 2012-12-19 中国科学院地质与地球物理研究所 一种基于超声测试获取煤样横波信号的控制方法及装置
CN102968404A (zh) * 2012-11-02 2013-03-13 中国石油天然气股份有限公司 地震资料的初至信息处理方法及装置
CN104166157A (zh) * 2013-05-16 2014-11-26 中国石油化工股份有限公司 基于特殊条件下炸药震源激发三分量接收的横波处理方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2008087505A2 (en) * 2007-01-20 2008-07-24 Spectraseis Ag Time reverse reservoir localization

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102830170A (zh) * 2012-07-23 2012-12-19 中国科学院地质与地球物理研究所 一种基于超声测试获取煤样横波信号的控制方法及装置
CN102968404A (zh) * 2012-11-02 2013-03-13 中国石油天然气股份有限公司 地震资料的初至信息处理方法及装置
CN104166157A (zh) * 2013-05-16 2014-11-26 中国石油化工股份有限公司 基于特殊条件下炸药震源激发三分量接收的横波处理方法

Also Published As

Publication number Publication date
CN104656135A (zh) 2015-05-27

Similar Documents

Publication Publication Date Title
CN104793264B (zh) 应用于钻机的地质状况实时反映与超前探测系统及方法
CN104216008B (zh) 一种井中压裂微地震事件识别方法
CN107272061B (zh) 一种次声信号与地震事件的自动关联方法
CN106353792A (zh) 一种适用于水力压裂微震震源定位的方法
US11789173B1 (en) Real-time microseismic magnitude calculation method and device based on deep learning
CN109765606B (zh) 一种基于反射槽波的回采工作面隐蔽断层性质的探测方法
CN104570110B (zh) 一种基于纵横波匹配的多分量资料联合速度分析方法
CN105719433B (zh) 一种基于孔内地震波的超前预报方法
CN110133715A (zh) 一种基于初至时差和波形叠加的微地震震源定位方法
CN108797662B (zh) 立柱下桩基长度无损检测方法及装置
CN107807382A (zh) 一种高密度面波勘探数据采集与处理方法
CN105735971B (zh) 一种基于弹性波的钻孔深度检测系统及其检测方法
CN103726836A (zh) 基于声波测井资料提取模式波慢度的方法
CN104569158A (zh) 基于爆破振动测试的岩体质量分类及动力参数估计方法
CN110058323A (zh) 一种致密砂岩地层脆性指数计算方法
CN106599377A (zh) 基于测井资料和煤体结构参量定量划分煤体结构的方法
CN106885846A (zh) 树木缺陷检测装置及检测方法
CN113034043A (zh) 一种基于地理特征处理技术的道路工程勘察测绘智能管理系统
CN105445793B (zh) 一种确定坏道数据的方法及装置
CN104656135B (zh) 一种地震剪切波的准确判别方法和系统
CN107703538A (zh) 地下不良地质勘测数据采集分析系统及方法
Dickmann et al. Towards the integration of smart techniques for tunnel seismic applications
CN106646610A (zh) 一种利用偏振约束aic算法自动拾取微震初至的算法
US20240134080A1 (en) Method and System for Real-Time Calculating a Microseismic Focal Mechanism Based on Deep Learning
CN105866834B (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