应力扰动对地震周期和地震矩影响的数值模拟

利用二维有限元数值模型,结合断层滑移弱化摩擦准则对断层滑动规律以及应力扰动对其影响进行了研究．数值计算结果表明,在均匀应力分布情况下, 平面断层滑动显示出典型的特征地震规律,断层面上的应力扰动对断层滑动规律产生影响,压应力增加明显延迟地震的发生时间,并增加地震释放的能量．应力扰动发生在地震破裂临界区时的影响比在震前滑移区时的影响显著．当发生在地震滑移区时,若应力扰动足够大,则压应力增大会造成地震发生时部分动力断层被暂时锁住,使得地震释放的能量变小,但可增加后续地震的能量; 而压应力减小则可导致地震规律产生更加复杂的变化,会即时触发地震．如果应力扰动发生在一个地震周期的早期,则触发的地震较小,但可导致随后的地震提前发生; 如果应力扰动发生在一个地震周期的后期,则会触发大地震．当应力扰动位于震前滑移区或破裂临界区时,小的扰动也可能产生类似的效果．应力扰动产生越晚,这种影响也越明显．应力扰动发生在破裂临界区的影响最明显．应力扰动的影响一般主要集中在应力发生扰动后的1—2个地震周期内．后续地震基本恢复无应力扰动时的特征地震规律．      

论速率—状态破裂模型和库仑破裂模型

考察了库仑破裂应力模型和速率－状态磨擦模型的预测效果。我们研究由于阶跃应力扰动（即地震引起的静态应力增加）以恒定速率在时刻t0时叠加到“背景”应力（如构造加载）上而引起的破裂时间的变化△t（时间提前）。△t的可预测性就意味着地震活动速率r(t)/r0的可预测性，r0是取决于构造应力的一恒定速率，地震活动速率可用地震目录来检验。符合余震序列一般特性的r(t)/r0模型必须能预测出符合大森定律的地震活动性衰减速率，序列的持续时间长度应少于主震轮回周期时间长度的百分之几，并能直接返回到背景活动速度。库仑模型要求一个断层在加载期间保持为闭锁状态，破裂瞬时发生，△t与t0无关。这些特性意味着地震活动速率会出现瞬时的无限增长，而持续时间为零。对不同的状态演变规则所作的r(t)/r0的数值计算表明，在主震发震时刻非常接近于破裂状态的断层上会发生余震，这些断层特性必定是“库仑式”的，从而使滑动演变规则可被排除。实际余震的数量特征可以制约速率－状态本构关系的参数；α也许会低于实验值，刚度也许会高，正应力也许会低于静岩压力。我们同时对库仑模型和速率－状态模型作了理论上的比较。当本构关系参量α相对于参量b减少时，速率－状态模型的断层特性就更接近于“库仑式”。这是因为滑动开始是减速的，表明断层接触状态初始是愈合的。当参量α数值较小时，减速更为明显，则断层表现为更加接近于闭锁状态。即使当速率－状态模型的△t具有库仑式特性时，它的大小仍然可能差－常数，常数大小与b有关。在这种情况下，速率－状态模型的特性类似于一个修正的库仑破坏模型，即弱化使破裂应力临界值降低，从而增加了时间的提前量。与非库仑式响应的偏离也与加载速率、弹性刚度、初始条件以及状态演变规则的假设等有关。     

断层演化过程中的周期和非周期地震滑移及非地震滑移一维断层模型分析

本文首先根据Dieterich和Ruina提出的含速率和状态的摩擦定律(Dieterich-Ruina定律), 基于一维弹簧-滑块模型推导了地震复发周期的解析表达式, 然后将该近似解与数值模拟结果以及Barbot等的相关研究进行了对比分析． 此外, 本文还利用数值模拟与理论分析研究了断层周期和非周期演化的力学成因机制以及非地震滑移形成的另类力学机制, 并讨论了一维弹簧-滑块模型的优点及其局限性． 结果表明: ① 震后滑移和自加速/成核阶段的持续时间在整个演化过程中不能被忽略; ② 在修正后的复发周期模型中, 复发周期的长短除了与断层特征尺度、 作用于断层面上的有效正应力和远场加载速率相关外, 还受Dieterich-Ruina定律中摩擦参数的取值以及临界滑移距离的影响; ③ 当给定各个物理参数和几何参数时, 目前所得到的解析近似解可以很好地估计地震的复发周期, 其相对误差可小于5%; ④ 在断层演化过程中, 施加剪切应力加载会产生非周期的地震滑移, 而在自加速/成核阶段后期或震后滑移阶段早期, 施加较大的剪切应力加载, 则会出现非地震滑移．      

