CN112083481A - 一种通过加速度获取速度和位移的方法 - Google Patents
一种通过加速度获取速度和位移的方法 Download PDFInfo
- Publication number
- CN112083481A CN112083481A CN202010951773.3A CN202010951773A CN112083481A CN 112083481 A CN112083481 A CN 112083481A CN 202010951773 A CN202010951773 A CN 202010951773A CN 112083481 A CN112083481 A CN 112083481A
- Authority
- CN
- China
- Prior art keywords
- acceleration
- record
- earthquake
- displacement
- time
- 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.)
- Pending
Links
- 230000001133 acceleration Effects 0.000 title claims abstract description 95
- 238000000034 method Methods 0.000 title claims abstract description 40
- 238000006073 displacement reaction Methods 0.000 title claims abstract description 33
- 238000001914 filtration Methods 0.000 claims abstract description 13
- 238000004364 calculation method Methods 0.000 claims abstract description 6
- 238000005070 sampling Methods 0.000 claims description 6
- 230000000694 effects Effects 0.000 abstract description 6
- 238000005457 optimization Methods 0.000 abstract description 3
- 238000010586 diagram Methods 0.000 description 4
- 230000010354 integration Effects 0.000 description 2
- 230000004075 alteration Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Images
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/24—Recording seismic data
- G01V1/242—Seismographs
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01P—MEASURING LINEAR OR ANGULAR SPEED, ACCELERATION, DECELERATION, OR SHOCK; INDICATING PRESENCE, ABSENCE, OR DIRECTION, OF MOVEMENT
- G01P15/00—Measuring acceleration; Measuring deceleration; Measuring shock, i.e. sudden change of acceleration
- G01P15/18—Measuring acceleration; Measuring deceleration; Measuring shock, i.e. sudden change of acceleration in two or more dimensions
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/16—Receiving elements for seismic signals; Arrangements or adaptations of receiving elements
- G01V1/18—Receiving elements, e.g. seismometer, geophone or torque detectors, for localised single point measurements
-
- 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
-
- 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
-
- 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/303—Analysis for determining velocity profiles or travel times
Abstract
本发明提供一种通过加速度获取速度和位移的方法,按照先后顺序分别包括以下5个部分:先获取加速度记录A的震前时长t(1),对加速度记录A进行基线校正处理得到加速度记录A1(2),对加速度记录A1进行巴特沃斯带通滤波处理得到加速度记录A2(3),对加速度记录A2积分计算速度时程V(4),以及对速度时程积分计算位移时程S(5)。采用本发明的技术方案,通过波形图对比,明显看出经过本发明方法计算之后,速度漂移和位移漂移效果不明显,基线漂移优化效果明显。
Description
技术领域
本发明涉及地震动数据处理技术领域,具体涉及一种通过加速度获取速度和位移的方法。
背景技术
强震仪是记录强烈地震近地面运动的自动触发式地震仪,被广泛的用于地震台网、流动台、核电站、大坝、水库等地方。强震仪记录的是三个方向(垂直、东西、南北)的加速度数据,记录时不仅包含了地震时地面运动信息,还包含有复杂的噪音,其中的低频噪音会导致加速度时程出现基线漂移。
基线漂移对加速度时程本身可能影响很小(一般不超过峰值加速度的2%),但通过积分求速度时程时,基线漂移被放大;而通过积分求位移时程时,基线漂移被进一步放大,往往与真实的位移时程相差甚远。因此,在使用加速度时程计算速度时程和位移时程,一般需要对基线漂移进行处理。
发明内容
有鉴于此,本发明目的是提供一种通过加速度获取速度和位移的方法,采用基线校正和巴特沃斯带通滤波加速度原始记录进行处理,然后对处理的加速度记录进行积分得到速度时程,再对速度时程进行积分得到位移时程。
为实现上述目的,本发明提供如下技术方案:一种通过加速度获取速度和位移的方法,按照先后顺序分别包括以下5个部分:先获取加速度记录A的震前时长t(1),对加速度记录A进行基线校正处理得到加速度记录A1(2),对加速度记录A1进行巴特沃斯带通滤波处理得到加速度记录A2(3),对加速度记录A2积分计算速度时程V(4),以及对速度时程积分计算位移时程S(5)。
优选地,获取加速度记录A的震前时长t(1)的具体步骤如下:一般而言,一台强震仪的记录下来的加速度记录的震前时长t是固定的,这样将每台强震仪的震前时长存储到配置中,该配置可以存储在文件中、数据中或者内存中。从配置中获取震前时长t,如果震前时长t不存在,则需要根据加速度记录A计算出震前时长,然后将震前时长t存储到配置中,以方便下次计算时使用。
优选地,对加速度记录A进行基线校正处理得到加速度记录A1(2)的具体步骤如下:根据震前时长t计算出震前记录平均值avg,再将加速度记录A减去该平均值avg,得到加速度记录A1;
式中,fs为加速度记录的采样率,n为震前记录长度。
优选地,对加速度记录A1进行巴特沃斯带通滤波处理得到加速度记录A2(3)的具体步骤如下:为了过滤掉加速度记录中的低频噪音,这里采用巴特沃斯滤波算法对加速度记录A1进行0.01Hz至20Hz的带通滤波,滤波的阶数为4,得到加速度记录A2。
优选地,对加速度记录A2积分计算速度时程V(4)的具体步骤如下:
加速度记录A2是离散的数据,基于时间上是等间隔,强震仪的采样频率为100Hz、200H z,时间间隔为0.01s、0.005s,相邻数据之间时间间隔足够小,所以这里采用二阶龙格-库塔法(Runge-Kutta methods),对加速度记录A2进行计算得到速度时程V;
其差分方程为:
优选地,对速度时程积分计算位移时程S(5)的具体步骤如下:同上述步骤(4)一样使用二阶龙格-库塔法,对速度时程V计算得到位移时程S。
与现有技术相比,本发明的有益效果为:
采用本发明的技术方案,通过加速度、速度以及位移的波形图对比,明显看出经过本发明方法计算之后,速度漂移和位移漂移效果不明显,基线漂移优化效果明显。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图;
图1本发明的整体方法中流程示意图;
图2采用本发明整体方法对应加速度时程波形图;
图3作为图2对照组的加速度时程波形图;
图4采用本发明整体方法对应速度时程波形图;
图5作为图4对照组的速度时程波形图;
图6采用本发明整体方法对位移时程波形图;
图7作为图6对照组的位移时程波形图。
具体实施方式
下面结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明的保护范围。
如图1~7所示,本发明提供一种通过加速度获取速度和位移的方法,按照先后顺序分别包括以下5个部分:先获取加速度记录A的震前时长t(1),对加速度记录A进行基线校正处理得到加速度记录A1(2),对加速度记录A1进行巴特沃斯带通滤波处理得到加速度记录A2(3),对加速度记录A2积分计算速度时程V(4),以及对速度时程积分计算位移时程S(5)。
其中,获取加速度记录A的震前时长t(1)的具体步骤如下:一般而言,一台强震仪的记录下来的加速度记录的震前时长t是固定的,这样将每台强震仪的震前时长存储到配置中,该配置可以存储在文件中、数据中或者内存中。从配置中获取震前时长t,如果震前时长t不存在,则需要根据加速度记录A计算出震前时长,然后将震前时长t存储到配置中,以方便下次计算时使用。
对加速度记录A进行基线校正处理得到加速度记录A1(2)的具体步骤如下:根据震前时长t计算出震前记录平均值avg,再将加速度记录A减去该平均值avg,得到加速度记录A1;
式中,fs为加速度记录的采样率,n为震前记录长度。
对加速度记录A1进行巴特沃斯带通滤波处理得到加速度记录A2(3)的具体步骤如下:为了过滤掉加速度记录中的低频噪音,这里采用巴特沃斯滤波算法对加速度记录A1进行0.01Hz至20Hz的带通滤波,滤波的阶数为4,得到加速度记录A2。
对加速度记录A2积分计算速度时程V(4)的具体步骤如下:
加速度记录A2是离散的数据,基于时间上是等间隔,强震仪的采样频率为100Hz、200H z,时间间隔为0.01s、0.005s,相邻数据之间时间间隔足够小,所以这里采用二阶龙格-库塔法(Runge-Kutta methods),对加速度记录A2进行计算得到速度时程V;
其差分方程为:
对速度时程积分计算位移时程S(5)的具体步骤如下:同上述步骤(4)一样使用二阶龙格-库塔法,对速度时程V计算得到位移时程S。
采用本发明的技术方案,通过波形图对比,明显看出经过本发明方法计算之后,速度漂移和位移漂移效果不明显,基线漂移优化效果明显。
需要说明的是,在本文中,诸如第一和第二等之类的关系术语仅仅用来将一个实体或者操作与另一个实体或操作区分开来,而不一定要求或者暗示这些实体或操作之间存在任何这种实际的关系或者顺序。而且,术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含,从而使得包括一系列要素的过程、方法、物品或者设备不仅包括那些要素,而且还包括没有明确列出的其他要素,或者是还包括为这种过程、方法、物品或者设备所固有的要素。
尽管已经示出和描述了本发明的实施例,对于本领域的普通技术人员而言,可以理解在不脱离本发明的原理和精神的情况下可以对这些实施例进行多种变化、修改、替换和变型,本发明的范围由所附权利要求及其等同物限定。
Claims (6)
1.一种通过加速度获取速度和位移的方法,其特征在于,按照先后顺序分别包括以下5个部分:先获取加速度记录A的震前时长t(1),对加速度记录A进行基线校正处理得到加速度记录A1(2),对加速度记录A1进行巴特沃斯带通滤波处理得到加速度记录A2(3),对加速度记录A2积分计算速度时程V(4),以及对速度时程积分计算位移时程S(5)。
2.根据权利要求1中所述的一种通过加速度获取速度和位移的方法,其特征在于,获取加速度记录A的震前时长t(1)的具体步骤如下:
一般而言,一台强震仪的记录下来的加速度记录的震前时长t是固定的,这样将每台强震仪的震前时长存储到配置中,该配置存储在文件中、数据中或者内存中;从配置中获取震前时长t,如果震前时长t不存在,则需要根据加速度记录A计算出震前时长,然后将震前时长t存储到配置中,以方便下次计算时使用。
4.根据权利要求3中所述的一种通过加速度获取速度和位移的方法,其特征在于,对加速度记录A1进行巴特沃斯带通滤波处理得到加速度记录A2(3)的具体步骤如下:为了过滤掉加速度记录中的低频噪音,这里采用巴特沃斯滤波算法对加速度记录A1进行0.01Hz至20Hz的带通滤波,滤波的阶数为4,得到加速度记录A2。
6.根据权利要求5中所述的一种通过加速度获取速度和位移的方法,其特征在于,
对速度时程积分计算位移时程S(5)的具体步骤如下:同上述步骤(4)一样使用二阶龙格-库塔法,对速度时程V计算得到位移时程S。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010951773.3A CN112083481A (zh) | 2020-09-11 | 2020-09-11 | 一种通过加速度获取速度和位移的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010951773.3A CN112083481A (zh) | 2020-09-11 | 2020-09-11 | 一种通过加速度获取速度和位移的方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN112083481A true CN112083481A (zh) | 2020-12-15 |
Family
ID=73736675
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010951773.3A Pending CN112083481A (zh) | 2020-09-11 | 2020-09-11 | 一种通过加速度获取速度和位移的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112083481A (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113009566A (zh) * | 2021-03-19 | 2021-06-22 | 中国地震局地球物理研究所 | 现地地震动预测模型及其构建方法 |
CN113049202A (zh) * | 2021-03-08 | 2021-06-29 | 中国地震局工程力学研究所 | 一种加速度积分位移的局部加权回归校正方法及系统 |
CN113465732A (zh) * | 2021-08-13 | 2021-10-01 | 重庆大学 | 基于emd-dnn的振动台结构位移响应预测方法及装置 |
CN115876413A (zh) * | 2023-02-14 | 2023-03-31 | 济南城建集团有限公司 | 基于加速度的移动车辆作用下桥梁位移估计方法 |
CN117607967A (zh) * | 2024-01-19 | 2024-02-27 | 中国建筑西南设计研究院有限公司 | 一种基于遗传算法的地震动基线校正方法及电子设备 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108181486A (zh) * | 2018-01-25 | 2018-06-19 | 中国科学院电子学研究所 | 加速度信号的处理方法及装置 |
CN109143370A (zh) * | 2018-07-25 | 2019-01-04 | 中国地震局地球物理研究所 | 地震动加速度记录基线漂移的校正方法 |
CN110657935A (zh) * | 2019-09-24 | 2020-01-07 | 浙江大学 | 一种地震波加速度积分处理方法及系统 |
-
2020
- 2020-09-11 CN CN202010951773.3A patent/CN112083481A/zh active Pending
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108181486A (zh) * | 2018-01-25 | 2018-06-19 | 中国科学院电子学研究所 | 加速度信号的处理方法及装置 |
CN109143370A (zh) * | 2018-07-25 | 2019-01-04 | 中国地震局地球物理研究所 | 地震动加速度记录基线漂移的校正方法 |
CN110657935A (zh) * | 2019-09-24 | 2020-01-07 | 浙江大学 | 一种地震波加速度积分处理方法及系统 |
Non-Patent Citations (2)
Title |
---|
曲中宪 等: "《数值分析》", 31 August 2004, 东北大学出版社 * |
朱明周: "岩土台阵的土动力学参数与场地反应研究", 《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》 * |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113049202A (zh) * | 2021-03-08 | 2021-06-29 | 中国地震局工程力学研究所 | 一种加速度积分位移的局部加权回归校正方法及系统 |
CN113049202B (zh) * | 2021-03-08 | 2022-07-12 | 中国地震局工程力学研究所 | 一种加速度积分位移的局部加权回归校正方法及系统 |
CN113009566A (zh) * | 2021-03-19 | 2021-06-22 | 中国地震局地球物理研究所 | 现地地震动预测模型及其构建方法 |
CN113009566B (zh) * | 2021-03-19 | 2024-01-30 | 中国地震局地球物理研究所 | 现地地震动预测模型及其构建方法 |
CN113465732A (zh) * | 2021-08-13 | 2021-10-01 | 重庆大学 | 基于emd-dnn的振动台结构位移响应预测方法及装置 |
CN115876413A (zh) * | 2023-02-14 | 2023-03-31 | 济南城建集团有限公司 | 基于加速度的移动车辆作用下桥梁位移估计方法 |
CN117607967A (zh) * | 2024-01-19 | 2024-02-27 | 中国建筑西南设计研究院有限公司 | 一种基于遗传算法的地震动基线校正方法及电子设备 |
CN117607967B (zh) * | 2024-01-19 | 2024-03-26 | 中国建筑西南设计研究院有限公司 | 一种基于遗传算法的地震动基线校正方法及电子设备 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112083481A (zh) | 一种通过加速度获取速度和位移的方法 | |
CN107272066B (zh) | 一种含噪地震信号初至走时拾取方法及装置 | |
CN104586384B (zh) | 一种心电信号实时心率检测方法及系统 | |
CN113009566B (zh) | 现地地震动预测模型及其构建方法 | |
CN108601543B (zh) | 一种ecg信号处理方法及装置 | |
CN109557587B (zh) | 一种vsp地震资料井筒波频率域滤波方法及装置 | |
CN105510976A (zh) | 一种多次波组合自适应衰减方法 | |
CN107479094A (zh) | 一种实现地震预警的方法 | |
CN107144879A (zh) | 一种基于自适应滤波与小波变换结合的地震波降噪方法 | |
CN112904414A (zh) | 一种地声事件定位方法及其失稳灾害预警方法、地声感知仪、监测系统及其可读存储介质 | |
CN109425897B (zh) | 消除地震资料野值干扰的方法及系统 | |
CN106073783A (zh) | 一种从光电容积脉搏波中提取呼吸率的方法 | |
CN104656140B (zh) | 一种基于垂直地震衰减规律的中值滤波方法以及系统 | |
CN110019175A (zh) | 一种基于行为轨迹库的区域碰撞分析方法 | |
CN102085094A (zh) | 柯氏音始发点与消失点的声强峰值斜率确定法 | |
CN109741762B (zh) | 声音活动检测方法及装置和计算机可读存储介质 | |
US11170760B2 (en) | Detecting speech activity in real-time in audio signal | |
CN109298451B (zh) | 一种改进偏斜度的自动拾取s波震相方法 | |
CN116383628A (zh) | 一种落石灾害微震信号的检测分类方法及系统 | |
CN109602414B (zh) | 一种多视角转换的心电信号数据增强方法 | |
CN106529425A (zh) | 一种心电信号的r波提取方法及系统 | |
CN111292748A (zh) | 一种可适应多种频率的语音录入系统 | |
CN106772582A (zh) | 一种基于图像轮廓分析的地震数据重构方法 | |
LU504061B1 (en) | Frequency domain filtering method and device of wellbore wave in vsp seismic data | |
CN111965702B (zh) | 一种基于地震初至的剩余静校正量处理方法及系统 |
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 | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20201215 |
|
RJ01 | Rejection of invention patent application after publication |