河西走廊西段及邻区主要断裂（一）-晚新生代逆断层与走滑断层的地震剖面解释
-
摘要: 基于地震剖面解释的构造分析表明，河西走廊西段及邻区发育有许多晚第三纪-第四纪的逆断层与走滑断层. 其主要类型有两种：一是与北祁连逆冲推覆作用有关的南倾逆断层，其走向北西，主要分布于北祁连山东北缘及河西走廊西南部地区；二是与阿尔金断裂活动有关的逆断层与走滑断层，主要分布于红柳峡、宽台山——合黎山地区. 这些断层与新生代始新世末期以来印-藏碰撞的远距离效应有关. 它们于晚第三纪开始活动，主要活动于第四纪时期，而且大多目前仍在活动，属于活动性断层，是河西走廊地区地震的潜发因素. 此外，它们还造成了地层的强烈变形，形成了一系列的地貌阶地和河流的水平错移.
-
引言
地震动是空间变化的，由于波列传播速度的有限性和相干性损失以及局部场地地质条件不同，结构各支承点的地震激励出现显著差异，这种差异可能对平面大尺度结构（例如桥梁、水坝）的地震反应产生重要影响。欧洲规范8 （European Committee for Standardization，2004）和城市轨道交通结构抗震设计规范GB 50909—2014 （中华人民共和国住房和城乡建设部，2014）也提到了地震动空间变化对结构的影响。地震动空间变化可以通过以下三种方法来分析：① 时间历程方法；② 随机振动方法；③ 反应谱方法。后两种方法的优点是应用了地震动的统计特征对结构响应进行分析。
对于随机振动分析，地震动输入由自功率谱密度矩阵（auto-power spectral density，缩写为APSD）确定；对于反应谱分析，地震动输入由地震反应谱确定，模态组合系数通常由结构随机振动理论推导得出（王君杰，1992；Berrah，Kausel，1992；Kiureghian，Neuenhofer，1992；王君杰等，1995）。一般情况下，随机振动法的数值积分运算量很大，基于随机振动理论的反应谱组合系数的计算也非常耗时。为减少积分时间，Harichandran和Vanmarcke （1986）提出了行波激励的解析公式，Loh和Ku （1995）基于简化的自功率谱密度函数提出了一种有效的计算方法，但采用了简单的空间相干函数。孙建梅等（2003）对Kiureghianh和Neuenhofer （1992）提出的多点输入反应谱方法进行了简化，即在相关系数的计算中忽略频率比介于0.7—1.2范围以外的部分，减少了计算量，但损失了精度。
为提高多点地震动激励下谱分析方法中相关系数的计算效率，本文拟基于对自功率谱和相关系数积分表达式的近似处理，提出相关系数的近似解析公式。
1. 结构地震响应
假定在空间某点r （
$r ＝ 1， 2， \cdots ， m$ ）处的地震激励是零均值平稳过程，其位移、速度和加速度函数分别记为${u_r} （ t ） $ ，${\dot u_r} （ t ） $ 和${\ddot u_r} （ t ） $ ，结构平稳响应的方差为（王君杰，1992；Kiureghian，Neuenhofer，1992）$$\begin{split} \sigma _k^2 ＝& \sum\limits_{r ＝ 1}^m {\sum\limits_{s ＝ 1}^m {{a_{kr}}{a_{ks}}} } {\rho _{{\rm{g}}r{\rm{g}}s}}{\sigma _{{\rm{g}}r}}{\sigma _{{\rm{g}}s}} ＋ 2\sum\limits_{r ＝ 1}^m {\sum\limits_{s ＝ 1}^m {\sum\limits_{j ＝ 1}^n {{a_{kr}}} } } {b_{kjs}}{\rho _{{\rm{g}}rjs}}{\sigma _{{\rm{g}}r}}{\sigma _{js}} ＋\\&\qquad\qquad \sum\limits_{r ＝ 1}^m {\sum\limits_{s ＝ 1}^m {\sum\limits_{i ＝ 1}^n {\sum\limits_{j ＝ 1}^n {{b_{kir}}} } } } {b_{kjs}}{\rho _{irjs}}{\sigma _{ir}}{\sigma _{js}} \text{，} \end{split}$$ (1) 式中：σk为第k个结构反应量的均方差；σgr和σgs分别为第r个空间点和第s个空间点的地面运动位移的均方差；σir和σjs分别为第r个地面运动激励下结构第i振型反应的均方差和第s个地面运动激励下结构第j振型反应的均方差；参数a和b均为与结构有关的系数；m为多点地震动的数目；n为参与计算的结构的振型数目；ρgrgs，ρgrjs和ρirjs定义为
$$ \left\{ \begin{array}{l}{\rho_{{\rm{g}}r{\rm{g}}s}} ＝ \dfrac{{{I_{{\rm{g}}r{\rm{g}}s}}}}{{\sqrt {{I_{{\rm{g}}r{\rm{g}}r}} {I_{{\rm{g}}s{\rm{g}}s}}} }} ， \\{\rho _{{\rm{g}}rjs}} ＝ \dfrac{{{I_{{\rm{g}}rjs}}}}{{\sqrt {{I_{{\rm{g}}r{\rm{g}}r}} {I_{jsjs}}} }} ，\\{\rho _{irjs}} ＝ \dfrac{{{I_{irjs}}}}{{\sqrt {{I_{irir}} {I_{jsjs}}} }} \text{，} \end{array}\right.$$ (2) 其中
$$ \left\{ \begin{array}{l}{I_{{\rm{g}}r{\rm{g}}s}} ＝ {{\rm{Re}}} \Bigg[{\displaystyle\int\nolimits_{ - \infty }^{ ＋ \infty } {\dfrac{1}{{{\omega ^4}}}S_{{\rm{g}}r{\rm{g}}s}^{\rm{A}} （ \omega ） {\rm{d}}\omega } } \Bigg] \text{，} \\ {I_{{\rm{g}}rjs}} ＝ {{\rm{Re}}} \Bigg[ { - \displaystyle\int\nolimits_{ - \infty }^{ ＋ \infty } {\dfrac{1}{{{\omega ^2}}}{H_j} （ \omega ） S_{{\rm{g}}r{\rm{g}}s}^{\rm{A}} （ \omega ） {\rm{d}}\omega } } \Bigg] \text{，} \\ {I_{irjs}} ＝ {{\rm{Re}}} \Bigg[ {\displaystyle\int\nolimits_{ - \infty }^{ ＋ \infty } {H_i^* （ \omega ） {H_j} （ \omega ） S_{{\rm{g}}r{\rm{g}}s}^{\rm{A}} （ \omega ） {\rm{d}}\omega } } \Bigg] \text{，} \end{array}\right.$$ (3) 式中，
${H_i} （ \omega ） $ 为结构自由振动第i阶振型的频响函数，$S_{{\rm{g}}r{\rm{g}}s}^{\rm{A}} （ \omega ） $ 为第r空间点地震加速度与第s空间点地震加速度的互功率谱密度函数，ω为圆频率。根据式（1），多点地震激励下结构反应的组合方法为$$ \begin{split}E[\max |{{\textit{z}}_k} （ t ） |] ＝& \Bigg[\sum\limits_{r ＝ 1}^m {\sum\limits_{s ＝ 1}^m {{a_{kr}}{a_{ks}}{\rho _{{\rm{g}}r{\rm{g}}s}}{u_{r， \max }}{u_{s， \max }}} } ＋ 2\sum\limits_{r ＝ 1}^m {\sum\limits_{s ＝ 1}^m {\sum\limits_{j ＝ 1}^n {{a_{kr}}} } } {b_{ksj}}{\rho _{{\rm{g}}rjs}}{u_{r， \max }} S {D_s} （ {\omega _j}， {\zeta _j} ） ＋ \\& \sum\limits_{r ＝ 1}^m {\sum\limits_{s ＝ 1}^m {\sum\limits_{i ＝ 1}^n {\sum\limits_{j ＝ 1}^n {{b_{kri}}} } } } {b_{ksj}}{\rho _{irjs}}S {D_r} （ {\omega _i}， {\zeta _i} ） S {D_s} （ {\omega _j}， {\zeta _j} ） {\Bigg]^{\tfrac{1}{2}}} ，\\[-20pt] \end{split}$$ (4) 式中，
$S {D_r} （ {\omega _i}，{\zeta _i} ） $ 和$S {D_s} （ {\omega _j}，{\zeta _j} ） $ 分别代表第r个地震激励下i阶振型最大位移和第s个地震激励下j阶振型最大位移，即平均位移反应谱。同时，${u_{r， \max }} ＝ S {D_r} （ 0， {\zeta _i} ） $ ，${u_{s， \max }} ＝ S {D_s} （ 0， {\zeta _j} ） $ 。对于式（3），根据复相干函数，互功率谱密度函数
${S_{rs}} （ \omega ） $ 可以写为$$ {\gamma _{rs}} （ {\rm{i}}\omega ） ＝ \frac{{{S_{rs}} （ {\rm{i}}\omega ） }}{{\sqrt {{S_{rr}} （ \omega ） {S_{ss}} （ \omega ） } }} ＝ {\rho _{rs}} （ \omega ， \Delta {{\boldsymbol{r}}_{rs}} ） \exp [ - {\rm{i}}{\theta _{rs}} （ \omega ， \Delta {{\boldsymbol{r}}_{rs}} ） ] {， } $$ (5) 式中，γ为相关因子，
$\;{\rho _{rs}} （ \omega ， \Delta {{\boldsymbol{r}}_{rs}} ） $ 为地震动的空间相干函数，${\theta _{rs}} （ \omega ， \Delta {{\boldsymbol{r}}_{rs}} ） $ 为地震动的相位差，$\Delta {{\boldsymbol{r}}_{rs}}$ 为空间距离的向量。国内外研究人员已经提出了多个
${\;\rho _{rs}} （ \omega ， \Delta {{\boldsymbol{r}}_{rs}} ） $ 表达式，例如Harichandran和Vanmarcke （1986）及Qu等（1996）提出的表达式等。Qu等（1996）提出的地震动空间相干函数模型如图1所示，其经验表达式为![]()
图 1 空间相干函数模型（Qu et al，1996）（a） ρrs-空间距离关系；（b） ρrs-频率关系Figure 1. Spatial coherence function model （Qu et al，1996）（a） ρrs-spatial distance relationship；（b） ρrs-frequency relationship$${\rho _{rs}} （ \omega ， \Delta {{\boldsymbol{r}}_{rs}} ） ＝ \exp ［ { - a （ \omega ） {{\left| {\Delta {{\boldsymbol{r}}_{rs}}} \right|}^{b （ \omega ） }}} ］ {， } a （ \omega ） ＝ {a_1}{\omega ^2} ＋ {a_2}， b （ \omega ） ＝ {b_1}\omega ＋ {b_2}， $$ (6) 式中：
${a_1} ＝ 1.678 {\text{×}} {10^{ - 5}}$ ，${a_2} ＝ 1.219 {\text{×}} {10^{ - 3}}$ ；${b_1} ＝ - 5.5 {\text{×}} {10^{ - 3}}$ ，${b_2} ＝ 0.7674$ 。Harichandran和Vanmarcke（1986）提出的地震动空间相干函数模型如图2所示，其经验表达式为
![]()
图 2 空间相干函数模型（Harichandran，Vanmarcke，1986）（a） ρrs-空间距离关系；（b） ρrs-频率关系Figure 2. Spatial coherence function model （Harichandran，Vanmarcke，1986）（a） ρrs-spatial distance relationship；（b） ρrs-frequency relationship$$ \left\{ \begin{array}{l}{\rho _{rs}} （ {\omega ， \Delta {{\boldsymbol{r}}_{rs}}} ） ＝ A\exp \left[ { - \dfrac{{2\left| {\Delta {{\boldsymbol{r}}_{rs}}} \right|B}}{{\alpha \theta （ \omega ） }}} \right] ＋ （ {1 - A} ） \exp \left[ { - \dfrac{{2\left| {\Delta {{\boldsymbol{r}}_{rs}}} \right|B}}{{\theta （ \omega ） }}} \right] ， \\ B ＝ 1 － A ＋ \alpha A，\\ \theta （ \omega ） ＝ \dfrac{K}{{\sqrt {1{\rm{ + }}{{\left( {\dfrac{\omega }{{{\omega _0}}}} \right)}^b}} }}， \end{array}\right.$$ (7) 式中
$A$ ，$\alpha $ ，$K$ ，$\theta $ ，${\omega _0}$ 和$b$ 为经验参数。对于台湾SMART-1台网第24次地震，$A ＝ 0.736$ ，$\alpha ＝ 0.147$ ，$K ＝ 5\;210$ ，${\omega _0} ＝ 6.85$ rad/s，$b ＝ 2.78$ 。本文中使用到的两个自功率谱函数经验模型见图3。第一个是Hu模型（胡聿贤，周锡元，1962，1965；王君杰，1992；王君杰，江近仁，1997），即
![]()
图 3 自功率谱密度函数${S_{rr}} （ \omega ） $ 的特征（ζg为场地土的阻尼比，下同）（a） 胡聿贤自功率谱密度模型；（b） 克拉夫−彭津自功率谱密度模型Figure 3. The characteristic for auto-power spectral density function （ζg is the damping ratio of site soil，the same below）（a） Hu’s auto-power spectral density model；（b） Clough-Penzien’s auto-power spectral density model$$ {S _{rr}} （ \omega ） ＝ {S _0}\frac{{\omega _{\rm{g}}^4 ＋ 4\zeta _{\rm{g}}^2\omega _{\rm{g}}^2{\omega ^2}}}{{{{ （ \omega _{\rm{g}}^2 － {\omega ^2} ） }^2} ＋ 4\zeta _{\rm{g}}^2\omega _{\rm{g}}^2{\omega ^2}}}\frac{{{\omega ^4}}}{{{\omega ^4} ＋ \omega _{\rm{c}}^4}}\text{；} $$ (8) 第二个是克拉夫−彭津模型（Clough，Penzien，1993），即
$$ {S _{rr}} （ \omega ） ＝ {S _0}\frac{{\omega _{\rm{g}}^4 ＋ 4\zeta _{\rm{g}}^2\omega _{\rm{g}}^2{\omega ^2}}}{{{{ （ \omega _{\rm{g}}^2 － {\omega ^2} ） }^2} ＋ 4\zeta _{\rm{g}}^2\omega _{\rm{g}}^2{\omega ^2}}}\frac{{{\omega ^4}}}{{{{ （ \omega _{\rm{f}}^2 － {\omega ^2} ） }^2} ＋ 4\zeta _{\rm{f}}^2\omega _{\rm{f}}^2{\omega ^2}}}{\text{.}} $$ (9) 本文在计算相关系数时，将以上两个模型分别简化为
$$ {S _{rr}} （ \omega ） ＝ {S _0}\frac{{{\omega ^4}}}{{{\omega ^4} ＋ \omega _{\rm{c}}^4}} \;\;\;\;\qquad\qquad 胡聿贤模型 $$ (10) $$ {S _{rr}} （ \omega ） ＝ {S _0}\frac{{{\omega ^4}}}{{{{ （ \omega _{\rm{f}}^2 － {\omega ^2} ） }^2} ＋ 4\zeta _{\rm{f}}^2\omega _{\rm{f}}^2{\omega ^2}}} \qquad 克拉夫-彭津模型 $$ (11) 式中，
${S_0}$ 为APSD的强度因子，${\omega _{\rm{g}}}$ 为场地土层的频率参数，${\zeta _{\rm{g}}}$ 为场地土层的阻尼参数，${\omega _{\rm{c}}}$ 为胡聿贤模型在低频段的滤波参数，${\omega _{\rm{f}}}$ 和${\zeta _{\rm{f}}}$ 为Clough-Penzien模型在低频段的滤波参数。计算中取${\omega _{\rm{g}}} ＝$ 1.5 Hz，${\omega _{\rm{c}}} ＝ $ 0.3 Hz，${\omega _{\rm{f}}} ＝ $ 0.25 Hz，${\zeta _{\rm{g}}} ＝$ 0.6，${\zeta _{\rm{f}}} ＝ $ 0.4，${S_0}$ ＝1。2. 近似自功率谱密度函数的精度
分别比较相干系数
$\;{\rho _{{\rm{g}}r{\rm{g}}s}}$ ，$\; {\rho _{{\rm{g}}rjs}} $ 和$\; {\rho _{irjs}} $ 的近似值 ［ 式（10） ］ 与精确值 ［ 式（8） ］ 之间的差异，结果如图4—6所示。其中，空间相干函数使用式（6），视波速向量$ {{\boldsymbol{v}}_{{\rm{app}}}} $ 的模取1 000 m/s。![]()
图 4 系数$ \;{\rho _{{\rm{g}}r{\text{g}}s}} $ 的精确值与近似值的比较Figure 4. The comparison between exact value and approximate value of coefficient$ \;{\rho _{{\rm{g}}r{\text{g}}s}} $ ![]()
图 6 系数$ \;{\rho _{irjs}} $ （$ {\zeta _i} = {\zeta_j} = 0.05 $ ）精确值与近似值的比较Figure 6. The comparison between exact value and approximate value of coefficient$\; {\rho _{irjs}} $ （$ {\zeta _i} = {\zeta_j} = 0.05 $ ）（a） $ \left| {\Delta {{\boldsymbol{r}}_{rs}}} \right| ＝ 300\;{\rm{m}} $，ωi＝0.1 Hz；（b） $ \left| {\Delta {{\boldsymbol{r}}_{rs}}} \right| ＝ 1\;000\;{\rm{m}} $，ωi＝0.1 Hz；（c） $ \left| {\Delta {{\boldsymbol{r}}_{rs}}} \right| ＝ 300\;{\rm{m}} $，ωi＝1.0 Hz； （d） $ \left| {\Delta {{\boldsymbol{r}}_{rs}}} \right| ＝ 1\;000\;{\rm{m}} $，ωi＝1.0 Hz；（e） $ \left| {\Delta {{\boldsymbol{r}}_{rs}}} \right| ＝ 300\;{\rm{m}} $，ωi＝4.0 Hz；（f） $ \left| {\Delta {{\boldsymbol{r}}_{rs}}} \right| ＝ 1\;000\;{\rm{m}} $，ωi＝4.0 Hz从图4可以看出，式（10）给出的系数
$\;{\rho _{{\rm{g}}r{\rm{g}}s}}$ 几乎是精确值。图5表明：当${\omega _j}$ 处于低频段时，系数$\; {\rho _{{\rm{g}}rjs}} $ 非常大；随着${\omega _j}$ 的增大及r点与s点之间空间距离的增大，$ \;{\rho _{{\rm{g}}rjs}} $ 逐渐减小至零，直至可以忽略，说明式（10）可以给出具有高精度的$ \;{\rho _{{\rm{g}}rjs}} $ 近似解答。![]()
图 5$\left| {\Delta {{\boldsymbol{r}}_{rs}}} \right| ＝ 300\;{\rm{m}}$ （a）和1 000 m （b）时系数$\; {\rho _{{\rm{g}}rjs}} $ 的精确值与近似值的比较（${\zeta _j} ＝ 0.05$ ）Figure 5. The comparison between exact value and approximate value of coefficient$\; {\rho _{{\rm{g}}rjs}} $ （${\zeta _j} ＝ 0.05$ ） with$ \left| {\Delta {{\boldsymbol{r}}_{rs}}} \right| ＝ 300\;{\rm{m}} $ （a） and 1 000 m （b）从图6a-d可以看出：当
${\omega _i}$ 和${\omega _j}$ 位于低频段时，由式（10）计算得出的相干系数$ \;{\rho _{irjs}} $ 近乎精确值，而近似的$ \;{\rho _{irjs}} $ 在${\omega _i}$ 和${\omega _j}$ 较大时会出现更大的误差；在高频段内，地震动自功率谱的值比低频段的值小得多，相比之下重要性不大。同时，系数$ \;{\rho _{irjs}} $ 的值随着${\omega _i}$ 和${\omega _j}$ 的增加而减小。以上结果表明，相干系数$ \;{\rho _{irjs}} $ 的近似值在工程实践中具有很好的准确性。3. 相关系数的近似公式
3.1 近似空间相干函数
$\;{\rho _{rs}} （ \omega ， \Delta {{\boldsymbol{r}}_{rs}} ） $ 如前所述，由于空间相干函数的经验表达式
$\; {\rho _{rs}} （ \omega ， \Delta {{\boldsymbol{r}}_{rs}} ） $ 的复杂性，需要进行数值积分来计算系数$\;{\rho _{{\rm{g}}r{\rm{g}}s}}$ ，$ \;{\rho _{{\rm{g}}rjs}} $ 和$\; {\rho _{irjs}} $ 。为了得到这三个系数的近似解析表达式，引入近似的$\;{\rho _{rs}} （ \omega ， {\Delta {\boldsymbol{r}}_{{{r}}s}} ） $ 。组合系数
$\;{\rho _{{\rm{g}}r{\rm{g}}s}}$ 与拟静态响应有关。拟静态响应是由位移贡献的，它由地震动中非常低的频率成分控制。因此，计算系数$\;{\rho _{{\rm{g}}r{\rm{g}}s}}$ 的相干函数近似表达式为$$ \begin{split}\\{\rho _{rs}} （ \omega ， \Delta {{\boldsymbol{r}}_{rs}} ） ＝ {\rho _{rs}} （ {\omega _{d， \max }}， \Delta {{\boldsymbol{r}}_{rs}} ） {， } \end{split}$$ (12) 式中，
${\omega _{d， \max }}$ 是地震动中位移的APSD达到最大值时对应的圆频率。组合系数
$ \;{\rho _{{\rm{g}}rjs}} $ 与式（1）或式（4）中的第二项耦合响应项有关，此耦合项由地震动中低频成分控制。因此，计算系数$ \;{\rho _{{\rm{g}}rjs}} $ 的相干函数近似表达式被定义为$$ {\rho _{rs}} （ \omega ， \Delta {{\boldsymbol{r}}_{rs}} ） ＝ {\rho _{rs}} （ {\omega _{{\rm{g}}rjs}}， \Delta {{\boldsymbol{r}}_{rs}} ） {， } $$ (13) 式中，
$ {\omega _{{\rm{g}}rjs}} $ 取$\sqrt {{\omega _j}{\omega _{d， \max }}}$ 或$（ {{\omega _j} ＋ {\omega _{d， \max }}} ） /2$ 。组合系数
$ \;{\rho _{irjs}} $ 与式（1）或式（4）中的第三项动力响应项有关，动力响应项由地震动中靠近以及介于${\omega _i}$ 和${\omega _j}$ 之间的频率成分控制。因此，计算系数$ \;{\rho _{irjs}} $ 的相干函数近似表达式为$$ {\rho _{rs}} （ \omega ， \Delta {{\boldsymbol{r}}_{rs}} ） ＝ {\rho _{rs}} （ {\omega _{irjs}}， \Delta {{\boldsymbol{r}}_{rs}} ） {， } $$ (14) 式中，
$ {\omega _{irjs}} $ 取$\sqrt {{\omega _i}{\omega _j}} $ 或$ （ {{\omega _i} + {\omega _j}} ） /2$ 。从图7—9所示结果可以看出，使用式（6）和式（8）计算得出的上述近似相干函数
$\;{\rho _{rs}} （ \omega ， {d_{rs}} ） $ 具有很好的准确性。由此可知，本文所提方法在$\; {\rho _{rs}} （ \omega ， \Delta {{\boldsymbol{r}}_{rs}} ） $ 的近似表达上具有良好的精度。![]()
图 7 系数$\;{\rho _{{\rm{g}}r{\rm{g}}s}}$ 的精确值与近似值的比较Figure 7. The comparison between exact value and approximate value of coefficient$\;{\rho _{{\rm{g}}r{\rm{g}}s}}$ ![]()
图 8$ \left| {\Delta {{\boldsymbol{r}}_{rs}}} \right| ＝ 300\;{\rm{m}} $ （a）和1 000 m （b）时系数$\; {\rho _{{\rm{g}}rjs}} $ （${\zeta _j} ＝ 0.05$ ）精确值与近似值的比较Figure 8. The comparison between exact value and approximate value of coefficient$ \;{\rho _{{\rm{g}}rjs}} $ （${\zeta _j} ＝ 0.05$ ） with$ \left| {\Delta {{\boldsymbol{r}}_{rs}}} \right| ＝ 300\;{\rm{m}} $ （a） and 1 000 m （b）![]()
图 9 系数$ \;{\rho _{irjs}} $ （${\zeta _i} ＝ {\zeta _j} ＝ 0.05$ ）精确值与近似值的比较Figure 9. The comparison between exact value and approximate value of coefficient$\; {\rho _{irjs}} $ （${\zeta _i} ＝ {\zeta _j} ＝ 0.05$ ）（a） $ \left| {\Delta {{\boldsymbol{r}}_{rs}}} \right| ＝ 300\;{\rm{m}} $，ωi＝0.1 Hz；（b） $ \left| {\Delta {{\boldsymbol{r}}_{rs}}} \right| ＝ 1\;000\;{\rm{m}} $，ωi＝0.1 Hz；（c） $\left| {\Delta {{\boldsymbol{r}}_{rs}}} \right| ＝ 300\;{\rm{m}}$，ωi＝1.0 Hz； （d） $ \left| {\Delta {{\boldsymbol{r}}_{rs}}} \right| ＝ 1\;000\;{\rm{m}} $，ωi＝1.0 Hz；（e） $\left| {\Delta {{\boldsymbol{r}}_{rs}}} \right| ＝ 300\;{\rm{m}}$，ωi＝4.0 Hz；（f） $ \left| {\Delta {{\boldsymbol{r}}_{rs}}} \right| ＝ 1\;000\;{\rm{m}} $，ωi＝4.0 Hz经过上述简化后，三个相关系数
$\;{\rho _{{\rm{g}}r{\rm{g}}s}}$ ，$\; {\rho _{{\rm{g}}rjs}} $ 和$ \;{\rho _{irjs}} $ 可以简化为$$ \left\{ \begin{array}{l}{\rho _{{\rm{g}}r{\rm{g}}s}} ＝ \dfrac{1}{{\sqrt {{I_{{\rm{g}}r{\rm{g}}r}}{I_{{\rm{g}}s{\rm{g}}s}}} }}{\rho _{rs}} （ {\omega _{d， \max }}， \Delta {{\boldsymbol{r}}_{rs}} ） {I_{{\rm{g}}r{\rm{g}}s}} ，\\ {\rho _{{\rm{g}}rjs}} ＝ \dfrac{1}{{\sqrt {{I_{{\rm{g}}r{\rm{g}}r}}{I_{jsjs}}} }}{\rho _{rs}} （ {\omega _{{\rm{g}}rjs}}， \Delta {{\boldsymbol{r}}_{rs}} ） {I_{{\rm{g}}rjs}} ，\\ {\rho _{irjs}} ＝ \dfrac{1}{{\sqrt {{I_{irir}}{I_{jsjs}}} }}{\rho _{rs}} （ {\omega _{irjs}}， \Delta {{\boldsymbol{r}}_{rs}} ） {I_{irjs}} . \end{array}\right.$$ (15) 对于胡聿贤自功率谱密度模型有
$$ \left\{ \begin{array}{l}{I_{{\rm{g}}r{\rm{g}}s}} ＝ {{\rm{Re}}} \Bigg({\displaystyle\int_{ - \infty }^\infty {\dfrac{1}{{{\omega ^4} ＋ \omega _{\rm{c}}^4}}{{\rm{e}}^{ - {\rm{i}}\omega \tau }}{\rm{d}}\omega } } \Bigg) { \text{，}} \\ {I_{{\rm{g}}rjs}} = {{\rm{Re}}}\Bigg[ { - \displaystyle\int_{ - \infty }^\infty {{H_j} （ {\rm{i}}\omega ） \dfrac{{{\omega ^2}}}{{{\omega ^4} ＋ \omega _{\rm{c}}^4}}{{\rm{e}}^{ - {\rm{i}}\omega \tau }}{\rm{d}}\omega } }\Bigg] {\text{，}} \\ {I_{irjs}} ＝ {{\rm{Re}}} \Bigg[{\displaystyle\int_{ - \infty }^\infty {H_{i}^* （ {\rm{i}}\omega ） {H_j} （ {\rm{i}}\omega ） \dfrac{{{\omega ^4}}}{{{\omega ^4} ＋ \omega _{\rm{c}}^4}}{{\rm{e}}^{ - {\rm{i}}\omega \tau }}{\rm{d}}\omega } }\Bigg] { \text{，}}\end{array} \right.$$ (16) 式中：
${\rm{Re}} $ 为复数的实部；$ \tau $ 为地震波从r点传播到s点的时间，$\tau ＝{\Delta {{\boldsymbol{r}}}_{rs}\cdot{{{\boldsymbol{v}}_{\rm{app}}}}}/{{\left|{{{\boldsymbol{v}}_{\rm{app}}}}\right|}^{2}}$ ，当r与s在相同位置时，$ \tau $ 等于零。对于克拉夫−彭津的自功率谱密度模型则有
$$\left\{ \begin{split}& {I_{{\rm{g}}r{\rm{g}}s}} ＝ {{\rm{Re}}} \Bigg[{\displaystyle\int_{ - \infty }^\infty {\dfrac{1}{{{{ （ \omega _{\rm{f}}^2 － {\omega ^2} ） }^2} ＋ 4\zeta _{\rm{f}}^2\omega _{\rm{f}}^2{\omega ^2}}}{{\rm{e}}^{ - {\rm{i}}\omega \tau }}{\rm{d}}\omega } } \Bigg] ， \\& {I_{{\rm{g}}rjs}} ＝ {{\rm{Re}}} \Bigg[ { - \displaystyle\int_{ - \infty }^\infty {{H_j} （ {\rm{i}}\omega ） \dfrac{{{\omega ^2}}}{{{{ （ \omega _{\rm{f}}^2 - {\omega ^2} ） }^2} ＋ 4\zeta _{\rm{f}}^2\omega _{\rm{f}}^2{\omega ^2}}}{{\rm{e}}^{ - {\rm{i}}\omega \tau }}{\rm{d}}\omega } } \Bigg] ， \\& {I_{irjs}} ＝ {{\rm{Re}}}\Bigg[ {\displaystyle\int_{ - \infty }^\infty {H_{i}^* （ {\rm{i}}\omega ） {H_j} （ {\rm{i}}\omega ） \dfrac{{{\omega ^4}}}{{{{ （ \omega _{\rm{f}}^2 － {\omega ^2} ） }^2} ＋ 4\zeta _{\rm{f}}^2\omega _{\rm{f}}^2{\omega ^2}}}{{\rm{e}}^{ - {\rm{i}}\omega \tau }}{\rm{d}}\omega } } \Bigg].\end{split} \right.\\[-15pt]$$ (17) 3.2 胡聿贤自功率谱密度模型下相关系数的解析表达式
根据复变函数的基本理论，式（16）的积分形式为
$$ I ＝ {{\rm{Re}}} \Bigg\{ { - 2\pi {\rm{i}} \sum\limits_{{{\rm{Im}}} {{\textit{z}}_{p}} ＜ 0} {{{\rm{Res}}} ［ f （ {{\textit{z}}_{p}} ） ］ \Bigg\}} } ＝ {{\rm{Re}}} \Bigg\{ { - 2\pi {\rm{i}} \sum\limits_{{{\rm{Im}}} {{\textit{z}}_{p}} ＜ 0} {\left[ {\frac{{B （ {{\textit{z}}_{p}} ） }}{{A （ {{\textit{z}}_{p}} ） }}} \right]} }\Bigg\} \text{，} $$ (18) 式中，
${{\rm{Res}}} ［ f （ {{\textit{z}}_{{p}}} ） ］$ 为复变函数的留数，$z_p$ 是函数f（z）的极点，Imzp是zp的虚部。根据式（18）得到
$$ {I_{{\rm{g}}r{\rm{g}}s}} ＝ {{\rm{Re}}} \Bigg\{ { - 2\pi {\rm{i}}\sum\limits_{p ＝ 1}^2 {\left[ {\frac{{B （ {{\textit{z}}_{{p}}} ） }}{{A （ {{\textit{z}}_{{p}}} ） }}} \right]} } \Bigg\} ；A （ {{\textit{z}}_{{p}}} ） ＝ \prod\limits_{q ＝ 1，q {\text{≠}} p }^4 { （ {{\textit{z}}_{{p}}} － {{\textit{z}}_q} ） } \text{，} B （ {{\textit{z}}_{{p}}} ） ＝{{\rm{e}}^{ - {\rm{i}}{{\textit{z}}_{{p}}}\left| \tau \right|}} \text{，} $$ (19) 其中，
${{\textit{z}}_1} = {\omega _{\rm{c}}}{{\rm{e}}^{{\rm{i}}5{\theta _{\rm{c}}}}}$ ，${{\textit{z}}_2} = {\omega _{\rm{c}}}{{\rm{e}}^{{\rm{i}}7{\theta _{\rm{c}}}}}$ ，${{\textit{z}}_3} = {\omega _{\rm{c}}}{{\rm{e}}^{{\rm{i}}{\theta _{\rm{c}}}}}$ ，${{\textit{z}}_4} = {\omega _{\rm{c}}}{{\rm{e}}^{{\rm{i}}3{\theta _{\rm{c}}}}}$ ；这里取下半平面的单极点留数$ {{\rm{Res}}} ［ f （ {{\textit{z}}_1} ） ］ $ 和${{\rm{Res}}} ［ f （ {{\textit{z}}_2} ） ］ $ 。如果
$\tau {\text{≥}} 0 $ ，则$$ {I_{{\rm{g}}rjs}} ＝ {{\rm{Re}}}\Bigg\{ { - 2\pi {\rm{i}}\sum\limits_{p ＝ 3}^4 {\left[ {\frac{{B （ {{\textit{z}}_{{p}}} ） }}{{A （ {{\textit{z}}_{{p}}} ） }}} \right]} } \Bigg\} ；A （ {{\textit{z}}_{{p}}} ） ＝ \prod\limits_{q ＝ 1，q {\text{≠}} p}^6 { （ {{\textit{z}}_{{p}}} － {{\textit{z}}_q} ） } \text{，} B （ {{\textit{z}}_{{p}}} ） ＝ {\textit{z}}_{{p}}^2{{\rm{e}}^{ - {\rm{i}}{{\textit{z}}_{{p}}}\tau }} \text{，} $$ (20) 其中，
${{\textit{z}}_1} ＝ {\omega _j}{{\rm{e}}^{{\rm{i}}{\theta _j}}}$ ，${{\textit{z}}_2} ＝ - {\omega _j}{{\rm{e}}^{ - {\rm{i}}{\theta _j}}}$ ，${{\textit{z}}_3} ＝ {\omega _{\rm{c}}}{{\rm{e}}^{{\rm{i}}5{\theta _{\rm{c}}}}}$ ，${{\textit{z}}_4} ＝ {\omega _{\rm{c}}}{{\rm{e}}^{{\rm{i}}7{\theta _{\rm{c}}}}}$ ；${{\textit{z}}_5} ＝ {\omega _{\rm{c}}}{{\rm{e}}^{{\rm{i}}{\theta _{\rm{c}}}}}$ ，${{\textit{z}}_6} ＝ {\omega _{\rm{c}}}{{\rm{e}}^{{\rm{i}}3{\theta _{\rm{c}}}}}$ ；这里取下半平面的单极点留数$ {{\rm{Res}}} ［ f （ {{\textit{z}}_3} ） ］ $ 和${{\rm{Res}}} ［ f （ {{\textit{z}}_4} ） ］ $ 。如果
$ \tau ＜ 0 $ ，则$$ {I_{{\rm{g}}rjs}} ＝ {{\rm{Re}}} \Bigg\{ { - 2\pi {\rm{i}}\sum\limits_{p ＝ 1}^4 {\left[ {\frac{{B （ {{\textit{z}}_{{p}}} ） }}{{A （ {{\textit{z}}_{{p}}} ） }}} \right]} } \Bigg\} ；A （ {{\textit{z}}_{{p}}} ） ＝ \prod\limits_{q ＝ 1，q {\text{≠}} p}^6 { （ {{\textit{z}}_{{p}}} － {{\textit{z}}_q} ） } \text{，} B （ {{\textit{z}}_{{p}}} ） ＝ {\textit{z}}_{{p}}^2{{\rm{e}}^{{\rm{i}}{{\textit{z}}_{{p}}}\tau }} \text{，} $$ (21) 其中，
${{\textit{z}}_1} ＝ {\omega _j}{{\rm{e}}^{ - {\rm{i}}{\theta _j}}}$ ，${{\textit{z}}_2} ＝ - {\omega _j}{{\rm{e}}^{{\rm{i}}{\theta _j}}}$ ，${{\textit{z}}_3} ＝ {\omega _{\rm{c}}}{{\rm{e}}^{{\rm{i}}5{\theta _{\rm{c}}}}}$ ，${{\textit{z}}_4} ＝ {\omega _{\rm{c}}}{{\rm{e}}^{{\rm{i}}7{\theta _{\rm{c}}}}}$ ；${{\textit{z}}_5} ＝ {\omega _{\rm{c}}}{{\rm{e}}^{{\rm{i}}{\theta _{\rm{c}}}}}$ ，${{\textit{z}}_6} ＝ {\omega _{\rm{c}}}{{\rm{e}}^{{\rm{i}}3{\theta _{\rm{c}}}}}$ ；这里取下半平面的单极点留数$ {{\rm{Res}}} ［ f （ {{\textit{z}}_1} ） ］ $ ，$ {{\rm{Res}}} ［ f （ {{\textit{z}}_2} ） ］ $ ，$ {{\rm{Res}}} ［ f （ {{\textit{z}}_3} ） ］ $ 和$ {{\rm{Res}}} ［ f （ {{\textit{z}}_4} ） ］ $ 。如果
$\tau {\text{≥}} 0 $ ，则$$ {I_{irjs}} ＝ {{\rm{Re}}}\Bigg\{ { - 2\pi {\rm{i}}\sum\limits_{p ＝ 3， 4， 7， 8} {\left[ {\frac{{B （ {{\textit{z}}_{{p}}} ） }}{{A （ {{\textit{z}}_{{p}}} ） }}} \right]} } \Bigg\} ； A （ {{\textit{z}}_{{p}}} ） ＝ \prod\limits_{q ＝ 1，q {\text{≠}} p}^8 { （ {{\textit{z}}_{{p}}} － {{\textit{z}}_q} ） } \text{，} B （ {{\textit{z}}_{{p}}} ） ＝ {\textit{z}}_{{p}}^4{{\rm{e}}^{ - {\rm{i}}{{\textit{z}}_{{p}}}\tau }} \text{，} $$ (22) 其中
$ {{\textit{z}}_1} ＝ {\omega _j}{{\rm{e}}^{{\rm{i}}{\theta _j}}} $ ，$ {{\textit{z}}_2} ＝ - {\omega _j}{{\rm{e}}^{ - {\rm{i}}{\theta _j}}} $ ，${{\textit{z}}_3} ＝ {\omega _{\rm{c}}}{{\rm{e}}^{{\rm{i}}5{\theta _{\rm{c}}}}}$ ，${{\textit{z}}_4} ＝ {\omega _{\rm{c}}}{{\rm{e}}^{{\rm{i}}7{\theta _{\rm{c}}}}}$ ；${{\textit{z}}_5} ＝ {\omega _{\rm{c}}}{{\rm{e}}^{{\rm{i}}{\theta _{\rm{c}}}}}$ ，${{\textit{z}}_6} ＝ {\omega _{\rm{c}}}{{\rm{e}}^{{\rm{i}}3{\theta _{\rm{c}}}}}$ ，${{\textit{z}}_7} ＝ {\omega _i}{{\rm{e}}^{ - {\rm{i}}{\theta _i}}}$ ，${{\textit{z}}_8} ＝ - {\omega _i}{{\rm{e}}^{{\rm{i}}{\theta _i}}}$ ；这里取下半平面的单极点留数$ {{\rm{Res}}} ［ f （ {{\textit{z}}_3} ） ］ $ ，$ {{\rm{Res}}} ［ f （ {{\textit{z}}_4} ） ］ $ ，$\ {{\rm{Res}}} ［ f （ {{\textit{z}}_7} ） ］ $ 和$ {{\rm{Res}}} ［ f （ {{\textit{z}}_8} ） ］ $ 。如果
$ \tau ＜ 0 $ ，则$$ {I_{irjs}} ＝ {{\rm{Re}}} \Bigg\{{ - 2\pi {\rm{i}}\sum\limits_{p ＝ 1}^4 {\left[ {\frac{{B （ {{\textit{z}}_{{p}}} ） }}{{A （ {{\textit{z}}_{{p}}} ） }}} \right]} } \Bigg\} ； A （ {{\textit{z}}_{{p}}} ） ＝ \prod\limits_{q ＝ 1，q {\text{≠}} p}^8 { （ {{\textit{z}}_{{p}}} － {{\textit{z}}_q} ） } \text{，} B （ {{\textit{z}}_{{p}}} ） ＝ {\textit{z}}_{{p}}^4{{\rm{e}}^{{\rm{i}}{{\textit{z}}_{{p}}}\tau }} \text{，} $$ (23) 其中，
${{\textit{z}}_1} ＝ {\omega _j}{{\rm{e}}^{ - {\rm{i}}{\theta _j}}}$ ，${{\textit{z}}_2} ＝ - {\omega _j}{{\rm{e}}^{{\rm{i}}{\theta _j}}}$ ，${{\textit{z}}_3} ＝ {\omega _{\rm{c}}}{{\rm{e}}^{{\rm{i}}5{\theta _{\rm{c}}}}}$ ，${{\textit{z}}_4} ＝ {\omega _{\rm{c}}}{{\rm{e}}^{{\rm{i}}7{\theta _{\rm{c}}}}}$ ；${{\textit{z}}_5} ＝ {\omega _{\rm{c}}}{{\rm{e}}^{{\rm{i}}{\theta _{\rm{c}}}}}$ ，${{\textit{z}}_6} ＝ {\omega _{\rm{c}}}{{\rm{e}}^{{\rm{i}}3{\theta _{\rm{c}}}}}$ ，${{\textit{z}}_7} ＝ {\omega _i}{{\rm{e}}^{{\rm{i}}{\theta _i}}}$ ，${{\textit{z}}_8} ＝ - {\omega _i}{{\rm{e}}^{ - {\rm{i}}{\theta _i}}}$ ；这里取下半平面的单极点留数$ {{\rm{Res}}} ［ f （ {{\textit{z}}_1} ） ］ $ ，$ {{\rm{Res}}} ［ f （ {{\textit{z}}_2} ） ］ $ ，$ {{\rm{Res}}} ［ f （ {{\textit{z}}_3} ） ］ $ 和$ {{\rm{Res}}} ［ f （ {{\textit{z}}_4} ） ］ $ 。以上诸式中，
${\theta _i} ＝ \arcsin {\zeta _i}$ ，${\theta _j} ＝ \arcsin {\zeta _j}$ ，${\theta _{\rm{c}}} ＝ \dfrac{\pi }{4}$ 。3.3 Penzien-Clough自功率谱密度下相关系数的解析表达式
式（16）的积分形式为
$$ {I_{{\rm{g}}r{\rm{g}}s}} ＝ {{\rm{Re}}}\Bigg\{ { - 2\pi {\rm{i}}\sum\limits_{p ＝ 1}^2 {\left[ {\frac{{B （ {{\textit{z}}_{{p}}} ） }}{{A （ {{\textit{z}}_{{p}}} ） }}} \right]} } \Bigg\} ；A （ {{\textit{z}}_{{p}}} ） ＝ \prod\limits_{q ＝ 1，q {\text{≠}} p}^4 { （ {{\textit{z}}_{{p}}} - {{\textit{z}}_q} ） } \text{，} B （ {{\textit{z}}_{{p}}} ） ＝ {{\rm{e}}^{ - {\rm{i}}{{\textit{z}}_{{p}}}\left| \tau \right|}} \text{，} $$ (24) 其中，
${{\textit{z}}_1} ＝ {\omega _{\rm{f}}}{{\rm{e}}^{ - {\rm{i}}{\theta _{\rm{f}}}}}$ ，${\textit{z}}{}_2 ＝ - {\omega _{\rm{f}}}{{\rm{e}}^{{\rm{i}}{\theta _{\rm{f}}}}}$ ，${{\textit{z}}_3} ＝ {\omega _{\rm{f}}}{{\rm{e}}^{{\rm{i}}{\theta _{\rm{f}}}}}$ ，${{\textit{z}}_4} ＝ - {\omega _{\rm{f}}}{{\rm{e}}^{ - {\rm{i}}{\theta _{\rm{f}}}}}$ ；这里取下半平面的单极点留数$ {{\rm{Res}}} ［ f （ {{\textit{z}}_1} ） ］ $ 和$ {{\rm{Res}}} ［ f （ {{\textit{z}}_2} ） ］ $ 。如果
$ \tau {\text{≥}} 0$ ，则$$ {I_{{\rm{g}}rjs}} ＝ {{\rm{Re}}} \Bigg\{{ - 2\pi {\rm{i}}\sum\limits_{p ＝ 3}^4 {\left[ {\frac{{B （ {{\textit{z}}_{{p}}} ） }}{{A （ {{\textit{z}}_{{p}}} ） }}} \right]} } \Bigg\} ；A （ {{\textit{z}}_{{p}}} ） ＝ \prod\limits_{q ＝ 1，q {\text{≠}} p}^6 { （ {{\textit{z}}_{{p}}} － {{\textit{z}}_q} ） } \text{，} B （ {{\textit{z}}_{{p}}} ） ＝ {\textit{z}}_{{p}}^2{{\rm{e}}^{ - {\rm{i}}{{\textit{z}}_{{p}}}\tau }} \text{，} $$ (25) 其中，
${{\textit{z}}_1} ＝ {\omega _j}{{\rm{e}}^{{\rm{i}}{\theta _j}}}$ ，$ {{\textit{z}}_2} ＝ - {\omega _j}{{\rm{e}}^{ - {\rm{i}}{\theta _j}}} $ ，${{\textit{z}}_3} ＝ {\omega _{\rm{f}}}{{\rm{e}}^{ - {\rm{i}}{\theta _{\rm{f}}}}}$ ，${\textit{z}}_4 ＝ - {\omega _{\rm{f}}}{{\rm{e}}^{{\rm{i}}{\theta _{\rm{f}}}}}$ ；${{\textit{z}}_5} ＝ {\omega _{\rm{f}}}{{\rm{e}}^{{\rm{i}}{\theta _{\rm{f}}}}}$ ，${{\textit{z}}_6} ＝ - {\omega _{\rm{f}}}{{\rm{e}}^{ - {\rm{i}}{\theta _{\rm{f}}}}}$ ；这里取下半平面的单极点留数$ {{\rm{Res}}} ［ f （ {{\textit{z}}_3} ） ］ $ 和$ {{\rm{Res}}} ［ f （ {{\textit{z}}_4} ） ］ $ 。如果
$\tau ＜ 0 $ ，则$$ {I_{{\rm{g}}rjs}} ＝ {{\rm{Re}}} \Bigg\{ { - 2\pi {\rm{i}}\sum\limits_{p ＝ 1}^4 {\left[ {\frac{{B （ {{\textit{z}}_{{p}}} ） }}{{A （ {{\textit{z}}_{{p}}} ） }}} \right]} } \Bigg\} ；A （ {{\textit{z}}_{{p}}} ） ＝ \prod\limits_{q ＝ 1，q {\text{≠}} p}^6 { （ {{\textit{z}}_{{p}}} － {{\textit{z}}_q} ） } \text{，} B （ {{\textit{z}}_{{p}}} ） ＝ {\textit{z}}_{{p}}^2{{\rm{e}}^{{\rm{i}}{{\textit{z}}_{{p}}}\tau }} \text{，} $$ (26) 其中，
${{\textit{z}}_1} ＝ {\omega _j}{{\rm{e}}^{ - {\rm{i}}{\theta _j}}}$ ，${{\textit{z}}_2} ＝ - {\omega _j}{{\rm{e}}^{{\rm{i}}{\theta _j}}}$ ，${{\textit{z}}_3} ＝ {\omega _{\rm{f}}}{{\rm{e}}^{ - {\rm{i}}{\theta _{\rm{f}}}}}$ ，${\textit{z}}{}_4 ＝ - {\omega _{\rm{f}}}{{\rm{e}}^{{\rm{i}}{\theta _{\rm{f}}}}}$ ；${{\textit{z}}_5} ＝ {\omega _{\rm{f}}}{{\rm{e}}^{{\rm{i}}{\theta _{\rm{f}}}}}$ ，${{\textit{z}}_6} ＝ - {\omega _{\rm{f}}}{{\rm{e}}^{ - {\rm{i}}{\theta _{\rm{f}}}}}$ ；这里取下半平面的单极点留数$ {{\rm{Res}}} ［ f （ {{\textit{z}}_1} ） ］ $ ，$ {{\rm{Res}}} ［ f （ {{\textit{z}}_2} ） ］ $ ，$ {{\rm{Res}}} ［ f （ {{\textit{z}}_3} ） ］ $ 和$ {{\rm{Res}}} ［ f （ {{\textit{z}}_4} ） ］ $ 。如果
$ \tau {\text{≥}} 0$ ，则$$ {I_{irjs}} ＝ {{\rm{Re}}} \Bigg\{ { - 2\pi {\rm{i}}\sum\limits_{p ＝ 3{{， 4， 7， 8}}} {\left[ {\frac{{B （ {{\textit{z}}_{{p}}} ） }}{{A （ {{\textit{z}}_{{p}}} ） }}} \right]} } \Bigg\} ；A （ {{\textit{z}}_{{p}}} ） ＝ \prod\limits_{q ＝ 1，q {\text{≠}} p}^8 { （ {{\textit{z}}_{{p}}} － {{\textit{z}}_q} ） } \text{，} B （ {{\textit{z}}_{{p}}} ） ＝ {\textit{z}}_{{p}}^4{{\rm{e}}^{ - {\rm{i}}{{\textit{z}}_{{p}}}\tau }} \text{，} $$ (27) 其中，
${{\textit{z}}_1} ＝ {\omega _j}{{\rm{e}}^{{\rm{i}}{\theta _j}}}$ ，$ {{\textit{z}}_2} ＝ - {\omega _j}{{\rm{e}}^{ - {\rm{i}}{\theta _j}}} $ ，$ {{\textit{z}}_3} ＝ {\omega _{\rm{f}}}{{\rm{e}}^{ - {\rm{i}}{\theta _{\rm{f}}}}} $ ，$ {{\textit{z}}_4} ＝ - {\omega _{\rm{f}}}{{\rm{e}}^{{\rm{i}}{\theta _{\rm{f}}}}} $ ；${{\textit{z}}_5} ＝ {\omega _{\rm{f}}}{{\rm{e}}^{{\rm{i}}{\theta _{\rm{f}}}}}$ ，${{\textit{z}}_6} ＝ - {\omega _{\rm{f}}}{{\rm{e}}^{ - {\rm{i}}{\theta _{\rm{f}}}}}$ ，${{\textit{z}}_7} ＝ {\omega _i}{{\rm{e}}^{ - {\rm{i}}{\theta _i}}}$ ，$ {{\textit{z}}_8} ＝ - {\omega _i}{{\rm{e}}^{{\rm{i}}{\theta _i}}} $ ；这里取下半平面的单极点留数$ {{\rm{Res}}} ［ f （ {{\textit{z}}_3} ） ］ $ ，$ {{\rm{Res}}} ［ f （ {{\textit{z}}_4} ） ］ $ ，$ {{\rm{Res}}} ［ f （ {{\textit{z}}_7} ） ］ $ 和$ {{\rm{Res}}} ［ f （ {{\textit{z}}_8} ） ］ $ 。如果
$\tau ＜ 0 $ ，则$$ {I_{irjs}} ＝ {{\rm{Re}}}\Bigg\{ { - 2\pi {\rm{i}}\sum\limits_{p ＝ 1}^4 {\left[ {\frac{{B （ {{\textit{z}}_{{p}}} ） }}{{A （ {{\textit{z}}_{{p}}} ） }}} \right]} } \Bigg\} ；A （ {{\textit{z}}_{{p}}} ） ＝ \prod\limits_{q ＝ 1，q {\text{≠}} p}^8 { （ {{\textit{z}}_{{p}}} － {{\textit{z}}_q} ） } \text{，} B （ {{\textit{z}}_{{p}}} ） ＝ {\textit{z}}_{{p}}^4{{\rm{e}}^{{\rm{i}}{{\textit{z}}_{{p}}}\tau }} \text{，} $$ (28) 其中，
${{\textit{z}}_1} ＝ {\omega _j}{{\rm{e}}^{ - {\rm{i}}{\theta _j}}}$ ，${{\textit{z}}_2} ＝ - {\omega _j}{{\rm{e}}^{{\rm{i}}{\theta _j}}}$ ，${{\textit{z}}_1} ＝ {\omega _{\rm{f}}}{{\rm{e}}^{ - {\rm{i}}{\theta _{\rm{f}}}}}$ ，${\textit{z}}_2 ＝ - {\omega _{\rm{f}}}{{\rm{e}}^{{\rm{i}}{\theta _{\rm{f}}}}}$ ；${{\textit{z}}_3} ＝ {\omega _{\rm{f}}}{{\rm{e}}^{{\rm{i}}{\theta _{\rm{f}}}}}$ ，${{\textit{z}}_6} ＝ - {\omega _{\rm{f}}}{{\rm{e}}^{ - {\rm{i}}{\theta _{\rm{f}}}}}$ ，${{\textit{z}}_7} ＝ {\omega _i}{{\rm{e}}^{{\rm{i}}{\theta _i}}}$ ，${{\textit{z}}_8} ＝ - {\omega _i}{{\rm{e}}^{ - {\rm{i}}{\theta _i}}}$ ；这里取下半平面的单极点留数$ {{\rm{Res}}} ［ f （ {{\textit{z}}_1} ） ］ $ ，$ {{\rm{Res}}} ［ f （ {{\textit{z}}_2} ） ］ $ ，$ {{\rm{Res}}} ［ f （ {{\textit{z}}_3} ） ］ $ 和$ {{\rm{Res}}} ［ f （ {{\textit{z}}_4} ） ］ $ 。以上诸式中，
${\theta _i} ＝ \arcsin {\zeta _i}$ ，${\theta _j} ＝ \arcsin {\zeta _j}$ ，${\theta _{\rm{f}}} ＝ \arcsin {\zeta _{\rm{f}}}$ 。4. 精度和计算效率的讨论
4.1 相关系数解析表达式的精度
本节以一座斜拉桥为例进行计算，其中跨240 m，边跨117.5 m，其有限元模型示意图如图10所示。假设地震波沿纵向传播，视波速取1 000 m/s。用式（9）或式（11）的自功率谱密度函数以及式（6）的地震动空间相干模型，计算该桥梁在多点地震激励下的反应。
计算中取前50阶振型。表1—4分别列出了桥梁不同构件中某些位置处的绝对位移、相对位移、轴力、剪力和弯矩。表中：Ln（n＝1—6）代表桥梁左侧的反应，Rn代表右侧的反应；n＝1时代表塔顶的绝对位移；n＝2时代表塔顶和塔底之间的相对位移；n＝3时代表边墩顶和梁端之间的相对位移，这项参数在落梁问题中很重要；n＝4时代表塔底的轴力；n＝5时代表塔底的剪力；n＝6时代表塔底的弯矩。
表 1 Qu等（1996）相干模型下的地震位移响应及相对误差Table 1. Seismic displacement response and its relative error under the coherence model of Qu et al （1996）胡聿贤的自功率谱密度模型 克拉夫−彭津的自功率谱密度模型 L1/cm R1/cm L2/cm R2/cm L3/cm R3/cm L1/cm R1/cm L2/cm R2/cm L3/cm R3/cm A 7.44 7.40 6.27 5.79 7.02 6.23 7.42 7.38 6.22 5.68 7.03 6.08 B 7.44 7.40 6.27 5.79 7.02 6.23 7.42 7.38 6.22 5.68 7.03 6.08 C 7.53 7.49 6.24 5.77 6.97 6.17 7.48 7.44 6.21 5.67 7.02 6.04 eBA 0 0 0 0 0 0 0 0 0 0 0 0 eCA 1.22% 1.22% −0.48% −0.35% −0.71% −0.96% 0.81% 0.81% −0.16% −0.18% −0.14% −0.66% 表 4 Harichandran和Vanmarcke （1986）相干模型下的地震力或弯矩响应及相对误差Table 4. Seismic force or moment response and its relative error under the coherence model of Harichandran and Vanmarcke （1986）胡聿贤的自功率谱密度模型 克拉夫−彭津的自功率谱密度模型 L4/kN R4/kN L5/kN R5/kN L6/kN·m R6/kN·m L4/kN R4/kN L5/kN R5/kN L6/kN·m R6/kN·m A 719 691 1111 1069 29656 27697 712 685 1109 1066 29377 27307 B 719 691 1107 1066 29617 27655 712 685 1106 1064 29355 27284 C 720 693 1111 1073 29688 27794 712 686 1106 1068 29393 27393 eBA 0 0 −0.36% −0.28% −0.13% −0.15% 0 0 −0.27% −0.19% −0.07% −0.08% eCA 0.14% 0.29% 0 0.37% 0.11% 0.35% 0 0.15% −0.27% 0.19% 0.05% 0.31% 定义两个相对误差：
$$ {e_{BA}} ＝ \frac{{B － A}}{A} ， {e_{CA}} ＝ \frac{{C － A}}{A} \text{，} $$ (29) 式中：A表示基于精确APSD模型和精确空间相干模型的响应的数值积分解；B表示基于精确APSD模型和近似空间相干模型的响应的数值积分解；C为3.2和3.3节中提出的解析解；
${e_{BA}} $ 和$ {e_{CA}} $ 是相对误差。在方法A和B中，使用了变步长的辛普森积分法。在积分极限从1 000，500，200到100 rad/s的情况下，将基于胡聿贤的功率谱密度函数与屈铁军、王君杰的空间相干函数的系数
$\;{\rho _{{\rm{g}}r{\rm{g}}s}}$ ，$ \;{\rho _{{\rm{g}}rjs}} $ ，$\; {\rho _{irjs}} $ 进行了比较，最终方法A采用200 rad/s。当相对积分精度分别为$ {10^{ - 8}} $ ，$ {10^{ - 5}} $ ，$ {10^{ - 4}} $ 和$ {10^{ - 3}} $ 时，分别计算出系数$\;{\rho _{{\rm{g}}r{\rm{g}}s}}$ ，$\; {\rho _{{\rm{g}}rjs}} $ ，$\; {\rho _{irjs}} $ ，并采用了精度为$ {10^{ - 4}} $ 时的积分步长。从表1—4可以看出，相对误差非常小，这说明近似空间相干函数对结构的反应影响非常小。从3.1节中的图7—9也可以得出这个结论。
表 2 Harichandran和Vanmarcke （1986）相干模型下的地震位移响应及相对误差Table 2. Seismic displacement response and its relative error under the coherence model of Harichandran and Vanmarcke （1986）胡聿贤的自功率谱密度模型 克拉夫−彭津的自功率谱密度模型 L1/cm R1/cm L2/cm R2/cm L3/cm R3/cm L1/cm R1/cm L2/cm R2/cm L3/cm R3/cm A 7.27 7.24 6.24 5.81 7.08 6.40 7.26 7.22 6.19 5.69 7.10 6.27 B 7.27 7.24 6.24 5.81 7.09 6.41 7.25 7.22 6.19 5.69 7.10 6.27 C 7.36 7.32 6.21 5.78 7.05 6.36 7.31 7.28 6.18 5.68 7.09 6.25 eBA 0 0 0 0 0.14% 0.16% −0.14% 0 0 0 0 0 eCA 1.24% 1.10% −0.48% −0.52% −0.42% −0.63% 0.69% 0.83% −0.16% −0.18% −0.14% −0.32% 表 3 Qu等（1996）相干模型下的地震力或弯矩响应及相对误差Table 3. Seismic force or moment response and its relative error under the coherence model of Qu et al （1996）胡聿贤的自功率谱密度模型 克拉夫−彭津的自功率谱密度模型 L4/kN R4/kN L5/kN R5/kN L6/kN·m R6/kN·m L4/kN R4/kN L5/kN R5/kN L6/kN·m R6/kN·m A 730 699 1058 1011 29192 27002 723 692 1057 1010 28931 26618 B 730 699 1054 1008 29159 26966 723 692 1055 1008 28915 26601 C 731 702 1057 1016 29219 27107 723 694 1054 1012 28944 26714 eBA 0 0 −0.38% −0.30% −0.11% −0.13% 0 0 −0.19% −0.20% −0.06% −0.06% eCA 0.14% 0.43% −0.09% 0.49% 0.09% 0.39% 0 0.29% −0.28% 0.20% 0.04% 0.36% 除了塔顶的绝对位移（即L1和R1）外，A和C方法的反应相对误差都小于1%。同时在地震分析中，绝对位移并不重要，而相对位移却很重要。这表明解析表达式可以与数值方法很好地吻合。
4.2 计算效率的比较
表5列出了数值方法A和解析方法C对该例桥梁的计算时间，计算中取前50阶振型。表6列出了计算中不同阶振型的计算时间。在本例中，只计算了以下反应：① 四个塔顶的绝对位移；② 四个塔顶与塔底之间的相对位移；③ 两个边墩顶部与梁端之间的相对位移；④ 四个塔底的轴力、剪力和弯矩。在本例中，使用英特尔Visual Fortran 11.1.067编写有限元程序，前50阶的特征值计算耗时50 s。
表 5 计算时间比较（单位：s）Table 5. Comparison of time consumption of computation （Unit：s）从表5可以发现，解析法求解的时间消耗只有数值求解的1/52左右。就绝对时间而言，解析法需要93—111 s，这在实际设计中是完全可以接受的。同时，数值解决方案需要5 833—6 273 s，比解析法要多耗时近98%。从表6中也可以得出同样的结论，计算时间随着所取振型阶数的增加而增加。对于大型桥梁，在计算中可能需要取数百阶振型。与用数值法相比，采用解析法计算大型结构响应可大大缩短计算用时。
表 6 计算时间比较（单位：s）Table 6. Comparison of time consumption of computation （Unit：s）5. 结论
本文在提出近似空间相干函数的基础上，建立了相关系数
$\;{\rho _{{\rm{g}}r{\rm{g}}s}}$ ，$ \;{\rho _{{\rm{g}}rjs}} $ 和$ \;{\rho _{irjs}} $ 的解析表达式，一方面保证了计算结果的精度，另一方面采用解析法可以更快速地获得三个相关系数的积分值。同时使用了仅需要一至两个地震动参数的APSD近似表达式。通过对数值计算结果的分析，本文所提相关系数$\;{\rho _{{\rm{g}}r{\rm{g}}s}}$ ，$ \;{\rho _{{\rm{g}}rjs}} $ 和$\; {\rho _{irjs}} $ 的解析式有足够的工程精度和较高的计算效率，可为平面大尺度结构多点激励下的地震响应计算提供高效方法。 -
-
期刊类型引用(2)
1. 吴堃，马杰，苏俊省，陈宇，黄信，齐麟，李长辉. 考虑局部场地效应的深水多跨连续梁桥纵桥向地震响应分析. 振动工程学报. 2025(01): 154-161 . 
百度学术
2. 李宇轩. 高层建筑结构抗震性能分析. 安徽建筑. 2025(04): 79-81 . 百度学术
其他类型引用(2)
下载:
计量
- 文章访问数: 1375
- HTML全文浏览量: 18
- PDF下载量: 80
- 被引次数: 4
