CN102721978A - 获得地震动样本的方法和装置 - Google Patents

Abstract

本发明实施例提供了一种获得地震动样本的方法和装置，涉及地震工程领域，能够获得与天然地震动具有相同非平稳特性的地震动样本，所述方法包括：获取预先存储的天然地震动记录的瞬时幅值和瞬时相位；对所述瞬时相位进行解卷；采用傅立叶（Fourier）变换获取所述解卷后的瞬时相位的幅值谱和相位谱；采用所述相位谱的高频分量构造随机频率过程；采用所述随机频率过程、所述幅值谱及所述相位谱的低频分量构造随机相位过程；根据所述瞬时幅值与所述随机相位过程获得地震动随机过程；根据所述地震动随机过程获得地震动样本。本发明可用于抗震设计领域。

Description

获得地震动样本的方法和装置

技术领域

本发明涉及地震工程领域，尤其涉及一种获得地震动样本的方法和装置。

背景技术

重要工程结构，如超高层建筑、核电站、海洋平台、大跨度桥梁、大型水坝和储液罐等的抗震设计需要进行动力时程分析，以确定工程结构在可能遭受的地震作用下动力响应的全过程，进而实现抗震设防目标。动力时程分析需要一定数量的、并且满足一定地震环境及场地条件的地震动加速度时程作为结构动力反应分析的输入。然而，天然的强震观测记录数量非常有限，无法满足实际工程的多种需求，因此需要利用数值方法进行地震动的模拟，以生成大量满足一定工程特性的样本地震动。

传统的非平稳地震动模拟方法可大致分为4大类：基于时域调制的方法、基于演变谱的方法、基于相位差谱的方法以及基于自回归滑移平均模型（ARMA）的方法。上述大部分地震动随机模拟方法主要是将给定天然地震动视为某一随机过程的一个样本，对该样本进行数据分析，从中提取出某一复杂变化的物理量（通常是相位或相位差）作为描述地震动随机性的随机变量并研究其频数分布，进而构造该随机变量的概率模型，并据此最终确定地震动的随机过程模型。然而在该过程中，需要对原始地震动进行非常繁琐的参数估计。

发明内容

本发明的实施例提供一种获得地震动样本的方法和装置，能够获得与天然地震动具有相同非平稳特性的地震动样本，无需繁琐的参数估计过程。

为达到上述目的，本发明的实施例采用如下技术方案：

一种获得地震动样本的方法，包括：

获取预先存储的天然地震动记录的瞬时幅值和瞬时相位；

对所述瞬时相位进行解卷；

采用傅立叶变换获取所述解卷后的瞬时相位的幅值谱和相位谱；

采用所述相位谱的高频分量构造随机频率过程；

采用所述随机频率过程、所述幅值谱及所述相位谱的低频分量构造随机相位过程；

根据所述瞬时幅值与所述随机相位过程获得地震动随机过程；

根据所述地震动随机过程获得地震动样本。

一种获得地震动样本的装置，包括：

第一获取单元，用于从预先存储的天然地震动记录获取所述天然地震动记录的瞬时幅值和瞬时相位；

解卷单元，用于对所述第一获取单元获取的瞬时相位进行解卷；

第二获取单元，用于采用傅立叶变换获取所述解卷单元解卷后的瞬时相位的幅值谱和相位谱；

第一构造单元，用于采用所述第二获取单元获取的相位谱的高频分量构造随机频率过程；

第二构造单元，用于采用所述第一构造单元构造的随机频率过程、所述第二获取单元获取的幅值谱及所述第二获取单元获取的相位谱的低频分量构造随机相位过程；

第一获得单元，用于根据所述第一获取单元获取的瞬时幅值与所述第二构造单元构造的随机相位过程获得地震动随机过程；

第二获得单元，用于根据所述地震动随机过程获得地震动样本。

本发明实施例提供的获得地震动样本的方法和装置，通过获取预先存储的天然地震动记录的瞬时幅值和瞬时相位，进而构造随机相位过程，瞬时幅值能够反映天然地震动记录的强度非平稳特性，而随机相位过程则能够反映天然地震动记录的频率非平稳特性，最后根据所述瞬时幅值与所述随机相位过程获得地震动随机过程，进而获得与原始的天然地震动记录具有相同非平稳特性的一系列地震动样本。

