引言

大量观测事实表明，地震事件之间存在一定的相关性，特别是同一构造区域内相继发生的主震与余震、以及空间上相邻的断层活动存在着显著的时空与力学关联。研究地震间的相关性，是理解地震序列演化规律、评估区域地震危险性、以及预测地震动强度参数的核心科学问题之一。

为量化这种相关性，赵银刚基于全球地震数据，通过线性回归方法分析了主震与最大余震之间的统计关系。研究发现，主震震级与最大余震震级之间存在显著的线性关系，且这一关系在不同断层类型中均成立，而这一关系可用于快速估计地震发生后的最大余震震级。王元杰基于国内外11次中、强地震的近场与远场强震动记录，采用方差主轴方法，对地震动三个平动分量之间的相关性进行了分析。结果表明：在强震段与下降段，最小方差主轴方向基本保持竖直；而在上升段，主轴方向变化较为复杂。

随着强震观测数据的积累和统计建模方法的进步，研究者开始运用如Copula函数等现代统计工具，定量刻画地震动参数之间（如峰值加速度、谱加速度、持时等）复杂的非线性、非对称关系。这类方法不仅能够描述单一地震中不同分量间强度参数的相关性，更可扩展到主震与余震参数间的联合概率建模，从而为概率地震危险性分析的发展奠定基础。

本文旨在对Copula函数在地震动强度参数相关性研究中的应用进行系统性综述。首先阐述其模型建立的大致流程；继而，分类评述其在主余震分析、竖向和水平地震相关性的典型案例；随后，深入讨论当前面临的主要挑战与局限性；最后，对未来可能的研究方向进行展望，以期为相关领域的研究者与工程师提供一个清晰的学术脉络与发展路径参考。

1 Copula函数理论基础与模型选择

1.1 核心原理：Sklar定理

Copula理论的核心是Sklar于1959年提出的著名定理。根据Sklar定理可知，若为一个具有边缘分布的维联合分布函数。那么，存在一个Copula函数，使得对于所有，有：

(1)

Copula函数的形式多样，能够描述变量之间各种的复杂相关结构，其所有一维边缘分布均服从[0,1]区间上的均匀分布。目前对地震间的相关性研究通常二维分布的Sklar定理。

假设)是和的联合累积分布函数，那么根据二维分布的Sklar定理可知，存在一个Copula函数将X1和X2的边缘分布函数连接起来得到X1和X2的联合分布函数，如下式所示：

(2)

(3)

(4)

式中：，,、为变量和的边缘累积分布函数，均为[0,1]间的均匀分布；和分别为和的边缘概率密度函数；为Copula分布函数，为Copula函数的密度函数；和分别为和的联合累积分布函数和联合概率密度函数。

2.2 二维Copula函数建立基本流程

构建二维Copula函数的联合分布模型通常包括三个步骤：(1)为拟合主震、余震各向地震动强度参数的边缘分布，本文将正态分布、Gumbel分布及广义极值分布三种分布函数作为备选模型。基于三者的K-S检验，经初步验证后，依据AIC准则和BIC准则信息准则，筛选出最优的边缘分布函数。（2）计算二者间的相关系数，一般计算pearson线性相关系数，Kendall秩相关系数以及Spearman秩相关系数。（3）确定最优的Copula函数族并估计Copula参数。参数估计常用极大似然法或矩估计法（一般通过Kendall‘s τ推断参数）。模型选择以AIC准则和BIC准则等信息准则为基础

2.3 地震相关性中常用Copula函数族

为了确定最优的Copula函数，通常选取以下五种常用Copula函数作为备选Copula函数，它们分别是Gaussian Copula函数、t-Copula函数、Frank Copula函数、Clayton Copula函数和Gumbel Copula函数。

