CN114966835B - 一种基于改进时域叠加法的近断层地震动拟合方法 - Google Patents

一种基于改进时域叠加法的近断层地震动拟合方法 Download PDF

Info

Publication number
CN114966835B
CN114966835B CN202210433801.1A CN202210433801A CN114966835B CN 114966835 B CN114966835 B CN 114966835B CN 202210433801 A CN202210433801 A CN 202210433801A CN 114966835 B CN114966835 B CN 114966835B
Authority
CN
China
Prior art keywords
earthquake motion
spectrum
fault
time
frequency
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
CN202210433801.1A
Other languages
English (en)
Other versions
CN114966835A (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.)
Fuzhou University
Original Assignee
Fuzhou University
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 Fuzhou University filed Critical Fuzhou University
Priority to CN202210433801.1A priority Critical patent/CN114966835B/zh
Publication of CN114966835A publication Critical patent/CN114966835A/zh
Application granted granted Critical
Publication of CN114966835B publication Critical patent/CN114966835B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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/28—Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/282—Application of seismic models, synthetic seismograms
    • 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. for interpretation or for event detection
    • 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. for interpretation or for event detection
    • G01V1/36—Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
    • G01V1/364—Seismic filtering
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G06F17/14—Fourier, Walsh or analogous domain transformations, e.g. Laplace, Hilbert, Karhunen-Loeve, transforms
    • G06F17/141—Discrete Fourier transforms
    • G06F17/142—Fast Fourier transforms, e.g. using a Cooley-Tukey type algorithm

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Remote Sensing (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Mathematical Physics (AREA)
  • Environmental & Geological Engineering (AREA)
  • Geology (AREA)
  • Geophysics (AREA)
  • Acoustics & Sound (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Data Mining & Analysis (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Discrete Mathematics (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Algebra (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

本发明涉及一种基于改进时域叠加法的近断层地震动拟合方法。包括：确定目标加速度反应谱；求解功率谱；由功率谱计算傅里叶幅值谱并与相位谱结合生成加速度时程；对生成加速度时程判断是否非平稳进而生成初始无脉冲地震动；对初始无脉冲地震动进行低频滤波生成近断层地震动时程的高频成分；使用等效速度脉冲模型生成近断层地震动时程的低频成分；对近断层地震动时程的高频成分和低频成分在时域叠加过程中的高低频峰值时刻提出了明确定义，生成新的近断层地震动低频成分并与近断层地震动高频成分在时域上进行叠加生成原始地震动反应谱；对原始地震动反应谱采用FFT法进行调整反应谱生成更符合实际情况近断层地震动。

Description

一种基于改进时域叠加法的近断层地震动拟合方法
技术领域
本发明涉及一种基于改进时域叠加法的近断层地震动拟合方法。
背景技术
随着国民经济和科技水平的不断发展，近年来越来越多大尺寸的建(构)筑作为生命线被修建在近断层区域甚至跨越断层。这种建(构)筑的自振周期往往都比一般的建(构)筑物大，容易受到低频率激励的影响。而近断层地震动是一种含有丰富低频成分的地震动，因此，考虑近断层地震动对这种结构的影响是很有必要的。到目前为止，近断层脉冲型地震动的地震记录还很匮乏，这对于缺少这种地震记录的区域来说，在进行大型建筑物的抗震设计时，是没有实测地震记录可以作为地震动输入的。但人造地震动可以在短时间内产生具有满足拟建场地条件的一系列地震动，这能有效弥补这个不足。所以，有必要研究拟合近断层脉冲型地震动的方法。人工合成地震动不仅可以得到满足各种场地条件的地震波，而且在一定程度上也促进了结构抗震设计的发展。因此，研究能够合成尽可能符合实际情况地震动的方法，显得尤为重要。
现有技术1【1】近断层区域划分及近断层速度脉冲型地震动模拟，其详细步骤有以下几步：
步骤S1：确定目标加速度反应谱Sa T(ξ,ω)：
目标加速度反应谱根据《公路工程抗震规范》(JTGB02-2013)中的5.2.1条确定，其形式如下式所示：
式中：
Tg—反应谱的特征周期；
T—结构自振周期；
Sa T(ξ,ω)—目标加速度反应谱；
Smax—反应谱最大值。
步骤S2：计算目标功率谱S(ω)；
根据当前采用提出的反应谱与功率谱转化的关系，其形式如下式所示：
式中：
Sa T(ξ,ω)—目标加速度反应谱；
S(ω)—目标功率谱；
ξ—阻尼比；
ω—结构自振频率；
T—结构自振周期；
γ—超越概率。
步骤S3：计算傅利叶幅值谱A(ω)；
功率谱和傅里叶幅值谱的数学关系如下式所示：
式中：A(ω)—傅里叶幅值谱
步骤S4：计算相位谱；
生成在-2π到2π上均匀分布的相位谱；
步骤S5：由IFFT得到加速度时程xa(t)；
采用下式所示的三角级数模型来模拟地震动X。
X＝g(t)·xa(t)
式中：
—相位角；
ωk—圆频率；
A(ωk)—幅值谱；
N—目标反应谱的的频域中频率的分隔点数；
g(t)—强度包络函数；
xa(t)—加速度时程；
X—地震动
利用快速傅里叶逆变换，计算得到平稳加速度时程；
步骤S6：乘以包络函数生成非平稳加速度时程；
为了将平稳加速度时程转化成非平稳加速度时程，通过乘以强度包络函数对生成的平稳加速度时程进行非均匀调制。包络函数的数学表达式如式所示
其中：
T1—平稳段的起始时刻；
T2—平稳段的结束时刻；
c—衰减阶段的变化速率。
步骤S7：与目标加速度反应谱拟合迭代，得到加速度时程；
为了让非平稳加速度时程的反应谱匹配目标加速度反应谱，需要对非平稳加速度时程进行迭代计算，使生成的非平稳加速度时程的反应谱逼近目标谱。具体的迭代计算步骤为：
a.计算出目标加速度反应谱与非平稳加速度时程的反应谱的比值；
b.对非平稳加速度时程作快速傅里叶变换计算得到相位谱和傅里叶幅值谱,将a中计算出的比值乘以傅里叶幅值谱；
c.将调整后的傅里叶幅值谱与相位谱结合，通过快速傅里叶逆变换计算得到平稳加速度时程；
d.乘以强度包络函数计算得到非平稳加速度时程；
e.计算非平稳加速度时程的反应谱与目标加速度反应谱之间的平均相对误差Em并判断其是否小于5％，若小于5％，则停止迭代计算，否则回到a。
经过n次迭代计算后，非平稳加速度时程的反应谱与目标加速度反应谱之间的平均相对误差为小于5％，即停止迭代计算。最终得到计算加速度反应谱曲线和目标加速度反应谱曲线.
步骤S8：将步骤S7最终生成的加速度时程经傅里叶变换得到傅里叶谱，将频率小于fr的幅值设为0
步骤S9：再将步骤S8最终生成的傅里叶幅值谱与步骤S4生成的相位谱结合，通过快速傅里叶逆变换计算生成高频加速度时程
步骤S10：采用等效速度脉冲模型来模拟低频速度脉冲时程：
式中ω(t)为包络函数，用式计算。
其中：
Vp—脉冲峰值；
fp—频率，取1/Tp，Tp为脉冲周期；
t1—余弦函数峰值发生时刻；
T—结构自振周期；
γ—衰减速率；
t0—包络函数峰值发生时刻；
将上述低频速度脉冲时程求导，得到低频加速度时程；
步骤S11：高，低频成分叠加生成近断层脉冲型地震动；
将步骤S9得到的高频加速度时程和步骤S10得到的低频脉冲加速度时程叠加生成近断层脉冲型地震动。
现有技术1【1】基于时域叠加法得到的近断层脉冲型地震动存在拟合得到的近断层脉冲型地震动的反应谱与目标反应谱不吻合的缺点。这是因为现有方法中，先进行高频成分加速度时程与目标反应谱的调整(步骤S7-S9)，然后进行高频成分与低频成分脉冲加速度时程在时域上的叠加(步骤S11)，得到近断层脉冲型地震动。由于高频成分与低频成分叠加后，引入了新的低频成分后，原来的拟合得到的反应谱又发生了改变，因而会导致拟合得到地震动与目标反应谱有很有明显的出入。因此，现有技术的不足会导致拟合的近断层脉冲型地震动的反应谱与目标反应谱有差别的明显不足。
发明内容
本发明的目的在于提供一种基于改进时域叠加法的近断层地震动拟合方法，解决基于时域叠加法生成近断层地震动反应谱与目标反应谱不吻合的难题。目前基于时域叠加法生成近断层脉冲型地震动是先调整高频成分加速度时程的反应谱，实现了高频成分加速度时程的反应谱与目标反应谱的吻合，再与低频成分脉冲加速度时程在时域上进行叠加，最终生成的近断层脉冲型地震动的反应谱在整个周期段与目标反应谱相差较大。所以，虽然基于时域叠加法生成近断层地震动在时域上具有速度脉冲特性，但是其生成的近断层地震动反应谱与目标反应谱不吻合。
为实现上述目的，本发明的技术方案是：一种基于改进时域叠加法的近断层地震动拟合方法，包括：确定目标加速度反应谱；求解目标功率谱；由功率谱计算傅里叶幅值谱并与相位谱结合生成加速度时程；对生成加速度时程判断是否非平稳加速度时程进而生成初始无脉冲地震动；对初始无脉冲地震动进行低频滤波生成近断层地震动时程的高频成分；使用等效速度脉冲模型生成近断层地震动时程的低频成分；对近断层地震动时程的高频成分和低频成分在时域叠加过程中的高低频峰值时刻提出了明确定义，生成新的近断层地震动低频成分并与近断层地震动高频成分在时域上进行叠加生成原始地震动反应谱；对原始地震动反应谱采用FFT法进行调整反应谱生成更符合实际情况近断层地震动。
在本发明一实施例中，所述确定目标加速度反应谱的具体实现步骤如下：
1)、根据式(1)求出速度放大系数设计谱βV(T)
其中βmm为各类场地上的拟速度均值谱峰值，根据式(2)确定：
βmm＝ΩCs (2)
式中：βV(T)为速度放大系数设计谱；Tg为反应谱的特征周期；Tp为脉冲周期；T为结构自振周期；Ω为各脉冲地震动记录速度放大系数谱最大值的平均值；Cs为场地系数；βmm为各类场地上的拟速度均值谱峰值；
2)、根据式(3)得等效加速度放大系数谱βVa(T)
式中：βVa(T)为等效加速度放大系数谱；ω为结构自振频率；PGV/PGA为峰值地面速度与峰值地面加速度比值；
3)、根据式(4)求得目标加速度反应谱Sa T(ξ,ω)，如下所示：
Sa T(ξ，ω)＝CRCdAβVa(T) (4)
式中：Sa T(ξ,ω)为目标加速度反应谱；βVa(T)为等效加速度放大系数谱；CR为风险系数；Cd为阻尼调整系数；A为设计基本地震动加速度峰值。
在本发明一实施例中，所述目标功率谱的计算公式如下：
式中：Sa T(ξ,ω)为目标加速度反应谱；S(ω)为目标功率谱；ξ为阻尼比；ω为结构自振频率；T为结构自振周期；λ为超越概率。
在本发明一实施例中，所述由功率谱计算傅里叶幅值谱并与相位谱结合生成加速度时程的具体实现步骤如下：
1)根据式(6)得到傅里叶幅值谱A(ω)：
Δω＝2π×fr×FFT长度
式中：A(ω)为傅里叶幅值谱；Δω为频率间隔；S(ω)为目标功率谱；
2)计算相位谱；生成在-2π到2π上均匀分布的相位谱；
3)根据式(7)生成加速度时程xa(t)：
式中：φk为相位角；ωk为圆频率；A(ωk)为幅值谱；N为目标反应谱的的频域中频率的分隔点数；g(t)为强度包络函数；xa(t)为加速度时程。
在本发明一实施例中，所述对生成加速度时程判断是否非平稳加速度时程进而生成初始无脉冲地震动的具体实现步骤如下：
判断生成的加速度时程是否为非平稳加速度时程，若是，即直接生成非平稳加速时程，即初始无脉冲地震动Xa；若为平稳加速度时程，则根据式(8)生成初始无脉冲地震动Xa：
Xa＝xa(t)×g(t) (8)
式中：Xa为初始无脉冲地震动；xa(t)为加速度时程；g(t)为强度包络函数；t为时间；T1为平稳段的起始时刻；T2为平稳段的结束时刻；c为衰减阶段的变化速率。
在本发明一实施例中，所述对初始无脉冲地震动进行低频滤波生成近断层地震动时程的高频成分的具体实现方式为：对生成的初始无脉冲地震动做傅里叶变换，将频率区间[0,fr]所对应的傅里叶幅值谱的值设为零，再经过傅里叶逆变换得到近断层地震动时程的高频成分Ah(t)。
在本发明一实施例中，所述使用等效速度脉冲模型生成近断层地震动时程的低频成分的具体实现步骤如下：
1)根据式(10)生成低频脉冲速度时程ν(t)
式中：ν(t)为低频脉冲速度时程；ω(t)为低频脉冲速度时程包络函数；Vp为脉冲峰值；t为时间；fp为脉冲频率；取1/Tp，Tp为脉冲周期；t1为余弦函数峰值发生时刻；T为结构自振周期；γ为衰减速率；t0为包络函数峰值发生时刻；
2)对生成的低频脉冲速度时程进行一次求导得到低频成分脉冲加速度时程α(t)；
3)得到近断层地震动时程的低频成分Al(t)：将得到的低频成分脉冲加速度时程在频率区间为[fr,∞)所对应的傅里叶幅值谱的值设为零，再经过傅里叶逆变换得到近断层地震动时程的低频成分Al(t)。
在本发明一实施例中，所述对近断层地震动时程的高频成分和低频成分在时域叠加过程中的高低频峰值时刻提出了明确定义，生成新的近断层地震动低频成分并与近断层地震动高频成分在时域上进行叠加生成原始地震动反应谱的具体实现步骤如下：
1)计算近断层地震动时程的高频成分Ah(t)峰值时刻tha，如式(12)所示；
ln(tha)＝1.35MW-6.88 (12)
其中：Mw为矩震级；tha为近断层地震动时程的高频成分Ah(t)峰值时刻；
2)计算近断层地震动时程的低频成分Al(t)峰值时刻tla，如式(13)所示；
ln(tla)＝1.35MW-6.88 (13)
其中：Mw为矩震级；tla为近断层地震动时程的低频成分Al(t)峰值时刻；
3)由高频峰值时刻tha和低频峰值时刻tla求得高低频峰值时刻差σt，如式(14)所示；
σt＝tha-tla (14)
4)把近断层地震动时程的低频成分Al(t)在时间轴上平移σt，得到新的近断层地震动时程的低频成分All(t)；
5)将新的近断层地震动时程的低频成分All(t)与近断层地震动时程的高频成分Ah(t)在时域进行叠加，得到原始地震动的反应谱Sa。
在本发明一实施例中，所述对原始地震动反应谱采用FFT法进行调整反应谱生成更符合实际情况近断层地震动的具体实现步骤如下：
1)求原始地震动的反应谱和目标加速度反应谱的平均相对误差；
2)判断得到的平均相对误差是否小于固定阈值；若平均相对误差大于固定阈值，则计算原始地震动的反应谱和目标加速度反应谱的比值，调整傅里叶幅值谱，进入步骤3)；若平均相对误差小于等于固定阈值，则输出近断层脉冲型地震动A(t)；
3)计算相位谱：重新生成在-2π到2π上均匀分布的相位谱；
4)根据式(15)生成新的加速度时程xp(t)
式中：φk为相位角；ωk为圆频率；A(ωk)为幅值谱；N为目标反应谱的的频域中频率的分隔点数；g(t)为强度包络函数；xp(t)为新的加速度时程；
5)判断生成的加速度时程是否为非平稳加速度时程，若是，即直接生成非平稳加速时程，即新的地震动Xp；若为平稳加速度时程，根据式(16)生成新的地震动Xp
Xp＝xp(t)×g(t) (16)
式中：Xp为新的地震动；g(t)为地震动包络函数；t为时间；T1为平稳段的起始时刻；T2为平稳段的结束时刻；c为衰减阶段的变化速率；
6)根据生成的新的地震动计算地震动反应谱；求地震动的反应谱和目标加速度反应谱的平均相对误差，并判断其是否小于固定阈值；若平均相对误差大于固定阈值，则重复步骤3)～步骤5)；若平均相对误差小于等于固定阈值，则输出近断层脉冲型地震动A(t)。
相较于现有技术，本发明具有以下有益效果：本发明提出的一种基于改进时域叠加法生成近断层地震动反应谱的方法，该方法生成的近断层地震动反应谱与目标反应谱相吻合。本发明先进行近断层脉冲型地震动高频成分拟合(步骤S1-步骤S8)，再进行近断层脉冲型地震动低频成分拟合(步骤S9-步骤S11)，最后利用FFT法对近断层脉冲型地震动的高频成分和低频成分在时域上叠加生成的地震动加速度时程的反应谱进行调整(步骤S16-步骤S21)，得到近断脉冲型地震动。通过对调整后得到的近断层脉冲型地震动的反应谱与目标反应谱进行对比，可以看出经过循环调整后的近断层脉冲型地震动的计算反应谱与目标反应谱在中长周期段基本完全重合，在短周期段误差较小。所以基于改进时域叠加法生成的近断层地震动拟合得到的近断层地震动反应谱不仅都保留了高频成分和低频成分的特性，而且都能够很好的拟合目标反应谱。
附图说明
图1为目标近断层地震动目标反应谱曲线。
图2为现有技术拟合的近断层地震动。
图3为本发明方法拟合的近断层地震动。
图4为本发明拟合近断层地震动的反应谱。
图5为本发明一实施例方法流程图。
具体实施方式
下面结合附图，对本发明的技术方案进行具体说明。
本发明提出了一种基于改进时域叠加法的近断层地震动拟合方法，包括：确定目标加速度反应谱；求解目标功率谱；由功率谱计算傅里叶幅值谱并与相位谱结合生成加速度时程；对生成加速度时程判断是否非平稳加速度时程进而生成初始无脉冲地震动；对初始无脉冲地震动进行低频滤波生成近断层地震动时程的高频成分；使用等效速度脉冲模型生成近断层地震动时程的低频成分；对近断层地震动时程的高频成分和低频成分在时域叠加过程中的高低频峰值时刻提出了明确定义，生成新的近断层地震动低频成分并与近断层地震动高频成分在时域上进行叠加生成原始地震动反应谱；对原始地震动反应谱采用FFT法进行调整反应谱生成更符合实际情况近断层地震动。
本发明拟合近断层脉冲型地震动时程，包括三个模块，即模块一：高频成分拟合；模块二：低频成分拟合；模块三：高低频叠加后基于FFT法进行目标反应谱调整。其中，如图5所示，具体过程步骤如下：
模块一：高频成分拟合；
步骤S1：确定目标加速度反应谱Sa T(ξ,ω)；
目标加速度反应谱采用现有技术2【2】中的设计反应谱，其形式如式(1)所示：
步骤S1.1：根据式(1)求出速度放大系数设计谱βV(T)
其中βmm为各类场地上的拟速度均值谱峰值，根据式(2)确定：
βmm＝ΩCs (2)
式中：
βV(T)—速度放大系数设计谱；
Tg—反应谱的特征周期；
Tp—脉冲周期；
T—结构自振周期；
Ω—各脉冲地震动记录速度放大系数谱最大值的平均值；
Cs—场地系数；
βmm—各类场地上的拟速度均值谱峰值；
步骤S1.2：根据式(3)得等效加速度放大系数谱βVa(T)
式中：
βVa(T)—等效加速度放大系数谱；
ω—结构自振频率；
PGV/PGA—峰值地面速度与峰值地面加速度比值；
步骤S1.3：根据式(4)求得目标加速度反应谱Sa T(ξ,ω)，如下所示：
Sa T(ξ，ω)＝CRCdAβVa(T) (4)
式中：
Sa T(ξ,ω)—目标加速度反应谱；
βVa(T)—等效加速度放大系数谱；
CR—风险系数；
Cd—阻尼调整系数；
A—设计基本地震动加速度峰值；
步骤S2：根据式(5)得到目标功率谱S(ω)：
式中：
Sa T(ξ,ω)—目标加速度反应谱；
S(ω)—目标功率谱；
ξ—阻尼比；
ω—结构自振频率；
T—结构自振周期；
λ—超越概率。
步骤S3：根据式(6)得到傅里叶幅值谱A(ω)：
Δω＝2π×fr×FFT长度
式中：
A(ω)—傅里叶幅值谱；
Δω—频率间隔；
步骤S4：计算相位谱；生成在-2π到2π上均匀分布的相位谱；
步骤S5：根据式(7)生成加速度时程xa(t)：
式中：
φk—相位角；
ωk—圆频率；
A(ωk)—幅值谱；
N—目标反应谱的的频域中频率的分隔点数；
g(t)—强度包络函数；
xa(t)—加速度时程；
步骤S6：判断步骤S5生成的加速度时程是否为非平稳加速度时程，若是，即直接生成非平稳加速时程，即初始无脉冲地震动Xa；若为平稳加速度时程，则进入步骤S7；
步骤S7：根据式(8)生成初始无脉冲地震动Xa
Xa＝xa(t)×g(t) (8)
式中：
Xa—初始无脉冲地震动；
g(t)—地震动包络函数；
t—时间；
T1—平稳段的起始时刻；
T2—平稳段的结束时刻；
c—衰减阶段的变化速率；
步骤S8：得到近断层地震动时程的高频成分Ah(t)：对步骤S7生成的初始无脉冲地震动做傅里叶变换，将频率区间[0,fr]所对应的傅里叶幅值谱的值设为零，再经过傅里叶逆变换得到近断层地震动时程的高频成分Ah(t)。
模块二：低频成分拟合；
步骤S9：根据式(10)生成低频脉冲速度时程ν(t)
式中：
ν(t)—低频脉冲速度时程；
ω(t)—低频脉冲速度时程包络函数；
Vp—脉冲峰值；
t—时间；
fp—脉冲频率；取1/Tp，Tp为脉冲周期；
t1—余弦函数峰值发生时刻；
T—结构自振周期；
γ—衰减速率；
t0—包络函数峰值发生时刻；
步骤S10：对步骤S9生成的低频脉冲速度时程进行一次求导得到低频成分脉冲加速度时程α(t)
步骤S11：得到近断层地震动时程的低频成分Al(t)：将步骤S9得到的低频成分脉冲加速度时程在频率区间为[fr,∞)所对应的傅里叶幅值谱的值设为零，再经过傅里叶逆变换得到近断层地震动时程的低频成分Al(t)。
模块三：FFT法调整近断层脉冲型地震动加速度时程的反应谱
步骤S12：计算由步骤S8得到的近断层地震动时程的高频成分Ah(t)峰值时刻tha，如式(12)所示；
ln(tha)＝1.35MW-6.88 (12)
其中：
Mw—矩震级；
tha—近断层地震动时程的高频成分Ah(t)峰值时刻；
步骤S13：计算由步骤S11得到近断层地震动时程的低频成分Al(t)峰值时刻tla，如式(13)所示；
ln(tla)＝1.35MW-6.88 (13)
其中：
Mw—矩震级；
tla—近断层地震动时程的低频成分Al(t)峰值时刻；
步骤S14：由步骤S12得到的高频峰值时刻tha和步骤S13得到的峰值时刻tla求得高低频峰值时刻差σt，如式(14)所示；
σt＝tha-tla (14)
步骤S15：把近断层地震动时程的低频成分Al(t)在时间轴上平移σt，得到新的近断层地震动时程的低频成分All(t)；
步骤S16：将由步骤S15得到新的近断层地震动时程的低频成分All(t)与步骤S8得到近断层地震动时程的高频成分Ah(t)在时域进行叠加，得到原始地震动的反应谱Sa。
步骤S17：求原始地震动的反应谱和目标加速度反应谱的平均相对误差。
步骤S18：判断步骤S17得到的平均相对误差是否小于固定阈值。若平均相对误差大于固定阈值，则计算原始地震动的反应谱和目标加速度反应谱的比值，调整傅里叶幅值谱，进入步骤S19；若平均相对误差小于等于固定阈值，则输出近断层脉冲型地震动A(t)。
步骤S19：计算相位谱：重新生成在-2π到2π上均匀分布的相位谱；
步骤S20：根据式(15)生成新的加速度时程xp(t)
式中：
φk—相位角；
ωk—圆频率；
A(ωk)—幅值谱；
N—目标反应谱的的频域中频率的分隔点数；
g(t)—强度包络函数；
xp(t)—新的加速度时程；
步骤S21：判断步骤S20生成的加速度时程是否为非平稳加速度时程，若是，即直接生成非平稳加速时程，即新的地震动Xp；若为平稳加速度时程，根据式(16)生成新的地震动Xp
Xp＝xp(t)×g(t) (16)
式中：
Xp—新的地震动；
g(t)—地震动包络函数；
t—时间；
T1—平稳段的起始时刻；
T2—平稳段的结束时刻；
c—衰减阶段的变化速率。
步骤S22：根据步骤S21生成的新的地震动计算地震动反应谱；求地震动的反应谱和目标加速度反应谱的平均相对误差，并判断其是否小于固定阈值。若平均相对误差大于固定阈值，则重复步骤S19～步骤S21；若平均相对误差小于等于固定阈值，则输出近断层脉冲型地震动A(t)。
本发明实施范例
1、地震动拟合目标：
采样频率fs＝100Hz，地震动总持时T＝20s。工程场地临近断层的潜在矩震级为Mw＝6.5级；工程场地震中距为R＝8km；场地类别为Ⅱ类；抗震设防烈度为7度，50年超越概率10％对应的设计基本地震动加速度峰值为0.15g。场地系数Cs＝1.0，场地特征周期Tg＝0.9s，速度脉冲周期Tp＝2.09s，风险系数CR＝1，阻尼调整系数Cd＝1。目标近断层地震动目标反应谱曲线如图1所示。
2、现有技术的拟合结果文献[王宇航，2015]
(1)经现有技术文献[王宇航，2015]【1】的步骤S1～S11，得到原有方法结果如图2(a)～2(d)所示；
3、本申请技术的拟合结果
(2)经本发明方法得到的结果如图3(a)～3(d)所示；
4、本申请技术与现有技术的效果比较
图4为本发明方法与现有技术的效果比较图。
参考文献：
【1】王宇航.近断层区域划分及近断层速度脉冲型地震动模拟[D].西南交通大学,2015.
【2】杨华平,钱永久,黎璟,等.近断层脉冲型地震设计谱研究[J].中国公路学报,2017,30(012):159-168.。
以上是本发明的较佳实施例，凡依本发明技术方案所作的改变，所产生的功能作用未超出本发明技术方案的范围时，均属于本发明的保护范围。

