CN108169794A - 近断层人工地震波与规范反应谱匹配的调整方法 - Google Patents

近断层人工地震波与规范反应谱匹配的调整方法 Download PDF

Info

Publication number
CN108169794A
CN108169794A CN201711341166.XA CN201711341166A CN108169794A CN 108169794 A CN108169794 A CN 108169794A CN 201711341166 A CN201711341166 A CN 201711341166A CN 108169794 A CN108169794 A CN 108169794A
Authority
CN
China
Prior art keywords
response spectrum
seismic wave
function
code response
spectrum
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
CN201711341166.XA
Other languages
English (en)
Other versions
CN108169794B (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 Railway Eryuan Engineering Group Co Ltd CREEC
Original Assignee
China Railway Eryuan Engineering Group Co Ltd CREEC
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 Railway Eryuan Engineering Group Co Ltd CREEC filed Critical China Railway Eryuan Engineering Group Co Ltd CREEC
Priority to CN201711341166.XA priority Critical patent/CN108169794B/zh
Publication of CN108169794A publication Critical patent/CN108169794A/zh
Application granted granted Critical
Publication of CN108169794B publication Critical patent/CN108169794B/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
    • G01V1/36Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00Details of seismic processing or analysis
    • G01V2210/50Corrections or adjustments related to wave propagation

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

本发明提供的近断层人工地震波与规范反应谱匹配的调整方法,通过将计算反应谱展开分解成若干层后,通过优化算法对需要调整的系数进行调整和重构小波,快速有效地实现与规范反应谱之间的匹配,使计算反应谱与规范反应谱的匹配程度更高,通过缩放Volterra级数能够提高匹配的精度和计算效率,有利于对工程结构进行抗震分析,提高抗震设计和结构动力分析的精确性和有效性,优化桥梁设计能力,提高桥梁安全系数,延长使用寿命,减轻震后破坏程度,降低震后维养成本,大幅降低匹配难度,从而能够达到对工程建筑结构、桥梁结构的抗震设计要求提供指导的技术效果,优化工程结构尤其是桥梁结构设计,提高设计效率。

Description

近断层人工地震波与规范反应谱匹配的调整方法
技术领域
本发明涉及人工地震波合成技术领域,特别涉及一种近断层人工地震波与规范反应谱匹配的调整方法。
背景技术
近断层地震动是地面运动中最为复杂的一种,其具有显著的竖向加速度分量,水平分量具有很强的断层破裂方向性效应,加速度、速度和位移峰值差别较大,竖向峰值加速度与水平峰值加速度之比远大于远场地震动的比值,其速度时程表现为大幅值、长周期、短持时的强脉冲形式,会对工程结构的产生较大的破坏,因此,在对结构进行抗震分析时,考虑近断层效应十分必要。
现有技术中,由于实际的近断层地震记录数量有限,尤其是竖向地震波的记录更少,需要采用人工合成地震波进行工程结构设计抗震模拟分析,根据远场地震动的功率谱与规范反应谱的对应关系,模拟生成的远场水平地震波很容易就能与规范反应谱实现匹配,竖向地震波再根据与水平加速度与竖向加速度之间的比值进行简化计算,也能够很容易的实现匹配,但具有近断层特性的地震波就难以与规范反应谱很好的匹配,并且匹配过程十分复杂繁琐,效率低下,同时在近断层地震作用下,竖向分量的影响较大,不容忽视,也需要更为精确的计算,因此,出于对桥梁工程结构设计分析的需要,以便更好对桥梁工程结构尤其指导对桥梁结构抗震设计要求,如何使具有近断层特点的人工地震波与规范反应谱匹配十分重要。
发明内容
本发明所要解决的技术问题在于针对现有生成的近断层特性的人工地震波与规范反应谱的匹配度较差,误差较大,匹配过程复杂,效率较低等上述不足,提供一种近断层人工地震波与规范反应谱匹配的调整方法,为了实现上述发明目的,本发明提供了以下技术方案:
一种近断层人工地震波与规范反应谱匹配的调整方法,包括如下步骤:
(1)选取速度冲击脉冲波函数和远场地震功率谱函数,得到时程函数,再乘以包络函数,生成具有近断层效应的第一地震波;
(2)对生成的第一地震波进行小波分解,通过Volterra级数展开得到第一地震波的每一层的系数和对应的小波成分,建立第一地震波的时程函数与其计算反应谱之间匹配的非线性系数关系;
(3)将得到的第一地震波的计算反应谱与规范反应谱相比较,得出计算反应谱与规范反应谱之间的误差指数;
(4)通过优化算法调整对应层的系数和重构小波来调整计算反应谱,使所述误差指数小于预定的容许误差;
(5)将根据满足要求的计算反应谱求得的时程函数乘以包络函数,即可获得与规范反应谱匹配的第二地震波。
本方法为了与场地的规范反应谱函数,即规范规定函数或者场地安评所得函数相匹配,在得到了人工生成的包含有近断层和远场人工波时程函数后,乘以包络函数获得非平稳随机过程的第一地震波,然后将第一地震波采用Volterra级数展开和小波分解,建立时程函数与计算反应谱之间的非线性系数的关系,再应用优化算法调整对应层的系数和重构小波,使第一地震波的计算反应谱与规范反应谱之间的误差指数满足要求,最终得出与规范谱匹配程度较好的且符合近断层场地特征的第二地震波,采用本方法通过将计算反应谱展开分解成若干层后,通过优化算法对需要调整的系数进行调整和重构小波,快速有效地实现与规范反应谱之间的匹配,使计算反应谱与规范反应谱的匹配程度更高,通过缩放Volterra级数能够提高匹配的精度和计算效率,有利于对工程结构进行抗震分析,提高抗震设计和结构动力分析的精确性和有效性,优化桥梁设计能力,提高桥梁安全系数,延长使用寿命,减轻震后破坏程度,降低震后维养成本,同时采用此方法模拟生成竖向地震波,模拟结果更加精确,使同时具有水平加速度分量和竖向加速度分量的人工地震波与规范反应谱更好的相互匹配,大幅降低匹配难度,从而能够达到对工程建筑结构、桥梁结构的抗震设计要求提供指导的技术效果,优化工程结构尤其是桥梁结构设计,提高设计效率。
进一步的,所述误差指数根据公式:
计算,其中, 表示第k个控制点的规范反应谱值,Sak) 表示第k个控制点的计算反应谱值,k=1,2,…Nc
进一步的,在所述步骤(2)中,通过所述Volterra级数得到的第一地震波的表达式如下:
其中,hk表示第k层的系数;u表示时间函数,τ为采样时间点。
进一步的,所述优化算法包含包括蚁群算法或粒子群算法。
与现有技术相比,本发明的有益效果:本发明通过将计算反应谱通过Volterra级数展开分解成若干层和小波分解后,通过优化算法对需要调整的系数进行调整和重构小波,快速有效地实现与规范反应谱之间的匹配,使计算反应谱与规范反应谱的匹配程度更高,通过缩放 Volterra级数能够提高匹配的精度和计算效率,有利于对工程结构进行抗震分析,提高抗震设计和结构动力分析的精确性和有效性,优化桥梁设计能力,提高桥梁安全系数,延长使用寿命,减轻震后破坏程度,降低震后维养成本,同时采用此方法模拟生成竖向地震波,模拟结果更加精确,使同时具有水平加速度分量和竖向加速度分量的人工地震波与规范反应谱更好的相互匹配,大幅降低匹配难度,从而能够达到对工程建筑结构、桥梁结构的抗震设计要求提供指导的技术效果,优化工程结构尤其是桥梁结构设计,提高设计效率。
附图说明
图1为本发明中一种近断层人工地震波与规范反应谱匹配的调整方法的流程图。
具体实施方式
下面结合实施例及具体实施方式对本发明作进一步的详细描述。但不应将此理解为本发明上述主题的范围仅限于以下的实施例,凡基于本发明内容所实现的技术均属于本发明的范围。
实施例1
一种近断层人工地震波与规范反应谱匹配的调整方法,包括如下步骤:
(1)选取速度冲击脉冲波函数和远场地震功率谱函数,得到时程函数,再乘以包络函数,生成具有近断层效应的第一地震波;
(2)对生成的第一地震波进行小波分解,通过Volterra级数展开得到第一地震波的每一层的系数和对应的小波成分,建立第一地震波的时程函数与其计算反应谱之间匹配的非线性系数关系;
(3)将得到的第一地震波的计算反应谱与规范反应谱相比较,得出计算反应谱与规范反应谱之间的误差指数;
(4)通过优化算法调整对应层的系数和重构小波来调整计算反应谱,使所述误差指数小于预定的容许误差;
(5)将根据满足要求的计算反应谱求得的时程函数乘以包络函数,即可获得与规范反应谱匹配的第二地震波。
先选取速度冲击脉冲波函数和远场地震功率谱函数,得到一个时程函数,用时程函数乘以包络函数得到一个非平稳的具有近断层效应的第一地震波,对生成的第一地震波进行小波分解,根据需要的精度,确定Volterra级数展开的级数,展开得到第一地震波的每一层的系数和对应的小波成分,得到的第一地震波的表达式如下:
其中,hk表示第k层的系数;u表示时间函数,τ为采样时间点,如0.02s、0.04s…。
然后建立第一地震波的时程函数与规范反应谱之间匹配的非线性系数关系,将得到的第一地震波的计算反应谱与规范反应谱相比较,得出计算反应谱与规范反应谱之间的误差指数,所述误差指数根据公式:
计算,其中, 表示第k个控制点的规范反应谱值,Sak) 表示第k个控制点的计算反应谱值,k=1,2,…Nc;ε为预设的容许误差,当误差指数εm大于ε时,通过优化算法,如蚁群算法或粒子群算法,调整需要调整的对应层的系数和重构小波,即缩放Volterra级数,直至误差指数εm小于或等于ε,根据满足误差的第一地震波的计算反应谱算出对应的时程函数,将此时程函数乘以包络函数,得到第二地震波,得到的所述第二地震波具有近断层非平稳特性,并与规范反应谱匹配程度更高,计算效率更快。