Gaussian Copula常被用来作为基准，评估其它复杂的 Copula 模型的性能，因此比较Gaussian Copula与t-Copula或Gumbel Copula的拟合优度就可以客观地量化数据中尾部相关性的显著程度，进而据此判断是否需要采用更复杂的模型。具体而言，t-Copula有明确而突出的优势：它最擅长刻画对称尾部相关性，因而很适合于建模强震-强余震这类极端事件的联合风险。因此，在实际应用中，研究者宜先借助数据分析方法（例如绘制经验尾部相关图）判断数据中是否存在显著的尾部相关性，若尾部相关特征明显且大致对称，t-Copula就是一个十分自然的候选模型。与此形成极好补充的是，Gumbel Copula是分析主余震强度参数相关性，尤其是强震-强余震极端风险传递机制时极其有力的工具，其上尾相关的特性本身与地震序列的物理认知高度契合，故在应用上也更可靠、更合理。Clayton Copula是阿基米德族Copula中又一重要成员，其最本质、最突出的数学特性是刻画非对称的下尾相关性，因此它在主余震相关性分析中有明确、独特的应用场景和解释价值，但也毋庸讳言，相较于Gumbel Copula,其适用范围更有限。FrankCopula 最核心的特点是能够描绘出变量之间对称而灵活的相互依存关系，而没有尾部的关联性。因此，当主震与余震的某一对强度参数存在较强的对称相关性时，Frank Copula相较于Gaussian Copula是更灵活、更恰当的选择。

2. 地震动水平和竖向强度参数的相关性研究

近年来，国内外研究学者对水平向地震动强度参数间的相关性进行了大量研究：BRADLEY^[12]以NGA数据库为基础，开发了水平向地震动峰值地面速度（PGV）与多种基于频谱的地面运动强度测量之间的经验相关方程，通过活动浅层地壳地震的地表运动和四个地表运动预测方程，得到了峰值地面加速度、5%阻尼伪光谱加速度、加速度谱强度和光谱强度间的相关性模型；BAKER等以太平洋地震工程研究中心数据库为基础，分析水平向地震动若干强度参数间的相关性，发现最大方向谱的相关性与其他谱加速度定义的相关性基本一致，且这些相关性对震级、距离或场地条件等因素无显著依赖。

目前的研究显示，从传统的横向到参数之间的研究，地动强度参数之间的相关性分析已经逐步扩展到水平和竖向强度参数共同相关的领域，从而提供了一个全新的视角来对地动三维效应进行全面评估。为深入认识近断层地震动的三维特征，李宁基于PEER强震数据库中的890组记录，对竖向与水平加速度峰值比随矩震级、断层距、场地条件及断层类型的变化规律进行了系统分析，并对近断层脉冲特性开展了对比研究。结果表明：竖向脉冲在出现数量、分布范围及其峰值、周期等衰减特性上均与水平脉冲存在明显差异。周正华以国内外多次地震近场自由地表加速度记录为基础，重点分析了竖向与水平向加速度峰值比及反应谱比值的演变规律，并特别以1999年我国台湾地区集集地震为例，揭示了竖向反应谱的特征。研究发现，V/H比值与震源距呈负相关，并随周期、场地条件及震源距的变化而呈现系统性差异。韩建平以9次地震事件中获取的86组三分量加速度记录为基础，重点分析了突出周期比值的统计规律，即纵向和横向加速度的峰值比值，加速度的反应谱比值和优异周期比值。

在PEER NGA-West2地震动数据库中，三分量地震波的方向定义为“传感器相对方向”，而非绝对的“地理方向”（如正北、正东）。数据本身并不包含固定的地理方位信息。因此，在分析所选取水平分量地震动记录时，先合成处理两个水平分量记录，从而使地震动强度参数的相关性分析结果更为可靠。王晓磊以975组地震动记录为基础，采用四种水平向地震动合成方法——几何平均值（Gm）、较大值（Lar）以及RotDnn，分别对主震和余震的两条正交水平地震动记录进行合成，确定了水平向和竖向地震动强度参数的相关系数和最优边缘分布函数，最终通过AIC和BIC准则确定全数据水平地震动RotD50强度参数和竖向地震动强度参数间的最优Copula函数。通过建立水平地震动和竖向地震动的联合分布函数，实现了水平向地震动强度参数条件下，竖向参数的条件均值预测，并对不同场地类别对竖向和水平向地震动相关性做了进一步探讨。研究结果表明：Copula函数可有效建立水平和竖向地震动强度参数的联合分布函数；基于四种方法合成的水平地震动强度参数与竖向地震动强度参数间的相关性基本相同.

