复杂介质地震定位中震源轨迹的计算

在地震定位中常常需要求解震源轨迹,但由于复杂介质中的震源轨迹较为复杂,难以给出其解析解,因此震源轨迹的计算通常仅限于简单介质模型.本文基于最小走时树射线追踪技术,提出了一种计算复杂介质中震源轨迹的方法.为回避发震时间问题,以观测到时差作为震源轨迹的约束条件.首先从模型节点中选出少量理论到时差与观测到时差之绝对差,即双重时差较小的点作为震源轨迹的代表点,然后以其中双重时差最小的点为初始点,在双重时差场中利用最小走时树射线追踪方法计算出初始点到其他震源轨迹代表点的射线路径作为震源轨迹.当选的震源轨迹代表点较多时,得到的震源轨迹较为粗略,此时可去掉射线经过次数较少的代表点的射线路径使震源轨迹更为精细.为减少计算量,对最小走时树射线追踪方法的终止条件做了修正.以一个复杂介质模型中的地震为例,计算了包括速度扰动、到时扰动等不同情况下的震源轨迹,结果表明所提出的震源轨迹计算方法切实可行.     

南海东部海盆扩张过程的数值模拟

南海的形成演化一直是国内外关注的热点之一,南海的扩张发生在早渐新世-早中新世,并且在南海的扩张中至少包括一次洋脊跃迁.本文采用上升离散地幔流和热对流耦合模型模拟了南海东部海盆扩张、洋中脊跃迁和扩张停止后岩石圈的热结构演化和岩浆熔融过程.数值模拟结果表明,洋中脊跃迁是南海扩张中的一个重要的过程,由于洋脊跃迁使得在两洋脊之...     

Linear Correlations between Peak Frequency of Gyrosynchrotron Spectrum and Photosphere Magnetic Fields

The gyrosynchrotron spectra are computed in a nonuniform magnetic field case, taking into account the self-and gyroresonanceabsorption. It is found that the peak frequency νp of the gyrosynchrotron spectrum systematically increases with the increasing photosphere magnetic field strength B0 and increasing viewing angle θ. It is also found for the first time that there are good positive linear correlations between νp and B0, and between log νp and log θ, with linear correlation coefficient 0.99 between νp and B0 and 0.95 between log νp and log θ. We apply the correlations to analyze two burst events observed with OVSA and find that the evolution tendencies of the photosphere magnetic field strength B0 estimated from the above expression are comparable with the observational results of SOHO/MDI. We also give a comparison of the diagnostic results of coronal magnetic field strength in both uniform and nonuniform source models.     

Maximum flux density of the gyrosynchrotron spectrum in a nonuniform source

The maximum flux density of a gyrosynchrotron radiation spectrum in a mag- netic dip|oe model with self absorption and gyroresonance is calculated. Our calculations show that the maximum flux density of the gyrosynchrotron spectrum increases with in- creasing low-energy cutoff, number density, input depth of energetic electrons, magnetic field strength and viewing angle, and with decreasing energy spectral index of energetic electrons, number density and temperature of thermal electrons. It is found that there are linear correlations between the logarithms of the maximum flux density and the above eight parameters with correlation coefficients higher than 0.91 and fit accuracies better than 10%. The maximum flux density could be a good indicator of the changes of these source parameters. In addition, we find that there are very good positive linear correla- tions between the logarithms of the maximum flux density and peak frequency when the above former five parameters vary respectively. Their linear correlation coefficients are higher than 0.90 and the fit accuracies are better than 0.5%.     

Modelling the response of shallow groundwater levels to combined climate and water-diversion scenarios in Beijing-Tianjin-Hebei Plain,China

A three-dimensional groundwater flow model was implemented to quantify the temporal variation of shallow groundwater levels in response to combined climate and water-diversion scenarios over the next 40 years (2011–2050) in Beijing-Tianjin-Hebei (Jing-Jin-Ji) Plain, China. Groundwater plays a key role in the water supply, but the Jing-Jin-Ji Plain is facing a water crisis. Groundwater levels have declined continuously over the last five decades (1961–2010) due to extensive pumping and climate change, which has resulted in decreased recharge. The implementation of the South-to-North Water Diversion Project (SNWDP) will provide an opportunity to restore the groundwater resources. The response of groundwater levels to combined climate and water-diversion scenarios has been quantified using a groundwater flow model. The impacts of climate change were based on the World Climate Research Programme’s (WCRP’s) Coupled Model Intercomparison Project phase 3 (CMIP3) multi-model dataset for future high (A2), medium (A1B), and low (B1) greenhouse gas scenarios; precipitation data from CMIP3 were applied in the model. The results show that climate change will slow the rate of decrease of the shallow groundwater levels under three climate-change scenarios over the next 40 years compared to the baseline scenario; however, the shallow groundwater levels will rise significantly (maximum of 6.71 m) when considering scenarios that combine climate change and restrictions on groundwater exploitation. Restrictions on groundwater exploitation for water resource management are imperative to control the decline of levels in the Jing-Jin-Ji area.     