附图说明

为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。

图1为本发明实施例提供的一种获取地震动样本的方法；

图2为本发明实施例提供的另一种获取地震动样本的方法；

图3为本发明实施例提供的一种获取地震动样本的装置。

具体实施方式

下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。

本发明实施例提供了一种获得地震动样本的方法，如图1所示，包括：

101、获取预先存储的天然地震动记录的瞬时幅值和瞬时相位。

所述预先存储的天然地震动记录指与要进行抗震设计的工程结构所在的场地具有相似地震环境和场地条件的天然地震动记录，其中天然地震动记录指利用强震仪记录到的天然地震产生的地震动加速度、速度或位移的时间过程。

可以采用多种方式获取预先存储的天然地震动记录的瞬时幅值和瞬时相位，其中优选的，可采用希尔伯特（Hilbert）变换获取所述预先存储的天然地震动记录的瞬时幅值和瞬时相位。当然本发明实施例对此是不作限定的，本领域技术人员可根据本领域公知常识或常用技术手段获取所述预先存储的天然地震动记录的瞬时幅值和瞬时相位。

具体但不限于此的，所述采用希尔伯特变换从预先存储的天然地震动记录获取所述天然地震动记录的瞬时幅值和瞬时相位包括：

使用下列公式获取所述天然地震动记录的瞬时幅值和瞬时相位，

$a_0\left(t\right)=\sqrt{x_0^2\left(t\right)+y_0^2\left(t\right)},$ ${\mathrm{\θ}}_0^{\′}\left(t\right)=\arctan \frac{y_0\left(t\right)}{x_0\left(t\right)},$ $y_0\left(t\right)=H\left[x_0\left(t\right)\right]=\frac{1}{\mathrm{\π}}P{\∫}_{-\∞}^{\∞}\frac{x_0\left(\mathrm{\τ}\right)}{t-\mathrm{\τ}}\mathrm{d\τ};$

其中，a₀(t)表示瞬时幅值，θ′₀(t)表示瞬时相位，x₀(t)表示预先存储的天然地震动记录，H表示希尔伯特变换，P表示取广义积分的柯西主值。

该步骤中，首先对预先存储的天然地震动记录进行希尔伯特变换，得到变换结果y₀(t)，然后再分别通过公式 $a_0\left(t\right)=\sqrt{x_0^2\left(t\right)+y_0^2\left(t\right)}$ 和 ${\mathrm{\θ}}_0^{\′}=\arctan \frac{y_0\left(t\right)}{x_0\left(t\right)}$ 获得瞬时幅值a₀(t)和瞬时相位θ′₀(t)。

102、对所述瞬时相位进行解卷；

对所述瞬时相位进行解卷以获得解卷后的瞬时相位。

103、采用傅立叶（Fourier）变换获取所述解卷后的瞬时相位的幅值谱和相位谱。

在本步骤中，可以基于傅立叶变换获取所述解卷后的瞬时相位的幅值谱和相位谱。

具体但不限于此的，所述傅立叶变换获取所述解卷后的瞬时相位的幅值谱和相位谱包括：

使用下列公式获取所述解卷后的瞬时相位的幅值谱和相位谱，

$A_{\mathrm{\θ},0}\left(f\right)=\sqrt{{\mathrm{Re}}^2\left[{\mathrm{\Ω}}_0\left(f\right)\right]+{\mathrm{Im}}^2\left[{\mathrm{\Ω}}_0\left(f\right)\right]},$ ${\mathrm{\Φ}}_{\mathrm{\θ},0}\left(f\right)=\arctan \frac{\mathrm{Im}\left[{\mathrm{\Ω}}_0\left(f\right)\right]}{\mathrm{Re}\left[{\mathrm{\Ω}}_0\left(f\right)\right]},$

${\mathrm{\Ω}}_0\left(f\right)={\∫}_{-\∞}^{\∞}{\mathrm{\θ}}_0\left(t\right)\exp (-i2\mathrm{\πft})\mathrm{dt};$

