基于FFT-MA的随机介质建模方法

随机介质是描述地球介质小尺度非均匀性的有效模型,在地震散射波场分析、储层描述等领域具有广泛应用,快速准确的随机介质建模方法是开展相关研究的基础与前提.本文首次将FFT-MA（Fast Fourier Transform Moving Average）算法引入到随机介质建模研究中,分析了该方法与传统随机介质建模方法相比具有的优势,并提出了基于FFT-MA的非平稳随机介质建模方法.建模实验表明,与传统的基于谱分解定理的随机介质建模方法相比,基于FFT-MA的方法在空间域产生随机数序列,使随机数序列与结构参数分离,因此,随机数序列与所建模型存在空间上的对应关系,可以分区域建模和局部修改模型.在非平稳随机介质模型建模时,滑动窗口的大小可以根据自相关长度变化而变化,避免了每个采样点都建立一次完整大小的模型,提高了建模效率.因此,FFT-MA随机介质建模方法能准确构建满足自相关函数要求的平稳及非平稳随机介质模型,具有建模效率高、灵活、实用的优点.     

Evaluation of nonlinear models for precipitation forecasting

Prediction of factors affecting water resources systems is important for their design and operation. In hydrology, wavelet analysis (WA) is known as a new method for time series analysis. In this study, WA was combined with an artificial neural network (ANN) for prediction of precipitation at Varayeneh station, western Iran. The results obtained were compared with the adaptive neural fuzzy inference system (ANFIS) and ANN. Moreover, data on relative humidity and temperature were employed in addition to rainfall data to examine their influence on precipitation forecasting. Overall, this study concluded that the hybrid WANN model outperformed the other models in the estimation of maxima and minima, and is the best at forecasting precipitation. Furthermore, training and transfer functions are recommended for similar studies of precipitation forecasting.     

宾川地区的地壳厚度和泊松比分布特征研究

This paper presents the changes of crust thickness and Poissons ratios distribution in the Binchuan region,where the first air-gun transmitting station and its a small dense array were deployed.From September 2011 to January 2014,more than 239 teleseismic events of M≥6.0 were recorded in 16 stations in the Binchuan region.Their P-wave receiver functions were analyzed respectively.The high spatial resolution result shows that the average crust thickness of Binchuan region is 45.3km,it follows the rule of"deeper in the north and east part,shallower in the south and west part."The deepest region is in Xiaoyindian Station;the crust thickness is 47.9km;the shallowest region is in Paiying Station,it has the thickness of 42.1km.It shows that the deeper Moho surface nearby the Chenghai fault and shallower nearby the Honghe fault;the isoline distribution of thickness changes greatly nearby the Chenghai fault and slowly nearby the Honghe fault.From the distribution of Poissons ratios,it is unevenly in the study area with a great difference from the north part to the south part,which shows a characteristic of"lower in the north,higher in the south".The Poissons ratio nearby the Honghe fault is medium too high(0.26≤σ≤0.29);lower nearby the Chenghai fault(σ≤0.26).This paper concludes the possible reason of different characteristic between Poissons ratio and crust thickness is thicker in the upper crust in the Binchuan region.     

接收函数确定Moho面速度和密度跃变的方法研究及应用

通过对单层模型反射和透射系数的推导,提出了利用接收函数一次转换波和多次波确定Moho面速度和密度跃变的速度-密度跃变(δβ-δρ)扫描叠加方法.利用反射率法计算了不同模型的远震理论地震图,按照与处理实际观测波形一致的方法和流程计算了理论接收函数;根据不同模型数值试验结果,深入分析了界面速度和密度跃变对接收函数震相幅度的影响.利用(δβ-δρ)扫描叠加方法,对理论接收函数进行了数值试验,结果证明了该方法的可行性.最后将该方法应用于位于青藏高原东北缘的高台(GTA)台和兰州(LZH)台,确定了两个台站下方Moho面的速度跃变分别约为(19±1)%和(20±1)%,密度跃变最小值为(4±2)%和(6±2)%.     

地壳形变分析的抗差最小二乘配置迭代解法

利用最小二乘配置进行地壳形变分析,其结果的合理性关键在于经验协方差函数的拟合.考虑到观测数据存在粗差的情况,提出基于观测值中位数初值的抗差最小二乘配置方法和基于中位参数法的抗差最小二乘配置方法.两种方法首先分别利用观测值中位数给出观测值初始权阵以及利用中位参数法给出最小二乘配置初始解,然后均在给定协方差函数参数初始值的情况下,应用合适的等价权进行抗差估计并通过迭代计算,最终获得稳健的协方差函数参数估值及最小二乘配置解.利用本文提出的两种方法以及传统方法分别对庐山地震的GPS垂直位移数据和意大利L'Aquila地震的InSAR同震位移数据进行处理分析.结果表明:相对传统方法,基于观测值中位数初值的抗差最小二乘配置方法效果更好,更具稳健性.     

利用近震高频接收函数研究四川理县西山村滑坡体结构