Claims (6)

1.一种基于改进时域叠加法的近断层地震动拟合方法，其特征在于，包括：确定目标加速度反应谱；求解目标功率谱；由功率谱计算傅里叶幅值谱并与相位谱结合生成加速度时程；对生成加速度时程判断是否非平稳加速度时程进而生成初始无脉冲地震动；对初始无脉冲地震动进行低频滤波生成近断层地震动时程的高频成分；使用等效速度脉冲模型生成近断层地震动时程的低频成分；对近断层地震动时程的高频成分和低频成分在时域叠加过程中的高低频峰值时刻提出了明确定义，生成新的近断层地震动低频成分并与近断层地震动高频成分在时域上进行叠加生成原始地震动反应谱；对原始地震动反应谱采用FFT法进行调整反应谱生成更符合实际情况近断层地震动；
所述确定目标加速度反应谱的具体实现步骤如下：
步骤S1.1、根据式(1)求出速度放大系数设计谱βV(T)
其中βmm为各类场地上的拟速度均值谱峰值，根据式(2)确定：
βmm＝ΩCs (2)
式中：βV(T)为速度放大系数设计谱；Tg为反应谱的特征周期；Tp为脉冲周期；T为结构自振周期；Ω为各脉冲地震动记录速度放大系数谱最大值的平均值；Cs为场地系数；βmm为各类场地上的拟速度均值谱峰值；
步骤S1.2、根据式(3)得等效加速度放大系数谱βVa(T)
式中：βVa(T)为等效加速度放大系数谱；ω为结构自振频率；PGV/PGA为峰值地面速度与峰值地面加速度比值；
步骤S1.3、根据式(4)求得目标加速度反应谱Sa T(ξ,ω)，如下所示：
Sa T(ξ，ω)＝CRCdAβVa(T) (4)
式中：Sa T(ξ,ω)为目标加速度反应谱；βVa(T)为等效加速度放大系数谱；CR为风险系数；Cd为阻尼调整系数；A为设计基本地震动加速度峰值；
所述目标功率谱的计算公式如下：
式中：Sa T(ξ,ω)为目标加速度反应谱；S(ω)为目标功率谱；ξ为阻尼比；ω为结构自振频率；T为结构自振周期；λ为超越概率；
所述由功率谱计算傅里叶幅值谱并与相位谱结合生成加速度时程的具体实现步骤如下：
1)根据式(6)得到傅里叶幅值谱A(ω)：
Δω＝2π×fr×FFT长度
式中：A(ω)为傅里叶幅值谱；Δω为频率间隔；S(ω)为目标功率谱；
2)计算相位谱；生成在-2π到2π上均匀分布的相位谱；
3)根据式(7)生成加速度时程xa(t)：
式中：φk为相位角；ωk为圆频率；A(ωk)为幅值谱；N为目标反应谱的频域中频率的分隔点数；g(t)为强度包络函数；xa(t)为加速度时程。
2.根据权利要求1所述的一种基于改进时域叠加法的近断层地震动拟合方法，其特征在于，所述对生成加速度时程判断是否非平稳加速度时程进而生成初始无脉冲地震动的具体实现步骤如下：
判断生成的加速度时程是否为非平稳加速度时程，若是，即直接生成非平稳加速时程，即初始无脉冲地震动Xa；若为平稳加速度时程，则根据式(8)生成初始无脉冲地震动Xa：
Xa＝xa(t)×g(t) (8)
式中：Xa为初始无脉冲地震动；xa(t)为加速度时程；g(t)为强度包络函数；t为时间；T1为平稳段的起始时刻；T2为平稳段的结束时刻；c为衰减阶段的变化速率。
3.根据权利要求1所述的一种基于改进时域叠加法的近断层地震动拟合方法，其特征在于，所述对初始无脉冲地震动进行低频滤波生成近断层地震动时程的高频成分的具体实现方式为：对生成的初始无脉冲地震动做傅里叶变换，将频率区间[0,fr]所对应的傅里叶幅值谱的值设为零，再经过傅里叶逆变换得到近断层地震动时程的高频成分Ah(t)。
4.根据权利要求1所述的一种基于改进时域叠加法的近断层地震动拟合方法，其特征在于，所述使用等效速度脉冲模型生成近断层地震动时程的低频成分的具体实现步骤如下：
1)根据式(10)生成低频脉冲速度时程ν(t)
式中：ν(t)为低频脉冲速度时程；ω(t)为低频脉冲速度时程包络函数；Vp为脉冲峰值；t为时间；fp为脉冲频率；取1/Tp，Tp为脉冲周期；t1为余弦函数峰值发生时刻；T为结构自振周期；γ为衰减速率；t0为包络函数峰值发生时刻；
2)对生成的低频脉冲速度时程进行一次求导得到低频成分脉冲加速度时程α(t)；
3)得到近断层地震动时程的低频成分Al(t)：将得到的低频成分脉冲加速度时程在频率区间为[fr,∞)所对应的傅里叶幅值谱的值设为零，再经过傅里叶逆变换得到近断层地震动时程的低频成分Al(t)。
5.根据权利要求1所述的一种基于改进时域叠加法的近断层地震动拟合方法，其特征在于，所述对近断层地震动时程的高频成分和低频成分在时域叠加过程中的高低频峰值时刻提出了明确定义，生成新的近断层地震动低频成分并与近断层地震动高频成分在时域上进行叠加生成原始地震动反应谱的具体实现步骤如下：
1)计算近断层地震动时程的高频成分Ah(t)峰值时刻tha，如式(12)所示；
ln(tha)＝1.35MW-6.88 (12)
其中：Mw为矩震级；tha为近断层地震动时程的高频成分Ah(t)峰值时刻；
2)计算近断层地震动时程的低频成分Al(t)峰值时刻tla，如式(13)所示；
ln(tla)＝1.35MW-6.88 (13)
其中：Mw为矩震级；tla为近断层地震动时程的低频成分Al(t)峰值时刻；
3)由高频峰值时刻tha和低频峰值时刻tla求得高低频峰值时刻差σt，如式(14)所示；
σt＝tha-tla (14)
4)把近断层地震动时程的低频成分Al(t)在时间轴上平移σt，得到新的近断层地震动时程的低频成分All(t)；
5)将新的近断层地震动时程的低频成分All(t)与近断层地震动时程的高频成分Ah(t)在时域进行叠加，得到原始地震动的反应谱Sa。
6.根据权利要求1所述的一种基于改进时域叠加法的近断层地震动拟合方法，其特征在于，对近断层地震动低频成分并与近断层地震动高频成分叠加后的原始地震动反应谱采用FFT法进行反应谱调整，生成与目标反应谱更加吻合的近断层地震动的具体实现步骤如下：
1)求原始地震动的反应谱和目标加速度反应谱的平均相对误差；
2)判断得到的平均相对误差是否小于固定阈值；若平均相对误差大于固定阈值，则计算原始地震动的反应谱和目标加速度反应谱的比值，调整傅里叶幅值谱，进入步骤3)；若平均相对误差小于等于固定阈值，则输出近断层脉冲型地震动A(t)；
3)计算相位谱：重新生成在-2π到2π上均匀分布的相位谱；
4)根据式(15)生成新的加速度时程xp(t)
式中：φk为相位角；ωk为圆频率；A(ωk)为幅值谱；N为目标反应谱的频域中频率的分隔点数；g(t)为强度包络函数；xp(t)为新的加速度时程；
5)判断生成的加速度时程是否为非平稳加速度时程，若是，即直接生成非平稳加速时程，即新的地震动Xp；若为平稳加速度时程，根据式(16)生成新的地震动Xp
Xp＝xp(t)×g(t) (16)
式中：Xp为新的地震动；g(t)为地震动包络函数；t为时间；T1为平稳段的起始时刻；T2为平稳段的结束时刻；c为衰减阶段的变化速率；
6)根据生成的新的地震动计算地震动反应谱；求地震动的反应谱和目标加速度反应谱的平均相对误差，并判断其是否小于固定阈值；若平均相对误差大于固定阈值，则重复步骤3)～步骤5)；若平均相对误差小于等于固定阈值，则输出近断层脉冲型地震动A(t)。
CN202210433801.1A 2022-04-24 2022-04-24 一种基于改进时域叠加法的近断层地震动拟合方法 Active CN114966835B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210433801.1A CN114966835B (zh) 2022-04-24 2022-04-24 一种基于改进时域叠加法的近断层地震动拟合方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210433801.1A CN114966835B (zh) 2022-04-24 2022-04-24 一种基于改进时域叠加法的近断层地震动拟合方法