其中，A_θ,0(f)表示解卷后的瞬时相位的幅值谱，Φ_θ,0(f)表示解卷后的瞬时相位的相位谱，θ₀(t)表示解卷后的瞬时相位，Re表示取复数的实部，Im表示取复数的虚部，i表示虚数单位，i²=-1。

该步骤中，首先对解卷后的瞬时相位进行傅立叶变换，得到傅立叶变换结果Ω₀(f)，然后分别获取傅立叶变换结果的实部和虚部，最后分别通过公式 $A_{\mathrm{\θ},0}\left(f\right)=\sqrt{{\mathrm{Re}}^2\left[{\mathrm{\Ω}}_0\left(f\right)\right]+{\mathrm{Im}}^2\left[{\mathrm{\Ω}}_0\left(f\right)\right]}$ 和 ${\mathrm{\Φ}}_{\mathrm{\θ},0}\left(f\right)=\arctan \frac{\mathrm{Im}\left[{\mathrm{\Ω}}_0\left(f\right)\right]}{\mathrm{Re}\left[{\mathrm{\Ω}}_0\left(f\right)\right]}$ 获得解卷后的瞬时相位的幅值谱A_θ,0(f)和相位谱Ω_θ,0(f)。

104、采用所述相位谱的高频分量构造随机频率过程。

其中，相位谱的高频分量可以通过多种滤波方法对所述相位谱进行滤波获得。

具体但不限于此的，本步骤中使用下列公式构造随机频率过程，

Ψ(f)＝A(f)Θ(f)， $A\left(f\right)=\sqrt{{\mathrm{\Φ}}_{\mathrm{\θ},n}^2\left(f\right)+{\widehat{\mathrm{\Φ}}}_{\mathrm{\θ},n}^2\left(f\right)},$ ${\widehat{\mathrm{\Φ}}}_{\mathrm{\θ},n}\left(f\right)=H\left[{\mathrm{\Φ}}_{\mathrm{\θ},n}\left(f\right)\right]=\frac{1}{\mathrm{\π}}P{\∫}_{-\∞}^{\∞}\frac{{\mathrm{\Φ}}_{\mathrm{\θ},n}\left(\mathrm{\τ}\right)}{f-\mathrm{\τ}}\mathrm{d\τ};$

其中，Ψ(f)表示随机频率过程，Φ_θ，n(f)表示相位谱的高频分量，Θ(f)表示随机频率函数，H表示希尔伯特变换，P表示取广义积分的柯西主值。

该步骤中，首先对相位谱的高频分量Φ_θ,d(f)进行希尔伯特变换，得到希尔伯特变换结果

然后通过公式 $A\left(f\right)=\sqrt{{\mathrm{\Φ}}_{\mathrm{\θ},n}^2\left(f\right)+{\widehat{\mathrm{\Φ}}}_{\mathrm{\θ},n}^2\left(f\right)}$ 得到确定性函数A(f)，最后通过公式Ψ(f)＝A(f)Θ(f)，即确定性函数A(f)与随机频率函数Θ(f)的乘积得到随机频率过程Ψ(f)。

105、采用所述随机频率过程、所述幅值谱及所述相位谱的低频分量构造随机相位过程。

其中，所述相位谱的低频分量可通过多种滤波方法对所述相位谱进行滤波获得。在本步骤中，所述随机频率过程、所述幅值谱及所述相位谱的低频分量均由预先存储的天然地震动的瞬时相位直接或间接得到，构造的随机相位过程反映预先存储的天然地震动记录的频率非平稳特性。

具体但不限于此的，使用下列公式构造随机相位过程，

$\mathrm{\θ}\left(t\right)={\∫}_{-\∞}^{\∞}\mathrm{\Ω}\left(f\right)\exp \left(i2\mathrm{\πft}\right)\mathrm{dt},$ Ω(f)＝A_θ,0(f)exp{i[Φ_θ,d(f)+Ψ(f)]}；