数值模拟静态应力扰动下的断层失稳:结果分析兼与Dieterich模型和Coulomb模型的对比

基于大量岩石力学实验,Dieterich和其他研究者(Dieterich,1978;Ruina,1983)首先提出了描述岩石摩擦过程的速率-状态摩擦定律(R-S摩擦定律).如今R-S摩擦定律已成为研究地震成核等地震演化机制的有效手段.Dieterich(1992,1994)最早提出了描述受静态剪应力扰动后断层失稳时间提前或推后的余震触发机制的解析模型.现在Dieterich模型已经成为解释余震随时间衰减规律的Omori定律等地震观测现象的有力工具.与之相对应,广泛使用的Coulomb应力失稳模型也可以给出断层受到静态剪应力扰动后,断层失稳时间的提前和推后量.Dieterich模型和Coulomb应力失稳模型基于不同的物理方法,所以在进行地震危险性评估时,二者均有各自的局限性.本文利用R-S摩擦定律控制的1-D弹簧-滑块模型,模拟计算了理论地震循环以及在不同静态剪切应力扰动下,失稳时间的提前和推后量的变化情况,然后将计算得到的时间提前和推后量分别与Dieterich模型和Coulomb模型的相应计算结果进行了定量化对比和差异性分析,并给出了相应的解释.数值模拟结果显示,对于R-S摩擦定律在参数不变的条件下,断层模型失稳时间的提前和推后量的大小强烈依赖于静态剪应力扰动的大小和作用时间,而且绝对值相同的正、负向静态剪应力扰动造成的失稳时间的提前和推后量的变化情况并不完全一致.在震后松弛/滑移阶段和闭锁阶段,时间提前和推后量是常数,且随静态剪应力扰动绝对值的增大而增大,两者的比值接近于1.0,这与Dieterich模型和Coulomb模型的结果是一致的,相应的差值小于两模型结果的10%.而在自加速阶段,模拟计算结果则存在与Dieterich模型和Coulomb模型结果不同的特征.首先,在自加速阶段模拟计算结果均偏离Coulomb模型,而且时间提前和推后量的比值小于1.0,相异于Coulomb模型的论断.不过当受到正向静态剪应力扰动后,Dieterich模型的结果和模拟计算结果是一致的,最大相差量不超过Dieterich模型结果的7%,可接近0.对于负向静态剪应力扰动,当其绝对值较小时,Dieterich模型的结果很接近模拟计算结果,相差量不超过该结果的14%.但对于绝对值较大的情况,模拟计算结果远大于Dieterich模型的结果,最大可达Dieterich模型结果的35倍,这是由于负向静态剪应力扰动后使得Ω=δ·θ/dc1的条件不再成立,进而使得Dieterich模型不再成立.总的来说,与模拟计算相比Dieterich模型可以很好地描述1-D断层受扰后失稳时间提前和推后量的变化情况,并且可以体现出正、负静态剪切应力扰动后失稳时间提前量和推后量变化的差异性,而Coulomb模型则不能完整地给出受到静态剪应力扰动后断层失稳时间提前或推后的估计值.     

基于Dieterich地震发生率模型分析前震成因的力学机制

前震是地震前兆观测中的一个重要物理量,前震的研究对地震预测的发展尤为重要,前震序列的典型特征包括加速发生的小震与古登堡-里克特定律(Gutenberg-Richter Law)中的b值变小.本文在Dieterich提出的地震发生率模型的基础上,探讨前震发生的力学成因机制及前震-主震-余震时域演化特征.模型分析结果表明,前震震源区断层及其断层周边剪应力加载速率的变化可能是前震发生及地震发生率变化的一个关键原因.由于前震震源区次级断层之间剪应力加载速率受来自主震成核过程中断层自加速滑移的影响升高,从而可导致这些次级断裂的加速失稳,即加速前震的发生.当震源区内由前震所产生的静态应力扰动不可忽略时,应力扰动和主震成核的共同作用也可对后续前震发生率产生影响.当正向应力扰动出现时,后续的前震序列的地震发生率会出现陡增,随后其地震发生率逐渐下降.而当负向应力扰动出现时,地震发生率会出现陡降,然后再次逐渐上升.基于Kostrov模型,本文得到了剪应力加载速率与古登堡-里克特定律中b值的关系式,结果表明前震序列中b值的减小与前震区内的剪应力加载速率的上升有关.     