3. 主余震强度参数相关性研究及余震模型预测

强震事件常伴随显著的地震序列活动，主震发生后往往会触发一系列余震，构成一个物理上连贯的灾害过程。传统概率地震危险性分析在建模时，通常将主震与余震视为统计独立事件，或仅采用经验性的简化关系进行描述，这种处理方式难以反映主震在强度、频谱及空间分布等方面对后续余震活动的影响。由于余震是由主震孕育而来的，从震源机制、震源位置、震级等方面看，余震和主震之间存在一定的联系。

目前，国内外研究学者对主余震地震动的相关性进行了大量研究：施雨捷在752组真实主余震序列的基础上，采用多元非线性回归方法，建立了主震与最大余震地震动参数间的数学关系。结果表明，余震相对地震动参数与主震矩震级及场地剪切波速呈正相关，与主余震震级比、剪切波速比呈负相关，与震源距相关性较弱。在车型比较方面，在VS30较低的情况下，预测结果接近CY2008车型；VS30较高时，与AS2008车型的误差不超过5%。这种模式可以更好地对余震地震动参数进行预测，提供理论基础，用于地震危险性分析和工程应用。丁佳伟基于NGA-West2数据库中的8次主余震事件（共560条记录），分析了主震与余震空间相关特性的差异。通过ASK14预测方程计算事件内残差，并采用地统计学半变异函数模型量化了谱加速度的空间自相关性。同时，结合Pearson相关系数与马尔可夫假设，建立了考虑空间距离的互相关模型。结果表明，余震在长周期段表现出比主震更强的空间相关性。温卫平^[21]通过建立与主震参数相关的余震地震动条件预测公式，揭示了主余震强度参数间的统计依赖关系。

已有研究证实，主震与余震的地震动强度参数间具有明显的相关性。朱瑞广以662条主余震地震动为基础，通过确定主余震强度参数间的相关系数和最优边缘分布函数，确定了主余震强度参数间的最优Copula函数。并通过基于Copula函数的主余震联合分布函数模型，给出了给定主震地震动参数条件下,余震地震动强度参数的条件分布和条件均值。研究结果表明:可以通过Copula条件均值对余震地震动强度参数的取值进行预测。该方法的核心逻辑在于，利用Copula函数构建主震与余震关键强度参数的联合概率分布。基于Sklar定理，该联合分布可分解为描述主震参数边缘分布、余震参数边缘分布以及二者间复杂相依结构的Copula函数。通过此模型，可精确推导出在给定主震强度参数观测值的条件下，余震相应参数超越某一阈值的条件概率。

上述基于Copula的条件概率建模方法，为深入分析余震自身参数间的相关性结构提供了有力工具。白治朋该研究基于NGA-West2数据库与CB14地震动预测方程，分析了2268条余震记录的谱加速度相关性，建立了余震谱加速度相关系数预测模型。该模型能有效表征不同周期（0.01~10 s）谱加速度间的相关结构，且预测误差较低。在系统把握余震参数间相关性的基础上，研究进一步向竖向地震动这一特定分量深化。王霈宇通过构建余震竖向地震动预测方程，建立主余震谱形参数相关模型，并提出余震竖向条件均值谱概念，为余震竖向地震动选取提供了科学合理的目标谱，对震后建筑安全评定等关键任务具有重要的理论指导意义。在实现参数层面的条件预测基础上，更为系统的概率风险分析框架得以发展。王泽龙提出了一种考虑主余震相关性的概率危险性分析方法：利用主震与余震在震级、断层信息等方面的统计关系，通过蒙特卡洛模拟生成余震场景；进而以余震条件均值谱为目标，计算给定主震条件下的余震危险性曲线及一致危险性谱。最终实现对区域余震地震动影响场的概率性划定。研究表明：该方法改进了传统分析中忽视序列地震相关性的不足，通过考虑主余震在多个方面的关联性，可以更合理地评估强震后余震对周围场地的威胁程度。

4. 当前挑战与未来展望

尽管Copula函数的应用取得了显著进展，但在理论深度和工程实用化方面仍面临一系列挑战，这也指明了未来的研究方向。

