专利名称：获得地震动样本的方法和装置的制作方法
技术领域：
本发明涉及地震工程领域，尤其涉及一种获得地震动样本的方法和装置。
背景技术：
重要工程结构，如超高层建筑、核电站、海洋平台、大跨度桥梁、大型水坝和储液罐等的抗震设计需要进行动力时程分析，以确定工程结构在可能遭受的地震作用下动力响应的全过程，进而实现抗震设防目标。动力时程分析需要一定数量的、并且满足一定地震环境及场地条件的地震动加速度时程作为结构动力反应分析的输入。然而，天然的强震观测记录数量非常有限，无法满足实际工程的多种需求，因此需要利用数值方法进行地震动的模拟，以生成大量满足一定工程特性的样本地震动。传统的非平稳地震动模拟方法可大致分为4大类基于时域调制的方法、基于演 变谱的方法、基于相位差谱的方法以及基于自回归滑移平均模型(ARMA)的方法。上述大部分地震动随机模拟方法主要是将给定天然地震动视为某一随机过程的一个样本，对该样本进行数据分析，从中提取出某一复杂变化的物理量(通常是相位或相位差)作为描述地震动随机性的随机变量并研究其频数分布，进而构造该随机变量的概率模型，并据此最终确定地震动的随机过程模型。然而在该过程中，需要对原始地震动进行非常繁琐的参数估计。
发明内容
本发明的实施例提供一种获得地震动样本的方法和装置，能够获得与天然地震动具有相同非平稳特性的地震动样本，无需繁琐的参数估计过程。为达到上述目的，本发明的实施例采用如下技术方案一种获得地震动样本的方法，包括获取预先存储的天然地震动记录的瞬时幅值和瞬时相位；对所述瞬时相位进行解卷；采用傅立叶变换获取所述解卷后的瞬时相位的幅值谱和相位谱；采用所述相位谱的高频分量构造随机频率过程；采用所述随机频率过程、所述幅值谱及所述相位谱的低频分量构造随机相位过程;根据所述瞬时幅值与所述随机相位过程获得地震动随机过程；根据所述地震动随机过程获得地震动样本。一种获得地震动样本的装置，包括第一获取单元，用于从预先存储的天然地震动记录获取所述天然地震动记录的瞬时幅值和瞬时相位；解卷单元，用于对所述第一获取单元获取的瞬时相位进行解卷；第二获取单元，用于采用傅立叶变换获取所述解卷单元解卷后的瞬时相位的幅值谱和相位谱；
第一构造单元，用于采用所述第二获取单元获取的相位谱的高频分量构造随机频率过程；第二构造单元，用于采用所述第一构造单元构造的随机频率过程、所述第二获取单元获取的幅值谱及所述第二获取单元获取的相位谱的低频分量构造随机相位过程；第一获得单元，用于根据所述第一获取单元获取的瞬时幅值与所述第二构造单元构造的随机相位过程获得地震动随机过程；第二获得单元，用于根据所述地震动随机过程获得地震动样本。本发明实施例提供的获得地震动样本的方法和装置，通过获取预先存储的天然地震动记录的瞬时幅值和瞬时相位，进而构造随机相位过程，瞬时幅值能够反映天然地震动记录的强度非平稳特性，而随机相位过程则能够反映天然地震动记录的频率非平稳特性，最后根据所述瞬时幅值与所述随机相位过程获得地震动随机过程，进而获得与原始的天然 地震动记录具有相同非平稳特性的一系列地震动样本。
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。图I为本发明实施例提供的一种获取地震动样本的方法；图2为本发明实施例提供的另一种获取地震动样本的方法；图3为本发明实施例提供的一种获取地震动样本的装置。
具体实施例方式下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。本发明实施例提供了一种获得地震动样本的方法，如图I所示，包括101、获取预先存储的天然地震动记录的瞬时幅值和瞬时相位。所述预先存储的天然地震动记录指与要进行抗震设计的工程结构所在的场地具有相似地震环境和场地条件的天然地震动记录，其中天然地震动记录指利用强震仪记录到的天然地震产生的地震动加速度、速度或位移的时间过程。可以采用多种方式获取预先存储的天然地震动记录的瞬时幅值和瞬时相位，其中优选的，可采用希尔伯特(Hilbert)变换获取所述预先存储的天然地震动记录的瞬时幅值和瞬时相位。当然本发明实施例对此是不作限定的，本领域技术人员可根据本领域公知常识或常用技术手段获取所述预先存储的天然地震动记录的瞬时幅值和瞬时相位。具体但不限于此的，所述采用希尔伯特变换从预先存储的天然地震动记录获取所述天然地震动记录的瞬时幅值和瞬时相位包括使用下列公式获取所述天然地震动记录的瞬时幅值和瞬时相位，
_=_=arctan载，少"(/)=啦(/)]=士pi^r;其中，a^t)表示瞬时幅值，9 / 0(t)表示瞬时相位，Xtl (t)表示预先存储的天然地震动记录，H表示希尔伯特变换，P表示取广义积分的柯西主值。该步骤中，首先对预先存储的天然地震动记录进行希尔伯特变换，得到变换结果Y0⑴，然后再分别通过公式&= jXo2(0 + w(0和砵(0 = arctan，获得瞬时幅值％⑴和瞬时相位0 ' o(t)。102、对所述瞬时相位进行解卷；对所述瞬时相位进行解卷以获得解卷后的瞬时相位。103、采用傅立叶(Fourier)变换获取所述解卷后的瞬时相位的幅值谱和相位谱。在本步骤中，可以基于傅立叶变换获取所述解卷后的瞬时相位的幅值谱和相位 P曰。具体但不限于此的，所述傅立叶变换获取所述解卷后的瞬时相位的幅值谱和相位谱包括使用下列公式获取所述解卷后的瞬时相位的幅值谱和相位谱，4,,(/) = VRe2[^(./)]Msn2[^(./')]，= 31013,1 Re[A(/)]，n0(f) = J■二0Q(t)exp(-HnfOdt ；其中，Ae；0(f)表示解卷后的瞬时相位的幅值谱，O "(f)表示解卷后的瞬时相位的相位谱，ha)表示解卷后的瞬时相位，Re表示取复数的实部，Im表示取复数的虚部，i表示虚数单位，I2=-Io该步骤中，首先对解卷后的瞬时相位进行傅立叶变换，得到傅立叶变换结果\(f)，然后分别获取傅立叶变换结果的实部和虚部，最后分别通过公式
^ w、/)1
4奶= >似。(/)] +如似0(/)]和‘(/) = ■她获得解卷后的瞬时雛的
幅值谱A "(f)和相位谱Q e’0(f)。104、采用所述相位谱的高频分量构造随机频率过程。其中，相位谱的高频分量可以通过多种滤波方法对所述相位谱进行滤波获得。具体但不限于此的，本步骤中使用下列公式构造随机频率过程，^ ( f ) = A ( f ) 0 ( f ) , AU、=权JJ) +奶JJ)，
-I (P.., (「、i
^,(/) = H[巧 (/)] = - P ,；
JT 厂丫 f-T其中，W (f)表示随机频率过程，Oo.Jf)表示相位谱的高频分量，O (f)表示随机频率函数，H表示希尔伯特变换，P表示取广义积分的柯西主值。该步骤中，首先对相位谱的高频分量O e，d(f)进行希尔伯特变换，得到希尔伯特变换结果,然后通过公式廣/) =(/)得到确定性函数A (f)，最后通过公式W (f) = A(f) 0 (f)，即确定性函数A(f)与随机频率函数 (f)的乘积得到随机频率过程W (f)。105、采用所述随机频率过程、所述幅值谱及所述相位谱的低频分量构造随机相位过程。其中，所述相位谱的低频分量可通过多种滤波方法对所述相位谱进行滤波获得。在本步骤中，所述随机频率过程、所述幅值谱及所述相位谱的低频分量均由预先存储的天然地震动的瞬时相位直接或间接得到，构造的随机相位过程反映预先存储的天然地震动记录的频率非平稳特性。具体但不限于此的，使用下列公式构造随机相位过程，^(/) = | Q(/)exp(i27u//)d/ Q (f) = A90 (f) exp {i [O 0，d(f)+W (f)]}；其中，0 (t)表示随机相位过程，Ae；0(f)表示解卷后的瞬时相位的幅值谱， O 0；d(f)表示相位谱的低频分量，(f)表示随机频率过程，i表示虚数单位，I2=-I ;106、根据所述瞬时幅值与所述随机相位过程获得地震动随机过程。根据所述瞬时幅值与所述随机相位过程获得地震动随机过程，具体但不限于此的，可采用将所述瞬时幅值与所述随机相位过程的余弦值相乘的方式获得地震动随机过程，也即采用下列公式获得地震动随机过程，X (t) = a0 (t) cos 0 (t)；其中，x(t)表示地震动随机过程,ajt)表示瞬时幅值，0 (t)表示随机相位过程。107、根据所述地震动随机过程获得地震动样本。为随机相位过程设定不同的随机数，利用伪随机数生成算法即可生成与所述给定随机数相对应的随机频率函数的一个样本，进而产生与所述随机数对应的地震动随机过程的一个样本，进而针对不同随机数可得到一系列的用于需要进行抗震设计的工程结构的输入地震动样本。可以理解的是，本发明实施例中的步骤101至106的内容及顺序是不作严格限定的。例如，步骤101可以为获取预先存储的天然地震动记录的瞬时相位，然后在步骤106之前可以增加获取预先存储的天然地震动记录的瞬时幅值的步骤。又如，步骤103可以为采用傅立叶变换获取所述解卷后的瞬时相位的相位谱，然后在步骤105之前可以增加获取采用傅立叶变换获取所述解卷后的瞬时相位的幅值谱的步骤。本发明实施例提供的获得地震动样本的方法，通过获取预先存储的天然地震动记录的瞬时幅值和瞬时相位，进而构造随机相位过程，瞬时幅值能够反映天然地震动记录的强度非平稳特性，而随机相位过程则能够反映天然地震动记录的频率非平稳特性，最后根据所述瞬时幅值与所述随机相位过程获得地震动随机过程，进而获得地震动样本，即获得与原始的天然地震动记录具有相同非平稳特性的一系列地震动样本。通过本发明实施例提供的方法获得地震动样本可作为输入应用于利用蒙特卡洛法对结构地震反应的数值模拟研究工作中，以使得工程结构能够更加合理的抵御可能发生的地震作用。进一步的，在本发明提供的一个实施例中，在步骤104之前还包括采用汉宁(Hanning)窗函数对所述相位谱分解,得到相位谱的高频分量和低频分量。具体但不限于此的，使用下列公式得到相位谱的高频分量和低频分量，O 0；d(f) = Hann[O 0’0(f) ], CVn(f) = O 0；o(f)-O 0；d(f)；
其中，O 0，d (f)表示相位谱的低频分量，①0，n (f)表示相位谱的高频分量，o e,0(f)表示解卷后的瞬时相位的相位谱，Hann表示利用汉宁窗函数对信号进行低通滤波处理。进一步的，在本发明提供的一个实施例中，步骤104中的所述随机频率函数 (f)为在区间[-1，I]均匀分布的随机频率函数，并满足对任意频率f，0 (f)为区间[-1，I]均匀分布的随机变量；及若f2,则0 (f:)与0 (f2)为相互独立的随机变量。针对任意给定的随机数，利用伪随机数生成算法即可生成与所述给定随机数相对应的随机频率函数 (f)的一个样本。 为了更好的说明本发明提供的获得地震动样本的方法，本发明提供了一个优选的具体实施例，该实施例提供了一种获得地震动样本的方法，如图2所示，包括201、采用希尔伯特变换获取所述预先存储的天然地震动记录XtlU)的瞬时幅值a0(t)和瞬时相位0 ' 0(t) o首先对预先存储的天然地震动记录进行希尔伯特变换，得到变换结果yci(t)V0(Z) 二 H[x0(/)j 二士 Pj' '。
TZt — T然后获取天然地震动记录X。(t)的瞬时幅值aQ(t)a0(t) = ylxl(t) + y；(t),和天然地震动记录X。(t)的瞬时相位0 ' 0(t)
,V (/I(/) = arctan —~ ,,
x0(0202、对所述瞬时相位0 ' 0(t)进行解卷，得到解卷后的瞬时相位0力)。203、采用傅立叶变换获取所述解卷后的瞬时相位0 Jt)的幅值谱Ae^f)和相位谱 O0 o(f)。首先对解卷后的瞬时相位ejt)进行傅立叶变换，得到傅立叶变换结果QJf):i20(/) = J ^ O0 (t)exp( -ilnf t)dt。然后获取解卷后的瞬时相位的幅值谱A" (f)AoxXf) = VRe:[A(/)] +Inr[A(./)]。获取解卷后的瞬时相位的相位谱O "(f)
^ . r、Im[ A( /')]巧(!(/)= arctan .1;；。
RefA (/)1其中，Re表示取复数的实部，Im表示取复数的虚部，i表示虚数单位，i2=_l。204、采用汉宁窗函数获取相位谱①"(f)的高频分量O0,n(f)和低频分量① e’d(f)。首先采用汉宁窗函数对相位谱0"(f)进行低通滤波，获取其低频分量① e’d(f)O0；d(f) = Hann
0然后通过相位谱O "(f)和低频分量①0，d(f)获取高频分量①0，n(f)
①0’n(f)=① 0’o(f)-C> 0’d(f)。205、采用所述相位谱的高频分量①e ,n(f)构造随机频率过程W (f)。首先对高频分量O e ,n(f)进行希尔伯特变换，得到希尔伯特变换结果為.n(/)
权利要求
1.一种获得地震动样本的方法，其特征在于，包括 获取预先存储的天然地震动记录的瞬时幅值和瞬时相位； 对所述瞬时相位进行解卷； 采用傅立叶变换获取所述解卷后的瞬时相位的幅值谱和相位谱； 采用所述相位谱的高频分量构造随机频率过程； 采用所述随机频率过程、所述幅值谱及所述相位谱的低频分量构造随机相位过程； 根据所述瞬时幅值与所述随机相位过程获得地震动随机过程； 根据所述地震动随机过程获得地震动样本。
2.根据权利要求I所述的获得地震动样本的方法，其特征在于，所述获取预先存储的天然地震动记录的瞬时幅值和瞬时相位包括 采用希尔伯特变换获取预先存储的天然地震动记录的瞬时幅值和瞬时相位。
3.根据权利要求I所述的获得地震动样本的方法，其特征在于，在所述采用所述相位谱的高频分量构造随机频率过程前，还包括 采用汉宁窗函数对所述相位谱分解，得到相位谱的高频分量和低频分量。
4.根据权利要求2所述的获得地震动样本的方法，其特征在于，所述采用希尔伯特变换从预先存储的天然地震动记录获取所述天然地震动记录的瞬时幅值和瞬时相位包括 使用下列公式获取所述天然地震动记录的瞬时幅值和瞬时相位， 。。(0=，=arctan^S ’刷=HM)]=M 二; 其中,ajt)表示瞬时幅值，0 ' C1 (t)表示瞬时相位,XtlU)表示预先存储的天然地震动记录，H表示希尔伯特变换，P表示取广义积分的柯西主值。
5.根据权利要求I所述的获得地震动样本的方法，其特征在于，所述采用傅立叶变换获取所述解卷后的瞬时相位的幅值谱和相位谱包括 使用下列公式获取所述解卷后的瞬时相位的幅值谱和相位谱，； 其中，AtJf)表示解卷后的瞬时相位的幅值谱，O "(f)表示解卷后的瞬时相位的相位谱，0。(0表示解卷后的瞬时相位，Re表示取复数的实部，Im表示取复数的虚部，i表示虚数单位，I2=-Io
6.根据权利要求3所述的获得地震动样本的方法，其特征在于，所述采用汉宁窗函数对所述相位谱分解，得到相位谱的高频分量和低频分量包括 使用下列公式得到相位谱的高频分量和低频分量， 其中，^ e,d(f)表示相位谱的低频分量，^ e,n(f)表示相位谱的高频分量，^ e,o (f)表示解卷后的瞬时相位的相位谱，Hann表示利用汉宁窗函数对信号进行低通滤波处理。
7.根据权利要求I所述的获得地震动样本的方法，其特征在于，所述采用所述相位谱的高频分量构造随机频率过程包括 使用下列公式构造随机频率过程，( f ) = A ( f ) 0 ( f ) , A(f) = pl(f) + 0ln(f)，-I [.OG 0n n (t)^oAf) = H[^,(/)] = -Pj ,； Tl j^J-T 其中，W (f)表示随机频率过程，表示相位谱的高频分量， (f)表示随机频率函数，H表示希尔伯特变换，P表示取广义积分的柯西主值。
8.根据权利要求7所述的获得地震动样本的方法，其特征在于，所述随机频率函数 (f)为在区间[-1，I]均匀分布的随机频率函数，并满足 对任意频率f， (f)为区间[-1，I]均匀分布的随机变量 '及 若古f2，则 (f\)与 (f2)为相互独立的随机变量。
9.根据权利要求I所述的获得地震动样本的方法，其特征在于，所述采用所述随机频率过程、所述幅值谱及所述相位谱的低频分量构造随机相位过程包括 使用下列公式构造随机相位过程， ^(0 = J x 0(./')exp(i2TC,/'/)d/ , Q (f) = A0；o(f) exp {i [O 0，d(f)+W (f)]}； 其中，0 (t)表示随机相位过程，A"(f)表示解卷后的瞬时相位的幅值谱，^Vd(f)表示相位谱的低频分量，(f)表示随机频率过程，i表示虚数单位，I2=-Io
10.根据权利要求I所述的获得地震动样本的方法，其特征在于，所述根据所述瞬时幅值与所述随机相位过程获得地震动随机过程包括 使用下列公式构造地震动随机过程X (t) = a0(t)cos 0 (t)； 其中，x(t)表示地震动随机过程，ajt)表示瞬时幅值，0⑴表示随机相位过程。
11.一种获得地震动样本的装置，其特征在于，包括 第一获取单元，用于从预先存储的天然地震动记录获取所述天然地震动记录的瞬时幅值和瞬时相位； 解卷单元，用于对所述第一获取单元获取的瞬时相位进行解卷； 第二获取单元，用于采用傅立叶变换获取所述解卷单元解卷后的瞬时相位的幅值谱和相位谱； 第一构造单元，用于采用所述第二获取单元获取的相位谱的高频分量构造随机频率过程; 第二构造单元，用于采用所述第一构造单元构造的随机频率过程、所述第二获取单元获取的幅值谱及所述第二获取单元获取的相位谱的低频分量构造随机相位过程； 第一获得单元，用于根据所述第一获取单元获取的瞬时幅值与所述第二构造单元构造的随机相位过程获得地震动随机过程； 第二获得单元，用于根据所述地震动随机过程获得地震动样本。
全文摘要
本发明实施例提供了一种获得地震动样本的方法和装置，涉及地震工程领域，能够获得与天然地震动具有相同非平稳特性的地震动样本，所述方法包括获取预先存储的天然地震动记录的瞬时幅值和瞬时相位；对所述瞬时相位进行解卷；采用傅立叶(Fourier)变换获取所述解卷后的瞬时相位的幅值谱和相位谱；采用所述相位谱的高频分量构造随机频率过程；采用所述随机频率过程、所述幅值谱及所述相位谱的低频分量构造随机相位过程；根据所述瞬时幅值与所述随机相位过程获得地震动随机过程；根据所述地震动随机过程获得地震动样本。本发明可用于抗震设计领域。
文档编号G01V1/28GK102721978SQ20121021029
公开日2012年10月10日 申请日期2012年6月19日 优先权日2012年6月19日
发明者张郁山 申请人:中国地震灾害防御中心