基于Dieterich地震发生率模型分析前震成因的力学机制

前震是地震前兆观测中的一个重要物理量,前震的研究对地震预测的发展尤为重要,前震序列的典型特征包括加速发生的小震与古登堡-里克特定律（Gutenberg-Richter Law）中的b值变小.本文在Dieterich提出的地震发生率模型的基础上,探讨前震发生的力学成因机制及前震-主震-余震时域演化特征.模型分析结果表明,前震震源区断层及其断层周边剪应力加载速率的变化可能是前震发生及地震发生率变化的一个关键原因.由于前震震源区次级断层之间剪应力加载速率受来自主震成核过程中断层自加速滑移的影响升高,从而可导致这些次级断裂的加速失稳,即加速前震的发生.当震源区内由前震所产生的静态应力扰动不可忽略时,应力扰动和主震成核的共同作用也可对后续前震发生率产生影响.当正向应力扰动出现时,后续的前震序列的地震发生率会出现陡增,随后其地震发生率逐渐下降.而当负向应力扰动出现时,地震发生率会出现陡降,然后再次逐渐上升.基于Kostrov模型,本文得到了剪应力加载速率与古登堡-里克特定律中b值的关系式,结果表明前震序列中b值的减小与前震区内的剪应力加载速率的上升有关.     

基于Chester-Higgs模型探讨摩擦生热对断层演化进程的影响

速率和状态相依赖的摩擦定律是本文采用的重要定律。结合Chester-Higgs摩擦模型和McKenzie-Brune摩擦生热模型,在一维弹簧-滑块-断层近似模型下,利用四阶变步长的Dormand-Prince算法,研究探讨了断层摩擦生热对断层演化的影响。结果表明:与忽略温度影响的情形相比,摩擦生热造成的温度上升可导致断层滑移时刻的略微提前,并伴随着摩擦系数和状态变量的下降,同时也使得断层的滑移量和应力降略有减小,而滑移速率有所增大;另外,在考虑温度影响时,有效正应力和临界滑移距离也会影响断层的演化过程,断层上的有效正应力越大,断层失稳时刻越提前,温度上升越明显;断层的临界滑移距离越大,断层失稳时刻则越迟,温度上升越显著,但当临界滑移距离超过5 cm时,具有不同临界滑移距离的断层,失稳时的温度则基本保持一致。      

汶川MW7.9地震余震序列触发机制研究

余震触发机制的Dieterich解析模型被广泛应用于区域地震活动性的定量分析以及依赖时间的概率地震预测模型的建立等方面. 基于滑移速率和状态相依赖的摩擦定律和弹簧-滑块模型, 从Dieterich断层滑移速率方程出发, 给出了静态应力扰动下触发地震的时钟提前或推后的近似解, 从而明确地阐明了触发地震的产生机制与断层的演化过程密切相关, 并与传统位错模型下库仑应力扰动时间提前或推后量作了比较. 采用对数线性拟合方法求得了汶川M_W7.9主震后余震序列持续时间, 符合Dieterich理论结果. 以汶川余震序列为例, 给出了两种不同的应力扰动模式在该余震序列中的应用. 结果表明, 经典Dieterich扰动解无法给出主震发生后即时余震数量的异常增加, 而考虑主震前后剪应力速率变化的Dieterich分段解则可反映出余震发生率及个数随时间的演化特征.      

汶川Mw7.9地震余震序列触发机制研究

余震触发机制的Dieterich解析模型被广泛应用于区域地震活动性的定量分析以及依赖时间的概率地震预测模型的建立等方面.基于滑移速率和状态相依赖的摩擦定律和弹簧-滑块模型,从Dieterich断层滑移速率方程出发,给出了静态应力扰动下触发地震的时钟提前或推后的近似解,从而明确地阐明了触发地震的产生机制与断层的演化过程密切相关,并与传统位错模型下库仑应力扰动时间提前或推后量作了比较.采用对数线性拟合方法求得了汶川Mw7.9主震后余震序列持续时间,符合Dieterich理论结果.以汶川余震序列为例,给出了两种不同的应力扰动模式在该余震序列中的应用.结果表明,经典Dieterich扰动解无法给出主震发生后即时余震数量的异常增加,而考虑主震前后剪应力速率变化的Dieterich分段解则可反映出余震发生率及个数随时间的演化特征.     