Publications (2)

Publication Number Publication Date
CN114966835A CN114966835A (zh) 2022-08-30
CN114966835B true CN114966835B (zh) 2024-06-04

Family

ID=82979322

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210433801.1A Active CN114966835B (zh) 2022-04-24 2022-04-24 一种基于改进时域叠加法的近断层地震动拟合方法

Country Status (1)

Country Link
CN (1) CN114966835B (zh)

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR20070115523A (ko) * 2006-06-02 2007-12-06 박은희 동적시스템의 실시간 해석을 위한 정확한 저주파수 및스펙트럼 위상 추정 방법
CN107480391A (zh) * 2017-08-24 2017-12-15 中铁二院工程集团有限责任公司 基于数据驱动的近断层非平稳地震动模拟方法
CN110389379A (zh) * 2019-07-12 2019-10-29 中国地震局地球物理研究所 能够表征地面永久位移的近断层地震加速度时程拟合方法
CN110794459A (zh) * 2019-11-14 2020-02-14 福州大学 一种海底近断层地震动的拟合方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110390169A (zh) * 2019-07-25 2019-10-29 西南交通大学 一种基于力学模型的地震滑坡危险性定量评价方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR20070115523A (ko) * 2006-06-02 2007-12-06 박은희 동적시스템의 실시간 해석을 위한 정확한 저주파수 및스펙트럼 위상 추정 방법
CN107480391A (zh) * 2017-08-24 2017-12-15 中铁二院工程集团有限责任公司 基于数据驱动的近断层非平稳地震动模拟方法
CN110389379A (zh) * 2019-07-12 2019-10-29 中国地震局地球物理研究所 能够表征地面永久位移的近断层地震加速度时程拟合方法
CN110794459A (zh) * 2019-11-14 2020-02-14 福州大学 一种海底近断层地震动的拟合方法

