地磁加卸载响应比与地震对应系关系研究

通过对浙、沪、苏邻近４个地磁观测点的近十年的资料,及邻近区域１０个地震进行了每日Ｚ分量日变加卸载响应比随时间变化和响应比月均值的相关研究。结果表明：加卸载响应比阈值ＰＯ（Ｚ）和响应比月均值ＰＭ（Ｚ）的大小与地震的震级、震中距及异常出现的早晚并非完全呈线性相关;当响应比值大于某一阈值时,有可能发生一定震级的地震;按该方法结果综合分析?载止１９９８年９月底,未来１年左右的时间内在该研究区域或更大范围（     

2021年漾濞MS 6.4地震云南数字化井水位同震响应特征分析

通过对漾濞M_S 6.4地震云南数字化井水位分钟值数据同震响应特征进行分析,分析观测仪器采样率对井水位映震能力的影响,并对水位同震响应机理进行初步分析。结果表明,对于漾濞M_S 6.4地震,13口同震响应观测井主要分布在震中距200 km范围内的红河断裂带北侧和小滇西地区;同震响应形态可分为阶升（上升）型、阶降（或缓降）型、阶变（阶升、阶降）—复原型等3种变化类型,阶升型同震响应幅度相对大于其他2种类型;由漾濞M_S 6.4地震烈度长轴和短轴共同划定的椭圆区域来划分此次地震静态应力和动态应力的影响范围,可以从各观测井受力状态初步分析水位同震响应变化形态的力学机理;在静态应力影响范围内,井水位同震响应主要表现为阶升变化;在静态应力影响范围以外,除以阶变（阶升、阶降）—复原变化为主外,还有阶降（缓降）、阶升变化。观测仪器采样率对井水位映震能力的影响较大,水位秒钟值数据的同震响应变化形态更加完整,但目前水位秒钟值数据的收集和使用都还在不便之处,这种情况有待改进。     

广域电磁法接收机数字信号FFT实时实现技术

针对广域电磁法（WFEM）接收机各频组数字信号快速傅里叶变换（FFT）点数不统一、低频组点数大不适合现场可编程门阵列（FPGA）实时实现的技术难题,本文采用插值和抽取相结合的抽样率转换方法,提出了满足广域电磁法接收机要求的等长度FFT变换技术,将各频组FFT变换点数统一变换到1024点,并以广域电磁法接收机实际采集到的电场信号进行了测试,最终在EP2C35F484C8 FPGA芯片上实现了1024点FFT变换,最高响应频率达到了93.84 MHz,完成1024点广域电磁信号FFT运算最快只需34 μs,并且能够有效提取广域电磁信号的振幅和相位信息.结果表明：在对广域电磁法接收机信号进行FIR低通滤波后,能够通过插值和抽取相结合的方法,将各个频组FFT变换的长度统一变换成1024点,实现广域电磁信号的FFT实时变换.本文提出的等长度FFT变换技术能够满足新一代高性能广域电磁法接收机FFT实时变换的需求.     

植被覆盖状况影响中国地表气温变化的观测事实

利用NOAA/AVHRR归一化植被指数（NDVI）及观测气温与再分析地表气温的差值（Observation Minus Reanalysis, OMR）分析了植被覆盖状况对中国地表气温变化的影响.结果表明,地表气温OMR趋势值与NDVI在空间上呈现出显著的负相关关系,植被覆盖状况差（NDVI小于0.1）的区域地表升温较为显著,气温OMR趋势值超过0.2℃/10a,而植被覆盖度高（NDVI大于0.5）的区域气温OMR趋势值则变化不大,甚至出现降温.气温OMR趋势值对植被的季节变化还有着敏感的响应.不同区域植被覆盖状况的差异可能导致中国地表气温变化对全球变暖的响应不同,预测中国未来气候变化需要考虑植被覆盖状况及其动态变化的影响.     

流域蒸散耗水率对气候和下垫面变异响应关系的稳定性研究