剪切载荷扰动对断层摩擦影响的实验研究

利用双轴伺服控制加载装置,采用3块花岗闪长岩标本组成的含有2个滑动面的直剪结构,开展了摩擦滑动实验。实验中通过在剪切方向上叠加正弦波状的位移扰动,研究了剪应力扰动对断层粘滑失稳的影响。研究表明,在恒定的正应力和位移速率下,标本表现出较为规则的粘滑;当在剪切方向叠加位移扰动后,随扰动振幅的增加,粘滑发生时间与扰动的相关性增大,粘滑的应力降和时间间隔的分布趋于离散,其中扰动能产生明显影响的临界振幅大致为0.05MPa;粘滑应力降和时间间隔随扰动振幅增大而逐渐离散的现象在较高正应力下更为明显,而相同扰动振幅下粘滑失稳发生时间与扰动的相关性也随正应力增加而增大;扰动周期对摩擦性状的影响不明显。研究结果意味着,地震引起的同震剪应力的变化不仅会引起邻近断层上地震发生时间的变化,也可能引起地震强度的变化     

三维速率—状态断层模型在地震轮回期内的非均匀库仑应力扰动

对遵从实验室得到的速率—状态摩擦定律的准动态连续三维断层模型进行了数值试验。我们做了一些测试,其中库仑应力（CS）或在整个断层增加,或仅在断层的局部增加。对于局部应力加载我们发现,如果触发几乎是瞬时的（1～2个月内）,成核发生在库仑应力增加区域的走向延伸部位。相反,如果在这几个月内地震没有发生,它可以在断层的任何地方成核,甚至比没有正的库仑应力加载要晚。这些特征意味着新的发现,其对三维模型而言是独特的,并且不能由一维弹簧一滑块模型来解释。该发现可应用在余震（依赖时间的）地震危险性评估中。     

关于Coulomb应力变化/扰动作用下的Dieterich余震触发机制的广义解

基于单自由度弹簧-滑块模型,滑移速率和状态相依赖的摩擦关系可用于对断层内部地震成核和断层失稳过程的定量化描述.Dieterich(1994)余震触发理论模型首次给出中强震后区域应力场受到静态应力扰动后所导致的区域地震活动性的时空变化特征,近期研究也表明Dieterich理论模型可进一步推广至依赖时间的地震预测模型的建立.本文从简单直观的断层群体化概念模型出发,推导出了包括静态剪应力和正应力扰动作用下广义Dieterich解.同Dieterich经典解相比,广义Coulomb应力变化:ΔCFF^G=Δτ－(μ₀－α)Δσ取代了Dieterich方程中原有的剪切应力扰动Δτ.从而表明余震发生率R同作用于断层上的正应力的变化(扰动)有着密切的相关性.进一步,我们讨论了传统Coulomb应力变化(扰动)模型在地震预测过程中可能存在的局限性.以1976年M_S7.8唐山大地震的主余震序列为例,采用本文中得到的结果,并结合视时窗分段方法,拟合了该地区地震活动性的时间演化过程.结果表明,除Coulomb应力变化(扰动)的影响外,主震前后加载于断层上的剪应力速率变化可对早期余震发生率产生很大影响.     

摩擦生热对断层演化过程的影响:Ruina-和Chester-Higgs-模型的对比分析