其中，θ(t)表示随机相位过程，A_θ,0(f)表示解卷后的瞬时相位的幅值谱，Φ_θ,d(f)表示相位谱的低频分量，Ψ(f)表示随机频率过程，i表示虚数单位，i²=-1；

106、根据所述瞬时幅值与所述随机相位过程获得地震动随机过程。

根据所述瞬时幅值与所述随机相位过程获得地震动随机过程，具体但不限于此的，可采用将所述瞬时幅值与所述随机相位过程的余弦值相乘的方式获得地震动随机过程，也即采用下列公式获得地震动随机过程，

x(t)＝a₀(t)cosθ(t)；

其中，x(t)表示地震动随机过程，a₀(t)表示瞬时幅值，θ(t)表示随机相位过程。

107、根据所述地震动随机过程获得地震动样本。

为随机相位过程设定不同的随机数，利用伪随机数生成算法即可生成与所述给定随机数相对应的随机频率函数的一个样本，进而产生与所述随机数对应的地震动随机过程的一个样本，进而针对不同随机数可得到一系列的用于需要进行抗震设计的工程结构的输入地震动样本。

可以理解的是，本发明实施例中的步骤101至106的内容及顺序是不作严格限定的。例如，步骤101可以为获取预先存储的天然地震动记录的瞬时相位，然后在步骤106之前可以增加获取预先存储的天然地震动记录的瞬时幅值的步骤。又如，步骤103可以为采用傅立叶变换获取所述解卷后的瞬时相位的相位谱，然后在步骤105之前可以增加获取采用傅立叶变换获取所述解卷后的瞬时相位的幅值谱的步骤。

本发明实施例提供的获得地震动样本的方法，通过获取预先存储的天然地震动记录的瞬时幅值和瞬时相位，进而构造随机相位过程，瞬时幅值能够反映天然地震动记录的强度非平稳特性，而随机相位过程则能够反映天然地震动记录的频率非平稳特性，最后根据所述瞬时幅值与所述随机相位过程获得地震动随机过程，进而获得地震动样本，即获得与原始的天然地震动记录具有相同非平稳特性的一系列地震动样本。通过本发明实施例提供的方法获得地震动样本可作为输入应用于利用蒙特卡洛法对结构地震反应的数值模拟研究工作中，以使得工程结构能够更加合理的抵御可能发生的地震作用。

进一步的，在本发明提供的一个实施例中，在步骤104之前还包括：

采用汉宁（Hanning）窗函数对所述相位谱分解，得到相位谱的高频分量和低频分量。

具体但不限于此的，使用下列公式得到相位谱的高频分量和低频分量，

Φ_θ,d(f)＝Hann[Φ_θ,0(f)]，Φ_θ，n(f)＝Φ_θ,0(f)-Φ_θ,d(f)；

其中，Φ_θ,d(f)表示相位谱的低频分量，Φ_θ,n(f)表示相位谱的高频分量，Φ_θ,0(f)表示解卷后的瞬时相位的相位谱，Hann表示利用汉宁窗函数对信号进行低通滤波处理。

进一步的，在本发明提供的一个实施例中，步骤104中的所述随机频率函数Θ(f)为在区间[-1,1]均匀分布的随机频率函数，并满足：

对任意频率f，Θ(f)为区间[-1,1]均匀分布的随机变量；及

若f₁≠f₂，则Θ(f₁)与Θ(f₂)为相互独立的随机变量。

针对任意给定的随机数，利用伪随机数生成算法即可生成与所述给定随机数相对应的随机频率函数Θ(f)的一个样本。

为了更好的说明本发明提供的获得地震动样本的方法，本发明提供了一个优选的具体实施例，该实施例提供了一种获得地震动样本的方法，如图2所示，包括：

201、采用希尔伯特变换获取所述预先存储的天然地震动记录x₀(t)的瞬时幅值a₀(t)和瞬时相位θ′₀(t)。

首先对预先存储的天然地震动记录进行希尔伯特变换，得到变换结果y₀(t)：

$y_0\left(t\right)=H\left[x_0\left(t\right)\right]=\frac{1}{\mathrm{\π}}P{\∫}_{-\∞}^{\∞}\frac{x_0\left(\mathrm{\τ}\right)}{t-\mathrm{\τ}}\mathrm{d\τ}.$