气候条件的变异和流域下垫面特征的改变是影响流域蒸散耗水的重要因素。本文聚焦于1900 2008年间全球83个典型流域数据,基于Budyko水热耦合平衡方程,探究100多年间不同条件下流域蒸散耗水率（AET/P）对气候和下垫面特征变异响应关系的稳定性。结果表明：(1)从长时间尺度看,大部分流域蒸散耗水率与气候干燥指数（PET/P）和流域特征参数（n值）变异的响应关系呈现较强的时间稳定性。从短时间尺度而言,半湿润流域内蒸散耗水率对干燥指数的响应系数?（AET/P）/?（PET/P）在20世纪内持续降低。不同气候条件下蒸散耗水率对流域特征参数的响应系数?（AET/P）/?（n）的变化差异显著。分不同下垫面特征来看,低n值（n<2）流域内AET/P对n值的变化更为敏感;(2)气候条件（PET/P）是大多数湿润区内蒸散耗水率的主导因素,在干旱与半干旱流域内,下垫面特征参数（n值）对AET/P贡献最大。在湿润区内,PET/P对AET/P的贡献程度随时间小幅提升;半湿润区内PET/P对AET/P的贡献度呈下降趋势。在低n值（n<2;流域持水能力较弱）流域内,n值对AET/P的贡献更多。在...     

高阶seislet变换及其在随机噪声消除中的应用

Seislet变换是一种小波类数学变换方法,主要根据不同小波级数上地震同相轴的局部倾角的不同来分析数据.一般意义上,测线方向上的离散小波变换（DWT）是一种特殊的零局部地震倾角的seislet变换.早期的工作基于低阶版本的离散小波变换来构建seislet变换,在本文中,通过使用Cohen-Daubechies-Feauveau (CDF) 9/7双正交小波变换（常用于JPEG2000压缩标准）作为框架,扩展高阶seislet变换方法.通过分析理论模型和实际数据的处理结果,并对比傅里叶变换、离散小波变换和低阶seislet变换,高阶seislet变换可以为地震数据提供更好的压缩比.因此更加适用于地震数据去噪处理.     

软土隧道基于地震响应的输入地震动排序

为探讨“最不利地震动”概念在地下结构抗震设计中的适用性,以软土地铁区间隧道为对象建立相应的地层-结构动力分析模型。以直径变形率为分析指标,基于动力时程方法研究18条不同输入地震动作用下隧道结构动力响应的分布及差异性,得出基于隧道地震响应的输入地震动排序,并通过调幅手段对比分析了地面峰值加速度（PGA）和隧道埋深变化对隧道结构地震动响应排序的影响规律。最后,评价了不同输入地震动参数,包括峰值加速度、峰值速度、峰值位移、绝对累积速度（CAV）和阿里亚斯（Arias）强度（I_A）与隧道地震响应之间的相关性。分析结果表明:① 随着PGA从0.5 m/s2增加到2 m/s2,地震动排序发生明显变化,并且不同输入地震动引起的隧道地震响应差异显著提高,最不利地震动引起的直径变形率与平均值的比值从1.1增加到1.9;② 隧道从浅埋到深埋的过程中,地震动排序结果基本保持不变;③ PGA为2 m/s2时,隧道地震响应与基岩面峰值速度（PBV）的相关性最好,相关系数达到0.94,其次是与基岩面峰值位移（PBD）和I_A,相关系数分别为0.62和0.48,相关性最差的是基岩面峰值加速度（PBA）和CAV,相关系数仅为0.37和0.22。研究结论可为今后软土隧道的输入地震动选择提供科学依据。      

全球地磁感应测深数据三维反演

全球地磁感应测深能获得地幔转换带及下地幔上部的导电结构.但目前稀疏的地磁台站分布及部分台站的观测数据稳定性较差,影响了三维反演对地下电性结构的分辨力和反演可靠性.为此,区别于传统的L2-范数反演方法,本文提出并实现了基于L1-范数的地磁测深响应三维反演技术.在反演中,利用L1-范数度量数据预测误差,降低"飞点"数据的影响,将相关系数较小的C-响应估计也纳入反演数据中.三维正演模拟采用球坐标系下的交错网格有限差分法,反演采用有限内存拟牛顿法.文中利用指数概率密度分布函数构造非高斯噪声的合成数据,并采用棋盘模型对反演方法的可靠性进行了验证.之后,我们将本文提出的三维反演方法用于全球129个地磁观测台站的C-响应数据反演,结果表明在地幔转换带深部,中国东北地区为高导电异常,南欧和北非则均为高阻;夏威夷在900km以下为高导;菲律宾海及以东地区在转换带表现为明显的高阻,这些结果与前人研究结果一致.由于采用了更多的台站数据,我们的反演结果还发现一些新的异常:南美洲南端,转换带表现为明显的高导;澳大利亚东南部,地幔转换带深部,也存在一个明显的高导异常,这些异常分布和地震层析成像的低速区一致.因此,L1-范数三维反演能够充分利用全球C-响应数据信息,提高地磁测深对地球深部电性结构的分辨能力,更好的研究全球地幔电性结构.     