高温下甲基丙烯酸正丁酯ATRP的合环失活

以2-溴代异丁酸乙酯[2-(E iB)-B r]为引发剂,CuB r/2,2’-bpy为催化体系,进行了甲基丙烯酸正丁酯(n-BMA)在异丙醇(IPA)/水中的原子转移自由基聚合,在16.5~100℃范围内研究了温度对聚合反应活性特征的影响。动力学结果表明,低温时反应有良好可控性,温度升至85℃以上时,聚合反应呈现出一系列非活性特征,聚合进行到一定程度后转化率不再增加,所有的链都成为死链。红外光谱的结果证明,较高温度时处于休眠态的增长链发生了合环反应,导致了增长链的永久失活。     

基于姿态参数的海洋可控源电磁数据一维反演

基于姿态信息的海洋可控源数据一维反演采用Occam方法进行,反演数据为不同收发距离观测到的电场振幅(MVO)。接收器和发射器姿态参数同时引入到反演中,正演过程引入发射源姿态参数,借助矢量分解将任意姿态的发射源分解为3个正交坐标轴分量源,分别计算接收点处电场响应,叠加得到总电场响应;根据接收器姿态信息,将计算坐标系中的电场矢量进行欧拉旋转获取任意姿态时测量到的电场响应。理论模型与实际数据反演结果表明,基于姿态信息进行海洋可控源勘探数据反演是可行的。     

内蒙古东乌旗地区崇根山岩块超基性岩地球化学特征

崇根山岩块是贺根山蛇绿岩带内超基性岩出露面积最大的岩块,岩性主要以二辉橄榄岩为主,方辉橄榄岩次之以及少量纯橄岩透镜体或团块。岩石学和地球化学特征表明这套超基性岩属于位于蛇绿岩套中下部的地幔橄榄岩序列,其二辉橄榄岩和方辉橄榄岩具有富Al的特征,分异指数很低(DI=0.1~2.79),而固结指数很高(SI=79.9~89.62),代表了原始地幔的残余。稀土元素表现为LREE弱亏损的特征,与阿尔卑斯型辉橄岩相似;纯橄岩稀土含量仅为球粒陨石的1/10,具有"U"型稀土元素分布模式,与阿尔卑斯橄榄岩和SSZ型蛇绿岩地幔橄榄岩相似。崇根山大面积出露的二辉橄榄岩可能代表了洋盆扩张早期的地幔橄榄岩的残余。整体上富Sr贫Ti的现象暗示了后期消减作用对地幔橄榄岩地球化学结构的影响。     

横向各向同性介质中地震波走时模拟

横向各向同性介质是地球内部广泛分布的一种各向异性介质.针对这种介质,我们对各向同性介质的最小走时树走时模拟方法进行了推广,推广后的方法可适用于非均匀、对称轴任意倾斜的横向各向同性介质模型.为保证计算效率,最小走时树的构建采用了一种子波传播区域随地震波传播动态变化的改进算法.对于弱各向异性介质,我们使用了一种新的地震波群速度近似表示方法,该方法基于用射线角近似表示相角的思想,对3种地震波（qP, qSV和qSH）均有较好的精度.应用本文地震波走时模拟方法对均匀介质、横向非均匀介质模型进行了计算,并将后者结果与弹性波方程有限元方法的模拟结果进行了对比,结果表明两者符合得很好.本文方法可用于横向各向同性介质的深度偏移及地震层析成像的深入研究.     

基于射线追踪技术计算地震定位中震源轨迹的改进方法

使用震源轨迹确定震源位置不仅稳健而且直观,但当介质复杂时震源轨迹难以给出解析解.基于最小走时树射线追踪技术计算震源轨迹的方法(以轨迹所在的残差场中残差最小的点(初始点)至残差较小的点(震源轨迹代表点)的射线路径表示震源轨迹)适用于复杂速度模型,但尚不能正确计算由多段组成的震源轨迹,同时兼顾计算轨迹的完整性和精细性较为困难,计算参数设置烦琐不适于大批量数据的自动处理.针对该方法存在的问题,本文对其进行了改进:(1)采用一种"削皮"算法选取震源轨迹所经过的模型单元的节点作为轨迹代表点;(2)将残差较小的区域作为震源轨迹计算区域(该区域依轨迹分布自适应地划分为若干个连通区域),从未计算的轨迹代表点中选取残差最小者作为射线路径初始点,利用最小走时树算法依次计算所有连通区域内的震源轨迹;(3)通过去掉较短的不再分叉的射线路径使震源轨迹更为精细.虚拟和真实事件的算例表明,改进方法有效克服了原方法的不足,可便捷地计算复杂速度模型中事件的震源轨迹,计算的轨迹精细且较完整.