CN109270575A - 一种基于建筑物地震响应等效的爆破地震波模型构造方法 - Google Patents

一种基于建筑物地震响应等效的爆破地震波模型构造方法 Download PDF

Info

Publication number
CN109270575A
CN109270575A CN201811310006.3A CN201811310006A CN109270575A CN 109270575 A CN109270575 A CN 109270575A CN 201811310006 A CN201811310006 A CN 201811310006A CN 109270575 A CN109270575 A CN 109270575A
Authority
CN
China
Prior art keywords
attenuation
signal
seismic waves
seismic
amplitude
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
CN201811310006.3A
Other languages
English (en)
Other versions
CN109270575B (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.)
Henan University of Technology
Original Assignee
Henan University of Technology
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 Henan University of Technology filed Critical Henan University of Technology
Priority to CN201811310006.3A priority Critical patent/CN109270575B/zh
Publication of CN109270575A publication Critical patent/CN109270575A/zh
Application granted granted Critical
Publication of CN109270575B publication Critical patent/CN109270575B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28Processing seismic data, e.g. for interpretation or for event detection
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/282Application of seismic models, synthetic seismograms
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/30Analysis

Landscapes

  • Engineering & Computer Science (AREA)
  • Remote Sensing (AREA)
  • Physics & Mathematics (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)
  • Buildings Adapted To Withstand Abnormal External Influences (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

本发明公开了一种基于建筑物地震响应等效的爆破地震波模型构造方法,包括如下步骤:选择实测爆破地震波速度信号作为待模拟信号;将爆破地震波速度信号转换为加速度信号;计算建筑结构分析对象在加速度荷载作用下的结构响应,并记录顶层速度响应的幅值;将爆破地震波速度信号进行多分辨小波分解和重构,确定能量比例达10%以上的重构子信号的主频和幅值;构造待模拟信号的爆破地震波模型:采用同一时程分析法,计算与顶层速度响应幅值结果相同时的a的具体取值;将求出的具体a值代入所构造的模型中,得到基于建筑物地震响应等效的爆破地震波模拟信号的具体表达式,得到的爆破地震波模型不仅具有实测波形的特性参数,还能适用于爆破地震效应机理研究。

Description

一种基于建筑物地震响应等效的爆破地震波模型构造方法
技术领域
本发明属于爆破地震效应评价和工程结构抗震领域,尤其涉及 一种基于建筑物地震响应等效的爆破地震波模型构造方法。
背景技术
爆破地震波信号作为动力荷载输入来研究爆破地震波对工程结 构的地震效应作用机理时,通常有两种形式:实测爆破地震波和人 工模拟爆破地震波。由于实测爆破地震波信号的离散性和随机性, 基于数值分析方法求解的爆破地震效应结果规律性不明显,且很难 进行爆破地震效应机理相关的结构动力方程推导,因而一般求不出 具体问题对应的解析解。爆破地震波模拟信号虽能够按照期望条件 人工合成具有不同主频、振幅、持续时间等主要特性参数的波形, 但即使合成的爆破地震波信号在特性参数和形态上高度相似,而加 载到特定结构后所输出的结构地震响应幅值差异较大。
发明内容
本申请旨在至少解决现有技术中存在的技术问题之一。为此, 本发明的目的之一在于提供一种基于建筑物地震响应等效的爆破地 震波模型构造方法,利用该方构造的爆破地震波模型不仅具有实测 波形的特性参数,还能适用于爆破地震效应机理研究。
为解决上述技术问题,本发明采用如下技术方案:
一种基于建筑物地震响应等效的爆破地震波模型构造方法,包 括如下步骤:
S1:选择实测爆破地震波速度信号v0(t)作为待模拟信号;
S2:将爆破地震波速度信号v0(t)转换为加速度信号a0(t);
S3:利用时程分析法计算建筑结构分析对象在加速度荷载作用 下的结构响应,并记录顶层速度响应的幅值;
S4:将爆破地震波速度信号v0(t)进行多分辨小波分解和重构, 确定能量比例达10%以上的重构子信号的主频fi和幅值Ai
S5:构造待模拟信号的爆破地震波模型,如下式所示:
其中:
ωi为主频fi对应的角频率,ωi=2πfi
Ag为地震波调幅系数,为保持模拟信号与待模拟爆破地震波信 号的幅值A0相等,取n为满足条件的重构子信号数目;
e-at为衰减函数,其中参数a为一常数;
S6:采用S3步骤中的同一时程分析法,计算参数a取不同值时 同一建筑结构分析对象在S5步骤所构造的爆破地震波模型作用下 的顶层速度响应幅值,记录与S3步骤中的顶层速度响应幅值结果 相同时的a的具体取值;
求解构造的爆破地震波模型作用下的结构响应时,需要先通过 对公式(3)一阶求导得出加速度荷载模型,然后通过某微小的时间 步长将模型离散成一系列的加速度信号序列。
S7:将S6步骤中求出的具体a值代入S5步骤所构造的模型中, 得到基于建筑物地震响应等效的爆破地震波模拟信号的具体表达 式。
进一步的,爆破地震波速度信号v0(t)通过如下公式转换为加速 度信号a0(t),
式中,Δt为采集速度列的时间步长,ti为第i个时间步长对应的时刻, a0(ti)为ti时刻对应的加速度值,k为采集速度列的数据点总数。
进一步的,步骤S4中主频fi和幅值Ai的具体求解过程为:分解 层数以分解后的频带宽度达20Hz以内来确定,假定爆破地震波信 号v0(t)分解层数取m时符合要求,则信号v0(t)在经小波分解和重构后 可以得到m+1个子频带,以vi(t)表示第i个子频带上的单支重构信号, 则总信号可表示为:
然后求出每个频带上重构子信号的主频、幅值及能量比例,从而确 定出能量比例达10%以上的重构子信号vvi的主频fi和幅值Ai,其中, i=1,2…..n,n为满足条件的重构子信号数目。
进一步的,利用时程分析法中的Wilson-θ法计算建筑结构分析 对象在加速度荷载作用下的结构响应。
进一步的,小波基函数取db8。
与现有技术相比,本发明的有益效果如下:
本发明构造出了与原爆破地震波实测信号具有相同幅值、持续 时间和多个优势频率的爆破地震波模拟信号,且模拟信号同时还能 一定程度上反映了被分析结构对象所在的场地条件特征,为爆破地 震效应的机理研究提供适应性强的爆破地震波模拟信号,为推导爆 破地震效应相关具体问题对应的解析解提供可行性。
综上所述,本发明对爆破地震效应安全评估、预测和控制具有 重要的理论意义,对建筑物附近的爆破工程方案设计和优化具有指 导作用。
附图说明
图1为本发明方法的流程示意图;
图2为本发明的实测爆破地震波信号;
图3为本发明的基于建筑物地震响应等效的爆破地震波模拟信 号。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术 方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明 一部分实施例,而不是全部的实施例。基于本发明中的实施例,本 领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他 实施例,都属于本发明保护的范围。
一种基于建筑物地震响应等效的爆破地震波模型构造方法,包 括如下步骤:
步骤一:为保证构造的爆破地震波模型能反映场地条件特征, 在爆区实测的爆破地震波中选择典型的待模拟地震波信号。实测的 爆破地震波信号一般是速度时程波,导出后是一系列离散的速度序 列数据,记为v0(t)。
步骤二:将实测的爆破地震波速度信号转换成加速度信号。将 爆破地震波速度信号v0(t)通过公式(1)转换为加速度信号(记为 a0(t))。公式(1)如下:
式中,Δt为采集速度列的时间步长,ti为第i个时间步长对应的时刻,a0(ti)为ti时刻对应的加速度值,k为采集速度列的数据点总数。
步骤三:利用时程分析法计算特定的建筑结构分析对象在加速 度荷载作用下的结构响应,并记录顶层速度响应的幅值。
步骤四:将待模拟爆破地震波信号进行多分辨小波分解和重构, 确定能量比例达10%以上的重构子信号的主频和幅值,分解层数以 分解后的频带宽度达20Hz以内来确定,假定爆破地震波信号v0(t)分 解层数取m时符合要求,则信号v0(t)在经小波分解和重构后可以得 到m+1个子频带。若以vi(t)表示第i个子频带上的单支重构信号,则 总信号可表示为:
然后求出每个频带上重构子信号的主频、幅值及能量比例,从而确 定出能量比例达10%以上的重构子信号vvi(i=1,2…..n,n为满足 条件的重构子信号数目)的主频(记为fi)和幅值(记为Ai)。
步骤五:将待模拟爆破地震波信号的模型构造成如下形式:
其中,ωi为能量比例达10%以上的重构子信号vvi的主频fi对应的角 频率,ωi=2πfi;Ag为地震波调幅系数,为保持模拟信号与待模拟爆 破地震波信号的幅值(记为A0)相等,取 为函数式 所构造模型的幅值;e-at为衰减函数,其中参数a为一常 数,具体值由如下所述的建筑物地震响应等效的方法来确定。
步骤六:采用步骤三中的同一时程分析法,计算参数a取不同 值时同一建筑结构分析对象在步骤五所构造的爆破地震波模型作用 下的顶层速度响应幅值,记录与步骤三中的顶层速度响应幅值结果 相同时的a的具体取值。注意,求解构造的爆破地震波模型作用下 的结构响应时,需要先通过对公式(3)一阶求导得出加速度荷载模 型,然后通过某微小的时间步长将模型离散成一系列的加速度信号 序列。
步骤七:最终确定出基于建筑物地震响应等效的爆破地震波模 拟信号的具体模型。将步骤六中求出的具体a值代入步骤五所构造 的模型,即可构造出基于建筑物地震响应等效的爆破地震波模拟信 号的具体函数表达式。
实施例
一种基于建筑物地震响应等效的爆破地震波模型构造方法,包 括如下步骤:
步骤一:为保证构造的爆破地震波模型能反映场地条件特征, 在爆区实测的爆破地震波中选择一条典型的地震波速度信号v0(t)作 为待模拟信号,其速度时程曲线见图2。实测的爆破地震波信号一 般是速度时程波,导出后是一系列离散的速度序列数据,记为v0(t)。
步骤二:将爆破地震波速度信号v0(t)通过公式(1)转换成加速度 信号a0(t)。
步骤三:取某三层钢筋混凝土建筑结构作为分析对象,结构的 特性参数为:各层质量m=[2762,2760,2300]kg,各层刚度 k=[2.485,1]×95N/m,采用时程分析法中的Wilson-θ法(取 θ=1.4)计算该结构在加速度信号a0(t)的作用下的地震响应,并记录 下顶层速度响应的幅值为2.73cm/s。
步骤四:将模拟信号v0(t)进行多分辨小波分解和重构,分解层 数取11层,小波基函数取db8,经多分辨小波分解和重构后的12 个子频带上子信号的主频、幅值及能量比例结果见表1,通过表1 确定出能量比例达10%以上的4个重构子信号的主频和幅值分别 为:63.52Hz,0.81cm/s;122.87Hz,1.46cm/s;213.45Hz,0.95cm/s; 582.05Hz,1.24cm/s。
步骤五:将待模拟爆破地震波信号的模型构造成如下形式:
V(t)=0.5224e-at[0.81cos(399t)+1.46cos(772t)+0.95cos(1341t)+1.24cos(3657t)]。
步骤六:采用Wilson-θ时程分析法,计算参数a取不同值时同 一建筑结构在步骤五所构造的爆破地震波模型作用下的顶层速度响 应幅值,最终确定出顶层速度响应幅值取2.73cm/s时a的取值为 0.05。
步骤七:将a=0.05代入步骤五所构造的模型,最终确定出基于 该建筑物地震响应等效的爆破地震波模拟信号的具体模型为:
V(t)=0.5224e-0.05t[0.81cos(399t)+1.46cos(772t)+0.95cos(1341t)+1.24cos(3657t)]。
该模拟信号的速度时程曲线见图3。
表1为本发明具体实施例所述记录的爆破地震波信号多分辨小 波分析结果。
子频带(Hz) 主频(Hz) 幅值(cm/s) 能量比例(%)
0~1.2208 1.04 0.05 0.19
1.2208~2.4415 2.08 0.03 0.17
2.4415~4.8829 4.16 0.13 0.68
4.8829~9.7657 7.29 0.17 1.07
9.7657~19.5313 16.66 0.19 1.12
19.5313~39.0625 22.91 0.52 7.16
39.0625~78.125 63.52 0.81 23.71
78.125~156.25 122.87 1.46 26.55
156.25~312.5 213.45 0.95 16.43
312.5~625 582.05 1.24 18.83
625~1250 667.43 0.52 4.05
1250~2500 1917.95 0.13 0.04
上述实施例仅仅是清楚地说明本发明所作的举例,而非对实施 方式的限定。对于所属领域的普通技术人员来说,在上述说明的基 础上还可以做出其它不同形式的变化或变动。这里也无需也无法对 所有的实施例予以穷举。而由此所引申出的显而易见的变化或变动 仍处于本发明的保护范围之中。

Claims (5)

1.一种基于建筑物地震响应等效的爆破地震波模型构造方法,其特征在于,包括如下步骤:
S1:选择实测爆破地震波速度信号v0(t)作为待模拟信号;
S2:将爆破地震波速度信号v0(t)转换为加速度信号a0(t);
S3:利用时程分析法计算建筑结构分析对象在加速度信号a0(t)作用下的结构响应,并记录顶层速度响应的幅值;
S4:将爆破地震波速度信号v0(t)进行多分辨小波分解和重构,确定能量比例达10%以上的重构子信号的主频fi和幅值Ai
S5:构造待模拟信号的爆破地震波模型,如下式所示:
其中:
ωi为主频fi对应的角频率,ωi=2πfi
Ag为地震波调幅系数,为保持模拟信号与待模拟爆破地震波信号的幅值A0相等,取n为满足条件的重构子信号数目;
e-at为衰减函数,其中参数a为一常数;
S6:采用S3步骤中的同一时程分析法,计算参数a取不同值时同一建筑结构分析对象在S5步骤所构造的爆破地震波模型作用下的顶层速度响应幅值,记录与S3步骤中的顶层速度响应幅值结果相同时的a的具体取值;
S7:将S6步骤中求出的具体a值代入S5步骤所构造的模型中,得到基于建筑物地震响应等效的爆破地震波模拟信号的具体表达式。
2.根据权利要求1所述的基于建筑物地震响应等效的爆破地震波模型构造方法,其特征在于,爆破地震波速度信号v0(t)通过如下公式转换为加速度信号a0(t):
式中,Δt为采集速度列的时间步长,ti为第i个时间步长对应的时刻,a0(ti)为ti时刻对应的加速度值,k为采集速度列的数据点总数。
3.根据权利要求1所述的基于建筑物地震响应等效的爆破地震波模型构造方法,其特征在于,步骤S4中主频fi和幅值Ai的具体求解过程为:分解层数以分解后的频带宽度达20Hz以内来确定,假定爆破地震波信号v0(t)分解层数取m时符合要求,则信号v0(t)在经小波分解和重构后可以得到m+1个子频带,以vi(t)表示第i个子频带上的单支重构信号,则总信号可表示为:
然后求出每个频带上重构子信号的主频、幅值及能量比例,从而确定出能量比例达10%以上的重构子信号vvi的主频fi和幅值Ai,其中,i=1,2…..n,n为满足条件的重构子信号数目。
4.根据权利要求1所述的基于建筑物地震响应等效的爆破地震波模型构造方法,其特征在于:利用时程分析法中的Wilson-θ法计算建筑结构分析对象在加速度荷载作用下的结构响应。
5.根据权利要求1所述的基于建筑物地震响应等效的爆破地震波模型构造方法,其特征在于:小波基函数取db8。
CN201811310006.3A 2018-11-02 2018-11-02 一种基于建筑物地震响应等效的爆破地震波模型构造方法 Active CN109270575B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811310006.3A CN109270575B (zh) 2018-11-02 2018-11-02 一种基于建筑物地震响应等效的爆破地震波模型构造方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811310006.3A CN109270575B (zh) 2018-11-02 2018-11-02 一种基于建筑物地震响应等效的爆破地震波模型构造方法

Publications (2)

Publication Number Publication Date
CN109270575A true CN109270575A (zh) 2019-01-25
CN109270575B CN109270575B (zh) 2019-11-26

Family

ID=65192125

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811310006.3A Active CN109270575B (zh) 2018-11-02 2018-11-02 一种基于建筑物地震响应等效的爆破地震波模型构造方法

Country Status (1)

Country Link
CN (1) CN109270575B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111829405A (zh) * 2020-07-16 2020-10-27 中铁十六局集团北京轨道交通工程建设有限公司 一种基于小波分析城区爆破安全控制的方法

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102183787A (zh) * 2011-03-07 2011-09-14 中国海洋石油总公司 一种基于地震记录变子波模型提高地震资料分辨率的方法
CN102253415A (zh) * 2011-04-19 2011-11-23 中国石油大学(华东) 基于裂缝等效介质模型的地震响应模式建立方法
WO2015014762A3 (en) * 2013-07-29 2015-04-02 Cgg Services Sa Method and device for the generation and application of anisotropic elastic parameters in horizontal transverse isotropic (hti) media
US20150134308A1 (en) * 2012-09-14 2015-05-14 Institute Of Geology And Geophysics, Chinese Academy Of Sciences Method and device for acquiring optimization coefficient, and related method and device for simulating wave field
CN105044771A (zh) * 2015-08-05 2015-11-11 北京多分量地震技术研究院 基于有限差分法的三维tti双相介质地震波场数值模拟方法
CN106154331A (zh) * 2016-06-29 2016-11-23 中国石油化工股份有限公司 正交介质地震波场模拟频散压制方法
CN106199697A (zh) * 2016-06-29 2016-12-07 中国石油化工股份有限公司 模拟微地震的弹性波正演方法
US20180210101A1 (en) * 2016-12-30 2018-07-26 China Petroleum & Chemical Corporation System and method for seismic inversion

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102183787A (zh) * 2011-03-07 2011-09-14 中国海洋石油总公司 一种基于地震记录变子波模型提高地震资料分辨率的方法
CN102253415A (zh) * 2011-04-19 2011-11-23 中国石油大学(华东) 基于裂缝等效介质模型的地震响应模式建立方法
US20150134308A1 (en) * 2012-09-14 2015-05-14 Institute Of Geology And Geophysics, Chinese Academy Of Sciences Method and device for acquiring optimization coefficient, and related method and device for simulating wave field
WO2015014762A3 (en) * 2013-07-29 2015-04-02 Cgg Services Sa Method and device for the generation and application of anisotropic elastic parameters in horizontal transverse isotropic (hti) media
CN105044771A (zh) * 2015-08-05 2015-11-11 北京多分量地震技术研究院 基于有限差分法的三维tti双相介质地震波场数值模拟方法
CN106154331A (zh) * 2016-06-29 2016-11-23 中国石油化工股份有限公司 正交介质地震波场模拟频散压制方法
CN106199697A (zh) * 2016-06-29 2016-12-07 中国石油化工股份有限公司 模拟微地震的弹性波正演方法
US20180210101A1 (en) * 2016-12-30 2018-07-26 China Petroleum & Chemical Corporation System and method for seismic inversion

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
刘文光 等: ""隔震结构动力弹塑性分析地震记录选择的波谱分类法研究"", 《建筑结构学报》 *
陈士海 等: ""基于小波包技术的爆破地震效应计算模型及安全判据研究"", 《爆炸与冲击》 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111829405A (zh) * 2020-07-16 2020-10-27 中铁十六局集团北京轨道交通工程建设有限公司 一种基于小波分析城区爆破安全控制的方法

Also Published As

Publication number Publication date
CN109270575B (zh) 2019-11-26

Similar Documents

Publication Publication Date Title
Yamamoto et al. Stochastic model for earthquake ground motion using wavelet packets
Bogdanoff et al. Response of a simple structure to a random earthquake-type disturbance
CN108387933B (zh) 一种确定地层品质因子的方法、装置及系统
CN109375252B (zh) 考虑不同发震构造最大可信地震的地震动参数评价方法
CN109375253B (zh) 基于全部发震构造最大可信地震的地震动参数评价方法
CN110083920A (zh) 一种地震作用下非比例阻尼结构随机响应的分析方法
CN105891888A (zh) 多域分频并行多尺度全波形反演方法
CN105205495A (zh) 基于emd-elm的非平稳脉动风速预测方法
CN107886078A (zh) 一种基于分层自适应阈值函数的小波阈值降噪方法
CN109270575B (zh) 一种基于建筑物地震响应等效的爆破地震波模型构造方法
Li et al. Simulation and generation of spectrum-compatible ground motions based on wavelet packet method
Bani-Hani et al. A multi-step approach to generate response-spectrum-compatible artificial earthquake accelerograms
CN114578417A (zh) 一种基于小波调整的脉冲型地震动拟合方法
Marano Non-stationary stochastic modulation function definition based on process energy release
CN103577877B (zh) 一种基于时频分析和bp神经网络的船舶运动预报方法
Basu et al. Comparative 1D ground response analysis of homogeneous sandy stratum using Linear, Equivalent Linear and Nonlinear Masing approaches
Majidi et al. Estimation of earthquake frequency content and its effect on dynamic analysis using continuous and discrete wavelet transform
CN107506521B (zh) 一种输电铁塔三维风荷载模拟方法
Sokolov et al. Instrumental intensity distribution for the Hector Mine, California, and the Chi-Chi, Taiwan, earthquakes: Comparison of two methods
Beneldjouzi et al. Numerical and random simulation procedure for preliminary local site characterization and site factor assessing
CN108416816B (zh) 多维核磁共振测井数据的压缩处理方法和装置
Quiroz‐Ramírez et al. Evaluation of the intensity measure approach in performance‐based earthquake engineering with simulated ground motions
CN110796036A (zh) 一种用于提高结构模态参数识别精度的方法
Hao Input seismic motions for use in the structural response analysis
Hole et al. Measurement and simulation of low frequency impulse noise and ground vibration from airblasts

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
CB03 Change of inventor or designer information
CB03 Change of inventor or designer information

Inventor after: Wei Haixia

Inventor after: Chen Jinfeng

Inventor after: Zhu Jie

Inventor after: Chu Huaibao

Inventor after: Qu Zheng

Inventor before: Wei Haixia

Inventor before: Chen Jinfeng

Inventor before: Zhu Jie

Inventor before: Chu Huaibao

Inventor before: Qu Zheng

GR01 Patent grant
GR01 Patent grant