Claims (4)

1.一种近断层人工地震波与规范反应谱匹配的调整方法,其特征在于,包括如下步骤:
(1)选取速度冲击脉冲波函数和远场地震功率谱函数,得到时程函数,再乘以包络函数,生成具有近断层效应的第一地震波;
(2)对生成的第一地震波进行小波分解,通过Volterra级数展开得到第一地震波的每一层的系数和对应的小波成分,建立第一地震波的时程函数与其计算反应谱之间匹配的非线性系数关系;
(3)将得到的第一地震波的计算反应谱与规范反应谱相比较,得出计算反应谱与规范反应谱之间的误差指数;
(4)通过优化算法调整对应层的系数和重构小波来调整计算反应谱,使所述误差指数小于预定的容许误差;
(5)将根据满足要求的计算反应谱求得的时程函数乘以包络函数,即可获得与规范反应谱匹配的第二地震波。
2.如权利要求1所述的近断层人工地震波与规范反应谱匹配的调整方法,其特征在于,所述误差指数根据公式:
计算,其中, 表示第k个控制点的规范反应谱值,Sak)表示第k个控制点的计算反应谱值,k=1,2,…Nc
3.如权利要求1或2所述的近断层人工地震波与规范反应谱匹配的调整方法,其特征在于,在所述步骤(2)中,通过所述Volterra级数展开得到的第一地震波的表达式如下:
其中,hk表示第k层的系数;u表示时间函数,τ为采样时间点。
4.如权利要求1或2所述的近断层人工地震波与规范反应谱匹配的调整方法,其特征在于,所述优化算法包括蚁群算法或粒子群算法。
CN201711341166.XA 2017-12-14 2017-12-14 近断层人工地震波与规范反应谱匹配的调整方法 Active CN108169794B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201711341166.XA CN108169794B (zh) 2017-12-14 2017-12-14 近断层人工地震波与规范反应谱匹配的调整方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201711341166.XA CN108169794B (zh) 2017-12-14 2017-12-14 近断层人工地震波与规范反应谱匹配的调整方法