远震接收函数已广泛用于反演台站下方的结构,然而由于地球的非弹性衰减作用,远震数据较难获得高频接收函数,对浅地表结构约束不足.为了克服这一问题,我们使用近震数据的高频接收函数来研究浅表速度结构,并应用于四川理县西山村滑坡体上3个宽频带地震仪记录到的近震事件.本文发展了接收函数V_P-k(V_P/V_S)叠加方法,结合接收函数H-k叠加和波形反演方法获得了台站下方滑坡体的厚度、S波速度和平均V_P/V_S比,并与钻孔得到的滑坡体厚度进行对比.结果表明,滑坡体具有小尺度的横向不均匀性,台站下方滑坡体的平均V_P/V_S比在2.4~3.1之间变化并且在底层存在78~143m·s~(-1)左右的S波低速层.本文观测到的高V_P/V_S比和底层低的S波速度结构,与电磁法获得的滑坡体底层低的电阻率和底部富水特征一致,表明滑坡体h1底界面的抗剪强度相对较弱,是潜在的滑坡危险区域.本文研究结果表明,利用近震接收函数能有效约束浅表的速度结构,进而能为滑坡灾害治理提供一定的地震学参考.     

羌塘盆地中部地区地壳S波速度结构及构造意义

羌塘盆地是我国最大的海相盆地.本文根据在羌塘盆地内布设的27个宽频带地震观测台站记录的远震波形数据,利用非线性复谱比反演算法得到各台站下方100 km深度范围内S波速度结构.结果表明.羌塘地区Moho深度较为平缓,平均深度为61 km;北羌塘地壳内低速层广泛分布;北羌塘具有两个较大的沉积盆地,龙尾错和白滩湖坳陷,沉积厚度分别有10 km和15 km.尽管北羌塘下地壳受到强烈的新生代火山岩作用改造,但是这种深部岩浆热作用会加速烃源岩中有机质的热演化历程.北羌塘两个盆地具有很好的油气前景.与北羌塘低速层分布相比,南羌塘下低速层更深,可能与班公怒江洋于中生代的俯冲消减及拉萨地体北向俯冲有关.     

利用远震接收函数揭示的喜马拉雅东构造结台阵下方地壳结构及其动力学意义

在喜马拉雅造山带的东缘,雅鲁藏布江缝合带在这里发生急剧转折,南迦巴瓦变质体快速隆起,然而关于东构造结的形成机制一直未有定论.利用围绕南迦巴瓦峰的48个宽频带地震台站记录的远震数据提取P波接收函数,采用改进的H-κ叠加方法和共转换点叠加方法综合研究了东构造结的地壳厚度、波速比分布和地壳结构特征.结果表明:研究区平均地壳厚度为64.03 km,大部分台站介于60.48~66.55 km范围;平均波速比为1.728,主要集中范围为1.696~1.742.东构造结地壳厚度横向变化剧烈,构造结西端和北端厚而中间薄,东构造结核部Moho面呈现上隆的构造形态,东西向上隆幅度约为6~7 km,南北向的上隆超过9~10 km.东构造结核部地壳上隆减薄可能由高密度、高波速的岩石圈撕裂残片拆沉到上地幔软流圈后重力失衡所致.平均波速比超过1.8的高值异常展布于东构造结的两侧,推测为环东构造结的壳内部分熔融体.东构造结地壳上隆减薄和壳内部分熔融的存在很可能均与幔源热物质的上涌有关,而软流圈地幔的上涌则可能由印度板片的撕裂引起.     

中国东北地区高分辨率地壳结构:远震接收函数

利用分布在东北地区的国家地震局台网、NECESSArray台网、吉林大学在长白山及其周边地区布设的26个临时台站总计259个台站接收到的16,070条高质量的P波接收函数,采用H-k和CCP(Common Conversion Point,共转换点)叠加成像方法,获得该区高分辨率的地壳结构.观测结果显示,东北地区莫霍界面深度和地表高程总体呈镜像关系;西部大兴安岭—太行山重力梯级带附近存在莫霍界面深度陡变带;中部的松辽盆地地区受晚中生代的地壳拉伸作用影响,地壳厚度较薄,北部的小兴安岭地区和南部的华北北缘造山带可能同样受拉伸运动影响,具有较小的地壳厚度;松辽盆地莫霍界面深度由西向东逐渐减小,推测这与太平洋板块俯冲作用有关;东部地区莫霍界面结构比较复杂,依兰—伊通断裂与敦化—密山断裂之间出现复杂震相,可能与该区存在地幔物质的底侵作用有关;长白山火山地区地壳厚度较大,对应较高的波速比,推测在该区地壳内存在岩浆囊.     

低频检波器频谱补偿算法建模及仿真

采用数值仿真的方法建立了动圈式检波器的仿真模型,推导检波器传递函数,通过对实际动圈式检波器频率特性的测试,验证了动圈式检波器数值仿真模型的正确性,分析动圈检波器频率响应特性,分别对动圈式检波器的自然频率、灵敏度、阻尼比进行仿真,对不同自然频率、灵敏度、阻尼比对检波器输出频率特性的影响进行了分析,使用二阶二次型反馈补偿方法,设计检波器低频补偿结构,结果表明:数值仿真模型与实际检波器的频率响应一致,补偿结构使检波器自然频率由2 Hz降低至0.2 Hz.