2000年姚安MS65地震震源及近邻区的低QS

利用震源距23 km范围内观测的2000年姚安MS65地震余震记录，计算了震源及近邻区域的尾波规一直达S波在频率15～20 Hz之间的傅里叶谱振幅.结果显示谱振幅随震源距增大而增大， 在对谱振幅进行了震源辐射方向性校正之后, 才出现谱振幅随震源距衰减的现象.由此获得了震源及近邻区域S波的Q(f)值,可表示为QS(f)=89f098其值比由尾波得出的姚安地区的平均QC(f)值低得多，表明了震源破裂带的强烈非均匀性对QS(f)的重大影响.     

昆仑山断层围陷波的分析和研究

对2001年昆仑山口西Ms81级地震产生的断层带,布设了沿断层和横跨断层的两条人工地震测线.通过对观测资料的定量分析和处理,求得了昆仑山断层带内部的细结构.分析工作包括从S波震相开始的振幅谱计算、速度频散计算、群速度测量,并用面波频散方法反演S波速度结构,用振幅谱比的方法估计断层带的Q值.野外试验结果表明,S波震相与围陷波组的时间差随炮点与台站之间距离增大而增加,在断层带外的测点上观测到与断层带相关的场地效应.最后得出昆仑山断层带宽度为250m、速度结构为断层内低速的分层结构和Q值为15（断层内）和30（围岩）.虽然昆仑山口西地震的震级比美国加州Landers地震的震级（Ms76）大,且地震产生的破裂带长度长得多,但是这两个地震断层带的宽度却相差不大.     

昆仑山MS81地震震源参数的多破裂段模拟研究

昆仑山MS81地震的已有研究结果在破裂带长度、破裂面方向、破裂面大小等震源破裂特征参数方面存在较大差异.本文采用D-InSAR技术首次获得昆仑山MS81地震干涉同震形变场，结合野外科学考察的实测值，进行了主破裂带InSAR视线向变化量的分解，通过对InSAR分解结果、野外科学考察、遥感解译等多源数据综合分析，重新划分了昆仑山地震的次级破裂段.进而通过对地震南北盘同震应变的分析，发现了昆仑山地震的南北两盘分别受挤压和拉张两种应力作用，研究表明多种岩石在拉张和压力作用下其最小主应力下的杨氏模量表现出非线性弹性特征，从而提出对昆仑山地震地表位移及震源特征参数分析时应考虑非线弹性介质导致的非线性弹性位移分布特征.基于上述原因，本文对Okada线弹性位错模型的算法进行了改进，提出了“多震源、非均一位错分量、多破裂段叠加”的线弹性模型，该模型模拟出的形变场干涉纹图较好地体现了地震形变场的分布特征，并由此获得了一套较为完整的地震发震断层的几何学特征参数，为破裂带长度、破裂面方向、破裂面大小等震源破裂特征参数研究提供了较好的解释.     

昆仑山MS81地震震源参数的多破裂段模拟研究

昆仑山MS81地震的已有研究结果在破裂带长度、破裂面方向、破裂面大小等震源破裂特征参数方面存在较大差异.本文采用D-InSAR技术首次获得昆仑山MS81地震干涉同震形变场，结合野外科学考察的实测值，进行了主破裂带InSAR视线向变化量的分解，通过对InSAR分解结果、野外科学考察、遥感解译等多源数据综合分析，重新划分了昆仑山地震的次级破裂段.进而通过对地震南北盘同震应变的分析，发现了昆仑山地震的南北两盘分别受挤压和拉张两种应力作用，研究表明多种岩石在拉张和压力作用下其最小主应力下的杨氏模量表现出非线性弹性特征，从而提出对昆仑山地震地表位移及震源特征参数分析时应考虑非线弹性介质导致的非线性弹性位移分布特征.基于上述原因，本文对Okada线弹性位错模型的算法进行了改进，提出了“多震源、非均一位错分量、多破裂段叠加”的线弹性模型，该模型模拟出的形变场干涉纹图较好地体现了地震形变场的分布特征，并由此获得了一套较为完整的地震发震断层的几何学特征参数，为破裂带长度、破裂面方向、破裂面大小等震源破裂特征参数研究提供了较好的解释.     

用随机模拟方法研究设定地震地面运动

用随机模拟方法研究设定地震地面运动可分两步进行：（1）采用中小地震的数字观测资料确定研究区路径、场地参数；（2）把确定的参数应用到研究区的强地面运动随机模拟中. 为此，本文首先采用了一种可行的方法，分离软基岩场地效应和非弹性衰减的影响，并把二者分别确定出来. 研究区北部中软基岩场地的平均效应在2～4Hz频率范围为15倍左右；研究区的S波品质因子为QS＝278f[KG*2]0.346. 把这些参数用于研究区的场地和路径模型中，并选择单拐角频率震源谱模型，随机模拟了研究区未来中强地震可能在北天山中段可能造成的地面运动，模拟加速度时程和反应谱可以服务于本地区的地震灾害预测和建筑物可靠性验算.     