利用速率-状态摩擦定律(Rate-and State-Dependent Friction Law:简称RSF定律),结合McKenzie-Brune摩擦生热模型,本文分别从Ruina提出的RSF定律和Chester-Higgs提出的RSF定律出发,通过一维弹簧-滑块模型,采用四阶变步长的Dormand-Prince算法,对断层演化过程进行了数值模拟,探讨了摩擦生热对断层演化进程的影响.模拟结果显示,与Ruina-模型相比,Chester-Higgs-模型在断层高速滑动时存有更大的摩擦强度,表明摩擦生热对断层具有一定的强化作用,且同临界滑移距离的取值相关.而且,Chester-Higgs-模型在失稳时的断层面温度远远低于Ruina-模型,表明摩擦生热在断层演化过程中能抑制断层面温度的剧烈升高,且正应力和临界滑移距离越大,两种模型的温差越为明显,而断层的刚度和尺度则对温度的影响很小.模拟两种模型周期演化过程的结果表明,在相同的初始条件下,Chester-Higgs-模型给出的断层失稳周期明显比Ruina-模型更短,说明摩擦生热对断层自身演化最显著的影响是较大地缩短了地震重复发生周期.当断层进入周期性演化后,Chester-Higgs-模型给出的摩擦强度大于Ruina-模型,且对第一次非周期性失稳的摩擦强度和剩余应力的继承性更好.另外,由ChesterHiggs-模型给出的静态应力降远小于Ruina-模型给出的结果,所对应的单个事件的滑移量也小于Ruina-模型.     

潮汐应力对2007年M_s6.4宁洱地震震源断层成核失稳过程的影响

在潮汐库仑破裂应力计算和相关性检验的基础上,采用速率-状态依赖性摩擦本构关系,进行潮汐应力对宁洱地震震源断层成核失稳过程影响的数值模拟研究,并讨论潮汐正应力和潮汐剪应力在潮汐应力触发中的作用,以及应力的振幅和周期振荡性对潮汐触发的影响.结果发现,宁洱地震与潮汐应力之间存在潮汐触发相关性.宁洱地震震源断层的成核过程存在一个转变周期,当潮汐应力周期接近和大于转变周期时,潮汐应力的周期性将直接影响潮汐触发作用的效果.潮汐应力振幅值的增加也有助于潮汐应力对断层地震成核失稳过程的影响,正应力和剪应力振幅值的影响效果差别不大,但是潮汐剪应力比正应力能更迅速的起到潮汐触发作用.潮汐应力触发地震的过程中,主要是正值的潮汐正应力和剪应力对地震起到触发作用,而负值的潮汐应力对断层的影响小于正值的潮汐应力,断层对负的潮汐应力不敏感.     

正应力扰动对断层滑动失稳影响的实验研究

利用双轴伺服控制加载装置,采用三块花岗闪长岩标本组成的含有两个滑动面的直剪结构,开展了摩擦滑动实验.实验中通过在垂直滑动面的载荷上叠加正弦波状和方波状的扰动,研究了正应力扰动对断层黏滑失稳的影响.研究表明,在恒定的正应力和位移速率下,标本表现为规则的黏滑,叠加正应力扰动后,随扰动振幅的增加黏滑发生时间与扰动的相关性增大,黏滑应力降和时间间隔的分布趋于离散.黏滑应力降和时间间隔的平均值随平均正应力的增加呈线性增长,扰动叠加后黏滑应力降的离散度随平均正应力的增加而增大;黏滑应力降和时间间隔主要受应力变化幅度的影响,而与应力变化的速率关系不大.剪应力和正应力扰动都会对断层黏滑失稳产生影响,而正应力扰动的影响更明显.这两种扰动对断层黏滑失稳影响的机制存在差异,剪应力扰动只是改变断层滑动的推动力,而正应力扰动则改变了断层面上凹凸体的接触状态.     

局部山体地形对断层动力学破裂过程的影响研究

基于曲线网格有限差分方法研究了垂直走滑断层在不同山体地形情况下的动力学破裂传播,分析并讨论了局部山体地形对断层破裂过程及相应地面地震动的影响,得到了各模型断层面的动力学破裂过程及相应的地表峰值速度特征。研究结果表明,山体地形尺度（山体高度及底部延展距离等）对断层动力学破裂过程影响较大,进而影响到相应的地面地震动分布。当山体地形处于自由地表上亚剪切向超剪切转换的位置附近时,山体地形会阻碍断层面上自由地表超剪切的产生。一般而言,对于具有一定埋深的断层,当山体地形底部延展距离一定时,山体高度越高,其对自由地表超剪切的阻碍程度越大;当山体高度一定时,地形底部延展距离越大,越会阻碍自由地表超剪切的产生,这种破裂过程的变化会导致相应地面地震动呈现不同特征的分布。此外,还探讨了断层破裂过程及相应地震动对成核区外初始剪切应力变化的响应,结果显示,当初始剪切应力较高时,高应力降引起的超剪切破裂会对断层破裂及相应的地震动分布起主导作用。      