Also Published As

Publication number Publication date
CN114966835A (zh) 2022-08-30

Similar Documents

Publication Publication Date Title
Tu et al. On the nature of compressive fluctuations in the solar wind
CN110794459B (zh) 一种海底近断层地震动的拟合方法
CN114442153B (zh) 一种近断层地震动拟合方法
CN107967240B (zh) 基于人工地震动生成过程的耐震时程加速度优化算法
Jeon et al. A narrowband active noise control system with frequency corrector
CN104616667A (zh) 一种用于汽车内的主动降噪方法
CN109002678B (zh) 基于水电工程水工建筑物抗震设计的人工地震模拟方法
CN114578417B (zh) 一种基于小波调整的脉冲型地震动拟合方法
CN104533732A (zh) 一种抑制风力发电机组塔架左右振动的控制方法及装置
CN111986644B (zh) 一种宽窄带混合前馈型主动噪声控制系统
US20240106360A1 (en) Method for reducing electromagnetic vibration of permanent magnet motor with fractional slot concentrated winding
CN103972912B (zh) 一种含风电电力系统频率响应的频域分析方法
CN114966835B (zh) 一种基于改进时域叠加法的近断层地震动拟合方法
CN108169794B (zh) 近断层人工地震波与规范反应谱匹配的调整方法
CN108736513B (zh) 考虑秒级风速相关性的年度的风电场频率风险评估方法
CN109828154A (zh) 一种分频段复合正交脉冲注入的三相电网阻抗测量方法
CN113484911B (zh) 一种基于地震动演化功率谱的地震位移反应计算方法
CN115236732B (zh) 一种基于相位谱重构的近断层地震动拟合方法
CN113255482A (zh) 一种基于hht脉冲参数识别的远场类谐和地震动合成方法
CN110826197B (zh) 一种基于改进Cholesky分解闭合解的风速场模拟方法
CN102928817A (zh) 一种利用时延估计进行转子碰摩声发射源定位的方法
Shama Spectrum compatible earthquake ground motions by morlet wavelet
CN115510381B (zh) 一种海上风机多元相干效应风场载荷构建方法
CN114459712B (zh) 一种基于自回归模型的地震模拟振动台实验补偿方法
CN104506188A (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
GR01 Patent grant
GR01 Patent grant