海南岛介质非弹性衰减及场地响应反演

利用海南数字地震台网18个地震台274条波形资料,运用遗传算法反演海南岛非弹性衰减系数及台站场地响应,计算得到海南岛非弹性衰减Q值随频率f的关系为Q（f）=343.7f^0.42,说明海南岛适中Q值区域,地壳相对稳定。18个地表台站的场地响应场与频率有关,其中：7个台站的场地响应接近1,平稳变化;2个台站在高频部分衰减;7个台站的场地响应小于1,变化平稳; 2个台站低频部分场地响应接近1,高频部分放大。     

基于运营环境和提升小波变换的桥梁损伤检测研究

根据损伤桥梁在车辆荷载作用下的动力响应特点,以及提升小波变换对信号突变信息的放大功能,提出了利用桥梁运营荷载作用下加速度响应提升小波变换系数的分布特性对结构损伤进行识别的方法。首先,采集桥梁在行车荷载作用下的加速度响应信号;然后,对加速度响应信号进行提升小波变换,分别利用加速度响应信号、加速度响应信号差,提升小波变换系数空间变化的峰值识别损伤位置;最后,对行车速度、损伤位置、损伤程度和测量噪声对损伤识别效果的影响进行了分析讨论。结果表明：在行车速度8m/s以下、测量噪声不高于5%情况下,利用运营荷载作用下桥梁单点动力响应信号提升小波变换,可以实现桥梁多处损伤的检测和识别。     

18000年前冰期气候的数值模拟

了解过去环境变化的事实是准确评估当前环境所处的历史位置和预测未来气候变化趋势的前提．本文利用CLIMAP建立的末次冰期的地球表面状况作为下边界条件,用GCM模拟了18000年前冰期全球气候状况．模拟结果显示：冰期全球平均地表气温比现代低约4.9℃,自由大气温度从低层到高层也均有所下降,全球平均降水比现在的少0．25mm／d,相对湿度小10％,平均海平面气压值上升．本模式的模拟结果和其他GCM的模拟结果作了比较,并和地质考古得到的有限的冰期气候资料作了比较,两者基本相符．     

北京重力地球潮参数和重力仪力学特性的测定

根据北京白家疃测站８台ＧＥＯ潮汐重力仪对比观测的资料，采用标准线性体流变模型模拟重力仪及记录系统的力学特性，计算各台仪器的格值改正系数、主潮汐波的振幅衰减系数和相位滞后；确定北京重力地球潮的振幅因子δ和相让滞后为δ（Ｍ_２）＝１．１７８４，δ（Ｏ_１）＝１．１８０４，δ（Ｍ_２）／δ（Ｏ_１）＝０．９９８３，（Ｍ_２）＝０．０７°，（Ｏ_１）＝０．１０°．     

地震活动性的统计分析:由过去推测将来的可能性研究

为了由过去的地震活动性推测将来的地震活动性,引入了地震（震级≥m）的期望年发生率v（≥m）来描述一个地区的地震活动性．根据全球地震目录（1964-1994年）以及南加州（1932-1995年）和华北（1970-1994年）两个区域地震目录资料,以统计样本量作为目录记录时间长短的相对量度,对由不同的统计样本量计算得出的地震实际年发生率v（≥m,T,t）进行了统计分析,得到三点结沦：①在统计样本量n足够大的情况下,地震实际年发生率表现出准平稳时间过程的特征,可近似地看作地震期望年发生率,本文给出了这种近似的误差（离差系数）与统计样本量之间的定量关系;②离差系数与统计样本量之间的关系与震级无关,表现出不同震级层次的相似性;③统计样本量相同时,不同震级的地震期望年发生率之间满足logv（≥）=a－bm的关系,形式上与G－R关系相似,但它给出了由小地震的统计特征估计大地震的期望年发生率及其统计误差的方法．基于上述结论,进一步讨论了地震活动性的统计特征在地震危险性分析中的潜在应用。     