然后获取天然地震动记录x₀(t)的瞬时幅值a₀(t)：

$a_0\left(t\right)=\sqrt{x_0^2\left(t\right)+y_0^2\left(t\right)},$

和天然地震动记录x₀(t)的瞬时相位θ′₀(t)：

${\mathrm{\θ}}_0^{\′}\left(t\right)=\arctan \frac{y_0\left(t\right)}{x_0\left(t\right)}.$

202、对所述瞬时相位θ′₀(t)进行解卷，得到解卷后的瞬时相位θ₀(t)。

203、采用傅立叶变换获取所述解卷后的瞬时相位θ₀(t)的幅值谱A_θ,0(f)和相位谱Φ_θ,0(f)。

首先对解卷后的瞬时相位θ₀(t)进行傅立叶变换，得到傅立叶变换结果Ω₀(f)：

${\mathrm{\Ω}}_0\left(f\right)={\∫}_{-\∞}^{\∞}{\mathrm{\θ}}_0\left(t\right)\exp (-i2\mathrm{\πft})\mathrm{dt}.$

然后获取解卷后的瞬时相位的幅值谱A_θ,0(f)：

$A_{\mathrm{\θ},0}\left(f\right)=\sqrt{{\mathrm{Re}}^2\left[{\mathrm{\Ω}}_0\left(f\right)\right]+I{m}^2\left[{\mathrm{\Ω}}_0\left(f\right)\right]}.$

获取解卷后的瞬时相位的相位谱Φ_θ,0(f)：

${\mathrm{\Φ}}_{\mathrm{\θ},0}\left(f\right)=\arctan \frac{\mathrm{Im}\left[{\mathrm{\Ω}}_0\left(f\right)\right]}{\mathrm{Re}\left[{\mathrm{\Ω}}_0\left(f\right)\right]}.$

其中，Re表示取复数的实部，Im表示取复数的虚部，i表示虚数单位，i²=-1。

204、采用汉宁窗函数获取相位谱Φ_θ,0(f)的高频分量Φ_θ，n(f)和低频分量Φ_θ,d(f)。

首先采用汉宁窗函数对相位谱Φ_θ,0(f)进行低通滤波，获取其低频分量Φ_θ,d(f)：

Φ_θ,d(f)＝Hann[Φ_θ,0(f)]。

然后通过相位谱Φ_θ,0(f)和低频分量Φ_θ,d(f)获取高频分量Φ_θ,n(f)：

Φ_θ,n(f)＝Φ_θ,0(f)-Φ_θ,d(f)。

205、采用所述相位谱的高频分量Φ_θ，n(f)构造随机频率过程Ψ(f)。

首先对高频分量Φ_θ，n(f)进行希尔伯特变换，得到希尔伯特变换结果

${\widehat{\mathrm{\Φ}}}_{\mathrm{\θ},n}\left(f\right)=H\left[{\mathrm{\Φ}}_{\mathrm{\θ},n}\left(f\right)\right]=\frac{1}{\mathrm{\π}}P{\∫}_{-\∞}^{\∞}\frac{{\mathrm{\Φ}}_{\mathrm{\θ},n}\left(\mathrm{\τ}\right)}{f-\mathrm{\τ}}\mathrm{d\τ}.$

然后构造确定性函数A(f)：

$A\left(f\right)=\sqrt{{\mathrm{\Φ}}_{\mathrm{\θ},n}^2\left(f\right)+{\widehat{\mathrm{\Φ}}}_{\mathrm{\θ},n}^2\left(f\right)}.$

根据确定性函数A(f)与随机频率函数Θ(f)获得随机频率过程Ψ(f)：

Ψ(f)＝A(f)Θ(f)。

在本步骤中，随机频率函数Θ(f)优选为在区间[-1,1]均匀分布的随机频率函数，并满足：

对任意频率f，Θ(f)为区间[-1,1]均匀分布的随机变量；及

若f₁≠f₂，则Θ(f₁)与Θ(f₂)为相互独立的随机变量。