电离层暴时的foF2预报技术研究

电离层是一个非常复杂的非线性系统.本文利用BP网络非线性输入-输出映射的特点预报电离层F2层临界频率.由于foF2在太阳活动高低年、不同经纬度、不同季节有着不同变化特点，本文仅以海口站、长春站的数据为例分别对1994年和2001年的foF2进行预报，预报结果表明发生电离层暴时用神经网络法预报foF2比国际参考电离层更有优势.此外，由于不同用户对预报精度的需求不同，本文以海口站2000、2001年的预报结果为例构建了残差模型.     

电离层暴时的foF2预报技术研究

电离层是一个非常复杂的非线性系统.本文利用BP网络非线性输入-输出映射的特点预报电离层F2层临界频率.由于foF2在太阳活动高低年、不同经纬度、不同季节有着不同变化特点，本文仅以海口站、长春站的数据为例分别对1994年和2001年的foF2进行预报，预报结果表明发生电离层暴时用神经网络法预报foF2比国际参考电离层更有优势.此外，由于不同用户对预报精度的需求不同，本文以海口站2000、2001年的预报结果为例构建了残差模型.     

2004年7月11日西藏Mw62级地震震源破裂过程研究

基于有限断层模型，利用远场体波波形数据研究了2004年7月11日西藏MW62级地震的震源破裂过程.结果表明该地震是一个以倾滑为主的浅源正断层型地震，震源深度为125km，断层面走向152°，倾角44°，平均倾滑角-117°.破裂在震中处成核，然后以28km/s的平均速度向两侧传播，在震中以东偏北5km处达到最大滑动43cm.该地震主张力轴近E W方向，受浅部NNW SSE或N S向裂谷带控制，青藏高原南部的逆冲运动是引发这次地震的直接原因.     

川滇及周边地区地壳横波衰减的成像研究

利用由149个台站记录的中国地震年报、四川和云南区域地震台网观测报告的5897个近震的27530条横波振幅与周期资料，反演得到川滇地区及周边的地壳衰减Q0值（频率为1 Hz时的Q值）空间分布图像、台站的场地响应以及由于介质衰减不同造成的震级偏差.结果表明，川滇及周边地区的地壳横波衰减Q0值在200～600范围内变化，平均值为400，横向变化显著.在川西-滇西北-保山地区、云贵交界-昆明-思茅地区以及海原-祁连、汾渭和东昆仑等活动块体边界带的Q0值明显较低，这些地区地震和构造活动相对活跃；而构造相对稳定的四川盆地、马尔康块体、桂西滇东断块以及鄂尔多斯、秦岭-大别等地区的Q0值则明显较高.此外，川滇地区大多数台站的场地响应都在1附近，仅有少数增益或衰减较大，达到2倍左右，且其分布无明显规律性.该区观测报告给出的大多数地震震级偏差不大于03，但也有少数震级偏差达04，其分布与介质品质因子Q0值分布相关，说明震级测定可能受地壳介质衰减横向变化的影响.     

川滇及周边地区地壳横波衰减的成像研究

利用由149个台站记录的中国地震年报、四川和云南区域地震台网观测报告的5897个近震的27530条横波振幅与周期资料，反演得到川滇地区及周边的地壳衰减Q0值（频率为1 Hz时的Q值）空间分布图像、台站的场地响应以及由于介质衰减不同造成的震级偏差.结果表明，川滇及周边地区的地壳横波衰减Q0值在200～600范围内变化，平均值为400，横向变化显著.在川西-滇西北-保山地区、云贵交界-昆明-思茅地区以及海原-祁连、汾渭和东昆仑等活动块体边界带的Q0值明显较低，这些地区地震和构造活动相对活跃；而构造相对稳定的四川盆地、马尔康块体、桂西滇东断块以及鄂尔多斯、秦岭-大别等地区的Q0值则明显较高.此外，川滇地区大多数台站的场地响应都在1附近，仅有少数增益或衰减较大，达到2倍左右，且其分布无明显规律性.该区观测报告给出的大多数地震震级偏差不大于03，但也有少数震级偏差达04，其分布与介质品质因子Q0值分布相关，说明震级测定可能受地壳介质衰减横向变化的影响.     

灵武地震尾波频率特性及介质参数的测定