等离子体层突起位置的变化与对流电场的某些特性

本文分析了四种影响等离子体层顶突起位置的因素:（1）对流电场E_x分量的存在;（2）白天和夜间的E_y分量的不对称性;（3）偏心极光椭圆带;（4）E_y分量的突然变化。文章着重讨论了在平静时和在磁活动性增强时正E_x分量存在的可能性。文章肯定了Volland电场模型,但指出它存在某些缺点,建议做进一步改进。     

行星际激波特性的数值研究

采用二维二分量 ＭＨＤ模型,应用 ＰＰＭ格式,研究内日球赤道面内 １ＡＵ附近激波波阵面参数的分布特性,以及太阳风速、扰动源宽度、扰动源初始传播方向和行星际电流片对分布的影响．结果指出：（１）均匀背景情况下激波阵面西侧为强磁场区,东侧的动压差比西侧大,而密度比和气压比基本对称;（２）背景太阳风流速和扰动源初始传播方向对参量分布影响较大;扰动源宽度影响不大;（３）当扰动源中心处于电流片东侧或西侧时,密度比、气压比和磁场比峰值均西偏,而动压差峰值和扰动源处于同侧;当扰动源中心与电流片重合时,激波阵面东西两侧都出现峰值,为双峰分布,东侧动力学参数大,西侧磁场强．上述结论与观测资料的统计结果趋势上是一致的．     

华北地台北缘地球物理场特征与金属矿床预测

对华北地台北缘的地球物理场特征进行了探讨,根据重、磁异常数据反演计算了该地区的莫霍界面、居里界面、磁性界面的起伏.利用地球物理场资料和反演计算结果对该地区的构造格架和断裂进行了推断,同时预测了9个成矿远景区:(1)集宁一呼和浩特金成矿区;(2)张家口一赤城金、银多金属成矿区;(3)密云一高岭金、铜成矿区;(4)青龙一马兰庄金成矿区;(5)秦皇岛、金多金属成矿区;(6)郝家营多金属成矿区;(7)承德地区金、多金属成矿区;(8)赤峰一喀喇沁旗金成矿区;(9)宁城东金、多金属成矿区.     

論用PKP波測定极远地震的震級

一、引言目前我国测定极远地震（震中距离Δ>110°）震级的方法有两种: 1.面波法计算震级的公式 M_LH=log（A/T）-log（A/T）^*,（1.1） 式中A是水平分向地动位移（μ）,T是相应的周期（秒）,取面波中A/T为最大者,-log（A/T）^*是索洛维约夫-谢巴林1957年发表的量规函数。     

利用重力资料反演三维密度界面的质面系数法

文中给出一深度为h的水平矩形质面在地面产生的重力公式,对比了此公式算出的重力异常与三维长方体公式算出的重力异常,得出:当三维地质体的厚度与其埋深之比小于0.2-0.4时,二者的重力效应几乎相等,它们之间最大偏差在1-3％以内。用给出的重力公式进一步导出了重力反演三维密度界面的质面系数法。用理想模式检验,结果表明,本文给出的方法收敛快,反演误差较小,计算效果较好。最后用一实例,验证了方法的适用性。     

起伏地形上的位场转换

根据起伏地形上观测的异常场分量（沿任一方向的分量）换算出垂直与水平分量,本文提供了二维情况下的转换方法。其算法是:（1）通过共形映照,将w平面上的地形曲线L的上部映照为z平面的上半平面,将曲线L映照成z平面上的实轴。（2）求解边界为实轴的Hilbert边值问题。 通过理论模型的计算,验证了该方法的正确性。     

ENSO和长江大水对天文因子的响应研究

ENSO（厄尔尼诺和南方涛动的总称）对世界尺度的天气气候有影响,本文用概率论统计检验方法,查明长江大水与ENSO以及ENSO与天文因子的相关关系,并系统分析了三者之间的相关关系。结果表明:ENSO事件对天文物理主周期和天文特征有明显的响应,长江大水年对太阳活动特征、节气日的月相年变化和ENSO也有显著的响应关系。这些关系对长江大水的准确预测有重要价值。文中对天文物理因子对ENSO和大水的影响机理进行了探讨。