针对任意给定的随机数，利用伪随机数生成算法即可生成与所述给定随机数相对应的随机频率函数Θ(f)的一个样本。

206、采用随机频率过程Ψ(f)、幅值谱A_θ,0(f)及相位谱的低频分量Φ_θ,d(f)构造随机相位过程θ(t)。

具体使用下列公式构造随机相位过程θ(t)：

$\mathrm{\θ}\left(t\right)={\∫}_{-\∞}^{\∞}\mathrm{\Ω}\left(f\right)\exp \left(i2\mathrm{\πft}\right)\mathrm{dt},$ Ω(f)＝A_θ,0(f)exp{i[Φ_θ,d(f)+Ψ(f)]}；

其中，i表示虚数单位，i²=-1。

207、根据瞬时幅值a₀(t)与随机相位过程θ(t)获得地震动随机过程x(t)。

具体采用以下公式获得地震动随机过程x(t)：

x(t)＝a₀(t)cosθ(t)。

208、根据地震动随机过程x(t)获得地震动样本。

为随机相位过程设定不同的随机数，利用伪随机数生成算法即可生成与所述给定随机数相对应的随机频率函数Θ(f)的一个样本，进而产生与所述随机数对应的地震动随机过程的一个样本，进而针对不同随机数可得到一系列的用于需要进行抗震设计的工程结构的输入地震动样本。

本发明实施例提供的获得地震动样本的方法，通过获取预先存储的天然地震动记录的瞬时幅值和瞬时相位，进而构造随机相位过程，瞬时幅值能够反映天然地震动记录的强度非平稳特性，而随机相位过程则能够反映天然地震动记录的频率非平稳特性，最后根据所述瞬时幅值与所述随机相位过程获得地震动随机过程，进而获得地震动样本，即获得与原始的天然地震动记录具有相同非平稳特性的一系列地震动样本。通过本发明实施例提供的方法获得的地震动样本可作为输入应用于利用蒙特卡洛法对结构地震反应的数值模拟研究工作中，以使得工程结构能够更加合理的抵御可能发生的地震作用。

与上述方法相对应的，本发明还提供了一种获得地震动样本的装置，如图3所示，包括：

第一获取单元31，用于从预先存储的天然地震动记录获取所述天然地震动记录的瞬时幅值和瞬时相位；

解卷单元32，用于对第一获取单元31获取的瞬时相位进行解卷；

第二获取单元33，用于采用傅立叶变换获取解卷单元32解卷后的瞬时相位的幅值谱和相位谱；

第一构造单元34，用于采用第二获取单元33获取的相位谱的高频分量构造随机频率过程；

第二构造单元35，用于采用第一构造单元34构造的随机频率过程、第二获取单元33获取的幅值谱及第二获取单元33获取的相位谱的低频分量构造随机相位过程；

第一获得单元36，用于根据第一获取单元31获取的瞬时幅值与第二构造单元35构造的随机相位过程获得地震动随机过程；

第二获得单元37，用于根据第一获得单元36获得地震动样本。

本发明实施例提供的获得地震动样本的装置，第一获取单元31获取预先存储的天然地震动记录的瞬时幅值和瞬时相位，进而第二构造单元35将构造随机相位过程，瞬时幅值能够反映天然地震动记录的强度非平稳特性，而随机相位过程则能够反映天然地震动记录的频率非平稳特性，最后第一获得单元36根据所述瞬时幅值与所述随机相位过程获得地震动随机过程，进而第二单元37获得地震动样本，即获得与原始的天然地震动记录具有相同非平稳特性的一系列地震动样本。

需要说明的是，本发明实施例提供的装置的各单元是与本发明实施例提供的方法对应设置的。

通过以上的实施方式的描述，所属领域的技术人员可以清楚地了解到本发明可借助软件加必需的通用硬件的方式来实现，当然也可以通过硬件，但很多情况下前者是更佳的实施方式。基于这样的理解，本发明的技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来，该计算机软件产品存储在可读取的存储介质中，如计算机的软盘，硬盘或光盘等，包括若干指令用以使得一台计算机设备（可以是个人计算机，服务器，或者网络设备等）执行本发明各个实施例所述的方法。