双向压缩条件下非均匀断层标本变形过程中的声发射活动特征

在双向压缩条件下,研究了一种非均匀断层模型标本在变形过程中的声发射(AE)时空特征.结果表明,随侧向应力sigma;2的增大,断层破裂强度提高,自加载至破坏的时间延长,破坏形式由突发失稳逐渐转变为渐进式破坏.预制断层及其非均匀性对标本变形过程中AE空间分布图象起明显的控制作用.预制断层,特别是预制断层上的强度不均匀部位及高强度段落,控制着微破裂空间分布的总体格局.破裂局部化均起始于断层上强度不均匀部位,随侧向应力sigma;2的增大,微破裂密集区域从预制断层强度不均匀部位逐渐扩展到整个高强度区段.侧向应力sigma;2的大小对AE时序特征具有显著影响.较低sigma;2条件下,断层表现为突发失稳,失稳发生在微破裂活动ldquo;增强-平静rdquo;的背景之上;而较高sigma;2条件下,断层表现为渐进式破坏,破坏前后AE频次持续增加,AER呈指数增长.sigma;2对b值变化也具有一定影响.当断层破坏表现为突发失稳时,b值在弹性阶段后期至弱化阶段显示出ldquo;前兆rdquo;性降低;而当断层表现为渐进式破坏时,b值变化平稳.      

断层上的地震成核过程与前兆模拟研究

为了把握断层上地震的前兆规律,对连续介质中断层面上的地震成核过程进行了数值模拟研究。摩擦采用了速率及状态依赖怀本构关系,并考虑了正应力变化的影响。结果表明,在成核的后期阶段：①最大位移速率单调加速增加;②滑动热点（最大滑动速率处）在后期阶段开始自发性迁移,且在空间上③当最大滑动速率达到可以明显探测的量级时（高于加载速率一个数量级以上）,倒刘时间为２０小时或更长一些,这时的速率变化可作为临震预测标志     

摩擦状态-速率依从的区域地震触发模型研究

设计一个具有摩擦状态-速率依从机制模拟大地震后区域触发地震的模型, 该模型由一定数量表示地震成核断层的函数组成。 大震对区域静应力场的影响由应力阶跃载荷表现, 利用摩擦状态-速率依从的一维自由度弹簧滑块模型结果计算地震发生时间提前。 结果表明, 施加幅度为地震应力降几十分之一至百分之一的应力阶跃, 可以引起地震发生速率的明显变化, 该变化通过统计检验, 与大地震后的区域触发地震活动相似。     

龙门山断层地震周期及其动力学过程模拟研究

在断层面上引入速率-状态相依摩擦本构关系、考虑铲形逆冲断层几何结构特征、断层下盘和上盘中下地壳及上地幔为黏弹性介质、上盘上地壳为弹塑性介质,本文用二维有限元动力学模型模拟了龙门山断层上大震准周期复发行为、分析了断层上地震孕育位置、地震周期不同阶段的应力/应变场演化特征.不同于近垂直走滑断层上的地震周期行为,大陆铲形逆冲断层上的构造应力的积累和释放过程更复杂、有其独特性.我们得到如下认识：（1）铲形逆冲断层上的地震复发是准周期行为.（2）龙门山断层最大库仑应力位于断层17~20 km深处,应力长期积累和同震释放都在此深度最大,说明地震会在此处孕育、发动.（3）在断层破裂的深部和浅部,同震滑动大小和构造应力释放大小并非同步,而是差异悬殊.（4）地震仅部分释放区域积累的应变能,断层上盘上地壳顶部和底部的褶皱、破裂等永久变形形式也是释放应变能的重要形式.（5）应变能密度增量的演化图像分为：震间、同震、震后期,清晰反应了龙门山断层附近的地震动力学过程.（6）地震发生除释放能量外,同时也对近断层的中地壳和断层底部有很大的应变能加载;这些加载,在震后期可能通过震后滑移、余震或中下地壳乃至上地幔的驰豫形变用几十年时间释放.以上对大陆内铲形逆冲断层上变形特征的了解,有助于我们在其地震周期行为中评估地震危险性.