Publications (2)

Publication Number Publication Date
CN108169794A true CN108169794A (zh) 2018-06-15
CN108169794B CN108169794B (zh) 2020-06-30

Family

ID=62525447

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201711341166.XA Active CN108169794B (zh) 2017-12-14 2017-12-14 近断层人工地震波与规范反应谱匹配的调整方法

Country Status (1)

Country Link
CN (1) CN108169794B (zh)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110020400A (zh) * 2019-01-07 2019-07-16 河海大学 一种调整地震波与目标反应谱精确匹配的影响矩阵方法
CN110069836A (zh) * 2019-04-03 2019-07-30 河海大学 一种高低频段交替与目标谱匹配的改进影响矩阵方法
CN110244355A (zh) * 2019-07-25 2019-09-17 西南交通大学 一种基于震源断层模型的脉冲地震动模拟方法
CN112285770A (zh) * 2020-10-28 2021-01-29 重庆大学 一种基于实数小波变换的地震波相位谱扰动方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2002168964A (ja) * 2000-11-30 2002-06-14 Tokyo Electric Power Co Inc:The 地震動の強さ推定方法、その装置及びそのプログラムを記録した記録媒体
CN101236256A (zh) * 2008-01-25 2008-08-06 北京工业大学 多阻尼比目标反应谱兼容的人工地震波合成方法
CN103837334A (zh) * 2012-11-20 2014-06-04 中国石油大学(北京) 一种机械故障的检测方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2002168964A (ja) * 2000-11-30 2002-06-14 Tokyo Electric Power Co Inc:The 地震動の強さ推定方法、その装置及びそのプログラムを記録した記録媒体
CN101236256A (zh) * 2008-01-25 2008-08-06 北京工业大学 多阻尼比目标反应谱兼容的人工地震波合成方法
CN103837334A (zh) * 2012-11-20 2014-06-04 中国石油大学(北京) 一种机械故障的检测方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
N.A.ALEXANDER 等: "Obtaining Spectrum Matching Time Series Using a Reweighted Volterra Series Algorithm(RVSA)", 《BULLETIN OF THE SEISMOLOGICAL SOCIETY OF AMERICA》 *

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110020400A (zh) * 2019-01-07 2019-07-16 河海大学 一种调整地震波与目标反应谱精确匹配的影响矩阵方法
CN110020400B (zh) * 2019-01-07 2023-05-23 河海大学 一种调整地震波与目标反应谱精确匹配的影响矩阵方法
CN110069836A (zh) * 2019-04-03 2019-07-30 河海大学 一种高低频段交替与目标谱匹配的改进影响矩阵方法
CN110244355A (zh) * 2019-07-25 2019-09-17 西南交通大学 一种基于震源断层模型的脉冲地震动模拟方法
CN110244355B (zh) * 2019-07-25 2021-06-08 西南交通大学 一种基于震源断层模型的脉冲地震动模拟方法
CN112285770A (zh) * 2020-10-28 2021-01-29 重庆大学 一种基于实数小波变换的地震波相位谱扰动方法