以上所述，仅为本发明的具体实施方式，但本发明的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本发明揭露的技术范围内，可轻易想到变化或替换，都应涵盖在本发明的保护范围之内。因此，本发明的保护范围应以所述权利要求的保护范围为准。

Claims (11)

1.一种获得地震动样本的方法，其特征在于，包括：

获取预先存储的天然地震动记录的瞬时幅值和瞬时相位；

对所述瞬时相位进行解卷；

采用傅立叶变换获取所述解卷后的瞬时相位的幅值谱和相位谱；

采用所述相位谱的高频分量构造随机频率过程；

采用所述随机频率过程、所述幅值谱及所述相位谱的低频分量构造随机相位过程；

根据所述瞬时幅值与所述随机相位过程获得地震动随机过程；

根据所述地震动随机过程获得地震动样本。

2.根据权利要求1所述的获得地震动样本的方法，其特征在于，所述获取预先存储的天然地震动记录的瞬时幅值和瞬时相位包括：

采用希尔伯特变换获取预先存储的天然地震动记录的瞬时幅值和瞬时相位。

3.根据权利要求1所述的获得地震动样本的方法，其特征在于，在所述采用所述相位谱的高频分量构造随机频率过程前，还包括：

采用汉宁窗函数对所述相位谱分解，得到相位谱的高频分量和低频分量。

4.根据权利要求2所述的获得地震动样本的方法，其特征在于，所述采用希尔伯特变换从预先存储的天然地震动记录获取所述天然地震动记录的瞬时幅值和瞬时相位包括：

使用下列公式获取所述天然地震动记录的瞬时幅值和瞬时相位，

$a_0\left(t\right)=\sqrt{x_0^2\left(t\right)+y_0^2\left(t\right)},$ ${\mathrm{\θ}}_0^{\′}\left(t\right)=\arctan \frac{y_0\left(t\right)}{x_0\left(t\right)},$ $y_0\left(t\right)=H\left[x_0\left(t\right)\right]=\frac{1}{\mathrm{\π}}P{\∫}_{-\∞}^{\∞}\frac{x_0\left(\mathrm{\τ}\right)}{t-\mathrm{\τ}}\mathrm{d\τ};$

其中，a₀(t)表示瞬时幅值，θ′₀(t)表示瞬时相位，x₀(t)表示预先存储的天然地震动记录，H表示希尔伯特变换，P表示取广义积分的柯西主值。

5.根据权利要求1所述的获得地震动样本的方法，其特征在于，所述采用傅立叶变换获取所述解卷后的瞬时相位的幅值谱和相位谱包括：

使用下列公式获取所述解卷后的瞬时相位的幅值谱和相位谱，

$A_{\mathrm{\θ},0}\left(f\right)=\sqrt{{\mathrm{Re}}^2\left[{\mathrm{\Ω}}_0\left(f\right)\right]+{\mathrm{Im}}^2\left[{\mathrm{\Ω}}_0\left(f\right)\right]},$ ${\mathrm{\Φ}}_{\mathrm{\θ},0}\left(f\right)=\arctan \frac{\mathrm{Im}\left[{\mathrm{\Ω}}_0\left(f\right)\right]}{\mathrm{Re}\left[{\mathrm{\Ω}}_0\left(f\right)\right]},$

${\mathrm{\Ω}}_0\left(f\right)={\∫}_{-\∞}^{\∞}{\mathrm{\θ}}_0\left(t\right)\exp (-i2\mathrm{\πft})\mathrm{dt};$

其中，A_θ,0(f)表示解卷后的瞬时相位的幅值谱，Φ_θ,0(f)表示解卷后的瞬时相位的相位谱，θ₀(t)表示解卷后的瞬时相位，Re表示取复数的实部，Im表示取复数的虚部，i表示虚数单位，i²=-1。

6.根据权利要求3所述的获得地震动样本的方法，其特征在于，所述采用汉宁窗函数对所述相位谱分解，得到相位谱的高频分量和低频分量包括：

使用下列公式得到相位谱的高频分量和低频分量，