本文提出了一种利用地方震尾波测定地壳介质相关长度和 Q 值的方法.该方法基于全面考虑了震源谱、介质散射特性和仪器响应诸因素对尾波优势频率的影响,得到公式t*=2/[1/f(1+12f2+1/f(1+22f2-1/f+1/2 d/df In I(f)] 其中t*=t/Q,f 为尾波优势频率,1,2分别为由介质相关长度和震源拐角频率确定的参数,I(f)为仪器响应.当给定震源参数(2)时,可通过对已测尾波 f-t 曲线反演计算求解相关长度和Q 值,应用这一方法,我们处理了灵武地区1982至1989年期间40余次地震的尾波资料.测定了相关长度和 Q 值,并对中强地震前后介质相关长度的时间变化作了初步研究.结果表明,震前介质相关长度表现出减小的迹象.      

唐山余震的尾波 Q 值研究r

利用唐山余震的数字化地震记录研究了唐山余震区的区域衰减特性,利用 Aki(1969)的单次散射理论,分析了滦县地区的尾波 Q 值及其随频率的变化.对于尾波的起始点,我们分别取1.5倍和2倍的 S 波到时进行了计算.两种取法都有强烈的 Q 值随频率的变化.在 Q=Q0(f/f0)的假定下,第一种情况 Q0分布在16-49之间,均值为28;值分布在1.04-1.51之间,均值为1.22.第二种情况 Q0在15——59之间,均值为34;值分布在0.96-1.43之间,均值为1.16.把第二种情况的关系绘入 Sato(1994)的 S 波和 S 尾波的 Q 值与频率的关系曲线中去,符合较好.通过比较不同作者在不同地区的结果,发现数据的离散不仅与所研究的地区有关,而且不同的尾波持续时间对结果亦有很大影响.我们发现低 Q0值对应于高值和较短的震中距,高 Q0值对应于低值和较长的震中距.这是 Q 值的深度依赖关系和频率依赖关系的综合效应.      

小近震三分向地震图S波记录振幅包络形状及Q值测定

本文采用三层地壳介质模型,讨论了地方台网记录的小近震三分向地震图S波段振幅包络形状问题.认为S波包络的规则形状主要是由于强反射界面上的一次和多次全反射波形成的.分析了地壳分层结构,震源辐射特性,介质衰减对n次全反射波振幅的影响.在此基础上,从理论上给出了从直达S波到尾波之间的S波段记录部分振幅包络形状函数和S波最大振幅.利用包络形状函数方法测定了宁夏地区介质品质因素Q值.      

青藏高原东北缘地壳S波速度结构与泊松比及其意义

利用甘肃地震台网16个台站记录的远震资料,采用最大熵谱反褶积方法,得到了各个台站的接收函数. 采用接收函数扫描法和线性反演方法对研究区的壳幔结构进行了研究,这两种接收函数方法得出的结果具有很好的一致性. 青藏高原东北缘地壳厚度变化剧烈,祁连块体为50～55 km、柴达木块体和河西走廊为45 km左右（合作台除外）,由北向南,Moho界面呈中央下凹的准对称状. 研究区地壳VP/VS介于166～185（σ=0215～0294,均值0254）,其均值接近或略低于全球平均值;S波速度结构可见壳幔过渡带具有明显的突跳,结合其他地球物理学证据,推断该区可能不存在岩浆底侵作用和地壳部分熔融现象. 该区地壳VP/VS值与地壳厚度呈反相关关系,推断该区地壳的主要组成成分以中酸性岩石为主,其45～55 km厚的地壳可能主要是通过上地壳的叠置形成的.     

不同散射模式下S波包络合成

为揭示地震波在地壳小尺度非均匀介质中的散射过程,更准确地描述地震波的包络展宽现象,本文基于多次各向异性散射理论,采用离散波数法求解能量密度积分方程,选取高斯型自相关函数表征的散射模式,得到S波能量密度包络。基于此,本文首先分析了单次散射和多次散射在形成S波能量密度包络中的贡献规律;然后探讨了吸收系数和总散射系数对合成S波能量密度包络的影响;最后对比了在不同散射模式下合成的S波能量密度包络的差异。结果显示:① 不同的散射模式下单次散射和多次散射对地震波散射过程的贡献规律是一致的,对于近震（震源距小于100 km）,单次散射模型可以近似合成S波能量密度包络;随着震源距增大,多次前散射模型可以更快地接近总能量密度包络;② 吸收系数增大会降低直达S波和尾波幅值,总散射系数增大会降低直达S波幅值,但使得S波尾波幅值升高;③ 前散射模式下S波能量密度包络随震源距的增大会导致峰值延迟,包络展宽,尾波衰减一致性更快等现象产生。      

青藏高原INDEPTH-Ⅲ剖面地壳厚度与泊松比：地质与地球物理含义

对INDEPTH Ⅲ台站的接收函数进行扫描，利用Moho界面产生的转换波和多次波的走时信息，估计台站下方的地壳平均波速比VP/VS和地壳厚度.结果显示：（1）沿着INDEPTH Ⅲ剖面，地壳厚度整体变化不大，均为65±5km，其中拉萨地块Moho界面埋深较羌塘地体要深约5～6km.结合其他研究资料，我们推断，在整个班公-怒江缝合带存在约10km的Moho错断，为拉萨地体北缘的地幔盖层向北俯冲到羌塘地体之下所致.（2）青藏高原地壳平均波速比整体都较高，可能与青藏高原地壳广泛存在的流体/部分熔融岩浆有关.拉萨地体北部异常高的地壳VP/VS可能与嘉黎-崩错右旋走滑断裂相关；而另一个泊松比异常区位于羌塘中北部（st36～st40），它可能是由热的地幔引起的壳内部分熔融所致.     

四川地震辐射能量和视应力的研究

利用成都数字地震台网的波形记录资料,研究了四川地区地震辐射能量〖WTHX〗E〖WTBZ〗s、σapp视应力等震源参数及相互关系.对时域波形数字记录首先进行数据处理,包括带通滤波,快速富氏变换,扣除仪器幅频特性;经反变换,回到时间域,给出数字波形时间序列,选取所有方向的整个p波和s波段数据;二次变换到频率域内计算速度和地动位移的功率谱积分,进而求得各台站记录的地震波辐射能量和视应力等震源参数.研究了2000年5月到2004年8月成都数字地震台网记录的224个3级以上地震的震源参数.四川地区地震辐射能量〖WTHX〗E〖WTBZ〗s与地震矩M0,应力降Δσ与视应力σapp,视应力σapp与震级ML大体成正比.其中地震视应力σapp与震级分布数据较宽,地震视应力与破裂尺度ra的关系也较分散.同时研究了四川地区近年视应力σapp的空间分布.方法是将计算的σapp值在一定空间邻域内作平均,从而绘出等值线.分析结果表明,σapp值的相对高值区分布在四川南部的川滇交界和北部的川甘交界一带.同时研究了四川地区地震视应力值的时间分布.最后对中小地震视应力σapp及其他震源谱参数的意义作了讨论.     

昆仑山断裂带围陷波的有限差分数值模拟解释

利用交错网格有限差分方法对昆仑山断裂带人工爆破产生的围陷波进行了三维数值模拟解释.为提高断裂带最终模型的可信度，在围陷波模拟的同时考虑了人工爆破记录的三个分量.对昆仑山断裂带围陷波的模拟结果表明，影响围陷波特性的断裂带深度主要在10 km以内.S波速度和断裂带宽度对围陷波的到时、频率、振幅和相位影响较大.数值模拟解释获得的昆仑山断裂带的细结构参数是：浅部断裂带宽度为300 m，深部为250 m；深度在400 m以上断裂带内S波速度为098 km/s，外部围岩S波速度为170 km/s，Q值为138；S波速度和Q值随着深度的增加而增加；1000 m以下断裂带内S波速度为280 km/s，围岩S波速度为33 km/s.     

新疆南天山西段地区介质衰减特征研究

利用新疆数字地震台网记录的15个台站的720条地震波形资料,用Atkinson方法联合反演得到南天山西段Q值随频率f的关系式为Q(f)=313.4f0.548;在400 km的震源距范围内,研究区地震波的几何衰减满足三段衰减模型;研究区的Q值小于天山中东段、北天山和柯坪地区的Q值,说明该区域构造活动强烈。其次根据S波衰减记录得到观测振幅谱,采用Moya遗传算法得到了新疆南天山西段15个台站的场地响应。结果显示,大多数台站的场地响应在1附近,西克尔台(XKR)、塔什库尔干台(TAG)、乌什台(WUS)和神木园台(SMY)的场地响应在1~4 Hz之间均大于2,表现出一定的放大效应,岳普湖台(YPH)呈现余弦波动,可能与台基岩性和仪器的工作状态有关。区域的选择对台站场地相应计算结果影响不大。     

2009年云南姚安M_S6.0地震余震序列尾波Q_c值研究

利用云南大姚地震台记录到的2009年姚安6.0级地震余震序列的数字化波形观测资料,基于尾波单次散射理论,估算该次地震震源区尾波Qc值。结果表明:中心频率为1.5 Hz时,振幅衰减率β(f)在-0.141 27～-0.043 05之间,平均值为-0.074 60;Qc值在33.3～99.8之间,平均值为65.5。并得出Qc值与频率的关系:尾波波源因子与震级成正比,与频率成反比。主震发生后头几天,尾波Qc值(1.5 Hz)波动较大,之后在±1.5倍方差线内趋于稳定,且没有超过5级的强余震发生。姚安地区尾波Qc值具有构造运动活跃、地震活动频繁地区的低Qc值特征。