Also Published As

Publication number Publication date
CN108169794B (zh) 2020-06-30

Similar Documents

Publication Publication Date Title
CN108169794A (zh) 近断层人工地震波与规范反应谱匹配的调整方法
Feng et al. Data‐driven modeling of the solar wind from 1 Rs to 1 AU
CN108073075A (zh) 基于ga优化bp神经网络的硅微加速度计温度补偿方法、系统
Spanos et al. An approach for non-linear stochastic analysis of U-shaped OWC wave energy converters
CN109991387B (zh) 一种模拟γ-TiAl合金在非等温条件下的相场分析法
CN107967240A (zh) 基于人工地震动生成过程的耐震时程加速度优化算法
CN110794459A (zh) 一种海底近断层地震动的拟合方法
CN106546400B (zh) 一种多输入多输出非高斯随机振动试验系统及试验方法
CN107894618B (zh) 一种基于模型平滑算法的全波形反演梯度预处理方法
CN113806845B (zh) 一种基于pod插值的非平稳风场高效模拟方法
CN109002678B (zh) 基于水电工程水工建筑物抗震设计的人工地震模拟方法
CN114578417B (zh) 一种基于小波调整的脉冲型地震动拟合方法
CN111651896B (zh) 一种基于实际风速及实际复杂地形的风电场流场计算方法
CN107422376B (zh) 一种非平稳非高斯地震动时程模拟方法
CN104392717A (zh) 一种基于声道谱高斯混合建模的快速语音转换系统及其方法
CN110244355B (zh) 一种基于震源断层模型的脉冲地震动模拟方法
CN114442153A (zh) 一种近断层地震动拟合方法
CN109446552B (zh) 多轴相关随机激励下结构疲劳寿命时域计算方法
Vasylkevych et al. A high‐accuracy global prognostic model for the simulation of Rossby and gravity wave dynamics
BENOIT et al. A nearshore wave atlas along the coasts of France based on the numerical modeling of wave climate over 25 years
CN110020400A (zh) 一种调整地震波与目标反应谱精确匹配的影响矩阵方法
Kaurkin et al. Assimilation of the AVISO altimetry data into the ocean dynamics model with a high spatial resolution using Ensemble Optimal Interpolation (EnOI)
Zalesny et al. Numerical model of nonhydrostatic ocean dynamics based on methods of artificial compressibility and multicomponent splitting
Pogarskii et al. Joint analysis of the wind and wave-field variability in the Indian Ocean area for 1998–2009
WO2022051030A3 (en) Generative modeling of quantum hardware

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