4.1 核心挑战

高维建模的复杂程度和计算成本：实际的工程问题往往会涉及到几十个相关变量，甚至是几百个相关变量。而Vine Copula（藤Copula）把高维联合分布自然、合理地分解为若干条件二元Copula对，可以进行更好的模块化建模，但其结构选择、参数估计在高维情况下变得异常复杂，计算量极大。

非平稳性与时变相关性：由于目前绝大多数Copula模型都假定相依结构是静态的，而地震动本身具有非平稳性，主余震序列中破裂过程又随时间演化，故变量间的相关性很自然会随时间或地震动强度变化。因此，发展能合理描述动态相依关系的时变Copula模型是一个重要而困难的课题。

由于目前Copula模型主要是从强度参数数据的角度做统计拟合，故其与地震震源物理、波动传播机理的直接联系十分薄弱，因此，如何将物理约束自然、合理地引入Copula参数之中，构造物理信息增强的Copula模型，是提高模型外推、预测可靠性的一大挑战。

数据稀缺性与模型不确定性：由于近断层脉冲型地震动、特大震级事件及某些构造环境下的观测记录都非常稀少，故用现有数据估计模型参数时必然伴随着较大的不确定性，而预测极端罕见事件时又会放大这种不确定性。因此，量化Copula模型本身的不确定性并合理传播到风险结果中，是当前研究的明显缺口。

4.2未来研究方向

目前基于Copula的文献大多从“点对点”参数对（如主震与余震强度参数）间的相关性来研究，因此未来一个有潜力的突破方向是把所识别的若干二元Copula关系自然、合理的整合为一个系统的参数影响力网络，从而揭示主震破坏潜力向余震传递的“关键路径”与内在机理。具体步骤可归纳如下： 首先把主震的各类强度参数（如PGA、PGV、不同周期的Sa、持时参数等）视为“源节点”，并将余震的对应以及其它强度参数视为“目标节点”，用最优Copula模型的相依强度（如Kendall’s τ）作为连接边的权重，进而构建一个加权的二分网络模型。在此基础上，借助网络科学中的方法（如最大权重匹配、关键路径分析）进行挖掘。

首先，这种方法实现了从描述孤立“二元关系”到理解“系统关联”的认知转变，将主余震系统视为一个整体网络，因而更有利于系统的讨论不同强度参数在风险传递中的相对重要性。其次，通过识别“关键路径”，能够揭示主震破坏如何通过特定的地震动特征维度（如特定频带能量、累积效应）被筛选、放大并传递给余震，从而为超越统计关联、探索其背后的物理机制提供切实可行结构化分析框架，所得结果直接对工程风险分析有极强的指导意义：明确在余震危险性评估中应优先关注的主震特征参数，据此发展更具物理基础、计算更高效的概率预测模型与地震动筛选方法。这一方向也彰显了Copula理论与复杂网络科学进行学科交叉的潜力，为解决地震工程中高维非线性系统的分析难题提供了新思想。这一研究方向将推动相关研究从静态的“相关性描述”迈向动态的“机理解释”，为基于物理机制和数据驱动相结合的主余震风险精细化评估开辟新途径。

5结论

本文系统梳理了Copula函数在地震动强度参数相关性研究中的应用进展。作为一种灵活的统计工具，Copula通过建立地震动强度参数的联合分布函数，可以得到给定一个强度参数的条件下，另一个强度参数的条件分布和条件均值。综述表明，该理论方法主要用于以下两个方面：

在水平与竖向地震动参数联合建模方面，通过构建竖向分量及其与水平分量的联合分布模型，实现了基于水平向参数对竖向参数的条件预测。在主余震序列相关性分析与风险量化方面，通过主余震强度参数的联合分布，可以得到在给定主震条件下的余震强度参数的概率预测，显著提升了余震危险性评估的合理性与精细化水平，为震后应急决策与韧性设计提供了更科学的依据。

然而，该领域仍面临高维扩展、非平稳性刻画、物理机制融合及数据不确定性等核心挑战。未来的研究应聚焦于发展高维高效建模技术、探索时变与非平稳相依结构，并致力于将数据驱动模型与地震物理规律相融合。特别值得关注的是，将离散的二元Copula关系整合为系统的参数影响力网络，以揭示主余震风险传递的关键路径，这有望推动研究从静态的“相关性描述”迈向动态的“机理解释”与“系统风险评估”。