Φ_θ,d(f)＝Hann[Φ_θ,0(f)]，Φ_θ，n(f)＝Φ_θ,0(f)-Φ_θ,d(f)；

其中，Φ_θ,d(f)表示相位谱的低频分量，Φ_θ，n(f)表示相位谱的高频分量，Φ_θ,0(f)表示解卷后的瞬时相位的相位谱，Hann表示利用汉宁窗函数对信号进行低通滤波处理。

7.根据权利要求1所述的获得地震动样本的方法，其特征在于，所述采用所述相位谱的高频分量构造随机频率过程包括：

使用下列公式构造随机频率过程，

Ψ(f)＝A(f)Θ(f)， $A\left(f\right)=\sqrt{{\mathrm{\Φ}}_{\mathrm{\θ},n}^2\left(f\right)+{\widehat{\mathrm{\Φ}}}_{\mathrm{\θ},n}^2\left(f\right)},$ ${\widehat{\mathrm{\Φ}}}_{\mathrm{\θ},n}\left(f\right)=H\left[{\mathrm{\Φ}}_{\mathrm{\θ},n}\left(f\right)\right]=\frac{1}{\mathrm{\π}}P{\∫}_{-\∞}^{\∞}\frac{{\mathrm{\Φ}}_{\mathrm{\θ},n}\left(\mathrm{\τ}\right)}{f-\mathrm{\τ}}\mathrm{d\τ};$

其中，Ψ(f)表示随机频率过程，Φ_θ,n(f)表示相位谱的高频分量，Θ(f)表示随机频率函数，H表示希尔伯特变换，P表示取广义积分的柯西主值。

8.根据权利要求7所述的获得地震动样本的方法，其特征在于，所述随机频率函数Θ(f)为在区间[-1,1]均匀分布的随机频率函数，并满足：

对任意频率f，Θ(f)为区间[-1,1]均匀分布的随机变量；及

若f₁≠f₂，则Θ(f₁)与Θ(f₂)为相互独立的随机变量。

9.根据权利要求1所述的获得地震动样本的方法，其特征在于，所述采用所述随机频率过程、所述幅值谱及所述相位谱的低频分量构造随机相位过程包括：

使用下列公式构造随机相位过程，

$\mathrm{\θ}\left(t\right)={\∫}_{-\∞}^{\∞}\mathrm{\Ω}\left(f\right)\exp \left(i2\mathrm{\πft}\right)\mathrm{dt},$ Ω(f)＝A_θ,0(f)exp{i[Φ_θ,d(f)+Ψ(f)]}；

其中，θ(t)表示随机相位过程，A_θ,0(f)表示解卷后的瞬时相位的幅值谱，Φ_θ,d(f)表示相位谱的低频分量，Ψ(f)表示随机频率过程，i表示虚数单位，i²=-1。

10.根据权利要求1所述的获得地震动样本的方法，其特征在于，所述根据所述瞬时幅值与所述随机相位过程获得地震动随机过程包括：

使用下列公式构造地震动随机过程：

x(t)＝a₀(t)cosθ(t)；

其中，x(t)表示地震动随机过程，a₀(t)表示瞬时幅值，θ(t)表示随机相位过程。

11.一种获得地震动样本的装置，其特征在于，包括：

第一获取单元，用于从预先存储的天然地震动记录获取所述天然地震动记录的瞬时幅值和瞬时相位；

解卷单元，用于对所述第一获取单元获取的瞬时相位进行解卷；

第二获取单元，用于采用傅立叶变换获取所述解卷单元解卷后的瞬时相位的幅值谱和相位谱；

第一构造单元，用于采用所述第二获取单元获取的相位谱的高频分量构造随机频率过程；

第二构造单元，用于采用所述第一构造单元构造的随机频率过程、所述第二获取单元获取的幅值谱及所述第二获取单元获取的相位谱的低频分量构造随机相位过程；

第一获得单元，用于根据所述第一获取单元获取的瞬时幅值与所述第二构造单元构造的随机相位过程获得地震动随机过程；

第二获得单元，用于根据所述地震动随机过程获得地震动样本。

Images