参考文献：

1. [1] 赵银刚,刘庆杰,王晨,等.基于线性回归分析的主余震相关关系[J].地震地磁观测与研究,2017,38(02):71-76.
2. [2] 王元杰,赵晓芬,王君杰,等.基于方差主轴方法的地震动三分量相关性研究[J].地震学报,2025,47(02):254-268.
3. [3] 丁国,陈隽.序列型地震动物理随机模型研究[J].工程力学,2017,34(09):125-138.
4. [4] Shin M, Kim B. Effects of frequency contents of aftershock ground motions on reinforced concrete (RC) bridge columns[J].Soil dynamics and earthquake engineering,2017,97:48-59.
5. [5] Xu Y, Tang X S, Wang J, et al. Copula-based joint probability function for PGA and CAV: A case study from Taiwan region[J]. Earthquake engineering & structural dynamics,2016,45(13):2123-2136.
6. [6] Tesfamariam S, Goda K. Seismic performance evaluation framework considering maximum and residual inter-story drift ratios: Application to non-code conforming reinforced concrete buildings in Victoria, British Columbia, Canada [J].Frontiers in built environment,2015,1:1-18.
7. [7] 宋帅,钱永久,吴刚.桥梁系统地震易损性分析的混合Copula函数方法[J].工程力学,2017,34(01):219-227.
8. [8] Nelsen R B. An introduction to copulas (Second Edition) [M]. New York: Springer,2006.
9. [9] 刘亭亭,于晓辉,吕大刚.地震动多元强度参数主成分与结构损伤的相关性分析[J].工程力学,2018,35(08):122-129+137.
10. [10] Akaike H. A new look at the statistical model identification[J].IEEE transactions on automatic control,1974,19(06):716-723.
11. [11] Aho K, Derryberry D, Peterson T. Model selection for ecologists: The worldviews of AIC and BIC[J].Ecology,2014,95(03):631-636.
12. [12] Bradley B A.Empirical correlations between peak ground velocity and spectrum-based intensity measures[J]. Earthquake spectra,2012,28(01):17-35.
13. [13] Baker J W, Bradley B A. Intensity measure correlations observed in the NGA-West2 database, and dependence of correlations on rupture and site parameters[J].Earthquake spectra,2017,33(01):145-156.
14. [14] 李宁,刘洪国,刘平,等.近断层竖向地震动特征统计分析[J].土木工程学报,2020,53(10):120-128.
15. [15] 周正华,周雍年,卢滔,等.竖向地震动特征研究[J].地震工程与工程振动,2003,23(03):25-29.
16. [16] 韩建平,魏宏亮.基于不同地震事件数据的竖向地震动特性统计分析[J].土木工程学报,2013,46(增刊1):127-133.
17. [17] 王晓磊,王浠铭,阎卫东,等.基于Copula函数的水平和竖向地震动强度参数相关性分析[J].工程力学,2023,40(05):79-92.
18. [18] Boore D M. Orientation-independent, nongeometric mean measures of seismic intensity from two horizontal components of motion[J].Bulletin of the seismological society of America,2010,100(04):1830-1835.
19. [19] 施雨捷,芮佳,黄小宁,等.基于主余震理论的余震地震动PGA衰减关系研究[J].青海大学学报,2022,40(02):78-85.
20. [20] 温卫平.主余震地震动参数特征及损伤谱研究[D].哈尔滨工业大学,2015.
21. [21] 朱瑞广,吕大刚.基于Copula函数的主余震地震动强度参数相关性分析[J].工程力学,2019,36(02):114-123.
22. [22] 白治朋,文勇波,贺秀儒,等.基于条件谱匹配的主余震序列型地震动选择方法研究[J].世界地震工程,2025,41(02):21-31.
23. [23] 王霈宇.基于Copula函数的余震竖向条件谱研究[D].燕山大学,2024.
24. [24] 王泽龙.考虑主余震序列相关性的余震地震危险性分析[D].燕山大学,2023.